Study on the Influence Law of Pre-cutting Groove on the Crack Expansion of Rock Breaking by Disc Cutter
-
摘要: 针对硬岩地质条件下传统滚刀破岩出现贯入度低、磨损严重和破岩效率低等问题，许多新型辅助破岩方法被提出，如通过激光、水射流等非接触式破岩方法在岩石上制造预切缝来辅助滚刀破岩。为了明晰预切缝条件下滚刀破岩机理，采用颗粒离散元方法建立预切缝在滚刀侧面（侧置式）和滚刀前面（前置式）两种辅助破岩仿真模型，探究了预切缝深度和间距对滚刀破岩裂纹拓展的影响规律。研究结果表明：1）侧置式辅助破岩模型中，预切缝主要影响裂纹的拓展，诱导破岩过程中侧向裂纹向预切缝底端扩展；前置式辅助破岩模型中，预切缝主要影响裂纹的萌生起裂，在贯入初期只产生侧向裂纹，并促使侧向裂纹向两侧扩展。2）滚刀垂直力和裂纹数量均随着预切缝深度的增加而减小，但滚刀垂直力随着刀间距的增大而增大，裂纹数量则随着刀间距的增加先增大后减小。3）同等条件下，侧置式辅助破岩模型的破岩效率要高于前置式辅助破岩模型，但前置式辅助破岩模型中滚刀的垂直力要远小于侧置式辅助破岩模型中的滚刀垂直力，可通过进一步增加前置式辅助破岩模型中滚刀贯入度来充分利用预切缝辅助破岩。4）滚刀贯入度一定时，侧置式辅助破岩模型适用于预切缝深度大于滚刀贯入度的情况，而前置式辅助破岩模型适用于预切缝深度小于滚刀贯入度的情况。Abstract: To solve the problems of low penetration, serious wear and low rock-breaking efficiency during traditional disc cutter rock breaking under the geological conditions of hard rock, many new auxiliary rock breaking methods have been proposed, such as creating pre-cutting grooves on the rock through laser, water jet, and other non-contact rock-breaking methods to assist disc cutter rock breaking. To clarify the rock breaking mechanism of the disc cutter under the condition of the pre-cutting grooves, the particle discrete element method is used to establish two auxiliary rock breaking simulation models for the pre-cutting grooves at the side of the disc cutter (side-type) and in front of the disc cutter (front-type). The influences of the depth of the pre-cutting grooves and the cutter spacing on the propagation of rock-breaking cracks by the disc cutter are explored. The research results show that: 1) The pre-cutting grooves in the side-type auxiliary rock-breaking model mainly affect the expansion of cracks, and the lateral cracks propagate to the bottom of the pre-cutting grooves during the process of inducing rock breaking. The pre-cutting grooves in the front-type auxiliary rock-breaking model mainly affect the initiation of cracks, and only lateral cracks are generated at the early stage of penetration, and the lateral cracks are promoted to propagate to both sides; 2) The vertical force of the disc cutter and the number of cracks decreases with the increase of the depth of the pre-cutting grooves, while the vertical force of the disc cutter increases with the increase of the cutter spacing, and the number of cracks first increases and then decrease with the increase of the cutter spacing; 3) Under the same condition, the side-type auxiliary rock-breaking model is higher than that of the front-type auxiliary rock-breaking model, but the vertical force of the disc cutter in the front-type auxiliary rock-breaking model is much smaller than the vertical force of the disc cutter in the side-type auxiliary rock-breaking model. The auxiliary effect of the pre-cutting grooves can be fully utilized by further increasing the penetration of the disc cutter in the front-type auxiliary rock-breaking model; 4) When the penetration of the disc cutter is constant, the side-type auxiliary rock-breaking model is suitable for the case where the depth of the pre-cutting grooves is greater than the penetration of the disc cutter, and the front-type auxiliary rock-breaking model is suitable for the case where the depth of the pre-cutting grooves is less than the penetration of the disc cutter.
-
Keywords:
- pre-cutting grooves /
- crack extension /
- disc cutter /
- particle discrete element /
- rock fragmentation
-
全断面隧道掘进机（tunnel boring machine，TBM）以其安全、可靠、环保等优点，被广泛应用于水利水电涵洞、采矿、铁路公路隧道施工等地下空间工程[1]。在硬岩应用中，最常用的刀具是滚刀，它已成为硬岩TBM的标准刀具[2]。然而，随着地下空间不断开发利用，面对全断面硬岩隧道，依然会遇到很多问题，如贯入度低、滚刀磨损率高、掘进效率低等[3]，而更换刀具意味着花费更多的时间和施工成本。据统计，在高磨损地层中更换和修复TBM刀具的时间成本为掘进总时间的1/3，而刀具成本约为挖掘施工成本的1/3[4]。
为了减少滚刀磨损，提高掘进效率，学者们从材料、磨损特性[5-6]、破岩机理以及滚刀的结构设计[7-8]等方面进行了大量研究，并提出了多种解决方案。虽然上述方法取得了较好的效果，但硬岩地质条件下TBM挖掘效率低、滚刀磨损率高的问题仍未得到有效解决。很多学者开始探究预切缝辅助破岩方法，即利用激光[9]、水射流[10-11]等非接触式破岩方法来制造预切缝，降低岩石围压从而提升滚刀破岩效率。2019年6月在中国开发了第1个集成高压水射流系统的隧道掘进机[12]，在掘进过程中，水射流和滚刀同时作用在掌子面上，水射流的主要作用是在掌子面切割形成预切缝，破坏掌子面岩石的完整性，从而辅助滚刀破岩。
关于预切缝辅助滚刀破岩效果，现已有学者进行了部分研究。Fenn[13]研究了一种典型的水射流预切缝辅助切割系统，使得盘形滚刀受到的垂直力和滚动力减少40%。Ciccu等[14]将水射流置于滚刀一侧进行辅助开挖，发现射流对滚刀一侧具有削弱作用，从而提高了破岩体积。Wang等[15]基于数值模型研究了水射流参数对辅助滚刀破岩性能的影响，其结果表明，最佳水射流喷嘴直径为1.5 mm，最佳水射流喷射速度为224.5 m/s。Cheng等[16]研究了预切缝参数对岩石破碎力学行为的影响，其结果表明，峰值穿透力随切口深度的增加而减小，随切口间距的增加而略有增加。徐福通等[17]通过滚刀压头贯入试验研究预切缝对TBM机械滚刀破岩效果的影响规律，其结果表明，与无切槽相比，有预切槽的岩样破碎时法向荷载降低 44.13%～53.05%。
虽然预切缝辅助滚刀破岩技术正在应用到TBM硬岩掘进工程中，但相应的理论研究相对较少，现有的研究大多是对辅助破岩的可行性和效果，以及破岩参数的影响进行探究，鲜有人注意到在硬岩地质条件下预切缝辅助滚刀破岩裂纹扩展机理。本文从岩石细观破坏出发，采用颗粒离散元方法[18]建立了预切缝下滚刀破岩数值模型，分析预切缝下滚刀破岩裂纹扩展机理，并研究不同预切缝深度和刀间距对预切缝辅助滚刀破岩裂纹扩展的影响。
1. 数值模型
1.1 颗粒离散元方法简介
近年来，颗粒离散元方法在岩土工程中得到了广泛的发展和应用[18]。它可以解决传统岩土理论无法解释的典型力学现象，揭示微观介质在复杂岩土介质条件下的累积损伤和破坏机理。颗粒离散元方法的计算原理主要是通过颗粒与颗粒之间运用牛顿运动定律进行反复计算，得到颗粒与颗粒之间的接触力，同时更新相应颗粒的速度、加速度以及位移；当颗粒得到了新参数后进行重新迭代，颗粒与颗粒之间的接触力将被重新赋值，通过迭代更新，不断进行计算，直到达到所要求的平衡条件为止，如图1所示。
图 1 颗粒离散元方法计算原理Fig. 1 Calculation principle of particle discrete element method
下载:
全尺寸图片
在颗粒离散元数值模拟中，岩石材料是由接触点黏结在一起的刚性圆盘组成。颗粒之间的接触模型一般有接触黏结模型、平行黏结模型和滑移模型。其中平行黏结模型可以传递力和力矩，黏结键作用于颗粒之间的圆形或矩形截面上，刚度由接触刚度和键合刚度共同决定。因此，黏结破坏立即导致刚度降低。从这个意义上说，平行黏结模型是一种更为真实的岩石类材料模型，在这种模型中，键可以在张力或剪切中断裂，从而降低刚度[19]。所以本文中选用线性平行连接模型定义颗粒间的接触。
1.2 岩石细观参数标定
在颗粒流离散元模型中，岩石的宏观力学性能是由微观颗粒的性质决定的，但两者之间并没有一个清晰的对应关系。为了在仿真中较为真实地模拟出实际岩石的宏观力学性能，除了选择合适的本构模型外，还需要对细观颗粒参数进行标定。本文中采用的岩样是从某项目现场开采的岩石强度较高、完整性较好的花岗岩。本文通过单轴压缩和巴西劈裂室内实验获得岩石的宏观力学参数，其力学参数如表1所示。
表 1 花岗岩试样宏观力学参数Table 1 Macro-mechanical parameters of granite samples抗压强度/MPa 抗拉强度/MPa 弹性模量/GPa 泊松比 密度/
(g·cm–3)156.2 26.5 62.5 0.27 2.97 为了匹配材料的宏观力学性能，建立单轴压缩和巴西劈裂模拟模型对岩石颗粒细观参数进行标定，采用试错法[20]不断地调试颗粒的微观力学参数，直至模拟得到岩石模型的宏观力学指标与实际岩石指标相近，如图2所示。经过反复的标定，岩石颗粒的细观参数如表2所示。岩石颗粒细观参数中平行键黏结部分和线性黏结部分与之相同。
图 2 颗粒流模型细观参数的确定Fig. 2 Determination of meso-parameters of particle flow model下载:
全尺寸图片
表 2 花岗岩试样微观参数Table 2 Microscopic parameters of granite samples颗粒接触参数 数值 颗粒最小半径/mm 0.5 颗粒粒径比 1.6 有效模量/GPa 35 刚度比 2.5 摩擦系数 0.5 黏结法向强度/MPa 50 黏结切向强度/MPa 100 内摩擦角/(°) 0 1.3 预切缝辅助滚刀破岩数值模型建立
通过在岩石模型上预先设置切缝，重点分析预切缝条件下滚刀垂直力所主导的裂纹产生与扩展过程，在此情况下，将滚刀的3维的滚压运动简化为2维的侵入问题是可行的[21-22]。本文以17寸盘形滚刀（即外径为432 mm的滚刀）为研究对象，确定其主要几何参数：滚刀刀刃顶部宽度为12 mm，过渡圆弧半径为5 mm，刀刃角为20°。将滚刀刀圈视为刚性体，以颗粒离散元中的刚性墙体模拟。目前，预切缝辅助滚刀破岩主要有两种方式[11]：一种是侧置式，即预切缝在滚刀侧面；另一种是前置式，即预切缝在滚刀的前方。采用上述标定后的微观参数分别建立了侧置式模型和前置式模型，如图3所示。整个模型宽400 mm，高200 mm，包含颗粒数目为54492。岩体除开挖掌子面外，使用刚性墙边界固定。滚刀的起始位置与岩石上表面的间距为0.1 mm，通过给滚刀施加垂直向下大小为0.1 mm/s的速度，使得滚刀匀速向下侵入岩石，通过设定滚刀最大竖向位移来控制滚刀贯入度。预切缝宽度设置为实验室实验常见平均宽度5 mm，切缝深度为h，滚刀刀间距为S，滚刀贯入度为P。
图 3 预切缝辅助滚刀破岩模型Fig. 3 Rock breaking model with disc cutter assisted by pre-cutting grooves下载:
全尺寸图片
1.4 预切缝辅助滚刀破岩数值模型实验验证
预切缝辅助滚刀破岩实验与模拟对比如图4所示。由图4可知，实验和模拟所得到的裂纹拓展方向较为相似，即当P为6 mm、S为120 mm时，滚刀破岩产生的裂纹朝着预切缝扩展并贯通，促使岩石开始起裂破碎，产生相应的岩渣。
图 4 预切缝辅助滚刀破岩实验与模拟对比Fig. 4 Experimental and simulation comparison of pre-cutting grooves assisted disc cutter rock breaking下载:
全尺寸图片
将仿真模拟中的滚刀垂直力转换为3维情况，从而与实验得到的滚刀垂直力进行对比。根据文献[23]得到如下转换公式：
$$ {F_{\rm{N}}} = 2\int_0^\beta {\rm{d}} {F_{\rm{N}}} = T_{\theta} \cos (\beta - \theta )r{\rm{d}}\theta = \frac{{2PT(\beta )}}{{\arccos (1 - P/r)}} $$ (1) 式中，FN为垂直力，β为滚刀刃半接岩角，θ为积分段与外侧半径夹角，Tθ为滚刀刀刃压力，T(β)为滚刀刀刃对称平面上的压力，r为滚刀半径。
根据式（1）将仿真结果转化后得到滚刀的垂直力平均值约为123.2 kN；实验中滚刀垂直力随时间变化曲线如图5所示。实验得到滚刀垂直力平均值约为158.6 kN；实验中滚刀垂直力和仿真垂直力比值大约为22.3%，这可能是由仿真模型和实验边界条件差别造成；实验中岩石厚度远大于仿真模型岩石厚度，从而提供了更强的约束，导致破岩载荷偏大。
图 5 预切缝辅助滚刀破岩垂直力变化曲线Fig. 5 Variation curve of normal force of pre-cutting grooves assisted disc cutter rock breaking下载:
全尺寸图片
从裂纹拓展方向和滚刀垂直力两方面来说，仿真结果和实验结果是相似的，说明了预切缝辅助滚刀破岩模型是可靠的。
2. 仿真结果和分析
2.1 不同预切缝下滚刀破岩裂纹扩展分析
为了分析两种破岩方式下预切缝对滚刀破岩的影响，以滚刀间距S为120 mm、切缝深度h为9 mm、贯入度P为9 mm的仿真结果为例，与相同参数设置下无预切缝仿真结果进行对比分析。
2.1.1 破岩过程分析
图6为无预切缝滚刀侵入岩体裂纹生成与扩展。由图6可知：无预切缝常规破岩时，刀具侵入岩体前期（贯入度0.5 mm以内），岩体主要发生弹性变形，几乎无裂纹产生。当贯入到0.5 mm时（图6（a）），滚刀下方形成微裂纹，其中，主要以张拉裂纹为主，伴有极少量剪切裂纹产生（红色裂纹表示张拉裂纹，绿色裂纹表示剪切裂纹）。当继续增加贯入度到3 mm时，产生大量的裂纹，出现宏观径向裂纹和侧向裂纹起裂并扩展，且此时的裂纹以张拉裂纹为主，即宏观裂纹的起裂和扩展以张拉破坏为主（图6（b））。随着贯入度的继续增加，径向裂纹缓慢扩展，侧向裂纹逐渐扩展到岩面（图6（c）），部分侧向裂纹与岩面贯通，产生宏观破碎。贯入度继续增加到9 mm时（图6（d）），径向裂纹和侧向裂纹进一步缓慢扩展，但两把滚刀之间依然没有裂纹交汇贯通，这是因为刀间距较大的原因。
图 6 无预切缝滚刀侵入岩体裂纹生成与扩展Fig. 6 Formation and propagation of cracks intruding into rock mass by the disc cutter without pre-cutting slot下载:
全尺寸图片
图7为侧置式模型滚刀侵入岩体裂纹生成与扩展。由图7可知：在刀具侵入岩体前期，裂纹的萌生起裂与常规破岩基本一致（图7（a）和（b））；随着贯入度的增加，径向裂纹依然沿原方向缓慢扩展，但因为两把滚刀中间预切缝的存在，部分侧向裂纹的扩展方向发生改变（图7（c）），两把滚刀产生的部分原本朝岩面扩展的侧向裂纹，朝着更为薄弱的预切缝底端扩展，并在预切缝的底端交汇贯通；并且，当贯入度逐步增加到9 mm时（图7（d）），新的侧向裂纹在预切缝的底端交汇贯通，且裂纹的弧度更大更深，形成的破碎区域也更大。
图 7 侧置式模型滚刀侵入岩体裂纹生成与扩展Fig. 7 Formation and propagation of cracks intruding into rock mass by the disc cutter of the side-mounted model下载:
全尺寸图片
图8为前置式模型滚刀侵入岩体裂纹生成与扩展。由图8可知：在刀具侵入岩体前期，只有滚刀刀刃两侧与岩体接触；在经历弹性变形阶段之后，贯入到0.5 mm时，只有少量张拉裂纹产生；当贯入度继续增加到3 mm时，出现宏观侧向裂纹起裂并扩展，无径向裂纹产生，且依然只有张拉裂纹；贯入度增加到6 mm时，有少量剪切裂纹产生，并出现宏观径向裂纹的起裂，但仍以张拉裂纹和侧向裂纹为主，部分侧向裂纹扩展至岩面；当进一步贯入到9 mm时，预切缝被填满并压实后，滚刀刀刃与岩体充分接触，径向裂纹得到进一步的起裂和扩展，而部分侧向裂纹在两把滚刀中间交汇贯通，岩石起裂破碎。
图 8 前置式模型滚刀侵入岩体裂纹生成与扩展Fig. 8 Formation and propagation of cracks intruding into rock mass by the disc cutter of the front-end model下载:
全尺寸图片
2.1.2 裂纹扩展特征分析
图9为常规破岩、侧置式、前置式3种情况下，贯入度从0到9 mm滚刀连续贯入仿真过程中，滚刀破岩载荷、裂纹数量与贯入度关系曲线（均以1#刀载荷为例）。由图9可以观察到：3种情况下滚刀所受到的侧向力变化趋势一致，均在0附近波动，且波动越来越明显；常规滚刀破岩和侧置式模型破岩垂直力变化趋势一致，且与Cheng等[16]研究中滚刀侵入实验载荷波动情况一致，在贯入初期均有一个最大峰值载荷出现，并伴随裂纹数量的快速增长，大部分裂纹在此阶段萌生起裂；随后，载荷降低到一个较低水平，载荷呈现上下波动的状态，体现了岩石的跃进破碎特性，且随着贯入度的增加，波动峰值越来越大。
图 9 切削力和裂纹数随贯入度变化曲线Fig. 9 Cutting force and crack number change curves with penetration下载:
全尺寸图片
前置式模型的垂直力与前两种情况有较大的差异，由于滚刀下方预切缝的存在，滚刀贯入初期并未出现一个较大峰值载荷，载荷从0开始上下波动，垂直力载荷一直处于缓慢增长较低的水平，这也是这阶段裂纹萌生和扩展较慢的原因。随着贯入度的增加，垂直力载荷增大，其波动幅度也越来越大，裂纹数量也呈逐步增长的趋势。由此可以说明：滚刀垂直力载荷特性与裂纹的生成和扩展有关；峰值载荷越大、波动越大时，裂纹的生成和扩展也越容易。
图10为这3种情况下滚刀的平均垂直力随着贯入度的变化曲线。由图10可知，3种情况下，滚刀垂直力均随着贯入度的增大而增大，但是侧置式模型和前置式模型的垂直力均明显小于无切缝的常规滚刀破岩。这也说明两种模型中，预切缝的存在能够有效减小滚刀垂直力，特别是前置式模型中，滚刀刀刃前期与预切缝直接接触，滚刀垂直力载荷大幅降低，但在滚刀下方难以形成密实核，裂纹数量较少。
图 10 不同破岩方式下垂直力随贯入度变化曲线Fig. 10 Vertical force varies with penetration under different rock breaking modes下载:
全尺寸图片
结合上述分析和第2.1.1节内容可知，常规破岩和侧置式模型滚刀贯入岩石过程中，岩石裂纹的生成和扩展可以分为4个阶段，如图9（a）和（b）所示：1）弹性变形阶段（O点至A点）。岩样发生了弹性变形，裂纹还未产生。2）裂纹快速萌生起裂、初步扩展阶段（A点至B点）。此阶段与滚刀的最大峰值载荷区间相对应，裂纹数量快速增长，剪切裂纹基本在这个阶段中产生。3）裂纹缓慢扩展、数量趋于稳定阶段（B点至C点）。4）裂纹贯通、产生宏观破坏阶段（C点以后）。侧置式模型中，预切缝的影响发生在第3）阶段和第4）阶段中，影响裂纹的扩展方向和宏观破坏程度。前置式模型只有3个阶段：①弹性变形阶段；②裂纹随贯入度增加缓慢产生并扩展阶段；③裂纹贯通、产生宏观破坏阶段。
2.2 预切缝深度对裂纹扩展的影响
在预切缝辅助机械刀具破岩的研究中，切缝深度一直被研究者视为决定辅助破岩能力的重要影响因素[24]。因此研究预切缝深度对两种模型中滚刀破岩裂纹扩展的影响。设置滚刀贯入度P为6 mm，刀间距S为60 mm，预切缝深度h分别为3、6、9、12 mm，进行仿真计算，获得各预切缝深度下平均垂直力和裂纹数量，如表3所示。
表 3 各预切缝深度下垂直力和裂纹数量统计Table 3 Statistics of vertical force and number of cracks at each pre-cut depth预切缝深度/mm 侧置式模型 前置式模型 垂直力/kN 裂纹数量 垂直力/kN 裂纹数量 3 135.1 2245 112.6 1872 6 116.3 1972 90.8 1653 9 111.7 1770 67.4 1139 12 103.4 1871 56.3 649 图11和12分别是滚刀垂直力和裂纹数量随预切缝深度变化曲线。由图11和12可知，两种模型的垂直力和裂纹数量均随着预切缝深度的增加而减小，但前置式模型的垂直力和裂纹数量减小幅度要远大于侧置式模型的减小幅度。
图 11 垂直力随预切缝深度变化曲线Fig. 11 Vertical force varies with the depth of the pre-cut seam下载:
全尺寸图片
图 12 裂纹数量随预切缝深度变化曲线Fig. 12 Number of cracks varies with the depth of the pre-cut joint下载:
全尺寸图片
图13和14分别展示了侧置式模型和前置式模型在各预切缝深度下岩石破碎情况。如图13（a）～（d）所示：在各预切缝深度下两把滚刀之间均产生了贯通裂纹；当预切缝深度大于滚刀贯入度6 mm时，贯通裂纹呈弧状，从滚刀刃部扩展到预切缝底端交汇，形成片状岩碴，且岩碴大小随着预切缝深度的增大而增大；当预切缝深度小于等于6 mm时，预切缝对裂纹的扩展影响较小，贯通裂纹避开预切缝在其下方交汇。如图14（a）～（d）所示：前置式模型中当预切缝深度小于等于滚刀贯入度6 mm时，会产生贯通裂纹；当其大于6 mm时，滚刀破岩载荷和裂纹数量大幅度下降，无贯通裂纹产生，破岩量也大幅度降低，不利于滚刀破岩。由此说明滚刀贯入度一定时，侧置式模型适用于预切缝深度大于贯入度情况，而前置式模型适用于预切缝深度小于贯入度情况。
图 13 侧置式模型中各预切缝深度下岩石破碎情况Fig. 13 Rock fragmentation at each pre-cut depth in the side-mounted model下载:
全尺寸图片
图 14 前置式模型中各预切缝深度下岩石破碎情况Fig. 14 Rock fragmentation at each pre-cut depth in the front-end model下载:
全尺寸图片
2.3 刀间距对裂纹扩展的影响
刀间距是滚刀在刀盘上一个重要布置参数，工程上刀盘的刀间距一般为40～120 mm[25]。设置滚刀贯入度P为6 mm，预切缝深度h为6 mm，刀间距S分别为70、80、90、100 mm，进行仿真计算，研究刀间距对两种模型中滚刀破岩裂纹扩展的影响，仿真结果如表4所示。
表 4 各刀间距下垂直力和裂纹数量统计Table 4 Statistics of vertical force and number of cracks at each cutter spacing刀间距/mm 侧置式模型 前置式模型 垂直力/kN 裂纹数量 垂直力/kN 裂纹数量 70 122.7 2248 96.2 1847 80 125.3 2324 105.3 1594 90 132.6 2412 112.7 1631 100 141.2 1980 109.4 1621 图15和16分别是滚刀垂直力和裂纹数量随刀间距的变化曲线。由图15和16可知：两种模型的垂直力均随着刀间距的增大而增大，且前置式模型的垂直力均要比侧置式模型的小很多；裂纹数量则随着刀间距的增加呈先增大后减小的趋势。
图 15 垂直力随刀间距变化曲线Fig. 15 Vertical force varies with the cutter spacing下载:
全尺寸图片
图 16 裂纹数量随刀间距变化曲线Fig. 16 Number of cracks varies with the cutter spacing下载:
全尺寸图片
图17和18分别展示了侧置式模型和前置式模型在各刀间距下岩石破碎情况。如图17（a）～（d）所示：对于侧置式辅助破岩，当刀间距小于100 mm时，均产生了交汇贯穿裂纹，形成片状岩碴，且岩碴大小随着刀间距的增大而增大；当刀间距为100 mm时，岩石未起裂破碎，所以此时临界刀间距为90 mm。如图18（a）～（d）所示，对于前置式辅助破岩，岩石产生贯通裂纹只发生在刀间距为70 mm时，也就是说当刀间距不大于70 mm时，前置式辅助破岩下裂纹贯通，临界刀间距为70 mm。由此可见，两种模型裂纹数量呈现出先增大后减小变化趋势，拐点出现在刀间距能够产生贯通裂纹的临界点上。
图 17 侧置式模型中各刀间距下岩石破碎情况Fig. 17 Rock fragmentation under each cutter spacing in the side-mounted model下载:
全尺寸图片
图 18 前置式模型中各刀间距下岩石破碎情况Fig. 18 Rock fragmentation under each cutter spacing in the front-end model下载:
全尺寸图片
3. 结 论
基于颗粒离散元的方法，建立了两种典型预切缝辅助滚刀破岩模型，分析了不同预切缝下滚刀破岩裂纹扩展机理，并研究了预切缝深度和刀间距这两个重要参数对两种模型破岩时裂纹扩展的影响，得出如下重要结论：
1）两种模型中预切缝辅助滚刀破岩机理不同，侧置式模型中预切缝的辅助作用发生在裂纹扩展阶段，诱导破岩过程中侧向裂纹向预切缝底端扩展，并交汇贯通；前置式模型中预切缝主要影响裂纹的萌生起裂，在贯入初期只产生侧向裂纹，并促使侧向裂纹向两侧扩展。
2）贯入度一定时，侧置式模型适用于预切缝深度大于滚刀贯入度的情况，而前置式模型适用于预切缝深度小于贯入度的情况。
3）两种模型破岩时垂直力随着预切缝深度的增加而减小，随着刀间距的增加而变大；产生的总裂纹数量随着预切缝深度的增加而减小，随着刀间距的增加先变大后减小，拐点出现在刀间距能够产生贯通裂纹的临界点上；切缝深度为6 mm时，侧置式和前置式辅助下岩石起裂破碎的临界刀间距分别为90 mm和70 mm。
4）在相同参数下，侧置式模型的辅助破岩效率要高于前置式模型的辅助破岩效率，但值得注意的是，前置式模型中滚刀的垂直力要远小于侧置式模型中的滚刀垂直力，可以通过进一步增加前置式模型中滚刀的贯入度来充分利用预切缝的辅助作用。
-
图 1 颗粒离散元方法计算原理
Fig. 1 Calculation principle of particle discrete element method
下载:
全尺寸图片
图 2 颗粒流模型细观参数的确定
Fig. 2 Determination of meso-parameters of particle flow model
下载:
全尺寸图片
图 3 预切缝辅助滚刀破岩模型
Fig. 3 Rock breaking model with disc cutter assisted by pre-cutting grooves
下载:
全尺寸图片
图 4 预切缝辅助滚刀破岩实验与模拟对比
Fig. 4 Experimental and simulation comparison of pre-cutting grooves assisted disc cutter rock breaking
下载:
全尺寸图片
图 5 预切缝辅助滚刀破岩垂直力变化曲线
Fig. 5 Variation curve of normal force of pre-cutting grooves assisted disc cutter rock breaking
下载:
全尺寸图片
图 6 无预切缝滚刀侵入岩体裂纹生成与扩展
Fig. 6 Formation and propagation of cracks intruding into rock mass by the disc cutter without pre-cutting slot
下载:
全尺寸图片
图 7 侧置式模型滚刀侵入岩体裂纹生成与扩展
Fig. 7 Formation and propagation of cracks intruding into rock mass by the disc cutter of the side-mounted model
下载:
全尺寸图片
图 8 前置式模型滚刀侵入岩体裂纹生成与扩展
Fig. 8 Formation and propagation of cracks intruding into rock mass by the disc cutter of the front-end model
下载:
全尺寸图片
图 9 切削力和裂纹数随贯入度变化曲线
Fig. 9 Cutting force and crack number change curves with penetration
下载:
全尺寸图片
图 10 不同破岩方式下垂直力随贯入度变化曲线
Fig. 10 Vertical force varies with penetration under different rock breaking modes
下载:
全尺寸图片
图 11 垂直力随预切缝深度变化曲线
Fig. 11 Vertical force varies with the depth of the pre-cut seam
下载:
全尺寸图片
图 12 裂纹数量随预切缝深度变化曲线
Fig. 12 Number of cracks varies with the depth of the pre-cut joint
下载:
全尺寸图片
图 13 侧置式模型中各预切缝深度下岩石破碎情况
Fig. 13 Rock fragmentation at each pre-cut depth in the side-mounted model
下载:
全尺寸图片
图 14 前置式模型中各预切缝深度下岩石破碎情况
Fig. 14 Rock fragmentation at each pre-cut depth in the front-end model
下载:
全尺寸图片
图 15 垂直力随刀间距变化曲线
Fig. 15 Vertical force varies with the cutter spacing
下载:
全尺寸图片
图 16 裂纹数量随刀间距变化曲线
Fig. 16 Number of cracks varies with the cutter spacing
下载:
全尺寸图片
图 17 侧置式模型中各刀间距下岩石破碎情况
Fig. 17 Rock fragmentation under each cutter spacing in the side-mounted model
下载:
全尺寸图片
图 18 前置式模型中各刀间距下岩石破碎情况
Fig. 18 Rock fragmentation under each cutter spacing in the front-end model
下载:
全尺寸图片
表 1 花岗岩试样宏观力学参数
Table 1 Macro-mechanical parameters of granite samples
抗压强度/MPa 抗拉强度/MPa 弹性模量/GPa 泊松比 密度/
(g·cm–3)156.2 26.5 62.5 0.27 2.97 表 2 花岗岩试样微观参数
Table 2 Microscopic parameters of granite samples
颗粒接触参数 数值 颗粒最小半径/mm 0.5 颗粒粒径比 1.6 有效模量/GPa 35 刚度比 2.5 摩擦系数 0.5 黏结法向强度/MPa 50 黏结切向强度/MPa 100 内摩擦角/(°) 0 表 3 各预切缝深度下垂直力和裂纹数量统计
Table 3 Statistics of vertical force and number of cracks at each pre-cut depth
预切缝深度/mm 侧置式模型 前置式模型 垂直力/kN 裂纹数量 垂直力/kN 裂纹数量 3 135.1 2245 112.6 1872 6 116.3 1972 90.8 1653 9 111.7 1770 67.4 1139 12 103.4 1871 56.3 649 表 4 各刀间距下垂直力和裂纹数量统计
Table 4 Statistics of vertical force and number of cracks at each cutter spacing
刀间距/mm 侧置式模型 前置式模型 垂直力/kN 裂纹数量 垂直力/kN 裂纹数量 70 122.7 2248 96.2 1847 80 125.3 2324 105.3 1594 90 132.6 2412 112.7 1631 100 141.2 1980 109.4 1621 -
[1] Deng Mingjiang.Challenges and thoughts on risk management and control for the group construction of a super-long tunnel by TBM[J].Engineering,2018,4(1):112–122. doi: 10.1016/j.eng.2017.07.001 [2] Huo Junzhou,Sun Xiaolong,Li Guangqing,et al.Multi-degree-of-freedom coupling dynamic characteristic of TBM disc cutter under shock excitation[J].Journal of Central South University,2015,22(9):3326–3337. doi: 10.1007/s11771-015-2873-3 [3] Lin Laikuang,Xia Yimin,Mao Qingsong,et al.Experimental study on wear behaviors of TBM disc cutter ring in hard rock conditions[J].Tribology Transactions,2018,61(5):920–929. doi: 10.1080/10402004.2018.1442895 [4] Su Weilin,Li Xinggao,Jin Dalong,et al.Analysis and prediction of TBM disc cutter wear when tunneling in hard rock strata:A case study of a metro tunnel excavation in Shenzhen,China[J].Wear,2020,446/447:203190. doi: 10.1016/j.wear.2020.203190 [5] Zhang Xuhui,Lin Laikuang,Xia Yimin,et al.Experimental study on wear of TBM disc cutter rings with different kinds of hardness[J].Tunnelling and Underground Space Technology,2018,82:346–357. doi: 10.1016/j.tust.2018.08.050 [6] 刘建琴,魏琛,郭伟.基于离散元方法的TBM刀盘磨损研究及优化设计[J].铁道工程学报,2020,37(9):85–90. doi: 10.3969/j.issn.1006-2106.2020.09.015 Liu Jianqin,Wei Chen,Guo Wei.Wear research and optimization design of TBM cutter based on discrete element method[J].Journal of Railway Engineering Society,2020,37(9):85–90 doi: 10.3969/j.issn.1006-2106.2020.09.015 [7] 孙伟,郭莉,周建军,等.TBM双滚刀破岩过程模拟及刀圈结构设计[J].煤炭学报,2015,40(6):1297–1302. doi: 10.13225/j.cnki.jccs.2014.3032 Sun Wei,Guo Li,Zhou Jianjun,et al.Rock fragmentation simulation under dual TBM disc cutter and design of cutter ring[J].Journal of China Coal Society,2015,40(6):1297–1302 doi: 10.13225/j.cnki.jccs.2014.3032 [8] Ning Bo,Xia Yimin,Lin Laikuang,et al.Experimental study on the adaptability of cutters with different blade widths under hard rock and extremely hard rock conditions[J].Acta Geotechnica,2020,15(11):3283–3294. doi: 10.1007/s11440-020-00958-0 [9] Erfan M R,Shahriar K,Sharifzadeh M,et al.Coupled T-H-M processes’ effect on specific energy in continuous wave fiber laser rock perforation[J].Journal of Laser Applications,2018,30(3):032005. doi: 10.2351/1.5018799 [10] Liu Songyong,Ji Huifu,Han Dongdong,et al.Experimental investigation and application on the cutting performance of cutting head for rock cutting assisted with multi-water jets[J].The International Journal of Advanced Manufacturing Technology,2018,94(5/6/7/8):2715–2728. doi: 10.1007/s00170-017-1072-9 [11] Zhang Jinliang,Li Yongchang,Zhang Yuansheng,et al.Using a high-pressure water jet-assisted tunnel boring machine to break rock[J].Advances in Mechanical Engineering,2020,12(10):168781402096229. doi: 10.1177/1687814020962290 [12] Zheng Yanlong,He Lei.TBM tunneling in extremely hard and abrasive rocks:Problems,solutions and assisting methods[J].Journal of Central South University,2021,28(2):454–480. doi: 10.1007/s11771-021-4615-z [13] Fenn O.The use of water jets to assist free-rolling cutters in the excavation of hard rock[J].Tunnelling and Underground Space Technology,1989,4(3):409–417. doi: 10.1016/0886-7798(89)90085-0 [14] Ciccu R,Grosso B.Improvement of disc cutter performance by water jet assistance[J].Rock Mechanics and Rock Engineering,2014,47(2):733–744. doi: 10.1007/s00603-013-0433-4 [15] Wang Fengchao,Zhou Dapeng,Zhou Xin,et al.Rock breaking performance of TBM disc cutter assisted by high-pressure water jet[J].Applied Sciences,2020,10(18):6294. doi: 10.3390/app10186294 [16] Cheng Jianlong,Jiang Zihao,Han Weifeng,et al.Breakage mechanism of hard-rock penetration by TBM disc cutter after high pressure water jet precutting[J].Engineering Fracture Mechanics,2020,240:107320. doi: 10.1016/j.engfracmech.2020.107320 [17] 徐福通,卢景景,周辉,等.预切槽和TBM机械滚刀的新型联合破岩模式研究[J].岩土力学,2021,42(5):1363–1372. doi: 10.16285/j.rsm.2020.1468 Xu Futong,Lu Jingjing,Zhou Hui,et al.Research on combined rock-breaking mode of pre-cutting groove and TBM mechanical cutter[J].Rock and Soil Mechanics,2021,42(5):1363–1372 doi: 10.16285/j.rsm.2020.1468 [18] 刘学伟,刘泉声,黄诗冰,等.裂隙岩体温度–渗流耦合数值流形方法[J].四川大学学报(工程科学版),2013,45(增刊2):77–83. Liu Xuewei,Liu Quansheng,Huang Shibing,et al.Study on numerical manifold method of coupled thermo-hydraulic of fractured rock masses[J].Journal of Sichuan University(Engineering Science Edition),2013,45(Supp2):77–83 [19] 李勇,蔡卫兵,朱维申,等.水力耦合作用下裂纹扩展演化机理的试验和颗粒流分析[J].工程科学与技术,2020,52(3):21–31. doi: 10.15961/j.jsuese.201900652 Li Yong,Cai Weibing,Zhu Weishen,et al.Experiment and particle flow analysis of crack propagation evolution mechanism under hydraulic coupling[J].Advanced Engineering Sciences,2020,52(3):21–31 doi: 10.15961/j.jsuese.201900652 [20] Lee H,Jeon S.An experimental and numerical study of fracture coalescence in pre-cracked specimens under uniaxial compression[J].International Journal of Solids and Structures,2011,48(6):979–999. doi: 10.1016/j.ijsolstr.2010.12.001 [21] Gong Q M,Jiao Y Y,Zhao J.Numerical modelling of the effects of joint spacing on rock fragmentation by TBM cutters[J].Tunnelling and Underground Space Technology,2006,21(1):46–55. doi: 10.1016/j.tust.2005.06.004 [22] 谭青,张魁,周子龙,等.球齿滚刀作用下岩石裂纹的数值模拟与试验观测[J].岩石力学与工程学报,2010,29(1):163–169. Tan Qing,Zhang Kui,Zhou Zilong,et al.Numerical simulation and experimental observation of rock cracks under action of spherical tooth hob cutter[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(1):163–169 [23] 汪珂.临空面辅助滚刀破岩的数值模拟研究[J].铁道工程学报,2020,37(10):96–102. doi: 10.3969/j.issn.1006-2106.2020.10.018 Wang Ke.Numerical simulation research on the rock fracturing by a TBM disc cutter assisted by free faces[J].Journal of Railway Engineering Society,2020,37(10):96–102 doi: 10.3969/j.issn.1006-2106.2020.10.018 [24] Liu Xiaohui,Liu Songyong,Ji Huifu.Mechanism of rock breaking by pick assisted with water jet of different modes[J].Journal of Mechanical Science and Technology,2015,29(12):5359–5368. doi: 10.1007/s12206-015-1137-3 [25] 谭青,易念恩,夏毅敏,等.TBM滚刀破岩动态特性与最优刀间距研究[J].岩石力学与工程学报,2012,31(12):2453–2464. doi: 10.3969/j.issn.1000-6915.2012.12.009 Tan Qing,Yi Nianen,Xia Yimin,et al.Research on rock dynamic fragmentation characteristics by TBM cutters and cutter spacing optimization[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(12):2453–2464 doi: 10.3969/j.issn.1000-6915.2012.12.009
