DEBUT(LANG='EN')

mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)

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

kou = DEFI_MATERIAU(ELAS=_F(E=206000.0,
                            NU=0.3))

moku = DEFI_MATERIAU(ELAS_ORTH=_F(E_L=240.0,
                                  E_N=6000.0,
                                  E_T=240.0,
                                  G_LN=400.0,
                                  G_LT=400.0,
                                  G_TN=400.0,
                                  NU_LN=0.016,
                                  NU_LT=0.016,
                                  NU_TN=0.016))

fieldmat = AFFE_MATERIAU(AFFE=(_F(GROUP_MA=('kou', ),
                                  MATER=(kou, )),
                               _F(GROUP_MA=('moku', ),
                                  MATER=(moku, ))),
                         MODELE=model)

load = AFFE_CHAR_MECA(DDL_IMPO=(_F(DX=0.0,
                                   DY=0.0,
                                   DZ=0.0,
                                   GROUP_MA=('kotei1', )),
                                _F(DX=0.0,
                                   DY=0.0,
                                   GROUP_MA=('kotei2', ))),
                      FORCE_ARETE=_F(FY=-10.0,
                                     GROUP_MA=('saika', )),
                      MODELE=model)

reslin = MECA_STATIQUE(CHAM_MATER=fieldmat,
                       EXCIT=_F(CHARGE=load),
                       MODELE=model)

unnamed = CALC_CHAMP(CONTRAINTE=('SIGM_NOEU', ),
                     CRITERES=('SIEQ_NOEU', ),
                     DEFORMATION=('EPSI_NOEU', ),
                     FORCE=('REAC_NODA', ),
                     RESULTAT=reslin)

IMPR_RESU(RESU=(_F(NOM_CHAM=('DEPL', ),
                   RESULTAT=reslin),
                _F(NOM_CHAM=('EPSI_NOEU', 'SIEQ_NOEU', 'SIGM_NOEU', 'REAC_NODA'),
                   RESULTAT=unnamed)),
          UNITE=80)

FIN()
