(51 intermediate revisions by 3 users not shown) | |||
Line 13: | Line 13: | ||
<span id='_Hlk82103515'></span><span id='_Hlk82103504'></span>Ping Li<sup>1</sup> [mailto:geopingli@163.com geopingli@163.com] | <span id='_Hlk82103515'></span><span id='_Hlk82103504'></span>Ping Li<sup>1</sup> [mailto:geopingli@163.com geopingli@163.com] | ||
− | Zhijian Chen<sup> | + | Zhijian Chen<sup>2</sup> [mailto:zhijianchen_geo@163.com zhijianchen_geo@163.com] |
− | Yi Ding<sup> | + | Yi Ding<sup>3</sup> [mailto:nsscm_yiding@163.com nsscm_yiding@163.com] |
==Affiliations:== | ==Affiliations:== | ||
Line 21: | Line 21: | ||
<sup>1</sup>College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China. | <sup>1</sup>College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China. | ||
− | <span id='OLE_LINK4'></span><sup> | + | <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/> | <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 [17-20]. 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; | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | (1) The soil around the pipe pile is elastic saturated soil, and the pipe pile is | + | |
(2) The soil around the pile only has disturbance effect in radial direction; | (2) The soil around the pile only has disturbance effect in radial direction; | ||
Line 60: | Line 57: | ||
(3) The soil parameters used in the model remain unchanged during the consolidation process. | (3) The soil parameters used in the model remain unchanged during the consolidation process. | ||
− | <div | + | <div id='img-1'></div> |
− | [[Image:Review_932178682507-image5.png|246px]] | + | {| 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 77: | Line 74: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum571777'></span> <math>\begin{array}{ | + | | <span id='ZEqnNum571777'></span> <math>\begin{array}{c} |
− | \frac | + | \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}\\ |
− | \frac | + | \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> | \end{array}</math> | ||
|} | |} | ||
Line 85: | Line 82: | ||
|} | |} | ||
+ | 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 | |
− | + | ||
− | + | ||
− | + | ||
− | The initial pore pressure distribution of soil around a pile after pile-driving can generally be obtained through numerical simulation or experimental test | + | |
{| class="formulaSCP" style="width: 100%; text-align: center;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 103: | Line 97: | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (2) | | 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 | 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 | ||
Line 112: | Line 105: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum743846'></span> <math> | + | | <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) | | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 127: | Line 119: | ||
| <span id='ZEqnNum201641'></span> <math>\begin{array}{c} | | <span id='ZEqnNum201641'></span> <math>\begin{array}{c} | ||
{u_1\vert }_{r=r_\mbox{d}}={u_2\vert }_{r=r_\mbox{d}}\\ | {u_1\vert }_{r=r_\mbox{d}}={u_2\vert }_{r=r_\mbox{d}}\\ | ||
− | k_\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> | \end{array}</math> | ||
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (4) | | 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 | The top of the ground is permeable, while the bottom is impermeable. Therefore, the vertical boundary can be expressed as | ||
Line 142: | Line 133: | ||
|- | |- | ||
| <span id='ZEqnNum798795'></span> <math>\begin{array}{c} | | <span id='ZEqnNum798795'></span> <math>\begin{array}{c} | ||
− | {u_1\vert }_{z=0}= | + | {u_1\vert }_{z=0}=\displaystyle\frac{\partial u_1}{\partial z}\Bigg|_{z=H}=0\\ |
− | {u_2\vert }_{z=0}= | + | {u_2\vert }_{z=0}=\displaystyle\frac{\partial u_2}{\partial z}\Bigg|_{z=H}=0 |
\end{array}</math> | \end{array}</math> | ||
|} | |} | ||
Line 149: | Line 140: | ||
|} | |} | ||
+ | ===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., Eqs.<span id='cite-ZEqnNum823648'></span>[[#ZEqnNum823648|(2)]]-<span id='cite-ZEqnNum201641'></span>[[#ZEqnNum201641|(4)]]), which yields | |
− | + | ||
− | 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., | + | |
{| class="formulaSCP" style="width: 100%; text-align: center;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 159: | Line 149: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum149568'></span> <math> | + | | <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> |
− | \frac | + | |- |
− | \frac | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (6) | ||
Line 172: | Line 161: | ||
{| style="text-align: center; margin:auto;" | {| 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} | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (7) | ||
Line 182: | Line 171: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum791600'></span> <math> | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (8) | ||
Line 193: | Line 182: | ||
|- | |- | ||
| <span id='ZEqnNum399916'></span> <math>\begin{array}{c} | | <span id='ZEqnNum399916'></span> <math>\begin{array}{c} | ||
− | {{\overline{u}}_1\vert }_{r=r_ | + | {{\overline{u}}_1\vert }_{r=r_{d}}={{\overline{u}}_2\vert }_{r=r_{d}}\\ |
− | k_ | + | 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> | \end{array}</math> | ||
|} | |} | ||
Line 200: | Line 189: | ||
|} | |} | ||
− | + | 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>. | |
− | where <math> | + | |
According to the separation of variables method, <math>{\overline{u}}_1</math> and <math>{\overline{u}}_2</math> can be written as | According to the separation of variables method, <math>{\overline{u}}_1</math> and <math>{\overline{u}}_2</math> can be written as | ||
Line 218: | Line 206: | ||
|} | |} | ||
+ | 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 | |
− | + | ||
− | Substituting Eq. (10) into Eq. (6), one can obtain | + | |
{| class="formulaSCP" style="width: 100%; text-align: center;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 229: | Line 216: | ||
|- | |- | ||
| <span id='ZEqnNum505543'></span> <math>\begin{array}{c} | | <span id='ZEqnNum505543'></span> <math>\begin{array}{c} | ||
− | \frac{\frac{ | + | \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\\ |
− | \frac{\frac{ | + | \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> | \end{array}</math> | ||
|} | |} | ||
Line 236: | Line 223: | ||
|} | |} | ||
+ | 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 | |
− | + | ||
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 246: | Line 232: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum188101'></span> <math>\begin{array}{ | + | | <span id='ZEqnNum188101'></span> <math>\begin{array}{l} |
− | R_{1mn}=a_{1mn} | + | R_{1mn}=a_{1mn}{J}_0\left({\beta }_{1mn}r\right)+b_{1mn}{Y}_0\left({\beta }_{1mn}r\right)\\ |
− | T_{1mn}=c_{1mn} | + | T_{1mn}=c_{1mn}{e}^{-\left(C_{d}{\beta }_{1mn}^2+C_{v}{\alpha }_n^2\right)t} |
\end{array}</math> | \end{array}</math> | ||
|} | |} | ||
Line 259: | Line 245: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum833305'></span> <math>\begin{array}{ | + | | <span id='ZEqnNum833305'></span> <math>\begin{array}{l} |
− | R_{2mn}=a_{2mn} | + | R_{2mn}=a_{2mn}{J}_0\left({\beta }_{2mn}r\right)+b_{2mn}{Y}_0\left({\beta }_{2mn}r\right)\\ |
− | T_{2mn}=c_{2mn} | + | T_{2mn}=c_{2mn}{e}^{-\left(C_{h}{\beta }_{2mn}^2+C_{v}{\alpha }_n^2\right)t} |
\end{array}</math> | \end{array}</math> | ||
|} | |} | ||
Line 267: | Line 253: | ||
|} | |} | ||
+ | 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 | |
− | + | ||
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 278: | Line 263: | ||
|- | |- | ||
| <span id='ZEqnNum354201'></span> <math>\begin{array}{c} | | <span id='ZEqnNum354201'></span> <math>\begin{array}{c} | ||
− | {\overline{u}}_1=\sum_{m=1}^{\infty }\left[A_{1mn} | + | {\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=\sum_{m=1}^{\infty }\left[A_{2mn} | + | {\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> | \end{array}</math> | ||
|} | |} | ||
Line 285: | Line 270: | ||
|} | |} | ||
− | + | Substituting Eq.<span id='cite-ZEqnNum354201'></span>[[#ZEqnNum354201|(14)]] into the boundary condition of Eq.<span id='cite-ZEqnNum791600'></span>[[#ZEqnNum791600|(8)]], it yields | |
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 294: | Line 278: | ||
|- | |- | ||
| <span id='ZEqnNum665918'></span> <math>\begin{array}{c} | | <span id='ZEqnNum665918'></span> <math>\begin{array}{c} | ||
− | {\overline{u}}_1=\sum_{m=1}^{\infty }A_{1mn}R_1\left({\beta }_{1mn}r\right) | + | {\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=\sum_{m=1}^{\infty }A_{2mn}R_2\left({\beta }_{2mn}r\right) | + | {\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> | \end{array}</math> | ||
|} | |} | ||
Line 301: | Line 285: | ||
|} | |} | ||
+ | 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 | |
− | + | ||
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 311: | Line 294: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum164985'></span> <math>R_1\left({\beta }_{1mn}r_ | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (16) | ||
|} | |} | ||
− | |||
and | and | ||
Line 324: | Line 306: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum741519'></span> <math>\ | + | | <span id='ZEqnNum741519'></span> <math>\left\{ \begin{array}{c} |
− | C_ | + | C_{d}{\beta }_{1mn}^2=C_{h}{\beta }_{2mn}^2\\ |
− | \frac{R_2\left({\beta }_{2mn}r_ | + | \displaystyle\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}=\displaystyle\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}</math> | + | \end{array}\right.</math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (17) | | 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 | |
− | + | ||
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 342: | Line 323: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum775505'></span> <math>\begin{array}{ | + | | <span id='ZEqnNum775505'></span> <math>\begin{array}{l} |
− | {\overline{u}}_1=\sum_{m=1}^{\infty }A_{2mn}\frac{R_2\left({\beta }_{2mn}r_ | + | {\overline{u}}_1=\displaystyle\sum_{m=1}^{\infty }A_{2mn}\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){e}^{-\left(C_{d}{\beta }_{1mn}^2+C_{v}{\alpha }_n^2\right)t}\\ |
− | {\overline{u}}_2=\sum_{m=1}^{\infty }A_{2mn}R_2\left({\beta }_{2mn}r\right) | + | {\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> | \end{array}</math> | ||
|} | |} | ||
Line 350: | Line 331: | ||
|} | |} | ||
− | + | <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 | |
− | <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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 359: | Line 339: | ||
|- | |- | ||
| <span id='ZEqnNum753186'></span> <math>\begin{array}{c} | | <span id='ZEqnNum753186'></span> <math>\begin{array}{c} | ||
− | A_{2mn} | + | 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> | \end{array}</math> | ||
|} | |} | ||
Line 366: | Line 346: | ||
|} | |} | ||
− | + | 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 | |
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 374: | Line 353: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>\begin{array}{ | + | | <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> | \end{array}</math> | ||
|} | |} | ||
Line 382: | Line 361: | ||
|} | |} | ||
− | + | 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 | |
− | 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 | + | |
{| class="formulaSCP" style="width: 100%; text-align: center;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 390: | Line 368: | ||
{| style="text-align: center; margin:auto;" | {| 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_ | + | | <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) | | 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>. | |
− | where <math>r_ | + | |
Then, the coefficient <math>A_{2mn}</math> can be specified as | Then, the coefficient <math>A_{2mn}</math> can be specified as | ||
Line 405: | Line 382: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>A_{2mn}=\frac{a_1\left(\ | + | | <math>A_{2mn}=\displaystyle\frac{a_1\left(\sin{\alpha }_nH-\sin{\alpha }_nh_0\right)}{{\alpha }_n{}^2}\displaystyle\frac{\displaystyle\frac{R_2\left({\beta }_{2mn}r_{d}\right)}{R_1\left({\beta }_{1mn}r_{d}\right)}{\delta }_1+{\delta }_2}{\displaystyle\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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (22) | ||
|} | |} | ||
− | |||
where | where | ||
− | <math>{\delta }_1=\frac{1}{{\beta }_{1mn}}\left\{\frac{1}{{\beta }_{1mn}}\left[R_1\left({\beta }_{1mn}r_ | + | <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_ | + | <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_ | + | <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_ | + | <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. | 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;" | |
− | + | |- | |
− | m_ | + | | |
− | r_0=0.25 | + | {| style="text-align: center; margin:auto;width: 100%;" |
− | H=20 | + | |- |
− | \end{array}</math> | + | | style="text-align: center;" | <math>\begin{array}{ccc} |
+ | m_{v}={10}^{-2}{m}^2/{kN} & {\gamma }_{w}=10{kN}/{m}^3 & k_{v}=2\times {10}^{-9}{m}/{s}\\ | ||
+ | r_0=0.25{ m} & r_{p}=1.2\times r_{d} & r_{e}=20\times r_{0}\\ | ||
+ | H=20{ m} & h_0=0.5{ m} & a_1=5{kPa}/{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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 447: | Line 427: | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (23) | | style="width: 5px;text-align: right;white-space: nowrap;" | (23) | ||
|} | |} | ||
− | |||
and | and | ||
Line 461: | Line 440: | ||
|} | |} | ||
− | + | <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 | |
− | <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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 469: | Line 447: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <span id='ZEqnNum350685'></span> <math>\ | + | | <span id='ZEqnNum350685'></span> <math>\left\{ \begin{array}{c} |
R_1\left({\beta }_{1mn}r\right)=R_2\left({\beta }_{2mn}r\right)\\ | 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) | P_1\left({\beta }_{1mn}r\right)=P_2\left({\beta }_{2mn}r\right) | ||
− | \end{array}</math> | + | \end{array}\right.</math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (25) | | 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 | |
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 490: | Line 467: | ||
|} | |} | ||
− | + | 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 | |
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 498: | Line 474: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>\overline{u}=\sum_{m=1}^{\infty }A_{mn}R\left({\beta }_{mn}r\right) | + | | <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) | | 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>. | |
− | 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 | After finite Fourier inverse transform, the final general solution without considering disturbance effect can be obtained as | ||
Line 513: | Line 488: | ||
{| style="text-align: center; margin:auto;" | {| 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 | + | | <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) | | 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 | |
− | + | ||
− | + | ||
− | + | ||
− | r_0\leq r\leq r_ | + | |
\end{array}\right\}</math> is discretized, as follow | \end{array}\right\}</math> is discretized, as follow | ||
Line 534: | Line 506: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>\ | + | | <math>\left\{\begin{array}{lll} |
− | r_i=r_0+i\Delta r & \Delta r=\frac{r_ | + | 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 | z_j=j\Delta z & \Delta z=\frac{H}{J} & j=0,1...,J | ||
− | \end{array}</math> | + | \end{array}\right.</math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (29) | | 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 | |
− | + | ||
− | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 552: | Line 523: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>\frac{u_{i,j}^k-u_{i,j}^{k-1}}{\Delta t}=C_{ | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (30) | ||
Line 572: | Line 545: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>\begin{array}{c} | + | | <math>\begin{array}{c} u_{-1,j}^k=u_{1,j}^k\\ u_{I,j}^k=0 \end{array}</math> |
− | 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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (32) | ||
Line 585: | Line 555: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>k_{ | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (33) | ||
Line 603: | Line 573: | ||
|} | |} | ||
− | + | where <math>k_{{h}i}=\left\{ \begin{array}{cc} | |
− | where <math>k_{ | + | k_{d}, & 0\leq i\leq I_{d}\\ |
− | k_ | + | k_{h}, & I_{d}<i\leq I |
− | k_ | + | \end{array}\right.</math> and <math>C_{{h}i}=\left\{ \begin{array}{cc} |
− | \end{array}</math> and <math>C_{ | + | C_{d}, & 0\leq i\leq I_{d}\\ |
− | C_ | + | C_{h}, & I_{d}<i\leq I |
− | C_ | + | \end{array}\right.</math>. The above FDM solution is programmed with MATLAB, and the specific programming steps are as follows: |
− | \end{array}</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> | For <math display="inline">i\in \left[0,I\right],j=0</math> | ||
Line 623: | Line 592: | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (35) | | style="width: 5px;text-align: right;white-space: nowrap;" | (35) | ||
|} | |} | ||
− | |||
For <math display="inline">i=I,j\in \left[1,J\right]</math> | For <math display="inline">i=I,j\in \left[1,J\right]</math> | ||
Line 636: | Line 604: | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (36) | | style="width: 5px;text-align: right;white-space: nowrap;" | (36) | ||
|} | |} | ||
− | |||
For <math display="inline">i=0,j\in \left[1,J\right]</math> | For <math display="inline">i=0,j\in \left[1,J\right]</math> | ||
Line 645: | Line 612: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>u_{0,j}^k=\left\{1-2\Delta t\left[\frac{C_{ | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (37) | ||
|} | |} | ||
− | |||
For <math display="inline">i\in \left[1,I-1\right],j=J</math> | For <math display="inline">i\in \left[1,I-1\right],j=J</math> | ||
Line 658: | Line 624: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>u_{i,J}^k=\frac{C_{ | + | | <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) | | 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> | |
− | For <math | + | |
{| class="formulaSCP" style="width: 100%; text-align: center;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 671: | Line 636: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>u_{i,j}^k=\frac{C_{ | + | | <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\\ |
− | I=1, | + | j=1,2,\cdots J-1 |
− | j=1, | + | |
\end{array}</math> | \end{array}</math> | ||
|} | |} | ||
Line 679: | Line 643: | ||
|} | |} | ||
+ | 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;" | |
− | + | ||
− | + | ||
− | + | ||
|- | |- | ||
− | | [[Image:Review_932178682507-image135. | + | | style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image135.jpeg|102px]] |
− | | [[Image:Review_932178682507-image136. | + | | 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_ | + | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 705: | Line 666: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>T_ | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (40) | ||
|} | |} | ||
− | ==4.1 Parameter analysis== | + | ===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= | + | <div id='img-3'></div> |
− | + | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | |
− | + | ||
|- | |- | ||
− | | [[Image:Review_932178682507-image144. | + | | style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image144.jpeg|300px]] |
− | | [[Image:Review_932178682507-image145. | + | | 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;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
Line 733: | Line 695: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | <math>u_ | + | | <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) | | style="width: 5px;text-align: right;white-space: nowrap;" | (41) | ||
Line 739: | Line 701: | ||
− | In | + | In [[#img-4|Figure 4]], the influence of permeability coefficient ratio <math>k_{d}/k_{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= | + | <div id='img-4'></div> |
− | [[Image:Review_932178682507-image164. | + | {| 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_ | + | 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= | + | <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='img-6'></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-image176.jpeg|240px]] | |
− | + | | style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image177.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 6'''. (a) Schematic diagram of laboratory test for a pipe pile. (b) Distribution diagram of initial pore pressure by pile-driving | ||
|} | |} | ||
− | |||
− | |||
− | |||
− | 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. | + | The calculated and measured values of the excess pore pressure are shown in [[#img-7|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= | + | <div id='img-7'></div> |
− | + | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | |
− | + | ||
|- | |- | ||
− | | [[Image:Review_932178682507-image178. | + | | style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image178.jpeg|300px]] |
− | | [[Image:Review_932178682507-image179. | + | | style="text-align: center;padding:10px;"| [[Image:Review_932178682507-image179.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 7'''. Data of 6 measured points and corresponding calculated values. (a) P4-P6.(b) P7-P9 | ||
|} | |} | ||
− | |||
− | + | ==5. Conclusions== | |
− | + | ||
− | < | + | 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. | Conceptualization: Ping Li; writing—original draft preparation: Ping Li and Zhijian Chen; writing—review and editing: Yi Ding. | ||
− | + | <!--===Informed Consent Statement=== | |
− | < | + | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | ===Informed Consent Statement=== | + | |
Not applicable.<br/> | Not applicable.<br/> | ||
Line 809: | Line 768: | ||
All data, models, and code generated or used during the study appear in the submitted article. | All data, models, and code generated or used during the study appear in the submitted article. | ||
− | + | ==Conflict of Interest== | |
− | + | ||
− | + | ||
The authors declared that they have no conflicts of interest to this work.<br/> | The authors declared that they have no conflicts of interest to this work.<br/> | ||
+ | <br/>--> | ||
− | + | ==References== | |
+ | <div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;"> | ||
− | [1] K T | + | [1] Tran K.T., McVay M., Herrera R., Lai P. 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, 31(9):1285-1295, 2011. |
− | [2] M E | + | [2] Mabsout M.E., Tassoulas J.L. A finite element model for the simulation of pile driving. International Journal for Numerical Methods in Engineering, 37(2):257-278, 1994. |
− | [3] J H | + | [3] Hwang J.H., Liang N., Chen C.H. Ground response during pile driving. Journal of Geotechnical and Geoenvironmental Engineering, 127(11):939-949, 2001. |
− | [4] M Y | + | [4] Fattah M.Y., Al-Soudani W.H., Omar M. Estimation of bearing capacity of open-ended model piles in sand. Arabian Journal of Geosciences, 9(3):1-14, 2016. |
− | [5] | + | [5] Azimi P., Karimpour-Fard M., Shariatmadari N., Tsuha C. A new approach to estimate the bearing capacity of driven piles. Arabian Journal of Geosciences, 14(13):1-12, 2021. |
− | [6] M R | + | [6] Mahmood M.R., Qadir S.J.A. 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, 11(21):1-15, 2018. |
− | [7] | + | [7] Dijkstra J., Broere W., Heeres O.M. Numerical simulation of pile installation. Computers and Geotechnics, 38(5):612-622, 2011. |
− | [8] X T | + | [8] Xu X.T., Liu H.L., Lehane B.M. Pipe pile installation effects in soft clay. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, 159(4):285-296, 2006. |
− | [9] | + | [9] Li L., Li J.P., Sun D.A., Zhang L.X. Time-dependent bearing capacity of a jacked pile: An analytical approach based on effective stress method. Ocean Engineering, 143:177-185, 2017. |
− | [10] M F | + | [10] Randolph M.F., Wroth C.P. An analytical solution for the consolidation around a driven pile. International Journal for Numerical Analytical Methods in Geomechanics, 3(3):217-229, 1979. |
− | [11] W D | + | [11] Guo W.D. Visco-elastic consolidation subsequent to pile installation. Computers and Geotechnics, 26(2):113-144, 2000. |
− | [12] Z K | + | [12] Gao Z.K., Shi J.Y. Consolidation solution of soil around single-pile after pile sinking. Rock and Soil Mechanics, 29(4):979-982, 2008. |
− | [13] S H | + | [13] Wang S.H., Ni P.P., Chen Z., Mei G.X. Consolidation solution of soil around a permeable pipe pile. Marine Georesources Geotechnology, 20(9):1097-1105, 2020. |
− | [14] | + | [14] Chen Z., Xiao T., Feng J.X., Ni P.P, Chen D.Q., Mei G.X., Chen Y.F. Mathematical characterization of pile-soil interface boundary for consolidation analysis of soil around permeable pipe pile. Canadian Geotechnical Journal, 58(9):1277-1288, 2020. |
− | [15] | + | [15] Bozozuk M., Fellenius B.H., Samson L. Soil disturbance from pile driving in sensitive clay. Canadian Geotechnical Journal, 15(3):346-361, 1978. |
− | [16] | + | [16] Fang T., Huang M., Tang K. Cross-section piles in transparent soil under different dimensional conditions subjected to vertical load: an experimental study. Arabian Journal of Geosciences, 13(21):1-8, 2020. |
− | [17] K R | + | [17] Massarsch K.R., Wersäll C. Cumulative lateral soil displacement due to pile driving in soft clay. Sound Geotechnical Research to Practice: Honoring Robert D. Holtz II, 462-479, 2013. |
− | [18] Z Y | + | [18] Li Z.Y., Wang K.H., Wu W.B., Leo C.J., Wang N. Vertical vibration of a large-diameter pipe pile considering the radial inhomogeneity of soil caused by the construction disturbance effect. Computers and Geotechnics, 85:90-102, 2017. |
− | [19] | + | [19] Zhou Y., Chai J.C. Equivalent ‘smear’ effect due to non-uniform consolidation surrounding a PVD. Géotechnique, 67(5):410-419, 2017. |
− | [20] | + | [20] Nazir R., Moayedi H., Subramaniam P., Ghareh S. Ground improvement using SPVD and RPE. Arabian Journal of Geosciences, 10(23):1-21, 2017. |
− | [21] R D | + | [21] Holtz R.D., Jamiolkowski M.B., Lancellotta R., Pedroni R. Prefabricated vertical drains: design and performance. Construction Industry Research & Information Assoc., Boston, Mass.: Butterworth, 1991. |
− | [22] | + | [22] Wang J., Yang Y.L., Fu H.T., Cai Y.Q., Hu X.Q., Lou X.M., Jin Y.W. Improving consolidation of dredged slurry by vacuum preloading using prefabricated vertical drains (PVDs) with varying filter pore sizes. Canadian Geotechnical Journal, 57(2):294-303, 2020. |
− | [23] Y B | + | [23] Deng Y.B., Liu G.B., Indraratna B., Rujikiatkamjorn C., Xie K.H. Model test and theoretical analysis for soft soil foundations improved by prefabricated vertical drains. International Journal of Geomechanics, 17(1):04016045, 2017. |
− | [24] | + | [24] Chen Z., Ni P.P., Mei G.X., Chen Y.F. Semi-analytical solution for consolidation of ground with partially penetrating PVDs under the free-strain condition. Journal of Engineering Mechanics, 147(2):04020148, 2020. |
− | [25] | + | [25] Tian Y., Wu W.B., Jiang G.S., El Naggar M.H., Mei G.X., Ni P.P. Analytical solutions for vacuum preloading consolidation with prefabricated vertical drain based on elliptical cylinder model. Computers Geotechnics, 116:103202, 2019. |
− | [26] L O | + | [26] Soderberg L.O. Consolidation theory applied to foundation pile time effects. Géotechnique, 12(3):217-225, 1962. |
− | [27] Z L | + | [27] Luo Z.L., Dong F.H. Statistical investigation of bearing capacity of pile foundation based on Bayesian reliability theory. Advances in Civil Engineering, 2019:9858617, 2019. |
− | [28] B M | + | [28] Lehane B.M., Gill D.R. Displacement fields induced by penetrometer installation in an artificial soil. International Journal of Physical Modelling in Geotechnics, 4(1):25-36, 2004. |
− | [29] | + | [29] de Chaunac H., Holeyman A. Numerical analysis of the set-up around the shaft of a closed-ended pile driven in clay. Géotechnique, 68(4):332-344, 2018. |
− | [30] Q M | + | [30] Chen Q.M., Haque M,N., Abu-Farsakh M., Fernandez B.A. Field investigation of pile setup in mixed soil. Geotechnical Testing Journal, 37(2):268-281, 2014. |
− | [31] P P | + | [31] Ni P.P., Mangalathu S., Mei G.X., Zhao Y.L. Laboratory investigation of pore pressure dissipation in clay around permeable piles. Canadian Geotechnical Journal, 55(9):1257-1267, 2018. |
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
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 [17-20]. 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.
As shown in Figure 1, all involved parameters include the depth of foundation, , the radius of pipe pile, , the radius of the disturbed zone, , and the influencing radius of pipe pile, . 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.
![]() |
Figure 1. Schematic diagram of consolidation model for the soil around a pipe pile |
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:
|
(1) |
where and are the excess pore pressure of disturbed zone (i.e., ) and undisturbed zone (i.e., ), respectively; (i.e., ) and (i.e., ) represent the consolidation coefficient in radial direction of undisturbed and disturbed zones, respectively; (i.e., ) represent the consolidation coefficient in vertical direction of the undisturbed zone; and represent the permeability coefficient in radial direction of disturbed and undisturbed zones, respectively; represent the permeability coefficient in vertical direction of the influencing zone; is the volume compression coefficient of soil; is the unit weight of water; and are the radial distance from the center of pile and the depth from the ground surface, respectively; and is consolidation time.
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
|
(2) |
The pile is impermeable and the excess pore pressure outside the influencing radius affected by pipe pile is ignored, it yields
|
(3) |
The interface between the disturbed and undisturbed zones needs to satisfy the flow continuity condition, that is
|
(4) |
The top of the ground is permeable, while the bottom is impermeable. Therefore, the vertical boundary can be expressed as
|
(5) |
According to the boundary condition of Eq.(5), the finite Fourier sine transform is performed on the governing equation (i.e., Eq.(1)) and its remaining solving conditions (i.e., Eqs.(2)-(4)), which yields
|
(6) |
|
(7) |
|
(8) |
|
(9) |
where ; and .
According to the separation of variables method, and can be written as
|
(10) |
where and are functions of ; and are functions of .
Substituting Eq.(10) into Eq.(6), one can obtain
|
(11) |
where ; and and are the separation constants.
It can be found that Eq.(11) is an ordinary differential equation, and the solution to Eq.(11) can be easily obtained as follows
|
(12) |
|
(13) |
where , are the zero-order Bessel functions of the first- and second- kind, respectively.
Combining Eqs.(12) and (13), and can be rewritten as
|
(14) |
Substituting Eq.(14) into the boundary condition of Eq.(8), it yields
|
(15) |
where ; ; ; and .
Substituting Eq.(15) into the continuous boundary condition (i.e., Eq.(9)), it yields
|
(16) |
and
|
(17) |
where ; and .
Combined with Eq.(16), Eq.(15) can be rewritten as
|
(18) |
and can be solved with Eq.(17). Then, using the initial condition of Eq.(7) and orthogonality of eigenfunctions, one can obtain
|
(19) |
According to Eq.(19), the coefficient can be obtained. Applying finite Fourier inverse transform to Eq.(18), the general solutions can be written as
|
(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
|
(21) |
where is the radius of plastic zone; and are constant to be determined. The distribution range of the initial excess pore pressure is .
Then, the coefficient can be specified as
|
(22) |
where
;
;
;
.
The verification and subsequent parameter sensitivity analysis using the parameters defined below. The parameters not mentioned can be obtained in the graphical analysis.
|
When and are equal, it can be seen from Eq.(17) that
|
(23) |
and
|
(24) |
and are obtained by solving Eq.(24). Substituting Eq. (23) and Eq.(24) back into , , and , it yields
|
(25) |
Combining Eq.(16) with Eq.(25), one can obtain
|
(26) |
So, it is found that and in Eq.(15) are equal. Write them together as
|
(27) |
where ; ; ; and .
After finite Fourier inverse transform, the final general solution without considering disturbance effect can be obtained as
|
(28) |
which is consistent with of the solution proposed by Gao and Shi [12].
In order to further verify the correctness of the proposed solution, the FDM is used for comparison. First, the solution domain is discretized, as follow
|
(29) |
where and are respectively the radial distance and depth corresponding to point or ; and are respectively radial mesh and depth mesh.
The governing equation (i.e., Eq.(1)) and solving conditions (i.e., Eqs.(2)-(5)) are discretized by using the forward-backward implicit difference scheme as shown below
|
(30) |
|
(31) |
|
(32) |
|
(33) |
|
(34) |
where and . The above FDM solution is programmed with MATLAB, and the specific programming steps are as follows:
For
|
(35) |
For
|
(36) |
For
|
(37) |
For
|
(38) |
For
|
(39) |
The time step is divided into 90 segments uniformly by the time of each order of magnitude (e.g., the time step is when between and ). The grid of the model is uniform, and the grid size is (0.25m) (0.25m) (Figure 2(a)). Three positions (i.e., Point 1 at , Point 2 at and Point 3 at ) of the excess pore pressure are selected to compare the FDM solution and the analytical solution. As shown in Figure 2(b), the proposed solution coincides with the FDM solution almost completely, which indicates the correctness of the proposed solution.
![]() |
|
(a) | (b) |
Figure 2. (a) Mesh of the model. (b) Curves of relation between time and excess pore pressure |
For convenience of analysis, the time factor, , and the dimensionless excess pore pressure, are defined as follows
|
(40) |
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., ), the excess pore pressure at different depths in the disturbed zone was still higher than that in the undisturbed zone (i.e., ). Figure 3(b) shows that the excess pore pressure increases linearly with depth direction under different . The growth rate of the excess pore pressure with depth is large in the disturbed zone (i.e., ), while the growth rate of the excess pore pressure with depth is small in the undisturbed zone (i.e., ).
![]() |
|
(a) | (b) |
Figure 3. (a) Distribution of along direction under different . (b) Distribution of along direction under different |
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 () for analysis:
|
(41) |
In Figure 4, the influence of permeability coefficient ratio on the excess pore pressure under deep and shallow foundations, where reflects the disturbance degree of pile installation. When , this solution can degenerate into the solution of homogeneous soil around a pile [12]. Whether it's shallow foundation (i.e., case of m) or deep foundation (i.e., case of m), the dissipation rate of the excess pore pressure gradually decreases with decreases in . That indicates that 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 , which indicated that soil disturbance has a significant hindering effect on . 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 .
![]() |
Figure 4. Influence of permeability coefficient ratio on the excess pore pressure |
Unlike which reflect the degree of disturbance, the radius of the disturbed zone, , reflects the disturbance range. The impact of on the excess pore pressure is introduced in Figure 5. This solution can degenerate into the solution of homogeneous soil around a pile when [12]. Additionally, compared with the case without disturbed zone (i.e., ), the dissipation rate of decreases gradually with the increase of 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 .
![]() |
Figure 5. Influence of the disturbance radius ratio on the excess pore pressure |
To demonstrate the validity of the proposed solution, a set of published laboratory test data is used for comparison [31]. As shown in 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 is equal to 20 times [31]. Because the Piezometer 1-3 (i.e., P1-P3, at distance of 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 mm and P7-P9 are located at distance of of mm. According to the initial excess pore pressures of P4-P9, the parameters of the initial excess pore pressure is determined as shown in Figure 6(b).
![]() |
|
(a) | (b) |
Figure 6. (a) Schematic diagram of laboratory test for a pipe pile. (b) Distribution diagram of initial pore pressure by pile-driving |
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.
![]() |
|
(a) | (b) |
Figure 7. Data of 6 measured points and corresponding calculated values. (a) P4-P6.(b) P7-P9 |
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, , on the whole dissipation period at pile-soil interface; The impact of radius ratio, , on the excess pore pressure is interesting. When is small, only impedes dissipation rate of excess pore pressure in the early and intermediate dissipation periods. But, with increases, can impedes dissipation rate of excess pore pressure on the whole dissipation period. Disturbance effect (e.g., and ) 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.
Conceptualization: Ping Li; writing—original draft preparation: Ping Li and Zhijian Chen; writing—review and editing: Yi Ding.
[1] Tran K.T., McVay M., Herrera R., Lai P. 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, 31(9):1285-1295, 2011.
[2] Mabsout M.E., Tassoulas J.L. A finite element model for the simulation of pile driving. International Journal for Numerical Methods in Engineering, 37(2):257-278, 1994.
[3] Hwang J.H., Liang N., Chen C.H. Ground response during pile driving. Journal of Geotechnical and Geoenvironmental Engineering, 127(11):939-949, 2001.
[4] Fattah M.Y., Al-Soudani W.H., Omar M. Estimation of bearing capacity of open-ended model piles in sand. Arabian Journal of Geosciences, 9(3):1-14, 2016.
[5] Azimi P., Karimpour-Fard M., Shariatmadari N., Tsuha C. A new approach to estimate the bearing capacity of driven piles. Arabian Journal of Geosciences, 14(13):1-12, 2021.
[6] Mahmood M.R., Qadir S.J.A. 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, 11(21):1-15, 2018.
[7] Dijkstra J., Broere W., Heeres O.M. Numerical simulation of pile installation. Computers and Geotechnics, 38(5):612-622, 2011.
[8] Xu X.T., Liu H.L., Lehane B.M. Pipe pile installation effects in soft clay. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, 159(4):285-296, 2006.
[9] Li L., Li J.P., Sun D.A., Zhang L.X. Time-dependent bearing capacity of a jacked pile: An analytical approach based on effective stress method. Ocean Engineering, 143:177-185, 2017.
[10] Randolph M.F., Wroth C.P. An analytical solution for the consolidation around a driven pile. International Journal for Numerical Analytical Methods in Geomechanics, 3(3):217-229, 1979.
[11] Guo W.D. Visco-elastic consolidation subsequent to pile installation. Computers and Geotechnics, 26(2):113-144, 2000.
[12] Gao Z.K., Shi J.Y. Consolidation solution of soil around single-pile after pile sinking. Rock and Soil Mechanics, 29(4):979-982, 2008.
[13] Wang S.H., Ni P.P., Chen Z., Mei G.X. Consolidation solution of soil around a permeable pipe pile. Marine Georesources Geotechnology, 20(9):1097-1105, 2020.
[14] Chen Z., Xiao T., Feng J.X., Ni P.P, Chen D.Q., Mei G.X., Chen Y.F. Mathematical characterization of pile-soil interface boundary for consolidation analysis of soil around permeable pipe pile. Canadian Geotechnical Journal, 58(9):1277-1288, 2020.
[15] Bozozuk M., Fellenius B.H., Samson L. Soil disturbance from pile driving in sensitive clay. Canadian Geotechnical Journal, 15(3):346-361, 1978.
[16] Fang T., Huang M., Tang K. Cross-section piles in transparent soil under different dimensional conditions subjected to vertical load: an experimental study. Arabian Journal of Geosciences, 13(21):1-8, 2020.
[17] Massarsch K.R., Wersäll C. Cumulative lateral soil displacement due to pile driving in soft clay. Sound Geotechnical Research to Practice: Honoring Robert D. Holtz II, 462-479, 2013.
[18] Li Z.Y., Wang K.H., Wu W.B., Leo C.J., Wang N. Vertical vibration of a large-diameter pipe pile considering the radial inhomogeneity of soil caused by the construction disturbance effect. Computers and Geotechnics, 85:90-102, 2017.
[19] Zhou Y., Chai J.C. Equivalent ‘smear’ effect due to non-uniform consolidation surrounding a PVD. Géotechnique, 67(5):410-419, 2017.
[20] Nazir R., Moayedi H., Subramaniam P., Ghareh S. Ground improvement using SPVD and RPE. Arabian Journal of Geosciences, 10(23):1-21, 2017.
[21] Holtz R.D., Jamiolkowski M.B., Lancellotta R., Pedroni R. Prefabricated vertical drains: design and performance. Construction Industry Research & Information Assoc., Boston, Mass.: Butterworth, 1991.
[22] Wang J., Yang Y.L., Fu H.T., Cai Y.Q., Hu X.Q., Lou X.M., Jin Y.W. Improving consolidation of dredged slurry by vacuum preloading using prefabricated vertical drains (PVDs) with varying filter pore sizes. Canadian Geotechnical Journal, 57(2):294-303, 2020.
[23] Deng Y.B., Liu G.B., Indraratna B., Rujikiatkamjorn C., Xie K.H. Model test and theoretical analysis for soft soil foundations improved by prefabricated vertical drains. International Journal of Geomechanics, 17(1):04016045, 2017.
[24] Chen Z., Ni P.P., Mei G.X., Chen Y.F. Semi-analytical solution for consolidation of ground with partially penetrating PVDs under the free-strain condition. Journal of Engineering Mechanics, 147(2):04020148, 2020.
[25] Tian Y., Wu W.B., Jiang G.S., El Naggar M.H., Mei G.X., Ni P.P. Analytical solutions for vacuum preloading consolidation with prefabricated vertical drain based on elliptical cylinder model. Computers Geotechnics, 116:103202, 2019.
[26] Soderberg L.O. Consolidation theory applied to foundation pile time effects. Géotechnique, 12(3):217-225, 1962.
[27] Luo Z.L., Dong F.H. Statistical investigation of bearing capacity of pile foundation based on Bayesian reliability theory. Advances in Civil Engineering, 2019:9858617, 2019.
[28] Lehane B.M., Gill D.R. Displacement fields induced by penetrometer installation in an artificial soil. International Journal of Physical Modelling in Geotechnics, 4(1):25-36, 2004.
[29] de Chaunac H., Holeyman A. Numerical analysis of the set-up around the shaft of a closed-ended pile driven in clay. Géotechnique, 68(4):332-344, 2018.
[30] Chen Q.M., Haque M,N., Abu-Farsakh M., Fernandez B.A. Field investigation of pile setup in mixed soil. Geotechnical Testing Journal, 37(2):268-281, 2014.
[31] Ni P.P., Mangalathu S., Mei G.X., Zhao Y.L. Laboratory investigation of pore pressure dissipation in clay around permeable piles. Canadian Geotechnical Journal, 55(9):1257-1267, 2018.
Published on 17/03/22
Accepted on 06/03/22
Submitted on 04/11/21
Volume 38, Issue 1, 2022
DOI: 10.23967/j.rimni.2022.03.008
Licence: CC BY-NC-SA license
Are you one of the authors of this document?