Télécharger la plaquette Cast3M


fr en

Télécharger ouvfiss2D.dgibi

Retour à la liste

Numérotation des lignes : 

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

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