Consolidation of a poroelastic half-space with anisotropic permeability and compressible constituents by axisymmetric surface loading

The fully coupled Biot quasi-static theory of linear poroelasticity is used to study the consolidation of a poroelastic half-space caused by axisymmetric surface loads. The fluid and solid constituents of the poroelastic medium are compressible and its permeability in the vertical direction is different from its permeability in the horizontal direction. An analytical solution of the governing equations is obtained by taking the displacements and the pore pressure as the basic state variables and using a combination of the Laplace and Hankel transforms. The problem of an axisymmetric normal load is discussed in detail. An explicit analytical solution is obtained for normal disc loading. Detailed numerical computations reveal that the anisotropy in permeability as well as the compressibilities of the fluid and solid constituents of the poroelastic medium have significant effects on the consolidation of the half-space. The anisotropy in permeability may accelerate the consolidation process and may lead to a dilution in the theoretical prediction of the Mandel-Cryer effect. The compressibility of the solid constituents may also accelerate the consolidation process.     

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.     

Dynamic response of pavements on poroelastic half-space soil medium to a moving traffic load

This paper presents an investigation of the steady-state response of pavement systems subjected to a moving traffic load. The traffic loads are simulated by four rectangular load pressures, and the rigid and flexible pavement systems are regarded as an infinite plate resting on a poroelastic half-space soil medium. The contact surface between the plate and the poroelastic half-space is assumed to be smooth and fully permeable. Kirchhoff small-deflection thin-plate theory is employed to analyze the plate, while Biot’s fully dynamic poroelastic theory is used to characterize the poroelastic half-space. The frequency wave-number domain solution of the pavement system is obtained by the compatibility condition between the plate and the poroelastic half-space. By applying the inverse fast Fourier transform, the time domain solution is obtained. Also, the influences of the load speed, the permeability of the soil, and the flexural rigidity of the plate on the response of the pavement system are investigated. The numerical results show that the influences of these parameters on the dynamic response of the pavement system are significant.     

简谐线源荷载作用下热流固耦合地基的动力响应

基于Biot波动理论和广义热弹性理论,对简谐线源荷载（力荷载和热荷载）作用下的热-流-固耦合地基的动力响应问题进行了研究。将地基看成是均质各向同性、完全饱和的多孔半空间介质,利用无量纲化和Fourier变换方法对热-流-固耦合控制方程进行简化,得到了变换域内应力分量、位移分量、温度分布及超孔隙水压力的一般解,并利用Fourier逆变换得到了相应的积分形式解答。通过数值计算对按热-流-固耦合理论、饱和多孔弹性理论和热弹性理论得到的结果进行了比较,同时分析了热-流-固耦合条件下热荷载激振频率对竖向应力、竖向位移、温度分布以及超孔隙水压力的影响。     

空沟对列车运行引起的地基振动隔振效果研究

为研究空沟对高速列车引起地基振动的隔振效果,在已有研究基础上改进出一种新的解析研究模型。模型中首次运用了饱和半空间模型来研究地基上隔振沟对高速列车的隔振效果;隔振沟则通过在饱和半空间土体上设置3个合适宽度、截面为矩形的弹性层来模拟;中间矩形弹性层为路堤,路堤上方放置了枕木与轨道。枕木与轨道分别通过纵向异性Kirchhoff薄板与Euler梁来模拟,饱和土地基采用Biot多孔饱和介质理论来描述。控制方程通过傅立叶变换与傅立叶级数展开,在变换域中进行求解。研究表明,随着列车运行速度的提高,空沟的隔振效果明显提高;饱和土体固-液相的耦合作用对隔振沟的隔振效果的影响明显,尤其当列车运行速度超过土体表面Rayleigh波速时,随着土体渗透系数的增加,空沟隔振效果显著降低。此外,列车运行速度超过土体表面波速时,饱和土地基上空沟的隔振效果明显优于相应单相弹性地基上空沟的隔振效果。     

Thermomechanical response of a poroelastic half-space soil medium subjected to time harmonic loads

The thermomechanical responses of a porous elastic medium subjected to time harmonic loads (normal force and thermal source) are investigated analytically in the context of generalised thermoelastic theory with one relaxation time. The material of the foundation, obeying Biot’s dynamic poroelastic theory, is idealised as a uniform, fully saturated poroelastic half-space stratum. The coupled governing equations are established based on Biot’s dynamic poroelastic theory and on generalised thermoelastic theory. Assuming the disturbances to be harmonically time dependent, the general solutions of stress, displacement, temperature distribution and excess pore water pressure are deduced using the Fourier transform, and the transformed solutions are numerically inverted. The differences among the coupled thermo-hydro-mechanical dynamic model (THMD), the hydro-mechanical dynamic model (HMD) and the thermo-elastic dynamic model (TMD) are discussed. In addition, the effects of the thermal loading frequency on the displacement, stress, temperature distribution and excess pore water pressure components are analysed in the numerical results.     

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.     

Quasi-static deformation due to two-dimensional seismic sources embedded in an elastic half-space in welded contact with a poroelastic half-space

The Biot linearized theory of fluid saturated porous materials is used to study the plane strain deformation of a two-phase medium consisting of a homogeneous, isotropic, poroelastic half-space in welded contact with a homogeneous, isotropic, perfectly elastic half-space caused by a two-dimensional source in the elastic half-space. The integral expressions for the displacements and stresses in the two half-spaces in welded contact are obtained from the corresponding expressions for an unbounded elastic medium by applying suitable boundary conditions at the interface. The case of a long dip-slip fault is discussed in detail. The integrals for this source are solved analytically for two limiting cases: (i) undrained conditions in the high frequency limit, and (ii) steady state drained conditions as the frequency approaches zero. It has been verified that the solution for the drained case (ω → 0) coincides with the known elastic solution. The drained and undrained displacements and stresses are compared graphically. Diffusion of the pore pressure with time is also studied.     

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

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

横观各向同性饱和地基三维瞬态Lamb问题

研究了横观各向同性饱和土地基在地表动力荷载作用下的三维瞬态响应。基于饱和多孔介质的三维Biot波动理论,利用Laplace变换,建立圆柱坐标系下横观各向同性饱和土的波动方程;解耦波动方程后,根据算子理论,并借助Fourier展开和Hankel变换技术,得到瞬态荷载作用下,饱和土介质的土骨架位移和应力、孔隙水相对位移和孔隙水压力的一般解;利用一般解,给出横观各向同性饱和地基在地表集中荷载激励下的瞬态Lamb问题的解答。数值算例结果表明,采用各向同性饱和介质的动力学模型,不能准确描述具有明显各向异性特性的饱和土地基的瞬态动力特性。     

层状地基中单桩负摩擦问题积分方程解法

利用Biot固结理论和积分方程方法研究了表面有堆载的层状地基中单桩负摩擦问题。根据层状饱和土的圆形载荷基本解得出了单桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化。再结合传递和刚度矩阵传递到各个层中去,对变换域内的积分方程采用Schapery 逆变换方法得到时域内单桩的近似积分方程。求解积分方程组并进行相应的数值逆变换,就可得出层状地基中的单桩在表面圆形均布载荷作用下的位移、轴力、孔压和桩侧摩阻力随时间的变化情况。计算结果表明,桩侧剪力和孔压分层明显。     

群桩基础引发饱和地基振动的近似解析解

建立了饱和半空间-群桩基础耦合动力近似解析模型,其中基桩考虑为欧拉梁,饱和地基采用Biot两相介质动力控制方程,二者在频率-波数域中利用桩-土相互作用点处的饱和地基柔度矩阵进行耦合,并采用快速傅里叶逆变换获得频域解答。分析了饱和地基中群桩基础的动力阻抗及刚性承台上施加简谐振动荷载时群桩基础引起的饱和地基振动响应。研究结果表明,荷载类型、激振频率以及地基渗透系数对饱和地基位移和孔压响应有明显影响。特别的,群桩基础动力阻抗峰值频率会随着地基渗透系数的增加而增大。     

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.     

层状饱和土体中排桩对简谐荷载隔振效果分析

根据Biot理论,利用已有的传递、透射矩阵法得到层状饱和土体内部受竖向圆形分布荷载作用下的基本解,再由Muki虚拟桩法,建立了频域内层状土-桩的第二类Fredholm积分方程。通过离散方法求解积分方程,得到了评价隔振效果的振幅比。与已知文献结果比较,验证了方法的正确性。数值结果表明：对于同种类型的振源,采用相同的隔振系统,在上软下硬的多层土体中的隔振效果比在上硬下软的土层中要好。采用较长的桩或刚度较大的桩,或桩之间的距离加密都可得到好的隔振效果。     

考虑热-水-力耦合效应多孔弹性地基的动力响应

通过对Biot波动方程的修正,得到考虑热-水-力学耦合效应的多孔弹性介质动力响应的控制方程,研究了简谐均布荷载作用下地基土体的热-水-力耦合动力响应问题。利用Fourier变换技术,得到地基中的应力、位移和孔隙水压力积分形式的解答。利用Fourier逆变换得到数值结果,分析了热-水-力学耦合条件下地基土体中温度增量、应力、位移和孔隙水压力响应的分布,并讨论了热源输入的影响, 结果表明：应力、位移和孔隙水压力随 的增大而有一定的减小。     

淤砂层模拟方法对上层理想流体动力反应的影响——平面P波入射

采用解析法,以弹性半空间-淤砂层-理想流体层为对象,分别将淤砂层模拟为固体介质和流体介质,分析了平面P波入射时上层理想流体的动力反应特性。通过算例分析,比较了将淤砂层模拟为弹性固体、黏弹性固体、理想流体和黏性流体介质时与将淤砂层模拟为两相多孔介质时计算结果的差异。研究表明当渗透系数很小时,多孔介质可以近似简化为黏弹性固体,甚至线弹性固体。这一分析研究对于在高坝-库水-淤砂-地基系统地震反应分析的数值模型中采用合理的淤砂层模型,提高计算效率具有重要参考意义,而且还可以用于校核数值模型的合理性。     

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.     

Transient dynamic response of multilayered saturated media subjected to impulsive loadings

This paper presents a stable and efficient method for calculating the transient solution of layered saturated media subjected to impulsive loadings by means of the analytical layer element method. Starting with the field equations based on Biot's linear theory for porous, fluid‐saturated media, and the seepage continuity equation, an analytical layer element for a single layer is established by applying Laplace‐Hankel integral transform. The global stiffness matrix in the transform domain for a layered saturated half‐space subjected to a transient circular patch loading is obtained by assembling the layer elements of each layer. The displacements in the time domain are derived by Laplace‐Hankel inverse transform of the global stiffness matrix. Numerical examples are conducted to verify the accuracy of the method and to demonstrate the influences of type of transient loading, buried depth of loading, permeability, and stratification of materials on the transient response of the multilayered saturated poroelastic media.     

Poroelastic behaviour of saturated brittle rock with anisotropic damage

A new anisotropic poroelastic damage model is proposed for saturated brittle porous materials. The model is formulated in the framework of the continuum damage mechanics. A second‐rank symmetric tensor is used to characterize material damage due to oriented microcracks. The classic Biot poroelastic theory is then extended to include poroelastic damage coupling. Both the deterioration of elastic properties and poroelastic coefficients is taken into account. A suitable procedure for determination of model parameters from standard laboratory tests is presented. The validity of the model is tested through comparison between numerical predictions and experimental data in various loading conditions. The overall performance of the model is evaluated. The choice of relevant effective stress for the microcrack propagation criterion in saturated cohesive geomaterials is discussed. Copyright © 2000 John Wiley & Sons, Ltd.     

Dynamics of unsaturated poroelastic solids at finite strain

We derive the governing equations for the dynamic response of unsaturated poroelastic solids at finite strain. We obtain simplified governing equations from the complete coupled formulation by neglecting the material time derivative of the relative velocities and the advection terms of the pore fluids relative to the solid skeleton, leading to a so‐called u^s ? p^w ? p^a formulation. We impose the weak forms of the momentum and mass balance equations at the current configuration and implement the framework numerically using a mixed finite element formulation. We verify the proposed method through comparison with analytical solutions and experiments of quasi‐static processes. We use a neo‐Hookean hyperelastic constitutive model for the solid matrix and demonstrate, through numerical examples, the impact of large deformation on the dynamic response of unsaturated poroelastic solids under a variety of loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.