Green functions for an incompressible linearly nonhomogeneous half-space

Summary Time-harmonic vibrations of an incompressible half-space having shear modulus linearly increasing with depth are studied. The half-space is subjected to a surface load which has vertical or hovizontal direction. The general solution of the time-harmonic, in the vertical direction nonhomogeneous problem is constructed for arbitrary angular distribution in the horizontal plane. Numerical results concerning surface displacements due to a point force are given for the case of nonzero shear modulus at the surface. These results show that nonhomogeneity can considerably increase amplitudes at large distances from the applied force.     

Integration of Singularities in FE/BE Analyses of Soil-foundation Interaction with Non-homogeneous Elastic Soils

This paper addresses the derivation of the boundary integral equations for a non-homogeneous elastic half-space subjected to constant surface tractions on an arbitrarily shaped area on the basis of the respective Green's functions. The type of non-homogeneity considered is a power law variation of Young's modulus with depth below the surface of the half-space. Two different methods, a contour integral and an integration-free approach are presented, applicable for arbitrarily and rectangular shaped boundary elements, respectively. In the first one the divergence theorem is used in order to reduce the integration of a two-dimensional surface element to an integration over the element's confining boundary only. In the second approach no integration at all is needed since the solution is found simply by evaluating functions to be determined at the boundaries of the loaded rectangle.     

Contact vibration analysis of the functionally graded material coated half-space under a rigid spherical punch

This paper studies the contact vibration problem of an elastic half-space coated with functionally graded materials (FGMs) subject to a rigid spherical punch. A static force superimposing a dynamic time-harmonic force acts on the rigid spherical punch. Firstly, we give the static contact problem of FGMs by a least-square fitting approach. Next, the dynamic contact pressure is solved by employing the perturbation method. Lastly, the dynamic contact stiffness with different dynamic contact displacement conditions is derived for the FGM coated half-space. The effects of the gradient index, coating thickness, internal friction, and punch radius on the dynamic contact stiffness factor are discussed in detail.     

The axisymmetric lamb’s problem for a finite prestrained half-space covered with a finite prestretched layer

The piecewise-homogeneous body model and the three-dimensional linearized theory of elastic waves in prestressed bodies are used to solve the axisymmetric time-harmonic Lamb’s problem for a finite prestrained half-space covered with a finite prestretched layer. It is assumed that the half-space and layer are incompressible and their deformation is described by the Treloar potential. The normal stress at the interface is calculated Published in Prikladnaya Mekhanika, Vol. 43, No. 3, pp. 132–143, March 2007.     

The influence of the third order elastic constants on the dynamical interface stress field in a half-space covered with a pre-stretched layer

The influence of the third order elastic constants on the dynamical (time-harmonic) axisymmetric interface stress field in the system which comprises the half-space and the pre-stretched covering layer is investigated within the framework of the piecewise homogeneous body model by employing the three-dimensional linearized theory of elastic waves in initially stressed bodies. The elasticity relations for the layer and half-space materials are given through the Murnaghan potential. It is assumed that the force acting on the free face plane of the covering layer is a time-harmonic point-located normal force. The influences due to both of the quantities of the pre-stretching of the layer and the third order elastic constants of the layer material on the interface normal stress are analysed. The numerical results are presented for concrete selected materials such as steel, aluminium and acrylic plastic.     

Surface displacements due to elastic wave scattering by buried planar and non-planar cracks

The scattering of time-harmonic plane longitudinal, shear, and Rayleigh waves by a crack in two dimensions embedded in a semi-infinite homogeneous isotropic elastic half-space has been studied in this paper. Two problems have been considered: a straight crack and a Y-shaped crack. A hybrid numerical technique combining a multipolar representation of the scattered field in the half-space with the finite element method has been used to obtain the far-field displacements as well as the stress-intensity factors for the crack tips. Results for vertical displacement on the free surface of the half-space are presented in this paper.     

Surface Green function for a soft elastic half-space: Influence of surface stress

Surface Green function for incompressible, elastically isotropic half-space coupled with surface stress is derived by using double Fourier transform technique. The result indicates that the surface displacement induced by a force tangential to the surface is the same as the usual solution for elastic half-spaces where the effect of surface stress is ignored. However, the displacement caused by a force normal to the surface involves an additional parameter, i.e. the ratio of specific surface stress to shear modulus. The parameter has the dimension of length, and may provide a means to introduce an intrinsic length scale for some related problems regarding the surface of an elastic half-space. This is extremely true for soft elastic media with very low shear modulus, because in that situation the magnitude of the parameter is relatively large. As an illustrative example, the proposed Green function is adopted to analyze the interaction between two molecules with circular section adsorbed on the surface of a soft elastic half-space. It is shown that surface stress remarkably affects the pair interaction potential when the distance between the molecules is not larger than several times of the intrinsic length scale.     

Anti-plane shear Green’s functions for an isotropic elastic half-space with a material surface

This paper presents analytical Green’s function solutions for an isotropic elastic half-space subject to anti-plane shear deformation. The boundary of the half-space is modeled as a material surface, for which the Gurtin–Murdoch theory for surface elasticity is employed. By using Fourier cosine transform, analytical solutions for a point force applied both in the interior or on the boundary of the half-space are derived in terms of two particular integrals. Through simple numerical examples, it is shown that the surface elasticity has an important influence on the elastic field in the half-space. The present Green’s functions can be used in boundary element method analysis of more complicated problems.     

Dynamic interaction of a pile with a transversely isotropic elastic half-space under transverse excitations

This investigation is concerned with a mathematical analysis of an elastic circular cylindrical pile embedded in a transversely isotropic half-space under lateral dynamic excitations. A combination of time-harmonic horizontal shear force and moment are applied at the top end of the pile. The boundary value problem is formulated by decomposing the pile-medium system into a fictitious pile and an extended transversely isotropic half-space. A Fredholm integral equation of the second kind governs the interaction problem, whose solution is then computed numerically. Selected results for dynamic compliance bending moment, displacement and slope profiles are presented for different transversely isotropic half-spaces to portray the influence of degree of anisotropy of the medium on various aspects of the solution.     

On the stress field singularity in an incompressible half-space weakened by two near-surface wedge-like cracks

We consider the equilibrium problem for an elastic incompressible half-space weakened by two near-surface wedge-like cracks, whose lie in the same plane perpendicular to the half-space surface and have a common vertex. We use the Papkovich-Neuber representation to reduce the problem to finding two harmonic functions satisfying the mixed boundary conditions. These functions are constructed in spherical coordinates by using a Mehler-Fock type integral representation in Legendre functions. The analytic solution thus obtained permits finding the character of the stress distribution near the common tip of the cracks.     

波数-频率域内地基土表面位移Green函数的理论分析

建立了柱面坐标系下分层弹性半空间地基土模型。利用钟阳刚度矩阵法和Haskell-Thomson传递矩阵法推导出所有分层土体之间的振动传递关系;根据Helmholtz定理将土体的位移向量分解成势函数的形式,推导出弹性半空间表面应力与位移之间的关系;再将分层土体和半空间地基土通过位移与应力之间的关系进行耦合,得到分层弹性半空间地基土模型表面位移与应力之间的关系。结合单位脉冲荷载作用下地基土表面的边界条件,推导出波数-频率域内地基土表面位移Green函数的解析解,用Matlab程序语言对理论进行实现并通过算例对地基土表面位移Green函数的特征进行了分析和总结。     

Green’s functions of a surface-stiffened transversely isotropic half-space

Green’s functions of a transversely isotropic half-space overlaid by a thin coating layer are analytically obtained. The surface coating is modeled by a Kirchhoff thin plate perfectly bonded to the half-space. With the aid of superposition technique and employing appropriate displacement potential functions, the Green’s functions are expressed in two parts; (i) a closed-form part corresponding to the transversely isotropic half-space with surface kinematic constraints, and (ii) a numerically evaluated part reflecting the interaction between the half-space and the plate in the form of semi-infinite integrals. Some limiting cases of the problem such as surface-stiffened isotropic half-space, Boussinesq and Cerruti loadings, and extremely flexible and rigid plates are also studied. For the classical Cerruti problem in transversely isotropic materials, the effects of incompressibility are highlighted. Numerical results are provided to show the effects of material anisotropy, relative stiffness factor, and load buried depth. The obtained Green’s functions play a key role in treating further mixed-boundary-value problems in surface stiffened transversely isotropic half-spaces.     

Elastic Displacement in a Half-Space Under the Action of a Tensor Force. General Solution for the Half-Space with Point Forces

The elastic displacement in an isotropic elastic half-space with free surface is calculated for a point tensor force which may arise from the seismic moment of seismic sources concentrated at an inner point of the half-space. The starting point of the calculation is the decomposition of the displacement by means of the Helmholtz potentials and a simplified version of the Grodskii-Neuber-Papkovitch procedure. The calculations are carried out by using generalized Poisson equations and in-plane Fourier transforms, which are convenient for treating boundary conditions. As a general result we compute the displacement in the isotropic elastic half-space with free surface caused by point forces with arbitrary structure and orientation, localized either beneath the surface (generalized Mindlin problem) or on the surface (generalized Boussinesq-Cerruti problems). The inverse Fourier transforms are carried out by means of Sommerfeld-type integrals. For forces buried in the half-space explicit results are given for the surface displacement, which may exhibit finite values at the origin, or at distances on the surface of the order of the depth of the source. The problem presented here may be viewed as an addition to the well-known static problems of elastic equilibrium of a half-space under the action of concentrated loads. The application of the method to similar problems and another approach to the starting point of the general solution are discussed.     

Finite-amplitude Love waves in a pre-stressed compressible elastic half-space with a double surface layer

The dispersive behavior of finite-amplitude time-harmonic Love waves propagating in a pre-stressed compressible elastic half-space overlaid with two compressible elastic surface layers of finite thickness is investigated. The half-space and layers are made of different pre-stressed compressible neo-Hookean materials. The dispersion relation which relates wave speed and wavenumber is obtained in explicit form. Results for the energy density and energy flux of the waves are also presented. The special case where the interfaces between the layers and the half-space are principal planes of the left Cauchy–Green deformation tensor is also investigated. Numerical results are presented showing the variation of the Love wave speed with the pre-stress and the propagation angle.     

Dynamics of a prestressed incompressible layered half-space under moving load

The paper addresses a plane problem: a concentrated force acts on a plate resting on an elastic half-space with homogeneous prestrain. The equations of motion of the plate incorporate shear and rotary inertia. The half-space is assumed to be incompressible and isotropic in the natural state. The elastic potential is given in general form and is only specified for numerical purposes. The dependence of the critical velocity of the load and the stress-strain state on the prestresses is analyzed for different ratios between the stiffnesses of the layer and half-space and different contact conditions. The calculations are carried out for a half-space with Bartenev-Khazanovich potential __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 3, pp. 36–54, March 2008.     

Dynamic Response of a Multilayered Poroelastic Half-Space to Harmonic Surface Tractions

The propagator matrix method is developed to study the dynamic response of a multilayered poroelastic half-space to time-harmonic surface tractions. In a cylindrical coordinate system, a method of displacement potentials is applied first to decouple the Biot’s wave equations into four scalar Helmholtz equations, and then, general solutions to those equations are obtained. After that, the propagator matrix method and the vector surface harmonics are employed to derive the solutions for a multilayered poroelastic half-space subjected to surface tractions. It is known that the original propagator algorithm has the loss-of-precision problem when the waves become evanescent. At present, an orthogonalization procedure is inserted into the matrix propagation loop to avoid the numerical difficulty of the original propagator algorithm. Finally, a high-order adaptive integration method with continued fraction expansions for accelerating the convergence of the truncated integral is adopted to numerically evaluate the integral solutions expressed in terms of semi-infinite Hankel-type integrals with respect to horizontal wavenumber. Furthermore, to validate the present approach, the response of a uniform poroelastic half-space is examined using the formulation proposed in this article. It is shown that the numerical results computed with this approach agree well with those computed with the analytical solution of a uniform half-space.     

Vertical vibrations of elastic foundation resting on saturated half-space

This paper is mainly concerned with the dynamic response of an elastic foun- dation of finite height bounded to the surface of a saturated half-space.The foundation is subjected to time-harmonic vertical loadings.First,the transform solutions for the governing equations of the saturated media are obtained.Then,based on the assumption that the contact between the foundation and the half-space is fully relaxed and the half- space is completely pervious or impervious,this dynamic mixed boundary-value problem can lead to dual integral equations,which can be further reduced to the Fredholm integral equations of the second kind and solved by numerical procedures.In the numerical exam- ples,the dynamic compliances,displacements and pore pressure are developed for a wide range of frequencies and material/geometrical properties of the saturated soil-foundation system.In most of the cases,the dynamic behavior of an elastic foundation resting on the saturated media significantly differs from that of a rigid disc on the saturated half-space. The solutions obtained can be used to study a variety of wave propagation problems and dynamic soil-structure interactions.     

Analytical solution of the spherical indentation problem for a half-space with gradients with the depth elastic properties

The contact problem for the impression of spherical indenter into a non-homogeneous (both layered and functionally graded) elastic half-space is considered. Analytical methods for solving this problem have been developed. It is assumed that the Lame coefficients vary arbitrarily with the half-space depth. The problem is reduced to dual integral equations. The developed methods make it possible to find the analytical asymptotically exact problem solution, suitable for a PC. The influence of the Lame coefficients variation upon the contact stresses and size of the contact zone with different radius of indenter as well as values of the impressing forces are studied. The effect of the non-homogeneity is examined. The developed method allows to construct analytical solutions with presupposed accuracy and gives the opportunity to do multiparametric and qualitative investigations of the problem with minimal computation time expenditure.     

Anti-plane time-harmonic Green's functions for a circular inhomogeneity in piezoelectric medium with a spring- or membrane-type interface

The present work focuses on the two-dimensional anti-plane time-harmonic dynamic Green's functions for a circular inhomogeneity immersed in an infinite matrix with an imperfect interface, where both the inhomogeneity and matrix are assumed to be piezoelectric and transversely isotropic. Two types of imperfect interface, the spring-type interface with electromechanical coupling and the membrane-type interface, are considered. The former type is often used to model the electromechanical damage along the interface while the latter is largely employed to simulate surface/interface effect of nano-sized inhomogeneity. By using the Bessel function expansions, explicit solutions for the electromechanical fields induced by a time-harmonic anti-plane line force and line charge located in an unbounded matrix as well as the circular inhomogeneity are respectively derived. The present solutions can recover the anti-plane Green's functions for some special cases, such as the dynamic or quasi-static Green's functions of piezoelectricity with perfect interface as well as the dynamic or quasi-static Green's functions of pure elasticity with imperfect interface. For detailed discussions, selected analytical results are presented. Dependence of the electromechanical fields on circular frequency as well as interface properties is examined. The size effect related to interfacial imperfection is also discussed.     

Anti-plane shear Lamb''s problem treated by gradient elasticity with surface energy

The consideration of higher-order gradient effects in a classical elastodynamic problem is explored in this paper. The problem is the anti-plane shear analogue of the well-known Lamb's problem. It involves the time-harmonic loading of a half-space by a single concentrated anti-plane shear line force applied on the half-space surface. The classical solution of this problem based on standard linear elasticity was first given by J.D. Achenbach and predicts a logarithmically unbounded displacement at the point of application of the load. The latter formulation involves a Helmholtz equation for the out-of-plane displacement subjected to a traction boundary condition. Here, the generalized continuum theory of gradient elasticity with surface energy leads to a fourth-order PDE under traction and double-traction boundary conditions. This theory assumes a form of the strain-energy density containing, in addition to the standard linear-elasticity terms, strain-gradient and surface-energy terms. The present solution, in some contrast with the classical one, predicts bounded displacements everywhere. This may have important implications for more general contact problems and the Boundary-Integral-Equation Method.