# PATRON
# ANALYSE MODALE
# 3D

DEBUT();
# MISE EN DONNEES

MAIL=LIRE_MAILLAGE();

MODELE=AFFE_MODELE(MAILLAGE=MAIL,
                   AFFE=_F(GROUP_MA='OBJET',
                           PHENOMENE='MECANIQUE',
                           MODELISATION='3D',),);
# MATERIAU

ACIER=DEFI_MATERIAU(ELAS=_F(E=2.1E11,
                            NU=0.3,
                            RHO=7800.0,),);

CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
                    AFFE=_F(GROUP_MA='OBJET',
                            MATER=ACIER,),);
# CONDITIONS AUX LIMITES ET CHARGEMENT
# CONDITIONS AUX LIMITES 

BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
                        DDL_IMPO=_F(GROUP_MA='ENCAST',
                                    LIAISON='ENCASTRE',),);
# CALCUL DES MATRICES ASSEMBLEES K ET M

MACRO_MATR_ASSE(MODELE=MODELE,
                CHAM_MATER=CHMAT,
                CHARGE=BLOQUAGE,
                NUME_DDL=CO('NUMEDDL'),
                MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
                              OPTION='RIGI_MECA',),
                           _F(MATRICE=CO('MASSE'),
                              OPTION='MASS_MECA',),),);
# CALCUL DES CINQ PREMIERS MODES

MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
                       MATR_B=MASSE,
                       CALC_FREQ=_F(
                       OPTION='PLUS_PETITE',
                       NMAX_FREQ=5,),);
#CHARGEMENT

EXCIT=AFFE_CHAR_MECA(MODELE=MODELE,
                     PRES_REP=_F(GROUP_MA='PRESSE',
                                 PRES=100.e5,),);

EXCEL=CALC_VECT_ELEM(OPTION='CHAR_MECA',
                     CHARGE=EXCIT,);

EXCASS=ASSE_VECTEUR(VECT_ELEM=EXCEL,
                    NUME_DDL=NUMEDDL,);
OMEGA = 10.0;


SINU = FORMULE(VALE='sin(OMEGA*INST)',
               NOM_PARA='INST',);
#CALCUL TRANSITOIRE

TRADIR=DYNA_LINE_TRAN(MATR_MASS=MASSE,
                      MATR_RIGI=RIGIDITE,
                      NEWMARK=_F(),
                      EXCIT=_F(VECT_ASSE=EXCASS,),
                      INCREMENT=_F(PAS=.01,
                                   INST_FIN=10.,),);

FIN();