Next Article in Journal
The Stability Analysis of a Double-X Queuing Network Occurring in the Banking Sector
Previous Article in Journal
Geometric Numerical Integration of Liénard Systems via a Contact Hamiltonian Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solution of Laterally Loaded Free-Head Long Piles in Elasto-Plastic Cohesive Soils

by
Ayman Abd-Elhamed
1,2
1
Physics and Engineering Mathematics Department, Faculty of Engineering-Mattaria, Helwan University, Cairo 11718, Egypt
2
Faculty of Engineering, King Salman International University, Ras-Sedr 46512, Egypt
Mathematics 2021, 9(16), 1961; https://doi.org/10.3390/math9161961
Submission received: 18 July 2021 / Revised: 11 August 2021 / Accepted: 13 August 2021 / Published: 16 August 2021
(This article belongs to the Section Engineering Mathematics)

Abstract

:
This research study presents a closed form solution of responses of laterally loaded long piles embedded on cohesive soils with a constant subgrade modulus. The surrounding soil medium is modelled as elastic-perfectly plastic. The closed form solution is derived by solving the governing differential equation of the pile–soil system. The most popular numerical computation software package MATLAB is utilized for the implementation of solutions. The provided analytical method reliably calculates the pile head deflection and bending moment required for engineering design purposes. Results are discussed and verified with solutions of an equivalent three-dimensional finite element (FE) model developed using ANSYS software. It was concluded that the proposed analytical model could efficiently provide the exact solution of embedded piles in elasto-plastic cohesive soil under lateral loads.

1. Introduction

When a shallow soil stratum is highly compressible and of low bearing capacity to support the load transmitted by the superstructure, piles are principally used to transmit the loads to the underlying stiffer soil strata of a higher bearing capacity. Various precious approaches have been proposed to analyze the laterally loaded piles, including the subgrade reaction approach [1,2] and p-y approach [3,4]. In the subgrade reaction approach, the supporting soil is idealized as a series of lateral elastic springs and the relationship between lateral applied force and lateral deflection pattern of a pile is treated as linear- elastic over the depth of the pile. Moreover, the subgrade reaction modulus is obtained through fitting with relevant numerical solutions [5,6]. In the p-y approach the soil is conveniently simulated by a series of independent non-linear springs varying with the depth, and the soil–pile interface responses are captured by representing the soil nonlinear springs with p-y curves along the pile length. Several researchers on the piles subjected to lateral loads continued to be unabated in which the soil is considered as a continuum-based model. In these studies, different approaches such as the finite element approach [7], boundary element approach [8], and finite difference approach [9], have been introduced to predict the response of laterally loaded piles. A large number of studies have examined the behaviour of the soil–pile interaction of laterally loaded piles considering that the relationship between soil lateral reaction and pile deflection is linear [10], and the lateral load behaviour of soil is non-linear. Therefore, the development of algebraic expressions via a closed form solution is proposed in this study to predict the non-linear responses of laterally loaded piles in clay strata. The proposed analytical solution in the current study provides a better approach for structural designers to simply solve for the lateral soil reaction, displacement, shear force, and bending moment responses of laterally loaded free-head long piles embedded in homogeneous cohesive soils with a constant lateral subgrade reaction with depth. Consequently, the techniques can be easily applied in practice as an alternative approach to analyse and design laterally loaded long piles. In addition, the proposed solution incorporates the contribution of the Winkler foundation model with the linear elastic-perfect-plastic p-y soil response. The proposed analytical solution is verified by comparing the results obtained for the pile’s deflection and induced bending moment with numerical results using the finite element software ANSYS. The accuracy of the location at which the maximum bending moment occurs, and the magnitude of that maximum bending moment along the pile are further verified analytically by comparing them to those obtained from the FE analysis. In addition, a numerical example is introduced to investigate the effect of the horizontal loads on the depth of the plastic zone.

2. Problem Definition and Formulation

A free-head cylindrical pile of length L, diameter D, and Young’s modulus E p embedded in a homogeneous clay soil and subjected to lateral load H 0 at an eccentricity, e above ground surface is schematically shown in Figure 1a. In order to simulate the pile–soil interaction, the pile shaft is assumed to be perfectly glued to the surrounding soils suggesting that there is no relative movement along the pile–soil interface. Furthermore, a series of springs distributed along the pile shaft are used to characterize that interaction for the pile subjected to the lateral load H 0 and bending moment, M 0 = H 0 e at the ground surface, as shown in Figure 1b.
Free-head and fixed-tip boundary conditions are assumed. Therefore, the pile head is unrestrained to rotate and move laterally while the pile base is fully fixed against translation and rotation. Unquestionably, the soil behaves non-linearly under higher lateral load. This nonlinearity is described by the well-known ideal elasto-plastic p-y relationship at any depth. This model is shown in Figure 2a in which k is the gradient of the p-y curve that offers the soil modulus at any point z below the surface along the pile, p u is the limiting soil horizontal resistance per unit depth, and y_y is the threshold of soil displacement above which the soil is initiated to yield. On the other hand, the soil near the ground surface presumably yields firstly and propagates downward to a depth called plastic depth, z p [11]. Below plastic depth, the soil reaction is linear and still in the elastic state. Accordingly, the lateral soil reaction along the pile shaft may be characterized by two zones: the plastic zone above which extends to a depth z p depending on both the applied lateral force and its point of application, and the bottom elastic zone below this depth (see Figure 1b and Figure 2b) [12]. It is worth noting that the limiting lateral force for clay soils increases with depth from 2 c u D to 3 c u D at ground surface up to 8 c u D to 12 c u D at a depth of about 3D, and beyond this depth it keeps constant [13]. Where c u is the undrained shear strength of the soil, D is the pile diameter. Consequently, the limiting lateral force is proposed to equal 9 c u D to the end of the plastic zone (see Figure 2b) [14].

3. Mathematical Formulation

In the case of the pile subjected to lateral forces, it tries to push the surrounding soil horizontally in the loading direction to its limit state. Subsequently, the surface soils tend to yield, and the cracks tend to be very shallow. The lateral soil reaction along the pile shaft may be characterized by two zones: the top plastic zone which extends downward to depth z p depending on both applied lateral force and its point of application, and the bottom elastic zone below this depth (see Figure 2).
The governing flexural equation for the response of the soil–pile system according to Winkler’s foundation concept can be written as:
E I d 4 y d z 4 + P ( z ) = 0 ,
In the above, E is the modulus of elasticity of concrete, I is the effective moment of inertia, y is the lateral deflection of the pile, and p(z) is the soil reaction at any point at a depth z along the axis of the pile which may be described for the elasto-plastic model as:
P ( z ) = {   k × y         z z p     p u                 z   z p ,
where k is the subgrade reaction modulus of the soil, z p is the depth of the plastic zone, and p u is the lateral limit force below which the soil reaction has a linear relationship with lateral pile deflection.
The integration of Equation (1) considering the boundary and continuity conditions presented in Table 1 and Table 2 leads to the analytical solution for piles subjected to lateral loads and embedded in linear-plastic soil.

3.1. Elastic Soil Behaviour

Based on the assumption that the applied lateral force does not exceed the ultimate load p u , the relationship between lateral soil resistance per unit pile length and pile lateral deflection under working loads is linear–elastic. Consequently, under this assumption the governing equation of motion can be expressed as:
E I d 4 y d z 4 + k × y = 0 ,
If one assumes that the modulus of subgrade reaction k in clay soil is independent of depth [15], thereby the general solution of Equation (3) is of the form:
y ( z ) = e β z { g 1 sin ( β z ) + g 2 cos ( β z ) } + e β z { g 3 sin ( β z ) + g 4 cos ( β z ) } ,
where β is the reciprocal of the characteristic length and g i ,   i = 1 , 2 , 3 , 4 are four independent constant parameters.
Assuming that the pile is infinite and the deflection at the tip of the long piles is negligible, then the lateral deflections at working loads can be simplified as
y ( z ) = e β z { g 1 sin ( β z ) + g 2 cos ( β z ) } ,
Further higher derivatives for Equation (5) yields
y ( z ) = β e β z { ( g 1 g 2 ) cos ( β z ) + ( g 1 + g 2 ) sin ( β z ) } ,
y ( z ) = 2 β 2 e β z { g 1 cos ( β z ) g 2 sin ( β z ) } ,
y ( z ) = 2 β 3 e β z { ( g 1 + g 2 ) cos ( β z ) + ( g 1 g 2 ) sin ( β z ) } ,
y ″″ ( z ) = 4 β 4 e β z { g 1 sin ( β z ) + g 2 cos ( β z ) } ,
By substituting Equations (5) and (9) into Equation (3), β for the free-load pile can be assigned as:
β = 1 / L = k D / 4 E I 4 ,
where EI is stiffness of pile section and D is width or diameter of pile.
Two constant parameters g 1   and   g 2 are obtained for piles subjected to applied horizontal load H 0 and applied bending moment M 0 by imposing the equilibrium conditions for both the bending moment and shear force at pile top using Table 1. Thereupon, by substituting the boundary conditions into Equations (7) and (8), two integration coefficients g 1   and   g 2 can be determined as:
g 1 = 1 2 β 2 E I { M 0 } ,
g 2 = 1 2 β 2 E I { M 0 + H 0 β } ,
The deflection along the laterally loaded pile can be obtained by substituting Equations (11) and (12) into Equation (5) as
y ( z ) = e β z 2 β 2 E I   { M 0 sin ( β z ) + ( M 0 + H 0 β ) cos ( β z ) } ,
Consequently, the bending moment M, and the shear force V are found by successive differentiation of Equation (13) as:
M ( z ) = e β z { M 0 cos ( β z ) + ( M 0 + H 0 β ) sin ( β z ) } ,
V ( z ) = e β z { H 0 cos ( β z ) ( 2 β M 0 + H 0 ) sin ( β z ) } ,
By setting V ( z ) = 0 , the point of zero shear z 0 can be written as:
z 0 = 1 β tan 1 H 0 2 β M 0 + H 0 ,
Therefore, the maximum bending moment can be obtained as follows:
M m a x = e β z 0 { M 0 cos ( β z 0 ) + ( M 0 + H 0 β ) sin ( β z 0 ) } ,

3.2. Elasto-Plastic Soil Behaviour

Admittedly, the soil behavior under higher lateral load levels is predominantly non-linear. A simple approach proposed by Madhav et al. [16] reasonably captures this nonlinearity through considering the elastic-plastic subgrade model. Consequently, the non-linear responses for long piles under lateral loading may be obtained through two zones as follows.

3.2.1. Plastic Zone z z p

p ( z ) = E I d 4 Y d z 4 = p u ,
Through several integrations for Equation (18), the key responses of the pile in the plastic zone can be expressed as follows:
V ( z ) = E I d 3 Y d z 3 = p u z + C 1 ,
M ( z ) = E I d 2 Y d z 2 = p u z 2 2 + C 1 z + C 2 ,
φ ( z ) = d Y d z = 1 E I [ p u z 3 6 + C 1 z 2 2 + C 2 z + C 3 ] ,
Y ( z ) = 1 E I [ p u z 4 24 + C 1 z 3 6 + C 2 z 2 2 + C 3 z + C 4 ] ,

3.2.2. Elastic Zone z z p

y ( z ) = e β z { C 5 sin ( β z ) + C 6 cos ( β z ) } ,
y ( z ) = β e β z   { ( C 5 C 6 ) cos ( β z ) ( C 5 + C 6 ) sin ( β z ) } ,
y ( z ) = 2 β 2 e β z { C 5 cos ( β z ) C 6 sin ( β z ) } ,
y ( z ) = 2 β 3 e β z { ( C 5 + C 6 ) cos ( β z ) + ( C 5 C 6 ) sin ( β z ) } ,
y ″″ ( z ) = 4 β 4 e β z { C 5 sin ( β z ) + C 6 cos ( β z ) } ,
The solution of the above equations requires boundary conditions to obtain the unknown parameters C i ,   i = 1 , 2 , 3 , 4 , 5 , 6 and z = z p . All boundary conditions associated with the above equations must be determined by recourse to the boundary conditions in Table 2. The equilibrium for the bending moment and shear force at z = 0 and compatibility for displacement and slope at z = z p is considered.
By considering equilibrium for shear at the pile head, the undetermined constants C 1 can be determined as:
C 1 = H 0 ,
Subsequently, by imposing the equilibrium for the bending moment at the pile head, the undetermined constants C 2 can be expressed as:
C 2 = M 0 ,
Then, the plastic zone develops at a depth z p of the form:
z p = ( H 0 2 p u 2 + 2 M 0 p u ) 0.5 + H 0 p u 1 β ,
Applying the compatibility for slope at z = z p by equating Equation (21) with Equation (24) yields
C 3 = E I β { y ( z p ) + M ( z p ) 2 β 2 E I } + p u z p 3 6 H 0 z p 2 2 M 0 z p ,
where
y ( z p ) = p u k ,
and
M ( z p ) = p u z p 2 2 + H 0 z p + M 0 ,
Similarly, applying the compatibility for deflection at z = z p by equating Equation (22) with Equation (23) yields,
C 4 = E I y ( z p ) + p u z p 4 24 H 0 z p 3 6 M 0 z p 2 2 C 3 z p ,
Based on Equations (21) and (22), the rotation φ t o p and horizontal deflection Y t o p of the pile at z = 0 can be expressed as
φ t o p = β { y ( z p ) + M ( z p ) 2 β 2 E I } + p u z p 3 6 H 0 z p 2 2 M 0 z p ,
Y t o p = y ( z p ) + 1 E I { p u z p 4 24 H 0 z p 3 6 M 0 z p 2 2 C 3 z p } ,
Furthermore, Equation (26) can be re-expressed as:
y ( z ) = 2 β 3 e β z { C 5 cos ( β z ) C 6 sin ( β z ) + C 5 sin ( β z ) + C 6 cos ( β z ) } , y ( z ) = ( β ) ( 2 β 2 e β z ) { C 5 cos ( β z ) C 6 sin ( β z ) } ( β ) y ( z )                                       ( 1 2 β ) ( 4 β 4 e β z ) { C 5 sin ( β z ) + C 6 cos ( β z ) } ( 1 2 β )   y ″″ ( z )
y ( z ) = β y ( z ) + 2 β 3 y ( z ) ,
Applying equilibrium principles at z = z p yields
C 5 = e β z p { y ( z p ) sin ( β z p ) M ( z p ) 2 β 2 E I cos ( β z p ) } ,
and C 6 is given by
C 6 = e β z p { y ( z p ) cos ( β z p ) + M ( z p ) 2 β 2 E I sin ( β z p ) } ,
There are two cases that may arise to obtain the maximum value of the bending moment and the corresponding position.
Case 1: the bending moment occurs in the elastic region.
In this case, the largest bending moment value is assumed to occur in the elastic region. Substituting into Equation (26) with z = z 0 and equating with zero yields
y ( z 0 ) = 2 β 3 E I e β z 0 { e β z p y ( z p ) [ sin ( β z p ) cos ( β z 0 ) + cos ( β z p ) cos ( β z 0 ) + e β z p M ( z p ) 2 β 2 E I [ sin ( β z p ) cos ( β z 0 ) cos ( β z p ) cos ( β z 0 ) ] + e β z p y ( z p ) [ sin ( β z p ) sin ( β z 0 ) cos ( β z p ) sin ( β z 0 ) ] e β z p M ( z p ) 2 β 2 E I [ sin ( β z p ) sin ( β z 0 ) + cos ( β z p ) sin ( β z 0 ) ] } = 0 ,
The depth of maximum bending moment which occurs in elastic zone can be expressed as
z 0 = 1 β tan 1 { sin ( β z p ) [ y ( z p ) + M ( z p ) 2 β 2 E I ] + cos ( β z p ) [ y ( z p ) M ( z p ) 2 β 2 E I ] cos ( β z p ) [ y ( z p ) + M ( z p ) 2 β 2 E I ] sin ( β z p ) [ y ( z p ) M ( z p ) 2 β 2 E I ] } ,
Furthermore,
M m a x = 2 β 2 E I e β z 0 { C 5 cos ( β z 0 ) C 6 sin ( β z 0 ) } ,
Case 2: the bending moment occurs in the plastic region.
This case assumes the occurrence of largest bending moment value is in the plastic region. Proceeding similarly as Case 1, the depth of maximum bending moment which occurs in elastic zone can be expressed as
z 0 = H 0 p u ,
and its value is
M m a x = H 0 2 2 p u + M 0 ,

4. Validation of the Proposed Methods

The capability of the proposed analytical solutions presented in this research to predict the lateral behaviour of single piles embedded in a cohesive soil and the corresponding pile responses have been demonstrated and compared with the FE results carried out using ANSYS software.
For this purpose, a concrete pile with length L = 15 m, diameter D = 0.4 m, Poisson’s ratio ν = 0.2, and Young’s modulus of elasticity Ep = 35,000 MPa is embedded in a clay layer with lateral reaction modulus ks = 50,000 kpa and shear strength Cu = 14.4 kpa. The pile head is unrestrained to rotate and move laterally while the pile tip is completely fixed.
Figure 3 presents the computed pile head deflections versus applied lateral forces for both the analytical solution and FE solution with ANSYS software. The figure indicates that the results of the proposed analytical solution are closer to the results from the FE method. The lateral deflection obtained by the ANSYS workbench 18.2 for different lateral loads, 2, 75.5, and 82 kN are shown in Figure 4. The captured results obtained from Figure 3 are compared with those obtained by proposed analytical method as shown in Table 3. The presented results clearly indicate excellent agreement with the FE solution.
Additionally, to ensure that the proposed solution produces a reliable pile deflection and bending moment not only at the pile head but also for the entire pile length, the lateral deflection and bending moment distribution along the pile shaft are presented and compared with the results of the FE analysis as shown in Figure 5.
The captured peak lateral deflection and bending moment obtained from the analytical solution are 16.710 mm and 140.459 kN·m, respectively, while the corresponding values obtained from the FE analysis are 17.809 mm and 141.053 kN·m, respectively. The results clearly indicate the agreement between the analytical and FE solutions in calculating peak value for design purposes.
Additionally, it is observed that the lateral deflection and bending moment profiles obtained using the proposed analytical method reach a good agreement with the FE solution.

5. Conclusions

An analytical solution of the lateral response of an infinite pile embedded in homogeneous elasto-plastic soil deposits has been proposed in the current study. The solution is based on conceptual assumptions of the Winkler model which treat the pile as a flexible beam and soil restraint surrounding the pile as a series of non-linear independent springs. The flexural differential equation describing the response of the soil–pile system was solved analytically. Additionally, the location at which the maximum bending moment occurs, and the magnitude of that maximum bending moment along the pile were obtained from the analysis. For verification of the proposed analytical procedure, the obtained results are compared with the FE results, carried out using ANSYS software, showing excellent agreements. The proposed exact solution could efficiently provide a better approach for structural designers to predict the non-linear responses of laterally loaded free-head long piles embedded in clay soil deposits with a uniform subgrade reaction modulus. Additionally, the proposed analytical solution has the advantage that it can provide a better approach for engineers to analyse and design the piles subjected to lateral loads without resorting to a 3D FE analysis, which is time consuming. Consequently, the mathematical forms of the theoretical solutions can be easily applied in practice as an alternative approach to analyse and design laterally loaded long piles.
It is emphasized that the simplified analytical solution obtained in this study is limited by the linearity of the soil–pile system, which often significantly affects the overall responses of the pile and the soil. Therefore, the proposed solution must be modified before applied to a non-linear case. Moreover, further investigation and analyses are needed to analyse the laterally loaded long piles in soil with variable subgrade reactions with depth, or for piles embedded in multilayered soil.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Nomenclature

DPile diameter or width H 0 Applied lateral load
LPile length M 0 Applied moment
IMoment of inertia of pile shaft y ( z ) Lateral deflection of pile
EYoung’s modulus of pile shaft φ ( z ) Pile rotation
kCoefficient of the subgrade reaction M ( z ) Pile bending moment
cuUndrained shear strength V ( z ) Pile shear force
puLimit lateral force p ( z ) Load per unit length
β Reciprocal of the characteristic length M m a Maximum moment
eDistance of pile above ground surface Y t o p Top horizontal deflection
g i Integration coefficients φ t o p Top pile rotation
C i Unknown parameters z 0 Point of zero shear
v Soil Poisson’s ratio z p Depth of the plastic zone

References

  1. Hetenyi, M. Beams on Elastic Foundations; University of Michigan Press: Ann Arbor, MI, USA, 1964. [Google Scholar]
  2. Matlock, H.; Reese, C. Generalized solutions for laterally loaded piles. J. Soil Mech. Found. Div. ASCE 1960, 86, 63–91. [Google Scholar] [CrossRef]
  3. Guo, W.D. Simple model for nonlinear response of 52 laterally loaded piles. J. Geotech. Geoenviron. Eng. 2013, 139, 234–252. [Google Scholar] [CrossRef]
  4. Su, D.; Zhou, Y.G. Effect of loading direction on the response of laterally loaded pile groups in sand. Int. J. Geomech. ASCE 2016, 16, 04015051. [Google Scholar] [CrossRef]
  5. Baguelin, F.; Frank, R.; Said, Y.H. Theoretical study of lateral reaction mechanism of piles. Geotechnique 1977, 27, 405–434. [Google Scholar] [CrossRef]
  6. Scott, R.F. Foundation Analysis; Prentice Hall: Englewood Cliffs, NJ, USA, 1981. [Google Scholar]
  7. Basu, D.; Salgado, R. Elastic analysis of laterally loaded pile in multi-layered soil. Geom. Geoeng. 2007, 2, 183–196. [Google Scholar] [CrossRef]
  8. Ai, Z.Y.; Feng, D.L.; Cheng, Y.C. BEM analysis of laterally loaded piles in multi-layered transversely isotropic soils. Eng. Anal. Bound. Elem. 2013, 37, 1095–1106. [Google Scholar] [CrossRef]
  9. Zhang, L.; Ernst, C.D.; Einstein, H.H. Nonlinear analysis of laterally loaded rock-socketed shafts. J. Geotech. Geoenviron. Eng. ASCE 2000, 126, 955–968. [Google Scholar] [CrossRef]
  10. William, H.; Celio, V.; Dipanjan, B.; Griffiths, D.V. Elastic solutions for laterally loaded piles. J. Geotech. Geoenviron. Eng. 2013, 139, 1096–1103. [Google Scholar]
  11. Hsiung, Y.M. Theoretical elastic-plastic solution for laterally loaded piles. J. Geotech. Geoenviron. Eng. 2003, 129, 475–480. [Google Scholar] [CrossRef]
  12. Fayun, L.; Yanchu, L.I.; Lei, L.L.; Jialai, W. Analytical solution for laterally loaded long piles based on Fourier-Laplace integral. Appl. Math. Model. 2014, 38, 5198–5216. [Google Scholar]
  13. Reese, L.L.; Welch, R.C. Lateral loading of deep foundations in stiff clay. J. Geotech. Eng. ASCE 1975, 101, 633–649. [Google Scholar]
  14. Broms, B.B. Lateral resistance of piles in cohesive soils. J. Soil Mech. Div. ASCE 1964, 90, 5198–5216. [Google Scholar]
  15. Banerjee, P.K.; Davies, T.G. The behaviour of axially and laterally loaded single piles embedded in nonhomogeneous soils. Geotechnique 1978, 28, 309–326. [Google Scholar] [CrossRef]
  16. Madhav, M.R.; Rao, N.S.; Madhavan, K. Laterally loaded pile in elasto-plastic soil. Soils Found. 1971, 11, 1–15. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Schematic diagram of pile under lateral load, (b) illustration of lateral soil reaction distribution along the pile shaft.
Figure 1. (a) Schematic diagram of pile under lateral load, (b) illustration of lateral soil reaction distribution along the pile shaft.
Mathematics 09 01961 g001
Figure 2. (a) Elasto-plastic behavior of the soil, (b) lateral soil reaction distribution [14].
Figure 2. (a) Elasto-plastic behavior of the soil, (b) lateral soil reaction distribution [14].
Mathematics 09 01961 g002
Figure 3. Comparison between the analytical and numerical lateral deflections of pile head versus applied lateral forces.
Figure 3. Comparison between the analytical and numerical lateral deflections of pile head versus applied lateral forces.
Mathematics 09 01961 g003
Figure 4. The lateral deflection profiles of piles for different lateral loads (a) 2 kN, (b) 75.5 kN, and (c) 80 kN applied at a distance of 1 m above ground surface.
Figure 4. The lateral deflection profiles of piles for different lateral loads (a) 2 kN, (b) 75.5 kN, and (c) 80 kN applied at a distance of 1 m above ground surface.
Mathematics 09 01961 g004aMathematics 09 01961 g004b
Figure 5. Comparison of the results with the FE results. (a) The distribution of lateral deflection along the pile. (b) The distribution of the bending moment along the pile.
Figure 5. Comparison of the results with the FE results. (a) The distribution of lateral deflection along the pile. (b) The distribution of the bending moment along the pile.
Mathematics 09 01961 g005
Table 1. Boundary conditions used for elastic solution.
Table 1. Boundary conditions used for elastic solution.
Boundary ConditionEquations
Shear force V = H 0 at z = 0 E I Y ( 0 ) = H 0
Bending moment M = M 0 at z = 0 E I Y ( 0 ) = M 0
Table 2. Boundary conditions used for elasto-plastic solution.
Table 2. Boundary conditions used for elasto-plastic solution.
Boundary ConditionEquations
Shear force V = H 0 at z = 0 E I Y ( 0 ) = H 0
Bending moment M = M 0 at z = 0 E I Y ( 0 ) = M 0
Slope from plastic zone equal slope from elastic zone at z = z p (compatibility for slope) Y ( z p ) p l a s t i c = y ( z p ) e l a s t i c
Deflection from plastic zone equal deflection from elastic zone at z = z p (compatibility for deflection) Y ( z p ) p l a s t i c = y ( z p ) e l a s t i c
Bending moment from plastic zone equal bending moment from elastic zone at z = z p y ( z p ) = M ( z p ) E I
Shear force from plastic zone equal shear force from elastic zone at z = z p y ( z p ) = V ( z p ) E I
For infinitely long piles, the pile is fixed at the bottom y ( L ) = 0
y ( L ) = 0
y ( L ) = 0
y ( L ) = 0
Table 3. Pile head deflections obtained from the analytical and FE methods.
Table 3. Pile head deflections obtained from the analytical and FE methods.
Applied Lateral Load (kN)Pile Head Deflections (mm)%
Analytical MethodFE, ANSYS Method
20.540.4319
75.514.3816.3413.7
8218.317.753
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abd-Elhamed, A. Analytical Solution of Laterally Loaded Free-Head Long Piles in Elasto-Plastic Cohesive Soils. Mathematics 2021, 9, 1961. https://doi.org/10.3390/math9161961

AMA Style

Abd-Elhamed A. Analytical Solution of Laterally Loaded Free-Head Long Piles in Elasto-Plastic Cohesive Soils. Mathematics. 2021; 9(16):1961. https://doi.org/10.3390/math9161961

Chicago/Turabian Style

Abd-Elhamed, Ayman. 2021. "Analytical Solution of Laterally Loaded Free-Head Long Piles in Elasto-Plastic Cohesive Soils" Mathematics 9, no. 16: 1961. https://doi.org/10.3390/math9161961

APA Style

Abd-Elhamed, A. (2021). Analytical Solution of Laterally Loaded Free-Head Long Piles in Elasto-Plastic Cohesive Soils. Mathematics, 9(16), 1961. https://doi.org/10.3390/math9161961

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop