Télécharger plas12.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas12.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Plas12.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. GRAPH='N';
  8. SAUT PAGE ;
  9. *
  10. *******************************************************
  11. * TEST PLAS12
  12. *
  13. * Un tuyauderie encastré en deux extremités
  14. * soumis à un choc thermique.
  15. *
  16. * GEOMETRIE : Longueur du tuyau l : 1000. mm
  17. * Rayon extérieure du tuyau r : 10. mm
  18. * Epaisseur du tuyau e : 100. mm
  19. * MATERIAU : Module d'Young E : dépendant
  20. * de la
  21. * température
  22. * Coeficient de Poisson Nu : 0.3
  23. * dilatation thermique Alpha : dépendant
  24. * de la
  25. * température
  26. * Courbe de traction : dépendant
  27. * de la
  28. * température
  29. * CHARGEMENT: Une tempérarure T uniforme passant
  30. * de 0°C à 800°C
  31. *
  32. * encastrement encastrement
  33. * |-------------------------------------|
  34. * p1 == CHOC THERMIQUE == p2
  35. *
  36. *******************************************************
  37. * Objectif du test:
  38. * Verifier la procédure PASAPAS avec un matériau
  39. * dépendant de la température.
  40. *******************************************************
  41. * Validation du calcul:
  42. * On compare la contrainte de PASAPAS avec la solution
  43. * exacte : EPSILLON = Alpha * (Delta T)
  44. * SIGMA = Interpolation de la valeur EPSILLON
  45. * sur la courbe de traction
  46. *******************************************************
  47. *******************************************************
  48. * procédure calculant le module D'Young avec une courbe
  49. * de traction
  50. DEBP CALYO EV1*EVOLUTION;
  51. X1=EXTR EV1 ABSC;
  52. Y1=EXTR EV1 ORDO;
  53. YOU1=(EXTR Y1 2)/(EXTR X1 2);
  54. FINP YOU1;
  55. *******************************************************
  56. *
  57. * 1) option de calcul
  58. opti dime 3 elem seg2;
  59. *
  60. * 2) données géométriques
  61. l=1000.; e = 10; r = 100.;
  62. *
  63. * 3) données de caractéristiques du matériau
  64. **** 3.1) courbe de traction en fonction de la
  65. **** température
  66. TITR 'Courbe de traction en fonction de la température';
  67. progx = prog 0. 0.001785 0.00386 0.00701 0.00984 0.0132 0.0167 0.0259 0.0491 0.0994;
  68. progy = prog 0. 350.23 470.05 553.73 610.14 635.94 663.59 717.97 795.39 877.42;
  69. ev1 = evol turq manu eps progx sig (progy*1.0);
  70. ev2 = evol turq manu eps progx sig (progy*0.8);
  71. ev3 = evol turq manu eps progx sig (progy*0.6);
  72. ev4 = evol turq manu eps progx sig (progy*0.4);
  73. ev5 = evol turq manu eps progx sig (progy*0.2);
  74. tra = nuage comp T 0. 200. 400. 600. 800. comp ECRO ev1 ev2 ev3 ev4 ev5;
  75. **** 3.2) Module D'Young en fonction de la température
  76. TITR 'Module d Young en fonction de la température';
  77. y1=calyo ev1;
  78. y2=calyo ev2;
  79. y3=calyo ev3;
  80. y4=calyo ev4;
  81. y5=calyo ev5;
  82. evy = evol manu T (prog 0. 200. 400. 600. 800.) YOUN (prog y1 y2 y3 y4 y5);
  83. SI (NEG GRAPH 'N'); DESS evy; FINSI;
  84. **** 3.3) Alpha en fonction de la température
  85. TITR 'Dilatation thermique en fonction de la température';
  86. eva = evol manu T (prog 0. 200. 400. 600. 800.) ALPH ((prog 10.4 11.9 13.0 13.8 14.25)*8E-6);
  87. SI (NEG GRAPH 'N'); DESS eva; FINSI;
  88. **** 3.4) la section du tuyau
  89. s = pi*((r*r) - ((r - e)*(r - e)));
  90. **** 3.5) Courbe d'ecrouissage en fonction de la temperature
  91. TITR 'Courbe d ecrouissage en fonction de la température';
  92. progy = progy enle 1 ;
  93. progx = (progx enle 1) - (progy / y1) ;
  94. ec1 = evol roug manu eps progx sig (progy*1.0);
  95. ec2 = evol roug manu eps progx sig (progy*0.8);
  96. ec3 = evol roug manu eps progx sig (progy*0.6);
  97. ec4 = evol roug manu eps progx sig (progy*0.4);
  98. ec5 = evol roug manu eps progx sig (progy*0.2);
  99. cec = nuage comp T 0. 200. 400. 600. 800. comp ECRO ec1 ec2 ec3 ec4 ec5;
  100. SI (NEG GRAPH 'N');
  101. TAB1=TABL;
  102. TAB1.1 = 'MARQ LOSA ';
  103. TAB1.2 = 'MARQ TRIB ';
  104. TAB1.3 = 'MARQ CROI ';
  105. TAB1.4 = 'MARQ TRIA ';
  106. TAB1.5 = 'MARQ CARR ';
  107. TAB1.6 = 'TIRR MARQ LOSA ';
  108. TAB1.7 = 'TIRR MARQ TRIB ';
  109. TAB1.8 = 'TIRR MARQ CROI ';
  110. TAB1.9 = 'TIRR MARQ TRIA ';
  111. TAB1.10 = 'TIRR MARQ CARR ';
  112. TAB1.'TITRE'=TABL;
  113. TAB1.'TITRE'. 1 = MOT ' Traction T=800.C';
  114. TAB1.'TITRE'. 2 = MOT ' Traction T=600.C';
  115. TAB1.'TITRE'. 3 = MOT ' Traction T=400.C';
  116. TAB1.'TITRE'. 4 = MOT ' Traction T=200.C';
  117. TAB1.'TITRE'. 5 = MOT ' Traction T=0.C';
  118. TAB1.'TITRE'. 6 = MOT 'Ecrouissage T=800.C';
  119. TAB1.'TITRE'. 7 = MOT 'Ecrouissage T=600.C';
  120. TAB1.'TITRE'. 8 = MOT 'Ecrouissage T=400.C';
  121. TAB1.'TITRE'. 9 = MOT 'Ecrouissage T=200.C';
  122. TAB1.'TITRE'. 10 = MOT 'Ecrouissage T=0.C';
  123. ev0 = ev5 et ev4 et ev3 et ev2 et ev1 ;
  124. ec0 = ec5 et ec4 et ec3 et ec2 et ec1 ;
  125. DESS (ev0 et ec0) lege tab1;
  126. FINSI;
  127. * 4) maillage
  128. TITR 'Tuyauderie encastré soumis à un choc thermique';
  129. p1 = 0 0 0; p2 = L 0 0;
  130. d1= p1 d 10 p2;
  131. SI (NEG GRAPH 'N'); TRAC D1; FINSI;
  132. *
  133. * 6) modele, materiaux et caracteristiques geometriques
  134. modti= MODE d1 mecanique elastique plastique isotrope tuya;
  135. mater= MATE modti 'YOUN' evy 'NU' 0.3 'ALPH' eva 'TALP' 0. 'TREF' 0. 'ECRO' cec;
  136. carac= 'CARA' modti 'EPAI' e 'RAYO' r 'CISA' 1. ;
  137. mc=mater et carac;
  138. *
  139. * 6) bloquage
  140. bloqtot=bloq depl rota (p1 et p2) ;
  141. *
  142. * 8) chargement thermique
  143. TBTEMPE = TABLE;
  144. TBTEMPE.0 =MANU CHPO d1 1 T 0.;
  145. TBTEMPE.1 =MANU CHPO d1 1 T 200.;
  146. TBTEMPE.2 =MANU CHPO d1 1 T 400.;
  147. TBTEMPE.3 =MANU CHPO d1 1 T 600.;
  148. TBTEMPE.4 =MANU CHPO d1 1 T 800.;
  149. TEMP = TABLE;
  150. TEMP.0 = 0.;
  151. TEMP.1 = 2.;
  152. TEMP.2 = 4.;
  153. TEMP.3 = 6.;
  154. TEMP.4 = 8.;
  155. EVT = EVOL MANU TEMPS (prog 0. 2. 4. 6. 8.) T (prog 0. 200. 400. 600. 800.);
  156. CH2 = CHAR T TEMP TBTEMPE;
  157. *
  158. * 9) chargement mecanique
  159. f1 = forc (0. 0. 0.) p1;
  160. ch1 = char meca f1 (evol manu (prog 0. 10.) (prog 0. 0.));
  161. chat = ch2 et ch1;
  162. *
  163. * 10) calcul nonlinéaire
  164. TAB1 = TABLE;
  165. TAB1.'MODELE' = modti;
  166. TAB1.'CARACTERISTIQUES' = MC;
  167. TAB1.'BLOCAGES_MECANIQUES' = BLOQTOT;
  168. TAB1.'TEMPS_CALCULES' = (prog 0. pas 1. 8.);
  169. TAB1.'CHARGEMENT' = CHAT;
  170. PASAPAS TAB1;
  171. *
  172. * 11) comparaison et message
  173. TBEVOL=TABL;
  174. TBEVOL.1=(EV1 + EV2)/2.;
  175. TBEVOL.2=EV2;
  176. TBEVOL.3=(EV2 + EV3)/2.;
  177. TBEVOL.4=EV3;
  178. TBEVOL.5=(EV3 + EV4)/2.;
  179. TBEVOL.6=EV4;
  180. TBEVOL.7=(EV4 + EV5)/2.;
  181. TBEVOL.8=EV5;
  182. TBCONT=TAB1.CONTRAINTES;
  183. PCALCUL=PROG; PANALYT=PROG; PEPSI=PROG;
  184. REPETER B1 ((DIME TBCONT) - 1);
  185. V1 = (EXTR (TBCONT.(&B1)) EFFX 1 1 1)/S;
  186. PCALCUL=PCALCUL ET (PROG (ABS V1));
  187. T1 = IPOL (EXTR EVT ABSC) (EXTR EVT ORDO) (TAB1.TEMPS.&B1);
  188. Y1 = IPOL (EXTR EVY ABSC) (EXTR EVY ORDO) T1;
  189. A1 = IPOL (EXTR EVA ABSC) (EXTR EVA ORDO) T1;
  190. PEPSI=PEPSI ET (PROG (A1*T1));
  191. EVC=TBEVOL.&B1;
  192. V2 = IPOL (EXTR EVC ABSC) (EXTR EVC ORDO) (A1*T1);
  193. PANALYT=PANALYT ET (PROG V2);
  194. FIN B1;
  195. ECALCUL=EVOL ROUG MANU 'Epsillon' PEPSI 'Sigma(Mpa)' PCALCUL;
  196. EANALYT=EVOL TURQ MANU 'Epsillon' PEPSI 'Sigma(Mpa)' PANALYT;
  197. TAB1 = TABL;
  198. TAB1.1 = 'MARQ LOSA ';
  199. TAB1.2 = 'MARQ TRIB ';
  200. TAB1.'TITRE' = TABL;
  201. TAB1.'TITRE'. 1 = MOT ' Cast3M';
  202. TAB1.'TITRE'. 2 = MOT 'Analytique';
  203. SI (NEG GRAPH 'N');
  204. DESS (ECALCUL et EANALYT) lege tab1;
  205. FINSI;
  206. ENEANA=EXTR 1 (INTG EANALYT);
  207. ENEK20=EXTR 1 (INTG ECALCUL);
  208. ERR1=ABS ((ENEANA - ENEK20)/ENEANA);
  209. opti echo 0 ;
  210. mess '***********************';
  211. mess 'Energie de déformation cumulée';
  212. mess 'k2000,theorie,erreur: ' ENEK20 ENEANA ERR1;
  213. mess '***********************';
  214. opti echo 1 ;
  215. si (ERR1 > 5.5d-3);
  216. erre 5;
  217. sinon;
  218. erre 0;
  219. finsi;
  220. fin;
  221.  
  222.  
  223.  
  224.  
  225.  

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