射流联合盾构切削钢筋混凝土仿真与试验

庄欠伟; 袁一翔; 徐天明; 张弛

doi:10.11779/CJGE202010006

射流联合盾构切削钢筋混凝土仿真与试验 English Version

庄欠伟^{1, 2,},
袁一翔^1, ,,
徐天明²,
张弛²

1.
上海大学土木工程系,上海　200444
2.
上海隧道工程有限公司,上海200238

基金项目:

上海市科委项目 18DZ2283200

上海市国资委科研项目 2017014

详细信息

作者简介:
庄欠伟（1978— ）,男,教授级高级工程师,硕士生导师,主要从事地下空间与工程方面教学和研究工作。E-mail：qianweizhuang@sohu.com

通讯作者:
袁一翔， E-mail：yyxyzdd@163.com

中图分类号: U455.43
计量
- 文章访问数: 312
- HTML全文浏览量: 17
- PDF下载量: 193
出版历程
- 收稿日期: 2020-02-03
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-09-30

Simulation and experiment on cutting reinforced concrete with jet combined shield method

1.
Department of Civil Engineering, Shanghai University, Shanghai 200444, China
2.
Shanghai Tunnel Engineering Co., Ltd., Shanghai 200238, China

摘要

摘要: 针对传统盾构刀具直接切削钢筋混凝土所产生的弊端,进行磨料水射流联合盾构刀具切削钢筋混凝土研究。通过仿真和试验探究了磨料水射流在不同横移速度下对钢筋和素混凝土的切削效果和切割机理并确定了切割钢筋混凝土的最佳横移速度。研究结果表明：磨料水射流对钢筋和混凝土的切割机理存在差异性,水射粒子对钢筋的切割贡献十分有限；混凝土的切口深度和宽度大于钢筋,对横移速度变化的敏感性小于钢筋；在高横移速度下混凝土切口损伤宽度得到扩大但并不显著、钢筋出现切割残余；切割深度随着横移速度的增加而减小,下降速度由陡变缓；10 mm/min的横移速度较好满足切割C55钢筋混凝土的要求。基于前期研究成果开展模拟盾构切削试验,分析了磨料水射流联合盾构刀具破除C55钢筋混凝土的效果、刀具损伤和参数特征。试验结果表明：磨料水射流联合盾构刀具切削能有效解决盾构刀具直接切削所产生的钢筋缠绕问题,减少了刀具损伤,控制了扭矩峰值。
- 隧道开挖 /
- 磨料水射流 /
- 盾构刀具 /
- 钢筋混凝土 /
- 仿真
Abstract: The research of abrasive water jet combined with shield cutting tools in the cutting of reinforced concrete is carried out because of the disadvantages of the traditional shield cutting tools in the direct cutting of reinforced concrete. The cutting effect and mechanism of abrasive water jet on reinforced concrete and plain concrete under different moving speed are explored, and the best moving speed of cutting reinforced concrete is determined through simulation and experiment. The results show that the cutting mechanism of abrasive water jet is different from that of concrete. The contribution of water jet particles to steel cutting is limited. The depth and width of concrete incision are greater than those of steel. The sensitivity of concrete to the change of transverse velocity is less than that of steel bar. The damage width of concrete notch is enlarged but not obvious, and the cutting residues occur in the reinforcement under high transverse velocity. The cutting depth decreases with the increase of the traverse speed. The descending speed slows down. The transverse velocity of 10 mm/min can meet the requirements of cutting C55 reinforced concrete. The cutting experiment of the simulation shield machine is carried out based on the previous research results. The effect, tool damage and parameter characteristics of the abrasive water jet combined shield machine tool to break C55 reinforced concrete are analyzed, and the results are compared with those of the shell knife direct cutting. The results show that the abrasive water jet combined with shield cutter cutting can effectively solve the problem of steel winding caused by the direct cutting of shield cutter, reduce the cutter damage and control the peak torque.
- tunnel excavation /
- abrasive water jet /
- shield cutter /
- reinforced concrete /
- simulation

HTML全文

0. 引言

隧洞涌水问题一直是地下工程建设中突出挑战之一,大量的涌水会造成隧洞失稳、人员伤亡、地面沉降及生态环境恶化等一系列恶劣后果^[1-2]。因此,隧洞涌水量的预测计算对防排水设计、衬砌措施选择以及环境影响评价等方面具有十分重要的意义^[3-4]。

目前,隧洞涌水量计算方法主要分为经验公式法、解析公式法及数值计算法3类,其中解析公式法由于其计算快速简洁、论推导严谨的优点而得到广泛地应用。解析公式法可分为镜像法、竖井法、保角变换法及其他方法^[5-8],但这些方法都基于不考虑地下水位的变化,隧洞处于地下水无限补给这一初始假定条件推导而出的。对于大量的实际工程而言,隧洞周围地下水往往不存在无限补给源,隧洞的开挖将疏干其周围地下水,进一步导致地下水位的降落。因此,采用传统解析公式计算疏干隧洞涌水量会造成计算结果明显偏大。

针对疏干条件下隧洞涌水量变化与地下水位降落等问题,已有不少学者开展了研究。Moon等^[9]利用离散单元法研究出疏干条件下裂隙岩体中隧洞地下水位的变化规律,提出将降低地下水位 $\bar{h}$ 代替初始地下水位 $h$ ,在Goodman公式基础上计算隧洞涌水量；苏凯等^[10]提出了利用渗透水头替代初始地下水位,并给出渗透水头的经验公式,简化了涌水量公式的使用。EI Tani等^[11]基于数值积分表达式与格林函数法,推导出单侧补给条件下隧洞的地下水位降深与涌水量解析公式。上述学者虽然采用多种方法推导出了疏干条件下隧洞涌水量解析公式,但是这些解析公式的推导是在假设稳定地下水水位位于隧洞顶部上方的情况下进行的,主要适用于深埋隧洞的地下水疏干问题。对于浅埋隧洞而言,其顶部上方地下水被完全疏干,上述解析公式无法真实地刻画浅埋隧洞的渗流特征,而且也无法计算出隧洞孔隙水压力的分布规律；此外,当前疏干条件下隧洞涌水量的研究仍着力在未衬砌隧洞上,对于衬砌隧洞,涌水量与衬砌水压力的分布规律仍未得到很好的解决。

本文针对浅埋、深埋及衬砌隧洞进行解析分析,然后利用数值模拟方法计算出不同隧洞半径、初始地下水位高度及围岩与衬砌相对渗透性下隧洞涌水量及水压力分布规律,将数值计算结果与解析分析相结合,提出了计算隧洞涌水量和孔隙水压力的解析公式。

1. 解析分析

1.1 浅埋隧洞

疏干条件下浅埋隧洞地下水位线快速降落至隧洞中心线两侧,形成稳定的降落漏斗,隧洞顶部地下水被完全疏干,见图1。将其周围地下水运动视为满足水平一维流,其数学模型为

$\frac{d}{d x} (u \frac{d u}{d x}) = 0$ ,

(1)

其边界条件为

$\begin{array}{l} u |_{x = r} = u_{0}, \\ u |_{x = L} = u_{L}, \end{array}}$

(2)

图 1 浅埋隧洞渗流场示意图

Figure 1. Schematic representation of shallow tunnel

下载: 全尺寸图片幻灯片

式中, $u$ 为隧洞右侧 $x$ 处的含水层厚度, $u_{0}$ , $u_{L}$ 分别为隧洞右侧边界与距中心L处的含水层厚度；自变量x取值范围为 $r \leq x \leq L$ 。

方程（1）对应边界条件式（2）的解为

$u^{2} = \frac{u_{L}^{2} - u_{0}^{2}}{L - r} x + (u_{0}^{2} - \frac{u_{L}^{2} - u_{0}^{2}}{L - r} r)$ 。

(3)

根据Darcy定律可得任意断面潜水流的单宽流量为

$q_{x} = k_{m} u \frac{d u}{d x}$ ,

(4)

式中, $q_{x}$ 为距隧洞右侧 $x$ 处任意断面上的单宽流量, $k_{m}$ 为围岩渗透系数。因此,隧洞中心处单宽流量为

$q_{x = 0} = k_{m} u \frac{d u}{d x} = k_{m} \frac{u_{L}^{2} - u_{0}^{2}}{2 (L - r)}$ 。

(5)

考虑到隧洞渗流场的对称性,可得浅埋隧洞涌水量公式：

$Q = \frac{k_{m}}{L - r} (u_{L}^{2} - u_{0}^{2}) = \frac{k_{m}}{L - r} (2 h D + h^{2})$ 。

(6)

实际上,浅埋隧洞底部地下水为径向流运动,这造成了地下水水头的额外损失,并不符合水平流的假定,导致了计算结果明显偏大。Hooghoudt^[12]于1940年在解决降水条件下暗管排水量的计算问题时,引入影响距离d来消除径向流所带来的的误差。

因此,本文借鉴相关思路,同样利用影响距离d来替代隧洞至底板距离“D”,进一步推导出浅埋隧洞涌水量计算公式,可表达为

$Q = \frac{k_{m}}{L - r} (2 h d + h^{2})$ 。

(7)

1.2 深埋隧洞

对于深埋隧洞而言,其周围地下水为径向流运动,满足拉普拉斯方程。Harr^[3]于1962年利用镜像法推导出无限补给条件下隧洞渗流场内任一点涌水量、水头与孔隙水压力分布的解析公式：

$Q = 2 π k_{m} \frac{h}{\ln {[1 + {(2 h / r)}^{2}]}^{0.5}} \approx \frac{2 π k_{m} h}{\ln (2 h / r)}$ ,

(8)

$H = [1 - \frac{\ln (\cos^{2} θ + {(2 h / x - \sin θ)}^{2})}{\ln (1 + {(2 h / r)}^{2})}] \times h$ ,

(9)

$p = [1 - \frac{\ln (\cos^{2} θ + {(2 h / x - \sin θ)}^{2})}{\ln (1 + {(2 h / r)}^{2})} - \frac{x \sin θ}{h}] \times h \times γ_{w} 。$

(10)

式中 $θ$ 为流场内任一点与隧洞中心线之间的逆时针夹角；x为流场内任一点至隧洞中心的距离； $γ_{w}$ 为水的重度；其他参数见上式。

疏干条件下深埋隧洞地下水位出现一定的降落,地下水未被完全疏干,见图2,其中 $h$ 为初始地下水位, $\bar{h}$ 为降落地下水位,s为地下水位降深。相对比无限补给状况,疏干条件下深埋隧洞的地下水水头与涌水量发生显著变化,简单地采用渗透水头替代初始地下水位是不合适的,也不符合镜像法的方法原理,同时会造成隧洞孔隙水压力计算结果的偏小。

图 2 深埋隧洞渗流场示意图

Figure 2. Schematic representation of deep tunnel

下载: 全尺寸图片幻灯片

本文将降落地下水位等效转化为未降落的等效地下水位 $h^{'}$ ,然后利用镜像法进而推导出疏干条件下深埋隧洞涌水量解析公式：

$Q^{'} = 2 π k_{m} \frac{h^{'}}{ln(2 h^{'} / r)}$ ,

(11)

式中, $h^{'}$ 为深埋疏干隧洞等效地下水位高度。

易知,深埋隧洞中心水平向与垂向上孔隙水压力表达式分别为

$p_{水平} = [1 - \frac{\ln (1 + {(2 h^{'} / x)}^{2})}{\ln (1 + {(2 h^{'} / r)}^{2})}] \times h^{'} \times γ_{w}$ ,

(12)

$p_{垂直} = [1 - \frac{\ln ({(2 h^{'} / x - 1)}^{2})}{\ln (1 + {(2 h^{'} / r)}^{2})} - \frac{x}{h^{'}}] \times h^{'} \times γ_{w}$ 。

(13)

1.3 衬砌隧洞

对于衬砌隧洞而言,衬砌圈的存在改变了隧洞周围渗流场,造成孔隙水压力与地下水位的重新分布,地下水水头的损失造成经衬砌圈进入隧洞的涌水量显著减少,见图3。Fernandez等^[13-14]给出了无限补给条件下衬砌隧洞涌水量的解析公式：

$Q = 2 π k_{m} \frac{h - h_{l}}{\ln (2 h / r_{l})}$ ,

(14)

$Q_{l} = 2 π k_{l} \frac{h_{l}}{\ln (r_{l} / r)}$ 。

(15)

图 3 衬砌隧洞渗流场示意图

Figure 3. Schematic representation of lined tunnel

下载: 全尺寸图片幻灯片

式中 $Q_{l}$ 为衬砌隧洞涌水量； $Q$ 为衬砌圈外侧涌水量； $k_{l}$ 为衬砌圈渗透系数； $r$ , $r_{l}$ 分别为隧洞半径与衬砌圈半径； $h_{l}$ 为衬砌圈水头。

疏干条件下衬砌隧洞的稳定地下水位会小于其初始地下水位,地下水降深取决于围岩与衬砌相对渗透性^[15]。同理,运用等效转化思想,将衬砌隧洞的降落地下水位转化为衬砌等效地下水位 ${h^{'}}_{l}$ ,进一步地,疏干条件下衬砌隧洞涌水量计算公式可表达为

$Q = 2 π k_{m} \frac{{h^{'}}_{l} - h_{l}}{\ln (2 {h^{'}}_{l} / r_{l})}$ ,

(16)

式中, ${h^{'}}_{l}$ 为衬砌隧洞等效地下水位,其他参数见上述公式。

为简化运算,令

$A = \frac{\ln (2 {h^{'}}_{l} / r_{l})}{\ln (r_{l} / r)}$ ,

根据流量连续性原理,令 $Q = Q_{l}$ ,可得

$h_{l} = \frac{{h^{'}}_{l}}{1 + A \times k_{l} / k_{m}}$ 。

(17)

而衬砌状态下隧洞中心水平向与垂向上衬砌水压力表达式为

$p_{水平} = [1 - \frac{\ln (1 + {(2 {h^{'}}_{l} / x)}^{2})}{\ln (1 + {(2 {h^{'}}_{l} / r_{l})}^{2})} \times \frac{A (k_{l} / k_{m})}{1 + A (k_{l} / k_{m})}] \times {h^{'}}_{l} \times γ_{w}$ ,

(18)

$p_{垂直} = [1 - \frac{\ln ({(2 {h^{'}}_{l} / x - 1)}^{2})}{\ln (1 + {(2 {h^{'}}_{l} / r_{l})}^{2})} \times \frac{A (k_{l} / k_{m})}{1 + A (k_{l} / k_{m})} - \frac{x}{{h^{'}}_{l}}] \times {h^{'}}_{l} \times γ_{w} 。$

(19)

对于浅埋、深埋及衬砌隧洞涌水量和衬砌水压力的解析公式中的关键参数（影响距离d,等效地下水位 $h^{'}$ ,衬砌等效地下水位 ${h^{'}}_{l}$ ）,可以基于数值模拟方法确定。

2. 数值模拟研究

2.1 计算方法与数值模型

疏干条件下隧洞周围渗流场处于不断变化之中,渗流场求解的难点在于隧洞边界条件的表达与溢出点奇异性的确定^[16]。郑宏等^[17]提出了一个新的Signorini型变分不等式提法,将潜在溢出面边界提为Signorini型互补边界条件,从理论上解决了溢出点的奇异性,已广泛地应用到水利水电工程中。本文基于Signorini型变分不等式法开发出有限元程序,计算疏干条件下隧洞渗流场。

假定取不同的隧洞半径r（2.5,5.0,7.5,10.0 m）、初始地下水位高度h（25～300 m）与围岩与衬砌相对渗透性（ $k_{m} / k_{l}$ =1,100,1000,10000, $k_{m}$ =1 m/d）,建立有限元模型,其中模型范围为两侧外部边界及底部取100倍隧洞半径^{[10, 15]},初始地下水位埋深取50 m,衬砌厚度取1.5 m,见图4。根据Signorini型变分不等式原理,模型边界条件设定：①左右两侧边界 $EB \cup FD$ 取定水头边界；②底边界BD取隔水边界；③隧洞边界与边界取Signorini型边界。

图 4 有限元模型网格剖分示意图

Figure 4. Finite element mesh of models

下载: 全尺寸图片幻灯片

2.2 隧洞涌水量

图5为各种数值案例下隧洞涌水量的变化曲线。可以看出：当隧洞半径与初始地下水位高度一定时,以隧洞h=100 m为例,随着围岩与衬砌相对渗透性 $k_{m} / k_{l}$ 的增大,隧洞涌水量由98.44 m³/d（ $k_{m} / k_{l}$ =1）显著减少至0.24 m³/d（ $k_{m} / k_{l}$ =10000）；当 $k_{m} / k_{l}$ 一定时,随着隧洞初始地下水位高度逐渐增大,疏干隧洞涌水量由17.11 m³/d（h=25 m）逐渐增大至329.53 m³/d（h=300 m）,呈指数形式增长,而衬砌隧洞涌水量增幅较小,4.66 m³/d增至60.59 m³/d（ $k_{m} / k_{l}$ =100）,说明衬砌圈具有良好的阻水作用。

图 5 隧洞涌水量数值计算结果

Figure 5. Water inflows of tunnels in different cases

下载: 全尺寸图片幻灯片

2.3 孔隙水压力分布

为了探讨初始地下水位高度对隧洞孔隙水压力分布的影响,以隧洞r=5.0 m为例,绘制出不同初始地下水位高度下疏干隧洞水平向和垂向的孔隙水压力分布曲线,见图6。

图 6 不同初始地下水位高度下疏干隧洞水平向与垂向的孔隙水压力分布

Figure 6. Distribution of pore water pressure along horizontal and vertical directions of drained tunnel under various initial groundwater levels

下载: 全尺寸图片幻灯片

由图6可知：在疏干隧洞水平向上,随着隧洞边界右侧距离的增加,孔隙水压力值逐渐增大,且在距离隧洞100 m范围内,孔压变化幅度最为剧烈,之后增幅变缓；在垂向上,隧洞边界处孔隙水压力为0 MPa,随着距隧洞顶拱的距离逐渐增大,孔隙水压力值先迅速增大到最大值,然后缓慢减小,在渗流自由面处重新降低为0 MPa。

以隧洞r=5.0 m,h=200 m为例,绘制出4种围岩与衬砌相对渗透性下隧洞水平向和垂向的孔隙水压力分布,如图7所示。可以看出：不同衬砌方案下隧洞孔隙水压力分布基本相似,在水平向上,在隧洞边界至衬砌圈范围内,孔隙水压力值迅速增加,且随着围岩与衬砌相对渗透性的增大,孔压由1.623 MPa ( $k_{m} / k_{l}$ =100)增至1.984 MPa( $k_{m} / k_{l}$ =10000)；在垂直向上,孔隙水压力值在衬砌圈表面快速增加至最大值,之后呈线性递减至0 MPa,这说明了衬砌圈的存在显著改变了隧洞孔隙水压力分布规律。

图 7 不同围岩与衬砌相对渗透性下衬砌隧洞水平向与垂向的.孔隙水压力分布

Figure 7. Distribution of pore water pressure along horizontal and vertical directions of lined tunnel with different k_m/k_l

下载: 全尺寸图片幻灯片

3. 讨论与验证

3.1 涌水量解析公式

（1）影响距离

将浅埋隧洞（ $r / h \geq 0.1$ ）计算工况结果汇总,绘制出浅埋隧洞d/D关于r/h之间的散点图,见图8。可以看出：浅埋隧洞d/D值分布在0.1～0.4内,关于呈线性分布趋势；在隧洞半径r一定时,影响距离随初始地下水位高度的增加而逐渐增大,以r=5.0 m为例,d/D由初始地下水位高度25 m（r/h=0.2）时0.314增加至初始地下水位高度50 m（r/h=0.1）时0.382；在初始地下水位高度一定时,影响距离随隧洞半径的增大而增大。

图 8 d/D关于r/h的散点分布及拟合曲线

Figure 8. Scatter distribution and fitting curves of d/D with respect to r/h

下载: 全尺寸图片幻灯片

针对浅埋隧洞d/D与r/h之间的分布规律,将定义为关于d/D的线性函数,表达式如下：

$d = (a r / h + b) \cdot D$ ,

(20)

式中,a与b分别为待定系数。

进一步地,根据图8中散点数据进行曲线拟合,可得到系数：

$\begin{array}{l} a = 0.0193 r^{2} - 0.3223 r + 0.3654, \\ b = - 0.0017 r^{2} + 0.0311 r + 0.3538, \end{array}}$

(21)

可知,浅埋隧洞涌水量半解析公式为

$\begin{array}{l} Q = \frac{k_{m}}{L - r} (2 h d + h^{2}), \\ d = (a r / h + b) \cdot D, \end{array}}$

(22)

式中,a与b见式（21）。

图8中相关性系数R²均大于0.95,说明半解析解与数值解具有良好的拟合度,验证了表达式具有可行性。

（2）等效地下水位

将深埋隧洞（ $r / h < 0.1$ ）计算工况结果汇总,绘制出深埋隧洞等效地下水位与初始地下水位之比 $h^{'} / h$ 关于 $r / h$ 之间的散点图,见图9。可以看出：深埋隧洞 $h^{'} / h$ 值处于0.3～1.0之间,关于 $r / h$ 呈线性分布；在隧洞半径r一定时,随着隧洞初始地下水位高度逐渐增大,等效地下水位与初始地下水位比值 $h^{'} / h$ 逐渐增大；在隧洞初始地下水位h一致时,洞径越大, $h^{'} / h$ 逐渐减小,这说明地下水位降深随洞径的增大而逐渐增大,从而导致隧洞上方的等效地下水位逐渐减少。

图 9

$h^{'} / h$ 关于

$r / h$ 的散点分布及拟合曲线

Figure 9. Scatter distribution and fitting curves of

$h^{'} / h$ with respect to

$r / h$

下载: 全尺寸图片幻灯片

为了便于计算,同理,将 $h^{'} / h$ 定义为关于 $r / h$ 的线性函数,表达式如下：

$h^{'} = (c r / h + d) \cdot h$ ,

(23)

式中,a,b分别为待定系数。

根据图9中数据进行曲线拟合,可得到系数：

$\begin{array}{l} c = - 0.1371 r^{2} + 2.5061 r - 18.2290, \\ d = 0.0060 r^{2} - 0.1124 r + 1.3456, \end{array}}$

(24)

因此,深埋隧洞涌水量解析公式表达为

$\begin{array}{l} Q = 2 π k_{m} \frac{h^{'}}{\ln (2 h^{'} / r)}, \\ h^{'} = (c r / h + d) h, \end{array}}$

(25)

式中,c,d见式（24）。

图9中相关性系数R²均大于0.95,显示了解析解与数值解具有很高的拟合度,说明了解析公式的可行性。

（3）衬砌等效地下水位

根据衬砌隧洞涌水量数值计算结果,反算出衬砌隧洞等效地下水位,绘制出衬砌隧洞等效地下水位与初始地下水位之比 ${h^{'}}_{l} / h$ 关于 $r / h$ 之间的散点图,见图10。可以看出：总体来看,衬砌隧洞等效地下水位与初始地下水位之比 ${h^{'}}_{l} / h$ 处于0.7～1.0之间,且随 $r / h$ 的增大而逐渐减少；在隧洞相对埋深 ${h^{'}}_{l} / h$ 一定时,值随围岩与衬砌相对渗透性 $k_{m} / k_{l}$ 的增大而逐渐增大,并逐渐趋向于1,这说明了地下水位降深随的增大而逐渐减少,逐渐趋向于0,从而导致衬砌隧洞上方的等效地下水位逐渐抬升。

图 10

$h^{'} / h$ 关于

$r / h$ 的散点分布及拟合曲线

Figure 10. Scatter distribution and fitting curves of

$h^{'} / h$ with respect to

$r / h$

下载: 全尺寸图片幻灯片

针对衬砌隧洞 ${h^{'}}_{l} / h$ 值的分布规律,将定义为关于 $r / h$ 的幂函数,考虑到 $k_{m} / k_{l}$ 工况下, ${h^{'}}_{l} / h$ 趋向于1,本文分别给出 $k_{m} / k_{l}$ =100与 $k_{m} / k_{l}$ =1000两种工况下衬砌隧洞等效地下水位的表达式：

${h^{'}}_{l} = {\begin{cases} e {(r / h)}^{f} \cdot h (k_{m} / k_{l} = 100) \\ w {(r / h)}^{t} \cdot h (k_{m} / k_{l} = 1000) \end{cases}$ ,

(26)

式中, $e$ , $f$ , $w$ 与 $t$ 分别为待定系数。

根据图10中数据进行曲线拟合,可得到系数：

$\begin{array}{l} e = - 0.0295 r + 0.8944, f = - 0.0104 r - 0.0163, \\ w = - 0.0052 r + 0.9980, t = - 0.0015 r + 0.0009 。 \end{array}}$

(27)

综合上式,衬砌隧洞涌水量解析公式表达为

$\begin{array}{l} Q = Q_{l} = 2 π k_{l} \frac{h_{l}}{\ln (r_{l} / r)}, \\ h_{l} = \frac{{h^{'}}_{l}}{1 + A \times k_{l} / k_{m}}, \end{array}}$

(28)

式中, ${h^{'}}_{l}$ 见式（26）与（27）。

图10中相关性系数R²均大于0.95,显示了解析解与数值解具有很高的拟合度,说明了解析公式的可行性。

3.2 孔压分布

利用等效地下水位 $h^{'}$ 与衬砌等效地下水位 ${h^{'}}_{l}$ 的表达式,计算出隧洞中心右侧水平向与顶拱垂向的衬砌水压力分布曲线,以隧洞r=5 m,h=200 m为例,将解析解、数值解及Femandez公式三者对比,见图11。可以看出：Femandez公式计算结果最大,解析解与数值解次之,其中水平向上,在距隧洞边界100 m范围内,两者拟合程度较高,之后解析解则逐渐小于数值解；垂直向上,解析解位于数值解与Femandez公式之间,且随着相对渗透性 $k_{m} / k_{l}$ 的增大,解析解与数值解的拟合程度增加。这是由于Femandez公式适用于无限补给条件下衬砌水压力的求解,不考虑水头损失,故计算值最大；解析解中等效地下水位是隧洞整体渗流场的综合体现,其值位于降落地下水位与初始地下水位之间,处于合理取值区间。

图 11 不同计算方法下隧洞中心右侧水平向与顶拱垂向的水压力分布曲线

Figure 11. Distribution of pore water pressure along horizontal and vertical directions of tunnel by different methods

下载: 全尺寸图片幻灯片

为了研究疏干条件下衬砌外水压力 $p_{l}$ 与远端孔隙水压力 $p_{h}$ 之间的变化关系,绘制出衬砌外水压力的折减系数（ $β = p_{l} / p_{h}$ ）与相对渗透性的关系曲线,见图12。可以看出：围岩与衬砌相对渗透性越大,衬砌水压力的折减系数越大,其中隧洞r=5 m, $k_{m} / k_{l}$ 由1增加至10000时, $β$ 由0.05增加至0.998；当相对渗透性 $k_{m} / k_{l}$ 和衬砌厚度一定时,折减系数随隧洞半径的增大而减小,以 $k_{m} / k_{l}$ =100为例, $β$ 由0.915（r=2.5 m）递减至0.691（r=10.0 m）。

图 12 折减系数与相对渗透性之间的关系曲线

Figure 12. Relationship between discount coefficient and k_m/k_l

下载: 全尺寸图片幻灯片

3.3 结果验证

为验证半解析公式的正确性,将半解析公式、数值解与解析公式（Goodman公式、苏凯公式、EI Tani公式等）三者计算结果相互对比,以r=5 m为例,不同计算方法下隧洞涌水量结果对比见图13。

图 13 不同计算方法隧洞涌水量结果对比

Figure 13. Comparison of water inflows into tunnel by different methods

下载: 全尺寸图片幻灯片

可以看出：隧洞涌水量随初始地下水位h的增大均逐渐增大,图13（a）中疏干隧洞Goodman公式计算结果最大,苏凯公式次之,EI Tani公式最小,而数值解与解析解处于苏凯公式与EI Tani公式之间,且两者拟合程度较高。这是由于,Goodman公式不考虑地下水位下降,因此计算结果值最大；苏凯公式与EI Tani公式均考虑到地下水位的下降,但苏凯公式仅仅假定隧洞周围渗流场处于似稳定状态,且其隧洞边界条件设为零孔隙水压力边界,造成计算结果相比较实际稳定状态偏大；而EI Tani公式在推导过程中仅考虑到单侧定水头补给,造成其计算值相对偏小。图13（b）中衬砌隧洞Femandez公式计算结果相对较大,而解析解与数值解两者结果略小,且围岩与衬砌相对渗透性越大,三者越接近。因此,经过与传统解析公式和数值解的对比,提出的解析公式更加真实刻画了隧洞渗流场特征,计算结果也更合理与准确。

4. 结论

针对浅埋、深埋及衬砌隧洞分别进行解析分析,引入关键参数（影响距离d、等效地下水位 $h^{'}$ 及衬砌等效地下水位 ${h^{'}}_{l}$ ）来表征渗流场,利用数值计算确定关键参数,给出了更切合实际的隧洞涌水量计算的半解析公式,同时探究了孔隙水压力的分布规律,得出以下结论：

（1）疏干条件下隧洞涌水量取决于隧洞半径、初始地下水位高度及围岩与衬砌相对渗透性。未衬砌隧洞涌水量随隧洞半径与初始地下水位高度之比的增大而逐渐减少；衬砌圈的存在显著改变了隧洞涌水量和孔隙水压力的分布规律,围岩与衬砌相对渗透性是决定涌水量大小的主要因素。

（2）浅埋隧洞d/D值介于0.1～0.4之间,关于呈线性分布；深埋隧洞 ${h^{'}}_{l} / h^{'}$ 值介于0.3～1.0之间,关于r/h呈线性分布；衬砌隧洞 ${h^{'}}_{l} / h^{'}$ 值介于0.7～1.0之间,与围岩与衬砌相对渗透性成正比,隧洞半径成反比,且关于r/h呈幂函数分布。

（3）利用数值计算结果和回归拟合方法,确定关键参数（影响距离d,等效地下水位 $h^{'}$ 及衬砌等效地下水位 ${h^{'}}_{l}$ ）的计算公式,提出了计算浅埋、深埋及衬砌隧洞涌水量的半解析公式,计算结果与数值解吻合较好。

（4）等效地下水位是隧洞整体渗流场的综合体现,其值位于降落地下水位与初始地下水位之间,且随着围岩与衬砌相对渗透性增大,孔压解析解与数值解的拟合程度逐渐提高。最后,将提出的解析公式、数值解与传统解析公式三者相互对比,验证了解析公式的合理性与可行性。

图 1 磨料水射流切割模型

Figure 1. Cutting model of abrasive water jet

下载: 全尺寸图片幻灯片

图 2 低横移速度下水射流冲击混凝土仿真

Figure 2. Simulation of water jet impacting concrete at low transverse velocity

下载: 全尺寸图片幻灯片

图 3 高横移速度下水射流冲击混凝土仿真

Figure 3. Simulation of water jet impacting concrete at high transverse velocity

下载: 全尺寸图片幻灯片

图 4 低横移速度下水射流冲击钢筋仿真

Figure 4. Simulation of water jet impacting steel bar at low transverse velocity

下载: 全尺寸图片幻灯片

图 5 高横移速度下水射流冲击钢筋仿真

Figure 5. Simulation of water jet impacting steel bar at high transverse velocity

下载: 全尺寸图片幻灯片

图 6 高压射流切割机床

Figure 6. High pressure jet cutting machine

下载: 全尺寸图片幻灯片

图 7 切割口损伤状态

Figure 7. Damage patterns of incision

下载: 全尺寸图片幻灯片

图 8 切割口宽度

Figure 8. Widths of incision

下载: 全尺寸图片幻灯片

图 9 不同横移速度下素混凝土切割深度

Figure 9. Cutting depths of concrete at different transverse velocities

下载: 全尺寸图片幻灯片

图 10 不同横移速度下钢筋切割深度

Figure 10. Cutting depths of steel bar under different traverse speeds

下载: 全尺寸图片幻灯片

图 11 钢筋混凝土内部切割形状

Figure 11. Internal cutting shape of reinforced concrete

下载: 全尺寸图片幻灯片

图 12 10 mm/min横移速度下钢筋混凝土切割深度

Figure 12. Cutting depths of reinforced concrete at traverse speed of 10 mm/min

下载: 全尺寸图片幻灯片

图 13 小型盾构机试验平台

Figure 13. Experimental platform of small shield machine

下载: 全尺寸图片幻灯片

图 14 刀盘

Figure 14. Cutter head

下载: 全尺寸图片幻灯片

图 15 射流切割后钢筋混凝土状态

Figure 15. Status of reinforced concrete cut by jet

下载: 全尺寸图片幻灯片

图 16 钢筋破坏形式

Figure 16. Failure forms of steel bar

下载: 全尺寸图片幻灯片

图 17 钢筋混凝土切割面

Figure 17. Cutting surface of reinforced concrete

下载: 全尺寸图片幻灯片

图 18 混凝土碎渣

Figure 18. Slag of concrete

下载: 全尺寸图片幻灯片

图 19 刀具损伤

Figure 19. Damage of cutter

下载: 全尺寸图片幻灯片

表 1 RHT模型关键参数

Table 1 Key parameters of RHT model

本构模型参数					失效参数
A	N	$P_{spall}^{*}$	Q₀	B_Q	D₁	D₂
1.6	0.61	0.1	0.6805	0.0105	0.04	1

下载: 导出CSV

表 2 J-C模型关键参数^[17]

Table 2 Key parameters of J-C model^[17]

本构模型参数					失效参数
A/MPa	B/MPa	n	C	m	D₁	D₂	D₃	D₄
507	320	0.28	0.064	1.06	0.1	0.76	1.57	-0.84

下载: 导出CSV

表 3 材料的力学参数

Table 3 Mechanical parameters of materials

类别	弹性模量/GPa	泊松比	密度/(g·cm^-3)
混凝土	35.5	0.2	2.75
钢筋	200	0.25	7.89

下载: 导出CSV

表 4 试验参数

Table 4 Experimental parameters

喷嘴材料	最大流量/(L·min^-1)	磨料直径/mm	粒子占比/%
蓝宝石	3.7	0.15~0.3	5
喷嘴直径/mm	靶距/mm	射流角度/(°)	射流速度/(m·s^-1)
0.33	20	90	800

下载: 导出CSV

参考文献(18)

[1]	徐前卫, 朱合华, 马险峰, 等. 地铁盾构隧道穿越桥梁下方群桩基础的托换与除桩技术研究[J]. 岩土工程学报, 2012, 34(7): 1217-1226. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201207008.htm XU Qian-wei, ZHU He-hua, MA Xian-feng, et al. Pile underpinning and removing technology of shield tunnels crossing through group pile foundations of road bridges[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(7): 1217-1226. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201207008.htm
[2]	孙波, 肖龙鸽, 孙正阳, 等. 深圳地铁盾构穿越建筑群及切削桩基施工[J]. 隧道建设, 2015, 35(6): 571-578. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201506021.htm SUN Bo, XIAO Long-ge, SUN Zheng-yang, et al. Shield tunneling and pile cutting construction in Shenzhen metro[J]. Tunnel Construction, 2015, 35(6): 571-578. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201506021.htm
[3]	傅德明. 盾构切削混凝土模拟试验和切削桩基施工技术[J]. 隧道建设, 2014, 34(5): 472-477. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201405020.htm FU De-ming. Model test on concrete cutting directly by shield and pile foundation cutting technology[J]. Tunnel Construction, 2014, 34(5): 472-477. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201405020.htm
[4]	袁大军, 王飞. 盾构切削大直径钢筋混凝土群桩的理论和实践[M]. 北京: 科学出版社, 2017. YUAN Da-jun, WANG Fei. Theory and Practice of Shield Cutting Large Diameter Reinforced Concrete Pile Group[M]. Beijing: Science Press, 2017. (in Chinese)
[5]	王飞, 袁大军, 董朝文, 等. 盾构直接切削大直径钢筋混凝土桩基试验研究[J]. 岩石力学与工程学报, 2013, 32(12): 2566-2574. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201312025.htm WANG Fei, YUAN Da-jun, DONG Chao-wen, et al. Test study of shiedcutting largediameter reinforced concrete piles directly[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(12): 2566-2574. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201312025.htm
[6]	弓永军. 磨料水射流切割技术研究现状及其发展趋势[J]. 液压与气动, 2016(10): 1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-YYYQ201610001.htm GONG Yong-jun. Research status and development trend of abrasive water jet cutting technology[J]. Hydraulic and Pneumatic, 2016(10): 1-5. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YYYQ201610001.htm
[7]	HOOD M. Cutting strong rock with a drag bit assisted by high-pressure water jets[J]. Africa Institute Mining Metal, 1976, 177(4): 43-54.
[8]	神山守. 地下建設における２液混合噴流の地盤掘削特性および地中支障物の切削性に関する研究[D]. 东京: 日本早稲田大学, 2018. KAMIYAMA MAMORU. Study on the Characteristics of Dual-Fluid Jets for the Soil Excavation and the Obstacles Cutting in Underground Construction[D]. Tokyo: Akita University, 2018. (in Japanese)
[9]	卢义玉, 连红军, 卢建中. 水射流辅助刀具破碎岩石机理及射流位置的研究[J]. 矿山机械, 2005, 33(11): 6-8. https://www.cnki.com.cn/Article/CJFDTOTAL-KSJX200511000.htm LU Yi-jun, LIAN Hong-jun, LU Jian-zhong. Study on the mechanism of rock breaking and the position of water jet assisted cutter[J]. Mining Machinery, 2005, 33(11): 6-8. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KSJX200511000.htm
[10]	卢义玉, 陆朝晖, 李晓红, 等. 水射流辅助PDC刀具切割岩石的力学分析[J]. 岩土力学, 2008, 29(11): 3037-3040. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200811030.htm LU Yi-jun, LU Zhao-hui, LI Xiao-hong, et al. Mechanical analysis of water jets assisting PDC bit to cut rocks[J]. Rock and Soil Mechanics, 2008, 29(11): 3037-3040. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200811030.htm
[11]	张文华, 汪志明, 于军泉, 等. 高压水射流-机械齿联合破岩数值模拟研究[J]. 岩石力学与工程学报, 2005, 24(23): 4373-4382. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200523031.htm WANG Fei, YUAN Da-jun, DONG Chao-wen, et al. Numerical simulation for combined breaking rock with high pressure water jet and mechanical bit[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(23): 4373-4382. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200523031.htm
[12]	曾锐. 滚筒式采煤机-高压水射流联合截割系统的设计与研究[D]. 徐州: 中国矿业大学, 2018. ZENG Rui. Design and Research on Combined Cutting System of Drum Shearer & High Pressure Water Jet[D]. Xuzhou: China Mining University, 2018. (in Chinese)
[13]	陈跃强. 磨料水射流-截齿联合破岩性能研究[D]. 徐州: 中国矿业大学, 2018. CHEN Yue-qiang. Research on Rock Breaking Performance of Pick Assisted with Abrasive Water Jet[D]. Xuzhou: China Mining University, 2018. (in Chinese)
[14]	杨清文, 王晓敏. 前混合磨料水射流切割钢板和混凝土的实验研究[J]. 兵工学报, 2005, 26(1): 133-135. https://www.cnki.com.cn/Article/CJFDTOTAL-BIGO200501031.htm YANG Wen-qing, WANG Xiao-min. Experimental study on cutting steel plate and concrete with premixed abrasive water jet[J]. Acta Armamentarii, 2005, 26(1): 133-135. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BIGO200501031.htm
[15]	王伟. 高压磨料水射流切割碳纤维复合材料的试验研究[D]. 哈尔滨: 哈尔滨理工大学, 2015. WANG Wei. Study on the Experiment of Carbon Fiber Composite Materials with High-Pressure Abrasive Water Jet Cutting[D]. Harbin: Harbin Institute of Ice Technology, 2015. (in Chinese)
[16]	赵宏伟. 磨料水射流切割微晶复合材料的试验研究[D]. 哈尔滨: 哈尔滨理工大学, 2017. ZHAO Hong-wei. Study on the Experiment of Microcrystalline Composite Materials with Abrasive Water Jet Cutting[D]. Harbin: Harbin Institute of Ice Technology, 2017. (in Chinese)
[17]	王哲, 吴淑伟, 姚王晶, 等. 盾构穿越既有桥梁桩基磨桩技术的研究[J]. 岩土工程学报, 2020, 42(1): 117-125. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202001019.htm WANG Zhe, WU Shu-wei, YAO Wang-jing. Grinding pile technology of shield tunnels crosssing pile foundation of existing bridges[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(1): 117-125. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202001019.htm
[18]	张文超, 武美萍, 任仲贺. 基于LS-DYNA仿真的射流加工参数分析[J]. 表面技术, 2017, 46(10): 268-276. https://www.cnki.com.cn/Article/CJFDTOTAL-BMJS201710042.htm ZHANG Wen-chao, WU Mei-ping, REN Zhong-he. Analysis of jet flow machining parameters based on LS-DYNA simulation[J]. Surface Technology, 2017, 46(10): 268-276. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BMJS201710042.htm

施引文献

资源附件(0)

图(19) / 表(4)

计量

文章访问数: 312
HTML全文浏览量: 17
PDF下载量: 193
被引次数: 0

0. 引言
1. 解析分析
1.1 浅埋隧洞
1.2 深埋隧洞
1.3 衬砌隧洞
2. 数值模拟研究
2.1 计算方法与数值模型
2.2 隧洞涌水量
2.3 孔隙水压力分布
3. 讨论与验证
3.1 涌水量解析公式
3.2 孔压分布
3.3 结果验证
4. 结论

射流联合盾构切削钢筋混凝土仿真与试验 English Version

作者简介: 庄欠伟（1978— ）,男,教授级高级工程师,硕士生导师,主要从事地下空间与工程方面教学和研究工作。E-mail：qianweizhuang@sohu.com

通讯作者: 袁一翔， E-mail：yyxyzdd@163.com

计量

出版历程

Simulation and experiment on cutting reinforced concrete with jet combined shield method

0. 引言

1. 解析分析

1.1 浅埋隧洞

1.2 深埋隧洞

1.3 衬砌隧洞

2. 数值模拟研究

2.1 计算方法与数值模型

2.2 隧洞涌水量

2.3 孔隙水压力分布

3. 讨论与验证

3.1 涌水量解析公式

3.2 孔压分布

3.3 结果验证

4. 结论

计量

出版历程

目录

0. 引言

1. 解析分析

1.1 浅埋隧洞

1.2 深埋隧洞

1.3 衬砌隧洞

2. 数值模拟研究

2.1 计算方法与数值模型

2.2 隧洞涌水量

2.3 孔隙水压力分布

3. 讨论与验证

3.1 涌水量解析公式

3.2 孔压分布

3.3 结果验证

4. 结论

作者简介:
庄欠伟（1978— ）,男,教授级高级工程师,硕士生导师,主要从事地下空间与工程方面教学和研究工作。E-mail：qianweizhuang@sohu.com

通讯作者:
袁一翔， E-mail：yyxyzdd@163.com