POURSUITE(identifier='0:2',
          LANG='EN')

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

mater1 = DEFI_MATERIAU(identifier='2:2',
                       ELAS=_F(ALPHA=1.1e-05,
                               E=200000.0,
                               NU=0.3))

mater3 = DEFI_MATERIAU(identifier='3:2',
                       ELAS=_F(ALPHA=0.0,
                               E=200000000.0,
                               NU=0.3))

mater2 = DEFI_MATERIAU(identifier='4:2',
                       ELAS=_F(ALPHA=2.3e-05,
                               E=70000.0,
                               NU=0.3))

fieldma1 = AFFE_MATERIAU(identifier='5:2',
                         AFFE=(_F(GROUP_MA=('innerring_Faces', ),
                                  MATER=(mater1, )),
                               _F(GROUP_MA=('top_Faces', ),
                                  MATER=(mater3, )),
                               _F(GROUP_MA=('outering_Faces', ),
                                  MATER=(mater2, ))),
                         AFFE_VARC=_F(EVOL=thload,
                                      NOM_VARC='TEMP',
                                      TOUT='OUI',
                                      VALE_REF=20.0),
                         MODELE=model1)

fix = AFFE_CHAR_MECA(identifier='6:2',
                     DDL_IMPO=(_F(DY=0.0,
                                  GROUP_NO=('outer_base', 'inner_base')),
                               _F(DX=0.0,
                                  GROUP_NO=('top_yslide', ))),
                     MODELE=model1)

reslin = MECA_STATIQUE(identifier='7:2',
                       CHAM_MATER=fieldma1,
                       EXCIT=_F(CHARGE=fix),
                       MODELE=model1)

reslin = CALC_CHAMP(identifier='8:2',
                    reuse=reslin,
                    CONTRAINTE=('SIGM_NOEU', 'SIEF_NOEU'),
                    CRITERES=('SIEQ_NOEU', ),
                    DEFORMATION=('EPSI_NOEU', ),
                    FORCE=('FORC_NODA', 'REAC_NODA'),
                    RESULTAT=reslin)

IMPR_RESU(identifier='9:2',
          FORMAT='MED',
          RESU=_F(NOM_CHAM=('EPSI_NOEU', 'SIEQ_NOEU', 'SIGM_NOEU', 'DEPL', 'SIEF_NOEU'),
                  RESULTAT=reslin),
          UNITE=80)

FIN(identifier='10:2',
    )
