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.     

Diffraction of plane SV waves by a shallow circular-arc canyon in a saturated poroelastic half-space

An analytical solution for the scattering and diffraction of incident plane SV waves by a shallow circular-arc canyon in a saturated poroelastic half-space is derived by the wave function expansion method. The solution is utilized to analyze the dependence of the computed surface motions on the incident frequencies, incident angles, porosity, boundary drainage and Poisson's ratio. It is shown that, depending on the incident angles, the surface displacement amplitudes around a canyon in a dry poroelastic half-space and saturated poroelastic half-space can be very different. The surface displacement amplitudes of an undrained saturated poroelastic half-space are close to those of a drained saturated poroelastic half-space. For low porosity, the surface displacement amplitudes of a saturated poroelastic half-space are almost identical to those of a dry poroelastic half-space, and drainage condition has little influence on the surface displacement amplitudes. But for high porosity, the effect of drainage condition becomes significant, and for the same porosity, the displacement amplitudes of an undrained saturated half-space will be larger than those of a drained saturated half-space. Poisson's ratio is also an important factor affecting the surface displacement amplitudes around the canyon, both in drained and undrained conditions, but leads to larger effects for an undrained saturated half-space than for a drained saturated half-space. Large pore pressures are found around the canyon and their amplitudes depend on the incident angles and frequencies. Below the surface, the amplitudes of pore pressures are less than they are at the surface, especially for high frequencies.     

Diffraction of plane SV waves by a cavity in poroelastic half-space

This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green’s functions of compressive and shear wave sources are derived based on Biot’s theory. The scattered waves are constructed using fi ctitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verifi cation of the accuracy is performed by: (1) checking the satisfaction extent of the boundary c...     

Dynamic responses of a saturated poroelastic half-space generated by a moving truck on the uneven pavement

A truck–pavement–ground coupling model was established to study the dynamic responses of a saturated poroelastic half-space generated by a moving heavy truck on the uneven pavement. The ground was simulated as a fully saturated poroelastic half-space governed by Biot’s theory. The overlying pavement was simplified as a Kirchhoff thin plate. With the assumption of a sinusoidal pavement surface, the dynamic wheel–pavement force was obtained through a linear Hertizian contact model. The numerical results showed that this dynamic load could make considerable contributions to the stress and excess pore water pressure responses in the ground. Furthermore, the effective stress path of the soil unit beneath the pavement caused by the moving truck was firstly calculated and presented. It was found that the differences between the total stress path and the effective stress path became significant as the truck speed increased, thus the effective stress path was more suitable than total stress path to reflect the stress history of soil elements in the saturated ground during the passage of high-speed traffics.     

Scattering of plane SV waves by cylindrical canyons in saturated porous medium

On the basis of Biot dynamic theory, an analytic solution of two-dimensional scattering and diffraction of plane SV waves by circular cylindrical canyons in a half space of saturated porous media is presented in this paper for the first time. The solution is obtained by employing the Fourier–Bessel series expansion technique. Parametric studies had been carried out, which includes: the angle of incidence, the frequency of the incident SV wave, the porosity of saturated porous medium and the stiffness and Poisson's ratio of the solid-skeleton. All the outcomes are useful for the seismic analysis of the surface topography 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.     

A viscous-spring transmitting boundary for cylindrical wave propagation in saturated poroelastic media

Based on the u–U formulation of Biot equation and the assumption of zero permeability coefficient, a viscous-spring transmitting boundary which is frequency independent is derived to simulate the cylindrical elastic wave propagation in unbounded saturated porous media. By this viscous-spring boundary the effective stress and pore fluid pressure on the truncated boundary of the numerical model are replaced by a set of spring, dashpot and mass elements, and its simplified form is also given. A u–U formulation FEA program is compiled and the proposed transmitting boundaries are incorporated therein. Numerical examples show that the proposed viscous-spring boundary and its simplified form can provide accurate results for cylindrical elastic wave propagation problems with low or intermediate values of permeability or frequency content. For general two dimensional wave propagation problems, spuriously reflected waves can be greatly suppressed and acceptable accuracy can still be achieved by placing the simplified boundary at relatively large distance from the wave source.     

Viscoelastic representation of surface waves in patchy saturated poroelastic media

Wave-induced flow is observed as the dominated factor for P wave propagation at seismic frequencies. This mechanism has a mesoscopic scale nature. The inhomogeneous unsaturated patches are regarded larger than the pore size, but smaller than the wavelength. Surface wave, e.g., Rayleigh wave, which propagates along the free surface, generated by the interfering of body waves is also affected by the mesoscopic loss mechanisms. Recent studies have reported that the effect of the wave-induced flow in wave propagation shows a relaxation behavior. Viscoelastic equivalent relaxation function associated with the wave mode can describe the kinetic nature of the attenuation. In this paper, the equivalent viscoelastic relaxation functions are extended to take into account the free surface for the Rayleigh surface wave propagation in patchy saturated poroelastic media. Numerical results for the frequency-dependent velocity and attenuation and the time-dependent dynamical responses for the equivalent Rayleigh surface wave propagation along an interface between vacuum and patchy saturated porous media are reported in the low-frequency range (0.1–1,000 Hz). The results show that the dispersion and attenuation and kinetic characteristics of the mesoscopic loss effect for the surface wave can be effectively represented in the equivalent viscoelastic media. The simulation of surface wave propagation within mesoscopic patches requires solving Biot’s differential equations in very small grid spaces, involving the conversion of the fast P wave energy diffusion into the Biot slow wave. This procedure requires a very large amount of computer consumption. An efficient equivalent approach for this patchy saturated poroelastic media shows a more convenient way to solve the single phase viscoelastic differential equations.     

Diffraction of plane P waves around an alluvial valley in poroelastic half-space

Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green’s functions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of fictitious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary...     

平面SV波在饱和半空间中沉积谷地周围的散射

采用一种特殊的间接边界积分方程法,求解了平面SV波在饱和半空间中任意形状沉积谷地周围的二维散射问题。结合饱和半空间中膨胀波源和剪切波源格林函数,由分布在沉积和半空间交界面附近两虚拟波源面上的波源分别构造沉积内外的散射波场,由交界面连续条件建立方程并求解确定虚拟波源密度,总波场反应即可由自由波场和散射波场叠加而得。然后通过边界条件验算、退化解答与现有结果的比较以及稳定性检验,验证了方法的计算精度。通过一组典型算例,研究了平面SV波在饱和半空间中沉积谷地周围散射的基本规律,详细给出了不同参数情况沉积谷地附近地表位移幅值和孔隙水压,着重分析了入射SV波频率和角度、边界渗透条件、沉积孔隙率等因素对场地反应的影响,得出了一些有益的结论。     

Scattering of plane P waves by circular-arc alluvial valleys with saturated soil deposits

This paper presents an analytical solution for two-dimensional scattering and diffraction of plane P waves by circular-arc alluvial valleys with shallow saturated soil deposits. The solution is based on Biot's dynamic theory for saturated porous media, and derived by employing Fourier–Bessel series expansion technique. In this analysis, soil deposits in the circular-arc valley are modeled as saturated porous media based on Biot's dynamic theory, and the circular-arc valley is assumed to be imbedded in an infinite half-space, filled with elastic single-phase media. Numerical results from this solution show that the amplitudes of displacement at the surface of an alluvial valley are mainly relative to the angle of incidence, the dimensionless frequency of incident P wave, the degree of saturation and porosity of soil deposits, and the stiffness and Poisson's ratio of the solid skeleton of the soil deposits. Furthermore, the proposed solution is compared with the previous solution, in which the soil deposit was modeled as an elastic single-phase solid.     

平面P波在饱和半空间中洞室周围的散射(II):数值结果

本文通过数值计算研究了入射平面P波在饱和半空间中洞室周围散射问题,分析了入射波频率和角度、边界渗透条件、孔隙率、泊松比等参数对散射的影响。研究表明,平面P波入射情况下,饱和半空间和单相(干土)半空间中洞室附近地表位移幅值的差别很大,干土情况的水平位移幅值相对较大,饱和情况的竖向地表位移幅值相对较大;由于波在洞室附近的干涉,饱和情况与干土情况的地表位移出现相位漂移。随着孔隙率的增大,洞室附近水平地表位移幅值逐渐减小,竖向地表位移幅值则逐渐增大;当孔隙率较低时,边界渗透条件对地表位移幅值的影响很小,而当孔隙率较大时,边界渗透条件的影响则不可忽视,不透水情况下,水平和竖向地表位移幅值的峰值均相对较大;随着入射频率的增加,孔隙率的影响逐渐增大,而且不透水情况下孔隙率的影响相对较大。随着泊松比的增大,洞室附近水平地表位移幅值逐渐降低,竖向地表位移幅值则逐渐增大;泊松比较小时,边界渗透条件对位移幅值的影响较大,泊松比较大时,边界渗透条件对位移幅值的影响则较小;随着入射频率的增加,泊松比的影响逐渐增大。当孔隙率较小时,半空间地表和洞室表面孔隙水压幅值较小,但空间变化比较剧烈,随着孔隙率的增大,孔隙水压逐渐增大但空间变化逐渐平缓;随着入射频率的增加,孔隙水压幅值逐渐增大,且孔隙水压的空间变化逐渐变得复杂。     

Propagation of partially coherent non-stationary random waves in a viscoelastic layered half-space

This paper presents a highly accurate method based on the precise integration method (PIM) and on the pseudo excitation method (PEM). The method computes the propagation behaviour of partially coherent non-stationary random waves in a viscoelastic, transversely isotropic solid, which consists of a multi-layered soil resting on a homogeneous semi-infinite space. The excitation source is a local rupture between two layers, which causes a partially coherent non-stationary random field. The analysis of non-stationary random wave propagation is transformed into that for deterministic waves by using PEM. The resulting governing equations in the frequency-wavenumber domain are linear ordinary differential equations, which are solved very precisely by using PIM. The evolutionary power spectral densities of the ground level responses are investigated and some typical earthquake phenomena are explained.     

Horizontal vibration of a rigid disk buried in a poroelastic half-space in contact with a fluid half-space

This paper addresses the horizontal vibration of a rigid disk embedded in a poroelastic half-space in contact with a fluid half-space using the poroelastic theory of potentials. The solution of this problem is expressed in terms of dual integral equations that are converted into Fredholm integral equations of the second kind and solved numerically. Selected numerical results for the horizontal dynamic impedance coefficient are examined based on different poroelastic materials, embedment depths, and excitation frequencies; furthermore, the results are analyzed for the cases in which there is and is no fluid overlying the poroelastic medium to examine the effect of fluid. The results of this study are helpful for designing a foundation embedded in the seabed due to dynamic horizontal forces.     

平面P波在饱和半空间中洞室周围的散射(I):解析解

利用波函数展开法给出了入射平面P波在饱和半空间中圆形洞室周围散射问题的一个解析解。半空间假定为无粘性流体饱和介质,满足Biot理论。采用一种基于实验数据的孔隙率和模量之间的线性关系来确定Biot模型中的介质参数。解答考虑了透水边界和非透水边界两种情况。对边界条件进行了数值检验,结果表明,随着级数截断项数的增大,边界残量衰减很快。解答为进一步研究入射波频率和角度、边界渗透条件、孔隙率、泊松比等参数对散射的影响奠定了基础。     

Dynamic response of a partially debonded pipeline embedded in a saturated poroelastic medium to harmonic plane waves

A practical model of partially debonded pipeline embedded in a saturated poroelastic medium is proposed, and the dynamic response of this model to harmonic plane waves is theoretical investigated. Biot׳s poroelastic theory is introduced to describe the dynamic equations of the saturated poroelastic medium, and the potentials obtained from Helmholtz decomposition theorem are expressed by wave function expansion method. The debonding areas around the pipeline are assumed to be filled with water. The disturbed solutions of basic field equations in these areas are expressed in terms of a scalar velocity potential. Different boundary conditions of bonded and debonded areas are adopted, and the expanded coefficients are obtained. An example of one partially debonded area is presented and analyzed. It is found that the stresses in the perfectly bonded and debonded areas show great difference, and the jump of dynamic stress at the connection points between these two areas is great in the case of low frequency. The effect of debonded areas on the dynamic stress under different thicknesses of lining is also examined.     

Dynamic response of a tunnel buried in a saturated poroelastic soil layer to a moving point load

The dynamic response of a tunnel buried in a two-dimensional poroelastic soil layer subjected to a moving point load was investigated theoretically. The tunnel was simplified as an infinite long Euler–Bernoulli beam, which was placed parallel to the traction-free ground surface. The saturated layer was governed by Biot’s theory. Combined with the specified boundary conditions along the beam and saturated poroelastic layer, the coupled equations of the system were solved analytically in the frequency–wavenumber domain based on Fourier transform. The time domain responses were obtained by the fast inverse Fourier transform. The critical velocity of the considered structure was determined from the dispersion curves. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium subjected to the underground moving load were investigated. It is concluded that, for coarse materials or fine materials subjected to the high-velocity loading, models ignoring the coupling effects between the pore fluid and the soil skeleton may cause errors. The shear modulus and the permeability coefficients of the saturated soil as well as the load moving velocity had significant influence on the displacement and pore pressure responses.     

Diffraction of plane sh waves in a half-space

Scattering of antiplane shear waves (SH) in two dimensions by surface and near-surface defects in a homogeneous, isotropic elastic semi-infinite medium has been studied. Attention has been focused here in the range of medium to long wavelengths. A combined finite element and analytical technique has been used to study the problems of scattering by semi-circular and triangular canyons. The results for the former case are compared with the known exact solution and those for the latter case are compared with some available approximate solutions. Finally a problem of multiple scattering by a triangular canyon and a nearby circular tunnel is studied. Numerical results are presented showing the effects of multiple scattering and different angles of incidence. These results are of interest in earthquake engineering.     

Love waves in an inhomogeneous fluid saturated porous layered half-space with linearly varying properties

The paper is concerned with the propagation of the Love waves in an inhomogeneous transversely isotropic fluid saturated porous layered half-space with linearly varying properties. The analysis is based on Biot's theory. Firstly, the dispersion equation in the complex form for the Love waves in an inhomogeneous porous layer is derived. Then the equation is solved by an iterative method. Detailed numerical calculation is presented for an inhomogeneous fluid saturated porous layer overlying a purely elastic half-space. The dispersion and attenuation of the Love waves are discussed. In addition, the upper and lower bounds of the Love wave speed are explored.     

Elastic waves in visco-elastic half-space generated by various vehicle loads

In this paper, the response of a visco-elastic half-space subjected to moving loads with static and dynamic components is investigated. Four types of vehicle loads are considered, including the moving point load, uniformly distributed wheel load, elastically distributed wheel load, and a train load simulated as a sequence of elastically distributed wheel loads. In each case, the influence of the moving loads traveling in the subsonic, transonic and supersonic ranges on the dynamic responses of the half-space is studied. The parametric study conducted herein enables us to grasp insight into the mechanism of wave propagation for a visco-elastic half-space under moving loads.