Code - Aster: Modélisation de L'amortissement en Dynamique Linéaire

Version 10
Code_Aster
Titre : Modélisation de l'amortissement en dynamique linéa[...] Date : 04/11/2011 Page : 1/13
Responsable : Emmanuel BOYERE Clé : R5.05.04 Révision : 7679
Modélisation de l'amortissement en dynamique

linéaire
Résumé :
Les analyses dynamiques linéaires des structures soumises à des forces ou des mouvements imposés
nécessitent d'ajouter des caractéristiques d'amortissement mécanique aux caractéristiques de rigidité et de
masse du modèle.
On dispose de plusieurs modélisations classiques, applicables à tous les types d'éléments finis disponibles :
• le modèle d'amortissement visqueux,

• le modèle d'amortissement hystérétique (dit aussi "amortissement structural")
pour l'analyse harmonique des matériaux viscoélastiques.
Pour les analyses utilisant les méthodes de réponse dynamique par recombinaison modale, avec une base
modale de modes propres réels, il est possible d'introduire des coefficients d'amortissement modaux.
Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Version 10
Code_Aster
Table des matières

1 Notion d'amortissement mécanique...................................................................................... 3
1.1 Modèles d'amortissement............................................................................................... 3
1.2 Définitions générales pour caractériser l'amortissement [bib1]...................................... 3
1.2.1 Coefficient de perte............................................................................................... 3
1.2.2 Amortissement réduit............................................................................................ 3
2 Modèle d'amortissement visqueux....................................................................................... 4
2.1 Définition physique de l'amortissement visqueux.......................................................... 4
2.2 Oscillateur harmonique avec amortissement visqueux.................................................. 4
2.2.1 Réponse à un lâcher d'excitation.......................................................................... 5
2.2.2 Réponse à une excitation harmonique.................................................................. 5
3 Modèle d'amortissement hystérétique.................................................................................. 5
3.1 Définition physique de l'amortissement hystérétique..................................................... 5
3.2 Oscillateur harmonique avec amortissement hystérétique............................................ 6
4 Autres modèles d'amortissement.......................................................................................... 7
5 Analyse de structure avec amortissement............................................................................ 7
5.1 Amortissement global de la structure............................................................................. 8
5.1.1 Amortissement visqueux proportionnel "global".................................................... 8
5.1.2 Influence des coefficients d'amortissement proportionnel.................................... 8
5.1.3 Amortissement hystérétique "global"..................................................................... 9
5.2 Amortissement localisé.................................................................................................. 10
5.2.1 Eléments amortisseurs.......................................................................................... 10
5.2.2 Amortissement affecté à tout type d'élément fini.................................................. 10
5.2.3 Construction de la matrice d'amortissement......................................................... 11
6 Utilisation de la matrice d'amortissement............................................................................. 11
6.1 Utilisation de la matrice d'amortissement visqueux....................................................... 11
6.1.1 Analyse dynamique linéaire directe...................................................................... 11
6.1.2 Analyse dynamique par recombinaison modale.................................................... 12
6.2 Utilisation de la matrice de rigidité complexe ................................................................ 12
6.3 Analyse modale complexe............................................................................................. 13
7 Bibliographie......................................................................................................................... 13
8 Description des versions du document................................................................................. 13

Version 10
Code_Aster
1 Notion d'amortissement mécanique

1.1 Modèles d'amortissement
Le mouvement des structures soumises à des forces ou des mouvements imposés, variables au cours
du temps, dépend, en particulier des propriétés d'amortissement, c'est-à-dire de la dissipation
d'énergie dans les matériaux constitutifs de la structure et dans les liaisons des différents éléments de
structure entre eux et avec le milieu environnant.
Les phénomènes physiques intervenant dans cette dissipation d'énergie sont nombreux frottements,
interaction fluide-structure dans une lame fluide, chocs, viscosité et plasticité, rayonnement vibratoire
aux appuis.
Les modèles de comportement représentant ces phénomènes sont souvent mal connus et il n'est pas
possible de les décrire explicitement au niveau élémentaire. C'est pourquoi les modèles les plus
utilisés sont les modèles simples qui permettent de reproduire à l'échelle macroscopique les
principaux effets sur les structures [bib1] [bib2]. Ceux disponibles actuellement dans Code_Aster
sont :
• l'amortissement visqueux : énergie dissipée proportionnelle à la vitesse du mouvement,

• l'amortissement hystérétique (dit aussi "amortissement structural") : énergie dissipée
proportionnelle au déplacement telle que la force d'amortissement est de signe opposée à
celui de la vitesse.
Notons que l'amortissement de Coulomb, qui correspond à un amortissement de frottement pour

lequel l'énergie dissipée est proportionnelle à la force de réaction normale à la direction de
déplacement n'est pas implanté actuellement dans Code_Aster.
Les valeurs des paramètres de ces modèles sont déduits de résultats expérimentaux. Au stade de la
conception, on se limite à l'utilisation de valeurs réglementaires.
1.2 Définitions générales pour caractériser l'amortissement [bib1]

1.2.1 Coefficient de perte
Le coefficient de perte  est un coefficient adimensionnel caractéristique de l'effet amortisseur défini

comme le rapport de l'énergie dissipée au cours d'un cycle à l'énergie potentielle maximum multipliée
par 2 :
E d par cycle
= éq 1.2-1
2  E p max
1.2.2 Amortissement réduit

Par définition l'amortissement réduit est égal à la moitié du coefficient de perte

= éq 1.2-2
2

Version 10
Code_Aster
2 Modèle d'amortissement visqueux

2.1 Définition physique de l'amortissement visqueux
Les dispositifs amortisseurs classiques (laminage d'un fluide visqueux à travers les orifices d'un piston
entraîné par le mouvement vibratoire) délivrent des forces proportionnelles à la vitesse du
mouvement et de signe opposé. Au cours d'un cycle, le travail de ces forces est positif : c'est
l'amortissement visqueux.
k c
m
f u
Pour un oscillateur simple de rigidité k , de masse m et d'amortissement visqueux c , la force
extérieure appliquée équilibre les trois composantes : force de rappel élastique ku , force
d'amortissement c u̇ et force d'inertie m ü d'où l'équation dynamique en mouvement absolu :
m üc u̇ku= f éq 2.1-1
Pour ce modèle d'amortissement visqueux l'énergie dissipée au cours d'un cycle de pulsation  est
proportionnelle à la vitesse vibratoire −u 0 sin  t  associée au déplacement u 0 cos   t  :
2
Ed par cycle=∫ 0
−c u0 sin td  u0 cos  t  =c u20
et l'énergie potentielle pour un déplacement sinusoïdal u0 cos t est :
0 1
E p max =∫  /2 ku 0 cos t d  u0 cos  t  = ku20
2
Pour un cycle de pulsation  et de déplacement sinusoïdal u 0 cos  t , le coefficient de perte est
proportionnel à la fréquence du mouvement :
c
= éq 2.1-2
k
2.2 Oscillateur harmonique avec amortissement visqueux

L'analyse classique du modèle non amorti associé à l'équation [éq. 2.1-1], mise sous la forme
 k−m   u=0
2
nous donne 0=
 k la pulsation propre.
m
L'amortissement critique à partir duquel l'équation différentielle [éq 2.1-1] n'a plus de solution
2k
oscillante est donné par les formules c critique =2  km=2m 0= ce qui permet de donner une
0
interprétation numérique de l'amortissement réduit, qui est souvent exprimé en pourcentage de
l'amortissement critique :
 c c
= = = éq 2.1-3
2 c critique 2m 0

Version 10
Code_Aster
2.2.1 Réponse à un lâcher d'excitation

f0
A partir d'une déformation statique u st = , un lâcher (libération du système) produit un
k
-  t
mouvement oscillatoire libre u  t =u e
l 0
cos  '0 t qui fait apparaître la pulsation propre du
0
système amorti ' 0=0   1−  .

2
Au cours du temps, l'amplitude extrémale u 1

, u2  diminue à chaque période de
− 0 T
e =e−2   =e− où  est le décrément logarithmique : =2 
2.2.2 Réponse à une excitation harmonique

j t
La réponse à une excitation harmonique de la forme f  t = f 0 e s’écrit avec une réponse forcée
 j  t− qui s’écrit avec la pulsation réduite λ= 

solution particulière permanente u t =u 0 e
0
ku 0 1
= 2
= H v  j   où H v  j  est la fonction de transfert complexe d’un oscillateur
f0 1− λ  j2  λ
simple avec amortissement visqueux.
u 0 ku 0 1
Le module de la réponse = =∣H v  j  ∣= fait apparaître une
u st f 0
  1−   2  %lambda 
2 2 2
amplification dynamique par rapport à la réponse statique u st .

'0
Cette amplification est maximale pour λ=
0

=  1−2  et donne la valeur du déplacement
u 0 max 1
maximal = . Si l'on observe la vitesse vibratoire u̇ t = j  u  t  , l'amplification
u st 2   1− 2
0
de la vitesse vibratoire est maximale pour λ= =1 et l'amplitude maximale de la vitesse est
0
1
u̇ 0 max = =Q , où Q est l'analogie mécanique du facteur de surtension des électriciens qui nous
2
font vivre. Ces propriétés sont à l'origine des méthodes de mesure des caractéristiques
d'amortissement des structures mécaniques.

Version 10
Code_Aster
3 Modèle d'amortissement hystérétique

3.1 Définition physique de l'amortissement hystérétique
Pour une excitation sinusoïdale appliquée à une structure élasto-plastique ou à une structure élastique
avec frottement, la courbe force-déplacement fait apparaître un travail positif de la force extérieure
qui correspond à une énergie dissipée dans la structure, que l'on peut en première approximation
représenter comme ci-dessous :
f f
u u
Elastoplastique Glissement
Dans les deux cas le coefficient de perte croit, en général avec l'amplitude du cycle. Pour des valeurs
faibles du coefficient de perte (< 0.2), la forme du cycle n'a pas d'influence sensible sur le mouvement
et on peut l'assimiler à une ellipse [bib1].
Dans le cas particulier d'une relation force-déplacement dont le cycle est de forme elliptique,
l'expression du coefficient de perte est simple. Pour une force appliquée F et un déplacement
u=u 0 cos  la force de rappel est ku 0 cos q et la force d'amortissement −hu0 sin  ce qui conduit
à la relation d'équilibre F =ku0 cos −hu 0 sin  .
f
k u0 θ=0
θ = −π → h u0
2
u0 u
k θ = π → -hu0
1 2
Les énergies dissipée au cours d'un cycle et potentielle maximum sont
2 0 ku02
 
E d par cycle=∫0 −hu0 sin  d u 0 cos  = hu et E pm ax =∫/2 ku 0 cos  d u 0 cos  =
2
0   2
d'où le coefficient de perte
 hu 20 h
= 2
=
ku 0
k éq 3.1-1
2
2

Version 10
Code_Aster
h
Pour un cycle sinusoïdal = t , Le coefficient d'amortissement hystérétique = est
k
indépendant de  . Il peut être déterminé à partir d'un essai sous chargement cyclique harmonique.
3.2 Oscillateur harmonique avec amortissement hystérétique

Le modèle d'amortissement hystérétique est utilisable pour traiter les réponses harmoniques de
structures avec des matériaux viscoélastiques.
2
L'énergie dissipée par cycle sous la forme E d par cycle=∫0  d  permet de mettre en évidence un
module d'YOUNG complexe E * à partir de la relation contrainte-déformation d'un matériau
j t j   t− 
viscoélastique = 0 e et =0 e où  0 et  0 sont les amplitudes et  la phase :
 
 

E *= = 0 e j = 0  cos  j sin 
 0 0
En notant E 1=
 0
0
cos  la partie réelle et E 2=
0
0
sin  la partie imaginaire on obtient
E1
E =E 1 jE 2= E 1  1 j  avec =
*
=tg  , où j est aussi appelé angle de perte.
E2
L'analyse classique de l'équation [éq 2.1-1] n'a de sens, avec un modèle d'amortissement
j t
hystérétique, que pour une excitation harmonique f  t = f 0 e qui conduit à l'équation
m ük  1 j   u=m ü  k  jh  u= f 0 e j  t éq 3.2-1
où la partie réelle du déplacement u représente le déplacement de la masse et h=k  .

Comme précédemment Cf. [§ 2.2], la réponse harmonique peut s'écrire, avec la pulsation réduite
 ku 0 1
λ= , sous la forme = = H h  j   où H h  j   est la fonction de transfert
0 f0 1− λ 2 j 
complexe d’un oscillateur simple avec amortissement hystérétique.
u0 ku0 1
Le module de la réponse = =∣H h  j ∣= fait apparaître une amplification
u st f0
  1− λ  
2
2
2
dynamique par rapport à la réponse statique, amplification qui est maximale pour λ=1 et donne la
u 0 max 1 1
valeur du déplacement maximal = = .
u st  2
En conclusion, l'amortissement réduit associé à l'amortissement hystérétique est :
 h h
= = = éq 3.2-2
2 2k 2m 02

Version 10
Code_Aster
4 Autres modèles d'amortissement

On ne traite pas ici des modèles représentant l'amortissement "ajouté" par les fluides immobiles
confinés ou les fluides en mouvements. On se reportera au fascicule [R4.07] couplage
fluide-structure.
5 Analyse de structure avec amortissement

Les modélisations présentées ne sont pas aisément généralisables aux différentes analyses de
structures Cf. [§1].
Remarque :
Les deux modélisations n'ont pas le même domaine d'analyse linéaire :
• l'amortissement visqueux est utilisable en analyse transitoire ou harmonique,
• l'amortissement hystérétique n'est utilisable qu'en analyse harmonique.
Les options de modélisations dans Code_Aster permettent la définition :
• d'un amortissement global pour la structure,
• d'amortissements localisés sur des mailles ou des groupes de mailles.
5.1 Amortissement global de la structure

En l'absence d'informations suffisantes sur les composants et liaisons créant une dissipation
d'énergie, une modélisation courante consiste à construire une matrice d'amortissement "global".
5.1.1 Amortissement visqueux proportionnel "global"

On se place dans le cadre des équations classiques de la dynamique des structures linéaires :
M Ü C U̇  KU =F  t  éq 5.1.1-1
La notion d'amortissement de RAYLEIGH permet de définir la matrice d'amortissement C comme

combinaison linéaire des matrices de rigidité et de masse :
C= K  M éq 5.1.1-2
Avantages :
• facile à mettre en œuvre en utilisant l'opérateur COMB_MATR_ASSE [U4.53.01], après avoir
assemblé les matrices de rigidité et de masse à coefficients réels.;
• utile pour la validation d'algorithmes de résolution ;
• historiquement, son succès est attaché aux méthodes d'analyse transitoire par recombinaison
modale à partir d'une base de modes propres réels.
Les propriétés d'orthogonalité des modes propres réels solution du problème aux valeurs
propres  K −M  =0 se traduisent par la diagonalisation simultanée dans le passage
2
en coordonnées modales généralisées de T K  et T M  .

L'amortissement de RAYLEIGH est une condition suffisante pour diagonaliser T C  .
T C  T
Le système d'équations modales q̈ T q̇ 2 q= T F  t  devient alors
 M  M
diagonal.
T
q̈2   q̇2 q= F t éq 5.1.1-3
T M 

Version 10
Code_Aster
Inconvénients :
• Cette modélisation ne permet pas de représenter l'hétérogénéité de la structure par rapport à
l'amortissement.
• L'amortissement effectivement introduit dans le modèle dépend fortement de l'identification
des coefficients  et  Cf. [§ 5.1.2].
5.1.2 Influence des coefficients d'amortissement proportionnel

Trois cas d'identification simples sont présentés ici pour illustrer, les effets induits par cette
modélisation :
• amortissement proportionnel aux caractéristiques d'inertie : α=0, β= β i
Ce cas a été très utilisé en résolution transitoire directe : si la matrice de masse est
diagonale, celle d'amortissement l'est encore et le gain en place mémoire est évident.
Le coefficient  peut être identifié à l'amortissement réduit expérimental x i du mode
propre  i , i  qui participe le plus à la réponse Cf. [éq. 2.1-1] d'où β i =2 i i . Pour
i
toute autre pulsation on obtient un amortissement modal réduit = β i . Les modes

élevés  >> i seront très peu amortis et les modes basse fréquence i trop
amortis.
ξ
2 ξ ω =β Μ
ξi
ωi ω
• amortissement proportionnel aux caractéristiques de rigidité : α=α j , β=0 .
 peut être identifié, comme précédemment à partir de x j associé au
Le coefficient
 
mode  j ,  j d'où α j =2 j  j . Pour toute autre pulsation on obtient un amortis-

sement modal réduit =α j . Les modes élevés  >>  j sont très amortis, etc....
j
ξ
2ξ ω= α Κ
ξj
ωj ω
• amortissement proportionnel complet : α=α j , β= β i

A partir d'une identification sur deux modes indépendants  j , w  et  j
i i j 
, w j , nous
obtiendrons pour toute autre pulsation un amortissement modal réduit
w wi
x=α j βi [ ]
. Dans l'intervalle w i , w j , la variation de l'amortissement réduit est
wj w
faible et en dehors on retrouve la combinaison des inconvénients précédents : les modes
extérieurs à l'intervalle sont trop amortis.

Version 10
Code_Aster
Dans aucun des cas précédents, on ne pourra reproduire une hypothèse d'amortissement modal égal
pour tous les modes. Des méthodes ont été imaginées pour tendre vers cet objectif [bib1].
5.1.3 Amortissement hystérétique "global"

La généralisation de l'équation de l'oscillateur simple avec amortissement hystérétique conduit au
 
système d'équations complexes ou F  est une excitation harmonique.
M Ü  K  1 j   U =F    Éq 5.1.3-1
Connaissant la matrice de rigidité réelle, il est possible de construire une matrice d'amortissement
hystérétique K h= j  K , avec un cœfficient de perte "global"  .
Comme précédemment en résolution par recombinaison modale, à partir d'une base de modes
propres réels, on obtient T M  q̈ j  T K h  qT K  q=T F  t  où la matrice
d'amortissement hystérétique généralisée est diagonale [ ]

T K h  = diag   i , comme la matrice
de rigidité généralisée [
T K = diag  i . ]
D'après la définition de l'amortissement réduit (Cf. [éq 1.2-2]), l'amortissement modal est constant

pour tous les modes d'où =
2
Avantages :
• facile à mettre en œuvre en utilisant l'opérateur COMB_MATR_ASSE [U4.53.01], après avoir
assemblé la matrices de rigidité.;
• très utile pour la validation d'algorithmes de résolution ;
• l'amortissement effectivement introduit dans le modèle est constant pour tous les modes de
la structure, comme le demande des règlements de construction.
Inconvénients :
• cette modélisation est mal adaptée pour les études industrielles, car elle ne permet pas de
représenter l'hétérogénéité de la structure par rapport à l'amortissement.
• seule l'analyse harmonique (en complexe) est possible.
•
5.2 Amortissement localisé
Pour les analyses nécessitant une modélisation représentant l'hétérogénéité de la structure, il est
possible d'affecter des caractéristiques d'amortissement localisées sur les mailles de la structure, en
fait sur des éléments du modèle.
5.2.1 Eléments amortisseurs

Il est possible d'appliquer des éléments amortisseurs discrets :

Version 10
Code_Aster
• sur des mailles POI1 : l'amortissement est lié au déplacement (respectivement la vitesse) du
nœud support,
• sur des mailles SEG2 : l'amortissement est lié au déplacement relatif (respectivement la
vitesse relative) des deux nœuds reliés.
L'opérateur AFFE_CARA_ELEM [U4.24.01] permet de définir pour chaque élément discret :
• une matrice d'amortissement de type visqueux a discret dont les termes sont affectés aux
différents degrés de liberté des nœuds concernés ; plusieurs modes de description de la
matrice sont disponibles.
• un coefficient de perte hystérétique discret multiplicateur de la matrice de rigidité de
l'élément discret affecté à la maille support.
5.2.2 Amortissement affecté à tout type d'élément fini

Le matériau élastique affecté à tout élément fini peut être défini avec des paramètres d'amortissement
par l'opérateur DEFI_MATERIAU [U4.23.01] :
• Amortissement visqueux proportionnel avec deux paramètres de RAYLEIGH  et  .
AMOR_ALPHA : 
AMOR_BETA : 
Pour tous les types d'éléments finis (de milieux continus, structuraux ou discrets), il est
possible de calculer les matrices élémentaires réelles correspondant à l'option de calcul
'AMOR_MECA', après avoir calculé les matrices élémentaires correspondant aux options de
calcul 'RIGI_MECA' et 'MASS_MECA'.
La matrice élémentaire de l'élément i affecté du matériau  j ,  j est alors de la forme :
• pour un élément fini
c elem i= j k elem i  j m elem i
• pour un élément discret
c elem i =a discret i
• Amortissement hystérétique avec un coefficient de e
AMOR_HYST : coef
Pour tous les types d'éléments finis (de milieux continus, structuraux ou discrets), il est
possible de calculer les matrices élémentaires complexes correspondant à l'option de calcul
'RIGI_MECA_HYST', après avoir calculé les matrices élémentaires correspondant aux
options de calcul 'RIGI_MECA'.
La matrice élémentaire de l'élément i affecté du matériau  j ,  j est alors de la forme :

• pour un élément fini

k *elem i =k elem i 1 j  j 
• pour un élément discret

k *elem i =k elem i 1 j  discret i 
5.2.3 Construction de la matrice d'amortissement
L'assemblage des matrices élémentaires d'amortissement est obtenu avec l'opérateur
ASSE_MATRICE habituel [U4.42.02] ou par la macro commande MACRO_MATR_ASSE [U4.31.02]. On
doit utiliser les mêmes numérotations et le même mode de stockage que pour les matrices de rigidité
et de masse (opérateur NUME_DDL [U4.42.01]).

Version 10
Code_Aster
Remarque :
La matrice d'amortissement obtenue est non proportionnelle
C≠ K  M ou K h ≠ j  K
6 Utilisation de la matrice d'amortissement

6.1 Utilisation de la matrice d'amortissement visqueux
6.1.1 Analyse dynamique linéaire directe
La matrice d'amortissement visqueux C , quel que soit son mode d'élaboration et son caractère
proportionnel ou non proportionnel, est utilisable pour l'analyse dynamique linéaire directe (mot clé
MATR_AMOR) avec les opérateurs :
• d'analyse transitoire DYNA_LINE_TRAN [R5.05.02] et [U4.54.01]
• d'analyse harmonique DYNA_LINE_HARM [R5.05.03] et [U4.54.02]
6.1.2 Analyse dynamique par recombinaison modale

Pour les analyses par recombinaison modale, on doit projeter cette matrice dans le sous-espace défini
par un ensemble  de modes propres réels, obtenus sur le problème non amorti associé
 K − 2
M  =0 .
Cette opération est possible avec la macro commande MACRO_PROJ_BASE [U4.55.11] ou avec
l'opérateur PROJ_MATR_BASE [U4.55.01].
Pour le calcul de la réponse dynamique en force ou en mouvement imposé dans l'espace modal, on
dispose de possibilités suivantes :
• utilisation de la matrice d'amortissement généralisée T C  :
• en analyse transitoire avec l'opérateur DYNA_TRAN_MODAL [R5.06.04] et [U4.54.03] et le
mot clé AMOR_GENE,
• en analyse sismique par méthode spectrale avec l'opérateur COMB_SISM_MODAL
[R4.05.03] et [U4.54.04] et le mot clé AMOR_GENE,
• en analyse harmonique avec l'opérateur DYNA_LINE_HARM [R5.05.03] et [U4.54.02] et
le mot clé MATR_AMOR.
Rappelons que dans le cas de l'amortissement hétérogène (utilisation des options
d'amortissement localisé), la matrice T C  n'est par diagonale.
• utilisation de l'amortissement modal visqueux en fournissant un amortissement modal réduit
constant pour tous les modes  ou une liste de valeurs  i .
Plusieurs méthodes d'identification de ces coefficients sont possibles mais il n'existe pas de
commande de construction automatique de la liste de valeurs : utilisation de l'hypothèse de
BASILE
 2 i i =diag
T C 
T M   , méthode du règlement RCC-G pour l'analyse sismique
avec amortissement du sol, exploitation de résultats expérimentaux, ...

• en analyse transitoire avec l'opérateur DYNA_TRAN_MODAL [R5.06.04] [U4.54.03] et le
mot clé AMOR_REDUIT.
• en analyse sismique par méthode spectrale avec l'opérateur COMB_SISM_MODAL
[R4.05.03] [U4.54.04] et les mots clé AMOR ou LIST_AMOR. Une évolution est demandée
pour généraliser le mot clé AMOR_REDUIT.

Version 10
Code_Aster
• en analyse harmonique une demande d'évolution avec l'opérateur DYNA_LINE_HARM

[R5.05.03] [U4.54.02] est déposée. Elle n'est pas traitée en version 3.6.
Pour les analyses par sous-structuration dynamique, avec l'utilisation d'une base modale (base de
RITZ) on se reportera à [R4.06.03] et [R4.06.04].
6.2 Utilisation de la matrice de rigidité complexe

La matrice de rigidité complexe K *= K  K h , où K h est une matrice imaginaire (au sens des
complexes !), est utilisable pour l'analyse harmonique directe avec l'opérateur DYNA_LINE_HARM
[R5.05.03] et [U4.54.02] et le mot clé MATR_RIGI.
Pour les analyses par recombinaison modale, aucune fonctionnalité n'est disponible actuellement pour
l'utilisation du modèle d'amortissement hystérétique.
6.3 Analyse modale complexe

La matrice d'amortissement visqueux C est indispensable pour l'analyse modale complexe avec les
opérateurs traitant le problème quadratique aux valeurs propres [R5.01.02]:
• par itérations inverses MODE_ITER_INV [U4.52.01]

• par itérations simultanées MODE_ITER_SIMULT [U4.52.02]
Rappelons que les modes propres complexes permettent une approche mieux adaptée à l'étude
dynamique des structures fortement amorties (amortissement réduit  > 20%). A ce jour aucun outil
de réponse dynamique par recombinaison modale utilisant une base de modes propres complexes
n'est disponible dans Code_Aster.
7 Bibliographie
1) "Amortissement" : Françoise GANTENBEIN et Michel LIVOLANT (CEA-DMT) in "Génie
Parasismique" Ouvrage collectif - Presses de l'E.N.P.C. (1985)
2) "Amortissement dans les calculs de structures" : Forum IPSI -  2 AS Volume XVIII N°2 (Juin
1994)
8 Description des versions du document

Version Auteur(s) Description des modifications
Aster Organisme(s)
3 J.R. LEVESQUE Texte initial
(EDF-R&D/IMA/MMN)

Code - Aster: Modélisation de L'amortissement en Dynamique Linéaire

Transféré par

Droits d'auteur :

Formats disponibles

Code - Aster: Modélisation de L'amortissement en Dynamique Linéaire

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Code - Aster: Modélisation de L'amortissement en Dynamique Linéaire

Transféré par

Droits d'auteur :

Formats disponibles

Version 10

Modélisation de l'amortissement en dynamique

• le modèle d'amortissement visqueux,

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Table des matières

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

1 Notion d'amortissement mécanique

• l'amortissement visqueux : énergie dissipée proportionnelle à la vitesse du mouvement,

Notons que l'amortissement de Coulomb, qui correspond à un amortissement de frottement pour

1.2 Définitions générales pour caractériser l'amortissement [bib1]

Le coefficient de perte  est un coefficient adimensionnel caractéristique de l'effet amortisseur défini

1.2.2 Amortissement réduit

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

2 Modèle d'amortissement visqueux

m üc u̇ku= f éq 2.1-1

2.2 Oscillateur harmonique avec amortissement visqueux

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

2.2.1 Réponse à un lâcher d'excitation

système amorti ' 0=0   1−  .

Au cours du temps, l'amplitude extrémale u 1

2.2.2 Réponse à une excitation harmonique

 j  t− qui s’écrit avec la pulsation réduite λ= 

amplification dynamique par rapport à la réponse statique u st .

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

3 Modèle d'amortissement hystérétique

Les énergies dissipée au cours d'un cycle et potentielle maximum sont

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

3.2 Oscillateur harmonique avec amortissement hystérétique

m ük  1 j   u=m ü  k  jh  u= f 0 e j  t éq 3.2-1

où la partie réelle du déplacement u représente le déplacement de la masse et h=k  .

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

4 Autres modèles d'amortissement

5 Analyse de structure avec amortissement

5.1 Amortissement global de la structure

5.1.1 Amortissement visqueux proportionnel "global"

La notion d'amortissement de RAYLEIGH permet de définir la matrice d'amortissement C comme

en coordonnées modales généralisées de T K  et T M  .

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

5.1.2 Influence des coefficients d'amortissement proportionnel

• amortissement proportionnel complet : α=α j , β= β i

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

5.1.3 Amortissement hystérétique "global"

d'amortissement hystérétique généralisée est diagonale [ ]

5.2.1 Eléments amortisseurs

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

5.2.2 Amortissement affecté à tout type d'élément fini

La matrice élémentaire de l'élément i affecté du matériau  j ,  j est alors de la forme :

Manuel de référence Fascicule r5.05 : Dynamique transitoire ou harmonique

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

6 Utilisation de la matrice d'amortissement