Télécharger la plaquette Cast3M


fr en

Télécharger pasapas.procedur

Retour à la liste

Numérotation des lignes : 

  1. * PASAPAS   PROCEDUR  MB234859  24/09/11    21:15:03     12008          
  2.  
  3. 'DEBP' PASAPAS  PRECED*'TABLE';
  4. 'SAUT' 1 'LIGNE';
  6. '------------------ DEBUT DE LA PROCEDURE PASAPAS --------------------';
  7. ************************************************************************
  8. *                                                                      *
  9. *   PRECED      objet de type TABLE                                    *
  10. *               C'est une table donnant , les arguments d'entrees, les *
  11. *               options, puis les  resultats des pas precedents        *
  12. *                                                                      *
  13. * Etat initial au debut du pas de temps decrit par le MCHELM WTAB.ETAT1*
  14. * champs de materiaux :                                                *
  15. *     WTAB.'MAT_REF' materiau etat a la temperature de REFERENCE       *
  16. *     WTAB.'MAT0'    materiau etat debut de calcul                     *
  17. *     WTAB.'MAT1'    materiau fin du pas de temps                      *
  18. *                                                                      *
  19. ************************************************************************
  20.  
  21. ****************   I PARTIE   ******************************************
  22. * Verifications des infos transmises et valeurs par defaut
  23. PAS_DEFA PRECED ;
  24. * Creation des tables de sortie et initialisations
  25. PAS_INIT PRECED ;
  26.  
  27. WTAB  = PRECED.'WTABLE'      ;
  28. conti = PRECED.'CONTINUATION';
  29. estim = PRECED.'ESTIMATION'  ;
  30. *-----------------------------------------------------------------------
  31. *              Appel aux procedures utilisateurs
  32. *-----------------------------------------------------------------------
  33. 'SI' WTAB.'PROCEDURE_PERSO1';
  34.   PRECED.'PERSO1_APPEL' = 1 ;
  35.   PERSO1 PRECED ;
  36.   'SI' ('EXIS' PRECED 'ARRET') ;
  37.     WTAB.'ARRET' = PRECED.'ARRET' ;
  38.   'FINSI' ;
  39. 'FINSI';
  40. *
  41. 'SI' WTAB.'PROCEDURE_REEV_MEC' ;
  42.   REEV_MEC PRECED 0 ;
  43. 'FINSI' ;
  44. *
  45. 'SI' WTAB.'PROCEDURE_REEV_THE' ;
  46.   REEV_THE PRECED 0 ;
  47. 'FINSI' ;
  48. *-----------------------------------------------------------------------
  49. *              Initialisation ETAT et MATERIAUX
  50. *-----------------------------------------------------------------------
  51. TEMP0 = conti.'TEMPS';
  52. WTAB.'ETAT0'   = PAS_ETAT PRECED TEMP0                    ;
  53. WTAB.'ETAT1'   = WTAB.'ETAT0'                             ;
  54. WTAB.'MAT0'    = PAS_MATE PRECED WTAB.'ETAT0'             ;
  55. WTAB.'MAT1'    = WTAB.'MAT0'                              ;
  56. *
  57. *---------   initialisation deformation thermique de reference
  58. 'SI' WTAB.'FOR_MECA' ;
  59.   'SI' WTAB.'ITHER' ;
  60.     'SI' ('EGA' WTAB.'ETREF' 'INCONNU');
  61.       WTAB.'ETAT_REF'= PAS_ETAT PRECED TEMP0 'REFERENCE';
  62.       WTAB.'MAT_REF' = PAS_MATE PRECED WTAB.'ETAT_REF'  ;
  63.       WTAB.'ETREF'   ='EPTH' WTAB.'MO_TOT' WTAB.'MAT_REF' WTAB.'TEMPERATURE_REFERENCE' ;
  64.     'FINS' ;
  65. *
  66.     'SI' ('EGA' ('TYPE' WTAB.'TET1') 'CHPOINT');
  67.       TET1 ='CHAN' 'CHAM' WTAB.'TET1' WTAB.'MO_TOT' 'STRESSES' 'TEMPERATURES' ;
  68.     'FINS';
  69.     ETT0  = 'EPTH' WTAB.'MO_TOT' WTAB.'MAT0' TET1 ;
  70.     WTAB.'ETHER1' = ETT0 - WTAB.'ETREF' ;
  71. *
  72.     'SI' WTAB.'POR1';
  73.       WTAB.'MSRTHER1'= PAS_EPTH PRECED WTAB.'MO_TOT' WTAB.'MAT0' TET1 ;
  74.     'FINS' ;
  75.   'FINS' ;
  76. 'FINS' ;
  77.  
  78. **************************   II PARTIE  *************************
  79. ******************************************************************
  80. ******************* boucle sur les pas de temps  *****************
  81. ******************************************************************
  82. * ------------quelques initialisations de l'algorithme
  83. IPAS  = 0 ;
  84. IAUTO = 1 ;
  85. ISAUV = VRAI;
  86. DTV   = ('EXTR' WTAB.'TEMPS_CALCULES' 1) - TEMP0 ;
  87. 'SI' ( DTV '<' 1.e-15 );
  88.   'SI' WTAB.'FREQUENTIEL' ;
  89.     DTV = 1.e-14 ;
  90.   'SINON' ;
  91.     DTV =('EXTR' WTAB.'TEMPS_CALCULES' 2) - TEMP0 ;
  92.   'FINSI';
  93. 'FINSI';
  94. * Pour l'option PAS_AJUSTE
  95. DT_AVANT = DTV ;
  96. IREREDU  = VRAI;
  97. PASFINAL = 0 ;
  98. *-----------------------------------------------------------------------
  99. *              Debut de la boucle sur les pas de temps
  100. *-----------------------------------------------------------------------
  101. 'REPETER' BEXTERN  ;
  102. *
  103. * ----- Determiner le nouvel instant a calculer TI
  104. *
  105.   TEMP0 = conti.'TEMPS' ;
  106.   NPAS  = 'DIME' WTAB.'TEMPS_CALCULES';
  107.   ICHG  = WTAB.'ICHG'; ICHG = ICHG + 1; WTAB.'ICHG' = ICHG;
  108.   'SI' ( ICHG '>' NPAS) ; 'QUITTER' BEXTERN; 'FINSI';
  109.   TI = 'EXTRAIRE'  WTAB.'TEMPS_CALCULES' ICHG ;
  110.   ITERTPS = (TI '<' TEMP0 ) ;
  111.   'SI' ('NON' WTAB.'FREQUENTIEL') ;
  112.     ITERTPS = ITERTPS 'OU' (TI 'EGA' PRECED.'TEMPS'. 0) ;
  113.     'SI' ('NON' WTAB.'REEQUILIBRAGE') ;
  114.       ITERTPS = ITERTPS 'OU' ('EGA' TI TEMP0) ;
  115.     'FINSI';
  116.   'FINSI';
  117.   'SI' ITERTPS ;
  118.     'ITERER' BEXTERN ;
  119.   'FINSI';
  120. *
  121. * ----- Ajustement du pas en fonction du nb d iteration du pas precedent
  122. *
  123.   'SI' WTAB.'PAS_AJUSTE';
  124.     PASFINAL=PASFINAL+1;
  125.     TIV=TI;
  126.     DTV=DT_AVANT;
  127.     'SI' WTAB.'CONV';
  128.       'SI' (WTAB.'NOMBRE_ITERATIONS' > 20);
  129.         DTV=DTV'/' (2. ** 0.5);
  130.       'FINSI';
  131.       'SI' (WTAB.'NOMBRE_ITERATIONS' <  7);
  132.         'SI' IREREDU;
  133.           DTV=DTV '*' (2. ** 0.5);
  134.         'FINSI';
  135.       'FINSI';
  136.       DT_AVANT= DTV;
  137.       IREREDU=VRAI;
  138.     'SINON';
  139.       DTV=0.0000000001D0* DT_AVANT;
  140.       DT_AVANT= DT_AVANT /2.;
  141.       IREREDU=FAUX;
  142.     'FINSI';
  143.     TTI = DTV* 1.0000000001 + TEMP0;
  144.     'SI' ( TTI '&lt;EG' TIV );
  145.       TI=TTI;
  146.       PASFINAL=0;
  147.     'FINSI';
  148.     'SI' (PASFINAL '>EG' 1.5) ; 
  149.       DT_AVANT = TI-TEMP0; 
  150.     'FINSI';
  151.   'FINSI';
  152. *
  153.   estim.'TEMPS' = TI;
  154.   WTAB.'DT'     = TI '-' TEMP0 ;
  155.   'SI'((WTAB.'CONSOLIDATION' 'OU' WTAB.'DYNAMIQUE') 'ET' (WTAB.'DT' 'EGA' 0.));
  156.     WTAB.'DT' = WTAB.'DTPREC' ;
  157.     'ITERER' BEXTERN;
  158.   'FINSI';
  159. *
  160. * ----- Compteur IPAS des pas de temps calcules
  161.   IPAS = IPAS + 1;
  162.   WTAB.'PAS'=IPAS;
  163.  
  164. * ----- Mise a jour du modele si chargement de type MODE :
  165.   'SI' WTAB.'CHAR_MODE' ;
  166.     MOD1P     =  WTAB.'MODELE' ;
  167.     MODEVOL1  = 'TIRE' PRECED.'MODELE' 'MODE' TI ;
  168.     WTAB.'MODELE' = MODEVOL1 ;
  169.  
  170.     'SI' WTAB.'CHAR_MATE' ;
  171.       MATEVOL1  = 'TIRE' PRECED.'CARACTERISTIQUES' 'MATE' TI ;
  172.       WTAB.'CARACTERISTIQUES' = MATEVOL1 ;
  173.     'FINSI' ;
  174.  
  175.     PAS_MODL PRECED ;
  176.  
  177. *   MaJ table rayonnement :
  178.     'SI' WTAB.'RAYO' ;
  179.       WTAB.'RAYONNEMENT' = 'PRRA' WTAB.'MOD_RAY' ;
  180.     'FINSI' ;
  181.  
  182. *   On sauvegarde les modeles successifs pour analyse resultats :
  183.     ESTIM.'MODELES' = WTAB.'MODELE' ;
  184.  
  185. *   Initialisation TINI si thermique :
  186.     'SI' WTAB.'THE1' ;
  187.       MODX      = WTAB.'MOD_COND' ;
  188.       MODXP     = MOD1P 'EXTR' 'MATE' 'CONDUCTION' ;
  189.       GEO1      = MODX  'EXTR' 'MAIL' ;
  190.       GEO1P     = MODXP 'EXTR' 'MAIL' ;
  191.       DGEO1     = GEO1  'DIFF' GEO1P ;
  192.  
  193. *     Adapation du modele : transport de la solution
  194.       'SI' ((WTAB.'ADAP_MODE') 'ET' ('NON' (TI '<' (WTAB.'TPS1_ADAP')))) ;
  195.         mess '*****  Transport solution adaptation' ;
  196.         IADAP1    = (WTAB.'INDI_ADAP') + 1 ;
  197.         'SI' (IADAP1 '&lt;EG' ('DIME' (WTAB.'LTPS_ADAP'))) ;
  198.           WTAB.'TPS1_ADAP' = WTAB.'LTPS_ADAP' 'EXTR' IADAP1 ;
  199.           WTAB.'INDI_ADAP' = IADAP1 ;
  200.         'SINON' ;
  201.           WTAB.'TPS1_ADAP' = 'MAXI' WTAB.'TEMPS_CALCULES' + 1. ;
  202.         'FINSI' ;
  203.  
  204. *       Projection du champ de temperature sur la partie adaptee :
  205.         DGEO1P    = DGEO1 'INTE' GEO1P ;
  206.         DGEO2P    = (DGEO1 'INTE' GEO1) 'INCLU' DGEO1P 'VOLU' 'BARY'  ;
  207.         CHT1P     = CONTI.'TEMPERATURES' 'EXCO' 'T' 'T' ;
  208.         CHAMT1P   = 'CHAN' 'CHAM' DGEO1P CHT1P ;
  209.         CHT2P     = 'PROI' CHAMT1P  DGEO2P ;
  210.         CHT1P     = ('REDU' CHT1P (GEO1P 'DIFF' DGEO1P)) 'ET' CHT2P ;
  211.         CONTI.'TEMPERATURES' = CHT1P ;
  212.  
  213. *       Calcul de la capacite courante sur la partie adaptee :
  214.         DMODX     = 'REDU' MODX DGEO2P ;
  215.         DMATX     = 'REDU' (WTAB.'CARACTERISTIQUES') DMODX ;
  216.         DMATX     = 'VARI' 'NUAG' DMODX DMATX CHT2P ;
  217.         CAPAX     = 'CAPA' DMODX DMATX ;
  218.         CAPA1P    = ('REDU' (WTAB.'CAPA_COURANTE') (GEO1P 'DIFF' DGEO1P)) 'ET' CAPAX ;
  219.         WTAB.'CAPA_COURANTE' = CAPA1P ;
  220.  
  221.         DGEO1     = (DGEO1 'INTE' GEO1) 'DIFF' DGEO2P ;
  222.       'FINSI' ;
  223.  
  224. *     Calcul temperature dans partie apportee du modele :
  225.       'SI' ('NON' (('VIDE' DGEO1) 'OU' ('VIDE' (DGEO1 'INTE' GEO1)))) ;
  226. *       Quantite de chaleur sur geometrie GEO1P :
  227.         CHT1P     = CONTI.'TEMPERATURES' 'EXCO' 'T' 'T' ;
  228.         CAPA1P    = WTAB.'CAPA_COURANTE' lump ;
  229.         CHQ1P     = CAPA1P * CHT1P ;
  230. *       Quantite de chaleur dans apport :
  231.         DMODX     = 'REDU' MODX DGEO1 ;
  232.         DMATX     = 'REDU' (WTAB.'CARACTERISTIQUES') DMODX ;
  233.         CHTINIX   = DMATX 'EXCO' 'TINI' 'T' ;
  234.         CHTINIX   = 'CHAN' 'CHPO' CHTINIX DMODX 'MOYE' ;
  235.         DMATX     = 'VARI' 'NUAG' DMODX DMATX CHTINIX ;
  236.         CAPAX     = ('CAPA' DMODX DMATX) lump ;
  237.         CHQX      = CAPAX * CHTINIX ;
  238. *       Calcul T sur GEO1 :
  239.         CAPA1     = CAPA1P et CAPAX ;
  240.         CHTX1     = 'MANU' 'CHPO' GEO1 1 'T' 1. 'NATURE' 'DIFFUS' ;
  241.         CHQX1     = CAPA1 * CHTX1;
  242.         CHT1      = (CHQ1P + CHQX) * ('INVE' CHQX1) ('MOTS' 'Q') ('MOTS' 'Q') ('MOTS' 'T') ;
  243. *        CHT1      = reso CAPA1 (CHQ1P + CHQX) ; comm plus lent avec RESO ;
  244.         CONTI.'TEMPERATURES' = CHT1 ;
  245.       'FINSI' ;
  246.  
  247. *      Sauvegarde TEMPERATURES au debut du pas apres APPORT :
  248.        IND1     = ('DIME' PRECED.'TEMPERATURES_APPORT') - 1 ;
  249.        CHTAPP1  = CONTI.'TEMPERATURES' 'EXCO' 'T' 'TAPP' ;
  250.        PRECED.'TEMPERATURES_APPORT' . IND1 = CHTAPP1 ;
  251.        ESTIM .'TEMPERATURES_APPORT'        = CHTAPP1 ;
  252.  
  253. *   Finsi formulation thermiquue
  254.     'FINSI' ;
  255.  
  256. *   Initialisation PROPORTIONS_PHASES si metallurgie :
  257.     'SI' WTAB.'FOR_METALLU' ;
  258.       MODX      = WTAB.'MODELE' 'EXTR' 'FORM' 'METALLURGIE' ;
  259.       MODXP     = MOD1P         'EXTR' 'FORM' 'METALLURGIE' ;
  260.       GEO1      = MODX  'EXTR' 'MAIL' ;
  261.       GEO1P     = MODXP 'EXTR' 'MAIL' ;
  262.       DGEO1     = GEO1 'DIFF' GEO1P ;
  263.       'SI' ('NON' (('VIDE' DGEO1) 'OU' ('VIDE' (DGEO1 'INTE' GEO1)))) ;
  264.         CHMET0    = CONTI.'PROPORTIONS_PHASES' ;
  265.         DMODX     = 'REDU' MODX DGEO1 ;
  266.         CHMET1    = 'ZERO' DMODX 'VARINTER' ;
  267.         LISCO     = 'EXTR' DMODX 'VARI' ;
  268.         CMP1      = 'EXCO' CHMET1 ('EXTR' LISCO 1) ('EXTR' LISCO 1) ;
  269.         CHMET1    = CHMET1 + (CMP1 + 1.D0);
  270.         CHMET1    = CHMET0 'ET' CHMET1 ;
  271.         CONTI.'PROPORTIONS_PHASES' = CHMET1 'REDU' MODX ;
  272. *trac CHMET1 MODX ;
  273.       'FINSI' ;
  274.     'FINS' ;
  275.  
  276. *   Initialisation MECA si besoin :
  277.     'SI' WTAB.'FOR_MECA' ;
  278.       MODM1    = 'VIDE' 'MMODEL' ;
  279.       MODM1P   = 'VIDE' 'MMODEL' ;
  280.       'SI'  WTAB.'MEC1' ;
  281.         MODM1    = MODM1 'ET' WTAB.'MOD_MEC' ;
  282.         MODM1P   = MODM1P 'ET' ('EXTR' MOD1P 'FORM' 'MECANIQUE') ;
  283.       'FINSI' ;
  284.       'SI'  WTAB.'POR1' ;
  285.         MODM1    = MODM1 'ET' WTAB.'MOD_POR' ;
  286.         MODM1P   = MODM1P 'ET' ('EXTR' MOD1P 'FORM' 'POREUX') ;
  287.       'FINSI' ;
  288.       'SI' ('NON' ('VIDE' MODM1)) ;
  289.         GEOM1    = 'EXTR' MODM1  'MAIL' ;
  290.         GEOM1P   = 'EXTR' MODM1P 'MAIL' ;
  291.         DGEOM1   = GEOM1 'DIFF' GEOM1P ;
  292.         'SI' ('NON' (('VIDE' DGEOM1) 'OU' ('VIDE' (DGEOM1 'INTE' GEOM1)))) ;
  293.           DMODM1   = 'REDU' MODM1 DGEOM1 ;
  294.  
  295. *         Mise a jour ETREF :
  296.           WTAB.'ETAT_REF'= PAS_ETAT PRECED TEMP0 'REFERENCE' ;
  297.           WTAB.'MAT_REF' = PAS_MATE PRECED WTAB.'ETAT_REF'   ;
  298.           TREF1    = 'EXCO' WTAB.'MAT_REF' 'TREF' 'T' ;
  299.           WTAB.'TEMPERATURE_REFERENCE' = 'CHAN' 'TYPE' TREF1 'TEMPERATURES' ;
  300.           DETREF1  = 'EPTH' DMODM1 WTAB.'MAT_REF' WTAB.'TEMPERATURE_REFERENCE' ;
  301.           WTAB.'ETREF' = WTAB.'ETREF' 'ET' DETREF1 ;
  302.  
  303. *         Mise a jour ETHER1 :
  304.           ETAT1    = PAS_ETAT PRECED TI ;
  305.           MAT1     = PAS_MATE PRECED ETAT1 ;
  306.           TET1     = ('TIRE' PRECED.'CHARGEMENT' 'TAPP' TI) 'EXCO' 'TAPP' 'T' ;
  307.           DETHER1  = ('EPTH' DMODM1 MAT1 TET1) - DETREF1 ;
  308.           WTAB.'ETHER1' = WTAB.'ETHER1' 'ET' DETHER1 ;
  309.  
  310. *         Mise a jour deplacements et deformations :
  311.           DEPDGEO1 = ('EXCO' ETAT1 ('EXTR' DMODM1 'DEPL')) 'CHAN' 'CHPO' DMODM1 ;
  312.           DEPSL1   = ('EPSI' DMODM1 DEPDGEO1) - DETHER1 ;
  313.  
  314. *         Mise a jour table estimation :
  315.           estim.DEPLACEMENTS               = conti.DEPLACEMENTS 'ET' DEPDGEO1 ;
  316.           estim.CONTRAINTES                = conti.CONTRAINTES  'ET' ('ZERO' DMODM1 'CONTRAIN') ;
  317.           estim.VARIABLES_INTERNES         = conti.VARIABLES_INTERNES  'ET' ('ZERO' DMODM1 'VARINTER') ;
  318.           estim.DEFORMATIONS_INELASTIQUES  = conti.DEFORMATIONS_INELASTIQUES  'ET' ('ZERO' DMODM1 'DEFINELA') ;
  319.           estim.DEFORMATIONS               = conti.DEFORMATIONS 'ET' ('ZERO' DMODM1 'DEFORMAT') ;
  320.  
  321.           'SI' ('EXISTE' PRECED 'DEFORMATIONS_APPORT');
  322.             estim.DEFORMATIONS_APPORT      = conti.DEFORMATIONS_APPORT 'ET' DEPSL1 ;
  323.           'FINSI' ;
  324.         'FINSI' ;
  325.  
  326.       'SINON' ;
  327.             estim.DEFORMATIONS_APPORT      = conti.DEFORMATIONS_APPORT ;
  328.       'FINSI' ;
  329.     'FINSI' ;
  330.  
  331. * Fin traitement CHAR_MODE
  332.   'FINSI' ;
  333. *
  334. * ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  335. * +++++++++++++++++ BOUCLE interne de bouclage mec_ther ++++++++++++++++
  336. * ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  337.   I_BOTH = 0 ;
  338.   'REPETER' BO_BOTH WTAB.'BOU_METH' ;
  339.     I_BOTH = I_BOTH + 1 ;
  340. *
  341. *   Re-initialisation de WTAB.'CHARGEMENT' en tete de boucle mec-ther
  342.     'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  343.       WTAB.'CHARGEMENT' = WTAB.'CHARGEMENT_I' ;
  344.     'FINSI';
  345. *
  346. *   ####################################################################
  347. *   ########### RESOLUTION D'UN PAS SI THERMIQUE/DIFFUSION #############
  348. *   ####################################################################
  349.     'SI' (WTAB.'FOR_THER' 'OU' WTAB.'FOR_DIFF') ;
  350. *
  351. *     Se placer dans la configuration adequate
  352.       'SI' WTAB.'FOR_MECA' ;
  353.         'SI' (I_BOTH > 1) ;
  354.           'SI' WTAB.'GRANDS_DEPLACEMENTS' ;
  355.             WTAB.'RECARI'   = VRAI ;
  356.             WTAB.'REA_GEOM' = VRAI ;
  357.             'FORM' WTAB.'FOR' ;
  358.           'FINSI';
  359.         'SINON' ;
  360.           'FORM' WTAB.'GE0_DEB' ;
  361.         'FINSI' ;
  362.       'FINSI' ;
  363. *
  364.       'SI' (TI 'NEG' TEMP0);
  365.         WTAB.'TEM_CALC' = 'PROG' TI;
  366.         'SI' ('EXISTE' PRECED 'PAS_MAX');
  367.           'SI' ((TEMP0 + PRECED.'PAS_MAX') '<' TI);
  368.             WTAB.'TEM_CALC' = 'PROG' TEMP0 'PAS' PRECED.'PAS_MAX' TI ;
  369.             WTAB.'TEM_CALC' = 'ENLE' WTAB.'TEM_CALC' 1;
  370.           'FINSI';
  371.         'FINSI';
  372.  
  373. *       Mise a jour blocages passes dans chargement
  374. *       -> ISO AVANT MAIS A FAIRE A CHAQUE FOIS SI DEPEND DE LA MECA ??
  375.         'SI' ('EGA' I_BOTH 1) ;
  376. *         Cas BLOCAGES_THERMIQUES (BLOT) :
  377.           'SI' WTAB.'CHAR_BLOT' ;
  378.             BLOT1 = 'TIRE' PRECED.'BLOCAGES_THERMIQUES' 'BLOT' TI ;
  379.             WTAB.'BLOCAGES_THERMIQUES' = BLOT1 ;
  380.           'FINSI' ;
  381. *
  382. *         Cas BLOCAGES_DIFFUSIONS (BLOD) :
  383.           'SI' WTAB.'CHAR_BLOD' ;
  384.             BLOD1 = 'TIRE' PRECED.'BLOCAGES_DIFFUSIONS' 'BLOD' TI ;
  385.             WTAB.'BLOCAGES_DIFFUSIONS' = BLOD1 ;
  386.           'FINSI' ;
  387.         'FINSI' ;
  388. *
  389. * ----- Appel a la procedure thermique (lineaire ou nonlineaire) ----
  390. *
  391.         'SI' ('NON' WTAB.'COUPLAGE'); 'MENAGE'; 'FINSI';
  392.         'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'NONLINEAIRE');
  393.           CHTER = TRANSNON PRECED;
  394.         'FINSI';
  395.         'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'LINEAIRE');
  396.           CHTER = TRANSLIN PRECED;
  397.         'FINSI';
  398.         'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'DUPONT');
  399.           WTAB.'NOUV_PAS' = ('EXTR' WTAB.'TEM_CALC' 1) - TEMP0;
  400.           CHTER = DUPONT2 PRECED;
  401.         'FINSI';
  402. *
  403. *       Verification d'une erreur en sortie
  404.         'SI' CHTER.'ERREUR';
  405.           'ERREUR' ' ';
  406.         'FINSI';
  407. *
  408. *       Se placer sur la configuration de reference (appel proc ut.)
  409.         'FORM' WTAB.'FOR0' ;
  410. *
  411. *        Reevaluation des donnees du calcul dependant de la temperature
  412.         'SI' WTAB.'PROCEDURE_REEV_THE' ;
  413.           REEV_THE PRECED  1 ;
  414.         'FINSI' ;
  415.       'FINS';
  416. *
  417.     'FINSI' ;
  418. *   ####################################################################
  419. *   ###### RESOLUTION D'UN PAS SI MECANIQUE/MECANIQUE DES FLUIDES ######
  420. *   ####################################################################
  421.     'SI' (WTAB.'FOR_MECA' 'OU' WTAB .'NVSTNL') ;
  422.       'SAUT' 1 'LIGN';
  423.       'MESSAGE' ' Numero du pas : ' IPAS '  Indice d evolution : ' ICHG ' -> temps :' TI;
  424. *
  425. *     Se placer sur la configuration de debut de pas
  426.       'FORM' WTAB.'GE0_DEB' ;
  427. *
  428.       'MENA' ;
  429.       IERRMEC = UNPAS PRECED;
  430. *
  431. *     Se placer sur la configuration de reference
  432.       'FORM' WTAB.'FOR0' ;
  433. *
  434. *     Test si on a converge
  435.       'SI' IERRMEC ;
  436.         'ERREUR' 997 ;
  437.       'FINSI';
  438. *
  439. *     Reevaluation des donnees du calcul dependant de la mecanique
  440.       'SI' WTAB.'PROCEDURE_REEV_MEC' ;
  441.         REEV_MEC PRECED 1;
  442.       'FINSI' ;
  443. *
  444.     'FINSI';
  445. *   ####################################################################
  446. *   ############ TEST DE LA CONVERGENCE MECANIQUE THERMIQUE ############
  447. *   ####################################################################
  448.     'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  449. *      je cherche la convergence thermo-mecanique : je dois au moins
  450. *      calculer deux fois mon probleme
  451.       WTAB.'CONV_THMEC' = FAUX ;
  452.       'SI' (I_BOTH > 1) ;
  453. *       Test sur toutes les composantes (notamment, si DIFFUSION)
  454.         CHT1     = WTAB.'THER_COURANT' ;
  455.         DUABS2   = 'ABS' (WTAB.'THER_COURANTP' - CHT1) ;
  456.         CRICOH1  = 0. ;
  457.         NBINCO1  = 'DIME' WTAB.'PRIM_TOT' ;
  458.         'REPE' BINCO1 NBINCO1 ;
  459.           LINCOI1  = 'MOTS' ('EXTR' WTAB.'PRIM_TOT' &BINCO1) ;
  460.           MDUABSI1 = 'MAXI' DUABS2 'AVEC' LINCOI1 ;
  461.           NORMI1   = 'MAXI' 'ABS' CHT1 'AVEC' LINCOI1 ;
  462.           'SI' (NORMI1 '&lt;EG' (1.E-14 * MDUABSI1)) ;
  463.             NORMI1 = 1. ;
  464.           'FINS' ;
  465.           CRICOHI1 = MDUABSI1 / NORMI1 ;
  466.           'SI' (CRICOHI1 '>' CRICOH1) ;
  467.             CRICOH1  = CRICOHI1 ;
  468.          'FINS' ;
  469.         'FIN' BINCO1 ;
  470.  
  471.         WTAB.'CONV_THMEC' = CRICOH1 '<' (WTAB . 'CRITERE_COHERENCE') ;
  472. *       Fin de boucle Si ne converge Pas
  473.         'SI'('NON' WTAB.'CONV_THMEC') ;
  474.           'SI' (I_BOTH >EG WTAB.'BOU_METH');
  475.             'MESS' 'PAS DE COHERENCE TROUVEE MECANIQUE THERMIQUE' ;
  476.             'ERREUR' 999;
  477.           'FINSI' ;
  478.         'FINSI' ;
  479.       'FINSI' ;
  480. *     THER_COURANTPrecedent dans WTAB pour test CONVERGENCE_MEC_THE
  481.       WTAB.'THER_COURANTP' = WTAB.'THER_COURANT' ;
  482.     'FINSI' ;
  483. *
  484.     'SI' WTAB.'CONV_THMEC' ;
  485.       'QUITTER' BO_BOTH ;
  486.     'FINSI' ;
  487. *
  488.   'FIN' BO_BOTH ;
  489. * ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  490. * +++++++++++++++++ BOUCLE interne de bouclage mec_ther ++++++++++++++++
  491. * ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  492. *
  493. * ----------------------------------------------------------------------
  494. *       Mise a jour de quantites de WTAB a partir de celles a TI
  495. * ----------------------------------------------------------------------
  496. *
  497. * Formulation THERMIQUE
  498.   'SI' WTAB.'FOR_THER';
  499.     'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'DUPONT');
  500.       WTAB.'TET0'              = estim.'TEMPERATURES' ;
  501.       WTAB.'VIEUXPAS'          = WTAB.'NOUV_PAS' ;
  502.       WTAB.'RELAXATION_DUPONT' = CHTER.'RELAXATION_DUPONT';
  503.     'FINSI';
  504.   'FINSI' ;
  505. *
  506. * Formulation MECANIQUE
  507.   'SI' WTAB.'FOR_MECA' ;
  508. *
  509.     'SI' WTAB.'ITHER' ;
  510.       WTAB.'TET1'   = WTAB.'TET2';
  511.       WTAB.'ETHER1' = WTAB.'ETHER2';
  512.       'SI' WTAB.'POR1';
  513.         WTAB.'MSRTHER1'= WTAB.'MSRTHER2';
  514.       'FINSI';
  515.     'FINSI' ;
  516. *
  517.     'SI' WTAB.'CHAR_DEFI' ; 
  518.       WTAB.'DEFOR1' = WTAB.'DEFOR2';
  519.     'FINSI';
  520. *
  521.     'SI' WTAB.'DYNAMIQUE' ;
  522.       'SI' ('NON' WTAB.'LIAISON_PERSISTANTE') ;
  523.         WTAB.'FREA1' = WTAB.'FOPL' ;
  524.       'FINSI';
  525.       WTAB.'REAPREC' = estim.'REACTIONS';
  526.     'FINSI';
  527. *
  528. *   Actualiser la configuration de debut de pas
  529.     'SI' ( WTAB.'GRANDS_DEPLACEMENTS' 'OU'
  530.           (WTAB.'FEFP_FORMULATION' 'ET' WTAB.'UPDATE_LAGRANGIAN')) ;
  531.       WTAB.'GE0_DEB' = WTAB.'FOR' ;
  532.     'FINSI';
  533. *
  534.   'FINSI';
  535. *
  536. * Option AUTOMATIQUE : temps effectivement calcule
  537.   'SI' (WTAB.'AUTOMATIQUE' 'OU' ('NON' WTAB.'CONV'));
  538.     TEMP0 = (TI '-' TEMP0) '*' WTAB.'AUTOCOEF' '+' TEMP0;
  539.     estim.'TEMPS' = TEMP0 ;
  540.     MSGZ = 'CHAI' ' Temps effectivement calcule';
  541.     'SI' ('NON' WTAB.'CONV');
  542.       MSGZ = 'CHAI' MSGZ ' et non converge';
  543.     'FINSI';
  544.     'MESS' MSGZ ' : 'TEMP0;
  545.   'FINSI';
  546.  
  547. * Procedure utilisateur
  548.   'SI' WTAB.'FOR_MECA' ;
  549.     'SI' WTAB.'PROCEDURE_PERSO1' ;
  550.       PRECED.'PERSO1_APPEL' = 2 ;
  551.       PERSO1 PRECED;
  552.       'SI' ('EXIS' PRECED 'ARRET') ;
  553.         WTAB.'ARRET' = PRECED.'ARRET' ;
  554.       'FINSI' ;
  555.     'FINSI';
  556. *
  557. *   Redefinir WTAB.'CHARGEMENT_I'
  558.     'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  559.       WTAB.'CHARGEMENT_I' = WTAB.'CHARGEMENT' ;
  560.     'FINSI';
  561.   'FINSI';
  562. *
  563. * pour ne pas sauvegarder la matrice de raideur si on doit la recalculer
  564.   'SI' ('EXIS' WTAB 'RRRR') ;
  565.     'SI' WTAB.'RECARI' ;
  566.       'OUBL' WTAB 'RRRR' ;
  567.     'FINSI';
  568.   'FINSI';
  569. *
  570. * ----------------------------------------------------------------------
  571. *        Remplir CONTINUATION et table de resultats + sauvegarde
  572. * ----------------------------------------------------------------------
  573.   ISAUV = PAS_RESU iauto PRECED ;
  574. *
  575. * Option AUTOMATIQUE
  576.   'SI' WTAB.'AUTOMATIQUE'; 
  577.     'SI' ((WTAB.'AUTOCOEF' 'EGA'  0.9999999D0) 'OU' (TEMP0 '>EG' TI)) ;
  578.       'MESS' ' Objectif atteint '; 'MESS' ' ';
  579.       WTAB.'AUTODEUX' = FAUX ;
  580.       WTAB.'AUTORED1' = 0 ;
  581.       WTAB.'AUTOREDU' = 1. ;
  582.       IAUTO = 1;
  583.     'SINON';
  584.       'MESS' ' Objectif non atteint: On recommence'; 'MESS' '  ';
  585.       IAUTO = IAUTO+1;
  586.       'SI' (IAUTO >EG WTAB.'AUTOPAS');
  587.         'ERREUR' 1000;
  588. *       'MESS' 'Pilotage non converge  apres ' IAUTO ' pas.';
  589.       'FINSI';
  590.       WTAB.'AUTODEUX'=VRAI;
  591.       ICHG=ICHG - 1; WTAB.'ICHG' = ICHG ;
  592.     'FINSI';
  593.   'FINSI';
  594. *
  595. * Option PAS_AJUSTE
  596.   'SI' WTAB.'PAS_AJUSTE';
  597.     'SI' ((TI < TIV) 'OU' ('NON' WTAB.'CONV')); 
  598.       ICHG=ICHG-1; WTAB.'ICHG' = ICHG ;
  599.     'FINSI';
  600.   'FINSI';
  601. *
  602.   'SI' WTAB.'ARRET';
  603.     'QUITTER' BEXTERN;
  604.   'FINSI';
  605. *
  606. 'FIN' BEXTERN;
  607. *-----------------------------------------------------------------------
  608. *              Fin de la boucle sur les pas de temps
  609. *-----------------------------------------------------------------------
  610. 'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  611.   WTAB.'AMORTISSEMENT' = TOTO ;
  612. 'FINSI';
  613. *
  614. 'SI' ('NEG' ('TYPE'  WTAB.'GE0_DEB') 'MOT') ;
  615.   'OUBL' WTAB 'GE0_DEB' ;
  616. 'FINSI' ;
  617. *
  618. * pour forcer le calcul de la rigidite si on  rappele  pasapas
  619. 'SI' ('EXIS' WTAB 'RRRR') ;
  620.   'OUBL' WTAB 'RRRR';
  621. 'FINSI';
  622. *
  623. * on tue WTAB.'MAT1' s'il existe car il peut etre tres gourmand en place
  624. WTAB.'MAT1' = 'MOT' 'INCONNU';
  625. *
  626. 'MESS' ' ';
  628. '------------------- FIN DE LA PROCEDURE PASAPAS ---------------------';
  629. 'SAUT' 1 'LIGNE';
  630. ************************************************************************
  631. 'FINPROC' PRECED ;
  632.  
  633.  
  634.  
  635.  

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