截齿破岩力学特性与碎片分离数值模拟研究

2022-10-20万丽荣高冠顺逯振国王志文马玉华

煤炭科学技术 2022年9期

万丽荣，高冠顺，逯振国，王志文，马玉华

(山东科技大学机械电子工程学院，山东青岛 266590)

0 引言

目前，煤炭开采工作面采高逐渐增高，开采技术及开采装备不断发展，巷岩掘进速度不断加快[1-2]。巷岩掘进设备主要以悬臂式掘进机为主，在截割过程中截齿承受较大冲击载荷，因此截齿是最容易磨损和损坏的部件，进而影响整个掘进机的工作性能[3-5]。为提高掘进机工作效率，国内外众多学者针对截齿截割机理及性能做了大量研究。王想等[6]研究了相关角度参数对截割力和比能耗的影响。杜鑫等[7]利用ANSYS/LS-DYNA，分析了不同截割线速度截齿受力，得出随着截割速度的增加，截齿所受最大峰值应力逐渐减小，随着线速度的继续增加，最大峰值应力减小的幅度变小。刘春生[8]探究了截割力谱的分形特征与安装角及切削厚度的关系。LIU等[9]利用PFC(2D)研究了岩石物理参数对岩石破坏和截割性能的影响。WANG等[10]研究了切削参数对截割比能耗及切削力的影响，并提出比能的一般模型。DAI等[11]和SU等[12]利用DEM(3D)来研究破岩过程，其研究的共同特点是碎片的分离为颗粒形式。LI等[13]利用LS-DYNA(3D)研究了不同截割头截割岩石力的变化情况，结果表明不同安装角下的截割力不同。毛君等[14]研究了截齿安装角对截割性能的影响得出安装角为45°时截割效率最高。王立平等[15]研究了基于Evans截割模型镐型截齿截割力的计算，得出了截割力计算公式。王杰等[16]基于单元破裂的方式，构建了连续-非连续单元法，并应用于岩石裂纹扩展问题。以上研究建立了多种煤岩破碎力学模型，对部分截割参数对截割性能的影响做了研究；对于镐型截齿破岩过程中截割参数对裂纹扩展和碎片分离的影响研究不足，对于裂纹和碎片分离，前人的研究多以试验为基础，从截齿受力规律与岩石破碎形式出发，对于运用三维有限元法数值模拟岩石的破碎研究较少。

数值模拟研究方面，虽然二维数值模拟方法在获得裂纹扩展和碎片分离方面具有优势，但获取的截割力与实验结果的一致性较差。截齿破岩实质上是三维的切削过程，利用三维数值模拟方法呈现截齿破岩中岩石的裂纹扩展和碎片分离过程，对揭示截齿破岩机理具有重要意义。以EBZ-160TY型掘进机截割头上的截齿为例，采用ANSYS/LS-DYNA(3D)建立截齿截割岩石有限元模型，通过结合侵彻失效和压缩损伤，研究在不同截割角及切削厚度下岩石裂纹扩展和碎片分离过程，并对截割力进行分析。

1 数值模拟模型的建立

1.1 模型建立及单元的选择

运用三维绘图软件SolidWorks建立截齿截割岩石的几何模型，然后导出.x_t格式文件，并导入到ANSYS软件中进行单元类型的设置和网格的划分[17]。本次数值模拟采用阶梯型镐型截齿，截齿的结构参数主要包括截齿长度L=160 mm、齿身长度Lj=80 mm、齿柄直径Djb=38 mm、齿尖夹角αjj=80°、齿尖合金头直径djh=16 mm，其结构如图1所示。图2为截齿截割岩石模型示意图，其中，岩石的长、宽、高分别为220、200、60 mm，V为截齿沿X轴正方向的截割速度，截齿截割角α为截齿中心线与X轴反方向间夹角，h为切削厚度。模型单元类型选取SOLID 164实体单元对截齿和岩石模型进行网格划分，采用单点积分以节省数值模拟时间。

图1 截齿结构示意Fig.1 Structure drawing of pick

图2 截齿截割岩石示意Fig.2 Schematic of rock cutting with pick

1.2 模型约束和截割参数设置

截齿沿X轴正方向运动，同时约束截齿的其他自由度。给岩石添加约束，岩石底面为全约束，X方向的2个侧面添加Y方向位移约束，Y方向2个侧面添加X方向位移约束。为缩短计算时间通过在岩块模型的非截割面添加非反射边界条件，采用小块岩石来模拟无限大岩块[18]。

1.3 定义接触和载荷施加

主要研究岩石裂纹的形成与碎片的剥落，在数值模拟过程中不需要考虑到截齿的变形问题且截齿弹性模量远远大于岩石材料弹性模量，为进一步提高数值模拟效率，将截齿材料定义为刚体[19]。岩块使用HJC本构模型，岩块及截齿材料物理力学参数见表1，选择岩块和截齿接触类型为ERODING_SUEFACE_TO_SURFACE。单元失效根据HJC材料中岩石材料属性和添加*MAT_ADD_EROSSION工作语句作为失效判别依据[20]。当单元模型中的应力、应变超过预设值时，单元失效并删除，多个删除的单元连接一起便形成了裂缝。

对截齿施加运动载荷，改变截割角及切削厚度，研究岩石的破碎及失效情况，建立的有限元模型如图3所示，截齿截割岩石过程中截齿受力如图4所示，在三向力中截割力整体变化趋势与合力相近[21]，本文研究三向力中的截割力，FN为牵引阻力，FS为截割力。

图3 截齿截割岩石模型Fig.3 Rock cutting model with pick

图4 截齿受力示意Fig.4 Force of pick

表1 截齿岩石材料物理力学参数

1.4 岩石本构模型

HJC模型是HOLMQUIST等[22]提出的材料损伤模型,该本构模型应用于解决岩石等材料的大变形问题。此模型主要由3个方程构成，分别为屈服面方程、压缩状态方程和损伤演化方程。HJC本构模型损伤演化方程如下：

(1)

HJC屈服面方程即强度方程，代表了无量纲静水压力和等效应力在不同损伤状态下的联系，方程式如下：

σ*=[A(1-D)+Bp*N][1+C(lnε*)]

(2)

式中：C为应变率敏感系数；D为损伤系数；A为准凝聚强度；B为准强度增大系数；p*为无量纲静水压力；ε*为无量纲应变率；σ*为等效应力；N为压力增大系数。

HJC压缩状态方程采用如图5所示曲线方程来表示静水压力和体积应变的关系；压力和体积的函数关系共分为3个阶段：线弹性阶段、过渡阶段、密实阶段。

图5 静水压力-体积应变曲线Fig.5 Hydrostatic pressure-volume strain curve

第1阶段：q

(3)

式中：q为压碎压力；u为体积应变；K为体积弹性模量。

第2阶段：qc

第3阶段：q1≤q为密实阶段，材料中的空隙被完全压实，此时压力和体积应变关系为

(4)

2 仿真结果分析

运行ANSYS/LS-DYNA数值模拟软件，截齿以恒定截割速度运动，当与岩石接触时，组成岩块的单元与截齿产生相互作用力，当岩石单元所受应力大于岩石强度时，部分岩石单元失效删除，岩石开始生成裂纹，随着截齿运动，裂纹在岩石表面和内部贯穿、连接。在截齿破岩过程中，得到裂纹的形成过程如图6所示，改变截割角和切削厚度得到裂纹扩展和碎片分离的数值模拟结果如图7、图8所示，不同截割角及切削厚度下，采集截齿破岩过程中截齿截割力的变化作为纵坐标，以截割时间、切削厚度作为横坐标，得到的仿真结果如图9—图11所示。

图6 裂纹及碎片分离过程Fig.6 Crack and fragment separation process

图7 不同截割角碎片分离情况Fig.7 Fragment separation at different cutting angles

为便于分析对裂纹的影响，将截齿隐去并选取截割角为50°，切削厚度为20 mm的数值模拟结果来研究岩石整个裂纹生成及碎片分离过程。如图6a所示，由于一些单元在拉伸应力和剪切应力作用下失效，在截齿周围的岩石上产生裂纹且裂纹是随机分布的。如图6b所示，随着截齿在X方向上的位移增加，裂纹在拉伸应力的作用下向外扩展。如图6c所示，更多的裂纹不规则地出现，并且形成一些断裂区域。裂纹在岩石的表面和内部传播并相互连接，碎片从基岩分离。如图6d所示，随着截割位移的增加，具有随机数量和形态的碎片从基岩中分离出来；同时，在碎片与基岩之间的大量单元由于拉应力而失效。

块岩数是衡量破岩效率的一个重要标准，块岩数越多，破岩效率越高，且块岩数越多，粉尘量越小，对改善工人井下作业环境和提高破岩时的安全性具有重要意义；以碎片含有的单元个数作为判别的标准，当碎片含有的单元个数大于10时作为块岩标准。为研究截割角、切削厚度和块岩数之间的关系，对不同截割角及切削厚度下截齿破岩做了数值模拟分析，结果如图7、图8所示。

图7为截割速度3 m/s，切削厚度20 mm，截割位移60 mm时，不同截割角下的裂纹扩展和碎片分离情况，其中图7a为35°截割角时裂纹及碎片生成情况，截割角越小，越有利于截齿侵入岩体，但几乎无碎片生成。图7b中截割角为40°，相对截齿的锋利程度下降，开始有微小碎片和裂纹生成。图7c为截割角45°时截齿截割岩石破碎情况，与35°相比有了明显提升，在截割相同位移时，有碎片生成。图7d截割角为50°时截割情况，与上述3种情况相比，裂纹数目增多，块岩数增加，岩石破碎效率变大。图7e截割角55°时裂纹数目和碎片数量减少，破岩效率有所下降。图7f截割角60°时，碎片数量持续减少，几乎无碎片生成。

在截割速度3 m/s，截割角50°时，不同切削厚度下的数值模拟结果如图8所示，图8a中切削厚度为10 mm时，裂纹生成不明显，无碎片生成。图8b中切削厚度由10 mm增加到15 mm，截齿周围的单元由于拉应力和剪切应力的作用而失效，开始生成裂纹，并伴随着碎片的生成。图8c切削厚度为20 mm时，裂纹扩展明显，碎片数目增多。

由图7、图8仿真结果统计得出不同截割角下块岩数目见表2，不同切削厚度下块岩数目见表3。从表2中可以看出，截割角50°时块岩数最多，切削厚度越大，破岩效率越高。

表2 不同截割角下块岩数目

表3 不同切削厚度下块岩数目

图9为当截割深度为20 m时，单个截齿截割角分别为35°、40°、45°、50°、55°、60°时截齿受力曲线，由图中6条受力曲线可知，截割力最大峰值出现在截割角为35°时，约32.7 kN，截割角为50°时，最大峰值截割力约25 kN。分析表4中的数据可以得出，截割力最大峰值随着截割角的增大，呈先减小后增加的趋势。