#### Abstract

The dynamic response of pile embedded in layered transversely isotropic soil and subjected to arbitrary vertical harmonic force is investigated. Based on the viscoelastic constitutive relations for a transversely isotropic medium, the dynamic governing equation of the transversely isotropic soil is obtained in cylindrical coordinates. By introducing the fictitious soil pile model and the distributed Voigt model, the governing equations of soil-pile system are also derived. Firstly, the vertical response of the soil layer is solved by using the Laplace transform technique and the separation of variables technique. Secondly, the analytical solution of velocity response in the frequency domain and its corresponding semianalytical solution of velocity response in the time domain are derived by means of inverse Fourier transform and convolution theorem. Finally, based on the obtained solutions, a parametric study has been conducted to investigate the influence of the soil anisotropy on the vertical dynamic response of pile. It can be seen that the influence of the shear modulus of soil in the vertical plane on the dynamic response of pile is more notable than the influence of the shear modulus of soil in the horizontal plane on the dynamic response of pile.

#### 1. Introduction

The dynamic interaction of soil-pile system is a complicated contact problem owing to the complexity of soil layer. Therefore, many investigators have paid their attention to develop theoretical models to analyze the soil layer for the dynamic interaction of soil-pile system. Many of these various theoretical models fall into four main categories. (1) The dynamic Winkler model [1–7], which represents the soil as a series of distributed springs attached to the pile, is the first category. The stiffness of the soil “springs” is calculated by assuming that wave energy can only propagate outwards under plane-strain conditions [8]. (2) The plane-strain model [9–16], which assumes that the soil is made up of an infinite number of infinitesimally thin, independent, horizontal, elastic layers that extend to infinity, is the second category. It can be noted that the plane-strain model neglects the strain in the vertical direction and the waves propagate only in the horizontal direction. (3) Various three-dimensional axisymmetric continuum models [17–22], in which the soil layer is three-dimensional and the wave effect of soil is considered, is the third category. Due to the complexity of the three-dimensional axisymmetric continuum models, it is very difficult to promote the use of these models in layered soil. (4) The radially inhomogeneous model [23–26], in which the soil region is assumed as a linear viscoelastic medium composed of two concentric regions, an inner annular region of disturbed medium, and an outer semi-infinite undisturbed region, is the fourth category. It can be seen that the model does not simulate the interaction of soil zones rigorously and it is only an approximate simplified model that may lead to the argument whether the results adequately reflect the real interaction of the soil zones [27].

Most of the previous investigations on the dynamic interaction of soil-pile system had treated soil as an isotropic medium. However, in practice soil deposits possess a certain degree of anisotropy owing to their deposition history resulting in properties that are different in the horizontal and vertical direction. The vertical and horizontal differences can be simulated by virtue of a transversely isotropic soil model [28] which a simple isotropic model cannot allow for. Tsai [29] studied the torsional dynamic response of a circular disk on a transversely isotropic half-space by utilizing integral transform technique. It is worth noting that the anisotropic material constants have obvious influence on the resonant amplitude and frequency of vibration. Liu and Novak [30] investigated the dynamic response of single pile embedded in transversely isotropic layered media by means of FEM combining with dynamic stiffness matrices of the soil derived from Green’s functions for ring loads. Chen et al. [31] proposed a model of transversely isotropic layered elastic media to study the dynamic horizontal impedances of double piles embedded in the ground and subjected to harmonic sway-rocking loadings at the head of the two piles. Chen et al. [32] studied the transient torsional dynamic response of a pile embedded in transversely isotropic saturated soil by utilizing the Laplace transform. Wang et al. [33] investigated the dynamic response of an end bearing pile embedded in transversely isotropic saturated soil when there was a time-harmonic torsional loading acting on the pile head. Notwithstanding the work above, it can be seen that no investigation corresponding to the vertical dynamic response of a pile embedded in layered transversely isotropic elastic soil has been reported in existing literatures until now.

Based on the above literature review, the main objective of this paper is to develop an extended soil-pile interaction model to simulate the vertical dynamic response of pile embedded in layered transversely isotropic soil. Utilizing the fictitious soil pile model to simulate the dynamic interaction between the pile and pile end soil [34], the governing equations of soil-pile system are established when the vertical wave effect of surrounding soil is taken into account. Then, the analytical solution of velocity response in the frequency domain and its corresponding semianalytical solution of velocity response in the time domain are derived by means of integral transform technique and separation of variables technique. Based on these solutions, a parametric study has been conducted to assess the influence of the soil anisotropy on the vibration behavior of the pile.

#### 2. Governing Equations

##### 2.1. Computational Model and Assumptions

The problem studied in this paper is the vertical vibration of viscoelastic pile embedded in layered transversely isotropic soil and the geometric model is shown in Figure 1. Based on the fictitious soil pile model, the soil-pile system is discretized into a total of segments along the vertical direction, which are numbered by . The thickness of the th () soil-pile segment is denoted by and the depth of the th () soil layer top is denoted by . If the numbers of the soil-pile segments are enough, the properties of soil and pile within each segment can be assumed to be homogeneous, respectively. The pile length is denoted by and the thickness of pile end soil is denoted by . is an arbitrary vertical harmonic force acting on the pile head.

The following assumptions are adopted during the analysis.(1)The surrounding soil of pile is layered, transversely isotropic and viscoelastic. The damping force of soil is proportional to the strain rate, and the proportional coefficient of the th soil layer is denoted by .(2)The top surface of soil layer is assumed to be free boundary without normal and shear stresses, and the bottom surface of pile end soil is assumed to be rigid boundary without displacements.(3)The dynamic interaction of the adjacent soil layers is simulated by using the distributed Voigt model. The spring constant and damping coefficient of the distributed Voigt model between the th soil layer and its upper adjacent soil layer are denoted by and , and the corresponding values between the th soil layer and its lower adjacent soil layer are denoted by and , respectively.(4)During vibration, the vertical wave effect of surrounding soil is taken into account, but the radial wave effect of surrounding soil is ignored.(5)Both the pile and the fictitious soil pile are vertical, viscoelastic, and circular in cross-section and have a perfect contact with the surrounding soil during vibration.(6)During vibration, the soil-pile system is subjected to small deformations and strains, and the conditions of displacement continuity and force equilibrium are satisfied at the interface of the adjacent pile (including fictitious soil pile) segments.

##### 2.2. Governing Equations of Soil-Pile System

###### 2.2.1. Dynamic Equation of Soil

Combining with the viscoelastic constitutive relations for a transversely isotropic medium proposed by Ding [35], the dynamic equilibrium equation of the transversely isotropic soil in cylindrical coordinates can be derived as follows where , , and denote the radial displacement, circumferential displacement, and vertical displacement, respectively. The coefficients of the above equations should satisfy the following equations where and denote the horizontal and vertical elastic modulus, respectively. denotes the shear modulus in the vertical plane. is the Poisson’s ration in the vertical direction caused by the horizontal stress and is the Poisson’s ratio in the horizontal direction caused by the vertical stress and and should satisfy . is the Poisson’s ratio in the orthogonal direction of the horizontal strain caused by the horizontal stress.

Owing to the assumption that only the vertical wave effect of surrounding soil is taken into consideration, the equilibrium equation of the transversely isotropic soil for the axisymmetric problem can be further rewritten as follows

Based on (3) and taking into account the viscosity of soil, the governing equation of the th transversely isotropic soil layer for the axisymmetric problem can be established as where is the vertical displacement of the th soil layer. , , and denote the viscous damping coefficient, density, and shear modulus in the vertical plane of the th soil layer. . denotes the ratio of the elastic modulus of the th soil layer in the horizontal direction to that in the vertical direction; , , and denote the shear modulus in the horizontal plane, the Poisson’s ratio in the horizontal direction caused by the vertical stress, and the Poisson’s ratio in the orthogonal direction of the horizontal strain caused by the horizontal stress of the th soil layer, respectively.

###### 2.2.2. Dynamic Equation of Pile

Denoting to be the vertical displacement of the th pile (including fictitious soil pile) segment and according to the Euler-Bernoulli rod theory, the dynamic equilibrium equation of pile can be established as where , , , , , , and denote the elastic modulus, cross-section area, radius, mass per unit length of pile (including fictitious soil pile), density, elastic longitudinal wave velocity, and viscous damping coefficient of the th pile (including fictitious soil pile) segment, respectively. is the frictional force of the th soil layer acting on the surface of the th pile shaft and can be expressed as

Combining with the assumptions, the boundary and initial conditions of soil-pile system can be established as follows.

*(**1) Boundary Conditions of Soil*. At the top surface of the th soil layer,

At the bottom surface of the th soil layer,

At an infinite radial distance of the th soil layer,

*(2) Boundary Conditions of Pile*. At the top surface of the th pile (including fictitious soil pile) segment,

At the bottom surface of the th pile (including fictitious soil pile) segment,

where and denote the displacement impedance function at the top and bottom surfaces of the th pile (including fictitious soil pile) segment. is the Laplace transform parameter.

*(3) Boundary Condition at the Interface of Soil-Pile System*

* (**4) Initial Conditions of Soil-Pile System*. Initial conditions of the th soil layer are as follows:

Initial conditions of the th pile (including fictitious soil pile) segment are as follows:

#### 3. Solutions of the Governing Equations

##### 3.1. Vibrations of the Soil Layer

Denoting to be the Laplace transform with respect to time of and associating with the initial condition (10a), (4) is transformed by using the Laplace transform technique as follows where is the shear wave velocity of the th soil layer in the vertical direction.

By virtue of the separation of variables technique and denoting , (11) can be decoupled as follows where constants and must satisfy the following equation

It can be seen that (12) is Bessel equation and (13) is ordinary differential equation of second order whose general solutions can be easily obtained. Associating with these general solutions, the vertical displacement of the th soil layer, can be derived as where and denote the modified Bessel functions of order zero of the first and second kind, respectively. , , , and are constants determined by the boundary conditions.

Converting and in the global coordinates into and in the local coordinates, (7a), (7b), and (7c) are transformed by using the Laplace transform technique and can be rewritten as follows

According to the properties of the modified Bessel functions, that is, if , then and , it can be obtained that from boundary conditions (16c). By means of boundary conditions (16a) and (16b), the following equation can be obtained where and denote the dimensionless complex stiffness of the upper surface and lower surface of the th soil layer. Substituting into (17) and solving it by using bisection method in the frequency domain, a series of eigenvalues, , can be derived. can also be derived by substituting into (14).

Then, the vertical displacement of the th soil layer can be rewritten as where and are a series of constants determined by the boundary conditions which reflect the dynamic interaction of soil and pile.

##### 3.2. Vibrations of the Pile

Denoting to be the Laplace transform with respect to time of , (5) is transformed by using the Laplace transform technique and can be rewritten combining with (6) and (18)

It is not difficult to obtain that the general solution of (19) can be expressed as where where , , , and are all dimensionless parameters. denotes the propagation time of elastic longitudinal wave in the th pile segment. is the modified Bessel functions of order one of the second kind.

Combining with the boundary conditions (8a) and (8b), the displacement impedance function at the head of the th pile segment can be derived in the local coordinates as follows where where denotes the displacement impedance function at the head of the ()th pile segment which can be obtained by using boundary conditions. Then, following the method of recursion typically used in the transfer function technique, the displacement impedance function at the head of pile can be derived as where where , , , and are all dimensionless parameters. denotes the propagation time of elastic longitudinal wave in the th pile segment. and can be obtained from the following equations where and denote the dimensionless complex stiffness of the upper surface and lower surface of the th soil layer.

Then, the velocity transfer function at the head of pile can be obtained as

Letting and substituting it into (27), the velocity response in the frequency domain at the head of pile can be obtained as where is the dimensionless velocity admittance at the pile head which can be expressed as

By means of the inverse Fourier transform and convolution theorem, the velocity response in the time domain at the head of pile can be expressed as , where denotes the Fourier transform of , which is the vertical excitation acting on the pile head.

In particular, the excitation can be regarded as a half-sine pulse in the nondestructive detection of pile foundation as follows where and denote the duration of the impulse or impulse width and the maximum amplitude of the vertical excitation, respectively. Then, the velocity response in the time domain at the head of pile can be expressed as where is the dimensionless velocity response which can be expressed as where denotes the dimensionless frequency. denotes the dimensionless pulse width which should satisfy . denotes the dimensional time variable which should satisfy .

#### 4. Analysis of Vibration Characteristics

According to the derivation process shown in the previous section, it can be seen that the difference between the shear modulus in the vertical plane and the shear modulus in the horizontal plane reflects the soil anisotropy. Therefore, based on the solutions, the influence of these two kinds of shear modulus of pile surrounding soil and pile end soil on the dynamic response of pile is studied in detail. Unless otherwise specified, the length, radius, density, and longitudinal wave velocity of pile are , , , and , respectively. The spring constant of the distributed Voigt model is equal to the elastic modulus of the lower soil layer, and the damping coefficient of the distributed Voigt model is ··.

##### 4.1. Influence of the Anisotropy of Pile Surrounding Soil on the Dynamic Response of Pile

Firstly, the influence of the shear modulus of pile surrounding soil in the vertical plane on the dynamic response of pile is investigated. Parameters of pile end soil are as follows: the thickness is three times that of pile diameter; the soil density is ; both the shear modulus in the vertical plane and the shear modulus in the horizontal plane are ; both the Poisson’s ratio in the horizontal direction caused by the vertical stress and the Poisson’s ratio in the orthogonal direction of the horizontal strain caused by the horizontal stress are 0.35; the damping coefficient is ··. Parameters of pile surrounding soil are as follows: the soil density is ; both the Poisson’s ratio in the horizontal direction caused by the vertical stress and the Poisson’s ratio in the orthogonal direction of the horizontal strain caused by the horizontal stress are 0.4; the damping coefficient is ·; the shear modulus in the horizontal plane is ; the shear modulus in the vertical plane is , respectively.

Figure 2 shows the influence of the shear modulus of pile surrounding soil in the vertical plane on the dynamic response of pile. As shown in Figure 2(a), it can be noted that the velocity admittance curves oscillate about a mean amplitude as the frequency increases. As the shear modulus of pile surrounding soil in the vertical plane increases, the amplitude of resonance peaks gradually decreases, but the resonance frequency of velocity admittance curves almost remains unchanged. As shown in Figure 2(b), it is observed that the amplitude of the incident pulses and reflective wave signals decreases with the increase of the shear modulus of pile surrounding soil in the vertical plane. As the shear modulus of pile surrounding soil in the vertical plane increases, the raising phenomenon between the incident pulses and the primary reflective wave signals will be gradually aggravated, and the declining phenomenon between the primary reflective wave signals and the secondary reflective wave signals will also be gradually intensified.

**(a) Velocity admittance curves**

**(b) Reflected wave signal curve**

After that, the influence of the shear modulus of pile surrounding soil in the horizontal plane on the dynamic response of pile is studied. Parameters of pile surrounding soil are as follows: the shear modulus in the vertical plane is , and the shear modulus in the horizontal plane is , respectively. The other parameters of soil-pile system are the same as those shown in the previous case.

Figure 3 shows the influence of the shear modulus of pile surrounding soil in the horizontal plane on the dynamic response of pile. As shown in Figure 3(a), it can be seen that the amplitude of resonance peaks gradually increases with the increase of the shear modulus of pile surrounding soil in the horizontal plane, but the change of resonance frequency can be ignored. As shown in Figure 3(b), it can be seen that the amplitude of the reflective wave signals decreases with the increase of the shear modulus of pile surrounding soil in the horizontal plane, but the decreasing ratio is small. The raising phenomenon between the primary reflective wave signals and the secondary reflective wave signals will be gradually intensified with the increase of the shear modulus of pile surrounding soil in the horizontal plane.

**(a) Velocity admittance curves**

**(b) Reflected wave signal curves**

##### 4.2. Influence of the Anisotropy of Pile End Soil on the Dynamic Response of Pile

In this section, the influence of the shear modulus of pile end soil in the vertical plane on the dynamic response of pile is firstly investigated. Parameters of pile surrounding soil are as follows: the soil density is ; both the shear modulus in the horizontal plane and the shear modulus in the vertical plane are ; both the Poisson’s ratio in the horizontal direction caused by the vertical stress and the Poisson’s ratio in the orthogonal direction of the horizontal strain caused by the horizontal stress are 0.4; the damping coefficient is ··. Parameters of pile end soil are as follows: the thickness is three times that of pile diameter; the soil density is ; both the Poisson’s ratio in the horizontal direction caused by the vertical stress and the Poisson’s ratio in the orthogonal direction of the horizontal strain caused by the horizontal stress are 0.35; the damping coefficient is ··; the shear modulus in the horizontal plane is ; the shear modulus in the vertical plane is , , , , , respectively.

Figure 4 shows the influence of the shear modulus of pile end soil in the vertical plane on the dynamic response of pile. As shown in Figure 4(a), as the shear modulus of pile end soil in the vertical plane increases, the amplitude of resonance peaks gradually decreases with the same resonance frequency, but the decreasing ratio is small. As shown in Figure 4(b), it can be seen that the amplitude of the reflective wave signals decreases with the increase of the shear modulus of pile surrounding soil in the vertical plane.

**(a) Velocity admittance curves**

**(b) Reflected wave signal curves**

Then, the influence of the shear modulus of pile end soil in the horizontal plane on the dynamic response of pile is studied. Parameters of pile end soil are as follows: the shear modulus in the vertical plane is , and the shear modulus in the horizontal plane is , , , , , respectively. The other parameters of soil-pile system are the same as those shown in the previous case.

Figure 5 shows the influence of the shear modulus of pile end soil in the horizontal plane on the dynamic response of pile. It can be seen that the influence of the shear modulus of pile end soil in the horizontal plane on the dynamic response of pile can be ignored.

**(a) Velocity admittance curves**

**(b) Reflected wave signal curves**

#### 5. Conclusions

By considering a pile embedded in layered transversely isotropic soil as a dynamic soil-pile interaction problem, the governing equations of soil-pile system are established when there is arbitrary vertical harmonic force acting on the pile head. Then, an analytical solution for the velocity response in the frequency domain and its corresponding semianalytical solution for the velocity response in the time domain have been derived by virtue of the transform technique and the separation of variables technique. An extensive parameter study has been undertaken to investigate the influence of the soil anisotropy on the vertical dynamic response of pile and the following conclusions have been obtained.(1)Whether for the pile surrounding soil or for the pile end soil, it can be seen that the influence of the shear modulus in the vertical plane on the dynamic response of pile is more notable than the influence of the shear modulus in the horizontal plane on the dynamic response of pile. Therefore, the shear modulus of soil in the vertical plane plays a leading role in the dynamic response of pile when only the vertical wave effect of soil is taken into account.(2)As the shear modulus of pile surrounding soil in the vertical plane increases, both the amplitude of the resonance peaks of velocity admittance curves and the reflective wave signals of reflected wave signal curves gradually decrease. As the shear modulus of pile surrounding soil in the horizontal plane increases, the amplitude of the resonance peaks of velocity admittance curves gradually increases, but the reflective wave signals of reflected wave signal curves gradually decrease.(3)As the shear modulus of pile end soil in the vertical plane increases, both the amplitude of the resonance peaks of velocity admittance curves and the reflective wave signals of reflected wave signal curves gradually decrease. The influence of the shear modulus of pile end soil in the horizontal plane on the dynamic response of pile can be ignored.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This research is supported by the National Natural Science Foundation of China (Grant no. 51309207), the China Postdoctoral Science Foundation Funded Project (Grant nos. 2012M521495 and 2013T60759), and the Fundamental Research Funds for the Central Universities (Grant no. CUG120821).