室内与现场反演固结系数差异特征与原因分析

2020-06-05赵俊明陈琦刘垚邓永锋

公路与汽运 2020年2期

赵俊明,陈琦,刘垚,邓永锋

(1.苏交科集团股份有限公司,江苏南京 210017；2.东南大学交通学院,江苏南京 210096)

连云港地区紧靠黄海,土层多为淤泥、淤泥质黏土、淤泥质亚黏土及淤泥混砂层,属于饱和的正常压密软黏土,具有压缩性高、沉降量大、排水固结慢、地基稳定性较差的特点。根据对已建道桥的调查,大部分道桥因地基沉降不均而产生严重的路面开裂和桥头跳车,这是由施工期间对工后沉降预估不足所致。在软土地基上建设公路面临复杂的软土固结沉降问题,固结系数是体现固结过程快慢的指标,是软基变形设计计算的关键参数。

目前常用的获取固结系数的方法主要有室内固结试验法、孔压静力触探法和现场反演分析计算法。对这3种方法之间的差异,不少学者开展了相应研究。如陈新彦等以珠江三角洲大量软土地区试验结果为依据,比较了地基沉降反演、室内固结试验和孔压静力触探试验所得固结系数的相互关系；胡荣华等基于深圳地区软土工程实例,对比了室内试验法、现场试验法、反演分析法、间接推算法结果的异同；顾绍付等以具体工程土样的固结试验和渗透试验数据为依据,用多种方法计算土样的固结系数,分析了两种参数之间的联系；孟高头等开展孔压静力触探和室内固结试验,对比了两种固结系数间的异同；邓永锋等根据Hong&Onitsuka的扰动度定义,给出了土样扰动对室内固结试验结果的影响。已有研究表明现场沉降反演法与孔压静力触探法的固结系数较接近,而室内试验所得固结系数比现场沉降反演法小1～2个数量级。但对于这种差异产生原因的研究较少,导致选取合理的固结系数作为施工依据比较困难。该文基于连云港临海高等级公路建设期沉降观测数据,从土样扰动、渗透系数各向异性及固结模式三方面分析差异产生原因。

1 固结系数获取方法

1.1 室内试验法

室内试验一般采用基于Terzaghi一维固结理论的解析法求竖向固结系数,固结方程为：

式中：Cv为固结系数,Cv＝k/(mvγw)；k为渗透系数；mv为体积压缩系数；γw为水的重度。

为便于求解,Terzaghi假定k和mv在土体压密固结过程中保持不变,即k、mv和Cv在固结过程中均为常数。

1.2 孔压静力触探法

探头贯入土层时,其周围土层中的应力、应变及孔隙水压力均产生重分布,探头贯入土层所产生的超静孔隙水压力的消散接近于水平径向扩散,故孔压静探在土层中所估算的为水平向固结系数Ch。按Terzaghi固结理论,孔压消散的轴对称方程为：

式中：u为孔隙水压力；r为任一点距探头中心的水平距离；t为消散时间。

将式(2)代入满足超孔压的初始条件和边界条件,可求解某一深度土层的水平向固结系数：

式中：T为时间因数；r0为探头半径。

目前国内外倾向于取固结度U为50%时的T50计算固结系数。

1.3 Asaoka沉降反演法

Asaoka法数学形式相对简单,参数易确定,仅需少量数据即可进行图解法求解固结系数。Asaoka基于Mikasa提出的一维大变形的连续性方程∂ε/∂t＝Cv(∂2ε/∂z2),构建了ε的级数解：

式中：T、F为时间t的函数；T′、T″,F′、F″分别为T、F对时间t的一阶和二阶导数。

根据边界条件,可得到Asaoka解析解：

式中：S i为在同一荷载作用下第i个时间段的沉降；β0和β1为参数,双面排水时lnβ1＝－6Cv/H2Δt,单面排水时lnβ1＝－2Cv/H2Δt。

2 结果对比分析

2.1 室内固结系数与反演法结果对比

根据汾灌(汾水—灌云)、沪苏浙(上海—江苏—浙江)高速公路的勘察和沉降资料,采用Asaoka法反演地基综合固结系数,并与现场原状土室内试验所得固结系数进行对比,结果见图1。

图1 Asaoka反演与室内平均固结系数对比

由图1可知：现场反演的平均固结系数比室内所得固结系数大10倍左右。

2.2 孔压静力触探法与沉降反演法结果对比

孟高头等基于珠江三角洲部分典型断面的现场沉降观测资料反算得出该地区软土土层的固结系数,并与孔压静力触探试验所得结果进行比较,结果见表1。

由表1可知：孔压静力触探反演结果与实测沉降资料的反算值很接近,处于同一数量级。该结果与蔡国军等对连盐高速公路沿线4个典型监测断面的反演结果一致。

表1 孔压静力触探消散试验值和现场沉降实测反算值对比

3 差异原因分析

从勘察取样扰动、固结系数各向异性及固结模式3个角度对几种方法所得固结系数差异的产生原因进行分析。

3.1 取样扰动对固结系数的影响

邓永锋等对连云港海相软土进行薄壁取样,讨论了土样扰动程度、应力水平与固结系数的关系,得到固结系数随扰动度SD和应力的变化(见图2)。

从图2可以看出：在扰动程度较低的情况下,屈服前的固结系数比屈服后的大1～2个数量级,且扰动度低的屈服前固结系数比扰动度高的屈服前固结系数大1～2个数量级。由于室内试验在取样运输过程中存在应力、扰动度变化的情况,土样应力状态会从屈服前变为屈服后,扰动程度增加,进而导致室内试验结果比现场试验求得的固结系数小1～2个数量级。

3.2 固结系数各向异性

根据已有研究成果,反演法所得固结系数与孔压静力触探法所得水平向固结系数Ch处于同一个数量级,差异很小。孟高头等以珠江三角洲大量钻探试验结果为依据,对孔压静力触探所得水平向固结系数Ch和室内固结试验所得Cv进行一元线性回归分析,结果见图3。

图2 固结系数随扰动度SD的变化

图3 Ch-Cv线性相关关系散点图及回归曲线

由图3可知：Ch、Cv间的关系为Cv＝0.086 5Ch＋0.685 9,相关系数r为0.905 4,回归方程的拟合程度很好,Ch和Cv具有明显的线性关系；Ch比Cv大1～2个数量级。

3.3 固结模式

Asaoka法是基于太沙基一维固结理论提出的反演法,求得的是包含水平、竖向固结的现场整体固结系数,而室内试验由于侧限仅能获得竖向固结系数。为进一步分析固结模式对固结系数反演结果的影响,以临海高等级公路连云港段某断面为例建立有限元数值模型,利用Asaoka法反演固结系数Cv1,同时根据压力不同设置等效室内试验固结系数按深度的加权平均值Cv2,通过渗透系数和固结系数的正相关关系建立水平、竖向渗透系数的比值kh/kv与固结系数的关系,分析kh/kv的改变对固结系数的影响。模型中,路基高度为3 m,硬壳层厚度为2 m,软土层深度为15 m,下卧层深度为10 m,地下水位在硬壳层下边界,水平影响区域为50 m(见图4)。土层基本参数见表2,堆载过程见图5。现场kh/kv一般为2～5,在这个范围内,设定竖向渗透系数kv不变,根据不同kh/kv及软土层厚度(15、20和25 m)下的沉降历时曲线进行反演分析。反演固结系数与竖向加权固结系数对比见表3。

图4 计算模型

图5 路基填土历时曲线

表2 土层计算参数

从表3可看出：1)Cv1/Cv2随着kh/kv的升高升高。kh/kv＝1时,Asaoka法平均反演固结系数为室内固结试验的5倍左右；kh/kv＝5时,Asaoka法平均反演固结系数为室内固结试验的10倍左右。2)Cv1/Cv2随软土层厚度的增加而减小。即使水平向渗透系数与竖向渗透系数相等,反演固结系数Cv1也为竖向固结系数Cv2的5倍,说明在现场水平和竖向二维固结模式中,水平向固结始终占主导地位,且随着kh/kv的增大而增大。3)在kh/kv的一定比值,反演固结系数Cv1比竖向固结系数Cv2大1～2个数量级,即反演法所得固结系数比室内试验所得固结系数大1～2个数量级。