Consolidation of soil induced by pile installation considering disturbance effect

<!-- metadata commented in wiki content


<span id='_Hlk51705514'></span><span id='_Hlk41492253'></span><span id='_Hlk55074535'></span>

==Consolidation of soil induced by pile installation considering disturbance effect==

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
by</div>

==Authors: E-mail:==

<span id='_Hlk82103515'></span><span id='_Hlk82103504'></span>Ping Li<sup>1</sup> [mailto:geopingli@163.com geopingli@163.com]

Zhijian Chen<sup>2</sup> [mailto:zhijianchen_geo@163.com zhijianchen_geo@163.com]

Yi Ding<sup>3</sup> [mailto:nsscm_yiding@163.com nsscm_yiding@163.com]

==Affiliations:==

<sup>1</sup>College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China.

<sup>2</sup>School of Earth Sciences and Engineering, Hohai University, Nanjing 210098, China.

<span id='OLE_LINK4'></span><sup>3</sup> Nanjing Supervision Station for Construction Market, Nanjing 210098, China.

<span id='_Hlk73473327'></span>'''Corresponding author: '''Zhijian Chen<br/>
-->
==Abstract==

In practice, the consolidation of soil around the pile has a great influence on the time-dependent bearing capacity of pile. However, most of the consolidation theory of soil around the pile neglects the disturbance effect of pile-driving on surrounding soil and regards the soil as homogeneous, which overestimates the consolidation efficiency of the soil, and obtains a higher pile bearing capacity. In view of this, a consolidation model of soil around a pipe pile considering soil disturbance effect is presented in this paper. Fourier transform and separation of variables are used to obtain the analytical solution, and then the solution is verified by degradation analysis and Finite Difference Method (FDM). Firstly, the radial and vertical distribution of excess pore pressure generated after pile-driving is analyzed. In the radial direction, the excess pore pressure decreases rapidly from the radius of pipe pile to the radius of the disturbed zone, then slowly decays to 0 from'' ''the radius of disturbed zone to the influencing radius of pipe pile. In the vertical direction, the excess pore pressure along the vertical direction grows linearly. Subsequently, the variation of the average excess pore pressure at the pile-soil interface with the permeability coefficient and radius of disturbance zone are analyzed. The permeability coefficient of disturbance zone has a remarkable negative effect on the excess pore pressure during the whole consolidation period. Increasing the radius of the disturbance zone will hinder the dissipation of the excess pore pressure in the intermediate and later stages. Finally, the validity of the proposed analytical solution is illustrated by comparing with the laboratory results.

'''Keywords''': Pipe pile, disturbance effect, consolidation theory, analytical solution<br/>

==1. Introduction==

Coastal areas are mostly soft soil foundation, which needs to be treated before construction. Static pile technique is a common method to deal with soft soil foundation, which has been widely used in the world [1-3]. Although the bearing capacity of the foundation can be improved by pressing the static pile into the soil by the static pile pressing machine, it is inevitable to squeeze the soil in the process of pile installation, resulting in disturbance of the surrounding soil around the static pile [4-6]. Further, pile installation could produce high excess pore pressure, which results in the slow development of pile bearing capacity [7]. High excess pore pressure could also cause some negative effects on the pile and surrounding soil, including deviation of pile and ground heave [8-10]. Therefore, it is necessary to put forward a set of prediction scheme to evaluate the excess pore pressure after pile installation, so that engineers can deal with the risks in advance.

Axisymmetric consolidation model is often used to evaluate the dissipation of excess pore pressure after pile-driving. Randolph and Wroth [10] first proposed an analytical solution of axisymmetric consolidation around an impermeable pile considering radial drainage through pile-driving test. Subsequently, in the theoretical study on the axisymmetric consolidation of the soil around the impermeable pile, Guo  [11] proposed the analytical solution considering viscoelastic soil, and Gao and Shi [12] proposed the analytical solution considering both radial and vertical drainage. Recently, a permeable pile technique has been developed to accelerate the consolidation rate of foundation soil. Therefore, some scholars have further studied axisymmetric consolidation around a permeable pile [13,14].

It should be emphasized that all above-mentioned studies are limited to undisturbed soils. In fact, the existing field data and test data show that the disturbance effect on the soil around the pile is inevitable in the process of pile installation [15,16]. In general, in the process of pile installation, radial compression will occur to the soil around the pile, leading to the decrease of the radial permeability coefficient of the disturbed soil, which will lead to the decrease of the radial drainage performance of the soil. When the vertical drainage distance is long, the soil drainage path is mainly radial. In this case, the disturbance effect of soil should not be ignored.

Similar to pile installation, the prefabricated vertical drains (i.e., PVDs) installation has been fully studied <sup>[17-20]</sup>. Compared with the smear effect by PVD installation [21-25], the disturbance effect induced by pile installation is much more significant due to its bigger radius. Thus, it is of significance to study the consolidation of soil around the pile considering disturbance effect. Additionally, the growth rate of pile bearing capacity could be overestimated by using the current consolidation model, which leads to some unknown risks (e.g., excessive settlement and insufficient bearing capacity) and reducing the accuracy of building’s reliability analysis [9,26,27]. Therefore, whether it is the study of consolidation theory of soil around pile or the subsequent analysis of pile bearing capacity, the study of consolidation theory considering disturbance effect is of great significance. Unfortunately, there is no an analytical solution to axisymmetric consolidation of soil considering the disturbance effect around a pipe pile.

In order to better predict the dissipation process of excess pore pressure and the increase process of pile bearing capacity after piling, it is necessary to further optimize the existing consolidation theory. On the condition that the disturbance effect cannot be ignored, this paper presents an analytical solution to axisymmetric consolidation of soil considering the disturbance effect around a pipe pile on the basis of previous studies. Fourier transform and separation of variables are used to solve the problem, and the solution is calibrated by degradation analysis and FEM. Subsequently, the distribution of excess pore pressure along vertical and radial directions are analyzed by graphical interpretation. Moreover, the concept of the dimensionless average excess pore pressure at pile-soil interface is put forward and the influence of disturbance parameters on it is analyzed. Most importantly, the validity of the solution is demonstrated by comparing with test data.

==2. Mathematical model==

===2.1 Model description===

As shown in [[#img-1|Figure 1]], all involved parameters include the depth of foundation,  <math>H</math>, the radius of pipe pile,  <math>r_{0}</math>, the radius of the disturbed zone,  <math>r_{d}</math>, and the influencing radius of pipe pile,  <math>r_{e}</math>. In this model, the following assumptions are given:

(1) The soil around the pipe pile is elastic saturated soil, and the pipe pile is impermeable pile;

(2) The soil around the pile only has disturbance effect in radial direction;

(3) The soil parameters used in the model remain unchanged during the consolidation process.

<div id='img-1'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[Image:Review_932178682507-image5.png|246px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figure 1'''. Schematic diagram of consolidation model for the soil around a pipe pile 
|}

===2.2 Governing equations===

According to the findings of Randolph and Wroth [10], the excess pore pressure generated during pile-driving will dissipate in both radial and vertical directions. Therefore, the following governing equations are established:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum571777'></span>  <math>\begin{array}{c}
\displaystyle\frac{\partial u_1}{\partial t}=C_{d}\left(\displaystyle\frac{{\partial }^2u_1}{\partial r^2}+\displaystyle\frac{1}{r}\frac{\partial u_1}{\partial r}\right){+}C_{v}\displaystyle\frac{{\partial }^2u_1}{\partial z^2}\\
\displaystyle\frac{\partial u_2}{\partial t}=C_{h}\left(\displaystyle\frac{{\partial }^2u_2}{\partial r^2}+\displaystyle\frac{1}{r}\displaystyle\frac{\partial u_2}{\partial r}\right){+}C_{v}\displaystyle\frac{{\partial }^2u_2}{\partial z^2}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}

where  <math>u_{1}</math> and  <math>u_{2}</math> are the excess pore pressure of disturbed zone (i.e.,  <math>r_0\leq r\leq r_{d}</math>) and undisturbed zone  (i.e.,  <math>r_{d}\leq r\leq r_{e}</math>), respectively;  <math>C_{h}</math> (i.e.,  <math>k_{h}/\left({\gamma }_{w}m_{v}\right)</math>) and  <math>C_{d}</math> (i.e.,  <math>k_{d}/\left({\gamma }_{w}m_{v}\right)</math>) represent the consolidation coefficient in radial direction of undisturbed and disturbed zones, respectively;  <math>C_{v}</math> (i.e.,  <math>k_{v}/\left({\gamma }_{w}m_{v}\right)</math>) represent the consolidation coefficient in vertical direction of the undisturbed zone;  <math>k_{h}</math> and  <math>k_{d}</math> represent the permeability coefficient in radial direction of disturbed and undisturbed zones, respectively;  <math>k_{v}</math> represent the permeability coefficient in vertical direction of the influencing zone;  <math>m_{v}</math> is the volume compression coefficient of soil;  <math>{\gamma }_{w}</math> is the unit weight of water;  <math>r</math> and  <math>z</math> are the radial distance from the center of pile and the depth from the ground surface, respectively; and  <math>t</math> is consolidation time.

===2.3 Initial and boundary conditions===

The initial pore pressure distribution of soil around a pile after pile-driving can generally be obtained through numerical simulation or experimental test [28-30]. Here, the initial excess pore pressure is generalized expressed as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum823648'></span>  <math>{u_1\vert }_{t=0}={u_2\vert }_{t=0}=u_0\left(r,z\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}

The pile is impermeable and the excess pore pressure outside the influencing radius  <math>r_\mbox{e}</math> affected by pipe pile is ignored, it yields

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum743846'></span>  <math>\frac{\partial u_1}{\partial r}\Bigg|_{r=r_0} ={u_2\vert }_{r=r_\mbox{e}}=</math><math>0</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

The interface between the disturbed and undisturbed zones needs to satisfy the flow continuity condition, that is

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum201641'></span>  <math>\begin{array}{c}
{u_1\vert }_{r=r_\mbox{d}}={u_2\vert }_{r=r_\mbox{d}}\\
k_\mbox{d}\displaystyle\frac{\partial u_1}{\partial r}\Bigg|_{r=r_\mbox{d}}=k_\mbox{h}\displaystyle\frac{\partial u_2}{\partial r}\Bigg|_{r=r_\mbox{d}}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

The top of the ground is permeable, while the bottom is impermeable. Therefore, the vertical boundary can be expressed as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum798795'></span>  <math>\begin{array}{c}
{u_1\vert }_{z=0}=\displaystyle\frac{\partial u_1}{\partial z}\Bigg|_{z=H}=0\\
{u_2\vert }_{z=0}=\displaystyle\frac{\partial u_2}{\partial z}\Bigg|_{z=H}=0
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}

===2.4 Solving procedures===

According to the boundary condition of Eq.<span id='cite-ZEqnNum798795'></span>[[#ZEqnNum798795|(5)]], the finite Fourier sine transform is performed on the governing equation (i.e., Eq.<span id='cite-ZEqnNum571777'></span>[[#ZEqnNum571777|(1)]]) and its remaining solving conditions (i.e., Eq.<span id='cite-ZEqnNum823648'></span>[[#ZEqnNum823648|(2)]]-<span id='cite-ZEqnNum201641'></span>[[#ZEqnNum201641|(4)]]), which yields

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum149568'></span>  <math> \frac{\partial \overline{u}_1}{\partial t}=C_d \left(\frac{\partial^2 \overline{u}_1}{\partial r^2}+ \frac{1}{r} \frac{\partial \overline{u}_1}{\partial r}\right) - C_v \left(\frac{N_n}{H} \right)^2 \overline{u}_1, \quad r_0\le r\le r_d ; </math>
|-
| <math> \frac{\partial {\overline{u}}_2}{\partial t}=C_h \left( \frac{\partial^2 {\overline{u}}_2}{\partial r^2}+ \frac{1}{r} \frac{\partial {\overline{u}}_2}{\partial r}\right) - C_v \left(\frac{N_n}{H} \right)^2\overline{u}_2, \quad r_d\le r\le r_e ; </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum969494'></span>  <math>{{\overline{u}}_1\vert }_{t=0}={{\overline{u}}_2\vert }_{t=0}=</math><math>{\int }_0^Hu_0\left(r,z\right)\sin\frac{N_nz}{H}{d}z</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum791600'></span>  <math>\displaystyle\frac{\partial {\overline{u}}_1}{\partial r}\Bigg|_{r=r_0}={{\overline{u}}_2\vert }_{r=r_{e}}=0</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum399916'></span>  <math>\begin{array}{c}
{{\overline{u}}_1\vert }_{r=r_{d}}={{\overline{u}}_2\vert }_{r=r_{d}}\\
k_{d}\displaystyle\frac{\partial {\overline{u}}_1}{\partial r}\Bigg|_{r=r_{d}}=k_{h}\displaystyle\frac{\partial {\overline{u}}_2}{\partial r}\Bigg|_{r=r_{d}}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}

where  <math>\overline{u}_i=\int_0^Hu_i \sin\frac{N_nz}{H}{d}z,\, i=1,2</math>; and  <math>N_n=\left(2n-1\right)\pi/2</math>.

According to the separation of variables method,  <math>{\overline{u}}_1</math> and  <math>{\overline{u}}_2</math> can be written as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum108209'></span>  <math>\begin{array}{c}
{\overline{u}}_1=R_{1mn}\left(r\right)T_{1mn}\left(t\right)\\
{\overline{u}}_2=R_{2mn}\left(r\right)T_{2mn}\left(t\right)
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
|}

where  <math>R_{1mn}</math> and  <math>R_{2mn}</math> are functions of  <math>r</math>;  <math>T_{1mn}</math> and  <math>T_{2mn}</math> are functions of  <math>t</math>.

Substituting Eq.(10) into Eq.(6), one can obtain

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum505543'></span>  <math>\begin{array}{c}
\displaystyle\frac{\displaystyle\frac{{d}T_{1mn}}{{d}t}}{C_{d}T_{1mn}}+\displaystyle\frac{C_{v}}{C_{d}}{\alpha }_n^2=\displaystyle\frac{1}{R_{1mn}}\left(\displaystyle\frac{{d}^2R_{1mn}}{{d}r^2}+\displaystyle\frac{1}{r}\displaystyle\frac{{d}R_{1mn}}{{d}r}\right)=-{\beta }_{1mn}^2\\
\displaystyle\frac{\displaystyle\frac{{d}T_{2mn}}{{d}t}}{C_{h}T_{2mn}}+\displaystyle\frac{C_{v}}{C_{h}}{\alpha }_n^2=\displaystyle\frac{1}{R_{2mn}}\left(\displaystyle\frac{{d}^2R_{2mn}}{{d}r^2}+\displaystyle\frac{1}{r}\displaystyle\frac{{d}R_{2mn}}{{d}r}\right)=-{\beta }_{2mn}^2
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (11)
|}

where  <math>{\alpha }_n=N_n/H</math>; and  <math>{\beta }_{1mn}</math> and  <math>{\beta }_{2mn}</math> are the separation constants.

It can be found that Eq.<span id='cite-ZEqnNum505543'></span>[[#ZEqnNum505543|(11)]] is an ordinary differential equation, and the solution to Eq.<span id='cite-ZEqnNum505543'></span>[[#ZEqnNum505543|(11)]] can be easily obtained as follows

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum188101'></span>  <math>\begin{array}{l}
R_{1mn}=a_{1mn}{J}_0\left({\beta }_{1mn}r\right)+b_{1mn}{Y}_0\left({\beta }_{1mn}r\right)\\
T_{1mn}=c_{1mn}{e}^{-\left(C_{d}{\beta }_{1mn}^2+C_{v}{\alpha }_n^2\right)t}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum833305'></span>  <math>\begin{array}{l}
R_{2mn}=a_{2mn}{J}_0\left({\beta }_{2mn}r\right)+b_{2mn}{Y}_0\left({\beta }_{2mn}r\right)\\
T_{2mn}=c_{2mn}{e}^{-\left(C_{h}{\beta }_{2mn}^2+C_{v}{\alpha }_n^2\right)t}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
|}

where  <math>{J}_0</math>,  <math>{Y}_0</math> are the zero-order Bessel functions of the first- and second- kind, respectively.

Combining Eqs.<span id='cite-ZEqnNum188101'></span>[[#ZEqnNum188101|(12)]] and <span id='cite-ZEqnNum833305'></span>[[#ZEqnNum833305|(13)]],  <math>{\overline{u}}_1</math> and  <math>{\overline{u}}_2</math> can be rewritten as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum354201'></span>  <math>\begin{array}{c}
{\overline{u}}_1=\displaystyle\sum_{m=1}^{\infty }\left[A_{1mn}{J}_0\left({\beta }_{1mn}r\right)+B_{1mn}{Y}_0\left({\beta }_{1mn}r\right)\right]{e}^{-\left(C_{d}{\beta }_{1mn}^2+C_{v}{\alpha }_n^2\right)t}\\
{\overline{u}}_2=\displaystyle\sum_{m=1}^{\infty }\left[A_{2mn}{J}_0\left({\beta }_{2mn}r\right)+B_{2mn}{Y}_0\left({\beta }_{2mn}r\right)\right]{e}^{-\left(C_{h}{\beta }_{2mn}^2+C_{v}{\alpha }_n^2\right)t}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (14)
|}

Substituting Eq.<span id='cite-ZEqnNum354201'></span>[[#ZEqnNum354201|(14)]] into the boundary condition of Eq.<span id='cite-ZEqnNum791600'></span>[[#ZEqnNum791600|(8)]], it yields

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum665918'></span>  <math>\begin{array}{c}
{\overline{u}}_1=\displaystyle\sum_{m=1}^{\infty }A_{1mn}R_1\left({\beta }_{1mn}r\right){e}^{-\left(C_{d}{\beta }_{1mn}^2+C_{v}{\alpha }_n^2\right)t}\\
{\overline{u}}_2=\displaystyle\sum_{m=1}^{\infty }A_{2mn}R_2\left({\beta }_{2mn}r\right){e}^{-\left(C_{h}{\beta }_{2mn}^2+C_{v}{\alpha }_n^2\right)t}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (15)
|}

where  <math>{\chi }_1\left({\beta }_{1mn}\right)=\frac{{J}_1\left({\beta }_{1mn}r_0\right)}{{Y}_1\left({\beta }_{1mn}r_0\right)}</math>;  <math>{\chi }_2\left({\beta }_{2mn}\right)=\frac{{J}_0\left({\beta }_{2mn}r_{e}\right)}{{Y}_0\left({\beta }_{2mn}r_{e}\right)}</math>;  <math>R_1\left({\beta }_{1mn}r\right)={J}_0\left({\beta }_{1mn}r\right)-</math><math>{\chi }_1\left({\beta }_{1mn}\right){Y}_0\left({\beta }_{1mn}r\right)</math>; and  <math>R_2\left({\beta }_{2mn}r\right)={J}_0\left({\beta }_{2mn}r\right)-</math><math>{\chi }_2\left({\beta }_{2mn}\right){Y}_0\left({\beta }_{2mn}r\right)</math>.

Substituting Eq.<span id='cite-ZEqnNum665918'></span>[[#ZEqnNum665918|(15)]] into the continuous boundary condition (i.e., Eq.<span id='cite-ZEqnNum399916'></span>[[#ZEqnNum399916|(9)]]), it yields

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum164985'></span>  <math>R_1\left({\beta }_{1mn}r_{d}\right)A_{1mn}=R_2\left({\beta }_{2mn}r_{d}\right)A_{2mn}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (16)
|}

and

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum741519'></span>  <math>\left\{ \begin{array}{c}
C_{d}{\beta }_{1mn}^2=C_{h}{\beta }_{2mn}^2\\
\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}=\frac{k_{h}{\beta }_{2mn}P_2\left({\beta }_{2mn}r_{d}\right)}{k_{d}{\beta }_{1mn}P_1\left({\beta }_{1mn}r_{d}\right)}
\end{array}\right.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (17)
|}

where  <math>P_1\left({\beta }_{1mn}r\right)={J}_1\left({\beta }_{1mn}r\right)-</math><math>{\chi }_1\left({\beta }_{1mn}\right){Y}_1\left({\beta }_{1mn}r\right)</math>; and  <math>P_2\left({\beta }_{2mn}r\right)={J}_1\left({\beta }_{2mn}r\right)-</math><math>{\chi }_2\left({\beta }_{2mn}\right){Y}_1\left({\beta }_{2mn}r\right)</math>.

Combined with Eq.<span id='cite-ZEqnNum164985'></span>[[#ZEqnNum164985|(16)]], Eq.<span id='cite-ZEqnNum665918'></span>[[#ZEqnNum665918|(15)]] can be rewritten as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum775505'></span>  <math>\begin{array}{l}
{\overline{u}}_1=\displaystyle\sum_{m=1}^{\infty }A_{2mn}\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}R_1\left({\beta }_{1mn}r\right){e}^{-\left(C_{d}{\beta }_{1mn}^2+C_{v}{\alpha }_n^2\right)t}\\
{\overline{u}}_2=\displaystyle\sum_{m=1}^{\infty }A_{2mn}R_2\left({\beta }_{2mn}r\right){e}^{-\left(C_{h}{\beta }_{2mn}^2+C_{v}{\alpha }_n^2\right)t}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (18)
|}

<math>{\beta }_{1mn}</math> and  <math>{\beta }_{2mn}</math> can be solved with Eq.<span id='cite-ZEqnNum741519'></span>[[#ZEqnNum741519|(17)]]. Then, using the initial condition of Eq.<span id='cite-ZEqnNum969494'></span>[[#ZEqnNum969494|(7)]] and orthogonality of eigenfunctions, one can obtain

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum753186'></span>  <math>\begin{array}{c}
A_{2mn}\displaystyle\int_{r_0}^{r_{d}}r{\left[\displaystyle\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}R_1\left({\beta }_{1mn}r\right)\right]}^2{d}r+A_{2mn}\displaystyle\int_{r_{d}}^{r_{e}}rR_2^2\left({\beta }_{2mn}r\right){d}r\\
=\displaystyle\int_{r_0}^{r_{d}}r\displaystyle\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}R_1\left({\beta }_{1mn}r\right)\displaystyle\int _0^Hu_0\left(r,z\right)\sin{\alpha }_nz{d}z{d}r+\displaystyle\int_{r_{d}}^{r_{e}}rR_2\left({\beta }_{2mn}r\right){\int }_0^Hu_0\left(r,z\right)\sin{\alpha }_nz{d}z{d}r
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (19)
|}

According to Eq.<span id='cite-ZEqnNum753186'></span>[[#ZEqnNum753186|(19)]], the coefficient  <math>A_{2mn}</math> can be obtained. Applying finite Fourier inverse transform to Eq.<span id='cite-ZEqnNum775505'></span>[[#ZEqnNum775505|(18)]], the general solutions can be written as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{l}
u_1=\frac{2}{H}\displaystyle\sum_{m=1}^{\infty }\displaystyle\sum_{n=1}^{\infty }A_{2mn}\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}R_1\left({\beta }_{1mn}r\right)\sin{\alpha }_nz{e}^{-\left(C_{d}{\beta }_{1mn}^2+C_{v}{\alpha }_n^2\right)t}\\
u_2=\frac{2}{H}\displaystyle\sum_{m=1}^{\infty }\displaystyle\sum_{n=1}^{\infty }A_{2mn}R_2\left({\beta }_{2mn}r\right)\sin{\alpha }_nz{e}^{-\left(C_{h}{\beta }_{2mn}^2+C_{v}{\alpha }_n^2\right)t}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (20)
|}

The initial excess pore pressure is mainly distributed in the plastic zone. It decays along the radius direction and increases along the depth. It is also worthy to highlight the fact that the initial excess pore pressure is calculated after pile-driving, at which time the excess pore pressure of the upper part of surface soil has been completely dissipated [14]. So, the specific expression form of the initial excess pore pressure is given as follows

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum163277'></span>  <math>u_0\left(r,z\right)=a_1\left(z-h_0\right)ln\frac{r_{p}}{r}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (21)
|}

where  <math>r_{p}</math> is the radius of plastic zone; <math>a_1</math> and  <math>h_0</math> are constant to be determined. The distribution range of the initial excess pore pressure is  <math>\left\{z\geq h_0;r_0\leq r\leq r_{p}\right\}</math>.

Then, the coefficient  <math>A_{2mn}</math> can be specified as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>A_{2mn}=\frac{a_1\left(\sin{\alpha }_nH-\sin{\alpha }_nh_0\right)}{{\alpha }_n{}^2}\frac{\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}{\delta }_1+{\delta }_2}{\frac{R_2^2\left({\beta }_{2mn}r_{d}\right)}{R_1^2\left({\beta }_{1mn}r_{d}\right)}{\delta }_3+{\delta }_4}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (22)
|}

where

<math>{\delta }_1=\frac{1}{{\beta }_{1mn}}\left\{\frac{1}{{\beta }_{1mn}}\left[R_1\left({\beta }_{1mn}r_{0}\right)-\right. \right. </math><math>\left. \left. R_1\left({\beta }_{1mn}r_{d}\right)\right]+\right. </math><math>\left. r_{d}ln\frac{r_{p}}{r_{d}}P_1\left({\beta }_{1mn}r_{d}\right)-\right. </math><math>\left. r_0ln\frac{r_{p}}{r_0}P_1\left({\beta }_{1mn}r_{0}\right)\right\}</math>;

<math>{\delta }_2=\frac{1}{{\beta }_{2mn}}\left\{\frac{1}{{\beta }_{2mn}}\left[R_2\left({\beta }_{2mn}r_{d}\right)-\right. \right. </math><math>\left. \left. R_2\left({\beta }_{2mn}r_{p}\right)\right]+\right. </math><math>\left. r_{p}ln\frac{r_{p}}{r_{p}}P_2\left({\beta }_{2mn}r_{p}\right)-\right. </math><math>\left. r_{d}ln\frac{r_{p}}{r_{d}}P_2\left({\beta }_{2mn}r_{d}\right)\right\}</math>;

<math>{\delta }_3=\frac{r_{d}{}^2}{2}\left[P_1^2\left({\beta }_{1mn}r_{d}\right)+\right. </math><math>\left. R_1^2\left({\beta }_{1mn}r_{d}\right)\right]-</math><math>\frac{r_{0}{}^2}{2}\left[P_1^2\left(m,n,r_{0}\right)+\right. </math><math>\left. R_1^2\left(m,n,r_{0}\right)\right]</math>;

<math>{\delta }_4=\frac{r_{e}{}^2}{2}\left[P_2{\left({\beta }_{2mn}r_{e}\right)}^2+\right. </math><math>\left. R_2^2\left({\beta }_{2mn}r_{e}\right)\right]-</math><math>\frac{r_{d}{}^2}{2}\left[P_2^2\left({\beta }_{2mn}r_{d}\right)+\right. </math><math>\left. R_2^2\left({\beta }_{2mn}r_{d}\right)\right]</math>.

==3. Degradation analysis and verification by FEM==

The verification and subsequent parameter sensitivity analysis using the parameters defined below. The parameters not mentioned can be obtained in the graphical analysis.

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\begin{array}{ccc}
m_\mbox{v}={10}^{-2}\mbox{ }\mbox{m}^2/\mbox{kN} & {\gamma }_\mbox{w}=10\mbox{kN}/\mbox{m}^3 & k_\mbox{v}=2\times {10}^{-9}\mbox{m}/\mbox{s}\\
r_0=0.25\mbox{ m} & r_\mbox{p}=1.2\times r_\mbox{d} & r_\mbox{e}=20\times r_\mbox{0}\\
H=20\mbox{ m} & h_0=0.5\mbox{ m} & a_1=5\mbox{kPa}/\mbox{m}
\end{array}</math>
|}
|}

===3.1 Degradation analysis===

When  <math>k_{h}</math> and  <math>k_{d}</math> are equal, it can be seen from Eq.<span id='cite-ZEqnNum741519'></span>[[#ZEqnNum741519|(17)]] that

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum989253'></span>  <math>{\beta }_{1mn}={\beta }_{2mn}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (23)
|}

and

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum808683'></span>  <math>{\chi }_1\left({\beta }_{1mn}\right)={\chi }_2\left({\beta }_{2mn}\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (24)
|}

<math>{\beta }_{1mn}</math> and  <math>{\beta }_{2mn}</math> are obtained by solving Eq.<span id='cite-ZEqnNum808683'></span>[[#ZEqnNum808683|(24)]]. Substituting Eq. <span id='cite-ZEqnNum989253'></span>[[#ZEqnNum989253|(23)]] and Eq.<span id='cite-ZEqnNum808683'></span>[[#ZEqnNum808683|(24)]] back into  <math>R_1\left({\beta }_{1mn}r\right)</math>,  <math>R_2\left({\beta }_{2mn}r\right)</math>,  <math>P_1\left({\beta }_{1mn}r\right)</math> and  <math>P_2\left({\beta }_{2mn}r\right)</math>, it yields

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <span id='ZEqnNum350685'></span>  <math>\left\{ \begin{array}{c}
R_1\left({\beta }_{1mn}r\right)=R_2\left({\beta }_{2mn}r\right)\\
P_1\left({\beta }_{1mn}r\right)=P_2\left({\beta }_{2mn}r\right)
\end{array}\right.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (25)
|}

Combining Eq.<span id='cite-ZEqnNum164985'></span>[[#ZEqnNum164985|(16)]] with Eq.<span id='cite-ZEqnNum350685'></span>[[#ZEqnNum350685|(25)]], one can obtain

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>A_{1mn}=A_{2mn}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (26)
|}

So, it is found that  <math>{\overline{u}}_1</math> and  <math>{\overline{u}}_2</math> in Eq.<span id='cite-ZEqnNum665918'></span>[[#ZEqnNum665918|(15)]] are equal. Write them together as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\overline{u}=\sum_{m=1}^{\infty }A_{mn}R\left({\beta }_{mn}r\right){e}^{-\left(C_{h}{\beta }_{mn}^2+C_{v}{\alpha }_n^2\right)t}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (27)
|}

where  <math>A_{mn}=A_{1mn}=A_{2mn}</math>;  <math>{\beta }_{mn}={\beta }_{1mn}={\beta }_{2mn}</math>;  <math>\chi \left({\beta }_{mn}\right)={\chi }_1\left({\beta }_{1mn}\right)=</math><math>{\chi }_2\left({\beta }_{2mn}\right)</math>; and <math>P\left({\beta }_{mn}r\right)=P_1\left({\beta }_{1mn}r\right)=</math><math>P_2\left({\beta }_{2mn}r\right)</math>.

After finite Fourier inverse transform, the final general solution without considering disturbance effect can be obtained as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u=\frac{2}{H}\sum_{m=1}^{\infty }\sum_{n=1}^{\infty }A_{mn}R\left({\beta }_{mn}r\right)\sin{\alpha }_nz{e}^{-\left(C_{h}{\beta }_{mn}^2+C_{v}{\alpha }_n^2\right)t}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (28)
|}

which is consistent with of the solution proposed by Gao and Shi [12].

===3.2 Verification by FDM===

In order to further verify the correctness of the proposed solution, the FDM is used for comparison. First, the solution domain  <math display="inline">\left\{\begin{array}{cc}
r_0\leq r\leq r_{e}; & 0\leq z\leq H
\end{array}\right\}</math> is discretized, as follow

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\left\{\begin{array}{lll}
r_i=r_0+i\Delta r & \Delta r=\frac{r_{e}-r_0}{I} & i=0,1...,I\\
z_j=j\Delta z & \Delta z=\frac{H}{J} & j=0,1...,J
\end{array}\right.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (29)
|}

where  <math display="inline">r_i</math> and  <math display="inline">z_j</math> are respectively the radial distance and depth corresponding to point <math display="inline">i</math> or <math display="inline">j</math>;  <math display="inline">\Delta r</math> and  <math display="inline">\Delta z</math> are respectively radial mesh and depth mesh.

The governing equation (i.e., Eq.<span id='cite-ZEqnNum571777'></span>[[#ZEqnNum571777|(1)]]) and solving conditions (i.e., Eqs.<span id='cite-ZEqnNum823648'></span>[[#ZEqnNum823648|(2)]]-<span id='cite-ZEqnNum798795'></span>[[#ZEqnNum798795|(5)]]) are discretized by using the forward-backward implicit difference scheme as shown below

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\frac{u_{i,j}^k-u_{i,j}^{k-1}}{\Delta t}=C_{{h}i}\left[\frac{u_{i-1,j}^{k-1}-2u_{i,j}^{k-1}+u_{i+1,j}^{k-1}}{{\left(\Delta r\right)}^2}+\right. </math><math>\left. \frac{1}{r_i}\frac{u_{i+1,j}^{k-1}-u_{i-1,j}^{k-1}}{2\Delta r}\right]+</math>
|-
| <math>C_{v}\frac{u_{i,j-1}^{k-1}-2u_{i,j}^{k-1}+u_{i,j+1}^{k-1}}{{\left(\Delta z\right)}^2}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (30)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u_{i,j}^0=u_0\left(r_i,z_j\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (31)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c} u_{-1,j}^k=u_{1,j}^k\\ u_{I,j}^k=0 \end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (32)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>k_{{h}I_{d}}\frac{u_{I_{d},j}^k-u_{I_{d}-1,j}^k}{\Delta r}=</math><math>k_{{h}I_{d}+1}\frac{u_{I_{d}+1,j}^k-u_{I_{d},j}^k}{\Delta r}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (33)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
u_{i,0}^k=0\\
u_{i,J-1}^k=u_{i,J+1}^k
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (34)
|}

where  <math>k_{{h}i}=\left\{ \begin{array}{cc}
k_{d}, & 0\leq i\leq I_{d}\\
k_{h}, & I_{d}<i\leq I
\end{array}\right.</math> and  <math>C_{{h}i}=\left\{ \begin{array}{cc}
C_{d}, & 0\leq i\leq I_{d}\\
C_{h}, & I_{d}<i\leq I
\end{array}\right.</math>. The above FDM solution is programmed with MATLAB, and the specific programming steps are as follows:

For  <math display="inline">i\in \left[0,I\right],j=0</math>

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u_{i,0}^k=0</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (35)
|}

For  <math display="inline">i=I,j\in \left[1,J\right]</math>

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u_{I,j}^k=0</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (36)
|}

For  <math display="inline">i=0,j\in \left[1,J\right]</math>

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u_{0,j}^k=\left\{1-2\Delta t\left[\frac{C_{{h}0}}{{\left(\Delta r\right)}^2}+\right. \right. </math><math>\left. \left. \frac{C_{v}}{{\left(\Delta z\right)}^2}\right]\right\}u_{0,j}^{k-1}+</math><math>\frac{2C_{{h}0}\Delta t}{{\left(\Delta r\right)}^2}u_{1,j}^{k-1}+</math><math>\frac{C_{v}\Delta t}{{\left(\Delta z\right)}^2}u_{0,j-1}^{k-1}+</math><math>\frac{C_{v}\Delta t}{{\left(\Delta z\right)}^2}u_{0,j+1}^{k-1}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (37)
|}

For  <math display="inline">i\in \left[1,I-1\right],j=J</math>

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u_{i,J}^k=\frac{C_{{h}i}\Delta t\left(2r_{i,j}-\Delta r\right)}{2{\left(\Delta r\right)}^2r_i}u_{i-1,J}^{k-1}+</math><math>\left\{1-2\Delta t\left[\frac{C_{{h}i}}{{\left(\Delta r\right)}^2}+\right. \right. </math><math>\left. \left. \frac{C_{v}}{{\left(\Delta z\right)}^2}\right]\right\}u_{i,J}^{k-1}+</math><math>\frac{C_{{h}i}\Delta t\left(2r_{i,j}+\Delta r\right)}{2{\left(\Delta r\right)}^2r_i}u_{i+1,J}^{k-1}+</math><math>\frac{2C_{v}\Delta t}{{\left(\Delta z\right)}^2}u_{i,J-1}^{k-1}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (38)
|}

For  <math>i\in \left[1,I-1\right],j\in \left[1,J-1\right]</math>

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u_{i,j}^k=\frac{C_{{h}i}\Delta t\left(2r_{i,j}-\Delta r\right)}{2{\left(\Delta r\right)}^2r_i}u_{i-1,j}^{k-1}+</math><math>\left\{1-2\Delta t\left[\frac{C_{{h}i}}{{\left(\Delta r\right)}^2}+\right. \right. </math><math>\left. \left. \frac{C_{v}}{{\left(\Delta z\right)}^2}\right]\right\}u_{i,j}^{k-1}+</math><math>\frac{C_{{h}i}\Delta t\left(2r_{i,j}+\Delta r\right)}{2{\left(\Delta r\right)}^2r_i}u_{i+1,j}^{k-1}+</math><math>\frac{C_{v}\Delta t}{{\left(\Delta z\right)}^2}u_{i,j-1}^{k-1}+\frac{C_{v}\Delta t}{{\left(\Delta z\right)}^2}u_{i,j+1}^{k-1},\, \begin{array}{c} I=1,2,\cdots I-1\\
j=1,2,\cdots J-1
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (39)
|}

The time step is divided into 90 segments uniformly by the time of each order of magnitude (e.g., the time step is  <math>{10}^1{ s}</math> when  <math>t</math> between  <math>{10}^2{ s}</math> and  <math>{10}^3{ s}</math>). The grid of the model is uniform, and the grid size is  <math display="inline">\Delta r</math> (0.25m) <math display="inline">\times \Delta z</math> (0.25m) ([[#img-2|Figure 2]](a)). Three positions (i.e., Point 1 at  <math>\left(r_0,H\right)</math>, Point 2 at  <math>\left(3r_0,{0}{.8}H\right)</math> and Point 3 at  <math>\left(5r_0,{0}{.6}H\right)</math>) of the excess pore pressure are selected to compare the FDM solution and the analytical solution. As shown in [[#img-2|Figure 2]](b), the proposed solution coincides with the FDM solution almost completely, which indicates the correctness of the proposed solution.

<div id='img-2'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|  style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image135.jpeg|102px]]
|  style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image136.jpeg|center|294px]]
|-
|  style="text-align: center;font-size: 75%;"|(a) 
|  style="text-align: center;font-size: 75%;"|(b) 
|- style="text-align: center; font-size: 75%;"
| colspan="2" style="padding:10px;"| '''Figure 2'''. (a) Mesh of the model. (b) Curves of relation between time and excess pore pressure
|}

==4. Results and Discussion==

For convenience of analysis, the time factor,  <math>T_{v}</math> , and the dimensionless excess pore pressure,  <math>\overline{u}</math> are defined as follows

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>T_{v}=C_{v}t/H^2,\overline{u}=u/u_0\left(r_0,H\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (40)
|}

===4.1 Parameter analysis===

[[#img-3|Figure 3]](a) illustrates the distribution of excess pore pressure along the depth. The excess pore pressure increases with the increase in depth. After consolidation for a period of time (i.e.,  <math>T_{v}={10}^{-3}</math> ), the excess pore pressure at different depths in the disturbed zone was still higher than that in the undisturbed zone (i.e.,  <math>r>r_{d}</math>). [[#img-3|Figure 3]](b) shows that the excess pore pressure increases linearly with depth direction under different <math>r</math>. The growth rate of the excess pore pressure with depth is large in the disturbed zone (i.e.,  <math>r\leq r_{d}</math>), while the growth rate of the excess pore pressure with depth is small in the undisturbed zone (i.e.,  <math>r > r_d</math>).

<div id='img-3'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|  style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image144.jpeg|300px]]
|  style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image145.jpeg|center|300px]]
|-
|  style="text-align: center;font-size: 75%;"|(a) 
|  style="text-align: center;font-size: 75%;"|(b) 
|- style="text-align: center; font-size: 75%;"
| colspan="2" style="padding:10px;"| '''Figure 3'''. (a) Distribution of  <math>\overline{u}</math> along  <math>r</math> direction under different  <math>z</math>. (b) Distribution of  <math>\overline{u}</math> along  <math>z</math> direction under different  <math>r</math>
|}


According to the research works proposed by Li et al.  [9], time-dependent bearing capacity is mainly affected by the dissipation of the excess pore pressure at the pile-soil interface. In order to study the influence of the soil disturbance induced by pile installation on the excess pore pressure at pile-soil interface, this paper introduces the dimensionless average excess pore pressure at pile-soil interface (<math>u_{{ave}}\left(t\right)</math>) for analysis:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>u_{ave}(t)=\frac{\displaystyle\int_{0}^H 2\pi r_0u_1\left(t,r_0,z\right){d}z}{\displaystyle\int_{0}^H 2\pi r_0 u_0\left(r_0,z\right){d}z}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (41)
|}


In [[#img-4|Figure 4]], the influence of permeability coefficient ratio  <math>k_{d}/k_\mbox{h}</math> on the excess pore pressure under deep and shallow foundations, where  <math>k_d/k_{h}</math> reflects the disturbance degree of pile installation. When  <math>k_{d}/k_h=1</math>, this solution can degenerate into the solution of homogeneous soil around a pile [12]. Whether it's shallow foundation (i.e., case of  <math>H=5</math>m) or deep foundation (i.e., case of  <math>H = 20</math>m), the dissipation rate of the excess pore pressure gradually decreases with decreases in  <math>k_{d}/k_{h}</math>. That indicates that  <math>u_{ave}\left(t\right)</math> could be overestimated at the anytime without considering the disturbance effect. The dissipation rate of the average excess pore pressure at pile-soil interface decreases obviously with the increase of  <math>k_{d}/k_{h}</math>, which indicated that soil disturbance has a significant hindering effect on  <math>u_{ave}\left(t\right)</math>. By comparing Case of shallow foundation and Case of deep foundation, it is clearly known that the disturbance effect has a greater influence on deep foundation than on shallow foundation. Moreover, the difference of drainage performance of deep and shallow foundations is gradually obvious with the increase of  <math>k_{d}/k_{h}</math>.

<div id='img-4'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[Image:Review_932178682507-image164.jpeg|300px]] 
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figure 4'''. Influence of permeability coefficient ratio  <math>k_{d}/k_{h}</math> on the excess pore pressure 
|}


Unlike  <math>k_{d}/k_{h}</math> which reflect the degree of disturbance, the radius of the disturbed zone,  <math>r_{d}</math>, reflects the disturbance range. The impact of  <math>r_{d}/r_{0}</math> on the excess pore pressure is introduced in [[#img-5|Figure 5]]. This solution can degenerate into the solution of homogeneous soil around a pile when  <math>r_{d}/r_{0}=1</math> [12]. Additionally, compared with the case without disturbed zone (i.e.,  <math>r_{d}/r_{0}=1</math>), the dissipation rate of  <math>u_{{ave}}\left(t\right)</math> decreases gradually with the increase of  <math>r_{d}</math> in the intermediate and later stages. Moreover, it is found that the consolidation rate of the deep foundation is smaller than that of the shallow foundation with the increase of  <math>r_{d}/r_{0}</math>.

<div id='img-5'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_932178682507-image174.jpeg|300px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figure 5'''. Influence of the disturbance radius ratio  <math>r_{d}/r_{0}</math> on the excess pore pressure
|}

===4.2 Test analysis===

To demonstrate the validity of the proposed solution, a set of published laboratory test data is used for comparison [31]. As shown in [[#img-6|Figure 6]](a), the test box is a cube with a side length of 1000 mm, and the test pile is a PVC pipe with a radius of 30 mm and length of 900 mm (length of 700 mm in the soil). The literature suggests that <math>r_e</math> is equal to 20 times <math>r_0</math> [31]. Because the Piezometer 1-3 (i.e., P1-P3, at distance of <math>r = 90</math> mm) are close to the pile, the measurement data fluctuates greatly, so only the data of P4 ~ P9 are used for analysis. P4-P6 are located at distance of <math>r = 210</math> mm and P7-P9 are located at distance of of <math>r  = 330</math> mm. According to the initial excess pore pressures of P4-P9, the parameters of the initial excess pore pressure is determined as shown in [[#img-6|Figure 6]](b).

<div id="_Hlk76471231" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
| [[Image:Review_932178682507-image176.jpeg|240px]]
| [[Image:Review_932178682507-image177.jpeg|center|300px]]
|}
</div>

<div id="_Hlk58349943" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 6. (a) Schematic diagram of laboratory test for a pipe pile and (b) distribution diagram of initial pore pressure by pile-driving.</div>

The calculated and measured values of the excess pore pressure are shown in Figure 7. It is observed that the calculated value is a good fit with the measured value. That indicates the proposed mathematical model has a good ability to predict the variation trend of the excess pore pressure.

<div id="_Hlk76470940" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
| [[Image:Review_932178682507-image178.jpeg|300px]]
| [[Image:Review_932178682507-image179.jpeg|center|300px]]
|}
</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 7. Data of 6 measured points and corresponding calculated values: (a) P4 ~ P6, and (b) P7 ~ P9.</div>

<span id='_Hlk30159828'></span>

===5. Conclusions===

<span id='_Hlk76423000'></span><span id='_Hlk76422546'></span>In this paper, an analytical solution to consolidation problem with disturbance effect in the soil around a pile is derived and its accuracy is verified by degradation analysis and the FDM. The influence of disturbance effect on the excess pore pressure dissipation of the soil around a pile under long and short foundation is analyzed. It was concluded that the excess pore pressure in the disturbed zone is more difficult to dissipate than that in the undisturbed zone; The dissipation rate of the average excess pore pressure gradually decreased obviously with the increase of permeability coefficient ratio,  <math>k_\mbox{d}/k_\mbox{h}</math> , on the whole dissipation period at pile-soil interface; The impact of radius ratio,  <math>r_\mbox{d}/r_\mbox{0}</math> , on the excess pore pressure is interesting. When  <math>r_\mbox{d}/r_\mbox{0}</math> is small,  <math>r_\mbox{d}/r_\mbox{0}</math> only impedes dissipation rate of excess pore pressure in the early and intermediate dissipation periods. But, with  <math>r_\mbox{d}/r_\mbox{0}</math> increases,  <math>r_\mbox{d}/r_\mbox{0}</math> can impedes dissipation rate of excess pore pressure on the whole dissipation period. Disturbance effect (e.g.,  <math>k_\mbox{d}/k_\mbox{h}</math> and  <math>r_\mbox{d}/r_\mbox{0}</math> ) has obvious negative effect on both shallow and deep foundation, and the negative effect on deep foundations is greater than that on shallow foundation. Finally, the proposed solution is used to predict a test case, which shows that the proposed solution has a good ability to predict the variation trend of the excess pore pressure. However, this model only assumes that the soil is elastic and ignores the effect of the external load, which will be the objective of subsequent research.

===Author Contributions===

Conceptualization: Ping Li; writing—original draft preparation: Ping Li and Zhijian Chen; writing—review and editing: Yi Ding.

<br/>

===Institutional Review Board Statement===

This study did not require ethical approval.<br/>

===Informed Consent Statement===

Not applicable.<br/>

===Data Availability Statement===

All data, models, and code generated or used during the study appear in the submitted article.

<br/>

===Conflict of Interest===

The authors declared that they have no conflicts of interest to this work.<br/>

===References===

[1] K T Tran, M McVay, R Herrera, P Lai. A new method for estimating driven pile static skin friction with instrumentation at the top and bottom of the pile. Soil Dynamics and Earthquake Engineering. 2011;31(9):1285-1295.

[2] M E Mabsout, J L Tassoulas. A finite element model for the simulation of pile driving. International Journal for Numerical Methods in Engineering. 1994;37(2):257-278.

[3] J H Hwang, N Liang, C H Chen. Ground response during pile driving. Journal of Geotechnical and Geoenvironmental Engineering. 2001;127(11):939-949.

[4] M Y Fattah, W H Al-Soudani, M Omar. Estimation of bearing capacity of open-ended model piles in sand. Arabian Journal of Geosciences. 2016;9(3):1-14.

[5] Pantea Azimi, M Karimpour-Fard, N Shariatmadari, C Tsuha. A new approach to estimate the bearing capacity of driven piles. Arabian Journal of Geosciences. 2021;14(13):1-12.

[6] M R Mahmood, S J A Qadir. Effect of particle size distribution of cohesionless soils on the ultimate carrying capacity of open ended pipe piles under different saturation conditions. Arabian Journal of Geosciences. 2018;11(21):1-15.

[7] J Dijkstra, W Broere, O M Heeres. Numerical simulation of pile installation. Computers and Geotechnics. 2011;38(5):612-622.

[8] X T Xu, H L Liu, B M Lehane. Pipe pile installation effects in soft clay. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering. 2006;159(4):285-296.

[9] L Li, J P Li, D A Sun, L X Zhang. Time-dependent bearing capacity of a jacked pile: An analytical approach based on effective stress method. Ocean Engineering. 2017;143:177-185.

[10] M F Randolph, C P Wroth. An analytical solution for the consolidation around a driven pile. International Journal for Numerical Analytical Methods in Geomechanics. 1979;3(3):217-229.

[11] W D Guo. Visco-elastic consolidation subsequent to pile installation. Computers and Geotechnics. 2000;26(2):113-144.

[12] Z K Gao, J Y Shi. Consolidation solution of soil around single-pile after pile sinking. Rock and Soil Mechanics. 2008;29(4):979-982.

[13] S H Wang, P P Ni, Z Chen, G X Mei. Consolidation solution of soil around a permeable pipe pile. Marine Georesources Geotechnology. 2019:1-9.

[14] Z Chen, T Xiao, J X Feng, P P Ni, D Q Chen, G X Mei, Y F Chen. Mathematical characterization of pile-soil interface boundary for consolidation analysis of soil around permeable pipe pile. Canadian Geotechnical Journal. 2020(ja).

[15] M Bozozuk, B H Fellenius, L Samson. Soil disturbance from pile driving in sensitive clay. Canadian Geotechnical Journal. 1978;15(3):346-361.

[16] T Fang, M Huang, K Tang. Cross-section piles in transparent soil under different dimensional conditions subjected to vertical load: an experimental study. Arabian Journal of Geosciences. 2020;13(21):1-8.

[17] K R Massarsch, C Wersäll. Cumulative lateral soil displacement due to pile driving in soft clay. Sound Geotechnical Research to Practice: Honoring Robert D. Holtz II; 2013: 462-479.

[18] Z Y Li, K H Wang, W B Wu, C J Leo, N Wang. Vertical vibration of a large-diameter pipe pile considering the radial inhomogeneity of soil caused by the construction disturbance effect. Computers and Geotechnics. 2017;85:90-102.

[19] Y Zhou, J C Chai. Equivalent ‘smear’effect due to non-uniform consolidation surrounding a PVD. Géotechnique. 2017;67(5):410-419.

[20] R Nazir, H Moayedi, P Subramaniam, S Ghareh. Ground improvement using SPVD and RPE. Arabian Journal of Geosciences. 2017;10(23):1-21.

[21] R D Holtz, M B Jamiolkowski, R Lancellotta, R Pedroni. Prefabricated vertical drains: design and performance; 1991.

[22] J Wang, Y L Yang, H T Fu, Y Q Cai, X Q Hu, X M Lou, Y W Jin. Improving consolidation of dredged slurry by vacuum preloading using prefabricated vertical drains (PVDs) with varying filter pore sizes. Canadian Geotechnical Journal. 2020;57(2):294-303.

[23] Y B Deng, G B Liu, B Indraratna, C Rujikiatkamjorn, K H Xie. Model test and theoretical analysis for soft soil foundations improved by prefabricated vertical drains. International Journal of Geomechanics. 2017;17(1):04016045.

[24] Z Chen, P P Ni, G X Mei, Y F Chen. Semi-analytical solution for consolidation of ground with partially penetrating PVDs under the free-strain condition. Journal of Engineering Mechanics. 2020;147(2):04020148.

[25] Y Tian, W B Wu, G S Jiang, M H El Naggar, G X Mei, P P Ni. Analytical solutions for vacuum preloading consolidation with prefabricated vertical drain based on elliptical cylinder model. Computers Geotechnics. 2019;116:103202.

[26] L O Soderberg. Consolidation theory applied to foundation pile time effects. Géotechnique. 1962;12(3):217-225.

[27] Z L Luo, F H Dong. Statistical investigation of bearing capacity of pile foundation based on Bayesian reliability theory. Advances in Civil Engineering. 2019;2019.

[28] B M Lehane, D R Gill. Displacement fields induced by penetrometer installation in an artificial soil. International Journal of Physical Modelling in Geotechnics. 2004;4(1):25-36.

[29] H de Chaunac, A Holeyman. Numerical analysis of the set-up around the shaft of a closed-ended pile driven in clay. Géotechnique. 2018;68(4):332-344.

[30] Q M Chen, M N Haque, M Abu-Farsakh, B A Fernandez. Field investigation of pile setup in mixed soil. Geotechnical Testing Journal. 2014;37(2):268-281.

[31] P P Ni, S Mangalathu, G X Mei, Y L Zhao. Laboratory investigation of pore pressure dissipation in clay around permeable piles. Canadian Geotechnical Journal. 2018;55(9):1257-1267.