The behaviour of an elastic non-homogeneous half-space. Part I-line and point loads

In the first part of this paper solutions are developed for the response of a non-homogeneous half-space subjected to either a surface point load or a surface line load. The non-homogeneity considered is a variation in Young's modulus (E) with depth (z) which takes the form E=m_EZ^α where m_E is a constant and α is referred to as the non-homogeneity parameter. The variation of these solutions as the non-homogeneity parameter α varies between the limits of zero (homogeneous soil) to unity (Gibson soil) gives some fresh insight into both these limiting cases.     

SURFACE DISPLACEMENTS OF A NON-HOMOGENEOUS ELASTIC HALF-SPACE SUBJECTED TO UNIFORM SURFACE TRACTIONS. PART I: LOADING ON ARBITRARILY SHAPED AREAS

A numerical technique is presented for the analysis of surface displacements of a non-homogeneous elastic half-space subjected to vertical and/or horizontal surface loads uniformly distributed over an arbitrarily shaped area. The non-homogeneity considered is a particular form of power variation of Young's modulus with depth. Since the exponent which determines the degree of non-homogeneity may vary from zero to unity, both the homogeneous half-space and the Gibson soil may be included as limiting cases in a single numerical scheme. In order to account for the arbitrary shape of the loading, the boundary of the loaded area is linearized piecemeal. This enables the modeling of any load pattern according to the desired degree of accuracy. Special attention is focused on the integration scheme, since the singularity associated with the Green's function becomes progressively more pronounced the greater the non-homogeneity parameter gets. The performance of the numerical procedure is studied using analytical solutions for rectangular shaped areas. Further comparisons with well-known solutions based on integral transform techniques for a uniformly distributed load acting on a circular area of the non-homogeneous soil mass show excellent agreement as well. © 1997 John Wiley & Sons, Ltd.     

The axial loading of a rigid circular anchor plate embedded in an elastic half-space

The present paper examines the axisymmetric problem related to the loading of a rigid circular anchor plate which is embedded in bonded contact with an isotropic elastic half-space. A Hankel transform development of the governing equations is used to reduce the associated mixed boundary value problem to a set of coupled Fredholm integral equations of the second kind. These equatons are solved in a numerical fashion to generate results of engineering interest. In particular, the results indicate the influence of the depth of embedment on the axial stiffness of the rigid anchor plate.     

竖向圆形荷载作用下弹性地基半空间问题的解析解

竖向圆形荷载作用下弹性半空间问题的位移和应力解是桩基分析的基础。利用Hankel积分变换,首先导出了弹性地基半空间位移与应力的积分形式的通解。通过适当地引入边界条件和界面位移和应力的连续条件,求得了内部作用竖向圆形荷载时弹性地基半空间位移与应力的积分形式解。在此基础上,给出了不同深度处荷载作用投影范围内竖向位移和竖向正应力的平均值。数值结果验证了解析解的正确性。     

Finite element modelling of thick plates on two‐parameter elastic foundation

This paper is intended to give some information about how to build a model necessary for bending analysis of rectangular and circular plates resting on a two‐parameter elastic foundation, subjected to combined loading and permitting various types of boundary conditions. The formulation of the problem takes into account the shear deformation of the plate and the surrounding interaction effect outside the plate. The numerical model based on an 18‐node zero‐thickness isoparametric interface element interacting with a thick Reissner–Mindlin plate element with three degrees of freedom at each of the nine nodes, which enforce C⁰ continuity requirements for the displacements and rotations of the midsurface, is proposed. Stiffness matrices of a special interface element are superimposed on the global stiffness matrix to represent the stiffening elastic foundation under and beyond the plate. Some numerical examples are given to illustrate the advantages of the method presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Vertical and horizontal deformation of an inhomogeneous elastic half-space

Vertical and horizontal deformations of surface footings have been studied for an inhomogeneous elastic half-space in which the shear modulus increases with an arbitrary power of depth, n, and Poisson's ratio is constant. A general solution for displacements has been obtained first for point loads applied in vertical and horizontal directions. These are then used in obtaining closed-form solutions for displacements of uniformly loaded circular and rectangular footings. Finally, a numerical method is described that can be used to analyse a rigid footing of an arbitrary shape, and results for rigid rectangular footings are given.     

Vertical deformation of rigid foundations of arbitrary shape on layered soil media

A numerical procedure is described for the analysis of vertical deformation of smooth, rigid foundations of arbitrary shape on homogeneous and layered soil media. The contact area at the interface of the foundation and soil medium is approximated by square subdivisions. The response of the system is then obtained from the superposition of the influence of the individual subdivisions. The flexibility influence coefficients are based on equivalent smooth, rigid circular areas with the same contact area as the square subdivisions. For foundations on a homogeneous, isotropic elastic half-space, the flexibility coefficients are given analytically by the integrated forms of the Boussinesq's solution. For a layered soil medium, the flexibility coefficients are determined from an axisymmetric finite element analysis which is essentially two dimensional. Thus, there is no necessity for a full three-dimensional finite element analysis. Comparison with solutions obtained using the integral transform technique for smooth, rigid rectangular foundations on a homogeneous, isotropic elastic half-space shows good agreement. Parametric solutions are presented for the response of rectangular foundations on some ‘typical’ soil profiles. The use of a simplified method to estimate the settlement of rectangular foundations on a layered soil medium by superposing solutions for homogeneous, elastic strata is discussed.     

饱和地基表面基础在Rayleigh波作用下的摇摆振动分析

考虑地基与基础的相互作用,运用Biot波动方程理论,研究了饱和地基表面有质量的刚性圆形基础在入射Rayleigh波作用下的摇摆振动问题。将饱和地基中的总波场分解为自由波场、刚体辐射波场及刚体散射波场,从而考虑基础对Rayleigh波的散射。采用Hankel变换求解土体控制方程,结合基础摇摆振动的混合边值条件以及Rayleigh波的动力作用获得基础摇摆振动振幅的表达式。计算结果分析了泊松比、基础质量、透水条件、土体渗透系数等对基础摇摆振动的影响。研究结果表明,与弹性介质相比,饱和地基中孔隙水的存在减弱了基础的摇摆振动,其作用效果与土体渗透系数及地基表面排水条件密切相关。     

Empirical Relationships Between the Elastic Settlement of Rigid Rectangular Foundations and the Settlement of the Respective Flexible Foundations

The problem of settlement of shallow foundations is among the most important ones in classical soil mechanics. And while for the settlement of flexible foundations elastic solutions are widely used, for rigid rectangular foundations where the actual contact pressure distribution is still unknown, the problem is approximated either analytically assuming a contact pressure distribution or semi‐empirically combining the theory of elasticity with experimental and/or numerical results. A third and often attractive choice is the use of simple empirical relationships or relevant tabulated values relating the elastic settlement of rigid foundations (ρ_R) with the settlement of the respective flexible foundations (e.g. at the center, ρ_Ce). Reviewing the relathionships of this third approach, the author revealed serious lack of consesous between the various sources; for example, according to the literature, ρ_R ranges between 68 and 125% of ρ_Ce, the time when it is well-known that ρ_R?<?ρ_Ce. In this paper, comparison of the settlement of 210 rigid foundation cases derived from 3D elastic finite element analysis, with the settlement of the respective flexible foundations derived from the theory of elasticity, led to simple empirical relationships between ρ_R and ρ_Ce as well as between ρ_R and ρ_Av (ρ_Av?=?average settlement of the flexible foundation) with coefficient of determination (R²) almost unity. The analysis showed that these relationships are largely independent of the aspect ratio of foundations and the thickness and Poisson’s ratio (ν) of the compressible medium, although separate relationships are given for ν?=?0.5, slightly increasing R². Finally, a correction factor for foundation rigidity is given exploting the known linear relationship that exists between the relative stiffness factor of foundations and settlement.

     

含圆柱形弹性夹杂地基沉降量的解析解

研究各向同性地基中含一个圆柱形弹性夹杂的地基沉陷问题。考虑全场的应力、变形特征,运用复变函数的罗朗级数展开技术,构造合适的应力函数,结合夹杂的边界条件和无穷远的性态,获得了全场复势的解析表达。数值结果显示：由于两相材料剪切模量之比的不同,在某个范围内沉陷量大于基底中点沉陷;在这个范围之外,则小于基底中点沉陷。当夹杂埋深为夹杂半径的5倍以上时,就可以运用所得结果,而不致引起太大误差;含圆形隧洞和刚性夹杂的问题可作为特殊情况得到。     

The Reissner–Sagoci problem for a transversely isotropic half‐space

A transversely isotropic linear elastic half‐space, z?0, with the isotropy axis parallel to the z‐axis is considered. The purpose of the paper is to determine displacements and stresses fields in the interior of the half‐space when a rigid circular disk of radius a completely bonded to the surface of the half‐space is rotated through a constant angle θ₀. The region of the surface lying out with the circle r?a, is free from stresses. This problem is a type of Reissner–Sagoci mixed boundary value problems. Using cylindrical co‐ordinate system and applying Hankel integral transform in the radial direction, the problem may be changed to a system of dual integral equations. The solution of the dual integral equations is obtained by an approach analogous to Sneddon's (J. Appl. Phys. 1947; 18 :130–132), so that the circumferential displacement and stress fields inside the medium are obtained analytically. The same problem has already been approached by Hanson and Puja (J. Appl. Mech. 1997; 64 :692–694) by the use of integrating the point force potential functions. It is analytically proved that the present solution, although of a quite different form, is equivalent to that given by Hanson and Puja. To illustrate the solution, a few plots are provided. The displacements and the stresses in a soil deposit due to a rotationally symmetric force or boundary displacement may be obtained using the results of this paper. Copyright © 2006 John Wiley & Sons, Ltd.     

PROPAGATION OF RAYLEIGH WAVES IN A HETEROGENEOUS INCOMPRESSIBLE SUBSTRATUM OVER A HOMOGENEOUS INCOMPRESSIBLE HALF-SPACE

The paper studies the propagation of Rayleigh waves in an incompressible layer with general variation of rigidity as μ=μ₀(1+bz)^m resting over a homogeneous incompressible elastic half-space. Instead of using the Whittaker function, the expansion formula proposed by Newlands has been used for better results at shallow depths. As a particular case for m=1, the results have been shown to coincide with those obtained by Newlands. The velocities have been computed for different values of m and the results are presented in graphs.     

The eccentric loading of a rigid circular foundation embedded in an isotropic elastic medium

This paper examines the problem of the eccentric loading of a rigid circular disc-shaped foundation embedded in bonded contact with an istropic elastic medium of infinite extent. The solution of this problem is achieved by using a Hankel integral transform technique which reduces the problem to two sets of dual integral equations. These systems of dual integral equations represent the equations which govern the axisymmetric and asymmetric states of deformation induced by the loaded foundation. Explicit results are derived both for the displacement and rotation of the circular foundation and for the contact stress at the interface.     

A note on the consolidation settlement of a rigid circular foundation on a poroelastic halfspace

This paper uses Biot's poroelasticity approach to examine the consolidation behaviour of a rigid foundation with a frictionless base in contact with a poroelastic halfspace. The mathematical development of the mixed boundary value problem involves a set of dual integral equations in the Laplace transform domain which cannot be conveniently solved by employing conventional procedures. In this paper, a numerical solution is developed using a scheme where the contact normal stress is approximated by a discretized equivalent. The influence of limiting drainage boundary conditions at the entire surface of the halfspace on the degree of consolidation of the rigid circular foundation is investigated. The results obtained in this study are compared with the corresponding results given in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Fundamental results concerning the settlement of a rigid foundation on an elastic medium due to an adjacent surface load

This paper examines the interaction between a rigid circular foundation resting in smooth contact with an isotropic elastic halfspace and a concentrated surface load which acts at a finite distance from the foundation. Owing to the action of the external load the rigid foundation experiences an extra settlement and a tilt. The expressions for the extra settlement and the tilt are evaluated in exact closed form. It is also shown that these deformations due to the external load satisfy Betti's reciprocal theorem. The auxiliary solution required for the application of the reciprocal theorem is derived from the analysis of the problem of a rigid circular foundation resting in smooth contact with an elastic medium and subjected to an eccentric load. The, results developed for the interaction between the rigid circular foundation and the external concentrated load are utilized to generate, among others, solutions for the settlement and tilt induced at a rigid foundation due to its interaction with uniformly or non-uniformly distributed loads with circular and square plan shapes.     

Semi‐analytical far field model for three‐dimensional finite‐element analysis

A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd.     

A theoretical estimate of the tilt of a surface foundation due to an inclined anchor load applied at the interior of the soil medium

This paper examines the problem of the interaction between a loaded rigid circular foundation located at the surface of an isotropic elastic halfspace and an inclined concentrated anchor load which is located at a finite depth along the axis of symmetry. Such inclined loads can be induced by, for example, anchor regions supporting earth retaining structures. The loaded rigid circular foundation resting in smooth contact with the elastic soil mass experiences a displacement and a tilt due to the action of the inclined anchor load. The magnitude of the rotational settlement is evaluated in exact closed form.     

The behaviour of an elastic non-homogeneous half-space. Part II–circular and strip footings

Solutions developed in the first part of this paper (i.e. describing the response of a non-homogeneous half-space subjected to surface point and line loads) are used in this part to obtain solutions for a variety of surface loadings. Consideration is given to a distributed load acting over a circular area or strip and a rigid disk or strip subjected to applied normal load and moment. It is established that the profiles of surface settlement due to uniformly distributed loads acting over a strip or circular area are strongly dependent on the degree of non-homogeneity. This dependency is reduced when the footing is rigid. When α = 1 the moduli variation is identical to the Gibson soil and the equivalence with the Winkler soil model is established.     

A simple model for inelastic footing response to transient loading

This paper presents a mechanical analogue which models the response of a rigid circular footing on an ideal elastoplastic half-space to transient loads. In the rational analysis of pile-driving dynamics, the response of soil at the base of a pile is often approximated by a footing on a semi-infinite half-space. Most existing base models employ the well-known Lysmer analogue to model the elastic response of the soil at the pile base, and account for the inelastic soil behaviour through the inclusion of a plastic slider with a slip load equal to the ultimate failure load of the footing. The improved model provides a force response which is significantly closer to the ideal response than existing models. The paper commences with a review of analytical solutions for the dynamic response of a rigid circular footing on an elastic half-space. Existing mechanical analogs for the system are reviewed, and an automatic matching process proposed which improves the accuracy of the analogs under transient loading. The inelastic response is then studied using the finite element method, and the mechanical analogs are modified to allow representation of the observed inelastic behaviour. Examples are presented illustrating close agreement between the proposed models and finite element analyses for a range of Poisson's ratio. The improved models have direct application for one-dimensional models of pile driving, particularly in the back-analysis of data from dynamic testing of piles. They are also applicable to studies of dynamic compaction.     

3D Elastic Solutions for Laterally Loaded Discs: Generalised Brazilian and Point Load Tests

This paper investigates the application of a double Fourier series technique to the construction of an elastic stress field in a cylindrical bar subject to lateral boundary loads. The lateral loads, including the constant load boundary conditions, are represented by two Fourier series: one on the perimeter of the circular section (r ₀, θ) and the other on the longitudinal curved surface parallel to the bar axis (z). The technique invokes acceptable potential functions of the Papkovich–Neuber displacement field, satisfying the governing partial differential equations, to assign appropriate odd and even trigonometric Fourier terms in cylindrical coordinates (r, θ, z). The generic solution decomposes the problem of interest to a state of stress caused by two independent boundary conditions along the z axis and θ-polar angle, both superimposed on a solution for which these potentials are the product of the trigonometric terms of the independent variables (θ, z). Constants appearing in the resultant second-order partial differential equations are determined from the generally mixed (tractions and/or displacements) boundary conditions. While the solutions are satisfied exactly at the ends of an infinite bar, they are satisfied weakly on average, in the light of Saint Venant’s approximation at the two ends of a finite bar. The application of the proposed analysis is verified against available elastic solutions for axisymmetric and non-axisymmetric engineering problems such as the indirect Brazilian Tensile Strength and Point Load Strength tests.