Télécharger la plaquette Cast3M


fr en

Télécharger pas_init.procedur

Retour à la liste

Numérotation des lignes : 

  1. * PAS_INIT  PROCEDUR  MB234859  24/02/29    21:15:04     11849          
  2. 'DEBPROC' PAS_INIT   PRECED*'TABLE'  ;
  3. ************************************************************************
  4. * Procedure qui renseigne les champs dans la table CONTINUATION.
  5. * et qui creee les tables ou seront sauves les resultats.
  6. * 
  7. * Trois cas sont a distinguer : 
  8. * - en cas d'un premier appel a PASAPAS, CONTINUATION n'existe pas 
  9. *   et elle est creee et completee de facon a contenir les champs 
  10. *   a l'instant TEMPS0 = 0.
  11. * - en cas de reprise d'un calcul, la procedure PAS_DEFA s'est chargee 
  12. *   de verifier que cet instant est bien present dans la liste des temps 
  13. *   calcules. Les champs resultats pour les instants posterieurs a cet 
  14. *   instant de reprise ont ete supprimes des tables de sauvegardes.
  15. *   TEMPS0 = PRECED.'TEMPS'. NSOR et il faut renseigner CONTINUATION de 
  16. *   facon a ce qu'elle contienne les champs a cet instant.
  17. * - en cas de poursuite d'un calcul, on repart du dernier pas calcule.
  18. *   La table CONTINUATION est redefinie dans le cas ou le modele 
  19. *   aurait change (cas de la fabrication additive ou des elements ont pu 
  20. *   etre ajoutes entre deux appels de PASAPAS).
  21. *  
  22. * Entree : Table PRECED transmises a PASAPAS  
  23. * Sortie : Aucune mais la table de PASAPAS contient desormais 
  24. *          - des tables pour stocker les resultats de calculs
  25. *          - la table CONTINUATION qui contient   les champs a TEMPS0
  26. *          - la table ESTIMATION   qui contiendra les champs a TI
  27. ************************************************************************
  28. WTAB    = PRECED.'WTABLE';
  29. TABSAUV = WTAB.'TABSAUV';
  30. inouv   = 'NON' ('EXIS' PRECED 'CONTINUATION') ;
  31. 'SI' inouv ;
  32.   PRECED.'CONTINUATION' = 'TABLE';
  33. 'FINSI';
  34. conti = PRECED.'CONTINUATION';
  35. ************************************************************************
  36. * Initialisations de l'instant de calcul
  37. ************************************************************************
  38. 'SI' ('EXIS' PRECED 'TEMPS') ;
  39.   NSOR = ('DIME' (PRECED.'TEMPS')) - 1;
  40.   'SI' inouv ;
  41.     TEMP0 = PRECED.'TEMPS'. NSOR;
  42.   'SINON';
  43.     TEMP0 = conti.'TEMPS';
  44.   'FINSI';
  45. 'SINON';
  46.   NSOR = 0; 
  47.   'SI' WTAB.'FREQUENTIEL' ;
  48.     TEMP0 = WTAB.'TEMPS0' ;
  49.   'SINON' ;
  50.     TEMP0 = 0.D0 ;
  51.   'FINSI' ;
  52.   PRECED.'TEMPS'    = 'TABLE';
  53.   PRECED.'TEMPS'. 0 = TEMP0 ;
  54. 'FINSI';
  55. WTAB.'TEMPS0' = TEMP0;
  56. WTAB.'NSOR'   = NSOR;
  57. conti.'TEMPS' = WTAB.'TEMPS0';
  58. ************************************************************************
  59. * Initialisation indices 'MODELES' de la table de sortie
  60. ************************************************************************
  61. 'SI' WTAB.'CHAR_MODE' ;
  62.  
  63. * Indices MODELES de la table de sortie :
  64.   'SI' ('NON' ('EXIS' PRECED 'MODELES')) ;
  65.     PRECED.'MODELES'     = 'TABLE' ;
  66.     PRECED.'MODELES' . 0 = WTAB.'MODELE' ;
  67. *    PRECED.'MAILLAGES'     = 'TABLE' ;
  68. *    PRECED.'MAILLAGES' . 0 = 'EXTR' (PRECED.'MODELES' . 0) 'MAIL' ;
  69.   'FINSI' ;
  70.   conti.'MODELES'   = PRECED.'MODELES'. NSOR ;
  71. * conti.'MAILLAGES' = PRECED.'MAILLAGES' . NSOR ;
  72.  
  73. 'FINSI' ;
  74. ************************************************************************
  75. * Initialisations pour la THERMIQUE
  76. ************************************************************************
  77. 'SI' WTAB.'FOR_THER';
  78. *
  79. * - Initialisation de la temperature
  80.   'SI' ('EXISTE' PRECED 'TEMPERATURES');
  81.     TH_COUR = PRECED.'TEMPERATURES'. NSOR;
  82.     'SI' (NSOR '>EG' 2) ;
  83.       WTAB.'TET0' = PRECED.'TEMPERATURES'. (NSOR - 2);
  84.     'FINSI';
  85.   'SINON';
  86.     'SI' ('EXIS' WTAB.'MAT_T' 'TINI') ;
  87.       MOTSTEMP = 'EXTR' WTAB.'MOD_T' 'DEPL' ;
  88.       DDDD = VIDE 'MCHAML' ;
  89.       'REPE' BMO ('DIME' MOTSTEMP) ;
  90.         MOTEMPI = 'EXTR' MOTSTEMP &BMO ;
  91.         DDDD = DDDD ET ('EXCO' WTAB.'MAT_T' 'TINI' MOTEMPI) ;
  92.       'FIN' BMO ;
  93.     'SINO' ;
  94.       DDDD = 'ZERO' WTAB.'MOD_T' 'TEMPERAT' ;
  95.     'FINS' ;
  96.     TH_COUR = 'CHANGER' 'CHPO' WTAB.'MOD_T' DDDD;
  97.     PRECED.'TEMPERATURES'    = 'TABLE';
  98.     PRECED.'TEMPERATURES'. 0 = TH_COUR;
  99.   'FINSI';
  100.   CONTI.'TEMPERATURES' = 'CHAN' 'ATTRIBUT' TH_COUR 'NATURE' 'DIFFUS';
  101. *
  102. * - Initialisation des reactions
  103.   'SI' ('NON' ('EXIS' PRECED 'REACTIONS_THERMIQUES')) ;
  104.     'SI' WTAB.'CHAR_BLOT' ;
  105.       BLOTHER = 'TIRE' WTAB.'BLOCAGES_THERMIQUES' 'BLOT' TEMP0 ;
  106.     'SINON' ;
  107.       BLOTHER = WTAB.'BLOCAGES_THERMIQUES' ;
  108.     'FINSI' ;
  109.     QTH_COUR = 'REAC' BLOTHER CONTI.'TEMPERATURES' ;
  110.     PRECED.'REACTIONS_THERMIQUES'     = 'TABLE'   ;
  111.     PRECED.'REACTIONS_THERMIQUES' . 0 = QTH_COUR ;
  112.   'FINS' ;
  113.   conti.'REACTIONS_THERMIQUES' = PRECED.'REACTIONS_THERMIQUES'. NSOR ;
  114. *
  115. * Indice 'TEMPERATURES_APPORT' de la table de sortie :
  116.   'SI' WTAB.'CHAR_MODE' ;
  117.     'SI' ('NON' ('EXIS' PRECED 'TEMPERATURES_APPORT')) ;
  118.       PRECED.'TEMPERATURES_APPORT'     = TABLE ;
  119.       NBT1     = 'DIME' PRECED.'TEMPERATURES' ;
  120.       'REPE' BTAPP1 NBT1 ;
  121.         IBT1     = &BTAPP1 - 1 ;
  122.         PRECED.'TEMPERATURES_APPORT' . IBT1 = PRECED.'TEMPERATURES' . IBT1 ;
  123.       'FIN' BTAPP1 ;
  124.     'FINSI' ;
  125.     conti.'TEMPERATURES_APPORT' = PRECED.'TEMPERATURES_APPORT'. NSOR ;
  126.   'FINSI' ;
  127. *
  128. 'FINS';
  129. ************************************************************************
  130. * Initialisations pour la MECANIQUE
  131. ************************************************************************
  132. 'SI' WTAB.'FOR_MECA' ;
  133.    MO_TOT = WTAB.'MO_TOT';
  134. *
  135. * - Initialisation du deplacement
  136.   'SI' ('NON' ('EXIS' PRECED 'DEPLACEMENTS'));
  137.     DDDD = 'ZERO' MO_TOT 'DEPLACEM';
  138.     DEP_COUR = 'CHANGER' 'CHPO' MO_TOT DDDD;
  139.     'DETRUIRE' DDDD;
  140.     PRECED.'DEPLACEMENTS'    = 'TABLE';
  141.     PRECED.'DEPLACEMENTS'. 0 = DEP_COUR;
  142.   'FINSI';
  143.   conti.'DEPLACEMENTS' = PRECED.'DEPLACEMENTS'. NSOR;
  144. *
  145. * - Initialisation des contraintes
  146.   'SI' ('NON' ('EXIS' PRECED 'CONTRAINTES'));
  147.     SIG_COUR = 'ZERO' MO_TOT 'CONTRAIN';
  148.     PRECED.'CONTRAINTES'    = 'TABLE';
  149.     PRECED.'CONTRAINTES'. 0 = SIG_COUR;
  150.   'FINSI';
  151.   conti.'CONTRAINTES' = PRECED.'CONTRAINTES'. NSOR;
  152. *
  153. * - Initialisation des vitesses
  154.   'SI' (WTAB.'DYNAMIQUE' 'OU' ('EXISTE' WTAB 'MOD_LIA'));
  155.     'SI' ('NON' ('EXISTE' PRECED 'VITESSES'))  ;
  156.       DDDD = 'ZERO' MO_TOT 'DEPLACEMENTS';
  157.       VIT_COUR = 'CHANGER' 'CHPO'  MO_TOT DDDD;
  158.       'DETRUIRE' DDDD;
  159.       PRECED.'VITESSES'    = 'TABLE';
  160.       PRECED.'VITESSES'. 0 = VIT_COUR;
  161.     'FINSI';
  162.     conti.'VITESSES' = PRECED.'VITESSES'. NSOR;
  163.   'FINSI';
  164. *
  165. * - Initialisation des accelerations
  166.   'SI' WTAB.'DYNAMIQUE';
  167.     'SI' ('NON' ('EXISTE' PRECED 'ACCELERATIONS')) ;
  168.       DDDD = 'ZERO' MO_TOT 'DEPLACEMENTS';
  169.       ACC_COUR = 'CHANGER' 'CHPO'  MO_TOT DDDD;
  170.       'DETRUIRE' DDDD;
  171.       PRECED.'ACCELERATIONS'    = 'TABLE';
  172.       PRECED.'ACCELERATIONS'. 0 = ACC_COUR;
  173.     'FINSI';
  174.     conti.'ACCELERATIONS' = PRECED.'ACCELERATIONS'. NSOR;
  175.   'FINSI';
  176. *
  177. * - Initialisation des reactions
  178.   'SI' ('NON' ('EXIS' PRECED 'REACTIONS'));
  179.     'SI' WTAB.'CHAR_BLOM' ;
  180.       BLOMECA = 'TIRE' WTAB.'BLOCAGES_MECANIQUES' 'BLOM' TEMP0 ;
  181.     'SINON' ;
  182.       BLOMECA = WTAB.'BLOCAGES_MECANIQUES' ;
  183.     'FINSI' ;
  184.     REA_COUR = 'REAC' BLOMECA conti.'DEPLACEMENTS';
  185.     PRECED.'REACTIONS'     = 'TABLE' ;
  186.     PRECED.'REACTIONS' . 0 = REA_COUR;
  187.   'FINSI';
  188.   conti.'REACTIONS' = PRECED.'REACTIONS'. NSOR;
  189. *
  190. * - Initialisation des variables internes
  191.   'SI' WTAB.'IPLAVI';
  192.     'SI' ('NON' ('EXIS' PRECED 'VARIABLES_INTERNES')) ;
  193.       VAR_COUR = 'ZERO' MO_TOT 'VARINTER';
  194.       PRECED.'VARIABLES_INTERNES'    = 'TABLE';
  195.       PRECED.'VARIABLES_INTERNES'. 0 = VAR_COUR;
  196.     'FINSI';
  197.     conti.'VARIABLES_INTERNES' = PRECED.'VARIABLES_INTERNES'. NSOR ;
  198.   'FINSI';
  199. *
  200. * - Initialisation des deformations inelastiques
  201.   'SI' ('EXISTE' PRECED 'DEFORMATIONS_INELASTIQUES');
  202.     DEI_COUR = PRECED.'DEFORMATIONS_INELASTIQUES' . NSOR;
  203.   'SINON';
  204.     'SI' inouv ;
  205.       DEI_COUR = 'ZERO' MO_TOT 'DEFINELA';
  206.       'SI' ('EXIS' TABSAUV  'DEFIN');
  207.         'SI' TABSAUV.'DEFIN';
  208.           PRECED.'DEFORMATIONS_INELASTIQUES'    = 'TABLE';
  209.           PRECED.'DEFORMATIONS_INELASTIQUES'. 0 = DEI_COUR;
  210.         'FINSI';
  211.       'FINSI';
  212.     'SINON' ;
  213.       DEI_COUR = conti.'DEFORMATIONS_INELASTIQUES' ;
  214.     'FINSI' ;
  215.   'FINSI';
  216.   conti.'DEFORMATIONS_INELASTIQUES' = DEI_COUR;
  217. *
  218. * - Initialisation des deformations totales
  219.   'SI' ('EXISTE' PRECED 'DEFORMATIONS');
  220.     DEF_COUR = PRECED.'DEFORMATIONS'. NSOR;
  221.   'SINON';
  222.     'SI' inouv ;
  223.       DEF_COUR = 'ZERO' MO_TOT 'DEFORMATIONS';
  224.       'SI' ('EXIS' TABSAUV  'DEFTO');
  225.         'SI' TABSAUV.'DEFTO';
  226.           PRECED.'DEFORMATIONS'     = 'TABLE';
  227.           PRECED.'DEFORMATIONS' . 0 = DEF_COUR;
  228.         'FINSI';
  229.       'FINSI';
  230.     'SINON' ;
  231.       DEF_COUR = conti.'DEFORMATIONS' ;
  232.     'FINSI' ;
  233.   'FINSI';
  234.   conti.'DEFORMATIONS' = DEF_COUR;
  235. * 
  236. * - Initialisation des deformations apport
  237.   'SI' ('EXISTE' PRECED 'DEFORMATIONS_APPORT');
  238.     DEF_APPO = PRECED.'DEFORMATIONS_APPORT'. NSOR ;
  239.   'SINON';
  240.     'SI' inouv;
  241.       DEF_APPO = 'ZERO' MO_TOT 'DEFORMATIONS';
  242.       'SI' ('EXIS' TABSAUV  'DEFAP');
  243.         'SI' TABSAUV.'DEFAP';
  244.           PRECED.'DEFORMATIONS_APPORT'     = 'TABLE';
  245.           PRECED.'DEFORMATIONS_APPORT' . 0 = DEF_APPO;
  246.         'FINSI';
  247.       'FINSI' ;
  248.     'SINON' ;
  249.       DEF_APPO = conti.'DEFORMATIONS_APPORT' ;
  250.     'FINSI';
  251.   'FINSI';
  252.   conti.'DEFORMATIONS_APPORT' = DEF_APPO;
  253.  
  254. * - Initialisation pour un tracking
  255.   'SI' ('EXIS' WTAB.'MODELE' 'MATE' 'EFEM');
  256.     'SI' ('EXIS' PRECED 'TRACKING');
  257.       'SI' PRECED.'TRACKING' ;
  258.         'SI' ('NON' ('EXIS' PRECED 'PROOTS'));
  259.           PRECED.'PROOTS' = 'TABLE';
  260.           PRECED.'PROOTS'. 0 = 0.;
  261.         'FINSI';
  262. *       'SI' inouv; conti.'PROOTS' = PRECED.'PROOTS'.NSOR;'FINS';
  263.         'SI' ('NON' ('EXIS' PRECED 'EROOTS'));
  264.           PRECED.'EROOTS' = 'TABLE';
  265.           PRECED.'EROOTS'. 0 = 0.;
  266.         'FINSI';
  267. *       'SI' inouv; conti.'PROOTS' = PRECED.'PROOTS'.NSOR;'FINS';
  268.         'SI' ('NON' ('EXIS' PRECED 'NROOTS'));
  269.           PRECED.'NROOTS' = 0;
  270.         'FINSI';
  271. *       'SI' inouv; conti.'EROOTS' = PRECED.'EROOTS'.NSOR;'FINS';
  272.         'SI' ('NON' ('EXIS' PRECED 'TABISO'));
  273.           PRECED.'TABISO' = 'TABLE';
  274.           PRECED.'TABISO'. 0 = 0.;
  275.         'FINSI';
  276. *       'SI' inouv; conti.'TABISO' = PRECED.'TABISO'.NSOR;'FINS';
  277.         'SI' ('NON' ('EXIS' PRECED 'RESTH'));
  278.           PRECED.'RESTH' = 'TABLE';
  279.           PRECED.'RESTH'. 0 = 0.;
  280.         'FINSI';
  281. *       'SI' inouv; conti.'RESTH' = PRECED.'RESTH'.NSOR;'FINS';
  282.         'SI' ('NON' ('EXIS' PRECED 'ISOTEST'));
  283.           PRECED.'ISOTEST' = 'TABLE';
  284.           PRECED.'ISOTEST'. 0 = 0.;
  285.         'FINSI';
  286. *       'SI' inouv; conti.'ISOTEST' = PRECED.'ISOTEST'.NSOR;'FINS';
  287.         'SI' ('NON' ('EXIS' PRECED 'ISOTOT'));
  288.           PRECED.'ISOTOT' = 'TABLE';
  289.           PRECED.'ISOTOT'. 0 = 0.;
  290.         'FINSI';
  291. *       'SI' inouv; conti.'ISOTOT' = PRECED.'ISOTOT'.NSOR;'FINS';
  292.         'SI' ('NON' ('EXIS' PRECED 'MESH_ISO'));
  293.           PRECED.'MESH_ISO' = 'TABLE';
  294.           PRECED.'MESH_ISO'. 0 = 0.;
  295.         'FINSI';
  296. *       'SI' inouv; conti.'MESH_ISO' = PRECED.'MESH_ISO'.NSOR;'FINS';
  297.         'SI' ('NON' ('EXIS' PRECED 'MESH'));
  298.           PRECED.'MESH' = 'TABLE';
  299.           PRECED.'MESH'. 0 = 0.;
  300.         'FINSI';
  301. *       'SI' inouv; conti.'MESH' = PRECED.'MESH'.NSOR;'FINS';
  302.       'FINSI';
  303.     'FINSI';
  304.   'FINSI';
  305. *
  306. * - Initialisation pour le non-local HELM
  307.   'SI' WTAB.'NLOC' ;
  308.     'SI' ('EGA' WTAB.'NON_LOCAL' 'HELM');
  309.       WTAB.'HELMHOLTZ' = 'TABLE' ;
  310.       MOLVIA = 'EXTR' WTAB.'NLOC_MODL' 'NLOC' ;
  311.       NHELM = 'DIME' MOLVIA ;
  312.       WTAB.'HELMHOLTZ' . 'N_VARI_NL'= NHELM ;
  313.       GEOHELM = 'EXTR' WTAB.'NLOC_MODL' 'MAIL' ;
  314.       MAHELM = 'REDU' WTAB.'CARACTERISTIQUES' WTAB.'NLOC_MODL' ;
  315. *     liste des parametres variables       
  316.       DDVA = 'EXTR' MAHELM 'COVA' ;
  317.       WTAB.'HELMHOLTZ' . 'NLOC_COVA'= DDVA ; 
  318.       'REPE' BH NHELM ;
  319.         WTAB.'HELMHOLTZ'. &BH  = 'TABLE' ;
  320.         LEMO  = 'EXTR' MOLVIA &BH ;
  321.         chaih1=chai 'pas_init: Helmholtz pour :' LEMO;
  322. *       mess chaih1;
  323.         WTAB.'HELMHOLTZ'. &BH. 'NOM' = LEMO ;
  324.         'SI' ('EXIS' DDVA LEMO) ;
  325.           WTAB.'HELMHOLTZ'. &BH. 'CTE' = FAUX ;
  326.         'SINON' ;
  327.           WTAB.'HELMHOLTZ'. &BH. 'CTE' = VRAI ;
  328.         'FINSI' ;
  329.         LEMOD = 'MODE' GEOHELM 'THERMIQUE' 'ANISOTROPE' ;
  330.         WTAB.'HELMHOLTZ'. &BH . 'H_MODELE' = LEMOD ;
  331. *       non linearite par defaut des helmholtz (sellier )
  332.         WTAB.'HELMHOLTZ'. &BH . 'LINEAIRE' = FAUX ;    
  333.       'FIN' BH ;
  334.     'FINSI';
  335.   'FINSI';
  336. *
  337. * En cas de chargement de deformations initialisation
  338.   'SI' WTAB.'CHAR_DEFI';
  339.      WTAB.'DEFOR1' = 'TIRE' WTAB.'CHARGEMENT' 'DEFI' TEMP0 ;
  340.   'FINSI';
  341. *
  342. * En presence de THERMIQUE (CHARGEMENT ou FORMULATION) initialisation de TEMPERATURE_REFERENCE
  343.   'SI' WTAB.'ITHER' ;
  344. *
  345. *   En cas de chargement thermique initialisation
  346.     'SI' WTAB.'CHAR_THE';
  347.       WTAB.'TET1' = 'TIRE' WTAB.'CHARGEMENT' 'T' TEMP0;
  348.     'SINON';
  349.       WTAB.'TET1' = CONTI.'TEMPERATURES' ;
  350.     'FINSI';
  351. *
  352. *    On prend la temperature de reference dans le MATERIAU
  353.      'SI' ('EXIS' WTAB.'MA_TOT' 'TREF');
  354.         TREF = 'ZERO' WTAB.'MO_TOT' 'TEMPERAT' ;
  355.         TSZ = 'EXTR' WTAB.'MO_TOT' 'ZONE' ;
  356.         ISZ = 1 ;
  357.         'REPE' BSZ (('DIME' TSZ) / 2) ;
  358.            MOSZI = TSZ . ISZ ;
  359.            MOTSTEMP = 'EXTR' MOSZI 'TEMP' ;
  360.            'REPE' BMO ('DIME' MOTSTEMP) ;
  361.               MOTEMPI = 'EXTR' MOTSTEMP &BMO ;
  362.               TREFI = 'CHAN' 'STRESSES' MOSZI ('EXCO' ('REDU' WTAB.'MA_TOT' MOSZI) 'TREF' MOTEMPI) ;
  363.               TREF = TREF + TREFI ;
  364.            'FIN' BMO ;
  365.            ISZ = ISZ + 2 ;
  366.         'FIN' BSZ ;
  367.         WTAB.'TEMPERATURE_REFERENCE' = 'CHAN' 'TYPE' TREF 'TEMPERATURES' ;
  368.      'SINO';
  369.         ERREUR 77 'AVEC' ('MOT' 'TREF') ('MOT' 'de MATERIAU');
  370.      'FINS';
  371. *
  372. ***  Desormais initialise dans PASAPAS, au meme endroit que ETREF
  373. ***  'SI' ('EGA' WTAB.'ETHER1' 'INCONNU');
  374. ***     On initialise les deformations thermiques a ZERO
  375. ***     WTAB.'ETHER1' = 'ZERO' MO_TOT 'DEFORMATIONS' ;
  376. ***     CB215821 : Pour une reprise propre, il faudrait mettre DEFORMATION_THERMIQUE dans CONTI
  377. ***  'FINS' ;
  378.   'FINS';
  379. 'FINS';
  380. ************************************************************************
  381. * Initialisation des differentes configurations (reference, debut pas)
  382. ************************************************************************
  383. WTAB.'FOR0'    = 'FORM' ;
  384. WTAB.'GE0_DEB' = WTAB.'FOR0';
  385. 'SI' WTAB.'FOR_MECA' ;
  386.   'SI'(WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  387.     WTAB.'GE0_DEB' = 'FORM' conti.'DEPLACEMENTS';
  388.   'FINS';
  389.   'FORM' WTAB.'FOR0' ;
  390. 'FINS';
  391. ************************************************************************
  392. * Initialisations pour la MECANIQUE DES FLUIDES
  393. ************************************************************************
  394. 'SI' WTAB.'NVSTNL'                          ;
  395.   'SI' ('NON' ('EXISTE' PRECED 'VITESSES_FLUIDE'))  ; 
  396.     DDDD = 'ZERO' WTAB.'MOD_NSL' 'DEPLACEMENTS'     ;
  397.     VFL_COUR = 'CHANGER' 'CHPO' WTAB.'MOD_NSL' DDDD ;
  398.     'DETRUIRE' DDDD                                 ;
  399.     PRECED.'VITESSES_FLUIDE'     = 'TABLE'          ;
  400.     PRECED.'VITESSES_FLUIDE'. 0  = VFL_COUR        ;
  401.   'FINSI'                                   ;
  402.   conti.'VITESSES_FLUIDE' = PRECED.'VITESSES_FLUIDE'. NSOR;
  403. 'FINSI'                                     ;
  404. ************************************************************************
  405. * Initialisations pour la METALLURGIE (T.L. en juin 2018)
  406. ************************************************************************
  407. 'SI' WTAB.'FOR_METALLU' ;
  408. *
  409. * - Initialisation des proportions de phase
  410.   'SI' ('NON' ('EXIS' PRECED 'PROPORTIONS_PHASES')) ;
  411. *   Valeurs par defaut : ZERO pour toutes les phases sauf la premiere : 1.D0
  412.     CHMET0  ='ZERO' WTAB.'MOD_MET' 'VARINTER';
  413.     LISCO   ='EXTR' WTAB.'MOD_MET' 'VARI';
  414.     CMP1    ='EXCO' CHMET0 ('EXTR' LISCO 1) ('EXTR' LISCO 1);
  415.     CHMET0  = CHMET0 + (CMP1 + 1.D0);
  416.     PRECED.'PROPORTIONS_PHASES'    = 'TABLE' ;
  417.     PRECED.'PROPORTIONS_PHASES'. 0 = CHMET0;
  418.   'FINS' ;
  419.   conti.'PROPORTIONS_PHASES' = PRECED.'PROPORTIONS_PHASES'. NSOR ;
  420.   CHTPOI0  = ('ZERO' WTAB.'MOD_MET' 'TEMPERATURES') 'NOMC' 'TPOI' ;
  421.   CHTPOI0  = 'CHAN' 'CHPO' WTAB.'MOD_MET' CHTPOI0 ;
  422.   WTAB . 'T_POINT0' = CHTPOI0 ;
  423. *  'SI' ('EXIS' PRECED 'T_POINT') ;
  424. *    TAB_TPOI = 'INDE' (PRECED.'T_POINT')   ;
  425. *    CHTPOI0  = PRECED.'T_POINT'.NSOR       ;
  426. *  'SINO' ;
  427. *    CHTPOI0  = ('ZERO' WTAB.'MOD_MET' 'TEMPERATURES') 'NOMC' 'TPOI' ;
  428. *    CHTPOI0  = 'CHAN' 'CHPO' WTAB.'MOD_MET' CHTPOI0 ;
  429. *    PRECED.'T_POINT'     = 'TABLE' ;
  430. *    PRECED.'T_POINT' . 0 = CHTPOI0 ;
  431. *  'FINS' ;
  432. * conti . 'T_POINT'            = CHTPOI0 ;
  433. 'FINSI' ;
  434. ************************************************************************
  435. * Initialisations pour la DIFFUSION
  436. ************************************************************************
  437. 'SI' WTAB.'FOR_DIFF' ;
  438. *
  439. * - Initialisation des concentrations
  440.   'SI' ('NON' ('EXIS' PRECED 'CONCENTRATIONS')) ;
  441.     CHCO0 = 'ZERO' WTAB. 'MOD_DIF' 'DEPLACEMENTS' ;
  442.     CHCO0 = 'CHAN' 'CHPO' WTAB. 'MOD_DIF' CHCO0   ;
  443.     CHCO0 = 'CHAN' 'ATTRIBUT' CHCO0  'NATURE' 'DIFFUS';
  444.     PRECED.'CONCENTRATIONS'    = 'TABLE' ;
  445.     PRECED.'CONCENTRATIONS'. 0 = CHCO0  ;
  446.   'FINS' ;
  447.   CONTI.'CONCENTRATIONS' = PRECED.'CONCENTRATIONS'. NSOR ;
  448. *
  449. * - Initialisation des reactions
  450.   'SI' ('NON' ('EXIS' PRECED 'REACTIONS_DIFFUSIONS')) ;
  451.     'SI' WTAB.'CHAR_BLOD' ;
  452.       BLODIFF = 'TIRE' WTAB.'BLOCAGES_DIFFUSIONS' 'BLOD' TEMP0 ;
  453.     'SINON' ;
  454.       BLODIFF = WTAB.'BLOCAGES_DIFFUSIONS' ;
  455.     'FINSI' ;
  456.     QDI_COUR = 'REAC' BLODIFF CONTI.'CONCENTRATIONS' ;
  457.     PRECED.'REACTIONS_DIFFUSIONS'     = 'TABLE' ;
  458.     PRECED.'REACTIONS_DIFFUSIONS' . 0 = QDI_COUR ;
  459.   'FINS' ;
  460.   CONTI.'REACTIONS_DIFFUSIONS' = PRECED.'REACTIONS_DIFFUSIONS'. NSOR ;
  461. 'FINS' ;
  462. ************************************************************************
  463. * Initialisations pour le CHANGEMENT DE PHASE (CB@!%*@! 2019)
  464. *-----------Doit etre place en dernier car peut necessiter les autres
  465. *-----------champs
  466. ************************************************************************
  467. 'SI' WTAB.'PHASE';
  468. *
  469.   MO_COVA ='EXTR' WTAB.'MAT_PHA' 'COVA';
  470.   'SI' (('DIME' MO_COVA) > 0);
  471.     CHAI1 ='CHAI' ' Le materiau du modele CHANGEMENT_PHASE ne peut pas varier';
  472.     'ERRE' CHAI1 ;
  473.   'FINS';
  474.  
  475.   'SI' ( 'EXISTE' PRECED 'PROPORTIONS_PHASE')  ;
  476.     CHPHA0 ='REDU' WTAB.'MOD_PHA' PRECED.'PROPORTIONS_PHASE'.NSOR ;
  477.   'SINON';
  478. *   On met un etat qui a du sens physique pour la proportion de phases
  479.     TABPHA ='EXTR' WTAB.'MOD_PHA' 'ZONE';
  480.     CHPHA0 ='VIDE' 'MCHAML';
  481.     'REPE' SURi (('DIME' TABPHA) / 2);
  482. *     Boucle sur les modeles CHANGEMENT_PHASE
  483.       MODEi  = TABPHA.((&SURi-1)*2 + 1) ;
  484.       MATEi  =('REDU' MODEi WTAB.'MAT_PHA' ) 'EXCO' 'PRIM';
  485.       LMOPRIM='EXTR' MODEi 'DEPL';
  486.       MOPRIM ='EXTR' LMOPRIM 1;
  487.       'SI' ('EXIS' WTAB.'CHARGEMENT' MOPRIM);
  488. *       On regarde dans le CHARGEMENT si la PRIMALE existe
  489.         CHPRIM ='TIRE' WTAB.'CHARGEMENT' MOPRIM TEMP0;
  490.         'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  491.           CHPRIM='CHAN'  'CHAM'  CHPRIM MODEi 'NOEUD';
  492.           CHPRIM= CHPRIM 'NOMC' 'PRIM';
  493.           CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  494.           CHPHA0= CHPHAi 'ET' CHPHA0;
  495.           'ITER' SURi ;
  496.         'FINS';
  497.       'FINS';
  498.  
  499. *      On regarde dans les indices connus de la TABLE PRECED
  500.       'SI' ('EXIS' PRECED 'TEMPERATURES');
  501.         CHPRIM = TH_COUR;
  502.         'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  503.           CHPRIM='CHAN'  'CHAM'  CHPRIM MODEi 'NOEUD';
  504.           CHPRIM= CHPRIM 'NOMC' 'PRIM';
  505.           CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  506.           CHPHA0= CHPHAi 'ET' CHPHA0;
  507.           'ITER' SURi ;
  508.         'FINS';
  509.       'FINS';
  510.  
  511.       'SI' ('EXIS' PRECED 'CONCENTRATIONS');
  512.         CHPRIM = CHCO0;
  513.         'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  514.           CHPRIM='CHAN'  'CHAM'  CHPRIM MODEi 'NOEUD';
  515.           CHPRIM= CHPRIM 'NOMC' 'PRIM';
  516.           CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  517.           CHPHA0= CHPHAi 'ET' CHPHA0;
  518.           'ITER' SURi ;
  519.         'FINS';
  520.       'FINS';
  521.  
  522. *      Si on arrive ici c'est qu'on n'a pas trouve le champ primal : on met zero
  523.       CHPHAi='ZERO' MODEi 'VARINTER';
  524.       CHPHA0= CHPHAi 'ET' CHPHA0;
  525.     'FIN' SURi;
  526.  
  527.     PRECED.'PROPORTIONS_PHASE'    = 'TABLE' ;
  528.     PRECED.'PROPORTIONS_PHASE'. 0 = CHPHA0;
  529.   'FINSI';
  530.   conti.'PROPORTIONS_PHASE' = CHPHA0;
  531.   WTAB.'BLO_PHASE'='PMAT' WTAB.'MOD_PHA' ;
  532. 'FINSI';
  533. *
  534. 'SI' WTAB.'SOLUBILITE';
  535.   WTAB.'BLO_SOL' = 'PMAT' WTAB.'MOD_SOL' ;
  536. 'FINSI';
  537. *
  538. PRECED.'ESTIMATION' = 'COPIER' conti;
  539. 'FINP';
  540.  
  541.  

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