Télécharger transnon.procedur

Retour à la liste

Numérotation des lignes :

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

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