"

6.26 Rigorous elasticity method for estimating pile deformations

As previously mentioned, Brom’s method for estimating pile head deflection requires estimating the coefficient of subgrade reaction Kh by means of empirical formulas. As such, the method is inherently approximate. A rigorous method for estimating the deflection but also the rotation of the head of a single pile has recently been developed by Zheng et al. (2024). Their formulation is based on modelling soil as elastic Tajimi-type medium, and ignoring possible slippage or separation at the soil-pile interface, therefore is applicable mainly to serviceability conditions. Zheng et al.’s method allows estimating the swaying Khh, rocking Krr and swaying-rocking Krh = Khr components of the head stiffness of a single floating pile embedded in elastic soil of finite thickness (Figure 6.109a).

It is reminded that the swaying stiffness Khh is the lateral load Hw that, when applied to the head of a pile results in head deflection y(0) = 1 (Figure 6.109b). The rocking stiffness Krr is the bending moment Mw that, when applied to the head of a pile results in head rotation ω(0) = 1 (Figure 6.109b). The swaying-rocking stiffness Krh is the bending moment that will develop to a pile with fixed head rotation ω(0) = 0 when the pile is subjected to unit displacement at its head y(0) = 1, considering the sign convention in Figure 6.109b. The swaying-rocking stiffness Khr is the lateral force that will develop to a pile with fixed head displacement y(0) = 0 when the pile is subjected to unit rotation at its head ω(0) = 1, considering the sign convention in Figure 6.109b. Moreover, it is Krh = Khr.

Using Zheng et al.’s method to calculate the stiffness components allows establishing the stiffness [K] and compliance [S] matrixes of the pile, as:

(6.128a) \left[ {\begin{array}{*{20}{c}}{{H_w}}\\{{M_w}}\end{array}} \right] = \left[ {\bf{K}} \right]\left[ {\begin{array}{*{20}{c}}{y(0)}\\{\omega (0)}\end{array}} \right] = \left[ {\begin{array}{*{20}{c}}{{K_{hh}}}&{{K_{hr}}}\\{{K_{rh}}}&{{K_{rr}}}\end{array}} \right]\left[ {\begin{array}{*{20}{c}}{y(0)}\\{\omega (0)}\end{array}} \right]

(6.128b) \left[ {\begin{array}{*{20}{c}}{y(0)}\\{\omega (0)}\end{array}} \right] = \left[ {\bf{S}} \right]\left[ {\begin{array}{*{20}{c}}{{H_w}}\\{{M_w}}\end{array}} \right] = {\left[ {\bf{K}} \right]^{ - 1}}\left[ {\begin{array}{*{20}{c}}{{H_w}}\\{{M_w}}\end{array}} \right]

It is reminded that the stiffness and compliance matrixes are independent of the load conditions, and of any restrains imposed to pile head deflection or rotation. The charts provided in Figure 6.110 allow calculating directly the components of the stiffness matrix for soil with Poisson’s ratio vs = 0.3 and different values of the pile-soil moduli ratio Ep/Es, the slenderness of the pile L/D and the dimensionless thickness of the foundation soil layer H/L. Having estimated these components one can calculate the compliance matrix by inverting the stiffness matrix, and subsequently the deflection y(0) and the rotation ω(0) of the pile head for any load combination acting at the pile head, as well as head fixity conditions.

It is interesting to observe in Figure 6.110 that there exists a critical pile slenderness ratio L/D, beyond which pile stiffness becomes independent of the length of the pile, and the bottom part of the pile no longer contributes to its response to lateral load. Therefore one can argue that, if the top part of the soil profile of dimensionless thickness equal to the critical pile slenderness comprises uniform soil conditions, this method (but also Broms method) will provide reasonably accurate results, even when the pile is installed in layered soil. This is not limited to serviceability conditions, but can be extrapolated to the estimation of the collapse lateral load. However, the critical pile slenderness is not the same for all components of stiffness, and on top of that it varies with the pile-soil moduli ratio Ep/Es.

Figure (a) on the left presents a schematic of a friction pile embedded in a soil layer of thickness H, that is underlaid by rigid base. The pile is subjected to a horizontal force Hw and a bending moment M=ehHw on its head. The pile's length is L<H and the pile's diameter is D. The vertical axis is denoted with z. The parameters of the soil are Young's modulus Es, Shear modulus Gs, Poisson's ratio vs. The parameters of the pile are Young's modulus Ep, moment of inertia Ip, area of cross section Ap. The soil reaction acting on the pile is denoted with p(z). Figure (b) on the right presents two free body diagrams of the pile, which are used to defined the swaying stiffness Khh, the swaying-rocking stiffness Krh=Khr, and the rocking stiffness Krr. In the first diagram the pile is subjected to unit deflection at its head and in the second diagram the pile is subjected to unit rotation at its head.
Figure 6.109. (a) Schematic of a floating pile embedded in a soil layer of finite thickness and subjected to lateral load and/or bending moment at its head; (b) Definition of the components of pile stiffness (Zheng et al. 2024).
Twelve charts used for the calculation of the dimensionless components of pile stiffness, as function of the problem parameters. The top row presents the variation of Khh/(EsD), Khr/(EsD^2) and Krr/(EsD^3) with L/D, for different values of the Ep/Es ratio. The plotted results correspond to H/D = 40 and vs = 0.3.The second row presents the variation of Khh/(EsD), Khr/(EsD^2) and Krr/(EsD^3) with Ep/Es, for different values of slenderness L/D. The plotted results correspond to H/D = 40 and vs = 0.3.The third row presents the variation of Khh/(EsD), Khr/(EsD^2) and Krr/(EsD^3) with Ep/Es, for different values of H/L. The plotted results correspond to L/D = 5 and vs = 0.3. The bottom row presents the variation of Khh/(EsD), Khr/(EsD^2) and Krr/(EsD^3) with H/L, for different values of the Ep/Es ratio. The plotted results correspond to L/D = 10 and vs = 0.3. All the parameters are defined in the previous Figure 6.109.
Figure 6.110. Variation of pile head stiffness components with pile-soil moduli ratio Ep/Es, pile slenderness L/D and dimensionless thickness of the foundation soil layer H/L. Results for soil Poisson’s ratio vs = 0.3 (Results from Zheng et al. 2024).

License

Icon for the Creative Commons Attribution 4.0 International License

Fundamentals of foundation engineering and their applications Copyright © 2025 by University of Newcastle & G. Kouretzis is licensed under a Creative Commons Attribution 4.0 International License, except where otherwise noted.