Télécharger la plaquette Cast3M


fr en

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