Télécharger soudage6.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier soudage6.dgibi
  2. * section : thermique conduction convection
  3. *----------------------------------------------------------------------*
  4. * S O U D A G E 6 . D G I B I *
  5. *----------------------------------------------------------------------*
  6. * Objet :
  7. * -------
  8. *
  9. * Exemple de simulation thermique du soudage d'un raidisseur sur une
  10. * plaque avec apport de matiere (4 passes). Les materiaux sont en acier
  11. * inoxydable (316 L).
  12.  
  13. * On modelise la conduction, la convection avec le milieu ambiant
  14. * et la source de chaleur selon une distribution gaussienne.
  15. *
  16. * L'apport de matiere est defini a l'aide d'un chargement de nom MODE
  17. * qui decrit l'evolution du modele au cours du temps. Un chargement de
  18. * nom MATE definit celle de ses caracteristiques.
  19. *
  20. * La procedure SOUDAGE permet de definir la sequence de soudage et
  21. * d'indexer les elements du maillage des cordons de soudure en fonction
  22. * de l'evolution de l'outil le long de la trajectoire de soudage.
  23. *
  24. * Description :
  25. * -------------
  26. *
  27. * Type de calcul : Thermique Transitoire
  28. * Mode de calcul : 3D
  29. * Type d'element : CUB8, PRY6
  30. * Chargement : Source de chaleur, Convection, Apport de matiere
  31. *
  32. *----------------------------------------------------------------------*
  33. *
  34. opti dime 3 elem cub8 mode trid ;
  35.  
  36. * Pour affichage traces, mettre IG1 a VRAI :
  37. IG1 = faux ;
  38.  
  39. * Pour calcul complet, mettre icomplet a VRAI :
  40. icomplet = faux ;
  41.  
  42. *opti trac psc ;
  43. ig1 = ig1 ou (ega (vale trac) 'PSC') ;
  44.  
  45. *----------------------------------------------------------------------*
  46. * Donnees du probleme *
  47. *----------------------------------------------------------------------*
  48. *
  49. *----------------------- Parametres geometriques ----------------------*
  50.  
  51. * Geometrie pieces a souder :
  52. * la1 : largeur plaque (m)
  53. * hr1 : hauteur du raidisseur (m)
  54. * lw1 : longueur a souder (m)
  55. * ep1 : epaisseur plaque (m)
  56. * er1 : epaisseur raidisseur (m)
  57. la1 = 60.e-3 ;
  58. hr1 = 40.e-3 ;
  59. lw1 = 30.e-3 ;
  60. ep1 = 30.e-3 ;
  61. er1 = 10.e-3 ;
  62. de1 = 1.5e-3 ;
  63. de2 = 5.e-3 ;
  64. de3 = 15.e-3 ;
  65.  
  66. * Geometrie cordons :
  67. * lc1 : largeur cordons (m)
  68. lc1 = 6.e-3 ;
  69. * Pour la geometrie modelisee, section de cordon identique :
  70. Sc1 = 0.5 * lc1 * lc1 ;
  71. * Positions cordons :
  72. Pc1 = (0.5*er1*(1 0 0)) plus ((lc1 lc1 0.) * 0.5) ;
  73. Pc2 = Pc1 plus (lc1 0. 0.) ;
  74. Pc3 = Pc1 plus ((lc1 lc1 0.) * 0.5) ;
  75. Pc4 = Pc1 plus (0. lc1 0.) ;
  76.  
  77. *--------------------- Donnees probleme thermique ---------------------*
  78.  
  79. * CI / CL :
  80. * Tini1 : temperature initiale :
  81. * Te1 : temperature milieu ambiant :
  82. * h1 : coefficient d'echange convectif (W/m2/K)
  83. * em1 : emissivite
  84. * tps0 : temps de maintien avant debut deplacement torche (s)
  85. Tini1 = 20. ;
  86. Te1 = 20. ;
  87. h1 = 20. ;
  88. em1 = 0.5 ;
  89. tps0 = 1. ;
  90.  
  91. * Parametres soudage :
  92. * Us1 : tension de soudage (V)
  93. * Is1 : Intensite de soudage (A)
  94. * Vs1 : Vitesse de soudage (m/s)
  95. * Vd1 : Vitesse de deplacement (m/s)
  96. * Vf1 : vitesse de debit de fil (m/s)
  97. * Df1 : diametre du fil (m) : pas utilisee dans le calcul
  98. * Xpas1 : pas de discretisation de l'apport de matiere (m)
  99. Us1 = 12. ;
  100. Is1 = 400. ;
  101. Vs1 = 0.3 / 60. ;
  102. Vd1 = 0.1 * Vs1 ;
  103. Df1 = 1.2e-3 ;
  104. Xpas1 = 5.e-3 ;
  105.  
  106. * L'intensite de soudage est deduite du debit volumique de matiere et
  107. * de l'energie thermique qu'il faut fournir pour la fondre (on ne tient
  108. * pas compte du rendement).
  109. * A titre indicatif, calcul effectue pour l'estimer (avec : rho = 7.6e3 kg/m3,
  110. * Cp = 500. J/kg/m3, dT = 1500. degC, Qlat = 1.e9 J/m3)
  111. Is0 = Sc1 * Vs1 * (7.6e3 * 500. * 1500. + 1.e9) / Us1 ;
  112. list Is0 ;
  113.  
  114. * La vitesse de fil est deduite de la section des cordons et de la vitesse
  115. * de soudage. Elle n'intervient pas dans le calcul car geometrie cordon fixee.
  116. * A titre indicatif, son calcul et affichage (m/min.) :
  117. Vf1 = Sc1 * Vs1 / (pi * Df1 * Df1 /4.) ;
  118. List (Vf1 * 60.) ;
  119.  
  120. * Parametres modele source Gaussienne 3D :
  121. * Qtot1 : puissance thermique de la source (J)
  122. * Rg1 : "rayon" de la source Gaussienne (m)
  123. * Zg1 : "profondeur" de la source Gaussienne (m)
  124. Qtot1 = Us1 * Is1 ;
  125. Rg1 = 6.0e-3 ;
  126. Zg1 = 4.5e-3 ;
  127.  
  128. *-------------------------- Donnees materiau --------------------------*
  129. * Ref. :
  130. * Camille Cambon, Issam Bendaoud, Sebastien Rouquette, Fabien Soulie.
  131. * “Influence of the first weld bead on strain and stress states in wire+arc additive manufacturing”.
  132. * The 12th International Seminar ”Numerical Analysis of Weldability”,
  133. * Institute for Materials Science, Joining and Forming (IMAT),Sep 2018, Seggau, Austria.
  134. * hal-01954354 https://hal.archives-ouvertes.fr/hal-01954354
  135. *
  136. * rho1 : masse volumique (kg/m3)
  137. * Tfus1 : temperature de changement de phase solide-liquide (degC)
  138. * Qlat1 : chaleur latente de fusion de l'acier (J/m3)
  139. rho1 = 7760. ;
  140. Tfus1 = 1450. ;
  141. Qlat1 = 2.e9 ;
  142.  
  143. * k1 : conductivite thermique (W/m/K) en fonction de la temperature
  144. lt1 = prog 20. 100. pas 100. 1500. ;
  145. 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. ;
  146. k1 = evol vert manu 'T' lt1 'K' lk1 ;
  147. si ig1 ;
  148. dess k1 titr 'Conductivite thermique 316L (W/m/K)' ;
  149. fins ;
  150.  
  151. * Cp1 : capacite thermique massique (J/m3) en fonction de la temperature
  152. * On derive l'enthalpie :
  153. lt1 = prog 0. pas 100. 1500. ;
  154. 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 ;
  155. 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 ;
  156. dhdt1 = (lh2 - lh1) / 100. ;
  157. lc0 = dhdt1 * 1.e9 / rho1 ;
  158. lt0 = prog 50. pas 100. 1450. ;
  159. cp0 = evol manu 'T' lt0 'C' lc0 ;
  160. idx0 = posi (maxi lc0) dans lc0 ;
  161. lc2 = lc0 enle idx0 ;
  162. lt2 = lt0 enle idx0 ;
  163. cp2 = evol oran manu 'T' lt2 'C' lc2 ;
  164. cp1 = cp2 ;
  165. si ig1 ;
  166. tleg1 = table ;
  167. tleg1 . titre = table ;
  168. tleg1 . 1 = 'Cp derivee publi.' ;
  169. tleg1 . 2 = 'Cp utilisee' ;
  170. dess (cp0 et cp1) titr ' Capacite thermique massique 316L (J/m3)' lege tleg1 ;;
  171. fins ;
  172.  
  173. *----------------------------------------------------------------------*
  174. * Maillage *
  175. *----------------------------------------------------------------------*
  176.  
  177. * Repere :
  178. O1 = 0 0 0 ;
  179. X1 = 1 0 0 ;
  180. Y1 = 0 1 0 ;
  181. Z1 = 0 0 -1 ;
  182. oeil0 = 100 +100 100 ;
  183.  
  184. *-------------------- Maillage Plaque et Raidisseur -------------------*
  185.  
  186. * Plaque : maillage surface dans plan (0,xy)
  187. Pp1 = 0.5*er1*X1 ;
  188. Pp2 = Pp1 plus (3.*lc1*X1) ;
  189. lp1 = ((0 0 0) droi Pp1 dini de1 dfin de1) et (Pp1 droi Pp2 dini de1 dfin de1) ;
  190. sp1 = lp1 tran (-0.4*ep1*Y1) dini de1 dfin de1 ;
  191. sp1 = sp1 coul gris ;
  192. lp2 = sp1 cote 2 ;
  193. lp3 = sp1 cote 3 ;
  194. Pp3 = lp2 poin proc (Pp2 moin Y1) ;
  195. PP4 = Pp3 moin (0.6*ep1*Y1) ;
  196. Pp5 = PP4 moin (0.5*er1*X1) moin (3.*lc1*X1) ;
  197. Pp6 = lp3 poin proc (-1 -1 0) ;
  198. CntSp2 = Pp3 droi Pp4 dini de1 dfin de2 ;
  199. CntSp2 = CntSp2 et (Pp4 droi Pp5 dini de2 dfin de2) ;
  200. CntSp2 = CntSp2 et (Pp5 droi Pp6 dini de2 dfin de1) ;
  201. CntSp2 = CntSp2 et (inve lp3) ;
  202. Sp2 = surf CntSp2 plan ;
  203. Sp2 = sp2 coul gris ;
  204. lp22 = CntSp2 elem comp Pp3 Pp4 ;
  205. Pp7 = la1*X1 ;
  206. Pp8 = Pp7 moin (ep1*Y1) ;
  207. Pp9 = Lp22 poin proc (Pp2 moin (ep1*Y1)) ;
  208. CntSp3 = Pp2 droi Pp7 dini de1 dfin de3 ;
  209. CntSp3 = CntSp3 et (Pp7 droi Pp8 dini de3 dfin de3) ;
  210. CntSp3 = CntSp3 et (Pp8 droi Pp9 dini de3 dfin de2) ;
  211. CntSp3 = CntSp3 et (inve (Lp2 et lp22)) ;
  212. Sp3 = surf CntSp3 plan ;
  213. Sp3 = Sp3 coul gris ;
  214. spz0 = sp1 et Sp2 et Sp3 ;
  215.  
  216. * Raidisseur : maillage surface dans plan (0,xy)
  217. * Pr1 = Pp1 : Noeuds l'arete a l'angle de la plaque et du raidisseur commun aux 2 pieces :
  218. Pr2 = Pp1 plus (2.*lc1*Y1) ;
  219. lr1 = (0 0 0) droi Pp1 dini de1 dfin de1 ;
  220. sr1 = lr1 tran (2.*lc1*Y1) dini de1 dfin de1 ;
  221. Pr2 = sr1 poin proc Pr2 ;
  222. sr1 = sr1 coul vert ;
  223. lr3 = sr1 cote 3 ;
  224. sr2 = lr3 tran ((hr1-lc1-lc1)*Y1) dini de1 dfin de2 ;
  225. sr2 = sr2 coul vert ;
  226. srz0 = sr1 et sr2 ;
  227.  
  228. * Extrusion maillage volumique :
  229. Sz0 = Spz0 et Srz0 ;
  230. Vp0 = spz0 volu tran (lw1*Z1) dini de1 dfin de1 ;
  231. Vr0 = srz0 volu tran (lw1*Z1) dini de1 dfin de1 ;
  232. Vs0 = Vp0 et Vr0 ;
  233.  
  234. *------------------------ Maillage des cordons ------------------------*
  235.  
  236. * Maillage surface cordons dans plan (O,x,y) :
  237. opti elem tet4 ;
  238.  
  239. * Cordon 1 :
  240. Pc12 = lp1 poin proc (Pp1 plus (lc1*X1)) ;
  241. Pc13 = (sr1 cote 2) poin proc (Pp1 plus (lc1*Y1)) ;
  242. CntC1 = lp1 elem comp Pp1 Pc12 ;
  243. CntC1 = CntC1 et (Pc12 droi Pc13 dini de1 dfin de1) ;
  244. CntC1 = CntC1 et ((sr1 cote 2) elem comp Pc13 Pp1) ;
  245. SuC1 = surf CntC1 plan ;
  246. Suc1 = Suc1 coul 1 ;
  247.  
  248. * Cordon 2 :
  249. Pc22 = lp1 poin proc (Pp1 plus (2.*lc1*X1)) ;
  250. Pc23 = Pc22 plus ((-0.5*lc1) (0.5*lc1) 0.) ;
  251. Pc24 = Pc23 moin (lc1*X1) ;
  252. Pc24 = Cntc1 poin proc Pc24 ;
  253. CntC2 = lp1 elem comp Pc12 Pc22 ;
  254. CntC2 = CntC2 et (Pc22 droi Pc23 dini de1 dfin de1) ;
  255. CntC2 = CntC2 et (Pc23 droi Pc24 dini de1 dfin de1) ;
  256. CntC2 = CntC2 et (((CntC1 inve) elem comp Pc24 Pc12) ) ;
  257. Suc2 = Cntc2 surf plan ;
  258. Suc2 = Suc2 coul 2 ;
  259.  
  260. * Cordon 3 :
  261. Pc33 = Pc23 plus ((-0.5*lc1) (0.5*lc1) 0.) ;
  262. CntC3 = (inve CntC2) elem comp Pc24 Pc23 ;
  263. CntC3 = CntC3 et (Pc23 droi Pc33 dini de1 dfin de1) ;
  264. CntC3 = CntC3 et (Pc33 droi Pc13 dini de1 dfin de1) ;
  265. CntC3 = CntC3 et ((inve CntC1) elem comp Pc13 Pc24) ;
  266. SuC3 = surf CntC3 plan ;
  267. SuC3 = SuC3 coul 3 ;
  268.  
  269. * Cordon 4 :
  270. CntC4 = (inve CntC3) elem comp Pc13 Pc33 ;
  271. CntC4 = CntC4 et (Pc33 droi Pr2 dini de1 dfin de1) ;
  272. CntC4 = CntC4 et ((sr1 cote 2) elem comp Pr2 Pc13) ;
  273. SuC4 = CntC4 surf plan ;
  274. SuC4 = SuC4 coul 4 ;
  275.  
  276. * Volume cordons :
  277. Suc0 = Suc1 et Suc2 et Suc3 et Suc4 ;
  278. *trac (0 0 1000) Suc0 ;
  279. Vc0 = Suc0 volu tran (lw1*Z1) dini de1 dfin de1 ;
  280.  
  281. * Maillage volumique de chaque cordon :
  282. Vc1 = Vc0 elem coul 1 ;
  283. Vc2 = Vc0 elem coul 2 ;
  284. Vc3 = Vc0 elem coul 3 ;
  285. Vc4 = Vc0 elem coul 4 ;
  286.  
  287. * Fusion noeuds communs cordons/plaque et cordons/raidisseur mais pas plaque/raidisseur
  288. elim ((enve Vc0) et (enve Vp0)) 1.e-5 ;
  289. elim ((enve Vc0) et (enve Vr0)) 1.e-5 ;
  290.  
  291. * Maillage total :
  292. V0 = Vc0 et Vs0 ;
  293.  
  294. * Affichages :
  295. si IG1 ;
  296. mot1 = chai ' Maillage plaque et raidisseur :' (nbno Vs0) ' noeuds,' (nbel Vs0) ' elem.' ;
  297. trac oeil0 face Vs0 titr mot1 ;
  298. mbox2 = boite (Sp1 et Sr1 et Vc0) ;
  299. mot1 = chai ' Maillage plaque, raidisseur et cordons :' (nbno V0) ' noeuds,' (nbel V0) ' elem.' ;
  300. trac oeil0 face V0 titr mot1 ;
  301. trac (0 1 10) face V0 boit mbox2 titr mot1 ;
  302. fins ;
  303.  
  304. *----------------------------------------------------------------------*
  305. * Sequence de soudage *
  306. *----------------------------------------------------------------------*
  307.  
  308. * Initialisation table de soudage :
  309. tsoud1 = table ;
  310. tsoud1 . vitesse_de_soudage = Vs1 ;
  311. tsoud1 . vitesse_de_deplacement = Vd1 ;
  312. tsoud1 . puissance_de_soudage = Qtot1 ;
  313. tsoud1 . diametre_de_fil = Df1 ;
  314. tsoud1 . vitesse_de_fil = Vf1 ;
  315. tsoud1 . point_de_depart = Pc1 ;
  316.  
  317. * Specification de la sequence :
  318. * Delai de 0,5 s avant de faire la 1ere passe :
  319. soudage tsoud1 point tps0 ;
  320. * 1ere passe :
  321. soudage tsoud1 passe droi (lw1*Z1) ;
  322. * 2e passe :
  323. soudage tsoud1 depla droi Pc2 abso ;
  324. soudage tsoud1 point 600. puis 0. debi 0. ;
  325. soudage tsoud1 passe droi (lw1*Z1) ;
  326. * 3e passe :
  327. soudage tsoud1 depla droi Pc3 abso ;
  328. soudage tsoud1 point 600. puis 0. debi 0. ;
  329. soudage tsoud1 passe droi (lw1*Z1) ;
  330. * 4e passe :
  331. soudage tsoud1 depla droi Pc4 abso ;
  332. soudage tsoud1 point 600. puis 0. debi 0. ;
  333. soudage tsoud1 passe droi (lw1*Z1) ;
  334. * Refroidissement
  335. soudage tsoud1 point 600. puis 0. debi 0. ;
  336.  
  337. * Sequencage maillage des cordons le long de la trajectoire :
  338. tab2 = soudage tsoud1 'MAIL' vc0 'PAS' Xpas1 'TEMP' ;
  339.  
  340. ttps1 = tab2.evolution_maillage.temps ;
  341.  
  342. si ig1 ;
  343. titr ' Trajectoire soudage des 4 passes (vert : depl. outil)' ;
  344. trac (1 2 10) (tsoud1.trajectoire et (aret vc1) et (aret vc2 coul defa) et (aret vc3 coul defa) et (aret vc4 coul defa)) ;
  345. dess tsoud1.evolution_deplacement titr 'Evolution temporelle du deplacement sur la trajectoire' ;
  346. dess tsoud1.evolution_puissance titr 'Evolution tempeorelle de la puissance thermique' ;
  347. *
  348. opti oeil (1. 2. 3.) ;
  349. waam tab2 visu cach ;
  350. *dess (evol manu tab2.temps_calcules (prog 1. pas 1. (dime tab2.temps_calcules))) titr 'Temps calcules' ;
  351. fins ;
  352.  
  353. *----------------------------------------------------------------------*
  354. * Modele Physique *
  355. *----------------------------------------------------------------------*
  356.  
  357. * Accrochage DDL thermique bord raidisseur sur plaque :
  358. clt1 = rela accro (enve Vr0) vp0 (mots T) ;
  359.  
  360. * sinte1 : surface interieure : a retirer de la surface qui convecte
  361. * partie du raidisseur en vis-a-vis de la plaque (surface y = ymin.)
  362. ptry0 = (vr0 coor 2) poin mini ;
  363. sury0 = (vr0 enve) elem appu stri ptry0 ;
  364. sinte1 = sury0 ;
  365. * partie de la plaque en vis-a-vis du raidisseur (maillage non-conforme)
  366. pclt1 = clt1 extr mail nomu ;
  367. supi1 = (enve Vp0) elem appu stri (pclt1 et sury0) ;
  368. sinte1 = sinte1 et supi1 ;
  369. * plan de symetrie en x=0
  370. ptx0 = (V0 coor 1) poin infe 1.e-6 ;
  371. sux0 = (enve V0) elem appu stri ptx0 ;
  372. sinte1 = sinte1 et sux0 ;
  373.  
  374. si ig1 ;
  375. trac (-5 1 3) cach ((enve V0) diff sinte1) boit (boite vc0) titr 'Surface de convection' ;
  376. fins ;
  377.  
  378. *---------------------- Modele / Carateristiques ----------------------*
  379.  
  380. * Conduction / Convection :
  381. mod1 = mode V0 thermique ;
  382. mat1 = mate mod1 rho rho1 k k1 'C' cp1 'TINI' Tini1 ;
  383.  
  384. * Source thermique :
  385. evqt1 = tsoud1.evolution_puissance ;
  386. evxs1 = tsoud1.evolution_deplacement ;
  387. chxs1 = tsoud1.trajectoire coor curv ;
  388. cgxs1 = char traj chxs1 evxs1 ;
  389. mod3 = mode V0 thermique source gaussienne isotrope_transverse ;
  390. mat3 = mate mod3 qtot evqt1 orig cgxs1 rgau Rg1 dire (1. 1. 0.) zgau Zg1 ;
  391.  
  392. * Tables d'evolution des modeles et caracteristiques :
  393. ttps1 = tab2.evolution_maillage.temps ;
  394. tmai1 = tab2.evolution_maillage.maillage ;
  395. nb1 = dime ttps1 ;
  396. tmod1 = table ;
  397. tmod2 = table ;
  398. tmod3 = table ;
  399. tmod4 = table ;
  400. tmat1 = table ;
  401. tmat2 = table ;
  402. tmat3 = table ;
  403. tmat4 = table ;
  404. repe b1 nb1 ;
  405. geoi1 = tmai1 . (&b1 - 1) ;
  406. geoi1 = geoi1 et Vp0 et Vr0 ;
  407. tmod1 . (&b1 - 1) = redu mod1 geoi1 ;
  408. tmat1 . (&b1 - 1) = redu mat1 (tmod1 . (&b1 - 1)) ;
  409. sconv1 = (enve geoi1) diff sinte1 ;
  410. *trac nclk cach sconv1 boit (boite vc0) ;
  411. tmod2 . (&b1 - 1) = mode sconv1 thermique convection ;
  412. tmat2 . (&b1 - 1) = mate (tmod2 . (&b1 - 1)) 'H' h1 'TC' Te1 ;
  413. tmod3 . (&b1 - 1) = redu mod3 geoi1 ;
  414. tmat3 . (&b1 - 1) = redu mat3 (tmod3 . (&b1 - 1)) ;
  415. tmod4 . (&b1 - 1) = mode sconv1 thermique rayonnement infini ;
  416. tmat4 . (&b1 - 1) = mate (tmod4 . (&b1 - 1)) 'EMIS' em1 'T_IN' Te1 ;
  417. fin b1 ;
  418.  
  419. * Chargements MODE / MATE :
  420. cgmod1 = char mode ttps1 tmod1 ;
  421. cgmod2 = char mode ttps1 tmod2 ;
  422. cgmod3 = char mode ttps1 tmod3 ;
  423. cgmod4 = char mode ttps1 tmod4 ;
  424. cgmod0 = cgmod1 et cgmod2 et cgmod3 et cgmod4 ;
  425.  
  426. cgmat1 = char mate ttps1 tmat1 ;
  427. cgmat2 = char mate ttps1 tmat2 ;
  428. cgmat3 = char mate ttps1 tmat3 ;
  429. cgmat4 = char mate ttps1 tmat4 ;
  430. cgmat0 = cgmat1 et cgmat2 et cgmat3 et cgmat4 ;
  431.  
  432. cg0 = cgmod0 et cgmat0 ;
  433.  
  434. *------------------------- Resolution PASAPAS -------------------------*
  435.  
  436. * temps calcules :
  437. ltca1 = tab2.temps_calcules ;
  438. it0 = posi tps0 dans ltca1 ;
  439. ltca1 = (prog 0. pas 0.2 tps0) et (ltca1 extr (lect (it0 + 1) pas 1 (dime ltca1))) ;
  440. si (non icomplet) ;
  441. ltca1 = ltca1 extr (lect 1 pas 1 20) ;
  442. fins ;
  443. list ltca1 ;
  444.  
  445. * temps des sauvegardes :
  446. si icomplet ;
  447. ltsg1 = ltca1 extr (lect 50 pas 50 ((dime ltca1) / 50 * 50)) ;
  448. si ((maxi ltsg1) neg (maxi ltca1)) ;
  449. ltsg1 = ltsg1 et (maxi ltca1) ;
  450. fins ;
  451. list ltsg1 ;
  452. fins ;
  453.  
  454. tpas1 = table ;
  455. tpas1.modele = tire cgmod0 mode 0. ;
  456. tpas1.caracteristiques = tire cgmat0 mate 0. ;
  457. tpas1.chargement = cg0 ;
  458. tpas1.blocages_thermiques = clt1 ;
  459. tpas1.temps_calcules = ltca1 ;
  460. tpas1.temps_sauves = ltca1 ;
  461. tpas1.celsius = vrai ;
  462.  
  463. si icomplet ;
  464. tpas1.temps_sauvegardes = ltsg1 ;
  465. opti sauv 'thermique_SoudageRaidisseur.sauv' ;
  466. sauv ;
  467. fins ;
  468.  
  469.  
  470. pasapas tpas1 ;
  471.  
  472. * Affichage temps de calcul :
  473. duree1 = temp enti horl ;
  474. duree1 = (duree1 / 1000) ;
  475. dmin1 = (duree1 / 60) ;
  476. dsec1 = duree1 - (60 * dmin1) ;
  477. mot1 = chai '***** DUREE DU CALCUL (s) :' dmin1 ' min' dsec1 ' s' ;
  478. mess mot1 ;
  479.  
  480. *------------------------ Petit post-traitement -----------------------*
  481.  
  482. si ig1 ;
  483.  
  484. * Animation du champ de temperature :
  485. if1 = (dime tpas1.temps) - 1 ;
  486. modf1 = tpas1.modeles.if1 ;
  487. mailf1 = (extr modf1 mate conduction) extr mail ;
  488. dep0 = manu chpo mailf1 3 ux 0. uy 0. uz 0. ;
  489. def1 = vide deforme ;
  490. liso1 = prog 0. 50. 100. PAS 150. 1450. ;
  491. vtra1 = vale trac ;
  492. si (ega vtra1 'PSC') ;
  493. opti ftra 'Temperature_soudage6_anim.ps' ;
  494. fins ;
  495. mbox1 = boite Vc0 ;
  496. repe b1 if1 ;
  497. modi1 = tpas1.modeles.(&b1-1) ;
  498. maili1 = (extr modi1 mate conduction) extr mail ;
  499. chti1 = tpas1.temperatures.(&b1-1) ;
  500. defoi1 = defo maili1 dep0 chti1 ;
  501. def1 = def1 et defoi1 ;
  502. si (ega vtra1 'PSC') ;
  503. mot1 = chai format '(F6.1)' 'Temperature (degC) au temps (s) :' (tpas1.temps.(&b1-1)) ;
  504. trac (1 2 3) chti1 maili1 liso1 titr mot1 boit mbox1 ;
  505. fins ;
  506. fin b1 ;
  507. si (neg vtra1 'PSC') ;
  508. trac (1 2 3) anim def1 liso1 boit mbox1 ;
  509. fins ;
  510.  
  511. fins ;
  512.  
  513. si icomplet ; sauv ; fins ;
  514.  
  515. *-------------- F I N S O U D A G E 6 . D G I B I --------------*
  516. fin ;
  517.  
  518.  
  519.  
  520.  
  521.  
  522.  
  523.  
  524.  
  525.  

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