Face stability of shield-driven tunnels using multi-tangent technique
-
摘要: 针对非线性条件下深埋隧道开挖面稳定性问题,采用多切线技术分段线性近似Power-Law(P-L)准则,基于极限分析上限定理,构建1种新的开挖面多曲线圆锥破坏机构,提出1种考虑拉应力截断Tension cut-off(T-C)效应的开挖面稳定性理论计算方法。通过与既有文献结果对比,验证所提方法的可靠性。根据反算的应力分布和破坏机构几何特点,探究不同参数条件下T-C对开挖面稳定性的影响。结果表明:①在线性条件下,黏聚力较大且内摩擦角较小时,T-C对开挖面临界支护力的影响较大;在非线性条件下,m的增加并不会放大T-C对开挖面稳定性的影响。②在抗拉强度较大时,T-C对开挖面稳定性影响较为显著,随着抗拉强度的折减,T-C控制的破坏机构尺寸逐渐减小,且顶部区域明显钝化。
-
关键词:
- 开挖面稳定性 /
- 多切线技术 /
- Power-Law准则 /
- 拉应力截断 /
- 极限分析
Abstract: To address the face stability problem of deep-buried tunnels under nonlinear conditions, using Power-Law (P-L) criterion the approximate piecewise and by multi-tangent technique, a novel multi-cone failure mechanism is developed for tunnel faces based on the kinematic limit analysis theorem. A theoretical framework considering the tension cut-off (T-C) is proposed to evaluate the stability of tunnel faces. The reliability of the proposed method is validated through comparison with that in the existing literatures. Based on the back-calculated stress distribution and geometric feature of the failure mechanism, the influences of T-C under different parameters on the stability of tunnel faces are investigated. The results indicate that: (1) under the linear condition, the effects of T-C on the critical support pressure are more significant with high cohesion and small internal friction angle. Under the nonlinear condition, an increase in the nonlinear coefficient will not exacerbate the impact of T-C on the face stability. (2) For higher tensile strength, the effects of T-C on the stability of tunnel faces are significant. As the tensile strength decreases, the failure mechanism governed by T-C gradually shrinks, and the top of the failure mechanism exhibits noticeable blunting and looks like a dome. -
0. 引言
随着环境温度周期性变化,季节性冻土在反复冻融循环过程中,内部孔隙水随之发生上下迁移、相变,土中孔隙数量和大小将随着孔隙水的冻结与融化而发生变化,从而引起土体内部微观结构发生改变,导致土体出现冻胀融沉。冻融循环是一种连续不断的强风化过程,不断破坏着土的结构性,是引起土体强度参数降低,产生劣化的重要原因[1]。为克服冻融循环对土体性质的劣化影响,工程中通常采用不同的方法对公路地基土进行改良。改良方式主要可以分为物理、化学和综合改良3种[2]。物理改良法主要是在土体中加入土工织物[3]以及聚合物纤维等[4];或是通过掺加炉渣、砂、碎石等粗颗粒来改变其颗粒级配以增强土体的物理、强度特性[5]。化学改良方法目前采用较多,在土体中外加水泥[6]、石灰[7]或是化学浆液[8]等,通过改良材料与土颗粒间发生化学反应,从而生成强度较高的新物质,以提高土体的强度和稳定性。综合法则是物理改良与化学改良的结合[9]。
可以发现,上述国内外学者关于冻融循环对改良土的研究虽多,但以往改良材料的选取大多数局限于水泥、石灰、聚丙烯高分子纤维等,并不全面,有待进一步去探索新型改良材料。众所周知,石灰、水泥通常表现出以拉伸开裂为主的脆性破坏模式;此外这种碱性化学材料的掺入,会对周围土壤环境造成严重的污染。随着国家的迅速发展,环境污染日益严重,环保问题已成为阻碍国家发展的重要因素之一,绿色、经济改良土的物理力学性质研究将成为国内外岩土工程领域的一个崭新课题。而木质素纤维作为造纸工业副产品之一,来源广泛、价格低廉,具有经济、环保以及良好耐久性的优点,因此岩土工程领域中木质素纤维改良土具有很好的应用前景。
此外研究表明,土的宏观力学性质从根本上受其微观结构的控制,通过对土微观结构的研究,可以从微观的角度上解释土的力学性质[10]。刘春等[11]基于PACS软件进行裂纹模式分析。PACS软件具有自动分割识别颗粒和孔隙的功能,可以有效地减少人为因素对处理结果的影响,因此本文选用PACS软件对SEM图片进行处理。
本文设计完成了多种掺量的木质素纤维改良土在不同冻融次数后的UU三轴试验,着重分析了不同冻融次数下改良土的应力应变特性、弹性模量、强度指标等力学性质的变化规律。探讨了改良土的抵抗冻融效果并得出最佳掺量。并对不同冻融次数后的改良土进行SEM电镜扫描试验,从微观孔隙、概率熵变化的角度解释冻融对土体宏观力学性质影响的机理。
1. 试验材料及试验方案
1.1 试验材料
试验所用土样为天津滨海新区软黏土,在105℃干燥箱中烘干24 h后,碾散并筛分,按《土工试验方法标准》(GB/T 50123—2019)测得土样基本物理指标见表1。
表 1 土样物理特性指标Table 1. Physical properties of soil sampleswop/% γ/(kN·m-3) Gs wL /% wP /% IP 19.2 18.7 2.67 42.4 22.5 20.4 试验用木质素纤维来自河北石家庄石粉厂。纤维呈絮状,长度1~3 mm,在水中具有良好的分散性。
1.2 土样制备
首先将土样晒干,然后粉碎并过筛,以除去大颗粒杂质。进行三轴试验的关键是确保制备的试样具有良好的一致性,为使纤维与土样充分均匀混合,本文试验的制样方法是:先称取一定质量干土,将不同配比(分别为干土质量的0%,0.25%,0.5%,0.75%,1%)的纤维加入一定量的水中(干土质量的20%)混合,再喷洒入干土中,充分搅拌30 min后,套上保鲜膜密封静置12 h待水分均匀分布后,再将土样分3层填入制样器,压实后削平土样表面。试样最终尺寸:ϕ39.1 mm×80 mm。
1.3 试验方案
冻融试验箱可控温度范围为-25℃~150℃,如图1(a)。将上述制得的土样用保鲜膜密封,形成封闭的冻融环境,防止冻融循环试验中水分的散失和外界补给。根据文献[1, 3, 4]可知,5次冻融后,冻融对土体的物理力学性质影响基本稳定。故本文的冻融次数分别选取0,1,3,5,7,9次,试样冻结、融化温度分别设为-20℃、20℃,冻结、融化时间均为12 h,以保证土样内部完全冻结、融化。
由于软黏土透水性较差,故本文对其进行不固结不排水三轴试验(UU),试验采用双联三轴应变控制式系统,如图1(b)所示,试验机轴向加载量程选择0~1 kN,三轴围压、反压压力范围0~1 MPa。对经过不同冻融次数后的土样进行UU三轴试验,围压分别选取100,200,300 kPa,剪切速率选取0.5 mm/min,剪切应变达到20%时停止试验。
微观试验采用型号为Quanta250的电镜扫描仪,如图1(c)。电子发射电压为5 kV,接口连接到计算机,可以将扫描图像直接传送到计算机。分别对不同冻融次数后木质素纤维掺量为0%,0.25%,0.75%的改良土进行电镜扫描,扫描倍数设为100,500,1000。
2. 试验结果与分析
2.1 应力-应变曲线
图2为0.75%掺量木质素纤维改良土的应力应变曲线。可以看出,在100 kPa围压下,改良土的应力应变曲线均表现为应变软化型,而在200,300 kPa的围压下,表现为应变硬化型。这体现了围压大小对冻融循环后土体应力应变曲线具有重要的影响。在较高围压下,由于侧限的存在,随机分布在土体中的纤维构成了三维网架结构,加强了土中潜在的薄弱面,减小土中裂隙的延伸和发展,提高土样塑性。同时高围压还会增强了土颗粒间的联锁力,增加了相对滑移面(即剪切面)上的正应力,从而减小了冻融过程中产生的裂缝和裂隙,围压的压实作用增强了土体的变形抗力。
此外未经过冻融的土体体现出更强的应变软化特性,随着冻融次数的增加逐渐由应变软化型向硬化型转变,应力应变曲线的平面位置随冻融循环次数的增加而不断下降。这主要是冻融对土体有一个结构重塑效应,土体在冻结时,内部水结冰形成的冰晶体在土体内部产生冻胀力,当冻胀力超过土颗粒间的黏结力时,增大土颗粒间的微孔隙,导致土颗粒结构发生破坏,重新排列,从而在宏观上应力幅值随冻融循环次数的增加表现出逐渐劣化衰减的趋势。
图3绘制了不同木质素纤维掺量改良土在100,300 kPa两种围压下的应力应变曲线。在100 kPa下,随着木质素纤维掺量的增加,应力应变曲线逐渐呈现出应变软化特性。这是由于絮状的木质素纤维在加入水中时,生成的浆体有一定的胶结固化特性,改良土产生了一定的脆性,表现出应变软化的特点。而在300 kPa围压下,由于高围压对土体的塑性提升大于胶结固化产生的脆性,所以表现出更强的应变硬化特性。
2.2 弹性模量
土体弹性模量的确定通常是取1%轴向应变时偏应力增量与轴向应变增量之比[12]:
E=ΔσΔε=σ1.0%−σ0ε1.0%−ε0, (1) 式中 E为弹性模量;
Δσ 为应力增量;Δε 为应变增量,σ0 和ε0 分别对应初始应力和初始应变。根据式(1),计算得到不同木质素纤维掺量在不同冻融次数后弹性模量的变化曲线,由于篇幅有限,仅列出100 kPa围压下的结果,如图4。冻融循环会使弹性模量大幅度减小,第1次冻融循环后的弹性模量降低幅度较大。而木质素纤维的掺入使弹性模量产生了不同程度的提高,其中,木质素纤维掺量为0.75%时改良土的弹性模量最大。
为探究在同一冻融次数下不同木质素纤维掺量抵抗冻融劣化的效果,定义弹性模量的衰减幅度
DE 为DE=EN−E0E0×100%, (2) 式中,
EN 为N次冻融后的弹性模量,E0为未冻融的弹性模量。经过上式计算,冻融9次后,不同木质素纤维掺量改良土弹性模量的衰减幅度见表2。在200,300 kPa围压下,0.75%掺量的改良土弹性模量衰减幅度分别为53.55%,47.33%,均为同一围压下的最低值。在100 kPa围压时,衰减幅度为58.74%,仅比0.25%掺量的改良土大,因此综合来看,可以说明0.75%木质素纤维改良土抵抗冻融对弹模劣化影响的效果最好。
表 2 改良土9次冻融后的弹性模量衰减幅度Table 2. Attenuation amplitudes of elastic modulus of reinforced soils after 9 freeze-thaw cycles(%) 纤维掺量/% 围压 100 kPa 200 kPa 300 kPa 0 -58.96 -54.38 -62.47 0.25 -52.13 -56.02 -58.48 0.50 -64.86 -60.64 -68.72 0.75 -58.74 -53.55 -47.33 1.00 -59.59 -56.73 -47.69 2.3 破坏强度
对于应变硬化型,取15%轴向应变对应的偏应力作为其破坏强度,对于应变软化型,取其峰值应力作为破坏强度。据此可以得到围压100 kPa下木质素纤维改良土的破坏强度,如图5所示,将改良土经9次冻融循环后的破坏强度衰减值列于表3。
表 3 改良土9次冻融后的强度衰减幅度Table 3. Attenuation amplitudes of strength of reinforced soils after 9 freeze-thaw cycles(%) 纤维掺量/% 围压 100 kPa 200 kPa 300 kPa 0 -34.90 -26.32 -33.02 0.25 -34.35 -21.67 -32.26 0.50 -25.50 -22.83 -22.96 0.75 -26.51 -16.92 -17.16 1.00 -26.33 -17.64 -19.82 由图5可以看出,在同一冻融次数下,当木质素纤维掺量小于0.75%时,强度随掺量的增加而增大,超过0.75%时,强度随掺量的增加反而大幅度减小,呈现出这种现象的原因主要是木质素纤维掺量较多,分布比较集中,在土样内部形成了薄弱面,从而降低了土样强度。0.5%和0.75%掺量的木质素纤维改良土性能最优,二者强度接近,而超过5次冻融后,0.75%掺量的改良土在强度方面的优势开始扩大。
另外,由表3可知,随着冻融次数的增加,土样强度均出现不同幅度的降低。9次冻融后,素土的强度衰减幅度在100 kPa围压下达到最大值34.9%,而0.75%木质素纤维改良土的衰减幅度在200,300 kPa围压下分别为16.92%,17.16%,均为同围压下的最小值,故综合来看,可以确定木质素纤维最佳掺量在0.75%附近。
2.4 抗剪强度指标
根据UU试验得到的应力应变曲线,对3种围压所形成的破坏包络线进行直线拟合,得到抗剪强度指标随冻融循环的变化规律,如图6所示。整体来看,改良土表现出更高的黏聚力和内摩擦角。素土的黏聚力随着冻融次数的增加,呈现出先增加后降低的趋势,内摩擦角随冻融次数的增加而减小7°,变化幅度较大。0.75%木质素纤维改良土的黏聚力和内摩擦角均为最大;且0.75%掺量的改良土随冻融次数的增加,抗剪强度指标变化幅度较小、相对更加稳定。因此也可以说明木质素纤维最佳掺量在0.75%附近。
3. 微观结构分析
100倍SEM图像如图7所示,其中黑色代表孔隙,白色代表土体。可以看出9次冻融后的孔隙明显大于1次冻融后的孔隙,且9次冻融后颗粒有明显的抱团聚拢现象。0.75%木质素纤维改良土的裂隙明显比素土更少。对改良土样局部放大1000倍后得到图8,可以发现纤维在土样内部主要起到了三维网架结构或是“桥梁”搭接作用,加强了土体之间的联结力,从而在宏观上增加了力学性能指标。
通过PACS图像处理软件对电镜扫描后的图像进行定量分析,可以得到孔隙度、概率熵等参数。
对100倍下的图像进行二值化处理,并降噪分割,可以得到不同冻融次数后的孔隙大小。如图9所示,随着冻融次数的增加,土样内部孔隙率明显增大,在宏观上体现出冻融对力学性质的劣化影响。另外,0.75%木质素纤维改良土相比于其他土样有更小的孔隙度,从而提高了土体的破坏强度,说明了改良效果最佳,这也与前文分析得到的数据结果一致。
概率熵是一个描述颗粒排列的结构参数,可用来分析冻融循环作用后的土颗粒排列情况。从图10中可以看出,随着冻融次数的增加,三种土样的概率熵都逐渐减小。这说明冻融作用对土颗粒排列的定向性有一定的影响,根据概率熵的定义可知,在冻融作用下,颗粒排列由较混乱状态向有序状态发展,也就是说颗粒排列的定向性越来越好。
这是因为,在冻融循环过程中,外界温度条件的改变使得孔隙水发生相变,形成一定体积的冰晶结构,相邻土颗粒在冰晶体的楔形力挤压作用下,颗粒发生相对位移。此外在低温冻结过程中,由于冻结锋面的存在,土体内部水分有向冻结锋面定向迁移的趋势,为土颗粒的位移提供条件,故而在冻融循环作用下,土颗粒间的相对位移改变了土颗粒的排列方式,使土颗粒排列的定向性越来越好。
此外,还可以发现,0.75%纤维改良土的概率熵要比另外两种土样的小。说明了0.75%纤维改良土中的土颗粒排列在经历冻融循环后更加有序,抵抗冻融的劣化影响更加明显。
4. 结论
本文以木质素纤维改良土为研究对象,开展了一系列冻融循环试验、UU三轴试验和SEM电镜扫描试验,分析总结了冻融次数和木质素纤维掺量对该改良土力学性质变化的影响规律,探讨了木质素纤维改良土的抵抗冻融效果以及微观机制。主要得到以下结论:
(1)随着围压的增大,冻融次数的增加,应力应变曲线逐渐由应变软化型过渡到应变硬化型。随着木质素纤维掺量的增加,应力应变曲线呈现出软化特性;高围压下,此现象不明显。
(2)首次冻融后,弹性模量降低幅度较大,之后降低幅度逐渐放缓。木质素纤维掺量为0.75%时,改良土的弹性模量、破坏强度和黏聚力均达到最大值,内摩擦角变化最为稳定,抵抗冻融劣化能力最强。在实际工程应用中,建议采用0.75%的纤维掺量来改良季节性冻土地区路基土。
(3)从微观图像上可以看出,冻融作用使土体内部微观孔隙增大,导致宏观上力学性能的衰减,概率熵随着冻融的增加而减小。纤维在改良土样内部主要起到了三维网架结构或是“桥梁”搭接作用,减小了冻融对土体孔隙损伤的影响,从而增强了土体的强度。
-
表 1 临界支护力与理论解对比
Table 1 Comparison between critical support pressures with theoretical solutions
单位: kPa 工况 c0/kPa φ/(°) σt/kPa 破坏准则 Li等[16] Mollon等[3] 本文结果 差异/% 情况1 7 17 22.9 T-C(ξ=0) 33.92 37.44 33.84 9.6 7 17 22.9 M-C 33.61 37.07 33.53 9.5 情况2 10 25 21.4 T-C(ξ=0) 11.34 14.23 11.32 20.4 10 25 21.4 M-C 10.63 13.71 10.74 21.7 情况3 15 15 55.0 T-C(ξ=0) 14.01 18.22 13.97 23.3 15 15 55.0 M-C 10.77 15.20 10.72 29.5 表 2 开挖面临界支护力对比
Table 2 Comparison of critical support pressures of tunnel faces for σt=20 kPa
单位: kPa c0 m = 1 m = 1.2 m = 1.4 m = 1.6 T-C M-C T-C P-L T-C P-L T-C P-L 4 75.29 74.88 123.81 123.23 192.81 191.97 291.53 291.01 6 38.03 37.81 62.30 61.86 95.35 94.57 138.84 137.99 8 19.96 19.64 33.10 32.67 51.46 50.73 75.22 74.11 10 9.80 9.23 17.01 16.42 27.33 26.57 41.25 40.13 12 3.65 2.57 7.56 6.56 13.05 12.09 20.91 19.75 14 — — 1.92 0.10 4.70 3.01 8.52 7.01 16 — — — — 0.03 — 1.69 — 注:—表示不需要支护。 -
[1] PERAZZELLI P, LEONE T, ANAGNOSTOU G. Tunnel face stability under seepage flow conditions[J]. Tunnelling and Underground Space Technology, 2014, 43: 459-469. doi: 10.1016/j.tust.2014.03.001
[2] SAADA Z, MAGHOUS S, GARNIER D. Stability analysis of rock slopes subjected to seepage forces using the modified Hoek-Brown criterion[J]. International Journal of Rock Mechanics and Mining Sciences, 2012, 55: 45-54. doi: 10.1016/j.ijrmms.2012.06.010
[3] MOLLON G, DIAS D, SOUBRA A H. Rotational failure mechanisms for the face stability analysis of tunnels driven by a pressurized shield[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2011, 35(12): 1363-1388. doi: 10.1002/nag.962
[4] 吕玺琳, 王浩然, 黄茂松. 盾构隧道开挖面稳定极限理论研究[J]. 岩土工程学报, 2011, 33(1): 57-62. LÜ Xilin, WANG Haoran, HUANG Maosong. Limit theoretical study on face stability of shield tunnels[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(1): 57-62. (in Chinese)
[5] 徐前卫, 唐卓华, 朱合华, 等. 盾构隧道开挖面极限支护压力研究[J]. 岩土工程学报, 2017, 39(7): 1234-1240. doi: 10.11779/CJGE201707009 XU Qianwei, TANG Zhuohua, ZHU Hehua, et al. Limit support pressure at excavation face of shield tunnels[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(7): 1234-1240. (in Chinese) doi: 10.11779/CJGE201707009
[6] 黄阜, 潘秋景, 凌同华. 地下水渗流作用下的盾构隧道开挖面安全系数上限分析[J]. 岩土工程学报, 2017, 39(8): 1461-1469. doi: 10.11779/CJGE201708013 HUANG Fu, PAN Qiujing, LING Tonghua. Upper bound analysis of factor of safety for shield tunnel face subjected to underground water seepage[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(8): 1461-1469. (in Chinese) doi: 10.11779/CJGE201708013
[7] 黄阜, 王芬, 张芝齐, 等. 基于Hoek-Brown破坏准则的盾构隧道开挖面支护力上限研究[J]. 铁道科学与工程学报, 2018, 15(11): 2892-2900. HUANG Fu, WANG Fen, ZHANG Zhiqi, et al. Upper bound solution of supporting pressure for shield tunnel face subjected to Hoek-Brown failure criterion[J]. Journal of Railway Science and Engineering, 2018, 15(11): 2892-2900. (in Chinese)
[8] 金大龙, 袁大军, 郑浩田, 等. 高水压条件下泥水盾构开挖面稳定离心模型试验研究[J]. 岩土工程学报, 2019, 41(9): 1653-1660. doi: 10.11779/CJGE201909009 JIN Dalong, YUAN Dajun, ZHENG Haotian, et al. Centrifugal model tests on face stability of slurry shield tunnels under high water pressures[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(9): 1653-1660. (in Chinese) doi: 10.11779/CJGE201909009
[9] 徐涛, 史庆锋, 章定文, 等. 泥水盾构开挖面泥膜渗透特性与压力传递机制[J]. 岩土工程学报, 2023, 45(9): 1878-1887. doi: 10.11779/CJGE20220866 XU Tao, SHI Qingfeng, ZHANG Dingwen, et al. Permeability characteristics of filter cake and pressure transfer on face during slurry shield tunnelling[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(9): 1878-1887. (in Chinese) doi: 10.11779/CJGE20220866
[10] 陈仁朋, 尹鑫晟, 李育超, 等. 泥水盾构泥膜渗透性及其对开挖面稳定性影响[J]. 岩土工程学报, 2017, 39(11): 2102-2108. doi: 10.11779/CJGE201711018 CHEN Renpeng, YIN Xinsheng, LI Yuchao, et al. Permeability of filter cake and its influence on face stability of slurry shield-driven tunnels[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(11): 2102-2108. (in Chinese) doi: 10.11779/CJGE201711018
[11] DRUCKER D C, PRAGER W. Soil mechanics and plastic analysis or limit design[J]. Quarterly of Applied Mathematics, 1952, 10(2): 157-165. doi: 10.1090/qam/48291
[12] PAUL B. A modification of the Coulomb-Mohr theory of fracture[J]. Journal of Applied Mechanics, 1961, 28(2): 259-268. doi: 10.1115/1.3641665
[13] HOEK E, MARTIN C D. Fracture initiation and propagation in intact rock: a review[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2014, 6(4): 287-300. doi: 10.1016/j.jrmge.2014.06.001
[14] CHARLES J A, SOARES M M. The stability of slopes in soils with nonlinear failure envelopes[J]. Canadian Geotechnical Journal, 1984, 21(3): 397-406. doi: 10.1139/t84-044
[15] CHEN W F. Limit Analysis and Soil Plasticity[M]. Amsterdam: Elsevier Scientific Pub Co, 1975.
[16] LI T Z, YANG X L. Three-dimensional face stability of shallow-buried tunnels with tensile strength cut-off[J]. Computers and Geotechnics, 2019, 110: 82-93. doi: 10.1016/j.compgeo.2019.02.014
[17] MICHALOWSKI R L. Failure potential of infinite slopes in bonded soils with tensile strength cut-off[J]. Canadian Geotechnical Journal, 2018, 55(4): 477-485. doi: 10.1139/cgj-2017-0041
[18] VERMEER P A, RUSE N, MARCHER T. Tunnel heading stability in drained ground[J]. Felsbau, 2002, 20(6): 8-18.
[19] PAN Q J, DIAS D. Upper-bound analysis on the face stability of a non-circular tunnel[J]. Tunnelling and Underground Space Technology, 2017, 62: 96-102. doi: 10.1016/j.tust.2016.11.010
-
其他相关附件