3D dynamic responses of a multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure

Dynamic Green's function plays an important role in the study of various wave radiation, scattering and soil-structure interaction problems. However, little research has been done on the response of transversely isotropic saturated layered media. In this paper, the 3D dynamic responses of a multi-layered transversely isotropic saturated half-space subjected to concentrated forces and pore pressure are investigated. First, utilizing Fourier expansion in circumferential direction accompanied by Hankel integral transform in radial direction, the wave equations for transversely isotropic saturated medium in cylindrical coordinate system are solved. Next, with the aid of the exact dynamic stiffness matrix for in-plane and out-of-plane motions, the solutions for multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure are obtained by direct stiffness method. A FORTRAN computer code is developed to achieve numerical evaluation of the proposed method, and its accuracy is validated through comparison with existing solutions that are special cases of the more general problems addressed. In addition, selected numerical results for a homogeneous and a layered material model are performed to illustrate the effects of material anisotropy, load frequency, drainage condition and layering on the dynamic responses. The presented solutions form a complete set of Green's functions for concentrated forces (including horizontal load in x(y)-direction, vertical load in z-direction) as well as pore pressure, which lays the foundation for further exploring wave propagation of complex local site in a layered transversely isotropic saturated half-space by using the BEMs.     

横观各向同性饱和地基的三维动力响应

首先引入位移函数,将直角坐标系下横观各向同性饱和土Biot波动方程转化为2个解耦的六阶和二阶控制方程;然后基于双重Fourier变换,求解了Biot波动方程,得到以土骨架位移和孔隙水压力为基本未知量的积分形式的一般解,并用一般解给出了饱和土总应力分量的表达式．在此基础上系统研究了横观各向同性饱和半空间体的稳态动力响应问题,考虑表面排水和不排水两种情况,得到了半空间体在任意分布的表面谐振荷载作用下,表面位移的稳态动力响应,文末给出了算例．     

Source mechanism model for ground motion simulation

A simple analytical model for computing ground motion in a layered half-space due to a buried seismic source is presented in this paper. The buried earthquake source is represented as a distribution of double couples varying in time as a ramp function on the fault plane. The analysis is simplified by first decoupling the governing equations into P-SV and SH problem by a coordinate transformation in the frequency-wave number domain. These two problems are solved separately and the final solution is obtained by the sum of solutions of these individual problems. Explicit expressions for ground motion in a layered half-space due to an impulsive double couple are derived. In the sequel, Green’s function for the displacement field in an infinite medium is also presented. The developed source mechanism model is also demonstrated by simulating ground motion for the Kucth earthquake (M_w = 7.7) of 26th January 2001.     

The interaction of punches on a transversely isotropic half-space

Three-dimensional contact problems on the interaction of two similar punches on an elastic transversely isotropic half-space (five elastic constants) are investigated, when the isotropy planes are perpendicular to the boundary of the half-space. In this connection the stiffness of the half-space boundary depends on the direction. The kernel of the integral equation of the contact problems is represented in a quadrature-free form using the theory of generalized functions. This form of the kernel enables it to be regularized at singular points and enables Galanov's method to be used to solve the contact problem with an unknown contact area.     

Stress field around an arbitrary thin inclusion in a transversely isotropic elastic half-space

We consider a thin flat inclusion of arbitrary shape located inside a transversely isotropic elastic half-space in the plane parallel to its boundary z = 0. An arbitrary tangential displacement is prescribed on the inclusion. The boundary of the half-space is stress-free. We need to find the complete field of stresses and displacements in this half-space. A governing integral equation is derived by the generalized method of images, introduced by the author. The case of circular inclusion is considered as an example. Two methods of solution of the governing integral equation are derived. A detailed solution is presented for the particular cases of radial expansion, torsion and lateral displacement of the inclusion. The solution is also valid for the case of isotropy. The governing integral equation for the case of isotropy is derived.     

Apparatus and applications of a generalized Fourier method for transversally isotropic bodies bounded by a plane and a paraboloid of rotation

We construct new basic solutions of equilibrium equations of a transversally isotropic paraboloid of rotation and obtain summation theorems that express them in terms of cylindrical solutions for a half-space and vice versa. The problem of the action of an axial concentrated force on an elastic transversally isotropic half-space with a fixed inclusion in the form of a paraboloid of rotation is investigated. The problem is solved by a generalized Fourier method and is reduced to a system of integral equations with a Fredholm operator on the condition that intersections of the boundaries of the half-space and inclusion are absent. We consider dependences of stresses on the shape of the paraboloidal inclusion and on the distance between the boundary surfaces and analyze the results of the calculation.     

3D dynamic Green’s functions in a multilayered poroelastic half-space

The complete 3D dynamic Green’s functions in the multilayered poroelastic media are presented in this study. A method of potentials in cylindrical coordinate system is applied first to decouple the Biot’s wave equations into four scalar Helmholtz equations, and then, general solutions to 3D wave propagation problems are obtained. After that, a three vector base and the propagator matrix method are introduced to treat 3D wave propagation problems in the stratified poroelastic half-space disturbed by buried sources. 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. At last, the validity of the present approach for accurate and efficient calculating 3D dynamic Green’s functions of a multilayered poroelastic half-space is confirmed by comparing the numerical results with the known exact analytical solutions of a uniform poroelastic half-space.     

横观各向同性饱和地基与中厚圆板的非轴对称动力相互作用

研究横观各向同性饱和土地基上中厚弹性圆板的非轴对称振动问题,即首先利用Fourier展开和Hankel变换技术,求解了简谐激励下横观各向同性饱和土地基的非轴对称Biot波动方程,然后按混合边值问题建立地基与弹性中厚圆板非轴对称动力相互作用的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程．文末给出了算例．数值结果表明,在一定频率范围内,地基表面的位移幅值随激振频率增加而增大,随距离的增大以振荡形式衰减变化．     

Assessing the mechanical responses for anisotropic multi-layered medium under harmonic moving load by Spectral Element Method (SEM)

The article is concerned with the mechanical responses of anisotropic multi-layered medium under harmonic moving load. An analytical solution for two-dimensional anisotropic multi-layered medium subjected to harmonic moving load is devoted via Spectral Element Method (SEM), while the anisotropic property is approximated as transverse isotropy. Starting with the constitutive equations of transversely isotropic body and the governing equations of motion based on the loading properties. The analytical spectral elements in the wavenumber domain are obtained according to the principle of wave superposition and Fourier transformation. Then, the spectral global stiffness matrix of the multi-layered medium is derived by assembling the nodded stiffness matrices of all layers depended on the different interlayer conditions between the adjacent layers, i.e. sliding and bonded. The corresponding analytical solutions are achieved by taking the Fourier series and Inverse Fast Fourier Transform (IFFT) algorithm. Finally, some examples are given to validate the accuracy of the proposed analytical solution, and to demonstrate the impact of both anisotropy, top layer thickness, interlayer conditions, and loading properties (velocity and natural frequency) on the mechanical response of the multi-layered medium.     

Boussinesq indentation of a transversely isotropic half-space reinforced by a buried inextensible membrane

The normal indentation of a rigid circular disk into the surface of a transversely isotropic half-space reinforced by a buried inextensible thin film is addressed. By virtue of a displacement potential function and the Hankel transform, the governing equations of this axisymmetric mixed boundary value problem are represented as a dual integral equation, which is subsequently reduced to a Fredholm integral equation of the second kind. Two important results of the contact stress distribution beneath the disk region as well as the equivalent stiffness of the system are expressed in terms of the solution of the Fredholm integral equation. When the membrane is located on the surface or at the remote boundary, exact closed-form solutions are presented. For the limiting case of an isotropic half-space the results are verified with those available in the literature. As a special case, the elastic fields of a reinforced transversely isotropic half-space under the action of surface axisymmetric patch loads are also given. The effects of anisotropy, embedment depth of the membrane, and material incompressibility on both the contact stress and the normal stiffness factor are depicted in some plots.     

Combined Fourier-wavelet transforms for studying dynamic response of anisotropic multi-layered flexible pavement with linear-gradual interlayers

The objective of this study is to develop an analytical solution for studying dynamic response of anisotropic multi-layered flexible pavement with linear-gradual interlayers under a harmonic moving load. In this study, the flexible pavement structure is simplified as one multi-layered flexible system, which is assumed to be a semi-infinite medium. A new approach combining the Fourier transform with the wavelet transform for solving the dynamic analytical solution. The wavelet transforms for solving inverse Fourier transform, in solving the solution in the physical domain, is superior to the conventional inverse Fourier transform. The linear-gradual interlayer between the adjacent layers is defined using the shear spring model, and the anisotropic property is simplified as transverse isotropy. Also, in the dynamic analytical solving processes, the motion in the transversely isotropic medium is decoupled into in-plane motion and out-of-plane motion because of the propagation of the waves in a transversely isotropic medium with coupling phenomena. The corresponding analytical solution is entered into a MATLAB-based computer program, which can compute the dynamic responses of an anisotropic multi-layered medium at different interlayer conditions. The accuracy of this program is confirmed through comparison with the results from the examples from the references. The influence analyses of linear-gradual interlayers and anisotropic properties of structural layers are illustrated. It is concluded that the proposed analytical solution-based computer program could be used in the multi-layered flexible pavement structural design and risk management in civil engineering.     

层状横观各向同性饱和地基上圆板的非轴对称振动

研究了层状横观各向同性饱和地基上弹性圆板的非轴对称振动问题．首先,通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的三维动力方程转化为一阶常微分方程组,基于径向Hankel变换,建立问题的状态方程,求解状态方程后得到传递矩阵;其次,利用传递矩阵,结合层状饱和地基的边界条件、排水条件及层间接触和连续条件,给出了任意简谐激振力作用下层状横观各向同性饱和地基动力响应的通解;然后,按混合边值问题建立层状饱和地基上弹性圆板非轴对称振动的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程,并给出了算例．     

Action of a shear stress pulse on the surface of a viscoelastic anisotropic half-space

The boundary value problem of the behavior of a viscoelastic half-space subjected to a surface shear stress pulse is solved in the linear formulation on the basis of dynamical correspondence principle. The medium occupying the half-space possesses the property of transverse isotropy. The exact solution of the problem is obtained using integral transforms. The solution is analyzed with reference to a Maxwell-Voigt model.Moscow Region. Translated from Mekhanika Polimerov, No. 5, pp. 933–937, September–October, 1969.     

Free and forced finite amplitude oscillation of a spherical shell of radially isotropic elastic material

A spherical elastic shell with radial transverse isotropy is considered. The periods of finite amplitude radial oscillation of the shell have been obtained in two cases, namely (i) when both the surfaces of the shell are free from traction, and (ii) when the shell boundaries are uniformly loaded in such a way that the pressure difference between inner and outer surfaces is constant with respect to time. It is observed that for free oscillation to take place it is necessary to impose a new restriction on the strain energy function in addition to those already obtained for finite amplitude oscillation of an isotropic elastic shell. In the forced oscillation case however the required conditions are of the same form as in the corresponding case for isotropic media.     

分析土壤半无限域的一种有限层方法(Ⅰ)——一般理论

本文研究了一种用于横观各向同性体动力学的有限层法.将土壤介质半无限域视为一个横观各向同性半空间体,介质的材料函数沿深度变化,将介质沿深度方向分成若干层,在每一层内材料函数用一个指数函数来模拟.这样,使求解问题的方程得到简化,利用Fourier变换我们得到了各层处“节点”力与位移的关系.这种有限层法实际上是一种半解析法.具有一般半解析法所具备的数据前后处理量少及计算量小等优点.本文所研究的这种有限层法为分析土壤与结构相互作用问题提供了一条新途径.     

Two-dimensional problem of a two-temperature generalized thermoelastic half-space subjected to ramp-type heating

In this work, we will consider a half-space filled with an elastic material, which has constant elastic parameters. The governing equations are taken in the context of the theory of two-temperature generalized thermoelasticity. A linear temperature ramping function is used to more realistically model thermal loading of the half-space surface. The medium is assumed initially quiescent. Laplace and Fourier transform techniques are used to obtain the general solution for any set of boundary conditions. The general solution obtained is applied to a specific problem of a half-space subjected to ramp-type heating. The inverse Fourier transforms are obtained analytically while the inverse Laplace transforms are computed numerically using a method based on Fourier expansion techniques. Some comparisons have been shown in figures to estimate the effect of the ramping parameter of heating.     

横观各向同性弹性层点力解

本文根据弹性层状结构的传递矩阵法思想，由横观各向同性弹性力学基本方程，导出了含应力和位移两类变量的混合方程，利用Ｆｏｕｒｉｅｒ变换和文献［７］的位移函数通解，以及计算机代数软件，得到了横观各向同性层的点力解，这个点力解可直接退化到各同性情形的解．     

Resultant forces and moments in mixed-mixed problems of the theory of elasticity

A problem is called mixed-mixed, when both normal and tangential displacements are prescribed on a part of the boundary, while the normal and tangential stresses are prescribed at the rest of the boundary. Exact closed form expressions have been derived for the resultant normal and tangential forces, tilting moment and torque, directly through the prescribed displacements, thus eliminating the need for determination of stresses. The problem solved treats a transversely isotropic elastic half-space, with arbitrary normal and tangential displacements prescribed inside a circle, and the rest of the boundary being stress-free. The interaction between an arbitrary force inside the half-space and a bonded punch is considered as an example. No similar result has ever been reported, even in the case of isotropy.     

Transient response of SH waves in a layered half-space with sub-surface and interface cracks

The transient scattering of SH waves by sub-surface and interface cracks parallel to the free surface in a layered elastic solid is investigated. The problem in frequency domain is solved by using a hybrid method which combines the finite element method of the near field with the boundary integral representation of the far field. The transient responses are then obtained by inverting the spectra via fast Fourier transform with the incident pulse Ricker of wavelet. Numerical results are presented for the surface displacements, dynamic stress intensity factors and wave motion in the layered half-space. Furthermore, the propagations of reflected, diffracted, and direct impact waves at any instant are clearly identified by the present method. To understand the mechanism of elastic wave interaction is very important in the field of ultrasonic non-destructive evaluation (NDE) and fracture mechanics studies.     

横观各向同性饱和弹性多孔介质非轴对称动力响应

应用Fourier展开和Hankel变换求解了简谐激励下横观各向同性饱和弹性多孔介质的非轴对称Biot波动方程,得到了一般解。用一般解给出了多孔介质总应力分量的表达式。最后对求解横观各向同性饱和弹性多孔介质非轴对称动力响应边值问题的方法作了系统说明,并且给出了数值分析特例。