宗 路，郭 晉，蔡改貧，肖賢煌，王 俊

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

磨礦作業(yè)是選礦流程中一道十分重要的工序，半自磨機(jī)具有磨礦流程短、生產(chǎn)率高、成本低的優(yōu)勢(shì)，在各類大型礦山中廣泛應(yīng)用，在實(shí)際生產(chǎn)中半自磨機(jī)存在著能耗太大的問(wèn)題。合理地設(shè)計(jì)磨機(jī)參數(shù)能夠有效地實(shí)現(xiàn)節(jié)能降耗，對(duì)半自磨機(jī)的能耗研究一直是國(guó)內(nèi)外學(xué)者研究的重要方向之一[1-4]。張永等[5]通過(guò)對(duì)半自磨機(jī)筒體襯板改進(jìn)及排礦格子板改造，提高了磨礦效率和襯板壽命。唐新民[6]用物理學(xué)方法對(duì)筒內(nèi)鋼球和礦石的運(yùn)動(dòng)規(guī)律展開(kāi)分析，得出高出襯板160、80 mm的壓條、鋼球大小和不同的充填率，與筒體襯板和鋼球的碎裂磨損和處理量、能耗的變化關(guān)系。楊樹(shù)新[7]、劉華[8]通過(guò)分析半自磨機(jī)襯板受力情況和襯板材質(zhì)方面指導(dǎo)設(shè)計(jì)襯板結(jié)構(gòu)。田秋娟等[9]、王繼生等[10]結(jié)合半自磨機(jī)設(shè)計(jì)參數(shù)和磨礦效果進(jìn)行分析研究。Toor[11]從半自磨機(jī)襯板的磨損影響磨礦性能的角度進(jìn)行量化統(tǒng)計(jì)，指出了襯板結(jié)構(gòu)變化對(duì)半自磨機(jī)整體性能的影響關(guān)系。賈曉帥等[12]分析了Mo、Si、Cu、Mn等不同合金元素及其含量對(duì)半自磨機(jī)襯板材料性能的影響，優(yōu)化后的襯板能夠獲得更好的力學(xué)性能。文獻(xiàn)[13]以筒體內(nèi)的單個(gè)介質(zhì)為研究對(duì)象，通過(guò)該介質(zhì)所受的合力分析襯板高度變化的影響。Djordjevic等[14]采用PFC軟件分析了半自磨機(jī)襯板結(jié)構(gòu)和轉(zhuǎn)速率對(duì)磨礦功耗的影響關(guān)系，為半自磨機(jī)的功耗研究方法提供了指導(dǎo)依據(jù)。Cleary[15]使用EDEM軟件模擬了半自磨機(jī)中顆粒的運(yùn)動(dòng)形態(tài)，研究了不同轉(zhuǎn)速率下筒內(nèi)物料的運(yùn)動(dòng)狀態(tài)、筒體轉(zhuǎn)矩和磨礦能耗，分析了相關(guān)參數(shù)變化對(duì)半自磨機(jī)磨礦能耗的影響。Collinao等[16]利用DEM軟件結(jié)合三維激光測(cè)量設(shè)備分析了半自磨機(jī)襯板的變形和磨損情況。

學(xué)者通過(guò)數(shù)理推導(dǎo)、樣機(jī)試驗(yàn)、模擬仿真在半自磨機(jī)設(shè)備研究上取得了一定的成果，但是側(cè)重于半自磨機(jī)工作參數(shù)的優(yōu)化，缺乏科學(xué)嚴(yán)謹(jǐn)?shù)陌胱阅C(jī)功率理論體系。本文中運(yùn)用動(dòng)力學(xué)分析建立以襯板結(jié)構(gòu)參數(shù)為變量的半自磨機(jī)破碎能耗數(shù)理模型，根據(jù)模型分析半自磨機(jī)襯板結(jié)構(gòu)參數(shù)變化對(duì)磨礦性能的影響趨勢(shì)和最優(yōu)參數(shù)范圍，利用仿真模擬和磨礦試驗(yàn)驗(yàn)證模型的可靠性，為改善磨礦效果、減少磨礦能耗提供一定的理論指導(dǎo)。

1 半自磨機(jī)工作原理

本文中所用半自磨機(jī)模型參考了某礦山φ10.37 m×5.19 m大型半自磨機(jī)的外型參數(shù)，并利用三維建模軟件建立該樣機(jī)的物理模型，如圖1所示。

圖1 半自磨機(jī)仿真物理模型Fig.1 Simulation physical model of semi-autogenous mill

半自磨機(jī)的工作原理是：磨機(jī)在電機(jī)的作用下以一定的轉(zhuǎn)速運(yùn)轉(zhuǎn)，利用襯板表面及帶有凹凸?fàn)钚问降奶嵘龡l襯板帶動(dòng)磨介和物料一同做周向運(yùn)動(dòng)，被提升的磨介和物料在達(dá)到拋落點(diǎn)后被拋出，通過(guò)這種周而復(fù)始的運(yùn)動(dòng)使得磨介與物料、物料與物料、物料與襯板之間相互沖擊來(lái)達(dá)到物料破碎的目的。在破碎過(guò)程中，物料存在明顯的分區(qū)現(xiàn)象，將半自磨機(jī)內(nèi)部物料分布區(qū)域分為提升區(qū)、死區(qū)（腎形蠕動(dòng)區(qū)）、瀉落區(qū)和拋落區(qū)4個(gè)區(qū)域，如圖2所示。

圖2 筒內(nèi)載荷分布區(qū)域劃分示意圖Fig.2 Sketch map of partition of load distribution in cylinder

2 混合運(yùn)動(dòng)狀態(tài)下有用功分析

2.1 提升區(qū)載荷有用功分析

設(shè)襯板軸向截面的結(jié)構(gòu)函數(shù)表達(dá)式為f（x），襯板的軸向截面面積S為

則單個(gè)襯板槽內(nèi)空間的軸向截面面積和體積為

式中：l為半自磨機(jī)筒體長(zhǎng)度，mm；N為半自磨機(jī)內(nèi)提升條的數(shù)量；R0為半自磨筒體內(nèi)徑，mm；R1為襯板頂端至筒心距離，mm。

在筒體旋轉(zhuǎn)一周的時(shí)間間隔內(nèi)，統(tǒng)計(jì)有N個(gè)槽體裝入鋼球送入提升區(qū)，則旋轉(zhuǎn)一周提升的鋼球總數(shù)量為

式中：P為鋼球自由堆放的孔隙率，%；D為鋼球直徑，mm。

當(dāng)鋼球磨介做拋物線運(yùn)動(dòng)下落到落回點(diǎn)時(shí)，其有用功為

式中：G為單個(gè)鋼球質(zhì)量，kg；v為鋼球脫離拋落點(diǎn)初速度，m/s；a為常系數(shù)，當(dāng)轉(zhuǎn)速n給定時(shí)

將式（6）代入式（5）中，得到區(qū)域內(nèi)旋轉(zhuǎn)一周時(shí)間內(nèi)球介所作功

則，完全提示區(qū)內(nèi)載荷的有用功率P1為

式中：R1為襯板頂端至筒心距離，R1=R0-H。

2.2 拋落區(qū)載荷有用功分析

磨介在向上運(yùn)動(dòng)過(guò)程中通常會(huì)伴隨有滑動(dòng)現(xiàn)象，故分析此球?qū)咏橘|(zhì)的運(yùn)動(dòng)速度v時(shí)必須引入與物料相關(guān)的滑動(dòng)摩擦系數(shù)即速度

將速度v代入式（5），得到做拋落運(yùn)動(dòng)的不完全提升區(qū)內(nèi)介質(zhì)運(yùn)動(dòng)的有用功為

在該區(qū)域內(nèi)設(shè)一半徑為R的圓軌跡，在軌跡上劃分出一厚度為dR的無(wú)限薄介質(zhì)層。該介質(zhì)層在半自磨機(jī)回轉(zhuǎn)體旋轉(zhuǎn)一周的時(shí)間范圍內(nèi)整體的質(zhì)量為

則該介質(zhì)層作用的有用功為

在內(nèi)層級(jí)半徑R2和外層級(jí)半徑R1之間積分上式，就可以得到在磨機(jī)轉(zhuǎn)一周時(shí)，整個(gè)落下球介所作之功

則做拋落運(yùn)動(dòng)的不完全提升區(qū)內(nèi)載荷的有用功率P2為

式中：R1為襯板頂端至筒心距離，R1=R0-H；R2為拋落

提升區(qū)載荷顆粒最小運(yùn)動(dòng)半徑

2.3 瀉落區(qū)載荷有用功分析

該區(qū)域載荷的重量G為

式中：φ'為R2半徑內(nèi)載荷的填充率，%；R2為瀉落提升區(qū)最大半徑，m；δ為載荷的堆密度，t/m3。

載荷重心點(diǎn)S位于摩擦角θ的夾角線上，其到圓心位置的距離為

式中Ω為瀉落區(qū)載荷的橫斷面對(duì)應(yīng)的圓心角，rad。根據(jù)力矩公式可得，瀉落區(qū)載荷對(duì)磨礦中心的力矩M為

磨礦過(guò)程中半自磨機(jī)需要由電動(dòng)機(jī)產(chǎn)生大小相同、方向相反的轉(zhuǎn)矩來(lái)克服力矩M對(duì)瀉落區(qū)域載荷作功，所以瀉落區(qū)內(nèi)載荷的有用功至少為

將上面得到的G值及X值代入公式，則

3 半自磨機(jī)破碎能耗模型的建立與分析

3.1 破碎能耗模型的建立

通過(guò)對(duì)混合運(yùn)動(dòng)狀態(tài)下的載荷運(yùn)動(dòng)形態(tài)劃分及有用功率分析，建立了各個(gè)區(qū)域內(nèi)載荷的磨礦有用功率，以此將各區(qū)域載荷有用功率進(jìn)行整合，即可得到半自磨機(jī)破碎能耗模型，即

式中：為半自磨機(jī)內(nèi)提升條高度，m；L為半自磨機(jī)內(nèi)提升條寬度，m；N為半自磨機(jī)內(nèi)提升條的數(shù)量，m。

3.2 襯板結(jié)構(gòu)參數(shù)對(duì)有用功率的影響分析

選取矩形襯板為研究對(duì)象，將試驗(yàn)測(cè)取的相關(guān)物理數(shù)據(jù)及選取的結(jié)構(gòu)參數(shù)代入式（19）中，分別得到以提升條高度、寬度和數(shù)量為變量的有用功率變化曲線。

圖3 有用功率變化曲線圖Fig.3 Useful power change curve

由圖3a可知，在提升條數(shù)量為0～32的范圍內(nèi)，半自磨機(jī)的有用功率隨著提升條數(shù)量的增加而急劇上升，但隨著提升條數(shù)量的增加，超過(guò)32個(gè)時(shí)，有用功率上升趨勢(shì)變緩。分析其原因是提升條數(shù)量的增加有利于提高筒內(nèi)載荷的拋出高度，增強(qiáng)了磨介對(duì)物料之間的自磨作用，提升了半自磨機(jī)的有用功率；但是由于半自磨機(jī)筒內(nèi)空間的限制，增加提升條數(shù)量則減少了載荷的破碎空間、襯板的槽內(nèi)空間和載荷運(yùn)動(dòng)的旋轉(zhuǎn)半徑，因此，隨著提升條的數(shù)量增加，提升載荷的增長(zhǎng)效果緩慢。分析圖3b可知，由在提升條高度為0～10 mm的范圍內(nèi)，半自磨機(jī)磨礦有用功率隨提升條高度的增加而急劇上升，但隨著提升條高度的進(jìn)一步增加，曲線出現(xiàn)拐點(diǎn)，在12～17 mm的高度范圍內(nèi)，有用功率近似達(dá)到最高值，高度高于17 mm之后功率開(kāi)始波動(dòng)下滑，由此判斷提升條高度為17 mm時(shí)有用功率為最優(yōu)值。由圖3c可知，隨著提升條高度的增加，有用功率呈現(xiàn)先增加后減小的趨勢(shì)，在20 mm左右處達(dá)到最大值。原因是當(dāng)寬度設(shè)計(jì)過(guò)大時(shí)，提升襯板之間的槽寬縮小，槽內(nèi)容量降低，可提升的載荷數(shù)量減少，有用功率隨之減?。划?dāng)襯板寬度過(guò)窄時(shí)，襯板之間的槽寬加大，槽內(nèi)載荷在提升區(qū)內(nèi)的滑動(dòng)和自傳加劇，難以達(dá)到理論拋落點(diǎn)拋出，導(dǎo)致磨礦過(guò)程中的有用功率減小。

4 仿真及實(shí)驗(yàn)

4.1 方案設(shè)計(jì)

將建立的半自磨機(jī)實(shí)驗(yàn)樣機(jī)模型導(dǎo)入EDEM軟件中，并設(shè)置入料物、襯板和磨介鋼球的材料屬性，不同材料之間的接觸參數(shù)，半自磨機(jī)筒體的轉(zhuǎn)速等仿真參數(shù)，設(shè)計(jì)3因素4水平的正交實(shí)驗(yàn)方案如表1所示，其中填充率為40%，轉(zhuǎn)速為48 r/min，磨礦時(shí)間為10 min。

表1 正交實(shí)驗(yàn)方案Tab.1 Orthogonal experiment scheme

4.2 結(jié)果

為了更好地量化描述半自磨機(jī)在不同襯板結(jié)構(gòu)參數(shù)設(shè)置下的磨礦性能，對(duì)各組仿真實(shí)驗(yàn)中地粘結(jié)鍵斷裂速率展開(kāi)分析，將各組試驗(yàn)中獲得的隨時(shí)間變化的粘結(jié)鍵斷裂值導(dǎo)出，繪制成曲線圖，不同提升條數(shù)量下的斷裂鍵數(shù)變化曲線如圖4所示。

圖4 不同提升條數(shù)量下的斷裂鍵數(shù)變化曲線Fig.4 Variation curves of fracture bond number under different lifting bars

從圖4中可以看出，提升條襯板條數(shù)的增加使得粘結(jié)鍵的斷裂速率逐漸增加。數(shù)量為32的曲線呈現(xiàn)出最大的變化斜率；提升條襯板數(shù)量從16增加到32的時(shí)候，粘結(jié)鍵的斷裂速度增幅趨勢(shì)相對(duì)減緩，變化幅度減弱，這是因?yàn)殡S著提升條襯板數(shù)量的增多襯板間的間距就相應(yīng)的減少，槽內(nèi)容量減小，導(dǎo)致被提升的載荷總量減少，影響了對(duì)載荷顆粒的提升效果，降低了半自磨機(jī)的破碎效率。

圖5為不同高度下斷裂鍵數(shù)變化曲線。從圖中可以明顯看出，提升條高度為17 mm的襯板，在仿真過(guò)程中粘結(jié)鍵的斷裂速度是相對(duì)較快的，高度超過(guò)17 mm的時(shí)候，粘結(jié)鍵的斷裂速率反而會(huì)隨著襯板高度的增加而下降，表明半自磨機(jī)提升襯板的高度設(shè)計(jì)有一個(gè)最佳范圍，超出范圍則磨機(jī)的磨礦性能會(huì)逐漸降低。

圖5 不同高度下斷裂鍵數(shù)變化曲線Fig.5 Variation curves of fracture bond number at different heights

圖6 不同寬度下斷裂鍵數(shù)變化曲線Fig.6 Variation curves of fracture bond number under different widths

圖6為不同寬度下斷裂鍵數(shù)變化曲線。從圖中可以明顯看出，寬度為20 mm的提升襯板，在仿真過(guò)程中粘結(jié)鍵的斷裂速度是最快的，當(dāng)寬度達(dá)到30、40 mm時(shí)，粘結(jié)鍵的斷裂速率反而會(huì)隨著提升襯板寬度的增加而下降，表明半自磨機(jī)提升襯板的寬度也有一個(gè)最佳范圍，當(dāng)大于或小于這一范圍時(shí)其磨礦性能開(kāi)始下降。

4.3 分析

課題組參照某大型礦山半自磨機(jī)回轉(zhuǎn)體主體結(jié)構(gòu)參數(shù)，縮小比例制作的半自磨機(jī)試驗(yàn)樣機(jī)如圖7所示，筒體有效內(nèi)徑和有效長(zhǎng)度分別為480、260 mm，并制作了不同規(guī)格的試驗(yàn)機(jī)提升襯板，利用DT96-E型的智能單相電能表采集半自磨機(jī)的磨礦電功耗數(shù)據(jù)，為方便觀察半自磨機(jī)內(nèi)部載荷顆粒的運(yùn)動(dòng)情況，將筒體一側(cè)端蓋替換為透明的亞克力材質(zhì)端蓋。以破碎能耗比為磨礦性能的評(píng)價(jià)指標(biāo)，即半自磨機(jī)在新生成單位質(zhì)量的新顆粒所消耗的能量，單位為kW/kg，按表1的試驗(yàn)方案進(jìn)行磨礦實(shí)驗(yàn)。

為了便于對(duì)試驗(yàn)結(jié)果的分析，對(duì)每大組磨礦后的物料進(jìn)行篩分，粒徑直徑區(qū)間為≥9、＜9～6、＜6～5、＜4～1的5個(gè)等級(jí)，并用電子稱稱出各個(gè)粒級(jí)的質(zhì)量和新生成顆粒的質(zhì)量（g）。將統(tǒng)計(jì)的各組破碎能耗比數(shù)據(jù)導(dǎo)入Design-Expert軟件中，分別繪制出提升條高度、數(shù)量、寬度對(duì)破碎比能耗產(chǎn)生影響的響應(yīng)曲面圖和等高線圖，分別如圖8、9所示。

圖7 半自磨機(jī)實(shí)驗(yàn)樣機(jī)Fig.7 Experimental prototype of semi-autogenous mill

圖8 提升條高度與提升條數(shù)量對(duì)破碎比能耗的響應(yīng)曲面圖和等高線圖Fig.8 Response surface and contour map of the height and number of lifting bars on energy consumption of crushing ratio

圖9 提升條高度與提升條寬度對(duì)破碎比能耗的響應(yīng)曲面圖和等高線圖Fig.9 Response surface and contour map of strip height and strip width to energy consumption of crushing ratio

從圖8a可知，當(dāng)提升條高度一定時(shí)，隨著提升條數(shù)量的增加，破碎比能耗逐漸降低，反映出比能耗值與提升條數(shù)量呈反比趨勢(shì)；當(dāng)提升條數(shù)量一定時(shí)，提升條高度從2 mm開(kāi)始增加，破碎比能耗呈現(xiàn)先減小后增大的趨勢(shì)，存在最低點(diǎn)。從圖8b可以發(fā)現(xiàn)，破碎比能耗的最低點(diǎn)出現(xiàn)于提升條高度為14～18 mm，提升條數(shù)量近似為32，故在半自磨機(jī)提升條襯板的寬度一定時(shí)，將提升條襯板的高度和數(shù)量參數(shù)分別設(shè)計(jì)為14～18 mm和32個(gè)的組合時(shí)，預(yù)計(jì)可以達(dá)到最佳的破碎能耗比。

圖9為提升條高度與提升條寬度對(duì)破碎比能耗的響應(yīng)曲面圖和等高線圖。從圖9a可知，當(dāng)給提升條高度一定時(shí)，隨著提升條寬度的增大，破碎比能耗呈現(xiàn)先減小后增大的趨勢(shì)；當(dāng)提升條寬度一定時(shí)，隨著提升條高度增加，破碎比能耗同樣也呈現(xiàn)先減小后增大的趨勢(shì)，所以必然存在一個(gè)最佳高度值和寬度值組合范圍。從圖9b可以看出，破碎比能耗的最低點(diǎn)出現(xiàn)于提升條高度為14～18 mm，提升條寬度為20～24 mm的組合范圍中，故在半自磨機(jī)提升條襯板的數(shù)量一定時(shí)，將提升條襯板的高度和寬度參數(shù)分別設(shè)計(jì)為14～18 mm和20～24 mm的組合時(shí)，預(yù)計(jì)可以達(dá)到最佳的破碎能耗比。

5 結(jié)論

1）根據(jù)半自磨機(jī)內(nèi)部載荷分布區(qū)域的劃分范圍，結(jié)合結(jié)構(gòu)參數(shù)和內(nèi)部載荷的動(dòng)力學(xué)分析，以區(qū)域載荷的運(yùn)動(dòng)入手累計(jì)求解整體載荷動(dòng)能，推導(dǎo)出能夠描述襯板結(jié)構(gòu)參數(shù)影響半自磨機(jī)磨礦有用功率的數(shù)學(xué)模型，分析了提升條襯板數(shù)量、高度、寬度對(duì)有用功率的影響規(guī)律。

2）基于離散元法對(duì)半自磨機(jī)內(nèi)部載荷顆粒的運(yùn)動(dòng)形態(tài)開(kāi)展可視化仿真實(shí)驗(yàn)，分析內(nèi)部載荷的運(yùn)動(dòng)狀態(tài)和粘結(jié)顆粒（物料）的破碎情況，得出半自磨機(jī)襯板結(jié)構(gòu)參數(shù)變化對(duì)磨礦性能的影響趨勢(shì)及最優(yōu)值：隨著半自磨機(jī)襯板數(shù)量增加，半自磨機(jī)磨礦效果越佳，隨著半自磨機(jī)襯板高度、寬度的增加，半自磨機(jī)磨礦效果曲線、載荷顆粒的拋落數(shù)量同呈先增大后減小的趨勢(shì)；當(dāng)半自磨機(jī)襯板數(shù)量為32，高度為17 mm，寬度為20 mm時(shí)，拋落的載荷顆粒數(shù)量相對(duì)最多，死區(qū)范圍相對(duì)最小，磨礦性能相對(duì)最佳，驗(yàn)證了半自磨機(jī)有用功率數(shù)學(xué)模型的準(zhǔn)確性。

參考文獻(xiàn)（References）：

[1]楊琳琳，文書(shū)明.自磨機(jī)和半自磨機(jī)的發(fā)展和應(yīng)用[J].國(guó)外金屬礦選礦，2013，4（7）：13-16.

[2]WEERASEKARA N S，POWELL L S，CLEARY P W，et al.The contribution of DEM to the science of comminution[J].Powder Technology，2013，248（248）：3-24.

[3]張光烈.我國(guó)自磨/半自磨技術(shù)新進(jìn)展[J].礦冶裝備，2011，42（3）：46-53.

[4]劉洪均，王彪，趙留成，等.自磨機(jī)/半自磨機(jī)功耗測(cè)試技術(shù)及其應(yīng)用[J].有色金屬（選礦部分），2013，35（2）：40-43.

[5]張永，孫東亮，韓寧剛.大型半自磨機(jī)襯板及格子板改造[J].有色金屬（選礦部分），2014，12（5）：75-76.

[6]唐新民.冬瓜山銅礦半自磨機(jī)筒體襯板和鋼球的研究[J].化工礦物與加工，2013，10（5）：20-25.

[7]楊樹(shù)新.自磨機(jī)筒體襯板對(duì)自磨工藝的影響[J].有色礦冶，2009，25（1）:50-52.

[8]劉華.φ5.5 m×1.8 m半自磨機(jī)襯板的改進(jìn)[J].礦山機(jī)械，2012，40（3）：133-135.

[9]田秋娟，郝萬(wàn)軍，鄧立營(yíng)，等.襯板的設(shè)計(jì)參數(shù)對(duì)球磨機(jī)磨礦效果的影響[J].礦山機(jī)械，2010，26（9）：73-76.

[10]王繼生，盧建坤，馬偉，等.襯板形狀對(duì)大型球磨機(jī)磨礦效率的影響[J].礦山機(jī)械，2013，41（8）：68-72.

[11]TOOR P.Quantifying the influence of liner wear on SAG mill performance[J].Sag Milling，2013，13（6）：2389-2408.

[12]賈曉帥，張志偉，黃雨，等.大型半自磨機(jī)襯板成分優(yōu)化模擬分析[J].金屬熱處理，2017，42（5）：1-8.

[13]MALEKIMM， YAHYAEIM， BANISI S.Converting AG to SAG mills:the Gol-E-Gohar iron ore company case[J].Powder Technology，2012，217（2）：100-106.

[14]DJORDJEVIC F，SHI N，MORRISON R.Determination of lifter design，speed and filling effects in AG mills by 3D DEM[J].Minerals Engineering，2004，39（17）：1135-1142.

[15]CLEARY P W.Predictinge hargemotion，powerdraw，segregation and wear inball mills using diserete element methods[J].Minerals Engineering，1998，11（11）：1061-1080.

[16]COLLINAO E，DAVILA P，IRARRAZABAL R，et al.Continuous improvement in SAG mill liner design using new technologies[J].XXVII International Mineral Processing Congress，2014，68（10）：2061-2080.