Editeur/Patrons/modes_pou.comm

   1
   2 # PATRON
   3 # ANALYSE MODALE
   4 # POUTRE
   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 CARELEM=AFFE_CARA_ELEM(MODELE=MODELE,
  26                        POUTRE=_F(GROUP_MA='OBJET',
  27                                  SECTION='CERCLE',
  28                                  CARA=('R','EP',),
  29                                  VALE=(1.0,0.01,),),);
  30 # CONDITIONS AUX LIMITES
  31
  32 BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
  33                         DDL_IMPO=_F(GROUP_MA='ENCAST',
  34                                     LIAISON='ENCASTRE',),);
  35 #MATRICES ASSEMBLEES K, M
  36 #
  37
  38 MACRO_MATR_ASSE(MODELE=MODELE,
  39                 CHAM_MATER=CHMAT,
  40                 CARA_ELEM=CARELEM,
  41                 CHARGE=BLOQUAGE,
  42                 NUME_DDL=CO('NUMEDDL'),
  43                 MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
  44                               OPTION='RIGI_MECA',),
  45                            _F(MATRICE=CO('MASSE'),
  46                               OPTION='MASS_MECA',),),);
  47 #5 PREMIERS MODES
  48
  49 MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
  50                        MATR_B=MASSE,
  51                        CALC_FREQ=_F(
  52                        OPTION='PLUS_PETITE',
  53                        NMAX_FREQ=5,),);
  54 # MODES AU FORMAT GMSH
  55
  56 IMPR_RESU(MODELE=MODELE,
  57           FORMAT='GMSH',
  58           UNITE=19,
  59           RESU=_F(MAILLAGE=MAIL,
  60                   RESULTAT=MODES,
  61                   NOM_CHAM='DEPL',),);
  62
  63 FIN();