
DEBUT();

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

model=AFFE_MODELE(MAILLAGE=mesh,
                  AFFE=_F(TOUT='OUI',
                          PHENOMENE='MECANIQUE',
                          MODELISATION='AXIS',),);

true=DEFI_FONCTION(
                   NOM_PARA='EPSI',
                   NOM_RESU='SIGMA',
                   VALE=(0.0023 ,34410,
                         0.0058 ,42780,
                         0.0082 ,45180,
                         0.0150 ,49320,
                         0.0220 ,51850,
                         0.0303 ,45020,
                         0.0519 ,55740,
                         0.0782 ,60690,
                         0.1024 ,62670,
                         0.1031 ,62720,
                         0.2033 ,67980,
                         0.3035 ,71260,
                         0.4036 ,73680,
                         0.5037 ,75620,
                         0.6038 ,77240,
                         0.7039 ,78630,
                         0.8039 ,79860,
                         0.9040 ,80960,
                         1.0040 ,81960,
                         ),
                   INTERPOL='LIN',
                   PROL_DROITE='LINEAIRE',
                   PROL_GAUCHE='CONSTANT',);

two25=DEFI_MATERIAU(ELAS=_F(E=20410000,
                            NU=0.3,),
                    TRACTION=_F(SIGM=true,),);

material=AFFE_MATERIAU(MAILLAGE=mesh,
                       AFFE=_F(TOUT='OUI',
                               MATER=two25,),);

load=AFFE_CHAR_MECA(MODELE=model,
                    DDL_IMPO=_F(GROUP_MA='Top',
                                DY=1,),);

fix=AFFE_CHAR_MECA(MODELE=model,
                   DDL_IMPO=_F(GROUP_MA='Center',
                               DY=0.0,),);

load1=DEFI_FONCTION(
                    NOM_PARA='INST',
                    VALE=(0.0,0.0,
                          1.0,1.0,
                          ),);

step1=DEFI_LIST_REEL(DEBUT=0.0,
                     INTERVALLE=_F(JUSQU_A=1.0,
                                   NOMBRE=100,),);

step=DEFI_LIST_INST(DEFI_LIST=_F(METHODE='MANUEL',
                                 LIST_INST=step1,),);

res_nl=STAT_NON_LINE(MODELE=model,
                     CHAM_MATER=material,
                     EXCIT=(_F(CHARGE=fix,),
                            _F(CHARGE=load,
                               FONC_MULT=load1,),),
                     COMP_INCR=_F(RELATION='VMIS_ISOT_TRAC',
                                  DEFORMATION='SIMO_MIEHE',),
                     INCREMENT=_F(LIST_INST=step,),
                     NEWTON=_F(REAC_INCR=1,
                               MATRICE='TANGENTE',
                               REAC_ITER=1,),
                     CONVERGENCE=_F(RESI_GLOB_MAXI=100,
                                    RESI_GLOB_RELA=1.0e-10,),);

res_nl=CALC_CHAMP(reuse =res_nl,
                  RESULTAT=res_nl,
                  TOUT='OUI',
                  CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',),
                  DEFORMATION=('EPSP_ELNO','EPSP_NOEU',),
                  CRITERES=('SIEQ_ELNO','SIEQ_NOEU','EPEQ_ELNO','EPEQ_NOEU',),
                  FORCE=('FORC_NODA','REAC_NODA',),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=8,
          RESU=_F(RESULTAT=res_nl,
                  NOM_CHAM='REAC_NODA',
                  GROUP_MA='Top',),);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(MAILLAGE=mesh,
                  RESULTAT=res_nl,
                  TOUT_CHAM='OUI',),);

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