Bsp. *comm für harmonische Analyse 3D Solid Modell

Date : 13/06/2012 Page : 9/10
Clé : U4.53.11
Révision : 9108
Example of use in imposed force
# -------------------------------------------------------------
#REFERENCE
: SFM/VPCS SDLD21
COMMANDS
# -------------------------------------------------------------
#SYSTEME
MASSE-RESSORT A 8 DEGREES OF FREEDOM AVEC DAMPING
#VISQUEUX PROPORTIONAL (REPONSE HARMONIC)
#9 SPRINGS 8 MASSES -- K = 1.E+5
M = 10
C = 50--
# -------------------------------------------------------------
DEBUT (CODE = _F (NOM = “SDLD21A”),);
# -------------------------------------------------------------
MY = LIRE_MAILLAGE ();
MO = AFFE_MODELE (
MESH = MY,
AFFE = ( _F (
_F (TOUT=
“OUI”,
PHENOMENE=
“MECANIQUE”,
MODELISATION= “DIS_T”,),
GROUP_NO= MASSES,
MODELISATION= “DIS_T”,),
),);
CARELEM = AFFE_CARA_ELEM ( MODELE = MO,
DISCRET = ( _F (
CARA=
VALE=
_F (CARA=
VALE=
_F (CARA=
VALE=
GROUP_MA= ` SPRING',
“K_T_D_L”,
(1.E+5,1., 1.),),
GROUP_NO= MASSES,
“M_T_D_N”,
10. ,),
GROUP_MA= DEADENED,
“A_T_D_L”,
(50. , 1. , 1.),),
),),);
CH = AFFE_CHAR_MECA (
MODELE = MO,
DDL_IMPO =
( _F (GROUP_NO = ` A_ET_B',
DX = 0. , DY = 0. , DZ = 0. ,),
_F (GROUP_NO = ` MASS',
DY = 0. , DZ = 0. ,),),
FORCE_NODALE =_F ( NODE = ` P4', FX = 1. ,),
);
MELR = CALC_MATR_ELEM (MODELE = MO, LOAD = CH,
OPTION = “RIGI_MECA”,
CARA_ELEM = CARELEM,);
MELM = CALC_MATR_ELEM (MODELE = MO, CHARGE= CH,
OPTION = “MASS_MECA”,
CARA_ELEM = CARELEM,);
MELC = CALC_MATR_ELEM (MODELE = MO, LOAD = CH,
OPTION = “AMOR_MECA”,
CARA_ELEM = CARELEM,);
VECT = CALC_VECT_ELEM (LOAD = CH, OPTION = “CHAR_MECA”,);
# -------------------------------------------------------------
NUM = NUME_DDL ( MATR_RIGI = MELR,);
MATASSR = ASSE_MATRICE (MATR_ELEM = MELR, NUME_DDL = NUM,);
MATASSM = ASSE_MATRICE (MATR_ELEM = MELM, NUME_DDL = NUM,);
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
default
Code_Aster
Titre : Opérateur DYNA_LINE_HARM
Responsable : Harinaivo ANDRIAMBOLOLONA
Date : 13/06/2012 Page : 10/10
Clé : U4.53.11
Révision : 9108
MATASSC = ASSE_MATRICE (MATR_ELEM = MELC, NUME_DDL = NUM,);
VECTASS = ASSE_VECTEUR (VECT_ELEM = VECT, NUME_DDL = NUM,);
# -------------------------------------------------------------
# EXCITATION FORCE SINUSOIDALE Of AMPLITUDE CRETE FX = 1.N
# WITH the POINT P4
BANDE OF FREQUENCIES 5 Hz - 40 Hz
# (AMPLITUDE INDEPENDANTE DE LA FREQUENCY)
LIFREQ = DEFI_LIST_REEL ( DEBUT = 5. ,
INTERVALLE =_F
(JUSQU_A=40.,
NOMBRE=70,),
);
DYNAHARM = DYNA_LINE_HARM (
MATR_MASS = MATASSM,
MATR_RIGI = MATASSR,
MATR_AMOR = MATASSC,
LIST_FREQ = LIFREQ,
EXCIT=_F
(VECT_ASSE = VECTASS,
COEF_MULT = 1. ,) ,
);
IMPR_RESU (MODELE = MO, RESU =_F
( RESULTAT = DYNAHARM,),);
FIN ();