0 引言
随着城市发展,地面交通拥堵已成为普遍现象,因此发展地下空间变得极为重要。矩形顶管法作为地下隧道的一种施工工法,有着空间利用率大、对周边环境影响小等优点; 但矩形顶管隧道在顶进结束后会继续产生地表沉降,且随时间增加而增大,进而危害到周围建( 构) 筑物。因此,对矩形顶管施工引起的工后地表沉降的研究具有重要意义。
目前关于矩形顶管施工引起土体变形的研究方法主要有: 经验法[1]、理论计算法[2]、数值模拟法[3 - 7]和实测数据分析法[8 - 11]等。经验法中: 林强强[1]证明Peck 公式可用于计算矩形顶管施工期间的地表沉降。理论计算法中: 王日东[2]采用随机介质理论得出矩形顶管推进引起的三维土体变形计算方法。数值模拟法中: 文献[3 - 7]均模拟了不同工况下矩形顶管施工引起的土体变形。实测数据分析法中: 温锁林[8]通过施工现场实测土体沉降、孔隙水压力和土体水平位移,得到施工对周边环境的影响规律; 施文捷等[9]通过现场实测数据,探讨了矩形顶管施工过程中周边水土压力及地表变形随时间的变化规律; 郭亮[10]通过施工现场监测土体变形,得到土体沉降的变化规律; 邓长茂等[11]通过对上海软土地层中3 个大截面矩形顶管施工实例分析,发现矩形顶管推进引起地表变形具有一般规律性。综上所述,目前对矩形顶管施工期间引起地表沉降的研究较多,但未见对工后地表沉降的研究,也不清楚Peck 公式是否适用于工后地表沉降计算。因此,有必要对矩形顶管引起的工后地表沉降作进一步研究。
本文借鉴盾构隧道工后地表沉降的计算方法,提出矩形顶管隧道施工引起的周围土体初始超孔隙水压力计算方法,采用分层总和法计算每层土体中由于超孔隙水压力消散引起的固结沉降量,叠加施工阶段的地表沉降量,得到工后总的地表沉降量。通过算例分析,研究工后地表横向和纵向沉降以及地表沉降速率的变化规律。
1 土体初始超孔隙水压力计算方法
1. 1 研究思路及假定
本文研究对象为软土地区矩形顶管施工引起的工后地表沉降( 包括施工期间沉降和固结沉降) 。在研究土体初始超孔隙水压力前,作以下假设: 1) 本文中出现的初始超孔隙水压力为土体中最大超孔隙水压力; 2) 土体开挖应力释放导致初始超孔隙水压力的产生,且隧道周围土体的应力释放率大小相等; 3) 矩形顶管隧道一侧视为挡土墙。
本文借鉴盾构隧道施工产生的土体超孔隙水压力计算方法: 通过计算出隧道水平线处起拱点( 图1 中A1点) 土体超孔隙水压力大小以及该点围压大小,二者相除得到应力释放率; 再根据应力释放理论,得出隧道周围土体初始超孔隙水压力大小及分布规律。
由于矩形和圆形在几何性质上存在不同,起拱点位置不一致。在矩形顶管中可假设矩形隧道一侧为挡土墙; 根据库仑土压力理论,图1 中A2点( 即墙趾点)为挡土墙滑动破坏面起拱点,因此选择A2点作为矩形顶管隧道初始超孔隙水压力的计算点( 以下简称隧道角点处) 。
1. 2 矩形顶管隧道周围土体初始超孔隙水压力计算
1. 2. 1 隧道围压的计算方法
矩形顶管隧道的围压受力模式见图2。
图1 盾构隧道与矩形顶管隧道初始超孔隙水压力计算点
图2 矩形顶管隧道受力示意图[12]
式( 4) —( 8) 中: U0 '为隧道角点处土体的初始超孔隙水压力; Δσz为平面应变条件下第二主应力的变化量; Δσr与Δσθ分别为径向和切向应力的变化量; Ps为隧道对土体的支护应力; ε、λ 为Henkel 孔隙水压力系数,ε = ( 3A - 1) /槡2 ,且饱和土中A = 0. 5,λ = 1; σ0 '为土体初始有效应力; u0为静止孔隙水压力; G 为土体Lame 弹性系数,G = E /[2( 1 + μ) ],E 为土体的弹性模量,μ 为土体泊松比; ηCF = 12 ( R' /tc) 2,其中R'定义为隧道的等效半径,即以隧道几何中心为圆心,根据矩形隧道超孔隙水压力计算点的位置作出半径为R'的外接圆,A2点的等效半径等于隧道外部尺寸对角线长度的一半; tc为隧道矩形顶管管节的厚度;D0 = Ec tc /( 1 - μ2c) ,Ec为隧道混凝土的弹性模量,μc为隧道混凝土的泊松比。
1. 2. 3 土体应力释放率计算方法
由式( 3) 和式( 4) 计算得到隧道角点( 图1 中A2点) 处的土体初始超孔隙水压力值( U0 ') 和水土压力值( σ3) ,则土体应力释放率α = U0 ' /σ3。
1. 2. 4 隧道周围土体初始超孔隙水压力计算
矩形顶管隧道周围的土体初始超孔隙水压力,可以通过隧道围压乘以土体应力释放率α 得到。由于矩形隧道存在转角,周围土体中隧道侧面超孔隙水压力与顶部及底部超孔隙水压力区域中间存在4 个过渡区域,这在圆形隧道中是没有的。矩形顶管隧道周边土体初始超孔隙水压力示意图见图3。
图3 矩形顶管隧道周边土体初始超孔隙水压力示意图
本文提出过渡区域可假定由一组对数螺旋线ρ = r0exp( a0θ') 组成,其中ρ 和θ'为对数螺旋线函数的变量,分别代表超孔隙水压力大小和该位置处与水平方向的夹角; r0、a0为对数螺旋线参数,可以通过曲线两端已知的初始超孔隙水压力大小联立计算得到。由此可计算得到4 个过渡区域( 图3 中虚线区域) 中土体的初始超孔隙水压力。
1. 3 土体扰动范围确定
本文在文献[14 - 15]的研究基础上提出了土体初始超孔隙水压力分布范围的具体计算方法。以隧道右半横截面为例,矩形顶管施工扰动范围如图4 所示,本文参考文献[15]的计算方法,可确定A、B、C 点的具体位置。
图4 矩形顶管隧道施工扰动范围
式( 9) —( 11) 中: a,b 为简便公式而设立的参数,a =H +h/2 +r',b = r'sin ( 45° +/2) +( H +h/2) cot ( 45° +/2)。
1. 4 分布范围内任一点土体初始超孔隙水压力计算
本文假定: 1) 由于矩形顶管开挖引起土体卸荷扰动,适用朗肯主动土压力理论,因此土体超孔隙水压力扩散面( 边界) 与水平面夹角为β = 45° + /2。超孔隙水压力由隧道开挖卸荷引起,而不是地表加载引起,因此超孔隙水压力是从隧道边界向外侧传递并扩散。2)由于沿隧道顶进方向距离较长,因此简化为平面应变问题,该方向在计算中长度均取1 m。3) 沿隧道边缘任一点垂直方向,土体内各点的初始超孔隙水压力值从隧道边界向外侧发生衰减。4) 土体超孔隙水压力扩散面上应力均匀分布。5) 周边土体各点的初始超孔隙水压力不受其余处影响,仅由该点相对应的与隧道相邻处初始超孔隙水压力U0传递。
以隧道横截面为例,本文超孔隙水压力传递模型见图5。
图5 隧道周围土体超孔隙水压力向四周传递示意图
1) 隧道拱顶上方土体初始超孔隙水压力计算。令隧道顶部向上传递应力的土体宽度为L。魏纲[17]在盾构隧道超孔隙水压力计算中提出令
( H0 + L) /( H0 +2R) = 0. 7,故L = 1. 4R - 0. 3H0。( 12)
式中: H0为盾构隧道的覆土埋深; R 为盾构隧道外半径,在矩形顶管隧道中,本文提出将等效半径R'( 1. 2. 2节中已给出取值方法) 代入式中计算更合适。
隧道顶部土体的初始超孔隙水压力值U0 =αPθ = 90°,其中Pθ = 90°为隧道顶部围压值。由竖向受力平衡推导得出侧面应力
式中d 为隧道到扰动边界的距离。
3) 隧道侧向土体超孔隙水压力计算。以隧道水平轴线上的侧向土体为例,推导得到隧道侧向距离k处的土体初始超孔隙水压力
式中j 为隧道到扰动边界的距离。
2 工后地表长期沉降计算
矩形顶管隧道施工引起的最终总的地表沉降量S由施工期间地表沉降S1和最终地表固结沉降量S2组成,则
图6 工后地表沉降计算模型简图
2. 1 施工引起地表沉降计算
林强强[1]证明Peck 公式[18]可用于计算矩形顶管施工期间的地表沉降S1,即:
式中: Vloss为隧道单位长度土体损失量,Vloss = hlη,其中η 为土体损失率; i 为施工阶段的地表沉降槽宽度系数; y 为计算点距隧道轴线的横向水平距离。
2. 2 施工结束后地表最终固结沉降计算
采用分层总和法单向压缩基本公式,可计算得到矩形顶管施工引起的地表最终固结沉降量S2。为简化计算,作出以下假设: 1) 由于各层土的附加应力不易于确定,假定土质均匀; 2) 土质参数( 压缩模量) 采用加权平均参数。
式中: Δp 为各层超孔隙水压力平均值; E—s为隧道上覆土层的压缩模量的加权平均值。
2. 3 地表沉降随时间的变化规律
土体固结t 时间之后的地表总沉降量S( y,t) 可表示为
式中: S2( y,t ) 为土体固结t 时间之后的固结沉降量;Uz为固结度,因为土体初始超孔隙水压力近似为三角形分布,可根据文献[19]中图2 - 52 中曲线( 3) 得到,其中土体固结时间因数( 竖向)
式中: H— 为y 处压缩土层的计算厚度; cv为土体固结系数( 竖向) ; t 为土体固结的时间。
通过研究沉降速率和时间之间的关系,可以得出单位时间内地表沉降量的变化。由不同时间t 所取的Uz计算得到地表沉降量S( y,t) 。令第i 个计算点的时间为ti,i从1 开始取值。则ti时间对应的地表沉降量为S( y,ti) ,t i + 1时间的地表沉降量为S( y,ti + 1) 。令( ti + 1 - ti) 时间段内的地表沉降量为ΔS( y,t i) ,则地表日沉降速率:
3 算例分析
3. 1 工程概况及计算过程
以上海市轨道交通六号线浦电路车站3 号出入口矩形顶管工程[1]作为实际案例,进行算例分析。顶管隧道长42 m,共28 个管节,每个管节长1. 5 m,内部尺寸为高3. 36 m、宽5. 24 m、厚度tc = 0. 5 m,即h = 4. 36m、l = 6. 24 m。覆土埋深H = 7. 2 m,地下水位0. 56m。土体各参数取值见表1。
表1 土体参数取值
地表沉降槽宽度系数i 的取值,参考文献[20]中的相关研究: i = k( H + 0. 5 h) ,黏性土中k 值分布范围为0. 37 ~ 0. 66,本文取k = 0. 37,得到i = 3. 47 m。土体初始有效应力可通过总应力减去土体中静止孔隙水压力得到。在该算例中,隧道角点处土体初始有效应力为110 kPa。
根据式( 4) —( 8) 可计算得到隧道角点处的土体初始超孔隙水压力为33. 01 kPa,由式( 3) 计算得到该点处隧道围压为106. 05 kPa,二者相除得到应力释放率α = 31. 12%。隧道顶部围压为124. 63 kPa,即可得到该点处初始超孔隙水压力为38. 79 kPa。
根据1. 4 节的计算方法,计算得到隧道周围土体扰动范围内的土体初始超孔隙水压力等值线图,见图7。根据分层总和法单向压缩基本公式,对地表横向最终固结沉降量进行计算。如图7 所示,土体共被划分为8 层,每层层厚为2 m,其中底层厚度为1. 95 m。y分别取0、1、2、4、6、10、14 m。
图7 土体初始超孔隙水压力等值线图( 单位: kPa)
根据式( 19) 可以计算得到横向地表最终固结沉降值,再利用式( 20) 可以计算得到地表总沉降量随时间的变化值,最后由式( 22) 可以计算得到隧道轴线上方地表总的日沉降速率随固结时间变化值。
3. 2 计算结果分析
图8 为隧道轴线上方地表工后总沉降计算值与实测数据[1]对比图。由于顶管全部贯通后实测工作结束,导致该断面工后仅实测25 d,没有后期数据。如图8 所示,采用本文方法所得计算结果与实测数据比较吻合且趋势一致,验证了本文方法的可靠性。
图8 隧道轴线上方地表总沉降计算值与实测数据对比
图9 为隧道轴线上方地表总沉降计算值随时间变化曲线。由图9 可知,初始阶段地表沉降增长速度较快,地表沉降随时间成凸曲线增长; 后期地表沉降随时间开始近似成线性增长,地表沉降速率逐渐减小至趋近不变。
图10 为横向地表总沉降随时间变化曲线。由图10 可知,随着时间增长,横向地表总沉降逐渐增大,但沉降曲线大致符合正态分布规律。
图11 为由式( 22) 计算得到的隧道轴线上方地表总的日沉降速率随固结时间变化曲线。由图11 可知,工后前3 个月地表沉降速率最大,随后半年内沉降速率的变化曲率迅速减小,并在接下来的数年里沉降速率保持稳定。
图9 隧道轴线上方地表总沉降计算值随时间变化曲线
图10 横向地表总沉降随时间变化曲线
图11 隧道轴线上方地表计算日沉降速率曲线
璩继立[21]的研究结果表明,虽然盾构施工引起的地面沉降槽宽度会随时间增长逐渐增大,但仍符合Peck 提出的Gauss 曲线分布规律。图12 为本文方法计算得到的地表最终固结沉降计算值与Peck 公式预测值比较。由图12 可知,二者变化规律较为吻合,表明矩形顶管的地表固结沉降也大致符合Gauss 曲线分布规律。
图12 地表最终固结沉降计算值与Peck 公式预测值比较
根据式( 18) 计算得到施工期间地表最大沉降量S1 = 20 mm,根据式( 19) 计算得到施工结束后地表最终最大固结沉降量S2 = 35. 45 mm,二者叠加得到工后总的最终最大地表沉降量S = 55. 45 mm。工后固结沉降占总沉降的63. 9%,不容忽视。但目前国内工程界对工后沉降还不够重视,仅在施工期间进行地表沉降监测,工后沉降几乎不考虑。本工程案例也是在顶管全部贯通后即停止实测,实测最大沉降量为27. 57mm,根据本文计算该断面后续还会产生27. 88 mm 的沉降量,这显然是不能忽略的。因此,建议在施工前先对工后沉降进行预测,根据计算结果来制定施工方案,同时延长监测时间,对工后监测频率可以适当降低。
4 结论与讨论
1) 提出矩形顶管隧道施工引起的周围土体初始超孔隙水压力计算方法,进而得到土体初始超孔隙水压力等值线图。运用分层总和法计算地表固结沉降量,叠加施工阶段的地表沉降量,得到工后地表总沉降量的计算方法,进而预估最终沉降量。
2) 本文方法计算得到的工后横向地表固结沉降计算值,其变化规律与Peck 公式预测值较为吻合; 横向地表总沉降值随着时间增长而逐渐增大,曲线基本符合正态分布规律。
3) 工后地表沉降主要发生在施工结束后半年内,沉降速率在施工结束后的3 个月内最大,随后半年内沉降速率的变化曲率迅速减小,并在较长时间沉降速率保持稳定。计算工后固结沉降占总沉降的60% 以上,不容忽视。
本文研究中作了较多假定,仅考虑土体开挖应力释放引起的超孔隙水压力,没有考虑其他因素; 土体扰动范围和超孔隙水压力的分布需作更深入的研究。另外需收集更多地表沉降实测数据,尤其是工后沉降数据,来进一步验证本文方法的可靠性。
摘自:隧道建设