姜 昆,王永龍,樊嘉榮,張夢(mèng)龍

(1.武漢科技大學(xué) 冶金裝備及其控制教育部重點(diǎn)實(shí)驗(yàn)室,武漢 430081;2.河南新農(nóng)人航空科技開發(fā)有限公司,河南 焦作 454950)

0 引言

農(nóng)作物病蟲害危及我國(guó)食品安全并制約農(nóng)產(chǎn)品的全面供給。傳統(tǒng)手動(dòng)、小型植保機(jī)械作業(yè)存在藥物殘留超標(biāo)、低效率及危害操作人員身體健康等負(fù)面問(wèn)題,不能滿足我國(guó)快速發(fā)展的專業(yè)防治需求[1-2]。多旋翼植保無(wú)人機(jī)具有噴灑高效性、立體性強(qiáng)、勞動(dòng)強(qiáng)度低等特點(diǎn),尤其適用于復(fù)雜地形作業(yè),是植保噴施的高效手段[3]。

由于風(fēng)場(chǎng)參與噴施作業(yè),無(wú)人機(jī)噴施霧滴運(yùn)動(dòng)規(guī)律與傳統(tǒng)噴施存在較大差異[4]。國(guó)內(nèi)外學(xué)者[5-7]對(duì)旋翼植保無(wú)人機(jī)的噴施特性及防治效果進(jìn)行了研究及探討。邱白晶等[8]指出:無(wú)人機(jī)作業(yè)速度及高度對(duì)沉積霧滴濃度、均勻度影響極顯著。秦維彩等[9]使用含示蹤劑的水溶液進(jìn)行噴霧試驗(yàn),得出噴施高度、噴幅對(duì)植株冠層霧滴沉積分布存在影響的結(jié)論。陳盛德等[10-11]以HY-B-10L無(wú)人機(jī)為實(shí)驗(yàn)平臺(tái),研究不同飛行參數(shù)對(duì)植株冠層霧滴沉積分布的影響,同時(shí)指出旋翼豎直風(fēng)場(chǎng)對(duì)霧滴沉積和穿透性影響最顯著。楊風(fēng)波等[12]運(yùn)用計(jì)算流體力學(xué)仿真,研究霧滴在風(fēng)場(chǎng)作用下的運(yùn)動(dòng)規(guī)律及漂移情況,并得出該方法能較準(zhǔn)確地模擬實(shí)際噴施作業(yè)的結(jié)論。

以往的研究多集中于飛行、噴施參數(shù)對(duì)沉積結(jié)果的影響,實(shí)際作業(yè)中,植保無(wú)人機(jī)飛行載荷是重要變量且不斷下降,不應(yīng)被忽視。藥液釋放致飛行載荷減小,間接改變旋翼轉(zhuǎn)數(shù)及機(jī)身仰俯角,導(dǎo)致旋翼風(fēng)場(chǎng)變化并影響噴施結(jié)果。為此,以紅旗n-10型植保無(wú)人機(jī)為例,結(jié)合數(shù)值分析與場(chǎng)地試驗(yàn),研究飛行載荷下降對(duì)霧滴沉積規(guī)律的影響。

1 旋翼轉(zhuǎn)速和仰俯角計(jì)算

多旋翼植保無(wú)人機(jī)巡航作業(yè)過(guò)程中,受空氣阻力F、飛行重力G、旋翼總拉力T的作用,如圖1所示。

圖1 多旋翼植保無(wú)人機(jī)平飛受力圖

定義:無(wú)人機(jī)中心點(diǎn)為坐標(biāo)原點(diǎn),行進(jìn)負(fù)方向?yàn)閄正向,機(jī)身左側(cè)至右側(cè)為Y正向,垂直地面向上為Z正向。

由圖1得無(wú)人機(jī)平飛受力平衡方程為

F=G·tanθ

(1)

(2)

其中

(3)

(4)

CD=CD1(1-sin3θ)+CD2(1-cos3θ)

(5)

式中F—無(wú)人機(jī)飛行阻力(N);

G—無(wú)人機(jī)飛行載荷(N);

T—旋翼總拉力(N);

θ—機(jī)身仰俯角(°);

ρ—飛行環(huán)境空氣密度(kg/m3);

V—平飛速度(m/s);

S—無(wú)人機(jī)最大橫截面積(m2);

N—旋翼轉(zhuǎn)數(shù)(r/min);

DP—旋翼直徑(m);

CT—旋翼拉力系數(shù);

CD—無(wú)人機(jī)阻力系數(shù);

CD1—仰俯角為0°時(shí),無(wú)人機(jī)阻力系數(shù);

CD2—仰俯角為90°時(shí),無(wú)人機(jī)阻力系數(shù)。

根據(jù)上式得飛行載荷G和仰俯角θ,飛行載荷G與旋翼轉(zhuǎn)數(shù)N的函數(shù)關(guān)系為

(6)

(7)

前后旋翼轉(zhuǎn)數(shù)差產(chǎn)生仰俯角,因轉(zhuǎn)數(shù)差不大,本文所求旋翼轉(zhuǎn)數(shù)為平均旋翼轉(zhuǎn)數(shù)。參考多旋翼飛行器設(shè)計(jì)與控制[13]提供的實(shí)驗(yàn)數(shù)據(jù),結(jié)合紅旗n-10所采用2880型旋翼具體參數(shù),求解旋翼拉力系數(shù)CT。帶入數(shù)據(jù)得出:n-10無(wú)人機(jī)在作業(yè)速度V為0、1、2、3m/s時(shí),旋翼轉(zhuǎn)數(shù)N與飛行負(fù)載G的函數(shù)關(guān)系如圖2所示;仰俯角θ與飛行負(fù)載G的函數(shù)關(guān)系如圖3所示。其中,空氣密度ρ取1.29kg/m3。

大面積作業(yè)時(shí),植保無(wú)人機(jī)按既定航線采取定高、定速巡航作業(yè)方式覆蓋噴施區(qū)域。伴隨藥液釋放,飛行載荷減小達(dá)12kg。為控制作業(yè)速度且保持力平衡,無(wú)人機(jī)旋翼轉(zhuǎn)數(shù)下降可達(dá)32%;同時(shí),飛行阻力有所增加致仰俯角增大明顯。旋翼流場(chǎng)是植保無(wú)人機(jī)噴施重要影響因素。根據(jù)上述函數(shù)關(guān)系求解旋翼轉(zhuǎn)速、仰俯角數(shù)據(jù)用于數(shù)值模擬參數(shù)設(shè)置。

圖2 飛行載荷G與旋翼轉(zhuǎn)數(shù)N函數(shù)圖

圖3 飛行載荷G與仰俯角θ函數(shù)圖

2 模型建立與設(shè)置

2.1 研究對(duì)象

以紅旗n-10型四旋翼小型植保無(wú)人機(jī)為研究對(duì)象(見圖4),主要作業(yè)參數(shù)如表1所示。其配備4個(gè)霧化噴頭,置于相應(yīng)旋翼電機(jī)下方,與旋翼垂直距離0.2m處。

圖4 n-10型四旋翼小型植保無(wú)人機(jī)

2.2 幾何模型建立

2880型旋翼為多曲面復(fù)雜零件,使用3D掃描儀對(duì)無(wú)人機(jī)正、反旋翼進(jìn)行掃描,Geomagic Design X軟件對(duì)掃描面逆向建模。模型導(dǎo)入前處理軟件Gambit中,結(jié)合n-10植保無(wú)人機(jī)實(shí)際噴施范圍及機(jī)身尺寸,在Gambit中建立,長(zhǎng)方體(6m×6m×3m)模擬計(jì)算區(qū)域,并進(jìn)行網(wǎng)格劃分,如圖5所示。機(jī)身結(jié)構(gòu)造型復(fù)雜且不是本文研究重點(diǎn),為減小計(jì)算成本暫不做考慮。

圖5 模擬計(jì)算區(qū)域

本文模擬計(jì)算域由4個(gè)包含旋翼部分的動(dòng)域和1個(gè)靜域組成,動(dòng)域與靜域交界面設(shè)置Interface連接,對(duì)區(qū)域交界處網(wǎng)格加密處理。經(jīng)網(wǎng)格無(wú)關(guān)性驗(yàn)證,網(wǎng)格劃分完畢后的單個(gè)旋翼網(wǎng)格數(shù)約為22.5萬(wàn),網(wǎng)格總數(shù)約為205.8萬(wàn),其處理結(jié)果如圖6所示。

(a) 正反旋翼網(wǎng)格化示意圖

(b) 計(jì)算域網(wǎng)格化示意圖

2.3 模擬參數(shù)設(shè)置

該無(wú)人機(jī)由4個(gè)無(wú)刷直流電機(jī)提供動(dòng)力,相鄰電機(jī)轉(zhuǎn)向相反以減小不必要的滾轉(zhuǎn)力矩,因此相鄰旋翼為反向旋轉(zhuǎn)設(shè)置。選用DPM(Discrete Phase Model)模型對(duì)離散相進(jìn)行噴施模擬,不考慮霧滴間的互相碰撞及霧滴體積變化對(duì)連續(xù)相的作用。

1)離散相噴施參數(shù):材料水密度997.0kg/m3,釋放霧滴數(shù)為500;霧滴噴射角度垂直四旋翼所在平面向下,噴施半錐角為30°;采用離散隨機(jī)游走模型。

2)邊界條件設(shè)置如圖5所示。計(jì)算區(qū)域右側(cè)面為速度入口邊界條件,速度設(shè)置為0(無(wú)人機(jī)懸停)、1、2、3m/s。有效噴幅區(qū)面域設(shè)置為粒子捕捉,其余進(jìn)出口離散邊界條件設(shè)置為粒子逃逸。

3)離散相采用非穩(wěn)態(tài)模擬,添加重力加速度g=9.81m/s2,模擬時(shí)長(zhǎng)10s,迭代時(shí)間步長(zhǎng)0.05s,時(shí)間步數(shù)為200。

3 旋翼風(fēng)場(chǎng)與有效沉積率的結(jié)果與分析

無(wú)人機(jī)滿載到空載的噴施過(guò)程中,霧滴沉積規(guī)律變化,本質(zhì)由旋翼風(fēng)場(chǎng)變化引起。本文先對(duì)不同飛行載荷下的旋翼風(fēng)場(chǎng)進(jìn)行探討,再模擬計(jì)算飛行載荷為22、20、18、16、14、12、10kg,作業(yè)速度為0(懸停)、1、2、3m/s,統(tǒng)計(jì)作業(yè)高度2m時(shí),無(wú)人機(jī)右后(行進(jìn)方向右后)、右前(行進(jìn)方向右前)噴頭噴施霧滴有效沉積率。

3.1 旋翼風(fēng)場(chǎng)變化

圖7為無(wú)外界環(huán)境風(fēng)場(chǎng),平飛速度為3m/s,高度為2m,飛行載荷為22、18、14、10kg時(shí),在觀測(cè)面y=0.442m(右側(cè)旋翼所在豎直平面),旋翼風(fēng)場(chǎng)速度云圖(顯示范圍為-25～10m/s,負(fù)號(hào)表示Z軸負(fù)方向)。

(a) 22kg

(c) 14kg

(d) 10kg

由圖7可知:旋翼風(fēng)場(chǎng)從前、后旋翼位置產(chǎn)生并向下流動(dòng),在此過(guò)程中逐漸減弱,接近地面時(shí)向周邊擴(kuò)散。由于來(lái)向風(fēng)場(chǎng)參與,前后旋翼風(fēng)場(chǎng)在機(jī)腹下約1m處有匯合現(xiàn)象。無(wú)人機(jī)飛行載荷從滿載22kg降至自重10kg過(guò)程中,旋翼轉(zhuǎn)數(shù)下降33%,在計(jì)算域y=0.442m平面,無(wú)人機(jī)正下方1m處風(fēng)速由圖7(a)約16m/s降至圖7(d)約11m/s,同時(shí)機(jī)身仰俯角增加約6°,使旋翼氣流漂移更明顯。

旋翼風(fēng)場(chǎng)流速隨飛行載荷下降而減弱顯著,穿透力減弱;受對(duì)向來(lái)流影響更小,行徑后方旋翼風(fēng)場(chǎng)對(duì)霧滴作用范圍大于前旋翼風(fēng)場(chǎng);機(jī)身仰俯角的增加,使旋翼氣流向后漂移加重;飛行速度越大,旋翼風(fēng)場(chǎng)向后漂移越嚴(yán)重。n-10型無(wú)人機(jī)的4個(gè)噴頭安裝于各旋翼下方0.2m處,該位置旋翼風(fēng)場(chǎng)穩(wěn)定、擾流小,能增強(qiáng)藥液穿透性,降低漂移。

3.2 霧滴沉積統(tǒng)計(jì)及分析

3.2.1 霧滴沉積率評(píng)價(jià)

n-10型無(wú)人機(jī)結(jié)構(gòu)與功能對(duì)稱設(shè)計(jì),為便于霧滴數(shù)量統(tǒng)計(jì),減少運(yùn)算時(shí)長(zhǎng),采用無(wú)人機(jī)右后、右前噴頭單獨(dú)噴施模擬的方式替代整體噴施,且不考慮霧滴破裂、蒸發(fā)等因素。使用Fluent 15.0軟件,計(jì)算作業(yè)速度為0、1、2、3m/s,飛行高度2m,飛行載荷為22、20、18、16、14、12、10kg時(shí),對(duì)應(yīng)旋翼轉(zhuǎn)數(shù)和仰俯角條件下無(wú)人機(jī)的噴施情況,獲取相應(yīng)工況下的霧滴有效沉積率(噴施霧滴在有效噴幅區(qū)內(nèi)的數(shù)量占噴施霧滴總數(shù)比),以評(píng)判各飛行載荷下的霧滴漂移情況。

霧滴有效沉積率計(jì)算見式(8)～式(10)。有效噴幅區(qū)域S是以噴嘴豎直所在直線為中心,r為半徑的圓形區(qū)域,S平面與噴嘴豎直距離為h。

S=π·r2

(8)

r=h·tanβ

(9)

(10)

式中β—噴頭半錐角(°);

θ—機(jī)身仰俯角(°);

h—噴嘴與地面垂直高度(m);

r—有效噴幅區(qū)半徑(m);

S—有效噴幅區(qū)面積(m2);

NS(n,θ)—旋翼轉(zhuǎn)數(shù)為n、機(jī)身仰俯角為θ時(shí)霧滴在有效噴幅區(qū)S的個(gè)數(shù);

N(n,θ)—旋翼轉(zhuǎn)數(shù)為n、機(jī)身仰俯角為θ時(shí)噴頭噴出霧滴總數(shù);

D(n,θ)—旋翼轉(zhuǎn)數(shù)為n、機(jī)身仰俯角為θ時(shí)霧滴有效沉積率(%)。

3.2.2 霧滴沉積分析

無(wú)人機(jī)右后噴頭霧滴有效沉積率模擬結(jié)果,如表2所示。

由表2可見:作業(yè)速度及飛行載荷不同,霧滴有效沉積率不同。滿載且懸停時(shí),右后噴頭霧滴有效沉積率可達(dá)85.1%;作業(yè)速度3m/s、飛行載荷10kg時(shí),霧滴有效沉積率低至31.2%。無(wú)人機(jī)飛行載荷從22kg降至10kg過(guò)程中,平飛速度為0、1、2、3m/s時(shí),作業(yè)初段與末段相較,右后噴頭噴施霧滴有效沉積率分別相差6%、9.1%、10.7%、11.7%。平飛速度相同,大飛行載荷作業(yè)下的右后噴頭霧滴有效沉積率優(yōu)于小飛行載荷作業(yè)時(shí)。安全作業(yè)速度內(nèi),這種差距隨作業(yè)速度提高更為顯著。

大飛行載荷作業(yè)時(shí),無(wú)人機(jī)強(qiáng)旋翼風(fēng)場(chǎng)使更多噴施霧滴被攜帶并快速向有效噴幅區(qū)運(yùn)動(dòng),一定程度上削弱了正面來(lái)流對(duì)霧滴的影響,霧滴群總體抗漂移性更好。同一平飛速度,無(wú)人機(jī)小飛行載荷作業(yè)擁有更大機(jī)身仰俯角,且速度越高、仰俯角越大。仰俯角的增加使旋翼所在平面及噴嘴相對(duì)無(wú)人機(jī)行徑方向逆向轉(zhuǎn)動(dòng),旋翼風(fēng)場(chǎng)隨之轉(zhuǎn)動(dòng),雖仰俯角增加不大,但也增加了霧滴群被旋翼風(fēng)場(chǎng)所攜帶而逃逸出有效噴幅區(qū)的幾率?？傮w而言,同等作業(yè)速度,右后噴頭在大飛行載荷作業(yè)時(shí)的霧滴有效沉積率優(yōu)于小飛行載荷。

如圖8所示:平飛速度為0、1、2、3m/s,植保無(wú)人機(jī)飛行載荷與右后噴頭噴施霧滴有效沉積率呈線性關(guān)系,決定系數(shù)R2分別為0.810 6、0.939 9、0.944 3、0.915 5,線性方程回歸效果顯著。在以上作業(yè)速度中,植保無(wú)人機(jī)飛行載荷下降對(duì)右后噴頭霧滴有效沉積率影響顯著。

圖8 右后噴頭在各飛行載荷的霧滴有效沉積率

無(wú)人機(jī)右前噴頭霧滴有效沉積率模擬結(jié)果,如表3所示。

由表3可見:滿載且懸停時(shí),右前噴頭霧滴有效沉積率為83.9%;作業(yè)速度為3m/s、飛行載荷10kg時(shí),霧滴有效沉積率低至27.3%。無(wú)人機(jī)飛行載荷從22kg降至10kg的過(guò)程中,平飛速度為0、1、2、3m/s時(shí),作業(yè)初段與末段相較,右前噴頭噴施霧滴有效沉積率分別相差3.2%、6.5%、16.1%、10.5%。平飛速度相同,大飛行載荷作業(yè)下的右前噴頭霧滴有效沉積率優(yōu)于小飛行載荷作業(yè)時(shí),安全作業(yè)速度內(nèi),這種差距隨作業(yè)速度提高更為顯著。

表3 右前噴頭各條件下霧滴有效沉積率Table 3 Effective deposition rate of droplets under the conditions of the right front nozzle %

低速作業(yè)時(shí),前噴頭有效沉積率與后噴頭大致相同。作業(yè)速度超過(guò)2m/s時(shí),受更強(qiáng)對(duì)向來(lái)流作用,前旋翼流場(chǎng)向后移,并與后旋翼流場(chǎng)匯合,后旋翼流場(chǎng)移動(dòng)幅度更小。在流場(chǎng)共同作用下,前噴頭噴施霧滴逃逸出有效噴幅區(qū)的比例相較于后噴頭更大,因此前噴頭的霧滴有效沉積率低于后噴頭。

如圖9所示:平飛速度為0、1、2、3m/s,植保無(wú)人機(jī)飛行載荷與右前噴頭噴施霧滴有效沉積率呈線性關(guān)系,決定系數(shù)R2分別為0.804 0、0.726 7、0.848 5、0.937 6,線性方程回歸效果顯著。在以上作業(yè)速度中,植保無(wú)人機(jī)飛行載荷下降對(duì)右前噴頭霧滴有效沉積率影響顯著。

綜上所述,植保無(wú)人機(jī)從滿載到空載的巡航過(guò)程中,前后噴頭噴施霧滴有效沉積率都隨飛行載荷減小而降低。作業(yè)速度為3m/s時(shí),無(wú)人機(jī)滿載與空載的霧滴有效沉積率相差約11%,差距隨作業(yè)速度提高更為明顯。

4 試驗(yàn)

4.1 無(wú)人機(jī)噴施試驗(yàn)

為驗(yàn)證飛行載荷對(duì)霧滴沉積的影響,對(duì)n-10型無(wú)人機(jī)與取樣水敏紙垂直高度2m,平飛速度1、2、3m/s,飛行載荷22、20、18、16、14、12、10kg的作業(yè)條件進(jìn)行藥液沉積試驗(yàn),地點(diǎn)選擇河南省武陡縣郭提村試驗(yàn)田,如圖10(a)所示。水敏紙(3cm×8cm)平放并固定于支撐桿頂端,每一試驗(yàn)架次掠過(guò)3張水敏紙,前后水敏紙相距5m,如圖10(b)所示。每架次噴施完成后統(tǒng)計(jì)3張水敏紙霧滴個(gè)數(shù)并記錄其平均數(shù)。試驗(yàn)時(shí),環(huán)境溫度為26°C,相對(duì)濕度53%,風(fēng)速1.0～2.0m/s。

采用Imagej圖形處理軟件統(tǒng)計(jì)水敏紙上的霧滴數(shù),處理過(guò)程如圖10(c)、(d)所示。

圖9 右前噴頭在各飛行載荷的霧滴有效沉積率

(a) 噴灑實(shí)驗(yàn)圖

(b) 布置示意圖

(c) Imagej霧滴處理

(d) Imagej霧滴處理結(jié)果

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

霧滴處理結(jié)果如圖11所示。

圖11 霧滴數(shù)量統(tǒng)計(jì)結(jié)果

由圖11可知:無(wú)人機(jī)平飛速度為1、2、3m/s時(shí),相較滿載與空載,取樣水敏紙上的霧滴數(shù)分別相差53、92、102,滿載噴灑時(shí)的水敏紙上的霧滴數(shù)多于空載噴灑時(shí)。試驗(yàn)霧滴沉積變化趨勢(shì)與仿真計(jì)算霧滴有效沉積率大致吻合,說(shuō)明作業(yè)過(guò)程中,飛行載荷對(duì)噴施霧滴有效沉積有直接影響。大飛行載荷時(shí),無(wú)人機(jī)下方水敏紙上的霧滴數(shù)多,小飛行載荷作業(yè)時(shí)噴施霧滴向四周擴(kuò)散較多,導(dǎo)致旋翼下方水敏紙上霧滴數(shù)更少。

5 結(jié)論

1)巡航作業(yè)中,飛行載荷隨藥液釋放而減小,植保無(wú)人機(jī)旋翼轉(zhuǎn)數(shù)下降可達(dá)32%,仰俯角增加可達(dá)10°,旋翼風(fēng)場(chǎng)隨飛行載荷下降而減弱并逐漸向后漂移。

2)豎直旋翼風(fēng)場(chǎng)會(huì)抑制霧滴漂移或擴(kuò)散,仰俯角增大會(huì)加重霧滴漂移。作業(yè)速度0～3m/s間,無(wú)人機(jī)滿載比空載時(shí)的霧滴沉積率高4%～11%左右,模擬與實(shí)測(cè)基本一致。植保無(wú)人機(jī)飛行載荷下降會(huì)降低霧滴有效沉積,漂移增加。