
DEBUT();

sig400=DEFI_FONCTION(
                     NOM_PARA='EPSI',
                     NOM_RESU='SIGM',
                     VALE=(0.00995 ,40.4 ,
                           0.0198 ,51.0 ,
                           0.02956 ,61.8 ,
                           0.03922 ,72.8 ,
                           0.04879 ,84.0 ,),
                     PROL_DROITE='LINEAIRE',
                     PROL_GAUCHE='LINEAIRE',);

MA=DEFI_MATERIAU(ELAS=_F(E=4000.0,
                         NU=0.3,),
                 TRACTION=_F(SIGM=sig400,),
                 PRAGER=_F(C=1000,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
                 AFFE=_F(TOUT='OUI',
                         PHENOMENE='MECANIQUE',
                         MODELISATION='3D',),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
                   AFFE=_F(TOUT='OUI',
                           MATER=MA,),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=(_F(GROUP_MA='Group_4',
                                 DX=0.0,
                                 DZ=0.0,),
                              _F(GROUP_MA='Group_5',
                                 DX=0.0,
                                 DY=0.0,),
                              _F(GROUP_MA='Group_2',
                                 DX=0.0,),
                              _F(GROUP_MA='Group_1',
                                 DX=6.0,),),);

inst=DEFI_LIST_REEL(DEBUT=0.0,
                    INTERVALLE=(_F(JUSQU_A=10.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=20.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=30.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=40.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=50.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=60.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=70.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=80.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=90.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=100.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=110.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=120.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=130.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=140.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=150.0,
                                   NOMBRE=12,),
                                _F(JUSQU_A=160.0,
                                   NOMBRE=12,),),);

fonc1=DEFI_FONCTION(
                    NOM_PARA='INST',
                    VALE=(0.0 ,0.0 ,
                          10.0 ,1.0 ,
                          20.0 ,0.0 ,
                          30.0 ,-1.0 ,
                          40.0 ,0.0 ,
                          50.0 ,1.0 ,
                          60.0 ,0.0 ,
                          70.0 ,-1.0 ,
                          80.0 ,0.0 ,
                          90.0 ,1.0 ,
                          100.0 ,0.0 ,
                          110.0 ,-1.0 ,
                          120.0 ,0.0 ,
                          130.0 ,1.0 ,
                          140.0 ,0.0 ,
                          150.0 ,-1.0 ,
                          160.0 ,0.0 ,),);

RESU=STAT_NON_LINE(MODELE=MODE,
                   CHAM_MATER=MATE,
                   EXCIT=_F(CHARGE=CHAR,
                            FONC_MULT=fonc1,),
                   COMPORTEMENT=_F(RELATION='VMIS_ECMI_TRAC',
                                   DEFORMATION='PETIT_REAC',),
                   INCREMENT=_F(LIST_INST=inst,),
                   NEWTON=_F(PREDICTION='TANGENTE',
                             MATRICE='TANGENTE',
                             REAC_ITER=1,),
                   CONVERGENCE=_F(RESI_GLOB_RELA=0.01,
                                  ITER_GLOB_MAXI=200,),
                   ARCHIVAGE=_F(LIST_INST=inst,
                                CHAM_EXCLU='VARI_ELGA',),);

RESU=CALC_CHAMP(reuse =RESU,
                MODELE=MODE,
                CHAM_MATER=MATE,
                RESULTAT=RESU,
                EXCIT=_F(CHARGE=CHAR,),
                CONTRAINTE='SIEF_ELNO',
                CRITERES='SIEQ_ELNO',);

RESU=CALC_CHAMP(reuse =RESU,
                RESULTAT=RESU,
                CONTRAINTE=('SIGM_NOEU','SIEF_NOEU',),
                CRITERES='SIEQ_NOEU',);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(MAILLAGE=MAIL,
                  RESULTAT=RESU,
                  NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();
#CHECKSUM:01dafe22424e902cb44b92e6c2656de2  -:FIN CHECKSUM
