基于EDEM的小型立式粉碎機粉碎刀具優化與研究

邢志中，張海東，王 孟，陳 騰，郭小軍，翟超男

(云南農業大學機電工程學院，云南昆明 650201)

?

基于EDEM的小型立式粉碎機粉碎刀具優化與研究

邢志中，張海東*，王 孟，陳 騰，郭小軍，翟超男

(云南農業大學機電工程學院，云南昆明 650201)

摘要[目的]提高小型立式粉碎機的粉碎性能。[方法]對小型立式粉碎機、大豆進行實體測量，建立其EDEM模型，提出粉碎刀具的優化方案，以仿真時間1.01～1.10 s內顆粒的旋轉動能均值、碰撞次數總和作為評價指標，對比優化前后的粉碎刀具性能，并重復3次仿真試驗，驗證優化的正確性。[結果]優化后粉碎機內顆粒的旋轉動能均值、碰撞次數總和均大于原始粉碎機內顆粒，其顆粒旋轉動能均值最大值為0.007 7 J，碰撞次數總和最大值為738次。[結論]優化后粉碎刀具性能明顯優于優化前粉碎刀具性能，該研究可為其他粉碎機粉碎刀具的改進提供理論參考。

關鍵詞EDEM；粉碎機；刀具；優化

隨著現代食品工業的不斷發展，人們對食品細化程度的要求越來越高[1-3]。例如，一些動植物體的不可食部分(如骨)、殼(如蛋殼)、蝦皮等可通過粉碎技術而成為易被人體吸收、利用的鈣源和甲殼素；巧克力配料粒度的細化使巧克力具有細膩滑潤的良好口感；還有調味品加工中微粉食品的巨大孔隙率造成集合孔腔，可吸收香氣經久不散[4]。提高粉碎效率、降低能耗、不斷改進超細粉碎設備對促進食品粉碎行業的蓬勃發展具有十分重要的意義。

EDEM作為全球首個多用途離散元素法建模軟件，已在礦山、化工、散粒體輸送、農業等領域發揮了重要作用[5-6]。李洪昌等[7]利用EDEM軟件仿真振動篩的工作過程，結果表明隨著振動篩振幅以及頻率的增大，所篩選的物料沿篩面后移的速度也隨之增大，并且在振幅為40 mm、頻率為6 Hz時出現了篩分損失。楊公波等[8]應用EDEM軟件對不同斗輪轉速下的斗輪機進行了數值仿真，結果表明隨著斗輪轉速的增加，挖掘阻力也隨之增加，并且在重力式卸料轉速范圍內，轉速為5 r/min卸料的效果最好。賈富國等[9]采用離散元軟件EDEM分析了糙米勻料盤的結構參數以及運行條件對料層厚度均勻性的影響，結果表明進料量對糙米的分布影響顯著，而勻料盤對其影響不顯著，并且勻料盤的最佳轉速為183 r/min。目前，有關EDEM軟件在立式粉碎機優化中的應用研究尚未見報道，鑒于此，筆者將小型立式粉碎機內的刀具進行優化，采用EDEM軟件對刀具優化前后粉碎機內的顆粒旋轉動能均值、碰撞次數總和進行比較，探討優化刀具對提高粉碎機性能的可行性，以期為立式粉碎機的優化提供一種新思路，并為今后的實際生產提供一定的參考依據。

1材料與方法

1.1立式粉碎機模型對小型立式粉碎機(圖1)進行測量，獲得其幾何尺寸，對其進行建模。由于EDEM軟件仿真當中只需要建立顆粒體與幾何體接觸的部分來進行仿真計算，并且為減輕計算機運算量和縮短計算時間，因此將立式粉碎機模型進行簡化。在UG軟件中建立其三維模型，另存為prt格式，用于后續仿真將其三維模型導入到EDEM軟件當中。其內部結構的模型截圖如圖2所示。

圖1　小型立式粉碎機Fig.1　Small vertical mill

圖2　粉碎機內部模型截圖 Fig.2　Internal screenshot of grinder model

1.2顆粒模型鑒于大豆顆粒的離散特性，該研究選用大豆品種“中黃39”為研究對象。采用精度為0.02 mm的游標卡尺測量其長、寬、高。共測量800顆，并根據公式計算每顆大豆的等效直徑與球形率。計算公式分別為：

(1)

φ=D/a×100%

(2)

式中，D為等效直徑，單位為mm；φ為球形率，單位為%；a、b、c分別為大豆長度、寬度、高度，單位均為mm。

大豆群體的三維尺寸、等效直徑及球形率結果見表1。由表1可以看出，大豆的球形率均值已接近98%，故采用EDEM軟件仿真分析小型立式粉碎機時，可將大豆顆粒看作球體。

表1　大豆群體的三維尺寸、等效直徑及球形率

1.3軟件仿真與粉碎刀具的優化設置顆粒與幾何體、顆粒與顆粒之間的接觸模型為Hertz-Mindlin(no slip)built-in。將材料庫導入到EDEM軟件中，根據材料庫及文獻[10-11]設置全局變量參數(表2)。

表2　全局變量參數

注：大豆密度采用排水法進行測定。

Note：Soybean density was detected by drainage method.

在顆粒模塊中，將大豆設置為球體，之后將UG軟件建立的粉碎機三維模型導入到EDEM中。添加顆粒工廠，命令其為虛擬，于1×10-12s內生成100顆中等大小的大豆，其大小采用正態分布生成，參照表1設置等效直徑平均值為3.442 5 mm，標準差為0.097 mm。在求解器當中設置固定時間步長百分比為20%，仿真時間總長4.00 s。圖3為仿真時間0.70 s時大豆顆粒在粉碎機內的狀態。由圖3可以看出，顆粒在粉碎機運行前，已經基本靜止穩定，因此，粉碎機于1.00 s開始工作，1.00～3.00 s為刀具粉碎顆粒時間，除此之外，設置數據寫出頻率為0.01 s，網格邊長為最小大豆半徑的2倍。

圖3　0.70 s粉碎機內顆粒狀態Fig.3　Particle state in shredder at 0.70 s

大豆在粉碎機內主要依靠刀具、粉碎內腔以及顆粒間的撞擊、摩擦、剪切等方式達到粉碎的效果，為加強刀具的粉碎性能，對粉碎機內的刀具進行優化，優化前后的刀具三維模型如圖4、5所示。

圖4　原始刀具模型Fig.4　The original tool model

圖5　優化后星形刀具模型Fig.5　The improved star model of tools

1.4仿真試驗在粉碎機內，顆粒的碰撞包括顆粒與刀具、顆粒與內壁以及顆粒之間的相互碰撞，其碰撞次數越多，粉碎效果越好。利用離散元軟件EDEM后處理器的數據導出功能對刀具優化前后粉碎機內顆粒在仿真時間為1.00～1.10 s的旋轉動能均值及顆粒碰撞次數總和進行對比，繪制相關圖表。

1.5仿真試驗驗證依照首次軟件仿真流程，所有參數設置不變，分別在優化前后粉碎機內生成100顆中等大小大豆，對比1.01～1.10 s時間段內優化前后粉碎機顆粒的旋轉動能均值以及碰撞次數總和。其中旋轉動能均值的計算方法為：對每個時間點所對應的顆粒旋轉動能均值的數值求和，再除以時間點個數。重復3次，輸出相應數據繪制圖表。2結果與分析

2.1仿真試驗結果與分析由圖6可以看出，在1.04 s這一時間點，星形刀具粉碎機內顆粒旋轉動能均值略大于原始粉碎機內顆粒的旋轉動能均值。其余各時間點，星形刀具粉碎機內顆粒旋轉動能均值均明顯大于原始粉碎機內顆粒的旋轉動能均值。

圖6　優化前后粉碎機內顆粒旋轉動能均值Fig.6　Mean of particle rotational kinetic energy in shredder before and after improvement

由圖7看以看出，在1.04與1.07 s 這2個時間點，星形刀具粉碎機內顆粒的碰撞次數小于原始粉碎機內顆粒碰撞次數，1.05 s這一時間點，兩者相等。其余各時間點，星形刀具粉碎機內顆粒碰撞次數均大于原始粉碎機內顆粒碰撞次數。對各時間點所對應的碰撞次數求和，結果表明，在1.00～1.10 s時間段內，星形刀具粉碎機內顆粒碰撞次數總和為635，在原始粉碎機內其值為541。可見，優化刀具后的粉碎機的顆粒旋轉動能均值及碰撞次數總和均明顯優于原始粉碎機。

圖7　優化前后粉碎機內顆粒碰撞次數Fig.7　Particle collision frequency in shredder before and after improvement

2.2仿真試驗驗證結果與分析由圖8可以看出，3次重復的星形刀具粉碎機內顆粒的旋轉動能均值均遠大于原始粉碎機內顆粒的旋轉動能均值。顆粒的旋轉動能均值最大值、最小值均出現在第1次重復，分別為0.007 7、0.004 0 J。

圖8　粉碎機內顆粒旋轉動能均值Fig.8　Mean of rotational kinetic energy in shredder

由圖9可以看出，3次重復的星形刀具粉碎機內顆粒的碰撞次數總和均遠大于原始粉碎機內顆粒碰撞次數總和。顆粒碰撞次數總和的最大值、最小值均出現在第1次重復，分別為738、560次。

圖9　粉碎機內顆粒碰撞次數總和Fig.9　The sum of total particle collision frequency in the mill

3結論

該研究通過利用離散元軟件EDEM對刀具優化前后的粉碎機內顆粒的旋轉動能均值及碰撞次數總和進行分析，得出星形粉碎刀具的粉碎性能優于原始刀具的結果；進行了3次重復仿真試驗來驗證優化結果，結果表明：優化刀具后粉碎機內顆粒的旋轉動能均值、碰撞次數總和均大于原始粉碎機內顆粒的旋轉動能均值、碰撞次數總和，其結果與首次仿真結果一致。旋轉動能均值最大值為0.007 7 J，碰撞次數總和最大值為738次。

參考文獻

[1] 譚向勇.中國食品工業的現狀及發展趨勢研究[J].北京工商大學學報(自然科學版)，2010，28(1)：1-7.

[2] 張艷，謝武.我國食品工業發展分析[J].市場論壇，2010(5)：40-41.

[3] 魏益民.我國農產品加工業發展現狀與趨勢分析[J].中國食物與營養，2006(10)：4-5.

[4] 高堯來，溫其標.超微粉體的制備及其在食品中的應用前景[J].食品科學，2002，23(5)：157-160.

[5] SPOGIS N.Multiphase modeling using EDEM-CFD coupling for fluent[Z].CFD OIL，2008：18-19.

[6] LI J，WEBB C，PANDIELLA S S，et al.Discrete particle motion on sieves-a numerical study using the DEM simulation[J].Powder technology，2003,133：190-202.

[7] 李洪昌，李耀明，唐忠，等.基于 EDEM 的振動篩分數值模擬與分析[J].農業工程學報，2011，27(5)：117-121.

[8] 楊公波，李郁，陳定方.基于 EDEM 仿真的斗輪堆取料機取料機理研究[J].武漢理工大學學報(交通科學與工程版)，2014，38(3)：680-683.

[9] 賈富國，姚麗娜，韓燕龍，等.基于離散元法的糙米勻料盤仿真優化設計[J].農業工程學報，2016，32(4)：235-241.

[10] VAN LIEDEKERKE P，TIJSKENS E，DINTWA E，et al.DEM simulations of the particle flow on a centrifugal fertilizer spreader[J].Powder technology，2009，190(3)：348-360.

[11] TAVAREZ F A，PLESHA M E.Discrete element method for modelling solid and particulate materials[J].International journal for numerical methods in engineering，2007，70(4)：379-404.

基金項目云南農業大學博士科研啟動基金；云南農業大學學生科技創新創業行動(2016ZKX121)。

作者簡介邢志中(1991- )，男，山西運城人，碩士研究生，研究方向：農業機械及設施。*通訊作者，副教授，博士，從事農業機械及其自動化研究。

收稿日期2016-04-18

中圖分類號S 22

文獻標識碼A

文章編號0517-6611(2016)15-267-04

Optimization and Research of Grinding Tool for Small Vertical Mill Based on EDEM

XING Zhi-zhong, ZHANG Hai-dong*, WANG Meng et al

(College of Mechanical & Electrical Engineering, Yunnan Agricultural University, Kunming, Yunnan 650201)

Abstract[Objective] To improve the crushing performance of small vertical crusher. [Method] Small vertical mill and soybean were physically measured so as to establish EDEM model. The optimization scheme of tool grinding was proposed. With the mean of rotational kinetic energy of 1.01-1.10 s particles and the collision number sum as the indexes, grinding tool performance was compared before and after optimization. Simulation experiment was carried out for three times; and the correctness of optimization was verified. [Result] After optimization, the mean of rotational kinetic energy of 1.01-1.10 s particles and the collision number sum were both greater than those of particles in the original mill. The mean maximum of particle rotational kinetic energy was 0.007 7 J; the maximum value of the sum of collision times reached 738. [Conclusion] After optimization, the performance of the grinding tool is obviously better than that of the former. This research provides theoretical references for the improvement of the grinding tool.

Key wordsEDEM; Pulverizer; Tool; Optimization