Processing math: 17%
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

基于空腔膨胀理论的TBM盘形滚刀侵岩载荷模型建模分析及试验验证

张魁, 张玉林, 郑学军, 周友行, 乔硕

张魁, 张玉林, 郑学军, 周友行, 乔硕. 基于空腔膨胀理论的TBM盘形滚刀侵岩载荷模型建模分析及试验验证[J]. 岩土工程学报, 2022, 44(12): 2263-2271. DOI: 10.11779/CJGE202212013
引用本文: 张魁, 张玉林, 郑学军, 周友行, 乔硕. 基于空腔膨胀理论的TBM盘形滚刀侵岩载荷模型建模分析及试验验证[J]. 岩土工程学报, 2022, 44(12): 2263-2271. DOI: 10.11779/CJGE202212013
ZHANG Kui, ZHANG Yu-lin, ZHENG Xue-jun, ZHOU You-hang, QIAO Shuo. Establishment and experimental verification of rock penetration load model for TBM disc cutters based on cavity expansion theory[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(12): 2263-2271. DOI: 10.11779/CJGE202212013
Citation: ZHANG Kui, ZHANG Yu-lin, ZHENG Xue-jun, ZHOU You-hang, QIAO Shuo. Establishment and experimental verification of rock penetration load model for TBM disc cutters based on cavity expansion theory[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(12): 2263-2271. DOI: 10.11779/CJGE202212013

基于空腔膨胀理论的TBM盘形滚刀侵岩载荷模型建模分析及试验验证  English Version

基金项目: 

国家自然科学基金项目 51704256

国家自然科学基金项目 11832016

国家自然科学基金项目 52175254

湖南省科技创新计划项目 2021RC2094

湖南省科技创新计划项目 2020RC2037

国家重点研发计划项目 YFB400800

合肥通用机械研究院有限公司项目 2021ZKKF043

详细信息
    作者简介:

    张魁(1985—),男,博士,副教授,硕士生导师,主要从事大型全断面隧道掘进装备破岩机理及其制造工艺研究。E-mail:zhangk@xtu.edu.cn

    通讯作者:

    郑学军, E-mail: zhengxuejun@xtu.edu.cn

  • 中图分类号: TU94

Establishment and experimental verification of rock penetration load model for TBM disc cutters based on cavity expansion theory

  • 摘要: 刃底密实核动态衍生现象严重影响了盘形滚刀的破岩机理,进而一定程度上影响了全断面岩石隧道掘进机(TBM)的掘进效能。参考H Alehossein钝楔形压头侵岩模型,将密实核衍生现象及密实核区周围复杂的应力场纳入考察,基于空腔膨胀理论和离散化建模思路,建立了常截面平刃滚刀侵岩模型。随后,分别以滚刀所受垂直力和刃底无因次平均接触应力为性能指标,将滚刀侵岩过程中刀圈刃部横截面几何尺寸参数和岩石物理力学性能参数作为考察变量,针对上述所建侵岩模型开展了十因素两水平的正交试验分析,探究了各考察变量的敏感程度,并分析了无因次平均接触应力的变化规律。最后,选取了2种具有不同粒径的岩石试样,制备了4种刀圈压头,并在TRW-3000三轴冲击破岩试验机上开展了侵岩试验验证。研究结果表明:对于垂直力和无因次平均接触应力而言,刀刃宽和侵深均为敏感因素;由于将密实核现象纳入考核,在滚刀压头侵岩作用下位于刃底的岩石应力水平明显大于单轴抗压强度;所建模型对于在第一次阶跃破碎之前的侵岩载荷速率的预测误差最大不超过10%,这表明模型具有良好的准确性。
    Abstract: The dynamic derivation of dense core at the cutter bottom seriously affects the rock-breaking mechanism of disc cutter (hereinafter referred to as cutter), and then affects the tunneling efficiency of full face rock tunnel boring machine (TBM) to a certain extent. Referring to the H. Alehossein blunt wedge indenter rock penetration model, the derivation phenomenon of dense core and the complex stress field around the dense core area is investigted, and a constant cross-section flat edge disc cutter (hereinafter referred to as cutter) rock penetration model is established based on the cavity expansion theory and discrete modeling idea. Then, taking the vertical force and dimensionless average contact stress as the performance evaluation indices, and selecting the geometric dimension parameters of the blade section of the cutter ring and the physical and mechanical property parameters of the rock as the variables, the ten-factor and two-level L12 (211) orthogonal tests are carried out. It is found that the blade width and penetration depth are the most sensitive factors. After that, the change characteristics of dimensionless contact stress at the cutter ring bottom are analyzed. It is found that because the dense core phenomenon is considered in the model, the stress of the rock under the cutter ring bottom is significantly greater than the uniaxial compressive strength, which is more in line with the engineering practice. Finally, two kinds of rock samples with different particle sizes are selected, four kinds of indenters are prepared, and the verification experiment of the theoretical model is carried out using the TRW-3000 rock breaking tester. The results show that both the cutter width and the penetration depth are sensitive factors for vertical force and dimensionless average contact stress. Because the dense core phenomenon is included in the assessment, the rock stress level at the cutter bottom is obviously greater than the uniaxial compressive strength under the intrusion of indenter. The maximum prediction error of the proposed model is less than 10%, which shows that it has good accuracy.
  • 压头侵岩机理及其理论建模一直是岩石力学领域的经典研究课题,受到众多学者的广泛关注。基于弹性力学理论,很多学者相继推导了压头侵岩过程中岩石的应力分布公式[1]。上述理论研究成果为进一步建立全断面岩石隧道掘进机(TBM)盘形滚刀(以下简称滚刀)破岩载荷模型奠定了基础。然而,由于上述模型对复杂刀岩作用过程进行了较大的假设简化,故其仅在分析压头静压破碎时具有较高的参考价值,尚不能较为准确全面地将岩石宏微观破碎特征纳入考察。

    在前述压头侵岩理论研究的基础上,许多学者基于滚刀线切割、滚压破岩实验数据,开展了滚刀载荷预测模型的搭建工作。IEvans等[2]和Roxborough等[3]认为,滚刀破岩时所受垂直力与滚刀压入岩石区域后在岩石表面上的投影面积成正比,并假定其比值为岩石的单轴抗压强度,进而基于该推论建立了滚刀垂直力预测公式。Rostami[4]根据线切割试验和压头压痕试验,建立了滚刀受力预测模型(业界广泛引述为CSM模型)。Rostami[5]随后进一步完善了滚刀线性切割试验方案,通过在滚刀刀圈侧面等周向间隔地安装多个应变片,以期对滚刀滚压破岩过程中的刀岩接触力进行实时监测,其研究结果表明:刀岩接触力存在前后卸荷破碎区,且应力集中程度较高。除了上述国外学者,国内多家科研院所也对滚刀侵岩载荷开展了理论建模研究。例如,上海交通大学640教研室[6]、东北工学院岩石破碎研究室[7]、中南大学[8]也分别提出了相应的受力预测模型,均具有一定的应用价值。也有很多学者基于现有模型(如CSM模型),对滚刀侵岩过程开展了诸多研究,虽然对侵岩理论的发展作出了贡献,但并未改变侵岩模型的机理。总的来说,由于学术界对于刀岩作用过程认识的局限性,加之岩石本构行为的复杂性,导致上述模型多基于试验数据回归得出,或基于大量简化假设推导得出,忽略了刀岩作用过程中客观存在的一些重要特征。

    正如文献研究所表明的,密实核现象就是滚刀破岩过程中必然发生的、且不可忽视的重要现象之一;侵岩试验中绝大部分(70%~85%)能量主要用于形成密实核区域[9]。也正是因为密实核的存在,使得在实际破岩过程中刀岩接触力水平要远大于现有理论预估值,进而导致刀具选型与更换缺乏理论指导,最终引起滚刀磨损失效等故障事件频发。据不完全统计,刀具费用约占掘进施工总费用的1/3[10],TBM停机后有2/3的时间用于刀具更换[11]。可见,将密实核现象纳入考察,并进一步完善发展滚刀侵岩载荷预测模型,具有一定的理论研究价值与工程应用前景。

    在现有科学技术水平的制约下,尽管密实核现象难以被直接加以在线观测与定量研究,但仍然有部分文献对其做出了前瞻性研究。Sanio等[12]在其空腔膨胀模型中指出,楔形滚刀下方存在应力状态类似于静水压力状态的致密岩粉团(即为密实核),并在此基础上给出了密实核区压力计算模型。Alehossein等[13]基于空腔膨胀模型,推导出了钝楔形压头垂直侵岩作用下,密实核区域周围的位移场、应力场和速度场的计算公式。张晓平等[14]对脆性岩石的破坏进行了大量试验研究,指出脆性材料的破坏分为:裂隙压密阶段、弹性阶段、裂纹稳定扩展阶段、裂纹加速扩展阶段和峰后阶段5个阶段。Liu等[15-16]开展了常截面平刃压头侵入板状花岗岩试验,观测到了致密的密实核。Liu等[17]采用红外热成像和声发射等无损检测技术,对密实核现象进行了观测,揭示了岩石损伤区的发展规律。李克金等[18]基于空腔膨胀理论,分析了各因素对密实核等区域边界尺寸的影响规律。此外,也有一些学者在建立常截面滚刀侵岩受力模型时尝试将密实核现象纳入考察。例如,孙伟等[19]和Huo等[20]将滚刀破岩过程分为弹性变形阶段、挤压破碎阶段、密实核破碎阶段与卸载阶段4个阶段,并采用分阶段建模思路,建立了滚刀受力预测模型。如何进一步考虑密实核区域周围的复杂场量及刀岩的三维耦合作用,建立出更加精确的侵岩载荷预测模型,依然是目前的研究难点。

    为了进一步完善滚刀侵岩载荷模型,特别是将密实核周围的复杂应力场纳入考察,本文基于空腔膨胀理论和离散化建模思路,建立了常截面平刃滚刀侵岩载荷预测模型;随后,以滚刀侵岩垂直力和无因次平均接触应力作为性能评价指标,选取刀圈刃形截面几何尺寸参数和岩石物理力学性能参数作为考察变量,开展了十因素两水平L12(211)正交试验分析,并根据正交试验所确定的敏感因素,分析了刃底无因次接触应力的变化特征。最后,为了验证所建理论模型的准确性,开展了滚刀侵岩试验。

    尽管常截面滚刀与H Alehossein侵岩研究中所述钝楔形压头在刃形尺寸上存在一定差异,但是二者所引入的岩石宏观破碎现象较为相似,故同样假定滚刀侵岩时存在:弹性变形阶段、弹塑性变形阶段、密实核生成与发育阶段、岩石裂纹衍生与损伤积累阶段和反复阶跃破碎阶段。刀具侵岩时,在岩石裂纹形成之前,同样认为岩石内部存在密实核、弹塑性失效区和弹性变形区3个应力区。

    不同于钝楔形压头,滚刀具有相对较为复杂的三维几何形状,导致滚刀侵岩理论建模困难。为此,尝试对滚刀进行如下离散化处理。如图 1所示,将刀岩接触段PRT上的刀圈按等角度间隔Δβ剖分成N份滚刀微元,并按照从左到右的顺序依次编号。N越大,则每份滚刀微元厚度Δβ就越小。当N趋近于无穷大时,每份滚刀微元厚度Δβ近似为零。定义OP为剖切角起始位置,逆时针为正,则第i份滚刀微元对应的剖切角βi

    βi=iΔβ (1)
    图  1  滚刀离散化及滚刀微元的定义
    Figure  1.  Discretization of cutter and definition of cutter element

    i=N时,称此时的剖切角为刀岩理论接触角ψ,由几何关系可得

    ψ=arccos(1h/r) (2)

    式中,h为滚刀侵深,r为滚刀半径。

    滚刀微元在其剖面上的运动可以视为垂直侵岩运动,由几何关系可得其等效侵深hβ

    hβ=r(rh)/cos(ψβi) (3)

    根据常截面滚刀的几何特征,如图 2所示,可以得到滚刀微元刀岩理论接触半径aβ的计算公式为

    aβ=2r0hβh2β+a0 (4)
    图  2  给定剖切角βi时滚刀微元侵岩理论模型示意图
    Figure  2.  Schematic diagram of theoretical model for rock penetration of cutter element under given cutting angleβi

    式中,r0为刃底圆弧段半径,2a0为刀刃宽度。

    为了便于建模,忽略裂纹衍生与损伤积累阶段中密实核物质的流失,仅考察首次阶跃破碎前的侵岩过程。在前述滚刀侵岩基本假设基础上,进一步假定在剖切角为βi的滚刀微元作用下,且在裂纹形成之前,其刃底岩体同样存在如图 2所示的3个应力区。图 2中,以岩面与滚刀中间对称面的交点O''为极坐标原点,半径坐标为R;密实核区边界的半径坐标为aβ,弹塑性失效区域外边界半径坐标为r。假设半圆形密实核区域处于静水压力状态,各点应力大小均为P

    为了建模方便,对滚刀微元的各项尺寸参数进行无因次化处理:定义无因次径向坐标为ξ=R/aβ,定义无因次侵深为δβ=hβ/2a0。则密实核边界的无因次半径为α=aβ/2a0,弹塑性失效区域外边界的无因次半径为ξ=r/aβ。如图 2所示,ξ为密实核区域, 1<\xi \le {\xi }_{*} 为弹塑性失效区域, \xi >{\xi }_{*} 为弹性变形区域。

    因任一滚刀微元的厚度远小于其刀刃宽,故可简化为平面应力问题。根据H Alehossein给出的空腔膨胀模型,滚刀下方存在的3个不同区域有着不同的本构模型:

    (1) 静水压力状态的半径为{a_\beta }的密实核区域。

    (2) 半径在{a_\beta }{r_*}之间的弹塑性失效区,应力满足Mohr-Coulomb屈服条件,在极坐标下可写为

    {K_\text{p}}({\sigma _\theta } - \zeta ) - ({\sigma _\text{r}} - \zeta ) = 0\text{,} (5)

    式中,\zeta = {\sigma _\text{c}}/({K_\text{p}} - 1){K_\text{p}} = (1 + \sin {\varphi _\text{b}})/(1 - \sin {\varphi _\text{b}}){\sigma _\text{c}}为岩石抗压强度,{\varphi _\text{b}}为岩石内摩擦角。

    (3) 半径在{r_*}之外的弹性区域,应力由Lamé解给出。在极坐标系下岩体弹性本构方程为

    \left. \begin{array}{l}{\epsilon }_{r}=({\sigma }_{r}-\nu {\sigma }_{\theta })/{E}_{r}\text{ }\text{,}\\ {\epsilon }_{\theta }=({\sigma }_{\theta }-\nu {\sigma }_{r})/{E}_{r}\text{ }\text{,} \end{array} \right\} (6)

    式中,\nu 为岩石泊松比,{\varepsilon _\text{r}}{\varepsilon _\theta }分别为岩体弹性变形区内某点的径向应变和环向应变,{\sigma _\text{r}}{\sigma _\theta }分别为径向应力和环向应力。

    基于H Alehossein钝楔形压头侵岩模型,进一步推导出了滚刀微元侵岩模型中应力场、位移场和速度场的基本计算公式。

    (1) 应力场

    密实核区应力场计算式为

    P=2{K}_\text{p}{S}_\text{o}^\text{l}{\xi }_{*}^{({K}_\text{p}-1)/{K}_\text{p}}/({K}_\text{p}-1)-\zeta \text{,} (7)

    式中, {S}_\text{o}^\text{l}={\sigma }_\text{c}/({K}_\text{p}+1)

    弹塑性失效区应力场计算公式为

    \left.\begin{array}{l} \sigma_r=\zeta-\left[\left(2 K_{\mathrm{p}} S_{\mathrm{o}}^1\right) /\left(K_{\mathrm{p}}-1\right)\right]\left(\xi_* / \xi\right)^{\left(K_{\mathrm{p}}-1\right) / K_{\mathrm{p}}}, \\ \sigma_\theta=\zeta-\left[\left(2 S_{\mathrm{o}}^1\right) /\left(K_{\mathrm{p}}-1\right)\right]\left(\xi_* / \xi\right)^{\left(K_{\mathrm{p}}-1\right) / K_{\mathrm{p}}} 。 \end{array}\right\} (8)

    \xi = {\xi _*},代入式(8),即可得到弹塑性失效区外边界上的应力状态。

    (2) 位移场

    在弹性变形区域,由拉美公式可得径向位移:

    {u_\text{r}} = \kappa a_\beta ^2\xi _*^2/2r = {a_\beta }\kappa \xi _*^2{\xi ^{ - 1}}/2\text{,} (9)

    式中,\kappa = {\sigma _\text{c}}(1 + \nu )(1 - {\text{sin}}{\varphi _\text{b}})/{E_\text{r}}

    (3) 速度场

    给定{\xi _*}{a_\beta }时,定义弹塑性失效区内任一点的位置变化率为{v_\text{r}},则{v_\text{r}}关于\xi 的一阶微分方程为

    {\text{d}}{v_\text{r}}/({\text{d}}\xi ) + K_\text{d}^{ - 1}{v_\text{r}}/\xi + K_\text{d}^{ - 1}A{\xi ^{ - \left( {{K_\text{p}} - 1} \right)/{K_\text{p}}}} = 0\text{,} (10)

    式中,{K_d} = (1 + \sin {\psi _0})/(1 - \sin {\psi _0}){\psi _0}为岩石剪胀角, A=\lambda \kappa \left[1+({a}_{\beta }\text{d}{\xi }_{\text{*}})/({\xi }_{\text{*}}\text{d}{a}_{\beta })\right]{\xi }_{\text{*}}{}^{\left({K}_\text{p}-1\right)/{K}_\text{p}},\lambda =\text{(}1-\nu {K}_\text{d}+ {K_\text{p}}{K_\text{d}} - \lambda = {K_\text{p}}\nu {\text{)}}/\left[ {(1 + \nu ){K_\text{p}}} \right]。

    解方程后可得速度场的表达式为

    {v_\text{r}}(\xi ) = {A_0}{\xi ^{ - 1/{K_\text{d}}}} - \mu \kappa {\xi _*} + {a_\beta }\frac{{{\text{d}}{\xi _*}}}{{{\text{d}}{a_\beta }}}{\left( {\frac{{{\xi _*}}}{\xi }} \right)^{\left( {{K_\text{p}} - 1} \right)/{K_\text{p}} - 1}}。 (11)

    如果获得弹塑性失效区域外边界的无因次半径{\xi _*}关于无因次侵深{\delta _\beta }的关系式,就可以利用该式计算出各个区域内任一点的应力、位移和速度值,以及密实核内应力值P。然后,根据刀岩接触区面积大小,即可求取滚刀微元的侵岩载荷。

    由于滚刀微元具有与所述钝楔形压头截然不同的刃形特征,因此需要重新修正速度场的边值条件。当刀岩接触半径增加{\text{d}}{a_\beta }时,在密实核边界上(\xi = 1),密实核膨胀体积与刀岩挤占体积的关系如图 3所示,边界上的速度满足如下关系式:

    {v_\text{r}}(1) = {\text{d}}{S_\text{c}}/{\text{d}}{S_\text{r}} = 2\tan \eta /\pi \text{,} (12)
    图  3  侵岩过程中刀岩体积协调关系示意图
    Figure  3.  Schematic diagram of cutter-rock volume coordination during rock penetration

    式中,{\text{d}}{S_r}为密实核膨胀面积{\text{d}}{S_r} = \pi {a_\beta }{\text{d}}{a_\beta }{\text{d}}{S_c}为刀刃挤压面积{\text{d}}{S_c} = 2{a_\beta }\tan \eta {\text{d}}{a_\beta }\eta 满足关系式:

    \tan \eta = ({a_\beta } - {a_0})/\sqrt {r_0^2 - {{({a_\beta } - {a_0})}^2}} 。 (13)

    \xi = {\xi _*},对式(9)求解拉格朗日时间导数可得

    {v_r}({\xi _*}) = \kappa \left( {{\xi _*} + {a_\beta }\frac{{{\text{d}}{\xi _*}}}{{{\text{d}}{a_\beta }}}} \right)。 (14)

    式(12),(14)为式(11)的边值条件。

    \xi = {\xi _*}代入式(11)可得

    {v_r}\left( {{\xi _*}} \right) = {A_0}\xi _*^{ - 1/{K_d}} - \mu \kappa \left( {{\xi _*} + {a_\beta }\frac{{{\text{d}}{\xi _*}}}{{{\text{d}}{a_\beta }}}} \right)。 (15)

    联立式(12),(14),(15)可得{\xi _*}关于{a_\beta }的一阶常微分方程:

    \frac{{{\text{d}}{\xi _{\text{*}}}}}{{{\text{d}}{a_\beta }}} = \frac{1}{{{a_\beta }}}\left[ { - {\xi _{\text{*}}} + \frac{{\gamma ({a_\beta })}}{{(1 + \mu ){\xi _{\text{*}}}^{1/{K_d}} - \mu {\xi _{\text{*}}}^{\left( {{K_\text{p}} - 1} \right)/{K_\text{p}} - 1}}}} \right]\text{,} (16)

    式中,

    \gamma ({a_\beta }) = 2\tan \eta /(\pi \kappa )。 (17)

    根据{\delta _\beta }的定义,可将式(4)改写为

    {a_\beta } = {a_0} + 2\sqrt {{a_0}{r_0}{\delta _\beta } - a_0^2\delta _\beta ^2} 。 (18)

    联立式(13),(18),可将式(18)改写为

    \gamma ({\delta _\beta }) = 4\sqrt {{a_0}{r_0}{\delta _\beta } - a_0^2\delta _\beta ^2} /\left[ {\pi \kappa ({r_0} - 2{a_0}{\delta _\beta })} \right]。 (19)

    对式(18)两边取微分得

    {\text{d}}{a_\beta } = \left[ {({a_0}{r_0} - 2a_0^2{\delta _\beta })/(\sqrt {{a_0}{r_0}{\delta _\beta } - a_0^2\delta _\beta ^2} )} \right]{\text{d}}{\delta _\beta }。 (20)

    将式(18)~(20)代入式(16),即可得到弹塑性失效区外边界无因次半径{\xi _*}关于无因次侵深{\delta _\beta }的方程:

    \frac{{{\text{d}}{\xi _*}}}{{{\text{d}}{\delta _\beta }}} = \frac{{{r_0} - 2{a_0}{\delta _\beta }}}{{\sqrt {{a_0}{r_0}{\delta _\beta } - a_0^2\delta _\beta ^2} + 2{\text{(}}{r_0}{\delta _\beta } - {a_0}\delta _\beta ^2{\text{)}}}} \cdot \\ \left[ { - {\xi _*} + \frac{{\gamma ({\delta _\beta })}}{{(1 + \mu ){\xi _*}^{1/{K_d}} - \mu {\xi _*}^{\left( {{K_\text{p}} - 1} \right)/{K_\text{p}} - 1}}}} \right]。 (21)

    \gamma ({\delta _\beta }) = 1即可得到所述方程的初始条件。

    给定滚刀侵深h以及滚刀微元任一剖面位置{\beta _i},参考3.1节基本场量公式,可获得岩体内任一位置处的应力、位移和速度大小。

    图 4所示,从图 1中任取一剖切角为{\beta _i}的滚刀微元,{F_i}为该滚刀微元所受合力,类比滚刀在滚压过程中的受力,可进一步将{F_i}分解为垂直于掌子面方向的垂直力\Delta {F_\text{v}}、平行于滚刀运动方向的切向力\Delta {F_\text{t}}和垂直于滚刀运动方向的侧向力\Delta {F_\text{s}}

    \left. \begin{array}{l}\Delta {F}_{\text{v}}=2P({\beta }_{i})a({\beta }_{i})r\mathrm{cos}(\psi -{\beta }_{i})\Delta \beta \text{ }\text{,}\\ \Delta {F}_{\text{t}}=2P({\beta }_{i})a({\beta }_{i})r\mathrm{sin}(\psi -{\beta }_{i})\Delta \beta \text{ }\text{,}\\ \Delta {F}_{text{s}}=P({\beta }_{i})r{h}_{\beta }\Delta \beta \text{ }。 \end{array} \right\} (22)
    图  4  滚刀及其微元受力分析示意图
    Figure  4.  Schematic diagram of force analysis of cutter and cutter element

    根据叠加原理,对滚刀微元各切削力分量进行累加求和,即可得到滚刀侵岩时总的切削力分量:

    \left. \begin{array}{l}{F}_\text{v}/2= \sum \limits_{i=1}^{N}\Delta {F}_\text{v}\text{,}\\ {F}_\text{t}=\sum \limits_{i=1}^{N}\Delta {F}_\text{t}\text{,}\\ {F}_\text{s}=\sum \limits_{i=1}^{N}\Delta {F}_\text{s}。 \end{array} \right\} (23)

    为了探究考虑密实核现象后刀岩接触应力水平的变化差异,进一步推导了刃底无因次平均接触应力{\sigma '_\text{m}}

    {\sigma '_\text{m}} = \frac{{\mathop \sum \nolimits_{i = 1}^N \Delta {F_\text{v}}}}{{\mathop \sum \nolimits_{i = 1}^N 2ra(\beta )\cos (\psi - \beta )\Delta \beta }}/{\sigma _\text{c}}。 (24)

    为了后续有针对性地开展滚刀侵岩过程的析因研究,现分别以侵岩垂直力{F_\text{v}}和无因次平均接触应力{\sigma '_\text{m}}为性能指标,将滚刀侵岩过程中刀圈刃部横截面几何尺寸参数(刀刃宽、刀刃角、过渡圆弧半径、刀刃半径)、岩石物理力学性能参数(抗压强度、抗拉强度、剪切强度、内摩擦角、剪胀角)以及侵深作为考察变量,各设计了一组十因素两水平的正交试验,标记为L12(211),其正交试验方案见表 1。如表 23所示,分别为以侵岩垂直力{F_\text{v}}和无因次平均接触应力{\sigma '_\text{m}}为性能指标,获得的极差分析表。

    表  1  L12(211)正交试验方案表
    Table  1.  List of L12(211) orthogonal level
    水平 侵深A/mm 刀刃宽B/mm 刀刃角C/(°) 过渡圆弧半径D/mm 滚刀半径E/mm 内摩擦角F/(°) 剪胀角G/(°) 抗压强度H/MPa 抗拉强度I/MPa 剪切强度J/MPa
    1 1.0 4.5 5 2.3 108 45 15 97 7.2 15
    2 1.6 8.5 10 6.9 216 55 25 137 8.2 25
    下载: 导出CSV 
    | 显示表格
    表  2  极差分析表(垂直力{F_\text{v}})
    Table  2.  Range analysis table (normal force{F_\text{v}})
    项目 A B C D E F G H I J 空列
    K1 496.59 609.07 676.11 657.22 547.92 697.65 602.71 631.23 650.54 661.67 673.91
    K2 824.04 711.56 644.51 663.40 772.70 622.97 717.91 870.77 670.09 658.96 646.71
    k1 82.76 101.51 112.69 109.54 91.32 116.27 100.45 105.21 108.42 110.28 112.32
    k2 137.34 118.59 107.42 110.57 128.78 103.83 119.65 145.13 111.68 109.83 107.79
    Rj 327.45 102.49 31.60 6.18 224.78 74.68 115.20 58.16 19.55 2.71 27.20
    注:AEGBFHCIDJ
    下载: 导出CSV 
    | 显示表格
    表  3  极差分析表(无因次接触应力{\sigma '_\text{m}})
    Table  3.  Range analysis table (dimensionless contact stress{\sigma '_\text{m}})
    项目 A B C D E F G H I J 空列
    K1 7.82 10.36 8.89 9.59 8.88 9.23 8.06 9.96 9.32 9.17 8.97
    K2 10.27 7.74 9.2 8.50 9.22 8.87 10.04 9.99 8.78 8.93 9.13
    k1 1.30 1.73 1.48 1.60 1.48 1.54 1.34 1.66 1.55 1.53 1.49
    k2 1.71 1.29 1.54 1.42 1.54 1.48 1.67 1.66 1.46 1.49 1.52
    Rj 2.45 2.62 0.33 1.09 0.34 0.35 1.98 1.83 0.53 0.24 0.16
    注:BAGHDIFECJ
    下载: 导出CSV 
    | 显示表格

    (1) 对于侵岩垂直力{F_\text{v}}而言,各因素的显著性顺序为侵深>刀刃半径>剪胀角>刀刃宽>内摩擦角>抗压强度>刀刃角>抗拉强度>过渡圆弧半径>剪切强度,其中,侵深、刀刃半径、剪胀角和刀刃宽的影响最为显著,内摩擦角和抗压强度次之,其他因素均为不显著。上述分析结论与孙剑萍等[21]基于四因素正交试验所获结果相类似,即侵深和刀刃半径对滚刀的破岩受力影响最为显著,进而对滚刀磨损的影响也最为显著。林赉贶等[22]开展了安装参数和掘进参数对滚刀破岩阻力影响的正交试验,同样指出侵深为影响滚刀切削阻力的最显著因素之一。

    (2) 对于无因次平均接触应力{\sigma '_\text{m}}而言,各因素的显著性顺序为刀刃宽>侵深>剪胀角>抗压强度>过渡圆弧半径>抗拉强度>内摩擦角>刀刃半径>刀刃角>剪切强度,其中,刀刃宽和侵深的影响较为显著,其他因素均为不显著因素。

    为了进一步确定刃底无因次平均接触应力的变化特征,给定花岗岩物理力学参数:{\sigma _\text{c}}=158 MPa,{E_\text{r}} =41.0 GPa,\nu =0.234,{\varphi _\text{b}}=60°,{\psi _0}=30°。给定平刃滚刀的尺寸参数:{a_0}=6.5 mm,{r_0}=2.3 mm。初始条件:{\xi _*}{\text{ = }}1a=6.5073 mm,{\delta _{\max }}=0.09。给定刀刃宽2{a_0}为9,13,17,21 mm,利用式(24)计算得到刃底无因次平均接触应力{\sigma '_\text{m}}随着侵深h的变化曲线图,如图 5所示。

    图  5  无因次平均接触应力随侵深的变化曲线
    Figure  5.  Variation of dimensionless average contact stresses with penetration depths

    图 5可知,刃底无因次平均接触应力{\sigma '_\text{m}}并非一个固定值,而是随着侵深h的增加而增加;当侵深h相同时,平均接触应力{\sigma '_\text{m}}随着刀刃宽2{a_0}的增加而减少;当侵深h较大时,由于侵岩过程中位于刃底的密实核受周围完整原岩的限制作用,平均接触应力{\sigma '_\text{m}}明显大于岩石抗压强度{\sigma _\text{c}}(图 5中表现为当侵深h增大时,无因次平均接触应力{\sigma '_\text{m}} > 1)。这明显与Evans等[2]给出的“刃底平均接触应力近似为{\sigma _\text{c}}”的假设不同。由于侵岩模型将密实核现象纳入了考察,因此更符合客观实际。

    为了验证前述侵岩理论模型的准确性,在图 6所示TRW-3000三轴破岩试验平台上开展了近似无围压条件下滚刀压头侵岩试验。该试验平台利用位于水平方向的XY向油缸提供侧向围压以夹持岩石试样(侧向围压设定为0.5 MPa,故可近似视作无围压条件),同时利用位于垂直方向的Z向油缸提供侵岩压力。

    图  6  TRW-3000三轴冲击破岩试验平台
    Figure  6.  TRW-3000 rock-breaking test machine

    试验采用4种编号为C1,C2,C3,C4的常截面滚刀压头,其中C1,C4为自制压头,其余压头均由工程用刀圈经电火花加工而成。各压头的截面尺寸参数如表 4所示。试验采用编号为R1,R2的两种不同粒径的中国花岗岩,即汨罗麻岩(细花)和金点隋麻(粗花)。如图 7所示,为制备的用于性能标定测试的标准圆柱形岩样,其在0.5 MPa围压下所测得的物理力学性能参数如表 5所示。为了降低岩样尺寸边界效应的不利影响,试验采用几何尺度远大于滚刀压头的块状岩样。

    表  4  滚刀压头截面尺寸参数表
    Table  4.  Cross-sectional size parameters of indenters
    编号 2{a_0}/mm 2{\theta _0}/(°) {r_0}/mm r/mm
    C1 6.5 20 1.75 108
    C2 11.3 13 2.80 216
    C3 11.7 18 5.70 216
    C4 13.0 20 3.50 216
    下载: 导出CSV 
    | 显示表格
    图  7  试验所用柱形花岗岩截面图
    Figure  7.  Sections of columnar granite used in tests
    表  5  0.5 MPa围压下岩石试样物理力学性能参数表
    Table  5.  Mechanical property parameters of rock samples under 0.5 MPa confining pressure
    编号 弹性模量/GPa 泊松比 抗压强度/MPa 抗拉强度/MPa 黏聚力/MPa 内摩擦角/(°) 剪胀角/(°)
    R1 41.94 0.18 136.81 8.19 22.25 52.73 20.08
    R2 41.74 0.15 97.43 7.25 16.59 54.12 21.75
    下载: 导出CSV 
    | 显示表格

    (1) 预压固定:XY方向按100 N/s的加载速度同时施加向紧固围压。随后,Z向加载单元按100 N/s的加载速度对压头施加10 kN垂直预压。

    (2) 正式加载及数据记录:保持侧向紧固围压恒定,静置10 min后,将Z向位移/载荷读数清零以消除载荷曲线调整期的影响。Z向加载单元以1 mm/min的加载速度侵入试样至给定侵深或产生第一次阶跃破碎,期间记录刀具垂直侵岩载荷曲线。

    (3) 卸载观察:一次试验完成后,按照先Z向后XY向的顺序卸载油缸压力,卸载速率不大于10 N/s以防止发生岩爆。待完全卸载后,取出岩样对压痕进行拍照观察,并记录压痕形状尺寸。

    (4) 重复试验:分别开展4种压头侵入R1和R2共计8组侵岩试验。其中,每组至少进行3次重复性试验,以获取相似度较高的3条载荷曲线样本。

    按照格里菲斯理论[23]可知,岩石试样等天然材料内部存在诸多微裂隙、气孔等缺陷,因此在外力作用下,所述缺陷周围极易产生应力集中现象[24]。在侵岩试验中,同样也观测到了上述现象,例如部分试样在加载末期沿压头中间对称面发生了脆断,如图 8所示。

    图  8  岩样脆断面的密实核形态图
    Figure  8.  Dense core morphology of brittle section of rock sample

    可以看出,这些脆断面上均存在充满致密岩粉的弧形密实核区,其特征为中间区域厚边缘薄,且左右基本对称。该形态特征与Yin等[25]在TBM滚刀侵岩试验中所观测到的密实核形态基本一致。

    针对每组试验所获侵岩载荷曲线,对其在第一次阶跃破碎前的曲线段进行拟合,获得实测拟合曲线及其拟合斜率,并求取拟合斜率的平均值。类似地,对相似条件下的理论侵岩载荷曲线,进行拟合,获得理论拟合曲线及其拟合斜率。为了便于对比分析,针对每一种岩石试样,将不同压头的实测/理论载荷曲线及其拟合曲线绘制于同一张图中,如图 9所示。

    图  9  0.5 MPa围压下各压头的实测、理论曲线
    Figure  9.  Measured and theoretical curves of each indenter under confining pressure of 0.5 MPa

    从整体上来看,左侧理论载荷拟合曲线基本上与右侧的实测载荷拟合曲线关于Y轴基本对称。这说明理论载荷曲线与实测曲线大致重合。如表 6所示,为不同条件下实测侵岩载荷曲线的拟合斜率平均值、理论拟合斜率及其误差值的汇总表;由表可知,最大预测误差不超过10%;这说明本论文所提出的滚刀侵岩理论模型在预测滚刀压头侵入特定花岗岩时垂直力变化速率方面具有一定的准确性。当然,由于所建理论模型是在上述特定简化假设的基础上推导完成,故所述模型仍存在如下局限性:

    表  6  侵岩载荷曲线斜率的误差分析
    Table  6.  Slope error analysis of cutting load curves
    组号 实测值
    /(kN·mm-1)
    理论值
    /(kN·mm-1)
    相对误差
    /%
    C1-R1 126.6 119 6
    C2-R1 212.6 200 6
    C3-R1 257.4 239 7
    C4-R1 231.3 219 5
    C1-R2 115.1 107 7
    C2-R2 215.3 201 7
    C3-R2 197.2 192 3
    C4-R2 214.9 203 6
    下载: 导出CSV 
    | 显示表格

    (1) 忽略了裂纹衍生与损伤积累阶段密实核物质的流失,仅考察了首次阶跃破碎前的侵岩过程,对于首次阶跃破碎之后的侵岩过程还需要进行深入研究。

    (2) 仅针对两种花岗岩试样进行了试验验证,因此对于其他种类岩石的适用性尚待进一步研究。

    此外,需要补充解释的是,部分试验曲线在侵岩初期出现了小幅波动(如图 9中虚线框所示),导致理论预测载荷曲线相对实测出现了一定的偏离。这可能是由于试验装置的局限性所导致的,而并非是理论建模预测精度的问题;未来可通过提高压头加工精度和增加紧固装置的方式来进一步提高试验过程稳定性。

    本文基于空腔膨胀理论,采用离散化建模思路,将滚刀侵岩过程中的密实核现象纳入考察,建立了可反映密实核区周围复杂应力关系的常截面平刃滚刀侵岩载荷预测模型,并以垂直力和无因次平均接触应力为指标,分别对所建模型开展了十因素两水平正交分析和试验验证,得到4点结论。

    (1) 正交试验结果表明,对于垂直力和无因次平均接触应力而言,其敏感因素大为不同,前者的敏感因素为侵深、刀刃半径、剪胀角、刀刃宽、内摩擦角和抗压强度;后者为刀刃宽和侵深。

    (2) 侵岩理论研究表明,刃底无因次平均接触应力并非一个固定值,而是随着侵深的增加而增加;最为重要的是,由于将密实核现象纳入考核,在滚刀压头侵岩作用下位于刃底的岩石应力水平明显大于单轴抗压强度。这与Evans I,Roxborough F F等给出的“刃底平均接触应力近似为{\sigma _\text{c}}”的假设明显不同,而更符合Rostami J提出的“刃底接触应力呈高度集中分布”的论断。

    (3) 在侵岩验证试验过程中,观察到了充满致密岩粉的弧形密实核区,其形态与理论假设以及其他学者得到的观察结果基本相符,但相关定量分析有待后续深入研究。

    (4) 侵岩验证试验表明,在给定条件下,所建模型对于在第一次阶跃破碎之前的侵岩载荷速率的预测误差最大不超过10%。这说明所建模型准确性良好。

  • 图  1   滚刀离散化及滚刀微元的定义

    Figure  1.   Discretization of cutter and definition of cutter element

    图  2   给定剖切角{\beta _i}时滚刀微元侵岩理论模型示意图

    Figure  2.   Schematic diagram of theoretical model for rock penetration of cutter element under given cutting angle{\beta _i}

    图  3   侵岩过程中刀岩体积协调关系示意图

    Figure  3.   Schematic diagram of cutter-rock volume coordination during rock penetration

    图  4   滚刀及其微元受力分析示意图

    Figure  4.   Schematic diagram of force analysis of cutter and cutter element

    图  5   无因次平均接触应力随侵深的变化曲线

    Figure  5.   Variation of dimensionless average contact stresses with penetration depths

    图  6   TRW-3000三轴冲击破岩试验平台

    Figure  6.   TRW-3000 rock-breaking test machine

    图  7   试验所用柱形花岗岩截面图

    Figure  7.   Sections of columnar granite used in tests

    图  8   岩样脆断面的密实核形态图

    Figure  8.   Dense core morphology of brittle section of rock sample

    图  9   0.5 MPa围压下各压头的实测、理论曲线

    Figure  9.   Measured and theoretical curves of each indenter under confining pressure of 0.5 MPa

    表  1   L12(211)正交试验方案表

    Table  1   List of L12(211) orthogonal level

    水平 侵深A/mm 刀刃宽B/mm 刀刃角C/(°) 过渡圆弧半径D/mm 滚刀半径E/mm 内摩擦角F/(°) 剪胀角G/(°) 抗压强度H/MPa 抗拉强度I/MPa 剪切强度J/MPa
    1 1.0 4.5 5 2.3 108 45 15 97 7.2 15
    2 1.6 8.5 10 6.9 216 55 25 137 8.2 25
    下载: 导出CSV

    表  2   极差分析表(垂直力{F_\text{v}})

    Table  2   Range analysis table (normal force{F_\text{v}})

    项目 A B C D E F G H I J 空列
    K1 496.59 609.07 676.11 657.22 547.92 697.65 602.71 631.23 650.54 661.67 673.91
    K2 824.04 711.56 644.51 663.40 772.70 622.97 717.91 870.77 670.09 658.96 646.71
    k1 82.76 101.51 112.69 109.54 91.32 116.27 100.45 105.21 108.42 110.28 112.32
    k2 137.34 118.59 107.42 110.57 128.78 103.83 119.65 145.13 111.68 109.83 107.79
    Rj 327.45 102.49 31.60 6.18 224.78 74.68 115.20 58.16 19.55 2.71 27.20
    注:AEGBFHCIDJ
    下载: 导出CSV

    表  3   极差分析表(无因次接触应力{\sigma '_\text{m}})

    Table  3   Range analysis table (dimensionless contact stress{\sigma '_\text{m}})

    项目 A B C D E F G H I J 空列
    K1 7.82 10.36 8.89 9.59 8.88 9.23 8.06 9.96 9.32 9.17 8.97
    K2 10.27 7.74 9.2 8.50 9.22 8.87 10.04 9.99 8.78 8.93 9.13
    k1 1.30 1.73 1.48 1.60 1.48 1.54 1.34 1.66 1.55 1.53 1.49
    k2 1.71 1.29 1.54 1.42 1.54 1.48 1.67 1.66 1.46 1.49 1.52
    Rj 2.45 2.62 0.33 1.09 0.34 0.35 1.98 1.83 0.53 0.24 0.16
    注:BAGHDIFECJ
    下载: 导出CSV

    表  4   滚刀压头截面尺寸参数表

    Table  4   Cross-sectional size parameters of indenters

    编号 2{a_0}/mm 2{\theta _0}/(°) {r_0}/mm r/mm
    C1 6.5 20 1.75 108
    C2 11.3 13 2.80 216
    C3 11.7 18 5.70 216
    C4 13.0 20 3.50 216
    下载: 导出CSV

    表  5   0.5 MPa围压下岩石试样物理力学性能参数表

    Table  5   Mechanical property parameters of rock samples under 0.5 MPa confining pressure

    编号 弹性模量/GPa 泊松比 抗压强度/MPa 抗拉强度/MPa 黏聚力/MPa 内摩擦角/(°) 剪胀角/(°)
    R1 41.94 0.18 136.81 8.19 22.25 52.73 20.08
    R2 41.74 0.15 97.43 7.25 16.59 54.12 21.75
    下载: 导出CSV

    表  6   侵岩载荷曲线斜率的误差分析

    Table  6   Slope error analysis of cutting load curves

    组号 实测值
    /(kN·mm-1)
    理论值
    /(kN·mm-1)
    相对误差
    /%
    C1-R1 126.6 119 6
    C2-R1 212.6 200 6
    C3-R1 257.4 239 7
    C4-R1 231.3 219 5
    C1-R2 115.1 107 7
    C2-R2 215.3 201 7
    C3-R2 197.2 192 3
    C4-R2 214.9 203 6
    下载: 导出CSV
  • [1] 徐小荷, 余静. 岩石破碎学[M]. 北京: 煤炭工业出版社, 1984.

    XU Xiao-he, YU Jing. Rock Fragmentation[M]. Beijing: China Coal Industry Publishing House, 1984. (in Chinese)

    [2]

    EVANS I, POMEROY C D. The Strength, Fracture and Workability of Coal[M]. London: Pergamon Press, 1966.

    [3]

    ROXBOROUGH F F, PHILLIPS H R. Rock excavation by disc cutter[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, , 1975, 12(12): 361–366.

    [4]

    ROSTAMI J. Design Optimization, Performance Predictions, and Economic Analysis of TBM Application in the Proposed Yucca Mountain Nuclear Waste Repository[D]. Golden: Colorado School of Mines, 1991.

    [5]

    ROSTAMI J. Development of a Force Estimation Model for Rock Fragmentation With Disc Cutters Through Theoretical Modeling and Physical Measurement of Crushed Zone Pressure [D]. Golden: Colorado School of Mines, 1997.

    [6] 孙鸿范, 陈健元, 陈刚. 掘进机盘形滚刀破岩力及计算载荷的研究[J]. 工程机械, 1980, 8: 1–7. https://www.cnki.com.cn/Article/CJFDTOTAL-GCJA198008000.htm

    SUN Hong-fan, CHEN Jian-yuan, CHEN Gang. Research on rock-breaking force and calculated load of disc cutter of cutterhead[J]. Construction Machinery and Equipment, 1980, 8: 1–7. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCJA198008000.htm

    [7] 余静. 岩石机械破碎规律和破岩机理模型[J]. 煤炭学报, 1982, 3: 10–18. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB198203001.htm

    YU Jing. Rules of rock fragmentation with mechanical methods and model of rock failure mechanism[J]. Journal of China Coal Society, 1982, 3: 10–18. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB198203001.htm

    [8]

    XIA Y M, TAO O Y, ZHANG X M, et al. Mechanical model of breaking rock and force characteristic of disc cutter[J]. Journal of Central South University, 2012, 19(7): 1846-1852. doi: 10.1007/s11771-012-1218-8

    [9] 刘泉声, 刘建平, 时凯, 等. 评价岩石脆性指标对滚刀破岩效率的影响[J]. 岩石力学与工程学报, 2016, 35(3): 498–510. doi: 10.13722/j.cnki.jrme.2015.0569

    LIU Quan-sheng, LIU Jian-ping, SHI Kai, et al. Evaluation of rock brittleness indexes on rock fragmentation efficiency by disc cutter[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(3): 498–510. (in Chinese) doi: 10.13722/j.cnki.jrme.2015.0569

    [10] 夏毅敏, 罗德志, 周喜温. 盾构地质适应性配刀规律研究[J]. 煤炭学报, 2011, 36(7): 1232-1236. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201107035.htm

    XIA Yi-min, LUO De-zhi, ZHOU Xi-wen. Study on the law of geology adaptability cutter selection for shield[J]. Journal of China Coal Society, 2011, 36(7): 1232-1236. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201107035.htm

    [11]

    ZHAO J, GONG Q M, EISENSTEN Z. Tunnelling through a frequently changing and mixed ground: a case history in Singapore[J]. Tunnelling and Underground Space Technology, 2007, 22(4): 388–400. doi: 10.1016/j.tust.2006.10.002

    [12]

    SANIO H P. Prediction of the performance of disc cutters in anisotropic rock[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1985, 22(3): 153–161.

    [13]

    ALEHOSSEIN H, DETOURNAY E, HUANG H. An analytical model for the indentation of rocks by blunt tools[J]. Rock Mechanics and Rock Engineering, 2000, 33(4): 267–284.

    [14] 张晓平, 王思敬, 韩庚友, 等. 岩石单轴压缩条件下裂纹扩展试验研究: 以片状岩石为例[J]. 岩石力学与工程学报, 2011, 30(9): 1772–1781. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201109008.htm

    ZHANG Xiao-ping, WANG Si-jing, HAN Geng-you, et al. Crack propagation study of rock based on uniaxial compressive test—a case study of schistose rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(9): 1772–1781. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201109008.htm

    [15]

    LIU J, CAO P. Study on rock fracture with TBM cutter under different confining stresses[J]. Indian Geotechnical Journal, 2016, 46(1): 104–114.

    [16]

    LIU J, CAO P, HAN D Y. Sequential indentation tests to investigate the influence of confining stress on rock breakage by tunnel boring machine cutter in a biaxial state[J]. Rock Mechanics and Rock Engineering, 2016, 49(4): 1479–1495.

    [17]

    LIU Q, LIU Q S, PAN Y C, et al. Experimental study on rock indentation using infrared thermography and acoustic emission techniques[J]. Journal of Geophysics and Engineering, 2018, 15(5): 1864–1877.

    [18] 李克金, 李文, 张德文, 等. 基于空腔膨胀理论的TBM滚刀破岩模型研究[J]. 武汉大学学报(工学版), 2020, 53(7): 583–590. https://www.cnki.com.cn/Article/CJFDTOTAL-WSDD202007003.htm

    LI Ke-jin, LI Wen, ZHANG De-wen, et al. Study on rock breakage model of TBM disc cutter based on cavity expansion theory[J]. Engineering Journal of Wuhan University, 2020, 53(7): 583–590. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-WSDD202007003.htm

    [19] 孙伟, 张旭, 赵奎山. 基于密实核理论的单滚刀多阶段受力预测模型[J]. 机械设计与制造, 2015(6): 9–12. https://www.cnki.com.cn/Article/CJFDTOTAL-JSYZ201506004.htm

    SUN Wei, ZHANG Xu, ZHAO Kui-shan. Multi-stage force prediction model of single disc cutter based on the dense nuclear theory[J]. Machinery Design & Manufacture, 2015(6): 9–12. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSYZ201506004.htm

    [20]

    HUO J Z, WANG W Z, SUN W, et al. The multi-stage rock fragmentation load prediction model of tunnel boring machine cutter group based on dense core theory[J]. The International Journal of Advanced Manufacturing Technology, 2017, 90(1/2/3/4): 277–289.

    [21] 孙剑萍, 胡瑜涛, 汤兆平, 等. 基于正交试验法的盾构机盘形滚刀磨损影响因素研究[J]. 中国科技论文, 2018, 13(10): 1158–1163. https://www.cnki.com.cn/Article/CJFDTOTAL-ZKZX201810012.htm

    SUN Jian-ping, HU Yu-tao, TANG Zhao-ping, et al. Influences on wear of disc cutter in the shield machine based on orthogonal test method[J]. China Sciencepaper, 2018, 13(10): 1158–1163. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZKZX201810012.htm

    [22] 林赉贶, 夏毅敏, 贾连辉, 等. 安装参数与掘进参数对滚刀破岩阻力的影响[J]. 浙江大学学报(工学版), 2018, 52(6): 1209–1215. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC201806021.htm

    LIN Lai-kuang, XIA Yi-min, JIA Lian-hui, et al. Influence of installation and tunnelling parameters on rock-breaking resistance of disc cutter[J]. Journal of Zhejiang University (Engineering Science), 2018, 52(6): 1209–1215. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC201806021.htm

    [23] 俞茂宏, 昝月稳, 范文, 等. 20世纪岩石强度理论的发展: 纪念Mohr-Coulomb强度理论100周年[J]. 岩石力学与工程学报, 2000, 19(5): 545–550. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200005000.htm

    YU Mao-hong, ZAN Yue-wen, FAN Wen, et al. Advances in strength theory of rock in 20 century—100 years in memory of the Mohr- coulomb strength theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2000, 19(5): 545–550. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200005000.htm

    [24] 谭青, 易念恩, 夏毅敏, 等. TBM滚刀破岩动态特性与最优刀间距研究[J]. 岩石力学与工程学报, 2012, 31(12): 2453–2464. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201212010.htm

    TAN Qing, YI Nian-en, XIA Yi-min, 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. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201212010.htm

    [25]

    YIN L J, GONG Q M, MA H S, et al. Use of indentation tests to study the influence of confining stress on rock fragmentation by a TBM cutter[J]. International Journal of Rock Mechanics and Mining Sciences, 2014, 72: 261–276.

  • 期刊类型引用(5)

    1. 黄仟. 铁路隧道TBM破碎岩石机理及破岩力影响因素研究. 江苏建筑职业技术学院学报. 2024(01): 26-31 . 百度学术
    2. 张魁,陈振宇,杨长,郑学军,张旭辉. 基于小孔应力集中的激光预钻孔条件下TBM滚刀侵岩模型建模与试验验证. 岩土工程学报. 2024(06): 1279-1288 . 本站查看
    3. 王帅,郭庚鑫,李阁强,董振乐,李东林,韩伟锋,焦雷,宋斌. 基于ADAMS的全断面竖井掘进机锥形刀盘载荷特性. 科学技术与工程. 2024(17): 7122-7130 . 百度学术
    4. 汪鼎华,江彪,徐国强,张魁. 移动激光预切槽辅助TBM滚刀破岩的仿真研究. 陕西科技大学学报. 2024(06): 151-156+179 . 百度学术
    5. 孙浩凯,高阳,朱光轩,徐飞,郑新雨. 隧道掘进机滚刀破岩动态荷载理论模型及试验研究. 岩土力学. 2023(06): 1657-1670 . 百度学术

    其他类型引用(3)

图(9)  /  表(6)
计量
  • 文章访问数:  134
  • HTML全文浏览量:  44
  • PDF下载量:  22
  • 被引次数: 8
出版历程
  • 收稿日期:  2021-10-10
  • 网络出版日期:  2022-12-13
  • 刊出日期:  2022-11-30

目录

/

返回文章
返回