4. Exemples#
4.1. Définition de la dépendance en fréquence des propriétés des matériaux viscoélastiques#
# fréquences pour lesquelles les paramètres du matériaux sont donnés
list_f=DEFI_LIST_REEL(VALE=(1,10,50,100,500,1000,1500,),);
# valeurs (de la partie réelle) du module de Young aux fréquences de list_f
list_E=DEFI_LIST_REEL(VALE=(23.2E6,58.E6,145.E6,203.E6,348.E6,435.E6,464.E6,),);
# valeurs du facteur de perte aux fréquences de list_f
list_eta=DEFI_LIST_REEL(VALE=(1.1,0.85,0.7,0.6,0.4,0.35,0.34,),);
fonc_E=DEFI_FONCTION(NOM_PARA=”FREQ”,
VALE_PARA=list_f,
VALE_FONC=list_E,
INTERPOL=(“LIN”,”LIN”,),
PROL_DROITE=”LINEAIRE”,
PROL_GAUCHE=”CONSTANT”,);
fonc_eta=DEFI_FONCTION(NOM_PARA=”FREQ”,
VALE_PARA=list_f,
VALE_FONC=list_eta,
INTERPOL=(“LIN”,”LIN”,),
PROL_DROITE=”LINEAIRE”,
PROL_GAUCHE=”CONSTANT”,);
4.2. Calcul des modes propres complexes#
modes=DYNA_VISCO(MODELE=modele,
CARA_ELEM=cara_ele,
# matériaux à propriétés constantes :
MATER_ELAS=_F(E=2.1e11,
NU=0.3,
RHO=7800.,
AMOR_HYST=0.002,
GROUP_MA=”SUPPORT”),
# matériaux à propriétés dépendantes de la fréquence :
MATER_ELAS_FO =_F(E=fonc_E,
AMOR_HYST=fonc_eta,
RHO=1200.,
NU=0.45,
GROUP_MA=”VISCO”),
TYPE_RESU=”MODE”,
TYPE_MODE=”BETA”,
# bande fréquentielle de recherche
FREQ=(1.,1500.),
EXCIT=_F(CHARGE=condlim),
);
4.3. Calcul de la réponse harmonique#
# DEFINITION DU CHARGEMENT
excit=AFFE_CHAR_MECA(MODELE=modele,
FORCE_NODALE=_F(GROUP_NO=”A”,
FZ=1.,),);
# DEFINITION DES FREQUENCES DE CALCUL DE LA REPONSE
listfr=DEFI_LIST_REEL(DEBUT=1.,
INTERVALLE=(_F(JUSQU_A=500.,
PAS=1.,),),);
# REPONSE HARMONIQUE
visco=DYNA_VISCO(MODELE=modele,
CARA_ELEM=cara_ele,
EXCIT=_F(CHARGE=(condlim, excit),),
MATER_ELAS=(_F(E=2.1E11,
NU=0.3,
RHO=7800.,
AMOR_HYST=0.002,
GROUP_MA=”DESSOUS”),
_F(E=7.0E10,
NU=0.3,
RHO=2700.,
AMOR_HYST=0.001,
GROUP_MA=”DESSUS”),),
MATER_ELAS_FO=(_F(E=fonc_E,
AMOR_HYST=fonc_eta,
RHO=1200.,
NU=0.45,
GROUP_MA=”VOLUME”,),),
TYPE_RESU=”HARM”,
TYPE_MODE=”REEL”,
LIST_FREQ=listfr,
# champs à sauvegardés
NOM_CHAM=(“DEPL”,”VITE”),
# sauvegarde des modes propres de la structure :
MODE_MECA=CO(“modes”),
);