含夹矸煤层滚筒磨损特性离散元模拟试验分析

赵丽娟,金 鑫,赵宇迪,王 斌

(辽宁工程技术大学 机械工程学院,辽宁 阜新 123000)

摘 要:复杂煤层赋存条件下工作的采煤机,其螺旋滚筒所受载荷具有非线性、时变性、冲击性和强耦合性,磨损问题不可避免,且严重影响其工作性能和生产效率。基于虚拟样机技术,采用离散元数值模拟方法研究含夹矸煤层条件下螺旋滚筒的磨损特性,以提高其截割性能、延长其使用寿命。以煤岩截割机理和离散元素法为基础,分析了滚筒与煤岩体间相互作用关系及滚筒磨损模型;以MG2×55/250-BWD型采煤机螺旋滚筒为工程对象,采用Hertz-Mindlin with bonding黏结模型为颗粒间接触模型、Hertz-Mindlin with Archard Wear模型为煤岩体与滚筒间接触模型,利用Proe/E与EDEM建立滚筒截割含夹矸煤层耦合模型;采用正交试验法,分析采煤机牵引速度、滚筒转速、螺旋升角和煤岩抗压强度对滚筒磨损问题的影响趋势,经过多元回归分析得到截齿、叶片磨损深度的评价模型;基于实际工况,综合考虑采煤机截割比能耗、载荷波动系数和生产率,构造多目标优化函数,利用MATLAB求解得到滚筒最优工作参数。基于EDEM离散元仿真分析结果,发现滚筒组件存在不同程度磨损,其中截齿齿尖局部接触区域磨损最为严重,叶片主要磨损区域为尾片齿座根部位置附近。煤岩抗压强度、采煤机牵引速度增加,截齿、叶片磨损加剧;滚筒转速增加,截齿磨损深度减小,叶片磨损深度先增加后减小;螺旋升角增加,截齿磨损深度先增加后减小,叶片磨损深度逐渐减小;其中,煤岩抗压强度是造成截齿、叶片磨损的主要因素。当截割煤岩抗压强度为30 MPa的煤层时,采煤机牵引速度为5.66 m/min、滚筒转速为80.40 r/min、螺旋升角为11.35°时,滚筒性能达到最优,此时采煤机截割比能耗为0.713 2 kW·h/m3、载荷波动系数为0.081 9、截齿、叶片和端盘的磨损深度分别为1.595 μm,1.061 7 nm和0.907 nm,生产率为315.465 5 t/h。离散元仿真结果中滚筒磨损部位分布与现场工况下滚筒实际磨损情况具有一致性,验证了建立的仿真模型与过程的可行性。基于此研究方法可实现滚筒磨损问题的定量分析,极大地降低试验成本,有助于高效、强力、耐磨滚筒的研发。

关键词:螺旋滚筒;含夹矸煤层;磨损;离散元数值模拟;正交试验

滚筒作为采煤机的工作机构,截割赋存条件复杂煤层时极易出现磨损失效,给整机工作可靠性和高效截割造成严重影响[1],因此研究滚筒与煤岩相互作用过程中磨损特性,对提高其使用寿命意义重大。目前,已有国内外专家学者展开相关研究:刘晓辉[2]建立了镐型截齿与煤岩相互作用接触模型,以截齿与煤岩的干涉量为分析指标,研究截割形式、结构参数和工作角度对截齿磨损特性的影响规律。张强等[3-4]采用小波包分析、神经网络等方法,建立截齿磨损识别模型,实现截齿磨损状态的在线监测和磨损程度的快速准确识别;赵洪盛[5]提出一种改进的聚焦形貌算法,用于测量计算截齿磨损体积,准确率高达95%;DEWANGAN等[6-7]通过采用扫描电子显微镜和能量弥散X射线探测器深入研究截齿磨损部位,对截齿磨损机理进行了分类。史秀宝等[8]在MATLAB环境下模拟了滚筒齿座与煤壁间的干涉,找出了末端齿座磨损原因,有助于滚筒的结构优化;卜庆辉等[9]分析了滚筒齿座、螺旋叶片尾端磨损原因,并给出了改进方案;赵丽娟等[10-11]对截齿磨损、折断等失效形式进行总结分析,并从材料、加工工艺、结构设计等方面提出了相应改进措施。罗晨旭[12]提出一种基于SUSAN边缘检测算法的图像处理方案来获取截齿磨损面轮廓尺寸信息,研究了煤岩界面位置形式对截齿磨损的影响规律。

以上研究主要侧重于滚筒磨损机理和磨损程度识别等方面,对滚筒磨损程度的定量分析、因素影响规律等尚未进行深入探讨。笔者采用EDEM离散元数值模拟方法,以MG2×55/250-BWD型采煤机螺旋滚筒为工程对象,创建滚筒截割含夹矸煤层耦合模型,分析了滚筒磨损机理,研究了采煤机牵引速度、滚筒转速、螺旋升角和煤岩抗压强度对滚筒磨损的影响规律,并建立截齿、螺旋叶片磨损数学模型。基于实际工况,综合考虑截割比能耗、载荷波动系数、生产率,确定了滚筒最佳工作参数,为滚筒磨损问题研究提供了方法指导,对实际工况下滚筒磨损程度预测及截割参数选取提供参考,具有一定工程应用价值。

1 理论背景

1.1 镐型截齿力学模型

采煤机截割、破碎含夹矸煤层时,镐型截齿受力如图1所示。

图1 滚筒单齿受力分析[13]

Fig.1 Stress analysis of drum pick[13]

截煤时,截割阻力和牵引阻力分别[14]

(1)

Yj=(0.5~0.8)Zj

(2)

式中,ZjYj分别为截割阻力和牵引阻力,为平均截割阻抗,N/mm;bp为截齿宽度,cm;Kψ为煤的脆性系数;hmax为截齿最大切削厚度,cm;θ为截齿所处位置角度,(°);tcp为截齿切削宽度,cm;Kz为自由表面系数;Ky为截角影响系数;Kφ为截齿前刃面形状影响系数;Kc为截齿排列方式系数;Kot为地压对煤壁影响系数;φ为截齿对于牵引方向的偏转角,(°);f′为截割阻抗系数;δcm为煤的单向抗压强度,MPa;Sj为截齿磨损面在截割平面上的投影面积,m2;Kδ为矿体应力状态体积系数。

截割夹矸时:截割阻力和牵引阻力分别为

Zj=PKkTkψkψkdky(0.25+1.8hmaxsin θtcp)+

PK0.1Sj

(3)

式中,PK为岩石的接触强度,MPa;kT为截齿类型系数;kψ为硬质合金头形状系数;kψ为刀头部形状系数;kd为硬质合金刀头直径系数;ky为截齿截角影响系数。

Yj=2.5Zj(0.15+0.000 56PK)/(10hmaxsin θ)0.4

(4)

侧向力可由截割阻力表示为

(5)

式中,Xj为侧向力,N;Ci为截齿排列方式影响系数(i=1,2,3)。

1.2 螺旋叶片与煤岩颗粒间力学模型

滚筒装煤过程中,破落的煤岩颗粒随叶片共同运动,两者间的力学耦合作用十分复杂。由于颗粒质量较小,忽略其自身重力的影响[15],为简化研究,依据散体力学理论,对单体煤岩颗粒进行受力分析,如图2所示。

图2 螺旋叶片与煤岩颗粒间力学分析

Fig.2 Mechanical analysis of blade and coal-rock particle

煤颗粒平衡方程为

(6)

式中,Pt,Px分别为落煤切向力和轴向抛煤力,N;Nα为抛煤时螺旋叶片对煤颗粒的正压力,N;fm为螺旋叶片与煤颗粒的摩擦因数;β为叶片螺旋升角,(°)。

根据落煤切向力与煤岩颗粒切向分速度求得装煤功率为

NZ=0.001Nαvt(sin β+fmcos β)

(7)

式中,NZ为装煤功率,kW;vtnDsin βsin(β+ρm)/cos ρmρm为摩擦角,(°);vt为煤岩颗粒切向分速度,m/min;D为颗粒所在叶片作用处的平均回转直径,m;n为滚筒转速,r/min。

装煤功率亦可表示[16]

NZ=0.1vqvjKzm/n

(8)

式中,vq,vj分别为牵引速度和截齿截割线速度,m/s;Kzm为装煤阻力系数,N/cm。

联立式(7),(8)可得到煤颗粒的正压力为

(9)

1.3 磨损模型

截割、装载过程中,煤岩体作为硬质颗粒对滚筒材料表面不断挤压,并存在相对运动,造成滚筒严重磨损,由磨损机理可知两者间磨损类型为磨粒磨损[17]。EDEM中Hertz-Mindlin with Archard Wear模型是基于J.F.Archard磨损理论[18],对Hertz-Mindlin(no slip)模型的扩展,可用于设备磨损区域表面磨损深度值的求解,因此本文用以研究滚筒磨损问题。

根据J.F.Archard磨损模型可知,滚筒磨损量与载荷及材料之间的关系为

V=KFL/H

(10)

式中,V为磨损体积,mm3;K为磨损常数;F为接触载荷,N;L为相对滑动距离,mm;H为材料硬度。

定义W=K/H,则EDEM中滚筒表面磨损深度值为

h=WFL/A

(11)

式中,h为磨损深度,mm;A为接触面积,mm2

2 EDEM离散元仿真模型建立

虚拟煤壁模型的构建是滚筒磨损离散元数值模拟试验的关键,为使虚拟煤壁更加接近实际,保证模拟结果的可信度,对兖州矿区杨村矿煤层进行取样,通过试验(图3)测定煤岩的物理、力学性质见表1。

图3 煤岩物理、力学性质测定试验

Fig.3 Physical and mechanical properties experiments of coal-rock

表1 煤岩试样物理力学性质

Table 1 Physical and mechanical characters of coal-rock samples

煤岩17煤夹矸1夹矸2夹矸3弹性模量/MPa2 0103 62012 10018 300密度/(kg·m-3)1 2802 4602 6302 610抗压强度/MPa12304252内摩擦角/(°)20.4833.2532.0038.72泊松比0.280.240.230.21

设置煤岩颗粒半径为12 mm,颗粒间接触模型为Hertz-Mindlin with bonding黏结模型,黏结参数可根据物理、力学性质测定结果,基于Hertz接触理论[19],采用项目组开发的采煤机工作机构优化设计及载荷计算软件[20]计算得到,如图4所示。

图4 颗粒间黏结参数

Fig.4 Bonding parameters between particles

采用Pro/ENGINEER族表功能建立采煤机滚筒三维模型,进行网格划分后导入EDEM中,并根据需要建立夹矸层厚度为250 mm,整体尺寸为2 000 mm×800 mm×1 000 mm的煤壁模型[21]。设置颗粒与滚筒间接触模型为Hertz-Mindlin with Archard Wear模型。在Geometry模块下,通过Dynamics选项添加滚筒运动参数,根据现场工作条件,设置滚筒转速为72 r/min,牵引速度为4.2 m/min,最终得到离散元仿真耦合模型如图5所示。

图5 EDEM仿真耦合模型

Fig.5 EDEM simulation coupling model

3 数值模拟结果分析

3.1 煤岩颗粒运动分析

利用EDEM后处理模块查看滚筒截割过程中颗粒运动状态,如图6所示。由图6可知,在滚筒截割作用下,外侧煤壁被截落的颗粒首先落入装煤区;内侧煤壁破碎后的颗粒首先在齿座根部附近形成堆积,在自身重力及螺旋叶片的推挤作用下,一部分沿滚筒轴线流向刮板输送机,一部分做抛物运动越过筒毂落入采空区,形成浮煤,且离滚筒越远,颗粒运动速度越大。

图6 滚筒动态截割过程颗粒运动状态

Fig.6 Particles motion state in dynamic cutting process of the drum

图7 滚筒组件磨损分布及载荷曲线

Fig.7 Wear distribution diagram and load curve of drum assembly

3.2 滚筒磨损分析

滚筒受煤壁和截落煤岩颗粒的共同作用,存在不同程度磨损。为便于观察,对滚筒主要磨损部位进行局部放大处理,滚筒组件磨损分布图及载荷曲线如图7所示。

由图7可知,滚筒组件中,截齿磨损最为严重,载荷峰值为45.15 kN,最大磨损深度为2.35 μm,磨损区域主要为合金头的齿尖局部接触区(受载区域B)及齿柄头部的轴肩处(受载区域A),且受载区域B磨损程度大于区域A;叶片载荷峰值为12.03 kN,最大磨损深度为1.06 nm,磨损主要集中在尾片齿座根部附近位置;端盘载荷峰值为5.02 kN,最大磨损深度为0.907 μm,磨损主要集中在端盘端面和端部相邻齿座的空隙部位。这是由于截割过程中,截齿受载区域B与煤壁直接接触,存在较大接触力,并伴有相对运动,使得区域B受摩擦力作用明显,磨损最严重;滚筒装煤时,叶片与煤岩颗粒间存在推挤作用力,是其产生磨损的主要原因。而齿根附近煤岩颗粒的堆积使得该部位叶片受力较大,造成局部磨损相对严重;端盘位于滚筒的靠近煤壁侧,其与煤壁间的容煤空间较小。由端部截齿截落的煤岩颗粒部分落入容煤空间后受端盘和煤壁的不断挤压而继续破碎,部分在端部相邻齿座空隙部位堆积,由滚筒的不断转动随叶片排出,因此端盘端面与端部相邻齿座空隙部位受摩擦力作用较大,磨损程度较为严重。此外,滚筒和煤岩间复杂力学耦合作用产生的较大交变冲击载荷会引起滚筒的剧烈振动[22],给滚筒组件带来不同程度的附加磨损。

3.3 可行性验证

由于现场工况下,滚筒磨损试验用时较长,且受EDEM计算速度的限制,采用离散元模拟方法进行定量化验证存在很大难度。因此,根据文献[23-24]中仿真模型的验证方法及文献[25]中截齿磨损的试验研究,文中通过对比仿真结果与现场磨损中(图8)滚筒磨损部位和特点,发现滚筒磨损情况具有一致性,证明了采用EDEM模拟滚筒截割复杂煤层进行力学分析及磨损研究是准确可靠的,具有可行性。

图8 螺旋滚筒现场磨损情况

Fig.8 Wear condition of drum on site

4 滚筒磨损特性多因素影响分析

滚筒落煤与装煤是多因素耦合作用过程,其载荷具有非线性、时变性和强耦合性[26]。单一的理论分析并不能直接反映各因素对滚筒磨损的具体影响,笔者采用正交试验法探究不同工况下滚筒磨损规律。

4.1 试验设计

采煤机运动学参数、滚筒结构参数及煤岩物理、力学性质等共同决定滚筒的载荷特性[27],与滚筒磨损情况密切相关。选取采煤机牵引速度、滚筒转速、螺旋升角、煤岩抗压强度为设计因素,按照L9(34)正交表进行模拟试验,分析其对滚筒磨损的影响趋势,因素水平表见表2。

以滚筒截齿、叶片磨损深度为试验结果,利用EDEM离散元软件进行仿真模拟,试验方案及结果(取5位有效数字)见表3,其中A,B,C和D分别为v,n,βδ编码。

表2 因素水平[2]

Table 2 Factor level[2]

水平因素牵引速度v/(m·min-1)滚筒转速n/(r·min-1)螺旋升角β/(°)抗压强度δ/MPa1270103024901442361101652

表3 试验配置方案及结果

Table 3 Experimental configuration and results

试验编号因素编码ABCD试验结果截齿磨损深度hj/μm叶片磨损深度hy/nm111110.646 50.117 3212221.384 31.166 8313332.079 20.685 6421233.507 91.541 5522310.970 70.955 5623121.630 00.725 0731322.244 31.090 3832133.867 03.545 1933211.611 10.191 3

4.2 结果分析

4.2.1 极差分析

依据表3中的试验配置方案,对截齿、叶片磨损深度进行极差分析,结果见表4。其中,Ti为因素水平之和为因素水平均值;R为极差;S为均方差。

对比表4中截齿磨损数据:R(D)>R(A)>R(C)>R(B),S(D)>S(A)>S(C)>S(B),根据极差分析原理[29]可知,煤岩抗压强度是其磨损的主要影响因素,影响高度显著,牵引速度次之,滚筒转速的影响最小;对比叶片磨损数据,R(D)>R(B)>R(A)>R(C),S(D)>S(B)>S(A)>S(C),可知其磨损主要影响因素为煤岩抗压强度,影响最为显著,滚筒转速次之,螺旋升角的影响最小。

表4 截齿、叶片磨损深度极差分析

Table 4 Analysis of extreme difference in wear depth of pick and valve

因素截齿磨损深度hj/μmABCD叶片磨损深度hy/nmABCDT14.110 06.398 76.143 53.228 31.969 72.749 14.387 41.264 2T26.108 66.222 06.503 35.258 63.222 05.667 42.899 62.982 1T37.722 45.320 35.294 29.454 14.826 71.601 92.731 45.772 2T11.370 02.132 92.047 81.076 10.656 60.916 41.462 50.421 4T22.036 22.074 02.167 81.752 91.074 01.889 10.966 50.994 0T32.574 11.773 41.764 73.151 41.608 90.534 00.910 51.924 1R1.204 10.359 50.403 12.075 30.952 31.355 10.552 01.502 7S0.603 20.192 80.207 01.058 40.477 40.698 70.303 80.758 4

4.2.2 因素影响趋势分析

以各因素水平为横坐标,截齿、叶片磨损深度均值为纵坐标,经过数据处理得到各因素水平趋势图如图9所示。

图9 截齿、螺旋叶片磨损因素水平趋势

Fig.9 Factors trend diagram of pick and spiral blade wear

分析图9可知:截齿、螺旋叶片磨损深度均随牵引速度、煤岩抗压强度增加而增加,但增加趋势有所不同。煤岩抗压强度增加,截齿、叶片磨损深度呈快速上升趋势,牵引速度增加则呈相对平缓上升趋势;随滚筒转速增加,截齿磨损深度逐渐减小,而叶片磨损深度先增加后减小;随螺旋升角增加,截齿磨损深度呈先增加后减小趋势,而叶片磨损深度逐渐减小。

呈现上述变化趋势的原因是:牵引速度增大,单位时间内截齿截割厚度、滚筒内煤流质量增加,导致单齿受力、叶片与煤流间摩擦力增加,从而加剧截齿、叶片磨损程度;滚筒转速增加时,煤岩的破碎程度较高,细碎、小块状煤体增多,使得煤体与叶片的接触面积增大,摩擦现象显著;当滚筒转速增大到一定值时,截落的煤岩颗粒大部分会越过筒毂落向采空区,使得叶片输送的煤流质量减小,从而减轻磨损程度;螺旋升角改变时,截齿相对位置、截割顺序与截割时间发生变化,导致大块煤随机崩落,截齿受力呈不规律变化;螺旋升角增加时,煤流颗粒轴向速度增大,叶片输煤能力增强,叶片与煤流颗粒间作用力减小,磨损程度降低;煤岩抗压强度越大,滚筒破碎煤岩体越困难,因此磨损越严重。

4.2.3 最优方案确定

由表4中各因素的均值及图9可得出,在截割含夹矸煤层时,为降低滚筒磨损,因素A,B,C的最优水平分别为A1,B3,C3。采煤机实际工作中还应综合考虑生产率、载荷波动系数、截割比能耗等问题,因此需对最优方案进行调整。

以滚筒截割夹矸抗压强度为30 MPa的工况为例,基于SPSS软件,采用数理统计方法对正交试验结果进行多元回归分析,分别得到截齿、叶片磨损深度评价模型:

hj=-7.546 6+0.429 3v+0.045 39n+0.956 0β-

0.016 03v2-0.000 302n2-0.038 58β2

hy=-18.440 7+0.120 6v+0.514 2n-0.507 8β+

0.014 68v2-0.002 910n2+0.015 99β2

以牵引速度、滚筒转速、螺旋升角为设计变量,截割比能耗Hw、载荷波动系数σ、截齿磨损深度hj、叶片磨损深度hy最小,生产率Q最大为多目标优化函数,考虑到各单目标的权重问题,对多目标函数进行规范化处理,转化后的目标函数为

(12)

式中,各单目标函数具有同等重要地位,因此w1=w2=w3=w4=w5=0.2[30];采用项目组基于MATLAB及Excel联合开发的采煤机工作机构优化设计及载荷计算软件(2014SR102903),在其中输入相关参数,根据计算结果进行拟合得到截割比能耗Hw、载荷波动系数σ和生产率Q函数模型[31];Hwmax,σmax,hjmax,hymax,Qmax分别为各单目标函数在论域中最大值。

采用MATLAB中Constrained nonlinear minimization求解器进行求解,设置初始解为[2 70 10],迭代误差为10-6,得到优化函数最优解为(v,n,β)=(5.66,80.40,11.35),即当牵引速度为5.66 m/min、滚筒转速为80.40 r/min、螺旋升角为11.35°时,滚筒性能达到最优,此时截割比能耗为0.713 2 kW·h/m3、载荷波动系数为0.081 9、截齿磨损深度为1.595 μm、叶片磨损深度为1.061 nm、生产率为315.465 5 t/h。

5 结 论

(1)基于煤岩物理、力学特性测试试验,利用EDEM离散元软件建立采煤机滚筒截割含夹矸煤层耦合模型,对滚筒磨损过程进行数值模拟,仿真结果与滚筒现场磨损情况具有一致性,验证了采用EDEM离散元法进行滚筒磨损模拟试验研究的可行性。

(2)利用EDEM数值模拟方法揭示了滚筒磨损机理,得出截齿齿尖主要磨损区域为受载区域B,叶片主要磨损区域为尾片齿座根部位置附近。

(3)基于EDEM滚筒磨损数值模拟,对采煤机牵引速度、滚筒转速、螺旋升角及煤岩抗压强度进行了四因素三水平正交试验,确定了因素影响显著性及影响趋势。基于实际工况,以截割比能耗、载荷波动系数、截齿磨损深度、叶片磨损深度、生产率为目标函数建立了多目标优化函数,得到了滚筒工作参数的最佳匹配,即当牵引速度为5.66 m/min、滚筒转速为80.40 r/min、螺旋升角为11.35°时,滚筒性能达到最优。

(4)离散元数值模拟方法实现了滚筒磨损的量化分析,为滚筒磨损问题研究提供了一种新的思路,极大地降低了试验成本,对研发高效、强力、耐磨滚筒具有借鉴意义。

参考文献(References):

[1] 张磊.采煤机磨损检测方法和减摩改进方法研究[D].北京:中国地质大学(北京),2018.

ZHANG Lei.A study of detection and anti-friction improvement method for shearer wear[D].Beijing:China University of Geosciences(Beijing),2018.

[2] 刘晓辉.镐型截齿与煤岩互作用力学与磨损特性研究[D].徐州:中国矿业大学,2016.

LIU Xiaohui.Research on mechanical and wear characteristic of conical pick interacted with coal-rock[D].Xuzhou:China University of Mining and Technology,2016.

[3] 张强,顾颉颖,刘峻铭,等.基于小波包与SOM 神经网络的截齿磨损状态识别[J].煤炭学报,2018,43(7):2077-2083.

ZHANG Qiang,GU Jieying,LIU Junming,et al.Pick wear condition identification based on wavelet packet and SOM neural network[J].Journal of China Coal Society,2018,43(7):2077-2083.

[4] 张强,刘志恒,王海舰,等.基于BP神经网络的截齿磨损程度在线监测[J].中国机械工程,2017,28(9):1062-1068.

ZHANA Qiang,LIU Zhiheng,WANG Haijian,et al.On-line monitoring of pick’s wear degrees based on BP neural network[J].China Mechanical Engineering,2017,28(9):1062-1068.

[5] 赵洪盛.基于聚焦形貌恢复的采煤机截齿磨损测量分析[D].太原:太原理工大学,2015.

ZHAO Hongsheng.Measurement and analysis of shearer pick wear based on shape from focus[D].Taiyuan:Taiyuan University of Technology,2015.

[6] DEWANGAN S,CHATTOPADHYAYA S,HLOCH S.Wear assessment of conical pick used in coal cutting operation[J].Rock Mechanics & Rock Engineering,2015,48(5):2129-2139.

[7] DEWANGAN S,CHATTOPADHYAYA S.Characterization of wear mechanisms in distorted conical picks after coal cutting[J].Rock Mechanics & Rock Engineering,2016,49(1):225-242.

[8] 史秀宝,曹志娟,袁泽清.滚筒末端齿座磨损原因的模拟分析[J].矿山机械,2017,45(8):9-12.

SHI Xiubao,CAO Zhijuan,YUAN Zeqing.Simulation on wear of pick holder at drum end and causal analysis[J].Mining & Processing Equipment,2017,45(8):9-12.

[9] 卜庆辉,史秀宝,刘洁,等.采煤机螺旋滚筒失效原因分析及改进[J].矿山机械,2018,46(2):1-5.

BU Qinghui,SHI Xiubao,LIU Jie,et al.Causal analysis on failure of spiral drum of shearer and improvement[J].Mining & Processing Equipment,2018,46(2):1-5.

[10] 赵丽娟,赵名扬,刘旭南.采煤机截齿失效形式及对策探讨[J].机械强度,2017,39(5):1224-1229.

ZHAO Lijuan,ZHAO Mingyang,LIU Xunan.Study on failure mode and countermeasures of shearer picks[J].Journal of Mechanical Strength,2017,39(5):1224-1229.

[11] 赵丽娟,田震,郭晨光.矿用截齿失效形式及对策[J].金属热处理,2015,40(6):194-198.

ZHAO Lijuan,TIAN Zhen,GUO Chenguang.Failure modes and countermeasures of a mining pick[J].Heat Treatment of Metals,2015,40(6):194-198.

[12] 罗晨旭.滚筒采煤机开采含煤岩界面煤层截割特性研究[D].徐州:中国矿业大学,2015.

LUO Chenxu.Study on cutting characteristics of shearer mining coal seam with coal-rock interface[D].Xuzhou:China University of Mining and Technology,2015.

[13] 张丹,刘春生,李德根.瑞利随机分布下滚筒截割载荷重构算法与数值模拟[J].煤炭学报,2017,42(8):2164-2172.

ZHANG Dan,LIU Chunsheng,LI Degen.Reconstruction algorithm and its numerical simulation for drum load based on Rayleigh stochastic theory[J].Journal of China Coal Society,2017,42(8):2164-2172.

[14] 赵丽娟,田震,刘旭南,等.薄煤层采煤机滚筒载荷特性仿真分析[J].系统仿真学报,2015,27(12):3102-3108.

ZHAO Lijuan,TIAN Zhen,LIU Xu’nan,et al.Simulation analysis of load characteristic of thin seam shearer drum[J].Journal of System Simulation,2015,27(12):3102-3108.

[15] 高魁东.薄煤层滚筒采煤机装煤性能研究[D].徐州:中国矿业大学,2014.

GAO Kuidong.Study on coal-loading performance of thin coal seam shearer[D].Xuzhou:China University of Mining and Technology,2014.

[16] 刘春生.滚筒式采煤机理论设计基础[M].徐州:中国矿业大学出版社,2003.

[17] 梅益铭.单磨粒磨削中的磨粒磨损过程及磨损状态检测方法研究[D].杭州:浙江大学,2017.

MEI Yiming.Research on grit wear process and monitoring technique of grain wear in single-grit grinding[D].Hangzhou:Zhejiang University,2017.

[18] 张延强.WK-75型矿用挖掘机斗齿的磨损分析及结构改进[D].太原:太原理工大学,2016.

ZHANG Yanqiang.Wear analysis and structure improvement on bucket tooth of WK-75 mining excavator[D].Taiyuan:Taiyuan University of Technology,2016.

[19] 田震,赵丽娟,刘旭南,等.基于离散元法的螺旋滚筒装煤性能研究[J].煤炭学报,2017,42(10):2758-2764.

TIAN Zhen,ZHAO Lijuan,LIU Xu’nan,et al.Loading performance of spiral drum based on discrete element method[J].Journal of China Coal Society,2017,42(10):2758-2764.

[20] 赵丽娟.采煤机工作机构优化设计及载荷计算软件[P].中国专利:2014SR102903.3,2014-07-22.

ZHAO Lijuan.Optimum design of shearer working mechanism and load calculating software[P].China Patent:2014SR102903.3,2014-07-22.

[21] 周文潮.复杂煤层条件下采煤机破煤过程的动力传递问题研究[D].阜新:辽宁工程技术大学,2017.

ZHOU Wenchao.Study on the power transmission of coal cutting process of shearer under complex coal seam condition[D].Fuxin:Liaoning Technical University,2017.

[22] 郝志勇,周正啟,袁智,等.基于实验测试的采煤机截割载荷的分形分布规律研究[J].应用力学学报,2019,36(2):417-423.

HAO Zhiyong,ZHOU Zhengqi,YUAN Zhi,et al.Research on fractal distribution law of shearer cutting load Based on experimental test[J].Chinese Journal of Applied Mechanics,2019,36(2):417-423.

[23] 赵丽娟,赵宇迪,金鑫,等.夹矸煤层条件下螺旋叶片磨损失效数值模拟研究[J].煤炭学报,2019,44(11):3589-3595.

ZHAO Lijuan,ZHAO Yudi,JIN Xin,et al.Numerical simulation study on wear failure of helical blade in coal seam with dirt band[J].Journal of China Coal Society,2019,44(11):3589-3595.

[24] XIA R,WANG X,LI B,et al.Discrete element method-(DEM-) based study on the wear mechanism and wear regularity in scraper conveyor chutes[J].Mathematical Problems in Engineering,2019:4191570.

[25] 王雁翔,杨道龙,王凌翔.镐形截齿截割磨损研究[J].煤矿机械,2014,35(2):47-49.

WANG Yanxiang,YANG Daolong,WANG Lingxiang.Study on wear performance of pick cutting[J].Coal Mine Machinery,2014,35(2):47-49.

[26] 赵丽娟,赵名扬.薄煤层采煤机装煤性能研究[J].煤炭学报,2017,42(7):1892-1898.

ZHAO Lijuan,ZHAO Mingyang.Loading performance of thin seam shearer[J].Journal of China Coal Society,2017,42(7):1892-1898.

[27] 葛帅帅,秦大同,胡明辉,等.不同截割工况下变速截割采煤机截割性能多目标优化[J].煤炭学报,2018,43(8):2338-2347.

GE Shuaishuai,QIN Datong,HU Minghui,et al.Multi-objective optimization for cutting performance of the variable speed cutting shearer under different conditions[J].Journal of China Coal Society,2018,43(8):2338-2347.

[28] 赵丽娟,王雅东,刘旭南.薄煤层采煤机强力螺旋滚筒设计研究[J].机械科学与技术,2019,38(11):1712-1719.

ZHAO Lijuan,WANG Yadong,LIU Xunan.Design research for powerful screw drum of thin coal seam shearer[J].Mechanical Science and Technology for Aerospace Engineering,2019,38(11):1712-1719.

[29] 燕浩,刘梅清,梁兴,等.基于正交试验大型轴流泵空化特性的数值模拟[J].华中科技大学学报(自然科学版),2014,42(12):35-40.

YAN Hao,LIU Meiqing,LIANG Xing,et al.Numerical simulation on cavitation characteristics of large scale axial-flow pumps based on orthogonal experiment[J].Journal of Huazhong University of Science and Technology (Natural Science Edition),2014,42(12):35-40.

[30] 赵丽娟,范佳艺.基于遗传算法的采煤机螺旋滚筒多目标优化设计[J].中国机械工程,2018,29(5):591-596.

ZHAO Lijuan,FAN Jiayi.Shearer’s helical drum multi-objective optimization design based on GA[J].China Mechanical Engineering,2018,29(5):591-596.

[31] 胡春梅.采煤机滚筒辅助设计软件的开发与应用[D].阜新:辽宁工程技术大学,2014.

HU Chunmei.Development and application of shearer drum-aided design software[D].Fuxin:Liaoning Technical University,2014.

Discrete element simulation analysis on the wear characteristics of drum in coal seam with gangue

ZHAO Lijuan,JIN Xin,ZHAO Yudi,WANG Bin

(School of Mechanical Engineering,Liaoning Technical University,Fuxin 123000,China)

Abstract:During its working process,the spiral drum completes two functions of coal-cutting and coal-loading.The load on the spiral drum of a shearer working under complex coal seam conditions has nonlinearity,time-varying,impact and strong coupling.Therefore,the wear problem of drum is inevitable,which seriously affects its working performance and production efficiency.Based on the virtual prototyping technology,the discrete element numerical simulation method was used to study the wear characteristics of drum in coal seam with gangue,so as to improve its cutting performance and prolong its service life.Based on the coal-rock cutting mechanism and the discrete element method,the interaction relationship between drum and coal-rock and the wear model of drum were analyzed.Taking the MG2×55/250-BWD shearer spiral drum as an example,and the Hertz-Mindlin with bonding model was used as the contact model between particles,the Hertz-Mindlin with Archard Wear model was used as the contact model between coal-rock and drum,the coupling model was established for drum cutting in the coal seam with gangue by Pro/E and EDEM.The influence trend of the traction speed,the drum rotation speed,the spiral angle and the compressive strength of coal-rock on the drum wear was analyzed by using the orthogonal test method,and the evaluation models of pick and blade wear depth were obtained through a multiple regression analysis.Based on the actual working condition,the cutting ratio energy consumption,load fluctuation coefficient and productivity were comprehensively considered,the multi-objective optimization function was constructed and solved by MATLAB to obtain the optimal working parameters of the drum.The results show that the drum assembly has different degrees of wear,and the wear of the local contact area of the pick tip is the most serious.The main wear area of the blade is near the pick root of the tail piece.The compressive strength of coal-rock and the traction speed increase,the wear of pick and blade is intensified.The drum rotation speed increases,the pick wear depth decreases,the blade wear depth increases first and then decreases.The spiral angle increases,the pick wear depth increases first and then decreases,the blade wear depth gradually decreases.The compressive strength of coal-rock is the main factor causing the wear of pick and blade.When the compressive strength of coal-rock is 30 MPa,and the traction speed is 5.66 m/min,the drum rotation speed is 80.40 r/min,and the spiral angle is 11.35°,the drum performance is optimal.The cutting ratio energy consumption is 0.713 2 kW·h/m3,the load fluctuation coefficient is 0.081 9,the wear depth of pick,blade and end plate are 1.595 μm,1.061 7 nm and 0.907 nm,respectively,and the productivity is 315.465 5 t/h.The distribution of the drum wear in the discrete element simulation results is consistent with the actual wear of the drum under the field condition,which verifies the feasibility of the simulation model and process.The quantitative analysis of the drum wear problem can be realized using the research method,which can greatly reduce the test cost and contribute to the research and development of efficient,powerful and wear-resistant drum.

Key words:spiral drum;coal seam with gangue;wear;discrete element numerical simulation;orthogonal test

中图分类号:TD421

文献标志码:A

文章编号:0253-9993(2020)09-3341-10

移动阅读

赵丽娟,金鑫,赵宇迪,等.含夹矸煤层滚筒磨损特性离散元模拟试验分析[J].煤炭学报,2020,45(9):3341-3350.

ZHAO Lijuan,JIN Xin,ZHAO Yudi,et al.Discrete element simulation analysis on the wear characteristics of drum in coal seam with gangue[J].Journal of China Coal Society,2020,45(9):3341-3350.

收稿日期:20190624

修回日期:20190903

责任编辑:陶 赛

DOI:10.13225/j.cnki.jccs.2019.0843

基金项目:国家自然科学基金资助项目(51674134);辽宁省教育厅重点资助项目(LJ2017ZL001)

作者简介:赵丽娟(1964—),女,辽宁阜新人,教授,博士生导师,博士。E-mail:zzz2120@126.com

通讯作者:金 鑫(1991—),女,辽宁葫芦岛人,博士研究生。E-mail:jx032615@163.com