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

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

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

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

DYNAMIC RESPONSE OF SATURATED POROUS ELASTIC FOUNDATION UNDER POROSITY ANISOTROPY$^{\bf 1)}$

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

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

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

轴对称热-力埋置载荷下半空间饱和孔隙介质的动力响应

基于修正的Biot热弹性本构理论,得到了饱和多孔介质热-力耦合的动力学控制方程.针对半空间孔隙介质在内置简谐热-力轴对称载荷作用下的动力问题,利用Hankel变换获得了响应的解析表达式,并利用Han-kel逆变换进行数值求解,分析了埋置深度、表面热边界条件等对响应的影响规律.结果表明:孔隙水压力在载荷作用处上方有负压出现;环向、径向及竖向应力在载荷作用处发生突变,且在载荷作用处上方均出现拉应力;在载荷作用处下方,孔隙水压力及竖向、径向和环向应力均随着深度的增大而减小.当内置热载荷仅设置温差,而无外热源输入时,温差对孔隙介质中的响应几乎不产生影响.当孔隙介质表面绝热时,孔隙水压力小于表面等温情况下的值.     

中心舱体与薄壁梁刚-柔耦合系统的热弹性-结构动力学分析

空间柔性结构受太阳热流冲击而诱发的振动是导致航天器失效的典型模式之一,准确预测结构热致振动的响应及稳定性是卫星设计的基础。针对常见的中心舱体与附属薄壁杆件组成的空间结构,提出了考虑刚-柔耦合、耦合热弹性和耦合热-结构三重耦合效应的热致振动分析理论模型。其中,刚-柔耦合是指舱体姿态角、顶端集中质量转动与柔性附件运动的耦合;耦合热弹性是指应变率与温度场的耦合;耦合热-结构是指舱体转动及结构变形与薄壁杆件吸收太阳热流的耦合。基于热弹性理论和Lagrange方程,推导了传热和运动的耦合方程;采用Laplace变换方法并使用Routh-Hurwitz稳定判据推导了稳定性边界方程。结果表明,该模型能够更为准确的给出热致振动响应及稳定性预测。     

半无限大体广义电磁热弹耦合的二维问题

基于Lord和Shulman广义热弹性理论,研究了热、电可导的半无限大体电磁热弹耦合的二维问题。半无限大体受热和外加恒定磁场的作用,文中建立了电磁热弹性耦合的控制方程,利用正则模态法求解得到了所考虑物理量的解析解,并用图形反映了各物理量的分布规律,从分布图上可以看出,介质中出现了电磁热弹耦合效应,各物理量的非零值仅在一个有限的区域内。     

层状分数阶黏弹性饱和地基与梁共同作用的时效研究

饱和地基与梁共同作用问题的研究在力学领域及工程界都具有重要意义. 采用分数阶Merchant模型研究饱和地基的流变固结, 该模型比常用整数阶黏弹性模型更能精确反映地基的时变特征. 基于层状正交各向异性黏弹性饱和地基的固结解答, 采用有限元法与边界元法耦合的方法, 研究梁与分数阶黏弹性饱和地基的共同作用问题. 依据Timoshenko梁理论将梁离散为若干单元, 进而得到梁的总刚度矩阵方程; 将黏弹性地基固结问题的精细积分解答作为边界积分的核函数, 采用边界元法建立地基柔度矩阵方程; 结合梁与地基接触面的位移协调条件以及力的平衡条件, 通过有限元法与边界元法的耦合, 最终求得层状分数阶黏弹性饱和地基与Timoshenko梁共同作用的解答. 将本文地基退化为Kelvin地基进行计算, 并与已有文献中的算例进行对比, 二者具有很好的一致性. 在此基础上, 探讨分数阶次和地基成层性对梁与黏弹性饱和地基共同作用的影响. 结果表明: 分数阶次高的黏弹性饱和地基的固结速率明显更快; 对于层状地基, 加固表层土体能有效控制地基整体沉降, 并减小差异沉降. 实际工程中, 应充分考虑饱和地基流变及土体分层性的影响, 以准确分析梁与地基的共同作用过程.      

非饱和多孔介质中混合元法的化学-热-渗流-力学耦合的本构模拟

在文献[1]中建立的多孔介质中化学-热-渗流-力学(CTHM)本构模型基础上,针对文献[2]建立的非饱和多孔介质中热-渗流-力学耦合分析的混合有限元方法,发展了非饱和多孔介质中混合元的化学-热-渗流-力学(CTHM)耦合本构模拟算法。采用非关联流动多重屈服准则模拟非饱和多孔介质的材料非线性行为。推导了u-pw-pa-T形式的包含了耦合率本构方程积分的向后欧拉映射算法和一致性弹塑性切线模量矩阵(单元刚度矩阵)的混合元一致性算法。本文给出了临界状态线(CSL)和状态边界面(SBS)两个屈服准则的一致性算法。数值结果显示了本文所发展的混合元耦合本构模拟算法在模拟由热、化学、力学荷载共同引起的多孔介质中化学-热-渗流-力学(CTHM)耦合行为的能力和有效性。     

多孔介质的热-电-化-力学耦合理论及应用

本文针对自然界中不同种类的多孔介质,评述了包含化学效应的多场耦合力学问题的国外内研究现状.介绍了研究多场耦合问题的理论架构,并基于唯象理论和热力学理论,建立了一般性的热-电-化-力学多场耦合理论,在此基础上简化为化学-力学耦合理论,利用相应的控制方程和本构关系, 给出了线性耦合系统的变分原理,并证明了化学-力学耦合理论架构的封闭性. 基于所提出的化-力耦合模型,通过数值算例解释了多孔介质中的化学-力学耦合现象.最后对多孔介质,特别是对活体生物软组织中的多场耦合研究中存在的问题进行了讨论,并展望了本领域的未来发展趋势.     

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.     

Effects of the dynamic wheel–rail interaction on the ground vibration generated by a moving train

Based on Biot’s fully dynamic poroelastic theory, the dynamic responses of the poroelastic half-space soil medium due to quasi-static and dynamic loads from a moving train are investigated semi-analytically. The dynamic loads are assumed to be generated from the rail surface irregularities. The vehicle dynamics model is used to simulate the axle loads (quasi-static loads) and the dynamic loads from a moving train. The compatibility of the displacements at wheel–rail contact points couple the vehicle and the track–ground subsystem, and yield equations for the dynamic wheel–rail loads. A linearized Hertzian contact spring between the wheel and rail is introduced to calculate the dynamic loads. Using the Fourier transform, the governing equations for the poroelastic half-space are then solved in the frequency–wavenumber domain. The time domain responses are evaluated by the fast inverse Fourier transform. Numerical results show that the dynamic loads can make important contribution to dynamic response of the poroelastic half-space for different train speed, and the dynamically induced responses lie in a higher frequency range. The ground vibrations caused by the moving train can be intensified as the primary suspension stiffness of the vehicle increases.     

Numerical analysis of the isolation of the vibration due to Rayleigh waves by using pile rows in the poroelastic medium

The isolation of the vibration due to harmonic Rayleigh waves using pile rows embedded in a saturated poroelastic half-space is investigated in this study. Based on Biot’s theory and the potential function method, the free field solution for Rayleigh waves along the surface of the poroelastic half-space is derived first. The fundamental solution for a harmonic circular patch load applied in the poroelastic half-space are obtained in terms of Biot’s theory and the integral transform method. Using Muki’s method and the fundamental solution for the circular patch load as well as the Rayleigh waves solution for the poroelastic half-space, the second kind of Fredholm integral equations in the frequency domain for pile rows are derived. Numerical solution of the integral equations yields the dynamic response of the pile–soil system to incident Rayleigh waves. Influences of various parameters on the vibration isolation effect of piles rows are investigated numerically. Numerical results suggest that for the same vibration source, the same pile rows will produce a better vibration isolation effect for the poroelastic medium than for a single phase elastic medium. Also, stiffer piles tend to have better vibration isolation effect than flexible piles. Moreover, the pile length and the spacing between neighboring piles in each pile row have significant influence on the vibration isolation effect of pile rows.     

直角坐标下横观各向同性饱和半空间体动力响应问题的解析解

基于双重Fourier变换技术，成功求解了直角坐标系下，横观各向同性弹性饱和多孔介质的三维Biot动力方程，得到了以固体骨架位移分量和孔隙流体压力为基本未知量的积分形式一般解．进而，用一般解给出了饱和介质的总应力分量表达式。在此基础上，研究了在任意分布的表面竖向和水平谐振力作用下，横观各向同性饱和半空间体的动力响应问题。数值结果表明，采用各向同性饱和介质的动力学模型，不能准确描述具有明显各向异性特性的饱和土地基的动力性能。     

横观各向同性饱和介质非轴对称动力响应分离变量解Separation of variables solution of non-axisymmetrical dynamic response of transversely isotropic saturated porous media

首先将横观各向同性饱和弹性多孔介质非轴对称问题的Bio t波动方程,变换为适宜于进行分离变量法求解的形式;然后在非轴对称简谐激励下,用分离变量法得到Bio t方程的一般解,即用分离变量法求得了多孔介质位移和应力分量的解析表达式;并给出了半空间横观各向同性饱和弹性多孔介质在表面竖向简谐荷载作用下表面竖向位移的数值分析结果,得出载荷对30倍受载半径以外的区域几乎无影响的结论。同时表明了本文的分析方法是切实可行的。     

移动荷载下上覆弹性板饱和地基的动力响应

用半解析法对移动条形荷载作用下上覆弹性板饱和两相弹性介质的动力响应问题进行了研究。由忽略土粒压缩和土体自重的Biot波动方程出发,对荷载进行Fourier展开。假设响应函数形式,利用待定系数法求解了考虑固液耦合作用的两相介质在移动荷载作用下的土体位移,有效应力及孔压表达式。求解过程中考虑了土体和弹性板之间的相互作用,并假设土体与弹性板的竖向位移相等。通过计算给出了各参数对土体位移和孔压响应和的影响。数值结果表明弹性梁刚度对位移和孔压响应有较大影响。     

饱和地基上基础动力阻抗函数的一种分析方法

提出了一种基础阻抗函数的时域求解方法, 并用于饱和地基情形. 通过给基础输入脉冲位移, 应用集中质量显式有限元方法结合局部透射人工边界, 得到地基施加于基础的反力时程, 然后根据阻抗函数的定义, 应用傅立叶变换求得基础阻抗函数. 通过算例, 与 Halpern 的结果进行了对比, 验证了该方法的有效性. 并以一方形基础为例, 分别讨论了泊松比、孔隙率、渗透系数和埋深对无量纲动力柔度的影响. 与以往方法相比, 该方法可以考虑复杂地基情形和不规则基础形式.      

Response of railway track system on poroelastic half-space soil medium subjected to a moving train load

Based on the dynamic poroelastic theory of Biot, dynamic responses of a track system and poroelastic half-space soil medium subjected to moving train passages are investigated by the substructure method. The whole system is divided into two separately formulated substructures, the track and the ground, and the rail is described by introducing the Green function for an infinitely long Euler beam subjected to the action of moving axle loads of the train and the reactions of the sleeper. Sleepers are represented by a continuous mass and the effect of the ballast is considered by introducing the Cosserat model for granular medium. Using the double Fourier transform, the governing equations of motion are then solved analytically in the frequency-wave-number domain. The time domain responses are evaluated by the inverse Fourier transform computation for a certain train speed. Computed results show that the shape of the rail displacements of the elastic and poroelastic soil medium are in good agreement with each other of the low train velocity, but the result of the poroelastic soil medium is significantly different to that of the elastic soil medium for the high train velocity which is higher than Rayleigh-wave speed in the soil. The influence of the soil intrinsic permeability on soil responses is discussed with great care in both time domain and frequency domain. The dynamic responses of the soil medium are considerably affected by the fluid phase as well as the load velocity.     

A numerical model for the isolation of moving-load induced vibrations by pile rows embedded in layered porous media

A numerical model is developed to analyse the isolation of moving-load induced vibrations using pile rows embedded in a layered poroelastic half-space. Based on Biot’s theory and the transmission and reflection matrices (TRM) method, the free wave field solution for a moving load applied on the surface of a layered poroelastic half-space and the fundamental solution for an harmonic circular patch load are determined. Using Muki and Sternberg’s method, the second kind of frequency domain Fredholm integral equations for the dynamic interaction between pile rows and the layered poroelastic half-space are derived. The time domain solution is recovered via inverse Fourier transform in order to obtain the amplitude reduction ratio and thus assess the vibration isolation efficiency of pile rows. A special case of the present model shows good agreement with an existing solution. Numerical results of this study show that the speed of moving loads has an important influence on the isolation of vibrations by pile rows: the same pile rows can achieve better isolation efficiencies for higher speed loads than for lower speed loads. Pile rows embedded in a two-layered poroelastic half-space with a softer overlying layer usually generate better vibration isolation effects than those with a stiffer overlying layer. Finally, better isolation vibration may be realized by increasing the pile length and decreasing the net spacing between neighboring piles in a pile row.