Télécharger waam4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier waam4.dgibi
  2. * section : thermique conduction convection
  3. *----------------------------------------------------------------------*
  4. * W A A M 4 . D G I B I *
  5. *----------------------------------------------------------------------*
  6. * Objet :
  7. * -------
  8. *
  9. * Exemple de simulation thermique d'une fabrication additive en WAAM.
  10. * L'exemple simule la realisation d'un piquage a 45 degres sur un tube
  11. * (tube A). Le depot et le tube A sont en acier 316L.
  12. *
  13. * On modelise la conduction, la convection et le rayonnement avec le
  14. * milieu ambiant. La source de chaleur suit 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 modelisation fait appel a la procedure SOUDAGE pour definir la
  21. * sequence de fabrication et indexer l'appartion des mailles au cours
  22. * du temps le long de la trajectoire de l'outil (sequancage de l'apport
  23. * de matiere).
  24. *
  25. * Cet exemple sert notamment a verifier le bon fonctionnement de
  26. * l'option MAIL de la procedure SOUDAGE utilisee pour le sequencage du
  27. * maillage dans le cas d'une trajectoire circulaire fermee.
  28. *
  29. * Un test de non regression est realise sur la temperature en un point
  30. * et un instant donne du calcul.
  31. *
  32. * Description :
  33. * -------------
  34. *
  35. * Type de calcul : Thermique Transitoire
  36. * Mode de calcul : 3D
  37. * Type d'element : CUB8, PRI6, TET4
  38. * Chargement : Source de chaleur, Convection, Rayonnement,
  39. * Apport de matiere
  40. *
  41. *----------------------------------------------------------------------*
  42. opti dime 3 elem cub8 ;
  43.  
  44. opti trac psc ;
  45.  
  46. * Pour affichages, mettre IG1 a VRAI :
  47. ig1 = faux ;
  48.  
  49. * Pour sortie PS sequencage apport de matiere, mettre ISEQMAIL a VRAI :
  50. iseqmail = faux ;
  51.  
  52. * Pour sortie PS champ temperature au cours du temps, mettre ITANIM a VRAI :
  53. iT_anim = faux ;
  54.  
  55. * Pour calcul complet, mettre ICOMPLET a VRAI :
  56. icomplet = faux ;
  57.  
  58. *----------------------------------------------------------------------*
  59. * Parametres du modele *
  60. *----------------------------------------------------------------------*
  61.  
  62. *------------------ Parametres geometrie et maillage ------------------*
  63.  
  64. * Rae1 : rayon exte. tube A (m)
  65. * Rai1 : rayon inte. tube A (m)
  66. * Rbe1 : rayon exte. tube B (m)
  67. * Rbi1 : rayon inte. tube B (m)
  68. * Rbm1 : rayon moyen tube B (m)
  69. * Angl1 : angle axes tubes A et B (degres)
  70. * Np1 : nombre de passe de fabrication
  71. * epp1 : epaisseur d'une passe
  72. * de1 : densite maillage passe (tube B)
  73. * de2 : densite maillage tube A
  74. * ne1 : nombre d'element sur la trajectoire circulaire d'une passe
  75. * ne2 : nombre d'element sur le perimetre du tube fabrique
  76. Rae1 = 150.e-3 ;
  77. Rai1 = 135.e-3 ;
  78. Rbe1 = 100.e-3 ;
  79. Epb1 = 10.e-3 ;
  80. Rbi1 = Rbe1 - Epb1 ;
  81. Rbm1 = 0.5 * (Rbe1 + Rbi1) ;
  82. Angl1 = 45. ;
  83. Np1 = 5 ;
  84. epp1 = 4.e-3 ;
  85. de1 = 2.e-3 ;
  86. de2 = 10. * de1 ;
  87. ne1 = 40 ;
  88. ne2 = enti (2. * pi * Rbm1 / de1 + 0.5) ;
  89. ne2 = ne2 / 2 ;
  90. ne2 = ne2 - (ne2 - (ne2 / 4 * 4)) ; comm ne2 doit etre un multiple de 4 ;
  91.  
  92. *----------------------- Parametres fabrication -----------------------*
  93.  
  94. * Us1/Us2 : tension 1ere/autres passes
  95. * Is1/Is2 : intensite 1ere/autres passes
  96. * eta1 : rendement
  97. * dtini1 : duree prechauffage 1ere passe (s)
  98. * dtpasse1: temps de pause entre 2 passes (s)
  99. * Dfil1 : diametre de fil (m)
  100. * Vfil1 : vitesse de debit de fil (m/s)
  101. * Vs1 : vitesse de soudage (m/s)
  102. * Rg1 : rayon source Gaussienne (m)
  103. * Zg1 : profondeur source Gaussienne (m)
  104. Us1 = 30. ;
  105. Us2 = 25. ;
  106. Is1 = 250. ;
  107. Is2 = 250. ;
  108. eta1 = 0.8 ;
  109. dtini1 = 1.0 ;
  110. dtpasse1 = 600. ;
  111. Dfil1 = 2.e-3 ;
  112. Vfil1 = 3.5 / 60. ;
  113. Vs1 = 300.e-3 / 60. ;
  114. Rg1 = 5.e-3 ;
  115. Zg1 = 3.e-3 ;
  116.  
  117. *----------------- Parametres modelisation thermique ------------------*
  118.  
  119. * T0 : temperature initiale et ambiante (degC)
  120. * H1 : coefficient de convection (W/m2)
  121. * Em1 : emissivite
  122. ZeroC1 = 273.16 ;
  123. T0 = 20. + ZeroC1 ;
  124. H1 = 20. ;
  125. Em1 = 0.5 ;
  126.  
  127. *----------------------------------------------------------------------*
  128. * MAILLAGE *
  129. *----------------------------------------------------------------------*
  130.  
  131. * Repere :
  132. O1 = 0 0 0 ;
  133. X1 = 1 0 0 ;
  134. y1 = 0 1 0 ;
  135. z1 = 0 0 1 ;
  136.  
  137. *------------------------- Maillage du depot --------------------------*
  138.  
  139. * Points et lignes de de construction :
  140. PCA1 = 1 0 0 ;
  141. PCB1 = 1 0 0 ;
  142. PA1 = PCA1 plus (Rae1 * Y1) ;
  143. PBE1 = PCB1 plus (Rbe1 * Y1) ;
  144. PBI1 = PCB1 plus (Rbi1 * Y1) ;
  145. RBM1 = (0.2 * Rbi1) + (0.8 * Rbe1) ;
  146. PBM1 = PCB1 plus (Rbm1 * Y1) ;
  147.  
  148. * Diametres ext., int. et moyen du tube B oriente suivant Ox ;
  149. LCBM1 = cerc ne1 rota 360. PBM1 O1 PCB1 ;
  150. LCBE1 = cerc ne2 rota 360. PBE1 O1 PCB1 elim ;
  151. LCBI1 = cerc ne2 rota 360. PBI1 O1 PCB1 elim ;
  152. LCBM1 = LCBM1 coul vert ;
  153. LCBE1 = LCBE1 coul roug ;
  154. LCBI1 = LCBI1 coul bleu ;
  155.  
  156. * Rotation lignes de contruction tube B a Angl1 de Ox :
  157. LCBE1 LCBI1 LCBM1 PCB1 = LCBE1 LCBI1 LCBM1 PCB1 tour Angl1 Y1 O1 ;
  158.  
  159. * Projections surfaces sur tube A :
  160. Sb1 = LCBE1 regl LCBI1 dini de1 dfin de1 ;
  161. Sb2 = Sb1 proj cyli (O1 moin PCB1) cyli O1 PCA1 PA1 ;
  162. LCBM2 = LCBM1 proj cyli (O1 moin PCB1) cyli O1 PCA1 PA1 ;
  163. PCB2 = PCB1 proj cyli (O1 moin PCB1) cyli O1 PCA1 PA1 ;
  164. *trac (LCBE1 et LCBI1 et LCBM1) (1000. * PCB1) titr mot1 ;;
  165.  
  166. * Construction maillage passe par passe &
  167. * Construction de la trajectoire :
  168. nep1 = enti (epp1/de1+0.5) ;
  169. VTB0 = vide maillage ;
  170. ltrj1 = LCBM2 coul roug ;
  171. ttrj1 = table ;
  172. Si1 = Sb2 ;
  173. Pi1 = Si1 poin proc O1 ;
  174. repe bpasse1 Np1 ;
  175. i1 = &bpasse1 ;
  176.  
  177. * Maillage passe :
  178. Vi1 = Si1 volu nep1 tran (epp1*PCB1) ;
  179. Vi1 = Vi1 coul i1 ;
  180. VTB0 = VTB0 et Vi1 ;
  181. Si1 = Vi1 face 2 ;
  182. *trac (Si1 et (aret VTB0)) ;
  183.  
  184. * Maillage trajectoire :
  185. Pi2 = Si1 poin proc O1 ;
  186. ltrj1 = ltrj1 plus (Pi2 moin Pi1) ;
  187. ttrj1 . i1 = ltrj1 ;
  188. Pi1 = Pi2 ;
  189. ltrj1 = inve ltrj1 ;
  190. *trac (ltrj1 et (aret VTB0)) ;
  191. fin bpasse1 ;
  192.  
  193. *------------------------- Maillage du tube A -------------------------*
  194.  
  195. * Maillage tranche tube A sous surface base tube B :
  196. Stai1 = Sb2 proj cyli (O1 moin PCB1) cyli O1 X1 (Rai1 * Y1) ;
  197. Vta1 = Sb2 volu 2 Stai1 ;
  198.  
  199. * Surface exterieure :
  200. dXA1 = Rbe1 / (cos 45.) + (Rae1 + Rbe1 * 0.2) ;
  201. PTA1 = PCB2 moin (dXA1 * X1) moin ((PCB2 coor 3) * Z1) plus (Rae1 * Y1) ;
  202. LTA1 = PTA1 cerc dini de2 dfin de2 rota O1 X1 180. ;
  203. PTA2 = LTA1 poin (nbno LTA1) ;
  204. LTA3 = LTA1 plus (2. * dXA1 * X1) ;
  205. LTA3 = inve LTA3 ;
  206. PTA3 = LTA3 poin 1 ;
  207. PTA4 = LTA3 poin (nbno LTA3) ;
  208. LTA2 = PTA2 droi PTA3 dini de2 dfin de2 ;
  209. LTA4 = PTA4 droi PTA1 dini de2 dfin de2 ;
  210.  
  211. * Contour du tube B :
  212. tsb2 = (cont sb2) part nesc conn ;
  213. csbe2 = tsb2 . 1 ;
  214. si ((mesu csbe2) < (mesu tsb2 . 2)) ;
  215. csbe2 = tsb2 . 2 ;
  216. fins ;
  217.  
  218. * Contour global et maillage surface exte. tube A :
  219. LTA0 = LTA1 et LTA2 et LTA3 et LTA4 et (inve csbe2) ;
  220. STA1 = surf LTA0 cyli O1 X1 ;
  221.  
  222. * Surace interieure :
  223. LTA5 = LTA1 proj coni ((PTA1 plus PTA2) * 0.5) cyli O1 X1 (Rai1 * Y1) ;
  224. PTA5 = LTA5 poin 1 ;
  225. PTA6 = LTA5 poin (nbno LTA5) ;
  226. LTA7 = LTA5 plus (2. * dXA1 * X1) ;
  227. PTA7 = LTA7 poin (nbno LTA7) ;
  228. PTA8 = LTA7 poin 1 ;
  229. LTA6 = PTA6 droi PTA7 dini de2 dfin de2 ;
  230. LTA8 = PTA8 droi PTA5 dini de2 dfin de2 ;
  231.  
  232. * Contour exte. Stai1 :
  233. tstai1 = (cont Stai1) part nesc conn ;
  234. cstai1 = tstai1 . 1 ;
  235. si ((mesu cstai1) < (mesu tstai1 . 2)) ;
  236. cstai1 = tstai1 . 2 ;
  237. fins ;
  238.  
  239. * Contour global et maillage surface inte. tube A :
  240. LTA02 = (LTA5 et LTA6 et LTA7 et LTA8) elem comp PTA5 PTA5 ;
  241. LTA02 = LTA02 et (inve cstai1) ;
  242. STA2 = surf LTA02 cyli O1 X1 ;
  243.  
  244. * Bords :
  245. STA3 = LTA0 regl LTA02 2 ;
  246.  
  247. * Volume partie sup. tube A :
  248. VTA2 = volu ((sta1 et sta2 et sta3) chan tri3) ;
  249.  
  250. * Extensions axiales :
  251. ptx0 = (VTA2 coor 1) poin mini ;
  252. ptx1 = (VTA2 coor 1) poin maxi ;
  253. STAx0 = (enve VTA2) elem appu stri ptx0 ;
  254. STAx1 = (enve VTA2) elem appu stri ptx1 ;
  255. VTAx0 = STAx0 volu tran 1 (-2. * Rae1 * X1) ;
  256. VTAx1 = STAx1 volu tran 1 (+2. * Rae1 * X1) ;
  257. VTA3 = VTAx0 et VTAx1 ;
  258.  
  259. * Volume partie infe. (opposee au piquage) :
  260. ptz0 = ((VTA2 et VTA3) coor 3) poin infe 1.e-6 ;
  261. ptz0 = (ptz0 coor 2) poin supe 0. ;
  262. STAz0 = ((VTA2 et VTA3) enve) elem appu stri ptz0 ;
  263. VTA4 = STAz0 volu 10 rota -180. O1 X1 ;
  264.  
  265. * Maillage total tube A :
  266. VTA0 = (VTA1 et VTA2 et VTA3 et VTA4) coul turq ;
  267. elim VTA0 1.e-6 ;
  268.  
  269. *------------------------ Affichages Maillages ------------------------*
  270.  
  271. oeil1 = -100 100 80 ;
  272. si ig1 ;
  273. trac oeil1 face (VTA0 et VTB0) titr 'Maillage tube (turquoise) et piquage par depot WAAM' ;
  274. trac (100*Y1) face (VTA0 et VTB0) boit (boite VTB0) titr 'Maillage tube et piquage : detail piquage vue cote' ;
  275. fins ;
  276.  
  277. *----------------------------------------------------------------------*
  278. * Sequence de fabrication *
  279. *----------------------------------------------------------------------*
  280.  
  281. * Puissances de soudage :
  282. Qs1 = eta1 * Us1 * Is1 ;
  283. Qs2 = eta1 * Us2 * Is2 ;
  284.  
  285. * Vitesse de deplacement entre 2 passes pour un temps de pause donnee :
  286. ldepl1 = ((ttrj1 . 1) poin 1) droi 1 ((ttrj1 . 2) poin 2) ;
  287. vdepl1 = (mesu ldepl1) / dtpasse1 ;
  288.  
  289. * Donnees d'entree :
  290. tfab1 = table ;
  291. tfab1 . vitesse_de_soudage = Vs1 ;
  292. tfab1 . vitesse_de_deplacement = vdepl1 ;
  293. tfab1 . puissance_de_soudage = Qs2 ;
  294. tfab1 . debit_de_fil = 0.25 * pi * Dfil1 * Dfil1 * Vfil1 ;
  295. tfab1 . point_de_depart = (ttrj1 . 1) poin 1 ;
  296. tfab1 . temps_de_coupure = 0.5 ;
  297.  
  298. * Sequence de fabrication :
  299. soudage tfab1 point dtini1 puis Qs1 debi 0. ;
  300. nfab1 = dime ttrj1 ;
  301. repe bfab1 nfab1 ;
  302. si (&bfab1 ega 1) ;
  303. soudage tfab1 passe mail (ttrj1 . &bfab1) puis Qs1 ;
  304. sino ;
  305. soudage tfab1 passe mail (ttrj1 . &bfab1) ;
  306. fins ;
  307. si (&bfab1 ega nfab1) ; quit bfab1 ; fins ;
  308. soudage tfab1 depla droi ((ttrj1 . (&bfab1+1)) poin 1) abso vite vdepl1 ;
  309. fin bfab1 ;
  310. soudage tfab1 point dtpasse1 puis 0. debi 0. ;
  311.  
  312. * Affichages sequence de fabrication :
  313. si ig1 ;
  314. opti oeil (PCB2 plus Z1) ;
  315. trac tfab1.trajectoire titr 'trajectoire realisation piquage' ;
  316. dess tfab1.evolution_deplacement titr 'evolution deplacement' ;
  317. dess tfab1.evolution_puissance titr 'evolution puissance thermique' ;
  318. dess tfab1.evolution_debit titr 'evolution debit apport de matiere' ;
  319. fins ;
  320.  
  321. *-------------------- Sequencage apport de matiere --------------------*
  322.  
  323. * Sequencage maillage le long de la trajectoire :
  324. tseq1 = soudage tfab1 mail VTB0 pas 10.e-3 temp ;
  325. ttps1 = tseq1 . evolution_maillage . temps ;
  326. tmai1 = tseq1 . evolution_maillage . maillage ;
  327.  
  328. oeil2 = 100 50 150. ;
  329. si iseqmail ;
  330. si (ega (vale trac) 'PSC') ; opti ftra 'waam4_sequ_mail.ps' ; fins ;
  331. mbox1 = boite VTB0 ;
  332. i1 = 0 ;
  333. repe b1 (dime tmai1) ;
  334. geoi1 = tmai1 . i1 ;
  335. trac oeil2 face (geoi1 et VTA0) titr 'Sequencage du maillage du depot' boit mbox1 ;
  336. i1 = i1 + 1 ;
  337. fin b1 ;
  338. fins ;
  339.  
  340. *----------------------------------------------------------------------*
  341. * Modelisation Thermique *
  342. *----------------------------------------------------------------------*
  343.  
  344. * Donnees materiau :
  345. t316L = biblio '316L' ;
  346. rho1 = ipol t316L.'RHO' (T0-Zeroc1) ;
  347. evk1 = t316L.'K' + Zeroc1 absc ;
  348. evc1 = t316L.'C' + Zeroc1 absc ;
  349.  
  350. * Modele de conduction :
  351. VT0 = VTA0 et VTB0 ;
  352. mod1 = mode VT0 thermique conduction ;
  353. mat1 = mate mod1 'RHO' rho1 'C' evc1 'K' evk1 'TINI' T0 ;
  354.  
  355. * Modele source thermique :
  356. * Evolution de la position de la source :
  357. evqt1 = tfab1.evolution_puissance ;
  358. evxs1 = tfab1.evolution_deplacement ;
  359. chxs1 = tfab1.trajectoire coor curv ;
  360. cgxs1 = char traj chxs1 evxs1 ;
  361.  
  362. * Evolution de la direction de la source :
  363. * On donne la normale au tube en chaque point de la trajectoire.
  364. * L'axe du tube est selon (O,x), la normale est donc selon la coord.
  365. * (y,z) de chaque point.
  366. ltraj1 = tfab1.trajectoire ;
  367. ldir1 = enum ;
  368. ltpx1 = prog ;
  369. evtxs1 = evol manu (evxs1 extr ordo) (evxs1 extr absc) ;
  370. nbdir1 = nbno ltraj1 ;
  371. tps0 = -1. ;
  372. ddir1 = vide maillage ;
  373. repe bdir1 nbdir1 ;
  374. pi1 = ltraj1 poin &bdir1 ;
  375. xsi1 = extr chxs1 scal pi1 ;
  376. txi1 = ipol evtxs1 xsi1 ;
  377. si (txi1 > tps0) ;
  378. pni1 = 0. (pi1 coor 2) (pi1 coor 3) ;
  379. ldir1 = ldir1 et (enum pni1) ;
  380. ltpx1 = ltpx1 et txi1 ;
  381. ddir1 = ddir1 et (pi1 droi 1 (pi1 plus pni1)) ;
  382. fins ;
  383. fin bdir1 ;
  384. cgdir1 = char dire ltpx1 ldir1 ;
  385.  
  386. * Visualisation :
  387. si ig1 ;
  388. mbox1 = boite VTB0 ;
  389. ddir1 = ddir1 coul oran ;
  390. trac (100 100 20) face (VTA0 et VTB0 et ddir1) boit mbox1 titr 'Directions d''orientation de la source de chaleur' ;
  391. fins ;
  392.  
  393. * Definition du modele :
  394. mod2 = mode VT0 thermique source gaussienne isotrope_transverse ;
  395. mat2 = mate mod2 'QTOT' evqt1 'ORIG' cgxs1 'DIRE' cgdir1 'RGAU' Rg1 'ZGAU' Zg1 ;
  396.  
  397. * Sequencage modeles :
  398. ltps1 = prog tabl ttps1 ;
  399. nb1 = dime tmai1 ;
  400. lgeo0 = enum nb1*vta0 ;
  401. lgeo1 = enum 'TABL' tmai1 ;
  402.  
  403. * Vectorisation des operations sur les listes de maillages et de modeles
  404. * en activant la parallelisation :
  405. opti para vrai ;
  406. lgeo1 = lgeo0 et lgeo1 ;
  407. lgeo2 = enve lgeo1 ;
  408. lmod1 = redu lgeo1 (mod1 et mod2) ;
  409. lmat1 = redu lmod1 (mat1 et mat2) ;
  410. lmod2 = mode lgeo2 thermique convection ;
  411. lmat2 = mate lmod2 'H' h1 'TC' T0 ;
  412. lmod3 = mode lgeo2 thermique rayonnement infini ;
  413. lmat3 = mate lmod3 'EMIS' Em1 'T_IN' T0 ;
  414.  
  415. * Chargements MODE / MATE :
  416. opti para faux ;
  417. * Retour en sequentiel : chaque commande sert de point d'arret
  418. cgmod1 = char mode ltps1 lmod1 ;
  419. cgmod2 = char mode ltps1 lmod2 ;
  420. cgmod3 = char mode ltps1 lmod3 ;
  421. cgmod0 = cgmod1 et cgmod2 et cgmod3 ;
  422.  
  423. cgmat1 = char mate ltps1 lmat1 ;
  424. cgmat2 = char mate ltps1 lmat2 ;
  425. cgmat3 = char mate ltps1 lmat3 ;
  426. cgmat0 = cgmat1 et cgmat2 et cgmat3 ;
  427.  
  428. *------------------------- Resolution PASAPAS -------------------------*
  429.  
  430. ltca0 = tseq1.temps_calcules ;
  431. si icomplet ;
  432. ltca1 = ltca0 ;
  433. sino ;
  434. ltca1 = ltca0 extr (lect 1 pas 1 11) ;
  435. fins ;
  436. list ltca1 ;
  437.  
  438. tpas1 = table ;
  439. tpas1.modele = cgmod0 ;
  440. tpas1.caracteristiques = cgmat0 ;
  441. tpas1.temps_calcules = ltca1 ;
  442.  
  443.  
  444. pasapas tpas1 ;
  445.  
  446. * Affichage temps de calcul :
  447. duree1 = temp horl ;
  448. duree1 = (duree1 / 1000) ;
  449. dmin1 = (duree1 / 60) ;
  450. dsec1 = duree1 - (60 * dmin1) ;
  451. mot1 = chai '***** DUREE DU CALCUL (s) :' dmin1 ' min' dsec1 ' s' ;
  452. mess mot1 ;
  453.  
  454. si icomplet ;
  455. opti sauv 'waam4.sauv' ;
  456. sauv ;
  457. fins ;
  458.  
  459. *------------------------ Petit post-traitement -----------------------*
  460.  
  461. si iT_anim ;
  462. si (ega (vale trac) 'PSC') ; opti ftra 'waam4_T_anim.ps' ; fins ;
  463.  
  464. * Sorties PS du champ de temperature :
  465. * Image toutes les :
  466. * - dtimag1 quand puiss. thermique non nulle
  467. * - FFwd1*dtimag1 sinon.
  468. dtimag1 = 2. ;
  469. si (ega (vale trac) 'PSC') ; dtimag1 = 0.5 ; fins ;
  470. FFwd1 = 10. ;
  471. liso1 = prog 50. 100. PAS 150. 1450. ;
  472. mbox1 = boite VTB0 ;
  473. tpsf1 = tpas1.temps . ((dime tpas1.temps) - 1) ;
  474. ifin1 = faux ;
  475. tpsi1 = 0. ;
  476. repe b1 ;
  477. modi1 = peche tpas1 modeles tpsi1 ipol ;
  478. geoi1 = (extr modi1 mate conduction) extr mail ;
  479. chti1 = peche tpas1 temperatures tpsi1 ipol ;
  480. chti1 = chti1 - ZeroC1 ;
  481. ptrj1 = tire cgxs1 traj tpsi1 ;
  482. Ri1 = ((ptrj1 coor 2) ** 2) + ((ptrj1 coor 3) ** 2) ** 0.5 ;
  483. yi1 = ptrj1 coor 2 ;
  484. thetai1 = yi1 / Ri1 * 180. / pi ;
  485. depl geoi1 tour thetai1 O1 X1 ;
  486. mot1 = chai format '(F6.1)' 'Temperature (degC) au temps (s) :' tpsi1 ;
  487. trac (100 0 100) chti1 geoi1 liso1 titr mot1 boit mbox1 ;
  488. depl geoi1 tour thetai1 X1 O1 ;
  489. tpsi1 = tpsi1 + dtimag1 ;
  490. qti1 = ipol evqt1 tpsi1 ;
  491. si (qti1 < 1.) ; comm si puissance nulle, on accelere le film (xFFwd1) ;
  492. tpsi1 = (FFwd1 - 1.) * dtimag1 + tpsi1 ;
  493. fins ;
  494. si (tpsi1 > tpsf1) ;
  495. si ifin1 ;
  496. quit b1 ;
  497. sino ;
  498. tpsi1 = tpsf1 ;
  499. ifin1 = vrai ;
  500. fins ;
  501. fins ;
  502. fin b1 ;
  503.  
  504. fins ;
  505.  
  506. *------------------------- Test non-regression ------------------------*
  507.  
  508. ttest1 = tpas1.temps.10 ;
  509. temp1 = peche tpas1 temperatures ttest1 ipol ;
  510. tmax1 = maxi abs temp1 ;
  511. tmaxref1 = 4272.9 ;
  512. si ((abs (tmax1 - tmaxref1)) > 1.) ;
  513. erre 5 ;
  514. fins ;
  515.  
  516. *------------------ F I N W A A M 4 . D G I B I -----------------*
  517. fin ;
  518.  
  519.  
  520.  
  521.  
  522.  

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