李 臣，姜志宏，蔡改貧

（江西理工大學(xué) 機(jī)電工程學(xué)院，江西 贛州 341000）

球磨機(jī)腔形結(jié)構(gòu)優(yōu)化及其數(shù)值模擬分析

李 臣，姜志宏，蔡改貧

（江西理工大學(xué) 機(jī)電工程學(xué)院，江西 贛州 341000）

球磨機(jī)作為應(yīng)用廣泛的粉磨機(jī)械，其磨礦效果一直是研究人員關(guān)注的問題。本研究從球磨機(jī)腔形結(jié)構(gòu)出發(fā)，提出了一種橢圓腔形結(jié)構(gòu)的新型球磨機(jī)。并針對物料在球磨機(jī)內(nèi)的運(yùn)動特性，建立了物料運(yùn)動的動力學(xué)模型。采用離散元軟件EDEM，從轉(zhuǎn)速差、粉磨效率、物料之間以及物料與磨料介質(zhì)之間的碰撞摩擦能三個方面，對所提出的橢圓腔形結(jié)構(gòu)球磨機(jī)與傳統(tǒng)的圓筒腔形球磨機(jī)進(jìn)行數(shù)值模擬對比分析。結(jié)果表明，所提出的橢圓腔形結(jié)構(gòu)球磨機(jī)的磨礦效果大幅提升。物料相對于球磨機(jī)筒體的轉(zhuǎn)速差增大；物料的粉磨效率提高；球磨機(jī)內(nèi)物料死區(qū)減小，物料之間以及物料與磨料介質(zhì)之間的碰撞摩擦增多，碰撞摩擦能增大。

球磨機(jī)腔形；離散元；磨礦效果；二次開發(fā)

球磨機(jī)在物料粉磨過程中具有廣泛的應(yīng)用，球磨機(jī)的研究在國內(nèi)外一直受到廣泛的關(guān)注和高度的重視[1-2]。國內(nèi)外學(xué)者和工程技術(shù)人員對球磨機(jī)的研究取得了一系列成果[3-4]，但對于球磨機(jī)腔形結(jié)構(gòu)的研究卻相對較少。相關(guān)研究發(fā)現(xiàn)球磨機(jī)腔形結(jié)構(gòu)對物料的粉磨效果有重要影響[5]。球磨機(jī)的基本工作原理是球磨機(jī)筒體的旋轉(zhuǎn)，帶動物料與鋼球在筒體內(nèi)作拋落運(yùn)動，利用物料之間以及物料與磨料介質(zhì)之間的碰撞摩擦使物料細(xì)化的過程。

本研究從球磨機(jī)的腔形結(jié)構(gòu)出發(fā)，提出了一種橢圓腔形結(jié)構(gòu)的新型球磨機(jī)。首先分析物料在球磨機(jī)筒體內(nèi)的運(yùn)動，建立物料運(yùn)動的動力學(xué)模型。利用離散元軟件EDEM及其二次開發(fā)，從轉(zhuǎn)速差、磨礦效率、物料之間以及物料與磨料介質(zhì)之間的碰撞摩擦能三個方面，對所提出的橢圓腔形結(jié)構(gòu)球磨機(jī)與傳統(tǒng)的圓筒形球磨機(jī)進(jìn)行數(shù)值模擬對比分析。為研發(fā)設(shè)計(jì)新型高效球磨機(jī)，提高球磨機(jī)的磨礦效果提供一定的參考。

1 球磨機(jī)物理模型

為便于仿真模擬，簡化球磨機(jī)筒體結(jié)構(gòu)模型由筒體、內(nèi)部襯板組成（圖1）。橢圓腔形截面尺寸如圖2所示。本研究所提出的橢圓腔形球磨機(jī)截面由兩個相同橢圓正交疊加組合而成。兩橢圓相交所得最外層曲線即為橢圓腔形截面的輪廓（圖2）。圖3為圓筒腔形球磨機(jī)腔形結(jié)構(gòu)示意圖。

圖1 橢圓腔形球磨機(jī)筒體Fig.1 Elliptical chamber ball mill cylinder

圖2 橢圓腔形截面尺寸Fig.2 Size of elliptical chamber’section

圖3 圓筒腔形球磨機(jī)筒體Fig.3 Cylindrical chamber ball mill cylinder

當(dāng)球磨機(jī)以一定的速度做旋轉(zhuǎn)運(yùn)動時，筒壁底層的物料在襯板、物料間的摩擦力和物料重力以及由于筒體旋轉(zhuǎn)而產(chǎn)生的離心力的作用下，將隨著筒體作旋轉(zhuǎn)上升運(yùn)動，被提升到一定的高度，當(dāng)顆粒的重力大于或等于離心力時，就開始脫離筒體內(nèi)壁，按照某一軌跡降落。這種周而復(fù)始的運(yùn)動就產(chǎn)生了連續(xù)的沖擊和研磨作用，從而達(dá)到粉磨的目的，并呈現(xiàn)出復(fù)雜的物料混合運(yùn)動[6]。且在運(yùn)動過程中，物料在筒體內(nèi)出現(xiàn)分區(qū)現(xiàn)象，通常將物料在筒體內(nèi)分3個運(yùn)動區(qū)域，分別為隨筒體上升區(qū)、瀉落面滑落區(qū)（此區(qū)域磨礦效果最差，又叫死區(qū)）、拋落區(qū)[7]，如圖4所示。

圖4 物料在筒體的分布Fig.4 Material distribution in the cylinder

2 物料運(yùn)動動力學(xué)模型

2.1 物料運(yùn)動方程

根據(jù)牛頓第二定律，建立物料在球磨機(jī)筒體內(nèi)的平動和轉(zhuǎn)動運(yùn)動方程[8]：

式中：v'為物料加速度，m/s2；w'為物料角加速度，rad/s2；F為物料所受的合力（包括重力、沖擊力、摩擦力），N；T為物料所受的合力矩，N·m；I為物料的轉(zhuǎn)動慣量，kg·m2；m為物料質(zhì)量，kg。

2.2 物料之間、物料與磨料介質(zhì)的受力分析

物料之間的碰撞，在接觸點(diǎn)處的受力Fc可分解為切向分力和法向分力；對于任一兩個碰撞物料i，j的接觸力計(jì)算公式如下[9]：

式中：Kτ，Kn分別為物料之間的切向和法向剛度系數(shù)，為常量；δτi，j，δni，j分別為物料在碰撞接觸點(diǎn)處的切向和法向變形量，mm；ητ，ηn分別為物料之間的切向和法向阻尼系數(shù)，為常量；ni，j為物料接觸點(diǎn)處的法向矢量；vri，j=vi-vj為物料碰撞時的相對速度，m/s；vsi，j=vri，j-（vri，j×ni，j）×ni，j+R（wi+wj）×ni，j為物料碰撞時接觸點(diǎn)處的滑動速度，m/s；wi、wj分別為兩碰撞物料角速度，rad/s；R為物料間接觸半徑，mm。

筒體內(nèi)所有物料之間、物料i與鋼球j以及物料i與筒體j間的合力統(tǒng)一表達(dá)式如式（5）所示：

2.3 碰撞摩擦能計(jì)算模型

離散元軟件EDEM中無法直接分析物料之間、物料與磨料介質(zhì)之間的碰撞摩擦能量。通過EDEM二次開發(fā)，將推導(dǎo)的碰撞摩擦能計(jì)算公式嵌入EDEM軟件中，進(jìn)行物料的碰撞摩擦分析。碰撞摩擦能主要包括碰撞摩擦引起的摩擦能和阻尼能[10]。

式中：N為碎磨過程中的相互摩擦對象，包括物料、鋼球和筒體；Ff為相互摩擦對象之間的摩擦力，N；dxslip為不考慮彈性變形的位移增量，mm。

式中：N為碎磨過程中相互碰撞的對象，包括物料、鋼球和筒體；cn為碰撞對象之間法向阻尼系數(shù)，為常量；x'n為碰撞對象之間法向速度，m/s；dxn為碰撞對象之間法向位移增量mm；ct為碰撞對象之間的切向阻尼系數(shù)，為常量；x't為碰撞對象之間切向速度，m/s；dxt為碰撞對象之間法向位移增量，mm。

3 球磨機(jī)腔形結(jié)構(gòu)數(shù)值模擬對比分析

通過在相同工況下的對比試驗(yàn)，分別從轉(zhuǎn)速差、粉磨效率、物料之間以及物料與磨料介質(zhì)之間的碰撞摩擦能三個方面，對所提出的橢圓腔形結(jié)構(gòu)球磨機(jī)與傳統(tǒng)的圓筒腔形球磨機(jī)進(jìn)行數(shù)值模擬對比分析。球磨機(jī)參數(shù)如表1、表2所示。

EDEM是使用離散元技術(shù)進(jìn)行顆粒系統(tǒng)仿真和分析的通用CAE軟件，可建立礦石顆粒的參數(shù)化模型[11]。本研究中物料的模擬采用多尺度內(nèi)聚物料模型[12-13]，物料內(nèi)部由黏結(jié)鍵連接，黏結(jié)鍵的斷裂過程就是物料粉磨過程，具體物料相關(guān)參數(shù)如表3所示。

按照試驗(yàn)參數(shù)，在EDEM軟件中進(jìn)行相關(guān)參數(shù)設(shè)置，設(shè)置仿真時間步長為1×10-6s仿真總時間為10 s。圖5、圖6分別為EDEM仿真過程中橢圓腔形和圓筒腔形球磨機(jī)內(nèi)不同時刻磨礦狀態(tài)。

表1 橢圓腔形球磨機(jī)參數(shù)Tab.1 Parameters of elliptical chamber ball mill

表2 圓筒腔形球磨機(jī)參數(shù)Tab.2 Parameters of cylindrical chamber ball mill

表3 物料參數(shù)Tab.3 Material parameters

圖5 橢圓腔形球磨機(jī)內(nèi)不同時刻磨礦狀態(tài)Fig.5 Grinding state of the Elliptical chamber ball mill at different time

圖6 圓筒腔形球磨機(jī)內(nèi)不同時刻磨礦狀態(tài)Fig.6 Grinding state of cylindrical chamber ball mill at different time

3.1 轉(zhuǎn)速差對比分析

轉(zhuǎn)速差為物料在筒體內(nèi)運(yùn)動，相對于筒體的速度差。轉(zhuǎn)速差越大，則物料在筒體內(nèi)的動能越大，物料受到的碰撞摩擦力越大，物料的粉磨效果越好。通過EDEM后處理，獲得橢圓腔形和圓筒腔形球磨機(jī)內(nèi)物料轉(zhuǎn)速差的變化曲線，如圖7所示。

圖7表明，本研究提出的改進(jìn)橢圓腔形球磨機(jī)，在相同的工作參數(shù)下，物料的整體轉(zhuǎn)速差高于傳統(tǒng)圓筒腔形球磨機(jī)，且峰值轉(zhuǎn)速差更大。說明改進(jìn)后的橢圓腔形球磨機(jī)能使物料在筒體內(nèi)獲得更大的動能，所受的碰撞摩擦力更大。

圖7 球磨機(jī)內(nèi)物料轉(zhuǎn)速差對比Fig.7 Comparisonofmaterialrotatingspeeddifferenceintheballmill

3.2 磨礦效率對比分析

本研究在物料粉磨數(shù)值模擬中采用的多尺度內(nèi)聚物料模型，可以通過物料內(nèi)部黏結(jié)鍵數(shù)目變化來描述物料的粉磨效果。通過EDEM后處理，分別獲取橢圓腔形和圓筒腔形球磨機(jī)磨礦過程中物料內(nèi)部黏結(jié)鍵數(shù)目的變化曲線，如圖8所示。

圖8 球磨機(jī)內(nèi)物料內(nèi)部黏結(jié)鍵斷裂趨勢對比Fig.8 Comparisonoffracturetendencyofinternalbondsinballmill

從圖8可以得出，相同工作參數(shù)下，傳統(tǒng)圓筒腔形球磨機(jī)中物料內(nèi)部黏結(jié)鍵斷裂呈現(xiàn)出線性化變化趨勢，且物料內(nèi)部黏結(jié)鍵斷裂個數(shù)少于改進(jìn)橢圓腔形球磨機(jī)。改進(jìn)橢圓腔形球磨機(jī)中物料內(nèi)部黏結(jié)鍵的斷裂呈現(xiàn)出非線性變化趨勢，物料內(nèi)部黏結(jié)鍵斷裂的速度整體高于傳統(tǒng)圓筒腔形球磨機(jī)，且在4～6 s時間段物料內(nèi)部黏結(jié)鍵斷裂速度最快。表明改進(jìn)橢圓腔形球磨機(jī)使得物料在球磨機(jī)內(nèi)的磨礦效率獲得了提升。

3.3 碰撞摩擦能對比分析

球磨機(jī)的磨礦效果取決于物料之間、磨料介質(zhì)與物料之間的碰撞摩擦，碰撞摩擦產(chǎn)生的能量越大，說明物料受到的碰撞、摩擦作用越多，越有利于物料的粉磨。通過EDEM二次開發(fā)，將碰撞摩擦能計(jì)算公式（6）、（7）、（8）內(nèi)嵌到EDEM中，通過后處理分別獲取物料在橢圓腔形和圓筒腔形球磨機(jī)內(nèi)碰撞摩擦能的變化趨勢，如圖9所示。

圖9 球磨機(jī)內(nèi)碰撞摩擦能對比Fig.9 Comparison of collision and friction energy in ball mill

圖9分析表明，傳統(tǒng)圓筒腔形球磨機(jī)內(nèi)物料碰撞摩擦能增長緩慢，且近似呈現(xiàn)出線性化增長趨勢。改進(jìn)橢圓腔形球磨機(jī)內(nèi)物料碰撞摩擦能增長迅速，且在4～6 s時間段增長最快。碰撞摩擦能的變化趨勢與物料內(nèi)部黏結(jié)鍵斷裂變化趨勢基本吻合。說明改進(jìn)橢圓腔形球磨機(jī)內(nèi)物料受到更多的碰撞、摩擦作用，磨礦效果更好。

3.4 球磨機(jī)內(nèi)死區(qū)對比

球磨機(jī)內(nèi)物料的運(yùn)動區(qū)域劃分為：隨筒體上升區(qū)、瀉落面滑落區(qū)（死區(qū)）、拋落區(qū)。增大拋落區(qū)的范圍，減小死區(qū)的范圍是提高磨礦效果的有效方式。對比數(shù)值模擬中橢圓腔形和圓筒腔形球磨機(jī)內(nèi)死區(qū)的大小，如圖10、圖11所示。

對比圖10與圖11得出，相比于傳統(tǒng)圓筒腔形，改進(jìn)橢圓腔形內(nèi)物料運(yùn)動死區(qū)明顯減少，更多的物料處于拋落區(qū)。說明改進(jìn)橢圓腔形球磨機(jī)能明顯減少死區(qū)的大小，提高粉磨效果。

圖10 圓筒腔形物料運(yùn)動矢量圖Fig.10 Motion vector in cylindrical chamber

圖11 橢圓腔形物料運(yùn)動矢量圖Fig.11 Motion vectors in elliptical chamber

4 結(jié)論

本研究從球磨機(jī)腔形結(jié)構(gòu)出發(fā)，提出了一種橢圓腔形結(jié)構(gòu)的新型球磨機(jī)。通過離散元軟件EDEM及其二次開發(fā)，對所提出的橢圓腔形結(jié)構(gòu)球磨機(jī)與傳統(tǒng)的圓筒腔形結(jié)構(gòu)球磨機(jī)進(jìn)行數(shù)值模擬對比分析，得出結(jié)論如下：

（1）橢圓腔形結(jié)構(gòu)球磨機(jī)內(nèi)物料相對于球磨機(jī)筒體的轉(zhuǎn)速差增大，物料在筒體內(nèi)獲得更大的動能，所受的碰撞摩擦力更大。

（2）橢圓腔形結(jié)構(gòu)球磨機(jī)內(nèi)物料的粉磨效率提高，物料內(nèi)部黏結(jié)鍵的斷裂呈現(xiàn)出非線性變化趨勢，物料內(nèi)部黏結(jié)鍵斷裂的速度更快、斷裂個數(shù)更多。

（3）橢圓腔形結(jié)構(gòu)球磨機(jī)內(nèi)物料死區(qū)減小，物料之間以及物料與磨料介質(zhì)之間的碰撞、摩擦作用增多，碰撞摩擦能增大。

[1] 王繼杰.球磨機(jī)的研究進(jìn)展及實(shí)用技術(shù)綜述[J].中國陶瓷工業(yè)，2007，14（5）：32-35. WANG Jijie.The research progress and the practical technology of the ball mill[J].Ceramic Industry of China，2007，14（5）：32-35.

[2] 程啟明，鄭 勇.球磨機(jī)多模型PID型神經(jīng)元網(wǎng)絡(luò)控制系統(tǒng)[J].中國電機(jī)工程學(xué)報(bào)，2008，28（2）：103-109. CHENG Qiming，ZHENG Yong.Control system of multi-model PID neuron network for ball mill[J].Proceedings of the CSEE，2008，28（2）：103-109.

[3] 李文亮，楊 濤，于向軍，等.國外大型球磨機(jī)發(fā)展現(xiàn)狀[J].礦山機(jī)械，2007，61（1）：13-15. LI Wenliang，YANG Tao，YU Xiangjun，et al.Development status of large scale ball mill in foreign countries[J].Mining&Processing Equipment，2007，61（1）：13-15.

[4] 耿 凡，袁竹林，孟德才，等.球磨機(jī)中顆?；旌线\(yùn)動的數(shù)值模擬[J].熱能動力工程，2009，24（5）：623-630. GENG Fan，YUAN Zhulin，MENG Decai，et al.Numerical simulation of the mixing movement of particle in a ball mill[J]. Journal of EngineeringforThermal Energy and Power，2009，24（5）：623-630.

[5] 畢秋實(shí).大型旋回破碎機(jī)腔形優(yōu)化設(shè)計(jì)及離散元仿真驗(yàn)證[D].長春：吉林大學(xué)，2014. BI Qiushi.Optimization design of large gyratory crusher chamber shape and DEM simulation for its performance[D].Changchun：Jilin University，2014.

[6] 孫利波.球磨過程的數(shù)學(xué)模型及其試驗(yàn)研究[D].濟(jì)南：山東大學(xué)，2006. SUN Libo.Model for ball mill and test investigation of its process [D].Jinan：Shandong University，2006.

[7] 樊狄鋒，孫 毅，毛亞郎，等.球磨機(jī)混合運(yùn)動狀態(tài)下介質(zhì)運(yùn)動形態(tài)的分析研究[J].礦山機(jī)械，2010，38（5）：81-84. FAN Difeng，SUN Yi，MAO Yalang，et al.Analysis of medium motions under the hybrid motion state in the ball mill[J].Mining& Processing Equipment，2010，38（5）：81-84.

[8] 黃理興.巖石動力學(xué)研究成就與趨勢 [J].巖土力學(xué)，2011，32（10）：2890-2900. HUANG Lixing.Development and new achievements of rock dynamics in China[J].Rock and Soil Mechanics，2011，32（10）：2890-2900.

[9] TSUJI Y，KAWAGUCHI T，TANAKA T.Discrete particle simulation of two-dimensional fluidized bed[J].Powder Technology，1993，77（1）：79-87.

[10] El Shamy U，DENISSEN C.Microscale characterization of energy dissipation mechanisms in liquefiable granular soils[J].Computers and Geotechnics，2010，37（7）：846-857.

[11]張 輝，張永震.顆粒力學(xué)仿真軟件EDEM簡要介紹[J].CADCAM與制造業(yè)信息化，2008，32（12）：48-49. ZHANG Hui，ZHANG Yongzhen.Brief introduction to particle mechanicssimulation softwareEDEM [J].CAD-CAM and Manufacturing Informatization，2008，32（12）：48-49.

[12]余世科.基于離散元法的多尺度內(nèi)聚顆粒模型構(gòu)建及沖擊破碎能耗分析[D].贛州：江西理工大學(xué)，2015. YU Shike.The establishment of multi-scale bonded particle model based on discrete element method and impact crushing energy consumption analysis[D].Ganzhou：Jiangxi University of Science and Technology，2015.

[13] 郭進(jìn)山，蔡改貧，盧小江，等.一種多尺度內(nèi)聚顆粒模型的構(gòu)建方法及數(shù)值模擬研究[J].科學(xué)技術(shù)與工程，2016，16（8）：72-78. GUO Jinshan，CAI Gaipin，LU Xiaojiang，et al.Investigation on construction method and numerical simulation of multi-scale cohesive particle model for rock [J].Science Technology and Engineering，2016，16（8）：72-78.

Optimization of Ball Mill's Chamber Structure and Its Numerical Simulation Analysis

LI Chen,JIANG Zhihong,CAI Gaipin
(School of Mechanical and Electrical Engineering,Jiangxi University of Science and Technology,Ganzhou 341000,Jiangxi,China)

The grinding effect for the ball mill has long been a research focus.A new type of ball mill with elliptical chamber structure is developed.The dynamic model of the material movement was established on the basis of the movement characteristics of material in the ball mill.Numerical simulation comparison analysis is performed on the elliptical-chamber structured ball mill and the traditional cylinder-chamber one from the aspects of rotational speed difference,grinding efficiency,impact friction energy between material and abrasive medium by using the discrete element software EDEM.The results show that the grinding effect of elliptical chamber structure of ball mill improves dramatically with increasing grinding efficiency,speed difference and friction energy.In addition,the dead zone of material in ball mill decreases while the friction and collision between the material and the abrasive medium increases.

ball mill chamber;discrete element;grinding effect;secondary development

TD453

A

10.3969/j.issn.1009－0622.2017.01.014

2016-11-11

國家自然科學(xué)基金項(xiàng)目（51464017）；江西省高?？萍悸涞赜?jì)劃項(xiàng)目（KJLD13045）

李 臣（1991-），男，山東臨沂人，碩士研究生，研究方向?yàn)榈V山機(jī)械及物料高效破碎。

姜志宏（1977-），男，新疆奎屯人，副教授，主要從事固體力學(xué)、物料高效破碎技術(shù)研究工作。