摘 ?要： 在如今的工業(yè)現(xiàn)場(chǎng)、食品業(yè)、制藥等行業(yè)，都存在著研磨這個(gè)加工流程。球磨機(jī)在研磨工序中是絕佳的選擇。由于其自身的優(yōu)勢(shì)很明顯，比如它的結(jié)構(gòu)相對(duì)比較簡(jiǎn)單，而且研磨的效率高，可靠性強(qiáng)等。使其在研磨工序中占有比較高的地位。但是球磨機(jī)的自身存在的缺陷也比較多，比如鑄造球磨機(jī)所需的鋼材較多，影響工作效率的因素較多等。本文針對(duì)球磨機(jī)存在的缺陷，對(duì)影響球磨機(jī)研磨效率的工作參數(shù)做出了相應(yīng)的研究，并利用正交實(shí)驗(yàn)的方法進(jìn)行分組實(shí)驗(yàn)，最后利用離散元的方法通過(guò)EDEM軟件建立仿真模型，利用仿真的結(jié)果對(duì)工作參數(shù)進(jìn)行調(diào)節(jié)，最終實(shí)現(xiàn)提高球磨機(jī)的工作效率的目的。

關(guān)鍵詞： 球磨機(jī);離散元;EDEM

中圖分類(lèi)號(hào)： TP3 ? ?文獻(xiàn)標(biāo)識(shí)碼： A ? ?DOI：10.3969/j.issn.1003-6970.2019.08.034

本文著錄格式：孫鵬. 基于離散元素法的巖土球磨機(jī)仿真與分析[J]. 軟件，2019，40（8）：145148+189

【Abstract】： In today's industrial field， food industry， pharmaceutical industry， there is a grinding process. Ball mill is the best choice in the grinding process. Because of its obvious advantages， such as its relatively simple structure， high grinding efficiency， strong reliability and so on. It occupies a relatively high position in the grinding process. But there are many defects in the ball mill itself， such as more steel for casting ball mill， more factors affecting work efficiency， and so on. Aiming at the shortcomings of ball mill， this paper studies the working parameters affecting the grinding efficiency of ball mill， and carries out grouping experiments by means of orthogonal experiment. Finally， the simulation model is established by EDEM software with the method of discrete element， and the working parameters are adjusted by the simulation results. Finally， the purpose of improving the working efficiency of ball mill is realized.

【Key words】： Ball mill; Discrete element; EDEM

0 ?引言

研磨工序在如今的各大行業(yè)都發(fā)揮著至關(guān)重要的作用，可以說(shuō)是在諸多行業(yè)成為了必不可少的一份工序，但是也存在著顯著的問(wèn)題，即研磨的工作效率相對(duì)太低，浪費(fèi)了很多時(shí)間成本和資源，在礦工廠的生產(chǎn)成本中，就研磨這道工序甚至占到了40%的生產(chǎn)成本[1]，所以如果能在研磨這道工序上效率有所提高，甚至只提高一小步都能為礦工廠的效益帶來(lái)巨大的提升，甚至擴(kuò)展到整個(gè)行業(yè)，為整個(gè)研磨行業(yè)帶來(lái)可觀的經(jīng)濟(jì)增長(zhǎng)[2]。

在研磨行業(yè)中，最常用的研磨設(shè)備即為球磨機(jī)，球磨機(jī)的分類(lèi)多種多樣，本文研究的主要是臥式球磨機(jī)。由于球磨機(jī)的優(yōu)勢(shì)比較突出[3]：自身的結(jié)構(gòu)以筒體為主，所以比較簡(jiǎn)單、工作過(guò)程中比較穩(wěn)定，不會(huì)出現(xiàn)大幅度抖動(dòng)破裂、工作場(chǎng)景密封性強(qiáng)，很好的減少了噪聲和灰塵、能研磨各種性質(zhì)的物料，甚至包括液態(tài)物料等等優(yōu)點(diǎn)。所以從球磨機(jī)被建造出來(lái)后，一直在研磨行業(yè)扮演著舉足輕重的地位，被廣泛的運(yùn)用到礦石、化工、食品等各種行業(yè)當(dāng)中，并且都為各行業(yè)帶來(lái)可觀的經(jīng)濟(jì)增長(zhǎng)。

據(jù)研究，影響球磨機(jī)工作效率的原因基本可以分為以下三種情況[4]：（1）球磨機(jī)自身的結(jié)構(gòu);球磨機(jī)主要是由筒體和支架組成，筒體的尺寸和性質(zhì)以及筒體內(nèi)襯板的形狀，尺寸都會(huì)對(duì)工作效率有所影響。（2）入磨物料的性質(zhì);包括固液物料分類(lèi)，物料的尺寸、數(shù)量、性質(zhì)等等。（3）球磨機(jī)運(yùn)轉(zhuǎn)的工作參數(shù);包括球磨機(jī)轉(zhuǎn)速的大小、料球比、筒體內(nèi)介質(zhì)填充率等。

根據(jù)以上分析，影響球磨機(jī)研磨效率的因素較多[5]，難以從全局角度優(yōu)化，具體分析如下。首先是球磨機(jī)自身的結(jié)構(gòu);一般在選型結(jié)束建造出來(lái)后整體尺寸就難以改變，只能在襯板上做文章，比如修改下襯板的尺寸，調(diào)整下襯板的角度和數(shù)量等，影響效果不是很大;再是入磨物料的性質(zhì)，只是作為外部因素，與球磨機(jī)自身的遠(yuǎn)轉(zhuǎn)情況關(guān)系不大，所以暫不予考慮;最后就是球磨機(jī)運(yùn)轉(zhuǎn)的工作參數(shù)，可以作為調(diào)整的和優(yōu)化方向進(jìn)行研究，也就是本文的研究?jī)?nèi)容。

1 ?離散元素法的基本理論與力學(xué)模型

1.1 ?離散元的基本理論

由于離散元最初的研究對(duì)象是巖土等一些非連續(xù)的介質(zhì)，這些介質(zhì)的運(yùn)動(dòng)狀況很難用肉眼去評(píng)判，所以離散元素法是通過(guò)將介質(zhì)劃分成相互獨(dú)立的顆粒單元，當(dāng)顆粒之間發(fā)生碰撞，也就會(huì)產(chǎn)生相應(yīng)的作用力[6]。通過(guò)牛頓第二運(yùn)動(dòng)定律，利用迭代的方法不斷迭代計(jì)算，最終可以求得在一個(gè)選定的時(shí)間步長(zhǎng)內(nèi)所有顆粒的運(yùn)動(dòng)狀況，跟蹤這些實(shí)時(shí)更新位置的顆粒，可得到一個(gè)宏觀的運(yùn)動(dòng)規(guī)律。運(yùn)用離散元素法求解問(wèn)題的時(shí)候有個(gè)前提條件：首先是選定的時(shí)間步長(zhǎng)要足夠小，在這個(gè)時(shí)間步長(zhǎng)內(nèi)，顆粒單元的碰撞不能干擾到其他單元，并且在這個(gè)時(shí)間不長(zhǎng)內(nèi)，要保持顆粒單元的加速度和速度不變。只有滿足這樣的前提條件，才能保證離散元素法的科學(xué)應(yīng)用。

在離散元素法中主要包含兩個(gè)方面的內(nèi)容[7]：一是顆粒單元之間的接觸碰撞，也就是力和位移的關(guān)系;用來(lái)求解顆粒單元碰撞時(shí)所受的力。二是牛頓第二運(yùn)動(dòng)定律的運(yùn)用;用來(lái)求解顆粒單元碰撞時(shí)經(jīng)過(guò)的位移，和顆粒單元的速度、加速度。由于離散元素法的應(yīng)用范圍很廣泛，所以接觸模型和牛頓第二定律的運(yùn)用也需要隨情況變化。一般是對(duì)于散體顆粒模型，是將單個(gè)的顆粒看做一個(gè)單元，對(duì)于巖土礦石，則是將每個(gè)不規(guī)則的巖土塊看做一個(gè)單元。

1.2 ?顆粒模型

顆粒模型根據(jù)研究的對(duì)象不同，建立的模型和計(jì)算方式也不同。常見(jiàn)的顆粒模型可以分為以下兩種：硬球模型和軟球模型[8]。這兩種模型各有優(yōu)缺點(diǎn)，適用于不同的場(chǎng)合。硬球模型主要適用的是顆粒單元運(yùn)動(dòng)較快的場(chǎng)景，在這樣的場(chǎng)景中，顆粒的碰撞是在一瞬間發(fā)生的，接觸時(shí)刻很短，顆粒單元所受的瞬時(shí)沖擊力很強(qiáng)，一般會(huì)直接破碎或者顆粒沒(méi)有明顯的變化，這種情況下，只需要考慮兩個(gè)單元顆粒之間的碰撞情況，不需要再考慮另外的碰撞干擾。而軟球模型適用的是顆粒單元運(yùn)動(dòng)速度較慢的情況，在這種情況下，顆粒之間的碰撞會(huì)持續(xù)在一段時(shí)間內(nèi)，會(huì)發(fā)生明顯的變形，通過(guò)計(jì)算顆粒之間的變形量可以得到顆粒所受到的彈性力。本文研究的是球磨機(jī)對(duì)巖土顆粒的研磨過(guò)程，所以選擇軟球模型來(lái)模擬顆粒之間的碰撞情況。

由于本文研究的顆粒對(duì)象是巖土，形狀不規(guī)則，受力不均勻，所以我們將巖土顆粒模擬成形狀相同，質(zhì)量相等的球形顆粒。利用球形顆粒之間的接觸碰撞模擬巖土研磨時(shí)發(fā)生破碎時(shí)的接觸碰撞。下圖1就是模擬的球形顆粒與顆粒之間的接觸和球星顆粒與球磨機(jī)邊界的接觸情況[9]。

將圖1的球形顆粒接觸情況中的顆粒受力情 況進(jìn)行詳細(xì)分析，可以得到如下圖2的顆粒受力形變圖：

球形顆粒之間的碰撞接觸再細(xì)化分，可以將之看作一個(gè)振動(dòng)模型[10]，顆粒碰撞時(shí)所產(chǎn)生的力與振動(dòng)模型中的力方向基本一致，這樣，我們就可以利用振動(dòng)模型中的求解方法取求解顆粒碰撞所受的力以及顆粒經(jīng)過(guò)的位移等變量。下圖3即為模擬成振動(dòng)模型的圖：

1.3 ?運(yùn)動(dòng)方程

3 ?球磨機(jī)的工作原理和參數(shù)

球磨機(jī)的工作原理就是以電能帶動(dòng)筒體高速旋轉(zhuǎn)運(yùn)動(dòng)，使得筒體內(nèi)的介質(zhì)和入磨物料混合充分，并通過(guò)相互之間的不斷碰撞破碎，達(dá)到物料被研磨的效果[13]。而影響球磨機(jī)研磨效果的因素也比較錯(cuò)綜復(fù)雜，難以全面的考量每種因素的影響效果，常見(jiàn)的分類(lèi)有幾下幾種：

（1）球磨機(jī)自身結(jié)構(gòu)的影響[14]：球磨機(jī)主要是由進(jìn)出料口、傳動(dòng)軸、筒體、支架等主要部分組成;其中筒體尺寸的大小對(duì)球磨機(jī)研磨效果的影響最深，因?yàn)橥搀w的直徑和筒長(zhǎng)決定了物料再研磨過(guò)程中的運(yùn)動(dòng)空間和運(yùn)動(dòng)路線，同時(shí)筒體的大小對(duì)球磨機(jī)的生產(chǎn)效率也有所影響，筒體越大代表著可入磨的物料數(shù)量越多，一次性生產(chǎn)量也就越多。筒體內(nèi)的襯板對(duì)物料的研磨效果也有所影響，一般不同尺寸、不同形狀、不同數(shù)目、不同材料等的襯板對(duì)物料的研磨效果都有所區(qū)別，同時(shí)還要考慮到襯板的磨損問(wèn)題，襯板如果被磨損的過(guò)多，就意味著筒體內(nèi)的直徑空間在逐漸變大，需要及時(shí)調(diào)整入磨物料的量，保證球磨機(jī)維持良好的研磨效果。

（2）工作參數(shù)影響：常見(jiàn)的球磨機(jī)的工作參數(shù)有料球比、球磨機(jī)的轉(zhuǎn)速、介質(zhì)填充率、介質(zhì)尺寸等。料球比指的是指的是入磨物料與介質(zhì)磨球的比值，如果料球比過(guò)大，意外著入磨物料過(guò)多，磨球太少，導(dǎo)致物料的破碎率降低，更多的是物料之間的無(wú)效接觸，研磨效率太低;而如果料球比過(guò)小，則意味著入磨物料過(guò)少，介質(zhì)磨球過(guò)多，會(huì)導(dǎo)致介質(zhì)磨球與筒壁和襯板的摩擦增大，有損磨機(jī)的壽命，而且使得球磨機(jī)的生產(chǎn)效率變低，更多處于空磨得狀態(tài)，浪費(fèi)生產(chǎn)資源。介質(zhì)填充率過(guò)大或過(guò)小也會(huì)影響球磨機(jī)的工作效率，如果過(guò)大，可磨的物料就會(huì)變少，球磨機(jī)處于過(guò)飽和狀態(tài)，而過(guò)小則會(huì)使得球磨機(jī)空轉(zhuǎn)的次數(shù)變多。球磨機(jī)轉(zhuǎn)速的大小對(duì)研磨效果的影響也比較大;轉(zhuǎn)速過(guò)大意外著物料受到的離心力過(guò)大，使得物料做圓周運(yùn)動(dòng)，達(dá)不到有效的研磨效果，轉(zhuǎn)速率過(guò)小的話，則會(huì)導(dǎo)致物料與介質(zhì)的碰撞破碎沖擊能量太小，不能有效的研磨物料。

（3）入磨物料的性質(zhì)：一般入磨物料的種類(lèi)比較多，很難精確評(píng)判某一種物料的可磨度，按常規(guī)來(lái)說(shuō)，如果物料的硬度較高，研磨效果就會(huì)相對(duì)較差，球磨機(jī)的生產(chǎn)效率就會(huì)較低，而如果物料的硬度較低，就比較容易發(fā)生碰撞破碎，也就意味著更容易被研磨[16]。同時(shí)，入磨物料的粒度對(duì)研磨效果也有影響，如果入磨物料的粒度太小，物料和磨礦介質(zhì)之間的碰撞產(chǎn)生的沖擊力就會(huì)比較小，研磨的效果會(huì)很不理想，而且可能還會(huì)導(dǎo)致磨機(jī)底部出現(xiàn)砂石沉積的現(xiàn)象，影響磨機(jī)的壽命;而如果粒度過(guò)大，會(huì)使得研磨的時(shí)間變長(zhǎng)，也就使得球磨機(jī)的生產(chǎn)效率變低。另外入磨物料自身的化學(xué)特性等也會(huì)對(duì)研磨效果有所影響[17]。

由以上分析可知，對(duì)球磨機(jī)研磨效果的影響因素較多，在本文中難以全面的進(jìn)行分析，所以本文主要是對(duì)影響球磨機(jī)研磨效果的工作參數(shù)做出了相關(guān)研究。

4 ?仿真實(shí)驗(yàn)分析

4.1 ?球磨機(jī)仿真參數(shù)選取

本文建立的實(shí)驗(yàn)球磨機(jī)仿真模型的尺寸為筒體直徑為800 mm，筒體長(zhǎng)度為600 mm。球磨機(jī)的結(jié)構(gòu)相對(duì)比較簡(jiǎn)單，我們主要研究的是筒體的內(nèi)容，筒體內(nèi)的襯板共十根，長(zhǎng)寬高分別為600 mm、10 mm、10 mm，如圖4所示。

本文通過(guò)離散元軟件EDEM建立起球磨機(jī)的仿真模型[18]。設(shè)定不同的轉(zhuǎn)速率、介質(zhì)尺寸、料球比和介質(zhì)填充率的大小依據(jù)正交法的法則進(jìn)行分組實(shí)驗(yàn)，根據(jù)實(shí)驗(yàn)結(jié)果可獲得每組參數(shù)的組合下球磨機(jī)的功率和入磨入料受到的沖擊能?？梢愿鶕?jù)單位質(zhì)量的功率和沖擊能進(jìn)行分析對(duì)比，最終可得出最佳的工作參數(shù)組合。

4.2 ?正交試驗(yàn)

根據(jù)以上選取的因素和水平，再結(jié)合第三章 對(duì)各因素水平的取值可以得到如下表1的因素水平表[19]。

在上表1中，我們?cè)O(shè)定了轉(zhuǎn)速率、介質(zhì)填充率、介質(zhì)尺寸、料球比四個(gè)因素，每個(gè)因素設(shè)置了三個(gè)水平值，所以可將正交表設(shè)為 ，也就是四因素三水平共九組試驗(yàn)[20]，如表2所示。

4.3 ?利用EDEM進(jìn)行仿真

對(duì)正交表九組試驗(yàn)利用EDEM仿真可得到如下九組結(jié)果：

5 ?結(jié)論

通過(guò)EDEM仿真結(jié)果，可以看出每一種參數(shù)組合導(dǎo)致的不同的顆粒運(yùn)動(dòng)效果，再利用極差分析法和方差分析法，分別評(píng)估了影響球磨機(jī)工作性能的四種參數(shù)的重要性，最終得出球磨機(jī)的最優(yōu)組合為：介質(zhì)填充率40%，介質(zhì)尺寸40 mm，料球比0.6，轉(zhuǎn)速率80%。在這樣的工作參數(shù)組合下，球磨機(jī)單位質(zhì)量損耗的功率和介質(zhì)對(duì)物料的沖擊能量都會(huì)最大。

參考文獻(xiàn)

[1] 周恩浦. 礦山機(jī)械：選礦機(jī)械部分[M]. 冶金工業(yè)出版社， 1979.

[2] 于雄， 張株瑞， 周應(yīng)軍. 離散制造企業(yè)設(shè)備管理系統(tǒng)與MES集成研究[J]. 軟件， 2016， 37（02）： 174-176

[3] Testing of a new dynamic Ergun equation for transport with positron emission particle tracking[J]. AIChE Journal， 2016， 62（3）： 939-946.

[4] 王連杰， 韋群. 三維云圖仿真系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)[J]. 軟件， 2018， 39（5）： 151-155.

[5] 何云. 球磨機(jī)磨碎介質(zhì)的合理選擇[J]. 硫磷設(shè)計(jì)與粉體工程， 1993（1）： 6-19.

[6] 段希祥. 我國(guó)粗磨球磨機(jī)鋼球尺寸狀況的分析[J]. 礦冶工程， 1998（1）： 23-26.

[7] 段希祥. 球磨機(jī)鋼球尺寸的理論計(jì)算研究[J]. 中國(guó)科學(xué)（A輯 數(shù)學(xué) 物理學(xué) 天文學(xué) 技術(shù)科學(xué)）， 1989， 32（8）： 856-863.

[8] 包敬海. 基于Matlab混合編程的海上溢油仿真系統(tǒng)的研究[J]. 軟件， 2012， 33（7）： 103-104.

[9] 孫軍鋒， 董為民， 姚輝， et al. 襯板提升條對(duì)磨機(jī)有用功率的影響[J]. 礦山機(jī)械， 2009（21）： 56-60.

[10] Cundall， P. A. （1971） A Computer Model for Simulating Progressive Large Scale Movements in Blocky Rock System. International Proceedings Symposium， 128-132.

[11] 陳媛媛， 黃善國(guó)， 郭嬋. 一種基于ALM的Modem自動(dòng)化測(cè)試系統(tǒng)[J]. 軟件， 2012， 33（6）： 104-106.

[12] R. D. Morrison， P. W. Cleary. Using DEM to Model ore breakage within a pilot scale SAG mill.Minerals Engineering. 2004， 17： 1117-1124.

[13] 王建紅. 基于零均值及零方差圖構(gòu)建的潤(rùn)滑油磨粒檢測(cè) [J]. 軟件， 2015， 36（6）： 120-124.

[14] Lemos J. V. A Hybrid Distinct Element Computational Model For The Half-Plane [M.S] University of Minnesota， 1983.

[15] CUNDALL P. A.， STRACK O D L.A discrete numerical model for granular assemble. Geotechnique， 1979， 29（1）： 47-65.

[16] M.S. Poowell. GN. Nurick. A. Study of Charge Motion in Rotary Mills Part 1-Extension of the Theory. Minerals engineering. 1996（9）： 259-268.

[17] M. S. Poowell. GN. Nurick. A. Study of Charge Motion in Rotary Mills Part-2-Experiment Work.Minerals Engineering. 1996（9）： 343-350

[18] CUNDALL P. A.， STRACK O. D. L. The distinct element meted as a toll for research in granular media-paet： department of civil and mineral engineering.University of Minnesota， national Science Foundation， 1979.

[19] P. W. Cleary. Charge behaviour and power consumption in ball mills： sensitivity to mill opereating conditions， liner geometry and charge composition， International Journal of Mineral Processing， 2001， 63： 79-114.

[20] R. Changramohan， M. S. Powell. Measurenent of particle interaction properties for incorporation in discrete element method simulation. Minerals.