平面P1波斜入射下海底洞室地震响应解析分析

基于无黏性可压缩理想流体介质波动理论和Biot流体饱和多孔介质波动理论,考虑水下饱和土的流固耦合,借助Hankel函数积分变换法(HFITM)给出入射平面P1波在海底洞室周围散射问题的解析解。相比传统研究中的"大圆弧假定",Hankel函数积分变换法可以较好地处理半空间表面边界条件。利用该解析解,计算分析了洞室表面透水条件、入射角度、入射频率、海水水深和饱和土的孔隙率等因素对水-土交界面处水平位移、竖向位移和洞室表面动水压力、环周总应力的影响。结果表明:洞室表面透水条件对水-土交界面处水平位移和竖向位移影响较小;随着斜入射角的增加,水-土交界面处竖向位移减小;随着入射频率的增加,水-土交界面处水平位移随之增加;海水水深为2.5倍SV波的波长时,水-土交界面处水平位移及洞室表面动水压力最大值最大;随着孔隙率的增加,水-土交界面处水平位移、竖向位移和洞室表面环周总应力减小,而洞室表面动水压力随之增加。     

简谐荷载作用下饱和土体中圆形衬砌隧道三维动力响应分析

研究了全空间饱和土体中圆形衬砌隧道在径向简谐点荷载作用下的三维动力响应,将衬砌用无限长圆柱壳来模拟,土体用Biot饱和多孔介质模型来模拟,引入两类势函数来表示土骨架的位移和孔隙水压力,并利用修正Bessel方程来求解各势函数,结合边界条件,得到频率-波数域内衬砌和土骨架位移、孔隙水压力的解答,最后进行Fourier逆变换得到时间-空间域内的响应。通过算例分析了荷载振动频率和土体渗透性对土体和衬砌位移响应及土体孔压的影响。结果表明,饱和土体和弹性土体的位移响应具有明显区别。随着荷载频率的增大,土体和隧道位移幅值减小,土体孔压幅值增大;随着土体渗透性增大,土体位移及孔压幅值减小。     

频域内饱和多孔介质的参数反演分析

基于Biot理论,假定固体颗粒和孔隙内流体均不可压缩,建立了以固体骨架位移表示的的控制方程.考虑单层饱和多孔介质在竖向简谐荷载作用下一维动力响应,通过理论推导获得了骨架位移、应力以及孔隙流体压力等物理量的解析表达式.基于饱和土的简谐动力模型试验数据,与所得到的理论解答相结合,将饱和多孔介质材料参数反演问题归结为非线性多峰函数的最优化问题.全局最优解的求解采用了遗传算法和模拟退火算法,并通过试验和数值算例验证了所得材料参数的正确性.     

Rayleigh波在饱和半空间中圆形洞室周围的散射

基于Biot的两相介质理论,采用间接边界积分方程法,求解了Rayleigh波在饱和半空间中圆形洞室周围的二维散射问题,分析了入射波频率、孔隙率、边界渗透条件和洞室埋深等参数对地表位移幅值、洞室表面和地表孔隙水压力的影响。研究结果表明,由于圆形洞室的存在,饱和半空间地表位移幅值和地表孔隙水压被明显放大:在透水条件下,水平和竖向位移最大值分别比自由场放大了10.1倍和11.2倍;在不透水条件下,水平和竖向位移最大值分别比自由场放大了12.0倍和9.6倍,地表孔隙水压力放大了2.1~3.0倍。地表位移和孔隙水压力最大值均出现在入射波近端洞室边界附近。随着Rayleigh波入射频率的增大和洞室埋深的增加,地表位移幅值的放大作用有所减小。在相同孔隙率条件下,当入射频率为1.0时,洞室表面孔隙水压力最大;当入射频率为2.0时,洞室表面孔隙水压力最小,洞室表面最大孔隙水压力出现在洞室顶部。     

具有球形空腔饱和土分数导数 黏弹性土的动力特性

将土骨架视为具有分数阶导数本构关系的黏弹性体,基于Biot两相饱和介质模型,建立具有球形空腔饱和分数导数黏弹性土体稳态振动的控制方程。通过引入势函数,得到球对称情形下具有球形空腔饱和分数导数黏弹性土体的位移、应力和孔隙流体压力等解析表达式。考察分数导数模型参数和饱和土参数等对土体振动特性的影响,结果表明,流体压缩性对饱和土体的动力特性有显著影响,而土骨架压缩性和流-固耦合系数的影响相对较小;分数导数阶数对土体动力特性的影响与材料参数比的取值有关。同时,边界不排水条件下饱和土体的动力响应大于排水条件下饱和土体的动力响应。     

爆炸荷载作用下深埋圆形隧洞的饱和黏弹性土-衬砌系统动力响应

假定土骨架服从标准线性固体黏弹性本构关系,研究了深埋圆形隧洞的饱和黏弹性土-弹性衬砌耦合系统在轴对称爆炸作用下的瞬态动力响应。首先,基于饱和土的Biot模型和衬砌的弹性理论,通过引入势函数和Laplace变换,利用弹性衬砌和饱和黏弹性土界面处的连续性条件以及边界条件,得到饱和黏弹性土体和弹性衬砌位移、应力和孔隙水压力等在Laplace变换域中的解析解。其次,利用Laplace数值Crump逆变换得到耦合系统在时间域的动力响应,数值分析了不同土体模型下土体-衬砌耦合系统的径向位移和环向应力以及土体孔隙水压力等。结果表明：对不同土体模型的土体-衬砌耦合系统,其在爆炸载荷作用下的动力响应性态基本一致,但动力响应的振动周期和幅值等具有明显的差异。同时,对于饱和黏弹性土-弹性衬砌系统,土体黏性参数对土体径向位移和孔隙水压力有明显的影响,但对土体环向应力影响较小。     

饱和冻土中弹性波的传播特性

用混合物连续介质理论,考虑了土颗粒骨架、冰、水三相介质,选取土颗粒位移、孔隙水位移、孔隙水压和孔隙冰压为基本变量,采用Bishop有效应力原理,建立了饱和冻土多孔介质的弹性波弥散方程。经理论推导,给出了饱和冻土中弹性波的传播速度及衰减的解析表达式。通过数值算例,探讨了饱和冻土中两种压缩波（Pl波和P2波）及剪切波（S波）的波速和衰减与频率和孔隙率、含冰量等土参数的关系。通过参数分析研究了饱和冻土中3种体波的传播特性。     

有限二维饱和多孔介质因载荷诱发Biot固结的解析解

推导了有限矩形区域饱和多孔介质因表面载荷诱发的Biot固结的一个解析解。假设多孔介质为均匀各向同性和线弹性,并被单相流体所饱和;控制方程组采用不可压缩多孔介质模型;孔隙压力场采用狄利克雷边界条件,上下表面位移场符合物理边界,而左右侧面位移场边界条件则由人为特别给定。利用有限正余弦变换和拉普拉斯变换及数值反演获得了物理空间孔隙压力场和位移场的半解析解,其体现为双重级数和的封闭形式。最后以某软黏土层平面应变固结为例,利用有限元分析软件ABAQUS对所给出的解析解进行了验证,同时基于该解析解考察了孔隙压力场和位移场的时空演化规律。所给出的解析解可用于深入分析有限二维饱和多孔介质的流-固耦合力学行为。     

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

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

非饱和土中圆柱形衬砌隧道的瞬态动力响应

在以往关于圆柱形衬砌隧道的瞬态动力响应中,衬砌周围土体大多假定为弹性介质或饱和介质。然而,自然界中的土体大多为非饱和介质。考虑土体与衬砌结构的动力相互作用及动荷载引起的附加质量密度的影响,研究了瞬态荷载作用下非饱和土中无限长深埋圆柱形衬砌隧道的动力响应。基于多孔介质混合物理论和连续介质力学理论,建立了非饱和土中圆柱形衬砌隧道受到瞬态荷载作用时衬砌及周围土体的控制方程,利用Durbin数值反演法得到了衬砌及土体在时间域的动力响应。数值分析了饱和度对瞬态荷载下径向位移、径向应力、环向应力和孔隙水压力的影响。结果表明：饱和度对衬砌及周围土体的瞬态响应影响显著;饱和度对径向位移沿径向的衰减影响较小,对环向应力和孔隙压力沿径向的衰减影响较大。     

移动环形荷载作用下饱和土中圆形衬砌隧洞动力响应研究

基于Biot理论,研究了饱和土中带有衬砌的圆形隧洞在移动环形荷载作用下的动力响应。假定衬砌为弹性体,土体为饱和多孔介质,引入两类势函数来表示土体、孔隙水和衬砌的位移,使隧洞的控制方程解耦。结合边界条件及连续条件,通过傅立叶变换得到频率-波数域中衬砌和土体的应力、位移和孔隙水压力解答,最后用傅立叶积分逆变换得到时-空域中的数值解。计算并比较了3种隧洞模型（弹性土体隧洞、饱和土体隧洞和饱和土衬砌隧洞）的动力响应分析。数值分析结果说明：（1）移动荷载速度对3种隧洞动力响应均具有较大影响;（2）弹性土体隧洞和饱和土体隧洞的动力响应具有明显区别,所以在富水地区的隧洞动力响应中土体应该视为饱和土体;（3）衬砌对隧洞动力响应有较大影响,故隧洞的动力分析中不能忽略衬砌作用。     

Modelling of fluid expulsion and deformation behaviour of dewatering sediments

A finite element model is developed for modelling coupled fluid expulsion/deformation behaviour of dewatering sediments subjected to external loadings under isothermal conditions. The non-linear deformation behaviour of the sediment (soil) skeleton is based on the force equilibrium equation in which the constitutive relationship of stress and strain is implemented by the modified Cam-Clay model in soil plasticity. The fluid flow behaviour in the model is described by the generalized porous media flow equation. The model allows temporal and spatial variations of porosity and permeability. The fluid viscosity and density are assumed to be temperature-dependent. The model also allows the development of single and multiple faults, depending upon the material (sediment and fluid) properties, loading and boundary conditions. Procedures are implemented for (1) updating the material properties such as porosity, permeability, fluid density and viscosity and (2) the development of faults which allow the formation of high-permeability conduits for fluid flow. The solution algorithm for displacements of the sediments and the excess pore (fluid) pressure is based on a residual load technique to handle the non-linear (elastic-plastic) deformation behaviour of the sediment skeleton. The model can be applied to one- and two-dimensional problems. Examples of a plane strain saturated sediment layer subjected to stepwise horizontal tractions versus time are given.     

Coupled simulation of wave propagation and water flow in soil induced by high‐speed trains

The purpose of this paper is to simulate the coupled dynamic deformation and water flow that occur in saturated soils when subjected to traffic loads, which is a problem with several practical applications. The wave propagation causes vibrations leading to discomfort for passengers and people in the surroundings and increase wear on both the vehicle and road structure. The water flow may cause internal erosion and material transport in the soil. Further, the increased pore water pressure could reduce the bearing capacity of embankments. The saturated soil is modelled as a water‐saturated porous medium. The traffic is modelled as a number of moving wheel contact loads. Dynamic effects are accounted for, which lead to a coupled problem with solid displacements, water velocity and pressure as primary unknowns. A finite element program has been developed to perform simulations. The simulations clearly demonstrate the induced wave propagation and water flow in the soil. The simulation technique is applicable to railway as well as road traffic. Copyright © 2007 John Wiley & Sons, Ltd.     

Elastic consolidation around a point sink embedded in a half-space with anisotropic permeability

The complete solution is presented for the transient effects of pumping fluid from a point sink embedded in a saturated, porous elastic half-space. It is assumed that the medium is homogeneous and isotropic with respect to its elastic properties and homogeneous but anisotropic with respect to the flow of pore fluid. The soil skeleton is modelled as a linear elastic material obeying Hooke's law, while the pore fluid is assumed to be incompressible with its flow governed by Darcy's law. The solution has been evaluated for a particular value of Poisson's ratio of the solid skeleton, i.e. 0.25, and the results have been presented graphically in the form of isochrones of excess pore pressure and surface profile for the half-space. The solutions presented may have application in practical problems such as dewatering operations in compressible soil and rock masses.     

Coupling of solid deformation and pore pressure for undrained deformation—a discrete element method approach

This paper presents a numerical scheme for fluid‐particle coupling that uses the discrete element method by taking into consideration solid deformation and pore pressure generation. A new water particle element is introduced to calculate pore water pressure due to porosity changes. The water particle element has the same size and shape as the solid element and experiences the same amount of deformation. On the basis of the effective stress principle at the element contact, the total force is equal to the sum of the force transmitted through the solid element contact and the water particle force due to pore water pressure. Analytical solutions of traditional soil mechanics problems, such as isotropic compression and consolidated triaxial undrained test, are used to quantitatively validate the proposed model. The numerical results show good agreement between the model and the analytical solutions. The model therefore provides an effective method to calculate pore pressure in a porous medium in discrete modeling.     

Hydro‐mechanical modeling of cohesive crack propagation in multiphase porous media using the extended finite element method

In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of fully saturated porous materials

Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd.