宋雷潔,李建平,楊 欣,王鵬飛

(河北農(nóng)業(yè)大學(xué) 機電工程學(xué)院,河北 保定 071001)

0 引言

隨著我國果園的規(guī)?；l(fā)展,病蟲害發(fā)生的面積也隨之增大,果樹的病蟲害防治工作是生產(chǎn)管理過程中一項重要環(huán)節(jié),施藥效果的好壞直接影響果實的品質(zhì)和產(chǎn)量[1]?，F(xiàn)階段我國果園中普遍應(yīng)用風(fēng)送式噴霧技術(shù),可以細化霧滴粒徑,帶動葉片間的擾動,從而增加樹膛中部霧滴的穿透性,提高樹膛內(nèi)部藥液沉積量[2-5],有效減少霧滴的飄移,農(nóng)藥利用率高[6-7]。

塔型噴霧機的噴霧出風(fēng)口由傳統(tǒng)的半圓形變?yōu)榇怪敝本€形,霧滴在軸流風(fēng)機的氣流作用下水平噴射至冠層,在一定程度上避免了向上無靶標(biāo)噴霧的藥液浪費[8]。由于塔型噴霧機立板高度較高,普遍存在頂部氣流流量較少的情況,導(dǎo)致噴霧機上方風(fēng)速較低;而氣流的導(dǎo)流過程在很大程度上影響著噴霧機頂部風(fēng)送效果,風(fēng)機提供的氣流通過噴霧機內(nèi)部導(dǎo)流板的導(dǎo)向作用,對噴霧機內(nèi)部氣流進行分配,改變氣流的分布,從而調(diào)整噴霧機頂部的風(fēng)送效果。不同尺寸、不同安裝角度的導(dǎo)流板會引起不同的導(dǎo)流效果,故可改變上方導(dǎo)流板的參數(shù),通過流體仿真軟件對導(dǎo)流板的不同安裝情況進行仿真分析,觀察因?qū)Я靼灏惭b參數(shù)的不同而引起的導(dǎo)流效果的變化,找出導(dǎo)流板的最佳安裝參數(shù)組合。

1 塔型風(fēng)送式果園噴霧機結(jié)構(gòu)與原理

塔型風(fēng)送式噴霧機主要由風(fēng)機、立板、導(dǎo)流裝置、噴頭、藥箱及藥泵等工作部件組成,通過AIP軟件繪制的塔型噴霧機的三維模型如圖1所示。

1.噴頭 2.導(dǎo)流板 3.藥箱 4.藥泵 5.風(fēng)機 6.立板圖1 塔型風(fēng)送式果園噴霧機三維模型圖Fig.1 Three-dimensional model diagram of tower type wind-driven orchard sprayer

噴霧機工作時,風(fēng)機提供強大的氣流,經(jīng)內(nèi)部導(dǎo)流板導(dǎo)向后輸送至出風(fēng)口,藥液在高速氣流的吹送下與空氣撞擊而霧化,被輸送至果樹。由于塔型噴霧機高度較高,若導(dǎo)流板安裝不合理,則會出現(xiàn)風(fēng)機附近氣流較強而頂端氣流少的情況,導(dǎo)致噴霧效果不均勻。因此,需通過改變導(dǎo)流板的安裝參數(shù)、調(diào)整導(dǎo)流板對腔體內(nèi)部氣流的導(dǎo)向作用,盡可能將氣流均勻分布于腔體內(nèi)各處,進而使得各出風(fēng)口所形成的氣流場基本相同,各出風(fēng)口的噴霧效果達到最佳狀態(tài)。

2 仿真模型的建立及仿真分析

2.1 技術(shù)路線

塔型風(fēng)送式果園噴霧機送風(fēng)裝置風(fēng)場流體分析的技術(shù)路線如圖2所示。針對研究內(nèi)容,查閱與之相關(guān)的文獻和資料,利用三維設(shè)計軟件建立仿真分析的流體域模型,通過仿真分析軟件對流體域模型進行分析計算,并通過觀察計算結(jié)果找出最佳組合參數(shù)。

圖2 技術(shù)路線圖Fig.2 Technical roadmap

2.2 仿真模型的建立

噴霧機立板高H=2 000mm、長L=1 000mm、寬W=85mm、風(fēng)機直徑D=810mm。由于導(dǎo)流板1、導(dǎo)流板2、導(dǎo)流板3(導(dǎo)流板自下而上設(shè)為導(dǎo)流板1、2、3、4、5)距離風(fēng)機位置較近,所以導(dǎo)流板1、導(dǎo)流板2、導(dǎo)流板3的長度不能過長,且安裝角度(導(dǎo)流板與水平方向夾角)設(shè)定為接近風(fēng)機的徑向方向效果最好,故設(shè)置β1=-30°、β2=0°、β3=30°,在導(dǎo)流板不會產(chǎn)生干涉的前提下設(shè)定l1=l2=l3=130mm不變。由于導(dǎo)流板4、導(dǎo)流板5距離風(fēng)機位置較遠,所以導(dǎo)流板的導(dǎo)向作用更加重要。若導(dǎo)流板4、導(dǎo)流板5安裝合理,則能增強噴霧機頂部氣流分布,從而提高噴霧機頂部的風(fēng)送效果。

利用正交試驗法對導(dǎo)流板4、導(dǎo)流板5不同安裝情況進行進一步仿真計算分析,通過正交試驗數(shù)據(jù)總結(jié)出導(dǎo)流板4、導(dǎo)流板5的安裝角度及長度對噴霧機頂部氣流流動分布的影響,并通過觀察氣流的分布效果及各監(jiān)測點處的氣流速度等數(shù)據(jù),找出最佳的安裝參數(shù)組合,使得噴霧機頂部的氣流分布達到最佳狀態(tài)。初定β4=β5=30°,并以30°為基準(zhǔn),上下各調(diào)整10°設(shè)為三水平。以模型之間不產(chǎn)生干涉為基準(zhǔn),設(shè)定l4=l5=250mm,并以250mm為基準(zhǔn),上下各調(diào)整50mm設(shè)為三水平。根據(jù)以上參數(shù)所設(shè)定的因素水平表,如表1所示。

表1 因素水平表Table 1 Factor level table

利用AIP軟件按照L9(34)正交試驗表的9種試驗參數(shù)建立對應(yīng)的流體域模型,利用STAR CCM+軟件分別對這9種情況進行內(nèi)部氣流分布的分析計算,通過觀察這9組仿真試驗的計算結(jié)果,對比導(dǎo)流板處氣流流動的效果、監(jiān)測點處氣流速度等參數(shù),總結(jié)出導(dǎo)流板的安裝參數(shù)對氣流流動的影響,從而進一步確定導(dǎo)流板的長度及安裝角度的最優(yōu)參數(shù)組合。

2.3 仿真分析

應(yīng)用STAR CCM+軟件進行仿真計算分析之前,首先需要對三維模型進行簡化,然后抽出仿真計算所需流體域。利用AIP軟件生成的流體域模型如圖3所示。

圖3 流體域模型Fig.3 Fluid domain model

仿真計算步驟如下:運行STAR CCM+軟件,將利用AIP軟件創(chuàng)建的流體域模型導(dǎo)入STAR CCM+軟件中,對其分割邊界并設(shè)定邊界名稱及類型。本次仿真設(shè)置網(wǎng)格大小為3mm,物理模型選擇Gas、Coupled Flow、Ideal Gas、Steady、Turbulent、K-Epsilon Turbulent;設(shè)置入口速度為15m/s,將最大迭代步數(shù)(Maximum Steps)設(shè)置為1 000步,新建Scalar Scene以檢測氣流的分布情況。在導(dǎo)流板所在的出風(fēng)口處,坐標(biāo)為(0.04,1.24,-0.02)、(0.04,1.54,-0.02)兩點分別設(shè)置監(jiān)測點,在Reports目錄中建立監(jiān)測點的速度和壓力reports,分別對每一監(jiān)測點的reports生成monitor及plot,來監(jiān)測出風(fēng)口處氣流速度及壓力;最后,開始進行計算,并觀察計算結(jié)果。

2.4 仿真結(jié)果

依照初定參數(shù)所設(shè)計的模型進行仿真分析,第1組仿真試驗(β4=20°、β5=20°、l4=200mm、l5=200mm)的仿真結(jié)果如圖4所示。

圖4 第1組仿真試驗結(jié)果Fig.4 Group 1 simulation test results

2.4.1 判斷模型

通過觀察云圖可以看出:導(dǎo)流板4、導(dǎo)流板5的周圍存在部分深藍色低速區(qū)域。深藍色區(qū)域的面積越大,則說明低速區(qū)域越大、氣流速度分布越不均勻。通過計算藍色低速區(qū)域的面積、觀察監(jiān)測點處的氣流速度等參數(shù)的變化來對比各項仿真試驗的氣流分布情況,從中找出最優(yōu)的結(jié)果。

將低速區(qū)圈出,繪制成二維圖形,建立直角坐標(biāo)系如圖5、圖6所示。利用定積分公式求出低速區(qū)域面積。由仿真云圖可觀察出左右兩側(cè)效果大致相同,呈近似對稱分布,故只對比單側(cè)氣流分布情況,以左側(cè)藍色低速區(qū)域面積為對象進行對比。

圖5 導(dǎo)流板5低速區(qū)域面積示意圖Fig.5 Schematic diagram of the area of the low velocity region of the deflector 5

圖6 導(dǎo)流板4低速區(qū)域面積示意圖Fig.6 Schematic diagram of the area of the low velocity region of the deflector 4

根據(jù)定積分公式,所求區(qū)域面積計算公式為

針對于第1組仿真試驗導(dǎo)流板5周圍低速區(qū)域面積(A區(qū)域)示意圖列出面積公式為

同理,列出導(dǎo)流板4周圍低速區(qū)域面積(B區(qū)域)公式為

其中,k為相應(yīng)線段的斜率。

分別求出A、B區(qū)域的面積。以下8組仿真試驗方案計算方法均與第1組仿真試驗相同,最終計算結(jié)果匯總?cè)绫?所示。

表2 實驗數(shù)據(jù)匯總表Table 2 Experimental data summary table

根據(jù)戴奮奮[9]研究所提出的末速度原則推導(dǎo)出公式為

其中,V1為風(fēng)機出口風(fēng)速(m/s);H1為風(fēng)機出口垂直高度(m);V2為氣流到達樹體末速度(m/s);K為考慮到風(fēng)量的沿程損失而設(shè)定的系數(shù),通常取K=1.3～1.8。

戴奮奮在蘋果園進行了試驗研究,結(jié)果表明:蘋果樹末速度取值為9～10m/s。故取V2=10m/s,K=1.5代入公式,計算得出V1=15m/s。這說明,當(dāng)出風(fēng)口風(fēng)速達到15m/s時,噴霧效果符合要求。

第2組仿真試驗(β4=20°、β5=30°、l4=250mm、l5=250mm)的仿真結(jié)果如圖7左圖所示。第3組仿真試驗(β4=20°、β5=40°、l4=300mm、l5=300mm)的仿真結(jié)果如圖7右圖所示。

圖7 第2、3組仿真試驗結(jié)果Fig.7 Group 2,3 simulation test results

由計算結(jié)果可以看出:第2組仿真試驗中,v4=7.947m/s,v5=3.665m/s,S4/S=2.109%,S5/S=1.279%;第3組仿真試驗中v4=7.305m/s,v5=5.401m/s,S4/S=2.977%,S5/S=1.107%。兩組仿真試驗中,導(dǎo)流板4、導(dǎo)流板5周圍均存在較大的深藍色低速區(qū)域,且所設(shè)監(jiān)測點處氣流速度較低,說明參數(shù)還需進行調(diào)整。

第4組仿真試驗(β4=30°、β5=20°、l4=250mm、l5=300mm)的仿真結(jié)果如圖8左圖所示,第5組仿真試驗(β4=30°、β5=30°、l4=300mm、l5=200mm)的仿真結(jié)果如圖8右圖所示。

圖8 第4、5組仿真試驗結(jié)果Fig.8 Group 4,5 simulation test results

由計算結(jié)果可以看出:第4組仿真試驗中,導(dǎo)流板5周圍存在較大面積的低速區(qū)域,S4/S=0.754%,S5/S=3.048%,v4=3.062m/s,v5=5.168m/s;第5組仿真試驗中,導(dǎo)流板4、導(dǎo)流板5周圍的低速區(qū)域面積均有所減少,S4/S=1.028%,S5/S=0.676%,v4=8.716m/s,v5=5.267m/s。兩組仿真實驗中,所設(shè)監(jiān)測點處的氣流速度較低,均未達到要求,參數(shù)還需進行調(diào)整。

第6組仿真試驗(β4=30°、β5=40°、l4=200mm、l5=250mm)的仿真結(jié)果如圖9左圖所示,第7組仿真試驗(β4=40°、β5=20°、l4=300mm、l5=250mm)的仿真結(jié)果如圖9右圖所示。

圖9 第6、7組仿真試驗結(jié)果Fig.9 Group 6,7 simulation test results

由計算結(jié)果可以看出:第6組仿真試驗中,導(dǎo)流板5處的低速區(qū)域基本消失,S4/S=0.401%,S5/S=0.238%。所設(shè)監(jiān)測點處的氣流速度也有所提高,v4=16.551m/s、v5=13.043m/s。第7組仿真試驗中導(dǎo)流板5周圍存在較大低速區(qū)域,S4/S=0.068%,S5/S=2.352%,v4=4.409m/s,v5=4.137m/s。兩導(dǎo)流板處出風(fēng)口氣流速度較低,效果不佳。

第8組仿真試驗(β4=40°、β5=30°、l4=200mm、l5=300mm)的仿真結(jié)果如圖10左圖所示,第9組仿真試驗(β4=40°、β5=40°、l4=250mm、l5=200mm)的仿真結(jié)果如圖10右圖所示。

圖10 第8、9組仿真試驗結(jié)果Fig.10 Group 8,9 simulation test results

由計算結(jié)果可以看出:第8組仿真試驗中,導(dǎo)流板5周圍存在部分低速區(qū)域,S4/S=0.046%,S5/S=1.784%,v4=8.442m/s,v5=4.441m/s,氣流速度較低,效果不佳,參數(shù)仍需調(diào)整;第9組仿真試驗中,導(dǎo)流板4、導(dǎo)流板5處的低速區(qū)域均明顯減小,S4/S=0.065%,S5/S=0.265%,v4=15.246m/s,v5=17.719m/s,出風(fēng)口處氣流速度整體提高,均可達到要求,效果最佳。

2.4.2 數(shù)據(jù)統(tǒng)計

將按L9(34)正交試驗表進行的9種仿真試驗計算結(jié)果匯總?cè)绫?所示。

2.4.3 仿真試驗結(jié)果

通過對比STAR CCM+對9種流體域分析的結(jié)果數(shù)據(jù),可以看出導(dǎo)流板的不同安裝參數(shù)影響著噴霧機內(nèi)部流體的分布。由計算結(jié)果可知:第9組仿真試驗的仿真結(jié)果最優(yōu),由云圖可明顯觀察到調(diào)整β4=β5=40°;導(dǎo)流板4、導(dǎo)流板5處的相對低速區(qū)域面積基本消失,各導(dǎo)流板處監(jiān)測點位置的氣流速度整體較高,效果得到明顯改善。

3 結(jié)論

1)利用AIP軟件對塔型噴霧機的流體域進行建模,并利用STAR CCM+軟件對流體域模型進行仿真分析,通過對比計算結(jié)果,總結(jié)出導(dǎo)流板的不同安裝參數(shù)對噴霧機內(nèi)氣流分布的影響,并從9組試驗中找出最合適的安裝參數(shù)。

2)在l1=l2=l3=130mm、β1=-30°、β2=0°、β3=30°條件下,通過對比9組仿真試驗的計算結(jié)果,得出結(jié)論:第9組仿真試驗中,即當(dāng)β4=40°、β5=40°、l4=250mm、l5=200mm時,v4=15.246m/s,S4/S=0.065%,v5=17.719m/s,S5/S=0.265%,導(dǎo)流板周圍氣流速度能達到相對較好的狀態(tài)。

3)通過對比不同安裝參數(shù)下導(dǎo)流板周圍氣流的分布、出風(fēng)口氣流速度等參數(shù),總結(jié)出導(dǎo)流板的不同安裝參數(shù)對塔型噴霧機出風(fēng)口處氣流變化的影響。基于此可進一步優(yōu)化噴霧機內(nèi)導(dǎo)流板的分布,對塔型噴霧機內(nèi)部結(jié)構(gòu)優(yōu)化的研究方向提供建議和參考。