Télécharger rupt11.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rupt11.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ****************************************************
  5. * *
  6. * VALIDATION DES PROCEDURES GTHETA ET T_PITETA *
  7. * PAR UNE PLAQUE EN TRACTION PURE. SOLUTION *
  8. * DE REFERENCE : ISIDA, On the tension of a *
  9. * strip with a central elliptical hole. *
  10. * Tran. Jap. Soc. Mec. Engng., vol. 21 (1955). *
  11. * *
  12. * GEOMETRIE : *
  13. * *
  14. * longueur totale de la plaque : 700 MM *
  15. * largeur totale de la plaque : 200 MM *
  16. * longueur totale de la fissure : 100 MM *
  17. * *
  18. * MATERIAU : *
  19. * *
  20. * Module d'Young : 19775 Kgf/mm² *
  21. * Coefficient de poisson : 0,3 *
  22. * *
  23. * CHARGEMENT : *
  24. * *
  25. * Traction pure : 100 Kgf/mm² *
  26. * *
  27. * SOLUTION ANALYTIQUE : *
  28. * *
  29. * Le taux de restitution d'énergie (G_THETA) : *
  30. * G = 102,23 Kgf/mm *
  31. * La dérivée seconde de l'énergie (T_PITETA) : *
  32. * dG/da = 3,650 Kgf/mm² *
  33. * *
  34. ****************************************************
  35. DEBP PRMAIL P0*POINT LI11*MAILLAGE LI12*MAILLAGE LI13*MAILLAGE
  36. LI14*MAILLAGE NVTAI*FLOTTANT NDC*ENTIER;
  37. P21=P0 MOIN (NVTAI 0.);P22=P21 PLUS (0. NVTAI);
  38. P23=P22 PLUS (NVTAI 0.);P24=P23 PLUS (NVTAI 0.);
  39. P25=P24 MOIN (0. NVTAI);
  40. LI21=P21 D NDC P22;LI22=P22 D NDC P23;
  41. LI23=P23 D NDC P24;LI24=P24 D NDC P25;
  42. SU21=COUT LI11 LI21;SU22=COUT LI12 LI22;
  43. SU23=COUT LI13 LI23;SU24=COUT LI14 LI24;
  44. SUT=SU21 ET SU22 ET SU23 ET SU24;
  45. FINP LI21 LI22 LI23 LI24 SUT;
  46. *
  47. OPTI DIME 2 ELEM QUA8 MODE PLAN DEFO;
  48. opti epsilon lineaire;
  49. PB=50. 0.; VTAI=0.5; TVTAI=TABLE;
  50. ***************************************
  51. * MAILLAGE CARRE EN POINTE DE FISSURE *
  52. ***************************************
  53. CVTAI=10. * VTAI; TVTAI=TABL; TNDC=TABL;
  54. TVTAI.1= 12. * VTAI; TNDC.1=8;
  55. TVTAI.2= 15. * VTAI; TNDC.2=6;
  56. TVTAI.3= 20. * VTAI; TNDC.3=5;
  57. TVTAI.4= 28. * VTAI; TNDC.4=4;
  58. TVTAI.5= 40. * VTAI; TNDC.5=3;
  59. TVTAI.6= 58. * VTAI; TNDC.6=3;
  60. TVTAI.7= 78. * VTAI; TNDC.7=3;
  61. TVTAI.8=100. * VTAI; TNDC.8=2;
  62. P11=PB MOIN (CVTAI 0.); P12=P11 PLUS (0. CVTAI);
  63. MCARR1=(P11 D 10 P12) TRAN 10 (CVTAI 0.);
  64. MCARR2=(INVE (COTE 3 MCARR1)) TRAN 10 (CVTAI 0.);
  65. MCARR0=MCARR1 ET MCARR2;
  66. PB=MCARR0 POIN PROCH PB;
  67. MFOND=MCARR0;
  68. *********************************************
  69. * ON ELARGIT ET ON MET A LA PLACE DE SUFISS *
  70. *********************************************
  71. LI11=COTE 1 MCARR1;LI12=COTE 2 MCARR1;
  72. LI13=COTE 2 MCARR2;LI14=COTE 3 MCARR2;
  73. I=0;
  74. REPETER BCOU (DIME TVTAI);
  75. I=I + 1;
  76. NVTAI=TVTAI.I;NDC=TNDC.I;
  77. LI11 LI12 LI13 LI14 SUT2=PRMAIL PB LI11 LI12 LI13 LI14 NVTAI NDC;
  78. MFOND=MFOND ET SUT2;
  79. FIN BCOU;
  80. ELIM 0.01 MFOND;
  81. SUT = (LI12 ET LI13) TRAN -12 DINI 8 DFIN 30 (0. (350. - 50.));
  82. SUT = SUT ET MFOND;
  83. CSUT = CONT SUT;
  84. ELIM 0.01 SUT;
  85. *******************************************
  86. * PREPARATION POUR CALCUL SUR MAILLES 1MM *
  87. *******************************************
  88. LICOT=CSUT ELEM COMP (SUT POIN PROC (0. 350.)) (SUT POIN PROC (0. 0.));
  89. LIGAM=CSUT ELEM COMP (SUT POIN PROC (50. 0.)) (SUT POIN PROC (100. 0.));
  90. LIFIS=CSUT ELEM COMP (SUT POIN PROC (0. 0.)) (SUT POIN PROC (50. 0.));
  91. LIHAU=CSUT ELEM COMP (SUT POIN PROC (100. 350.))
  92. (SUT POIN PROC (0. 350.));
  93. OB1= modeli sut mecanique elastique;
  94. MA1=MATE OB1 YOUN 19775. NU 0.3 ALPHA 1.5E-5;
  95. RI1=RIGI MA1 OB1;
  96. CDL1=(BLOQ UX LICOT) ET (BLOQ UY LIGAM);
  97. FO1=PRES MASS OB1 -100. LIHAU;
  98. U=RESO (RI1 ET CDL1) FO1;
  99. SIG=SIGM MA1 U OB1;
  100. titr 'Fissure deformee sous la pression';
  101. vec1=vect fo1 fx fy 0.01 roug;
  102. *trac (defo sut u vec1) face (cont sut);
  103. ***** Solution analytique ******
  104. COE1 = (0.91/19775.)*100.*100.*pi;
  105. COE2 = pi*50./200;
  106. COS1 = cos (COE2*180./pi); SIN1 = sin (COE2*180./pi);
  107. G_ANA = COE1*50.*(1./COS1);
  108. DGDA_ANA = COE1*((1./COS1) + (COE2*SIN1/(COS1*COS1)));
  109. ********************************
  110. * Appel à la procédure G_THETA *
  111. ********************************
  112. TAB = TABLE ;
  113. TAB.'MAILLAGE'= SUT;
  114. TAB.'FISSURE' = LIFIS;
  115. TAB.'FONDFISS' = PB;
  116. AFF = TABLE; AFF.1 = OB1;
  117. TAB2=TABLE;
  118. TAB2.'DEPLACEMENT'=U;
  119. TAB2.'CONTRAINTE'=SIG;
  120.  
  121. SUPTAB = TABLE ;
  122. SUPTAB.'OBJECTIF' = MOT 'J';
  123. SUPTAB.'LEVRE_SUPERIEURE' = lifis;
  124. SUPTAB.'FRONT_FISSURE' = PB ;
  125. SUPTAB.'MODELE' = ob1;
  126. SUPTAB.'CARACTERISTIQUES' = ma1;
  127. SUPTAB.'SOLUTION_RESO' = u;
  128. SUPTAB.'CHARGEMENTS_MECANIQUES' = FO1;
  129. i = 5;
  130. repeter bbbb 1;
  131. i = i + 1;
  132. SUPTAB.'COUCHE' = i;
  133. G_THETA SUPTAB;
  134. G_NUM = SUPTAB.'RESULTATS';
  135. ERR1 = (ABS ((G_NUM - G_ANA)/G_ANA))*100.;
  136. MESS 'G_NUM =' G_NUM 'Erreur =' (CHAI ERR1 ' %');
  137. fin bbbb;
  138. *********************************
  139. * Appel à la procédure T_PITETA *
  140. *********************************
  141. SUPTAB.'OBJECTIF' = MOT 'DJ/DA';
  142. SUPTAB.'BLOCAGES_MECANIQUES'=CDL1 ;
  143. i = 5;
  144. repeter bbbb 1;
  145. i = i + 1;
  146. SUPTAB.'COUCHE' = i;
  147. G_THETA SUPTAB;
  148. DGDA_NUM = SUPTAB.'RESULTATS';
  149. ERR2 = (ABS ((DGDA_NUM - DGDA_ANA)/DGDA_ANA))*100.;
  150. MESS 'DGDA_NUM =' DGDA_NUM 'Erreur =' (CHAI ERR2 ' %');
  151. fin bbbb;
  152. **************************************
  153. ************ Arrêt du calcul *********
  154. **************************************
  155. SI ((ERR1 > 0.65) OU (ERR2 > 0.5));
  156. ERRE 5;
  157. SINO;
  158. ERRE 0;
  159. FINSI;
  160. FIN;
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  

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