MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2017/1794950 1794950 Research Article Vertical Vibration Characteristics of a Variable Impedance Pile Embedded in Layered Soil http://orcid.org/0000-0001-5473-1560 Wu Wenbing 1 2 Xu Xuelian 1 Liu Hao 1 http://orcid.org/0000-0001-8485-6364 Fang Changliang 1 3 http://orcid.org/0000-0002-2067-5350 Dou Bin 1 Liang Rongzhu 1 Quinn Dane 1 Engineering Research Centre of Rock-Soil Drilling & Excavation and Protection Ministry of Education China University of Geosciences Wuhan Hubei 430074 China cugb.edu.cn 2 College of Civil Engineering and Architecture Guangxi University Nanning Guangxi 530004 China gxu.edu.cn 3 Institute of Drilling Engineering and Fluid Mining TU Bergakademie Freiberg Freiberg Germany tu-freiberg.de 2017 342017 2017 30 09 2016 21 02 2017 342017 2017 Copyright © 2017 Wenbing Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In engineering applications, various defects such as bulging, necking, slurry crappy, and weak concrete are always observed during pile integrity testing. To provide more reasonable basis for assessing the above defects, this paper proposed simple and computationally efficient solutions to investigate the vertical vibration characteristics of a variable impedance pile embedded in layered soil. The governing equations of pile-soil system undergoing a vertical dynamic loading are built based on the plane strain model and fictitious soil pile model. By employing the Laplace transform method and impedance function transfer method, the analytical solution of the velocity response at the pile head is derived in the frequency domain. Then, the corresponding semianalytical solution in the time domain for the velocity response of a pile subjected to a semisinusoidal force applied at the pile head is obtained by adopting inverse Fourier transform and convolution theorem. Based on the presented solutions, a parametric study is conducted to study the vertical vibration characteristics of variable cross-section pile and variable modulus pile. The study gives an important insight into the evaluation of the construction quality of pile.

National Natural Science Foundation of China 51309207 51322807 51678547 China Postdoctoral Science Foundation 2016M600711 MOE Engineering Research Center of Rock-Soil Drilling & Excavation and Protection 201402 Central Universities-Cradle CUGL150411
1. Introduction

With the wide application of piles for supporting various upper structures all over the world, the construction quality of pile foundations has attracted increasing attention in recent decades. In practice, whatever the construction method is, various defects in pile are often found during pile integrity testing, such as bulging, necking, slurry scrappy, and weak concrete. The common methods for assessing the pile quality are static load testing and dynamic testing. Static load testing is inconvenient and time-consuming for a great deal of surcharge load is required to test the ultimate bearing capacity of pile foundation. Dynamic testing, including low-strain and high-strain testing, has been widely utilized in the nondestructive examination of pile due to its simplicity and reliability .

Since the pile dynamics is the basis of various methods of dynamic pile testing, many investigators have paid their attention to develop various dynamic interaction models for the dynamic response of pile subjected to vertical loads. Van Koten et al.  primarily utilized a distributed system of single Voigt unit (a linear spring and a dashpot connected in parallel) to simulate the soil reactions acting on the pile and derived an analytical solution for the dynamic interaction of pile-soil system. Following Van Koten et al.’s research, a large number of researchers applied this method to investigate the vertical vibration of pile foundation for different engineering conditions (e.g., Randolph and Deeks ; Michaelides et al. ; Yesilce and Catal ; Wang et al. ; Zhang et al. [10, 11]). Although the Voigt model is of great simplicity and high efficiency in practical engineering, it cannot allow for the wave effect of stress propagating in pile surrounding soil. Allowing more rigorous treatment of pile surrounding soil, Novak  developed the plane strain model assuming an elastic soil layer consisting of independent infinitesimally thin horizontal layers extending to infinity. The plane strain model has attracted numerous followings for lots of engineering cases, such as Rajapakse and Shah , Mamoon and Banerjee , El Naggar and Novak , Militano and Rajapakse , and Wu et al. . It is worth noting that the plane strain model is one of the most important methods to simulate the interaction between pile and surrounding soil due to its compact and elegant properties; that is, there is zero strain in the vertical direction and the stress waves propagate only in the horizontal direction. In a different approach, Nogami and Novak  first adopted a three-dimensional continuum model (i.e., the soil layer is three-dimensional and its wave effect is considered) to investigate the vertical vibration of pile foundation. Afterwards, Nogami and Konagai , Senjuntichai et al. , Hu et al. , Wang et al. , Wu et al. [23, 24], and Zheng et al. [25, 26] also developed various three-dimensional continuum models extensively in their investigations on the dynamic response of pile-soil system. It can be found that the three-dimensional continuum model is theoretically more rigorous than the previously mentioned Voigt models and plane strain models. Nevertheless, it is very difficult to promote the use of these three-dimensional continuum models in layered soil because of its complexity.

Based on the above literature review, the main objective of this study is to develop simple and computationally efficient solutions for the analysis of vertical vibration characteristics of a variable impedance pile embedded in layered soil. In this paper, the dynamic interaction between pile and its surrounding soil is described by the plane strain model, and the dynamic interaction between pile and its underlying soil is simulated by the fictitious soil pile model . Laplace transform and transfer function technique have been utilized to establish the impedance function of pile head in the frequency domain. By virtue of Fourier inverse transform, the semianalytical solution of the velocity response of pile head subjected to a semisinusoidal force is also obtained, this constituting the theoretical basis of the low-strain testing method of pile integrity. A parametric study is conducted to investigate the influence of pile properties on the vertical dynamic response of the variable impedance pile.

2. Mathematical Model and Assumptions 2.1. Computational Model

Considering the defects of pile and the inhomogeneity of soil, the vertical vibration of a variable impedance pile embedded in layered soil is theoretically studied based on the fictitious soil pile model. Figure 1 shows a schematic diagram of the pile-soil system, in which the finite soil layers underlying the pile end are assumed to be soil column, namely, fictitious soil pile whose cross-section area is the same as the pile. The fictitious soil pile is in perfect contact with the pile end and satisfies the plane section of deformation assumption. Taking into account the variation of modulus or cross-sectional dimension of the pile or soil properties, the pile-soil system is divided into a total number of m along the vertical direction, which are numbered by 1,2,,jm from the bottom of fictitious soil pile to pile head. The properties of pile (including fictitious soil pile) and surrounding soil layer are assumed to be homogeneous within each segment but may vary from segment to segment. lj and hj denote the thickness and top surface depth of the jth pile-soil segment, respectively. ρpj and rj represent the density and radius of the pile, respectively. Hp and Hs indicate the length of pile and the thickness of pile end soil, respectively. q(t) defines the vertical exciting force acting on the pile head.

Schematic of pile-soil interaction model.

2.2. Underlying Assumptions

The pile-soil model is developed based on the following assumptions:

The pile-soil system is subjected to small deformations and strains, and the pile and soil remain in perfect contact during the vibration. The stress and displacement of the interface of pile-soil system are continuous, and only vertical displacements are considered.

The pile is vertical, viscoelastic, and cylindrical. The support conditions are based on the fictitious soil pile model that overlies rigid bedrock.

The surrounding soil is isotropic linear viscoelastic medium with a hysteretic-type damping regardless of frequency. The surrounding soil is infinite in the radial direction and its free top surface has no normal or shear stress.

The dynamic stress of soil transfers to the pile shaft through the complex stiffness on the contact surface of the every segment of pile-soil system.

The conditions of displacement continuity and force equilibrium are satisfied at the interface of fictitious soil pile and pile end.

3. Governing Equations and Their Solutions 3.1. Dynamic Equation of Soil and Its Solution

Denoting Wj=Wj(r) as the vertical displacement of the jth soil layer, the dynamic equilibrium equation of the jth soil layer undergoing vertical dynamic loading can be expressed as [13, 15](1)r2d2Wjdr2+rdWjdr-βj2r2Wj=0,where βj=iω/vsj1+iDsj; vsj=Gj/ρj represents the shear wave velocity of the jth soil layer, and ρj, Gj, and Dsj indicate the density, shear modulus, and material damping of the jth soil layer, respectively. ω denotes the circular frequency and i=-1 depicts the imaginary unit.

Based on the theory of Bessel functions, the general solution of (1) is derived as(2)Wjr=AjK0βjr+BjI0βjr,where I0(·) and K0(·) represent the modified Bessel functions of order zero of the first and second kind, respectively. Aj and Bj are integration constants determined by boundary conditions. Based on assumption (3), it can be obtained that Bj=0 from (2).

The vertical shear stress at arbitrary point in the soil layer is written as(3)τjr=GjdWjrdr=-GjAjβjK1βjr,where Gj=Gj(1+iDsj).

Then, the vertical shear complex stiffness around the pile shaft of jth soil layer is obtained as(4)KKj=-2πrjτjrjWjrj=2πrjGjβjK1βjrjK0βjrj,where I1(·) and K1(·) represent the modified Bessel functions of order one of the first and second kind, respectively.

3.2. Dynamic Equation of Pile (Including Fictitious Soil Pile) and Its Solution

Denoting uj(z,t) as the vertical displacement of the jth pile (including fictitious soil pile) segment, the dynamic equation of pile can be expressed as(5)EpjApj2ujz,tz2+Apjδpj3ujz,tz2t-mpj2ujz,tt2-fjz,t=0,where Epj, mpj, δpj, and Apj=πrj2 represent the elastic modulus, mass per unit length, viscous damping coefficient, and cross-section area, respectively. fj(z,t)=KKjuj(z,t) depicts the frictional force of the jth soil layer acting on the surface of the pile shaft. KKj denotes the vertical shear complex stiffness around the jth pile segment.

Boundary condition at the top of pile can be expressed as(6)EpmApmumz,tz+Apmδpm2umz,tztz=0=-qt.

Boundary condition at the bottom of fictitious soil pile can be written as(7)u1z,tz=h=0.

Continuity conditions of the interface of the adjacent pile or fictitious soil pile segments can be founded as(8)ujz,t=uj+1z,t;EpjApjuiz,tz=Epj+1Apj+1upj+1z,tz.

Initial conditions of pile and fictitious soil pile segments can be derived as(9)ujz,tt=0=0ujz,ttt=0=0.

Denoting Uj(z,s)=0uj(z,t)e-stdt as the Laplace transform of uj(z,t) with respect to time and combining with the initial conditions (9), (5) is transformed by applying the Laplace transform (two-sided) as follows:(10)EpjApj2Ujz,sz2+Apjδpjs2Ujz,sz2-mpjs2Ujz,s-KKjUjz,s=0.

Equation (10) can be further reduced as(11)Vpj21+δpjEpjs2Ujz,s2z-s2+1ρpjApjKKjUjz,s=0, where Vpj, ρpj, and Epj denote the elastic longitudinal wave velocity, density, and elastic modulus of the jth pile (including fictitious soil pile) segment, respectively.

The general solution of (11) can be expressed as(12)Ujz,s=Cjcosλ-jzlj+Djsinλ-jzlj, where λ-j=-(s2+KKj/ρpjApj)tj2/(1+(δpj/Epj)s) illustrates dimensionless eigenvalue. Cj and Dj are constants determined by the boundary conditions. tj=lj/Vpj depicts the propagation time of elastic longitudinal wave in the jth pile (including fictitious soil pile) segment. Tc denotes the propagation time of elastic longitudinal wave propagating from the pile head to the pile end.

Combining with (7), the displacement impedance function at the top of the first fictitious soil pile segment (z=h1) can be obtained as(13)Zp1z=h1=-Ep1Ap1U1z,s/z+Ap1δp1sU1z,s/zz=h1U1z,sz=h1=-ρp1Ap1Vp11+δp1/Ep1sλ-1tanλ-1-ϕ1t1,where ϕ1=arctanZp0t1/ρp1Ap1Vp1λ-11+δp1/Ep1s. Zp0 denotes the displacement function at the bottom of the fictitious soil pile.

Then, following the method of recursion typically utilized in the transfer function technique , the displacement impedance function of the top of the jth pile segment or fictitious soil pile segment can be founded as(14)Zpjz=hj=-EpjApjUjz,s/z+ApjδpjsUjz,s/zz=hjUjz,sz=hj=-ρpjApjVpj1+δpj/Epjsλ-jtanλ-j-ϕjtj, where ϕj=arctanZpj-1tj/ρpjApjVpjλ-j1+δpj/Epjs.

Through further recursion, the displacement impedance function at the pile head can be expressed as(15)Zpmz=0=-EpmApmUmz,s/z+ApmδpmsUmz,s/zz=0Umz,sz=0=-ρpmApmVpm1+δpm/Epmsλ-mtanλ-m-ϕmtm,where ϕm=arctanZpm-1tm/ρpmApmVpmλ-m1+δpm/Epms.

Following the definition of Wang et al. , the velocity transit function at the pile head can be expressed as(16)Hv=s×1Zpm=-stmρpmApmVpm1+δpm/Epmsλ-mtanλ-m-ϕm.

Letting s=iω and substituting it into (16), the velocity response in the frequency domain at the pile head can be derived as(17)Hviω=-1ρpmApmVpmiω-t-m1+δpm/Epmiω-λ-mtanλ-m-ϕm.

Nondimensionalizing expression (17) yields the dimensionless velocity admittance of the pile head:(18)Hv=iω-t-m1+δpm/Epmiω-λ-mtanλ-m-ϕm,where ω-=Tcω and t-m=tm/Tc depict the dimensionless frequency and dimensionless propagation time of elastic longitudinal wave in the mth pile segment, respectively.

By virtue of the inverse Fourier transform, the velocity response at the pile head is derived as V(t)=IFT[Q(ω)Hv(iω)] in the time domain, where Q(ω) represents the Fourier transform of q(t) with respect to time. In particular, for the low-strain integrity testing of pile, the vertical excitation can be regarded as a semisinusoidal force as follows:(19)qt=QmaxsinπT0t,t0,T00,tT0, where T0 and Qmax denote the pulse width and maximum amplitude of the excitation, respectively. The velocity response of pile head in the time domain can be derived as(20)Vt=-12QmaxρpmApmVpm-iω-t-m1+δpm/EpmTciω-λ-mtanλ-m-ϕmT-0π2-T-02ω-21+e-iω-T-0eiω-t-dω-.

The dimensionless velocity response of pile head in the time domain is expressed as(21)V=-12-iω-t-m1+δpm/EpmTciω-λ-mtanλ-m-ϕmT-0π2-T-02ω-21+e-iω-T-0eiω-t-dω-, where T-0=T0/Tc, t-=t/Tc, t, and δ-pm=δpm/(EpmTc) denote the dimensionless pulse width, dimensionless time, time, and dimensionless viscous damping coefficient, respectively.

4. Parametric Study and Discussion

In this section, the vertical dynamic response of variable cross-section and variable modulus pile is systematically investigated, which can provide the theoretical basis for various methods of dynamic pile testing. In order to focus on the influence of variable impedance on the vertical dynamic response of pile, the surrounding soil is set to be homogeneous in the following analysis. For practical engineering, it only needs to divide the surrounding soil into finite layers according to the actual soil conditions. Unless otherwise specified, the mass density and shear wave velocity of pile surrounding soil are given as 1800kg/m3 and 180m/s, respectively. The mass density, shear wave velocity, Poisson’s ratio, viscoelastic damping coefficient, and thickness of pile end soil are 2000kg/m3, 220m/s, 0.35, 1000 N·m−3·s, and 3 times pile diameter, respectively.

4.1. Vertical Vibration Characteristics of Variable Cross-Section Pile

According to engineering experience, the variation of cross-section at different locations will result in different vertical vibration characteristics of pile. Therefore, the influence of bulging and necking at different locations along the pile shaft on the vertical vibration characteristics at pile head is investigated in the following analyses. The length, radius, mass density, and longitudinal wave velocity of pile are given as 20m, 0.5m, 2500kg/m3, and 3800m/s, respectively.

Case 1 (variable cross-section segment at the top of pile shaft).

For this engineering condition, the pile shaft is divided into two segments, which are marked by 1 and 2 from pile end to pile top. Segment 1 has a uniform cross-section and satisfies the design requirements with a length of 19m, and segment 2 has a variable cross-section with a length of 1m. Denote the variable cross-section ratio to describe the variable cross-section degree of pile shaft; that is, RA represents the ratio of actual radius to design radius. RA>1, RA<1, and RA=1 represent the bulging segment, necking segment, and uniform segment, respectively.

Figure 2 shows the influence of the variation degree of variable cross-section segment at pile head on the vertical dynamic response of pile head. As illustrated in Figure 2(a), compared with the uniform cross-section pile, the formant amplitude of velocity admittance curves of pile with necking segment at pile head becomes smaller within low frequency range, but the formant amplitude of velocity admittance curves becomes larger within high frequency range. In contrast, compared with the uniform cross-section pile, the formant amplitude of velocity admittance curves of pile with bulging segment at pile head becomes larger within low frequency range, but the formant amplitude of velocity admittance curves becomes smaller within high frequency range. It can also be seen that the larger the necking degree or bulging degree of variable cross-section segment, the more obvious the variation of the formant amplitude of velocity admittance curves of necking pile or bulging pile. From Figure 2(b), compared with the uniform cross-section pile, the peak values of incident wave signals and preliminary orthokinetic reflected signals from pile end in the reflected wave signal curves of pile with necking segment at pile head become smaller, and the larger the necking degree of variable cross-section segment, the larger the reducing degree of peak value. It can also be seen that the signals between the incident wave and the preliminary orthokinetic reflected wave from pile end tend to inclining downward, and the larger the necking degree of variable cross-section segment at pile head, the more obvious the inclining phenomena. By contrast, for the pile with bulging segment at pile head, the existence of bulging segment at pile head will enlarge the peak values of incident wave signals and preliminary orthokinetic reflected signals from pile end in the reflected wave signal curves. The larger the bulging degree of variable cross-section segment, the larger the increasing range of peak value. It is also observed that the signals between the incident wave and the preliminary orthokinetic reflected wave from pile end tend to inclining upward, and the larger the bulging degree of variable cross-section segment at pile head, the more obvious the inclining phenomena.

Influence of the variation degree of variable cross-section segment on the vertical dynamic response of pile head for Case 1.

Reflected wave signal curves

Case 2 (variable cross-section segment at the end of pile shaft).

Owing to the variation of cross-section, the pile shaft is also divided into two segments, which are marked by 1 and 2 from pile end to pile top. Segment 1 has a variable cross-section with a length of 1m, and segment 2 has a uniform cross-section and satisfies the design requirements with a length of 19m. The other parameters are the same as those shown in Case 1.

Figure 3 illustrates the influence of the variation degree of variable cross-section segment at pile end on the vertical dynamic response of pile head. As shown in Figure 3(a), compared with the uniform cross-section pile, the resonant frequency of the velocity admittance curves of pile with necking segment at pile end gradually increases with the increase of necking degree of variable cross-section segment. The formant amplitude of velocity admittance curves remains basically invariable within low frequency range, but the formant amplitude of velocity admittance curves within high frequency range gradually decreases as the necking degree of variable cross-section segment increases. For the case of pile with bulging segment at pile end, the resonant frequency of the velocity admittance curves gradually decreases as the bulging degree of variable cross-section segment increases. The formant amplitude of velocity admittance curves within high frequency range gradually increases with the increase of the necking degree of variable cross-section segment, but the formant amplitude of velocity admittance curves within low frequency range remains basically invariable. As seen in Figure 3(b), compared with the uniform cross-section pile, the peak values of preliminary orthokinetic reflected signals from pile end in the reflected wave signal curves of pile with necking segment at pile end become smaller, and the arrival time of reflected signals from pile end will be ahead. For the case of pile with bulging segment at pile end, reverse reflected signals with small amplitude will occur before the preliminary orthokinetic reflected signals from pile end. The peak values of preliminary orthokinetic reflected signals from pile end become larger, and the arrival time of reflected signals from pile end will be delayed. It can also be observed that the larger the necking degree or bulging degree of variable cross-section segment, the more obvious the above phenomena.

Influence of the variation degree of variable cross-section segment on the vertical dynamic response of pile head for Case 2.

Reflected wave signal curves

Case 3 (variable cross-section segment at the middle part of pile shaft).

For Case 3, the pile shaft is divided into three segments, which are marked by 1, 2, and 3 from pile end to pile top. Segments 1 and 3 have a uniform cross-section and satisfy the design requirements with a length of 9.5m, and segment 2 has a variable cross-section with a length of 1m. The other parameters are the same as those given in Case 1.

Figure 4 depicts the influence of the variation degree of variable cross-section segment at the middle part of pile shaft on the vertical dynamic response of pile head. From Figure 4(a), compared with the uniform cross-section pile, the oscillatory part of the velocity admittance curves displays increasing complex characteristics of nonuniform peaks with irregular distribution for the pile with necking or bulging segment at the middle part. Therefore, it is difficult to intuitively assess the pile quality from the velocity admittance curves for these conditions. As illustrated in Figure 4(b), for the pile with necking segment at the middle part of pile shaft, the orthokinetic reflected signals at the middle part of pile shaft are first received between the incident wave signals and preliminary orthokinetic reflected signals from pile end, and the reverse reflected signals will occur after the orthokinetic reflected signals. The peak values of preliminary orthokinetic reflected signals from pile end are reduced because of the orthokinetic reflected signals at the middle part of pile shaft. By contrast, for the pile with bulging segment at the middle part of pile shaft, the reverse reflected signals with orthokinetic reflected signals at the middle part of pile shaft are first received between the incident wave signals and preliminary orthokinetic reflected signals from pile end, and the peak values of preliminary orthokinetic reflected signals from pile end are also reduced owing to the orthokinetic reflected signals at the middle part of pile shaft. It can be seen that the larger the necking degree or bulging degree of variable cross-section segment, the more obvious the above phenomena.

Influence of the variation degree of variable cross-section segment on the vertical dynamic response of pile head for Case 3.

Reflected wave signal curves

Case 4 (variable cross-section segment located between the middle and the top of pile shaft).

Similar to Case 3, the pile shaft is also divided into three segments, which are marked by 1, 2, and 3 from pile end to pile top. Segment 1 and segment 3 have a uniform cross-section and satisfy the design requirements, and segment 2 has a variable cross-section. The lengths of each segment are l1=13 m, l2=1 m, and l3=6 m, respectively. The other parameters are the same as those shown in Case 1.

Figure 5 shows the influence of variation degree of variable cross-section segment located between the middle and the top of pile shaft on the vertical dynamic response of pile head. Similar to Case 3, it is also difficult to estimate the pile quality from the velocity admittance curves for Case 4 because of the complex characteristics of nonuniform peaks with irregular distribution. As depicted in Figure 5(b), for the pile with necking segment located between the middle and the top of pile shaft, the orthokinetic reflected signals with reverse reflected signals will first appear at a distance of 6 m to the pile top, and the second-order orthokinetic reflected signals with reverse reflected signals caused by the necking segment will be received again at a distance of 12 m to the pile top. The peak values of the above-mentioned orthokinetic and reverse reflected signals increase as the necking degree of pile segment increases. For the pile with bulging segment located between the middle and the top of pile shaft, the reverse reflected signals with orthokinetic reflected signals will first occur at a distance of 6 m to the pile top, and the second-order reverse reflected signals with orthokinetic reflected signals caused by the bulging segment will be received again at a distance of 12 m to the pile top. The peak values of the above-mentioned reverse and orthokinetic reflected signals increase with the increase of the bulging degree of pile segment. It can be found that the larger the necking degree or bulging degree of variable cross-section segment, the larger the peak values of the above-reflected signals.

Influence of the variation degree of variable cross-section segment on the vertical dynamic response of pile head for Case 4.

Reflected wave signal curves

Case 5 (variable cross-section segment located between the middle and the end of pile shaft).

In this case, the pile shaft is also divided into three segments, which are marked by 1, 2, and 3 from pile end to pile top. Segment 1 and segment 3 have a uniform cross-section and satisfy the design requirements, and segment 2 has a variable cross-section. The lengths of each segment are l1=6m, l2=1m, and l3=13m, respectively. The other parameters are the same as those shown in Case 1.

Figure 6 illustrates the influence of variable degree of variable cross-section segment located between the middle and the end of pile shaft on the vertical dynamic response of pile head. Similar to Cases 3 and 4, it is difficult to evaluate the pile quality from the velocity admittance curves for Case 5 due to the complex characteristics of nonuniform peaks with irregular distribution. As shown in Figure 6(b), for the pile with necking segment located between the middle and the end of pile shaft, the orthokinetic reflected signals with reverse reflected signals will occur at a distance of 13 m to the pile top. The peak values of above-mentioned orthokinetic and reverse reflected signals increase with the increase of the necking degree of pile. For the pile with bulging segment located between the middle and the end of pile shaft, the reverse reflected signals with orthokinetic reflected signals will appear at a distance of 13 m to the pile top. The peak values of the above-mentioned reverse and orthokinetic reflected signals increase as the bulging degree of pile segment increases. It is worth noting that the larger the necking degree or bulging degree of variable cross-section segment, the larger the peak values of the above-reflected signals.

Influence of the variation degree of variable cross-section segment on the vertical dynamic response of pile head for Case 5.

Reflected wave signal curves

4.2. Vertical Vibration Characteristics of Variable Modulus Pile

During the construction of bored pile, the modulus of pile may be nonuniform along the pile shaft owing to the immature construction craft, improper construction management, or other adverse factors. Therefore, this section discusses the influence of the location and degree of nonuniform pile segment on the vertical dynamic response of pile head. First, the influence of the location of nonuniform pile segment on the vertical dynamic response of pile is investigated. The location and length of nonuniform pile segment are referred to the five cases shown in Section 4.1. The length, radius, and mass density of pile are 20 m, 0.5 m, and 2500kg/m3, respectively. The nonuniform characteristic of pile is reflected by varying the longitudinal wave velocity of normal pile segment and nonuniform pile segment, and the longitudinal wave velocity of normal pile segment and nonuniform pile segment are 3800m/s and 3000m/s, respectively.

Figure 7 presents the influence of the location of nonuniform pile segment on the vertical dynamic response of pile head. As shown in Figure 7(a), when the nonuniform pile segment is located at the top of pile (Case 1), the formant amplitude of velocity admittance curves gradually increases with the increase of frequency. For the other four cases, the oscillatory part of the velocity admittance curves displays increasing complex characteristics of nonuniform peaks with irregular distribution for the pile with variable modulus segment. From Figure 7(b), the influence of the location of nonuniform pile segment on the vertical dynamic response of pile head for the five cases is basically similar to that of pile with necking segment at corresponding location.

Influence of the location of nonuniform pile segment on the vertical dynamic response of pile head.

Reflected wave signal curves

Then, the influence of the degree of nonuniform pile segment on the vertical dynamic response of pile head is further discussed. The pile shaft is divided into three segments, which are marked by 1, 2, and 3 from pile end to pile top. Segment 1 and segment 3 are normal pile segments and have normal modulus satisfying the design requirements; segment 2 has nonnormal modulus. Denote the variable wave impedance ratio to describe the nonuniform degree of pile segment; that is, RV represents the ratio of wave impedance of normal pile segment to that of nonuniform pile segment. RV>1 and RV<1 mean the increase and decrease of wave impedance in the nonuniform pile segment, respectively. RA=1 indicates that the pile is uniform as a whole. The lengths of each segment are l1=6m, l2=1m, and l3=13m, respectively.

Figure 8 indicates the influence of the degree of nonuniform pile segment on the vertical dynamic response of pile head. As shown in Figure 8(a), similar to the regulation of pile with necking segment or bulging segment at the corresponding location in the pile shaft, there are complex characteristics of nonuniform peaks with irregular distribution in the velocity admittance curves when the wave impedance of nonuniform pile segment increases or decreases. The larger the variation range of wave impedance, the more obvious the complex oscillation. Therefore, it is also difficult to intuitively assess the nonuniform characteristics of pile shaft only from the velocity admittance curves. As depicted in Figure 8(b), as the wave impedance of nonuniform pile segment decreases, the variation regulation of reflected wave curves is basically similar to that of pile with necking segment at corresponding location. Meanwhile, as the wave impedance of nonuniform pile segment increases, the variation regulation of reflected wave curves is basically similar to that of pile with bulging at corresponding location.

Influence of the degree of nonuniform pile segment on the vertical dynamic response of pile head.

Reflected wave signal curves

5. Conclusions

Based on the plane strain model and fictitious soil pile model, this paper derives the analytical solution of velocity admittance in the frequency domain and the corresponding semianalytical solution of the velocity response of pile in the time domain. Based on the solutions, a parametric study is conducted to investigate the vertical dynamic response of pile with variable cross-section segment and variable modulus segment, and the following conclusions are obtained:

When there are necking or bulging segments located between the top and the end of pile shaft, the oscillatory part of the velocity admittance curves displays increasing complex characteristics of nonuniform peaks with irregular distribution. Therefore, it is difficult to intuitively assess the construction quality of pile from the velocity admittance curves for these conditions.

For the pile with necking segment or decreasing modulus segment located between the top and the end of pile shaft, the orthokinetic reflected signals are first received between the incident wave signals and preliminary orthokinetic reflected signals from pile end, and the reverse reflected signals will occur after the orthokinetic reflected signals.

For the pile with bulging segment or increasing modulus segment located between the top and the end of pile shaft, the reverse reflected signals will first appear between the incident wave signals and preliminary orthokinetic reflected signals from pile end, and the orthokinetic reflected signals are received after the reverse reflected signals.

The larger the necking degree or bulging degree of variable cross-section segment, the larger the peak values of the reflected signals between the incident wave signals and preliminary orthokinetic reflected signals from pile end.

The larger the variation range of the modulus of variable modulus pile segment, the larger the peak values of the reflected signals between the incident wave signals and preliminary orthokinetic reflected signals from pile end.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research is funded by the National Natural Science Foundation of China (Grants nos. 51309207, 51322807, and 51678547) and the China Postdoctoral Science Foundation Funded Project (Grant no. 2016M600711). The authors also wish to acknowledge the support of the Research Funds provided by MOE Engineering Research Center of Rock-Soil Drilling & Excavation and Protection (Grant no. 201402) and the Fundamental Research Funds for the Central Universities-Cradle Plan for 2015 (Grant no. CUGL150411).

Liao S. T. Roesset J. M. Dynamic response of intact piles to impulse loads International Journal for Numerical and Analytical Methods in Geomechanics 1997 21 4 255 275 10.1002/(sici)1096-9853(199704)21:460;255::aid-nag86962;3.0.co;2-j 2-s2.0-0031126285 Yanlong L. Shouyi L. Yang Y. Xing T. Temperature stress and surface insulation measures of concrete face slabs during cold wave period International Journal of Civil Engineering 2015 13 4 501 507 2-s2.0-84949434289 Li Y. Sun Y. Li B. Xu Z. Penalty function-based method for obtaining a reliability indicator of gravity dam stability Computers and Geotechnics 2017 81 19 25 10.1016/j.compgeo.2016.07.003 Tian H. Ziegler M. Kempka T. Physical and mechanical behavior of claystone exposed to temperatures up to 1000°C International Journal of Rock Mechanics and Mining Sciences 2014 70 144 153 10.1016/j.ijrmms.2014.04.014 2-s2.0-84901197204 Van Koten H. Middendorp P. Van Brederode P. An analysis of dissipative wave propagation in a pile International Seminar on the Application of Stress-Wave Theory on Piles, Stockholm 1980 Rotterdam, The Netherlands A.A. Balkema Randolph M. F. Deeks A. J. Dynamic and static soil models for axial response Proceedings of the 4th International Conference on the Application of Stress Wave Theory to Piles September 1992 The Hague, The Netherlands 3 14 Michaelides O. Gazetas G. Bouckovalas G. Chrysikou E. Approximate non-linear dynamic axial response of piles Geotechnique 1998 48 1 33 53 10.1680/geot.1998.48.1.33 2-s2.0-0032007409 Yesilce Y. Catal H. H. Free vibration of piles embedded in soil having different modulus of subgrade reaction Applied Mathematical Modelling 2008 32 5 889 900 10.1016/j.apm.2007.02.015 ZBL1130.74032 2-s2.0-38349110268 Wang K. Wu W. Zhang Z. Leo C. J. Vertical dynamic response of an inhomogeneous viscoelastic pile Computers and Geotechnics 2010 37 4 536 544 10.1016/j.compgeo.2010.03.001 2-s2.0-77952547564 Zhang Q.-Q. Zhang Z.-M. A simplified nonlinear approach for single pile settlement analysis Canadian Geotechnical Journal 2012 49 11 1256 1266 10.1139/T11-110 2-s2.0-84865401517 Zhang Q.-Q. Liu S.-W. Zhang S.-M. Zhang J. Wang K. Simplified non-linear approaches for response of a single pile and pile groups considering progressive deformation of pile-soil system Soils and Foundations 2016 56 3 473 484 10.1016/j.sandf.2016.04.013 2-s2.0-84965081505 Novak M. Dynamic stiffness and damping of piles Canadian Geotechnical Journal 1974 11 4 574 598 10.1139/t74-059 Rajapakse R. K. N. D. Shah A. H. On the longitudinal harmonic motion of an elastic bar embedded in an elastic half-space International Journal of Solids and Structures 1987 23 2 267 285 10.1016/0020-7683(87)90060-6 ZBL0604.73060 2-s2.0-0023271405 Mamoon S. M. Banerjee P. K. Time-domain analysis of dynamically loaded single piles Journal of Engineering Mechanics 1992 118 1 140 160 10.1061/(ASCE)0733-9399(1992)118:1(140) 2-s2.0-84904994637 El Naggar M. H. Novak M. Nonlinear analysis for dynamic lateral pile response Soil Dynamics and Earthquake Engineering 1996 15 4 233 244 10.1016/0267-7261(95)00049-6 2-s2.0-0030158838 Militano G. Rajapakse R. K. N. D. Dynamic response of a pile in a multi-layered soil to transient torsional and axial loading Geotechnique 1999 49 1 91 109 10.1680/geot.1999.49.1.91 2-s2.0-0033081188 Wu W. Jiang G. Dou B. Leo C. J. Vertical dynamic impedance of tapered pile considering compacting effect Mathematical Problems in Engineering 2013 2013 9 304856 10.1155/2013/304856 2-s2.0-84893859543 Nogami T. Novak M. Soil‐pile interaction in vertical vibration Earthquake Engineering & Structural Dynamics 1976 4 3 277 293 10.1002/eqe.4290040308 2-s2.0-0016892789 Nogami T. Konagai K. Time domain axial response of dynamically loaded single piles Journal of Engineering Mechanics 1986 112 11 1241 1252 10.1061/(ASCE)0733-9399(1986)112:11(1241) 2-s2.0-0022807309 Senjuntichai T. Mani S. Rajapakse R. K. N. D. Vertical vibration of an embedded rigid foundation in a poroelastic soil Soil Dynamics and Earthquake Engineering 2006 26 6-7 626 636 10.1016/j.soildyn.2006.01.013 2-s2.0-33746014096 Hu C. B. Wang K. H. Xie K. H. Time domain analysis of vertical dynamic response of a pile considering the effect of soil-pile interaction Chinese Journal of Computational Mechanics 2004 21 8 392 399 2-s2.0-7244236692 Wang K. Zhang Z. Leo C. J. Xie K. Dynamic torsional response of an end bearing pile in saturated poroelastic medium Computers and Geotechnics 2008 35 3 450 458 10.1016/j.compgeo.2007.06.013 2-s2.0-42649112088 Wu W. B. Wang K. H. Zhang Z. Q. Leo C. J. Soil-pile interaction in the pile vertical vibration considering true three-dimensional wave effect of soil International Journal for Numerical and Analytical Methods in Geomechanics 2013 37 17 2860 2876 10.1002/nag.2164 2-s2.0-84887824861 Wu W. Jiang G. Huang S. Leo C. J. Vertical dynamic response of pile embedded in layered transversely isotropic soil Mathematical Problems in Engineering 2014 2014 12 126916 10.1155/2014/126916 MR3240658 2-s2.0-84934992833 Zheng C. Liu H. Ding X. Kouretzis G. P. Sloan S. W. Poulos H. G. Non-axisymmetric response of piles in low-strain integrity testing Géotechnique 2017 67 2 181 186 10.1680/jgeot.15.P.193 Zheng C. Kouretzis G. P. Ding X. Liu H. Poulos H. G. Three-dimensional effects in low-strain integrity testing of piles: analytical solution Canadian Geotechnical Journal 2015 53 2 225 235 10.1139/cgj-2015-0231 2-s2.0-84938905670