矩形移动荷载作用下饱和-非饱和土双层地基的动力响应分析1)

2022-03-12李奎奎赵建昌王立安

力学与实践 2022年1期

李奎奎赵建昌王立安

(兰州交通大学土木工程学院，兰州 730070)

地基在移动荷载作用下的动力响应一直是岩土工程和交通工程领域的重点关注课题。随着运输业的发展，车辆的行驶速度越来越快，载重也越来越大，在此背景下研究地基土在移动荷载作用下动力响应问题对于道路结构设计和环境振动影响分析都具有重要意义。地基土体是一种具有多相组成和骨架孔隙结构的复杂介质，在动力作用下由于土体内部存在多场耦合效应，使其动力响应分析变得非常复杂。将地基土体视为固体骨架和孔隙水组成的两相饱和介质，利用Biot 多孔介质波动理论进行动力响应分析已经相当成熟[1-4]。徐斌等[5]根据Biot 波动理论，采用传递、反射矩阵方法研究了移动荷载作用下层状饱和土动力响应问题。王立安等[6]对非均匀饱和地基的动力特性做了研究。徐长节等[7]和周凤玺等[8]基于Vardoulakis 等[9]的理论研究了非饱和土中波的传播规律。李绍毅[10]采用半解析法研究了移动荷载作用下多层非饱和铁路地基的共振问题，分析了土体饱和度对非饱和地基振动的影响。Lu等[11]研究了矩形移动荷载作用下非饱和地基的动力响应。Fang 等[12]研究了列车荷载作用下非饱和地基的振动响应。Song 等[13]对埋置荷载作用下的非饱和地基动力响应问题做了研究。

上述关于非饱和地基的动力响应研究中，学者们都将地基考虑为单层非饱和介质，而实际自然界中的地基由于地下水的存在，地基底部土体往往是完全饱和的，只有顶部一定深度范围内受到环境气候的影响表现为非饱和状态。因此，对地基进行更精确的理论分析时，应将地基考虑为由下层饱和土和上层非饱和土层叠而成的半空间。

基于以上分析，文中将地基考虑为饱和土和非饱和土层叠而成的半空间模型，利用三重Fourier 变换和降阶法对饱和土和非饱和土的控制方程进行耦合求解，推导出地基在矩形移动荷载作用下的动力响应解。通过算例分析，对模型验证计算并给出了地表位移峰值沿上层非饱和土体饱和度和土层厚度的变化规律曲线，同时也给出了矩形移动荷载作用下地表(z= 0) 竖向位移的空间分布及竖向振动峰值与速度的关系曲线。

1 问题模型

如图1 所示，由非饱和土和饱和土层叠而成的半空间表面作用一矩形移动荷载Q，荷载移动速度为c，荷载分布长、宽分别为2a和2b。上层非饱和土厚度为H，下层饱和土为无限厚度。图中，λ∗,µ∗为饱和土对应的Lamb 常数，其表达式分别为

图1 计算模型Fig.1 Computational model

式中，ρs，ρw和ρa分别为非饱和土中固体、水和气体的质量密度。

2 非饱和土控制方程及求解

2.1 运动方程

根据混合物理论，考虑固-液-气之间的惯性耦合及毛细管压力，非饱和土的运动方程写为[12-13]

式中，u为非饱和土体中固体骨架的位移矢量(ui，i=x，y，z)；w，v分别为孔隙中水相和气相相对于固体骨架的相对位移矢量；θ为固体骨架的体应变。“(·)”、“(··)” 表示对时间t的一、二阶导数。pw，pa分别为孔隙水压力和孔隙气压力；α为Biot 系数，α= 1-Ksf/Ks，其中Ksf和Ks是土骨架和土颗粒的体积模量，λ+ 2µ/3 =Ksf；χ为有效应力参数，取值与有效饱和度Se相同，Se的表达式为：Se= (Sr-Sw0)/(1-Sw0)，其中Sw0为束缚水饱和度 (或残余饱和度)；g为重力加速度。∇和∇2分别为Hamilton 算子和Laplace 算子，即：∇=∂/x+∂/y+∂/z，∇2=∂2/x2+∂2/y2+∂2/z2。对于非饱和土的总孔隙压力p，有

2.2 本构方程

根据有效应力原理，非饱和土本构关系写为[13]

式中，σij(i,j=x,y,z)为总应力，εij是应变张量，δij为Kronecker 符号。

2.3 质量守恒方程

式中，A11～A14,A21～A24见参考文献[13]。

2.4 非饱和土控制方程求解

对于Hamilton 算子∇，有如下散度计算式

式中，θw和θa分别为孔隙水和气体的体应变。

对空间坐标x，y和时间坐标t引入三重Fourier变换，其变换关系为

利用式(10) 对式(2)～式(9) 进行三重Fourier变换和算子运算，并通过降阶法求解常微分方程组得到非饱和土的位移、孔压及应力通解，详细解法参见文献[11-12，14-15]。

式中，Θ为特征向量矩阵；±γ1，±γ2，...，±γ5为矩阵Θ对应的特征根；D1～D10为积分常数，通过边界条件确定。

3 饱和土控制方程及求解

3.1 饱和土控制方程

由Biot 波动理论，饱和土动力控制方程为[16]

(1) 运动方程

式中各参数定义参见文献[16]，其中上标“*” 表示饱和土对应的物理量。

3.2 饱和土控制方程求解

采用与非饱和土相同的求解方法，同样利用Fourier 变换和散度运算求解饱和土控制方程，再采用降阶法求解常微分方程组可得出饱和土的通解。此处，直接给出饱和土在变换域中的位移、孔隙压力和应力通解。

4 边界问题

4.1 矩形移动荷载的数学描述

对于半空间顶面作用的移动矩形分布荷载Q，利用Heaviside 阶跃函数将其写为

式中，δ( ) 表示Dirac-delta 函数。

4.2 边界条件

地表为透水、透气边界，则半空间表面(z= 0)处的边界条件为

4.3 边界方程求解

通过边界条件和连续条件组成的14 个方程可确定出积分常数D1～D14。对式(20) 和式(21) 做三重Fourier 变换，再将饱和土和非饱和土的位移、应力和孔压的通解代入，并分别取z=0，z=H，得到关于D1～D14的线性方程组

至此，系统中所有参数均已确定。将D1～D14回代到式(11)、式(12) 和式(16)、式(17)，即可计算出非饱和土层和饱和土层中的位移、应力和孔隙压力。对计算结果做Fourier 逆变换，可得到时间-空间域中各场量的计算结果。

5 验证与分析

5.1 算法验证

对前文推导结果进行编程计算，关于式(11)、式(12)和式(16)、式(17)的Fourier 逆变换可通过数值积分或离散Fourier 逆变换实现，本文中采用计算效率较高的离散Fourier 逆变换实现。引入移动坐标轴Xt，Xt=x-ct，使计算结果能以振源为中心分布。取上层非饱和土的饱和度Sr= 1，则半空间地基退化为单一饱和半空间，此时，有效应力参数χ= 1。则上、下层土体的特性参数值相同，土体参数取值与文献[16] 相同，按文献[16] 给定计算参数进行取值计算，其中主要参数的取值有：Lamb 常数µ取30 MPa，n取0.3；非饱和土中固体和水的质量密度ρs和ρw的取值分别为2700 kg/m3和1000 kg/m3。图2 为本文退化解与单层饱和半空间模型[16] 的结果对比，图中显示，退化解与文献解能够很好地吻合。(图中¯uz为无量纲竖向位移，¯uz=uzµ∗/(aq0)。)

图2 退化计算及结果对比(q0 =60 kN/m2，a=2 m，b=1 m)Fig.2 Calculation of degradation and comparison of results(q0 =60 kN/m2，a=2 m，b=1 m)

5.2 算例分析

图3 为移动荷载作用下地表(z= 0) 竖向位移的空间分布。图4 为地表位移峰值沿上层非饱和土体饱和度和土层厚度的变化规律曲线。图4(a)显示：当非饱和土层饱和度小于0.15 时，地表位移峰值随饱和度的增大而急速减小，当非饱和土层饱和度大于0.15 时则随之增大，当非饱和土层饱和度接近0.9 时(趋于完全饱和) 则再次急速减小，曲线的拐点位置随非饱和土层的厚度发生变化。图4(b) 显示：当非饱和土层的厚度小于0.5 m 时，地表振动位移发生剧烈波动。当非饱和土层的厚度大于0.5 m 时，随着非饱和土层的厚度增大，地表位移也随之增大。但当非饱和土层达到一定厚度后，振动位移将收敛于恒定值，该收敛值与单层非饱和地基的结果一致，说明地基完全体现为单层非饱和土地基，不再受下部饱和土层的影响，该临界厚度与土体饱和度相关，还能看出当地基上层土体为中度饱和时，层状土对地表位移峰值的影响较小。