DEBUT(LANG='FR')


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

mesh = MODI_MAILLAGE(identifier='1:1',
                     reuse=mesh,
                     MAILLAGE=mesh,
                     ORIE_PEAU_2D=_F(GROUP_MA=('kotei', 'saika')))

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

elemprop = AFFE_CARA_ELEM(identifier='3:1',
                          COQUE=_F(COQUE_NCOU=5,
                                   EPAIS=0.1,
                                   GROUP_MA=('zentai', ),
                                   INER_ROTA='OUI'),
                          MODELE=model)

mater = DEFI_MATERIAU(identifier='4:1',
                      ELAS=_F(E=206000.0,
                              NU=0.3))

fieldmat = AFFE_MATERIAU(identifier='5:1',
                         AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MODELE=model)

load = AFFE_CHAR_MECA(identifier='6:1',
                      DDL_IMPO=_F(DX=0.0,
                                  DY=0.0,
                                  DZ=0.0,
                                  GROUP_MA=('kotei', )),
                      FORCE_ARETE=_F(FZ=-10.0,
                                     GROUP_MA=('saika', )),
                      MODELE=model)

reslin = MECA_STATIQUE(identifier='7:1',
                       CARA_ELEM=elemprop,
                       CHAM_MATER=fieldmat,
                       EXCIT=_F(CHARGE=load),
                       MODELE=model,
                       SOLVEUR=_F(NPREC=-1,
                                  RESI_RELA=140.0))

unnamed = CALC_CHAMP(identifier='8:1',
                     CARA_ELEM=elemprop,
                     CONTRAINTE=('SIGM_ELNO', ),
                     CRITERES=('SIEQ_ELNO', ),
                     DEFORMATION=('EPSI_ELNO', ),
                     RESULTAT=reslin)

unnamed2 = POST_CHAMP(identifier='9:1',
                      EXTR_COQUE=_F(NIVE_COUCHE='SUP',
                                    NOM_CHAM=('EPSI_ELNO', 'SIEQ_ELNO', 'SIGM_ELNO'),
                                    NUME_COUCHE=1),
                      RESULTAT=unnamed)

unnamed3 = CALC_CHAMP(identifier='10:1',
                      CARA_ELEM=elemprop,
                      CONTRAINTE=('SIGM_NOEU', ),
                      CRITERES=('SIEQ_NOEU', ),
                      DEFORMATION=('EPSI_NOEU', ),
                      MODELE=model,
                      RESULTAT=unnamed2)

IMPR_RESU(identifier='11:1',
          FORMAT='MED',
          RESU=(_F(CARA_ELEM=elemprop,
                   NOM_CHAM=('DEPL', ),
                   RESULTAT=reslin),
                _F(CARA_ELEM=elemprop,
                   NOM_CHAM=('SIGM_NOEU', 'EPSI_NOEU', 'SIEQ_NOEU'),
                   RESULTAT=unnamed),
                _F(CARA_ELEM=elemprop,
                   NOM_CHAM=('SIEQ_ELNO', 'EPSI_ELNO', 'SIGM_ELNO'),
                   RESULTAT=unnamed)),
          UNITE=80)

FIN()