基于EDEM的雙齒輥破碎機(jī)破碎效率影響因素仿真分析

張瑞新 劉 煜 鄭群飛 孫健東 武海龍 劉文文

(1.中國(guó)礦業(yè)大學(xué)(北京)資源與安全工程學(xué)院，北京 100083；2.華北科技學(xué)院，河北 廊坊 065201)

隨著我國(guó)露天采礦生產(chǎn)能力的提高，露天采礦中巖石剝離物剝離量也隨之增大，為滿足剝離物膠帶運(yùn)輸對(duì)粒度的需求，大多物料需要進(jìn)行破碎，破碎設(shè)備能力的發(fā)揮直接關(guān)系到整個(gè)礦山的生產(chǎn)效能和綜合效益[1-2]。我國(guó)礦山地質(zhì)條件錯(cuò)綜復(fù)雜，面對(duì)的破碎對(duì)象特征多變，對(duì)雙齒輥破碎機(jī)破碎效率有較大影響。離散元素法在仿真顆粒破碎方面的引入為研究雙齒輥破碎機(jī)提供了新的理論依據(jù)和創(chuàng)新途徑[3]。近年來，計(jì)算機(jī)仿真技術(shù)取得了巨大進(jìn)步，利用模擬仿真對(duì)雙齒輥破碎機(jī)進(jìn)行關(guān)鍵技術(shù)研究突破了以往研究中的觀測(cè)難、統(tǒng)計(jì)難、計(jì)算量大、誤差大和耗時(shí)長(zhǎng)等弊端，因此對(duì)破碎機(jī)建立計(jì)算機(jī)模擬仿真技術(shù)平臺(tái)逐漸成為一種常用方法[4-6]。本研究以MMD1000系列雙齒輥篩分破碎機(jī)為研究對(duì)象，應(yīng)用離散元軟件EDEM對(duì)雙齒輥破碎機(jī)破碎效率進(jìn)行仿真研究，考察物料性質(zhì)、物料粒度分布、齒輥轉(zhuǎn)速等變量對(duì)破碎機(jī)破碎效率的影響規(guī)律，為提高雙齒輥破碎機(jī)生產(chǎn)效能和破碎機(jī)選型設(shè)計(jì)提供參考。

1 模型構(gòu)建

1.1 破碎機(jī)理和破碎機(jī)模型

雙齒輥破碎機(jī)的破碎是利用安裝在輥?zhàn)由系母邚?qiáng)度破碎齒進(jìn)行剪切、扭彎物料，迫使物料產(chǎn)生新的裂隙，裂隙之間發(fā)生貫通，從而產(chǎn)生破壞[6]。雙齒輥破碎機(jī)兼有一定的篩分功能，對(duì)物料的篩分是通過交錯(cuò)排列的破碎齒使合格粒級(jí)物料通過相對(duì)緩慢運(yùn)動(dòng)的兩根齒輥之間的空隙和破碎齒輥與側(cè)梳板之間的空隙直接排出，達(dá)到篩分的目的[7]。雙齒輥破碎機(jī)的破碎過程可分為初始破碎、深度破碎和第三級(jí)破碎。各破碎過程機(jī)理與特點(diǎn)如表1所示，具體破碎過程如圖1所示：

表1 各級(jí)破碎過程的破碎機(jī)理與特點(diǎn)Table 1 Crushing mechanism and characteristics of crushing process at all levels

圖1 各級(jí)破碎過程示意Fig.1 Sketch map of crushing process at all levels

結(jié)合雙齒輥破碎機(jī)的破碎機(jī)理和MMD1000系列雙齒輥篩分破碎機(jī)結(jié)構(gòu)參數(shù)(部分參數(shù)如表2所示)，并綜合考慮到研究對(duì)象傳動(dòng)結(jié)構(gòu)的復(fù)雜性以及模擬仿真時(shí)EDEM軟件的特點(diǎn)等諸多因素，在不影響仿真的前提下，利用三維建模軟件SolidWorks對(duì)該系列雙齒輥破碎機(jī)的傳動(dòng)結(jié)構(gòu)等進(jìn)行簡(jiǎn)化處理，建立了雙齒輥破碎機(jī)的簡(jiǎn)化三維模型，如圖2所示。

表2 破碎機(jī)部分參數(shù)Table 2 Parameters of crusher

1.2 顆粒黏結(jié)模型建立

顆粒接觸模型是離散單元法的重要基礎(chǔ)，仿真中根據(jù)礦石物料破碎特點(diǎn)，需要考慮多個(gè)顆粒間的碰撞，因此本文采用Hertz-Mindlin with bonding built-in黏結(jié)接觸模型，該黏結(jié)接觸模型對(duì)于礦巖結(jié)構(gòu)建模具有普遍適用性[8]。黏結(jié)模型中小顆粒間通過設(shè)定的Bonding鍵黏結(jié)在一起，以達(dá)到礦石模型的內(nèi)部力學(xué)結(jié)構(gòu)特征，Bonding鍵可以阻止顆粒的相對(duì)運(yùn)動(dòng)，當(dāng)施加在物料上的外力達(dá)到物料所能承受的最大應(yīng)力時(shí)這種結(jié)合就被破壞，從而實(shí)現(xiàn)物料破碎[9]。該黏結(jié)接觸模型設(shè)定顆粒在某一時(shí)刻ts被黏結(jié)起來，在此之前，顆粒通過默認(rèn)的Hertz-Mindlin接觸模型產(chǎn)生相互作用，法向和切向黏結(jié)力Fn、Ft及法向和切向黏結(jié)力矩Tn、Tt隨著時(shí)間步從零開始逐漸增加：

圖2 雙齒輥破碎機(jī)三維模型Fig.2 3D model of Double-Geared Roller Crusher

δFn=-υnSnAδt，

(1)

δFt=-υtStAδt，

(2)

δTn=-ωnStJδt，

(3)

(4)

當(dāng)法向和切向應(yīng)力超過某個(gè)特定的值時(shí)，黏結(jié)就被破壞。因此，定義法向和切向應(yīng)力的最大值分別為

(5)

(6)

礦石模型的尺寸主要根據(jù)雙齒輥破碎機(jī)的進(jìn)料尺寸和不同尺寸物料占比情況確定，按照一定的系數(shù)進(jìn)行顆粒填充，填充公式為

α×VRel=N×VFraction

(7)

式中：α為填充體積分?jǐn)?shù)，其值不宜過大，一般取值為0.56；VRel為實(shí)際物料顆粒的體積；N為小顆粒的體積；VFraction為填充所需的小顆粒的數(shù)量[10]。

由于本次仿真時(shí)的破碎機(jī)排料尺寸在300 mm以下，經(jīng)多次前期仿真分析，當(dāng)小顆粒的半徑設(shè)為60 mm時(shí)，既能保證仿真的真實(shí)性又能提高仿真效率。

2 仿真方案確定及參數(shù)設(shè)定

結(jié)合影響破碎機(jī)破碎效率的因素，對(duì)MMD1000系列雙齒輥破碎機(jī)在不同影響因素條件下的破碎效率進(jìn)行仿真模擬，仿真方案如圖3所示。

圖3 仿真方案示意Fig.3 Schematic diagram of simulation program

通過在EDEM中設(shè)定物料的單位面積法向剛度、單位面積切向剛度、法向極限強(qiáng)度、切向極限強(qiáng)度定義物料性質(zhì)，物料性質(zhì)參數(shù)如表3所示。物料粒度分布情況根據(jù)黑岱溝露天礦5#破碎站的物料粒度統(tǒng)計(jì)數(shù)據(jù)而定。最大入料粒度尺寸為1 500 mm，粒度在0～300 mm的物料占比大于75%，粒度在300～ 1 000 mm的物料占比小于20%，粒度大于 1 000 mm的物料占比低于5%。將物料粒度大于300 mm的定義為大尺寸物料，通過調(diào)節(jié)不同粒度顆粒數(shù)目實(shí)現(xiàn)不同物料粒度分布，將仿真物料粒度分布劃分為A、B、C 3組，其中A組物料大塊率為10%，B組物料大塊率為15%，C組物料大塊率為20%。

表3 物料性質(zhì)參數(shù)Table 3 Material properties parameters

雙齒輥破碎機(jī)為了適應(yīng)不同的來料情況，齒輥轉(zhuǎn)速可進(jìn)行適當(dāng)調(diào)節(jié)，調(diào)節(jié)范圍為設(shè)計(jì)時(shí)齒輥轉(zhuǎn)速前后富余量的5%～6%。針對(duì)MMD1000系列雙齒輥破碎機(jī)相關(guān)技術(shù)數(shù)據(jù)，確定破碎機(jī)齒輥轉(zhuǎn)速，仿真中將齒輥轉(zhuǎn)速設(shè)為31 r/min、33 r/min、35 r/min 3級(jí)。

時(shí)間步長(zhǎng)是仿真計(jì)算模塊里每2次運(yùn)算之間的時(shí)間差值。如果時(shí)間步長(zhǎng)設(shè)置不合理，將會(huì)導(dǎo)致失真，使仿真不準(zhǔn)確，而且可能出現(xiàn)替換后的顆粒“炸開”。為了保證仿真的連續(xù)性，固定時(shí)間步長(zhǎng)為Reyleigh時(shí)間步長(zhǎng)的5%～40%，該時(shí)間步長(zhǎng)為EDEM仿真準(zhǔn)靜態(tài)微粒集合中時(shí)間步長(zhǎng)在理論上的最大值，由式(8)決定：

(8)

式中：TR為瑞利時(shí)間步長(zhǎng)；R為顆粒半徑；ρ是顆粒材料密度；G是顆粒材料的剪切模量；σ是泊松比。計(jì)算可得，本仿真的時(shí)間步長(zhǎng)為1.65×10-6s，此時(shí)固定時(shí)間步長(zhǎng)為Reyleigh時(shí)間步的18.37%。

3 仿真結(jié)果分析

經(jīng)確定仿真方案和相關(guān)參數(shù)后，導(dǎo)入雙齒輥破碎機(jī)三維模型并進(jìn)行仿真，待仿真結(jié)束后，在EDEM后處理模塊輸出仿真數(shù)據(jù)并進(jìn)行數(shù)據(jù)處理和分析。由于仿真中Bonding鍵的斷裂是隨機(jī)發(fā)生的，為了清晰表達(dá)出Bonding鍵斷裂數(shù)目隨時(shí)間的變化關(guān)系，采用Bonding鍵數(shù)目在不同破碎時(shí)間內(nèi)的疊加圖表示各因素對(duì)破碎機(jī)破碎效率的影響。圖4為不同粒度分布的軟巖物料在不同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目隨時(shí)間變化的曲線，圖5為不同粒度分布的硬巖物料在不同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目隨時(shí)間變化的曲線。

圖4 不同粒度分布的軟巖物料(煤炭)在不同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目隨破碎時(shí)間的變化曲線Fig.4 Number of breaking Bonding bonds changed with breaking time of soft rock material (coal) with different particle size distribution at different gear roll speeds

圖5 不同粒度分布硬巖物料(巖石)在不同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目隨破碎時(shí)間的變化曲線Fig.5 Number of breaking Bonding bonds changed with breaking time of hard rock material (rock) with different particle size distribution at different gear roll speeds

對(duì)比分析圖4、圖5可知：①物料大塊率越低，單位時(shí)間內(nèi)斷裂的Bonding鍵數(shù)目越均勻，即物料出料越均勻，物料大塊率增加時(shí)，單位時(shí)間內(nèi)斷裂的Bonding鍵數(shù)目波動(dòng)越大，這與破碎機(jī)咬合物料的角度、物料尺寸有關(guān)；②針對(duì)圖5中硬巖物料，當(dāng)物料硬度較大時(shí)，破碎時(shí)會(huì)出現(xiàn)突然破碎，物料伴有一定的“跳蹦”現(xiàn)象，所以圖5(a)、圖5(b)、圖5(c)中Bonding鍵斷裂數(shù)目曲線出現(xiàn)一定波動(dòng)，而軟巖物料易于破碎，齒輥轉(zhuǎn)速較低時(shí)不會(huì)出現(xiàn)所謂的“跳蹦”，但齒輥轉(zhuǎn)速增大時(shí)也會(huì)伴有該現(xiàn)象的發(fā)生；③對(duì)于同種粒度分布的物料，破碎相同量的硬巖物料所花費(fèi)的時(shí)間普遍大于破碎軟巖物料所需時(shí)間，說明物料硬度的增加將會(huì)降低破碎機(jī)的破碎效率，另外，軟巖物料破碎完成時(shí)Bonding鍵斷裂數(shù)目較硬巖物料Bonding鍵斷裂數(shù)目大，軟巖物料Bonding鍵斷裂數(shù)目約為硬巖物料Bonding鍵斷裂數(shù)目的2倍左右，而Bonding鍵是組合物料塊度所設(shè)定的黏結(jié)鍵，斷裂數(shù)目越多，說明破碎后的物料成塊率越低，相反，斷裂數(shù)目越少，破碎后的物料成塊率越高。所以，上述對(duì)比結(jié)果說明：硬巖物料破碎后出料較軟巖出料成塊率大，即硬巖物料被強(qiáng)制排出的幾率大于軟巖物料。

圖6為軟巖物料不同物料粒度分布在相同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目隨時(shí)間的變化關(guān)系，圖7為硬巖物料不同物料粒度分布在相同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目隨時(shí)間的變化關(guān)系。

由圖6、圖7可知：①圖6中軟巖物料大塊率由A組(10%)增大到B組(15%)時(shí)，在相同齒輥轉(zhuǎn)速下，曲線斜率相似，即破碎機(jī)破碎效率基本相同，但破碎過程中的Bonding鍵斷裂數(shù)目A組略大于B組，當(dāng)大塊率由B組(15%)增大到C組(20%)時(shí)，破碎時(shí)間較短時(shí)物料破碎效率相似，破碎時(shí)間較長(zhǎng)時(shí)C組物料破碎效率明顯下降，這說明軟巖物料大塊率在一定范圍內(nèi)時(shí)對(duì)破碎機(jī)破碎效率影響較小，當(dāng)大塊率超過該范圍時(shí)將導(dǎo)致破碎機(jī)的破碎效率降低；②比較圖7可知，對(duì)于硬巖物料，大塊率的增加會(huì)降低破碎機(jī)的破碎效率，且降低的幅度逐漸減少；③通過分析圖6、圖7中Bonding鍵斷裂數(shù)目可知，A組物料粒度分布下的Bonding鍵斷裂數(shù)目高于B組和C組物料粒度分布下的Bonding鍵斷裂數(shù)目，即物料粒度越低，破碎機(jī)破碎越徹底，強(qiáng)制排出越少，破碎效率越高，特別是隨著物料硬度的增加，這種現(xiàn)象表現(xiàn)得越明顯。

4 結(jié) 論

(1)物料性質(zhì)對(duì)破碎機(jī)的破碎效率影響較大。物料粒度分布相同，硬度不同時(shí)，Bonding鍵斷裂數(shù)目差距很大，硬度大時(shí)破碎相同量的物料耗時(shí)長(zhǎng)，破碎機(jī)破碎效率低,Bonding鍵斷裂數(shù)目少。

圖6 不同物料粒度分布軟巖物料在相同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目Fig.6 Number of Bonding bonds with different material particle sizes of soft rock material distributed at same gear speed

圖7 不同物料粒度分布硬巖物料在相同齒輥轉(zhuǎn)速下的Bonding鍵斷裂數(shù)目Fig.7 Number of Bonding bonds with different material particle sizes of hard rock material distributed at same gear speed

(2)物料大塊率的增加會(huì)降低破碎機(jī)的破碎效率。對(duì)于軟巖物料，大塊率在一定范圍內(nèi)時(shí)對(duì)破碎機(jī)破碎效率影響較小，當(dāng)大塊率超過該范圍時(shí)將導(dǎo)致破碎機(jī)的破碎效率降低。對(duì)于硬巖物料，大塊率的增加會(huì)降低破碎機(jī)的破碎效率，且與大塊率增加的范圍關(guān)系不大。因此，嚴(yán)格控制破碎機(jī)的入料粒度，特別是物料大塊率，有助于提高破碎機(jī)的破碎效率。

(3)對(duì)于軟巖物料，當(dāng)物料粒度分布一定時(shí)，適當(dāng)提高雙齒輥破碎機(jī)齒輥轉(zhuǎn)速能在一定程度上增大破碎機(jī)破碎效率，但破碎機(jī)破碎效率增加的速度逐步減小。對(duì)于硬巖物料，當(dāng)物料粒度分布一定時(shí)，齒輥轉(zhuǎn)速為31 r/min時(shí)破碎機(jī)的破碎效率略高于齒輥轉(zhuǎn)速為33 r/min和35 r/min時(shí)破碎機(jī)的破碎效率。齒輥轉(zhuǎn)速的變化還會(huì)造成破碎后物料的粒度發(fā)生改變：低齒輥轉(zhuǎn)速下會(huì)使破碎機(jī)破碎后物料粒度分布較為均勻，強(qiáng)制排出大塊較少；而高齒輥轉(zhuǎn)速易造成強(qiáng)制排出大塊物料，物料粒度分布不均。

[1] 羅 銘，王興容，郝其昌．固定式液壓破碎機(jī)的開發(fā)研制及其應(yīng)用[J]．金屬礦山，2007(10)：92-96．

Luo Ming，Wang Xingrong，Hao Qichang．Development of fixed type hydraulic crusher and its application in mines[J]．Metal Mine，2007(10)：92-96．

[2] 貢凱軍，何有鈞．MMD雙齒輥破碎機(jī)生產(chǎn)能力的計(jì)算方法[J]．金屬礦山，1997(5)：37-38.

Gong Kaijun，He Youjun．Efficiency calculation of MMD double tooth-roller crusher[J]．Metal Mine，1997(5)：37-38．

[3] 徐 泳，孫其誠(chéng)，張 凌，等．顆粒離散元法研究進(jìn)展[J]．力學(xué)進(jìn)展，2003，33(2)：251-260．

Xu Yong，Sun Qicheng，Zhang Ling，et al．Advance in discrete element methods for particulate materials[J]．Advances in Mechanics，2003，33(2)：251-260．

[4] 徐紅玉，許占民，艾爾肯·艾海提，等．BP-1100/700型顎式破碎機(jī)動(dòng)顎仿真優(yōu)化[J]．金屬礦山，2010，40(7)：116-118．

Xu Hongyu，Xu Zhanmin，Arken Ahat，et al．Simulation optimization analysis for BP-1100/700 jaw type crusher[J]．Metal Mine，2010，40(7)：116-118．

[5] 黃鵬鵬，肖觀發(fā)，李 成，等.基于EDEM的物料破碎效果仿真分析[J]．礦山機(jī)械，2014，41(10)：76-80．

Huang Pengpeng， Xiao Guanfa，Li Cheng，et al．Simulation and analysis on material crushing effects based on EDEM[J]．Mining & Processing Equipment,2014，41(10)：76-80．

[6] 崔 斌，郝秀娟，陳希紅，等．雙齒輥破碎機(jī)生產(chǎn)能力的計(jì)算與分析[J]．機(jī)械工程師，2010，41(5)：147-148．

Cui Bin，Hao Xiujuan，Chen Xihong，et al．Efficiency calculation and analysis of double tooth-roller crusher[J]．Mechanical Engineer，2010，41(5)：147-148．

[7] 梁克敏，吳曉紅，牛婷志．淺析新型雙齒輥破碎機(jī)的生產(chǎn)能力[J]．礦山機(jī)械，2009，37(15)：87-89．

Liang Kemin，Wu Xiaohong，Niu Tingzhi．Simple analysis on the production capacity of the new Bi-Toothed Roll Crusher[J]．Mining & Processing Equipment，2009，37(15)：87-89．

[8] 王國(guó)強(qiáng)，郝萬(wàn)軍，王繼新．離散單元法及其在EDEM上的實(shí)踐[M]．西安：西安工業(yè)大學(xué)出版社，2010．

Wang Guoqiang，Hao Wanjun，Wang Jixin．Discrete Element Method and Its Practice on EDEM[M]．Xi′an： Xi′an Technological University Press，2010．

[9] 胡國(guó)明．顆粒系統(tǒng)的離散元素法分析仿真[M]．武漢：武漢理工大學(xué)出版社，2010．

Hu Guoming．Analysis and Simulation of Granular System by Discrete Element Method Using EDEM[M]．Wuhan：Wuhan University of Technology Press，2010．

[10] 余世科．基于離散元法的多尺度內(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.