路基沉降预测的拓展指数曲线模型

2010-01-26秦尚林陈善雄许锡昌

铁道标准设计 2010年2期

秦尚林,陈善雄,许锡昌

(中国科学院武汉岩土力学研究所岩土力学与工程国家重点实验室,武汉 430071)

1 概述

随着国民经济的快速发展,高速铁路在经济发展中将具有非常重要的战略意义,我国高速铁路建设目前已全面展开,为保证其运行过程中的快速、舒适和安全,路基工后沉降量控制是一个关键的环节,因此,科学、准确地预估线下结构物沉降变形成为高速铁路路基工程中非常重要的一个课题[1～3]。

目前沉降预测方法主要有三类,一是经典的分层总和法,二是基于Biot固结理论的数值计算方法,三是通过现场实测沉降资料,采用曲线拟合的方法进行估算。前两种方法由于假定条件和确定计算指标的试验技术上的问题,常常导致计算的沉降量与实测值相差甚远,因此,第三种方法成为工程中最为常用的方法[4]。指数曲线法作为一种常用的曲线拟合方法,具有简单实用的优点,但是其前提是假定荷载是一次施加,这与实际情况不符[5],因此,该模型往往后期沉降拟合较好,而前期沉降拟合不够稳定,与实际测量值相差较大[6]；另外,高速铁路路基沉降变形量级较小,且相对波动较大,指数曲线法对沉降变形观测数据的单调性有着严格的要求,局部数据的小幅度起伏变化都可能导致无法进行预测计算,在武广铁路客运专线沉降评估实践过程中,也发现大部分沉降观测断面直接用指数曲线法无法计算[7]。

本文主要针对常规的指数曲线模型的不足,进行适当的改进,以期达到不仅可以预测小变形情形,而且可以有效改善其初期阶段预测误差大的缺点。

2 指数曲线法的改进

2.1 常规指数曲线法

根据工程实践经验,从沉降曲线的“拐点”B开始,路基沉降量s与时间t的关系规律可近似表示为指数曲线模型(图1)。在不考虑次固结沉降的情况下,其方程可由式(1)表示

(1)

式中，s0为t0时刻的沉降量；st为t时刻的沉降量；s∞为最终沉降量；η为待定常数。

(2)

对式(2)取对数得

(3)

图1 指数曲线法推算沉降示意

2.2 拓展指数曲线模型

前面已经指出,指数曲线法的主要缺点：一是填筑阶段预测误差大,二是对小变形情况不适用。为改善其预测性能,并能将填筑期观测数据纳入分析时间段以内,为此,在常规的指数曲线法基础上引入荷载系数的概念,对模型进行必要的修正,以期克服上述缺点,并且适用于沉降量小、数据起伏波动大的高速铁路路基沉降的预测分析。

假定荷载增量、加载速率变化不大的情况下,沉降变形的增量与荷载增量成正比,从而引入荷载系数ξ,将指数曲线模型式(1)修正为

(4)

3 拓展指数曲线模型求解方法

在实测数据中取两点：(t1,s1),(t2,s2),要求s1

(5)

(6)

式(5)除以式(6)可得

(7)

式(7)中只有u一个未知数,通过Excel可求得u的值,然后将其代入式(4)得

(8)

由式(8)可知,st与1-e-ut成线性关系,由此经过拟合可求得s∞。将求得的u与s∞代入式(4)即可求得任意时刻的沉降量。

4 拓展指数曲线模型适用性探讨

改进后的指数曲线法究竟能否达到预期要求,适用于高速铁路沉降预测,下面以武广铁路客运专线DK1 524+140.46～DK1 524+481.57段路基沉降板观测断面为例来进行说明。

本段路基全长341.11 m,位于丘谷区，地势平坦开阔,丘坡表层为粉质黏土,下伏呈砂土状的全风化花岗岩、强风化花岗岩和弱风化花岗岩。采用了换填、搅拌桩、CFG桩等多种方式进行地基加固处理。填筑于2007年11月7日完成,填土层厚2.1～2.7 m,共设置了10个沉降板观测断面,沉降观测从2007年7月10日开始,至2008年5月24日为止,沉降观测312 d,由沉降板测得的总沉降量在5.40～11.30 mm。

限于篇幅,仅以断面DK1 524+300沉降板观测数据为例进行计算分析,而直接列出其他断面的计算结果。DK1 524+300断面沉降板的填土高度-时间-沉降(H-t-s)曲线见图2。

图2 DK1 524+300断面沉降板H-t-s曲线

在沉降曲线上选取两点(t1,s1)=(61,4.30),(t2,s2)=(312,11.30),代入式(7)可得

(9)

运用Excel可求得u=0.006 63,将其代入式(8)可得

该式通过直线拟合可得s∞=13.878 mm。

将u及s∞的值代入式(4)可得

st=13.878(1-e-0.006 63t)

(10)

根据式(10),新模型预测曲线如图3所示。

图3 DK1 524+300断面实测值与模型预测曲线

拓展指数曲线模型计算结果见表1。

表1 拓展指数曲线模型与常规指数曲线预测结果对比

常规的指数曲线法不能对该断面进行预测,而拓展指数曲线模型可进行预测。拓展指数曲线模型拟合曲线与实测曲线吻合程度很好,相关系数高。

该区段10个断面计算结果见表2。

表2 拓展指数曲线模型沉降预测结果汇总

从表2可以看出：拓展指数曲线模型能够计算DK1 524+140.46～DK1 524+481.57段的10个路基断面的数据,并且计算所得的相关系数高,均在0.99以上,拟合数据与实测数据的误差平方和范围为5.828～13.878 mm2,吻合程度很好；改进前,指数曲线法不能直接计算这些断面,改进后,全部能够计算,而且能够得到很好的预测效果。

5 拓展指数曲线模型与传统方法对比

为进一步考察改进指数曲线方法的预测效果,针对上述10个断面,采用双曲线法、三点法计算结果与之进行对比。双曲线法、三点法的具体分析方法参见文献[5]。

运用双曲线法及三点法对DK1 524+140.46～DK1 524+481.57区段10个沉降板观测断面工后沉降进行了拟合分析与预测,分析结果汇总见表3。

表3 实测结果与预测结果汇总

由表3可以看出,3种方法对该段10个观测断面的预测结果相差不是很大,都可以较好地进行沉降预测分析。相比较而言,拓展指数曲线模型与双曲线模型拟合的相关系数非常接近,均优于三点法的相关系数；而拓展指数曲线模型计算的误差平方和大部分介于其他两种方法之间,优于三点法而比双曲线法略大；用拓展指数曲线模型计算的最终沉降量也基本上介于双曲线法和三点法之间,比三点法大而比双曲线法小,可见拓展指数曲线模型拟合精度是比较高的,能够得到很好的预测效果。