謝海軍,李志偉,張 靜
(華南農(nóng)業(yè)大學(xué) 工程學(xué)院,廣州 510642)
小粒種子是指平均直徑小于Φ3mm的種子[1],主要包括蔬菜、油菜、苜蓿、谷子、芝麻、花卉、林業(yè)和包衣煙草等。目前,針對小粒徑種子播種主要有育苗移栽和直播[2],我國超2/3蔬菜栽培采用穴盤育苗移栽,工廠化穴盤育苗是實現(xiàn)種苗商品化供應(yīng)主要方式。精量播種機是工廠化穴盤育苗的核心設(shè)備,前期研究表明,排種器轉(zhuǎn)速、吸孔結(jié)構(gòu)及正負氣壓對播種機的空穴率、單粒率、多粒率3項指標有顯著影響[3-5]。近年來,國內(nèi)外學(xué)者對氣力排種器的研究,有通過試驗的手段尋求各因素水平間的最佳組合,也有利用CDF軟件對排種器氣流場進行仿真[6-13]。本文運用ANSYS/FLUNET軟件氣力滾筒式排種器氣流場動力學(xué)模型,在吸附種子的條件下,仿真分析了吸孔結(jié)構(gòu)分別為沉孔、錐孔和直孔對排種器吸種效果的影響。
在薄板厚度相同的條件下,設(shè)置有沉孔吸孔、錐孔吸孔和直孔吸孔。利用ANSYS軟件FLUNET模塊建立吸孔氣流場區(qū)域的CFD結(jié)構(gòu)模型如圖1所示。圖1中,1為氣流入口邊界,2為氣流出口邊界,3為吸孔,4為種子,其余邊界為固定壁面邊界。采用自由網(wǎng)格劃分方式,單元截面形狀為三角形,單位長度間隔個數(shù)為100,在吸孔附近對網(wǎng)格進行進一步細化,更有利于計算結(jié)果。
圖1 吸種孔結(jié)構(gòu)問題區(qū)域
ANSYS/FLUENT是一個用于分析二維和三維流體流動場常用工具[14],由于在排種器上的吸孔中的氣流場是軸對稱流動,所以氣流場分析可以簡化為二維軸對稱流體流動問題進行求解。仿真分析時對模型進行如下假設(shè):①流場氣體為不可壓縮氣體;②室溫恒為25℃;③吸孔入口壓力恒為1.01×105Pa,出口壓力為0Pa;④排種器吸孔進氣口氣流速度分布是均勻;⑤在所有壁面上施加無滑移邊界條件。
仿真試驗方法為在吸孔直徑分別為Φ0.6mm、Φ0.8mm和Φ1.0mm分別以負壓為-3.0kPa對孔口結(jié)構(gòu)為沉孔、錐孔、直孔的吸孔進行氣流場仿真,得到3種孔口結(jié)構(gòu)下的吸孔入口中心氣流速度及種子表面平均氣壓,探討沉孔、錐孔、直孔的吸種和攜種效果。
負壓為-3.0kPa時,孔徑為Φ0.6mm的沉孔、錐孔和直孔氣流場的速度云圖和吸孔入口處壓力云圖如圖2所示。
(a) 沉孔速度云圖 (b) 沉孔壓力云圖
(c) 錐孔速度云圖 (d) 錐孔壓力云圖
(e) 直孔速度云圖 (f) 直孔壓力云圖
由圖2可知:負壓為-3.0kPa、吸孔直徑為Φ0.6mm時,沉孔型吸孔入口氣流平均速度為22.79m/s,被吸附種子表面平均氣壓為-0.844kPa;錐孔型吸孔入口氣流平均速度為18.83m/s,被吸附種子表面平均氣壓為-1.976kPa;沉孔型吸孔入口氣流平均速度為36.40m/s,被吸附種子表面平均氣壓為-0.406kPa。
負壓為-3.0kPa時,孔徑為Φ0.8mm的沉孔、錐孔和直孔氣流場的速度云圖和吸孔入口處壓力云圖如圖3所示。
由圖3可知:負壓為-3.0kPa、吸孔直徑為Φ0.8mm時,沉孔型吸孔入口氣流平均速度為24.84m/s,被吸附種子表面平均氣壓為-0.819kPa;錐孔型吸孔入口氣流平均速度為22.65m/s,被吸附種子表面平均氣壓為-2.651kPa;沉孔型吸孔入口氣流平均速度為33.09m/s,被吸附種子表面平均氣壓為-0.469kPa;
負壓為-3.0kPa時,孔徑為Φ1.0mm的沉孔、錐孔和直孔氣流場的速度云圖和吸孔入口處壓力云圖如圖4所示。
(a) 沉孔速度云圖 (b) 沉孔壓力云圖
(c) 錐孔速度云圖 (d) 錐孔壓力云圖
(e) 直孔速度云圖 (f) 直孔壓力云圖
(a) 沉孔速度云圖 (b) 沉孔壓力云圖
(c) 錐孔速度云圖 (d) 錐孔壓力云圖
(e) 直孔速度云圖 (f) 直孔壓力云圖
由圖4可知:負壓為-3.0kPa、吸孔直徑為Φ1.0mm時,沉孔型吸孔入口氣流平均速度為24.59m/s,被吸附種子表面平均氣壓為-0.973kPa;錐孔型吸孔入口氣流平均速度為21.59m/s,被吸附種子表面平均氣壓為-2.602kPa;沉孔型吸孔入口氣流平均速度為31.70m/s,被吸附種子表面平均氣壓為-0.393kPa。
由仿真結(jié)果可知:將吸孔入口氣流速度和種子表面負壓統(tǒng)計如表1所示。
表1 吸孔入口氣流速度和種子表面負壓統(tǒng)計
根據(jù)流體力學(xué)理論,氣流對種子的吸附力為種子的繞流阻力[15]。繞流阻力的計算公式為
(1)
式中C—無因次系數(shù),與物體形狀、表面狀態(tài)和雷偌系數(shù)有關(guān);
ρk—空氣密度(g/cm3);
A—種子在垂直于運動方向平面上的投影面積(m2);
vx—吸孔附近氣體流速(m/s)。
由表1可知:在負壓和吸孔直徑相同時,直孔吸孔入口的氣流速度最大,沉孔與錐孔吸孔入口氣流速度相差不大;錐孔吸孔種子表面負壓最大,其次是沉孔,直孔吸孔種子表明負壓最小。因此,由氣體擾流公式可知:直孔吸孔入口氣力產(chǎn)生的擾流阻力大,種子的吸附效果好;而沉孔吸孔入口氣流速度與錐孔吸孔入口氣流速度相差不大,所以沉孔吸孔在吸孔入口處產(chǎn)生的擾流阻力相差不大,沉孔與錐孔對種子的吸附效果相近。
在吸孔附近,種子處于具有一定氣體流速的流場中,且氣體流場的分布呈現(xiàn)出放射狀的圓錐體,設(shè)錐角為γ。假設(shè)種子是球形且滾筒上以吸孔為中心的球面上氣流速度相等,由此可知距吸孔中心x處,種子在垂直于運動方向平面的投影面積A[15-16]為
A=2πx2(1-cosγ)
(2)
假定氣體不可壓縮,根據(jù)質(zhì)量守恒定律可得,吸孔為中心的球面上氣流速度vx為
(3)
式中Q—氣體流量(m3/h)。
由文獻[17]可知,滾筒吸孔處的氣流速度為
(4)
式中R—氣體常量;
Te—氣體絕對溫度(K);
k—氣體比熱容(J/kg·℃);
p—氣室負壓(Pa);
P0—大氣壓(Pa)。
假設(shè)氣體不可壓縮,則直徑為d的吸孔內(nèi)氣體流量Q為
(5)
吸孔附近的種子在氣體繞流作用下被吸起時的平衡方程為
(6)
式中dz—種子直徑(m2);
ρz—種子密度(g/cm3)。
由式(1)~ 式(5)聯(lián)立可得
(7)
由式(7)可知:吸孔對種子的作用范圍與種子的物料特性與吸孔結(jié)構(gòu)和負壓有關(guān)。
由表1可知:沉孔截面突變,有氣體回流現(xiàn)象,產(chǎn)生明顯渦流損耗;而錐孔截面漸變,無回流現(xiàn)象,負壓損耗小,所以錐孔吸孔對種子表面負壓大于沉孔吸孔對種子的負壓。錐孔對種子的吸附范圍大于沉孔吸孔,對種子吸附后的穩(wěn)定性優(yōu)于沉孔吸孔。因此,錐孔吸孔吸種效果優(yōu)于沉孔吸孔。
在負壓和吸孔直徑相同時,直孔吸孔入口處氣流速度大與錐孔吸孔入口氣流速度,但由于直孔吸孔在種子在垂直于運動方向平面上的投影面積較小,所以在吸種階段,直孔吸孔和錐孔吸孔在吸孔入口產(chǎn)生擾流阻力大小有待試驗驗證。但是,錐孔吸孔對種子表面的負壓遠大于直孔吸孔對種子負壓,且錐孔吸孔與種子接觸面較大,所以對于小粒徑種子,錐孔吸孔吸附多粒種子的概率增大。
氣力滾筒排種器試驗平臺如圖5所示,播種對象為甘藍種子,直徑為1.4mm。
1.穴盤 2.種箱 3.清種裝置 4.滾筒 5.輸送裝置 6.機架 7.控制箱
利用氣力滾筒排種器試驗平臺,對影響單粒率的吸孔口結(jié)構(gòu)形式,在滾筒轉(zhuǎn)速和負壓相同的條件下進行試驗。
試驗中,滾筒轉(zhuǎn)速10r/min,負壓為-3.0kPa,吸孔直徑分別為Φ0.6、Φ0.8、Φ1.0mm,以甘藍種子為播種對象,選擇孔徑分別為錐孔和直孔進行試驗,統(tǒng)計吸孔單粒率。
進行吸孔直徑對單粒率影響的試驗,相同的試驗條件下,吸孔直徑分別為0.6、0.8、1.0mm時,錐孔吸孔單粒率分別為88.2%、90.3%、91.2%;直孔吸孔的單粒率分別為91.3%、93.2%和94.2%。
試驗表明:對于小粒徑種子,直孔吸孔吸種效果優(yōu)于錐孔吸孔。錐孔吸孔的對種子的吸附作用范圍大,增大吸孔吸附多粒種子的概率。因此,對于小粒徑種子播種,直孔吸孔對種子的吸種效果較好。
1)在吸孔直徑相同時,錐孔吸孔對種子表面負壓大于沉孔吸孔對種子的負壓, 對種子吸附后的穩(wěn)定性優(yōu)于沉孔吸孔。因此,錐孔吸孔吸種效果優(yōu)于沉孔吸孔。
2)在吸孔直徑相同時,直孔吸孔入口氣流速度大于錐孔吸孔入口氣流速度,但是直孔吸孔對種子的作用范圍小于錐孔吸孔。
3)驗證試驗表明:對于小粒徑種子播種,直孔吸孔對種子的吸種效果優(yōu)于錐孔吸孔。
參考文獻:
[1] 華南農(nóng)業(yè)大學(xué).滾筒氣力式蔬菜播種機:中國,ZL20052 0023056.5[P].2006-07-19.
[2] 牛菊菊.滾筒氣力式蔬菜育苗穴盤排種器的研究[D].廣州:華南農(nóng)業(yè)大學(xué),2004.
[3] 陳佳琦.滾筒氣力式蔬菜播種機的設(shè)計與試驗[D].廣州:華南農(nóng)業(yè)大學(xué),2006.
[4] 林敘彬.氣力滾筒式蔬菜播種機的改進設(shè)計與試驗[D].廣州:華南農(nóng)業(yè)大學(xué),2009.
[5] 王廣源.氣力滾筒式蔬菜播種機的結(jié)構(gòu)設(shè)計與研究[D].廣州:華南農(nóng)業(yè)大學(xué),2012.
[6] 陳進,李耀明,王希強,等.氣吸式排種器吸孔氣流場的有限元分析[J].農(nóng)業(yè)機械學(xué)報,2007,38(9):59-62.
[7] 袁月明,馬旭,金漢學(xué),等.氣吸式水稻芽排種器氣室流場研究[J].農(nóng)業(yè)機械學(xué)報,2005,36(6):42-45.
[8] 廖慶喜,李繼波,覃國良.氣力式油菜精量排種器氣流場仿真分析[J].農(nóng)業(yè)機械學(xué)報,2009,40(7):78-82.
[9] 廖慶喜,楊波,李旭,等.內(nèi)充氣吹式油菜精量排種器氣室流場仿真與試驗[J].農(nóng)業(yè)機械學(xué)報,2012,43(4):51-54.
[10] 龔智強.氣吸振動盤式精密排種裝置理論與試驗研究[D].鎮(zhèn)江:江蘇大學(xué),2013.
[11] 余佳佳.氣力式油菜精量排種器結(jié)構(gòu)解析于排種過程仿真研究[D].武漢:華中農(nóng)業(yè)大學(xué),2013.
[12] 祁兵.中央集排氣送式精量排種器設(shè)計與試驗研究[D].北京:中國農(nóng)業(yè)大學(xué),2014.
[13] 叢錦玲.油菜小麥兼用型氣力式精量排種系統(tǒng)及其機理研究[D].武漢:華中農(nóng)業(yè)大學(xué),2014.
[14] 朱紅鈞,林元華,謝龍漢. FLUENT流體分析工程案例精講[M].北京:電子工業(yè)出版社,2013.
[15] Fallak S S, Sverker P E P.Vacuum nozzle design for seed metering[J].Transactions of the ASAE,1984,27(3):688-696.
[16] Shafi S, Holmes R G. Air-jet seed metering a theoretical and experimental study[J].Transactions of the SAE,1990,33(5):1432-1438.