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

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