Télécharger waam2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier waam2.dgibi
  2. * section : thermique conduction convection mecanique plastique
  3. *----------------------------------------------------------------------*
  4. * W A A M 2 . D G I B I *
  5. *----------------------------------------------------------------------*
  6. * Objet :
  7. * -------
  8. *
  9. * Exemple de simulation thermomecanique d'un depot de matiere WAAM.
  10. * L'exemple simule la realisation d'un "mur" sur la tranche d'une tole
  11. * en acier inox (donnees 316L issues de la litterature, ref. fournie).
  12. *
  13. * On modelise la conduction, la convection avec le milieu ambiant
  14. * et la source de chaleur selon une distribution gaussienne.
  15. *
  16. * Le comportement mecanique est elastoplastique avec un ecrouissage
  17. * isotrope.
  18. *
  19. * L'apport de matiere est defini a l'aide d'un chargement de nom MODE
  20. * qui decrit l'evolution du modele au cours du temps. Un chargement de
  21. * nom MATE definit celle de ses caracteristiques.
  22. *
  23. * La modelisation fait appel aux procedures SOUDAGE et WAAM pour
  24. * definir la sequence de soudage et la mailler.
  25. *
  26. * On resout la thermique et la mecanique separemment, par deux appels
  27. * distincts a la procedure PASAPAS.
  28. *
  29. * Enfin, ce cas-test sert de validation a la prise en compte de la
  30. * FUSION du materiau dans le calcul du comportement. Pour cela, on verifie
  31. * que les def. inelastiques et variables internes sont nulles aux pts
  32. * d'integration ou T > Tfusion a l'aide de la procedure PERSO1.
  33. *
  34. * Description :
  35. * -------------
  36. *
  37. * Type de calcul : Thermique, Mecanique
  38. * Mode de calcul : 3D
  39. * Type d'element : CUB8
  40. * Chargement : Source de chaleur, Convection, Apport de matiere
  41. *
  42. *----------------------------------------------------------------------*
  43. *
  44. opti dime 3 elem cub8 ;
  45.  
  46. * icomplet = vrai : simulation de 10 passes de couches WAAM
  47. * icomplet = faux : 10 premiers pas de temps
  48. * ivalid : mettre a vrai pour verifier mise a zero EPIN et VARINT avec PERSO1
  49. icomplet = faux ;
  50. ivalid = vrai ;
  51.  
  52. * ig1 : mettre a vrai pour affichages
  53. ig1 = faux ;
  54. *opti trac psc eptr 5 ;
  55. ig1 = ig1 ou (ega (vale trac) 'PSC') ;
  56.  
  57. *----------------------- Parametres du probleme -----------------------*
  58.  
  59. * Parametres geometrie :
  60. * lw1 : longueur faite en WAAM
  61. * ls1 : longueur support
  62. * hs1 : hauteur support
  63. * e1 : epaisseur de la plaque
  64. lw1 = 80.e-3 ;
  65. ls1 = 124.e-3 ;
  66. hs1 = 50.e-3 ;
  67. e1 = 6.e-3 ;
  68.  
  69. * Position thermocouples :
  70. * Bord deport WAAM initial en (0 0 0) :
  71. TC1 = 40.e-3 0 -3.e-3 ;
  72. TC2 = 40.e-3 0 -5.e-3 ;
  73. TC3 = 0.e-3 0 -25.e-3 ;
  74. TC4 = 80.e-3 0 -10.e-3 ;
  75.  
  76. * Parametres apport de chaleur :
  77. * Us1 : tension electrique de soudage (V) 1ere couche
  78. * Us2 : tension electrique de soudage (V) couches suivantes
  79. * Is1 : intensite electrique soudage (A) 1ere couche
  80. * Is2 : intensite electrique soudage (A) couches suivantes
  81. * Vs1 : vitesse de soudage (m/s)
  82. * Eta1 : rendement de la source
  83. * R0 : rayon de la distribution Gaussienne
  84. Us1 = 15.0 ;
  85. Us2 = 13.0 ;
  86. Is1 = 120. ;
  87. Is2 = 100. ;
  88. Vs1 = 0.30 / 60. ;
  89. Eta1 = 0.8 ;
  90. R0 = 3.e-3 ;
  91.  
  92. * Parametres apport de matiere :
  93. * dfil1 : diametre file (m)
  94. * vfil1 : vitesse defilement fil (m/s) 1ere passe
  95. * vfil2 : vitesse defilement fil (m/s) autres passes
  96. dfil1 = 1.2e-3 ;
  97. vfil1 = 3.2 / 60. ;
  98. vfil2 = 2.5 / 60. ;
  99.  
  100. * Parametres sequence soudage :
  101. * dtini1 : delai initial avant debut deplacement torche
  102. * delai1 : delai entre 2 passes de WAAM
  103. * nbpass1: nombre de passes
  104. dtini1 = 0.8 ;
  105. delai1 = 30. ;
  106. nbpass1 = 10 ;
  107.  
  108. * Parametres conditions thermiques initiales et aux limites :
  109. * Tini1 : temperature initiale (degC)
  110. * Te1 : temperature de convection (degC)
  111. * h1 : coefficient d'echange convectif (W/m2)
  112. Tini1 = 20. ;
  113. Te1 = 20. ;
  114. h1 = 20. ;
  115.  
  116. *-------------------------- Proprietes 316 L --------------------------*
  117. * Ref. :
  118. * Camille Cambon, Issam Bendaoud, Sebastien Rouquette, Fabien Soulie.
  119. * "Influence of the first weld bead on strain and stress states in wire+arc additive manufacturing”.
  120. * The 12th International Seminar ”Numerical Analysis of Weldability",
  121. * Institute for Materials Science, Joining and Forming (IMAT),Sep 2018, Seggau, Austria.
  122. * hal-01954354 https://hal.archives-ouvertes.fr/hal-01954354
  123. *
  124. * rho1 : masse volumique (kg/m3)
  125. * Tfus1 : temperature de changement de phase solide-liquide (degC)
  126. * Qlat1 : chaleur latente de fusion de l'acier (J/m3)
  127. rho1 = 7760. ;
  128. Tfus1 = 1450. ;
  129. Qlat1 = 2.e9 ;
  130.  
  131. * k1 : conductivite thermique (W/m/K) en fonction de la temperature
  132. lt1 = prog 20. 100. pas 100. 1500. ;
  133. lk1 = prog 14. 15.2 16.6 17.9 19. 20.6 21.8 23.1 24.3 26. 27.3 28.6 29.9 45. 60. 60. ;
  134. k1 = evol vert manu 'T' lt1 'K' lk1 ;
  135. si ig1 ;
  136. dess k1 titr 'Conductivite thermique 316L (W/m/K)' ;
  137. fins ;
  138.  
  139. * Cp1 : capacite thermique massique (J/m3) en fonction de la temperature
  140. * On derive l'enthalpie :
  141. lt1 = prog 0. pas 100. 1500. ;
  142. lh1 = prog 0.983 1.36 1.76 2.19 2.62 3.06 3.51 3.96 4.43 4.92 5.42 5.92 6.44 6.97 9.51 ;
  143. lh2 = prog 1.36 1.76 2.19 2.62 3.06 3.51 3.96 4.43 4.92 5.42 5.92 6.44 6.97 9.51 10.1 ;
  144. dhdt1 = (lh2 - lh1) / 100. ;
  145. lc0 = dhdt1 * 1.e9 / rho1 ;
  146. lt0 = prog 50. pas 100. 1450. ;
  147. cp0 = evol manu 'T' lt0 'C' lc0 ;
  148. idx0 = posi (maxi lc0) dans lc0 ;
  149. lc2 = lc0 enle idx0 ;
  150. lt2 = lt0 enle idx0 ;
  151. cp2 = evol oran manu 'T' lt2 'C' lc2 ;
  152. cp1 = cp2 ;
  153. si ig1 ;
  154. tleg1 = table ;
  155. tleg1 . titre = table ;
  156. tleg1 . titre . 1 = mot 'Cp derivee publi.' ;
  157. tleg1 . titre . 2 = mot 'Cp utilisee' ;
  158. dess (cp0 et cp1) lege tleg1 titr ' Capacite thermique massique 316L (J/m3)' ;
  159. fins ;
  160.  
  161. * Proprietes mecaniques 316L (dependantes de T) :
  162. * nu1 : coefficient de Poisson
  163. nu1 = 0.3 ;
  164.  
  165. * Ym1 : module de Young (Pa)
  166. lt1 = prog 20. 100. pas 100. 1200. ;
  167. lym1 = prog 197000. 191500. 184000. 176500. 168000. 160000. 151500. 142500. 130000. 108000. 81500. 32000. 7400. ;
  168. lym1 = 1.e6 * lym1 ;
  169. ym1 = evol vert manu T lt1 youn lym1 ;
  170.  
  171. * Sy1 : limite d'elasticite (Pa)
  172. lt1 = prog 20. pas 100. 1400. ;
  173. lsy1 = prog 287. 237. 198. 172. 157. 151. 145. 136. 127. 115. 79. 38. 24. 18. 2. ;
  174. lsy1 = 1.e6 * lsy1 ;
  175. sy1 = evol vert manu T lt1 SIGY lsy1 ;
  176.  
  177. * Alph1 : dilatation thermique (/K)
  178. lt1 = prog 20. 100. 300. 400. 600. 700. 800. 1000. 1400. ;
  179. lalph1 = prog 1.55 1.60 1.71 1.75 1.84 1.87 1.90 1.94 1.96 ;
  180. lalph1 = 1.e-5 * lalph1 ;
  181. alph1 = evol vert manu T lt1 ALPH lalph1 ;
  182.  
  183. * Hx : module d'ecrouissage cinematique (Pa)
  184. lt1 = prog 20. 700. pas 100. 1300. ;
  185. lh1 = prog 2400. 2400. 2350. 1500. 800. 725. 150. 10. ;
  186. lh1 = 1.e6 * lh1 ;
  187. Hx1 = evol vert manu T lt1 H lh1 ;
  188.  
  189. * Courbes d'ecrouissage :
  190. nb1 = dime lt1 ;
  191. lep1 = prog 0. 1. ;
  192. nua1 = vide nuage ;
  193. ecro1 = vide evolution ;
  194. tleg1 = table ;
  195. tleg1 . titre = table ;
  196. repe b1 nb1 ;
  197. ti1 = extr lt1 &b1 ;
  198. syi1 = ipol sy1 ti1 ;
  199. hi1 = ipol Hx1 ti1 ;
  200. lmsi1 = prog syi1 (syi1+hi1) ;
  201. ecroi1 = evol oran manu eps lep1 sig lmsi1 ;
  202. ecroi1 = ecroi1 coul &b1 ;
  203. nua1 = nua1 et (nuag comp T ti1 comp ECRO ecroi1) ;
  204. ecro1 = ecro1 et ecroi1 ;
  205. tleg1.titre.&b1 = chai 'T = ' (enti ti1) ;
  206. fin b1 ;
  207.  
  208. si ig1 ;
  209. dess ym1 titr ' Module de Young (Pa)' ;
  210. dess sy1 titr ' Limite elastique (Pa)' ;
  211. dess alph1 titr ' Coef. dilatation thermique (/K)' ;
  212. dess (1.e-6*ecro1) xbor 0. 0.2 ybor 0. 500.
  213. titr 'Courbes d ecrouissage a differentes temperatures (MPa)' lege tleg1 ;
  214. fins ;
  215.  
  216. *------------------- Sequence de soudage / Maillage -------------------*
  217.  
  218. * Parametres discretisation apport de matiere :
  219. * xp1 : "pas" de discetisation de l'apport de matiere en espace (m)
  220. * dz1 : increment WAAM selon (O,z) / hauteur cordon
  221. * debi1 : debit volumique de fil 1ere passe
  222. * debi2 : debit volumique de fil autres passes
  223. xp1 = e1 / 2. ;
  224. debi1 = 0.25 * dfil1 * dfil1 * pi * vfil1 ;
  225. debi2 = 0.25 * dfil1 * dfil1 * pi * vfil2 ;
  226. dz1 = debi1 / e1 / Vs1 ;
  227. dz2 = debi2 / e1 / Vs1 ;
  228.  
  229. * Parametres modele source Gaussienne 3D :
  230. Qtot1 = Eta1 * Us1 * Is1 ;
  231. Qtot2 = Eta1 * Us2 * Is2 ;
  232. *Qtot1 = Qtot1 - (Qlat1 * e1 * dz1 * Vs1) ;
  233. *Qtot2 = Qtot2 - (Qlat1 * e1 * dz2 * Vs1) ;
  234. *list Qtot1 ;
  235. *list Qtot2 ;
  236. Rg1 = ((2. / 3.) ** 0.5) * R0 ;
  237.  
  238. * Parametre geometrie / trajectoire :
  239. e1s2 = 0.5 * e1 ;
  240. dz1s2 = 0.5 * dz1 ;
  241.  
  242. * table SOUDAGE :
  243. tso1 = tabl ;
  244. tso1.vitesse_de_soudage = Vs1 ;
  245. tso1.puissance_de_soudage = Qtot2 ;
  246. tso1.diametre_de_fil = dfil1 ;
  247. tso1.vitesse_de_fil = vfil2 ;
  248. tso1.point_de_depart = (0 e1s2 dz1) ;
  249. tso1.largeur_de_passe = e1 ;
  250.  
  251. * Défintion de la sequence de soudage :
  252. * Repetition sequence de 2 passes en AR :
  253. si (nbpass1 mult 2) ;
  254. nb1 = nbpass1 / 2 ;
  255. sino ;
  256. nb1 = (nbpass1 - 1) / 2 ;
  257. fins ;
  258. soudage tso1 point dtini1 puis Qtot1 debi debi1 ;
  259. repe b1 nb1 ;
  260. si (&b1 ega 1) ;
  261. soudage tso1 passe droi (lw1 0 0) puis Qtot1 debi debi1 ;
  262. sino ;
  263. soudage tso1 passe droi (lw1 0 0) ;
  264. fins ;
  265. soudage tso1 depla couche pause delai1 ;
  266. soudage tso1 passe droi ((-1.*lw1) 0 0) ;
  267. soudage tso1 depla couche pause delai1 ;
  268. fin b1 ;
  269. si ((nbpass1-1) mult 2) ;
  270. soudage tso1 passe droi (lw1 0 0) ;
  271. fins ;
  272. soudage tso1 point 60. puis 0. debi 0. ;
  273.  
  274. si ig1 ;
  275. trac tso1.trajectoire titr 'trajectoire depot WAAM mur' ;
  276. dess tso1.evolution_deplacement titr 'evolution deplacement depot WAAM mur' ;
  277. dess tso1.evolution_puissance titr 'evolution puissance thermique depot WAAM mur' ;
  278. dess tso1.evolution_debit titr 'evolution debit apport de matiere depot WAAM mur' ;
  279. fins ;
  280.  
  281. * Maillage du Mur :
  282. * On teste la donnee d'un LISTREEL pour la discretisation en espcace
  283. ne1 = 5 ;
  284. tab2 = waam tso1 mail pas (prog 2.e-3 xp1 4.e-3 xp1 2.e-3 xp1) larg e1 dens 1.e-3 'TEMP' 0.2 ;
  285. elim tab2.maillage 1.e-6 ;
  286. mur1 = tab2.maillage coul vert ;
  287.  
  288. si ig1 ;
  289. trac cach mur1 titr ' Maillage global du mur' ;
  290. tmai1 = tab2.evolution_maillage.maillage ;
  291. waam tab2 visu cach ((tso1.trajectoire) plus (0 0 1.e-5)) ;
  292. fins ;
  293.  
  294. * Maillage du support sous-jacent :
  295. me1 = -1. * e1 ;
  296. mdz1 = -1. * dz1 ;
  297. zmin1 = (mur1 coor 3) mini ;
  298. pz0 = (mur1 coor 3) poin infe (zmin1 + 1.e-5) ;
  299. sz0 = (enve mur1) elem appu stri pz0 ;
  300. sup1 = sz0 volu tran ne1 (0 0 me1) ;
  301.  
  302. xmin1 = (sup1 coor 1) mini ;
  303. px0 = (sup1 coor 1) poin infe (xmin1 + 1.e-5) ;
  304. px1 = (sup1 coor 1) poin supe (xmin1 + lw1 - 1.e-5) ;
  305. sx0 = (enve sup1) elem appu stri px0 ;
  306. sx1 = (enve sup1) elem appu stri px1 ;
  307. sup1 = sup1 et (sx0 volu tran ne1 (me1 0 0)) et (sx1 volu tran ne1 (e1 0 0)) ;
  308. sup1 = sup1 coul gris ;
  309.  
  310. dx1 = (ls1 - lw1) * 0.5 - e1 ;
  311. me1 = -1. * e1 ;
  312. xmin2 = xmin1 + me1 ;
  313. zmin2 = zmin1 + me1 ;
  314. p1 = sup1 poin proc (xmin2 0 zmin1) ;
  315. p2 = sup1 poin proc (xmin2 0 zmin2) ;
  316. p3 = sup1 poin proc ((xmin1 + lw1 + e1) 0 zmin2) ;
  317. p4 = sup1 poin proc ((xmin1 + lw1 + e1) 0 zmin1) ;
  318. p5 = p4 plus (dx1 0 0) ;
  319. p6 = p5 moin (0 0 hs1) ;
  320. p7 = p6 moin (ls1 0 0) ;
  321. p8 = p7 plus (0 0 hs1) ;
  322. de1 = 0.5 * e1 ;
  323. de2 = 6.2e-3 ;
  324. cs2 = (p5 droi p6 dini de1 dfin de2) et
  325. (p6 droi p7 dini de2 dfin de2) et
  326. (p7 droi p8 dini de2 dfin de1) et
  327. (p8 droi p1 dini de1 dfin de1) et
  328. (p1 droi p2 dini de1 dfin de1) et
  329. (p2 droi p3 dini de1 dfin de1) et
  330. (p3 droi p4 dini de1 dfin de1) et
  331. (p4 droi p5 dini de1 dfin de1) ;
  332. fs2 = surf cs2 plan ;
  333. sup2 = fs2 volu tran (0 e1 0) 2 ;
  334. sup2 = sup2 coul turq ;
  335.  
  336. sup0 = sup1 et sup2 ;
  337.  
  338. mail1 = sup0 et mur1 ;
  339. si ig1 ;
  340. mot1 = chai 'Maillage total "Mur WAAM" et support :' (nbno mail1) ;
  341. mot1 = chai mot1 ' /' (nbel mail1) ' noeuds/elem.' ;
  342. trac face mail1 titr mot1 ;
  343. fins ;
  344.  
  345. *------------- Accrochage DDL thermique maillage support --------------*
  346.  
  347. * Relation cinematique DDL mecanique entre maillages non conformes :
  348. clu1 = sup1 rela accro sup2 ;
  349.  
  350. * Encastrement bord plaque support :
  351. ptx0 = (sup0 coor 1) poin infe ((sup0 coor 1 mini) + 1.e-5) ;
  352. clu2 = bloq depl ptx0 ;
  353.  
  354. * Relation cinematique DDL thermique entre maillages non conformes :
  355. clt1 = sup1 rela accro sup2 (mots 'T') ;
  356. clt2 = (TC1 et TC2 et TC3 et TC4) rela accro sup0 (mots 'T') ;
  357.  
  358. * Surface interface maillages support non-conformes :
  359. pinte1 = clt1 extr mail nomu ;
  360. sinte1 = (enve sup0) elem appu stri pinte1 ;
  361.  
  362. si ig1 ;
  363. trac cach (-1 -1 0.5) ((ptx0 coul roug) et mail1) titr 'Points encastrement support (rouge)' ;
  364. fins ;
  365.  
  366. *--------------------- Modele / Caracteristique -----------------------*
  367.  
  368. * Conduction / Convection :
  369. mod1 = mode mail1 thermique ;
  370. mat1 = mate mod1 rho rho1 k k1 'C' cp1 'TINI' Tini1 ;
  371.  
  372. * Source thermique :
  373. evqt1 = tso1.evolution_puissance ;
  374. evxs1 = tso1.evolution_deplacement ;
  375. chxs1 = tso1.trajectoire coor curv ;
  376. cgxs1 = char traj chxs1 evxs1 ;
  377. mod3 = mode mail1 thermique source gaussienne ;
  378. mat3 = mate mod3 qtot evqt1 orig cgxs1 rgau Rg1 ;
  379.  
  380. * Mecanique :
  381. mod4 = mode mail1 mecanique elastique plastique isotrope fusion cons 316L ;
  382. mat4 = mate mod4 youn ym1 nu nu1 alph alph1 ecro ecro1 tref tini1 talp tini1 tfus Tfus1 ;
  383.  
  384. * Tables d'evolution des modeles et caracteristiques :
  385. ttps1 = tab2.evolution_maillage.temps ;
  386. tmai1 = tab2.evolution_maillage.maillage ;
  387. nb1 = dime ttps1 ;
  388. tmod1 = table ;
  389. tmod2 = table ;
  390. tmod3 = table ;
  391. tmod4 = table ;
  392. tmat1 = table ;
  393. tmat2 = table ;
  394. tmat3 = table ;
  395. tmat4 = table ;
  396. repe b1 nb1 ;
  397. geoi1 = tmai1 . (&b1 - 1) ;
  398. geoi1 = geoi1 et sup0 ;
  399. tmod1 . (&b1 - 1) = redu mod1 geoi1 ;
  400. tmat1 . (&b1 - 1) = redu mat1 (tmod1 . (&b1 - 1)) ;
  401. sconv1 = (enve geoi1) diff sinte1 ;
  402. tmod2 . (&b1 - 1) = mode sconv1 thermique convection ;
  403. tmat2 . (&b1 - 1) = mate (tmod2 . (&b1 - 1)) 'H' h1 'TE' Te1 ;
  404. tmod3 . (&b1 - 1) = redu mod3 geoi1 ;
  405. tmat3 . (&b1 - 1) = redu mat3 (tmod3 . (&b1 - 1)) ;
  406. tmod4 . (&b1 - 1) = redu mod4 geoi1 ;
  407. tmat4 . (&b1 - 1) = redu mat4 (tmod4 . (&b1 - 1)) ;
  408. fin b1 ;
  409.  
  410. * Chargements MODE / MATE :
  411. cgmod1 = char mode ttps1 tmod1 ;
  412. cgmod2 = char mode ttps1 tmod2 ;
  413. cgmod3 = char mode ttps1 tmod3 ;
  414. cgmod4 = char mode ttps1 tmod4 ;
  415. cgmod0 = cgmod1 et cgmod2 et cgmod3 ;
  416.  
  417. cgmat1 = char mate ttps1 tmat1 ;
  418. cgmat2 = char mate ttps1 tmat2 ;
  419. cgmat3 = char mate ttps1 tmat3 ;
  420. cgmat4 = char mate ttps1 tmat4 ;
  421. cgmat0 = cgmat1 et cgmat2 et cgmat3 ;
  422.  
  423. cg0 = cgmod0 et cgmat0 ;
  424.  
  425. *----------------------------------------------------------------------*
  426. * THERMIQUE *
  427. *----------------------------------------------------------------------*
  428. *
  429. * Liste des temps de calcul :
  430. ltca1 = tab2.temps_calcules ;
  431. si (non icomplet) ;
  432. ltca1 = ltca1 extr (lect 1 pas 1 10) ;
  433. fins ;
  434.  
  435. *-------------------- Resolution Thermique PASAPAS --------------------*
  436.  
  437. * Table PASAPAS :
  438. tab1 = table ;
  439. tab1.modele = tire cgmod0 mode 0. ;
  440. tab1.caracteristiques = tire cgmat0 mate 0. ;
  441. tab1.chargement = cg0 ;
  442. tab1.blocages_thermiques = clt1 et clt2 ;
  443. tab1.temps_calcules = ltca1 ;
  444.  
  445.  
  446. pasapas tab1 ;
  447.  
  448. * Affichage temps de calcul :
  449. duree1 = temp enti horl ;
  450. duree1 = (duree1 / 1000) ;
  451. dmin1 = (duree1 / 60) ;
  452. dsec1 = duree1 - (60 * dmin1) ;
  453. mot1 = chai '***** DUREE DU CALCUL THERMIQUE (s) :' dmin1 ' min' dsec1 ' s' ;
  454. mess mot1 ;
  455.  
  456. *----------------------------------------------------------------------*
  457. * MECANIQUE *
  458. *----------------------------------------------------------------------*
  459. *
  460. *-------------------------- Procedure PERSO1 --------------------------*
  461. * Verification mise a zero variables internes et def. plastique si T>Tfusion :
  462. debp perso1 tu1*table ;
  463.  
  464. si (tu1.perso1_appel ega 1) ;
  465. quit perso1 ;
  466. fins ;
  467.  
  468. tfus1 = tu1.data_perso1.tfusio ;
  469. cht2 = (tu1.wtable.tet2) exco 'T' 'T' ;
  470. mod1 = tu1.wtable.modele extr form mecanique ;
  471. chamt2 = chan cham mod1 cht2 stresses ;
  472. si ((maxi chamt2) >EG tfus1) ;
  473. mess ' ****** PERSO1 : verif. mise a zero Epin & Var_Int' ;
  474. masq1 = chamt2 masq supe tfus1 ;
  475. masq1 = masq1 nomc 'SCAL' ;
  476. *trac (0 1 1) 13 masq1 mod1 titr ' Elem T < tfus.' ;
  477. *
  478. * Def. Inelastiques :
  479. epi0 = tu1.estimation.deformations_inelastiques redu mod1 ;
  480. lmc0 = extr epi0 comp ;
  481. lmc1 = mots (dime lmc0) * 'SCAL' ;
  482. epi1 = masq1 * epi0 lmc1 lmc0 lmc0 ;
  483. err1 = maxi abs epi1 ;
  484. *trac 13 epi1 mod1 titr ' MAJ def. inelas.' ;
  485. *
  486. * Var. Internes :
  487. var0 = tu1.estimation.variables_internes redu mod1 ;
  488. lmc0 = extr var0 comp ;
  489. lmc1 = mots (dime lmc0) * 'SCAL' ;
  490. var1 = masq1 * var0 lmc1 lmc0 lmc0 ;
  491. err1 = (maxi abs var1) + err1 ;
  492. si (err1 > ((vale prec)*10.)) ; erre 308 ; fins ;
  493. *trac (0 1 1) 13 var1 mod1 titr ' MAJ var. inter.' ;
  494. fins ;
  495.  
  496. finp ;
  497. *-------------------- Resolution Mecanique PASAPAS --------------------*
  498.  
  499. * Chargement thermique :
  500. cgt1 = char 'T' (tab1.temps) (tab1.temperatures) ;
  501.  
  502. * Table PASAPAS :
  503. tab2 = table ;
  504. tab2.modele = tire cgmod4 mode 0. ;
  505. tab2.caracteristiques = tire cgmat4 mate 0. ;
  506. tab2.chargement = cgmod4 et cgmat4 et cgt1 ;
  507. tab2.blocages_mecaniques = clu1 et clu2 ;
  508. tab2.temps_calcules = ltca1 ;
  509. tab2.processeurs = mot comportement ;
  510.  
  511. si icomplet ;
  512. ntca1 = dime ltca1 ;
  513. tab2.temps_sauvegardes = ltca1 extr (lect 50 pas 50 (ntca1 / 50 * 50) ntca1) ;
  514. opti sauv 'waam2.sauv' ;
  515. fins ;
  516.  
  517. tab2.procedure_perso1 = ivalid ;
  518. tab2.data_perso1 = table ;
  519. tab2.data_perso1.tfusio = tfus1 ;
  520.  
  521.  
  522. pasapas tab2 ;
  523.  
  524. * Affichage temps de calcul :
  525. duree1 = temp enti horl ;
  526. duree1 = (duree1 / 1000) ;
  527. dmin1 = (duree1 / 60) ;
  528. dsec1 = duree1 - (60 * dmin1) ;
  529. mot1 = chai '***** DUREE DU CALCUL MECANIQUE (s) :' dmin1 ' min' dsec1 ' s' ;
  530. mess mot1 ;
  531.  
  532. *------------------------ Petit post-traitement -----------------------*
  533.  
  534. * Animation contrainte de von Mises :
  535. vtra1 = vale trac ;
  536.  
  537. si ig1 ;
  538. if1 = (dime tab2.temps) - 1 ;
  539. def1 = vide deforme ;
  540. si (ega vtra1 'PSC') ;
  541. opti ftra 'Contrainte_Waam2_Anim.ps' ;
  542. fins ;
  543. liso1 = prog 25. pas 25. 350. ;
  544. mbox1 = boite mail1 ;
  545. repe b1 if1 ;
  546. i1 = &b1-1 ;
  547. modi1 = tab2.modeles.i1 ;
  548. maili1 = (extr modi1 form mecanique) extr mail ;
  549. depi1 = tab2.deplacements.i1 ;
  550. sigi1 = tab2.contraintes.i1 ;
  551. sigi1 = 1.e-6 * sigi1 ;
  552. sigi1 = vmis modi1 sigi1 ;
  553. defoi1 = defo maili1 depi1 1. sigi1 modi1 ;
  554. def1 = def1 et defoi1 ;
  555. si (ega vtra1 'PSC') ;
  556. mot1 = chai format '(F6.1)' 'Contrainte de von Mises (MPa) au temps (s) :' (tab2.temps.i1) ;
  557. trac (1 -1.2 1) defoi1 liso1 titr mot1 boit mbox1 ;
  558. fins ;
  559. fin b1 ;
  560. si (neg vtra1 'PSC') ;
  561. trac (1 -1.2 1) anim def1 liso1 ;
  562. fins ;
  563. fins ;
  564.  
  565. si icomplet ; sauv ; fins ;
  566.  
  567. *------------------ F I N W A A M 2 . D G I B I -----------------*
  568. fin ;
  569.  
  570.  
  571.  
  572.  
  573.  
  574.  
  575.  
  576.  
  577.  
  578.  
  579.  

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