Soil Structure Interaction For Building
Soil Structure Interaction For Building
Soil Structure Interaction For Building
Soil-Structure
Interaction for
Building Structures
This report was produced by the NEHRP Consultants Joint Venture, a joint venture of the Applied Technology Council (ATC)
and the Consortium of Universities for Research in Earthquake Engineering (CUREE). While endeavoring to provide practical
and accurate information, the NEHRP Consultants Joint Venture, the authors, and the reviewers assume no liability for, nor
express or imply any warranty with regard to, the information contained herein. Users of information contained in this report
assume all liability arising from such use.
Unless otherwise noted, photos, figures, and data presented in this report have been developed or provided by NEHRP
Consultants Joint Venture staff or consultants engaged under contract to provide information as works for hire. Any similarity
with other published information is coincidental. Photos and figures cited from outside sources have been reproduced in this
report with permission. Any other use requires additional permission from the copyright holders.
Certain commercial software, equipment, instruments, or materials may have been used in the preparation of information
contributing to this report. Identification in this report is not intended to imply recommendation or endorsement by NIST, nor is
it intended to imply that such software, equipment, instruments, or materials are necessarily the best available for the purpose.
NIST policy is to use the International System of Units (metric units) in all its publications. In this report, information is
presented in both metric units and U.S. Customary Units (inch-pound), as the inch-pound system is the preferred system of
units in the U.S. earthquake engineering industry.
Cover illustration – Rendering of the Sherman Oaks building structural and foundation system, used in an example application
of soil-structure interaction principles (courtesy of C. Haselton).
NIST GCR 12-917-21
Prepared for
U.S. Department of Commerce
National Institute of Standards and Technology
Engineering Laboratory
Gaithersburg, MD 20899
By
NEHRP Consultants Joint Venture
A partnership of the Applied Technology Council and the
Consortium of Universities for Research in Earthquake Engineering
September 2012
Participants
National Institute of Standards and Technology
John (Jack) R. Hayes, Jr., Director, National Earthquake Hazards Reduction Program
Steven L. McCabe, Deputy Director, National Earthquake Hazards Reduction Program
John (Jay) L. Harris III, Project Manager
This work is the result of an extensive literature search and collection of available
information on soil-structure interaction, discussions with researchers and
practitioners on the state of SSI knowledge and practice, conduct of problem-focused
investigations, and analytical parametric studies. A workshop of invited experts and
other stakeholders was convened to receive feedback on the developing report and
preliminary recommendations. Input from this group was instrumental in shaping the
final product.
Jon A. Heintz
Program Manager
=0.33: (a) geometry; (b) x-direction; and (c) y-direction .................. 2-15
of a homogeneous half-space, with zero hysteretic damping, and
Figure 2-5 Plot of dynamic stiffness modifiers and damping ratios versus
dimensionless frequency, for square footings embedded in a
=0.33: (a) geometry; and (b) x-direction (y-direction similar) .......... 2-16
homogeneous half-space, with zero hysteretic damping, and
Figure 2-6 Plot of dynamic stiffness modifiers and radiation damping ratios versus
dimensionless frequency comparing uniform half-space and non-uniform
profiles in which G varies with depth: (a) translation for circular
foundations; and (b) rotation for square foundations ........................... 2-19
with =0.4 and p/s =1.3: (a) geometric parameters; (b) lateral
dimensionless frequency, for single piles in a homogeneous half-space,
Table 2-1 Values of Shear Wave Velocity and Shear Modulus Reduction for
Various Site Classes and Shaking Amplitudes .................................... 2-10
Table 2-2a Elastic Solutions for Static Stiffness of Rigid Footings at the
Ground Surface .................................................................................... 2-11
Table 2-2b Embedment Correction Factors for Static Stiffness of Rigid
Footings ............................................................................................... 2-12
Table 2-3a Dynamic Stiffness Modifiers and Radiation Damping Ratios for
Rigid Footings...................................................................................... 2-13
Table 2-3b Dynamic Stiffness Modifiers and Radiation Damping Ratios for
Embedded Footings ............................................................................. 2-14
Table 2-4a Equations for Static Stiffness of Single Piles ...................................... 2-29
Table 2-4b Equations for Dynamic Stiffness Modifiers and Damping Ratios for
Single Piles .......................................................................................... 2-30
Table 5-1 Summary of Field Forced-Vibration Tests Used to Evaluate
Foundation Impedance Functions .......................................................... 5-2
Table 5-2 Summary of Laboratory-Scale Tests of Dynamic Response of
Footings ................................................................................................. 5-7
Table 7-1 Summary of Modeling Approaches Considered in Example
Applications ........................................................................................... 7-6
Table 7-2 Summary of Earthquake Events Recorded at the Sherman Oaks
Building ............................................................................................... 7-10
Table 7-3 Summary of Effective Profile Depths and Average Effective Profile
Velocities for the Sherman Oaks Building .......................................... 7-13
Table 7-4 Calculation of Shallow Foundation Stiffness and Damping
Parameters for the Sherman Oaks Building ......................................... 7-15
Table 7-5 Calculation of Pile Stiffness and Damping Parameters for the
Sherman Oaks Building ....................................................................... 7-17
Table 7-6 Comparison of Fundamental Periods for Alternative Foundation
Modeling Configurations for the Sherman Oaks Building .................. 7-20
Table 7-7 Summary of Sherman Oaks Building Parametric Stick Models .......... 7-28
Table 7-8 Comparison of Results for Moment Frame, Core Shear Wall, and
Perimeter Shear Wall Stick Models and Alternative Foundation
Modeling Configurations ..................................................................... 7-30
This report presents a synthesis of the body of knowledge contained in SSI literature,
which has been distilled into a concise narrative and harmonized under a consistent
set of variables and units. Specific techniques are described by which SSI
phenomena can be simulated in engineering practice, and recommendations for
modeling seismic soil-structure interaction effects on building structures are
provided.
Problems associated with the practical application of SSI for building structures are
rooted in a poor understanding of fundamental SSI principles. Soil-structure
interaction topics are generally not taught in graduate earthquake engineering
courses, so most engineers attempting SSI in practice must learn the subject on their
own. Unfortunately, practice is hindered by a literature that is often difficult to
understand, and codes and standards that contain limited guidance. Most articles rely
heavily on the use of wave equations in several dimensions and complex arithmetic
to formulate solutions and express results. Moreover, nomenclature is often
inconsistent, and practical examples of SSI applications are sparse. This gives rise to
the present situation in which soil-structure interaction is seldom applied, and when it
is, modeling protocols vary widely and are not always well conceived.
Although the NEHRP Recommended Provisions for Seismic Regulations for New
Buildings and Other Structures has provided guidance for consideration of SSI
effects in forced-based procedures for several decades, these procedures have not
found significant use in practice. Practical application of SSI gained momentum
following publication of FEMA 440, Improvement of Inelastic Seismic Analysis
Procedures (FEMA, 2005), which provided the design community with procedures
for incorporating the effects of soil-structure interaction in nonlinear static pushover-
type analyses. These procedures were eventually adopted into ASCE/SEI 41-06,
In 2009, the National Institute of Standards and Technology (NIST) initiated the Task
Order 69221 Project entitled “Improved Procedures for Characterizing and Modeling
Soil-Structure Interaction for Performance-Based Seismic Engineering.” The
purpose of this project was to develop consensus guidance for implementing soil-
structure interaction in response history analyses, such that input ground motions
accurately reflect the input at the base of structures, and that structural models
include elements that account for the geotechnical and foundation conditions
associated with the building under consideration. Work also included an extensive
review of available research on soil-structure interaction, evaluation of existing SSI
guidelines for static-type analyses, and development of recommendations for
improvement where necessary.
1. When is the use of foundation springs and dashpots important, and which
structural response parameters are affected?
Once the decision to implement SSI has been made, a basic level of understanding of
the physical phenomenon and a practical analysis methodology for simulating its
effects are needed. This report describes the principal components of SSI in a clear
and concise way, and consistent nomenclature is used throughout. Explicit
computational tools that can be used in engineering practice are provided, and
applications of SSI to force-based analysis procedures, pushover (displacement-
based) procedures, and response history analysis procedures are described.
Methods that can be used to evaluate the above effects can be categorized as direct
and substructure approaches. In a direct analysis, the soil and structure are included
within the same model and analyzed as a complete system. In a substructure
approach, the SSI problem is partitioned into distinct parts that are combined to
formulate the complete solution.
Evaluation of site response using wave propagation analysis through the soil is
important to this approach. Such analyses are most often performed using an
equivalent linear representation of soil properties in finite element, finite difference,
Chapter 2 describes inertial SSI effects, beginning with the behavior of simple
structure-soil systems to provide insight into the conditions for which inertial SSI
effects are most critical. Also presented are detailed procedures for computing
foundation stiffness and damping, both for idealized conditions in classical solutions
and for more realistic conditions that may include flexible foundation systems, non-
uniform soil, and material nonlinearity. Both shallow foundation systems (e.g.,
footings, mats) and deep foundation systems (e.g., piles) are discussed.
Chapter 3 describes the manner by which transfer functions can be computed for
various foundation configurations considering kinematic interaction effects, and how
transfer function amplitudes can be used to modify response spectral ordinates.
Chapter 4 describes how SSI procedures described in Chapter 2 and Chapter 3 are
implemented in seismic design provisions contained in currently available consensus
standards and design guidelines.
Chapter 5 describes the relatively limited number of SSI experiments and case
history studies available in the literature, and summarizes some of the lessons learned
from that work. Efforts to calibrate or verify analysis procedures against laboratory
data or field performance data are noted, where applicable.
Chapter 8 answers three fundamental questions that arise when contemplating the use
of SSI in the seismic analysis of a building structure, summarizes SSI analysis
procedures that would be followed on typical project in a step-by-step manner, and
provides short-term and long-term research needs.
A list of Symbols defining key notation, and a list of References cited from the body
of SSI literature are provided at the end of this report.
Section 2.1 discusses system behavior and highlights some of the principal effects of
inertial interaction and the conditions for which its effects are significant. The
methods focus on single degree-of-freedom systems, but they can be extrapolated to
multi-degree-of-freedom systems with a dominant first mode. Section 2.2 provides a
relatively detailed description of how foundation springs and dashpots can be
specified to represent the flexibility and damping associated with soil-foundation
interaction in translational and rotational vibration modes for shallow foundations
(e.g., footings and mats). Section 2.3 provides corresponding solutions for the
stiffness and damping characteristics of deep foundations (e.g., piles and drilled
shafts). Some of the procedures given in Section 2.2 are coded into available
computer programs such as DYNA6 (Western Engineering, 2011). This program can
also be used for pile foundations, although the results are relatively approximate.
Section 2.4 presents several models that can be used to evaluate shallow foundation
response for conditions involving nonlinear material behavior or geometric
nonlinearities (i.e. gapping).
A rigid base refers to soil supports with infinite stiffness (i.e., without soil springs).
A rigid foundation refers to foundation elements with infinite stiffness (i.e., not
deformable). A fixed base refers to a combination of a rigid foundation elements on
a rigid base. A flexible base analysis considers the compliance (i.e., deformability) of
both the foundation elements and the soil.
F
(2-1)
k
From structural dynamics, the undamped natural vibration frequency, , and period,
T, of the structure are given by Clough and Penzien (1993) as:
By substituting Equation 2-1 into Equation 2-2, an expression for the square of
period is obtained as:
2 m
T 2 2 2
F
2 m
(2-3)
F
(a) (b)
Figure 2-1 Schematic illustration of deflections caused by force applied to: (a)
fixed-base structure; and (b) structure with vertical, horizontal, and
rotational flexibility at its base.
Now consider the same structure with vertical, horizontal, and rotational springs at its
base, representing the effects of soil flexibility against a rigid foundation, as depicted
in Figure 2-1b. The vertical spring stiffness in the z direction is denoted kz, the
horizontal spring stiffness in the x direction is denoted kx, and the rotational spring is
denoted kyy, representing rotation in the x-z plane (about the y-y axis). If a force, F, is
applied to the mass in the x direction, the structure deflects, as it does in the fixed-
u f h
F
k
F F F h
(2-4)
h
k k x k yy
If Equation 2-4 is substituted into Equation 2-3, an expression for flexible base
period, T , is obtained as:
(2-5)
F
T 1 1 h2
m
2
k
T k k x k yy
(2-6)
m
Equation 2-6 simplifies into a classical period lengthening expression (Veletsos and
Meek, 1974):
1
T k kh 2
(2-7)
T k x k yy
In previous work by Veletsos and Nair (1975) and Bielak (1975), it has been shown
that the dimensionless parameters controlling period lengthening are:
, and
s 4 BLh
h h B m
, , , (2-8)
V sT B L
where h is the structure height (or height to the center of mass of the first mode
mass (or effective modal mass),s is the soil mass density, and is the Poisson’s
shape), B and L refer to the half-width and half-length of the foundation,m is the
ratio of the soil. Previous work was applicable to circular foundations, and has been
adapted here for rectangular shapes considering the ratio, B/L.
To the extent that h/T quantifies the stiffness of the superstructure, the term h/(VsT) in
Equation 2-8 represents the structure-to-soil stiffness ratio. The term h/T has units of
velocity, and will be larger for stiff lateral force resisting systems, such as shear
walls, and smaller for flexible systems, such as moment frames. The shear wave
velocity, Vs, is closely related to soil shear modulus, G, computed as:
Vs G / s (2-9)
For typical building structures on soil and weathered rock sites, h/(VsT) is less than
0.1 for moment frame structures, and between approximately 0.1 and 0.5 for shear
wall and braced frame structures (Stewart et al., 1999b). Period lengthening
structure system. The mass ratio, m/s4BLh, is the ratio of structure mass to the mass
ratio, B/L, in Equation 2-8 are aspect ratios describing the geometry of the soil-
of soil in a volume extending to a depth equal to the structure height, h, below the
foundation. In Equation 2-7, it can be seen that period lengthening has no
fundamental dependence on mass. The mass ratio term was introduced so that period
lengthening could be related to easily recognizable characteristics such as structural
first mode period, T, and soil shear wave velocity, Vs, rather than structural stiffness,
taken as 0.15 (Veletsos and Meek, 1974). The Poisson’s ratio of the soil, , affects
k, and soil shear modulus, G. The effect of mass ratio is modest, and it is commonly
Using models for the stiffness of rectangular foundations (of half-width, B; half-
length, L; and L ≥ B) resting on a homogeneous isotropic half-space with shear wave
velocity, Vs, period lengthening ratios can be calculated with the results shown in
Figure 2-2a, which is plotted for the special case of a square footing (L = B).
and varying ratios of h/B. In this plot, = 0.33, B/L = 1.0, hysteretic
versus structure-to-soil-stiffness ratio for square foundations (L = B)
rotation, . This implies that inertial SSI effects would be more significant in tall
to-foundation-width ratio, h/B, due to increased overturning moment and foundation
buildings, but this is not the case. Tall buildings typically have low h/(VsT) ratios,
which is more important for controlling inertial SSI effects. Hence period
lengthening in tall buildings is near unity (i.e., little or no period lengthening). For a
fixed ratio of h/B, period lengthening is observed to decrease modestly with
foundation-width-to-length ratio, B/L, due to increased foundation size (and therefore
stiffness) normal to the direction of loading.
This damping is composed of two parts: (1) contributions from soil hysteresis
(hysteretic damping); and (2) radiation of energy away, in the form of stress waves,
0 f i
T T
1
n
(2-10)
where i is the structural damping in the superstructure assuming a fixed base, which
is generally taken as 5% for typical structural systems. More refined estimates of i
are possible based on structural system type and configuration, as described in
PEER/ATC-72-1, Modeling and Acceptance Criteria for Seismic Design and
al., 1999b) have shown that f ranges from approximately 0% to 25%. The exponent,
Analysis of Tall Buildings (ATC, 2010). Observations from case studies (Stewart et
n, on the period lengthening term in Equation 2-10 is taken as 3 for linearly viscous
structural damping, and 2 otherwise (e.g., for hysteretic damping) (Givens, 2013).
Analytical models for foundation damping have been presented by Veletsos and Nair
(1975), Bielak (1975 and 1976), Roesset (1980), Wolf (1985), Aviles and Perez-
Rocha (1996), Maravas et al. (2007), and Givens (2013), among others. The classical
solution of Veletsos and Nair accounts for the frequency dependence of foundation
and others is complex-valued (i.e., composed of real plus imaginary values), which
complicates the interpretation of its physical meaning. Bielak’s work utilizes the
same conditions except that the foundation is assumed to be a cylinder penetrating a
half-space to an embedment depth, D, and the resulting expressions are real-valued.
The value of exponent n in Equation 2-10 is taken as 3 for the Veletsos and Bielak
solutions because structural damping is assumed to be viscous.
T T n s 1
be re-written as:
f s x yy
T T T T n yy
1 1
(2-11a)
ns nx
T x T yy
Where s is soil hysteretic damping, x and yy are damping ratios related to radiation
damping from translational and rotational modes (described further in Section 2.2),
and Tx and Tyy are fictitious vibration periods, calculated as if the only source of the
vibration was foundation translation or rotation, as follows:
T x 2 T yy 2
m mh 2
(2-11b)
kx k yy
Exponents ns, nx, and nyy depend on the specific form of damping associated with the
respective components of the foundation damping, and all other terms are as
previously defined. However, because none of these terms would be expected to be
Note that for n = ns, the period lengthening terms in front of the i term in Equation
linearly viscous, it is recommended to take these exponents as 2 (Givens, 2013).
2-10 and the s term in Equation 2-11a are weight factors that together sum to unity.
Accordingly, Equation 2-11a can be viewed as a “mixing rule” for damping in
different vibration modes and sources. Because Wolf’s results were produced
neglecting the frequency dependence of foundation stiffness terms, Equation 2-11a
can provide more accurate results if those effects are included in the period
lengthening calculation.
Soil hysteretic damping, s, is strain-dependent, and can typically be evaluated from
information in the literature. Classical models are summarized in Kramer (1996).
More contemporary empirical models by Darendeli (2001) and Menq (2003) account
for overburden pressure and shear strain in a consistent manner across multiple soil
types.
The Wolf solution for foundation damping in Equation 2-11a, along with the classical
Veletsos, Bielak, and Roesset solutions, neglect contributions from terms involving
the product of two damping ratios. Maravas et al. (2007) presents exact solutions in
which those terms are included. Like Wolf, Maravas et al. (2007) utilizes hysteretic
damping so exponents n = 2, and if terms involving the product of two damping
ratios are excluded, Equation 2-11a is recovered.
h/(VsT). In Figure 2-2b, all exponents were taken as 2. Damping f decreases with
increasing values of h/B, indicating that lateral movements of the foundation (which
rocking (which dominates at high h/B ratios). Radiation damping terms (x and yy)
dominate at low h/B ratios) dissipate energy into soil more efficiently than foundation
Analysis procedures for T / T and f similar to those described above have been
validated relative to observations from instrumented buildings shaken by earthquakes
(Stewart et al., 1999a; 1999b). These studies show that the single most important
parameter controlling the significance of inertial interaction is h/(VsT), and that
inertial SSI effects are generally negligible for h/(VsT) < 0.1, which occurs in flexible
structures (e.g., moment frame buildings) located on competent soil or rock.
Conversely, inertial SSI effects tend to be significant for stiff structures, such as shear
wall or braced frame buildings, located on softer soils.
The effect of inertial SSI on the base shear of a building is illustrated in Figure 2-3.
Because base shear for elastic response is commonly computed based on pseudo-
spectral acceleration in the first mode, the figure depicts the variation in pseudo-
spectral acceleration versus period in both linear and log scales. The pseudo-spectral
acceleration for a flexible-base structure, Sa , is obtained by entering the spectrum
drawn for effective damping ratio, 0, at the corresponding elongated period, T .
Figure 2-3 Illustration of inertial SSI effects on spectral acceleration (base shear)
associated with period lengthening and change in damping.
The effect of SSI on base shear is related to the slope of the spectrum. Base shear
tends to increase when the slope is positive and decrease when the slope is negative.
The period at which the spectral peak occurs, referred to as the predominant period of
ground motion, Tp, is generally controlled by the tectonic regime, earthquake
magnitude, and site-source distance (Rathje et al., 2004), and will only match the site
period in certain cases involving large impedance contrasts due to soil layering. In
the absence of unusual site effects, typical values of Tp range from approximately 0.2
to 0.5 seconds for shallow crustal earthquakes in tectonically active regions, such as
California.
for Equation 2-12a is:
k j k j 1 2i j (2-12b)
where:
c j
j (defined for kj > 0) (2-13a)
2k j
harmonic excitation and response at a given frequency. The phase difference, j,
The imaginary part of the complex impedance represents a phase difference between
between force and (lagged) displacement is (Clough and Penzien 1993; Wolf 1985):
Angle j is also known as a loss angle. For example, if j is 10%, peak harmonic
displacement will lag peak force by 0.197 radians (11.3 degrees). When j goes to
infinity, j is bounded by /2.
Many impedance function solutions are available for rigid circular or rectangular
foundations located on the surface of, or embedded within, a uniform, elastic, or
visco-elastic half-space. In the case of a rigid rectangular foundation resting on the
surface of a half-space with shear wave velocity Vs, Pais and Kausel (1988), Gazetas
(1991), and Mylonakis et al. (2006) review impedance solutions in the literature and
present equations for computing the stiffness and damping terms in Equation 2-12.
Solutions describe translational stiffness and damping along axes x, y, and z, and
rotational stiffness and damping about those axes (denoted xx, yy, and zz). Stiffness
Poisson’s ratio of the soil, , dynamic stiffness modifiers, j, and embedment
is denoted kj, and is a function of foundation dimensions, soil shear modulus, G,
modifiers, j:
k j K j j j (2-14a)
K j GB m f B L , , j f B L , a 0 (2-14b)
j f B L , D B , d w B , Aw B L (2-14c)
where Kj is the static foundation stiffness at zero frequency for mode j, and m = 1 for
translation, and m = 3 for rotation. Shear modulus, G, should reflect the effects of
modulus reduction with increasing shear strain amplitude. ASCE/SEI 7-10,
Minimum Design Loads for Buildings and Other Structures (ASCE, 2010), and
FEMA P-750, NEHRP Recommended Seismic Provisions for New Buildings and
Other Structures (FEMA, 2009), provide the information presented in Table 2-1 for
adjusting the shear modulus and shear wave velocity for large strain levels.
the soil mass density. An average effective value of Vs is generally computed across
B
a0 (2-15)
Vs
essentially unique for half-space conditions, but may not be so in presence of a stiff
spring and dashpot coefficients that depend on a0. This can be taken as the frequency
corresponding to the period associated with the dominant response of the structure.
In most cases, this will be the first-mode, flexible-base period.
Notes: (1) SDS is the short period spectral response acceleration parameter defined in
ASCE/SEI 7-10; use straight line interpolation for intermediate values of SDS/2.5.
(2) Value should be evaluated from site-specific analysis.
Table 2-2a lists expressions for static foundation stiffness, Kj, for three translational
and three rotational degrees of freedom for a rigid rectangular footing located at the
ground surface. These equations are similar for Pais and Kausel (1988), Gazetas
(1991), and Mylonakis et al. (2006).
stiffness. Factors, j, to increase Kj for the effects of embedment are provided in
Embedment of foundations below the ground surface increases static foundation
Table 2-2b. The Pais and Kausel (1988) equations are most often used in practice.
The equations by Gazetas (1991) and Mylonakis et al. (2006) are more general,
accounting for embedment effects resulting from gapping between the soil and
foundation side walls.
Equations for dynamic stiffness modifiers, j, and radiation damping ratios, j, for
rigid footings located at the ground surface are provided in Table 2-3a. Dynamic
stiffness modifiers and radiation damping ratios for embedded footings are provided
in Table 2-3b.
1 1
GB
B L
Translation along z-axis
2 2
GB 2GL
B B L
Translation along y-axis 2 2.5
L 0.65 B
K x , sur 6.8 2.4 K x , sur K y , sur GL 1
2 0.75
GB 0.2
B L
Translation along x-axis
L B
K zz , sur GB 3 4.25 4.06 K zz , sur GJ t0.75 4 111
2.45 10
B L
Torsion about z-axis
I y 3
L 0.75 L
3.73 0.27
2.4 0.15
GB 3
1 1
G
B B
Rocking about y-axis K yy , sur K yy , sur
L B
K xx , sur 3.2 B 0.8 I x 2.4 0.5 L
0.25
GB 3
1 1
G L
B
0.75
Rocking about x-axis K xx , sur
0.25 D D B Aw
z 1.0 0.25 z 1 1 1.3 1 0.2
0.8 2/3
L / B B 21B L 4 BL
Translation along z-axis
1.34 D D z A
y 1.0 0.33 y 1 0.15 1 0.52 2
0.8 0.4
1 L / B B
w w
Translation along y-axis
B BL
Translation along x-axis x y Same equation as for y, but Aw term changes for B L
1.32 D B d w
zz 1 1.3 zz 1 1.4 1
L B
0.9 0.9
L / B B
Torsion about z-axis
D D 2 dw dw B
0.6
B
1.6
B
Rocking about y-axis 1.5
D D dw dw dw B
0.2
xx 1.0 xx 1 1.26 1
B D
2
1.6
B 0.35 L / B B B L
Rocking about x-axis
4 L / B a 0
K y , sur / GB 2 y
Translation along y-axis y 1.0 y
4 L / B a
x 1.0 x
K x , sur / GB 2 x
Translation along x-axis 0
0.33 0.03 L / B 1 a 0
2
4 / 3 L / B L / B a 02 a
zz 1.0 zz 0
3
2 zz
K zz , sur / GB
2
Torsion about z-axis
1 0.33 L / B 1
0
1 3 L / B 1
0.8
1.4
a 0
3
a 2
0.7
4 / 3 L / B a a0
yy 1.0 yy
3 2
0.55a 02
1.4 2 K yy , sur 2 yy
2
0
Rocking about y-axis
0.6 a0 a 0
1 1.75 L / B 1
L / B
1.8
GB
3 3
0.55 0.01 L / B 1 a 0 4 / 3 L / B a 02 a
xx 1.0 xx 0
2
2 xx
K xx , sur / GB 2.2
0.4 0.4
Rocking about x-axis
2.4 a0
0
L / B
L / B
2
3 2
a
3 3
4 L / B D / B 1 L / B a
z
K z ,emb / GB
0
2 z
Translation along z-axis
4 L / B D / B 1 L / B a
y
K y,emb / GB
2 y
0
Translation along y-axis
4 L / B D / B L / B a
x
K x,emb / GB
0
2 x
Translation along x-axis
4 / 3 3 L / B D / B L / B D / B 3 L / B D / B D / B L / B ( L / B) a 02 a
zz 0
3 2 3
zz ,emb 2 zz
Torsion about z-axis
3
a 02
GB 1 3 L / B 1
K 1.4
0.7
L 3 D L 3
4 / 3 3 a 02 4 L D
D L D D L
3 3 2 3
B B B B B B B B a
yy B
3 B
yy
0
Rocking about y-axis
3 1 1.75 L / B 1 a 0 3
K 2
K 1.8
yy , emb
yy , emb
2
GB GB
D D 3 L 3
4 / 3 3 a 02 4 L 1 D
L D D L
3
B B B B B B B
B B a0
xx
2 xx
3
K xx , emb K xx , emb
Rocking about x-axis
GB 1 1.75 L / B 1
3
1.8
GB
2
a 0 3
Notes: Soil hysteretic damping, s, is additive to foundation radiation damping, j.
; from Table 2-3a
a 0 B / V s ; 2 1 / (1 2 ) ; 2.5
Figure 2-4 shows the variation in dynamic stiffness modifiers versus frequency for
rigid footings located at the ground surface. In the case of translational stiffness,
dynamic stiffness modifiers (x, y) are essentially unity, regardless of frequency or
foundation aspect ratio. For rotational stiffness, however, dynamic stiffness
modifiers for rocking (xx, yy) degrade markedly with frequency, but are relatively
insensitive to aspect ratio.
a 0 B / V s a 0 B / V s
Figure 2-4 Plot of dynamic stiffness modifiers and damping ratios versus dimensionless frequency,
hysteretic damping, and =0.33: (a) geometry; (b) x-direction; and (c) y-direction.
for rectangular footings resting on the surface of a homogeneous half-space, with zero
Because soil hysteretic damping, s, is taken as zero, Figure 2-4 also shows the
variation in radiation damping ratios for translation (x, y) and rotation (xx, yy)
versus frequency. Translational radiation damping is only modestly affected by the
direction of shaking or the aspect ratio of the foundation. The modest increase of
translational damping with aspect ratio is a result of the increased foundation size
(i.e., larger wave source).
j term should be taken as the sum of radiation damping for the appropriate vibration
source of foundation damping. When used to calculate the dashpot coefficient, cj, the
mode (from Equation 2-13a and Table 2-3) and soil hysteretic damping, s, provided
by a geotechnical engineer.
Figure 2-5 shows the variation in dynamic stiffness modifiers and radiation damping
ratios versus frequency for embedded foundations. In equations provided by Pais and
Kausel (1988), dynamic stiffness modifiers are unaffected by embedment, and this
a 0 B / V s
Figure 2-5 Plot of dynamic stiffness modifiers and damping ratios versus
dimensionless frequency, for square footings embedded in a
The elasto-dynamic analyses upon which Figure 2-5 is based assume perfect contact
between soil and basement walls. Accordingly, the solutions indicate much higher
damping levels than those for shallow foundations (Gazetas, 1991). These damping
levels may not be reliable when gaps form between foundations and the adjacent soil,
which reduces the potential for radiation damping from basement walls. In studies
performed by Stewart et al. (1999b), buildings shaken by earthquakes generally do
not exhibit damping levels consistent with such models. As a result, the impedance
of embedded foundations can be conservatively estimated from the equations for
static stiffness in Table 2-2a and adjusted by dynamic stiffness modifiers for surface
foundations from Table 2-3b.
In most cases, Vs profiles are evaluated away from foundations (i.e., in the free-field)
and reflect a variation of shear modulus with depth. Variation in soil shear modulus
with depth, and the presence of additional weight from a structure, complicates the
selection of an appropriate shear wave velocity in the calculation of static foundation
stiffnesses. To evaluate a single effective Vs value for use in computations, it is
necessary to: (1) correct Vs values measured in the free-field to account for
overburden pressures associated with the added weight of the structure; and (2)
calculate an average effective Vs value over an appropriate depth range.
m , as follows:
Soil shear modulus, G0, is known to increase with mean effective confining stress,
G0 G a m
n
pa
(2-16)
where Ga is the shear modulus at atmospheric pressure, m is the effective confining
stress, pa is taken as approximately 100 kPa, and n varies from approximately 0.5 for
granular soils (Hardin and Black, 1968; Marcuson and Wahls, 1972) to 1.0 for
cohesive soils with plasticity index (PI) greater than 6.5 (Yamada et al., 2008).
Recognizing that Vs is proportional to the square root of shear modulus, free-field
measurements of shear wave velocity (at depth) can be corrected to account for
overburden pressures due to the added weight of the structure as follows:
(z ) v (z )
Vs , F (z ) Vs (z ) v
n/ 2
v (z )
(2-17)
effective stress from the self-weight of the soil at depth z, and v (z) is the increment
of vertical stress due to the weight of the structure at depth z, which can be computed
using classical Boussinesq stress distribution theory (e.g., Fadum, 1948). The
overburden correction in Equation 2-17 is typically significant only at shallow depths
(i.e., 50% to 100% of the foundation dimension).
Vs , avg
zp
z
(2-18a)
i
Vs ,F (z ) i
n
i 1
When soil shear modulus increases with depth, some of the seismic energy radiating
from the foundation reflects upward towards the foundation, hence it is not “lost” as
occurs in a uniform half-space. Impedance solutions for this phenomenon are
Figure 2-6 shows a plot of dynamic stiffness modifiers and radiation damping ratios
comparing results for a uniform half-space and non-uniform profiles in which G
varies with depth, as shown. The effect on radiation damping is more pronounced in
rotational stiffness modifier, yy, for square foundations (Figure 2-6b) is modest.
rotation (Figure 2-6b) than in translation (Figure 2-6a). Also, the effect on the
Hence, the effect of variation in soil shear modulus with depth is most critical for
static stiffness and radiation damping associated with foundation rocking. Because
rocking is often an insignificant contributor to overall foundation damping, the
practical impact of soil non-homogeneity is primarily related to its effect on static
stiffness.
a 0 r / v s 0 a 0 B / v s 0
Figure 2-6 Plot of dynamic stiffness modifiers and radiation damping ratios
versus dimensionless frequency comparing uniform half-space and
non-uniform profiles in which G varies with depth: (a) translation for
circular foundations (from Gazetas, 1991); and (b) rotation for square
foundations (from Vrettos, 1999).
In the extreme case of a rigid material at depth in a soil profile, radiation damping
from body wave propagation disappears at frequencies lower than the fundamental
frequency of the soil column. While no geologic materials are actually rigid, this can
Classical impedance function solutions, such as those presented in Table 2-2 and
Table 2-3, strictly apply for rigid foundations. As illustrated in Figure 2-1, soil-
foundation interaction for rigid foundations can be represented by individual springs
for each foundation degree of freedom. Actual foundation slabs and basement walls,
however, are non-rigid structural elements. The few theoretical solutions that exist
apply to circular foundations supporting a rigid core (Iguchi and Luco, 1982),
flexible perimeter walls (Liou and Huang, 1994), or rigid concentric walls (Riggs and
rotational stiffness, krr, and rotational radiation damping ratio, rr, for the cases of a
Waas, 1985). Figure 2-7 shows the effect of flexible foundation elements on
a 0 r f / V s a0
rigid core (Iguchi and Luco, 1982) and flexible perimeter walls (Liou and
Huang, 1994).
where rf is the foundation radius, tf is the foundation thickness, and Ef and f are the
Young’s modulus and Poisson’s ratio of the foundation concrete. The case of = 0
corresponds to a rigid foundation slab.
Liou and Huang (1994) showed that foundation flexibility does not significantly
affect translational stiffness and damping terms for the case of flexible perimeter
walls. For rotational stiffness and radiation damping, Figure 2-7 shows that
foundation flexibility effects are relatively modest for the case of flexible perimeter
walls, and most significant for the case of a rigid core.
Typical practice does not adjust the impedance function for non-rigid foundations as
shown in Figure 2-7. Instead, foundations springs are distributed across the extent of
the foundation, as illustrated in Figure 1-2c. Distributed springs allow the foundation
to deform in a natural manner given the loads imposed by the superstructure and the
spring reactions. For vertical springs, this can be accomplished by calculating the
vertical translational impedance, as described above, and normalizing it by the
foundation area to compute stiffness intensity, k zi (also known as coefficient of
subgrade reaction), with dimensions of force per cubic length:
k zi
kz
(2-20a)
4 BL
c zi
cz
(2-20b)
4 BL
To correct for underestimation of rotational stiffness, strips along the foundation edge
(of length ReL) are assigned stiffer springs. When combined with springs in the
interior, the total rotational stiffness of the foundation is reproduced. Harden and
Hutchinson (2009) present expressions for end length ratios and spring stiffness
increases as a function of L/B using static stiffnesses from Gazetas (1991).
More generally, the increase in spring stiffness, Rk, can be calculated as a function of
foundation end length ratio, Re, as:
3k yy
i 3 1 R e
z
3
1 1 R e
4 k BL
Rocking ( yy ) : R k , yy 3
(2-21a)
3k xx
i 3 1 R e
z
4k B L
3
1 1 R e
Rocking ( xx) : R k , xx 3
(2-21b)
Equations 2-21 were derived by matching the moment produced by the springs for a
unit foundation rotation to the rotational stiffness kyy or kxx. In these equations, a
value of Re can be selected (typically in the range of 0.3 to 0.5), which then provides
a unique Rk. This correction for rotational stiffness, however, does not preserve the
original vertical stiffness kz. This is considered an acceptable approximation, in
general, because rocking is the more critical foundation vibration mode in most
structures.
1 R
Rocking ( yy) : R c , yy
R k , yy 1 1 R e
4c zi BL3
3 3
(2-21c)
e
1 R
Rocking ( xx) : R c , xx
R k , xx 1 1 R e
4c zi B 3 L
3 3
(2-21d)
e
Use of the above procedures for modifying vertical spring impedances will reproduce
the theoretical rotational stiffness and damping through distributed vertical springs
and dashpots. While this allows foundation flexibility to be accounted for, in the
sense that foundation structural elements connected to springs and dashpots are non-
rigid, a question that remains is whether or not the rotational impedance computed
using a rigid foundation impedance function is an appropriate target for calibration.
For the case of a rigid core illustrated in Figure 2-7 it is not, but solutions for more
practical situations are not available.
the embedment modifier, kx/x) can be applied to the spring at the base slab
the embedded stiffness attributable to the base slab (i.e., the stiffness without
Previous sections have discussed linear spring stiffness, but have not addressed
foundation capacity. Pushover procedures for SSI analysis typically utilize elastic-
perfectly-plastic force-deflection relationships for springs. Hence, limiting spring
forces (i.e., capacities) are needed for vertical and horizontal springs.
In the case of vertical springs, the capacity is the unfactored bearing capacity of the
foundation distributed over the tributary area of the spring (dA). Bearing capacity
should be calculated considering the foundation geometry, drained or undrained shear
strength parameters as appropriate, soil unit weight, and simultaneous presence of
both horizontal and vertical loads on the foundations. These concepts are discussed
In the case of horizontal springs, the capacity of springs located at the level of a
footing or mat should reflect the unfactored sliding resistance at the slab-soil
interface. The capacity of springs along basement walls should reflect the unfactored
passive earth pressure.
Limiting lateral and vertical capacities of foundations are usually not simultaneously
realizable. This is especially important in the presence of geometric nonlinearities
such as soil-foundation gapping, described further in Section 2.4.
Buildings founded on soft soils may have pile-supported footings or mats, especially
when the foundation is not embedded (i.e., no basement levels). This section
discusses the effective stiffness and damping of pile-supported foundations. Only the
case of vertical piles is considered, as battered piles are seldom used in building
structures.
The impedance of pile groups for lateral and rotational vibration modes using elasto-
dynamic solutions is also covered. Pile-to-pile interaction effects are considered and
the manner by which vertical responses of piles are combined to develop rotational
impedance is discussed. Pile stiffnesses from elasto-dynamic solutions are contrasted
with the discrete element modeling typically performed in practice.
An important consideration when piles are combined with shallow spread footings or
a mat foundation is whether or not lateral resistance is provided by the shallow
foundation elements in combination with the piles. Soil might be expected to settle
Springs and dashpots effectively replace a single pile in the numerical modeling of a
pile-supported foundation, as schematically illustrated in Figure 2-9. The impedance
of single piles, represented by k jp and c jp , can be described in the notation used for
shallow foundations (Equations 2-12 to 2-14a). The dynamic stiffness for a
particular vibration mode is denoted k jp , and the corresponding dashpot, representing
the effects of damping, is denoted c jp . Subscript j represents the vibration mode,
which is taken as x (horizontal) and z (vertical).
k jp K jp p
j (2-22a)
where:
K jp j E s d
j w pj w sj w bj f E p E s , L p d
jp f E p E s , p s , w sj , , a 0p
(2-22b)
s and p are the mass densities for soil and pile materials, respectively; is the
diameter; Es and Ep are the Young’s moduli for soil and pile materials, respectively;
long pile does not deflect over its entire length, but only to a certain depth, termed the
active pile length, La (shown in Figure 2-10).
a 0p d / V s a 0p d / V s
Figure 2-10 Plot of dynamic stiffness modifiers and damping ratios versus
The active pile length is on the order of 10 to 20 pile diameters, depending on pile-
soil stiffness contrast, soil non-homogeneity, and fixity conditions at the pile head
(Randolph, 1981; Gazetas, 1991; Pender 1993; Syngros 2004; Karatzia and
Mylonakis 2012). Piles with lengths Lp > La essentially behave as infinitely long
beams, and the actual length does not affect flexural response. Active lengths tend to
be greater for dynamic loading than for static loading, due to the ability of elastic
waves to travel further down the pile than a static stress field.
Expressions for active pile length for lateral deformations can be cast in the form:
Ep
La d
Es
(2-24)
where and μ are dimensionless constants, and all other terms are as previously
(1981) and Fleming et al. (1992) recommend = 1.8 and μ = 0.25. For dynamic
defined. For fixed-head piles in homogeneous soil under static loading, Randolph
G0, is then evaluated from shear wave velocity, Vs, as G0 sVs2 (from Equation
values of Vs for use in the analysis of pile impedance. Small-strain shear modulus,
2-9). Soil shear modulus, G, should be reduced relative to G0 for large strain effects
E s 2 1 G
(e.g., Table 2-1). Young’s modulus for soil can then be related to shear modulus as:
(2-25)
As an alternative to the use of Table 2-1 for evaluating modulus reduction effects,
soil strains adjacent to a horizontally loaded pile can be evaluated as (Kagawa and
1 u z
Kraft, 1981):
z (2-26)
2.5d
where γ(z) denotes an average soil shear strain at depth z, and u(z) the corresponding
horizontal pile displacement. On the basis of this equation, a strain-compatible soil
shear modulus can be obtained through conventional modulus reduction curves (e.g.,
Darendeli, 2001;Vucetic and Dobry, 1991).
Weight factors in Equation 2-22b (wpj, wsj, and wbj) represent the relative
contributions of the pile structural stiffness, pile-soil interaction through side-load
transfer, and pile-soil interaction through toe resistance for vibration mode j. These
weight factors always sum to unity (i.e., wpj + wsj + wbj = 1.0), and are not required
The equations are used to determine dimensionless parameters j and weight factors
Equations related to the static stiffness of single piles are provided in Table 2-4a.
(wpj, wsj, and wbj). Note that the equations for j depend, in turn, on a series of
are and Lp (both related to Ep/Es), and the corresponding modulus, z. Equations
Ep/Es as indicated in the table. For vertical vibration (j=z), the additional variables
for dynamic stiffness modifiers ( jp terms) and damping ratios ( jp terms) for single
piles are provided in Table 2-4b.
Damping ratios reflect material damping in the pile and soil materials (p and s,
respectively) as well as radiation damping (rj terms). The dynamic stiffness
modifiers and radiation damping ratios for single piles obtained from these
expressions are plotted in Figure 2-10.
When piles are used as part of a building foundation system, they are usually
configured in groups to support continuous mat foundations or discrete pile caps for
individual load-bearing elements. The impedance of a pile group cannot be
determined by simple addition of individual pile impedances because grouped piles
interact through the soil by “pushing” or “pulling” each other through waves emitted
from their periphery. This is called a group effect, and it can significantly affect the
impedance of a pile group as well as the distribution of head loads among individual
piles in the group. Group effects depend primarily on pile spacing, frequency, and
number of piles. They are more pronounced in the elastic range, and dynamic group
effects decrease in the presence of material nonlinearity.
The ratio of the pile group impedance in any oscillation mode, k Gj , to the sum of the
individual static pile impedances in the same oscillation mode, N piles k Pj , static , is the
efficiency factor of the pile group (Kaynia and Kausel, 1982). Efficiency factors are
generally less than unity for low frequencies, but can increase significantly at higher
frequencies under low strain conditions. Negative efficiency factors are also
possible, which suggests a phase difference of over 90 degrees between oscillations
of a single pile and oscillations of the pile group at the same frequency. Note that
these factors strictly refer to dynamic compliance of the group and are different from
the familiar efficiency factors for group bearing capacities in foundation engineering.
E
Poulos and Davis (1980)
x 1/ 4 x 3/ 4 p
1/ 4
1
Es
Translation along x-axis Scott (1981), Mylonakis (1995)
2
E
3/ 40
x 2 p
Dobry et al. (1982)
Es
Syngros (2004)
w px 1/ 4
Dobry et al. (1982)
w sx 3 / 4 Long pile ( L p / d L active ) Mylonakis (1995)
w bx 0 Mylonakis and Roumbas (2001)
Ep Ω tanh( L p )
Randolph and Wroth (1978)
z z
1/ 2 1/ 2
2 Es 1 Ω tanh( L p )
Translation along z-axis Scott (1981)
Ep
1/ 2
1 E s
Ω
2
z 2
Mylonakis and Gatezas (1998),
Randolph (2003), Salgado (2008)
4 2 E Lp
L p z p
1 1
Es
2
d
w pz 1 w sz w bz
Thomas (1980)
2
2 cosh L p 1 2 sinh L p
w bz
a o
3 p / s p 2
xp 1
32 x 1
Translation along x-axis Mylonakis and Roumbas (2001)
a o a o
p / s p 2 1 p 1/ 2
zp 1 w sz
8 x 1
Translation along z-axis Mylonakis (2011)
2
Degree of Freedom Damping Ratio Reference
xp p s rx
1 3 3 Dobry et al. (1982)
Translation along x-axis
4 4 4 Mylonakis and Roumbas (2001)
a o
p 3/ 4
rx
2 x 1 x
3
Gazetas and Dobry (1984a,b)
1
a op w bz 0.21a op
1.2
rz p sz
z 4 1 z
3/ 4
w Gazetas and Dobry (1984b)
Results for horizontal and rocking oscillations are provided in Figure 2-11 for pile
groups in square configurations computed using the solution by Mylonakis and
Gazetas (1999). Peaks and valleys observed in the plots are due to destructive and
constructive interference of the waves between piles, which tend to increase and
decrease in dynamic impedance, as first identified by Wolf and Von Arx (1978), and
explained by Kaynia and Kausel (1982), Nogami (1983), and Dobry and Gazetas
(1988).
The above effects tend to decrease with non-homogeneity and nonlinearity in the soil,
as the waves emitted from the periphery of the piles become less coherent (El-Naggar
and Novak, 1994 and 1995; Michaelides et al., 1998). At the low normalized
frequencies of interest in most practical problems ( a0p 0.4) , efficiency factors are
less than one and saturate as the number of piles, Npiles, increases. Hence, for large
piles groups, low-frequency efficiencies would not be significantly lower than those
shown for the 44 pile group in Figure 2-11. Efficiency factors above unity for the
rocking mode are due to the intrinsic out-of-phase movement of the piles located on
opposite sides of the rocking axis.
Figure 2-11 Plots of pile group efficiency factors and damping ratios versus dimensionless
frequency for square pile groups for: (a) lateral loading at head of pile group
under zero cap rotation; (b) moment at head of pile group, introducing rocking
and rocking results are for Ep/Es = 1000, Lp/d = 20, p/s = 1.3, =0.4, (pile
under zero cap translation; and (c) vertical loading at head of pile group. Lateral
spacing)/d = 5, p = 0, and s = 0.05. Vertical resutls are for Ep/Es = 100, Lp/d =
15, p/s = 1.4, = 0.4, (pile spacing)/d = 5, p = 0, and s = 0.05.
In Figure 2-12, the group impedance of a 33 pile group is compared to that of a
footing of equivalent dimensions. Results are presented relative to a common
normalized frequency using the conversion:
a 0 a 0P
B
(2-27)
d
Results show that the effect of the piles is dramatic for rotational stiffness, increasing
kyy by factors up to 50 relative to a shallow footing alone. Translational stiffness is
also increased, but to a lesser degree. Interestingly, damping ratios are on the same
order of magnitude, with somewhat greater variation with frequency caused by pile-
to-pile interaction effects.
a 0 a 0p (B / d ) a 0 a 0p (B / d )
Figure 2-12 Comparison between the impedance (stiffness and damping ratio) of
a 3x3 pile group and the impedance of a footing with equivalent
dimensions for: (a) lateral loading at head of pile group; and (b)
s = 0.05.
In engineering practice, piles are typically modeled using discrete element methods,
as prepared in common commercial computer programs such as APILE (Ensoft, Inc.
program for analysis of axial capacity of piles) and LPILE (Ensoft, Inc. program for
analysis of piles under lateral loads). In such programs, it is straightforward to
extract head load-deflection relationships from the analysis, which are often used in
lieu of elasto-dynamic solutions described above.
Unfortunately, the p-y (nonlinear lateral spring) and t-z (nonlinear vertical spring)
relations used in these programs were developed to represent large-deformation
response in static or cyclic loading, and they do not accurately capture stiffness or
damping. In fact, many spring formulations have an initial stiffness that is infinite
(i.e., the backbone curve is vertical at the origin). When the pile-soil system is not
expected to yield, the elasto-dynamic solutions presented above provide a superior
Most of the research performed on nonlinear SSI has been related to structural
yielding with linear, or equivalent-linear, soil (Case 1, above) or soil
yielding/gapping with a linear structure (Case 2/Case 3, above). A brief overview of
this work is presented in the following sections.
global or system ductility demand, s; or (2) the conventional member ductility
The performance of yielding structures is typically represented in terms of: (1) the
demand, (Priestley and Park, 1987; Paulay and Priestley, 1992). For the structure
in Figure 2-1, these factors are defined as:
s max
y
(2-28)
max
y
(2-29)
A key difference between these factors is that member ductility demand, , refers
exclusively to structural deformations, whereas global or system ductility demand, s,
encompasses rigid body movements associated with translation and rotation of the
foundation, which do not reflect strains in the superstructure. Equation 2-28 and
Equation 2-29 are geometric relations, and the former always provides smaller
numerical values than the latter for a given set of structural response values
(Ciampoli and Pinto, 1995; Mylonakis and Gazetas, 2000). Because of this, the use
of Equation 2-28 for assessing SSI effects on structural response might be of limited
value.
There is mounting analytical and experimental evidence that material and geometric
nonlinearities in the soil may be beneficial to the seismic response of a structure.
This has led some authors (e.g., Gazetas, 2006; Gajan and Kutter, 2008) to propose
revising the foundation design philosophy by allowing significant yielding in the soil
close to the foundation, or the foundation itself, to dissipate energy and protect the
superstructure. This requires control of settlement and tilting of the structure.
Hence, the analysis and design process considering soil nonlinearity involves
optimization of the trade-offs between the potentially beneficial effects of soil
Soil-structure interaction studies with nonlinear soil and foundation behavior can be
classified into three approaches: (1) continuum models, (2) beam-on-nonlinear
winkler foundation (BNWF) models, and (3) plasticity-based macro-models. The
first approach is by far the most computationally demanding, and has been employed
to a limited extent (Borja and Wu, 1994; Jeremic et al., 2009). Available findings
suggest the creation of stress-induced inhomogeneities under the foundation, which
may limit wave radiation away from the structure and cause wave reflections leading
to resonance effects.
The second and third approaches for nonlinear soil modeling are briefly described
below. Although both approaches can consider material nonlinearities, only macro-
element approaches are currently configured for material nonlinearities (gapping).
The basic description of the models, their input parameters, and a brief comparison
with experimental data are provided. The reader is referred to related literature for
additional information. The emphasis in the following discussion is on two specific
models that are implemented in OpenSees, Open System for Earthquake Engineering
Simulation (McKenna, 1997; OpenSees, 2011). Their implementation could be
reasonably extended to other computational platforms in the future. Much of the
content is adapted from Gajan et al. (2010).
Starting with the pioneering work of McClelland and Focht (1958), beam-on-
nonlinear winkler foundation (BNWF) models have been used for many years for
analyzing the response of foundations, most notably piles, for static loads (Matlock,
1970; Cox et al., 1974) and dynamic loads (Penzien, 1970; Nogami et al., 1992;
Boulanger et al., 1999). Key advantages of these models over continuum
formulations lies in their ability to describe soil-structure interaction phenomena by
one-dimensional nonlinear springs distributed along the soil-foundation interface. It
is well-known that the modulus of the springs (also known as modulus of sub-grade
reaction) is not uniquely a soil property, but also depends on foundation stiffness,
geometry, frequency, response mode, and level of strain. A limitation of the
approach relates to its one-dimensional nature. A spring responds only to loads
acting parallel to its axis, so loads acting in a perpendicular direction have no effect
on the response of the spring. Accordingly, the concept of plastic potential and flow
rule cannot be explicitly incorporated. Nevertheless, the BNWF approach is popular
because of its simplicity and predictive abilities on a variety of problems.
The material models are mechanistic, based on an arrangement of various linear and
nonlinear springs, gap elements, and dashpots. Radiation damping can be accounted
for using a dashpot that is placed in parallel with the far-field elastic component. The
backbone curves are thus characterized by a linear-elastic region, followed by an
increasingly growing nonlinear region. The QzSimple2 material has an asymmetric
hysteretic response, with a backbone curve defined by an ultimate load on the
compression side and a reduced strength in tension to account for the low strength of
soil in tension. The PxSimple1 material is envisioned to capture the passive
resistance, associated stiffness, and potential gapping of embedded shallow footings
subjected to lateral loads. This material model is characterized by a pinched
hysteretic behavior, which can more suitably account for the phenomena of gapping
during unloading on the opposite side of a footing. The TxSimple1 material is
intended to capture the frictional resistance along the base of a shallow foundation.
This material is characterized by a large initial stiffness and a broad hysteresis, as
anticipated for frictional behavior associated with foundation sliding.
The functional forms and parameters describing the p-x, t-x, and q-z springs are
similar, so only the q-z model is described here. The backbone curve has linear and
nonlinear regions. The linear-elastic portion of the backbone curve is described by
the initial stiffness kz:
q k zs (2-30)
where q represents the spring force, and s represents the spring deflection. The upper
limit of the linear-elastic region, defined as q0, is taken as a fraction of the ultimate
load qult as follows:
q 0 C r q ult (2-31)
q q ult q ult q0 (for |s| > s0)
n
cs 50
cs 50 s s 0
(2-32)
The unload-reload rules that operate with the backbone curve are relatively simple,
generally consisting of the familiar Masing rules (i.e., the shape of the unload and
reload portion of the cyclic loop matches twice the backbone curve). The drag and
gap component is parameterized by a bilinear closure spring in parallel with a
nonlinear drag spring. The cyclic response of each of the material models, when
subjected to a sinusoidal displacement, is demonstrated in Figure 2-14.
User-defined parameters for the q-z element can be synthesized based on two
physical parameters obtained from the results of a typical high-quality geotechnical
site investigation (i.e., bearing capacity, qult, and elastic stiffness, kz) and several
relatively subtle parameters defining the details of the elements described above. The
need for several relatively unfamiliar parameters presents a barrier to implementation
of this type of model in practice. These parameters include:
Radiation damping (cz). This dashpot coefficient is considered to be a physical
parameter that is well documented in the literature (e.g., Section 2.2). The
parameter is sensitive to soil stiffness, footing shape, aspect ratio and
embedment.
Tension capacity (TP). The tension capacity parameter, TP, determines the
maximum magnitude of the drag force in Component 1 of the nonlinear springs.
It is the ratio of tension capacity to bearing capacity with typical selected values
of 0 to 0.10 (as suggested in Boulanger et al., 1999), although, more recently
some experts (e.g., Kutter) have recommended using a TP value of zero.
Distribution and magnitude of vertical stiffness. As illustrated in Figure 2-8
and Figure 2-13, two parameters are necessary to account for the distribution and
building, with a vertical factor of safety, FSv = 2.3, resting on dense sand with Dr =
80%. In these tests, M is the moment applied by the horizontal force, V, at the level
of the soil-foundation interface. Figures 2-15a through 2-15d show moment versus
rotation, settlement versus rotation, shear force versus sliding, and settlement versus
sliding histories, as predicted by the model and compared with the experiment. The
model compares reasonably well with the experimental results in terms of capturing
the shapes of the hysteresis loops and rotational and lateral stiffnesses. Additional
comparisons are provided in Raychowdhury (2008).
plastic approach or a simple G- and s- correction based on a characteristic shear
(3) a mechanism for describing plastic flow is incorporated in the form of a hypo-
strain; and (4) uplift behavior can be modeled by means of a nonlinear model
allowing for separation between the footing and soil (Pecker and Chatzigogos, 2010).
While PBM models are rational and can capture plastic effects, available
formulations possess a number of drawbacks, notably an inability to incorporate
flexible foundation behavior, the effects of stress-induced inhomogeneity on
radiation damping, and failure modes other than general shear failure. In addition,
there is limited experimental validation of the models. More details are given in
Cremer et al. (2001), Houlsby and Cassidy (2002), Chatzigogos et al. (2009) and
Pecker and Chatzigogos (2010). A variant of these models, focusing on uplift
phenomena in compliant soil, is examined in some detail below.
The model differs from other macro-element models in that its constitutive relations
are obtained by tracking the geometry of gaps and contacts of the soil-foundation
interface. To this end, the CIM provides nonlinear relations between cyclic loads and
displacements of the soil-foundation system during combined cyclic loading (i.e.,
vertical, shear, and moment).
Soil-foundation contact is tracked in the CIM using a parameter called the critical
contact area ratio, A/Ac, where A is the area of the footing, and Ac is the area of the
footing required to have contact with the soil to support the vertical and shear loads.
The ratio A/Ac can be considered to be an alternate definition of the factor of safety
with respect to bearing capacity. For a two-dimensional shear wall structure loaded
in the plane of the wall, the ratio A/Ac equals the footing length ratio 2L/Lc, as shown
in Figure 2-17a. In Figure 2-17b, the foundation position is tracked relative to the
underlying soil surface, which is pushed as far as the “soil_max” surface, but which
rebounds to the “soil_min” surface when unloaded. Zero stress transfer between the
soil and foundation occurs in the gap region.
With seven user-defined input parameters, the CIM is intended to capture the
essential features (load capacities, stiffness degradation, energy dissipation, and
permanent deformations) of the cyclic load-deformation behavior of shallow
foundations. One advantage of the CIM relative to the BNWF model is that the
(a) (b)
Figure 2-17 Contact interface model (CIM): (a) definition of critical contact length; and (b)
tracking of foundation position relative to soil pressures (Gajan and Kutter, 2009,
with permission from ASCE).
Input parameters for the CIM are the ultimate vertical load, QULT, the length of
elastic rotation limit, elastic, the rebound ratio, Rv, and the internal node spacing, DL.
footing, 2L, the initial vertical stiffness, Kv, the initial horizontal stiffness, Kh, the
stiffness and footing geometry. Of those parameters, elastic, Rv, and DL are relatively
The initial rotation stiffness is calculated by the CIM based on the specified vertical
Figure 2-18 shows a comparison of the CIM model to results from centrifuge test
SSG02–05 conducted at the University of California at Davis (summarized in Gajan
resting on dense sand with Dr = 80%. Figures 2-18a through 2-18d show moment
versus rotation, settlement versus rotation, shear force versus sliding, and settlement
versus sliding histories, as predicted by the model and compared with the experiment.
The CIM demonstrates good comparison with the experimental results in terms of
capturing the shapes of the hysteresis loops and rotational and lateral stiffnesses.
Additional comparisons are presented in Gajan (2006) and Gajan et al. (2008).
800 150
Shear force (KN)
Moment (KNm)
100
400
50
0
0
-400
-50
-800 -100
Experiment
40 40
CIM Simulation
Settlement (mm)
Settlement (mm)
0 0
-40 -40
-80 -80
-120 -120
-160 -160
-0.04 -0.02 0 0.02 0.04 -40 -20 0 20 40
Rotation (rad) Sliding (mm)
Figure 2-18 Comparison of CIM simulation to centrifuge tests for a medium aspect ratio
building on dense sand, with Dr = 80%, and FSv = 2.6: (a) moment-rotation;
(b) settlement-rotation; (c) shear-sliding; and (d) settlement-sliding (Gajan
and Kutter, 2009, with permission from ASCE).
Gajan et al. (2010) compared the performance of the BNWF and CIM nonlinear SSI
models for hypothetical site conditions and experimental datasets. Based on those
comparisons, recommendations for model selection are as follows:
If the simulations are to be used for structural design of footing elements, or the
footing flexibility is anticipated to contribute to the foundation response, the
BNWF model should be chosen. This model can be used to more directly
evaluate internal moments and shears used for section design.
If the normalized moment to shear ratio M / (V L) is less than approximately
3.0, and sliding is not restrained by slabs and grade beams, then the moment
capacity of the footing will be sensitive to shear load, and vice versa. In this
between the moment, shear, and axial responses. For cases with M / (V L)
case, the CIM model is preferred because of its ability to account for coupling
greater than 3.0, rocking will tend to dominate, and both models should produce
may also be neglected for very small M / (V L) ratios where sliding is known
similar results if the parameter selection protocols herein are followed. Coupling
Base-slab averaging results from adjustment of spatially variable ground motions that
would be present within the envelope of the foundation, which are averaged within
the foundation footprint due to the stiffness and strength of the foundation system.
Base-slab averaging can be understood by recognizing that the motion that would
have occurred in the absence of the structure is spatially variable. Placement of a
foundation slab across these variations produces an averaging effect in which the
foundation motion is less than the localized maxima that would have occurred in the
free-field. Torsional rotations, referred to as the “tau effect” (Newmark, 1969), can
also be introduced.
Motions of surface foundations are modified relative to the free-field when seismic
waves are incoherent. Incoherence of the incident waves at two different points
means that they have variations in their phase angle. Some incoherence is
deterministic (i.e., predictable), because it results from wave passage. For example,
as illustrated in Figure 3-1a, the presence of a non-zero vertical angle causes waves to
arrive at different points along the foundation of a building at different times. This is
Incoherence that remains when waves are aligned to have common arrival times is
stochastic, and is quantified by lagged coherency models. Stochastic incoherence
results from source-to-site heterogeneities in the seismic path of travel, which scatters
seismic waves. Lagged coherency is also well-documented in array studies (e.g.,
Abrahamson et al., 1991; Ancheta et al., 2011). As a practical matter, incoherence
from wave passage and lagged coherency is always present in earthquake ground
motions to some degree.
Figure 3-1 Illustration of foundation subjected to inclined shear waves: (a) schematic
In the presence of incoherent wave fields, translational base-slab motions are reduced
relative to the free-field, and rotational motions are introduced. The reduction in
translational motion is generally the more important result. Reductions of base-slab
translation and the introduction of torsion and rotation in the vertical plane are effects
that tend to become more significant with increasing frequency. The frequency-
dependence of these effects is primarily associated with: (1) the increased effective
size of the foundation relative to the seismic wavelengths at high frequencies; and (2)
significant reductions in lagged coherency with increasing frequency (Abrahamson et
al., 1991).
There are numerous theoretical models for predicting the relationship between
foundation input and free field ground motions in the presence of inclined, but
V
sin a 0k s
V app V app
Hu
, a 0k
V
(3-1b)
a 0k s
V
2 Vs
app
V app
Hu , a 0k
2
(3-1c)
2 Vs
In the above expressions, a 0k is similar to a0 defined in Equation 2-15, except that the
foundation dimension is related to base contact area, as follows:
B eA
a k
0 (3-2)
Vs
If, as indicated in array studies, Vapp ranges from approximately 2.0 km/s to 3.5 km/s,
then for a typical soil site, a reasonable estimate of the velocity ratio, Vapp/Vs, is
approximately 10. In Figure 3-1b, the result labeled “wave passage only” shows the
transfer function between uFIM and ug based on Equations 3-1. Using this model,
wave passage alone causes relatively modest base-slab reductions in ground motion
across the frequency range of engineering interest.
Transfer functions of recorded foundation input and free-field motions are generally
significantly lower at high frequencies than predicted by wave passage models. This
occurs because wave passage is a relatively modest contributor to the spatial
variation in ground motion that drives base-slab averaging. Additional sources of
variability include stochastic phase variability (quantified by lagged coherency) and
stochastic variations in Fourier amplitudes. Two approaches for capturing these
effects in the analysis of transfer functions are: (1) continuum modeling of the soil
and foundation system subject to input motions with a defined coherency function
(Computer Program SASSI2000, Lysmer et al., 1999; Ostadan et al., 2005); and (2)
application of a semi-empirical simplified model (Veletsos et al., 1997; Kim and
Stewart, 2003).
either vertically or at an angle v to the vertical, as shown in Figure 3-1a. The results
foundations on the surface of an elastic half-space to incoherent S-waves propagating
parameter, a, related to lagged coherency and wave inclination that scales the
essentially independent of foundation shape, but is strongly dependent upon a
frequency axis of the theoretical transfer function. For vertically propagating waves
this transfer function (adapted from Veletsos and Prasad, 1989) can be written as:
b0
(3-3)
where b0 4 a a0k and I0 and I1 are modified Bessel functions, zero and first
order, respectively. Equation 3-3 was developed for circular foundations. The
4 term adapts a 0k (for rectangles) to a0 defined for an equivalent radius that
preserves foundation area. For small and large values of the argument (2b02), the
Bessel function summation in Equation 3-3 can be written in terms of power series
and exponential functions, respectively (Watson, 1995); for routine application, these
approximations can be expressed as:
for b0 1
I 0 (2b02 ) I1 (2b02 )
b0 6 b08 b010
2 4
1 b b
b0 16b0
0 2 0
Note that the exponential terms in Equation 3-3 and Equation 3-4 cancel for b0 > 1.
The two functions in Equation 3-4 have a misfit of 0.0073 at b0 = 1, which is accurate
enough for practical purposes.
where Vs is a representative small-strain shear wave velocity for the soil beneath the
foundation, which can be calculated as described in Section 2.2.2.
semi-empirical approach near the upper and lower limits of a. Two key
amplitude variability. Figure 3-1b shows transfer function Hu calculated using the
observations from this figure are: (1) as a increases (indicating increasingly spatially
variable motions), Hu decreases significantly; and (2) for the range of a supported by
case history data, Hu from the semi-empirical procedure is much lower than from
wave passage models. This can be attributed to significant contributions of stochastic
phase and amplitude variability to base-slab averaging.
The data set considered by Kim and Stewart (2003) consists of buildings with mat
foundations, footing and grade beam foundations, and grade beam and friction pile
foundations, generally with base dimensions, BeA , in the range of 15 m to 40 m.
Although the Veletsos models strictly apply to rigid foundations, the semi-empirical
model applies to the more realistic foundation conditions present in the calibration
data set.
Errors could occur, however, when the model is applied to conditions beyond the
calibration data set. In particular, the effects of incoherence in the Veletsos models is
taken as proportional to wavelength, thus implying strong scaling with frequency and
distance. Array data indicate that distance scaling is much weaker than the frequency
scaling (Abrahamson et al., 1991; Ancheta et al., 2011), so the model would be
expected to over-predict the effects of incoherence (under-predict Hu) for very large
foundations. The opposite would be true for small foundations. Even within the
parameter range of the calibration data set, it should be recognized that the empirical
model fits the data in an average sense, and should not be expected to match any
particular observation.
If the base slab of a building is embedded below the ground surface (i.e., the structure
has a basement), foundation-level motions are further reduced as a result of ground
motion reduction with depth below the free surface. The available solutions apply to
rigid cylinders embedded in a uniform soil of finite or infinite thickness (half-space).
Analytical solutions by Kausel et al. (1978) and Day (1978) describe foundation
input motions at the base of embedded cylinders as a function of free-field surface
ground motion ug. When subjected to vertically propagating coherent shear waves,
embedded cylinders experience a reduction in base translational motion, relative to
ug, due to ground motion reductions with depth and wave scattering effects.
Rotations in the vertical plane are also introduced as a result of differential
displacements imposed upon the cylinders over their embedded depth. These transfer
functions can be adapted to rectangular foundation shapes as:
D
H u 0.45, 1.1 (3-6b)
Vs
L D D
H yy 0.26 1 cos ,
Vs
(3-6c)
ug Vs 2
D
H yy 0.26, (3-6d)
Vs 2
where D is the embedment depth, as shown in Figure 3-2a. Velocity, Vs, in this case
should be interpreted as the average effective profile velocity, Vs, avg, defined in
Section 2.2.2. The transfer functions for rotation in the xx direction have the same
B/ug.
form as Equation 3-6(c) and Equation 3-6(d), except that it expresses the rotation
Figure 3-2b plots the transfer functions for translation and rotation. The reduction of
translation is substantial at high frequencies, saturating at about 70% of fE, which is
the fundamental frequency of the soil column between the surface and depth, D.
With a high-frequency de-amplification level of 0.45, the embedment effect is often
more important than the base-slab averaging effect for building structures.
Flores-Berones and Whitman (1982), Fan et al. (1991), and Nikolaou et al. (2001)
describe the kinematic response of vertical piles and pile groups in elastic soil
subjected to vertically propagating coherent shear waves. Similar solutions for
inclined (coherent) waves have been presented by Barghouthi (1984), Mamoon and
Banerjee (1990), and Kaynia and Novak (1992). Because the incident motions
assumed in the development of these models were coherent, these models do not
adequately incorporate base-slab averaging effects. Consequently, model predictions
do not compare favorably to data (Kim and Stewart, 2003). In particular, transfer
function ordinates, Hu, are significantly over-predicted by the Fan et al. (1991)
models at high frequencies.
Kim and Stewart (2003) found that observed variations between foundation input and
free-field ground motions at building sites in California could be adequately
represented with the models for shallow foundations in Section 3.1 and Section 3.2.
The pile-supported buildings considered were generally founded on alluvial soils, and
the piles were not end-bearing. At building sites in Japan, however, Mikami et al.
(2006) found that model predictions for base-slab averaging (SASSI, and the semi-
empirical method in Section 3.1) overestimated Hu relative to data. The pile
foundations in the Japanese data set have relatively high flexural stiffness, and were
more nearly end-bearing, in comparison to the California data set. Development of
analytical solutions for the kinematic interaction problem for pile-supported
foundations of varying flexural rigidity subjected to realistic (incoherent) input
motion remains an important research need.
H u ( f ), f fL
S a FIM ( f )
(3-7a)
Sa( f )
At higher frequencies, the approximation from Equation 3-7 does not hold because
high frequency spectral ordinates and PGA are controlled by lower-frequency
components of the ground motion. At high frequencies, a conservative
approximation of spectral ordinates can generally be obtained from:
H u ( f L ), f fL
S a FIM ( f )
(3-7b)
Sa( f )
The limiting frequency, fL, depends on the frequency content of ug. For typical stiff
soil or rock ground motions having mean periods in the range of 0.2 sec to 0.5 sec, fL
has been found to be approximately 5 Hz, as documented in Appendix E of FEMA
440, Improvement of Nonlinear Static Seismic Analysis Procedures (FEMA, 2005).
Long-period ground motions resulting from near-fault directivity pulses or soft soil
effects, however, can have much lower limiting frequencies, fL, such that no
significant spectral ordinate reductions from kinematic interaction are realizable.
Further research is needed to identify factors controlling fL, and to develop more
reliable recommendations for engineering application.
Differences between transfer functions and ratios of response spectra are illustrated in
Figure 3-4 using the Rancho Cucamonga data presented in Figure 3-3. East-west
(E-W) response spectra are shown in Figure 3-4a. Reductions in foundation motions
relative to free-field motions are apparent for periods less than approximately 0.7 sec
(of frequencies greater than 1.4 Hz). Figure 3-4b shows the ratio of response spectra
(RRS) for the E-W direction along with: (1) the transfer function model, Hu, for both
base-slab averaging and embedment effects; and (2) the RRS model derived from Hu
using Equation 3-7. The RRS model captures the general trends of the data, although
there are significant period-to-period variations in the data (and even some RRS
ordinates greater than unity). The RRS data does not show the strong decrease in
spectral ratios with decreasing period that is evident in the Hu function, which
illustrates the saturation effect described above.
Acceleration histories representing the FIM can be modified from free-field motions
using the following procedure:
1. Calculate the Fourier transforms of ug.
For most practical situations, this procedure could be avoided by merely selecting
and modifying ground motions for compatibility with Sa-FIM in lieu of Sa, in which
case no further modification is needed.
(a)
(b)
However, trends in practice are tending toward increased use of SSI. This has been
driven principally by seismic retrofit projects in which SSI analysis is used to gain
better insight into structural performance and to improve accuracy in the analytical
simulation of important structural response quantities. This chapter describes the
implementation of SSI procedures in currently available engineering standards and
guidelines. Limitations of those procedures are also briefly discussed.
Figure 4-1 Schematic illustration of the shape of the design response spectrum in the
NEHRP Recommended Provisions.
0
(4-2)
The change in base shear is related to the change in seismic coefficient (or spectral
obtained from the design spectrum at an elongated period, T . The term 0.05 0
acceleration). The C s term in Equation 4-2 represents the seismic coefficient
0.4
from the fixed-base structural damping value of i = 0.05, to the flexible-base value
represents the reduction in spectral ordinate associated with a change in damping
of 0.
slope with respect to period. Coupled with the requirement that 0 must exceed i
It is important to note that the shape of the design spectrum is flat, or has a negative
Modification of design base shear for SSI effects in equivalent lateral force
procedures has a potentially significant shortcoming. There is no link between base
shear reduction factors intended to represent structural ductility (i.e., R factors) and
soil-structure interaction. Crouse (2001) noted that existing R factors may already
reflect the beneficial effects of soil-structure interaction, and modifying the base
shear to account for both SSI and ductility may be unconservative in some cases.
Accordingly, there is a need to revisit the definition of R factors with respect to SSI
effects, and define values that represent structural ductility effects alone.
1
T k kh 2
(4-3)
T k x k yy
however, ASCE/SEI 7-10 does not specify how lateral stiffness, kx, or rotational
stiffness, kyy, are to be evaluated. The Commentary to the NEHRP Recommended
Provisions provides some guidance related to circular foundations, but the
rectangular foundation models contained in Chapter 2, Section 2.2 are more useful
for practical engineering applications.
The shear modulus, G, used in conjunction with equations for static foundation
stiffness, must be reduced from the shear modulus at small strain levels, G0, to
account for large strain effects. ASCE/SEI 7-10 (and the NEHRP Recommended
Provisions) provides values for adjusting the shear modulus and shear wave velocity
for large strain levels based on site class and spectral response acceleration levels.
These values are provided in Table 2-1.
0 f i
T T
1
n
(4-4)
with fixed-base structural damping, i = 0.05, and exponent, n = 3 (for ideally viscous
the foundation damping factor, f, is not evaluated directly from controlling variables
material damping). In ASCE/SEI 7-10 (and the NEHRP Recommended Provisions),
function of structure aspect ratio, h/r. Note that r is an equivalent foundation radius,
which is calculated to match the foundation area for squat structures and the
The principal limitations of force-based procedures in ASCE/SEI 7-10 are: (1) use of
simplified spectra that can only result in a decrease in base shear as period lengthens;
A potentially important consideration associated with the use of the SSI procedures
in Chapter 19 of ASCE/SEI 7-10 is the value of the fixed-base fundamental period, T.
Chapter 12 of ASCE/SEI 7-10 contains approximate methods for evaluation of T, and
limiting values (i.e., CuTa), which bias the estimate of T to intentionally produce
conservative values of design base shear. In SSI procedures, T should be taken as the
best estimate value of period, without deliberate bias.
Figure 4-3 Schematic illustration of a pushover analysis and development of a pushover curve
for a structure with a flexible base.
Powell (2006) describes common ways by which the pushover curve is combined
with a design response spectrum to estimate the seismic displacement in a structure.
Three such methods are known as the Capacity Spectrum Method (ATC, 1996), the
Coefficient Method (FEMA, 1997; FEMA, 2000; and ASCE, 2007), and Equivalent
Linearization (FEMA, 2005). These methods are illustrated in Figure 4-4.
response spectrum for flexible-base damping ratios, 0, that are greater than the
field response spectrum for kinematic interaction effects; and (3) reduction of the
fixed-base structural damping ratio, i. The manner by which these components are
evaluated in displacement-based analysis procedures is described below. Slight
modifications to the notation contained in reference engineering standards and
guidelines have been made for consistency with the notation adopted in this report.
1 2 B e 0.3048
adapted from FEMA 440 as follows:
1.2
1 T
A
14100
1
RRS bsa (4-5a)
T fL
RRS bsa 1
1
2 B eA 0.3048 f L T
1
1.2
(4-5b)
14100 fL
2 Df L
RRS emb cos T
1
V sr
(4-6b)
fL
where Vsr is the strain-reduced shear wave velocity evaluated using the reduction
factors in Table 2-1. In Equations 4-5, the equivalent foundation dimension BeA is
expressed in units of meters. These equations are a curve-fit of the semi-empirical
base-slab averaging transfer function described in Equation 3-3 (in Chapter 3). A
shear wave velocity term does not appear in Equations 4-5 because the Vs terms
cancel in the expression for b0 in Equation 3-3. The resulting RRS curves for base-
slab averaging are shown in Figure 4-5.
Figure 4-5 Ratios of response spectra (uFIM/ug) for base slab averaging using
the semi-empirical formulation adopted in FEMA 440.
limiting period, fL, is taken as 5 Hz (0.2 sec). As of this writing, these equations are
in the process of being revised in the next version of ASCE/SEI 41 (ASCE, 2013) to
reflect the recommendations contained in Chapter 3 of this report.
damping ratio, f, which is then combined with the fixed-base structural damping
In FEMA 440, the objective of the damping analysis is to estimate the foundation
ratio, i, to estimate 0 using Equation 4-4 (with n = 3). The principal challenge is to
extract f from the results of the pushover analysis of the structure in both its fixed-
foundation flexibility can significantly reduce radiation damping (yy) from rotational
base and flexible-base condition. As described in Chapter 2, Section 2.2.3,
Table 2-2 (the dynamic stiffness modifier, x, is assumed as unity). The effective
Assuming shaking in the x-direction, stiffness Kx is then evaluated using equations in
K yy
K *fixed h 2
T
(4-7)
2
K *fixed
T
1
Kx
*
where K fixed is the equivalent fixed-base stiffness of the structure evaluated from:
2
M
2
T
*
K fixed (4-8)
Note that dynamic stiffness modifier, yy, is also taken as unity. The value of Kyy
estimated from Equation 4-7 reflects the stiffness of the foundation structural
elements as implemented in the pushover analysis, so no assumptions of foundation
rigidity are required.
The next step is to reduce the period lengthening ratio from the small-displacement
periods). Taking as the expected ductility demand for the system (including
condition to the large-displacement (i.e., post-yield) condition (with elongated
structure and soil effects), the effective period lengthening in the post-yield state is
computed as:
T 1 T 2
0.5
1 1
T eff T
(4-9)
Most of the resources listed at the beginning of this chapter (e.g., ATC-40, FEMA
440, FEMA P-750, ASCE/SEI 41-06, and ASCE/SEI 7-10) are silent on the
implementation of SSI effects in response history analyses. Similar to ASCE/SEI
7-10 (Section 12.13), they permit the use of soil springs in principal, but offer no
specific guidance on how the springs should be selected or utilized in a response
history analysis.
Subterranean levels are modeled in both the SLE and MCE analyses, including the
mass, stiffness, and structural capacities of structural elements such as walls,
columns, and slabs. Response history analysis for the SLE (Figure 4-6b), is
performed with a relatively simple model that omits the surrounding soil and does not
include soil springs. Response history analysis for the MCE (Figure 4-6c), is
performed with springs and dashpots representing soil-foundation interaction along
basement walls and below the base slab. In this case, ground motions are applied to a
rigid “bathtub” surrounding the subterranean portions of the structure. In both the
SLE and MCE analyses, the motion applied at the base of the model can be either the
free-field motion (ug) or the foundation input motion (uFIM). These recommendations
are derived largely from the recommendations of Naeim et al. (2008).
Procedures for calculating spring stiffnesses and capacities are not specified in the
PEER Guidelines, but can be taken as those from Chapter 2 of this report. Similarly,
foundation input motions can be modified from free-field motions using the
procedures in Chapter 3 of this report.
Guidance for the seismic analysis of nuclear safety-related facilities in the United
States is provided in ASCE 4-98, Seismic Analysis of Safety-Related Nuclear
Both the direct analysis and substructure approaches for SSI described in Chapter 1
are permitted under ASCE 4-98. Provisions related to response history analyses as
well as equivalent lateral force-based analyses are included. Referring to Figure 1-2
(in Chapter 1), consideration of kinematic interaction effects, foundation flexibility,
and damping are included. ASCE 4-98 does not consider base slab averaging effects,
but it does consider embedment effects. The foundation input motion adjusted for
embedment effects is referred to as the Foundation Input Response Spectra (FIRS),
which is computed for the free-field conditions at the foundation level using wave
propagation analysis (de-convolution). The procedure in the upcoming version of
ASCE 4 (ASCE, in preparation) will also allow for base-slab averaging through
response history analysis, in which the input motion is specified with a defined
coherency function, as described in Section 3.1, using software such as SASSI
(Ostadan et al., 2005).
Testing has played a relatively minor role in the research and evolution of soil-
structure interaction. For the most part, research has been dominated by numerical
analysis exercises. The models developed from numerical and theoretical studies,
however, apply for idealized conditions, so testing and seismic monitoring play a
vital role in guiding the implementation of idealized models in practice.
This chapter identifies and reviews experiments and case studies available in the
literature. Information is separated into field-scale and laboratory-scale tests, which
typically have different objectives and different applications. Tests involving
shallow foundations and dynamic loading are emphasized. Results from tests
identified in this chapter, where applicable, have been interpreted for use in SSI
modeling in other chapters of this report.
Seismic monitoring and field testing of structures is generally performed with one of
two objectives: (1) evaluation of system properties such as the fundamental mode
period and damping of an SSI system; or (2) estimation of foundation impedance
ordinates representing the stiffness and damping characteristics of soil-foundation
interaction.
5.1.1 System Studies
In system studies, seismic data from well-instrumented buildings are used in system
identification analyses to estimate modal vibration parameters. To evaluate SSI
effects, input-output pairs must be selected that isolate the system behavior
associated with the structure alone (fixed-base properties) and the full system
(flexible-base properties), as described in Stewart and Fenves (1998). Results can
then be compared to period lengthening and foundation damping models of the type
described in Section 2.1. Stewart et al. (1999a; 1999b), for example, describe the
results of such analyses and lessons learned regarding the practical application of
impedance functions for rigid circular foundations.
System studies can also be undertaken using forced-vibration tests (e.g., Yu et al.,
2005). The input-output pairs that should be used for evaluation of fixed- and
flexible-base properties are given by Tileylioglu (2008). No major study applying
these techniques to evaluate period lengthening and foundation damping has been
completed to date.
Table 5-1 Summary of Field Forced-Vibration Tests Used to Evaluate Foundation Impedance
Functions
Excitation Results
Lin and
shaker on 1.3 (D=1.5); modal freq.
33m 0 to 1.5 17.5 305 7 to 70 kx, cx, kyy, cyy Jennings
ground 1.5 (D = 0) only
(1984)
Luco et al.
NS: 2.16; shaker on NS: 0.8 to 2.5; NS: 1.06; NS: 0.8-2.5;
2525m 4 to 5.5 300 kx, cx, kyy, cyy (1988); Wong
EW: 1.26 roof EW: 0.8 to 1.75 EW: 1.1 EW: 0.8-1.75
et al. (1988)
1.33m; shaker on Crouse et al.
0 n/a 120; 75 10 to 60 n/a kx, cx, kyy, cyy 0 to 60
1.21.1m fndn. (1990)
cyykx-yy
shaker on kx, cx, kyy,
de Barros and
diam. = 10.8m 5.2 9.37 300 roof and 2 to 20 2 5 to 14
Luco (1995)
fndn. cx‐yy, kz, cz
shaker on Tileylioglu et
4.14.1m 0 6.0, 13 198 5 to 15 1.15, 1.29 kx, cx, kyy, cyy 5 to 15
roof al. (2011)
Notes: (1) Foundation embedment depth; (2) Fundamental mode, fixed-base frequency; (3) Vs=Shear wave velocity of soil;
(4) Frequency range; and (5) Period lengthening.
Symbols: D=embedment depth; n/a=not available; NS=North-South building axis; EW=East-West building axis;
diam=diameter; and fndn=foundation.
The first field investigations of foundation impedance provided results over a limited
range of frequencies (Lin and Jennings, 1984; Luco et al., 1988; Wong et al., 1988)
or for very small structures representative of strong motion instrument huts (Crouse
et al., 1990). More recently, de Barros and Luco (1995) tested a relatively large
model structure (of a nuclear reactor) and provided impedance ordinates over a
relatively wide frequency range (approximately 4 Hz to 20 Hz). Figure 5-1 shows
impedance ordinates evaluated by de Barros and Luco. Results are shown in non-
normalized form due to uncertainty in the shear modulus of the foundation soils, and
illustrate the noisy character of the data, especially at frequencies under 4 Hz or
greater than 14 Hz. Also shown in Figure 5-1 are three model predictions for
Figure 5-1 Non-normalized impedance values from experimental data compared with theoretical
predictions for a nuclear containment structure at Hualien, Taiwan for translational (top)
and rotational (bottom) modes (adapted from de Barros and Luco, 1995).
Two practical difficulties associated with field testing for impedance ordinates and
comparison to model predictions have been encountered in the previous work. First,
limited resolution of the data acquisition system with respect to analogue-to-digital
signal conversion and time-stamping contribute significantly to noise in the results.
Most previous studies have not formally evaluated noise effects, which can lead to
spurious results (e.g., impedance ordinates in Figure 5-1 for frequencies outside the
4 Hz to 14 Hz range). Second, shear wave velocity profiles have generally been
established using downhole or suspension logging methods in the free-field.
Use of downhole or suspension logging methods can have limited resolution very
near the ground surface (e.g., Andrus et al., 2004). Because the soil materials
immediately below the foundation exert the greatest influence on foundation
stiffness, this introduces uncertainty in the selection of an appropriate value of Vs for
The final study listed in Table 5-1 overcame many of the practical difficulties
encountered in previous work through the use of high-fidelity, modern sensors and
data acquisition equipment available through the NEES@UCSB equipment site
(http://nees.ucsb.edu/). Seismic velocities were measured immediately adjacent to
the foundation so that they incorporated overburden effects. Sample results are
shown in Figure 5-2, and the data are observed to be numerically stable across the
range of tested frequencies. Stronger damping was evident in the translational versus
the rotational vibration modes, and radiation damping was a significant contributor in
each case.
(a) (b)
Figure 5-2 Normalized impedance values from experimental data compared with
theoretical predictions for the Garner Valley site for: (a) translational; and
(b) rotational modes (adapted from Tileylioglu et al., 2011). Velocity, Vsm,
and modulus, Gm, denote median values from test data.
Field tests have been supported by nuclear regulatory agencies in the United States
and abroad to provide data that can be used to validate analysis procedures, such as
CLASSI (Wong, 1979) and SASSI (Lysmer et al, 1999), commonly used for nuclear
reactor structures. In one case (Hualien site in Taiwan), these experiments can be
used to infer impedance functions. Such cases were presented above and are not
repeated here (de Barros and Luco, 1995). Following is a discussion of three similar
experiments:
Lotung, Taiwan. This experiment utilized a ¼-scale containment model
constructed in Lotung, Taiwan. The objective of the study was to validate SSI
analysis methodologies. The Lotung model was a cylindrical structure with a
total height of 15 m and embedment depth of 4.5 m. The free-field, structure,
and internal components were instrumented to record motions and pressures at
the soil-foundation interface. Earthquakes recorded in 1986 were used for
subsequent analyses. The model reactor was installed in a soft soil site with
measured shear wave velocities, shear moduli, and damping curves with shear
strain. The data were used in a blind prediction exercise (Ostadan et al., 1991).
Because the data consisted of earthquake recordings and not controlled
experiments, it was not practical to directly identify impedance functions across a
significant frequency range, as described in Section 5.1.2.
However, laboratory tests are limited in their ability to reproduce certain field
conditions (e.g., Novak, 1987). For example, the finite size of a laboratory test
container precludes radiation damping of waves with quarter-wavelengths on the
order of the container dimension. This is illustrated in Figure 5-3, which shows
theoretical model-based damping ratios increasing with foundation size, whereas the
experimental data show essentially constant damping ratio at the hysteretic (material)
damping level. The difference is due to radiation damping that is present in the
theoretical model, but not in the experiments, due to the small size of the laboratory
container.
h=45.7cm
square/rec.
applied static vertical, Dobry et al.
=0.33;
B=5.6 to moist sand; level ground
force at top of torsion, (1986)
=121 pcf;
13.3 cm concrete D=152.4
0 1 block and coupled Report Erden (1974)
s=0.025
L=B-6B block 2B=304.8
suddenly swaying- Stokoe and
circle 2L=304.8
released rocking Erden (1985)
radius=14.9cm
h=45.9 to
46.7cm applied static Gazetas and
=0.33;
moist sand; level ground vertical,
square/rec. force at top of Stokoe (1991)
=121 pcf;
concrete ~7.5 to D=152.4 coupled
B=5.6-13.3 cm 1 block and Report Erden (1974)
s=0.025
block 18.8 2B=304.8 swaying-
L=B-6B suddenly Stokoe and
2L=304.8 rocking
circle released Erden (1985)
radius=14.9cm
vertical,
aluminum level ground base/sidewall
h=2.8cm Dry, No. 120 torsion,
block D=35.6 shear and Downloadable Gadre and
B=1.9cm 3.8 30 Nevada sand coupled
1cm below 2B=61.0 passive/active data; report Dobry (1998)
L=1.9cm Dr=75% swaying-
surface 2L=91.4 force
rocking
Nevada sand applied slow
level ground
h=50.8cm double Dr=60 to 80% sinusoidal Rosebrook and
D=53.0 Downloadable
B=1.73cm aluminum 0 to 1.5 20 and cyclic force rocking Kutter (2001a,
2B=90.0 data; report
L=6.68cm shear walls bay mud and dynamic 2001b, 2001c)
2L=175.0
Cu=100 kPa base shaking
applied slow
h=48.6cm level ground
sinusoidal
B=2.5 to single wall Nevada sand D=53.0 Downloadable Chang et al.
2.5 20 cyclic force rocking
1.73cm with frame Dr=80% 2B=90.0 data; report (2007)
and dynamic
L=2.5 to 7.1cm 2L=175.0
base shaking
applied slow
h=276.4cm level ground
sinusoidal
B=13.5 to Nevada sand D=53.0 Downloadable
bridge pier 8.6 20 cyclic force rocking Ugalde (2007)
17.75cm Dr=80% 2B=90.0 data; report
and dynamic
L=13.5 to 2L=175.0
base shaking
17.75cm
saturated; applied vertical,
level ground
steel base dense sand sinusoidal torsion, Ghosh and
D=22.0
radius=1.5cm plate with 0 50 (Dr=85%) over cyclic force coupled Report Madabhushi
2B=23.5
dual dome loose sand (freq=50Hz swaying- (2007)
2L=56.0
(Dr=45%) for 500ms) rocking
single and Nevada sand applied slow
h=50.75cm level ground
double rigid Dr=60 to 80% sinusoidal
B=1.63cm D=53.0 Downloadable Gajan and
steel or 0 to 7 20 and cyclic force rocking
L=6.75 to 2B=90.0 data; report Kutter (2008)
aluminum bay mud and dynamic
7.0cm 2L=175.0
shear wall Cu=100 kPa base shaking
2L=total length; 2B=total width; h=height; g-level=gravity load multiplier applied while testing; D=embedment depth; =Poisson's ratio;
=unit weight of soil; s=soil damping ratio; Dr = relative density of sand; and Cu = undrained shear strength of clay.
Symbols:
Substantial research has been available for many years, yet there is relatively limited
implementation of soil-structure interaction in engineering practice. There is often a
lag between the state of knowledge and the state of practice, and this has proven to be
particularly true in the case of soil-structure interaction. Even among relatively
experienced practitioners, it is clear that there is room for improvement in the
technical approaches used in SSI modeling, as well as the manner in which structural
engineers, geotechnical engineers, engineering seismologists, and other design and
construction professionals interact with each other on projects involving soil-
structure interaction.
This chapter summarizes general observations on the state of the practice with regard
to modeling soil-structure interaction effects on building structures, and makes
specific recommendations for the process of communication and collaboration
between design professionals. As a result, the intended audience for this chapter is
wider than it is for other chapters in this report, and includes owners, project
managers, architects, construction managers, and others who are involved in
managing, designing, and constructing the built environment.
6.1 Overview
The state of the practice with regard to SSI was discussed informally and anecdotally
with selected structural engineers, geotechnical engineers, and members of the
project team. Design professionals engaged in these discussions represented small
and large firms with different backgrounds, levels of experience, and geographic
locations. Collectively their experience covered a wide range of building sizes,
occupancy types, and structural materials. Comments attributed to geotechnical
engineers may also be related to engineering seismologists that are providing
recommendations as part of the geotechnical engineering scope of services, either as
a member of the geotechnical engineering firm, as a subconsultant, or as an
independent consultant on projects.
6.2 Observations
This section synthesizes observations of the state of practice related to general issues,
collaboration between design professionals, information needed and shared,
understanding of SSI principles, and implementation in terms of SSI analyses, soil
springs, and common approaches to design problems.
On most projects, structural engineers are engaged by the architect, but geotechnical
engineers are engaged by the building owner. This arrangement is primarily the
result of a perceived increase in liability for geotechnical engineering, and the
reluctance of architects, and their professional liability insurers, to engage
geotechnical engineers as subconsultants. As a result, geotechnical engineers are not
directly managed by the architect as lead design professional, and geotechnical
Traditionally, base-slab averaging and embedment effects have not been considered,
and this information is not regularly requested. With the publication of FEMA 440,
Improvement of Nonlinear Static Seismic Analysis Procedures (FEMA, 2005) and
other subsequent engineering resources, structural engineers are beginning to take
advantage of these effects on new projects, and are beginning to request related
design information.
Although there are consistent needs on many projects, few structural engineers report
the availability of checklists, or other standardized lists, to assist in organizing and
requesting the geotechnical information necessary for foundation design or SSI
modeling. A sample checklist of geotechnical information needed by the structural
engineer would be considered a useful aid.
Geotechnical engineers can provide better recommendations when they have more
detailed information on which to base them. In addition to information from soil
Understanding of SSI principles varies widely across both the structural and
geotechnical engineering disciplines. Most structural engineers can appreciate that
SSI effects are more pronounced in soft soils, and many are aware that foundation
input motions can differ from free-field ground motions. Although SSI effects are
known to be significant on stiff, squat, shorter-period buildings, many practitioners
mistakenly believe that SSI would have a larger effect on taller, flexible, longer-
period buildings.
Analysis procedures used to address seismic loading include equivalent lateral force
procedures, modal response spectrum analysis procedures, nonlinear static
(pushover) procedures, and nonlinear response history procedures. Although there is
increasing use of nonlinear procedures in design, response history analyses are
performed by a minority of structural engineers, and among that minority, response
history analyses are performed on a relatively small percentage of projects. Currently
available codes, standards, and guidelines provide SSI provisions for force-based
procedures (e.g., equivalent lateral force and response spectrum analysis) and
Most software packages commonly used in building design have the ability to
include uniaxial springs representing support flexibility. However, some are not able
to model compression-only properties that might be appropriate for non-embedded
(i.e., surface) foundations. More sophisticated programs can vary the compression
and tension properties of the springs, and permit bilinear or tri-linear force-
displacement relationships. Multi-support excitation is difficult to implement and is
not widely available in commonly used structural engineering software, though it is
available in programs such as OpenSees (McKenna, 1997).
Vertical springs are almost always used in the design of mat foundations in order to
properly understand the distribution of design forces in the mat. This is typically
done with a single modulus of subgrade reaction representing long-term settlement
properties of the supporting soil. Use of springs under gravity load-carrying grade
beams is also relatively frequent. Springs are occasionally placed beneath spread
footings, strip footings, grade beams, piles, and piers. To a lesser extent, they are
used to model the combination of soil and foundation vertical flexibility (e.g.,
individual piles or pile groups).
In comparison, horizontal springs are much less frequently used. In rare cases, they
are used on piers and piles to understand the distribution of lateral forces to the
foundation elements, particularly if the foundation consists of different systems with
different stiffnesses. Horizontal springs are infrequently used to represent the passive
pressure developed on the sides of pile caps, grade beams, or retaining walls.
Often the modulus of subgrade reaction, given for long-term gravity loading, is used
as the basis for dynamic loading conditions. This can lead to significant
underestimation of soil stiffness during dynamic loading, and significant
overestimation of the amount of displacement and rotation caused by the deflecting
soil. The cause for this appears to be a lack of understanding about the differences
between long-term and short-term loading effects on soil, a lack of clarity on what is
needed by the structural engineer, and a lack of understanding of the limits of the
information the geotechnical engineer is providing.
When foundation springs are used, practice varies on the extent to which uncertainty
is used to bound expected soil properties. Resources like FEMA 356 and ASCE/SEI
41-06 specify a range of two times the expected values to one-half of the expected
values, but this is dependent on the variability in site conditions and scatter of
material properties. Some structural engineers use a simple, conservative upper
bound assuming a fixed base (or an essentially infinitely rigid spring). Others modify
the spring values provided by the geotechnical engineer in some way. The best result
is obtained when the geotechnical engineer provides upper- and lower-bounds around
expected values. This usually narrows the bounds recommended by FEMA 356 and
ASCE/SEI 41-06, but is dependent on the variability in site conditions and scatter
present in measured soil properties.
For deep foundations such as piles or piers, different stiffnesses are occasionally
provided for tension and for compression. Lateral springs are typically developed
using common commercial software programs that provide the displacement, shear,
and moment profile considering the variation in soil properties over the depth of the
element.
Figure 6-3 illustrates a building with a basement that does not have retained earth on
one or more sides. In this situation, the loading caused by the soil pressure on one
side of the building is typically analyzed separately and added to the building inertial
loading by linear superposition. Structural engineers report the use of Models 2A,
2B, 3A, and 3B, as shown in Figure 6-2, for modeling partially embedded buildings.
Figure 6-4 shows a building without a basement. Structural engineers report the use
of different modeling approaches for buildings without a basement, depending on the
type of foundation system:
6.3 Recommendations
Better clarity on the part of structural engineers with regard to the modeling approach
for the soil and foundation will help geotechnical engineers provide better
recommendations. Having a defined number of meetings between the geotechnical
and structural engineers will improve communication and help ensure proper
development and implementation of foundation recommendations.
Owners can level the playing field in the competitive market place by clearly
requesting these additional scope items in the proposal. Owners can also coordinate
the timing of consultant selection and engagement so that geotechnical and structural
engineers are on a concurrent project delivery schedule. Some projects might benefit
from a foundation design charrette involving a foundation contractor to discuss
project specific installation issues that can help in refining and improving the design.
Soil springs are often not properly, or adequately, implemented in practice. In some
cases, the modulus of subgrade reaction developed for long-term settlement is being
used for dynamic loading situations. In other cases, a single value or single force-
displacement relationship is being used for the soil properties. Figure 6-5 provides a
sample format for presentation of soil spring data.
25
qc(UB) = 22.5 ksf (Short‐term: Upper Bound)
kUB=1920 kcf
qc(UB) = 20 ksf (Short‐term: Expected)
20
k=960 kcf qc(UB) = 18 ksf (Short‐term: Lower Bound)
Bearing Pressure (ksf)
15
kLB=480 kcf
qall = 10 ksf (Short‐term D+L+EQ at ASD)
10
qall = 7.5 ksf (Long‐term D+L at ASD)
qall = 5 ksf (Long‐term D at ASD)
5
kLT= 180 kcf
0
0 0.5 1 1.5 2
Settlement (inches)
The sample format includes both the relationship for long-term settlement, which
needs only a lower-bound, and recommendations for the expected (or target)
properties under dynamic loading, including both an upper-bound and lower-bound.
Case studies illustrating the level of effort, relative benefits, and potential drawbacks
associated with simplified SSI modeling approaches are provided in example
applications presented in Chapter 7 for a variety of different building configurations.
Methods for developing soil springs and dashpots described in this report are
implemented on case-study buildings using several models with different
idealizations at the soil-foundation interface. Examples are used to illustrate the
corresponding level of effort for different modeling approaches. Overall results and
selected structural response quantities from the models are compared and contrasted.
Analyses are performed in the context of a substructure approach to modeling SSI
effects.
The overall approach for development of example applications involved: (1) the
selection of suitable instrumented buildings that have experienced (and recorded)
earthquake ground motions; (2) development of baseline models for comprehensive
substructure-based analysis of seismic response; (3) calibration of baseline models to
approximately match the recorded response of the buildings through variation of
structural parameters; and (4) systematic variation in the idealization of the soil-
foundation interface to evaluate the impact of different modeling approaches on the
predicted response of the buildings.
Work included a review of the results from similar studies performed by others and
published in the literature. Information from these studies was used to guide the
example applications presented herein, and corroborate the resulting observations.
Prior studies reported by Naeim et al. (2008) and Tileylioglu et al. (2010) utilized
typical structural engineering software packages, such as ETABS, Extended Three
Dimensional Analysis of Building Systems (Computers and Structures, Inc.), and
SAP2000, Integrated Software for Structural Analysis and Design (Computers and
Analyses initially employed the full substructure modeling approach, designated the
Baseline Model (or MB) as described in Chapter 6, except that kinematic base
rocking was applied to the base of vertical foundation springs along with depth-
variable ground motions. This results in double-counting of kinematic rotation
effects. Fortunately, these effects were minor for the structures investigated, and
double-counting did not have a significant impact on results.
Springs were elastic, with no compression capacity limit, and zero tension capacity.
Kinematic effects were evaluated in a manner similar to that described in Chapter 3
(except that kinematic base rocking was applied). Significant difficulties were
encountered in implementing multi-support excitation using typical engineering
production software packages. Results were considered reliable for displacement
response, but not reliable for forces. In general, however, good matches between
computed and observed responses were reported using the baseline modeling (MB)
approach.
Selected elements of the full substructure modeling approach were then omitted from
the models to investigate their impact on the computed response. The following
factors did not have a significant impact on the results: (1) consideration of multi-
support excitation along basement walls; and (2) application of a zero-tension
condition in the foundation springs. Consideration of kinematic interaction effects
had a significant impact on the distribution of interstory drifts, particularly below
grade. Consideration of foundation springs had a significant impact on building
vibration periods and distribution of interstory drifts.
Two modeling approaches commonly used in practice were shown to provide poor
results: (1) fixing the structure at the ground surface, truncating the embedded
portions of the structure, and applying the free-field translation as the input motion
(Chapter 6, Model 1); and (2) modeling the embedded base of the structure, using
horizontal and vertical springs with ends fixed against translation, and applying free-
field motions as input at the base slab level (Chapter 6, Model 3).
Two case-study buildings were selected from CSMIP, consisting of: (1) a 13-story,
reinforced concrete moment frame structure (designed in 1964; retrofitted in 1994),
with two basement levels, located in Sherman Oaks, California; and (2) a 10-story,
reinforced concrete shear wall core and perimeter moment frame structure (designed
in 1970), without any basement levels, located in Walnut Creek, California. The
exterior elevations of the buildings are shown in Figure 7-1.
(a) (b)
Figure 7-1 Exterior elevations of two case-study buildings: (a) 13-story Sherman
Oaks building; and (b) 10-story Walnut Creek building (CESMD,
2011, with permission).
No single building was found to meet all of the selection criteria. The Sherman Oaks
building satisfied criteria (1) and (3) above, while the Walnut Creek building satisfied
criteria (2), (3), and (4) above.
Once baseline models have been successfully calibrated to match recorded data, the
soil-foundation interface is idealized using one of the following simplified modeling
approaches from Chapter 6:
Model 1. In Model 1, only the above-ground portion of the structure is modeled,
the base is fixed at the ground surface, and the free-field ground motion, ug, is
applied at the base of the model.
Model 2. In Model 2, the above-ground and subterranean portions of the
structure are modeled, the based is fixed at the bottom of the embedded portion
of the structure, the soil surrounding the embedded portion is ignored (i.e., no
horizontal foundation springs are used), and the free-field ground motion, ug, is
applied at the base of the model.
Model 3. In Model 3, the above-ground and subterranean portions of the
structure are modeled, horizontal and vertical soil springs are included, the far
end of each spring is fixed against translation, and the free-field ground motion,
ug, is applied at the base slab while horizontal springs at other levels remain
fixed. In the example applications, Model 3 is investigated using response
history analysis, although in practice, it is typically used in nonlinear static
(pushover) analysis applications.
Model 4. In Model 4, the above-ground and subterranean portions of the
structure are modeled, horizontal and vertical soil springs are included, the far
end of each spring is fixed against a rigid “bathtub” element, and the horizontal
foundation input motion, uFIM, or free-field motion, ug, is applied to the rigid
element. The key difference between Model 4 and Model MB is the manner in
which the seismic demand is applied. In Model 4, the effect of kinematic loading
on basement walls associated with depth-variable displacement histories is
neglected. Further, the use of free-field motions in lieu of foundation input
motions neglects kinematic interaction altogether by replacing the recorded
Table 7-1 summarizes the properties of the modeling approaches considered in the
example applications. Variations in Model MB (denoted MB.1 and MB.2) are
considered as follows:
In Model MB.1, the embedded portion of the building is assumed to be rigid.
The specification of seismic demand is not modified. The objective of this model
is to investigate the effects of flexibility in the subterranean structural elements.
In Model MB.2, there is a change in the way that soil flexibility is modeled. In
this model, springs are not allowed to develop tension to investigate the effects of
nonlinearities in the foundation springs.
Table 7-1 Summary of Modeling Approaches Considered in Example Applications
Model Variations
Parameters No. 1 No. 2 No. 3 No. 4 MB MB.1 MB.2
Structural Foundation
N/A Flexible Flexible Flexible Flexible Rigid Flexible
Elements
4a: uFIM
Input Motion ug ug ug uFIM uFIM uFIM
4b: ug
Depth-variable Ground
No No No No Yes Yes Yes
Motion
Not
Spring Tension n/a n/a Permitted Permitted Permitted Permitted
Permitted
All modeling schemes, except Model MB.2, permit tension to develop in the soil
springs. The rationale for this approach is that, prior to an earthquake, actual
basement wall-to-soil contact pressures range between ‘at rest’ and ‘active’ earth
pressures. Springs have an initial condition of zero force. Since earthquake shaking
will impose alternating cycles of increased and decreased pressures relative to the
initial state, wall-to-soil contact can be represented by a spring that develops tension,
provided that the level of deformation does not lead to gap formation. Given the
modest levels of demand imposed on the Sherman Oaks and Walnut Creek buildings,
this was expected to be a reasonable assumption, and this assumption was tested
using Model MB.2 on the Sherman Oaks building.
The Sherman Oaks building is a 13-story structure above grade with two basement
levels below grade. The seismic force-resisting system consists of reinforced
concrete moment-resisting frames that extend from the roof to the foundation,
supplemented by perimeter concrete shear walls in the subterranean levels. The
gravity system consists of cast-in-place concrete slabs supported on concrete beams
Construction drawings of the Sherman Oaks building were made available for
inspection through the auspices of the California Strong Motion Instrumentation
Program. The building measures 50 m (164 ft) tall from the ground surface to the
roof. The plan dimensions of the superstructure are 21.9 m (72 ft) wide by 57.6 m
(189 ft) long, although it widens at the first floor to match the foundation dimensions.
The height of the first story is 7.0 m (23 ft), and other above-grade stories are 3.6 m
(11.75 ft). The basement levels are embedded approximately 6.2 m (20.5 ft) below
the ground surface. The height of the first basement level is 3.5 m (11.5 ft), and the
height of the second basement level is 2.7 m (9 ft). The characteristics of the
foundation are described in Section 7.2.3 below.
A plan view of the Sherman Oaks site is shown in Figure 7-3. The ground surface
elevation is approximately 216 m (709 ft) above mean sea level. Boring information
and geophysical logs were obtained from geotechnical investigations conducted in
the vicinity of the site (LeRoy, Crandall and Associates, 1978 and 1982).
Figure 7-3 Plan view of the Sherman Oaks site showing locations of borehole
and geophysical logs used for site characterization (adapted from
LeRoy, Crandall and Associates, 1978 and 1982).
Figure 7-4 Subsurface characteristics of the Sherman Oaks site: (a) shear wave velocity
profile; and (b) material profile (adapted from LeRoy, Crandall and Associates,
1978 and 1982).
The average moist unit weight was taken to be 20 kN/m3 (124 lb/ft3). To estimate
strength properties of the foundation soils, direct shear test results from samples were
utilized. In the shallow soils, where most of the soil-foundation load transfer will
occur, the soil is unsaturated, so drained shear strengths were used. For the range of
surcharge pressures over the soil profile extending to the bottom of piles at a depth of
parameters inferred from available data are: (1) c′ = 20.1 kN/m2, and ′ 30°, from
approximately 16.2 m (53 ft) below grade, the soil Mohr Coulomb strength
0.0 m to 6.1 m; and (2) c′ = 12.9 kN/m2, and ′ = 38°, from 6.1 m to 16.2 m.
The Sherman Oaks building was instrumented in 1977 by the California Strong
Motion Instrumentation Program and designated CSMIP Station No. 24322.
Figure 7-5 Sherman Oaks building (CSMIP Station No. 24322) instrument
locations (CESMD, 2011, with permission).
Since 1977, six earthquake events have been recorded and processed by CSMIP at
this station. The events, along with peak accelerations for the second basement level,
ground floor, and roof are presented in Table 7-2. Note that in the 1987 Whittier and
1994 Northridge events, the peak accelerations followed an unusual pattern in which
the largest recorded motions occurred at the ground floor (i.e., the recorded motions
at the foundation and roof levels were smaller). As a result, subsequent calculations
are based on peak accelerations measured at the foundation level.
Horizontal translations recorded at the second basement (foundation) level were used
as foundation input motions, uFIM. The foundation input motion, uFIM, is the modified
free-field, ug, response due to base-slab averaging and embedment effects. Typically,
ug is known, and uFIM must be calculated based on transfer functions. In this case,
uFIM was measured and ug was inferred using transfer functions to remove the base-
Figure 7-6 Response spectra for recorded foundation input motion, uFIM, and inferred free-
field motion, ug, at the Sherman Oaks building, 1994 Northridge earthquake.
An equivalent linear ground response analysis was performed using DEEPSOIL v4.0
(University of Illinois, 2011) to estimate the depth-variable ground response adjacent
to the embedded portion of the structure, taking the recording of uFIM as input at
6.7 m with an elastic half-space. Input motions in ground response analyses can be
specified as “outcropping” or “within,” the former requiring an elastic half-space, and
the latter requiring a rigid base (Kwok et al., 2007). For the Sherman Oaks building,
the motion was recorded “within” the profile, but on a large foundation slab that
could be interpreted as representing an outcropping condition. Analyses were
performed for both conditions, with “outcropping” results appearing to be more
realistic.
Figure 7-7 Computed variation of peak ground acceleration (PGA) and peak
ground velocity (PGV) with depth at the Sherman Oaks site, 1994
Northridge earthquake.
The foundation is rectangular in plan, measuring 36.3 m (119 ft) wide by 57.6 m
(189 ft) long. The foundation consists of bored pile groups of varying configurations,
interconnected by pile caps and grade beams. The cast-in-place (bored) concrete
piles measure 51 cm (20 in) in diameter, and extend to a depth of approximately
9.9 m (32.5 ft) below the lowest basement level (i.e., approximately 16.2 m below
grade). Typical pile spacing varies between 0.9 m (3 ft) to 1.2 m (4 ft), center to
center. Concrete basement walls surround all embedded levels below grade. A
foundation plan developed based on information contained in construction drawings
provided by CSMIP is shown in Figure 7-8.
For the Sherman Oaks building, resistance is provided by both shallow and deep
foundation elements. It is assumed that the shallow foundation elements (principally
grade beams) remain in contact with the soil, which is justified by the sandy nature of
the soils at the site and the presence of friction piles rather than end-bearing piles.
Horizontal translation (x and y), base spring BL 17.8 6.2 to 24.0 296.2
Notes: (1) Calculated using overall foundation half-width, B=11 m; half-length, L=28.8 m.
In Table 7-3, the maximum depth considered is greater than zp by the amount of the
foundation embedment, D (i.e., the depth range extends to D + zp). For the specific
case of the horizontal base spring, the depth interval begins at the base of the
foundation.
density, s, of 20 kN/m3/g, and values of average effective profile velocity, Vs, avg,
In Table 7-4, shear modulus, G, was evaluated from Equation 2-9, using a soil mass
s
Horizontal, overall
k x ,total x K x , sur x cx,total 2k x ,total emb
k , cx ,total
x-direction
65.2 6.48e6 1.45 0.13 1 0.076 0.114
x ,total
k y ,total y K y , sur y s
Horizontal, overall
cy ,total 2k y ,total emb
k , c y ,total
y-direction 65.2 7.04e6 1.45 0.13 1 0.070 0.106
y ,total
s
Horizontal, base spring
k x ,base x K x , sur cx,base 2k x ,base sur
k , cx ,base
x-direction 88.6 8.81e6 --- 0.11 1 0.065 ---
x ,base
k y ,base y K y , sur s
Horizontal, base spring
cy ,base 2k y ,base sur
k , c y ,base
y-direction 88.6 9.56e6 --- 0.11 1 0.060 ---
y ,base
s
k z z K z , sur z cz 2kz emb
kz , cz
Vertical, z
65.2 8.55e6 1.22 0.13 1 --- 0.132
s
k xx xx K xx , sur xx cxx 2kxx emb
kxx , cxx
Rocking about x-axis
58.0 1.05e9 1.74 0.14 1 --- 0.010
s
k , c yy
k yy yy K yy , sur yy cyy 2k yy emb
Rocking about y-axis
73.9 5.55e9 1.58 0.12 0.99 --- 0.003
yy
(3) Calculated at a frequency corresponding to the first-mode period of the flexible-base structure.
hysteretic damping ratio, s. For the Sherman Oaks site, the soil hysteretic damping
radiation damping ratio (determined using equations in Chapter 2) and the soil
ratio was evaluated using Darendeli (2001), and determined to be 0.088 for the
Northridge earthquake, and 0.009 for the other events.
Vertical Stiffness and Damping. Calculations for overall vertical and rotational
stiffness and damping ratios are shown in Table 7-4, using values of average
effective profile velocity from Table 7-3. Vertical springs were distributed over the
footprint of the foundation, as shown in Figure 7-9, using a vertical stiffness intensity
that is normalized by area (Equation 2-20). Stiffness intensities were adjusted near
the edges of the foundation to match the overall rocking stiffness values given in
Table 7-4 (using Equation 2-21 and an end length ratio of Re = 0.4). The stiffness of
any individual spring in the model was then computed as the product of stiffness
intensity and the tributary area for the spring.
Vertical stiffness intensities in the central zone of the foundation were computed
using Equation 2-20a. Edge intensities were increased by factors of Rk = 5.33 for the
xx-direction and Rk = 3.58 for the yy-direction to correct for underestimation of
rotational stiffness (Equations 2-21a, 2-21b). Corner intensities were evaluated as the
average of the intensities in the xx- and yy-directions. Dashpot intensities were
calculated based on stiffness intensities using Equation 2-20b, and reduced by Rc
(Equations 2-21c, 2-21d), to correct for overestimation of rotational damping.
Results are shown in Table 7-5. Results are provided for 3x3 and 4x4 pile group
configurations and applied to individual pile groups shown in Figure 7-8. Results for
3x3 pile groups were applied to groups with 12 piles, and results for 4x4 pile groups
were applied to groups with more than 13 piles. The modeled stiffness for each pile
group is the product of the individual pile stiffness, k jp , the pile group efficiency
factor, k Gj , and the number of piles, Npile, in each group.
Table 7-5 Calculation of Pile Stiffness and Damping Parameters for the Sherman Oaks
Building
Static Pile Pile Group
Active Pile Stiffness, Efficiency
Pile Group
Length, La Vs,avg G(1) K jp (2) Factor,
j (3)
Spring/ Damping Ratio,
Dashpot (m) (m/s) (MPa) (MN/m) k Gj (3)
Horizontal, x, y
4.1 180.1 63.4 237.3 0.47 0.2
(3x3 group)
Vertical, z
9.9 183.7 67.6 859.3 0.43 0.25
(3x3 group)
Horizontal, x, y
4.1 180.1 63.4 237.3 0.36 0.05
(4x4 group)
Vertical, z
9.9 183.7 67.6 859.3 0.30 0.29
(4x4 group)
Notes: (1) Calculatedusing a shear modulus reduction factor, G/G0=0.5 for the Northridge earthquake.
(2) Calculatedfor Ep=474.8 MPa, and d=0.6 m.
(3) Approximated using Figure 2-11.
k y ,total k y ,base
k1 k2 k1
Translation in the k y ,base k y , piles
y-direction 48 16 16
k x ,total k x ,base
k4 k5 k 4
Translation in the k x ,base k x , piles
x-direction 27 6 6
k3 k zG k z
Translation in the
z-direction Note: k3 is not a constant value; is a function of Npiles, and kz is a function
of spring intensities shown in Figure 7-9.
Figure 7-10 Plan view, transverse section, and longitudinal section of foundation nodes
for the Sherman Oaks building, and calculation of associated spring
stiffnesses.
The overall horizontal stiffness, including contributions from base shear, passive
pressure resistance against the basement walls, and pile groups, was distributed
around the foundation perimeter (i.e., interior nodes have no horizontal springs).
Fundamental periods of vibration for each of the model variants are shown in Table
7-6. The resulting modeled periods were only modestly affected by different
idealizations of the soil-foundation interface.
Table 7-6 Comparison of Fundamental Periods for Alternative Foundation
Modeling Configurations for the Sherman Oaks Building
Fundamental Period
(sec)
Model Longitudinal Transverse
MB (Baseline Model) 2.67 2.72
MB.1 (rigid subterranean structure) 2.35 2.68
MB.2 (no tension in foundation springs) 2.65 2.73
Model 4 (bathtub) 2.67 2.72
Model 3 (fixed horizontal springs) 2.34 2.65
Model 2 (fixed at foundation) 2.67 2.71
Model 1 (fixed at grade) 2.34 2.67
Model MB.1 and Model MB.2. Results for Model MB.1 in Figure 7-11 showed that
rigid subterranean structural elements caused an increase structural response,
particularly in the NS direction, likely due to the change in period. In the case of
Model MB.2, results in Figure 7-11 showed that allowing geometric nonlinearities
(i.e., gap formation) had no discernible impact on response. This suggests that gaps
would not be expected to form in the foundation springs between the soil and the
basement walls, and supports the equivalent-linear soil-foundation modeling
assumptions suggested in Section 7.1.3.
Model 4. Results for Model 4 in Figure 7-12 showed that the bathtub model
introduces negligible changes in displacement response over the height of the
structure (i.e., less than 3.6% different at the roof level). Changes in story drift and
story shear force profiles were also negligible. Peak floor accelerations were most
sensitive to the change in modeling configuration. Values above the ground level
were relatively unaffected, but values in the subterranean levels were sensitive to the
Model 3. Results for Model 3 in Figure 7-13 showed the least agreement with Model
MB. Fixing the upper-level foundation springs against displacement, and applying
input motions to the base slab, caused large differences in all response quantities
including building vibration periods, displacement histories, drift ratios, and story
shears. Given the significant discrepancies observed, use of this modeling
configuration is not recommended.
Model 2. Results for Model 2 in Figure 7-14 showed that modeling the subterranean
levels, even while ignoring the effects of the surrounding soil, can provide good
results for some response quantities. Model 2 exhibited good agreement for building
vibration periods and displacement histories. Reasonable agreement was observed
for maximum displacement and drift ratios, but story shears and peak floor
accelerations differed more significantly, particularly in the subterranean levels.
Model 1. Results for Model 1 indicated that ignoring the subterranean levels
significantly alters the period of vibration. As a result, displacement histories were
more out-of-phase than most other modeling configurations, as shown in Figure 7-15.
Differences in story drifts, story shears, and peak floor accelerations were relatively
large (up to 50% different) in some cases.
In Figure 7-16, peak displacement, drift, and story shear response quantities from all
modeling configurations were synthesized and plotted in a single figure. Results for
Model 3 are clear outliers for each of the parameters considered. Results for Model 4
are closest to Model MB, followed by Model 2, and then Model 1. Differences in
response quantities, when they occurred, were generally greater in the subterranean
levels than in the levels above grade.
Analysis of the Sherman Oaks building showed that key building response quantities
can be affected by the foundation modeling assumptions, even when the structure-to-
soil stiffness ratio, h/(VsT), does not indicate a strong dependence on SSI effects.
Additional parametric studies were performed to investigate how different structural
system characteristics might influence the magnitude of SSI effects on modeled
response. A series of simplified stick models of the Sherman Oaks building were
developed, and the resulting response quantities for each model were compared.
An idealized stick model for the Sherman Oaks building is shown in Figure 7-17.
Although the superstructure has been simplified into a single equivalent stick, each
model includes detailed modeling at the foundation level, essentially maintaining
foundation geometry, and spring and dashpot configurations. Details on stick model
development and calibration with the full-building model are provided in
Appendix B.
Figure 7-17 Elevation an idealized stick model of the Sherman Oaks building.
Detailed analytical results for all parametric stick models are provided in
Appendix B. Results are compared between the moment frame stick model and the
Table 7-8 summarizes relative comparisons between the results for each response
quantity for the moment frame (MF), core shear wall (SW), and perimeter shear wall
(SW3) models. Information is based on plots provided in Appendix B. An
explanation for the nomenclature used in the table is provided in the footnotes.
Table 7-8 Comparison of Results for Moment Frame, Core Shear Wall, and Perimeter Shear Wall
Stick Models and Alternative Foundation Modeling Configurations
Moment Core Shear Perimeter Shear Wall
Frame (MF) Wall (SW) (SW3)
Response X-Direct. Y-Direct. X-Direct. Y-Direct. X-Direct. Y-Direct.
Quantity Location (E-W) (N-S) (E-W) (N-S) (E-W) (N-S)
Period --- 1 < 2, 4a, MB Similar 1 < 2, 4a, MB 1 < 2, 4a, MB 1 < 2 < 4a, MB 1 < 2, 4a, MB
Displacement Superstructure 4a, MB <1, 2 4a, MB <2<1 1<< 4a, MB <2 1<<2<4a, MB 1< 4a, MB <2 1<<2, 4a < MB
Basement 4a, MB <1, 2 Similar 4a,MB < 2 4a,MB < 2 4a,MB < 2 4a,MB < 2
Story Drift Superstructure 2, 4a, MB <1 4a<MB <2<1 Varies 1 < 2, 4a, MB Similar 1, 2< 4a < MB
Story Shear Superstructure 2, 4a, MB <1 2, 4a, MB <1 4a,MB <2<1 1< 4a, MB <2 4a< MB< 2<<1 2< 4a< MB<<1
Peak Floor
Superstructure Similar 4a< 1< MB< 2 4a < MB <1, 2 4a < MB <1, 2 4a < MB <1< 2 2, 4a, MB <<1
Acceleration
Basement 4a, MB<<2 4a, MB<<2 4a, MB < 2 4a, MB < 2 4a, MB<<2 4a, MB<<2
Notes: (1) 1, 2, 4a, and MB refer to Model 1, Model 2, Model 4a, and Model MB foundation modeling configurations.
(2) “1<<2, 4a < MB” indicates, for example, that the response of Model 1 is less than the response of Models 2 and 4a; the
response of Models 2 and 4a are similar; and the response of all models is less than the response of Model MB.
(3) “<<” indicates a comparatively larger difference in response than “<”.
(4) “Similar” indicates that the results do not vary significantly across the model types.
(5) “Varies” indicates that the trend in the results is not clear; some models have larger values in some cases and smaller
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
2 1-MF-Stick 2 1-MF-Stick
2-MF-Stick 2-MF-Stick
0 4a-MF-Stick 0 4a-MF-Stick
MB-MF-Stick MB-MF-Stick
-2 -2
-40 -30 -20 -10 0 10 20 30 40 -40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm] Max. Displacement [cm]
Northridge, X-Direction Northridge, Y-Direction
14 14
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
Figure 7-18 Comparison of maximum displacement envelopes in each direction for foundation Models
1, 2, 4a, and MB, for each structural system variant (MF, SW, SW1, SW2, and SW3).
Story Number
8 8
6 6
4 4
2 2 1-MF-Stick
2-MF-Stick
0 0 4a-MF-Stick
MB-MF-Stick
-2 -2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015 -0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm] Max. Interstory Drifts [cm/cm]
Northridge, Interstory Drifts, X-Direction Northridge, Interstory Drifts, Y-Direction
14 14
12 12
10 10
Story Number
Story Number
8 8
6 6
4 4
12 12
10 10
Story Number
Story Number
8 8
6 6
4 4
12 12
10 10
Story Number
Story Number
8 8
6 6
4 4
12 12
10 10
Story Number
Story Number
8 8
6 6
4 4
Figure 7-19 Comparison of maximum story drift envelopes in each direction for foundation Models 1, 2,
4a, and MB, for each structural system variant (MF, SW, SW1, SW2, and SW3).
Floor Number
MB-MF-Stick
Floor Number
8 8
6
6
4
4
1-MF-Stick 2
2 2-MF-Stick
0
0 4a-MF-Stick
MB-MF-Stick -2
-2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Normalized Story Shear
Normalized Story Shear
Northridge, X-Direction Northridge, Y-Direction
14 14
1-SW -Stick
12 12
2-SW -Stick
10 10 4a-SW -Stick
Floor Number
Floor Number
MB-SW -Stick
8 8
6 6
4 4
1-SW -Stick
2 2-SW -Stick 2
4a-SW -Stick
0 0
MB-SW -Stick
-2 -2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear Normalized Story Shear
Northridge, X-Direction Northridge, Y-Direction
14 14
12 12 1-SW 1-Stick
2-SW 1-Stick
10 10
4a-SW 1-Stick
Floor Number
Floor Number
8 8 MB-SW 1-Stick
6 6
4 4
1-SW 1-Stick
2 2-SW 1-Stick 2
4a-SW 1-Stick
0 0
MB-SW 1-Stick
-2 -2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear Normalized Story Shear
Northridge, X-Direction Northridge, Y-Direction
14 14
1-SW 2-Stick
12 12 2-SW 2-Stick
10 10 4a-SW 2-Stick
Floor Number
Floor Number
MB-SW 2-Stick
8 8
6 6
4 4
1-SW 2-Stick
2 2-SW 2-Stick 2
4a-SW 2-Stick
0 0
MB-SW 2-Stick
-2 -2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear Normalized Story Shear
Northridge, X-Direction Northridge, Y-Direction
14 14
1-SW 3-Stick
12 12 2-SW 3-Stick
10 10 4a-SW 3-Stick
Floor Number
Floor Number
MB-SW 3-Stick
8 8
6 6
4 1-SW 3-Stick 4
2 2-SW 3-Stick 2
4a-SW 3-Stick
0 0
MB-SW 3-Stick
-2 -2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear Normalized Story Shear
Figure 7-20 Comparison of maximum story shear envelopes in each direction for foundation Models
1, 2, 4a, and MB, for each structural system variant (MF, SW, SW1, SW2, and SW3).
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
2 1-MF-Stick 2 1-MF-Stick
2-MF-Stick 2-MF-Stick
0 4a-MF-Stick 0 4a-MF-Stick
MB-MF-Stick MB-MF-Stick
-2 -2
-1.5 -1 -0.5 0 0.5 1 1.5 -1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g] Peak Floor Accel. [g]
Northridge, X-Direction Northridge, Y-Direction
14 14
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
12 12
10 10
Floor Number
Floor Number
8 8
6 6
4 4
Figure 7-21 Comparison of peak floor acceleration envelopes in each direction for foundation Models
1, 2, 4a, and MB, for each structural system variant (MF, SW, SW1, SW2, and SW3).
The following general observations from parametric studies of the Sherman Oaks
building were made:
Observed SSI effects became larger as the building stiffness increased. They
were relatively minor for the moment frame (MF) building, more discernible for
the core shear wall (SW) building, and readily apparent for the stiffest perimeter
shear wall building (SW3).
Observed SSI effects correlated well with the structure-to-soil stiffness ratio,
h/(VsT); the larger the ratio, the greater the observed SSI effect.
The following additional observations were made, relative to the specific case of
these building variants and these input motions:
For all building variants in this study, fundamental periods are past the peak
spectral response and located on the descending (i.e., velocity-controlled) branch
of the response spectrum for the Northridge earthquake. As such, any stiffening
in the structure could be assumed to lead to an increase in response.
Observed differences in modeled periods between the different foundation
models were relatively small for a given structure. A possible exception
occurred in the case of shear wall buildings, in which Model 1 periods were
noticeably less than the other models.
Large differences were observed between model results in the basement levels.
Different trends were observed in comparisons between superstructure results
and comparisons between basement results.
As the superstructure stiffened, model-to-model variations in superstructure
response increased such that they approached the magnitude of observed
variations in basement level response. This was particularly noticeable in the
case of peak floor acceleration.
In the superstructure, Model 4a overall results were the most similar to Model
MB results for all building types. Model 1 and Model 2 were usually
conservative for story shear and peak floor acceleration for all building types.
Model 1 was observed to result in conservative displacements and drifts in the
moment frame building, but unconservative displacements and drifts in the shear
wall buildings.
In the basement levels, Model 2 results were higher than Models 4a results across
all response quantities. In the case of displacements and peak floor accelerations,
Model 2 results exceeded Model MB results, but in the case of story drift and
story shear, Model MB often had larger values than Model 2.
The Walnut Creek building is a10-story structure with no subterranean levels. The
seismic force-resisting system consists of a reinforced concrete shear wall core with
perimeter precast and cast-in-place concrete frames. Construction drawings of the
Walnut Creek building were made available for inspection through the auspices of
the California Strong Motion Instrumentation Program. The building was designed
in 1970. The plan dimensions of the building are 31.8 m (104 ft) wide by 45.2 m
(148 ft) long. A typical floor plan is shown in Figure 7-22.
Figure 7-22 Walnut Creek building typical floor plan, based on construction
drawings provided by CSMIP.
A plan view of the Walnut Creek site is shown in Figure 7-24. The ground surface
elevation is approximately 51.8 m (170 ft) above mean sea level. Soil conditions for
the site were obtained from portions of a geotechnical report, prepared by Harding,
Miller, Lawson & Associates (1970), and from a geotechnical report with seismic
refraction investigation of the subsurface soils, prepared by Raney Geotechnical
(1983), in the vicinity of the site. Information from a total of 19 borings ranging in
depth from 1.5 m (5 ft) to 17.4 m (57 ft), and four seismic refraction surveys, were
obtained from within and around the footprint of the Walnut Creek building.
Borings indicated predominantly west-dipping contacts of sandy clays and silts with
variable thicknesses of 0.6 m (2 ft) to 5.5 m (18 ft) overlaying siltstone and sandstone
of the Orinda Formation. The depth of the water table was measured between 6.1 m
(20 ft) and 12.2 m (40 ft) below the ground surface. The average moist unit weight
Figure 7-24 Plan view of the Walnut Creek site showing borehole and refraction
survey locations used for site characterization (adapted from
Harding, Miller, Lawson & Associates, 1970; Raney Geotechnical,
1983).
The Walnut Creek building was instrumented in 1979 by the California Strong
Motion Instrumentation Program and designated CSMIP Station No. 58364.
Instrumentation includes 16 accelerometers at the locations including the ground
floor, third floor, eighth floor, and roof levels. There are two vertical sensors at the
ground level, allowing the base rocking effects to be measured. There are no free-
field instruments in the vicinity of the site.
Since 1979, five earthquake events have been recorded and processed by CSMIP at
this station. The events, along with peak accelerations for the ground floor and roof
are presented in Table 7-9.
Table 7-9 Summary of Earthquake Events Recorded at the Walnut Creek Building
Ground Floor Roof
Earthquake (g) (g)
N-S E-W V1 V2 N-S E-W V
Livermore 80A 0.030 0.033 0.022 0.023 0.116 0.108 ----
Livermore 80B 0.061 0.048 0.024 0.018 0.137 0.192 ----
1989 Loma Prieta 0.102 0.046 0.053 0.047 0.208 0.164 ----
2007 AlumRock 0.005 0.003 NR NR 0.018 0.015 ----
2008 Alamo 0.029 0.018 NR NR 0.057 0.034 ----
Symbols: N-S=north-south; E-W=east-west; V=vertical; V1=Chan.12 (west); V2=Chan.13 (East); and
NR=no record available.
Figure 7-26 CSMIP Station No. 58364: Walnut Creek 10-Story commercial building,
sensor location sketch (CESMD, 2011).
Since there were no free-field instruments in the vicinity of the site, free-field ground
motions, ug, were obtained from the recorded motions using transfer functions to
remove the base-slab averaging effects. The uFIM and ug response spectra for the
1989 Loma Prieta earthquake are shown in Figure 7-27.
Figure 7-27 Response spectra for foundation input motion, uFIM, and free-field motion, ug, at the
Sherman Oaks building, 1989 Loma Prieta earthquake.
The foundation is rectangular in plan, measuring 31.8 m (104 ft) wide by 45.2 m
(148 ft) long. The foundation consists of concrete spread footings, drilled shafts, and
mat elements. A foundation plan developed based on information contained in
construction drawings provided by CSMIP is shown in Figure 7-28.
Values Vs, avg in Table 7-10 have not been corrected for overburden pressures due to
the added weight of the structure. Before averaging, measured free-field velocities,
Vs, should be increased to account for the presence of structural overburden.
Horizontal translation (x and y), base spring BL 8.6 4.3 to 12.9 434.3
Notes: (1) Calculated using mat foundation half-width, B=5.8 m; half-length, L=12.8 m; depth D=4.3 m.
structure, v, was evaluated at two-foot intervals below the foundation assuming a
For the Walnut Creek building, the increase in vertical stress due to the weight of the
Mat Foundation Flexibility. The Walnut Creek building has a central core
consisting of shear walls supported by an underlying reinforced concrete mat
foundation that is larger in plan and extends beyond the dimensions of the core.
To evaluate the effects of this flexibility on the resulting foundation impedance, plots
similar to Figure 2-7 (from Iguchi and Luco, 1982) were utilized with a ratio of core
to foundation radius equal to 0.75. Results are shown in Table 7-12, which indicate
that the foundation stiffness ratios are close to one, meaning there is little deviation
from the rigid foundation case. Therefore, the full mat dimensions of 11.6 m (38 ft)
by 25.6 m (84 ft) were used to calculate foundation impedance below the core walls
in the Walnut Creek building.
Table 7-12 Evaluation of Soil-to-Foundation Stiffness Ratios for Flexible Mat Foundation
below Shear Wall Core in the Walnut Creek Building
(1)
Ratio
Basis for Core Mat rc k (2)
Equivalent Radius Radius
Vibration Mode Radius, rf rc (m) rf (m) rf (Eq. 2-19) k 0
Assessed based on Iguchi and Luco (1982), using rc/rf =0.75 and =100.
Notes: (1) Calculated using Ef =22,894.7 MPa; f =0.2, tf =1.22 m, a0 =0.13, and G based on Vs,F,avg from Table 7-11.
(2)
density, s, of 18.1 kN/m3/g, and values of average effective profile velocity
In Table 7-13, shear modulus, G, was evaluated from Equation 2-9, using a soil mass
s
Horizontal, base spring
k x ,base x K x , sur cx,base 2k x ,base sur
k , cx ,base
x-direction 469.4 2.24e7 --- 0.11 1 0.057 ---
x ,base
s
Horizontal, base spring
k y ,base y K y , sur cy ,base 2k y ,base sur
k , c y ,base
y-direction 469.4 2.40e7 --- 0.11 1 0.053 ---
y ,base
s
k z z K z , sur z cz 2kz emb
kz , cz
Vertical, z
376.4 2.35e7 1.22 0.12 1 --- 0.111
s
k xx xx K xx , sur xx cxx 2kxx emb
kxx , cxx
Rocking about x-axis
371.6 8.49e8 1.70 0.12 1 --- 0.005
s
k , c yy
k yy yy K yy , sur yy cyy 2k yy emb
Rocking about y-axis
379.3 2.78e9 1.54 0.12 0.99 --- 0.002
yy
(3) Calculated at a frequency corresponding to the first-mode period of the flexible-base structure (from system identification, Stewart and Stewart, 1997).
Dashpot coefficients were determined using Equation 2-13a and the total foundation
damping ratio. The total foundation damping ratio is equal to the sum of the
hysteretic damping ratio, s. For the Walnut Creek site, the soil hysteretic damping
radiation damping ratio (determined using equations in Chapter 2) and the soil
Vertical stiffness intensities in the central zone of the mat foundation were computed
using Equation 2-20a. Edge intensities were increased by factors of Rk = 5.48 for the
xx-direction and Rk = 3.21 for the yy-direction to correct for underestimation of
rotational stiffness (Equations 2-21a, 2-21b). Corner intensities were evaluated as the
Footing and Drilled Shaft Stiffness and Damping. Individual springs and dashpots
were calculated for horizontal and vertical translation modes for the spread footings
and drilled shafts supporting the columns of the perimeter frames. For simplicity,
spread footings and drilled shafts were treated as individual shallow square
foundations with dimensions of 3.05 m (10 ft) by 3.05 m (10 ft). Calculations were
similar to those presented for the mat foundation, using the same soil parameters.
For reference, calculated values were 3.91 x 106 kN/m for horizontal springs (kx and
ky); 4.97 x 106 kN/m for vertical springs (kz); 1.68 x 103 kN-s/m for horizontal
dashpots (cx and cy); and 2.75 x 104 kN-s/m for vertical dashpots (cz).
Because the Walnut Creek core wall foundation has limited embedment, and there
are no basement levels, multi-support excitation along embedded portion of the
building was not considered. As a result, Model MB converges to Model 4 (bathtub
model). Using this model as a basis, alternative modeling configurations were
developed and studied, and the resulting response quantities were compared for the
following variants:
Model 2, which included explicit modeling of the subterranean foundation
elements, assuming a fixed base at the foundation level and omitting the
surrounding soil; and
Model 1, which ignored the response of the subterranean foundation elements by
assuming a rigid base at the ground level.
Implementation of Model 4 (bathtub model), Model 2, and Model 1 for the Walnut
Creek building is illustrated in Figure 7-30. Because of the minimal embedment of
the Walnut Creek building structure, Model 1 and Model 2 are nearly identical.
Fundamental periods of vibration for each of the model variants are shown in Table
7-14. The resulting modeled periods were only modestly affected by different
idealizations of the soil-foundation interface (i.e., Model 1 and Model 2 have the
same fundamental period; Model 4 is slightly more flexible with a longer
fundamental period).
Table 7-14 Comparison of Fundamental Periods for Alternative Foundation
Modeling Configurations for the Walnut Creek Building
Fundamental Period
(sec)
Model E-W Direction
Model 4 (bathtub) 0.83
Model 2 (fixed at foundation) 0.78
Model 1 (fixed at grade) 0.78
In Figure 7-31, the fixed-base models (Model 1 and Model 2) capture roof
acceleration and displacement histories reasonably well, but the phasing was quite
different. Observed differences in phasing are likely due to differences in higher
mode response.
The presence of two vertical instruments at the ground level allowed for the
evaluation of base rotation, and computation of lateral roof displacement due to base
rocking. High-pass filtering of the base rocking data, which was needed to produce a
physically meaningful result, is described in Appendix C. Comparison of modeled
versus recorded roof displacement histories due to base rotation are shown in Figure
7-33.
Based upon previous findings from Naeim et al. (2008), Tileylioglu et al. (2010), and
observations from example applications and parametric studies herein, the following
conclusions can be made:
MB-type models have an encouraging ability to match observed building
responses from recorded motions with modest tuning of structural parameters
(e.g., damping ratios, building masses, or element stiffnesses). These results
suggest that the relatively simple equivalent-linear spring and dashpot approach
for SSI modeling presented in Chapter 2 can provide a satisfactory representation
of foundation impedance.
Of the model types studied, Model 4, including foundation springs and dashpots
within a rigid bathtub element, provided the best match to MB-type models, and
is a reasonable and practical simplification to variable support excitation.
Model 2, which omits foundation springs, but explicitly includes modeling of the
subterranean structure, was the next best alternative. Model 1 results were
observed to vary significantly from MB-type model results, but generally
provided a conservative estimate of force-based demands for design use.
Model 3 results were highly variable, and use of Model 3 approaches to
foundation modeling is not recommended.
Above-ground building responses, as measured by envelopes of peak response
parameters such as displacement, drift, story shear, and floor acceleration, were
only modestly affected by soil-structure interaction effects in the buildings
studied. The effects became more significant as the stiffness of the
superstructure increased. In contrast, below-ground responses were very
sensitive to soil-structure interaction effects and foundation modeling
assumptions (i.e., kinematic ground motion descriptions and spring/dashpot
configurations). Sensitivity of below-ground response was observed for all
structures, across the full parametric range of stiffnesses studied.
In this chapter, conclusions and recommendations from prior chapters are distilled
into specific, step-by-step procedures to guide the implementation of soil-structure
interaction modeling in a design setting. Section 8.1 answers the three questions
posed in Chapter 1, which are intended to help guide decisions regarding when SSI is
likely to be important and what geotechnical information is needed to model SSI
effects. Section 8.2 provides guidance on implementation of inertial and kinematic
interaction effects. Section 8.3 describes further study and future research needed to
improve the state of knowledge for soil-structure interaction.
Question 1: When is the use of foundation springs and dashpots important, and
which structural response parameters are affected?
The structure-to-soil stiffness ratio, h/(VsT), can be used as a relative measure for
determining when SSI effects will become significant. In this expression, h is the
structure height, Vs is the soil shear wave velocity; and T is the fixed-base building
period. In applying the structure-to-soil stiffness ratio, values of h, Vs, and T for a
given soil-foundation-structure system are not unique. The following values should
be used:
Height. Height, h, is the effective height to the center of mass for the first mode
shape, taken as approximately two-thirds of the modeled building height.
Shear wave velocity. Shear wave velocity should be taken as the average
effective profile velocity, Vs, avg, calculated based on overburden-corrected shear
wave velocities below the foundation, Vs, F (z). Guidance on calculating average
effective profile velocities is provided in Section 2.2.2 and illustrated in
Chapter 7.
Period. Period should be taken as the best estimate value of the fixed-based
building period in the direction under consideration. The structure-to-soil
stiffness ratio should be evaluated separately in each direction.
When using the structure-to-soil stiffness ratio, it is important to recognize that the
ratio is an approximate relative measure, and not an absolute criterion. Even when
h/(VsT) < 0.1, relative distributions of moments and shear forces in a building can be
modified relative to the fixed-base condition, particularly in dual systems, structures
with significant higher-mode responses, and subterranean levels of structures.
Additional information is provided in Chapter 2, and SSI effects on other building
response quantities are investigated in Chapter 7.
Base-slab averaging effects become important within the period range of engineering
interest for foundation sizes (measured as an equivalent foundation half width) of
about 20 m (66 ft) or larger for typical soil or weathered rock sites in California.
Embedment effects are sensitive to the depth of embedment, and typically become
important when a structure has two or more subterranean levels. Both base-slab
averaging and embedment effects principally impact short period spectral ordinates
(at periods less than approximately 1.0 sec). Base-slab averaging and embedment
effects are introduced in Chapter 3. Implementation in standards and guidelines is
described in Chapter 4.
In this section, procedures for developing springs and dashpots, modifying ground
motions due kinematic interaction, and including soil-structure interaction in
response history analyses are summarized, along with references to more detailed
information.
The steps for developing springs and dashpots in SSI analyses are summarized as
follows:
1. Develop the required input parameters for analysis:
a. Geotechnical and shear wave velocity profiles. Examples of the type of
information that is required are given in Figure 7-4 and Figure 7-25.
Uncertainty in the shear wave velocity profile should be considered on the
basis of scatter in data from measurements at the site.
b. Shear strength parameters and their variation with depth. Below the ground
water table, undrained strength parameters are required. Drained strength
parameters are generally acceptable above the ground water table.
Variability in shear strengths should be considered based on the range
observed in material-specific testing.
c. Poisson’s ratio, (. Can generally be taken as 0.3 for sands and 0.45 for
clays.
d. Soil hysteretic damping ratio, (s). Strain-dependent soil damping can be
measured using site-specific dynamic material testing, but for most projects
can be taken from existing empirical relationships. Several such
taken as the average of the first mode periods in each orthogonal direction.
b. Calculate the dimensionless frequency, a0, using Equation 2-15.
c. Calculate the dynamic stiffness modifier, j, using equations in Table 2-3a.
d. Calculate radiation damping ratios, j, using equations in Table 2-3a for
surface foundations, or Table 2-3b for embedded foundations.
5. Evaluate frequency-dependent spring and dashpot (impedance) coefficients:
modifier for embedment, j, and the dynamic stiffness modifier, j, using
a. Stiffness, kj, can be taken as the product of the static stiffness, Kj,sur, the
Equation 2-14a.
b. Dashpot coefficient, cj, can be evaluated from radiation damping ratios, j,
using Equation 2-13a.
c. Repeat as needed for the range of static stiffnesses that account for
uncertainty in the soil profile.
6. Evaluate limiting spring forces (i.e., capacities), which depend on the strength of
the foundation soils, as described in Section 2.2.4:
a. For vertical springs, the capacity is the unfactored foundation bearing
capacity.
b. For horizontal springs at the base-slab level, the capacity is derived from
shear-sliding resistance at the soil-foundation interface.
c. For basement walls, the capacity is derived using passive earth pressure
theory.
d. Repeat as needed to account for uncertainty in shear strength parameters.
7. Distribute vertical springs and dashpots around the foundation:
damping).
Procedures for deep (e.g., pile) foundations, described in Section 2.3, are similar to
those presented above for shallow foundations. Differences include the depth range
over which dynamic properties are evaluated; the equations used to calculate static
stiffness, dynamic modifiers, and damping ratios; and the need for group modifiers in
the case of pile groups. Computer programs (e.g., DYNA6, Western Engineering,
2011; SASSI, Lysmer et al., 1999) can perform Step 2 through Step 5 for finite
element analyses, utilizing procedures similar to those discussed in this report. The
DYNA6 program can also be used for analysis of pile impedance as well.
An important consideration when deep foundation elements (e.g., piles) are combined
with shallow foundation elements (e.g., spread footings or mats) is whether or not
resistance from both shallow and deep foundation elements can be combined. If the
soil is expected to settle away from the shallow foundation elements (e.g., the case of
consolidating soils and end-bearing piles), then lateral load resistance should be
derived on the basis of the piles, pile caps, and basement walls only, and the
resistance provided by shallow foundation elements should be ignored.
The steps for modifying ground motions due to kinematic interaction are summarized
as follows:
1. Collect the required input for the analysis:
a. Specification of seismic demand. Determine if seismic demands for design
are to be specified in the form of an acceleration response spectrum a set of
ground motion acceleration time histories.
b. Foundation Dimensions. Determine the area of foundation, as represented
by B eA , and the embedment depth, D.
c. Shear wave velocity profile. Determine shear wave velocities to depth, D.
2. Calculate the transfer function for base-slab averaging:
a. Transfer function, Hu(), is calculated using Equation 3-3 and Equation 3-4.
Input parameters include B eA and a.
b. Parameter, a, is evaluated using Equation 3-5. This base-slab averaging
model is calibrated for B eA in the range of 15 m to 40 m.
Hu(), for embedment is calculated using Equations 3-6a through 3-6d. Input
3. Calculate the transfer function for embedment (if D > 0). The transfer function,
parameters include depth, D, and the average effective profile velocity, Vs, avg,
over that depth.
4. For each frequency, , the combined transfer function ordinate for base-slab
averaging and embedment is taken as the product of the results from (2) and (3).
5. Modify free-field ground motions to foundation input motions:
a. Calculate response spectrum modifiers using Equation 3-7a or Equation
3-7b, which relate the ratio of foundation to free-field response spectral
ordinates to the transfer function at the corresponding frequency.
b. Modify acceleration time histories, as needed, using the procedures in
Section 3.4. Note the limitations described in Section 3.4 are for ground
motions dominated by long-period energy (e.g., soft soils and near-fault
effects).
Future research needs in soil-structure interaction have been organized into two
general thematic areas. The first theme involves relatively short-term
recommendations expanding on current studies to: (1) provide tangible insights into
the benefits of SSI analysis for owners and practicing engineers; and (2) further
explore the benefits and limitations of SSI response history analysis procedures,
possibly leading to improved procedures. The second theme involves relatively long-
term recommendations intended to address fundamental limitations in the state of SSI
knowledge, which limit the accuracy and reliability of SSI models available for use
in engineering practice.
Research needs for expanding the state of knowledge for soil-structure interaction are
as follows:
The foundation damping model of Veletsos et al. (various) produces different
results than similar models by others. A critical examination of the derivation of
that model is needed, followed by the development of equations for foundation
damping that properly consider hysteretic damping from soil response, radiation
damping from rotational and translational vibration modes, and the sensitivity of
radiation damping to different soil stiffness profiles. This problem is discussed
in Section 2.1.
The rotational stiffness of shallow foundation systems with non-rigid structural
foundation elements is poorly understood. In particular, the effects of coupled
versus uncoupled rotations at the base of lateral-load bearing elements on
radiation damping and overall system impedance need to be investigated. This
problem is discussed in Section 2.2.3.
The impedance of pile-supported foundations is poorly understood for realistic
pile and soil conditions, especially for pile groups. Elasto-dynamic solutions for
piles in idealized soil profiles exist in the literature, but are not used in practice,
in part because they only apply at low displacement levels. The discrete element
models that are used (e.g., LPILE, APILE) are poorly constrained for stiffness,
and are intended for non-seismic problems. Next-generation element models for
dynamic loading of piles are needed that accurately capture the stiffness from
Renderings of the Sherman Oaks building configuration are provided in Figure A-1.
The renderings include cut-away views showing details in the first story, basement,
and foundation regions. The structural configuration in the first story and basement
levels were modeled in a simplified manner to keep the structural modeling aspects
tractable, and to focus on comparisons between variations in the idealization of the
soil-foundation interface.
In the first story of the building, the one-story extension of the structure along the
south longitudinal face of the building (shown in Figure A-1) was not included in the
structural model (i.e., both the mass and the stiffness of the one-story extension were
excluded). In the basement region, the geometry was simplified such that the shear
walls were placed in-line with the perimeter frames of the superstructure. This kept
the footprint of the model at 21.9 m wide by 57.6 m long (72 ft by 189 ft) from
foundation to roof. Although this geometric simplification was made, the flexibility
inherent in the real configuration was captured through elastic connecting springs
modeled with the stiffness properties of the horizontal slab components that
interconnect the frame lines of the superstructure with the basement wall lines. This
is shown schematically in Figure A-2.
Figure A-2 Schematic illustration of elastic springs connecting the framing lines of the
superstructure with the wall lines in the basement levels.
Fiber elements are composed of both core and cover concrete materials (Concrete04
in OpenSees) and reinforcing steel materials (Steel02 in OpenSees). Fiber elements
model flexural behavior using nonlinear concrete and steel material models. Shear
and torsional flexibilities are modeled to be linear-elastic, and are combined together
into a single component model using the SectionAggregator approach in OpenSees.
An expected yield strength of 462 MPa (67 ksi) was used for the grade 60 steel
(Melchers, 1999). This value is slightly lower than, but still comparable to, an
expected yield strength of 517 MPa (75 ksi), which is suggested in ASCE/SEI 41-06
(ASCE, 2007). An expected initial stiffness of 20,000 GPa (29,000 ksi) was used,
along with a post-yield hardening stiffness of 2% of the initial stiffness. A nominal
concrete strength of 35 MPa (5.0 ksi) was included in the model, without provision
for expected strength. Use of nominal concrete strength in lieu of expected strength
was judged to have minimal impact on the structural response predictions because of
the mild nonlinearity experienced in the Northridge earthquake.
The calculated building mass included the mass of all structural elements (beams,
columns, and slabs); 0.5 kPa (10 psf) for partitions; 0.6 kPa (12 psf) for mechanical,
electrical, and plumbing components; and 25% of the design live load taken as 0.6
kPa (12 psf) for a design live load of 2.5 kPa (50 psf). These masses are accounted
for both above and below grade.
Damping was modeled as 4.5% Rayleigh damping, anchored to the first and second
mode periods of the building (2.9 sec and 1.0 sec, respectively). This level of
damping is a calibrated value, which is discussed below. In developing the Rayleigh
damping matrix, degrees-of-freedom associated with foundation springs were
excluded. This was necessary to avoid double counting of foundation damping
because the soil-foundation model included dashpots at these degrees-of-freedom.
Basement shear walls are 30 cm (1 ft) thick, modeled in the simplified manner
described above (and shown in Figure A-2). Figure A-2 shows how the moment
resisting frame, simplified shear wall, and slab elements connect to the framing nodes
in the two basement levels of the building. At each connection location, the nodes
are placed at the same coordinates in the model (but are shown offset in Figure A-2
to illustrate the connectivity). Vertical and rotational degrees of freedom are
constrained together for all nodes at each connection location. Figure A-2 also
depicts the rigid beams needed to support the simplified shear wall models and the
soil springs at the base of the building (soil dashpots are not shown). Stiffness in the
The stiffness and damping properties of the soil are modeled using vertical and
horizontal springs and dashpots, as described in Chapter 7. Soil springs are linear-
elastic (with model variant MB.2 considering the no-tension gap springs). Dashpots
are linear in all model variations.
The depth-variable nature of input ground motions over the height of the basement
walls were specified using the MultipleSupportExcitation approach in OpenSees. In
this approach, the acceleration, velocity, and displacement acceleration-histories are
all specified at each subterranean level. Specifying each ground motion history
removes the need to integrate motions within OpenSees.
The Sherman Oaks Baseline Model (MB) was calibrated against recordings from the
1994 Northridge earthquake. The calibrated model was then used to predict the
response of the building for the 1992 Landers and 1987 Whittier earthquakes to
assess the stability of the calibration.
The first stage of calibration targeted the shear stiffness of the reinforced concrete
walls and slabs in the subterranean levels. A multiplier on the theoretical shear
modulus of the uncracked concrete, Gc, was used as the calibration parameter to
match near-ground response. A definitive ratio of cracked to uncracked stiffness in
shear walls is not available in the literature. For example, ASCE/SEI 41-06 suggests
that the full (unreduced) value of Gc be used for both cracked and uncracked walls.
PEER/ATC-72-1 states that the cracked shear stiffness should be “substantially
lower” than the uncracked shear stiffness, but there are limited test data available for
use in quantifying this ratio. It was found that, over the range of 0.25 to 0.40
considered, near-ground response was not highly sensitive to this parameter, and
acceptable results were achieved with a stiffness multiplier of 0.25.
The second stage of calibration targeted the building period. Adjustments to the
structural mass and stiffness were used as the calibration parameters for the modeled
building period. The structure mass was modified through the application of a scale
factor. In fiber element models, stiffness is not an input parameter, but is computed
from component dimensions and material properties. Component stiffness was
adjusted through modification of the stiffness of reinforcing steel within the fiber
element model. This approach is not meant to suggest that the rebar stiffness is
The initial modeled building period was found to be too low. To increase the period,
final calibration factors included multipliers of 0.7 on rebar stiffness and 1.0 on mass.
Table A-1 compares the first mode periods of the building in each direction with the
periods predicted by the Baseline Model (MB). Measured building periods were
computed using non-parametric system identification through the calculation of
transfer functions between the base-level and roof acceleration histories. Although
not an exact match in all cases, Table A-1 shows that modeled periods agree
reasonably well with the measured periods in the three events studied.
Table A-1 Comparison of Measured and Modeled Periods for the Sherman Oaks Building
Identified Period Identified Period MB Period MB Period
Earthquake Mode 1 (sec) Mode 2 (sec) Mode 1 (sec) Mode 2 (sec)
Long Trans Long Trans Long Trans Long Trans
1994 Northridge 2.92 3.28 0.84 0.94 2.67 2.72 0.84 0.94
1992 Landers 2.56 2.72 0.73 0.80 2.14 2.68 0.74 0.80
1987 Whittier 2.33 2.29 0.82 0.83 2.33 2.29 0.82 0.83
The third stage of calibration targeted Rayleigh damping. Adjustments to the level of
damping were used to match the amplitude of displacements in the superstructure.
The OpenSees model used full Rayleigh damping with two matching frequencies,
taken as the first- and second-mode vibration frequencies. A target level of 4.5% of
critical damping was found to provide reasonable results.
Soil spring capacities were compared to demands computed for the 1994 Northridge
using the calibrated Model MB to investigate the potential for failure of the
foundation during earthquake shaking. Limiting spring capacities for the Sherman
Oaks building (presented in Chapter 7) exceeded spring force demands by a factor of
three or more. This margin of safety indicates that foundation failure is unlikely, and
limiting spring capacities were not needed in the present simulations.
The calibrated Model MB was then used to predict the response of the building for
the 1992 Landers and 1987 Whittier earthquakes to assess the stability of the
calibration. Results for displacement histories and peak displacements are shown in
Figure A-5 through Figure A-8. The quality of the match for Whittier is better than
for Northridge, likely due to a lack of nonlinearity in the response. In the case of
Landers, the predicted response is consistently weaker than the recorded response.
The Sherman Oaks building moment frame is a 13-story structure above grade, with
a total of seven bays in the east-west (E-W) direction (X-direction) and a total of two
bays in the north-south (N-S) direction (Y-direction). The overall configuration of
the full-building Baseline Model (MB) is shown in Figures B-1 through B-3. These
figures also show the location of the master nodes used in the stick model. The stick
model configuration is shown in Figure B-4.
In the stick model, the mass of each floor, the effective lateral, axial, and rotational
stiffness, and the effective yield strength, are all lumped at the master node for each
level. Geometric and material definitions from the full-building OpenSees model
code were extracted and the rotational and axial stiffness at each floor level were
calculated. A static pushover analysis of the full-building Model MB was used to
determine the lateral stiffness and yield strength at each level.
First, a series of static pushover analyses were performed involving the application of
a single point load at the master node of the floor under consideration, while
constraining all other master nodes of the structure. The pushover analysis was
displacement-controlled, and the recorded force-displacement data at the master node
Figure B-2 Longitudinal elevation (X-direction) of the full-building model of the Sherman
Oaks building.
under consideration was used to determine the effective lateral stiffness at each floor
level. Adjustments to the effective lateral stiffness were made as needed to account
for the stiffness contribution of adjacent floors. Once all contributions to stiffness
had been accounted for, force-displacement curves were created. These force-
displacement curves were then used to create idealized force-displacement curves per
ASCE/SEI 41-06 (ASCE, 2007). An example force-displacement curve is shown in
Figure B-5.
This methodology, however, resulted in overall stiffness and yield strength values for
the stick model that exceeded the full-building Model MB values. Another pushover
analyses was performed involving the application of a single point load at the master
node of the roof, pushed to a defined displacement of 20 cm (8 in). This
5000
Force (kips)
4000
3000
2000
1000
0
0 1 2 3 4 5 6 7
Displacement (in)
Figure B-5 Idealized force-displacement curve for the Sherman Oaks building.
Figure B-6 Shear wall and basement level node definitions, constraints, and
connectivity for the Sherman Oaks stick model.
The first stage of calibration targeted the hysteretic behavior observed at each floor
during the response histories. Hysteresis curves were used to compare the predefined
lateral stiffness of the stick model to the full-building Model MB. Softening of the
lateral stiffness of the stick model was needed in both lateral directions at all floor
levels. Plots of absolute displacement and peak acceleration response showed greater
convergence after this calibration.
The second stage of calibration targeted the mild bilinear behavior observed in story
levels above the first story. This stage consisted of defining the yield force at each
level where bilinear behavior occurred in an attempt to match maximum floor
displacements and peak floor accelerations. Each iteration caused divergence of
results with no clear trend toward better calibration, so further adjustment of this
parameter was abandoned.
Final Stick Model Properties. The final stick model was created using
twoNodeLink element objects. This allowed the stick elements spanning from floor
to floor to be modeled as axial, shear, and rotational springs. Final element
properties for the moment frame stick model are given in Tables B-1 through B-3.
The axial and rotational spring stiffness values in Table B-1 are based on length
(story height), area (total column area per story level), and elastic modulus
information taken directly from the full-building Model MB. Rotational stiffness
about global X- and Y-directions were determined by solving a system of linear
equations for a unit rotation. The lateral force-deformation properties in Table B-2
are based on lateral pushover analyses described above. Masses (weights) used in the
stick model are provided in Table B-3.
Table B-2 Lateral Force-Deformation Properties of the Moment Frame Stick Model
Story Yield Strength (kips) Initial Stiffness (kips/in) Strain-Hardening Ratio
E-W N-S E-W N-S E-W N-S
1 5500 6500 2566.8 2521.9 0.0371 0.0142
2 12000 13050 5342.1 4415.0 -0.0541 -0.0176
3 10171 11605 4893.8 3755.8 -0.0348 -0.0478
4 9600 10200 4724.2 3521.5 -0.0318 0.0000
5 8000 8214 4584.9 3429.5 -0.0643 -0.0234
6 6700 7300 4452.4 3389.5 -0.0460 -0.0255
7 5286 6000 4375.8 3351.1 -0.0345 -0.0312
8 5000 5733 4261.2 3293.0 -0.0210 -0.0370
9 4343 5086 4169.9 3242.0 0.0183 -0.0225
10 3771 4571 4141.1 3205.5 0.0040 -0.0090
11 3000 3743 4142.6 3215.3 0.0245 0.0339
12 3000 3600 4235.6 3346.7 0.0263 0.0115
13 5500 5000 4920.4 3864.7 0.0820 0.2642
Table B-4 compares the first-, second-, and third-mode periods of the moment frame
(MF) stick model and the full-building Model MB. Development and calibration of
the stick model was performed using a relatively early version of the Sherman Oaks
OpenSees Model MB, which changed over time as model components and ground
motion inputs were refined. For this reason, Model MB results presented in this
appendix are not an exact match to those presented in Section 7.2.
Table B-4 Comparison of Modeled Periods for the Moment Frame Stick
Model and the Full-Building Model MB
Identified Period Identified Period Identified Period
Mode 1 (sec) Mode 2 (sec) Mode 3 (sec)
Translational Torsional Translational
Model (Global Y-axis) (about Global Z-axis) (Global X-axis)
MB-MF-Full 2.45 2.31 2.21
MB-MF-Stick 2.87 N.A. 2.62
Figure B-7 compares the roof displacement response history for the moment frame
stick model to the center of the full-building Model MB in the X-direction and the
Y-direction. Responses in the Y-direction are more in phase than in the X-direction.
Displacement amplitudes match fairly well. In some cases, larger amplitudes are
observed in the stick model, but in other cases, larger amplitudes are observed in the
full-building model. Similar results were obtained for nodes located at the mid-
height of the building.
10
-10
-20
-30
-40
0 10 20 30 40 50 60
Time (sec)
10
-10
-20
-30
0 10 20 30 40 50 60
Time (sec)
Figure B-7 Comparison of roof displacement histories for the moment frame
stick model and the full-building Model MB.
Figures B-8 through B-11 compare response envelopes for maximum displacement,
story drift, story shear, and peak floor acceleration for each direction. Results for the
moment frame stick model are typically within about 10% of those for the full-
building Model MB. Maxima in the positive and negative directions are reported. In
some cases, results for the moment frame stick model are larger in one direction, but
smaller in the perpendicular direction. It was found that calibration attempts often
improved results in one direction at the expense of the other direction.
Based on the following comparisons, the calibration process for the moment frame
stick model was judged to be sufficient for the purposes of the parametric study.
12
10
Floor Number
8
0 MB-MF-Stick
MB-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
0 MB-MF-Stick
MB-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
12
10
Story Number
0 MB-MF-Stick
MB-MF-Full
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-9 Comparison of maximum story drift envelopes for the moment frame
stick Model MB and the full-building Model MB.
12
10
Floor Number
8
0 MB-MF-Stick
MB-MF-Full
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Northridge, Y-Direction
14
12 MB-MF-Stick
MB-MF-Full
10
Floor Number
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Figure B-10 Comparison of maximum story shear envelopes for the moment
frame stick Model MB and the full-building Model MB.
As mentioned earlier, development and calibration of the stick model was performed
using a relatively early version of the Sherman Oaks OpenSees Model MB, which
changed over time. In Figure B-10, a spike in story shear force can be observed in
the basement level for the stick model. This is because early versions of the full-
building Model MB included an artificial phase shift in the ground motions used for
multi-support excitation, which caused a similar spike in shear demand. This phase
shift was removed in subsequent versions of the full-building Model MB, and the
spike no longer appears in the results for the final full-building model depicted in the
figure.
12
10
Floor Number
0 MB-MF-Stick
MB-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
0 MB-MF-Stick
MB-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-11 Comparison of peak floor acceleration envelopes for the moment
frame stick Model MB and the full-building Model MB.
Calibrations were applied to Baseline Model (MB) results only. The calibrated
moment frame stick model was imported into models capturing the foundation
modeling configurations associated with Model 1, Model 2, and Model 4a. Table
B-5 compares the first-, second-, and third-mode periods of the moment frame stick
model and the full-building model for each foundation configuration.
Table B-5 Comparison of Modeled Periods for the Moment Frame Stick
Model and the Full-Building Model for each Alternative
Foundation Configuration
Identified Period Identified Period Identified Period
Mode 1 (sec) Mode 2 (sec) Mode 3 (sec)
Translational Torsional Translational
Model (Global Y-axis) (about Global Z-axis) (Global X-axis)
Periods in each simplified stick model are longer than the periods observed in the
corresponding full-building model. This was attributed to the stick elements being
modeled as linear elements and calibrated with an elastic stiffness that was softer
than the full-building model.
Figures B-12 through B-23 compare response envelopes for maximum displacement,
story drift, story shear, and peak floor acceleration for each alternative foundation
model in each direction. In general, comparisons are similar to trends observed for
Model MB.
12
10
Floor Number
0 1-MF-Stick
1-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
0 1-MF-Stick
1-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Figure B-12 Comparison of maximum displacement envelopes for the moment
frame stick Model 1 and the full-building Model 1.
12 1-MF-Stick
1-MF-Full
10
Story Number
8
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Northridge, Interstory Drifts, Y-Direction
14
12
10
Story Number
0 1-MF-Stick
1-MF-Full
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-13 Comparison of maximum story drift envelopes for the moment frame
stick Model 1 and the full-building Model 1.
12
10
Floor Number
0 1-MF-Stick
1-MF-Full
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Northridge, Y-Direction
14
12
10
Floor Number
0 1-MF-Stick
1-MF-Full
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Figure B-14 Comparison of normalized story shear envelopes for the moment
frame stick Model 1 and the full-building Model 1.
12
10
Floor Number
8
0 1-MF-Stick
1-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
0 1-MF-Stick
1-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-15 Comparison of peak floor acceleration envelopes for the moment
frame stick Model 1 and the full-building Model 1.
12
10
Floor Number
0 2-MF-Stick
2-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
0 2-MF-Stick
2-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Figure B-16 Comparison of maximum displacement envelopes for the moment
frame stick Model 2 and the full-building Model 2.
12 2-MF-Stick
2-MF-Full
10
Story Number
8
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Northridge, Interstory Drifts, Y-Direction
14
12
10
Story Number
0 2-MF-Stick
2-MF-Full
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-17 Comparison of maximum story drift envelopes for the moment frame
stick Model 2 and the full-building Model 2.
12
10
Floor Number
0 2-MF-Stick
2-MF-Full
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Northridge, Y-Direction
14
12
10
Floor Number
0 2-MF-Stick
2-MF-Full
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Figure B-18 Comparison of normalized story shear envelopes for the moment
frame stick Model 2 and the full-building Model 2.
12
10
Floor Number
8
0 2-MF-Stick
2-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
0 2-MF-Stick
2-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-19 Comparison of peak floor acceleration envelopes for the moment
frame stick Model 2 and the full-building Model 2.
12
10
Floor Number
0 4a-MF-Stick
4a-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
0 4a-MF-Stick
4a-MF-Full
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Figure B-20 Comparison of maximum displacement envelopes for the moment
frame stick Model 4a and the full-building Model 4a.
Story Number
8
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Northridge, Interstory Drifts, Y-Direction
14
12
10
Story Number
0 4a-MF-Stick
4a-MF-Full
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-21 Comparison of maximum story drift envelopes for the moment frame
stick Model 4a and the full-building Model 4a.
12
10
Floor Number
0 4a-MF-Stick
4a-MF-Full
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Northridge, Y-Direction
14
12
10
Floor Number
0 4a-MF-Stick
4a-MF-Full
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Figure B-22 Comparison of normalized story shear envelopes for the moment
frame stick Model 4a and the full-building Model 4a.
12
10
Floor Number
8
0 4a-MF-Stick
4a-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
0 4a-MF-Stick
4a-MF-Full
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-23 Comparison of peak floor acceleration envelopes for the moment
frame stick Model 4a and the full-building Model 4a.
Figures B-24 through B-27 overlay response envelopes for maximum displacement,
story drift, story shear, and peak floor acceleration for each alternative foundation
configuration of the moment frame models in each direction. The following trends
were observed:
Displacements for Model 1 are typically the largest, followed by Model 2.
Displacements for Model 4a are very similar to Model MB.
Drifts trends are similar to displacement trends, except that Model MB has larger
drifts than the others models in the basement levels.
Story shear results for Model 1 are the largest. Story shears for Model 2, Model
4a and Model MB are similar.
Peak floor acceleration results for all models are similar in the upper floors.
Model 2 produces significantly larger results in the ground story and basement
levels.
12
10
Floor Number
8
2 1-MF-Stick
2-MF-Stick
0 4a-MF-Stick
MB-MF-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-MF-Stick
2-MF-Stick
0 4a-MF-Stick
MB-MF-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
12
10
Story Number
2 1-MF-Stick
2-MF-Stick
0 4a-MF-Stick
MB-MF-Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-25 Comparison of maximum story drift envelopes for the moment frame
stick Models 1, 2, 4a, and MB.
12
10
Floor Number
8
2 1-MF-Stick
2-MF-Stick
0 4a-MF-Stick
MB-MF-Stick
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Northridge, Y-Direction
14
12 1-MF-Stick
2-MF-Stick
10 4a-MF-Stick
MB-MF-Stick
Floor Number
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Figure B-26 Comparison of normalized story shear envelopes for the moment
frame stick Models 1, 2, 4a, and MB.
12
10
Floor Number
2 1-MF-Stick
2-MF-Stick
0 4a-MF-Stick
MB-MF-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-MF-Stick
2-MF-Stick
0 4a-MF-Stick
MB-MF-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-27 Comparison of peak floor acceleration envelopes for the moment
frame stick Models 1, 2, 4a, and MB.
To investigate the effects of structural stiffness on SSI, moment frame stick model
parameters were scaled to represent the stiffness associated with shear wall systems.
ASCE/SEI 7-10 (ASCE, 2010) and ACI 318-08 (ACI, 2008) were used as references
for determining appropriate scale factors for stiffness. A shear wall stick model was
created by converting the base shear of the full-building Model MB into a design
base shear for an equivalent shear wall using the ratio of response modification
coefficients (R factors) for the special reinforced concrete moment frame and special
reinforced concrete shear wall systems, respectively. Assuming a bearing wall
system, this ratio is 1.6. Based on an amplified base shear, and material properties
predefined in the full-building Model MB, a required shear wall area was determined.
In the case of the core wall configuration, a geometry using two separate interior core
walls, 16 inches thick, was assumed. The cores were configured to match the size
and dimensions of typical elevator, stair, and restroom core areas in typical buildings.
The stiffness was then determined by summing the flexural and shear deflections of a
slender cantilever wall using the following equation:
b b
Vh3 1.2Vh
(B-1)
12 EI Gc A
Where the shear modulus of concrete, Gc, was assumed to be 40% of the elastic
modulus of concrete. Total deflection was primarily controlled by the flexural term,
and the stiffness of the shear wall system was approximately three times the stiffness
of the moment frame system. The stiffness contribution from the gravity system was
also considered. When combined with the stiffness of the shear wall system, a total
stiffness scale factor of 5 was obtained. The core shear wall building model, with a
scale factor of 5, was designated SW. Table B-6 compares the first-, second-, and
third-mode periods of the shear wall stick models for each foundation configuration.
Table B-6 Comparison of Modeled Periods for the Shear Wall Stick Models
for each Alternative Foundation Configuration
Identified Period Identified Period Identified Period
Mode 1 (sec) Mode 2 (sec) Mode 3 (sec)
Translational Torsional Translational
Model (Global Y-axis) (about Global Z-axis) (Global X-axis)
1 Stick (SW) 1.32 N.A. 1.17
2 Stick (SW) 1.38 N.A. 1.19
4a Stick (SW) 1.40 N.A. 1.20
MB Stick (SW) 1.40 N.A. 1.20
Figures B-28 through B-43 overlay response envelopes for maximum displacement,
story drift, story shear, and peak floor acceleration for each alternative foundation
configuration, for each shear wall model, in each direction. The following trends
were observed:
Differences between the models are much more significant in the stiffer shear
wall systems than the moment frame system, demonstrating that SSI effects have
a larger influence on response in stiffer structural systems.
Displacements in Model 1 are typically the smallest. Displacements in Model 4a
are very similar to Model MB.
Drift results are similar to displacement results, except that the Model MB has
larger drifts in the basement levels.
Story shear results vary between the X-direction and Y-direction for all models.
Peak floor accelerations in the superstructure are larger in Model 1 and Model 2.
Model 2 produces significantly larger floor accelerations in the ground story and
basement levels.
12
10
Floor Number
8
2 1-SW -Stick
2-SW -Stick
0 4a-SW -Stick
MB-SW -Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW -Stick
2-SW -Stick
0 4a-SW -Stick
MB-SW -Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
12
10
Story Number
2 1-SW -Stick
2-SW -Stick
0 4a-SW -Stick
MB-SW -Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
12
10
Story Number
2 1-SW -Stick
2-SW -Stick
0 4a-SW -Stick
MB-SW -Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-29 Comparison of maximum story drift envelopes for the SW stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
8
2 1-SW -Stick
2-SW -Stick
0 4a-SW -Stick
MB-SW -Stick
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Northridge, Y-Direction
14
12 1-SW -Stick
2-SW -Stick
10 4a-SW -Stick
MB-SW -Stick
Floor Number
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Figure B-30 Comparison of normalized story shear envelopes for the SW stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
2 1-SW -Stick
2-SW -Stick
0 4a-SW -Stick
MB-SW -Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW -Stick
2-SW -Stick
0 4a-SW -Stick
MB-SW -Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-31 Comparison of peak floor acceleration envelopes for the SW stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
8
2 1-SW 1-Stick
2-SW 1-Stick
0 4a-SW 1-Stick
MB-SW 1-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW 1-Stick
2-SW 1-Stick
0 4a-SW 1-Stick
MB-SW 1-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Figure B-32 Comparison of maximum displacement envelopes for the SW1 stick
Models 1, 2, 4a, and MB.
12
10
Story Number
2 1-SW 1-Stick
2-SW 1-Stick
0 4a-SW 1-Stick
MB-SW 1-Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
12
10
Story Number
2 1-SW 1-Stick
2-SW 1-Stick
0 4a-SW 1-Stick
MB-SW 1-Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-33 Comparison of maximum story drift envelopes for the SW1 stick
Models 1, 2, 4a, and MB.
Floor Number
8
6
4 1-SW 1-Stick
2 2-SW 1-Stick
4a-SW 1-Stick
0
MB-SW 1-Stick
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear
Northridge, Y-Direction
14
12 1-SW 1-Stick
2-SW 1-Stick
10
4a-SW 1-Stick
Floor Number
8 MB-SW 1-Stick
6
4
2
0
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear
Figure B-34 Comparison of normalized story shear envelopes for the SW1 stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
2 1-SW 1-Stick
2-SW 1-Stick
0 4a-SW 1-Stick
MB-SW 1-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW 1-Stick
2-SW 1-Stick
0 4a-SW 1-Stick
MB-SW 1-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-35 Comparison of peak floor acceleration envelopes for the SW1 stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
8
2 1-SW 2-Stick
2-SW 2-Stick
0 4a-SW 2-Stick
MB-SW 2-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW 2-Stick
2-SW 2-Stick
0 4a-SW 2-Stick
MB-SW 2-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Figure B-36 Comparison of maximum displacement envelopes for the SW2 stick
Models 1, 2, 4a, and MB.
12
10
Story Number
2 1-SW 2-Stick
2-SW 2-Stick
0 4a-SW 2-Stick
MB-SW 2-Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
12
10
Story Number
2 1-SW 2-Stick
2-SW 2-Stick
0 4a-SW 2-Stick
MB-SW 2-Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-37 Comparison of maximum story drift envelopes for the SW2 stick
Models 1, 2, 4a, and MB.
Floor Number
8 MB-SW 2-Stick
6
4
2
0
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear
Northridge, Y-Direction
14
1-SW 2-Stick
12 2-SW 2-Stick
10 4a-SW 2-Stick
Floor Number
MB-SW 2-Stick
8
6
4
2
0
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear
Figure B-38 Comparison of normalized story shear envelopes for the SW2 stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
2 1-SW 2-Stick
2-SW 2-Stick
0 4a-SW 2-Stick
MB-SW 2-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW 2-Stick
2-SW 2-Stick
0 4a-SW 2-Stick
MB-SW 2-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-39 Comparison of peak floor acceleration envelopes for the SW2 stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
8
2 1-SW 3-Stick
2-SW 3-Stick
0 4a-SW 3-Stick
MB-SW 3-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW 3-Stick
2-SW 3-Stick
0 4a-SW 3-Stick
MB-SW 3-Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Figure B-40 Comparison of maximum displacement envelopes for the SW3 stick
Models 1, 2, 4a, and MB.
12
10
Story Number
2 1-SW 3-Stick
2-SW 3-Stick
0 4a-SW 3-Stick
MB-SW 3-Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
12
10
Story Number
2 1-SW 3-Stick
2-SW 3-Stick
0 4a-SW 3-Stick
MB-SW 3-Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-41 Comparison of maximum story drift envelopes for the SW3 stick
Models 1, 2, 4a, and MB.
Floor Number
8
6
4 1-SW 3-Stick
2 2-SW 3-Stick
4a-SW 3-Stick
0
MB-SW 3-Stick
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear
Northridge, Y-Direction
14
1-SW 3-Stick
12 2-SW 3-Stick
10 4a-SW 3-Stick
MB-SW 3-Stick
Floor Number
8
6
4
2
0
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Normalized Story Shear
Figure B-42 Comparison of normalized story shear envelopes for the SW3 stick
Models 1, 2, 4a, and MB.
12
10
Floor Number
2 1-SW 3-Stick
2-SW 3-Stick
0 4a-SW 3-Stick
MB-SW 3-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge,X-Direction
Northridge, Y-Direction
14
12
10
Floor Number
2 1-SW 3-Stick
2-SW 3-Stick
0 4a-SW 3-Stick
MB-SW 3-Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-43 Comparison of peak floor acceleration envelopes for the SW3 stick
Models 1, 2, 4a, and MB.
Figures B-44 through B-47 show comparisons between moment frame stick Model
MB and shear wall stick Model MB results for maximum displacement, story drift,
story shear, and peak floor acceleration in each direction. Results clearly indicate a
significant difference in building response between stiffer and more flexible
structural systems.
Northridge, X-Direction
14
12
10
Floor Number
0 MB-MF-Stick
MB-SW -Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
Northridge, Y-Direction
14
12
10
Floor Number
0 MB-MF-Stick
MB-SW -Stick
-2
-40 -30 -20 -10 0 10 20 30 40
Max. Displacement [cm]
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
12
10
Story Number
0 MB-MF-Stick
MB-SW -Stick
-2
-0.015 -0.01 -0.005 0 0.005 0.01 0.015
Max. Interstory Drifts [cm/cm]
Figure B-45 Comparison of maximum story drift envelopes for the moment frame
stick Model MB and the shear wall stick Model MB.
12
10
Floor Number
8
0 MB-MF-Stick
MB-SW -Stick
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Northridge, Y-Direction
14
12 MB-MF-Stick
MB-SW -Stick
10
Floor Number
-2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Normalized Story Shear
Figure B-46 Comparison of normalized story shear envelopes for the moment
frame stick Model MB and the shear wall stick Model MB.
12
10
Floor Number
0 MB-MF-Stick
MB-SW -Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Northridge, Y-Direction
14
12
10
Floor Number
0 MB-MF-Stick
MB-SW -Stick
-2
-1.5 -1 -0.5 0 0.5 1 1.5
Peak Floor Accel. [g]
Figure B-47 Comparison of peak floor acceleration envelopes for the moment
frame stick Model MB and the shear wall stick Model MB.
Renderings of the Walnut Creek building configuration are provided in Figure C-1.
The renderings include cut-away views showing details of the core wall and
perimeter frame configurations, as well as the interior layout of a typical floor level.
Figure C-2 provides plan views of the Walnut Creek building shear wall core at the
foundation, first floor, and typical floor levels.
The building model is composed of the core walls and the perimeter frame. The
framing elements were modeled with force-based beam-column elements using fiber
sections (nonlinearBeamColumn elements in OpenSees). Two methods of modeling
the core wall were considered: (1) use of four-node quadrilateral elements (Quad or
Shell elements in OpenSees); and (2) use of horizontal and vertical boundary frame
elements in combination with diagonal strut elements.
The core walls are 30.5 cm (12 in) thick. To account for the cracked section
properties of the concrete walls, a multiplier on the theoretical uncracked stiffness
must be applied. This multiplier should be dependent on the level of in-situ cracking
of the shear wall and on the intensity of ground shaking, which would dictate the
level of cracking that the wall will experience during an earthquake. A definitive
ratio of cracked to uncracked stiffness in shear walls is not available in the literature.
In the case of flexural stiffness, ASCE/SEI 41-06 (ASCE, 2007) suggests a multiplier
of 0.8 for nominally uncracked walls, and a multiplier of 0.5 for cracked walls.
When using plane-stress Quad elements, the stiffness in the plane of each wall
Poisson’s ratio, (taken as 0.25). The material chosen to represent the concrete was
element is based on the input values of wall dimensions, material behavior, and
The calculated building mass included the mass of all structural elements (beams,
columns, and slabs); 0.5 kPa (10 psf) for partitions; 0.6 kPa (12 psf) for mechanical,
electrical, and plumbing components; and 25% of the design live load taken as 0.6
kPa (12 psf) for a design live load of 2.5 kPa (50 psf). The base of the core wall
extends 4.3 m (14 feet) below the ground level, which includes a 3 m (10 feet)
extension of the core wall, 1.2 m (4 feet) of mat foundation, and sand backfill. The
mass of the core wall (including the sand) is also included in the calculated building
mass. In the model, the mass of the slab and framing elements is assigned to the
nodes, and the mass of each shear wall is distributed to the lower boundary beam.
Damping was modeled as 5% of critical Rayleigh damping, anchored to the first and
third mode periods of the building. The damping was uniformly distributed to all
structural elements, but degrees-of-freedom associated with foundation springs were
excluded. This was necessary to avoid double counting of foundation damping
because the soil-foundation model included dashpots at these degrees-of-freedom.
Recorded motions from the 1989 Loma Prieta earthquake were used to validate the
numerical model. The first step in the validation was a system identification to
Figure C-6 5%-damped elastic response spectra of recorded motions for the
Walnut Creek building in the 1989 Loma Prieta earthquake.
The first two periods identified for the E-W direction were 0.81sec and 0.18sec,
respectively. The effective stiffness of the wall columns and diagonal struts were
then modified so that the numerical model gave an accurate representation of the first
two periods as well as the relative displacement history between the ground floor and
the roof. An effective stiffness equal to 30% of the gross section properties resulted
in the first two modeled periods equal to 0.83sec and 0.23sec, respectively.
Because two vertical instruments were present at the ground floor level, it is possible
to compute base rocking (i.e., the difference between the two vertical recordings
divided by the horizontal distance between them). As shown in Figure C-10,
however, low frequency noise causes base rocking contributions evaluated in this
manner to have unrealistic features. Accordingly, the vertical records were filtered
using high-pass Butterworth filters with a corner frequency of 0.4 Hz.
Figure C-9 Maximum story drift profile for the Walnut Creek building model
subjected to the 1989 Loma Prieta earthquake.
Project Manager
David A. Hutchinson
Buehler & Buehler, Structural Engineers
600 Q Street, Suite 200
Sacramento, California 95811
Workshop Participants
Peter Behnam Robert D. Hanson
KPFF Consulting Engineers 2926 Saklan Indian Drive
6080 Center Drive, Suite 300 Walnut Creek, California 94595
Los Angeles, California 90045
Doug Hohbach
Lauren Carpenter Hohbach-Lewin
WHL Consulting Engineers, Inc. 260 Sheridan Avenue, Suite 150
350 South Figueroa Street, Suite 400 Palo Alto, California 94306
Los Angeles, California 90071
Ben Hushmand
Michael Cochrane Hushmand Associates, Inc.
Weidlinger Associates, Inc. 250 Goddard Road
4551 Glencoe Avenue, Suite 350 Irvine, California 92618
Marina del Rey, California 90292
Peter Lee
John Gavan Skidmore, Owings & Merrill LLP
KPFF Consulting Engineers One Front Street, Suite 2500
6080 Center Drive, Suite 300 San Francisco, California 94111
Los Angeles, California 90045
Anne Lemnitzer
Tony Ghodsi University of California, Los Angeles
Englekirk Structural Engineers Dept. of Civil and Environmental Engineering
925 Fort Stockton Drive, Suite 202 5731 Boelter Hall
San Diego, California 92103 Los Angeles, California 90095