Topic 3 (Fema) Sdof

Structural Dynamics of
Linear Elastic Single-Degree-of-Freedom

(SDOF) Systems
Instructional Material Complementing FEMA 451, Design Examples
SDOF Dynamics 3 - 1
Structural Dynamics
Equations of motion for SDOF structures

Structural frequency and period of vibration
Behavior under dynamic load
Dynamic magnification and resonance
Effect of damping on behavior
Linear elastic response spectra
SDOF Dynamics 3 - 2
Importance in Relation to ASCE 7-05
Ground motion maps provide ground

accelerations in terms of response spectrum
coordinates.
Equivalent lateral force procedure gives base

shear in terms of design spectrum and period
of vibration.
Response spectrum is based on 5% critical

damping in system.
Modal superposition analysis uses design

response spectrum as basic ground motion
input.
SDOF Dynamics 3 - 3
Idealized SDOF Structure

F(t)
Mass
F ( t ), u ( t )
Damping
Stiffness
u(t)
SDOF Dynamics 3 - 4
Equation of Dynamic Equilibrium

f I (t )
0 .5 f S ( t )
fD (t )
F (t )
0 .5 f S ( t )
F (t ) f I (t ) fD (t ) fS (t ) = 0
fI (t ) + fD (t ) + fS (t ) = F (t )
SDOF Dynamics 3 - 5
Observed Response of Linear SDOF

Applied Force, kips
40
0
-40
0.00
0.20
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
Displacement, in
0.50
0.00
-0.50
0.00
15.00
0.20
Velocity, in/sec
0.00
-15.00
0.00
0.20
Acceleration, in/sec2
400.00
0.00
-400.00
0.00
0.20
Time, sec
SDOF Dynamics 3 - 6

(Development of Equilibrium Equation)
Spring Force, kips
Inertial Force, kips
Damping Force, Kips
30.00
4.00
50.00
15.00
2.00
25.00
0.00
0.00
0.00
-15.00
-2.00
-25.00
-30.00
-0.60
-4.00
-20.00
-50.00
-500
-0.30
0.00
0.30
Displacement, inches
Slope = k
= 50 kip/in
f S ( t ) = k u( t )
0.60
-10.00
0.00
10.00
20.00
Velocity, In/sec
Slope = c
= 0.254 kip-sec/in
f D ( t ) = c u& ( t )
-250
250
500
Slope = m
= 0.130 kip-sec2/in
f I ( t ) = m u&&( t )
SDOF Dynamics 3 - 7

f I (t )
0 .5 f S ( t )
fD (t )
F (t )
0 .5 f S ( t )
fI (t ) + fD (t ) + fS (t ) = F (t )
m u&&( t ) + c u& ( t ) + k u ( t ) = F ( t )
SDOF Dynamics 3 - 8
Properties of Structural Mass

Internal Force
Mass
M
1.0
Acceleration
Includes all dead weight of structure

May include some live load
Has units of force/acceleration
SDOF Dynamics 3 - 9
Damping
Damping Force
Properties of Structural Damping
C
1.0
Velocity
In absence of dampers, is called inherent damping

Usually represented by linear viscous dashpot
Has units of force/velocity
SDOF Dynamics 3 - 10
Damping
Damping Force
Properties of Structural Damping (2)
AREA =
ENERGY
DISSIPATED
Displacement
Damping vs displacement response is

elliptical for linear viscous damper.
Stiffness
Spring Force
Properties of Structural Stiffness
K
1.0
Displacement
Includes all structural members

May include some seismically nonstructural members
Requires careful mathematical modelling
Has units of force/displacement
Spring Force
Stiffness
Properties of Structural Stiffness (2)
AREA =
ENERGY
DISSIPATED
Displacement
Is almost always nonlinear in real seismic response

Nonlinearity is implicitly handled by codes
Explicit modelling of nonlinear effects is possible
Undamped Free Vibration

Equation of motion:
Initial conditions:
Assume:
u ( t ) = A sin( t ) + B cos( t )
A=
u& 0
u (t ) =
u& 0
Solution:
m u&&( t ) + k u( t ) = 0
u& 0 u 0
B = u0
k
m
sin( t ) + u 0 cos( t )
Undamped Free Vibration (2)

u& 0
T = 0.5 sec
1.0
3
2
1
0
-1
-2
-3
u0
0.0
0.5
1.0
1.5
2.0
Time, seconds
Circular Frequency
(radians/sec)
k
m
Cyclic Frequency
(cycles/sec, Hertz)
f =
2
Period of Vibration
(sec/cycle)
1
2
T =
=
f
Approximate Periods of Vibration

(ASCE 7-05)
Ta = C h
x
t n
Ct =
Ct =
Ct =
Ct =
0.028, x = 0.8
0.016, x = 0.9
0.030, x = 0.75
0.020, x = 0.75
for steel moment frames

for concrete moment frames
for eccentrically braced frames
for all other systems
Note: This applies ONLY to building structures!
T = 0.1N
a
For moment frames < 12 stories in height, minimum

story height of 10 feet. N = number of stories.
Empirical Data for Determination

of Approximate Period for Steel Moment Frames
Ta = 0.028hn0.8
Periods of Vibration of Common Structures

20-story moment resisting frame
T = 1.9 sec
T = 1.1 sec
T = 0.15 sec
20-story braced frame

T = 1.3 sec
T = 0.8 sec
T = 0.1 sec
Gravity dam
Suspension bridge
T = 0.2 sec
T = 20 sec
Adjustment Factor on Approximate Period

(Table 12.8-1 of ASCE 7-05)
T = Ta Cu Tcomputed
SD1
> 0.40g
0.30g
0.20g
0.15g
< 0.1g
Cu
1.4
1.4
1.5
1.6
1.7
Applicable ONLY if Tcomputed comes from a properly

substantiated analysis.
Which Period of Vibration to Use

in ELF Analysis?
If you do not have a more accurate period
(from a computer analysis), you must use T = Ta.
If you have a more accurate period from a computer
analysis (call this Tc), then:
if Tc > CuTa
use T = CuTa
if Ta < Tc < TuCa use T = Tc

if Tc < Ta
use T = Ta
Damped Free Vibration
m u&&( t ) + c u& ( t ) + k u ( t ) = 0
Initial conditions: u0
u& 0
st
Assume: u ( t ) = e
Equation of motion:
Solution:
u( t ) = e
u& 0 + u 0
sin( D t )
u 0 cos( D t ) +
D
c
c
=
=
2m
cc
D = 1
Damping in Structures
c
c
=
=
2m
cc
cc is the critical damping constant.
is expressed as a ratio (0.0 < < 1.0) in computations.

Sometimes is expressed as a% (0 < < 100%).
Displacement, in
Time, sec
Response of Critically Damped System, =1.0 or 100% critical

True damping in structures is NOT viscous. However, for low
damping values, viscous damping allows for linear equations
and vastly simplifies the solution.
Spring Force, kips
Damping Force, Kips
30.00
4.00
50.00
15.00
2.00
25.00
0.00
0.00
0.00
-15.00
-2.00
-25.00
-30.00
-0.60
-4.00
-20.00
-50.00
-500
-0.30
0.00
0.30
0.60
-10.00
0.00
10.00
20.00
Velocity, In/sec
-250
250
2
Acceleration, in/sec
500
Damped Free Vibration (2)

3
2
1
0
-1
-2
-3
0% Damping
10% Damping
20% Damping
0.0
0.5
1.0
1.5
2.0
Time, seconds
Damping in Structures (2)

Welded steel frame
Bolted steel frame
= 0.010
= 0.020
Uncracked prestressed concrete

Uncracked reinforced concrete
Cracked reinforced concrete
= 0.015
= 0.020
= 0.035
Glued plywood shear wall

Nailed plywood shear wall
= 0.100
= 0.150
Damaged steel structure

Damaged concrete structure
= 0.050
= 0.075
Structure with added damping
= 0.250

Inherent damping
is a structural (material) property

independent of mass and stiffness
Inherent = 0.5 to 7.0% critical

Added damping
C
is a structural property dependent on

mass and stiffness and
damping constant C of device
Added = 10 to 30% critical

Measuring Damping from Free Vibration Test

For all
damping values
u1
=
ln
u2
u1
Amplitude
0.5
u2
u3
2
1 2
u0 e
-0.5
-1
0.00
0.50
1.00
For very low

damping values
1.50
2.00
2.50
3.00
u1 u 2
2 u2
Time, Seconds
Undamped Harmonic Loading

Equation of motion:
= frequency of the forcing function
Force, Kips
T =
150
100
50
0
-50
-100
-150
0.00
m u&&( t ) + k u ( t ) = p 0 sin( t )
0.25
= 0.25 sec
po=100 kips
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Time, Seconds
Undamped Harmonic Loading (2)

&&( t ) + k u ( t ) = p 0 s in ( t )
Equation of motion: m u
Assume system is initially at rest:
Particular solution: u ( t ) = C s i n ( t )
Complimentary solution: u( t ) = A sin(t ) + B cos(t )
Solution:
p0
u (t ) =
sin( t ) sin( t )
2
k 1 ( / )

Define
Loading frequency
Structures natural frequency
Dynamic magnifier
Transient response
(at structures frequency)
p0
1
u( t ) =
sin( t ) sin( t ) )
2 (
k 1
Steady state
Static displacement
response
(at loading frequency)
= 4 rad / sec = 2 rad / sec
Steady state
response (in.)
Transient
response (in.)
Total response
(in.)
200
100
0
-100
-200
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
10
5
0
-5
-10
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
10
5
0
-5
-10
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
10
5
0
-5
-10
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
sp ace e t,
Loading (kips)
= 0.5 uS = 5.0 in.
Time, seconds
4 rad / sec = 4 rad / sec

Loading
(kips)
150
100
50
0
-5 0
-1 0 0
-1 5 0
0 .0 0
= 0.99 u S = 5.0 in.
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0.25
0.50
0.75
1.00
1.25
1.50
1.7 5
2.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
500
Steady state
response (in.)
250
0
-250
-500
0.00
500
250
Transient
response (in.)
0
-250
-500
0.00
80
0
p
Total response
(in.)
40
-40
-80
0.00
T ime, seconds
Undamped Resonant Response Curve

80
2 uS
Displacement, in.
40
-40
Linear envelope
-80
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Time, seconds

Loading (kips)
150
100
50
0
-5 0
-1 0 0
-1 5 0
0 .0 0
= 1.01
u S = 5.0 in.
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
500
Steady state
response (in.)
250
0
-2 5 0
-5 0 0
0 .0 0
500
Transient
response (in.)
250
0
-2 5 0
-5 0 0
0 .0 0
80
40
0
p
Total response
(in.)
-4 0
-8 0
0 .0 0
2 .0 0
T im e , s e c o n d s
= 4 rad / sec
Loading (kips)
150
100
50
0
-50
-100
-150
0.00
= 2.0 uS = 5.0 in.
= 8 rad / sec
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1.0 0
1.2 5
1 .5 0
1 .75
2 .00
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
Steady state
response (in.)
3
0
-3
-6
0 .0 0
6
3
Transient response
(in.)
0
-3
-6
0 .0 0
6
3
0
p
Total response
(in.)
-3
-6
0 .0 0
Response Ratio: Steady State to Static

(Signs Retained)
Magnification Factor 1/(1- 2)
12.00
8.00
In phase
4.00
Resonance
0.00
-4.00
180 degrees out of phase
-8.00
-12.00
0.00
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio

(Absolute Values)
12.00
Resonance
10.00
8.00
6.00
4.00
Slowly
loaded
Rapidly
loaded
2.00
1.00
0.00
0.00
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio
Damped Harmonic Loading

Equation of motion:
m u&&( t ) + cu& ( t ) + k u ( t ) = p 0 sin( t )
Force, Kips
T =
150
100
50
0
-50
-100
-150
0.00
= 0.25 sec
po=100 kips
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Time, Seconds

Equation of motion:
m u&&( t ) + cu& ( t ) + k u ( t ) = p 0 sin( t )

Assume system is initially at rest
Particular solution:
u ( t ) = C sin( t ) + D cos( t )
Complimentary solution:
u ( t ) = e t [ A sin( D t ) + B cos( D t ) ]
Solution:
+ C sin( t ) + D cos( t )
c
2 m
D = 1 2

Transient response at structures frequency
(eventually damps out)
u(t ) = e
[ A sin( D t ) + B cos( D t ) ]
C sin(t ) + D cos( t )
Steady state response,
at loading frequency
po
1 2
C=
k (1 2 ) 2 + (2 ) 2
po
2
D=
k (1 2 ) 2 + (2 ) 2
Damped Harmonic Loading (5% Damping)

BETA=1 (Resonance)
Beta=0.5
Beta=2.0
Displacement Amplitude, Inches
50
40
30
20
10
0
-10
-20
-30
-40
-50
0.00
1.00
2.00
3.00
4.00
5.00
Time, Seconds

50
40
1
Static
2
30
20
10
0
-10
-20
-30
-40
-50
0.00
1.00
2.00
3.00
4.00
5.00
Time, Seconds
Harmonic Loading at Resonance

Effects of Damping
200
150
100
50
0
-50
-100
-150
-200
0.00
1.00
2.00
3.00
4.00
5.00
Time, Seconds
0% Damping
%5 Damping
14.00
Resonance
0.0% Damping
5.0 % Damping
10.0% Damping
25.0 % Damping
Dynamic Response Amplifier
12.00
10.00
8.00
RD =
6.00
1
(1 2 ) 2 + ( 2 ) 2
4.00
2.00
Slowly
loaded
0.00
0.00
0.50
Rapidly
loaded
1.00
1.50
2.00
2.50
3.00
Frequency Ratio,
Summary Regarding Viscous Damping

in Harmonically Loaded Systems
For systems loaded at a frequency near their

natural frequency, the dynamic response
exceeds the static response. This is referred to
as dynamic amplification.
An undamped system, loaded at resonance, will

have an unbounded increase in displacement
over time.

Damping is an effective means for dissipating energy in
the system. Unlike strain energy, which is recoverable,
dissipated energy is not recoverable.
A damped system, loaded at resonance, will have a

limited displacement over time with the limit being (1/2)
times the static displacement.
Damping is most effective for systems loaded at or near

resonance.
CONCEPT of ENERGY STORED

and Energy DISSIPATED
F
Energy
Stored
Energy
Dissipated
2
1
LOADING
YIELDING
Energy
Recovered
2
u
UNLOADING
Total
Energy
Dissipated
3
u
UNLOADED
General Dynamic Loading
F(t)
Time, T

Solution Techniques
Fourier transform
Duhamel integration
Piecewise exact
Newmark techniques
All techniques are carried out numerically.
Piecewise Exact Method
dF
F ( ) = Fo +
dt
Fo
dF
dt
dt

Initial conditions
u o , 0 = 0 u&o , 0 = 0
Determine exact solution for 1st time step
u1 = u ( )
u&1 = u& ( )
u&&1 = u&&( )
Establish new initial conditions
u o ,1 = u ( )
u&0 ,1 = u& ( )
LOOP
Obtain exact solution for next time step
u 2 = u ( )
u& 2 = u& ( )
u&&2 = u&&( )

Advantages:
Exact if load increment is linear

Very computationally efficient
Disadvantages:
Not generally applicable for inelastic behavior

Note: NONLIN uses the piecewise exact method for
response spectrum calculations.
Newmark Techniques
Proposed by Nathan Newmark

General method that encompasses a family of different
integration schemes
Derived by:
Development of incremental equations of motion
Assuming acceleration response over short time step
Newmark Method
Advantages:
Works for inelastic response

Disadvantages:
Potential numerical error

Note: NONLIN uses the Newmark method for
general response history calculations
Development of Effective Earthquake Force
GROUND ACC, g
0.40
0.20
0.00
-0.20
-0.40
0.00
1.00
2.00
3.00
4.00
5.00
6.00
TIME, SECONDS
Earthquake Ground Motion, 1940 El Centro

0.3
0.2
0.1
40
0
-0.1
-0.2
-0.3
0
10
20
30
40
50
60
Time (sec)
Ground Displacement (cm)
15
10
Ground Velocity (cm/sec)
Ground Acceleration (g's)
0.4
30
20
10
0
-10
-20
-30
0
10
20
30
40
50
Time (sec)
0
-5
-10
-15
0
10
20
30
40
50
60
Many ground motions now

are available via the
Internet.
Time (sec)
60
u&&r
Ground Acceleration Response History
0.40
GROUND ACC, g
u&&g
u&&t
0.20
0.00
-0.20
-0.40
0.00
1.00
2.00
3.00
4.00
5.00
TIME, SECONDS
m[u&&g ( t ) + u&&r ( t )] + c u&r ( t ) + k ur ( t ) = 0

mu&&r ( t ) + c u&r ( t ) + k ur ( t ) = mu&&g ( t )
6.00
Simplified form of Equation of Motion:

mu&&r (t ) + cu&r (t ) + kur (t ) = mu&&g (t )
Divide through by m:
c
k
u&&r (t ) + u&r (t ) + ur (t ) = u&&g (t )
m
m
Make substitutions:
c
= 2
m
k
=2
m
Simplified form:
u&&r (t ) + 2 u&r (t ) + 2ur (t ) = u&&g (t )

For a given ground motion, the response

history ur(t) is function of the structures
frequency and damping ratio .
Structural frequency
u&&r (t ) + 2 u& r (t ) + u r (t ) = u&&g (t )

2
Damping ratio
Ground motion acceleration history
Response to Ground Motion (1940 El Centro)

0.4
Excitation applied to structure

with given and
0.3
0.2
0.1
0
-0.1
SOLVER
-0.2
-0.3
0
10
20
30
40
Time (sec)
60
6
Structural Displacement (in)
Change in ground motion

or structural parameters
and requires recalculation of structural
response
50
Computed response
4
2
0
-2
-4
Peak displacement
-6
0
10
20
30
Time (sec)
40
50
60
The Elastic Displacement Response Spectrum

An elastic displacement response spectrum is a plot
of the peak computed relative displacement, ur, for an
elastic structure with a constant damping , a varying
fundamental frequency (or period T = 2/ ), responding
to a given ground motion.
5% damped response spectrum for structure
responding to 1940 El Centro ground motion
DISPLACEMENT, inches
16
12
0
0
10
PERIOD, Seconds
Computation of Response Spectrum for

El Centro Ground Motion
Displacement, Inches
0.08
0.06
Computed response
0.04
0.02
0.00
-0.02
-0.04
-0.06
-0.08
10
11
12
Time, Seconds
10.00
Elastic response spectrum
= 0.05
T = 0.10 sec
Umax= 0.0543 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
Computation of Response Spectrum

for El Centro Ground Motion
0.40
0.30
Computed response
0.20
0.10
0.00
-0.10
-0.20
-0.30
-0.40
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.20 sec
Umax = 0.254 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds

0.80
0.60
Computed response
0.40
0.20
0.00
-0.20
-0.40
-0.60
-0.80
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.30 sec
Umax = 0.622 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds

1.20
0.90
Computed response
0.60
0.30
0.00
-0.30
-0.60
-0.90
-1.20
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.40 sec
Umax = 0.956 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds

2.40
1.80
Computed response
1.20
0.60
0.00
-0.60
-1.20
-1.80
-2.40
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.50 sec
Umax = 2.02 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds

3.20
2.40
1.60
Computed response
0.80
0.00
-0.80
-1.60
-2.40
-3.20
0
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.60 sec
Umax= -3.00 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
Complete 5% Damped Elastic Displacement

Response Spectrum for El Centro
Ground Motion
12.00
10.00
8.00
6.00
4.00
2.00
0.00
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Period, Seconds
Development of Pseudovelocity
Response Spectrum
35.00
5% damping
Pseudovelocity, in/sec
30.00
25.00
20.00
15.00
10.00
PSV (T ) D
5.00
0.00
0.0
1.0
2.0
3.0
4.0
Period, Seconds
Development of Pseudoacceleration
Response Spectrum
400.0
5% damping
Pseudoacceleration, in/sec
350.0
300.0
250.0
200.0
PSA (T ) 2 D
150.0
100.0
50.0
0.0
0.0
1.0
2.0
3.0
4.0
Period, Seconds
Note About the Pseudoacceleration Response Spectrum

The pseudoacceleration response spectrum represents the total
acceleration of the system, not the relative acceleration. It is nearly
identical to the true total acceleration response spectrum for lightly
damped structures.
400.0
5% damping
350.0
Peak ground
acceleration
300.0
250.0
200.0
150.0
100.0
50.0
0.0
0.0
1.0
2.0
3.0
4.0
Period, Seconds
u&&r
0.40
GROUND ACC, g
u&&g
u&&t
PSA is TOTAL Acceleration!
0.20
0.00
-0.20
-0.40
0.00
1.00
2.00
3.00
4.00
5.00
TIME, SECONDS

6.00
Difference Between Pseudo-Acceleration

and Total Acceleration
(System with 5% Damping)
300.00
Acceleration (in/sec )
350.00
250.00
200.00
150.00
100.00
50.00
0.00
0.1
1
Period (sec)
Total Acceleration
10
Pseudo-Acceleration
Difference Between Pseudovelocity

and Relative Velocity
40
Ve locity (in/se c)
35
30
25
20
15
10
5
0
0.1
1
Period (sec)
Relative Velocity
10
Pseudo-Velocity
Displacement Response Spectra

for Different Damping Values
Damping
25.00
0%
5%
10%
20%
20.00
15.00
10.00
5.00
0.00
0.0
1.0
2.0
3.0
4.0
5.0
Period, Seconds
Pseudoacceleration Response Spectra

Damping
Pseudoacceleration, g
4.00
0%
5%
10%
20%
3.00
2.00
1.00
0.00
Peak ground
acceleration
0.0
1.0
2.0
3.0
4.0
5.0
Period, Seconds
Damping Is Effective in Reducing the

Response for (Almost) Any Given Period
of Vibration
An earthquake record can be considered to be the
combination of a large number of harmonic components.
Any SDOF structure will be in near resonance with one
of these harmonic components.
Damping is most effective at or near resonance.
Hence, a response spectrum will show reductions due to
damping at all period ranges (except T = 0).
Amplitude

Response for Any Given Period of
Vibration
4.00
2.00
0.00
-2.00
-4.00
0.0
6.0
12.0 18.0 24.0 30.0 36.0 42.0 48.0 54.0 60.0 66.0 72.0 78.0 84.0 90.0
Time (sec)
Example of an artificially generated wave to

resemble a real time ground motion
accelerogram.
Generated wave obtained by combining five
different harmonic signals, each having equal
amplitude of 1.0.
The Artificial Wave Is the Sum of Five Harmonics

1
0.5
0
-0.5
-1
0.0
Amplitud
e
T = 5.0 s
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 4.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 3.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
Time (sec)

1
0.5
0
-0.5
-1
0.0
Amplitude
T = 2.0 s
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 1.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
Summation
4.00
2.00
0.00
-2.00
-4.00
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
Time (sec)
Damping Reduces the Response

at Each Resonant Frequency
14.00
0.0% Damping
5.0 % Damping
10.0% Damping
25.0 % Damping
Fourier amplitude
12.00
10.00
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio,
Frequency (Hz)
FFT curve for the combined wave
Use of an Elastic Response Spectrum

Example Structure
12.00
K = 500 k/in
M = 2000/386.4 = 5.18 k-sec2/in
= (K/M)0.5 =9.82 rad/sec
T = 2/ = 0.64 sec
5% critical damping
10.00
W = 2,000 k
8.00
6.00
4.00
2.00
0.00
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Period, Seconds
At T = 0.64 sec, displacement = 3.03 in.


400.0
K = 500 k/in
W = 2,000 k
M = 2000/386.4 = 5.18 k-sec2/in
= (K/M)0.5 =9.82 rad/sec
T = 2/ = 0.64 sec
5% critical damping
Pseudoacceleration, in/sec2
Example Structure
350.0
300.0
250.0
200.0
150.0
100.0
50.0
0.0
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
Period, Seconds
At T = 0.64 sec, pseudoacceleration = 301 in./sec2

Base shear = M x PSA = 5.18(301) = 1559 kips
4.0
Response Spectrum, ADRS Space

1.00
Diagonal lines represent
period values, T
0.80
T = 0.64s
0.60
0.40
0.20
0.00
0.00
2.00
4.00
6.00
8.00
10.00
12.00
Four-Way Log Plot of Response Spectrum

Line of increasing
displacement
100
PSEUDOVELOCITY, in/sec
D=10.0
1.0
0.1
Line of constant
displacement
10
1.0
0.01
D=
1
.01
0.1
0.001
PSV
0.1
0.1
10
100
1000
Circular Frequency
per Second
, Radiand
Circular Frequency
(radians/sec)

Line of constant
acceleration
100
PSA=1000
10000
100000
Line of increasing
acceleration
10
100
10000
PSA = PSV
10
100
1000
0.1
0.1
10
100
1000
Circular Frequency
Radiand per
Circular ,
Frequency
Second
(radians/sec)

2
100
DI
S
ec
s
n/
i
N,
IO
T
A 100
R
00
LE
PL
AC
0 EM
10
E
10
NT
,i
n
10
E
C
AC 10
0
10
0
0
1.
10
1
0.
01
0.
0.
1
01
0
0.
0.1
0.1
10
100
1000
Circular Frequency
, Radiand per
Circular Frequency
Second
(radians/sec)

Plotted vs Period
D
is
pl
ac
em
.0
10
1
00
0.
0.
00
1
01
0.
0.
01
1.00
10
0.
0.
1
0
1.
1.
0
g
n,
io
10.00
at
10
.0
er
el
cc
en
t,
in
100.00
0.10
0.01
0.10
1.00
10.00
PERIOD, Seconds
Development of an Elastic
Response Spectrum
Problems with Current Spectrum:
100.00
en
t
10
.0
D
is
.0
10
n,
pl
ac
em
tio
ra
10.00
0.
1
1.
0
0
1.
10
0.
0.
01
1.00
01
0.
1
00
0.
0.10
0.01
0.
00
1
e
el
,i
n.
c
Ac
For a given earthquake,

small variations in structural
frequency (period) can produce
significantly different results.
0.10
1.00
PERIOD, Seconds
10.00
It is for a single earthquake; other

earthquakes will have different
Characteristics.
1940 El Centro, 0.35 g, N-S

A
en
t,
in
100
10
.0
0% Damping
5% Damping
10% Damping
20* Damping
10
0.
0.
1
0
1.
1.
0
10
1
00
0.
0.
00
1
0.
01
01
0.
Pseudo Velocity, In/Sec
g
n,
io
.0
10
at
D
is
pl
ac
em
er
el
cc

0.1
0.01
0.1
1
Period, Seconds
10
5% Damped Spectra for Four California Earthquakes

Scaled to 0.40 g (PGA)
Pseuso Velocity, in/sec
Different earthquakes
100.0
will have different spectra.
El Centro
Loma Prieta
North Ridge
San Fernando
Average
10.0
1.0
0.1
0.01
0.10
1.00
10.00
Period, seconds
Smoothed Elastic Response Spectra

(Elastic DESIGN Response Spectra)
Newmark-Hall spectrum
ASCE 7 spectrum
Newmark-Hall Elastic Spectrum

100
0% Dam ping
5% Dam ping
Displacem ent (in)
10% Dam ping
10
max u&
Observations
v&& max v&&

at short T
v0
max u
max u&&
v max v
v&& 0
0.1
0.01
0.1
Period1 (sec)
10
at long T
100
Very Stiff Structure (T < 0.01 sec)
Relative displacement
Total acceleration
Zero
Ground acceleration
Very Flexible Structure (T > 10 sec)
Total acceleration
Ground displacement
Zero

100
10
.0
.0
10
g
n,
io
1.
0
0
1.
0% Damping
5% Damping
10% Damping
20* Damping
4.25 in.
0.
01
0.
1
10
0.
12.7 in/s
01
0.
0.35g
1
00
0.
0.1
0.01
Ground Maxima
0.
00
1
at
D
is
pl
ac
em
en
t,
er
el
cc
in
10
0.1
1
Period, Seconds
10
Newmarks Spectrum Amplification Factors

for Horizontal Elastic Response
Damping
% Critical
.05
1
2
3
5
7
10
20
One Sigma (84.1%)

av
ad
aa
5.10 3.84 3.04
4.38 3.38 2.73
3.66 2.92 2.42
3.24 2.64 2.24
2.71 2.30 2.01
2.36 2.08 1.85
1.99 1.84 1.69
1.26 1.37 1.38
aa
3.68
3.21
2.74
2.46
2.12
1.89
1.64
1.17
Median (50%)
av
ad
2.59 2.01
2.31 1.82
2.03 1.63
1.86 1.52
1.65 1.39
1.51 1.29
1.37 1.20
1.08 1.01

1) Draw the lines
corresponding to max
v&& , v& , v
2) Draw line max v
&&
3
4
from Tb to Tc
6
3) Draw line V
from Tc to Td
4) Draw line D
from Td to Te
max v&
max v
5) Draw connecting line

from Ta to Tb
Ta
Tb
Tc
Td
Te Tf

from Te to Tf
Spectral Response
Acceleration, Sa
ASCE 7
Uses a Smoothed Design Acceleration Spectrum
Short
period
acceleration
2
SDS
SD1
Long period
acceleration
SDS
Sa = 0.6
T + 0.4 SDS
T0
Sa = SDS
3
4
SD1
T
TS
Sa = L 2D1
T
Sa =
TS
Period, T
T = 1.0
TL
Note exceptions at larger periods

The ASCE 7 Response Spectrum

is a uniform hazard spectrum based on
probabilistic and deterministic seismic
hazard analysis.
(SDOF) Systems
SDOF Dynamics 3 - 1
This set of slides covers the fundamental concepts of structural dynamics of

linear elastic single-degree-of-freedom (SDOF) structures. A separate topic
covers the analysis of linear elastic multiple-degree-of-freedom (MDOF)
systems. A separate topic also addresses inelastic behavior of structures.
Proficiency in earthquake engineering requires a thorough understanding of
each of these topics.
FEMA 451B Topic 3 Notes
Slide 1
Structural Dynamics

SDOF Dynamics 3 - 2
This slide lists the scope of the present topic. In a sense, the majority of the
material in the topic provides background on the very important subject of
response spectra.
Slide 2

coordinates.

of vibration.

damping in system.

input.
SDOF Dynamics 3 - 3
The relevance of the current topic to the ASCE 7-05 document is provided
here. Detailed referencing to numbered sections in ASCE 7-05 is provided
in many of the slides. Note that ASCE 7-05 is directly based on the 2003
NEHRP Recommended Provisions for Seismic Regulations for New
Buildings and Other Structures, FEMA 450, which is available at no charge
from the FEMA Publications Center, 1-800-480-2520 (order by FEMA
publication number).
Slide 3

F(t)
Mass
F ( t ), u ( t )
Damping
Stiffness
u(t)
SDOF Dynamics 3 - 4
The simple frame is idealized as a SDOF mass-spring-dashpot model with a

time-varying applied load. The function u(t) defines the displacement
response of the system under the loading F(t). The properties of the
structure can be completely defined by the mass, damping, and stiffness as
shown.
The idealization assumes that all of the mass of the structure can be lumped
into a single point and that all of the deformation in the frame occurs in the
columns with the beam staying rigid. Represent damping as a simple
viscous dashpot common as it allows for a linear dynamic analysis. Other
types of damping models (e.g., friction damping) are more realistic but
require nonlinear analysis.
Slide 4

f I (t )
0 .5 f S ( t )
fD (t )
F (t )
0 .5 f S ( t )
F (t ) f I (t ) fD (t ) fS (t ) = 0
fI (t ) + fD (t ) + fS (t ) = F (t )
SDOF Dynamics 3 - 5
Here the equations of motion are shown as a force-balance. At any point in

time, the inertial, damping, and elastic resisting forces do not necessarily act
in the same direction. However, at each point in time, dynamic equilibrium
must be maintained.
Slide 5

Applied Force, kips
40
0
-40
0.00
0.20
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
Displacement, in
0.50
0.00
-0.50
0.00
15.00
0.20
Velocity, in/sec
0.00
-15.00
0.00
0.20
Acceleration,
in/sec2
400.00
0.00
-400.00
0.00
0.20
Time, sec
SDOF Dynamics 3 - 6
This slide (from NONLIN) shows a series of response histories for a SDOF
system subjected to a saw-tooth loading. As a result of the loading, the
mass will undergo displacement, velocity, and acceleration. Each of these
quantities are measured with respect to the fixed base of the structure.
Note that although the loading is discontinuous, the response is relatively
smooth. Also, the vertical lines show that velocity is zero when displacement
is maximum and acceleration is zero when velocity is maximum.
NONLIN is an educational program for dynamic analysis of simple linear and
nonlinear structures. Version 7 is included on the CD containing these
instructional materials.
Slide 6

Spring Force, kips
Damping Force, Kips
30.00
4.00
50.00
15.00
2.00
25.00
0.00
0.00
0.00
-15.00
-2.00
-25.00
-30.00
-0.60
-0.30
0.00
0.30
Slope = k
= 50 kip/in
f S ( t ) = k u( t )
0.60
-4.00
-20.00
-10.00
0.00
10.00
20.00
-50.00
-500
Velocity, In/sec
Slope = c
= 0.254 kip-sec/in
f D ( t ) = c u& ( t )
-250
250
500
Slope = m
= 0.130 kip-sec2/in
f I ( t ) = m u&&( t )
SDOF Dynamics 3 - 7
These X-Y curves are taken from the same analysis that produced the
response histories of the previous slide. For a linear system, the resisting
forces are proportional to the motion. The slope of the inertial-force vs
acceleration curve is equal to the mass. Similar relationships exist for
damping force vs velocity (slope = damping) and elastic force vs
displacement (slope = stiffness).
The importance of understanding and correct use of units cannot be over
emphasized.
Slide 7

f I (t )
0 .5 f S ( t )
fD (t )
F (t )
0 .5 f S ( t )
f I (t ) + fD (t ) + fS (t ) = F (t )
m u&&( t ) + c u& ( t ) + k u ( t ) = F ( t )
SDOF Dynamics 3 - 8
Here the equations of motion are shown in terms of the displacement,

velocity, acceleration, and force relationships presented in the previous slide.
Given the forcing function, F(t), the goal is to determine the response history
of the system.
Slide 8

Internal Force
Mass
M
1.0
Acceleration

SDOF Dynamics 3 - 9
Mass is always assumed constant throughout the response. Section 12.7.2

of ASCE 7-05 defines this mass in terms of the effective weight of the
structure. The effective weight includes 25% of the floor live load in areas
used for storage, 10 psf partition allowance, operating weight of all
permanent equipment, and 20% of the flat roof snow load when that load
exceeds 30 psf.
Slide 9
Damping
Damping Force
C
1.0
Velocity

Except for the case of added damping, real structures do not have discrete
dampers as shown. Real or inherent damping arises from friction in the
material. For cracked concrete structures, damping is higher because of the
rubbing together of jagged surfaces on either side of a crack.
In analysis, we use an equivalent viscous damper primarily because of the
mathematical convenience. (Damping force is proportional to velocity.)
Slide 10
Damping
Damping Force
AREA =
ENERGY
DISSIPATED
Displacement

The force-displacement relationship for a linear viscous damper is an ellipse.

The area within the ellipse is the energy dissipated by the damper. The
greater the energy dissipated by damping, the lower the potential for
damage in structures. This is the primary motivation for the use of added
damping systems. Energy that is dissipated is irrecoverable.
Slide 11
Stiffness
Spring Force
K
1.0
Displacement

In this topic, it is assumed that the force-displacement relationship in the

spring is linear elastic. Real structures, especially those designed according
to current seismic code provisions, will not remain elastic and, hence, the
force-deformation relationship is not linear. However, linear analysis is often
(almost exclusively) used in practice. This apparent contradiction will be
explained as this discussion progresses.
The modeling of the structure for stiffness has very significant uncertainties.
Section 12.7.3 of ASCE 7-05 provides some guidelines for modeling the
structure for stiffness.
Slide 12
Stiffness
Spring Force
AREA =
ENERGY
DISSIPATED
Displacement

This is an idealized response of a simple inelastic structure. The area within

the curve is the inelastic hysteretic energy dissipated by the yielding
material. The larger hysteretic energy in relation to the damping energy, the
greater the damage.
In this topic, it is assumed that the material does not yield. Nonlinear
inelastic response is explicitly included in a separate topic.
Slide 13

Equation of motion:
Initial conditions:
Assume:
u ( t ) = A sin( t ) + B cos( t )
A=
u& 0
u (t ) =
u& 0
Solution:
m u&&( t ) + k u( t ) = 0
u& 0 u 0
B = u0
k
m
sin( t ) + u 0 cos( t )
In this unit, we work through a hierarchy of increasingly difficult problems.

The simplest problem to solve is undamped free vibration. Usually, this type
of response is invoked by imposing a static displacement and then releasing
the structure with zero initial velocity. The equation of motion is a second
order differential equation with constant coefficients. The displacement term
is treated as the primary unknown.
The assumed response is in terms of a sine wave and a cosine wave. It is
easy to see that the cosine wave would be generated by imposing an initial
displacement on the structure and then releasing. The sine wave would be
imposed by initially shoving the structure with an initial velocity. The
computed solution is a combination of the two effects.
The quantity is the circular frequency of free vibration of the structure
(radians/sec). The higher the stiffness relative to mass, the higher the
frequency. The higher the mass with respect to stiffness, the lower the
frequency.
Slide 14

u& 0
T = 0.5 sec
1.0
3
2
1
0
-1
-2
-3
u0
0.0
0.5
1.0
1.5
2.0
Time, seconds
Circular Frequency
(radians/sec)
k
m
Cyclic Frequency
(cycles/sec, Hertz)
f =
Period of Vibration
(sec/cycle)
T =
1
2
=
f
This slide shows a computed response history for a system with an initial
displacement and velocity. Note that the slope of the initial response curve
is equal to the initial velocity (v = du/dt). If this term is zero, the free vibration
response is a simple cosine wave. Note also that the undamped motion
shown will continue forever if uninhibited. In real structures, damping will
eventually reduce the free vibration response to zero.
The relationship between circular frequency, cyclic frequency, and period of
vibration is emphasized. The period of vibration is probably the easiest to
visualize and is therefore used in the development of seismic code
provisions. The higher the mass relative to stiffness, the longer the period of
vibration. The higher the stiffness relative to mass, the lower the period of
vibration.
Slide 15

(ASCE 7-05)
Ta = Ct hnx
Ct =
Ct =
Ct =
Ct =
0.028, x = 0.8
0.016, x = 0.9
0.030, x = 0.75
0.020, x = 0.75

T = 0.1N
a

One of the first tasks in any seismic design project is to estimate the period
of vibration of the structure. For preliminary design (and often for final
design), an empirical period of vibration is used. Section 12.8.2 of ASCE 705 provides equations for estimating the period. These equations are listed
here.
Slide 16

Ta = 0.028hn0.8
Ta is based on curve-fitting of data obtained from measured response of

California buildings after small earthquakes. As will be seen later, the
smaller the period, the larger the earthquake force that must be designed for.
Hence, a lower bound empirical relationship is used.
Because the empirical period formula is based on measured response of
buildings, it should not be used to estimate the period for other types of
structure (bridges, dams, towers).
Slide 17

T = 1.9 sec
T = 1.1 sec
T = 0.15 sec

T = 1.3 sec
T = 0.8 sec
T = 0.1 sec
Gravity dam
Suspension bridge
T = 0.2 sec
T = 20 sec
This slide shows typical periods of vibration for several simple structures.
Engineers should develop a feel for what an appropriate period of vibration
is for simple building structures.
For building structures, the formula T = 0.1 in is the simplest reality check.
The period for a 10-story building should be approximately 1 sec. If a
computer analysis gives a period of 0.2 sec or 3.0 sec for a 10-story building,
something is probably amiss in the analysis.
Slide 18

T = Ta Cu Tcomputed
SD1
> 0.40g
0.30g
0.20g
0.15g
< 0.1g
Cu
1.4
1.4
1.5
1.6
1.7

In some cases, it is appropriate to remove the conservatism from the

empirical period formulas. This is done through use of the Cu coefficient.
This conservatism arises from two sources:
1. The lower bound period was used in the development of the period
formula.
2. This lower bound period is about 1/1.4 times the best-fit period.
The empirical formula was developed on the basis of data from California
buildings. Buildings in other parts of the country (e.g., Chicago) where
seismic forces are not so high will likely be larger than those for the same
building in California.
It is important to note that the larger period cannot be used without the
benefit of a properly substantiated analysis, which is likely performed on
a computer.
Slide 19

in ELF Analysis?
if Tc > CuTa
use T = CuTa

if Tc < Ta
use T = Ta
This slide shows the limitations on the use of CuTa. ASCE-7-05 will not allow
the use of a period larger than CuTa regardless of what the computer
analysis says. Similarly, the NEHRP Recommended Provisions does not
require that you use a period less than Ta.
Slide 20
m u&&( t ) + c u& ( t ) + k u ( t ) = 0
u& 0
st
Assume: u ( t ) = e
Equation of motion:
Solution:
u& + u 0
u ( t ) = e t u 0 cos( D t ) + 0
sin( D t )
D
c
c
=
2m
cc
D = 1 2
This slide shows the equation of motion and the response in damped free
vibration. Note the similarity with the undamped solution. In particular, note
the exponential decay term that serves as a multiplier on the whole
response.
Critical damping (cc) is defined as the amount of damping that will produce
no oscillation. See next slide.
The damped circular frequency is computed as shown. Note that in many
practical cases (x < 0.10), it will be effectively the same as the undamped
frequency. The exception is very highly damped systems.
Note that the damping ratio is often given in terms of % critical.
Slide 21
=
c
c
=
2m
cc

Displacement, in
Time, sec

The concept of critical damping is defined here. A good example of a

critically damped response can be found in heavy doors that are fitted with
dampers to keep the door from slamming when closing.
Slide 22
Spring Force, kips
Damping Force, Kips
30.00
4.00
50.00
15.00
2.00
25.00
0.00
0.00
0.00
-15.00
-2.00
-25.00
-30.00
-0.60
-0.30
0.00
0.30
0.60
-4.00
-20.00
-10.00
0.00
10.00
20.00
-50.00
-500
Velocity, In/sec
-250
250
500
An earlier slide is repeated here to emphasize that damping in real

structures is NOT viscous. It is frictional or hysteretic. Viscous damping is
used simply because it linearizes the equations of motion. Use of viscous
damping is acceptable for the modeling of inherent damping but should be
used with extreme caution when representing added damping or energy loss
associated with yielding in the primary structural system.
Slide 23

3
2
1
0
-1
-2
-3
0% Damping
10% Damping
20% Damping
0.0
0.5
1.0
1.5
2.0
Time, seconds
This slide shows some simple damped free vibration responses. When the
damping is zero, the vibration goes on forever. When the damping is 20%
critical, very few cycles are required for the free vibration to be effectively
damped out. For 10% damping, peak is approximately of the amplitude of
the previous peak.
Slide 24

Welded steel frame
Bolted steel frame
= 0.010
= 0.020

= 0.015
= 0.020
= 0.035

= 0.100
= 0.150

= 0.050
= 0.075
= 0.250
Some realistic damping values are listed for structures comprised of different
materials. The values for undamaged steel and concrete (upper five lines of
table) may be considered as working stress values.
Slide 25

Inherent damping


Added damping
C


The distinction between inherent damping and added damping should be

clearly understood.
Slide 26

For all
damping values
u1
ln
Amplitude
0.5
u2
u3
u1
=
u2
2
1 2
-1
0.00
For very low

damping values
u0 e t
-0.5
0.50
1.00
1.50
2.00
2.50
3.00
u1 u 2
2 u2
Time, Seconds
One of the simplest methods to measure damping is a free vibration test.

The structure is subjected to an initial displacement and is suddenly
released. Damping is determined from the formulas given. The second
formula should be used only when the damping is expect to be less than
about 10% critical.
Slide 27

Equation of motion:
Force, Kips
T =
150
100
50
0
-50
-100
-150
0.00
m u&&( t ) + k u ( t ) = p 0 sin( t )
0.25
= 0.25 sec
po=100 kips
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Time, Seconds
The next series of slides covers the response of undamped SDOF systems
to simple harmonic loading. Note that the loading frequency is given by the
omega term with the overbar. The loading period is designated in a similar
fashion.
Slide 28

&&( t ) + k u ( t ) = p 0 s in ( t )
Particular solution: u ( t ) = C s in ( t )
Solution:
u (t ) =
p0
1
k 1 ( / ) 2
sin( t ) sin( t )
This slide sets up the equation of motion for undamped harmonic loading
and gives the solution. We have assumed the system is initially at rest.
Slide 29

Define
Loading frequency
Dynamic magnifier
u( t ) =
Transient response
1
p0
(sin( t ) sin( t ) )
k 1 2
Static displacement
Steady state
response
Here we break up the response into the steady state response (at the
frequency of loading) and the transient response (at the structures own
natural frequency). Note that the term po/k is the static displacement. The
dynamic magnifier shows how the dynamic effects may increase (or
decrease) the response. This magnifier is a function of the frequency ratio
. Note that the magnifier goes to infinity if the frequency ratio is 1.0. This
defines the resonant condition.
In other words, the response is equal to the static response, times a
multiplier, times the sum of two sine waves, one in phase with the load and
the other in phase with the structures undamped natural frequency.
Slide 30
= 0.5 uS = 5.0 in.
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
10
5
0
-5
-10
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Transient
response (in.)
10
5
0
-5
-10
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Total response
(in.)
10
5
0
-5
-10
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Loading (kips)
Steady state
response (in.)
sp ace e t,
200
100
0
-100
-200
0.00
Time, seconds
This is a time-history response of a structure with a natural frequency of 4

rad/sec (f = 2 Hz, T = 0.5 sec), and a loading frequency of 2 rad/sec (f = 1
Hz, T = 1 sec), giving a frequency ratio of 0.5. The harmonic load
amplitude is 100 kips. The static displacement is 5.0 inches. Note how the
steady state response is at the frequency of loading, is in phase with the
loading, and has an amplitude greater than the static displacement. The
transient response is at the structures own frequency. In real structures,
damping would cause this component to disappear after a few cycles of
vibration.
Slide 31

Loading
(kips)
150
100
50
0
-5 0
-1 0 0
-1 5 0
0 .0 0
= 0.99 u S = 5.0 in.
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0.25
0.50
0 .7 5
1.00
1.25
1 .5 0
1.7 5
2.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
5 00
Steady state
response (in.)
2 50
0
-2 50
-5 00
0 .0 0
500
250
Transient
response (in.)
0
-250
-500
0.00
80
40
0
-40
Total response
(in.)
-80
0.00
T ime, seconds
In this slide, has been increased to 4 rad/sec, and the structure is almost
at resonance. The steady state response is still in phase with the loading,
but note the huge magnification in response. The transient response is
practically equal to and opposite the steady state response. The total
response increases with time.
If one looks casually at the steady state and transient response curves, it
appears that they should cancel out. Note, however, that the two responses
are not exactly in phase due to the slight difference in the loading and
natural frequencies. This can be seen most clearly at the time 1.75 sec into
the response. The steady state response crosses the horizontal axis to the
right of the vertical 1.75 sec line while the transient response crosses exactly
at 1.75 sec.
In real structures, the observed increased amplitude could occur only to
some limit and then yielding would occur. This yielding would introduce
hysteretic energy dissipation (apparent damping), causing the transient
response to disappear and leading to a constant, damped, steady state
response.
Slide 32

80
2 uS
Displacement, in.
40
-40
Linear envelope
-80
0.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Time, seconds
This is an enlarged view of the total response curve from the previous slide.
Note that the response is bounded within a linear increasing envelope with
the increase in displacement per cycle being 2 times the static
displacement.
Slide 33

Loading (kips)
150
100
50
0
-5 0
-1 0 0
-1 5 0
0 .0 0
= 1.01
u S = 5.0 in.
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
500
Steady state
response (in.)
250
0
-2 5 0
-5 0 0
0 .0 0
500
Transient
response (in.)
250
0
-2 5 0
-5 0 0
0 .0 0
Total response
(in.)
40
80
-4 0
-8 0
0 .0 0
2 .0 0
In this slide, the loading frequency has been slightly increased, but the
structure is still nearly at resonance. Note, however, that the steady state
response is 180 degrees out of phase with the loading and the transient
response is in phase. The resulting total displacement is effectively identical
to that shown two slides back.
Slide 34
= 4 rad / sec
Loading (kips)
150
100
50
0
-50
-100
-150
0.00
= 2.0 uS = 5.0 in.
= 8 rad / sec
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
Steady state
response (in.)
3
0
-3
-6
0 .0 0
6
3
Transient response
(in.)
0
-3
-6
0 .0 0
Total response
(in.)
-3
-6
0 .0 0
The loading frequency is now twice the structures frequency. The important
point here is that the steady state response amplitude is now less than the
static displacement.
Slide 35

(Signs Retained)
12.00
8.00
In phase
4.00
Resonance
0.00
-4.00
-8.00
-12.00
0.00
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio
This plot shows the ratio of the steady state response to the static
displacement for the structure loaded at different frequencies. At low loading
frequencies, the ratio is 1.0, indicating a nearly static response (as
expected). At very high frequency loading, the structure effectively does not
have time to respond to the loading so the displacement is small and
approaches zero at very high frequency. The resonance phenomena is very
clearly shown. The change in sign at resonance is associated with the inphase/out-of-phase behavior that occurs through resonance.
Slide 36

(Absolute Values)
12.00
Resonance
10.00
8.00
6.00
4.00
Slowly
loaded
Rapidly
loaded
2.00
1.00
0.00
0.00
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio
This is the same as the previous slide but absolute values are plotted. This
clearly shows the resonance phenomena.
Slide 37

Equation of motion:
m u&&( t ) + cu& ( t ) + k u ( t ) = p 0 sin( t )
Force, Kips
T =
150
100
50
0
-50
-100
-150
0.00
= 0.25 sec
po=100 kips
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Time, Seconds
We now introduce damping into the behavior. Note the addition of the
appropriate term in the equation of motion.
Slide 38

Equation of motion:
m u&&( t ) + cu& ( t ) + k u ( t ) = p 0 sin( t )

u ( t ) = C sin( t ) + D cos( t )
Solution:
c
2 m
D = 1 2
This slide shows how the solution to the differential equation is obtained.
The transient response (as indicated by the A and B coefficients) will damp
out and is excluded from further discussion.
Slide 39

u ( t ) = e t [ A sin( D t ) + B cos( D t ) ] +
C sin(t ) + D cos(t )
C=
po
1 2
k (1 2 ) 2 + (2 ) 2
D=
po
2
k (1 2 ) 2 + (2 ) 2
This slide shows the C and D coefficients of the steady state response.
Note that there is a component in phase with the loading (the sine term) and
a component out of phase with the loading (the cosine term). The actual
phase difference between the loading and the response depends on the
damping and frequency ratios.
Note the exponential decay term causes the transient response to damp out
in time.
Slide 40

BETA=1 (Resonance)
Beta=0.5
Beta=2.0
50
40
30
20
10
0
-10
-20
-30
-40
-50
0.00
1.00
2.00
3.00
4.00
5.00
Time, Seconds
This plot shows the response of a structure at three different loading

frequencies. Of significant interest is the resonant response, which is now
limited. (The undamped response increases indefinitely.)
Slide 41

50
40
1
Static
2
30
20
10
0
-10
-20
-30
-40
-50
0.00
1.00
2.00
3.00
4.00
5.00
Time, Seconds
For viscously damped structures, the resonance amplitude will always be

limited as shown.
Slide 42

Effects of Damping
200
150
100
50
0
-50
-100
-150
-200
0.00
1.00
2.00
3.00
4.00
5.00
Time, Seconds
0% Damping
%5 Damping
A comparison of damped and undamped responses is shown here. The

undamped response has a linear increasing envelope; the damped curve will
reach a constant steady state response after a few cycles.
Slide 43
14.00
Resonance
0.0% Damping
5.0 % Damping
10.0% Damping
25.0 % Damping
12.00
10.00
8.00
RD =
6.00
1
(1 2 ) 2 + ( 2 ) 2
4.00
2.00
Slowly
loaded
0.00
0.00
0.50
Rapidly
loaded
1.00
1.50
2.00
2.50
3.00
Frequency Ratio,
This plot shows the dynamic magnification for various damping ratios. For
increased damping, the resonant response decreases significantly. Note
that for slowly loaded structures, the dynamic amplification is 1.0 (effectively
static). For high frequency loading, the magnifier is zero.
Note also that damping is most effective at or near resonance (0.5 < <
2.0).
Slide 44



over time.
A summary of some of the previous points is provided.
Slide 45



resonance.
Summary continued.
Slide 46

F
Energy
Stored
Energy
Dissipated
2
1
LOADING
YIELDING
Energy
Recovered
2
u
UNLOADING
Total
Energy
Dissipated
3
u
UNLOADED
It is very important that the distinction between stored energy and dissipated
energy be made clear. (Note that some texts use the term absorbed
energy in lieu of stored energy.)
In the first diagram, the system remains elastic and all of the strain energy is
stored. If the bar were released, all of the energy would be recovered.
In the second diagram, the applied deformation is greater than the elastic
deformation and, hence, the system yields. The energy shown in green is
stored, but the energy shown in red is dissipated. If the bar is unloaded, the
stored energy is recovered, but the dissipated energy is lost. This is shown
in Diagrams 3 and 4.
Slide 47
F(t)
Time, T
The discussion will now proceed to general dynamic loading. By general

loading, it is meant that no simple mathematical function defines the entire
loading history.
Slide 48

Solution Techniques
Fourier transform
Duhamel integration
Piecewise exact
Newmark techniques
There are a variety of ways to solve the general loading problem and all are
carried out numerically on the computer. The Fourier transform and
Duhamel integral approaches are not particularly efficient (or easy to
explain) and, hence, these are not covered here. Any text on structural
dynamics will provide the required details.
The piecewise exact method is used primarily in the analysis of linear
systems. The Newmark method is useful for both linear and nonlinear
systems. Only the basic principles underlying of each of these approaches
are presented.
Slide 49
F ( ) = Fo +
Fo
dF
dt
dF
dt
dt
In the piecewise exact method, the loading function is broken into a number
of straight-line segments. In a sense, the name of the method is a misnomer
because the method is not exact when the actual loading is smooth (like a
sine wave) because the straight line load segments are only an
approximation of the actual load. When the actual load is smooth, the
accuracy of the method depends on the level of discretization when defining
the loading function.
For earthquake loads, the load is almost always represented by a recorded
accelerogram, which does consist of straight line segments. (There would
be little use in trying to interpolate the ground motion with smooth curves.)
Hence, for the earthquake problem, the piecewise exact method is truly
exact.
Slide 50

Initial conditions
u o , 0 = 0 u&o , 0 = 0
u1 = u ( )
u&1 = u& ( )
u&&1 = u&&( )
u o ,1 = u ( )
u&0 ,1 = u& ( )
LOOP
u 2 = u ( )
u& 2 = u& ( )
u&&2 = u&&( )
The basic idea of the piecewise exact method is to develop a solution for a
straight line loading segment knowing the initial conditions. Given the initial
conditions and the load segment, the solution at the end of the load step is
determined and this is then used as the initial condition for the next step of
the analysis. The analysis then proceeds step by step until all load
segments have been processed.
Slide 51

Advantages:

Disadvantages:

It should be noted that the piecewise exact method may be used for
nonlinear analysis in certain circumstances. For example, the fast
nonlinear analysis (FNA) method developed by Ed Wilson and used in SAP
2000 utilizes the piecewise exact method. In FNA, the nonlinearities are
right-hand sided, leaving only linear terms in the left-hand side of the
equations of motion.
Slide 52
Newmark Techniques

integration schemes
Derived by:
The Newmark method is one of the most popular methods for solving the
general dynamic loading problem. It is applicable to both linear and
nonlinear systems. It is equally applicable to both SDOF and MDOF
systems.
The Newmark method is described in more detail in the topic on inelastic
behavior of structures.
Slide 53
Newmark Method
Advantages:

Disadvantages:

The advantages and disadvantages of the Newmark method are listed. The
principal advantage is that the method may be applied to inelastic systems.
The method also may be used (without decoupling) for multiple-degree-offreedom systems.
Slide 54
GROUND ACC, g
0.40
0.20
0.00
-0.20
-0.40
0.00
1.00
2.00
3.00
4.00
5.00
6.00
TIME, SECONDS
In an earthquake, no actual force is applied to the building. Instead, the

ground moves back and forth (and up and down) and this movement induces
inertial forces that then deform the structure. It is the displacements in the
structure, relative to the moving base, that impose deformations on the
structure. Through the elastic properties, these deformations cause elastic
forces to develop in the individual members and connections.
Slide 55

0.3
0.2
0.1
40
0
-0.1
-0.2
-0.3
0
10
20
30
40
50
60
Time (sec)
15
10
0.4
30
20
10
0
-10
-20
-30
0
10
20
30
40
50
60
Time (sec)
0
-5
-10
-15
0
10
20
30
40
50
60

Internet.
Time (sec)
Earthquake ground motions usually are imposed through the use of the
ground acceleration record or accelerogram. Some programs (like Abaqus)
may require instead that the ground displacement records be used as input.
Slide 56
u&&g
u&&t
u&&r
GROUND ACC, g
0.40
0.20
0.00
-0.20
-0.40
0.00
1.00
2.00
3.00
4.00
5.00
6.00
TIME, SECONDS

In this slide, it is assumed that the ground acceleration record is used as

input. The total acceleration at the center of mass is equal to the ground
acceleration plus the acceleration of the center of mass relative to the
moving base. The inertial force developed at the center of mass is equal to
the mass times the total acceleration.
The damping force in the system is a function of the velocity of the top of the
structure relative to the moving base. Similarly, the spring force is a function
of the displacement at the top of the structure relative to the moving base.
The equilibrium equation with the zero on the response history spectrum
(RHS) represents the state of the system at any point in time. The zero on
the RHS reflects the fact that there is no applied load.
If that part of the total inertial force due to the ground acceleration is moved
to the right-hand side (the lower equation), all of the forces on the left-hand
side are in terms of the relative acceleration, velocity, and displacement.
This equation is essentially the same as that for an applied load (see Slide 8)
but the effective earthquake force is simply the negative of the mass times
the ground acceleration. The equation is then solved for the response
history of the relative displacement.
Slide 57

mu&&r (t ) + cu& r (t ) + kur (t ) = mu&&g (t )
u&&r (t ) +
k
c
u&r (t ) + ur (t ) = u&&g (t )
m
m
Make substitutions:
c
= 2
m
k
=2
m
Simplified form:
u&&r (t ) + 2 u&r (t ) + 2ur (t ) = u&&g (t )

In preparation for the development of response spectra, it is convenient to

simplify the equation of motion by dividing through by the mass. When the
substitutions are made as indicated, it may be seen that the response is
uniquely defined by the damping ratio, the undamped circular frequency of
vibration, and the ground acceleration record.
Slide 58

u&&r (t ) + 2 u& r (t ) + 2 u r (t ) = u&&g (t )

Damping ratio
This restates the point made in the previous slide. A response spectrum is
created for a particular ground motion and for a structure with a constant
level of damping. The spectrum is obtained by repeatedly solving the
equilibrium equations for structures with varying frequencies of vibration and
then plotting the peak displacement obtained for that frequency versus the
frequency for which the displacement was obtained.
Slide 59

0.4

with given and
0.3
0.2
0.1
0
-0.1
SOLVER
-0.2
-0.3
0
10
20
30
40
Time (sec)
60
6

response
50
Computed response
4
2
0
-2
-4
Peak displacement
-6
0
10
20
30
Time (sec)
40
50
60
The next several slides treat the development of the 5% damped response
spectrum for the 1940 El Centro ground motion record. The solver
indicated in the slide is a routine, such as the Newmark method, that takes
the ground motion record, the damping ratio, and the system frequency as
input and reports as output only the maximum absolute value of the relative
displacement that occurred over the duration of the ground motion. It is
important to note that by taking the absolute value, the sign of the peak
response is lost. The time at which the peak response occurred is also lost
(simply because it is not recorded).
Slide 60

16
12
0
0
10
PERIOD, Seconds
This slide is a restatement of the previous point.
Slide 61

0.08
0.06
Computed response
0.04
0.02
0.00
-0.02
-0.04
-0.06
-0.08
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.10 sec
Umax= 0.0543 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
Here, the first point in the response spectrum is computed. For this and all
subsequent steps, the ground motion record is the same and the damping
ratio is set as 5% critical. Only the frequency of vibration, represented by
period T, is changed.
When T = 0.10 sec (circular frequency = 62.8 radians/sec), the peak
computed relative displacement was 0.0543 inches. The response history
from which the peak was obtained is shown at the top of the slide. This
peak occurred at about 5 sec into the response, but this time is not recorded.
Note the high frequency content of the response.
The first point on the displacement response spectrum is simply the
displacement (0.0543 inches) plotted against the structural period (0.1 sec)
for which the displacement was obtained.
Slide 62

0.40
0.30
Computed response
0.20
0.10
0.00
-0.10
-0.20
-0.30
-0.40
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.20 sec
Umax = 0.254 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
Here the whole procedure is repeated, but the system period is changed to
0.2 sec. The computed displacement history is shown at the top of the slide,
which shows that the peak displacement was 0.254 inches. This peak
occurred at about 2.5 sec into the response but, as before, this time is not
recorded. Note that the response history is somewhat smoother than that in
the previous slide.
The second point on the response spectrum is the peak displacement (0.254
inch) plotted against the system period, which was 0.2 sec.
Slide 63

0.80
0.60
Computed response
0.40
0.20
0.00
-0.20
-0.40
-0.60
-0.80
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.30 sec
Umax = 0.622 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
The third point on the response spectrum is the peak displacement (0.622
inch) plotted against the system period, which was 0.3 sec. Again, the
response is somewhat smoother than before.
Slide 64

1.20
0.90
Computed response
0.60
0.30
0.00
-0.30
-0.60
-0.90
-1.20
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.40 sec
Umax = 0.956 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
The fourth point on the response spectrum is the peak displacement (0.956
inch) plotted against the system period, which was 0.40 sec.
Slide 65

2.40
1.80
Computed response
1.20
0.60
0.00
-0.60
-1.20
-1.80
-2.40
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.50 sec
Umax = 2.02 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
The next point on the response spectrum is the peak displacement (2.02
inches) plotted against the system period, which was 0.50 sec.
Slide 66

3.20
2.40
Computed response
1.60
0.80
0.00
-0.80
-1.60
-2.40
-3.20
0
10
11
12
Time, Seconds
10.00
= 0.05
T = 0.60 sec
Umax= -3.00 in.
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds
The next point on the response spectrum is the peak displacement (3.03
inches) plotted against the system period, which was 0.60 sec. Note that
only the absolute value of the displacement is recorded.
The complete spectrum is obtained by repeating the process for all
remaining periods in the range of 0.7 through 2.0 sec. For this response
spectrum, 2/0.1 or 20 individual points are calculated, requiring 20 full
response history analyses. A real response spectrum would likely be run at
a period resolution of about 0.01 sec, requiring 200 response history
analyses.
Slide 67

Ground Motion
12.00
10.00
8.00
6.00
4.00
2.00
0.00
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Period, Seconds
This is the full 5% damped elastic displacement response spectrum for the
1940 El Centro ground motion. Note that the spectrum was run for periods
up to 4.0 sec. This spectrum was generated using NONLIN.
Note also that the displacement is nearly zero when T is near zero. This is
expected because the relative displacement of a very stiff structure (with T
near zero) should be very small. The displacement then generally increases
with period, although this trend is not consistent. The reductions in
displacement at certain periods indicate that the ground motion has little
energy at these periods. As shown later, a different earthquake will have an
entirely different response spectrum.
Slide 68
Response Spectrum
35.00
5% damping
30.00
25.00
20.00
15.00
10.00
PSV (T ) D
5.00
0.00
0.0
1.0
2.0
3.0
4.0
Period, Seconds
If desired, an elastic (relative) velocity response spectrum could be obtained

in the same way as the displacement spectrum. The only difference in the
procedure would be that the peak velocity computed at each period would be
recorded and plotted.
Instead of doing this, the velocity spectrum is obtained in an approximate
manner by assuming that the displacement response is harmonic and,
hence, that the velocity at each (circular) frequency is equal to the frequency
times the displacement. This comes from the rules for differentiating a
harmonic function.
Because the velocity spectrum so obtained is not exact, it is called the
pseudovelocity response spectrum.
Note that it appears that the pseudovelocity at low (near zero) periods is also
near zero (but not exactly zero).
Slide 69
Response Spectrum
400.0
5% damping
350.0
300.0
250.0
200.0
PSA (T ) 2 D
150.0
100.0
50.0
0.0
0.0
1.0
2.0
3.0
4.0
Period, Seconds
The pseudoacceleration spectrum is obtained from the displacement

spectrum by multiplying by the circular frequencies squared. Note that the
acceleration at a near zero period is not near zero (as was the case for
velocity and displacement). In fact, the pseudoacceleration represents the
total acceleration in the system while the pseudovelocity and the
displacement are relative quantities.
Slide 70

damped structures.
400.0
5% damping
350.0
Peak ground
acceleration
300.0
250.0
200.0
150.0
100.0
50.0
0.0
0.0
1.0
2.0
3.0
4.0
Period, Seconds
For very rigid systems (with near zero periods of vibration), the relative
acceleration will be nearly zero and, hence, the pseudoacceleration, which is
the total acceleration, will be equal to the peak ground acceleration.
Slide 71
u&&g
u&&t
u&&r
GROUND ACC, g
0.40
0.20
0.00
-0.20
-0.40
0.00
1.00
2.00
3.00
4.00
5.00
6.00
TIME, SECONDS

This slide explains why the pseudoacceleration is equal to the total

acceleration. The relative displacement is multiplied by omega to get
pseudovelocity. The pseudovelocity then is multiplied by omega to get the
total acceleration.
Slide 72

300.00
350.00
250.00
200.00
150.00
100.00
50.00
0.00
0.1
1
Period (sec)
Total Acceleration
10
Pseudo-Acceleration
This plot shows total acceleration and pseudoacceleration for a 5% damped

system subject to the El Centro ground motion. Note the similarity in the two
quantities. The difference in the two quantities is only apparent at low
periods.
The difference can be much greater when the damping is set to 10%, 20%,
or 30% critical, and the differences can appear in a wider range of periods.
Slide 73

40
Ve lo city (in/se c)
35
30
25
20
15
10
5
0
0.1
1
Period (sec)
Relative Velocity
10
Pseudo-Velocity
This plot shows relative velocity and pseudovelocity for a 5% damped

system subject to the El Centro ground motion. Here, the differences are
much more apparent than for pseudoacceleration, and the larger differences
occur at the higher periods. The differences will be greater for systems with
larger amounts of damping.
Slide 74

Damping
25.00
0%
5%
10%
20%
20.00
15.00
10.00
5.00
0.00
0.0
1.0
2.0
3.0
4.0
5.0
Period, Seconds
The higher the damping, the lower the relative displacement. At a period of
2 sec, for example, going from zero to 5% damping reduces the
displacement amplitude by a factor of two. While higher damping produces
further decreases in displacement, there is a diminishing return. The %
reduction in displacement by going from 5 to 10% damping is much less that
that for 0 to 5% damping.
Slide 75

Damping
4.00
0%
5%
10%
20%
3.00
2.00
1.00
0.00
Peak ground
acceleration
0.0
1.0
2.0
3.0
4.0
5.0
Period, Seconds
Damping has a similar effect on pseudoacceleration. Note, however, that

the pseudoacceleration at a (near) zero period is the same for all damping
values. This value is always equal to the peak ground acceleration for the
ground motion in question.
Slide 76

of Vibration
Damping is generally effective at all periods (except at T = 0). The reason

for this is that ground motions consist of a large number of harmonics, each
at a different frequency. When a response spectrum analysis is run for a
particular period, there will be a near resonant response at that period.
Damping is most effective at resonance and, hence, damping will be
effective over the full range of periods for which the response spectrum is
generated.
Slide 77

Vibration
Amplitude
4.00
2.00
0.00
-2.00
-4.00
0.0
6.0
12.0 18.0 24.0 30.0 36.0 42.0 48.0 54.0 60.0 66.0 72.0 78.0 84.0 90.0
Time (sec)

accelerogram.
amplitude of 1.0.
To demonstrate the point made in the previous slide, an artificial ground

motion is made up from the sum of five simple harmonics.
Slide 78

1
0.5
0
-0.5
-1
0.0
Amplitud
e
T = 5.0 s
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 4.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 3.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
Time (sec)
Each of the harmonics has an amplitude of 1.0. The first three of the
harmonics with T = 5, 4, and 3 sec are shown.
Slide 79

1
0.5
0
-0.5
-1
0.0
Amplitude
T = 2.0 s
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 1.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
Summation
4.00
2.00
0.00
-2.00
-4.00
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
Time (sec)
The remaining two harmonics (at T = 2 and 1 sec) and the sum are shown.
Slide 80

14.00
0.0% Damping
5.0 % Damping
10.0% Damping
25.0 % Damping
Fourier amplitude
12.00
10.00
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio,
Frequency (Hz)
The Fourier amplitude spectrum of the artificial ground motion is shown at

the left. This spectrum shows the five discrete harmonics that are in the
artificial motion. If the response spectrum is run at intervals of 0.2 sec, there
will be resonant response at each of these frequencies. Damping will be
very effective in reducing the response at each of the frequencies.
Slide 81

Example Structure
12.00
K = 500 k/in
M = 2000/386.4 = 5.18 k-sec2/in
= (K/M)0.5 =9.82 rad/sec
T = 2/ = 0.64 sec
5% critical damping
10.00
W = 2,000 k
8.00
6.00
4.00
2.00
0.00
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Period, Seconds

This is a simple example of the use of an elastic displacement response

spectrum. If the system is assumed to have 5% damping (matching the
spectrum) and the system period is known, the peak displacement may be
easily computed. Note that the sign of the displacement (positive or
negative) and the time that the displacement occurred is not known as this
information was discarded when the spectrum was generated.
Slide 82

400.0
K = 500 k/in
W = 2,000 k
M = 2000/386.4 = 5.18 k-sec2/in
=
(K/M)0.5
=9.82 rad/sec
T = 2/ = 0.64 sec
5% critical damping
Pseudoacceleration, in/sec 2
Example Structure
350.0
300.0
250.0
200.0
150.0
100.0
50.0
0.0
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Period, Seconds

This is a simple example of the use of an elastic pseudoacceleration

response spectrum. If the system is assumed to have 5% damping
(matching the spectrum) and the system period and mass are known, the
peak base shear may be easily computed. Note that the sign of the shear
(positive or negative) and the time that the shear occurred is not known as
this information (related to pseudoacceleration) was discarded when the
spectrum was generated.
Slide 83

1.00
period values, T
0.80
T = 0.64s
0.60
0.40
0.20
0.00
0.00
2.00
4.00
6.00
8.00
10.00
12.00
Another type of spectrum plot is the acceleration-displacement response

spectrum (ADRS), which is also called a demand spectrum. Here,
displacement is plotted on the x-axis and pseudoacceleration is plotted on
the y-axis. Periods of vibration are represented as radial lines.
This kind of spectrum is most commonly used in association with capacity
spectra developed from nonlinear static pushover analysis. A demand
spectrum is also useful in assessing stiffness and damping requirements of
base-isolated systems.
Slide 84

Line of increasing
displacement
100
D=10.0
1.0
0.1
Line of constant
displacement
10
1.0
0.01
D=
1
.01
0.1
0.001
PSV
0.1
0.1
10
100
1000
Circular Frequency
per Second
, Radiand
Circular Frequency
(radians/sec)
Response spectra often are plotted on four-way log paper. This type of
spectrum is often called a tripartite spectrum because the displacement,
pseudovelocity, and pseudoacceleration are all shown on the same plot.
On the plot, pseudovelocity is plotted on the vertical axis. Lines of constant
and logarithmically increasing displacement are generated as shown. The
use of circular frequency on the horizontal axis is rarely used in practice but
is convenient for illustrating the development of the plot.
Slide 85

Line of constant
acceleration
100
PSA=1000
10000
100000
Line of increasing
acceleration
10
100
10000
PSA = PSV
10
100
1000
0.1
0.1
10
100
1000
Circular Frequency
Radiand per
Circular ,
Frequency
Second
(radians/sec)
Lines of constant and logarithmically increasing pseudoacceleration are

obtained in a similar manner.
Slide 86

2
100
DI
S
PL
AC
0 EM
10
EN
s
n/
,i
T,
10
10
N
IO
AT 100
00
ER
EL
C
AC 10
00
in
1.
10
0
0.
0.
ec
10
0.
01
0.
1
00
0.1
0.1
10
100
1000
Circular Frequency
, Radiand per
Circular Frequency
Second
(radians/sec)
This is a completed spectrum for the 5% damped 1940 El Centro earthquake

with maximum acceleration = 0.35g.
Slide 87

Plotted vs Period
n.
,i
ac
em
en
t
D
is
pl
10
.0
g
1
00
0.
0.
00
1
01
0.
0.
01
1.00
10
0.
0.
1
0
1.
1.
0
10.00
0.10
0.01
.0
10
n,
ti o
ra
le
ce
Ac
100.00
0.10
1.00
10.00
PERIOD, Seconds
Response spectra usually are plotted versus structural period or structural

cyclic frequency. This is the same spectrum as shown in the previous slide,
but it is plotted versus period.
Slide 88
Response Spectrum
100.00
ce
is
pl
a
10
.
,g
.
10
n
tio
10.00
1.
0
0
1.
0
0.
1
1
0.
0.
0
1.00
0
0.
a
er
m
en
el
t,
i
n.
c
Ac

0.
0
01
01
0
0.
0.10
0.01
0.10
1.00
PERIOD, Seconds
10.00

Characteristics.
The use of a single earthquake spectrum in structural design is not

recommended for the reasons shown on this slide. The same site
experiencing different earthquakes (or different components of the same
earthquake) often will have dissimilar spectra.
Slide 89

ac
em
en
t
D
is
pl
10
.0

.0
10
n,
ti o
ra
le
ce
Ac
g
0% Damping
5% Damping
10% Damping
20* Damping
10
0.
0.
1
0
1.
1.
0
10
1
00
0.
0.
00
1
0.
01
01
0.
,i
n.
100
0.1
0.01
0.1
1
Period, Seconds
10
Note the significant changes (for any given damping value) in the 1.5 sec
period range.
Slide 90

100.0
El Centro
Loma Prieta
North Ridge
San Fernando
Average
10.0
1.0
0.1
0.01
0.10
1.00
10.00
Period, seconds
The spectra are scaled to 0.4 g with 5% damping. Note the differences.
Slide 91

ASCE 7 spectrum
Because real ground motion spectra are difficult to work with in a design
office, a variety of empirical spectra have been generated. One of the
earliest of these empirical spectra was developed by Nathan Newmark. The
next several slides describe this in detail.
The spectrum used by ASCE 7-05 is simpler than the Newmark spectrum,
but explanation of the background of the ASCE 7 spectrum is more difficult.
Certain key aspects of the ASCE 7 spectrum are presented in the topic on
seismic load analysis.
Slide 92

100
0% Dam ping
5% Dam ping
Displacem ent (in)
10% Dam ping
10
max u&
Observations
v&& max v&&

at short T
v0
max u
max u&&
v max v
v&& 0
0.1
0.01
0.1
Period1 (sec)
10
at long T
100
The Newmark spectrum is based on the following observations:
The pseudoacceleration at very low periods is exactly equal to the peak

ground acceleration.
The relative displacement at very long periods is exactly equal to the

peak ground displacement.
At intermediate periods, the displacement, pseudovelocity, and

pseudoacceleration are equal to the ground values times some empirical
constant.
Slide 93
Total acceleration
Zero
Ground acceleration
For very low period (high frequency) buildings, the maximum relative
displacement will be zero. The maximum acceleration will approach the
ground acceleration.
Slide 94
Total acceleration
Ground displacement
Zero
For very high period (low frequency) buildings, the maximum relative
displacement will be equal to the maximum ground displacement. The
maximum total acceleration will approach zero.
Slide 95

100
n.
,i
ac
em
en
t
D
is
pl
.0
10
10
.0
n,
ti o
ra
le
ce
Ac
10
1.
0
0
1.
0% Damping
5% Damping
10% Damping
20* Damping
4.25 in.
0.
01
0.
1
01
0.
10
0.
12.7 in/s
0.35g
0.
00
1
1
00
0.
0.1
0.01
Ground Maxima
0.1
1
Period, Seconds
10
The yellow line shows the maximum recorded ground displacement, velocity,
and acceleration from the 1940 El Centro earthquake. These lines clearly
form a lower bound to the elastic response spectra. Note how the building
response displacements, velocities, and accelerations are amplifications of
the ground values. Note also how the amplifications decrease with
increased damping.
Slide 96

Damping
% Critical
.05
1
2
3
5
7
10
20
One Sigma (84.1%)

av
ad
aa
5.10 3.84 3.04
4.38 3.38 2.73
3.66 2.92 2.42
3.24 2.64 2.24
2.71 2.30 2.01
2.36 2.08 1.85
1.99 1.84 1.69
1.26 1.37 1.38
aa
3.68
3.21
2.74
2.46
2.12
1.89
1.64
1.17
Median (50%)
av
ad
2.59 2.01
2.31 1.82
2.03 1.63
1.86 1.52
1.65 1.39
1.51 1.29
1.37 1.20
1.08 1.01
Newmark has developed a series of amplification factors to be used in the

development of design spectra. These are based on the average of dozens
of spectra recorded on firm soil sites for the western United States. Values
are shown for the median and median plus one standard deviation.
Slide 97

1) Draw the lines
v&& , v& , v
2) Draw line max v
&&
3
4
from Tb to Tc
6
3) Draw line V
from Tc to Td
4) Draw line D
from Td to Te
max v&
max v

from Ta to Tb
Ta
Tb
Tc
Td
Te Tf

from Te to Tf
These are the steps in the development of the Newmark spectrum. Note
that actual values are not present.
Slide 98
Spectral Response
Acceleration, Sa
ASCE 7
Short
period
acceleration
2
SDS
SD1
Long period
acceleration
SDS
T + 0.4 SDS
T0
Sa = 0.6
Sa = SDS
3
4
SD1
T
TLSD1
Sa = 2
T
Sa =
TS
Period, T
T = 1.0
TL

This plot shows the basic relationships used for the ASCE 7 spectrum. Note
that the vertical axis is pseudoacceleration. The spectrum is derived from a
series of maps giving spectral acceleration values for short period (T = 0.2
sec) or long period (T = 1 sec) buildings. Note that the part of the spectrum
to the right of TL (Curve 4) was introduced in the 2003 NEHRP
Recommended Provisions and in ASCE 7-05.
The maps are based on very stiff soils. For design purposes, the
acceleration spectra is not reduced to the ground acceleration at low periods
(Line 1 on the plot). Damping is assumed to be 5% critical.
Slide 99

hazard analysis.
This slide notes that the ASCE 7 spectrum is a uniform hazard spectrum.
This concept is covered in detail in the topic on seismic hazard analysis.
The main purpose of this side is a transition into the hazards topic.
Slide 100
(SDOF) Systems
Structural Dynamics

SDOF Dynamics 3 - 1
SDOF Dynamics 3 - 2

coordinates.
F(t)
F ( t ), u ( t )
Mass

of vibration.
Damping
Stiffness
u(t)
damping in system.

t

input.
SDOF Dynamics 3 - 3
SDOF Dynamics 3 - 4

Applied Force, kips
40
f I (t )
F (t )
0
-40
0.00
0 .5 f S ( t )
fD (t )
0.20
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
0.40
0.60
0.80
1.00
Displacement, in
0.50
0 .5 f S ( t )
0.00
-0.50
0.00
0.20
Velocity, in/sec
15.00
0.00
F (t ) f I (t ) fD (t ) fS (t ) = 0
fI (t ) + fD (t ) + fS (t ) = F (t )
-15.00
0.00
0.20
400.00
0.00
-400.00
0.00
0.20
Time, sec
FEMA 451B Topic 3 Handouts
SDOF Dynamics 3 - 5
SDOF Dynamics 3 - 6
Structural Dynamics 1

Spring Force, kips
4.00
50.00
15.00
2.00
25.00
0.00
0.00
0.00
-15.00
-2.00
-25.00
-30.00
-0.60
-4.00
-20.00
0.00
0.30
0.60
Slope = k
= 50 kip/in
f S ( t ) = k u( t )
-10.00
0.00
10.00
f I (t )
F (t )
Damping Force, Kips
30.00
-0.30
20.00
0 .5 f S ( t )
-50.00
-500
-250
250
fD (t )
0 .5 f S ( t )
500
Velocity, In/sec
fI (t ) + fD (t ) + fS (t ) = F (t )
Slope = m
= 0.130 kip-sec2/in
Slope = c
= 0.254 kip-sec/in
m u&&( t ) + c u& ( t ) + k u ( t ) = F ( t )
f I ( t ) = m u&&( t )
f D ( t ) = c u&( t )
SDOF Dynamics 3 - 7
SDOF Dynamics 3 - 8
Damping
M
Damping Force
Internal Force
Mass
1.0
1.0
Acceleration
Velocity


SDOF Dynamics 3 - 9
AREA =
ENERGY
DISSIPATED
Spring Force
K
1.0
Displacement
Displacement

Stiffness
Damping Force
Damping


Spring Force
Stiffness
Initial conditions:
Assume:
u& 0
u (t ) =
u& 0
Displacement

1.0
1.5
2.0
Cyclic Frequency
(cycles/sec, Hertz)
f =
2
Period of Vibration
(sec/cycle)
1
2
T =
=

Ta = 0.028h
0.8
n
Ct =
Ct =
Ct =
Ct =
0.028, x = 0.8
0.016, x = 0.9
0.030, x = 0.75
0.020, x = 0.75

Time, seconds
k
m
Ta = Ct hnx
0.5
Circular Frequency
(radians/sec)
sin( t ) + u 0 cos( t )
(ASCE 7-05)
T = 0.5 sec
u0
0.0
k
m
1.0
3
2
1
0
-1
-2
-3
B = u0

u ( t ) = A sin( t ) + B cos( t )
A=
Solution:
AREA =
ENERGY
DISSIPATED
u& 0
m u&&( t ) + k u( t ) = 0
u& 0 u 0
Equation of motion:
T = 0.1N
a


T = 1.9 sec
T = 1.1 sec
T = 0.15 sec

T = 1.3 sec
T = 0.8 sec
T = 0.1 sec
Gravity dam
Suspension bridge
T = 0.2 sec
T = 20 sec

in ELF Analysis?
T = Ta Cu Tcomputed

SD1
> 0.40g
0.30g
0.20g
0.15g
< 0.1g

Cu
1.4
1.4
1.5
1.6
1.7
if Tc > CuTa

if Tc < Ta
m u&&( t ) + c u& ( t ) + k u ( t ) = 0
u& 0
st
Assume: u ( t ) = e

Displacement, in
u& + u 0
u ( t ) = e t u 0 cos( D t ) + 0
sin( D t )
D
Time, sec
D = 1 2
Damping Force, Kips

50.00
15.00
2.00
25.00
0.00
0.00
0.00
-15.00
-2.00
-25.00

4.00
3
2
1
0
-1
-2
-3
0% Damping
10% Damping
20% Damping
0.0
-30.00
-0.60
-0.30
0.00
0.30
0.60
-4.00
-20.00
-10.00
0.00
10.00
20.00
-50.00
-500
Velocity, In/sec
Spring Force, kips
Solution:
30.00
c
c
=
=
cc
2m
Equation of motion:
c
c
=
cc
2m
use T = Ta
use T = CuTa
0.5
1.0
1.5
2.0
Time, seconds
-250
250
500
Welded steel frame

Bolted steel frame
= 0.010
= 0.020

= 0.015
= 0.020
= 0.035

= 0.100
= 0.150

= 0.050
= 0.075
= 0.250
Inherent damping

Added damping
For all
damping values
u1
u2
u3

Equation of motion:
1 2
u0 e
-0.5
-1
0.00
For very low

damping values
0.50
1.00
1.50
2.00
2.50
3.00
u1 u 2
2 u2
Time, Seconds
m u&&( t ) + k u ( t ) = p 0 sin( t )
T =
Force, Kips
Amplitude
0.5

u
ln 1 =
u2

1

150
100
50
0
-50
-100
-150
0.00
= 0.25 sec
po=100 kips
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
Time, Seconds

&&( t ) + k u ( t ) = p 0 s in ( t )
Define

Particular solution: u ( t ) = C s i n ( t )
u (t ) =
u( t ) =
1
p0
sin( t ) sin( t )
2
k 1 ( / )
Loading frequency
Dynamic magnifier
Solution:
Transient response
1
p0
(sin( t ) sin( t ))
k 1 2
Static displacement
Steady state
response

200
100
0
-100
-200
0.00
Loading (kips)
0.25
150
100
50
0
-5 0
-1 0 0
-1 5 0
0 .0 0
0.50
0.75
1.00
1.25
1.50
1.75
2.00
0.50
0.75
1.00
1.25
1.50
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
1.75
0 .25
0 .5 0
0 .7 5
1 .0 0
1.2 5
1.5 0
1.75
2.00
0.2 5
0.50
0 .75
1.0 0
1.25
1 .50
1.7 5
2.00
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
-5 00
0.00
2.00
500
250
0
Transient
response (in.)
0.25
10
5
0
-5
-10
0.00
0.50
0.75
1.00
1.25
1.50
1.75
-250
-500
0 .00
2.00
80
40
0.25
0.50
0.75
1.00
1.25
1.50
1.75
0
-40
Total response
(in.)
-80
0.00
2.00
Time, seconds
T ime, seconds

80
150
100
50
0
-5 0
-1 0 0
-1 5 0
0 .0 0
Loading (kips)
2 uS
40
Displacement, in.
0 .2 5
2 50
-2 50
0.25
Total response
(in.)
= 0.99 u S = 5.0 in.
5 00
Steady state
response (in.)
10
5
0
-5
-10
0.00
Transient
response (in.)

Loading
(kips)
10
5
0
-5
-10
0.00
Steady state
response (in.)
= 0.5 uS = 5.0 in.
= 1.01
u S = 5.0 in.
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1.2 5
1.5 0
1 .7 5
500
Steady state
response (in.)
250
0
-2 5 0
-5 0 0
0 .0 0
500
250
Transient
response (in.)
-40
0
-2 5 0
-5 0 0
0 .0 0
Linear envelope
0.25
0.50
0.75
1.00
1.25
1.50
1.75
40
Total response
(in.)
2.00
0
p
-80
0.00
80
Time, seconds
-4 0
-8 0
0 .0 0
2 .0 0
T im e , s ec o n d s
Loading (kips)
150
100
50
0
-50
-100
-150
0.00
= 2.0 uS = 5.0 in.
=8 rad / sec
0.25
0.50
0.75
1.00
1.25
1.50
1.75
2.00
3
0
-3
-6
0 .0 0
0 .2 5
0 .5 0
0.7 5
1 .00
1 .2 5
1 .5 0
1 .7 5
2 .0 0
6
3
Transient response
(in.)
0
-3
-6
0 .0 0
(Signs Retained)
12.00
Steady state
response (in.)
= 4 rad / sec
8.00
In phase
4.00
Resonance
0.00
-4.00
-8.00
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .75
2.00
Total response
(in.)
-3
-12.00
0.00
-6
0 .0 0
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio
0 .2 5
0 .5 0
0 .7 5
1 .0 0
1 .2 5
1 .5 0
1 .7 5
2 .0 0
(Absolute Values)
Equation of motion:
12.00
m u&&( t ) + cu&( t ) + k u ( t ) = p 0 sin( t )
10.00
T =
8.00
= 0.25 sec
6.00
Slowly
loaded
4.00
Rapidly
loaded
2.00
1.00
0.00
0.00
Force, Kips
Resonance
0.50
1.00
1.50
2.00
2.50
150
100
50
0
-50
-100
-150
0.00
po=100 kips
0.25
0.50
0.75
3.00
1.00
1.25
1.50
1.75
2.00
Time, Seconds
Frequency Ratio
Equation of motion:
m u&&( t ) + cu& ( t ) + k u ( t ) = p 0 sin( t )

u ( t ) = e t [ A sin( D t ) + B cos( D t ) ] +

u ( t ) = C sin( t ) + D cos( t )
C sin( t ) + D cos( t )
u(t ) = e
[ A sin( D t ) + B cos( D t ) ]
Solution:

c
2 m
D = 1 2
C=

BETA=1 (Resonance)
Beta=0.5
Beta=2.0
2
po
k (1 2 ) 2 + (2 ) 2

40
D=
50
50
40
30
20
10
0
-10
-20
1
Static
2
30
20
10
0
-10
-20
-30
-40
-30
-50
0.00
-40
-50
0.00
1 2
po
2 2
k (1 ) + (2 ) 2
1.00
2.00
3.00
4.00
5.00
Time, Seconds
1.00
2.00
3.00
4.00
5.00
Time, Seconds
14.00
Resonance
Effects of Damping
0.0% Damping
5.0 % Damping
10.0% Damping
25.0 % Damping
12.00
200
150
100
50
0
-50
-100
10.00
8.00
6.00
2.00
1.00
2.00
3.00
4.00
(1 ) + ( 2 ) 2
Slowly
loaded
Rapidly
loaded
5.00
Time, Seconds
0% Damping
1
2 2
4.00
-150
-200
0.00
RD =
0.00
0.00
0.50
%5 Damping
1.00
1.50
2.00
2.50
3.00
Frequency Ratio,





over time.


resonance.

F
Energy
Stored
Energy
Dissipated
2
1
F(t)
LOADING
YIELDING
Energy
Recovered
2
u
Total
Energy
Dissipated
UNLOADING
Time, T
UNLOADED

Solution Techniques
Fourier transform
Duhamel integration
Piecewise exact
Newmark techniques
F ( ) = Fo +
Fo
dF
dt
dF
dt
dt

Initial conditions
Advantages:
u&1 = u& ( )
u&&1 = u&&( )

Disadvantages:
u o ,1 = u ( )
u o , 0 = 0 u&o , 0 = 0
u1 = u ( )
u& 0,1 = u& ( )

LOOP
u 2 = u ( )
u& 2 = u& ( )
u&&2 = u&&( )

Newmark Method
Newmark Techniques
Advantages:
integration schemes
Derived by:

Disadvantages:

0.3
0.2
0.1
40
0
-0.1
-0.2
-0.3
0
10
20
30
40
50
60
Time (sec)
15
GROUND ACC, g
0.40
0.20
0.00
-0.20
-0.40
0.00
10
20
10
0
-10
-20
-30
0
1.00
2.00
3.00
4.00
5.00
20
30
40
50
60
Time (sec)
-5

Internet.
-10
6.00
10
-15
10
20
30
40
50
60
Time (sec)
u&&t
30
TIME, SECONDS
u&&g
0.4
mu&&r (t ) + cu& r (t ) + kur (t ) = mu&&g (t )
u&&r
u&&r (t ) +
GROUND ACC, g
0.40
0.20
k
c
u&r (t ) + ur (t ) = u&&g (t )
m
m
0.00
Make substitutions:
-0.20
-0.40
0.00
1.00
2.00
3.00
4.00
5.00
m[&&
ug ( t ) + u&&r ( t )] + c u&r ( t ) + k ur ( t ) = 0
Simplified form:
u&&r (t ) + 2 u&r (t ) + 2u r (t ) = u&&g (t )



0.4
k
=2
m
c
= 2
m
6.00
TIME, SECONDS

with given and
0.3
0.2
0.1
0
-0.1
SOLVER
-0.2
-0.3
10
20
30
Time (sec)

response
Damping ratio
40
50
60
6
u&&r (t ) + 2 u& r (t ) + 2 u r (t ) = u&&g (t )
Computed response
4
2
0
-2
-4
Peak displacement
-6
0
10
20
30
40
50
60
Time (sec)


0.08

0.06
Computed response
0.04
0.02
0.00
-0.02
-0.04
-0.06
-0.08
10
11
12
Time, Seconds

10.00

8.00
= 0.05
12
16
T = 0.10 sec
8
Umax= 0.0543 in.
4.00
2.00
0.00
0.00
0
0
6.00
10
0.50
1.00
1.50
2.00
Period, Seconds
PERIOD, Seconds


0.80
0.30
0.40
Computed response
0.20
0.10
0.00
-0.10
-0.20
-0.30
0.60
Computed response
0.40
0.20
0.00
-0.20
-0.40
-0.60
-0.80
-0.40
10
11
12
10
11
12
8.00
Umax = 0.254 in.

8.00
T = 0.20 sec
10.00
10.00
= 0.05
Time, Seconds
Time, Seconds
= 0.05
6.00
T = 0.30 sec
4.00
Umax = 0.622 in.
2.00
0.00
0.00
0.50
1.00
1.50
6.00
4.00
2.00
0.00
0.00
2.00
0.50
1.00
1.50
2.00
Period, Seconds
Period, Seconds


1.20
2.40
0.90
Computed response
0.60
0.30
0.00
-0.30
-0.60
-0.90
-1.20
1.80
Computed response
1.20
0.60
0.00
-0.60
-1.20
-1.80
-2.40
10
11
12
Time, Seconds
10
11
12
6.00
= 0.05
4.00
T = 0.50 sec
2.00
Umax = 2.02 in.
0.00
0.00
8.00
0.50
1.00
1.50
2.00
Period, Seconds
Umax = 0.956 in.
T = 0.40 sec
10.00
8.00
= 0.05
Time, Seconds
10.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
Period, Seconds

Ground Motion

3.20
2.40
12.00
Computed response
1.60
0.80
0.00
10.00
-0.80
-1.60
-2.40
0
10
11
-3.20
12
Time, Seconds
10.00
= 0.05
T = 0.60 sec
Umax= -3.00 in.
8.00
8.00
6.00
4.00
6.00
2.00
4.00
2.00
0.00
0.0
0.00
0.00
0.50
1.00
1.50
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Period, Seconds
2.00
Period, Seconds
Response Spectrum
Response Spectrum
400.0
35.00
30.00
5% damping
350.0
5% damping
25.00
20.00
15.00
300.0
250.0
200.0
PSA (T ) 2 D
150.0
100.0
10.00
PSV (T ) D
5.00
50.0
0.0
0.0
0.00
0.0
1.0
2.0
3.0
1.0
2.0
3.0
4.0
Period, Seconds
4.0
Period, Seconds

damped structures.
u&&g
u&&t
u&&r
400.0
0.40
GROUND ACC, g
5% damping
350.0
Peak ground
acceleration
300.0
250.0
0.20
0.00
-0.20
-0.40
0.00
200.0
1.00
2.00
3.00
4.00
5.00
6.00
TIME, SECONDS
150.0
m[&&
ug ( t ) + u&&r ( t )] + c u&r ( t ) + k ur ( t ) = 0
100.0
50.0
0.0
0.0
1.0
2.0
3.0
4.0
Period, Seconds


40
35
300.00
Ve locity (in/se c)
350.00
250.00
200.00
150.00
100.00
30
25
20
15
10
50.00
0.00
0.1
1
Period (sec)
Total Acceleration
10
0.1
1
Period (sec)
Pseudo-Acceleration
Relative Velocity
Pseudo-Velocity

Damping
4.00
0%
5%
10%
20%
20.00

Damping
25.00
10
15.00
10.00
5.00
0%
5%
10%
20%
3.00
2.00
1.00
0.00
0.00
0.0
1.0
2.0
3.0
4.0
5.0
Period, Seconds
Peak ground
acceleration

2.0
3.0
4.0
5.0
Period, Seconds

Vibration
Amplitude

1.0

of Vibration
0.0
4.00
2.00
0.00
-2.00
-4.00
0.0
6.0
12.0 18.0 24.0 30.0 36.0 42.0 48.0 54.0 60.0 66.0 72.0 78.0 84.0 90.0
Time (sec)

accelerogram.
amplitude of 1.0.
Amplitud
e
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 2.0 s
0.0
T = 4.0 s
0.0
1
0.5
0
-0.5
-1
T = 5.0 s
Amplitude
1
0.5
0
-0.5
-1
6.0
12.0
18.0
24.0
30.0
36.0
42.0
1
0.5
0
-0.5
-1
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
T = 1.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
Summation
1
0.5
0
-0.5
-1
T = 3.0 s
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
54.0
60.0
66.0
72.0
78.0
84.0
90.0
4.00
2.00
0.00
-2.00
-4.00
0.0
6.0
12.0
18.0
24.0
30.0
36.0
42.0
48.0
72.0
78.0
84.0
90.0
12.00
K = 500 k/in
10.00
10.00
W = 2,000 k
M = 2000/386.4 = 5.18 k-sec2/in
8.00
6.00
= (K/M)0.5 =9.82 rad/sec
4.00
T = 2/ = 0.64 sec
2.00
5% critical damping
0.0% Damping
5.0 % Damping
10.0% Damping
25.0 % Damping
12.00
66.0

Example Structure
14.00
Fourier amplitude
60.0

0.00
0.00
54.0
Time (sec)
Time (sec)
8.00
6.00
4.00
2.00
0.50
1.00
1.50
2.00
2.50
3.00
Frequency Ratio,
0.00
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Period, Seconds
Frequency (Hz)

M = 2000/386.4 = 5.18 k-sec2/in

= (K/M)0.5 =9.82 rad/sec
T = 2/ = 0.64 sec
5% critical damping
1.00
350.0
300.0

period values, T
0.80
W = 2,000 k
Pseudoacceleration, in/sec 2
K = 500 k/in
400.0
Example Structure
250.0
200.0
150.0
100.0
50.0
0.0
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
T = 0.64s
0.60
0.40
0.20
Period, Seconds
0.00
0.00
2.00
4.00
6.00
8.00
10.00
12.00
Line of increasing
displacement
Line of constant
acceleration
100
1.0
100
0.1
10
1.0
0.01
D=
1
.01
0.1
10000
PSA=1000
Line of constant
displacement
D=10.0
0.001
PSV
100000
Line of increasing
acceleration
10
100
PSA = PSV
10000
10
1000
100
0.1
0.1
0.1
10
100
0.1
1000
10
100
1000
Circular Frequency
, Radiand per
Circular Frequency
Second
Circular Frequency
per Second
, Radiand
Circular Frequency
(radians/sec)
(radians/sec)
Plotted vs Period
1000
0.
1
0.10
0.01
0.
01
0.
00
1
10
.0
1.
0
0.
1
1.
0
0.
1
.
, in
en
t
ce
m
Di
sp
la
10
.0
01
0.
0.
01
0.
00
1
1
00
0.
0% Damping
5% Damping
10% Damping
20* Damping
10
0.
10
0.
1
00
0.

Characteristics.
10
0
1.
0
1.
01
0.
PERIOD, Seconds
10.00
,g
.0
10
n
io
1.00
n,
tio
t
ra
0.10

.0
10
ra
le
ce
Ac
le
ce
Ac

1.00
0.10
0.01

100
10.00
10.00
Response Spectrum
100.00
1.00
PERIOD, Seconds
is
pl
ac
em
en
t,
in
.
0.10
100
(radians/sec)
0.
00
1
10
Circular Frequency
, Radiand per
Circular Frequency
Second
0.
01
is
pl
ac
em
en
t,
in
.
1.
0
0.1
0.1
10
.0
1
00
0.
1
00
0.
1.00
01
0.
01
0.
0.
1
10
0.
10
1
0.
10.00
0
1.
10
0
0
1.
.0
10
in
10
n
,i
N
IO
AT 100
00
ER
EL
C
AC 10
00
n,
tio
EN
T,
10
100.00
DI
SP
LA
C
0 EM
10
ec
/s
ra
le
ce
Ac
100
0.1
0.01
0.1
1
Period, Seconds
10


100.0
El Centro
Loma Prieta
North Ridge
San Fernando
Average
10.0
1.0
0.1
0.01
0.10
1.00
ASCE 7 spectrum
10.00
Period, seconds

100
0% Dam ping
5% Dam ping
Displacem ent (in)
10% Dam ping
10
max u&
Observations
v&& max v&&

at short T
v0
max u
max u&&
v max v
v&& 0
0.1
0.01
0.1
Period1 (sec)
10
at long T
Zero
Total acceleration
100
Ground acceleration
100
is
pl
ac
em
en
t,
in
.
D
10
.0
.0
10
n,
tio
1.
0
0
1.
0% Damping
5% Damping
10% Damping
20* Damping
4.25 in.
0.
01
0.
1
10
0.
12.7 in/s
01
0.
ra
le
ce
Ac
10
0.35g
Zero
0.1
0.01
Ground Maxima
0.
00
1
Ground displacement
1
00
0.
Total acceleration
0.1
1
Period, Seconds
10

Damping
% Critical
.05
1
2
3
5
7
10
20
One Sigma (84.1%)

aa
av
ad
5.10 3.84 3.04
4.38 3.38 2.73
3.66 2.92 2.42
3.24 2.64 2.24
2.71 2.30 2.01
2.36 2.08 1.85
1.99 1.84 1.69
1.26 1.37 1.38
aa
3.68
3.21
2.74
2.46
2.12
1.89
1.64
1.17
Median (50%)
av
ad
2.59 2.01
2.31 1.82
2.03 1.63
1.86 1.52
1.65 1.39
1.51 1.29
1.37 1.20
1.08 1.01

1) Draw the lines
2) Draw line A
from Tb to Tc
3) Draw line V
from Tc to Td
4) Draw line D
from Td to Te
v&& , v& , v
max v&&
g
max v&
max v

from Ta to Tb
Ta
Tb
Tc
Td
Te Tf

from Te to Tf
ASCE 7
Spectral Response
Acceleration, Sa

Short
period
acceleration
2
SDS
SD1
Long period
acceleration
SDS
Sa = 0.6
T + 0.4 SDS
T0
Sa = SDS
3
4
SD1
T
TS
Sa = L 2D1
T
Sa =

hazard analysis.
TS
Period, T
T = 1.0
TL


Topic 3 (Fema) Sdof

Uploaded by

Document Informationclick to expand document information

Copyright:

Available Formats

Topic 3 (Fema) Sdof

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Topic 3 (Fema) Sdof

Uploaded by

Copyright:

Available Formats

Structural Dynamics of

Linear Elastic Single-Degree-of-Freedom

Instructional Material Complementing FEMA 451, Design Examples

Equations of motion for SDOF structures

Instructional Material Complementing FEMA 451, Design Examples

Importance in Relation to ASCE 7-05

Ground motion maps provide ground

Equivalent lateral force procedure gives base

Response spectrum is based on 5% critical

Modal superposition analysis uses design

Idealized SDOF Structure

Instructional Material Complementing FEMA 451, Design Examples

Equation of Dynamic Equilibrium

Observed Response of Linear SDOF

Observed Response of Linear SDOF

Inertial Force, kips

Damping Force, Kips

Instructional Material Complementing FEMA 451, Design Examples

Equation of Dynamic Equilibrium

Properties of Structural Mass

Includes all dead weight of structure

Properties of Structural Damping

In absence of dampers, is called inherent damping

Properties of Structural Damping (2)

Damping vs displacement response is

Instructional Material Complementing FEMA 451, Design Examples

Properties of Structural Stiffness

Includes all structural members

Properties of Structural Stiffness (2)

Is almost always nonlinear in real seismic response

Undamped Free Vibration

Instructional Material Complementing FEMA 451, Design Examples

Undamped Free Vibration (2)

Instructional Material Complementing FEMA 451, Design Examples

Approximate Periods of Vibration

for steel moment frames

Note: This applies ONLY to building structures!

For moment frames < 12 stories in height, minimum

Empirical Data for Determination

Instructional Material Complementing FEMA 451, Design Examples

Periods of Vibration of Common Structures

20-story braced frame

Instructional Material Complementing FEMA 451, Design Examples

Adjustment Factor on Approximate Period

Applicable ONLY if Tcomputed comes from a properly

Which Period of Vibration to Use

if Ta < Tc < TuCa use T = Tc

Instructional Material Complementing FEMA 451, Design Examples

Damped Free Vibration

Instructional Material Complementing FEMA 451, Design Examples

cc is the critical damping constant.

is expressed as a ratio (0.0 < < 1.0) in computations.

Response of Critically Damped System, =1.0 or 100% critical

Inertial Force, kips

Damping Force, Kips

Instructional Material Complementing FEMA 451, Design Examples

Damped Free Vibration (2)