Télécharger ouvfiss2D.dgibi

Retour à la liste

Numérotation des lignes :

  1. *
  2. *
  3. * TEST DE l'OPERATEUR OUVFISS EN 2D
  4. *
  5. * ON TIRE SUR UN BARREAU ENDOMMAGEABLE
  6. *
  7. * UN DEFAUT EST CREE POUR LOCALISER L'ENDOMMAGEMENT
  8. *
  9. * UNE SEULE FISSURE EST GENEREE, DONC L'OUVERTURE EST Delta(L) - Sigma/E0
  10. *
  11. * On UTILISE DU MAZARS A ECROUISSAGE LINEAIRE
  12. *
  13. OPTI DIME 2 ELEM QUA4 MODE PLAN CONT;
  14. h=0.02;
  15. L=0.2;
  16. *
  17. * caracteristiques
  18. *
  19. YG=3.E10;
  20. Ft=3.E6;
  21. Kt=Ft/YG;
  22. GF=300.;
  23. *
  24. * AT<-10 pour écrouissage linéaire
  25. *
  26. AT=-20.;
  27. *
  28. * BT est la déformation à laquelle la contrainte s'annule
  29. *
  30. BT=2*GF/Ft/h;
  31. DENSITE h;
  32. P1=0. 0.;
  33. P2=L 0.;
  34. P3=0. H;
  35. D1=D P1 P3;
  36. S1=D1 TRANS P2;
  37. D2=S1 COTE 3;
  38. MOD1=MODE S1 MECANIQUE ELASTIQUE ISOTROPE ENDOMMAGEMENT MAZARS;
  39. MAT1=MATE MOD1 YOUN YG NU 0.2 KTR0 (KT/0.9) ATRA
  40. AT BTRA BT ACOM 1.4 BCOM 1900.;
  41. CH1=MANU CHAM MOD1 'TYPE'
  42. 'CARACTERISTIQUES' 'POSI' 'RIGIDITE' KTR0 1 1 1 (0.1*KT);
  43. MAT1=MAT1 - CH1;
  44. CL1=BLOQ D1 UX;
  45. CL2=BLOQ P1 UY;
  46. CL3=BLOQ D2 UX;
  47. CLT=CL1 ET CL2 ET CL3;
  48. F1=DEPI CL3 1.;
  49. PROG1=PROG 0. 1.;
  50. EVOL1=EVOL MANU PROG1 PROG1;
  51. CHAR1=CHARGEMENT F1 EVOL1 'DIMP';
  52. T0=Kt*L;
  53. DELTAT=T0/10.;
  54. TF=BT*h;
  55. LT1=PROG (0.9*T0) PAS (T0/10.) TF;
  56. TAB1=TABLE;
  57. TAB1.MODELE=MOD1;
  58. TAB1.CARACTERISTIQUES=MAT1;
  59. TAB1.BLOCAGES_MECANIQUES=CLT;
  60. TAB1.CHARGEMENT=CHAR1;
  61. TAB1.TEMPS_CALCULES=LT1;
  62. PASAPAS TAB1;
  63. *
  64. * courbe globale
  65. *
  66. ev1=@global tab1 evol1 cl3 fx;
  67. *
  68. * ouverture de fissure théorique
  69. *
  70. progd=extr ev1 absc;
  71. progf=extr ev1 ordo;
  72. *
  73. * contrainte = F/h, deformation=F/(E*h),
  74. * deplacement élastique = F*L/(E*h)
  75. *
  76. progouv=progd - (progf * L / YG / h);
  77. evouv=evol manu progd 'deplacement' progouv 'ouverture';
  78. *Dess evouv;
  79. ouvfiss tab1;
  80. *
  81. * ouverture calculee
  82. *
  83. progouv2=prog;
  84. n1=dime tab1.temps;
  85. repeter bou1 n1;
  86. ouv1=extr tab1 . OUV . (&BOU1 - 1) EPXX 1 1 1;
  87. PROGOUV2=INSE PROGOUV2 &BOU1 OUV1;
  88. FIN BOU1;
  89. evouv2=evol manu progd 'deplacement' progouv2 'ouverture';
  90. *Dess evouv2;
  91. EVERR=ABS (EVOUV - EVOUV2);
  92. SOM1=EXTR 1 (INTG EVERR);
  93. SOM2=EXTR 1 (INTG EVOUV);
  94. ERR_REL=SOM1/SOM2;
  95. SI (ERR_REL > 1.E-2) ;
  96. ERREUR 5;
  97. FINSI;
  98.  
  99. FIN;
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales