謝海軍，李志偉，張 靜

(華南農(nóng)業(yè)大學(xué) 工程學(xué)院，廣州 510642)

0 引言

小粒種子是指平均直徑小于Φ3mm的種子[1]，主要包括蔬菜、油菜、苜蓿、谷子、芝麻、花卉、林業(yè)和包衣煙草等。目前，針對(duì)小粒徑種子播種主要有育苗移栽和直播[2]，我國超2/3蔬菜栽培采用穴盤育苗移栽，工廠化穴盤育苗是實(shí)現(xiàn)種苗商品化供應(yīng)主要方式。精量播種機(jī)是工廠化穴盤育苗的核心設(shè)備，前期研究表明，排種器轉(zhuǎn)速、吸孔結(jié)構(gòu)及正負(fù)氣壓對(duì)播種機(jī)的空穴率、單粒率、多粒率3項(xiàng)指標(biāo)有顯著影響[3-5]。近年來，國內(nèi)外學(xué)者對(duì)氣力排種器的研究，有通過試驗(yàn)的手段尋求各因素水平間的最佳組合，也有利用CDF軟件對(duì)排種器氣流場(chǎng)進(jìn)行仿真[6-13]。本文運(yùn)用ANSYS/FLUNET軟件氣力滾筒式排種器氣流場(chǎng)動(dòng)力學(xué)模型，在吸附種子的條件下，仿真分析了吸孔結(jié)構(gòu)分別為沉孔、錐孔和直孔對(duì)排種器吸種效果的影響。

1 吸孔結(jié)構(gòu)模型

在薄板厚度相同的條件下，設(shè)置有沉孔吸孔、錐孔吸孔和直孔吸孔。利用ANSYS軟件FLUNET模塊建立吸孔氣流場(chǎng)區(qū)域的CFD結(jié)構(gòu)模型如圖1所示。圖1中，1為氣流入口邊界，2為氣流出口邊界，3為吸孔，4為種子，其余邊界為固定壁面邊界。采用自由網(wǎng)格劃分方式，單元截面形狀為三角形，單位長度間隔個(gè)數(shù)為100，在吸孔附近對(duì)網(wǎng)格進(jìn)行進(jìn)一步細(xì)化，更有利于計(jì)算結(jié)果。

圖1 吸種孔結(jié)構(gòu)問題區(qū)域

2 吸孔氣流場(chǎng)仿真

2.1 模型假設(shè)

ANSYS/FLUENT是一個(gè)用于分析二維和三維流體流動(dòng)場(chǎng)常用工具[14]，由于在排種器上的吸孔中的氣流場(chǎng)是軸對(duì)稱流動(dòng)，所以氣流場(chǎng)分析可以簡化為二維軸對(duì)稱流體流動(dòng)問題進(jìn)行求解。仿真分析時(shí)對(duì)模型進(jìn)行如下假設(shè)：①流場(chǎng)氣體為不可壓縮氣體；②室溫恒為25℃;③吸孔入口壓力恒為1.01×105Pa，出口壓力為0Pa；④排種器吸孔進(jìn)氣口氣流速度分布是均勻；⑤在所有壁面上施加無滑移邊界條件。

2.2 仿真方法

仿真試驗(yàn)方法為在吸孔直徑分別為Φ0.6mm、Φ0.8mm和Φ1.0mm分別以負(fù)壓為-3.0kPa對(duì)孔口結(jié)構(gòu)為沉孔、錐孔、直孔的吸孔進(jìn)行氣流場(chǎng)仿真，得到3種孔口結(jié)構(gòu)下的吸孔入口中心氣流速度及種子表面平均氣壓，探討沉孔、錐孔、直孔的吸種和攜種效果。

2.3 仿真結(jié)果

負(fù)壓為-3.0kPa時(shí)，孔徑為Φ0.6mm的沉孔、錐孔和直孔氣流場(chǎng)的速度云圖和吸孔入口處壓力云圖如圖2所示。

(a) 沉孔速度云圖 (b) 沉孔壓力云圖

(c) 錐孔速度云圖 (d) 錐孔壓力云圖

(e) 直孔速度云圖 (f) 直孔壓力云圖

由圖2可知：負(fù)壓為-3.0kPa、吸孔直徑為Φ0.6mm時(shí)，沉孔型吸孔入口氣流平均速度為22.79m/s，被吸附種子表面平均氣壓為-0.844kPa；錐孔型吸孔入口氣流平均速度為18.83m/s，被吸附種子表面平均氣壓為-1.976kPa；沉孔型吸孔入口氣流平均速度為36.40m/s，被吸附種子表面平均氣壓為-0.406kPa。

負(fù)壓為-3.0kPa時(shí)，孔徑為Φ0.8mm的沉孔、錐孔和直孔氣流場(chǎng)的速度云圖和吸孔入口處壓力云圖如圖3所示。

由圖3可知：負(fù)壓為-3.0kPa、吸孔直徑為Φ0.8mm時(shí)，沉孔型吸孔入口氣流平均速度為24.84m/s，被吸附種子表面平均氣壓為-0.819kPa；錐孔型吸孔入口氣流平均速度為22.65m/s，被吸附種子表面平均氣壓為-2.651kPa；沉孔型吸孔入口氣流平均速度為33.09m/s，被吸附種子表面平均氣壓為-0.469kPa；

負(fù)壓為-3.0kPa時(shí)，孔徑為Φ1.0mm的沉孔、錐孔和直孔氣流場(chǎng)的速度云圖和吸孔入口處壓力云圖如圖4所示。

(a) 沉孔速度云圖 (b) 沉孔壓力云圖

(c) 錐孔速度云圖 (d) 錐孔壓力云圖

(e) 直孔速度云圖 (f) 直孔壓力云圖

(a) 沉孔速度云圖 (b) 沉孔壓力云圖

(c) 錐孔速度云圖 (d) 錐孔壓力云圖

(e) 直孔速度云圖 (f) 直孔壓力云圖

由圖4可知：負(fù)壓為-3.0kPa、吸孔直徑為Φ1.0mm時(shí)，沉孔型吸孔入口氣流平均速度為24.59m/s，被吸附種子表面平均氣壓為-0.973kPa；錐孔型吸孔入口氣流平均速度為21.59m/s，被吸附種子表面平均氣壓為-2.602kPa；沉孔型吸孔入口氣流平均速度為31.70m/s，被吸附種子表面平均氣壓為-0.393kPa。

由仿真結(jié)果可知：將吸孔入口氣流速度和種子表面負(fù)壓統(tǒng)計(jì)如表1所示。

表1 吸孔入口氣流速度和種子表面負(fù)壓統(tǒng)計(jì)

3 仿真結(jié)果分析

根據(jù)流體力學(xué)理論，氣流對(duì)種子的吸附力為種子的繞流阻力[15]。繞流阻力的計(jì)算公式為

(1)

式中C—無因次系數(shù)，與物體形狀、表面狀態(tài)和雷偌系數(shù)有關(guān)；

ρk—空氣密度(g/cm3)；

A—種子在垂直于運(yùn)動(dòng)方向平面上的投影面積(m2)；

vx—吸孔附近氣體流速(m/s)。

由表1可知：在負(fù)壓和吸孔直徑相同時(shí)，直孔吸孔入口的氣流速度最大，沉孔與錐孔吸孔入口氣流速度相差不大；錐孔吸孔種子表面負(fù)壓最大，其次是沉孔，直孔吸孔種子表明負(fù)壓最小。因此，由氣體擾流公式可知：直孔吸孔入口氣力產(chǎn)生的擾流阻力大，種子的吸附效果好；而沉孔吸孔入口氣流速度與錐孔吸孔入口氣流速度相差不大，所以沉孔吸孔在吸孔入口處產(chǎn)生的擾流阻力相差不大，沉孔與錐孔對(duì)種子的吸附效果相近。

在吸孔附近，種子處于具有一定氣體流速的流場(chǎng)中，且氣體流場(chǎng)的分布呈現(xiàn)出放射狀的圓錐體，設(shè)錐角為γ。假設(shè)種子是球形且滾筒上以吸孔為中心的球面上氣流速度相等，由此可知距吸孔中心x處，種子在垂直于運(yùn)動(dòng)方向平面的投影面積A[15-16]為

A=2πx2(1-cosγ)

(2)

假定氣體不可壓縮，根據(jù)質(zhì)量守恒定律可得，吸孔為中心的球面上氣流速度vx為

(3)

式中Q—?dú)怏w流量(m3/h)。

由文獻(xiàn)[17]可知，滾筒吸孔處的氣流速度為

(4)

式中R—?dú)怏w常量；

Te—?dú)怏w絕對(duì)溫度(K)；

k—?dú)怏w比熱容(J/kg·℃)；

p—?dú)馐邑?fù)壓(Pa)；

P0—大氣壓(Pa)。

假設(shè)氣體不可壓縮，則直徑為d的吸孔內(nèi)氣體流量Q為

(5)

吸孔附近的種子在氣體繞流作用下被吸起時(shí)的平衡方程為

(6)

式中dz—種子直徑(m2)；

ρz—種子密度(g/cm3)。

由式(1)～ 式(5)聯(lián)立可得

(7)

由式(7)可知：吸孔對(duì)種子的作用范圍與種子的物料特性與吸孔結(jié)構(gòu)和負(fù)壓有關(guān)。

由表1可知：沉孔截面突變，有氣體回流現(xiàn)象，產(chǎn)生明顯渦流損耗；而錐孔截面漸變，無回流現(xiàn)象，負(fù)壓損耗小，所以錐孔吸孔對(duì)種子表面負(fù)壓大于沉孔吸孔對(duì)種子的負(fù)壓。錐孔對(duì)種子的吸附范圍大于沉孔吸孔，對(duì)種子吸附后的穩(wěn)定性優(yōu)于沉孔吸孔。因此，錐孔吸孔吸種效果優(yōu)于沉孔吸孔。

在負(fù)壓和吸孔直徑相同時(shí)，直孔吸孔入口處氣流速度大與錐孔吸孔入口氣流速度，但由于直孔吸孔在種子在垂直于運(yùn)動(dòng)方向平面上的投影面積較小，所以在吸種階段，直孔吸孔和錐孔吸孔在吸孔入口產(chǎn)生擾流阻力大小有待試驗(yàn)驗(yàn)證。但是，錐孔吸孔對(duì)種子表面的負(fù)壓遠(yuǎn)大于直孔吸孔對(duì)種子負(fù)壓，且錐孔吸孔與種子接觸面較大，所以對(duì)于小粒徑種子，錐孔吸孔吸附多粒種子的概率增大。

4 試驗(yàn)

4.1 試驗(yàn)材料

氣力滾筒排種器試驗(yàn)平臺(tái)如圖5所示，播種對(duì)象為甘藍(lán)種子，直徑為1.4mm。

1.穴盤 2.種箱 3.清種裝置 4.滾筒 5.輸送裝置 6.機(jī)架 7.控制箱

4.2 試驗(yàn)方法

利用氣力滾筒排種器試驗(yàn)平臺(tái)，對(duì)影響單粒率的吸孔口結(jié)構(gòu)形式，在滾筒轉(zhuǎn)速和負(fù)壓相同的條件下進(jìn)行試驗(yàn)。

試驗(yàn)中，滾筒轉(zhuǎn)速10r/min，負(fù)壓為-3.0kPa，吸孔直徑分別為Φ0.6、Φ0.8、Φ1.0mm，以甘藍(lán)種子為播種對(duì)象，選擇孔徑分別為錐孔和直孔進(jìn)行試驗(yàn)，統(tǒng)計(jì)吸孔單粒率。

4.3 試驗(yàn)結(jié)果

進(jìn)行吸孔直徑對(duì)單粒率影響的試驗(yàn)，相同的試驗(yàn)條件下，吸孔直徑分別為0.6、0.8、1.0mm時(shí)，錐孔吸孔單粒率分別為88.2%、90.3%、91.2%；直孔吸孔的單粒率分別為91.3%、93.2%和94.2%。

試驗(yàn)表明：對(duì)于小粒徑種子，直孔吸孔吸種效果優(yōu)于錐孔吸孔。錐孔吸孔的對(duì)種子的吸附作用范圍大，增大吸孔吸附多粒種子的概率。因此，對(duì)于小粒徑種子播種，直孔吸孔對(duì)種子的吸種效果較好。

5 結(jié)論與討論

1)在吸孔直徑相同時(shí)，錐孔吸孔對(duì)種子表面負(fù)壓大于沉孔吸孔對(duì)種子的負(fù)壓， 對(duì)種子吸附后的穩(wěn)定性優(yōu)于沉孔吸孔。因此，錐孔吸孔吸種效果優(yōu)于沉孔吸孔。

2)在吸孔直徑相同時(shí)，直孔吸孔入口氣流速度大于錐孔吸孔入口氣流速度，但是直孔吸孔對(duì)種子的作用范圍小于錐孔吸孔。

3)驗(yàn)證試驗(yàn)表明：對(duì)于小粒徑種子播種，直孔吸孔對(duì)種子的吸種效果優(yōu)于錐孔吸孔。

參考文獻(xiàn)：

[1] 華南農(nóng)業(yè)大學(xué).滾筒氣力式蔬菜播種機(jī)：中國，ZL20052 0023056.5[P].2006-07-19.

[2] 牛菊菊.滾筒氣力式蔬菜育苗穴盤排種器的研究[D].廣州:華南農(nóng)業(yè)大學(xué),2004.

[3] 陳佳琦.滾筒氣力式蔬菜播種機(jī)的設(shè)計(jì)與試驗(yàn)[D].廣州:華南農(nóng)業(yè)大學(xué),2006.

[4] 林?jǐn)⒈?氣力滾筒式蔬菜播種機(jī)的改進(jìn)設(shè)計(jì)與試驗(yàn)[D].廣州:華南農(nóng)業(yè)大學(xué),2009.

[5] 王廣源.氣力滾筒式蔬菜播種機(jī)的結(jié)構(gòu)設(shè)計(jì)與研究[D].廣州:華南農(nóng)業(yè)大學(xué),2012.

[6] 陳進(jìn)，李耀明，王希強(qiáng)，等.氣吸式排種器吸孔氣流場(chǎng)的有限元分析[J]．農(nóng)業(yè)機(jī)械學(xué)報(bào)，2007，38(9)：59-62．

[7] 袁月明，馬旭，金漢學(xué)，等.氣吸式水稻芽排種器氣室流場(chǎng)研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2005，36(6)：42-45.

[8] 廖慶喜，李繼波，覃國良.氣力式油菜精量排種器氣流場(chǎng)仿真分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2009，40(7)：78-82.

[9] 廖慶喜，楊波，李旭，等.內(nèi)充氣吹式油菜精量排種器氣室流場(chǎng)仿真與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2012，43(4)：51-54.

[10] 龔智強(qiáng).氣吸振動(dòng)盤式精密排種裝置理論與試驗(yàn)研究[D].鎮(zhèn)江：江蘇大學(xué)，2013.

[11] 余佳佳.氣力式油菜精量排種器結(jié)構(gòu)解析于排種過程仿真研究[D].武漢：華中農(nóng)業(yè)大學(xué)，2013.

[12] 祁兵.中央集排氣送式精量排種器設(shè)計(jì)與試驗(yàn)研究[D].北京：中國農(nóng)業(yè)大學(xué)，2014.

[13] 叢錦玲.油菜小麥兼用型氣力式精量排種系統(tǒng)及其機(jī)理研究[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.