DEBUT(LANG='EN')

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

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

elemprop = AFFE_CARA_ELEM(MASSIF=_F(ANGL_AXE=(0.0, 90.0),
                                    GROUP_MA=('Partition_3', ),
                                    ORIG_AXE=(0.0, 14500.0, 10000.0)),
                          MODELE=model)

mater = DEFI_MATERIAU(ELAS_ORTH=_F(E_L=375.0,
                                   E_N=7500.0,
                                   E_T=375.0,
                                   G_LN=500.0,
                                   G_LT=500.0,
                                   G_TN=500.0,
                                   NU_LN=0.016,
                                   NU_LT=0.016,
                                   NU_TN=0.016))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh,
                         MODELE=model)

load = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
                                  DY=0.0,
                                  DZ=0.0,
                                  GROUP_MA=('kotei', )),
                      FORCE_ARETE=_F(FY=50.0,
                                     GROUP_MA=('saica', )),
                      MODELE=model)

reslin = MECA_STATIQUE(CARA_ELEM=elemprop,
                       CHAM_MATER=fieldmat,
                       EXCIT=_F(CHARGE=load),
                       MODELE=model,
                       SOLVEUR=_F(METHODE='MUMPS',
                                  RESI_RELA=0.01))

unnamed = CALC_CHAMP(CARA_ELEM=elemprop,
                     CHAM_MATER=fieldmat,
                     CONTRAINTE=('SIGM_NOEU', ),
                     CRITERES=('SIEQ_NOEU', ),
                     DEFORMATION=('EPSI_NOEU', ),
                     MODELE=model,
                     RESULTAT=reslin)

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

FIN()
