# Dynamic Response of an Inhomogeneous Viscoelastic Pile in a Multilayered Soil to Transient Axial Loading

Authors: Zhiqing Zhang; Jian Zhou; Kuihua Wang; Qiang Li; Kaifu Liu
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

A quasi-analytical solution is developed in this paper to investigate the mechanism of one-dimensional longitudinal wave propagating in inhomogeneous viscoelastic pile embedded in layered soil and subjected to a transient axial loading. At first, the pile-soil system is subdivided into several layers along the depth direction in consideration of the variation of cross-sectional acoustic impedance of the pile or differences in soil properties. Then, the dynamic governing equation of arbitrary soil layer is established in cylindrical coordinates and arbitrary viscoelastic pile segment is modeled using a single Voigt model. By using the Laplace transform and boundary conditions of the pile-soil system, the vertical impedance at the top of arbitrary pile segment is defined in a closed form in the frequency domain. Then by utilizing the method of recursion typically used in the Transfer Function technique, the vertical impedance at the pile top can be derived in the frequency domain and the velocity response of an inhomogeneous viscoelastic pile subjected to a semi-sine wave exciting force is obtained in a semi-analytical form in the time domain. Selected numerical results are obtained to study the mechanism of longitudinal wave propagating in a pile with a single defect or double defects.


---

## Body

## 1. Introduction

Pile vibration theory can provide valuable guidance for both the dynamic design of embedded foundations and dynamic nondestructive integrity testing of piles. For the dynamic design of pile foundations, the most concerned problem is the study of vibration characteristics of embedded piles in the low frequency range, and accordingly the dynamic reaction from the surrounding soil is an essential component in the relevant study. In light of this, many soil models, such as simplified continuum model [1–3], continuum model [4–6], and finite element model [7], have been developed to investigate the pile-soil dynamic interaction. For the dynamic nondestructive integrity testing of piles, low strain integrity testing technique has received wide application in assessing the construction quality of piles due to its relatively low-cost and simplicity. This testing technique is mainly based on the stress wave propagation theory through a bar. As a result, the theoretical study of a pile embedded in the soil and subjected to a dynamic vertical load under small deformation condition is increasingly becoming more important. For instance, Davis and Dunn [8] firstly proposed the mechanical admittance (or mobility function) method to determine the length and cross-sectional area of a pile by applying a steady-state harmonic excitation at a specified set of frequencies on the pile top. In subsequent studies, Davis and Robertson [9] extended the mechanical admittance method to determine the pile head stiffness. Accounting for the convenience of transient excitation, Higgs [10] modified the mechanical admittance method by applying a vertical impact load on the pile top instead of steady-state harmonic excitation. After that, Lin et al. [11] introduced the impact-echo method to assess the integrity of piles by using the amplitude Fourier spectrum of the displacement record at the pile head instead of the mobility function. Watson et al. [12] developed a wavelet transform signal processing method instead of the conventional Fourier based methods to locate the position of the pile tip. It is noted that, in the low strain integrity pile testing where the hammer is relatively small compared to the pile dimension, Rayleigh and shear waves will radiate from the impact loading and the effects of three-dimensional (3D) waves on the near field responses are obvious. Subsequently, several researchers conducted relevant studies and proposed several methods to diminish the effects of 3D waves on the dynamic response of the pile top. Liao and Roesset [13, 14] investigated the influence of 3D waves on the dynamic response at the top of intact and defective pile by comparing one-dimensional (1D) wave theory and 3D axisymmetric finite element simulation results. It is shown from their studies that 3D effects are mainly influenced by the frequency and are more strongly manifested at high frequencies. Chow et al. [15] found that the velocity response curves resemble that of a pile with a defect near the pile head when considering 3D effects and further proposed that the potential source of error can be removed by maintaining a distance between hammer and receiver that is greater than 50% of the pile radius. Chai et al. [16] found that when the ratio of the characteristic length of an impact pulse to the pile radius is large enough, the components of Rayleigh waves in the wave field at the pile top are diminished. In this study, Chai et al. still proposed that the receiver should be placed at positions between 0.5R and 0.75R (R = pile radius) from the pile axis to diminish the influence of the multireflections. Lu et al. [17] investigated the 3D characteristics of wave propagation in pipe-pile using elastodynamic finite integration technique and found that the interferences of Rayleigh waves are weakest at an angle of 90° from where hammer hits. Furthermore, for the drilled piles with high slenderness ratio, it is difficult to detect the pile length and deep flaw from the traditional low strain pile integrity testing technique due to insufficient impact energy, testing signal decay, and soil-pile interaction. To solve this problem, Ni et al. [18] adjusted the testing devices for acquiring a lower frequency signal and developed a new numerical signal process method to enhance the reflection signals from the pile tip. It is also shown from the experimental results that the testing signal identification abilities can be improved by the modified method.Most of the previous studies on the low strain pile integrity testing did not consider the effect of pile material damping, abrupt variation of surrounding soil properties, and multidefects in pile on the dynamic response. It is worth noting that the material damping indeed exists in a pile, the properties of the surrounding soil may change greatly in certain embedment depth, and a pile may contain several defects. Wang et al. [19] investigated the vertical dynamic response of an inhomogeneous viscoelastic pile and analyzed the effect of pile material damping and soil properties on the mechanical admittance and velocity response of the pile top. However, in this study, the surrounding soil reaction on the pile is approximately simulated by a general Voigt model which cannot veritably and accurately reflect the pile-soil interaction. Therefore, the objective of this paper is to develop a practical solution to evaluate the theoretical capabilities of the nondestructive dynamic response method in detecting the existence and location of single or double defects in a viscoelastic pile embedded in a multilayered soil. Using the solution developed, a parametric study has been undertaken to investigate the mechanism of one-dimensional elastic longitudinal wave propagating in a defective pile. Finally, the theoretical model developed in the present paper is validated by comparison of the theoretical fitted curve and field measured curve of velocity response.


## 2. Formulation of the Problem

### 2.1. Geometry and Basic Assumption

The system examined is an inhomogeneous viscoelastic pile embedded in a multilayered soil and the geometric model is shown in Figure1. To portray the variation of cross-sectional acoustic impedance (the product of density, cross-sectional area, and the one-dimensional elastic longitudinal wave velocity) of a pile or differences in soil properties, the pile-soil system is subdivided into a total of N segments (layers) numbered by 1,2,…,j,…,N from pile tip to pile head. The thickness of the jth (1≤j≤N) soil layer is equal to the length of the jth (1≤j≤N) pile segment and is denoted by hj. In order to derive an analytical or quasi-analytical solution for this problem, the assumptions are made as follows: (1) the surrounding soil is a linearly viscoelastic layer and the pile is vertical, elastic, and circular in cross-section. The pile and soil layer properties are assumed to be homogeneous within each segment or layer, respectively, but may change from segment to segment or layer to layer; (2) the pile-soil system is subjected to small deformations and strains during the vibration; the pile has a perfect contact with the surrounding soil during the vibration; (3) the free surface of the soil has no normal and shear stresses and the soil is infinite in the radial direction; (4) the soil at the base of the pile is modeled using a spring with elastic constant kpb and a dashpot with damping coefficient cpb; (5) the contact traction acting at the jth soil layer due to j-1th and j+1th soil layers is treated as the distributed Winkler subgrade model independent of the radial distance (see Figure 2).Figure 1
Model of pile-soil interaction.Figure 2
Pile segment used in derivation of impedance (in the local coordinate system).


### 2.2. Basic Equations and Solutions

The axisymmetric dynamic response of a linear viscoelastic soil is considered here and the governing equation of thejth soil layer can be expressed in terms of vertical displacement as follows [20]:(1)ηj2+iDvjηj2-2+2Dsj∂2wjr,z,t∂z2+1+iDsj1r∂∂r+∂2∂r2wjr,z,t=1Vj2∂2wjr,z,t∂t2,where ηj=(λj+2Gj)/Gj=2(1-μj)/(1-2μj); Vj=Gj/ρj is the shear wave velocity in the jth soil layer; μj is Poisson ratio; λj and Gj are Lame’s elastic constant and shear modulus of the jth soil layer, respectively; λj′ and Gj′ are corresponding viscosity coefficients about Lame’s elastic constant and shear modulus, respectively; Dsj=Gj′/Gj and Dvj=λj′/λj are hysteretic type damping which is irrespective with frequency; wj(r,z,t) is the vertical displacement of the jth soil layer; i=-1.If the pile material is modeled as a single Voigt model, the axial displacement of thejth pile segment is governed by the following one-dimensional equation of motion:(2)Epjρpj∂2ujz,t∂z2+ηpjρpj∂3ujz,t∂z2∂t+2πrjρpjApjfjz,t=∂2ujz,t∂t2,where rj, ρpj, Epj, Apj, ηpj, and uj(z,t) are the radius, mass density, Young’s modulus, cross-sectional area, viscoelastic damping coefficient, and axial displacement of the jth pile segment, respectively; fj(z,t) denotes the contact traction along the jth pile-soil interface due to the surrounding soil.


### 2.3. Initial Conditions of the Pile-Soil System

The pile-soil system is stationary in the initial state(t=0). Accordingly, the corresponding initial conditions can be written as(3)wjr,z,tt=0=0;∂wjr,z,t∂tt=0=0;j=1,2,…,N,(4)ujz,tt=0=0∂ujz,t∂tt=0=0;j=1,2,…,N.


### 2.4. Boundary and Continuity Conditions of the Pile-Soil System

Equations (1) and (2) cannot be solved analytically in the time domain; thus they are solved by using Laplace transform. In order to solve the pile-soil interaction problem, the boundary and continuity conditions of the pile-soil system will be formulated in the Laplace transform space and expressed in the local coordinate system (as shown in Figure 2). The Laplace transform of a function f(r,z,t) with respect to t is defined as(5)f∗r,z,s=∫0∞fr,z,te-stdt,where s is the Laplace transform parameter.


#### 2.4.1. The Boundary and Continuity Conditions of the Soil Layer

The boundary conditions of thejth (1≤j≤N) soil layer can be written in the local coordinate system as(6)wj∗r⟶∞,z,s=0,(7)∂wj∗r,z,s∂z-kstjEsjwj∗r,z,sz=0=0(8)∂wj∗r,z,s∂z+ksbjEsjwj∗r,z,sz=hj=0,where kstj and ksbj denote the distributed reaction coefficient of the j+1th and j-1th soil layer acting on the jth soil layer, respectively. For the Nth soil layer, kstN is equal to zero due to the free surface of the soil.The continuity condition of the interface between thejth and j+1th (1≤j≤N-1) soil layer can be written as(9)kstj=ksbj+1.


#### 2.4.2. The Boundary Conditions of the Pile Segment

The boundary condition at the top of thejth (1≤j≤N) pile segment can be expressed in the local coordinate system as(10)duj∗z,sdz+ηpjsEpjduj∗z,sdzz=0=-Pj∗sEpjApj,where Pj∗(s) denotes the axial force of j+1th pile segment acting on the top of the jth pile segment. It is noted that PN∗(s)=P∗(s) for the Nth pile segment.Allowing for the continuity conditions of the vertical displacement and axial force at the interface of adjacent pile segments, the boundary condition at the base of thejth (1≤j≤N) pile segment can then be expressed in the local coordinate system as(11)duj∗z,sdz+ηpjsEpjduj∗z,sdz+uj∗z,sZj-1EpjApjz=hj=0,where Zj-1 denotes the vertical impedance at the top of the j-1th pile segment. It is worth noting that Z0=kpb+scpb for the 1st pile segment. Moreover, soil response at the pile toe is approximated by the response of a vertically vibrating rigid disk on the surface of an elastic half-space. The soil spring and damping coefficient at the pile toe can be given as follows [21]:(12)kpb=4G0r11-μ0,cpb=3.4r12ρ0G01-μ0,where G0 and μ0 denote the shear modulus and Poisson’s ratio of the soil underlying the pile toe.


#### 2.4.3. The Continuity Conditions of the Displacement and Stress of thejth (1≤j≤N) Pile-Soil Interface

Consider the following:(13)wj∗r,z,sr=rj=uj∗z,s,(14)fj∗z,s=τrzj∗r,z,s=Gj1+iDsj∂wj∗r,z,s∂rr=rj.


### 2.5. Solution of thejth Soil Layer

Combining (1) with the initial conditions given in (3) and applying the Laplace transform yield(15)ηj2+iDvjηj2-2+2Dsj∂2wj∗r,z,s∂z2+1+iDsj1r∂∂r+∂2∂r2wj∗r,z,s=sVj2wj∗r,z,s.To solve (15), a single-variable function wj∗(r,z,s)=Rj(r)Zj(z) is introduced and substituting wj∗(r,z,s) into (15) yields(16)d2Zjzdz2+Jj2Zjz=0,d2Rjrdr2+1rdRjrdr-qj2Rjr=0,where qj2=(ηj2+iDvj(ηj2-2)+2DsjJj2+(s/Vj)2)/(1+iDsj).The general solutions of (16) can be written as(17)Zjz=Cjsin⁡Jjz+Djcos⁡Jjz,(18)Rjr=AjK0qjr+BjI0qjr,where I0(qjr) and K0(qjr) are the modified Bessel functions of the first and second kinds of zero order, respectively; Aj, Bj, Cj, and Dj are the constants which remain to be determined later from the boundary conditions.Since the functionI0(qjr)→∞ when r→∞, the constant Bj in (18) should vanish to zero to satisfy the boundary condition given in (6). Then, the substitution of boundary conditions given in (7) and (8) into (17) results in(19)tan⁡hjJj=k¯sbj+k¯stjhjJjhjJj2-k¯sbjk¯stj,where k¯sbj=ksbjhj/Esj and k¯stj=kstjhj/Esj denote the dimensionless reaction coefficients at the base and top of the jth soil layer, respectively.Then, the solution of (15) can be written in a series expansion as(20)wj∗r,z,s=∑m=1∞AmjK0qmjrsin⁡Jmjz+φmj,where φmj=arctan⁡Jmjhj/k¯stj, qmj2=([ηj2+i(Dvj(ηj2-2)+2Dsj)]Jmj2+(s/Vj)2)/(1+iDsj).Based on continuity condition given in (14), the contact traction along the 1st pile-soil interface can be expressed as(21)fj∗r,z,sr=rj=Gj1+iDsj∂wj∗r,z,s∂rr=rj=-Gj1+iDsj·∑m=1∞AmjqmjK1qmjrjsin⁡Jmjz+φmj,where K1(qmjr) denotes the modified Bessel functions of the second kind of the first order.


### 2.6. Solution of thejth Pile Segment

Based on the initial conditions given in (4), applying the Laplace transform into (2) yields(22)Epjρpj+ηpjsρpjd2uj∗z,sdz2+2πrjρpjApjfj∗z,s=s2uj∗z,s.For thejth pile segment, by utilizing the stress continuity condition that is given in (14), substituting (21) into (22) gives(23)Epjρpj+ηpjsρpjd2uj∗z,sdz2-s2uj∗z,s=2πrjGj1+iDsjρpjApj·∑m=1∞AmjqmjK1qmjrjsin⁡Jmjz+φmj.The solution of (23) can be written as(24)uj∗z,s=α1jeγjz+α2je-γjz+∑m=1∞ψmjsin⁡Jmjz+φmj,where ψmj=-2Gj(1+iDsj)AmjqmjK1(qmjrj)/(rj[Epj(1+η¯pjstc)Jmj2+ρpjs2]); γj=ρpjs2/[Epj(1+η¯pjstc)]; α1j and α2j are undetermined constants; η¯pj=ηpj/(Epjtc) is the dimensionless damping coefficient of jth pile segment; tc=∑j=1N(hj/Vpj) is the time of the elastic longitudinal wave propagating from the pile top to pile tip; Vpj=Epj/ρpj is the 1D elastic longitudinal wave velocity in the jth pile segment.By using the displacement continuity condition that is given in (13), substituting (20) and (24) into (13) results in(25)α1jeγjz+α2je-γjz+∑m=1∞ψmjsin⁡Jmjz+φmj=∑m=1∞AmjK0qmjrjsin⁡Jmjz+φmj.By invoking the orthogonality of eigenfunctionssin⁡(Jmjz+φmj) (m=1,2,3,…), multiplying sin⁡(Jmjz+φmj) on both sides of (25), and then integrating over the interval z=[0,hj], the undetermined coefficient Amj is found to be(26)Amj=1LmjEmj∫0hjQjsin⁡Jmjz+φmjdz,where Lmj=hj/2-(sin⁡(2Jmjhj+2φmj)-sin⁡(2φmj))/(4Jmj), Qj=α1jeγjz+α2je-γjz, (27)Emj=2Gj1+iDsjqmjK1qmjrjrjEpj1+η¯pjstcJmj2+ρpjs2+K0qmjrj.It is mathematically convenient at this stage to introduce the following dimensionless variables:(28)L¯mj=Lmjhj;q¯mj=qmjhj;r¯j=rjhj;J¯mj=Jmjhj;γ¯j=γjhj.The amplitude of the vertical displacement of thejth pile segment is then given by(29)uj∗z,s=α1jeγjz+∑m=1∞ξmjsin⁡Jmjz+φmj+α2je-γjz+∑m=1∞ζmjsin⁡Jmjz+φmj,where (30)ξmj=vmjγ¯jeγ¯jsin⁡J¯mj+φmj-sin⁡φmj-J¯mjeγ¯jcos⁡J¯mj+φmj-cos⁡φmjJ¯mj2+γ¯j2,ζmj=vmj-γ¯je-γ¯jsin⁡J¯mj+φmj-sin⁡φmj-J¯mje-γ¯jcos⁡J¯mj+φmj-cos⁡φmjJ¯mj2+γ¯j2,vmj=-2Gj1+iDsjq¯mjK1q¯mjr¯jEpj1+η¯pjstcJ¯mj2+ρpjhj2s2r¯jL¯mjEmj.Based on the boundary conditions of thejth pile segment, then substituting (10) and (11) into (29), the variables α1j, α2j are obtained and the vertical impedance function at the top end of the jth pile segment can be written as(31)Zjs=Pjsujz=0,s=-EpjApj+η¯pjstcApjα1j/α2jγj+∑m=1∞ξmjJmjcos⁡φmj+-γj+∑m=1∞ζmjJmjcos⁡φmjα1j/α2j1+∑m=1∞ξmjsin⁡φmj+1+∑m=1∞ζmjsin⁡φmj,where(32)α1jα2j=-1+η¯pjstc-γ¯je-γ¯j+∑m=1∞ζmjJmjcos⁡J¯mj+φmj+Z¯j-1e-γ¯j+∑m=1∞ζmjsin⁡J¯mj+φmj1+η¯pjstcγ¯jeγ¯j+∑m=1∞ξmjJmjcos⁡J¯mj+φmj+Z¯j-1eγ¯j+∑m=1∞ξmjsin⁡J¯mj+φmj,Z¯j-1=Zj-1hjEpjApj.For the 1st pile segment, Z¯0=(kpb+scpb)h1/(Ep1Ap1) denotes the dimensionless reaction coefficient at the pile base.Through recursion impedance function from the 1st pile segment to theNth pile segment, the vertical impedance function at the pile head can be expressed as(33)ZNs=PNsuNz=0,s=-EpNApN+η¯pNstcApNα1N/α2NγN+∑m=1∞ξmNJmNcos⁡φmN+-γN+∑m=1∞ζmNJmNcos⁡φmNα1N/α2N1+∑m=1∞ξmNsin⁡φmN+1+∑m=1∞ζmNsin⁡φmN.Lettings=iω, the frequency response function of displacement (admittance function of displacement) at the pile head can be written as(34)Huω=1ZN=-1EpNApN+η¯pNstcApNα1N/α2N1+∑m=1∞ξmNsin⁡φmN+1+∑m=1∞ζmNsin⁡φmNα1N/α2NγN+∑m=1∞ξmNJmNcos⁡φmN+-γN+∑m=1∞ζmNJmNcos⁡φmN.Then, the response function of velocity (or admittance function of velocity) at the pile head can be written as(35)Hvω=iωHuω=1ρpNApNVpNHv′ω,where Hv′(ω) is the dimensionless response function of velocity at the pile head and(36)Hv′ω=-iωρpNApNVpNEpNApN+η¯pNstcApNα1N/α2N1+∑m=1∞ξmNsin⁡φmN+1+∑m=1∞ζmNsin⁡φmNα1N/α2NγN+∑m=1∞ξmNJmNcos⁡φmN+-γN+∑m=1∞ζmNJmNcos⁡φmN.If the Fourier transform of the longitudinal exciting forceP(t) acting on the pile top is denoted by P(ω), the velocity response of the pile top in the time domain can be calculated by convolution theorem as(37)Vt=IFTPωHω.When the force at the top of the pile is a half-sine pulse,(38)Pt=Qmax⁡sin⁡πtT0t<T00t≥T0,where T0 and Qmax denote the duration of the impulse or impulse width and the maximum amplitude of the exciting force, respectively. Then, the velocity response of the pile top in the time domain can be further expressed as (39)Vt=12π∫-∞+∞Qmax⁡VpNApNρpNHvω·πT0π2-T02ω21+e-iωT0eiωT0dω.To facilitate analysis, it is useful to introduce the normalized velocity response:(40)V′t=Vtmax⁡Vt,where max⁡[V(t)] denotes the maximum value of V(t).


## 2.1. Geometry and Basic Assumption

The system examined is an inhomogeneous viscoelastic pile embedded in a multilayered soil and the geometric model is shown in Figure1. To portray the variation of cross-sectional acoustic impedance (the product of density, cross-sectional area, and the one-dimensional elastic longitudinal wave velocity) of a pile or differences in soil properties, the pile-soil system is subdivided into a total of N segments (layers) numbered by 1,2,…,j,…,N from pile tip to pile head. The thickness of the jth (1≤j≤N) soil layer is equal to the length of the jth (1≤j≤N) pile segment and is denoted by hj. In order to derive an analytical or quasi-analytical solution for this problem, the assumptions are made as follows: (1) the surrounding soil is a linearly viscoelastic layer and the pile is vertical, elastic, and circular in cross-section. The pile and soil layer properties are assumed to be homogeneous within each segment or layer, respectively, but may change from segment to segment or layer to layer; (2) the pile-soil system is subjected to small deformations and strains during the vibration; the pile has a perfect contact with the surrounding soil during the vibration; (3) the free surface of the soil has no normal and shear stresses and the soil is infinite in the radial direction; (4) the soil at the base of the pile is modeled using a spring with elastic constant kpb and a dashpot with damping coefficient cpb; (5) the contact traction acting at the jth soil layer due to j-1th and j+1th soil layers is treated as the distributed Winkler subgrade model independent of the radial distance (see Figure 2).Figure 1
Model of pile-soil interaction.Figure 2
Pile segment used in derivation of impedance (in the local coordinate system).


## 2.2. Basic Equations and Solutions

The axisymmetric dynamic response of a linear viscoelastic soil is considered here and the governing equation of thejth soil layer can be expressed in terms of vertical displacement as follows [20]:(1)ηj2+iDvjηj2-2+2Dsj∂2wjr,z,t∂z2+1+iDsj1r∂∂r+∂2∂r2wjr,z,t=1Vj2∂2wjr,z,t∂t2,where ηj=(λj+2Gj)/Gj=2(1-μj)/(1-2μj); Vj=Gj/ρj is the shear wave velocity in the jth soil layer; μj is Poisson ratio; λj and Gj are Lame’s elastic constant and shear modulus of the jth soil layer, respectively; λj′ and Gj′ are corresponding viscosity coefficients about Lame’s elastic constant and shear modulus, respectively; Dsj=Gj′/Gj and Dvj=λj′/λj are hysteretic type damping which is irrespective with frequency; wj(r,z,t) is the vertical displacement of the jth soil layer; i=-1.If the pile material is modeled as a single Voigt model, the axial displacement of thejth pile segment is governed by the following one-dimensional equation of motion:(2)Epjρpj∂2ujz,t∂z2+ηpjρpj∂3ujz,t∂z2∂t+2πrjρpjApjfjz,t=∂2ujz,t∂t2,where rj, ρpj, Epj, Apj, ηpj, and uj(z,t) are the radius, mass density, Young’s modulus, cross-sectional area, viscoelastic damping coefficient, and axial displacement of the jth pile segment, respectively; fj(z,t) denotes the contact traction along the jth pile-soil interface due to the surrounding soil.


## 2.3. Initial Conditions of the Pile-Soil System

The pile-soil system is stationary in the initial state(t=0). Accordingly, the corresponding initial conditions can be written as(3)wjr,z,tt=0=0;∂wjr,z,t∂tt=0=0;j=1,2,…,N,(4)ujz,tt=0=0∂ujz,t∂tt=0=0;j=1,2,…,N.


## 2.4. Boundary and Continuity Conditions of the Pile-Soil System

Equations (1) and (2) cannot be solved analytically in the time domain; thus they are solved by using Laplace transform. In order to solve the pile-soil interaction problem, the boundary and continuity conditions of the pile-soil system will be formulated in the Laplace transform space and expressed in the local coordinate system (as shown in Figure 2). The Laplace transform of a function f(r,z,t) with respect to t is defined as(5)f∗r,z,s=∫0∞fr,z,te-stdt,where s is the Laplace transform parameter.


### 2.4.1. The Boundary and Continuity Conditions of the Soil Layer

The boundary conditions of thejth (1≤j≤N) soil layer can be written in the local coordinate system as(6)wj∗r⟶∞,z,s=0,(7)∂wj∗r,z,s∂z-kstjEsjwj∗r,z,sz=0=0(8)∂wj∗r,z,s∂z+ksbjEsjwj∗r,z,sz=hj=0,where kstj and ksbj denote the distributed reaction coefficient of the j+1th and j-1th soil layer acting on the jth soil layer, respectively. For the Nth soil layer, kstN is equal to zero due to the free surface of the soil.The continuity condition of the interface between thejth and j+1th (1≤j≤N-1) soil layer can be written as(9)kstj=ksbj+1.


### 2.4.2. The Boundary Conditions of the Pile Segment

The boundary condition at the top of thejth (1≤j≤N) pile segment can be expressed in the local coordinate system as(10)duj∗z,sdz+ηpjsEpjduj∗z,sdzz=0=-Pj∗sEpjApj,where Pj∗(s) denotes the axial force of j+1th pile segment acting on the top of the jth pile segment. It is noted that PN∗(s)=P∗(s) for the Nth pile segment.Allowing for the continuity conditions of the vertical displacement and axial force at the interface of adjacent pile segments, the boundary condition at the base of thejth (1≤j≤N) pile segment can then be expressed in the local coordinate system as(11)duj∗z,sdz+ηpjsEpjduj∗z,sdz+uj∗z,sZj-1EpjApjz=hj=0,where Zj-1 denotes the vertical impedance at the top of the j-1th pile segment. It is worth noting that Z0=kpb+scpb for the 1st pile segment. Moreover, soil response at the pile toe is approximated by the response of a vertically vibrating rigid disk on the surface of an elastic half-space. The soil spring and damping coefficient at the pile toe can be given as follows [21]:(12)kpb=4G0r11-μ0,cpb=3.4r12ρ0G01-μ0,where G0 and μ0 denote the shear modulus and Poisson’s ratio of the soil underlying the pile toe.


### 2.4.3. The Continuity Conditions of the Displacement and Stress of thejth (1≤j≤N) Pile-Soil Interface

Consider the following:(13)wj∗r,z,sr=rj=uj∗z,s,(14)fj∗z,s=τrzj∗r,z,s=Gj1+iDsj∂wj∗r,z,s∂rr=rj.


## 2.4.1. The Boundary and Continuity Conditions of the Soil Layer

The boundary conditions of thejth (1≤j≤N) soil layer can be written in the local coordinate system as(6)wj∗r⟶∞,z,s=0,(7)∂wj∗r,z,s∂z-kstjEsjwj∗r,z,sz=0=0(8)∂wj∗r,z,s∂z+ksbjEsjwj∗r,z,sz=hj=0,where kstj and ksbj denote the distributed reaction coefficient of the j+1th and j-1th soil layer acting on the jth soil layer, respectively. For the Nth soil layer, kstN is equal to zero due to the free surface of the soil.The continuity condition of the interface between thejth and j+1th (1≤j≤N-1) soil layer can be written as(9)kstj=ksbj+1.


## 2.4.2. The Boundary Conditions of the Pile Segment

The boundary condition at the top of thejth (1≤j≤N) pile segment can be expressed in the local coordinate system as(10)duj∗z,sdz+ηpjsEpjduj∗z,sdzz=0=-Pj∗sEpjApj,where Pj∗(s) denotes the axial force of j+1th pile segment acting on the top of the jth pile segment. It is noted that PN∗(s)=P∗(s) for the Nth pile segment.Allowing for the continuity conditions of the vertical displacement and axial force at the interface of adjacent pile segments, the boundary condition at the base of thejth (1≤j≤N) pile segment can then be expressed in the local coordinate system as(11)duj∗z,sdz+ηpjsEpjduj∗z,sdz+uj∗z,sZj-1EpjApjz=hj=0,where Zj-1 denotes the vertical impedance at the top of the j-1th pile segment. It is worth noting that Z0=kpb+scpb for the 1st pile segment. Moreover, soil response at the pile toe is approximated by the response of a vertically vibrating rigid disk on the surface of an elastic half-space. The soil spring and damping coefficient at the pile toe can be given as follows [21]:(12)kpb=4G0r11-μ0,cpb=3.4r12ρ0G01-μ0,where G0 and μ0 denote the shear modulus and Poisson’s ratio of the soil underlying the pile toe.


## 2.4.3. The Continuity Conditions of the Displacement and Stress of thejth (1≤j≤N) Pile-Soil Interface

Consider the following:(13)wj∗r,z,sr=rj=uj∗z,s,(14)fj∗z,s=τrzj∗r,z,s=Gj1+iDsj∂wj∗r,z,s∂rr=rj.


## 2.5. Solution of thejth Soil Layer

Combining (1) with the initial conditions given in (3) and applying the Laplace transform yield(15)ηj2+iDvjηj2-2+2Dsj∂2wj∗r,z,s∂z2+1+iDsj1r∂∂r+∂2∂r2wj∗r,z,s=sVj2wj∗r,z,s.To solve (15), a single-variable function wj∗(r,z,s)=Rj(r)Zj(z) is introduced and substituting wj∗(r,z,s) into (15) yields(16)d2Zjzdz2+Jj2Zjz=0,d2Rjrdr2+1rdRjrdr-qj2Rjr=0,where qj2=(ηj2+iDvj(ηj2-2)+2DsjJj2+(s/Vj)2)/(1+iDsj).The general solutions of (16) can be written as(17)Zjz=Cjsin⁡Jjz+Djcos⁡Jjz,(18)Rjr=AjK0qjr+BjI0qjr,where I0(qjr) and K0(qjr) are the modified Bessel functions of the first and second kinds of zero order, respectively; Aj, Bj, Cj, and Dj are the constants which remain to be determined later from the boundary conditions.Since the functionI0(qjr)→∞ when r→∞, the constant Bj in (18) should vanish to zero to satisfy the boundary condition given in (6). Then, the substitution of boundary conditions given in (7) and (8) into (17) results in(19)tan⁡hjJj=k¯sbj+k¯stjhjJjhjJj2-k¯sbjk¯stj,where k¯sbj=ksbjhj/Esj and k¯stj=kstjhj/Esj denote the dimensionless reaction coefficients at the base and top of the jth soil layer, respectively.Then, the solution of (15) can be written in a series expansion as(20)wj∗r,z,s=∑m=1∞AmjK0qmjrsin⁡Jmjz+φmj,where φmj=arctan⁡Jmjhj/k¯stj, qmj2=([ηj2+i(Dvj(ηj2-2)+2Dsj)]Jmj2+(s/Vj)2)/(1+iDsj).Based on continuity condition given in (14), the contact traction along the 1st pile-soil interface can be expressed as(21)fj∗r,z,sr=rj=Gj1+iDsj∂wj∗r,z,s∂rr=rj=-Gj1+iDsj·∑m=1∞AmjqmjK1qmjrjsin⁡Jmjz+φmj,where K1(qmjr) denotes the modified Bessel functions of the second kind of the first order.


## 2.6. Solution of thejth Pile Segment

Based on the initial conditions given in (4), applying the Laplace transform into (2) yields(22)Epjρpj+ηpjsρpjd2uj∗z,sdz2+2πrjρpjApjfj∗z,s=s2uj∗z,s.For thejth pile segment, by utilizing the stress continuity condition that is given in (14), substituting (21) into (22) gives(23)Epjρpj+ηpjsρpjd2uj∗z,sdz2-s2uj∗z,s=2πrjGj1+iDsjρpjApj·∑m=1∞AmjqmjK1qmjrjsin⁡Jmjz+φmj.The solution of (23) can be written as(24)uj∗z,s=α1jeγjz+α2je-γjz+∑m=1∞ψmjsin⁡Jmjz+φmj,where ψmj=-2Gj(1+iDsj)AmjqmjK1(qmjrj)/(rj[Epj(1+η¯pjstc)Jmj2+ρpjs2]); γj=ρpjs2/[Epj(1+η¯pjstc)]; α1j and α2j are undetermined constants; η¯pj=ηpj/(Epjtc) is the dimensionless damping coefficient of jth pile segment; tc=∑j=1N(hj/Vpj) is the time of the elastic longitudinal wave propagating from the pile top to pile tip; Vpj=Epj/ρpj is the 1D elastic longitudinal wave velocity in the jth pile segment.By using the displacement continuity condition that is given in (13), substituting (20) and (24) into (13) results in(25)α1jeγjz+α2je-γjz+∑m=1∞ψmjsin⁡Jmjz+φmj=∑m=1∞AmjK0qmjrjsin⁡Jmjz+φmj.By invoking the orthogonality of eigenfunctionssin⁡(Jmjz+φmj) (m=1,2,3,…), multiplying sin⁡(Jmjz+φmj) on both sides of (25), and then integrating over the interval z=[0,hj], the undetermined coefficient Amj is found to be(26)Amj=1LmjEmj∫0hjQjsin⁡Jmjz+φmjdz,where Lmj=hj/2-(sin⁡(2Jmjhj+2φmj)-sin⁡(2φmj))/(4Jmj), Qj=α1jeγjz+α2je-γjz, (27)Emj=2Gj1+iDsjqmjK1qmjrjrjEpj1+η¯pjstcJmj2+ρpjs2+K0qmjrj.It is mathematically convenient at this stage to introduce the following dimensionless variables:(28)L¯mj=Lmjhj;q¯mj=qmjhj;r¯j=rjhj;J¯mj=Jmjhj;γ¯j=γjhj.The amplitude of the vertical displacement of thejth pile segment is then given by(29)uj∗z,s=α1jeγjz+∑m=1∞ξmjsin⁡Jmjz+φmj+α2je-γjz+∑m=1∞ζmjsin⁡Jmjz+φmj,where (30)ξmj=vmjγ¯jeγ¯jsin⁡J¯mj+φmj-sin⁡φmj-J¯mjeγ¯jcos⁡J¯mj+φmj-cos⁡φmjJ¯mj2+γ¯j2,ζmj=vmj-γ¯je-γ¯jsin⁡J¯mj+φmj-sin⁡φmj-J¯mje-γ¯jcos⁡J¯mj+φmj-cos⁡φmjJ¯mj2+γ¯j2,vmj=-2Gj1+iDsjq¯mjK1q¯mjr¯jEpj1+η¯pjstcJ¯mj2+ρpjhj2s2r¯jL¯mjEmj.Based on the boundary conditions of thejth pile segment, then substituting (10) and (11) into (29), the variables α1j, α2j are obtained and the vertical impedance function at the top end of the jth pile segment can be written as(31)Zjs=Pjsujz=0,s=-EpjApj+η¯pjstcApjα1j/α2jγj+∑m=1∞ξmjJmjcos⁡φmj+-γj+∑m=1∞ζmjJmjcos⁡φmjα1j/α2j1+∑m=1∞ξmjsin⁡φmj+1+∑m=1∞ζmjsin⁡φmj,where(32)α1jα2j=-1+η¯pjstc-γ¯je-γ¯j+∑m=1∞ζmjJmjcos⁡J¯mj+φmj+Z¯j-1e-γ¯j+∑m=1∞ζmjsin⁡J¯mj+φmj1+η¯pjstcγ¯jeγ¯j+∑m=1∞ξmjJmjcos⁡J¯mj+φmj+Z¯j-1eγ¯j+∑m=1∞ξmjsin⁡J¯mj+φmj,Z¯j-1=Zj-1hjEpjApj.For the 1st pile segment, Z¯0=(kpb+scpb)h1/(Ep1Ap1) denotes the dimensionless reaction coefficient at the pile base.Through recursion impedance function from the 1st pile segment to theNth pile segment, the vertical impedance function at the pile head can be expressed as(33)ZNs=PNsuNz=0,s=-EpNApN+η¯pNstcApNα1N/α2NγN+∑m=1∞ξmNJmNcos⁡φmN+-γN+∑m=1∞ζmNJmNcos⁡φmNα1N/α2N1+∑m=1∞ξmNsin⁡φmN+1+∑m=1∞ζmNsin⁡φmN.Lettings=iω, the frequency response function of displacement (admittance function of displacement) at the pile head can be written as(34)Huω=1ZN=-1EpNApN+η¯pNstcApNα1N/α2N1+∑m=1∞ξmNsin⁡φmN+1+∑m=1∞ζmNsin⁡φmNα1N/α2NγN+∑m=1∞ξmNJmNcos⁡φmN+-γN+∑m=1∞ζmNJmNcos⁡φmN.Then, the response function of velocity (or admittance function of velocity) at the pile head can be written as(35)Hvω=iωHuω=1ρpNApNVpNHv′ω,where Hv′(ω) is the dimensionless response function of velocity at the pile head and(36)Hv′ω=-iωρpNApNVpNEpNApN+η¯pNstcApNα1N/α2N1+∑m=1∞ξmNsin⁡φmN+1+∑m=1∞ζmNsin⁡φmNα1N/α2NγN+∑m=1∞ξmNJmNcos⁡φmN+-γN+∑m=1∞ζmNJmNcos⁡φmN.If the Fourier transform of the longitudinal exciting forceP(t) acting on the pile top is denoted by P(ω), the velocity response of the pile top in the time domain can be calculated by convolution theorem as(37)Vt=IFTPωHω.When the force at the top of the pile is a half-sine pulse,(38)Pt=Qmax⁡sin⁡πtT0t<T00t≥T0,where T0 and Qmax denote the duration of the impulse or impulse width and the maximum amplitude of the exciting force, respectively. Then, the velocity response of the pile top in the time domain can be further expressed as (39)Vt=12π∫-∞+∞Qmax⁡VpNApNρpNHvω·πT0π2-T02ω21+e-iωT0eiωT0dω.To facilitate analysis, it is useful to introduce the normalized velocity response:(40)V′t=Vtmax⁡Vt,where max⁡[V(t)] denotes the maximum value of V(t).


## 3. Problem Definition

The possible defects due to voids or inclusions in a pile after construction can generally be categorized as neck, bulb, and weak concrete. The geometry of typical defective piles is plotted in Figure3. WRN (neck width) and WRB (bulb width) denote the radius of the neck and bulb, respectively. FH and FL denote the embedment depth and length of the defect, respectively. If a pile contains more than one defect (i.e., two defects), the geometry of defective pile can be the combination of Figures 3(a), 3(b), and 3(c).Figure 3
Geometry of typical defective piles: (a) pile with a neck, (b) pile with a bulb, and (c) pile with a weak concrete.


## 4. Results and Discussion

In the past, the mechanical admittance method has received wide application in the field pile integrity test. However, the mechanical admittance method is difficult to accurately identify defects when a pile contains more than one flaw. Furthermore, quick attenuation of the amplitude of admittance curve in the high frequency range caused by the pile material damping can also make the location of the defect hard to be accurately detected. As a result, in previous few years, the method by virtue of time history of velocity response of the pile top has then been widely used instead of the mechanical admittance method in China due to its simplicity of subsequent parameter analysis. Accordingly, the influence of geometric and physical characteristics of the defects and pile material damping on the velocity response of the pile top will be investigated in the following section. Unless otherwise specified, the pile properties employed in the following analysis arej=1,2,…,N, ρpj=2500 kg/m3, Vpj=3790 m/s, η¯pj=0.0015, H=19 m, rj=0.5 m, t0=1.5×10-3 s and Poisson’s ratio = 0.2; the soil properties are ρj=1800 kg/m3, Vsj=120 m/s, Dsj=Dvj=0.02, k¯sbj=k¯stj=0.01 and Poisson’s ratio = 0.3.


### 4.1. Pile Containing a Single Defect

#### 4.1.1. Effect of Width of Defect

Figure4 shows the influence of the neck width (WRN) on the velocity response of the pile top. It can be seen from Figure 4 that the first reflective signal from the defect (RS1) arrives at t=tRS1=0.00300 s (tRS1 is the time of arrival of the first reflective signal from the defect). Accordingly, the embedment depth of the defect can be calculated as FH=Vrod×tRS1/2=5.68 m (Vrod is the 1D elastic longitudinal wave velocity in an intact pile, error = 5.3%). The reflective signal from the pile tip (RST) arrives at t=tRST=0.00967 s (tRST is the time of arrival of the reflective signal from the pile tip), and the length of the pile H=Vrod×tRST/2=18.32 m (error = 3.6%). When the degree of defect is relatively great (WRN=0.3 m) (namely, the peak of RS1 is above zero axis), the second reflective signal from the defect (RS2) appears and arrives at t≈2tRS1  (t=0.00598 s), and the amplitude of RS2 is much lower than that of RS1 because of quick dissipation of the energy. It is worth noting that this signal (RS2) can be used to further check the location of the defect. For instance, the embedment of the defect can be calculated as FH=Vrod×tRS2/4=5.67 m (error = 5.5%). It is also observed from Figure 4 that the amplitude of RS1 increases greatly with the decrease of WRN. This result indicates that the higher the amplitude of the reflective signal is, the greater the degree of defect will be. Furthermore, for the pile with a neck, the phase of RS1 is the same as that of the incident pulse due to the abrupt decrease of the cross-sectional area.Figure 4
Variation of the normalized velocity response of the pile top with the neck width (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure5 shows the influence of the bulb width (WRB) on the velocity response of the pile top. It can be seen from Figure 5 that RS1 and RST arrive at t=tRS1=0.00296 s and t=tRST=0.00973 s, respectively. Therefore, FH=Vrod×tRS1/2=5.61 m (error = 6.5%) and H=Vrod×tRST/2=18.32 m (error = 3.6%). It is also observed from Figure 6 that the amplitude of the RS1 increases greatly with the increase of WRB. Moreover, the phase of RS1 is the opposite to that of the incident pulse due to the abrupt increase of the cross-sectional area. However, the phase of RS2 is the same as that of the incident pulse, which is similar to the reflective signal of the pile with a neck. As a result, RS2 is easily identified as a defect by mistake.Figure 5
Variation of the normalized velocity response of the pile top with the bulb width (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure 6
Variation of the normalized velocity response of the pile top with the longitudinal wave velocity in weak concrete (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure6 shows the influence of the longitudinal wave velocity in a weak concrete (Vp2) on the velocity response of the pile top. It can be seen from Figure 6 that the velocity response curve of the pile with a weak concrete is similar to that of the pile with a neck. The phase of RS1 is the same as that of the incident pulse due to the abrupt decrease of Vp2 and the amplitude of RS1 increases greatly with the decrease of Vp2. Furthermore, the time of arrival of RST increases obviously with the decrease of Vp2. The reason for this result is that RST will need much more time to arrive at the pile top as Vp2 decreases.


#### 4.1.2. Effect of Material Damping of the Pile

Figure7 shows the influence of material damping of the pile on the velocity response of the pile top. It can be seen from Figure 7 that the material damping of the pile has marked influence on the velocity response of the pile top. The amplitude of RS1 and RST decreases markedly with the increasing material damping of the pile, which results in the fact that the degree of defect seems to be smaller than the actual degree. Moreover, the width of RS1 and RST shows obvious increase with the increase of pile material damping, resulting in the marginal of the reflective signal becoming more and more ambiguous. As a result, it is difficult to accurately obtain the time of arrival of the signal reflected from the defect and determine the location of the defect. For instance, RS1 arrives at t=tRS1=0.00311 s and FH=Vrod×tRS1/2=5.89 m (error = 1.8%), when pile material damping is not considered (η¯j=0). However, when pile material damping is considered (i.e., η¯j=0.006), RS1 arrives at t=tRS1=0.00265 s, and FH=Vrod×tRS1/2=5.02 m, the error can be reached to 16.3%. Therefore, it is needed to increase the calculated embedment depth of the defect when determining the location of the defect in the field test.Figure 7
Variation of the normalized velocity response of the pile top with material damping of the pile (N=3, WRN=0.30 m, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).


#### 4.1.3. Effect of Length of Defect

Figure8 shows the influence of the neck length on the velocity response of the pile top. It can be seen from Figure 8 that the neck length has marked influence on the velocity response. As the neck length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase. Moreover, the time of arrival of RS2 and RST is difficult to accurately identify with the increase of neck length.Figure 8
Variation of the normalized velocity response of the pile top with neck length (N=3, WRN=0.30 m, FH=h1=6.0 m).Figure9 shows the influence of the bulb length on the velocity response of the pile top. It can be seen from Figure 9 that the bulb length has marked influence on the velocity response of the pile top. When the bulb length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase, and the time of arrival of RS2 is difficult to obtain. However, the increase of bulb length has negligible influence on the amplitude and the time of arrival of RST.Figure 9
Variation of the normalized velocity response of the pile top with blub length (N=3, WRB=0.70 m, FH=h1=6.0 m).Figure10 shows the influence of the bulb length on the velocity response of the pile top. It can be seen from Figure 10 that the weak concrete length has marked influence on the velocity response of the pile top. As the weak concrete length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase. Furthermore, the time of arrival of RST increases with the increase of weak concrete length.Figure 10
Variation of the normalized velocity response of the pile top with weak concrete length (N=3, VP2=2000 m/s, FH=h1=6.0 m).


#### 4.1.4. Effect of Surrounding Soil Properties

In some cases, the soil properties may change greatly in certain embedment depth. In order to explicitly identify the influence of the variation of the surrounding soil properties, the pile is assumed to be homogeneous and four cases of soil properties are investigated. Case 1: the shear velocity of soil layers decreases along the depth direction (namely,V1=60 m/s, V2=100 m/s, V3=140 m/s); Case 2: the shear velocity of soil layers increases along the depth direction (namely, V1=140 m/s, V2=100 m/s, V3=60 m/s); Case 3: the soil has a soft interlayer (namely, V1=120 m/s, V2=60 m/s, V3=140 m/s); Case 4: the soil has a hard interlayer (namely, V1=120 m/s, V2=240 m/s, V3=140 m/s). It can be seen from Figure 11(a) that the velocity curve shows minor upward inclination, as the shear velocity of the soil gradually decreases along the depth direction. However, as the shear velocity of soil gradually increases along the depth direction, the velocity curve shows minor downward inclination. Based on these results, it can be concluded that the gradual variation of the surrounding soil properties along the depth has negligible influence on the identification of the defect in a pile. It can be seen from Figure 11(b) that the abrupt variation of the surrounding soil properties along the depth has certain influence on the velocity response of the pile. The phase of the reflective signal of the interface of adjacent soil layers is, respectively, the same as and opposite to that of the incident pulse for Case 3 and Case 4. The shape of the reflective signal of the soil with a soft or hard interlayer is in a half-sine form, which is obviously different from that of the defective pile. Therefore, this result should be noted in the field test in case of identifying this kind of signal as a defect by mistake.Figure 11
Effect of interface of adjacent soil layers on the normalized velocity response of the pile (N=3, h1=12.0 m, h2=1.0 m, h3=6.0 m).
(a)
(b)


### 4.2. Pile Containing Multidefects

In Section4.1, the mechanism of longitudinal wave propagating in a pile with a single defect has received detailed investigation. However, in practical engineering, a pile may contain more than one flaw along the pile body. It is noted that the reflected signals from different defects in the velocity curve may interact due to the existence of the potential second reflected signal from the defect. Therefore, the following section will take a pile with double defects as an example to illustrate the influence of the signal reflected from one defect near the pile top on the signal reflected from the other defect near the pile tip.Figure12 shows the time history of the velocity response of the pile with a bulb near the pile top and a neck near the pile tip. It can be seen from Figure 12 that the characteristics of the reflected signals from bulb and neck are obvious for r4=0.60 m when considering the relative amplitude of the signals. However, the time of arrival of the reflected signal from the neck (t≈0.00300 s) is very close to that of the second reflected signal from the bulb (t≈0.00591 s) for r4=0.70 m; accordingly the neck near the pile tip is easy to be identified as the second reflected signal from the bulb (referred to in Figure 9). The reason for this result is that the second reflected signal from the bulb and reflected signal from the neck overlap in the time domain. Therefore, if the amplitude and width of the first reflected signal from the bulb are great, and the second reflected signal is obvious, other testing methods are suggested to further check the integrity of the pile.Figure 12
Time history of the normalized velocity response of the pile with bulb near the pile top and neck near the pile tip (N=5, h1=7.0 m, h2=0.5 m, h3=4.5 m, h4=1.0 m, h5=6.0 m).Figure13 shows the time history of the velocity response of the pile with a neck near the pile top and a bulb near the pile tip. It can be seen from Figure 13 that the characteristics of the reflected signals from the neck and bulb are obvious for r4=0.40 m. However, the reflected signal from the neck cannot be identified for r4=0.30 m. The reason for this result is that the reflected signal from the bulb is overlapped and covered up by the second reflected signal from the neck. It is worth noting that this case will not influence the result of pile integrity test because of the existence of the neck near the pile top.Figure 13
Time history of the normalized velocity response of the pile with neck near the pile top and bulb near the pile tip (N=5, h1=7.0 m, h2=0.5 m, h3=4.5 m, h4=1.0 m, h5=6.0 m).


## 4.1. Pile Containing a Single Defect

### 4.1.1. Effect of Width of Defect

Figure4 shows the influence of the neck width (WRN) on the velocity response of the pile top. It can be seen from Figure 4 that the first reflective signal from the defect (RS1) arrives at t=tRS1=0.00300 s (tRS1 is the time of arrival of the first reflective signal from the defect). Accordingly, the embedment depth of the defect can be calculated as FH=Vrod×tRS1/2=5.68 m (Vrod is the 1D elastic longitudinal wave velocity in an intact pile, error = 5.3%). The reflective signal from the pile tip (RST) arrives at t=tRST=0.00967 s (tRST is the time of arrival of the reflective signal from the pile tip), and the length of the pile H=Vrod×tRST/2=18.32 m (error = 3.6%). When the degree of defect is relatively great (WRN=0.3 m) (namely, the peak of RS1 is above zero axis), the second reflective signal from the defect (RS2) appears and arrives at t≈2tRS1  (t=0.00598 s), and the amplitude of RS2 is much lower than that of RS1 because of quick dissipation of the energy. It is worth noting that this signal (RS2) can be used to further check the location of the defect. For instance, the embedment of the defect can be calculated as FH=Vrod×tRS2/4=5.67 m (error = 5.5%). It is also observed from Figure 4 that the amplitude of RS1 increases greatly with the decrease of WRN. This result indicates that the higher the amplitude of the reflective signal is, the greater the degree of defect will be. Furthermore, for the pile with a neck, the phase of RS1 is the same as that of the incident pulse due to the abrupt decrease of the cross-sectional area.Figure 4
Variation of the normalized velocity response of the pile top with the neck width (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure5 shows the influence of the bulb width (WRB) on the velocity response of the pile top. It can be seen from Figure 5 that RS1 and RST arrive at t=tRS1=0.00296 s and t=tRST=0.00973 s, respectively. Therefore, FH=Vrod×tRS1/2=5.61 m (error = 6.5%) and H=Vrod×tRST/2=18.32 m (error = 3.6%). It is also observed from Figure 6 that the amplitude of the RS1 increases greatly with the increase of WRB. Moreover, the phase of RS1 is the opposite to that of the incident pulse due to the abrupt increase of the cross-sectional area. However, the phase of RS2 is the same as that of the incident pulse, which is similar to the reflective signal of the pile with a neck. As a result, RS2 is easily identified as a defect by mistake.Figure 5
Variation of the normalized velocity response of the pile top with the bulb width (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure 6
Variation of the normalized velocity response of the pile top with the longitudinal wave velocity in weak concrete (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure6 shows the influence of the longitudinal wave velocity in a weak concrete (Vp2) on the velocity response of the pile top. It can be seen from Figure 6 that the velocity response curve of the pile with a weak concrete is similar to that of the pile with a neck. The phase of RS1 is the same as that of the incident pulse due to the abrupt decrease of Vp2 and the amplitude of RS1 increases greatly with the decrease of Vp2. Furthermore, the time of arrival of RST increases obviously with the decrease of Vp2. The reason for this result is that RST will need much more time to arrive at the pile top as Vp2 decreases.


### 4.1.2. Effect of Material Damping of the Pile

Figure7 shows the influence of material damping of the pile on the velocity response of the pile top. It can be seen from Figure 7 that the material damping of the pile has marked influence on the velocity response of the pile top. The amplitude of RS1 and RST decreases markedly with the increasing material damping of the pile, which results in the fact that the degree of defect seems to be smaller than the actual degree. Moreover, the width of RS1 and RST shows obvious increase with the increase of pile material damping, resulting in the marginal of the reflective signal becoming more and more ambiguous. As a result, it is difficult to accurately obtain the time of arrival of the signal reflected from the defect and determine the location of the defect. For instance, RS1 arrives at t=tRS1=0.00311 s and FH=Vrod×tRS1/2=5.89 m (error = 1.8%), when pile material damping is not considered (η¯j=0). However, when pile material damping is considered (i.e., η¯j=0.006), RS1 arrives at t=tRS1=0.00265 s, and FH=Vrod×tRS1/2=5.02 m, the error can be reached to 16.3%. Therefore, it is needed to increase the calculated embedment depth of the defect when determining the location of the defect in the field test.Figure 7
Variation of the normalized velocity response of the pile top with material damping of the pile (N=3, WRN=0.30 m, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).


### 4.1.3. Effect of Length of Defect

Figure8 shows the influence of the neck length on the velocity response of the pile top. It can be seen from Figure 8 that the neck length has marked influence on the velocity response. As the neck length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase. Moreover, the time of arrival of RS2 and RST is difficult to accurately identify with the increase of neck length.Figure 8
Variation of the normalized velocity response of the pile top with neck length (N=3, WRN=0.30 m, FH=h1=6.0 m).Figure9 shows the influence of the bulb length on the velocity response of the pile top. It can be seen from Figure 9 that the bulb length has marked influence on the velocity response of the pile top. When the bulb length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase, and the time of arrival of RS2 is difficult to obtain. However, the increase of bulb length has negligible influence on the amplitude and the time of arrival of RST.Figure 9
Variation of the normalized velocity response of the pile top with blub length (N=3, WRB=0.70 m, FH=h1=6.0 m).Figure10 shows the influence of the bulb length on the velocity response of the pile top. It can be seen from Figure 10 that the weak concrete length has marked influence on the velocity response of the pile top. As the weak concrete length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase. Furthermore, the time of arrival of RST increases with the increase of weak concrete length.Figure 10
Variation of the normalized velocity response of the pile top with weak concrete length (N=3, VP2=2000 m/s, FH=h1=6.0 m).


### 4.1.4. Effect of Surrounding Soil Properties

In some cases, the soil properties may change greatly in certain embedment depth. In order to explicitly identify the influence of the variation of the surrounding soil properties, the pile is assumed to be homogeneous and four cases of soil properties are investigated. Case 1: the shear velocity of soil layers decreases along the depth direction (namely,V1=60 m/s, V2=100 m/s, V3=140 m/s); Case 2: the shear velocity of soil layers increases along the depth direction (namely, V1=140 m/s, V2=100 m/s, V3=60 m/s); Case 3: the soil has a soft interlayer (namely, V1=120 m/s, V2=60 m/s, V3=140 m/s); Case 4: the soil has a hard interlayer (namely, V1=120 m/s, V2=240 m/s, V3=140 m/s). It can be seen from Figure 11(a) that the velocity curve shows minor upward inclination, as the shear velocity of the soil gradually decreases along the depth direction. However, as the shear velocity of soil gradually increases along the depth direction, the velocity curve shows minor downward inclination. Based on these results, it can be concluded that the gradual variation of the surrounding soil properties along the depth has negligible influence on the identification of the defect in a pile. It can be seen from Figure 11(b) that the abrupt variation of the surrounding soil properties along the depth has certain influence on the velocity response of the pile. The phase of the reflective signal of the interface of adjacent soil layers is, respectively, the same as and opposite to that of the incident pulse for Case 3 and Case 4. The shape of the reflective signal of the soil with a soft or hard interlayer is in a half-sine form, which is obviously different from that of the defective pile. Therefore, this result should be noted in the field test in case of identifying this kind of signal as a defect by mistake.Figure 11
Effect of interface of adjacent soil layers on the normalized velocity response of the pile (N=3, h1=12.0 m, h2=1.0 m, h3=6.0 m).
(a)
(b)


## 4.1.1. Effect of Width of Defect

Figure4 shows the influence of the neck width (WRN) on the velocity response of the pile top. It can be seen from Figure 4 that the first reflective signal from the defect (RS1) arrives at t=tRS1=0.00300 s (tRS1 is the time of arrival of the first reflective signal from the defect). Accordingly, the embedment depth of the defect can be calculated as FH=Vrod×tRS1/2=5.68 m (Vrod is the 1D elastic longitudinal wave velocity in an intact pile, error = 5.3%). The reflective signal from the pile tip (RST) arrives at t=tRST=0.00967 s (tRST is the time of arrival of the reflective signal from the pile tip), and the length of the pile H=Vrod×tRST/2=18.32 m (error = 3.6%). When the degree of defect is relatively great (WRN=0.3 m) (namely, the peak of RS1 is above zero axis), the second reflective signal from the defect (RS2) appears and arrives at t≈2tRS1  (t=0.00598 s), and the amplitude of RS2 is much lower than that of RS1 because of quick dissipation of the energy. It is worth noting that this signal (RS2) can be used to further check the location of the defect. For instance, the embedment of the defect can be calculated as FH=Vrod×tRS2/4=5.67 m (error = 5.5%). It is also observed from Figure 4 that the amplitude of RS1 increases greatly with the decrease of WRN. This result indicates that the higher the amplitude of the reflective signal is, the greater the degree of defect will be. Furthermore, for the pile with a neck, the phase of RS1 is the same as that of the incident pulse due to the abrupt decrease of the cross-sectional area.Figure 4
Variation of the normalized velocity response of the pile top with the neck width (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure5 shows the influence of the bulb width (WRB) on the velocity response of the pile top. It can be seen from Figure 5 that RS1 and RST arrive at t=tRS1=0.00296 s and t=tRST=0.00973 s, respectively. Therefore, FH=Vrod×tRS1/2=5.61 m (error = 6.5%) and H=Vrod×tRST/2=18.32 m (error = 3.6%). It is also observed from Figure 6 that the amplitude of the RS1 increases greatly with the increase of WRB. Moreover, the phase of RS1 is the opposite to that of the incident pulse due to the abrupt increase of the cross-sectional area. However, the phase of RS2 is the same as that of the incident pulse, which is similar to the reflective signal of the pile with a neck. As a result, RS2 is easily identified as a defect by mistake.Figure 5
Variation of the normalized velocity response of the pile top with the bulb width (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure 6
Variation of the normalized velocity response of the pile top with the longitudinal wave velocity in weak concrete (N=3, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).Figure6 shows the influence of the longitudinal wave velocity in a weak concrete (Vp2) on the velocity response of the pile top. It can be seen from Figure 6 that the velocity response curve of the pile with a weak concrete is similar to that of the pile with a neck. The phase of RS1 is the same as that of the incident pulse due to the abrupt decrease of Vp2 and the amplitude of RS1 increases greatly with the decrease of Vp2. Furthermore, the time of arrival of RST increases obviously with the decrease of Vp2. The reason for this result is that RST will need much more time to arrive at the pile top as Vp2 decreases.


## 4.1.2. Effect of Material Damping of the Pile

Figure7 shows the influence of material damping of the pile on the velocity response of the pile top. It can be seen from Figure 7 that the material damping of the pile has marked influence on the velocity response of the pile top. The amplitude of RS1 and RST decreases markedly with the increasing material damping of the pile, which results in the fact that the degree of defect seems to be smaller than the actual degree. Moreover, the width of RS1 and RST shows obvious increase with the increase of pile material damping, resulting in the marginal of the reflective signal becoming more and more ambiguous. As a result, it is difficult to accurately obtain the time of arrival of the signal reflected from the defect and determine the location of the defect. For instance, RS1 arrives at t=tRS1=0.00311 s and FH=Vrod×tRS1/2=5.89 m (error = 1.8%), when pile material damping is not considered (η¯j=0). However, when pile material damping is considered (i.e., η¯j=0.006), RS1 arrives at t=tRS1=0.00265 s, and FH=Vrod×tRS1/2=5.02 m, the error can be reached to 16.3%. Therefore, it is needed to increase the calculated embedment depth of the defect when determining the location of the defect in the field test.Figure 7
Variation of the normalized velocity response of the pile top with material damping of the pile (N=3, WRN=0.30 m, h1=12.5 m, FL=h2=0.5 m, FH=h1=6.0 m).


## 4.1.3. Effect of Length of Defect

Figure8 shows the influence of the neck length on the velocity response of the pile top. It can be seen from Figure 8 that the neck length has marked influence on the velocity response. As the neck length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase. Moreover, the time of arrival of RS2 and RST is difficult to accurately identify with the increase of neck length.Figure 8
Variation of the normalized velocity response of the pile top with neck length (N=3, WRN=0.30 m, FH=h1=6.0 m).Figure9 shows the influence of the bulb length on the velocity response of the pile top. It can be seen from Figure 9 that the bulb length has marked influence on the velocity response of the pile top. When the bulb length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase, and the time of arrival of RS2 is difficult to obtain. However, the increase of bulb length has negligible influence on the amplitude and the time of arrival of RST.Figure 9
Variation of the normalized velocity response of the pile top with blub length (N=3, WRB=0.70 m, FH=h1=6.0 m).Figure10 shows the influence of the bulb length on the velocity response of the pile top. It can be seen from Figure 10 that the weak concrete length has marked influence on the velocity response of the pile top. As the weak concrete length increases from 0.5 m to 2.0 m, the width and amplitude of RS1 and RS2 show marked increase. Furthermore, the time of arrival of RST increases with the increase of weak concrete length.Figure 10
Variation of the normalized velocity response of the pile top with weak concrete length (N=3, VP2=2000 m/s, FH=h1=6.0 m).


## 4.1.4. Effect of Surrounding Soil Properties

In some cases, the soil properties may change greatly in certain embedment depth. In order to explicitly identify the influence of the variation of the surrounding soil properties, the pile is assumed to be homogeneous and four cases of soil properties are investigated. Case 1: the shear velocity of soil layers decreases along the depth direction (namely,V1=60 m/s, V2=100 m/s, V3=140 m/s); Case 2: the shear velocity of soil layers increases along the depth direction (namely, V1=140 m/s, V2=100 m/s, V3=60 m/s); Case 3: the soil has a soft interlayer (namely, V1=120 m/s, V2=60 m/s, V3=140 m/s); Case 4: the soil has a hard interlayer (namely, V1=120 m/s, V2=240 m/s, V3=140 m/s). It can be seen from Figure 11(a) that the velocity curve shows minor upward inclination, as the shear velocity of the soil gradually decreases along the depth direction. However, as the shear velocity of soil gradually increases along the depth direction, the velocity curve shows minor downward inclination. Based on these results, it can be concluded that the gradual variation of the surrounding soil properties along the depth has negligible influence on the identification of the defect in a pile. It can be seen from Figure 11(b) that the abrupt variation of the surrounding soil properties along the depth has certain influence on the velocity response of the pile. The phase of the reflective signal of the interface of adjacent soil layers is, respectively, the same as and opposite to that of the incident pulse for Case 3 and Case 4. The shape of the reflective signal of the soil with a soft or hard interlayer is in a half-sine form, which is obviously different from that of the defective pile. Therefore, this result should be noted in the field test in case of identifying this kind of signal as a defect by mistake.Figure 11
Effect of interface of adjacent soil layers on the normalized velocity response of the pile (N=3, h1=12.0 m, h2=1.0 m, h3=6.0 m).
(a)
(b)


## 4.2. Pile Containing Multidefects

In Section4.1, the mechanism of longitudinal wave propagating in a pile with a single defect has received detailed investigation. However, in practical engineering, a pile may contain more than one flaw along the pile body. It is noted that the reflected signals from different defects in the velocity curve may interact due to the existence of the potential second reflected signal from the defect. Therefore, the following section will take a pile with double defects as an example to illustrate the influence of the signal reflected from one defect near the pile top on the signal reflected from the other defect near the pile tip.Figure12 shows the time history of the velocity response of the pile with a bulb near the pile top and a neck near the pile tip. It can be seen from Figure 12 that the characteristics of the reflected signals from bulb and neck are obvious for r4=0.60 m when considering the relative amplitude of the signals. However, the time of arrival of the reflected signal from the neck (t≈0.00300 s) is very close to that of the second reflected signal from the bulb (t≈0.00591 s) for r4=0.70 m; accordingly the neck near the pile tip is easy to be identified as the second reflected signal from the bulb (referred to in Figure 9). The reason for this result is that the second reflected signal from the bulb and reflected signal from the neck overlap in the time domain. Therefore, if the amplitude and width of the first reflected signal from the bulb are great, and the second reflected signal is obvious, other testing methods are suggested to further check the integrity of the pile.Figure 12
Time history of the normalized velocity response of the pile with bulb near the pile top and neck near the pile tip (N=5, h1=7.0 m, h2=0.5 m, h3=4.5 m, h4=1.0 m, h5=6.0 m).Figure13 shows the time history of the velocity response of the pile with a neck near the pile top and a bulb near the pile tip. It can be seen from Figure 13 that the characteristics of the reflected signals from the neck and bulb are obvious for r4=0.40 m. However, the reflected signal from the neck cannot be identified for r4=0.30 m. The reason for this result is that the reflected signal from the bulb is overlapped and covered up by the second reflected signal from the neck. It is worth noting that this case will not influence the result of pile integrity test because of the existence of the neck near the pile top.Figure 13
Time history of the normalized velocity response of the pile with neck near the pile top and bulb near the pile tip (N=5, h1=7.0 m, h2=0.5 m, h3=4.5 m, h4=1.0 m, h5=6.0 m).


## 5. Application in Engineering

Figure14 shows the comparison of the fitted theoretical curves and the measured curve from the model pile. The pile is an embedded, prefabricated, nonuniform square pile with a bulb near the pile top and a neck near the pile tip, which was installed at a site in China. The length and cross-sectional dimension of the pile are 10 m and 400 mm × 400 mm, respectively. Field measurements recorded the propagation time of the impulse traveling from the pile top to pile tip to pile top as 4.6 ms. 1D elastic longitudinal wave velocity of the pile is estimated to be 4347.8 m/s. The surrounding soil is soft plastic viscoelastic backfill soil with mass density of 1700 kg/m3 and shear wave velocity of 90 m/s. For the purpose of analysis, the square pile is considered as an equivalent circular pile with radius r=225.5 mm. The theoretical velocity curve has been derived through adjusting the radius, depths, and material damping of the pile segments as well as the duration of impulse to give a good fit to the measured curve. The fitted curves require the following input: (1) the embedment of the bulb is 4.9 m, and the cross-sectional area of the bulb to the normal shaft area is 210%; (2) the embedment of the neck is 7.20 m, and the cross-sectional area of the neck decreased to a ratio of 54% of the normal shaft area; (3) the impulse width of the exciting force is 0.59 ms, and the dimensionless damping coefficient of the pile η¯pj=0.0015. By comparison, the actual condition of the pile is listed as follows: at the depth of 4.8 m, the enlarged area ratio is 200%; at the depth of 7.16 m, the cross-sectional area decreased to a ratio of 55%. Moreover, it is noted from Figure 14 that fitted curves 1 and 2 match well with the measured curve when t<0.0045 s. When t>0.0045 s, in the case of the damping coefficient at the pile toe cpb taken in terms of (14), the amplitude of RST in fitted curve 1 is much higher than that in measured curve. However, RST in fitted curve 2 matches well with that in the measured curve, when cpb is taken as a relatively great value. The reason for this result is that the amplitude of RST decreases with the increase of soil damping at the pile toe and the corresponding soil damping may be underestimated by virtue of (14) for the square pile.Figure 14
Comparison of the fitted theoretical curves and measured curve from the model pile.


## 6. Conclusions

(1)
For a pile with a single neck, the phase of RS1 and RS2 is the same as that of the incident pulse. The amplitude of the signals reflected from the neck increases markedly with the decrease of the neck width.(2)
For a pile with a single bulb, the phase of RS1 is the opposite to that of the incident pulse. However, the phase of RS2 is the same as that of the incident pulse, which is easily identified as a defect by mistake. The amplitude of the signals reflected from the bulb increases greatly with the increase of bulb width.(3)
For a pile with a single weak concrete, the phase of RS1 and RS2 is the same as that of the incident pulse due to the abrupt decrease ofVp2 and the amplitude of RS1 increases greatly with the decrease of Vp2. The time of arrival of RST increases obviously with the decrease of Vp2.(4)
The amplitude and width of RS1 and RST, respectively, decrease and increase with the increase of pile material damping, which makes the degree of the defect seem to be smaller than the actual degree and the time of arrival of the signal reflected from the defect difficult to accurately identify.(5)
As the length of single defect (such as neck, bulb, or weak concrete) increases, the width and amplitude of RS1 and RS2 from the defect show marked increase. For a pile containing a neck, it is noted that the time of arrival of RST is difficult to identify with the increase of neck length.(6)
The gradual variation of the surrounding soil properties along the depth has negligible influence on identification of the defect in a pile. The abrupt change in adjacent soil layer properties has certain influence on the velocity response and this result should be noted in the field test in case of identifying this kind of signal as a defect by mistake.(7)
For a pile with a bulb near the pile top and a neck near the pile tip, when the amplitude and width of the first reflected signal from the bulb are great, and the second reflected signal is obvious, other testing methods are suggested to further check the pile integrity.



---

*Source: 495253-2015-06-08.xml*