，， ，

(1. 湖北工業(yè)大學(xué)機(jī)械工程學(xué)院，湖北武漢430068； 2. 湖北省地質(zhì)勘查裝備中心，湖北武漢430022)

在對(duì)巖石成分進(jìn)行光譜分析時(shí)，需要將采集的樣品粉碎研磨后，裝入密閉容腔中進(jìn)行混料處理，其目的是通過混合過程得到組成和性質(zhì)均勻的礦粉，以保證光譜分析的效果[1]?；旌闲Ч暮脡闹苯佑绊懼庾V分析的精度，所以對(duì)礦粉混合均勻性提出了更高的要求[2-3]。所謂混合，即礦粉微粒在外力的作用下發(fā)生運(yùn)動(dòng)速度和方向的改變，使粉體微粒均勻分布的操作過程[4-5]。由于混合過程中微粒的運(yùn)動(dòng)十分復(fù)雜，在如何提高微?；旌暇鶆蚨确矫嬉恢笔侨缃裱芯康臒狳c(diǎn)和難點(diǎn)，許多學(xué)者已進(jìn)行了相關(guān)研究。梁增華等[6]針對(duì)鈦粉在不同振動(dòng)方向下的振動(dòng)狀態(tài)，通過仿真分析了鈦粉填充密度的影響。任繼良等[7]研究了顆粒運(yùn)動(dòng)軌跡對(duì)火藥運(yùn)動(dòng)接觸的影響，通過仿真分析其混合均勻性。張紫薇等[8]模擬了典型轉(zhuǎn)速下三維滾筒內(nèi)的碳球顆?；旌线^程，分析研究滾筒內(nèi)顆粒流態(tài)隨時(shí)間的演化過程，展示不同顆粒流態(tài)下的顆粒混合結(jié)構(gòu)特性。當(dāng)前，在礦粉微粒的混合方法、混合均勻度及其混合效率等方面的關(guān)聯(lián)性還有待深入研究。

為了提高微粒的混合均勻度和效率，本文在前人研究的基礎(chǔ)上提出了一種混料方法，即在振動(dòng)混料的條件下，向混料容器中加入一種攪拌介質(zhì)，使其在混合過程中起到攪拌礦粉的作用，以獲得更好的混料效果。通過EDEM軟件建立微粒和攪拌介質(zhì)仿真模型，模擬微粒在混料容器中的運(yùn)動(dòng)狀態(tài)，將混合到一定程度的微粒進(jìn)行試樣分析，應(yīng)用標(biāo)準(zhǔn)偏差和離散度來判斷混料的效果，分析攪拌介質(zhì)對(duì)微?；旌暇鶆蚨鹊挠绊?，并通過振動(dòng)混料實(shí)驗(yàn)對(duì)仿真結(jié)果進(jìn)行驗(yàn)證。

1 固體介質(zhì)的選用

針對(duì)礦粉微粒的混合特性，對(duì)其施加一定的機(jī)械外力，能夠?qū)ξ⒘＿M(jìn)行有效的重新排列來提高混合均勻性，這種機(jī)械外力可以統(tǒng)稱為物料混合力，通過混合運(yùn)動(dòng)作用在礦粉微粒上[9-10]。

由于加入了攪拌介質(zhì)，在混合過程中，它與礦粉微粒之間相互摩擦碰撞會(huì)產(chǎn)生磨損現(xiàn)象，從而影響光譜分析的效果,因此，原則上應(yīng)選用耐磨且不與礦粉微粒發(fā)生反應(yīng)的介質(zhì)。

為了充分表達(dá)攪拌效果，該介質(zhì)應(yīng)懸浮在礦粉物料中，與礦粉一起受到物料混合力。當(dāng)采用振動(dòng)作為物料混合的驅(qū)動(dòng)力，所需的混料運(yùn)動(dòng)將受到振幅和頻率的影響，通過振動(dòng)參數(shù)定性分析可知，密度太大的介質(zhì)，會(huì)使振動(dòng)幅度有減小的趨勢(shì)，影響微粒在容腔中的運(yùn)動(dòng)；密度太小的介質(zhì)，與微粒碰撞效果不明顯，起不到攪拌的作用，因此，所選介質(zhì)的密度應(yīng)與實(shí)驗(yàn)礦粉的密度接近[11-12]。

同時(shí)還應(yīng)考慮攪拌介質(zhì)的幾何尺寸，尺寸太大的介質(zhì)因占用的空間大，會(huì)阻礙微?；旌系倪\(yùn)動(dòng)軌跡，微粒在運(yùn)動(dòng)中的相對(duì)位移就會(huì)減小，從而影響微粒混合的均勻性，因此，應(yīng)避免選用尺寸較大的介質(zhì)。

綜合來看，選用瑪瑙球作為攪拌介質(zhì)較為合適[13]。主要是因?yàn)楝旇У拿芏韧V粉的密度相近，且其主要成分是二氧化硅，即使在攪拌過程中發(fā)生元素滲透，也不影響光譜分析的效果。如實(shí)驗(yàn)需要將粒度為50～150 μm的礦料及粉狀催化劑均勻混合，宜選用半徑為2～3 mm的瑪瑙球作為攪拌介質(zhì)。實(shí)驗(yàn)時(shí)，可將2個(gè)坩堝密閉緊扣形成一個(gè)密閉容腔，模型如圖1所示，呈杯狀結(jié)構(gòu)，底面半徑為22 mm，高度為12 mm，將其密閉緊扣。把瑪瑙球和被混物料放入其中，混合均勻后將瑪瑙球取出以供實(shí)驗(yàn)用。

1—坩堝；2—瑪瑙球；3—粉狀光譜分析試樣。圖1 微粉混合示意圖Fig.1 Mixed schematic diagram

2 建立顆粒接觸模型

采用EDEM軟件中的Hertz-Mindlin無(wú)滑動(dòng)接觸模型，其結(jié)構(gòu)參數(shù)和物理意義如下：設(shè)半徑分別為R1、R2的2個(gè)球形顆粒發(fā)生彈性接觸，法向重疊的計(jì)算公式為

α=R1+R2-|r1-r2|，

(1)

式中：r1、r2是2個(gè)顆粒球心位置矢量。

顆粒間的接觸面為圓形，接觸半徑a為

(2)

式中,R*為等效粒子半徑，可由下式求出

(3)

顆粒間的法向力Fn可由下式求得：

(4)

式中,E*為等效彈性模量，由下式求出

(5)

式中：E1、υ1、E2、υ2分別為顆粒1和顆粒2的彈性模量和泊松比。

(6)

顆粒間切向力Ft可由下式求出：

Ft=-Stδ，

(7)

式中:δ為切向重疊量；St為切向剛度，由下式求出：

(8)

式中，G*為等效剪切模量，由下式求出：

(9)

式中G1和G2為2個(gè)顆粒的剪切模量。

顆粒間的切向阻尼力Ft可由下式求出：

(10)

3 幾何模型和混合評(píng)價(jià)方法

3.1 幾何模型

用Solidworks建立密閉坩堝三維模型導(dǎo)入到EDEM。由于礦粉微粒形狀的復(fù)雜性，還原真實(shí)形態(tài)的仿真比較困難，因此在EDEM中采用簡(jiǎn)化模型，用3種不同的顆粒模型分別代表礦粉、粉狀催化劑和瑪瑙球，如圖2、3所示。

3.2 微粒的生成方法

在Factory模塊下建立3個(gè)顆粒場(chǎng)，F(xiàn)actory1中選擇動(dòng)態(tài)生成方式將礦粉微粒無(wú)限生成，直至距離坩堝底面大約5 mm的位置，從0.1 s開始生成顆粒。Factory2的顆粒場(chǎng)平面建立在距底面9.5 mm的位置，在礦粉微粒的基礎(chǔ)上生成4 mm的粉狀催化劑。在Factory3中選擇static方式生成顆粒，在顆粒類型(type)選擇瑪瑙球。根據(jù)粉體微粒及振動(dòng)的相關(guān)條件，設(shè)置微粒和瑪瑙球的相關(guān)參數(shù)，如表1所示。

圖2 礦粉或催化劑模型(粒徑/75 μm)Fig.2 Mineral powder or catalyst model (particle size 75 μm)

圖3 瑪瑙球模型(半徑為2.5 mm)Fig.3 Agate ball model (radius 2.5 mm)

表1 礦粉、催化劑和瑪瑙球相關(guān)參數(shù)

3.3 3個(gè)指標(biāo)的評(píng)價(jià)

利用標(biāo)準(zhǔn)偏差、混合度、混合指數(shù)對(duì)混合程度的高低進(jìn)行評(píng)價(jià)[14]。在EDEM軟件中，將密閉坩堝容器進(jìn)行網(wǎng)絡(luò)空間格劃分，保證顆粒所占區(qū)域的網(wǎng)格類型近似相同。首先對(duì)坩堝建立4×4×4的空間格，導(dǎo)出各個(gè)時(shí)間段中每個(gè)空間格所占各種顆粒的數(shù)量，并計(jì)算出空間格中顆粒的總量，設(shè)xi為單元格中催化劑微粒所占的比例數(shù)，表示所有比例數(shù)的均值,本文選取催化劑為研究對(duì)象[15]。通過得到的數(shù)據(jù)對(duì)混合效果進(jìn)行分析。標(biāo)準(zhǔn)偏差也稱為均方差根，表示數(shù)據(jù)波動(dòng)的幅度，其計(jì)算方法是為

(11)

標(biāo)準(zhǔn)偏差S值越小，說明數(shù)據(jù)波動(dòng)的幅度越小，混合的效果就越好，但是，S值只與測(cè)定值和標(biāo)準(zhǔn)平均值的差值有關(guān)，而與各測(cè)定值的本身大小無(wú)關(guān)[16]。根據(jù)光譜分析實(shí)驗(yàn)的需要，在實(shí)驗(yàn)前抽樣實(shí)測(cè)礦粉混合的均勻程度對(duì)光譜分析精度的影響，礦粉混合均勻度對(duì)應(yīng)的S值應(yīng)為0.05～0.06較合適。

實(shí)際上，當(dāng)混合物的組分比例相差懸殊時(shí)，標(biāo)準(zhǔn)偏差就很難反映出混合效果了，還需引入離散度R來反應(yīng)混合均勻程度。離散度的公式為

(12)

由于本文中所述的混料物的粒徑為50～150 μm的礦料及粉狀催化劑的組成比較接近，并無(wú)太大懸殊，因此，本文中僅以標(biāo)準(zhǔn)偏差S值作為混合均勻度的評(píng)價(jià)指標(biāo)。

4 模擬結(jié)果分析與實(shí)際應(yīng)用

利用EDEM軟件模擬，對(duì)瑪瑙球?qū)ΦV粉微?；旌暇鶆蚨鹊挠绊戇M(jìn)行分析。

4.1 瑪瑙球數(shù)量對(duì)礦粉微粒均勻度的影響

在水平面上對(duì)密閉坩堝施加一個(gè)水平振動(dòng)的機(jī)械外力。 通過對(duì)振動(dòng)特性的研究，查閱參考相關(guān)文獻(xiàn)[17]，設(shè)定振動(dòng)頻率和振幅分別為2 kHz和0.5 mm時(shí)，坩堝混料較為合理。 在此基礎(chǔ)上，瑪瑙球的個(gè)數(shù)對(duì)坩堝內(nèi)礦粉微粒的均勻性影響如圖4所示。

圖4 瑪瑙球數(shù)量對(duì)微?；旌暇鶆蚨鹊挠绊慒ig.4 Effect of agate bulb number on particle mixing uniformity

由圖可知：1 min內(nèi)，標(biāo)準(zhǔn)偏差會(huì)出現(xiàn)上下起伏現(xiàn)象，這是由于微粒之間發(fā)生碰撞有一個(gè)過渡的過程，之后微粒均勻性逐漸清晰，標(biāo)準(zhǔn)偏差隨時(shí)間增加呈逐漸下降趨勢(shì)。

2 min時(shí)，存放1個(gè)瑪瑙球的坩堝仿真環(huán)境下，標(biāo)準(zhǔn)偏差為0.059 8，最先達(dá)到光譜分析的實(shí)驗(yàn)條件。不加瑪瑙球的條件下，其標(biāo)準(zhǔn)偏差在2.5 min后才達(dá)到實(shí)驗(yàn)要求，而加入2個(gè)瑪瑙球在3 min后才達(dá)到效果，相比不加瑪瑙球的條件，混合速率更低，這是由于多個(gè)瑪瑙球的體積過大，阻礙了微粒的運(yùn)動(dòng)軌跡，減少礦粉微粒的相對(duì)運(yùn)動(dòng)，礦粉微粒容易積成塊狀難以混合，即瑪瑙球的數(shù)量應(yīng)避免選用太多。

4 min之后，混合處于穩(wěn)定，微?；旌匣具_(dá)到均勻，這時(shí)瑪瑙球?qū)炝暇鶆蚨扔绊懖淮螅?，瑪瑙球有助于礦粉微粒的混合，且時(shí)間在2.5 min左右，坩堝中加入1個(gè)瑪瑙球較為優(yōu)化。

4.2 翻面傾角對(duì)礦粉混合均勻度的影響

通過改變坩堝與水平面的傾斜角度，沿坩堝平面方向施加一維振動(dòng)，并對(duì)坩堝施加一定的轉(zhuǎn)速，這樣可以把積壓在坩堝底部的礦粉微粒實(shí)現(xiàn)翻面運(yùn)動(dòng)，從而判斷其對(duì)混合均勻度的影響。

設(shè)定實(shí)驗(yàn)的振幅為0.5 mm，頻率為2 kHz，轉(zhuǎn)速為120 r/min,坩堝中放入1個(gè)瑪瑙球。不同的傾角下，礦粉微粒混合均勻度的影響如圖5所示。

圖5 翻面傾角對(duì)微粒混合均勻度的影響Fig.5 Effect of face dip angle on mixing uniformity of particles

從圖中可知：2 min時(shí)，傾角為25 °的坩堝標(biāo)準(zhǔn)偏差為0.059 2，最先達(dá)到試驗(yàn)要求,且2～3 min的混合效率明顯優(yōu)于其他狀態(tài)下的混合效率，這是由于瑪瑙球在上升過程中受到重力的作用而降低了上升幅度，在下降過程中受重力的作用加速下降，減少了微粒的振動(dòng)幅度，傾斜一定角度能夠使坩堝運(yùn)動(dòng)的速度方向與重力方向形成傾斜夾角，從而可以削弱重力的影響；當(dāng)傾角為10 °較小時(shí)，微?；旌闲Чc水平振動(dòng)基本相同，并未到達(dá)最佳混合效果；當(dāng)傾角為45 °較大時(shí)，微粒的混合效果差，雖然削弱了重力的影響，但坩堝主要由底部的支撐轉(zhuǎn)換為側(cè)壁的支撐，瑪瑙球的上升幅度會(huì)下降，微粒下降到一定的高度也無(wú)法下降，從而影響了微粒的混合均勻性。

所以，合理的偏轉(zhuǎn)角度能夠?qū)⒅亓ψ饔棉D(zhuǎn)換到側(cè)壁上面，使側(cè)壁提供少量支撐力，從而有利于瑪瑙球在坩堝中更好的攪拌混合微粒。綜上所述，坩堝與水平面傾角為25 °左右時(shí)，能達(dá)到最佳混合效果。

4.3 振動(dòng)混料試驗(yàn)

依據(jù)上述的模擬結(jié)果，制造了一臺(tái)礦粉混料裝置(如圖6所示)，混料盤的盤面與水平面成一定的傾角 (0～35 °)。

1—偏心振動(dòng)軸;2—控制面板;3—混料盤。圖6 振動(dòng)混料實(shí)驗(yàn)臺(tái)Fig.6 Vibration mixing test rig

在偏心軸的作用下實(shí)現(xiàn)沿盤面方向的一維振動(dòng)(頻率0～3 kHz、振幅0.3～0.8 mm)實(shí)現(xiàn)振動(dòng)混料?；炝媳P還可圍繞偏心軸作旋轉(zhuǎn)運(yùn)動(dòng)(轉(zhuǎn)速0～300 r/min)，配合上述傾角，可實(shí)現(xiàn)混合料的翻面，從而提高混合效率。以上參數(shù)可根據(jù)實(shí)驗(yàn)需要進(jìn)選取和調(diào)節(jié)。實(shí)驗(yàn)時(shí)，可將圖1中的坩堝(裝有需混合礦粉、催化劑和瑪瑙球)安放在混料盤的圓孔中，此裝置可存放多種實(shí)驗(yàn)礦粉的混料。啟動(dòng)機(jī)器可按設(shè)定參數(shù)實(shí)現(xiàn)自動(dòng)混料，物料混合均勻后將瑪瑙球取出以供實(shí)驗(yàn)用。

設(shè)定基本參數(shù)，以混料盤的安裝孔全部裝有坩堝后，且每個(gè)坩堝的礦粉都能快速混合均勻?yàn)榍疤?。?yōu)選的實(shí)驗(yàn)參數(shù)為振動(dòng)頻率2 kHz、振幅0.3 mm，圓周轉(zhuǎn)速120 r/min。

通過改變瑪瑙球的數(shù)量和調(diào)節(jié)混料盤的傾角，將混合到各個(gè)時(shí)間點(diǎn)的微粒各選取10組樣品，按照上述模擬空間格的方法，將樣品盡量抹平在篩孔尺寸為0.074 mm的標(biāo)準(zhǔn)篩中，使微粒竟可能落在標(biāo)準(zhǔn)篩的空隙中, 把模擬空間格平面化。在計(jì)算機(jī)成像平面中，由于礦粉與催化劑微粒灰度相差較大，可以利用像素格分析法對(duì)標(biāo)準(zhǔn)篩進(jìn)行拍照，使用計(jì)算機(jī)將照片分割成像素格再進(jìn)行對(duì)比。由此方法得到各個(gè)時(shí)間段10組樣品中催化劑微粒所占比例，通過公式(11)得到混合均勻度的標(biāo)準(zhǔn)偏差S值。

瑪瑙球個(gè)數(shù)和傾角對(duì)礦粉混合均勻度的影響如圖7—8所示。

圖7 瑪瑙球數(shù)量對(duì)礦粉混合均勻度影響的對(duì)比Fig.7 Comparison of influence of agate ball quantity on mixing uniformity of mineral powder

圖8 礦粉混合均勻度的樣品值對(duì)比下傾角翻面不同F(xiàn)ig.8 Comparison of sample values of mineral powder mixing uniformity under different turning angles

由圖可以清楚地看出，坩堝中放入1個(gè)瑪瑙球，且混料盤的水平傾角為25 °時(shí)，最先滿足實(shí)驗(yàn)標(biāo)準(zhǔn)條件(參照礦粉混合均勻度對(duì)應(yīng)的S值)，且時(shí)間為2.5 min左右時(shí)的混合效率最為明顯；4 min后，混合趨于穩(wěn)定，微粒均勻度基本達(dá)到飽和，且在實(shí)際操作時(shí)該裝置運(yùn)行平穩(wěn)，實(shí)驗(yàn)效果良好。

通過實(shí)驗(yàn)還發(fā)現(xiàn)，采用該振動(dòng)混料實(shí)驗(yàn)臺(tái)，適當(dāng)增加混料盤轉(zhuǎn)速，混料效果更好。

綜上所述：實(shí)驗(yàn)結(jié)果與虛擬仿真結(jié)果基本吻合，即瑪瑙球在坩堝中起攪拌的作用，能夠更好地混合礦粉微粒。

5 結(jié)論

1)瑪瑙球的加入能充分發(fā)揮混料攪拌作用，可以高效地實(shí)現(xiàn)礦粉及催化劑的混合，優(yōu)選的振動(dòng)時(shí)間在2.5 min時(shí)只需在坩堝中加入1個(gè)瑪瑙球。

2)合適的振動(dòng)混料翻面傾角，能更好地發(fā)揮瑪瑙球的攪拌作用，加快微粉混合均勻的速度、混合效率更高，優(yōu)選的混料傾角為20～30 °。

3)實(shí)驗(yàn)結(jié)果反映了仿真模型的正確性，可為大批量混料作業(yè)的生產(chǎn)應(yīng)用提供借鑒和參考。