* fichier : mrcframe_test.dgibi
*
**********************************************************
*  Test sur la procedure MRCFRAME, fonction pour determiner 
*  la marge de securité d'un element type POUT et TIMO soumis à un 
*  chargement sismique. Les effort sont calculées selon la methode
*  spectrale avec et sans la prise en compte des termes de covariances
*
*
*  Section
*
*  b=0.3
*  h=0.4
*  c=0.025
*  armature long 4f12
*  cadre phi 8
*  At = f8 mm
*
*  4 Cast test:
*    - V_TEST1 = 1 Poteau court         ly=1.0 lz=1.0 s = 0.95 m
*    - V_TEST1 = 2 Poutre courte        ly=1.0 lz=1.0 s = 0.95 m
*    - V_TEST1 = 3 Poteau long          ly=2.0 lz=1.5 s = 0.1 m
*    - V_TEST1 = 4 Poutre longe         ly=2.0 lz=1.5 s = 0.1 m
*
*  Pour V_TEST1 = 1, 3 type de chargement:
*     - Cha1: N0 = -1600.0 kN
*             Nmax = 10 kN
*             Vymax = 200 kNm
*             Vzmax = 10 kNm
*     - Cha2: N0 = -1600.0 kN
*             Nmax = 10 kN
*             Vymax = 10 kNm
*             Vzmax = 200 kNm
*     - Cha3: N0 = -1600.0 kN
*             Nmax = 700 kN
*             Vymax = 10 kNm
*             Vzmax = 10 kNm    
*  Pour V_TEST1 = 2, 3 type de chargement:   
*     - Cha1: N0 = -1600.0 kN
*             Nmax = 10 kN
*             Vymax = 10 kNm
*             Vzmax = 200 kNm
*     - Cha3: N0 = -1600.0 kN
*             Nmax = 700 kN
*             Vymax = 10 kNm
*             Vzmax = 10 kNm       
*  Pour V_TEST1 = 3, 3 type de chargement:
*     - Cha1: N0 = -1600.0 kN
*             Nmax = 1 kN
*             Mymax = 10 kNm
*             Mzmax = 1 kNm
*     - Cha2: N0 = -1600.0 kN
*             Nmax = 1 kN
*             Mymax = 1 kNm
*             Mzmax = 10 kNm
*     - Cha3: N0 = -1600.0 kN
*             Nmax = 700 kN
*             Mymax = 1 kNm
*             Mzmax = 1 kNm  
*  Pour V_TEST1 = 4, 2 type de chargement:
*     - Cha1: N0 = -1600.0 kN
*             Nmax = 1 kN
*             Mymax = 10 kNm
*             Mzmax = 1 kNm
*     - Cha3: N0 = -1600.0 kN
*             Nmax = 700 kN
*             Mymax = 1 kNm
*             Mzmax = 1 kNm  
*
*
*
*  Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
*              et Nicolas ILE  /DEN/DANS/DM2S/SEMT/EMSI
*
**********************************************************
*
OPTI DIME 3 MODE TRID ELEM SEG2;
OPTI TRAC PSC;
*
* dimensions de la section 
B_Y1 = 0.3;
B_Z1 = 0.4;
* barres longitudianles
PHI_LON1 = 12.0;
* enrobage
ENR1 = 0.025;
* Armature transversales
A_CADRE1 = (2.0)*((PI)*((0.004)**(2.0)));
A_CADRE2 = (2.0)*((PI)*((0.004)**(2.0)));
* barres longitudianles
PHI_LON1 = PROG 12.0 12.0 12.0 12.0;
* Position y barres
PHI_Y = PROG ((-1.0)*((B_Y1 - (2.0*ENR1))/(2.0))) 
             ((+1.0)*((B_Y1 - (2.0*ENR1))/(2.0)))
             ((+1.0)*((B_Y1 - (2.0*ENR1))/(2.0)))
             ((-1.0)*((B_Y1 - (2.0*ENR1))/(2.0)));
* Position z barres
PHI_Z = PROG ((-1.0)*((B_Z1 - (2.0*ENR1))/(2.0))) 
             ((-1.0)*((B_Z1 - (2.0*ENR1))/(2.0)))
             ((+1.0)*((B_Z1 - (2.0*ENR1))/(2.0)))
             ((+1.0)*((B_Z1 - (2.0*ENR1))/(2.0)));
*
*
V_TEST1 = 1;
*V_TEST1 = 2;
*V_TEST1 = 3;
*V_TEST1 = 4;
*
*
SI (V_TEST1 EGA 1);
  V_CHA1 = 1;
*  V_CHA1 = 2;
*  V_CHA1 = 3;
FINSI;
SI (V_TEST1 EGA 2);
  V_CHA1 = 1;
*  V_CHA1 = 3;
FINSI;
SI (V_TEST1 EGA 3);
  V_CHA1 = 1;
*  V_CHA1 = 2;
*  V_CHA1 = 3;
FINSI;
SI (V_TEST1 EGA 4);
  V_CHA1 = 1;
*  V_CHA1 = 3;
FINSI;
*
*
*
P1 = 0. 0. 0.;
P2 = 1. 0. 0.;
L1 = D 1 P1 P2;
*
*
*------------ Cas Pouteau et Poutre courte ----------------------
SI ((V_TEST1 EGA 1) OU  (V_TEST1 EGA 2));
*
* longuer ly lz
  LTOT1Y = 1.0;
  LTOT1Z = 1.5;
* espacement cadres
  S_CAD1 = 0.95;
* calcul aire armature longitudinales
  ALL1 = 0.0;
  I = 1;
  NN1 = DIME PHI_LON1;
*
  REPE IND1 (NN1);
    PHI_I1 = EXTR PHI_LON1 I;
    AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
    AA1 = AA1/(1.E6);
    ALL1 = ALL1 + AA1;
    I = I + 1;
  FIN IND1;
* calcul omega_s et omega_w
  OM_SY1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  OM_SZ1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  OM_WZ1 = ((A_CADRE1)*(500.E6))/(((B_Y1)*(S_CAD1))*(25.E6));
  OM_WY1 = ((A_CADRE2)*(500.E6))/(((B_Z1)*(S_CAD1))*(25.E6));
  LAMB_Y1 = ((LTOT1Z)/(B_Y1));
  LAMB_Z1 = ((LTOT1Y)/(B_Z1));
  CHI_Z1 = ((B_Z1 - ((2.0)*(ENR1)))/(B_Z1));
  CHI_Y1 = ((B_Y1 - ((2.0)*(ENR1)))/(B_Y1));
* test 
  VAL_TT1 = LAMB_Y1 <EG (OM_SY1/OM_WY1);
  VAL_TT2 = LAMB_Z1 <EG (OM_SZ1/OM_WZ1);
* vymax et vzmax
  TAN_Y1 = (((LAMB_Y1**2.0) + 1.0)**(0.5)) - LAMB_Y1;
  TAN_Z1 = (((LAMB_Z1**2.0) + 1.0)**(0.5)) - LAMB_Z1;
  VY_MAX1 = OM_WY1*CHI_Y1 + ((0.5 - OM_WY1)*(TAN_Y1));
  VZ_MAX1 = OM_WZ1*CHI_Z1 + ((0.5 - OM_WZ1)*(TAN_Z1));
  VY_MAX1 = ((VY_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  VZ_MAX1 = ((VZ_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
* n traction et n compression plus nlim
  N_TRAC1 = ((OM_SY1)*(((B_Y1)*(B_Z1))*(25.e6)));
  N_COMP1 = ((1.0 + OM_SY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMY1 = 0.5 + OM_SY1 - ((OM_WY1)*(LAMB_Y1 + 1.0 - CHI_Y1));
  N_LIMY2 = 0.5 - OM_SY1 + ((OM_WY1)*(LAMB_Y1 - 1.0 + CHI_Y1));
  N_LIMZ1 = 0.5 + OM_SZ1 - ((OM_WZ1)*(LAMB_Z1 + 1.0 - CHI_Z1));
  N_LIMZ2 = 0.5 - OM_SZ1 + ((OM_WZ1)*(LAMB_Z1 - 1.0 + CHI_Z1));
  N_LIMY1 = ((N_LIMY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMY2 = ((N_LIMY2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMZ1 = ((N_LIMZ1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMZ2 = ((N_LIMZ2)*(((B_Y1)*(B_Z1))*(-25.e6)));
*
  VY_MAX1 = VY_MAX1/1000.;
  VZ_MAX1 = VZ_MAX1/1000.;
  N_TRAC1 = N_TRAC1/1000.;
  N_COMP1 = N_COMP1/1000.;
  N_LIMY1 = N_LIMY1/1000.;
  N_LIMY2 = N_LIMY2/1000.;
  N_LIMZ1 = N_LIMZ1/1000.;
  N_LIMZ2 = N_LIMZ2/1000.;
*
  SI (V_TEST1 EGA 1);
    TOL_1 = 1.E-8;
    TOL_2 = 1.E-6;
  FINSI;
  SI (V_TEST1 EGA 1);
    TOL_1 = 2.E-8;
    TOL_2 = 1.E-6;
  FINSI;
*
FINSI;
*
*
*
*
*------------ Cas Pouteau et Poutre long ----------------------
SI ((V_TEST1 EGA 3) OU  (V_TEST1 EGA 4));
*
* longuer ly lz
  LTOT1Y = 2.0;
  LTOT1Z = 1.5;
* espacement cadres
  S_CAD1 = 0.1;
* calcul aire armature longitudinales
  ALL1 = 0.0;
  I = 1;
  NN1 = DIME PHI_LON1;
*
  REPE IND1 (NN1);
    PHI_I1 = EXTR PHI_LON1 I;
    AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
    AA1 = AA1/(1.E6);
    ALL1 = ALL1 + AA1;
    I = I + 1;
  FIN IND1;
* calcul omega_s et omega_w
  OM_SY1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  OM_SZ1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  OM_WZ1 = ((A_CADRE1)*(500.E6))/(((B_Y1)*(S_CAD1))*(25.E6));
  OM_WY1 = ((A_CADRE2)*(500.E6))/(((B_Z1)*(S_CAD1))*(25.E6));
  LAMB_Y1 = ((LTOT1Z)/(B_Y1));
  LAMB_Z1 = ((LTOT1Y)/(B_Z1));
  CHI_Z1 = ((B_Z1 - ((2.0)*(ENR1)))/(B_Z1));
  CHI_Y1 = ((B_Y1 - ((2.0)*(ENR1)))/(B_Y1));
* test 
  VAL_TT1 = LAMB_Y1 <EG (OM_SY1/OM_WY1);
  VAL_TT2 = LAMB_Z1 <EG (OM_SZ1/OM_WZ1);
* vymax et vzmax
  TAN_Y1 = (((LAMB_Y1**2.0) + 1.0)**(0.5)) - LAMB_Y1;
  TAN_Z1 = (((LAMB_Z1**2.0) + 1.0)**(0.5)) - LAMB_Z1;
  VY_MAX1 = OM_WY1*CHI_Y1 + ((0.5 - OM_WY1)*(TAN_Y1));
  VZ_MAX1 = OM_WZ1*CHI_Z1 + ((0.5 - OM_WZ1)*(TAN_Z1));
  VY_MAX1 = ((VY_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  VZ_MAX1 = ((VZ_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
* n traction et n compression plus nlim
  N_TRAC1 = ((OM_SY1)*(((B_Y1)*(B_Z1))*(25.e6)));
  N_COMP1 = ((1.0 + OM_SY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMY1 = 0.5 + OM_SY1 - ((OM_WY1)*(LAMB_Y1 + 1.0 - CHI_Y1));
  N_LIMY2 = 0.5 - OM_SY1 + ((OM_WY1)*(LAMB_Y1 - 1.0 + CHI_Y1));
  N_LIMZ1 = 0.5 + OM_SZ1 - ((OM_WZ1)*(LAMB_Z1 + 1.0 - CHI_Z1));
  N_LIMZ2 = 0.5 - OM_SZ1 + ((OM_WZ1)*(LAMB_Z1 - 1.0 + CHI_Z1));
  N_LIMY1 = ((N_LIMY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMY2 = ((N_LIMY2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMZ1 = ((N_LIMZ1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  N_LIMZ2 = ((N_LIMZ2)*(((B_Y1)*(B_Z1))*(-25.e6)));
*
  VY_MAX1 = VY_MAX1/1000.;
  VZ_MAX1 = VZ_MAX1/1000.;
  N_TRAC1 = N_TRAC1/1000.;
  N_COMP1 = N_COMP1/1000.;
  N_LIMY1 = N_LIMY1/1000.;
  N_LIMY2 = N_LIMY2/1000.;
  N_LIMZ1 = N_LIMZ1/1000.;
  N_LIMZ2 = N_LIMZ2/1000.;
*N traction
  N_TRAC1 =((ALL1)*(500.E6));
*N compression
  N_COMP1 = (((B_Y1)*(B_Z1))*(25.E6))+((ALL1)*(500.E6));
  N_COMP1 = ((-1.0)*(N_COMP1));
*
*DEFORMATION DU BETON
  EPSI_BET = 0.0035;
*DEFORMATION DE L'ACIER
  EPSI_ACI = 0.045;
*
*HAUTEUR NETTE
  D_YP = B_Z1-ENR1;
  D_ZP = B_Y1-ENR1;
  D_YN = (B_Z1-ENR1);
  D_ZN = (B_Y1-ENR1);
*
*position de L'AXE NEUTRE
  AXEN_YP = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_YP;
  AXEN_ZP = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_ZP;
  AXEN_YN = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_YN;
  AXEN_ZN = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_ZN;
*
  FF_YP = ((EPSI_BET)/(AXEN_YP));
  FF_ZP = ((EPSI_BET)/(AXEN_ZP));
  FF_YN = (-1.0)*((EPSI_BET)/(AXEN_YN));
  FF_ZN = (-1.0)*((EPSI_BET)/(AXEN_ZN));
  EPS0_YP = (FF_YP)*((B_Z1/2.0) - AXEN_YP);
  EPS0_ZP = (FF_ZP)*((B_Y1/2.0) - AXEN_ZP);
  EPS0_YN = (FF_YN)*(AXEN_YN - (B_Z1/2.0));
  EPS0_ZN = (FF_ZN)*(AXEN_ZN - (B_Y1/2.0));
*
*N BALANCED
  NBIL_YP = ((-0.8)*(25.E6))*((B_Y1)*(AXEN_YP));
  NBIL_YN = ((-0.8)*(25.E6))*((B_Y1)*(AXEN_YN));
  NBIL_ZP = ((-0.8)*(25.E6))*((B_Z1)*(AXEN_ZP));
  NBIL_ZN = ((-0.8)*(25.E6))*((B_Z1)*(AXEN_ZN));
*M BALANCED
  MBIL_YP = ((1.0)*(NBIL_YP))*((B_Z1/2.0) - ((0.4)*(AXEN_YP)));
  MBIL_ZP = ((1.0)*(NBIL_ZP))*((B_Y1/2.0) - ((0.4)*(AXEN_ZP)));
  MBIL_YN = ((-1.0)*(NBIL_YN))*((B_Z1/2.0) - ((0.4)*(AXEN_YN)));
  MBIL_ZN = ((-1.0)*(NBIL_ZN))*((B_Y1/2.0) - ((0.4)*(AXEN_ZN)));
*
  I = 1;
  REPE IND1 (NN1);
* Calcolo di epsi
    EPSI_YP = EPS0_YP - ((FF_YP)*(EXTR PHI_Z I));
    EPSI_ZP = EPS0_ZP - ((FF_ZP)*(EXTR PHI_Y I));
    EPSI_YN = EPS0_YN - ((FF_YN)*(EXTR PHI_Z I));
    EPSI_ZN = EPS0_ZN - ((FF_ZN)*(EXTR PHI_Y I));
*calcolo di sigma
    SIGI_YP = (210000.0E6)*(EPSI_YP);
    SIGI_YN = (210000.0E6)*(EPSI_YN);
    SIGI_ZP = (210000.0E6)*(EPSI_ZP);
    SIGI_ZN = (210000.0E6)*(EPSI_ZN);
    SI ((ABS(SIGI_YP)) > 500.E6);
      SIGI_YP = (SIGI_YP/(ABS(SIGI_YP)))*(500.E6);
    FINSI;
    SI ((ABS(SIGI_YN)) > 500.E6);
      SIGI_YN = (SIGI_YN/(ABS(SIGI_YN)))*(500.E6);
    FINSI;
    SI ((ABS(SIGI_ZP)) > 500.E6);
      SIGI_ZP = (SIGI_ZP/(ABS(SIGI_ZP)))*(500.E6);
    FINSI;
    SI ((ABS(SIGI_ZN)) > 500.E6);
      SIGI_ZN = (SIGI_ZN/(ABS(SIGI_ZN)))*(500.E6);
    FINSI;
*
    PHI_I1 = EXTR PHI_LON1 I;
    AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
    AA1 = AA1/(1.E6);
*
    NBIL_YP = NBIL_YP + ((AA1)*(SIGI_YP));
    NBIL_YN = NBIL_YN + ((AA1)*(SIGI_YN));
    NBIL_ZP = NBIL_ZP + ((AA1)*(SIGI_ZP));
    NBIL_ZN = NBIL_ZN + ((AA1)*(SIGI_ZN));
*
    MBIL_YP = MBIL_YP + (((AA1)*(SIGI_YP))*(EXTR PHI_Z I));
    MBIL_YN = MBIL_YN + (((AA1)*(SIGI_YN))*(EXTR PHI_Z I));
    MBIL_ZP = MBIL_ZP + (((AA1)*(SIGI_ZP))*(EXTR PHI_Y I));
    MBIL_ZN = MBIL_ZN + (((AA1)*(SIGI_ZN))*(EXTR PHI_Y I));
*
  I = I + 1;
  FIN IND1;
*
  MBIL_YP = MBIL_YP/1000.;
  MBIL_YN = MBIL_YN/1000.;
  MBIL_ZP = MBIL_ZP/1000.;
  MBIL_ZN = MBIL_ZN/1000.;
  NBIL_YP = NBIL_YP/1000.;
  NBIL_YN = NBIL_YN/1000.;
  NBIL_ZP = NBIL_ZP/1000.;
  NBIL_ZN = NBIL_ZN/1000.;
  N_TRAC1 = N_TRAC1/1000.;
  N_COMP1 = N_COMP1/1000.;
*
  SI (V_TEST1 EGA 3);
    TOL_1 = 1.E-8;
    TOL_2 = 1.E-6;
  FINSI;
  SI (V_TEST1 EGA 4);
    TOL_1 = 1.E-6;
    TOL_2 = 1.E-6;
  FINSI;
*
FINSI;
*
*
*
*
* -------------------- Cas Pouteau courte ----------------------
SI (V_TEST1 EGA 1);
* Caracteristiques
  CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
                         'LIBZ' LTOT1Z 
                        'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
                        'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
                        'ASWZ' A_CADRE1    'YACI' 210000.0E6  
                        'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POTEAU';
*
*
* Chargement CAS1
  SI (V_CHA1 EGA 1);
    CH_SEI1 = MANU CHML L1 
        'C11' (10.E3*10.E3) 'C22' (10.E3*10.E3) 'C33' (200.E3*200.E3)
              'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
              'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
              'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
              'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
              'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
              'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                         'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
* Chargement CAS2
  SI (V_CHA1 EGA 2);
    CH_SEI1 = MANU CHML L1 
        'C11' (10.E3*10.E3) 'C22' (200.E3*200.E3) 'C33' (10.E3*10.E3)
              'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
              'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
              'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
              'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
              'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
              'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                         'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
* Chargement CAS3
  SI (V_CHA1 EGA 3);
    CH_SEI1 = MANU CHML L1 
        'C11' (700.E3*700.E3) 'C22' (10.E3*50.E3) 'C33' (10.E3*50.E3)
              'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
              'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
              'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
              'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
              'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
              'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                       'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
FINSI;
*
*
* -------------------- Cas Poutre courte ----------------------
SI (V_TEST1 EGA 2);
* Caracteristiques
  CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
                         'LIBZ' LTOT1Z 
                       'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
                       'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
                       'ASWZ' A_CADRE1    'YACI' 210000.0E6  
                       'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POUTRE';
*
* Chargement CAS1
  SI (V_CHA1 EGA 1);
    CH_SEI1 = MANU CHML L1 
        'C11' (10.E3*10.E3) 'C22' (10.E3*10.E3) 'C33' (200.E3*200.E3)
              'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
              'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
              'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
              'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
              'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
              'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                         'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
* Chargement CAS3
  SI (V_CHA1 EGA 3);
    CH_SEI1 = MANU CHML L1 
        'C11' (700.E3*700.E3) 'C22' (10.E3*50.E3) 'C33' (10.E3*50.E3)
              'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
              'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
              'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
              'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
              'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
              'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                       'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
FINSI;
*
*
*
* -------------------- Cas Pouteau long ----------------------
SI (V_TEST1 EGA 3);
* Caracteristiques
  CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
                         'LIBZ' LTOT1Z 
                         'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
                         'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
                         'ASWZ' A_CADRE1    'YACI' 210000.0E6  
                         'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POTEAU';
*
*
* Chargement CAS1
  SI (V_CHA1 EGA 1);
    CH_SEI1 = MANU CHML L1 
        'C11' (1.E3*1.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
              'C44' (00.0*00.0) 'C55' (10.E3*10.E3) 'C66' (1.E3*1.E3)
              'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
              'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
              'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
              'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
              'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                       'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' (NBIL_YP*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
* Chargement CAS2
  SI (V_CHA1 EGA 2);
    CH_SEI1 = MANU CHML L1 
          'C11' (1.E3*1.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
          'C44' (00.0*00.0) 'C55' (1.E3*1.E3) 'C66' (10.E3*10.E3)
          'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
          'C15' (00.0*00.0) 'C16' (00.0*00.0);
     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
          'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
          'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
          'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                       'C56' (00.0*00.0)) ;
     CH_STA1 = MANU CHML L1 'EFFX' (NBIL_ZP*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
* Chargement CAS3
  SI (V_CHA1 EGA 3);
     CH_SEI1 = MANU CHML L1 
          'C11' (700.E3*700.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
          'C44' (00.0*00.0) 'C55' (1.E3*1.E3) 'C66' (1.E3*1.E3)
          'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
          'C15' (00.0*00.0) 'C16' (00.0*00.0);
     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
          'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
          'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
          'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                       'C56' (00.0*00.0)) ;
     CH_STA1 = MANU CHML L1 'EFFX' (-1600*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
FINSI;
*
*
*
* -------------------- Cas Poutre long ----------------------
SI (V_TEST1 EGA 4);
* Caracteristiques
  CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
                         'LIBZ' LTOT1Z 
                         'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
                         'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
                         'ASWZ' A_CADRE1    'YACI' 210000.0E6  
                         'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POUTRE';
*
*
* Chargement CAS1
  SI (V_CHA1 EGA 1);
    CH_SEI1 = MANU CHML L1 
          'C11' (1.E3*1.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
          'C44' (00.0*00.0) 'C55' (10.E3*10.E3) 'C66' (1.E3*1.E3)
          'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
          'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
          'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
          'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
          'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                       'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' (NBIL_YP*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
* Chargement CAS3
  SI (V_CHA1 EGA 3);
    CH_SEI1 = MANU CHML L1 
          'C11' (700.E3*700.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
          'C44' (00.0*00.0) 'C55' (1.E3*1.E3) 'C66' (1.E3*1.E3)
          'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
          'C15' (00.0*00.0) 'C16' (00.0*00.0);
    CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
          'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
          'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
          'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
                       'C56' (00.0*00.0)) ;
    CH_STA1 = MANU CHML L1 'EFFX' (-1600*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
                           'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  FINSI;
FINSI;
*
*
*
*
TAB_INP                       = TABLE;
TAB_INP.'MAILLAGE'            = L1;
TAB_INP.'EFFORT_SEISME'       = CH_SEI1;
TAB_INP.'EFFORT_STATIQUE'     = CH_STA1;
TAB_INP.'CARACTERISTIQUES'    = CH_CAR1;
TAB_INP.'FC_BET'              = 30.E6;
TAB_INP.'GAM_C'               = 1.20;
TAB_INP.'ALP_C'               = 1.00;
TAB_INP.'FS_ACI'              = 500.E6;
TAB_INP.'GAM_S'               = 1.0;
*
LLL1 = LECT 1;

TAB22 = MRCFRAME 'ELLIPSOIDE' TAB_INP 2.E-8 LLL1 ;
TAB33 = MRCFRAME 'RECTANGLE' TAB_INP 1.E-6 LLL1 ;
*
*
LAMB1 = ('EXTR' TAB22.CH_LAMBDA 'LAMB' 1 1 1);
LAMB2 = ('EXTR' TAB33.CH_LAMBDA 'LAMB' 1 1 1);
*
*-------------------------------------------------------------------
SI (V_TEST1 EGA 1);
  SI (V_CHA1 EGA 1);
    LAMB3 = (VZ_MAX1/200.);
  FINSI;
  SI (V_CHA1 EGA 2);
    LAMB3 = (VY_MAX1/200.);
  FINSI;
  SI (V_CHA1 EGA 3);
    LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  FINSI;
FINSI;
*
*-------------------------------------------------------------------
SI (V_TEST1 EGA 2);
  SI (V_CHA1 EGA 1);
    LAMB3 = (VZ_MAX1/200.);
  FINSI;
  SI (V_CHA1 EGA 3);
    LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  FINSI;
FINSI;
*
*-------------------------------------------------------------------
SI (V_TEST1 EGA 3);
  SI (V_CHA1 EGA 1);
    LAMB3 = (MBIL_YN/10.);
  FINSI;
  SI (V_CHA1 EGA 2);
    LAMB3 = (MBIL_ZN/10.);
  FINSI;
  SI (V_CHA1 EGA 3);
    LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  FINSI;
FINSI;
*
*-------------------------------------------------------------------
SI (V_TEST1 EGA 4);
  SI (V_CHA1 EGA 1);
    LAMB3 = (MBIL_YN/10.);
  FINSI;
  SI (V_CHA1 EGA 3);
    LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  FINSI;
FINSI;
*
*
*
ERR_11 = (ABS(LAMB1 - LAMB3))/LAMB3;
SI (ERR_11 > 1.E-4);
  ERRE 5;
FINSI;
SI (LAMB2 > LAMB1);
  ERRE 5;
FINSI;
LIST LAMB1;
LIST LAMB2;
LIST LAMB3;
FIN;