固结系数为非常数时的砂井固结方程差分解.pdf

固结系数为非常数时的砂井固结方程差分解 程江涛，刘莉 中国地质大学(武汉)工程学院，武汉（430074） E-mail：cjtlaotou@163.com 摘 要:砂井主要用以增加饱和软土排水途径，使地基在荷载下加速固结。目前对砂井地基 排水固结问题多采用常固结系数的太沙基—伦杜立克理论，然而软土属于高压缩性土，随着 软土被压密，孔隙比将会发生较大变化，土的性质也会随之发生变化，因此，软土固结系数 的变化不容忽视。基于以上考虑，根据固结系数为非常数时的太沙基—伦杜立克方程，采用 (Crank—Nicolson)六点对称差分格式对其初始条件和边界条件进行有限差分离散，建立了有 限差分方程组。在此基础上，运用 Matlab 编制了砂井固结方程差分程序，并以佛山市和顺 —北潜公路干线软基处理为例，研究了其在不同荷载不同时间段的孔隙水压力以及固结度的 变化。 关键词: 砂井；固结方程；固结系数；固结度；孔隙水压力；有限差分法 1. 引言 砂井主要用以增加饱和软土的排水途径，使地基在荷载下加速固结。砂井问题属于空间 问题，其变形属性（1～3 维)取决于荷载的分布形式，而其中渗流属于三维问题，与荷载分 布无关。要较为精确的分析砂井地基的固结，必须考虑荷载的分布和渗流的空间属性。然而， 要得到这样的解答是极其困难的。为得到简单实用的解答，通常假设变形只发生在竖向，而 渗流呈空间轴对称性，从而将整个砂井地基的固结问题转化为单井固结问题。 目前对砂井地基的排水固结问题多采用固结系数为常数的太沙基—伦杜立克理论[1],[2]， 然而软土属于高压缩性土，随着软土压密，孔隙比将会发生较大变化，土的性质也会随之发 生变化，因此，软土固结系数的变化不容忽视。基于以上考虑，本文针对固结系数为非常数 时的太沙基—伦杜立克方程[3],[4]，采用(Crank—Nicolson)六点对称差分格式对其初始条件和 边界条件进行了有限差分离散，建立了有限差分方程组。在此基础上，运用 Matlab 编制了 砂井固结方程差分程序，并以佛山市和顺—北潜公路干线软基处理为例，研究了其在不同荷 载不同时间段的孔隙水压力以及固结度的变化。 2. 太沙基—伦杜立克方程差分格式 2.1 太沙基—伦杜立克方程 假设土体只发生竖向变形，径向渗流，瞬时加载下的径向固结方程式为: ? 2u 1 ?u ?u C ( r + r ) = r （1） h ?r 2 r ?r ?t 式中: Ch —径向固结系数，为有效应力的函数: C h = cons tan t （2.a） Ch = α vσ ′ + β v （2.b） Ch = Cv (σ ′) （2.c） 求解条件为： 1. u = u （3.a） r t=0 r0 - 1 - 2. ur = urw （3.b） r=rw ?u 3. r = 0 （3.c） ?r r=rh 式中：u0 为 t = 0 时，仅考虑径向渗流时，地基中任一点孔压: urw 为仅考虑径向渗流时， 砂井内任一点