孔隙率各向异性下饱和多孔弹性地基动力响应

由于不同的沉积条件和应力状态, 天然土体通常表现出一定的各向异性特征. 文章研究地基上表面受温度载荷和机械载荷时, 孔隙率各向异性参数变化对饱和多孔弹性地基热-水-力耦合动力响应问题的影响. 基于Lord-Shulman广义热弹性理论, 结合孔隙率各向异性基本假设, 建立了孔隙率各向异性饱和多孔弹性地基热-水-力耦合动力响应模型, 利用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度分布的解析表达式并加以图示. 正则模态法是一种利用加权残差求得解析解的方法, 相较于其他方法能快速求解偏微分方程. 当孔隙率各向异性参数为1时, 可将该各向异性耦合动力响应模型退化为热-水-力耦合动力响应模型验证该地基模型的合理性. 着重分析了孔隙率各向异性参数变化对不同物理量的影响. 结果表明: 孔隙率各向异性参数变化对物理量均有一定影响. 在地基上表面受温度载荷作用时, 对超孔隙水压力和竖向应力影响最为明显; 在地基上表面受机械载荷作用时, 对超孔隙水压力和温度影响明显. 整体而言, 无论地基上表面受何种载荷, 随着各向异性参数增大, 峰值逐渐减小, 在地基深度增加方向峰值所在位置向靠近地基上表面方向移动.      

饱和多孔黏弹地基热-水-力耦合动力响应分析

天然土体由于沉积条件和应力状态不同, 往往会表现出一定的流变性. 本文研究地基上表面受外载荷作用时, 渗透系数和孔隙率变化对饱和多孔黏弹性地基热-水-力耦合动力响应问题的影响. 基于Biot波动方程、达西定律和Lord-Shulman广义热弹性理论, 并引入了考虑黏弹性松弛时间因子的Kelvin-Voigt黏弹性模型研究地基上表面受热/力源作用时, 孔隙率和渗透系数变化对均质各向同性饱和多孔黏弹性地基中所考虑的各无量纲量的影响. 根据不同的边界条件采用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度的解析表达式, 结合算例分析了不同变量对各物理量的影响. 正则模态法是一种加权残差法, 可不经正、反积分变换将方程快速解耦并消除数值反变换的局限性. 结果表明: 无论何种载荷作用时, 载荷频率变化对所有考虑的物理量均有明显的影响; 孔隙率和渗透系数均对无量纲超孔隙水压力有明显的影响, 当仅考虑热载荷作用时, 孔隙率和渗透系数变化对无量纲温度均无影响. 正则模态法可广泛应用于岩土工程领域, 尤其适用于商业建筑、高速铁路和公路能源基础的热、力学特性研究中. 该研究结果可为工程施工奠定一定的理论基础, 具有一定的指导性意义.      

Research on coupled thermo-hydro-mechanical dynamic response characteristics of saturated porous deep-sea sediments under vibration of mining vehicle

The excessive deformation of deep-sea sediments caused by the vibration of the mining machine will adversely affect the efficiency and safety of mining. Combined with the deep-sea environment, the coupled thermo-hydro-mechanical problem for saturated porous deep-sea sediments subject to the vibration of the mining vehicle is investigated. Based on the Green-Lindsay(G-L) generalized thermoelastic theory and Darcy's law, the model of thermo-hydro-mechanical dynamic responses for saturated porous deep-sea sediments under the vibration of the mining vehicle is established. We obtain the analytical solutions of non-dimensional vertical displacement, excess pore water pressure, vertical stress, temperature, and change in the volume fraction field with the normal mode analysis method, and depict them graphically. The normal mode analysis method uses the canonical coordinate transformation to solve the equation, which can quickly decouple the equation by ignoring the modal coupling effect on the basis of the canonical mode. The results indicate that the vibration frequency has obvious influence on the vertical displacement, excess pore water pressure, vertical stress, and change in volume fraction field. The loading amplitude has a great effect on the physical quantities in the foundation, and the changes of the physical quantities increase with the increase in loading amplitude.     

荷载作用下多孔饱和地基的热-水-力耦合动力响应分析

基于广义热弹性理论,结合达西定律,对Biot波动方程进行修正,研究了一个受到荷载作用的多孔饱和地基的热-水-力多场耦合动态响应问题。建立了多孔饱和地基在荷载作用下的热-水-力耦合模型及控制方程,该模型可退化为热弹性耦合模型。采用正则模态法求解,得到了问题的解析解,讨论了热-水-力耦合模型和热弹性耦合模型的区别,分析了荷载频率变化对地基中各物理量的影响。最终给出了无量纲的竖向位移、超孔隙水压力、竖向应力和温度等物理量的分布规律。     

条形荷载下梯度非均匀非饱和土的动力响应分析

基于多孔介质混合物理论, 建立了梯度非均匀非饱和土地基模型, 研究了条形荷载作用下梯度非均匀非饱和土地基的动力响应问题. 通过傅里叶积分变换和Helmholtz矢量分解原理, 获得频域内非饱和土地基动力响应问题的通解, 结合回传射线矩阵法和边界条件, 求解获得了非均匀非饱和土层中位移、应力以及孔隙压力的计算列式. 假设沿深度方向梯度非均匀非饱和土的物理力学性质按幂函数连续变化, 通过数值傅里叶逆变换得到了非均匀非饱和土地基中的应力、位移以及孔隙压力等物理量的数值解, 分析讨论了土体非均匀性对非饱和土介质动力响应的影响规律. 结果表明: 土体非均匀性显著改变了非饱和土中竖向位移、正应力和孔隙压力在其深度方向上的振动模态, 其中孔隙气压在其深度方向的振动频率随着梯度因子的增加而不断增大, 波峰值不断靠近地表处附近; 竖向位移随着梯度因子的增大不断减小; 正应力和孔隙水压随着梯度因子的增大先增大后减小, 并且土体非均匀程度越高, 正应力与孔隙水压的幅值越大.      

非饱和土-隧道系统动力响应计算的波函数法

针对非饱和地基土中埋置隧道的三维动力响应计算问题, 提出了波函数法.采用无限长的Flügge薄壁圆柱壳模拟圆形隧道衬砌,采用流、固、气组成的三相介质模拟非饱和地基土体.分别采用分离变量法以及Helmholtz矢量分解定理求解薄壁圆柱壳的振动控制方程与非饱和土的波动方程.根据隧-土交界面与地表面处的应力、位移以及孔隙流体压力等边界条件,利用平面波与柱面波的转换性质,实现了隧道内作用单位简谐载荷时隧道衬砌与土体系统动力响应的耦合求解.通过与既有单相弹性介质2.5维有限元-边界元法、两相饱和多孔介质2.5维有限元-边界元法以及三相非饱和介质Pip in Pip半解析法的计算结果进行对比, 验证了本文计算方法的可靠性. 最后,基于该方法, 通过算例分析了不同饱和度下非饱和土-隧道系统的动力响应特征.结果表明, 饱和度对土体动位移与超孔隙水压力的幅值响应有较大影响.该方法的非饱和地基土参数退化后,也可用来计算和分析饱和地基土或单相弹性地基土与隧道系统的动力响应.      

半空间饱和土内置点载荷作用下的热弹性波动

基于Biot波动理论及热弹性动力理论,利用已建立的饱和多孔弹性介质热流固耦合控制方程,研究半无限地基在内置点热力源作用下的动力响应问题. 求解过程引用Hankel变换技术,得到了热力源作用下土体中温度增量、应力、位移和孔隙水压力的积分形式解答.利用Hankel数值逆变换得到计算结果,分析了热流固耦合条件下激振频率对竖向位移和孔隙水压力响应的影响. 对热流固耦合、热弹性和多孔弹性模型计算结果进行了比较.      

黏弹性准饱和土中球空腔动力特性

在频率域内研究了内水压力作用下分数导数型黏弹性准饱和土中球空腔的稳态动力响应. 通 过引入与孔隙率有关的应力系数合理地确定了介质和孔隙水共同承担的内水压力值. 将土骨 架和衬砌分别视为具有分数导数本构模型的黏弹性体和多孔柔性材料, 基于Biot两相介质模 型, 通过引入位移势函数解耦得到了内水压力作用下分数导数型黏弹性准饱和土中半封闭球 形空腔的位移、应力和孔隙水压力解析表达式. 考察了物性和几何各参数对球形空腔动力响 应的影响, 结果表明: 分数导数本构模型更合理地描述了土体的动力学行为; 饱和度对应力 和孔隙水压力影响较大, 而对位移影响较小.     

砂岩岩石力学参数各向异性研究

对胜利油田不同深度的四口井的砂岩岩芯进行了三轴压缩试验和动态岩石力学参数试验。根据获得的垂直地层方向和平行地层方向的试验数据,计算了各岩石力学参数的各向异性系数。通过分析围压以及深度对各向异性系数的影响,给出了砂岩峰值强度各项异性系数、静泊松比各向异性系数、纵波波速各向异性系数、纵横波速比和动泊松比各向异性系数与深度的拟合公式。发现峰值强度、纵波波速和纵横波速比的各向异性系数都随着深度的增加而减小,而泊松比各向异性系数随着深度的增加而增大。动、静弹性模量各向异性系数与深度相关性都很差,动弹性模量各向异性较弱,静弹性模量相对较强。     

成桩后黏性土中超孔压和有效应力试验研究

针对预制桩沉桩挤土效应问题,采用现场原位试验和数值模拟方法,分析了桩土界面超孔压和竖向有效应力随时间的变化规律,并应用水力压裂理论、孔穴扩张理论和修正剑桥(modified Cam-clay, MCC) 模型,揭示了成桩后桩周土体固结过程的超孔压沿桩长的分布形式。结果表明：桩周超孔压大小同孔压与传感器入土深度成正比;在外载荷作用下,黏性土中超孔压和竖向有效应力分别随时间呈“增–减–平”和“平–增–平”的三段式变化趋势;成桩后约25天,超孔压基本消散完毕,竖向有效应力趋于稳定最大值;基于修正理论公式的数值模型得出,不同载荷下,桩侧中下部出现负超孔压,距桩的水平距离和距桩端竖向距离越小,超孔压越大。可为同类静压沉桩挤土效应问题提供参考依据。

     

爆炸荷载作用下饱和半空间中衬砌隧道弹性动力响应IBIEM模拟

基于Biot两相介质理论,采用一种高精度间接边界积分方程法（IBIEM）研究了饱和半空间中浅埋衬砌隧道在内部爆炸荷载作用下的瞬态弹性动力反应。通过典型算例,给出了爆炸荷载作用下隧道附近地表位移、衬砌动应力、围岩径向位移和衬砌表面孔隙水压的时程响应,并对比分析了饱和半空间和全空间中隧道动力响应的区别。研究表明：覆土层厚度对浅埋隧道-围岩整体动力响应特征具有明显影响;衬砌表面透水状态对爆炸荷载的时程响应的影响不显著;随半空间饱和介质孔隙率增加,围岩受隧道内部爆炸影响程度降低,衬砌承担的爆炸作用增大;当和直达波、衬砌内部反射波的峰值叠加作用时,半空间表面反射波对衬砌隧道拱顶附近响应影响显著,使得衬砌动应力幅值、径向位移相比深埋情况大幅度增加。     

Discrete simulation and micromechanical analysis of two-dimensional saturated granular media

In this study, a novel approach to incorporate the pore water pressure in the discrete element method (DEM) to comprehensively model saturated granular media was developed. A numerical model was constructed based on the DEM by implanting additional routines in the basic DEM code; pore water pressure calculations were used with a two-dimensional (2D) model to simulate the undrained behavior of saturated granular media. This model coupled the interaction of solid particles and the pore fluid in saturated granular media. Finally, several 2D undrained shear tests were simulated. The test results showed that the model could predict the response of the saturated granular soil to shear loading. The effect of initial compaction was investigated. Biaxial tests on dense and loose specimens were conducted, and the effect of the initial density on the change in shear strength and the volume change of the system was investigated. The overall behavior of loose and dense specimens was phenomenologically similar to the real granular material. Constant volume tests were simulated, and the results were compared to those from the coupled model. Induced anisotropy was micromechanically investigated by studying the contact force orientation. The change in anisotropy depended on the modeling scheme. However, the overall responses of the media obtained using the coupled and constant volume methods were similar.     

Discrete simulation and micromechanical analysis of two-dimensional saturated granular media

In this study, a novel approach to incorporate the pore water pressure in the discrete element method （DEM） to comprehensively model saturated granular media was developed. A numerical model was constructed based on the DEM by implanting additional routines in the basic DEM code; pore water pressure calculations were used with a two-dimensional （2D） model to simulate the undrained behavior of satu- rated granular media. This model coupled the interaction of solid particles and the pore fluid in saturated granular media. Finally, several 2D undrained shear tests were simulated. The test results showed that the model could predict the response of the saturated granular soil to shear loading. The effect of initial compaction was investigated. Biaxial tests on dense and loose specimens were conducted, and the effect of the initial density on the change in shear strength and the volume change of the system was inves- tigated. The overall behavior of loose and dense specimens was phenomenologically similar to the real granular material. Constant volume tests were simulated, and the results were compared to those from the coupled model. Induced anisotropy was micromechanically investigated by studying the contact force orientation. The change in anisotropy depended on the modeling scheme. However, the overall responses of the media obtained usinz the couoled and constant volume methods were similar.     

A Water Retention Curve Model for the Simulation of Coupled Thermo-Hydro-Mechanical Processes in Geological Porous Media

This paper presents a new water retention curve (WRC) model for the simulation of coupled thermo-hydro-mechanical processes in geological porous media. The model simultaneously considers the impact of porosity and temperature on suction, for both wetting processes and drying processes. The model is based on an idealization of porous geological media as having an isotropic and homogeneous microscopic pore structure. Suction is expressed as a function of degree of saturation, porosity, surface tension of the water–air interface, and the length of air bubble perimeter of the pores per unit area on a random 2D cross-section of the medium. The tension of water–air interface is written as a function of temperature, and the length of perimeter of the water–air interface of the pores becomes a function of porosity and degree of saturation. The final equation of the new WRC is a function of suction, effective degree of saturation, temperature, porosity, pore-gas pressure, and the rate of degree of saturation change with time for both wetting and drying processes. The model was used to fit experimental data of the FEBEX bentonite, with good agreements between measured and calculated results.     

Material instabilities of anisotropic saturated multiphase porous media

This paper analyses the material instability of fully saturated multiphase porous media. On account of the fact that anisotropic mechanical behaviours are widely observed in saturated and partially saturated geomaterials, the anisotropic constitutive model developed by Rudnicki for geomaterials is used to model the anisotropic mechanical behaviour of the solid skeleton of saturated porous geomaterials in axisymmetric compression test. The inertial coupling effect between solid skeleton and pore fluid is also taken into account in dynamic cases. Conditions for static instability (strain localisation) and dynamic instability (stationary discontinuity and flutter instability) of fully saturated porous media are derived. The critical modulus, shear band angle for strain localisation, and the bound within which flutter instability may occur are given in explicit forms. The effects of material parameters on material instability are investigated in detail by numerical computations.     

白洋淀内水体的流动特征研究

研究了风对白洋淀内水体循环流动的影响，研究指出：(1)风应力可看作是水体循环 的主要驱动力. (2)在风应力的作用下，水体表层的流动主方向与风的方向一致，下层水体 的流动方向与表层方向相逆或近似相逆. (3)表层以下水体流动复杂，次表层水平流动有明 显的双涡结构. (4)淀底地形对水体垂向流动影响很大，垂向速度变化复杂，垂向上有双涡 流动现象. 从动力学角度研究了淀内水体流动的情况.     

VERTICAL VIBRATION ANALYSIS OF AXISYMMETRIC SATURATED SOIL

Based on Biot’s dynamic consolidation equations,by means of Laplace-Hankel transform technology,the integral solutions of stress and displacement in saturated soil with subjacent rock-stratum under axisymmetric arbitrary excitations were derived.In- fluence of the reflected wave generated by the boundary was revealed.Numerical results indicate that the vibration frequency has some effect on the vertical displacement of satu- rated soil.The vertical displacement at the surface of saturated soil lags in phase with the load.Furthermore,the dynamic permeability coefficient of saturated soil has significant effect on the vertical displacement at the initial stage of load applied,but when the load becomes stable,the effect is inapparent.     

非饱和土半空间Lamb问题及能量传输特性

地球表面绝大多数土层处于非饱和状态, 故采用传统饱和两相介质理论进行动力学分析时, 结果往往与实际情况不符. 针对这一问题, 本文以非饱和半空间作为研究对象, 基于连续介质力学和多孔介质理论, 考虑非饱和多孔介质中各相的质量守恒方程、动量守恒方程、本构方程以及有效应力原理等基本方程, 建立了以骨架位移、孔隙水压力和孔隙气压力为基本未知量的动力学控制方程. 针对非饱和半空间表面在竖向集中简谐荷载作用下的动力学响应及能量传输问题, 建立了频域内经典Lamb问题的轴对称计算模型, 采用Helmholtz分解法, 通过引入势函数Φ和Ψ表示骨架的位移分量, 结合本构方程获得了不同边界条件下半空间表面位移场和能量场等物理量的解析解答, 并通过数值算例对荷载参数(激振频率)、材料参数(饱和度、渗透系数)等影响因素进行了分析与讨论. 结果表明: (1)饱和度的升高或者激振频率下降, 都会提高非饱和半空间的表面位移幅值; (2)当渗透系数下降至一临界值时, 地表位移幅值会趋于一极限值, 并且透水(气)边界与不透水(气)边界条件下渗透系数的影响表现出明显的差异性.      

饱和土地基在简谐荷载作用下的解析分析

基于多孔介质混合物理论,用解析的方法研究了不可压饱和土地基受到简谐荷载作用下的动力响应问题。利用Fourier积分变换求解耦合方程组,得到了二维饱和土介质在简谐荷载作用下的通解。针对表面透水的具有下卧基岩的饱和土层以及半无限饱和土地基的边界条件,获得了固体骨架位移、孔隙流体位移、固体骨架有效应力以及孔隙流体压力的积分形式解答,并通过数值算例分析了饱和土地基在简谐荷载作用下的响应。     

Steady state responses of poroelastic half-space soil medium to a moving rectangular load

The dynamic steady state responses of a poroelastic half-space soil medium subjected to a moving rectangular load are investigated analytically/numerically. The full dynamic poroelastic theory of Biot is employed, under the assumption of an incompressible solid grain and neglecting the apparent mass density. Using the Fourier transform, the governing equations of motion are then reduced to a system of four coupled ordinary differential equations which are solved semi-analytically. Soil vertical displacements, accelerations and pore water pressures induced by moving load are calculated. Computed result shows that load velocity and intrinsic permeability of the soil medium shows an apparent effect on its dynamic responses and pore water pressures.