基于正交試驗(yàn)的制粒機(jī)成型仿真及結(jié)構(gòu)參數(shù)優(yōu)化

楊濤濤,周星宇,陸曉麗,錢(qián)善華,俞經(jīng)虎

(江南大學(xué)機(jī)械工程學(xué)院，江蘇省食品先進(jìn)制造裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室，江蘇 無(wú)錫 214122)

環(huán)模制粒機(jī)由于制粒過(guò)程穩(wěn)定、生產(chǎn)效率高、能耗低、成型率高等一系列優(yōu)點(diǎn)[1]成為飼料生產(chǎn)的主要機(jī)械之一。近年來(lái)，國(guó)內(nèi)外學(xué)者在環(huán)模制粒成型的功耗、顆粒成型質(zhì)量和機(jī)理、力學(xué)特性、模輥磨損以及振動(dòng)特性等方面進(jìn)行了大量的研究，為環(huán)模制粒機(jī)的結(jié)構(gòu)優(yōu)化和工藝參數(shù)最佳匹配提供了寶貴的參考。Rolfe等[2]研究顆粒尺寸大小、環(huán)模轉(zhuǎn)速和物料含水率三者對(duì)環(huán)模成型機(jī)物料成型的影響,結(jié)果表明，增大物料含水率、提高環(huán)模轉(zhuǎn)速、降低物料尺寸可以提高物料成型壓力和密度，但在成型過(guò)程也會(huì)增加產(chǎn)品溫度。沈江飛[3]建立制粒機(jī)參數(shù)化多體力學(xué)模型，并設(shè)計(jì)正交試驗(yàn)分析多個(gè)因素對(duì)機(jī)器振動(dòng)的影響規(guī)律，結(jié)果表明，增大主軸直徑和支撐架剛度或減小壓輥底座與主軸的連接剛度可有效降低轉(zhuǎn)軸振動(dòng)和擠壓力變化幅值，而轉(zhuǎn)軸軸承間距存在與軸承剛度相關(guān)的最優(yōu)值。Kaliyan等[4]針對(duì)秸稈物料的壓縮特性和致密過(guò)程物料的本構(gòu)模型兩方面進(jìn)行研究，建立了其彈塑性固體模型，通過(guò)對(duì)物料的壓縮試驗(yàn)，分析得出彈性模量和粘性系數(shù)與原料的耐磨性和抗壓強(qiáng)度有關(guān)。王詠梅等[5]應(yīng)用POLYFLOW軟件對(duì)環(huán)模制粒機(jī)擠壓區(qū)流場(chǎng)進(jìn)行模擬，研究不同工藝參數(shù)對(duì)流場(chǎng)壓力、速度、剪切速率和黏度分布特性的影響規(guī)律，發(fā)現(xiàn)喂料量為6 t·h-1、環(huán)模線(xiàn)速度為6.5 m·s-1和物料含水率為15%時(shí)所形成的流場(chǎng)有助于苜蓿草粉的制粒成型。李震等[6]采取流固耦合法對(duì)模輥關(guān)鍵部位進(jìn)行疲勞壽命分析，結(jié)果顯示，環(huán)模在靠近?？孜恢锰幾钊菀装l(fā)生疲勞破壞，其最小壽命為3.12×106次，并通過(guò)傳統(tǒng)應(yīng)力分析辦法驗(yàn)證耦合法的可行性及有效性。制粒機(jī)結(jié)構(gòu)參數(shù)的改變定然導(dǎo)致內(nèi)部物料流動(dòng)的變化，從而影響制粒機(jī)的成型效果。但是，由于物料的擠壓過(guò)程復(fù)雜，目前較少通過(guò)模擬制粒機(jī)成型數(shù)值來(lái)研究結(jié)構(gòu)參數(shù)對(duì)制粒效果的影響。因此，本文以谷物飼料為對(duì)象，應(yīng)用COMSOL軟件模擬分析物料在制粒機(jī)內(nèi)擠壓成型的流動(dòng)過(guò)程，并基于正交試驗(yàn)，研究制粒機(jī)的結(jié)構(gòu)參數(shù)對(duì)成型速率及能耗的影響，結(jié)合權(quán)矩陣綜合分析方法提出最佳的結(jié)構(gòu)參數(shù)組合。

1 原理與方法

1.1 環(huán)模制粒機(jī)成型原理

環(huán)模制粒機(jī)主要由喂料機(jī)、攪拌室和制粒機(jī)三部分組成。機(jī)器工作時(shí)，物料先是通過(guò)喂料機(jī)輸送到攪拌室內(nèi)，與室內(nèi)噴嘴灑出的調(diào)質(zhì)原料進(jìn)行混合攪拌，提高含水率、黏度等特性。當(dāng)完全攪拌后，物料便被送入制粒機(jī)，其內(nèi)部運(yùn)作截面如圖1所示。根據(jù)所受壓力不同，可將制粒機(jī)內(nèi)物料運(yùn)動(dòng)區(qū)域分為三個(gè)區(qū)：送料區(qū)、變形壓緊區(qū)、擠壓成型區(qū)。

注：1—壓輥；2—環(huán)模；3—切刀；4—擠壓成型區(qū)；5—變形壓緊區(qū)；6—送料區(qū)。Note:1—Pressure roller；2—Ring die；3—Cutter；4—Extrusion zone；5—Deformation compression zone；6—Feeding area.圖1 制粒機(jī)結(jié)構(gòu)Fig.1 Pellet mill structure

圖中，R為環(huán)模內(nèi)半徑，r為壓輥半徑，x為模輥間隙，α為擠壓成型區(qū)角度，β為變形壓緊區(qū)角度，根據(jù)圖1所示的幾何關(guān)系推導(dǎo)α、β數(shù)值。

(1)

(2)

成型過(guò)程中，環(huán)模以一定角速度順時(shí)針轉(zhuǎn)動(dòng)。物料首先進(jìn)入到送料區(qū)，在環(huán)模的帶動(dòng)下，慢慢地向變形壓緊區(qū)運(yùn)動(dòng)。在此期間，環(huán)模與壓輥的間距變小，物料不斷地受到擠壓，密度增大，與輥面產(chǎn)生的摩擦推動(dòng)壓輥同樣作順時(shí)針轉(zhuǎn)動(dòng)，但轉(zhuǎn)速比環(huán)模小。當(dāng)運(yùn)動(dòng)到變形壓緊區(qū)，物料繼續(xù)受壓擠，所受壓強(qiáng)隨弧長(zhǎng)呈二次曲線(xiàn)增長(zhǎng)[7](圖2)，少量物料也被擠出?？?。物料在擠壓成型區(qū)所受壓強(qiáng)基本為恒定值。在此區(qū)間，與壓輥接觸面的減少使得擠壓力急劇增大，同時(shí)顆粒間的接觸面積增大，物料產(chǎn)生較好地黏結(jié)，密度便穩(wěn)定不變。當(dāng)擠壓力能夠克服?？變?nèi)物料與孔壁摩擦?xí)r，物料便被擠出?？?，再通過(guò)切刀的切割，最終形成飼料顆粒。

圖2 物料所受壓強(qiáng)變化趨勢(shì)Fig.2 Change trend of pressure on materials

1.2 數(shù)值仿真

1.2.1數(shù)學(xué)模型 為了簡(jiǎn)化仿真計(jì)算，對(duì)數(shù)學(xué)模型作出以下假設(shè)：①物料運(yùn)動(dòng)為無(wú)彈性非牛頓流體流動(dòng)，且考慮壓縮性；②制粒機(jī)內(nèi)各截面物料分布均勻；③忽略黏性生熱及質(zhì)量力的影響。本文依據(jù)SZLH420制粒機(jī)主要技術(shù)參數(shù)(表1)在COMSOL中建立了制粒機(jī)中心截面的流域模型(圖3)，流域由環(huán)模的旋轉(zhuǎn)域和內(nèi)部的固定流域組成，中間空白圓形區(qū)域是壓輥所在位置。

表1 SZLH420制粒機(jī)主要技術(shù)參數(shù)Table 1 Main technical parameter of SZLH420 pellet mill

圖3 制粒機(jī)中心截面流域Fig.3 Mathematical model of central section of pellet mill

1.2.2仿真參數(shù)及邊界條件 流體流動(dòng)的控制方程主要包括質(zhì)量守恒方程、動(dòng)量守恒方程和描述流體黏性的流變學(xué)本構(gòu)方程[8]。在擠壓過(guò)程中，物料以塑性存在于流域中；當(dāng)剪切應(yīng)力超過(guò)臨界應(yīng)力時(shí)，則表現(xiàn)為流體的流動(dòng)特性。因此，本文采用HBP (Herschel-Bulkley-Papanastasiou)模型[9]對(duì)物料流變特性進(jìn)行表征,其模型及等效黏度系數(shù)表達(dá)式如下。

ταβ=τβα=μB(εαβ)n+τy(1-e-mγ)

(3)

(4)

式中，μB為流動(dòng)一致性系數(shù)，單位Pa·s；τy是臨界屈服應(yīng)力，單位N·m-2；m是模型參數(shù)，單位s，參數(shù)值越大，塑性存在時(shí)間越短；n是流動(dòng)特性指數(shù)，n<1為假塑性流體，n=1為牛頓流體，n>1為膨脹性流體。

考慮到制粒前的調(diào)質(zhì)攪拌過(guò)程，依據(jù)谷物飼料的初始屬性(表2)設(shè)置仿真中進(jìn)入制粒機(jī)前物料密度為845 kg·m-3，粘度流動(dòng)一致性系數(shù)μB為650 Pa·s，臨界屈服應(yīng)力τy圓整為6 N·m-2。在環(huán)模內(nèi)做擠壓流動(dòng)的物料可近似認(rèn)為是假塑性非牛頓流體且處于塑性時(shí)間較短，因此流動(dòng)特性指數(shù)n可設(shè)為0.95，模型參數(shù)m設(shè)為100。

表2 谷物飼料參數(shù)Table 2 Grain feed parameter

入口邊界設(shè)為質(zhì)量流，流率0.4 kg·s-1，出口邊界設(shè)為壓力出口，壓力為大氣壓，壓輥處圓形壁面設(shè)為動(dòng)網(wǎng)格，其余為無(wú)滑移壁面，摩擦系數(shù)為0.31。設(shè)置環(huán)模旋轉(zhuǎn)域轉(zhuǎn)速20 rad·s-1，并通過(guò)一致邊界實(shí)現(xiàn)旋轉(zhuǎn)域與固定域物料流動(dòng)的連續(xù)性。

1.3 正交試驗(yàn)

1.3.1正交試驗(yàn)設(shè)計(jì) 制粒機(jī)的結(jié)構(gòu)變化必然會(huì)造成成型速率及模輥能耗的變化。其中制粒機(jī)的環(huán)模內(nèi)徑、?？组L(zhǎng)度、壓輥直徑以及模輥間隙是影響成型效果及功耗的主要參數(shù)，依據(jù)SZLH420制粒機(jī)的原始參數(shù)設(shè)計(jì)了一組L9(43)正交試驗(yàn)，參數(shù)水平如表3所示。并運(yùn)用COMSOL軟件對(duì)各個(gè)試驗(yàn)方案進(jìn)行模擬仿真分析。

表3 正交試驗(yàn)參數(shù)Table 3 Pramater for orthogonal test

1.3.2試驗(yàn)指標(biāo) 成型速率和能耗是評(píng)價(jià)制粒機(jī)工作的重要指標(biāo)。通過(guò)仿真并不能直接反應(yīng)這兩個(gè)重要指標(biāo)，需要進(jìn)行數(shù)據(jù)轉(zhuǎn)換。本文將單位時(shí)間內(nèi)擠壓成型區(qū)及變形壓緊區(qū)的模孔平均產(chǎn)量Q(g·s-1)作為速率評(píng)價(jià)指標(biāo)，其表達(dá)式如下。

(5)

式中，vi為單個(gè)?？啄┒藦较蛩俾剩瑔挝?m·s-1；rm為?？装霃剑瑸?×10-3m；n為區(qū)域模孔個(gè)數(shù)；ρ為物料密度，為845 kg·m-3。

能耗(W)指標(biāo)主要由兩部分組成[10]，一是壓輥碾軋功耗(W1)，二是變形區(qū)、壓緊區(qū)內(nèi)物料受壓擠出?？椎钠骄龉?W2)，其表達(dá)式如下。

(6)

式中，P模輥間平均壓強(qiáng)，單位Pa，可從仿真中獲??；f為摩擦系數(shù)，為0.31；R為環(huán)模內(nèi)徑，單位m；B為模輥實(shí)際工作厚度，取0.1 m；α為擠壓成型區(qū)角度，單位rad；β為變形壓緊區(qū)角度，單位rad；Pi為單個(gè)?？鬃畲髩簭?qiáng)，單位Pa。

2 結(jié)果與分析

2.1 中心截面流域網(wǎng)格

依據(jù)物理場(chǎng)以及所設(shè)定的邊界條件對(duì)制粒機(jī)中心截面流域進(jìn)行較細(xì)化網(wǎng)格劃分，局部效果如圖4所示。模型邊界處采用自由四邊形劃分，其余區(qū)域則是自由三角形劃分，網(wǎng)格質(zhì)量良好，并沒(méi)有出現(xiàn)較大網(wǎng)格，可進(jìn)行下一步仿真分析。

圖4 局部網(wǎng)格劃分Fig.4 Local meshing

2.2 仿真結(jié)果分析

從圖5可以看出，速度場(chǎng)和溫度場(chǎng)分布基本呈中心對(duì)稱(chēng)。?？變?nèi)物料由于環(huán)模順時(shí)針轉(zhuǎn)動(dòng)具有較高的絕對(duì)速度，且離中心越遠(yuǎn)處速度越大，瞬時(shí)速度最高可達(dá)5.4 m·s-1。近環(huán)模及壓輥壁面物料在摩擦力的作用下產(chǎn)生較快的流動(dòng)。而模輥間物料運(yùn)動(dòng)速度普遍較慢，隨著模輥間隙的減小，速度慢慢降低，甚至出現(xiàn)回流及渦旋的現(xiàn)象。這是因?yàn)樵搮^(qū)域流動(dòng)速度相對(duì)較低，無(wú)法克服沿流動(dòng)方向正壓強(qiáng)梯度的影響而造成的。

圖5 仿真結(jié)果Fig.5 Simulation result

在壓力場(chǎng)中，主要分為正壓區(qū)及負(fù)壓區(qū)。隨著模輥間距減小，兩區(qū)域物料逐漸被壓縮，產(chǎn)生的反作用力愈來(lái)愈大，在接近最小間隙處正壓區(qū)達(dá)到79.1 MPa，負(fù)壓區(qū)達(dá)到78 MPa。而最小間隙處由于正負(fù)壓的抵消，壓強(qiáng)相對(duì)較低。模孔內(nèi)，物料只受壁面摩擦力作用不再受壓擠，壓強(qiáng)沿?？诇p小?？傮w來(lái)說(shuō)，仿真結(jié)果符合實(shí)際物料運(yùn)動(dòng)過(guò)程，能夠明顯反映不同區(qū)域的流動(dòng)狀態(tài)。

2.3 正交結(jié)果分析

基于正交試驗(yàn)的仿真結(jié)果如表4及圖6所示。單以產(chǎn)量為優(yōu)化目標(biāo)時(shí)，最佳方案為A2B2C1D1，即環(huán)模內(nèi)徑420 mm、模孔長(zhǎng)度54 mm、壓輥直徑178 mm、模輥間隙1.5 mm，且產(chǎn)量主次順序?yàn)镈(模輥間隙)>B(?？组L(zhǎng)度)>C(壓輥直徑)>A(環(huán)模內(nèi)徑)；若以能耗小為優(yōu)化目標(biāo)時(shí)，最佳方案為A3B1C3D3，即環(huán)模內(nèi)徑425 mm、?？组L(zhǎng)度52 mm、壓輥直徑182 mm、模輥間隙2.5 mm,且產(chǎn)量主次順序?yàn)镈(模輥間隙)>B(模孔長(zhǎng)度)>C(壓輥直徑)>A(環(huán)模內(nèi)徑)。

表4 正交試驗(yàn)結(jié)果Table 4 Result of orthogonal test

圖6 各因素對(duì)指標(biāo)影響Fig.6 Influence of various factors on indicators

2.4 權(quán)矩陣分析

2.4.1正交結(jié)果分析 從正交試驗(yàn)結(jié)果上看，產(chǎn)量高的方案明顯能耗比較大。因此，為了獲得產(chǎn)量高且能耗較少的綜合方案，本文采用權(quán)矩陣方法優(yōu)化選出最優(yōu)方案。根據(jù)權(quán)矩陣原理[11]，在計(jì)算產(chǎn)量的權(quán)矩陣時(shí)，指標(biāo)以較大值為優(yōu)化目標(biāo)；在計(jì)算能耗的權(quán)矩陣時(shí)，指標(biāo)以較小值為優(yōu)化目標(biāo)。最終計(jì)算出不同指標(biāo)的權(quán)矩陣ω1、ω2，并由權(quán)矩陣平均計(jì)算可得到正交試驗(yàn)的綜合評(píng)價(jià)的權(quán)矩陣ωa，其計(jì)算結(jié)果如下所示。

ω1=[0.012 9，0.012 9，0.012 9，0.062 0，0.062 8，0.061 2，0.052 4，0.051 7，0.051 2，0.216 6，0.204 6，0.198 8]T

(7)

ω2=[0.010 0，0.010 0，0.010 2，0.068 3，0.061 1，0.063 4，0.032 0，0.031 2，0.032 9，0.156 6，0.189 0，0.214 3]T×10-4

(8)

ωa=[0.006 5，0.006 5，0.006 5，0.031 0，0.031 4，0.030 6，0.026 2，0.025 9，0.025 6，0.108 3，0.102 3，0.099 4]T

(9)

從權(quán)矩陣中各因素水平的權(quán)重可知，因素A中A3的權(quán)重最大，因素B中B2的權(quán)重最大，因素C中C1的權(quán)重最大，因素D中D1的權(quán)重最大，因此以較高的產(chǎn)量、較低的能耗為優(yōu)化目標(biāo)的最佳參數(shù)為A3B2C1D1，即環(huán)模內(nèi)徑設(shè)置成425 mm，模孔長(zhǎng)度設(shè)置成54 mm，壓輥直徑設(shè)置成178 mm，模輥間隙設(shè)置成1.5 mm。

2.4.2優(yōu)化結(jié)果分析 再次運(yùn)用COMSOL對(duì)優(yōu)化后方案進(jìn)行模擬，優(yōu)化前后對(duì)比結(jié)果如表5及圖7所示。對(duì)比各項(xiàng)指標(biāo)得出，優(yōu)化前后擠壓區(qū)?？啄┒似骄俣葀提高了0.022 m·s-1，產(chǎn)量指標(biāo)增加了1.634 g·s-1，能耗指標(biāo)增加了1.591 kJ。由圖7可以看出，優(yōu)化前后距中軸線(xiàn)小角度內(nèi)(擠壓變形區(qū))的壓力差距較大，最高可達(dá)約15 MPa；當(dāng)角度大于15°時(shí)，兩者的壓力基本相同?？傮w來(lái)說(shuō)，優(yōu)化后的方案較之于優(yōu)化前在環(huán)模的內(nèi)表面壓力增大，增大的壓強(qiáng)使得模孔內(nèi)物料擠出流動(dòng)的速率加快，產(chǎn)量也相應(yīng)提高，但同時(shí)能耗也加大。和其他正交試驗(yàn)方案相比，優(yōu)化后的方案是在能耗增長(zhǎng)不太大的情況下，制粒機(jī)的產(chǎn)量得到一定提升，從而使得整體制粒效率進(jìn)一步提高。

圖7 正壓區(qū)環(huán)模內(nèi)表面壓力變化Fig.7 Pressure change of inner surface of ring die in positive pressure zone

表5 優(yōu)化前后對(duì)比Table 5 Comparison result before and after optimization

3 討論

目前，由于制粒機(jī)擠壓過(guò)程的復(fù)雜性以及物料性質(zhì)的特殊性，很少通過(guò)模擬制粒機(jī)內(nèi)部流動(dòng)過(guò)程來(lái)對(duì)其結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化，大部分都是以擠壓試驗(yàn)進(jìn)行研究。叢宏斌等[12]依據(jù)成型機(jī)理建立壓輥受力和物料成型能耗模型，繪制不同參數(shù)下壓輥載荷分布曲線(xiàn)和能耗曲線(xiàn)，并以此提出設(shè)備的技術(shù)參數(shù)優(yōu)化方法。但物料被壓縮制粒過(guò)程很難通過(guò)擠壓試驗(yàn)得出普遍性規(guī)律，這類(lèi)方法相對(duì)于模擬流動(dòng)仿真有一定的局限性且成本較高。因此，本文利用COMSOL仿真研究了制粒機(jī)環(huán)模內(nèi)徑、?？组L(zhǎng)度、壓輥直徑和模輥間隙對(duì)成型性能的影響，發(fā)現(xiàn)通過(guò)對(duì)制粒機(jī)內(nèi)部模型的仿真模擬，能夠有效地反映物料在實(shí)際工作時(shí)的運(yùn)動(dòng)狀態(tài)，對(duì)性能影響由大到小的結(jié)構(gòu)因素排序?yàn)?模輥間隙、?？组L(zhǎng)度、壓輥直徑、環(huán)模內(nèi)徑。采用權(quán)矩陣法確定優(yōu)化后的最佳參數(shù)組合為環(huán)模內(nèi)徑425 mm、?？组L(zhǎng)度54 mm、壓輥直徑178 mm、模輥間隙1.5 mm。相比于優(yōu)化前，能耗雖然增大，但成型速率增加，能夠有效地提高工作效率及性能。后期還可以改變COMSOL軟件中材料參數(shù)、邊界條件來(lái)研究不同物料和工藝參數(shù)對(duì)制粒機(jī)成型性能的影響，從而實(shí)現(xiàn)制粒機(jī)成型性能進(jìn)一步優(yōu)化。