Exemples
========

.. code-block:: text




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',);




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),

);




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'),

);