Télécharger la plaquette Cast3M


fr en

Télécharger transnon.procedur

Retour à la liste

Numérotation des lignes : 

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

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