Télécharger la plaquette Cast3M


fr en

Télécharger pasapas.procedur

Retour à la liste

Numérotation des lignes : 

  1. * PASAPAS   PROCEDUR  MB234859  23/01/09    21:15:05     11548          
  2.  
  3. 'DEBP' PASAPAS  PRECED*'TABLE';
  4. 'SAUT' 1 'LIGNE';
  6. '------------------ DEBUT DE LA PROCEDURE PASAPAS --------------------';
  7. ************************************************************************
  8. *                                                                      *
  9. *                                                                      *
  10. *** PRECED      objet de type TABLE                                    *
  11. *               C'est une table donnant , les arguments d'entrees, les *
  12. *               options, puis les  resultats des pas precedents        *
  13. *                                                                      *
  14. * Etat initial au debut du pas de temps decrit par le MCHELM WTAB.ETAT1*
  15. * champs de materiaux :                                                *
  16. *     WTAB.'MAT_REF' materiau etat a la temperature de REFERENCE       *
  17. *     WTAB.'MAT0'    materiau etat debut de calcul                     *
  18. *     WTAB.'MAT1'    materiau fin du pas de temps                      *
  19. *                                                                      *
  20. ************************************************************************
  21.  
  22. ****************   I PARTIE   ******************************************
  23.  
  24. *---------- attribution des valeurs par defaut et verifs elementaires
  25. PAS_DEFA PRECED ;
  26. *--------  initialisation du pas precedent si besoin et creation des tables de sorties.
  27. PAS_INIT PRECED ;
  28.  
  29. WTAB  = PRECED.'WTABLE'      ;
  30. conti = PRECED.'CONTINUATION';
  31. estim = PRECED.'ESTIMATION'  ;
  32. *
  33. WTAB.'GE0_DEB'=WTAB.'FOR0';
  34. 'SI'( WTAB.'GRANDS_DEPLACEMENTS' 'OU'  WTAB.'FEFP_FORMULATION');
  35.   'SI' ( 'EXIS' conti 'DEPLACEMENTS');
  36.     WTAB.'GE0_DEB'= 'FORME' conti . 'DEPLACEMENTS';
  37.   'FINSI';
  38. 'FINSI';
  39. 'FORM' WTAB.'FOR0';
  40.  
  41. *--------- quel est le delta-temps du dernier pas calcule ? ---------
  42. 'SI' ('EGA' WTAB.'DT' ('MOT' 'INCONNU'));
  43.   DELTAN = 0.D0 ;
  44.   WTAB.'DT' = DELTAN ;
  45. 'SINON';
  46.   DELTAN= WTAB.'DT';
  47. 'FINSI';
  48.  
  49. *------------Appel a une procedure utilisateur PERSO1
  50. 'SI' WTAB.'PROCEDURE_PERSO1';
  51.   PRECED.'PERSO1_APPEL' = 1 ;
  52.   PERSO1 PRECED ;
  53.   'SI' ('EXIS' PRECED 'ARRET') ;
  54.     WTAB.'ARRET' = PRECED.'ARRET' ;
  55.   'FINSI' ;
  56. 'FINSI';
  57.  
  58. * ------------Appel a une procedure utilisateur REEV_MEC
  59. 'SI' WTAB.'PROCEDURE_REEV_MEC' ;
  60.   REEV_MEC PRECED 0 ;
  61. 'FINSI' ;
  62.  
  63. * ------------Appel a une procedure utilisateur REEV_THE
  64. 'SI' WTAB.'PROCEDURE_REEV_THE' ;
  65.   REEV_THE PRECED 0 ;
  66. 'FINSI' ;
  67.  
  68. *---------   initialisation ETAT ET MATERIAUX
  69. WTAB.'ETAT0'   = PAS_ETAT PRECED WTAB.'TEMPS0'            ;
  70. WTAB.'ETAT1'   = WTAB.'ETAT0'                             ;
  71. WTAB.'MAT0'    = PAS_MATE PRECED WTAB.'ETAT0'             ;
  72. WTAB.'MAT1'    = WTAB.'MAT0'                              ;
  73.  
  74. *---------   initialisation deformation thermique de reference
  75. 'SI' WTAB.'FOR_MECA' ;
  76.   'SI' (WTAB.'CHAR_THE' 'OU'  WTAB.'FOR_THER' ) ;
  77.     'SI' ('EGA' WTAB.'ETREF' 'INCONNU');
  78.       WTAB.'ETAT_REF'= PAS_ETAT PRECED WTAB.'TEMPS0' 'REFERENCE';
  79.       WTAB.'MAT_REF' = PAS_MATE PRECED WTAB.'ETAT_REF'          ;
  80.       WTAB.'ETREF'   ='EPTH' WTAB.'MO_TOT' WTAB.'MAT_REF' WTAB.'TEMPERATURE_REFERENCE' ;
  81.     'FINS' ;
  82. *
  83.     'SI' ('EGA' ('TYPE' WTAB.'TET1') 'CHPOINT');
  84.       TET1 ='CHAN' 'CHAM' WTAB.'TET1' WTAB.'MO_TOT' 'STRESSES' 'TEMPERATURES' ;
  85.     'FINS';
  86.     ETT0  = 'EPTH' WTAB.'MO_TOT' WTAB.'MAT0' TET1 ;
  87.     WTAB.'ETHER1' = ETT0 - WTAB.'ETREF' ;
  88. *
  89.     'SI' WTAB.'POR1';
  90.       WTAB.'MSRTHER1'= PAS_EPTH PRECED WTAB.'MO_TOT' WTAB.'MAT0' TET1 ;
  91.     'FINS' ;
  92.   'FINS' ;
  93. 'FINS' ;
  94.  
  95. **************************   II PARTIE  *************************
  96. ******************************************************************
  97. ******************* boucle sur les pas de temps  *****************
  98. ******************************************************************
  99.  
  100. * ------------quelques initialisations de l'algorithme
  101. PASFINAL = 0 ;
  102. IPAS     = 0 ;
  103. ICHG     = 0 ;
  104. IAUTO    = 1 ;
  105. IREREDU  = VRAI;
  106. TEMP0    = WTAB.'TEMPS0';
  107. ISAUV    = VRAI;
  108. DTV      =('EXTR' WTAB.'TEMPS_CALCULES' 1) - TEMP0 ;
  109. 'SI' ( DTV '<' 1.e-15 );
  110.   'SI' WTAB.'FREQUENTIEL' ;
  111.     DTV = 1.e-14 ;
  112.   'SINON' ;
  113.     DTV =('EXTR' WTAB.'TEMPS_CALCULES' 2) - TEMP0 ;
  114.   'FINSI';
  115. 'FINSI';
  116. WTAB.'ARRET'   = FAUX;
  117. WTAB.'AUTODEUX'= FAUX;
  118. WTAB.'CONV'    = VRAI;
  119. DT_AVANT       = DTV ;
  120. WTAB.'NOMBRE_ITERATIONS'=12;
  121.  
  122. * Attention Bizarre : RFROT est initialise ici et utilise dans UNPAS ????
  123. RFROT=0;
  124.  
  125. *---------- debut boucle sur pas de temps
  126. 'REPETER' BEXTERN  ;
  127. *
  128. * ----- Determiner le nouvel instant a calculer TI
  129. *
  130.   NPAS='DIME' WTAB.'TEMPS_CALCULES';
  131.   ICHG=ICHG + 1;  WTAB.'ICHG'= ICHG;
  132.   'SI' ( ICHG '>' NPAS) ; 'QUITTER' BEXTERN; 'FINSI';
  133.   TI = 'EXTRAIRE'  WTAB.'TEMPS_CALCULES' ICHG ;
  134.   TIV=TI;
  135.   'SI' (WTAB.'FREQUENTIEL') ;
  136.     'SI' (TI '<' TEMP0 ) ;  'ITERER'  BEXTERN ; 'FINSI' ;
  137.   'SINON' ;
  138.     ITERTPS = (TI '<' TEMP0 ) 'OU' (TI 'EGA' PRECED.'TEMPS'. 0) ;
  139.     'SI' ('NON' WTAB.'REEQUILIBRAGE') ;
  140.       ITERTPS = ITERTPS 'OU' ('EGA' TI TEMP0) ;
  141.     'FINSI';
  142.     'SI' ITERTPS ;
  143.       'ITERER' BEXTERN ;
  144.     'FINSI';
  145.   'FINSI';
  146. *
  147. * ----- Ajustement du pas en fonction du nb d iteration du pas precedent
  148. *
  149.   PASFINAL=PASFINAL+1;
  150.   'SI' WTAB.'PAS_AJUSTE';
  151.     DTV=DT_AVANT;
  152.     'SI' WTAB.'CONV';
  153.       'SI' (WTAB.'NOMBRE_ITERATIONS' > 20);
  154.         DTV=DTV'/' (2. ** 0.5);
  155.       'FINSI';
  156.       'SI' (WTAB.'NOMBRE_ITERATIONS' <  7);
  157.         'SI' IREREDU;
  158.           DTV=DTV '*' (2. ** 0.5);
  159.         'FINSI';
  160.       'FINSI';
  161.       DT_AVANT= DTV;
  162.       IREREDU=VRAI;
  163.     'SINON';
  164.       DTV=0.0000000001D0* DT_AVANT;
  165.       DT_AVANT= DT_AVANT /2.;
  166.       IREREDU=FAUX;
  167.     'FINSI';
  168.     TTI = DTV* 1.0000000001 + TEMP0;
  169.     'SI' ( TTI '&lt;EG' TIV );
  170.       TI=TTI;
  171.       PASFINAL=0;
  172.     'FINSI';
  173.   'FINSI';
  174. *
  175.   estim.'TEMPS'=TI;
  176.   'SI' (PASFINAL '>EG' 1.5) ; DT_AVANT= TI-TEMP0; 'FINSI';
  177.   WTAB.'T_FINAL'= TI;
  178.   WTAB.'DT_INIT' = TI-TEMP0;
  179.   'SI'((WTAB.'CONSOLIDATION' 'OU' WTAB.'DYNAMIQUE')
  180.        'ET' (WTAB.'DT_INIT' 'EGA' 0.));
  181.     'ITERER' BEXTERN;
  182.   'FINSI';
  183.   WTAB.'DT'= WTAB.'DT_INIT';
  184. *
  185. * ----- Compteur IPAS des pas de temps calcules
  186.   IPAS = IPAS + 1;
  187.   WTAB.'PAS'=IPAS;
  188.  
  189. * ----- Mise a jour du modele si chargement de type MODE :
  190.   'SI' (WTAB.'CHAR_MODE') ;
  191.     MOD1P     =  WTAB.'MODELE' ;
  192.     MODEVOL1  = 'TIRE' PRECED.'MODELE' 'MODE' TI ;
  193.     WTAB.'MODELE' = MODEVOL1 ;
  194.  
  195.     'SI' (WTAB.'CHAR_MATE') ;
  196.       MATEVOL1  = 'TIRE' PRECED.'CARACTERISTIQUES' 'MATE' TI ;
  197.       WTAB.'CARACTERISTIQUES' = MATEVOL1 ;
  198.     'FINSI' ;
  199.  
  200.     PAS_MODL PRECED ;
  201.  
  202. *   MaJ table rayonnement :
  203.     'SI' WTAB.'RAYO' ;
  204.       WTAB.'RAYONNEMENT' = 'PRRA' WTAB.'MOD_RAY' ;
  205.     'FINSI' ;
  206.  
  207. *   On sauvegarde les modeles successifs pour analyse resultats :
  208.     ESTIM.'MODELES'   = WTAB.'MODELE' ;
  209.  
  210. *   Initialisation TINI si thermique :
  211.     'SI' (WTAB.'THE1') ;
  212.       MODX      = WTAB.'MOD_COND' ;
  213.       MODXP     = MOD1P 'EXTR' 'MATE' 'CONDUCTION' ;
  214.       GEO1      = MODX  'EXTR' 'MAIL' ;
  215.       GEO1P     = MODXP 'EXTR' 'MAIL' ;
  216.       DGEO1     = GEO1  'DIFF' GEO1P ;
  217.  
  218. *     Adapation du modele : transport de la solution
  219.       'SI' ((WTAB.'ADAP_MODE') 'ET' ('NON' (TI '<' (WTAB.'TPS1_ADAP')))) ;
  220.         mess '*****  Transport solution adaptation' ;
  221.         IADAP1    = (WTAB.'INDI_ADAP') + 1 ;
  222.         'SI' (IADAP1 '&lt;EG' ('DIME' (WTAB.'LTPS_ADAP'))) ;
  223.           WTAB.'TPS1_ADAP' = WTAB.'LTPS_ADAP' 'EXTR' IADAP1 ;
  224.           WTAB.'INDI_ADAP' = IADAP1 ;
  225.         'SINON' ;
  226.           WTAB.'TPS1_ADAP' = 'MAXI' WTAB.'TEMPS_CALCULES' + 1. ;
  227.         'FINSI' ;
  228. *erre 21 ;
  229. *       Projection du champ de temperature sur la partie adaptee :
  230.         DGEO1P    = DGEO1 'INTE' GEO1P ;
  231.         DGEO2P    = (DGEO1 'INTE' GEO1) 'INCLU' DGEO1P 'VOLU' 'BARY'  ;
  232.         CHT1P     = (WTAB.'TET1') 'EXCO' 'T' 'T' ;
  233.         CHAMT1P   = 'CHAN' 'CHAM' DGEO1P CHT1P ;
  234.         CHT2P     = 'PROI' CHAMT1P  DGEO2P ;
  235.         CHT1P     = ('REDU' CHT1P (GEO1P 'DIFF' DGEO1P)) 'ET' CHT2P ;
  236.         WTAB.'TET1' = CHT1P ;
  237.  
  238. *       Calcul de la capacite courante sur la partie adaptee :
  239.         DMODX     = 'REDU' MODX DGEO2P ;
  240.         DMATX     = 'REDU' (WTAB.'CARACTERISTIQUES') DMODX ;
  241.         DMATX     = 'VARI' 'NUAG' DMODX DMATX CHT2P ;
  242.         CAPAX     = 'CAPA' DMODX DMATX ;
  243.         CAPA1P    = ('REDU' (WTAB.'CAPA_COURANTE') (GEO1P 'DIFF' DGEO1P)) 'ET' CAPAX ;
  244.         WTAB.'CAPA_COURANTE' = CAPA1P ;
  245.  
  246.         DGEO1     = (DGEO1 'INTE' GEO1) 'DIFF' DGEO2P ;
  247.       'FINSI' ;
  248.  
  249. *     Calcul temperature dans partie ajoutee du modele :
  250.       'SI' ('NON' (('VIDE' DGEO1) 'OU' ('VIDE' (DGEO1 'INTE' GEO1)))) ;
  251. *       Quantite de chaleur sur geometrie GEO1P :
  252.         CHT1P     = (WTAB.'TET1') 'EXCO' 'T' 'T' ;
  253.         CHQ1P     = (WTAB.'CAPA_COURANTE' lump) * CHT1P ;
  254. *       Quantite de chaleur dans ajout :
  255.         DMODX     = 'REDU' MODX DGEO1 ;
  256.         DMATX     = 'REDU' (WTAB.'CARACTERISTIQUES') DMODX ;
  257.         CHTINIX   = DMATX 'EXCO' 'TINI' 'T' ;
  258.         CHTINIX   = 'CHAN' 'CHPO' CHTINIX DMODX 'MOYE' ;
  259.         DMATX     = 'VARI' 'NUAG' DMODX DMATX CHTINIX ;
  260.         CAPAX     = 'CAPA' DMODX DMATX ;
  261.         CHQX      = (CAPAX lump) * CHTINIX ;
  262. *       Calcul T sur GEO1 :
  263.         CAPA1     = (WTAB.'CAPA_COURANTE') et CAPAX ;
  264.         CHTX1     = 'MANU' 'CHPO' GEO1 1 'T' 1. 'NATURE' 'DIFFUS' ;
  265.         CHQX1     = CAPA1 * CHTX1;
  266.         CHT1      = (CHQ1P + CHQX) * ('INVE' CHQX1) ('MOTS' 'Q') ('MOTS' 'Q') ('MOTS' 'T') ;
  267. *       CHT2      = reso (lump CAPA1) (CHQ1P + CHQX) ; mess '############ avec RESO' ;
  268. *si (ipas > 10) ; erre 21 ; fins ;
  269.         WTAB.'TET1' = CHT1 ;
  270. *List ((resu (CHQ1P + CHQX)) maxi abs) ;
  271. *List ((resu (CAPA1*CHT1)) maxi abs) ;
  272. *trac CHT1 GEO1 ;
  273.       'FINSI' ;
  274.     'FINSI' ;
  275.  
  276. *   Initialisation PROPORTIONS_PHASES si metallurgie :
  277.     'SI' WTAB.'FOR_METALLU' ;
  278.       MODX      = (WTAB.'MODELE') 'EXTR' 'FORM' 'METALLURGIE' ;
  279.       MODXP     = MOD1P           'EXTR' 'FORM' 'METALLURGIE' ;
  280.       GEO1      = MODX  'EXTR' 'MAIL' ;
  281.       GEO1P     = MODXP 'EXTR' 'MAIL' ;
  282.       DGEO1     = GEO1 'DIFF' GEO1P ;
  283.       'SI' ('NON' (('VIDE' DGEO1) 'OU' ('VIDE' (DGEO1 'INTE' GEO1)))) ;
  284.         CHMET0    = ESTIM.'PROPORTIONS_PHASES' ;
  285.         DMODX     = 'REDU' MODX DGEO1 ;
  286.         CHMET1    = 'ZERO' DMODX 'VARINTER' ;
  287.         LISCO     = 'EXTR' DMODX 'VARI' ;
  288.         CMP1      = 'EXCO' CHMET1 ('EXTR' LISCO 1) ('EXTR' LISCO 1) ;
  289.         CHMET1    = CHMET1 + (CMP1 + 1.D0);
  290.         CHMET1    = CHMET0 'ET' CHMET1 ;
  291.         ESTIM.'PROPORTIONS_PHASES' = CHMET1 'REDU' MODX ;
  292. *trac CHMET1 MODX ;
  293.       'FINSI' ;
  294.     'FINS' ;
  295.  
  296.   'FINSI' ;
  297.  
  298. * ----- Mise a jour blocages passes dans chargement
  299. * Cas BLOCAGES_MECANIQUES (BLOM) :
  300.   'SI' (WTAB.'CHAR_BLOM') ;
  301.     BLOM1    = 'TIRE' PRECED.'BLOCAGES_MECANIQUES' 'BLOM' TI ;
  302.     WTAB.'BLOCAGES_MECANIQUES' = BLOM1 ;
  303.  
  304. *   MAJ Operateur Resolution :
  305.     'SI' ('EXIS' WTAB 'RRRR') ;
  306.       RRRR = WTAB.'RRRR' 'EXTR' 'RIGI' 'NOMU' ;
  307.       WTAB.'RRRR' = RRRR 'ET' BLOM1 ;
  308.     'FINSI' ;
  309.   'FINSI' ;
  310.  
  311. * Cas BLOCAGES_THERMIQUES (BLOT) :
  312.   'SI' (WTAB.'CHAR_BLOT') ;
  313.     BLOT1    = 'TIRE' PRECED.'BLOCAGES_THERMIQUES' 'BLOT' TI ;
  314.     WTAB.'BLOCAGES_THERMIQUES' = BLOT1 ;
  315.   'FINSI' ;
  316.  
  317. * Cas BLOCAGES_DIFFUSIONS (BLOD) :
  318.   'SI' (WTAB.'CHAR_BLOD') ;
  319.     BLOD1    = 'TIRE' PRECED.'BLOCAGES_DIFFUSIONS' 'BLOD' TI ;
  320.     WTAB.'BLOCAGES_DIFFUSIONS' = BLOD1 ;
  321.   'FINSI' ;
  322.  
  323. * ----- On se place dans la configuration deformee
  324.   'SI'( WTAB.'GRANDS_DEPLACEMENTS' 'OU'  WTAB.'FEFP_FORMULATION');
  325.     'SI' (WTAB.'RECALCUL');
  326.       'SI' ('EGA'  ( 'TYPE' WTAB.'GE0_DEB')  'CONFIGUR');
  327.         'FORM'  WTAB.'GE0_DEB';
  328.       'FINSI';
  329.     'FINSI';
  330.   'FINSI';
  331.  
  332. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  333. *++++++++++++++++++ BOUCLE interne de bouclage mec_ther ++++++++++++++
  334. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  335.   I_BOTH = 0 ;
  336.   'REPETER' BO_BOTH WTAB.'BOU_METH' ;
  337.     I_BOTH = I_BOTH + 1 ;
  338. *
  339. *  re-initialisation de WTAB.'CHARGEMENT' en tete de boucle mec-ther
  340. *
  341.     'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  342.       WTAB.'CHARGEMENT' = WTAB . 'CHARGEMENT_I' ;
  343.     'FINSI';
  344. *
  345. *#######################################################################
  346. *#################### CALCUL D'UN PAS SI THERMIQUE #####################
  347. *#######################################################################
  348.     'SI' (WTAB.'FOR_THER' 'OU' WTAB.'FOR_DIFF') ;
  349.       'SI' WTAB.'FOR_MECA' ;
  350. *---------- en cas de convergence mecanique thermique
  351.         'SI' (I_BOTH > 1) ;
  352. *----------   reprise de la configuration meca de la derniere estimation
  353. *             du pas present si coherence mecanique/thermique
  354.           'SI' ('EXIS' estim 'FOR' );
  355.             WTAB .'RECARI'  = VRAI ;
  356.             WTAB .'RECADET' = VRAI ;
  357.             WTAB.'REA_GEOM' = VRAI ;
  358.             'FORM' estim.'FOR'      ;
  359.           'FINSI';
  360.         'FINSI' ;
  361.       'FINSI' ;
  362.  
  363.       'SI' (TI 'NEG' TEMP0);
  364.         'SI' ('NON' ('EXISTE' PRECED 'PAS_MAX'));
  365.           WTAB.'TEM_CALC' = 'PROG' TI;
  366.         'SINON';
  367.           'SI' ((TEMP0 + PRECED.'PAS_MAX') '>EG' TI);
  368.             WTAB.'TEM_CALC' = 'PROG' TI;
  369.           'SINON';
  370.             WTAB.'TEM_CALC' = 'PROG' TEMP0
  371.                                    'PAS'  PRECED.'PAS_MAX'  TI;
  372.             WTAB.'TEM_CALC' = 'ENLE'  WTAB.'TEM_CALC' 1;
  373.           'FINSI';
  374.         'FINSI';
  375. *
  376. * ----- Appel a la procedure thermique (lineaire ou nonlineaire) ----
  377. *
  378.         'SI' ('NON' WTAB.'COUPLAGE'); 'MENAGE'; 'FINSI';
  379.         'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'NONLINEAIRE');
  380.           CHTER = TRANSNON PRECED;
  381.         'FINSI';
  382.         'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'LINEAIRE');
  383.           CHTER = TRANSLIN PRECED;
  384.         'FINSI';
  385.         'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'DUPONT');
  386.           WTAB.'NOUV_PAS' = ('EXTR' WTAB.'TEM_CALC' 1) - TEMP0;
  387.           CHTER = DUPONT2 PRECED;
  388.         'FINSI';
  389.  
  390. *       Verification d'une erreur en sortie
  391.         'SI' CHTER.'ERREUR';
  392.           'ERREUR' ' ';
  393.         'FINSI';
  394.  
  395. *       Modele Metallurgique (T.L. en juin 2018)
  396.         'SI' WTAB.'FOR_METALLU' ;
  397.           estim.'PROPORTIONS_PHASES' = CHTER.'PROPORTIONS_PHASES' ;
  398.         'FINSI';
  399.  
  400.         'SI' WTAB.'PHASE';
  401.           WTAB.'PROPPHAS'          = CHTER.'PROPPHAS' ;
  402.           estim.'PROPORTIONS_PHASE'= CHTER.'PROPPHAS' ;
  403.         'FINSI' ;
  404.  
  405.         'SI' WTAB.'FOR_THER' ;
  406.           WTAB.'TET2'=CHTER.'INITIAL(2)' ;
  407.           WTAB.'CLIT'=CHTER.'CLIT' ;
  408.           'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'NONLINEAIRE') ;
  409.             'SI' ('NEG' WTAB.'RIBLO_T' 'INCONNU');
  410.               WTAB.'RIBLO_T'       = CHTER.'RIBLO_T'       ;
  411.               WTAB.'LISEA_T'       = CHTER.'LISEA_T'       ;
  412.             'FINSI' ;
  413.           'FINSI' ;
  414.           'SI' WTAB.'RAYO' ;
  415.             WTAB.'RAYONNEMENT' = CHTER.'RAYONNEMENT' ;
  416.           'FINSI' ;
  417.         'FINSI' ;
  418.         'SI' WTAB.'FOR_DIFF' ;
  419.           WTAB.'CO2' =CHTER.'CO2';
  420.           WTAB.'CLID'=CHTER.'CLID' ;
  421.         'FINSI' ;
  422.  
  423.         'SI' WTAB.'FOR_THER' ;
  424.           estim.'TEMPERATURES' = WTAB.'TET2' ;
  425.           estim.'REACTIONS_THERMIQUES'='REAC' WTAB.'CLIT' WTAB.'TET2';
  426.         'FINSI';
  427.  
  428.         'SI' WTAB.'FOR_DIFF' ;
  429.           estim.'CONCENTRATIONS' = WTAB.'CO2' ;
  430.           estim.'REACTIONS_DIFFUSIONS' ='REAC' WTAB.'CLID' WTAB.'CO2';
  431.         'FINSI' ;
  432. *** Reevaluation des donnees du calcul dependant de la temperature
  433.         'SI' WTAB.'PROCEDURE_REEV_THE' ;
  434.           REEV_THE PRECED  1 ;
  435.         'FINSI' ;
  436.       'FINS';
  437.  
  438. *-----convergence mecanique thermique : reprise configuration  ------
  439. *---------- du pas precedent pour refaire la mecanique  -------------
  440.       'SI'( WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  441.         'FORM' WTAB.'GE0_DEB';
  442.       'FINSI' ;
  443.     'FINSI' ;
  444. *################## fin du calcul d'un pas thermique #################
  445.  
  446. *#######################################################################
  447. *#################### CALCUL D'UN PAS SI MECANIQUE #####################
  448. *#######################################################################
  449.     'SI' WTAB.'FOR_MECA' ;
  450.       'SAUT' 1 'LIGN';
  451.       'MESSAGE' ' Numero du pas : ' IPAS '  Indice d evolution : ' ICHG ' -> temps :' TI;
  452.  
  453. *--------------------- calcul de la temperature ---------------------
  454. * CB215821 : On met la thermique dans le Chargement pour 'TIRER' le chargement
  455. *            au bon instant dans UNPAS pour la MECANIQUE
  456. *            Il faudra faire de meme avec la MECANIQUE si la THERMIQUE en depend
  457. *        'SI' WTAB.'FOR_THER';
  458. *           THTA1='TABLE';
  459. *           THTA2='TABLE';
  460. *           THTA1.0=TEMP0;
  461. *           THTA1.1=TI;
  462. *           THTA2.0= conti.'TEMPERATURES' 'ENLE' 'LX' ;
  463. *           THTA2.1= estim.'TEMPERATURES' 'ENLE' 'LX' ;
  464. *           CHTH='CHAR' 'T' THTA1 THTA2;
  465. *           CHASANST= WTAB.'CHARGEMENT' 'ENLE' 'T';
  466. *           WTAB.'CHARGEMENT'= CHASANST 'ET' CHTH;
  467. *        'FINSI';
  468.       'SI' WTAB.'CHAR_THE';
  469.         WTAB.'TET2' ='TIRE' WTAB.'CHARGEMENT' 'T   ' TI;
  470.       'FINSI';
  471.  
  472. *------- Caracteristiques initiales en cas de grands deplacements -----
  473.       'SI' WTAB.'ITCAR';
  474.         'SI'  ('EGA' WTAB.'CARA'  'INCONNU');
  475.           WTAB.'CARA' = WTAB.'MAT0';
  476.         'FINSI';
  477.       'FINSI';
  478.  
  479. ************************************************************************
  480. ****************    APPEL A UNPAS       ********************************
  481. ************************************************************************
  482.       'MENA' ;
  483.       TT = UNPAS  PRECED;
  484.  
  485. *---------- Mise a jour du MODELE du PAS precedent ---------------
  486.       WTAB.'MO_TOT_PREC' = TT.'MO_TOT_PREC' ;
  487.  
  488. *---------- retour a la configuration initiale -------------------
  489.       'SI' (WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  490.         'FORM'  WTAB.'FOR0';
  491.       'FINSI';
  492.  
  493. *------- test si on a converge
  494.       WTAB.'CONV' = TT.'CONV';
  495.       'SI' ('NON' TT.'CONV' 'ET' ('NON' WTAB.'CONVERGENCE_FORCEE'));
  496. *          en dynamique et en poreux on garde le PRECED.FREA1 du dernier pas
  497. *          converge (dans le cas ou on voudrait continuer le calcul)
  498. *          'SI' WTAB.'DYNAMIQUE' ; WTAB.'FREA1' = F1 ; 'FINSI';
  499. *          'SI' WTAB.'CONSOLIDATION' ;WTAB.'FREA1' = FF1 ; 'FINSI';
  500.         'ERREUR' 997 ;
  501.       'FINSI';
  502.  
  503. *     Remplissage de la table ESTIMATION : champs en fin de pas
  504.       DEP0= conti.'DEPLACEMENTS';
  505. *     DEP = ENLEVER DEP0 'LX';
  506.       estim.'DEPLACEMENTS'= TT.'ZU1' ;
  507. *     'DETRUIRE' DEP;
  508.       estim.'DEFORMATIONS'= TT.'DEFF';
  509.       estim.'CONTRAINTES' = TT.'SIGF'  ;
  510.       estim.'DFGRAD'      = TT.'DFGRAD';
  511.       'SI' WTAB.'IPLAVI' ;
  512.         estim.'VARIABLES_INTERNES'        = TT.'VARF';
  513.         estim.'DEFORMATIONS_INELASTIQUES' = TT.'DFPF';
  514.       'FINSI' ;
  515.       'SI' WTAB.'GRANDS_DEPLACEMENTS' ;
  516.         estim.'FOR' = TT.'FOR' ;
  517.       'FINSI' ;
  518.       'SI'  WTAB . 'DYNAMIQUE' ;
  519.         estim.'VITESSES'      = TT.'VITI' ;
  520.         estim.'ACCELERATIONS' = TT.'ACCEI' ;
  521.       'FINSI';
  522. *       'SI'  WTAB .'ITCAR';
  523. *          estim.'CARA_ACTU' = TT.'CARA' ;
  524. *       'FINSI';
  525. ***      Reevaluation des donnees du calcul dependant de la mecanique
  526.       'SI' WTAB.'PROCEDURE_REEV_MEC' ;
  527.         REEV_MEC PRECED 1;
  528.       'FINSI' ;
  529.     'FINSI';
  530. *#######################################################################
  531. *################### FIN DU CALCUL D'UN PAS MECANIQUE ##################
  532. *#######################################################################
  533.  
  534.  
  535.  
  536. ***
  537. *** Test de la convergence mecanique thermique
  538. ***
  539.     'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  540. *      je cherche la convergence thermo-mecanique : je dois au moins
  541. *      calculer deux fois mon probleme
  542.       'SI' (I_BOTH > 1) ;
  543. *       Test sur toutes les composantes (notamment, si DIFFUSION)
  544.         CHT1     = WTAB . 'THER_COURANT' ;
  545.         DUABS2   = 'ABS' (WTAB.'THER_COURANTP' - CHT1) ;
  546.         CRICOH1  = 0. ;
  547.         NBINCO1  = 'DIME' WTAB . 'PRIM_TOT' ;
  548.         'REPE' BINCO1 NBINCO1 ;
  549.           LINCOI1  = 'MOTS' ('EXTR' WTAB.'PRIM_TOT' &BINCO1) ;
  550.           MDUABSI1 = 'MAXI' DUABS2 'AVEC' LINCOI1 ;
  551.           NORMI1   = 'MAXI' 'ABS' CHT1 'AVEC' LINCOI1 ;
  552.           'SI' (NORMI1 '&lt;EG' (1.E-14 * MDUABSI1)) ;
  553.             NORMI1   = 1. ;
  554.           'FINS' ;
  555.           CRICOHI1 = MDUABSI1 / NORMI1 ;
  556.           'SI' (CRICOHI1 '>' CRICOH1) ;
  557.             CRICOH1  = CRICOHI1 ;
  558.          'FINS' ;
  559.         'FIN' BINCO1 ;
  560.  
  561.         WTAB.'CONV_THMEC' = CRICOH1 '<' (WTAB . 'CRITERE_COHERENCE') ;
  562. *       Fin de boucle Si ne converge Pas
  563.         'SI'('NON' WTAB.'CONV_THMEC') ;
  564.           'SI' (I_BOTH >EG WTAB.'BOU_METH');
  565.             'MESS' 'PAS DE COHERENCE TROUVEE MECANIQUE THERMIQUE' ;
  566.             'ERREUR' 999;
  567.           'FINSI' ;
  568.         'FINSI' ;
  569.       'SINON' ;
  570. *       je suis sur la premiere boucle de convergence thermomecanique
  571.         WTAB.'CONV_THMEC' = FAUX ;
  572.       'FINSI' ;
  573. *     THER_COURANTPrecedent dans WTAB pour test CONVERGENCE_MEC_THE
  574.        WTAB . 'THER_COURANTP' = WTAB . 'THER_COURANT' ;
  575.     'SINON' ;
  576. *     ici je ne cherche pas la convergence thermomecanique
  577.       WTAB.'CONV_THMEC' = VRAI ;
  578.     'FINSI' ;
  579.  
  580.     'SI' (WTAB.'CONV_THMEC') ;
  581.        'QUITTER' BO_BOTH ;
  582.     'FINSI' ;
  583.  
  584.     'SI' ('NEG' ('TYPE'  WTAB.'GE0_DEB') 'MOT') ;
  585.       FORME WTAB.'GE0_DEB';
  586.     'FINSI' ;
  587.  
  588.   'FIN' BO_BOTH ;
  589.  
  590. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  591. *++++++++++++++++ fin de boucle interne de convergence mec_ther ++++++++
  592. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  593.  
  594. *----- Preparation du pas suivant : Thermique et Diffusion
  595.   'SI' WTAB.'FOR_THER';
  596.     WTAB.'TET2' = estim.'TEMPERATURES' ;
  597.     WTAB.'TET1'=WTAB.'TET2';
  598.     'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'DUPONT');
  599.       WTAB.'TET0'              = WTAB.'TET1'              ;
  600.       WTAB.'VIEUXPAS'          = WTAB.'NOUV_PAS'          ;
  601.       WTAB.'RELAXATION_DUPONT' = CHTER.'RELAXATION_DUPONT';
  602.       'SI' ('EXIS' CHTER  'INITIAL(1)') ;
  603.         WTAB.'TET1' = CHTER.'INITIAL(1)';
  604.       'FINSI' ;
  605.     'FINSI';
  606.   'FINSI' ;
  607.  
  608.   'SI' WTAB.'FOR_DIFF';
  609.     WTAB.'CO1' = WTAB.'CO2'             ;
  610.     WTAB.'CO2' = estim.'CONCENTRATIONS' ;
  611.   'FINSI' ;
  612.  
  613. *----- Preparation du pas suivant : Mecanique
  614.   'SI' WTAB.'FOR_MECA' ;
  615.     'SI' ('EXISTE' TT 'SECOND_MEMBRE');
  616.       WTAB . 'SECOND_MEMBRE' = TT . 'SECOND_MEMBRE' ;
  617.     'FINSI' ;
  618.     'SI' ('EXISTE' TT 'SUCCES');
  619.       WTAB . 'SUCCES' = TT . 'SUCCES' ;
  620.     'FINSI' ;
  621.     'SI' ('EXISTE' TT 'AUTOCOEF');
  622.       WTAB . 'AUTOCOEF' = TT . 'AUTOCOEF' ;
  623.     'FINSI' ;
  624.     'SI' ('EXISTE' TT 'AUTOREDU');
  625.       WTAB . 'AUTOREDU' = TT . 'AUTOREDU' ;
  626.     'FINSI' ;
  627.     'SI' ('EXISTE' TT 'AUTORED1');
  628.       WTAB . 'AUTORED1' = TT . 'AUTORED1' ;
  629.     'FINSI' ;
  630.     'SI' ('EXISTE' TT 'RIBLO_M');
  631.       WTAB . 'RIBLO_M' = TT . 'RIBLO_M' ;
  632.       WTAB . 'LISEA_M' = TT . 'LISEA_M' ;
  633.     'FINSI' ;
  634.     'SI' ('EXISTE' TT 'INCREMENT');
  635.       WTAB . 'INCREMENT' = TT . 'INCREMENT' ;
  636.     'FINSI' ;
  637.     'SI' ('EXISTE' TT 'INITEMPS');
  638.       WTAB . 'INITEMPS' = TT . 'INITEMPS' ;
  639.     'FINSI' ;
  640.     'SI' ('EXISTE' TT 'DT');
  641.       WTAB . 'DT' = TT . 'DT' ;
  642.     'FINSI' ;
  643.     'SI' ('EXISTE' TT 'NBPLAS');
  644.       WTAB . 'NBPLAS' = TT . 'NBPLAS' ;
  645.     'FINSI' ;
  646.     'SI' ('EXISTE' TT 'NOMBRE_ITERATIONS');
  647.       WTAB . 'NOMBRE_ITERATIONS' = TT . 'NOMBRE_ITERATIONS' ;
  648.     'FINSI' ;
  649.     'SI' ('EXISTE' TT 'FFROT');
  650.       WTAB . 'FFROT' = TT . 'FFROT' ;
  651.     'FINSI' ;
  652.     'SI' ('EXISTE' TT 'RESIDU');
  653.       WTAB . 'RESIDU' = TT . 'RESIDU' ;
  654.     'FINSI' ;
  655.     'SI' ('EXISTE' TT 'FNONL');
  656.       WTAB . 'FNONL' = TT . 'FNONL' ;
  657.     'FINSI' ;
  658.     'SI' ('EXISTE' TT 'DTPREC');
  659.       WTAB . 'DTPREC' = TT . 'DTPREC' ;
  660.     'FINSI' ;
  661. *  SP : Modification de WTAB.'TET2' suite a appel a UNPAS
  662. *  CB : SI (IPILOT) ==> TET2 est TIRE entre TEMPS0 et TI
  663.     'SI' ('EXISTE' TT 'TET2');
  664.       WTAB . 'TET2' = TT . 'TET2' ;
  665.     'FINSI' ;
  666.     'SI' ('EXISTE' TT 'XDENO');
  667.       WTAB . 'XDENO'  = TT . 'XDENO'  ;
  668.     'FINSI' ;
  669.     'SI' ('EXISTE' TT 'XDENOM');
  670.       WTAB . 'XDENOM' = TT . 'XDENOM' ;
  671.     'FINSI' ;
  672.     'SI' ('EXISTE' TT 'LASTKTAN');
  673.       WTAB . 'LASTKTAN' = TT . 'LASTKTAN' ;
  674.     'FINSI' ;
  675.     'SI' ('EXISTE' TT 'DFGRAD');
  676.       WTAB . 'DFGRAD'=  TT .'DFGRAD';
  677.     'FINSI' ;
  678.     'SI' ('EXIS' TT 'ETHER2');
  679.       WTAB.'ETHER1'= TT.'ETHER2';
  680.       WTAB.'ETHER2'= TT.'ETHER2';
  681.     'FINSI';
  682.     'SI' ('EXIS' TT 'MSRTHER2');
  683.       WTAB.'MSRTHER1'= TT.'MSRTHER2';
  684.       WTAB.'MSRTHER2'= TT.'MSRTHER2';
  685.     'FINSI';
  686.     'SI' ('EXIS' TT 'DEFOR2');
  687.       WTAB.'DEFOR1'= TT.'DEFOR2';
  688.       WTAB.'DEFOR2'= TT.'DEFOR2';
  689.     'FINSI';
  690.     'SI'('EXIS' TT 'ETAT2');
  691.       WTAB.'ETAT1'= TT.'ETAT2';
  692.     'FINSI';
  693.  
  694. * fin des reactualisations.
  695. *--------------------------- deplacement ------------------------------
  696.     TT.'DEPTOT' = estim.'DEPLACEMENTS' ;
  697. *-------------------- forces de reaction ------------------------------
  698.     TT.'REAC'= 'REAC' TT.'DEPTOT' WTAB.'CLIM';
  699.     estim.'REACTIONS' = TT.'REAC' ;
  700. *------------------ deformation d'origine thermique ------------------
  701.     'SI' WTAB.'CHAR_THE';
  702.       WTAB.'ETHER2' = TT.'ETHER2' ;
  703.       'SI' WTAB.'POR1' ;
  704.         WTAB.'MSRTHER2' = TT.'MSRTHER2';
  705.       'FINSI';
  706.     'FINSI';
  707. *------- en dynamique : vitesse, forces de reactions, operateur -------
  708.  
  709.     'SI' WTAB . 'DYNAMIQUE' ;
  710.       'SI' WTAB . 'LIAISON_PERSISTANTE' ;
  711. * corection des vitesses pour avoir des vitesses relatives nulles
  712. * aux points qui sont en contact
  713.         TT.'VITI' =  TT.'VITI' +  TT.'VAUX' ;
  714.       'SINON' ;
  715. * forces d'acceleration et d'amortissement a la fin du pas    *
  716.         WTAB.'FREA1'  = TT.'FOPL' - TT.'REAC';
  717. *     en cas d'appuis unilateraux avec choc elastique essai de corriger
  718. *     les vitesses fournies par le schema
  719.  
  720.         TT.'VITI'= VITEUNIL TT.'ZRAIDV' WTAB.'MASSE' TT.'VITI'
  721.         tt.'DEPT'   DEP0 WTAB.'DT'  WTAB.'SECOND_MEMBRE' WTAB;
  722.         'SI' ( 'EXIS' WTAB 'RATE_VITEUNIL');
  723.            'ERREUR' 'Probleme dans viteunil' ; 'FINSI';
  724.       'FINSI' ;
  725.       'SI' ( 'NEG' WTAB.'REAPREC' 'INCONNU');
  726.         'SI' ('EGA' ('TYPE' (EXTR TT.'REAC' 'MAIL')) 'MAILLAGE');
  727.           reacdif='REDU' WTAB.'REAPREC' (EXTR TT.'REAC' 'MAIL');
  728.           TT.'REAC'=  TT.'REAC' - REACDIF;
  729.         'FINSI' ;
  730.         WTAB.'REAPREC'=TT.'REAC';
  731.       'SINON';
  732.         WTAB . 'REAPREC'= TT.'REAC';
  733.       'FINSI';
  734.       estim.'VITESSES'  = TT.'VITI' ;
  735.       estim.'REACTIONS' = TT.'REAC' ;
  736.     'FINSI';
  737.  
  738. *---------- en consolidation : operateur ------------------------------
  739.     'SI' WTAB . 'CONSOLIDATION' ;
  740.       WTAB.'OPERATEUR'  = ROP;
  741.     'FINSI';
  742.  
  743. *-------------------------- configuration -----------------------------
  744.  
  745.     'SI'   WTAB.'GRANDS_DEPLACEMENTS' ;
  746.       WTAB.'RECALCUL'= VRAI     ;
  747.       WTAB.'GE0_DEB' = TT.'FOR' ;
  748.       'SI' WTAB.'ITCAR' ;
  749.         WTAB.'CARA' = TT.'CARA' ;
  750.       'FINSI' ;
  751.     'FINSI';
  752.     'FORM' WTAB.'FOR0';
  753.  
  754. ***********************************************************************
  755.     'SI' ('ET' (WTAB.'FEFP_FORMULATION')
  756.                 (WTAB.'UPDATE_LAGRANGIAN')) ;
  757. * tc            'SI' ( 'NEG' WTAB.'FOR'  'INCONNU'); DETR GEOM1; 'FINSI';
  758.       WTAB.'GE0_DEB' = TT.'FOR' ; WTAB.'RECALCUL'=VRAI;
  759.     'FINSI';
  760. ************************************************************************
  761.  
  762. *--------------------- destruction de l operateur ----------------------
  763.  
  764.     'SI' (WTAB . 'RECADET') ;
  765. *     'DETR' RH 'ELEM';
  766.       'OUBLIER'   ROP;
  767.     'FINSI';
  768.  
  769. *----------------------------- temps initial --------------------------
  770.     DT2 =( TI - TEMP0)/10000.;
  771.     'SI' (WTAB . 'AUTOMATIQUE' 'OU' ('NON' WTAB.CONV));
  772.       TEMP0 = ( TI - TEMP0) *  WTAB.'AUTOCOEF' + TEMP0;
  773.       'MESS' ' Temps effectivement calcule et non converge : ' TEMP0;
  774.       WTAB.'TEMPS0' = TEMP0;
  775.     'SINON';
  776.       TEMP0=TI;
  777.       WTAB.'TEMPS0' = TEMP0;
  778.     'FINSI';
  779.  
  780. *----------------- Appel a une procedure personnelle -------------------
  781.  
  782.     'SI' (WTAB.'PROCEDURE_PERSO1');
  783.       PRECED.'PERSO1_APPEL' = 2 ;
  784.       PERSO1 PRECED;
  785.       'SI' ('EXIS' PRECED 'ARRET') ;
  786.         WTAB . 'ARRET' = PRECED . 'ARRET' ;
  787.       'FINSI' ;
  788.       TEMP0 =  WTAB.'TEMPS0' ;
  789.     'FINSI';
  790. *----------------- On redefinit WTAB.'CHARGEMENT_I' ------------------
  791.     'SI'   WTAB.'CONVERGENCE_MEC_THE' ;
  792.       WTAB.'CHARGEMENT_I' = WTAB.'CHARGEMENT' ;
  793.     'FINSI';
  794.   'FINSI';
  795.  
  796. *++++++++++++++++++ fin preparation pas suivant mecanique +++++
  797. *##################  Sauvegarde ###############################
  798.   'SI' ('NON' WTAB.'FOR_MECA') ;
  799.     DT2 =( TI - TEMP0)/10000.;
  800.     TEMP0=TI;
  801.     WTAB.'TEMPS0' = TEMP0;
  802.     TT = 'TABLE' ;
  803.   'FINSI';
  804. * pour ne pas sauvegarder la matrice de raideur si on doit la recalculer
  805.   'SI' ('EXIS' WTAB 'RRRR')        ;
  806.     'SI'  WTAB.'RECARI'            ;
  807.        WTAB ='ENLE' WTAB 'RRRR'    ;
  808.        PRECED.'WTABLE'= WTAB       ;
  809.       'SI' ('EXIS' estim 'FOR')    ;
  810.          estim ='ENLE' estim 'FOR' ;
  811.          PRECED.'ESTIMATION'=estim ;
  812.       'FINSI';
  813.     'FINSI';
  814.   'FINSI';
  815.  
  816. *-- remplissage de la table de continuation et des tables de resultats
  817. *-- s'il faut garder les resultats, sauvegarde incrementale et
  818. *-- fantomisation si economique
  819.   ISAUV = PAS_RESU TEMP0 iauto PRECED ;
  820.  
  821. *******************************************
  822. *****   cas du pilotage automatique *******
  823. *******************************************
  824.  
  825.   'SI' WTAB . 'AUTOMATIQUE'; IAUTO=IAUTO+1;
  826.     'SI' ((WTAB . 'AUTOCOEF'   'EGA'  0.9999999D0) 'OU'
  827.           (TEMP0 '>EG' TI)) ;
  828.           TEMP0=TI;
  829.         'MESS' ' Objectif atteint ';
  830.         'MESS' ' ';
  831.         WTAB.'AUTODEUX' = FAUX ;
  832.         WTAB.AUTORED1 = 0 ;
  833.         WTAB.AUTOREDU = 1. ;
  834.         IAUTO=1;
  835.     'SINON';
  836.         'MESS' ' Objectif non atteint: On recommence';
  837.         'MESS' '  ';
  838.         'SI' ( IAUTO >EG WTAB.'AUTOPAS');
  839.          'ERREUR' 1000;
  840. *       'MESS' 'Pilotage non converge  apres ' IAUTO ' pas.';
  841.         'FINSI';
  842.         WTAB.'AUTODEUX'=VRAI;ICHG=ICHG - 1;
  843.     'FINSI';
  844.   'FINSI';
  845. *
  846.   'SI'(WTAB.'CHAR_THE' 'OU' WTAB.'FOR_THER');
  847.     WTAB.'TET1'=WTAB.'TET2';
  848.     'SI' WTAB.'FOR_MECA';
  849.       WTAB.'ETHER1' =WTAB.'ETHER2';
  850.     'FINSI';
  851.     'SI' WTAB.'POR1';
  852.       WTAB.'MSRTHER1'=WTAB.'MSRTHER2';
  853.     'FINSI';
  854.   'FINSI' ;
  855.   'SI' WTAB.'CHAR_DEFI' ; WTAB.'DEFOR1'=WTAB.'DEFOR2';'FINSI';
  856.   'SI' WTAB.'FOR_DIFF'  ; WTAB.'CO1'   =WTAB.'CO2'   ;'FINSI' ;
  857. ************************************************************************
  858. ******  cas de la convergence forcee
  859. ************************************************************************
  860.   'SI' WTAB.'PAS_AJUSTE';
  861.     'SI' ((TI < TIV) 'OU' ('NON' WTAB.'CONV')); ICHG=ICHG-1;'FINSI';
  862.   'FINSI';
  863.   'SI'  WTAB.'ARRET';
  864.     'QUITTER' BEXTERN;
  865.   'FINSI';
  866. ************************************************************************
  867. *************** FIN DE LA BOUCLE SUR LES PAS DE TEMPS ******************
  868. ************************************************************************
  869. 'FIN'  BEXTERN;
  870. ************************************************************************
  871.  
  872. 'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  873.   WTAB.'AMORTISSEMENT' = TOTO ;
  874. 'FINSI';
  875.  
  876. 'SI' ('NEG' ('TYPE'  WTAB.'GE0_DEB') 'MOT') ;
  877.   'OUBL' WTAB 'GE0_DEB' ;
  878. 'FINSI' ;
  879.  
  880. * on tue WTAB.'MAT1' s'il existe car il peut etre tres gourmand en place
  881. WTAB.'MAT1' = 'MOT' 'INCONNU';
  882.  
  883. 'MESS' ' ';
  885. '------------------- FIN DE LA PROCEDURE PASAPAS ---------------------';
  886. 'SAUT' 1 'LIGNE';
  887. *
  888. ************************************************************************
  889. *
  890. 'SI' ('EXIS' WTAB 'RRRR') ;
  891. * pour forcer le calcul de la rigidite si on  rappele  pasapas
  892.   WTAB1= 'ENLEVER' WTAB 'RRRR';
  893.   PRECED .'WTABLE'= WTAB1;
  894. 'FINSI';
  895.  
  896. * Typage de la table en sortie de PASAPAS ;
  897. PRECED.'SOUSTYPE' = 'MOT' 'PASAPAS' ;
  898.  
  899. 'FINPROC' PRECED ;
  900.  
  901.  

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