Télécharger la plaquette Cast3M


fr en

Télécharger unpas.procedur

Retour à la liste

Numérotation des lignes : 

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