DEBUT(identifier=u'0:1',
      LANG='EN')

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

model1 = AFFE_MODELE(identifier=u'2:1',
                     AFFE=(_F(GROUP_MA=('GrMesh_2_Volumes', ),
                              MODELISATION=('3D', ),
                              PHENOMENE='MECANIQUE'),
                           _F(GROUP_MA=('GrMesh_3_Volumes', ),
                              MODELISATION=('3D', ),
                              PHENOMENE='MECANIQUE')),
                     MAILLAGE=mesh)

mater2 = DEFI_MATERIAU(identifier=u'3:1',
                       ELAS=_F(E=21000000000.0,
                               NU=0.3))

mater1 = DEFI_MATERIAU(identifier=u'4:1',
                       ELAS=_F(E=11000000000.0,
                               NU=0.3))

fieldma1 = AFFE_MATERIAU(identifier=u'5:1',
                         AFFE=(_F(GROUP_MA=('GrMesh_3_Volumes', ),
                                  MATER=(mater2, )),
                               _F(GROUP_MA=('GrMesh_2_Volumes', ),
                                  MATER=(mater1, ))),
                         MODELE=model1)

fix = AFFE_CHAR_MECA(identifier=u'6:1',
                     DDL_IMPO=_F(DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,
                                 GROUP_MA=('base', )),
                     MODELE=model1)

load = AFFE_CHAR_MECA(identifier=u'7:1',
                      DDL_IMPO=_F(DZ=-0.001,
                                  GROUP_MA=('load', )),
                      MODELE=model1)

liason = AFFE_CHAR_MECA(identifier=u'8:1',
                        LIAISON_MAIL=_F(GROUP_MA_ESCL=('fit_b', ),
                                        GROUP_MA_MAIT=('GrMesh_3_Volumes', ),
                                        TYPE_RACCORD='MASSIF'),
                        MODELE=model1)

reslin = MECA_STATIQUE(identifier=u'9:1',
                       CHAM_MATER=fieldma1,
                       EXCIT=(_F(CHARGE=load),
                              _F(CHARGE=fix),
                              _F(CHARGE=liason)),
                       MODELE=model1)

reslin = CALC_CHAMP(identifier=u'10:1',
                    reuse=reslin,
                    CONTRAINTE=('SIGM_NOEU', ),
                    CRITERES=('SIEQ_NOEU', ),
                    DEFORMATION=('EPSI_NOEU', ),
                    FORCE=('FORC_NODA', 'REAC_NODA'),
                    RESULTAT=reslin)

IMPR_RESU(identifier=u'11:1',
          FORMAT='MED',
          RESU=_F(MAILLAGE=mesh,
                  NOM_CHAM=('EPSI_NOEU', 'SIEQ_NOEU', 'SIGM_NOEU', 'DEPL'),
                  RESULTAT=reslin),
          UNITE=3)

FIN(identifier=u'12:1',
    )
