DEBUT(LANG='EN')

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

model = AFFE_MODELE(AFFE=(_F(GROUP_MA=('sessyokub2', 'sessyokub', 'men'),
                             MODELISATION=('DKT', ),
                             PHENOMENE='MECANIQUE'),
                          _F(GROUP_MA=('Line_8', 'Line_5', 'Line_7', 'Line_6'),
                             MODELISATION=('POU_D_E', ),
                             PHENOMENE='MECANIQUE')),
                    MAILLAGE=mesh)

elemprop = AFFE_CARA_ELEM(COQUE=_F(ANGL_REP=(0.0, 0.0),
                                   COQUE_NCOU=1,
                                   EPAIS=10.0,
                                   GROUP_MA=('sessyokub2', 'sessyokub', 'men'),
                                   INER_ROTA='OUI'),
                          MODELE=model,
                          POUTRE=_F(CARA=('A', 'IY', 'IZ', 'EY', 'EZ', 'JX', 'RY', 'RZ'),
                                    GROUP_MA=('Line_8', 'Line_5', 'Line_7', 'Line_6'),
                                    SECTION='GENERALE',
                                    VALE=(100.0, 833.3, 833.3, 0.000217, -0.000206, 1347.87, 5.0, 5.0),
                                    VARI_SECT='CONSTANT'))

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

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

load = AFFE_CHAR_MECA(DDL_IMPO=_F(DRX=0.0,
                                  DRY=0.0,
                                  DRZ=0.0,
                                  DX=0.0,
                                  DY=0.0,
                                  DZ=0.0,
                                  GROUP_NO=('kotei', )),
                      LIAISON_ELEM=(_F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                       CARA_ELEM=elemprop,
                                       GROUP_MA_1=('sessyokub', ),
                                       GROUP_NO_2=('setti4', ),
                                       OPTION='COQ_POU'),
                                    _F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                       CARA_ELEM=elemprop,
                                       GROUP_MA_1=('sessyokub', ),
                                       GROUP_NO_2=('setti5', ),
                                       OPTION='COQ_POU'),
                                    _F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                       CARA_ELEM=elemprop,
                                       GROUP_MA_1=('sessyokub2', ),
                                       GROUP_NO_2=('setti6', ),
                                       OPTION='COQ_POU'),
                                    _F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                       CARA_ELEM=elemprop,
                                       GROUP_MA_1=('sessyokub2', ),
                                       GROUP_NO_2=('setti7', ),
                                       OPTION='COQ_POU')),
                      MODELE=model)

load0 = AFFE_CHAR_MECA(LIAISON_ELEM=(_F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                        CARA_ELEM=elemprop,
                                        GROUP_MA_1=('sessyokub', ),
                                        GROUP_NO_2=('setti4', ),
                                        OPTION='COQ_POU'),
                                     _F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                        CARA_ELEM=elemprop,
                                        GROUP_MA_1=('sessyokub', ),
                                        GROUP_NO_2=('setti5', ),
                                        OPTION='COQ_POU'),
                                     _F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                        CARA_ELEM=elemprop,
                                        GROUP_MA_1=('sessyokub2', ),
                                        GROUP_NO_2=('setti6', ),
                                        OPTION='COQ_POU'),
                                     _F(AXE_POUTRE=(0.0, 1.0, 0.0),
                                        CARA_ELEM=elemprop,
                                        GROUP_MA_1=('sessyokub2', ),
                                        GROUP_NO_2=('setti7', ),
                                        OPTION='COQ_POU')),
                       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=80)

FIN()
