Seismic Design of Structures Overview
Seismic Design of Structures Overview
PG 2019
Spring 2020 Semester
Fawad A. Najam
Department of Structural Engineering
NUST Institute of Civil Engineering (NICE)
National University of Sciences and Technology (NUST)
H-12 Islamabad, Pakistan
Cell: 92-334-5192533, Email: fawad@[Link]
Acknowledgement
• The material for the preparation of these lectures slides are taken from different sources.
• The primary source for these lecture slides are the lectures of Dr. Naveed
Anwar and Dr. Pennung Warnitchai at Asian Institute of Technology (AIT),
Thailand
• Class Notes of Prof. Dr. Worsak Kanok-Nukulchai at Asian Institute of Technology (AIT), Thailand
• The material is taken solely for educational purposes. All sources are duly acknowledged.
3
Full 3D, Nonlinear, Equations, Charts,
Partial Differential Closed Form with 2D/3D Linear Static
Inelastic Dynamic Tables, Rules,
Equations Approximations FEA/Matrix
FEA Limits
• Fairly General, Unified (FEM, BEM ...) • Structural Material (RC, PSC, HRS, CFS, timber ... )
• Input: Structural Geometry, Element/Member • Design Code (ACI, BS Codes, EuroCode, JIS ... )
Cross-sections, Reinforcement details
• Design Approach (working stress, ultimate strength, limit
• Output: Element/ Member Actions, state ... )
Displacements ... etc.
• Structural Members (beams, columns, slabs, footings ... )
Structural
Model Linear Nonlinear
E, A, I, L, G etc. = Constant, K= Constant E ≠ Constant, EI ≠ Constant, K≠ Constant
Seismic
Loading
Equivalent static analysis method is commonly used for the seismic design of ordinary and "regular"
buildings and structures.
Dynamic analysis methods (Response spectrum and Time history) are required for "irregular"
buildings and structures, very important structures, and structures that seismic response is not
dominated by the fundamental vibration mode. ... more tedious, more difficult.
• Definition
• Natural vibration of a structure released from initial condition and subjected to no external load
or damping
• Main governing equation -Eigenvalue Problem
𝑴 𝑼ሷ 𝑡 + 𝑪 𝑼ሶ 𝑡 + 𝑲 𝑼 𝑡 = 𝑷(𝑡)
• Solution gives
• Natural Frequencies
• Associated mode shapes
• An insight into the dynamic behavior and response of the structure
• A mode shape is a set of relative (not absolute) nodal displacements for a particular mode
of free vibration for a specific natural frequency
• The Time Period and Mode Shapes are inherent properties of the structure and do not
depend on the applied loads
• The Modal Analysis should be run before applying loads any other analysis to check the
model and to understand the response of the structure.
• Modal analysis is precursor to most types of analysis including Gravity Load Analysis,
Response Spectrum Analysis, Time History Analysis, Push-over Analysis, etc.
• Modal analysis is a useful tool even if full Dynamic Analysis is not performed.
• The Time Period and Mode Shapes, together with animation immediately exhibit the strengths
and weaknesses of the structure.
• Modal analysis can be used to check the accuracy of the structural model
– The Time Period should be within reasonable range,
– The disconnected members are identified
– Local modes are identified that may need suppression
• T2 = 4.96 sec
Translation in
• Structure
Mass • Attachments
• Occupants
• Tower
• Podium
Stiffness • Basement
• Foundation
• Soil
Mode 1
Mode 2
Mode 3
≅
+ + +…
A Detailed 3D Elastic
Structural Model
• Mass source
• DL = Dead load
• LL MEP = MEP and other permanent equipment live load
• LL STO = Storage live load
𝑚 𝑢ሷ 𝑡 + 𝑐 𝑢ሶ 𝑡 + 𝑘 𝑢 𝑡 = −𝑚 𝑢ሷ 𝑔 (𝑡)
𝑘 2
𝑐 =2𝑚𝜉𝜔 , 𝜔 =
𝑚
Solution
𝑢𝑜 𝑇, 𝜉 = 𝑚𝑎𝑥 𝑢(𝑡, 𝑇, 𝜉)
𝑢ሶ 𝑜 𝑇, 𝜉 = 𝑚𝑎𝑥 𝑢(𝑡,
ሶ 𝑇, 𝜉)
𝑢ሷ 𝑜 𝑇, 𝜉 = 𝑚𝑎𝑥 𝑢(𝑡,
ሷ 𝑇, 𝜉)
Equation of motion: 𝑚 𝑢ሷ + 𝑐 𝑢ሶ + 𝑘 𝑢 = −𝑚 𝑢ሷ 𝑔 𝑡
It turns out that, for a given ground acceleration 𝑢ሷ 𝑔 𝑡 , the deformation response depends
only on 𝜔𝑛 (or 𝑇𝑛 ) and 𝜉 of the SDOF system.
(b) (c)
Deformation response of three Deformation response spectrum
SDOF systems with 𝜉 = 0.02 for 𝜉 = 0.02
Once 𝑢0 is obtained from the deformation response spectrum, the corresponding peak
internal forces 𝑓𝑠0 can be determined by:
𝑓𝑠0 = 𝑘 𝑢0
𝐴 is not the real peak acceleration response but it has units of acceleration.
Dynamic Static
Stationary
at the instant where 𝑢 𝑡 = 𝑢0
𝐴
It can be written in the form 𝑉𝑏0 = 𝓌 (33)
𝑔
*It is used in the building codes to represent the coefficient by which the structural weight is
multiplied to obtain the base shear.
1 2 1 2 1 1
𝐸𝑠0 = 𝑘𝑢 0 = 𝑘𝐷 = 𝑚𝜔 𝑛 𝐷 = 𝑚𝒱 2
2 2 (34)
2 2 2 2
Solution 𝑓𝑆𝑜
The unit of mass is therefore the unit of force divided by the unit of acceleration.
𝐹 Free-end
The lateral stiffness 𝑲 in this case is determined from
𝛿
𝐹
𝐾 =
𝛿
Fixed-End
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 41
The lateral stiffness of this SDF system is
The total weight of the pipe is 10.79 × 12 = 129.5 𝑙𝑏, which may be
neglected relative to the lumped weight of 5200 𝑙𝑏 . Thus
𝑤 5.20
𝑚= = = 0.01347 𝑘𝑖𝑝 − 𝑠𝑒𝑐 2 Τ𝑖𝑛.
𝑔 386
The natural vibration frequency and period of the system are
𝑘 0.211
𝜔𝑛 = = = 3.958 𝑟𝑎𝑑Τ𝑠𝑒𝑐 𝑇𝑛 = 1.59 𝑠𝑒𝑐
𝑚 0.01374
𝐷 = 5.0 𝑖𝑛. = 𝑢𝑜
𝐴 = 0.20g
𝑀𝑐 (12.48 × 12)(4.5Τ2)
𝜎𝑚𝑎𝑥 = = = 46.5 𝑘𝑠𝑖
𝐼 7.23
Deformation response of a
SDOF system with very
long natural period
A very long period system is extremely flexible. The mass would be expected to remain
essentially stationary while the ground below moves.
• For Park Filed 66 earthquake where ground motion is very short and shock like, the effect
of damping would be small, as in the case of half cycle sine pulse excitation.
• The construction of response spectra plots requires the solution of single degree of
freedom systems for a sequence of natural frequency and of the damping ratio in the
range of interest.
• Every solution provides only one point (the maximum value) of the response spectrum.
• Since a large number of systems must be analyzed in order to fully plot each response
spectrum, the task is lengthy and time consuming even with the use of computer.
• Once these curves are constructed and are available for the excitation of interests, the
analysis for the design of structures subjected to dynamic loading is reduced to a
simple calculation of natural frequency of the system and the use of response spectra.
• Dynamic analysis of a system with 𝒏 degree of freedom can be transformed to the problem of
solving n systems in which each one is a single degree of freedom system
• The understanding and mastery of the concepts and methods of solutions for a single degree
of freedom system is quit important.
• Consider a quantity 𝑉 for an SDF system with natural frequency 𝜔 related to its peak
deformation 𝐷 ≡ 𝑢𝑜 due to earthquake ground motion:
2𝜋
𝑉 = 𝜔𝐷 = 𝐷
𝑇
• Using the above expression, the displacement response spectrum can be converted to
the Pseudo-velocity response spectrum.
• Consider a quantity 𝐴 for an SDF system with natural frequency 𝜔 related to its peak
deformation 𝐷 ≡ 𝑢𝑜 due to earthquake ground motion:
2
2𝜋 2
𝐴=𝜔 𝐷= 𝐷
𝑇
• Using the above expression, the displacement response spectrum can be converted to
the Pseudo-acceleration response spectrum.
• The quantity 𝐴 is related to the peak value of base shear 𝑉𝑏𝑜 or the peak value of the
equivalent static force 𝑓𝑠𝑜 .
𝑉𝑏𝑜 = 𝑓𝑠𝑜 = 𝑚𝐴
𝐴 𝐴
𝑉𝑏𝑜 = 𝑤 , 𝑓𝑠𝑜 = 𝑤
𝑔 𝑔
• The three spectra (deformation, pseudo-velocity, and pseudo-acceleration) are simply different
ways of presenting the same information on structural response for a given ground motion.
• Knowing one of the spectra, the other two can be obtained by algebraic operations mentioned
earlier.
The response spectrum for a given ground motion component 𝑢ሷ 𝑔 (𝑡) can be developed by
implementation of the following steps:
1) Numerically define the ground acceleration 𝑢ሷ 𝑔 (𝑡); typically, the ground motion
ordinates are defined every 0.02 sec.
2) Select the natural vibration period 𝑇 and damping ratio 𝜉 of an SDF system.
3) Compute the deformation response 𝑢(𝑡) of this SDF system due to the ground
motion 𝑢ሷ 𝑔 (𝑡) by any of the numerical methods.
4) Determine 𝑢𝑜 , the peak value of 𝑢(𝑡).
Arturo Danusso (1880 - 1968) Theodore von Karman (1881 - 1963) Maurice Anthony Biot (1905 - 1985)
Mathematical Formulation:
Summary:
Notes:
Notes:
• Static methods specified in building codes are based on single mode response and appropriate
for simple and regular structures
• Dynamic analysis should be used for complex buildings to determine significant response
characteristics
• Effects of structure’s dynamic characteristics on vertical distribution of lateral forces
• Increase in dynamic loads due to torsional motions
• Influence of higher modes, resulting in an increase in story shear and deformations
The experimental results from cyclic loading conditions indicate that the cyclic force-deformation
behavior of a structure depends on
o the structural material (concrete, steel)
o the type of structural members (beam, shear member, axial member)
o how members are assembled into the structural system.
Force 𝑭
Deformation 𝒖
Hysteresis Loop
Deformation
Force
Deformation
Deformation
𝜇 = 𝑑𝑢𝑐𝑡𝑖𝑙𝑖𝑡𝑦 𝑓𝑎𝑐𝑡𝑜𝑟 = 𝑢𝑚 /𝑢𝑦
In order to understand the effects of the inelastic force-deformation relation on the earthquake
response, it is necessary to evaluate the peak deformation of an elastroplastic system and
compare this deformation to the peak deformation of the corresponding linear system.
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 89
RESPONSE OF ELASTOPLATIC SYSTEM TO EARTHQUAKE GROUND MOTION
𝑑2 𝑢 𝑑𝑢 𝑑 2 𝑢𝑔
Equation of motion: 𝑚 2 + 𝑐 + 𝑓𝑠 𝑢 = −𝑚
𝑑𝑡 𝑑𝑡 𝑑𝑡 2
The nonlinear resisting force 𝑓𝑠 can also be obtained from the elastoplastic force-deformation
relation.
For a low-level ground motion, the system may behave like a linear elastic system with
1 𝑘
𝑓= 𝑇 = 1ൗ𝑓 ,
2𝜋 𝑚
where 𝑘 = initial stiffness in the elastic range.
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 90
EFFECTS OF INELASTICITY ON EARTHQUAKE RESPONSES
Response of a linearly elastic system with 𝑇=0.5 sec and 𝜉 = 0 to El Centro ground motion:
Deformation
𝑢𝑡
𝑓
Force 𝑠ൗ𝑊
Time intervals
of yielding
𝑢𝑡
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 92
EFFECTS OF INELASTICITY ON EARTHQUAKE RESPONSES
• For each system, the computer ductility factor 𝜇 is the “ductility demand” imposed on
elastoplastic system by the ground motion.
• The system should be designed such that its “ductility capacity” (i.e. the ability to
deform beyond the elastic limit) exceeds the “ductility demand”.
Ductility
demand For systems with long natural periods,
𝜇 The ductility demand 𝜇 ≈ 𝑓0 /𝑓𝑦
Note:
𝐷𝑦 = 𝑢𝑦
the yield deformation
𝑉𝑦 = 𝜔𝑛 𝑢𝑦
the yield pseudo-velocity
𝐴𝑦 = 𝜔𝑛2 𝑢𝑦
the yield pseudo-acceleration
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 101
DESIGN OPTIONS
𝑓𝑦
If the structure is designed for strength 𝑓0 = ൗ
𝑊
0.919 𝑊 or larger, it will remain within its linearly
elastic range during the El Centro ground motion;
therefore, it needs not to be ductile.
𝑓𝑠 3% post-yielding
stiffness
(Miranda, ASCE [Link]. Eng., 119, No5, May, 1993) 𝑢
INELASTIC DESIGN RESPONSE SPECTRA
The behavior of the structure under a monotonically increasing load would be:
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 109
UNIFORM BUILDING CODE (1994)
𝑍𝐼𝐶
Base Shear Design Base Shear 𝑉𝑏 = 𝑊
𝑅𝑤
𝑅𝑤 = Reduction Factor = 𝑅Ω 𝑅𝑦
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 110
REDUCTION FACTOR
𝑅𝑦
𝑅Ω 𝑅𝛺 = factor for reducing “yield strength demand” (𝑓𝑦 )
to “design base shear” (𝑉𝑏 ).
𝑅𝛺 is associated with the overstrength of the system.
Natural period, 𝑇
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 111
REDUCTION FACTOR 𝑹𝒘
𝑯
𝑹𝒘
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 112
REDUCTION FACTOR 𝑹𝒘 TABLE 16N-STRUCTURAL SYSTEMS
𝑹𝒘 𝑯
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 𝑁𝑜𝑡𝑒: 𝑵. 𝑳. = 𝑁𝑜 𝑙𝑖𝑚𝑖𝑡 113
REDUCTION FACTOR 𝑹𝒘 TABLE 16N-STRUCTURAL SYSTEMS
𝑹𝒘 𝑯
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 114
The Code-based Response Spectrum Analysis (RSA) Procedure
Performance-based Seismic Design of Buildings – Semester: Spring 2018 (Fawad A. Najam) 115
The Code-based Response Spectrum Analysis (RSA) Procedure
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 116
The Standard RSA Procedure (ASCE 7-10, IBS 2012, EC 8)
Linear Elastic Model Determine Spectral Acceleration Determine Elastic Base Shears Reduce Elastic Base Shear to account for
for each Significant Mode inelasticity
Base Shear
Spectral Acceleration (SA)
N
𝑉𝑒𝑙
Stories 𝑆𝐴3
𝑅 = Response Modification
For Initial
Viscous Damping
… Factor, ASCE 7 (or
Behavior Factor, EC 8)
𝑆𝐴2 𝑉𝑖𝑛
𝜙1 𝜙2 𝜙3 𝑉𝑒𝑙
𝑉𝑖𝑛 =
𝑉𝑏1 𝑉𝑏2 𝑉𝑏3 𝑅/𝐼
𝑆𝐴1
Eigen-value Analysis 𝑁
𝕂 − 𝜔2𝕄 Φ = 𝕆 ∆𝑒𝑙 = ∆𝑖𝑛
𝑉𝑏𝑛 = 𝛤𝑛 . 𝑚𝑖 . 𝜙𝑖,𝑛 . 𝑆𝐴𝑛
Determine Modal Properties 𝑇3 𝑇2 𝑇1 𝑖=1
𝑇𝑖 , 𝜙𝑖 , 𝛤𝑖 Time Period (sec) 𝑉𝑒𝑙 = (𝑉𝑏1 )2 + (𝑉𝑏2 )2 + (𝑉𝑏3 )2 + … 𝑥𝑟 Roof Displacement
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 118
Modal Combination Rules
• ABSSUM Rule N
• Add the absolute maximum value from each mode. Not so ro rn 0
popular and not used in practice n =1
• SRSS
N
• Square Root of Sum of Squares of the peak response from
each mode. Suitable for well separated natural frequencies.
ro n0
r 2
n =1
• CQC
N N
• Complete Quadric Combination is applicable to large range of
structural response and gives better results than SRSS.
ro
i =1 n =1
r r
in i 0 n 0
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 119
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 120
The Input – Output Summary
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 121
The RSA Procedure in UBC 97
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 122
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 123
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 124
Dynamic Analysis Procedure in IBC 2000
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 125
The RSA Procedure in IBC 2000
Design Spectral Values
• Adjust Maximum Considered Earthquake
(MCE) values of 𝑆𝑠 and 𝑆1 for local site effects
𝑆𝑀𝑆 = 𝐹𝑎 × 𝑆𝑠
𝑆𝑀1 = 𝐹𝑣 × 𝑆1
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 126
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 127
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 128
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 129
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 130
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 131
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 132
Routine Design Office Practice
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 133
Principal Directions
• Lack of definitions of the principal directions in code
• The direction of the base reaction of the mode shape associated with the fundamental frequency of the
system is used to define the principal direction.
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 134
Directional and Orthogonal Effects
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 135
Accidental Torsion
• Arise from several factors • Structures assigned to Seismic Design Category C,
• Rotational components of ground motions D, E, or F, where Type 1a or 1b torsional irregularity
• Effects of non-structural elements exists, accidental torsional moment needs to be
• Actual distribution of dead and live loads amplified.
• Uncertainties in defining building’s material
properties for dynamic analysis
• Generally 5% of eccentricity from center of mass is
considered
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 136
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 137
Seismic Design as per
ASCE 7-10
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 138
The RSA Procedure (ASCE 7-10)
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 139
The RSA Procedure
(ASCE 7-10)
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 140
The RSA Procedure (ASCE 7-10)
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 141
The RSA Procedure (ASCE 7-10)
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 142
The RSA Procedure (ASCE 7-10)
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 143
Strength Design Load Combinations (ASCE 7-16)
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 144
Response Spectrum Analysis Procedure in BCP (2007)
Performance-based Seismic Design of Buildings – Semester: Spring 2018 (Fawad A. Najam) 145
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 146
Design Response Spectrum (UBC 1997, p 2-38)
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 148
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 149
The Response Spectrum Analysis (RSA) Procedure
Performance-based Seismic Design of Buildings – Semester: Spring 2018 (Fawad A. Najam) 150
The RSA Procedure as prescribed
in ASCE 7-16
Software Demonstration
ETABS v 2016
Performance-based Seismic Design of Buildings – Semester: Spring 2020 (Fawad A. Najam) 151
Thank you