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',
25 # CONDITIONS AUX LIMITES ET CHARGEMENT
26 # CONDITIONS AUX LIMITES
28 BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
29 DDL_IMPO=_F(GROUP_MA='ENCAST',
30 LIAISON='ENCASTRE',),);
31 # CALCUL DES MATRICES ASSEMBLEES K ET M
33 MACRO_MATR_ASSE(MODELE=MODELE,
36 NUME_DDL=CO('NUMEDDL'),
37 MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
39 _F(MATRICE=CO('MASSE'),
40 OPTION='MASS_MECA',),),);
41 # CALCUL DES CINQ PREMIERS MODES
43 MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
48 # ECRITURE DES MODES AU FORMAT GMSH
50 IMPR_RESU(MODELE=MODELE,
53 RESU=_F(MAILLAGE=MAIL,