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. * 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. BETA 1.06 ;
  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=(INTG EVERR);
  93. SOM2=(INTG EVOUV);
  94. ERR_REL=SOM1/SOM2;
  95. SI (ERR_REL > 1.E-2) ;
  96.   ERREUR 5;
  97. FINSI;
  98.  
  99. ouvfiss tab1 'PRIN' ;
  100. repeter bou1 n1;
  101.   ouv1 = tab1 . OUV . (&BOU1 - 1) ;
  102.   vec1 = 'VECTEUR' ouv1 mod1 mat1 ;
  103. *  trac vect1 s1 'NCLK' ;
  104. FIN BOU1;
  105.  
  106. FIN;
  107.  
  108.  
  109.  
  110.  

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