Module 06: Response Spectrum Analysis
1 © 2016 ANSYS, Inc. September 29, 2016
ANSYS Mechanical Linear and Nonlinear Dynamics
Release 17.0
Module 06 Topics
© 2016 ANSYS, Inc.
September 29, 2016
2
Analysis
A. What is Response Spectrum Analysis?
© 2016 ANSYS, Inc.
September 29, 2016
3
... What is Response Spectrum Analysis?
© 2016 ANSYS, Inc.
September 29, 2016
4
... What is Response Spectrum Analysis?
© 2016 ANSYS, Inc.
September 29, 2016
Large model
Long transient
Large model Mode extraction
↓
Mode shapes
Small model Long transient
↓
Response spectrum
Combined solution Fast, approximate
Full solution Slow, accurate
Large model
Long transient
5
B. Generating the Response Spectrum
© 2016 ANSYS, Inc.
September 29, 2016
6
... Generating the Response Spectrum
© 2016 ANSYS, Inc.
September 29, 2016
process can be described as follows:
•
•
note the maximum absolute amplitude over time
m
k
c
k m
ω ≈
e.g., for oscillator with frequency ω = 30 Hz
max absolute value over time
S = 95 m/s2
7
... Generating the Response Spectrum
© 2016 ANSYS, Inc.
September 29, 2016
ω = 30 Hz, S = 95 m/s2
ω = 50 Hz, S = 138 m/s2
ω = 70 Hz, S = 86 m/s2
8
... Generating the Response Spectrum
© 2016 ANSYS, Inc.
September 29, 2016
9
... Generating the Response Spectrum
© 2016 ANSYS, Inc.
September 29, 2016
10
... Generating the Response Spectrum
© 2016 ANSYS, Inc.
11
September 29, 2016
by multiplying or dividing by the frequency.
– remember to convert frequency units; ω rad/s = 2πf Hz
a
= S /(2πf )2
Sd = Sv /(2πf )
a
v d
= S /(2πf )
S = S (2πf )
v
a d
= S (2πf )
S = S (2πf )2
C. Types of Analyses
© 2016 ANSYS, Inc.
12
September 29, 2016
D. Single Point Response Spectrum (SPRS) Analysis
© 2016 ANSYS, Inc.
13
September 29, 2016
© 2016 ANSYS, Inc.
14
September 29, 2016
stresses in a particular direction.
mode | frequency | mode shape | | participation factor | | |
1 | ω1 | {φ}1 | | γ1 | | |
2 | ω2 | {φ}2 | | γ2 | | |
3 | ω3 | {φ}3 | | γ3 | | |
… | … | … | | … | | |
( )
D
M
T
i
i
i
i
γ = {φ} [M ]{ }
−ω [ ]+ [K ]
{φ} = {0}
2
... Participation Factor, γ
Modal Analysis
... Spectrum Values, S
© 2016 ANSYS, Inc.
15
September 29, 2016
mode | frequency | mode shape | spectrum value | participation factor | | |
1 | ω1 | {φ}1 | S1 | γ1 | | |
2 | ω2 | {φ}2 | S2 | γ2 | | |
3 | ω3 | {φ}3 | S3 | γ3 | | |
… | … | … | … | … | | |
between spectrum points.
... Mode Coefficients, A
© 2016 ANSYS, Inc.
16
September 29, 2016
depending on the type of spectrum input.
mode | frequency | mode shape | spectrum value | participation factor | mode coefficient | |
1 | ω1 | {φ}1 | S1 | γ1 | A1 | |
2 | ω2 | {φ}2 | S2 | γ2 | A2 | |
3 | ω3 | {φ}3 | S3 | γ3 | A3 | |
| | | | | | |
…
…
…
…
…
…
displacement velocity acceleration
i i
i
i i
i i
ω2
ω
S γ S γ
A =
Ai = Siγi Ai =
displacement velocity acceleration
... Response, R
© 2016 ANSYS, Inc.
17
September 29, 2016
combined using some method.
mode | frequency | mode shape | spectrum value | participation factor | mode coefficient | response |
1 | ω1 | {φ}1 | S1 | γ1 | A1 | {R}1 |
2 | ω2 | {φ}2 | S2 | γ2 | A2 | {R}2 |
3 | ω3 | {φ}3 | S3 | γ3 | A3 | {R}3 |
… | … | … | … | … | … | … |
for displacement response for velocity response
for acceleration response
2
i i i i
{R}i
{R}i
ω A {φ}
= Ai {φ}i
= ωi Ai {φ}i
{R} =
E. Mode Combination Methods
© 2016 ANSYS, Inc.
18
September 29, 2016
same sign.
… Mode Combination Methods
© 2016 ANSYS, Inc.
19
September 29, 2016
… Square Root of the Sum of the Squares (SRSS) Method
© 2016 ANSYS, Inc.
September 29, 2016
∑
=
N
i
{R}
{R}= {R} + {
i =1
2
2
N
2
2
1 2
R} + + {R}
mode | frequency | mode shape | spectrum value | participation factor | mode coefficient | response |
1 | ω1 | {φ}1 | S1 | γ1 | A1 | {R}1 |
2 | ω2 | {φ}2 | S2 | γ2 | A2 | {R}2 |
3 | ω3 | {φ}3 | S3 | γ3 | A3 | {R}3 |
| | | | | | |
20
…
…
…
…
…
…
…
… Square Root of the Sum of the Squares (SRSS) Method
© 2016 ANSYS, Inc.
21
September 29, 2016
frequencies.
all modes.
… Sufficiently-Spaced Modes
© 2016 ANSYS, Inc.
22
September 29, 2016
ω1 ω2 ω3 ω4 ω5 ω6 ω7 ω8
N
i
{R}
{R}=
∑
i =1
2
… Closely-Spaced Modes (Correlated)
© 2016 ANSYS, Inc.
23
September 29, 2016
method
is not
applicable.
ω1
ω2 ω3
ω4 ω5 ω6
ω7 ω8
… Closely-Spaced Modes (Correlated)
© 2016 ANSYS, Inc.
24
September 29, 2016
critical damping ratio:
other
… Closely-Spaced Modes (Correlated)
© 2016 ANSYS, Inc.
25
September 29, 2016
between modes.
0.0 ≤ ε ≤ 1.0
ε = 0:
ε = 1:
0.0 < ε < 1.0:
uncorrelated
fully correlated
partially correlated
… Mode Combination Methods
© 2016 ANSYS, Inc.
26
September 29, 2016
2. Complete Quadratic Combination (CQC) method
3. Rosenblueth (ROSE) method
1
⎟
⎟
⎠
⎞ 2
⎜
⎜
⎝
⎛
N N
j
i
∑ ∑ ij i =1 j =i
{R}=
kε {R} {R}
2
1. Square Root of the Sum of Squares (SRSS) Method
1
2
⎞
⎝ i =1 ⎠
i
{R}= ⎛⎜ ∑N {R} ⎟
1
⎠
⎞ 2
⎛
N N
⎝ i =1 j =1
{R}= ⎜⎜ ∑ ∑εij {R}i {R}j ⎟⎟
Additional guidance on choosing a mode combination method will be provided in the “Recommendations” section.
F. Rigid Response
© 2016 ANSYS, Inc.
27
September 29, 2016
mid frequency
high frequency
low frequency
fSP
frequency at peak response
(spectral peak)
fZPA
frequency at rigid response
(zero period acceleration)
ZPA
u
... Spectral Regions
© 2016 ANSYS, Inc.
28
September 29, 2016
... Spectral Regions
© 2016 ANSYS, Inc.
29
September 29, 2016
low frequency
1
⎟
⎠
⎞ 2
⎜
⎜
⎛
N N
⎝ i =1 j =1
j ⎟
i
∑ ∑ ij p p
p
{R }=
ε {R } {R }
Low-Frequency Range
... Spectral Regions
High-Frequency Range
© 2016 ANSYS, Inc.
September 29, 2016
N
with themselves) and can therefore be combined algebraically, as follows:
{Rr }= ∑{Rr }i
i =1
high frequency
30
... Spectral Regions
Mid-Frequency Range
© 2016 ANSYS, Inc.
31
September 29, 2016
mid frequency
... Spectral Regions
Mid-Frequency Range
© 2016 ANSYS, Inc.
32
September 29, 2016
0.0 ≤ α ≤ 1.0
α = 0:
α = 1:
0.0 < α < 1.0:
periodic rigid
part periodic, part rigid
… Lindley-Yow method
© 2016 ANSYS, Inc.
33
September 29, 2016
ai
i
S
α = ZPA
α=0
α=1
ZPA
ZPA: acceleration at zero period
Sai: acceleration at the ith frequency
… Gupta method
© 2016 ANSYS, Inc.
34
September 29, 2016
fZPA
f1 f2
⎪
⎪⎩
⎪
⎧
2
1
2 1
2
1
1
0
3
for
for f ≤ f ≤ f
for
⎨ln( f / f )
ln( f / f )
f
f =
fi ≥ f2
i
i 1
fi ≤ f1
i
= ( f1 + 2 fZPA )
v,max
Sa,max
α =
2πS
α=0
α=1
ZPA
Rigid Response
© 2016 ANSYS, Inc.
35
September 29, 2016
– Affects all modes with f > f1.
Rigid Response Calculations
© 2016 ANSYS, Inc.
36
September 29, 2016
mode | frequency | spectrum value | response | | | |
1 | ω1 | S1 | {R}1 | | | |
2 | ω2 | S2 | {R}2 | | | |
3 | ω3 | S3 | {R}3 | | | |
… | … | … | … | | | |
… Rigid Response Calculations
© 2016 ANSYS, Inc.
37
September 29, 2016
according to the method selected.
mode | frequency | spectrum value | response | rigid response coefficient | | |
1 | ω1 | S1 | {R}1 | α1 | | |
2 | ω2 | S2 | {R}2 | α2 | | |
3 | ω3 | S3 | {R}3 | α3 | | |
… | … | … | … | … | | |
Lindley - Yow :
i
i
i
ln( f / f )
/ f1 )
ln( f
S
ai
2 1
0.0 ≤ αi ≤ 1.0
Gupta : α =
α = ZPA
… Rigid Response Calculations
© 2016 ANSYS, Inc.
38
September 29, 2016
coefficient.
Now, the periodic modes and rigid modes will be treated separately.
mode | frequency | spectrum value | response | rigid response coefficient | periodic component | rigid compone nt |
1 | ω1 | S1 | {R}1 | α1 | {Rp}1 | {Rr}1 |
2 | ω2 | S2 | {R}2 | α2 | {Rp}2 | {Rr}2 |
3 | ω3 | S3 | {R}3 | α3 | {Rp}3 | {Rr}3 |
… | … | … | … | … | … | … |
periodic component rigid component
2
i i
i
p
R
{Rr }i α{R}ii
{R } = 1−α { }
=
… Rigid Response Calculations
© 2016 ANSYS, Inc.
39
September 29, 2016
ω1
ω2
ω3 ω4 ω5 ω6
mode | frequency | response | rigid response coefficient | periodic component | rigid component |
1 | 0.5 | {R}1 | 0.0 | {R}1 | {0} |
2 | 1.4 | {R}2 | 0.0 | {R}2 | {0} |
3 | 3.6 | {R}3 | 0.18 | 0.98 {R}3 | 0.18 {R}3 |
4 | 6.4 | {R}4 | 0.61 | 0.79 {R}4 | 0.61 {R}4 |
5 | 12 | {R}5 | 1.0 | {0} | {R}5 |
6 | 25 | {R}6 | 1.0 | {0} | {R}6 |
periodic
transition
rigid
… Rigid Response Calculations
© 2016 ANSYS, Inc.
September 29, 2016
mode | frequency | spectrum value | response | rigid response coefficient | periodic component | rigid component |
1 | ω1 | S1 | {R}1 | α1 | {Rp}1 | {Rr}1 |
2 | ω2 | S2 | {R}2 | α2 | {Rp}2 | {Rr}2 |
3 | ω3 | S3 | {R}3 | α3 | {Rp}3 | {Rr}3 |
… | … | … | … | … | … | … |
{ }
1
(SRSS, CQC, or ROSE).
– Recall: CQC or ROSE are used if closely spaced modes are present.
SRSS CQC ROSE
⎟
⎟
⎠
1
⎞ 2
⎜
⎜
⎛
⎟
⎟
⎠
1
⎞ 2
⎜
⎜
⎝
⎛
2 ⎞ 2
⎛
N N
⎝ i =1 j =1
p
i j
∑ ∑ ij p
p
N N
i =1 j =1
p
i j
∑ ∑ ij p
p
N
⎝ i =1
i
∑ p
p
{R }=
{R }=
R ⎟
⎠
{R }= ⎜
ε {R } {R }
kε {R } {R }
40
… Rigid Response Calculations
{Rr }= ∑{Rr }i
i =1
© 2016 ANSYS, Inc.
41
September 29, 2016
mode | frequency | spectrum value | response | rigid response coefficient | periodic component | rigid component |
1 | ω1 | S1 | {R}1 | α1 | {Rp}1 | {Rr}1 |
2 | ω2 | S2 | {R}2 | α2 | {Rp}2 | {Rr}2 |
3 | ω3 | S3 | {R}3 | α3 | {Rp}3 | {Rr}3 |
… | … | … | … | … | … | … |
N
© 2016 ANSYS, Inc.
42
September 29, 2016
mode | frequency | spectrum value | response | rigid response coefficient | periodic component | rigid component |
1 | ω1 | S1 | {R}1 | α1 | {Rp}1 | {Rr}1 |
2 | ω2 | S2 | {R}2 | α2 | {Rp}2 | {Rr}2 |
3 | ω3 | S3 | {R}3 | α3 | {Rp}3 | {Rr}3 |
… | … | … | … | … | … | … |
2 2
r p
t
R } + {R }
{R }= {
… Rigid Response Calculations
G. Missing Mass Response
© 2016 ANSYS, Inc.
43
September 29, 2016
eff ,i i
T
i
i
D M
= γ2
γ = {φ} [M ]{ }
Missing Mass Response
© 2016 ANSYS, Inc.
44
September 29, 2016
mid frequency
high frequency
low frequency
R1
R2 R3
R4 R5 R6
R7 R8
R9
R10
R11
– We could have many modes extend far beyond fZPA.
ZPA
Neglected modes
… Missing Mass Response
© 2016 ANSYS, Inc.
45
September 29, 2016
then we don’t have to extract all modes above fZPA.
© 2016 ANSYS, Inc.
46
September 29, 2016
… Missing Mass Response Calculation
accurately represented by a static acceleration analysis.
{FT }= −[M ]{D}SZPA
{F }= −[M ]{φ} γ S
j j j ZPA
∑
N
j j ZPA
j
{F }= −
j =1
N
∑
j =1
[M ]{φ} γ S
© 2016 ANSYS, Inc.
47
September 29, 2016
excitation.
{
N
j j
N
M T
⎜
⎜
⎝ ⎠
⎛ ⎞
{F }= [M ]
{F }= {F }−
∑
j =1
∑ j j =1
φ} γ − {D}⎟⎟SZPA
M M
−1
{R }= [K ] {F }
… Missing Mass Response Calculation
© 2016 ANSYS, Inc.
48
September 29, 2016
components.
N
{Rr }= ∑{Rr }i + {RM }
i =1
2
2
r
p
{R } + {R }
{R}=
… Missing Mass Response Calculation
H. Multi-Point Response Spectrum (MPRS) Analysis
© 2016 ANSYS, Inc.
49
September 29, 2016
and mass).
… Multi-Point Response Spectrum
© 2016 ANSYS, Inc.
September 29, 2016
{RMPRS} is the total response of the MPRS analysis
{RSPRS}1 is the total response of the SPRS analysis for spectrum 1
{RSPRS}2 is the total response of the SPRS analysis for spectrum 2
etc
2
50
2
1 2
SPRS
SPRS
R
{R } + { } +
{RMPRS }=
I. Recommendations
© 2016 ANSYS, Inc.
51
September 29, 2016
response spectrum system at the solution level.
Recommended Solution Procedure
© 2016 ANSYS, Inc.
52
September 29, 2016
ω1 ω2 ω3
ω4 ω5 ω6
ω7 ω8
Recommended Solution Procedure
© 2016 ANSYS, Inc.
53
September 29, 2016
Analysis Settings
© 2016 ANSYS, Inc.
54
September 29, 2016
analysis is therefore a prerequisite.
– This will happen automatically when the Modal and Response Spectrum analysis systems are linked in the Project Schematic as shown previously; however, depending on the sequence of events leading up to combined solution of the linked systems, you may end up having to solve the Modal analysis twice. This situation can be prevented by setting “Future Analysis” to “MSUP Analyses” in the Analysis Settings Details of the Modal system before its first solution.
… Analysis Settings
© 2016 ANSYS, Inc.
55
September 29, 2016
… Analysis Settings
© 2016 ANSYS, Inc.
56
September 29, 2016
… Loads and Supports
© 2016 ANSYS, Inc.
57
September 29, 2016
… Loads and Supports
© 2016 ANSYS, Inc.
58
September 29, 2016
… Loads and Supports
© 2016 ANSYS, Inc.
59
September 29, 2016
… Loads and Supports
© 2016 ANSYS, Inc.
September 29, 2016
missing mass responses, but is higher than that of periodic responses.
RS Acceleration RS Velocity RS Displacement
Gupta: requires f1 and f2
Lindley: requires ZPA
60
Gupta: requires f1 and f2
Gupta: requires f1 and f2
… Results
© 2016 ANSYS, Inc.
61
September 29, 2016
shear strain and stress, and equivalent stress.
… Results
© 2016 ANSYS, Inc.
62
September 29, 2016
Workshop 06.1: Suspension Bridge
© 2016 ANSYS, Inc.
63
September 29, 2016
Module 06: Response Spectrum Analysis
END
64 © 2016 ANSYS, Inc. September 29, 2016