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.     

Scattering of plane P,SV or Rayleigh waves by a shallow lined tunnel in an elastic half space

Based on the plane complex variable theory and the image technique, an analytical solution is presented for scattering of plane harmonic P, SV or Rayleigh waves by a shallow lined circular tunnel in an elastic half space. The major contribution of this study is the treatment of the orthogonality of the boundary conditions along the half surface and the cavity wall. In terms of the image technique, the scattered waves by the half surface are simulated as transmitting from the image source of the origin of the tunnel. Using two different conformal mapping functions, we obtained the complex-valued stresses and displacements of the elastic medium and the liner in the image domain, respectively. The boundary value problem results in a set of infinite algebraic equations. The accuracy of the present approach is verified by comparing the present solution results with the available published data. Parametric study indicates that the embedment depth, the shear modulus and the thickness of the liner have significant influences on the dynamic response of the liner and the medium.     

Application of Biot's poroelasticity to some soil dynamics problems in civil engineering

This review type of paper shows how the poroelastodynamic theory of Biot can be applied to some soil dynamics problems encountered in transportation engineering, which have been solved by the present authors. These problems involve rigid walls retaining poroelastic soil and subjected to harmonic seismic waves and moving loads on poroelastic soil. Both classes of problems involve a soil layer over bedrock, are of the plane strain type and are solved analytically by two methods: a direct (almost exact and exact for the above two classes of problems) method and an approximate method. The effects of shear modulus, porosity, permeability and hysteretic damping of the soil medium as well as the seismic frequency for retaining walls and velocity for moving loads on the dynamic response are numerically evaluated in order to assess their relative importance on that response.     

Vertical vibration of an embedded rigid foundation in a poroelastic soil

This paper considers time-harmonic vertical vibration of an axisymmetric rigid foundation embedded in a homogeneous poroelastic soil. The soil domain is represented by a homogeneous poroelastic half space that is governed by Biot's theory of poroelastodynamics. The foundation is subjected to a time-harmonic vertical load and is perfectly bonded to the surrounding half space. The contact surface can be either fully permeable or impermeable. The dynamic interaction problem is solved by employing an indirect boundary integral equation method. The kernel functions of the integral equation are the influence functions corresponding to vertical and radial ring loads, and a ring fluid source applied in the interior of a homogeneous poroelastic half space. Analytical techniques are used to derive the solution for influence functions. The indirect boundary integral equation is solved by using numerical quadrature. Selected numerical results for vertical impedance of rigid foundations are presented to demonstrate the influence of poroelastic effect, foundation geometry, hydraulic boundary condition along the contact surface and frequency of excitation.     

饱水介质波动分析的谐波输入初始间断的处理方法

本文以饱水两相介质的土力学模型为研究对象,在假定两相介质为弹性介质条件下,采用了显式有限元法和透射边界进行了饱和弹性半空间动力响应问题的研究。为避免谐波输入初始间断的影响,文中提出了一个处理函数,并以弹性半空间为算例,对饱水介质和单相介质分别进行了在底边界P波垂直入射时的动力响应分析,验证了该处理函数的有效性和实用性。     

Analytical and numerical analysis of pressure drawdown in a poroelastic reservoir with complete overburden effect considered

Starting from an analytical reservoir model that incorporates full interaction with an elastic overburden, a new hybrid mathematical approach is developed by combining two numerical discretization methods. A tabular reservoir (petroleum reservoir or an aquifer) in an infinite or semi-infinite domain is viewed as a macroscopic displacement discontinuity, allowing use of the efficient displacement discontinuity mathematical method to calculate stresses and displacements that arise because of pressure changes. A 3-D finite element method using a poroelastic formulation is used to discretize the reservoir itself. By coupling the displacement discontinuity and finite element methods, a 3-D large-scale poroelastic reservoir can be simulated within an infinite or semi-infinite domain. The numerical model has been verified through comparison to known solutions, and some time-dependent pressure drawdown problems are analyzed. Results indicate that including the complete overburden (reservoir surroundings) response has a significant effect on pressure drawdown in a poroelastic reservoir during pumping, and should be incorporated in appropriate applications such as well test equations and subsidence analyses.     

Dynamic response of a piecewise circular tunnel embedded in a poroelastic medium

Recently, considerable efforts have been devoted to evaluation of seismic dynamic response of a circular tunnel. Conventional approaches have considered integral liners embedded in an elastic medium. In this study, we re-examine the problem with piecewise liners embedded in a porous medium. Surrounding saturated porous medium of tunnels is described by Biot's poroelastic theory, while the liner pieces and the connecting joints are treated as curved beams and characterized by curved beam theories. The scattered wave field in the porous medium is obtained by the wave function expansion method. The differential equations governing the vibration of a curved beam is discretized by the General Differential Quadrature (GDQ) method. The domain decomposition method is used to establish the global discrete dynamic equations for the piecewise tunnel. The surrounding soil and the tunnel are coupled together via the stress and the displacement continuation conditions which are implemented by the boundary collocation method. Numerical results demonstrate that the stiffness difference between the liner piece and the connecting joints has a considerable influence on the internal forces of the liner piece.     

A plane strain soil-structure interaction model

A plane strain model for dynamic soil-structure interaction problems under harmonic state is presented. The boundary element method is used to study the response of a homogeneous isotropic linear elastic soil. The far field displacement at the free surface is approximated by an outgoing Rayleigh wave. The finite element method is used to describe the response of the building, of the foundation and possibly of a finite part of the inhomogeneous non-linear soil. Two coupling procedures are described. The model is applied to a problem previously studied in the antiplane case. Incident P, SV and Rayleigh waves are considered. The results show an amplification and an attenuation of the structure motion with frequency when incident Rayleigh waves and P, SV body waves are respectively considered.     

Steady-state harmonic response of a rigid plate bearing on a liquid-saturated poroelastic halfspace

Soil–structure interaction problems are typically modelled by assuming subgrade behaviour to be either elastic or viscoelastic. Herein, compliance functions that may be used to solve soil–structure interaction problems are evaluated by treating the subgrade as a liquid-saturated poroelastic material whose behaviour is governed by Biot's theory. The compliances are evaluated for the harmonic rocking and vertical motions of rigid permeable and impermeable plates bearing on a poroelastic halfspace. Comparisons are made with elastic solutions which assume the subgrade to be either completely drained or undrained. Also, solid and fluid contact stresses are reported for the poroelastic case and compared to the solid contact stresses for the elastic cases.     

A seismic free field input model for FE-SBFE coupling in time domain

A seismic free field input formulation of the coupling procedure of the finite element (FE) and the scaled boundary finite-element (SBFE) is proposed to perform the unbounded soil-structure interaction analysis in time domain. Based on the substructure technique, seismic excitation of the soil-structure system is represented by the free-field motion of an elastic half-space. To reduce the computational effort, the acceleration unit-impulse response function of the unbounded soil is decomposed into two functions; linear and residual. The latter converges to zero and can be truncated as required. With the prescribed tolerance parameter, the balance between accuracy and efficiency of the procedure can be controlled. The validity of the model is verified by the scattering analysis of a hemi-spherical canyon subjected to plane harmonic P, SV and SH wave incidence. Numerical results show that the new procedure is very efficient for seismic problems within a normal range of frequency. The coupling procedure presented herein can be applied to linear and nonlinear earthquake response analysis of practical structures which are built on unbounded soil. Supproted by: the National Key Basic Research and Development Program under Grant No. 2002CB412709     

Uniform half-plane elastodynamic problems by an approximate boundary element method

In this paper the applicability of an approximate Boundary Element Method to uniform half-plane elastodynamic problems is investigated. This method employs the concept of images to construct approximate fundamental solutions for the half-plane and does not require any half-plane surface discretization. The method is formulated in the frequency domain for the case of harmonic disturbances or the Laplace transform domain for the case of transient disturbances. In the latter case a numerical inversion of the transformed solution is necessary to obtain the time domain response. The proposed method can be used as an alternative to boundary element methods that either utilize the infinite plane fundamental solution and thus require a half-plane surface discretization, or employ the exact half-plane fundamental solution, which even though leads to no surface discretization, is of a very lengthy and complicated form. Two characteristics numerical examples are used to illustrate the proposed method and study its advantages and disadvantages.     

A discontinuous Galerkin method for seismic wave propagation in coupled elastic and poroelastic media

Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.     

3D scattering of obliquely incident plane SV waves by an alluvial valley embedded in a fluid-saturated,poroelastic layered half-space

The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obliquely incident SV waves. A wave-number transform is first applied along the valley’s axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is confirmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are discussed in detail resulting in some conclusions.     

Three‐dimensional models of reservoir sediment and effects on the seismic response of arch dams

The important effects of bottom sediments on the seismic response of arch dams are studied in this paper. To do so, a three‐dimensional boundary element model is used. It includes the water reservoir as a compressible fluid, the dam and unbounded foundation rock as viscoelastic solids, and the bottom sediment as a two‐phase poroelastic domain with dynamic behaviour described by Biot's equations. Dynamic interaction among all those regions, local topography and travelling wave effects are taken into account. The results obtained show the important influence of sediment compressibility and permeability on the seismic response. The former is associated with a general change of the system response whereas the permeability has a significant influence on damping at resonance peaks. The analysis is carried out in the frequency domain considering time harmonic excitation due to P and S plane waves. The time‐domain results obtained by using the Fourier transform for a given earthquake accelerogram are also shown. The possibility of using simplified models to represent the bottom sediment effects is discussed in the paper. Two alternative models for porous sediment are tested. Simplified models are shown to be able to reproduce the effects of porous sediments except for very high permeability values. Copyright © 2004 John Wiley & Sons, Ltd.     

A model for the seismic analysis of a periodic viaduct when considering the pile–soil–structure interaction

In this study, a new model is developed for the aseismic design of a periodic viaduct when the pile–soil–structure interaction is considered. To account for the influence of the pile–soil–structure interaction, a wavenumber domain boundary element method (WDBEM) model for the periodic pile row supporting the viaduct is developed using the sequence Fourier transform as well as the boundary element method for the elastic medium. By using the WDBEM model for the pile row, the transfer matrices for the beams and piers, the joint conditions at the beam–beam–pier (BBP) junction as well as the periodicity condition for the viaduct, the wavenumber domain response of the periodic viaduct to spatially harmonic waves is determined. Based on the wavenumber domain response of the viaduct, the space-domain response of the viaduct to an arbitrary seismic wave can be obtained by invoking the inverse sequence Fourier transform method. Numerical results show that when the periodic viaduct is exposed to the spatially harmonic wave, resonances may occur at the bounding frequencies of the passbands of the characteristic waves of the viaduct. Also, it is found that the coincidence between the traveling seismic wave and characteristic waves of the viaduct will generate additional resonant frequencies located in passbands of the characteristic waves.     

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.     

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 effects of moving loads on road pavements: A review

This review paper deals with the dynamic response of road pavements to moving loads on their surface. The road pavement can be modeled as a beam, a plate, or the top layer of a layered soil medium. The foundation soil can be modeled as a system of elastic springs and dashpots or a homogeneous or layered half-space. The material behavior of the pavement can be elastic or viscoelastic, while that of the foundation layers elastic, viscoelastic, water-saturated poroelastic or even inelastic. The loads are concentrated or distributed of finite extent, may vary with time and move with constant or variable speed. The analysis is done by analytical, analytical/numerical and purely numerical methods, such as finite element and boundary element methods, under conditions of plane strain or full three-dimensionality. A number of representative examples is presented in order to illustrate the problem and the methods of analysis, demonstrate the dynamic effects of moving loads on the layered soil medium and indicate the implications of the results on road and airport pavement design.     

Combined paraxial-consistent boundary conditions finite element model for simulating wave propagation in elastic half-space media

In this paper, a finite element model of a soil island is coupled to both a consistent transmitting boundary and a paraxial boundary, which are then used to model the propagation of waves in semi-infinite elastic layered media. The formulation is carried out in the frequency domain while assuming plane strain conditions. It is known that a discrete model of this type, while providing excellent results for a wide range of physical parameters in the context of a half-space problem, may deteriorate rapidly at low frequencies of excitation. This is so because at low frequencies the various waves in the model eventually attain characteristic wavelengths which exceed the distance of the bottom boundary, which then causes that boundary to fail. Also, the paraxial boundaries themselves break down at very low frequencies. In this paper, this difficulty is overcome and the model׳s performance is improved upon dramatically by incorporating an artificial buffer layer sandwiched between the bottom of the soil medium and the underlying elastic half-space. Applications dealing with rigid foundations resting on homogenous or layered half-space media are shown to exhibit significant improvement. Following extensive simulations, clear guidelines are provided on the performance of the coupled model and an interpretation is given on the engineering significance of the findings. Finally, clear recommendations are provided for the practical use of the proposed modelling strategy.     

流体饱和半空间二维地形三分量弹性波散射间接边界元模拟

无限长局部地形地震波斜入射响应问题称为二维三分量问题,在计算量远小于三维的情况下,一定程度上反映了近地表场地的三维动力响应特征.基于天然土体的成层性及固液两相耦合特性,以层状多孔介质内部移动线荷载（孔隙水压）动力格林函数作为基本解,开展流体饱和半空间二维地形三分量弹性波散射的2.5维间接边界元模拟研究.总场响应由自由波场和散射波场叠加构成,前者可由直接刚度法求得,后者则通过施加移动虚拟均布荷载和移动虚拟孔隙水压所产生的动力响应来模拟.该方法优势在于离散仅限于地形底边界（无须离散自由地表）,格林函数计算不存在奇异性（荷载可直接加在边界上）,容易控制计算精度,对复杂边界条件具有很强的适应性.在退化验证和精度比较的基础上,以梯形凹陷和半椭圆沉积地形为例,模拟了时域和频域的流体饱和半空间三维弹性波散射响应.研究表明：局部地形的地震动响应依赖于入射频率、入射角度、边界透水条件、土层刚度和土层厚度等,入射波、反射波和散射波相互干涉,极大延长了位移的振动持续时间.