Procédures | Cast3M

* MRCFRAM2  PROCEDUR  BP208322  20/05/05    21:15:02     10601          
DEBPROC MRCFRAM2 TYP_CAL2*'MOT     '
   TAB_IN2*'TABLE' TOL_2*'FLOTTANT';
*
* --------------------------------------------------------
*  Procedure FRCFRAM2
*    Procedure appele par MRCFRAME pour le calcul de la marge sismique 
*    d'elements frame type POUTEAU LONG
*
*  Developpers:
*    Alberto FRAU
*    DEN/DANS/DM2S/SEMT/EMSI
*    Nicolas ILE
*    DEN/DANS/DM2S/SEMT/EMSI
*    Giulia DI STEFANO
*
*
*  LAM1 = MRCFRAM2 MOT1 TAB1 TOL1;
*
*
*  INPUT:
*
*    MOT1 [MOT] afin de specifier le type de calcul:
*              'ELLIPSOIDE' calcul des marges avec l'enveloppe 
*                           ellipsiodal
*              'RECTANGLE' calcul des marges avec l'enveloppe 
*                           prismatique
*
*    TAB1 type TABLE pour definir les donnes suivantes:
*         TAB1.'BASE_Y1'   dimension section direction Y en m [FLOTTANT]
*         TAB1.'BASE_Z1'   dimension section direction Z en m [FLOTTANT]
*         TAB1.'LONG_Y1'   longuer ly [FLOTTANT]
*         TAB1.'LONG_Z1'   longuer lz [FLOTTANT]
*         TAB1.'ENROBA1'   enrobage en m [FLOTTANT]
*         TAB1.'ACADRY1'   aire de l'armature transversale 
*                                   direction y en m2 [FLOTTANT]
*         TAB1.'ACADRZ1'   aire de l'armature transversale 
*                                   direction z en m2 [FLOTTANT]
*         TAB1.'S_CADR1'   espacement des cadres en 
*                                                    m  [FLOTTANT]
*         TAB1.'LI_PHI1'   LISTREEL des diam barres d'acier en mm
*         TAB1.'PY_PHI1'   LISTREEL de la position y en m de 
*                                   chaque barre par rapport au 
*                                   baricentre de la section
*         TAB1.'PZ_PHI1'  LISTREEL de la position z en m de 
*                                   chaque barre par rapport au 
*                                   baricentre de la section 
*         TAB1.'E_ACIE1'   Module Young acier [FLOTTANT]
*         TAB1.'EP_BET1'   deformation ultime beton [FLOTTANT]
*         TAB1.'EP_ACI1'   deformation ultime acier [FLOTTANT]
*         TAB1.'OM_SY1'    omega_s_y [FLOTTANT]
*         TAB1.'OM_SZ1'    omega_s_z [FLOTTANT]
*         TAB1.'OM_WZ1'    omega_w_y [FLOTTANT]
*         TAB1.'OM_WY1'    omega_w_z [FLOTTANT]
*         TAB1.'LAMB_Y1'   lambda_y [FLOTTANT]
*         TAB1.'LAMB_Z1'   lambda_y [FLOTTANT]
*         TAB1.'CHI_Z1'    chi_z [FLOTTANT]
*         TAB1.'CHI_Y1'    chi_y [FLOTTANT]
*         TAB1.'EFFX0'     EFFX statique [FLOTTANT]
*         TAB1.'EFFY0'     EFFY statique [FLOTTANT]
*         TAB1.'EFFZ0'     EFFZ statique [FLOTTANT]
*         TAB1.'MOMX0'     MOMX statique [FLOTTANT]
*         TAB1.'MOMY0'     MOMY statique [FLOTTANT]
*         TAB1.'MOMZ0'     MOMZ statique [FLOTTANT]
*         TAB1.'NNNN1'     (EFFX)2 - matrice X seisme [FLOTTANT]
*         TAB1.'NNVY1'     Cov EFFX et EFFY -matrice X seisme [FLOTTANT]
*         TAB1.'NNVZ1'     Cov EFFX et EFFZ -matrice X seisme [FLOTTANT]
*         TAB1.'NNMT1'     Cov EFFX et MOMX -matrice X seisme [FLOTTANT]
*         TAB1.'NNMY1'     Cov EFFX et MOMY -matrice X seisme [FLOTTANT]
*         TAB1.'NNMZ1'     Cov EFFX et MOMZ -matrice X seisme [FLOTTANT]
*         TAB1.'VYVY1'     (EFFY)2 - matrice X seisme [FLOTTANT]
*         TAB1.'VYVZ1'     Cov EFFY et EFFZ -matrice X seisme [FLOTTANT]
*         TAB1.'VYMT1'     Cov EFFY et MOMX -matrice X seisme [FLOTTANT]
*         TAB1.'VYMY1'     Cov EFFY et MOMY -matrice X seisme [FLOTTANT]
*         TAB1.'VYMZ1'     Cov EFFY et MOMZ -matrice X seisme [FLOTTANT]
*         TAB1.'VZVZ1'     (EFFZ)2 - matrice X seisme [FLOTTANT]
*         TAB1.'VZMT1'     Cov EFFZ et MOMX -matrice X seisme [FLOTTANT]
*         TAB1.'VZMY1'     Cov EFFZ et MOMY -matrice X seisme [FLOTTANT]
*         TAB1.'VZMZ1'     Cov EFFZ et MOMZ -matrice X seisme [FLOTTANT]
*         TAB1.'MTMT1'     (MOMX)2 - matrice X seisme [FLOTTANT]
*         TAB1.'MTMY1'     Cov MOMX et MOMY -matrice X seisme [FLOTTANT]
*         TAB1.'MTMZ1'     Cov MOMX et MOMZ -matrice X seisme [FLOTTANT]
*         TAB1.'MYMY1'     (MOMY)2 - matrice X seisme [FLOTTANT]
*         TAB1.'MYMZ1'     Cov MOMY et MOMZ -matrice X seisme [FLOTTANT]
*         TAB1.'MZMZ1'     (MOMZ)2 - matrice X seisme [FLOTTANT]
*         TAB1.'FCD1'      fcd du beton [FLOTTANT]
*         TAB1.'FSD1'      fsd du acier [FLOTTANT]
*
*    TOL1     Tolerance du calcul (1.e-4)
*
*
*  OUTPUT:
*
*
*    LAM1     Marge sismique
*
*
*                   
* --------------------------------------------------------
*
*
* Input Phase
BASE_Y1 = TAB_IN2.'BASE_Y1';
BASE_Z1 = TAB_IN2.'BASE_Z1';
LONG_Y1 = TAB_IN2.'LONG_Y1';
LONG_Z1 = TAB_IN2.'LONG_Z1';
ENROBA1 = TAB_IN2.'ENROBA1';
ACADRY1 = TAB_IN2.'ACADRY1';
ACADRZ1 = TAB_IN2.'ACADRZ1';
S_CADR1 = TAB_IN2.'S_CADR1';
LI_PHI1 = TAB_IN2.'LI_PHI1';
PY_PHI1 = TAB_IN2.'PY_PHI1';
PZ_PHI1 = TAB_IN2.'PZ_PHI1';
E_ACIE1 = TAB_IN2.'E_ACIE1';
EP_BET1 = TAB_IN2.'EP_BET1';
EP_ACI1 = TAB_IN2.'EP_ACI1';
OM_SY1  = TAB_IN2.'OM_SY1';
OM_SZ1  = TAB_IN2.'OM_SZ1';
OM_WZ1  = TAB_IN2.'OM_WZ1';
OM_WY1  = TAB_IN2.'OM_WY1';
LAMB_Y1 = TAB_IN2.'LAMB_Y1';
LAMB_Z1 = TAB_IN2.'LAMB_Z1';
CHI_Z1  = TAB_IN2.'CHI_Z1';
CHI_Y1  = TAB_IN2.'CHI_Y1';
EFFX0   = TAB_IN2.'EFFX0';
EFFY0   = TAB_IN2.'EFFY0';
EFFZ0   = TAB_IN2.'EFFZ0';
MOMX0   = TAB_IN2.'MOMX0';
MOMY0   = TAB_IN2.'MOMY0';
MOMZ0   = TAB_IN2.'MOMZ0';
NNNN1   = TAB_IN2.'NNNN1';
NNVY1   = TAB_IN2.'NNVY1';
NNVZ1   = TAB_IN2.'NNVZ1';
NNMT1   = TAB_IN2.'NNMT1';
NNMY1   = TAB_IN2.'NNMY1';
NNMZ1   = TAB_IN2.'NNMZ1';
VYVY1   = TAB_IN2.'VYVY1';
VYVZ1   = TAB_IN2.'VYVZ1';
VYMT1   = TAB_IN2.'VYMT1';
VYMY1   = TAB_IN2.'VYMY1';
VYMZ1   = TAB_IN2.'VYMZ1';
VZVZ1   = TAB_IN2.'VZVZ1';
VZMT1   = TAB_IN2.'VZMT1';
VZMY1   = TAB_IN2.'VZMY1';
VZMZ1   = TAB_IN2.'VZMZ1';
MTMT1   = TAB_IN2.'MTMT1';
MTMY1   = TAB_IN2.'MTMY1';
MTMZ1   = TAB_IN2.'MTMZ1';
MYMY1   = TAB_IN2.'MYMY1';
MYMZ1   = TAB_IN2.'MYMZ1';
MZMZ1   = TAB_IN2.'MZMZ1';
FCD1    = TAB_IN2.'FCD1';
FSD1    = TAB_IN2.'FSD1';
*
ENR1 = ENROBA1;
*
NLIMC1 = (-1.0 - OM_SY1)*(((BASE_Y1)*(BASE_Z1))*(FCD1));
NLIMT1 = (OM_SY1)*(((BASE_Y1)*(BASE_Z1))*(FCD1));
*
*  Static Point 
TINI1                     = TABLE;
TINI1.'TYPE'              = CHAINE 'P_LONG';
TINI1.'NN'                = EFFX0;
TINI1.'VY'                = EFFY0;
TINI1.'VZ'                = EFFZ0;
TINI1.'MT'                = MOMX0;
TINI1.'MY'                = MOMY0;
TINI1.'MZ'                = MOMZ0;
TINI1.'BY'                = BASE_Y1;
TINI1.'BZ'                = BASE_Z1;
TINI1.'FCD'               = FCD1;
TINI1.'FSD'               = FSD1;
TINI1.'WSY'               = OM_SY1;
TINI1.'WSZ'               = OM_SZ1;
TINI1.'WWY'               = OM_WY1;
TINI1.'WWZ'               = OM_WZ1;
TINI1.'LY'                = LAMB_Y1;
TINI1.'LZ'                = LAMB_Z1;
TINI1.'XIY'               = CHI_Y1;
TINI1.'XIZ'               = CHI_Z1;
TINI1.'PHI_LON'           = LI_PHI1;
TINI1.'PHI_Y'             = PY_PHI1;
TINI1.'PHI_Z'             = PZ_PHI1;
TINI1.'EPSI_BET'          = EP_BET1;
TINI1.'EPSI_ACI'          = EP_ACI1;
TINI1.'YOUNG_ACIER'       = E_ACIE1;
*
'SI' ('EGA' TYP_CAL1 'ELLIPSOIDE');
*
*-------------------------------------------------------------
*   Methode Ellipse
*
* Test about the static state stress
  T_ST1 = G_ULTIFR TINI1;
*
  'SI' (T_ST1 '&lt;EG' 0.0);
    LAM_FIN1 = 0.0;
    TIT1 = 'CHAINE' '  Element ' I_IE1 ' Lambda=' 
                    LAM_FIN1 '.....Sig0 externe';
    'MESS' TIT1;
  'SINON';
    MAI1 = 0. 0. 0.;
*   Matrix X
    MATX1 = VIDE 'RIGIDITE'/'RIGIDITE';
    MATX1 = MATX1 'ET' ('MANU' 'RIGIDITE' MAI1 (MOTS 'UX' 'UY' 'UZ')
          ('PROG' MYMY1 MYMZ1 NNMY1
                  MYMZ1 MZMZ1 NNMZ1
                  NNMY1 NNMZ1 NNNN1));
*
*   Static state stress
    EF_ST0 = 'MANU' 'CHPO' MAI1 3 'FX' MOMY0 'FY' MOMZ0 'FZ' EFFX0;
    EF_ST0 = 'CHAN' 'ATTRIBUT' EF_ST0 'NATURE' 'DISCRET';
*
*
*   Initial Vecteur
    ALP_0 =  'MANU' 'CHPO' MAI1 3 'UX' 2.0 'UY' 2.0 'UZ' 2.0;     
    XALP_0 = (MATX1*ALP_0)/(('XTMX' MATX1 ALP_0)**(0.5));
    XALP_0 = ('CHAN' ATTRIBUT XALP_0 NATURE DISCRET) + EF_ST0;
*
*       direction used for the ellispoide methode
    TB_DIR = 'TABLE';
    TB_DIR. 1 = 'MANU' 'CHPO' MAI1 3 'UX'  1.0 'UY'  0.0 'UZ'  0.0;
    TB_DIR. 2 = 'MANU' 'CHPO' MAI1 3 'UX'  0.0 'UY'  1.0 'UZ'  0.0;
    TB_DIR. 3 = 'MANU' 'CHPO' MAI1 3 'UX'  0.0 'UY'  0.0 'UZ'  1.0;
    TB_DIR. 4 = 'MANU' 'CHPO' MAI1 3 'UX' -1.0 'UY'  0.0 'UZ'  0.0;
    TB_DIR. 5 = 'MANU' 'CHPO' MAI1 3 'UX'  0.0 'UY' -1.0 'UZ'  0.0;
    TB_DIR. 6 = 'MANU' 'CHPO' MAI1 3 'UX'  0.0 'UY'  0.0 'UZ' -1.0;
*
    SB0 =    ('XTX' (XALP_0 - EF_ST0))**(0.5);
    RB0 =     SB0;
    LAMB0 =   RB0/SB0;
*
*   incremental step
    H_VAL1 = (SB0)/(100000.0);
    H_VAL2 = (SB0)/(100000.0);
*  
*   variable for the convergence test
    ID1_CON = 0;
*
* loop for each direction
    I_K1 = 1;
    'REPE' IND3 (6);
*
      TAB_ITE = 'TABLE';
      TAB_ITE. 0 = 'TABLE';
      TAB_ITE. 0 . 'ALPA' = ALP_0;
      TAB_ITE. 0 . 'LAMBDA' = LAMB0;
      TAB_ITE. 0 . 'RB' = RB0;
*
      ALP_K = (TB_DIR. I_K1);
      I_M1 = 1;
*     first loop in order to have n_elllips=n_limite
      'REPE' IND2 (20);
*       sauving Lambda for j iteraction
        TAB_ITE. I_M1 = 'TABLE';
        TAB_ITE. I_M1 . 'ALPA' = ALP_K;
*       determination of the stress state for the chosen direction
        XALP_J = (MATX1*ALP_K)/(('XTMX' MATX1 ALP_K)**(0.5));
        XALP_J = ('CHAN' ATTRIBUT XALP_J NATURE DISCRET) + EF_ST0;
*       nomalisez vector of hte rayon of the ellipsoide for the 
*       chosen direction
        SBJ = ('XTX' (XALP_J - EF_ST0))**(0.5);
        BJ = (XALP_J - EF_ST0)/(SBJ);
*       rayon of the ellipsoide for the chosen direction
        RBVAL1 = (TAB_ITE.(I_M1 - 1). 'RB');
*       second loop in order to have f=0
        I_I1 = 1;
        'REPE' IND1 (200);
*         strees state x0, x0+h, x0-h
          VV1 = (RBVAL1 + H_VAL1)*(BJ);
          VV2 = (RBVAL1 - H_VAL1)*(BJ);
          VV3 = (RBVAL1)*(BJ);
          VV1 = ('CHAN' 'ATTRIBUT' VV1 'NATURE' 'DISCRET') + EF_ST0;
          VV2 = ('CHAN' 'ATTRIBUT' VV2 'NATURE' 'DISCRET') + EF_ST0;
          VV3 = ('CHAN' 'ATTRIBUT' VV3 'NATURE' 'DISCRET') + EF_ST0;
          TINIA1 = 'COPI' TINI1;
          TINIB1 = 'COPI' TINI1;
          TINIC1 = 'COPI' TINI1;
          TINIA1.'NN'                = EXTR VV1 'FZ' MAI1;
          TINIA1.'VY'                = 0.0;
          TINIA1.'VZ'                = 0.0;
          TINIA1.'MT'                = 0.0;
          TINIA1.'MY'                = EXTR VV1 'FX' MAI1;
          TINIA1.'MZ'                = EXTR VV1 'FY' MAI1;
          TINIB1.'NN'                = EXTR VV2 'FZ' MAI1;
          TINIB1.'VY'                = 0.0;
          TINIB1.'VZ'                = 0.0;
          TINIB1.'MT'                = 0.0;
          TINIB1.'MY'                = EXTR VV2 'FX' MAI1;
          TINIB1.'MZ'                = EXTR VV2 'FY' MAI1;
          TINIC1.'NN'                = EXTR VV3 'FZ' MAI1;
          TINIC1.'VY'                = 0.0;
          TINIC1.'VZ'                = 0.0;
          TINIC1.'MT'                = 0.0;
          TINIC1.'MY'                = EXTR VV3 'FX' MAI1;
          TINIC1.'MZ'                = EXTR VV3 'FY' MAI1;
          'OUBL' VV1;
          'OUBL' VV2;
          'OUBL' VV3;
          F_VAL1 = G_ULTIFR TINIC1;
          D_VAL1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/
                                               ((2.0)*(H_VAL1));
*          MESS '-----';
*          LIST F_VAL1;
*          LIST D_VAL1;
*          SI (I_I1 EGA 6);
*            'QUIT' IND1;
*          FINSI;
*         controle
          'SI' ((('ABS'(F_VAL1)) '<' ((TOL_2)/(100.0))) 
                                  'OU' (I_I1 'EGA' 200));
            RBJ = RBVAL1;
            'QUIT' IND1;
          'SINON';
            'SI' (('ABS'(D_VAL1)) < 1.E-12);
              RBVAL2 = RBVAL1;
              RBVAL1 = (RBVAL1)*(1.1);
            'SINON';
*           otherwise
              RBVAL2 = RBVAL1;
              RBVAL1 = RBVAL1 - ((F_VAL1)/(D_VAL1));
            'FINSI';
          'FINSI';
          I_I1 = I_I1 + 1;
        'FIN' IND1;
*       computation of lambda
        LAMBJ = RBJ/SBJ;
        TAB_ITE. I_M1 . 'LAMBDA' = LAMBJ;
        TAB_ITE. I_M1 . 'RB' = RBJ;
        XCOUB1 = (((TAB_ITE. I_M1 .'LAMBDA')*(TAB_ITE. I_M1 .'LAMBDA'))*
                  (MATX1));
*       amplified stress state using the lambda coef
        XALP_L = (XCOUB1*ALP_K)/(('XTMX' XCOUB1 ALP_K)**(0.5));
        XALP_L = ('CHAN' ATTRIBUT XALP_L  NATURE DISCRET) + EF_ST0;
*       Dh*ex  ; Dh*ey ; Dh*ez
        DE_HX = 'MANU' 'CHPO' MAI1 3 'FX' H_VAL2 'FY' 0.0 'FZ' 0.0;
        DE_HY = 'MANU' 'CHPO' MAI1 3 'FX' 0.0 'FY' H_VAL2 'FZ' 0.0;
        DE_HZ = 'MANU' 'CHPO' MAI1 3 'FX' 0.0 'FY' 0.0 'FZ' H_VAL2;
*       x0*lam + Dh*ex  ; x0*lam + Dh*ey ; x0*lam + Dh*ez
        VV_HXP = XALP_L + DE_HX;
        VV_HYP = XALP_L + DE_HY;
        VV_HZP = XALP_L + DE_HZ;
*       x0*lam + Dh*ex  ; x0*lam + Dh*ey ; x0*lam + Dh*ez
        VV_HXN = XALP_L - DE_HX;
        VV_HYN = XALP_L - DE_HY;
        VV_HZN = XALP_L - DE_HZ;
*       x0*lam
        VV_0 = XALP_L;
*
        VV_HXP1 = 'EXTR' VV_HXP 'FX' MAI1;
        VV_HXP2 = 'EXTR' VV_HXP 'FY' MAI1;
        VV_HXP3 = 'EXTR' VV_HXP 'FZ' MAI1;
        VV_HYP1 = 'EXTR' VV_HYP 'FX' MAI1;
        VV_HYP2 = 'EXTR' VV_HYP 'FY' MAI1;
        VV_HYP3 = 'EXTR' VV_HYP 'FZ' MAI1;
        VV_HZP1 = 'EXTR' VV_HZP 'FX' MAI1;
        VV_HZP2 = 'EXTR' VV_HZP 'FY' MAI1;
        VV_HZP3 = 'EXTR' VV_HZP 'FZ' MAI1;
        VV_HXN1 = 'EXTR' VV_HXN 'FX' MAI1;
        VV_HXN2 = 'EXTR' VV_HXN 'FY' MAI1;
        VV_HXN3 = 'EXTR' VV_HXN 'FZ' MAI1;
        VV_HYN1 = 'EXTR' VV_HYN 'FX' MAI1;
        VV_HYN2 = 'EXTR' VV_HYN 'FY' MAI1;
        VV_HYN3 = 'EXTR' VV_HYN 'FZ' MAI1;
        VV_HZN1 = 'EXTR' VV_HZN 'FX' MAI1;
        VV_HZN2 = 'EXTR' VV_HZN 'FY' MAI1;
        VV_HZN3 = 'EXTR' VV_HZN 'FZ' MAI1;
*
        VV_01 = 'EXTR' VV_0 'FX' MAI1;
        VV_02 = 'EXTR' VV_0 'FY' MAI1;
        VV_03 = 'EXTR' VV_0 'FZ' MAI1;
*
 
        TINIA1 = 'COPI' TINI1;
        TINIB1 = 'COPI' TINI1;
        TINIA1.'NN'                = VV_HXP3;
        TINIA1.'VY'                = 0.0;
        TINIA1.'VZ'                = 0.0;
        TINIA1.'MT'                = 0.0;
        TINIA1.'MY'                = VV_HXP1;
        TINIA1.'MZ'                = VV_HXP2;
        TINIB1.'NN'                = VV_HXN3;
        TINIB1.'VY'                = 0.0;
        TINIB1.'VZ'                = 0.0;
        TINIB1.'MT'                = 0.0;
        TINIB1.'MY'                = VV_HXN1;
        TINIB1.'MZ'                = VV_HXN2;
        VAL_X1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/
                                                ((2.0)*(H_VAL2));
*
        TINIA1 = 'COPI' TINI1;
        TINIB1 = 'COPI' TINI1;
        TINIA1.'NN'                = VV_HYP3;
        TINIA1.'VY'                = 0.0;
        TINIA1.'VZ'                = 0.0;
        TINIA1.'MT'                = 0.0;
        TINIA1.'MY'                = VV_HYP1;
        TINIA1.'MZ'                = VV_HYP2;
        TINIB1.'NN'                = VV_HYN3;
        TINIB1.'VY'                = 0.0;
        TINIB1.'VZ'                = 0.0;
        TINIB1.'MT'                = 0.0;
        TINIB1.'MY'                = VV_HYN1;
        TINIB1.'MZ'                = VV_HYN2;
        VAL_Y1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/
                                                ((2.0)*(H_VAL2));
*
        TINIA1 = 'COPI' TINI1;
        TINIB1 = 'COPI' TINI1;
        TINIA1.'NN'                = VV_HZP3;
        TINIA1.'VY'                = 0.0;
        TINIA1.'VZ'                = 0.0;
        TINIA1.'MT'                = 0.0;
        TINIA1.'MY'                = VV_HZP1;
        TINIA1.'MZ'                = VV_HZP2;
        TINIB1.'NN'                = VV_HZN3;
        TINIB1.'VY'                = 0.0;
        TINIB1.'VZ'                = 0.0;
        TINIB1.'MT'                = 0.0;
        TINIB1.'MY'                = VV_HZN1;
        TINIB1.'MZ'                = VV_HZN2;
        VAL_Z1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/
                                                ((2.0)*(H_VAL2));
*
*       costruction n_limit
        DELTG1 = 'MANU' 'CHPO' MAI1 'UX' ((-1.0)*(VAL_X1))
                                    'UY' ((-1.0)*(VAL_Y1))
                                  'UZ' ((-1.0)*(VAL_Z1));
        DELTG1 = DELTG1/((XTX DELTG1)**(0.5));
*       check on the n_limit
        'SI' ((('XTX' DELTG1)**(0.5)) < 1.E-7);
*         if |n_limit|=0 rotate the chosen direction by 10 degrees
          MATR1 = 'VIDE' 'RIGIDITE'/'RIGIDITE';
          MATR1 = MATR1 'ET' 
             ('MANU' 'RIGIDITE' MAI1 ('MOTS' 'UX' 'UY' 'UZ') 'QUEL'
                ('PROG'  (COS(10.0)) ((-1.0)*(SIN(10.0))) 0.0
                         (SIN(10.0))      (COS(10.0))     0.0
                             0.0             0.0          1.0));
          AAA1 = MATR1*ALP_K;
          DELTG1 = MANU CHPO MAI1 'UX' (EXTR AAA1 MAI1 'FX')
                                  'UY' (EXTR AAA1 MAI1 'FY')
                                  'UZ' (EXTR AAA1 MAI1 'FZ');
        'FINSI';
*       convergence test
        VET_TEST = (TAB_ITE. I_M1 . 'ALPA') 
                 - (TAB_ITE. (I_M1 - 1) . 'ALPA');  
        'SI' ((('XTX' VET_TEST)**(0.5)) '<' TOL_2);
*         deduction of the final value of Lambda
          'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0);
            'SI' ('NEG' ('TYPE' LAM_FIN1) 'FLOTTANT');
              LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA');
            'SINON';
              'SI' (LAM_FIN1 '>' (TAB_ITE. I_M1 . 'LAMBDA'));
                LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA');
              'FINSI';
            'FINSI';
          'FINSI';
*          MESS '-----------------------------';
*          LIST F_VAL1;
*          LIST I_I1;
*          LIST LAMBJ;
*          LIST I_M1;
*         MESS '*';
*          LIST NLIMT1;
*          LIST NLIMC1;
*          LIST TINIC1.'MY';
*          LIST TINIC1.'MZ';
*          LIST TINIC1.'NN';
*          MESS '*';
*          LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UX' MAI1);
*          LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UY' MAI1);
*          LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UZ' MAI1);
*          LIST ('EXTR' DELTG1 'UX' MAI1);
*          LIST ('EXTR' DELTG1 'UY' MAI1);
*          LIST ('EXTR' DELTG1 'UZ' MAI1);
*          MESS '-----------------------------';
          'QUIT' IND2;
        'SINON';
          'SI' (I_M1 'EGA' 50);
            ID1_CON = 1;
            'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0);
              'SI' ('NEG' ('TYPE' LAM_FIN1) 'FLOTTANT');
                LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA');
              'SINON';
                'SI' (LAM_FIN1 '>' (TAB_ITE. I_M1 . 'LAMBDA'));
                  LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA');
                'FINSI';
              'FINSI';
            'FINSI';
          'SINON';
*           otherwise updating the n_limit
            ALP_K = DELTG1;
          'FINSI';
        'FINSI';
*       Display
*        MESS '-----------------------------';
*        LIST F_VAL1;
*        LIST I_I1;
*        LIST LAMBJ;
*        LIST I_M1;
*       MESS '*';
*        LIST NLIMT1;
*        LIST NLIMC1;
*        LIST TINIC1.'MY';
*        LIST TINIC1.'MZ';
*        LIST TINIC1.'NN';
*        MESS '*';
*        LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UX' MAI1);
*        LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UY' MAI1);
*        LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UZ' MAI1);
*        LIST ('EXTR' DELTG1 'UX' MAI1);
*        LIST ('EXTR' DELTG1 'UY' MAI1);
*        LIST ('EXTR' DELTG1 'UZ' MAI1);
*        MESS '/////';
*        LIST (('XTX' VET_TEST)**(0.5));
*        MESS '-----------------------------';
*        SI (I_K1 > 1);
*          ERRE 5;
*        FINSI;
        SI (I_M1 EGA 30);
          ERRE 5;
        FINSI;
        I_M1 = I_M1 + 1;
      'FIN' IND2;
      I_K1 = I_K1 + 1;
    'FIN' IND3;
*
*    screen message
    'SI' (ID1_CON 'EGA' 0);
      TIT1 = 'CHAINE' '  Element ' I_ELE1 ' Lambda=' 
                      LAM_FIN1;
      'MESS' TIT1;
    'FINSI';
    'SI' (ID1_CON 'EGA' 1);
      TIT1 = 'CHAINE' '  Element ' I_ELE1 ' Lambda=' 
                       LAM_FIN1 '.....Attention non convergence';
      'MESS' TIT1;
    'FINSI';
  'FINSI';
'SINON';
*
*-------------------------------------------------------------
*   Methode Rectangle
*
*
* Test about the static state stress
  T_ST1 = G_ULTIFR TINI1;
*
  ID1_CON = 0;
  'SI' (T_ST1 '&lt;EG' 0.0);
    LAM_FIN1 = 0.0;
    TIT1 = 'CHAINE' '  Element ' I_IE1 ' Lambda=' 
                    LAM_FIN1 '.....Sig0 externe';
    'MESS' TIT1;
  'SINON';
    LAMB0 =  1.0; 
    H_VAL3 = LAMB0/100;
    TOL_2 = FCBET1/1000;
    LL_LAMX = PROG -1.0 1.0 1.0 -1.0 -1.0 1.0 1.0 -1.0;
    LL_LAMY = PROG -1.0 -1.0 1.0 1.0 -1.0 -1.0 1.0 1.0;
    LL_LAMZ = PROG -1.0 -1.0 -1.0 -1.0 1.0 1.0 1.0 1.0;
*
*   reference list for the 8 point for rectangle solid
    LL_I = PROG 1. 2. 3. 4. 5. 6. 7. 8.;
*
*   creation of the resuming table
    TAB_ITE = 'TABLE';
    TAB_ITE. 0 = 'TABLE';
    TAB_ITE. 0 .'LAMBDA' = LAMB0 ;
*
    I_MM1 = 0;
    ID1_CON = 0;
*
*   loop for in order to find f=0
    'REPE' IND4 (200);
*
*     stress state for the 8 points
      X1_0 = (((MYMY1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMY0;
      Y1_0 = (((MZMZ1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMZ0;
      Z1_0 = (((NNNN1)**(0.5))*(((-1.0)*(LAMB0)))) + EFFX0;
      X2_0 = (((MYMY1)**(0.5))*(((1.0)*(LAMB0)))) + MOMY0;
      Y2_0 = (((MZMZ1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMZ0;
      Z2_0 = (((NNNN1)**(0.5))*(((-1.0)*(LAMB0)))) + EFFX0;
      X3_0 = (((MYMY1)**(0.5))*(((1.0)*(LAMB0)))) + MOMY0;
      Y3_0 = (((MZMZ1)**(0.5))*(((1.0)*(LAMB0)))) + MOMZ0;
      Z3_0 = (((NNNN1)**(0.5))*(((-1.0)*(LAMB0)))) + EFFX0;
      X4_0 = (((MYMY1)**(0.5))*(((1.0)*(LAMB0)))) + MOMY0;
      Y4_0 = (((MZMZ1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMZ0;
      Z4_0 = (((NNNN1)**(0.5))*(((-1.0)*(LAMB0)))) + EFFX0;
      X5_0 = (((MYMY1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMY0;
      Y5_0 = (((MZMZ1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMZ0;
      Z5_0 = (((NNNN1)**(0.5))*(((1.0)*(LAMB0)))) + EFFX0;
      X6_0 = (((MYMY1)**(0.5))*(((1.0)*(LAMB0)))) + MOMY0;
      Y6_0 = (((MZMZ1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMZ0;
      Z6_0 = (((NNNN1)**(0.5))*(((1.0)*(LAMB0)))) + EFFX0;
      X7_0 = (((MYMY1)**(0.5))*(((1.0)*(LAMB0)))) + MOMY0;
      Y7_0 = (((MZMZ1)**(0.5))*(((1.0)*(LAMB0)))) + MOMZ0;
      Z7_0 = (((NNNN1)**(0.5))*(((1.0)*(LAMB0)))) + EFFX0;
      X8_0 = (((MYMY1)**(0.5))*(((1.0)*(LAMB0)))) + MOMY0;
      Y8_0 = (((MZMZ1)**(0.5))*(((-1.0)*(LAMB0)))) + MOMZ0;
      Z8_0 = (((NNNN1)**(0.5))*(((1.0)*(LAMB0)))) + EFFX0;
*
*     put the N11, N22 et N12 stress state in the lists
      LL_X0 = 'PROG' (X1_0) (X2_0) (X3_0) (X4_0)
                     (X5_0) (X6_0) (X7_0) (X8_0);
      LL_Y0 = 'PROG' (Y1_0) (Y2_0) (Y3_0) (Y4_0)
                     (Y5_0) (Y6_0) (Y7_0) (Y8_0);
      LL_Z0 = 'PROG' (Z1_0) (Z2_0) (Z3_0) (Z4_0)
                     (Z5_0) (Z6_0) (Z7_0) (Z8_0);
      LL_F0 = 'PROG';
*
*     deduction of f for each point
      III1 = 1;
      'REPE' IND2 (8);
        TINIC1 = 'COPI' TINI1;
        TINIC1.'NN'                = ('EXTR' LL_Z0 III1);
        TINIC1.'VY'                = 0.0;
        TINIC1.'VZ'                = 0.0;
        TINIC1.'MT'                = 0.0;
        TINIC1.'MY'                = ('EXTR' LL_X0 III1);
        TINIC1.'MZ'                = ('EXTR' LL_Y0 III1);
        LL_F0 = LL_F0 ET 
          ('PROG' (G_ULTIFR TINIC1));
        III1 = III1 + 1;
      'FIN' IND2;
*     saving 
      TAB_ITE. I_MM1 .'L_X' = LL_X0;
      TAB_ITE. I_MM1 .'L_Y' = LL_Y0;
      TAB_ITE. I_MM1 .'L_Z' = LL_Z0;
      TAB_ITE. I_MM1 .'L_F' = LL_F0;
*     finding the minimum value of f and the corresponding point
      EV1 = 'EVOL' 'MANU' LL_I (LL_F0);
      II1 II2 VAL_F = 'MINI' EV1;
*     computation of D (gradiant) in case of Point_1=Argmin(fmin)
      X1_DP = 
           (((MYMY1)**(0.5))*((('EXTR' LL_LAMX II1)*(LAMB0 + H_VAL3)))) 
                                                          + MOMY0;
      Y1_DP = 
           (((MZMZ1)**(0.5))*((('EXTR' LL_LAMY II1)*(LAMB0 + H_VAL3)))) 
                                                          + MOMZ0;
      Z1_DP = 
           (((NNNN1)**(0.5))*((('EXTR' LL_LAMZ II1)*(LAMB0 + H_VAL3)))) 
                                                          + EFFX0;
      X1_0 = 
           (((MYMY1)**(0.5))*((('EXTR' LL_LAMX II1)*(LAMB0)))) 
                                                          + MOMY0;
      Y1_0 = 
           (((MZMZ1)**(0.5))*((('EXTR' LL_LAMY II1)*(LAMB0)))) 
                                                          + MOMZ0;
      Z1_0 = 
           (((NNNN1)**(0.5))*((('EXTR' LL_LAMZ II1)*(LAMB0)))) 
                                                          + EFFX0;
      X1_DN = 
           (((MYMY1)**(0.5))*((('EXTR' LL_LAMX II1)*(LAMB0 - H_VAL3)))) 
                                                          + MOMY0;
      Y1_DN = 
           (((MZMZ1)**(0.5))*((('EXTR' LL_LAMY II1)*(LAMB0 - H_VAL3)))) 
                                                          + MOMZ0;
      Z1_DN = 
           (((NNNN1)**(0.5))*((('EXTR' LL_LAMZ II1)*(LAMB0 - H_VAL3)))) 
                                                          + EFFX0;
      TINIA1 = 'COPI' TINI1;
      TINIB1 = 'COPI' TINI1;
      TINIA1.'NN'                = Z1_DP;
      TINIA1.'VY'                = 0.0;
      TINIA1.'VZ'                = 0.0;
      TINIA1.'MT'                = 0.0;
      TINIA1.'MY'                = X1_DP;
      TINIA1.'MZ'                = Y1_DP;
      TINIB1.'NN'                = Z1_DN;
      TINIB1.'VY'                = 0.0;
      TINIB1.'VZ'                = 0.0;
      TINIB1.'MT'                = 0.0;
      TINIB1.'MY'                = X1_DN;
      TINIB1.'MZ'                = Y1_DN;
      D_VAL1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/
                                              ((2.0)*(H_VAL3));
 
      'OUBL' X1_DP;
      'OUBL' Y1_DP;
      'OUBL' Z1_DP;
      'OUBL' X1_DN;
      'OUBL' Y1_DN;
      'OUBL' Z1_DN;
*      'OUBL' X1_0;
*      'OUBL' Y1_0;
*      'OUBL' Z1_0; 
*
*     Display
*      MESS '-----------------------';
*      LIST II1;
*      LIST LAMB0;
*      LIST VAL_F;
*      LIST D_VAL1;
*      LIST X1_0;
*      LIST Y1_0;
*      LIST Z1_0;
*      SI (I_MM1 EGA 5);
*        ERRE 5;
*      FINSI;
*
*     test convergence
      'SI' (('ABS'(VAL_F)) '<' TOL_1);
*       convergence!!!
        LAM_FIN1 = LAMB0;
        'QUIT' IND4;
      'SINON';
        'SI' (I_MM1 'EGA' 199);
          LAM_FIN1 = LAMB0;
          ID1_CON = 1;
          'QUIT' IND4;
        'FINSI';
        'SI' (('ABS'(D_VAL1)) < 1.E-12);
          LAMB1 = LAMB0;
          LAMB0 = LAMB0*1.10;
        'SINON';
*         otherwise
          LAMB1 = LAMB0;
          LAMB0 = LAMB0 - ((VAL_F)/(D_VAL1));
          'SI' (((ABS(LAMB1 - LAMB0))/(LAMB0)) '<' (1.E-8));
             LAMB0 = ((LAMB1)*(0.90));
          'FINSI';
        'FINSI';
      'FINSI';
      'OUBL' VAL_F;
      'OUBL' D_VAL1;
*
      TAB_ITE. (I_MM1 + 1) = 'TABLE';
      TAB_ITE. (I_MM1 + 1) .'LAMBDA' = ABS(LAMB0) ;
*
      I_MM1 = I_MM1 + 1;
    'FIN' IND4;
 
  'FINSI';
*
*    screen message
    'SI' (ID1_CON 'EGA' 0);
      TIT1 = 'CHAINE' '  Element ' I_ELE1 ' Lambda=' 
                      LAM_FIN1;
      'MESS' TIT1;
    'FINSI';
    'SI' (ID1_CON 'EGA' 1);
      TIT1 = 'CHAINE' '  Element ' I_ELE1 ' Lambda=' 
                       LAM_FIN1 '.....Attention non convergence';
    'FINSI';
*    LIST TAB_ITE;
*    LIST I_MM1;
*    LIST TAB_ITE.I_MM1;
*ERRE 5;
*
'FINSI';
FINPROC LAM_FIN1;