Télécharger rupt17.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rupt17.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **************************************************************
  5. *
  6. * rupt17.dgibi
  7. *
  8. * CAS TEST SUR LE CALCUL DE J EN THERMOPLASTICITE
  9. * POUR FISSURE PROCHE (ou SUR) INTERFACE LIAISON BIMETALLIQUE
  10. *
  11. * CREATION : BP, 11 Aout 2011
  12. * MODIFS :
  13. *
  14. * REM : - la presence de decharges induit une legere dependance au contour
  15. * - fissure parallele a linterface => pas besoin du 'MODELES_COMPOSITES'
  16. *
  17. **************************************************************
  18.  
  19.  
  20. *******************************************************
  21. *** Options de calcul et de tracés
  22.  
  23. OPTI DIME 2 ELEM QUA8 MODE plan defo ;
  24.  
  25. *deformation lineaire
  26. OPTI 'EPSI' 'LINEAIRE';
  27.  
  28. *calcul elastique?
  29. * calel = vrai ;
  30. calel = faux ;
  31.  
  32. * Complet => modif de la finesse du maillage + du pas de chargement
  33. COMPLET = faux;
  34. * COMPLET = vrai;
  35.  
  36. * options de tracé
  37. GRAPH = faux ;
  38. *opti 'TRAC' 'PSC' 'EPTR' 5;
  39.  
  40.  
  41. *******************************************************
  42. *** Géométrie
  43.  
  44. dxelim = 1e-6 ;
  45. h = 100. ;
  46. W = 100. ;
  47. hint = 50. ;
  48. hfis = 34. ;
  49. * hfis = 50. ;
  50. * avec hfis = 50., on a la fissure dans l'interface (à tester + tard)
  51. asW = 0.5 ;
  52.  
  53.  
  54. **********************************************************************
  55. *** Données matériau
  56.  
  57. *--- Données de l'acier inox 304L (20 et 350°C)
  58. lT11 = prog 20. 350. ;
  59. lE1 = prog 197. 172. ;
  60. evE1 = evol manu 'T' lT11 'YOUN' (lE1*1e3) ;
  61. *
  62. lep201 = prog 0. 0.0006935 0.00104 0.00140
  63. 0.00198 0.00305 0.00511 0.00716 0.0112 0.0163 0.0315
  64. 0.0516 0.1020 0.152 0.303 0.503 1.004 5.006 ;
  65. lsi201 = prog 0. 136.6 165.6 178.0 192.5 207.0
  66. 219.4 227.7 240.1 254.6 287.7 323.2 391.8 456.2 545.3
  67. 630.4 767.5 1212.0 ;
  68. * Conversion courbe traction > ecrouissage
  69. lsi201 = lsi201 enle 1 ;
  70. lep201 = (lep201 enle 1) - (lsi201 / (eve1 ipol 20.)) born mini 0. ;
  71. ta201 = evol manu lep201 lsi201 ;
  72. *
  73. le3501 = prog 0. 0.0006347 0.00087 0.00119
  74. 0.00171 0.00274 0.00478 0.00682 0.0109 0.0159 0.0311
  75. 0.0513 0.1017 0.152 0.303 0.503 1.004 5.006 ;
  76. ls3501 = prog 0. 109.2 115.6 118.1 121.9 127.0
  77. 134.6 141.0 151.1 161.3 189.2 221.8 290.5 354.8 443.9
  78. 529.0 666.1 1110.6 ;
  79. * Conversion courbe traction > ecrouissage
  80. ls3501 = ls3501 enle 1 ;
  81. le3501 = (le3501 enle 1) - (ls3501 / (eve1 ipol 350.)) born mini 0. ;
  82. ta3501 = evol manu le3501 ls3501 ;
  83. *
  84. CTrac1 = NUAGE 'COMP' 'T' 20. 350.
  85. 'COMP' 'ECRO' ta201 ta3501 ;
  86.  
  87. *--- Données des aciers 16 et 18MND5 (20 - 350°C)
  88. lT12 = prog 20. 350. ;
  89. lE2 = prog 204. 180. ;
  90. evE2 = evol manu 'T' lT12 'YOUN' (lE2*1e3) ;
  91. * dans le cas COMPLET pb avec la pente de lsi202
  92. * lep202 = prog 0. 0.001691 0.001891 0.00219 0.00269
  93. * 0.00369 0.00569 0.00769 0.00969 0.0117 0.0168 0.0219
  94. * 0.0320 0.0421 0.0522 0.0623 0.0723 0.0823 0.0924
  95. * 0.102 0.303 1.00 5.00 ;
  96. * lsi202 = prog 0. 345 345 345 345 345 345 345 345 355.4
  97. * 372.6 386.4 410.6 434.2 448.1 461.2 470.8 477.7
  98. * 484.6 491.6 570.9 671.5 834.2 ;
  99. * lep202 = prog 0. 0.001691 0.00969 0.0117 0.0168 0.0219
  100. * 0.0320 0.0421 0.0522 0.0623 0.0723 0.0823 0.0924
  101. * 0.102 0.303 1.00 5.00 ;
  102. * lsi202 = prog 0. 345 346 355.4
  103. * 372.6 386.4 410.6 434.2 448.1 461.2 470.8 477.7
  104. * 484.6 491.6 570.9 671.5 834.2 ;
  105. lep202 = prog 0. (0.300/204.) 0.002 0.0025 0.008 0.0117
  106. 0.0168 0.0219
  107. 0.0320 0.0421 0.0522 0.0623 0.0723 0.0823 0.0924
  108. 0.102 0.303 1.00 5.00 ;
  109. lsi202 = prog 0. 300. 320. 330. 345. 355.4
  110. 372.6 386.4 410.6 434.2 448.1 461.2 470.8 477.7
  111. 484.6 491.6 570.9 671.5 834.2 ;
  112. * Conversion courbe traction > ecrouissage
  113. lsi202 = lsi202 enle 1 ;
  114. lep202 = (lep202 enle 1) - (lsi202 / (eve1 ipol 20.)) born mini 0. ;
  115. ta202 = evol manu lep202 lsi202 ;
  116. *
  117. le3502 = prog 0. 0.001423 0.001679 0.00206 0.0026
  118. 0.00366 0.00573 0.00778 0.00983 0.0119 0.0169 0.0220
  119. 0.0321 0.0422 0.0523 0.0623 0.0724 0.0824 0.0925
  120. 0.102 0.303 1.00 5.00 ;
  121. ls3502 = prog 0. 256.2 266.1 281.0 290.0 299.0 310.9
  122. 319.9 328.9 334.9 349.8 364.8 385.7 400.2 412.3 422.9
  123. 429.0 435.0 441.0 447.0 509.0 586.7 709.3 ;
  124. * Conversion courbe traction > ecrouissage
  125. ls3502 = ls3502 enle 1 ;
  126. le3502 = (le3502 enle 1) - (ls3502 / (eve1 ipol 350.)) born mini 0. ;
  127. ta3502 = evol manu le3502 ls3502 ;
  128. *
  129. CTrac2 = NUAGE 'COMP' 'T' 20. 350.
  130. 'COMP' 'ECRO' ta202 ta3502 ;
  131.  
  132. si(GRAPH);
  133. titr ' Courbes d ecrouissage ' ;
  134. tleg1 = table ;
  135. tleg1 . titre = table ;
  136. tleg1 . titre . 1 = ' 304L a 20 degC' ;
  137. tleg1 . titre . 2 = ' 304L a 350 degC' ;
  138. tleg1 . titre . 3 = ' 18MND5 a 20 degC' ;
  139. tleg1 . titre . 4 = ' 18MND5 a 20 degC' ;
  140. dess ((coul ta201 bleu) et (ta3501 coul roug)
  141. et (coul ta202 turq) et (ta3502 coul rose)) lege tleg1 ;
  142. dess ((coul ta201 bleu) et (ta3501 coul roug)
  143. et (coul ta202 turq) et (ta3502 coul rose))
  144. 'XBOR' 0. 0.02 lege tleg1 ;
  145. fins;
  146.  
  147.  
  148. *******************************************************
  149. *** Maillage
  150.  
  151. si(COMPLET); nW = 50 ;
  152. sino; nW = 24;
  153. fins;
  154. nhF = enti (nW * (h - hint) / h) ;
  155.  
  156. P1 = 0. (-0.5 * h) ;
  157. P2 = W (-0.5 * h) ;
  158. l1 = P1 d nW p2 ;
  159.  
  160. si (hfis ega hint) ;
  161. nhA = enti (nW * hint / h) ;
  162. SA = l1 tran nhA (0. hint) ;
  163. l2 = l1 plus (0. hint) ;
  164. SF = l2 tran nhF (0. (h - hint)) ;
  165. pfis = l2 point proc ((asW * W) (hint - (0.5 * h))) ;
  166. p3 = l2 point proc (W (hint - (0.5 * h))) ;
  167. llig = l2 elem comp pfis p3 ;
  168. elim dxelim (llig et SA) ;
  169. l3 = l1 plus (0. h) ;
  170. elim dxelim (SF et l3) ;
  171. levF = diff l2 llig ;
  172. levA = levF plus (0. 0.) ;
  173. elim dxelim (SA et levA) ;
  174. sinon ;
  175. nhA = enti (nW * hfis / h) ;
  176. SA1 = l1 tran nhA (0. hfis) ;
  177. l2 = l1 plus (0. hfis) ;
  178. SA2 = l2 tran (nhF - nhA) (0. (hint - hfis)) ;
  179. l4 = l1 plus (0. hint) ;
  180. SF = l4 tran nhF (0. (h - hint)) ;
  181. pfis = l2 point proc ((asW * W) (hint - (0.5 * h))) ;
  182. p3 = l2 point proc (W (hint - (0.5 * h))) ;
  183. llig = l2 elem comp pfis p3 ;
  184. elim dxelim (llig et SA1) ;
  185. l3 = l1 plus (0. h) ;
  186. elim dxelim (SF et l3) ;
  187. levF = diff l2 llig ;
  188. levA = levF plus (0. 0.) ;
  189. elim dxelim (SA1 et levA) ;
  190. elim dxelim (SA2 et SF) ;
  191. SA = SA1 et SA2 ;
  192. finsi ;
  193.  
  194. plaq = (SA coul bleu) et (SF coul roug) ;
  195. si(GRAPH);trac plaq ;fins;
  196.  
  197.  
  198. *******************************************************
  199. *** MODELE MECANIQUE
  200.  
  201.  
  202. *--- Modèle elastique ou elasto plastique ?
  203. si calel ;
  204. mo316 = MODELE sa MECANIQUE ELASTIQUE ISOTROPE ;
  205. mo508 = MODELE sf MECANIQUE ELASTIQUE ISOTROPE ;
  206. ma316 = materiau mo316 young evE1 nu 0.3 alph 1E-5 ;
  207. ma508 = materiau mo508 young evE2 nu 0.3 alph 1E-5 ;
  208. sinon ;
  209. mo316 = MODELE sa MECANIQUE ELASTIQUE plastique ;
  210. mo508 = MODELE sf MECANIQUE ELASTIQUE plastique ;
  211. ma316 = materiau mo316 young evE1 nu 0.3 alph 1E-5 ecro CTrac1 ;
  212. ma508 = materiau mo508 young evE2 nu 0.3 alph 1E-5 ecro CTrac2 ;
  213. finsi ;
  214. mo = mo316 et mo508 ;
  215. ma = ma316 et ma508 ;
  216.  
  217. * ----- CONDITIONS LIMITES
  218. cd1 = bloq l1 UY ;
  219. cd2 = RELA ENSE UY l3 ;
  220. cd3 = bloq pfis UX ;
  221. cdl = cd1 et cd2 et cd3 ;
  222.  
  223. *----- Chargement
  224. lx1 = prog 0. 1. ;
  225. ly1 = prog 0. 1. ;
  226. ev1 = evol manu lx1 ly1 ;
  227. PRES1 = PRES MASS mo (-100.) l3 ;
  228. *PRES1 = PRES MASS mo dxelim l3 ;
  229. CHAP = CHAR MECA PRES1 ev1 ;
  230.  
  231. * ----- Définition de la temperature
  232. Ti = 20. ; Tf = 350. ;
  233. chT1 = manu chpo plaq 1 'T' Ti ;
  234. chx = plaq coor 1 ;
  235. chT2 = chT1 - (nomc 'T' (chx / W * (Ti - Tf))) ;
  236. TT1 = table ; TT1 . 0 = 0. ; TT1 . 1 = 1. ;
  237. TT2 = table ; TT2 . 0 = chT1 ; TT2 . 1 = chT2 ;
  238. CHAT = CHAR 'T' TT1 TT2 ;
  239.  
  240. * si (COMPLET); tcha = prog 0. 0.1 pas 0.05 1. ;
  241. si (COMPLET); tcha = prog 0. pas 0.05 1. ;
  242. sino; tcha = prog 0. pas 0.20 1. ;
  243. fins;
  244.  
  245.  
  246. *******************************************************
  247. *** PROCEDURE PAS A PAS
  248.  
  249. TAB2 = TABLE;
  250. TAB2.'MODELE' = mo ;
  251. TAB2.'CARACTERISTIQUES' = ma ;
  252. TAB2.'BLOCAGES_MECANIQUES' = cdl ;
  253. TAB2.'CHARGEMENT'= CHAP et CHAT ;
  254. TAB2.'TEMPS_CALCULES' = tcha ;
  255.  
  256. PASAPAS TAB2 ;
  257.  
  258. si(GRAPH);
  259.  
  260. ii = (dime tcha) - 1 ;
  261. def = defo tab2.deplacements.ii plaq ;
  262. trac chT2 def 'TITR' 'T et deformee finale';
  263.  
  264. si (non calel);
  265. var2 = TAB2 . 'VARIABLES_INTERNES' . ii;
  266. var2prg = (prog 0.) et (exp ((prog -3. PAS 0.1 -1.) * (log 10.)));
  267. trac var2 mo def var2prg 'TITR' 'variables_internes finales';
  268. fins;
  269.  
  270. fins;
  271.  
  272.  
  273. *******************************************************
  274. *** CALCUL DU TAUX DE RESTITUTION D'ENERGIE
  275.  
  276. *--- boucle sur les contours lnc ----------------------
  277. si(COMPLET);
  278. lnc = lect 1 2 3 5 7 8 9 11 ;
  279. lcoul = mots VIOL BLEU TURQ VERT OLIV JAUN ORAN ROUG;
  280. sino;
  281. lnc = lect 1 2 3 5 7 ;
  282. lcoul = mots BLEU TURQ VERT JAUN ROUG;
  283. fins;
  284. nlc = dime lnc ; i = 0 ;
  285. tath = table ;
  286. lJM = prog ;
  287. tath.0 = tcha ;
  288.  
  289. repeter bloc1 nlc ;
  290. i = i + 1 ;
  291. ncou = extr lnc i;
  292.  
  293. SUPTAB = TABLE ;
  294. SUPTAB.'OBJECTIF'= 'J';
  295. SUPTAB.'LEVRE_SUPERIEURE' = levF ;
  296. SUPTAB.'LEVRE_INFERIEURE' = levA;
  297. SUPTAB.'COUCHE' = ncou ;
  298. SUPTAB.'FRONT_FISSURE' = pfis ;
  299. SUPTAB.'SOLUTION_PASAPAS' = TAB2;
  300. G_THETA SUPTAB ;
  301.  
  302. si graph ;
  303. dess suptab.EVOLUTION_RESULTATS
  304. 'TITR' (chai i 'couche' ncou'(G_THETA)');
  305. fins ;
  306.  
  307. mocoul1 = extr lcoul i;
  308. si(ega i 1);
  309. evotot = coul suptab.EVOLUTION_RESULTATS mocoul1;
  310. sino;
  311. evotot = evotot et (coul suptab.EVOLUTION_RESULTATS mocoul1);
  312. fins;
  313.  
  314. tath.i = extr (suptab.EVOLUTION_RESULTATS) ordo ;
  315. Jmaxi = extr tath.i (dime tath.i);
  316. lJM = lJM et Jmaxi ;
  317.  
  318. fin bloc1 ;
  319. *--- fin de boucle sur les contours lnc ----------------------
  320.  
  321.  
  322. lel = W / nW ;
  323. evJ = evol 'BLEU' manu 'dimcouche' (lel*lnc) 'J' lJM ;
  324. si graph ;
  325. dess (evJ) 'TITR' (chai 'stabilite de J');
  326.  
  327. ybor0 = 0.;
  328. ybor1 = maxi (extr evJ 'ORDO');
  329. ybor1 = (1.1 * ybor1);
  330. tdess1 = tabl;
  331. tdess1 . 1 = mot 'MARQ PLUS';
  332. dess evJ 'YBOR' ybor0 ybor1 tdess1
  333. 'TITR' (chai 'stabilite de J');
  334.  
  335. dess evotot 'YBOR' ybor0 ybor1 'TITR' (chai 'evolution de J');
  336. fins ;
  337.  
  338. * opti donn 5 trac X ;
  339.  
  340.  
  341. *******************************************************
  342. *** TEST DE BON FONCTIONNEMENT
  343.  
  344. * valeurs testees
  345. Jmax1 = maxi lJM; Jmin1 = mini lJM;
  346. mess 'Jmax1 Jmin1=' Jmax1 Jmin1;
  347.  
  348. * bornes de references (au 11 aout 2011)
  349. si(COMPLET); Jmax0 = 75.; Jmin0 = 65.;
  350. sino; Jmax0 = 80.; Jmin0 = 70.;
  351. fins;
  352.  
  353. *test
  354. SI ((Jmax1 < Jmax0) et (Jmin1 > Jmin0));
  355. ERRE 0;
  356. SINON;
  357. ERRE 5;
  358. FINSI;
  359.  
  360.  
  361.  
  362. FIN ;
  363.  
  364.  
  365.  
  366.  
  367.  
  368.  
  369.  
  370.  
  371.  
  372.  
  373.  

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