Performance evaluation of shield tail brushes based on compression and grease escape tests
-
摘要: 盾尾密封失效事故时有发生,严重时易造成了重大涌水涌砂安全事故,已经引起工程界与学术界的高度关注。基于盾尾刷动态压缩试验和盾尾密封油脂逃逸试验结果及分析,结果表明:①动态压缩试验测试结果能较好反映盾尾刷弹塑性力学性能,大幅度节约试验时间;②随着盾尾间隙不断减小,盾尾刷贴合力、弹性系数不断增大,塑性变形占比不断增大、弹性变形占比不断减小。盾尾间隙为110 mm时弹性变形占比为90%,盾尾间隙压缩至50 mm时弹性变形占比降低到60%;③建立了盾尾间隙、盾尾油脂逃逸量与盾尾刷贴合力关系图,获得了盾尾密封系统的工作安全区域,可对盾尾密封状态有效性进行评估。建立的盾尾刷密封性能评价方法,为盾尾刷性能评价及检测提供了参考。Abstract: The failures of shield tail sealing have occurred from time to time, which may result in water and sand-gushing accidents in severe cases. This phenomenon has attracted great attention from the engineering and academic circles. According to the tests of dynamic compression and sealing grease escape of shield tail brushes, the results show that: (1) The dynamic compression tests can reflect the mechanical properties of the shield tail brushes, and save the test time greatly. (2) With the decrease of the shield tail gap, the adhesion force and elastic coefficient of the shield tail brushes increase continuously, the proportion of plastic deformation increases, and that of elastic deformation decreases. When the shield tail gap is 110 mm, the proportion of elastic deformation is 90%, and when it is compressed to 50 mm, the proportion of elastic deformation reduces to 60%. (3) The relationship among the gap, grease escape and adhesion force of shield tail is established, and the working safety area of the shield tail sealing system is obtained so as to evaluate the effectiveness of the status of shield tail sealing. The evaluation method for the performance of shield tail sealing is established, and can provide a reference for the performance evaluation and detection of the shield tail brushes.
-
0. 引言
潜水面上升可能诱发边坡失稳,艾伯凡矿灾[1](Aberfan Disaster)、沃楚西特水坝[2](Wachusett Dam)等就属于这一破坏模式。Brand[3]建议研究边坡的破坏机理时,应该考虑边坡内土体单元的应力应变响应。边坡内潜水面上升的应力路径可以简化为等剪应力排水路径[3-6]。针对等剪应力排水路径下土体的力学响应,研究者们进行了一系列工作[2, 6-7]。Chu等[2]通过试验证明了等剪应力排水路径下土体的非稳定态与边坡失稳存在相关性。Lashkari等[7]利用离散元数值模拟等剪应力排水路径,发现松砂和密砂的非稳定态存在区别。此外,研究者还发现各向异性对等剪应力排水路径的非稳定态存在显著影响。Fanni等[8]通过试验发现大主应力方向从垂直于沉积面旋转到45°时,非稳定态的应力比不断减小,指出在评估边坡稳定时不能忽略各向异性的影响。
一些研究者也在利用本构模型模拟该应力路径。Lashkari等[7]采用了各向异性临界状态塑性本构模型[9]模拟了不同孔隙比和剪应力值下的等剪应力排水试验,表明了模型具有模拟该应力路径的能力。但本构模型能否模拟各向异性条件下的等剪应力排水试验,目前仍然未知。
各向异性对砂土的力学响应有显著影响[10-16]。如何用本构模型描述饱和砂土的力学响应是目前的研究热点之一。为总结各向异性对于砂土力学响应影响的规律,研究者们进行了大量物理试验和离散元数值试验。物理试验方面如Yoshimine等[17]和Nakata等[18]研究了不同加载方向和中主应力系数对丰浦砂应力应变关系的影响,发现当大主应力平行于沉积面、中主应力系数为1时,土样表现出更强的剪缩特性,更易发生静态液化。离散元数值试验方面如Fu等[19]再现了直剪试验,研究各向异性砂土的峰值强度。这些试验为本构模型的建立和验证提供了基础和依据,各类模型通过模拟试验证明了自身的有效性[20-23]。如Wang等[22]通过离散元数值试验总结各向异性对塑性模量和剪胀的影响,在CycLiq本构模型[24]的基础上加入了组构各向异性,并通过离散元试验和物理试验对新模型进行了验证。
本文利用Wang等[22]提出的各向异性CycLiq本构模型,对Fanni等[8]进行的一系列定轴不排水试验与定轴等剪应力排水试验进行了模拟,验证了该模型对于各向异性砂土在不同主应力方向下、不同应力路径力学响应的模拟能力,并研究了初始组构强度对于模型模拟定轴等剪应力排水试验结果的影响。
1. 各向异性CycLiq模型
1.1 本构模型
基于排水与不排水循环扭剪试验现象研究,张建民等[25-26]将砂土的体应变εv分解为有效球应力变化引起的弹性体应变εvc和剪切引起的塑性体应变εvd,又将剪切引起的塑性体应变εvd进一步划分为可逆性剪切体应变εvd, re与不可逆性剪切体应变εv, ir。在此基础上,王刚等[27]和Wang等[24]建立了一个能够模拟饱和砂土液化前后应力应变响应的弹塑性循环本构模型,并在模型中引入了状态变量,使得模型能够统一描述砂土在不同状态下的静动力响应。为考虑各向异性,Wang等[22]基于试验和离散元数值模拟总结的规律,对模型进行了改进,引入了组构张量,在塑性模量与剪胀率的计算中引入了各向异性的影响,使得模型能够描述不同主应力方向下的响应。
CycLiq模型的应力应变关系可以用下式来描述:
˙ε=p˙r2G+(r2G+I3K)˙p+(n+D3I)⟨L⟩。 (1) 式中:p为有效球应力;r为偏应力张量s与球应力p之比;G和K为弹性剪切模量和体积模量;I为单位张量;n为塑性偏应变流动方向与加载方向;D为剪胀率;L为塑性加载指数。参数上加(·)表示增量;〈〉为Macaulay计算符,当x > 0时,〈x〉= x,当x ≤ 0时,〈x〉= 0。
弹性剪切模量和体积模量采用Richart等[28]提出的线性弹性模量公式,在此不再赘述。塑性加载指数L通过下式计算:
L=p˙r:nH。 (2) 式中:H为塑性模量。
砂土临界状态采用Been等[29]提出的状态变量Ψ =e−ec,但ec与p的关系采用了与Wang等[22]所使用的指数关系[30]不同的对数关系:
ec=Γ−λeln(p)。 (3) 这一改动是因为Fanni等[8]给出了本文模拟试样临界状态下试样的e-p曲线,并采用式(3)的对数关系进行了拟合,且e-p曲线的对数关系被广泛采用,许多模型采用了这种关系[31-32],因此对模型中的公式进行了替换。不过需要指出的是,在液化状态时,对数关系不再适用,因此砂土的动力本构模型大都采用指数关系[33-34]。本文中不涉及液化内容,因此可以选取对数关系进行模拟。
为了考虑各向异性的影响,Wang等[22]在模型中引入了组构张量Fn,并考虑了组构张量对于塑性模量与剪胀率的影响。组构张量Fn的演化采用了如下公式:
dFn=⟨L⟩c[n−(1+D)Fn]。 (4) 式中:c为模量参数,代表组构张量演化的速率。通过总结离散元数值试验的规律,Wang等[22]发现组构张量的模||Fn||与组构张量和加载方向的联合不变量An =Fn: n之差||Fn||-An可以作为一个良好的指标,以反映塑性模量和剪胀率随加载角度的变化。以水平沉积的试样在三轴压缩中的加载角度为0°,随着加载角度的增加,指标||Fn||-An不断增大,同时塑性模量不断减小,剪胀率从不断增加。因此,将||Fn||-An引入到了塑性模量与剪胀率的计算公式中:
C=exp[−Δ1(‖Fn‖−An)1/122], (5) H=23hg(θ)Gexp(−npΨ)⋅(Mexp(−npΨ)Mmˉρρ−1C), (6) Dir=dir+Δ2(‖Fn‖−An)dirDir, iso。 (7) 式中:h,Δ1,Δ2,M,dir,np为模型参数;g(θ)为形状函数,这里采用的是Zhang等[35]建议的函数;Mm为历史最大应力比;Dir, iso为各向同性条件下的不可逆剪胀率,详细公式可见Wang等[24]。
1.2 参数标定
Fanni等[8]试验中采用的砂为一种细硅砂,颗粒相对质量密度为2.64,平均粒径D50为0.21 mm,均匀系数Cu为2.3,颗粒呈圆形至亚圆形。因缺乏相关试验,大部分参数参照了Wang等[22]标定的丰浦砂参数,仅对部分可确定的临界状态参数(np,Γ,λe)进行了重新标定。标定结果如表 1所示。
其中,Γ和λe由Fanni等[8]根据试验拟合得到,在更换了Wang等[22]模型中ec与p的关系式后可直接使用。图 1展示了np的标定方法,以Fanni等[8]进行的一组定轴等剪应力排水试验(中主应力系数b = 0.2,加载角度α = 30°)结果参照进行标定。图中点代表试验所得试样到达非稳定态时的有效球应力与孔隙比,虚线和实线分别为np = 6和12时的模拟结果。可以看到当其他参数不变时,np值越大,模拟得到的非稳定态时有效球应力随孔隙比变化越快,曲线更加陡峭。这是由于等剪应力排水路径中的非稳定态代表试样处于峰值应力比状态,而峰值应力比可由Mexp(-npΨ)决定。
除了表 1中列出的模型参数外,初始组构强度,即初始组构张量的模||Fnin||,也是一个重要的参数。它代表了试样在试验开始时的各向异性程度,||Fnin||越大,试样的各向异性越强。不同的制样方式制得的试样往往有着不同的各向异性程度[36],因此需要取不同的初始组构强度。本次模拟中涉及两种不同的制样方式,湿装夯实法(moist tamped, MT)和干砂雨法(dry pluviated, DP),其中干砂雨法得到的试样孔隙比要小于湿装夯实法,结合Yang等[37]、叶斌等[38]通过物理试验与离散元模拟均发现干砂雨法的各向异性强于湿装夯实法,可以推测干砂雨法在制样中,由于受到重力影响,土颗粒的排列具有定向性,而湿装夯实法在制样中,土颗粒更倾向于随机排列,导致其各向异性相对较弱。最终,湿装夯实法制得的试样初始组构强度||Fnin||取为0.5,干砂雨法制得的试样初始组构强度||Fnin||取为0.7。
2. 不同主应力方向下模型的响应
2.1 定轴不排水剪切试验
Fanni等[8]进行了一系列空心圆柱定轴不排水剪切试验。试验中,初始有效球应力均为171 kPa,大小主应力之差均为64 kPa,中主应力系数b分别为0,0.2,0.5,大主应力与沉积方向的夹角α分别为0°,22.5°,30°和45°,两种制样方法,即湿装夯实法和干砂雨法,各5组试验,如图 2,3所示。其中,湿装夯实法所得试样孔隙比在0.761~0.785,干砂雨法所得试样孔隙比在0.674~0.697。采用表 1中的模型参数对试验进行模拟,模拟结果如图 2,3所示。图 2(a),(b)中横纵坐标为Fanni等[8]给出的八面体剪应变和八面体偏应力,图 3同理。
表 1 各向异性CycLiq模型参数Table 1. Model parameters of anisotropic CycLiq modelGo κ h dre, 1 dre, 2 dir α γd, r np nd M λe Γ Δ1 Δ2 c 100 0.008 4.0 0.7 30 0.3 30 0.05 12 2 1.2 0.008 0.761 1.35 0.8 5 图 2展示了湿装夯实法试样定轴不排水试验结果与模拟。湿装夯实法得到的试样均为松砂,各试验和模拟均出现静态液化。总体上来看,随着b和α的增加,试验与模拟中试样的峰值偏应力均不断减小。从图 2(a),(c)中可以看出,试验中b = 0.2的3组曲线极为接近,甚至α = 30°的偏应力要高于α = 22.5°,这主要是试验中不同试样的孔隙比相差较大造成的,在模拟中该现象并不明显。此外,b = 0.5的模拟结果在加载初期剪应力出现了明显的减小,这主要是由于本构模型中组构对塑性模量的影响描述不够准确,在大孔隙比下放大了组构对塑性模量的影响。
图 3展示了干砂雨法试样定轴不排水试验结果与模拟。相比于湿装夯实法,干砂雨法得到的试样孔隙比更小,可认为密砂。各组模拟的曲线与试验结果均有较好的吻合。对于b = 0.2的3组试验和模拟,α = 22.5°时均表现出明显的剪胀,而在α = 45°时则表现出强烈的剪缩。对于b = 0,α = 0°和b = 0.2,α = 22.5°两组试验,试验和模拟的曲线在图 3(a),(b)均出现了交叉,这是由于两者的孔隙比相差较大,在孔隙比、中主应力系数、主应力方向三者的共同作用下出现了图中的效果。
2.2 定轴等剪应力排水试验
在定轴不排水剪切试验的基础上,Fanni等[8]进行了一系列定轴等剪应力排水试验,以探究各向异性对等剪应力排水应力路径的影响。试验中采用了湿装夯实法(MT)和干砂雨法(DP)两种制样方法,初始有效球应力均为171 kPa,大小主应力之差均为64 kPa,中主应力系数b分别为0,0.2,主应力与沉积方向的夹角α分别为0°,22.5°和45°。每个夹角均用不同的孔隙比进行试验,对于湿装夯实法,b = 0,α = 0°组采用了6个不同的孔隙比,b = 0.2,α = 22.5°组采用了4个,b = 0.2,α = 45°组采用了10个;对于干砂雨法, b = 0.2,α = 22.5°组1个,b = 0.2,α = 45°组2个。采用表 1中的模型参数进行模拟,模拟中保持偏应力不变,当有效球应力无法再减小时,即试样达到峰值应力比时,认为此时试样处于非稳定态,模拟结束。
图 4给出了湿装夯实法试样在不同b和α下的定轴等剪应力排水试验结果与模拟。图中圆代表初始状态,三角形代表试样达到非稳定态,即处于峰值应力比状态。当相同b和α时,非稳定态的有效球应力p随孔隙比e的增加而增加,3组中模拟与试验均吻合较好。当孔隙比e不变、α增加时,模拟与试验得到的有效球应力均增加。
图 5给出了干砂雨法试样在不同b和α下的定轴等剪应力排水试验结果与模拟。b = 0.2,α = 45°组试验与模拟结果相近,规律也与湿装夯实法中一致;模拟结果在定量上与试验结果存在差距,推测这是由于模型参数或者初始组构强度不准确导致的。当α增加时,模拟与试验得到的有效球应力均增加,这与湿装夯实法中的规律是一致的。
图 6给出了试验和模拟中非稳定态的结果与趋势线。可以看到试验得到的非稳定态在e-p空间内具有良好的线性关系,对模拟得到的非稳定态进行线性回归,得到R2均在0.99以上,表明模型对等剪应力比试验具有良好的模拟能力。对比湿装夯实法的3组试验,可以发现模拟得到的3条趋势线之间的间距要略小于试验所得,这主要受初始组构强度的影响,具体情况在下文中阐释。
为了说明初始组构强度对于定轴等剪应力排水条件下土体响应的影响,在数值模拟中对不同初始组构强度条件进行了分析。图 7中给出了不同初始组构强度对非稳定态模拟结果的影响。从图中可以发现初始组构越强,不同主应力方向下的非稳定态差异越大,即峰值应力比差异越大。式(5),(6)中可以解释这一现象。试样达到峰值应力比时,有塑性模量H = 0。当初始组构强度增加,即||Fn||增加时,指标||Fn||-An = ||Fn||(1 – nF: n)(其中nF指组构张量的方向)的值域[0, 2||Fn||]随之扩大,式(5)所得C的值域也相应扩大。对于定轴试验,初始组构强度增加,导致了不同的b和α计算得到C值之间的差异增加,因H = 0,进而使得式(6)中历史最大峰值应力比Mm(即非稳定态时的应力比)的差异增加,在剪应力q不变的情况下,非稳定态时的球应力p之间的差异也就增加了。这些反映在图 7中即为趋势线之间的间距增加。
3. 结论
本文利用各向异性CycLiq本构模型,模拟了Fanni等[8]进行的定轴不排水试验与定轴等剪应力排水试验,对比试验结果与模拟结果,可以得到以下3点结论。
(1)对于定轴不排水剪切试验,该模型与试验均表现为随着中主应力系数和加载角度的增加,土体剪缩增强,更易发生静态液化。在定量上,模型模拟的应力路径和应力应变曲线与试验也有较好的吻合。
(2)对于定轴等剪应力排水试验,模型与试验均表现为非稳定态时的平均有效应力随孔隙比和加载角度的增加而增加,且相同加载角度的非稳定态在e-p空间内呈线性分布。总体来说,CycLiq模型对各向异性下砂土等剪应力排水响应的模拟效果较好,具有运用到实际工程中的潜力。
(3)不同初始组构强度对于等剪应力排水试验的结果有显著的影响。初始组构强度越高,不同主应力方向作用下的非稳定态时平均有效应力相差越大。
-
表 1 盾尾刷力学性能指标
Table 1 Mechanical parameters of shield tail brush
盾尾间隙/mm 盾尾刷贴合力/N 压紧板残余高度/mm 压紧板残余角度/(°) 弹性系数/(N·mm-1) 50 1362.2 135 35.6 16 75 997.2 159 44.5 11.9 90 392 174 48 4.7 110 294 183 51.4 4.02 -
[1] 何 川, 封 坤, 方 勇. 盾构法修建地铁隧道的技术现状与展望[J]. 西南交通大学学报, 2015, 50(1): 97-109. doi: 10.3969/j.issn.0258-2724.2015.01.015 HE Chuan, FENG Kun, FANG Yong. Review and Prospects on Constructing Technologies of Metro Tunnels Using Shield Tunnelling Method [J]. Journal of Southwest Jiaotong University, 2015, 50(1): 97-109. (in Chinese) doi: 10.3969/j.issn.0258-2724.2015.01.015
[2] 朱 伟, 陈仁俊. 盾构隧道施工技术现状及展望(第3讲)——盾构隧道应用前景及发展方向[J]. 岩土工程界, 2002, 5(1): 18-20, 52. doi: 10.3969/j.issn.1674-7801.2002.01.040 ZHU Wei, CHEN Ren-jun. Current status and prospects of shield tunnel construction technology (lecture3)-application prospects and development directions of shield tunnels[J]. Geotechnical Engineering, 2002, 5(1): 18-20, 52. (in Chinese) doi: 10.3969/j.issn.1674-7801.2002.01.040
[3] 李艳春. 盾构法隧道施工中盾尾刷的优化与保护[J]. 湖南城市学院学报(自然科学版), 2015, 24(3): 33-34. doi: 10.3969/j.issn.1672-7304.2015.03.012 LI Yanchun. Shield tunneling optimization and protection shield tail brush[J]. Journal of Hunan City University (Natural Science), 2015, 24(3): 33-34. (in Chinese) doi: 10.3969/j.issn.1672-7304.2015.03.012
[4] 任广艳. 透水涌砂冒险堵漏隧道坍塌撤人不及: 广东省佛山市轨道交通2号线"2.7"透水坍塌重大事故分析[J]. 吉林劳动保护, 2019(8): 40-42. REN Guangyan. Seepage sand gushing, risk plugging, tunnel collapse, evacuation failure—analysis of the "2.7" serious accident of seepage collapse in Foshan Rail Transit Line 2, Guangdong Province[J]. Jilin Labour Protection, 2019(8): 40-42. (in Chinese)
[5] YU C, ZHOU A N, CHEN J, et al. Analysis of a tunnel failure caused by leakage of the shield tail seal system[J]. Underground Space, 2020, 5(2): 105-114. doi: 10.1016/j.undsp.2018.11.003
[6] 李奕, 钟志全. 一种新型盾尾刷的设计与应用[J]. 建筑机械化, 2011, 32(1): 82-84. doi: 10.3969/j.issn.1001-1366.2011.01.020 LI Yi, ZHONG Zhiquan. Design and application of a new-type shield tail brush[J]. Construction Mechanization, 2011, 32(1): 82-84. (in Chinese) doi: 10.3969/j.issn.1001-1366.2011.01.020
[7] 霍志光. 新型盾尾刷设计与应用分析[J]. 价值工程, 2011, 30(8): 87. doi: 10.3969/j.issn.1006-4311.2011.08.059 HUO Zhiguang. Design and application analysis of new shield tail brush[J]. Value Engineering, 2011, 30(8): 87. (in Chinese) doi: 10.3969/j.issn.1006-4311.2011.08.059
[8] 沈桂丽, 刘金祥, 李剑雄, 等. 盾尾密封刷实验平台的设计与研究[J]. 隧道建设, 2015, 35(8): 851-854. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201508023.htm SHEN Guili, LIU Jinxiang, LI Jianxiong, et al. Design of and study on tail brush experimental platform[J]. Tunnel Construction, 2015, 35(8): 851-854. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201508023.htm
[9] WEI Linchun. Experimental study on mechanical bahavior of wire brushes on shield tail[J]. Tunnel Construction, 2021, 41(2): 206-211.
[10] YE Guan-lin, HAN Lei, SANTOSH Kumar Yadav, et al. Investigation on the tail brush induced loads upon segmental lining of a shield tunnel with small overburden[J]. Tunnelling and Underground Space Technology, 2020, 97: 103283. doi: 10.1016/j.tust.2020.103283
[11] 钟小春, 莫暖娇, 余明学, 等. 盾尾刷环形密封系统单元试验及水密性机制研究[J]. 岩土工程学报, 2023, 45(2): 354-361. doi: 10.11779/CJGE20211464 ZHONG Xiaochun, MO Nuanjiao, YU Mingxue, et al. Study on the unit test of shield tail brush annular sealing system and its watertightness mechanism[J]. Journal of Geotechnical Engineering, 2023, 45(2): 354-361. (in Chinese) doi: 10.11779/CJGE20211464
[12] 黄旭民, 黄林冲, 梁禹. 施工期同步注浆影响下盾构隧道管片纵向上浮特征分析与应用[J]. 岩土工程学报, 2021, 43(9): 1700-1707. doi: 10.11779/CJGE202109015 HUANG Xumin, HUANG Linchong, LIANG Yu. Analysis and application of longitudinal uplift characteristics of segments of shield tunnels affected by synchronous grouting during construction period[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(9): 1700-1707. (in Chinese) doi: 10.11779/CJGE202109015
[13] 王德乾, 郑筱彦, 斯芳芳, 等. 盾尾密封油脂与盾构机及施工地层的适应性研究[J]. 铁道建筑技术, 2021(7): 8-12, 69. https://www.cnki.com.cn/Article/CJFDTOTAL-TDJS202107002.htm WANG Deqian, ZHENG Xiaoyan, SI Fangfang, et al. Study on adaptability of shield tail sealing paste, shield machine and shield construction strata[J]. Railway Construction Technology, 2021(7): 8-12, 69. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TDJS202107002.htm