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

2021-11-09周凤玺张雅森曹小林牟占霖

力学学报 2021年7期

关键词：非饱和渗透系数边界条件

周凤玺张雅森曹小林牟占霖

* (兰州理工大学土木工程学院,兰州 730050)

† (西部土木工程防灾减灾教育部工程研究中心,兰州 730050)

** (东南大学土木工程学院,南京 211189)

引言

弹性半空间的波动响应一直是弹性动力学领域的研究热点,其中Lamb 问题最具代表性.随着研究的深入,不同集中荷载形式(点源和线源、表面和内部等)作用下单相弹性介质的Lamb 问题已经取得了比较完备的体系.近年来,有关多相多孔介质的Lamb 问题逐渐受到人们的重视.自Biot[1-2]建立了两相介质的波动方程后,国内外学者针对饱和半空间Lamb 问题已经取得了一系列研究成果,主要包括荷载作用于半空间表面[3-5]和半空间内部[6-9]以及层状地基[10-14]等不同方面的动力响应研究.相对于饱和土,在工程建设中大量涉及到的是处于地下水位以上的非饱和土体,而已有研究表明,介质中孔隙气体的存在对其动力响应行为、弹性波传播特性以及能量传输产生巨大影响,因此研究非饱和半空间的动力学行为在岩土工程、地震工程等领域有着重要的理论和应用价值.

由于非饱和多孔介质物理力学特性的复杂性,使得对非饱和半空间Lamb 问题及能量传输特性的研究成果较少.王春玲等[15-17]采用积分变换法和消元法求得了非饱和地基受竖向简谐荷载作用下的稳态响应积分变换解,但其最终解的形式十分复杂,不便于应用.徐明江等[18-20]以三相多孔介质模型为基础,通过引入双变量本构关系,采用解析法研究了简谐荷载作用下非饱和土地基的动力响应问题,给出了积分形式的解答,但并未考虑颗粒间吸应力对非饱和半空间动力学特征的影响.Zhang 等[21]在假定土骨架为多孔弹性连续介质,且具有均匀性和各向同性的基础上,通过应用Fiourier 展开技术和Hankel 积分的方法得到了在内部激励作用下的非饱和土中动态格林函数解,但上述研究成果均未讨论不同边界条件对非饱和半空间表面的动力学响应特征及能量传输特性的影响规律.

本文在已有研究成果的基础之上,考虑非饱和土中粒间吸应力的作用,结合质量守恒方程、动量平衡方程及有效应力原理等基本方程,运用Helmholtz 分解法,在柱坐标系下建立了非饱和半空间的动力学控制方程.分别考虑透水(气)和不透水(气)两种边界条件,对频域内的轴对称问题进行求解,得到了非饱和半空间表面受到竖向集中简谐荷载作用下的解析解答,并且通过参数分析讨论了在不同饱和度、不同振动频率以及不同渗透系数下的动力响应特性和能量传输的变化及其影响规律,以期为不同边界条件下非饱和半空间的表面振动问题提供参考依据.

1 基本方程

如图1 所示的非饱和半空间表面受到频率为ω,幅度为q0的垂直简谐荷载作用,考虑轴对称性,问题的基本方程包括如下6 类.

图1 半空间计算模型Fig.1 Half-space calculation mode

(1) 时域内固体骨架动量平衡方程

(2) 时域内孔隙流体运动平衡方程

式(1)和式(2)中,σr,σθ,σz,τrz分别为代表单元体上的总应力和剪应力;ur,uθ,uz分别为径向位移、环向位移和竖向位移;wr,wθ,wz分别为孔隙中液体的径向位移、环向位移和竖向位移;vr,vθ,vz分别为孔隙中气体的径向位移、环向位移和竖向位移;pl和pa分别为孔隙水压力和孔隙气压力;分别为固相、液相和气相的表观密度.

考虑孔隙中的流体渗透形式符合达西定律,用固有渗透系数K表征气相和液相的渗透系数kl和ka分别为

则式(2a)～式(2d)中系数bl和ba可表示为如下形式

式中,n为孔隙率;Sr为饱和度;g为重力加速度;ηl和ηa分别为液相和气相的动力黏度系数.

(3) 有效应力原理

有效应力原理是土力学的核心,目前关于非饱和土中有效应力原理大致可分为单变量理论[22]、双变量理论[23]和复合变量理论[24].其中,Lu 等[25]在考虑微观颗粒间作用力和有效应力的基础上提出了吸应力表示的有效应力公式

式中,σs=−(pa−pl)Se表示粒间吸应力,Se为有效饱和度,Se=(Sr−Sw0)/(1−Sw0),Sw0为残余饱和度.吸应力表示的有效应力公式是非饱和介质中两相流体压力的函数,主要受饱和度变化的影响,消除了Bishop 有效应力原理中对参数χ的依赖性.

(4) 质量守衡方程

忽略各相之间的质量交换,非饱和多孔介质各相的质量守衡方程可描述为[26-28]

(5) 本构方程

考虑固相颗粒的可压缩性,弹性各向同性非饱和介质的本构方程为

式中,α=1−Kb/Ks,Ks和Kb分别为固相颗粒以及骨架的体积压缩模量; λ 和 µ为 Lame 常数.

(6) 渗流连续方程

式中,系数aij的具体表达式详见附录A.本构方程(7)和渗流连续方程(8) 的具体推导过程分别见附录B 和附录C.

2 问题的求解

2.1 控制方程

考虑简谐荷载作用下各位移分量的形式如下

将式(9)代入孔隙流体运动平衡方程(2)中,并结合式(1)、式(7)和式(8),整理后可得频域内非饱和半空间的动力控制方程

根据Helmholtz 分解定理,引入柱坐标系下的两个势函数Φ和Ψ表示ur和uz如下

将式(12a)分别代入式(10a)和式(11)中,可得

考虑到因变量的时间导数同二阶空间导数的乘积与其时间导数同一阶空间导数的乘积相比,前者是高阶小量[21].因此,利用式(8)可以得到

将式(16)代入式(13),可得

将式(17)代入式(15),整理后可得

式中

方程(14)和方程(18)即为以势函数表示的频域内控制方程.在得到势函数解答的基础上,可通过基本方程得到各位移分量、应力分量等物理量的解答.

2.2 解析解答

利用分离变量法求解微分方程(14)和(18),可得

式中,I0(ζ11r),K0(ζ11r),I0(ζ12r),K0(ζ12r),I0(ζ2r),K0(ζ2r)分别为第一类和第二类零阶贝塞尔函数;

待定系数A,B,C,D由具体的边界条件确定.

考虑如下边界条件

根据边界条件(21)～ (23)可得

因此,对于势函数可重写成以下形式

将式(26)和式(27)代入式(12)和式(17)中,最终可获得非饱和多孔介质中孔隙压力及位移分量的表达式为

将孔隙压力和位移的结果代入本构方程(7)中,可得各应力分量的表达式为

对于本文所讨论问题的描述,在非饱和半空间表面(z=0)分别考虑透水(气)和不透水(气)两种边界条件.

(1) 当半空间表面排水(气)时,土体表面孔隙气压力和孔隙水压力为0,即

此时,式(28)～式(30)中各参数分别为

(2)当半空间表面不排水(气)时,土体表面孔隙流体和土骨架之间的相对位移为0,即

此时孔隙中流体和土骨架间的位移关系由式(2)得

此时,式(28)～式(30)中各参数可表示为

2.3 能量场特性

通常情况下,半空间表面单位面积的能量传播情况可由其表面牵引力和质点运动速度的积表示[29].因此,对于本文所考虑的非饱和多孔介质材料,单位面积上的能量,可由式(34)表示,结合两类边界条件(1)和(2)可求得不同透水(气)条件下,非饱和半空间表面受到简谐荷载作用时能量的传输性质.

3 解答的有效性验证

在经典Lamb 问题的分析研究中,学者们往往采用位移解的形式进行描述.为验证本文计算结果的准确性,利用式(29)得出的非饱和半空间表面的竖向位移uz和水平位移ur的解析解答,同文献[29]在饱和半空间中的计算结果进行对比,其分析结果十分接近,说明本文所得结果可以和经典饱和半空间理论很好地衔接,进一步证明本文计算结果的有效性,如图2 所示.

图2 有效性验证Fig.2 Validity verification

4 数值分析与讨论

为讨论不同边界条件下相关参数对非饱和半空间动力响应和能量传输特性的影响规律,本文将通过数值算例分析在不透水(气)条件、透水(气)条件下土体表面位移及能量变化受到饱和度、振动频率、渗透系数的影响情况.数值算例中所选取的计算参数如表1 所示.

表1 计算参数[17]Table 1 Calculation parameters[17]

在激振圆频率ω=1 rad/s 时,图3 绘制了非饱和半空间表面处为不透水(气)条件、透水(气)条件时,饱和度变化对其表面位移及能量传输特性的影响曲线.由图3 可见,非饱和半空间表面位移幅值会随着饱和度的增大而增大.这主要是由于随着饱和度的升高,非饱和介质中基质吸力降低从而引起粒间吸应力降低,使半空间抵抗外力变形的能力减弱,导致半空间表面位移幅值会呈现出逐渐增大的趋势.在不透水(气)条件下,整体的位移幅值低于透水(气)条件下的位移幅值,当饱和度较低时,孔隙内部存在大量气体,由于气体本身有很强的可压缩性,因此当饱和度较低时,两种不同边界条件下的位移幅值相差很小;当饱和度较高时,土中孔隙水的含量明显升高,非饱和介质抵抗变形的能力也会随之提升,因此当饱和度较高时,不透水(气)条件下的位移幅值会较为明显的低于透水(气)条件下的位移幅值.同位移幅值的变化情况类似,半空间中的能量同样呈现出随着距离振源位置的增大而振荡下降的趋势,且当表面不透水(气)条件下,孔隙流体压力占比相对更大,但由于总应力没有变化,导致有效应力占比相对较小,因此在不透水(气)条件下半空间表面受外荷载作用时的总能量值依然小于透水(气)条件下的能量.

图3 饱和度对位移和能量的影响曲线Fig.3 Influence curve of saturation on displacement and energy

为了分析荷载激振圆频率对非饱和半空间表面位移及能量传输特性的影响,图4 绘出了不同频率下相关物理量的变化曲线.从图4 可以看出,随着激振频率的逐渐增加,不论是径向还是竖向位移的幅值均逐渐减小.因为在荷载振动频率较小时,透水(气)条件下的地表孔隙水(气)压力更容易消散,所以在透水(气)和不透水(气)两种不同边界条件下的位移幅值显现出一定的差异,荷载振动频率较大时的现象与之相反.且由于施加的外力水平不变,因此非饱和半空间中的能量变化会呈现出相似的变化趋势.

图4 激振频率对位移和能量的影响曲线Fig.4 Influence curve of excitation frequency on displacement and energy

在ω=1 rad/s,饱和度Sr=0.6 时,图5 给出了固有渗透系数的变化对半空间表面位移幅值和能量传输特性的影响曲线.由图5 可见,随着固有渗透系数的逐渐降低,骨架位移也随之减小.随着距离振源位置的逐渐增加,位移幅值呈现出振荡下降的现象.当渗透系数很低时,两种边界条件的性质趋于一致,因为在条件(1)的情况下,表面虽透水(气),但由于半空间内部的孔隙水(气)压力难以快速消散,依然会对土骨架产生一定的支持作用,所以位移幅值会随着渗透系数的降低而呈现出一定的下降趋势.在不透水(气)条件下,孔隙流体会持续影响半空间的表面位移,且孔隙中流体和土骨架之间没有相对位移,因此不透水(气)边界条件下非饱和半空间表面的位移幅值会略低于透水(气)边界条件下的位移幅值,且当土骨架位移幅值减小时,总能量也会呈现出减小的现象,当渗透系数下降至K<1×10−13m/s 时,地表位移幅值受渗透系数影响趋于一极限值,并且随着渗透系数的逐渐降低,这两种不同边界条件所产生的宏观现象会逐渐趋于一致,但其差异性将一直存在.