横观各向同性饱和地基中埋置荷载的非轴对称瞬态响应

基于孔隙介质的Biot理论,首先利用Laplace变换,给出圆柱坐标系下横观各向同性饱和弹性多孔介质在变换域上的波动方程;将波动方程解耦后,根据方位角的Fourier展开和径向Hankel变换,求解了Biot波动方程,得到以土骨架位移、孔隙水压力和土介质总应力分量的积分形式的一般解;借助一般解,建立了有限厚度饱和土层和饱和半空间的精确动力刚度矩阵,并由土层的层间界面连续条件建立三维非轴对称层状饱和地基的总刚度方程;在此基础上,系统研究了横观各向同性饱和半空间体在内部集中荷载激励下的动力响应,并给出了问题的瞬态解答.该研究为运用边界元法求解饱和地基动力响应奠定了理论基础.     

A practical numerical method for large strain liquefaction analysis of saturated soils

Accurate prediction of the liquefaction of saturated soils is based on strong coupling between the pore fluid phase and soil skeleton. A practical numerical method for large strain dynamic analysis of saturated soils is presented. The u–p formulation is used for the governing equations that describe the coupled problem in terms of soil skeleton displacement and excess pore pressure. A mixed finite element and finite difference scheme related to large strain analysis of saturated soils based on the updated Lagrangian method is given. The equilibrium equation of fluid-saturated soils is spatially discretized by the finite element method, whereas terms associated with excess pore pressure in the continuity equation are spatially discretized by the finite difference method. An effective cyclic elasto-plastic constitutive model is adopted to simulate the non-linear behavior of saturated soils under dynamic loading. Several numerical examples that include a saturated soil column and caisson-type quay wall are presented to verify the accuracy of the method and its usefulness and applicability to solutions of large strain liquefaction analysis of saturated soils in practical problems.     

Propagation of seismic waves through liquefied soils

To predict the earthquake response of saturated porous media it is essential to correctly simulate the generation, redistribution, and dissipation of excess pore water pressure during and after earthquake shaking. To this end, a reliable numerical tool requires a dynamic, fully coupled formulation for solid–fluid interaction and a versatile constitutive model. Presented in this paper is a 3D finite element framework that has been developed and utilized for this purpose. The framework employs fully coupled dynamic field equations with a u–p–U formulation for simulation of pore fluid and solid skeleton interaction and a SANISAND constitutive model for response of solid skeleton. After a detailed verification and validation of the formulation and implementation of the developed numerical tool, it is employed in the seismic response of saturated porous media. The study includes examination of the mechanism of propagation of the earthquake-induced shear waves and liquefaction phenomenon in uniform and layered profiles of saturated sand deposits.     

二维土坝的两相动力分析

提出了一种分析饱和土坝动力反应的方法,考虑了土坝的两相介质特性,在固液耦联动力方程的基础上,选取固相位移,液相位移、孔隙水压作为场变量,采用伽辽金加权残数法进行有限元空间离散化,然后在时域上采用Ｗｉｌｓｏｎ－θ法进行逐步积分。该方法不仅能计算出固相位移和液相位移,而且能直接得到孔隙水压的反应过程。文中以一饱和土坝模型进行算例分析,并与将其作为单相介质时的结果进行了比较。该法可用于分析饱和介质的地震     

Three‐dimensional dynamic response of underground openings in saturated rock masses

In this paper a three‐dimensional dynamic analysis program for saturated porous rocks and soils is presented. The theoretical formulations incorporated in the proposed computer program are the extension of Biot's two‐phase theory to non‐linear region. Numerical study using the program is carried out to demonstrate the three‐dimensional features of the program and to show the effects of pore water on the dynamic response of underground openings in saturated rock masses. It is shown that underground openings in saturated porous media could be significantly more vulnerable to the potential damages associated with high motions and shear failure than those in dry media. Copyright © 2001 John Wiley & Sons, Ltd.     

Decoupling of the coupled poroelastic equations for quasistatic flow in deformable porous media containing two immiscible fluids

The study of the poroelastic behavior of sedimentary materials containing two immiscible fluids in response to either applied stress or pore pressure change in a quasistatic limit, i.e., negligible second time-derivatives, is of great importance to many hydrogelogical problems, e.g., land subsidence caused by withdrawal of subsurface fluids. The poroelasticity models developed for analyzing these problems feature partial differential equations that are coupled in the terms describing viscous damping and strain field. To determine closed-form analytical solutions for induced volumetric strain (dilatation) of the solid framework and its interaction with fluid flows, the choice of normal coordinates whose transformation can be performed to decouple these poroelastic equations is highly desirable. In this paper, we show that normal coordinates for decoupling these equations are real-valued and equal to three different linear combinations of the dilatations of the solid and the fluids (or equivalently, three different linear combinations of two individual fluid pressures and solid dilatation). In contrast to fully saturated porous media, it is found that the viscous damping effect must be represented in normal coordinates in the presence of the second fluid. The resulting decoupled equations representing independent motional modes are a Laplace equation and two diffusion equations, which can be solved analytically under a variety of initial and boundary conditions. Thus, after inverse transformation of normal coordinates is performed, the closed-form analytical solutions for induced solid volumetric strain and excess pore fluid pressures can be obtained simultaneously from our decoupled partial differential equations.     

基于Comsol的饱和多孔介质动力方程的数值模拟及应用

两相饱和多孔介质的动力响应问题在地震工程领域具有重要的研究意义,由于涉及到固相和液相的动力耦合,使得该问题的求解尤为复杂。本文利用Comsol在求解多场耦合问题上的优点,针对Biot饱和多孔介质u-U耦合形式下的波动方程特征,经过一系列微分算子运算和矩阵变换得到导数形式下的波动方程,基于Comsol Multiphysics提供的广义偏微分方程模式对变形后的波动方程进行求解,并把改进后的无限元边界应用到无限域动力问题的模拟中。通过与饱和多孔介质动力响应的解析解进行对比,验证模型求解技术的可行性和正确性,并在此基础上讨论饱和土地基中空沟隔振效果与饱和土体参数孔隙率、泊松比的关系。通过研究分析,可以为饱和土地基中空沟隔振设计提供一些有价值的参考。     

含地下水的土层对地震波传播的影响

本文综合地介绍了含地下水的土层对地震波传播的研究成果,并对一些普遍性的原则问题进行了讨论.在不变形介质的前提下详细研究了二维地下水的波动问题,应用和开发了线单元-特征线法、交替方向隐式差分法和伽辽金有限单元法,求解了一些边值问题.还研究了介质的不均匀性、地层的成层性等对地下水波动的影响以及地下水波在界面上的透射、反射现象,发现了一些新的规律.从普遍方程出发,论证了地下水波在可变形介质和不变形介质中的差别,导出了两者的换算关系.开辟了由不变形介质的解求可变形介质的解的一个间接途径. 在可变形介质的前提下,直接从两相介质动力方程组出发,对含水土层中的波动问题进行了详尽研究,考虑了弹塑性、非线性本构关系、剪胀性、刚度退化等各种因素的影响.通过计算对各种土壤刚度、渗阻条件、边界透水条件、荷载形式、本构关系和不均匀性等情况做了比较研究,发现一些新规律.在垂直地震荷载作用下发现了两种不同形态的孔隙水压反应--阶变形和振荡形的机制,并用间断波理论解释了孔压阶跃的原因.最后对发现的新现象做了一些讨论和总结.     

饱和土体瞬态响应有限元分析

给出基于Biot多孔介质理论分析饱和土体在动载荷作用下瞬态响应的有限元公式,数值计算部分采用本文有限元法分别计算一维饱和土柱在两种不同类型动载荷作用下的瞬态响应,并将数值计算结果与文献中的解析解进行比较,二者结果十分吻合,从而验证本文方法的可行性。     

Low-frequency dilatational wave propagation through fully-saturated poroelastic media

The strong coupling of applied stress and pore fluid pressure, known as poroelasticity, is relevant to a number of applied problems arising in hydrogeology and reservoir engineering. The standard theory of poroelastic behavior in a homogeneous, isotropic, elastic porous medium saturated by a viscous, compressible fluid is due to Biot, who derived a pair of coupled partial differential equations that accurately predict the existence of two independent dilatational (compressional) wave motions, corresponding to in-phase and out-of-phase displacements of the solid and fluid phases, respectively. The Biot equations can be decoupled exactly after Fourier transformation to the frequency domain, but the resulting pair of Helmholtz equations cannot be converted to partial differential equations in the time domain and, therefore, closed-form analytical solutions of these equations in space and time variables cannot be obtained. In this paper we show that the decoupled Helmholtz equations can in fact be transformed to two independent partial differential equations in the time domain if the wave excitation frequency is very small as compared to a critical frequency equal to the kinematic viscosity of the pore fluid divided by the permeability of the porous medium. The partial differential equations found are a propagating wave equation and a dissipative wave equation, for which closed-form solutions are known under a variety of initial and boundary conditions. Numerical calculations indicate that the magnitude of the critical frequency for representative sedimentary materials containing either water or a nonaqueous phase liquid is in the kHz–MHz range, which is generally above the seismic band of frequencies. Therefore, the two partial differential equations obtained should be accurate for modeling elastic wave phenomena in fluid-saturated porous media under typical low-frequency conditions applicable to hydrogeological problems.     

Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

二维非均匀饱和土体的地震反应分析

将饱和土体视为由弹性骨架和不可压缩流体组成的两相孔隙介质,取固相位移,液相位移,孔隙水压力为场变量。从两相孔介质的动力耦合微分方程出发,利用伽辽金原理和Ｗｉｌｓｏｎ－θ法导出了三场有限元积分方法,该方法可进行二维非均匀饱和土体的地震反应分析,文中通过算例讨论了非均匀饱和土体地震反应的一些特征。     

Finite-element simulation of hot-water-type geothermal reservoirs

A fully coupled finite-element model is presented to simulate geothermal reservoirs in terms of the unknown variables of displacement, pore pressure and temperature using an elasto-plastic stress-strain relationship. The mathematical model combines the simultaneous mechanisms of subsidence, mass flow and heat flow through porous media by utilizing the equilibrium, fluid flow and energy flow equations. The governing partial differential equations are based upon constant physical parameters (except fluid density) and formulated for hot-water-type geothermal reservoirs. A simultaneous solution algorithm is applied to the modelling equations yielding transient and spatial distributions of the unknown variables. Application of the simulator to actual field conditions is performed by simulating the production history of a well located in the Denizli-Kızıldere geothermal field in Turkey. The results show an excellent agreement between the computed and the observed values of temperature and pore pressure variations throughout the reservoir.     

Boundary integral formulation and two-dimensional fundamental solutions for dynamic behavior analysis of unsaturated soils

In this paper the coupled equations governing the dynamic behavior of unsaturated soils are derived based on the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) [Gatmiri B. Analysis of fully coupled behavior of unsaturated porous medium under stress, suction and temperature gradient. Final report of CERMES-EDF, 1997] and Gatmiri et al. (1998) [Gatmiri B, Delage P, Cerrolaza M, UDAM: a powerful finite element software for the analysis of unsaturated porous media. Adv Eng Software 1998; 29(1): 29–43]. In this formulation, the solid skeleton displacements, water pressure and air pressure are presumed to be independent variables. The Boundary Integral formulations as well as fundamental solutions for such a dynamic u–p_w–p_a theory are presented in this paper for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated two dimensional (2D) fundamental solutions for such deformable porous medium with linear elastic behavior are derived in Laplace transform domain using the method of Hörmander. Finally, some numerical results are presented to show the accuracy of the proposed solutions. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case.     

饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法

本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.     

基于流固耦合动力模型的饱和土体-隧道体系地震反应研究

基于ABAQUS软件平台,应用自行开发的流固耦合动力模型孔压单元模拟场地土体,并通过黏弹性人工边界方法实现地震动的输入,对饱和土体场地中的双孔隧道结构在地震荷载作用下的动力反应进行研究。计算结果表明:在地震反应结束时刻,场地土体位移幅值在两隧道之间以及两隧道的附近区域较大,而远离隧道的区域则较小;场地底部区域土体的孔压幅值较大,而场地顶部区域土体则较小;隧道左右两侧拱腰部位的衬砌的应力较大,而拱顶部位则较小。计算结果同时表明了流固耦合动力模型孔压单元在饱和土体-隧道体系地震反应研究中的适用性。     

饱和土中浅埋隧道耦合二维动力响应解答探讨

针对平面波耦合作用下饱和半空间中浅埋隧道二维动力响应的典型问题的解答进行探讨。首先基于Biot饱和孔隙介质理论建立饱和土场方程,并通过中间势函数的代入推导得到饱和土体中的快纵波、慢纵波及横波的Helmholtz势函数方程,随后引入柱坐标系下饱和土体和弹性衬砌的势函数通解,并代入边界条件中,建立以第n组参数为变量的多元一次齐次线性方程组。采用本文的推导思路,可以解决饱和土体中的浅埋隧道在地震作用时P波和S波耦合作用下的响应问题。     

基于Kelvin-Voigt黏弹性骨架的含非饱和流体孔隙介质BISQ模型(英文)

本文综合考虑了在波传播过程中孔隙介质的三种重要力学机制——"Biot流动机制一squirt流动机制-固体骨架黏弹性机制",借鉴等效介质思想,将含水饱和度引入波动力学控制方程,并考虑了不同波频率下孔隙流体分布模式对其等效体积模量的影响,给出了能处理含粘滞性非饱和流体孔隙介质中波传播问题的黏弹性Biot/squirt(BISQ)模型。推导了时间-空间域的波动力学方程组,由一组平面谐波解假设,给出频率-波数域黏弹性BISQ模型的相速度和衰减系数表达式。基于数值算例分析了含水饱和度、渗透率与频率对纵波速度和衰减的影响,并结合致密砂岩和碳酸盐岩的实测数据,对非饱和情况下的储层纵波速度进行了外推,碳酸盐岩储层中纵波速度对含气饱和度的敏感性明显低于砂岩储层。     

Time harmonic point load and dynamic contact problem of contacting fluid and poroelastic half-spaces

The dynamic response of contacting fluid and fluid-saturated poroelastic half- spaces to a time-harmonic vertical point force or a point pore pressure is investigated. The solutions are formulated using the boundary conditions at the fluid-porous medium interface. The point load solutions are then used to solve the dynamic problem of the vertical vibration of a rigid disc (both permeable and impermeable discs are included) on the surface of the poroelastic half-space. The contact problems are solved by integrating the point force and point pore pressure solutions over the contact area with unknown discontinuous force and pore pressure distributions, which are determined from the boundary conditions. The solutions are expressed in terms of dual integral equations, which are converted to Fredholm integral equations of the second kind and solved numerically. Selected numerical results for the vertical dynamic compliance coefficient for the cases with or without fluid overlying the poroelastic half-space are presented to show the effects of the fluid. The influence of the permeability condition of the disc on the compliance of the poroelastic half-space is investigated. The displacement, vertical stress, pore pressure in the poroelastic half-space and water pressure in the fluid half-space are also examined for different poroelastic materials and frequencies of excitation. The present results are helpful in the study of the dynamic response of foundations on the seabed under seawater.     

二相介质饱和土中群桩动力阻抗分析

用流体饱和多孔介质材料描述土体，由饱和土和群桩及承台系统的位移协调条件和力平衡条件建立饱和土和群桩及承台系统动力相互作用的控制方程，分析饱和土中群桩动力阻抗。结果表明：孔隙流体对饱和土中桩基础动力阻抗有一定的影响；在饱和土具有不同的流体渗透系数时，饱和土中群桩动力阻抗也有一定差别。在地基上与基础结构动力相互作用研究中应该考虑地基土中孔隙流体的影响。