The Green￠s function of the harmonic horizontal force applied at the interior of the saturated half-space soil

By using integral transform methods, the Green's functions ofhorizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous con-ditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green's functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper.     

The Green’s function of the harmonic horizontal force applied at the interior of the saturated half-space soil

By using integral transform methods, the Green’s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above-mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green’s functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper. Foundation item: State Natural Science Foundation (59879012) and Doctoral Foundation from State Education Commission (98024832).     

半空间饱和土在内部简谐水平力作用下的Ge函数

根据积分变换方法得出了半空间内部作用简谐水平力时的Gree函数．首先，利用Hankel积分变换方法，直接对频域内的Biot波动方程进行求解，得出Biot波动方程的通解；利用通解和半空间内部作用水平力时边界上的边界条件，以及力作用面上的连续性条件，可以得出上述边值问题的解；对于边值问题在变换域内的解进行相应的逆变换，就可以得出频域内的Gree函数．本文得到的线弹性退化解与文献中的结果吻合．最后，文中给出了两个算例．      

Solution of dynamic Green׳s function for unsaturated soil under internal excitation

The closed form three-dimensional Green׳s function of a semi-infinite unsaturated poroelastic medium subjected to an arbitrary internal harmonic loading is derived, with consideration of capillary pressure and dynamic shear modulus varying with saturation. By applying the Fourier expansion techniques and Hankel integral transforms to the circumferential and radial coordinates, respectively, the general solution for the governing partial differential equations is obtained in the transformed domain. A corresponding boundary value problem is formulated. The integral solutions for the induced displacements, pore pressure and net stress are then determined considering the continuity conditions. The formulas are compared with the degenerated solution of saturated soils and confirmed. Numerical results reveal that the response of the unsaturated half-space depends significantly on the saturation by altering dynamic shear modulus to account for the effects of matric suction on soil stiffness. Slight differences between the results occur if only the saturation is taken into account. Moreover, a large source-depth results in a pronounced contribution to the reduction of surface displacement amplitudes. The analytical solutions concluded in the study offer a broader application to dynamic response associated with axi-symmetric and asymmetric conditions.     

Horizontal vibrations of a disk on a poroelastic half-space

This paper analytically examines the horizontal vibration of a rigid disk on a saturated poroelastic half-space. The pressure-solid displacement form of the harmonic equations of motion for asymmetric dynamic problem are developed from the form of the equations originally presented by Biot. Making use of a new method the solution of the above equations is obtained. According to the mixed boundary -value conditions, the dual integral equations of the horizontal vibration of a rigid disk on a saturated poroelastic half-space are established. By appropriate transforms, it is shown that the dual integral equations can be reduced to a pair of Fredholm integral equations of the second kind, whose solutions are then computed. Numerical results for the horizontal dynamic compliance coefficient are given at the end of this paper.     

Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

Vertical dynamic response of a disk on a saturated poroelastic half-space

This paper considers the vertical dynamic response of a disk on a saturated poroelastic half-space. Firstly the pressure-solid displacement form of the harmonic equations of motion for a poroelastic solid are developed from the form of the equations originally presented by Biot. These equations are solved by a new method. Then the mixed boundary value problem for the vertical harmonic vibration of a disk on a poroelastic half-space is studied. The two types of drainage conditions at the surface of the poroelastic half-space are considered: (a) the surface of the poroelastic half-space is assumed to be completely pervious both within and exterior to the plate; (b) The interface between the plate and the poroelastic half-space is assumed to be impervious and the exterior region is assumed to be pervious. By using the Hankel transform techniques, the paper develops the governing dual integral equations. These governing integral equations are further reduced to systems of standard Fredholm integral equations of the second kind by Abel transform.     

Zero-stress,cylindrical wave functions around a circular underground tunnel in a flat,elastic half-space: Incident P-waves

The zero-stress boundary conditions at the surface of the half-space in the presence of surface and sub-surface cavities for in-plane, incident cylindrical P- and SV-waves have always posed challenging problems. The outgoing cylindrical P- and SV-waves can be represented by Hankel functions of radial distance coupled with the sine and cosine functions of angle. Together, at the half-space surface the P- and SV-wave functions are not orthogonal over the semi-infinite radial distance from 0 to infinity. Thus, to simultaneously satisfy the zero in-plane, normal, and shear stresses, an approximation of the geometry is often made. This paper presents an analytical formulation of the boundary-valued problem, where the Hankel wave functions are expressed in integral form, changing the representation from cylindrical to rectangular coordinates, so that the zero-stress boundary conditions at the half-space surface can be applied in a more straightforward way.     

Dynamic response of a flexible plate on saturated soil layer

An analytical approach is developed to study the dynamic response of a flexible plate on single-layered saturated soil. The analysis is based on Biot's two-phased theory of poroelasticity and also on the classical thin-plate theory. First, the governing differential equations for saturated soil are solved by the use of Hankel transform. The general solutions of the skeleton displacements, stresses, and pore pressures, derived in the transformed domain, are subsequently incorporated into the imposed boundary conditions, which leads to a set of dual integral equations describing the corresponding mixed boundary value problem. These governing integral equations are finally reduced to the Fredholm integral equations of the second kind and solved by standard numerical procedures. The accuracy of the present solution is validated via comparisons with existing solutions for an ideal elastic half-space. Furthermore, some numerical results are presented to show the influences of the layer depth, the plate flexibility, and the soil porosity on the dynamic compliances.     

In-plane dynamic Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic half-space

The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic(TI)half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf(1985),and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.     

Time harmonic point load and dynamic contact problem of contacting fluid and poroelastic half-spaces

The dynamic response of contacting fluid and fluid-saturated poroelastic half- spaces to a time-harmonic vertical point force or a point pore pressure is investigated. The solutions are formulated using the boundary conditions at the fluid-porous medium interface. The point load solutions are then used to solve the dynamic problem of the vertical vibration of a rigid disc (both permeable and impermeable discs are included) on the surface of the poroelastic half-space. The contact problems are solved by integrating the point force and point pore pressure solutions over the contact area with unknown discontinuous force and pore pressure distributions, which are determined from the boundary conditions. The solutions are expressed in terms of dual integral equations, which are converted to Fredholm integral equations of the second kind and solved numerically. Selected numerical results for the vertical dynamic compliance coefficient for the cases with or without fluid overlying the poroelastic half-space are presented to show the effects of the fluid. The influence of the permeability condition of the disc on the compliance of the poroelastic half-space is investigated. The displacement, vertical stress, pore pressure in the poroelastic half-space and water pressure in the fluid half-space are also examined for different poroelastic materials and frequencies of excitation. The present results are helpful in the study of the dynamic response of foundations on the seabed under seawater.     

A new series solution method for two-dimensional elastic wave scattering along a canyon in half-space

This paper presents a semi-analytical method for studying the two-dimensional problem of elastic wave scattering by surface irregularities in a half-space. The new method makes use of the member of a c-completeness family of wave functions to construct the scattering fields, and then applies equal but opposite tractions to those of the foregoing constructed scattering fields on the horizontal surface of the half-space to produce additional scattering fields. These additional scattering fields are a series of Lamb's solutions. Thus the whole scattering field constructed in the series automatically satisfies the Navier equations, the condition of zero traction on the half-space surface, and the radiation boundary conditions at infinity. Using the traction-free conditions along the canyon surface, the coefficients of the series solutions are determined via a least-squares method. For incident P, SV, and Rayleigh waves, the numerical results are presented for the scattering displacements in the vicinity of a semi-circular canyon in the half-space.     

Vertical and horizontal vibrations of a rigid disc on a multilayered transversely isotropic half-space

A half-space containing horizontally multilayered regions of different transversely isotropic elastic materials as well as a homogeneous half-space as the lowest layer is considered such that the axes of material symmetries of different layers and the lowest half-space to be as depth-wise. A rigid circular disc rested on the free surface of the whole half-space is considered to be under a forced either vertical or horizontal vibration of constant amplitudes. Because of the involved integral transforms, the mixed boundary value problems due to mixed condition at the surface of the half-space are changed to some dual integral equations, which are reduced to Fredholm integral equations of second kind. With the help of contour integration, the governing Fredholm integral equations are numerically solved. Some numerical evaluations are given for different combinations of transversely isotropic layers to show the effect of degree of anisotropy of different layers on the response of the inhomogeneous half-space.     

Dynamic displacements and stresses in the vicinity of a cylindrical cavity embedded in a half-space

The two-dimensional response of a viscoelastic half-space containing a buried, unlined, infinitely long cylindrical cavity of circular cross-section subjected to harmonic plane SH, P, SV and Rayleigh waves is obtained by an indirect boundary integral method based on the two-dimensional Green's functions for a viscoelastic half-space. An extensive critical review of the existing numerical results obtained by other techniques is presented together with some new numerical results describing the motion on the ground surface and the motion and stresses on the wall of the cavity for P, SV, SH and Rayleigh waves.     

Two dimensional transient fundamental solution due to suddenly applied load in a half-space

A transient Green function due to suddenly applied line loads in an isotropic and homogeneous half-space is reported in this paper. The derivation of the half-space Green function in the Laplace and the Fourier transform spaces is first reviewed. Following an explicit inversion of the Fourier transform, the inverse Laplace transform is implemented along the contour integral on the p-complex plane in an integral form. The half-space Green function consists of full-space Green functions and a singularity-free complementary term. It can be easily incorporated into current transient boundary elements using the transient full-space Green function. Combined with finite elements, the half-space Green function can be used in a hybrid procedure to solve transient half-space problems without discretization of the free surface. Numerical results are presented to illustrate transient wave propagation in a half-space.     

Seismic response of a cylindrical shell embedded in a layered viscoelastic half-space. I: Formulation

A method to obtain the three-dimensional harmonic response of a infinitely long cylindrical shell of circular cross-section embedded in a layered viscoelastic half-space and subjected to harmonic plane waves impinging at an oblique angle with respect to the axis of the shell is presented. The procedure combines an indirect integral representation for the field in the exterior half-space with a model of the pipeline or tunnel based on Donnell shell theory. The integral representation for the soil is based on the use of moving Green's functions for the layered viscoelastic half-space. The accuracy of the formulation is tested by comparison of results obtained by using different discretizations. Extensive comparisons with previous two- and three-dimensional results for the case of a shell embedded in a uniform half-space and some new numerical results for a shell embedded in a multilayered half-space are presented in a companion paper.     

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅰ): Formulation

This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.     

Diffraction of obliquely incident waves by a cylindrical cavity embedded in a layered viscoelastic half-space

The three-dimensional harmonic response in the vicinity of an infinitely long, cylindrical cavity of circular cross-section buried in a layered, viscoelastic half-space is obtained when the half-space is subjected to homogeneous plane waves and surface waves impinging at an oblique angle with respect to the axis of the cavity. The solution is obtained by an indirect boundary integral method based on the use of moving Green's functions for the viscoelastic half-space. Numerical results describing the motion on the ground surface and the motion and stresses on the wall of the cavity are presented for obliquely incident P-, SV-, SH- and Rayleigh waves with different horizontal angles of incidence.     

Impedances of embedded rigid strip foundations

This paper is concerned with the dynamic response of rigid strip foundations of arbitrary geometry embedded in a homogeneous elastic half-space. The embedded rigid foundation is modelled by an equivalent domain in a uniform half-space which is subjected to an appropriate body force field. The components of the impedance matrix are determined through the solution of a linear simultaneous equation system which is established by invoking rigid body displacements of discrete locations within the equivalent domain and appropriate equilibrium consideration. It is found that high numerical efficiency and flexibility can be achieved using the body force model when compared to boundary integral formulations through the selection of appropriate displacement influence functions and a ‘parent domain’ in the analysis. Numerical results are presented to illustrate the influence of the embedment ratio, frequency of excitation, foundation geometry and Poisson's ratio on the vertical, horizontal, rocking and coupled impedances of a single embedded foundation. The effect on the impedance due to the presence of an adjacent embedment is investigated for various distances between foundations and embedment ratios.     

横观各向同性饱和地基中埋置荷载的非轴对称瞬态响应

基于孔隙介质的Biot理论,首先利用Laplace变换,给出圆柱坐标系下横观各向同性饱和弹性多孔介质在变换域上的波动方程;将波动方程解耦后,根据方位角的Fourier展开和径向Hankel变换,求解了Biot波动方程,得到以土骨架位移、孔隙水压力和土介质总应力分量的积分形式的一般解;借助一般解,建立了有限厚度饱和土层和饱和半空间的精确动力刚度矩阵,并由土层的层间界面连续条件建立三维非轴对称层状饱和地基的总刚度方程;在此基础上,系统研究了横观各向同性饱和半空间体在内部集中荷载激励下的动力响应,并给出了问题的瞬态解答.该研究为运用边界元法求解饱和地基动力响应奠定了理论基础.