Acessibilidade / Reportar erro

Hankel transformation method for solving the Westergaard problem for point, line and distributed loads on elastic half-space

Abstract

The Hankel transformation method was used in this work to determine the normal and shear stress distributions due to point, line and distributed loads applied to the surface of an elastic media. The elastic media considered in this study was assumed to be inextensible in the horizontal directions, and the only non-vanishing displacement component is the vertical component. Such materials were first considered by Westergaard as models of elastic half-space with alternating layers of soft and stiff materials with the stiff materials of negligible thickness and so closely spaced that the composite material characteristics is idealised as isotropic and homogeneous. The study adopted a displacement formulation. The Hankel transformation was applied to the governing Cauchy - Navier differential equation of equilibrium to reduce the problem to a second order ordinary differential equation (ODE) in terms of the deflection function in the Hankel transform space. Solution of the ODE subject to the boundedness condition yielded bounded deflection function in the Hankel transform space. Equilibrium of the internal vertical forces and the external applied load, was used to obtain the constant of integration, for the deflection function in the Hankel transform space. Inversion yielded the deflection in the physical domain variable. The stress displacement equations were then used to determine the Cauchy stresses. The vertical stress distributions due to line and distributed load over a rectangular area were also determined by using the point load solution for vertical stresses as Green functions, and then performing integration along the line and over the rectangular area of the load. The results obtained for vertical stresses due to point, line and distributed loads were determined in terms of dimensionless influence coefficients which were presented. The results obtained for the deflection, normal and shear stresses due to point, line and distributed load agreed with the solutions originally presented by Westergaard who used a stress function method.

Keywords
Hankel transform method; Westergaard problem; Cauchy-Navier displacement equation of equilibrium; elastic half-space

1. INTRODUCTION

1.1 Background

The determination of the normal and shear stress distributions due to point, line and distributed loads applied to the surface of linear elastic media are problems of the classical mathematical theory of elasticity (Ike, 2006Ike C.C. (2006). Principles of Soil Mechanics. De-Adroit Innovation, Enugu., 2018aIke C.C. (2018a). Hankel transform method for solving axisymmetric elasticity problems of circular foundation on semi-infinite soils. International Journal of Engineering and Technology (IJET) Vol 10. No 2 April - May 2018 pp 549 - 564. DOI: 10.21817/ijet/2018/v10 i2/181002111.
https://doi.org/10.21817/ijet/2018/v10 i...
, 2018bIke C.C. (2018b). General Solutions for Axisymmetric Elasticity Problems of Elastic Half Space using Hankel Transform Method. International Journal of Engineering and Technology (IJET) Vol 10. No 2 April - May 2018 pp 565 - 580. DOI: 10.21817/ijet/2018/v10i2/181002112.
https://doi.org/10.21817/ijet/2018/v10i2...
). Such problems are extensively encountered and applied in the elastic stress and settlement analysis of solids and structural footings (Rocscience, 2018Rocscience (2018). Westergaard stress solution method. https://www.rocscience.com/ documents/../Westergaard_Stress_Solution_Method pdf. Retrieved 02/09/2018.
https://www.rocscience.com/ documents/.....
; Westergaard, 1938Westergaard H.M. (1938). A problem of elasticity suggested by a problem in soil mechanics: soft soil reinforced by numerous strong horizontal sheets Contributions to the mechanics of solids. Stephen Timoshenko 60th Birthday Anniversary Volume Macmillan New York.). The foundational principles for the formulation of such problems, as indeed all elasticity problems are the simultaneous consideration of a system of fifteen differential equations in fifteen unknowns, namely: three differential equations of equilibrium, six equations representing the generalized Hooke’s stress - strain relations and six equations expressing strains in terms of the displacement components. The system of fifteen equations are required to additionally satisfy the compatibility conditions as well as the traction and deformation boundary conditions. The governing equations are more complex when anisotropy and heterogeneity of the three dimensional elastic media as well as the non-linear considerations of large deformations are taken into account. However, even for the cases of simplifications in the material properties of the media introduced by the assumptions of isotropy and homogeneity and the consideration of small deformations, the governing equations of the mathematical theory of elasticity are unwieldy in their numbers and demand extensively intensive analytical demanding techniques.

Three methods used in the formulation of the mathematical theory of elasticity problems are the displacement method, stress method and mixed (hybrid) methods (Ojedokun and Olutoge, 2012Ojedokun O.Y. and Olutoge F.A. (2012). Application of Boussinesq’s and Westergaard’s formulae in analysing foundation stress distribution for a failed telecommunication mast. African Journal of Mathematics and Computer Science Research Vol 5(4) pp 71 - 77 15 February 2012. Available on line at http://www.academicjournals.org/AJMCSR DOI:10.5897/AJMCSR 11.155 ISSN:2006-9731 Academic Journals.
http://www.academicjournals.org/AJMCSR D...
; Westergaard, 1964Westergaard H.M. (1964). Theory of Elasticity and Plasticity. Dover New York.; Barbar, 2010Barbar J.R. (2010). Elasticity Third Revised Edition. Springer Science and Bussiness Media. Dordrecht The Netherlands. DOI: https/doi.org/10.007/978-90-481-38098-4.
https/doi.org/10.007/978-90-481-38098-4...
; Sadd, 2014Sadd M.H. (2014). Elasticity Theory, Application and Numerics, Third Edition. University of Rhode Island. Elsevier Academic Press Amsterdam https://personnel egr.uri/edu/teggart/ courses/ce571/chapter 7. ppx.
https://personnel egr.uri/edu/teggart/ c...
; Bowles, 1997Bowles J.E. (1997). Foundation Analysis and Design, International Edition. McGraw Hill International Book Company, Tokyo (1997) ISBN 0-07-912247-7.). The displacement methods of the theory of elasticity presented by Navier and Lamé involve a reformulation of the system of fifteen governing equations of equilibrium, stress - strain and kinematics such that only the three components of displacement (in a three dimensional problem) become the primary unknown variables (Kachanov et al., 2003Kachanov M.L., Shafiro B. and Tsukrov I. (2003). Handbook of Elasticity Solutions. Springer Science and Bussiness Media Kluwer Academic Publishers Dordrecht The Netherlands.; Sitharam and Govinda Reju, 2017Sitharam T.G. and Govinda Reju L. (2017). Applied Elasticity for Engineers Module: Elastic Solutions with Applications in Geomechanics 14.139.172.204/nptel/1/CSE/web/ 105108070/module 8/lecture 17.pdf.; Palaniappan, 2011Palaniappan D. (2011). A general solution of equations of equilibrium in linear elasticity. Applied Mathematical Modelling 35 (2011) pp 5494 - 5499 Elsevier.; Hazel, 2015Hazel A. MATH 35021 Elasticity www.maths.manchester.ac.lc/_ahazel/MATHS. Nov 30, 2015.
www.maths.manchester.ac.lc/_ahazel/MATHS...
). The advantage of the displacement formulation is the obvious reduction in the number of equations to be solved from fifteen to three coupled equations in terms of the three displacement components. The other unknowns Cauchy stresses and strains are then found using the strain displacement, and stress - strain relations after the Navier - Lamé equations are solved for the Cartesian displacement components.

Stress methods, introduced by Michell and Beltrami involve the reformulation of the system of fifteen governing equations such that only the six Cauchy stress components (in a three dimensional formulation) become the primary variables (Ike, 2006Ike C.C. (2006). Principles of Soil Mechanics. De-Adroit Innovation, Enugu., 2018aIke C.C. (2018a). Hankel transform method for solving axisymmetric elasticity problems of circular foundation on semi-infinite soils. International Journal of Engineering and Technology (IJET) Vol 10. No 2 April - May 2018 pp 549 - 564. DOI: 10.21817/ijet/2018/v10 i2/181002111.
https://doi.org/10.21817/ijet/2018/v10 i...
, 2018bIke C.C. (2018b). General Solutions for Axisymmetric Elasticity Problems of Elastic Half Space using Hankel Transform Method. International Journal of Engineering and Technology (IJET) Vol 10. No 2 April - May 2018 pp 565 - 580. DOI: 10.21817/ijet/2018/v10i2/181002112.
https://doi.org/10.21817/ijet/2018/v10i2...
; Rocscience, 2018Rocscience (2018). Westergaard stress solution method. https://www.rocscience.com/ documents/../Westergaard_Stress_Solution_Method pdf. Retrieved 02/09/2018.
https://www.rocscience.com/ documents/.....
; Westergaard, 1938Westergaard H.M. (1938). A problem of elasticity suggested by a problem in soil mechanics: soft soil reinforced by numerous strong horizontal sheets Contributions to the mechanics of solids. Stephen Timoshenko 60th Birthday Anniversary Volume Macmillan New York.; Ojedokun and Olutoge, 2012Ojedokun O.Y. and Olutoge F.A. (2012). Application of Boussinesq’s and Westergaard’s formulae in analysing foundation stress distribution for a failed telecommunication mast. African Journal of Mathematics and Computer Science Research Vol 5(4) pp 71 - 77 15 February 2012. Available on line at http://www.academicjournals.org/AJMCSR DOI:10.5897/AJMCSR 11.155 ISSN:2006-9731 Academic Journals.
http://www.academicjournals.org/AJMCSR D...
; Westergaard, 1964; Barbar, 2010Barbar J.R. (2010). Elasticity Third Revised Edition. Springer Science and Bussiness Media. Dordrecht The Netherlands. DOI: https/doi.org/10.007/978-90-481-38098-4.
https/doi.org/10.007/978-90-481-38098-4...
; Sokolnikoff, 1956Sokolnikoff I.S. (1956). Mathematical Theory of Elasticity Second Edition. Data McGraw Hill Publishing Company Ltd Bombay New Delhi.; Timoshenko and Goodier, 1970Timoshenko S.P. and Goodier J.N. (1970). Theory of Elasticity Third Edition, McGraw Hill New York.; Green and Zerna, 1954Green A.E. and Zerna W. (1954). Theoretical Elasticity, Oxford University Press, London.). The advantage is the simplification inherent in the reduction of the number of governing equations from fifteen to six. The system of coupled Beltrami - Michell equations in terms of stresses are then solved to find the six stress components. Thereafter, the strains are found using the stress - strain relations, and the displacements obtained by integration of the strain - displacement equations subject to the displacement boundary conditions.

The mixed (hybrid) method, which is not commonly used in the theory of elasticity involve a reformulation of the governing equations in terms of some Cauchy stress components and some displacement components as the primary unknown variables.

Many researchers have derived some mathematical solutions to the Navier - Lamé displacement based formulations and the Beltrami - Michell stress - based formulations of the mathematical theory of elasticity; and presented stress and displacement functions that are solutions of the Beltrami - Michell stress and Navier - Lamé displacement equations (Nwoji et al., 2017aNwoji C.U., Onah H.N., Mama B.O., Ike C.C. (2017a). Solution of elastic half space problem using Boussinesq displacement potential functions. Asian Journal of Applied Sciences (AJAS) Vol 5, No 5 October 2017 pp. 1100 - 1106., 2017bNwoji C.U., Onah H.N., Mama B.O., Ike C.C. (2017b). Solution of the Boussinesq problem using Green and Zerna displacement potential function method. The Electronic Journal of Geotechnical Engineering (EJGE) 2017 Volume 22 Bundle 11 (22.11) pp. 4305 - 4314. Available at ejge.com.; Ike, 2017Ike C.C. (2017). First principles derivation of a stress function for axially symmetric elasticity problems, and application to Boussinesq problem. Nigerian Journal of Technology (NIJOTECH) Vol 36 No 3 July 2017 pp 767 - 772. http://dx-doi/10.4314/nijt/v.36i3.15. (www.nijotech.com).
http://dx-doi/10.4314/nijt/v.36i3.15...
; Ike et al., 2017aIke C.C., Onah H.N. and Nwoji C.U. (2017a). Bessel functions for axisymmetric elasticity problems of the elastic half space soil, a potential function method. Nigerian Journal of Technology (NIJOTECH) Vol 36 No 3 July 2017 pp. 773 - 781. http://dx.doi.org/10.4314/nijt/v.36i3.16.
http://dx.doi.org/10.4314/nijt/v.36i3.16...
, 2017bIke C.C., Mama B.O., Onah H.N. and Nwoji C.U. (2017b). Trefftz Harmonic function method for solving Boussinesq problem. Electronic Journal of Geotechnical Engineering (EJGE) 2017 (22.12) pp 4589 - 4601. Available at ejge.com.; Apostol, 2017Apostol B.F. (2017). Elastic Displacement in Half Space under the Action of a Tension force, General Solution for the Half Space with Point Forces. Journal of Elasticity Vol 126 pp 231 - 244. DOI:10.1007/S10659-016-9592-3.
https://doi.org/10.1007/S10659-016-9592-...
; Lurie and Vasilev, 1995Lurie S.A. and Vasilev V.V. (1995). The Biharmonic Problem in the Theory of Elasticity Gordon and Breach Publishers United States.; Chan, 2013Chan K.T. (2013). Analytical methods in Geomechanics CRC Press Taylor and Francis Group New York.; Teodorescu, 2013Teodorescu P.P. (2013). Treatise on Classical Elasticity Theory and Related Problems. Mathematical and Analytical Techniques with Applications to Engineering. Springer Dordrecht. DOI:10.1007/978-94-007-2616-1.
https://doi.org/10.1007/978-94-007-2616-...
; Abeyartne, 2012Abeyartne R. (2012). Continuum Mechanics, Volume 11 of Lecture notes on the mechanics of Elastic Solids. Cambridge. http://web.mit.edu/abeyartne/lecture notes/html 11th may 2012. Retrieved 8th April 2016.
http://web.mit.edu/abeyartne/lecture not...
; Davis and Salvadurai, 1996Davis R.O. and Salvadurai A.P.S. (1996). Elasticity and Goemechanics. Cambridge University Press.; Padio-Guidugli and Favata, 2014Padio-Guidugli P. and Favata A. (2014). Elasticity for Geotechnicians. A modern exposition of Kelvin, Boussinesq, Flammant, Cerrutti, Melan and Mindlin problems. Solid mechanics and its applications. Springer.). Some of the stress functions are: Boussinesq - Papkovich stress functions, Love stress functions, Airy stress functions, Morera stress function and Maxwell stress function. Displacement potential functions in the literature are: Boussinesq displacement potential functions, Papkovich - Neuber potential functions, Trefftz functions, Green and Zerna functions and Cerrutti functions.

The problem of point load applied at the origin (0, 0, 0) of a homogeneous, semi-infinite, linear elastic half-space is the Boussinesq problem, found in the theory of elasticity applied to soil/solid mechanics. The Boussinesq problem is a fundamental problem whose solutions are fundamental Green functions used in deriving solutions for distributed loads acting on given lines or given areas on the boundary of the semi-infinite linear elastic soil/solid media.

The Boussinesq stress distribution theory is very extensively used in the determination of the distributions of stresses in elastic media considered semi-infinite and three dimensional. The theory assumes that the semi-infinite medium is linear elastic, homogeneous and isotropic. However, most soils are not isotropic neither are they homogeneous.

Westergaard sought to solve the theory of elasticity problem of finding stress distributions at any point under a foundation that would cater for alternating horizontal layers of soft and stiff materials (Ike, 2006Ike C.C. (2006). Principles of Soil Mechanics. De-Adroit Innovation, Enugu., 2018a, 2018b; Rocscience, 2018; Anyaegbunam, et al., 2011Anyaegbunam A.J., Osadebe N.N., Eze-Uzoamaka O.J. (2011). Non-Existence of Solution for Horizontally Rigid Half-Space ASCE, Journal of Geotechnical and Geoenvironmental Engineering 137(4) April, pp 431 - 434. http://ascelibrary.org/gto/ resource/1/jgefk/v.137/i4.doi.10.1061/(ASCE)GT 1943- 5606.0000444.
http://ascelibrary.org/gto/ resource/1/j...
). For such materials with alternating soft and stiff layers, the Boussinesq stress distribution theory cannot correctly give the stresses, since the fundamental assumptions for constructing the Boussinesq theory have been violated by the material properties of the half-space.

Westergaard presented solutions for the stress fields for such problems by making the following fundamental assumptions (Ike, 2006Ike C.C. (2006). Principles of Soil Mechanics. De-Adroit Innovation, Enugu., 2018a, 2018b; Rocscience, 2018):

  • (i) The stiff layers are very closely spaced such that the composite material characteristics could be approximated or idealized as isotropic and homogeneous.

  • (ii) The stiff layers are of negligible thickness as compared to the thickness of the soft materials.

  • (iii) The stiff layer is inextensible, in the horizontal directions. Thus the horizontal components of displacement ux=uy=0 for both the x and the y - Cartesian coordinate directions.

  • (iv) (iv) The inextensible stiff layer restrains any horizontal strains in both the stiff and the soft layers. The horizontal normal strains in the x and y Cartesian coordinate directions are zero, (εxx=εyy=0). The horizontal shear strain γxy is zero (Anyaegbunam et al., 2011Anyaegbunam A.J., Osadebe N.N., Eze-Uzoamaka O.J. (2011). Non-Existence of Solution for Horizontally Rigid Half-Space ASCE, Journal of Geotechnical and Geoenvironmental Engineering 137(4) April, pp 431 - 434. http://ascelibrary.org/gto/ resource/1/jgefk/v.137/i4.doi.10.1061/(ASCE)GT 1943- 5606.0000444.
    http://ascelibrary.org/gto/ resource/1/j...
    ).

Generally, sedimentary soils such as natural clay strata highlight the anisotropic characteristics of soil masses (Bowles, 1997Bowles J.E. (1997). Foundation Analysis and Design, International Edition. McGraw Hill International Book Company, Tokyo (1997) ISBN 0-07-912247-7.). Consequently, the Westergaard theory is a better representation of the elasticity problem of three dimensional soil media (Rocscience, 2018Rocscience (2018). Westergaard stress solution method. https://www.rocscience.com/ documents/../Westergaard_Stress_Solution_Method pdf. Retrieved 02/09/2018.
https://www.rocscience.com/ documents/.....
). However many geotechnical engineers prefer the use of Boussinesq theory to model elastic soils primarily because it gives higher values of vertical stresses and could be considered more conservative from the perspective of design.

1.2 The Hankel transformation method and why it is used/adopted in this study

Hankel transformations are integral transformations whose kernels are Bessel functions, (Everitt and Kalf, 2007Everitt W.N. and Kalf H. (2007). The Bessel Differential Equation and the Hankel transform. Journal of Computational and Applied Mathematics, Elsevier Volume 208, Issue 1, November 2007, pp 3 - 19. https://doi.org/10.1016/j.cam.2006.10.029.; Duffy, 2004Duffy D.G. (2004). Transform Methods for Solving Partial Differential Equations Second Edition. Chapman and Hall/CRC ISBN 9781584884514.; Negero, 2016Negero N.T. (2016). Zero Order Hankel Transform Method for Partial Differential Equations. International Journal of Modern Sciences and Engineering Technology (IJMSET) Volume 3 Issue 10, 2016 pp 24 - 36. ISSN 2349 - 3755. Available at www.ijmset.com). They can also be referred to as Fourier - Bessel transforms. Hankel transforms are integral transforms applied to physical domains represented by infinite two dimensional plane with polar/cylindrical coordinate (r, θ, z) where 0r, (Tuteja et al., 2014Tuteja R., Jaloree S., and Goyal A. (2014). Application of Hankel transform of I-function of one variable for solving axisymmetric Dirichlet potential problem. Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 14, Issue 4, Version 1.0 pg 11 - 16. Online ISSN: 2249-4626.; Negero, 2016; Malacka Zuzana, 2018Malacka Zuzana S.K. (2018). Hankel Transform and Free Vibration of a Large Circular Membrane. 17th Conference on Applied Mathematics APLIMAT 2018 Faculty of Mechanical Engineering, Slovak University of Technology in Bratislava pp 717 - 723. evim.stuba.sk/APLIMAT2018/proceedings/papers/0717_Malacka.pdf.). The radial coordinate term of the Laplacian is the object of the Hankel transform. The Hankel transforms arise naturally in solving boundary value problems formulated and expressed in terms of the cylindrical polar coordinates system (Tuteja et al., 2014; Andrews and Shivamoggi, 1999Andrews L.C. and Shivamoggi B.K. (1999). Integral Transforms for Engineers. SPIE Optical Engineering Press. A publication of SPIE - The International Society for Optical Engineering, Bellingham Washington USA ISBN: 0-8194-3232-6.; Voegtle, 2017Voegtle I.C. (2017). The Bessel Function, the Hankel Transform and an Application to Differential Equations. Electronic Thesis and Dissertations Jack N. Averity College of Graduate Studies, Georgia Southern University Summer 2017. https://digitalcommons.georgiasouthern.edu/etd.
https://digitalcommons.georgiasouthern.e...
).

They are thus particularly useful in dealing with boundary value problems that are characterised by circular symmetry (Piessens, 2000Piessens R. (2000). “The Hankel Transform” In “The Transforms and Applications Handbook: Second Edition” Ed. Alexander D. Poularikas Boca Raton: CRC Press LLC.; Yokoyama, 2014Yokoyama Takehito (2014). The Hankel transform. www.stat.phys.titech.ac.jp/~yokayama /memo4.pdf. Accessed 1/7/2018.
www.stat.phys.titech.ac.jp/~yokayama /me...
). The Hankel transform of order, v equal to zero is the two-dimensional Fourier transform of a circularly symmetric function of three Cartesian coordinate variables. The most important special cases of the Hankel transform correspond to the order zero and one (i.e. v = 0, and v = 1).

The conditions for the existence of the Hankel transformation for a function and its derivatives are similar to the conditions for the existence of the Fourier transform for the function and its derivatives, and these conditions are usually satisfied for the functions that/which describe most physical and engineering models. The Hankel transform and its inverse are linear operators.

Similar to the Laplace and Fourier transforms, it is required that for the function being transformed in the Hankel transform space, the function and its derivatives with respect to the radial coordinate variable must vanish as r approaches infinity (r → ∞). (Tuteja et al., 2014Tuteja R., Jaloree S., and Goyal A. (2014). Application of Hankel transform of I-function of one variable for solving axisymmetric Dirichlet potential problem. Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 14, Issue 4, Version 1.0 pg 11 - 16. Online ISSN: 2249-4626.; Malacka Zuzana, 2018Malacka Zuzana S.K. (2018). Hankel Transform and Free Vibration of a Large Circular Membrane. 17th Conference on Applied Mathematics APLIMAT 2018 Faculty of Mechanical Engineering, Slovak University of Technology in Bratislava pp 717 - 723. evim.stuba.sk/APLIMAT2018/proceedings/papers/0717_Malacka.pdf.; Negero, 2016Negero N.T. (2016). Zero Order Hankel Transform Method for Partial Differential Equations. International Journal of Modern Sciences and Engineering Technology (IJMSET) Volume 3 Issue 10, 2016 pp 24 - 36. ISSN 2349 - 3755. Available at www.ijmset.com; Andrews and Shivamoggi, 1999Andrews L.C. and Shivamoggi B.K. (1999). Integral Transforms for Engineers. SPIE Optical Engineering Press. A publication of SPIE - The International Society for Optical Engineering, Bellingham Washington USA ISBN: 0-8194-3232-6.).

In the elastic half-space problem, the stresses and displacements are required to vanish as r → ∞, and this satisfies a key requirement for the application of the Hankel transformation. Again elastic half-space problems involving point load acting at the origin are axisymmetric and the stresses and displacement fields are expected to show circular symmetry about the axis of symmetry which is the axis of the application of the point load. Being boundary value problems with circular symmetry, the elastic half-space problem can be solved using the Hankel transform method.

Liao and Wang (1998Liao J.J. and Wang C.D. (1998). Elastic solutions for a transversely isotropic half-space subjected to a point load. International Journal for Numerical and Analytical Methods in Geomechanics (Int. J. Numer. Anal. Meth. Geomech) 22 pp. 425 - 447.) derived and presented closed-form solutions for the stresses and displacement fields in transversely isotropic, elastic half-space under the action of a point load. They used the Fourier and Hankel transforms in a cylindrical coordinate system in deriving the closed-form solutions. They assumed the plane of transverse isotropy to be parallel to the horizontal surface, and the half-space to be bounded by the horizontal surface (z = 0). They obtained their closed-form solution by using the superposition of two infinite spaces; one with a point load acting in the interior, and the other being free loading. Their solutions were identical with the Mindlin and Boussinesq solutions if the transversely isotropic half-space becomes homogeneous, linear elastic and isotropic.

Tarn and Wang (1987Tarn J.Q. and Wang Y.M. (1987). A fundamental solution for a transversely isotropic elastic space. J. Chin. Inst. Eng Volume 10 Issue 1. pp 13 - 21.) used the Fourier and Hankel transforms to solve the three dimensional (3D) elasticity problem of transversely isotropic, infinite space, and obtained solutions which reduced to the Kelvin solution for the case of isotropic media. Closed-form solutions for the stresses and displacement fields in transversely isotropic media due to the action of a point load have also been presented by Barden (1963Barden L. (1963). Stresses and displacements in a cross-anisotropic soil. Geotechnique 13 pp. 198 - 210.), and Misra and Sen (1975Misra B. and Sen B.R. (1975). Stresses and displacements in granular materials due to surface load. Int. J. Engng Sci. 13 pp. 743 - 761.).

The Hankel transform method was used by Ike (2018aIke C.C. (2018a). Hankel transform method for solving axisymmetric elasticity problems of circular foundation on semi-infinite soils. International Journal of Engineering and Technology (IJET) Vol 10. No 2 April - May 2018 pp 549 - 564. DOI: 10.21817/ijet/2018/v10 i2/181002111.
https://doi.org/10.21817/ijet/2018/v10 i...
) in a stress-based formulation to obtain general solutions for stress and displacement fields in homogeneous, isotropic, linear elastic semi-infinite soil subject to uniformly distributed axisymmetric load over a circular area on the surface. Ike (2018a) applied the Hankel transformation to the biharmonic stress compatibility equation to obtain bounded stress functions for the elastic half-space problem. He also applied Hankel transform to the Love stress functions to obtain the stresses and displacements in the Hankel transform space; and enforced the boundary conditions to determine the unknown constants of the Love stress function in the Hankel transform space. By inversion of the Hankel transform expressions, the stresses and displacements for the elastic half-space under uniformly distributed load on a circular area were found by Ike (2018a) in the physical domain variables. His solutions were identical with solutions in literature that used Boussinesq stress functions. Ike (2018b) used the Hankel transform method to obtain general solutions for stresses and displacements in semi-infinite, linear elastic isotropic media under axisymmetric load. He used Hankel transformation of the governing equations of the 3D axisymmetric elasticity problem in a stress-based formulation to obtain the Love stress function in the Hankel transform space. Ike (2018b) also applied the Hankel transformation to the Love stresses and displacement fields to obtain general solutions for the stresses and displacements in the Hankel transform space. The general solutions obtained were used to solve the particular axisymmetric elasticity problem originally solved using Boussinesq stress potential functions by Boussinesq. Ike (2018b) obtained the same expressions for stresses and displacements fields as obtained by Boussinesq even with the Hankel transform method in a Love stress function formulation.

1.3 Research aim and objectives

The research aim is to use the Hankel transformation method to solve the Westergaard problem for point, line and distributed loads on the surface of linear elastic, isotropic half-space. The specific objectives include:

  • (i) to use the Hankel transformation method to obtain bounded solutions to the Cauchy - Navier displacement equations of equilibrium for the Westergaard problem for the case of point load Q 0 acting at the origin of the semi-infinite elastic space.

  • (ii) to use the condition of equilibrium of the internal vertical normal stresses and vertical point load to determine the unknown constant of integration in the integration problem of the Cauchy - Navier displacement equation of equilibrium.

  • (iii) to apply the Hankel transformation to the vertical stress - displacement equation and obtain the vertical normal stresses σ¯zz(β,z) in the Hankel transform space variable, β.

  • (iv) to obtain the vertical stress distribution w(r, z) in the physical domain space variables by Hankel inverse transformation of the solution obtained in the Hankel transform space.

  • (v) to obtain the vertical normal stresses for the Westergaard point load problem in the physical domain variables using the vertical normal stress - displacement relations; and express the stresses in terms of non-dimensional Westergaard vertical stress influence coefficients.

  • (vi) to use the solution for a point load as Green function and obtain by superposition theory the solution for vertical stresses in the elastic half-space due to uniformly distributed line load of intensity q 0 and finite length.

  • (vii) to use the point load solution as Green’s function and obtain by superposition principles, the vertical stress distribution at any arbitrary depth under a corner of rectangular foundation areas subject to uniformly distributed load of intensity, p 0.

2. THEORETICAL FRAMEWORK

In the Westergaard problem, the horizontal components of the displacement field are assumed to be zero, and the only non-vanishing displacement component is the vertical displacement. The displacement field components in terms of the three dimensional Cartesian coordinates x, y, z components are

u x = u = 0 (1)

u y = v = 0 (2)

u z = w ( x , y , z ) (3)

where ux=u, uy=v, and uz=w(x,y,z) are respectively the x, y, and z Cartesian coordinate components of the displacement.

Cauchy - Navier displacement formulation of the differential equation of equilibrium in the absence of body force then reduce to the partial differential equation (PDE):

G 2 w ( x , y , z ) + ( λ + 2 G ) 2 w z 2 ( x , y , z ) = 0 (4)

where ∇2 is the Laplacian operator, given by:

2 = 2 x 2 + 2 y 2 (5)

G is the shear modulus of elasticity given in terms of the Young’s modulus of elasticity E, and the Poisson’s ratio, μ as:

G = E 2 ( 1 + μ ) (6)

and λ is the Lamé’s constant (parameter or modulus) given by:

λ = μ E ( 1 + μ ) ( 1 2 μ ) = 2 μ G ( 1 2 μ ) (7)

Thus, Navier displacement equation of equilibrium can be expressed as:

G ( 2 w x 2 + 2 w y 2 ) + ( λ + 2 G ) 2 w z 2 = 0 (8)

The strain - displacement equations from the small deformation assumptions of the three dimensional (3D) theory of elasticity become, for the Westergaard problem:

ε x x = u x = 0 (9)

ε y y = v y = 0 (10)

ε z z = w z = 0 (11)

γ x y = u y + v x = 0 (12)

γ x z = u z + w x = w x (13)

γ y z = v z + w y = w y (14)

where εxx , εyy and εzz are normal strains, γxy , γxz and γyz are shear strains.

The stress - displacement relations for the Westergaard problem are found form the stress - strain laws using the strain - displacement relations as:

σ z z = ( λ + 2 G ) ε z z = ( λ + 2 G ) w z (15)

τ x z = σ x z = G w x (16)

τ y z = σ y z = G w y (17)

where σzz is the vertical normal stress τxz and τyz are shear stresses.

3. METHODOLOGY

3.1 Hankel transformation of the governing equations

For the Westergaard problem involving a point load Q 0 the cylindrical polar coordinate system can be adopted to express the Cauchy - Navier displacement equation of equilibrium as:

G ( r 2 w ) + ( λ + 2 G ) 2 w z 2 = 0 (18)

where

r 2 = 2 r 2 + 1 r r (19)

Thus,

G ( 2 w r 2 + 1 r w r ) + ( λ + 2 G ) 2 w z 2 = 0 (20)

or,

G λ + 2 G ( 2 w r 2 + 1 r w r ) + 2 w z 2 = 0 (21)

where

G λ + 2 G = 1 2 μ 2 ( 1 μ ) = α ( μ ) (22)

where α(μ) is a parameter, defined to depend on the Poisson’s ratio, μ.

Applying the Hankel transformation to the Cauchy-Navier differential equation of equilibrium, we have:

0 { G λ + 2 G ( 2 w r 2 + 1 r w r ) + 2 w z 2 } r J 0 ( β r ) d r = 0 (23)

where rJ 0r) is the kernel (nucleus) of the Hankel transformation, β is the Hankel transformation parameter, and J 0r) is the Bessel function of the first kind of order zero.

Using the linearity property of the Hankel transformation, we have:

0 G λ + 2 G ( 2 w r 2 + 1 r w r ) r J 0 ( β r ) d r + 0 2 w z 2 r J 0 ( β r ) d r = 0 (24)

Simplifying,

d 2 d z 2 0 w ( r , z ) r J 0 ( β r ) d r + G λ + 2 G 0 ( 2 w r 2 + 1 r w r ) r J 0 ( β r ) d r = 0 (25)

d 2 d z 2 0 r w ( r , z ) r J 0 ( β r ) d r β 2 G λ + 2 G 0 r w ( r , z ) J 0 ( β r ) d r = 0 (26)

Let

W ¯ ( β , z ) = 0 r w ( r , z ) J 0 ( β r ) d r (27)

W¯(β,z) is the Hankel transform of w(r, z).

Then,

d 2 W ¯ ( β , z ) d z 2 G λ + 2 G β 2 W ¯ β ( z ) = 0 (28)

d 2 W ¯ ( β , z ) d z 2 α 2 β 2 W ¯ ( β , z ) = 0 (29)

Using differential (D) operator methods or the method of trial functions, the general solution for W¯(β,z) is:

W ¯ ( β , z ) = c 1 ( β ) exp ( α β z ) + c 2 ( β ) exp ( α β z ) (30)

where c 1(β) and c 2(β) are the two integration constants.

For bounded solutions, W¯(β,z)is required to be finite and defined as z → ∞.

Therefore,

c 2 ( β ) = 0 (31)

Then for bounded solution,

W ¯ ( β , z ) = c 1 ( β ) exp ( α β z ) (32)

The constant of integration c1(β) is found from the boundary condition.

3.2 Enforcement of boundary condition

For a point load Q0 acting at the origin of the half-space as shown in Figure 1, the Hankel transformation of the load is q¯(β) where:

q ¯ ( β ) = 0 Q 0 δ ( r = 0, z = 0 ) 2 π r r J 0 ( β r ) d r (33)

where in Equation (33), δ is the Dirac delta function. Evaluation of Equation (33) yields:

q ¯ ( β ) = Q 0 2 π (34)

Figure 1
Point load at the origin O of an elastic half space made of Westergaard material

The boundary condition is established from the equilibrium of internal vertical normal stresses and external forces as:

σ z z ( r = 0, z = 0 ) + Q 0 δ ( r = 0, z = 0 ) 2 π r = 0 (35)

σ z z | z = 0, r = 0 = Q 0 δ ( r = 0, z = 0 ) 2 π r (36)

In Hankel transform space,

σ ¯ z z ( β , z ) = ( λ + 2 G ) W ¯ ( β , z ) z = ( λ + 2 G ) z c 1 ( β ) exp ( α β z ) (37)

σ ¯ z z ( β , z ) = ( λ + 2 G ) c 1 ( β ) ( α β exp ( α β z ) ) = α β ( λ + 2 G ) c 1 ( β ) exp ( α β z ) (38)

Using the boundary condition,

σ ¯ z z ( β , z = 0 ) | = α β ( λ + 2 G ) c 1 ( β ) = Q 0 2 π (39)

c 1 ( β ) = Q 0 2 π α β ( λ + 2 G ) (40)

Then,

W ¯ ( β , z ) = Q 0 2 π α β ( λ + 2 G ) exp ( α β z ) = Q 0 2 π α β G α 2 exp ( α β z )

= Q 0 α 2 π β G exp ( α β z ) = Q 0 α exp ( α β z ) 2 π G β (41)

4. RESULTS

4.1 Vertical displacement in the physical (real) space coordinates

The uniformly distributed line load of intensity q0 and finite length is considered to act on the surface of the half-space as shown in Figure 2. By inversion, the vertical displacement is obtained in the domain space variables as:

w ( r , z ) = 0 β W ¯ ( β , z ) J 0 ( β r ) d β = 0 β Q 0 2 π α β ( λ + 2 G ) exp ( α β z ) J 0 ( β r ) d β w ( r , z ) = Q 0 2 π α ( λ + 2 G ) 0 exp ( α β z ) J 0 ( β r ) d β = Q 0 α 2 π G 0 exp ( α β z ) J 0 ( β r ) d β (42)

w ( r , z ) = Q 0 α 2 π G 1 R (43)

where

R = ( x 2 + y 2 + α 2 z 2 ) 1 / 2 (44)

4.2 Vertical stresses in the physical domain variables

This is found from Equation (44) as:

σ z z ( r , z ) = ( λ + 2 G ) z { Q 0 2 π α ( λ + 2 G ) 0 exp ( α β z ) J 0 ( β r ) d β } (45)

σ z z ( r , z ) = ( λ + 2 G ) Q 0 2 π α ( λ + 2 G ) { z 0 exp ( α β z ) J 0 ( β r ) d β } (46)

σ z z ( r , z ) = Q 0 2 π α 0 α β exp ( α β z ) J 0 ( β r ) d β (47)

σ z z ( r , z ) = Q 0 2 π 0 β exp ( α β z ) J 0 ( β r ) d β (48)

σ z z ( r , z ) = Q 0 2 π α z ( r 2 + α 2 z 2 ) 3 / 2 (49)

σ z z ( r , z ) = Q 0 2 π z 2 α z 3 ( r 2 + α 2 z 2 ) 3 / 2 (50)

σ z z ( r , z ) = Q 0 z 2 1 2 π ( α z 2 ) 3 / 2 ( r 2 + α 2 z 2 ) 3 / 2 (51)

σ z z ( r , z ) = Q 0 z 2 α 2 π ( z 2 r 2 + α 2 z 2 ) 3 / 2 (52)

σ z z ( r , z ) = Q 0 z 2 α 2 π ( 1 r 2 / z 2 + α 2 ) 3 / 2 (53)

σ z z ( r , z ) = Q 0 z 2 α 2 π ( α 2 + r 2 z 2 ) 3 / 2 (54)

σ z z ( r , z ) = Q 0 z 2 I w ( r , z , α ( μ ) ) (55)

where

I w ( r , z , α ( μ ) ) = α 2 π ( α 2 + r 2 z 2 ) 3 / 2 (56)

Iw(r,z,α(μ)) is the Westergaard vertical stress influence coefficient at the point (r, z) in soil with Poisson’s ratio, μ.

But

r 2 = x 2 + y 2 (57)

I w ( x , y , z , α ( μ ) ) = α 2 π ( α 2 + x 2 + y 2 z 2 ) 3 / 2 (58)

The other normal stresses are computed from the stress - strain relation as:

σ x x = σ y y = μ 1 μ σ z z = k 0 σ z z = k 0 Q 0 z 2 I w ( r , z , α ( μ ) ) (59)

where k 0 is the coefficient of lateral stress at rest.

The shear stresses are found from the shear stress - displacement relations as:

τ x y = 0 = τ y x (60)

τ y z = τ z y = G w y = G y ( Q 0 α ( μ ) 2 π G R ) (61)

= G Q 0 2 π G y 1 R = G Q 0 2 π G y ( x 2 + y 2 + α 2 z 2 ) 1 / 2 (62)

τ z y = Q 0 α y 2 π R 3 (63)

Similarly,

τ z x = τ x z = G w x = G x ( Q 0 α ( μ ) 2 π G R ) (64)

τ x z = G Q 0 α 2 π G x ( x 2 + y 2 + α 2 z 2 ) 1 / 2 (65)

τ x z = Q 0 α x 2 π R 3 (66)

4.3 Uniformly distributed line load of intensity q 0 and finite length

Here, the vertical stress is found by using the solution for the point load as a Green function, as:

σ z z = 0 L 1 z 2 α 2 π ( α 2 + x 2 + y 2 z 2 ) 3 / 2 q 0 d y (67)

where

x = B (68)

Figure 2
Uniformly distributed line load of finite length acting on the surface of an elastic half-space

σ z z = q 0 α 2 π z 2 0 L ( α 2 + B 2 z 2 + y 2 z 2 ) 3 / 2 d y (69)

σ z z = q 0 z α 2 π ( n m 2 + α 2 1 ( m 2 + n 2 + α 2 ) 1 / 2 ) (70)

where

m = B z (71)

n = L z (72)

σ z z = q 0 z I W L ( m , n , α ( μ ) ) (73)

I W L ( m , n , α ( μ ) ) = α 2 π ( n m 2 + α 2 1 ( m 2 + n 2 + α 2 ) 1 / 2 ) (74)

IWL(m,n,α(μ)) is the non-dimensional vertical stress influence coefficient for uniformly distributed line load of intensity q 0 and finite length for a Westergaard elastic half-space.

4.4 Vertical stress distribution below a corner point of a rectangular area due to uniformly distributed load on a rectangular foundation ( B × L )

A uniformly distributed load is considered to act on a rectangular area B × L on the surface of the elastic half-space as shown in Figure 3. The vertical stress for the Westergaard problem is obtained in this case as the double integration problem over the loaded rectangular foundation area:

σ z z ( x , y , z ) = 0 L 0 B p 0 d x d y z 2 α 2 π ( α 2 + x 2 + y 2 z 2 ) 3 / 2 (75)

σ z z = p 0 0 L 0 B α 2 π z 2 ( α 2 + x 2 + y 2 z 2 ) 3 / 2 d x d y (76)

σ z z = p 0 2 π 0 L 0 B α z 2 ( α 2 + x 2 + y 2 z 2 ) 3 / 2 d x d y (77)

σ z z = p 0 2 π tan 1 ( L B α z 1 ( L 2 + B 2 + α 2 z 2 ) 1 / 2 ) (78)

σ z z = p 0 2 π tan 1 ( m 1 n 1 α ( m 1 2 + n 1 2 + α 2 ) 1 / 2 ) (79)

where

m 1 = L z (80)

n 1 = B z (81)

σ z z = p 0 I c w (82)

I c w = 1 2 π tan 1 ( m 1 n 1 α ( m 1 2 + n 1 2 + α 2 ) 1 / 2 ) (83)

Icw(m1,n1,α(μ)) is the non-dimensional vertical stress influence coefficient for uniformly distributed load on a rectangular foundation area for points at any depth z under a corner point in a Westergaard elastic half-space.

Figure 3
Uniformly distributed load of intensity p0 on a rectangular foundation (B × L)

For μ = 0,

α 2 ( μ = 0 ) = 1 2 μ 2 2 μ | μ = 0 = 1 2 = 0.5 (84)

α ( μ = 0 ) = 0.5 (85)

Then,

I w ( 0, z ; α = 0.5 ) = 0.5 2 π ( 0.5 + r 2 z 2 ) 3 / 2 (86)

Values of Iw(r,z,α=0.5) for various value of r/z can be calculated and presented in Tables. Table of values of the Westergaard vertical stress influence coefficient for μ = 0 for various values of r/z is shown in Table 1. The results agree with results obtained by Fadum (1948Fadum Ralph E. (1948). Influence values for estimating stresses in elastic foundations. Proceedings Second International Conference on Soil Mechanics and Foundation Engineering Vol 3 pp 77 - 84. https://www.issmge.org/uploads/ publications/1/43/ 194803-0020.pdf.
https://www.issmge.org/uploads/ publicat...
). The influence values are also calculated for various values of μ and r/z and presented in Table 2.

Table 1
Influence values for vertical, normal stresses due to point load at the origin (Hankel transform solution of the Westergaard problem) σzz(r,z)=QoIwz2(r,z,α(μ)) for μ=0

Table 2
Westergaard vertical stress influence coefficient for point load applied at the origin on the surface of an elastic half-space for various values of μ and (r/z)

5. DISCUSSION

The Hankel transformation method has been successfully used in this work to solve the Westergaard problem in the theory of elasticity involving finding vertical stresses due to a point load acting at the origin of an elastic half-space for which all the horizontal components of the displacement are assumed to vanish and the only non-vanishing displacement component is the vertical displacement. A displacement based approach was adopted in the study. Accordingly, the governing partial differential equation (PDE) of equilibrium is the Cauchy - Navier (PDE) given as Equation (4), which is solved by the Hankel transformation technique. Other considerations used in the solution were the strain displacement relations of the small - displacement assumptions of the 3D theory of elasticity given for the Westergaard problem as Equations (9 - 14); and the stress displacement equations given for the Westergaard problem as Equations (15 - 17).

Due to the symmetrical nature of the point load problem about the axis of application of the point load, the Cauchy - Navier displacement equation of equilibrium to be solved was expressed using the cylindrical polar coordinates system as Equation (18) or in expanded form, Equations (20) or (21). The Hankel transformation was applied to obtain Equation (23). The linearity property of the Hankel transformation was used to express the Hankel transformed Cauchy - Navier equation as the second order ordinary differential equation (ODE) given as Equation (28) in terms of W¯(β,z), the vertical displacement in the Hankel transform space. Equation (28) was solved using methods for solving ODEs to obtain the solution as Equation (30) which contained two integration constants. The requirements for the boundedness of the vertical displacement was used to obtain the integration constant c 2(β) as Equation (31); yielding the bounded solution for W¯(β,z) as Equation (32) which still contained one constant of integration. The equilibrium of internal vertical normal stresses and external load was used to obtain the boundary condition as Equation (36). Enforcement of the vertical stress boundary condition gives the unknown constant of integration c 1(β) as Equation (41); and the fully determined bounded solution for W¯(β,z)as Equation (42). By inversion of Equation (42), the vertical displacement in the physical domain space variables was found as Equation (44). The vertical stress distribution was obtained from Equation (44) using the vertical stress - displacement relation as Equation (48). Evaluation of the integral yielded the solution as Equation (49). The vertical stress was expressed in dimensionless form as Equation (54) or (55) where the Westergaard vertical stress influence coefficient is defined as Equation (56) or (58). The normal stress distribution in the horizontal directions σxx and σyy are found as Equation (59), while the shear stress distributions τxy , τzy , and τzx are found as Equations (60), (63) and (66).

The point load solution was used to express the vertical stresses due to uniformly distributed line load of intensity q 0 and finite length L as Equation (67). Evaluation of the integral gives the vertical stress for the distributed line load as Equation (70).

Similarly, the point load solution was used to express the vertical stress at points below the corner of rectangular foundations (B × L) carrying uniformly distributed load of intensity p 0 as Equation (67). The evaluation of Equation (75) gives the solution in dimensionless form as Equation (79) where the non-dimensional coefficient was defined as Equation (83).

For value of the Poisson ratio, μ = 0, the Westergaard vertical stress influence coefficient for point load was found as Equation (86); which was calculated and presented in tabulated form for values of r/z ranging from r/z =0 to r/z = 10 in Table 1. Values of the Westergaard vertical stress influence coefficients for various values of r/z and μ = 0, and μ = 0.40 are also shown in Table 2.

The deflection of points (x, y, 0) on the xy coordinate plane is given by substitution of z =0 into Equation (44) by:

w = ( x , y , z = 0 ) = Q 0 α 2 π G 1 ( x 2 + y 2 ) (87)

w = ( x , y , z = 0 ) = Q 0 α 2 π 2 ( 1 + μ ) E ( x 2 + y 2 ) = Q 0 α ( 1 + μ ) E π ( x 2 + y 2 ) (88)

However, at x = y = z =0 which is the origin and the point of application of the load, the deflection function becomes undefined and indeterminate, and this is expectedly due to the singular property of the point load acting at the origin, which introduces a singularity in the deflection function at the point of application (0,0,0). The normal and shear stress fields are similarly undefined and singular at the origin. However, on the xy coordinate plane at points (x, y, z = 0) the normal and shear stress fields are found by substitution of z = 0 into the expressions for the stresses. Thus,

σ z z ( x , y , z = 0 ) = σ z z ( r , z = 0 ) = (89)

σ x x ( x , y , z = 0 ) = σ x x ( r , z = 0 ) = (90)

σ y y ( x , y , z = 0 ) = σ y y ( r , z = 0 ) = (91)

τ x y ( x , y , z = 0 ) = τ x y ( r , z = 0 ) = 0 (92)

τ y z ( r , z = 0 ) = τ z y ( x , y , z = 0 ) = Q 0 α y 2 π ( x 2 + y 2 ) 3 / 2 = Q 0 α y 2 π r 3 (93)

τ z x ( x , y , z = 0 ) = τ z x ( r , z = 0 ) = Q 0 α x 2 π ( x 2 + y 2 ) 3 / 2 = Q 0 α x 2 π r 3 (94)

It is thus observed that at points on the xy Cartesian coordinate plane, the shear stresses τzy and τzx do not vanish but have definite values that depend on the coordinates of the point; and this is in line with observations made by Anyaegbunam et al. (2011Anyaegbunam A.J., Osadebe N.N., Eze-Uzoamaka O.J. (2011). Non-Existence of Solution for Horizontally Rigid Half-Space ASCE, Journal of Geotechnical and Geoenvironmental Engineering 137(4) April, pp 431 - 434. http://ascelibrary.org/gto/ resource/1/jgefk/v.137/i4.doi.10.1061/(ASCE)GT 1943- 5606.0000444.
http://ascelibrary.org/gto/ resource/1/j...
). However, the non-vanishing of the shear stresses τzy and τzx on the xy Cartesian coordinate plane do not in any way vitiate the results obtained, since the shear stress free boundaries were never imposed as a condition for the formulation of the governing equations unlike the classical Boussinesq problem where the displacement functions were derived to satisfy apriori the shear stress free boundary conditions on the xy plane.

6. CONCLUSIONS

The following conclusions are made from the study:

  • (i) The Hankel transformation method has been successfully used to solve the Westergaard problem of finding vertical stresses in the elastic half-space due to a point load at the origin, and the solution used as Green functions to obtain solutions for uniformly distributed loads on line and on a rectangular area.

  • (ii) The Hankel transformation of the governing Cauchy - Navier displacement equation of equilibrium formulated in cylindrical polar coordinates simplified the boundary value problem (BVP) to a second order ordinary differential equation which could be solved using methods for solving ODEs, to obtain bounded solutions for the displacement.

  • (iii) The displacement function obtained was found to be defined and single valued at all points in the elastic half-space except at the origin which is the point of application of the concentrated load at which the displacement function becomes singular, (undefined or indeterminate).

  • (iv) The normal stresses at any point on the xy coordinate plane are singular, (undefined or indeterminate).

  • (v) The shear stress τxy on the xy plane, as well as on the entire elastic half-space is zero. This is due to the foundational assumption of longitudinal inextensibility which made u = v = 0, and thus the shear stress γxy would naturally vanish from the strain - displacement relations for all points in the elastic half-space.

  • (vi) The shear stresses τzy and τzx do not vanish at points on the xy plane for which z = 0, but are singular (undefined or indeterminate) at the origin.

REFERENCES

  • Abeyartne R. (2012). Continuum Mechanics, Volume 11 of Lecture notes on the mechanics of Elastic Solids. Cambridge. http://web.mit.edu/abeyartne/lecture notes/html 11th may 2012. Retrieved 8th April 2016.
    » http://web.mit.edu/abeyartne/lecture notes/html
  • Andrews L.C. and Shivamoggi B.K. (1999). Integral Transforms for Engineers. SPIE Optical Engineering Press. A publication of SPIE - The International Society for Optical Engineering, Bellingham Washington USA ISBN: 0-8194-3232-6.
  • Anyaegbunam A.J., Osadebe N.N., Eze-Uzoamaka O.J. (2011). Non-Existence of Solution for Horizontally Rigid Half-Space ASCE, Journal of Geotechnical and Geoenvironmental Engineering 137(4) April, pp 431 - 434. http://ascelibrary.org/gto/ resource/1/jgefk/v.137/i4.doi.10.1061/(ASCE)GT 1943- 5606.0000444
    » http://ascelibrary.org/gto/ resource/1/jgefk/v.137/i4.doi.10.1061/(ASCE)GT 1943- 5606.0000444
  • Apostol B.F. (2017). Elastic Displacement in Half Space under the Action of a Tension force, General Solution for the Half Space with Point Forces. Journal of Elasticity Vol 126 pp 231 - 244. DOI:10.1007/S10659-016-9592-3.
    » https://doi.org/10.1007/S10659-016-9592-3
  • Barbar J.R. (2010). Elasticity Third Revised Edition. Springer Science and Bussiness Media. Dordrecht The Netherlands. DOI: https/doi.org/10.007/978-90-481-38098-4
    » https/doi.org/10.007/978-90-481-38098-4
  • Barden L. (1963). Stresses and displacements in a cross-anisotropic soil. Geotechnique 13 pp. 198 - 210.
  • Bowles J.E. (1997). Foundation Analysis and Design, International Edition. McGraw Hill International Book Company, Tokyo (1997) ISBN 0-07-912247-7.
  • Chan K.T. (2013). Analytical methods in Geomechanics CRC Press Taylor and Francis Group New York.
  • Davis R.O. and Salvadurai A.P.S. (1996). Elasticity and Goemechanics. Cambridge University Press.
  • Duffy D.G. (2004). Transform Methods for Solving Partial Differential Equations Second Edition. Chapman and Hall/CRC ISBN 9781584884514.
  • Everitt W.N. and Kalf H. (2007). The Bessel Differential Equation and the Hankel transform. Journal of Computational and Applied Mathematics, Elsevier Volume 208, Issue 1, November 2007, pp 3 - 19. https://doi.org/10.1016/j.cam.2006.10.029.
  • Fadum Ralph E. (1948). Influence values for estimating stresses in elastic foundations. Proceedings Second International Conference on Soil Mechanics and Foundation Engineering Vol 3 pp 77 - 84. https://www.issmge.org/uploads/ publications/1/43/ 194803-0020.pdf
    » https://www.issmge.org/uploads/ publications/1/43/ 194803-0020.pdf
  • Green A.E. and Zerna W. (1954). Theoretical Elasticity, Oxford University Press, London.
  • Hazel A. MATH 35021 Elasticity www.maths.manchester.ac.lc/_ahazel/MATHS Nov 30, 2015.
    » www.maths.manchester.ac.lc/_ahazel/MATHS
  • Ike C.C. (2006). Principles of Soil Mechanics. De-Adroit Innovation, Enugu.
  • Ike C.C. (2017). First principles derivation of a stress function for axially symmetric elasticity problems, and application to Boussinesq problem. Nigerian Journal of Technology (NIJOTECH) Vol 36 No 3 July 2017 pp 767 - 772. http://dx-doi/10.4314/nijt/v.36i3.15 (www.nijotech.com).
    » http://dx-doi/10.4314/nijt/v.36i3.15
  • Ike C.C. (2018a). Hankel transform method for solving axisymmetric elasticity problems of circular foundation on semi-infinite soils. International Journal of Engineering and Technology (IJET) Vol 10. No 2 April - May 2018 pp 549 - 564. DOI: 10.21817/ijet/2018/v10 i2/181002111.
    » https://doi.org/10.21817/ijet/2018/v10 i2/181002111
  • Ike C.C. (2018b). General Solutions for Axisymmetric Elasticity Problems of Elastic Half Space using Hankel Transform Method. International Journal of Engineering and Technology (IJET) Vol 10. No 2 April - May 2018 pp 565 - 580. DOI: 10.21817/ijet/2018/v10i2/181002112.
    » https://doi.org/10.21817/ijet/2018/v10i2/181002112
  • Ike C.C., Onah H.N. and Nwoji C.U. (2017a). Bessel functions for axisymmetric elasticity problems of the elastic half space soil, a potential function method. Nigerian Journal of Technology (NIJOTECH) Vol 36 No 3 July 2017 pp. 773 - 781. http://dx.doi.org/10.4314/nijt/v.36i3.16
    » http://dx.doi.org/10.4314/nijt/v.36i3.16
  • Ike C.C., Mama B.O., Onah H.N. and Nwoji C.U. (2017b). Trefftz Harmonic function method for solving Boussinesq problem. Electronic Journal of Geotechnical Engineering (EJGE) 2017 (22.12) pp 4589 - 4601. Available at ejge.com.
  • Kachanov M.L., Shafiro B. and Tsukrov I. (2003). Handbook of Elasticity Solutions. Springer Science and Bussiness Media Kluwer Academic Publishers Dordrecht The Netherlands.
  • Liao J.J. and Wang C.D. (1998). Elastic solutions for a transversely isotropic half-space subjected to a point load. International Journal for Numerical and Analytical Methods in Geomechanics (Int. J. Numer. Anal. Meth. Geomech) 22 pp. 425 - 447.
  • Lurie S.A. and Vasilev V.V. (1995). The Biharmonic Problem in the Theory of Elasticity Gordon and Breach Publishers United States.
  • Malacka Zuzana S.K. (2018). Hankel Transform and Free Vibration of a Large Circular Membrane. 17th Conference on Applied Mathematics APLIMAT 2018 Faculty of Mechanical Engineering, Slovak University of Technology in Bratislava pp 717 - 723. evim.stuba.sk/APLIMAT2018/proceedings/papers/0717_Malacka.pdf.
  • Misra B. and Sen B.R. (1975). Stresses and displacements in granular materials due to surface load. Int. J. Engng Sci. 13 pp. 743 - 761.
  • Negero N.T. (2016). Zero Order Hankel Transform Method for Partial Differential Equations. International Journal of Modern Sciences and Engineering Technology (IJMSET) Volume 3 Issue 10, 2016 pp 24 - 36. ISSN 2349 - 3755. Available at www.ijmset.com
  • Nwoji C.U., Onah H.N., Mama B.O., Ike C.C. (2017a). Solution of elastic half space problem using Boussinesq displacement potential functions. Asian Journal of Applied Sciences (AJAS) Vol 5, No 5 October 2017 pp. 1100 - 1106.
  • Nwoji C.U., Onah H.N., Mama B.O., Ike C.C. (2017b). Solution of the Boussinesq problem using Green and Zerna displacement potential function method. The Electronic Journal of Geotechnical Engineering (EJGE) 2017 Volume 22 Bundle 11 (22.11) pp. 4305 - 4314. Available at ejge.com.
  • Ojedokun O.Y. and Olutoge F.A. (2012). Application of Boussinesq’s and Westergaard’s formulae in analysing foundation stress distribution for a failed telecommunication mast. African Journal of Mathematics and Computer Science Research Vol 5(4) pp 71 - 77 15 February 2012. Available on line at http://www.academicjournals.org/AJMCSR DOI:10.5897/AJMCSR 11.155 ISSN:2006-9731 Academic Journals.
    » http://www.academicjournals.org/AJMCSR DOI:10.5897/AJMCSR
  • Padio-Guidugli P. and Favata A. (2014). Elasticity for Geotechnicians. A modern exposition of Kelvin, Boussinesq, Flammant, Cerrutti, Melan and Mindlin problems. Solid mechanics and its applications. Springer.
  • Palaniappan D. (2011). A general solution of equations of equilibrium in linear elasticity. Applied Mathematical Modelling 35 (2011) pp 5494 - 5499 Elsevier.
  • Piessens R. (2000). “The Hankel Transform” In “The Transforms and Applications Handbook: Second Edition” Ed. Alexander D. Poularikas Boca Raton: CRC Press LLC.
  • Rocscience (2018). Westergaard stress solution method. https://www.rocscience.com/ documents/../Westergaard_Stress_Solution_Method pdf Retrieved 02/09/2018.
    » https://www.rocscience.com/ documents/../Westergaard_Stress_Solution_Method pdf
  • Sadd M.H. (2014). Elasticity Theory, Application and Numerics, Third Edition. University of Rhode Island. Elsevier Academic Press Amsterdam https://personnel egr.uri/edu/teggart/ courses/ce571/chapter 7. ppx
    » https://personnel egr.uri/edu/teggart/ courses/ce571/chapter 7. ppx
  • Sitharam T.G. and Govinda Reju L. (2017). Applied Elasticity for Engineers Module: Elastic Solutions with Applications in Geomechanics 14.139.172.204/nptel/1/CSE/web/ 105108070/module 8/lecture 17.pdf.
  • Sokolnikoff I.S. (1956). Mathematical Theory of Elasticity Second Edition. Data McGraw Hill Publishing Company Ltd Bombay New Delhi.
  • Tarn J.Q. and Wang Y.M. (1987). A fundamental solution for a transversely isotropic elastic space. J. Chin. Inst. Eng Volume 10 Issue 1. pp 13 - 21.
  • Teodorescu P.P. (2013). Treatise on Classical Elasticity Theory and Related Problems. Mathematical and Analytical Techniques with Applications to Engineering. Springer Dordrecht. DOI:10.1007/978-94-007-2616-1.
    » https://doi.org/10.1007/978-94-007-2616-1
  • Timoshenko S.P. and Goodier J.N. (1970). Theory of Elasticity Third Edition, McGraw Hill New York.
  • Tuteja R., Jaloree S., and Goyal A. (2014). Application of Hankel transform of I-function of one variable for solving axisymmetric Dirichlet potential problem. Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 14, Issue 4, Version 1.0 pg 11 - 16. Online ISSN: 2249-4626.
  • Voegtle I.C. (2017). The Bessel Function, the Hankel Transform and an Application to Differential Equations. Electronic Thesis and Dissertations Jack N. Averity College of Graduate Studies, Georgia Southern University Summer 2017. https://digitalcommons.georgiasouthern.edu/etd
    » https://digitalcommons.georgiasouthern.edu/etd
  • Westergaard H.M. (1938). A problem of elasticity suggested by a problem in soil mechanics: soft soil reinforced by numerous strong horizontal sheets Contributions to the mechanics of solids. Stephen Timoshenko 60th Birthday Anniversary Volume Macmillan New York.
  • Westergaard H.M. (1964). Theory of Elasticity and Plasticity. Dover New York.
  • Yokoyama Takehito (2014). The Hankel transform. www.stat.phys.titech.ac.jp/~yokayama /memo4.pdf Accessed 1/7/2018.
    » www.stat.phys.titech.ac.jp/~yokayama /memo4.pdf
  • Available online: November 21, 2018.

Publication Dates

  • Publication in this collection
    2019

History

  • Received
    03 Oct 2018
  • Reviewed
    20 Nov 2018
  • Accepted
    21 Nov 2018
Individual owner www.lajss.org - São Paulo - SP - Brazil
E-mail: lajsssecretary@gmsie.usp.br