Qiang Luo•Ming Wei•Qingyuan Lu•Tengfei Wang

Abstract Piled embankments have been extensively used for high-speed rail over soft soils because of their effectiveness in minimizing differential settlement and shortening the construction period.Stress concentration ratio,defined as the ratio of vertical stress carried by pile heads(or pile caps if applicable)to that by adjacent soils,is a fundamental parameter in the design of piled embankments.In view of the complicated load transfer mechanism in the framework of embankment system,this paper presents a simplified analytical solution for the stress concentration ratio of rigid pile-supported embankments.In the derivation,the effects of cushion stiffness,pile–soil interaction,and pile penetration behavior are considered and examined.A modified linearly elastic-perfectly plastic model was used to analyze the mechanical response of a rigid pile–soil system.The analytical model was verified against field data and the results of numerical simulations from the literature.According to the proposed method,the skin friction distribution,pile–soil relative displacement,location of neural point,and differential settlement between the pile head(or cap)and adjacent soils can be determined.This work serves as a fast algorithm for initial and reasonable approximation of stress concentration ratio on the design aspects of piled embankments.

Keywords Piled embankments·Pile–soil interaction·Pile penetration·Cushion·Rigid pile·High-speed railway

1 Introduction

High-speed railway network in China is the world’s longest and most extensively used—with a total length of 37,900 km by the end of 2020.Owing to the continuing development of high-speed rail projects in China and other countries,track substructures require careful attention to ensure the safe delivery of passengers[1].In the case of a high-speed rail embankment,two major challenges may arise when weak foundation soil occurs:low strength and high compressibility of weak soils that poses geotechnical problems to infrastructures and control of post-construction settlement[2,3].Embankment reinforced with rigid piles has been universally recognized as an effective solution to minimize the settlement of weak subsoil,which contributes greatly to deformation[4–7].

Typical piled embankment primarily comprises piles,pile caps,an angular platform consisting of gravels and geosynthetics(if applicable),and an overlying embankment fill,as shown in Fig.1.The law of stress redistribution among all these components can be complex owing to the interactions within the earthwork system[8].Extensive investigations have been performed to elucidate the load transfer mechanism of piled embankments(e.g.,[9–11])for design purposes.Some analytical and numerical models[12–20]have been developed based on diversified assumptions to account for the soil arching effect at the base of embankments,classified into categories such as limit equilibrium,frictional,and rigid models.However,some fundamental factors may be overlooked in these models(e.g.,cushion layer,geosynthetic reinforcement,and load transfer efficiency of pile–soil system)when evaluating the mechanical behavior of piled embankments.Hence,some theoretical methods[8,21–23]have been developed to analyze the pile–soil interaction below the base of embankments,which can be divided into two main categories:semi-analytical method and closed-form solution.Although these approaches can be used to reveal the deformation characteristics of piles and the surrounding soil,they are not applicable when considerable loads are anticipated under high fill[24].In addition,the role of the cushion layer,which is critical in interpreting the mechanical behavior of foundation reinforced by rigid piles,is ignored in some of the mentioned methods.

Fig.1 Profile of a rigid pile-supported embankment system for high-speed rail

Stress concentration ratio is taken as a fundamental parameter in design of piled embankments and is a key indicator of the deformation and bearing capacity of underlying reinforced foundation;however,few methods are available for reasonably estimating its value.Combined with recent advances in elucidating the load transfer mechanism of piled embankments[25,26],this paper presents a novel analytical solution for predicting the stress concentration ratio.The main advantages of the proposed method over conventional methodologies include proper consideration of cushion stiffness and pile penetration behavior.Besides,a modified linearly elastic-perfectly plastic model was employed to determine the shear resistance and deformation at the pile–soil interface.Based on the improved interface model,negative skin friction along pile shafts can be captured,and more rational skin friction distribution is obtained considering the variation of soil stress level.The proposed method was verified against the results of numerical analysis and field measurements from two case histories.

2 Skin friction distribution

2.1 Pile–soil interface model

The mobilization of skin friction along the pile shaft primarily depends on the shear characteristics of the pile–soil interface.Numerous investigations have been conducted to obtain the stress–strain relationship for the soil–structure interface,providing theoretical descriptions such as exponential,linearly elastic-perfectly plastic,hyperboloid,and bilinear models.In this study,the linearly elastic-perfectly plastic model was selected as the pile–soil interface model owing to its computational simplicity and reasonable accuracy[8,21,22,27,28].

Within the theoretical framework of the linearly elasticperfectly plastic model,two parameters,i.e.,shear stiffness(Ks)and ultimate shear strength of the soil–pile interface(τu),must be determined.Chen et al.[8]illustrated that both parameters would increase with depth and lateral earth pressure;however,an empirical formula[29](see Eq.(1))has been more widely used in many studies(e.g.,[21,22,30])to calculateKswith the shear modulus of surrounding soil(Gs)being a constant.

whereUz，cis the critical shear displacement;landrpare the values of pile length and radius,respectively;rmis the influential radius of the pile in terms of soil domain.As specified in Fig.2a,Ksis independent of depth.

Fig.2 Linearly elastic-perfectly plastic model for pile–soil interaction in two forms:a constant shear stiffness;b normalized critical shear displacement as a constant

The results of shear tests on the sand–steel interface,performed by Evgin and Fakharian[31],indicate that the normal stress level is slightly related toUz，cbut is significantly related to theKs.It appears more appropriate to adopt a constant critical shear displacement in the formulation of the elastic-perfectly plastic model[32–34].In addition,the sample size was reported to significantly affect the critical shear displacement of a shear stress–displacement curve obtained from direct shear tests[35].Consequently,the concept of normalized shear displacement,Uz/L(Lrepresents the box dimension in the shear direction),was introduced to eliminate the effect of sample size in the determination of a representative critical shear displacement,as shown in Fig.2b.In general,the skin friction between construction material and soil is governed by material type,soil properties,and the roughness[36].The parameters for describing the shear behavior at the soil–structure interface are available in literature.

Considering the discussions above,a load transfer model was established based on normalized shear displacement.Hence,the modified shear stiffnessat the pile–soil interface can be calculated by

whereK0is the lateral earth pressure coefficient;γis the unit weight of adjacent soil;psrepresents the vertical stress applied to the soil mass;caandφacorrespond to the adhesion and friction angle of the soil–pile interface,respectively,which can be estimated using the strength reduction method as follows(in case experimental data are not accessible):

whereRinis the strength reduction factor;csandφsare the effective cohesion component and effective friction angle of the adjacent soil,respectively.In general,the value ofRinis 2/3 for sand–steel interface and 1/2 for the clay–steel interface,and a rough interface results in a largerRin.

2.2 Mobilization of skin friction

For a rigid pile-supported embankment,stress concentration occurs at the pile cap(or pile head)under an embankment load.The pile cap would penetrate upward into the cushion layer in extreme conditions,and a negative skin friction can be anticipated along the upper portion of the pile shaft.Below the neural plane,positive skin friction is mobilized to resist the downward movement of the rigid pile.Assuming that the pile–soil interaction follows the linearly elastic-perfectly plastic model,pile tip penetration occurs once the lower pile–soil interface reaches a critical state.The pile basal resistance is then mobilized to counteract the loads from the upper structure and the positive skin friction.In general,the skin friction along the pile shaft can be divided into three segments,corresponding to I,II,and III in Fig.3a.Within Segments I and III,the shear behavior of the pile–soil interface has reached the plastic stage,while the pile–soil interface undergoes an elastic deformation in Segment II.It was assumed thatvaried linearly with depth in Segment II,which was consistent with the findings of Chen et al.[8].Therefore,the mathematical expression forcan be obtained as

Fig.3 Profiles of relative displacement per unit length(a),relative displacement(b),and skin friction(c)on the pile shaft

wherezis the depth(the elevation of the pile head as the reference level);are the critical relative displacements per unit length related to the soil type and properties of the pile surface,respectively;lis the pile length;l1andl2are the locations of the demarcation point.It is noteworthy thatif the subsoil is homogeneous.

According to the geometrical relationship in Fig.3a,we have

Integrating for the area under the curve in Fig.3a results in an expression for the distribution of pile–soil relative displacement(Uz),as graphically represented in Fig.3b.The pile–soil relative displacement at the pile head(Uz，H)can be inferred as

Combining Eqs.(2),(3),and(5)yields(see Fig.3c)

whereA=ca+K0pstanφaandB=K0γtanφa.

3 Stress concentration ratio

3.1 Element selection from a piled embankment system

It is typical to use geosynthetic basal reinforcement in piled embankments for high-speed rail.Several studies have suggested that geosynthetic reinforcement has limited effect on both the stress concentration ratio and maximum settlement of an embankment when the tensile stiffness of the geosynthetic material is less than 860 kN/m[9,37]or when the pile spacing is sufficiently close[38].However,some researchers[12,13,39]argued that geosynthetic reinforcement is crucial for determining the stress concentration ratio.Given the abovementioned controversy and China’s TB 10001–2016 design code,this paper focuses on the pile–soil interaction in the following sections;therefore,necessary simplifications will be performed regarding soil arching and the effect of geosynthetic reinforcement.

The pile groups are generally arranged in a regular spacing pattern,and the pile spacing is significantly smaller compared with the span of an embankment.The performances of a pile–soil system under an embankment load can be treated as the same,except for those beneath the slope[40,41].Therefore,a cylindrical unit cell was introduced,as shown in Fig.4.The equivalent diameterde=2brepresenting the affected region can be calculated by[42]

wherebis the equivalent radius of a cylindrical unit cell;Sdenotes the pile spacing;cgrepresents the shape factor of 1.05 and 1.13 for triangular and square arrangements,respectively.

The following assumptions and simplifications were made in the derivation of the analytical solution:

(a) A one-dimensional compression problem is considered for the system.

(b) The embankment load is treated as a uniformly distributed vertical load,p.

(c) The pressures acting on the pile cap and soil mass are uniformly distributed(closely spaced plies).

(d) The pile configurations are identical,and a linear elastic and isotropic material is used for the pile.

(e) The penetration resistance at the pile cap and pile tip varies linearly with the penetration level.Nonlinear behavior and punching shear failure of cushion,and retraining effect of geosynthetics reinforcement in high fill conditions are not involved.

(f) The foundation is regarded as a perfectly elastoplastic material,and the self-weight of the rigid pile is ignored.

Based on the assumptions,a relationship can be established amongp,pp,andps(see Fig.4b):

whereppandpsdenote the pressure on the pile cap and soil surface beneath the cushion layer,respectively;mcp=defines the area percentage for the pile cap;AcpandAcsare the cross-sectional areas of the pile cap and projected area of soil,respectively.

3.2 Stress and deformation of pile and soil

As demonstrated in Fig.5,elements of the pile and soil mass having a thickness of dzwere considered for mechanical analysis.The force equilibrium equations yields

Fig.5 Schematic of force equilibrium analysis(a)for a given pile element(b)and soil element(c)under one-dimensional compression

wherearepresents the pile shaft diameter;bdenotes the outer diameter of the soil element;σpandσsare the vertical stresses of the pile and soil,respectively.

Substituting Eq.(8)in Eq.(11)and then integrating Eq.(11)on both sides with respect to the variablezyields

whereσpiandσsiare the vertical stresses of the pile and soil,respectively;C1-C8are integral constants;subscripti=1，2，3，4 corresponds to the depthz∈(0，l1)，(l1，l0)，(l0，l2)，(l2，l),respectively.

Integrating Eq.(12)yields

wherewpuandwpbrepresent the compressive deformation of the pile shaft above and below the neutral plane,respectively;Epis the elastic modulus of the pile material.

Similarly,integrating Eq.(13)yields the expression for the compressive deformation of soil mass above and below the neutral plane,wsuandwsb,respectively:

whereEsis the constrained modulus of adjacent soil.For layered soil,Escan be estimated by[43]

whereEsiandhiare the constrained modulus and thickness of the corresponding soil layer,respectively.

3.3 Boundary and compatibility conditions

Following the previous assumptions,the penetration of the pile head/cap(xpu)and pile tip(xpb)can be obtained as

whereCcandCrare the flexibility factors(m/kPa)for the cushion layer and substratum,respectively.They reflect the effect of the cushion layer on the mechanical performance of the pile–soil system.According to Eqs.(12)and(13),σp(l)andσs(l)can be determined as follows:

The stress boundary conditions for the pile and soil mass at the elevation level of the pile cap(z=0)can be expressed as(referring to Fig.4)

The pile penetration at the elevation of the cushion layer(z=0)is equal to the relative displacement between pile and the adjacent soil ofz=0.Combining Eqs.(7)and(17)yields

The following equations can be obtained as per compatibility conditions for Eqs.(12)and(13):

Given the dimensions and physical properties of the cushion layer,the compressive deformation of the cushion layer is not evaluated in deriving the displacement compatibility equations.However,it is noteworthy that the deformation effect of the cushion layer on the behavior of the pile–soil system can be observed in Eq.(17),in the form of a counter force acting on the pile cap.Based on the displacement compatibility condition(Fig.6),the following can be inferred:

Fig.6 Schematic of displacement compatibility given the pile penetration under applied loads

wheress1andsp1are the settlements of the soil surface and pile cap,respectively;sp2andss2denote the settlements of the pile tip and upper surface of the substratum,respectively.Because no differential settlement occurs at the neural plane,the transformation of Eq.(22)yields

3.4 Solution for stress concentration ratio

In summary,13 variables must be determined,i.e.,integral constantsC1-C8from Eqs.(12)and(13);andl1,l0,l2,ps,andpp.They can be solved analytically based on Eqs.(6),(10),(19)–(21),and(23)because the system is consistent.These 13 equations are highly nonlinear,rendering it challenging to solve using an explicit method.Consequently,reducing the number of unknown variables becomes top priority;in fact,C1-C8,ps,andppcan be expressed byl1,l0,andl2.Subsequently,a numerical technique known as the Newton–Raphson method is used to solve the system of equations.The general procedure is as follows:

i. Combining Eqs.(10),(19),and(20)yields

ii. Substituting Eqs.(12)and(13)in Eqs.(19)and(21),respectively,and then combining them with Eq.(24)yields the expressions forC1-C8in terms ofl1,l0,andl2.

iii. Substituting Eqs.(14a),(15a),and(17a)in Eq.(23a)yields the first equation of system.

iv. Substituting Eqs.(14b),(15b),and(17b)in Eq.(23b)yields the second equation of system.

v. Eq.(6)is the third equation of system;the combination of the three equations is then solved using the commercial software MATHEMATICA.

Stress concentration ratio is defined as follows[9]

Substituting Eq.(24)in Eq.(25)yields

Some observations from Eq.(26)are made:

(a) A small embankment load(p)results inl0+l1→0,consequently,napproaches 1.0.This is consistent with field observations because the differential settlement between the pile and soil is limited under a light load owing to the cushion layer.

(b) When the embankment load is extremely heavy(p→∞),napproaches 1.0.This can be attributed to that the foundation settlement becomes extremely large as the pile group reaches a failure state,and the effect of stiffness difference between the pile and soil on the mechanical response of the pile–soil system is less pronounced.

(c) When the cushion layer is completely flexible(Cc→∞),napproaches 1.0 without pile cap.

4 Case study

Two cases obtained from the literature were used to validate the proposed analytical solution.Case 1 was reported by Zhang et al.[44],and Case 2 was numerically performed by Chen et al.[8].The parameters required in the analysis of the two cases are summarized in Table 1.The load share ratio was used in abovementioned documents,defined as the proportion of total embankment load carried by the piles.A relationship is established between and stress concentration ratio(n)as follows:

Table 1 Parameters of interest in the analysis for Cases 1 and 2

An alternative approach was used for identifying parametersandCr,which are unavailable from the literature.The value of was estimated based on[45]or can be obtained through a back analysis of direct shear test results[31,46–49],ranging from 0.1%to 2%for various soil types.CcandCrwere estimated using the empirical formulae reported by Wu et al.[50]as follows:

wherehcandEcare the thickness and constrained modulus of the cushion layer,respectively;ξis an empirical coefficient(0.8 for soft soil and 1.2 for stiff soil);Aris the cross-sectional area of the pile tip;Erandvrrepresent the constrained modulus and Poisson’s ratio for the stratum,respectively.

4.1 Case 1

Zhang et al.[44]reported a field test performed on piled embankments over silty clay of medium compressibility.The embankment has been reinforced by cement-fly ashgravel(CFG)piles and a geogrid layer over pile caps.Field measurements presented little tension in the geogrid layer,approximately 8–10 kN/m.The geosynthetic reinforcement contributed little to the load redistribution on the pile cap and soil surface beneath the cushion layer.The results of two test sections,G1 and G3,were used in this study for validation.It is noteworthy that all piles were floating and arranged in an equilateral triangular pattern with a circular cap.

Figure 7 demonstrates the comparisons between the computed and measured values of the load share ratio from Sections G1 and G3.Three specifications were involved,namely Britain’s BS8006 design code,Nordic handbook(NGG)[19],and EBGEO[51].It is clear that the proposed method yielded satisfactory predictions, while the remaining methods generally overestimated the value ofRp.The primary causes were that some fundamental factors for the mechanism of load transfer,such as pile penetration behavior and skin friction distribution,have been included in the analytical solution.In addition,different approaches to account for the function of the cushion layer have significantly affected the computed results,at least in the case of a high embankment.Generally,the present method demonstrates promise in determining stress concentration ratio in engineering practice.

Fig.7 Measured values vs.computed values for the load share ratio(R p)in Sections G1(a)and G3(b)

4.2 Case 2

In the second scenario,a numerical simulation[8]was performed on a rigid pile-supported embankment without geosynthetic reinforcement over soft soil based on the commercial software PLAXIS.An axisymmetric model was established to represent the unit cell introduced in Fig.4,i.e.,a cylindrical soil column incorporating a single pile.Other details regarding the development of the numerical model can be found in[8].

The results of numerical simulation,Chen’s method[8],and the analytical solution in this study were compared in terms of the skin friction distribution along the pile shaft(Fig.8).The skin friction varied similarly for all three methods,and the distribution pattern was consistent with the concept of three segments introduced in Fig.3.Compared with Chen’s method,the present method agreed well with the simulation results when the ultimate state reached the pile–soil interface.

Fig.8 Comparison of calculated distributions of skin friction along the pile shaft[8]

Figure 9 demonstrates the profile of relative displacement between the pile and adjacent soil from the analytical solution.The pile penetrations at the head and tip are 9.52 and 14.28 mm,respectively.Based on Eqs.(22)and(23),the differential settlement at the pile cap can be obtained as 23.8 mm,similar to the results of numerical analysis(20.0 mm)and Chen’s method(24.0 mm).The load share ratio(Rp)and elevation level of the neutral plane are presented in Table 2,which allows the present solution to be validated.

Table 2 Comparison between three methods in Case 2 regarding the results of load share ratio and neutral plane location

Fig.9 Relative displacement between surrounding soil and pile varying with depth in Case 2 using the present method

5 Concluding remarks

A simple analytical framework was presented for estimating the stress concentration ratio of piled embankment considering the mechanical response of pile–soil system.Two application scenarios reported in existing literature were used to validate the proposed method.The following conclusions were obtained:

(1) The fast algorithm provided rational estimations for the load share ratio,skin friction distribution,location of neural plane,and differential settlement between the pile cap and the adjacent soils.

(2) The pile–soil interaction,stiffness of cushion layer,and pile penetration behavior were crucial for the identification of stress concentration ratio of piled embankments.

(3) The current method offers benefits such as easy access of parameters with clear physical meanings,affording well understanding of the impact of cushion stiffness on stress concentration ratio,and the consideration of pile–soil interaction.

Although the present algorithm yielded reliable predictions,it should only be treated as a simplified approach and a general guide to the field response.Future work is to deal with arching and an advanced model for geosynthetic reinforcement,and reduce the assumptions and simplifications to handle more complex situations.Despite such limitations,this study serves as a fast algorithm for the determination of stress concentration ratio on the design aspects.

AcknowledgmentsThis work was supported by the National Natural Science Foundation of China(Grant Nos.52078435 and 41901073),Shanghai Key Laboratory of Rail Infrastructure Durability and System Safety(Grant No.R202003),and China Postdoctoral Science Foundation(Grant No.2019M663556).

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License,which permits use,sharing,adaptation,distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s)and the source,provide a link to the Creative Commons licence,and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence,unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use,you will need to obtain permission directly from the copyright holder.To view a copy of this licence,visit http://creativecommons.org/licenses/by/4.0/.