[BACK]
 Computers, Materials & ContinuaDOI:10.32604/cmc.2022.021899 Article

Numerical Analysis of Laterally Loaded Long Piles in Cohesionless Soil

1Physics and Engineering Mathematics Department, Faculty of Engineering-Mattaria, Helwan University, Cairo, 11718, Egypt
2Faculty of Engineering, King Salman International University, South Sinai, Ras Sedr, Egypt
3Basic and Applied Science Department, College of Engineering and Technology, Arab Academy for Science, Technology and Maritime Transport, Cairo, 11799, Egypt
*Corresponding Author: Ayman Abd-Elhamed. Email: moh_fathy79_6@aast.edu
Received: 19 July 2021; Accepted: 30 August 2021

Abstract: The capability of piles to withstand horizontal loads is a major design issue. The current research work aims to investigate numerically the responses of laterally loaded piles at working load employing the concept of a beam-on-Winkler-foundation model. The governing differential equation for a laterally loaded pile on elastic subgrade is derived. Based on Legendre-Galerkin method and Runge-Kutta formulas of order four and five, the flexural equation of long piles embedded in homogeneous sandy soils with modulus of subgrade reaction linearly variable with depth is solved for both free- and fixed-headed piles. Mathematica, as one of the world's leading computational software, was employed for the implementation of solutions. The proposed numerical techniques provide the responses for the entire pile length under the applied lateral load. The utilized numerical approaches are validated against experimental and analytical results of previously published works showing a more accurate estimation of the response of laterally loaded piles. Therefore, the proposed approaches can maintain both mathematical simplicity and comparable accuracy with the experimental results.

Keywords: Numerical solution; laterally loaded pile; cohesionless soil; Legendre-Galerkin; Runge-Kutta

1  Introduction

Pile foundations are frequently used, especially in weaker soils, to support various structures subjected to lateral loads such as high-rise buildings, communication towers, wind turbines, earth-retaining structures, bridges, tanks and offshore structures. Lateral loads owing to wind, wave, dredging, traffic and seismic events are considered significant on these structures since they are eventually transmitted to the piles [1,2]. As a result, the piles have been analyzed by considering a concentrated force and/or moment acting on the top of the pile. Over many decades, several methods have been proposed for designing and analysis of piles subjected to lateral loads including the subgrade reaction approach [3], the p-y approach [4,5], the finite element approach [6], the finite difference approach [7], and the analytical method [8]. Of these approaches, elastic solutions based on beam-on-Winkler-foundation model, albeit approximate, are probably the most widely used in engineering practice due to their simplicity, as well as they provide satisfactory results. Winkler's model is a particularly attractive approach used to reliably capture the soil-pile interaction. In this model the pile is simulated as a flexural beam connected to a series of narrowly spaced independent and continuous Winkler springs and dashpots distributed along the pile shaft. To investigate the mechanical behavior of laterally loaded piles in clay, several numerical investigations using the LPILE software were performed by Moayed et al. [9]. Moreover, Chang [10] derived an analytic solution to get the responses of laterally loaded long piles in cohesive soil considering constant subgrade reaction modulus. Furthermore, Different techniques have been proposed to analyze the behaviour of piles frequently subjected to lateral loads in sandy soil with different boundary conditions at their ends, including power series solution [11], finite element method [12], and finite difference method [3,13]. Further full-scale tests were performed to explore the behavior of laterally loaded piles either in sand [14] or in clay [15]. Numerical simulations considering both theoretical predictions and experimental validations are the common powerful tools of analysis in the field of geotechnical engineering particularly complex applications.

Recently, many researchers used Legendre polynomials in different methods to construct various mathematical models. These methods can solve Lane-Emden type of differential equation [16], differential equation with second and fourth order [17], the equation of Cahn-Hilliard [18], integral equation of Fredholm type [19], Helmholtz equation [20], Volterra integral equations in the second kind [21], integral-differential of Fredholm type in linear form [22] and Abels integral equation [23]. Based on Legendre-Galerkin method, the pile flexural equation can be written as Ay(η)=f, where A is a differential operator. The solution of our problem can be approximated in Legendre series as y=j=0ncjPj(η). By applying the Galerkin method to minimize the residual yields j=0ncjAPj(η)f,Pr(η)=0. The differential equation is converted into discrete linear system. This system is solved, and the coefficients are determined for our approximate solution of the differential equation. In addition to the above method, the Bogacki–Shampine method [24] is a powerful numerical solution used for solving ordinary differential equations of the investigated problem. It is a Runge–Kutta method of two successive orders, for example 4th and 5th, with multi stages. An adaptive step size is implemented based on the error estimated between the solutions of the two successive orders. Elbashbeshy et al. [25,26] used this method in flow and heat transfer problems of fluid/nanofluid over different stretching surfaces and found that it is accurate by comparing it with methods used to solve the same problems.

This paper aims to present simple numerical methods to capture the behaviour of single piles under lateral working loads. First, mathematical formulation of the problem and the derivation of the governing differential equations are presented. Thereafter, two techniques are developed for numerical approximations of the derived equations. Then the effect of boundary conditions of the pile head on the behaviors of piles is studied as well. In order to examine the validity of the proposed numerical techniques, the obtained numerical results for deflection and bending moment along the laterally loaded piles are compared with results of previous studies. The proposed numerical techniques in the current study provide a better approach for structural designers to simply solve for the displacement and bending moment responses of laterally loaded piles. Consequently, the techniques can be easily applied in practice as an alternative approach to analyze and design laterally loaded long piles. It is worth noting that, the proposed solutions are applicable only for homogeneous soil with reaction modulus that varies linearly with depth. So, further investigation and analyses are needed for multi-layered soil with the modulus of subgrade reaction varies with any functions with the depth.

2  Mathematical Formulation

The model under examination consists of a pile that is assumed to be perfectly glued to the surrounding soils suggesting that there is no relative movement between the soil and the pile [27]. It is worth noting that the horizontal and vertical loads are rarely coupled. Thereupon, it is common among structural designers to neglect the vertical load effect during handling the laterally loaded pile problem. Accordingly, the pile under consideration is subjected to applied lateral force H0 at distance e above the soil surface, as shown in Fig. 1. Furthermore, the thickness of soil strata, L, coincides with the embedment length of pile. The pile elastic modulus and moment of inertia are Ep and Ip, respectively. The soil is characterized by a subgrade modulus k which linearly increases with the depth as a good approximation for granular soils [28]. For a laterally loaded pile on elastic homogeneous foundation, the governing differential equation employing the subgrade reaction theory can be derived by considering the equilibrium of transverse forces over a differential segment of the pile, as shown in Fig. 2.

Figure 1: Illustration of pile–soil system using Winkler model

Figure 2: Forces acting in differential pile element

Assuming all forces acting horizontally and applying equilibrium conditions leads to

V(V+dVdzdz)k(z)y(z)dz=0(1)

According to that,

dVdz+k(z)y(z)=0(2)

Summation of moments about axis at end of element yields

M(M+dMdzdz)+Vdz=0(3)

The relation between bending moment and the shear force is given by

V=dMdz(4)

Differentiating Eq. (4) with respect to z and substituting it into Eq. (2) provide

d2Mdz2+k(z)y(z)=0(5)

The basic moment-curvature relationship of elementary pile can be written as

M=EpIpd2ydz2(6)

The differential equation of a laterally loaded pile can take the form

EpIpd4ydz4+k(z)y(z)=0,(7)

where y(z) is the horizontal pile deflection at depth z and EpIp is the flexural stiffness of pile.

Assuming a linear variation of soil stiffness with soil depth, k can be written as follows [29]

k=nhz,(8)

where here nh is the coefficient of subgrade reaction that represents the rate of increase of subgrade reaction modulus with depth z.

Eqs. (7) and (8) yield

EpIpd4ydz4+nhzy(z)=0.(9)

The solution of Eq. (9) requires an iterative procedure to achieve convergence of the relative stiffness factor T as follows [30]

T=(EpIpnh)0.2.(10)

Defining a dimensionless variable x=z/T and introducing it into Eq. (9) yields

EpIpnhT5d4ydx4+xy(x)=0.(11)

By substituting Eq. (10) into Eq. (11), the governing equation for the soil-pile system is rewritten the form

d4ydx4+xy(x)=0.(12)

3  Boundary Conditions

As shown in Fig. 3, the solution of Eq. (12) requires boundary conditions to balance the number of equations and the number of unknowns.

Figure 3: Boundary conditions for the pile under consideration

Boundary conditions of infinitely long piles fixed at the bottom can be written as

y(LT)=0,y(LT)=0,y(LT)=0,y(LT)=0.(13)

The boundary conditions at the top of the pile depend on the circumstances of the lateral deflection, slope, bending moment and shear force. These are generalized into the following two categories.

The pile head is not restrained against rotation and translation. Therefore, the boundary conditions at the pile head can take the form

y(0)=HOT3EpIp,y(0)=MOT2EpIp.(14)

The pile head is completely restrained against rotation. The boundary conditions at the pile head are given by

y(0)=HOT3EpIp,y(0)=0.(15)

4  Legendre-Galerkin Method

4.1 Preliminaries

Orthogonal systems play a vital role in performing mathematical analysis. This can be due to functions belonging to very general classes can be expanded in series of orthogonal functions e.g., Fourier-Bessel series, Fourier series, etc. Orthogonal polynomials of degree n, Pn(x);n=0,1,2,, represent an important class of orthogonal systems where many of the special functions encountered in the applications, such as Jacobi, Legendre, Laguerre, Hermite and Chebyshev polynomials, are part of that class. The Legendre polynomial is defined as [31].

Pn(x)=2nk=0n(nk)(n+k12n)xk,(16)

with the orthogonality on the interval [1,1], i.e.,

11Pn(x)Pm(x)dx={22n+1,ifm=n0,ifmn.(17)

Lemma [32]: If the arbitrary constants q,k,m and j are positive integers then

Pk(q)(x)=i=0kqCq(k,i)Pi(x),(k+1)iseven(18)

and

xmPj(x)=n=02mΘm,n(j)Pj+mn(x),Pr(x)=0forr1,(19)

where

Cq(k,i)=2q1(2i+1)Γ(q+ki2)Γ(q+ki+12)Γ(q)Γ(2q+ki2)Γ(3q+ki2),(20)

Θm,n(j)=(1)n2j+mnm!(2j+2m2n+1)Γ(j+1)Γ(j+mn+1)×k=max(0,jn)j+mnΓ(j+k+1)2k(n+kj)!Γ(3j+2m2nk+2),×=0jk(1)Γ(2j+mnk+1)Γ(j+m+n+1)!(jk)!Γ(j+1)Γ(k++1)×2F1(jkn,j+m+;3j+2m2nk+2;2),(21)

and

2F1(N,a;2a;2)={0,forNodd(12)N2(a+12)N2,forNeven.

Theorem: If the arbitrary constants k, l, q and m are positive integers, then

11Pk(q)(x)Pl(x)dx={2i=0kqCq(k,i)2i+1,ifl=i,(k+1)iseven0,ifli(22)

and

11xmPk(x)Pl(x)dx={2n=02mΘm,n(j)2(j+mn)+1ifl=j+mn0,iflj+mn,(23)

where Cq(k,i) and Θm,n(j) are defined in Eqs. (20) and (21), respectively.

Proof. By using the given lemma and Eq. (17), the theorem is proved.

4.2 Legendre-Galerkin Method

To solve Eq. (12), the solution domain of our problem must be changed from the interval [0,L] where L to the interval [1,1] using the linear transformation

x=L2(η+1).(24)

So, Eq. (12) could be rewritten as

y(η)+(2L)3(η+1)y(η)=0,(25)

subject to the conditions of free-head pile

y(1)=(2L)3HT3EpIp,y(1)=(2L)2MT2EpIp,y(1)=y(1)=y(1)=y(1)=0,(26)

Moreover, the conditions of fixed-head pile are:

y(1)=(2L)3HT3EpIp,y(1)=y(1)=y(1)=y(1)=y(1)=0,(27)

The solution of Eq. (25) subject to the conditions given in Eqs. (26) or (27) is written in finite expansion of shifted Legendre function as

yn(η)j=0ncjPj(η).(28)

Reducing Eq. (25) is performed by orthogonalizing the residual with respect to the basic functions as follows

y(η),Pr(η)+(2L)3(η+1)y(η),Pr(η)=0,(29)

where the inner product , is defined as

ζ(η),χ(η)=11ζ(η)χ(η)dη.(30)

By substituting Eq. (28) into Eq. (29), a discrete system of n+1 unknowns is given in the form

Ac=b,(31)

whereA=[a0,0a0,1a0,3a0,n1a0,na1,0a1,1a1,3a1,n1a2,na2,0a2,1a2,3a2,n1a2,nP0(1)P1(1)P2(1)Pn1(1)Pn(1)P0(1)P1(1)P2(1)Pn1(1)Pn(1)P0(1)P1(1)P2(1)Pn1(1)Pn(1)P0(1)P1(1)P2(1)Pn1(1)Pn(1)P0(1)P1(1)P2(1)Pn1(1)Pn(1)P0(1)P1(1)P2(1)Pn1(1)Pn(1)],

c=[c0c1c3cn1cn]T and aj,r=Pj(η),Pr(η)+(2L)3(η+1)Pj(η),Pr(η). All values of aj,r could be calculated using the given theorem. The vector b based on the boundary conditions as follows.

In this case, the vector b is written as

b=[000(2L)3HT3EpIp(2L)2MT2EpIp0000]T.(32)

The linear system in Eq. (31) is solved using one of the available numerical techniques. By obtaining the coefficients values of {cj}0n, the approximate solution given in Eq. (28) becomes

y(η)=3.99HEI2/5nh3/5(1.9×106+8.02×105η8.4×104η2+2.9×103η3+9.6×103η40.109η5+0.284η60.012η71.13η8+2.06η93.17η10+5.04η11+6.35η1235.29η13+24.47η14+52.42η1574.02η1624.81η17+82.27η189.96η1947.01η20+16.04η21+13.97η226.69η231.73η24+η25)+2.58MEI2/5nh3/5(9.44×107+3.615η8.314η2+6.1213η30.010η40.079η5+0.376η60.373η70.591η8+1.02η92.26η10+11.16η119.75η1233.95η13+59.65η14+23.97η15103.32η16+18.82η17+88.61η1840.16η1941.67η20+26.40η21+10.35η228.13η231.06η24+η25).(33)

By using the inverse linear transformation η=(2L)x1andx=zT, the solutions of the governing differential equation for a laterally loaded pile on elastic subgrade, Eq. (12), are given for deflection and moment as

y(z)=2.429Hnh3/5(EI50.67EI24/5znh1/5+0.07EI22/5z3nh3/53.05×107EI21/5z4nh4/58.33×103EI4z5nh+1.84×103EI19/5z6nh6/5+1.7×105EI18/5z7nh7/55.7×105EI17/5z8nh8/5+1.04×105EI16/5z9nh9/52.99×106EI3z10nh2+1.28×106EI14/5z11nh11/53.73×107EI13/5z12nh12/5+7.28×108EI12/5z13nh13/51.02×108EI11/5z14nh14/5+1.07×109EI2z15nh15/58.76×1011EI9/5z16nh16/5+5.62×1012EI8/5z17nh17/52.84×1013EI7/5z18nh18/5+1.13×1014EI6/5z19nh19/53.52×1016EIz20nh20/5+8.35×1018EI4/5z21nh21/51.46×1019EI3/5z22nh22/5+1.79×1021EI2/5z23nh23/51.36×1023EI1/5z24nh24/5+4.85×1026z25nh5)+1.619Mnh3/5(EI51.08nh1/5EI24/5z+0.31nh2/5EI23/5z2+1.1×106nh4/5EI21/5z48.34×103nhEI4z5+3.01×103nh6/5EI19/5z63.8×104nh7/5EI18/5z7+1.2×105nh8/5EI17/5z84.8×106nh9/5EI16/5z9+2.6×106nh2EI3z103.97×107nh11/5EI14/5z112.81×108nh12/5EI13/5z12+2.15×108nh13/5EI12/5z133.97×107nh11/5EI14/5z112.81×108nh12/5EI13/5z12+2.15×108nh13/5EI12/5z134.47×109nh14/5EI11/5z14+5.82×1010nh3EI2z155.44×1011nh16/5EI9/5z16+3.84×1012nh17/5EI8/5z172.09×1013nh18/5EI7/5z18+8.86×1015nh19/5EI6/5z192.89×1016nh4EIz20+7.17×1018nh21/5EI4/5z211.31×1019nh22/5EI3/5z22+1.65×1021nh23/5EI2/5z231.29×1023nh24/5EI1/5z24+4.71×1026nh5z25)(34)

and

M(z)=H(EI27/5z8.8×106nh1/5EI26/5z20.41nh2/5EI5z3+0.134nh3/5EI24/5z4+1.8×103nh4/5EI23/5z57.8×103nhEI22/5z6+1.8×103nh6/5EI21/5z76.5×104nh7/5EI4z8+3.4×104nh8/5EI19/5z91.19×104nh9/5EI18/5z10+2.7×105nh2EI17/5z114.5×106nh11/5EI16/5z12+5.48×107nh12/5EI3z135.11×108nh13/5EI14/5z14+3.71×109nh14/5EI13/5z152.11×1010nh3EI12/5z16+9.41×1012nh16/5EI11/5z173.24×1013nh17/5EI2z18+8.52×1015nh18/5EI9/5z191.64×1016nh19/5EI8/5z20+2.20×1018nh4EI7/5z211.82×1020nh21/5EI6/5z22+7.07×1023nh22/5EIz23)+Mnh1/5(EI28/5+2.3×105nh2/5EI26/5z20.27nh3/5EI5z3+0.15nh4/5EI24/5z42.64×102nhEI23/5z5+1.1×103nh6/5EI22/5z65.6×104nh7/5EI21/5z7+3.8×104nh8/5EI4z87×105nh9/5EI19/5z96.01×106nh2EI18/5z10+5.4×106nh11/5EI17/5z111.3×106nh12/5EI16/5z12+1.978×107nh13/5EI3z132.11×108nh14/5EI14/5z14+1.69×109nh3EI13/5z151.038×1010nh16/5EI12/5z16+4.91×1012nh17/5EI11/5z171.781×1013nh18/5EI2z18+4.878×1015nh19/5EI9/5z199.76×1017nh4EI8/5z20+1.349×1018nh21/5EI7/5z211.151×1020nh22/5EI6/5z22+4.57×1023nh23/5EIz23).(35)

In this case, the vector b is written as b=[000(2L)3HT3EpIp00000]T. The approximate solution given in Eq. (28) becomes

y(η)=1.599HEI2/5nh3/5(6.15×106+1.45×104η8.656×104η21.82×103η3+3.95×102η40.153η5+0.147η6+0.527η71.94η8+3.62η94.53η104.10η11+30.45η1237.298η1328.146η14+94.987η1530.193η1690.098η17+72.802η18+35.205η1954.997η20+0.553η21+19.385η224.528η232.727+η25).(36)

By using the inverse linear transformation η=(2L)x1andx=zT, the solutions of the governing differential equation for a laterally loaded pile on elastic subgrade, Eq. (12), are given for deflection and moment as

y(z)=0.928Hnh3/5(EI50.50EI23/5z2nh2/5+0.179EI22/5z3nh3/52.72×106EI21/5z4nh4/58.308×103EI4z5nh6.2×105EI19/5z6nh6/5+6.75×104EI18/5z7nh7/51.7×104EI17/5z8nh8/5+3.5×105EI16/5z9nh9/51.21×105EI3z10nh2+4×106EI14/5z11nh11/59.336×107EI13/5z12nh12/5+1.557×107EI12/5z13nh13/51.944×108EI11/5z14nh14/5+1.870×109EI2z15nh31.414×1010EI9/5z16nh16/5+8.484×1012EI8/5z17nh17/54.048×1013EI7/5z18nh18/5+1.529×1014EI6/5z19nh19/54.526×1016EIz20nh4+1.026×1017EI4/5z21nh21/51.725×1019EI3/5z22nh22/5+2.025×1021EI2/5z23nh23/51.482×1023EI1/5z24nh24/5+5.089×1026z25nh5)(37)

and

M(z)=0.927HEInh4/5(EI23/51.078EI22/5znh1/5+3.2×105EI21/5z2nh2/5+0.166EI4z3nh3/5+1.8×103EI19/5z4nh4/50.028EI18/5z5nh+9.624×103EI17/5z6nh6/52.5×103EI16/5z7nh7/5+1.09×103EI3z8nh8/54.4×104EI14/5z9nh9/5+1.2×104EI13/5z10nh22.43×105EI12/5z11nh11/5+3.5×106EI11/5z12nh12/53.93×107EI2z13nh13/5+3.397×108EI9/5z14nh14/52.31×109EI8/5z15nh3+1.24×1010EI7/5z16nh16/55.236×1012EI6/5z17nh17/5+1.72×1013EIz18nh18/54.316×1015EI4/5z19nh19/5+7.979×1017EI3/5z20nh41.026×1018EI2/5z21nh21/5+8.186×1021EI1/5z22nh22/53.056×1023z23nh23/5).(38)

5  Runge-Kutta Formulas of Order 4 and 5 (RKBS45)

The differential equation of laterally loaded pile, Eq. (12), subjected to the boundary conditions in Eqs. (14) or (15) is converted into a system of linear first-order differential equations as follows

y=y1,y1=y2,y2=y3,y3=xy.(39)

If the pile top condition is a free-head one, then the system of equations given in Eq. (39) is subjected to the initial conditions

y(0)=s1,y1(0)=s2,y2(0)=MT2EpIp,y3(0)=HT3EpIp,(40)

where the constants s1 and s2 are unknowns to be evaluated during the process of solution.

On the other hand, if the pile top condition is a fixed-head one, then the system of equations given in Eq. (39) is subjected to the initial conditions

y(0)=s3,y1(0)=0,y2(0)=s4,y3(0)=HT3EpIp,(41)

where the constants s3 and s4 are unknowns to be determined as a part of the numerical solution.

The system of equations Eq. (39) can be rewritten in the vector form as

Y(x)=F(x,Y(x)),0xxmax,(42)

where Y(x)=[y(x)y1(x)y2(x)y3(x)] and xmax is used as an adaptation for x.

To obtain local error estimates for adaptive step-size control effectively, consider two Runge-Kutta formulas of different orders p and p^=p+1. A Runge-Kutta process generates a sequence Yn as an approximation of Y(xn) for 0=x0<x1<<xn=xmax. In the interval from xn to xn+1=xn+hn, there are two approximations of Y(xn+1) called Yn+1 and Y^n+1 for p and p^=p+1, respectively. Their forms are

Yn+1=Yn+hni=0sbiKiandY^n+1=Yn+hni=0sb^iKi,(43)

where s is the number of stages and,

K0=F(xn,Yn)andKi=F(xn+cihn,Yn+hnj=0i1(ai,jKj)),i=1,2,,s,(44)

Considering p=4 and s=7, the coefficients ai,j, bi, b^i and ci=j=0i1ai,j, i=1,2,,s, can be evaluated by Bogacki and Shampine [24] to produce the efficient pair of formulas Yn+1 (4th order formula) and Y^n+1 (5th order formula).

The solution Y(xmax)=[y(xmax)y1(xmax)y2(xmax)y3(xmax)] is generated based on the unknowns s1 and s2 for free-head pile and the unknowns s3 and s4 for fixed-head pile cases. The values of all these unknowns are determined upon solving the equations y(xmax)=y1(xmax)=y2(xmax)=y3(xmax)=0 to satisfy the given boundary conditions given in Eqs. (14) and (15). The value of xmax is adapted by gradual increase till the maximum difference between any values of s1, s2, s3 or s4 in succession is less than the required accuracy τ=106.

The error between the two numerical solutions Yn+1 and Y^n+1 is calculated as

en+1=max|Yn+1Y^n+1|.(45)

In case of en+1τ, one can use Yn+1 or Y^n+1 as the final approximate value of Y(xn+1). On the other hand, if en+1>τ, the error en+1 is used to adapt the step size hn to h^n as follows

h^n=hn(τen+1)1/p.(46)

The adapted step size is used to estimate the new values of Yn+1 and Y^n+1 until achieving en+1τ.

6  Validation of the Proposed Methods

The performance and capability of the proposed methods to predict the behavior of laterally loaded piles in cohesion1ess soil have been demonstrated by comparing the obtained numerical results and the observed results from field experiments in full-scale lateral load tests reported by Cox. et al. [33]. In these tests, the flexible free-head steel tube pipe pile of 0.61 m in diameter, 21 m in length, 9.525mm in wall thickness and 163000 kN⋅m2 bending rigidity, was embedded in a deposit of submerged sand. The soil profile at the site is composed of a uniformly fine-graded sand with internal angle of friction, ϕ=39 and submerged density of ɣ=10.4 kN/m3. To investigate the soil properties below the ground surface, the Standard Penetration Test, SPT, was performed by [33] and reported as 18 blows per 30 cm. Accordingly, the recommended value of the horizontal subgrade reaction constantnh for submerged sand is estimated to be 15 MN/m3, [34]. Moreover, an analytical power series solution is also used for comparison, including the solution proposed by Fayun et al. [11], which reported a simplified analytical solution for laterally loaded long piles based on Fourier-Laplace integral. The computed pile head deflections are plotted and compared with both field experiments results and results from Fayun et al. as shown in Fig. 4. The figure indicates that the proposed numerical solutions are consistent with both original analytic results and experimental test results, which further validates the proposed numerical solutions. The figure clearly indicates that the presented numerical solutions are more suitable for analyzing the response of laterally loaded long piles in sand layers.

Figure 4: Comparison between the numerical, analytical and experimental lateral deflections of pile head vs. applied lateral force

7  Results and Discussion

The design of pile foundation under lateral loads is extensively bound to study both the pile head deflection and the maximum bending moment. Legendre-Galerkin method and Runge-Kutta formulas of order four and five were employed to solve the flexural equation of long piles embedded in homogeneous cohesionless soil with a modulus of subgrade reaction increases linearly with depth. In the numerical simulation, the pile is subjected to horizontal force H0= 100 kN at distance e=0.305m above-soil level, as shown in Fig. 1. The considered parameters in terms of pile length L, pile bending rigidity EPIP and coefficient of subgrade reaction nh, are set to be 21 m, 163000 kN/m2 and 15000 kN/m3 respectively. In order to ensure that our results produce reliable pile deflection and bending moment not only at the pile head but also for the entire pile length, and the analysis is equally applicable for different pile head conditions, the obtained results are compared with the results of simplified solution proposed by Fayun et al. [11]. The responses of free-head pile foundation in terms of lateral deflection and bending moment distributions are shown in Fig. 5. In addition, the corresponding responses in the case of fixed-head pile foundation are also presented in Fig. 6.

Figure 5: Pile deflection and bending moment profiles of free-head long pile

The plotted curves reveal complete overlaps between Legendre-Galerkin solution and Runge-Kutta solution. Furthermore, the comparison charts demonstrate a very good agreement between the numerical results estimated via Legendre-Galerkin and Runge-Kutta and the corresponding procedure introduced by Fayun et al. [11]. These observations suggest that the lateral deflection and bending moment profiles can be represented accurately by the proposed method. However, insignificant differences between the location of peak deflection and bending moment from proposed methods and solution proposed by Fayun et al. can also be observed regardless the type of head condition.

Figure 6: Pile deflection and bending moment profiles of fixed-head long pile

8  Conclusion

Funding Statement: The authors received no specific funding for this study.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

## References

1.  A. Chore, R. Ingle and V. Sawant, “Parametric study of laterally loaded pile groups using simplified FE models,” Coupled Systems Mechanics, vol. 1, no. 1, pp. 1–7, 2012. [Google Scholar]

2.  A. Mostofi and K. Bargi, “Analytical and numerical evaluation of flexible response of floating piers to ship berthing impact,” International Journal of Structural and Civil Engineering Research, vol. 2, no. 1, pp. 249–259, 2011. [Google Scholar]

3.  H. Matlock and C. Reese, “Generalized solutions for laterally loaded piles,” Journal of the Soil Mechanics and Foundations Division, vol. 86, no. 5, pp. 63–91, 1960. [Google Scholar]

4.  M. Heidari, M. Jahanandish, H. El Naggar and A. Ghahramani, “Nonlinear cyclic behavior of laterally loaded pile in cohesive soil,” Canadian Geotechnical Journal, vol. 51, no. 2, pp. 129–143, 2014. [Google Scholar]

5.  D. Su and Y. Zhou, “Effect of loading direction on the response of laterally loaded pile groups in sand,” International Journal of Geomechanics, vol. 16, no. 2, pp. 1–12, 2016. [Google Scholar]

6.  D. Su and J. Li, “Three-dimensional finite element study of a single pile response to multidirectional lateral loadings incorporating the simplified state-dependent dilatancy model,” Computers and Geotechnics, vol. 50, no. 5, pp. 129–142, 2013. [Google Scholar]

7.  S. Gleser, “Lateral load tests on vertical fixed-head and free-head piles,” Symp. on Lateral Load Tests on Piles, ASTM STP, vol. 154, pp. 75–101, 1953. [Google Scholar]

8.  M. Ashour and G. Norris, “Modelling lateral soil-pile response based on soil-pile interaction,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 126, no. 5, pp. 420–428, 2000. [Google Scholar]

9.  Z. Moayed, A. Judi and B. Rabe, “Lateral bearing capacity of piles in cohesive soils based on soil's failure strength control,” Electronic Journal of Geotechnical Engineering, vol. 13(Dpp. 1–11, 2008. [Google Scholar]

10. Y. Chang, “Discussion on lateral pile loading tests,” L.B. Feagin. Trans. ASCE, vol. 102, pp. 272–278, 1937. [Google Scholar]

11. L. Fayun, L. Yanchu, L. Lei and W. Jialai, “Analytical solution for laterally loaded long piles based on Fourier-laplace integral,” Applied Mathematical Modelling, vol. 38, pp. 5198–5216, 2014. [Google Scholar]

12. E. Conte, A. Troncone and M. Vena, “Nonlinear three-dimensional analysis of reinforced concrete piles subjected to horizontal loading,” Computers and Geotechnics, vol. 49, pp. 123–133, 2013. [Google Scholar]

13. S. Gleser, “Lateral load tests on vertical fixed-head and free-head piles,” Symp. on Lateral Load Tests on Piles, ASTM STP, vol. 154, pp. 75–101, 1953. [Google Scholar]

14. P. Ruesta and F. Townsend, “Evaluation of laterally loaded pile group at roosevelt bridge,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 123, no. 12, pp. 1153–1161, 1997. [Google Scholar]

15. D. Brown, L. Resse and M. O'Neill, “Cyclic lateral loading of a large-scale pile group,” Journal of Geotechnical Engineering, vol. 113, no. 11, pp. 1326–1343, 1987. [Google Scholar]

16. S. Yousefi, “Legendre wavelets method for solving differential equations of lane-emden type,” Applied Mathematics and Computation, vol. 181, pp. 1417–1422, 2006. [Google Scholar]

17. J. Shen, “Efficient spectral-galerkin method i. direct solvers for second- and fourth-order equations by using legendre polynomials,” SIAM Journal on Scientific Computing, vol. 15, pp. 1489–1505, 1994. [Google Scholar]

18. X. Ye, “The legendre collocation method for the cahn-hilliard equation,” Applied Mathematical Modelling, vol. 150, pp. 87–108, 2003. [Google Scholar]

19. A. Khater, A. Shamardan, D. Callebaut and M. Sakran, “Numerical solutions of integral and integro-differential equations using legendre polynomials,” Numerical Algorithms, vol. 46, pp. 195–218, 2007. [Google Scholar]

20. B. Bialecki and A. Karageorghis, “Legendre gauss spectral collocation for the helmholtz equation on a rectangle,” Numerical Algorithms, vol. 36, pp. 203–227, 2004. [Google Scholar]

21. W. Zhengsu, C. Yanping and H. Yunqing, “Legendre spectral galerkin method for second-kind volterra integral equations,” Frontiers of Mathematics in China, vol. 4, pp. 181–193, 2009. [Google Scholar]

22. M. Fathy, M. El-Gamel and M. El-Azab, “Legendre-galerkin method for the linear fredholm integrodifferential equations,” Applied Mathematics and Computation, vol. 243, pp. 789–800, 2014. [Google Scholar]

23. S. Yousefi, “Numerical solution of abel's integral equation by using legendre wavelets,” Applied Mathematics and Computation, vol. 175, pp. 574–580, 2006. [Google Scholar]

24. P. Bogacki and L. Shampine, “An efficient runge-kutta (4, 5) pair,” Computers & Mathematics with Applications, vol. 32, pp. 15–28, 1996. [Google Scholar]

25. E. Elbashbeshy, T. Emam, M. El-Azab and K. Abdelgaber, “Effect of thermal radiation on flow, heat and mass transfer of a nanofluid over a stretching horizontal cylinder embedded in a porous medium with suction/injection,” Journal of Porous Media, vol. 18, no. 3, pp. 215–229, 2015. [Google Scholar]

26. E. Elbashbeshy, T. Emam, M. El-Azab and K. Abdelgaber, “Slip effect on flow, heat, and mass transfer of a nanofluid over a stretching horizontal cylinder in the presence of suction/injection,” Thermal Science, vol. 20, no. 6, pp. 1813–1824, 2016. [Google Scholar]

27. A. Hanna and A. Sharif, “Drag force on a single pile in clay subjected to surcharge loading,” International Journal of Geo-Mechanics, ASCE, vol. 6, no. 2, pp. 89–99, 2006. [Google Scholar]

28. L. Xu, F. Cai, G. Wang and K. Ugai, “Nonlinear analysis of laterally loaded single piles in sand using modified strain wedge model,” Computers and Geotechnics, vol. 51, pp. 60–71, 2013. [Google Scholar]

29. H. Poulos and E. Davis, Pile Foundation Analysis and Design, New York: Wiley, 1980. [Google Scholar]

30. L. Reese and H. Matlock, “Non-dimensional solutions for laterally loaded piles with soil modulus assumed proportional to depth,” in  Proc. of the 8th Texas Conf. on Soil Mechanics and Foundation Engineering, Austin, TX, pp. 1–23, 1956. [Google Scholar]

31. N. Lebedev and R. Silverman, Special Functions and Their Applications, New York: Dover Publications, 1972. [Google Scholar]

32. E. Doha, “On the construction of recurrence relations for the expansion and connection coefficients in series of jacobi polynomials,” Journal of Physics A: Mathematical and Theoretical, vol. 37, no. 3, pp. 657, 2004. [Google Scholar]

33. W. Cox, L. Reese and B. R. Grubbs, “Field testing of laterally loaded piles in sand,” in Proc. of the 6th Annual Offshore Technology Conf., Houston, TX, Paper No. OTC 2079, vol. 2, pp. 459–472, 1974. [Google Scholar]

34. R. Liang, “Drilled shaft foundations for noise barrier walls and slope stabilization,” Final Report, FHWA/OH-2002/038, Ohio Department of Transportation, OH, 2002. https://www.semanticscholar.org/paper/Drilled-Shaft-Foundations-for-Noise-Barrier-Walls-Liang/0399dd8a66a13b23de4e984975569f89b17a5e4d [Google Scholar]

 This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.