DEBUT()
YOUNG = DEFI_PARA_SENSI(VALE=1.)
FORC = DEFI_PARA_SENSI(VALE=1.)
MA=LIRE_MAILLAGE()
MA=DEFI_GROUP( reuse=MA,   MAILLAGE=MA,
                           CREA_GROUP_NO=_F(
                         GROUP_MA = ( 'LP9P10',  'LP9P10B',
                                      'LP8P9',   'LP8P9B',
                                         ))
                       )

MO=AFFE_MODELE( MAILLAGE=MA,
                             AFFE=_F(  TOUT = 'OUI',
                                    PHENOMENE = 'MECANIQUE',
                                    MODELISATION = 'C_PLAN')
                       )
ZERO =    DEFI_CONSTANTE(VALE=0.)
NU =    DEFI_CONSTANTE(VALE=0.3)
ALPHA =    DEFI_CONSTANTE(VALE=0.)

MAT=DEFI_MATERIAU(ELAS_FO=_F(  E = YOUNG, NU = NU,
                                TEMP_DEF_ALPHA =20.,
                               ALPHA = ALPHA, RHO=7800.,))

CHMAT=AFFE_MATERIAU(    MAILLAGE=MA,
                                  AFFE=_F(  TOUT = 'OUI',
                                         MATER = MAT)
                          )

CH=AFFE_CHAR_MECA_F(           MODELE=MO,
                       FORCE_INTERNE=_F(TOUT='OUI',FY=FORC,),
                               DDL_IMPO=(
                                    _F(  GROUP_NO = 'LP8P9',
                                               DX = ZERO),
                                             _F(  GROUP_NO = 'LP8P9B',
                                               DX = ZERO),
                                             _F(  GROUP_NO = 'P11',
                                               DY = ZERO))
                            )

DEP=MECA_STATIQUE(   MODELE=MO,   CHAM_MATER=CHMAT,
                      SENSIBILITE = (FORC),
                              EXCIT=_F( CHARGE = CH),
                     )

DEP=CALC_ELEM(reuse=DEP,RESULTAT=DEP,
                          SENSIBILITE = (FORC),
              OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)

FOND=DEFI_FOND_FISS(    MAILLAGE=MA,
                              FOND_FISS=_F( GROUP_NO = ('P0',)),
                           NORMALE=(0., 1., 0.,)
                      )

THETA0=CALC_THETA(       MODELE=MO,
                         THETA_2D=_F(  GROUP_NO = ('P0',),
                                    MODULE = 1.,
                                    R_INF = 3.75E-5,
                                    R_SUP = 7.50E-5),
                        DIRECTION=(1., 0., 0.,)
                     )


FIN()
MA1=LIRE_MAILLAGE()
MO1=AFFE_MODELE( MAILLAGE=MA1, AFFE=_F(  TOUT = 'OUI',
                 PHENOMENE = 'MECANIQUE', MODELISATION = 'C_PLAN'))

FIN()