随着人类对地下空间利用规模加大,竖井被广泛应用于地下矿产资源开采、交通隧道和废弃物质封存等工程中(
刘翔宇等,2024;
杨小聪等,2024)。然而,竖井施工常见的方法包括爆破法和机械法,目前多采用爆破法进行施工。爆破法是以钻眼爆破为基础,需大量下井人员,工作条件艰苦且作业环境危险系数大,由于工序机械化程度低,只能通过人海战术提高掘进速率(
Jin et al,2025)。随着金属矿山开采向深部和高强度围岩条件推进,在高地应力、硬岩及复杂构造条件下传统爆破施工面临着安全性差和成井质量难以控制等问题(
唐宇等,2021)。因此,寻找合适的替代方法已成为岩石工程界面临的主要问题,为此具有破岩安全性好、效率高和投入小等优点的机械法施工成为竖井施工技术的发展方向(
王少锋等,2021;
刘志强等,2022;
付鸿鑫等,2023)。
机械法中的竖井掘进机作为用于竖井施工的现代化复杂掘进装备具有明显优势,作业时由上而下在地层中钻进,为矿体开采和地下空间建设提供了通道。竖井掘进机采用机械破岩的方式将岩石从岩体上连续分离出来,依靠机械或自重掉落方法排出岩渣,利用辅助吊盘等配套设备完成井筒支护,依靠地面上的提绞设备完成人员、材料和破碎岩石的提升(
吕旦等,2023)。相对于爆破法和传统机械法,竖井掘进机施工方法可极大地减少竖井内作业人员数量,降低工人劳动强度,改善作业环境。针对目前竖井施工现有设备及工艺存在的问题,自主研制了一种重力式竖井掘进机,该设备对地层适应能力强,可实现连续掘进,提高施工效率。在重力式竖井掘进机中,截齿侵入角度是制约重力式竖井掘进机截齿滚动破岩效率和能耗水平的关键参数,直接影响截齿—岩体接触状态及岩石破裂方式。因此,系统研究其破岩机理对优化掘进参数具有重要意义。
张鑫(2022)通过可控试验研究系统揭示了切入角度等参数对盘形滚刀截割性能的影响规律,明确了侵入角度变化对截割效率与比能耗的主导作用。
武茂超等(2023)基于多工况现场数据分析指出,截齿破岩效率对截割角度变化高度敏感,且法向力与切向力比值随截割角度的增大呈现先增大后减小的非线性关系。
刘先珊等(2023)从理论力学角度建立了截齿破岩受力模型,阐明了侵入相关角度对破岩力大小及分布特征的控制作用,为截齿结构参数优化提供了理论依据。
肖建成等(2022)采用数值模拟方法揭示了贯入角度对截齿破岩效果的影响规律,表明较小贯入角度有利于增强岩体破碎程度。以上研究均是对截齿以不同角度垂直侵入岩进行研究,对于截齿基于设备自身重力及掘进滚筒旋转产生的旋转力与惯性力作用下以不同角度旋转破岩的研究鲜有涉及。
因此,针对重力式竖井掘进机独特的破岩方式,本研究基于MultiFracS有限—离散元(FDEM)数值模拟软件和相似理论,采用数值模拟与相似试验相结合的方法,对滚动不同截齿侵入角度下掘进滚筒旋转破岩展开一系列研究。通过上述研究,为实际竖井施工提供有效参考,为重力式竖井掘进机的研发提供理论基础。
1 滚动截齿破岩FDEM分析
1.1 截齿破岩模型建立
重力式竖井掘进机由前置掘进系统、后置驱动系统、液压系统和控制系统组成。其利用设备的自身重力及掘进滚筒旋转产生的旋转力和惯性力,通过掘进滚筒传递给截齿,即随着掘进滚筒高速旋转,截齿在机身自重力及掘进滚筒旋转产生的旋转力和惯性力共同作用下以不同角度旋转破碎岩石。截齿在岩石中往复旋转破碎,由于插入的路径与拔出的路径存在一定偏角,使得作用区域的岩石在滚动截齿作用下发生剪切和拉压破坏。
滚动截齿破坏岩石是一个极为复杂的过程,通过探究岩石在单个截齿不同作用下产生不同程度的损伤,科学描述了单个截齿对岩体的破坏,进一步探究了单截齿滚动作用下岩石破坏特性。在构建单截齿滚动破岩有限元模型时,将截齿简化为刚性截齿。本文采用MultiFracS有限—离散元(FDEM)数值模拟软件,基于滚动截齿破岩连续性建立了相应的滚动截齿作用下岩石破坏FDEM数值模型。通过建模软件Gmsh将截齿破岩模型划分为三节点三角形单元实体有限元网格,并在实体单元的边界上插入初始零厚度的四节点节理单元,从而实现模型的离散化,所构建的单截齿破岩模型如
图1所示。考虑到模拟掘进滚筒旋转360°破岩的工作存在耗时过长的问题,为兼顾计算效率和精度(
Rostami et al,2020;
Guo et al,2023),对不同区域岩石采用不同密集程度网格进行划分,即截齿与岩石接触区域的网格进行加密细化。为了真实有效地模拟竖井中的施工环境,对岩石试件左右下表面施加约束条件,上表面为自由面。
1.2 材料标定
截齿采用新型钛合金材料,密度为7 900 kg/m
3,弹性模量为270 GPa,泊松比为0.3。选取新乡某矿西风井地层较为常见的砂岩作为试验对象;根据一维条件、均匀性假设和试验设备的尺寸要求,将岩石试样的尺寸(直径×高度)设计为50 mm×100 mm和50 mm×25 mm的标准尺寸的圆柱体,如
图2所示。
如
图3所示,分别开展岩石单轴压缩试验和巴西劈裂试验,测得砂岩力学特性参数用于岩石有限—离散元模型的细观参数匹配。最终测定岩石的抗压强度为45.4 MPa,抗拉强度为7.4 MPa。
参数标定过程中有诸多调整策略,合适的方法能够大大减轻标定工作量。为此,学者们普遍认为可以通过单轴压缩和巴西劈裂试验结果反演进行标定,当某一组参数条件下模拟得出的荷载—位移曲线与试验结果一致时,认为该组参数是合理的,该试验方法比较成熟,已成功应用于数值仿真中参数标定环节(
任义等,2022;
张开雨等,2022)。对岩石试样开展单轴压缩和巴西劈裂试验模拟,如
图3所示。对比分析裂纹情况可知断裂位置基本一致,因此仿真对应的细观参数能够很好地反映真实砂岩的宏观参数,说明标定的颗粒细观参数是合理的。
对单轴压缩试验和巴西劈裂试验的数值仿真与室内试验数据进行统计,结果如
表1所示。
从
表1可以看出,单轴压缩强度的数值仿真结果与试验结果相近,数值试验所得值为45.2 MPa,二者相对误差仅为0.4%;岩石抗拉强度的数值试验值为7.3 MPa,与室内试验值的相对误差为1.35%。单轴压缩试验与巴西劈裂试验结果的数值试验均与室内试验参数差距小,验证了FDEM方法的准确性,说明该组数据可作为FDEM模拟中岩石模型参数的依据。最终获得离散元仿真模型中所采用的细观参数如
表2所示。
1.3 不同侵入角度对截齿滚动破岩机理的影响
为验证截齿侵入角度对单截齿破岩行为的影响,通过控制变量法,即保持截齿侵入深度为5 mm,滚筒旋转速度为35 r/min,对不同侵入角度下单截齿破岩进行仿真模拟分析,侵入角度为滚筒滚动过程中截齿与岩石接触角度,因此对侵入角度为30°、45°、60°、75°和90°条件下进行数值仿真研究。
(1)不同截齿侵入角度下岩石破坏状况
不同滚动截齿侵入角度条件下岩石最终破坏状况和裂纹演变如
图4所示,随着滚动截齿侵入角度的增加,岩体内部裂纹数量逐渐增加,随着微裂纹的产生与扩展,与滚动接触的区域破碎范围逐渐增大但内部裂纹未完全贯通,所以破碎体积逐渐变小。其中滚动截齿在一定时间内侵入角度为30°、45°、60°、75°和90°时分别对应截齿滚动破岩距离为32.62,35.63,39.72,45.17,56.25 mm,说明随着截齿侵入角度的增加,裂纹发育得越远越深,截齿滚动破岩距离越大。同时,岩石破碎范围呈现出非对称状态,这是由于随着滚动截齿侵入角度的增加,截齿与岩石接触面积不断减小,旋转产生的旋转力部分转变为摩擦力,截齿对岩石的挤压与摩擦力随侵入角度的增大而不断减小,岩石凹陷深度和破坏面积增加,破坏效果显著。
为了进一步探究不同侵入角度下单截齿滚动破岩的效果,对岩体内部微裂缝数量进行统计分析,结果如
图5所示。从
图5可以看出,脆性材料的岩石随着截齿侵入角度的增大,裂纹数逐渐增加,但是增加幅度先增大后减小,拉伸裂纹数量增加幅度较小,主要是剪切裂纹数量增加。
(2)不同侵入角度下截齿破碎比能变化
不同截齿侵入角度下截齿作用力变化趋势如
图6所示,随着截齿侵入角度的增加,截齿作用力呈现出上升趋势,截齿破碎岩石之后产生大量岩石碎块,由于岩体已发生干扰破坏,作用力逐渐下降,随着截齿侵入深度的不断增加,密实核产生时间无法确定,因此截齿破岩作用力呈现出波浪式浮动变化趋势。
对不同截齿侵入角度下单截齿破岩作用力和岩体破碎体积进行统计,结果如
图7所示。由
图7可知,随着截齿侵入角度的不断增加,截齿作用力逐渐减小,岩体破碎体积呈现出小幅上升趋势,由此可见,在不考虑其他因素的条件下,重力式竖井掘进机滚动截齿与岩体的接触面积随着侵入角度的减小而不断提高,截齿对岩体的挤压与摩擦力随侵入角度的减小而不断增大。
破碎比能是评判截齿破岩效率的重要参数之一,其定义为截齿破碎单位体积岩石所消耗的能量,破碎比能越低,效率越高(
陈新明等,2022;
张将令等,2022;Chen et al,2025)。其计算公式为
式中:SE为破碎比能(J/mm3);F为截齿作用(kN);h为截齿侵入深度(mm);V为岩石体积(mm3)。
根据
式(1)计算数值模拟下的破碎比能,对不同截齿侵入角度下破岩数据进行统计,结果如
表3所示。由
表3可知,随着截齿侵入角度的增加,单截齿破碎比能呈现下降趋势,较小的破碎比能对应较高的破岩效率,因此截齿侵入角度为90°时截齿破岩效率最高。
综上所述,在其他条件不变的前提下,随着截齿侵入角度的增加,拉伸裂纹起主要作用,剪切裂纹的数量随之增加,岩石破碎体积随之增大的幅度较小,截齿破岩作用力随之减小,破碎比能逐渐降低,截齿破岩效率越高。在重力式竖井掘进机施工中,截齿侵入角度对掘进机破岩性能的影响机制对于揭示前置掘进系统的几何排布规律与设计方法具有重要意义,因此尽量选取最优角度,提高竖井掘进效率。
2 试验验证
为了进一步探究重力式竖井掘进机不同滚动截齿侵入角度作用下岩石破坏机理,对滚动截齿破岩展开一系列研究。由于重力式竖井掘进机目前仍处于研发阶段,因此采用简单可控的室内试验进行研究,有效避免了其他无关因素的干扰,更加直观地再现出滚动截齿作用下岩石破坏特性。根据相似性原理,试验模型的尺寸缩放并不影响最终得到的结果,其研究对象的行为与真实情况存在一定相似性(
Chen et al,2023;
刘大兵等,2024)。因此,通过相似试验,对滚动截齿侵入岩石过程进行分析,有助于全面系统地研究不同滚动截齿侵入角度作用下岩石破坏机理问题。
2.1 试验试样及力学性能
在试验前应选取与试样相似的材料,相似材料的选择对模型试验的成功起着至关重要的作用(
邓嵘等2023;
孙超等,2024)。选取的材料应具有与原材料相似的物理力学性质,综合考虑下选取符合相似试验理论的混凝土材料作为试验对象;混凝土制作工序简单,在凝固过程中变形较小,且与岩石材料构造相似,在一定程度上能够真实模拟岩石动态破坏特性(
欧阳涛,2014;
潘文龙等,2023;
史勇等,2024)。水泥选用普通硅酸盐水泥P·O 42.5,由水泥、水、砂按照1.0∶0.5∶3.8进行配比,养护28 d得到,采用标准试块模制作完成混凝土制样,并经WAW-600型材料试验机如
图8所示获得试件的力学参数,其中测试试样所得力学参数如
表4所示。
2.2 截齿滚动破岩试验
(1)试验目的及内容
截齿滚动破岩过程是截齿与岩石之间角度不断变化的过程,为满足试验要求,本试验采用WAW-600型材料试验机对不同角度岩石试样下单截齿破岩规律进行研究。在试验过程中,观察岩石裂纹演变过程,分析截齿破岩力和岩石破碎块度,进而得出单截齿滚动破岩的一般性结论。
(2)试验方案
为了有效地模拟出截齿滚动破岩过程,采用截齿对不同角度岩石试样进行侵入试验,为提高试验效率,选出具有代表性的截齿和岩石角度进行研究,选取角度为30°、45°、60°、75°和90°,采用钢板对岩石四周施加相同约束条件。试验时将装有不同角度试块的岩石盒放置在试验台上,通过压力装置将固定于加载装置上部钢板之上的截齿逐步加载至给定侵入深度后卸载,如
图9所示。在该过程中利用高速摄像仪观察与截齿接触区域处岩石瞬时破坏过程,捕捉岩石破碎块碎裂过程,并记录载荷变化曲线。
(3)试验结果与分析
对不同截齿侵入角度下岩石破坏状况进行统计,结果如
图10所示。可以看出,随着截齿侵入角度的增加,岩石凹陷深度和破坏面积增加,裂纹发育得越远越深,同时发现岩石破碎范围呈现出非对称状态。这是由于随着滚动截齿侵入角度的增加,截齿与岩石的接触面积不断减小,掘进滚筒旋转产生的旋转力与惯性力部分转变为摩擦力,截齿对岩石的挤压与摩擦力随侵入角度的增大而不断减小,破坏效果显著。
通过计算机采集仪对不同角度下截齿破岩所需作用力进行统计并绘制变化曲线,如
图11所示。由
图11可以看出,当截齿侵入角度为30°时,截齿作用力为23.02 kN;当截齿侵入角度为45°时,截齿作用力为19.73 kN;当截齿侵入角度为60°时,截齿作用力为17.77 kN;当截齿侵入角度增大到75°时,截齿作用力为16.10 kN;当截齿侵入角度增大90°时,截齿作用力为15.06 kN。可以看出,随着滚动截齿侵入角度的增加,截齿作用力降低幅度逐渐减小。
随着截齿侵入角度的增加,截齿作用力逐渐减小。由此可见,在不考虑其他因素的条件下,重力式竖井掘进机滚动截齿与岩石接触面积随着侵入角度的减小而不断提高,滚动截齿对岩石的挤压与摩擦力随侵入角度的减小而不断增大。上述结果与数值模拟中所得结论一致,验证了数值模拟结果的准确性。
3 结论
(1)由FDEM方法标定的单轴压缩和巴西劈裂试验模拟结果与室内试验数据一致性较高,验证了FDEM方法的正确性。同时,也验证了滚动截齿破岩数值模拟结果的可靠性。
(2)岩石在滚动截齿作用下产生大量由拉伸和剪切破坏引起的裂纹,随着截齿侵入角度的增加,拉伸裂纹起主导作用,剪切裂纹数量随之增加;岩石破碎体积随之增大的幅度较小,截齿破岩作用力则随之减小;由破碎比能计算公式得出,当截齿侵入角度为90°时,破碎比能最低(0.05 J/mm3),截齿破岩效率最高。
(3)相似试验与模拟试验均具有较好的模拟精度,表明基于MultiFracS有限—离散元数值模拟软件开展截齿滚动破岩相关模拟试验具有可行性。基于相似理论的滚动截齿破岩试验探究了岩石在滚动截齿作用下的破坏演化过程,研究了截齿侵入角度对截齿破岩的影响。研究表明截齿破岩所需作用力随截齿侵入角度的增大而减小,且岩石破碎范围呈现出非对称状态。
国家重点研发计划项目“金属矿深部规模化安全高效连续智能开采技术与装备”(2023YFC2907203)