Code_Aster ®
Version
6.4
Titre :
Notice de modélisation de l'amortissement mécanique
Date :
08/11/02
Auteur(s) :
J. PELLET Clé
:
U2.06.03-B Page
: 1/6
Organisme(s) : EDF-R&D/AMA
Manuel d'Utilisation
Fascicule U2.06 : Dynamique
Document : U2.06.03
Notice de modélisation de l'amortissement
mécanique
Résumé
Les analyses dynamiques linéaires et non-linéaires, pour l'étude de la réponse vibratoire avec une excitation en
force ou en mouvement imposé ou pour l'analyse modale complexe, nécessitent d'ajouter des caractéristiques
d'amortissement mécanique aux caractéristiques de rigidité et de masse.
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 une base modale de modes propres réels, il est possible d'introduire des coefficients
d'amortissement modaux.
Manuel d'Utilisation
Fascicule U2.06 : Dynamique
HT-66/02/003/A
Code_Aster ®
Version
6.4
Titre :
Notice de modélisation de l'amortissement mécanique
Date :
08/11/02
Auteur(s) :
J. PELLET Clé
:
U2.06.03-B Page
: 2/6
1
Modèle d'amortissement visqueux
Le modèle d'amortissement visqueux est le plus couramment utilisé. Il correspond à la modélisation
d'une énergie dissipée proportionnelle à la vitesse vibratoire :
1
1
Ed = vT Cv = uT Cu
éq 1-1
2
2
où C est la matrice d'amortissement visqueux, à coefficients réels.
Il conduit aux équations classiques de la dynamique des structures :
Mu + Cu + Ku = f (t)
éq 1-2
avec K matrice de rigidité et M matrice de masse.
1.1
Amortissement visqueux proportionnel "global"
Cette modélisation, facile à mettre en oeuvre, correspond à :
C = K + M
éq 1.1-1
Elle est disponible actuellement, en utilisant l'opérateur COMB_MATR_ASSE [U4.72.01], après avoir
assemblé les matrices de rigidité et de masse à coefficients réels, mais elle est d'une utilité faible :
· validation d'algorithmes de résolution,
· inutile 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 (dissipation aux appuis ou aux assemblages). De
plus l'identification globale des coefficients et n'est possible, en analyse modale
expérimentale, que pour deux fréquences propres [ f1 f2 ] distinctes; elle donne, pour les
fréquences propres fi [ f , f
i
2 ] avec i = 2 f
i , une loi d'évolution de l'amortissement
réduit de la forme :
i = i +
i
1.2
Amortissement visqueux proportionnel des éléments du modèle
1.2.1 Caractéristiques
d'amortissement
Il est possible de construire une matrice d'amortissement à partir de chaque élément du modèle,
comme pour la rigidité et la masse.
Deux fonctionnalités sont utilisables :
· l'affectation d'éléments discrets, sur des mailles POI1 ou SEG2, par l'opérateur AFFE_CARA_ELEM
[U4.42.01]. Celui-ci permet de définir, avec plusieurs modes de description possibles, une matrice
d'amortissement pour chaque degré de liberté.
· la définition d'une caractéristique d'amortissement pour tout matériau élastique par l'opérateur
DEFI_MATERIAU [U4.43.01] par :
AMOR_ALPHA
:
[R]
AMOR_BETA
:
ce matériau étant ensuite affecté aux mailles concernées.
Manuel d'Utilisation
Fascicule U2.06 : Dynamique
HT-66/02/003/A
Code_Aster ®
Version
6.4
Titre :
Notice de modélisation de l'amortissement mécanique
Date :
08/11/02
Auteur(s) :
J. PELLET Clé
:
U2.06.03-B Page
: 3/6
1.2.2 Calcul des matrices d'amortissement
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' ou 'MASS_MECA_DIAG'. Chaque matrice élémentaire est alors de la forme :
· quand le matériau i, de caractéristiques d'amortissement visqueux proportionnel (i i), est
affecté à l'élément elem
celem = i kelem + i melem
· pour un élément discret
celem = adiscret
Cette opération est possible avec :
mel
[matr_elem_DEPL_R] = CALC_MATR_ELEM
(
/
OPTION:
'AMOR_MECA'
MODELE:
mo
[modele]
CHAM_MATER:
chmat
[cham_mater]
CARA_ELEM:
cara
[cara_elem]
);
L'assemblage de toutes les matrices élémentaires d'amortissement est obtenu avec l'opérateur
ASSE_MATRICE habituel [U4.61.22]. On notera que l'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.61.11]).
On remarque que la matrice d'amortissement obtenue est, en général, non proportionnelle :
C ° = K + M
1.2.3 Utilisation de la matrice d'amortissement visqueux
La matrice C est utilisable pour l'analyse dynamique linéaire directe (mot-clé MATR_AMOR) avec les
opérateurs de réponse dynamique linéaire :
· transitoire
DYNA_LINE_TRAN
[U4.53.02]
· harmonique
DYNA_LINE_HARM
[U4.53.11]
Elle est indispensable pour l'analyse modale complexe avec les opérateurs de recherche des valeurs
propres :
· par itérations inverses
MODE_ITER_INV
[U4.52.04]
· par itérations simultanées
MODE_ITER_SIMULT
[U4.52.03]
Pour les analyses en base modale, on doit projeter cette matrice dans le sous-espace défini par un
ensemble de modes propres réels. Cette opération est possible avec l'opérateur
PROJ_MATR_BASE [U4.63.12]. Notons que dans le cas général ( C non proportionnelle), la matrice
projetée n'est pas diagonale. Elle reste néanmoins utilisable (mot-clé AMOR_GENE) pour le calcul de la
réponse dynamique en force ou en mouvement imposé dans l'espace modal, avec l'opérateur de
réponse dynamique linéaire :
· transitoire
DYNA_TRAN_MODAL
[U4.53.21]
Manuel d'Utilisation
Fascicule U2.06 : Dynamique
HT-66/02/003/A
Code_Aster ®
Version
6.4
Titre :
Notice de modélisation de l'amortissement mécanique
Date :
08/11/02
Auteur(s) :
J. PELLET Clé
:
U2.06.03-B Page
: 4/6
1.2.4 Utilisation de l'amortissement modal visqueux
Pour les analyses en base modale de modes propres réels, l'équation différentielle dynamique en
coordonnées généralisées :
T
2
q& + 2 q& + q
i
=
f
i
i
i
i
i
i
(t)
éq
1.2.4-1
µi
fait apparaître un coefficient d'amortissement modal i exprimé comme une fraction de
l'amortissement critique et la masse généralisée du mode µi , qui dépend du mode de normalisation
du mode propre.
Dans le cas d'une matrice d'amortissement C strictement proportionnelle, les coefficients i se
déduisent des termes diagonaux de la matrice d'amortissement généralisée T C par :
T
i C i
2 i i =
T
i M i
et, dans le cas de modes propres normés à la masse modale unitaire,
T
2 i i = i C i
On peut utiliser cette relation dans le cas d'une matrice d'amortissement C non proportionnelle, en
appliquant l'hypothèse de BASILE, qui est acceptable pour des amortissements faibles (notamment s'il
n'y a pas d'amortissement localisé dominant) et des modes propres réels suffisamment découplés.
Les coefficients d'amortissement modaux peuvent être fournis par commande (mot clé AMOR_REDUIT)
à deux opérateurs pour :
· l'analyse transitoire dans l'espace modal
DYNA_TRAN_MODAL
[U4.53.21]
· l'analyse sismique par spectre d'oscillateur
COMB_SISM_MODAL
[U4.84.01]
Notons qu'il n'existe aucun outil d'extraction automatique de ces coefficients, à partir de la matrice
d'amortissement généralisée T C ), concept produit par l'opérateur PROJ_MATR_BASE [U4.63.12].
2
Modèle d'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. Le coefficient d'amortissement hystérétique est
déterminé à partir d'un essai sous chargement cyclique harmonique à la pulsation pour lequel on
obtient une relation contrainte-déformation qui permet de définir :
· l'énergie dissipée par cycle sous la forme :
E = d
d
cycle
· le module d'YOUNG complexe E* à partir de la relation contraintes-déformations :
jt
j(t
- )
= e
0
et = e
0
avec 0 et
amplitudes
les
,
0
phase
la
0
j
E* =
=
e
=
0
(cos + j sin )
0
0
où E* = E1 + j E2 = E1 (1+ j )
Manuel d'Utilisation
Fascicule U2.06 : Dynamique
HT-66/02/003/A
Code_Aster ®
Version
6.4
Titre :
Notice de modélisation de l'amortissement mécanique
Date :
08/11/02
Auteur(s) :
J. PELLET Clé
:
U2.06.03-B Page
: 5/6
avec
0
0
1
E =
(cos )
= partie réelle et E2 =
(sin )
= partie imaginaire
0
0
E
= 1 = tg = facteur de dissipation
E2
Ceci conduit aux équations de la dynamique des structures :
M u&+ K *(1+ j ) u = f ()
éq 2-1
avec K matrice de rigidité élastique réelle, M matrice de masse et le coefficient d'amortissement
hystérétique. Notons que l'on parle souvent de matrice de rigidité complexe.
2.1
Amortissement hystérétique "global"
Cette modélisation, facile à mettre en oeuvre, correspond à :
(-M 2 + j K + K)u = f()
éq
2.1-1
Elle est disponible actuellement, en utilisant l'opérateur COMB_MATR_ASSE [U4.72.01], après avoir
assemblé la matrice de rigidité à coefficients réels, mais elle est d'une utilité faible :
· validation d'algorithmes de résolution,
· inutile 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 (dissipation localisée dans des zones particulières
de la structure traitées avec des matériaux visco-élastiques).
2.2
Amortissement hystérétique des éléments du modèle
2.2.1 Caractéristiques
d'amortissement
Il est possible de construire une matrice de rigidité complexe à partir de chaque élément du modèle,
comme pour la rigidité réelle et la masse.
Deux fonctionnalités sont utilisables :
· l'affectation d'éléments discrets, sur des mailles POI1 ou SEG2, par l'opérateur
AFFE_CARA_ELEM [U4.42.01]. Celui-ci permet de définir, avec plusieurs modes de
description possibles, une matrice de rigidité réelle pour chaque degré de liberté et un
coefficient d'amortissement hystérétique à appliquer à cette matrice.
AMOR_HYST :
éta
[R]
· la définition d'une caractéristique d'amortissement pour tout matériau élastique par l'opérateur
DEFI_MATERIAU [U4.43.01] par le mot clé :
AMOR_HYST :
éta
[R]
ce matériau étant ensuite affecté aux mailles concernées.
Manuel d'Utilisation
Fascicule U2.06 : Dynamique
HT-66/02/003/A
Code_Aster ®
Version
6.4
Titre :
Notice de modélisation de l'amortissement mécanique
Date :
08/11/02
Auteur(s) :
J. PELLET Clé
:
U2.06.03-B Page
: 6/6
2.2.2 Calcul des matrices d'amortissement
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'.
Chaque matrice élémentaire est alors de la forme :
· quand le matériau i, de caractéristiques d'amortissement hystérétique i, est affecté à
l'élément elem
k*elem = kelem ( 1 + j i )
· pour un élément discret défini par une matrice de rigidité kdiscret et un coefficient
d'amortissement hystérétique
k*elem = kdiscret ( 1 + j )
Cette opération est possible avec :
mel
[matr_elem_DEPL_C] = CALC_MATR_ELEM
(
/
OPTION:
'RIGI_MECA_HYST'
MODELE:
mo
[modele]
CHAM_MATER:
chmat
[cham_mater]
CARA_ELEM:
cara
[cara_elem]
RIGI_MECA:
rigi
[matr_elem_*]
CHARGE
:
l_char
[l_char_meca]
);
L'assemblage de la matrice de rigidité complexe K *, à partir des matrices élémentaires est obtenu
avec l'opérateur ASSE_MATRICE habituel [U4.61.22]. On notera que l'on doit utiliser la même
numérotation et le même mode de stockage que pour la matrice de masse (opérateur NUME_DDL
[U4.61.11]).
Le chargement utilisé pour le calcul de la matrice de rigidité réelle (OPTION `RIGI_MECA') doit être
renseigné par le mot clé `CHARGE' pour le calcul de la matrice de rigidité élémentaire complexe.
2.2.3 Utilisation de la matrice de rigidité complexe
La matrice de rigidité complexe K * est utilisable pour l'analyse dynamique linéaire directe (mot-clé
MATR_RIGI) avec l'opérateur de réponse dynamique linéaire :
· harmonique
DYNA_LINE_HARM
[U4.53.11]
Pour la recherche de valeurs propres, aucune fonctionnalité n'est disponible actuellement pour
l'utilisation du modèle d'amortissement hystérique.
Pour les analyses en base modale, aucune fonctionnalité n'est disponible actuellement pour l'utilisation
du modèle d'amortissement hystérétique.
Manuel d'Utilisation
Fascicule U2.06 : Dynamique
HT-66/02/003/A
Document Outline