Module 5
Module 5
1
“If a civil engineer is to acquire fruitful experience in a brief span of time,
expose him to the concepts of earthquake engineering, no matter if he is
later not to work in earthquake country.”
2
Seismic Response of Linear Elastic SDOF systems
In this lecture, we will study the seismic response of linear SDOF systems subjected
to earthquake excitations.
By definition, all linear systems are elastic systems. They are also referred to as
linearly elastic systems in these lectures to emphasize both properties.
Permanent deformation due to damages
as a result of energy dissipation
Linear elastic system Elastic-perfectly plastic
fs system (Elasto plastic system)
fs
u u
No energy is dissipated by system and system returns Energy is dissipated by system when system is stressed
to its original position on deloading without permanent beyond elastic range. System returns to its original
deformation position on deloading but with permanent deformation 3
Strong motion record
Earthquake shaking capable of causing damage to structures is called Strong Ground Motion*
1 2 3 4 5
* There is no specific PGA to designate an earthquake as strong ground motion earthquake. Usually, an
earthquake with PGA≈0.05g is considered strong ground motion 4
Strong ground motions Accelerograms
g
u
5
Accelerogram used in lectures
g
Ground acceleration, u
Ground velocity, u g
Ground displacement, u g
𝑐 𝑘
⇒ 𝑢ሷ + 𝑢ሶ + 𝑢 = −𝑢ሷ 𝑔 (𝑡)
𝑚 𝑚
2𝜁𝑚𝑢𝜔
ሶ 𝑛
⇒ 𝑢ሷ + + 𝜔𝑛 2 𝑢 = −𝑢ሷ 𝑔 (𝑡)
𝑚
⇒ 𝑢ሷ + 2𝜁𝜔𝑛 𝑢ሶ + 𝜔𝑛 2 𝑢 = −𝑢ሷ 𝑔 (𝑡)
The time variation of ground displacement, from the given time variation of ground acceleration, can be
determined by using any appropriate time stepping numerical method.
Closer the time interval, more accurate will be solution. Typically, the time interval is chosen to be 1/100 to 1/50
of a second, requiring 1500 to 3000 ordinates to describe the round motion of above given El- Centro ,1940, ground
acceleration record having a duration of 30 sec.
Response quantities
𝑡
Sometime, calculating the total acceleration, 𝑢ሷ 𝑜 , in the mass due to base excitations would be
needed if the structure is supporting sensitive equipment. 𝑢ሷ 𝑜 𝑡 can be reduced by using a proper
vibration isolation technique discussed in Module 4
Response History
Response history is a graph between
Acceleration Response history
structural response to ground acceleration
history (Accelerogram) and time
If the response quantity of interest is
Velocity Response history
displacement, u , the Response History is called
Displacement or deformation response history
(a)
(b)
Response histories of building recorded at roof level (Figure a) to ground motion recorded (Figure b) at the
Holiday Inn building’s ground level (during the M6.4 1971 San Fernando EQ)
Displacement Response, u, due to ground acceleration, 𝑢ሷ 𝑔
𝑢, in
𝑢, in
Relative displacement, u, at the
mass level of the SDOF system Relative displacement, u, at mass level
𝑢 = 𝑢𝑡 − 𝑢𝑔
𝑢 = 𝑢𝑡 − 𝑢𝑔 m
k c
𝑢ሷ 𝑔
Base 𝑇𝑛 = 0.5𝑠𝑒𝑐, 𝜁 = 2%
𝑇𝑛 = 0.5𝑠𝑒𝑐, 𝜁 = 2%
𝑢ሷ 𝑔
𝑢ሷ 𝑔 , g
𝑢ሷ 𝑔𝑜 = 0.319𝑔 𝑢ሷ 𝑔 , g
El Centro,1940, ground acceleration applied 𝑢ሷ 𝑔𝑜 = 0.319𝑔
at the base of SDOF system 11
Dynamic Amplification
Ratio of the Response of the structure in the mass at roof level to the motion of
the ground itself is known as Dynamic Amplification. DA ≈ 350/225=1.56
DA ≈ 65/30=2.17
DA ≈ 25/15=2.17
Holiday Inn ground and building roof motion during the M6.4 1971 San Fernando earthquake: (a) north-south
ground acceleration, velocity, and displacement and (b) north-south roof acceleration, velocity, and displacement. 12
Influence of Tn and ζ on Peak displacement, uo , in a linear elastic
SDOF system
𝑢ሷ + 2𝜁𝜔𝑛 𝑢ሶ + 𝜔𝑛 2 𝑢 = −𝑢ሷ 𝑔 (𝑡)
The above given equation indicates that Structural displacement, u, depends upon Tn and ζ
Thus any two systems having the same values of Tn and ζ will have the same deformation
response u(t) even though one system may be more massive than the other or one may be
stiffer than the other
m1 m2= 8m1 𝑘2 = 3𝐸𝐼/ℎ23 = 3𝐸𝐼/ ℎ1 /2 3
𝐸𝐼
𝑘2 = 24 = 8𝑘1
𝑘1 = 3𝐸𝐼/ℎ13 ℎ1 3
h1 𝑚2 8𝑚1
h2=h1/2 𝑇𝑛2 = 2𝜋 ൗ𝑘 = 2𝜋
𝑚1 2 8𝑘1
𝑇𝑛1 = 2𝜋 ൗ𝑘
1
𝑚1
EI= Constt 𝑇𝑛2 = 2𝜋 ൗ𝑘 = 𝑇𝑛1
1
13
Effect of Tn on Deformation Response history
𝑢ሷ 𝑔 , 𝑔
𝑢ሷ 𝑔𝑜 = 0.319𝑔
ug , g
go = 0.319g
u
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).
16
SDOF or MDOF system ?
The concept of response spectrum was introduced by M.A. Biot in his Doctoral thesis in
1932 to characterize the effects of earthquakes on structures.
G.W. Housner popularized its use in earthquake engineering.
The importance of the concept is the fact that, for design purposes, we are only interested
in the maximum values of a structure’s response to an earthquake.
A plot of the peak values of a Response quantity as a function of the natural vibration
period Tn of the system, or a related parameter such as circular frequency ωn or cyclic
frequency fn , is called the Response Spectrum for that quantity.
18
Deformation Response Spectrum
The deformation response spectrum is a plot of uo against Tn for fixed ζ.
A similar plot for 𝑢ሶ 𝑜 is the relative velocity response spectrum and for 𝑢ሷ 𝑜 𝑡 is the total
acceleration response spectrum.
Figure on next slide summarizes the procedure to develop the deformation response
spectrum. The spectrum is developed for 2% damped SDOF systems when subjected to El
Centro ground motions as shown in part (a) of the figure.
19
Deformation Response spectrum
(a) El-centro ground acceleration; (b) Deformation response of three SDOF systems
with ζ=2% and Tn=0.5,1, and 2 sec; (c) Deformation response spectrum for ζ=2%
20
Deformation Response Spectrum
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 and thus require more
computational effort.
23
Pseudo-Acceleration Response Spectrum
The pseudo acceleration spectrum is obtained from the displacement spectrum by multiplying by
2
the circular frequencies squared. 2𝜋
𝜔𝑛 2 𝐷 = 𝐴 = 𝐷
𝑇𝑛
𝐷 ≡ 𝑢𝑜 2
Tn D-in. A=D*(2π/Tn)2 𝐴 = 2𝜋ൗ𝑇 𝐷
𝑛
Note that the acceleration at a near zero period is not near zero (as was the case for velocity and
displacement). In fact, the pseudo acceleration represents the total acceleration in the system while
the pseudo velocity and the displacement are relative quantities.
Pseudo-Acceleration Response Spectrum
The pseudo acceleration response spectrum represents the total acceleration of the
system, not the relative acceleration (?).
For very rigid systems (with near zero periods of vibration), the relative acceleration
will be nearly zero and, hence, the pseudo acceleration, which is the total acceleration, will
be equal to the peak ground acceleration.
At Tn=0, PGA, 𝑢ሷ 𝑔𝑜 = A
A caution about Pseudo responses
1. uo is same as D by definition.
𝑡
3. Similarly, 𝑢ሷ 𝑜 is not taken as A which by definition= ωn2D
26
Difference Between Pseudo-Acceleration , A, and
Total Acceleration, 𝑢ሷ 𝑡 𝑜
ζ=0.05
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% of critical
29
Pseudo Acceleration Response Spectra for Different Damping Values
Damping has a similar effect on pseudo acceleration. Note, however, that the pseudo
acceleration 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.
Pseudo Acceleration response spectra for 1940 El-centro earthquake for different values of ζ 30
Four-Way Log Plot of Response Spectrum (Combined D-V-A spectrum)
≈ 1.5 sec
Combined D-V-A response spectrum for N-S component of El Centro ground motion; ζ = 0,5,10 and 20%
32
Relation between Base shear and Pseudo acceleration
𝑓𝑠𝑜 = 𝑘. 𝑢𝑜
since 𝑘 = 𝜔𝑛 2 . 𝑚 ⇒ 𝑓𝑠𝑜 = 𝑘. 𝑢𝑜 = (𝜔𝑛 2 . 𝑚). 𝑢𝑜
⇒ 𝑓𝑠𝑜 = 𝑚. 𝐴
Base shear, Vb , is an estimate of the
𝐹𝑥 = 0 ⇒ 𝑉𝑏𝑜 −𝑓𝑠𝑜 = 0 maximum expected lateral force on the
base of the structure due to seismic activity
𝑉𝑏𝑜 = 𝑚. 𝐴
𝑊 𝐴
𝑉𝑏𝑜 = .𝐴 = .𝑊
𝑔 𝑔
Where W is the weight of the structure and g is the gravitational acceleration. When written in this form,
A/g may be interpreted as the Base shear coefficient or lateral force coefficient . It is used in building
codes to represent the coefficient by which the structural weight is multiplied to obtain the base shear
33
Problem 5.1
Consider two water tanks , each supported on R.C cantilever hollow circular tower with different
heights. Height of one tower is 60 ft and weight of supported water tank= 10 k. Height of second
tower is 30 ft with weight of water tank is 80 k.
If E and I for both towers = 3000 ksi and I= 5*105 in4, Calculate:
The peak value of displacement response , base shear and bending moment for both systems using
EL-centro 1940 response spectra for ζ=0.05 and comment on the results.
Ignore self weight of tower, P-Δ effect
𝑊1 10 k s2 𝑊2 80 k s2
𝑚1 = = = 0.026 k. 𝑚2 = = 2 = 0.208 k.
𝑔 32.2 ∗ 12 in/s 2 in. 𝑔 32.2 ∗ 12 in/s in.
k k
3 3000 5 ∗ 105 in4 3 3000 2 5 ∗ 105 in4
3𝐸𝐼 in2 in
𝑘1 = 3 = = 12.1 k/in 𝑘2 = = 48.4 k/in
ℎ1 (60 ∗ 12 in)3 (30 ∗ 12 in)3
𝑚1 0.026
𝑇𝑛1 =2𝜋 = 2𝜋 = 0.29 s
𝑘1 12.1
𝑚2 0.208
𝑇𝑛2 =2𝜋 = 2𝜋 = 0.29 s
𝑘2 48.4
𝑢𝑜 = 𝐷 = 0.67″
2𝜋 2 2𝜋 2 0.67
𝐴= 𝐷 = ft = 314.8 in/s2 0.67″
𝑇𝑛 0.29 12
𝑉𝑏𝑜 = 𝑘𝑢𝑜 = 𝑚𝐴
k
𝑉𝑏𝑜 for tank no. 1 = 12.1 ∗ 0.67 in = 8.1 k 𝑀𝑏𝑜 for tank no. 1 = 8.1 ∗ 60 = 486 ft. k
in
k 𝑀𝑏𝑜 for tank no. 1 = 32.2 ∗ 30 = 966 ft. k
𝑉𝑏𝑜 for tank no. 2 = 48.4 ∗ 0.67 in = 32.4 k
in
Conclusion: Since both water tank have same Tn ,therefore, they respond by same amount
when subjected to same EQ (i.e., same uo and A). However, tank having higher stiffness and
mass experience higher value of base shear
Spectral regions in Response Spectrum
𝐴
≈ const. for 𝑇𝑛 < 𝑇𝑐
𝑢ሷ 𝑔𝑜
36
Design Spectrum
Response spectrum cannot be used for the design of new structures, or the seismic safety
evaluation of existing structures due to the following reasons:
Response spectrum for a ground motion recorded during the past is inappropriate for
future design or evaluation.
The response spectrum is not smooth and jagged, specially for lightly damped structures
and, therefore, for a given earthquake, small variations in structural frequency (period) can
produce significantly different results. .
The response spectrum for different ground motions recorded in the past at the same site
are not only jagged but the peaks and valley are not necessarily at the same periods. This
can be seen from the graph given on next slide where the response spectra for ground
motions recorded at the same site during past three earthquakes are plotted
37
1956: 4.8
1940: 3.0
1968: 1.75
Tn = 0.4 sec
Figure: Response spectra for the N-S component of ground motions recorded at the imperial valley Irrigation
38ζ = 2%.
district substation, El Centro, California, during earthquakes of May 18,1940;Feb 9,1956;and April 8,1968;
Design Spectrum
If nothing have been recoded at the site, the design spectrum should be based on
ground motions recorded at other sites under similar conditions such as magnitude of the
earthquake, the distance of the site from causative fault, the fault mechanism, the
geology of the travel path of seismic waves from the source to the site, and the local soil
conditions at the site.
39
Design Spectrum
For practical applications, design spectra are presented as smooth curves or straight lines.
Smoothing is carried , using statistical analysis, out to eliminate the peaks and valleys in
the response spectra that are not desirable for design. For this purpose statistical analysis of
response spectra is carried out for the ensemble of ground motions.
Each ground motion, for statistical analysis is normalized (scaled up or down) so that all
ground motions have the same peak ground acceleration, say 𝑢ሷ 𝑔𝑜 ;other basis for
normalization can be chosen.
40
Design Spectrum (Newark-Hall)
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 Newmark spectrum is based on the following
observations:
The pseudo acceleration 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,
pseudo velocity, and pseudo acceleration are equal
to the ground values times some empirical constant.
41
Structural Response as a function of Ground motion
43
Construction of Design Spectrum (Newmark-Hall)
Newmark and Hall have developed procedures to construct such design spectra from
ground motion parameters. One such procedure is illustrated in given figure.
𝛼𝑉 𝑢ሶ 𝑔𝑜
𝛼𝐷 𝑢𝑔𝑜
The recommended period values Ta = 1/33 sec, Tb = 1/8 sec, Te = 10 sec, and Tf = 33 sec. 44
Construction of Design Spectrum (Newmark-Hall)
Construction of Design Spectrum for firm soil (Newmark-Hall)
𝛼𝑣 = 2.3
Using the values on previous slide
𝛼𝐷 = 2.01
and values given in table 6.9.2 (slide 𝛼𝐴 = 2.71
43) for 84.1 percentile and ζ =5%, the
Pseudo-velocity design spectrum can be
dawn as shown in Fig u re 6 .9 .4
46
Construction of Design Spectrum for firm soil (Newmark-Hall)
The Pseudo-acceleration and displacement design spectra drawn by using above given
equation are drawn in Figures 6.9.5 and 6.9.6 on next two slides
47
Construction of Design Spectrum for firm soil (Newmark-Hall)
Acceleration sensitive region
2𝜋
𝐴 = 𝑉𝜔𝑛 = 𝑉
𝑇𝑛
48
Construction of Design Spectrum (firm soil)
𝑉 𝑇𝑛
𝐷= = 𝑉
𝜔𝑛 2𝜋
49
Design Spectrum for various values of ζ
Pseudo- acceleration design spectrum (84.1 th percentile) drawn on log scale for ground
motions with u go = 1g , u go = 48 in/sec, and u go = 36 in. ; ζ = 1,2,5,10 and 20 %.
51 51
Envelope Design spectrum
For some sites a design spectrum is the envelope of two different elastic design spectra as
shown below
Site
52
Problem 5.2
(a) A full water tank is supported on an 80-ft-high cantilever tower. It is idealized as an SDF system
with weight w = 100 kips, lateral stiffness k = 4 kips/in., and damping ratio ζ = 5%. The tower
supporting the tank is to be designed for ground motion characterized by the design spectrum of Fig.
6.9.5 scaled to 0.5g peak ground acceleration. Determine the design values of lateral deformation
and base shear.
(b) The deformation computed for the system in part (a) seemed excessive to the structural designer,
who decided to stiffen the tower by increasing its size. Determine the design values of deformation
and base shear for the modified system if its lateral stiffness is 8 kips/in.; assume that the damping
ratio is still 5%. Comment on how stiffening the system has affected the design requirements. What
is the disadvantage of stiffening the system?
53
Design Spectra (Building code of Pakistan 2007)
T= natural time period, Ca and Cv= seismic zoning coefficients that depends on
soil type and Zoning factor, Z 54
Soil profile types (Building code of Pakistan 2007)
55
Seismic Zone factor (Building code of Pakistan 2007)
56
Seismic coefficients (Building code of Pakistan 2007)
57
Design Response Spectra (Building code of Pakistan 2007)
58
Base shear determination by static force method
(Building code of Pakistan 2007)
The total design base shear in a given direction shall be
determined from the following formula:
The total design base shear need not exceed the following:
The total design base shear shall not be less than the following:
In addition, for Seismic Zone 4, the total base shear shall also
not be less than the following:
Determine the design values of lateral deformation and base shear for the
data mentioned in problem M 6.2. The tank is required to be constructed in
Zone 2B (as per BCP 07) on soil type SE. Consider elastic design and take
Importance factor, I= 1.0
60
Problems for practice
61