A Damage Mechanics Model for Porous Media

When considering saturated poroelastic geomaterials, their consolidation response can be influenced by the evolution of damage in the porous skeleton. This paper introduces the concept of effective stress and the modified Mohr–Coulomb criterion as a means for examining the problem of consolidation response of a damage-susceptible poroelastic geomaterial, which is considered to be a novel development in the application of damage mechanics to the study of poroelastic phenomena. Then, a damage mechanics model for porous media is established and applied to analyze the Biot’s consolidation problem of damaged soils. The comparison between the numerical predictions and the classical solutions shows that the computational damage model developed in this paper is effective and feasible in analyzing the consolidation problem of damaged porous media.     

Dynamic responses of a poroelastic half‐space from moving trains caused by vertical track irregularities

The vibrations of railway tracks on a poroelastic half‐space generated by moving trains are investigated through a vehicle–track–ground coupling model. The theoretical model incorporates a vehicle, a track, and a fully saturated poroelastic half‐space soil medium. The source of vibration excitation is divided into two components: the quasi‐static loads and the dynamic loads. The quasi‐static loads are related to the static component of the axle loads, whereas the dynamic loads are due to the dynamic wheel–rail interaction. A linear Hertizian contact spring is introduced between each wheelset and the rail to consider the dynamic loads. Biot's dynamic theory is used to characterize the poroelastic half‐space soil medium. Using the Fourier transform, the governing equations for the track–ground system are solved and the numerical results are presented for a single axle vehicle model. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium are investigated. In addition, the different roles of the moving axle loads and the roughness‐induced dynamic loads are identified. It is concluded that the vibration level of the free field off the track predicted by the poroelastic soil medium is smaller than that predicted by the elastic soil medium for vehicle speed below the Rayleigh wave speed of the poroelastic half‐space, whereas it is larger for vehicle speed above the Rayleigh wave speed. The dynamic loads play an important role in the dynamic responses of the track–ground system. Copyright © 2010 John Wiley & Sons, Ltd.     

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

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

爆轰荷载作用下球空腔热流固耦合动力响应

基于饱和多孔弹性介质热流固耦合动力模型(THMD)及控制方程,研究了爆轰荷载作用下球空腔内壁面受到随时间变化热、力冲击作用下的热流固耦合动力响应。利用Laplace变换技术,对控制方程进行解耦,获得了温度、位移、孔隙水压力和应力积分形式解。采用Laplace逆变换得到数值结果,分析了热流固耦合参数的影响特性,并将热流固耦合结果与热弹性条件下的结果进行比较,以验证热流固耦合条件下结果的正确性。     

横观各向同性饱和地基竖向振动的衰减特性

基于提出的横观各向同性饱和多孔介质Biot波动方程的一般解,研究了饱和半空间地基在竖向点源简谐激振荷载作用下地表振动的衰减特征,分析了激振频率以及横观各向同性饱和土介质的各向异性参数和孔隙渗透系数对地表振动特征的影响。计算结果表明,低频和高频激振时,地表位移衰减特性存在明显差异;在饱和土的各向异性参数中,纵向和水平方向动态渗透系数比值和刚度系数比值对地表位移衰减影响最大,这也说明采用各向同性饱和介质的动力学模型不能准确地描述具有明显各向异性特性的饱和土地基的动力特性。     

Mode of a spherical cavity’s thermo-elastodynamic response in a saturated porous medium for non-torsional loads

The thermo-hydro-elastodynamic model (THED), in which the thermo-osmosis and thermal-filtration phenomenon in a two-phase porous thermoelastic medium can be considered, was previously presented by the authors [Ganbin Liu, Kanghe Xie, Rongyue Zheng. Model of nonlinear coupled thermo-hydro-elastodynamics response for a saturated poroelastic medium. Sci China (Ser E) 2009;52(8):2373–83] and is used in this paper to investigate the thermo-elastodynamic response of a spherical cavity in a saturated poroelastic medium when subjected to a time-dependent non-torsional thermal/mechanical source. The Separated Variable Method is introduced, and the non-zero displacement potentials are expanded in terms of the Legendre polynomial. Solutions for the displacement, temperature increment, pore pressure and stress are obtained in the domain of the Laplace transform. Numerical results are also performed for different modes and are compared with the results of the thermoelastic model (TED) to ascertain the validity and the difference between these two models.     

A 2.5-D coupled FE–BE model for the dynamic interaction between saturated soil and longitudinally invariant structures

A two-and-a-half-dimensional (2.5-D) coupled finite element–boundary element (FE–BE) model is presented to simulate the three-dimensional dynamic interaction between saturated soils and structures with longitudinally invariant geometries. A regularized 2.5-D boundary integral equation for saturated porous media is derived that avoids the evaluation of singular traction integrals. The 2.5-D coupled FE–BE model is established by using the continuity conditions on the soil–structure interface. The developed model is verified through comparison with an existing semi-analytical method. Two case studies of a tunnel embedded in a poroelastic half-space and the efficiency of a vibration isolating screen are presented.     

Wave propagation in thermoelastic saturated porous medium

Biot’s theory for wave propagation in saturated porous solid is modified to study the propagation of thermoelastic waves in poroelastic medium. Propagation of plane harmonic waves is considered in isotropic poroelastic medium. Relations are derived among the wave-induced temperature in the medium and the displacements of fluid and solid particles. Christoffel equations obtained are modified with the thermal as well as thermoelastic coupling parameters. These equations explain the existence and propagation of four waves in the medium. Three of the waves are attenuating longitudinal waves and one is a non-attenuating transverse wave. Thermal properties of the medium have no effect on the transverse wave. The velocities and attenuation of the longitudinal waves are computed for a numerical model of liquid-saturated sandstone. Their variations with thermal as well as poroelastic parameters are exhibited through numerical examples.     

衬砌隧道中移动轴向激励作用下两相多孔介质动力响应

为了研究衬砌隧道中地铁振动在饱和地层中的传播情况,采用移动轴向激励模拟技术,建立了隧道-衬砌-两相介质的动力分析模型;利用波函数展开法、傅里叶变换法等,推导了频域内衬砌隧道移动轴向激励作用下两相多孔介质动力响应的解析解,并给出了两相多孔介质临界速度的经验公式;通过离散快速傅里叶逆变换得到了时-空域内两相介质的动力响应。结果表明：无衬砌隧道中移动的轴向常激励作用下,两相介质临界速度只与介质的剪切模量和密度有关,且数值接近其剪切波速的1.1倍;对于衬砌隧道,介质的临界速度随着衬砌剪切模量的增大而增大,随着衬砌密度的增大而减小,衬砌对振动向介质中的传播有一定的削弱作用,衬砌剪切模量与介质剪切模量相差越大,削弱越明显;动力响应频率越接近激振频率,其幅值越大,所对应的临界速度越小。     

Dynamic response of a non-circular lined tunnel with visco-elastic imperfect interface in the saturated poroelastic medium

A semi-analytical method is proposed to investigate the dynamic response of a non-circular tunnel with visco-elastic imperfect interface in poroelastic medium. Biot’s dynamic theory is used to simulate the saturated poroelastic medium, and the governing equations are solved by reducing them into three Helmholtz equations that the potential functions satisfy. The visco-elastic interface model with elastic and viscosity coefficients is adopted to analyze the interface effect around the non-circular lined tunnel. The analytic solutions of displacements and stresses are expanded in terms of wave functions. Some numerical examples are provided to analyze the effect of visco-elastic interface on the dynamic stress around the tunnel.     

The fundamental solution of poroelastic plate saturated by fluid and its applications

In this paper, the numerical model of the transverse vibrations of a thin poroelastic plate saturated by a fluid was proposed. Two coupled dynamic equations of equilibrium related to the plate deflection and the equivalent moment were established for an isotropic porous medium with uniform porosity. The fundamental solutions for a porous plate were derived both in the Laplace transform domain and in the time domain. A meshless method was developed and demonstrated in the Laplace transform domain for solving two coupled dynamic equations. Numerical examples demonstrated the accuracy of the method of the fundamental solutions and comparisons were made with analytical solutions. The proposed meshless method was shown to be simple to implement and gave satisfactory results for a poroelastic plate dynamic analysis. Copyright © 2009 John Wiley & Sons, Ltd.     

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

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

An anisotropic damage–plasticity model for saturated quasi‐brittle materials

This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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

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

Numerical analysis of 1‐D compression wave propagation in saturated poroelastic media

A numerical procedure is presented for the simulation of 1‐D compression wave propagation in saturated poroelastic media. The media are modelled as a two‐phase system consisting of compressible fluid and solids. Viscous coupling forces resulting from the relative motion between phases are characterized as Darcy type. The numerical procedure can account for effects of axial strain, nonlinear material behaviour, and various drained and undrained boundary conditions. Time integration is carried out explicitly and isothermal conditions are assumed. The method is capable of modelling shock wave fronts without introducing artificial viscosity. Numerical results are in close agreement with analytical solutions for several simplified cases and indicate that mass coupling may have important effects on fluid velocity and wave speed. Corresponding effects on solid velocity and wave speed are much smaller. Numerical results also indicate that damping occurs in a saturated poroelastic column and is dependent on the value of hydraulic conductivity. Copyright © 2007 John Wiley & Sons, Ltd.     

层状饱和地基上刚性圆板的扭转振动及参数分析

用解析的方法首次研究了层状饱和地基上刚性圆板的扭转振动特性.运用Hankel变换求解饱和介质动力问题的控制方程,按混合边值条件建立了层状饱和地基上刚性圆板扭转振动的第二类Fredholm积分方程.数值算例给出了动力柔度系数和扭转角幅值随无量纲频率的变化曲线,与单相弹性及匀质饱和介质情况进行了对比分析,并进行了参数研究.     

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

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

列车荷载下上覆弹性层横观各向同性饱和地基的动力响应

基于Biot多孔弹性介质的波动理论,构造了轨道、道碴、枕木及弹性层的横观各向同性饱和地基在列车荷载下的动力计算模型。利用Fourier变换,得到了列车荷载作用下横观各向同性软土地基上弹性层动力响应的解析结果。利用离散Fourier逆变换得到数值计算结果,分析了荷载速度、地基的各向异性参数、弹性层刚度系数及厚度对位移和孔压响应的影响。分析结果表明：弹性层对控制地基振动作用显著,地表振动幅值随荷载速度的增加而增大,软土的横向弹性模量对地表振动及土中孔压有较大影响。     

Three‐dimensional time‐harmonic fundamental solutions for a fluid‐saturated poroelastic half‐space with partially permeable free surface

By virtue of a pair of scalar potentials for the displacement of the solid skeleton and the pore fluid pressure field of a saturated poroelastic medium, an alternative solution method to the Helmholtz decomposition is developed for the wave propagation problems in the framework of Biot's theory. As an application, a comprehensive solution for three‐dimensional response of an isotropic poroelastic half‐space with a partially permeable hydraulic free surface under an arbitrarily distributed time‐harmonic internal force field and fluid sources is developed. The Green's functions for the poroelastic fields, corresponding to point, ring, and disk loads, are reduced to semi‐infinite complex‐valued integrals that can be evaluated numerically by an appropriate quadrature scheme. Analytical and numerical comparisons are made with existing elastic and poroelastic solutions to illustrate the quality and features of the solution. Copyright © 2016 John Wiley & Sons, Ltd.