Modal Analysis
Modal Analysis
Correlation
Structural Dynamics Agenda Topics
structure-
structure-borne air-
air-borne
structure-
structure-borne
Validate
Cost of Change
Engineer
Natural frequency is the frequency at which a system naturally vibrates once it has
been forced into motion
f(t)
x(t)
k c
ground
k
n natural frequency (rad/sec)
Unrestricted © Siemens AG 2019 m
Page 10 Siemens PLM Software
Natural Frequency
Resonance is the buildup of large amplitude that occurs when a structure is excited
at its natural frequency
Amplitude
ω f = 0.4 ω f = 1.01 ω f =1.6
3 Single Degree of Freedom Systems with ω n = 1.01 Frequency
same mass, stiffness and damping
c
d n 1 2
FORCE
!
Unrestricted © Siemens AG 2019
Page 16 Siemens PLM Software
Frequency Response Functions
FORCE
-50.00
N2/Hz
dB
-100.00
0.00 Hz 1024.00
Amplitude
1.0
g/N
Damping 0.7
0.5
0.3
0.0
618
Mode Shape Phase
°
1.6
1.4
1.3
1.2
1.1
1.0
0.9
Imag
0.8
g/N
0.7
0.6
0.5
0.4
0.3
0.2
0.1
-0.1
1
Real
g/N
-700e-3
Unrestricted © Siemens AG 2019 10 50 100 150 200 250 300 350 400 450 500 550 600 651
Hz
Page 19 Siemens PLM Software
Quality Factor
2776.18 2722.33
2800.0
2400.0 2228.38
2355.74
1000.0
400.0
0.0
618
Phase
°
1st mode
2nd mode
3rd mode
Unrestricted © Siemens AG 2019
Page 23 Siemens PLM Software
Mode Shape:
3rd mode
1.6
1.0
Imag
g/N
0.5
-0.1
1
g/N
-700e-3
Unrestricted © Siemens AG 2019 10 100 200 300 400 500 651
Page 28 Siemens PLM Software
Hz
FRFs determine mode shapes
Modal parameters
• Frequency ▪ Troubleshooting
• Damping ▪ Simulation and prediction
▪ Optimization
• Mode Shape ▪ Diagnostics and health monitoring
f(t)
x(t)
m
k c
ground
14
10
8
0
0 25 50 100 150 200 250 300 350 400 450 500
Hz
Measure the Frequency Response Functions
Response
Frequency
Damping
Mode Shapes
Input Input
Unrestricted © Siemens AG 2019
Page 32 Siemens PLM Software
Linearity of the FRF
1.00 0.10
)
g/N
Log
(
Log
N2
FRF DRV:1:+X/FOR:1:+X
FRF DRV:1:+X/FOR:1:+X (1)
FRF DRV:1:+X/FOR:1:+X (2)
F AutoPow er FOR:1:+X
F AutoPow er FOR:1:+X
Hz
F AutoPow er FOR:1:+X
10.0e-6 0.00 Linear 100.00
180.00 0.00 Linear 100.00
Hz
Phase
°
1.00e-6 -180.00
0.00 Hz 100.00 0.00 Hz 100.00
Excitation
• Laboratory
(shakers, hammer, force cell, …)
• Operational excitations
(road simulation, flight simulation, wind
excitation, …)
• Unusual excitations
(loudspeaker, gun shot, explosion, …)
Response
• (Accelerometers, Laser,…)
Measurement system
• Data acquisition system: SCADAS
(minimum 2 channels)
• Processing on PC (Test.Lab)
Input
Good for wide range of structures
Limited frequency range
Response
FRF
Which tip
should I
use?
Amplitude
/Hz
Log
g/N
dB
/
N2
Metal Tip
Amplitude
/Hz
Log
g/N
dB
/
2
N
Nylon Tip
Amplitude
/Hz
Log
g/N
dB
/
N2
Rubber Tip
Any Difference?
F1x
F1y
F1z
F2x
F2y
F2z
…
F9x
F9y
F9z
Unrestricted © Siemens AG 2019
Page 43 Siemens PLM Software
Modal Test
Nine Point Modal x Three Directions
• Force In the ‘Z’ Direction. (One Direction)
• Triaxial Accelerometer
F1x
F1y
F1z
F2x
F2y One complete Row or one
F2z complete column required
to calculate mode shape
…
F9x
F9y
F9z
Unrestricted © Siemens AG 2019
Page 44 Siemens PLM Software
Modal Test
Nine Point Modal x Three Directions
• Force In the ‘Z’ Direction. (One Direction)
• Triaxial Accelerometer
F1x
F1y
F1z A1x/F1z A1x/F1z A1x/F1z A1x/F1z A1x/F1z A1x/F1z … A1x/F1z A1x/F1z A1x/F1z
F2x
F2y
F2z
…
F9x
F9y
F9z
Unrestricted © Siemens AG 2019
Page 45 Siemens PLM Software
Modal Test
Nine Point Modal x Three Directions
• Force In the ‘Z’ Direction. (One Direction)
• Triaxial Accelerometer
F1x
F1y
F1z A1x/F1z A1x/F1z A1x/F1z A1x/F1z A1x/F1z A1x/F1z … A1x/F1z A1x/F1z A1x/F1z
F2x
F2y
F2z
Complete row possible to
calculate mode shape ☺
…
F9x
F9y
F9z
Unrestricted © Siemens AG 2019
Page 46 Siemens PLM Software
Modal Test
Nine Point Modal x Three Directions
• Force In the ‘Z’ Direction. (One Direction)
• Triaxial Accelerometer
F1x
F1y
F1z A1x/F1z A1y/F1z A1z/F1z
F2x
F2y
F2z A1x/F2z A1y/F2z A1z/F2z
…
F9x
F9y
F9z A1x/F9z A1y/F9z A1z/F9z
F1x
F1y
Incomplete!
F1z A1x/F1z A1y/F1z A1z/F1z
F2x
F2y
F2z A1x/F2z A1y/F2z A1z/F2z
…
F9x
F9y
F9z A1x/F9z A1y/F9z A1z/F9z
Time Frequency
Time Consuming to setup
Control frequency range
Control Force Amplitude
Input
Better for larger structures
Typically fixed excitation point, multiple
response points - measured in batches
Response
FRF
1.00
F Coherence DRV:1:+X
F Coherence DRV:2:+X
F Coherence ENG:1:+Y
F Coherence FUSL:5:+X
0.00
0.00 Hz 100.00
-1
-2
Joseph Fourier
-3 (º1768 - †1830)
-4 Théorie analytique de la chaleur (1822)
2.5
• Fourier’s law of heat conduction
2
u 2u 2u
1.5
k 2 2
0.5
t x y
-0.5
-1
-1.5
-2
-2.5
Amplitude
Amplitude
Amplitude
Frequency
Time
(Hz)
(seconds)
When the spectral content of your signal does not correspond to an available spectral line.
f 1 Hz
5 V Sine Wave - 3 Hz 5 V Sine Wave - 2.5 Hz
5V 5V
1 2 3 4 5 6 1 2 3 4 5 6
Hz Hz
63 Unrestricted © Siemens AG 2019
Page 63 Siemens PLM Software
Periodic Signals
T T = N t
FFT stiches time over infinity Are these signals the same?
YES!
T T = N t
!
Unrestricted © Siemens AG 2019
Page 66 Siemens PLM Software
DSP Error: Leakage
Smearing of
spectral content across
ALL FREQUENCIES
T T = N t
Must be removed
Original signal
Window function
1
x
=
Windowed signal
AKA No Window
Freq. domain
1 2 1 2 3
3
Mode 1
Response DOF
1
Mode 2
3
Mode 3
Driving point FRFs
Bad
excitation
point
FRF Point:1
FRF Point:1
FRF Point:2
FRF Point:1
FRF Point:2
FRF Point:3
FRF Point:1
FRF Point:2
FRF Point:3
FRF Point:4
FRF Point:1
FRF Point:2
FRF Point:3
FRF Point:4
FRF Point:5
FRF Point:1
FRF Point:2
Which is the
appropriate
FRF Point:3
excitation point?
FRF Point:4
FRF Point:5
FRF Point:1
FRF Point:2
FRF Point:3
FRF Point:4
FRF Point:5
Pneumatic suspension
Elastic cords
No influence from
surrounding
modes
Frequency
Large influence
from surrounding
modes
Frequency
H ()
n
vi liT
v l
*
i i
H
i 1 j i j *
i
(m/s2)/N
• How many modes to curve-fit?
Log
Solutions
• Stabilization diagram
• Mode indicator functions
f1 = 50 Hz
d1 = 20 %
Amplitude
25 50 75 Frequency
f1 = 25 Hz f2 = 75 Hz
d1 = 10 % d2 = 10 %
Amplitude
25 50 75 Frequency
f1 = 25 Hz f2 = 50 Hz f3 = 75 Hz
d1 = 5 % d2 = 5 % d3 = 5 %
Amplitude
25 50 75 Frequency
Stability
o : new pole
Amplitude
Model Order
f : frequency
d : damping
s : all
0
Frequency
Mode Indicator Function (MIF) helps identify the modes for a system where multiple
reference FRFs were measured
▪ commonly used to detect the presence of repeated roots
Double “dip”
indicates two
modes at same
frequency
Modal Assurance Criterion (MAC) describes how similar the shapes are for a given mode
pair using a scale of 0 to 1 (e.g. 0% to 100%)
.5
6 Points
15 Points
.5
.5
• 0 – 60 Hz band
• PolyMAX: 90/105 modes
• LSCE: 70/105 modes
1.00 1.00
)
(m/s2)/N
Log
(
Amplitude
/
Hz
10.0e-6 4.00 Linear 20.00
180.00 4.00 Linear 20.00
Hz
Phase
°
F Coherence w ing:vvd:+Z/Multiple
F Coherence back:vde:+Y/Multiple
0.05 -180.00
Hz Hz
LSCE
Unrestricted © Siemens AG 2019
Page 131 Siemens PLM Software
PolyMax Validation – “Noisy” FRFs
LMS PolyMAX
Unrestricted © Siemens AG 2019
Page 132 Siemens PLM Software
PolyMax Validation – FRF Synthesis with PolyMax
0.10 10.0e-3
)
g/N
Log
g/N
Log
(
(
Hz Hz
100e-6 4.00 Linear 20.00 1.00e-6 4.00 Linear 20.00
180.00 4.00 Linear 20.00 180.00 4.00 Linear 20.00
Hz Hz
Phase
Phase
°
°
-180.00 -180.00
Hz Hz
1.00
)
(m/s2)/N
Log
(
Hz
10.0e-6 10.00 Linear 64.00
180.00 10.00 Linear 64.00
Hz
Phase
°
-180.00
10.00 Hz 64.00
Test Mesh
STL File
MOVE 1
MOVE 1
MOVE 1
MOVE 2
MOVE 1
MOVE 2
MOVE 1
MOVE 2
MOVE 3
MOVE 1
MOVE 2
MOVE 3
MOVE 1
A/F
MOVE 2
A/F
MOVE 3
A/F
Hz Hz
MOVE 1
Amplitude
Amplitude
g/N
ALL
EVEN THIS ONE IS
MOVES
Amplitude
Amplitude
EFFECTED
g/N
MOVE 1
A/F
Partial Shape 1
MOVE 2
A/F
Partial Shape 2
MOVE 3
A/F
Partial Shape 3
Hz Hz
Complete Shape
Unrestricted © Siemens AG 2019
Page 156 Siemens PLM Software
Multi-Run Modal: Transmission
TIME
Functional Activities
▪ Simulate product
performance before
prototypes are
available
TIME
▪ Use single
prototype & testing
for validation
Functional Activities
Product 1 Product 2
TIME TIME
Correlation Correlation
Pretest Pretest
Analysis Analysis
FE Modeling
Process Improvement
Unrestricted © Siemens AG 2019
Page 163 Siemens PLM Software
Pretest & Correlation: Motivation
COST !
▪ Minimize Failures
▪ FEA accuracy degrades as mode order increases
▪ Simulation results are used for design decisions in Acoustics, NVH,
Durability, Loads, etc…
▪ Modal Assurance Criterion (MAC) describes how similar the shapes are for a given
mode pair using a scale of 0 to 1
.5
MAC = 100%
Unrestricted © Siemens AG 2019
Page 166 Siemens PLM Software
MAC Example
MAC = 1.5%
Unrestricted © Siemens AG 2019
Page 167 Siemens PLM Software
MAC Example
6 Points
15 Points
Parameters:
▪ MAC threshold value for matching
of FE/Test mode pairs
• Small frequency
differences
In this case:
• 9 DOFs can be removed
to improve MAC from 85%
to 95% for Mode Pair 1, 1
TEST FE
TEST FE
TEST FE
MAC WILL
NOT CHANGE
MUCH
TEST FE
MAC WILL
CHANGE A
LOT
TEST FE
Sensitivity Analysis:
•Consideration of physical
exhaust system leads
engineer to consider effect of
weld on this junction (ignored
in the FE model)
FE Model Updating:
• Manual updating
• Element thickness increased
locally in weld location
• Amount of thickness
increase guided by MAC and
Frequency correlation
MAC Correlation
• Improved from 0.69 to 0.8
Frequency Correlation
• Improved from max 15%
error to 6%
Simcenter3D Pretest & Correlation Increases Reliability of FE Models and Test Results
• FE Models accurate over entire operating frequency range
• Test results capturing all modes uniquely
• Design decisions based on the complete and correct information
16.0
(m/s2)/N
Amplitude
9.0
8.0
7.0
6.0
5.0
Damping 4.0
3.0
2.0
Hz
0.0 10.00 Linear 950.00
180 10.00 Linear 950.00
Hz
Phase
°
Mode Shape
-180
➢ Requires Multiple FRFs
10 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 950
Hz
ss sv s s s s s s s f s d 32
ss sv s s s o s s d v f s d 31
ss sd s s s s s s s f s d 30
ss sv s s s s s s s v s d 29
ss sv s s s s s d s v s d 28
ss sf s s s s s df f s o 27
ss os s s s s s ss o s 26
ss s s s s s s s s s 25
ss fs s s s s s s f s 24
ss s s s s s s ss s 23
ss s s s s s s dv s 22
ss s s s s s s df s 21
ss s s s s s s ss s 20
Amplitude
ss s s s s s s sv s 19
g/lbf
ss s s s s s s f o s 18
ss s s s s s s d s 17
ss s s s s s s s d 16
df s v v s s s d f 15
sd s ov s s s v d 14
dd s o v s s s d f 13
df s v s s s f f 12
f s s s s s s f 11
df s o v s o o o 10
f s o o d o o 9
o o f o o 8
f s o 7
o o 6
o 5
63.3e-3
85.1 Linear 837
Hz