DEBUT(LANG='FR')

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

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

mater = DEFI_MATERIAU(ELAS=_F(E=6000.0,
                              NU=0.4))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         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', ))),
                      MODELE=model)

load0 = AFFE_CHAR_MECA(FORCE_ARETE=_F(FY=10.0,
                                      GROUP_MA=('saika', )),
                       MODELE=model)

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

unnamed = CALC_CHAMP(CHAM_MATER=fieldmat,
                     CONTRAINTE=('SIGM_NOEU', ),
                     CRITERES=('SIEQ_NOEU', ),
                     DEFORMATION=('EPSI_NOEU', ),
                     FORCE=('FORC_NODA', 'REAC_NODA'),
                     MODELE=model,
                     RESULTAT=reslin)

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

FIN()
