2-D deformation analysis of a half-space due to a long dip-slip fault at finite depth

Closed form analytical expressions of stresses and displacements at any field point due to a very long dip-slip fault of finite width buried in a homogeneous, isotropic elastic half-space, are presented. Airy stress function is used to derive the expressions of stresses and displacements which depend on the dip angle and depth of the upper edge of the fault. The effect of dip angle and depth of the upper edge of the fault on stresses and displacements is studied numerically and the results obtained are presented graphically. Contour maps for stresses and displacements are also presented. The results of Rani and Singh (1992b) and Freund and Barnett (1976) have been reproduced.     

Elastic solutions for a transversely isotropic half-space subjected to a point load

We rederive and present the complete closed-form solutions of the displacements and stresses subjected to a point load in a transversely isotropic elastic half-space. The half-space is bounded by a horizontal surface, and the plane of transverse isotropy of the medium is parallel to the horizontal surface. The solutions are obtained by superposing the solutions of two infinite spaces, one acting a point load in its interior and the other being free loading. The Fourier and Hankel transforms in a cylindrical co-ordinate system are employed for deriving the analytical solutions. These solutions are identical with the Mindlin and Boussinesq solutions if the half-space is homogeneous, linear elastic, and isotropic. Also, the Lekhnitskii solution for a transversely isotropic half-space subjected to a vertical point load on its horizontal surface is one of these solutions. Furthermore, an illustrative example is given to show the effect of degree of rock anisotropy on the vertical surface displacement and vertical stress that are induced by a single vertical concentrated force acting on the surface. The results indicate that the displacement and stress accounted for rock anisotropy are quite different for the displacement and stress calculated from isotropic solutions. © 1998 John Wiley & Sons, Ltd.     

半空间饱和土中圆形衬砌对弹性压缩波的散射

借助于Biot 波动理论和弹性波的传播理论,采用复变函数和多级坐标法,对半空间饱和土中圆形衬砌结构对弹性稳态压缩波的散射问题进行求解和分析。利用一个半径很大的圆弧来逼近半空间直边界,将待解问题转化为稳态弹性压缩波在一个大圆孔和一个弹性衬砌结构的散射问题。通过引入势函数,将饱和土的Biot波动方程和衬砌的弹性波动方程解耦成Helmholtz 方程,借助复变函数级数展开便可以预先写出该组Helmholtz方程的通解。然后,通过引用复变量,把饱和土和衬砌结构中的应力、位移及孔压用设定的势函数表示出来,再利用半空间饱和土和衬砌结构的连续性条件和近似直边界的圆弧边界和衬砌内边界的边界条件求解出该组势函数的特解。最后,利用势函数的特解,得到饱和土中的位移,应力和孔压及衬砌结构的位移和应力;变换不同的参数求解衬砌结构内外边界的动应力和孔压的集中系数,通过对算例结果的分析得出一系列有益的结论。     

Static elastic deformation in an orthotropic half-space with rigid boundary model due to non-uniform long strike slip fault

The solution of static elastic deformation of a homogeneous, orthotropic elastic uniform half-space with rigid boundary due to a non-uniform slip along a vertical strike-slip fault of infinite length and finite width has been studied. The results obtained here are the generalisation of the results for an isotropic medium having rigid boundary in the sense that medium of the present work is orthotropic with rigid boundary which is more realistic than isotropic and the results for an isotropic case can be derived from our results. The variations of displacement with distance from the fault due to various slip profiles have been studied to examine the effect of anisotropy on the deformation. Numerically it has been found that for parabolic slip profile, the displacement in magnitude for isotropic elastic medium is greater than that for an orthotropic elastic half-space, but, in case of linear slip, the displacements in magnitude for an orthotropic medium is greater than that for the isotropic medium.     

Static deformation of an orthotropic multilayered elastic half-space by two-dimensional surface loads

The transfer matrix approach is used to solve the problem of static deformation of an orthotropic multilayered elastic half-space by two-dimensional surface loads. The general problem is decoupled into two independent problems. The antiplane strain problem and the plane strain problem are considered in detail. Integral expressions for displacements and stresses at any point of the medium due to a normal line load and a shear line load, acting parallel to a symmetry axis, are obtained. In the case of a uniform half-space, closed form analytic expressions for displacements and stresses are derived. The procedure developed is quite easy and convenient for numerical computations.     

层状半空间中洞室对入射平面SV波的散射

利用间接边界元法,在频域内求解了层状弹性半空间中洞室对入射平面SV波的散射问题。通过自由场反应分析,求得假想洞室边界上各点位移和各单元应力响应。在洞室边界各个单元上施加虚拟分布荷载,求得位移和应力的格林函数。根据应力边界条件确定虚拟分布荷载,将自由场位移响应和虚拟分布荷载产生的位移响应叠加起来,即得到问题的解答。比较了层状半空间和均匀半空间中洞室对入射平面SV波的放大作用。结果表明,层状半空间情况有可能导致较大的地表位移幅值,尤其是对于较低频率入射波。     

Plane strain deformation of a multi-layered poroelastic half-space by surface loads

The Biot linearized quasi-static theory of fluid-infiltrated porous materials is used to formulate the problem of the two-dimensional plane strain deformation of a multi-layered poroelastic half-space by surface loads. The Fourier-Laplace transforms of the stresses, displacements, pore pressure and fluid flux in each homogeneous layer of the multi-layered half-space are expressed in terms of six arbitrary constants. Generalized Thomson-Haskell matrix method is used to obtain the deformation field. Simplified explicit expressions for the elements of the 6 × 6 propagator matrix for the poroelastic medium are obtained. As an example of the possible applications of the analytical formulation developed, formal solution is given for normal strip loading, normal line loading and shear line loading.     

Surface load dynamic solution of saturated transversely isotropic multilayer half-space

This paper treats the dynamic response of a multilayered transversely isotropic fluid saturated poroelastic half-space under surface time-harmonic traction. The governing system of partial differential equations is uncoupled with the use of a set of physically meaningful and complete potential functions that decompose different body waves in a saturated poroelastic transversely isotropic medium. After expressing the equations in the Hankel-Fourier domain, a proper algebraic factorization is applied to generate reflection and transmission matrices for decomposed waves. All responses including displacements, stresses, and pore fluid pressure for both general patch load and point load are presented in the form of semi-infinite line integrals. The verification of the method is confirmed with the degeneration of the solutions presented here to the existing solutions for dried both homogeneous and multilayered elastic half-spaces as well as poroelastic half-space. Selected numerical results are depicted to investigate the effects of layering and pore pressure on responses of a transversely isotropic poroelastic medium. The load distribution effects are studied by comparison of the patch and point load responses. Also, resonance notion and effective parameters on this phenomenon such as layering system and anisotropy contrast are discussed. Significant influence of materials and layering configuration on number and amplitude of resonances depicted through the numerical evaluation.     

Deformation of a layered half-space due to a very long tensile fault

The problem of the coseismic deformation of an earth model consisting of an elastic layer of uniform thickness overlying an elastic half-space due to a very long tensile fault in the layer is solved analytically. Integral expressions for the surface displacements are obtained for a vertical tensile fault and a horizontal tensile fault. The integrals involved are evaluated approximately by using Sneddon’s method of replacing the integrand by a finite sum of exponential terms. Detailed numerical results showing the variation of the displacements with epicentral distance for various source locations in the layer are presented graphically. The displacement field in the layered half-space is compared with the corresponding field in a uniform half-space to demonstrate the effect of the internal boundary. Relaxed rigidity method is used for computing the postseismic deformation of an earth model consisting of an elastic layer of uniform thickness overlying a viscoelastic half-space.     

Elastic consolidation around a point sink embedded in a half-space with anisotropic permeability

The complete solution is presented for the transient effects of pumping fluid from a point sink embedded in a saturated, porous elastic half-space. It is assumed that the medium is homogeneous and isotropic with respect to its elastic properties and homogeneous but anisotropic with respect to the flow of pore fluid. The soil skeleton is modelled as a linear elastic material obeying Hooke's law, while the pore fluid is assumed to be incompressible with its flow governed by Darcy's law. The solution has been evaluated for a particular value of Poisson's ratio of the solid skeleton, i.e. 0.25, and the results have been presented graphically in the form of isochrones of excess pore pressure and surface profile for the half-space. The solutions presented may have application in practical problems such as dewatering operations in compressible soil and rock masses.     

Deformation of an elastic layer coupling in different ways to a base due to a very long vertical strike-slip dislocation

The closed-form analytic expressions for the displacement and stresses at any point of an elastic layer lying over a base due to a very long vertical strike-slip dislocation are obtained. The interface between the layer and the base is assumed to be either ‘smooth-rigid’ or ‘rough-rigid’ or ‘welded’. The variations of displacement and stresses with the horizontal distance from the fault for different types of coupling of the layer with the base have been studied. It is found that the displacement for ‘welded interface’ lies between the displacements due to ‘smooth rigid’ and ‘rough-rigid’ interfaces for different positions of the observer and different values of the ratio of rigidities of the layer and half-space.     

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.     

Elastic solutions for a transversely isotropic half-space subjected to buried asymmetric-loads

Elastic closed-form solutions for the displacements and stresses in a transversely isotropic half-space subjected to various buried loading types are presented. The loading types include finite line loads and asymmetric loads (such as uniform and linearly varying rectangular loads, or trapezoidal loads). The planes of transverse isotropy are assumed to be parallel to its horizontal surface. These solutions are directly obtained from integrating the point load solutions in a transversely isotropic half-space, which were derived using the principle of superposition, Fourier and Hankel transformation techniques. The solutions for the displacements and stresses in transversely isotropic half-spaces subjected to linearly variable loads on a rectangular region are never mentioned in literature. These exact solutions indicate that the displacements and stresses are influenced by several factors, such as the buried depth, the loading types, and the degree and type of rock anisotropy. Two illustrative examples, a vertical uniform and a vertical linearly varying rectangular load acting on the surface of transversely isotropic rock masses, are presented to show the effect of various parameters on the vertical surface displacement and vertical stress. The results indicate that the displacement and stress distributions accounted for rock anisotropy are quite different for those calculated from isotropic solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

Vertical vibration of a rigid strip footing on viscoelastic half-space

This paper presents an analytical method for modeling the dynamic response of a rigid strip footing subjected to vertical-only loads. The footing is assumed to rest on the surface of a viscoelastic half-space; therefore, effects of hysteretic soil damping on the impedance of the foundation and the generated ground vibrations are considered in the solution. To solve the mixed boundary value problem, we use the Fourier transform to cast a pair of dual integral equations providing contact stresses, which are solved by means of Jacobi orthogonal polynomials. The resulting soil and footing displacements and stresses are obtained by means of the Fourier inverse transform. The solution provides more realistic estimates of footing impedance, compared to existing solutions for elastic soil, as well as of the attenuation of ground vibrations with distance of the footing. The latter is important for the estimation of machine vibration effects on nearby structures and installations.     

Static deformation of two welded half-spaces due to very long strike-slip dislocations

Closed form analytic expressions for displacement and stresses at any point of either of the two homogeneous, isotropic, perfectly elastic half-spaces in welded contact due to very long strike-slip dislocations are obtained. Both cases of vertical and horizontal strike-slip dislocations are discussed in detail. Variation of the displacement with horizontal distance from the fault and with vertical distance from the interface for a vertical strike-slip fault is studied numerically.     

弹性矩形板下横观各向同性多层地基分析

利用弹性矩形板与多层地基表面的竖向位移协调条件与光滑接触条件,由横观各向同性多层地基应力与位移非耦合的传递矩阵解,推导出弹性矩形板下竖向应力和位移的解析解。在此基础上,编制了相应的程序,并进行了数值计算。计算结果表明：矩形板刚度对板底竖向位移及板中心下的竖向应力有着较为显著的影响;板底竖向位移及板中心下的竖向应力随着板刚度的增加而减小,相同荷载作用下横观各向同性地基与均匀各向同性地基模型的计算结果差异较大,实际工程中很有必要采用更符合土体性质的横观各向同性地基模型     

Deformation of a poroelastic layer overlying an elastic half‐space due to dip‐slip faulting

An analytical solution of the plane strain problem of the deformation of a homogeneous, isotropic, poroelastic layer of uniform thickness overlying a homogeneous, isotropic, elastic half‐space due to two‐dimensional seismic sources buried in the elastic half‐space has been obtained. The integral expressions for the displacements, stresses and pore pressure have been obtained using the stress function approach by applying suitable boundary conditions at the free surface and the interface. The solution obtained is in the Laplace–Fourier transform domain. The case of a vertical dip‐slip line dislocation for the oceanic crust model of Earth is studied in detail. Schapery's formula is used for the Laplace inversion and the extended Simpson's formula for the Fourier inversion. Diffusion of pore pressure in the layer is studied numerically. Contour maps showing the pore pressure in the poroelastic layer have been plotted. The effect of the compressibility of the solid and fluid constituents on pore pressure has also been studied. Copyright © 2015 John Wiley & Sons, Ltd.     

层状场地中局部不均体对平面P波的散射

利用间接边界元方法,在频域内求解了层状场地中局部不均体对平面P波的散射。利用精确的土层动力刚度矩阵进行自由场反应分析,求得位移和应力响应。通过计算虚拟分布荷载的格林影响函数,求得相应位移和应力;根据边界条件确定虚拟分布荷载,最终得到问题的解答。研究了入射P波时,不均体宽度、埋深、厚度、入射角、入射频率度等参数对地表位移幅值的影响,并与相应自由场的结果进行了比较。不均体对P波散射有重要影响,在工程场地地震安全性评价中,应合理考虑局部不均体对场地设计地震动参数确定的影响。     

Static deformation of two welded monoclinic elastic half-spaces due to a long inclined strike-slip fault

Static deformation of two monoclinic elastic half-spaces in welded contact due to a long inclined strike-slip fault situated in one of the half-spaces is studied analytically and numerically. Closed-form algebraic expressions for the displacement at any point of the medium are obtained. The variation of the displacement at the interface with the horizontal distance from the fault is studied. The effect of anisotropy on the displacement field is examined. It is found that while the anisotropy of the source half-space has a significant effect on the displacement at the interface, the anisotropy of the other half-space has only a marginal effect.     

Dynamic Green's functions for an anisotropic poroelastic half-space

The dynamic responses of an anisotropic poroelastic half-space under an internal point load and fluid source are investigated in the frequency domain in this paper. By virtue of Fourier transform and Stroh formalism, the three-dimensional (3D) general solutions of the anisotropic Biot's coupling dynamics equations are derived in the frequency domain. Considering the two surface conditions, permeable and impermeable, the analytical solutions for displacement fields and pore pressure in half-space under a point source (point load or a fluid source) are obtained. When the material properties are isotropic, the numerical results of the poroelastic half-space are in excellent agreement with the existing analytical solutions. For anisotropic half-space cases, numerical results show the strong dependence of the dynamic Green's functions on the material properties.