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',),),);
40 # CALCUL DES CINQ PREMIERS MODES
42 MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
49 EXCIT=AFFE_CHAR_MECA(MODELE=MODELE,
50 PRES_REP=_F(GROUP_MA='PRESSE',
53 EXCEL=CALC_VECT_ELEM(OPTION='CHAR_MECA',
56 EXCASS=ASSE_VECTEUR(VECT_ELEM=EXCEL,
61 SINU = FORMULE(VALE='sin(OMEGA*INST)',
65 TRADIR=DYNA_LINE_TRAN(MATR_MASS=MASSE,
68 EXCIT=_F(VECT_ASSE=EXCASS,),