11 MODELE=AFFE_MODELE(MAILLAGE=MAIL,
12 AFFE=_F(GROUP_MA='OBJET',
13 PHENOMENE='MECANIQUE',
14 MODELISATION='3D',),);
17 ACIER=DEFI_MATERIAU(ELAS=_F(E=2.1E11,
21 CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
22 AFFE=_F(GROUP_MA='OBJET',
24 # CONDITIONS AUX LIMITES ET CHARGEMENT
25 # CONDITIONS AUX LIMITES
27 BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
28 DDL_IMPO=_F(GROUP_MA='ENCAST',
29 LIAISON='ENCASTRE',),);
30 # CALCUL DES MATRICES ASSEMBLEES K ET M
32 MACRO_MATR_ASSE(MODELE=MODELE,
35 NUME_DDL=CO('NUMEDDL'),
36 MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
38 _F(MATRICE=CO('MASSE'),
39 OPTION='MASS_MECA',),),);
41 P=POST_ELEM(MASS_INER=_F(TOUT='OUI',
43 # CALCUL DES CINQ PREMIERS MODES
45 MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
52 EXCIT=AFFE_CHAR_MECA(MODELE=MODELE,
53 PRES_REP=_F(GROUP_MA='PRESSE',
56 EXCEL=CALC_VECT_ELEM(OPTION='CHAR_MECA',
59 EXCASS=ASSE_VECTEUR(VECT_ELEM=EXCEL,
64 SINU = FORMULE(VALE='sin(OMEGA*INST)',
68 MACRO_PROJ_BASE(BASE=MODES,
69 MATR_ASSE_GENE=(_F(MATRICE=CO('MASSPROJ'),
71 _F(MATRICE=CO('RIGIPROJ'),
72 MATR_ASSE=RIGIDITE,),),
73 VECT_ASSE_GENE=_F(VECTEUR=CO('EXCPROJ'),
77 TRAMOD=DYNA_TRAN_MODAL(MASS_GENE=MASSPROJ,
79 AMOR_REDUIT=(.02,.01,),
80 INCREMENT=_F(INST_INIT=0.0,
83 EXCIT=_F(VECT_GENE=EXCPROJ,