# Soil-Structure Interaction in Transversely Isotropic Layered Media Subjected to Incident Plane SH Waves

Authors: Zhenning Ba; Xi Gao
Journal: Shock and Vibration
 (2017)

---

## Abstract

The dynamic soil-structure interaction (SSI) for incident plane SH waves is analyzed for a two-dimensional (2D) model of a shear wall on a rigid foundation by using the indirect boundary element method (IBEM). The rigid foundation utilized in this study is embedded in transversely isotropic (TI) soil layers over bedrock. The accuracy of the IBEM method is verified and analyzed by setting a semicylindrical, rigid foundation-shear wall structure system in the single TI soil layer and multiple TI soil layers over bedrock. This study shows that the TI characteristics of the site have a significant impact on the effective input motion and the superstructure response. In a single soil layer, the increase in the shear modulus ratio in the vertical and horizontal directions has a certain degree of amplified action on the effective input motion and the superstructure response. Simultaneously, the corresponding peak frequency of the response increases. In multiple soil layers, the changes in the effective input motion and the superstructure response are also affected by the TI characteristics of the soil layers, and the impact of this effect is related to the sequence of the layers.


---

## Body

## 1. Introduction

Dynamic soil-structure interactions (SSI) under the effect of seismic waves are important problems in the fields of soil dynamics, seismology, and earthquake engineering, because of their particular significance to practical engineering applications. Methods for solving soil-structure interaction problems can be divided into analytical and numerical methods. Early analytical methods utilized wave function expansion to obtain the out-of-plane response of a two-dimensional (2D), semicircular, rigid foundation (with a shear wall representing the superstructure) when excited by plane SH waves [1, 2]. Subsequently, many scholars adopted analytical methods to study dynamic soil-structure interactions, including the dynamic response of the interaction of a structure system with different foundations. Semielliptical and hemispherical foundations are some examples of differing foundations [3–6]. In addition, the interaction of a soil-structure system under different soil conditions (dry, saturated, etc.) was studied [7–9].The analytical method is limited to the dynamic response of foundations with a single shape in a homogeneous half-space. However, in recent years, the development of computer technology has promoted the application of numerical methods, such as the finite element method and the boundary element method. The key to study the dynamic interaction between soil and structure using the finite element method is to apply the appropriate artificial transmitting boundary. Lysmer [10] was the first to utilize the finite element method in this manner and calculated the dynamic response of a three-dimensional foundation embedded in a homogeneous half-space by setting a viscous absorbing boundary condition. Tassoulas and Kausel [11] calculated all of the stiffness coefficients of foundations embedded in a layered half-space using the semidiscrete method and setting a rigid transmitting boundary. They also discussed the influence of the embedded depth of the foundation in detail. Subsequently, many other scholars proposed other boundary conditions, including the rigid transmitting boundary [12], the iterative transmitting boundary [13], the Smith boundary [14], the paraxial boundary [15, 16], and the Higdon boundary [17]. These researchers conducted in-depth studies on the dynamic interactions between soil and structure.The boundary element method (BEM) does not require artificial boundaries because it uses Green’s functions as a basic solution, which automatically satisfies the radiation conditions at infinity. Green’s functions have also been used to calculate the impedance functions of embedded foundations in [18]. Karabalis and Beskos [19, 20] studied the stiffness coefficients of three-dimensional and two-dimensional embedded foundations and their dynamic response under seismic action using the time integral method. de Barros and Luco [21] calculated the dynamic response of a two-dimensional semicircular foundation embedded in a layered half-space using Green’s function. Fu et al. [22–25] studied the dynamic response of the interaction of a two-dimensional, semicylindrical, rigid foundation-shear wall system in single and multiple soil layers using the boundary element method combined with Green’s function distributed on a slope.In order to simplify the model, most scholars assumed the subject site to be a homogenous or isotropic half-space, which was different from the actual site. Studies have shown that the soil in an actual site often is anisotropic in different degrees, and the anisotropy is mainly reflected in the difference in elastic properties in the horizontal and vertical directions. Therefore, it is more realistic to consider the ground soil as TI soil [26]. Compared to isotropic ground soil, there is much less research available on the dynamic interaction of a soil-structure in TI soil. In existing studies, Wang and Rajapakse [27] calculated the dynamic response of a two-dimensional foundation under low and medium frequency using Green’s function [28, 29] for anisotropic soil; Kirkner [30] studied the dynamic interaction of TI half-space and a circular foundation for the first time; and Wu et al. [31] used an analytical method to investigate the vibration problem of a rigid round plate on TI-saturated half-space under the action of harmonic torsional load. Although these studies have adopted the site model of a TI soil layer, they did not have access to a systematical analysis or research on the influence of the characteristic parameters of TI ground soil on the dynamic response of soil-structure.In this study, the layered characteristic of natural soil has been considered, Green’s function of distributed loads on inclined lines with the indirect boundary element method (IBEM) is adopted as an elementary solution, and the solution of the effective input motion and the superstructure response in the dynamic interaction of soil-structure in TI ground soil has been given. By comparing with the effective input motion and the superstructure response in isotropic ground soil, the method has been verified; the effective input motion and the superstructure response are calculated by setting the single TI soil layer and the multiple TI soil layers as examples. In addition, the influence of the TI characteristics of the ground soil on the effective input motion and the superstructure response has been analyzed and discussed.


## 2. Methodology

### 2.1. Model

The model shown in Figure1 is embedded in TI-layered ground soil and consists of a rigid semicylinder foundation that supports an elastic shear wall. The TI ground is formed using N layers of horizontal TI soil and TI bedrock (half-space) underneath. In TI soil, the axis of symmetry of the medium is assumed to be normal to the horizontal surface. The TI soil layer is determined by the shear modulus in plane normal to axis of symmetry GhiL, horizontal shear modulus, the shear modulus in planes normal to plane of transverse isotropy GviL, vertical shear modulus [32, 33], and the damping ratio ζiL(i=1~N). The TI bedrock (half-space) is determined by the shear modulus in the horizontal direction GhR; the shear modulus in the vertical direction GvR; and the damping ratio ζR. The superscripts L and R represent the soil layer and the bedrock (half-space). The cross section of the foundation is semicircular, but it extends along the y-axis with radius of a and a mass of M0. The shear wall has the same width as the foundation with a shear wave velocity of βb, a height of H, a mass of Mb per unit length, and a damping ratio ζb. The shear wall has a rigid connection (no slip) between the rigid foundation and the TI-layered ground soil, and the interface between the foundation and the ground soil isS. Incident SH wave inputs from angle ψ in the bedrock surface have a circular frequency of ω. In the case of incoming SH waves, the foundation response (i.e., the effective input motion) and the superstructure response are sought after.Figure 1
Foundation-shear wall model.The main steps for determining the foundation response (i.e., the effective input motion) and the superstructure response include (1) the direct stiffness method which was introduced in literature [34] is adopted to determine the response of the free field in TI ground soil; (2) the indirect boundary element method (IBEM) from literature [35] is used to determine the dynamic stiffness coefficient outside of the embedded foundation plane; (3) the effective input motion is solved for, in other words, the sum of the following two parts of the response which is determined: (a) the response when the mass M0 of foundation and the mass Mb of superstructure are not taken into account; (b) the additional displacement response when the mass inertia forces F0 and Fb of the foundation and superstructure are taken into account; (4) the method in [2] is used to determine the superstructure response; (5) the relative effective input motion and superstructure response with respect to the free field are determined, namely, the specific value of the displacement amplitude of that point before and after the existence of superstructure.


### 2.2. The Response of the Free Field in Layered TI Ground Soil

The response of the free field refers to the site response when there is no foundation or superstructure. Based on the accurate out-of-plane dynamic stiffness matrix for TI-layered medium given in the literature [36], the direct stiffness method introduced in literature [34] is adopted to determine the response of the free field in TI ground soil. For the stiffness matrix of the TI soil layer, we assume that there are ascending and descending SH waves in the layer, and then the displacement at any point in the layer can be expressed as (1)vx,z=ASHexp⁡iktz+BSHexp⁡-iktzexp⁡-ikx=vzexp⁡-ikx,where v(z) can be considered as the amplitude transmitting along the x-direction(2)vz=ASHexp⁡iktz+BSHexp⁡-iktz,where ASH and BSH are the amplitudes of the ascending and descending SH waves in the soil layer, respectively, k=ωcos⁡ψSHL/csL∗ is the wave number along the x-direction, and t=-i1-1/cos2ψSHL and ψSHL are the included angles of the ascending or descending SH waves in the soil layer in the horizontal direction(3)csL∗=cshL∗2cos2ψSHL+csvL∗2sin2⁡ψSHL,which is the complex wave velocity of SH waves transmitting in the ψSHL direction. cshL∗=GhL∗/ρ and csvL∗=GvL∗/ρ are the complex wave velocities of SH waves in a soil layer transmitting in the horizontal and vertical directions, GhL∗=GhL1+2iζL and GvL∗=GvL1+2iζL are the complex shear moduli in the horizontal and vertical directions, GhL and GvL are the horizontal and vertical shear moduli when there is no damping in the soil layer, and ζL is the damping ratio.By substitutingz=0 and z=d (d is the thickness of soil layer) into (2), the displacement amplitudes v1 and v2 of the top and bottom of the soil layer can be obtained as(4)v1v2=11exp⁡iktdexp⁡-iktdASHBSH.Through the relational expression of stress and displacement, the shear stressesτyx and τyz can be obtained as(5)τyx=GhL∗∂v∂x=iktGhL∗ASHexp⁡iktz-BSHexp⁡-iktzexp⁡-ikx,(6)τyz=GvL∗∂v∂z=iktGvL∗ASHexp⁡iktz-BSHexp⁡-iktzexp⁡-ikx.By substitutingz=0 and z=d into (6) and introducing the amplitudes of the external load, Q1=τyz(0) and Q2=τyz(d), the amplitudes of the external load on the top and at the bottom of the soil layer can be obtained as(7)Q1Q2=iktGvL∗-11exp⁡iktd-exp⁡-iktdASHBSH.Eliminating coefficientsASH and BSH in (4) and (7), the layered stiffness matrix SSHL of the TI soil layer is (8)Q1Q2=SSHLv1v2=ωsin⁡ψSHLGvL∗csL∗sinωsin⁡ψSHL/csL∗cosωsin⁡ψSHLcsL∗d-1-1cosωsin⁡ψSHLcsL∗dv1v2.Under an infinite radiation condition, the outgoing wave with an amplitude ofBSH will be generated if load is arranged on the surface of TI half-space. The subscript 0 represents the free surface of the TI half-space in (4) and (7). Setting ASH=0, v0=v1, and Q0=Q1=-τyz0, the stiffness coefficient SSHR of the TI half-space will be obtained after the elimination of BSH(9)Q0=SSHRv0=iGvR∗ωsin⁡ψSHRcsR∗v0.After determining all of the stiffness matrices in the TI soil layer and the TI half-space, the accurate dynamic stiffness matrixSSH of the TI-layered site can be obtained, by combining SSHL and SSHR(10)SSH=k11L1k12L1k21L1k22L1+k11L2k12L2k21L2k22L2+k11L3⋱k22LN-1+k11LNk12LNk21LNk22LN+SSHRN+1×N+1,where kijLn (n=1~N,ij=1~2) denotes the ijth submatrix of the nth TI layer’s dynamic stiffness matrix SSHL. Therefore, the discrete dynamic equilibrium equation in the TI-layered site can be expressed as(11)Q1,Q2,…,QnT=SSHv1,v2,…,vnT,where Q1,Q2,…,QnT is the vector of the amplitude of an external load acting on the interface of the soil layer and v1,v2,…,vnT is the vector of the displacement amplitude on the interface of the soil layer.After obtaining the overall dynamic stiffness matrix of TI-layered site, the control point for SH wave transmitting from the bedrock surface can be selected on the outcrop of the bedrock, so that the other elements are zero in the vectorQ1,Q2,…,QnT of amplitude of external load, only(12)Qn=SSHRv0=iGvR∗ωsin⁡ψSHRcsR∗v0,where GvR∗ is the complex shear modulus of the bedrock half-space in the vertical direction, ψSHR is the incident angle of SH waves in the bedrock half-space, and csR∗ is the complex wave velocity of SH waves of the bedrock transmitting in the incident direction. In (3), v0 represents the movement of the bedrock outcrop. Substituting (12) into (11), the displacement in the interface of the TI soil layer can be obtained, and the amplitude coefficient of the ascending and descending wave in any soil layer can be obtained using (4). Finally, the amplitude of displacement and stress in arbitrary point of the TI-layered half-space can be obtained using (1), (5), and (6), namely, the free field response of SH waves.


### 2.3. Green’s Function of the Uniformly Distributed Inclined Load Acting in the TI Medium

In this section, the Green’s function of the uniformly distributed inclined load acting in the TI medium is derived, including both the displacement and stress, which serve the elementary solution of IBEM to obtain the stiffness coefficient of the foundation. The specific deviation of Green’s function is described as follows.Firstly, the load uniformly distributed on inclined lines in the spatial domain is expanded to the wave number domain. Secondly, two fictitious interfaces are introduced at the upper and lower boundaries of the loaded layer. In order to derive the dynamic response of the loaded layer, the two fictitious interfaces have to be fixed to obtain the counterforces of the two ends. Actually, the dynamic response of the loaded layer includes the particular solutions induced by the load applied, homogeneous solutions caused by the reaction loads and the solution derived from the reverse counterforces at the two ends. However, for those unloaded TI layers, their solutions only contain those homogeneous ones induced by the reverse counterforces at the two ends of the loaded layer. Finally, the solutions for every TI layer in the spatial domain are obtained by superposing a Fourier inversion on those solutions in the wave number domain, that is, Green’s function of uniformly distributed load on inclined lines [37]. Green’s function of displacement and tractive force (assuming that normal direction is known) is presented below:(13)gux,ξl=∫-∞∞gupx,ξl,k+guhx,ξl,k+gurx,ξl,ke-ikxdk,gtx,ξl=∫-∞∞gtpx,ξl,k+gthx,ξl,k+gtrx,ξl,ke-ikxdk,and Green’s function of displacement and tractive force for those unloaded TI layers is described as follows:(14)gux,ξl=∫-∞∞gurx,ξl,ke-ikxdk,gtx,ξl=∫-∞∞gtrx,ξl,ke-ikxdk,where gux,ξl and gtx,ξl are Green’s functions of displacement and tractive force for distributed loads. This applies when the loads are distributed on inclined lines with density q(ξl) in the unit (midpoint ξl). The superscript “p” represents the particular solution in the fixed layer, “h” the homogeneous solution in the fixed layer, “r” the solution to the counterforce in fixed end, and “k” the wave number along the horizontal direction.After determining Green’s function, the total displacement and tractive force of any pointx=(x,z) in the TI-layered ground soil can be expressed as(15)vx=∑l=1Kgux,ξlqξl,(16)tyx=∑l=1Kgtx,ξlqξl,where K is the total unit number. The total number of units (element size) is determined by the wavelength of the elastic wave in the TI-layered ground soil. This is necessary in order to meet the requirement of discrete convergence. Equations (15) and (16) can be understood as the dynamic response (scattered wave field) inside the layered TI half-space that is generated by the movement of the foundation. In addition, they can be simulated as the total dynamic response generated by applying a uniformly distributed load on all of the boundary units of the foundation.


### 2.4. Out-of-Plane Dynamic Stiffness Coefficient

Because the connection between the ground soil and the foundation is completely rigid, the displacements of the points in boundaryS of the foundation are the same. This displacement is equal to the out-of-plane displacement Δ of the rigid foundation (the foundation can only generate out-of-plane displacement under the excitation of out-of-plane load), and the displacements of points in the boundary S of the foundation can be expressed as(17)vx=∑l=1Kgux,ξlqξl=Δ.Equation (17) can also be written as(18)∑l=1Kgux,ξlqξlΔ=1=∑l=1Ngux,ξlΛl,where Λl (l=1~K) represents the density of the uniformly distributed load applied on the lth unit when the foundation has unit displacement. Assuming that the rigid contact boundary conditions in S of the foundation can be independently met at ξl of each unit, an equation set formed byK equations about Λ1,Λ2,…,ΛK will be obtained by combining (18). Λl (l=1~K) can be obtained by solving the linear system of equations. The following will be obtained by substituting Λl into (16) (19)tyx=∑l=1Kgtx,ξlΛlΔ.The resultant force acting on the foundation will be determined by integrating (19) along the boundary(20)Fy=∫Stxds=∫S∑l=1Kgtx,ξlΛlΔds=KyyΔ.Equation (20) represents the relational expression of force and displacement acting on the foundation, while Kyy is the dynamic stiffness coefficient of the foundation(21)Kyy=∫S∑l=1Kgtx,ξlΛlds.


### 2.5. The Effective Input Motion and the Superstructure Response

The effective input motionΔ is the displacement response of the foundation under harmonic wave excitation. The displacement response Δ includes two parts Δ1 and Δ2: the first one is the response without considering the mass M0 of the foundation and the mass Mb of the superstructure and the second one is the additional displacement response that considers the mass inertia force F0 of the foundation and the mass inertia force Fb of the superstructure. In the first part above, the virtual work principle is applied to the problem of displacement field and stress field of scattering and radiation of incident wave to obtain the displacement response Δ1 which has been used in [38]. The displacement response Δ1 is applied to the embedded foundation in [39], identified as(22)Δ1=∫Svfx,z∑l=1Kgtx,ξlΛl-tfx,ξldsKyy,in which vf and tf are the out-of-plane displacement and stress of the free field, respectively, alongS.The displacement response in the second part,Δ2, can be obtained using (23)Δ2=F0+FbKyy.The foundation is a rigid body and the following can be obtained under dynamic load(24)Fb=ω2M0Δ.The shear wall is an elastomer and the inertial force can be obtained by elastic dynamics(25)Fb=ω2Mbtan⁡kbHkbHΔ,in which kb=ω/βb represents the shear wave number of the shear wall.Substituting (23), (24), and (25) into Δ=Δ1+Δ2, the effective input motion can be obtained(26)Δ=Δ11-ω2/KyyM0+Mbtan⁡kbH/KbH.The superstructure responseΔb is the relative displacement between the top of the shear wall and the foundation and can be obtained from [10](27)Δb=Δ1cos⁡kbH-1.The relative effective input motion and superstructure response with respect to the free field (i.e., the specific value of the displacement amplitude of that point before and after the existence of the superstructure) can be represented by the following:(28)Δ¯=Δvf,Δ¯b=Δbvf.


## 2.1. Model

The model shown in Figure1 is embedded in TI-layered ground soil and consists of a rigid semicylinder foundation that supports an elastic shear wall. The TI ground is formed using N layers of horizontal TI soil and TI bedrock (half-space) underneath. In TI soil, the axis of symmetry of the medium is assumed to be normal to the horizontal surface. The TI soil layer is determined by the shear modulus in plane normal to axis of symmetry GhiL, horizontal shear modulus, the shear modulus in planes normal to plane of transverse isotropy GviL, vertical shear modulus [32, 33], and the damping ratio ζiL(i=1~N). The TI bedrock (half-space) is determined by the shear modulus in the horizontal direction GhR; the shear modulus in the vertical direction GvR; and the damping ratio ζR. The superscripts L and R represent the soil layer and the bedrock (half-space). The cross section of the foundation is semicircular, but it extends along the y-axis with radius of a and a mass of M0. The shear wall has the same width as the foundation with a shear wave velocity of βb, a height of H, a mass of Mb per unit length, and a damping ratio ζb. The shear wall has a rigid connection (no slip) between the rigid foundation and the TI-layered ground soil, and the interface between the foundation and the ground soil isS. Incident SH wave inputs from angle ψ in the bedrock surface have a circular frequency of ω. In the case of incoming SH waves, the foundation response (i.e., the effective input motion) and the superstructure response are sought after.Figure 1
Foundation-shear wall model.The main steps for determining the foundation response (i.e., the effective input motion) and the superstructure response include (1) the direct stiffness method which was introduced in literature [34] is adopted to determine the response of the free field in TI ground soil; (2) the indirect boundary element method (IBEM) from literature [35] is used to determine the dynamic stiffness coefficient outside of the embedded foundation plane; (3) the effective input motion is solved for, in other words, the sum of the following two parts of the response which is determined: (a) the response when the mass M0 of foundation and the mass Mb of superstructure are not taken into account; (b) the additional displacement response when the mass inertia forces F0 and Fb of the foundation and superstructure are taken into account; (4) the method in [2] is used to determine the superstructure response; (5) the relative effective input motion and superstructure response with respect to the free field are determined, namely, the specific value of the displacement amplitude of that point before and after the existence of superstructure.


## 2.2. The Response of the Free Field in Layered TI Ground Soil

The response of the free field refers to the site response when there is no foundation or superstructure. Based on the accurate out-of-plane dynamic stiffness matrix for TI-layered medium given in the literature [36], the direct stiffness method introduced in literature [34] is adopted to determine the response of the free field in TI ground soil. For the stiffness matrix of the TI soil layer, we assume that there are ascending and descending SH waves in the layer, and then the displacement at any point in the layer can be expressed as (1)vx,z=ASHexp⁡iktz+BSHexp⁡-iktzexp⁡-ikx=vzexp⁡-ikx,where v(z) can be considered as the amplitude transmitting along the x-direction(2)vz=ASHexp⁡iktz+BSHexp⁡-iktz,where ASH and BSH are the amplitudes of the ascending and descending SH waves in the soil layer, respectively, k=ωcos⁡ψSHL/csL∗ is the wave number along the x-direction, and t=-i1-1/cos2ψSHL and ψSHL are the included angles of the ascending or descending SH waves in the soil layer in the horizontal direction(3)csL∗=cshL∗2cos2ψSHL+csvL∗2sin2⁡ψSHL,which is the complex wave velocity of SH waves transmitting in the ψSHL direction. cshL∗=GhL∗/ρ and csvL∗=GvL∗/ρ are the complex wave velocities of SH waves in a soil layer transmitting in the horizontal and vertical directions, GhL∗=GhL1+2iζL and GvL∗=GvL1+2iζL are the complex shear moduli in the horizontal and vertical directions, GhL and GvL are the horizontal and vertical shear moduli when there is no damping in the soil layer, and ζL is the damping ratio.By substitutingz=0 and z=d (d is the thickness of soil layer) into (2), the displacement amplitudes v1 and v2 of the top and bottom of the soil layer can be obtained as(4)v1v2=11exp⁡iktdexp⁡-iktdASHBSH.Through the relational expression of stress and displacement, the shear stressesτyx and τyz can be obtained as(5)τyx=GhL∗∂v∂x=iktGhL∗ASHexp⁡iktz-BSHexp⁡-iktzexp⁡-ikx,(6)τyz=GvL∗∂v∂z=iktGvL∗ASHexp⁡iktz-BSHexp⁡-iktzexp⁡-ikx.By substitutingz=0 and z=d into (6) and introducing the amplitudes of the external load, Q1=τyz(0) and Q2=τyz(d), the amplitudes of the external load on the top and at the bottom of the soil layer can be obtained as(7)Q1Q2=iktGvL∗-11exp⁡iktd-exp⁡-iktdASHBSH.Eliminating coefficientsASH and BSH in (4) and (7), the layered stiffness matrix SSHL of the TI soil layer is (8)Q1Q2=SSHLv1v2=ωsin⁡ψSHLGvL∗csL∗sinωsin⁡ψSHL/csL∗cosωsin⁡ψSHLcsL∗d-1-1cosωsin⁡ψSHLcsL∗dv1v2.Under an infinite radiation condition, the outgoing wave with an amplitude ofBSH will be generated if load is arranged on the surface of TI half-space. The subscript 0 represents the free surface of the TI half-space in (4) and (7). Setting ASH=0, v0=v1, and Q0=Q1=-τyz0, the stiffness coefficient SSHR of the TI half-space will be obtained after the elimination of BSH(9)Q0=SSHRv0=iGvR∗ωsin⁡ψSHRcsR∗v0.After determining all of the stiffness matrices in the TI soil layer and the TI half-space, the accurate dynamic stiffness matrixSSH of the TI-layered site can be obtained, by combining SSHL and SSHR(10)SSH=k11L1k12L1k21L1k22L1+k11L2k12L2k21L2k22L2+k11L3⋱k22LN-1+k11LNk12LNk21LNk22LN+SSHRN+1×N+1,where kijLn (n=1~N,ij=1~2) denotes the ijth submatrix of the nth TI layer’s dynamic stiffness matrix SSHL. Therefore, the discrete dynamic equilibrium equation in the TI-layered site can be expressed as(11)Q1,Q2,…,QnT=SSHv1,v2,…,vnT,where Q1,Q2,…,QnT is the vector of the amplitude of an external load acting on the interface of the soil layer and v1,v2,…,vnT is the vector of the displacement amplitude on the interface of the soil layer.After obtaining the overall dynamic stiffness matrix of TI-layered site, the control point for SH wave transmitting from the bedrock surface can be selected on the outcrop of the bedrock, so that the other elements are zero in the vectorQ1,Q2,…,QnT of amplitude of external load, only(12)Qn=SSHRv0=iGvR∗ωsin⁡ψSHRcsR∗v0,where GvR∗ is the complex shear modulus of the bedrock half-space in the vertical direction, ψSHR is the incident angle of SH waves in the bedrock half-space, and csR∗ is the complex wave velocity of SH waves of the bedrock transmitting in the incident direction. In (3), v0 represents the movement of the bedrock outcrop. Substituting (12) into (11), the displacement in the interface of the TI soil layer can be obtained, and the amplitude coefficient of the ascending and descending wave in any soil layer can be obtained using (4). Finally, the amplitude of displacement and stress in arbitrary point of the TI-layered half-space can be obtained using (1), (5), and (6), namely, the free field response of SH waves.


## 2.3. Green’s Function of the Uniformly Distributed Inclined Load Acting in the TI Medium

In this section, the Green’s function of the uniformly distributed inclined load acting in the TI medium is derived, including both the displacement and stress, which serve the elementary solution of IBEM to obtain the stiffness coefficient of the foundation. The specific deviation of Green’s function is described as follows.Firstly, the load uniformly distributed on inclined lines in the spatial domain is expanded to the wave number domain. Secondly, two fictitious interfaces are introduced at the upper and lower boundaries of the loaded layer. In order to derive the dynamic response of the loaded layer, the two fictitious interfaces have to be fixed to obtain the counterforces of the two ends. Actually, the dynamic response of the loaded layer includes the particular solutions induced by the load applied, homogeneous solutions caused by the reaction loads and the solution derived from the reverse counterforces at the two ends. However, for those unloaded TI layers, their solutions only contain those homogeneous ones induced by the reverse counterforces at the two ends of the loaded layer. Finally, the solutions for every TI layer in the spatial domain are obtained by superposing a Fourier inversion on those solutions in the wave number domain, that is, Green’s function of uniformly distributed load on inclined lines [37]. Green’s function of displacement and tractive force (assuming that normal direction is known) is presented below:(13)gux,ξl=∫-∞∞gupx,ξl,k+guhx,ξl,k+gurx,ξl,ke-ikxdk,gtx,ξl=∫-∞∞gtpx,ξl,k+gthx,ξl,k+gtrx,ξl,ke-ikxdk,and Green’s function of displacement and tractive force for those unloaded TI layers is described as follows:(14)gux,ξl=∫-∞∞gurx,ξl,ke-ikxdk,gtx,ξl=∫-∞∞gtrx,ξl,ke-ikxdk,where gux,ξl and gtx,ξl are Green’s functions of displacement and tractive force for distributed loads. This applies when the loads are distributed on inclined lines with density q(ξl) in the unit (midpoint ξl). The superscript “p” represents the particular solution in the fixed layer, “h” the homogeneous solution in the fixed layer, “r” the solution to the counterforce in fixed end, and “k” the wave number along the horizontal direction.After determining Green’s function, the total displacement and tractive force of any pointx=(x,z) in the TI-layered ground soil can be expressed as(15)vx=∑l=1Kgux,ξlqξl,(16)tyx=∑l=1Kgtx,ξlqξl,where K is the total unit number. The total number of units (element size) is determined by the wavelength of the elastic wave in the TI-layered ground soil. This is necessary in order to meet the requirement of discrete convergence. Equations (15) and (16) can be understood as the dynamic response (scattered wave field) inside the layered TI half-space that is generated by the movement of the foundation. In addition, they can be simulated as the total dynamic response generated by applying a uniformly distributed load on all of the boundary units of the foundation.


## 2.4. Out-of-Plane Dynamic Stiffness Coefficient

Because the connection between the ground soil and the foundation is completely rigid, the displacements of the points in boundaryS of the foundation are the same. This displacement is equal to the out-of-plane displacement Δ of the rigid foundation (the foundation can only generate out-of-plane displacement under the excitation of out-of-plane load), and the displacements of points in the boundary S of the foundation can be expressed as(17)vx=∑l=1Kgux,ξlqξl=Δ.Equation (17) can also be written as(18)∑l=1Kgux,ξlqξlΔ=1=∑l=1Ngux,ξlΛl,where Λl (l=1~K) represents the density of the uniformly distributed load applied on the lth unit when the foundation has unit displacement. Assuming that the rigid contact boundary conditions in S of the foundation can be independently met at ξl of each unit, an equation set formed byK equations about Λ1,Λ2,…,ΛK will be obtained by combining (18). Λl (l=1~K) can be obtained by solving the linear system of equations. The following will be obtained by substituting Λl into (16) (19)tyx=∑l=1Kgtx,ξlΛlΔ.The resultant force acting on the foundation will be determined by integrating (19) along the boundary(20)Fy=∫Stxds=∫S∑l=1Kgtx,ξlΛlΔds=KyyΔ.Equation (20) represents the relational expression of force and displacement acting on the foundation, while Kyy is the dynamic stiffness coefficient of the foundation(21)Kyy=∫S∑l=1Kgtx,ξlΛlds.


## 2.5. The Effective Input Motion and the Superstructure Response

The effective input motionΔ is the displacement response of the foundation under harmonic wave excitation. The displacement response Δ includes two parts Δ1 and Δ2: the first one is the response without considering the mass M0 of the foundation and the mass Mb of the superstructure and the second one is the additional displacement response that considers the mass inertia force F0 of the foundation and the mass inertia force Fb of the superstructure. In the first part above, the virtual work principle is applied to the problem of displacement field and stress field of scattering and radiation of incident wave to obtain the displacement response Δ1 which has been used in [38]. The displacement response Δ1 is applied to the embedded foundation in [39], identified as(22)Δ1=∫Svfx,z∑l=1Kgtx,ξlΛl-tfx,ξldsKyy,in which vf and tf are the out-of-plane displacement and stress of the free field, respectively, alongS.The displacement response in the second part,Δ2, can be obtained using (23)Δ2=F0+FbKyy.The foundation is a rigid body and the following can be obtained under dynamic load(24)Fb=ω2M0Δ.The shear wall is an elastomer and the inertial force can be obtained by elastic dynamics(25)Fb=ω2Mbtan⁡kbHkbHΔ,in which kb=ω/βb represents the shear wave number of the shear wall.Substituting (23), (24), and (25) into Δ=Δ1+Δ2, the effective input motion can be obtained(26)Δ=Δ11-ω2/KyyM0+Mbtan⁡kbH/KbH.The superstructure responseΔb is the relative displacement between the top of the shear wall and the foundation and can be obtained from [10](27)Δb=Δ1cos⁡kbH-1.The relative effective input motion and superstructure response with respect to the free field (i.e., the specific value of the displacement amplitude of that point before and after the existence of the superstructure) can be represented by the following:(28)Δ¯=Δvf,Δ¯b=Δbvf.


## 3. Verification of Accuracy

The method used in this study is verified through comparison with the effective input motion of isotropic media and the superstructure response given in [2, 23]. According to our method, the results in isotropic ground soil can be calculated by using the same shear modulus in the horizontal and vertical directions (Gv=Gh). In order to describe the stiffness of the shear wall, a dimensionless parameter, ε=βLH/βba, is defined in which H and βb are the parameters of the shear wall. We set βR/βL=2,ε=2 in Figures 2 and 3 to perform the numerical calculation and the results obtained in this study are compared to existing ones. It can be seen from the figures that our results correlate well with the results obtained by Trifunac [2] and Liang et al. [23].Figure 2
Comparison between study results and literature results [23].
(a)
(b)Figure 3
Comparison between study results and literature results [2].


## 4. Numerical Results and Analysis

### 4.1. Single TI Soil Layer

Taking the semicylindrical rigid foundation embedded in single TI soil layer on isotropic bedrock as an example, the foundation extends infinitely in they-axis direction with the same cross section. The superstructure is simplified to be a piece of shear wall with the same width, and there is a rigid connection between the rigid foundation and TI soil (no slip). Figures 4, 5, 6, and 7 show the effective input motion and the superstructure response when the shear modulus ratio in the vertical and horizontal directions of TI soil is different. During the calculation, the shear modulus ratio is taken as Gv/Gh = 0.5, 1.0, and 2.0. The bedrock is considered to be isotropic with a density ratio of ρR/ρL=1.0, and the damping ratio of the soil layer and bedrock is ζL=0.05 and ζR=0.02, respectively. The thickness of TI soil layers is D/a=2, 3, 4, 5, the mass of the upper shear wall is Mb/M0=1, 2, 4, and the rigid of the shear wall is ε= 0, 2, 4. To maintain consistency in the vibration frequency of the soil under three different situations, the dimensionless frequency is still defined as ωa/(GV+Gh)/2ρ.Figure 4
The effective input motion for different flexibility and mass of superstructure and forD/a=2.Figure 5
The effective input motion for different thicknesses of soil layer.Figure 6
The superstructure response for different flexibility and mass of superstructure and forD/a=2.Figure 7
The superstructure response for different thicknesses of soil layer.


#### 4.1.1. The Effective Input Motion

As shown in Figure4, it can be found that the effective input motion |Δ-| in single TI soil layer (Gv/Gh = 0.5 and 2.0) is completely different from that in single isotropic layer (Gv/Gh = 1.0) over bedrock. The effective input motion in single TI soil layer fluctuates along with the changing effective input motion in single isotropic layer, and an oscillation phenomenon occurs. According to [23], for ε=2 and D/a=2, the resonant frequencies ωa/(GV+Gh)/2ρ of the shear wall are 0.785,2.356,3.927,…; for ε=4 and D/a=2, the resonant frequencies ωa/(GV+Gh)/2ρ of the shear wall are 0.393,1.178,1.964,…; and so on. As the flexibility (ε) of the shear wall increases, the resonant frequency points of the shear wall concentrate gradually. The effective input motion at a resonant frequency point of the shear wall is zero, and it reaches its maximum value at the first inherent frequency point of the soil layer. This is consistent with research on the effective input motion of a semicylindrical embedded foundation in [2]. In addition, with the increase inε, the peak value of the effective input motion increases and the corresponding frequency of its peak point decreases, and the peak value tends to concentrate at the same frequency point. Some numerical results for D/a=2 and Mb/M0=1 are selected as shown in Table 1. In Table 1, “first frequency” is the first peak frequency of effective input motion and “peak value” is the corresponding peak value.Table 1
Numerical results of effective input motion forD/a=2 and Mb/M0=1.
ε
G
v
/
G
h
First frequency
Peak value
0
0.5
0.541
1.159
1.0
0.636
1.219
2.0
0.711
1.286
2
0.5
0.536
1.410
1.0
0.606
1.697
2.0
0.636
1.988
4
0.5
0.370
3.081
1.0
0.375
3.425
2.0
0.380
3.477Meanwhile, with the increase in the shear modulus ratio of the soil layer in the vertical and horizontal directions, the corresponding frequency of peak value of the effective input motion tends to be higher. This trend is identical to the research on the resonance characteristics of SH waves traveling to a TI-layered site [40]. In addition, the corresponding frequency of the peak value of the effective input motion moves to a higher frequency as the shear modulus ratio in the vertical and horizontal directions increases. In addition, the peak value of the effective input motion also increases, which illustrates that the increase in the shear modulus ratio of the TI soil layer in the vertical and horizontal directions enlarges the response of the effective input motion. For example, for ε=2, D/a=2, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak value of effective input motion is 1.928, 2.216, and 2.470. Moreover, with an increase in the mass of the upper shear wall structure, the peak value of the effective input motion increases gradually, and the maximum peak value appears at the first inherent frequency point of the TI soil layer. At this time, the parameters of the TI soil layer also have a significant impact on the effective input motion. With the increase in shear modulus in the vertical and horizontal directions, the peak value of the effective input motion soars up and the corresponding frequency moves to be higher. The author believes that the change in the parameters of the TI soil layer caused the change in the inherent frequency of the TI soil layer, which further leads to the peak change in the effective input motion.In addition, as the thickness of the TI soil layer changes (D/a=2, 3, 4, and 5), the effective input motion fluctuates in the frequency domain. When the thickness of the soil layer increases, the effective input motion at the resonant frequency point of the shear wall is still zero, but the peak value at the inherent frequency point of the soil layer decreases significantly, and its value tends to be the same. In summary, the greater the thickness of the TI soil layer is, the more significant the fluctuation of the effective input motion will be. Furthermore, as the shear modulus ratio of the soil layer in the vertical and horizontal directions increases, the corresponding frequency of peak point of the effective input motion still moves to be higher. For example, for ε=2, D/a=1, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak frequency of effective input motion is 0.536, 0.606, and 0.636.


#### 4.1.2. The Superstructure Response

As shown in Figure6, the superstructure response |Δ-b| in single TI soil layer (Gv/Gh = 0.5 and 2.0) that can be found via calculation is quite different from that of the superstructure in single isotropic soil layer (Gv/Gh = 1.0). In addition, there is a certain vibration phenomenon that occurs in the low frequency range. The impact of soil TI parameters on the superstructure response in single TI soil layer is obviously greater than that in single isotropic soil layer. There are significant differences between the corresponding structural responses of foundations with different TI parameters, and the impact mainly occurs in the low frequency range. As shown in Figure 6, as the shear modulus ratio of the TI soil layer rises in the vertical and horizontal directions, the corresponding frequency of the peak point of the superstructure response is enhanced. Meanwhile, the value of the peak point gradually increases, which illustrates that the increase in the shear modulus ratio in the vertical and horizontal directions of the TI soil layer increases the response of the superstructure. For example, for ε=2, D/a=4, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak value of the superstructure response is 1.460, 2.045, and 2.506. In addition, when an increasing frequency, the impact of the parameters of the TI medium on the superstructure response weakens but tends to stabilize in the end.The characteristics of the shear wall also have an influence on its relative displacement with respect to the foundation. As the mass of the upper shear wallMb increases, the effective input motion increases gradually, but the superstructure response weakens significantly. As a result, the corresponding frequency of its peak point augments. Therefore, it is clear that a change in the mass of the superstructure will lead to a change in the seismic energy distribution between the foundation and the superstructure. In addition, the superstructure response decreases when the effective input motion increases. Also, the flexibility (ε) of the shear wall itself has a significant impact on the superstructure response. As the flexibility of the shear wall increases, the peak value of the superstructure response significantly increases, and the corresponding frequency of its peak point decreases. Moreover, as the thickness of the TI soil layer changes, the superstructure response fluctuates violently.


### 4.2. Multiple TI Soil Layers

Utilizing a semicylindrical foundation (embedded in TI soil) on an isotropic bedrock as an example, the foundation extends infinitely in they-axis direction with the same cross section. The superstructure is simplified to be a piece of shear wall structure with the same width, and it establishes a rigid connection between the rigid foundation and the TI soil (no slip). Figures 8 and 9 show the effective input motion and the superstructure response. Considering the characteristics of natural soil, there are two conditions selected in multiple TI soil layers: normal-sequence and mixed-sequence. For normal-sequence ground soil, the shear modulus ratio of TI soil layers in the vertical (horizontal) direction is Gv1L(Gh1L) : Gv2L(Gh2L) : Gv3L(Gh3L) : Gv4L(Gh4L) = 1 : 2 : 3 : 4; for mixed-sequence ground soil, the shear modulus ratio of TI soil layers in the vertical (horizontal) direction is Gv1L(Gh1L) : Gv2L(Gh2L) : Gv3L(Gh3L) : Gv4L(Gh4L) = 2 : 1 : 3 : 4 and 1 : 3 : 2 : 4. The medium density of soil layers is the same, and the shear modulus ratio of soil layers in a normal-sequence and mixed-sequence ground soil in the vertical and horizontal directions is GviL/GhiL=0.5 (i=1~4). The shear modulus of equivalent single TI soil layer (Gv1L(Gh1L) : Gv2L(Gh2L) : Gv3L(Gh3L) : Gv4L(Gh4L) = 1 : 1 : 1 : 1) in the vertical and horizontal directions shall be solved according to the equivalent shear wave velocity in the corresponding direction.Figure 8
The effective input motion in multiple TI soil layers (ε=2, Mb/M0=1).Figure 9
The superstructure response in multiple TI soil layers (ε=2, Mb/M0=1).The bedrock is isotropic, the mass of shear wall isMb/M0=1 with rigidity of ε=2, and the density of the bedrock is the same as that of each soil layer. The damping ratio of each soil layer is taken as ζiL(ζ-L)=0.05 (i=1~4), and that of the bedrock is ζR=0.02. The dimensionless frequency is defined as ωa/(G-vL+G-hL)/2ρL.


#### 4.2.1. The Effective Input Motion

From Figure8, we can see that there is a difference between the effective input motion in multiple TI soil layers and that in equivalent single TI soil layer. As the soil layer count increases, the superstructure response will exhibit more wave phenomenon. Furthermore, the sequence of soil layers has a significant impact on the effective input motion. The effective input motion in mixed-sequence layered soil is obviously larger than that in normal-sequence layered soil. This illustrates that the effective input motion increases because of a softer layer.


#### 4.2.2. Superstructure Response

It can be seen from Figure9 that there is a distinction between the superstructure response in multiple TI soil layers and that in equivalent single TI soil layer. In addition, the soil layer sequence has a significant impact on the superstructure response. This indicates that when calculating the superstructure response, the sequence of the actual soil layer shall be considered and shall not be solved simply based on normal-sequence ground soil and single TI soil layer. In general, the peak value of the superstructure response in mixed-sequence TI-layered soil, in single TI soil layer, and in normal-sequence TI-layered soil is successively reduced, while the frequency at the peak value point increases.


## 4.1. Single TI Soil Layer

Taking the semicylindrical rigid foundation embedded in single TI soil layer on isotropic bedrock as an example, the foundation extends infinitely in they-axis direction with the same cross section. The superstructure is simplified to be a piece of shear wall with the same width, and there is a rigid connection between the rigid foundation and TI soil (no slip). Figures 4, 5, 6, and 7 show the effective input motion and the superstructure response when the shear modulus ratio in the vertical and horizontal directions of TI soil is different. During the calculation, the shear modulus ratio is taken as Gv/Gh = 0.5, 1.0, and 2.0. The bedrock is considered to be isotropic with a density ratio of ρR/ρL=1.0, and the damping ratio of the soil layer and bedrock is ζL=0.05 and ζR=0.02, respectively. The thickness of TI soil layers is D/a=2, 3, 4, 5, the mass of the upper shear wall is Mb/M0=1, 2, 4, and the rigid of the shear wall is ε= 0, 2, 4. To maintain consistency in the vibration frequency of the soil under three different situations, the dimensionless frequency is still defined as ωa/(GV+Gh)/2ρ.Figure 4
The effective input motion for different flexibility and mass of superstructure and forD/a=2.Figure 5
The effective input motion for different thicknesses of soil layer.Figure 6
The superstructure response for different flexibility and mass of superstructure and forD/a=2.Figure 7
The superstructure response for different thicknesses of soil layer.


### 4.1.1. The Effective Input Motion

As shown in Figure4, it can be found that the effective input motion |Δ-| in single TI soil layer (Gv/Gh = 0.5 and 2.0) is completely different from that in single isotropic layer (Gv/Gh = 1.0) over bedrock. The effective input motion in single TI soil layer fluctuates along with the changing effective input motion in single isotropic layer, and an oscillation phenomenon occurs. According to [23], for ε=2 and D/a=2, the resonant frequencies ωa/(GV+Gh)/2ρ of the shear wall are 0.785,2.356,3.927,…; for ε=4 and D/a=2, the resonant frequencies ωa/(GV+Gh)/2ρ of the shear wall are 0.393,1.178,1.964,…; and so on. As the flexibility (ε) of the shear wall increases, the resonant frequency points of the shear wall concentrate gradually. The effective input motion at a resonant frequency point of the shear wall is zero, and it reaches its maximum value at the first inherent frequency point of the soil layer. This is consistent with research on the effective input motion of a semicylindrical embedded foundation in [2]. In addition, with the increase inε, the peak value of the effective input motion increases and the corresponding frequency of its peak point decreases, and the peak value tends to concentrate at the same frequency point. Some numerical results for D/a=2 and Mb/M0=1 are selected as shown in Table 1. In Table 1, “first frequency” is the first peak frequency of effective input motion and “peak value” is the corresponding peak value.Table 1
Numerical results of effective input motion forD/a=2 and Mb/M0=1.
ε
G
v
/
G
h
First frequency
Peak value
0
0.5
0.541
1.159
1.0
0.636
1.219
2.0
0.711
1.286
2
0.5
0.536
1.410
1.0
0.606
1.697
2.0
0.636
1.988
4
0.5
0.370
3.081
1.0
0.375
3.425
2.0
0.380
3.477Meanwhile, with the increase in the shear modulus ratio of the soil layer in the vertical and horizontal directions, the corresponding frequency of peak value of the effective input motion tends to be higher. This trend is identical to the research on the resonance characteristics of SH waves traveling to a TI-layered site [40]. In addition, the corresponding frequency of the peak value of the effective input motion moves to a higher frequency as the shear modulus ratio in the vertical and horizontal directions increases. In addition, the peak value of the effective input motion also increases, which illustrates that the increase in the shear modulus ratio of the TI soil layer in the vertical and horizontal directions enlarges the response of the effective input motion. For example, for ε=2, D/a=2, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak value of effective input motion is 1.928, 2.216, and 2.470. Moreover, with an increase in the mass of the upper shear wall structure, the peak value of the effective input motion increases gradually, and the maximum peak value appears at the first inherent frequency point of the TI soil layer. At this time, the parameters of the TI soil layer also have a significant impact on the effective input motion. With the increase in shear modulus in the vertical and horizontal directions, the peak value of the effective input motion soars up and the corresponding frequency moves to be higher. The author believes that the change in the parameters of the TI soil layer caused the change in the inherent frequency of the TI soil layer, which further leads to the peak change in the effective input motion.In addition, as the thickness of the TI soil layer changes (D/a=2, 3, 4, and 5), the effective input motion fluctuates in the frequency domain. When the thickness of the soil layer increases, the effective input motion at the resonant frequency point of the shear wall is still zero, but the peak value at the inherent frequency point of the soil layer decreases significantly, and its value tends to be the same. In summary, the greater the thickness of the TI soil layer is, the more significant the fluctuation of the effective input motion will be. Furthermore, as the shear modulus ratio of the soil layer in the vertical and horizontal directions increases, the corresponding frequency of peak point of the effective input motion still moves to be higher. For example, for ε=2, D/a=1, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak frequency of effective input motion is 0.536, 0.606, and 0.636.


### 4.1.2. The Superstructure Response

As shown in Figure6, the superstructure response |Δ-b| in single TI soil layer (Gv/Gh = 0.5 and 2.0) that can be found via calculation is quite different from that of the superstructure in single isotropic soil layer (Gv/Gh = 1.0). In addition, there is a certain vibration phenomenon that occurs in the low frequency range. The impact of soil TI parameters on the superstructure response in single TI soil layer is obviously greater than that in single isotropic soil layer. There are significant differences between the corresponding structural responses of foundations with different TI parameters, and the impact mainly occurs in the low frequency range. As shown in Figure 6, as the shear modulus ratio of the TI soil layer rises in the vertical and horizontal directions, the corresponding frequency of the peak point of the superstructure response is enhanced. Meanwhile, the value of the peak point gradually increases, which illustrates that the increase in the shear modulus ratio in the vertical and horizontal directions of the TI soil layer increases the response of the superstructure. For example, for ε=2, D/a=4, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak value of the superstructure response is 1.460, 2.045, and 2.506. In addition, when an increasing frequency, the impact of the parameters of the TI medium on the superstructure response weakens but tends to stabilize in the end.The characteristics of the shear wall also have an influence on its relative displacement with respect to the foundation. As the mass of the upper shear wallMb increases, the effective input motion increases gradually, but the superstructure response weakens significantly. As a result, the corresponding frequency of its peak point augments. Therefore, it is clear that a change in the mass of the superstructure will lead to a change in the seismic energy distribution between the foundation and the superstructure. In addition, the superstructure response decreases when the effective input motion increases. Also, the flexibility (ε) of the shear wall itself has a significant impact on the superstructure response. As the flexibility of the shear wall increases, the peak value of the superstructure response significantly increases, and the corresponding frequency of its peak point decreases. Moreover, as the thickness of the TI soil layer changes, the superstructure response fluctuates violently.


## 4.1.1. The Effective Input Motion

As shown in Figure4, it can be found that the effective input motion |Δ-| in single TI soil layer (Gv/Gh = 0.5 and 2.0) is completely different from that in single isotropic layer (Gv/Gh = 1.0) over bedrock. The effective input motion in single TI soil layer fluctuates along with the changing effective input motion in single isotropic layer, and an oscillation phenomenon occurs. According to [23], for ε=2 and D/a=2, the resonant frequencies ωa/(GV+Gh)/2ρ of the shear wall are 0.785,2.356,3.927,…; for ε=4 and D/a=2, the resonant frequencies ωa/(GV+Gh)/2ρ of the shear wall are 0.393,1.178,1.964,…; and so on. As the flexibility (ε) of the shear wall increases, the resonant frequency points of the shear wall concentrate gradually. The effective input motion at a resonant frequency point of the shear wall is zero, and it reaches its maximum value at the first inherent frequency point of the soil layer. This is consistent with research on the effective input motion of a semicylindrical embedded foundation in [2]. In addition, with the increase inε, the peak value of the effective input motion increases and the corresponding frequency of its peak point decreases, and the peak value tends to concentrate at the same frequency point. Some numerical results for D/a=2 and Mb/M0=1 are selected as shown in Table 1. In Table 1, “first frequency” is the first peak frequency of effective input motion and “peak value” is the corresponding peak value.Table 1
Numerical results of effective input motion forD/a=2 and Mb/M0=1.
ε
G
v
/
G
h
First frequency
Peak value
0
0.5
0.541
1.159
1.0
0.636
1.219
2.0
0.711
1.286
2
0.5
0.536
1.410
1.0
0.606
1.697
2.0
0.636
1.988
4
0.5
0.370
3.081
1.0
0.375
3.425
2.0
0.380
3.477Meanwhile, with the increase in the shear modulus ratio of the soil layer in the vertical and horizontal directions, the corresponding frequency of peak value of the effective input motion tends to be higher. This trend is identical to the research on the resonance characteristics of SH waves traveling to a TI-layered site [40]. In addition, the corresponding frequency of the peak value of the effective input motion moves to a higher frequency as the shear modulus ratio in the vertical and horizontal directions increases. In addition, the peak value of the effective input motion also increases, which illustrates that the increase in the shear modulus ratio of the TI soil layer in the vertical and horizontal directions enlarges the response of the effective input motion. For example, for ε=2, D/a=2, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak value of effective input motion is 1.928, 2.216, and 2.470. Moreover, with an increase in the mass of the upper shear wall structure, the peak value of the effective input motion increases gradually, and the maximum peak value appears at the first inherent frequency point of the TI soil layer. At this time, the parameters of the TI soil layer also have a significant impact on the effective input motion. With the increase in shear modulus in the vertical and horizontal directions, the peak value of the effective input motion soars up and the corresponding frequency moves to be higher. The author believes that the change in the parameters of the TI soil layer caused the change in the inherent frequency of the TI soil layer, which further leads to the peak change in the effective input motion.In addition, as the thickness of the TI soil layer changes (D/a=2, 3, 4, and 5), the effective input motion fluctuates in the frequency domain. When the thickness of the soil layer increases, the effective input motion at the resonant frequency point of the shear wall is still zero, but the peak value at the inherent frequency point of the soil layer decreases significantly, and its value tends to be the same. In summary, the greater the thickness of the TI soil layer is, the more significant the fluctuation of the effective input motion will be. Furthermore, as the shear modulus ratio of the soil layer in the vertical and horizontal directions increases, the corresponding frequency of peak point of the effective input motion still moves to be higher. For example, for ε=2, D/a=1, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak frequency of effective input motion is 0.536, 0.606, and 0.636.


## 4.1.2. The Superstructure Response

As shown in Figure6, the superstructure response |Δ-b| in single TI soil layer (Gv/Gh = 0.5 and 2.0) that can be found via calculation is quite different from that of the superstructure in single isotropic soil layer (Gv/Gh = 1.0). In addition, there is a certain vibration phenomenon that occurs in the low frequency range. The impact of soil TI parameters on the superstructure response in single TI soil layer is obviously greater than that in single isotropic soil layer. There are significant differences between the corresponding structural responses of foundations with different TI parameters, and the impact mainly occurs in the low frequency range. As shown in Figure 6, as the shear modulus ratio of the TI soil layer rises in the vertical and horizontal directions, the corresponding frequency of the peak point of the superstructure response is enhanced. Meanwhile, the value of the peak point gradually increases, which illustrates that the increase in the shear modulus ratio in the vertical and horizontal directions of the TI soil layer increases the response of the superstructure. For example, for ε=2, D/a=4, Mb/M0=2, and Gv/Gh = 0.5, 1.0, and 2.0, the first peak value of the superstructure response is 1.460, 2.045, and 2.506. In addition, when an increasing frequency, the impact of the parameters of the TI medium on the superstructure response weakens but tends to stabilize in the end.The characteristics of the shear wall also have an influence on its relative displacement with respect to the foundation. As the mass of the upper shear wallMb increases, the effective input motion increases gradually, but the superstructure response weakens significantly. As a result, the corresponding frequency of its peak point augments. Therefore, it is clear that a change in the mass of the superstructure will lead to a change in the seismic energy distribution between the foundation and the superstructure. In addition, the superstructure response decreases when the effective input motion increases. Also, the flexibility (ε) of the shear wall itself has a significant impact on the superstructure response. As the flexibility of the shear wall increases, the peak value of the superstructure response significantly increases, and the corresponding frequency of its peak point decreases. Moreover, as the thickness of the TI soil layer changes, the superstructure response fluctuates violently.


## 4.2. Multiple TI Soil Layers

Utilizing a semicylindrical foundation (embedded in TI soil) on an isotropic bedrock as an example, the foundation extends infinitely in they-axis direction with the same cross section. The superstructure is simplified to be a piece of shear wall structure with the same width, and it establishes a rigid connection between the rigid foundation and the TI soil (no slip). Figures 8 and 9 show the effective input motion and the superstructure response. Considering the characteristics of natural soil, there are two conditions selected in multiple TI soil layers: normal-sequence and mixed-sequence. For normal-sequence ground soil, the shear modulus ratio of TI soil layers in the vertical (horizontal) direction is Gv1L(Gh1L) : Gv2L(Gh2L) : Gv3L(Gh3L) : Gv4L(Gh4L) = 1 : 2 : 3 : 4; for mixed-sequence ground soil, the shear modulus ratio of TI soil layers in the vertical (horizontal) direction is Gv1L(Gh1L) : Gv2L(Gh2L) : Gv3L(Gh3L) : Gv4L(Gh4L) = 2 : 1 : 3 : 4 and 1 : 3 : 2 : 4. The medium density of soil layers is the same, and the shear modulus ratio of soil layers in a normal-sequence and mixed-sequence ground soil in the vertical and horizontal directions is GviL/GhiL=0.5 (i=1~4). The shear modulus of equivalent single TI soil layer (Gv1L(Gh1L) : Gv2L(Gh2L) : Gv3L(Gh3L) : Gv4L(Gh4L) = 1 : 1 : 1 : 1) in the vertical and horizontal directions shall be solved according to the equivalent shear wave velocity in the corresponding direction.Figure 8
The effective input motion in multiple TI soil layers (ε=2, Mb/M0=1).Figure 9
The superstructure response in multiple TI soil layers (ε=2, Mb/M0=1).The bedrock is isotropic, the mass of shear wall isMb/M0=1 with rigidity of ε=2, and the density of the bedrock is the same as that of each soil layer. The damping ratio of each soil layer is taken as ζiL(ζ-L)=0.05 (i=1~4), and that of the bedrock is ζR=0.02. The dimensionless frequency is defined as ωa/(G-vL+G-hL)/2ρL.


### 4.2.1. The Effective Input Motion

From Figure8, we can see that there is a difference between the effective input motion in multiple TI soil layers and that in equivalent single TI soil layer. As the soil layer count increases, the superstructure response will exhibit more wave phenomenon. Furthermore, the sequence of soil layers has a significant impact on the effective input motion. The effective input motion in mixed-sequence layered soil is obviously larger than that in normal-sequence layered soil. This illustrates that the effective input motion increases because of a softer layer.


### 4.2.2. Superstructure Response

It can be seen from Figure9 that there is a distinction between the superstructure response in multiple TI soil layers and that in equivalent single TI soil layer. In addition, the soil layer sequence has a significant impact on the superstructure response. This indicates that when calculating the superstructure response, the sequence of the actual soil layer shall be considered and shall not be solved simply based on normal-sequence ground soil and single TI soil layer. In general, the peak value of the superstructure response in mixed-sequence TI-layered soil, in single TI soil layer, and in normal-sequence TI-layered soil is successively reduced, while the frequency at the peak value point increases.


## 4.2.1. The Effective Input Motion

From Figure8, we can see that there is a difference between the effective input motion in multiple TI soil layers and that in equivalent single TI soil layer. As the soil layer count increases, the superstructure response will exhibit more wave phenomenon. Furthermore, the sequence of soil layers has a significant impact on the effective input motion. The effective input motion in mixed-sequence layered soil is obviously larger than that in normal-sequence layered soil. This illustrates that the effective input motion increases because of a softer layer.


## 4.2.2. Superstructure Response

It can be seen from Figure9 that there is a distinction between the superstructure response in multiple TI soil layers and that in equivalent single TI soil layer. In addition, the soil layer sequence has a significant impact on the superstructure response. This indicates that when calculating the superstructure response, the sequence of the actual soil layer shall be considered and shall not be solved simply based on normal-sequence ground soil and single TI soil layer. In general, the peak value of the superstructure response in mixed-sequence TI-layered soil, in single TI soil layer, and in normal-sequence TI-layered soil is successively reduced, while the frequency at the peak value point increases.


## 5. Conclusions

In this study, the IBEM method is applied to study the effective input motion and the response of the shear wall under the incidence of SH waves transmitting into layered TI soil. The accuracy of this method is verified by comparing the obtained results with existing research. The semicylindrical foundation-shear wall structure system in single TI soil layer and multiple TI soil layers is utilized as experimental examples. The impact of TI characteristics of soil layers and a shear wall structure on the effective input motion and the superstructure response is studied. The following conclusions are made:(1)
The study on the effective input motion and the superstructure response in single TI soil layer shows that the TI characteristics of the ground soil have a significant impact on the effective input motion and the superstructure response. As the shear modulus ratio of the TI soil layer in the vertical and horizontal directions increases, the peak value of the effective input motion and the superstructure response also increases to a certain degree, and their corresponding peak frequency increases. The increasing flexibility of the upper shear wall structure also increases the effective input motion and the peak value of the superstructure response. In addition, as the mass of superstructure increases, the effective input motion increases while the superstructure response significantly decreases. It is clear that a change in the mass of the superstructure will lead to a change in the seismic energy distribution between the foundation and the superstructure.(2)
The study on the effective input motion and the superstructure response in multiple TI soil layers shows that there is a significant difference between the effective input motion and the superstructure response in multiple TI soil layers and that in single TI soil layer. The difference is related to the sequence of the soil layers, and mixed TI soil layers significantly increase the effective input motion and the superstructure response.



---

*Source: 2834274-2017-06-21.xml*