離心式投飼機(jī)拋料盤對(duì)飼料破碎率的影響

張麗珍， 朱倓

（上海海洋大學(xué)工程學(xué)院，上海海洋可再生能源工程技術(shù)研究中心，上海 201306）

投喂飼料是水產(chǎn)養(yǎng)殖中一項(xiàng)十分重要的工作內(nèi)容，其完成度的好壞直接關(guān)系到所養(yǎng)殖生命體的生長(zhǎng)發(fā)育情況，從而影響到養(yǎng)殖戶的經(jīng)濟(jì)效益[1]。我國(guó)的水產(chǎn)養(yǎng)殖模式主要有淡水養(yǎng)殖、海水養(yǎng)殖和循環(huán)水養(yǎng)殖，其中淡水養(yǎng)殖占比較大，而池塘養(yǎng)殖則是淡水養(yǎng)殖中最重要的部分[2]。在池塘養(yǎng)殖的飼料投喂中，離心式投飼機(jī)由于擁有較大的拋撒范圍以及較好的穩(wěn)定性等優(yōu)點(diǎn)而獲得了廣泛的應(yīng)用[3-6]。但其普遍存在飼料破碎率較高的問(wèn)題，不僅會(huì)污染水體，影響生物體的攝食行為和生長(zhǎng)環(huán)境，而且會(huì)導(dǎo)致飼料的利用率下降、成本上升。因此，亟需對(duì)降低破碎率的影響因素進(jìn)行研究。

目前，對(duì)于離心式投飼機(jī)的研究大多集中于投餌量精準(zhǔn)控制、投飼均勻性等方面，對(duì)破碎率的研究較少。馬迪紅[7]將飼料顆粒用質(zhì)量點(diǎn)代替，通過(guò)Adams軟件進(jìn)行仿真，根據(jù)其動(dòng)能變化發(fā)現(xiàn)飼料顆粒碰撞后跳出拋料盤從而導(dǎo)致破碎率增大的情況，設(shè)計(jì)將拋料盤添加上蓋板防止飼料跳出。但是因?yàn)閭鹘y(tǒng)計(jì)算機(jī)輔助工程（computer aided engineering，CAE）軟件（如Adams等）無(wú)法準(zhǔn)確建立飼料顆粒的真實(shí)模型用以模擬破碎，存在一定局限性。為了探究拋料盤加速肋板彎折角度對(duì)破碎率的影響，胡慶松等[8]在仿真軟件中通過(guò)不斷改變肋板彎折角度，觀察運(yùn)動(dòng)軌跡并進(jìn)行數(shù)據(jù)分析得到破碎率最小情況下的最佳肋板彎折角度。但其研究單個(gè)顆粒的運(yùn)動(dòng)情況，忽略了真實(shí)拋料情況下大量顆粒的綜合影響，受限于飼料顆粒的計(jì)算數(shù)量，無(wú)法準(zhǔn)確地對(duì)影響因素進(jìn)行定量分析。離散元仿真（discrete element modeling，DEM）軟件EDEM具有高數(shù)量級(jí)顆粒碰撞與破碎的仿真模擬以及數(shù)值計(jì)算能力[9-11]，廣泛應(yīng)用于礦物、谷物等物料運(yùn)動(dòng)的研究中[12-18]，可以用來(lái)研究離心式投飼機(jī)飼料落到拋料盤，與拋料盤碰撞引起破碎的情況。

本文通過(guò)對(duì)飼料顆粒在拋料盤上的動(dòng)力學(xué)分析，確定影響破碎率的主要因素，利用EDEM軟件構(gòu)建單個(gè)飼料顆粒黏結(jié)模型，以飼料小顆粒間黏結(jié)鍵斷裂比率作為破碎率的衡量標(biāo)準(zhǔn)，進(jìn)行單因素仿真實(shí)驗(yàn)及三因素三水平正交實(shí)驗(yàn)，建立各項(xiàng)參數(shù)與飼料破碎率之間的二次回歸模型，得到各因素對(duì)破碎率的影響規(guī)律，為拋料盤的優(yōu)化設(shè)計(jì)提供參考。

1 材料與方法

1.1 離心式投飼機(jī)工作原理

離心式投飼機(jī)結(jié)構(gòu)如圖1所示，飼料由料箱經(jīng)螺旋輸送裝置，從落料口落到電機(jī)驅(qū)動(dòng)的拋料盤上被拋撒出去。拋料盤結(jié)構(gòu)如圖2所示，圖中R為肋板的曲率半徑，α是拋料盤工作時(shí)盤面與水平面之間的夾角，即拋料盤的傾角。

圖1 投飼機(jī)結(jié)構(gòu)Fig.1 Structure of feeding machine

圖2 拋料盤結(jié)構(gòu)Fig.2 Structure of centrifugal disc

1.2 飼料在拋料盤上的運(yùn)動(dòng)分析

飼料顆粒從落料口落到拋料盤上，隨著拋料盤的轉(zhuǎn)動(dòng)，沿著肋板MA加速運(yùn)動(dòng)，直至離開(kāi)拋料盤（圖3）。

圖3 飼料在盤上運(yùn)動(dòng)Fig.3 Feed movement on centrifugal disc

飼料顆粒在拋料盤上運(yùn)動(dòng)期間，受到自身重力G、盤底面作用力Fd、與盤面的摩擦力F1、與肋板側(cè)面的摩擦力F2、離心力Fce、科氏力Fv，受力分析如圖4所示。

圖4 顆粒受力分析Fig.4 Particle stress analysis

當(dāng)顆粒沿肋板由M運(yùn)動(dòng)到A（圖3）時(shí)所受到的盤面摩擦力計(jì)算如下。

式中，μ為飼料與拋料盤表面之間的動(dòng)摩擦系數(shù)；m為顆粒質(zhì)量，kg；α為盤面與水平面的夾角，rad；g為重力加速度，m·s-2。

顆粒所受離心力如公式（2）所示。

式中，l為MA圓弧長(zhǎng)度，mm；v為顆粒做圓周運(yùn)動(dòng)的切向速度，m·s-1；v1為顆粒在A點(diǎn)的切向速度，m·s-1；θ為v1與OA連線的夾角，rad；ω為拋料盤轉(zhuǎn)動(dòng)角速度，rad·s-1。

飼料顆粒在拋料盤側(cè)面所受到的摩擦力如公式（3）所示。

式中，故沿MA飼料顆粒運(yùn)動(dòng)方程如公式（4）所示。

安全檢查表法簡(jiǎn)便、易于掌握，是常用的分析風(fēng)險(xiǎn)的方法之一，目前也有打分的檢查表法?？蓪?shí)現(xiàn)半定量評(píng)價(jià)。該檢查方法的計(jì)分方法是根據(jù)實(shí)際檢查結(jié)果對(duì)安全檢查表所有的評(píng)價(jià)項(xiàng)目分別賦予“優(yōu)秀”“良好”“中等”“較差”等不同的定性等級(jí)，同時(shí)對(duì)相關(guān)等級(jí)賦予相應(yīng)的權(quán)重，對(duì)權(quán)重進(jìn)行累計(jì)求和，得到實(shí)際評(píng)價(jià)值。然查詢預(yù)先設(shè)定好的安全級(jí)別表，確定評(píng)價(jià)值在安全級(jí)別表中的位置，得到安全等級(jí)。檢查表的內(nèi)容在符合相關(guān)安全要求的前提下，還應(yīng)該簡(jiǎn)明扼要，切合實(shí)際，層次分明，重點(diǎn)突出，在實(shí)際操作過(guò)程中，應(yīng)依據(jù)以下三點(diǎn)制定安全檢查表[5]：

式中，α為盤面與水平面的夾角，rad。

飼料在拋料盤內(nèi)運(yùn)動(dòng)方程如公式（5）所示。

通過(guò)對(duì)飼料顆粒的運(yùn)動(dòng)分析可知，從單個(gè)顆粒的運(yùn)動(dòng)微分方程（5）來(lái)看，影響飼料破碎的主要因素有拋料盤轉(zhuǎn)速ω、拋料盤傾角α和速度分量夾角θ，而速度分量夾角由肋板曲率半徑R決定。從顆粒碰撞能量損耗的角度來(lái)看，飼料的破碎情況也與拋撒過(guò)程中碰撞總次數(shù)有關(guān)，即與肋板的個(gè)數(shù)有關(guān)。因此，以拋料盤轉(zhuǎn)速ω、拋料盤傾角α、速度分量夾角θ以及肋板個(gè)數(shù)作為影響因素進(jìn)行研究。

1.3 飼料黏結(jié)模型的建立

1.3.1 模型選擇 為了能夠更加準(zhǔn)確且真實(shí)地模擬破碎情況，需要建立單個(gè)飼料顆粒的黏結(jié)模型。選擇2號(hào)蝦飼料，顆粒為近似圓柱體，長(zhǎng)度約為3 mm，底面直徑約為1.5 mm。使用多個(gè)球體小顆粒及小顆粒之間產(chǎn)生的黏結(jié)鍵[19-20]填充組成單個(gè)圓柱體形狀的飼料顆粒，在確定球體小顆粒的物理半徑后，合理地設(shè)置接觸半徑可以讓黏結(jié)鍵同時(shí)承擔(dān)切向與法向應(yīng)力，當(dāng)顆粒黏結(jié)鍵受到超過(guò)臨界法向或切向應(yīng)力時(shí)，黏結(jié)鍵斷開(kāi)，表明顆粒飼料破碎[21]。模型如圖5所示。

圖5 小球體顆粒間的黏結(jié)模型Fig.5 Adhesion model between small spherical particles

1.3.2 飼料黏結(jié)參數(shù)的確定 單個(gè)圓柱體飼料模型需用球形小顆粒進(jìn)行填充，用單個(gè)直徑為0.2 mm的球形小顆粒進(jìn)行填充，填充個(gè)數(shù)N的計(jì)算公式如下。

式中，α為填充體積分?jǐn)?shù)，取0.56；Vr為單個(gè)圓柱形飼料顆粒體積，mm3；N為球形小顆粒個(gè)數(shù)；Vf為球形小顆粒體積，mm3；R為球形小顆粒半徑，為0.2 mm；r為圓柱形飼料顆粒底面半徑，為0.75 mm；h為圓柱形飼料顆粒高度，為3 mm。將參數(shù)值代入公式（6）得填充的球形小顆粒個(gè)數(shù)N為88個(gè)。

確定顆粒數(shù)后，以單個(gè)飼料圓柱體大小的容器作為顆粒工廠（particle factory）產(chǎn)生小顆粒進(jìn)行填充。球形小顆粒填充滿顆粒飼料后，利用bonding模型產(chǎn)生小顆粒之間的黏結(jié)鍵，黏結(jié)過(guò)程如圖6所示。

圖6 顆粒黏結(jié)過(guò)程Fig.6 Small particle bonding process

式中，ν為小顆粒材料泊松比，為0.3；E為小顆粒彈性模量，為9.1×107Pa；R為小顆粒半徑，取0.2 mm；B為經(jīng)驗(yàn)系數(shù)，取0.5[23]。得到單位面積法向剛度Kt為9.8×106N·m-3；單位面積切向剛度Kn為4.9×106N·m-3。

1.4 拋料盤性能仿真

1.4.1 仿真模型構(gòu)建與物性參數(shù)確定 將圖2所示的拋料盤幾何體模型用SolidWorks軟件進(jìn)行建模后導(dǎo)入EDEM，其直徑為200 mm。為了確定仿真所需的可靠接觸參數(shù)，采用提筒試驗(yàn)法[24-26]進(jìn)行參數(shù)標(biāo)定，如圖7所示。仿真所需物性參數(shù)見(jiàn)表1。

表1 物性參數(shù)Table.1 Material properties of particles

圖7 拋料盤仿真模型Fig.7 Simulation model of centrifugal disc

1.4.2 單因素試驗(yàn) 為了研究各個(gè)因素對(duì)飼料破碎率的影響規(guī)律，并為后續(xù)的正交試驗(yàn)提供合適的參數(shù)設(shè)定范圍，選擇3個(gè)因素進(jìn)行單因素試驗(yàn)，其他因素置于范圍中間值。工作中螺旋送料裝置在螺旋軸轉(zhuǎn)速為5 r·s-1的情況下落料口每秒流出顆粒體積為75 cm3（約為1 000顆飼料顆粒總體積），所以設(shè)置仿真條件參數(shù)為：流量1 000顆·s-1、落料方向?yàn)榇怪毕蛳?、落料口距離盤心垂直距離H=15 cm；水平距離h=5 cm、接觸參數(shù)設(shè)置等均保持不變。試驗(yàn)參數(shù)變化范圍：拋料盤轉(zhuǎn)速5、10、15、20、25、30，35 r·s-1；肋板個(gè)數(shù)3、4、5、6、7、8、9個(gè)；肋板曲率半徑60、70、80、90、100、110、120 mm；拋料盤傾角0°、5°、10°、15°、20°、25°，30°。以飼料在拋撒過(guò)程中斷裂的黏結(jié)鍵占比作為破碎率，拋料仿真模型如圖7所示。

1.4.3 正交回歸試驗(yàn) 為了明確拋料盤結(jié)構(gòu)參數(shù)對(duì)破碎率的影響規(guī)律，將單因素試驗(yàn)中破碎率達(dá)到最低時(shí)3個(gè)因素的參數(shù)組合作為中間水平，進(jìn)行三因素三水平正交回歸試驗(yàn)并進(jìn)行分析。仍以飼料在拋灑過(guò)程中斷裂黏結(jié)鍵占比作為破碎率。試驗(yàn)因素水平見(jiàn)表2。

1.5 驗(yàn)證試驗(yàn)

為了驗(yàn)證破碎率回歸方程優(yōu)化結(jié)果的可靠性，將優(yōu)化后的參數(shù)值加工生產(chǎn)出拋料盤以及普通拋料盤（直徑200 mm、無(wú)傾斜角度、肋板個(gè)數(shù)為4、肋板平直無(wú)曲率）分別安裝后進(jìn)行對(duì)比試驗(yàn)。根據(jù)水產(chǎn)行業(yè)相關(guān)試驗(yàn)標(biāo)準(zhǔn)[27]，在空曠平地上放置好可調(diào)速離心式投飼機(jī)(轉(zhuǎn)速在5～35 r·s-1可調(diào))，料箱中裝滿飼料顆粒，以固定電機(jī)轉(zhuǎn)速啟動(dòng)，待其工作狀態(tài)穩(wěn)定后，將其整個(gè)拋料過(guò)程分為前、中、后3個(gè)階段，分別用收納袋在出口處收集3個(gè)階段拋出的飼料顆粒，每個(gè)階段收集到的飼料重量不少于3 kg用網(wǎng)孔邊長(zhǎng)尺寸為顆粒直徑0.8倍的18目篩篩出破碎顆粒并稱重，按式（7）計(jì)算破碎率，重復(fù)3次取平均值。

式中，y為飼料破碎率；m1為篩出的破碎飼料質(zhì)量，kg；m2為收集到的飼料質(zhì)量，kg。

2 結(jié)果與分析

2.1 單因素試驗(yàn)結(jié)果分析

由圖8可知，拋料盤傾角、肋板個(gè)數(shù)、肋板曲率半徑、盤的轉(zhuǎn)速都會(huì)影響飼料在拋撒過(guò)程中的破碎率。轉(zhuǎn)速與破碎率成正比關(guān)系，因此盤的具體轉(zhuǎn)速在滿足拋灑范圍需求的情況下取轉(zhuǎn)速為20 r·s-1。在拋料盤傾角、肋板個(gè)數(shù)、肋板曲率半徑變化時(shí)，飼料破碎率變化均呈現(xiàn)出先變小再變大的規(guī)律，當(dāng)肋板個(gè)數(shù)為5、拋料盤傾角為15°、肋板曲率半徑為70 mm時(shí)，破碎率表現(xiàn)最小，分別為2.14%、2.18%、1.29%。

圖8 單因素試驗(yàn)結(jié)果Fig.8 Single factor test results

2.2 正交試驗(yàn)結(jié)果分析

正交試驗(yàn)結(jié)果如表3所示，破碎率范圍在1.53%～4.58%之間，說(shuō)明不同因素水平對(duì)破碎率影響較大，需要進(jìn)行優(yōu)化，確定最佳參數(shù)組合。

表3 三因素三水平正交試驗(yàn)結(jié)果Table 3 Results of three factor and three level orthogonal test

采用Design-expert軟件對(duì)正交試驗(yàn)結(jié)果進(jìn)行數(shù)據(jù)分析，二次多元回歸擬合得到方差分析結(jié)果如表4所示，模型極顯著（P＜0.001），影響破碎率因素由大到小分別為C、A2、BC、C2、B、A、B2、AB、AC，其中C2、B、A對(duì)破碎率影響顯著（0.001＜P＜0.05）；B2、AB、AC對(duì)破碎率影響不顯著（P＞0.05）；C、A2、BC對(duì)破碎率的影響呈現(xiàn)極顯著（P＜0.001）。由此可見(jiàn)，結(jié)構(gòu)參數(shù)對(duì)于破碎率的影響由大到小依次為：曲率半徑＞肋板個(gè)數(shù)＞拋料盤傾角。已知各個(gè)系數(shù)間不存在線性相關(guān)性，剔除掉不顯著因素后，優(yōu)化得到破碎率y與各個(gè)影響因素之間的回歸方程為如下。

表4 方差分析結(jié)果Table 4 Variance analysis

2.3 響應(yīng)面分析

根據(jù)表4的方差分析結(jié)果，在交互項(xiàng)中BC對(duì)于破碎率影響呈現(xiàn)極顯著（P＜0.001），對(duì)其影響破碎率的情況需要進(jìn)行分析。肋板個(gè)數(shù)(B)與曲率半徑（C）交互作用對(duì)破碎率的影響如圖9所示，當(dāng)肋板個(gè)數(shù)不變時(shí)，破碎率隨著曲率半徑的增大呈現(xiàn)先減小再增大的趨勢(shì)。當(dāng)曲率半徑不變時(shí)，破碎率隨著肋板個(gè)數(shù)的增大呈現(xiàn)先減小后增大的趨勢(shì)。當(dāng)曲率半徑在65～75 mm，肋板個(gè)數(shù)在4.5～5.5個(gè)時(shí)，破碎率較小。肋板曲率半徑對(duì)于破碎率的影響相較于肋板個(gè)數(shù)更為明顯，說(shuō)明肋板曲率半徑是這一組交互作用中的更為重要的影響因素。

圖9 因素交互作用曲面Fig.9 Response surface curve

2.4 驗(yàn)證結(jié)果分析

為了得到最優(yōu)工作參數(shù)組合，根據(jù)回歸方程式（8），將破碎率回歸模型[28]作為目標(biāo)函數(shù)，各因素范圍作為約束條件，求解最小值后得到最優(yōu)參數(shù)組合為：拋料盤傾角12.76°（取12.8°），肋板個(gè)數(shù)5.80個(gè)（取6個(gè)），肋板曲率半徑為71.49 mm(取71.5 mm)，此時(shí)破碎率為2.13%。將所得到的參數(shù)組合用于仿真驗(yàn)證，得到破碎率為2.29%，與優(yōu)化結(jié)果基本吻合。由圖10可知，在安裝了改進(jìn)后的拋料盤情況下，離心式投飼機(jī)在各轉(zhuǎn)速飼料破碎率均低于相同條件下未改進(jìn)的離心式投飼機(jī)，平均破碎率降低38.8%。且在高轉(zhuǎn)速（＞15 r·s-1）情況下飼料破碎率降低得更為顯著。

圖10 改進(jìn)前后破碎率對(duì)比Fig.10 Comparison of crushing rate before and after improvement

3 討論

目前，對(duì)離心式投飼機(jī)的飼料破碎率尚未展開(kāi)深入研究，但在其他類型的投飼設(shè)備的飼料破碎研究領(lǐng)域有著較多研究成果。如邵愷懌等[29]設(shè)計(jì)氣動(dòng)投飼裝置的關(guān)鍵機(jī)械結(jié)構(gòu)，通過(guò)試驗(yàn)調(diào)試在一定程度上減少了飼料的破碎率。王志勇等[30]設(shè)計(jì)了網(wǎng)箱自動(dòng)投飼平臺(tái)，通過(guò)自動(dòng)控制系統(tǒng)設(shè)計(jì)試驗(yàn)方案，得出飼料顆粒粒徑、管道布置、輸送速率對(duì)飼料破碎率的影響。本研究采用仿真模擬與實(shí)物驗(yàn)證試驗(yàn)相結(jié)合的方法確定了影響因素與飼料顆粒破碎率之間的二次回歸模型，獲得拋料盤最優(yōu)工作參數(shù)組合：肋板個(gè)數(shù)為6個(gè)、拋料盤傾角為12.8°、肋板曲率半徑為71.5 mm。仿真結(jié)果與實(shí)物試驗(yàn)值較為吻合，優(yōu)化后破碎率降低明顯，可以作為今后離心式投飼機(jī)拋料機(jī)構(gòu)優(yōu)化的參考。

本研究針對(duì)飼料顆粒破碎仿真中所存在的問(wèn)題，采用離散元仿真方法對(duì)單個(gè)飼料顆粒內(nèi)部進(jìn)行模型建立，以小顆粒間黏結(jié)鍵的斷裂表征飼料破碎，解決了以往只能通過(guò)間接獲得飼料顆粒破碎數(shù)據(jù)的弊端，一定程度上避免了試驗(yàn)的偶然性。此外，本研究中所涉及到的影響因素暫時(shí)僅包含拋料盤的結(jié)構(gòu)參數(shù)，后期通過(guò)增加試驗(yàn)因素，并同時(shí)將均勻性、投飼效率等其他離心式投飼機(jī)關(guān)鍵性能指標(biāo)納入考量，從而使得優(yōu)化進(jìn)一步貼合生產(chǎn)使用實(shí)際。