基于离散元法的粉碎机冲击力测试

邢志中 郭卫 张武刚 王渊 路正雄 李实军

摘要：为分析前后立式粉碎机粉碎内腔所受冲击力的情况及其影响因素，对改进前后立式粉碎机、中黄39大豆进行实体测量，测定其物理特性与力学特性，建立其EDEM模型并进行模拟仿真，得知：改进前后粉碎机粉碎下盖内腔所受冲击力均大于各自上盖内腔，改进后粉碎机下盖内腔所受冲击力大于改进前，改进前后粉碎机上盖内腔所受冲击力大小无明显区别。以1. 01 -1. 20 s改进前后粉碎机粉碎下盖内腔所受冲击力为评价指标，研究不同大豆数目以及不同粉碎刀具的转速对下盖内腔所受冲击力的影响。结果表明：随着大豆颗粒从100粒增至800粒，改进前后粉碎机下盖内腔所受冲击力增大值出现在大豆数目为600左右。随着粉碎刀具转速从6200增至7 000 rpm改进前后粉碎机下盖内腔所受冲击力逐渐增大。且在相同转速以及相同大豆颗数的情况下，改进后粉碎机下盖内腔所受冲击力均大于改进前。改进前后粉碎机下盖内腔所受冲击力最大值出现大豆数目为600，刀具转速为7 000 rpm，其值分别为64. 97与69. 56 N，改进前后粉碎机下盖内腔所受冲击力最小值出现大豆数目为400，刀具转速为6 200 rpm，其值分别为39. 35与41. 27 N.试验与软件仿真结果所得结论一致，该结果为研究立式粉碎机的粉碎内腔冲击力提供理论参考与依据。

关键词：离散元;中黄39大豆;立式粉碎机;冲击力;粉碎内腔

中图分类号：TD 402

文献标志码：A

文章编号：1672 -9315（2019）06 -1041 -08

DOI：10. 13800/j.cnki.xakjdxxb. 2019. 0617

开放科学（资源服务）标识码（OSID）：

收稿日期：2019 -04 -18

责任编辑：高佳

基金项目：国家重点研发计划子课题（ 2017 YFC0804310）

第一作者：邢志中（1991 -]，男，山西运城人，博士研究生，E-mail：xustxingzz@ 163.com

通信作者：郭卫（1955 -），男，陕西黄陵人，教授，博士生导师，E-mail：guow@ xust.edu.cn

0 引言

粉碎机的更新推动着各行各业的不断进步，矿业、食品等行业的发展无一不需要粉碎机[1]。而粉碎机的磨损已经是粉碎设备行业所常见的现象，同时它也是各制造商及学者非常重视的问题。粉碎内腔所受冲击力为粉碎机磨损的一个重要因素，因此，针对粉碎内腔所受冲击力展开研究，分析其影响因素，对延长粉碎机的使用寿命，避免不必要的损害，具有十分重要的意义[2-3]。

离散元法是求解与分析复杂离散系统的运动规律与力学特性的一种新型数值方法，而EDEM作为全球首个现代化离散元模型模拟的软件，已在制造业、矿山、食品行业等领域发挥了重要作用[4-5]。中国矿业大学的梁艳坤通过离散元法对神东矿区哈拉沟煤矿垮落带中破碎岩体溃砂的速度进行了深度研究，最终得知松散砂在溃砂通道内流动的速度分布呈现出纺锤形[6]。清华大学的李睿采用离散元法对球床的堆积结构开产研究，从仿真与试验验证结果得知：在球床的大空间区域颗粒堆积具有随机性，边缘部位则呈现有序堆积[7]。中国农业大学的彭飞应用EDEM软件研究喂料器的喂料过程，得到了主轴直径、螺距、物料转速与喂料稳定性的关系，并确定了这3个因素的最佳参数组合[8]。本研究应用离散元软件EDEM研究立式粉碎机粉碎内腔所受冲击力的大小与分布情况。同时由笔者已发表学术论文以及试验研究，得知改进后的圆锥上盖粉碎机粉碎效果优于圆柱上盖粉碎机（改进前粉碎机）[9]。在此基础上，分析改進前后粉碎机粉碎内腔所受冲击力的差别，大豆颗数、粉碎刀具转速对改进前后粉碎机内腔所受冲击力的影响，所得结论为今后的研究与实际生产、延长其使用寿命提供一定的参考依据。

1 研究对象及粉碎机模型的建立

1.1 立式粉碎机及其模型

研究对象为小型立式粉碎机与中黄39大豆，利用EDEM软件模拟分析大豆对小型立式粉碎机内腔的冲击力。研究所用的立式粉碎机型号为HC-700（图1），改进前后粉碎机上盖内腔如图2所示，对立式粉碎机几何尺寸进行测量，依此采用UG软件建立其三维模型。在EDEM软件仿真过程中，只需建立颗粒所接触的几何体，因此对其模型进行简化[1O -12]，改进前后粉碎机内腔三维模型如图3（a）、（b）所示。

1.2 颗粒模型

本研究以中黄39大豆为样本（图4）。随机选取800粒，测量其长、宽、高[13]。计算出每颗大豆的等效直径D与球形率φ，从而求出大豆长、宽、

将前文UG软件建立的粉碎机三维模型导入到EDEM中。根据粉碎机的工作转速，设置粉碎刀具转速为6 200 rpm.添加颗粒工厂于1×10 s内生成100颗中黄39大豆，其大小采用正态分布生成，参照表1设置中黄39大豆直径平均值为3. 452 mm.标准差为0.087.

在求解器中设置仿真时间总长4s，1-3s为粉碎机运行时间，除此之外，设置数据写出频率为0.01s，网格边长为最小大豆颗粒半径的2倍[22 -25]。

2.2 粉碎机内腔所受冲击力分析

在EDEM后处理的coloring模块中选择部件为粉碎机内腔，为便于观察，将粉碎机内所有颗粒隐藏，分析改进前后粉碎机内腔所受颗粒冲击力的大小及分布情况。图4（a）、（b）、（c）、（d）分别为1. 02及1.12 s时，改进前后粉碎机内腔所受颗粒冲击力的情况。

从图4（a）、（b）、（c）、（d）上下腔体所受冲击力云图可以看出，颜色越趋近于深灰色冲击力越小，颜色越趋近于红色冲击力越大，在粉碎机工作的过程中，改进前后粉碎机内腔所受颗粒的冲击力多集中于下腔体。通过EDEM后处理中的数据导出功能，将1.01 -1. 20 s时间段改进前后粉碎机上下盖内腔所受冲击力导出进行比较，重复9次仿真实验，求出各时间点上下盖内腔所受冲击力均值。其中A代表改进前的粉碎机，B代表改进后的粉碎机（图5）。

从图5可以看出，1. 01 -1.20 s时间段内，改进前后粉碎机下盖内腔所受冲击力分别大于各自上盖内腔所受冲击力。其中，改进前粉碎机上盖内腔的最大受力为10. 32 N，最小受力为1.17 N，受力波动较平缓;下盖内腔的最大受力为34. 06 N，最小受力为18. 69 N，受力波动较大。改进后粉碎机上盖内腔的最大受力为11. 74 N，最小受力为1.16N，受力波动较平缓，下盖内腔的最大受力为43. 98N，最小受力为25. 02 N，受力波动较大。因此在改进前后粉碎机上下盖设计的过程中，应着重考虑下盖的强度。除此之外得知，改进前后粉碎机上盖内腔所受冲击力无明显大小区别，而改进后粉碎机下盖内腔所受冲击力除在1. 12与1.17 s时小于改进前粉碎机，其余各时间点改进后粉碎机下盖所受冲击力均大于改进前。进而可以得出虽然改进后粉碎机粉碎效果有所提高，但其下盖内腔所受冲击力也明显增大。

3 大豆数量、刀具转速对粉碎内腔所受颗粒冲击力的影响

3.1 大豆数量对粉碎机内腔所受冲击力的影响

按照上诉仿真步骤，分别在粉碎机内生成100，200，300，400，500，600，700，800粒中黄39大豆进行仿真，其他仿真条件保持不变，每种大豆数目均进行9次重复试验。同样通过EDEM软件将1.01-1. 20 s时间段内，改进前后粉碎机下盖内腔所受冲击力导出，并求出每次仿真不同大豆數目在各时间点所对应冲击力数值之和的平均值，再对不同大豆数目所对应的9次重复仿真结果求均值，绘制相关折线图（图6）。其中A代表改进前的粉碎机，B代表改进后的粉碎机。

从图6可以看出，1. 01 -1.20 s时间段内，粉碎机内生成不同大豆数量进行粉碎，下盖内腔所受到冲击力均不相等。改进前后粉碎机下盖内腔所受冲击力最小值，均出现在大豆数量为100颗时，其值分别为33. 83与37. 61 N;改进前后粉碎机下盖内腔所受冲击力最大值，均出现在大豆数量为600颗时，其值分别为74. 59与83. 21 N.改进前后粉碎机下盖内腔所受冲击力与大豆颗粒数目有关，并且当大豆数量小于600时，随着大豆颗数的增加，下盖内腔所受冲击力逐渐增大。

3.2 刀具转速对粉碎机内腔所受冲击力的影响

同样按照首次仿真步骤，分别设置粉碎刀具转速为：6 200，6 400，6 600，6 800，7 000 rpm进行仿真，其他仿真条件不变，每种转速下均进行9次重复试验。将1. 01 -1. 20 s时间段内，改进前后粉碎机下盖内腔所受冲击力导出，并求出每次仿真不同转速在各时间点所对应冲击力数值之和的平均值，再对不同刀具转速所对应的9次重复仿真结果求出均值，绘制相关折线图（图7）。其中A代表改进前的粉碎机，B代表改进后的粉碎机。

从图7可以看出，在1. 01 -1.20 s时间段内，当粉碎刀具转速为6 200 rpm时，改进前后粉碎机的下盖内腔所受冲击力最小，分别为38. 36和39. 01N.当粉碎刀具转速为7 000 rpm时，改进前后粉碎机的下盖内腔所受冲击力最大，分别为47. 03和47. 13 N.随着粉碎刀具转速的增大，改进前后粉碎机下盖内腔所受冲击力均明显增大。

3.3 大豆颗数、刀具转速对粉碎机内腔所受冲击力的综合影响

分别在中黄39大豆颗数为400，600，800时，设置改进前后粉碎机粉碎刀具转速为6 200，6 400，6 600，6 800，7 000 rpm，其余条件不变，每种情况9次重复仿真，共进行270次仿真。分别求出1. 01 -1. 20 s时间段内，各时间点对应下盖所受冲击力数值之和的平均值，绘制表格如图8所示。

从图8可以看出，在相同大豆数目时，改进前后粉碎机随着粉碎刀具转速增加，下盖内腔所受冲击力同样明显增大。其中，改进前粉碎机下盖内腔所受冲击力最大值出现在粉碎刀具为7 000 rpm，大豆数目为600时，其值为：64. 97 N;最小值出现在粉碎刀具为6 200 rpm，大豆数目为400时，其值为：39. 35 N.改进后粉碎机下盖内腔所受冲击力最大值出现在粉碎刀具为7 000 rpm，大豆数目为600时，其值为：69. 56 N;最小值出现在粉碎刀具为6 200 rpm，大豆数目为400时，其值为：41. 27 N.除此之外，在相同的转速以及相同大豆颗数的情况下，改进后粉碎机下盖内腔所受冲击力均大于改进前。

4 试验验证

试验所用大豆为山东祥丰种业有限公司所生产的中黄39大豆，变频器为上海本川自动化科技有限公司所生产的XYS-AT2型变频器（图9）。采用BSF120-3AA-T常温应变片（南京聚航科技有限公司生产）布置在改进前后粉碎机上下盖内腔，并在粉碎机上下盖进行打孔处理，应变片的连接导线通过孔洞与外部ASMCl-9型电阻应变仪连接（每通道最高采样率为2 000 S/s，南京聚航科技有限公司生产），最终由计算机对所收集的数据进行处理（图10）。在中黄39大豆颗数为400.600，800时，通过变频器调至改进前后粉碎机粉碎刀具转速为6200，6 400，6 600，6 800，7 000 rpm，其余条件不变，每种情况重复9次试验。结果如图11，图12所示。

通过对比图5与图11，图8与图12.可以看出，在相同的转速以及相同大豆颗数的情况下，改进前后粉碎机上下盖内腔所受冲击力的仿真与试验值有略微偏差。但粉碎机上下盖内腔受力试验结果与仿真结果的总体改变趋势一致。

5 结论

1）改进前后粉碎机下盖内腔所受冲击力均大于各自上盖内腔所受冲击力，即下盖内腔比上盖内腔的磨损更为严重，在粉碎机设计过程中，应着重考虑下盖的强度。改进后粉碎机下盖内腔所受冲击力大于改进前粉碎机，改进后粉碎机上盖内腔所受冲击力与改进前上盖内腔所受冲击力无明显大小差别。

2）在相同粉碎机内，当中黄39大豆数目小于600时，随着中黄39大豆数目的增加，下盖内腔所受冲击力明显增大。即中黄39大豆数目为600左右时，下盖内腔的磨损随之越严重。

3）在相同粉碎机内，随着粉碎刀具转速的增加，下盖内腔所受冲击力明显增大。即粉碎刀具转速越高，下盖内腔磨损随之越严重。

4）在相同的转速以及相同大豆颗数的情况下，改进后粉碎机下盖内腔所受冲击力均大于改进前。即虽然改进后粉碎机提高了粉碎效果，但其下盖内腔磨损程度比改进前粉碎机更为严重。

5） EDEM软件仿真所得改进前后粉碎机内腔受力结果与试验验证结果基本一致，表明了采用EDEM软件分析粉碎机内腔受力的可行性。

参考文献（ References）：

[1] 陈树召.大型露天煤矿他移式破碎站半连续工艺系统优化与应用研究[J].煤炭学报，2012，37（3）：533 - 534.

CHEN Shu-zhao. Study on optimization and application ofsemi-continuous process system of other moving crushingstation in large open pit coal mine[J].Journal of CoalMine，2012 ，37 （3）：533 - 534.

[2]李潤军，单仁亮，李润圣，等，盾构机主驱动密封维修改造关键技术[J].西安科技大学学报，2014，34（5）：

579 -584.

LI Run-jun， SHAN Ren-Iiang， LI Run-sheng， et al.Keytechnologies of shield machine main drive seal repair andtransformation[J].Journal of Xi' an University of Sci-ence and Technology ，2014 ，34（5）：579 - 584.

[3]刘瑜，杜长龙，付林，等，煤块冲击破碎速度研究[J].振动与冲击，2011，30（3）：18 -21.

LIU Yu， DU Chang-long， FU Lin，et al.Study on the im-pact crushing speed of coal blocks[J].Study on the im-pact crushing speed of coal blocks[J].Vibration and im-pact，2011 ，30（3）：18 - 21.

[4] 李守巨，李德，于申.压头作用下岩石破碎过程的细观模拟[J].西安科技大学学报，2016，36（6）：769 - 774.

LI Shou-ju， LI De， YU Shen. Microscopic simulation ofrock fragmentation process under pressure head[ J].Journal of Xi' an University of Science and Technology，2016，36（6）：769 -774.

[5]孙学阳，安孝会，苗霖田，等.煤矿井工开采对上覆反向滑坡扰动的模拟研究[J].西安科技大学学报，2017 ，37 （1）：71 - 77.

SUN Xue-yang， AN Xiao-hui， MIAO Li-tian， et al.A sim-ulation study on the disturbance of overlying reverselandslide caused by coal mining[J].Journal of Xi' an U—niversitv of Science and Technology， 2017， 37（1）：71 - 77.

[6] 梁艳坤，隋旺华，朱涛，等.哈拉沟煤矿垮落带破碎岩体溃砂的离散元数值模拟研究[J].煤炭学报，2017 ，42（2）：470 - 476.

LIANC Yan-kun， SUI Wang-hua， ZHU Tao， et al.Dis-crete element numerical simulation of sand burst in frac-tured rock mass of collapse zone in Haragou Coal Mine[J]. Journal of Coal Mine ，2017 ，42（2）：470 - 476.

[7]李睿，任成，杨星团，等.基于单元体理论的球床堆积结构与传热算法研究[J].核动力工程，2016，37（3）：39 -42.

LI Rui， REN Cheng ，YANC Xing-tuan，et al.An analysisof structures and heat transfer for packed beds[J].Nu-clear Power Engineering， 2016 ，37（3）：39 - 42.

[8]彭飞，李腾飞，康宏彬，等.小型制粒机喂料器参数优化与试验[J].农业机械学报，2016，47（2）：51 -58.

PENG Fei， LI Teng-fei， KANC Hong-bin， et al.Parameteroptimization and test of feeder for small granulator[J].Journal of Agricultural Machinery， 2016， 47（ 2）： 51 - 58.

[9]邢志中，李贵荣，张海东，等.基于EDEM的小型立式粉碎机数值模拟与优化研究[J].食品工业，2017 （1）：206 - 209.

XING Zhi-zhong， LI Gui-rong， ZHANG Hai-dong， et al.Study on numerical simulation and optimization of smallvertical mill based on EDEM[J].Food IndustW， 2017（1）：206 - 209.

[10]毛君，刘歆妍，陈洪月，等.基于EDEM的采煤机滚筒工作性能的仿真研究[J].煤炭学报，2017，42（4）：1069 - 1077.

MAO Jun， LIU Xin-yan， CHEN Hong-yue，et al.Simula-tion study on working performance of shearer drum basedon EDEM[J].Journal of Coal Science， 2017， 42（ 4）：1069 - 1077.

[11]

Ouyang Y， Feng G，Read J J，et al.Estimating the ratio ofpond size to irrigated soybean land in Mississippi：a casestudy[J].Water Science and Technology： Water Supply，2016 ，16（6）：1639 - 1647.

[12] Thakur S C，Morrissey J P，Sun J， et al.Micromechanicalanalysis of cohesive granular materials using the discreteelement method with an adhesive elasto-plastic contactmodel[J].Cranular Matter ，2014， 16（3）：383 - 400.

[13] Woo S M， Uyeh D D， Sagong M S，et al.Development ofseeder for mixed planting of corn and soybeans[J].In-ternational Journal of Agricultural and Biological Engi-neering，2017 ，10（3）：95 - 101.

[14]陈海涛，李桐辉，王洪飞，等.气吸滚筒式垄上三行大豆密植排种器设计与参数优化[J].农业工程学报，2018 ，34（ 17）：16 - 24.

CHEN Hai-tao， LI Tong-hui， WANG Hong-fei，et al.De-sign and parameter optimization of three rows of soybeanseed metering device on air-suction drum ridge[J].Jour-nal of Agricultwal Engineering ，2018 ，34（ 17）： 16 - 24.

[15]常杰，王菊，刘勇，等，低共熔离子液体催化大豆油环氧化制备增塑剂[J].华南理工大学学报（自然科学版），2018 ，46（2）：93 - 101.

CHANG Jie， WANG Ju， LIU Yong， et al.Preparation ofplasticizers by epoxidation of soybean oil catalyzed bylow eutectic ionic liquids[J].Joumal of South China U一niversity of Technology（ Natural Science Edition）， 2018，46（2）：93 - 101.

[16] Constantine D A，Wang Y， Terrell E J.Effect of recipro-cation frequency on friction and wear of vibrating con-tacts lubricated with soybean-based BlOO biodiesel[J].Tribology Letters ，2013 ，50（2）：279 - 285.

[17]

Lisjak A，Liu Q，Zhao Q，et al.Numerical simulation of a_coustic emission in brittle rocks by two-dimensional fi—nite-discrete element analysis[J].Ceophysical JournalIntemational ，2013， 195（1）：423 - 443.

[18] Mukhopadhyay G，Palit P，Bhattacharvya S. Developmentof AISI A2 tool steel beater head for an impact crusher ina sinter plant[J].Metallography， Microstructure， and A-nalysis，2015 ，4（2）：114 - 121.

[19]王德福，王沫，李利橋.锤片式粉碎机粉碎玉米秸秆机理分析与参数优化[J].农业机械学报，2017 ，48（ 11）：165 - 171.

WANG De-fu， WANG Mo， LI Li-qiao. Mechanism analy-sis and parameter optimization of com straw crushing byhammer mill[ J]. Jomnal of Agricultural Machinery，2017 ，48（ 11）：165 - 171.

[20]金磊，曾亚武，叶阳，等.不规则颗粒及其集合体三维离散元建模方法的改进[J].岩土I程学报，2017 ，39（7）：1273 - 1281.

JIN Lei， ZENG Ya-wu， YE Yang， et al.Improvement of3D discrete element modeling method for irregular parti-cles and their aggregates[J].Journal of Geotechnical En-gineering ，2017 ，39（7）：1273 - 1281.

[21]王扬，吕凤妍，徐天月，等，大豆籽粒形状和尺寸分析及其建模[J].吉林大学学报（工学版），2018，48（2）：507 -517.

WANG Yang， LV Feng-yan， XU Tian-yue， et al.Shapeand size analysis and modeling of soybean seeds[J].Journal of Jilin University（ Engineering Edition）， 2018，48（2）：507 - 517.

[22]孙逸飞，宋顺翔，高玉峰.级配散粒体循环荷载下变形特性的离散元模拟[J].中国矿业大学学报，2018 ，47（4）：874 - 878.

SUN Yi-fei， SONC Shun-xiang， CAO Yu-feng. DEM sim-ulation of deformation characteristics of graded granularmaterials under cyclic loading[J].Journal of China Uni-versity of Mining and Technology， 2018， 47（ 4）： 874 -878.

[23]

Celigueta M A，Latorre S， Arrufat F， et al.Accurate mod-elling of the elastic behavior of a continuum with the Dis-crete Element Method[J].Computational Mechanics，2017 ，60（6）：1-14.

[24] Lee S H， Choi M， Kim T T，et al_Switching terahertz waves with gate-controlled active graphene metamaterials[J]. Nature Materials ，2012 ，11（ 11）：936 - 941.

[25]龍雪，刘社文，季顺迎，水位变化对正倒锥体冰载荷影响的离散元分析[J].力学学报，2019 ，51（1）：74 - 84.

LONC Xue， LIU She-wen， JI Shun-ying. Discrete elementanalysis of the impact of water level change on ice load offorward and inverted cones[J]. Joumal of Mechanics，2019 ，51（1）：74 - 84.