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/11/04    21:15:04     12064          
  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. NSOR = WTAB.'NSOR' ;
  39. 'SI' ('EXIS' PRECED 'TEMPS') ;
  40.   TEMP0 = PRECED.'TEMPS'. NSOR;
  41. 'SINON';
  42.   'SI' WTAB.'FREQUENTIEL' ;
  43.     'REPETER' BTFRE ;
  44.       TI = 'EXTRAIRE' WTAB.'TEMPS_CALCULES' &BTFRE ;
  45.       'SI' (TI '>' 0.D0) ;
  46.          TEMP0 = TI ;
  47.          'QUITTER' BTFRE ;
  48.       'FINSI' ;
  49.     'FIN' BTFRE ;
  50.   'SINON' ;
  51.     TEMP0 = 0.D0 ;
  52.   'FINSI' ;
  53.   PRECED.'TEMPS'    = 'TABLE';
  54.   PRECED.'TEMPS'. 0 = TEMP0 ;
  55. 'FINSI';
  56. conti.'TEMPS' = TEMP0 ;
  57. ************************************************************************
  58. * Initialisation indices 'MODELES' de la table de sortie
  59. ************************************************************************
  60. 'SI' WTAB.'CHAR_MODE' ;
  61.  
  62. * Indices MODELES de la table de sortie :
  63.   'SI' ('NON' ('EXIS' PRECED 'MODELES')) ;
  64.     PRECED.'MODELES'     = 'TABLE' ;
  65.     PRECED.'MODELES' . 0 = WTAB.'MODELE' ;
  66. *    PRECED.'MAILLAGES'     = 'TABLE' ;
  67. *    PRECED.'MAILLAGES' . 0 = 'EXTR' (PRECED.'MODELES' . 0) 'MAIL' ;
  68.   'FINSI' ;
  69.   conti.'MODELES'   = PRECED.'MODELES'. NSOR ;
  70. * conti.'MAILLAGES' = PRECED.'MAILLAGES' . NSOR ;
  71.  
  72. 'FINSI' ;
  73. ************************************************************************
  74. * Initialisations pour la THERMIQUE
  75. ************************************************************************
  76. 'SI' WTAB.'FOR_THER';
  77. *
  78. * - Initialisation de la temperature
  79.   'SI' ('EXISTE' PRECED 'TEMPERATURES');
  80.     TH_COUR = PRECED.'TEMPERATURES'. NSOR;
  81.     'SI' (NSOR '>EG' 2) ;
  82.       WTAB.'TET0' = PRECED.'TEMPERATURES'. (NSOR - 2);
  83.     'FINSI';
  84.   'SINON';
  85.     'SI' ('EXIS' WTAB.'MAT_T' 'TINI') ;
  86.       MOTSTEMP = 'EXTR' WTAB.'MOD_T' 'DEPL' ;
  87.       DDDD = VIDE 'MCHAML' ;
  88.       'REPE' BMO ('DIME' MOTSTEMP) ;
  89.         MOTEMPI = 'EXTR' MOTSTEMP &BMO ;
  90.         DDDD = DDDD ET ('EXCO' WTAB.'MAT_T' 'TINI' MOTEMPI) ;
  91.       'FIN' BMO ;
  92.     'SINO' ;
  93.       DDDD = 'ZERO' WTAB.'MOD_T' 'TEMPERAT' ;
  94.     'FINS' ;
  95.     TH_COUR = 'CHANGER' 'CHPO' WTAB.'MOD_T' DDDD;
  96.     PRECED.'TEMPERATURES'    = 'TABLE';
  97.     PRECED.'TEMPERATURES'. 0 = TH_COUR;
  98.   'FINSI';
  99.   conti.'TEMPERATURES' = 'CHAN' 'ATTRIBUT' TH_COUR 'NATURE' 'DIFFUS';
  100. *
  101. * - Initialisation des reactions
  102.   'SI' ('NON' ('EXIS' PRECED 'REACTIONS_THERMIQUES')) ;
  103.     'SI' WTAB.'CHAR_BLOT' ;
  104.       BLOTHER = 'TIRE' WTAB.'BLOCAGES_THERMIQUES' 'BLOT' TEMP0 ;
  105.     'SINON' ;
  106.       BLOTHER = WTAB.'BLOCAGES_THERMIQUES' ;
  107.     'FINSI' ;
  108.     QTH_COUR = 'REAC' BLOTHER conti.'TEMPERATURES' ;
  109.     PRECED.'REACTIONS_THERMIQUES'     = 'TABLE'   ;
  110.     PRECED.'REACTIONS_THERMIQUES' . 0 = QTH_COUR ;
  111.   'FINS' ;
  112.   conti.'REACTIONS_THERMIQUES' = PRECED.'REACTIONS_THERMIQUES'. NSOR ;
  113. *
  114. * Indice 'TEMPERATURES_APPORT' de la table de sortie :
  115.   'SI' WTAB.'CHAR_MODE' ;
  116.     'SI' ('NON' ('EXIS' PRECED 'TEMPERATURES_APPORT')) ;
  117.       PRECED.'TEMPERATURES_APPORT'     = TABLE ;
  118.       NBT1     = 'DIME' PRECED.'TEMPERATURES' ;
  119.       'REPE' BTAPP1 NBT1 ;
  120.         IBT1     = &BTAPP1 - 1 ;
  121.         PRECED.'TEMPERATURES_APPORT' . IBT1 = PRECED.'TEMPERATURES' . IBT1 ;
  122.       'FIN' BTAPP1 ;
  123.     'FINSI' ;
  124.     conti.'TEMPERATURES_APPORT' = PRECED.'TEMPERATURES_APPORT'. NSOR ;
  125.   'FINSI' ;
  126. *
  127. 'FINS';
  128. ************************************************************************
  129. * Initialisations pour la MECANIQUE
  130. ************************************************************************
  131. 'SI' WTAB.'FOR_MECA' ;
  132.    MO_TOT = WTAB.'MO_TOT';
  133. *
  134. * - Initialisation du deplacement
  135.   'SI' ('NON' ('EXIS' PRECED 'DEPLACEMENTS'));
  136.     DDDD = 'ZERO' MO_TOT 'DEPLACEM';
  137.     DEP_COUR = 'CHANGER' 'CHPO' MO_TOT DDDD;
  138.     'DETRUIRE' DDDD;
  139.     PRECED.'DEPLACEMENTS'    = 'TABLE';
  140.     PRECED.'DEPLACEMENTS'. 0 = DEP_COUR;
  141.   'FINSI';
  142.   conti.'DEPLACEMENTS' = PRECED.'DEPLACEMENTS'. NSOR;
  143. *
  144. * - Initialisation des contraintes
  145.   'SI' ('NON' ('EXIS' PRECED 'CONTRAINTES'));
  146.     SIG_COUR = 'ZERO' MO_TOT 'CONTRAIN';
  147.     PRECED.'CONTRAINTES'    = 'TABLE';
  148.     PRECED.'CONTRAINTES'. 0 = SIG_COUR;
  149.   'SINON';
  150.     SIG_COUR = PRECED.'CONTRAINTES'. NSOR ;
  151.     'SI' WTAB.'GRANDS_DEPLACEMENTS' ;
  152.       SIG_COUR = 'CAPI' WTAB.'HYPOTHESE_DEFORMATIONS' SIG_COUR MO_TOT conti.'DEPLACEMENTS' ;
  153.     'FINSI';
  154.   'FINSI';
  155.   conti.'CONTRAINTES' = SIG_COUR ;
  156. *
  157. * - Initialisation des vitesses
  158.   'SI' (WTAB.'DYNAMIQUE' 'OU' ('EXISTE' WTAB 'MOD_LIA'));
  159.     'SI' ('NON' ('EXISTE' PRECED 'VITESSES'))  ;
  160.       DDDD = 'ZERO' MO_TOT 'DEPLACEMENTS';
  161.       VIT_COUR = 'CHANGER' 'CHPO'  MO_TOT DDDD;
  162.       'DETRUIRE' DDDD;
  163.       PRECED.'VITESSES'    = 'TABLE';
  164.       PRECED.'VITESSES'. 0 = VIT_COUR;
  165.     'FINSI';
  166.     conti.'VITESSES' = PRECED.'VITESSES'. NSOR;
  167.   'FINSI';
  168. *
  169. * - Initialisation des accelerations
  170.   'SI' WTAB.'DYNAMIQUE';
  171.     'SI' ('NON' ('EXISTE' PRECED 'ACCELERATIONS')) ;
  172.       DDDD = 'ZERO' MO_TOT 'DEPLACEMENTS';
  173.       ACC_COUR = 'CHANGER' 'CHPO'  MO_TOT DDDD;
  174.       'DETRUIRE' DDDD;
  175.       PRECED.'ACCELERATIONS'    = 'TABLE';
  176.       PRECED.'ACCELERATIONS'. 0 = ACC_COUR;
  177.     'FINSI';
  178.     conti.'ACCELERATIONS' = PRECED.'ACCELERATIONS'. NSOR;
  179.   'FINSI';
  180. *
  181. * - Initialisation des reactions
  182.   'SI' ('NON' ('EXIS' PRECED 'REACTIONS'));
  183.     'SI' WTAB.'CHAR_BLOM' ;
  184.       BLOMECA = 'TIRE' WTAB.'BLOCAGES_MECANIQUES' 'BLOM' TEMP0 ;
  185.     'SINON' ;
  186.       BLOMECA = WTAB.'BLOCAGES_MECANIQUES' ;
  187.     'FINSI' ;
  188.     REA_COUR = 'REAC' BLOMECA conti.'DEPLACEMENTS';
  189.     PRECED.'REACTIONS'     = 'TABLE' ;
  190.     PRECED.'REACTIONS' . 0 = REA_COUR;
  191.   'FINSI';
  192.   conti.'REACTIONS' = PRECED.'REACTIONS'. NSOR;
  193. *
  194.   'SI' WTAB.'IPLAVI';
  195. *
  196. *   - Initialisation des variables internes
  197.     'SI' ('NON' ('EXIS' PRECED 'VARIABLES_INTERNES')) ;
  198.       VAR_COUR = 'ZERO' MO_TOT 'VARINTER';
  199.       PRECED.'VARIABLES_INTERNES'    = 'TABLE';
  200.       PRECED.'VARIABLES_INTERNES'. 0 = VAR_COUR;
  201.     'FINSI';
  202.     conti.'VARIABLES_INTERNES' = PRECED.'VARIABLES_INTERNES'. NSOR ;
  203. *
  204. *   - Initialisation des deformations inelastiques
  205.     'SI' ('EXISTE' PRECED 'DEFORMATIONS_INELASTIQUES');
  206.       DEI_COUR = PRECED.'DEFORMATIONS_INELASTIQUES' . NSOR;
  207.     'SINON';
  208.       'SI' inouv ;
  209.         DEI_COUR = 'ZERO' MO_TOT 'DEFINELA';
  210.         'SI' ('EXIS' TABSAUV  'DEFIN');
  211.           'SI' TABSAUV.'DEFIN';
  212.             PRECED.'DEFORMATIONS_INELASTIQUES'    = 'TABLE';
  213.             PRECED.'DEFORMATIONS_INELASTIQUES'. 0 = DEI_COUR;
  214.           'FINSI';
  215.         'FINSI';
  216.       'SINON' ;
  217.         DEI_COUR = conti.'DEFORMATIONS_INELASTIQUES' ;
  218.       'FINSI' ;
  219.     'FINSI';
  220.     conti.'DEFORMATIONS_INELASTIQUES' = DEI_COUR;
  221. *
  222.   'FINSI';
  223. *
  224. * - Initialisation des deformations totales
  225.   'SI' ('EXISTE' PRECED 'DEFORMATIONS');
  226.     DEF_COUR = PRECED.'DEFORMATIONS'. NSOR ;
  227.     'SI' WTAB.'GRANDS_DEPLACEMENTS' ;
  228.       DEF_COUR = 'CAPI' WTAB.'HYPOTHESE_DEFORMATIONS' DEF_COUR MO_TOT conti.'DEPLACEMENTS' ;
  229.     'FINSI';
  230.   'SINON';
  231.     'SI' inouv ;
  232.       DEF_COUR = 'ZERO' MO_TOT 'DEFORMATIONS';
  233.       'SI' ('EXIS' TABSAUV  'DEFTO');
  234.         'SI' TABSAUV.'DEFTO';
  235.           PRECED.'DEFORMATIONS'     = 'TABLE';
  236.           PRECED.'DEFORMATIONS' . 0 = DEF_COUR;
  237.         'FINSI';
  238.       'FINSI';
  239.     'SINON' ;
  240.       DEF_COUR = conti.'DEFORMATIONS' ;
  241.     'FINSI' ;
  242.   'FINSI';
  243.   conti.'DEFORMATIONS' = DEF_COUR;
  244. * 
  245. * - Initialisation des deformations apport
  246.   'SI' ('EXISTE' PRECED 'DEFORMATIONS_APPORT');
  247.     DEF_APPO = PRECED.'DEFORMATIONS_APPORT'. NSOR ;
  248.   'SINON';
  249.     'SI' inouv;
  250.       DEF_APPO = 'ZERO' MO_TOT 'DEFORMATIONS';
  251.       'SI' ('EXIS' TABSAUV  'DEFAP');
  252.         'SI' TABSAUV.'DEFAP';
  253.           PRECED.'DEFORMATIONS_APPORT'     = 'TABLE';
  254.           PRECED.'DEFORMATIONS_APPORT' . 0 = DEF_APPO;
  255.         'FINSI';
  256.       'FINSI' ;
  257.     'SINON' ;
  258.       DEF_APPO = conti.'DEFORMATIONS_APPORT' ;
  259.     'FINSI';
  260.   'FINSI';
  261.   conti.'DEFORMATIONS_APPORT' = DEF_APPO;
  262.  
  263. * - Initialisation pour un tracking
  264.   'SI' ('EXIS' WTAB.'MODELE' 'MATE' 'EFEM');
  265.     'SI' ('EXIS' PRECED 'TRACKING');
  266.       'SI' PRECED.'TRACKING' ;
  267.         'SI' ('NON' ('EXIS' PRECED 'PROOTS'));
  268.           PRECED.'PROOTS' = 'TABLE';
  269.           PRECED.'PROOTS'. 0 = 0.;
  270.         'FINSI';
  271.         'SI' ('NON' ('EXIS' PRECED 'EROOTS'));
  272.           PRECED.'EROOTS' = 'TABLE';
  273.           PRECED.'EROOTS'. 0 = 0.;
  274.         'FINSI';
  275.         'SI' ('NON' ('EXIS' PRECED 'NROOTS'));
  276.           PRECED.'NROOTS' = 0;
  277.         'FINSI';
  278.         'SI' ('NON' ('EXIS' PRECED 'TABISO'));
  279.           PRECED.'TABISO' = 'TABLE';
  280.           PRECED.'TABISO'. 0 = 0.;
  281.         'FINSI';
  282.         'SI' ('NON' ('EXIS' PRECED 'RESTH'));
  283.           PRECED.'RESTH' = 'TABLE';
  284.           PRECED.'RESTH'. 0 = 0.;
  285.         'FINSI';
  286.         'SI' ('NON' ('EXIS' PRECED 'ISOTEST'));
  287.           PRECED.'ISOTEST' = 'TABLE';
  288.           PRECED.'ISOTEST'. 0 = 0.;
  289.         'FINSI';
  290.         'SI' ('NON' ('EXIS' PRECED 'ISOTOT'));
  291.           PRECED.'ISOTOT' = 'TABLE';
  292.           PRECED.'ISOTOT'. 0 = 0.;
  293.         'FINSI';
  294.         'SI' ('NON' ('EXIS' PRECED 'MESH_ISO'));
  295.           PRECED.'MESH_ISO' = 'TABLE';
  296.           PRECED.'MESH_ISO'. 0 = 0.;
  297.         'FINSI';
  298.         'SI' ('NON' ('EXIS' PRECED 'MESH'));
  299.           PRECED.'MESH' = 'TABLE';
  300.           PRECED.'MESH'. 0 = 0.;
  301.         'FINSI';
  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.   conti.'VITESSES_FLUIDE_0' = PRECED.'VITESSES_FLUIDE'. NSOR;
  404. 'FINSI'                                     ;
  405. ************************************************************************
  406. * Initialisations pour la METALLURGIE (T.L. en juin 2018)
  407. ************************************************************************
  408. 'SI' WTAB.'FOR_METALLU' ;
  409. *
  410. * - Initialisation des proportions de phase
  411.   'SI' ('NON' ('EXIS' PRECED 'PROPORTIONS_PHASES')) ;
  412. *   Valeurs par defaut : ZERO pour toutes les phases sauf la premiere : 1.D0
  413.     CHMET0  ='ZERO' WTAB.'MOD_MET' 'VARINTER';
  414.     LISCO   ='EXTR' WTAB.'MOD_MET' 'VARI';
  415.     CMP1    ='EXCO' CHMET0 ('EXTR' LISCO 1) ('EXTR' LISCO 1);
  416.     CHMET0  = CHMET0 + (CMP1 + 1.D0);
  417.     PRECED.'PROPORTIONS_PHASES'    = 'TABLE' ;
  418.     PRECED.'PROPORTIONS_PHASES'. 0 = CHMET0;
  419.   'FINS' ;
  420.   conti.'PROPORTIONS_PHASES' = PRECED.'PROPORTIONS_PHASES'. NSOR ;
  421.   CHTPOI0  = ('ZERO' WTAB.'MOD_MET' 'TEMPERATURES') 'NOMC' 'TPOI' ;
  422.   CHTPOI0  = 'CHAN' 'CHPO' WTAB.'MOD_MET' CHTPOI0 ;
  423.   WTAB . 'T_POINT0' = CHTPOI0 ;
  424. *  'SI' ('EXIS' PRECED 'T_POINT') ;
  425. *    TAB_TPOI = 'INDE' (PRECED.'T_POINT')   ;
  426. *    CHTPOI0  = PRECED.'T_POINT'.NSOR       ;
  427. *  'SINO' ;
  428. *    CHTPOI0  = ('ZERO' WTAB.'MOD_MET' 'TEMPERATURES') 'NOMC' 'TPOI' ;
  429. *    CHTPOI0  = 'CHAN' 'CHPO' WTAB.'MOD_MET' CHTPOI0 ;
  430. *    PRECED.'T_POINT'     = 'TABLE' ;
  431. *    PRECED.'T_POINT' . 0 = CHTPOI0 ;
  432. *  'FINS' ;
  433. * conti . 'T_POINT'            = CHTPOI0 ;
  434. 'FINSI' ;
  435. ************************************************************************
  436. * Initialisations pour la DIFFUSION
  437. ************************************************************************
  438. 'SI' WTAB.'FOR_DIFF' ;
  439. *
  440. * - Initialisation des concentrations
  441.   'SI' ('NON' ('EXIS' PRECED 'CONCENTRATIONS')) ;
  442.     CHCO0 = 'ZERO' WTAB. 'MOD_DIF' 'DEPLACEMENTS' ;
  443.     CHCO0 = 'CHAN' 'CHPO' WTAB. 'MOD_DIF' CHCO0   ;
  444.     CHCO0 = 'CHAN' 'ATTRIBUT' CHCO0  'NATURE' 'DIFFUS';
  445.     PRECED.'CONCENTRATIONS'    = 'TABLE' ;
  446.     PRECED.'CONCENTRATIONS'. 0 = CHCO0  ;
  447.   'FINS' ;
  448.   conti.'CONCENTRATIONS' = PRECED.'CONCENTRATIONS'. NSOR ;
  449. *
  450. * - Initialisation des reactions
  451.   'SI' ('NON' ('EXIS' PRECED 'REACTIONS_DIFFUSIONS')) ;
  452.     'SI' WTAB.'CHAR_BLOD' ;
  453.       BLODIFF = 'TIRE' WTAB.'BLOCAGES_DIFFUSIONS' 'BLOD' TEMP0 ;
  454.     'SINON' ;
  455.       BLODIFF = WTAB.'BLOCAGES_DIFFUSIONS' ;
  456.     'FINSI' ;
  457.     QDI_COUR = 'REAC' BLODIFF conti.'CONCENTRATIONS' ;
  458.     PRECED.'REACTIONS_DIFFUSIONS'     = 'TABLE' ;
  459.     PRECED.'REACTIONS_DIFFUSIONS' . 0 = QDI_COUR ;
  460.   'FINS' ;
  461.   conti.'REACTIONS_DIFFUSIONS' = PRECED.'REACTIONS_DIFFUSIONS'. NSOR ;
  462. 'FINS' ;
  463. ************************************************************************
  464. * Initialisations pour le CHANGEMENT DE PHASE (CB@!%*@! 2019)
  465. *-----------Doit etre place en dernier car peut necessiter les autres
  466. *-----------champs
  467. ************************************************************************
  468. 'SI' WTAB.'PHASE';
  469. *
  470.   MO_COVA ='EXTR' WTAB.'MAT_PHA' 'COVA';
  471.   'SI' (('DIME' MO_COVA) > 0);
  472.     CHAI1 ='CHAI' ' Le materiau du modele CHANGEMENT_PHASE ne peut pas varier';
  473.     'ERRE' CHAI1 ;
  474.   'FINS';
  475.  
  476.   'SI' ( 'EXISTE' PRECED 'PROPORTIONS_PHASE')  ;
  477.     CHPHA0 ='REDU' WTAB.'MOD_PHA' PRECED.'PROPORTIONS_PHASE'.NSOR ;
  478.   'SINON';
  479. *   On met un etat qui a du sens physique pour la proportion de phases
  480.     TABPHA ='EXTR' WTAB.'MOD_PHA' 'ZONE';
  481.     CHPHA0 ='VIDE' 'MCHAML';
  482.     'REPE' SURi (('DIME' TABPHA) / 2);
  483. *     Boucle sur les modeles CHANGEMENT_PHASE
  484.       MODEi  = TABPHA.((&SURi-1)*2 + 1) ;
  485.       MATEi  =('REDU' MODEi WTAB.'MAT_PHA' ) 'EXCO' 'PRIM';
  486.       LMOPRIM='EXTR' MODEi 'DEPL';
  487.       MOPRIM ='EXTR' LMOPRIM 1;
  488.       'SI' ('EXIS' WTAB.'CHARGEMENT' MOPRIM);
  489. *       On regarde dans le CHARGEMENT si la PRIMALE existe
  490.         CHPRIM ='TIRE' WTAB.'CHARGEMENT' MOPRIM TEMP0;
  491.         'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  492.           CHPRIM='CHAN'  'CHAM'  CHPRIM MODEi 'NOEUD';
  493.           CHPRIM= CHPRIM 'NOMC' 'PRIM';
  494.           CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  495.           CHPHA0= CHPHAi 'ET' CHPHA0;
  496.           'ITER' SURi ;
  497.         'FINS';
  498.       'FINS';
  499.  
  500. *      On regarde dans les indices connus de la TABLE PRECED
  501.       'SI' ('EXIS' PRECED 'TEMPERATURES');
  502.         CHPRIM = TH_COUR;
  503.         'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  504.           CHPRIM='CHAN'  'CHAM'  CHPRIM MODEi 'NOEUD';
  505.           CHPRIM= CHPRIM 'NOMC' 'PRIM';
  506.           CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  507.           CHPHA0= CHPHAi 'ET' CHPHA0;
  508.           'ITER' SURi ;
  509.         'FINS';
  510.       'FINS';
  511.  
  512.       'SI' ('EXIS' PRECED 'CONCENTRATIONS');
  513.         CHPRIM = CHCO0;
  514.         'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  515.           CHPRIM='CHAN'  'CHAM'  CHPRIM MODEi 'NOEUD';
  516.           CHPRIM= CHPRIM 'NOMC' 'PRIM';
  517.           CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  518.           CHPHA0= CHPHAi 'ET' CHPHA0;
  519.           'ITER' SURi ;
  520.         'FINS';
  521.       'FINS';
  522.  
  523. *      Si on arrive ici c'est qu'on n'a pas trouve le champ primal : on met zero
  524.       CHPHAi='ZERO' MODEi 'VARINTER';
  525.       CHPHA0= CHPHAi 'ET' CHPHA0;
  526.     'FIN' SURi;
  527.  
  528.     PRECED.'PROPORTIONS_PHASE'    = 'TABLE' ;
  529.     PRECED.'PROPORTIONS_PHASE'. 0 = CHPHA0;
  530.   'FINSI';
  531.   conti.'PROPORTIONS_PHASE' = CHPHA0;
  532.   WTAB.'BLO_PHASE' = 'PMAT' WTAB.'MOD_PHA' ;
  533. 'FINSI';
  534. *
  535. 'SI' WTAB.'SOLUBILITE';
  536.   WTAB.'BLO_SOL' = 'PMAT' WTAB.'MOD_SOL' ;
  537. 'FINSI';
  538. *
  539. PRECED.'ESTIMATION' = 'COPIER' conti;
  540. 'FINP';
  541.  
  542.  

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