DEBUT(LANG='EN')

mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=2)

model = AFFE_MODELE(AFFE=(_F(MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE',
                             TOUT='OUI'),
                          _F(GROUP_MA=('ta4', 'ta5', 'ta6', 'ta1', 'ta2', 'ta3'),
                             MODELISATION=('POU_D_E', ),
                             PHENOMENE='MECANIQUE'),
                          _F(GROUP_MA=('sss', ),
                             MODELISATION=('DKT', ),
                             PHENOMENE='MECANIQUE')),
                    MAILLAGE=mesh)

elemprop = AFFE_CARA_ELEM(COQUE=_F(ANGL_REP=(0.0, 0.0),
                                   COQUE_NCOU=1,
                                   EPAIS=8.0,
                                   GROUP_MA=('sss', ),
                                   INER_ROTA='OUI'),
                          MODELE=model,
                          POUTRE=_F(CARA=('A', 'IY', 'IZ', 'JX'),
                                    GROUP_MA=('ta4', 'ta5', 'ta6', 'ta1', 'ta2', 'ta3'),
                                    SECTION='GENERALE',
                                    VALE=(100.0, 833.0, 833.0, 1347.87),
                                    VARI_SECT='CONSTANT'))

mater = DEFI_MATERIAU(ELAS=_F(E=7500.0,
                              NU=0.4,
                              RHO=0.0459))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MODELE=model)

load = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
                                  DY=0.0,
                                  DZ=0.0,
                                  GROUP_MA=('kotei', )),
                      LIAISON_SOLIDE=(_F(GROUP_NO=('ke1', )),
                                      _F(GROUP_NO=('ke2', )),
                                      _F(GROUP_NO=('ke3', )),
                                      _F(GROUP_NO=('ke4', )),
                                      _F(GROUP_NO=('ke5', )),
                                      _F(GROUP_NO=('ke6', )),
                                      _F(GROUP_NO=('ke7', )),
                                      _F(GROUP_NO=('ke8', ))),
                      MODELE=model)

ASSEMBLAGE(CARA_ELEM=elemprop,
           CHAM_MATER=fieldmat,
           CHARGE=(load, ),
           MATR_ASSE=(_F(MATRICE=CO('mass'),
                         OPTION='MASS_MECA'),
                      _F(MATRICE=CO('stifness'),
                         OPTION='RIGI_MECA')),
           MODELE=model,
           NUME_DDL=CO('nddl'))

modes = CALC_MODES(CALC_FREQ=_F(NMAX_FREQ=10),
                   MATR_MASS=mass,
                   MATR_RIGI=stifness,
                   OPTION='PLUS_PETITE',
                   VERI_MODE=_F(STOP_ERREUR='NON'))

IMPR_RESU(FORMAT='MED',
          RESU=_F(NOM_CHAM=('DEPL', ),
                  RESULTAT=modes),
          UNITE=3)

FIN()
