
DEBUT();

MA=DEFI_MATERIAU(ELAS=_F(E=210000000000.0,
                         NU=0.3,),
                 ECRO_LINE=_F(D_SIGM_EPSI=3.283e+9,
                              SY=155.8e+6,),);

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

MAIL=MODI_MAILLAGE(reuse =MAIL,
                   MAILLAGE=MAIL,
                   ORIE_PEAU_3D=_F(GROUP_MA='load',),);

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='hidari',
                                 DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,),
                              _F(GROUP_MA='migi',
                                 DX=0.0,
                                 DY=0.0,),),);

loadP=AFFE_CHAR_MECA(MODELE=MODE,
                     FORCE_FACE=_F(GROUP_MA='load',
                                   FX=1600000000,),);

ramp=DEFI_FONCTION(
                   NOM_PARA='INST',
                   VALE=(0.0,0.0,
                         1.0,1.0,
                         ),
                   INTERPOL='LIN',);

inst=DEFI_LIST_REEL(DEBUT=0.0,
                    INTERVALLE=_F(JUSQU_A=1.0,
                                  NOMBRE=40,),);

RESU=STAT_NON_LINE(MODELE=MODE,
                   CHAM_MATER=MATE,
                   EXCIT=(_F(CHARGE=CHAR,),
                          _F(CHARGE=loadP,
                             FONC_MULT=ramp,),),
                   COMP_INCR=_F(RELATION='VMIS_ISOT_LINE',
                                TOUT='OUI',),
                   INCREMENT=_F(LIST_INST=inst,),
                   METHODE='IMPL_EX',
                   CONVERGENCE=_F(RESI_GLOB_RELA=1e-6,
                                  ITER_GLOB_MAXI=10,),
                   ARCHIVAGE=_F(LIST_INST=inst,
                                CHAM_EXCLU="""VARI_ELGA""",),);

RESU=CALC_ELEM(reuse =RESU,
               MODELE=MODE,
               CHAM_MATER=MATE,
               RESULTAT=RESU,
               OPTION=('SIEF_ELNO','EPSI_ELNO','SIEQ_ELNO',),);

reac1=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
                              INTITULE='DEPL',
                              RESULTAT=RESU,
                              NOM_CHAM='DEPL',
                              GROUP_NO='load',
                              RESULTANTE='DX',),);

IMPR_TABLE(TABLE=reac1,);

RESU=CALC_NO(reuse =RESU,
             RESULTAT=RESU,
             OPTION=('SIGM_NOEU','SIEQ_NOEU',),);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=(_F(MAILLAGE=MAIL,
                   RESULTAT=RESU,
                   NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),
                   NOM_CMP=('DX','DY','DZ',),),
                _F(MAILLAGE=MAIL,
                   RESULTAT=RESU,
                   NOM_CHAM=('SIEQ_NOEU','SIEF_NOEU',),),
                _F(MAILLAGE=MAIL,
                   RESULTAT=RESU,
                   NOM_CHAM='EPSG_NOEU',),),);

FIN();
#CHECKSUM:d8acfcef4885d1b1f9225d0a8440e14a  -:FIN CHECKSUM
