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

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

mesh = MODI_MAILLAGE(identifier='2:1',
                     reuse=mesh,
                     MAILLAGE=mesh,
                     ORIE_PEAU_2D=(_F(GROUP_MA=('outer_fit', ),
                                      GROUP_MA_SURF=('outerring', )),
                                   _F(GROUP_MA=('inner_fit', ),
                                      GROUP_MA_SURF=('innterring', ))))

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

mater2 = DEFI_MATERIAU(identifier='4:1',
                       ELAS=_F(E=206000.0,
                               NU=0.3))

mater1 = DEFI_MATERIAU(identifier='5:1',
                       ELAS=_F(E=200000.0,
                               NU=0.3))

fieldma1 = AFFE_MATERIAU(identifier='6:1',
                         AFFE=_F(MATER=(mater1, ),
                                 TOUT='OUI'),
                         MODELE=model1)

step = DEFI_LIST_REEL(identifier='7:1',
                      DEBUT=0.0,
                      INTERVALLE=_F(JUSQU_A=1.0,
                                    NOMBRE=20))

load2 = DEFI_FONCTION(identifier='8:1',
                      NOM_PARA='INST',
                      VALE=(0.0, 0.0, 1.0, 1.0))

fix = AFFE_CHAR_MECA(identifier='9:1',
                     DDL_IMPO=_F(DY=0.0,
                                 GROUP_MA=('inner_base', 'outer_base')),
                     MODELE=model1)

contact = DEFI_CONTACT(identifier='10:1',
                       FORMULATION='CONTINUE',
                       MODELE=model1,
                       ZONE=_F(CONTACT_INIT='INTERPENETRE',
                               GROUP_MA_ESCL=('inner_fit', ),
                               GROUP_MA_MAIT=('outer_fit', )))

resnonl = STAT_NON_LINE(identifier='11:1',
                        CHAM_MATER=fieldma1,
                        COMPORTEMENT=_F(DEFORMATION='PETIT',
                                        RELATION='ELAS',
                                        TOUT='OUI'),
                        CONTACT=contact,
                        CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                       RESI_GLOB_RELA=1e-07),
                        EXCIT=_F(CHARGE=fix),
                        INCREMENT=_F(LIST_INST=step),
                        METHODE='NEWTON',
                        MODELE=model1,
                        NEWTON=_F(MATRICE='TANGENTE',
                                  REAC_INCR=1,
                                  REAC_ITER=1),
                        SOLVEUR=_F(METHODE='MUMPS'))

resnonl = CALC_CHAMP(identifier='12:1',
                     reuse=resnonl,
                     CONTRAINTE=('SIGM_ELNO', 'SIGM_NOEU', 'SIEF_NOEU'),
                     CRITERES=('SIEQ_NOEU', 'SIEQ_ELNO'),
                     RESULTAT=resnonl)

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

FIN(identifier='14:1',
    )
