Editeur/Patrons/modes_3D.comm

   1
   2 # PATRON
   3 # ANALYSE MODALE
   4 # 3D
   5
   6 DEBUT();
   7 # MISE EN DONNEES
   8
   9 MAIL=LIRE_MAILLAGE();
  10
  11 MODELE=AFFE_MODELE(MAILLAGE=MAIL,
  12                    AFFE=_F(GROUP_MA='OBJET',
  13                            PHENOMENE='MECANIQUE',
  14                            MODELISATION='3D',),);
  15 # MATERIAU
  16
  17 ACIER=DEFI_MATERIAU(ELAS=_F(E=2.1E11,
  18                             NU=0.3,
  19                             RHO=7800.0,),);
  20
  21 CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
  22                     AFFE=_F(GROUP_MA='OBJET',
  23                             MATER=ACIER,),);
  24
  25 # CONDITIONS AUX LIMITES ET CHARGEMENT
  26 # CONDITIONS AUX LIMITES
  27
  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
  32
  33 MACRO_MATR_ASSE(MODELE=MODELE,
  34                 CHAM_MATER=CHMAT,
  35                 CHARGE=BLOQUAGE,
  36                 NUME_DDL=CO('NUMEDDL'),
  37                 MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
  38                               OPTION='RIGI_MECA',),
  39                            _F(MATRICE=CO('MASSE'),
  40                               OPTION='MASS_MECA',),),);
  41 # CALCUL DES CINQ PREMIERS MODES
  42
  43 MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
  44                        MATR_B=MASSE,
  45                        CALC_FREQ=_F(
  46                        OPTION='PLUS_PETITE',
  47                        NMAX_FREQ=5,),);
  48 # ECRITURE DES MODES AU FORMAT GMSH
  49
  50 IMPR_RESU(MODELE=MODELE,
  51           FORMAT='GMSH',
  52           UNITE=19,
  53           RESU=_F(MAILLAGE=MAIL,
  54                   RESULTAT=MODES,
  55                   NOM_CHAM='DEPL',),);
  56
  57 FIN();