A Semianalytical Solution for Passively Loaded Piles Adjacent to Surcharge Load

Zhang, Hao; Shi, Minglei; Yang, Lei; Guo, Yuancheng

doi:https://doi.org/10.1155/2020/2398389

Advances in Civil Engineering

On this page

Abstract Introduction Conclusions Appendix Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 2398389 | https://doi.org/10.1155/2020/2398389

A Semianalytical Solution for Passively Loaded Piles Adjacent to Surcharge Load

Hao Zhang,¹Minglei Shi,²Lei Yang,³and Yuancheng Guo¹

Academic Editor: Claudio Mazzotti

Received29 Jul 2019

Accepted17 Feb 2020

Published10 Jun 2020

Abstract

Piles adjacent to a surcharge load commonly support not only active loads from superstructures but also the passive loads caused by soil lateral movement. To investigate the influence of passive load and the response along pile shafts of existing actively loaded piles, a load transfer model for analyzing the soil-pile interaction was developed based on plastic deformation theory and the triparameter soil model. An analytical solution for the deformation and internal force of such piles was proposed using the transfer matrix method, in which the transfer matrix coefficients for piles in free, plastic, and elastic zones were analytically obtained by considering the second-order axial force effect caused by lateral loading and soil yielding based on the triparameter soil model. The proposed methodology was validated by comparing its predictions with field measurements and previously published results. A good match between model predictions, field measurements, and previously published results implies that the proposed method can be used to evaluate the response of passive piles adjacent to a surcharge load. Parametric studies were also carried out to investigate the influence of surcharge pressure, soil resistance, and boundary conditions on the behavior of passively loaded piles adjacent to a surcharge load.

1. Introduction

Piles below existing buildings are commonly subjected to additional passive loading from moving soil caused by, for example, the construction of road embankments or fills adjacent to buildings. These piles are commonly classified as passive piles [1, 2]. Passive loading can be significant, particularly for soft soil deposits and large lateral soil displacement which can induce changes in pile bearing and deformation characteristics. Vertically loaded piles originally designed according to the pile bearing capacity and settlement control may become piles controlled by lateral passive loading, which is commonly neglected and often causes engineering accidents, such as the collapse of the Shanghai BaoGang Steel Warehouse caused by a steel slag surcharge in 1998 [3, 4], and the toppling of a nearly finished 13-story residential building in Shanghai due to an excess surface load to the north and an excavation to the south in 2009 [5].

Proper consideration of the reaction of passive piles to surcharge load is the key to avoid such accidents. Numerous numerical and theoretical methods have been developed to analyze passive piles [6–8]. Numerical methods usually utilize the finite-element method [9, 10] or the finite difference method [11], which can consider complex boundary conditions and geological profiles but requires significant computational capacity and longer time. Theoretical methods can be categorized as displacement-based methods and pressure-based methods [12]. The displacement-based methods can better reflect the field behavior of piles due to soil movements by using free-field soil displacements to determine the soil-pile interaction forces [9, 13, 14]. However, a complex numerical simulation or a large number of field data is needed to determine the free-field soil displacements in advance, which increases the difficulty of applying the displacement-based method in engineering practice [1, 15]. The pressure-based methods assume the distribution of passive loading induced by lateral soil displacements based on laboratory tests or theoretical analyses and calculate the deflection and bending moment of the laterally loaded pile [16–19]. The availability of analytical solutions for laterally loaded piles makes the pressure-based methods more practical.

The capacity of the pressure-based methods to reflect the real passive pile behavior is closely related to the load transfer model employed for predicting the soil-pile interactions. Different soil-pile interaction models have been proposed to evaluate the responses of laterally loaded pile [20–24]. For passively loaded piles, Ito et al. proposed a theoretical method for analyzing the growth mechanism of lateral forces acting on stabilizing piles based on plastic deformation theory [25]. A special boundary condition (the lateral pressure between piles was assumed to be the active Earth pressure) was used to derive the passive load, which may result in the lateral soil pressure in the plastic zone different from the actual soil pressure [26]. Shen presented a limit analysis method in which the soil-pile interaction was analyzed using a plane-strain slip-line theory to consider the soil plastic flow around the piles [27]. Ashour et al. extended the strain wedge (SW) model and deduced the mobilized nonuniformly distributed passive load based on soil-pile interaction in an incremental fashion [16]. However, these studies are primarily about stabilizing piles which have no active loading and cannot be applied to the passively loaded piles adjacent to surcharge load.

Therefore, this paper aims to develop a semianalytical solution for passively loaded piles adjacent to surcharge load (i.e., the piles subjected to both active loading from superstructure and passive loading from adjacent surcharge load). Specifically, a load transfer model for soil-pile interactions is proposed using plastic deformation theory and the triparameter soil model. The passive load induced by a neighboring surcharge load was determined in an incremental fashion. Then, the deformation and internal force solutions for piles under combined active and passive loads were derived using the transfer matrix method. The proposed methodology was validated by comparing its predictions with field measurements and previously published results. Finally, parametric studies were carried out to investigate the influence of surcharge pressure, soil resistance, and boundary conditions.

2. New Method for Analysis of Piles under Both Active and Passive Loadings

2.1. Problem Definition

The problem addressed in this study is schematically shown in Figure 1(a) where the pile is subjected to both the active loading from the superstructure and the passive loading from the adjacent surcharge load, such as an embankment. The pile is embedded in a multilayered soil deposit with n sublayers and has length L, diameter d, and bending stiffness EI. The pile has length L_a and L_b above and below the ground surface, respectively. The embedded part is divided into n segments according to the n soil layers. The active loads from the structure include the horizontal force V₀, moment M₀, and axial load N₀ acting on the pile head. A distributed load acting on the pile above the ground surface (called the free zone) is considered to simulate the extra external loading from, for example, wind or wave. The pile is also subjected to a passive load below the ground surface induced by nearby surcharge loading, which will be determined in Section 2.2.

(a)

(b)

Under the combined active and passive loads, the pile deforms continuously and may push the surrounding soil to its limit state. The soil begins to yield at the ground surface and extends downward, although the rest of the soil along the pile is still in an elastic state, and thus the soil state along the pile shaft can be divided into plastic and elastic zones (Figure 1(b)) [22, 28]. The reaction of the soil is denoted by a series of spring-slider elements (which are deactivated in the elastic state) along the shaft and each element is described by an elastic-perfectly plastic curve.

The influence of axial loads on the response of this passively loaded pile is also considered. An empirical axial force distribution on the pile shaft [29–31] is used to simulate the distribution of axial forces along the pile shaft in a multilayered soil deposit,

For the above ground portion of the pile, the axial force is expressed bywhere η (= A_pγ_G) is a coefficient describing the increase in axial force along the pile in the free zone due to the weight of the pile and A_p and γ_p represent the cross-sectional area and the unit weight of the pile, respectively.

For the pile below the ground, the axial force is assumed to vary with depth as follows:where is the axial force at the bottom of the pile section in the (i − 1)th soil layer and is a coefficient dependent on skin friction mobilization in the soil.

2.2. Load Transfer Model of Soil-Pile Interaction

The model for analyzing soil-pile interactions of existing piles adjacent to a surcharge load is developed following the model of Ito et al. [25]. The soil-pile interaction growth mechanism is shown in Figure 2. With increasing adjacent surcharge pressure, the lateral stress P_BB′ acting on surface BB′ in front of pile increases and causes plastic deformation of the soil. For the load transfer between moving soil and piles, the plastic soil is assumed to start from surface BB′ and extends to surface OO′ where the lateral stress may equal the static Earth pressure P₀ in considering long-term soil consolidation [26], as shown in Figure 2. The plastic deformation zone in the soil continues to develop with increasing lateral stress P_BB′ acting on surface BB′ (Figures 2(a) and 2(c)). When the lateral stress acting on surface BB′ increases sufficiently, soil is plastically deformed at the front and back of the piles (Figure 2(d)).

(a)

(b)

(c)

(d)

In order to analyze the load transfer between the moving soil and pile, the following assumptions are made (Figure 2(d)):(a)The soil layer operates under plane-strain conditions in regard to depth. When the soil moves laterally, two sliding surfaces occur along lines BEDG and B′E′D′G′, in which lines EB and E′B′ make an angle π/4 + φ/2 with the x-axis, and lines DG and D′G′ make an angle π/4 − φ/2 with the x-axis.(b)Soil in BEDGG′D′E′B′ becomes plastic and obeys the Mohr-Coulomb yield criterion. The soil layer is represented as a plastic solid with an angle of internal friction φ and cohesion c.(c)Even if frictional forces act on surfaces BEDG and B′E′D′G′, the stress distribution in soil BEDGG′D′E′B′ is almost the same as if no frictional forces on those surfaces exist.(d)Piles are rigid.

The lateral stress at any point in BB′E′E, EE′A′A, AA′D′D, and DD′G′G can be obtained from “Appendix A” and expressed as follows:

In the BB′E′E zone,

In the EE′A′A zone,

In the AA′D′D zone,

In the DD′G′G zone,where , , , and are the lateral stress at any point x in the BB′E′E zone, EE′A′A zone, AA′D′D zone, and DD′G′G zone, respectively; D₁ is pile spacing; D₂ is net spacing between piles; ; ; ; ; ; ; ; ; ; .

is the total soil pressure acting on surface BB′ which is given bywhere is the additional lateral stress acting on surface BB′ induced by the surcharge load and can be estimated using Boussinesq’s solution and the Flamant integral solution according to the distribution characteristics of surcharge load [26, 32]; K₀ is the static Earth pressure coefficient.

is the lateral stress acting on the surface AA′ which is given by (Appendix A)

The passive load P_s (z) acting on piles induced by moving soil can be obtained by the following steps:(a)Using (7), the total soil pressure P_BB′ (z) acting on surface BB′ can be obtained. Then, the lateral stress P_AA′ (z) acting on surface AA′ can be obtained using (8).(b)When P_AA′ (z) ≤ P₀ = K₀γz, substituting into (3), x = x_o can be obtained. If x_o is in zone BB′E′E (), as given by the condition in (3), the passive load acting on piles can be calculated by

Otherwise, x = x_o can be obtained by substituting into (4). The passive load acting on piles can be calculated by(c)When P_AA′ (z) > P₀ = K₀γz, substituting into (5), it gives x = x_o ()). If x_o is in zone AA′D′D (), as given by the condition in (5), the passive load acting on piles can be calculated by

Otherwise, x = x_o can be obtained by substituting into (6). If , as given by the condition in (6), the passive load acting on piles can be calculated by

If , the passive load acting on piles can be calculated by

The soil resistance behind the piles can be represented by springs for evaluating pile response, which is simulated using a perfectly elastic-plastic model. The reaction of the soil behind piles can be described aswhere k is the soil reaction modulus. k can be described using a triparameter soil model , in which m and n are soil reaction modulus coefficients; z₀ is the initial height of the ground surface; b₁ is the calculated pile width; p_u is ultimate soil resistance; y is pile deflection; is ultimate soil displacement. According to Matlock, for clay and for sand. ε_c is the soil strain, and ε_c = 0.005–0.02 [20].

δ is a correction coefficient for considering the influence of plastic soil deformation on the reaction calculated width of the soil behind the piles. When plastic soil deformation enters the DD′G′G zone, the expansion of plastic deformation may decrease the reaction calculated width of the soil behind the piles (Figures 2(c) and 2(d)):

2.3. Governing Differential Equations

Assuming the pile behaves as a Euler-Bernoulli beam, the differential equations governing pile deflection in the free, plastic, and elastic zones (Figure 1(b)) can be obtained by discretizing each pile segment into subslices. The accuracy of the solution is improved with finer pile slices.

For the pile in free zone above the ground surface, each pile slice has an individual coordinate system marked with the subscript a (Figure 3(a)). The pile length L_a is spilt into n_a subslices with length l_a = L_a/n_a. The distributed lateral load acting on the κth pile slice can be calculated as the average of the load acting on the slice:

(a)

(b)

(c)

The average axial force of the κth slice for the pile in the free zone can be calculated using equation (1):

For the pile beneath the ground surface, the pile length L_b is divided into n segments by the soil layers. Each pile segment has its individual coordinate system marked with the subscript i (Figures 3(b) and 3(c)) and is split into n_i subslices. The average axial force on the pile slice can be calculated using equation (2):where the subscripts i and refer to the pile slice of the ith pile segment in the ith soil layer.

Based on equations (9)∼(13), the average passive load () acting on the pile slice of the ith pile segment in the ith soil layer can be calculated as the average load acting on the corresponding slice as follows:

In the plastic zone, the average ultimate soil resistance behind the pile slice of the ith pile segment, , can be regarded as the average ultimate soil resistance of the top and bottom of the corresponding slice:

Similarly in the elastic zone, linear elastoplastic spring models are used, and the average soil reaction modulus () behind the pile slice of the ith pile segment can be calculated as the average soil reaction modulus of the top and bottom of the corresponding slice as

For the κth pile slice in the free zone (Figure 3(a)), the governing differential equation is

For the pile slice of the ith pile segment in the plastic zone (Figure 3(b)), the differential equation iswhere the subscript p indicates the plastic state.

When the soil resistance is in the elastic stage (Figure 3(c)), the differential equation can be expressed aswhere the subscript e indicates the plastic state.

2.4. Analysis for the Free Zone

Solving (22), the deflection of the ith pile slice is given bywhere is the pile deflection on the κth pile slice at depth ; ; ; C_o1, C_o2, C_o3, and C_o4 are integration constants.

Using equation (25), pile slope , bending moment , and shear force at depth are given by

The transfer matrix function of the ith pile slice in the free zone can be expressed aswhere ; , , , and are the pile deflection, slope, bending moment, and shear force at the bottom of the κth pile slice in the free zone, respectively; ; , , , and are the pile deflection, slope, bending moment, and shear force at the top of the κth pile slice, respectively.

is the transfer matrix coefficient of the κth pile slice:where , , , , and (j = 1–4) are 20 dimensionless coefficients, which are a function of .

As a result, based on the deflection and stress continuity of the pile, the transfer matrix function for the pile segment in the free zone can be expressed aswhere and are 5-by-1 transposed matrices and represent pile deflection, slope, bending moment, and shear force for the head and bottom of the pile segment in the free zone, respectively; is the overall transfer matrix coefficient of the pile in the free zone.

2.5. Analysis for the Plastic Zone

In the plastic soil zone, as the passive load and ultimate soil resistance can be determined using equations (19) and (20), respectively; the form of equation (23) is similar to equation (22), and the solution methods for the two differential equations are similar. The transfer matrix function of equation (23) of the pile slice of the ith pile segment can be expressed aswhere ; , , , and are the pile deflection, slope, bending moment, and shear force at the bottom of the pile slice of the ith pile segment; ; , , , and are the pile deflection, slope, bending moment, and shear force at the head of pile slice.

is the transfer matrix coefficient of the pile slice of the ith pile segment when the pile section is in the plastic zone. When ,where , , , , and (j = 1–4) are 20 dimensionless coefficients, which are function of .

Assuming the pile in the plastic zone extends from the sth slice of the tth pile segment to the uth pile slice of the pile segment, the matrix transfer function for the pile in the plastic zone is given bywhere is the overall transfer matrix coefficient of the pile in the plastic zone.

2.6. Analysis of the Elastic Zone

Similarly for piles in the elastic soil zone, the deflection of the pile slice of the ith pile segment in the elastic zone can be obtained by solving equation (24).

When or ,

When or ,where ; ; . ; ; ; ; , , , and are integration constants.

By substituting equation (33) or (34) into (26), the transfer matrix function of equation (24) for the pile slice of the ith pile segment can be expressed aswhere .

is the transfer matrix coefficient of the pile slice of the ith pile segment when the pile section is in the elastic zone:where , , , , and (j = 1–4) are 20 dimensionless coefficients, which are a function of .

Assuming that the pile range in the elastic zone is from the sth slice of the tth pile segment to the uth pile slice of the pile segment, the matrix transfer function for the pile in the elastic zone can be expressed aswhere is the overall transfer matrix coefficient of the pile in the elastic zone.

2.7. General Solutions

By combing equations (27), (30), and (35) with the transfer matrix approach [33–35], the transfer matrix function for passively loaded piles adjacent to an embankment can be expressed aswhere and are pile deflection vectors, slope, bending moment, and shear force for the pile tip and pile head, respectively. is the transfer matrix coefficient, which is calculated by equations (31) and (36) for piles in the plastic zone and elastic zone, respectively. is the overall transfer matrix coefficient of the pile in the free zone as deduced in equation (29). is defined as the overall pile transfer matrix coefficient, which is the product of the transfer matrix coefficient of each pile slice in the sequence from pile head to pile toe.

Equation (38) is composed of four equations with eight variables including the pile deflection, slope, bending moment, and shear force at the pile head and pile toe. The boundary conditions at the pile head and pile toe must be specified to solve this equation. Five types of boundary conditions are considered in this study: free-head, slip-head, fixed-head, free-base, and fixed-base.

For the free-head case, the boundary conditions at the pile head are

For the slip-head case, the boundary conditions at the pile head are

For the fixed-head case, the boundary conditions at the pile head are

For the free-base case, the boundary conditions at the pile head are

For the fixed-base case, the boundary conditions at the pile head are

Substituting one of the pile head boundary conditions and one of the pile base boundary conditions into equation (38), the remaining pile head variables can be calculated. Thus, based on the response of the pile head , the solution for the pile in the free zone at any depth can be calculated by

Similarly, the pile response at any depth beneath the ground surface can be calculated by

The algorithm for calculating and is summarized in Figure 4, in which the pile deflection convergence criteria ε = 0.003 of the (j − 1)th and jth iteration guarantee iteration accuracy.

3. Validation

In this section, the proposed analytical method is validated by analyzing two cases, one from existing literature and a new case.

3.1. Case from [9]

Based on the displacement-based method, a simplified numerical procedure for analyzing the pile response to lateral soil displacement was described in [9], in which free-field displacements were directly used to solve for the soil-pile load transfer mechanism. This method was applied and verified using several full-scale and centrifuge model tests. Figure 5 shows the calculation model [9], in which the pile embedded in layered clay-sand was analyzed. Clay and sand physical properties were summarized in [9]. The pile was assumed to have been installed prior to construction of the embankment and was connected to a rigid cap such that the pile head could be deflected without rotating (φ_c0 = 0°; V_c0 = 0 kN). The pile has a length of L = 24 m, diameter d = 1.0 m, spacing s = 2.0 m, and modulus E = 28 GPa. The soil resistance was assumed to increase linearly with depth (n = 1; z₀ = 4.5 m), and the coefficients of the soil reaction modulus (m = 2 MN/m⁴ for clay and 30 MN/m⁴ for the sand) were evaluated using Skempton’s empirical method [36].

Figure 6 shows the bending moment and pile deflection predicted by the numerical results [9] and the proposed method. The variation of bending moment and pile deflection with pile depth by the two methods agreed well. There is an underprediction of the pile deflection of the proposed method in comparison to the numerical results [9]. The discrepancy may be because of the different soil models used by the two methods. The pile-soil interaction was modeled by hyperbolic soil springs and the soil resistance was assumed constant with depth in [9], while elastic-perfectly plastic soil model was used and the soil resistance was assumed to increase linearly with depth in this paper. The soil resistance of the proposed method is larger than [9] under 9 m, which result in smaller pile deflection of the proposed method.

(a)

(b)

3.2. Bridge Pile Foundation near an Embankment

The method described above is also employed to analyze the response of offset piles columns in the SX Bridge at Lianyungang, China. SX Bridge is one of the three bridges of the XinG Road interchange, a semicloverleaf semidirectional compounded interchange (Figure 7(a)). XinG Road passes beneath the SX Bridge between the #26 and #28 piers (Figure 7(b)). XinG Road has a width of 42 m and a slope of 3%. The road shoulder height adjacent to the #26 pier is 1.1 m, and the road shoulder height adjacent to the #28 pier is 2.3 m. The #26 and #28 piers are all composed of 8 sets of piles columns. The length and diameter of piles are 62 m and 1.6 m, respectively. The length and diameter of columns are 8 m and 1.2 m, respectively. The space between the #28-1, #28-2, #28-3, and #28-4 piles and the space between the #26-1, #26-2, #26-3, and #26-4 piles are all 6.06 m. The space between the #28-5, #28-6, #28-7, and #28-8 piles and the space between the #26-5, #26-6, #26-7, and #26-8 are all 4.85 m. The distance from the #26 pier (#26-1∼#26-8) to the XinG Road slope is 3.5 m. The distance from the #28 piers to the XinG Road slope increased from 2.37 m at #28-1 to 4.20 m at #28-8 (Figure 8).

(a)

(b)

The 16 columns of #26 and #28 piers were offset away from the embankment after the filling of the XinG was completed. The lateral displacement of column heads in the #26 pier range from 3.2 to 6.6 cm and the lateral displacement of column heads in the #28 piers are all greater than 10.0 cm. The maximum lateral displacement is approximately 13.3 cm located at the head of the 2th column of the #28 pier (Figure 8).

The average XinG Road embankment height is 1.7 m; however, due to the presence of a 1.0 m–1.5 m side trench outside the #28 pier, the effective embankment height in the calculation is 3.0 m. The embankment fill unit weight is 17 kN/m³. Soil mechanical indicators are based on the data from the investigation site and associated triaxial tests (Table 1). The soft soil resistance coefficient is recommended as 0.55 MN/m⁴ based on the resistance coefficient reference value of the rock and soil provided by CRSSDIG [37]. In addition, sliding rubber bearings are used on the top of the columns, and the forces on the column heads are N₀ = 5379 kN, M₀ = 0 kN·m, and V₀ = 120 kN. The piles and columns are all composed of reinforce concrete, with a unit weight and modulus of 25 kN/m³ and 28.5 GPa, respectively.

Figure 9 shows the calculated column and pile lateral displacement compared to the measured data. It is noted that the calculated results of #28 piers are mainly smaller than the measured results, while the calculated results of #26 piers are mainly larger than the measured results. The reason is that the average embankment height is used in the proposed method, which is lower than the actual height at #28 piers and higher than the actual heights at #26 piers (Figure 7). Figure 10 shows the calculated bending moment along the column and pile for #28-2 and #28-6 from the proposed method and the improved strain wedge (SW) method [19]. The predictions of the improved SW method [19] are obtained from back analysis of measured pile head displacements which reflected real field behavior; however, it does not consider the effect of soil arching between piles and the rotation of pile head. This is the reason for the underprediction of pile bending moment of the improved SW method [19] in comparison to the proposed method.

(a)

(b)

3.3. Factors Influencing Pile Response

Taking #28-2 as an example, the effect of embankment pressure, soil resistance, and column boundary conditions on pile-column response is calculated. The influence of embankment pressure on the performance of neighboring piles columns is investigated by varying the embankment height between 1.0 m, 2.0 m, and 3.0 m. The other parameters are held constant as in previous sections. Increasing embankment height (or embankment pressure) causes an increase in bending moment and deflection of neighboring piles columns (Figure 11). The maximum pile bending moment increased from 3333.7 kN·m for H = 1.0 m to 5254.8 kN·m for H = 3.0 m, a 57.6% increase. The column head displacement varied from 5.92 cm for H = 1.0 m to 13.85 cm for H = 3.0 m, an increase of 1.34 times. With increasing embankment height (or embankment pressure), more passive loads are transferred from moving soil to the piles, and the surcharge loading of the soft soil has a significant impact on nearby piles.

(a)

(b)

The influence of soil resistance on pile-column performance is investigated by varying the mud clay reaction modulus coefficient between 2 MN/m⁴, 3 MN/m⁴, and 4 MN/m⁴. Figure 12 shows the effect of the clay reaction modulus coefficient on pile-column response. With increasing clay reaction modulus coefficient, the pile-column bending moment and deflection decrease significantly (Figure 12). The maximum pile bending moment and column head displacement are reduced by 30.2% and 40.8% from m = 2 MN/m⁴ to m = 3 MN/m⁴. The reductions increase to 46.2% and 61.7% for m = 4 MN/m⁴. These results clearly show that soil resistance has a significant impact on pile response; hence, the geological conditions around the piles should be carefully investigated, and reinforcement or protective measures should be employed in soft soils to ensure the stability and safety of piles adjacent to the embankment.

(a)

(b)

The influence of boundary conditions is investigated by comparing the free-head case (M_c0 = 0, V_c0 = 0), slip-head case (φ_c0 = 0, V_c0 = 0), and fixed-head case (y_c0 = 0, φ_c0 = 0). Figure 13 shows the effect of boundary conditions on the piles-columns responses. Bending moment and deflection along the pile and column for the fixed-head case are significantly different from the other two cases. By comparing this case to the free-head and slip-head cases, there is a negative moment in the upper part of the pile, and the maximum pile displacement occurs at 1/2 of the soft clay thickness. Comparing the free-head and slip-head cases, the difference in bending moment primarily occurs in the upper 2/3 s of the soft clay, and the difference in pile deflection primarily occurs within 7d depth beneath the pile head.

(a)

(b)

4. Conclusions

Based on plastic deformation theory, a load transfer model was presented to analyze the soil-pile interaction of passive piles adjacent to a surcharge load. A semianalytical solution was developed to assess the response of piles under active and passive loads using the transfer matrix approach.

The proposed method was validated against the displacement-based method proposed by Goh et al. [9] and field measurements at SX Bridge, Lianyungang, China. Results showed that the proposed method can provide good response predictions (i.e., deflection, bending moment, and shear forces) for piles adjacent to a surcharge load.

A parametric study was conducted to investigate the influence of the neighboring surcharge pressure, soil resistance, and boundary conditions of pile/column heads on piles adjacent to a surcharge load. Results indicate that the maximum bending moment and deflection along piles increase with increasing neighboring surcharge pressure and decreasing soil resistance. The boundary conditions of pile/column heads have a significant influence of the response of piles/columns.

Appendix

A. Lateral Stress in Soils Surrounding Piles

The equilibrium of forces acting on a small soil element in zones BB′E′E, EE′A′A, AA′D′D, and DD′G′G is, respectively, considered in the x-direction (Figure 14).

(a)

(b)

(c)

(d)

In the BB′E′E zone,

In the EE′A′A zone,

In the AA′D′D zone,

In the DD′G′G zone,where ; ; ; ; ; ; ; .

Equations (A.1)∼(A.4) are differential equations with separated variables and can be integrated as follows:

When ,

When ,where , , , and are integration constants.

Assuming the additional lateral stress acting on surface BB′ induced by the surcharge load is which can be determined from Boussinesq’s solution and the Flamant integral solution according to the surcharge load distribution characteristics [26, 32], the total soil pressure acting on surface BB′, , can be obtained using (7).

The constant in (A.5) can be obtained considering (see (7)) as the boundary condition (where ):

Using (A.5) and (A.9), the lateral stress at any point x in the BB′E′E can be expressed using (3).

In regard to the EE′A′A zone, considering (3) for as the boundary condition (surface EE′), the constant in (A.6) can be obtained:

Using (A.6) and (A.10), the lateral stress at any point x in EE′A′A can be expressed using (4).

Substituting x = 0 into (4), the lateral stress acting on the surface AA′, , can be obtained from (8).

In regard to the AA′D′D zone, considering as the boundary condition of the surface AA′, the constant in (A.7) can be obtained:

Using (A.7) and (A.11), the lateral stress at any point x in the AA′D′D can be expressed using (5).

In regard to the DD′G′G zone, considering (5) for as the boundary condition (surface DD′), the constant in (A.8) can be obtained:

Using (A.8) and (A.12), the lateral stress at any point x in the DD′G′G can be expressed using (6).

Notations

, , :	Calculation coefficients
A_p:	Cross-sectional area of the pile
, , :	Calculation coefficients
c:	Cohesion of the soil
, , :	Calculation coefficients
d:	Pile diameter
D₁:	Pile spacing
D₂:	Net spacing between piles
, , :	Calculation coefficients
EI:	Bending stiffness; the pile has length
, , :	Calculation coefficients
k:	Reaction modulus of soil mass
:	Average value of soil reaction behind the pile in the pile slice of the ith pile segment below the ground
K₀:	Coefficient of the static Earth pressure
L:	Pile length
L_a:	Pile length above the ground surface
L_b:	Pile length below the ground surface
l_a:	Length of each pile slice in free zone
M₀:	Moment acting on the pile head
, ,:	Moment for pile in the free zone, plastic zone, and elastic zone
N₀:	Axial load acting on the pile head
:	Average axial force of the κth slice in free zone
:	Average axial force for the pile in pile slice of the ith pile segment below the ground
p:	Soil resistance
p_u:	Ultimate soil resistance
:	Average value of ultimate soil resistance for the pile in the pile slice of the ith pile segment
P₀:	Static Earth pressure
P_s (z):	Distributed passive load acting on the pile below ground surface
:	Average value of passive load acting on the pile in the pile slice of the ith pile segment
:	Average distributed force acting on the κth pile slice in free zone
:	Overall transfer matrix coefficient of the entire pile
:	Overall transfer matrix coefficient of the pile in the free zone
:	Transfer matrix coefficient of the κth pile slice in the free zone
:	Overall transfer matrix coefficient of the pile in the plastic zone
:	Transfer matrix coefficient of the pile slice of the ith pile segment in the plastic zone
:	Overall transfer matrix coefficient of the pile in the elastic zone
:	Transfer matrix coefficient of the pile slice of the ith pile segment in the elastic zone
, , :	Matrix variables
V₀:	Horizontal force acting on the pile head
, , :	Shear force for pile in the free zone, plastic zone, and elastic zone
y:	Pile deflection
, , :	Deflection for the pile in the free zone, plastic zone, and elastic zone
z:	Pile depth
, ,:	Pile depth for the pile in the free zone, elastic zone, and plastic zone
γ, γ_p:	Unit weight of the soil and pile
φ:	Internal friction of the soil
, , :	Slope for the pile in the free zone, plastic zone, and elastic zone
, , :	Calculation coefficients
η:	Increase coefficient of the axial force in the free zone
δ:	Correction coefficient of soil reaction.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work wass supported by the National Natural Science Foundation of China (Grant no. 51608490) and the Key Scientific Research Project of Colleges and Universities in Henan Province (Grant no. 16A560009), for which the authors are grateful.

References

H. G. Poulos, “Ground movements-a hidden source of loading on deep foundations,” DFI Journal-The Journal of the Deep Foundations Institute, vol. 1, no. 1, pp. 37–53, 2007.
View at: Publisher Site | Google Scholar
E. Beer, “The effects of horizontal loads on piles, due to surcharge or seismic effects,” in Proceedings of 9th International Conference on Soil Mechanics and Foundation Engineering, Tokyo, Japan, July 1977.
View at: Google Scholar
M. Yang, B. T. Zhu, and F. Q. Chen, “Pilot study on collapse of an industrial building due to adjacent surcharge loads,” Chinese Journal of Geotechnical Engineering, vol. 24, no. 4, pp. 446–450, 2002.
View at: Google Scholar
B. T. Zhu and M. Yang, “Collapse of a piled industrial building due to surcharge loads,” in Proceedings of the 6th International Conference on Tall Buildings (ICTBVI), Singapore, December 2005.
View at: Google Scholar
W. D. Wang, Q. Li, Y. Hu, J. W. Shi et al., “Field investigation of collapse of a 13-story high-rise residential building in Shanghai,” Journal of Performance of Constructed Facilities, vol. 31, no. 4, Article ID 04017012, 2017.
View at: Publisher Site | Google Scholar
E. A. Ellis and S. M. Springman, “Modelling of soil-structure interaction for a piled bridge abutment in plane strain FEM analyses,” Computers and Geotechnics, vol. 28, no. 2, pp. 79–98, 2001.
View at: Publisher Site | Google Scholar
M. R. Karim, S.-C. R. Lo, and C. T. Gnanendran, “Behaviour of piles subjected to passive loading due to embankment construction,” Canadian Geotechnical Journal, vol. 51, no. 3, pp. 303–310, 2014.
View at: Publisher Site | Google Scholar
J. L. Zheng, J. H. Zhang, and X. F. Li, “Effects of embankment load on piled bridge abutments on soft clay,” China Journal of Highway and Transport, vol. 26, no. 2, pp. 48–55, 2013.
View at: Google Scholar
A. T. C. Goh, C. I. Teh, and K. S. Wong, “Analysis of piles subjected to embankment induced lateral soil movements,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 123, no. 9, pp. 792–801, 1997.
View at: Publisher Site | Google Scholar
D. P. Stewart, R. J. Jewell, and M. F. Randolph, “Numerical modelling of piled bridge abutments on soft ground,” Computers and Geotechnics, vol. 15, no. 1, pp. 21–46, 1993.
View at: Publisher Site | Google Scholar
H. G. Poulos and E. H. Davis, Pile Foundation Analysis and Design, John Wiley and Sons, New York, NY, USA, 1980.
D. P. Stewart, R. J. Jewell, and M. F. Randolph, “Design of piled bridge abutments on soft clay for loading from lateral soil movements,” Géotechnique, vol. 44, no. 2, pp. 277–296, 1994.
View at: Publisher Site | Google Scholar
Y. K. Chow, “Analysis of piles used for slope stabilization,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 20, no. 9, pp. 635–646, 1996.
View at: Publisher Site | Google Scholar
F. Y. Liang, F. Yu, J. B. Li et al., “Analysis of bearing capacity of a single pile under adjacent building subjected to lateral soil movements,” Rock and Soil Mechanics, vol. 31, no. 2, pp. 449–454, 2010.
View at: Google Scholar
M. K. Kelesoglu and S. F. Cinicioglu, “Free-field measurements to disclose lateral reaction mechanism of piles subjected to soil movements,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 136, no. 2, pp. 331–343, 2010.
View at: Publisher Site | Google Scholar
M. Ashour and H. Ardalan, “Analysis of pile stabilized slopes based on soil-pile interaction,” Computers and Geotechnics, vol. 39, pp. 85–97, 2012.
View at: Publisher Site | Google Scholar
W. D. Guo, “Elastic models for nonlinear response of rigid passive piles,” International Journal for Numeriacl and Analytical Methods in Geomechanics, vol. 38, no. 18, pp. 1969–1989, 2014.
View at: Publisher Site | Google Scholar
S. Muraro, A. Madaschi, and A. Gajo, “On the reliability of 3D numerical analyses on passive piles used for slope stabilisation in frictional soils,” Géotechnique, vol. 64, no. 6, pp. 486–492, 2014.
View at: Publisher Site | Google Scholar
H. Zhang, M. L. Shi, W. S. Hu et al., “Analysis of sloping pier nearby embankment of overpass bridge in interchange,” Journal of Southeast University (Natural Science Edition), vol. 43, no. 3, pp. 617–623, 2013.
View at: Google Scholar
H. Matlock, “Correlations for design of laterally loaded piles in soft clay,” in Proceedings of the 2nd Annual Offshore Technology Conference, Houston, TX, USA, April 1970.
View at: Google Scholar
Y.-m. Hsiung, “Theoretical elastic-plastic solution for laterally loaded piles,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 129, no. 5, pp. 475–480, 2003.
View at: Publisher Site | Google Scholar
W. D. Guo, “Simple model for nonlinear response of 52 laterally loaded piles,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 139, no. 2, pp. 234–252, 2012.
View at: Google Scholar
L. Zhang, “Nonlinear analysis of laterally loaded rigid piles in cohesionless soil,” Computers and Geotechnics, vol. 36, no. 5, pp. 718–724, 2009.
View at: Publisher Site | Google Scholar
L. Zhang and S. Ahmari, “Nonlinear analysis of laterally loaded rigid piles in cohesive soil,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 37, no. 2, pp. 201–220, 2013.
View at: Publisher Site | Google Scholar
T. Ito and T. Matsui, “Methods to estimate lateral force acting on stabilizing piles,” Soils and Foundations, vol. 15, no. 4, pp. 43–59, 1975.
View at: Publisher Site | Google Scholar
M. X. Zhu, W. M. Gong, G. P. Xu et al., “Analysis of bearing mechanism of axially loaded isolation piles under large-scale surcharge loading,” Chinese Journal of Rock Mechanics and Engineering, vol. 33, no. 2, pp. 421–432, 2014.
View at: Google Scholar
Z. J. Shen, “Stabilizing force and design method for anti-slide piles,” Chinese Journal of Geotechnical Engineering, vol. 14, no. 1, pp. 51–56, 1992.
View at: Google Scholar
W. D. Guo, “On limiting force profile, slip depth and response of lateral piles,” Computers and Geotechnics, vol. 33, no. 1, pp. 47–67, 2006.
View at: Publisher Site | Google Scholar
L. Zhang, X. N. Gong, and J. L. Yu, “Analysis of deformation and internal force for single pile under combined vertical and lateral loads considering yielding of soil,” Rock and Soil Mechanics, vol. 32, no. 8, pp. 2441–2446, 2011.
View at: Google Scholar
M. H. Zhao, L. G. Wu, and J. H. Liu, “Finite difference solution for loadbearing and anti-slide pile with consideration of P-Δ effect,” Engineering Mechanics, vol. 25, no. 3, pp. 102–106, 2008.
View at: Google Scholar
M. Zhu, Y. Zhang, W. Gong et al., “Generalized solutions for axially and laterally loaded piles in multilayered soil deposits with transfer matrix method,” International Journal of Geomechanic, vol. 17, no. 4, Article ID 04016104, 2016.
View at: Google Scholar
C.-D. Wang, “Lateral stress caused by horizontal and vertical surcharge strip loads on a cross-anisotropic backfill,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 29, no. 14, pp. 1341–1361, 2005.
View at: Publisher Site | Google Scholar
T. Nogami and S. K. Paulson, “Transfer matrix approach for nonlinear pile group response analysis,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 9, no. 4, pp. 299–316, 1985.
View at: Publisher Site | Google Scholar
Q. T. Liu and G. R. Ni, Transfer Matrix Approach for Structure Analysis, Railway Publishing House Co., Ltd., Beijing, China, 1997.
L. Zhang, M. H. Zhao, and X. J. Zou, “Behavior of laterally loaded piles in multilayered soils,” International Journal of Geomechanic, vol. 15, no. 2, Article ID 06014017, 2015.
View at: Publisher Site | Google Scholar
A. W. Skempton, “The bearing capacity of clays,” in Proceedings of the Building Research Congress, London, England, January 1951.
View at: Google Scholar
A. J. Zhang, “Model of interaction between pile and soil & simplified calculating method of pile used to stability unstable slope,” South China University of Technology, Guangzhou, China, 2013, PhD. thesis.
View at: Google Scholar

Copyright

Copyright © 2020 Hao Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

427

Downloads

975

Citations