DEBUT(LANG='EN')


mesh = LIRE_MAILLAGE(identifier='0:1',
                     FORMAT='MED',
                     UNITE=20)

modelth = AFFE_MODELE(identifier='1:1',
                      AFFE=_F(MODELISATION=('AXIS', ),
                              PHENOMENE='THERMIQUE',
                              TOUT='OUI'),
                      MAILLAGE=mesh)

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

mater1 = DEFI_MATERIAU(identifier='3:1',
                       ELAS=_F(ALPHA=1.1e-05,
                               E=200000.0,
                               NU=0.3),
                       THER=_F(LAMBDA=0.5,
                               RHO_CP=1.0))

mater2 = DEFI_MATERIAU(identifier='4:1',
                       ELAS=_F(ALPHA=2.3e-05,
                               E=70000.0,
                               NU=0.3),
                       THER=_F(LAMBDA=0.5,
                               RHO_CP=1.0))

mater3 = DEFI_MATERIAU(identifier='5:1',
                       ELAS=_F(ALPHA=0.0,
                               E=200000000.0,
                               NU=0.3),
                       THER=_F(LAMBDA=0.5,
                               RHO_CP=1.0))

matth = AFFE_MATERIAU(identifier='6:1',
                      AFFE=(_F(GROUP_MA=('innerring_Faces', ),
                               MATER=(mater1, )),
                            _F(GROUP_MA=('outering_Faces', ),
                               MATER=(mater2, )),
                            _F(GROUP_MA=('top_Faces', ),
                               MATER=(mater3, ))),
                      MAILLAGE=mesh)

linst = DEFI_LIST_REEL(identifier='7:1',
                       VALE=(0.0, 1.0))

loadst = AFFE_CHAR_THER(identifier='8:1',
                        MODELE=modelth,
                        TEMP_IMPO=_F(GROUP_MA=('innerring_Faces', 'top_Faces', 'outering_Faces'),
                                     TEMP=80.0))

mecabc = AFFE_CHAR_MECA(identifier='9:1',
                        DDL_IMPO=(_F(DY=0.0,
                                     GROUP_NO=('outer_base', 'inner_base')),
                                  _F(DX=0.0,
                                     GROUP_NO=('top_yslide', ))),
                        MODELE=modmeca)

resuth = THER_LINEAIRE(identifier='10:1',
                       CHAM_MATER=matth,
                       ETAT_INIT=_F(VALE=0.0),
                       EXCIT=_F(CHARGE=loadst),
                       INCREMENT=_F(LIST_INST=linst),
                       MODELE=modelth)

matmeca = AFFE_MATERIAU(identifier='11:1',
                        AFFE=(_F(GROUP_MA=('innerring_Faces', ),
                                 MATER=(mater1, )),
                              _F(GROUP_MA=('outering_Faces', ),
                                 MATER=(mater2, )),
                              _F(GROUP_MA=('top_Faces', ),
                                 MATER=(mater3, ))),
                        AFFE_VARC=_F(EVOL=resuth,
                                     NOM_VARC='TEMP',
                                     TOUT='OUI',
                                     VALE_REF=0.0),
                        MODELE=modmeca)

resumeca = MECA_STATIQUE(identifier='12:1',
                         CHAM_MATER=matmeca,
                         EXCIT=_F(CHARGE=mecabc),
                         LIST_INST=linst,
                         MODELE=modmeca)

IMPR_RESU(identifier='13:1',
          FORMAT='MED',
          RESU=(_F(RESULTAT=resuth),
                _F(RESULTAT=resumeca)),
          UNITE=80)

FIN()