桩周饱和黏土非齐次渗流固结解

2021-10-20高子坤傅长荣

水利水电科技进展 2021年5期

高子坤, 傅长荣

(1.丽水学院土木工程系，浙江丽水 323000; 2.浙西南地质灾害防治研究所，浙江丽水 323000)

超静孔隙水压力的消散速度及程度影响基础的承载力，并反映饱和土层固结度及压缩量随时间的变化规律。该问题相关的研究主要包括初始超静孔隙水压力的确定，反映超静孔隙水压力时空变化规律理论或数值解答的推导[1-2]及其在工程与环境问题解决中的应用等内容。

Sagaseta[3-4]探讨了不排水条件下桩-土作用引起的地基土位移、应力求解等问题；高子坤等[1,5-6]对半空间中压桩挤土造成的附加位移、应力、超静孔隙水压力及其消散过程展开研究，得到相应问题的解答，但针对土体超静孔隙水压力消散固结问题的解答局限于齐次边界；陈文[7]分析了饱和黏土中沉桩挤土机理及效应；刘时鹏等[8]采用K0固结排水模型，推导了饱和黏土中压桩挤土问题的应力场并得到初始超静孔隙水压力解答；曹权等[9]针对饱和软黏土中压桩挤土问题，引入对数应变参数描述大变形及软化等特征，推导土中孔扩张造成的应力场，并分析了超静孔隙水压力分布规律；姚笑青等[10-11]对沉桩在饱和黏土中引起的超静孔隙水压力进行了估算。

本文针对空间轴对称非齐次固结问题，建立非齐次固结定解条件并推导解答，计算得到超静孔隙水压力随时间、空间变化的模拟数据，绘制时间-平均固结度曲线、水头等值线和流速矢量图。结合工程实例，应用实测的桩基承载力随时间变化数据换算平均固结度，以验证理论推导的正确性。

1 解答推导

空间轴对称问题固结微分方程为[12]

(1)

式中：u为超静孔隙水压力；Ch、Cv分别为水平和竖直方向的固结系数；z、r分别为竖直和水平(径)向坐标。

建立固结模型如图1所示，包括桩端平面(地下或地表)、桩尖平面、桩土接触界面以及初始超静孔隙水压力影响范围等。

图1 固结模型

考虑固结开始时，场地附近的河流或水库水位影响，建立桩端平面水头边界；考虑基桩以承压含水砂砾石为持力层，或为加速固结排水设计砂井抽水降压等措施的影响，建立桩尖平面水力梯度边界，即式(1)的边界条件为

(2)

初始条件为

u|t=0=f(r,z)

(3)

式中：φ(t)、ψ(t)分别为桩端边界处水头压力函数和桩尖平面边界处水力梯度函数；f(r,z)为任意函数；rw为桩孔半径；re为f(r,z)影响半径；H为桩长。

对定解方程式(1)进行齐次化运算，令：

u(r,z,t)=V(r,z,t)+W(z,t)

(4)

其中W(z,t)=zψ(t)+φ(t)

依据式(2)及W(z,t)表达式可得桩尖平面水力梯度Wz(H,t)=∂W(H,t)/∂z=ψ(t)，桩端平面水头W(0,t)=φ(t)。

将式(4)代入式(1)可得：

(5)

RZT′-ChDZT-CvTRZ″=0

(6)

引入参数μ、λ，对式(6)分离变量可得：

(7)

求解式(7)可得：

(8)

(9)

T(t)=Ae-λk,iCht

(10)

结合初始条件V(r,z,0)=f(r,z)-W(z,0)，将式(8)～(10)代入V(r,z,t)可得：

(11)

其中λk,i=αi2+nμk

Mi=Y0(αir)-Y0(αire)J0(αir)/J0(αire)

式中：αi(i=1,2,3,…)为特征方程J1(αirw)·Y0(αire)-J0(αire)Y1(αirw)=0的特征值；J0、J1分别为0阶和1阶1类Bessel函数；Y0、Y1分别为0阶和1阶2类Bessel函数。

依据u(r,z,t)=V(r,z,t)+W(z,t)可得式(1)的理论解：

u(r,z,t)=V(r,z,t)+zψ+φ

(12)

由式(12)可得t时刻平均固结度理论计算式为

(13)

2 工程实例

某试桩承载力实测值及其换算的平均固结度如表1所示[1,12]。

表1 承载力及平均固结度

f1(r,z)=βΔσ0+αfΔτ0

(14)

其中Δσ0=(Δσr+Δσθ+Δσz)/3

式中：Δσr、Δσθ、Δσz、Δτrz为沉桩引起的桩周土体的应力增量，由文献[1]确定；A为Skempton孔压参数；β为Henkel孔隙压力系数，对于饱和黏土β=1。

为便于比较及验证，也对式(13)所示的初始条件进行固结渗流问题计算：

f2(r,z)=γ′z

(15)

2.1 固结渗流计算

依据式(12),对水头和水力梯度为φ=3 m、ψ=-0.4和φ=3 m、ψ=0.4两种工况进行计算，绘制超静孔隙水压力等值线和渗流速度矢量如图2～5所示，图2和图4基于初始条件式(14)，图3和图5基于初始条件式(15)。由图2～5比较分析可知：固结开始时刻(0 d)，等值线及流速矢量的空间分布与初始条件吻合，验证了级数解答的正确性与收敛性；

图2 超静孔隙水压力等值线(初始条件式(14)，φ=3 m，ψ=-0.4;单位：kPa)

图3 超静孔隙水压力等值线(初始条件式(15)，φ=3 m，ψ=-0.4;单位：kPa)

图4 超静孔隙水压力等值线(初始条件式(14)，φ=3 m，ψ=0.4;单位：kPa)

图5 超静孔隙水压力等值线(初始条件式(15)，φ=3 m，ψ=0.4;单位：kPa)

在固结稳定阶段(110 d)，等值线及流速矢量的空间分布由边界取值确定，为稳定渗流；固结渗流中间阶段(10～50 d)，渗流的空间分布和流速矢量的大小与方向受边界条件和初始条件共同作用，随时间变化，为非稳定渗流；负水力梯度将提高排水固结速度，并最终在桩周土中产生负超静孔隙水压力,即超固结，如图2和图3所示；正水力梯度的作用与负水力梯度相反，使土层的超静孔隙水压力无法完全消散，如图4和图5所示。