Effect of Soil–Pile–Structure Interaction on Seismic Design of Tall and Massive Buildings Through Case Studies

In this study, the effects of Soil-Pile-Structure Interaction (SPSI) forces on the seismic response of two tall and massive building projects were investigated. The first case study is an emergency hospital building located in south of Switzerland which is an existing reinforced concrete construction with concrete columns. Due to the importance of the building and its potential danger in failure cases of the load-bearing structure, seismic analysis was first carried out by using a 3D finite element model for the entire structure and a new FORTRAN program based on the enhanced cone model for the soil under the structure. The implementation of simulation was conducted in two different conditions, namely fixed-base behavior and soil–structure interaction, which are considered. Based on the results obtained by using the cone frustum approach, the period of vibration in the case of concentrated piles with the correction factor for pile-group action, exhibiting the value of (1.2 s), is increased up to 42% compared with the case of fixed-base (0.7 s). In the second phase, the dynamic responses of the Green Building at the MIT campus in Cambridge during three ground motions were studied through both approaches. It has shown that the first period of vibration in the case of the MIT Science building increased by about 11% (from 1.37 to 1.53 s) with respect to the SPSI effect. Besides, a comparison was made between the structural responses of the obtained cone method and 3D finite element simulation in ABAQUS which indicated that these results were in good agreement with the direct results. It was concluded that the enhanced cone model as a convenient, fast, and rather accurate method can be applied for foundation vibration and dynamic soil–structure interaction analysis in practical engineering projects whenever possible.


Introduction
The Soil-Pile-Structure Interaction (SPSI) plays an important role in assessing the dynamic response to seismic ground motions or vibrations internally generated within the structural system. It is very important to note that the SPSI can alter the superstructure's performance by influencing the dynamic properties of the structure. A fundamental understanding of this complex phenomenon is essential for improvements of the structural design and construction practices, and, as a result, for reduction of loss of life and monetary damages caused by earthquakes. Extensive field testing and experimental investigations on different aspects of pile behavior have resulted in a number of empirical and approximate analytical methods for the pile-foundation design (Craig and Kurdila 2006;Han 2002;Guin and Banerjee 1998;Kaynia 1982;Dobry and Gazetas 1988;Chore et al. 2010). In addition, other studies have resulted in more rigorous schemes for pile analysis. The results of research from various studies also back up the point that dynamic impedance of foundation is affected by several geometric and material factors, such as foundation shape (i.e., disk, rectangular, strip), type of soil profile (i.e., deep uniform or multi-layer deposit, shallow stratum on rock), foundation embedment (i.e., surface foundation, embedded foundation, pile foundation), shear wave velocity and type of loading (static, dynamic) (Fan et al. 1991;Pitilakis et al. 2008).
Generally, three direct numerical approaches have so far been evolved for predicting the dynamic pile behavior. The first approach represents soil as an elastic continuum, the second represents the soil-pile system by a set of discrete masses, springs, and dashpots, and the third one uses either boundary element or finite element methods based on the actual soil properties Massumi and Tabatabaiefar 2007).
The continuum model was considered to be a generalized Winkler model. The major contribution regarding this approach was carried out by Novak (1987), who introduced an approximate method based on the analytical solution of Baranov (1967). The soil was considered to be made of a set of independent infinitesimally thin horizontal layers in this model. Therefore, it automatically considered wave propagation in the horizontal direction. A set of equations for stiffness and damping ratio based on frequencyindependent parameters was provided for a single pile under vertical vibration. Thereafter, further progress on this model was made by Novak and Nogami (1977) and Dobry and Gazetas (1988). The reduction in soil stiffness and the increase in damping associated with strong shaking are sometimes modeled crudely in these analyses by making arbitrary reductions in shear moduli and arbitrary increases in viscous damping. Moreover, several researchers developed dynamic nonlinear finite element analysis in the time domain using the full 3-dimensional wave equations to investigate the seismic soil-structure interaction (Ismail et al. 2018), but it is not feasible for engineering practice at present because of the time needed for the computations. A rigorous finiteelement analysis of a pile foundation subjected to lateral loading was done by Ülker-Kaustell et al. (2010). These numerical methods have an advantage in that they can analyze the soil-pile system as a whole. Therefore, no separate analysis is required to calculate the dynamic interaction factors for pile-group analysis (Sadeghi Hokmabadi et al. 2012).
The major drawbacks of the analytical and numerical models in prior research studies were overcome with the introduction of the cone model which applies for shallow foundations on a half space. Afterwards, Wolf has rearranged and extended the concept of the cone model to cover a complete range of dynamic excitations and physical situations. The enhanced model was based on one-dimensional wave propagation in the cone segments with reflections and refractions at layer interfaces. Due to the simplification of the physical problem, the mathematics of the cone models can be solved rigorously (Pal and Baidya 2018;Jaya and Prasad 2004;Meek and Wolf 1994).
The present study investigates the effectiveness of applying the enhanced cone model to analyze the dynamic response of pile groups in two tall and massive building projects. For this purpose, a new computer program (FCONE) was developed based on the cone frustum model. This program utilizes the equivalent dynamic modeling concept in the analysis of soil dynamics and soil-pile-structure interaction.

The Concept of Cone Model
A simplified one-dimensional cone method was first introduced by Meek and Wolf (1994) to evaluate the dynamic response of a surface foundation on a half space (Wolf 1985). Linear behavior of the site is assumed; hence, the soil is assumed to remain linearly elastic with hysteretic material damping during dynamic excitation. Figure 1 shows the loaded disk on the surface of a homogeneous half space, modeled in the onedimensional strength-of-material approach as a truncated semi-infinite bar with a vertical axis (Wolf and Deeks 2004;Wolf 1994).
An enhanced method based on the cone frustum approach has been developed to model an embedded foundation that is an order of magnitude larger than the disk, such as a pile foundation in soil media. The soil region is divided into several segments by massless rigid disks covering the volume of foundation, as shown in Fig. 2. The disks will coincide with the interfaces of the half space. In addition, further disks are selected to adequately represent the dynamic behavior. After assemblage, the dynamic stiffness matrices of all the cone frustum segments are then calculated and subsequently assembled together with the underlying half space to form the dynamic stiffness matrix of the corresponding backbone cone. Applying a unit load at the disk and solving for the displacement amplitudes at all disks lead to a column in the flexibility matrix [G(w)] of the free field condition. Thus, the force-displacement relationship is obtained as Eq.
(1), where {u(w)} is the displacement amplitudes and {P(w)} is the force amplitudes: This procedure is repeated for all disks to obtain all columns of the dynamic flexibility matrix [G(w)]. The dynamic stiffness of a single cone S 1 (w) modeling an incident wave induced by a disk load of radius r 0 on a half space is given by Eq. (2), where K s is the static translational stiffness coefficient of the cone modeling the incident wave, c is the relevant wave velocity, ω is the load circular frequency, and ρ is the soil density. The value of the aspect ratios of a particular cone is shown in Table 1 for each type of motion.
The dynamic stiffness matrix of the cone frustum segments is obtained by inverting the dynamic flexibility matrix [S f (w)] = [G −1 (w)], and, then, the forces acting on the top and bottom disks and their corresponding displacements are related as follows: It is well-known that the behavior of a pile in a compact group is quite different from that of a single pile because of pile-soil-pile interaction. These interaction effects depend mainly on pile spacing, frequency of excitation, and number of piles. These dynamic interactions are incorporated into the analysis by a factor which is defined as the ratio of fractional increase in displacement amplitude at receiver pile u r f (w) to the displacement amplitude at the loaded source pile u s (w). The dynamic normal-axial interaction factor α v (w) can be obtained from Eq. (4), where C s is shear-wave velocity, ζ g is hysteretic damping ratio, ω is frequency of excitation, r 0 is radius of the pile, and d is distance between two piles. The following equation can also be written in polar coordinates (θ, r) to obtain good estimates for any arbitrary angle, where u θ f (w) and u rad f (w) are tangential and radial displacement amplitudes at receiver pile, respectively: When analyzing soil-structure interaction, two types of damping have to be taken into account, namely the geometric attenuation (radiation attenuation) that results from the distance of the excitation source to the object and the material damping due to the friction between the individual soil particles. The latter can be taken into account by substitution of G and G (1 + 2iζ) and E by E (1 + 2iζ), with ζ as the degree of attenuation of the soil. The material damping is independent of the excitation frequency (Meek and Wolf 1994;Wolf 1985). When the distance between the piles is large (greater than about 20 diameters), the piles do not affect each other and the group stiffness and damping are simply the vector sums of the individual pile stiffness and damping. The dynamic stiffness matrix of the pile embedded in the half space can be obtained by replacing the properties of the soil between two consecutive disks with the mechanical properties of the pile section. This can be performed by adding the difference between the dynamic stiffness matrix of the pile and that of the corresponding soil cylinder [ΔS(w)], leading to: where P(w) is force amplitude and ΔK and ΔM are the difference between stiffness matrices and mass matrices of the pile and soil cylinder, respectively. Equation (6) can be written as: P z P x p y m z m x m y 2 6 6 6 6 6 6 4 3 7 7 7 7 7 7 5 where v, h, r, and t denote vertical, horizontal, rocking, and torsional motions, respectively. It can be noted that coupling occurs between horizontal and rocking motions. Each coefficient of the impedance matrix is a complex function of the excitation frequency which can be expressed in the form of real and imaginary components. Furthermore, the dynamic stiffness coefficient obtained from the analysis is then decomposed into two parts: In Eq. (8), K st is the static stiffness, i the imaginary unit, and k ij and c ij , respectively called stiffness and damping coefficients, are dimensionless coefficients depending on the Poisson's coefficient of the soil and the frequency parameter.
where C s ¼ ffiffiffiffiffiffiffiffi ffi G=ρ p is the shear wave velocity in the soil. Finally, the total displacement of receiver pile (i) due to lateral load of source pile ( j) can be written as: or In order to achieve a sufficient accuracy, the number of disks n has to be selected so that the length L of a pile element does not exceed about one tenth of the shortest wavelength propagating from the pile. Because the shear wave velocity is less than the dilatational wave velocity in the soil, the wavelength of shear waves is lesser than that of the longitudinal ones and the number of disks can be chosen as (Wolf and Deeks 2004;Wolf 1994): Transportation Infrastructure Geotechnology (2020) 1:13-45

Case Study 1: (an Emergency Hospital in South of Switzerland)
The building is of a general rectangular shape with the overall dimensions of 73 m in the longitudinal direction, 25 m in the width direction, and a height of 35.6 m. The building is a new reinforced concrete construction with concrete columns. The primary lateral load-carrying system consists of three reinforced cores and an integrated reinforced concrete exterior frame system. For the earthquake resistance evaluation, a special micro-zoning map and a response spectrum developed by Resonance Ltd. for Brig were used. The building is classified as a category III construction with an importance factor of 1.4. This is the highest class of building categories in Switzerland. The building was modeled in three dimensions (3D) and calculated for non-ductile behavior according to the SIA Standard 262 (2013 edition). The summary of the building characteristics and modeling assumptions is given in the succeeding texts.

Longitudinal Direction
The earthquake loads are carried by the three reinforced concrete cores and two shear walls together with the reinforced concrete framework. As the geometry of all floors is similar, continuous transfer of earthquake loads is guaranteed from the roof to the foundation. Additionally, due to the layout and stiffness differences, the mass center and the stiffness center of the building do not match resulting in torsional stresses due to lateral earthquake effects.

Lateral (Short) Direction
Three reinforced concrete core and two reinforced concrete frames transfer the horizontal equivalent loads corresponding to the earthquake effects from the higher levels through the lobby level to the basements of the building. There is a concern of soft story in the building. As a result, high shear stress concentrations occur in the building due to this irregularity. Due to the stiffness differences, the resistance of the building in the short direction is lower than that in the longitudinal direction.

Numerical Investigation
The earthquake loads are determined according to the Swiss Code SIA 260-262; Merkblatt SIA 2018. The following material characteristics have been used: concrete: Ec = 20,000 Mpa (elastic modulus); f ′c = 19.2 Mpa (compressive strength).
For the seismic structural analysis, the following four computer programs were used: -ETABS, which is a well-known structural analysis and design software from CSI, Berkeley, USA. -OpenSees or Open System for earthquake engineering Simulation. This is an earthquake simulation software based on the displacement-based structural analysis considering material non-linearities. The program is very suitable for assessing the behavior of buildings in seismic situations. The purpose of modeling with OpenSees is to perform a non-linear static pushover (NSP) analysis of the reinforced concrete system with concrete cores, shear walls, and the heavy concrete columns. The modeling is based on ASCE 41-13 (2014). Material non-linearities and hysteretic behavior are considered for the elements under normal forces, bending, and shear effects. For the calculations, both force-based and displacement-based methods have been used. For the displacement-based method, ASCE 41 procedures have been implemented in the calculation of target displacements. -FCONE program was a new FORTRAN program developed at Smteam based on the enhanced cone model concept for soil modeling. This program utilizes the equivalent dynamic modeling concept in the analysis of soil dynamics and pilefoundation interaction. -Abaqus/Cae was used to investigate the structure response according to seismic soilpile-structure interaction effect. The finite element simulation phase of the present study along with a description of boundary conditions, the meshing and analysis techniques, and the materials' constitutive laws is also presented in the subsequent texts.

Pile Group Effect
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 non-linearity. Numerical analysis methods of soil-pile-structure interaction (SPSI) problems can be classified into the following two approaches: direct method and substructure method. In the direct method for analysis of SPSI, it is necessary to model underlying soil, foundation, and structure with best precision and accuracy since they are supposed to be analyzed together. In the substructure method, soil-structure system is divided into the following two parts: the first part is the structure and the second part is the soil that has common boundary with the foundation and in which each part can be modeled in detail (Wolf 1985). In order to consider the influence of the piles on the dynamic response of the structure, the substructure method was used in this study. This consists of three consecutive steps: (1) Calculation of foundation input motion in the absence of superstructure inertia (kinematic interaction) includes translational displacements and rotational components.
(2) Calculation of the dynamic impedances (spring and damping coefficients) for horizontal rotation and the combination of base foundation and pile system.
(3) Calculation of the seismic response of the building supported on the springs and dashpots of step 2 and subjected at its base to the foundation input motion of step 1.
The method of Fan et al. (1991), based on the BDWF model (Beam on Dynamic Winkler Foundation), is used to determine the foundation input motion. In this method, the displacement of the pile head (U11) as a function of the free-field motion (Uff) is calculated as follows: Here, Γ is a function of spring and damping constants, frequency, and mechanical characteristics of the piles, such as Ip, Ap, Mp, and Ep.
A detailed study of gazetas shows that the kinematic factor Iu = U 11 / U ff for small frequencies, like the first frequency of the hospital Brig (0.7 Hz), is equal to 1. This means no increase in the amplitude of the free-field earthquake motions. Therefore, in the case of the Brig hospital, the kinematic interaction effect was neglected.
The arrangement of pile-group in rectangular foundation is shown in Fig. 3. In this case, both rigid and flexible foundations are considered. In the fixedbase model, the structure is being analyzed without soil and foundation with fixed boundary condition to resist all displacements or rotations at the bottom.
In the cone model, the structure is built on a rigid foundation with springs and dashpots in the bottom in place of the soil. The response spectrum method was used for determining the displacements and accelerations.
To determine the vertical axial load in the upper part of the piles, a finite element model of the building was created using the ETABS program. In the modeling of the core, 3D shell elements were used. The elastic modulus E of the concrete has been reduced to 60% of the design value to include the effects of cracking in the concrete. The (q) factor (response factor) is set to 2, and the importance factor is assumed as 1.4. In 2015, a new earthquake hazard map of Switzerland ( Fig. 4) was issued with a corresponding design spectrum, which was used for the present study. The modal analysis was performed to determine natural modes of frequencies for the structural model in two different conditions, namely fixed-base behavior and behavior including SPSI (containing embedded piles with pile-group effect).
The results of ETABS analysis are given in Table 2, and the first vibration mode shape of the building with dynamic pile group effect obtained from the ETABS is shown in Fig. 5.

Target Displacements
The calculated target displacement in long and short directions according to the Euronorm 8 and ASCE 41-13, 2014, is found to be 12 and 15 cm, respectively. These values represented the demand on the structure to have a sufficient ductility level under earthquake loads before a collapse of the building can occur.

Non-Linear Pushover Analysis
A pushover analysis involves the application of static lateral loads distributed over the height of the structure and calculation of the resulting displacements in a # Pile Type Diameter (Cm) Length (m) 1 100 22 2 100 28 3 120 30 Fig. 3 The arrangement of pile group in rectangular foundation model of the SPSI system. A pushover analysis of a structure with a flexible base is schematically illustrated in Fig. 6. A non-linear pushover analysis was reported for the reinforced concrete system model of the building with the parameters defined in ASCE 41-13. Behavior of each wall was determined based on the wall thickness, the wall length, and the material properties as provided by the ETABS model. At each point on the pushover curve, the deformations of all components in the structural system are related to the reference displacement. The ETABS model has been read into OpenSees for conducting the pushover analysis. The pushover analysis resulted in displacement capacities of 20 cm in the transverse direction and 22 cm in the longitudinal direction, as shown in Fig. 7. This is compared with the target displacements, found from ASCE 41-13, of 12 cm in the longitudinal direction and 15 cm in the transverse direction.

Dynamic Stiffness of Pile Groups
Dynamic stiffnesses (spring and damping coefficients) for a pile group are strongly frequency-dependent. In the FCONE program reported in the article, the impedance matrix for different blocks (pile groups) was calculated as a function of frequency, using the cone model concepts according to John P. Wolf (Meek and Wolf 1994). For the determination of stiffness, a new method of solving dynamic interaction problems was developed wherein the seismic behavior of pile-soilpile-structure interaction (SPSI) system is investigated as a whole considering both pile group effect and soil constitutive model. As mentioned previously, modifications of the flexibilities have to be performed for modeling an embedded Fig. 4 The response spectrum for the Brig ground motion and the seismity map of Switzerland Table 2 The results of two analyses obtained from the ETABS mode

Period of vibration (s)
The roof displacement (cm) Fixed-base foundation (no piles) 0.7 6.5 Flexible foundation (including piles) 1.2 11.5 foundation since the double cone models (especially for the vertical degree of freedom) are too flexible at large distances from the loaded disk. The pile-soilpile interaction effect in the pile-group has been evaluated through the FCONE program by using the matrix of interaction factors (α υ ), which is frequencydependent (Stacul and Squeglia 2018;El Naggar and El Naggar 2007). Results of the program for one 3 × 3 block under a core (i.e., spring and damper constants as a function of frequencies) are summarized in the following section.
The spectral values have been reduced to take account of the radiation and material damping of the soil. It turns out that in the case of the flexible-base structure, in addition to the added rocking component to the horizontal motion of the structure, a part of the structure vibration energy will transmit to the soil layer and can be dissipated. This results in a reduction of the natural frequency in comparison to the fixed-base model and an increase of the damping (mainly due to radiation damping resulting from the wave propagation and hysteresis damping of the soil materials).
Detailed mechanical characteristics of the soil have been explained in Table 3, where E is elastic modulus, Gs is shear modolus, ρ is unit weight, Cs is shear wave velocity, υ is Poisson's ratio, and ξ is material damping ratio. The soil region was assumed to be homogeneous with uniform shear modulus and density during Fig. 5 The first vibration mode shape of the Brig structure (period: 1.195) Fig. 6 Schematic illustration of a pushover analysis and development of a pushover curve for a structure with a flexible base (Stewart et al. 2012) analysis (Stewart et al. 1999a;Stewart et al. 1999b). The horizontal, vertical, and rotational dynamic stiffness and damping ratios of a soil layer captured by using FCONE program are shown as a function of dimensionless frequency in Fig. 8.
The dynamic stiffness coefficients of an embedded foundation, such as a pile foundation, in a homogeneous half space are sensitive to variations in the soil material damping ratio, which is in turn a function of strain amplitude. In the lower frequency range (a 0 < 1), considerable deviations exist for certain damping coefficient (Wolf and Deeks 2004;Meek and Wolf 1991). Figure 8b presents the dynamic behavior below and above the so-called cut-off frequency, where an abrupt change in response occurs. As can be seen in Fig. 8, the variation in damping coefficient at a lower frequency is much greater than that at a higher frequency.

The Finite-Element Simulation
The finite element code ABAQUS/Cae model as an accurate direct way was used to investigate the seismic behavior of the structure under the SPSI effect. For this purpose, the infinite soil medium is reduced to a finite region using infinite element boundaries. In this case, an elastic perfectly plastic model was adopted for modeling

Displacement Capacity
Target Displacement Fig. 7 The pushover curve obtained by OpenSees analysis the concrete piles with a modulus of elasticity of 30,000 (N/mm 2 ). In the 3D finite element model developed using Abaqus/Cae, both the pile and the soil were modeled using eight-nodded solid continuum elements (C3D8R) to account for the continuum nature of the soil. This element has a quadratic shape function which is well-suited to model the medium with bending-dominated deformation. A too a) Horizontal Stiffness for a 3×3 Pile group b) Vertical Stiffness for a 3×3 Pile group Fig. 8 The vertical, horizontal, and rotational dynamic stiffness for a 3 × 3 pile group coarse mesh may deviate from the expected response, whereas a too fine mesh requires very long computational time. Hence, a relatively fine mesh was adopted for the pile and a coarser mesh was adopted for the soil. Overall, 11,089 elements for the piles and 10,205 elements for 10-story SPSI model were used. The soil was modeled using 98,400 elements, and the soil-pile interaction was modeled by defining the tangential and normal behavior of the contact surfaces in the FE model (Ghandil and Behnamfar 2015;Bao et al. 2012).
A strain-hardening model using Mohr-Coulomb failure criterion was adopted for the soil which uses the five significant soil parameters, such as shear modulus (G), Poisson's ratio (ν′), shearing resistance angle (φ′), effective cohesion (c′), and angle of dilatancy (ψ), as shown in Table 4. An angle of internal friction of 35°was used in the definition of the Mohr-Coulomb failure criterion. The geometrical details adopted for the SPSI analysis are shown in Fig. 9.
Results of 3D numerical prediction of first mode of vibration in eigenvalue analysis for both cases (fixed-base and SPSI interaction) and maximum lateral displacements in short direction are shown in Figs. 10 and 11, respectively. Table 5 provides a good comparison between the results of SPSI analysis obtained by the cone method and FEM solution, which indicate that the results of the cone method are in close agreement with the FEM results. As the structural stiffness in comparison with soil or building height increases, the period rate increases as well. For the first mode of vibration in modal analysis, a good tendency with negligible difference was observed. Therefore, the cone model provided reasonable accuracy of modal analysis for the high-rise building responses because the deviation was limited to 3%. From the comparison between the maximum lateral displacements, it can be observed that the foundation displacement results obtained from the cone method give good agreement with the FEM analysis results with about 6% difference, and the most difference was observed in the roof displacement with around 10%. The non-linearity of the soil medium results in increased lateral displacements of stories; however, as can be seen by comparing the results, it does not have a significant impact on the lateral displacements in this case. It can be also concluded that the maximum lateral displacement of structures increases with the increasing number of stories, and the upper stories' displacements are more affected with SPSI than the lower stories.

Dynamic Analysis of an Emergency Hospital in South of Switzerland
Specification of the free-field ground motion is one of the most important factors in SPSI analysis. A range of intensities, durations, and frequency contents was covered by selecting the ground motions to enable a comprehensive evaluation of the SPSI effect. In this research, the used earthquake excitation was the 1940 El Centro earthquake which has peak ground acceleration (PGA) of 0.349 g, commonly used in structural dynamics. The ground motion records aligned in the two horizontal directions (shown in Fig. 12) were selected from the strong motion database of the Pacific Earthquake Engineering Research Center (PEER 2012).
The time history analysis was about 50 s duration consisting of 2500 steps under an actual earthquake accelerogram and a free vibration segment of 0.02 s  stories of fixed-base and flexible models were measured by using the ETABS program. As mentioned, the significance of the displacement mechanisms was expected to depend on the excitation, soil, and structure properties. In particular, lateral displacements and story drifts of the structures resting on rather soft soils may be significantly amplified by SPSI (Jennings and Kuroiwa 1968). The time histories of roof displacements under earthquake excitation with SPSI effect along X and Y directions are demonstrated in Fig. 13. In addition, the maximum lateral displacements and story drifts in both cases are shown versus the building floors in Fig. 14 (parts a and b) and (parts c and d), respectively.
From the comparison, it can be observed that the top deflection of the 9th-floor building (roof) with SSI effect increased by 1.77 times (and 2.12 times) than the case without taking into consideration the SSI effect, along the X direction (and the Y direction), respectively.
It should be noted that soil-structure interaction analysis typically involves dealing with several uncertainties. The results of such dynamic analysis may require determining the limitation of material properties and models. In real engineering projects, the Fig. 11 The maximum lateral displacements in short direction limitation ranges are determined in such a way that they stay on the conservative side.
The uncertainties have remarkable impact on the structure response, and, hence, it is essential to be cautious in this matter.

Case Study 2: Dynamic Analysis of the Green Building at MIT Campus
The Green building is Building 54 at the MIT campus in Cambridge, Massachusetts, USA. This building is the tallest structure in Cambridge, designed and It has 20 stories with an overall height of 83 m plus a one-story basement below ground level. The height of the first story is 10 m, and that of the others is 3.5 m. The ground surface is approximately 6.1 m above sea level, and depth to bedrock is approximately 30.55-40.0 m below sea level. Building 54 includes administrative sections and classrooms for Natural Sciences department. As shown in Fig. 15, the building orientation in the plan has a deviation of 25°to the north. Two 25-cm-thick shear walls extend from the basement-foundation level to the roof level at the two narrow ends.
The underground floor has a link to the MIT colleges' connecting tunnel. The structure was made of cast-in-place reinforced concrete, and the thickness of the concrete slab at each story is 101 mm. In Fig. 16, the position of the shear walls, the dimensions, and the sizes are specified. The building is placed on a concrete mat foundation with concrete pile caps that have thickness of 350 mm and loading capacity of 50 tons. The foundation dimension is 34 × 16.5 m. Despite the existence of doors, windows, and entrances, from the shear walls up to the roof level, the building still could be considered as a symmetric structure. Figure 17 shows a typical depth versus shear wave velocity log from a site approximately 50 m (150 ft) from the building (Jen 2011;Haley and Aldrich 1983).

Site Condition
Building 54 is located near the Charles River. The structure is placed on the alluvial deposits of the river. Geotechnical studies of the site condition have been thoroughly performed, and the soil properties, especially shear wave velocity of the soil, were determined.
The site was 6.1 m above mean sea level, and, according to geotechnical studies, the shear wave velocity varied between 100 and 600 m/s. The depth of bedrock was estimated to be around 30.55 to 34 m from the ground surface. Based on a shear wave velocity boring log from a site approximately 50 m from the building (Jen 2011;Haley Fig. 18 The soil profile and the wave propagation in the soil layers and Aldrich 1983), the fundamental site frequency is estimated to be ∼ 1.5 Hz. Figure 18 displays the soil profile and the wave propagation in the soil layers. Later, Trocha (2013) studied the dynamic behavior of this building, again using ambient data. The natural frequencies were identified-with some differences from those identified by Çelebi et al. (2014) as 0.75 Hz (1st NS), 2.85 Hz (2nd NS), 0.68 Hz (1st EW), 2.45 Hz (2nd EW), and 1.45 Hz (1st Tor). The fundamental frequency of the building (0.68 Hz) is less than the frequency of the site; hence, no radiation damping is  considered. In addition to geotechnical data, the site seismicity data were also available. In 2012, the magnitude 4.0 earthquake occurred near Hollis Center, Maine, where data recorded with accelerometers were available. In the present study, three ground motions including recorded ground motion in Boston resulting from (October 16, 2012) Hollis Center Earthquake, Imperial Valley (El Centro in 1940) earthquake and employed ground motion in the rehabilitation of Hotel Azadi, Tehran, Iran are considered as the input motions for the dynamic analysis and structure evaluation.

Modeling and Dynamic Analysis
The modeling of Building 54 has been performed in detail, meaning that the structure was simulated along with the surrounding soil. Computation of maximum lateral displacement for all stories of the building corresponding to the soft soil during the three mentioned earthquakes was done based on both conditions with and without SSI. Although the structure data were clear and definite, there was considerable uncertainty in soil properties which caused complex challenges in the analysis. The "springs and dashpots" method proposed by the NEHRP Seismic Provisions (Stewart et al. 2012) allows evaluating the mechanical and dynamical characteristics of soil that have been impacted by a combination of three mechanisms, including footing conditions,  The MIT building has a mat foundation with 34.0 m in length and 20.5 m in width, which is supported by pile groups. Pile groups are fixed with pile caps located at both longitudinal directions of the building (i.e., EW side). The length of the piles is 15 m, and their diameter is 0.3556 m, and each group contains 30 piles. In order to capture the soil-structure interaction effects, it would be needed to consider combined response of mat foundation and pile groups.
Kinematic interaction due to existence of stiff elements causes recorded ground motions in the free field to be significantly less than what is imposed to the foundation. The main reasons associated with this phenomenon are the base-slab averaging (BSA), Period Ҫelebi et al.    embedment effects, and kinematic pile response. In this section, foundation input motion (FIM) is calculated based on the proposed method in the NIST report (Stewart et al. 2012). In this regard, the following steps are followed: & Fourier transform of free-field ground motion is calculated. & Amplitude of Fourier transform is modified by consideration of base-slab averaging and embedment effects. & Inverse of Fourier transform is performed to obtain foundation input motion. Figure 19 depicts the acceleration records of the Hollis Center earthquake in the following three situations: free field, modified by base-slab averaging (BSA) consideration, and modified with both BSA embedment and embedment effect. As shown in Fig. 20, the Fourier transform of these cases is multiplied to consider kinematic interaction effects.
Dynamic soil stiffness parameters, i.e., stiffness and damping parts as well as effective foundation input for vertical, horizontal, and rocking degrees of freedom, are calculated using the FCONE program. Outcomes are depicted against the dimensionless frequency parameter, a 0 = w2r 0 /c s (Figs. 21, 22, 23).
The mentioned characterization and structure evaluation method is normally used for seismic analysis of infrastructures as it was also applied to Building 54. The first three identified mode shapes obtained from the ETABS analysis are shown in Fig. 24, where it can be seen that the application of SSI could cause approximately 10% increase in the captured periods of the building. A comparison was also made between the results of modal analysis by Celebi et al. (2014) and these 3D numerical predictions captured by using ETABS and FCONE programs. Good agreements have been achieved in the comparison of modal analysis results.
In the following section, time-history analyses are conducted under previously considered ground motion records which were modified by consideration of kinematic interaction. In general, SSI tended to amplify the lateral deflections of the superstructures. Experienced roof displacements under excitation of considered ground motions with and without SSI effects are summarized in Table 6. Besides, the results of time-history analysis for all cases are displayed in more detail, as shown in Fig. 25. It can be seen that the displacements increased by consideration of the SSI effects.
To determine the sensitivity in the structure's response from an analytical perspective, we computed both minimum and maximum horizontal accelerations under excitation of considered ground motions with and without SSI effects, which is presented in Table 7.
Moreover, the horizontal roof acceleration response history of the building during the considered earthquakes is plotted in Fig. 26, for fixed-base conditions. It can be seen that the experienced roof accelerations are decreased by consideration of the SSI effect. Figure 27 shows a schematic of the building depicting important dimensions as well as locations and orientations of the 36 accelerometers deployed throughout the building. Each accelerometer channel is connected via a cable to a 36-channel central recording system. The accelerometers are distributed as three per level to

Conclusion
The present study investigates the effectiveness of applying the enhanced cone model to analyze the dynamic response of pile groups in two tall and massive building projects. For this purpose, we have developed a new computer program in which the seismic behavior of pile-soil-pile-structure interaction on a large building is investigated as a whole considering pile-group effect and soil constitutive model according to the cone theory. This program utilizes the equivalent dynamic modeling concept in the analysis of soil dynamics and pile-foundation interaction. The main conclusions drawn from the results are given as follows: & The enhanced cone model allows the non-linearity of the near-field soil (plastic zone) to be represented by an equivalent linear manner, thereby reducing the degrees of freedom compared to the FEM method and consequently reducing computational effort. The results presented in this article are obtained from substructure approach in the frequency domain which is formulated on the basis of the cone method derived from the fundamental solution for a homogeneous, isotropic, and linear-elastic continuum. The presented model is also capable of incorporating the angular wave incidence from the far-field. & By comparing the obtained results, it can be concluded that considering the effects of the SPSI can alter the dynamic characteristics of the structural system. It is shown that the effects of soil-structure interaction are demonstrated as increasing the period of vibrations and the lateral deflections as well as decreasing the experienced roof accelerations. Consequently, the in-story drifts of the structures were amplified in comparison with the fixed-base model. Therefore, ignoring the SPSI effects in such cases could affect the predicted damage level of structural and non-structural elements during an earthquake, resulting in considerable deviations. & Linear (frequency domain) SSI analysis can be employed for prediction of dynamic response of structures in low-to-moderate-intensity shaking earthquakes, such as these real cases that were located in a moderate-level seismic zone. Besides, the consequences might differ for extreme shaking that results in highly non-linear soil where the non-linear effects, including gapping, sliding, and uplift, are greatest in the immediate vicinity of the soil-structure boundary, and these cannot be captured using equivalent linear techniques. & The capability of the cone method to predict the dynamic SPSI effects gives it an advantage over other time-consuming rigorous methods. This approach leads to physical insight with conceptual clarity, which is simple to use and solve, as the mathematical solution is simplified, and provides sufficient generality (layered site, embedment, all frequencies) and acceptable engineering accuracy. The accuracy of any analysis is limited anyway, because of the many uncertainties, some of which can never be eliminated (for instance, the definition of the dynamic loads and the values of the dynamic soil properties in the analysis of a structure). The analysis with cones fits the size and economics of practical engineering projects. Thus, the results presented in this article can contribute to relevant engineers in designing safe and economical structures based on foundation vibration and dynamic SPSI analysis in practical projects whenever possible.