Télécharger transnon.procedur

Retour à la liste

Numérotation des lignes :

  1. * TRANSNON PROCEDUR FANDEUR 22/01/19 21:15:17 11256
  2. 'DEBPROC' TRANSNON PRECED*'TABLE';
  3. *--------------------------------------------------------------------------------*
  4. * *
  5. * T R A N S N O N *
  6. * --------------- *
  7. * *
  8. * RESOLUTION D'UN PROBLEME DE THERMIQUE TRANSITOIRE NON-LINEAIRE *
  9. * METHODE A UN PAS DE TEMPS ( THETA SCHEMA ) *
  10. * *
  11. * ETAB, TABLE CONTENANT EN ENTREE : *
  12. * *
  13. * INDICE 'TET1' CHAMP DE TEMPERATURE INITIAL AU PAS 0 *
  14. * INDICE 'RAYO' LOGIQUE VALANT VRAI POUR UNE CONDITION *
  15. * DE RAYONNEMENT *
  16. * INDICE 'EMISSIVITE' MCHAML DECRIVANT LES FACTEURS D'EMISSIVITE *
  17. * NOM DE LA COMPOSANTE : EMIS *
  18. * INDICE 'CELSIUS' LOGIQUE VALANT VRAI SI L'UNITE EST LE *
  19. * DEGRE CELSIUS (CAPITAL SI RAYONNEMENT) *
  20. * INDICE 'MOD_THE' OBJET MODELE THERMIQUE *
  21. * INDICE 'MOD_CON' OBJET MODELE CONVECTION *
  22. * INDICE 'BLOCAGES_THERMIQUES' MATRICE DE BLOCAGE *
  23. * INDICE 'MAT_THE' OBJET MATERIAU THERMIQUE. *
  24. * CE CHAMP PEUT AVOIR DES COMPOSANTES DE *
  25. * TYPE FLOTTANT OU EVOLUTION (ABS-ORD). *
  26. * ABS : VARIABLE QUELCONQUE *
  27. * ORD : VALEUR DE LA COMPOSANTE CONCERNE POUR *
  28. * LA VALEUR ABS. *
  29. * INDICE 'MAT_CON' OBJET MATERIAU CONVECTION *
  30. * INDICE 'CHARGEMENT' CHARGEMENT DECRIVANT LES *
  31. * VALEURS DES VARIABLES EXTERNES (EX: TE, *
  32. * FLUX,TEMPERATURES IMPOSEES ,...) *
  33. * INDICE 'PHASE' LOGIQUE VALANT VRAI SI CHGT DE PHASE 'PARFAIT' *
  34. * INDICE 'SOLUBILITE' LOGIQUE VALANT VRAI SI CHGT DE PHASE 'SOLUBILITE' *
  35. * INDICE 'TEMPS0' TEMPS INITIAL (CORRESPOND A TET1) *
  36. * INDICE 'TEM_CALC' LISTREEL : TEMPS DES RESULTATS A CALCULER *
  37. * INDICE 'RELAXATION_THETA' VALEUR DU COEFFICIENT DE RELAXATION *
  38. * (VALEUR PAR DEFAUT 0.5) *
  39. * INDICE 'SOUS_RELAXATION' VALEUR DU COEFF. DE SOUS-RELAXATION *
  40. * (VALEUR PAR DEFAUT 0.5) *
  41. * INDICE 'PRECISION' VALEUR DU CRITERE DE CONVERGENCE *
  42. * (VALEUR PAR DEFAUT 1.E-4) *
  43. * INDICE 'MAXITERATION' NOMBRE MAXIMUM D'ITERATIONS AUTORISEES *
  44. * INDICE 'PROJECTION' LOGIQUE VALANT VRAI SI COUPLAGE ET SI LE *
  45. * MAILLAGE DE LA MECANIQUE ET DE LA THERMIQUE *
  46. * EST DIFFERENT *
  47. * INDICE 'TTOL' PRECISION ABSOLUE pour la THERMIQUE et la DIFFUSION *
  48. * *
  49. * ETAB CONTIENT EN SORTIE *
  50. * *
  51. * INDICE INITIAL(2) DERNIER CHAMP DE TEMPERATURE CALCULE *
  52. * *
  53. * REMARQUE : CERTAINES VARIABLES ONT DES NOMS IMPOSES. *
  54. * LE TEMPS : TEMPS *
  55. * LA TEMPERATURE : T *
  56. * LA TEMPERATURE EXTERIEURE : TE *
  57. * LES TEMPERATURES IMPOSEES : TIMP *
  58. * LES FLUX : Q *
  59. * *
  60. * LE NOM DES AUTRES VARIABLES EST LAISSE AU CHOIX DE *
  61. * L'UTILISATEUR *
  62. * *
  63. **********************************************************************************
  64. * Recuperation de l'information contenue dans ETAB *
  65. **********************************************************************************
  66. ERR_VAL = FAUX ;
  67. STAB ='TABL' ;
  68. ETAB = PRECED.'WTABLE';
  69.  
  70. NBPART = ETAB.'NBPART';
  71. LPRIM1 = ETAB.'PRIM_TOT' ;
  72. LDUAL1 = ETAB.'DUAL_TOT' ;
  73. LPRIM2 =('MOTS' 'TEMP') 'ET' LPRIM1 ;
  74. NBDUA1 = 'DIME' LDUAL1 ;
  75.  
  76. * Recuperation des tolerence absolues (Tout est dans l'indice TTOL)
  77. 'SI' ('NEG' ETAB.'TTOL' ('MOT' 'INCONNU')) ;
  78. LOG_TOL = VRAI ;
  79. TTOL = ETAB.'TTOL' ;
  80. 'SINO' ;
  81. LOG_TOL = FAUX ;
  82. 'FINS' ;
  83. XPETI = ('VALE' 'PETIT') / ('VALE' 'PREC') ;
  84.  
  85. * Appel a PAS_MODL :
  86. * Mise a jour des indices de PRECED.WTABLE relatifs aux modeles
  87. * si PRECED.WTABLE.MODELE a ete modifie
  88. PAS_MODL PRECED ;
  89.  
  90. MODMEC = ETAB.'MO_TOT' ; 'COMM' 'Formulation MECANIQUE & POREUX' ;
  91.  
  92. MODTOT = ETAB.'MOD_TOT' ; 'COMM' 'Formulation THERMIQUE & THERMOHYDRIQUE & DIFFUSION & CHANGEMENT_PHASE' ;
  93. MATTOT0 ='REDU' MODTOT ETAB.'MAT_TOT' ;
  94. MAILTOT ='EXTR' MODTOT 'MAIL';
  95.  
  96. * -------------- presence de conduction
  97. ICONDU = ETAB.'CONDUCTION';
  98. 'SI' ICONDU;
  99. MATCOND0=ETAB.'MAT_COND';
  100. MODCOND =ETAB.'MOD_COND';
  101. 'FINS';
  102.  
  103. * -------------- presence de convection
  104. ICONVE= ETAB.'CONVECTION' ;
  105. 'SI' ICONVE;
  106. MATCONV0=ETAB.'MAT_CON';
  107. MODCONV =ETAB.'MOD_CON';
  108. 'FINS';
  109.  
  110. * -------------- presence d'advection
  111. IADVEC=ETAB.'ADVECTION';
  112. 'SI' IADVEC;
  113. MATADV0 = ETAB.'MAT_ADV';
  114. MODADV = ETAB.'MOD_ADV';
  115. 'FINS';
  116.  
  117. * -------------- presence de rayonnement
  118. IRADI = ETAB.'RAYO';
  119. 'SI' IRADI ;
  120. MATRAY0 = ETAB.'MAT_RAY';
  121. MODRAY = ETAB.'MOD_RAY';
  122. 'SI' ('NON' ('EXIS' ETAB 'CTE_STEFAN_BOLTZMANN')) ;
  123. CTE_SB = 5.67037321E-8 ;
  124. 'MESS' 'ATTENTION ! Constante de Stefan-Boltzmann en unites SI' ;
  125. 'MESS' ' soit ' CTE_SB ;
  126. ETAB.'CTE_STEFAN_BOLTZMANN' = CTE_SB ;
  127. 'FINS' ;
  128. 'FINS' ;
  129.  
  130. * -------------- presence de sources de chaleur
  131. ISOURQ = ETAB.'SOURCE_Q';
  132. 'SI' ISOURQ ;
  133. MATSOQ0 = ETAB.'MAT_SOQ';
  134. MODSOQ = ETAB.'MOD_SOQ';
  135. 'FINS';
  136.  
  137. * -------------- presence de changement de phase 'PARFAIT'
  138. IPHASE=ETAB.'PHASE';
  139. 'SI' IPHASE;
  140. MODPHA = ETAB.'MOD_PHA' ;
  141. MATPHA0 = ETAB.'MAT_PHA' ;
  142. BLOPHA = ETAB.'BLO_PHASE';
  143. PROPHA0 = ETAB.'PROPPHAS' ;
  144. PROPHA1 = PROPHA0 ;
  145. 'FINS';
  146.  
  147. * -------------- presence de changement de phase 'SOLUBILITE'
  148. ISOLU=ETAB.'SOLUBILITE';
  149. 'SI' ISOLU;
  150. BLOSOL = ETAB.'BLO_SOL' ;
  151. 'FINS';
  152.  
  153. * -------------- presence Modele Metallurgie (T.L. en juin 2018)
  154. IMETAL = ETAB.'FOR_METALLU' ;
  155. 'SI' IMETAL;
  156. MODMET = ETAB.'MOD_MET' ;
  157. PHA_MET0= PRECED.'ESTIMATION'.'PROPORTIONS_PHASES' ;
  158. 'FINS';
  159.  
  160. * -------------- presence de thermohydrohydrique
  161. ITHM1=ETAB.'THM1';
  162. 'SI' ITHM1;
  163. MATTHM0 = ETAB.'MAT_THM';
  164. MODTHM = ETAB.'MOD_THM';
  165. 'FINS';
  166.  
  167. * -------------- presence de diffusion
  168. IDIFF=ETAB.'FOR_DIFF';
  169. 'SI' IDIFF;
  170. MATDIF0 = ETAB.'MAT_DIF' ;
  171. MODDIF = ETAB.'MOD_DIF' ;
  172. 'FINS';
  173.  
  174. * CB215821 : Je crois que 'RIBLO_T' et 'LISEA_T' ne sont jamais remplis... (Voir pasapas ; pas_defa ; transnon)
  175. IRIB = 'NEG' ETAB.'RIBLO_T' 'INCONNU' ;
  176. 'SI' IRIB;
  177. RIBL_VAL = ETAB.'RIBLO_T' ;
  178. LISE_VAL = ETAB.'LISEA_T' ;
  179. 'FINS' ;
  180.  
  181. * Parallelisation automatique des calculs ?
  182. IPARAL = 'EGA' ETAB.'PROCESSEURS' 'AUTOMATIQUE' ;
  183. PARTLOCA = 'EGA' ETAB.'PROCESSEURS' 'COMPORTEMENT' ;
  184. 'OPTI' 'PARA' IPARAL ;
  185.  
  186. 'SI' IPARAL ;
  187.  
  188. MODRTOT ='PART' 'ARLE' MODTOT NBPART ;
  189. MATRTOT0 ='REDU' MATTOT0 MODRTOT ;
  190. 'SI' ICONDU;
  191. MODRCOND='PART' 'ARLE' MODCOND NBPART ;
  192. MATRCOND='REDU' MATCOND0 MODRCOND ;
  193. 'FINS';
  194. 'SI' ICONVE;
  195. MODRCON ='PART' 'ARLE' MODCONV NBPART ;
  196. MATRCON ='REDU' MATCONV0 MODRCON ;
  197. 'FINS';
  198. 'SI' IADVEC;
  199. MODRADV ='PART' 'ARLE' MODADV NBPART ;
  200. MATRADV ='REDU' MATADV0 MODRADV ;
  201. 'FINS' ;
  202. 'SI' IRADI;
  203. MODRRAY ='PART' 'ARLE' MODRAY NBPART ;
  204. MATRRAY ='REDU' MATRAY0 MODRRAY ;
  205. 'FINS';
  206. 'SI' ISOURQ;
  207. MODRSOQ ='PART' 'ARLE' MODSOQ NBPART ;
  208. MATRSOQ ='REDU' MATSOQ0 MODRSOQ ;
  209. 'FINS';
  210. 'SI' IPHASE;
  211. MODRPHA ='PART' 'ARLE' MODPHA NBPART ;
  212. MATRPHA ='REDU' MATPHA0 MODRPHA ;
  213. 'FINS';
  214. 'SI' IMETAL;
  215. MODRMET ='PART' 'ARLE' MODMET NBPART ;
  216. PHA_MET0='REDU' PHA_MET0 MODRMET ;
  217. 'FINS';
  218. 'SI' ITHM1;
  219. MODRTHM ='PART' 'ARLE' MODTHM NBPART ;
  220. MATRTHM ='REDU' MATTHM0 MODRTHM ;
  221. 'FINS' ;
  222. 'SI' IDIFF ;
  223. MODRDIFF='PART' 'ARLE' MODDIF NBPART ;
  224. MATRDIFF='REDU' MATDIF0 MODRDIFF ;
  225. 'FINS' ;
  226.  
  227. 'SINO' ;
  228. MODRTOT = MODTOT ;
  229. MATRTOT0 = MATTOT0 ;
  230. 'SI' ICONDU;
  231. MODRCOND= MODCOND ;
  232. MATRCOND= MATCOND0;
  233. 'FINS';
  234. 'SI' ICONVE;
  235. MODRCON = MODCONV ;
  236. MATRCON = MATCONV0;
  237. 'FINS';
  238. 'SI' IADVEC;
  239. MODRADV = MODADV ;
  240. MATRADV = MATADV0 ;
  241. 'FINS';
  242. 'SI' IRADI;
  243. MODRRAY = MODRAY ;
  244. MATRRAY = MATRAY0 ;
  245. 'FINS';
  246. 'SI' ISOURQ;
  247. MODRSOQ = MODSOQ ;
  248. MATRSOQ = MATSOQ0 ;
  249. 'FINS';
  250. 'SI' IPHASE;
  251. MODRPHA = MODPHA ;
  252. MATRPHA = MATPHA0 ;
  253. 'FINS';
  254. 'SI' IMETAL;
  255. MODRMET = MODMET ;
  256. 'SI' PARTLOCA ;
  257. MODCOMP1 = 'PART' 'ARLE' MODRMET (ETAB.'NBPART') ;
  258. 'FINSI' ;
  259. 'FINS';
  260. 'SI' ITHM1;
  261. MODRTHM = MODTHM ;
  262. MATRTHM = MATTHM0 ;
  263. 'FINS';
  264. 'SI' IDIFF;
  265. MODRDIFF = MODDIF ;
  266. MATRDIFF = MATDIF0 ;
  267. 'FINS' ;
  268. 'FINS' ;
  269.  
  270. *------- Initialisation de la liste des temps a calculer ---------
  271. TE_CALC = ETAB.'TEM_CALC' ;
  272. DCALC ='DIME' TE_CALC ;
  273. ITE_IT = 0 ;
  274.  
  275. *-------------- Calcul du premier pas de calcul -------------------
  276. TEM0 = ETAB.'TEMPS0' ;
  277. TEM1 ='EXTR' TE_CALC 1;
  278. T_PAS = TEM1 - TEM0 ;
  279.  
  280. *-------------- modif pour recuperer le temps du dernier pas de calcul
  281. TEM_VAL = TEM1 ;
  282. *-------------- fin de modif -----------------------------------------
  283.  
  284. *---- Initialisation de la valeur du critere de fin d'iteration ---
  285. ZPREC1 = ETAB.'PRECISION' ;
  286.  
  287. *--------- Initialisation du nombre d'iterations autorisees -------
  288. 'SI' ('EXIS' ETAB 'MAXITERATION');
  289. MAX1 = ETAB.'MAXITERATION';
  290. 'SINO';
  291. MAX1 = 150;
  292. 'FINS';
  293.  
  294. *-------- Recherche des variables selon lesquelles varie ----------
  295. *---------------------- le champ de materiau ----------------------
  296. * Si couplage, on regarde les caracteristiques variables du modele complet
  297. si ETAB.'COUPLAGE' ;
  298. LVAREX = 'EXTR' (ETAB.'CARACTERISTIQUES') 'DEVA';
  299. sino ;
  300. LVAREX = 'EXTR' MATTOT0 'DEVA';
  301. fins ;
  302. ILVA = 'DIME' LVAREX;
  303. IMATCHG = 'EXIS' LVAREX 'TEMP' ;
  304. NONLIN1 = (('EGA' ILVA 1) 'ET' ('NON' IMATCHG))
  305. 'OU' ('>EG' ILVA 2)
  306. 'OU' ETAB.'PROCEDURE_PARATHER' ;
  307. 'SI' NONLIN1 ;
  308. 'SI' ETAB.'COUPLAGE';
  309. CONTI=PRECED.'CONTINUATION';
  310. *--- Creation du MCHAML contenant tous les parametres mecaniques --
  311. *-------- Ces parametres sont determines au temps TEM0 ------------
  312. *----- Ces parametres sont projetes sur le maillage thermique -----
  313. DEP1 ='CHAN' 'CHAM' CONTI.'DEPLACEMENTS' MODMEC 'NOEUD';
  314. CON1 ='CHAN' 'NOEUD' CONTI.'CONTRAINTES' MODMEC ;
  315. CHAR1_ME = DEP1 'ET' CON1 ;
  316. * 'SI' ('EXIS' CONTI 'VITESSES');
  317. * VIT1 = 'CHAN' 'CHAM' CONTI.'VITESSES' MODMEC 'NOEUD';
  318. * CHAR1_ME = CHAR1_ME 'ET' VIT1;
  319. * 'FINS';
  320. 'SI' ('EXIS' CONTI 'VARIABLES_INTERNES');
  321. VAR1 ='CHAN' 'NOEUD' CONTI.'VARIABLES_INTERNES' MODMEC ;
  322. CHAR1_ME = CHAR1_ME 'ET' VAR1 ;
  323. * DEF1 ='CHAN' 'NOEUD' ETAB.'DEI_EST' MODMEC ;
  324. * CHAR1_ME = CHAR1_ME 'ET' VAR1 'ET' DEF1;
  325. 'FINS';
  326. 'SI' ETAB.'PROJECTION';
  327. CHAR_ME ='CHAN' 'CHAM' ('PROI' MAILTOT CHAR1_ME) MODRTOT 'RIGIDITE';
  328. 'SINO';
  329. CHAR_ME ='CHAN' 'CHAM' ('CHAN' 'CHPO' CHAR1_ME MODMEC) MODRTOT 'RIGIDITE';
  330. 'FINS';
  331. 'FINS';
  332. 'FINS';
  333.  
  334. *-----------recherche des composantes de materiau qui varient et initialisation
  335. * --------- des logiques de recalcul des matrices.-------------------
  336. LCOREX = 'EXTR' MATTOT0 'COVA' ;
  337.  
  338. IKDEP =('EXIS' LCOREX 'K ') 'OU' ('EXIS' LCOREX 'K1 ') 'OU' ('EXIS' LCOREX 'K11 ') ;
  339. ICDEP =('EXIS' LCOREX 'RHO ') 'OU' ('EXIS' LCOREX 'C ') 'OU' ('EXIS' LCOREX 'M ') ;
  340. IECDEP = 'EXIS' LCOREX 'H ' ;
  341. IKDDEP =('EXIS' LCOREX 'KD ') 'OU' ('EXIS' LCOREX 'KD1 ') 'OU' ('EXIS' LCOREX 'KD11') ;
  342. ICDDEP = 'EXIS' LCOREX 'CDIF' ;
  343. IKVITE = 'EXIS' LCOREX 'VITE' ;
  344.  
  345. ICPCONS = ETAB.'CAPACONST' ;
  346. IKCONS = ETAB.'CONDCONST' ;
  347.  
  348. 'SI' ITHM1 ;
  349. IKDEP= VRAI ;
  350. ICDEP= VRAI ;
  351. 'FINS' ;
  352.  
  353. *------- Cas particulier ou il n'y a pas de CAPACITE ------------------*
  354. *- Ceci correspond au cas ou, au moins, une des 2 composantes 'C ' -*
  355. *- et/ou 'RHO ' n'est pas definie dans le materiau (facultatives). -*
  356. LCOMPM = 'EXTR' MATTOT0 'COMP' ;
  357.  
  358. ICAPATH = ETAB.'FOR_THER' 'ET' ('EXIS' LCOMPM 'C ') 'ET' (('EXIS' LCOMPM 'RHO ') 'OU' ('EXIS' LCOMPM 'M '));
  359. ICAPADI = IDIFF 'ET' ('EXIS' LCOMPM 'CDIF') ;
  360. ICAPA = ICAPATH 'OU' ICAPADI 'OU' ICPCONS ;
  361. 'SI' ('NON' ICAPATH) ;
  362. ICDEP = FAUX ;
  363. 'FINS' ;
  364. 'SI' ('NON' ICAPADI) ;
  365. ICDDEP = FAUX ;
  366. 'FINS' ;
  367.  
  368. *---------- Initialisation de l'unite de temperature --------------
  369. 'SI' ('EXIS' ETAB 'CELSIUS');
  370. ICELS = ETAB.'CELSIUS' 'ET' ETAB.'FOR_THER' ;
  371. 'SINO';
  372. ICELS = FAUX;
  373. 'FINS';
  374.  
  375. *---------- Initialisation de la temperature de reference----------
  376. 'SI' ICELS;
  377. TREF = 273.15 ;
  378. 'SINO';
  379. TREF = 0. ;
  380. 'FINS';
  381. ETAB.'TREF'=TREF;
  382.  
  383. *----------- Initialisation du champ de temperatures ------------
  384. CHTREF ='VIDE' 'CHPOINT' / 'DIFFUS' ;
  385. U0 ='VIDE' 'CHPOINT' / 'DIFFUS' ;
  386. 'SI' ETAB.'FOR_THER' ;
  387. U0 = ETAB.'TET1';
  388. 'SI' ICELS ;
  389. CHTREF =(0. * U0) + TREF ;
  390. 'FINS' ;
  391. 'FINS' ;
  392.  
  393. 'SI' IDIFF ;
  394. U0 = U0 'ET' ETAB.'CO1';
  395. 'FINS' ;
  396.  
  397. ETAB.'THER_COURANT' = U0 ;
  398. U0MTREF = U0 ;
  399. U0 = U0 + CHTREF ;
  400.  
  401. *---------- y-a-t-il definition des coefficients lambda ----------
  402. * En regime stationnaire (sans CAPACITE), ou en presence d'advection on force LAMBDA1 a 1.
  403. * SINO on ne resout pas le bon probleme !
  404. 'SI' (IADVEC 'OU' ('NON' ICAPA)) ;
  405. LAMBDA1 = 1. ;
  406. 'SINO';
  407. 'SI' ( 'EXIS' ETAB 'RELAXATION_THETA' ) ;
  408. LAMBDA1 = ETAB.'RELAXATION_THETA' ;
  409. 'SINO' ;
  410. LAMBDA1 = 1. ;
  411. 'FINS' ;
  412. 'FINS' ;
  413. ETAB.'RELAXATION_THETA' = LAMBDA1 ;
  414.  
  415.  
  416. 'SI' ('<EG' (LAMBDA1 - 1.0) 1.E-6) ;
  417. I_LAMBD1 = 3 ;
  418. 'SINO' ;
  419. 'SI' ('<EG' (LAMBDA1 - 0.5) 1.E-6) ;
  420. I_LAMBD1 = 2 ;
  421. 'SINO' ;
  422. 'SI' ('<EG' LAMBDA1 1.E-6) ;
  423. I_LAMBD1 = 1 ;
  424. 'SINO' ;
  425. I_LAMBD1 = 0 ;
  426. 'FINS' ;
  427. 'FINS' ;
  428. 'FINS' ;
  429.  
  430. *---- En cas de rayonnement : recalcul des facteurs de forme -----------
  431. *---- Si geometrie modifiee ou s'il n'existent pas, --------------
  432. 'SI' IRADI ;
  433. ETAB.'CHMAT' = MATTOT ;
  434. IAPPEL = 2 ;
  435. t_rayo = PAS_RAYO PRECED tem1 IAPPEL ;
  436. t_rayo ='MOT' 'INCONNU' ;
  437. 'FINS' ;
  438.  
  439. ********************************************************************
  440. * Rappel des donnees du probleme *
  441. ********************************************************************
  442. 'SI' ETAB.'MESSTHER' ;
  443. ETAB.'MESSTHER' = FAUX ;
  444. 'SI' ICAPA ;
  445. 'MESS' ' *** Probleme TRANSITOIRE ' ;
  446. 'SINO' ;
  447. 'MESS' ' *** Probleme STATIONNAIRE' ;
  448. 'FINS' ;
  449.  
  450. 'SI' ('EGA' ETAB.'PROCEDURE_THERMIQUE' 'NONLINEAIRE');
  451. 'MESS' ' *** Methode a un pas de temps (theta methode) ' ;
  452. 'SINO';
  453. 'MESS' ' *** Methode a deux pas de temps (DUPONT2) ' ;
  454. 'MESS' ' > Calcul du premier pas par la theta-methode ' ;
  455. 'FINS';
  456.  
  457. 'MESS' ' *** Donnees utilisees du probleme :';
  458. 'MESS' ' > Champ de temperatures a l instant : ' TEM0;
  459. 'SI' NONLIN1;
  460. 'MESS' ' > K,C ou H variables.';
  461. 'SINO';
  462. 'MESS' ' > K,C et H constants.';
  463. 'FINS';
  464.  
  465. 'SI' ICELS;
  466. 'MESS' ' > Unite de temperature utilisee : Celsius.';
  467. 'SINO';
  468. 'MESS' ' > Unite de temperature utilisee : Kelvin.';
  469. 'FINS';
  470.  
  471. 'SI' ('EGA' ETAB.'PROCEDURE_THERMIQUE' 'NONLINEAIRE');
  472. 'MESS' ' > Coefficient de relaxation : ' LAMBDA1;
  473. 'SINO';
  474. 'MESS' ' > Coefficient de relaxation pour la theta-methode : ' LAMBDA1 ;
  475. 'MESS' ' > Coefficient de relaxation pour la methode Dupont2 : ' ETAB.'RELAXATION_DUPONT' ;
  476. 'MESS' ' > Coefficient de sous-relaxation : ' LAMBDA2;
  477. 'FINS';
  478.  
  479. 'SI' (IPHASE 'OU' ISOLU);
  480. 'SI' (IPHASE 'ET' ISOLU) ;
  481. 'MESS' ' > Changement de phase ''PARFAIT'' et ''SOLUBILITE'' prevus.';
  482. 'SINO' ;
  483. 'SI' IPHASE ;
  484. 'MESS' ' > Changement de phase ''PARFAIT'' prevu.';
  485. 'SINO';
  486. 'MESS' ' > Changement de phase ''SOLUBILITE'' prevu.';
  487. 'FINS';
  488. 'FINS' ;
  489. 'SINO';
  490. 'MESS' ' > Pas de changement de phase prevu.';
  491. 'FINS';
  492.  
  493. 'SI' ('NEG' ETAB.'BLOCAGES_THERMIQUES' 'INCONNU');
  494. 'MESS' ' > Presence de temperatures imposees.';
  495. 'SINO';
  496. 'MESS' ' > Pas de temperatures imposees.';
  497. 'FINS';
  498.  
  499. 'SI' ('EXIS' ETAB.'CHARGEMENT' 'Q ');
  500. 'MESS' ' > Presence de flux.';
  501. 'SINO';
  502. 'MESS' ' > Pas de presence de flux.';
  503. 'FINS';
  504.  
  505. 'SI' ICONDU;
  506. 'MESS' ' > Presence de conduction';
  507. 'SINO';
  508. 'MESS' ' > Pas de presence de conduction ';
  509. 'FINS';
  510.  
  511. 'SI' ICONVE;
  512. 'MESS' ' > Presence de convection.';
  513. 'SINO';
  514. 'MESS' ' > Pas de terme de convection.';
  515. 'FINS';
  516.  
  517. 'SI' IRADI;
  518. 'MESS' ' > Presence de rayonnement.';
  519. 'SINO';
  520. 'MESS' ' > Pas de presence de rayonnement.';
  521. 'FINS';
  522.  
  523. 'SI' IADVEC;
  524. 'MESS' ' > Presence d advection.';
  525. 'SINO';
  526. 'MESS' ' > Pas de presence d advection';
  527. 'FINS';
  528.  
  529. 'SI' ICPCONS ;
  530. 'MESS' ' > Presence d une capacite constante.';
  531. 'FINS';
  532.  
  533. 'SI' IKCONS ;
  534. 'MESS' ' > Presence d une conductivite constante.';
  535. 'FINS';
  536. 'MESS' ' ';
  537. 'FINS';
  538.  
  539. ********************************************************************
  540. * Boucle sur les pas de temps *
  541. ********************************************************************
  542. 'REPETER' BOUC_SO ;
  543. ITE3 = &BOUC_SO;
  544. TEMM = (LAMBDA1 * TEM1) + ((1 - LAMBDA1) * TEM0);
  545. TH_COUR = U0MTREF ;
  546. ETAB.'THER_COURANT' = TH_COUR ;
  547.  
  548. 'MESS' ' ';
  549. 'MESS' '---- Calcul du champ de temperature au temps :'TEM1;
  550. 'MESS' ' ';
  551. *---- initialisation du materiau
  552. 'SI' (NONLIN1 'OU' IMATCHG) ;
  553. *-------------- Cas de la temperature et du temps ---------------
  554. CHVAR ='CHAN' 'CHAM' TH_COUR MODRTOT 'RIGIDITE' ;
  555. CHVAR2 ='MANU' 'CHML' MODRTOT 'TEMP' TEMM 'RIGIDITE' ;
  556. CHVAR = CHVAR 'ET' CHVAR2 ;
  557. *-------------- Cas des variables mecaniques --------------------
  558. *---- si il y a couplage on les "rentre" toutes dans CHVAR ------
  559. 'SI' ETAB.'COUPLAGE';
  560. CHVAR = CHVAR 'ET' CHAR_ME;
  561. 'FINS';
  562. *-------------- Mise a jour des variables externes --------------
  563. 'SI' ETAB.'PROCEDURE_PARATHER' ;
  564. PARATHER PRECED TEMM ;
  565.  
  566. * Remise en place des MCHAML sur les bons MMODELS car on suppose
  567. * que l'utilisateur ne connais pas tous les indices de ETAB.
  568. MATRTOT0 ='REDU' ETAB.'MAT_TOT' MODRTOT ;
  569. 'SI' ICONDU;
  570. MATRCOND='REDU' ETAB.'MAT_COND' MODRCOND;
  571. 'FINS';
  572. 'SI' ICONVE;
  573. MATRCON ='REDU' ETAB.'MAT_CON' MODRCON ;
  574. 'FINS';
  575. 'SI' IADVEC;
  576. MATRADV ='REDU' ETAB.'MAT_ADV' MODRADV ;
  577. 'FINS' ;
  578. 'SI' IRADI;
  579. MATRRAY ='REDU' ETAB.'MAT_RAY' MODRRAY ;
  580. 'FINS';
  581. 'SI' IPHASE;
  582. MATRPHA ='REDU' ETAB.'MAT_PHA' MODRPHA ;
  583. 'FINS';
  584. 'SI' ITHM1;
  585. MATRTHM ='REDU' ETAB.'MAT_THM' MODRTHM ;
  586. 'FINS' ;
  587. 'SI' IDIFF ;
  588. MATRDIFF='REDU' ETAB.'MAT_DIF' MODRDIFF;
  589. 'FINS' ;
  590. 'FINS' ;
  591.  
  592. *-------------- Cas des autres variables --------------------
  593. 'REPETER' BOU_NU2 ILVA;
  594. ITE2 = &BOU_NU2;
  595. MOTI = 'EXTR' LVAREX ITE2;
  596. 'SI' ('EXIS' LPRIM2 MOTI) ;
  597. 'ITER' BOU_NU2 ;
  598. 'FINS' ;
  599. 'SI' ETAB.'COUPLAGE';
  600. 'SI' ('EXIS' CHAR_ME MOTI );
  601. 'SI' ('NEG' LAMBDA1 0.);
  602. 'MESS'
  603. '*** Le materiau thermique depend de parametres mecaniques';
  604. 'ERRE' '*** Il faut initialiser SOUS-RELAXATION a 0.';
  605. 'FINS';
  606. 'ITERER' BOU_NU2;
  607. 'FINS';
  608. 'FINS';
  609. 'SI' ('EXIS' ETAB.'CHARGEMENT' MOTI);
  610. CH_1 = 'TIRE' ETAB.'CHARGEMENT' MOTI TEMM;
  611. TYP_1 = 'TYPE' CH_1;
  612.  
  613. * on suppose que les chargements qui ne sont pas d'origine mecanique
  614. * sont appuyes sur le maillage thermique.
  615. 'SI' ('EGA' TYP_1 'CHPOINT');
  616. CHVAR2 = 'CHAN' 'CHAM' CH_1 MODRTOT 'RIGIDITE';
  617. 'SINO';
  618. CHVAR2 = 'CHAN' 'RIGIDITE' CH_1 MODRTOT;
  619. 'FINS';
  620. CHVAR = CHVAR + CHVAR2;
  621. 'SINO';
  622. 'MESS' '*** Il manque une donnee du chargement ***';
  623. 'MESS' '*** Variable concernee : 'MOTI;
  624. 'ERREUR' ' ';
  625. 'FINS';
  626. 'FIN' BOU_NU2;
  627. 'FINS';
  628.  
  629. *----------------- Calcul du champ de materiau a T= U0 -----------------
  630. 'SI' (NONLIN1 'OU' IMATCHG) ;
  631. CHVAR1 ='REDU' CHVAR MODRTOT;
  632. CHMAT ='VARI' 'NUAG' MODRTOT MATRTOT0 CHVAR1 'RIGIDITE';
  633. 'SI' (IDIFF 'OU' ETAB.'THE1') ;
  634. CHMAT1 = 'REDU' CHMAT ETAB.'MOD_THE' ;
  635. 'SI' ITHM1 ;
  636. ma ff0 = @MATETHM MODRTHM TH_COUR ;
  637. CHMAT1 = CHMAT1 'ET' ma ;
  638. 'FINS' ;
  639. 'SINO' ;
  640. 'SI' ITHM1 ;
  641. CHMAT1 ff0 = @MATETHM MODRTHM TH_COUR ;
  642. 'FINS' ;
  643. 'FINS' ;
  644.  
  645. 'SINO';
  646. CHMAT = MATRTOT0 ;
  647. 'SI' (IDIFF 'OU' ETAB.'THE1') ;
  648. CHMAT1 = 'REDU' CHMAT ETAB.'MOD_THE' ;
  649. 'SI' ITHM1 ;
  650. ma ff0 = @MATETHM MODRTHM TH_COUR ;
  651. CHMAT1 = CHMAT1 'ET' ma ;
  652. 'FINS' ;
  653. 'SINO' ;
  654. 'SI' ITHM1 ;
  655. CHMAT1 ff0 = @MATETHM MODRTHM TH_COUR ;
  656. 'FINS' ;
  657. 'FINS' ;
  658. 'FINS';
  659.  
  660. * On conserve le materiau instancie CHMAT dans ETAB (Pour PAS_RAYO...)
  661. ETAB.'CHMAT' ='REDU' CHMAT MODTOT ;
  662.  
  663. MATRTOT ='REDU' CHMAT1 MODRTOT ;
  664. 'SI' ICONDU;
  665. MATRCOND='REDU' CHMAT1 MODRCOND;
  666. 'FINS';
  667. 'SI' ICONVE;
  668. MATRCON ='REDU' CHMAT1 MODRCON;
  669. 'FINS';
  670. 'SI' IRADI;
  671. MATRRAY ='REDU' CHMAT1 MODRRAY;
  672. 'FINS';
  673. 'SI' IADVEC;
  674. MATRADV ='REDU' CHMAT1 MODRADV;
  675. 'FINS' ;
  676. 'SI' ISOURQ ;
  677. MATRSOQ ='REDU' CHMAT1 MODRSOQ ;
  678. 'FINS' ;
  679. 'SI' IPHASE;
  680. MATRPHA ='REDU' CHMAT1 MODRPHA;
  681. 'FINS';
  682. 'SI' ITHM1;
  683. MATRTHM ='REDU' CHMAT1 MODRTHM ;
  684. 'FINS' ;
  685. 'SI' IDIFF ;
  686. MATRDIFF='REDU' CHMAT1 MODRDIFF;
  687. 'FINS' ;
  688.  
  689. *--- Calcul des matrices d'iteration On les recalcule a cause de perso
  690. IKDETO=FAUX;IBCONS=FAUX;
  691. 'SI' ICONDU ;
  692. KDETO ='COND' MODRCOND MATRCOND;
  693. IKDETO= VRAI;
  694. 'FINS';
  695. 'SI' ITHM1;
  696. KDCC ='COND' MODRTHM MATRTHM;
  697. 'SI' IKDETO;
  698. KDETO=KDETO 'ET' KDCC;
  699. 'SINO';
  700. KDETO=KDCC;
  701. IKDETO=VRAI;
  702. 'FINS';
  703. 'FINS';
  704. 'SI' (('NON' IKDEP) 'ET' IKDETO);
  705. BCONS= KDETO;
  706. IBCONS=VRAI;
  707. 'FINS';
  708.  
  709. 'SI' ICONVE;
  710. KDCCG= 'COND' MODRCON MATRCON;
  711. 'SI' IKDETO;
  712. KDETO=KDETO 'ET' KDCCG;
  713. 'SINO';
  714. KDETO=KDCCG;
  715. IKDETO=VRAI;
  716. 'FINS';
  717. 'SI' ('NON' IECDEP);
  718. 'SI' IBCONS ;
  719. BCONS=BCONS 'ET' KDCCG;
  720. 'SINO';
  721. BCONS=KDCCG;
  722. IBCONS=VRAI;
  723. 'FINS';
  724. 'FINS';
  725. 'FINS';
  726.  
  727. 'SI' IADVEC;
  728. * la conduction de du model d'advection est deja calculer grace a icondu
  729. 'SI' ('NON' IKVITE);
  730. KDCC2= 'ADVE' MODRADV MATRADV;
  731. 'SI' IKDETO;
  732. KDETO=KDETO 'ET' KDCC2;
  733. 'SINO';
  734. KDETO= KDCC2 ;
  735. IKDETO=VRAI;
  736. 'FINS';
  737. 'SI' IBCONS;
  738. BCONS = BCONS 'ET' KDCC2;
  739. 'SINO';
  740. BCONS= KDCC2;
  741. IBCONS=VRAI;
  742. 'FINS';
  743. 'FINS';
  744. 'FINS';
  745.  
  746. 'SI' IDIFF;
  747. KDEDI = 'COND' MODRDIFF MATRDIFF ;
  748. 'SI' IKDETO;
  749. KDETO=KDETO 'ET' KDEDI;
  750. 'SINO';
  751. KDETO=KDEDI;
  752. IKDETO=VRAI;
  753. 'FINS';
  754. 'SI' ('NON' IKDDEP);
  755. 'SI' IBCONS;
  756. BCONS=BCONS 'ET' KDEDI;
  757. 'SINO';
  758. BCONS = KDEDI;
  759. IBCONS= VRAI ;
  760. 'FINS';
  761. 'FINS';
  762. 'FINS' ;
  763.  
  764. 'SI' IKCONS ;
  765. 'SI' IKDETO ;
  766. KDETO = KDETO 'ET' ETAB.'CONDUCTIVITE_CONSTANTE' ;
  767. 'SINO' ;
  768. KDETO = ETAB.'CONDUCTIVITE_CONSTANTE' ; KDETO = VRAI ;
  769. 'FINS' ;
  770. 'SI' IBCONS;
  771. BCONS=BCONS 'ET' ETAB.'CONDUCTIVITE_CONSTANTE' ;
  772. 'SINO';
  773. BCONS= ETAB.'CONDUCTIVITE_CONSTANTE' ; IBCONS=VRAI ;
  774. 'FINS';
  775. 'FINS' ;
  776. *
  777. AA = KDETO * LAMBDA1 ;
  778. CCONS = AA ;
  779. BB = KDETO ;
  780. 'SI' ('NON' IBCONS) ;
  781. BCONS = EXTR KDETO 'RIGI' 'MULT'; IBCONS=VRAI;
  782. 'FINS';
  783. *
  784. 'SI' ICAPA ;
  785. ETAB.'CTHE_COURANTE' = 'VIDE' 'RIGIDITE'/'RIGIDITE' ;
  786. ETAB.'CDIF_COURANTE' = 'VIDE' 'RIGIDITE'/'RIGIDITE' ;
  787. 'SI' ICAPATH ;
  788. 'SI' ICONDU;
  789. CDET = 'CAPA' MODRCOND MATRCOND ;
  790. ETAB.'CTHE_COURANTE' = ETAB.'CTHE_COURANTE' et CDET ;
  791. CDET = CDET / T_PAS;
  792. AA = CDET 'ET' AA ;
  793. 'SI' ICDEP ; CCONS = CDET 'ET' CCONS ; 'FINS' ;
  794. 'FINS';
  795. 'SI' ITHM1;
  796. CDET = 'CAPA' MODRTHM MATRTHM ;
  797. ETAB.'CTHE_COURANTE' = ETAB.'CTHE_COURANTE' et CDET ;
  798. CDET = CDET / T_PAS;
  799. AA = CDET 'ET' AA ;
  800. 'SI' ICDEP; CCONS = CDET 'ET' CCONS ; 'FINS' ;
  801. 'FINS';
  802. 'FINS';
  803. 'SI' ICAPADI ;
  804. CDET = 'CAPA' MODRDIFF MATRDIFF ;
  805. ETAB.'CDIF_COURANTE' = ETAB.'CDIF_COURANTE' et CDET ;
  806. CDET = CDET / T_PAS;
  807. AA = CDET 'ET' AA ;
  808. 'SI' ICDDEP ; CCONS = CDET 'ET' CCONS ; 'FINS' ;
  809. 'FINS' ;
  810. ETAB.'CAPA_COURANTE' = ETAB.'CTHE_COURANTE' et ETAB.'CDIF_COURANTE' ;
  811. 'SI' ICPCONS ;
  812. ETAB.'CAPA_COURANTE' = ETAB.'CAPA_COURANTE' et ETAB.'CAPACITE_CONSTANTE' ;
  813. CDET = ETAB.'CAPACITE_CONSTANTE' / T_PAS ;
  814. AA = CDET 'ET' AA ;
  815. 'SI' ICDDEP ; CCONS = CDET 'ET' CCONS ; 'FINS' ;
  816. 'FINS' ;
  817. 'FINS' ;
  818.  
  819. *----------------- y-a-t-il des blocages ---------------------------
  820. * on le laisse ici car perso peut l'avoir modifie de memepour les autres
  821. * matrices
  822. IKBLOC=FAUX;
  823. 'SI' ETAB.'FOR_THER' ;
  824. IKBLOC=VRAI;
  825. MAT_BLOC = ETAB.'BLOCAGES_THERMIQUES' ;
  826. 'SI' ('EGA' MAT_BLOC 'INCONNU') ;
  827. MAT_BLOC = 'EXTR' KDETO 'RIGI' 'MULT' ;
  828. 'FINS' ;
  829. MAT_BLTH = MAT_BLOC ;
  830. 'FINS' ;
  831.  
  832. 'SI' IDIFF ;
  833. MAT_BLDI = ETAB.'BLOCAGES_DIFFUSIONS' ;
  834. 'SI' ('EGA' MAT_BLDI 'INCONNU') ;
  835. MAT_BLDI = 'EXTR' KDEDI 'RIGI' 'MULT' ;
  836. 'FINS' ;
  837. 'SI' IKBLOC;
  838. MAT_BLOC = MAT_BLOC 'ET' MAT_BLDI;
  839. 'SINO';
  840. MAT_BLOC = MAT_BLDI;
  841. 'FINS';
  842. 'FINS' ;
  843.  
  844. AA = AA 'ET' MAT_BLOC ;
  845. MAT_CHPO = MAT_BLOC * U0 ;
  846. SECON = MAT_CHPO ;
  847.  
  848. *------------------- y-a-t-il des temperatures imposees--------
  849. ICHUIMP=FAUX;
  850. 'SI' ETAB.'FOR_THER' ;
  851. 'SI' ('EXIS' ETAB.'CHARGEMENT' 'TIMP') ;
  852. CH_UIMP = 'TIRE' ETAB.'CHARGEMENT' 'TIMP' TEM1 ;
  853. 'SI' ICELS ;
  854. CH_UIMP = CH_UIMP + TREF ;
  855. 'FINS' ;
  856. ICHUIMP=VRAI;
  857. 'FINS';
  858. 'FINS';
  859. 'SI' IDIFF;
  860. 'SI' ('EXIS' ETAB.'CHARGEMENT' 'CIMP') ;
  861. CH_CI = 'TIRE' ETAB.'CHARGEMENT' 'CIMP' TEM1 ;
  862. 'SI' ICHUIMP;
  863. CH_UIMP= CH_UIMP 'ET' CH_CI;
  864. 'SINO';
  865. CH_UIMP=CH_CI;
  866. 'FINS';
  867. ICHUIMP=VRAI;
  868. 'FINS';
  869. 'FINS';
  870. 'SI' ('NON' ICHUIMP);
  871. CH_UIMP='DEPI' MAT_BLOC 0.;
  872. 'FINS';
  873. MAT_CHPO = CH_UIMP - MAT_CHPO ;
  874. SECON = CH_UIMP - SECON ;
  875.  
  876. *---------- Si convection on cherche TE ------------
  877. *----------calcul convection second membre --------------
  878. 'SI' ICONVE;
  879. 'SI' ('EXIS' ETAB.'CHARGEMENT' 'TECO' ) ;
  880. CH_TEC ='TIRE' ETAB.'CHARGEMENT' 'TECO' TEMM ;
  881. 'SI' ICELS ;
  882. CH_TEC = CH_TEC + TREF ;
  883. 'FINS' ;
  884. CHAL_TEI ='CONV' MODRCON MATRCON CH_TEC ;
  885.  
  886. MAT_CHPO = MAT_CHPO + CHAL_TEI;
  887. 'SI' ('NON' IECDEP) ;
  888. SECON = SECON + CHAL_TEI;
  889. 'FINS' ;
  890. 'SINO' ;
  891. 'SI' ('EXIS' MATRCON 'TC') ;
  892. CH_TEC = 'EXCO' MATRCON 'TC' 'T' ;
  893. CH_TEC = 'CHAN' 'CHPO' MODRCON CH_TEC ;
  894. 'SI' ICELS ;
  895. CH_TEC = CH_TEC + TREF ;
  896. 'FINS' ;
  897. CHAL_TEI ='CONV' MODRCON MATRCON CH_TEC ;
  898. MAT_CHPO = MAT_CHPO + CHAL_TEI;
  899. 'SI' ('NON' IECDEP) ;
  900. SECON = SECON + CHAL_TEI;
  901. 'FINS' ;
  902. 'FINS' ;
  903. 'FINS';
  904. 'FINS';
  905.  
  906. *--------- SI terme de flux -------------------------------
  907.  
  908. 'SI' ITHM1 ;
  909. ma FTGC = @MATETHM MODRTHM TH_COUR ;
  910. MAT_CHPO = MAT_CHPO - FTGC;
  911. SECON = SECON - FTGC;
  912. 'FINS';
  913.  
  914. 'REPE' Boucle NBDUA1 ;
  915. MOTCOMP = 'EXTR' LDUAL1 &Boucle ;
  916. 'SI' ('EXIS' ETAB.'CHARGEMENT' MOTCOMP) ;
  917. CH_DUA = 'TIRE' ETAB.'CHARGEMENT' MOTCOMP TEMM ;
  918. MAT_CHPO = MAT_CHPO + CH_DUA ;
  919. SECON = SECON + CH_DUA ;
  920. 'FINS' ;
  921. 'FIN' Boucle ;
  922.  
  923. *------------- calcul des termes rayonnement --------------------------
  924.  
  925. 'SI' IRADI ;
  926. IAPPEL = 3 ;
  927. t_rayo = PAS_RAYO PRECED TEMM IAPPEL ;
  928. 'SI' ('EXIS' t_rayo 'ADDI_MATRICE') ;
  929. RIG_RAYO = LAMBDA1 '*' t_rayo.'ADDI_MATRICE';
  930. ccons = RIG_RAYO 'ET' ccons ;
  931. aa = RIG_RAYO 'ET' aa ;
  932. bb = t_rayo.'ADDI_MATRICE' 'ET' bb;
  933. 'FINS';
  934. 'SI' ('EXIS' t_rayo 'ADDI_SECOND') ;
  935. MAT_CHPO = MAT_CHPO 'ET' t_rayo.'ADDI_SECOND' ;
  936. 'FINS';
  937. 'SI' ('EXIS' t_rayo 'RAYO_VALEUR');
  938. STAB.'RAYONNEMENT' = t_rayo.'RAYO_VALEUR' ;
  939. 'FINS';
  940. t_rayo = 'MOT' 'INCONNU' ;
  941. 'FINS';
  942.  
  943. *------------- calcul des sources de chaleur --------------------
  944.  
  945. 'SI' ISOURQ ;
  946. CHQ1 = SOUR MODRSOQ MATRSOQ ;
  947. MAT_CHPO = MAT_CHPO 'ET' CHQ1 ;
  948. SECON = SECON 'ET' CHQ1 ;
  949. 'FINS';
  950.  
  951. *------------- en cas de procedure utilisateur CHARTHER --------------
  952.  
  953. 'SI' etab.'PROCEDURE_CHARTHER' ;
  954. taat = charther preced temm ;
  955. 'SI' ('EXIS' taat 'ADDI_MATRICE') ;
  956. mat_z = lambda1 '*' taat.'ADDI_MATRICE' ;
  957. ccons = mat_z 'ET' ccons ;
  958. aa = mat_z 'ET' aa ;
  959. bb = taat.'ADDI_MATRICE' 'ET' bb ;
  960. 'FINS';
  961. 'SI' ('EXIS' taat 'ADDI_SECOND') ;
  962. MAT_CHPO = MAT_CHPO 'ET' taat.'ADDI_SECOND' ;
  963. 'FINS';
  964. oublier taat ;
  965. 'FINS';
  966.  
  967. * Pour le Changement de phase 'PARFAIT' : Blocages et Jeux
  968. 'SI' IPHASE;
  969. AA = AA 'ET' BLOPHA ;
  970. Q_LAT_INT=('SOUR' MODPHA ('EXCO' MATTOT0 'DUAL') MATTOT0 'ELEM') / T_PAS ;
  971. JEU_PHAS = 'PHAJ' MODTOT MATTOT0 U0MTREF ;
  972. MAT_CHPO = MAT_CHPO + JEU_PHAS ;
  973. SECON = SECON + JEU_PHAS ;
  974. 'FINS';
  975.  
  976. * Pour le Changement de phase 'SOLUBILITE' : Blocages et Jeux
  977. 'SI' ISOLU;
  978. AA = AA 'ET' BLOSOL ;
  979. JEU_SOL = 'PHAJ' MODTOT MATRTOT0 U0MTREF ;
  980. MAT_CHPO = MAT_CHPO + JEU_SOL ;
  981. SECON = SECON + JEU_SOL ;
  982. 'FINS';
  983. *--------------- Calcul du premier residu --------------------------
  984.  
  985. B3 = BB * U0 ;
  986. MAT_CHPO = MAT_CHPO - B3 ;
  987. ********************************************************************
  988. * Boucle sur les iterations *
  989. ********************************************************************
  990. * on peut mettre n'importe quoi c'est pour
  991. * ne pas faire de tests dans la boucle
  992. ZNACCE = 2 ;
  993. ZITAC = 0 ;
  994. IT = 0 ;
  995. ITDEP = 3 ;
  996. 'SI' IPHASE;
  997. ITDEP=5;
  998. 'FINS';
  999. ITACC = 0 ;
  1000. ACFP1 = MAT_CHPO * 0. ;
  1001. ACFP2 = ACFP1 ;
  1002. MAT_CHPP= ACFP1 ;
  1003. ACFP3 = ACFP1 ;
  1004. ACFEP1 = ACFP1 ;
  1005. ACFEP2 = ACFP1 ;
  1006. CORREC = 0.;
  1007. U0 ='ENLE' U0 'LX';
  1008. U1 = U0 ;
  1009. DU0 ='VIDE' 'CHPOINT' / 'DIFFUS' ;
  1010. DU1 ='VIDE' 'CHPOINT' / 'DIFFUS' ;
  1011. EXCPHA ='VIDE' 'CHPOINT' / 'DISCRET';
  1012.  
  1013. 'REPETER' BOU_IT1;
  1014. * IT est le compteur de BOU_IT1, ITACC doit etre =< 0 pour qu'on accelere
  1015. IT = IT + 1 ;
  1016. ITACC = ITACC - 1 ;
  1017. ZITAC = ZITAC + 1 ;
  1018. ITE_IT = &BOU_IT1;
  1019.  
  1020. * Nombre max d'iteration autorise est-il depasse?
  1021. 'SI' (ITE_IT '>' MAX1);
  1022. 'SAUTER' 1 'LIGNE' ;
  1023. 'MESS' '*** LE CALCUL N A PAS CONVERGE ***';
  1024. ERR_VAL= VRAI;
  1025. 'QUIT' BOUC_SO ;
  1026. 'FINS';
  1027. AUTORISA= VRAI;
  1028. BB = BCONS;CC=CCONS;
  1029. *
  1030. *---------------------------------------------------------------------
  1031. * La force motrice de l'iteration est fixee: RESIDU
  1032. * on va calculer un nouveau champ de deplacement
  1033. *
  1034. * calcul de l'increment de l'increment de deplacement zdep1
  1035. * par resolution lineaire
  1036. *
  1037. *
  1038. * acceleration de convergence
  1039. *
  1040. CORRECP = CORREC;
  1041. CORREC = 0;
  1042. ACFP0 = (MAT_CHPO - MAT_CHPP) 'ENLE' FLX ;
  1043. ACFEP0 = ACFP0;
  1044. ACFEP0 = ACFEP0 - CORRECP ;
  1045. 'SI' ('MULT' IT ZNACCE) ;
  1046. 'SI' (IT > ITDEP);
  1047. CORREC = 'ACT3' ACFEP2 ACFEP1 ACFEP0 ACFP3 ACFP2 ACFP1 ACFP0 ;
  1048. MAT_CHPO = MAT_CHPO - CORREC ;
  1049. 'FINS';
  1050. 'FINS';
  1051. MAT_CHPP=MAT_CHPO ;
  1052.  
  1053. SI (IT > 3);
  1054. 'DETR' ACFP3;
  1055. 'DETR' ACFEP2;
  1056. 'FINS';
  1057. ACFP3 = ACFP2 ; ACFP2 = ACFP1 ; ACFP1 = ACFP0 ;
  1058. ACFEP2 = ACFEP1 ; ACFEP1 = ACFEP0 ;
  1059.  
  1060. MENAGE;
  1061.  
  1062. 'SOUC' 0;
  1063. 'SI' IRIB ;
  1064. 'SI' IPHASE;
  1065. DU1 BID BID BID BID = 'RESO' AA MAT_CHPO 'INIB' RIBL_VAL LISE_VAL EXCPHA ;
  1066. 'SINO';
  1067. DU1 BID BID BID = 'RESO' AA MAT_CHPO 'INIB' RIBL_VAL LISE_VAL ;
  1068. 'FINS';
  1069. 'SINO';
  1070. 'SI' IPHASE;
  1071. DU1 = 'RESO' AA MAT_CHPO EXCPHA ;
  1072. 'SINO';
  1073. DU1 = 'RESO' AA MAT_CHPO ;
  1074. * mess ' ITeration ' it;
  1075. * aaa= extr aa rigi mult ; mess ' aaa ' ; list aaa;
  1076. * mess ' mat_chpo' ; list mat_chpo;
  1077. * mess ' du1 ' ; list du1;
  1078. 'FINS';
  1079. 'FINS';
  1080.  
  1081. ZRAID_T = 'EXTR' AA 'CONT';
  1082. 'SI' ('EXIS' ZRAID_T 'NITER');
  1083. RIBL_VAL = ZRAID_T.7;
  1084. LISE_VAL = ZRAID_T.6;
  1085. IRIB=VRAI;
  1086. 'FINS';
  1087.  
  1088.  
  1089. 'SI' IPHASE;
  1090. * Flux generalises ('LX')
  1091. EXCPHA ='EXCP' MODTOT MATPHA0 PROPHA0 Q_LAT_INT U0MTREF DU1 ;
  1092. 'SOUC' 0;
  1093. PROPHA1 ='PROP' MODTOT PROPHA0 Q_LAT_INT DU1 ;
  1094. AUTORISA=(ITE_IT '>' 1) 'ET' ('NON' ('SOUC')) ;
  1095. 'FINS';
  1096.  
  1097.  
  1098. * 'SI' ISOLU;
  1099. ** Flux generalises ('LX')
  1100. * EXCPHA =('EXCS' MODTOT CHMAT U0MTREF DU1) 'ET' EXCPHA ;
  1101. * AUTORISA=(ITE_IT '>' 1) 'ET' ('NON' ('SOUC')) ;
  1102. * 'FINS';
  1103.  
  1104. U1 = U0 + DU1 ;
  1105.  
  1106. 'SI' ('EGA' I_LAMBD1 0) ;
  1107. U_BOU1 = (LAMBDA1 * U1) + ((1 - LAMBDA1) * U0);
  1108. 'SINO' ;
  1109. 'SI' ('EGA' I_LAMBD1 1) ; U_BOU1 = U0 ; 'FINS' ;
  1110. 'SI' ('EGA' I_LAMBD1 2) ; U_BOU1 = 0.5 * (U0 + U1) ; 'FINS' ;
  1111. 'SI' ('EGA' I_LAMBD1 3) ; U_BOU1 = U1 ; 'FINS' ;
  1112. 'FINS' ;
  1113.  
  1114. TH_COUR = U_BOU1 - CHTREF;
  1115. ETAB.'THER_COURANT' = TH_COUR ;
  1116.  
  1117.  
  1118. * Calcul des Proportions de Phases (Modele Metallurgique par T.L. en mai 2018)
  1119. 'SI' IMETAL ;
  1120. CHT0 ='CHAN' 'CHAM' ('EXCO' U0 'T' 'T') MODRMET ;
  1121. CHT1 ='CHAN' 'CHAM' ('EXCO' U1 'T' 'T') MODRMET ;
  1122. CHETPS0 ='MANU' 'CHML' MODRMET 'TEMP' tem0 'STRESSES' ;
  1123. CHETPS1 ='MANU' 'CHML' MODRMET 'TEMP' tem1 'STRESSES' ;
  1124. si PARTLOCA ;
  1125. che11 = assi tous redu (PHA_MET0 'ET' CHETPS0 'ET' CHT0 'ET' CHMAT) MODCOMP1 ;
  1126. che22 = assi tous redu ( CHETPS1 'ET' CHT1 'ET' CHMAT) MODCOMP1 ;
  1127. PHA_MET1 = assi tous 'COMP' MODCOMP1 che11 che22 ;
  1128. PHA_MET1 = (etg PHA_MET1) redu MODRMET ;
  1129. che11 che22 = 0 0 ;
  1130. sinon ;
  1131. PHA_MET1 = 'COMP' MODRMET (PHA_MET0 'ET' CHETPS0 'ET' CHT0 'ET' CHMAT)
  1132. ( CHETPS1 'ET' CHT1 'ET' CHMAT ) ;
  1133. finsi ;
  1134. 'FINS' ;
  1135.  
  1136.  
  1137. *-------- Calcul du critere de fin d'iteration -------------------
  1138. DDU_SANS ='EXCO' ('ABS' (DU1 - DU0)) LPRIM1 ;
  1139. 'SI' LOG_TOL ;
  1140. ZPREC3 ='MAXI' DDU_SANS ;
  1141. 'MESS' ' Iteration numero :' ITE_IT/23 ' Critere absolu :' ZPREC3 ;
  1142. CRITABS ='MAXI' ('MASQ' (DDU_SANS - TTOL) 'SUPE' 0.D0) ;
  1143. 'SI' ((CRITABS '&lt;EG' 0.5) 'ET' AUTORISA) ;
  1144. 'QUITTER' BOU_IT1 ;
  1145. 'FINS' ;
  1146.  
  1147. 'SINO' ;
  1148. CHP_CRIT_REL = DDU_SANS * ((('ABS' U1) + XPETI) ** -1) LPRIM1 LPRIM1 LPRIM1 ;
  1149. ZPREC2 = 'MAXI' 'ABS' CHP_CRIT_REL ;
  1150. 'MESS' ' Iteration numero :' ITE_IT/23 ' Critere relatif :'ZPREC2 ;
  1151. 'SI' ((ZPREC2 '&lt;EG' ZPREC1) 'ET' AUTORISA) ;
  1152. 'QUITTER' BOU_IT1 ;
  1153. 'FINS' ;
  1154. 'FINS' ;
  1155.  
  1156. DU0 = DU1 ;
  1157.  
  1158.  
  1159. *****************************************************************
  1160. * CALCUL DU NOUVEAU RESIDU *
  1161. *****************************************************************
  1162. MAT_CHPO = SECON ;
  1163. * Gestion des variables externes *
  1164. *****************************************************************
  1165.  
  1166. 'SI' NONLIN1 ;
  1167.  
  1168. *-------------- Cas de la temperature et du temps ---------------
  1169.  
  1170.  
  1171. CHVAR ='CHAN' 'CHAM' TH_COUR MODRTOT 'RIGIDITE' ;
  1172. CHVAR2 ='MANU' 'CHML' MODRTOT 'TEMP' TEMM 'RIGIDITE' ;
  1173. CHVAR = CHVAR 'ET' CHVAR2 ;
  1174.  
  1175. *-------------- Cas des variables mecaniques --------------------
  1176. *---- si il y a couplage on les "rentre" toutes dans CHVAR ------
  1177.  
  1178. 'SI' ETAB.'COUPLAGE';
  1179. CHVAR = CHVAR 'ET' CHAR_ME;
  1180. 'FINS';
  1181.  
  1182. *-------------- Mise a jour des variables externes --------------
  1183. 'SI' ETAB.'PROCEDURE_PARATHER' ;
  1184. PARATHER PRECED TEMM ;
  1185.  
  1186. * Remise en place des MCHAML sur les bons MMODELS car on suppose
  1187. * que l'utilisateur ne connais pas tous les indices de ETAB.
  1188. MATRTOT0 ='REDU' ETAB.'MAT_TOT' MODRTOT ;
  1189. 'SI' ICONDU;
  1190. MATRCOND='REDU' ETAB.'MAT_COND' MODRCOND;
  1191. 'FINS';
  1192. 'SI' ICONVE;
  1193. MATRCON ='REDU' ETAB.'MAT_CON' MODRCON ;
  1194. 'FINS';
  1195. 'SI' IADVEC;
  1196. MATRADV ='REDU' ETAB.'MAT_ADV' MODRADV ;
  1197. 'FINS' ;
  1198. 'SI' IRADI;
  1199. MATRRAY ='REDU' ETAB.'MAT_RAY' MODRRAY ;
  1200. 'FINS';
  1201. 'SI' IPHASE;
  1202. MATRPHA ='REDU' ETAB.'MAT_PHA' MODRPHA ;
  1203. MATPHA ='REDU' ETAB.'MAT_PHA' MODPHA ;
  1204. 'FINS';
  1205. 'SI' ITHM1;
  1206. MATRTHM ='REDU' ETAB.'MAT_THM' MODRTHM ;
  1207. 'FINS' ;
  1208. 'SI' IDIFF ;
  1209. MATRDIFF='REDU' ETAB.'MAT_DIF' MODRDIFF;
  1210. 'FINS' ;
  1211. 'FINS' ;
  1212.  
  1213. *-------------- Cas des autres variables --------------------
  1214. 'REPETER' BOU_NU2 ILVA;
  1215. ITE2 = &BOU_NU2;
  1216. MOTI = 'EXTR' LVAREX ITE2;
  1217. 'SI' ('EXIS' LPRIM2 MOTI) ;
  1218. 'ITER' BOU_NU2 ;
  1219. 'FINS' ;
  1220. 'SI' ETAB.'COUPLAGE';
  1221. 'SI' ('EXIS' CHAR_ME MOTI );
  1222. 'SI' ('NEG' LAMBDA1 0.);
  1223. 'MESS' '*** Le materiau thermique depend de parametres mecaniques';
  1224. 'ERRE' '*** Il faut initialiser SOUS-RELAXATION a 0.';
  1225. 'FINS';
  1226. 'ITERER' BOU_NU2;
  1227. 'FINS';
  1228. 'FINS';
  1229. 'SI' ('EXIS' ETAB.'CHARGEMENT' MOTI);
  1230. CH_1 = 'TIRE' ETAB.'CHARGEMENT' MOTI TEMM;
  1231. TYP_1 = 'TYPE' CH_1;
  1232.  
  1233. * on suppose que les chargements qui ne sont pas d'origine mecanique
  1234. * sont appuyes sur le maillage thermique.
  1235. 'SI' ('EGA' TYP_1 'CHPOINT ');
  1236. CHVAR2 = 'CHAN' 'CHAM' CH_1 MODRTOT 'RIGIDITE';
  1237. 'SINO';
  1238. CHVAR2 = 'CHAN' 'RIGIDITE' CH_1 MODRTOT;
  1239. 'FINS';
  1240. CHVAR = CHVAR + CHVAR2;
  1241. 'SINO';
  1242. 'MESS' '*** Il manque une donnee du chargement ***';
  1243. 'MESS' '*** Variable concernee : 'MOTI;
  1244. 'ERREUR' ' ';
  1245. 'FINS';
  1246. 'FIN' BOU_NU2;
  1247. 'FINS';
  1248.  
  1249. *----------------- Calcul du champ de materiau a T= U_BOU1--------------
  1250. 'SI' NONLIN1 ;
  1251. CHVAR1 = 'REDU' CHVAR MODRTOT;
  1252. CHMAT = 'VARI' 'NUAG' MODRTOT MATRTOT0 CHVAR1 'RIGIDITE';
  1253.  
  1254. 'SI' (IDIFF 'OU' ETAB.'THE1') ;
  1255. CHMAT1 = 'REDU' CHMAT ETAB.'MOD_THE' ;
  1256. 'SI' ITHM1 ;
  1257. ma ff0 = @MATETHM MODRTHM TH_COUR ;
  1258. CHMAT1 = CHMAT1 'ET' ma ;
  1259. 'FINS' ;
  1260. 'SINO' ;
  1261. 'SI' ITHM1 ;
  1262. CHMAT1 ff0 = @MATETHM MODRTHM TH_COUR ;
  1263. 'FINS' ;
  1264. 'FINS' ;
  1265.  
  1266. 'SINO';
  1267. CHMAT = MATRTOT0 ;
  1268. 'SI' (IDIFF 'OU' ETAB.'THE1') ;
  1269. CHMAT1 = 'REDU' MATRTOT ETAB.'MOD_THE' ;
  1270. 'SI' ITHM1 ;
  1271. ma ff0 = @MATETHM MODRTHM TH_COUR ;
  1272. CHMAT1=CHMAT1 'ET' ma ;
  1273. 'FINS' ;
  1274. 'SINO' ;
  1275. 'SI' ITHM1 ;
  1276. CHMAT1 ff0 = @MATETHM MODRTHM TH_COUR ;
  1277. 'FINS' ;
  1278. 'FINS' ;
  1279. 'FINS';
  1280.  
  1281. ETAB.'CHMAT' ='REDU' CHMAT MODTOT ;
  1282.  
  1283. MATRTOT ='REDU' CHMAT1 MODRTOT ;
  1284. 'SI' ICONDU;
  1285. MATRCOND='REDU' CHMAT1 MODRCOND;
  1286. 'FINS';
  1287. 'SI' ICONVE;
  1288. MATRCON ='REDU' CHMAT1 MODRCON ;
  1289. 'FINS';
  1290. 'SI' IRADI;
  1291. MATRRAY ='REDU' CHMAT1 MODRRAY ;
  1292. 'FINS';
  1293. 'SI' IADVEC;
  1294. MATRADV ='REDU' CHMAT1 MODRADV ;
  1295. 'FINS' ;
  1296. 'SI' IPHASE;
  1297. MATRPHA ='REDU' CHMAT1 MODRPHA ;
  1298. MATPHA ='REDU' CHMAT1 MODPHA ;
  1299. 'FINS';
  1300. 'SI' ITHM1;
  1301. MATRTHM ='REDU' CHMAT1 MODRTHM ;
  1302. 'FINS' ;
  1303. 'SI' IDIFF ;
  1304. MATRDIFF='REDU' CHMAT1 MODRDIFF;
  1305. 'FINS' ;
  1306.  
  1307. *-------------- CALCUL matrice de conductivite si besoin --------------
  1308. 'SI' IKDEP ;
  1309. 'SI' ICONDU; BB = ('COND' MODRCOND MATRCOND) 'ET' BB;'FINS';
  1310. 'SI' ITHM1 ; BB = ('COND' MODRTHM MATRTHM ) 'ET' BB;'FINS';
  1311. 'FINS' ;
  1312. 'SI' IECDEP;
  1313. 'SI' ICONVE; BB = ('COND' MODRCON MATRCON ) 'ET' BB;'FINS';
  1314. 'FINS';
  1315. 'SI' IKVITE; BB = ('ADVE' MODRADV MATRADV ) 'ET' BB;'FINS';
  1316. 'SI' IKDDEP ; BB = ('COND' MODRDIFF MATRDIFF) 'ET' BB;'FINS';
  1317.  
  1318. *-------------- CALCUL matrice de capacite si besoin------------------
  1319. XDIV = -1.d0 * T_PAS ;
  1320. 'SI' ICDEP ;
  1321. CDETX = 'VIDE' 'RIGIDITE'/'RIGIDITE' ;
  1322. 'SI' ICONDU ;
  1323. CDETY = 'CAPA' MODRCOND MATRCOND ;
  1324. CC = (CDETY/XDIV) 'ET' CC ;
  1325. CDETX = CDETX et CDETY ;
  1326. 'FINS';
  1327. 'SI' ITHM1 ;
  1328. CDETY = 'CAPA' MODRTHM MATRTHM ;
  1329. CC = (CDETY/XDIV) 'ET' CC ;
  1330. CDETX = CDETX et CDETY ;
  1331. 'FINS';
  1332. ETAB.'CTHE_COURANTE' = CDETX ;
  1333. 'FINS';
  1334. 'SI' ICDDEP ;
  1335. CDETY = 'CAPA' MODRDIFF MATRDIFF ;
  1336. CC = (CDETY/XDIV) 'ET' CC ;
  1337. ETAB.'CDIF_COURANTE' = CDETY ;
  1338. 'FINS';
  1339.  
  1340. *-------------- CALCUL des termes de convection si besoin ---------
  1341. 'SI' ICONVE;
  1342. 'SI' IECDEP;
  1343. 'SI' (('EXIS' ETAB.'CHARGEMENT' 'TECO') 'OU' ('EXIS' MATRCON 'TC')) ;
  1344. CHAL_TE = 'CONV' MODRCON MATRCON CH_TEC ;
  1345. MAT_CHPO = MAT_CHPO + CHAL_TE ;
  1346. 'FINS';
  1347. 'FINS';
  1348. 'FINS';
  1349.  
  1350. *
  1351. *------------- calcul des termes rayonnement --------------------------
  1352. 'SI' IRADI ;
  1353. IAPPEL = 3 ;
  1354. t_rayo = PAS_RAYO PRECED TEM1 IAPPEL ;
  1355. 'SI' ('EXIS' t_rayo 'ADDI_MATRICE') ;
  1356. bb = bb 'ET' t_rayo.'ADDI_MATRICE' ;
  1357. 'FINS';
  1358. 'SI' ('EXIS' t_rayo 'ADDI_SECOND') ;
  1359. MAT_CHPO = MAT_CHPO 'ET' t_rayo.'ADDI_SECOND' ;
  1360. 'FINS';
  1361. 'SI' ('EXIS' t_rayo 'RAYO_VALEUR');
  1362. STAB.'RAYONNEMENT' = t_rayo.'RAYO_VALEUR';
  1363. 'FINS';
  1364. t_rayo = 'MOT' 'INCONNU' ;
  1365. 'FINS';
  1366.  
  1367. *------------- en cas de procedure utilisateur CHARTHER --------------
  1368.  
  1369. 'SI' etab.'PROCEDURE_CHARTHER' ;
  1370. TAAT = charther preced TEM1 ;
  1371. 'SI' ('EXIS' TAAT 'ADDI_MATRICE') ;
  1372. bb = bb 'ET' taat.'ADDI_MATRICE' ;
  1373. 'FINS';
  1374. 'SI' ('EXIS' taat 'ADDI_SECOND') ;
  1375. MAT_CHPO = MAT_CHPO 'ET' taat.'ADDI_SECOND' ;
  1376. 'FINS';
  1377. oublier taat ;
  1378. 'FINS' ;
  1379.  
  1380. *--------------- calcul K*u1 et (lambda1*k + C/DT) * DDU1
  1381. 'SI' ('EGA' I_LAMBD1 1) ;
  1382. DS = BB * U0 ;
  1383. 'SINO' ;
  1384. 'SI' ('NEG' I_LAMBD1 3) ;
  1385. chp_z = LAMBDA1 * DU1 ;
  1386. 'SINO' ;
  1387. chp_z = DU1 ;
  1388. 'FINS' ;
  1389. DS = BB * ( chp_z + U0 ) ;
  1390. 'FINS' ;
  1391. DS2 = CC * DU1;
  1392. MAT_CHPO = MAT_CHPO + DS2 - DS ;
  1393. 'SI' ITHM1 ;
  1394. ma FTGC = @MATETHM MODRTHM TH_COUR ;
  1395. MAT_CHPO = MAT_CHPO - FTGC ;
  1396. 'FINS';
  1397. 'FIN' BOU_IT1;
  1398.  
  1399. 'SI' IPHASE ;
  1400. * Evolution de la proportion de phase (impossible actuellement en parallele car Non-Local)
  1401. PROPHA1 ='PROP' MODTOT PROPHA0 Q_LAT_INT DU1 ;
  1402. PROPHA0 = PROPHA1 ;
  1403. 'FINS' ;
  1404.  
  1405. TH_COUR = U1 - CHTREF ;
  1406. ETAB.'THER_COURANT' = TH_COUR ;
  1407.  
  1408. *--------------- Appel a la procedure utilisateur PERSO2 ---------------
  1409. 'SI' ETAB.'PROCEDURE_PERSO2';
  1410. * CB215821 : Il serait bien de mettre dans ESTIMATION la solution COURANTE
  1411. * 'SI' ('EXIS' PRECED 'TEMPERATURES' ) ;
  1412. * tab_z = 'INDEX' PRECED.'TEMPERATURES' ;
  1413. * IND0 = tab_z. ('DIME' tab_z) ;
  1414. * ETAB.'TEMI' = PRECED.'TEMPERATURES' . IND0 ;
  1415. * 'SINO';
  1416. * tab_z = 'INDEX' PRECED.'CONCENTRATIONS' ;
  1417. * IND0 = tab_z. ('DIME' tab_z) ;
  1418. * ETAB.'TEMI' = PRECED.'CONCENTRATIONS' . IND0 ;
  1419. * 'FINS' ;
  1420. PERSO2 PRECED ;
  1421. * ETAB = 'ENLEVER' ETAB 'TEMI' ;
  1422. 'SI' ('EXIS' PRECED 'ARRET') ;
  1423. ETAB.'ARRET' = PRECED.'ARRET' ;
  1424. 'FINS' ;
  1425. 'FINS';
  1426.  
  1427. *******************************************************************
  1428. * Reinitialisation des variables en vue du pas suivant *
  1429. *******************************************************************
  1430. 'SI' ((ITE3+1) '&lt;EG' DCALC);
  1431. TEM0 = TEM1;
  1432. TEM1 ='EXTR' TE_CALC (ITE3 + 1);
  1433. T_PAS = TEM1 - TEM0 ;
  1434. U0 = U1 ;
  1435. U0MTREF= U0 - CHTREF ;
  1436.  
  1437. 'SINO';
  1438. 'QUIT' BOUC_SO;
  1439. 'FINS';
  1440.  
  1441. *------------------------ On fait le menage -----------------------
  1442. 'MENA' ;
  1443. 'FIN' BOUC_SO ;
  1444.  
  1445. 'SAUT' 1 'LIGN' ;
  1446.  
  1447. STAB.'ERREUR' = ERR_VAL ;
  1448. 'SI' ETAB.'FOR_THER' ;
  1449. CHTMP ='EXCO' TH_COUR ETAB.'PRIM_T' 'NOID' ;
  1450. STAB.'INITIAL(2)' ='CHAN' 'ATTRIBUT' CHTMP 'NATURE' 'DIFFUS';
  1451. STAB.'CLIT' = MAT_BLTH ;
  1452. 'FINS' ;
  1453. 'SI' IDIFF ;
  1454. CHTMP ='EXCO' TH_COUR ETAB.'PRIM_DIF' 'NOID' ;
  1455. STAB.'CO2' ='CHAN' 'ATTRIBUT' CHTMP 'NATURE' 'DIFFUS';
  1456. STAB.'CLID' = MAT_BLDI ;
  1457. 'FINS' ;
  1458.  
  1459. 'SI' IRIB ;
  1460. STAB.'RIBLO_T' = RIBL_VAL ;
  1461. STAB.'LISEA_T' = LISE_VAL ;
  1462. 'FINS' ;
  1463.  
  1464. 'SI' IPHASE ;
  1465. STAB.'PROPPHAS' = PROPHA1 ;
  1466. 'FINS' ;
  1467.  
  1468. *Ecriture des Proportions de Phases en sortie (Modele Metallurgique par T.L. en mai 2018)
  1469. 'SI' IMETAL ;
  1470. STAB.'PROPORTIONS_PHASES' ='REDU' PHA_MET1 ETAB.'MOD_MET' ;
  1471. 'FINS' ;
  1472.  
  1473. * Sauvegarde capacite courante :
  1474. 'SI' (ICDEP ou ICDDEP) ;
  1475. ETAB.'CAPA_COURANTE' = ETAB.'CTHE_COURANTE' et ETAB.'CDIF_COURANTE' ;
  1476. 'SI' ICPCONS ;
  1477. ETAB.'CAPA_COURANTE' = ETAB.'CAPA_COURANTE' et ETAB.'CAPACITE_CONSTANTE' ;
  1478. 'FINSI' ;
  1479. 'FINSI' ;
  1480.  
  1481. 'OPTI' 'PARA' FAUX ;
  1482. 'FINPROC' STAB ;
  1483.  

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