Lecture Notes

LECTURE NOTES
Prepared by Dr. Ashraf A El Damatty, Ph.D., P. Eng. Professor Department of Civil & Environmental Engineering The University of Western Ontario London, Ontario, Canada
Notes can not be copied or reproduced without approval of Dr. Ashraf El Damatty
Acknowledgement
The author would like to acknowledge the contribution of the following people in preparing this set of notes.
Dr. Maged Youssef in developing the section describing the seismic requirements of reinforced concrete structures, Mr. Mohamed Semelawy in preparing the example for seismic design of RC moment resisting frames, Mr. Mahmoud Hassan in preparing the example for seismic design of steel structures, Mr. Ahmed Hamada in preparing the modeling section, and Mr. Mohamed Darwish and Mr. Mohammad Siddique in typing the notes.
List of Proposed References
1) General, good coverage of code Title: Elements of Earthquake Engineering and Structural Dynamics By: Andre Filiatrault Publisher: Polytechnic international press ISBN: 2-553-01021-4
2) Structural Dynamics Author Title Publisher Chopra, Anil K. Dynamics of structures: theory and applications to earthquake engineering / Anil K. Chopra. Englewood Cliffs, N.J. : Prentice Hall, 1995.
3) Steel Structures Author Title Publisher Author Title Publisher Bruneau, Michel, Ph.D. Ductile design of steel structures / Michel Bruneau, Chia-Ming Uang, Andrew Whittaker. New York : McGraw-Hill, 1998. Canadian Institute of Steel Construction Handbook of Steel Construction Universal Offset Limited, 9th Edition, 2007.
4) Concrete and Masonry Structures Authors Title Publisher Author Title Publisher Thomas Paulay, and M. J. N. Priestley Seismic Design of Reinforced Concrete and Masonry Buildings/ Thomas Paulay, and M. J. N. Priestley Willey- Interscience, 1992 Cement Association of Canada Concrete Design Handbook Ottawa, 3rd Edition, 2006
CHAPTER 1 EARTHQUAKE GROUND MOTIONS CHARACTERISTICS
Seismic Analysis and Design of Buildings
By: Prof. A. A. El Damatty
Page: 1-1 .
1- Causes and Effects of Earthquakes The earth is divided to six major tectonic plates; these are the Eurasian, Pacific, American, African, Indian, and Antarctic. The vast majorities of damaging earthquakes originate at, or adjacent to, the boundaries of tectonic plates, due to relative deformations at the boundaries. Because of the nature of the rough interface between adjacent plates, a stick-slip phenomenon generally occurs, rather than smooth continuous relative deformation, or creep. The relative deformation of the adjacent plates is resisted at the rough interface by friction. When the induced stresses exceed the frictional capacity of the interface, slip occurs releasing the elastic energy stored in the rock in the form of shock waves propagating through the medium at the ground-wave velocity. The inertial response of structures to the ground accelerations resulting from the energy released during fault slip is the primary interest to structural engineers. 2- Seismic Waves The rupture point within the earths crust represents the source of emission of energy. It is known as hypocenter, focus, or source The epicenter is the point on the earths surface immediately above the hypocenter. The Focal depth is the depth of the hypocenter below the epicenter. The Focal distance is the distance from the hypocenter to a given reference point.
Page: 1-2 .
The energy released by earthquakes is propagated by different types of waves. The body waves, originate at the rupture zone and include P waves (primary or dilatation waves) and S waves (secondary or shear waves). The P waves involve particles moving parallel to the direction of propagation of the wave. The S waves involve particles moving perpendicular to the direction of propagation of the wave. P and S waves have different characteristic velocities Vp and Vs. The velocities usually satisfy the following relation:
Vp Vs
As a consequence, the time interval T between the arrival of P and S waves at a given site is thus proportional to the focal distance Xf, hence: Xf = Vp T/(3-1) Recording the P-S time interval at three or more sites enables the epicentral position to be estimated as shown in Fig. 1
X1 T1V p 3 1 T2V p 3 1 X2
T3V p 3 1 X3
X1, X2, and X3 are recording stations
Fig. 1 Determination of Epicenter
Page: 1-3 .
3-Characteristics of Earthquakes The effect of earthquakes on buildings depends on the seismic ground motion that the structure experiences. In general, there are six components of seismic ground motions; three translation and three rotational components. Current strong motion measuring instruments only record three translation components. There is no reliable measurement on the rotational components at present. Seismic ground motion is very complex. There is a continuous effort to understand what are the relevant parameters that denote the damage potential of such motions to buildings. Factors that affect the damage potential are: Peak ground acceleration (PGA) Peak ground velocity (PGV) Peak ground displacement (PGD) Duration of Strong Shaking Frequency Content The current opinion is that one needs at least two parameters to describe the damage potential; one parameter denotes the intensity of shaking and frequency content; and the other denotes the duration of shaking. The combination of level of shaking and frequency content is reflected in a single quantity called Response Spectrum which will be described later.
Page: 1-4 .
a) Magnitude The accepted measure of magnitude is the Richter scale which is related to the maximum deformation of the surface-wave motion recorded by a standard Wood-Anderson seismograph located at a distance of 100 km from the epicenter. The accepted relationship between energy released E, and Richter magnitude M is: log10E = 11.4 + 1.5M It is clear that an increase in one Richter unit leads to about 30 times increase in the energy released.
(b) Intensity Earthquake intensity is a subjective estimate of the effects of an earthquake and is dependent on peak acceleration, velocity, duration, and frequency content. The most widely used scale is the modified Mercalli, MM, scale which was refined by Richter in 1958. The effective range is from MM2, which is felt by persons at rest on upper floor of building to MM12 where damage is nearly total. A listing of the complete scale is given below:
Page: 1-5 .
Masonry A : Good workmanship, mortar, and design; reinforced, especially laterally, and bound together by using steel, concrete, etc.; designed to resist lateral forces. Masonry B : Good workmanship and mortar; reinforced, but not designed in detail to resist lateral forces. Masonry C: Ordinary workmanship and mortar; no extreme weaknesses like failing to tie in at corners, but neither reinforced nor designed against horizontal forces. Masonry D: Weak materials, such as adobe; poor mortar; low standards of workmanship; weak horizontally. Modified Mercalli Intensity Scale of 1931 (abridged and rewritten by C. F. Richter)
Page: 1-6 .
Page: 1-7 .
4- Characteristics of Earthquake Accelerograms (a) Accelerograms Accelerographs record ground accelerations in optical or digital forms as time history records. Integration of the records enables velocities and displacements to be estimated. Three different time history accelerations are shown in Fig. 2
Fig. 2 Time History Earthquake Records
Page: 1-8 .
(b) Vertical Acceleration Vertical accelerations recorded by accelerographs are generally lower than corresponding horizontal components and frequently are richer in low-period components. It is often assumed that peak vertical accelerations are approximately two-thirds of peak horizontal values. However, there is an increasing evidence that peak horizontal and vertical components are similar for near-field records. (C) Influence of Soil Stiffness It is generally accepted that soft soils modifies the characteristics of strong ground motion transmitted to the surface from the underlying bed rock. Soft soils usually result in amplification of long-period components and reduction of the short-period ones.
Fig. 3 Comparison of lake bed (1-3) and rock (4-6) Accelerographs, Mexico City 1985, Paulay and Priestley (1992)
Page: 1-9 .
Fig. 3 shows accelerations recorded at adjacent sites on rock and on medium depth lake deposits in the 1985 Mexico earthquake. The epicenter of this earthquake was approximately 400 km from Mexico City and the peak bedrock acceleration was about 0.05g. It can be noticed from Fig. 3 that this acceleration was amplified about five times by the elastic characteristics of the lake bed deposits. A modified ground motion, having predominant energy in the period range of 2 to 5 sec was generated. As a consequence, buildings with natural periods in this range were subjected to extremely strong motions with many failures. (d) Geographical Amplification Geographic features may have a significant influence on local intensity of ground motion. In particular, steep ridges may amplify the base rock accelerations by resonance effects. A structure built on top of ridge may thus be subjected to intensified shaking. Extensive damage Buildings abandoned
No damage 50m 200 m Fig. 4 1000 m
Page: 1-10 .
5- Attenuation Relationships The attenuation relationship represents the reduction in peak ground acceleration with distance from epicenter. Three major factors contribute to the attenuation: 1) The energy released from an earthquake may consider to be radiated away from the source as a combination of spherical and cylindrical waves. The increase in surface area of the waves as they move away from the source implies that accelerations will decrease with distance as the sum of a number of terms proportional to Re-1/2, Re-1, Re-2 and ln Re ; where Re is the distance to the point source or cylindrical axis. 2) The energy is reduced due to material attenuation or damping of the transmitting medium. 3) Attenuation may result from wave scattering at interfaces between different layers of material. For small-to-moderate earthquakes, the source could reasonably be considered as a point source, and spherically radiating waves would characterize attenuation. For large earthquakes, with fault movement over several kilometers, cylindrical waves might seem more appropriate. Hence attenuation relationships for small and large earthquakes might be expected to be different. Typical attenuation relationships taken the form:
a0 =C1e C 2 M ( Re +C3 ) C 4
Page: 1-11 .
where: a0 is the peak ground acceleration M is the Richter magnitude of the earthquake Re is the epicentral distance C1, C2, C3, and C4 are constants Many relationships of the general form is given by above equations have been proposed by various researchers. 6- Relation between Intensity and Ground Acceleration Fig. 5 shows the relationship between MM intensity and peak acceleration resulting from a number of studies.
Fig. 5 Relationship between Intensity and Peak Acceleration (Paulay and Priestley (1992))
To some extent the scatter exhibited by Fig. 5 can be reduced by use of effective ground acceleration (EPA) which is related to the peak response
Page: 1-12 .
acceleration of short period elastic oscillation rather than the actual maximum ground acceleration. Despite the wide scatter in Fig. 5, there is an approximate linear relationship between the logarithm of PGA and MM intensity, I. An average relationship may be written as: (PGA)Ave = 10-2.4 + 0.34I m/sec Meanwhile, a conservative estimate for (PGA) can be evaluated using the relationship: (PGA) = 10-1.95 + 0.32I m/sec
Page: 1-13 .
7- Return Periods: Probability of Occurrence To assess the seismic risk for a given site, it is necessary to know not only the characteristics of strong ground motion that are feasible for the site, but also frequency with which such events are expected. It is common to express this by the return period of an earthquake of given magnitude, which is the average recurrence interval for earthquakes of equal or longer magnitude. Generally, larger earthquakes occur less frequently than small ones. The probability of occurrence (inverse of the return period) of earthquakes of different magnitude M can be expressed by the following distribution: (M) = Ve-M (M) is the probability of an earthquake of magnitude M or greater occurring in a given volume V per unit volume. and are constants related to the location of a given volume V. Fig. 6 shows data for different tectonic zones compared with predictions of the above equation.
Page: 1-14 .
Fig. 6 Magnitude-Probability Relationship (Paulay and Priestley (1992)) 8- Design intensity The intensity of ground motion adopted for design depends on the seismicity of the area and also on the level of structural response contemplated and the acceptable risk associated with that response. Three levels of response are usually identified: 1) Serviceability limit state where building operations are not disrupted by the level of ground shaking. 2) A damage control limit state where repairable damage to the building may occur. 3) A survival limit state under an extreme event earthquake, where severe and possible irrepairable damage might occur but collapse and loss of life are avoided.
CHAPTER 2 RESPONSE OF A SINGLE DEGREE OF FREEDOM SYSTEM
Page: 2-1 .
2-1 Equation of motion of a SDOF

Basic components of a dynamic problem: 1. 2. 3. 4. 5. Mass (m = W/g) Applied load P(t) Inertia force = FI Damping force = FD Elastic force = FS
Single Degree of Freedom Model: C M K Using DAlemberts principle; a dynamic system can be considered to be in equilibrium if the inertia force is included in the free body diagram. FD P(t) FI FS FI + FD + FS = P(t) Where:
d2x x FI = mass x acceleration = m 2 = m&& dt
x(t) P(t)
& FD = Dashpot constant x velocity = cx FS = Spring stiffness x displacement = kx
Equation of motion:
m&& + cx + kx = P(t) x &
Page: 2-2 .
Elastic Stiffness for Elements: k = P/ Force Required for Unit displacement
1) Springs in Series: K1 P1 P1 = P2 = P P1 = K11 ;

K=
K2 P2
P,
P2 = K22 ;
= 1 + 2
KK P P P = = = 1 2 1 + 2 P P2 K1 + K2 1 + K1 K2 KK K = 1 2 K1 + K2
2) Springs in Parallel: P1 , 1 , K 1 P, , K P2 , 2 , K 2 1 = 2 =
K= P P + P K11 + K22 = 1 2=
K = K1 + K2
Page: 2-3 .
3) Axial Stiffness of a Bar: P L 4) Torsion: T L 5) Beams in Flexure: a) Cantilever: P

PL3 3 EI = K = L3 3 EI
PL EA
EA L
K =
TL GJ
GJ L
K =
b) Simple Beam: P
PL3 48 EI = K = L3 48 EI
Page: 2-4 .
c)
PL3 12 EI = K = L3 12 EI
d)
PL3 3 EI = K = 3 EI L3
e)
K1
K2
K1 =
12 EI 1 3 L1
K2 =
12 EI L3 2
K1 and K2 are two springs in parallel (have the same displacement)

K = K1 + K 2 = 12E ( I1 + I 2 ) 3 L1
Page: 2-5 .
Types of Dynamic Loading: a) Free Vibration: X(t)

& X0
X0
0 2 4 6
Time
b) Forced Harmonic Vibration: P(t)

P (t ) = P0 sin t
P0
2 4 6 Time
(Machine Vibration)
c) Impulse Loading: (Driving of Piles)
P(t)
P0
Time d) Random Vibration: P(t)

0 200 400 600 800 1000 1200 1400 1600 1800
Time
Page: 2-6 .
2-2 Free Vibration Response of a SDOF

x Equation of motion: m&& + kx = 0
&& + 0 x = 0 x 2
2 0 =
k m
Solution of the above D.E.: X(t) = Asin(0t) + Bcos(0t)
& & Substituting with the initial conditions: X(0) = X0 and X (0) = X 0 & X0 we obtain: X(t) = sin(0t ) + X 0cos(0t )
0
The above equation can be written in a compact form as:
& 2 X0 2 X(t) = + X 0 [sin( 0 t + )] 0
= tan1 & 0 X / 0 0
& X0
X(t) X0 T = 2/0
Page: 2-7 .
2-3 Undamped Free Vibration of a Damped System

Equation of motion:
m&& + cx + kx = 0 x & k c && + x + x = 0 & x m m c && + x + 0 x = 0 & 2 x m

(damping ratio)
Define: = 2m 0
2 & x Equation of motion: && + 2 0 x + 0 x = 0
For 0 < < 1 The solution of the above equation becomes: X(t) = e
0t
[A1sin(Dt ) + A2cos(Dt )]
2
Where: D = 1 0 (Frequency of vibration of a damped system)
A1 =
& X 0 + 0 X 0
A2 = X 0
The above equation can be written in a compact form as: X(t) =
Xe0t [sin(Dt + )]
Page: 2-8 .
1 0.8 0.6 0.4 0.2 X(t) 0 -0.2 0 -0.4 -0.6 -0.8 -1 Time 5 10 15 20
TD = 2/D
Response to Harmonic Motion C M K

2 & x Equation of motion: && + 2 0 x + 0 x =
x(t) P(t) = P0cost or P0sint
P0 cos t _____________(*) m The solution consists of a homogeneous solution and a particular one:
Homogeneous: XC(t) =
e 0t [A1sin(Dt ) + A2 cos(Dt )]
Try a particular solution in the form: XP(t) = Csin(t) + Dcos(t)
Page: 2-9 .
Substituting for the particular solution into (*) and solving for the constants C and D we obtain: XP(t) = XP cos(t ) Where:
XP =
P0 2 m0 2 2 1 + 2 0 0
2
Amplitude of forced response
2 1 = tan 0 2 1 0
The angle at which the response lags the applied load
Adding the complementary and particular solutions, the total response is: P0 cos(t ) 2 m0 X (t ) = e0t [A1sin(Dt ) + A2cos(Dt )] + 2 2 2 1 + 2 0 0
transient term Decays exponentially and occurs at the damped frequency of the structure
steady state term Takes place at the forced frequency. Lags behind the force by phase angle .
Page: 2-10 .
Practically in design the steady state response is only considered. Define DMF Dynamic Magnification Factor:
XP =
P0
2 m0
1 2 2 1 + 2 0 0
2
P P X St = 0 = 0 2 i.e. K m0
X P = X St .DMF
1
DMF =
2 2 1 + 2 0 0
Dynamic Magnification Factor
12 11 Dynamic Magnification Factor 10 9 8 7 6 5 4 3 2 1 0 0 0.5 1 1.5 /0 2 2.5 3 =0 = 0.05 = 0.1 = 0.2
Page: 2-11 .
Example 2-1 A simply supported beam of length L carries an electric motor of weight W = 12 KN at its mid span. The flexural rigidity EI of the beam is such that the static deflection under the weight of the machine is 3 mm. The equivalent viscous damping is such that after 10 cycles, the free vibration amplitude due to initial displacement X0 is reduced to 60 % of X0. The motor operates at 600 RPM and causes a harmonic force due to unbalance F0= 2.5 KN. Neglecting the mass of the beam, determine: a) The value of the undamped frequency 0. b) The value of the damped frequency D. c) The magnitude of the maximum dynamic displacement X for the given values of the parameters due to the operation of the machine. d) The maximum total deflection including gravity and dynamic loading.
Page: 2-12 .
Page: 2-13 .
Page: 2-14 .
Page: 2-15 .
2-4 Equation of Motion of a SDOF under Earthquake Motion

C M K g(t) X(t)
The figure above shows a simulation for a SDOF system subjected to an earthquake ground motion described by a time history ground acceleration g(t). X(t) represents the displacement response relative to ground of this SDOF system. The equation of motion of the SDOF is given as:
& m( &&g (t) + &&(t)) + cx + kx = 0 y x & m&&(t) + cx + kx = m&&g (t) x y

As such, due to earthquake, the SDOF is subjected to a dynamic y force P(t) = m&&g (t) Dividing the previous equation by the mass m, we get:
Page: 2-16 .
&&(t) + x
c k & x + x = &&g (t) y m m
k Defining the natural frequency of the system 0 = m and its c = damping ratio 2m0 , the equation of motion can be written as: y && + 2 0 x + 02 x = &&g (t) & x
The solution of the above 2nd order D.E. can be obtained by Time History Analysis Response Spectrum Analysis
Provides the entire time response.
Provides only maximum responses.
Can be applied to linear and non-linear systems. Involves significant computations.
Applicable only to linear systems.
Simple
Page: 2-17 .
2-5 Time History Analysis

Closed form mathematical solutions exist for the free vibration response of a SDOF system as well as the responses to harmonic and impulse loading. Dynamic problems which involve forces that can not be expressed in simple mathematical forms (example earthquake forces) can not be solved in a closed form. For such problems, the response is obtained by applying numerical integration techniques to the equation of motion. In general, the numerical integration techniques can be categorised to: 1) Explicit Algorithms. 2) Implicit Algorithms. One of the most widely applied implicit algorithms will be discussed in details in this section. Implicit Algorithms Newmark method U Wilson method Newmark method We will assume:
& & x n+1 = x n + (1 ) &&n + &&n+1 t x x
& x n+1 = x n + x n t + 1 &&n + &&n+1 t 2 _______________ (2) x x 2
[(
_______________ (1)
Page: 2-18 .
and are implicit parameters that can be determined to obtain integration accuracy and stability. = 1/2 and = 1/6 correspond to linear acceleration. = 1/2 and = 1/4 correspond to constant average acceleration. Equation of motion at t = (n+1)t:
& m&&n+1 + cxn+1 + kxn+1 = F n+1(t) x
_______________ (3)
For constant average acceleration scheme ( = 1/2 and = 1/4):
& & x n+1 = x n + &&n + &&n+1 (t / 2) x x
_______________ (4)
& x n +1 = x n + x n t + &&n + &&n +1 t 2 / 4 x x
](
_______________ (5)
From (5):
&&n +1 = x n +1 x n x n t 4 / t 2 &&n & x x
](
_______________ (6)
Substituting (6) into (4), we get:

& & & x n +1 = x n + &&n (t / 2 ) + x n +1 x n x n t 4 / t 2 &&n (t / 2) ___ (7) x x
[[
](
Equation (7) can be simplified to give: & & x n +1 = (x n +1 x n )(2 / t ) x n Substituting (6) and (8) into (3), we get:
_______________ (8)
4xn 4xn n 2xn n & 4m 2c n+1 n+1 & x + x ___ (9) + && + c 2 + + k x = F (t) + m 2 + t t t t t
Page: 2-19 .
StepByStep Solution Using Newmark method ( = , = ) A- Initial Calculations: 1. Evaluate k, m and c.

0 &0 x0 2. Initialize x , x ,and &&
3. Select a time step t. 4. Calculate: a0 = 4/t2, a1 = 2/t, a2 = 4/t, a3 = t/2 .

5. Evaluate k = k + a0 m + a1c
B- For Each time step:

1. Calculate effective load at time t = (n+1)t; F n +1 & x & F n+1 = F n+1 + m a0 x n + a2 x n + &&n + c a1 x n + x n
) (
2. Solve kx n +1 = F n +1
3. Calculate acceleration and velocity at t = (n+1)t

&&n+1 = a0 x n+1 x n a2 x n &&n & x x
& & x n +1 = x n + a3 &&n +1 + &&n x x
Page: 2-20 .
Example 2-2 The frame structure shown below is subjected to a strong ground motion pulse with a peak acceleration A0 = 1.2g. The latter is modeled as a triangular pulse load of duration t1 equal to half the period of vibration of the structure. Assuming zero damping, determine by numerical integration the maximum displacement of the structure, Xmax, and also the time, tmax, at which it occurs. Use a time step t=T/10 sec. W = 200 KN I= 3.0 m W 250x33 4.5 m g(t) t1= 0.5T A0= 1.2g t1/2 t1 t X
g(t) A0
Page: 2-21 .
Page: 2-22 .
Page: 2-23 .
Page: 2-24 .
Page: 2-25 .
2-6 Concept of Elastic Response Spectrum

As a central concept of earthquake engineering, the response spectrum provides a convenient means to summarize the peak response of all possible linear SDOF systems to a particular component of ground motion. Construction of an Elastic Response Spectrum The following steps are conducted in order to construct the response spectrum of a given ground motion. 1. Consider a SDOF subjected to a ground motion given by the time history analysis acceleration g(t). The SDOF system
2 m = 2 T= has a certain period 0 k and a damping ratio . M x

C K/2 2. Equation of Motion: K/2 g(t) x is the relative displacement between the mass and the ground
& m( &&g (t) + &&(t)) + cx + kx = 0 y x & m&&(t) + cx + kx = m&&g (t) x y

2 y & x Or && + 2 0 x + 0 x = &&g (t)
Page: 2-26 .
Knowing the time history record g(t), the response of the structure x(t) can be evaluated using any of the numerical methods.
The maximum value Sd for the response displacement is recorded. This maximum displacement Sd(T,) depends on the period and the damping ratio of the SDOF. Sd(T,) is called the spectral displacement. 3. The maximum velocity response Sv(T,) can be obtained by the relation: Sv(T,) = 0Sd(T,) 4. The maximum acceleration response Sa(T,) can be obtained by the relation: Sa(T,) = 0Sv(T,) = 02Sd(T,) 5. Steps (1) to (4) are repeated for different values for the natural period of the system T as well as for different damping ratios .
Page: 2-27 .
Page: 2-28 .
The following observations could be noticed: 1. Systems with very short period, say T < Ta = 0.025 s have Sa= g0 and Sd is very small. Why?? 2. Systems with very long period, say T > 15 s have Sd approaches g0 and Sa is very small. Why?? 3. For shortperiod systems, Sa can be idealised as constant. 4. For mediumperiod systems, Sv can be idealised as constant. 5. For longperiod systems, Sd can be idealised as constant. 6. Sa, Sv and Sd are plotted together on a log-log scale as function of T and .
Figure 1 Response Spectrum (=0, 2, 5 and 10%) for El Centro earthquake
N.B. Each earthquake has its own response spectrum
Page: 2-29 .
Based on the above observations one can conclude the following: The long period region is called the displacement sensitive region as the structural response is most directly related to ground displacement. The medium period region is called the velocity sensitive region as the structural response is most directly related to ground velocity. The short period region is called the acceleration sensitive region as the structural response is most directly related to ground acceleration. Elastic Design Spectrum Based on the analysis of a large ensemble of ground motions recorded on firm ground an elastic design spectrum was developed by researchers as shown.
Figure 2 Elastic Design Spectrum
Page: 2-30 .
The amplification factors for two different non exceedance probabilities, 50% (median) and 84.1% (median plus one standard deviation) are given in Table 1 for several values of damping. Table 1 Amplification Factors: Elastic Design Spectra Median One sigma Damping, (50 Percentile) (84.1 Percentile) (%) A V D A V D 1 3.21 2.31 1.82 4.38 3.38 2.73 2 2.74 2.03 1.63 3.66 2.92 2.42 5 2.12 1.65 1.59 2.71 2.30 2.01 10 1.64 1.37 1.2 1.99 1.84 1.69 20 1.17 1.08 1.01 1.26 1.37 1.38
Page: 2-31 .
2-7 Seismic Response of Structures Using Response Spectrum for a SDOF System
The elevated tank shown in the opposite figure can be treated as a SDOF system having a period T given by: H K m = 2 T= 0 k The maximum response of the tank to an earthquake (e.g. El M
Centro earthquake) can be obtained by applying the following steps: 1. A damping ratio is assumed for the structure. 2. Using T and , the response spectrum of El Centro can be used to evaluate the maximum relative displacement Sd experienced by the mass M during the earthquake. 3. The maximum velocity Sv(T,) = 0Sd(T,) 4. The maximum acceleration Sa(T,) = 0Sv(T,) 5. The maximum base shear Vmax = MSa = M02Sd 6. The maximum overturning moment Mmax= MHSa= MH02Sd
Page: 2-32 .
Example 2-3 Consider a steel spherical elevated water tank supported by a single circular pipe. The diameter of the tank is 5 m and it is made of steel 0.01 m thick. The circular steel pipe is 1 m in diameter with a wall thickness of 0.02 m wall thickness. The height of the water tank (center of tank to ground level) is 22.5 m. Assuming the structure can be treated as a SDOF system with the total mass concentrated at the water tank level and the foundation condition is taken as rigid, a) Calculate the maximum displacement, base shear force and the overturning moment caused by the 1940 El Centro earthquake. (Assume 2% damping ratio) b) How does the maximum displacement and base shear change if the thickness of the supporting pipe is doubled?
Page: 2-33 .
Page: 2-34 .
Page: 2-35 .
CHAPTER 3 SEISMIC ANALYSIS OF MULTI DEGREES OF FREEDOM STRUCTURES
Page: 3-1
3-1 Multi Degree of Freedom Systems

a) Single Storey Building
12 EI K= 3 L
3 EI K = 3 L
3 EI K = 3 L
3 EI K = 3 L
K=
12 EI L3
Since all columns have the same displacement, they can be treated as springs in parallel:
K = Ki = 12 EI 3EI 12 EI + 3 3 + 3 3 L L L
By neglecting the mass of the columns compared to the mass of the floor, the frame can be treated as a SDOF C M K = Ki
Page: 3-2
b) 2-Storey Shear Building M2 U2 K2 M1 U1 K1
P2(t)
P1(t)
K1 is the stiffness of the 1st storey. K2 is the stiffness of the 2nd storey. 2 DOF System U1 C1 K1 K1U1 C11 M1 M11 M1 C2 K2 K2(U2 U1) C2(2 1) P1(t) K2(U2 U1) C2(2 1) M2 M22 P2(t) M2 U2
Page: 3-3
Equations of motion: M11 + K1U1 + C11 - K2(U2 U1) - C2(2 1) = P1(t) M22 + K2(U2 U1) + C2(2 1) = P2(t) The two above equations can be written in matrix form as:
&& & M1 0 U1 K1 + K2 K2 U1 C1 + C2 C2 U1 P (t) 1 + + & = && 0 M K2 U2 C2 C2 U2 P2 (t) U 2 2 K2

Or:
&& & [M ]{U}+ [K]{U}+ [C]{U}= {F(t)}

For earthquake excitations: P1(t)= -M1g(t) P2(t)= -M2g(t)
Page: 3-4
3-2 Undamped Free Vibration of MDOF Systems (Determination of Natural Frequencies and Mode Shapes)
&& [M ]{U}+ [K]{U} = {0} _____________________________(1) Assuming that: {U} = {}sin( t + )
Or:
{U} = ei(t+ ) {}
{&&} = 2 {U} U
Now equation (1) becomes:

2
[[K] [M]]{} = {0} (Eigen Value Problem) ________________(2)
The non-trivial solution of the above equation is obtained when:

2 det [K] [M ] = 0 (characteristic equation)
For n degrees of freedom, the above will lead to n values of 1, 2, n. Each value of r can be back-substituted into equation (2) to obtain the equivalent mode shape r.
Page: 3-5
Example 3-3 Find the natural frequencies and storey reinforced concrete structure. Where: M1 = M2 = 30,000 Kg ES= 200,000 MPa Ic= 0.005 m4
M2
M1
10 m
4m
4m
mode shapes for the following two
Page: 3-6
Page: 3-7
Page: 3-8
Page: 3-9
3-3 Mathematical Properties of the Mode Shapes

1. If [K] and [M ] are positive definite matrices, therefore for n degrees of freedom, 12, 22, n2 are real roots. 2. If
{i }
T
and
{ }
j
are
two
different
mode
shapes
{i } [M ]{ j } = {i } [K]{ j } = 0
T
3.
{ } form a complete set of vectors; i.e. any function can be

j
expressed by these shape functions. 4. Normalization of {}: Evaluate 5.
1 {j } j T mj = { j } [M ]{ j } j = j m
{}
{ } [M]{ }=1
T j
{i }T [M ]{i } = mi*
1 = If
{i } [K]{i } = i2mi*
T
6.
[]
{ } { }
2
........ n Ortho-normalized modes nxn
{}
12 2 T 2 [K] = T [M ] = [I ] and
[] []
[] []
2 n nxn
By: Prof. A. A. El Damatty Page: 3-10
3-4 Modal Analysis Procedure to Evaluate Seismic Response of MDOF Systems

&& & [M ]nxn{U}nx1 + [K]nxn{U}nx1 + [C]nxn{U}nx1 = {F(t)}nx1
Assume that:
{U}nx1 = X1(t){1}+ X2 (t){2}+......X n (t){n }

i.e. {U}nx1 = Xi (t){i } = []nxn{X}nx1
i=1 n
By pre-multiplying the equation of motion by {} , the following

T
equation is obtained:
& & [ ]T [M ][ ]{X&}+ [ ]T [K ][ ]{X } + [ ]T [C][ ]{X } = [ ]T {F (t)}

The product [ ] [M ][ ]will lead to the following diagonal matrix:
T
m1 0 0 0 m2 0 0 T [ ] [M ][ ] = m3 0 mn
Where mi = {i } [M ]{i }
T
Meanwhile, the product [ ] [K ][ ] will also lead to the following

T
diagonal matrix:
m112 0 0 0 2 0 0 m22 T [ ] [K ][ ] = 2 0 m33 2 mnn

Assuming that the product [ ] [C][ ] also lead to a diagonal matrix,
T
a total number of n decoupled equations of motion are obtained. Each equation of motion corresponds to a certain mode shape i and is given as:
&&i + 2 i i x i + i2 x i = & x Fi ( t ) T where Fi (t ) = [i ]1xn {F (t )}nx1 mi
For earthquake excitations:
& {F (t)}T = { M1Y&g (t)
& & M 2Y&g (t ) ........ M nY&g (t )}
As such, the equation of motion of the ith mode can be written as:
&&i + 2 i i x i + i2 x i = & x Li && Yg (t ) mi
Where:
mi = Li =
M
j =1
2 ji
M
j =1
ji
ji is the component j of the mode i

i = Li / m i is called the modal participation factor for
mode i
Using the response spectrum technique:
(x i )max
& (x i )max
= i S d = i S v = i S a
= i S d { ji }
(&&i )max x
{U }
ji max
& {U } & {U& } {F }
ji max
= i S v { ji } = i S a { ji }
= i S a { ji M j }
ji max
ji max
{Vi }max
L2 = i S a ji M j = i S a mi j =1
n
Since the maximum responses of the ith modes do not occur at the same time, a certain rule for modal combination has to be adopted. Modal Combination: Let the modal maximum responses be: Q1, Q2,.Qj
2 2 2 1. Square root of sum of squares (SRSS): Q = Q1 + Q2 + ...Q j + ..
2. Absolute Sum: Q = Q1 + Q2 + .... Q j + ... 3. Complete Quadratic Combination (CQC):

Q = Qi2 +
2 i =1 i j n
Qi Q j
2 1 + ij
1 i j ij = where i + j
Q1Q2 2 1 + 12
2 2 2 e.g. for 2 modes: Q = Q1 + Q2 + 2
Sequence of Modal Combination: Mode 1 Mode 2 Mode 3 Combination
Forces
Shear
Moment
Example 3-4 Use the modal analysis procedure to evaluate the response of the frame shown in example 3-3 to El Centro earthquake.
3-5 Earthquake Response of MDOF Using Time History Analysis

The equation of motion for each vibration mode is given as:
&&i + 2 i i x i + i2 x i = & x
Li && Yg (t ) mi
The above equation can be solved using the Newmarks method to obtain xi(t). The nodal displacements associated with mode i can be then obtained as follows:
{U } = X (t){ }
ji i ji
The time history variation of the nodal displacements due to all mode shapes is then given as:
{U }
j nx1
= X i (t){ ji }
n i=1
&& Also: U j
{ }
nx1
&& = X i (t){ ji }
n i =1
Knowing the time history variation of the displacements all other quantities such as the internal forces can be evaluated.
Typical time history responses at the top level of a three storey building due to each one of the three modes of vibration together with the total time history response at the same level are shown below.
The time history variations of the base shear force associated with each one of the three modes of vibration as well as the time history variation of the total base shear force are given below.
CHAPTER 4 CODE PROCEDURES FOR EARTHQUAKE RESISTANT STRUCTURES
Page: 4-1 By: Prof. A. A. El Damatty_____
Inelastic Behaviour and ductility Most of the modern seismic design codes allow inelastic behaviour to occur in structures when they are subjected to strong earthquake motions. A load reduction factor R is usually applied to the earthquake load intensity. Structures are designed under earthquake loads which are less than what are expected if the structures remained elastic. Inelastic behaviour is expected to occur when the structure is subjected to a strong ground motion. If the structure possesses enough ductility, the energy of the earthquake will dissipate through inelastic hysterisis loops that yielded members undergo.
Accordingly, ductility is a very important property that governs the earthquake response of the structure.
Local ductility l = m/y
Global ductility
The global ductility g of a building is usually defined by the following relation:

g y
g =
where : g is the maximum displacement that can be sustained by the structure before collapse. y is the displacement at which the first plastic hinge forms in the structure.
The relation between the load factor R and the global ductility g can be established using one of the following concepts: 1) Equal Displacement F R = d/c g = a/b d c a b From similarity: g = R
2) Equal Energy F R = d/c g = a/b (b-a) c =1/2 d (ad/c) (d/c)2 = 2(b-a)/a d c a b R2 = 2(g -1) R = 2(g-1)
Load Reduction Factor (Rd)
Global Ductility g
Local Ductility, l
V=
S (Ta ) M v I EW R d R0
S (2.0) M v I EW Rd R0
Ta S (Ta) Mv IE Rd Ro W
: Fundamental period of the structure : Design spectral response acceleration : Factor to account for higher mode effect on base shear : Importance factor : Load reduction factor : Calibration factor : Dead Load + 25% Snow Load + 60% Storage Load
According to the NBCC 2005 Code, the seismic hazard at various locations in Canada is defined by five ground motion properties: Sa(0.2), Sa(0.5), Sa(1.0), Sa(2.0), and PGA Sa(T) represents the 5% damped spectral acceleration for period T in seconds, expressed as a ratio of the gravitational acceleration. Seismic Data for Toronto and Montreal is given below: Location Toronto Montreal Seismic Data Sa (0.2) 0.26 0.69 Sa (0.5) 0.13 0.34 Sa (1.0) 0.055 0.14 Sa (2.0) 0.015 0.048 PGA 0.17 0.43
Depending on the soil conditions at the site, the design spectral acceleration S(T) is determined using the set of equations provided on the following page.
Design Spectral Acceleration
Definition of Soil Type
Mv depends on the period of the structure Ta. Mv increases with an increase in Ta (for flexible buildings). It reaches a value of 2.5 for Ta 2.0 seconds
The values of Rd, Ro are given in Table 4.1.8.9 for different structure systems.
Notice, that there is a limit on the height of the building for some of the lateral resisting systems depending on the seismic zone, which is defined by the parameters IEFaSa(0.2) and IEFvSa(1.0)
Estimate of Ta
Moment Resisting Frames => Ta = 0.085(hn)3/4 steel = 0.075(hn)3/4 concrete = 0.1N others
------ (1)
Braced Frame => Ta = 0.025hn Shear Walls => Ta = 0.05(hn)3/4
----------- (2) ---------- (3)
The period can be calculated using computer model leading to a value Ta. However, Ta should not be taken greater than (1) For Moment resisting Frame, 1.5Ta as calculated in (1) (2) For Braced Frame, (3) For Shear Walls, 2.0 Ta as calculated in (2) 2.0 Ta as calculated in (3)
Force distribution Ft Fx =
(V Ft )W x h x
n (Wi hi ) i =1
Ft = 0.07Ta V 0.25V For Ta 0.7 sec. => Ft = 0
Fx hx
hn
Reduction in Moment Mx = Jx Fi (hi hx) Jx = 1.0 for hx 0.6hn

(1 J )h x 0.6hn
Jx = J +
for hx < 0.6hn
J depends on Ta and the type of lateral resisting system.
(Table 4.1.8.11)
Definition of Irregularities
A structure is considered to be irregular if it has any of the type of irregularities described in the above table. Equivalent Static Forces Analysis may be used if any of the following criteria is met: a) IEFaSa(0.2) <0.35 b) Regular structure, with hn 60.0 m and Ta 2.0 sec. (in both two orthogonal directions)
c) Structure with structural irregularities of Type 1, 2, 3, 4, 5, 6, or 8 with hn 20 m and Ta 0.5 sec. (in both two orthogonal directions)
*Otherwise dynamic analysis has to be carried.
Accidental Torsion
Bx = max/ave max = maximum storey displacement at the extreme points induced by equivalent static force applied at distances 0.1Dnx ave = average displacement at the extreme points induced by equivalent static force applied at distances 0.1Dnx When only equivalent static analysis is required (i.e. no dynamic analysis is required) Tx = Fx (ex + 0.1 Dnx) Tx = Fx (ex - 0.1 Dnx)
Direction of Loading If components of the lateral system are oriented along a set of orthogonal axes, 1) Independent analyses about each of the principal axis shall be performed.
100% 100% Independent
y x
2) Case of no orthogonal axis and IEFaSa (0.2) < 0.35, Independent analyses about any two orthogonal axes are permitted.
permitted 100% 100%
permitted 100% 100%
y x
Independent
Independent
3) Case of no orthogonal axis and IEFaSa (0.2) 0.35,
y 100% together 30% 30% together 100% x
30%
100%
100% 30% together
together
Dynamic Analysis Procedure
Linear Analysis
Non-linear Analysis
Modal Response Spectra
Time History
Time History
Ground motion history Spectral acceleration values compatible with S(T) S(T) defined by the Code
Scaling used in Dynamic Analysis 1) Evaluate the inelastic base shear force V from the Code equation (taking RdR0 into account) 2) Evaluate Ve' from Linear Dynamic Analysis (based on S(T)) 3) Ve (elastic) = Ve'.IE (Dynamic analysis)
4) Vd' (inelastic) =
Ve R d R0
(Dynamic analysis)
5) Compare Vd' to V
Vd ' V
Vd' < V
Vd = V d ' Bx < 1.7 Vd = 0.8V Bx 1.7 Vd = V
6) Redo the dynamic analysis after scaling dynamic excitation by a ratio = Vd/Ve' to obtain the members internal forces. 7) To estimate deflection and interstorey drift, multiply the dynamic analysis results obtained from (6) by the ratio RdR0
*The dynamic analysis should include the effect of accidental torsion
Bx 1.7 ( 0.1Dnx)Fx
Bx < 1.7 center of mass shifted by ( 0.05Dnx)Fx
Time History Analysis
Selection of a Ground Input Acceleration ag(t) Clause 4.1.8.12(3) ag(t) is such that its response spectra is compatible with S(T) Its response spectrum should be equal or exceed the target spectra throughout of the period range of interest.
This can be done by scaling or modifying actual records of similar magnitude and similar distances that contribute to seismic hazard at the site.
Halchuk and Adams (2004) have conducted a deaggregation of seismic hazard for selected Canadian cities.
Halchuk K, S. and Adams, J. Deaggregation of Seismic Hazard for Selected Canadian Cities 13th World Conference on Earthquake Engineering, Vancouver, B. C., Canada, August 1-16,2004, paper No. 2470.
Deaggregation: Dividing the total hazard into contributions based on distance and magnitude.
Halchuk and Adams (2004)
Halchuk and Adams (2004)
Atkinson and Bersnev (1998) have developed various ground motion time histories which are compatible with the uniform hazard spectra used in NBCC 2005. Atkinson, G., and Bersnev, I., (1998), Compatible ground-motion time histories for new national seismic hazard maps, Canadian Journal for Civil Engineering 25: 305-318. Records for both eastern and Western Canada were developed in this study. o For long period structures Selected earthquake records should (a) Have the same deaggregation parameters (D,M) as those specified for S(1.0) of the location. (b) Be scaled to match S (T1), where T1 is the fundamental period of the structure. o For short period structures Selected earthquake records should (a) Have the same deaggregation parameters (D,M) as those specified for S(0.2) of the location. (b) Be scaled to match S(T1), where T1 is the fundamental period of the structure.
CHAPTER 5 SEISMIC ANALYSIS USING COMPUTER MODELING
By: Prof. A. A El Damatty Page: 5-1 .
Two- Dimensional Seismic Analysis Using SAP2000
m3 = 57290.4kg W 690 x 419 W 460 x 286
Given: 2-D moment resistant steel frame with dimensions shown in Figure 1. E = 2x10 kN/m
8
W 460 x 286
4m
m2 = 42967.8 kg W 690 x 419
W 530 x 447
W 530 x 447
4m
m1 = 28645.2 kg
W 610 x 551
10 m
Figure 1
Required: 1- Free vibration: period, & mode shapes. 2- Time history Analysis. 3- Response Spectrum analysis.
W 610 x 551
W 690 x 419
4m
Solution: 1- Select the units from the drop box
2- From File menu select New Model From Template, this will display different modal templates, choose portal frame type and fill the properties of the portal frame: o Number of stories: 3 o Number of bays: 1 o Story height: 4 o Bay width: 10 3- The screen will be titled vertically showing the frame in 2-dimension (x-z plane) & 3- dimension, close the 3 dimension window. 4- Click on Set Elements on the main toolbar & check the labels box
in the joints & frames area to see the joint & frame numbering. 5- Click on joints 1 and 6. From the Assign menu select Joint, from joints select Restraints then choose fixed support from the fast restraints. 6- From the Define menu select Materials, the define material dialog box will show up, highlight STEEL then click Modify / Show Material to check the properties of the material. 7- Verify that modulus of elasticity is 2.0 E+11 then click OK twice. 8- From the Define menu select Frame Sectionsto display the define frame sections dialog box. 9- In the Click to: choose Import I / Wide Flange, another box will display for Section Properties File, from Sap2000n folder in your hard drive choose Cisc.pro (see Fig. 2). Double click on the Cisc.pro and choose from the given sections: W410x100, W610x551, W530x447, &
W460x286. (you can press the Ctrl key on the keyboard to choose all of them at one step). 10- Select the column of the first story (element 1, 4), then from Assign menu select Frame, Section, choose W610x551, make the same for the column of the second, third stories (element 2, 5, & element 3,6). After each selection for beam or column section you can click on the undeformed shape icon to see the element numbering again. 11- From Set Elements click on the frame labels to remove it. 12-The final Figure will be as shown in Fig. 3.
Figure 2
Fig. 2
Fig. 3
13-Select joints 2 and 6. From Assign menu select Joint then select Masses, a dialog will display. In the dialog box of direction 1 write 14322.6. 14-Repeat step 16 for joints 3, 7 and write 21483.9 in the dialog box of direction 1, also joints 4, 8 and write 28645.2 in the dialog box of direction 1. 15-Select Nodes 2 & 6 then from the Assign menu choose Joint, Constraints this will display Constraints box.
16-In the Click to: choose Add Rod, accept the default Rod 1, Click Ok twice to finish this part. 17-Repeat step 17, 18 for Nodes 3 , 7 & Nodes 4 , 8 to make Rod 2 & Rod 3. 18-Form the Define menu select FunctionsResponse Spectrum to display define response spectrum functions box. 19-In the Click to: choose User Spectrum, a dialog box for response spectrum function definition will display. 20-In the Function Name write montreal-est-C, input the values of the periods and spectral accelerations corresponding to the response spectrum of Montraleast with soil type C and damping of 0.05, see Fig. 4. 21-Form the Define menu select Functions Time History to display define time history functions box.
Fig. 4
22- In the Click to: choose Add Function From File, a dialog box for time history function definition will display. 23-In the Function Name write ELCENTRO, click Open File to choose the file having the data for your time history from its location in your computer. 24-Write in the box of number of points per line the number of reading per line, in our example we have 1 per line. Then choose either Function at Equal Period Step of or Time and Function Values depending on how your data file is organized. In our case we choose Time and Function values, see Fig. 5. 25-Form the Define menu select Analysis CasesResponse Spectrum to display define response spectra box. 26-Choose Add New Case, a dialog box for response spectrum case data will display. 27-In the Response Spectrum Case Data box: o Accept the default Spectrum Case Name: ACASE1 o Analysis Case Type: Response Spectrum o Accept the default Excitation Angle: 0 o Accept the default Modal Combination Option, CQC.
Fig. 5
o In the Damping box: 0.05 o Accept the default Directional Combination option, SRSS. o In the Function area select montreal-est-C in U1 direction o For U1 direction write 1 at the scale factor. o Click ADD button. o Click Ok. 28- Choose Add New Case, a dialog box for response spectrum case data will display. 29-In the Time History Case Data box: o Accept the default History Case Name: ACASE2 o For modal damping: click Modify / Show and write 0.05 in Damping for all modes. o For Load type: choose acceleration. o In Number of output time steps box: 1560. o In Output Time Step Size box: 0.02 o In the Load Assignments area select ELCENTRO from Functions for U1 direction o For U1 direction write 1 at the scale factor and leave the arrival time, angle 0. o Click ADD button. 30-Click Ok twice to finish this part. 31-From the Analyze menu select Set Options to display the analysis option dialog box and choose the x-z plane.
Output: We will divide the output file into: 1- Output can be obtained from the screen for the time history of a certain quantity such as: o Time history for a joint displacement, rotation, acceleration or velocity. o Time history for bending moment, shear force or axial force at a specific location of an element. o Time history for base shear force and overturning moment. 2- Output obtained from the 2Dframe.out output file.
Time History 1. Time history for displacement ux of Joint 4: o Click on Joint 4, and then from Display menu click on Show Time History Traces this will display Time History Display Definition. From the choose function part double click on joint 4 at List of Functions, joint 4 will now appear at Plot Function. o The default is drawing displacement ux for this joint, if you want to change click on Define Function icon, another box will display for Time History Functions. o Highlight joint 4, click on Modify/Show TH Function, another box for Time History joint function will appear. Choose the function you want then click Ok twice. o Now Highlight Joint 4 then click the Display icon, another display window will appear as in Fig. 4. You can either plot this figure or
save it by clicking in file the menu and choose Print Graphic or Print Tables To File o The same procedure for joint 4 can be applied for frame element by choosing the required frame element. 2. Time history for Base Shear & Base Moment: o Click on the Define Function icon, the box for Time History Functions will display. o In the Click to: choose Add Base Functions, another window will display for Base functions, Click on Base Shear X and Base Moment Y. Click Ok Twice. o Now Base Shear X and Base Moment Y will be at List of Functions, use the same procedure described for joint 4 to plot the Base Shear X and Base Moment Y.
Fig.6
2Dframe.out output file: In this file one can find information about: 1. Natural frequencies and mode shapes 2. Modal Participation Factors. 3. Mass Participation Factors. 4. Joint Displacements due to SPEC1 Spectra. 5. Frame Element Internal forces due to SPEC1 Spectra. N.B. as shown in Fig. 8, some specific can be selected for printing.
Fig. 8
Three-Dimensional Seismic Analysis

The building model can be developed using many graphical interface structural programs like ETABS. ETABS is a special purpose analysis and design program developed specifically for building structures. An example for a 30storey reinforced concrete building will be illustrated in this section. The building is assumed to be located in Toronto. The total height of the building is 92 m. Modeling Main Steps: 1- Develop the building geometry and axes using CAD or ETABS
2- Define materials properties 3- Define the building elements properties Shell elements are used to model the slabs and the shear walls, and frame elements are used to model the columns and the rigid frames. The stiffness of the pile foundations or shallow foundations are usually accounted for through the use of equivalent vertical, horizontal and rotational springs.
To account for the cracked section properties, the section stiffness is reduced using stiffness modifiers. The reducing factors can be evaluated from the CSAA23.3-04 Clause 21.2.5.2.1 Table 21.1 and Clause 21.2.5.2.2.
4- Construct the Model
5- Preliminary Dynamic Analysis The building is assumed to be located in Toronto. The soil profile is type D. NBCC 2005 spectral values for Toronto Sa (0.2) = 0.26 Sa (0.5) = 0.13 Sa(1.0) = 0.055 Sa(2.0) = 0.015 Using Soil Type D Clause 4.1.8.4.6 S (0.2) = 3.315 m/sec2 S (0.5) = 1.785 m/sec2 S(1.0) = 0.7553 m/sec2 S(2.0) = 0.206 m/sec2 S(4.0) = 0.103 m/sec2 Notice that for dynamic analysis S has units of acceleration. It is obtained by multiply Sa by Fa or Fv and the ground acceleration g ------------ Fa = 1.30 and Fv = 1.40 The Design spectral acceleration values can be determined using NBCC 2005
Conduct Analysis to obtain the following results 1- The building modal participating mass ratios, mode shapes and periods are obtained from the analysis:
Mode 1 2 3 4 5 6 Period 4.235024 3.50052 1.108758 1.067125 0.953012 0.529784 UX 28.7694 23.0926 10.3055 17.4395 1.7659 3.843 UY 40.5759 9.3154 3.987 15.8869 12.1331 0.7729 UZ 0 0 0 0 0 0 S umUX 28.7694 51.862 62.1674 79.6069 81.3728 85.2159 S umUY 40.5759 49.8913 53.8783 69.7652 81.8983 82.6712 S umUZ 0 0 0 0 0 0 RX 58.0654 13.4977 0.0702 26.0066 1.8824 0.0372 RY 40.2014 30.6581 0.2505 27.8469 0.5209 0.179 RZ 0.5623 40.7607 0.1777 34.5346 4.1669 0.95 S umRX 58.0654 71.5631 71.6334 97.6399 99.5223 99.5595 S umRY 40.2014 70.8595 71.11 98.9568 99.4778 99.6568 S umRZ 0.5623 41.3231 41.5007 76.0354 80.2023 81.1524
T = 4.235 sec (Combined mode X-Y - dominant Y ) T= 3.50 sec (Combined mode X-Y-Torsion - dominant X and Torsion) 2- Estimate the periods using NBCC 2005 (hn) Ta (MRF) = Shear Walls Ta = = 2.23 sec = 1.50 sec
Maximum allowable periods according to NBCC 2005 Ta Ta (1.50)(2.23) = 3.345 sec for (MRF) (2.0)(1.50) = 3.0 sec for (Shear walls)
Use Ta = 3.0 sec to evaluate the base shear force based on the NBCC 2005 6- Code calculations Mv = 2.50 (Table 4.1.8.9) Assume moderately ductile shear wall building Rd = 2.0 IE = 1.0 W = 143883 kN =The dead load of the building including the self weight of all elements, the super imposed load, and 0.25 of snow loads. = 0.206 Ro = 1.50
The base shear force is calculated using the NBCC 2005 equation: Clause (4.1.8.11)
and should not be less than
V = 24,700 KN Conduct static analysis using the equivalent static earthquake load to obtain Bx. The static analysis involves the cases of eccentricity. Bx = max/ave max =0.25 m Bx = 2.70 7- Evaluation of spectrum scaling factor and re-analysis ave = 0.09 m
Along X-Direction From initial analysis V'xex = 19,073 kN V'xey = 12,619 kN V'xe = = 22,869 KN
Vxe = V'xe * IE = 22,869 *1.0 = 22,869 KN V'xd = = 7,623 KN
Compare V'xd to V
V'xd V
V'xd < V
Vd = V'xd Regular Vd = 0.8V Irregular Vd = V
Bx 1.7 Vd = V = 24,700 KN Scaling Factor = Vd/V'xe = 1.08 Multiply all values of S(T) by 1.08 and redo the analysis to obtain the members forces.
Repeat these steps along the Y-Direction From initial analysis V'yex = 12,619 kN V'yey = 20,158 kN V'ye = = 23,782 KN
Vye = V'ye * IE = 23,782 *1.0 = 23,782 KN V'yd = = 7,927 KN
Compare V'yd to V
V'yd V
V'yd < V
Vd = V'yd Bx < 1.7 Vd = 0.8V Bx 1.7 Vd = V = 24,700 KN Scaling Factor = Vd/V'ye = 1.04 Multiply all the S(T) by 1.04 and redo the analysis to obtain the members forces. Bx 1.7 Vd = V
Earthquake analysis cases -Specx : response spectrum in X- Direction -Specy : response spectrum in Y- Direction -Specxeccp : response spectrum in X- Direction with eccentricity +0.10 -Specxeccn : response spectrum in X- Direction with eccentricity -0.10 -Specyeccp : response spectrum in Y- Direction with eccentricity +0.10 -Specyeccn : response spectrum in Y- Direction with eccentricity -0.10 -Specxudp : The storey shear force obtained from the case Spec X and applied at a distance of +0.10 Dnx from the shear center.
Dnx = Plan dimension of the building at level x perbendicular to the direction of seismic loading being considered -Specxudn : The storey shear force obtained from the case Spec X and applied at a distance of -0.10 Dnx from the shear center. -Specyudp : The storey shear force obtained from the case Spec Y and applied at a distance of +0.10 Dnx from the shear center. -Specyudn : The storey shear force obtained from the case Spec Y and applied at a distance of -0.10 Dnx from the shear center. Earthquake X-Direction is an envelope to all the X-direction cases. Earthquake Y-Direction is an envelope to all the Y-direction cases 8- Load Combinations Based on NBCC 2005. Clause 4.1.3.2. Case 1 (E1): 1.40 Dead loads Case 2 (E2): (1.25 Dead Loads or 0.90 Dead loads) + 1.50 Live Load + (0.25 Snow or 0.40 Wind Load) Case 3 (E3): (1.25 Dead loads or 0.90 Dead loads) + 1.50 Snow + (0.50 live loads or 0.40 Wind) Case 4 (E4): (1.25 Dead loads or 0.90 Dead loads) + 1.40 Wind X-Direction + (0.50 live loads or 0.40 Snow) Case 5 (E5): (1.25 Dead loads or 0.90 Dead loads) + 1.40 Wind Y-Direction + (0.50 live loads or 0.40 Snow)
Earth Quake Combinations Case 6 (E6): 1.0 Dead loads +( 1.0 Earthquake X-Direction + 0.30 Earthquake Y-Direction) + 0.50 live loads + 0.25 Snow Case 7 (E7): 1.0 Dead loads + (1.0 Earthquake Y-Direction + 0.30 Earthquake X-Direction) + 0.50 live loads + 0.25 Snow 9- Results a) Shear wall Design -Magnify the moments obtained from the final analysis with the P- factor U2
U2 = 1 ( 1 f pf hVf
f = inter storey drift due to factored lateral load. pf = factored axial loads in the story. Vf = factored total story shear. h = story height. - Magnify the shear forces obtained from the final analysis by a factor = (Rd). (Ro) b) Columns Design - Magnify the moments obtained from the final analysis with the P- factor U2 - Magnify the shear forces obtained from the final analysis by a factor = (Rd). (Ro)
c) Frame Girders Design (Ductile MRF) - Magnify the moments obtained from final analysis by a factor = (
Rd R0 ( Moderatelyductilewalls ) Rd ' R0' ( DuctileMRF )
). (U2)
) ------- Reduction factors for Ductile MRF
d) Deflection and interstorey drift - Magnify the deflection and interstorey forces obtained from the final analysis by a factor = (Rd). (Ro)
CHAPTER 6 SEISMIC DESIGN OF STEEL STRUCTURES
By: Prof. A A. El Damatty
Page: 6-1 .
6-1 Concept of Capacity Design Capacity design was developed in the late 1960s in New Zealand as an approach to resist the effect of severe earthquakes. The concept of capacity design is based on the following: Inelastic action is unavoidable during severe earthquakes. The designer dictates where inelastic response should occur. Zones of possible inelastic actions (called plastic zones) are selected to be regions where large plastic deformation can develop without significant loss of strength (regions which possess enough ductility) These regions are detailed to suppress premature failure modes, such as local buckling or member instability in the case of steel structures, or shear failure in the case of concrete structures. The surrounding members are designed in such a way that their capacities are greater than that the values corresponding to the maximum capacity of the plastic zone. These members are intended to remain elastic during an earthquake.
Page: 6-2 .
An illustration for the concept of plastic design is shown in Fig. 6.1. This example involves a cantilever beam of total length L consisting of a brittle segment (such as a fiber composite material) of length a and a ductile steel segment of length b.
P
brittle material
steel
a L (L/b) Mp (steel) Mp (steel)
Fig. 6.1 Capacity design approach would aim at making the brittle material stronger than needed to ensure that plastic hinging occurs first in the steel segment of the cantilever. Therefore, the moment capacity of the brittle material should satisfy: Mbrittle L/b Mp (steel) is a number greater than 1.0 to account for the possible reserve strength of the cantilever beyond its nominal yield strength.
Page: 6-3 .
Steel is in general a ductile material However, brittle failure of a steel member might result due to any of the following effects: (a) Overall buckling (b) Local buckling (c) Failure of connections
Stress-strain curve for a typical steel coupon
Two examples emphasizing the capacity design of moment resisting and concentric bracing ductile frame structures are presented in these notes.
Page: 6-4 .
6-2 Strength Design Before performing a seismic (capacity) design of a lateral resisting system, the structure is analyzed under the effects of various load combinations (DL, LL, W, S, E) defined by the NBCC and the maximum factored forces and moments resulting from various load combinations are determined. Members are designed to satisfy the strength and serviceability criteria defined in the design Code.
The members and connections will be then re-designed according to the capacity design provisions (clauses 27 ..)
6-3 Types of Resistance Ultimate (Mr) r = 0.9 Nominal (Mp) p = 1.0 Probable (Mp) p = 1.1 Ry
Page: 6-5 .
Moment Resisting Frames

(Energy dissipates through yielding of beams)
Page: 6-6 .
6-4 Seismic Design Provisions for Ductile Moment Resisting Frames, Rd = 5.0, R0 = 1.5 (S16-01) No Limit on the height of the building Plastic hinges can form only in beams (flexural members). Plastic hinges in columns are permitted only at the base, except in single storey buildings. Beams (27.2.2) o Class 1 sections o Laterally braced according to the requirements of clause 13.7 (b)
Lr 17250 + 15500k = ry Fy
k is the ratio between the smaller factored moment to the larger factored moment at opposite ends of the beam; k is positive for double curvature, and shall be based on the both gravity and seismic loads. Bending moment due to seismic loads may be taken as linear variable from maximum value to zero at both ends of the beam.
Page: 6-7 .
Columns (27.2.3.2) Class 1 or Class 2 Mrc' [1.1 Ry Mpb + Vh (x+ dc/2)] where: Mrc' = Sum of the column factored flexural resistances at the intersection of the beam and column Mpb = Nominal plastic moment resistance of the beam Vh = Shear acting at plastic hinge when 1.18 Ry Mpb is reached x = Distance from centre of the beam plastic hinge to the column face
Cf
Mrc' = 1.18 Mpc [1-
C y
] Mpc
Mpc = nominal plastic moment resistance of the column Cf = results from summation of Vh at the considered level. More requirements are specified in Clause 27.2.3.1 if the plastic hinge develops in the column. ( 27.2.4)
Panel Zone
When plastic hinges form in adjacent beams, the panel zone shall resist forces arising from beam moments of: [1.1 Ry Mpb + Vh (x+ dc/2)]
Page: 6-8 .
Resistance of Panel Zone

2 3bc t c
(a) Vr = 0.55 dc w Fyc [1+
'
dc db w
'
] 0.66 dc w' Fyc
The above equation can be used if the following conditions are satisfied 1) When IEFaSa(0.2) 0.55, the sum of the panel zone depth and width divided by the panel zone thickness does not exceed 90. 2) Satisfy the width-to-thickness limit of Clause 13.4.1.1(a) 3) Doubler plates are groove or fillet welded. Doubler plate can be included in calculating the thickness to width ratio if it is connected to the column web near the centre of the panel.
Page: 6-9 .
Moderately Ductile Moment-Resisting Frames (Rd = 3.5, R0 = 1.5)
Same provisions as ductile Moment Resisting Frames except that: (i) The beams shall be Class 1 or 2 sections (ii) The bracing beam shall meet the requirements of Clause 1.7 (a) instead of 13.7 (b)
Lcr 25000 + 15000k = ry Fy
Page: 6-10 .
Evaluation of Forces in Panel Zone
M3 M1 V1 C1 V4 V3
C3 Vh V2 M2 M4 C4 V1 V3 h3/2 h4/2 V4 l1/2 l2/2 0.95 dc
C2
0.95 db
V2
Equilibrium Equations
M3 M1 M2 M4 + =0 (1 1 ) (1 2 ) (1 3 ) (1 4 )
(1)
Vh = [ M1 + M2 -
3M 3 M 4 4 ]/0.95db (1 3 ) (1 4 )
(2)
1 =
0.95d c 0.95d c 0.95d b 0.95d b , 2 = , 3 = , 4 = l1 l2 h3 h4
Page: 6-11 .
Concentric Braced Frames (Energy dissipates through yielding of bracing members)
Page: 6-12 .
Type Tension-Compression
Moderately Ductile If IEFaSa(0.2) < 0.35 No limit Otherwise H 40 m
Limited Ductility If IEFaSa(0.2) < 0.35 No limit Otherwise H 60 m
(a)
When H > 32 m, seismic When H > 48 m, seismic force shall be increased force shall be increased by 3% per meter above by 2% per meter above 32 m 48 m
Chevron
(b)
Same as (a)
Same as (a)
Tension only
20 m
40 m
If H > 16 m, seismic force shall be increased by 3% per meter of height above 16 m
If H > 32 m, seismic force shall be increased by 3% per meter of height above 32 m
K-truss N.A. N.A.
Page: 6-13 .
Moderately ductile concentrically braced frames (Rd = 3.0, R0 =1.3) Diagonal Bracing Limit on the width-to-thickness ratio o For IEFaSa(0.2) 0.35 Sections Rectangular section and HSS Legs of angle and flanges of Channel Other elements Class 1 Class 2 145/Fy 170/Fy 300/Fy Class 1 kL/r 100 kL/r = 200
o For IEFaSa(0.2) < 0.35 HSS => Class 1 Other sections => Class 2 Legs of angle = 170/Fy
Page: 6-14 .
Bracing Connection Factored resistance of the bracing connections shall be equal or exceed both AgRyFy and 1.2 the times the nominal compressive resistance of the bracing members (1.2Cp) The tensile force should not exceed the combined effect of gravity in the bracing and the effect of seismic loads corresponding to RdR0 = 1.3 (Redistribution due to buckling should be taken into account when calculating the seismic loads corresponding to RdR0 = 1.3)
Beams The capacity design of the beam should consider two conditions: 1- Compression bracing = 1.2Cp Tension bracing = AgRyFy
2- Compression bracing buckled and retain a force = 0.2 AgRyFy Tension bracing = AgRyFy
Columns Factored resistance of columns shall equal or exceed the effects of gravity and the brace forces described above for connections design Class 1 or 2 Bending resistance 0.2 ZFy
Page: 6-15 .
Seismic Design of a Moderately Ductile MRF and Moderately Ductile braced frame
1. Eight storeys 2. Located Montreal, Quebec. 3. Loading according to NBCC 05
Moment Resisting Frame
Braced Frame
Moment Resisting Frame
6x6000mms Braced Frame
6000
6000
8000
6000
6000
Plan of Example Building Loads: Dead Load: (NBCC 4.1.5) 1. self weight = 0.10 (thickness of slab) * 25 (Concrete unit weight)= 2.5 kN/m2 2. partition = 1.0 kN/m2 3. mechanical service = 0.5 kN/m2 4. flooring = 0.5 kN/m2 5. structural elements (beams & columns) = 0.25 kN/m2 D.L. = = 4.75 kN/m2 Live load: (NBCC 4.1.6) 1. L.L. on typical office floor = 2.4 kN/m2
Page: 6-16 .
Roof load 1. self weight = 0.10 (thickness of slab) x 25 (Concrete unit weight) = 2.5kN/m2 2. mechanical loading = 1.6 kN/m2 3. roofing = 0.5 kN/m2 4. structural elements (beams & columns) = 0.25 kN/m2 D.L. = = 4.85 kN/m2 5. snow load = 2.2 kN/m2
Analysis under Factored Loads The following load cases were considered: i. 1.25 DL + 1.5 LL + 0.5 Snow + P- effect ii. 1.25 DL + 0.5 LL + 1.5 Snow + P- effect iii. DL + E (No eccen.) + 0.5LL + 0.25SL + P- effect (both directions) iv. DL + E (Eccen. e) + 0.5LL + 0.25SL + P- effect (both directions)
Mass center
Earthquake load (N-S) (No Eccen.)
Earthquake load (N-S) (+e Eccen.)
Earthquake load (N-S) (-e Eccen.)
Earthquake load (E-W) (No Eccen.)
Earthquake load (E-W) (+e Eccen.)
Earthquake load (E-W) (-e Eccen.)
Page: 6-17 .
Seismic Loads Calculations Using the Equivalent Static Approach (NBCC 05) The total weight of the floors are given in the table below, Floor Height (m) 8 7 6 5 4 3 2 1 29.7 26.1 22.5 18.9 15.3 11.7 8.1 4.5
Weight ( kN) 6220.8 5472 5472 5472 5472 5472 5472 5472 W = 44524.8 kN
Due to symmetry, the centre of mass in each storey coincides with the centre of geometry. Lateral Earthquake force V V = S(Ta) . Mv . IE . W / (Rd . Ro) S(2.0) . Mv . IE . W / (Rd . Ro) Time period Ta North-south direction: Steel moment fames, T = 0.085 (hn)3/4 = 1.08 sec. East-west direction: Braced frames, T = 0.025 (hn) = 0.74 sec
NBCC 4.1.8.11 NBCC 4.1.8.11
For Montral, the 5% damped spectral response acceleration ratios, Sa(T), are provided in Table-C2 of NBCC. Sa(0.2) 0.69 Sa(0.5) 0.34 Sa(1.0) 0.14 Sa(2.0) 0.048 PGA 0.14
Site class: C (very dense soil and soft rock) is assumed. From table 4.1.8.4.b, the value of acceleration-based site coefficient Fa = 1
Page: 6-18 .
From table 4.1.8.4.c, the value of the velocity based coefficient Fv = 1.0 The design spectral acceleration value S(T) T =< 0.2 sec.; S(T) = Fa . Sa(0.2) = 1 * 0.69 = 0.69 T = 0.5 sec.; S(T) = Fv . Sa(0.5) = 1 * 0.34 = 0.34 take the smallest = Fa . Sa(0.2) = 1 * 0.69 = 0.69 T = 1.0 sec.; S(T) = Fv . Sa(1.0) = 1 * 0.14 = 0.14 T = 2.0 sec.; S(T) = Fv . Sa(2.0) / 2 = 1 * 0.048 = 0.024 Using linear interpolation North-south direction S(1.08) = 0.13 East-west direction S(0.74) = 0.24 From table 4.1.8.11., Sa(0.2)/Sa(2.0) = 14.375 > 8.0 , so higher mode factor Mv = 1, o Base overturning reduction factor Jy = 0.884 (N-S direction) o Base overturning reduction factor Jx = 0.92 (E-W direction) Normal importance category for this office building is assumed, IE = 1.0. From table 4.1.8.9. North-south direction, moderately ductile moment-resisting frames, Rd = 3.5, Ro = 1.5, East-west direction, moderately ductile concentrically braced fames, Rd = 3.0, Ro = 1.3 Lateral Earthquake force V North-south direction V = S(Ta).Mv.IE.W/(Rd.Ro) = 1102.6 kN S(2.0). Mv. IE.W / (Rd.Ro) East-west direction V = S(Ta).Mv.IE.W/(Rd.Ro) = 2740.2 kN S(2.0). Mv. IE.W / (Rd.Ro) Concentrated force at the top of the building Ft = 0.07.Ta .V North-south direction = 83.4 kN East-west direction = 142.4 kN Wi . hi = 770808 kN Fi = (V - Ft). Wi . hi / ( Wi . hi)
Page: 6-19 .
Floor Level 8 7 6 5 4 3 2 1
North-south direction Force (kN) 244.2 + 83.4 188.8 162.7 136.7 110.6 84.6 58.5 32.5
East-west direction Force (kN) 622.6 + 142.4 481.2 414.9 348.5 282.1 215.7 149.3 82.9
The structure shall be designed to resist overturning effects caused by the earthquake forces determined in Sentence (6) 4.1.8.11 and the overturning moment at level x, M, shall be determined using the following equation: Mx = Jx Fi ( hi hx) where Jx = 1.0 for hx 0.6 hn, and Jx = J + (1-J) (hx / 0.6 hn) for hx < 0.6 hn, and (N S) direction For floors 5 to 8, Jy = 1.0 Floor 4 : Jy = 0.884 + (1-0. 884) (15.3 / 17.82) Floor 3 : Jy = 0. 884 + (1-0. 884) (11.7 / 17.82) Floor 2 : Jy = 0. 884 + (1-0. 884) (8.1 / 17.82) Floor 1 : Jy = 0. 884 + (1-0. 884) (4.5 / 17.82) Ground : Jy = 0. 884 + (1-0. 884) (0 / 17.82) (E W) direction For floors 5 to 8, Jx = 1.0 Floor 4 : Jx = 0.92 + (1-0.92) (15.3 / 17.82) Floor 3 : Jx = 0.92 + (1-0.92) (11.7 / 17.82) Floor 2 : Jx = 0.92 + (1-0.92) (8.1 / 17.82) Floor 1 : Jx = 0.92 + (1-0.92) (4.5 / 17.82) Ground : Jx = 0.92 + (1-0.92) (0 / 17.82)
= 0.983 = 0.960 = 0.937 = 0.913 = 0.884
= 0.988 = 0.973 = 0.956 = 0.940 = 0.920
Page: 6-20 .
Overturning Floor moment (kN.m) (N S) direction 0 1179.7 3039.3 5484.8 8422.5 11758.7 15399.7 19251.5 24212.9
Overturning moment (kN.m) (E W) direction 0 2754.1 7240.8 13221.3 20456.4 28707.2 37734.7 47299.9 59629.8
Reduced overturning moment (kN.m) (N S) direction 0 1179.7 3039.3 5484.8 8284.4 11290.3 14425.3 17582.3 21404.2
Reduced overturning moment (kN.m) (E W) direction 0 2754.1 7240.8 13221.3 20224.9 27918.4 36088.1 44471.5 54859.5
8 7 6 5 4 3 2 1 Ground
Evaluation of Seismic Loads Using Response Spectrum Analysis The natural frequencies and mode shapes of the building were evaluated using the program SAP 2000 and the results are shown below: Ve' (N-S) = 3868kN Ve' (E-W) = 5150 kN Ve (N-S) = Ve'.IE = 3868 kN Ve (E-W) = Ve'.IE = 5150 kN Vd' (N-S) = Vd' (E-W) =
Ve R d R0 Ve R d R0
= 737 kN < 1103kN = 1321 kN < 2740 kN
Page: 6-21 .
First mode shape (Period =1.46075 sec)
Second mode shape (Period = 1.32939 sec)
Page: 6-22 .
Modal Participating Mass Ratios OutputCase StepType StepNum Text modal modal modal modal modal modal modal modal modal modal modal modal Text Mode Mode Mode Mode Mode Mode Mode Mode Mode Mode Mode Mode Unitless 1 2 3 4 5 6 7 8 9 10 11 12 Period Sec 1.460747 1.329386 1.269189 0.534707 0.463849 0.450689 0.300466 0.24875 0.248726 0.241159 0.233314 UX Unitless 0.68881 1.543E-08 0.00006013 0.14181 0.00001184 6.655E-09 0.0645 7.493E-10 1.234E-09 3.006E-07 6.89E-08 UY Unitless 9.718E-14 0.70367 0.00024 8.891E-13 0.00014 0.2019 2.248E-13 4.123E-10 0.04821 2.695E-11 2.802E-10 UZ Unitless 2.785E-09 5.495E-12 2.891E-09 3.087E-11 2.35E-10 1.728E-09 5.162E-08 5.093E-08 0.09393 8.926E-08 0.00004203 0.0009
0.266385 0.000004504 0.000001035
Page: 6-23 .
3-D view for the structure Torsional Sensitivity Level 8 (East-west direction) Torsional sensitivity can be determined by calculating the ratio, Bx = max / ave max = 0.01456 m ave = (0.01456 + 0.01290) / 2 = 0.01373 Bx = 0.01456/ 0.01373 = 1.06045 Level 8 (North-south direction) Torsional sensitivity can be determined by calculating the ratio, Bx = max / ave max = 0.02080 m ave = (0.02080+ 0.01807) / 2 = 0.019435 Bx = 0.02080 / 0.019435 = 1.0702
Page: 6-24 .
Interstorey Drift East-west direction The maximum interstorey drift = 0.00228 * Rd Ro = 0.008268 m < 0.025hs (floor height) < 0.09 m North-south direction The maximum interstorey drift = 0.00724 * Rd Ro = 0.003801m < 0.025hs (floor height) < 0.09 m
3-D view for the structure
Page: 6-25 .
Strength design according to CISC 95 North-south direction: Two moment resisting frames along the outer edges
Elevation of Moment Resisting Frame
Page: 6-26 .
East-west direction: Two braced frames

F8
S8
F7
S7
F6
S6
F5
S5
F4
S4
F3
S3
F2
S2
F1
S1
Elevation of Braced Frame
Page: 6-27 .
Capacity design of the Rigid Frames An interior joint of the moment resisting frame is considered to illustrate the concept of capacity design. (Floor F4)
W 690 x 140
W 690 x 140
Level 4
W 760 x 147
d b t w Beam length = L1 = L2 = 6000 mms Column height = h1 = h2 = 3600 mms 1 = 2 = 0.95dc/L1 = 0.1192 3 = 4 = 0.95db/ h1 = 0.1805 Beams W 690x140
Beam W690x140 684 254 18.9 12.4
Column W760x147 753 265 17.0 13.2
Mp
Mp/2
25000 (15000)(0.5) Lcr = 55.65) > ( = 50) Not good 350 ry
Page: 6-28 .
Assume the cross beam will be repeated every 2.0 m
Mp
2Mp/3
Mp/3
25000 (15000)(0.667) Lcr = 37.105) < ( = 42.8) OK 350 ry
Mr1 = Mr2 = MP = 1410 kN.m; Mp1 = 1410/0.9 = 1566.6 kN Mp1 = Mp2 = (1.1). (1.1). 1410 / 0.9 = 1895 kN.m Assume beam reach their ultimate strengths, Mrc 1.1. Ry . Mpb + Vh (x + dc/2)
w = 35.7 kN/m
dc/2
x L
Lh
Example for F4 Vh = VhE + VG VhE = 2Mpb' /Lh Lh = L-2S S = x + dc/2 VG = wL/2 => w = (DL +0.5 LL) x = db / 2 for welded connection
Page: 6-29 .
Example for F4 (see joint location on Page 6-25) S = (x + dc/2) = 684/2 + 753/2 = 718.5 mm Lh = 6000 (2) (718.5) = 4563 mm VhE = (2) (1895) / 4.563 = 830 kN w = [4.75 + (0.5)(2.4)](6) = 35.7 kN/m VG = (35.7). (6) / 2.0 = 107 kN
305 kN 107 107 107 107 kN 107 107 376 529 1228 107 107 466 1083 107 107 830 1895 755 kN.m 918
Unbalanced moment = (1.1 Ry . Mp + Vh (x + dc/2) ) = (1895) . (2) + (830) . (0.7185) . (2) = 4982 kN.m Mrc should be greater than 4982 / 2 = 2491 kN.m Mrc = (1.18). ( ) . Mpc [ 1 Cf Cy
Take Cf = (107)(10) = 1070 kN . Mpc = 1580 kN.m Cy = A . fy = (18700 ) x (350) / 1000 = 6545 kN Mrc = (1.18). (1580). [1 1070 ] = 1830 kN.m < 2491 kN.m Not good (0.9)(6545)
Page: 6-30 .
Try W920 x 381 . Mpc = 5280 kN.m Cy = 17010 kN Mrc = (1.18). (5280). [1 1070 ] = 5794 kN.m > 2491 kN.m OK (0.9)(17010)
Panel Zone
M3 M1 M2 M [ + 4 4 ] = 0 .0 1 1 1 2 1 3 1 4
Vh = [ M 1 + M 2
M 3 . 3 M 4 . 4 ] / 0.95d b 1 3 1 4
(0.95)(951) = 0.20 6000 (0.95)(684) = 0.18 3600
1 = 2 = 0.95dc/L1 = 3 = 4 = 0.95db/ h1 =
2M 3 ( 2491)(2) = M3 = M4 = (1 0.20) (1 0.18)
2552 kN.m
Vh = [ 2491 + 2491
0.18(2552 ) 0.18( 2552 ) ] /[(0.95).(0.684)] = 5943 kN 1 0.18 1 0.18
(3).(310)(43.6) 2 ] = 4468 < 5943 Not good Vr = (0.55).(0.9).(951).(24.4).(0.35)[1 + (951) x(684)(24.4)
Try: W 920 x 585 d b t w W920x585 960 427 55.9 31
(3).(427)(55.9) 2 ] = 6170 > 5943 ok Vr = (0.55).(0.9).(960).(31).(0.35)[1 + (960) x(684)(50)
Vr should not exceed (0.66dcw'Fy), i.e. should not exceed
Page: 6-31 .
Vr
(0.66) (960) (0.9) (31) (0.35)
6188 kN
Check clause 13.4.1.1 a = 960 (2) .(55.9) = 848.2 h = [ 684 (2) (18.9)] = 647
h / w = 647 / 31 = 20.87 a / h = 848.2 / 647 = 1.31 > 1.0 kv = 5.34 + (4) / (a/h)2 = 7.57 OK
439
kv = 64.9 Fy
k h 439 v w Fy
Fs = 0.66 Fy
OK
Page: 6-32 .
Revision of braced frame design for ductility Consider S8 HSS 178 x 178 x 6 L = 5380 mm b = 177.8 mm t = 6.35 mm r = 69.6 mm A = 4250 mm2 y = 350 MPa Cr = 789 kN Tr = 1339 kN Cf = 248 kN Tf = 99 kN Slenderness ratio (Clause 27.5.3.1) kL/r = 5380/69.6 = 77.3 < 200 Wall slenderness (Clause 27.5.3.2) (b 4t) / t = Limit =
330
177.8 (4).(6.35) = 24.0 6.35
ok
= 17.63 not good
The brace does not satisfy local buckling limit on wall slenderness Revised section HSS 152x152x10 A = 5210 mm2 b = 152.4 mm2 t = 9.53 mm2 r = 57.6 mm2 Tr = 1641 kN
Page: 6-33 .
Cr = 916 kN kL/r = 93.4 < 200 ok (b 4t) / t = 11.99 <17.63 ok Clause (27.5.4.2), connection between the brace and the column should be designed to resist a force Fc which is the largest of: a) Tu = Ag. Ry. y = (1.1). (5210) .(350) = 2005 kN b) Cu = 1.2C p =
1 .2 .(916) = 1222 kN ( 0 .9 )
Fc should not exceed the effect of gravity load + effect of seismic loads with Rd. Ro = 1.3 Fc should be less than (65 + Fc 1151 kN
(362).(3).(1.3) ) (1.3)
Connection should be designed to resist a force = 1151 kN
Beam design Verification The capacity design for the beam should consider two conditions: 1. Compression bracing = 1.2 Cp Tension bracing = Ag. Ry. Fy 2. Compression bracing buckled and retain a force = 0.2 Ag. Ry .Fy Tension bracing = Ag. Ry .Fy The beam need not to resist load effects exceeding to Rd. Ro = 1.3
Page: 6-34 .
Capacity design of strut beam

HSS 152x152x10 After capacity design
T7b F6 C6b p
C7b T6b
HSS 203x203x13 After capacity design
Example beam of level 7 T7b = Ag . Ry . Fy = (5210). (1.1). (0.35) = 2005 kN C7b = 0.2 Ag . Ry . Fy = 401 kN T6b = (9260)(1.1)(0.35) = 3565 kN C6b = 713 kN - P = (T6b C7b) cos P = - (3565 - 401) 0.743 = -2350 kN Mf = (wDL + 0.5 wLL). L2 / (8) = 272 kN.m Try W 460x 89 Cr = 3591 kN, Mr = 624 kN.m (laterally supported)
2350 (0.85) 272 P (0.85) M f + = + = 1.0 ok 3591 624 Cr Mr
Page: 6-35 .
Intersection beam Example beam at F7

HSS 152x152x10 C8b F7 p T8b p
C7b
T7b HSS 152x152x10
2.0 m
T8 = T7 = Ag . Ry . Fy = 2005 kN C8 = C7 = 1.2 Cp = 1222 kN (effective length 5.38 m) P = 0.5 (Ti + Ci) cosi - 0.5 ( Ti+1 + Ci+1) cosi+1 , = 0.75
P = 0.5 (2005 + 1222) (0.743) - (0.75).(0.5).(2005+1222).(0.743) = 300 kN Mf = MDL + 0.5 MLL
Page: 6-36 .
Try W 310x 24, Mf = 45.7 kN.m ( From SAP2000 analysis) Cr = Tr = 958 kN, Mr = 78.3 kN.m ( 2.0 m unbraced length)
P Mf 396 45.7 + = + = 0.99 ok Cr M r 958 78.3
Capacity design of column Example storey (6)
HSS 152x152x10 HSS 152x152x10
2005 1222
F8
F6
HSS 203x203x13 HSS 203x203x13 3565 3040
F4
HSS 203x203x13 3565
(1) Forces Based on Brace capacities
T = Ag Fy Ry C = 1.2 Cp
Page: 6-37 .
2 Cr or Tu Sin 1341 817
3 SRSS
4 Earthquake (RdRo) = 1.3 358
5 EQ (Smallest of 3 and 4) 358
6 D+0.5L+P
7 EQ + D +0.5L+P 764
2005 1222
1341
406
4543
3565 3040
2080
2080
596
2676
2385 2033 7273 2385 4660 4660 1003 5663
3565
SRSS combines the maximum brace contribution at any level above the considered level with the square root of the sum of the squares of all other brace contributions about that level. Example force at storey S4: SRSS = 2033 + 2385 +
( 2385 ) 2 + (817 ) 2 + (1341 ) 2 = 7273 kN
The column at this level should be designed to resist an axial load = 5663 kN resulting from the combination of the EQ and the (DL+0.5LL+P). Column should be Class 1 or 2 and should resist Mf = 0.2 Mp (S16-01-27.5.5.2)
CHAPTER 7 SEISMIC DESIGN OF REINFORCED CONCRETE STRUCTURES
Page: 7-1 By: Prof. A. A. El Damatty____
MOMENT RESISTING FRAMES (MRF)
Ductile MRF
=> Rd = 4.0 R0 = 1.7
Moderately Ductile MRF
=> Rd = 2.0 R0 = 1.4
Conventional MRF
=> Rd = 1.5 R0 = 1.3
Factored, Nominal and Probable Resistances: Mr = Factored Flexural Resistance c = 0.65 / s = 0.85
Mn = Nominal Flexural Resistance c = 1.00 / s = 1.00
Mp = Probable Flexural Resistance c = 1.00 / s = 1.25
Mn Mr
1.20
For Beams,
Mp Mr
1.47
For Columns,
Mp Mr
1.57
Seismic Provisions for Ductile MRF
Beams -Ductile Flexural Members
(21.3.1.1)
- Axial compressive Force - Lclear
Ag f 'c 10
4d
b d
0.3
-b
250 mm
-b
width of supporting member + 2 x 3/4d
d max.
b
d max. Plan View
Beams (cont..)
M-ve joint
M+ve joint
min = 1.4bwd/fy max = 0.025 min. 2 continuous bars (21.3.2.1)
M + ve jo int
M ve jo int 2
+ ve M any sec tion
M jo int 4
ve M any sec tion
M jo int 4
(21.3.2.2)
Beams (cont..)
50 mm
S d/4 8 db 24dst. 300 mm
2d S d/2 S d/4 100 mm
*Lap splices shall not be used:

-within the joints -within a distance 2d from the face of the joint -within a distance d from any plastic hinge caused by inelastic lateral displacements (21.3.2.3)
Beams (cont..)
Hoops Detailing Cross -tie
All extensions > 6db 60 mm
Beams (cont..) Shear Strength Requirement: Vr Vf Vr The smaller of { Vp { V determined from analysis with RdR0 = 1
Vp => shear corresponding to the probable moment resistance at the forces of the joints while the member is loaded with the tributary transverse load
Vr can be calculated using Clause 11 with = 0 and = 450 (Vr = Vc + Vs) Vc = cf'c bwdv Vs = sAvfydvcot / S Av = area of shear reinforcement within a distance s dv = effective shear depth = maximum of 0.9d or 0.72h = factor to account for density of concrete = angle of inclination of compressive stress with longitudinal bar = factor to account for the shear resistance of cracked concrete
Columns
(21.4.1)
-Ductile Frame members subjected to Flexural and Axial Loads - width -b/h 0.4 300 mm
-P >
Ag f 'c 10
M rct
M nbL
(include slab RFT)
M nbR
(include slab RFT)
M rcb
Mrc 1.1Mnb Strong Column Weak Beam
Columns (cont..) *Longitudinal RFT: 0.01
0.06
ld
1.3ld
Tension Lap splice within the centre half of the member length
Tension Lap splice at any section
Columns (cont..) Transverse RFT: Ash 0.2kn kp (Ag/Ach) (f'c/fyh) s hc , Ash 0.09 ( f'c/fyh) s hc * Stirrups: Ach = Core area hc = dimension of concrete core perpendicular to the direction of hoop bars kn =
nl nl 2
Pf P0
kp =
l0
Spacing > 6 db 150 mm
Spacing > b 100+(350-hx)/3 6db
IF Pf 0.5c f'c Ag
l0 < 1.5 h
clear height /6
IF Pf > 0.5c f'c Ag
l0 < 2h
clear height /6
Column Cross sections:
hx
hx
hx
(21.4.4.4)
hx > max. of 200 mm core dimension/3
hx = Distance between laterally supported
Shear Strength Requirements
K2
+ ve pr
M ve pr
M colt =
k1 ( M +ve + M ve ) pr pr k1 + k 2
lu K1
M colb
Vcol =
M colt + M colb lu
Vrcol Vf
Vrcol The smaller of { Vcol
{ V determined from analysis using RdR0 = 1
Vr can be calculated using Clause 11 with 0.1 and 450
Beam-Column Joints: Column stirrups should continue through the joint except if the joint is confined by Structural members as shown;
b < 3/4h < 3/4b
Applied shear on Joints:
Vcolumn
1.25 As fy
Vfb = 1.25 Asfy - Vcolumn
Vcolumn
C2 =1.25As2fy 1.25As2fy C1 =1.25As1fy
Vfb = 1.25 As1fy + 1.25As2fy - Vcolumn
Allowable shear in beam-column joints (a) Continued joints: 2.2 cf'c AJ
(21.5.4.1)
(b) Joints confined on 3 faces or two opposite faces: 1.6 cf'c AJ
or
(c) All others: 1.3 cf'c AJ
AJ is lesser of Ag of column or bw hcol
ldh
12db
(21.5.5.2)
* ldh = max of { 8db {150 mm { 0.2
fy f c'
db
* without the standard 900 hook => use 2.5ldh if depth of concrete below bar 300 mm => use 3.5 ldh otherwise (21.5.5.4)
Seismic Design of a Ductile Moment Resisting Frame
2.20
3.00 7.40 m
2.20 7.40 m 7.40 m
Analysis of the Moment Resisting Frame (MRF) was carried out using plane frame finite element analysis program. Three load cases were considered: o Dead Load (D) o Live Load (L) o Earthquake Load (E) To account for reduction of member stiffness due to concrete cracking (Clause 21.2.5.2.1) o IBeam = 0.4 Ig o IColumn=
fs 0.5 + 0.6* * Ig 1.0 f ' c * Ag (0.6 0.7 Ig)
Design of the First Story Beam

A a b B c
C L
B.M. values are recorded in Table 1
Table 1: Bending moment along the beam [kN.m] (no redistribution) AB a b BA BC D -117.45 113.43 88.26 -174.40 -159.10 L -55.77 54.28 42.10 -83.45 -76.02 E () 120.60 48.28 48.13 119.00 119.50 Load Cases 1.25D+1.5L -230.47 223.21 173.48 -343.18 -312.91 1.0D+1.0E 3.15 161.71 136.39 -55.40 -39.60 1.0D-1.0E -238.05 65.15 40.13 -293.40 -278.60 1.0D+0.5L+1.0E -24.74 188.85 157.44 -97.13 -77.61 1.0D+0.5L-1.0E -265.94 92.29 61.18 -335.13 -316.61 Design moments Max. Moment 3.15 223.21 173.48 -55.40 -39.60 Min. Moment -265.94 65.15 40.13 -343.18 -316.61
c 85.10 40.58 48.49 167.25 133.59 36.61 153.88 56.90 167.25 36.61
Negative moments are calculated at the face of the columns - For better, more economic design; moment distribution will be applied to reduce ve B.M. - For a ductile MRF, CSA allows for a maximum redistribution of 20% (clause 9.4.2) - For simplicity only D & L B.M. will be redistributed B.M. values after redistribution are recorded in Table 2
Table 2: Bending moment along the beam [kN.m] (after redistribution) AB a b BA BC D -93.96 140.35 117.62 -139.52 -127.28 L -44.62 67.11 59.16 -66.76 -60.82 E () 120.60 48.28 48.13 119.00 119.50 Load Cases 1.25D+1.5L -184.37 276.10 235.76 -274.54 -250.32 1.0D+1.0E 26.64 188.63 165.75 -20.52 -7.78 1.0D-1.0E -214.56 92.07 69.49 -258.52 -246.78 1.0D+0.5L+1.0E 4.33 222.18 195.33 -53.90 -38.19 1.0D+0.5L-1.0E -236.87 125.62 99.07 -291.90 -277.19 Design moments Max. Moment 26.64 276.10 235.76 -20.52 -7.78 Min. Moment -236.87 92.07 69.49 -291.90 -277.19
c 116.92 55.78 48.49 229.83 165.41 68.43 193.30 96.32 229.83 68.43
Design for flexural reinforcement

db max = lj 600 = = 25 24 24
(Clause 21.5.5.6)
Maximum bar size that can be used is 25M (25.2 mm in diameter) Section at joint B From Table 2 Mf = -291.9 kN.m. Assuming flexural lever arm of 0.75 h = 0.75*600 = 525 mm
As = 291.9*106 = 1635.3 mm 2 0.85* 400*525
Slab reinforcement within (3*hslab) from the side of the beam is effective 4-10 M bars are effective (As = 400 mm2) As required = 1635.3 400 = 1235.3 mm2 Try 4- 20 M (As = 1200 mm2) Note: it would be unwise to be too conservative when designing top reinforcement. Since beam shear, joint shear, column moments and column shears will be increased.

+ Clause 21.3.2.2: M R ve 0.5* M R ve
M R min 0.25* M R max
at column face at any section
Assume As+ve = 4-15M As = 1200 +400 = 1600 mm2 d = 700 30 11.3 19.2/2 = 649 mm As = 4 * 200 = 800 mm2 d = 30 11.3 16/2 = 49.3 mm fc = 30 MPa 1 = 0.85 0.0015 * fc = 0.805 1 = 0.97 0.0025 * fc = 0.895
0.72 h 0.72*700 504 dv = = = 0.9 d 0.9*649 584.1
4 - 20 M 110
10 M @ 300 mm top&bottom
700
4- 15 M 0.0035 400
Ts = As* s*fy = 1600 * 0.85 * 400 = 544*103 N Cs = As* ( s*fs 1 * c*fc ) = 800 * (0.85 * fs 0.805 * 0.65 * 30) N C s = 1 * 1* c*fc * b * C = 0.805 * 0.895 * 0.65 *30 * 400 * C
-Ts + Cs + C s = 0 544*10 680 * (fs 18.47) 5619.7 * C = 0

x
=0
. (1)
d' s' = cc 1 = 0.0035 Compatibility of strains c 49.3 fs= E s * s' = 700 1 c
49.3 1 c
. (2)
Solving (1) and (2) C = 69.74 mm ,and fs = 205.1 MPa a = 0.895*69.74 = 62.42 mm
a M R ve = Cc d + Cs' ( d d ' ) = 333.3 kN.m. > M f = 291.9 kN.m 2
Clause 21.3.2.1
1.4* bw * d 1.4* 400*649 = = 908.6 < As o.k fy 400 As max = 0.025* bw * d = 0.025* 400*649 = 6400 > As o.k As min =
Bottom Bars for +ve reinforcement Beffective of the T-Section = 2*7.4/10+0.4 = 1.88 m For Section a Try As = 8-15 M d 635 mm As * s * fy 1600*0.85* 400 a= = = 18.43 mm < tflange=110 mm ' c * 1 * f c * b 0.65*0.805*30*1880
a 18.34 M R ve = Ts *(d ) = 1600*0.85* 400*(635 ) = 340.42 kN.m. > Mf = 276.1 kN.m. 2 2
O.K. MR+ve at column face d 650 mm a = 8.72 mm

MR
+ve
As= 4-15 M =800 mm2 = 175.6 > 0.5 * 333.3 = 166.6 kN.m
MR for section with As-ve = 2- 20 M
(needed later for bar cut-off)
for simplicity ignore bottom reinforcement a = 33.77 mm

MR = 129.5 kN.m
2- 20 M
400
Design of transverse reinforcement In order to avoid brittle shear failure of beams; shear deformation should always be in the elastic range. Shear design is based on the probable moment of resistance of the beam, so we assure that when plastic hinges are formed at the ends of the beam, the beam is still in the elastic range regarding shear deformation. Mp-ve = 1.47 * 333.33 = 490.00 kN.m Mp+ve = 1.47 * 175.60 = 258.00 kN.m
A
490.0 kN.m
137.5 kN
137.5 kN
B
258.0 kN.m
2.2 m
3m
2.2 m
238.6 kN
36.4 kN
P = 1.0 D + 0.5 L = 111 + 0.5* 53 = 137.5 kN

238.6 kN 101.1 kN
36.4 kN 490.0 kN.m
34.9 kN.m 258.0 kN.m 338.2 kN.m
Region of Plastic Hinging L = 3.0 m
Note: when reversing the direction of the lateral load B.M.D. & S.F.D. will be mirrored i.e. MB = 490 kN.m and VB = 238.6 kN

Section at column face
Vf = 238.6 kN Clause 21.3.4.2 =0 (excessive cracking is expected at high lateral drift due to E.Q.) o = 45
Vs =
s * Av * fy * dv *cot( )
s Try 4-legged 10 M, Av=400 mm2

S= 0.85* 400* 400*584.1*1 = 332.9 mm 238.6*103
Check shear design requirements (i) Maximum shear (clause 11.3.3)

Vr , max = 0.25* c * f c' * bw * dv = 0.25*0.65*30*400*584.1=1139.0 kN < Vf = 238.6 O.K.
(ii)
Minimum amount of stirrups (11.2.8.2)

S Av * fy 0.06* f c' * bw 400* 400 1217 mm 0.06* 30 * 400 O.K.
(iii)
Spacing limit (11.3.8.3)

0.125* c * f c' * bw * dv = 0.125*0.65*30* 400*584.1 = 569.5 kN > Vf 600 Smax = 0.7*584.1 = 408.87 mm
( No reduction in Smax is needed )

> S O.K.
Check Anti-buckling requirements
(Clause 21.3.3.2)
Hoops are required for a distance of (2d) from the face of the column. The maximum hoop spacing (S) is (i) (ii) (iii) (iv) d/4 = 649/4 = 162.25 8*d bar, longitudinal = 8*19.5 = 156 24* d bar, hoops = 24*11.3 = 271 300
governs
Use 4-legged 10 M hoops @150 mm for a distance of (2d) = 1.25 m Section at mid span Vf = 101.1 kN Use 2-legged U-shape stirrups
S= 0.85* 200* 400*584.1*1 = 392.86 mm 101.1*103
Check plastic hinging region (Clause 21.3.3.1) Hoops have to be provided for the regions where plastic hinging may occur and for a distance (d) beyond it from previoiusly drawn B.M.D. LPlastic Hinging 3.0 m hoops has to be provided for a distance of 3.0 + 0.649 = 3.65 m (from both column faces as the B.M.D. is mirrored case of reverse lateral loading) 7.4 2*3.65 = 0.1 m (insignificant length. Provide hoops for the whole span) Reinforcement cut-off Table 3.1 (Page 3-15 Part II of the CSA A23.3) For M 20 Ld= 530 mm (for bottom bars k1=1.0) = 1.3 * 530 = 689 mm (for top bars k1=1.3) For M 15 Ld= 390 mm (for bottom bars k1=1.0) = 1.3 * 390 = 507 mm (for top bars k1=1.3)
The theoretical cut-off location of the top 2-M 20 is @ 1.51 m from the column face. Clause 12.10.4 states that an embedment length of at least (d) or (12 dbar) is provided
beyond the theoretical cut-off point.

Lr + dv cot( ) = 1.51 + 0.584*cot(35) = 2.34 m L= d ld + max 12* dbar = 0.69 + 0.584 = 1.28 m
Bar Splice (Clause 12.15.1) For the 2-20 M bars LSplice = 1.3 * Ld = 1.3* 689 900 mm
1.25 D + 1.50 L (Cut-off of bottom renforcement) 276.1 235.8 Mr Mf 340.7 1.71 theoretical cut-off Point ( M f = 175.6 ) 175.6
175.6
Formation of Plastic Hinges at bends of the beam (Cut-off of top renforcement)
1.51 1.51
490.0 M f = 129.5 258.0 338.2 338.2 258.0 M r = 129.5
490.0
C L
2- Legged 10 M hoops @ 150 mm 2
0.9 2-20M L 2 =3.65 m 2-20M L 1=3.1 m 2-10M L=3.0 m 4-15M L 4 =5.7 m 4-15M L=8.15 m 12 d b = 240 mm L dh = 300 mm (taken as 500 - column dimension) L 1 = 2.34+0.50+0.24 3.1 m L 2 = 3.45+0.50+0.24+0.45 4.65 m L 3 = 2*( 3.45+0.45)+0.50 8.30 m L 4 = 2*( 2.00+0.584*cot(35)) 5.7 m
2-20M L 3 =8.3 m 2-20M L=5.5 m
2-20M L=8.3 m 2-20M L=5.5 m 2-10M L=3.0 m 4-15M L 4 =5.7 m 4-15M L=7.9 m
4 - 20 M 110
4- 15 M 400
Section 1-1
2- 20 M + 2- 10 M 110
700
8- 15 M 400
Section 2-2
Design of the first story interior column From F.E. analysis, the straining actions at top & bottom of the column in consideration are depicted in Table 3. Factored axial load and B.M. are shown in Table 4.
Table 3: Straining actions at the column ends PE PL () MD ML Location PD Units kN kN kN kN.m kN.m Bottom -1329.68 -540.89 2.68 2.62 1.37 Top -1329.68 -540.89 2.68 6.35 3.27
ME () kN.m 213.81 112.84
VD kN 2.56 2.56
VL kN 1.33 1.33
VE () kN.m 93.33 93.33
Table 4: Factored Axial loads and Moments 1 2 3 1.25D+1.5L 1.0D+1.0E 1.0D-1.0E

Pf,Bottom [kN] Pf,top [kN] Mf,bottom [kN.m] Mf,top [kN.m]
4 1.0D+0.5L+1.0E
5 1.0D+0.5L-1.0E
-2473.44
-2473.44 5.33
-1327.00
-1327.00
-1332.36
-1332.36
-1597.45
-1597.45
-1602.81
-1602.81
216.43
119.20
-211.19
-106.49
217.11
120.83
-210.51
-104.85
12.85
Preliminary selection of column Try 500 X 500 mm

500
8 - 25 M
As = 8 25 M bars = 4000 mm2 As,min = 0.01*500*500 = 2500 < As O.K. As,max = 0.06*500*500 = 15000 > As O.K. (Clause 21.4.3.1)
500
Check column capacity Draw the section interaction diagram and check if all load combinations is within the limits
5000 4500 4000 3500 3000
P [kN]
1
2500 2000
5
1500 1000 500 0 0 50 100 150 200
4 2
250
300
350
400
450
500
B.M [kN.m]
Strong column-weak beam requirements (clause 21.4.2.2) To ensure the formation of plastic hinges in the beams CSA A23.3 requires that: Mnc MPb For case 2 Mr = 411.0 kN.m Mn= 1.2* Mr = 493.2 kN.m For case 3 Mr = 415.0 kN.m Mn= 1.2* Mr = 498.0 kN.m For case 4 Mr = 425.5kN.m Mn= 1.2* Mr = 510.6 kN.m For case 5 Mr = 429.0kN.m Mn= 1.2* Mr = 514.8 kN.m Case 1 does not involve lateral load MPb = 490 + 258 = 748.0 kN.m Mnc = 2*493.2 = 986.4 kN.m Mnc > MPb O.K.
Design of transverse reinforcement of the column (Clause 21.4.5.1) Vmax = 2.56 + 0.5*1.33 + 93.33 = 96.55 kN Shear in columns when plastic hinges are formed in beams
+ M coltop = ( M Pr + M Pr ) *
kc 1 = (490 + 258) * = 374 kN 2 kc
The first story column is connected to a strong foundation, hence it is expected that the column will hinge at the bottom
M PColumn ,max = 1.57 * 432 = 678.24 kN .m
VColumn =
374 + 678 = 277 kN 3.8
Vdesign = 277 kN Shear Design (clause 21.4.5.2)

Vc = 0.65*0.1* 30*500*(0.72*500) = 64.08 kN Vs = 277 - 64.08 = 212.9 kN
Using transverse reinforcement as depicted in figure Av = (2+2*cos(45)) *100 = 341 mm2

0.85*341* 400*0.72*500*cot(45) = 196.0 mm 212.9*103
S=
0.125* c * f c' * bw * dv = 0.125*0.65*30*500*0.72*500 = 438.75 kN > Vf 600 mm Smax = 0.7*0.72*500 = 252 mm > S O.K.
Av * fy 0.06* f * bw
' c
341* 400 830 mm 0.06* 30 *500
O.K.
Confinement requirements (Clause 21.4.4.2) nl 8 kn = = = 1.33 nl 2 6
Po = 1 * f c' ( Ag Ast ) + f y * Ast = 0.805*30*(500*500-4000)+400*4000=7541 kN Pf 1602 = = 0.212 kp = Po 7541
Hence the total area of reinforcement is :

Ash = 0.2* kn * kp * Ag f c' * * S * hc Ach fyh
5002 30 = 0.2*1.33*0.212* * * 420* S 4202 400 = 2.5175 * S f c' 30 = 0.09* 420* * S = 2.835* S Ash , min = 0.09* S * hc * 400 fyh 341 for Ash = 341 mm 2 S= = 120.3 mm 2 2.835
Anti-buckling requirements Maximum spacing is : (i) (ii) (iii) b/4 = 500/4 = 125 mm 6*d bar * 25 = 150 mm
Sx = 100 +
hx
350 hx 350 190 = 100 + = 153 mm 3 3
Use 10M hoops @ 120 mm As this is the first story column, hoops shall be provided for the whole clear height (clause 21.4.4.6).
Splice details Columns bars should be spliced at mid-height with a tension lap splice Ld=820 mm (Table 3.1 A23.3 page 3-15) Lap length = 1.3*820 = 1066 mm
8 - 25 M
500
500
Column Reinforcement Details
1066
Check Joint Shear
VCol.
F2 f2
F1 f1
Vjoint = F1+F2-Vcol. F 1 = 1.25* As * fy = 1.25*800* 400 = 400 kN F 2 = 1.25*1200* 400 = 600 kN

490
374 258 374 374
Vcol =
374* 2 = 197 kN 3.8
Vjoint = 600 + 400 197 = 803 kN
Vr , Joint = 2.2 * c * f c' * Aj = 2.2*1.0*0.65* 30 *500* 400 = 1566.5 kN > VJoint OK

374 374 490 374 258
Shear Walls
Ductile Shear Walls Rd = 3.5, R0 = 1.6 Moderately Ductile Shear Walls Rd = 2.0, R0 = 1.4 Conventional Constructional Shear Walls Rd = 1.5, R0 = 1.3
Check of Ductility of Shear Walls o Walls effectively continuous in cross section o Plastic hinge at the base
(21.6.7)
id =
f Rd R0 f w
(hw lw / 2)
0.004
ic = (
cu l w
2c
0.002) 0.025
id = inelastic rotational demand on a wall ic = inelastic rotational capacity of a wall cu = 0.0035 lw = length of wall hw = vertical height of wall w = wall overstrength factor =
Ratio of the load corresponds to nominal resistance factored load
1.3
-ve +ve C 0.0035
Ductile Flexural walls
(21.6.1.1)
For hw/lw > 2.0 hw For hw/lw 2.0
Rd = 3.5 Rd = 2.0
Plastic hinge Region = 1.5 lw
lw
( 21.6.3.1)
Effective flange width
> distance to adjacent wall

25% of the wall height above the Section under consideration
For all elevations above the plastic hinge region the design SF and BM shall be increased by a ratio of the Factored Moment resistance to the factored moment (21.6.2.2.C)
BM* factored resistance/factored moment
Original design BM
Factored resistance
Factored moment
Rectangular wall (Plastic hinge Region)
Concentrated RFT
Spacing of distributed RFT 300 mm (21.6.5.2)
bw
lw
Strain diagram
C
0.0035
If C> min. of {4 bw => bw > lu/10 {0.3lw lu : Clear distance between floors 0.0025 distributed 0.0600 (21.6.4.3)
(21.6.3.4)
0.0015 bwlw Asconcentrated (min. 4 bars) 0.06 x Area of concentrated RFT (21.6.4.3) C 0.55lw If C> 0.14 wlw => confine the compression region as a column w = 1/s = 1.18 (over strength factor)
db
1/10 bw 25 mm 20 mm
(21.6.4.4)
25 cm wall => db 20 cm wall => db
*Tie concentrated RFT as columns: Spacing of Ties least of { 6db {24 dst { bw 0.0025 (21.6.6.9)
*Horizontal RFT ratio *Splices:-
(21.6.5.1)
Concentrated RFT: * max. 50% at the same location *1.5 ld * half the height of the each floor must be clear of splices Distributed RFT: * 1.5ld
Rectangular wall (Outside Plastic hinge Region)
Spacing of distributed RFT 450 mm (21.6.5.1)
bw lw
0.0025 distributed 0.0600
(21.6.4.3)
0.0010 bwlw Asconcentrated (min. 4 bars) 0.06 x Area of (21.6.4.4) concentrated RFT
(21.6.4.3)
*db
1/10 bw
(21.6.4.4)
*Ties for concentrated RFT are the same as those for columns.
*Horizontal RFT ratio * Splices: 1.5ld
0.0025 (21.6.4.1)
Shear Strength Requirements for Ductile Walls Vr > V f Vr the smallest of { Vp -> shear corresponding to the development of the probable nominal capacity of the wall { V determined from the analysis with RdR0 = 1
For id 0.015 => shear demand =0
0.10cf'cbwdv
For id = 0.015 => shear demand = 0.18 For Ps = 0.1f'cAs For Ps = 0.2f'cAs
0.15cf'cbwdv
= 450 = 350
CHAPTER 8 LESSONS LEARNED FROM THE PERFORMANCE OF BUILDING DURING PAST EARTHQUAKES
Post earthquake reconnaissance reports have provided overviews of damage patterns and detailed descriptions of specific failure modes and structural weaknesses. In addition, many reports have attempted to explain the behaviour of individual buildings by correlating analytical and experimental predictions with observed damage. The major lessons learned from past earthquakes have repeated themselves over and over and have confirmed the principles applied in modern seismic codes.
Lessons learned 1) Designing to Code does not always provide safeguard against excessive damage in severe earthquakes. 2) Well-designed, well detailed and well constructed buildings usually resist earthquake loads without excessive damage. 3) Poor construction practice can lead to severe damage and collapse 4) Ground failure can cause severe damage and collapse. 5) Buildings subjected to successive earthquakes may suffer progressive weakening. 6) Ductility and redundancy provide safety against collapse. 7) Stiff elements, not considered in the design, affect the seismic response of a building.
Primary examples are moment resisting frames that are filled with stiff masonry bricks as shown in next page.
Typical damage Inflicted to Infill Walls and Concrete frame Elements (From EERI, Pub. No. 86-02)
Framed bays with infill walls become stiffer once the frame is in contact with infill wall. This increases the stiffness of the frame and consequently the dynamic properties of the structural system.
If the infill walls are distributed unsymmetrically in plan, torsional modes of vibration can be excited. The infill walls will crack if they are not designed for the interaction forces occurring between frame and wall. The interaction forces may also damage the frame.
8) Problems with soft stories: A soft story can be generated when a shear system is replaced at the first story level with a number of columns. Earthquakes have shown that the stiffness discontinuity with height often directs forces to places where Structures strength perform is minimal. better in Soft story
earthquakes when they have uniform, or gradually changing, stiffness and strength over the height. Abrupt changes in stiffness or strength are responsible for some of the dramatic earthquake failure.
Soft Story Effect (From EERI, Pub. No. 86-02)
9) Problems with vertical geometric irregularities
10) Importance of Horizontal Diaphragms: Horizontal diaphragms (floor system) are responsible of carrying seismic forces to the vertical-load resisting building elements. Failure of a horizontal diaphragm leads to concentration of forces on certain elements of the lateral resisting system.
11) Problems with horizontal (plan) irregularities: a) Torsion irregularity b) Reentrant corners c) Diaphragms discontinuity d) Non-perpendicular systems
Torsion irregularity opening Reentrant corner
Non-perpendicular systems
12) Problems with corner columns 13) Problems with exterior panels and parapets 14) Unreinforced masonry buildings usually perform very poorly 15) Precast concrete elements must be well tied together. 16) Steel buildings generally perform well(except the extensive number of failures in welded beam-column connections observed during the 1994 Northridge earthquake) 17) Damage caused by fires and damage caused to non-structural components and building contents often exceed the consequences of inadequate structural performance.

Lecture Notes

Uploaded by

Document Informationclick to expand document information

Copyright:

Available Formats

Lecture Notes

Uploaded by

Document Information

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

Lecture Notes

Uploaded by

Copyright:

Available Formats

LECTURE NOTES

List of Proposed References

CHAPTER 1 EARTHQUAKE GROUND MOTIONS CHARACTERISTICS

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

X1, X2, and X3 are recording stations

Fig. 1 Determination of Epicenter

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Fig. 2 Time History Earthquake Records

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

No damage 50m 200 m Fig. 4 1000 m

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

CHAPTER 2 RESPONSE OF A SINGLE DEGREE OF FREEDOM SYSTEM

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

2-1 Equation of motion of a SDOF

& FD = Dashpot constant x velocity = cx FS = Spring stiffness x displacement = kx

m&& + cx + kx = P(t) x &

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Elastic Stiffness for Elements: k = P/ Force Required for Unit displacement

1) Springs in Series: K1 P1 P1 = P2 = P P1 = K11 ;

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

3) Axial Stiffness of a Bar: P L 4) Torsion: T L 5) Beams in Flexure: a) Cantilever: P

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

K1 and K2 are two springs in parallel (have the same displacement)

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty

Types of Dynamic Loading: a) Free Vibration: X(t)

b) Forced Harmonic Vibration: P(t)

c) Impulse Loading: (Driving of Piles)

Time d) Random Vibration: P(t)

Seismic Analysis and Design of Buildings

By: Prof. A. A. El Damatty