Télécharger la plaquette Cast3M


fr en

Télécharger unpas.procedur

Retour à la liste

Numérotation des lignes : 

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