Next Article in Journal
Numerical Analysis of a Self-Acting Gas Bearing Lubricated with a Low-Boiling-Point Medium Using an Advanced Model Based on the Finite Difference Methods and Universal Computational Fluid Dynamics Software
Next Article in Special Issue
Settlement Foundations by Exploring the Collapse of Unsaturated Soils
Previous Article in Journal
EMG Activity of Lower Limb Muscles during Anti-Gravity Treadmill Running with Different Loads and Speeds
Previous Article in Special Issue
Stress State Analysis of the Soil Plug of Open-Ended Piles during Impact Driving Using Particle Flow Code (PFC)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

The Possibility of Detrimental Effects on Soil–Structure Interaction in Seismic Design Due to a Shift in System Frequency

Department of Civil Engineering, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(17), 7519; https://doi.org/10.3390/app14177519
Submission received: 27 April 2024 / Revised: 15 August 2024 / Accepted: 16 August 2024 / Published: 26 August 2024
(This article belongs to the Special Issue Soil-Structure Interaction in Structural and Geotechnical Engineering)

Abstract

:
Soil–structure interaction (SSI) leads to a modification in the dynamic properties of structure, but due to the complexity of analysis, it is traditionally assumed in seismic designs that the structure is fixed-supported on the ground, which brings about potential risks to the seismic performances of structure. The study works on the possibility of SSI having detrimental effects by comparing the dynamic responses of the SSI system to a fixed-base structure, and presents charts for an evaluation of the system frequency of SSI for the purpose of engineering practice. In order to reveal the physical nature, the SSI model is reduced to its simplest form, consisting of a SDOF oscillator, a three-dimensional rectangular foundation, and a multi-layered half-space. The energy dissipation in the soil is achieved by foundation impedances and the substructure method. Previously, the foundation impedances are usually acquired by two-dimensional or axisymmetric three-dimensional models in uniform half-space to avoid the high cost of the more realistic, fully 3D models, while a high-precision indirect boundary element method is employed, combined with the non-singular Green’s functions of distributed loads to calculate the foundation impedances. Although SSI dampens the peak amplitude of structure response in the frequency domain, case studies on four buildings’ responses to 42 earthquakes in the time history show a possibility of 15–20% that SSI amplifies the dynamic responses of structures, such as the maximum and the mean values in the time history, depending on the properties of the structures and the site, as well as the frequency component of incident waves.

1. Introduction

It is assumed in coded seismic designs that structures are fixed-supported at their foundations. In reality, soil experiences obvious flexibility during dynamic excitations, which is responsible for the mechanism of soil–structure interaction (SSI). Since the soil usually acts as a cushion under the structure foundation to decrease the structure responses, SSI was first considered beneficial for the seismic performances of structures several decades ago. Nevertheless, scholars realized, especially after the 2000s, that SSI may sometimes amplify the structure responses compared with a fixed-base assumption, or have detrimental effects on the seismic performances of structures, and that ignoring this mechanism may result in potential risks to seismic designs [1].
In recent years, progress on seismological numerical simulations added more proof regarding this recognition. It was shown, using a probabilistic model based on damage and repair cost, that buildings on very soft soil are extremely likely to incur smaller losses due to SSI, but SSI has an adverse effect on the buildings on moderately soft soil, with a probability up to 0.4 [2]. It was concluded using a statistical study on SSI that the underlying soil may either increase or decrease the structural response depending on the particular site–structure–earthquake scenario [3]. An analysis using damage spectra on SSI reveals that SSI increases the damage index before a threshold period, and the conventional fixed-base model underestimates the damages of buildings with periods less than the threshold [4]. A study using fragility curves shows that neglecting SSI results in nonconservative predicted damage levels of structural and non-structural elements [5].
Although the analysis of SSI is necessary for structural seismic designs, due to its complexity and cost, it currently only applies to a few very important infrastructures, such as nuclear power plants [6], huge dams [7] or long bridges [8,9]. Actually, it was the development of the nuclear power plants in the 1970s that promoted the extensive research on SSI. The analysis of SSI applies to a very small number of civil structures used only for some special purposes; for example, the influence of adjacent topography on buildings [10], the influence of a poroelastic site on deep piled foundations [11,12], and so on. For most civil structures, however, little attention is paid to seismic designs involving SSI. The conventional fixed-base assumption is usually adopted for convenience, and it is only possible to very roughly determine the potentially detrimental effects of SSI by the shear wave velocity, or the period, of the site. Th method, although widely used among civil engineers, is not accurate and cannot be parametrically quantified. For example, many studies using numerical simulations [13,14] and earthquake damage investigations [1,15] have revealed that structures with some of the same features in a certain local area exhibit much more serious damage than other structures without those features.
Nevertheless, there is an ongoing effort to integrate SSI analysis into the framework of coded seismic designs in engineering practice. In ASCE/SCI-7 [16] and NEHRP (FEMA P-1050-1/2015) [17], there are coded provisions in recommended chapters in all versions after the 1990s (v1997, v2000, v2003, v2009, v2015) to involve SSI analysis, which is based on the substructure method and foundation impedance functions. In CECS 160 (general rules for the performance-based seismic design of buildings) [18], similar provisions are also recommended, and the formulas on the evaluation of structure modes are improved. A research report by the National Institute of Standards and Technology (NIST) provides modified coefficient curves, considering SSI based on the concept of traditional coefficient curves in seismic designs [19]. However, these methods are rarely adopted in practice for the following two reasons: The analysis involves the calculation of foundation impedance functions, which are too complex for engineering designs. In addition, the analysis considers only the beneficial effects of SSI and reduces base shear force by at most 30%, which is unsafe or lack of safety reverse.
The foundation of a structure is usually comprised of various shapes, and it forms a group of multi-scatters in infinite half-space, which creates a complex scattering field under seismic excitations. In the past few decades, intensive effort has been put into its numerical simulations, such as raising the incident wave to a more shallow layer to reduce the range of the near-field [20], setting incident waves inside the artificial boundary to reduce the data storage on free-field motion and seismic loads [21], or carrying out a large-scale calculation, from earthquake sources to structure sites, through multi-scale simulation [22], and so on. However, up to now, the numerical simulations on SSI are still very costly and not widely available among civil engineers [23].
This paper investigates the possibility of SSI having detrimental effects using case studies on four buildings whose structural data, site data and earthquake records are available. The scattering field of the half-space is modeled by foundation impedances, which are obtained by a high-precision indirect boundary element method that we proposed in our previous research [24]. Under the earthquake excitations, kinematic interaction and inertial interaction are coupled for SSI, so some studies focus on the former in a half-space with a foundation [25,26], or on the latter with a structure in a uniform half-space [27]. The SSI model, considering the two types of interactions, consists of an oscillator, a 3D rectangular foundation, and a multi-layered half-space. It is shown in the study that the shift in peak frequency between SSI system and fixed-base oscillator leads to the different seismic performances between them, so the paper also proposes a corresponding chart of peak frequency, or system frequency, between SSI and the fixed-base oscillator, to shed some light on SSI analysis in engineering seismic designs. The target parameters in the chart consist of the foundation embedment, soil-layer thickness, length-to-width ratio, and so on, to cater a variety of scenarios in coded seismic designs.

2. Methodology

2.1. Model

The model of the soil–structure interaction system (SSI) consists of a superstructure, a three-dimensional rectangular foundation, and a horizontally multi-layered half-space, as shown in Figure 1. The superstructure is simplified to a single-degree-of-freedom (SDOF) oscillator, and its dynamic properties are characterized by mass Mb, spring coefficient kb, damping coefficient cb and height H, so the fixed-base frequency is ω b = k b / M b and the damping ratio is ξ b = c b / 2 M b ω b . The foundation is represented by a rigid rectangular foundation of mass M0, width a, length b and embedment c. The half-space consists of a soil layer of thickness D on a bedrock, and each sub-layer, as well as the bedrock, is assumed to be elastic and isotropic. The material properties of the jth sub-layer are characterized by shear-wave velocity βj, Poisson’s ratio νj, mass density ρj and damping ratio ξj (j = 1, 2, 3, …), and those of the bedrock are characterized by shear-wave velocity βR, Poisson’s ratio νR, mass density ρR and damping ratio ξR. The SSI system is excited by a train of harmonic SV waves coming from the bedrock, with motion in the XOZ-plane. The circular frequency of the excitation is ω, and the incident angle is θ, measured horizontally.
The deformation of the structure is represented by the relative horizontal translation ub with respect to the z’-axis, which is fixed on the foundation, and the foundation motion is represented by the generalized vector U = { u , φ } T e i ω t , with u of the horizontal translation and φ of the clockwise rotation about the reference point (0, 0, c), and eiωt is the harmonic time factor, omitted hereafter. The vertical displacement, as well as the torsional motion of the SSI system, are not considered in our problem.

2.2. Dynamic Equilibrium

The dynamic equilibrium equations of the SSI are as follows [28,29,30]:
M V ¨ + C V ˙ + K V = M U ¨ *
where
M = [ M b M b M b M b M 0 + M b S 0 + M b M b S 0 + M b I 0 + I b ] ,   C = [ c b 0 0 0 c hh c hm 0 c mh c mm ] ,   K = [ k b 0 0 0 k hh k hm 0 k mh k mm ] ,   M = [ M b M b M 0 + M b S 0 + M b S 0 + M b I 0 + M b ] ,   V = { u b u 0 φ 0 H / a }
The symbol S0 is the area moment of the foundation with respect to the axis z = c, I0 and Ib are the inertia moments of the foundation and the superstructure, respectively, and we have Ib = MbH’ for an oscillator here, for simplicity. The parameter H’ = H + c is the structure height from the structure top to the foundation bottom. Vector U * = { u * , φ * H } T is the foundation input motion, corresponding to the dynamic response of the excavation where the foundation is placed, and U 0 = { u 0 , φ 0 H } T is the relative displacement of the foundation with respect to the foundation input motion [28,29,30], so we have U = U * + U 0 . The physical meaning of U , U * and U 0 in SSI is illustrated in Figure 2.
Foundation impedances are complex-valued and are shown by
K 0 ( ω ) = [ K hh ( ω ) K hm ( ω ) K mh ( ω ) K mm ( ω ) ] = [ k hh ( ω ) k hm ( ω ) k mh ( ω ) k mm ( ω ) ] + i ω [ c hh ( ω ) c hm ( ω ) c mh ( ω ) c mm ( ω ) ]
The real and the imaginary parts act as the spring and dashpot coefficients of the site, with the elements Khh of horizontal impedance, Khm = Kmh of coupled horizontal-rocking impedance and Kmm of rocking impedance. Foundation impedances are obtained using a high-precision indirect boundary element method (IBEM), combined with the non-singular Green’s functions of distributed loads presented in [24]. The foundation impedances in Equation (1) involve both translation and rotation, and the elements are already normalized by H’/a to be unified to the same units.
The boundary element method (BEM) needs element meshing only at boundaries, and it does not need artificial boundaries like the Finite Element Method (FEM) because it uses Green’s functions to automatically satisfy the radiation conditions at infinity. The calculation on Green’s functions is carried out in the frequency domain by Fourier transformation, and its logic is similar to the “displacement method” in the Structural Mechanics textbook. The Green’s functions correspond to the displacement or traction responses of some point in the half-space, subjected to some unit-distributed load, then the foundation impedances are obtained by a series of fictitious distributed loads acting on the foundation-soil boundary. The value of the fictitious loads, as well as the foundation impedances, are all found by the continuous conditions on the boundary.
The dimensionless parameters necessary for the following analysis are:
(1)
The dimensionless mass of the foundation and the structure are M0/Ms and Mb/M0, respectively, where Ms is the mass of soil replaced by the foundation.
(2)
The dimensionless frequency of incident wave is defined as η = λ 1 / 2 a = ω a / π β 1 , where λ 1 is the shear wavelength of the first sublayer. Correspondingly, the dimensionless fixed-base frequency is η b = ω b a / π β 1 , and a general range of ηb in engineering practice is 0–0.4 [28,29,30], and the dimensionless system frequency is η ˜ b = ω ˜ b a / π β 1 .
(3)
The dimensionless size of the foundation and the dimensionless thickness of the soil layer are a/a, b/a, c/a and D/a, respectively.
(4)
The dimensionless mass density and the shear wave velocity of the soil layer are ρR/ρj and βR/βj, respectively.

3. SSI System versus Fixed-Base Oscillator

3.1. Dynamic Properties of SSI

The simplest layered half-space, a single soil layer on bedrock, is taken for the analysis of the dynamic properties of the SSI system. The following parameters of the site are used: ρR/ρj = 1, νR = νj = 1/3, ξR = ξj = 0.02, βR/βj= 2, D/c = 2, (j = 1, 2, 3, …). The foundation embedment is c/a = 0.5, and the foundation mass is M0/Ms = 0.2. Meanwhile, for the structure, we have a damping ratio ξb = 0.02, a height H/a = 2, and mass Mb/M0 = 10 following the recommended relationship Mb/M0 = 4.9H/a in [31]. The length-to-width ratio of the foundation varies, as b/a = 1/4, 1/2, 1, 2 and 4.
Figure 3a shows the spectrum amplitude of structure response | u b | for a fixed-base oscillator and SSI systems with a varying b/a ratio. The x-axis is the dimensionless incident frequency η = 0–0.3, which covers the predominant frequency domain of energy in earthquakes. Three structure frequency values are concerned: ηb = 0.1, 0.15, 0.2. For a fixed-base oscillator, a rigid connection between the structure and the ground is assumed without any energy dissipation through the soil, while for an SSI system, the soil behaves like a cushion under the structure and adds additional flexibility to the system, so the system frequency of the SSI is obviously lower than the resonant frequency of the corresponding fixed-base oscillator. Further, as the length-to-width ratio b/a decreases, the peak has a higher frequency and a higher amplitude, approaching the fixed-base case.
Figure 3b shows the dynamic properties (system frequency, system damping, peak amplitude) of the structure response | u b | , as well as the varying b/a ratios. The x-axis is the dimensionless structure frequency ηb = 0–3. The system frequency η ˜ b corresponds to the peak frequency of | u b | , or the eigen-value of Equation (1) [28,29,30], and here we use the ratio η ˜ b / η b to indicate a shift in system frequency with respect to the structure frequency η b . Although the spectrum shape is similar between the fixed-base oscillator and the SSI system with various length-to-width ratios, for the values of b/a and ηb concerned in the figure, the system frequency shifts up to 50% from the structure frequency. For example, it shifts by 12% for a framed structure of width 2a = 40 m, length 2b = 40 m and fixed-base frequency ω b = 6.28 Hz on a moderate rigid site (β1 = 300 m/s), and it shifts by 16% for the same structure on a moderate soft site (β1 = 250 m/s). In fact, the shift in peak frequency between the SSI system and the fixed-base oscillator leads to a completely different correspondence with respect to the frequency of excitations, which is further investigated in Figure 4.
The second sub-figure of Figure 3b describes the system damping ξ ˜ b which is obtained by the half-power method from | u b | , and we use the ratio ξ ˜ b / ξ b to show a comparison between the SSI system and the fixed-base oscillator. The ratio ξ ˜ b / ξ b > 1 indicates that the system damping is obviously larger than the structure damping, while the third sub-figure is the peak amplitude | u b max | between the SSI system and the fixed-base oscillator. Correspondingly, the peak amplitude in the frequency domain of the SSI system is smaller than the fixed-base oscillator, which gives | u b_SSI max / u b_fixed-base max | < 1 .

3.2. Shift in System Frequency

Figure 4 shows a spectrum correspondence between the structure response | u b | (the line style follows Figure 3) and the incident wave. Due to the shift in peak frequency, the SSI system demonstrates different seismic performances compared to the fixed-base oscillator. When the predominant frequency of the incident wave matches the peak frequency of the SSI system (acceleration record 1 in Figure 4), the dynamic responses of the SSI system are probably larger than the fixed-base oscillator, which probably has detrimental SSI effects on the seismic performances of a structure, so the seismic designs on the assumption of a fixed-base oscillator are probably unsafe. Whereas, for a case in which the predominant frequency matches that of the fixed-base oscillator (acceleration record 2 in Figure 4), the responses of the former are probably smaller than the latter, which probably has beneficial SSI effects, so the seismic designs on the assumption are probably too conservative. Certainly, for a case in which the predominant frequency of the incident wave is not obvious, it also probably results in beneficial SSI effects, since the peak amplitude of the SSI system is always smaller than that of the corresponding fixed-base oscillator. In summary, a smaller peak amplitude of SSI system in the frequency domain does not necessarily mean that the dynamic responses of SSI in the time domain is always smaller than the fixed-base structure. SSI effects are sometimes detrimental to the seismic performances of structures, depending on the dynamic properties of both the structure itself and the soil.
Figure 5 presents the ratio η ˜ b / η b (system frequency/fixed-base frequency) as that which is shown in Figure 3b, but for different structure and site parameters. The structure mass is for three values Mb/M0 = 5, 10, 15, and two kinds of half-space are concerned, βR/βj = 1 and 2, with the other parameters the same as Figure 3. βR/βj = 1 corresponds to a uniform half-space. The shift in system frequency in uniform half-space is larger than that in layered half-space, since the bedrock adds rigidity to the site. Also, the shift in system frequency for a high structure is larger than that for a short structure, since a high structure is more flexible according to our parameters.

4. Detrimental Effects of SSI

4.1. Data on Structures and Sites

A further study on the time domain was performed on four buildings in California, USA, to investigate the possibility of SSI having detrimental effects. The parameters were selected from the published literature: the Hollywood Storage Building (California Strong Motion Instrumentation Program (CSMIP) station No. 24236) [32], the Millikan Library (National Strong Motion Project (NSMP) No. 5407) [33], a building in Sherman Oaks (CSMIP station No. 24322) and a building in Walnut Creek (CSMIP station No. 58364) [19]. These buildings have been used for strong-motion observations for decades, and their structural and site data are available, as well as earthquake records, which is helpful to establish proper SSI models. The details of the four buildings and the stratifications of the site can be found in our previous paper [34], and we list key information on the buildings in Table 1 and on the sites in Table 2, which are used to establish the SSI model. The influence of SSI is probably not taken into account in the seismic designs of the buildings. The fixed-base frequency ωb/2π for the four buildings ranges from 0.83 Hz to 2.16 Hz, and the shear-wave velocity of the soil at the top of the site ranges from 184.7 m/s (soft soil) to 336 m/s (moderate rigid soil). A rough formula H e = H / 3 suggested in [35] is used as the effective height of the structures, since the structures are simplified to the equivalent SDOF oscillator.

4.2. Case Studies

Figure 6a,b are the dynamic responses of the Hollywood Storage Building for an assumed SSI system and an assumed fixed-base structure, respectively. The excitation, as shown at the top of the figure, is the acceleration record obtained at the ground floor of the building in the NS direction (a component of the horizontal direction) during the Northridge earthquake in 1994, and it is normalized to an acceleration of a peak value of 0.1 g. The record can be found on the USGS website https://strongmotioncenter.org/ (accessed on 28 January 2019). It is also assumed that the excitation is a train of SV waves coming from the bedrock with an incident angle θ = 90°. The acceleration responses of both the SSI system and the fixed-base structure are shown in the time history u ¨ b ( t ) (left column) and the response spectrum of the time history R S ( u ¨ b ) (right column).
The responses in the time history in Figure 6 were gathered according to the following steps: First, the time history of an earthquake acceleration was converted into the frequency domain at 4096 + 1 equally spaced points from 0 to 50 Hz using Fourier transformation to acquire the spectrum of earthquake acceleration. The corresponding values on the spectrum of the structure response ub were calculated at these frequency points as well. Thereafter, the spectrum of the structure response was multiplied by the spectrum of the earthquake acceleration at the corresponding frequency points to produce the amplified spectrum of the structure. Finally, the response in the frequency domain was converted to the time domain using inverse fast Fourier transformation.
In the time history graph (left column), for the time window 0–40 s, the value denoted by an asterisk “*” is the maximum value during the time history, and the value presented in the bottom right corner is the mean value. In the response spectrum graph, the value denoted by an asterisk shows the same. It can be noticed that the maximum value of the SSI system (0.33 g) is 57% larger than that of the fixed-base oscillator (0.21 g), and the mean value of the SSI system (4.33 × 10−2 g) is 48% larger than that of the fixed-base oscillator (2.92 × 10−2 g). Also, the maximum value in the response spectrum of the SSI system (0.96 g) is 25% larger than that of the fixed-base oscillator (0.80 g). We list these three values (the maximum and the mean values in the time history, and the maximum value in the response spectrum) for the four buildings to 42 earthquake events in Table 3 (XOZ plane) and Table 4 (YOZ plane) for both the SSI system and the fixed-base oscillator. The excitations are the acceleration records from the ground floor of each building for each earthquake event (https://strongmotioncenter.org/, accessed on 28 January 2019).
Figure 7 illustrates the ratio of the three responses (maximum response in the time history, mean response in the time history, and maximum response spectrum of structure acceleration response) of the SSI system to those of the fixed-base oscillator. The x-axis corresponds to each event in Table 3 and Table 4, and the y-axis is the percentage value. The symbol “FB” represents a “fixed-base oscillator”. A value larger than 100% on the y-axis indicates that the dynamic response of the SSI system is larger than that of the corresponding fixed-base oscillator, or that SSI amplifies the dynamic response of the structure and has a detrimental effect on the seismic designs. As shown in the figure, about 15–20% of the 42 selected earthquakes on four buildings’ events had detrimental effects. The SSI mechanism amplifies the structure responses by as much as 1.5 times for a few cases. This justifies the idea that a shift on the spectrum may lead to considerably different responses in the time domain, depending on the frequency component of the incident waves and the parameters of the structure as well as the soil.

4.3. Evaluation of System Frequency

Involving SSI in the seismic designs could result in neutral/negligible, beneficial or detrimental effects on structure safety. Due to the uncertainties of seismic motions, there is little necessity for the seismic analysis on SSI having beneficial effects, while it is necessary to develop a method that determines the detrimental effects to ensure the seismic safety of structures. However, seismic hazard assessments on SSI are costly, even for linear problems. For a seismic engineering analysis, an efficient method that is applicable in coded standard procedures is needed. For an engineering seismic analysis, it needs efficient method that is applicable in coded standard procedures. Since the structure period (frequency) is the most important parameter in the seismic designs, it makes sense to predict on the shift of structure period (frequency) in engineering practice. In the following, we describe a method, using Figure 5 to predict the system frequency of the Millikan Library.
The parameters of the structure are a = 11.65 m, b = 12.55 m, c = 4.3 m, He = 29.3 m, M0 = 1.43 × 106 kg, Mb = 1.07 × 107 kg, so we have b/a = 1.08, c/a = 0.37, He/a = 2.52, Mb/M0 = 7.48. The fifth layer (β5 = 609.6 m/s) is assumed as bedrock, and we simplify the multi-layered half-space to a uniform half-space by the equivalent shear wave velocity of the site
β eq = 5.49 + 4.26 + 3.66 + 6.71 5.49 298.7 + 4.26 387.1 + 3.66 454.2 + 6.71 487.7 = 393.0   m / s  
so we have the dimensionless frequency in the NS direction in the equivalent uniform half-space as
η b = ( 2.16 Hz × 2 π ) × 11.65 m 393.0 m / s × π = 0.128
We find the curve in parameters b/a = 1, He/a = 2, Mb/M0 = 10 in uniform half-space at the x-axis of ηb = 0.128, then we measure the value of y-axis of η ˜ b / η b = 0.85 , to evaluate the system frequency ω ˜ b / 2 π = 1.836 Hz . The evaluation of the system frequency using Figure 5 is very close to the true value 1.84 Hz. While in an EW direction, it gives η b = 0.075 , and we measure η ˜ b / η b = 0.94 from y-axis, to evaluate the system frequency ω ˜ b / 2 π = 1.18 Hz , which is also very close to the true value 1.20 Hz.
The parameters of the other three buildings are a little far from the parameter combinations described in Figure 5, so it is difficult to use Figure 5 to predict their system frequency. It is necessary to find a general method capable of covering a variety of structure and site parameters. However, the SSI system is multi-parametric, and dozens of target parameters are necessary to obtain the global function, which maps subdomain parameters to the system frequency. At last, we describe a possible way, based on an artificial neural network, to realize the evaluation of system frequency. The artificial neural network is an algorithm that fits a complex mapping of multiple inputs and outputs. Through appropriate machine learning, this algorithm could theoretically converge to any function.
The system frequency of SSI depends on the properties of the structure itself as well as the underlying soil. The soil that can determine the system frequency is the surface layer across dozens or hundreds of meters, the constructions and the parameters of which are statistically consistent in adjacent areas with the same geological environment. Moreover, the structural parameters follow general principles in structural designs and engineering practices. Therefore, it is feasible to create a neural network mapping the parameters of the two substructures (structure and site) to the SSI system, and the method is applicable to cities whose geological stratum data have been accumulated plentifully over the past decades, in order to obtain a numerical model on seismic hazard assessment involving SSI effects. This will be the focus of our future work.

5. Conclusions

To date, it has been assumed in most seismic designs that the structures are fixed-supported on the foundation, or that the underlying soil is completely rigid. In fact, the mechanism of the soil–structure interaction introduces considerable differences to the seismic performance of a structure. This study investigated the possibility of SSI amplifying the dynamic responses of a structure, or SSI having detrimental effects on the seismic performance of a structure, by comparing the dynamic responses of the SSI system and a fixed-base structure. The SSI model consists of an oscillator, a 3D rectangular foundation, and a multi-layered half-space. The foundation impedance functions are obtained by the high-precision indirect boundary element method from our previous studies.
The soil behaves like a cushion under the structure and adds additional flexibility, so from the viewpoint of the spectrum in the frequency domain, the system frequency of SSI is lower than the resonant frequency of the fixed-base structure, the system damping of SSI is larger than the structure damping, and the system amplitude of SSI is smaller than the peak amplitude of fixed-base structure. However, a smaller peak amplitude of the SSI system in the frequency domain does not necessarily mean that the dynamic responses of SSI in the time domain are always smaller than the fixed-base structure. Case studies on four buildings to 42 earthquakes in the time domain show a 15–20% possibility that SSI would have detrimental effects on the seismic performances of structure, or SSI amplifying the dynamic responses of structure.
This is because there is a shift in peak frequency between the SSI system and the fixed-base structure. When the predominant frequency of the incident wave matches the peak frequency of the SSI system, the dynamic responses of the SSI system are probably larger than the fixed-base oscillator, which probably has detrimental SSI effects on the seismic performance of a structure. The influence of SSI on the seismic performance of a structure depends on the dynamic properties of both the structure itself and the soil.
It should be noted that the possibility of SSI having detrimental effects is based on the four selected buildings as well as the earthquake records obtained from the buildings, so the database is not very plentiful, and the result is not very generalized. More cases are needed in the future to enrich the topic and the conclusion. It should also be noted that both the soil and the structure are likely to experience some degree of nonlinearity during strong earthquake events, and the detrimental effects of SSI under nonlinear conditions are the focus of our future work.

Author Contributions

Conceptualization, J.F.; Methodology, J.F.; Software, J.F.; Formal analysis, J.F.; Investigation, J.F. and Y.L.; Writing—original draft, W.T.; Writing—review & editing, W.T.; Funding acquisition, W.T., J.F. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Provincial Natural Science Foundation of Shaanxi (2022JM-226, 2022JQ-513), National Natural Science Foundation of China (52278218), Fundamental Research Funds for Central Universities (G2021KY05103).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors are grateful to one of the anonymous reviewers for his/her constructive comments, which lead to improvements in the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mylonakis, G.; Gazetas, C. Seismic soil-structure interaction: Beneficial or detrimental? J. Earthq. Eng. 2000, 4, 277–301. [Google Scholar] [CrossRef]
  2. Khosravikia, F.; Mahsuli, M.; Ghannad, M.A. The effect of soil–structure interaction on the seismic risk to buildings. Bull. Earthq. Eng. 2018, 16, 3653–3673. [Google Scholar] [CrossRef]
  3. Moghaddasi, M.; Cubrinovski, M.; Chase, J.G.; Pampanin, S.; Carr, A. Probabilistic evaluation of soil- foundation-structure interaction effects on seismic structural response. Earthq. Eng. Struct. Dyn. 2011, 40, 135–154. [Google Scholar] [CrossRef]
  4. Nakhaei, M.; Ghannad, M.A. The effect of soil-structure interaction on damage index of buildings. Eng. Struct. 2008, 30, 1491–1499. [Google Scholar] [CrossRef]
  5. Forcellini, D. Seismic fragility of tall buildings considering soil structure interaction (SSI) effects. Structures 2022, 45, 999–1011. [Google Scholar] [CrossRef]
  6. Dong, V.N.; Dookie, K.; Duan, D.N. Nonlinear seismic soil-structure interaction analysis of nuclear reactor building considering the effect of earthquake frequency content. Structures 2020, 26, 901–914. [Google Scholar]
  7. Dong, V.N.; Dookie, K.; Chunse, P.; Byounghan, C. Seismic Soil–Structure Interaction Analysis of Concrete Gravity Dam Using Perfectly Matched Discrete Layers with Analytical Wavelengths. J. Earthq. Eng. 2019, 25, 1657–1678. [Google Scholar]
  8. Zhang, N.; Wang, L.; Zhang, Y.; Lu, H.; Dai, D.; Cai, G.; Gao, Y. Effect of a V-shaped canyon on the seismic response of a bridge under oblique incident SH waves. Earthq. Eng. Struct. Dyn. 2023, 53, 496–514. [Google Scholar] [CrossRef]
  9. Zheng, S.-X.; Jia, H.-Y.; Zhao, C.-H.; Qu, H.-L.; Shi, X.-L. Seismic response analysis of long-span and asymmetrical suspension bridges subjected to near-fault ground motion. Eng. Fail. Anal. 2020, 115, 104615. [Google Scholar] [CrossRef]
  10. Oztoprak, S.; Cinicioglu, S.F.; Oztorun, N.K.; Alhan, C. Impact of neighbouring deep excavation on high-rise sun plaza building and its surrounding. Eng. Fail. Anal. 2020, 111, 104495. [Google Scholar] [CrossRef]
  11. Cui, C.; Zhang, S.; Chapman, D.; Meng, K. Dynamic impedance of a floating pile embedded in poro-visco-elastic soils subjected to vertical harmonic loads. Geomech. Eng. 2018, 15, 793–803. [Google Scholar]
  12. Cui, C.Y.; Meng, K.; Wu, Y.J.; Chapman, D.; Liang, Z.M. Dynamic response of pipe pile embedded in layered visco-elastic media with radial inhomogeneity under vertical excitation. Geomech. Eng. 2018, 16, 609–618. [Google Scholar]
  13. Renzi, S.; Madiai, C.; Vannucchi, G. A simplified empirical method for assessing seismic soil-structure interaction effects on ordinary shear-type buildings. Soil Dyn. Earthq. Eng. 2013, 55, 100–107. [Google Scholar] [CrossRef]
  14. Fu, J.; Todorovska, M.I.; Liang, J. Correction factors for SSI effects predicted by simplified models: 2D versus 3D rectangular embedded foundations. Earthq. Eng. Struct. Dyn. 2018, 47, 1963–1983. [Google Scholar] [CrossRef]
  15. Cinna, L. Mexico, San Francisco, Los Angeles and Kobe: What next? Nat. Hazards 1997, 16, 287–296. [Google Scholar]
  16. ASCE/SEI 7-16; Minimum Design Loads for Buildings and Other Structures. American Society of Civil Engineers: Reston, VA, USA, 2016.
  17. FEMA P-1050-1/2015; National Earthquake Hazard Reduction Program (NEHRP) Recommended Seismic Provisions for New Buildings and Other Structures: Part 2 Commentary. Building Seismic Safety Council (BSSC): Washington, DC, USA, 2015.
  18. CECS 160: 2004; General Rules for Performance-Based Seismic Design of Buildings (Trial Version). China Association for Engineering Construction Standardization (CECS): Beijing, China, 2004.
  19. National Institute of Standards and Technology. Soil-Structure Interaction for Building Structures; Report No. NIST GCR 12-917-21; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2012.
  20. Zhao, M.; Liu, P.; Zhang, J.; Zhang, G.; Gao, Z.; Du, X. A direct time-domain FEM-SBFEM-SBPML method for soil-structure interaction analysis using quadtree mesh. Comput. Geotech. 2023, 161, 105597. [Google Scholar] [CrossRef]
  21. Liu, J.; Bao, X.; Wang, D.; Tan, H.; Li, S. The internal substructue method for seismic wave input in 3D dynamic soil-structure interaction analysis. Soil Dyn. Earthq. Eng. 2019, 127, 105847. [Google Scholar] [CrossRef]
  22. Zhang, L.; Wang, J.T.; Xu, Y.J.; He, C.H.; Zhang, C.H. A Procedure for 3D Seismic Simulation from Rupture to Structures by Coupling SEM and FEM. Bull. Seismol. Soc. Am. 2020, 110, 1134–1148. [Google Scholar] [CrossRef]
  23. Zhang, X.; Far, H. Seismic behaviour of high-rise frame-core tube structures considering dynamic soil-structure interaction. Bull. Seismol. Soc. Am. 2022, 20, 5073–5105. [Google Scholar] [CrossRef]
  24. Fu, J.; Liang, J.; Ba, Z. Non-singular boundary element method on impedances of three-dimensional rectangular foundations. Eng. Anal. Bound. Elem. 2019, 99, 100–110. [Google Scholar] [CrossRef]
  25. Luis, A.P.; Suarez, A.; Aznarez, J.J.; Maeso, O. Kinematic internal forces in deep foundations with inclined piles. Earthq. Eng. Struct. Dyn. 2015, 44, 2129–2135. [Google Scholar]
  26. Chen, S.S.; Liao, K.H.; Shi, J.Y. A dimensionless parametric study for forced vibrations of foundation-soil systems. Comput. Geotech. 2016, 76, 184–193. [Google Scholar] [CrossRef]
  27. Ghannad, M.A.; Jahankhah, H. Site-dependent strength reduction factors for soil-structure systems. Soil Dyn. Earthq. Eng. 2007, 27, 99–110. [Google Scholar] [CrossRef]
  28. Luco, J.E.; Lanzi, A. Approximate soil-structure interaction analysis by a perturbation approach: The case of stiff soils. Soil Dyn. Earthq. Eng. 2013, 51, 97–110. [Google Scholar] [CrossRef]
  29. Lanzi, A.; Luco, J.E. Approximate soil-structure interaction analysis by a perturbation approach: The case of soft soils. Soil Dyn. Earthq. Eng. 2014, 66, 415–428. [Google Scholar] [CrossRef]
  30. Fu, J.; Liang, J.; Todorovska, M.I.; Trifunac, M.D. Soil-structure system frequency and damping: Estimation from eigenvalues and results for a 2D model in layered half-space. Earthq. Eng. Struct. Dyn. 2018, 47, 2055–2075. [Google Scholar] [CrossRef]
  31. Luco, J.E. Soil-structure interaction and identification of structural models. In Proceedings of the 2nd ASCE Conference on Civil Rngineering and Nuclear Power, Knoxville, TN, USA, 15–17 September 1980; Geotechnical Topics. ASCE: New York, NY, USA, 1980; Volume II, pp. 10.1–10.30. [Google Scholar]
  32. Trifunac, M.D.; Hao, T.Y.; Todorovska, M.I. Response of a 14-Story Reinforced Concrete Structure to Nine Earthquakes: 61 Years of Observation in the Hollywood Storage Building; Report CE 01-02; Department of Civil Engineering, University of Southern California: Los Angeles, CA, USA, 2001. [Google Scholar]
  33. Luco, J.E.; Wong, H.L.; Trifunac, M.D. Soil-Structure Interaction Effects on Forced Vibration Tests; Report 86-05; University of Southern California, Department of Civil Engineering: Los Angeles, CA, USA, 1986. [Google Scholar]
  34. Fu, J.; Liang, J. Performance degradation of tuned-mass-dampers arising from ignoring soil-structure interaction effects. Soil Dyn. Earthq. Eng. 2019, 125, 105701. [Google Scholar]
  35. Todorovska, M.I.; Trifunac, M.D. The system damping, the system frequency and the system response peak amplitudes during in-plane building-soil interaction. Earthq. Eng. Structrual Dyn. 1992, 21, 127–144. [Google Scholar] [CrossRef]
Figure 1. The SSI model consisting of a superstructure, a three-dimensional rectangular foundation, and a horizontally multi-layered half-space.
Figure 1. The SSI model consisting of a superstructure, a three-dimensional rectangular foundation, and a horizontally multi-layered half-space.
Applsci 14 07519 g001
Figure 2. The foundation motion U , foundation input motion U * and foundation relative motion U 0 .
Figure 2. The foundation motion U , foundation input motion U * and foundation relative motion U 0 .
Applsci 14 07519 g002
Figure 3. Dynamic responses (a) and dynamic properties (b) of the SSI system (βR/βj = 2, D/a = 2, c/a = 0.5, Mb/M0 = 10, H/a = 2).
Figure 3. Dynamic responses (a) and dynamic properties (b) of the SSI system (βR/βj = 2, D/a = 2, c/a = 0.5, Mb/M0 = 10, H/a = 2).
Applsci 14 07519 g003
Figure 4. Peak frequency of the structure response and predominant frequency of excitations.
Figure 4. Peak frequency of the structure response and predominant frequency of excitations.
Applsci 14 07519 g004
Figure 5. Shift in system frequency in uniform half-space and a layered half-space (D/c = 2, c/a = 0.5, ρR/ρj = 1, νR = νj = 1/3, ξR = ξj = 0.02).
Figure 5. Shift in system frequency in uniform half-space and a layered half-space (D/c = 2, c/a = 0.5, ρR/ρj = 1, νR = νj = 1/3, ξR = ξj = 0.02).
Applsci 14 07519 g005
Figure 6. (Top) The acceleration record obtained at the ground floor of the Hollywood Storage Building in a NS direction during the Northridge earthquake in 1994; (a) the structural response of an SSI system for time history (left) and for response spectrum (right); (b) the structural response of a fixed-base oscillator for time history (left) and for response spectrum (right).
Figure 6. (Top) The acceleration record obtained at the ground floor of the Hollywood Storage Building in a NS direction during the Northridge earthquake in 1994; (a) the structural response of an SSI system for time history (left) and for response spectrum (right); (b) the structural response of a fixed-base oscillator for time history (left) and for response spectrum (right).
Applsci 14 07519 g006
Figure 7. Ratio of dynamic responses between the SSI system and the fixed-base oscillator (shown by “FB”) (red circle: Hollywood Storage Building, gray star: Millikan Library, blue triangle: Sherman Oaks Building, green rectangle: Walnut Creek Building).
Figure 7. Ratio of dynamic responses between the SSI system and the fixed-base oscillator (shown by “FB”) (red circle: Hollywood Storage Building, gray star: Millikan Library, blue triangle: Sherman Oaks Building, green rectangle: Walnut Creek Building).
Applsci 14 07519 g007
Table 1. Buildings for case studies.
Table 1. Buildings for case studies.
ParametersHollywood Storage BuildingMillikan LibrarySherman Oaks BuildingWalnut Creek Building
a (short side)7.77 m11.65 m18.15 m15.9 m
b (long side)33.15 m12.55 m28.8 m22.6 m
c8.74 m4.30 m16.10 m4.20 m
H45.6 m43.9 m50.0 m39.2 m
He26.3 m25.3 m28.9 m22.6 m
M0 (kg)1.38 × 107 kg1.43 × 106 kg4.57 × 107 kg1.43 × 106 kg
Mb1.17 × 107 kg1.07 × 107 kg1.66 × 107 kg1.07 × 107 kg
Mb/M00.857.480.247.48
ξ b 2%2%2%2%
NSb/a4.270.931.591.42
c/a1.120.340.890.26
He/a3.382.171.591.42
ω b a / π β 1 0.0700.1681.100.14
ω b / 2 π 0.83 Hz2.16 Hz1.54 Hz1.52 Hz
ω ˜ b / 2 π 0.78 Hz1.86 Hz1.38 Hz1.40 Hz
ω ˜ b / ω b ( η ˜ b / η b ) 0.930.840.900.92
EWa/b0.231.080.630.70
c/b0.260.370.560.186
He/b0.792.021.001.00
ω b a / π β 1 0.720.110.610.28
ω b / 2 π 2.00 Hz1.26 Hz1.69 Hz2.10 Hz
ω ˜ b / 2 π 1.78 Hz1.20 Hz1.55 Hz1.88 Hz
ω ˜ b / ω b ( η ˜ b / η b ) 0.890.940.920.89
Table 2. Site parameters for case studies.
Table 2. Site parameters for case studies.
S-Wave Velocity
(m/s)
Thickness
(m)
Mass Density
(kg/m3)
Damping Ratio
Hollywood Storage Building
1184.715.22050.42%
2362.715.22082.42%
3554.730.52082.42%
4624.841.52082.42%
Bedrock1045.52082.42%
Millikan Library
1298.75.491846.92%
2387.14.261846.92%
3454.23.661846.92%
4487.76.711846.92%
5609.682.291846.92%
6762.016.161846.92%
Bedrock944.81846.92%
Sherman Oaks Building
11602.7 2040.82%
22056.72040.82%
326013.72040.82%
433023.82040.82%
Bedrock5142040.82%
Walnut Creek Building
13365.01846.92%
Bedrock4341846.92%
Table 3. Dynamic responses of structures in the time domain of target buildings in the XOZ-plane.
Table 3. Dynamic responses of structures in the time domain of target buildings in the XOZ-plane.
EarthquakeYear Time   History   u ¨ b ( t ) R S ( u ¨ b )
Maximum (g)Mean (g) ×10−2Maximum (g)
SSI SystemFixed-BaseSSI SystemFixed-BaseSSI SystemFixed-Base
Hollywood Storage Building (NS)
Kern County19520.430.287.725.251.671.40
San Fernando 19710.370.155.523.121.170.85
Whittier Narrow 19870.320.154.924.340.970.90
Northridge 19940.330.214.332.920.960.80
Chino Hills 20080.410.187.242.751.730.68
Encino20140.270.111.490.610.570.23
Millikan Library (NS)
Lytle Creek19700.410.3311.409.833.442.56
San Fernando19710.130.213.415.610.571.10
Whittier Narrow19870.160.373.127.100.902.66
Yorba Linda20020.130.262.915.760.691.96
San Simeon20030.120.133.453.490.871.03
Sherman Oaks Building (NS)
Lander19920.210.285.646.551.301.66
Whittier Narrow19870.130.151.381.520.340.33
Northridge19940.200.222.433.700.781.51
Chatsworth20070.0950.1021.001.040.330.47
Chino Hills20080.200.162.442.570.731.01
Encino20140.0610.120.0530.060.180.33
Walnut Creek Building (EW)
Livermore19800.330.344.946.201.321.48
Livermore Aftershock19800.270.394.337.841.502.37
Loma Prieta19890.240.334.545.831.191.68
Alamo20080.100.121.251.770.300.45
Table 4. Dynamic responses of structures in the time domain of target buildings in the YOZ-plane.
Table 4. Dynamic responses of structures in the time domain of target buildings in the YOZ-plane.
EarthquakeYear Time   History   u ¨ b ( t ) R S ( u ¨ b )
Maximum (g)Mean (g) ×10−2Maximum (g)
SSI SystemFixed-BaseSSI SystemFixed-BaseSSI SystemFixed-Base
Hollywood Storage Building (EW)
Kern County19520.440.327.847.022.392.50
San Fernando 19710.380.194.013.171.451.36
Whittier Narrow 19870.490.538.6015.422.094.69
Northridge 19940.450.326.604.153.111.95
Chino Hills 20080.350.415.526.451.943.14
Encino20140.190.252.883.670.881.48
Millikan Library (EW)
Lytle Creek19700.390.499.4715.442.373.77
San Fernando19710.200.196.657.741.351.53
Whittier Narrow19870.270.345.536.361.301.80
Yorba Linda20020.1070.1052.623.540.560.70
San Simeon20030.440.5511.7414.903.134.19
Sherman Oaks Building (EW)
Lander19920.190.275.208.101.201.99
Whittier Narrow19870.200.262.955.500.871.64
Northridge19940.1870.1902.913.111.131.15
Chatsworth20070.0550.1060.0600.0870.190.43
Chino Hills20080.200.243.585.181.461.63
Encino20140.0610.120.0530.060.180.33
Walnut Creek Building (NS)
Livermore19800.280.376.258.172.072.42
Livermore Aftershock19800.140.172.753.270.851.01
Loma Prieta19890.140.203.092.940.961.33
Alamo20080.130.191.112.150.360.59
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tao, W.; Fu, J.; Li, Y. The Possibility of Detrimental Effects on Soil–Structure Interaction in Seismic Design Due to a Shift in System Frequency. Appl. Sci. 2024, 14, 7519. https://doi.org/10.3390/app14177519

AMA Style

Tao W, Fu J, Li Y. The Possibility of Detrimental Effects on Soil–Structure Interaction in Seismic Design Due to a Shift in System Frequency. Applied Sciences. 2024; 14(17):7519. https://doi.org/10.3390/app14177519

Chicago/Turabian Style

Tao, Weifeng, Jia Fu, and Yugang Li. 2024. "The Possibility of Detrimental Effects on Soil–Structure Interaction in Seismic Design Due to a Shift in System Frequency" Applied Sciences 14, no. 17: 7519. https://doi.org/10.3390/app14177519

APA Style

Tao, W., Fu, J., & Li, Y. (2024). The Possibility of Detrimental Effects on Soil–Structure Interaction in Seismic Design Due to a Shift in System Frequency. Applied Sciences, 14(17), 7519. https://doi.org/10.3390/app14177519

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop