Télécharger rupt8.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rupt8.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Rupt8.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. ****************************************************
  8. * *
  9. * VALIDATION DE LA METHODE G_THETA EN CAS *
  10. * D'UNE PLAQUE EN TRACTION PURE. SOLUTION *
  11. * DE REFERENCE : Compendium of STRESS INTENSITY *
  12. * FACTORS by Rooke & Cartwright. Section Plates *
  13. * and Shells "Central crack in a finite width *
  14. * plate : uniform bending moment" *
  15. * *
  16. * *
  17. * GEOMETRIE : *
  18. * *
  19. * longueur de la fissure : 2a = 2 x 0.30 M *
  20. * longueur de la plaque : 2b = 2 x 1.00 M *
  21. * Hauteur de la plaque : 2h = 2 x 5.25 M *
  22. * Epaisseur de la plaque : e = 0.001 M *
  23. * *
  24. * MATERIAU : *
  25. * *
  26. * Module d'Young : E = 2.E11 *
  27. * Coefficient de poisson : NU = 0.3 *
  28. * *
  29. * CHARGEMENT : *
  30. * *
  31. * Traction pure : SIG = 5.E6 N/M2 *
  32. * *
  33. * SOLUTION ANALYTIQUE : *
  34. * *
  35. * Le taux de restitution d'energie : G = 129.5 *
  36. * *
  37. ****************************************************
  38. opti echo 1;
  39. OPTI DIME 3 ELEM TRI3 DENS 1.;
  40. N4 = 6; N8 = 13;
  41. N12 = (N4 + N8)*(-1);
  42. O1 = 0. 0. 1000.;
  43. ADIM = 0.3; BDIM = 1.; TAFIS = ADIM / 33.;
  44. HH1 = BDIM / 2.0;
  45. P1 = 0. 0. 0.; P2 = BDIM 0. 0.; P3 = ADIM 0. 0.;
  46. P3G = (ADIM - TAFIS) 0. 0.; P3D = (ADIM + TAFIS) 0. 0.;
  47. ALF1 = ATG HH1 ADIM; ALF2 = ATG HH1 (BDIM - ADIM);
  48. X1 = TAFIS * (COS ALF1); Y1 = TAFIS * (SIN ALF1);
  49. P31 = (ADIM - X1) Y1 0.;
  50. X1 = TAFIS * (COS ALF2); Y1 = TAFIS * (SIN ALF2);
  51. P32 = (ADIM + X1) Y1 0.;
  52. C1 = C N4 P3G P3 P31; C2 = C N8 P31 P3 P32;
  53. C3 = C N4 P32 P3 P3D;
  54. SU1 = COUT (C1 ET C2 ET C3) P3;
  55. P4 = 0. HH1 0.; P5 = BDIM HH1 0.;
  56. L1 = D N4 P1 P4 DINI (2.0*TAFIS) DFIN (8.0*TAFIS);
  57. L2 = D N12 P4 P31 DINI (9.0*TAFIS) DFIN (0.3*TAFIS);
  58. L3 = D N12 P1 P3G DINI (3.0*TAFIS) DFIN (1.5*TAFIS);
  59. L4 = D N8 P4 P5 ;
  60. L5 = D N12 P5 P32 DINI (14.0*TAFIS) DFIN (0.3*TAFIS);
  61. L6 = D N4 P2 P5 DINI (2.0*TAFIS) DFIN (10.0*TAFIS);
  62. L7 = D N12 P2 P3D DINI (10.0*TAFIS) DFIN (1.5*TAFIS);
  63. SU2 = DALL L1 L3 C1 L2;
  64. SU3 = DALL L2 L4 L5 C2;
  65. SU4 = DALL L6 L5 C3 L7;
  66. SU5 = L4 TRAN 7 (0. (1.0*BDIM) 0.);
  67. L8 = COTE 3 SU5; Y1 = COOR 2 (POINT L8 FINA);
  68. P6 = 0. (Y1 + (BDIM/4.)) 0.;
  69. P7 =BDIM (Y1 + (BDIM/4.)) 0.;
  70. L9 = D N4 P7 P6;
  71. SU6 = COUT L8 L9;
  72. SU7 = L9 TRAN 8 (0. (3.5*BDIM) 0.);
  73. L10 = COTE 3 SU7; Y1 = COOR 2 (POINT L10 FINA);
  74. L11 = D 1 P3 P3D;
  75. SUT = COUL TURQ (SU1 ET SU2 ET SU3 ET SU4 ET SU5 ET
  76. SU6 ET SU7);
  77. ELIM 1.E-8 SUT;
  78. L1 = (CONT SUT) ELEM APPU (SUT POIN
  79. DROI (0. 0. 0.) (BDIM 0. 0.) 1.E-8);
  80. L2 = (CONT SUT) ELEM APPU (SUT POIN
  81. DROI (BDIM 0. 0.) (BDIM Y1 0.) 1.E-8);
  82. L3 = (CONT SUT) ELEM APPU (SUT POIN
  83. DROI (BDIM Y1 0.) (0. Y1 0.) 1.E-8);
  84. L4 = (CONT SUT) ELEM APPU (SUT POIN
  85. DROI (0. Y1 0.) (0. 0. 0.) 1.E-8);
  86. PINI3 = POIN L3 INIT; PFIN3 = POIN L3 FINA;
  87. PMIL3 = (CHAN POI1 L3) DIFF (PINI3 ET PFIN3);
  88. MOEXT = 1. / ((NBNO L3) - 1);
  89. MOMIL = MOEXT * (NBNO PMIL3);
  90. ******
  91. ****** SOLUTION CONTRAINTE ET DEPLACEMENT
  92. ******
  93. EPAI_T = 0.001;
  94. NB_MUL = 21;
  95. EPAI_C = EPAI_T * 1.E-6;
  96. EPAI_I = (EPAI_T - (EPAI_C*NB_MUL))/(NB_MUL - 1);
  97. ELEFIS = SUT ELEM APPU LARG P3;
  98. ELEFIS = SUT ELEM APPU LARG ELEFIS;
  99. ELEFIS = SUT ELEM APPU LARG ELEFIS;
  100. TBEPAI = TABLE; TBEXCE = TABLE;
  101. REPETER BC_MUL1 ((2*NB_MUL) - 1);
  102. III = (&BC_MUL1 - NB_MUL) / 2.;
  103. TBEXCE.&BC_MUL1 = (EPAI_I + EPAI_C)*III;
  104. SI (EGA 0. ((&BC_MUL1/2.) - (&BC_MUL1/2)) 1.E-10);
  105. TBEPAI.&BC_MUL1 = EPAI_I;
  106. SINON;
  107. TBEPAI.&BC_MUL1 = EPAI_C;
  108. FINSI;
  109. FIN BC_MUL1;
  110. REPETER BC_MUL2 ((2*NB_MUL) - 1);
  111. SI (EGA &BC_MUL2 1);
  112. MOD_MUL = MODE ELEFIS MECANIQUE ELASTIQUE
  113. DKT CONS (CHAI &BC_MUL2);
  114. MAT_MUL = MATE MOD_MUL YOUN 2.E11 NU 0.3
  115. EPAI TBEPAI.&BC_MUL2 EXCE TBEXCE.&BC_MUL2;
  116. SINON;
  117. MOD_MU1 = MODE ELEFIS MECANIQUE ELASTIQUE
  118. DKT CONS (CHAI &BC_MUL2);
  119. MAT_MU1 = MATE MOD_MU1 YOUN 2.E11 NU 0.3
  120. EPAI TBEPAI.&BC_MUL2 EXCE TBEXCE.&BC_MUL2;
  121. MOD_MUL = MOD_MUL ET MOD_MU1;
  122. MAT_MUL = MAT_MUL ET MAT_MU1;
  123. FINSI;
  124. FIN BC_MUL2;
  125. MOD_R = MODE (SUT DIFF ELEFIS) MECANIQUE ELASTIQUE DKT;
  126. MAT_R = MATE MOD_R YOUN 2.E11 NU 0.3 EPAI EPAI_T;
  127. MO1 = MOD_MUL ET MOD_R;
  128. MA1 = MAT_MUL ET MAT_R;
  129. BL1 = BLOQ UY RX (L7 ET L11);
  130. BL2 = BLOQ UX RY L4;
  131. BL3 = BLOQ UZ RZ L3;
  132. RI = BL1 ET BL2 ET BL3 ET (RIGI MA1 MO1);
  133. FO2 = ((FORC FY MOEXT (PINI3 ET PFIN3)) ET
  134. (FORC FY MOMIL PMIL3))*5.e3;
  135. FO3 = FO2;
  136. DEPTO = RESO RI FO3 ;
  137. REATO = REAC DEPTO RI;
  138. SIGTO = SIGMA MA1 MO1 DEPTO ;
  139. ******************************
  140. SUPTAB = TABLE ;
  141. SUPTAB.'OBJECTIF' = MOT 'J';
  142. SUPTAB.'LEVRE_SUPERIEURE' = (cont sut) ELEM comp p3 p1;
  143. SUPTAB.'FRONT_FISSURE' = P3 ;
  144. SUPTAB.'MODELE' = MO1;
  145. SUPTAB.'CARACTERISTIQUES' = MA1;
  146. SUPTAB.'SOLUTION_RESO' = DEPTO;
  147. SUPTAB.'CHARGEMENTS_MECANIQUES' = FO3;
  148. SUPTAB.'ELEMENT_MULTICOUCHE' = ELEFIS ELEM 1;
  149. SUPTAB.'COUCHE' = 7;
  150. G_THETA SUPTAB;
  151. * RECUPERATION DU RESULTAT ET COMPARAISONS ELEMENT 25 ;
  152. * ;
  153. G = SUPTAB.'RESULTATS'.'GLOBAL' ;
  154. Kth = 50.1E5 ;
  155. K = (2.E11 * G)**0.5 ;
  156. mess ' ' ;
  157. mess ' ' ;
  158. mess 'taux de restitution energetique G :' G ;
  159. mess 'solution calculee Kcal :' K ;
  160. mess 'solution analytique Kthe :' Kth ;
  161. * CODE FONCTIONNEMENT ;
  162. RESI = abs (( k - Kth ) / Kth ) ;
  163. SI (RESI < 5E-2) ;
  164. ERRE 0 ;
  165. SINO ;
  166. ERRE 5 ;
  167. FINSI ;
  168. *RECUPERATION DU RESULTAT ET COMPARAISONS ELEMENT 117 ;
  169. * ;
  170. SUPTAB.'ELEMENT_MULTICOUCHE' = ELEFIS ELEM 117;
  171. G_THETA SUPTAB;
  172. G = SUPTAB.'RESULTATS'.'GLOBAL' ;
  173. Kth = 50.1E5 ;
  174. K = (2.E11 * G)**0.5 ;
  175. mess ' ' ;
  176. mess ' ' ;
  177. mess 'taux de restitution energetique G :' G ;
  178. mess 'solution calculee Kcal :' K ;
  179. mess 'solution analytique Kthe :' Kth ;
  180. * CODE FONCTIONNEMENT ;
  181. RESI = abs (( k - Kth ) / Kth ) ;
  182. SI (RESI < 5E-2) ;
  183. ERRE 0 ;
  184. SINO ;
  185. ERRE 5 ;
  186. FINSI ;
  187. *RECUPERATION DU RESULTAT ET COMPARAISONS ELEMENT 118 ;
  188. * ;
  189. SUPTAB.'ELEMENT_MULTICOUCHE' = ELEFIS ELEM 118;
  190. G_THETA SUPTAB;
  191. G = SUPTAB.'RESULTATS'.'GLOBAL' ;
  192. Kth = 50.1E5 ;
  193. K = (2.E11 * G)**0.5 ;
  194. mess ' ' ;
  195. mess ' ' ;
  196. mess 'taux de restitution energetique G :' G ;
  197. mess 'solution calculee Kcal :' K ;
  198. mess 'solution analytique Kthe :' Kth ;
  199. * CODE FONCTIONNEMENT ;
  200. RESI = abs (( k - Kth ) / Kth ) ;
  201. SI (RESI < 5E-2) ;
  202. ERRE 0 ;
  203. SINO ;
  204. ERRE 5 ;
  205. FINSI ;
  206. FIN;
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  

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