王立宗，范吉軍，余南輝，張藝凡，劉澤軍，張永林

（武漢輕工大學，武漢 430023）

0 引言

碾米機主要利用碾白輥高速旋轉(zhuǎn)去除糙米糠層獲得白米，傳統(tǒng)臥式碾米機在碾米過程中產(chǎn)生的米糠碎米的混合料在下吸風口被收集［1-3］。為此，需要專門的糠粞分離器加以分離［4-8］，同時若下吸風量過大，會使大米沉淀于米篩底部，造成大米碾白不均，碎米率增大。國內(nèi)一些企業(yè)提出采用上吸風式碾米機，可實現(xiàn)對米糠的高效收集，并杜絕了米粒存積于下部的缺陷［9-11］。目前的研究為碾米機的米糠碎米分離問題提供了新的解決思路，但針對上吸風式碾米機風道設計及試驗鮮有報道。

近年來，隨著計算流體動力學（Computational Fluid Dynamics，CFD）技術(shù)的發(fā)展，已成為對農(nóng)業(yè)流體機械進行流場分析的重要手段［12-15］。張雙等［16］利用Fluent軟件對MNML30型立式砂輥碾米機內(nèi)風道進行有限元分析，研究了碾白輥轉(zhuǎn)速對碾米機內(nèi)氣流場的影響；李祖吉等［17］通過EDEM—Fluent耦合對傳統(tǒng)臥式碾米機進行氣固兩相流耦合仿真求解，獲得了內(nèi)部顆粒的運動狀態(tài)以及碾米機內(nèi)流場分布；江濤等［18］利用Fluent研究了聯(lián)合收割機風篩選裝置吹風口結(jié)構(gòu)對流場分布規(guī)律的影響，通過對比仿真結(jié)果得出最有利于清選的結(jié)構(gòu)分布；劉帥等［19］運用Fluent中DPM模型對花生清選室流場進行模擬，得出花生運動軌跡。

運用空氣動力學知識，對傳統(tǒng)碾米機內(nèi)風道進行改進。設計一種采用獨特環(huán)形分離風道的上吸風式碾米機，以實現(xiàn)在碾米機內(nèi)部直接對米糠碎米進行高效分離收集；利用Fluent軟件對碾米機風道內(nèi)流場進行仿真分析；制作試驗樣機，對分離效果進行試驗，驗證設計的有效性，對碾米工藝的改進具有重要意義。

1 結(jié)構(gòu)與工作原理

如圖1、圖2所示，一種帶有環(huán)形風道的上吸風碾米機，主要由碾白系統(tǒng)和上吸風式排糠系統(tǒng)組成。碾白系統(tǒng)主要包括：螺旋推進器、空心軸、碾白輥、米篩網(wǎng)、拔料輥等，上吸風排糠系統(tǒng)主要包括：環(huán)形風道、噴風管、氣壓隔板、上吸風口、下進風口等。

圖1 主軸系統(tǒng)結(jié)構(gòu)圖Fig.1 Structural diagram of the main shaft system

上吸風機置于碾米機頂部，米篩網(wǎng)與外殼之間的空間形成環(huán)形風道，環(huán)形風道下部為下進風口，下進風口處設置有氣壓隔板，將風道與外界大氣隔開，使環(huán)形風道下部維持負壓狀態(tài)。同時，氣壓隔板使下進風口上端變窄，形成狹管效應，上升風力更加集中，防止米糠從下進風口排出。氣壓隔板同時還具有折流板的作用，使上升氣流形成環(huán)形繞流。噴風機向空心軸內(nèi)噴風，噴風管連通空心軸和碾白輥的噴風孔，氣流直接從空心軸經(jīng)噴風管流入碾白室，減少壓力損失（見圖2）。

圖2 主軸系統(tǒng)橫截面圖Fig.2 Section view of the main shaft system

2 模型建立及仿真參數(shù)

2.1 簡化建模

針對環(huán)形風道內(nèi)的氣體流場進行分析，建立簡化模型，對噴風槽附近區(qū)域進行分析，圖3中箭頭示意氣流流向。對機械部件建立三維模型，并對流體域模型進行網(wǎng)格劃分。

圖3 仿真模型Fig.3 Simulation model

2.2 物理模型及參數(shù)

設定噴風機風速為10 m/s，上吸風口風速為6 m/s，下進風口為壓力入口，與大氣相通，氣壓為0 Pa。

依據(jù)雷諾數(shù)Re判定碾米機內(nèi)氣流流動狀態(tài)

式中 d——特征長度，m；

v——流速，m/s；

ρ——流體密度，kg/m3；

μ——動力黏度系數(shù)，Pa·s。

空心軸內(nèi)徑d為80 mm，進風口風速為10 m/s，空氣密度ρ為1.177 kg/m3，空氣動力黏度系數(shù)μ為1.8×10-5Pa·s。得空心軸處雷諾數(shù)Re=52 311>4 000，流場為湍流狀態(tài)，故選用k-ε湍流模型進行有限元分析。

許多研究中均采用多孔介質(zhì)模型模擬糧食堆內(nèi)氣流流動問題［20-22］，故采用此方法模擬碾白室內(nèi)糙米對氣流產(chǎn)生的阻力。氣流穿過多孔介質(zhì)區(qū)內(nèi)部孔隙會產(chǎn)生壓力損失，損失系數(shù)采用Ergun方程項進行計算：

式中 1/α——黏性阻力系數(shù)，m-2；

ε——孔隙率，%；

ds——顆粒等效球直徑，m；

ρs——顆粒密度，kg/m3；

C2——慣性阻力系數(shù)，m-1。

依據(jù)曾勇的研究可得碾白室內(nèi)物料的孔隙率約為0.5，糙米的等效球直徑約為3 mm，帶入式中計算出黏性阻力系數(shù)為3.3×107m-2，慣性阻力系數(shù)為 4 666 m-1。

2.3 顆粒懸浮速度

風選機械一般依據(jù)目標顆粒懸浮速度進行設計。顆粒的自由懸浮速度公式如下所示：

式中 v0——懸浮速度，m/s；

g——重力加速度，m/s2；

ρs——顆粒密度，kg/m3；

C——阻力系數(shù)。

碎米密度約1 600 kg/m3，米糠密實密度約為1 500 kg/m3；空氣密度 1.17 kg/m3。

由于阻力系數(shù)C是顆粒雷諾數(shù)Rep的函數(shù)，采用分區(qū)懸浮速度公式，根據(jù)待測物料的等效球直徑來確定阻力系數(shù)：

碎米顆粒ds約為1.2 mm，計算出碎米懸浮速度約為 6.4 m/s［23-25］；米糠顆粒 ds為 0.22 mm［26-27］，懸浮速度約為1.6 m/s，遠小于碎米。

3 結(jié)果與分析

對噴風管轉(zhuǎn)至水平位置以及傾斜45°的情形進行仿真分析。圖4為速度云圖，可以看出噴風管處流速明顯加快、風力集中，可提高米糠排出效率。環(huán)形風道內(nèi)流速分布均勻，無停滯區(qū)。氣壓隔板間空隙處的上升氣流流速從1.4 m/s上升至3～4 m/s（見表1），狹管效應明顯，有效防止米糠下落。

圖4 速度云圖Fig.4 Velocity nephogram

如圖5所示，由于上吸風機的引風作用，不論噴風管轉(zhuǎn)至水平或斜45°位置，環(huán)形風道內(nèi)氣流均保持向上流動，可以將米糠粉塵高效地向上運輸。氣壓隔板產(chǎn)生折流作用，上升氣流改變流向，形成環(huán)形繞流。

圖5 氣流流線Fig.5 Air streamline

對目標區(qū)域內(nèi)流速進行等距取樣測量，得到具有代表性的平均流速。由表1可知，噴風管旋轉(zhuǎn)位置對流場影響較??；噴風槽附近的米篩網(wǎng)網(wǎng)眼風速極快，達到21～25 m/s；由于物料的阻礙作用，遠離噴風槽處的米篩網(wǎng)網(wǎng)眼氣流接近于0 m/s；環(huán)形風道內(nèi)平均風速為4.1～4.6 m/s。圖6為兩種情形的垂直X方向風速分量等值面，僅噴風槽和上吸風口附近垂直方向流速分量大于6.5 m/s，環(huán)形風道內(nèi)流速大于2 m/s而小于6.5 m/s，介于米糠與碎米懸浮速度之間，可實現(xiàn)顆粒分離；碾白室內(nèi)噴風槽處平均風速為10.5～11.7 m/s，較高的流速易于大米流態(tài)化，能提高碾白效果。

圖6 X軸方向流速分量等值面Fig.6 Isosurface of flow velocity component in X-axis direction

表1 風速對照表Tab.1 Comparison of air velocity 單位：m·s-1

圖7為氣壓隔板處壓力云圖，噴風槽處的米篩網(wǎng)內(nèi)側(cè)的氣壓高達5 500 Pa，環(huán)形風道內(nèi)為負壓，約為-50 Pa，米篩網(wǎng)內(nèi)外壓差最大約5 550 Pa，遠大于常規(guī)碾米機噴風槽處氣壓差最低要求［28］，較高的內(nèi)外壓差可以提高米糠碎米排出速率，并防止米糠從下進風口排出。

圖7 壓力云圖Fig.7 Pressure nephogram

4 樣機試驗

為檢驗該新型碾米機對米糠碎米的分離效果，按1：2比例搭建樣機。選擇合適的風機使上吸風口和噴風機出口風速與仿真一致，為區(qū)分米糠與碎米，將碎米染成藍色，利用高速相機拍攝顆粒運動過程，并利用軟件PIVlab測量顆粒速度?？梢杂^察到：絕大部分米糠從上吸風口排出，僅氣壓隔板上沉積少量米糠，碎米均從下進風口排出，圖8為利用PIVlab測量得到的環(huán)形風道處米糠以及下進風口處碎米的速度矢量圖，米糠整體向上運動，最大上升速度為1.15 m/s；碎米向下掉落，最大掉落速度為0.25 m/s，分離效果基本與仿真預測結(jié)果一致，滿足實際生產(chǎn)需求。

圖8 風道內(nèi)顆粒速度矢量圖Fig.8 Vector of particle velocity in duct

5 結(jié)論

（1）運用Fluent軟件對上吸風式環(huán)形風道碾米機的流場進行分析。結(jié)果表明：環(huán)形風道內(nèi)氣流速度為4.1～4.6 m/s，垂直方向分速度為2.0～6.5 m/s，均介于米糠與碎米懸浮速度之間，沒有產(chǎn)生較大旋渦區(qū)域，有利于米糠從上吸風口快速排出，碾白輥的轉(zhuǎn)動對流場影響較小，工作狀況穩(wěn)定。

（2）仿真結(jié)果證實氣壓隔板形成狹管效應，下進風口處上升氣流由1.4 m/s加速至3～4 m/s，同時具有折流效應，有助于形成環(huán)形繞流。噴風槽處米篩網(wǎng)內(nèi)外壓差最大約5 550 Pa，空心軸內(nèi)的氣流直接從噴風管噴出，減少了風能損失。

（3）制作了試驗樣機，測得米糠最大上升速度為1.15 m/s，碎米最大下落速度為0.25 m/s，可實現(xiàn)米糠與碎米的有效分離。