Calculation of ultimate supporting forces of shield tunnels in unsaturated sandy soils considering soil arching effects
-
摘要: 因施工中受注浆及掘进参数影响,地层损失引起的土拱效应逐渐发挥,滑面上剪应力逐渐增大至抗剪强度,同时应考虑表观黏聚力对剪应力的贡献作用,而Terzaghi公式并不能体现这些因素。针对该问题,基于主应力偏转理论,推导了非饱和砂土不同位置的侧压力系数公式,结合表观黏聚力公式和试验数据,得到与基质吸力相关的分段形式的表观黏聚力表达式,进而推导出不同土拱效应发挥程度下的松动土压力及盾构极限支护力的计算公式。算例结果表明同一埋深下,不同饱和度土层中拱顶处松动土压力随土拱效应发挥而逐渐减小且均远小于自重应力。随着含水率的增大松动土压力存在一个阈值;不同土拱效应发挥程度下,埋深越大,不同饱和度地层拱顶处松动土压力趋向不同的定值;盾构极限支护力公式计算结果亦存在同样规律,该结论将为非饱和砂土层盾构施工过程中土仓压力的合理设置提供理论指导,保证施工安全。Abstract: Due to the influences of grouting parameters and tunnelling parameters in the actual construction, the soil arching effects caused by the ground loss gradually in-crease, and the shear stress on the sliding surface increases to the shear strength gradually. The contribution of the apparent cohesion of unsaturated sandy soils to the shear stress should be also considered. However, the above factors are not considered in Terzaghi's formula. In order to solve the problem, based on the principal stress rotation theory, the expressions for the lateral pressure coefficient of unsaturated sandy soils at different positions are deduced. Based on with the theoretical formula for the apparent cohesion and the experimental data, a piecewise function for the apparent cohesion related to the matric suction is obtained, and the formula for loosening soil pressure as well as the formula for ultimate supporting force is derived in consideration of different soil arching effects. The results of the example show that with the soil arching effects, at the same depth, the loosening soil pressure decreases gradually and it is always less than the self-weight stress. With the increase of water content, there is a threshold value for the loosing earth pressure. Under different soil arching effects, with the increase of the depth, the loose earth pressure at the arch crown of strata with different saturations will tend to different constant values. Meanwhile, the formula for the limit supporting forces of shield tunnels has the same law. This research results will provide theoretical guidance for controlling reasonable chamber pressure of shield tunnelling in the unsaturated sandy soils, and ensure the safety of construction.
-
0. 引言
近年来的地震灾害表明,与地上结构的地震惯性力破坏不同,地下结构的破坏主要取决于场地土的变形,场地条件比地震动强度对管道破坏的影响更大[1-3]。长输埋地管道在敷设过程中不可避免要穿越多变的场地,如沉积河谷场地、开挖回填场地以及地质运动直接形成的多种土介质场地等[4]。灾后现场调查表明[5-8],相比于均匀场地,埋设在非均匀场地中管道的破坏率明显更高。因此,研究管-土结构体系在非均匀场地中的地震反应,探讨非均匀场地中埋地管道的抗震性能,具有重要的学术意义和工程价值。
不同的地震动入射方向对长线型地下结构的地震反应影响很大。早期Hindy等[9]利用管道的集中质量模型研究了地震动沿管道横向、纵向以及不同入射角下管道的响应。结果显示管道中的轴向应力远大于弯曲应力。因此,纵向地震反应通常被认为是更关键的,埋地管道的地震响应分析和设计也集中在对结构轴向应变的评估,忽略了由横向地震动作用引起的结构弯曲应变。然而,Ariman等[10]的研究表明大直径管道存在显著的弯曲应变,同时指出对大直径管道的地震响应分析仅考虑纵向地震动作用的影响,结果可能是不适用的。Zerva等[11]研究了随机地震作用下管道的轴向和水平横向的响应。尽管也观察到管道的轴向应变超过了弯曲应变,但同时也指出对大直径管道,应考虑纵向、横向地震作用下管道的地震响应。闫孔明等[12]开展了非一致激励下埋地管道的振动台试验,将弯曲应变结果与国内和美国管道规范结果对比,表明规范设计时弱化甚至忽略了管线弯曲应变,但是这种做法难以解释管线弯折破坏。虽然许多研究人员考虑了地震动的输入方向问题,但对于纵向和横向地震动作用下管道地震响应所占的主导地位仍存在争议,尚未达成共识。且目前的研究多为经验关系、理论和数值探索性分析,试验研究相对较少。
近年来,国内外学者相继开展了埋地管道振动台模型试验的开拓性研究。Kawabata等[13]对埋地柔性管道进行了振动台试验,研究了在环刚度相似的情况下,厚度对柔性管道动力特性的影响,利用不同的管道基础考虑场地的变化,但是不同场地是分开进行的,无法研究场地变化处管道的动力特性。Jafarzadeh等[14]利用振动台试验研究了不同的加速度幅值、地震动频率、土壤密度,管道基础条件和振动方向,对加速度放大系数和管道变形的影响。该试验通过在部分箱体底部放入混凝土基础模拟非均匀场地,采用简谐波作为输入地震动,未考虑其非一致激励的影响。
由于地震动传播速度有限,不可能同时到达长输埋地管道上所有位置,考虑行波效应的非一致激励是更合理的地震动输入方法。Newmark [15]提出地下管道在地震波作用下的分析方法,考虑行波效应,略去惯性力的影响和管土间相对位移。Wang等[16]不考虑惯性力和阻尼影响,提出拟静力方法,指出地震波相位差对地下管道的影响较大。此外,已开展的埋地管道研究表明非一致地震激励引起的管道沿线不同的土体运动,可能是造成埋地管道失效破坏的重要原因[17-18]。鉴于此,国内学者开展了非一致地震激励下埋地管道的振动台试验。孟海[19]首次采用两个独立剪切模型箱,开展了埋地管道纵向、横向非一致激励振动台试验,对管-土结构体系的加速度、位移、应变等结果进行分析,结果表明,横向地震动作用下,管道产生了纵向弯曲变形,应变反应比纵向地震动作用下更大,管-土之间的滑移更小。但试验中两独立模型箱中间无覆盖土层,模型土体不连续,造成了中间部分的管道裸露。Yan等[20]采用独立控制的双台面振动台,开展了埋地管道三向非一致激励双台阵振动台试验研究,探讨了管道的加速度、应变以及土体加速度等参数的变化规律,模型土体亦存在不连续的情况。上述虽考虑了非一致激励的影响,但均未考虑场地非均匀性的影响。
本试验在充分吸取现有振动台试验和理论研究成果对埋地管道结构地震反应认识的基础上,考虑场地条件、模型结构材料及输入地震动特性等,针对埋设在非均匀场地中的管道开展了横向一致、非一致激励振动台模型试验。本文对不同场地条件、不同激励方式、不同地震强度作用下埋地管道的加速度和应变等试验数据进行分析,研究了横向非一致激励下非均匀场地中埋地管道的地震响应特性及其变化规律。
1. 振动台模型试验介绍
1.1 试验模型设计
本试验采用自主研制的悬挂式连续体模型箱(SCMB),在北京工业大学振动台台阵系统上开展了埋地管道三台阵振动台试验,如图1所示。模型箱的制作考虑3个方面。
(1)结构坚固,连接可靠,在最大加载等级下模型箱可以稳定地进行试验。
(2)在考虑振动台系统的承载能力的情况下,模型箱尺寸可以保证场地土体的连续性和有利于模型管道的布置。
(3)模型箱与场地土之间的满足动力试验要求,模型箱效应不明显;有利于多点非一致激励加载。
对连续体模型箱空箱基频试验和自由场试验进行了有限元数值模拟,结果表明所设计连续体模型箱的动力特性和边界效应基本满足试验要求。SCMB的设计及适用性验证分析详见文献[21]。
试验中管道和土体的相似比基于Backingham定理[22]设计。相似关系的确定详见文献[23]。管道和土体的几何相似比Sl=1/10,弹性模量相似比SE分别为1/57,1/4,密度相似比分别为1/6.5,1,应变相似比Sε=1,频率相似比Sf=5,加速度相似比Sa=2.5。
管道原型为大直径管道,长60 m,直径1.5 m,厚20 mm,根据几何相似关系,模型管道长6.0 m,直径0.15 m,厚2 mm。为了更好地反映管道在地震动作用下的变形特性,便于测得管道的应变,采用模量较低的有机玻璃管。模型管道的敷设如图2所示。管道的材料参数如表1所示。
表 1 管道的材料参数Table 1. Material parameters of pipelines类型 材质 弹性模量/GPa 密度/(kg·m-3) 长度/m 管径/mm 壁厚/mm 原型 钢管 210 7800 60 1500 20 模型 有机玻璃管 3.2 1200 6 150 2 试验模型土取自北京地铁14号线北工大站10 m深的砂土和14 m深的粉质黏土。土体采用直径20 cm的铁锤夯实,将场地土划分为3个区域,均匀场地土体均为砂土,非均匀场地两端区域为砂土,中间区域为粉质黏土如图3(a)所示。
1.2 试验输入地震动及加载方案
为研究输入地震动特性对埋地管道结构地震反应的影响,本试验选择El-El-Centro地震动、Kobe地震动及BJ-北京人工合成地震动3条地震动记录作为输入地震动,其频率和峰值根据模型土的相似关系按5倍和2.5倍进行调整,如图4所示,为考虑地震动行波效应的影响,相邻子台的地震动输入延迟设置为1 s,在每级加载前进行WN-白噪声激振,获取场地的动力特性。振动台试验加载工况如表2所示。
表 2 振动台试验加载工况表Table 2. Loading events of shaking table tests加载等级 地震动 横向加速度峰值/g 原型 1/10模型 扫频 WN 0.05 0.05 一级加载 EL,Kobe,BJ 0.10 0.25 二级加载 EL,Kobe,BJ 0.20 0.50 三级加载 EL,Kobe,BJ 0.40 1.00 四级加载 EL,Kobe,BJ 0.62 1.55 1.3 试验传感器的布置
非均匀场地中传感器的布置如图3,5所示,均匀场地中的布置与之一致。图3(a)展示了布置在土体表面的加速度计(A)和激光位移传感器(D);图3(b)展示了模型管道的9个应变监测截面S1-1—S5-1;图3(c)展示了布置在#3,#4,#5箱体垂直中心线的X和Y向加速度计。图5(a)展示了#3,#4箱中管道环向布置的薄膜土压力计(MP);图5(b)展示了#2,#3,#4箱中管道环向布置的应变片(S);图5(c)展示了#1和#5箱中管道环向布置的应变片;图5(d)展示了布置在#2,#3,#4箱中管道上的微型加速度计(MA)。
2. 试验结果和分析
2.1 加速度时程曲线
图6,7展示了在横向一致和非一致激励下不同场地下埋地管道(MAY1)和周围土体(AY17)在监测断面2(#4箱)的加速度时程和傅氏谱。由图可知,埋地管道和及其周围土体在不同激励方式,不同场地条件下的加速度时程曲线都吻合较好。说明横向地震动作用下,管道主要受周围土体的约束,其加速度反应基本与周围土体相同。在一致激励作用下,管道和周围土体加速度傅式谱的频率分布较为一致,在较低加载等级下,非均匀场地中的管道和土体的加速度幅值比均匀场地大35%左右。在较高加载等级下,非均匀场地和均匀场地中管道和土体的加速度反应几乎相同。在非一致激励下,管道的加速度傅式谱的频带相比于管周土稍宽,低频分量放大幅度稍大,同一加载等级下,非均匀场地下管道和土体的加速度反应稍小。
2.2 管周应变时程曲线
图8展示了0.5g加载时,横向一致和非一致激励下均匀场地和非均匀场地中管道应变截面S4-2监测点0°和180°,90°和270°的应变时程曲线。结果表明,在所有工况中管道截面0°和180°测点的应变时程的相位和波形基本吻合,而管道截面90°和270°测点的应变时程相位相反。与管道截面0°和180°测点的应变响应相比,管道截面90°和270°测点的应变相对较大,非均匀场地中管道的应变比均匀场地中稍大。与横向一致激励相比,在加载等级0.5g非一致激励作用下管道的应变没有明显改变。说明在横向地震动0.5g作用下,均匀和非均匀场地中的管道被周围土体夹裹随其主要产生了纵向弯曲变形。
2.3 管道轴向应变峰值
图9,10展示了均匀和非均匀场地中,在横向一致、非一致地震动作用下,管道的纵向弯曲拉、压应变峰值(每个截面多个测点的拉伸、压缩应变峰值平均值)。
由图可知,管道拉、压应变峰值曲线在均匀和非均匀场地中沿管道轴向的变化趋势相似,中间截面S3-1、S4-1处管道拉伸、压缩应变峰值较大,随着加载等级的增加,管道拉、压应变峰值随之增大,且增大的速率逐渐变大。在一致激励作用下,在非均匀场地土质变化过渡区及粉质黏土中(管道截面应变S3-1、S4-1)的拉伸应变峰值比均匀场地中更大,在1.0g,1.55g加载等级时,与均匀场地相比,最大拉应变峰值高30%左右,但其最大压缩应变峰值比均匀土中略小。此外,在较高加载等级1.0g,1.55g时,均匀土中管道的应变响应非一致激励比一致激励大两倍左右。但非均匀土中管道的应变响应变化不明显,且与一致激励下不同,在非一致激励较高加载等级1.0g,1.55g时,最大拉应变峰值比均匀场地低30%左右,且出现在粉质黏土中管道应变截面S3-1处。说明在非一致激励和非均匀场地共同作用下管道的地震响应比单一影响因素作用下有所减小。
2.4 管道弯曲应变与加载等级的关系
图11,12是管道横截面S2-3,S3-1,S4-1纵向弯曲峰值应变随加载等级变化趋势图。
由图11,12可知,均匀和非均匀场地中埋地管道的纵向拉伸和压缩应变在一致、非一致激励作用下,随着加载等级的增大基本呈线性增加的趋势。对有机玻璃材料的拉伸、压缩材性试验显示,在±10000 με范围内其应力-应变大致呈线性关系。结合振动台试验结果可知,在一致、非一致激励作用下,动应力-应变基本呈线性关系,埋地管道刚度基本不变,也说明管道反应基本服从于周围土体的反应,与2.1节中管道加速度反应分析的结论相互印证。在一致激励作用下非均匀场地中的管道的拉伸应力-应变关系曲线的斜率相比于均匀场地明显增大,说明非均匀场地中埋地管道的拉伸应变对地震等级的增大更加敏感。
3. 结论
本文对横向非一致激励下砂土均匀场地和砂土、粉质黏土非均匀场地中的埋地管道进行了多台阵振动台试验研究,得到3点结论。
(1)埋地管道在横向一致、非一致激励下管道与周围土体的加速度反应基本一致,但在非一致激励下,埋地管道的加速度傅式谱频带稍宽,低频分量放大幅度稍大。一致激励下,非均匀场地中的管道和土体的加速度幅值比均匀场地大35%左右,非一致激励下,非均匀场地下管道和土体的加速度反应反而稍小。
(2)均匀场地和非均匀场地中,埋地管道在横向一致、非一致激励下随土体纵向弯曲变形,产生了弯曲应变,与均匀场地相比,非均匀场地中管道的应变明显增大,但在加载等级0.5g非一致激励作用下管道的应变没有明显改变,非一致激励对管道的应变响应影响较小。
(3)在较低加载等级下,地震动激励方式和场地条件变化对管道应变响应的影响较小,但随加载等级增加二者的影响较为显著,在只考虑非一致激励或非均匀场地相比于一致激励和均匀场地对管道响应的影响时,管道的应变都显著增大,而同时考虑非一致激励和非均匀场地时,管道的应变响应反而较只考虑单一因素时减小。
-
-
[1] BROMS B B, BENNERMARK H. Stability of clay at vertical openings[J]. Journal of the Soil Mechanics and Foundations Division, ASCE, 1967, 96(1): 71-94.
[2] LECA E, DORMIEUX L. Upper and lower bound solutions for the face stability of shallow circular tunnels in fric-tionalmaterial[J]. Géotechnique, 1990, 40(4): 581-606. doi: 10.1680/geot.1990.40.4.581
[3] LEE IM, NAM SW. The study of seepage forces actingon the tunnel lining and tunnel face in shallow tunnels[J]. Tunnelling and Underground Space Technology, 2001, 16(1): 31-40. doi: 10.1016/S0886-7798(01)00028-1
[4] 吕玺琳, 王浩然, 黄茂松. 盾构隧道开挖面稳定极限理论研究[J]. 岩土工程学报, 2011, 33(1): 57-62. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201101012.htm LÜ Xi-lin, WANG Hao-ran, HANG Mao-song. Limit theoretical study on face stability of shield tunnels[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(1): 57-62. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201101012.htm
[5] ANAGNOSTOU G, KOVARI K. The face stability of slurry -shield-driven tunnels[J]. Tunnelling and Underground Space Technology, 1994, 9(2): 165-174. doi: 10.1016/0886-7798(94)90028-0
[6] 吕玺琳, 李冯缔, 黄茂松, 等. 三维盾构隧道开挖面极限支护压力数值及理论解[J]. 同济大学学报(自然科学版), 2012, 40(10): 1469-1473. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201210005.htm LÜ Xi-lin, LI Feng-di, HUANG Mao-song, et al. Three-dimensional numerical and analytical solutions of limit support pressure at shield tunnel face[J]. Journal of Tongji University(Natural Science), 2012, 40(10): 1469-1473. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ201210005.htm
[7] 汤旅军, 陈仁朋, 尹鑫晟, 等. 密实砂土地层盾构隧道开挖面失稳离心模型试验研究[J]. 岩土工程学报, 2013, 35(10): 1830-1838. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201310010.htm TANG Lü-jun, CHEN Ren-peng, YIN Xin-sheng, et al. Centrifugal model tests on face stability ofshield tunnels in dense sand[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(10): 1830-1838. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201310010.htm
[8] TERZAGHI K. Theoretical Soil Mechanics[M]. New York: John Wiley and Sons, 1943: 37-42.
[9] 黎春林. 盾构隧道施工松动土压力计算方法研究[J]. 岩土工程学报, 2014, 36(9): 1714-1720. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201409024.htm LI Chun-lin. Method for calculating loosening earth pressure during construction of shield tunnels[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(9): 1714-1720. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201409024.htm
[10] 陈若曦, 朱斌, 陈云敏, 等. 基于主应力轴旋转理论的修正Terzaghi松动土压力[J]. 岩土力学, 2010, 31(5): 1402-1406. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201005013.htm CHEN Ruo-xi, ZHU Bin, CHEN Yun-min, et al. Modified Terzaghi loozening earth pressure based on theory of main stress axes rotation[J]. Rock and Soil Mechanics, 2010, 31(5): 1402-1406. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201005013.htm
[11] MARSTON A, ANDERSON A O. The theory of loads on pipes in ditches and tests of cement and clay drain tile and sewer pipe[R]. Ames: Iowa Engineering Experiment Station, Iowa State Collage, 1913.
[12] HANDY R L. The arch in soil arching[J]. Journal of Geotechnical Engineering, 1985, 111(3): 302-318.
[13] 蔺港, 孔令刚, 詹良通, 等. 基于太沙基土拱效应考虑基质吸力影响的松动土压力计算模型[J]. 岩土力学, 2015, 36(7): 2095-2104. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507045.htm LIN Gang, KONG Ling-gang, ZHAN Liang-tong, et al. An analytical model for loosening earth pressure considering matric suction based on Terzaghi soil arch effect[J]. Rock and Soil Mechanics, 2015, 36(7): 2095-2104. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201507045.htm
[14] BISHOP A W, ALPAN I, BLIGHT G E, et al. Factors controlling the shear strength of partly saturated cohesive soils[C]//ASCE Research Conference on Shear Strength of Cohesive Soils, 1960, Boulder: 503-532.
[15] FREDLUND D G, MORGENSTERN N R, WIDGER R A. The shear strength of unsaturated soils[J]. Canadian Geotechnical Journal, 1978, 15: 313-321.
[16] 林鸿州, 李广信, 于玉贞, 等. 基质吸力对非饱和土抗剪强度的影响[J]. 岩土力学, 2007, 28(9): 1931-1936. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200709030.htm LIN Hong-zhou, LI Guang-xin, YU Yu-zhen, et al. Influence of matric suction on shear strength behavior of unsaturated soils[J]. Rock and Soil Mechanics, 2007, 28(9): 1931-1936. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200709030.htm
[17] 董倩, 侯龙, 赵宝云. 基质吸力对非饱和粉质砂土抗剪强度的影响[J]. 中南大学学报, 2012, 43(10): 4017-4021. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201210041.htm DONG Qian, HOU Long, ZHAO Bao-yun. Influence of matric suction on shear strength of unsaturated silty sand[J]. Journal of Central South University (Science and Technology), 2012, 43(10): 4017-4021. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201210041.htm
[18] 加瑞. 盾构隧道垂直土压力松动效应的研究[D]. 南京: 河海大学, 2007. JIA Rui. Study on Relaxation Effect of Vertical Soil Pressure for Shield Tunnel[D]. Nanjing: Hohai University, 2007. (in Chinese)
[19] 秦建设. 盾构施工开挖面变形与破坏机理研究[D]. 南京: 河海大学, 2005. QIN Jian-she. Study on Face deformation and Collapse of Earth Pressure Shield Tunnel[D]. Nanjing: Hohai University, 2005. (in Chinese)
[20] LEE I M, LEE J S, NAM S W. Effect of seepage force ontunnel face stability reinforced with multistep pipegrouting[J]. Tunneling and Underground Space Technology, 2004, 19(6): 551-565.
[21] ANAGNOSTOU G, KOVARI K. The face stability of slurry-shield-driven tunnels[J]. Tunneling and Underground Space Technology, 1994, 9(2): 165-174.