L’Analyse Modale Expérimentale

- Méthodes et Moyens -

ENST 2013/20014 ATMANI Y.

 Introduction

Actuellement :
z Une vitesse de plus en plus grande dans les transports
z Une plus grande économie dans le carburant

z Ces deux contraintes sont réalisée en réduisant la masse des structures

Conséquences
z Les structures deviennent plus légère
z Les résonances se déplacent aux domaines des fréquences des forces d’excitation
z Les structures cassent à cause des chargements dynamiques
z La statique a fait l’objet d’étude de plus d’un siècle
Les essais dynamiques

Les essais en dynamique des structures

z Terme qui la différentie des essais statique et de l’analyse statique
z L’analyse modale et les essais modaux
z Les déformées opérationnelles (Operational Deflection Shapes)

Objectifs
z Qu’est ce qu’est un Essai Modal
z Pourquoi faire un Essai Modal ?
z Comment faire un Essai Modal ?
 Les essais dynamiques

Analyse Modale contre Déformé Opérationnelle Sortie

z Analyse Modale
–excitation artificielle (marteau, 1 or plusieurs
vibrateurs)
–la force appliqué à l’entrée est bien connue
(“sous notre contrôle”) entrée

z Déformé Opérationnelle
–investigation (visualisation) du comportement de
la structure dans des “conditions de
fonctionnements”
–les forces d’entrée ne sont pas connus et
compliqués à évaluer
C’est quoi?
Pourquoi?

Comment?
Définition de l’Essai Modal

C’est construire un modèle mathématique

décrivant le comportement et les propriétés
vibratoires d’une structure par des moyens
expérimentaux

Fréquence Naturelle
Amortissement Modal
Résidus
Paramètres Modaux

Théoricien Expérimentateur
z Valeur Propre z Fréquence Naturelle
z Pourcentage d’Amortissement z Facteur de Perte
z Vecteur Propre z Forme Modale

Amplitude Premier
e
tanc Mode Deuxième
Di s
Mode
Troisième
Mode
Fréq
uenc
e
Poutre
Vue e
sur le Vu le
r
D omai
ne F su aine
m l
réque
ntiel Do oda
Force M
Le Modèle Modal – Paramètres Globaux
La Fréquence Naturelle et l’Amortissement Modal

Domaine Fréquentiel Domaine Temporel

1
Taux de décroissance σ =
la bande à 3 dB τ
3dB
=2σ

Temps
Fréquence τ
ω0 T

1
Fréquence Naturelle :ω0 = 2πf0 Fréquence Naturelle : f0 =
T
σ σ
Le Taux d’Amortissement: ζ = Le Taux d’Amortissement: ζ =
ω0 ω0
 Le Modèle Modal – Paramètres Locaux
Résidus
Résidu: la “rigidité” du mode

H(ω0)

Amplitude Premier 2σ
n ce Mode
ta Deuxième
Dis Mode
Troisième
Mode ω0 Fréquence
Fréq Résidu: R = H(ω0) · σ
uenc
e

Vue u e Résidu du PE : Riir = a · ψir2

sur le V le
r
D omai
ne F su aine
m l
réque
ntiel Do oda Résidu général : Rijr = a · ψir · ψjr
Force M

Poutre Le Résidu du Point

Entraînant normalise la
Déformé Modale
C’est quoi

Pourquoi?
Comment
 Résolution des Problèmes

Fonction de Réponse en Fréquence

Réponses élevée

Réponse Vibratoire en fonctionnement

SDM et FRS (en supposant que le Modèle Modale est validé)
Scénarios du type
Il Arrive Quoi Si
Modification de la Dynamique Structurelle
Structural Dynamics Modification
z Masse
z Rigidité
z Amortisseur réglé
z Déplacer la fréquence de résonance

Simulation de la Réponse Forcée Réponse élevée Réponse basse

Réponse
Forced Response Simulation précédente
z Comment une structure va se
comporter après qu’une ou plusieurs
forces soit appliquées?
 Pourquoi faire des essais Modaux?
z Résoudre les problèmes
–Pour réduire les niveaux excessives de vibration
z Modélisation-Éléments Finis
–Pour s’assurer que les fréquences de résonances sont éloignés de celle de
l’excitation
–Validation par tests sur prototypes
–Raffinement du model mathématique en introduisant l’amortissement
–Condition préalable dans l’industrie aéronautique
–Aujourd'hui aussi couramment utilisé dans l’industrie automobile
z Analyse des assemblage des structure
–Prédire le comportement dynamique de sous composants assemblés
z Simulation des scénarios type “Il Arrive Quoi Si”
–Détermination des forces
–Réponse au excitations complexes
C’est quoi

Pourquoi

Comment?
Comment conduire un Essai Modal ?
Modélisation
z Géométrie
1 z Définition des Degrés De Libertés
z Direction X ou XYZ
Mesures
2 z Excitation par marteau ou vibrateur
z Fonctions de Réponse en Fréquence
z Fonction de Cohérence pour validation
Lissage des courbes
z Fréquence
3 z Amortissement
z Résidus
Validation
z MAC (Modal Assurance Criteria)
4 z MCF (Modal Confidence Factor)
z Dispersion de la phase (Phase Scatter)
z Diagramme de stabilité
Modélisation

f(t)
Les modèles SDOF et MDOF
m x(t)

c k Différentes techniques d’Analyse

Modale

Exciter une structure

Mesurer correctement les données

= + + +⋅⋅⋅+
Le post-processing des données de
l’Analyse Modale
 Analyse du Système

Excitation a(t)

Réponse b(t)

h(τ) Fonction de Réponse Fonction de Réponse

Impulsionelle |H(f) |
en Fréquence

temps

fréquence
 Les Descripteurs d’un Système

a(t) h(τ ) b(t) A( f ) H( f ) B( f )

Fonction de Réponse Impulsionelle Fonction de Réponse en Fréquence

Convolution Multiplication
+∞
b(t) = ∫ h(τ )a(t −τ )dτ = h(t) ∗ a(t) B( f ) = H( f ) ⋅ A( f )
−∞

h(τ ) et H( f ) Sont des descripteurs de système

indépendant des signaux introduits
Analyse du Système
Détermination des propriétés intrinsèques du système

z Excitation du système par une force connue

Sortie z Mesure de la Sortie
z Faire correspondre la Sortie à l’Entrée

Fonction de Réponse en Fréquence :

Sortie Mouvement Réponse

H(ω) = = =
Entrée Force Excitation

z La FRF montre les propriétés

intrinsèques d’un système
Entrée dynamique, indépendamment du
type de l’excitation et de la force de
l’excitation
Définition de la Fonction de Réponse en Fréquence
F(f) H(f) X(f)

F H X

f ∠H f f

H(f ) = X(f )
F(f )
H(f) est la Fonction de Réponse en Fréquence du système
F(f) est la Transformé de Fourier de l’entrée f(t)
X(f) est la Transformé de Fourier de la sortie x(t)
 Transformations

• Transformé de Laplace
∞

∫
− st
G (s) = g (t )e dt
0

s e s t u n e v a ra ib le c o m p le x e
s = σ + jω

• Transformé de Fourier
+∞
G (ω ) = ∫
− jω t
g (t )e dt
−∞

p o u r g (t ) c a u s a le i.e . g (t ) = 0 , t < 0
G ( jω ) = G ( s ) s = jω

• Ne pas confondre la variable σ ( la partie réelle de s) avec la constante

du système σ appelé constante de décroissance !
 FRF contre Fonction de Transfert

Fonction de Réponse de Fréquence = Fonction de transfert le long de l’axe jω

R R∗
H (s) s = jω = +
jω − p jω − p ∗
ou
R R∗
H (s) = +
jω − ( − σ + jω d ) jω − ( − σ − jω d )
s = jω
Modèle de Laplace pour un système SDOF
f(t)
m x(t)
M&
x&(t ) + Cx&(t ) + Kx (t ) = f (t )
c k

⎡⎣ ms 2 + cs + k ⎤⎦ X ( s ) = F ( s )
z La Fonction de Transfert
1
X (s) 1 m
H (s) = = =
F ( s ) ms 2 + cs + k ( s − r ) ( s − r ∗ ) 1
r , r ∗ = −σ ± jω d sont les pôles p de H(s) I.e. H ( s ) = m
( s − p ) ( s − p∗ )
∗
H(s) peut être développé dans une forme de fraction R R
partielle: H (s) = +
s − p s − p∗
Avec R est le résidu du pole p: 1 1
R = H ( s ) ( s − p )s = p = m = m
p − p ∗ 2 jω d
Lieux des Pôles – Plan de Laplace
jω

ω 0 = La fréquence naturelle non-amortie

p = −σ + jω d ω d = La fréquence naturelle amortie

θ
σ σ
taux de décroissance σ
cos(θ ) = ζ =
p * = −σ − jω d ω0
le ratio d'amortissement
 La forme simplifié d’un système SDOF
Déplacement
Déplacement
d = D sinωnt
D

Time

T Fréquence
1
T

m
Période, Tn in [sec]
k 1
Fréquence, fn= T in [Hz = 1/sec]
n

k
ωn= 2 π fn =
m
Système Masse et Ressort

temps

m1
k
ωn = 2πfn =
m + m1 m

Augmenter la masse
réduit la fréquence
Système Masse, Ressort et Amortisseur

temps

Augmenter
l’amortissement
m
réduit l’amplitude

k c1 + c2
 Le Modèle basique SDOF

f(t)

m x(t)

c k

M&
x&(t ) + Cx&(t ) + Kx (t ) = f (t )
&
x&( t ) = Vecteur Accélération
M = masse (force/accé.)
x&( t ) = Vecteur Vitesse
C = amortissement (force/vite.)
K = rigidité (force/dépl.) x( t ) = Vecteur Déplacement
f (t ) = Vecteur force appliqué
 Les Modèles SDOF:
Domaine Temporel et Fréquentiel
F(ω) H(ω) X(ω)

|H(ω)| 1 1
ω2m ωc
f(t) 1
k

m x(t)
ω
c k ∠ H(ω) ω0 = √ k/m
0º ω
– 90º
– 180º

X (ω ) 1
H (ω ) = =
&+ cx&+ kx = f
mx& F (ω ) − ω 2m + jωc + k
 Le Modèle SDOF: Domaine Fréquentiel
f(t) |H(ω)| 1 1
ω2m ωc
x(t) 1
m k

c ω
k ∠ H(ω) ω0 = √ k/m
0º ω
– 90º
– 180º

m a = f − cv − kx ⇔ m &
x&+ cx&+ kx = f
X (ω ) 1 R R∗
H (ω ) = = 2 = +
F (ω ) −ω m + jωc + k jω − ( −σ + jωd ) jω − ( −σ − jωd )
1 k c2 c
avec R= ωd = − 2 σ=
j 2mωd m 4m 2m
Le Modèle SDOF: Domaine Fréquentiel

R R∗
H (ω ) = − avec p = −σ + jωd
jω − p jω − p ∗

R
A la résonance : H (ω ) ≈
σ
Réel et Imaginaire [Link] & Magnetidue &
Nyquist
Phase Phase

Im
|H| |H|
Re
Re

θ θ
Im
 La Largeur de Bande de la Résonance
( estimation de R, ωd et σ)
R R∗
H ( jω ) = +
jω − ( −σ + jω d ) jω − ( − σ − jω d )
R
FRF dans la proximité de ωd : H ≈
jω − ( −σ + jω d )
z

R
z FRF à ω =ωd : H ≈
σ -3dB
R
z FRF à ω =ωd+σ : H ≈
jσ + σ
R
z FRF à ω =ωd-σ : H ≈
− jσ + σ
ωd-σ ωd ωd+σ

R H ( jω d )
H ( jω ± jσ ) = ≈ alors ∆ ω 3 dB = 2σ
2σ 2
 Résidu
( estimation de R, ωd et σ)
z Sens du Résidu, R = rigidité du pole
z Une estimation de R nécessite une estimation précise de σ
Log Magnitude
100 ∆ ω 3 dB
R = H ( jω d )σ = H ( jω d ) ⋅
10−2 2
10−4

10−6

10−8 Résidu ~ la surface sous la FRF

10−10

∆ω3dB/2 Fréquence
0
σ
Phase le ratio d'am ortissem ent ζ =
ω0
−180
Fréquence
 Matrice Modale
⎧ X 1 (ω )⎫ ⎧ ⋅ ⎫
⎪ X (ω )⎪ ⎡ H11 (ω ) H 21 (ω ) ⋅ ⋅ H1n (ω )⎤ ⎪
⋅ ⎪
⎪ 2 ⎪ ⎢ ⋅ H 22 (ω ) ⋅ ⎥ ⎪ ⎪
Modèle Modale (
⎪X3 ω ⎪ ) ⎢ ⎥ ⎪ F3 (ω )⎪
⎨ ⎬=⎢ ⋅ H 23 (ω ) ⋅ ⎥⎨ ⎬
(Domaine Fréquentiel) ⎪ ⋅ ⎪ ⎢ ⋅ ⎥⎪ ⋅ ⎪
⋅ ⎥
⎪ ⋅ ⎪ ⎢ ⎪ ⋅ ⎪
⎪ ⎪ ⎢⎣ H n1 (ω ) ⋅ ⋅ ⋅ H nn (ω )⎥⎦ ⎪ ⎪
⎩ n ⎭
X (ω ) ⎩ ⋅ ⎭

⋅
⋅
X4
X3
X2
X1
H22
F3
H21
Modèle MDOF
Amplitude

1+2

d1+ d2
2 1
Fréquence
m
d1
Phase
Fréquence
dF 0°
1 2
-90°
1+2

-180°
Quel est le problème avec les Modèles Pysiques ?

Coordonnées Physiques = CHAOS Espace Modal = Simplicité

Γ1

Rotor 1 q1

2σ1 ω01
2

Roulement Roulement Γ2

1 q2

Bâti 2σ 2 ω02
2
Γ3

1 q3

2σ3 ω03
2
Modèle MDOF: Domaine Fréquentiel
f(t)

m x(t) • sous la forme d’une fraction partielle:

k c

X (ω ) 1 R R∗
H (ω ) = = = +
F (ω ) −ω m + jω c + k
2
jω − ( −σ + jω d ) jω − ( −σ − jω d )
1 k c2 c
avec R= ωd = − σ=
j 2 mω d m 4m 2 2m

⎛ R ijr ∗
R ijr ⎞
H ij (ω ) = ∑H ijr =∑⎜
⎜
r ⎝ jω − p r
+
jω − p r
⎟⎟
r ⎠
Le Comportement Modal d’un système MDOF
• La réponse dynamique d’une structure peut être décomposé en un ensemble
discret de mouvement particulier indépendant. Ces mouvements sont
appelées les Modes de Vibrations.
• Un mode est décrit par les paramètres modaux :
− la Fréquence naturelle et l’Amortissement ( Valeurs Propres)
− Déformé Modale ( Vecteurs Propres)

= + + +⋅⋅⋅+

{x ( t )} = q 1 ( t ) {φ }1 + q 2 ( t ) {φ }2 + q 3 ( t ) {φ }3 + . .. + q n ( t ) {φ }n
 Le Comportement Modal d’un système MDOF
m
{ x } = [ Φ ]{q } = ∑ {φ }r q r
r =1

Cette relation montre que la déflection dynamique {x} est une combinaison
linéaire des déformés modales, les coefficients {φ} sont les déplacements
modaux

= + + +⋅⋅⋅+
Le Modèle Dynamique Modal du MDOF
le Modèle de la Réponse
{ x } = [ H ]{ F } [H] peut être représenté en terme de Propriétés Modales

⎛ a {φ } {φ } T a ∗ {φ } ∗ {φ } T ⎞
∗
m
[ ] ∑
H = ⎜ r r r
+
r r r ⎟
r =1 ⎜
⎝
( jω − p r ) ( jω − p r ∗ ) ⎟
⎠
m a r φ ir φ jr a r ∗φ ir φ jr
H ij = ∑ +
r =1 ( jω − p r ) ( jω − p r )
∗

• Le dénominateur dans toutes les FRF est le même et il contient la fréquence

Modale et l’Amortissement
− la Fréquence Modale et l’Amortissement sont des paramètres globaux
• Le numérateur est constant , et pour chaque mode il dépend de la déflection
modale de la réponse et du point d’excitation
− l’Amplitude est un paramètre local et dépend de la déformé modale
• Chaque ligne et colonne dans |H| contient la totalité de la déformé modale
Les Formes Alternatives du Modèle Modal
m
1 mr
H ij ( f ) = ∑
1 − ( f f r ) + 2 jζ r ( f f r )
2
r =1

m ar Ψ ir Ψ jr ar ∗ Ψ ir ∗ Ψ jr ∗
H ij (ω ) = ∑ +
r =1 jω − pr jω − pr ∗
m aijr Rijr ∗
H ij (ω ) = ∑ +
r =1 ( jω − pr ) ( jω − p ) r
∗

m Rijr Rijr ∗
H ij (ω ) = ∑ +
r =1 ( jω − pr ) ( jω − p ) r
∗

m uijr + jvijr uijr + jvijr

H ij (ω ) = ∑ +
r =1 −σ r + j (ω − ω dr ) −σ r + j (ω + ω dr )
m Ψ ir Ψ jr
H ij (ω ) = ∑
r =1 −ω 2 + 2 jσ rω + ω 0 r 2
m 2 RijrR jω + 22 RijrR σ r − 2 RijrI ω d
H ij (ω ) = ∑
r =1 −ω 2 + 2 jσ rω + ω 0 r 2
avec Rijr = RijrR + RijrI est le résidu comlexe et pr = −σ r + jω dr est le lieu du pole
Conclusion sur les systèmes MDOF

Un modèle MDOF = Σ modèles SDOF

Résidus et Déflections Modales
⎛ R ijr ∗
R ijr ⎞
H ij (ω ) = ∑H ijr =∑⎜
⎜
r ⎝ jω − p r
+
jω − p r
⎟⎟
r ⎠
|Hij|

ω
R ijr = a r φ ir φ jr

R éponse en i
H ij =
Force en j

1 2 3 4 5 6
 Résidus et Déflections Modales
m
[ R ]r [ R ∗ ]r
{ X (ω )} = [ H (ω ) ]{ F (ω )} H (ω ) = ∑
r =1 jω − p r
+
jω − p r ∗
[ R ] = a r {φ }r {φ }r
T
Hypothèse :
⎡ φ 1φ 1 φ 1φ 2 ... φ 1φ j ... φ 1φ n ⎤
⎢φ φ φ 2φ 2 ... φ 2φ j ... φ 2φ n ⎥⎥
⎢ 2 1
⎢ ... ... ... ⎥
i.e. [ R ]r = ar ⎢ ⎥
⎢ ... ... ... ⎥
⎢ ... ... ... ⎥
⎢ ⎥
⎢⎣φ nφ 1 φ 1φ 2 ... φ nφ j ... φ nφ n ⎥⎦
Chaque colonne et ligne dans [R]r contient la déformé modale du
mode r multiplié par un élément lui appartenant
Toute l’information dans [H(ω)] est contenue dans
une SEULE colonne ou ligne
Pics en quadrature

Pour des modes légèrement couplées la FRF ( la compliance) à la

résonance est aproximé par un SDOF:

R ijr a r φ iφ j ⎧⎪ ∝ φ i pour une colonne

H ij (ω dr ) ≅ = =⎨
σr σr ⎪⎩ ∝ φ j pour une ligne
A la résonance le déplacement( ou l’accélération) est déphasé de
90° par rapport de l’excitation, i.e. Hij(ωdr) est imaginaire pur
( et Rijr est imaginaire pur aussi)
L’ Amplitude de la partie imaginaire de la FRF ( compliance ou
accélérance) est proportionnelle à la déflection modale

Ιm [Hij(ω)]

ω
Déformée Modale à partir des Pics en Quadrature

#5
#4
#3
#2
DDL # 1
Fonction de Réponse en Fréquence
[m/s²] Time(Response) - Input
Working : Input : Input : FFT Analyzer

80

40

-40

-80

0 40m 80m 120m

[s]
160m 200m 240m
FFT
Sortie Domaine Fréquentiel Domaine Temporel
[m/s²] Autospectrum(Response) - Input
Working : Input : Input : FFT Analyzer
10

100m
1
FFT
[(m/s²)/N] (Magnitude))
Frequency Response H1(Response,Excitation) - Input (Magnitude

Inverse
Impulse Response h1(Response,Excitation) - Input (Real Part)
[(m/s²)/N/s]
Working : Input : Input : FFT Analyzer Working : Input : Input : FFT Analyzer

⇒
10m 2k
10

1m 1k
0 200 400 600 800 1k 1,2k 1,4k 1,6k
[Hz]
0
100m

-1k
[N] Autospectrum(Excitation) - Input

Entrée
Working : Input : Input : FFT Analyzer -2k
1 0 200 400 600 800 1k 1,2k 1,4k 1,6k 0 40m 80m 120m 160m 200m 240m
[Hz] [s]
100m

10m

Fonction de Fonction de
1m

100u
0 200 400 600 800 1k 1,2k 1,4k 1,6k

Réponse en Réponse
[Hz]

Fréquence Impulsionelle

[N] Time(Excitation) - Input

Working : Input : Input : FFT Analyzer
FFT Sortie Mouvement Réponse
= =
200

100
H(ω) =
Entrée Force Excitation
0

-100

-200
0 40m 80m 120m 160m 200m 240m
[s]
 Les avantages de la
Fonction de Réponse en Fréquence

F(f) H(f) X(f)

z Les Fonctions de Réponse en Fréquence sont des propriétés des systèmes

dynamiques linéaires

z Elles sont indépendantes de la Fonction d’Excitation

z L’Excitation peut être une fonction du temps Périodique, Aléatoire ou

Transitoire

z Les résultat obtenu avec un seul type d’excitation peuvent être utilisé pour
prédire la réponse du système pour n’importe quel autre type d’excitation
Les différentes Formes de la FRF

Compliance Rigidité Dynamique

(déplacement / force) (force / déplacement)

Mobilité Impédance
(vitesse / force) (force / vitesse)

Inertance ou Réceptance Masse Dynamique

(accélération / force) (force /accélération)
Les Différents Estimateurs des FRF

F(f) H(f) X(f) H( f )= X ( f )

F( f )

H1( f ) = GFX ( f )
GFF ( f )

H 2 ( f ) = GXX ( f )
GXF ( f )

G XX GFX
H3( f ) = ⋅ = H1 ⋅ H 2
GFF GFX

2
GFX G G* H
γ ( f )=
2
= FX ⋅ FX = 1
GFF ⋅GXX GFF GXX H 2
 Leakage dans la FRF
• Biais de résolution

Erreur de biais de résolution

Bs largeur de bande du système

τs temps de réponse du système ≈ 1 Bs
B la largeur de bande de mesure
B τ longeur de l'enregistrement ≈ 1 B
B B

si B >> Bs ⇒ leakage dans la fréquence

coresspondant à T << τ s ⇒ troncation dans le temps
Comment éliminer le leakage dans les FRF:
• Elimination par zoom:
B << Bs
co rresp o n d an t à
T >> τ s

• Minimisation à la résonance par l’emploi de H2 :

H > H 2 > H1

• Minimisation dans les anti-résonance par l’emploi de H1 :

H 2 > H1 > H
Quel Estimateur de FRF Doit On Utiliser?
Précision
G XX GFX
Définitions : H1( f ) = GFX ( f ) H 2 ( f ) = GXX ( f ) H3( f ) = ⋅
GFF ( f ) GXF ( f ) GFF GFX

Précision pour des systèmes avec: H1 H2 H3

Bruit à l’entrée - Meilleur -
Bruit à la sortie Meilleure - -
Bruit à la entrée+sortie - - Meilleure
Pics (leakage) - Meilleure -
Vallées (leakage) Meilleure - -

L’utilisateur peut choisir H1, H2 or H3 après la mesure

Mesures 1

f(t)
Les modèles SDOF et MDOF
m x(t)

c k Différentes techniques d’Analyse

Modale

Exciter une structure

Mesurer correctement les données

= + + +⋅⋅⋅+
Le post-processing des données de
l’Analyse Modale
Trois Types d’Analyse Modale

1. Essai par Marteau

– Coups par impact du marteau...mesures sériels ou parallèles
– Excite un large plage de fréquence rapidement
– La technique la plus utilisé

2. Essai par Vibrateur

– Vibrateur pour excitation modale... mesures sériels ou parallèles
– Plusieurs types d’excitation
– Utilisé plutôt pour les structures complexes

3. L’Analyse Modale Opérationnelle

– Utilise l’excitation naturelles des structures... mesures sériels ou
parallèles
– Une technique très avancée!
Les Différents Types d’Analyse Modale
(Avantages)

1. Essai par Marteau

– Rapide et facile
– Typiquement pas cher
– Peut être réalisé par un personnel pas très qualifié
2. Essai par Vibrateur
– Plus répétable qu’un essai par marteau
– Plusieurs types d’entrées sont disponibles
– Peut être utilisé pour des analyses type MIMO
3. Analyse Modale Opérationnelle
– Pas de besoin pour des conditions de frontières particulières
– Les mesures peuvent être effectué in-situ
– Utilise l’excitation naturelle de la structure en fonctionnement
– On peut effectuer d’autre tests lorsque on collecte les données pour un essai OMA
 Les Différents Types d’Analyse Modale
(Inconvénients)
1. Essai par Marteau
–Les mesures impulsionelles introduisent des facteurs de crête important
–La force d’entrée peut être différente d’une mesure à l’autre ( des opérateurs différents, des
endroits difficiles, etc.)
–Des coups ‘Calibré’ nécessite beaucoup de précautions ( coups doubles, etc.)
–La tête du marteau peut causer parfois une décroissance de la performance ( un niveau de
signal faible)
2. Essai par Vibrateur
–Un réglage de l’essai plus difficile (fixateurs, excitateur, etc.)
–Usuellement plus d’équipements et de canaux sont requis
–Des opérateurs plus qualifié sont requis
3. L’Analyse Modale Opérationnelle
–Le Modèle Modal n’est plus à l’échelle
–L’excitation est supposé couvrir tous le domaine fréquentiel auquel on s’intéresse
–Des données temporelles très grandes sont parfois exigée
–Très coûteuse en ressource informatique ( puissance et temps processeur)
Mesures 2

f(t)
Les modèles SDOF et MDOF
m x(t)

c k Différentes techniques d’Analyse

Modale

Exciter une structure

Mesurer correctement les données

= + + +⋅⋅⋅+
Le post-processing des données de
l’Analyse Modale
 Les Techniques d’Excitation
Excitation par Marteau Excitation par Vibrateur
Excitation en déplacement Réponse Multicanal ou les
points de réponse peuvent se déplacer
In Out

In
Out

z Les petites structures homogènes z Les structures grosses ou complexes

z Une technique polyréférence z Plusieurs types de signaux d’excitations
z Méthode rapide – pas de fixations sont possible
z Consommation de temps – temps
⎡H11 H12.........H1n ⎤ d’installation ⎡H11.......... ... ⎤
⎢ ⎥ ⎢ ⎥
.......... .......... .. H .......... ...
[H] = ⎢ ⎥ [H] = ⎢⎢ 21 ⎥
⎢.......... .......... .. ⎥ .......... .........⎥
⎢ ⎥ ⎢ ⎥
⎣.......... .......... .. ⎦ ⎣Hn1.......... ... ⎦
Essai par Marteau sur une poutre libre-libre
Méthode du marteau itinérant :
z Réponse mesuré dans un seul point
z Excitation de la structure par un marteau instrumenté avec un
capteur de force dans plusieurs points
z Les FRF’s sont calculé entre le point de mesure et les points d’excitation
z Les Modes de la structure sont identifié
Amplitude Premier
e Mode
tanc Second
Di s Mode
Troisième
Mode
Fréq Poutre
uenc
e Accélération
Force
Force
Vue Force le e
sur l Force
Force r
eD oma Force
Force e s ain al
u
ine Force m
Fréq
uent Force
Force Vu Do Mod
iel Force
Force

Press anywhere
to advance animation
 Mesure de la Matrice FRF (SISO)
Simple Input Simple Output

Une ligne
z Une seule Excitation intérante
z Une seule Réponse Fixe (référence)
SISO

X1 H11 H12 H13 ...H1n F1

X2 H21 H22 H23 ...H2n F2
X3 = H31 H32 H33...H 3n F3
: : :
Xn Hn1 Hn2 Hn3...Hnn Fn
 Mesure de la Matrice FRF (SIMO)
Simple Input Multiple Output

Plusieurs lignes
z Une seule Excitation intérante
z Plusieurs Réponses Fixes (référence)
SIMO

X1 H11 H12 H13 ...H1n F1

X2 H21 H22 H23 ...H2n F2
X3 = H31 H32 H33...H 3n F3
: : :
Xn Hn1 Hn2 Hn3...Hnn Fn
Essai par Vibrateur sur une poutre libre-libre
La méthode par Vibrateur:
z Excitation de la en un seul point par vibrateur instrumenté par un
capteur de force
z Réponse mesuré dans de nombreux of points
z Les FRF’s sont calculé entre le point de mesure et les points d’excitation
z Les Modes de la structure sont identifié
Amplitude Premier
e
tanc Mode Second
Di s Mode
Troisième
Mode
Fréq Poutre
uenc
e Accélération
Vue
sur l
eD ur leine
oma
ine ue soma dal
Fréq
uent
V D Mo
iel Excitation
par ∝ Force
Bruit Blanc

Press anywhere
to advance animation
Mesure de la Matrice FRF (Vibrateur en SIMO)

Une colonne
z Une seule Excitation Fixe (référence)
z Une seule Réponse itinérante SISO
ou
z Plusieurs Réponses(Itinérantes) SIMO
Multiple-Output:
Optimiser la consistance des données

X1 H11 H12 H13 ...H1n F1

X2 H21 H22 H23 ...H2n F2
X3 = H31 H32 H33...H 3n F3
: : :
Xn Hn1 Hn2 Hn3...Hnn Fn
Pourquoi le ” Multiple-Input and Multiple-Output” ?

z Plusieurs-Entrées: Pour les structures grandes et/ou complex plusieurs

vibrateurs sont néçessaire pour:
–distibuer efficacement l’energie d’excitation
et
–éviter le comportement non-linéaire

z Plusieurs-Sorties: Mesurer les sorties au meme temps pour optimiser la

consistence des données
i.e. MIMO
Les situations dans lesquels le MIMO est nécessaire
z Une seul ligne ou une seule colonne n’est pas suffisante pour
déterminer tous les modes dans les situations suivantes:

–Plusieurs modes à la meme fréquence (racines multiples),

e.g. structures symmetriques

–Structures complexes ayant des modes locaux, i.e. Une

DDL référence avec une déformé modale qui ne correspond
pas à tous les modes
z Dans tous les cas plus de colonnes ou plus de lignes
hdoivent etre mésurée - i.e. polyréférence.
Solutions:
–Une excitation par Impact par le Marteau avec plus de DDLs
de réponse
–Un vibrateur déplacé dans de differents DDLs de référence
–MIMO
Mesure de la Matrice FRF( Vibrateur MIMO)
Plusieurs Colonnes
z Plusieurs Excitations Fixes (références)
z Une Seule Réponse Itinérante MISO
ou
z Plusieurs Réponses Itinérantes MIMO

X1 H11 H12 H13 ...H1n F1

X2 H21 H22 H23 ...H2n F2
X3 = H31 H32 H33...H 3n F3
: : :
Xn Hn1 Hn2 Hn3...Hnn Fn
La Question Éternelle dans l’Analyse Modale …

F1
a

F2
 Excitation par Impact
Mesurer une ligne de la matrice FRF en
changeant la position du point d’impact

H11(ω) H12(ω) ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ H15(ω)

⋅⋅⋅⋅⋅
#5 ⋅⋅⋅⋅⋅
Accéléromètre #4 ⋅⋅⋅⋅⋅
#3 ⋅⋅⋅⋅⋅
#2 ⋅⋅⋅⋅⋅
#1

LAN

Capteur
de Force Marteau
d’Impact
 Excitation par Impact
a(t)

z La durée de l’impulsion et son amplitude dépendent du:

– Poids du marteau
– La tête du marteau (acier, plastique ou caoutchouc)
– Les caractéristiques dynamiques de la surface
– La vitesse à l’impact

z La largeur de la bade de fréquence est inversement proportionnelle à la

durée de l’impulsion
a(t) 1 GAA(f) 1
2 2

t f
 Les Fonctions de Fenêtrage pour
l’Excitation par Impact
Critère
z Comment sélectionner la
(shift) et la longueur des
fenêtres transitoire et
exponentielle :
Fenêtrage Transitoire du signal d’entrée

Fenêtrage Exponentiel
du signa de sorite

z Leakage due to exponential time

weighting on response signal is well
defined and therefore correction of the
measured damping value is often possible
Compensation du Fenêtrage Exponentiel
b(t)
Avec une fenêtre Fonction de fenêtrage
1
exponentielle Signal Original
appliqué au signal de
sortie, la constante
de temps mesurée
est très petite et la
constante de Temps
décroissance (par
Signal fenêtré
conséquent le taux
d’amortissement) est shift Longueur τ W
trop grand Longueur de l’enregistrement, T

Correction de la constante de décroissance σ et du taux d’amortissement ζ:

σ σm σ W
σ = σm − σ W ζ= = − = ζm − ζ W
ω0 ω0 ω0
Valeur Valeur 1
corrigé
σW =
Mesurée τW
Gamme des Marteaux
Déscription Application

grand marteau Bâtiments et

ponts
(4 kg)

marteau de Grand arbres et

main (1kg) grands
outillages
Marteau (0.5 Structure
kg) d’automobile et
outillages

Utilisation Composants
générale( 0.1
lb)
Mini marteau Disques durs,
circuits
imprimés, aubes
de turbines
Excitation par Marteau d’Impact

z Avantages:
–rapide
–Pas de problème de fixation
–Pas de masse de chargement
variable
–Portatif et idéal pour les essais sur
terrains
–La moins cher relativement aux
autres
z Inconvénients:
–Des facteurs de crête très grand
peuvent entraîner la structure a avoir
un comportement non-linaire
z Conclusion –Une force avec des pics très
important nécessaire pour exciter les
– très adapté au terrain grandes structures peuvent
– Utile pour déterminer le vibrateur occasionner des dommages locaux!
et les endroits de fixation –Un signal hautement déterministe
conduit a une approximation non
linéaire
 Excitation par Vibrateur
Mesurer une colonne de la matrice FRF en
changeant la position du capteur de la réponse

H11(ω) ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅
Accéléromètre H21(ω) ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅
#5 ⋅
#4 ⋅
#3 ⋅
#2 H51(ω) ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅
#1
Capteur de
Force

LAN

Vibrateur Amplificateur
de
Puissance
Fixation des Capteurs et du Vibrateur
Capteur de
a F Force Vibrateur
Montage de l’accéléromètre :
z Goujon fileté Accéléromètre

z Ciment
z Colle
z Base magnétique Propriétés de la Tige
Rigidité Axiale : Grande
Rigidité de Flexion : Petite
Capteur de Force et Vibrateur:
z Goujon fileté Avantages de la Tige:
z Pas de Moment d’Excitation
z Tige (Tige de connexion)
z Pas de Charge par Inertie de Rotation
z Protection du Vibrateur
z Protection du Capteur
z Aide a positionner le Vibrateur
 Fixation de l’Excitateur et de la Tige
Capteur de Force
Structure a mesurer Tige mince Excitateur

Accéléromètre

Les mesures de la force et de l’accélération ne sont pas affecté par la compliance de la tige,mais ...
Une correction Mineure de la masse est requise pour déterminer l’excitation effective
Fm
Fs
&
&
Fs = (m + M) X
Structure

Masse de la pointe &

&
Fm = M X
m
Matériau Masse du m+M
Piézoélectrique Vibrateur/Marteau, M Fs = Fm
M
La Calibration par une Masse Libre
&
X&(f ) &
X&(f ) 1
A(f ) = =
&
=
&(f ) m
F (f ) mX
Capteur de Force
Excitateur

Accéléromètre
Mass
e

⎛X&&(f ) ⎞ ⎛X&&(f ) ⎞
Mag ⎜ ⎟ ∠⎜ ⎟
⎝ F (f ) ⎠ ⎝ F ( f ) ⎠
Mag=1/masse Phase=0°

fréqeunce fréqeunce
 Tiges de fixations :
la liaison maîtresse pour la précision
• Utilité:
− Pour transmettre uniquement la force dans la direction sensible du capteur de force / tête
d’impédance
− Fusible mécanique-protection de la structure et de l’excitateur
− Rend l’alignement structure/excitateur plus aisé

• Tige idéale: rigidité de flexion nulle , rigidité axiale infinie, masse nulle
• Réalisation pratique:
− câble de tension (caractéristiques : rigidité de flexion nulles, masse faible)
− Push/pull (caractéristiques : une certaine rigidité de flexion, plus de masse)

• Conclusion:utiliser les Câbles de tensions de préférence !

• Limitations du câble de tension:
− on ne doit pas bouger la structure lors de la pré-charge!
− la pré-charge doit être effectué et vérifié
−La mise au point nécessite plus de temps
− Force limité (pic en sinus de < 400 N) et en fréquence ( <400 Hz)
 Comment obtenir des FRF de bonnes qualité:
quelques aspects sur le signal d’entrée
Minimiser l’interaction entre Choisir la meilleure position de
Excitateur/Structure à Tester l’excitateur
Solution : une armature légère en Solution : minimiser au maximum le
magnésium poids total de l’excitateur et des
fixations
Obtenir des mesures de force valides
Test modal basse fréquence performant
Solution : tige avec la plus petite
valeur possible de rigidité de Solution : augmenter au maximum le
flexion déplacement pic à pic
Distribuer et minimiser les niveaux de Capteur
force dynamiques Force/Impédance
Excitateur
Solution : un réglage flexible Structure Modal
a mesurer
Capacité de l’excitateur à encaisser Tige
de grand efforts latéraux
solution : construction renforcée
Mesure
avec 2×4 ressorts
au PE
 La Force de Réaction du Vibrateur

Réaction par Réaction par Exemple d’un

support externe l’inertie de l’excitateur mauvais
montage
Suspension de la
Structure a Tester

Suspension de la
Structure a Tester

Exciter
Support de Suspension
l’Excitateur
Choix de la Technique d’Excitation:

Il est fonction de:

• Application
• Non-linéarités
• La vitesse du test
• Leakage (analyse)
• Plage dynamique
• Facteur de Crête
• Bruit à l’entrée et/ou la sortie
• L’équipement disponible
Excitation Sinusoïdale
a(t)
A
RMS
A
Facteur de Crête = RMS = 2
Time

B(f1)
z Pour l’étude des non-linéarités, par
exemple: la distorsion harmonique

z Pour une excitation sur une large bande il faut que: A(f1)
– L’onde Sinus balaye lentement la plage
fréquentielle d’intérêt
– La condition de quasi-stationnarité soit satisfaite
Excitation par Sinus Balayé

Avantages
z Facteur de Crête faible
z Rapport Signal/Bruit grand
z La force d’entrée mieux contrôlé
z possibilité d’étudier les non-linéarités

Inconvénients
z Très lente
z Pas d’approximation linéaire des systèmes non-linéaires
Excitation Aléatoire
a(t)

Temps

Variation Aléatoire de l’amplitude et de la phase

⇒ le Moyennage va nous donner une estimation linéaire
optimale dans le cas des non-linéarités

Sortie du
Système
GAA(f), N = 1 GAA(f), N = 10 B(f1)

A(f1)
Fréq. Fréq. Entrée du
Système
Excitation Aléatoire
z Signal aléatoire :
– Caractérisé par sa densité de puissance spectrale (GAA) et ça densité de
probabilité d’amplitude (p(a))
a(t)

p(a) Temps

z Peut être limité dans une bande suivant la plage de fréquence qui nous
intéresse
GAA(f) Bande de base GAA(f) Zoom

Fréq. Fréq.
Bande de Bande de
Fréquence Fréquence
z Le signal n’est pas périodique dans le temps d’analyse time ⇒ problème de
Leakage dans les estimateurs spectraux
Excitation Aléatoire

Avantages
z La meilleure approximation linéaire du système
z Zoom
z Facteur de Crête acceptable
z Rapport Signal/Bruit acceptable

Inconvénients
z Leakage
z Nécessité de Moyennage (lenteur)
Aléatoire Impulsive
z Caractéristiques du signal Aléatoire Impulsive :
– Donne la meilleure approximation linéaire d’un système non-linéaire
– Travaille avec le zoom

a(t)

Temps

Avantages
z La meilleure approximation linéaire d’un système
z Pas de leakage (si un fenêtrage temporel rectangulaire est utilisé)
z Relativement rapide

Inconvénients
z Signal/Bruit et Facteur de Crête pas optimum
z Une pondération temporel spéciale est parfois nécessaire
Excitation Pseudo Aléatoire
z Signal Pseudo Aléatoire :
– Blocks de signaux aléatoires répété chaque T
a(t)

Time
T T T T

z La période de temps est égale à la longueur de l’enregistrement T

– Le spectre de ligne coïncide avec les lignes de l’analyseur
– Pas de pondération des non linéarités Sortie du
Système
GAA(f), N = 1 GAA(f), N = 10
B(f1)

A(f1)
Fréq. Fréq. Sortie du Système
 Excitation Pseudo Aléatoire
z Signal Pseudo Aléatoire :
– Caractérisé par Puissance/RMS (GAA) et densité de probabilité d’amplitude (p(a))
a(t)

p(a)
T T T T Temps

z Peut être limité dans une bande suivant la plage de fréquence qui nous intéresse
GAA(f) Bande de Base GAA(f) Zoom

Fréq. Fréq.
Plage de Fréq. Plage de Fréq.

La période de temps égal à T

⇒ Pas de leakage si on utilise une fenêtre Rectangulaire
Excitation Pseudo Aléatoire

Avantages
z Pas de leakage
z Rapide
z Zoom
z Facteur de Crête acceptable
z Rapport Signal/Bruit acceptable
Inconvénients
z Pas d’approximation linéaire des systèmes non-linéaires
 Excitation Aléatoire contre Pseudo Aléatoire
Flanc de la résonance Pic de la résonance (leakage)
A=50 γ2(f1)=0.999 A=50 γ2(f2)=0.999
B(f1) B(f1)

Aléatoire

A(f1) A(f2)

Pic de résonance Anti-résonance

A=50 γ2(f3)=1 A=50 γ2(f4)=0.864
B(f1) B(f1)

Pseudo-Aléatoire

A(f3) A(f4)
Multisinusoïdal (Chirp)
Pour un sinus balayé répété chaque temps d’enregistrement, Tr

Temps
Tr
C’est un type spécial d’un signal pseudo aléatoire dans
lequel le facteur de crête est minimisé (< 2)

Il a les avantages et les inconvénients du signal pseudo

aléatoire normal mais avec un facteur de crête plus faible
Autres Avantages:
– Forme de spectre idéale: le spectre est constitué d’un
spectre d’amplitude plat et d’un spectre de la phase lisse

Applications:
– Essai des structures ayant un comportement non-linéaire
 Aléatoire Périodique
Un signal aléatoire et pseudo-périodique combiné donne un
signal d’excitation qui présente les caractéristiques suivantes :
–Pas de leakage dans l’analyse
–La meilleure approximation linéaire du système
A A A B B B C C C

Signal pseudo-aléatoire
changeant avec le
temps: T T T
réponse transitoire
réponse de l’état
permanent
Les données analysé sur le
temps: A B C
(réponse de l’état permanent)
Inconvénients:
z La durée de l’essai est plus long que celui de l’essai pseudo aléatoire ou
aléatoire
 Impulsion Périodique
Cas spécial du signal pseudo aléatoire
Impulsion Rectangulaire, de Hanning, ou Gaussienne sur un temps définie
par l’utilisateur ∆ t répété toute les ∆ T définie aussi par l'utilisateur

∆t ∆T Temps

Le spectre en ligne pour une

sin x
impulsion Rectangulaire avec
x
une enveloppe
1/∆ t Fréquence

z Leakage peut être diminué en utilisant une pondération rectangulaire temporelle

z La pondération transitoire et exponentielle temporelle peut être utilisé pour
augmenter le rapport Signal/Bruit
z Gating of reflections with transient time weighting
z Les effets de non-linéarité ne sont pas pondérées
z Le signal a un grand facteur de crête
Impulsion Périodique

Avantages Inconvénients

z Rapide z Pas d’approximation linéaire

z Pas de leakage d’un système non-linéaire
(seulement avec une pondération z Facteur de Crête élevée
rectangulaire)
z Un niveau de pics élevée peut
z Gating of reflections causer l’excitation des non-
(pondération transitoire temporelle) linéarités
z Le spectre de l’excitation suit z Pas de Zoom
dans la bande de base z Une pondération temporelle
spectrum follows frequency spéciale peut être requise
span in baseband pour augmenter le rapport
z Facile à implémenter Signal/Bruit . Ce qui peut
introduire un leakage
Éléments pour le choix de la Technique d’Excitation

z Pour l’étude des non-linéarités: Excitation par Sinus Balayé

z Pour un système légèrement Excitation Aléatoire

non-linéaire:

z Pour un système parfaitement Excitation Pseudo Aléatoire

linéaire:

z Pour les mesures sur Excitation par Impact

terrains:

z Pour les mesures sur terrain Excitation par Impact Aléatoire

hautes résolutions:
Mesures 3

f(t)
m x(t)
Les modèles SDOF et MDOF
c k
Différentes techniques d’Analyse
Modale

Exciter la structure

+ + +⋅⋅⋅+ Mesurer correctement les données

=

Le post-processing des données de

l’Analyse Modale
 Anatomie d’un Accéléromètre
Électronique hybride de
Masse sismique conditionnement
Vis élastique
Isolation électrique
Élément piézoélectrique

Connecteur électrique
Protection thermique et mécanique

Vis de fixation

Masse Sismique F=MA Schéma Simplifié

Cristal Piézoélectrique Cristal FET
FET Piézoélectrique
Signal

Structure à Tester
 Types d’Accéléromètres
PMC: piézo-électrique en mode charge, PMVEI: piézo-électrique en mode voltage à électronique intégré, PR: piézo-résistif,CV: capacité variable

PMC PMVEI PR CV
Poids relatif trés petit trés petit petit moyenne
Prix relatif moyen moyen petit moyenne
Design miniaturisé oui oui oui non
Résistance aux g élevée basse trés basse Trés grande (200000g) Trés grande (10000g)
Plage de fréquence usuelle trés élvée trés élvée moyenne basse
Limitation mécanique des surcharges non non oui parfois oui
Dérive du zéro aux chocs élevée mauvaise mauvaise bonne temps de recouvrement
Réponse en phase plate oui oui non-amorti oui/amorti non oui
Réponse DC non non oui oui
S ensibilité aux vibrations élvée élvée basse élvée
S ensibilité à l'environnement basse trés basse basse basse
Utilisation aux hautes températures oui non non non
Utilisation pour la mesure des chocs moyenne moyenne elvée basse
Impulsion longue ( test de crash) non non oui oui
Mesure de la distance non non non non
Générateur autonome oui non non non
Câbles longs ( base impédance de sortie) non oui oui oui
Température cryogénique oui non non non
Calibration par shunt non non oui oui
Types d’Accéléromètres
L’Environnement de la Mesure

Humidité Champs électromagnétiques

Fluctuations des températures Haut niveaux sonores

Hautes et basses températures Montage non conforme

Radiations Vibrations transversales

(rayons X, rayons γ)

Sortie en relation
Flexion de la base uniquement avec la
vibrations

Flexion de la base

Entrée
Vibration
Montage d’un Accéléromètre
Techniques de montage:
-Goujon Goujon
-Ciment
Ciment
-Base magnétique
Colle d’abeille max 40°C
-Base de Swivel
Niveau -Colle

dB
30

20 Base magnétique

10

200 500 1k 2k 5k 10k 20k 30k 50k

Sensibilité et Domaine de Fréquence
Chargement de l’Objet a Tester:
Masse Dynamique
-2
0.1 pC /m s
1
0.65 g ⇒ M > 7g < M
10
10 pC /m s -2
54 g ⇒ M > 600g
M

1000 pC/ms -2
470 g ⇒ M > 5kg
Test de Résonance de Montage des Accéléromètres
dB
m
Montage idéal
M
Hz k/m
fm = 2 π
dB
m m
Suspension libre k
M
Hz

m
dB
Câble? OK
M Montage non conforme Montage? OK
Accéléromètre? OK
Hz

dB Montage sur plaque légère

m
effet “ Chargement Masse”
M
Hz
 Analyseur Bi-Voies

Schéma de principe du fonctionnement d ’un analyseur de spectre =

Effectue des mesures temporelles et fréquentielles

a ( ω)
cmplx
Préamplificateur

a(t) Filtre
passe bas
Echantilloneur
Fenètre
Temporelle FFT
moyenne
FRFi
1/n Σ (FRFi )
Filtre Fenètre
b(t) passe bas
Echantilloneur
Temporelle FFT
Echantillonnage Traitement de Signal b ( ω)
cmplx

Données de Sortie et Affichage

Niveaux - Fréq. Echantillonnage Type et paramètres des

d'entrée - Nbre d'échantillons fenètres temporelles
- Méthode et paramètres de Moyennage
- Paramètres d'affichage
 Moyennage de Spectre Bi-Voies
Spectre Auto-
temps Auto-Spectre
corrélation
F
de Fourier F-1
a (t ) A( f )
GAA ( f ) Raa (τ )

Canal A Réponse en Spectre

de Fourier
Fréquence F-1 h(τ )
H1 ( f )
F-1 Auto-
Spectre
corrélation
Croisé Rab (τ )

GAB ( f ) Puissance

Canal B Cohérence de Sortie Cohérente

γ 2 ( f ) ⋅ GBB ( f )
F γ (f)
2

Spectre Auto-
temps de Fourier Auto-Spectre corrélation
b(t )
B( f ) GBB ( f )
F-1 Rbb (τ )

Enregistrement Analyse Moyennage Post-Processing

Planifier le Test Avant de Commencer !

1. Choisir l’Excitation Approprié

– Marteau, Vibrateur, ou OMA?

2. Mise au Point Correcte de l’Analyseur FFT

– Plage de Fréquence, Résolution, Pondération, Fenêtrage
– Attention: l’Analyseur FFT est un ANALYZEUR BLOCK !

3. Une Bonne Distribution des Points de Mesures

– S’assurer q’un nombre de points suffisants est mesurées pour voir
tous les modes d’intérêt
– Méfiez vous du ’spatial aliasing’

4. La Mise au Point Physique

– Le montage de l’accéléromètre est CRITIQUE!
– Uni-axial contre Triaxial
– S’assurer que l’orientation du DDL est correcte
– Monter l’objet à tester ...le montage peut affecter la mesure!
– Calibrer le système
Où est-ce que l'Excitation devrait être appliquée?

Réponse "i"
H ij = X i =
F j Excitation " j"

Mesure en 1 j {X} = [H]{F}

Point Entraînant 2
i
i=j ⎧ X 1 ⎫ ⎡H11 H12 ⎤ ⎧F1 ⎫
⎨ ⎬=⎢ ⎥ ⋅⎨ ⎬
⎩ X 2 ⎭ ⎣H21 H22 ⎦ ⎩F2 ⎭

Mesure de j
Transfert X 1= H11 ⋅ F1 + H12 ⋅ F2
i≠j X 2 = H21 ⋅ F1 + H22 ⋅ F2
i
Vérification de la Mesure dans un Point Entraînant
Im [Hij]
z Tous les pics doivent être dans

⎡ X (f) ⎤ ⎡ X&(f) ⎤ ⎡ X&

&(f) ⎤
Im⎢ ,Re et Im
⎣ F(f) ⎥⎦ ⎢⎣ F(f) ⎥⎦ ⎢ F(f) ⎥
⎣ ⎦

Mag [Hij]

z Chaque anti-résonance dans

Mag [Hij] doit être entouré par
une paire de résonance

Phase [Hij]
z Les variations de la Phase
doivent être autour de 180°
La mesure dans un Point Entraînant
La qualité de la mesure dans un PE est très importante , car
les résidus PE sont utilisé pour normer le Modèle Modal

Postulat au PE : Problèmes au PE :
– Les Résidus de tout les modes – Un couplage modal très élevée , avec
doivent être estimer précisément à tous les modes qui sont en phase
partir d’une seule mesure – Des effets résiduels très élevée venant
des modes du corps rigide

Ιm ⏐Aij⏐ Log Mag⏐Aij⏐ Mesuré

ω
=
Re ⏐Aij⏐ Sans les modes
du corps rigide
ω
ω
Tests pour la validité des données mesurées :
la Cohérence
2
GFX ( f )
Cohérence γ 2( f ) =
GFF ( f ) G XX ( f )

– Mesure combien l’énergie introduite dans le système cause la réponse

– la mesure est plus fiable lorsque la cohérence avoisine le ‘1’
– Moins de ‘0.75’ est la plage de mauvaise cohérence
z Son interprétation est que pour chaque fréquence , elle montre le
degré de la relation de linéarité entre l’entrée et la sortie .
 Raisons pour une Faible Cohérence

Mesures Difficiles :
z Bruit dans le signal de sortie mesurée
z Bruit dans le signal d’entrée mesurée
z Les autres entrées ne sont pas en corrélation avec le signal d’entrée mesuré

Mauvaises Mesures:
z Leakage
z Systèmes variants avec le temps
z Les non-linéarités du système
z Instabilité des DDL
z Le temps de Propagation n’est pas compensé
Comment augmenter la Cohérence?
Une mauvaise cohérence peut indiquer la présence de plusieurs problèmes
• Si la cohérence est faible dans les pics modaux, le leakage probablement
affecte la mesure ⇒ Augmenter le temps d’échantillonnage ou changer la
nature du signal d’excitation
• Le capteur est dans un noeud ⇒ Changer la place du point de mesure
• Des cohérences faibles aux petites fréquences, problème fréquent pour les
capteurs piézoélectriques ⇒ Opter pour un vibrométre laser pour la mesure
de la réponse
• La présence de non-linéarités dans le système ⇒ Il faut les identifié et les
éliminer
• Bruit dans la mesure ⇒ Vérifier le signal temporel et s’assurer que
l’acquisition se fait dans un niveau de voltage convenable
 Tests pour la validité des données mesurées:
la Linéarité
X1 = H·F1 X1+X2 = H·(F1 + F2)
Linéarité ⇒
X2 = H·F2 a·X1 = H·(a· F1)

H(ω)
– Plus qu’on entre de force dans le
système autant de réponse égale
F(ω) X(ω) à cette force en sortira
– Puisque la FRF est un rapport ,
l’amplitude de la FRF doit être la
même quelque soit la force
d’entrée

La Linéarité peut être détecté en appliquant différents niveaux

d’amplitude et en observant la proportionnalité entre l’entrée et
la sortie. Ceci est fait usuellement en utilisant une excitation
sinusoïdale
Tests pour la validité des données mesurées:
la Réciprocité

H ij( f )= H ji( f )
j

La Réciprocité est le phénomène qui fait que la réponse mesuré dans le

DDL j et excité dans le DDL i donne le même résultat comme une
réponse mesurée dans le DDL i et excité dans DDL j
 Astuces pour des Obtenir des
Résultats de Meilleure Qualité
z Vérifier l’intégrité de la chaîne de mesure avent de commencer le test :
–La calibration des capteurs
–La calibration du rapport de masse
z Vérifier la convenance des capteurs d’entrées et de sorties :
–Plages Opérationnelles (fréquence, plage dynamique, réponse en phase)
–La masse additionnelle des accéléromètres
–Le montage des accéléromètre
–La sensibilité aux facteurs de l’environnement
–La Stabilité
z Vérifier la convenance des pré-réglage de l’essai :
–La postions des capteurs et leurs alignement
–Pré-test: rattling, CAF, modes du corps rigide, rapport Signal/Bruit, approximation linéaire,
signal d’excitation, racines multiples, réciprocité de Maxwell, mesure de la force,
connexion excitateur-capteur d’entrée-tige-structure

Des Mesure de Qualité pour la FRF sont la Base de

L’Analyse Modale Expérimentale!
Lissage des Courbes et Validation

f(t) Les modèles SDOF et MDOF

m x(t)

Différentes techniques d’Analyse

c k
Modale

Exciter la structure

Mesurer correctement les données

= + + +⋅⋅⋅+
Le post-processing des données de
l’Analyse Modale
De l’Essai vers l’Analyse
H( f )
FRF
Mesurée

Fréquence
Lissage de Courbe
(Pattern Recognition)

H( f )

Analyse Modale

Fréquence
De l’Essai vers l’Analyse

H( f )

Analyse Modale

Fréquence

f(t) f(t) f(t)

m x(t) m x(t) m x(t)
Modèles SDOF
c k c k c k
 Caractérisations des Modes

Tout Les Modes Peuvent Être Caractérisé Par:

1. Fréquence de Résonance
2. Amortissement Modal
3. Déformée Modale
 Analyse Modale – Étape par Étape
1. Inspection Visuelle des Données
–Cherchez les modes évidents dans la FRF
–Inspecter TOUTES LES FRFs…parfois des modes sont visibles dans une seule
FRF et pas dans toutes (noeuds )
–Utilisez partie Imaginaire et cohésion pour vérification
–Sommer l’amplitude de toutes les mesures pour les indices
2. Choisir le Lisseur de Courbe Adéquat
–Modes légèrement couplés: les techniques SDOF
–Modes fortement couplés: les techniques MDOF
–Mesures stable: les techniques globale
–Mesures instables: les techniques locales
–Mesures MIMO : les techniques Poly-référence
3. Analyse
–Utiliser plus qu’un lisseur de courbe pour voir s’il se concorde
–Faire attention aux calculs de Résidus
–Est ce que les déformées modales ont un sens ?
 Analyse Modale – Inspection des Données
1. Inspection Visuelle des Données
–Cherchez les modes évidents dans la FRF
–Inspecter TOUTES LES FRFs…parfois des modes sont visibles dans une seule
FRF et pas dans toutes (noeuds )
–Utilisez partie Imaginaire et cohésion pour vérification
–Sommer l’amplitude de toutes les mesures pour les indices
2. Choisir le Lisseur de Courbe Adéquat
–Modes légèrement couplés: les techniques SDOF
–Modes fortement couplés: les techniques MDOF
–Mesures stable: les techniques globale
–Mesures instables: les techniques locales
–Mesures MIMO : les techniques Poly-référence
3. Analyse
–Utiliser plus qu’un lisseur de courbe pour voir s’il se concorde
–Faire attention aux calculs de Résidus
–Est ce que les déformées modales ont un sens ?
 Analyse Modale – Lissage de Courbe
1. Inspection Visuelle des Données
–Cherchez les modes évidents dans la FRF
–Inspecter TOUTES LES FRFs…parfois des modes sont visibles dans une seule
FRF et pas dans toutes (noeuds )
–Utilisez partie Imaginaire et cohésion pour vérification
–Sommer l’amplitude de toutes les mesures pour les indices
2. Choisir le Lisseur de Courbe Adéquat
–Modes légèrement couplés: les techniques SDOF
–Modes fortement couplés: les techniques MDOF
–Mesures stable: les techniques globale
–Mesures instables: les techniques locales
–Mesures MIMO : les techniques Poly-référence
3. Analyse
–Utiliser plus qu’un lisseur de courbe pour voir s’il se concorde
–Faire attention aux calculs de Résidus
–Est ce que les déformées modales ont un sens ?
 Comment Marche le Lissage de Courbe?
z Le Lissage de Courbe est le processus avec lequel on
estime les Paramètres Modaux à partir des mesures

2σ
z Chercher la fréquence de
|H|
résonance
– Une fréquence qui avec
R/σ une petite excitation cause
une grande réponse

z Chercher l’amortissement
ω – Quel est le Q du pic?
ωd
z Chercher le résidu
0 – Essentiellement ‘la surface
Phase sous la courbe’
-180
Fréquence
Comment Marche le Lissage de Courbe?

2σ R
H ≅
|H| σ − ( jω d − jω )

avec
R/σ σ = ω 0ζ
ωd = ω0 1−ζ 2

ω
ωd

Le lissage des courbes est une méthode de réduction des données

- Les informations peuvent être perdus
- L’information perdue est un bruit utile!
Le lissage des Courbes

• Transformation
• Lissage par moindres carrées
• Les méthodes de lissage des courbes
•Polynomiale
•Exponentielle complexe

• Utilisation des termes résiduels

Transformation
Temps
(t) L ⎡⎣ x ( t ) ⎤⎦ = X ( s )
F ⎡⎣ x ( t ) ⎤⎦ = X (ω )

“Choisir le domaine des

variables indépendantes F −1 ⎡⎣ x ( t ) ⎤⎦ L−1 ⎡⎣ x ( t ) ⎤⎦
pour simplifier les calculs
et l’intrepretation” Synthèse
Fréquence Laplace
(ω) (s)

Lissage

• Les transformées font la conversion d’un domaine vers un

autre sans perte ou gain d’information
• Les transformées sont généralement réversibles
• Le processus de lissage des courbes est irréversible et
l’information peut être perdue- l’information perdue est le bruit
Lissage des courbes

• Transformation
• Lissage par moindres carrées
• Les méthodes de lissage des courbes
•Polynomiale
•Exponentielle complexe

• Utilisation des termes résiduels

 Lisseur de courbes SDOF

|H|
H4
H1 H2 1
X ( jω ) m
H (s) jω = =
H3
F ( jω ) ( jω ) 2 + jω 2ζ + ω 0

ω
1
• Équations avec trois inconnues ζ ω0
m
• La solution est trouvée par les trois observations H1 H2 H3
• Prendre d’autres points sur la courbe va engendré un problème surdéterminé
Méthode des Moindres Carrées
Modéle Mathématique:
Yk = bX k + N k
Y
avec N k est l'erreur aléatoire
Erreur Aléatoire Xk,Yk (mesures individuelles)

N k = Yk − bX k

∑ Nk = ∑ (Yk − bX k )
2 2

X
∂∑ Nk
2

= −2∑ (Yk − bX k )X k
∂b
∂∑ Nk
2

= 0 donne l'estimation bˆ de b qui minimise ∑ N k

2

∂b

bˆ =
∑ Yk X k
∑ Xk2
La MMC est une méthode qui élimine le bruit aléatoire
 Processus Paramétrique
Y
Les paramètres modaux sont estimés
à partir d’un modèle paramétrique
y = bx lissé au données du signal traité

Y X Données expérimentales

y = bx + c Modèle Paramétrique

Traitement du Signal

Y X
Estimation des Paramètres

y = ax 2 + bx + c
Paramètres Modaux

X
Lissage des Courbes

• Transformation
• Lissage par moindres carrées
• Les méthodes de lissage des courbes
•Polynomiale
•Exponentielle complexe

• Utilisation des termes résiduels

Techniques de Lissage des Courbes

Système à un degré de liberté (SDOF)

•Polynôme SDOF
•Lissage de Cercle

Système à plusieurs degrés de liberté (SDOF)

•Polynôme SDOF
•Exponentielle Complexe
•Paramètre direct
•Polyréférence
• ………………..
 Méthodes de Lissage des Courbes

Basée sur la mesure d’1 ligne ou 1 colonne de la matrice de réponse en

fréquence
Méthodes avec Mode Unique
• Rapide
• Simple à utiliser
• Les erreurs on lieu lorsque les modes sont fortement couplées

Méthodes avec Modes Multiples

• Traite le Couplage Modal
• Demande un opérateur plus qualifié
 Lisseurs de Courbes Évoluées

• Lisseur de Courbe Global

¾ la fréquence et l’amortissement sont estimé globalement de la totalité

des mesures en utilisant les méthodes des moindres carrées et le modèle
polynomial lissé

¾les résidus sont calculées à partir de la totalité des mesures en utilisant

les valeurs globales de l’amortissement

• Lisseur de courbe polyréférence

¾utilisé pour l’identification des racines répétés

¾deux ou plus de colonnes ou de lignes de la matrice FRF sont mesurées

¾un lisseur de courbe global utilisant le modèle de lissage polynomial

Lisseur de Courbe Polynomial SDOF

R R*
H (ω ) = + + A + A ω + A ω 2

jω jω * 0 1 2
p p
R = R1 + jR2 est le résdiu
p = σ + jωd
σ est le taux de décroissance
ωd est la fréqunce naturelle amortie
A0 , A1 , A2 sont les fonctions de constantes résiduelles
( effet des modes hors bande)
Lisseur de Courbe Polynomial

1 1
H (s) = m H (s) = m
SDOF ( s − p )( s − p * ) s 2 + js 2ζ + ω 0

Forme en Fraction Partielle

m
[ R ]r ⎡⎣ R * ⎤⎦
H (ω ) = ∑ + r

r =1 jω − pr jω − pr *

Forme en Fraction Rationnelle

am s m + am −1 s m −1 + .... + a0
H (ω ) =
bn s n + bn −1 s n −1 + .... + b0
Lisseur de courbe polynomial

am s m + am −1 s m −1 + .... + a0
|H| H (ω ) =
bn s n + bn −1 s n −1 + .... + b0

e1 ei Inconnues: am ........a0 , bn .........b0

en ei = FRFanalytique − FRFmesurée
Erreur:

Problème:
Minimiser le vecteur erreur
N
{ J } = ∑ ei*e j
i =1
 Lisseur de courbe polynomial

+ + + …….. =

Travaille mieux dans:

• Les bandes limités en fréquence
• Modes limités en nombre

La fonction de lissage contient des termes extra polynomiaux qui sont utilisé lors le lissage des
courbes pour comptabiliser les effets des modes hors bande.
Exponentielle Complexe

+ + +…=

Méthode:
• MDOF h ( t ) = 2 R e − σ t s i n (ω d t )
• Domaine temporel
Avantages:
• Prend en compte les modes multiples
• Précision du domaine temporel
• Dans Me’scope , son implémentation est du type Globale
Exponentielle Complexe

+ + +…=

Travaille mieux dans:

• les bandes larges en fréquence
• un ou plusieurs Modes extra (calculées) en plus que les pics modaux
sont présent dans la bande du curseur

Si le bouton Count Peak est utilisé pour déterminer le nombre des

modes, augmenter le nombre dans la fenêtre Modes par un ou
plusieurs modes extra avant l’utilisation de cette méthode
 Exponentielle Complexe

+ + +…=

Diagramme de Stabilisation
Nombre de Modes
Modes
7 Stables
6
5
4

Fréquence
Les modes physiques sont des modes qui se répètent pour de multiples
ordres de modèles
Les modes stables sont des modes physiques qui se répètent
consécutivement dans deux modèles ou plus et qui satisfassent à un critère
définit pat l’utilisateur ( MAC, Phase Scatter…etc)
Lissage des Courbes

• Transformation
• Lissage par moindres carrées
• Les méthodes de lissage des courbes
•Polynomiale
•Exponentielle complexe

• Utilisation des termes résiduels

Termes Résiduels
|A|

ω
|A|

|A|

ω
 Les Résidus sont Directement
Liées à la Déformée Modale !
* z Les Résidus exprime la
Rijr Rijr rigidité d’un mode pour
Hij (ω) = ∑ Hijr = ∑ +
r r jω − p r jω − pr* chaque FRF mesurée

z Par conséquent ils sont

Second
fonction de la déformée
Amplitude Premier Mode modale en chaque point
anc
e Mode mesurée
t
Dis
Troisième
Mode
Fréq Poutre
uenc
e Accélération
Force
Force
Force
Force
Force
Force
Force
Force
Force
Force
Force
Force
SDOF contre MDOF pour le Lissage des Courbes
z Utiliser les méthodes SDOF pour les modes FAIBLEMENT COUPLEES

z Utiliser les méthodes MDOF pour les modes FORTEMENT COUPLEES

z On peut combiner les deux techniques SDOF et MDOF !

SDOF MDOF
(couplage faible) (couplage fort)
|H|

ω
Lissage de Courbes Local contre Global
z Local : les résonances, amortissement, et les résidus sont
calculée pour chaque FRF en premier lieu …et après ils
sont combinées pour le lissage des courbes

z Global : les résonances, l’amortissement, et les résidus sont

calculées a travers toutes les FRFs

|H|

ω
 Résumé sur le Lissage des Courbes
Système à un Degré De Liberté (SDOF)
z Structure simple
z Peu de modes
z Des modes largement espacés
Système à Plusieurs Degrés De Libertés (MDOF)
z Modes multiples
z Beaucoup de modes
z Des modes proches entre eux
[(m/s²)/N] (Magnitude))
Frequency Response H1(Response,Excitation) - Input (Magnitude
Working : Input : Input : FFT Analyzer

10
Locale
z Basé sur un seul DDL
100m

Global
Basé sur un plusieurs DDLs
0 200 400 600 800 1k 1,2k 1,4k 1,6k
[Hz]
z

Monoréférence
z Basé sur une seule colonne ou une seule ligne de la matrice FRF
Polyréférence
z Basé sur des lignes ou des colonnes multiples de matrice FRF
z Les structures symétriques montre plusieurs modes à la même fréquence
z Un seul pique ne veut pas dire nécessairement un seul mode
 Analyse Modale – Analyse des Résultas
1. Inspection Visuelle des Données
–Cherchez les modes évidents dans la FRF
–Inspecter TOUTES LES FRFs…parfois des modes sont visibles dans une seule
FRF et pas dans toutes (noeuds )
–Utilisez partie Imaginaire et cohésion pour vérification
–Sommer l’amplitude de toutes les mesures pour les indices
2. Choisir le Lisseur de Courbe Adéquat
–Modes légèrement couplés: les techniques SDOF
–Modes fortement couplés: les techniques MDOF
–Mesures stable: les techniques globale
–Mesures instables: les techniques locales
–Mesures MIMO : les techniques Poly-référence
3. Analyse
–Utiliser plus qu’un lisseur de courbe pour voir s’il se concorde
–Faire attention aux calculs de Résidus
–Est ce que les déformées modales ont un sens ?
Quel Lisseur de Courbe doit on utilisé?
Fonction de Réponse en Fréquence
Hij (ω) pour couplage faible

Réel Imaginaire

Nyquist Log Magnitude

Quel Lisseur de Courbe doit on utilisé?
Fonction de Réponse en Fréquence
Hij (ω) pour couplage fort

Réel Imaginaire

Nyquist Log Magnitude

 La Synthèse du Modèle Modal
Théorique Expérimental

Hˆ ij

continu distribué structure

ou
Calculer: [ m ],[ k ], ( [c ] ) Mesurer:
{H1}Ligne ou Colonne
Résoudre:
Domaine-t: [m ]{x&&} + [c ]{x&} + [k ]{x} = {f } Estimer:
Domaine-s: ⎡[ m] s 2 + [c ] s + [ k ]⎤ { X (s )} = {f }
mat diag [ p] et [Φ]
⎣ ⎦ m
ar
[B ]; matrice carctéristique: det [B ] = 0 ⇒ [ λ] [H ] = ∑ jω − p {Φ}r + {Φ}r
T

r =1
⎡⎣B ( λ ) ⎤⎦ { x} = 0 ⇒ {ψ}λ ⇒ {Ψ} r

Modéles de la Réponse: [B]{x} = {F} Relation { X} = [H]{F}

[B] = [H ]
−1
Au domaine-ω :
 Le Modèle de la Réponse
{ X (ω )} = [ H (ω )]{F (ω )}
Avec:
m
[ R ]r ⎡⎣ R ∗ ⎤⎦
H (ω ) = ∑ + r

r =1 jω − p r jω − p r ∗
m Rijr Rijr ∗
H ij = ∑ +
r =1 jω − p r jω − p r ∗
On démontre que : [ R ]r = a r {φ }r {φ }r ⇒ Rij = a r φ ir φ jr
T

⎡ φ1φ1 φ1φ 2 ⋅ ⋅ φ1φ j ⋅ ⋅ φ1φ m ⎤

⎢φ φ φ 2φ j ⋅ ⎥⎥
On voit que chaque
⎢ 2 1 colonne/ligne de
⎢ ⋅ ⋅ ⋅ ⎥ [R] contient la
[ R ]r =⎢ ⎥ ⋅ ar déformé modale
⎢ φ iφ 1 φ iφ j φ iφ m ⎥ multiplié par une
élément
⎢ ⋅ ⋅ ⋅ ⎥ appartenant a cette
⎢ ⎥
⎢⎣φ nφ1 ⋅ ⋅ ⋅ φ nφ j ⋅ ⋅ φ nφ m ⎥⎦ n = m colonne/ligne
 Estimation de la Déformé Modale

Ligne (i ) = { Rij } = {φ iφ1 . φ iφ m }r ⋅ a r

T

r
φ iφ 2 . φ iφ i

En général:
⎧ φ1φ j ⎫
⎪φ φ ⎪ Rijr = a rφ irφ jr
⎪ 2 j⎪
⎪⎪ . ⎪⎪ Dans le Point Entraînant:
Colonne ( j ) = { Rij } =⎨ ⎬ ⋅ ar
⎪φ jφ j ⎪ Rijr = a rφ jr 2
r

⎪ . ⎪
⎪ ⎪ Rijr
⎩⎪φ nφ j ⎭⎪ φ jr =
ar

Le vecteur de la déformé modale complet est calculé à partir du Point

Entraînant
 La Normalisation (la mise à l’échelle)
• La Normalisation de la Déformée Modale • La Normalisation pat unité de masse
Pour un systéme SDOF: Modale :
1 1
M =1 a=
R= m [1] donne
2 jω d
Générallement: 2 jω d
Rijr = a r φ iφ j MDOF:
2 jω dr
[ Φ ] [ m ][ Φ ] = mat diag [ M r ] {φ }r = { R }r
T

Rijr
SDOF: R = aφ 2 [2]
φ mφ = M [3]
Remplaçons [2] et [3] dans [1]:
φ 2
aφ 2 =
2 jω d M
1
a=
M 2 jω d
Le Processus de Normalisation

• Rijr est déterminé a partir des mesures du point entraînant

• Dans le cas d’une Masse Modale unitaire Mr=1
− ar est calculé a partir de
1 1
ar = =
2 jω d M r 2 jω d
− Les déformées modales normalisées au Point Entraînant sont
déterminé a partir de Rijr et ar

φ jjr = R jjr 2 jω dr
− Le vecteur de déformée modale complet est calculé à partir des
mesures des fonctions de réponse
2 jω dr
{φ }r = { R }r
Rijr
La synthèse de Hij

Les Paramètres Modaux de base

m a r φ ir φ jr a r ∗φ ir ∗φ jr ∗
H ij (ω ) = ∑ +
r =1 jω − p r jω − p r ∗
( p r = −σ r + jω dr )
N.B: quelques logiciels d’analyse modale utilisent les résidus “réels”
1
R= m ⇒ masse modale unitaire
ωd
1
a=
ωd
m a r φ irφ jr a r ∗φ ir ∗φ jr ∗
H ij (ω ) = ∑ +
r =1 2 j ( jω − p r ) 2 j ( jω − p r ∗ )
Les Mesures Point Entraînant (PE)
La qualité de la mesure dans un PE est très importante , car
les résidus PE sont utilisé pour normer le Modèle Modal

Postulat au PE : Problèmes au PE :
– Les Résidus de tout les modes – Un couplage modal très élevée , avec
doivent être estimer précisément à tous les modes qui sont en phase
partir d’une seule mesure – Des effets résiduels très élevée venant
des modes du corps rigide

Ιm ⏐Aij⏐ Log Mag⏐Aij⏐ Mesuré

ω
=
Re ⏐Aij⏐ Sans les modes
du corps rigide
ω
ω
 Les Modèles Tronqués (modes)
Modèle réel
Log ⏐H⏐

Modèle tronqué

Modes dans le Modèle Rijr Rijr ∗

Modèle tronqué H ij = ∑
r jω − p r ∗
+
jω − p r ∗
1 1
Modèle “amélioré” H ij ' = + H +
m resω 2
ij
k res
Termes Résiduels
Conclusion
• Les termes résiduels compensent le manque de la banse en mode, mais ne peut être utilisé
pour la synthèse des FRFs non mesurées
• Le modèle modal ne peut être amélioré en ajoutant les termes résiduels
• Les termes résiduels ne sont utilisé que lors du processus de lissage des courbes
La validation d’un Modèle Modal

Outils de Validation
z MAC (Modal Assurance Criteria)
z MCF (Modal Confidence Factor)
z Dispersion de la phase (Phase Scatter)
z Diagramme de stabilité
Modal Assurance Criterion, MAC
z Synthesized FRF
z Modal Assurance Criteria MAC
T 2
ψi ψ j
MAC(i, j) =
(ψiT ψi ) ⋅ (ψ Tj ψ j )

[ψ ]T [ M ][ψ ] = m at diag { m1, m 2 ,K m n }

[ψ ]T [ K ][ψ ] = m at diag { k 1, k 2 ,K k n }
 Au-Delà de l’ Analyse Modale
Modèle L’Essai Modal
Prédire le modèle de l’essai
Élément Finis
Les données sur
l’amortissement

Comparaison et Lissage
Validation du Validation Validation
Modèle Modal des courbes

Mise à jour

Modèle Final

Simulation type
Résolution des problèmes
“Il Arrive Quoi Si?”
Au-Delà de l’ Analyse Modale

Analyse Modale
Expérimentale

Modèle Dynamique
basée sur les
Paramètres Modaux

F
Modification Simulation
Structurale de la Réponse

Modification Matérielle Simuler la Réponse Réelle

Spécification de la Résonance
 Conclusion

z Toutes les Structures Physiques peuvent être caractérisé par un modèle

SDOF simple

z Les Fonctions de Réponse en Fréquence sont la meilleure moyen a utiliser

pour mesurer les résonances

z Ils existe trois techniques d’analyse modale disponible aujourd'hui : Essai

par Marteau, Essai par Vibrateur, et l’analyse Opérationnelle

z La planification et réglage correct avant de commencer le test peut

diminuer le coût et les efforts…et assure la précision en minimisant les
mauvais résultats

z Il existe beaucoup de techniques de lissage de courbes disponible, il faut

utiliser la technique la plus adapté au problème traité
 Littérature complémentaire

z Frequency Analysis par [Link]

(Brüel & Kjær Theory and Application Handbook BT 0007-11)

z Modal Analysis of Large Structures - Multiple Exciter Systems par K. Zaveri

(Brüel & Kjær Theory and Application Handbook BT 0001-12)

z Modal Testing: Theory and Practice par D.J. Ewins

(Brüel & Kjær Theory and Application Handbook BT 0015-12)

z Dual Channel FFT Analysis par H. Herlufsen

(Brüel & Kjær Technical Review No. 1 & 2, 1984)

z Structural Testing, Part 1: Mechanical Mobility Measurement par O. Døssing

(Brüel & Kjær Theory and Application Booklet BR 0458-12)

z Structural Testing, Part 2: Modal Analysis and Simulation par O. Døssing

(Brüel & Kjær Theory and Application Booklet BR 0507-11)
Appendice: Paramètres d’Amortissement

2σ
Bande à 3 dB ∆f3 dB = , ∆ω3 dB = 2σ
2π

1 ∆f3 db ∆ω3 dB
Facteur de Perte η= = =
Q f0 ω0

η ∆f3 dB ∆ω3 dB
Pourcentage ζ= = =
2 2f0 2ω0
d’Amortissement
∆ω3 dB
Constante σ = ζ ω0 = π ∆f3 dB =
de Décroissance 2

f0 ω0
Facteur de Qualité Q= =
∆f3 dB ∆ω3 dB
 Appendice: Paramètres d’Amortissement

h(t ) = 2⋅ R ⋅ e−σt ⋅ sin( ωdt ) , où la constante de Décroissance est donnée par e-σt
∇
~
L’Enveloppe est donnée par l’amplitude analytique h(t):h( t ) = h2 ( t ) + h 2 ( t )
e − σt
h(t)
1
σ=
constante de décroissance τ
1 Temps
τ=
Constante de Temps σ La Constante de Temps, τ, est
σ 1 déterminé par le temps que va
Pourcentage d’Amortissement ζ= =
ω0 2πf0 τ prendre l’amplitude pour
1 décroître par un facteur de e =
Facteur de perte η = 2⋅ζ =
πf0 τ 2,72…
Qualité 1 ou
Q= = πf0 τ
η 10 log (e2) = 8.7 dB