E Kalkan

Practical Guidelines to Select and
Scale Earthquake Records for

Nonlinear Response History Analysis of
Structures
By Erol Kalkan and Anil K. Chopra
Open-File Report 2009–####

U.S. Department of the Interior
Practical Guidelines to Select and Scale Earthquake
Records for Nonlinear Analysis of Structures
U.S. Geological Survey

Practical Guidelines to Select and
Scale Earthquake Records for
Nonlinear Response History Analysis of
Structures
By Erol Kalkan and Anil K. Chopra
Open-File Report 2009–####
Suggested citation:
Kalkan E. and Chopra A.K., 2009, Practical Guidelines to Select and Scale Earthquake
Records
for Nonlinear Response History Analysis of Structures: U.S. Geological Survey Open-File
Report 2009–####, 113 p.
U.S. Department of the Interior

Records for Nonlinear Response History Analysis of Structures
U.S. Geological Survey
ii
Abstract
Earthquake engineering practice is increasingly using nonlinear response history analysis
(RHA) to demonstrate performance of structures. This rigorous method of analysis requires
selection and scaling of ground motions appropriate to design hazard levels. Presented herein is
a modal-pushover-based scaling (MPS) method to scale ground motions for use in nonlinear
RHA of buildings and bridges. In the MPS method, the ground motions are scaled to match (to
a specified tolerance) a target value of the inelastic deformation of the first-mode inelastic SDF
system whose properties are determined by first-mode pushover analysis. Appropriate for first-
mode dominated structures, this approach is extended for structures with significant
contributions of higher modes by considering elastic deformation of higher-mode SDF systems
in selecting a subset of the scaled ground motions. Based on results presented for two bridges
and six actual buildings, covering low-, mid-, and high-rise building types in California, the
accuracy and efficiency of the MPS procedure are established and its superiority over the
ASCE 7-05 scaling procedure is demonstrated.
1
2
Acknowledgements
Dr. Kalkan wishes to acknowledge the generous support of the Earthquake Engineering Research
Institute in particular Prof. Thalia Anagnos, Mr. Marshall Lew and Prof. Jonathan Bray for
providing him the 2008 EERI/FEMA NEHRP Professional Fellowship in Earthquake Hazard
Reduction. Special thanks are extended to Prof. Rakesh Goel, Prof. Sashi K. Kunnath, Prof.
Carlos Ventura and Dr. Emrah Erduran for their support in development of the OpenSEES
models.
3
4
Contents
Page
ABSTRACT........................................................................................................................1
ACKNOWLEDGEMENTS.................................................................................................3
1. INTRODUCTION.........................................................................................................7
1.1. RESEARCH OBJECTIVE.............................................................................................9
2. MODAL-PUSHOVER-BASED SCALING................................................................10
2.1. MPS PROCEDURE: SUMMARY................................................................................10
2.2. ESTIMATING SDF-SYSTEM INELASTIC DEFORMATION........................................13
3. CODE-BASED SCALING PROCEDURE................................................................15

4. GROUND MOTIONS SELECTION PROCEDURE..................................................16
4.1 GROUND MOTION ENSEMBLE.................................................................................17
5. EVALUATION OF MPS PROCEDURE: LOW- AND MID-RISE BUILDINGS.........24

5.1 SELECTED BUILDINGS.............................................................................................24
5.2. SYSTEMS ANALYZED...............................................................................................26
5.3 FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................27
5.4. EVALUATION OF MPS CONCEPT.............................................................................27
5.5. EVALUATION OF MPS AND CODE-BASED SCALING PROCEDURES...................28
5.6. MULTI-MODE CONSIDERATIONS.............................................................................30
6. EVALUATION OF MPS PROCEDURE: HIGH-RISE BUILDINGS..........................52

6.1 BUILDINGS SELECTED.............................................................................................52
6.2 ANALYTICAL MODELS AND CALIBRATION STUDIES............................................54
6.3. FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................55
6.4. EVALUATION OF MPS PROCEDURE.......................................................................56
6.5. COMPARING MPS AND CODE-BASED SCALING PROCEDURES.........................57
6.6. MULTI-MODE CONSIDERATIONS.............................................................................58
7. EVALUATION OF MPS PROCEDURE: SHORT - PERIOD BUILDING.................80

7.1 BUILDING DETAILS AND ANALYTICAL MODEL......................................................80
5
7.2. FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................81
8. EVALUATION OF MPS PROCEDURE: BRIDGES.................................................90

8.1 DESCRIPTION OF BRIDGES AND ANALYTICAL MODELS.....................................90
8.2 FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................93
9. CONCLUSIONS.....................................................................................................105
10. APPENDIX-A: COMPUTATION OF MINIMUM SCALE FACTOR.....................107
11. REFERENCES....................................................................................................108
NOTATION.....................................................................................................................112
6
1.INTRODUCTION
Seismic evaluation of existing structures and of proposed design of new structures is usually
based on nonlinear static (or pushover) analysis procedures, but nonlinear response history
analysis (RHA) is now being increasingly used. In the latter approach, the seismic demands are
determined by nonlinear RHA of the structure for several ground motions. Procedures for
selecting and scaling ground motion records for a site-specific hazard are described in building
codes [e.g. IBC 2006 (ICBO 2006) and CBC 2007 (ICBO 2007)] and have been the subject of
much research in recent years.
Current performance-based design and evaluation methodologies prefer intensity-based
methods to scale ground motions over spectral matching techniques that modify the frequency
content or phasing of the record to match its response spectrum to the target spectrum. In
contrast, intensity-based scaling methods preserve the original non-stationary content and only
modify its amplitude. The primary objective of intensity-based scaling methods is to provide
scale factors for a small number of ground motion records so that nonlinear RHA of the structure
for these scaled records is accurate, i.e., it provides an accurate estimate in the median value of
the engineering demand parameters (EDPs), and is efficient, i.e., it minimizes the record-to-
record variations in the EDP. Scaling ground motions to match a target value of peak ground
acceleration (PGA) is the earliest approach to the problem, which produces inaccurate estimates
with large dispersion in EDP values [Nau and Hall 1984; Miranda 1993; Vidic et al. 1994;
Shome and Cornell 1998]. Other scalar intensity measures (IMs) such as: effective peak
acceleration, Arias intensity and effective peak velocity have also been found to be inaccurate
and inefficient [Kurama and Farrow 2003]. None of the preceding IMs consider any property of
the structure to be analyzed.
Including a vibration property of the structure led to improved methods to scale ground
motions, e.g., scaling records to a target value of the elastic spectral acceleration, A ( T1 ) from the
code-based design spectrum or PSHA-based uniform hazard spectrum at the fundamental
vibration period of the structure, T1, provides improved results for structures whose response is
dominated by their first-mode [Shome et al. 1998]. However, this scaling method becomes less
7
accurate and less efficient for structures responding significantly in their higher vibration modes
or far into the inelastic range [Mehanny 1999; Alavi and Krawinkler 2000; Kurama and Farrow
2003]. To consider higher mode response, a scalar IM that combines the spectral accelerations
A ( T1 ) and A ( T2 ) at the first two periods and vector IM comprised of A ( T1 ) and the ratio of
A ( T1 ) / A ( T2 ) have been developed [Bazzurro 1998; Shome and Cornell 1999]. Although this
vector IM improves accuracy, it remains inefficient for near-fault records with a dominant
velocity pulse [Baker and Cornell 2006].
To recognize the lengthening of the apparent period of vibration due to yielding of the
structure, a scalar IM defined as a combination of A ( T1 ) and A ( cT1 ) where c > 1 , has been
considered [Mehanny 1999; Cordova et al. 2000]; alternatively, scaling earthquake records to
minimize the difference between its elastic response spectrum and the target spectrum has been
proposed [Kennedy et al. 1984; Malhotra 2003; Alavi and Krawinkler 2004; Naeim et al. 2004;
Youngs et al. 2007].
International Building Code (IBC) [ICBO 2006] and California Building Code (CBC)
[ICBO 2007] require that earthquake records be scaled according to the ASCE 7-05 provisions
[ASCE 2005]. For two dimensional analyses of regular structures, ground motions are scaled
such that the average value of the 5%-damped elastic response spectra for a set of scaled motions
is not less than the design response spectrum over the period range from 0.2T1 to 1.5T1. For
structures having plan irregularities or structures without independent orthogonal lateral load
resisting systems where three-dimensional analyses need to be performed, ground motions
should consist of appropriate horizontal components.
All the preceding scaling methods utilize IMs based on elastic response of the structure,
but do not explicitly consider its inelastic response. They lead to scale factors that depend only
on the structural period(s), independent of the structural strength. The elastic-response-based IMs
may not be appropriate for near-fault sites where the inelastic spectral deformation can be
significantly larger than corresponding elastic spectral deformation [Bozorgnia and Mahin 1998;
Alavi and Krawinkler 2000; Baez and Miranda 2000; Chopra and Chintanapakdee 2004]. This
limitation has been overcome in recently proposed IMs based on the inelastic deformation
spectrum, leading to improved estimate of the median EDPs and reduced dispersion of EDPs
[Bazzurro and Luco 2004; Luco and Cornell 2007]. Through incremental dynamic analyses,
response of generic frames to different intensity levels of near-fault ground motions
8
demonstrated that scaling records with the IM defined as the inelastic deformation of the first-
mode inelastic SDF system is accurate, efficient and sufficient compared to elastic-response-
based IMs [Tothong and Luco 2007; Tothong and Cornell 2008]. Required in this approach are
attenuation relationships for the inelastic deformation with given ground motion properties
(magnitude, fault distance, site condition etc.) and mean rate of occurrence of the hazard level
considered [Tothong and Cornell 2008].
1.1. RESEARCH OBJECTIVE
The objective of this report is to develop a new practical method for selecting and scaling
earthquake ground motion records in a form convenient for evaluating existing structures or
proposed designs for new structures. The selection procedure presented considers the important
characteristics of the ground motions (e.g., pulse, directivity, fling, basin, duration etc.) consistent
with the hazard conditions. The scaling procedure presented explicitly considers structural strength
and is based on the standard IM of spectral acceleration that is available from the USGS seismic
hazard maps, where it is mapped for periods of 0.2 sec and 1.0 sec for the entire U.S. to facilitate
construction of site-specific design spectrum [Petersen et al. 2008] or it can be computed from
the uniform hazard spectrum obtained by probabilistic seismic hazard analysis (PSHA) for the
site.
Based on modal pushover analysis, the procedure presented herein explicitly considers
the strength of the structure, obtained from the first-mode pushover curve and determines scaling
factors for each record to match a target value of the deformation of the first-mode inelastic SDF
system estimated by established procedures. Appropriate for first-mode dominated structures,
this approach is extended for structures with significant contributions of higher modes. Based on
results presented for computer models of six actual buildings [4-story reinforced concrete (RC),
4-, 6-, 13-, 19- and 52-story steel special moment resisting frame (SMRF)] and two bridges [two
span and multi-span], the effectiveness of this scaling procedure is established and its superiority
over the ASCE 7-05 scaling procedure is demonstrated.
9
2. MODAL-PUSHOVER-BASED SCALING
In the modal pushover-based scaling (MPS) procedure, each ground motion record is scaled by a
scale factor selected to ensure that the peak deformation of the first-mode inelastic SDF system
due to the scaled record is close enough to a target value of the inelastic deformation. The force-
deformation relation for the first-mode inelastic SDF system is determined from the first-mode
pushover curve. The target value of the inelastic deformation is the median deformation of the
inelastic SDF system for a large ensemble of (unscaled) earthquake records compatible with the
site-specific seismic hazard conditions. Nonlinear RHA of the inelastic SDF system provides the
peak deformation of the system to each record in the ensemble, and the median of the data set
provides the target value. Alternatively, the median deformation of the inelastic SDF system can
be estimated as the deformation of the corresponding linearly elastic system, known directly
from the target spectrum, multiplied by the inelastic deformation ratio; empirical equations for
this ratio are available for systems with known yield-strength reduction factor [Chopra and
Chintanapakdee 2004].
For first-mode dominated structures, scaling earthquake records to the same target value
of the inelastic deformation is expected be sufficient. Because higher vibration modes are known
to contribute significantly to the seismic response of mid-rise and high-rise buildings, the MPS
procedure checks for higher-mode compatibility of each record by comparing its scaled elastic
spectral displacement response values at higher-mode vibration periods of the structure against
the target spectrum. This approach ensures that each scaled earthquake record satisfies two
requirements: (1) the peak deformation of the first-mode inelastic SDF system is close enough to
the target value of the inelastic deformation; and (2) the peak deformation of the higher-mode
elastic SDF system is not far from the target spectrum.
2.1.MPS PROCEDURE: SUMMARY
The MPS procedure is summarized below in a step-by-step form:
1. For the given site, define the target pseudo-acceleration response spectrum either as the
PSHA-based uniform hazard spectrum, or code-based design spectrum, or the median
10
pseudo-acceleration spectrum for a large ensemble of (unscaled) earthquake records

compatible with the site-specific seismic hazard conditions.
2. Compute the frequencies n (periods Tn ) and f n of the first few modes of elastic
vibration of the structure.
First-mode Dominated Structures
3. Develop the base shear-roof displacement Vb1-ur1 relation or pushover curve by nonlinear
static analysis of the structure subjected to gradually increasing lateral forces with an
invariant force distribution s1 = mf1, associated with the first-mode, where m is the
*
structural mass matrix. Gravity loads, including those present on the interior (gravity)
frames, are applied before starting the pushover analysis.
4. Idealize the pushover curve and select a hysteretic model for cyclic deformations, both
appropriate for the structural system and materials [Han and Chopra 2005; Bobadilla and
Chopra 2007]. Determine the yield-strength reduction factor R y (equals strength required
for the structure to remain elastic divided by the yield strength of the structure) from:
R y = M1* A1 Vb1 y , where M1* is the effective modal mass and Vb1 y is the yield point
value of base shear determined from the idealized pushover curve.

5. Convert the idealized pushover curve to the force-deformation Fs1 L1 - D1 relation of the
first-mode inelastic SDF system by utilizing Fs1/L1= Vb1 / M 1 and D1 = ur1 / 1fr1 in which
*
L1 = f1mi , fr1 is the value of f 1 at the roof, 1 = (f1T mi )/(f1T mf1 ) and each element of
the influence vector i is equal to unity.

I
6. For the first-mode inelastic SDF system, establish the target value of deformation D1
from D1 = CR D1 , where D1 = ( T1 2p ) A1 and A1 is the target pseudo-spectral

I 2
acceleration at period T1 and C R is determined from an empirical equation (shown in the

next section) for the inelastic deformation ratio corresponding to the yield-strength
reduction factor R y , determined in Step 4.

7. Compute the peak deformation D1 = max D1 (t ) of the first-mode inelastic SDF system
I
defined by the force deformation relation developed in Steps 4 and 5 and damping ratio
11
1. The initial elastic vibration period of the system is T1 = 2p(L1D1y / Fs1y)1/2 . For a SDF
I
system with known T1 and 1, D1 can be computed by nonlinear RHA due to one of the
selected ground motions u&

&g (t ) multiplied by a scale factor SF, to be determined to satisfy
Step 8, by solving
&
D&+ 2  D& + F [ D ] / L = -(SF)u&
&
1 1 1 1 s1 1 1 g (t ) (2.1)
I
8. Compare the normalized difference between the target value of the deformation D1 of
I
the first-mode inelastic SDF system (Step 6) and the peak deformation D1 , determined in
Step 7 against a specified tolerance, e
D1 = D1I - D1I / D1I < e (2.2)
9. Determine the scale factor SF such that the scaled record (SF) u&
&g ( t ) satisfies the criterion
of Eq. (2). Because Eq. (1) is nonlinear, SF cannot be determined a priori, but requires an
iterative procedure starting with an initial guess. Starting with SF = 1, Steps 7 and 8 are
implemented and repeated with modified values of SF until Eq. (2) is satisfied.
Successive values of SF are chosen by trial and error or by a convergence algorithm, e.g.,
Newton-Raphson iteration procedure. For a given ground motion, if the Eq. 2 is satisfied
by more than one SF, the SF closest to one should be taken.
Repeat Steps 7 and 8 for as many records as deemed necessary; obviously the scaling
factor SF will be different for each record. These scaling factors will be shown to be
appropriate for structures that respond dominantly in the first-mode.
Higher-mode Considerations
10. Establish target values of deformation of higher-mode SDF systems, treated as elastic
systems, from the target spectrum Dn = ( Tn 2p ) An , where the mode number n = 2. We

2
have found that the second-mode is mostly adequate for buildings susceptible to higher-
mode effects.
11. By linear RHA, calculate the peak deformation D2 �max t D2 (t ) of the nth-mode elastic
SDF system with known T2 and 2 due to a selected ground motion u&
&g (t ) multiplied by its
scale factor SF determined in Step 9.

12. Compute the normalized difference between the target value of deformation D2 (Step 10)
and the peak deformation determined in Step 11.
12
D 2 = D2 - D2 / D2 (2.3)
and rank the scaled records based on their D 2 value; the record with the lowest average D 2
is ranked the highest.
13. From the ranked list, select the final set of records with their scale factors determined in
Step 9 to be used in nonlinear RHA of the structure.
2.2.ESTIMATING SDF-SYSTEM INELASTIC DEFORMATION
The inelastic deformation ratio C R is required in Step 6 to estimate the deformation of the
inelastic SDF system. Such equations were first developed by Veletsos and Newmark [1960] as a
function of elastic vibration Tn and ductility factor m . However, in selecting and scaling ground
motion records for nonlinear RHA of an existing building or of a proposed design of a new
building, the inelastic deformation ratio should be expressed as a function of Tn and the yield-
strength reduction factor R y ; these quantities are determined in Steps 7 and 4, respectively. The
inelastic deformation ratio can be expressed as a function of elastic vibration period and yield-
strength reduction factor Ry. Response data for 216 ground motions recorded on NEHRP site
classes B, C, and D demonstrated that the mean inelastic deformation ratio is influenced little by
soil condition, by magnitude if Ry < 4 (but significantly for larger Ry), or by site-to-fault distance
so long as it exceeds 10 km [Ruiz-Garcia and Miranda 2002]. Regression analysis of these data
led to an equation for the inelastic deformation ratio as a function of Tn and Ry; this equation is
restricted to elastoplastic systems.
Median values of C R have been presented for non-degrading bilinear hysteretic systems
subjected to seven ensembles of far-fault ground motions (each with 20 records), representing
large or small earthquake magnitude and distance, and NEHRP site classes, B, C, or D; and for
two ensembles of near-fault ground motions. Regression analysis of these data led to the
empirical equation [Chopra and Chintanapakdee 2004]:
-1
� �a �T1 ��
�
d
-1
CR = 1 + �
( LR - 1) + � b + c �
� �� (2.4)
�R �
�
� � y �Tc ��
� �
in which, the limiting value of CR at Tn = 0 is:
13
1 � Ry - 1 �
LR = 1+
� � (2.5)
Ry � a �
where a is the post-yield stiffness ratio and Tc is the period separating the acceleration and
velocity-sensitive regions of the target spectrum; the parameters in Eq. (2.4) are: a=61, b=2.4,
c=1.5, and d=2.4.
Equations (2.4) and (2.5) and values of their parameters are valid for far-fault ground
motions, independent of (1) earthquake magnitude and distance; and (2) NEHRP site class B, C,
and D; and also for near-fault ground motions.
3.CODE-BASED SCALING PROCEDURE

The procedures and criteria in the 2006 IBC and 2007 CBC for the selection and scaling of
ground motions for use in nonlinear RHA of structures are based on the ASCE 7-05 provisions
14
[ASCE 2005]. According to ASCE 7-05, earthquake records should be selected from events of
magnitudes, fault distance and source mechanisms that comply with the maximum considered
earthquake. If the required number of appropriate records is not available, appropriate simulated
ground motions may be included to make up the total number required.
For two-dimensional analysis of symmetric-plan buildings, ASCE 7-05 requires intensity-
based scaling of ground motion records using appropriate scale factors so that the average value
of the 5%-damped response spectra for the set of scaled records is not less than the design
response spectrum over the period range from 0.2T1 to 1.5T1. The design value of an engineering
demand parameter (EDP)—member forces, member deformations or story drifts—is taken as
the average value of the EDP over seven (or more) ground motions, or its maximum value over
all ground motions, if the system is analyzed for fewer than seven ground motions.
The ASCE 7-05 scaling procedure does not insure a unique scaling factor for each record;
obviously, various combinations of scaling factors can be defined to insure that the average
spectrum of scaled records remains above the design spectrum (or amplified spectrum incase of
3-D analyses) over the specified period range. Because it is desirable to scale each record by the
smallest possible factor, an algorithm is developed (Appendix-A) and used in applying the code-
scaling procedure in Section 4.
4.GROUND MOTIONS SELECTION

PROCEDURE
15
Both record selection and scaling are equally important processes for success of any nonlinear
RHA applied. Appositely selecting records considering the hazard conditions for a given site
helps to reduce the dispersion of EDPs and increase accuracy by achieving better estimates of the
“true median”. Before scaling ground motions, one needs to define the hazard conditions
associated with a given site either through deterministic or probabilistic site-specific hazard
analysis or alternatively from the USGS seismic hazard maps. The parameters that need to be
considered in identifying the scenario conditions are those that have the most influence on
ground motion spectral shape (Graizer and Kalkan 2009):
- Magnitude range of anticipated significant event(s)
- Distance range of the site from the causative fault(s)
- Site-condition (site-geology generally described by average shear-wave velocity

within 30 m)
- Basin effect (if basin exists)
- Directivity effect
Spectral shape (i.e., response spectrum normalized by PGA) defines ground motion demand
characteristic on MDF systems. Therefore in selecting candidate records for nonlinear RHAs,
one needs to carefully identify records whose spectral shapes are close to each other. The
dependence of ground motion spectral shape on the first three parameters coloring seismic
hazard conditions is explained in detail in the following:
Magnitude Dependence
In general, events with larger magnitude yield wider response spectra. In order to find the degree
of magnitude influence on response spectral shape, average spectral shapes of earthquakes
ranging from magnitude 4.9 to 7.9 in the extended NGA database (Graizer and Kalkan 2009) are
plotted in Figure 4.1. As shown, the spectral peak gradually shifts from ~0.15 sec for the lowest
magnitude earthquake (M4.9) to ~0.5 sec for the largest events (M7.6 - 7.9). Maximum
amplitudes of the average spectral shape are relatively stable varying from 2.3 to 2.6, with higher
amplitudes at smaller magnitudes.
16
Distance Dependence
As reported in previous studies (e.g., Abrahamson and Silva 1997), predominant period shifts to
higher values with increase in distance from the fault for a given earthquake. Figure 4.2 depicts
such a distance dependence on the spectral shape whereby variations of maximum period for
different distance bins are plotted for the 1999 M7.6 Chi-Chi earthquake. For this particular
event, predominant period shifts from about 0.35 sec at the closest distances (0-20 km bin) to 1.2
sec at farthest fault distances (120-140 km bin). Similar observations are valid for the other
major events investigated.
Soil Condition Dependence

In addition to magnitude and distance dependence, spectral shape also depends on site
conditions. Predominant period of spectral shape from a rock site is generally lower than that of
a soil site. This tendency is shown in Figure 4.3 which compares the average spectral shape in
VS30 bin of 180-360 m/sec with that of 540-720 m/sec. At short periods, spectral shape of average
rock site remains above that of average soft-soil, while the reverse takes place beyond the
spectral period of 0.3 sec. Unlike predominant period shifting to higher values with reduction in
VS30, peak intensity of spectral shapes remains similar.
4.1 GROUND MOTION ENSEMBLE
In order to test the MPS procedure, a total of twenty one near-fault strong earthquake ground
motions were compiled from the Next Generation Attenuation project earthquake ground motion
database [Power et al. 2006]. These motions were recorded during seismic events with moment
magnitude, M ≥ 6.5 at closest fault distances, Rcl ≤ 12 km and belonging to NEHRP site
classification C and D. The selected ground motion records and their characteristic parameters
are listed in Table 4.1.
17
Shown in Figure 4.4 are the pseudo-acceleration response spectrum for each ground
motion and the median of the 21 response spectra. The median spectrum is taken to be the design
spectrum for purposes of evaluating the MPS procedure.
The median spectrum of the ground motion ensemble is presented next in Figure 4.5 as a
four-way logarithmic plot, together with its idealized version (dashed-line). The idealized
spectrum is divided logically into three period ranges: the long-period region to the right of point
d, Tn > Td, is called the displacement-sensitive region; the short-period region to the left of point
c, Tn < Tc , is called the acceleration-sensitive region; and the intermediate-period region between
points c and d, Tc < Tn < Td, is called the velocity-sensitive region [Chopra 2007; Section 6.8].
Note that the velocity-sensitive region is unusually narrow, which is typical of near-fault ground
motions.
Table 4.1 Selected earthquake ground motions
18
19
Figure 4.1 Comparison of average spectral shape of earthquakes in magnitude range

of 4.9 to 7.9 (Increase in magnitude shifts the predominant period to higher
values).
Figure 4.2 Shift of predominant period of average spectral shape to higher values with
increase in average distance within each 20 km distance bin (Data values
correspond to the 1999 M7.6 Chi-Chi earthquake).
20
Figure 4.3 Comparison of average “rock” and “soil” spectral shapes and their transfer
function (Rock/Soil).
21
SPECTRAL REGIONS
Acceleration Velocity Displacement
sensitive sensitive sensitive
)
m
(c
p.
is
2
10 10
D
10
2
Pseudo Velocity (cm/s)
0
10
10
1
1
-1
10
10
10
0
-2
10
Ps
eu
10
-1
do
0
A
-3
10
10
cc
.(
10
-2
cm
/s
²)
Ta Tb Tc Td Te Tf
-1 0.04 0.15 1 3 8 33
10
-1 0 1
10 10 10
Period (s),Tn
Figure 4.4 Median elastic response spectrum for the selected ensemble of ground
motions shown by a solid line, together with its idealized version in dashed
line; spectral regions are identified;  = 5%.
22
2.5
Median
ith GM
2
Pseudo Acceleration (g)
1.5
1
Design Spectrum
0.5
0
0 1 2 3 4 5
Period (s), Tn
Figure 4.5 Individual response spectra for 21 ground motions and their median
response spectrum; = 5%.
23
5.EVALUATION OF MPS PROCEDURE:

LOW- and MID-RISE BUILDINGS
The efficiency and accuracy of the MPS and ASCE 7-05 scaling procedures will be first
evaluated based on the three steel SMRF buildings representing low- and mid-rise building types
in California.
A scaling procedure is considered efficient if the dispersion of EDPs due to the scaled
records are small; it is accurate if the median value of the EDPs due to scaled ground motions is
close to the benchmark results, defined as the median values of EDPs, determined by nonlinear
RHA of the building to each of the 21 unscaled ground motions (Chapter 4). In this section, the
median values of EDPs determined from a set of 7 ground motions, scaled according to MPS and
ASCE 7-05 scaling procedures, will be compared. The median value, x̂ defined as the geometric
mean and the dispersion measure,  of n observed values of xi are calculated from
1/ 2
�n � �n �
�� ln xi � �� (ln xi - ln xˆ ) 2 �
xˆ = exp �i =1 �;  = �i =1 � (5.1)
� n � � n -1 �
�
� �
� �
� �
�
The EDPs selected are peak values of story drift ratio, i.e., peak relative displacement
between two consecutive floors normalized by story height; floor displacements normalized by
building height; column and beam plastic rotations.
5.1 SELECTED BUILDINGS
The first set of buildings selected to evaluate the efficiency and accuracy of the MPS method are
existing 4-, 6-, and 13-story steel SMRF buildings. The six and thirteen-story buildings are
instrumented, and their motions have been recorded during past earthquakes. A brief description
of these buildings and their analytical models follow.
24
5.1.1 Four-story Building

This building located in Northridge, California was designed in compliance with 1988 Uniform
Building Code (UBC). The structural system is composed of perimeter SMRFs to resist lateral
loads and interior gravity frames. The floor plan and elevation of the building together with beam
and column sizes are shown in Figure 5.1. In this and subsequent figures, gravity frames are
excluded from the plan. The columns are embedded into grade beams and anchored to the top of
pile caps; thus, displacements and rotations are essentially restrained in all directions. All
columns are made of A-572 grade 50 steel. The girders and beams are made of A-36 steel. The
floor systems are composed of 15.9 cm thick slab (8.3 cm light weight concrete and 7.6 cm
composite metal deck). The total seismic weight of the building was estimated to be
approximately 10,880 kN.
5.1.2. Six-story Building

The six-story building located in Burbank, California was designed in 1976 in accordance with
the 1973 UBC requirements (Figure 5.2). It has a square plan measuring 36.6 m by 36.6 m with a
8.3 cm thick light weight concrete slab over 7.6 cm metal decking. Shear studs between the slab
and beams were provided on the interior beams in the North-South direction only. The structural
system is essentially symmetrical in plan. Moment continuity of each of the perimeter frames is
interrupted at the ends where a simple shear connection is used to connect to the transverse-
frame column along its weak axis. The plan view of the building and the elevation of a typical
frame together with member sizes are shown in Figure 5.3. The interior frames of the building
were designed as gravity frames and consist of simple shear connections only. All columns are
supported by base plates anchored on foundation beams, which in turn are supported on a pair of
9.75 - 0.75 m diameter concrete piles. Section properties were computed for A-36 steel with an
assumed yield stress of 303 MPa as established from coupon tests [Anderson and Bertero 1991].
The minimum concrete compressive strength at 28 days was 20.7 MPa, except it was 13.8 MPa
for the slab on grade. The total seismic weight of the building was estimated to be approximately
34,650 kN. The building was instrumented with a total of 13 strong motion sensors at the
ground, 2nd, 3rd and roof levels.
25
5.1.3. Thirteen-story Building

Located in South San Fernando Valley, 5 km southwest of the epicenter of the 1994 Northridge
earthquake, this 13-story building (with one basement) was designed according to the 1973 UBC
(Figure 5.4). The footprint of the building is 53.3 by 53.3 m. The exterior frames of the building
are moment-resisting frames and interior frames were designed for gravity load bearing. The
foundation consists of piles, pile caps and grade beams. The floor plan of the perimeter frames
and a typical elevation of one of these frames are shown in Figure 5.5; note that member sizes
are indicated and that the corner columns are composed of box sections. Floor systems consist of
6.4 cm of 20.7 MPa-concrete fill over 7.5 cm steel decking. The roof system is lighter with 5.7
cm of vermiculite fill on 7.6 cm steel decking. The total seismic weight of the building was
estimated to be 68,950 kN. Motions of this building during the 1994 Northridge earthquake were
recorded by 7 sensors: three each in the North-South and East-West directions, and one in the
vertical direction.
5.2.SYSTEMS ANALYZED
These symmetric-plan buildings may be analyzed as two-dimensional systems for ground

motions applied independently along each axis. The frames selected for modeling are: SMRF
along Line 1 for the 4-story building, SMRF long Line 1 for the 6-story building, and the SMRF
along Line G for the 13-story building. The analytical models were generated in the open source
finite element platform, OpenSEES using a force-based nonlinear beam-column element
[Neuenhofer and Filippou 1998]. This element utilizes a layered ‘fiber’ section at each
integration point, which in turn is associated with uni-axial material models and enforces
Bernoulli beam assumptions for axial force and bending. Centerline dimensions were used in
element modeling. Included in the frame model was one half of the building mass distributed
proportionally to the floor nodes. Special features such as local connection fracture were not
simulated; consequently, modeling of members and connections was based on the assumption of
stable hysteresis loops derived from a bilinear stress-strain model. The columns were assumed to
be fixed at the base level. Complete details of the analytical models and calibration studies for
the six- and thirteen-story buildings against the recorded data are reported in Kalkan (2006). The
26
first three natural vibration periods and modes of each building are shown in Figure 5.6 and the
first-mode pushover curves in Figure 5.7, where P-D effects are included.
5.3 FIRST-MODE SDF-SYSTEM PARAMETERS
In order to establish first-mode SDF system properties for each building, the global hysteretic
behavior of buildings is described by the cyclic pushover curves [Han et al., 2004; Chopra and
Bobadilla, 2007] as shown in Figure 5.8. These were determined by nonlinear static analysis of
the buildings subjected to the modal force distribution ( s1 = mφ1 ) with its magnitude varied and
reversed to cause the cyclic roof displacement described by Figure 5.9. Comparison of the cyclic
and monotonic pushover curve in Figure 5.10 indicates cycle-to-cycle deterioration of stiffness
due to P-Δ effects.
The force-deformation relation for the first-mode SDF system is determined from the
base shear – roof displacement relation defined by modal pushover curve by utilizing Fs1/L1= Vb1
/ M 1 and D1 = ur1 / 1fr1. The hysteretic force-deformation relation is idealized by the peak-
*
oriented model [Ibarra and Krawinkler, 2005; Ibarra et al., 2005], with the monotonic curve
idealized as tri-linear (Fig. 5.10). In the hysteretic model, the cyclic behavior is described by a
series of deterioration rules. These deterioration parameters are determined from the cyclic force
deformation relation FS1 / L1 − D1 obtained from the cyclic pushover curve (Fig. 5.8) by an
iterative trial-and-error procedure to obtain a best fit of the model to the actual cyclic curve. The
associated hysteretic model provides reasonable representation of the global cyclic behavior of
the buildings as shown in Figure 5.11.
5.4.EVALUATION OF MPS CONCEPT
As a first step in evaluating the concept underlying the MPS procedure, the target value of
I
deformation D1 is computed not as described in Step 6 of the procedure (Chapter 2) but as the
median value of peak deformation of the first-mode inelastic SDF system due to 21 ground
I
motions determined by nonlinear RHA. The MPS method utilizing this D1 value is denoted
27
henceforth as MPS*. The 21 ground motions are divided into 3 sets each containing 7 records
(Table 5.1). The records in each set are selected randomly from at least 3 different earthquakes to
avoid any dominant influence of a single event on the ground motion set. An appropriate scale
factor for each record is determined by implementing Steps 1-8 of the MPS procedure.
Efficiency and accuracy of the MPS* procedure are evaluated for each ground motion set
separately by comparing the median values of EDPs determined by nonlinear RHA of the
building due to the 7 scaled records against the benchmark EDPs. Figure 5.12 shows the
benchmark EDPs for all three buildings; results from individual records are also included to
demonstrate the large dispersion. Almost all of the excitations drive all three buildings well into
the inelastic range as shown in Figure 5.13 where the roof displacement values due to 21 ground
motions are identified on the first-mode pushover curve. Also shown is the median value. The
post-yield branch of the first-mode pushover curve exhibits negative slope because of P-Δ
effects.
Comparisons of the EPDs due to the MPS procedure with the benchmark results are
presented in Figures 5.14 through 5.16 for the three buildings. Included are the EDPs due to each
of the 7 scaled ground motions to show the dispersion of the data. The results are organized for
each building in three parts corresponding to the three ground motion sets. These results
demonstrate that the MPS procedure is accurate; the median values of EDPs due to every small
(7) subset of scaled ground motion closely match the benchmark results, which were determined
from a large (21) set of ground motions. The dispersion of the EDP values due to the 7 scaled
records about their median value is much smaller compared to the data for the 21 unscaled
records in Figure 5.17. These results collectively demonstrate that the concept underlying the
MPS procedure is accurate and efficient in scaling records for nonlinear RHA of buildings.
5.5.EVALUATION OF MPS AND CODE-BASED SCALING

PROCEDURES
The preceding implementation of the MPS concept is the same as the MPS procedure described
I I
in Section 2 except for how D1 was computed. Previously, the exact value of D1 was
determined by nonlinear RHA of the first-mode inelastic SDF system, but it will now be
28
estimated according to Step 6, using an empirical equation for C R , in accordance with the MPS
procedure. In utilizing CR equation, zero post-yield stiffness is assumed, although the idealized
first-mode SDF systems (Fig. 5.10) have negative post-yield stiffness. This choice is dictated by
the fact that the original CR equation was determined for stable systems with non-negative post-
yield stiffness ratio [Chopra and Chintanapakdee 2004].
In CR equation, using zero post-yield stiffness seems to be plausible because the
variability in the peak displacement demand is not affected significantly by the hysteretic
behavior [Kurama and Farrow, 2003; Gupta and Kunnath, 1998]. Figure 5.17 compares the
I
“exact” target value of deformation, D1 (continuous horizontal line) with estimated target value
I
of deformation, D1 (dashed horizontal line) using the CR equation with zero post-yield stiffness;
D1I values from individual records for each of the three buildings are also included. Notably,
I I
“exact” target value of deformation, D1 is determined as the median value of D1 from 21 ground
motions based on first-mode inelastic SDF system having peak-oriented hysteretic behavior as
I
shown in Figure 5.11. As Figure 5.17 indicates, “exact” and estimated values of D1 are close to
each other; the discrepancy between them becomes less as the initial period, Tn of the inelastic
SDF system prolongs.
I
Once D1 is estimated, an appropriate scale factor for each record is determined based on
the inelastic first-mode SDF systems (Fig. 5.11) by implementing Steps 7-8 of the MPS
procedure. Table 5.2 lists the scaled factors computed for each set and for each building. The
EDPs determined by nonlinear RHAs of the structure due to a set of 7 ground motions scaled
according to MPS procedure are compared against the benchmark EDPs. Figures 5.18 - 5.26
present such comparisons for the three buildings and for the three sets of ground motions. These
results demonstrate that the MPS procedure is much superior compared to the ASCE 7-05
procedure for scaling ground motion records. This superiority is apparent in two respects: First,
for each building and each ground motion set, the ground motions scaled according to the MPS
procedure lead to median values of EDPs that are much closer to the benchmark values than the
corresponding results based on the ASCE 7-05 procedure. The only exception is the combination
of 6-story building and Ground Motion Set 3. In this case, the MPS procedure leads to estimates
of the EDPs that are only slightly better than results from the ASCE 7-05 scaling procedure.
29
Second, the dispersion in the EDP values due to the 7 scaled records around the median value is
much smaller when the records are scaled according to the MPS procedure compared to the
ASCE 7-05 scaling procedure. However, even with the MPS scaling, the dispersion of EDPs for
the upper stories of 6- and 13-story buildings is noticeable, particularly for Ground Motion Set 2,
indicating that the higher-mode contributions to the seismic demands are significant. These
factors are considered for 6- and 13-story buildings in the next chapter.
An alternative way of comparing MPS and ASCE 7-05 scaling methods is based on the
ratio of the EDP value due to a scaled record and the benchmark value. The deviation of the
median, D of this ratio from unity is an indication of the error or bias in estimating the median
EDP value, and the dispersion,  of this ratio is an indication of the scatter in the individual
EDPs, determined from the scaled ground motions. Included also in the comparison is the MPS*
I
procedure based on “exact” values of D1 instead of Step 6. Figure 5.27 presents the median, D
of the EDP ratio for story drifts determined from records scaled according to the MPS*, MPS,
and ASCE 7-05 scaling methods. Comparing these D values against 1.0, it is apparent that the
MPS* method is most accurate (least biased), the MPS method is only slightly less accurate. The
bias in the MPS methods is generally less than 20%, except in the case of the 6-story building
and Ground Motion Set 3. The ASCE 7-05 method is least accurate generally overestimates the
EDPs, with the overestimation exceeding 50% in some cases.
Figure 5.28 presents the dispersion of the EDP ratio for story drifts determined from
records scaled according to the MPS*, MPS, and ASCE 7-05 scaling methods. It is apparent that
the MPS* scaling method leads to the smallest dispersion, and it becomes only slightly larger in
the MPS method. Dispersion is largest in the ASCE 7-05 scaling method, becoming unacceptably
large for some combinations of buildings and ground motion sets.
5.6.MULTI-MODE CONSIDERATIONS
As demonstrate here, the MPS method based solely on the first-mode inelastic SDF system
(Steps 1-9 of the method) is superior over the ASCE 7-05 scaling method. Considering the higher
modes of vibration is expected to improve the method further for mid-rise and high-rise
buildings.
30
To consider the 2nd mode effect into account, the 21 records scaled based on the first-
mode response only (Steps 1-9 of the method) are ranked by accounting for elastic second mode
response according to Steps 10-13 of the method. The seven records with the highest ranks (see
Step 12) were defined as Ground Motion Set 4. Note that this set is different for each building
(Table 5.3).
Figure 5.29 compares the median EDPs from ground motions scaled by the MPS method
with the benchmark values for the 6- and 13-story building; results of individual ground motions
are also shown. It is apparent that considering higher modes in selecting ground motions in the
MPS method provides accurate estimates of the median EDPs and reduces the record-to-record
variability (compared to the results achieved by Ground Motion Sets 1-3). This improved
accuracy and efficiency is demonstrated in Figure 5.30, where the D and —the median value of
the ratio of the estimated story drift to its benchmark value, and dispersion of this ratio—are
plotted for the four set of ground motions. It is apparent that Ground Motion Set 4 is more
accurate and efficient than Ground Motion Sets 1 though 3.
The improvement achieved for these buildings is modest because the higher-mode
computations are not especially significant in the response of selected buildings. Such
improvement is expected to be more pronounced in the case of taller buildings responding
significantly in their higher modes.
Table 5.1 List of records for three ground motion sets
31
Table 5.2 Scale factors computed for three buildings and for three sets of seven
ground motions
Table 5.3 Scale factors for 6-, and 13-story buildings considering higher-modes
32
(a) (b)
Figure 5.1 4-story building: (a) plan; and (b) elevation of perimeter SMRF.
33
Figure 5.2 Overview of the 6-story steel building in Burbank, CA.
Moment resisting connection (b)

(a)
34
Figure 5.4 Overview of the 13-story steel building in South San Fernando Valley, CA.
Moment resisting connection
(a)
(b)
35
Figure 5.6 Natural vibration periods and modes of 4-, 6-, and 13-story buildings.
4-Story 6-Story 13-Story

0.3 0.3 0.15
Base Shear / Weight
0.2 0.2
0.09
0.1 0.1
0.03
0 0 0
0 1 2 3 0 1 2 0 1 2 3 4
Roof Disp. / Bldg. Height (%)
Figure 5.7 First-mode pushover curves for 4-, 6-, and 13-story buildings.
36

0.4 0.3 0.15
Base Shear / Weight
0 0 0
-0.4 -0.3 -0.15

-4 -2 0 2 4 -2 -1 0 1 2 -4 -2 0 2 4
Figure 5.8 First-mode cyclic pushover curve (solid line) and monotonic pushover
curve (dashed line), for 4-, 6-, and 13-story buildings.
4
Displacement D  D y
-2
-4
-6
0 1 2 3 4 5 6 7 8 9 10 11
Cycle Number
Figure 5.9 Description of roof-displacement amplitudes for cyclic pushover analysis.
37

0.3 0.12
0.3
0.2 0.08
FS1 / L1 (g)
0.2
0.1 0.04
0.1
0 0 0
0 10 20 30 40 0 10 20 30 40 50 0 50 100 150
D1 (cm)
Figure 5.10 Comparison of first-mode pushover curve (solid line) and its idealized
trilinear model (dashed line) for 4-, 6, and 13-story buildings.

0.5 0.4 0.15
FS1 / L1 (g)
0 0 0
-0.5 -0.4 -0.15

-50 0 50 -40 0 40 -150 0 150
D1 (cm)
Figure 5.11 Comparison of first-mode cyclic pushover curve (solid line) and its
hysteretic model (dashed lines), for 4-, 6-, and 13-story buildings.
38
4
4-Story
3
Floor
1
Median
ith GM
0
6
6-Story
5
4
Floor
2
1
0
14
12 13-Story
10
8
Floor
6
4
2
0
0 1 2 0 2 4 6 80 0.02 0.04 0 0.01 0.02 0.03
Floor Disp. / Bldg. Height (%) Story Drift Ratio (%) Col. Pl. Rot. (rad) Beam Pl. Rot. (rad)
Figure 5.12 Median values of EDPs determined by nonlinear RHA of three buildings for
21 ground motions; results for individual ground motions are also
included.

0.5 0.3 0.15
Base Shear / Weight
0.4
0.2 0.10
0.3
Median
Median
Median
0.2
0.1 0.05
0.1
0 0 0
0 1 2 3 0 1 2 0 1 2
Figure 5.13 Roof displacements determined by nonlinear RHA of three buildings for 21
ground motions identified on first-mode pushover curves.
39
4
GM Set-1
3
Floor
2
o Benchmark
1 MPS
Median
ith GM
0
4
GM Set-2
3
Floor
0
4
GM Set-3
3
Floor
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
Fl. Disp. / Bldg. Height (%) Story Drift Ratio (%) Col. Pl. Rot. (rad) Beam Pl. Rot. (rad)
Figure 5.14 Comparison of median EDPs based on the MPS concept with benchmark
EDPs for the 4-story building; individual results for each of the seven
scaled ground motions are also presented.
40
6
5 GM Set-1
4
Floor
3
o Benchmark
2
MPS
1 Median
ith GM
0
6
5 GM Set-2
4
Floor
3
2
0
6
5 GM Set-3
4
Floor
2
1
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
41
14
12 GM Set-1
10
8
Floor
6 o Benchmark
4 MPS
2 Median
ith GM
0
14
12 GM Set-2
10
8
Floor
6
4
2
0
14
12 GM Set-3
10
8
Floor
6
4
2
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
42

60 80 150
50
120
60
40
D1 (cm)
90
30
I
40
60
20
10 20 30
0 4 8 12 16 20 0 4 8 12 16 20 0 4 8 12 16 20
Ground Motion Number
I
Figure 5.17 Peak deformation D1 values of the first-mode inelastic SDF system for 21
ground motions for 4-, 6-, and 13-story buildings; “exact” target value of
I
deformation D1 is identified by horizontal continuous line; horizontal
I
dashed line indicates target value of deformation D1 established by CR
equation.
4
GM Set-1
3
Floor
2
o Benchmark
1 MPS
Median
ith GM
0
4
3
Floor
2
o Benchmark
1 ASCE-7
Median
ith GM
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
Figure 5.18 Comparison of median EDPs for Ground Motion Set 1 scaled according to
MPS (top row) and ASCE 7-05 (bottom row) scaling procedures with
benchmark EDPs; individual results for each of seven scaled ground
motions are also presented. Results are for the 4-story building.
43
4
GM Set-2
3
Floor
2
o Benchmark
1 MPS
Median
ith GM
0
4
3
Floor
2
o Benchmark
1 ASCE-7
Median
ith GM
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
4
GM Set-3
3
Floor
2
o Benchmark
1 MPS
Median
ith GM
0
4
3
Floor
2
o Benchmark
1 ASCE-7
Median
ith GM
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
Figure 5.20 Comparison of median EDPs for Ground Motion Set 3 computed for the 4-
story building based on the MPS and ASCE 7-05 earthquake record scaling
procedures with benchmark EDPs; individual results for each of seven
scaled ground motions are also presented. Results are for the 4-story
building.
44
6
GM Set-1
5
4
Floor
3
o Benchmark
2
MPS
1 Median
ith GM
0
6
4
Floor
3
o Benchmark
2
ASCE-7
1 Median
ith GM
0
0 1 2 0 1 2 3 4 0 0.02 0 0.01 0.02 0.03
6
GM Set-2
5
4
Floor
3
o Benchmark
2
MPS
1 Median
ith GM
0
6
4
Floor
3
o Benchmark
2
ASCE-7
1 Median
ith GM
0
0 1 2 0 1 2 3 4 0 0.02 0 0.01 0.02 0.03
45
6
GM Set-3
5
4
Floor
3
o Benchmark
2
MPS
1 Median
ith GM
0
6
4
Floor
3
o Benchmark
2
ASCE-7
1 Median
ith GM
0
0 1 2 0 1 2 3 4 0 0.02 0 0.01 0.02 0.03
Figure 5.23 Comparisons of median EDPs for Ground Motion Set 3 scaled according to
14
12 GM Set-1
10
8
Floor
6
o Benchmark
4 MPS
2 Median
ith GM
0
14
12
10
8
Floor
6
o Benchmark
4 ASCE-7
2 Median
ith GM
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
46
14
12 GM Set-2
10
8
Floor
6
o Benchmark
4 MPS
2 Median
ith GM
0
14
12
10
8
Floor
6
o Benchmark
4 ASCE-7
2 Median
ith GM
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
14
12 GM Set-3
10
8
Floor
6
o Benchmark
4 MPS
2 Median
ith GM
0
14
12
10
8
Floor
6
o Benchmark
4 ASCE-7
2 Median
ith GM
0
0 1 2 0 2 4 0 0.02 0.04 0 0.01 0.02 0.03
47
GM Set-1 GM Set-2 GM Set-3

4
MPS*
*
MPS
3
ASCE-7
Floor
4-story
0
6
4
Floor
1
6-story
0
14
12
10
8
Floor
2
13-story
0
0.5 1 0.5 1 0.5 1 1.5
, or
Figure 5.27 Median story drift ratios D MPS , D MPS * and D ASCE - 7 for three ground
motions sets and for three buildings.
48

4
MPS*
*
MPS
3
ASCE-7
Floor
4-story
0
6
4
Floor
1
6-story
0
14
12
10
8
Floor
2
13-story
0
0 0.2 0.4 0.6 0 0.2 0.4 0.6 0 0.2 0.4 0.6 0.8
, or
Figure 5.28 Dispersion of story drift ratios  MPS ,  MPS * and  ASCE - 7 for three ground
motions sets and for three buildings.
49
6
6-story
5 GM Set-4
4
Floor
3
o Benchmark
2
MPS
1 Median
ith GM
0
0 1 2 0 1 2 3 0 0.01 0.02 0 0.01 0.02 0.03
13-story
14
12 GM Set-4
10
8
Floor
6 o Benchmark
4 MPS
2 Median
ith GM
0
0 1 2 0 1 2 3 0 0.01 0.02 0 0.01 0.02 0.03
MPS procedure with benchmark EDPs; individual results for each of seven
scaled ground motions are also presented. Results are for the 6- (top row),
and 13-story building (bottom row).
50
6
GM Set-1 6-Story
5 GM Set-2
GM Set-3
4 GM Set-4
Floor
0
0.5 1 1.5 0 0.2 0.4 0.6 0.8
14
13-Story
12
10
8
Floor
0
0.5 1 1.5 0 0.2 0.4 0.6 0.8
 MPS
Figure 5.30 Median D MPS and dispersion  MPS of story drift ratios for four ground
motions sets and for two buildings.
51
6. EVALUATION OF MPS PROCEDURE:

HIGH-RISE BUILDINGS
Nonlinear RHA is being increasingly used to estimate seismic demands for high-rise buildings in
the U.S. This chapter tests the accuracy and effectiveness of the MPS procedure in scaling
ground motions to estimate the median seismic demands (with small dispersion) using two high-
rise moment-frame steel buildings.
6.1 BUILDINGS SELECTED
6.1.1 19-story Building

The building shown in Figure 6.1 is located in Century City – Los Angeles, designed in 1966-67
and constructed in 1967. It has 19 stories above ground and 4 stories of parking below the
ground level. The vertical load carrying system consists of 11.4 cm thick reinforced concrete
slabs supported on steel beams. There is no composite action between the slab and steel beams
due to lack of shear studs. The lateral load resisting system consists of four ductile steel moment
frames in the longitudinal direction and five X-braced frames in the transverse direction.
Moment resisting connections are used at the intersection of beams and columns. Perimeter
columns are standard I-sections except at the first story where columns are built-up box-sections.
The foundation system consists of 22 m long driven steel I-beam piles, capped in groups and
connected by 61 cm square reinforced concrete tie beams [Hart, 1973; Naeim 1998].
The building was initially instrumented with only three sensors at the time of the 1971
San Fernando earthquake. 15 sensors were in place during the 1994 Northridge earthquake
52
(Figure 6.2). During the Northridge event (its epicenter was 20 km away from the site), the
recorded peak horizontal accelerations were 0.32 g at the basement, 0.53 g at the ground floor
and 0.65 g at the roof; this intensity of shaking resulted in moderate damage in the building in
the form of buckling in some braces at upper floor levels in the transverse direction [Naeim,
1997]. There was no damage in the perimeter moment frames.

The second building selected is one of the tallest buildings in downtown Los Angeles, designed
in 1988 and constructed in 1988–90. This building comprises a 52-story steel frame office tower
and five levels of basement as underground parking. The floor plans of the tower are not
perfectly square; the tip of every corner is clipped and the middle third of each side is notched. In
groups of about five stories, above the 36 th story, the corners of the floors are clipped further to
provide a setback (Fig. 6.3).
The structural system of the building is composed of a braced-core, twelve columns

(eight on the perimeter and four in the core), and eight 91.4 cm deep outrigger beams at each
floor connecting the inner and outer columns. The core, which is about 17 m by 21 m, is
concentrically braced between the level-A (the level below the ground-level) and the 50 th storey.
Moment resisting connections are used at the intersection of beams and columns. The outrigger
beams, about 12 m long, link the four core columns to the eight perimeter columns to form a
ductile moment resisting frame. The outrigger beams are laterally braced to prevent lateral
torsional buckling and are effectively connected to the floor diaphragm by shear studs to transmit
the horizontal shear force to the frame. Perimeter columns are standard I-sections, while the core
columns are built-up sections with square cross section at the lower floor and crucifix section at
the upper floors. The interior core is concentrically braced. The building foundation is concrete
spread footings (2.74 to 3.35 m thick) supporting the steel columns with 13 cm thick concrete
slab on grade.
The building is instrumented with twenty accelerometers to record its translational and
torsional motions (Figure 6.4). During the Northridge earthquake (its epicenter was 30 km away
from the site), the recorded values of peak horizontal accelerations were 0.15 g at the basement,
0.17 g at the level-A, and 0.41 g at the roof; no structural damage was observed. The latest
53
recorded event, the 2008 Chino-Hills earthquake (its epicenter is 47 km away from the site),
generated a PGA of 0.06 g at the ground and 0.263 g at the roof level.
6.2 ANALYTICAL MODELS AND CALIBRATION STUDIES

Three dimensional (3-D) computer model of the building was developed. Steel columns, beams,
and braces were modeled by force-based nonlinear beam-column element in the open source
finite element platform, OpenSEES. To properly model the buckling response of the brace, a
slight imperfection is modeled at the brace’s mid-span. Centerline dimensions were used in
element modeling. The building weight, including estimates of non-structural elements such as
partition walls and the mechanical equipment in the roof, was estimated to be 102,309 kN. Nodes
at each floor were constrained to have the same lateral displacement to simulate rigid diaphragm
behavior. The estimated floor mass and mass moment of inertia were lumped at the centre of
mass at each floor. Panel zone deformations and local connection fracture were not considered
therefore modeling of members and connections was based on the assumption of stable
hysteresis loops derived from a bilinear stress-strain model with 3% strain hardening. The
expected yield stress for steel members equal to ~250 MPa was used. The columns were
assumed to be fixed at the base level. The P-D effects were included in the global system level.
As shown in Figure 6.5, the first six natural vibration periods of the building are
identified by frequency domain analysis of the motion of the roof relative to base motion
recorded during the Northridge event. The computer model was able to match the measured
periods (Table 6.1). Figure 6.6 plots the computed mode shapes for the first six modes; although
the motion is dominantly in the transverse direction, a slight torsional component exists.
Rayleigh damping was selected to be 3% of critical for the first and sixth modes (Fig. 6.7)
Nonlinear RHA of the building subjected to the three components of the motion recorded at the
base level during the Northridge event leads to the relative displacement response in two
horizontal directions at the roof, eight and second floors shown in Figure 6.8, where it is
54
compared with the motions derived from records. The excellent agreement between the
computed and recorded displacements indicates that the computer model is adequate.

The 3-D model of the 52-story building in OpenSEES included 58 separate column types and 23
different beam types. The building weight, including estimates of non-structural elements such as
partition walls and the mechanical equipment in the roof, was estimated to be 235,760 kN. The
material model is based on stable hysteresis loops derived from a bilinear stress-strain model
with 2% strain hardening. All steel framing including columns are ASTM A-572 (grade 50) with
a nominal yielding strength of 345 MPa.
The first nine natural vibration periods of the building identified from the recorded
relative roof motions with respect to base during the Northridge and Chino-Hills events are
shown in Figure 6.9. the computer model was able to match the measured periods reasonable
well (Table 6.2). Figure 6.10 plots the mode shapes for the first nine modes; torsional modes are
well separated from translational modes. Rayleigh damping was selected to be 4% of critical for
the first and ninth modes (Fig. 6.11). Nonlinear RHA of the building subjected to the three
components of the motion recorded at the base level during the Northridge event leads to the
relative displacement response in two horizontal directions at the roof, thirty five, twenty second
and fourteenth floor shown in Figure 6.12. The excellent agreement between the computed and
recorded displacements indicates that the computer model is adequate.
6.3.FIRST-MODE SDF-SYSTEM PARAMETERS
The force-deformation relation FS1 / L1 − D1 for the first-mode SDF system for each building is
determined from the base shear – roof displacement relation defined by first-mode pushover
curve by utilizing Fs1/L1= Vb1 / M 1 and D1 = ur1 / 1fr1. For both buildings, only the E-W direction
*
is considered in the pushover analysis. The resulting force-deformation relations for the first-
mode SDF system are shown in Figure 6.13. For the 19-story building, a bilinear hysteretic
force-deformation relation is found to be adequate, while for the 52-story building, the hysteretic
55
force-deformation relation is idealized by the peak-oriented model [Ibarra and Krawinkler, 2005;
Ibarra et al., 2005] with the monotonic curve idealized as tri-linear.
6.4.EVALUATION OF MPS PROCEDURE
The efficiency and accuracy of the MPS procedure will be evaluated for each of three ground
motion sets (Table 5.1), separately by comparing the median values of EDPs determined by
nonlinear RHA of each building due to the 7 scaled records against the benchmark EDPs.
Although the 21 records selected are the most intense records in the NGA database
consistent with the seismic hazard defined in Chapter 4.3, the 19- and 52-story buildings remain
essentially within their linearly elastic range. To impose significant inelastic deformations, these
21 ground motion records were amplified by a factor of two and their median response spectrum
(Figure 6.14) is treated as the target spectrum. Nonlinear RHAs under scaled records indicates
significant nonlinearity from both buildings as shown in Figure 6.15, where the resultant peak
roof displacement values are identified on the first-mode pushover curve; and their median
values are marked as vertical dashed line. Benchmark EDPs for two buildings are presented in
Figure 6.16; also included are the results from individual records to demonstrate the large
dispersion.
I
In the MPS procedure, a target value of deformation, D1 is estimated through the
empirical CR equation (Step 6 of the procedure - Chapter 2.1). Figure 6.17 compares the “exact”
I
target value of deformation, D1 (dashed horizontal line) with estimated target value of
I I
deformation, D1 (continuous horizontal line) using the C R equation; D1 values from individual
records are also included to show its record-to-record variability. Notably, “exact” target value of
I I
deformation, D1 is determined as the median value of D1 from 21 ground motions based on
I I
first-mode inelastic SDF system. Estimated value of D1 reasonably matches “exact” value D1
with slight overestimation for both buildings.
56
I
Once D1 is estimated, an appropriate scale factor for each record is determined based on
the inelastic first-mode SDF systems (Fig. 6.13) by implementing Steps 7-8 of the MPS
procedure. The scale factors established for each record in each of the three sets are presented in
Table 6.3 for both buildings. The EDPs determined by nonlinear RHAs due to three sets of 7
ground motions scaled according to the MPS procedure are compared against the benchmark
EDPs. Figures 6.18 and 6.19 present such comparisons respectively for the 19- and 52-story
buildings considering the three sets of ground motions. As shown, the median values of the peak
floor displacement, story drift ratio and column plastic rotations are well-estimated (difference is
less than 20%), whereas beam plastic rotations for both buildings are slightly overestimated at
intermediate floors. The dispersion of the EDP values due to the 7 scaled records is reduced as
compared to the dispersion associated with the original records (Fig. 6.16), but this reduction is
less at intermediate and upper floors, indicating that higher-mode contributions to structural
response should also be considered in identifying the most appropriate scaled records.
6.5.COMPARING MPS AND CODE-BASED SCALING PROCEDURES
In addition to the MPS, scale factors for each of 7 records in three sets are determined in
accordance with the ASCE -7 procedure (described in Section 3). The EDPs determined by
nonlinear RHA of the structure due to a set of 7 ground motions scaled according to the ASCE 7-
05 scaling procedure are compared against the benchmark EDPs. Figures 6.20 - 6.21 present
such comparisons for the two buildings and for the three sets of 7 ground motions.
As opposed to the low- and mid-rise buildings (Chapter 5), the ground motions scaled
according to the ASCE 7-05 scaling method reasonably estimate the median benchmark EDPs
for both buildings; however the dispersion of EDPs is noticeably large for all three sets for the
19-story building and for the first two sets for the 52-story building. The dispersion in the EDP
values due to the 7 scaled records around the median value is much smaller when the records are
scaled according to the MPS procedure compared to the ASCE 7-05 scaling procedure.
An alternative way of comparing MPS and ASCE 7-05 scaling methods is based on the
ratio of the EDP value due to a scaled record and the benchmark value. The deviation of the
median D of this ratio from unity is an indication of the error or bias in estimating the median
EDP value, and the dispersion  of this ratio is an indication of the scatter in the individual
57
EDPs, determined from the scaled ground motions. Figure 6.22 presents the median D of the
EDP ratio for story drifts determined from records scaled according to the MPS and ASCE 7-05
scaling methods. Comparing these D values against 1.0, it is apparent that the MPS method is
more accurate. The bias in the MPS methods is generally less than 25%. The bias associated with
the ASCE 7-05 method is generally much higher in all cases. As it is evident in Figure 6.23 that
the MPS scaling method leads to the smallest dispersion; dispersion is much larger in the ASCE
7-05 scaling method for 19-story building and unacceptably high for 52-story building for the
first two sets of records.
6.6.MULTI-MODE CONSIDERATIONS
It is well known that the contribution of higher-modes is more significant in the structural
response of longer-period buildings. These higher vibration modes are therefore considered in
selecting the most appropriate set of 7 ground motions out of the 21 records that were scaled
initially based on the first-mode response only. This is achieved by implementing Steps 10-13 of
the MPS method. The seven records with the highest ranks (see Step 12) were defined as Ground
Motion Set 4. Note that this set is different for each building (Table 6.4).
Figure 6.24 compares the median EDPs from ground motions scaled by the MPS method
with the benchmark values for the 19- and 52-story buildings; results of individual ground
motions are also shown. Considering higher modes in selecting ground motions in the MPS
method provides accurate estimate of the median EDPs and reduces record-to-record variability
(compared to the results achieved by Ground Motion Sets 1-3). These enhanced accuracy and
efficiency are demonstrated in Figure 6.25, where the D and —the median value of the ratio of
the estimated story drift to its benchmark value, and dispersion of this ratio—are plotted for the
four set of ground motions. Ground Motion Set 4 is more accurate and efficient than Ground
Motion Sets 1 though 3; this improvement is more pronounced for the 52-story building.
Table 6.1 Measured and computed natural periods for 19-story building
58
Table 6.2 Measured and computed natural periods for 52-story building
Table 6.3 Scale factors for two buildings and for three sets of seven ground motions
59
Table 6.4 Scale factors for two buildings considering higher-modes
60
Figure 6.1 Overview of the 19-story building in Century City, Los Angeles, CA.
61
Figure 6.2 Instrumentation lay-out of the 19-story building.
62
Figure 6.3 Overview of the 52-story building in downtown, Los Angeles, CA.
63
Figure 6.4 Instrumentation lay-out of the 52-story building.
E-W Direction N-S Direction Torsional Direction

0.2
Fourier Amp. of Rel. Acceleration
1994 Northridge Eq.

3.7 sec
3.4 sec
0.8 sec
1.4 sec
0.6 sec
1.1 sec
3.7 sec
0.8 sec
1.1 sec
0
0 1 2 3 4 50 1 2 3 4 50 1 2 3 4 5
Period (sec)
Figure 6.5 Identification of natural periods for the 19-story building using recorded
relative motion at the roof level during the 1994 Northridge earthquake.
64
E-W Component N-S Component Torsional Component
18 Mode 3 Mode 5 Mode 2
16
[1,3,5]
14
Mode 1 [1,3,6]
12 Mode 5
10
Floor
8 Mode 4
[2,4,6]
6 Mode 4
4 Mode 2
2 Mode 6
0
-0.2 -0.1 0 0.1 0.2 -0.2 -0.1 0 0.1 0.2 -5 0 5
   -5
xj yj j x10
Figure 6.6 Natural vibration modes of the 19-story building.
6
Rayleigh damping
Modal Damping Ratio (%)
2
1st-Mode
6th-Mode
0
0 5 10
Natural Frequency, (wn)
Figure 6.7 Modal damping ratios for the 19-story building.
65
E-W Direction N-S Direction

50
th
19 Floor
0
OpenSees
Observed
-50
50
Displacement (cm)
th
8 Floor
-50
50
nd
2 Floor
-50
0 10 20 30 40 00 10 20 30 40 50
Duration (sec)
Figure 6.8 Comparison of observed and computed floor displacements in two

horizontal directions of the 19-story building at different floor levels.
(Recorded data is from the M6.7 1994 Northridge earthquake).
66
E-W Direction N-S Direction Torsional Direction

0.1
1.7 sec
(x2)
1.8 sec
1994 Northridge Eq.
0.9 sec
1.7 sec
0.9 sec
5.6 sec
Fourier Amp. of Rel. Acceleration
5.9 sec
0.9 sec
4.7 sec
0
0.04
(x2)
1.7 sec
2008 Chino-Hills Eq.

0.9 sec
1.7 sec
1.7 sec
0.9 sec
0.9 sec
5.3 sec
5.6 sec
0
0 1 2 3 4 5 6 7 1 2 3 4 5 6 7 1 2 3 4 5 6 7
Period (sec)
Figure 6.9 Identification of natural periods for the 52-story building using recorded
relative motion at the roof level (for E-W and N-S directions) and 49 th floor
level (for torsional direction) during the 1994 Northridge (upper-panels)
and 2008 Chino-Hills (lower-panels) earthquakes.
67
E-W Component N-S Component Torsional Component
50
[1,3,5,6,8,9]
45
40 Mode 1 Mode 3
35 9
30 [1,2,4,5,7,9] Mode 2
Floor
25
5 6
4
20
15 8 7
[2,3,4,6,7,8]
10
0
-0.2 -0.1 0 0.1 0.2 -0.2 -0.1 0 0.1 0.2 -2 -1 0 1 2
   -4
xj yj j x10
Figure 6.10 Natural vibration modes of the 52-story building.
Rayleigh damping
6
Modal Damping Ratio (%)
2
1st-Mode
9th-Mode
0
0 2 4 6 8 10
Natural Frequency, (wn)
Figure 6.11 Modal damping ratios for the 52-story building.
68
E-W Direction N-S Direction

3
Roof OpenSees
Observed
0
-3
3
35th Floor
0
Displacement (cm)
-3
3
22nd Floor
-3
3
14th Floor
-3
15 20 25 30 35 40 4515 20 25 30 35 40 45
Duration (sec)
Figure 6.12 Comparison of recorded and computed floor displacements in two

horizontal directions of the 52-story building at different floor levels.
(Recorded data is from the M5.4 2008 Chino-Hills earthquake).
69
0.3 0.05
19-Story Bldg. 52-Story Bldg.
0.04
0.2
0.03
FS1 / L1 (g)
FS1 / L1 (g)
0.02
0.1
0.01
0 0
0 50 100 150 200 250 0 100 200 300
D1 (cm) D1 (cm)
model (dashed line) for the 19- and 52-story buildings.
7
Median
6 ith GM
Pseudo Acceleration (g)
0
0 2 4 6 8 10
Period (s), Tn
Figure 6.14 Individual response spectra for 21 ground motions (up-scaled by 2) and
their median response spectrum; = 5%.
70
0.08 0.03
19-story Bldg. 52-story Bldg.
0.06
Base Shear / Weight
Base Shear / Weight

0.02
0.04
Median
Median 0.01
0.02
0 0
0 1 2 3 4 5 0 1 2
Roof Disp. / Bldg. Height (%) Roof Disp. / Bldg. Height (%)
Figure 6.15 Roof displacements determined by nonlinear RHA of the 19- and 52-story
buildings for 21 ground motions identified on first-mode pushover curves.
15
Floor
10
5
Median
ith GM
0
0 2 4 0 2 4 6 0 0.02 0.04 0 0.02 0.04 0.06
50Floor Disp. / Bldg. Height (%) Story Drift Ratio (%) Col. Pl. Rot. (rad) Beam Pl. Rot. (rad)
40
30
Floor
20
10 Median
ith GM
0
0 1 0 2 0 0.01 0 0.01 0.02
Figure 6.16 Median values of EDPs determined by nonlinear RHA of 19- (top row) and
52-story (bottom row) buildings for 21 ground motions; results for
individual ground motions are also included.
71
19-Story Building 52-Story Building

500 500
CR Equation
Median
400 400
ith GM
300 300
D1 (cm)
D1 (cm)
I
I
200 200
100 100
5 10 15 20 5 10 15 20
Ground Motion Number Ground Motion Number
I
ground motions for the 19- and 52-story buildings; “exact” target value of
I
deformation D1 is identified by horizontal dashed line; horizontal
I
continuous line indicates target value of deformation D1 established by CR
equation.
72
GM Set-1
15
Floor
10
o Benchmark
5 MPS
Median
ith GM
0
GM Set-2
15
Floor
10
GM Set-3
15
Floor
10
0
0 2 4 0 5 0 0.02 0 0.02 0.04
Figure 6.18 Comparison of median EDPs based on the MPS procedure with benchmark
73
50
GMSet-1
40
30
Floor
20 o Benchmark
MPS
10 Median
ithGM
0
50
GMSet-2
40
30
Floor
20
10
0
50
GMSet-3
40
30
Floor
20
10
0
0 1 2 0 5 0 0.02 0 0.02 0.04
74
GM Set-1
15
Floor
10
o Benchmark
5 ASCE-7
Median
ith GM
0
GM Set-2
15
Floor
10
GM Set-3
15
Floor
10
0
0 2 4 0 5 0 0.02 0 0.02 0.04
Figure 6.20 Comparison of median EDPs based on the ASCE 7-05 ground motion
scaling procedure with benchmark EDPs for the 19-story building;
individual results for each of the seven scaled ground motions are also
presented.
75
50
GM Set-1
40
30
Floor
20 o Benchmark
ASCE-7
10 Median
ith GM
0
50
GM Set-2
40
30
Floor
20
10
0
50
GM Set-3
40
30
Floor
20
10
0
0 1 2 0 5 0 0.02 0 0.02 0.04
scaling procedure with benchmark EDPs for the 52-story building;
individual results for each of the seven scaled ground motions are also
presented.
76
ASCE-7
15 MPS
Floor
10
19-story
0 GM Set-1 GM Set-2 GM Set-3
0 0.5 1 1.5 0 0.5 1 1.5 0 0.5 1 1.5 2
50
40
30
Floor
20
10
52-story
0
0 0.5 1 1.5 0 0.5 1 1.5 0 0.5 1 1.5 2
or
Figure 6.22 Median story drift ratios D MPS and D ASCE - 7 for three ground motions sets
and for two buildings.
77
ASCE-7
15 MPS
Floor
10
19-story
0 GM Set-1 GM Set-2 GM Set-3
0 0.5 1 0 0.5 1 0 0.5 1 1.5
50
40
30
Floor
20
10
52-story
0
0 0.5 1 0 0.5 1 0 0.5 1 1.5
 MPS or  ASCE -7
Figure 6.23 Dispersion of story drift ratios  ASCE -7 and  MPS for three ground
motions sets and for two buildings.
78
GM Set-4
15
Floor
10
o Benchmark
5 MPS
Median
ith GM
0
0 2 4 0 5 0 0.02 0 0.01 0.02
50 Fl. Disp. / Bldg. Height (%) Story Drift Ratio (%) Col. Pl. Rot. (rad) Beam Pl. Rot. (rad)
GM Set-4
40
30
Floor
20
10
0
0 1 2 0 5 0 0.02 0 0.02 0.04
MPS procedure with benchmark EDPs; individual results for each of seven
scaled ground motions are also presented. Results are for the 19- (top
row), and 52-story building (bottom row).
79
15
Floor
10
GM Set-1
5
GM Set-2
GM Set-3
GM Set-4 19-Story Bldg.
0
0 0.5 1 1.5 2 0 0.2 0.4 0.6 0.8 1
50
40
30
Floor
20
10
52-Story Bldg.
0
0 0.5 1 1.5 2 0 0.2 0.4 0.6 0.8 1
D MPS  MPS
Figure 6.25 Median D MPS and dispersion  MPS of story drift ratios for four ground
motions sets and for 19- and 52-story buildings.
80
7. EVALUATION OF MPS PROCEDURE:

SHORT - PERIOD BUILDING
The efficiency and accuracy of the MPS procedure were established in previous sections
considering structures with first-mode vibration period ranging from 1.0 to 5.6 sec. This chapter
further assesses the MPS procedure considering a building with a shorter first-mode vibration
period (T1 = 0.6 s). The short period range (T 1 << 1.0 s) has a particular importance in predicting
the inelastic target displacement values because the equal displacement rule is often violated in
this range meaning that inelastic displacement values become larger than the elastic ones.
7.1 BUILDING DETAILS AND ANALYTICAL MODEL
The 4-story building (Fig. 7.1), located in Watsonville, California is selected for the assessment
of the MPS procedure. It was originally designed and constructed in 1948 as a three-story
building and the fourth story was added in 1955. Its vertical carrying system consists of concrete
slabs supported by composite concrete-steel columns. The lateral load system consists of
concrete shear-walls in both directions. The foundation system consists of spread footings below
shear walls.
The building was instrumented in 1983 with 13 strong motion sensors on three levels; the
instrumentation lay-out is shown in Figure 7.2. During the 1989 Loma Prieta earthquake, peak
horizontal accelerations of 0.66 g at the ground level and 1.24 g at the roof level were recorded
in this building. This recorded data is used to calibrate the 3-D computer model developed in
OpenSEES. Centerline dimensions were used in element modeling. The beams, columns, and
shear walls were modeled with nonlinear-beam-column element based on fiber sections
containing confined concrete, unconfined concrete and steel reinforcing bars. This material
model enabled capturing of the rapid strength loss after the building’s peak strength. The
crushing strain of the unconfined concrete was selected to be equal to 0.004 and that for confined
concrete was selected to be that corresponding to the rupture of confining steel based on the
81
Mander model. The strength of concrete and steel was selected based on the values specified in
the structural drawings. Nodes at each floor were constrained to have the same lateral
displacement to simulate rigid diaphragm behavior. The estimated floor mass and mass moment
of inertia were lumped at the centre of mass at each floor. The P-D effects in the global system
level were included in modeling. Also considered in the model is the foundation flexibility at the
base as per the FEMA-356 [ASCE 2000] recommendations by attaching six linear springs —three
along the x-, y-, and z-translation, two about the x- and y-rocking, and one about the z-torsion
because the lateral load resisting system of the building consists of shear walls in both directions.
7.2.FIRST-MODE SDF-SYSTEM PARAMETERS
The force-deformation relation FS1 / L1 − D1 for the first-mode SDF system is determined from
the base shear – roof displacement relation defined by first-mode pushover curve by utilizing
Fs1/L1= Vb1 / M 1 and D1 = ur1 / 1fr1. Only the E-W direction is considered in the pushover
*
analyses. The resulting force-deformation relation for the first-mode SDF system is shown in
Figure 7.3 where the hysteretic force-deformation relation is idealized as tri-linear.
Similar to the other buildings, the efficiency and accuracy of the MPS procedure will be
evaluated for each of three ground motion sets (Table 5.1), separately by comparing the median
values of EDPs determined by nonlinear RHA of 4-story building under 7 scaled records against
the benchmark EDPs. Nonlinear RHAs under unscaled records indicates significant nonlinearity
as shown in Figure 7.4, where the resultant peak roof displacement values are identified on the
first-mode pushover curve; and their median values are marked as vertical dashed line.
Benchmark EDPs are presented in Figure 7.5; also included are the results from individual
records to demonstrate the large dispersion.
I
empirical CR equation (Step 6 of the procedure - Chapter 2.1). Figure 7.6 compares the “exact”
82
I
target value of deformation, D1 (dashed horizontal line) with estimated target value of
I I
deformation, D1 (continuous horizontal line) using the C R equation; D1 values from individual
records are also included to show its record-to-record variability. Notably, “exact” target value of
I I
deformation, D1 is determined as the median value of D1 from 21 ground motions based on
I
first-mode inelastic SDF system. Estimated value of D1 underestimates the “exact” target value
I
of deformation by only 6%. Once D1 is estimated, an appropriate scale factor for each record is
determined based on the inelastic first-mode SDF systems (Fig. 7.7) by implementing Steps 7-8
of the MPS procedure. The EDPs determined by nonlinear RHAs due to three sets of 7 ground
motions scaled according to the MPS and ASCE 7-05 methods are compared against the
benchmark EDPs. Figures 7.7 and 7.8 present such comparisons considering the three sets of
ground motions. As shown, the median values of the peak floor displacement, story drift ratio
and column plastic rotations are better estimated (difference is less than 10%) using the MPS
procedure, same applies for beam and column plastic rotations. The dispersion of the EDP values
due to the 7 scaled records is reduced as compared to the dispersion associated with the original
records (Fig. 7.6), but this reduction is less at the first floor.
Similar to the previous chapters, an alternative way of comparing MPS and ASCE 7-05
scaling methods is based on the ratio of the EDP value due to a scaled record and the benchmark
value. The deviation of the median, D of this ratio from unity is an indication of the error or bias
in estimating the median EDP value, and the dispersion,  of this ratio is an indication of the
scatter in the individual EDPs, determined from the scaled ground motions. Figure 7.9 presents
the median, D of the EDP ratio for story drifts determined from records scaled according to the
MPS and ASCE 7-05 scaling methods. Comparing these D values against 1.0, it is apparent that
the MPS and ASCE 7-05 methods have the same accuracy for the 4-story (short-period) building.
The bias in the MPS methods is generally less than 30%, except in the case the Ground Motion
Set-2 at upper level and Set-3 at the first floor. Figure 7.10 presents the dispersion of the EDP
ratio for story drifts determined from records scaled according to the MPS and ASCE 7-05
scaling methods. It is apparent that the dispersion associated with both methods is similar.
Comparison of there sets implies that dispersion varies significantly from one set to other; this
83
observed variation at short period building was not as large for the other buildings, which have
T1 >1.0 sec.
Table 7.1 Scale factors computed for three sets of seven ground motions
84
Figure 7.1 Overview of the 4-story RC building in Watsonville, CA.
Figure 7.2 Instrumentation lay-out of the 4-story RC building.
85
0.5
0.4
FS1 / L1 (g)
0.3
0.2
0.1
0
0 5 10 15 20
D1 (cm)
model (dashed line) for the 4-story RC building.
0.5
0.4
Base Shear / Weight
0.3
Median
0.2
0.1
0
0 1 2
Figure 7.4 Roof displacements determined by nonlinear RHA of the 4-story RC

building for 21 ground motions identified on first-mode pushover curve.
86
4
4-Story
3
Floor
1 Median
ith GM
0
0 1 2 0 1 2 3 0 0.01 0.02 0 0.01 0.02 0.03
Figure 7.5 Median values of EDPs determined by nonlinear RHA of 4-story RC building
for 21 ground motions; results for individual ground motions are also
included.
50
CR Equation
40 Median
ith GM
30
D1 (cm)
I
20
10
0
0 4 8 12 16 20
I
I
ground motions; “exact” target value of deformation D1 is identified by
horizontal dashed line; horizontal continuous line indicates target value of
I
deformation D1 established by CR equation.
87
4
GM Set-1
3
Floor
2
o Benchmark
1 MPS
Median
ith GM
0
4
GM Set-2
3
Floor
0
4
GM Set-3
3
Floor
0
0 1 20 2 40 0.01 0.02 0 0.01 0.02 0.03
88
4
GM Set-1
3
Floor
2
o Benchmark
1 ASCE-7
Median
ith GM
0
4
GM Set-2
3
Floor
0
4
GM Set-3
3
Floor
0
0 1 20 2 40 0.01 0.02 0 0.01 0.02 0.03
scaling procedure with benchmark EDPs for the 4-story building; individual
results for each of the seven scaled ground motions are also presented.
89

4
ASCE-7
MPS
3
Floor
4-story
0
0 0.5 1 1.5 0 0.5 1 1.5 0 0.5 1 1.5 2
or
Figure 7.9 Median story drift ratios D MPS and D ASCE - 7 for three ground motions sets.
3
Floor
0
0 0.5 1 0 0.5 1 0 0.5 1 1.5
or
Figure 7.10 Dispersion of story drift ratios  MPS and  ASCE - 7 for three ground motions
sets.
90
8.EVALUATION OF MPS PROCEDURE:

BRIDGES
The seismic design of typical highway bridges in the U.S. is governed by the American
Association of Highway and Transportation Officials bridge design specifications (AASHTO-
2004) or Seismic Design Criteria (SDC-2006) in California. These specifications permit the
designer to utilize a variety of methods for seismic analysis. According to the SDC-2006,
equivalent static analysis and linear elastic dynamic analysis are the appropriate analytical tools
for estimating the displacement demands, and pushover analysis is the appropriate analytical tool
for establishing the displacement capacities. Nonlinear RHA is yet to be included in the SDC.
Next generation of the SDC is expected to incorporate nonlinear RHA similar to current building
codes for validation of new design and performance assessment of existing bridges.
Ground motion records for nonlinear RHA need to be appropriately selected and scaled
for accurate estimates of EDPs with small dispersion. In this chapter, the accuracy and
effectiveness of the MPS procedure in estimating median EDPs will be demonstrated on two
ordinary standard bridges. “Ordinary Standard Bridge” is a term used by the California
Department of Transportation (Caltrans) to identify bridges which can be designed using the
direct and simplified approach outlined in the SDC-2006. For a bridge to be considered as an
Ordinary Standard Bridge, SDC-2006 sets forth the following basic requirements: (1) the span
length should be less than 90 m, (2) the bridge should be constructed with normal weight
concrete, (3) foundations must be supported on spread footings, pile caps with piles or pile
shafts, and (4) the soil is not susceptible to liquefaction or lateral spreading during strong
shaking.
8.1 DESCRIPTION OF BRIDGES AND ANALYTICAL MODELS
Two types of existing ordinary bridges located in California have been considered: single bent
two-span overpass bridge and multi-span single-frame bridge.
91
8.1.1 Single-bent Two-span Overpass Bridge
Selected single-bent two-span overpass bridge is a representative of an overcrossing designed

according to post-Northridge Caltrans specifications. This bridge consists of a two-span
continuous deck and single-bent composed of two octagonal columns with spiral reinforcement.
Elevation, deck and column dimensions and column section detailing are shown in Figure 8.1.
The bridge has stub wall abutments which are restrained in the longitudinal and transverse
directions as a result of end diaphragm and wing wall interaction with the soil.
The computer model of the bridge was created in OpenSEES finite element platform.
Due the abutments, the bridge was assumed to be fixed against vertical, translation and rotation
about the longitudinal axis of the superstructure but longitudinal and transverse displacements
were permitted. The column bent footings were supported on translational and rotational springs
in each orthogonal direction. The stiffness of the translational springs to model the abutments
was determined in accordance with the Caltrans bridge design specifications (SDC-2006). The
abutment stiffness in the longitudinal direction was computed as:
K = (Ki.w.h) / 1.7 (8.1)
where Ki is the initial stiffness of the abutment and was taken to be equal to 11.5 kN/mm per
meter. w and h are the width and height of the abutment’s diaphragm. For the abutment stiffness
in the transverse direction and foundation stiffness in both translational directions, an empirical
value of 7,000 kN/m per pile was used to account for soil flexibility.
The finite element model of the bridge is represented by 3-D frame elements passing
through the mid-depth of the superstructure and 3-D frame elements passing through the
geometric centre and mid-depth of the columns and the cap beam. Two frame elements were
used for modeling columns (Fig. 8.2). Since the cross-section of the superstructure is uniform, it
was deemed sufficient to locate 8 nodes within each span. The mass of each member was equally
distributed to its end nodes. The seismic weight of the bridge was approximated to be 11.5 MN.
The box-girder was assumed to be integral with the bent, thus full continuity was employed at
'
the superstructure-bent intersection. The characteristic strength of unconfined concrete ( f c ) was
assumed to be 28 MPa with an ultimate strain of 0.005, while the yield strength of both
longitudinal and spiral reinforcement was specified as 413 MPa. The characteristic strength and
92
ultimate strain of confined concrete were computed to be 41.6 MPa and 0.0169 using the model
proposed by Mander et al. (1988). A bilinear model with a postyield stiffness of 1% of the initial
stiffness was used to model reinforcing steel. For simplicity, the columns were modeled as
circular columns with a diameter of 1.65 m. The diameter of the equivalent circular column was
selected such that both core and cover concrete areas of the original octagonal column are equal
to those of the equivalent circular column. The columns were modeled using 3-D, fiber-
discretized, nonlinear beam-column elements. The deck elements were assumed to remain
elastic.
8.1.2 Multi-span Single-frame Bridge
The multi-span single-frame bridge selected is representative of typical multi-span, single-frame

prestressed concrete bridges built according to post-Northridge Caltrans design specifications.
Figure 8.3 shows the bridge elevation, deck and column dimensions, and column section
detailing. It has four spans with 12.34 m and 16.45 m in length; bents are 19.5 m, 28 m and 25.25
m in height. The bents of the bridge consist of single double-spiral columns.
The multi-span single-frame bridge was also modeled as an elastic superstructure sitting
on nonlinear columns on elastic foundation. The assumption of elastic superstructure was based
on the capacity design approach employed by SDC-2006. An overview of the computer model of
the bridge is shown in Figure 8.4. The compressive strength of unconfined concrete and the yield
strength of longitudinal reinforcement were taken as 28 MPa and 413 MPa, respectively. The
compressive strength and ultimate strain of confined concrete were computed as 40.2 MPa and
0.0157 using Mander’s model. A bilinear model with a post-yield stiffness of 1% of the initial
stiffness was used to model reinforcing steel.
The columns of the bridge rest on shallow foundations. Six elastic springs in 3
translational and 3 rotational directions were used to model the soil effect. The approximate
expressions proposed in FEMA-356 (FEMA 2000) were used to compute properties of the
corresponding spring constants.
Seat type abutments are used at both ends of the bridge. Spring systems were used to
model the dynamic stiffness of the abutments. In the transverse direction, shear keys are
93
designed to break off during a strong ground motion. Hence, seat type abutments do not possess
any stiffness in the transverse direction. In the vertical direction, the movement of the bridge is
prevented at the abutments in both upward and downward directions. Thus, the abutments were
modeled as restrained supports in the vertical direction. In the longitudinal direction, the bridge
is free to move in the opposite direction of the abutment at each end. Towards the abutment there
is a certain amount of gap before the deck makes contact with the abutment. When the deck and
the abutment are in contact, the stiffness of the abutment is computed as:
K = (Ki.w.h) / 5.5 (8.2)
where Ki is taken to be equal to 3.5 kN/mm per meter. A spring which has no stiffness in tension
and remains elastic in compression with an initial stiffness of 10,060 kN/m and with a 10 cm gap
was used to model the abutment behavior in the longitudinal direction.
8.2 FIRST-MODE SDF-SYSTEM PARAMETERS
Figure 8.5 shows the first two modes of vibration and their periods for both bridges. For the
single-bent two-span overpass bridge, the first mode with a period 0.54 sec involves a transverse
translation of the deck; the second mode with a period 0.52 sec involves a longitudinal
translation of the superstructure. The multi-span single-frame bridge has the first mode period of
2.47 sec in the translational direction and 1.06 sec in the longitudinal direction.
Modal pushover curves for the bridge were developed in the transverse and longitudinal
directions separately. Similar to the modal pushover analyses procedure for buildings (Chopra,
2001), the distribution of lateral forces was determined from the shape of the fundamental
transverse mode, and fundamental longitudinal mode, multiplied by tributary mass. For each
direction, the resultant pushover curves for both bridges are demonstrated in Figure 8.6 with a
thick solid line. In the same figure, bilinear idealization of pushover cures is shown in thick
dashed line. For the single-bent two-span overpass bridge the pushover curves are quite similar
in the longitudinal and transverse directions.
94
As in previous chapters, the efficiency and accuracy of the MPS procedure are evaluated for each
of the three sets of 7 ground motions (Table 5.1), separately by comparing the median values of
EDPs determined by nonlinear RHA of the bridges due to 7 scaled records against the
benchmark EDPs. The MPS procedure is applied based only on first-mode of vibration (Section
2.1) and only for estimating response to transverse excitation. EDPs adopted are global response
parameters: peak value of deck drift ratio and bent-column plastic rotation. The benchmark EDPs
are computed first via nonlinear RHA due to 21 unscaled ground motions; the resultant median
EDPs are presented in Figure 8.7 together with results from individual records. As expected,
unscaled records led to large dispersion of EDPs. The peak values of deck drifts due to 21
ground motions range from 1% to over 4%; these values are identified on the first-mode
pushover curves together with their median values in Figure 8.8. These comparisons indicate that
all of the excitations drive both bridges well into the inelastic range.
I
empirical CR equation (Step 6 of the procedure - Chapter 2.1) using post-yield stiffness ratio, and
yield-strength ratio. Yield-strength ratio, R y is determined (Step 4 of the MPS procedure) from:
R y = M1* A1 Vb1 y , where M1* is the effective modal mass and Vb1 y is the yield point value of
base shear determined from the idealized pushover curve (Fig. 8.6). R y is computed as 3.10 and
3.46, respectively for single-span and multi-span bridges. Figure 8.9 compares the estimated
I
target value of deformation, D1 using the empirical CR equation against the exact median value.
Exact value is computed from nonlinear RHAs of inelastic first-mode SDF systems with force-
I
deformation relation shown in Figure 8.6 due to 21 unscaled records. D1 values from individual
records are also included to show its record-to-record variability. As Figure 8.9 demonstrates the
I I
empirical value of D1 overestimates the “exact” value D1 by 12% to 14%.
I
Once D1 is estimated, Steps 7 and 8 of the MPS procedure are implemented to determine
an appropriate scale factor for each record. The scale factors established for each record in each
of the three sets are presented in Table 8.1 for both bridges. The EDPs determined by nonlinear
RHAs of bridges due to the three sets of 7 ground motions scaled according to MPS procedure
95
are compared against the benchmark EDPs. Figures 8.10 and 8.11 exhibit these comparisons for
the single-bent overpass bridge and multi-span bridge, respectively. The estimate of median
value of EDPs determined by the MPS procedure are within 90% to 112% of the benchmark
values for the single-bent overpass. For the multi-span bridge, median EDPs are overestimated as
much as 18% for the ground motion Set-2 and 3; the overestimation is less for the Set-1. For
both bridges, the dispersion of EDPs due to scaled records is reduced significantly as compared
to the dispersion due to the unscaled records (Fig. 8.7). These results demonstrate that the MPS
procedure leads to scaled ground motions that yield accurate estimate of median EDPs with
significantly reduced dispersion compared to the unscaled ground motions.
Table 8.1 Scale factors computer for ground motion records
96
30.95 m 30.52 m
A A
Figure 8.1 Elevation, deck and column dimensions, and column section details of
single-bent overpass bridge.
97
L2
L1
h
L1 = 30.95 m
L2 = 30.52 m
V
S = 8.80 m
T L h = 8.50 m
Figure 8.2 Computer model of single-bent overpass bridge.
98
A A
Figure 8.3 Elevation, deck and column dimensions, and column section details of
multi-span bridge.
99
h3
h2
L1 = 40.50 m
h1 L2 = 54.00 m
h1 = 19.75 m
h2 = 28.00 m
h3 = 25.25 m
Figure 8.4 Computer model of multi-span bridge.
T1 = 0.54s T2 = 0.52s
T1 = 2.47s T2 = 1.06s
Figure 8.5 First two modes of vibration and their periods for single-bent overpass
bridge [top] and multi-span bridge [bottom].
100
0.5 0.5
0.4 Ti = 0.54 s 0.4 Ti = 0.52 s

a = 0.005 a = 0.005
0.3 0.3
FS1 / L1 (g)
FS1 / L1 (g)
0.2 0.2
Teff = 0.86 s Teff = 0.82 s
0.1 0.1
Transverse Dir. Longitudinal Dir.
0 0
0 5 10 15 20 0 5 10 15 20
D1 (cm) D1 (cm)
0.25 0.7
Teff = Ti = 1.06 s
Teff = Ti = 2.47 s
0.6
0.2 a = 0.07
0.5
0.15
FS1 / L1 (g)
= 0.18
FS1 / L1 (g)
a
0.4
0.1 0.3
0.2
0.05
0.1
Transverse Dir.
Longitudinal Dir.
0 0
0 20 40 60 80 100 0 20 40 60 80 100
D1 (cm) D1 (cm)
bilinear model (dashed line) for single-bent overpass bridge [top-panels];
multi-span bridge [bottom-panels].
101
5 0.04
Median
ith GM
Column Plastic Rot. (rad)

Deck Drift Ratio (%) 4
0.03
3
0.02
2
0.01
1
0 0
5 10 15 20 5 10 15 20
5 0.04
Median
ith GM Column Plastic Rot. (rad)
4
Deck Drift Ratio (%)
0.03
3
0.02
2
0.01
1
0 0
5 10 15 20 5 10 15 20
Figure 8.7 Median values of EDPs in transverse direction determined by nonlinear

RHA of single-bent overpass bridge [top-panels], multi-span bridge
[bottom-panels] due to 21 ground motions; results for individual ground
motions are also included.
102
0.4 0.25
0.2
0.3
Base Shear / Weight

Base Shear / Weight
0.15
0.2
Median
Median
0.1
0.1
0.05
0 0
0 1 2 3 4 5 0 1 2 3 4
Deck Drift Ratio (%) Deck Drift Ratio (%)
Figure 8.8 Deck drift ratios in transverse direction determined by nonlinear RHA of
single-bent overpass bridge [right] and multi-span single frame bridge [left]
for 21 ground motions identified on first-mode pushover curve.
60 120
CR Equation
50 Median
100
ith GM
40
80
D1 (cm)
D1 (cm)
30
I
I
60
20
40
10
0 20
5 10 15 20 5 10 15 20
I
ground motions for single-bent overpass bridge [right] and multi-span
I
single frame bridge [left]; “exact” target value of deformation D1 is
identified by horizontal dashed line; horizontal continuous line indicates
I
target value of deformation D1 established by CR equation.
103

5
Benchmark
MPS
4
Median
ith GM
3
0
0.04
0.03
0.02
0.01
0
1 2 3 4 5 6 71 2 3 4 5 6 71 2 3 4 5 6 7
EDPs in transverse direction for the single-bent overpass bridge; individual
results for each of the seven scaled ground motions are also presented.
104

7
Benchmark
6 MPS
Median
5
ith GM
4
1
0.05
0.04
0.03
0.02
0.01
1 2 3 4 5 6 71 2 3 4 5 6 71 2 3 4 5 6 7
EDPs in transverse direction for the multi-span bridge; individual results
for each of the seven scaled ground motions are also presented.
105
9.CONCLUSIONS
A modal-pushover-based scaling (MPS) method has been developed to scale ground motions for
use in nonlinear response history analysis (RHA) of buildings. In the MPS method, the ground
motions are scaled to match (to a specified tolerance) a target value of the inelastic deformation
of the first-mode inelastic SDF system—its properties determined by first-mode pushover
analysis—and the elastic deformation of 2nd-mode SDF systems are considered in selecting a
subset of the scaled ground motions to account for higher-mode effects.
Six building structures with a fundamental period ranging from 0.6 to 5.6 sec
representing low-, mid-, and high-rise building types in California, as well as two bridges were
utilized for testing the MPS procedure. For each structure, the median values of engineering
demand parameters (EDPs) due to three sets of 7 ground motions scaled by the MPS method
were computed by nonlinear RHA and compared against the benchmark values of EDPs,
determined by nonlinear RHA of the structure for 21 unscaled records. The ASCE 7-05
procedure was also applied to buildings in order to compare its performance in predicting the
median EDPs and reducing their dispersion with the MPS procedure. Such comparison led to the
following conclusions:
1. Even for the most intense near-fault ground motions, which represent a severe test, the
MPS method estimates the median value of seismic demands to a good degree of
accuracy (within 25% of the benchmark value). In contrast, the ASCE 7-05 scaling
method overestimates the demand by 20 to 50% for 4-, and 6-story steel buildings, and its
overestimation exceeds 50% for 13-story building. For the 19-, and 52-story buildings,
bias (overestimation or underestimation) associated with the ASCE 7-05 scaling method
reaches over 25%; for the MPS, the bias remained under 20% except for the 52-story
building for ground motion Set-3 which resulted in 25% of overestimation at certain
floors. The estimates of median value of EDPs determined by the MPS procedure are
within 90 to 112% of the benchmark values for the single-bent overpass. For the multi-
span bridge, median EDPs are overestimated as much as 18% for the ground motion Set-
106
2 and 3; the overestimation is less for the Set-1. These results demonstrate that the MPS
procedure leads to scaled ground motions that yield accurate estimate of median EDPs.
2. The dispersion of responses due to ground motion scaled by the MPS method is much
smaller compared to the ASCE 7-05 scaling method; in the latter method, dispersion is
unacceptably large for some combinations of buildings and bridges and ground motions
sets. Significant reduction of dispersion compared to the unscaled ground motions and
ASCE 7-05 procedure indicates that the MPS procedure is efficient.
3. Using the exact value of target deformation, defined as the median deformation of the
first-mode inelastic SDF system for a large ensemble of unscaled records determined by
nonlinear RHA, leads to the most accurate and efficient version of the MPS method.
Because this rigorous approach is not suitable for practical application, the target
deformation may be estimated from the deformation of the corresponding linear system,
available from the design spectrum, and empirical equations for the inelastic deformation
ratio. The increase in bias and dispersion resulting from this approximation is small. The
resulting practical version of the MPS method uses attenuation relations for elastic
spectral ordinates that are currently available; new attenuation relations for inelastic
spectral deformation are not required.
4. For first-mode dominated structures, scaling earthquake records to the target value of the
inelastic deformation is sufficient in producing accurate estimates of median EDPs and in
reducing the dispersion of EDPs due to individual ground motions. For mid-rise and
high-rise buildings where higher vibration modes are known to contribute significantly to
the seismic response, the MPS method requires an additional step to rank the scaled
ground motions based on the closeness of the elastic deformation of 2nd-mode elastic
SDF systems to their target values. Selecting a subset of highest-ranked ground motions
leads to a method that is more accurate and efficient for estimating seismic demands for
taller buildings.
This study has focused on developing the MPS method for scaling ground motions and its
evaluation based on low-, mid-, and high-rise buildings in regular plan, and two ordinary bridges.
In all models stable force deformation relations were considered.
107
10.APPENDIX-A: COMPUTATION OF
MINIMUM SCALE FACTOR
Although not explicitly stated in the ASCE 7-05 scaling procedure, it is desired to scale records
with small scale factors. To obtain a minimum scale factor for each of seven records, residuals
between the record’s scaled spectrum and target spectrum is minimized between 0.2T1 and 1.5T1
through a method of least square whereby square of sum of residuals is expressed as
n
2
l = ��
Ai - ( SF �
� Ai ) �
� (A.1)
i =1
in which Ai and Ai are respectively the target spectral acceleration and record’s (unscaled)
spectral acceleration at ith spectral period. n is the number of periods (log spaced) covered
between 0.2T1 and 1.5T1. SF is the scaling factor minimizing l, i.e., d l / dSF @ 0 ; and it is
computed from
�n ��n �
SF = ��
�i =1
( Ai �
A )
i ��/ �( Ai �
��i =1
Ai ) �
�
(A.2)
Note that Eq. (A.2) yields an optimal scaling factor to ensure that the record’s scaled
spectrum match closely the target spectrum between 0.2T1 and 1.5T1. Occasionally, the average
spectrum of seven scaled records may fall below the target spectrum in this period range and
violates the ASCE 7-05 requirement. If this situation is encountered, each record needs to be
further amplified with the same factor of the ratio of maximum difference.
108
11.REFERENCES
Abrahamson, N. A. and Silva, W. J., 1997. Empirical response spectral attenuation relations for
shallow crustal earthquakes. Seismol. Res. Lett. Vol. 68, pp. 94-127.
American Association of State Highway and Transportation Officials (2004). AASHTO - LRFD
Bridge design specifications. 3rd ed. Washington (DC).
Alavi, B., and Krawinkler, H. (2000). “Consideration of near-fault ground motion effects in
seismic design,” Proc. of the 12th World Conference on Earthquake Engineering, Paper No.
2665, Auckland, New Zealand.
Alavi, B., and Krawinkler, H. (2004). “Behavior of moment-resisting frame structures subjected
to near-fault ground motions,” Eq. Eng. and Str. Dyn., Vol. 33, No. 6, pp. 687-706.
American Society of Civil Engineers (2005), ASCE 7-05 Minimum Design Loads for Buildings,
Reston, VA.
Anderson, J., and Bertero, V. V. (1991). “Seismic performance of an instrumented six-story steel
building,” Earthquake Engineering Research Center, University of California, Berkeley,
Report No. 1991-11, 134 p.
Báez, J. I., and Miranda, E. (2000). “Amplification factors to estimate inelastic displacement
demands for the design of structures in the near field,” Proc. of the Twelfth World Conf. on
Eq. Eng., Paper No. 1561, Auckland, New Zealand.
Baker, J. W., and Cornell, A. C. (2006). “Spectral shape, epsilon and record selection,”
Earthquake Engineering & Structural Dynamics, Vol. 35, No. 9, pp. 1077-1095.
Bazzurro, P., 1998. Probabilistic Seismic Demand Analysis, Ph.D. thesis, Dept. of Civil and Env.
Eng., Stanford University, CA. (Available online at
http://www.stanford.edu/group/rms/Thesis/index.html; last accessed on 06/2008).
Bazzurro, P., and Luco, N. (2004). “Parameterization of Non-Stationary Acceleration Time
Histories. Lifelines Program Project 1G00 Addenda. Pacific Earthquake Engineering
Research (PEER) Center, University of California: Berkeley, CA, 83 p.
Bazzurro, P., and Luco, N. (2006). “Do scaled and spectrum-matched near-source records
produce biased nonlinear structural responses,” Proc. of the 8th National Conf. on Eq. Eng.,
Paper No. 1029, San Francisco, CA.
109
Bolt B. A. and Gregor N. J. Synthesized strong ground motions for the seismic condition
assessment of the eastern portion of the San Francisco bay bridge. Report UCB/EERC-93/12,
University of California, Earthquake Engineering Research Center, Berkeley, CA.
Bozorgnia, Y., and Mahin, S. A. (1998). “Ductility and strength demands of near-fault ground
structure-specific scalar intensity measures 389 motions of the Northridge earthquake,”
Proceedings of the 6th U.S. National Conf. on Eq. Eng., Seattle, WA.
Carballo J. E. and Cornell C. A. Probabilistic seismic demand analysis: Spectrum matching and
design. Department of Civil and Environmental Engineering, Stanford University 2000,
Report NO. RMS-41.
Chopra, A.K. (2001). Dynamics of Structures: Theory and Applications to Eq. Eng., 2nd Ed.,
Prentice Hall, Englewood Cliffs, N.J.
Chopra, A. K., and Chinatanapakdee, C. (2004). “Inelastic Deformation Ratios for Design and
Evaluation of Structures: Single-Degree-of-Freedom Bilinear Systems,” J. of Str. Eng.
(ASCE), Vol. 130, No. 9, pp. 1304-1319.
Cordova, P. P., Deierlein, G. G., Mehanny, S. S. F., and Cornell, C. A., 2000. Development of a
two-parameter seismic intensity measure and probabilistic assessment procedure,
Proceedings of the 2nd U.S.-Japan Workshop on Performance-Based Seismic Design
Methodology for Reinforced Concrete Building Structures, PEER Report 2000/10, Pacific
Earthquake Engineering Research Center, University of California, Berkeley, CA.
Graizer V. Kalkan E. 2009. Prediction of Response Spectral Acceleration Ordinates based on
PGA Attenuation, Earthquake Spectra, Vol. 25, No. 1, pp. 36 – 69.
International Conference of Building Officials (2006). International Building Code, Whittier,
CA.
International Conference of Building Officials (2007). California Building Code, Whittier, CA.
Krawinkler, H., and Al-Ali, A. (1996). “Seismic demand evaluation for a 4-story steel frame
structure damaged in the Northridge earthquake,” The Structural Design of Tall Buildings,
Vol. 5, No. 1, pp. 1-27.
Kennedy, R. P., Short, S. A., Merz, K. L., Tokarz, F. J., Idriss, I. M., Power, M. S., and Sadigh, K.
(1984). “Engineering characterization of ground motion-task 1: Effects of characteristics of
free-field motion on structural response,” NUREG/CR-3805, U.S. Regulatory Commission,
Washington, D.C.
110
Kurama, Y., and Farrow, K. (2003). “Ground motion scaling methods for different site conditions
and structure characteristics,” Eq. Eng. and Str. Dyn., Vol. 32, No. 15, pp. 2425-2450.
Lilhhand, K., and Tseng, W. S. (1989). “Development and application of realistic earthquake
time histories compatible with multiple damping design spectra,” Proc. of the 9th World
Conf. on Eq. Eng. Tokyo-Kyoto, Japan, Vol. 2, pp. 819-830.
Luco, N., and Cornell, A. C. (2007). “Structure-Specific Scalar Intensity Measures for Near-
Source and Ordinary Earthquake Ground Motions,” Earthquake Spectra, Vol. 23, No. 2, pp.
357–392.
Malhotra, P. K. (2003). “Strong-Motion Records for Site-Specific Analysis,” Earthquake Spectra,
Vol. 19, No. 3, pp. 557-578.
Mehanny, S. S. F. (1999). “Modeling and Assessment of Seismic Performance of Composite
Frames with Reinforced Concrete Columns and Steel Beams,” Ph.D. thesis, Dept. of Civil
and Env. Eng., Stanford University, California.
Miranda, E. (1993). “Evaluation of site-dependent inelastic seismic design spectra,” J. of Str.
Eng. (ASCE), Vol. 119, No. 5, pp. 1319–1338.
Naeim, F., Alimoradi, A., and Pezeshk, S. (2004). “Selection and scaling of ground motion time
histories for structural design using genetic algorithms,” Earthquake Spectra, Vol. 20, No. 2,
pp. 413-426.
Nau, J., and Hall, W. (1984). “Scaling methods for earthquake response spectra,” J. of Str. Eng.
(ASCE), Vol. 110, No. 91-109.
Neuenhofer, A., and Filippou, F.C. (1998). “Geometrically nonlinear flexibility-based frame
finite element,” J. of Str. Eng. (ASCE), Vol. 124, No. 6, pp. 704-711.
OpenSEES. Open Source finite element platform for earthquake engineering simulations 2006.
University of California Berkeley, Pacific Earthquake Engineering Center. (Available online
at http://OpenSEES.berkeley.edu/; last accessed on 06/2008).
Petersen, M. D., Frankel, A. D., Harmsen, S. C., Mueller, C. S., Haller, K. M., Wheeler, R. L.,
Wesson, R. L., Zeng, Y., Boyd, O. S., Perkins, D. M., Luco, N., Field, E. H., Wills, C. J., and
Rukstales, K. S., (2008). “Documentation for the 2008 Update of the United States National
Seismic Hazard Maps,” U.S. Geological Survey Open-File Report 2008–1128, 61 p.
(Available online at: http://pubs.usgs.gov/of/2008/1128/; last accessed on 06/2008).
111
Power, M., Chiou, B., Abrahamson, N., and Roblee, C. 2006. “The Next Generation of Ground
Motion Attenuation Models” (NGA) project: An overview,” In Proc. of the Eighth National
Conf. on Eq. Eng., Paper No. 2022, San Francisco, CA.
Ruiz-Garcia, J., and Miranda, E. (2002). “Inelastic displacement ratios for evaluation of existing
structures,” Eq. Eng. and Str. Dyn., Vol. 31, No. 3. pp. 539-560.
Seismic Design Criteria (2006), available online at
http://www.dot.ca.gov/hq/esc/techpubs/manual/othermanual/other-engin-manual/seismic-
design-criteria/sdc.html
Silva W. and Lee K. State-of-the-art for assessing earthquake hazards in the United States;
Report 24, WES RASCAL Code for Synthesizing Earthquake Ground Motions 1987.
Miscellaneous Paper S-73-1, U.S. Army Engineer Waterways Experiment Station, Vicksburg,
MS.
Shome, N., and Cornell, A. C. (1998). “Normalization and scaling accelerograms for nonlinear
structural analysis,” Prof. of the 6th U.S. National Conf. on Earthquake Engineering, Seattle,
WA.
Shome, N., Cornell, C. A., Bazzurro, P., and Carballo, J. E., 1998. Earthquakes, records, and
nonlinear responses, Earthquake Spectra, Vol. 14, No.3, pp. 469–500.
Shome, N., and Cornell, C. A., 1999. Probabilistic Seismic Demand Analysis of Nonlinear
Structures, Reliability of Marine Structures Program Report No. RMS-35, Dept. of Civil and
Env. Eng., Stanford University, CA. (Available online at
http://www.stanford.edu/group/rms/Thesis/index.html; last accessed on 06/2008).
Tothong, P., and Cornell, A. C. (2008). “Structural performance assessment under near-source
pulse-like ground motions using advanced ground motion intensity measures,” Eq. Eng. and
Str. Dyn. Vol. 37, No. 7, pp. 1013-1037.
Veletsos, A. S., and Newmark, N. M. (1960). “Effect of inelastic behavior on the response of
simple systems to earthquake motions,” Proc. of the Second World Conf. on Eq. Eng., Tokyo,
Japan, Vol. II, pp. 895–912.
Vidic, T. Fajfar, P., and Fischinger, M. (1994). “Consistent inelastic design spectra: Strength and
displacement,” Eq. Eng. and Str. Dyn., Vol. 23, No. 5, pp. 507-521.
Youngs, R. Power, M., Wang, G., Makdisi, F., and Chin, C. C. (2007). “Design Ground Motion
Library (DGML) – Tool for Selecting Time History Records for Specific Engineering
112
Applications (Abstract),” SMIP07 Seminar on Utilization of Strong-Motion Data, pp. 109 –

110. (Available online at http://www.conservation.ca.gov/cgs/smip/docs/seminar
/SMIP07/Pages/Paper8_Youngs.aspx; last accessed on 06/2008).
NOTATION
The following symbols are used in this report:
a = Parameter in empirical formula for CR

An = Target pseudo-spectral acceleration
b = Parameter in empirical formula for CR
c = Parameter in empirical formula for CR
Ratio of peak deformations of inelastic and corresponding elastic SDF
CR = systems for systems with known yield-strength reduction factor
d = Parameter in empirical formula for CR
Dn = Target value of nth mode elastic deformation
D1I = First-mode target value of inelastic spectral displacement
Dn(t) = Deformation of a SDF system
D1I = Peak deformation of inelastic SDF system
Dn = Peak deformation of elastic SDF system
D1,y = Yield deformation of inelastic SDF system
Fs1 = System resisting force under first-mode pushover
LR = Deformation ratio CR for zero-period system
m = Mass matrix of a MDF system
M = Moment magnitude of earthquake
M* Effective modal mass
n = Mode sequence number
Ry = Yield-strength reduction factor
Rcl = Closest distance to co-seismic rupture plane
sn* = Load vector of modal pushover analysis
SF = Ground motion scaling factor
Ta = Period defined in Newmark-Hall smooth design spectrum (Fig. 4.6)
Tb = Period defined in Newmark-Hall smooth design spectrum (Fig. 4.6)
113
Tc = Period separating the acceleration- and velocity-sensitive regions

Td = Period separating the velocity- and displacement-sensitive regions
Te = Period defined in Newmark-Hall smooth design spectrum (Fig. 4.6)
Tf = Period defined in Newmark-Hall smooth design spectrum (Fig. 4.6)
Tn = Elastic natural vibration period
ur1 = Roof displacement of a MDF system under first-mode pushover
u&
&g = Earthquake ground acceleration
Vb1 = Base shear under first-mode pushover
Vb1y = Global yield strength under first-mode pushover
VS30 = Average shear-wave velocity within 30m depth from surface
a = Ratio of post-yield and initial stiffness
 = Damping ratio
 = Modal participation factor
f = Mode shape
i = Influence vector
114

E Kalkan

Uploaded by

Copyright:

Available Formats

E Kalkan

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

E Kalkan

Uploaded by

Copyright:

Available Formats

Practical Guidelines to Select and

Scale Earthquake Records for

By Erol Kalkan and Anil K. Chopra

Open-File Report 2009–####

U.S. Geological Survey

By Erol Kalkan and Anil K. Chopra

Open-File Report 2009–####

U.S. Department of the Interior

U.S. Geological Survey

2.2. ESTIMATING SDF-SYSTEM INELASTIC DEFORMATION........................................13

3. CODE-BASED SCALING PROCEDURE................................................................15

5. EVALUATION OF MPS PROCEDURE: LOW- AND MID-RISE BUILDINGS.........24

5.2. SYSTEMS ANALYZED...............................................................................................26

5.3 FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................27

5.4. EVALUATION OF MPS CONCEPT.............................................................................27

5.5. EVALUATION OF MPS AND CODE-BASED SCALING PROCEDURES...................28

5.6. MULTI-MODE CONSIDERATIONS.............................................................................30

6. EVALUATION OF MPS PROCEDURE: HIGH-RISE BUILDINGS..........................52

6.2 ANALYTICAL MODELS AND CALIBRATION STUDIES............................................54

6.3. FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................55

6.4. EVALUATION OF MPS PROCEDURE.......................................................................56

6.5. COMPARING MPS AND CODE-BASED SCALING PROCEDURES.........................57

6.6. MULTI-MODE CONSIDERATIONS.............................................................................58

7. EVALUATION OF MPS PROCEDURE: SHORT - PERIOD BUILDING.................80

7.2. FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................81

7.3. EVALUATION OF MPS PROCEDURE.......................................................................81

8. EVALUATION OF MPS PROCEDURE: BRIDGES.................................................90

8.2 FIRST-MODE SDF-SYSTEM PARAMETERS.............................................................93

8.3. EVALUATION OF MPS PROCEDURE.......................................................................94

1.1. RESEARCH OBJECTIVE

2.1.MPS PROCEDURE: SUMMARY

The MPS procedure is summarized below in a step-by-step form:

pseudo-acceleration spectrum for a large ensemble of (unscaled) earthquake records

vibration of the structure.

First-mode Dominated Structures

value of base shear determined from the idealized pushover curve.

the influence vector i is equal to unity.

from D1 = CR D1 , where D1 = ( T1 2p ) A1 and A1 is the target pseudo-spectral

acceleration at period T1 and C R is determined from an empirical equation (shown in the

reduction factor R y , determined in Step 4.

selected ground motions u&

systems, from the target spectrum Dn = ( Tn 2p ) An , where the mode number n = 2. We

scale factor SF determined in Step 9.

and the peak deformation determined in Step 11.

2.2.ESTIMATING SDF-SYSTEM INELASTIC DEFORMATION

3.CODE-BASED SCALING PROCEDURE

4.GROUND MOTIONS SELECTION

- Magnitude range of anticipated significant event(s)

- Distance range of the site from the causative fault(s)

- Site-condition (site-geology generally described by average shear-wave velocity

- Basin effect (if basin exists)

Soil Condition Dependence

4.1 GROUND MOTION ENSEMBLE

Table 4.1 Selected earthquake ground motions

Figure 4.1 Comparison of average spectral shape of earthquakes in magnitude range

5.EVALUATION OF MPS PROCEDURE:

5.1 SELECTED BUILDINGS

5.1.1 Four-story Building