DEBUT(LANG='EN')


mesh = LIRE_MAILLAGE(identifier='0:1',
                     FORMAT='MED',
                     UNITE=2)

model = AFFE_MODELE(identifier='1:1',
                    AFFE=_F(MODELISATION=('3D', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

elast_pl = DEFI_FONCTION(identifier='2:1',
                         INTERPOL=('LIN', ),
                         NOM_PARA='EPSI',
                         NOM_RESU='TOUTRESU',
                         PROL_DROITE='LINEAIRE',
                         PROL_GAUCHE='CONSTANT',
                         VALE=(0.0015, 105.0, 0.05, 200.0, 0.2, 300.0))

mater = DEFI_MATERIAU(identifier='3:1',
                      ELAS=_F(E=70000.0,
                              NU=0.3),
                      TRACTION=_F(SIGM=elast_pl))

fieldmat = AFFE_MATERIAU(identifier='4:1',
                         AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh)

depl_imp = DEFI_FONCTION(identifier='5:1',
                         NOM_PARA='INST',
                         VALE=(0.0, 0.0, 1.0, 1.0))

listr = DEFI_LIST_REEL(identifier='6:1',
                       DEBUT=0.0,
                       INTERVALLE=_F(JUSQU_A=1.0,
                                     PAS=0.1))

times = DEFI_LIST_INST(identifier='7:1',
                       DEFI_LIST=_F(LIST_INST=listr),
                       METHODE='MANUEL')

load = AFFE_CHAR_MECA(identifier='8:1',
                      FORCE_FACE=_F(FY=10.0,
                                    GROUP_MA=('saika', )),
                      MODELE=model)

fix = AFFE_CHAR_MECA(identifier='9:1',
                     DDL_IMPO=_F(DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,
                                 GROUP_MA=('kotei', )),
                     MODELE=model)

resnonl = STAT_NON_LINE(identifier='10:1',
                        CHAM_MATER=fieldmat,
                        COMPORTEMENT=_F(DEFORMATION='SIMO_MIEHE',
                                        RELATION='VMIS_ISOT_TRAC'),
                        EXCIT=(_F(CHARGE=load,
                                  FONC_MULT=depl_imp),
                               _F(CHARGE=fix)),
                        INCREMENT=_F(LIST_INST=listr),
                        METHODE='NEWTON',
                        MODELE=model,
                        NEWTON=_F(MATRICE='TANGENTE',
                                  REAC_ITER=1,
                                  REAC_ITER_ELAS=1))

unnamed = CALC_CHAMP(identifier='11:1',
                     CONTRAINTE=('SIGM_NOEU', ),
                     CRITERES=('SIEQ_NOEU', ),
                     DEFORMATION=('EPSI_NOEU', ),
                     FORCE=('FORC_NODA', ),
                     MODELE=model,
                     RESULTAT=resnonl)

IMPR_RESU(identifier='12:1',
          FORMAT='MED',
          RESU=(_F(MAILLAGE=mesh,
                   NOM_CHAM=('DEPL', ),
                   RESULTAT=resnonl),
                _F(MAILLAGE=mesh,
                   NOM_CHAM=('EPSI_NOEU', 'SIGM_NOEU', 'SIEQ_NOEU'),
                   RESULTAT=unnamed)),
          UNITE=80)

FIN()