Télécharger la plaquette Cast3M


fr en

Télécharger unpas.procedur

Retour à la liste

Numérotation des lignes : 

  1. * UNPAS     PROCEDUR  FANDEUR   22/01/19    21:15:18     11256          
  2. 'DEBPROC' UNPAS PRECED*'TABLE';
  3. *
  4. *----------------------------------------------------------------------*
  5. *         PROCEDURE UNPAS                                              *
  6. *                                                                      *
  7. *    calcul d'un increment de solution en grand deplacement plastique  *
  8. *    par la methode des residus                                        *
  9. *----------------------------------------------------------------------*
  10. *    En entree
  11. *     PRECED  la table passee à PASAPAS
  12. *
  13. *  kich nota champ de materiau : etablit au debut de procedure UNPAS
  14. *  estimation des caracteristiques fin de pas (TI     ): ZMAT22  -->> ZMATXX
  15. *  caracteristiques debut de pas              (TEMPS0) : ZMAT11
  16. *  puis dans la boucle de non-convergence
  17. *      ZMAT1 initialise avec ZMAT11
  18. *  materiau fin de pas de temps : ZMAT2 sorti de COMP,
  19. *    range dans WTAB.'MAT1' en sortie unpas
  20. ************************************************************************
  21. *    sortie STAB12  indice :                                           *
  22. *                                                                      *
  23. *      DEPT   increment de deplacement  sur le pas                     *
  24. *      SIGF   contraintes a la fin du pas                              *
  25. *      VARF   variables internes a la fin du pas                       *
  26. *      DFPF   deformation inelastique a la fin du pas                  *
  27. *      CONV   logique valant vrai si pas de probleme de convergence    *
  28. *      DEFF   deformations a la fin du pas si grandes deformations     *
  29. *----------------------------------------------------------------------*
  30. *
  31. * Appel a PAS_MODL :
  32. * Mise a jour des indices de PRECED.WTABLE relatifs aux modeles
  33. * si PRECED.WTABLE.MODELE a ete modifie
  34. PAS_MODL PRECED ;
  35.  
  36. *CB215821 : Recuperation de XPETIT (07/12/2016)
  37. XPETIT   = 'VALE' 'PETI'    ;
  38. XGRAND   = 'VALE' 'GRAND'   ;
  39.  
  40. *  declanchement du recalcul de la matrice
  41.   ITRCLC =  1  * WTAB.'DELTAITER';
  42.   'SI' (ITRCLC <  20) ; ITRCLC =  20; 'FINSI';
  43.  
  44. WTAB = PRECED.'WTABLE'      ;
  45. LAG_TOT=VRAI;
  46. 'SI' ('EXIS' WTAB 'LAG_TOT'); LAG_TOT=WTAB.'LAG_TOT'; 'FINS';
  47. conti= PRECED.'CONTINUATION';
  48. estim= PRECED.'ESTIMATION'  ;
  49. *
  50. *  initialisation et reprise des valeurs des tables ****************
  51. *
  52. MXMYMZ = 'MOTS' 'MX'  'MY'  'MZ'  'MT'  'FP'  'FPQ'  'FTP'
  53.                 'IMX' 'IMY' 'IMZ' 'IMT' 'IFP' 'IFPQ' 'IFTP' ;
  54. MXMFLX = 'MOTS' 'MX'  'MY'  'MZ'  'MT'  'FLX' 'FP'  'FPQ'  'FTP'
  55.                 'IMX' 'IMY' 'IMZ' 'IMT' 'FLX' 'IFP' 'IFPQ' 'IFTP' ;
  56. MNPRIM = 'MOTS' 'UX'  'UY'  'UZ'  'UR'  'UT'   'RX'   'RY'  'RZ'  'RT'
  57.                       'P'   'PQ'  'TP'   'ALFA' 'BETA'
  58.                 'IUX' 'IUY' 'IUZ' 'IUR'  'IUT'  'IRX'  'IRY' 'IRZ' 'IRT'
  59.                 'IP' 'IPQ'  'ITP' 'IALF' 'IBET' ;
  60. MLPRIM = MNPRIM 'ET' ('MOTS' 'LX');
  61. MNDUAL = 'MOTS' 'FX'  'FY'  'FZ'  'FR'   'FT'   'MX'   'MY'  'MZ'  'MT'
  62.                       'FP'  'FPQ' 'FTP'  'FALF' 'FBET'
  63.                 'IFX' 'IFY' 'IFZ' 'IFR'  'IFT'  'IMX'  'IMY' 'IMZ' 'IMT'
  64.                 'IFP' 'IFPQ' 'IFTP' 'IFAL' 'IFBE' ;
  65. MLDUAL = MNDUAL 'ET' ('MOTS' 'FLX');
  66. MLDEPL = 'MOTS' 'UX' 'UY' 'UZ' 'UR' 'UT'  'IUX' 'IUY' 'IUZ' 'IUR' 'IUT'
  67.                 'ALFA' 'BETA' 'IALF' 'IBET' ;
  68. MLROTA = 'MOTS' 'RX' 'RY' 'RZ' 'RT' 'P' 'PQ' 'TP'
  69.                 'IRX' 'IRY' 'IRZ' 'IRT' 'IP' 'IPQ' 'ITP';
  70. MLDEFOR = 'MOTS'  'EPXX' 'EPYY' 'EPZZ' 'EPSS' 'EPTT' 'EPRR'
  71.                   'GAXY' 'GAXZ' 'GAYZ' 'GAST' 'GASN' 'GATN'
  72.                   'GARZ' 'GART' 'GAZT' 'GXY '
  73. *                 'CX  ' 'CY  ' 'CZ  '
  74.                   'EPSE' ;
  75. MVPRIM = 'MOTS' 'VTX' 'VTY' 'VTZ' 'VTR' 'VTT' 'VWX' 'VWY' 'VWZ' 'VWT'
  76.                 'VVP' 'VVPQ' 'VVTP' 'VALF' 'VBET'
  77.           'IVTX' 'IVTY' 'IVTZ' 'IVTR' 'IVTT' 'IVWX' 'IVWY' 'IVWZ' 'IVWT'
  78.                 'IVVP' 'IVPQ' 'IVTP' 'IVAL' 'IVBE' 'VLX';
  79. MOCA = 'MOTS' 'VECT' 'VX  ' 'VY  ' 'VZ  ' 'VXF ' 'VYF ' 'VZF '
  80.               'V1X ' 'V1Y ' 'V1Z ' 'V2X ' 'V2Y ' 'V2Z ' ;
  81.  
  82. *  definition de variables locales
  83. PASUNIL = FAUX;
  84. *'SI' ( 'EXIS' PRECED 'ECRIT' ) ; list wtab; 'FINS';
  85. ICERAM = WTAB.'CERAMIQUE' ;
  86. IDYN   = WTAB.'DYNAMIQUE';
  87. IELANL = WTAB.'NON_LINEAIRE';
  88. IENDOM = WTAB.'ENDOMMAGEMENT';
  89. IFEFP  = WTAB.'FEFP_FORMULATION' ;
  90. IFEFPUL= WTAB.'UPDATE_LAGRANGIAN';
  91. IFTOL  = 'NEG' WTAB.'FTOL' 'INCONNU' ;
  92. HYPDEF = 'MOT' WTAB.'HYPOTHESE_DEFORMATIONS';
  93. IGRD   = WTAB.'GRANDS_DEPLACEMENTS';
  94. IKSIA  = WTAB.'K_SIGMA';
  95. IKTAN  = WTAB.'K_TANGENT' ;
  96. IMPLP  = WTAB.'LIAISON_PERSISTANTE';
  97. IMTOL  = 'NEG' WTAB.'MTOL' 'INCONNU';
  98. IPILOT = WTAB.'AUTOMATIQUE';
  99. IPLAST = WTAB.'PLASTIQUE';
  100. IPLAVI = WTAB.'IPLAVI';
  101. IPREDIC = WTAB.'PREDICTEUR';
  102. IRCON  = WTAB.'RAIDCONST';
  103. IRAUGLU  = WTAB.'RAIDAUGM';
  104. IRAUG = IRAUGLU ;
  105. AUTAUG = WTAB.'AUTOAUGM';
  106. ISOL   = WTAB.'CONSOLIDATION';
  107. ISSTE  = WTAB.'SUBSTEPPING';
  108. ITHER  = WTAB.'CHAR_THE' 'OU' WTAB.'FOR_THER';
  109. IVIEXT = WTAB.'VISCO_EXTERNE';
  110. IVIDOM = WTAB.'VISCODOMMAGE';
  111. IVISCO = WTAB.'VISCOPLASTIQUE';
  112. LNLOC  = WTAB.'NLOC';
  113. LOGDEF = WTAB.'CHAR_DEFI';
  114. LOGPRE = WTAB.'CHAR_PRES' ;
  115. LOGPIL = WTAB.'CHAR_PILO' ;
  116. NITMA  = WTAB.'NITERINTER_MAX';
  117. NSSTE  = WTAB.'NMAXSUBSTEPS';
  118. POR1   = WTAB.'POR1' ;
  119. TI     = WTAB.'T_FINAL';
  120. EKREAC = WTAB.'REAC_GRANDS';
  121. ekreac = 1e50;
  122. ZMAXIT = WTAB.'MAXITERATION' ;
  123. ZNACCE = VRAI;
  124. XCONV = 0.;
  125. DEPSTDM = 0.;
  126. ZNCONS = WTAB.'NITER_KTANGENT'     ;
  127. ZPREC  = WTAB.'PRECISION'          ;
  128. ZPREK  = WTAB.'PRECISINTER'        ;
  129. ZPRECD = WTAB.'PRECDECHARGE'       ;
  130. ZPRECM = WTAB.'PRECFLEX'           ;
  131. ZCLIM0 = WTAB.'BLOCAGES_MECANIQUES';
  132. *-- Autres initialisations en non-local helmhoktz --
  133. 'SI' (LNLOC 'ET' ('EGA' WTAB.'NON_LOCAL' 'HELM'));
  134.   TAHELM = WTAB.'HELMHOLTZ' ;
  135.   NHELM =  TAHELM . 'N_VARI_NL' ;
  136. 'FINSI' ;
  137. *-- On initialise ZCLIM qui va contenir l'ensemble des C.L.
  138. ZCLIM  = ZCLIM0    ;
  139. DT     = WTAB.'DT' ;
  140. TEMPS0=WTAB.'TEMPS0';
  141. *---  doit on reactualiser la geometrie temp0 ne 0 et grand depl ? ---
  142. 'SI'  WTAB.'RECALCUL'  ;
  143.     WTAB . 'RECARI'   = VRAI ;
  144.     WTAB . 'RECADET'  = VRAI ;
  145.     WTAB . 'REA_GEOM' = VRAI ;
  146. *   on suppose que l'on est sur la bonne configuration
  147.     GEOM1 = 'FORM';
  148. *  'FORM'  GEOM1;
  149. 'SINON';
  150.     GEOM1 = WTAB.'FOR0' ;
  151.     WTAB . 'RECARI'  = FAUX ;
  152.     WTAB . 'RECADET' = FAUX ;
  153. 'FINS';
  154.  
  155. 'SI' ('OU' ('OU' IENDOM  IVIDOM) ICERAM);
  156.      WTAB.'RECARI'= VRAI;
  157.      LAG_TOT=0;
  158. 'FINS';
  159.  
  160. 'SI' IRAUG;
  161.    RIG_AUG  = WTAB.'RIGIDITE_AUGMENTEE';
  162. 'FINS';
  163.  
  164. 'SI' IRCON;
  165.    RIG_CONS = WTAB.'RIGIDITE_CONSTANTE';
  166.    MAI_CONS ='EXTR' RIG_CONS 'MAIL'    ;
  167. 'FINS';
  168.  
  169. *---------- chamelem etat estime pour la fin du pas de temps ------
  170.  ZETAT2 = PAS_ETAT PRECED TI    ;
  171.  ZMAT22 = PAS_MATE PRECED ZETAT2;
  172.  ZMATFI = ZMAT22                ;
  173.  
  174. *----------------------- nouveau chargement ----------------------
  175. 'SI' ('EXIS' (WTAB.'CHARGEMENT') 'MECA');
  176.       ZFEXT2 = 'TIRE' WTAB.'CHARGEMENT' 'MECA' TI;
  177.       TYP_2  = 'TYPE' ZFEXT2;
  178.       'SI' ('NEG' TYP_2 'CHPOINT ');
  179.  'MESS' '*** Erreur dans la definition du chargement ( MECA ) ***';
  180.  'ERRE' '*** Le type du champ n est pas CHPOINT ***';
  181.       'FINS';
  182. 'SINON';
  183.     ZFEXT2= 'VIDE' 'CHPOINT' / 'DISCRET';
  184. 'FINS';
  185.  
  186. 'SI' ('EXIS' (WTAB.'CHARGEMENT') 'FORC');
  187.       F2_FOR = 'TIRE' WTAB.'CHARGEMENT' 'FORC' TI;
  188.        TYP_2 = 'TYPE' F2_FOR;
  189.        'SI' ('NEG' TYP_2 'CHPOINT ');
  190.  'MESS' '*** Erreur dans la definition du chargement ( FORC ) ***';
  191.  'ERRE' '*** Le type du champ n est pas CHPOINT ***';
  192.       'FINS';
  193.        ZFEXT2 = ZFEXT2 '+' F2_FOR ;
  194. 'FINS';
  195.  
  196. 'SI' (LOGPRE 'ET' ('NON' IGRD)) ;
  197.    MOP   = WTAB.'MOD_PRE' ;
  198.    ZPEXT ='TIRE' WTAB.'CHARGEMENT' 'PRES' TI ;
  199.   'SI' ('EXIS' WTAB 'MAT_PRE') ;
  200.      ZFPEXT = 'BSIG' MOP ZPEXT ('REDU' ZMATFI MOP) ;
  201.   'SINON' ;
  202.      ZFPEXT = 'BSIG' MOP ZPEXT ;
  203.   'FINS' ;
  204.    ZFEXT2 = ZFEXT2 '+' ZFPEXT ;
  205. 'FINS' ;
  206.  
  207.  
  208. *--------------- Si il existe des deplacements imposes --------------
  209.          'SI' ('EXIS' (WTAB.'CHARGEMENT') 'DIMP');
  210.              F2_DEP = 'TIRE' WTAB.'CHARGEMENT' 'DIMP' TI;
  211.              ZFEXT2 = ZFEXT2 + F2_DEP;
  212.          'FINS';
  213. *--------------- Si il un increment de deplacements imposes -----------------
  214.          'SI' ('EXIS' (WTAB.'CHARGEMENT') 'DINC');
  215.              F2_INC = 'TIRE' WTAB.'CHARGEMENT' 'DINC' TI;
  216. *              list f2_inc;
  217.              F2_mail = extr f2_inc maillage;
  218.              zclim_mail = zclim extrai 'MAILLAGE';
  219.              F2_Mcmp = zclim_mail elem appu larg f2_mail;
  220.              zclim_inc = zclim redu f2_mcmp;
  221.  
  222.              f2_base = (zclim_inc * conti.'DEPLACEMENTS') exco 'FLX' 'FLX';
  223. *              list  f2_base;
  224.              ZFEXT2 = ZFEXT2 + F2_INC + F2_base; ;
  225.          'FINS';
  226. *---------------- si chargement deformation actualisation DEFOR ------
  227.          'SI' WTAB.'CHAR_DEFI' ;
  228.              WTAB.'DEFOR2' = 'TIRE' WTAB.'CHARGEMENT' 'DEFI' TI;
  229.          'FINS';
  230.  
  231. *---------- dynamique : preparation du second membre --------------
  232.          'SI' ( WTAB.'DYNAMIQUE' ) ;
  233.             'SI'  ('EGA' WTAB.'FREA1' 'INCONNU');
  234.                'SI' ('EXIS' WTAB.'CHARGEMENT' 'MECA');
  235.                    F1 = 'TIRE' WTAB.'CHARGEMENT' 'MECA' TEMP0;
  236.                'FINS';
  237.                'SI' ('EXIS' WTAB.'CHARGEMENT' 'FORC');
  238.                   F1F = 'TIRE' WTAB.'CHARGEMENT' 'FORC' TEMP0;
  239.                 'SI' ('EGA' ('TYPE' F1) 'CHPOINT ') ;
  240.                   F1 = F1 + F1F ;
  241.                 'SINON';
  242.                   F1 = F1F ;
  243.                 'FINS' ;
  244.                'FINS';
  245.  
  246.                'SI' WTAB.'PROCEDURE_CHARMECA';
  247.                   TFF1 = CHARMECA PRECED TEMPS0 ;
  248.                  'SI' ('EXIS' TFF1 'ADDI_SECOND');
  249.                     FF1=TFF1 .'ADDI_SECOND';
  250.                  'FINS';
  251.                  'SI' ('EGA' ('TYPE' F1) 'CHPOINT ');
  252.                     F1=F1 + FF1;
  253.                  'SINON';
  254.                     F1 = FF1;
  255.                  'FINS';
  256.                'FINS';
  257.  
  258.                'SI' LOGPRE ;
  259.                  MOP = WTAB.'MOD_PRE' ;
  260.                  ZPEXT0 = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TEMPS0 ;
  261.                  'SI' ('EXIS' WTAB 'MAT_PRE') ;
  262.                    FF1 = 'BSIG' MOP ZPEXT0 ('REDU' WTAB.'MAT1' MOP) ;
  263.                  'SINON' ;
  264.                    FF1 = 'BSIG' MOP ZPEXT0 ;
  265.                  'FINS' ;
  266.                  'SI' ('EGA' ('TYPE' F1) 'CHPOINT ');
  267.                    F1=F1 + FF1;
  268.                  'SINON';
  269.                    F1 = FF1;
  270.                  'FINS';
  271.                'FINS';
  272. *
  273.                  LAF0 = 'BSIG' WTAB.'MO_TOT' conti.'CONTRAINTES' ZMAT22;
  274.                'SI'  IRCON  ;
  275.                   LAF0 = LAF0 'ET'
  276.                   ('REDU' MAI_CONS (RIG_CONS '*' conti.'DEPLACEMENTS'));
  277.                'FINS';
  278.                'SI' (('EXIS' WTAB.'CHARGEMENT' 'MECA') 'OU'
  279.                     ('EXIS' WTAB.'CHARGEMENT' 'FORC') 'OU'
  280.                    wtab.'PROCEDURE_CHARMECA');
  281.                   LAF1  = F1  - LAF0 ;
  282.                'SINON';
  283.                   LAF1 = -1 * LAF0;
  284.                'FINS';
  285.                LAF2 = 'ENLEVER' (ZCLIM0 '*' conti.'DEPLACEMENTS') 'FLX';
  286. * forces exterieures + reactions - forces interieures au debut du calcul
  287. * c.a.d. (masse*acceleration initiale)+(amortissement*vitesse initiale)
  288.                 WTAB.'FREA1' = LAF1 - LAF2 ;
  289.             'FINS';
  290.             'SI' wtab.'LIAISON_PERSISTANTE' ;
  291. *  forces d'acceleration au debut du pas
  292.               'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  293.                  FF4 = WTAB.'AMORTISSEMENT'* conti.'VITESSES';
  294.                  WTAB.'FMAN'= WTAB.'FREA1' - FF4 ;
  295.               'SINON' ;
  296.                  WTAB.'FMAN'= WTAB.'FREA1' ;
  297.               'FINS';
  298.             'FINS';
  299. *------------  il faut calculer la matrice de masse tout de suite -------------
  300.             'SI' ( ('NON' ('EXIS' WTAB MASSE )) 'OU'
  301.                                WTAB.'GRANDS_DEPLACEMENTS');
  302.                 WTAB.'MASSE' = 'MASS' WTAB.'MO_TOT' ('CHAN' 'MASSE'
  303.                          ZMAT22 WTAB.'MO_TOT' ) ;
  304.                'SI' WTAB.'MASSCONST';
  305.                   WTAB.'MASSE'=WTAB.'MASSE' 'ET' WTAB.'MASSE_CONSTANTE';
  306.                'FINS';
  307.             'FINS';
  308.             FF = WTAB.'MASSE' *conti.'VITESSES';
  309.             FF4 = 4. / DT * FF; 'DETR' FF;
  310. * partie du second membre qui ne depend que des informations du pas prec
  311.             WTAB.'FREA1' = FF4 + WTAB.'FREA1';
  312.          'FINS';
  313. *--------- consolidation  : preparation du second membre ---------
  314.           'SI' WTAB.'CONSOLIDATION' ;
  315.             FF  = 'BSIG'  WTAB.'MOD_POR' conti.'CONTRAINTES' ;
  316.             FF4 = 'EXCO'  WTAB.'MOT_POR' FF
  317.                     WTAB.'MOT_POR' 'NOID' 'NATURE' 'DISCRET' ;
  318.             'SI' WTAB.'DYNAMIQUE';
  319.                 WTAB.'FREA1' = WTAB.'FREA1' +   FF4   ;
  320.             'SINON';
  321.                  WTAB.'FREA1' =  FF4   ;
  322.             'FINS';
  323.              ZFEXT2 = ZFEXT2 '+' FF4  ;  'DETR' FF ;
  324. *           ---- traitement des flux si besoin ----
  325.             'SI' ( 'EXIS' WTAB.'CHARGEMENT' 'FLUX' ) ;
  326.                FLUXT0= 'TIRE' WTAB.'CHARGEMENT' 'FLUX' TEMPS0;
  327.                FLUXTI= 'TIRE' WTAB.'CHARGEMENT' 'FLUX' TI  ;
  328.                FACFLU =  -1. * WTAB.'DT';
  329.                FLUXT =   ( FACFLU * (1 - WTAB.'TETA') * FLUXT0 )
  330.                        + ( FACFLU * WTAB.'TETA' *FLUXTI ) ;
  331.                ZFEXT2 = ZFEXT2 '+' FLUXT ;
  332.               'DETR' FLUXT ; 'DETR' FLUXT0; 'DETR' FLUXTI;
  333.            'FINS' ;
  334.         'FINS';
  335. *----------------- non-local type HELM : preparation --------------
  336.         'SI' (LNLOC 'ET' ('EGA' WTAB.'NON_LOCAL' 'HELM'));
  337.              PAS_HELM PRECED ;
  338.         'FINS' ;
  339. *----------------- calcul de la masse si frequentiel -------------
  340.         'SI' ( WTAB.'FREQUENTIEL' 'ET' ('NON' ('EXIS' WTAB 'MASSE' )));
  341.              WTAB.'MASSE' = 'MASS' WTAB.'MO_TOT'('CHAN' 'MASSE'
  342.                          ZMAT22 WTAB.'MO_TOT' ) ;
  343.         'FINS';
  344.  
  345. *----------------- pilotage indirect : preparation --------------
  346.         'SI' (LOGPIL);
  347.              ZFPILIN = WTAB.'FORCES_PILOTEES' 'ET' WTAB.'DEPLACEMENTS_PILOTES' ;
  348.                          DTAU = ('IPOL' WTAB.'PARAMETRE_DE_PILOTAGE' TI) -
  349.                                 ('IPOL' WTAB.'PARAMETRE_DE_PILOTAGE' TEMPS0) ;
  350.              ETA0 = 'EXTR' PRECED.'COEFFICIENT_DE_PILOTAGE' 
  351.                              ('DIME' PRECED.'COEFFICIENT_DE_PILOTAGE') ;
  352.              D_ETA=0.;                                                   
  353.         'FINS' ;
  354.  
  355. *-      Second membre
  356.         ZFCONSTA = ZFEXT2 ;
  357.  
  358.  
  359. ************************************************************************
  360. *             Parallélisation du GIBIANE via les ASSISTANTS            *
  361. ************************************************************************
  362. ZMODLI = WTAB.'MO_TOT'     ;
  363. ZMODLP = WTAB.'MO_TOT_PREC';
  364. NBPART = WTAB.'NBPART'     ;
  365.  
  366. PARALLEL = FAUX ;
  367. PARTLOCA = FAUX ;
  368. ZMODL  = ZMODLI ;
  369. 'SI' ('EGA' WTAB.'PROCESSEURS' 'COMPORTEMENT') ;
  370.   PARALLEL = VRAI ;
  371.   PARTLOCA = VRAI ;
  372.   MODRELOC = 'PART' 'ARLE' ZMODLI NBPART ;
  373. 'FINS' ;
  374. 'SI' ('EGA' WTAB.'PROCESSEURS' 'AUTOMATIQUE') ;
  375.   PARALLEL = VRAI ;
  376.   PARTLOCA = FAUX ;
  377.   ZMODL = 'PART' 'ARLE' ZMODLI NBPART ;
  378.   'OPTI' 'PARA' VRAI ;
  379. 'FINS' ;
  380.  
  381. * Test sur un MODELE qui aurait changé
  382. 'SI' ('NEG' ZMODLI ZMODLP)                                   ;
  383.    CONTI.'CONTRAINTES' =('REDU' CONTI. 'CONTRAINTES' ZMODLI) '+' ('ZERO' ZMODLI 'CONTRAINTES') ;
  384.    ESTIM.'CONTRAINTES' =('REDU' ESTIM. 'CONTRAINTES' ZMODLI) '+' ('ZERO' ZMODLI 'CONTRAINTES') ;
  385.    CONTI.'DEFORMATIONS'=('REDU' CONTI.'DEFORMATIONS' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT'   ) ;
  386.    ESTIM.'DEFORMATIONS'=('REDU' ESTIM.'DEFORMATIONS' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT'   ) ;
  387.   'SI' IPLAVI ;
  388.      CONTI.'DEFORMATIONS_INELASTIQUES'=('REDU' CONTI.'DEFORMATIONS_INELASTIQUES' ZMODLI) '+' ('ZERO' ZMODLI 'DEFINELA') ;
  389.      ESTIM.'DEFORMATIONS_INELASTIQUES'=('REDU' ESTIM.'DEFORMATIONS_INELASTIQUES' ZMODLI) '+' ('ZERO' ZMODLI 'DEFINELA') ;
  390.      CONTI.'VARIABLES_INTERNES'       =('REDU' CONTI.'VARIABLES_INTERNES'        ZMODLI) '+' ('ZERO' ZMODLI 'VARINTER') ;
  391.      ESTIM.'VARIABLES_INTERNES'       =('REDU' ESTIM.'VARIABLES_INTERNES'        ZMODLI) '+' ('ZERO' ZMODLI 'VARINTER') ;
  392.   'FINS';
  393.  
  394.    WTAB.'ETAT1' = PAS_ETAT PRECED TEMPS0      ;
  395.    WTAB.'MAT1'  = PAS_MATE PRECED WTAB.'ETAT1';
  396.  
  397.   'SI' ITHER ;
  398.      WTAB.'ETHER1'    =('REDU' WTAB.'ETHER1'   ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT') ;
  399.     'SI' POR1;
  400.        WTAB.'MSRTHER1'=('REDU' WTAB.'MSRTHER1' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT') ;
  401.     'FINS';
  402.   'FINS';
  403.  
  404.   'SI' ('EXIS' WTAB 'DEFOR1' );
  405.      WTAB.'DEFOR1' =
  406.                    ('REDU'  WTAB.'DEFOR1' ZMODLI) '+'
  407.                    ('ZERO' ZMODLI 'DEFORMAT'    ) ;
  408.   'FINS';
  409.  
  410.   'SI' ('EXIS' WTAB 'DEFOR2' );
  411.      WTAB.'DEFOR2' =
  412.                    ('REDU'  WTAB.'DEFOR2' ZMODLI) '+'
  413.                    ('ZERO' ZMODLI 'DEFORMAT'    ) ;
  414.   'FINS';
  415. 'FINS';
  416.  
  417. *-----------materiau au debut du pas ------------------------------
  418. ZMAT11 = WTAB.'MAT1'        ;
  419. 'SI' ('EGA' WTAB.'PROCESSEURS' 'AUTOMATIQUE') ;
  420.    ZMAT1 ='REDU' ZMAT11 ZMODL ; COMM 'REDU au cas ou le modele est parallele';
  421. 'SINO';
  422.    ZMAT1 = ZMAT11;
  423. 'FINS';
  424.  ZMAT1I  = ZMAT1; COMM 'Version Initiale de ZMAT1 pour BONOCONV';
  425.  
  426.  
  427. ************************************************************************
  428. *                       Quelques initialisations                       *
  429. ************************************************************************
  430. STAB12             = 'TABL'                    ;
  431. STAB12.'ZU1'       =  CONTI.'DEPLACEMENTS'     ;
  432. STAB12.'SIGF'      =  CONTI.'CONTRAINTES'      ;
  433. STAB12.'DEFF'      =  CONTI.'DEFORMATIONS'     ;
  434. STAB12.'FNONL'     =  WTAB.'FNONL'             ;
  435. STAB12.'RESIDU'    =  WTAB.'RESIDU'            ;
  436. STAB12.'XDENO'     =  WTAB.'XDENO'             ;
  437. STAB12.'XDENOM'    =  WTAB.'XDENOM'            ;
  438. 'SI' ( 'EXIS' WTAB 'ETAT1' );
  439.   STAB12.'ETAT1'   = WTAB.'ETAT1';
  440. 'FINS';
  441. STAB12.'ETAT2'= ZETAT2;
  442. 'SI' ('EXIS' WTAB 'DEFOR1' );
  443.   STAB12.'DEFOR1'  = WTAB.'DEFOR1'  ;
  444.   STAB12.'DEFOR2'  = WTAB.'DEFOR2'  ;
  445. 'FINS';
  446. 'SI' ( 'EXIS' WTAB 'FNONL');
  447.   STAB12.'FNONL'   = WTAB.'FNONL'   ;
  448. 'FINS';
  449. 'SI' ('EXIS' WTAB  'TET1') ;
  450.   STAB12.'TET1'    = WTAB.'TET1'    ;
  451.   STAB12.'TET2'    = WTAB.'TET2'    ;
  452. 'FINS';
  453. * STAB12.'SUCCES' = VRAI ;
  454. 'SI' ('NEG' WTAB.'AUTOCOEF' 'INCONNU') ;
  455.   STAB12.'AUTOCOEF' = WTAB.'AUTOCOEF' ;
  456. 'FINS' ;
  457. 'SI' ('NEG' WTAB.'AUTOREDU' 'INCONNU') ;
  458.   STAB12.'AUTOREDU' = WTAB.'AUTOREDU' ;
  459. 'FINS' ;
  460. 'SI' ('NEG' WTAB.'SECOND_MEMBRE' 'INCONNU') ;
  461.   STAB12.'SECOND_MEMBRE' = WTAB.'SECOND_MEMBRE' ;
  462. 'FINS' ;
  463. 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  464.   STAB12.'LASTKTAN' = WTAB.'LASTKTAN' ;
  465. 'FINS' ;
  466. 'SI' ('NEG' WTAB.'AUTORED1' 'INCONNU') ;
  467.   STAB12.'AUTORED1' = WTAB.'AUTORED1' ;
  468. 'FINS' ;
  469. 'SI' ('NEG' WTAB.'LISEA_M' 'INCONNU') ;
  470.   STAB12.'LISEA_M' = WTAB.'LISEA_M' ;
  471.   STAB12.'RIBLO_M' = WTAB.'RIBLO_M' ;
  472. 'FINS' ;
  473. 'SI' ('NEG' WTAB.'INCREMENT' 'INCONNU' );
  474.   STAB12.'INCREMENT' = WTAB.'INCREMENT';
  475.   INCRPREC = STAB12.'INCREMENT'        ;
  476. 'FINS';
  477. STAB12.'FFROT'    = WTAB.'FFROT'    ;
  478. STAB12.'INITEMPS' = WTAB.'INITEMPS' ;
  479. STAB12.'DT'       = WTAB.'DT'       ;
  480. STAB12.'DTPREC'   = WTAB.'DTPREC'   ;
  481.  
  482. 'SI'  ITHER  ;
  483.    STAB12.'TETA1'  = WTAB.'TET1';
  484.    STAB12.'TETA2'  = WTAB.'TET2';
  485.    ETT0 ='REDU' WTAB.'ETHER1' ZMODL;     'COMM' 'REDU au cas ou le modele soit parallele';
  486.   'SI' POR1;
  487.     MSRTT0='REDU' WTAB.'MSRTHER1' ZMODL; 'COMM' 'REDU au cas ou le modele soit parallele';
  488.   'FINS';
  489.  
  490. 'SINON' ;
  491.    che3 ='MANU' 'CHML' ZMODL 'T' 20. 'STRESSES' ;
  492.    che4 ='MANU' 'CHML' ZMODL 'T' 20. 'STRESSES' ;
  493. 'FINS' ;
  494.  
  495. ZMATI   ='REDU' ZMAT22 ZMODLI ;
  496. ZMAT    ='REDU' ZMATI  ZMODL  ;
  497.  
  498. * SP : initialisation DFGRAD en presence d'un modele MECANIQUE
  499. *      (DFGRAD mis a INCONNU par defaut dans PAS_INIT)
  500. 'SI' ('NEG' WTAB.'MOD_MEC' 'INCONNU') ;
  501.   'SI' (('NEG' WTAB.'DFGRAD' 'INCONNU') 'ET' IGRD) ;
  502.     geot = 'FORM' ; 'FORM' WTAB.'FOR0' ;
  503.     STAB12.'DFGRAD' ='GRAD' ZMODL ZMAT STAB12.'ZU1' ;
  504.     'FORM' geot ;
  505. **'SINO';
  506. **  STAB12.'DFGRAD' ='REDU' WTAB.'DFGRAD' ZMODL;
  507.   'FINS' ;
  508. 'SINON' ;
  509.   STAB12.'DFGRAD' = WTAB.'DFGRAD' ;
  510. 'FINS' ;
  511.  
  512. com_sig  = 'EXTR' ZMODLI  'CONT';
  513. HPP_EPS  =  FAUX                ;
  514. EPS_EPS  = 'TEXT' '        '    ;
  515. * Option a n'utiliser que par les utilisateurs avertis
  516. 'SI' ('EXIS' PRECED 'ACCELERATION') ;
  517.   III = PRECED.'ACCELERATION' ;
  518.   'SI' ('EGA' ('TYPE' III) 'ENTIER') ; ZNACCE = III ; 'FINS' ;
  519. 'FINS' ;
  520. * Matrice tangente : non utilisee si IPLAVI a FAUX
  521. IKTAN  = IKTAN 'ET' IPLAVI ;
  522. 'SI' (WTAB.'K_TANGENT' 'ET' ('NON' IPLAVI)) ;
  523.   'MESS' 'IPLAVI faux : pas de matrice tangente ->'
  524.          ' on utilise la rigidite elastique' ;
  525. 'FINS' ;
  526. * Matrice tangente par perturbation :
  527. * Option non disponible si non local ou si IPLAVI a FAUX
  528. IPERT = WTAB.'K_TANGENT_PERT' 'ET' ('NON' LNLOC) 'ET' IPLAVI ;
  529. ZPERC1 = WTAB.'K_TANG_PERT_C1' ;  ZPERC2 = WTAB.'K_TANG_PERT_C2' ;
  530. * Matrice tangente : partie symetrique utilisee
  531. 'SI' WTAB.'K_TANGENT_SYME' ;
  532.   ZKTASYM = 'MOT' 'SYME' ;
  533. 'SINON' ;
  534.   ZKTASYM = 'TEXTE' '    ' ;
  535. 'FINS';
  536. * Matrice tangente : pas d'acceleration en cas de modele FEFP ou SSTE
  537. 'SI' IKTAN ;
  538.   'SI' (IFEFP 'OU' ISSTE) ; ZNACCE = FAUX; 'FINS';
  539. 'FINS' ;
  540. *
  541. 'SI' IFTOL ;
  542.     ZFTOL = 'ABS' WTAB.'FTOL' ;
  543. 'FINS';
  544. 'SI'  IMTOL  ;
  545.     ZMTOL = 'ABS' WTAB.'MTOL' ;
  546. 'FINS';
  547. ITCAR = ( 'EXIS' ZMATI  'EPAI') 'OU' ('EXIS' ZMATI 'INRY') 'OU'
  548.         ( 'EXIS' ZMATI  'MODS') 'OU'
  549.        (( 'EXIS' ZMODLI 'ELEM' 'JOI1') 'ET' ('EXIS' ZMATI 'V1X ')
  550.                                        'ET' ('EXIS' ZMATI 'V1Y '))  ;
  551. *
  552. * CB215821 : Devrait t-on mutualiser ITCAR et WTAB.ITCAR? ==>
  553. *            Ce ne sont pas tout a fait les memes tests aujourd'hui
  554. 'SI' (ITCAR 'EGA' FAUX);
  555.    ZMAT2 = ZMATI ;
  556.    ZMAT2R= ZMAT  ;
  557. 'FINS';
  558.  
  559. 'SI' WTAB.'ITCAR';
  560.   'SI' ('EGA' WTAB.'CARA' ('MOT' 'INCONNU'));
  561.      WTAB.'CARA' = ZMAT11 ;
  562.   'FINS';
  563.    CARA1 = ZMAT1;
  564. 'FINS';
  565.  
  566. 'SI' ('OU' ('OU' IVISCO IVIDOM) IVIEXT);  ZPREK = 5.E-7 ;  'FINS';
  567. 'SI' IENDOM;                              ZPREK = ZPREC ;  'FINS';
  568.  
  569. * on fait ici la séparation poreux .. pour l'avoir sur
  570. *les modèles partitionnes
  571. 'SI' POR1;
  572.     MO_PORI = 'EXTR' ZMODLI 'FORM' 'POREUX';
  573.     MO_POR  = 'EXTR' ZMODL  'FORM' 'POREUX';
  574.  
  575.     MA_POR  = 'REDU' ZMAT22 MO_POR      ;
  576. ** kich ma_por0 intervient si ISOL
  577. ** initialement     MA_POR0= 'REDU' MO_POR (STAB12.'MAT1');
  578.     MA_POR0 = 'REDU' ZMAT11 MO_POR      ;
  579.  
  580.     MAI_POR = 'EXTR' MO_POR  'MAILLAGE' ;
  581.     MAI_PORI= 'EXTR' MO_PORI 'MAILLAGE' ;
  582. 'FINS';
  583.  
  584. * recuperation de certains champs, si nbpart>1 zmodl est partitionné
  585. * sinon c'est le modele initial
  586. *                                                                      *
  587. DEFT0        = 'REDU' conti.'DEFORMATIONS' ZMODL              ;
  588. ZSIGF        = 'REDU' STAB12.'SIGF' ZMODL                     ;
  589. 'SI' IPLAVI ;
  590.      ZDEIF   = 'REDU' CONTI.'DEFORMATIONS_INELASTIQUES' ZMODL ;
  591.      ZVARF   = 'REDU' CONTI.'VARIABLES_INTERNES'        ZMODL ;
  592.      com_var = 'EXTR' ZMODLI 'VARI' ;
  593.      com_dei = 'EXTR' ZMODLI 'DEIN' ;
  594.      lnom    =  com_var             ;
  595.     'SI' ISOL  ;
  596.        com_maa = 'EXTR' ZMODLI 'MATE' ;
  597.     'FINS' ;
  598. 'FINS' ;
  599. *
  600. *   teste t'on les moments  ?
  601. TSTMOM = ITCAR 'OU' ('EGA' ('VALE' 'MODE') 'PLANGENE') ;
  602. *   teste t'on les POREUX   ?
  603. 'SI' POR1  ; TSTMOM=VRAI ; 'FINS';
  604. *
  605. IKLFFF=VRAI;
  606. 'SI'TSTMOM;    'SI' IFTOL;     'SI' IMTOL;
  607.        IKLFFF=FAUX;
  608. 'FINS';        'FINS';        'FINS';
  609. *
  610. 'SI'('NON' TSTMOM);
  611.     'SI' IFTOL;
  612.        IKLFFF=FAUX;
  613.     'FINS';
  614. 'FINS';
  615. *
  616. GEOREF0  = WTAB.'FOR0' ;
  617. WTAB.'CONV'=VRAI;
  618. ISOUSPPP=0;
  619. 'SI' ('EXIS' wtab 'ISOUSPAS'); ISOUSPPP = wtab.'ISOUSPAS'; 'SINON'; isousppp=0; 'FINSI';
  620. WTAB . 'ISOUSPAS' = 0;
  621. NSOUSPAS = WTAB . 'MAXSOUSPAS';
  622. ZCCONV = VRAI ;
  623. KNOCONV = 0 ;
  624.  
  625. * Initialisation CHArgement SANS T :
  626. CHASANST= WTAB.'CHARGEMENT' 'ENLE' 'T' ;
  627.  
  628. RED_URG = 0;
  629.    residc = 'VIDE' 'CHPOINT' / 'DISCRET';
  630. 'SI' autaug; iraug = faux; 'FINSI';
  631. ************************************************************************
  632. ****** boucle de non convergence
  633. ************************************************************************
  634. * isoucomp : Indicateur SOUci dans COMPortement, FAUX par defaut 
  635. isoucomp = faux ;
  636. resmul = 1;
  637.     augmult = 0.60000000;
  638.      RFNS  = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  639.      RTRSF = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  640.      NONCONV = FAUX;
  641.      excfconv = 0;
  642. 'REPETER'  BONOCONV NSOUSPAS ;
  643. 'SI' ('EGA' RED_URG 0);
  644.  'SI' ((XCONV < ZPREC) 'ET' (DEPSTDM < ZPREC) 'ET' (resmul > 0.99) 'ET' (AUGMULT < 100.));
  645.     augmult = 0.60000000;
  646. *pv resmul = 1;
  647.     augauto = augmult;
  648.     augk    = augmult;
  649.     'SI' autaug ; iraug = faux; 'FINSI';
  650.     znacce=VRAI;
  651.   'SI' (IFEFP 'OU' ISSTE) ; ZNACCE = FAUX; 'FINS';
  652.  'SINON';
  653.     augmult = augmult / 1.5;
  654. **  si (augmult < 0.6); augmult = 0.6; finsi;
  655.  'FINSI';
  656. 'SINON';
  657.    augmult = augmult * 1.5;
  658.  'SI' autaug;   IRAUG = VRAI; 'FINSI';
  659. 'FINSI';
  660. ***   wtab.'RECARI' = vrai;
  661.  
  662.   KNOCONV       = KNOCONV+1                 ;
  663.   STAB12.'CONV' = FAUX                      ;
  664.   DT_INIT       = STAB12.'DT'               ;
  665.   DTINI         = STAB12.'DT'               ;
  666.  
  667.   ZSIG0         = ZSIGF                     ;
  668.   zsig0_1       = zsig0;
  669.  'SI' IPLAVI ;
  670.     ZEPS0 =  ZDEIF                          ;
  671.     ZVAR0 =  ZVARF                          ;
  672.  'FINS';
  673.  
  674.   ZU1           = STAB12.'ZU1'   ;
  675.   ZDETOT = 'ENLE' ZU1 'LX';
  676.  GR_U_DEB = 'MOT' 'INCONNU';
  677.  'SI' ('EXIS' STAB12 'DFGRAD' );
  678.   GR_U_DEB      = STAB12.'DFGRAD';
  679.  'FINSI';
  680.   'SI'  ITHER  ;
  681.      TETA1 = STAB12.'TETA1' ;
  682.      TETA2 = STAB12.'TETA2' ;
  683.      DTETD = TETA2 '-' TETA1;
  684.   'FINS' ;
  685.  
  686. * materiau au debut du pas en cas de non convergence etat1=etat2
  687.   'SI' (knoconv > 1) ;
  688.      ZMAT1 = 'REDU' (PAS_MATE PRECED STAB12.'ETAT1') ZMODL ;
  689.   'SINON';
  690.      ZMAT1 = ZMAT1I ;
  691.   'FINS' ;
  692. *------ caracteristiques initiales en cas de grands deplacements ------
  693.   'SI' WTAB.'ITCAR';
  694.      MECAR1 = 'EXCO'  MOCA  ZMAT1  'NOID'            ;
  695.      MECAR2 = 'EXCO'  MOCA  CARA1  'NOID'            ;
  696.      MECAR1 = 'CHAN' 'TYPE' MECAR1 'CARACTERISTIQUES';
  697.      MECAR2 = 'CHAN' 'TYPE' MECAR2 'CARACTERISTIQUES';
  698.      ZMAT1  =  ZMAT1 '-' MECAR1 '+' MECAR2           ;
  699.   'FINS';
  700.  
  701. *----------- Calcul du champ de materiau a la fin du pas ------------
  702.   'SI' (knoconv '>' 1) ;
  703.      MCHC = PAS_ETAT PRECED TI  ;
  704.      MMMM = PAS_MATE PRECED MCHC;
  705.   'SINON' ;
  706.      MMMM = ZMAT22              ;
  707.   'FINS' ;
  708.    MMMM = 'REDU' MMMM ZMODL;
  709.  
  710. *----- Caracteristiques initiales en cas de grands deplacements -----
  711.   'SI' WTAB.'ITCAR';
  712.      MECAR1 = 'EXCO'  MOCA  MMMM   'NOID'            ;
  713.      MECAR2 = 'EXCO'  MOCA  CARA1  'NOID'            ;
  714.      MECAR1 = 'CHAN' 'TYPE' MECAR1 'CARACTERISTIQUES';
  715.      MECAR2 = 'CHAN' 'TYPE' MECAR2 'CARACTERISTIQUES';
  716.      MMMM = MMMM '-' MECAR1 '+' MECAR2               ;
  717. *   (fdp) on reporte les variations du materiau sur le pas
  718. *    dans les autres instances du champ materiau
  719. *   (ZMAT et ZMATI semblent suffirent)
  720.      MECAR1 = 'EXCO'  MOCA  ZMAT  'NOID'             ;
  721.      MECAR2 = 'EXCO'  MOCA  MMMM  'NOID'             ;
  722.      MECAR1 = 'CHAN' 'TYPE' MECAR1 'CARACTERISTIQUES';
  723.      MECAR2 = 'CHAN' 'TYPE' MECAR2 'CARACTERISTIQUES';
  724.      ZMAT   =  ZMAT  '-' MECAR1 '+' MECAR2           ;
  725.      ZMATI  = 'REDU' ZMAT ZMODLI                     ;
  726.   'FINS';
  727.  
  728. *------------ Calcul de la rigidite a la fin du pas ----------------
  729. *  Test ci-dessous est inutile, WTAB.'BLOCAGES_MECANIQUES' et ZCLIM0
  730. *  correspondent aux memes rigidites ... (cf. ligne 117)
  731.    AA1 ='EXTR' WTAB.'BLOCAGES_MECANIQUES' 'MAIL'     ;
  732.    AA2 ='EXTR' ZCLIM0 'MAIL' ;
  733.    AA3 ='DIFF' AA1 AA2       ;
  734.    AA4 = NBNO  AA3           ;
  735.   'SI' ( (WTAB.'RECARI' 'OU' IRAUG      ) 'OU'
  736.          ('NON' ('EXIS' WTAB 'RRRR')    ) 'OU' ('NEG' AA4 0));
  737.     'SI' ('OU' ('OU' IENDOM  IVIDOM) ICERAM);
  738.        HOOKENDO ='HOOK' ZMODL MMMM     ZVAR0;
  739.        RH       ='RIGI' ZMODL HOOKENDO MMMM ;
  740.       'DETR' HOOKENDO;
  741.     'SINON';
  742.        'SI' ((LAG_TOT 'EGA' 1) 'ET' vrai);
  743.           'FORM' GEOREF0;
  744.           HOOKRH = 'HOOK' ZMODL MMMM;
  745.           HOOKRH2 = 'PICA' HYPDEF HOOKRH ZMODL ZU1;
  746.           'FORM' GEOM1;
  747.           RH = 'RIGI' HOOKRH2 ZMODL MMMM;
  748.           'DETRUI' HOOKRH;
  749.           'DETRUI' HOOKRH2;
  750.        'SINON';
  751.           RH  = 'RIGI' ZMODL MMMM ;
  752.        'FINSI';
  753.     'FINS';
  754. *   RH peut contenir des CL
  755.     ZCL = 'EXTR' RH 'RIGI' 'MULT' ;
  756.     'SI' ('NEG' ('DIME' ZCL) 0) ;
  757.       ZCLIM = ZCLIM0 'ET' ZCL ;
  758.     'FINSI' ;
  759.     RRRR = RH 'ET' ZCLIM0 ;
  760. *   Prise en compte d'eventuelles RIGIDITE_CONSTANTE
  761.     'SI' IRCON;
  762.       RRRR = RRRR 'ET' RIG_CONS;
  763.     'FINS';
  764. *
  765.     ZRAIDNA = RRRR;
  766.     'SI' (IRAUG 'ET' ('NON' AUTAUG));
  767.       RRRR = RRRR 'ET' RIG_AUG  ;
  768.     'FINSI';
  769. *
  770.     'SI' (AUTAUG 'ET' ('NON' IRAUGLU));
  771.        RIG_AUG  = 'MASSE' ZMODL MMMM ;
  772. ***    mess 'actualisation rig_aug';
  773.     'FINSI';
  774.     'SI' (AUTAUG 'ET' IRAUG);
  775. *  recalcul masse pour l'augmentation
  776.       RRRR = RRRR 'ET' (RIG_AUG * augauto) 'ET' (RH   * augk) ;
  777.     'FINS';
  778.     'SI' ('NON' IRAUG);
  779. *  recalcul masse pour l'augmentation
  780.       'SI' IRAUGLU;
  781. ***     'MESS' ' augmentation residuelle ';
  782.         RRRR = RRRR 'ET' RIG_AUG ;
  783.     'FINS';
  784.     'FINSI';
  785.  
  786.     'SI' ('EGA' ('DIME' ZCL) 0) ;
  787. *   Stockage de la rigidite pour eviter de la recalculer
  788.     WTAB.'RRRR'=RRRR;
  789.     'FINSI';
  790. *
  791.   'SINON';
  792.      RRRR=WTAB.'RRRR';
  793.   'FINS';
  794.   ZRAID=RRRR;
  795.  
  796. *------------ consolidation ou dynamique faut-il recalculer l'operateur?
  797. * -----: preparation du pas de temps ------------
  798.     'SI' ( WTAB.'CONSOLIDATION' 'OU' WTAB.'DYNAMIQUE');
  799.        DT = TI '-' TEMP0;
  800.       'SI' ( '>' (DELTAN '*' 0.9999) DT) ;
  801.          WTAB . 'RECAOP' = VRAI;
  802.       'FINS';
  803.       'SI' ( '<' (DELTAN '*' 1.0001) DT) ;
  804.          WTAB . 'RECAOP' = VRAI;
  805.       'FINS';
  806.       'SI' WTAB.'MATVAR';
  807.          WTAB . 'RECAOP' = VRAI;
  808.       'FINS';
  809.        DELTAN=WTAB.'DT';
  810.     'FINS';
  811. *---------------------- Formation de l operateur -----------------------
  812.      'SI' ( WTAB.'DYNAMIQUE' 'OU' WTAB.'CONSOLIDATION');
  813.        'SI' ('NEG' WTAB.'OPERATEUR' 'INCONNU');
  814.           ZRAID = WTAB.'OPERATEUR';
  815.          'SI' (WTAB . 'RECAOP') ;
  816.             ZRAID = RRRR ;
  817.          'FINS';
  818.        'FINS';
  819.      'FINS';
  820.  
  821. *------------ operateur frequentiel -----------------------
  822. *------------ operateur amortissement en frequentiel
  823.      'SI' WTAB.'FREQUENTIEL' ;
  824.         RRR2 = 'AMOR'  ZMODL ZMAT ;
  825.         RR2  = 'CHAN' 'INCO' RRR2
  826.               ('MOTS' 'ALFA' 'BETA') ('MOTS' 'IALF' 'IBET')
  827.               ('MOTS' 'FALF' 'FBET') ('MOTS' 'FALF' 'FBET') 'QUEL' ;
  828.         RR3  = 'CHAN' 'INCO' RRR2
  829.               ('MOTS' 'ALFA' 'BETA') ('MOTS' 'ALFA' 'BETA')
  830.               ('MOTS' 'FALF' 'FBET') ('MOTS' 'IFAL' 'IFBE') 'QUEL' ;
  831.         RRR2 = RR2 'ET' RR3 ;
  832.  
  833.         RR1  = ZRAID ;
  834.         OMEGI= 2.* PI * TI ;
  835.         RRR1 = OMEGI * OMEGI * (-1.) * WTAB.'MASSE' ;
  836.         RR1  = ZRAID 'ET' RRR1 ;
  837.         RR4  = 'CHAN' 'INCO' (RR1 '*' (-1.D0))
  838.               ('MOTS' 'ALFA' 'BETA') ('MOTS' 'IALF' 'IBET')
  839.               ('MOTS' 'FALF' 'FBET') ('MOTS' 'IFAL' 'IFBE') 'QUEL';
  840.         ZRAID= RR1   'ET' RR4  ;
  841.         RR5  = OMEGI '*'  RRR2 ;
  842.         ZRAID= ZRAID 'ET' RR5  ;
  843.      'FINS' ;
  844.  
  845. *--------------- et la perméabilité  ----------------------------------
  846.      'SI' (WTAB.'CONSOLIDATION') ;
  847.        'SI' (WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'MATVAR') ;
  848.           WTAB .'PERMEABILITE'= 'PERM' WTAB.'MOD_POR' MMMM ;
  849.           WTAB .'RECAOP'      =  VRAI                      ;
  850.        'FINS';
  851.      'FINS';
  852.  
  853. *------------- Cas de la consolidation ou de la dynamique -------------
  854. *------------- il faut recalculer l'operateur d'iteration -------------
  855.      'SI' (WTAB . 'RECAOP') ;
  856.        'SI' ( WTAB.'DYNAMIQUE');
  857.           ZRAID =  4.D0 '/'( DT ** 2) '*' WTAB.'MASSE'     'ET' ZRAID;
  858.          'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  859.             ZRAID = WTAB.'AMORTISSEMENT' '*' (2.D0 '/' DT) 'ET' ZRAID;
  860.          'FINS';
  861.        'FINS' ;
  862.  
  863.        'SI' ( WTAB.'CONSOLIDATION');
  864.           ZRAID =-1.* DT* WTAB.'TETA'* WTAB.'PERMEABILITE' 'ET' ZRAID ;
  865.        'FINS' ;
  866.         WTAB . 'OPERATEUR'= ZRAID ;
  867.         WTAB . 'RECAOP'   = FAUX  ;
  868.      'FINS';
  869.      ZRAIDNA = ZRAID;
  870.  
  871.  
  872. *-------------- traitement des contacts frottements automatiques -------
  873.       CDEP    = STAB12.'ZU1'               ;
  874.       CDEPSLX ='ENLE' CDEP 'LX'            ;
  875.  
  876.       BFCONT  = FAUX                       ;
  877.       BCLIM2  = FAUX                       ;
  878.      'SI' WTAB.'CONTACT';
  879.         MODCON= WTAB.'MODCONTA';
  880.        'SI' ('EXIS' WTAB 'MATCONTA') ;
  881.           WTAB.'MATCONTA'= 'REDU' MODCON WTAB.'MATCONTA';
  882. **       list resu wtab.'MATCONTA';
  883.           MODFRO matfro  CJEU CRR RFROT='RFCO' MODCON  WTAB.'CONV' WTAB.'MATCONTA' ;
  884.        'SINON' ;
  885.           CJEU CRR RFROT='RFCO' MODCON  WTAB.'CONV'                ;
  886.        'FINS';
  887.  
  888. *
  889.        'SI' ('NEG' 0 CJEU) ;
  890.          FADHE = 'EXCO' CJEU 'FADH' 'NOID' 'FADH' ;
  891.          CDAP  = 'EXCO' CJEU 'FLX'  'NOID' 'FLX' ;
  892.        'FINSI' ;
  893.  
  894.        'SI' (WTAB . 'MODAL') ;
  895.           CRR  = 'PJBA' CRR   ZMODL  MMMM  ;
  896.           MCRR = 'EXTR' CRR  'MAIL' 'MULT' ;
  897.           MCDAP= 'EXTR' CDAP 'MAIL'        ;
  898.          'REPETER' BCDA ('NBNO' MCDAP) ;
  899.             PBCDA ='POIN' MCDAP &BCDA ;
  900.             PCRR  ='POIN' MCRR  &BCDA ;
  901.             CHCR  ='MANU' 'CHPO' PCRR 1 'FLX' ('EXTR' CDAP 'FLX' PBCDA)
  902.                    'NATURE' 'DISCRETE' ;
  903.            'SI' ('EGA' 1 &BCDA) ;
  904.               CCDA = CHCR ;
  905.            'SINON' ;
  906.               CCDA = CHCR 'ET' CCDA ;
  907.            'FINS' ;
  908.          'FIN' BCDA ;
  909.           CDAP = CCDA ;
  910.        'FINS' ;
  911.  
  912. *       ATTENTION : mettre les conditions de frottement en premier pour
  913. *                   les numeroter en dernier
  914.         ZCLIM2 = 'VIDE' 'RIGIDITE' / 'RIGIDITE';
  915.         ZFCONT  ='VIDE' 'CHPOINT' / 'DISCRET';
  916.        'SI'  ('NEG' CRR 0);
  917.           BCLIM2 = VRAI                ;
  918.           BFCONT = VRAI                ;
  919.           ZCLIM2 = CRR    'ET' ZCLIM2  ;
  920.           CCOR   = CRR    '*'  CDEPSLX 'EXCO' 'FLX' 'FLX';
  921.           ZFCONT = ZFCONT '+'  CCOR    ;
  922.        'FINS';
  923.  
  924.        'SI' ( 'NEG' 0 CJEU);
  925.           BFCONT = VRAI                ;
  926.           ZFCONT = ZFCONT '+'  CDAP    ;
  927.        'FINS';
  928.  
  929.        'SI' ('NEG' 0 RFROT);
  930.           BCLIM2 = VRAI                ;
  931.           BFCONT = VRAI                ;
  932.           ZCLIM2 = RFROT  'ET' ZCLIM2  ;
  933.           CCOR   = RFROT  '*'  CDEPSLX 'EXCO' 'FLX' 'FLX' ;
  934.           ZFCONT = ZFCONT '+'  CCOR    ;
  935.        'FINS';
  936. *
  937. *      Mise a jour de ZCLIM et ZFCONSTA si necessaire
  938.        'SI' BCLIM2;
  939.          ZCLIM  = ZCLIM2 'ET' ZCLIM ;
  940.          ZRAID  = ZCLIM2 'ET' ZRAID;
  941.        'FINS';
  942. *
  943.        'SI' BFCONT;
  944.          ZFCONSTA = ZFEXT2 '+' ZFCONT ;
  945.        'FINS';
  946.      'FINS';
  947. *
  948.       ZFCONST1 = ZFCONSTA;
  949.       zclimp = zclim;
  950.  
  951.  
  952. * ---------------  pilotage automatique ********************************
  953.      ISNPB = FAUX ;
  954.      AL1 = 1. ;COEPI = 1.d0; COEINC=0.d0;COEPI0=1.d0;DAL1=100.D0;
  955.      CORPREC = 1. ;
  956. *    CORPREC = 10.;
  957.      'SI' (('EGA' ipredic  'HPP') 'ET' WTAB.'CONV');
  958.           EPS_EPS = 'TEXTE' 'LINEAIRE';
  959.           HPP_EPS = VRAI;
  960.      'FINS';
  961. *
  962.      'SI' IPILOT ;
  963.           'SI' ( WTAB.'AUTODEUX' ) ;
  964.                 COEPI = 'ABS' ( STAB12.'AUTOCOEF');
  965.                 COEPI = COEPI / (1.-COEPI);
  966.                 'SI' (COEPI > 1.D0) ; COEPI=1.D0;'FINS';
  967.                 COEPI0=COEPI;STAB12.'AUTOCOEF'=COEPI;
  968.            'SINON';
  969.                 STAB12.'AUTOCOEF' = 1.D0;
  970.            'FINS' ;
  971.            RED1 = 1. ;
  972.            RED2 = 0 ;
  973.            'SI' ('NEG' WTAB.'AUTORED1' 'INCONNU') ;
  974.                 'SI' (STAB12.'AUTORED1' > 0);
  975.                      STAB12.'AUTORED1' = STAB12.'AUTORED1' - 1;
  976.                      'SI' (STAB12.'AUTORED1' 'EGA' 0) ;
  977.                            COEPI = 3 *  COEPI ;
  978.                            STAB12.'AUTOREDU' = STAB12.'AUTOREDU' / 3.;
  979.                            'SI' (STAB12.'AUTOREDU' > 1.1 );
  980. *                      on travaille encore avec un critere reduit
  981.                                 STAB12.'AUTORED1' = 4 ;
  982.                                 RED1 = 3. ;
  983.                            'FINS';
  984.                      'MESS' 'On multiplie le critere de pilotage par 3';
  985.                      'FINS';
  986.                 'FINS';
  987.            'FINS';
  988.            'SI' (COEPI > 1d0);
  989.                  RED1 = RED1 / COEPI ;COEPI =1d0;
  990.            'FINS';
  991.            'SI' ( 'NEG' WTAB.'NBPLAS' 'INCONNU') ;
  992.                  'SI' ( WTAB.'NBPLAS' 'EGA' 0) ;
  993.                     'SI' (COEPI < 0.) ;COEPI = COEPI * -2.;'FINS';
  994.                  'FINS';
  995.            'FINS';
  996.            COEPI = 'ABS' COEPI ;
  997.            COEPI0 = COEPI;
  998. *    sans pilotage
  999.      'SINON';
  1000.            STAB12.'AUTOCOEF'= 1.D0;
  1001.      'FINS';
  1002.  
  1003. ************************************************************************
  1004. *-------------    quelques initialisations  pour la boucle ETIQ ********
  1005. ************************************************************************
  1006.      URG    = FAUX; RED_URG = 0 ; IT = 0 ; c_zdepr = faux;ITACC = 0;
  1007.      ZICONV = VRAI; MMC     = 0 ; MMCMAX = 0 ; EPSM = 0.; DPSMAX = xpetit; dpsmaxp = dpsmax;
  1008.      DEPSTDM  = 0. ;  DEKREAC1 = 0. ;XCONVNOR = 0. ; ITNORM1  = 0 ;
  1009.      DEPSTREF = 100. * (WTAB . 'MAXDEFOR') ;
  1010.      ZDEPL='EXCO' conti.'DEPLACEMENTS' 'LX' 'NOID' 'LX';
  1011.      GR_U_K   = GR_U_DEB ;
  1012.      DITNORM1 = 0  ;NBCYCLE1 = 0  ;
  1013.      zdept = zu1 - (zu1 'ENLE' 'LX');
  1014.      zdeptq = zdept; zdeptp = zdept;
  1015.      zprecnc=1e-5;  FDEF = 0. ; ITNV = -5 ;
  1016.      TABCONV = 'TABL';
  1017.      'SI' IFEFPUL;
  1018.         XUPDA = 1;
  1019.      'FINS';
  1020.  
  1021. *********** Corrections en cas de materiaux variables ****************************
  1022.       DEPST0=0 ;
  1023.      'SI' WTAB.'MATVAR' ;
  1024.        'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  1025.           XXX3 ='ELAS' ZMODL ZSIG0 ZMAT  ZVAR0;
  1026.           XXX4 ='ELAS' ZMODL ZSIG0 ZMAT1 ZVAR0;
  1027.        'SINON';
  1028.           XXX3 ='ELAS' ZMODL ZSIG0 ZMAT       ;
  1029.           XXX4 ='ELAS' ZMODL ZSIG0 ZMAT1      ;
  1030.        'FINS';
  1031.         DDEF0 = XXX4 - XXX3;
  1032.        'DETR' XXX3;
  1033.        'DETR' XXX4;
  1034.         DEPST0=-1.* DDEF0;
  1035.      'FINS';
  1036.  
  1037. ***** En cas de chargement thermiques *********************************
  1038.       DMSRT0=0;
  1039.      'SI' ITHER ;
  1040.        'SI' (WTAB.'MATVAR' 'ET' IPILOT);
  1041.          'MESS' 'Le pilotage n est pas possible avec un materiau qui depend de la temperature' ;
  1042.          'ERREUR' 19 ;
  1043.        'FINS';
  1044.  
  1045.        'SI' ('EGA' ('TYPE' TETA2) 'CHPOINT');
  1046.           MCHTETA2 ='CHAN' 'CHAM' TETA2 ZMODL 'STRESSES' 'TEMPERATURES' ;
  1047.        'SINON';
  1048.           MCHTETA2 = TETA2;
  1049.        'FINS';
  1050.         ETREF= 'REDU' WTAB.'ETREF' ZMODL ; 'COMM' 'REDU au cas ou le modele soit parallele';
  1051.         ETT  =('EPTH' ZMODL ZMAT MCHTETA2 ) - WTAB.'ETREF' ;
  1052.         DTT  =  ETT '-' ETT0 ;
  1053.  
  1054.        'SI' WTAB.'POR1';
  1055. *------   Cas du milieu poreux avec chargement thermique ----------
  1056. *         cas isotrope seulement pour le moment
  1057. *         et on  ne s'occupe pas du alpha-reference !!
  1058.           MSRTT = PAS_EPTH PRECED ZMODL ZMAT MCHTETA2 ;
  1059.           DMSRT0= MSRTT '-' MSRTT0 ;
  1060.        'FINS' ;
  1061.  
  1062.         DEPST0 = DTT '+' DEPST0 ;
  1063.      'FINS';
  1064.  
  1065. *     calcul de DSIGT0 et FTHE en cas de chargement si necessaire ***********
  1066.       FTHE = 0 ;
  1067.      'SI' (DEPST0 'NEG' 0);
  1068.        'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  1069.           DSIGT0 = 'ELAS' ZMODL DEPST0 ZMAT ZVAR0;
  1070.        'SINON';
  1071.           DSIGT0 = 'ELAS' ZMODL DEPST0 ZMAT      ;
  1072.        'FINS' ;
  1073.  
  1074.        'SI' (DMSRT0 'NEG' 0);
  1075.           DSIGT0 = DSIGT0 '+' DMSRT0 ;
  1076.        'FINS';
  1077.         FTHE ='BSIG' ZMODL DSIGT0 ZMAT ;
  1078.      'FINS';
  1079.  
  1080. *-------------- deformations imposes ==> contraintes imposees **********************************
  1081.      'SI' LOGDEF;
  1082.            DDEFOR0 ='REDU' (STAB12.'DEFOR2' - STAB12.'DEFOR1') ZMODL;
  1083.           'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  1084.              DSI1 = 'ELAS' ZMODL DDEFOR0 ZMAT ZVAR0;
  1085.           'SINON';
  1086.              DSI1 = 'ELAS' ZMODL DDEFOR0 ZMAT      ;
  1087.           'FINS';
  1088.            FDEF = 'BSIG' ZMODL DSI1 ZMAT ;
  1089.      'FINS';
  1090.  
  1091. *-------------- pression imposee et grands deplacements ?***************
  1092.      'SI' (LOGPRE 'ET' IGRD) ;
  1093.        MOP = WTAB.'MOD_PRE' ;
  1094.        ZPEXT0 = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TEMPS0 ;
  1095.        ZPEXTF = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TI     ;
  1096. *
  1097.        'SI' ('EXIS' WTAB 'MAT_PRE') ;
  1098.           ZFPEXT0 = 'BSIG' MOP ZPEXT0 ('REDU' ZMAT11 MOP) ;
  1099.           ZFPEXTF = 'BSIG' MOP ZPEXTF ('REDU' ZMAT22 MOP) ;
  1100.        'SINON' ;
  1101.           ZFPEXT0 = 'BSIG' MOP ZPEXT0 ;
  1102.           ZFPEXTF = 'BSIG' MOP ZPEXTF ;
  1103.        'FINS' ;
  1104.        ZFP0F = 'COPIER' ZFPEXTF ; COEFP=1.D0;
  1105.      'FINS';
  1106.  
  1107. * ---------  ktangent et fefp******************************************
  1108.      'SI' IKTAN ;
  1109.            'SI' IFEFP ;
  1110.                 IKT_SAUV = VRAI ;
  1111.                 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  1112.                        'MESS' 'FEFP: Start with LASTKTAN' ;
  1113.                      ZRIKTA = STAB12.'LASTKTAN' ;
  1114.                      ZRAID = ZCLIM 'ET' ZRIKTA ;
  1115.                 'SINON' ;
  1116.                      'MESS' 'FEFP: Previous KTAN not available' ;
  1117.                      ZRAID = ZRAID 'ET' ('KSIGMA' ZMODL ZSIG0 ZMAT) ;
  1118.                 'FINS' ;
  1119.            'SINON' ;
  1120.                 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  1121.                      IKT_SAUV = VRAI ;
  1122.                      'SI' IPERT ;
  1123.                          'MESS' 'Matrice tangente par perturbation - '
  1124.                          'Demarrage avec KTAN = LASTKTAN' ;
  1125.                      'SINON' ;
  1126.                          'MESS' 'Matrice tangente "coherente" - '
  1127.                          'Demarrage avec KTAN = LASTKTAN' ;
  1128.                      'FINS' ;
  1129.                       ZRIKTA = STAB12.'LASTKTAN' ;
  1130.                       ZRAID = ZCLIM 'ET' ZRIKTA ;
  1131.                 'SINON' ;
  1132.                      'SI' IPERT ;
  1133.                           IKT_SAUV='NEG' WTAB.'K_TANGENT_ITER0'
  1134.                                    'MAT_ELASTIQUE';
  1135.                          'MESS' 'Matrice tangente par perturbation - '
  1136.                          'Demarrage avec KTAN = rigidite elastique' ;
  1137.                      'SINON' ;
  1138.                IKT_SAUV = ('NEG' WTAB.'K_TANGENT_ITER0' 'MAT_ELASTIQUE')
  1139.                'ET' ('NEG' WTAB.'K_TANGENT_ITER0' 'MAT_TANGENTE') ;
  1140.                          'SI' ('EGA' WTAB.'K_TANGENT_ITER0'
  1141.                               'MAT_ELASTIQUE') ;
  1142.                                'MESS' 'Matrice tangente "coherente" - '
  1143.                             'Demarrage avec KTAN = rigidite elastique' ;
  1144.                          'SINON' ;
  1145.                                'MESS' 'Matrice tangente "coherente" - '
  1146.                                'Demarrage avec KTAN (DTTAN = 0.)' ;
  1147.                                DTTAN = 0. ;
  1148.                                ZRIKTA = 'KTAN' ZMODL ZSIG0 ZVAR0 ZMAT
  1149.                                     'PREC' ZPREK 'DT  ' DTTAN ZKTASYM ;
  1150.                                ZRAID = ZCLIM 'ET' ZRIKTA ;
  1151.                          'FINS' ;
  1152.                      'FINS' ;
  1153.                 'FINS' ;
  1154.            'FINS' ;
  1155.      'FINS' ;
  1156.  
  1157. *--------  en grands deplacements option K_SIGMA ***********************
  1158. *
  1159.      'SI' (IGRD 'ET' ('NON' HPP_EPS));
  1160.        'SI' (IKSIA 'ET' ('NON' IFEFP)) ;
  1161.           KSI1    ='KSIGMA' ZMODL ZSIG0 ZMAT;
  1162.           ZRAIDINI= ZRAID           ;
  1163.           ZRAID   = ZRAID 'ET' KSI1 ;
  1164.        'FINS' ;
  1165.      'FINS' ;
  1166.  
  1167. * Y a-t-il des forces non conservatives ( forces suiveuses)? ***********
  1168.       ADDISEC0 = FAUX;
  1169.       ADDISEC2 = FAUX;
  1170.      'SI' WTAB.'PROCEDURE_CHARMECA';
  1171. *       on ajoute l indice ADDI_MATRICE pour signaler a charmeca qu on
  1172. *       souhaite aussi l operateur linearisé des Forces NL de charmeca
  1173.         PRECED . 'ADDI_MATRICE' = vrai;
  1174.         TFP22  = CHARMECA  PRECED   WTAB.'T_FINAL';
  1175.         PRECED . 'ADDI_MATRICE' = faux;
  1176. *   FP22 = F^suiv_n+1
  1177.          DMZPRES = 0.D0 ;
  1178.        'SI' ('EXIS' TFP22 'ADDI_SECOND') ;
  1179.           FP22    = TFP22.'ADDI_SECOND'  ;
  1180.           FP022   = 'COPIER' FP22        ;
  1181.           MZPRES  = 'MAXI' 'ABS' FP22    ;
  1182.           DMZPRES = MZPRES               ;
  1183.           ADDISEC2= VRAI                 ;
  1184.        'FINS';
  1185.        'SI' ('EXIS' TFP22 'ADDI_MATRICE');
  1186.           ZRAID = ZRAID 'ET' TFP22.'ADDI_MATRICE';
  1187.        'FINS';
  1188.  
  1189. *  FP0 = F^suiv_n
  1190.         TFP0 = CHARMECA  PRECED   TEMPS0 ;
  1191.        'SI'('EXIS' TFP0 'ADDI_SECOND')   ;
  1192.           FP0     = TFP0.'ADDI_SECOND'   ;
  1193.           MZPRES0 = 'MAXI' 'ABS' FP0     ;
  1194.           DMZPRES = DMZPRES '-' MZPRES0  ;
  1195.           ADDISEC0= VRAI                 ;
  1196.        'FINS';
  1197. *      'SI' ('EXIS' TFP0 'ADDI_MATRICE');
  1198. *         ZRAID = ZRAID 'ET' TFP0.'ADDI_MATRICE';
  1199. *      'FINS';
  1200.         COEFP   = 1.D0              ;
  1201.      'FINS';
  1202.  
  1203.  
  1204. * -----------calcul de la partie constante du second membre  **********
  1205. *                          en consolidation
  1206. *     ZFP1 est cense contenir : - B0*SIG0 et
  1207. *                                 DT*(1-TETA)*FI0 + DT*H*P
  1208. *
  1209. *     dans ZFCONSTA on met le second membre de u ***********************
  1210.  
  1211. *  en dynamique ********************************************************
  1212. * ZFP1 est cense contenir : F0 + 4/DT*M*V0 - B0*SIG0
  1213.      'SI' IDYN  ;
  1214.         UNSURH  = 1.D0 '/' STAB12.'DT'    ;
  1215.         ZFP1    = WTAB.'FREA1'            ;
  1216.         ZDYFEXT = ZFCONSTA 'ENLEVER' 'FLX';
  1217.         ZFCONSTA= ZFCONSTA '+' ZFP1       ;
  1218.      'FINS';
  1219.       ZFEXT     = ZFCONSTA 'ENLEVER' 'FLX';
  1220.  
  1221. *---------deplacement (ou jeu) e imposer e la fin du pas *************
  1222. * on separe les efforts ZFEXT (=F^ext_n+1)   ***************************
  1223. * et deplacement (ou jeu) ZFLX1 (u^imp_n+1) a imposer a la fin du pas
  1224.      ZFLX1 = 'EXCO' ZFCONSTA 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET';
  1225.      'SI' ('NEG' STAB12.'FFROT' 'INCONNU'); FFROT = STAB12.'FFROT';
  1226.      'SINON' ; FFROT=ZFEXT * 0; 'FINS';
  1227.  
  1228. *   calcul des forces externes deja equilibrees au debut du pas ********
  1229. *   par B*SIGMA : ZF1 = F^int_n = B*sigma_n + K^cst*u_n
  1230.      ZF1 = 'BSIG' ZMODL ZSIG0 ZMAT1;
  1231.      'SI' IRCON;
  1232.            ZF1 = ZF1 '+' ('REDU' MAI_CONS (RIG_CONS '*' ZU1));
  1233.      'FINS';
  1234.      'SI'  IDYN  ;  FFDYN = 'COPIER' ZF1; 'FINS';
  1235.      'SI' ISOL ;
  1236.            GRAP0= 'GRAD' MO_POR ZU1 MA_POR0 'CONS' ;
  1237.            XXX1 = 'GRAD' MO_POR ZU1 MA_POR  'CONS' ;
  1238.            XXXS =((1.- WTAB. 'TETA' )*GRAP0)+ (WTAB. 'TETA' * XXX1);
  1239.            XXX2 = STAB12.'DT' * ('GNFL' MO_POR XXXS) ;
  1240.            XXX3 = ZF1  ;
  1241.            ZF1  = XXX3 - XXX2;'DETR' XXX3;
  1242.            'DETR' XXX2 ;
  1243.      'FINS';
  1244.  
  1245. * initialisation des variables forces et deplacement*******************
  1246. * zzd est le deplacement au pas precedent   (=u_n)  et ZLX=lambda_n
  1247.      ZZD ='ENLE' ZU1 'LX';
  1248.      ZLX = ZU1 'EXCO' 'LX' 'NOID' 'LX' 'NATURE' 'DIFFUS';
  1249. * ---------  flxini est la partie des FLX deja realisee au debut du pas
  1250. **pv     FLXINI= ZZD '*' ZCLIM;
  1251. * FREAP  : -1*reactions du pas precedent    (=F^reac_n) transportees sur les nou
  1252. **pv     FREAP = ZLX '*' ZCLIM;
  1253.     FZU1 = ZU1 '*' ZCLIM;
  1254.     FLXINI = FZU1 'EXCO' 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET';
  1255.     FREAP  = FZU1 'ENLE' 'FLX';
  1256.  
  1257.  
  1258.      FEXT0 = ZF1 '+' FREAP;
  1259. * FEXT0 est le chargement externe (sans reactions vu par la
  1260. * structure le pas d'avant)   = F^ext_n = F^int_n - F^reac_n
  1261. *
  1262. *  on va calculer le premier residu c'est a dire le desequilibre *******
  1263. *  entre les forces externes et le calcul B*SIGMA.
  1264. *  le sigma qui sert est celui qui existerait si le champ de
  1265. *  deplacement ne changeait pas (ZU1). faire attention aux FLX
  1266. *  En pilotage on reprend ce residu que l'on multiplie par COEPI
  1267. *  XXX1 = [F^ext_n+1 ; Du^imp]
  1268. *  DFEXT0 = increment des forces et des FLX a imposer (le residu
  1269. *  du pas precedent) = [DF^ext ; Du^imp]
  1270.      XXX1   = ZFCONSTA '-' FLXINI ;
  1271.      DFEXT0 = XXX1     '-' FEXT0  ;
  1272.  
  1273. * si pression suiveuse dfext0 contient en plus l'increment des forces
  1274. * de pression du uniquement a la reactualisation de la geometrie (sans
  1275. * augmentation du module)
  1276. * mais comme F^int_n equilibre deja F^ext_n + F^suiv_n +... et qu'on est
  1277. * toujours sur config_n, on doit avoir : DFEXT0= [DF^ext ; Du^imp]
  1278. * avec DF^ext qui ne contient pas de forces suiveuses (...a vérifier)
  1279.     'SI' ADDISEC0;
  1280.        DFEXT0 = DFEXT0 '+' FP0     ;
  1281.     'FINS';
  1282.     'SI' (LOGPRE 'ET' IGRD) ;
  1283.        DFEXT0 = DFEXT0 '+' ZFPEXT0 ;
  1284.     'FINS' ;
  1285.      DFEXT0F = DFEXT0 'ENLEVER' 'FLX';
  1286.      DFEXT0L = DFEXT0 'EXCO' 'FLX' 'NOID' 'FLX';
  1287.  
  1288. *    DFEXT0 = [DFEXT0F ; DFEXT0L]
  1289. *           = [F^ext_n+1 - (F^int_n - F^reac_n - F^suiv_n) ; Du^imp]
  1290. *    RESIDU = forces exterieures sans reactions
  1291. *            (avec des termes supplementaires le cas echeant p.ex. en
  1292. *             dynamique ou en poreux) - forces interieures
  1293. *    et increment des relations imposees
  1294. *             ---> la resolution fournira dU et dR
  1295.  
  1296. *     RESIDU = [F^ext_n+1 + F^ther_n+1 + F^defi_n+1 - F^int_n ; Du^imp]
  1297.       RESIDU =  XXX1     '+'FTHE      '+'FDEF      '-'ZF1;
  1298.       ZDFINI ='COPIER' DFEXT0          ;
  1299.       ZFPLO  = ZF1 '-' FTHE '-' FDEF   ;
  1300.      'SI' (ITHER 'OU' WTAB.'MATVAR');
  1301.         ZDFINI=ZDFINI '+' FTHE;
  1302.      'FINS';
  1303.      'SI' LOGDEF;
  1304.         ZDFINI= ZDFINI '+' FDEF;
  1305.      'FINS';
  1306.      'SI' (LOGPRE 'ET' IGRD) ;
  1307.         RESIDU = RESIDU '+' ZFPEXTF ;
  1308.         ZDFINI = ZDFINI '+' ZFPEXTF '-' ZFPEXT0 ;
  1309.      'FINS';
  1310.  
  1311.      'SI' ADDISEC0 ;
  1312.         ZDFINI = ZDFINI '-' FP0 ;
  1313.      'FINS';
  1314.  
  1315.      'SI' ADDISEC2 ;
  1316. *       mess ' fp22 ' ; list resu fp22;
  1317.         RESIDU = RESIDU '+' FP22 ;
  1318. *       on tient compte de l'increment des forces suiveuses en direction
  1319. *       et en module
  1320.         ZDFINI = ZDFINI '+' FP22 ;
  1321.      'FINS';
  1322.  
  1323. * ici, on a :
  1324. * RESIDU = [F^ext_n+1 + F^ther_n+1 + F^defi_n+1 + F^suiv_n+1
  1325. *           - F^int_n  ;  Du^imp]
  1326. *        = [ DF^tot ; Du^imp]
  1327. * ZDFINI = [F^ext_n+1 + F^ther_n+1 + F^defi_n+1 + F^suiv_n+1
  1328. *           - (F^int_n - F^reac_n)  ;  Du^imp]
  1329. *
  1330.       IMPO12= FAUX;
  1331.       IMPO12U = FAUX;
  1332.       XXX1  ='EXTR' ZCLIM 'MAIL' 'UNIL';
  1333.      'SI' (('NBEL' XXX1) '>' 0);
  1334.         IMPO12 = VRAI;
  1335.         IMPO12U = VRAI;
  1336.         DIMPO12='REDU' DFEXT0L XXX1 ;
  1337.         DIMPOV = DFEXT0L '-' DIMPO12;
  1338.      'FINS';
  1339.       stab12.'SECOND_MEMBRE' = RESIDU '*' 1.D0;
  1340.  
  1341.  
  1342. ************************************************************************
  1343. ***               1ERE RESOLUTION                                    ***
  1344. ************************************************************************
  1345.      FEXCI = 'VIDE' 'CHPOINT' / 'DISCRET' ;
  1346.  
  1347.  
  1348.      'SI' (WTAB.'ADHERENCE') ;
  1349.        'SI' ('EXIS' WTAB 'MAIL_BLOM') ;
  1350.          FADRE = 'REDU' FADHE WTAB.'MAIL_BLOM' ;
  1351.        'SINON' ;
  1352.          FADRE = FADHE ;
  1353.        'FINSI' ;
  1354.        FEXCI = FEXCI 'ET' FADRE ;
  1355.      'FINSI' ;
  1356. *
  1357. *
  1358.  
  1359.      'SI' LOGPIL ;      
  1360. *    Modification du residu      
  1361.         RESIDU = RESIDU 'ET' (ETA0 '*' ZFPILIN);        
  1362.     'FINSI' ;
  1363. *  recalcul frottement
  1364. 'SI' (WTAB.'CAFROTTE' 'ET' IMPO12);
  1365.     'SI' (WTAB.'FROCABL');
  1366.        zsigfT = 'REDU' zsigf ZMODLI;
  1367.        FFROT = 'EXCC' ZCLIM ZDEPT MODFRO MATFRO ZSIGFT;
  1368.        RTRSF = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  1369.        RFNS  = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  1370.        FEXCI = FEXCI '+' FFROT ;
  1371.      'FINSI';
  1372.     'SI' (WTAB.'FROCOUL');
  1373.         excfconv = excfconv + 1 ;  
  1374.        RTRSF*'RIGIDITE' RFNS*'RIGIDITE'  FFROT = 'EXCF' ZCLIM (ZDEPT          ) MODFRO MATFRO excfconv;
  1375.        FEXCI = FEXCI '+' FFROT ;
  1376.      'FINSI';
  1377.  
  1378. 'FINS';
  1379.  
  1380.  
  1381.      'SI' ('EXIS' STAB12 'RIBLO_M' ) ;
  1382.          'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12);
  1383.             ZDEP1 BID BID BID BID = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  1384.                                             STAB12.'RIBLO_M'
  1385. **                                           STAB12.'LISEA_M' FEXCI ;
  1386.                                             STAB12.'LISEA_M' FEXCI RFNS RTRSF ;
  1387.          'SINON';
  1388.             ZDEP1 BID BID BID     = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  1389.                                             STAB12.'RIBLO_M'
  1390.                                             STAB12.'LISEA_M'      ;
  1391.          'FINS';
  1392.          'OUBLIER' STAB12 'RIBLO_M';
  1393.      'SINON';
  1394.          'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12U);
  1395. **               ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU FEXCI           ;
  1396.                  ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU FEXCI RFNS RTRSF          ;
  1397.          'SINON';
  1398.                ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU;
  1399.          'FINS';
  1400.      'FINS';
  1401. **   mess 'zdep1 apres resou' ;list (zdep1 exco 'LX' 'LX');
  1402.  
  1403.      'SI' LOGPIL ; 
  1404.        ZDEPII= 'RESO' ZRAID ZFPILIN ; 
  1405.        ZDEPILO = 0. * ZDEPII ;     
  1406. *    Calcul de D_eta par appel a la procedure PILOINDI
  1407.            D_ETA = PILOINDI PRECED ZU1  ZDEPILO ZDEP1 ZDEPII DTAU; 
  1408. *    Mise à jour  
  1409.            ZDEP1 = ZDEP1 '+' (D_ETA '*' ZDEPII) ;       
  1410.        ETA = ETA0 ;  
  1411.     'FINSI' ;
  1412.  
  1413.  
  1414.       ZRAID_T = 'EXTR' ZRAID 'CONT';
  1415.  
  1416.      'SI'  IDYN;
  1417.         STAB12.'ZRAIDV'= ZRAID;
  1418.      'FINS';
  1419.  
  1420.      'SI' ('EXIS' ZRAID_T 'NITER');
  1421.           STAB12.'RIBLO_M' = ZRAID_T. 7 ;
  1422.           STAB12.'LISEA_M' = ZRAID_T. 6 ;
  1423.      'FINS';
  1424. *  pour eventuel calcul de augauto
  1425.  zdep1d = zdep1 'ENLE' 'LX';
  1426. *  et initialisation zdepl (sert pour xnum)
  1427.  ZDEPL = 'EXCO' ZDEP1 'LX' 'NOID' 'LX' ;     
  1428. * On sauve le deplacement initial pour la convergence forcee ***********
  1429.      'SI' ('EXIS' WTAB 'DEPI');
  1430.      'SINON';
  1431.      'FINS';
  1432. *
  1433. *        calcul d'une norme pour la convergence**************************
  1434. *
  1435.       XXX1= ZFEXT;
  1436.      'SI' ADDISEC2;
  1437.        XXX1=ZFEXT + FP22;
  1438.      'FINS';
  1439.      'SI' (LOGPRE 'ET' IGRD) ;
  1440.        XXX1 = ZFEXT + ZFPEXTF ;
  1441.      'FINS';
  1442.  
  1443.  
  1444.      'SI' LOGPIL ;      
  1445. *    Modification du chargement
  1446.         XXX1 = XXX1 '+' ((ETA '+' D_ETA) '*' ZFPILIN);  
  1447.     'FINSI' ;
  1448.  
  1449.  
  1450.      ZDEP1P50 = ZDEP1       ;
  1451.      XDENO='XTY' ZDEP1P50  ( XXX1 -( RESIDU  'EXCO'
  1452.         'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET')) MLPRIM MLDUAL;
  1453.      'SI' ('VERI' xdeno) ; 'SINON'; xdeno = 1; 'FINSI';
  1454.      MZDEP1M = 'MAXI' ZDEP1P50 'ABS' 'AVEC' MLDEPL;
  1455.      MZFM = 'MAXI' (FTHE + FDEF + ZF1) 'ABS' 'AVEC' MLDUAL;
  1456.      XDENO1 = 'ABS' XDENO  + (MZFM * MZDEP1M);
  1457.      MZDEP1M = MZDEP1M + XPETIT;
  1458.      XDENO=XDENO1/MZDEP1M;
  1459.      XDENO = XDENO + MZFM;
  1460.      'SI' (XDENO < XPETIT); XDENO = 1.; 'FINSI';
  1461.      XDENOM=XDENO;
  1462.      'SI' TSTMOM ;
  1463.         ZDEP1P50 = ZDEP1 + XPETIT;          ;
  1464.         XDENOM = XDENO1/('MAXI' ZDEP1P50 'ABS' 'AVEC' MLROTA);
  1465.      XDENOM = XDENOM + XPETIT       ;
  1466.      'FINS' ;
  1467.      'SI' IPILOT ;
  1468.            'SI' ( WTAB.'AUTODEUX' ) ;
  1469.               XDENO = STAB12.'XDENO';
  1470.               XDENOM = STAB12.'XDENOM';
  1471.            'FINS';
  1472.      'FINS';
  1473.      'SI' ('NON' WTAB.'CONV');
  1474. *        'MESS' 'non convergence on garde xdeno xdenom ' xdeno xdenom;
  1475.             'SI' (XDENO < STAB12.'XDENO'); XDENO = STAB12.'XDENO'; 'FINSI';
  1476.             'SI' (XDENOM < STAB12.'XDENOM'); XDENOM = STAB12.'XDENOM'; 'FINSI';
  1477.      'FINS';
  1478.  
  1479.      IAFAIR=FAUX;
  1480.      'SI'  WTAB.'CONV';
  1481.      STAB12.'INCREMENT' = 'COPIER' ZDFINI ; 'FINS';
  1482.      RESIDNOR = 'COPIER' RESIDU ;
  1483.      'SI' IPILOT;
  1484.             XXX3=DFEXT0F * ( 1-COEPI) ;
  1485.             XXX2 = ZFEXT - XXX3;'DETR' XXX3;
  1486.            'DETR' ZFEXT; ;ZFEXT = XXX2;
  1487.             XXX1 = RESIDU * COEPI; 'DETR' RESIDU;
  1488.             XXX2 = 1.D0 -COEPI * FREAP;RESIDU = XXX1 + XXX2;
  1489.            'SI' IMPO12;
  1490.                 RESIDU = 1.D0 - COEPI * DIMPO12 + RESIDU;
  1491.            'FINS';
  1492.            IAFAIR=VRAI;
  1493.      'FINS';
  1494. *
  1495. *  petite correction du residu pour esperer gagner du temps ************
  1496. *
  1497.      INIT = FAUX ;
  1498.      'SI'(('NEG' STAB12.'FNONL' 'INCONNU') 'ET'
  1499.            (WTAB.'INITIALISATION'));
  1500.         'SI' IPILOT;
  1501.               'SI' ( WTAB.'AUTODEUX' 'ET' (COEPI 'NEG' 1.D0)) ;
  1502.               'MESS' 'Initialisation a partir du pas precedent '  COEPI;
  1503.                     IAFAIR=VRAI; INIT = VRAI;
  1504.                      XXX1= RED1 * STAB12.'FNONL';
  1505.                      XXX2= XXX1 + RESIDU ;
  1506.                     'DETR' RESIDU ;'DETR' XXX1 ;RESIDU =XXX2;
  1507.                'FINS' ;
  1508.         'SINON';
  1509. * on fait la correction si le pas precedent a converge sans non convergence
  1510.             'SI' (WTAB.'CONV' 'ET' (ISOUSPPP 'EGA' 0));
  1511. * on fait la correction si le pas precedent etait non lineaire.
  1512.               'SI' (('MAXI' 'ABS' STAB12.'FNONL') > (ZPREC * XDENO)) ;
  1513. * on enleve le residu du pas precedent pour recuperer l'increment
  1514. * nominal du second membre e imposer f1.
  1515. * f2 et l'increment du second membre du pas precedent
  1516.                 STAB12.'INCREMENT'=STAB12.'INCREMENT' - STAB12.'RESIDU';
  1517.                     zdeps = WTAB.'ZDEP1' + zdep1;
  1518.                     FFNO= 'XTY' STAB12.'FNONL' zdeps  MNDUAL MNPRIM;   
  1519.                     F1 = STAB12.'INCREMENT' ;
  1520.                     F2 = INCRPREC ;
  1521.                     f12 = 'XTY' f1 zdeps MNDUAL MNPRIM;                
  1522.                     f22 = 'XTY' f2 zdeps MNDUAL MNPRIM;                
  1523.                     AMPL = f12/(f22 + XPETIT);
  1524.                     DTPREC=1E30; AMPLT = 0;
  1525.                     'SI' ('NEG' STAB12.'DTPREC' 'INCONNU');
  1526.                          DTPREC = STAB12.'DTPREC';
  1527.                          'SI' (DTPREC > XPETIT);
  1528.                              AMPLT=WTAB.'DT_INIT' /DTPREC;
  1529.                          'FINS';
  1530.                      'FINS';
  1531. *  Le chargement n'est il pas de fluage ou de thermique ?';
  1532.                     XDCOMP = ('XTY' ZDEP1 F1 MNPRIM MNDUAL) ;  
  1533.                    'SI' (('ABS' XDCOMP) < (ZPREC * XDENO * mzdep1m)
  1534.                         'OU' (('ABS' FFNO) > (('ABS' F22) * 2.e2   )));
  1535.                         DTPREC = STAB12.'DTPREC';
  1536.                         AMPL=AMPLT;
  1537. * la decharge est-elle significative
  1538.                       'SI'((F12/(f22 + XPETIT)) < -0.05);ampl=0.;'FINS';
  1539.                         STAB12.'INITEMPS'=VRAI;LOGTEMP=FAUX;;
  1540.                         'MESS'
  1541. 'Pas d increment de charge, initialisation calculee avec le temps';
  1542.                    'SINON';
  1543.                         AMPL = F12 / (F22 + XPETIT);
  1544.                         LOGTEMP=STAB12.'INITEMPS';
  1545.                         STAB12.'INITEMPS'=FAUX;
  1546.                    'FINS';
  1547. * changement de modele on n'initialise pas
  1548.                    'SI' ('NEG' ZMODLI ZMODLP);  AMPL = 0; 'FINSI';
  1549.  
  1550.                     AMPL = MINI (prog AMPL AMPLT);
  1551.                    'SI' ((AMPL > 0) 'ET' (AMPL < 2e1)) ;
  1552.       'MESS' 'Initialisation a partir de la solution precedente Coeff'
  1553.                          AMPL;
  1554.                          XXX1 = AMPL * STAB12.'FNONL';
  1555.                          XXX2 =RESIDU+ XXX1;
  1556.                          'DETR' XXX1;'DETR' RESIDU;
  1557.                          RESIDU = XXX2;
  1558.                          IAFAIR=VRAI; INIT = VRAI;
  1559.                    'FINS';
  1560.               'FINS';
  1561.             'FINS';
  1562.         'FINS';
  1563.      'FINS';
  1564.      WTAB.'ZDEP1'=zdep1;
  1565. *
  1566. *  initialisation en plastique *****************************************
  1567. *
  1568.      'SI' IPLAVI ;
  1569.           'SI' ('NEG' WTAB.'MOVA' 'RIEN') ;
  1570.                   ACC0 = 'EXCO' (WTAB.'MOVA') ZVAR0 ;
  1571.           'FINS' ;
  1572.      'FINS' ;
  1573. *
  1574. *    debut des iterations internes boucle etiquette ********************
  1575. *
  1576.      IKT =  FAUX ;
  1577.      IPREM =  VRAI ; RECA_K = FAUX; RECA_N = 0;
  1578.      DEPST  = 'ZERO' ZMODL 'DEFORMATIONS' ;
  1579.      DEPSTP = depst ;
  1580.      DEPSTK = DEPSTP ;
  1581. *
  1582. *       initialisation acceleration de convergence ********************
  1583. *
  1584.      iafair = vrai;
  1585.      PASTEST = FAUX;
  1586.      ZITAC= 0 ;
  1587. *     on peut mettre n'importe quoi c'est pour
  1588. *     ne pas faire de tests dans la boucle
  1589.      ACFP1  = 'COPIER' ZFEXT2 *0. ;
  1590.      ACFP2  = ACFP1 ;
  1591.      ZDEPLD = ACFP1 ;
  1592.      ACFP3  = ACFP1 ;
  1593.      ACFEP1 = ACFP1 ;
  1594.      ACFEP2 = ACFP1 ;
  1595.      FCORF  = 'COPIER' FREAP;
  1596. * initialisation du meilleur critere
  1597.      XCONVMIN = 1e20;
  1598.      DPSMREF  = 0  ;
  1599.      XCONV    = 0. ;
  1600.      XCONVP   = 1. ;
  1601.      NSOINCR  = 1  ;
  1602.      'SI' ('NEG' WTAB.'SOUS_INCREMENT'  'INCONNU') ;
  1603.          NSOINCR = WTAB.'SOUS_INCREMENT' ;
  1604.      'FINS';
  1605.      NONCONV = FAUX;
  1606.      CORREC=0;
  1607.      PASREINI=VRAI;
  1608.      coefmt=1.;
  1609. *
  1610. * initialisation des messages pour l'iteration en cours ****************
  1611. *
  1612.      'SI' IKLFFF ;
  1613.           'MESSAGE'
  1614. 'Iter'*13   'Nplas'*26    'Critere'*39        'Deps.max'*52       'Eps.max'*65      'Crit.flex'*78
  1615. ;
  1616.      'SINON';
  1617.          'MESSAGE'
  1618. 'Iter'*13   'Nplas'*26     'Fresidu'*39        'Deps.max'*52       'Eps.max'*65      'Mresidu'*78
  1619. ;
  1620.      'FINS';
  1621.  INSTAB = FAUX;
  1622. RESIDIN = RESIDU * 1.;
  1623. hpp_exit = faux; 'SI' hpp_eps; hpp_exit=vrai; 'FINSI';
  1624.  
  1625. *=======================================================================
  1626. *======= DEBUT DE LA BOUCLE DE CONVERGENCE                       =======
  1627. *=======================================================================
  1628. *  trac cach v1 nclk;
  1629. *  trac (surf1 et surf2 et surf3 et _rext et _rint) nclk;
  1630. *  trac cach (engr1 et engr2) nclk;
  1631. dpsmaxp = 0;
  1632.  
  1633. IT_RECA = 0;
  1634. depst = depst * 0;
  1635. depstp = depst;
  1636. 'REPETER'  ETIQ ;
  1637. resmul = 2. * resmul;
  1638. 'SI' (resmul > 1.0); resmul = 1.0; 'SINON'; itacc = 4; 'FINSI';
  1639. **resmul = 1.;
  1640. *  sauvegarde des etats sur la configuration debut de sous-increment: geom1
  1641. *  mise a jour apres la sortie de etiq
  1642. depstp = depst;
  1643. zsig0 = zsig0_1;
  1644.  
  1645. INSTAB = FAUX;
  1646.  
  1647. nconvr = faux;
  1648. *IT est le compteur de ETIQ, ITACC doit etre =< 0 pour qu'on accelere
  1649. IT= IT + 1 ;
  1650. ITACC = ITACC - 1;
  1651. ZITAC = ZITAC + 1 ;
  1652. *
  1653. *---------------------------------------------------------------------
  1654. *           La force motrice de l'iteration est fixee: RESIDU
  1655. *           on va calculer un nouveau champ de deplacement
  1656. *
  1657. *   calcul de l'increment de l'increment de deplacement zdep1
  1658. *   par resolution lineaire
  1659. *-----------------------------------------------------------------------
  1660. *
  1661. *  petits travaux pour acceleration de convergence
  1662. *
  1663. CORRECP = CORREC;
  1664. CORREC = 0;
  1665. PASTEST=FAUX;
  1666. ACFEP0 = RESIDU - FCORF;
  1667. ACFP0 = ACFEP0 'ENLE' 'FLX';
  1668. ACFEP0 = ACFP0;
  1669. ACFEP0 = ACFEP0 - CORRECP ;
  1670. 'SI' IGRD ;
  1671.   'FORM' GEOM1;
  1672. *  rem : 0.1 est la valeur par defaut de EKREAC (= . 'REAC_GRANDS')
  1673.   'SI' ((IT > 1) 'ET' (URG 'OU' (ITACC > 3) 'OU' INSTAB 'OU'
  1674.   (('MULT' (IT - IT_RECA)  ITRCLC) 'ET' ('NON' HPP_EPS))));
  1675.   URG = FAUX;
  1676.       IT_RECA = IT;
  1677.       dekreac1 = 0;
  1678.       INSTAB = FAUX;
  1679.       URG=FAUX;
  1680.         PASREINI=FAUX;
  1681.       GEOR ZMATTEMP = 'FORM' ZDEPT   ZMODLI  ZMATI ;
  1682. *     A cause de FORM, ZMATTEMP n'est plus parallele...
  1683.       ZMATTEMP = 'REDU' ZMATTEMP ZMODL;
  1684.       zdepr = zdept;
  1685.       c_zdepr = faux;
  1686.       RECA_K = VRAI;
  1687.       RECA_N = RECA_N + 1;
  1688.  
  1689.      'SI' (RECA_N > 20) ;
  1690.         nonconv = vrai;
  1691.      'FINS';
  1692.  
  1693.       txt_k ='CHAI' ' Recalcul de K (= K^el';
  1694.        'SOUC' 0;
  1695.        'SI' ((LAG_TOT 'EGA' 1) 'ET' vrai);
  1696.           'FORM' GEOREF0;
  1697.           HOOKRH = 'HOOK' ZMODL ZMATI;
  1698.           HOOKRH2 = 'PICA' HYPDEF HOOKRH ZMODL (ZU1 + ZDEPT);
  1699.           'FORM' GEOR;
  1700.           RITC = 'RIGI' HOOKRH2 ZMODL zmattemp 'NOER';
  1701.           'DETRUI' HOOKRH;
  1702.           'DETRUI' HOOKRH2;
  1703.           'FORM' GEOM1;
  1704.        'SINON';
  1705.          ritc  ='RIGI'   zmodl zmattemp 'NOER' ;
  1706.        'FINSI';
  1707.       PASOK = 'SOUCI';
  1708.       'SI' PASOK;
  1709.          'MESS' 'rigi rate. on reessaye avec la configuration geom1';
  1710.          'FORM' GEOM1;
  1711.          ritc  ='RIGI'   zmodl zmatI ;
  1712.       'FINSI';
  1713.  
  1714.  
  1715.  
  1716.       ZMATTEMP = 0 ;
  1717.  
  1718. *     RH peut contenir des CL
  1719.       ZRI = 'EXTR' RITC 'RIGI' 'NOMU' ;
  1720.       ZCL = 'EXTR' RITC 'RIGI' 'MULT' ;
  1721.       'SI' ('NEG' ('DIME' ZCL) 0) ;
  1722.         ZCLIM = ZCLIM0 'ET' ZCL ;
  1723.       'SINON' ;
  1724.         ZCLIM = ZCLIM0 ;
  1725.       'FINSI' ;
  1726.  
  1727.       BFCONT  = FAUX                       ;
  1728.       BCLIM2  = FAUX                       ;
  1729.      'SI' WTAB.'CONTACT';
  1730. *       recalcul raideur de contact et jeu
  1731. *   Maintenant, on veut la configuration finale pour les directions et les jeux
  1732.   'FORM' GEOM1;
  1733.   'FORM' ZDEPT ;
  1734.         MODCON  = wtab.'MODCONTA';
  1735.        'SI' ('EXIS' WTAB 'MATCONTA');
  1736.           MODFRO matfro CJEU CRR RFROT='RFCO' MODCON  WTAB.'CONV'  WTAB.'MATCONTA';
  1737.        'SINON' ;
  1738.           CJEU CRR RFROT='RFCO' MODCON  WTAB.'CONV'                ;
  1739.        'FINS' ;
  1740.  
  1741.        'SI' ('NEG' 0 CJEU) ;
  1742.          FADHE = 'EXCO' CJEU 'FADH' 'NOID' 'FADH' ;
  1743.          CDAP  = 'EXCO' CJEU 'FLX'  'NOID' 'FLX' ;
  1744.        'FINSI' ;
  1745.  
  1746.        'SI' (WTAB . 'MODAL') ;
  1747.           CRR  = 'PJBA' CRR  ZMODL MMMM ;
  1748.           MCRR = 'EXTR' CRR 'MAIL' 'MULT' ;
  1749.           MCDAP = 'EXTR' CDAP 'MAIL' ;
  1750.  
  1751.          'REPETER' BCDA ('NBNO' MCDAP) ;
  1752.             PBCDA = MCDAP 'POINT' &BCDA ;
  1753.             PCRR  ='POINT' MCRR &BCDA   ;
  1754.             CHCR  ='MANU' 'CHPO' PCRR 1 'FLX' ('EXTR' CDAP 'FLX' PBCDA)
  1755.                    'NATURE' 'DISCRETE' ;
  1756.            'SI' ('EGA' 1 &BCDA) ;
  1757.               CCDA = CHCR ;
  1758.            'SINON' ;
  1759.               CCDA = CCDA 'ET' CHCR ;
  1760.            'FINS' ;
  1761.          'FIN' BCDA ;
  1762.  
  1763.           CDAP = CCDA ;
  1764.        'FINS' ;
  1765.  
  1766. *       ATTENTION : mettre les conditions de frottement en premier pour
  1767. *                   les numeroter en dernier
  1768.         ZCLIM2 ='VIDE' 'RIGIDITE' / 'RIGIDITE';
  1769.         ZFCONT ='VIDE' 'CHPOINT' / 'DISCRET';
  1770.         CDEP = ZU1   + ZDEPT  ;
  1771.         CDEPSLX ='ENLE' CDEP 'LX'              ;
  1772.        'SI'  ('NEG' CRR 0);
  1773.           BCLIM2 = VRAI                ;
  1774.           BFCONT = VRAI                ;
  1775.           ZCLIM2 = CRR    'ET' ZCLIM2  ;
  1776.           CCOR   = CRR '*' CDEPSLX   'EXCO' 'FLX' 'FLX'  ;
  1777.           ZFCONT = ZFCONT '+'  CCOR    ;
  1778.        'FINS';
  1779.  
  1780.        'SI' ( 'NEG' 0 CJEU);
  1781.           BFCONT = VRAI                ;
  1782.           CDAP_M = CDAP EXCO FLX FLX EXTR MAILLAGE;
  1783.           ZFCONT = ZFCONT '+'  CDAP    ;
  1784.        'FINS';
  1785.  
  1786.        'SI' ('NEG' 0 RFROT);
  1787.           BCLIM2 = VRAI                ;
  1788.           BFCONT = VRAI                ;
  1789.           ZCLIM2 = RFROT  'ET' ZCLIM2  ;
  1790.         CDEPSLX ='ENLE' CDEP 'LX'              ;
  1791.           CCOR   = RFROT  '*'  CDEPSLX 'EXCO' 'FLX' 'FLX' ;
  1792.           ZFCONT = ZFCONT  +  CCOR    ;
  1793.        'FINS';
  1794.  
  1795. *      Mise a jour de ZCLIM, ZFCONSTA et RESIDU si necessaire
  1796.        'SI' BCLIM2;
  1797.          ZCLIM  = ZCLIM2 'ET' ZCLIM             ;
  1798.  
  1799. *  pas de correction a faire sur residu car il est hors reactions pour qu'elles
  1800. *  sortent en total a la resolution incrementale
  1801.  
  1802.        'FINS';
  1803.  
  1804.      'FINS';
  1805.  
  1806.      'SI' BFCONT;
  1807.         ZFCONSTA = ZFEXT2 '+' ZFCONT;
  1808.      'SINO';
  1809.         ZFCONSTA = ZFEXT2           ;
  1810.      'FINS';
  1811.       ZFCONST1 = ZFCONSTA;
  1812.  
  1813.          zclimp = zclim;
  1814.  
  1815.      'SI' IDYN;
  1816.         ZFCONSTA = ZFCONSTA '+' ZFP1  ;
  1817.      'FINS';
  1818.  
  1819.       ZFLX1 = 'EXCO' ZFCONSTA 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET';
  1820.       ZRAID =  ZCLIM 'ET' ZRI;
  1821. *  repasser sur la configuration de calcul des raideurs
  1822.      'FORM' geor;
  1823.      'SI' iksia ;
  1824.         sigttcca ='PICA' HYPDEF ZMODL ZSIGF ZDEPT ;
  1825.         ksigtc= 'KSIGM' sigttcca zmodl zmat;
  1826.         ZRAID= ZRAID 'ET' ksigtc;
  1827.      'FINS';
  1828.      'SI' IRCON;
  1829.         ZRAID = ZRAID 'ET' RIG_CONS;
  1830.         txt_k = 'CHAI' txt_k ' + K^cst';
  1831.      'FINS';
  1832. *  remise a jour impo12
  1833.       IMPO12= FAUX;
  1834.       IMPO12U = FAUX;
  1835.       XXX1  ='EXTR' ZCLIM 'MAIL' 'UNIL';
  1836.      'SI' (('NBEL' XXX1) '>' 0);
  1837.         IMPO12 = VRAI;
  1838.         IMPO12U = VRAI;
  1839.         DIMPO12='REDU' DFEXT0L XXX1 ;
  1840.         DIMPOV = DFEXT0L '-' DIMPO12;
  1841.      'FINS';
  1842. * pv
  1843.      'SI' IRAUG;
  1844. * xksx = xtmx zdep1d  ksigtc ;
  1845. *'SI' ((abs xksx) < xpetit );
  1846. *      zdep1d = (zu1 + zdept) 'ENLE' 'LX';
  1847. *      xksx = xtmx zdep1d  ksigtc ;
  1848. *'FINSI';
  1849.  
  1850.  'SI'  iprem; zdep1d = zu1 'ENLE' 'LX'; 'FINSI';
  1851.   xkx = (xtmx zdep1d rh)  + xpetit;
  1852.   xmx = (xtmx zdep1d rig_aug) + xpetit;
  1853.   'SI' (('VERI' xkx) 'ET' ('VERI' xmx) 'NON');
  1854.     zu1l = zu1 'ENLE' 'LX';
  1855.     xkx = (xtmx zu1l rh)  + xpetit;
  1856.     xmx = (xtmx zu1l rig_aug) + xpetit;
  1857.   'FINSI';
  1858. 'SI' faux;
  1859.   xktx = xty zdep1d acfp0  mnprim mndual;  
  1860.   xksx1 = xkx - xktx;
  1861. *  'SI' (xksx1 > xksx); xksx = xksx1 ; 'FINSI';
  1862.    xksx = xksx1;
  1863.   'SI' (('VERI' xksx) 'ET' ('VERI' xkx) 'ET' ('VERI' xmx) 'NON');
  1864.     zu1l = zu1 'ENLE' 'LX';
  1865.     xkx = xty zu1l (rh * zu1l) mnprim mndual + xpetit;    
  1866.     xmx = xty zu1l (rig_aug  * zu1l) mnprim mndual + xpetit; 
  1867.     xksx = xkx;
  1868. **   xksx = xty zu1l (ksigtc * zu1l) mnprim mndual;
  1869.   'FINSI';
  1870. * si les coef sont negatif, c'est qu'on n'est pas instable sur ce mode
  1871.  'SI' ((abs xksx) > xpetit );
  1872.     augauto = xksx / xmx  ;
  1873.     augk= xksx /xkx  ;
  1874.  'SINON';
  1875.   mess ' abs xksx ' (abs xksx) ' max zdep1d ' (maxi abs zdep1d);
  1876.  
  1877.  'FINSI';
  1878. 'FINSI';
  1879. *'SI' iprem;
  1880.       augauto = xkx / xmx; augk = 1.;
  1881. *'FINSI';
  1882.  
  1883.  
  1884.       augauto = abs augauto;
  1885.       augk = abs augk;
  1886.  
  1887.       augmult = augmult * 1.02;
  1888.       'SI' (augmult > 1D4 ); augmult=1D4 ; 'FINSI';
  1889.  
  1890.       augauto = augauto * augmult ;
  1891.       augk = augk * augmult;
  1892.  
  1893.  
  1894. 'SI' (iraug 'ET' autaug) ; 'MESS' 'multiplicateur d augmentation masse' augauto ' raideur' augk ; 'FINSI';
  1895.         ZRAIDNA=  ZRAID ;
  1896.        'SI' AUTAUG;
  1897.         ZRAID = ZRAID 'ET' (RIG_AUG * augauto) 'ET' (RH * augk) ;
  1898.        'SINON';
  1899.         ZRAID = ZRAID 'ET' RIG_AUG ;
  1900. ***      'MESS' ' augmentation residuelle ';
  1901.        'FINSI';
  1902.  
  1903.         txt_k = 'CHAI' txt_k ' + K^aug';
  1904. ***     ZNACCE = FAUX;
  1905.  
  1906.     'SINON';
  1907.       augmult = augmult * 0.55 ;
  1908. 'SI' iraug; 'MESS' 'nouveau augauto augmult' ' ' augauto augmult;  'FINSI';
  1909.      'SI' IRAUGLU;
  1910.         zraid = zraid 'ET' rig_aug;
  1911.      'FINSI';
  1912.     'FINSI';
  1913. ***     'FINS';
  1914.  
  1915.      'SI' WTAB.'PROCEDURE_CHARMECA';
  1916.         PRECED . 'ADDI_MATRICE' = vrai;
  1917.         TFP22= CHARMECA  PRECED   WTAB.'T_FINAL';
  1918.         PRECED . 'ADDI_MATRICE' = faux;
  1919.        'SI' (EXIS TFP22 'ADDI_MATRICE');
  1920.           zraid = zraid 'ET' TFP22.'ADDI_MATRICE';
  1921.           txt_k = 'CHAI' txt_k ' + K^cent';
  1922.        'FINS';
  1923.      'FINS';
  1924.  
  1925.      'SI' IDYN;
  1926. *      bp : en toute rigueur, il faudrait aussi recalculer la MASSE ...
  1927. *           et ajouter l'amortissement le cas échéant ...
  1928.         MMA  = WTAB.'MASSE' '*' ( 4.D0 '/' WTAB.'DT' '/' WTAB.'DT');
  1929.         ZRAID= ZRAID 'ET' MMA;
  1930.       'FINS';
  1931.       'MESS' ( 'CHAI' txt_k ' ) dans config deformee ' DEKREAC1 );
  1932.       'FORM' GEOM1;
  1933. *      on impose le recalcul de K a la prochaine iteration si it=2
  1934.       'SI' (IT 'EGA' 2) ;
  1935.          ITACC = 5 ;
  1936.       'SINON';
  1937. *  pas de raison d'interferer avec les accelerations de convergence
  1938. *****    ITACC = 1 ;
  1939.       'FINS' ;
  1940.        GR_U_K  = GR_U_FIN ;
  1941.        zdeptm = zdept ;
  1942.   'FINS';
  1943. 'FINS';
  1944.  
  1945. *
  1946. *  acceleration de convergence effective ------------------------------
  1947. *
  1948. ACCEL = FAUX;
  1949. 'SI' (('NON' instab) 'ET' (ITACC '&lt;EG' 0) 'ET' (IT '>' 3) 'ET' ZNACCE );
  1950.  ITACC = 2;
  1951.  ACCEL = VRAI;
  1952.    CORREC = 'ACT3' ACFEP2  ACFEP1 ACFEP0
  1953.                    ACFP3   ACFP2  ACFP1  ACFP0 ;
  1954.   'SI' (wtab.'STABILITE' 'ET' ('NON' HPP_EPS));
  1955. *    verif que l'acceleration ne renvoie pas en arriere
  1956.  
  1957.  
  1958.  
  1959.      acc_ref = xty (acfp0         ) zdep1d MNDUAL MNPRIM;    
  1960.      acc_ref = acc_ref + xpetit;
  1961.      acc_dir = xty (acfp0 - correc) zdep1d MNDUAL MNPRIM;    
  1962.      acc_rap = acc_dir/acc_ref;
  1963.     'SI' (acc_rap '>' 256.) ;
  1964. * en cas d'acceleration trop grande, on la limite
  1965.        'MESS' 'Limitation acceleration ' ' 'acc_rap;  CORREC = correc * (256./acc_rap);
  1966.      'FINSI';
  1967. *    'SI' ((acc_rap '&lt;EG' 0.) 'OU' (acc_dir '&lt;EG' 0.)) ;
  1968.      'SI' (acc_rap '&lt;EG' 0.) ;
  1969. * en cas d'acceleration retrograde, on n'accelere pas
  1970.        'MESS' 'Annulation acceleration: instable ' ' ' acc_dir;
  1971.        correc =  0.;
  1972.        itacc=2;
  1973.      'FINSI';
  1974.   'FINSI';
  1975.   RESIDU = RESIDU '-' CORREC;
  1976.   RESIDC = RESIDC '-' CORREC;
  1977. 'FINS';
  1978.  
  1979.  ACFP3  = ACFP2  ;
  1980.  ACFP2  = ACFP1  ;
  1981.  ACFP1  = ACFP0  ;
  1982.  ACFEP2 = ACFEP1 ;
  1983.  ACFEP1 = ACFEP0 ;
  1984. *
  1985. * Resolution ---------------------------------------------------------
  1986. * on obtient ZDEP1 = [ du ; Dlx ]
  1987. 'SI' (resmul < 0.99);
  1988. 'MESS' 'Reduction du chargement. Coefficient: ' resmul;
  1989. 'SI' (resmul < 1.d-8); 'ERREUR' 996; 'FINSI';
  1990. **ZDETOT = ZZD   '+' ZDEPT ;
  1991. **XXX1   = ZCLIM '*' ZDETOT;
  1992. * -1*reactions a l'iteration it
  1993. **FCORF = 'ENLEVER'  XXX1 'FLX';
  1994.   residu = ((residu -fcorf) * resmul) +  fcorf;
  1995. **  residu = (residu 'ENLE' 'FLX') 'ET' ((RESIDU 'EXCO' 'FLX' 'NOID' 'FLX')  * resmul);
  1996.     ACFP0 = (RESIDU  - FCORF) ;
  1997. **  ACFP0=ACFP0 + ((FFROT- FFROTP) * ZCLIM) ;
  1998. 'FINSI';
  1999.   'SOUCI' 0;
  2000. 'SI' IAFAIR;
  2001. *
  2002.  
  2003.   FEXCI = 'VIDE' 'CHPOINT' / 'DISCRET' ;
  2004.   'SI' (WTAB.'ADHERENCE') ;
  2005.     'SI' ('EXIS' WTAB 'MAIL_BLOM') ;
  2006.       FADRE = 'REDU' FADHE WTAB.'MAIL_BLOM' ;
  2007.     'SINON' ;
  2008.       FADRE = FADHE ;
  2009.     'FINSI' ;
  2010.     FEXCI = FEXCI 'ET' FADRE ;
  2011.   'FINSI' ;
  2012. *
  2013. *  recalcul frottement
  2014. 'SI' (WTAB.'CAFROTTE' 'ET' IMPO12);
  2015.     'SI' (WTAB.'FROCABL');
  2016.        zsigfT = 'REDU' zsigf ZMODLI;
  2017.        FFROT = 'EXCC' ZCLIM ZDEPT MODFRO MATFRO ZSIGFT;
  2018.        RTRSF = 'VIDE' 'RIGIDITE' / 'RIGIDITE';
  2019.        RFNS  = 'VIDE' 'RIGIDITE' / 'RIGIDITE';
  2020.        FEXCI = FEXCI '+' FFROT ;
  2021.      'FINSI';
  2022.     'SI' (WTAB.'FROCOUL');
  2023.         excfconv = excfconv + 1 ;  
  2024.        RTRSF*'RIGIDITE' RFNS*'RIGIDITE'  FFROT = 'EXCF' ZCLIM (ZDEPT           )  MODFRO MATFRO excfconv;
  2025. si faux;
  2026.   rfr = (zdept exco 'LX' 'NOID' 'LX') * rfrot REDU LOP1;
  2027.   rde = (zdept enle 'UZ')  REDU lop1;
  2028.   VR1      = 'VECT' RDE 'DEPL'  'ROUG'   ;
  2029.   VR2      = 'VECT' (ffrot *rfrot redu lop1)  'FORC' 'JAUN'  ;
  2030.   OEIL1    = (1.E5 -1.E5 1.E5) ;
  2031.   'TRAC' NCLK OEIL1 (VR1 et vr2) ('ARET' LOP1 ) 'TITR' 'Forces de reaction.' ;
  2032. finsi;
  2033.        FEXCI = FEXCI '+' FFROT ;
  2034.      'FINSI';
  2035.  
  2036. 'FINS';
  2037. zdeptp = zdept;
  2038.  
  2039. *
  2040.   'SI' ( ('EXIS' STAB12 'RIBLO_M') 'ET' (IPREM 'OU' RECA_K) ) ;
  2041.     'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12U);
  2042.        ZDEP1 BID BID BID BID = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  2043.                               STAB12.'RIBLO_M'   STAB12.'LISEA_M'
  2044.                               FEXCI RFNS RTRSF;
  2045. **                            FEXCI ;
  2046.     'SINON';
  2047.        ZDEP1 BID BID BID     = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  2048.                               STAB12.'RIBLO_M'   STAB12.'LISEA_M' ;
  2049.     'FINS';
  2050.  
  2051.   'SINON';
  2052.     'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12U);
  2053. **    ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU FEXCI ;
  2054. **   mess 'fexci 2 avant reso'; list fexci;
  2055.       ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU FEXCI       RFNS RTRSF;
  2056.     'SINON';
  2057.  
  2058.       ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU 'NOID';
  2059.       'SI' (WTAB.'MAN' 'ET' IPREM);
  2060.           ORDRE      = WTAB.'ORDRE' ;
  2061.           ZDEP2 IOUT = CORMAN ZRAIDINI ZMODL ZMAT  ORDRE
  2062.                        ZU1 ZSIG0      RESIDNOR  WTAB ;
  2063.           'SI' (IOUT 'EGA' 1) ;
  2064.              ZDEP1=ZDEP2;
  2065.           'FINS';
  2066.       'FINS';
  2067.     'FINS';
  2068.   'FINS';
  2069.   monsouc ='SOUCI';
  2070.   'SI' monsouc; 'SI' ('NON' autaug);
  2071.            'ERREUR' ('VALE' 'SOUCI');
  2072.    'FINSI'; 'FINSI';
  2073.  
  2074.  
  2075.      'SI' LOGPIL ; 
  2076. *      Calcul de D_eta par appel a la procedure PILOINDI
  2077.        D_ETA = PILOINDI PRECED ZU1 ZDEPT ZDEP1 ZDEPII DTAU ;       
  2078. *      Mise à jour  
  2079.            ZDEP1 = ZDEP1 '+' (D_ETA '*' ZDEPII) ;
  2080.            ETA = ETA '+' D_ETA ;
  2081.          'FINSI' ;
  2082.  
  2083.  
  2084. tabconv . it = 1;
  2085. *
  2086. *  verif sens ----------------------------------------------------------
  2087. *
  2088. si (resmul < 1d-2) ; souci 0; finsi;
  2089. zdep1s = zdep1 'ENLE' 'LX';
  2090. monsouc = souci;
  2091. pasunil = monsouc;
  2092. 'SI' (WTAB.'STABILITE' 'ET' (('NON' HPP_EPS) 'OU' monsouc) 'ET' autaug );
  2093.     sens = xtmx zraid zdep1s;
  2094.     'SI' ((sens < 0) 'OU' monsouc);
  2095.      'MESS' 'Annulation sous-increment: instable ' ' ' sens;
  2096.      itacc = 3  ;
  2097.      'SI' monsouc ;
  2098.      'MESS' 'souci: ' (vale 'SOUC'); 'FINSI';
  2099.                             resmul = resmul * 0.25;
  2100.   'SI' autaug;   iraug = vrai; 'FINSI';
  2101.      HPP_EPS = FAUX;
  2102.      pastest = vrai;
  2103.      instab = vrai;
  2104.      urg = vrai;
  2105.      augmult = augmult * 1.4;
  2106.     'SI' (iraug 'ET' (augmult < 1d3));
  2107.       ZDEPT = zdeptp ;
  2108.       zdep1 = (zdept - zdeptp) + (zdept 'EXCO' 'LX' 'NOID' 'LX');
  2109.       'ITERER' etiq;
  2110.      'FINSI';
  2111.     'FINS';
  2112. 'FINS';
  2113. *
  2114.   ZRAID_T = 'EXTR' ZRAID 'CONT';
  2115.   'SI' ('EXIS' ZRAID_T 'NITER');
  2116.        'SI' IDYN ; STAB12.'ZRAIDV'= ZRAID; 'FINS';
  2117.         STAB12.'RIBLO_M' = ZRAID_T.7;
  2118.         STAB12.'LISEA_M' = ZRAID_T.6;
  2119.         WTAB.'MAIL_BLOM' = 'EXTR' ZRAID_T.7 'MAIL' 'MULT' ;
  2120.   'FINS';
  2121. 'FINS';
  2122. RECA_K = FAUX;
  2123. *
  2124. * 1ere iteration ------------------------------------------------------
  2125. 'SI'   IPREM  ;
  2126.           ZDEPT  = 'COPIER' ZDEP1 ;
  2127.           ZDEPTM = ZDEPT ;
  2128. *         ZDEPf  = ZDEPT ;
  2129.           ZDELA  =   'COPIER' ZDEPT ;
  2130.           'SI' (WTAB . 'MODAL') ;
  2131.              'SI' ('EXIS' WTAB 'MODCONTA') ;
  2132.                 mamoco1 = 'EXTR' (ZDEPT 'ENLEVER' 'LX') 'MAIL' ;
  2133.              'FINS';
  2134.           'FINS' ;
  2135. *   iprem est faux: on est apres la premiere operation---------------
  2136. *
  2137. 'SINON';
  2138. *     zdept est l'increment de deplacement total avec les lagrangiens
  2139. *     de la solution complete
  2140.         XXX1   = ZDEPT 'ENLEVER' 'LX' ;
  2141. *       ZDEPf  = (ZDEPT + ZDEPTP) * 0.5;
  2142. *       on cumule les deplacements mais pas les lx
  2143. *       Du^(i)  = Du^(i-1) + du   ;   Dlx^(i) = 0 + Dlx
  2144.  
  2145.         ZDEPT  = XXX1 '+' ZDEP1 ;
  2146.        'DETR' XXX1 ;
  2147. 'FINS' ;
  2148. *
  2149. 'SI' (WTAB . 'MODAL') ;
  2150.      'SI' ('EXIS' WTAB 'MODCONTA') ;
  2151. * mettre les point materiels dans zdept
  2152.          zdeptu1 =  'REDU' zdept   mamoco1      ;
  2153.          ch_dco  =  'RECO' zdeptu1 ZMODL ZMATFI ;
  2154.          ZDEPT   = ('EXCO' zdept 'LX' 'LX') 'ET' zdeptu1 'ET' ch_dco ;
  2155.      'FINS';
  2156. 'FINS' ;
  2157. *
  2158. * ATTENTION ON EST SUR LA CONFIGURATION GEOM1?
  2159.  
  2160.  
  2161. *
  2162.  
  2163. *
  2164. *--- CAS 2 ------------------------------------------------------------*
  2165. *  si on part dans les decors on redemarre a 0
  2166. 'SI' ((XCONV > 1E8) 'ET' PASREINI 'ET' ('NON' IPILOT)) ;
  2167.         'MESS' 'Reinitialisation du schema';
  2168.       ZDEPT = zdeptp ;
  2169.       zdep1 = (zdept - zdeptp) + (zdept 'EXCO' 'LX' 'NOID' 'LX');
  2170.         PASREINI=FAUX;
  2171.         ITACC=3;
  2172.         'ITERER' ETIQ;
  2173. 'FINS';
  2174.  
  2175. *
  2176. *--- CAS 3 ------------------------------------------------------------*
  2177. 'SI'  ACCEL ;
  2178.  
  2179. *   pilotage automatique
  2180.     'SI' IPILOT;
  2181.               RED1 = STAB12.'AUTOREDU' ;
  2182.  
  2183. *       reduction du critere de pilotage red2 sert d'incateur si
  2184. *       pas accelere tous les 2 pas
  2185.  
  2186. * MODIFICATION CB215821 : 18/06/2015
  2187.            'SI' (RED2 < (IT '/' 20));
  2188.                 STAB12.'AUTOREDU' = STAB12.'AUTOREDU' '*' 3.D0 ;
  2189.                 RED2 = RED2 '+' 1 ;
  2190.                 ITACC= 3;
  2191.            'FINS';
  2192.  
  2193.             RED1 = RED1 '/' (STAB12.'AUTOREDU') ;
  2194.             'SI' (RED1 '<' 0.9) ;  STAB12.AUTORED1 = 10 ; 'FINS' ;
  2195.  
  2196.             'SI' ( STAB12.'AUTOREDU' '>' 1.D0 ) ;
  2197.                'MESS' 'On divise le critere de pilotage par '
  2198.                 STAB12.'AUTOREDU';
  2199.             'FINS' ;
  2200.  
  2201. *      test si snap back et si refus de l'acceleration
  2202.             'SI' IPLAVI ;
  2203.                   XXX2 = 'EXCO' ZDEPT  'LX' 'NOID' 'LX' ;
  2204.                   XXX1 = 2.D0  '*' XXX2;
  2205.                   XXX3 = ZDEPT '-' XXX1;
  2206.             'SINON' ;
  2207.                   XXX1 = COEPI '*' ZDELA;
  2208.                   XXX3 = ZDEPT '-' XXX1 ;
  2209.                   XXX2 = 'EXCO' XXX3 'LX' 'NOID' 'LX';
  2210.                   XXX1 = XXX2  '*' 2.D0 ;
  2211.                   XXX4 = COEPI '*' ZDELA;
  2212.                   XXX3 = ZDEPT '-' XXX4 '-' XXX1;
  2213.            'FINS' ;
  2214.  
  2215.            'DETR' XXX2 ;'DETR' XXX1; 'DETR' XXX3;
  2216.             XXX1 ='EXCO' ZDEPT 'LX' 'NOID' 'LX' ;
  2217.             XXX2 = XXX1 '*' 2.D0 ;
  2218.             XXX3 = ZDEPT '-' XXX2;
  2219.             SRT  = 'XTY' ZDFINI XXX3 MLDUAL MLPRIM;            pvpv      ;
  2220.            'DETR' XXX1 ;  'DETR' XXX2; 'DETR' XXX3;
  2221.  
  2222.             ISNPB = FAUX;
  2223.            'SI'  (SRT '<' 0) ;
  2224.                 ISNPB = VRAI ;
  2225.            'FINS';
  2226.  
  2227.             OO = WTAB.AUTOCRIT ;
  2228.             OO = OO '/' STAB12.'AUTOREDU' ;
  2229.             U1MA = AUTOPILO ZDEPT (COEPI'*'ZDELA) ZMODLI ZMATFI WTAB;
  2230.             AL1 = OO '/' U1MA ;
  2231. *        al1 est le coefficient de normalisation
  2232.            'SI' (( al1 '>EG' 1.d0) 'ET' (coepi '>EG' 1.d0)) ;
  2233. *     pour eviter d'aller au dela de alpha=1 on ignore le critere
  2234.                  AL1 = 1.d0 ;
  2235.            'FINS' ;
  2236.             'SI'((AL1 '>' 1.D0) 'ET' (COEPI '>' 0.D0));
  2237.               'SI' ( AL1 '>' (1.D0 '/' COEPI));
  2238.                    AL1 = 1.D0 '/' COEPI;
  2239.               'FINS';
  2240.             'FINS';
  2241. *  normalisation
  2242.  
  2243.               PASTEST=VRAI;
  2244.               XXX1   = ZLX   '*' ( 1.d0 '-' AL1);
  2245.               XXX3   = ZDEPT '*' AL1 ;
  2246.               ZDEPT   = XXX3  '+' XXX1;
  2247.              'DETR' XXX3;
  2248.     'FINS';
  2249. 'SINON';
  2250.       'SI' IPILOT ;
  2251. **
  2252. *   le calcul d'al1 est fait pour eviter la convergence
  2253. *   s'il est  externe a l'interval 0.5 --- 2.
  2254.             OO = WTAB.'AUTOCRIT' ;
  2255.             OO = OO '/' STAB12.'AUTOREDU' ;
  2256.             U1MA = AUTOPILO ZDEPT (COEPI'*'ZDELA) ZMODLI ZMATFI WTAB;
  2257.             AL1 = OO '/' U1MA ;
  2258.            'SI' (( AL1 '>EG' 1.d0) 'ET' (COEPI '>EG' 1.d0)) ;
  2259.                  AL1 = 1.d0 ;
  2260.            'FINS' ;
  2261.       'FINS' ;
  2262. 'FINS';
  2263. *
  2264. *  garder les reactions pour le test de convergence
  2265. *
  2266. ZDEPLP = ZDEPL                         ;
  2267. ZDEPL  = ZDEPT 'EXCO' 'LX' 'NOID' 'LX' ;
  2268.  
  2269. *----------------------------------------------------------------------
  2270. *     le nouveau champ est fixe on va tester l'equilibre(convergence)
  2271. *     et calculer la force motrice pour l'iteration suivante
  2272. *----------------------------------------------------------------------
  2273. *
  2274. *  calcul de  fcorf = lambda * m   force de reaction -----------------
  2275. *             fcoru = m * u     depimp (flx)
  2276.  
  2277. 'DETRUIRE' FCORF;
  2278. *       calcul du deplacement total u
  2279. zdetotp = zdetot;
  2280. ZDETOT = ZZD   '+' ZDEPT ;
  2281. XXX1   = ZCLIM '*' ZDETOT;
  2282. * -1*reactions a l'iteration it
  2283. FCORF = 'ENLEVER'  XXX1 'FLX';
  2284. * valeur des relations imposees a l'iteration it
  2285. FCORU = 'EXCO'  XXX1 'FLX'  'NOID' 'FLX';
  2286. 'DETR' XXX1 ;
  2287. *
  2288. *  dans le cas des modeles endommageables de Lemaitre, on ecoule
  2289. *  en tenant compte, dans les iterations internes, de la variation du
  2290. *  materiau avec la temperature
  2291. *
  2292. ZMATT = ZMAT    ;
  2293. 'SI' ('ET' ('ET' ITHER ('OU' IENDOM IVIDOM ) ) WTAB.'MATVAR')  ;
  2294. * on recupere certain materiau avec les parametres fct de la temperatue
  2295. * voir PAS_mate  (il ne faut que la dependance thermique)
  2296.       ZMATT = 'REDU' WTAB.'MA_COMP' ZMODL;
  2297. 'FINS';
  2298.  
  2299. 'SI' IFEFP;
  2300. * Update or total lagrangian -----------------------------------------
  2301.   'SI' IFEFPUL;
  2302. *  mess ' update lagrangian ZRIKTA';
  2303.     GEOM2 = 'FORM' ZDEPT ;
  2304.     ZRIKTA ZSIGF ZVARF ZDEIF = 'ECFEFP'
  2305.          ZMODL ZEPS0 ZVAR0 ZDEPT ZMAT ZPREK NITMA XUPDA;
  2306.   'SINON';
  2307. *  mess ' total lagrangian ZRIKTA';
  2308.     chp_z = ZDEPT + ZU1 ;
  2309.     GEOM2 = 'FORM' chp_z ;
  2310.     ZRIKTA ZSIGF ZVARF ZDEIF = 'ECFEFP'
  2311.          ZMODL ZEPS0 ZVAR0 chp_z ZMAT ZPREK NITMA ;
  2312.     chp_z = 1 ;
  2313.   'FINS';
  2314.   FEQU2 = 'BSIG' ZMODL ZSIGF ZMAT ;
  2315.   ZRAID = ZRIKTA 'ET' ZCLIM ;
  2316.    'FORM' GEOM1 ;
  2317. **
  2318.   XXX1 = 'EXCO' (WTAB.'MOVA') ZVARF ;
  2319.   EPSM = 'MAXI' XXX1 'AVEC' MLDEFOR ;
  2320.   ACC  = 'ABS' ( XXX1 - ACC0 ) ;
  2321. *'DETR' XXX1 ;
  2322.   MMC = 'MASQUE' ACC 'SUPERIEUR' 1.D-10 'SOMME' ;
  2323.   'SI' (MMC '>' MMCMAX) ; MMCMAX = MMC ; 'FINS' ;
  2324.   DPSMAX  = 'MAXI' ACC 'AVEC'  MLDEFOR ;
  2325.   DEPST   = 'CHAN' 'TYPE' ZDEIF 'DEFORMATIONS' ;
  2326.  
  2327.   GR_U_FIN= 'MOT' 'INCONNU';
  2328.  
  2329. 'SINON';
  2330. * cas standard --------------------------------------------------------
  2331.   ZMAT05=ZMAT;
  2332.   GEOM2 = GEOM1;
  2333.   'SI' IGRD;
  2334.     zdepm  = zdept '*' 0.5;
  2335.     'FORM' GEOM1 ;
  2336.      GEOM2m = 'FORM' zdepm  ;
  2337.      GEOM2  = 'FORM' zdepm  ;
  2338.     'FORM' GEOM1 ;
  2339.   'FINS';
  2340. *
  2341. *     calcul de l'increment de deformation elast totale  DEPST--------
  2342. *
  2343. * Pour les hypotheses de deformations non lineaires,
  2344. * on calcule l'increment de deformation sur le pas avec EPSI 'LINE'
  2345. * mais en se placant sur une configuration intermediaire (mi pas)
  2346.   'SI' (('NEG' HYPDEF 'LINEAIRE') 'ET' IGRD);
  2347.     'SOUC' 0;
  2348.     'SI' ITCAR;
  2349.       GEOMIL  ZMATTEMP = 'FORM' ZDEPM ZMODLI (ETG ZMAT05) ;
  2350. *     A cause de FORM, ZMATTEMP n'est plus parallele...
  2351.       ZMATTEMP         = 'REDU' ZMATTEMP ZMODL ;
  2352.       DEPST            = 'EPSI' 'LINE' 'NOER' ZMODL ZDEPT ZMATTEMP ;
  2353.       ZMATTEMP = 0 ;
  2354.     'SINON';
  2355.       GEOMIL           = 'FORM' ZDEPM    ;
  2356.       DEPST            = 'EPSI' 'LINE' 'NOER' ZMODL ZDEPT ZMAT05   ;
  2357.     'FINS';
  2358.     PASOK = 'SOUCI';
  2359.  
  2360.     'SI' ('NON' PASOK);
  2361.       'FORM'  geom1;
  2362.       DEPST = 'CAPI' HYPDEF DEPST ZMODL ZDEPM;
  2363.       'SI' ('MAXI' 'ABS' DEPST '>' 10);
  2364.         PASOK = VRAI;
  2365.       'FINS';
  2366.     'FINS';
  2367.  
  2368.     'SI' PASOK;
  2369.  mess 'annulation sous-increment: deformation non lineaire non calculable';
  2370.       itacc = 3;
  2371.       resmul = resmul * 0.25;
  2372.       'SI' autaug;   iraug = vrai; 'FINSI';
  2373.       HPP_EPS = FAUX;
  2374.       pastest = vrai;
  2375.       urg = vrai;
  2376.       itacc = 3;
  2377.       augmult = augmult * 1.4;
  2378.       ZDEPT = zdeptp ;
  2379.       zdep1 = (zdept - zdeptp) + (zdept 'EXCO' 'LX' 'NOID' 'LX');
  2380.       'SI' IGRD;
  2381.         'FORM' GEOM1 ;
  2382.       'FINS';
  2383.       'ITERER' etiq;
  2384.     'FINSI';
  2385.   'SINON';
  2386.     PASOK = vrai;
  2387.   'FINS';
  2388.  
  2389. ** on n'a pas pu calcule DEPST : on utilise les deformations lineaires
  2390.   'SI' PASOK;
  2391. *    'MESS' 'deformation lineaire';
  2392.     'SI' IGRD;
  2393.       'FORM' GEOM1;
  2394.     'FINSI';
  2395.     DEPST = 'EPSI' 'LINE' ZMODL ZDEPT ZMAT05;
  2396.   'FINS';
  2397.  
  2398. * calcul des increments de deformation et contrainte sur la configuration initiale
  2399. * en lagrangien total
  2400.   'SI' IGRD ;
  2401.     'SI' (LAG_TOT 'EGA' 1);
  2402. *  passage de depst et sig0 sur for0
  2403.       'FORM' GEOREF0;
  2404.       DEPST = 'CAPI' HYPDEF DEPST ZMODL ZU1 ;
  2405.       ZSIG0 = 'CAPI' HYPDEF ZSIG0 ZMODL ZU1 ;
  2406.     'FINSI';
  2407. *  passage de depst et sig0 sur geom2
  2408.     'SI' (LAG_TOT 'EGA' 2);
  2409.       DEPST = 'PICA' HYPDEF DEPST ZMODL ZDEPT ;
  2410.       ZSIG0 = 'PICA' HYPDEF ZSIG0 ZMODL ZDEPT ;
  2411.       'FORM'  geom2;
  2412.     'FINS';
  2413.     'SI' (LAG_TOT 'EGA' 3);
  2414.       DEPST = 'PICA' HYPDEF DEPST ZMODL ZDEPM ;
  2415.       ZSIG0 = 'PICA' HYPDEF ZSIG0 ZMODL ZDEPM ;
  2416.       'FORM'  geom2m;
  2417.     'FINS';
  2418.   'FINS';
  2419.  
  2420. *  ICI depst et zsig0 sont sur la configuration initiale en lagrangien total et
  2421. *      il faut soustraire les parties thermique et deformations imposees
  2422.   'SI' IPLAVI;
  2423.     'SI' ('OU' ('OU' IENDOM IVIDOM)  ICERAM);
  2424.       HOOKENDO = 'HOOK'          ZMODL ZMAT05 ZVAR0          ;
  2425.       'SI' IGRD;
  2426.          DEPST   = 'EPSI' 'QUAD' ZMODL ZDEPT  HOOKENDO ZMAT05;
  2427.       'SINON';
  2428.         DEPST    = 'EPSI' 'LINE' ZMODL ZDEPT  HOOKENDO ZMAT05;
  2429.       'FINS';
  2430.     'FINS';
  2431.  
  2432. *    si thermoplastique on enleve e alpha *dt et d'autres termes
  2433. *    si le materiau depend de la temperature
  2434.     'SI' (ITHER 'OU' WTAB.'MATVAR') ;
  2435.        XXX1  = DEPST0 '*' COEPI;
  2436.        XXX2  = DEPST  '-' XXX1 ;
  2437. *       'DETR' DEPST ; 'DETR' XXX1;
  2438.        DEPST = XXX2 ;
  2439.     'FINS' ;
  2440.  
  2441.     'SI' LOGDEF  ;
  2442.        XXX1 = DDEFOR0 '*' COEPI;
  2443.        XXX2 = DEPST   '-' XXX1 ;
  2444. *       'DETR' DEPST ; 'DETR' XXX1;
  2445.        DEPST = XXX2 ;
  2446.     'FINS' ;
  2447.  
  2448.   'SINON';
  2449. *   en elasticite, on a directement l'increment de contrainte
  2450.     DSIGT= 'ELAS'  ZMODL DEPST ZMAT05 ;
  2451.     'SI' (('NON' ITHER) 'ET' WTAB.'MATVAR') ;
  2452.       XXXXX3 = 'ELAS' ZMODL ZSIG0 ZMAT05 ;
  2453.       XXXXX4 = 'ELAS' ZMODL ZSIG0 ZMAT1  ;
  2454.       XXDEFO = XXXXX4 '-' XXXXX3 ;
  2455.       DSIGTV = 'ELAS' ZMODL XXDEFO ZMAT05 ;
  2456.       DSIGT  =  DSIGT '+' DSIGTV ;
  2457. *      'DETR' XXXXX3 ; 'DETR' XXXXX4 ;
  2458. *      'DETR' XXDEFO ; 'DETR' DSIGTV ;
  2459.     'FINS' ;
  2460. *   si thermoplastique on enleve alpha *dT et d'autres termes
  2461. *   si le materiau depend de la temperature
  2462.     'SI' ITHER  ;
  2463.       XXX1 = DSIGT0 '*' COEPI;
  2464.       XXX2 = DSIGT  '-' XXX1;
  2465.       'DETR' DSIGT ; 'DETR' XXX1;
  2466.       DSIGT =  XXX2  ;
  2467.     'FINS';
  2468.  
  2469.     'SI' LOGDEF  ;
  2470.       XXX1 = DSI1  '*' COEPI;
  2471.       XXX2 = DSIGT '-' XXX1;
  2472. *      'DETR' DSIGT ; 'DETR' XXX1;
  2473.       DSIGT =  XXX2 ;
  2474.     'FINS';
  2475.   'FINS';
  2476. *
  2477. *  Calcul du gradient des deplacements (total) a la fin du pas de ------
  2478. *  temps par rapport a la configuration de reference GEOREF0 (initiale)
  2479. *  rappel : On a : ZDETOT = ZZD + ZDEPT ;
  2480. 'SI' igrd;
  2481. 'SI' ('NEG' GR_U_DEB 'INCONNU');
  2482.   'SI' (LAG_TOT 'NEG' 1); geot = 'FORM';
  2483.      'FORM' WTAB.FOR0; 'FINS' ;
  2484.    GR_U_FIN = 'GRAD' ZMODL ZMAT ZDETOT ;
  2485.    D_GR_U   =  GR_U_FIN '-' GR_U_DEB   ;
  2486.   'SI' (LAG_TOT 'NEG' 1); 'FORM' GEOT; 'FINS';
  2487. 'SINON';
  2488.    GR_U_FIN = 'MOT' 'INCONNU';
  2489. 'FINS';
  2490. 'FINS';
  2491. *
  2492.   'SI' ZCCONV ;
  2493.     STAB12.'DT' = DTINI '*' COEPI ;
  2494.   'SINON' ;
  2495.     STAB12.'DT' = 0. ;
  2496.   'FINS' ;
  2497. *
  2498. *  si on est plastique ou viscoplastique on ecoule pour avoir-----------
  2499. *  le nouveau champ de contraintes
  2500. *
  2501. MMC = 0 ;
  2502. 'SI'  IPLAVI ;
  2503.  
  2504. *...cas SSTE...
  2505.    'SI' ISSTE;
  2506.       ZRIKTA ZSIGF ZVARF ZDEIF =
  2507.                         'SSTE' ZMODL ZSIG0 ZVAR0 DEPST ZMATT
  2508.                                ZPREK NSSTE NITMA;
  2509.       zvar0 = ('EXCO' zvar0 ('ENLE' lnom ('DIME' lnom))) 'ET'
  2510.               ('EXCO' zvarf ssii);
  2511.       zvar0 = 'CHANGER' 'TYPE' zvar0 'VARIABLES INTERNES';
  2512.  
  2513. *...cas non SSTE...
  2514.    'SINON';
  2515.       ZSIG01=ZSIG0;
  2516.       ZVAR01=ZVAR0;
  2517.       ZEPS01=ZEPS0;
  2518. **     defineto = DEPIN0 ; TABCONT='TABLE';
  2519.  
  2520.      'SI' (nsoincr 'NEG' 1);
  2521.         ZSIGM=ZSIG0 '/' 2;
  2522.      'FINS';
  2523.  
  2524.      'REPETER'  sousinc nsoincr;
  2525.  
  2526.        tsodeb='FLOT' stab12 .'DT'/nsoincr*( &sousinc - 1)+conti.'TEMPS';
  2527.        tsofin='FLOT' stab12 .'DT'/nsoincr*  &sousinc     +conti.'TEMPS';
  2528.       'SI' ('NON' ZCCONV) ;
  2529.          tsodeb = tsofin ;
  2530.       'FINS' ;
  2531.  
  2532.        che1= 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'TEMP' tsodeb;
  2533.        che2= 'ADET' 'NOUV' ZMODL CHASANST tsofin 'TEMP' tsofin;
  2534.  
  2535.        'SI' ('EXIS' WTAB 'MOD_LIA') ;
  2536.           vite1 =  'CHAN' 'COMP' conti.'VITESSES' MLPRIM MVPRIM ;
  2537.           che1  = che1
  2538.             'ET' ('CHAN' 'CHAM' ZMODL vite1                'STRESSES')
  2539.             'ET' ('CHAN' 'CHAM' ZMODL conti.'DEPLACEMENTS' 'STRESSES');
  2540.           che2  = che2
  2541.             'ET' ('CHAN' 'CHAM' ZMODL  ZDETOT              'STRESSES')
  2542.             'ET' ('CHAN' 'CHAM' ZMODL  ZFCONST1            'STRESSES');
  2543.        'FINS';
  2544.  
  2545.        'SI' ither ;
  2546.          TETDE=STAB12.'TET1' + (DTETD *('FLOT' (&sousinc - 1)/nsoincr));
  2547.          TETDF=STAB12.'TET1' + (DTETD *('FLOT'  &sousinc     /nsoincr));
  2548.          che3 = 'CHAN' 'CHAM' ZMODL TETDE 'STRESSES' ;
  2549.          che4 = 'CHAN' 'CHAM' ZMODL TETDF 'STRESSES' ;
  2550.        'FINS' ;
  2551.  
  2552.        che1 = che1 'ET' che3 ;
  2553.        che2 = che2 'ET' che4 ;
  2554.        aa5  = DEPST '*' ( (&sousinc '-' 1.) '/' nsoincr);
  2555.        che5 = DEFT0 '+'   aa5;
  2556.        cc5  = DEPST '*' ('FLOT' &sousinc '/' nsoincr);
  2557.        che6 = DEFT0 '+'  CC5;
  2558. *On met les deformations en tete des champs pour COMP ('KTAN' 'PERT')
  2559.        che1 = che5 'ET' che1 ;
  2560.        che2 = che6 'ET' che2 ;
  2561. * pour les materiaux on garde toujours la valeur a la fin du pas
  2562. * et au debut du pas
  2563. *       'SI' ('NEG'  ZMAT1 'INCONNU') ;
  2564.           che1 = che1 'ET' ZMAT1 ;
  2565. *       'FINS' ;
  2566. *
  2567. *  Modele NON LINEAIRE UTILISATEUR + GRANDES DEFORMATIONS
  2568. *  On ajoute les gradients de deplacements qui seront transformes en
  2569. *  gradients de transformation avant appel a UMAT (dans WKUMA1.ESO)
  2570.        'SI' (( 'NEG' GR_U_DEB 'INCONNU') 'ET' igrd) ;
  2571.          gru1 = GR_U_DEB + ( D_GR_U * ('FLOT' (&sousinc-1) / nsoincr));
  2572.          che1 = che1 'ET' gru1 ;
  2573.          gru2 = GR_U_DEB + ( D_GR_U * ('FLOT'  &sousinc    / nsoincr));
  2574.          che2 = che2 'ET' gru2 ;
  2575.        'FINS';
  2576.  
  2577.        'SI' LNLOC ;
  2578.           CHE1A = che1 ;
  2579.           CHE2A = che2 ;
  2580.  
  2581.           ISTEP = 1 ;
  2582.           chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'STEP' ISTEP ;
  2583.           che11 = CHE1A 'ET' chm_z ;
  2584.           che11 = che11 'ET' ZSIG01 'ET' ZVAR01 'ET' ZEPS01 ;
  2585.  
  2586.           chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsofin 'STEP' ISTEP ;
  2587.           che22 = CHE2A 'ET' chm_z ;
  2588.           che22 = che22 'ET' ZMATT ;
  2589.  
  2590.           'SI' PARTLOCA ;
  2591.             cha11 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che11 ;
  2592.             cha22 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che22 ;
  2593.             souci 0 ;
  2594.             cho22 = 'ASSIS' 'TOUS' 'COMP' MODRELOC cha11 cha22 ;
  2595.             isoucomp = souci ;
  2596.             cho2  = 'REDU' ZMODL ('ETG' cho22) ;
  2597.             cha11 = 0 ; cha22 = 0 ; cho22 = 0 ;
  2598.           'SINON';
  2599.              souci 0 ;
  2600.              cho2  = 'COMP' ZMODL che11 che22 ;
  2601.              isoucomp = souci ;
  2602.           'FINS' ;
  2603.  
  2604.           ZVARF = 'CHANGER' ('EXCO' cho2 com_var 'NOID')
  2605.                             'TYPE' 'VARIABLES INTERNES' ;
  2606. * NLOC ne traitant pas des champs //, on reduit tout au modele initial.
  2607.           'SI' (PARALLEL 'ET' ('NON' PARTLOCA)) ;
  2608.             ZVARF = 'REDU'  ZVARF ZMODLI ;
  2609.           'FINS' ;
  2610. *   cas non-local MOYE ou SB
  2611.           'SI' ('NEG'  WTAB.'NON_LOCAL' 'HELM');
  2612.             'SI' ('EGA' WTAB.'NON_LOCAL' 'SB') ;
  2613.                MOD_SB =  WTAB.'NLOC_MODL' ;
  2614.                CONTP = 'PRIN' ('REDU' ZSIG01 MOD_SB) MOD_SB ;
  2615.                ZVARF = ZVARF '+' CONTP '+' WTAB.'NLOC_SB_REGU' ;
  2616.             'FINS' ;
  2617.             ZVARN = 'NLOC'  ZVARF  WTAB.'CONN'  ;
  2618.             ZVARN = 'CHANGER' ('EXCO' ZVARN com_var)
  2619.                           'TYPE' 'VARIABLES INTERNES' ;
  2620.           'SINON' ;
  2621. *   cas non-local HELM
  2622.             MOD_HELM =  WTAB.'NLOC_MODL' ;
  2623.             ZVARN = 'COPI' ZVARF ;
  2624.             ZVARF1 = 'REDU' ZVARN MOD_HELM ;
  2625.             LICOHELM = 'EXTR' ZVARF1 'COMP' ;
  2626.             'REPE' BH NHELM ;
  2627.               LEMO = TAHELM . &BH. 'NOM' ;
  2628.               MOPRE = 'EXTR' LEMO 1 3 ;
  2629.               ZVAUX = 'EXCO' ZVARF1 LEMO 'SCAL';
  2630.               ZVAUX2 = 'CHAN' 'NOEUD' MOD_HELM ZVAUX;
  2631.               ZVAUX2 = 'PROI' TAHELM . &BH . 'H_MODELE'  ZVAUX2 'MINI';
  2632.               FSOUR = 'SOUR' TAHELM . &BH . 'H_MODELE' ZVAUX2 ;
  2633.               ZVNEW = 'RESO' TAHELM . &BH . 'H_OPER' FSOUR ;
  2634.               ZVNEW = 'NOMC' ZVNEW 'SCAL' ;
  2635.               ZVNEW = 'CHAN' 'CHAM' ZVNEW MOD_HELM 'STRESSES'
  2636.                                              'VARIABLES INTERNES';
  2637.               DZVN = ZVNEW '-' ZVAUX ;
  2638.               DZVN = 'NOMC' DZVN LEMO ;
  2639.               ZVARN = ZVARN '+' DZVN ;
  2640.  
  2641. *             cas particulier SMAX
  2642.               'SI' (('EXIS' LICOHELM 'SMAX') 'ET' ('NEG' MOPRE 'ERF'));
  2643.                 ZSMAX = 'EXCO' ZVARF1 'SMAX' 'SMAX' ;
  2644.                 ZSMAN = 'MANU' 'CHML' MOD_HELM 'SMAX' ('MAXI' ZSMAX)
  2645.                     'TYPE' 'VARIABLES INTERNES' 'STRESSES' ;
  2646.                 ZVARN = ZVARN '+' ZSMAN '-' ZSMAX ;
  2647.               'FINS' ;
  2648.             'FIN' BH ;
  2649.           'FINS' ;
  2650. * On rend ZVARN // si besoin :
  2651.           'SI' (PARALLEL 'ET' ('NON' PARTLOCA)) ;
  2652.             ZVARN = 'REDU' ZVARN ZMODL ;
  2653.           'FINS' ;
  2654.  
  2655.           ISTEP = 2 ;
  2656.           chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'STEP' ISTEP ;
  2657.           che11 = CHE1A 'ET' chm_z ;
  2658.           che11 = che11 'ET' ZSIG01 'ET' ZVARN 'ET' ZEPS01 ;
  2659.  
  2660.           chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsofin 'STEP' ISTEP ;
  2661.           che22 = CHE2A 'ET' chm_z ;
  2662.           che22 = che22 'ET' ZMATT ;
  2663.  
  2664.           'SI' PARTLOCA ;
  2665.             cha11 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che11 ;
  2666.             cha22 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che22 ;
  2667.             souci 0 ;
  2668.             cho22 = 'ASSIS' 'TOUS' 'COMP' MODRELOC cha11 cha22 ;
  2669.             isoucomp = souci ;
  2670.             cho2  = 'REDU' ZMODL ('ETG' cho22) ;
  2671.             cha11 = 0 ; cha22 = 0 ; cho22 = 0 ;
  2672.           'SINON';
  2673.             souci 0 ;
  2674.             cho2  = 'COMP' ZMODL che11 che22 ;
  2675.             isoucomp = souci ;
  2676.           'FINS' ;
  2677.  
  2678. * Un peu de menage
  2679.           CHE1A = 1 ;
  2680.           CHE2A = 1 ;
  2681.           ZVARN = 1 ;
  2682.           ZVARF = 1 ;
  2683.           chm_z = 1 ;
  2684.  
  2685.        'SINON' ;
  2686.            che11 = che1 'ET' ZSIG01 'ET' ZVAR01 'ET' ZEPS01 ;
  2687.            che22 = che2 'ET' ZMATT ;
  2688. *           si ( 'EXIST'  PRECED 'ECRIT' ) ;
  2689. *           'FINS';
  2690.            'SI' PARTLOCA;
  2691.                cha11 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che11 ;
  2692.                cha22 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che22 ;
  2693.                souci 0 ;
  2694.                cho22 = 'ASSIS' 'TOUS' 'COMP' MODRELOC cha11 cha22 ;
  2695.                isoucomp = souci ;
  2696.                cho2  = 'REDU' ZMODL ('ETG' cho22) ;
  2697.                cha11 = 0 ; cha22 = 0 ; cho22 = 0 ;
  2698.            'SINON';
  2699. *          'SI' ( 'EXIS' PRECED 'ECRIT' );
  2700. *          'FINS';
  2701.               souci 0 ;
  2702.               cho2  = 'COMP' ZMODL che11 che22 ;
  2703.               isoucomp = souci ;
  2704.            'FINS';
  2705.        'FINS' ;
  2706.        si isoucomp ;
  2707.          mess ' Souci dans le comportement' ;
  2708.        fins ;
  2709.    ZSIGF = 'CHANGER' ('EXCO' cho2 com_sig 'NOID')'TYPE' 'CONTRAINTES' ;
  2710.    ZVARF = 'CHANGER' ('EXCO' cho2 com_var 'NOID')'TYPE' 'VARIABLES INTERNES' ;
  2711.    ZDEIF = 'CHANGER' ('EXCO' cho2 com_dei 'NOID')'TYPE' 'DEFORMATIONS INELASTIQUES' ;
  2712.  
  2713.        'SI' ('EXIS' WTAB 'MOD_LIA') ;
  2714.          ZFLIA = 'EXCO' cho2 'FLIA' 'NOID';
  2715.        'FINS' ;
  2716.        'SI' (nsoincr 'NEG' 1);
  2717.          'SI' (&sousinc 'EGA' nsoincr);
  2718.            aaa1=zsigf/2;
  2719.            aaa2=aaa1+ zsigm;
  2720.            'DETR' aaa1; 'DETR' zsigm;
  2721.            zsigm=aaa2;
  2722.          'SINON';
  2723.            aaa2=zsigf+ zsigm;
  2724.            'DETR' zsigm;
  2725.            zsigm=aaa2;
  2726.          'FINS';
  2727.        'FINS';
  2728.        zsig01 = zsigf ;
  2729.        ZVAR01 = ZVARF ;
  2730.        ZEPS01 = ZDEIF ;
  2731. **       defineTO=ZDEIF+defineTO;
  2732.      'FIN' sousinc;
  2733. *
  2734. * Matrice tangente par perturbation evaluee pour la derniere iteration calculee
  2735. * A voir : cas grand deplacement ZSIGF et ZMAT, cas poreux et thermique ZSIGF
  2736.      'SI' (IKTAN 'ET' IPERT) ;
  2737.        Z1COMP = che11 ;
  2738.        Z2COMP = che22 'ET' ZSIGF ;
  2739.      'FINS' ;
  2740. *
  2741.      'SI' (nsoincr 'NEG' 1);
  2742. * pour tenir compte de ce que le travail de la correction est 1/2 FU,
  2743. * on la multiplie par 2
  2744.           CONT  =(ZSIGF + ZSIG0)          ;
  2745.           ZSIGM = ZSIGM*(2. /nsoincr)     ;
  2746.           ZSIGM = ZSIGM - CONT            ;
  2747.           BZSIGM='BSIG' ZMODL ZSIGM ZMAT;
  2748. **        ZDEPSPL=defineTO - ZEPS0 ;
  2749.      'FINS';
  2750. * on enleve toutes traces de champs inutiles
  2751.      zsig01=1;ZVAR01=1;ZEPS01=1;
  2752. **     defineTO=1; tabcont =1;
  2753.  
  2754. *     le fin d'en dessous est le fin de "si isste ... sinon ... finsi"
  2755. *     on est encore dans "si iplavi ... finsi"
  2756.    'FINS';
  2757. *...fin du si ISSTE sinon ...
  2758. *
  2759. *        cas particulier poreux et thermique
  2760. *
  2761.    'SI' ITHER;
  2762.      'SI' POR1 ;
  2763.         ZSIGF = ZSIGF - ( COEPI * DMSRT0 ) ;
  2764.      'FINS';
  2765.    'FINS';
  2766. *
  2767. *  max de epse pendant l'iteration
  2768. *  nbre de points qui ont une evolution non lineaire
  2769. *
  2770.    'SI' ('EGA' WTAB.'MOVA' 'RIEN') ;
  2771. **    EPSM = 0.; DPSMAX=0.;
  2772.       MMC=0;
  2773.    'SINON';
  2774.       XXX1 = 'EXCO' (WTAB.'MOVA') ZVARF ;
  2775. **    EPSM = 'MAXI' XXX1 ;
  2776.       ACC = 'ABS' ( XXX1 - ACC0 ) ;
  2777. **    'DETR' XXX1 ;
  2778.       MMC = 'MASQUE' ACC 'SUPERIEUR' 1.D-10 'SOMME' ;
  2779.       'SI' (MMC > MMCMAX) ; MMCMAX = MMC ; 'FINS' ;
  2780. **    DPSMAX = 'MAXI' ACC ;
  2781.    'FINS';
  2782. *
  2783. *  dans le cas SANS  plasticite,viscoplas,endommagement
  2784. *  on calcule directement le nouveau champ de contrainte
  2785. *
  2786. 'SINON' ;
  2787.     ZSIGF = ZSIG0 + DSIGT ;
  2788. 'FINS';
  2789. * si il y a lieu transport sur configuration debut de pas
  2790.     'SI' IGRD ;
  2791. **      ZSIG0 = ZSIG0_1;
  2792.       'SI' (LAG_TOT 'EGA' 1);
  2793. **        ZSIG0 =  'PICA' HYPDEF ZSIG0 ZMODL ZU1;
  2794.         ZSIGF =  'PICA' HYPDEF ZSIGF ZMODL ZU1;
  2795.         DEPST =  'PICA' HYPDEF DEPST ZMODL ZU1;
  2796.         'FORM' GEOM1;
  2797.       'FINS';
  2798.       'SI' (LAG_TOT 'EGA' 2);
  2799.         'FORM' GEOM1;
  2800. **        ZSIG0 =  'CAPI' HYPDEF ZSIG0 ZMODL ZDEPT;
  2801.         ZSIGF =  'CAPI' HYPDEF ZSIGF ZMODL ZDEPT;
  2802.         DEPST =  'CAPI' HYPDEF DEPST ZMODL ZDEPT;
  2803.       'FINS';
  2804.       'SI' (LAG_TOT 'EGA' 3);
  2805.         'FORM' GEOM1;
  2806. **        ZSIG0 =  'CAPI' HYPDEF ZSIG0 ZMODL ZDEPM;
  2807.         ZSIGF =  'CAPI' HYPDEF ZSIGF ZMODL ZDEPM;
  2808.         DEPST =  'CAPI' HYPDEF DEPST ZMODL ZDEPM;
  2809.       'FINS';
  2810.     'FINS';
  2811. *
  2812.  DPSMAX = 'MAXI' 'ABS' DEPST 'AVEC' MLDEFOR ;
  2813.     DEPSTI = DEPST;
  2814. 'SI' (IGRD 'ET' ('NON' HPP_EPS) 'ET' ('EGA' LAG_TOT 1));
  2815.     'FORM' georef0;
  2816.      DEPSTI = 'CAPI'  ZU1 HYPDEF DEPST ZMODL;
  2817. 'FINSI';
  2818. EPSMX   = DEFT0 '+' DEPSTI;
  2819. EPSM   ='MAXI' 'ABS' EPSMX 'AVEC' MLDEFOR ;
  2820. 'DETR' EPSMX;
  2821. 'SI' (IGRD 'ET' ('NON' HPP_EPS) 'ET' ('EGA' LAG_TOT 1));
  2822.      'FORM' geom1; 
  2823.      'DETR' DEPSTI;
  2824. 'FINSI';
  2825. *
  2826. *--- CAS 4 -------------------------------------------------------------
  2827. *  on  determine les forces resultantes de l'increment------------------
  2828. *
  2829. *   en grand deplacements on transforme pi en sigma
  2830. *   on met a jour les coordonnees
  2831. *
  2832.  ZMAT2=ZMAT;
  2833.  
  2834. 'SI' (IGRD 'ET' (('NON' HPP_EPS) 'OU' (IT > 0 )));
  2835.  
  2836.   ref = 50;
  2837.   coefmul = 1.d0 / (DPSMAX + XPETIT);
  2838.  
  2839.   'SI' (((coefmul   < ref) 'ET' AUTAUG) 'OU' PASUNIL);
  2840. *   il faut reduire d'urgence les deplacements ;
  2841.     zdeptlx = zdept 'EXCO' 'LX' 'NOID' 'LX';
  2842.     coefmul = 0.0;
  2843.     'SI' ((red_urg > 3) 'ET' ('NON' instab)); resmul = resmul * 0.25;
  2844.            'FINSI';
  2845.  
  2846.       augmult = augmult * 1.5;
  2847.       urg = vrai;
  2848.      'SI' autaug;    iraug = vrai; 'FINSI';
  2849. * on garde zdep1d pour le calcul des coeff d'acceleration
  2850.        zdep1d = zdep1 'ENLE' 'LX';
  2851.      'SI' AUTAUG;
  2852. *  reprise du vieux zdept. Tout le reste est recalcule a chaque pas de etiq
  2853.       ZDEPT = zdeptp ;
  2854.       zdep1 = (zdept - zdeptp) + (zdept 'EXCO' 'LX' 'NOID' 'LX');
  2855.      itacc = 3;
  2856.      'FINSI';
  2857.      'SI' ((dpsmaxp < 5e-4) 'OU' (dpsmaxp > 2e-2 )) ;
  2858.          PASTEST=VRAI;
  2859.      'SINON';
  2860.      'SI' (red_urg > 4);   nconvr = vrai; 'FINSI';
  2861.          tabconv . it = 1;
  2862.      'FINSI';
  2863.     'SI' (RED_URG > 2 )  ;
  2864.      'SI' ('NON' nonconv);  '