{"id":625,"date":"2025-03-14T04:54:37","date_gmt":"2025-03-14T04:54:37","guid":{"rendered":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/chapter\/6-26-rigorous-elasticity-method-for-estimating-pile-deformations\/"},"modified":"2026-03-16T14:15:15","modified_gmt":"2026-03-16T14:15:15","slug":"6-26-rigorous-elasticity-method-for-estimating-pile-deformations","status":"publish","type":"chapter","link":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/chapter\/6-26-rigorous-elasticity-method-for-estimating-pile-deformations\/","title":{"raw":"6.26 Rigorous elasticity method for estimating pile deformations","rendered":"6.26 Rigorous elasticity method for estimating pile deformations"},"content":{"raw":"As previously mentioned, Brom\u2019s method for estimating pile head deflection requires estimating the coefficient of subgrade reaction <em>K<sub>h<\/sub><\/em> 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 <em>et al.<\/em> (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 <em>et al.<\/em>\u2019s method allows estimating the swaying <em>K<sub>hh<\/sub><\/em>, rocking <em>K<sub>rr<\/sub><\/em> and swaying-rocking <em>K<sub>rh <\/sub><\/em>= <em>K<sub>hr<\/sub><\/em> components of the head stiffness of a single floating pile embedded in elastic soil of finite thickness (<a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109a<\/a>).\n\nIt is reminded that the swaying stiffness <em>K<sub>hh<\/sub><\/em> is the lateral load <em>H<sub>w<\/sub><\/em> that, when applied to the head of a pile results in head deflection <em>y<\/em>(0) = 1 (<a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>). The rocking stiffness <em>K<sub>rr<\/sub><\/em> is the bending moment <em>M<sub>w<\/sub><\/em> that, when applied to the head of a pile results in head rotation <em>\u03c9<\/em>(0) = 1 (<a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>). The swaying-rocking stiffness <em>K<sub>rh<\/sub><\/em> is the bending moment that will develop to a pile with fixed head rotation <em>\u03c9<\/em>(0) = 0 when the pile is subjected to unit displacement at its head <em>y<\/em>(0) = 1, considering the sign convention in <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>. The swaying-rocking stiffness <em>K<sub>hr<\/sub><\/em> is the lateral force that will develop to a pile with fixed head displacement <em>y<\/em>(0) = 0 when the pile is subjected to unit rotation at its head <em>\u03c9<\/em>(0) = 1, considering the sign convention in <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>. Moreover, it is <em>K<sub>rh<\/sub><\/em> = <em>K<sub>hr<\/sub><\/em>.\n\nUsing Zheng <em>et al.<\/em>\u2019s method to calculate the stiffness components allows establishing the stiffness [<strong>K<\/strong>] and compliance [<strong>S<\/strong>] matrixes of the pile, as:\n\n<strong>(6.128a)<\/strong> [latex]\\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}}}&amp;{{K_{hr}}}\\\\{{K_{rh}}}&amp;{{K_{rr}}}\\end{array}} \\right]\\left[ {\\begin{array}{*{20}{c}}{y(0)}\\\\{\\omega (0)}\\end{array}} \\right][\/latex]\n\n<strong>(6.128b) <\/strong>[latex]\\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][\/latex]\n\nIt 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 <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.110-hr.png\">Figure 6.110<\/a> allow calculating directly the components of the stiffness matrix for soil with Poisson\u2019s ratio <em>v<sub>s<\/sub><\/em> = 0.3 and different values of the pile-soil moduli ratio <em>E<sub>p<\/sub><\/em>\/<em>E<sub>s<\/sub><\/em>, the slenderness of the pile <em>L<\/em>\/<em>D<\/em> and the dimensionless thickness of the foundation soil layer <em>H<\/em>\/<em>L<\/em>. Having estimated these components one can calculate the compliance matrix by inverting the stiffness matrix, and subsequently the deflection <em>y<\/em>(0) and the rotation <em>\u03c9<\/em>(0) of the pile head for any load combination acting at the pile head, as well as head fixity conditions.\n\nIt is interesting to observe in <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.110-hr.png\">Figure 6.110<\/a> that there exists a critical pile slenderness ratio <em>L<\/em>\/<em>D<\/em>, 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\u00a0extrapolated 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 <em>E<sub>p<\/sub><\/em>\/<em>E<sub>s<\/sub><\/em>.\n\n[caption id=\"attachment_623\" align=\"aligncenter\" width=\"800\"]<img class=\"wp-image-623 size-full\" src=\"https:\/\/libraryresources.nse.org.ng\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341.png\" alt=\"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.\" width=\"800\" height=\"482\"> 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 <em>et al.<\/em> 2024).[\/caption]\n\n[caption id=\"attachment_623\" align=\"aligncenter\" width=\"936\"]<img class=\"wp-image-2830 size-large\" src=\"https:\/\/libraryresources.nse.org.ng\/wp-content\/uploads\/sites\/9\/2026\/03\/6.110-hr-scaled.png\" alt=\"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.\" width=\"936\" height=\"1024\"> Figure 6.110. Variation of pile head stiffness components with pile-soil moduli ratio <em>E<sub>p<\/sub>\/E<sub>s<\/sub><\/em>, pile slenderness <em>L\/D<\/em> and dimensionless thickness of the foundation soil layer <em>H\/L<\/em>. Results for soil Poisson\u2019s ratio <em>v<sub>s<\/sub><\/em> = 0.3 (Results from Zheng <em>et al.<\/em> 2024).[\/caption]","rendered":"<p>As previously mentioned, Brom\u2019s method for estimating pile head deflection requires estimating the coefficient of subgrade reaction <em>K<sub>h<\/sub><\/em> 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 <em>et al.<\/em> (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 <em>et al.<\/em>\u2019s method allows estimating the swaying <em>K<sub>hh<\/sub><\/em>, rocking <em>K<sub>rr<\/sub><\/em> and swaying-rocking <em>K<sub>rh <\/sub><\/em>= <em>K<sub>hr<\/sub><\/em> components of the head stiffness of a single floating pile embedded in elastic soil of finite thickness (<a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109a<\/a>).<\/p>\n<p>It is reminded that the swaying stiffness <em>K<sub>hh<\/sub><\/em> is the lateral load <em>H<sub>w<\/sub><\/em> that, when applied to the head of a pile results in head deflection <em>y<\/em>(0) = 1 (<a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>). The rocking stiffness <em>K<sub>rr<\/sub><\/em> is the bending moment <em>M<sub>w<\/sub><\/em> that, when applied to the head of a pile results in head rotation <em>\u03c9<\/em>(0) = 1 (<a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>). The swaying-rocking stiffness <em>K<sub>rh<\/sub><\/em> is the bending moment that will develop to a pile with fixed head rotation <em>\u03c9<\/em>(0) = 0 when the pile is subjected to unit displacement at its head <em>y<\/em>(0) = 1, considering the sign convention in <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>. The swaying-rocking stiffness <em>K<sub>hr<\/sub><\/em> is the lateral force that will develop to a pile with fixed head displacement <em>y<\/em>(0) = 0 when the pile is subjected to unit rotation at its head <em>\u03c9<\/em>(0) = 1, considering the sign convention in <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.109-hr.png\">Figure 6.109b<\/a>. Moreover, it is <em>K<sub>rh<\/sub><\/em> = <em>K<sub>hr<\/sub><\/em>.<\/p>\n<p>Using Zheng <em>et al.<\/em>\u2019s method to calculate the stiffness components allows establishing the stiffness [<strong>K<\/strong>] and compliance [<strong>S<\/strong>] matrixes of the pile, as:<\/p>\n<p><strong>(6.128a)<\/strong> [latex]\\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][\/latex]<\/p>\n<p><strong>(6.128b) <\/strong>[latex]\\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][\/latex]<\/p>\n<p>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 <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.110-hr.png\">Figure 6.110<\/a> allow calculating directly the components of the stiffness matrix for soil with Poisson\u2019s ratio <em>v<sub>s<\/sub><\/em> = 0.3 and different values of the pile-soil moduli ratio <em>E<sub>p<\/sub><\/em>\/<em>E<sub>s<\/sub><\/em>, the slenderness of the pile <em>L<\/em>\/<em>D<\/em> and the dimensionless thickness of the foundation soil layer <em>H<\/em>\/<em>L<\/em>. Having estimated these components one can calculate the compliance matrix by inverting the stiffness matrix, and subsequently the deflection <em>y<\/em>(0) and the rotation <em>\u03c9<\/em>(0) of the pile head for any load combination acting at the pile head, as well as head fixity conditions.<\/p>\n<p>It is interesting to observe in <a href=\"https:\/\/oercollective.caul.edu.au\/app\/uploads\/sites\/143\/2025\/03\/6.110-hr.png\">Figure 6.110<\/a> that there exists a critical pile slenderness ratio <em>L<\/em>\/<em>D<\/em>, 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\u00a0extrapolated 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 <em>E<sub>p<\/sub><\/em>\/<em>E<sub>s<\/sub><\/em>.<\/p>\n<figure id=\"attachment_623\" aria-describedby=\"caption-attachment-623\" style=\"width: 800px\" class=\"wp-caption aligncenter\"><img decoding=\"async\" class=\"wp-image-623 size-full\" src=\"https:\/\/libraryresources.nse.org.ng\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341.png\" alt=\"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&lt;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.\" width=\"800\" height=\"482\" srcset=\"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341.png 800w, https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341-300x181.png 300w, https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341-768x463.png 768w, https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341-65x39.png 65w, https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341-225x136.png 225w, https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-content\/uploads\/sites\/9\/2025\/03\/6.109-hr-e1743555074341-350x211.png 350w\" sizes=\"(max-width: 800px) 100vw, 800px\" \/><figcaption id=\"caption-attachment-623\" class=\"wp-caption-text\">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 <em>et al.<\/em> 2024).<\/figcaption><\/figure>\n<figure id=\"attachment_623\" aria-describedby=\"caption-attachment-623\" style=\"width: 936px\" class=\"wp-caption aligncenter\"><img decoding=\"async\" class=\"wp-image-2830 size-large\" src=\"https:\/\/libraryresources.nse.org.ng\/wp-content\/uploads\/sites\/9\/2026\/03\/6.110-hr-scaled.png\" alt=\"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.\" width=\"936\" height=\"1024\" \/><figcaption id=\"caption-attachment-623\" class=\"wp-caption-text\">Figure 6.110. Variation of pile head stiffness components with pile-soil moduli ratio <em>E<sub>p<\/sub>\/E<sub>s<\/sub><\/em>, pile slenderness <em>L\/D<\/em> and dimensionless thickness of the foundation soil layer <em>H\/L<\/em>. Results for soil Poisson\u2019s ratio <em>v<sub>s<\/sub><\/em> = 0.3 (Results from Zheng <em>et al.<\/em> 2024).<\/figcaption><\/figure>\n","protected":false},"author":1,"menu_order":39,"template":"","meta":{"pb_show_title":"","pb_short_title":"6.26 Rigorous elasticity method for estimating pile deformations","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-625","chapter","type-chapter","status-publish","hentry"],"part":421,"_links":{"self":[{"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/pressbooks\/v2\/chapters\/625","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/wp\/v2\/users\/1"}],"version-history":[{"count":1,"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/pressbooks\/v2\/chapters\/625\/revisions"}],"predecessor-version":[{"id":626,"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/pressbooks\/v2\/chapters\/625\/revisions\/626"}],"part":[{"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/pressbooks\/v2\/parts\/421"}],"metadata":[{"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/pressbooks\/v2\/chapters\/625\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/wp\/v2\/media?parent=625"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/pressbooks\/v2\/chapter-type?post=625"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/wp\/v2\/contributor?post=625"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/wp-json\/wp\/v2\/license?post=625"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}