不同因素下掘锚机滚筒截割性能仿真分析.pdf

基：D-02)收#Jul.JOURNALOFMACHINEDESIGN20232023年月No.7Vol.40机设第4 0 卷第7 期计械不同因素下掘锚机滚筒截割性能仿真分析陈洪月1.2，海秋晗，王春刚3，张德生4（1.辽宁工程技术大学机械工程学院，辽宁阜新123000;2.辽宁工程技术大学矿产资源开发利用技术与装备研究院，辽宁阜新123000;3.陕西延长石油巴拉素煤业有限公司，陕西西榆林719000;4.天地科技股份有限公司，北京100013)摘要：为研究不同因素下掘锚机滚筒简的截割性能，分别对不同转速、煤岩硬度和振动现象下滚筒截割的动态过程进行仿真分析。通过对截割阻力、波动系数、截割比能耗、底板平整性和力在截割面后方煤壁的传导进行分析，得到不同因素对滚筒截割性能的影响。研究结果表明：滚筒受到的截割阻力随煤岩硬度的升高而增大；振动对滚筒截割阻力也产生一定影响，相较于煤岩硬度，振动对滚筒截割阻力的影响要小；波动系数与截割比能耗随着煤岩硬度的升高及振动的产生而增大。衡量底板平整性的3项指标为最大间隙值、煤岩标准差与平整度，其均随滚筒转速的提高而增加；滚筒在截割过程中，位于煤岩两侧的颗粒受到的力要比位于中间的稳定，当颗粒未受到约束时，记录点载荷由上至下逐渐减小，由前至后逐渐减小。关键词：EDEM；煤岩硬度；振动；截割阻力；底板平整性中图分类号：TD421文献标识码：A文章编号：10 0 1-2 354(2 0 2 3)0 7-0 0 4 5-0 9Simulation and analysis of cutting performance of windlass excavator sdrum under different working conditionsCHEN Hongyue-2,HAI Qiuhan,WANG Chungang,ZHANG Desheng*4(1.School of Mechanical Engineering,Liaoning Technical University,Fuxin 123000;2.Institute of Mineral Resources Development and Utilization Technology and Equipment,Liaoning Technical University,Fuxin 123000;3.Shaanxi Yanchang Petroleum Balasu Coal Industry Co.,Ltd.,Yulin 719000;4.Tian Di Science&Technology Co.,Ltd.,Beijing 100013)Abstract:In this article,in order to explore the cutting performance of the windlass excavators drum under different fac-tors,the dynamic cutting process with different rotational speeds,coal-and-rock hardness and vibration is simulated and analyzed.Through the analysis on cutting resistance,wave coefficient,cutting-rotio energy consumption,floor flatness,force conduction inthe coal wall behind the cutting face,the influence of different factors on the drums cutting performance is worked out.The re-sults show that the drums cutting resistance increases with the ever-growing coal-and-rock hardness.The vibration also has cer-tain influence on the drums cutting resistance.Compared with the coal-and-rock hardness,the vibration has less influence onthe drums cutting resistance.The wave coefficient and the cutting-ratio energy consumption increase with the ever-growing coal-and-rock hardness and vibration.The three indexes to measure the floor flatness,such the maximum clearance value,the coal-and-rock standard deviation,and the floor flatness,increase with the drums rotating speed.In the cutting process,the force on日期：2 0 2 1-0 6-2 1；修订日期：2 0 2 3-0 1-2 4金项目：国家自然基金面上项目（518 7 4 157）；辽宁省教育厅基础研究项目（LN2020JCL026）；辽宁工程技术大学创新领军团队（LNTU20T46机计设第4 0 卷第7 期械the particles located on both sides of the coal and rock is more stable than that located in the middle.When the particles are notconstrained,the load at the recording point decreases from top to bottom and from front to back.Key words:EDEM;coal-and-rock hardness;vibration;cutting resistance;floor flatness掘锚机作为一种新型采掘设备，可同时进行掘进、锚杆支护和运煤等作业，其工作性能与滚筒有直接的联系。滚筒作为掘锚机主要截割部件之一，其截割性能直接决定了掘锚机的工作效率。不同于采煤机滚筒，掘锚机滚筒没有螺旋叶片且齿座是直接焊在滚筒上的，这使掘锚机滚筒与采煤机滚筒有着很大的区别。国内外对掘锚机滚筒的研究有：麻成标基于单镐齿和采掘头的切削机制，采用ANSYS软件进行采掘头破岩过程的仿真模拟，得出了采掘头在转速和布齿方式双变量条件下对切削性能的影响规律；刘在政2 采用LS-DYNA软件对掘锚机滚筒截齿破煤过程进行了仿真研究，分析破煤过程中截齿的应力状态并引用截割阻力的计算公式对仿真进行验证，证明了建模方法的可行性与正确性，研究工作为深人探讨破煤机理，并为改进截齿的设计提供可靠的依据；魏晓华等3 通过MATLAB软件进行仿真，得到了工作状态下掘进机截割头的转速和转矩的变化，为掘进机优化设计提供了理论依据等。文中以MB670-1型号掘锚机为研究对象，通过Pro/E软件建立滚筒简化模型，采用EDEM仿真软件建立煤岩模型，模拟掘锚机滚筒的截割过程。综合考虑煤岩硬度、滚筒转速和振动现象等因素对掘锚机滚筒截割性能的影响，采用单因素试验法对截割阻力、波动系数、截割比能耗、底板平整性和煤岩压缩力进行分析总结，为研究不同因素滚筒截割性能的影响提供了新的方向。1滚筒截割过程仿真模型的建立1.1颗粒模型的建立掘锚机滚筒截落煤块的尺寸通常在0 50 mm之间4 。在EDEM软件中定义6 种不同尺寸的颗粒来模拟煤岩，如图1所示其中，图1a为岩石颗粒，图1b图lf为煤颗粒。在EDEM软件中添加煤、岩石和钢3种材料。各参数如表1和表2 所示(a)1 mm(b)3mm(c)5mm(d)10 mm(e）30mm(f)50mm图1颗粒模型表1材料参数材料泊松比剪切模量G/Pa密度p/(kg/m)J4煤0.159.561081 420J5煤0.171.62101420J6煤0.221.721091 420岩石0.251072600钢0.3710107 800表2接触参数接触参数恢复系数静摩擦因数动摩擦因数煤一煤0.50.60.05煤一岩石0.50.70.08煤一钢0.50.40.05岩石-岩石0.60.80.1岩石-钢0.60.50.07钢一钢0.70.20.01为更好地模拟煤岩的物理状态，选择Hertz-Mind-linBonding模型。如图2 所示，模型通过“黏结键 将颗粒黏结。模型在相互接触的颗粒间传递法向应力F、切向应力力F,和力矩M,并且承受一定的切向运动和法向运动5-8 M2图2颗粒黏结模型.m:472023年7 月陈洪月，等不同因素下掘锚机滚筒截割性能仿真分析通过设置颗粒间的法向和切向黏结强度使颗粒能够承受一定的拉伸与剪切作用，当颗粒间的作用力超过黏结强度时，黏结会发生破坏。此后，颗粒将视作刚性球体进行接触求解。定义法向应力与切向应力的最大值如下9-10 F2T.n+RmaxAJ(1)-FTTRminAJ式中：R-颗粒黏结半径，mm,R=VR,R,；R,R2-球1,2 的颗粒半径,mm；F.,F,颗粒的法向、切向力，N；T.,T.果颗粒的法向、切向力矩，Nm；A单位接触面积,mm，A=R;1J一一极惯性矩,mmt,JTR42颗粒间的黏结力和力矩随时步变化的数学表达式如下：8F=-U.S.AStSF,=-0,S,Aotn=-W,S.J8t(2)ST,=St2式中:8 F.,8F,相应时间点颗粒的法向黏结力、切向黏结力，N；ST,8T,相应时间点颗粒间的法向力矩、切向力矩，St-时间步长，s；U,颗粒的法向、切向速度，m/s；颗粒的法向、切向角速度，rad/s；S.,S.颗粒的法向、切向接触刚度，N/m。通过计算确定表1中不同硬度煤岩颗粒间的黏结参数，如表3所示表3颗粒黏结参数数值颗粒黏结参数J4煤/5煤./6煤单位面积法向刚度/（N/m）4.01096.91097.6109单位面积切向刚度/（N/m）1.71092.91093.1109单位面积法向应力/Pa4.01065.01066.0106单位面积切向应力/Pa1.61062.01062.41061.2仿真模型的建立参考山特维克MB670-1型号掘锚机滚筒，确定主要参数,并采用Pro/E软件建立滚筒的三维模型。滚筒主要参数如表4 所示。表4滚筒主要参数滚筒宽度B/mm滚筒直径D,/mm掘进速度u/(m/min)52001 1503在EDEM软件中建立一个5.5m2.0m3.0m的箱体用于填充煤岩颗粒，使填充后的颗粒生成黏结键并将顶板接触面和滚筒截割面进行隐藏。在煤层上方建立一个箱体并通过MCU工具给予一定的力来模拟煤岩受到的顶板压力，单位面积顶板压力计算公式为11-14 DLLok,P。=z h k z +-(3)2Lk式中：P一-顶板压力,N;hi一动压系数;k2一悬顶片帮系数；-直接顶、基本顶岩石密度，kg/m；L控顶距,m;Lo基本顶初次来压步距，m；h,D-直接顶、基本顶岩层厚度，m通过计算得出顶板压力值为6.6 10 7 N。滚筒截割模型如图3所示，X向为滚筒截割方向，Y向为竖直方向，Z向为滚筒轴向方向。同时，顶板上方受到一定的顶板压力，生成的煤岩两侧受到箱体的固定约束。顶板顶板压力掘锚机滚筒固定固定约束约束巷道面V煤壁图3滚筒截割模型1.3设定仿真参数设置时间储存间隔为0.0 5s，网格尺寸为2 R，mino48计机设第4 0 卷第7 期械2滚筒截割性能的衡量指标2.1波动系数波动系数指的是滚筒受到的载荷的波动大小，其表达式为：(4)1=1式中：F一滚筒瞬时载荷，N；F滚筒载荷均值，N，F=F滚筒载荷记录时间点总数该公式用于计算表5与表6 中的波动系数，具体计算结果如表5与表6 所示表5不同硬度系数下滚筒截割的波动系数与截割比能耗载荷转矩落煤波动硬度截割比能耗方向均值均值质量系数系数Hw/(kW.h/m)F/NTM/(Nm)Mm/kg8X-2.21050.41F4Y1.31051.810573080.620.61Z-7.61036.19X-2.31050.47Y1.41052.01057 4720.740.66Z-8.51037.51X-2.81050.53J6Y1.61052.41057.4120.770.80Z1.710*12.05表6不同振动下滚筒仿真结果载荷转矩均值落煤 波动截割比振动 方向均值TM/质量系数能耗H/F/N(Nm)Mm/kg8(kWh/m)mX-2.21050.41无振动Y1.31051.8x1057.3080.620.61Z-7.61036.19X-2.31050.50轴向Y1.41051.91057.3370.620.64振动Z4.310316.69X-2.31050.46竖直Y1.41051.91057.3410.650.64振动Z-5.910311.822.2截割比能耗截割比能耗指的是滚筒截割单位体积的煤岩时所消耗的能量大小，表达式如下：PmtnTMH(5)W95503600Mm式中：H一截割比能耗,kWh/m；Pm煤岩密度，kg/m；截割时间，s;n滚筒转速，r/min;TM-滚筒转矩均值，Nm；Mm截落煤岩质量，kg。式（5)用于计算表5与表6 中的截割比能耗,具体计算结果如表5与表6 所示3煤岩硬度对滚筒截割性能的影响为研究不同煤岩硬度对掘锚机滚筒截割性能的影响，取滚筒转速为30 r/min，掘进速度为3m/min，只改变煤岩硬度这一种因素，分别对f4,J5,f6这3种硬度系数的煤岩进行截割仿真。通过对仿真结果进行分析，绘制不同硬度系数下滚筒受到的截割阻力，如图4所示。6X105-f4f5f6N/（弹205101520时间/s图4不同硬度系数下滚筒所受截割阻力通过EDEM后处理模块，导出滚筒三向载荷与转矩数据,求出不同硬度系数下滚筒截割的波动系数与截割比能耗，如表5所示，结合图4 和表5可知，滚筒截割阻力、波动系数与截割比能耗的大小与煤岩硬度成正比，随着煤岩硬度增大而增加。滚筒在截割煤岩过程中，截割阻力整体是先增大后逐渐稳定的趋势，但由于煤岩构成是复杂492023年7 月陈洪月，等同因素下掘锚机滚筒截割性能仿真分析和不规律的，有时会骤增然后恢复到稳定大小。当煤岩硬度系数为f4时，滚筒受到的截割阻力最大值约为400kN,波动系数与截割比能耗较小,故滚筒载荷的波动与能量的消耗较小；当煤岩硬度系数为6 时，滚筒受到的截割阻力最大值约为6 0 0 kN,波动系数与截割比能耗较大，故滚筒载荷的波动与能量的消耗较大。为了让掘锚机具有较好的截割性能，滚筒可以选择在硬度较低的工作面上进行工作。4振动现象对滚筒截割性能的影响掘锚机在工作过程中由于受到载荷冲击，使滚筒产生振动现象【15，为研究振动现象对掘锚机滚筒截割性能的影响,通过EDEM与RecurDyn软件耦合，对掘锚机滚筒截割过程进行联合仿真。由于RecurDyn软件中无法直接导出滚筒振动曲线，故采用间接观测法确定滚筒振动频率和振幅。在RecurDyn后处理界面导出滚筒的位移曲线与加速度曲线，确定滚筒振幅和频率。将负载状态下滚筒位移曲线与空载状态下滚筒位移曲线相减，将其最大值作为仿真振幅值。图5所示为滚筒振幅曲线。10最大值5-1005101520时间/s图5滚筒振幅曲线由于加速度曲线无法直接得到滚筒的振动频率，因此，对加速度曲线进行傅里叶变换，并绘制振幅-频率曲线。将振幅最大值所对应的频率作为仿真频率。振幅-频率曲线如图6 所示。通过对图5和图6 分析可知，滚筒振幅为7.5mm，振动频率为0.6 3Hz。因此,在EDEM前处理界面设置滚筒转速为30 r/min，掘进速度为3m/min，只改变滚筒振动方向这一因素，分别对无振动、轴向振动和竖直振动3种不同振动情况进行仿真分析。不同振动下滚筒受到的截割阻力，如图7 所示1510(0.63,7.5)50-5-10-150.00.51.01.52.0频率/Hz图6滚筒振幅频率曲线6105无振动轴向振动竖直振动4205101520时间/s图7不同振动下滚筒所受截割阻力通过EDEM后处理模块，导出滚筒三向载荷与转矩数据，求出不同振动情况下滚筒截割的波动系数与截割比能耗。通过对图7 和表6 分析可知，振动现象使滚筒受到的截割阻力、波动系数与截割比能耗增大。无振动时，滚筒受到的截割阻力最大值约为4 0 0 kN，波动系数与截割比能耗较小；当受到振动现象影响后，滚筒截割阻力最大值约为50 0 kN，波动系数与截割比能耗较大。5滚筒转速对底板平整性的影响掘锚机在工作过程中随着滚筒不断截割，截割槽内底板会不断地发生变化。截割工作结束后，截割槽内底板会变得参差不齐，槽内底板的具体情况即为50机计设械第4 0 卷第7 期底板平整性。若底板平整性过差，会严重影响掘锚机后续的工作，而滚筒转速即为底板平整性的首要影响条件。为研究滚筒转速对底板平整性的影响，选取掘进速度为3m/min，煤岩硬度系数为f4作为固定参数，只改变滚筒转速这一因素。分别取转速为20,30，4 0,50 r/m i n 进行仿真，通过计算不同转速下的底板最大差值h与煤岩标准差对底板平整性进行分析5.1底板最大间隙值h在截割槽处建立尺寸为0.4 m0.2m5.0m的统计区I，统计区I如图8 所示。为了方便观察底板最大间隙值，将统计区1设定为单独显示，并采用不同颜色表示颗粒的不同高度。将不同转速时的底板正视图汇总在一起并标记出最大间隙值，如图9所示图8统计区1最高点Position/mm100mm1.55e-03最低点20r/min1.51e-03最高点M105mm1.47e-03最低点30r/min最高点1.43e-03.117mm最低点1.39e-0340 r/min最高点1.35e-03t120mm最低点50 r/min图9底板最大间隙值如图9所示，底板最大间隙值随滚筒转速的提高而增大，当转速从2 0 r/min提高至50 r/min时，底板最大间隙值从10 0 mm增大至12 0 mm，且转速为30r/min及4 0 r/min时，底板最大间隙值相差较大。由图9还可以看出，转速为2 0,30 r/min时，底板最大间隙值相近；转速为4 0,50 r/min时底板最大间隙值相近。因此，为保证底板具有较小的最大间隙值，可选择在较低的转速时进行工作5.2底板煤岩质量标准差底板煤岩质量标准差能客观准确反映统计区内煤岩质量的离散程度，将统计区等分成一定数量，求出对应的标准差，反映出底板平整性。不同情况的底板平整性与煤岩质量标准差的关系如图10 所示4kg4 kg4 kg5kgkgkg4kg4kg4kg433kg0-00=0.89（a）底板平整（b）底板较平整7kgkgkg4kg2kg0-1.67（c）底板不平整图10平整性与标准差关系示意图将图8 所示统计区1沿Z轴划分成50 份并对各网格中煤岩的质量进行统计，通过公式计算出标准差并绘制曲线图，计算式如下所示：(m/-Em)2(6)N=1式中：N-网格数；第k个网格中煤岩质量，kg；E.m-N个网格中煤岩质量的平均值，kg。不同转速下煤岩质量标准差统计图如图11所示1.21.11.00.920304050滚筒转速/（r/min）图11不同滚筒转速时的煤岩质量标准差曲线由图11可知，统计区的煤岩质量标准差随滚筒转速的提高而增大；转速为2 0,30 和4 0 r/min时，统计区内煤岩质量标准差有较大的差异，而转速为4 0 r/min512023年7 月陈洪月，等：不同因素下掘锚机滚筒截割性能仿真分析及50 r/min的煤岩质量标准差相差较小。为保证底板煤岩质量具有较小的标准差，可选择较低的滚筒转速。5.3平整度P平整度可直接反映底板的平整性，其值越小表示底板的平整性越好。具体计算式如下所示16 P=E(7)将计算后的具体数值统计在一起绘制出平整度曲线图，如图12 所示2520%/辉木151020304050滚筒转速/（T/min）图12平整度随滚筒转速变化曲线由图12 可知，底板的平整度与滚筒转速成正比，在滚筒转速较低时底板具有较好的平整性。通过对底板最大间隙值、底板煤岩质量标准差和平整度的分析可知，这3个指标均与滚筒转速成正比。因此，可选择较低的滚筒转速，以保证截割槽底板拥有较好的平整性。虽然转速为2 0 r/min时底板平整性最佳，但考虑到转速对工作效率的影响可适当调高滚筒的转速。因此，可选转速为30 r/min，既可保证滚筒具有较高的工作效率，也可使底板具有较好的平整性。6滚筒截割过程中力在煤岩中的传导滚筒在截割过程中产生的力会传向综掘面后方的煤岩，随着滚筒的不断截割，传导的力也不断地发生变化。为研究滚筒截割过程中，力在后方煤岩传递的情况，以转速为30 r/min、煤岩硬度系数为f4和无振动的工况为例建立记录区进行分析。由于滚筒掘进深度为1.0m，因此，只对滚筒未截割的区域进行分析，0 s时颗粒应力云图如图13所示。Total Force/N3.00e+022.40e+02滚简截制区城载1.80e+02荷1.20e+026.00e+010.00e+00图130s时颗粒应力云图分别在煤壁左侧、中间、右侧创建记录点，在每侧的顶部、中部和底部由前至后建立9个尺寸为0.1m0.1mx0.1m的网格作为记录点，如图14 所示，顶部记录点载荷记录区域3煤壁右侧煤壁舒456中部记录点左侧底部记录点滚筒截割区域9后前图14载荷记录点示意图在EDEM后处理界面分别导出滚筒截割0,5,10，15,20s时各记录点颗粒的载荷数据，由于记录点中包含的颗粒数量不等，故选取平均力作为单颗粒受力的参考数据。将数据导出后分别绘制每侧记录点颗粒载荷随时间变化的点线图，如图15所示。200顶部记录点：3中部记录点：46180底部记录点：8916014012010080604020005101520时间/s(a）左侧52机计第4 0 卷第7 期设械80顶部记录点：3中部记录点：一6底部记录点：897060504005101520时间/s(b)中间顶部记录点：3100中部记录点：46底部记录点：789806040420005101520时间/s（c）右侧图15不同时间点颗粒载荷变化曲线图由图15可看出，中间各记录点数据波动较大，左右两侧各记录点数据较为稳定，这主要是因为滚筒长度比仿真生成的煤壁长度小，中间各记录点所在平面处于滚筒中心面，而左右两侧各记录点所在平面在滚筒所在平面之外，这使得中间颗粒比两侧颗粒受到的载荷更为复杂。通过对图15a和图15c两图分析可知：由上到下，位于顶部的颗粒受到的载荷最大,其次是中间颗粒,最后是底部颗粒。这是因为顶部颗粒同时受到顶板的压力与滚筒截割向后方煤岩传递的力，中间颗粒主要受到滚筒截割向后方煤岩传递的力，底部颗粒虽然也受到滚筒截割向后方煤岩传递的力，但由于位置距离滚筒较远，颗粒受到的力较小；由前至后，位于后方的颗粒载荷最大，位于前端与中间的颗粒受到的载荷较小且相差不大。这是因为处于前端与中间部分的颗粒在受力后会产生细微的移动从而将力向其他方向传递，而后方颗粒由于同时受到箱体后平面与左右平面的约束无法产生移动,因此,受到的力相对较大。通过对图15b分析可知：由于中间记录区位于滚筒中心的平面，故颗粒受到的载荷相比于两侧较为复杂,但仍存在一些规律。由上至下,与两侧记录点情况相同，由于顶板压力与滚筒距离的影响，顶部的颗粒受到的载荷最大，其次是中间颗粒，最后是底部颗粒；由前至后，位于前端的颗粒受力最大,其次是中间部分的颗粒,最后是位于后方的颗粒。这是因为，前端的颗粒与滚筒距离较近，后方的颗粒与滚筒距离较远，虽然后方颗粒受到箱体后平面的约束，但由于两侧为可活动颗粒，所以受到的力会向其他方向传递。7结论以MB670-1型号掘锚机研究对象，通过EDEM软件分别对滚筒转速、振动和煤岩硬度进行仿真，分析截割阻力、波动系数、截割比能耗和底板平整性等因素的变化趋势。通过分析得出以下结论：(1)滚筒在截割煤岩过程中，截割阻力呈先增大后逐渐稳定的趋势。煤岩硬度对截割阻力有较大的影响，煤岩硬度系数从4 增大至f6时，最大截割阻力从400kN增大至6 0 0 kN；振动相比于煤岩硬度对截割阻力的影响要小，无振动时，滚筒受到的最大截割阻力约为4 0 0 kN；在振动的影响下，滚筒受到的截割阻力增大至50 0 kN；波动系数和截割比能耗随着煤岩硬度的升高和振动的产生而增大。（2)底板平整性的3项指标：最大底板间隙值、底板煤岩质量标准差和平整度均随滚筒转速的升高而增加。考虑到滚筒转速对工作效率的影响，可选取滚筒转速为30 r/min，使滚筒具有较高工作效率的同时，底板的平整性也较好（3）滚筒截割过程中，由于位置原因，两侧颗粒受到的力要比中间的稳定。当颗粒未受到约束时，顶部颗粒受到的力相较于底部颗粒受到的力要小，前端颗粒受到的力相较于后端颗粒受到的力要小。参考文献1麻成标.掘锚机采掘头镐齿布置方式优化设计J.矿业研究与开发，2 0 2 0,4 0(1）：130-134.2刘在政.掘锚机单截齿破煤特性研究J.铁道建筑技术，2016(10):90-93.532023年7 月陈洪月，等不同因素下掘锚机滚筒截割性能仿真分析3魏晓华，师建国.纵轴式掘进机截割过程动力学建模与仿真J.辽宁工程技术大学学报，2 0 13，（1）：8 9-92.4张建军，郑杨，崔正龙.刨煤机刨削试验用模拟煤壁配比及结构层研究J.煤矿机电，2 0 16（1)：16-18.5王卫华，李夕兵.离散元法及其在岩土工程中的应用综述J.岩土工程技术,2 0 0 5,19(4)：17 7-18 1.6Su O,Ali Akcin.Numerical simulation of rock cutting usingthe discrete element method J.International Journal of RockMechanics and Mining Sciences,2011,48(3):434-442.7徐泳，孙其诚，张凌，等.颗粒离散元法研究进展J.力学进展,2 0 0 3,33(2);2 51-2 6 0.8徐宝鑫.截割头截齿安装参数的离散元仿真分析D.沈阳：沈阳理工大学，2 0 15.9LI Zhanfu.A study of particles penetration in sieving processon a linear vibrationscreen J.International Journal of CoalScience&Technology,2015,2(4):299-305.10】王国强,郝万军，王继新.离散单元法及其在EDEM上的实践M.西安：西北工业大学出版社，2 0 10.11孙广义,潘启新，黄占龙.单体支柱工作面顶板压力计算的探讨J.黑龙江矿业学院学报，1998,8（1）：6-9.12孙广义,马云东，李东.顶板压力计算的盲数辨识研究J.辽宁工程技术大学学报,2 0 0 3,2 2（S1）：7-8.13 Wang Jinhua.Development and prospect on fully mechanizedmining in Chinese cial mines J.International Jourmal ofCoal Science&Technology.2014,1(3):153-160.14 Wang Jinhua,Yu Bin,Kang Hongpu,et al.Key technologiesand equipment for a fully mechanized top-coal caving opera-tion with a large mining height at ultra-thick coal seams J.International Journal of Coal Sciene&Technology,2015,2(2):97-101.15谢苗，闫江龙，毛君.采煤机截割部振动特性分析J.机械强度,2 0 17,39(2）：2 54-2 6 0.16袁建明，郭城，叶方平.新型卸料斗及落料平整堆积仿真分析J.机械设计，2 0 2 0,37（7）：4 5-50.作者简介：陈洪月（198 2 一），男，教授，博士，研究方向：机械设计及理论。E-mail:海秋哈（1997 一），男，蒙古族，硕士研究生，研究方向：机械设计及理论。E-mail: