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
  67. température';
  68. progx = prog 0. 0.001785 0.00386 0.00701 0.00984
  69. 0.0132 0.0167 0.0259 0.0491 0.0994;
  70. progy = prog 0. 350.23 470.05 553.73 610.14
  71. 635.94 663.59 717.97 795.39 877.42;
  72. ev1 = evol manu eps progx sig (progy*1.0);
  73. ev2 = evol manu eps progx sig (progy*0.8);
  74. ev3 = evol manu eps progx sig (progy*0.6);
  75. ev4 = evol manu eps progx sig (progy*0.4);
  76. ev5 = evol manu eps progx sig (progy*0.2);
  77. tra = nuage comp T 0. 200. 400. 600. 800.
  78. comp TRAC ev1 ev2 ev3 ev4 ev5;
  79. SI (NEG GRAPH 'N');
  80. TAB1=TABL;
  81. TAB1.1 = 'MARQ LOSA ';
  82. TAB1.2 = 'MARQ TRIB ';
  83. TAB1.3 = 'MARQ CROI ';
  84. TAB1.4 = 'MARQ TRIA ';
  85. TAB1.5 = 'MARQ CARR ';
  86. TAB1.'TITRE'=TABL;
  87. TAB1.'TITRE'. 1 = MOT 'T=800.C';
  88. TAB1.'TITRE'. 2 = MOT 'T=600.C';
  89. TAB1.'TITRE'. 3 = MOT 'T=400.C';
  90. TAB1.'TITRE'. 4 = MOT 'T=200.C';
  91. TAB1.'TITRE'. 5 = MOT 'T=0.C';
  92. DESS (ev5 et ev4 et ev3 et ev2 et ev1) lege tab1;
  93. FINSI;
  94. **** 3.2) Module D'Young en fonction de la température
  95. TITR 'Module d Young en fonction de la
  96. température';
  97. y1=calyo ev1;
  98. y2=calyo ev2;
  99. y3=calyo ev3;
  100. y4=calyo ev4;
  101. y5=calyo ev5;
  102. evy = evol manu T (prog 0. 200. 400. 600. 800.)
  103. YOUN (prog y1 y2 y3 y4 y5);
  104. SI (NEG GRAPH 'N'); DESS evy; FINSI;
  105. **** 3.3) Alpha en fonction de la température
  106. TITR 'Dilatation thermique en fonction de la
  107. température';
  108. eva = evol manu T (prog 0. 200. 400. 600. 800.)
  109. ALPH ((prog 10.4 11.9 13.0 13.8 14.25)*8E-6);
  110. SI (NEG GRAPH 'N'); DESS eva; FINSI;
  111. **** 3.4) la section du tuyau
  112. s = pi*((r*r) - ((r - e)*(r - e)));
  113. * 4) maillage
  114. TITR 'Tuyauderie encastré soumis à un choc
  115. thermique';
  116. p1 = 0 0 0; p2 = L 0 0;
  117. d1= p1 d 10 p2;
  118. SI (NEG GRAPH 'N'); TRAC D1; FINSI;
  119. *
  120. * 5) Température de référence
  121. TREFE = MANU CHPO d1 1 T 0.;
  122. *
  123. * 6) modele, materiaux et caracteristiques geometriques
  124. modti=MODE d1 mecanique elastique plastique
  125. isotrope tuya;
  126. mater=MATE modti 'YOUN' evy 'NU' 0.3 'ALPH' eva
  127. 'TRAC' tra;
  128. carac= 'CARA' modti 'EPAI' e 'RAYO' r 'CISA' 1. ;
  129. mc=mater et carac;
  130. *
  131. * 6) bloquage
  132. bloqtot=bloq depl rota (p1 et p2) ;
  133. *
  134. * 8) chargement thermique
  135. TBTEMPE = TABLE;
  136. TBTEMPE.0 =TREFE;
  137. TBTEMPE.1 =MANU CHPO d1 1 T 200.;
  138. TBTEMPE.2 =MANU CHPO d1 1 T 400.;
  139. TBTEMPE.3 =MANU CHPO d1 1 T 600.;
  140. TBTEMPE.4 =MANU CHPO d1 1 T 800.;
  141. TEMP = TABLE;
  142. TEMP.0 = 0.;
  143. TEMP.1 = 2.;
  144. TEMP.2 = 4.;
  145. TEMP.3 = 6.;
  146. TEMP.4 = 8.;
  147. EVT = EVOL MANU TEMPS (prog 0. 2. 4. 6. 8.)
  148. T (prog 0. 200. 400. 600. 800.);
  149. CH2 = CHAR T TEMP TBTEMPE;
  150. *
  151. * 9) chargement mecanique
  152. f1 = forc (0. 0. 0.) p1;
  153. ch1 = char meca f1 (evol manu (prog 0. 10.)
  154. (prog 0. 0.));
  155. chat = ch2 et ch1;
  156. *
  157. * 10) calcul nonlinéaire
  158. TAB1 = TABLE;
  159. TAB1.'TEMPERATURE_REFERENCE'=TREFE;
  160. TAB1.'MODELE' = modti;
  161. TAB1.'CARACTERISTIQUES' = MC;
  162. TAB1.'BLOCAGES_MECANIQUES' = BLOQTOT;
  163. TAB1.'TEMPS_CALCULES' = (prog 0. pas 1. 8.);
  164. TAB1.'CHARGEMENT' = CHAT;
  165. PASAPAS TAB1;
  166. *
  167. * 11) comparaison et message
  168. TBEVOL=TABL;
  169. TBEVOL.1=(EV1 + EV2)/2.;
  170. TBEVOL.2=EV2;
  171. TBEVOL.3=(EV2 + EV3)/2.;
  172. TBEVOL.4=EV3;
  173. TBEVOL.5=(EV3 + EV4)/2.;
  174. TBEVOL.6=EV4;
  175. TBEVOL.7=(EV4 + EV5)/2.;
  176. TBEVOL.8=EV5;
  177. TBCONT=TAB1.CONTRAINTES;
  178. PCALCUL=PROG; PANALYT=PROG; PEPSI=PROG;
  179. REPETER B1 ((DIME TBCONT) - 1);
  180. V1 = (EXTR (TBCONT.(&B1)) EFFX 1 1 1)/S;
  181. PCALCUL=PCALCUL ET (PROG (ABS V1));
  182. T1 = IPOL (EXTR EVT ABSC) (EXTR EVT ORDO)
  183. (TAB1.TEMPS.&B1);
  184. Y1 = IPOL (EXTR EVY ABSC) (EXTR EVY ORDO) T1;
  185. A1 = IPOL (EXTR EVA ABSC) (EXTR EVA ORDO) T1;
  186. PEPSI=PEPSI ET (PROG (A1*T1));
  187. EVC=TBEVOL.&B1;
  188. V2 = IPOL (EXTR EVC ABSC) (EXTR EVC ORDO) (A1*T1);
  189. PANALYT=PANALYT ET (PROG V2);
  190. FIN B1;
  191. ECALCUL=EVOL ROUG MANU 'Epsillon' PEPSI
  192. 'Sigma(Mpa)' PCALCUL;
  193. EANALYT=EVOL TURQ MANU 'Epsillon' PEPSI
  194. 'Sigma(Mpa)' PANALYT;
  195. TAB1 = TABL;
  196. TAB1.1 = 'MARQ LOSA ';
  197. TAB1.2 = 'MARQ TRIB ';
  198. TAB1.'TITRE' = TABL;
  199. TAB1.'TITRE'. 1 = MOT 'CASTEM2000';
  200. TAB1.'TITRE'. 2 = MOT 'Analytique';
  201. SI (NEG GRAPH 'N');
  202. DESS (ECALCUL et EANALYT) lege tab1;
  203. FINSI;
  204. ENEANA=EXTR 1 (INTG EANALYT);
  205. ENEK20=EXTR 1 (INTG ECALCUL);
  206. ERR1=ABS ((ENEANA - ENEK20)/ENEANA);
  207. mess '***********************';
  208. mess 'Energie de déformation cumulée';
  209. mess 'k2000,theorie,erreur: ' ENEK20 ENEANA ERR1;
  210. mess '***********************';
  211. si (ERR1 > 5.5d-3);
  212. erre 5;
  213. sinon;
  214. erre 0;
  215. finsi;
  216. fin;
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  

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