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 = vrai;
  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. ta201 = evol manu lep201 lsi201 ;
  69. *
  70. le3501 = prog 0. 0.0006347 0.00087 0.00119
  71. 0.00171 0.00274 0.00478 0.00682 0.0109 0.0159 0.0311
  72. 0.0513 0.1017 0.152 0.303 0.503 1.004 5.006 ;
  73. ls3501 = prog 0. 109.2 115.6 118.1 121.9 127.0
  74. 134.6 141.0 151.1 161.3 189.2 221.8 290.5 354.8 443.9
  75. 529.0 666.1 1110.6 ;
  76. ta3501 = evol manu le3501 ls3501 ;
  77. *
  78. CTrac1 = NUAGE 'COMP' 'T' 20. 350.
  79. 'COMP' 'TRAC' ta201 ta3501 ;
  80.  
  81. *--- Données des aciers 16 et 18MND5 (20 - 350°C)
  82. lT12 = prog 20. 350. ;
  83. lE2 = prog 204. 180. ;
  84. evE2 = evol manu 'T' lT12 'YOUN' (lE2*1e3) ;
  85. * dans le cas COMPLET pb avec la pente de lsi202
  86. * lep202 = prog 0. 0.001691 0.001891 0.00219 0.00269
  87. * 0.00369 0.00569 0.00769 0.00969 0.0117 0.0168 0.0219
  88. * 0.0320 0.0421 0.0522 0.0623 0.0723 0.0823 0.0924
  89. * 0.102 0.303 1.00 5.00 ;
  90. * lsi202 = prog 0. 345 345 345 345 345 345 345 345 355.4
  91. * 372.6 386.4 410.6 434.2 448.1 461.2 470.8 477.7
  92. * 484.6 491.6 570.9 671.5 834.2 ;
  93. * lep202 = prog 0. 0.001691 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 346 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.300/204.) 0.002 0.0025 0.008 0.0117
  100. 0.0168 0.0219
  101. 0.0320 0.0421 0.0522 0.0623 0.0723 0.0823 0.0924
  102. 0.102 0.303 1.00 5.00 ;
  103. lsi202 = prog 0. 300. 320. 330. 345. 355.4
  104. 372.6 386.4 410.6 434.2 448.1 461.2 470.8 477.7
  105. 484.6 491.6 570.9 671.5 834.2 ;
  106. ta202 = evol manu lep202 lsi202 ;
  107. *
  108. le3502 = prog 0. 0.001423 0.001679 0.00206 0.0026
  109. 0.00366 0.00573 0.00778 0.00983 0.0119 0.0169 0.0220
  110. 0.0321 0.0422 0.0523 0.0623 0.0724 0.0824 0.0925
  111. 0.102 0.303 1.00 5.00 ;
  112. ls3502 = prog 0. 256.2 266.1 281.0 290.0 299.0 310.9
  113. 319.9 328.9 334.9 349.8 364.8 385.7 400.2 412.3 422.9
  114. 429.0 435.0 441.0 447.0 509.0 586.7 709.3 ;
  115. ta3502 = evol manu le3502 ls3502 ;
  116. *
  117. CTrac2 = NUAGE 'COMP' 'T' 20. 350.
  118. 'COMP' 'TRAC' ta202 ta3502 ;
  119.  
  120. si(GRAPH);
  121. dess ((coul ta201 bleu) et (ta3501 coul roug)
  122. et (coul ta202 turq) et (ta3502 coul rose));
  123. dess ((coul ta201 bleu) et (ta3501 coul roug)
  124. et (coul ta202 turq) et (ta3502 coul rose))
  125. 'XBOR' 0. 0.02;
  126. fins;
  127.  
  128.  
  129. *******************************************************
  130. *** Maillage
  131.  
  132. si(COMPLET); nW = 50 ;
  133. sino; nW = 24;
  134. fins;
  135. nhF = enti (nW * (h - hint) / h) ;
  136.  
  137. P1 = 0. (-0.5 * h) ;
  138. P2 = W (-0.5 * h) ;
  139. l1 = P1 d nW p2 ;
  140.  
  141. si (hfis ega hint) ;
  142. nhA = enti (nW * hint / h) ;
  143. SA = l1 tran nhA (0. hint) ;
  144. l2 = l1 plus (0. hint) ;
  145. SF = l2 tran nhF (0. (h - hint)) ;
  146. pfis = l2 point proc ((asW * W) (hint - (0.5 * h))) ;
  147. p3 = l2 point proc (W (hint - (0.5 * h))) ;
  148. llig = l2 elem comp pfis p3 ;
  149. elim dxelim (llig et SA) ;
  150. l3 = l1 plus (0. h) ;
  151. elim dxelim (SF et l3) ;
  152. levF = diff l2 llig ;
  153. levA = levF plus (0. 0.) ;
  154. elim dxelim (SA et levA) ;
  155. sinon ;
  156. nhA = enti (nW * hfis / h) ;
  157. SA1 = l1 tran nhA (0. hfis) ;
  158. l2 = l1 plus (0. hfis) ;
  159. SA2 = l2 tran (nhF - nhA) (0. (hint - hfis)) ;
  160. l4 = l1 plus (0. hint) ;
  161. SF = l4 tran nhF (0. (h - hint)) ;
  162. pfis = l2 point proc ((asW * W) (hint - (0.5 * h))) ;
  163. p3 = l2 point proc (W (hint - (0.5 * h))) ;
  164. llig = l2 elem comp pfis p3 ;
  165. elim dxelim (llig et SA1) ;
  166. l3 = l1 plus (0. h) ;
  167. elim dxelim (SF et l3) ;
  168. levF = diff l2 llig ;
  169. levA = levF plus (0. 0.) ;
  170. elim dxelim (SA1 et levA) ;
  171. elim dxelim (SA2 et SF) ;
  172. SA = SA1 et SA2 ;
  173. finsi ;
  174.  
  175. plaq = (SA coul bleu) et (SF coul roug) ;
  176. si(GRAPH);trac plaq ;fins;
  177.  
  178.  
  179. *******************************************************
  180. *** MODELE MECANIQUE
  181.  
  182.  
  183. *--- Modèle elastique ou elasto plastique ?
  184. si calel ;
  185. mo316 = MODELE sa MECANIQUE ELASTIQUE ISOTROPE ;
  186. mo508 = MODELE sf MECANIQUE ELASTIQUE ISOTROPE ;
  187. ma316 = materiau mo316 young evE1 nu 0.3 alph 1E-5 ;
  188. ma508 = materiau mo508 young evE2 nu 0.3 alph 1E-5 ;
  189. sinon ;
  190. mo316 = MODELE sa MECANIQUE ELASTIQUE plastique ;
  191. mo508 = MODELE sf MECANIQUE ELASTIQUE plastique ;
  192. ma316 = materiau mo316 young evE1 nu 0.3 alph 1E-5 trac CTrac1 ;
  193. ma508 = materiau mo508 young evE2 nu 0.3 alph 1E-5 trac CTrac2 ;
  194. finsi ;
  195. mo = mo316 et mo508 ;
  196. ma = ma316 et ma508 ;
  197.  
  198. * ----- CONDITIONS LIMITES
  199. cd1 = bloq l1 UY ;
  200. cd2 = RELA ENSE UY l3 ;
  201. cd3 = bloq pfis UX ;
  202. cdl = cd1 et cd2 et cd3 ;
  203.  
  204. *----- Chargement
  205. lx1 = prog 0. 1. ;
  206. ly1 = prog 0. 1. ;
  207. ev1 = evol manu lx1 ly1 ;
  208. PRES1 = PRES MASS mo (-100.) l3 ;
  209. *PRES1 = PRES MASS mo dxelim l3 ;
  210. CHAP = CHAR MECA PRES1 ev1 ;
  211.  
  212. * ----- Définition de la temperature
  213. Ti = 20. ; Tf = 350. ;
  214. chT1 = manu chpo plaq 1 'T' Ti ;
  215. chx = plaq coor 1 ;
  216. chT2 = chT1 - (nomc 'T' (chx / W * (Ti - Tf))) ;
  217. TT1 = table ; TT1 . 0 = 0. ; TT1 . 1 = 1. ;
  218. TT2 = table ; TT2 . 0 = chT1 ; TT2 . 1 = chT2 ;
  219. CHAT = CHAR 'T' TT1 TT2 ;
  220.  
  221. * si (COMPLET); tcha = prog 0. 0.1 pas 0.05 1. ;
  222. si (COMPLET); tcha = prog 0. pas 0.05 1. ;
  223. sino; tcha = prog 0. pas 0.20 1. ;
  224. fins;
  225.  
  226.  
  227. *******************************************************
  228. *** PROCEDURE PAS A PAS
  229.  
  230. TAB2 = TABLE;
  231. TAB2.'MODELE' = mo ;
  232. TAB2.'CARACTERISTIQUES' = ma ;
  233. TAB2.'BLOCAGES_MECANIQUES' = cdl ;
  234. TAB2.'CHARGEMENT'= CHAP et CHAT ;
  235. TAB2.'TEMPS_CALCULES' = tcha ;
  236.  
  237. PASAPAS TAB2 ;
  238.  
  239. si(GRAPH);
  240.  
  241. ii = (dime tcha) - 1 ;
  242. def = defo tab2.deplacements.ii plaq ;
  243. trac chT2 def 'TITR' 'T et deformee finale';
  244.  
  245. si (non calel);
  246. var2 = TAB2 . 'VARIABLES_INTERNES' . ii;
  247. var2prg = (prog 0.) et (exp ((prog -3. PAS 0.1 -1.) * (log 10.)));
  248. trac var2 mo def var2prg 'TITR' 'variables_internes finales';
  249. fins;
  250.  
  251. fins;
  252.  
  253.  
  254. *******************************************************
  255. *** CALCUL DU TAUX DE RESTITUTION D'ENERGIE
  256.  
  257. *--- boucle sur les contours lnc ----------------------
  258. si(COMPLET);
  259. lnc = lect 1 2 3 5 7 8 9 11 ;
  260. lcoul = mots VIOL BLEU TURQ VERT OLIV JAUN ORAN ROUG;
  261. sino;
  262. lnc = lect 1 2 3 5 7 ;
  263. lcoul = mots BLEU TURQ VERT JAUN ROUG;
  264. fins;
  265. nlc = dime lnc ; i = 0 ;
  266. tath = table ;
  267. lJM = prog ;
  268. tath.0 = tcha ;
  269.  
  270. repeter bloc1 nlc ;
  271. i = i + 1 ;
  272. ncou = extr lnc i;
  273.  
  274. SUPTAB = TABLE ;
  275. SUPTAB.'OBJECTIF'= 'J';
  276. SUPTAB.'LEVRE_SUPERIEURE' = levF ;
  277. SUPTAB.'LEVRE_INFERIEURE' = levA;
  278. SUPTAB.'COUCHE' = ncou ;
  279. SUPTAB.'FRONT_FISSURE' = pfis ;
  280. SUPTAB.'SOLUTION_PASAPAS' = TAB2;
  281. G_THETA SUPTAB ;
  282.  
  283. dess suptab.EVOLUTION_RESULTATS
  284. 'TITR' (chai i 'couche' ncou'(G_THETA)');
  285.  
  286. mocoul1 = extr lcoul i;
  287. si(ega i 1);
  288. evotot = coul suptab.EVOLUTION_RESULTATS mocoul1;
  289. sino;
  290. evotot = evotot et (coul suptab.EVOLUTION_RESULTATS mocoul1);
  291. fins;
  292.  
  293. tath.i = extr (suptab.EVOLUTION_RESULTATS) ordo ;
  294. Jmaxi = extr tath.i (dime tath.i);
  295. lJM = lJM et Jmaxi ;
  296.  
  297. fin bloc1 ;
  298. *--- fin de boucle sur les contours lnc ----------------------
  299.  
  300.  
  301. lel = W / nW ;
  302. evJ = evol 'BLEU' manu 'dimcouche' (lel*lnc) 'J' lJM ;
  303. dess (evJ) 'TITR' (chai 'stabilite de J');
  304.  
  305. ybor0 = 0.;
  306. ybor1 = maxi (extr evJ 'ORDO');
  307. ybor1 = (1.1 * ybor1);
  308. tdess1 = tabl;
  309. tdess1 . 1 = mot 'MARQ PLUS';
  310. dess evJ 'YBOR' ybor0 ybor1 tdess1
  311. 'TITR' (chai 'stabilite de J');
  312.  
  313. dess evotot 'YBOR' ybor0 ybor1 'TITR' (chai 'evolution de J');
  314.  
  315. * opti donn 5 trac X ;
  316.  
  317.  
  318. *******************************************************
  319. *** TEST DE BON FONCTIONNEMENT
  320.  
  321. * valeurs testees
  322. Jmax1 = maxi lJM; Jmin1 = mini lJM;
  323. mess 'Jmax1 Jmin1=' Jmax1 Jmin1;
  324.  
  325. * bornes de references (au 11 aout 2011)
  326. si(COMPLET); Jmax0 = 75.; Jmin0 = 65.;
  327. sino; Jmax0 = 80.; Jmin0 = 70.;
  328. fins;
  329.  
  330. *test
  331. SI ((Jmax1 < Jmax0) et (Jmin1 > Jmin0));
  332. ERRE 0;
  333. SINON;
  334. ERRE 5;
  335. FINSI;
  336.  
  337.  
  338.  
  339. FIN ;
  340.  
  341.  
  342.  
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  

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