劉德柱, 周勇, 張國忠, 張夢(mèng)月, 柯燴彬, 楊全軍

(華中農(nóng)業(yè)大學(xué)工學(xué)院， 農(nóng)業(yè)農(nóng)村部長江中下游農(nóng)業(yè)裝備重點(diǎn)實(shí)驗(yàn)室， 武漢 430070)

施肥作業(yè)是實(shí)現(xiàn)農(nóng)業(yè)高產(chǎn)、高效與糧食安全的重要保證[1]，但不合理的施肥會(huì)導(dǎo)致作物減產(chǎn)和環(huán)境污染。再生稻可提高稻田復(fù)種指數(shù)，增加農(nóng)民收入，其追肥以人工撒肥和離心圓盤式撒肥機(jī)撒肥為主[2]。人工撒肥效率低、施肥均勻性差；離心圓盤式撒肥機(jī)作業(yè)幅度寬、效率高，但容易受到風(fēng)力等因素的影響，導(dǎo)致作業(yè)范圍內(nèi)的撒肥不均，且存在施肥過程中施肥量調(diào)節(jié)不準(zhǔn)確[3-4]的問題。因此，研究再生稻追肥機(jī)械，提高施肥均勻性和實(shí)現(xiàn)施肥量準(zhǔn)確調(diào)節(jié)，對(duì)于提高肥料的利用率、減少環(huán)境污染具有重要意義[5-6]。

施肥機(jī)械是施肥作業(yè)的三大支柱之一，按照施肥方式可分為撒肥機(jī)、種肥施肥機(jī)與追肥施肥機(jī)三種[7]。近年來，學(xué)者對(duì)插秧施肥機(jī)械和中耕施肥機(jī)械研究較多。王金峰等[8]對(duì)水田側(cè)深施肥裝置關(guān)鍵部件排肥器和氣力輸送系統(tǒng)進(jìn)行設(shè)計(jì)與分析，通過風(fēng)力將肥料輸送在水稻根部側(cè)深位置，實(shí)現(xiàn)了定位施肥。齊興源等[9]設(shè)計(jì)了一種稻田氣力式變量施肥機(jī)，并采用圓錐擋板噴撒器噴撒肥料顆粒，在滿足幅寬要求的基礎(chǔ)上，還能夠在幅寬方向上變量控制施肥。施印炎等[10]設(shè)計(jì)了一種新型軸分段式變量追肥機(jī)執(zhí)行機(jī)構(gòu)，分析并確定了其關(guān)鍵部件結(jié)構(gòu)和參數(shù)，實(shí)現(xiàn)了對(duì)常規(guī)槽輪排肥器的精準(zhǔn)控制。以上研究實(shí)現(xiàn)了施肥機(jī)具定位深施肥和施肥量的準(zhǔn)確調(diào)節(jié)，但針對(duì)地表追肥的定位施肥研究較少。由于再生稻分蘗期后逐漸封行，植株較高[11-12]，使用撒肥機(jī)存在施肥不均和施肥量調(diào)節(jié)不準(zhǔn)確的缺點(diǎn)，而采用定位施肥的方式可提高肥料利用率，降低農(nóng)業(yè)成本，追肥效果較好。

為了解決上述問題，本文提出了定位雙側(cè)施肥的再生稻追肥方案，即將肥料顆粒施撒于靠近水稻根系兩側(cè)的土壤；還設(shè)計(jì)了一種以交錯(cuò)斜槽可調(diào)式外槽輪排肥器為基礎(chǔ)的再生稻氣送式雙側(cè)施肥裝置，并針對(duì)影響施肥裝置的主要因素進(jìn)行了試驗(yàn)，得到了施肥裝置最佳工作參數(shù)，為再生稻追肥機(jī)整機(jī)的設(shè)計(jì)提供了參考。

1 材料與方法

1.1 施肥裝置設(shè)計(jì)

1.1.1再生稻追肥方案 根據(jù)再生稻的生長特點(diǎn)，采用對(duì)行雙側(cè)定位施肥的方式進(jìn)行再生稻地表追肥，施肥位置如圖1所示。將肥料施入秧苗兩側(cè)，肥料顆粒散布在靠近秧苗根系分布的范圍可以實(shí)現(xiàn)肥料的充分利用，從而避免撒肥機(jī)拋撒不均導(dǎo)致的肥料浪費(fèi)。

再生稻對(duì)行雙側(cè)施肥機(jī)主要由插秧機(jī)底盤和氣送式雙側(cè)排肥裝置組成，如圖2所示。氣送式雙側(cè)排肥裝置主要由風(fēng)機(jī)、電機(jī)、肥箱、排肥器、混肥管、分肥器等部分組成。

注：1—施肥裝置；2—雙排肥口；3—植株；4—顆粒肥料；5—土壤。Note: 1—Fertilizer device; 2—Two Fertilizer outlets; 3—Rice plant; 4—Granulated fertilizer; 5— Soil.圖1 施肥位置Fig.1 Location of fertilizer

注：1—風(fēng)機(jī)；2—電機(jī)；3—肥箱；4—排肥器；5—混肥管；6—分肥器。Note: 1—Fan; 2—Motor; 3—Box of fertilizer; 4—Fertilizer apparatus; 5—Mixed fertilizer tube; 6—Fertilizer distributor.圖2 再生稻對(duì)行雙側(cè)施肥機(jī)Fig.2 Structue of fertilizer applicator

1.1.2施肥機(jī)關(guān)鍵部件與工作原理 氣送式雙側(cè)施排肥裝置試驗(yàn)臺(tái)架如圖3所示，主要由電機(jī)、風(fēng)機(jī)、交錯(cuò)斜槽可調(diào)式外槽輪排肥器、肥箱、混合管、分肥器、電機(jī)調(diào)速器等部分組成。工作時(shí)肥箱內(nèi)的顆粒肥料經(jīng)過交錯(cuò)斜槽可調(diào)式外槽輪排肥器定量排出進(jìn)入混肥管，風(fēng)機(jī)將高速氣流吹入混肥管，混肥管利用文丘里原理[14]在混肥管顆粒肥料入口處產(chǎn)生氣壓差，肥料依靠重力和氣壓差進(jìn)入混合管，被高速氣流吹至分肥器，分肥器分成均勻的兩行經(jīng)過排肥口排至再生稻植株兩側(cè)土壤，實(shí)現(xiàn)再生稻植株兩側(cè)定量施肥。

1.1.3交錯(cuò)斜槽可調(diào)式外槽輪排肥器設(shè)計(jì) 傳統(tǒng)的外槽輪式排肥器常用的槽輪有直槽、斜槽、交錯(cuò)槽等形式，其中直槽結(jié)構(gòu)簡單、生產(chǎn)制作方便，應(yīng)用最為廣泛，但其齒和槽的結(jié)構(gòu)導(dǎo)致其存在排肥的瞬時(shí)斷條及周期性脈動(dòng)[13]，本課題組設(shè)計(jì)的交錯(cuò)斜槽可調(diào)式外槽輪排肥器有效克服了直槽輪式外槽輪排肥器的缺點(diǎn)。其由排肥殼、交錯(cuò)斜槽槽輪、阻塞輪、調(diào)量輪、排肥舌、出肥口等部分組成，如圖4所示。工作時(shí)，肥料隨槽輪的轉(zhuǎn)動(dòng)，兩側(cè)的交錯(cuò)斜槽槽輪同時(shí)工作，通過左右兩邊的調(diào)量輪壓縮交錯(cuò)斜槽的工作長度實(shí)現(xiàn)對(duì)排肥量的調(diào)節(jié)，從而實(shí)現(xiàn)排肥量可調(diào)的穩(wěn)定施肥。

注：1—電機(jī)；2—風(fēng)機(jī)；3—排肥器；4—肥箱；5—混肥器；6—分肥器；7—電機(jī)調(diào)速器。Note: 1—Motor;2—Fan;3—Fertilizer apparatus;4—Box of fertilizer； 5—Mixed fertilizer tube; 6—Fertilizer distributor; 7—Motor governor.圖3 施肥裝置試驗(yàn)臺(tái)架Fig.3 Structure of test bench

注：1—排肥殼；2—交錯(cuò)斜槽槽輪；3—阻塞輪；4—調(diào)量輪；5—排肥舌；6—出肥口。Note: 1—Fertilizer shell; 2—Staggered spiral adjustable outer grooved wheel; 3—Shutoff block; 4—Adjusting wheel; 5—Fertilizer tongue; 6—Fertilizer outlet.圖4 排肥器結(jié)構(gòu)Fig.4 Structure of fertilizer apparatus

1.1.4混肥管設(shè)計(jì) 混肥管是使空氣和肥料混合均勻并運(yùn)送肥料的部件，其結(jié)構(gòu)如圖5所示，主要包括氣流入口、收縮管、肥料入口、氣肥混合室、擴(kuò)張管、輸肥管。

根據(jù)文丘里原理[14]，氣流從混肥管氣流入口進(jìn)入，隨著收縮管管道截面漸減小，流速變大，壓強(qiáng)變小，在肥料入口處產(chǎn)生壓差[15]。交錯(cuò)斜槽排肥器定量排出顆粒肥料，通過重力和壓差的作用進(jìn)入氣肥混合室，并與氣流混合均勻，經(jīng)擴(kuò)張管進(jìn)入輸肥管道。其結(jié)構(gòu)尺寸參照經(jīng)典文丘里管設(shè)計(jì)，混肥管氣流入口和風(fēng)機(jī)出風(fēng)口相連，因此D1與風(fēng)機(jī)出風(fēng)口直徑一致為45 mm，α取21°，氣肥混合室直徑取0.7D1，即31.5 mm，其長度與直徑相等為31.5 mm,β取8°[16]，根據(jù)管道內(nèi)顆粒輸送原理[17]，入口風(fēng)速計(jì)算公式如下。

注：1—?dú)饬魅肟冢?—收縮管； 3—肥料入口；4—?dú)夥驶旌鲜遥?—擴(kuò)張管；6—輸肥管。Note: 1—Flow inlet; 2—Constricted tube; 3—Fertilizer inlet; 4—Mixed room; 5—Expansion tube; 6—Fertilizer delivery tube.圖5 混肥管結(jié)構(gòu)Fig.5 Structure of fertilizer mixing tube

(1)

式中，va為入口風(fēng)速，m·s-1；KL為物料粒度系數(shù)，16；Kd為物料特性系數(shù)，(2～5)×10-5；ρf為顆粒肥密度，為1 430 kg·m-3；L為輸肥管長度，m。

由于輸肥管的長度(L)一般不超過2 m，因此,輸肥管長度對(duì)入口風(fēng)速(va)的影響可忽略不計(jì)。為了保證氣流速度能夠滿足肥料輸送要求，應(yīng)保留30%的空余量。計(jì)算后可得入口風(fēng)速(va)為24.87 m·s-1,取整后可得入口風(fēng)速為25 m·s-1。

顆粒肥料在輸送管道中的料氣輸送比計(jì)算如下。

(2)

(3)

式中μ為料氣輸送比；Gf為單位時(shí)間內(nèi)顆粒肥料輸送質(zhì)量，kg·s-1；Ga為單位時(shí)間氣流質(zhì)量，kg·s-1；ρa(bǔ)為空氣密度，標(biāo)準(zhǔn)狀況下取1.293 kg·m-3。

綜合(1)(2)(3)計(jì)算輸肥管直徑D2。

(4)

根據(jù)再生稻追肥要求，選取追肥量最大為600 kg·hm-2。料氣輸送比μ取值可在0～2[18]，本研究取值為1，計(jì)算圓整后得輸肥管直徑D2為38 mm。

1.1.5分肥器的設(shè)計(jì) 分肥器是氣送式雙側(cè)施肥裝置的核心部件之一，各行施肥量一致性很大程度取決于分肥器的結(jié)構(gòu)和形狀。如圖6所示，輸肥管輸送的空氣和肥料混合流通過分肥器實(shí)現(xiàn)兩行均勻分肥，分肥管的角度θ會(huì)影響氣流運(yùn)動(dòng)狀態(tài)[19]，進(jìn)而影響施肥裝置的施肥性能。根據(jù)再生稻生長特性， 60%以上的植株根系分布在離禾蔸中心 10 cm 土層內(nèi)[20]，為了使肥料顆粒施撒在植株根系附近，兩排肥口的中心距應(yīng)小于20 cm。

注：1—輸肥管；2—排肥口1； 3—排肥口2。Note: 1—Fertilizer delivery tube; 2—Fertilizer outlet 1; 3—Fertilizer outlet 2.圖6 分肥器結(jié)構(gòu)Fig.6 Structure of fertilizer distributor

1.2 臺(tái)架試驗(yàn)

1.2.1試驗(yàn)臺(tái)架及材料 本研究選用肥料為長江中下游常用的顆粒復(fù)合肥(總養(yǎng)分大于等于45%，中化化肥有限公司)，千粒重26.7 g，含水率0.45%，球形率90%。試驗(yàn)臺(tái)架自制，在排肥試驗(yàn)前應(yīng)對(duì)臺(tái)架進(jìn)行調(diào)試，以保證各施肥機(jī)構(gòu)各部件運(yùn)行平穩(wěn)可靠。

1.2.2試驗(yàn)設(shè)計(jì) 根據(jù)已有研究和前期的理論分析，確定影響氣送式雙側(cè)施肥機(jī)構(gòu)排肥性能的主要因素及其工作范圍為：排肥軸轉(zhuǎn)速10～50 r·min-1；槽輪工作長度為10～50 mm；分肥器角度90°～180°。結(jié)合 Design-expert 8.0 軟件進(jìn)行三因素二次回歸旋轉(zhuǎn)組合試驗(yàn)。試驗(yàn)因素水平編碼見表1。

1.2.3評(píng)價(jià)指標(biāo) 依據(jù)施肥機(jī)械質(zhì)量評(píng)價(jià)技術(shù)規(guī)范[21]，施肥機(jī)械的性能應(yīng)滿足總排量穩(wěn)定性變異系數(shù)≤7.8%，各行排肥量一致性變異系數(shù)≤13%，所以本文選取它們?yōu)樵u(píng)價(jià)指標(biāo)。

①總排肥量穩(wěn)定性變異系數(shù)。

(5)

(6)

(7)

(8)

②各行排肥量一致性變異系數(shù)。

(9)

(10)

(11)

(12)

1.2.4參數(shù)優(yōu)化 為得到施肥裝置最佳工作參數(shù)組合，建立目標(biāo)函數(shù)與約束條件。

(13)

采用Design-Expert 8.0中的優(yōu)化設(shè)計(jì)功能，對(duì)上述目標(biāo)函數(shù)進(jìn)行求解。

2 結(jié)果與分析

2.1 施肥結(jié)果分析

根據(jù)三因素二次回歸旋轉(zhuǎn)組合試驗(yàn)設(shè)計(jì)了23組試驗(yàn)，包括15組析因試驗(yàn)和8組誤差檢驗(yàn)試驗(yàn)，結(jié)果見表2。

由表2可知，本試驗(yàn)中施肥裝置總排肥量穩(wěn)定性變異系數(shù)均低于7.8%，各行施肥量一致性變異系數(shù)均低于13%，符合施肥裝置設(shè)計(jì)要求。

表2 施肥裝置試驗(yàn)設(shè)計(jì)與結(jié)果Table 2 Design and results of fertilizer device experiment

2.2 回歸模型建立

采用Design-expert 8.0軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行回歸擬合，建立總排肥量穩(wěn)定性變異系數(shù)(αa)、各行排肥量一致性變異系數(shù)(αb)與排肥軸轉(zhuǎn)速、槽輪工作長度、分肥器角度的二次回歸模型，并進(jìn)行方差分析，結(jié)果見表3、表4。進(jìn)而得出顯著項(xiàng)，并根據(jù)因素編碼公式進(jìn)行回代，得出回歸方程如下。

表4 模型αb方差分析及顯著性檢驗(yàn)Table 4 Variance analysis and significant test of model αb

αa=24.54-0.03X1-0.24X2-0.33X3-0.001X1X2-0.002X1X3-0.000 2X2X3+0.006X12+0.006X22+0.001X32

(14)

αb=10.51+0.22X1-0.28X2-0.18X3-0.004X1X2-0.001X1X3+0.000 4X2X3+0.004X12+0.006X22+0.001X32

(15)

2.3 模型方差分析

由上述回歸方程和方差分析(表3、4)可知，總排肥量穩(wěn)定性變異系數(shù)和各行排肥量一致性變異系數(shù)回歸模型檢驗(yàn)均極顯著(P<0.01)，失擬不顯著(P>0.1)，方程與實(shí)際擬合情況較好。

表3 模型αa方差分析及顯著性檢驗(yàn)Table 3 Variance analysis and significant test of model αa

總排肥量穩(wěn)定性變異系數(shù)回歸模型中，一次項(xiàng)排肥軸轉(zhuǎn)速、槽輪工作長度對(duì)總排肥量穩(wěn)定性變異系數(shù)的影響極顯著，分肥器角度對(duì)總排肥量穩(wěn)定性變異系數(shù)的影響顯著；交互項(xiàng)對(duì)總排肥量穩(wěn)定性變異系數(shù)的影響均不顯著；二次項(xiàng)對(duì)總排肥量穩(wěn)定性變異系數(shù)的影響均極顯著。各因素對(duì)總排肥量穩(wěn)定性變異系數(shù)的影響主次順序?yàn)榕欧瘦S轉(zhuǎn)速、槽輪工作長度、分肥器角度。

各行排肥量一致性變異系數(shù)回歸模型中，一次項(xiàng)排肥軸轉(zhuǎn)速和分肥器角度對(duì)各行排肥量一致性變異系數(shù)影響極顯著，槽輪工作長度對(duì)各行排肥量一致性變異系數(shù)的影響顯著；交互項(xiàng)對(duì)各行排肥量一致性變異系數(shù)的影響均不顯著；二次項(xiàng)中X22、X32對(duì)各行排肥量一致性變異系數(shù)的影響均顯著，X12對(duì)各行排肥量一致性變異系數(shù)的影響不顯著。各因素對(duì)各行排肥量一致性變異系數(shù)的影響主次順序?yàn)榕欧瘦S轉(zhuǎn)速、分肥器角度、槽輪工作長度。

2.4 響應(yīng)面分析

2.4.1排肥軸轉(zhuǎn)速和槽輪工作長度的交互作用

圖7為分肥器角度為135°時(shí)，排肥軸轉(zhuǎn)速和槽輪工作長度的對(duì)總排肥量穩(wěn)定性變異系數(shù)和各行排肥量一致性變異系數(shù)的響應(yīng)面。排肥軸轉(zhuǎn)速在10～50 r·min-1、槽輪工作長度在10～50 mm時(shí)，隨著交錯(cuò)斜槽槽輪工作長度的增加，總排肥量穩(wěn)定性變異系數(shù)先減小后增大，各行排肥量一致性變異系數(shù)也先減小后增大；隨著轉(zhuǎn)速的增加，總排肥量穩(wěn)定性變異系數(shù)先減小后增大，各行排肥量一致性變異系數(shù)逐漸增大。

圖7 交互因素的響應(yīng)面Fig.7 Interaction factors response surface

2.4.2排肥軸轉(zhuǎn)速和分肥器角度的交互作用

圖8為排肥器交錯(cuò)斜槽工作長度為30 mm時(shí)，排肥軸轉(zhuǎn)速和分肥器角度對(duì)總排肥量穩(wěn)定性變異系數(shù)和各行排肥量一致性變異系數(shù)的響應(yīng)面。當(dāng)分肥器角度90～180°、排肥軸轉(zhuǎn)速10～50 r·min-1時(shí)，隨著分肥器角度的增加，總排肥量穩(wěn)定性變異系數(shù)先減小后增大，各行排肥量一致性變異系數(shù)也先減小后增大；隨著排肥軸轉(zhuǎn)速的增加，總排肥量穩(wěn)定性變異系數(shù)逐漸增大，各行排肥量一致性變異系數(shù)也逐漸增大。

圖8 交互因素的響應(yīng)面Fig.8 Interaction factors response surface

2.4.3槽輪工作長度和分肥器角度的交互作用

圖9為排肥軸轉(zhuǎn)速為30 r·min-1時(shí)，槽輪工作長度和分肥器角度對(duì)總排肥量穩(wěn)定性變異系數(shù)和各行排肥量一致性變異系數(shù)的響應(yīng)面。當(dāng)分肥器角度90～180°、槽輪工作長度為10～50 mm時(shí)，隨著分肥器角度的增加，總排肥量穩(wěn)定性變異系數(shù)先減小后增大，各行排肥量一致性變異系數(shù)逐漸增大；隨著槽輪工作長度的增加，總排肥量穩(wěn)定性變異系數(shù)逐漸增大，各行排肥量一致性變異系數(shù)也逐漸增大。

圖9 交互因素的響應(yīng)面Fig.9 Interaction factors response surface

2.5 優(yōu)化試驗(yàn)驗(yàn)證分析

采用Design-Expert 8.0中的優(yōu)化設(shè)計(jì)功能對(duì)目標(biāo)函數(shù)進(jìn)行求解，得到再生稻氣送式雙側(cè)施肥機(jī)構(gòu)各因素最佳參數(shù)組合：排肥軸轉(zhuǎn)速16 r·min-1，槽輪工作長度22 mm，分肥器角度114°，此時(shí)總排肥量穩(wěn)定性變異系數(shù)為1%，各行排肥量一致性變異系數(shù)為2.45%。

將得到的最佳工作參數(shù)進(jìn)行臺(tái)架試驗(yàn)驗(yàn)證，最終得到總排肥量穩(wěn)定性變異系數(shù)為1.34%，各行排肥量一致性變異系數(shù)為2.71%，與優(yōu)化結(jié)果基本一致，均滿足再生稻追肥的要求。

3 討論

以往研究中，定位施肥多在插秧施肥等作業(yè)中，對(duì)再生稻地表追肥研究的定位施肥研究相對(duì)較少，同時(shí)為了減少撒肥機(jī)撒肥不均和施肥量調(diào)節(jié)不準(zhǔn)確，本文根據(jù)再生稻生長特點(diǎn)和追肥要求，提出了再生稻的對(duì)行雙側(cè)定位施肥方案，設(shè)計(jì)了一種再生稻氣送式雙側(cè)施肥裝置，并進(jìn)行了三元二次回歸正交旋轉(zhuǎn)組合設(shè)計(jì)臺(tái)架試驗(yàn)。從響應(yīng)面圖結(jié)合臺(tái)架試驗(yàn)的實(shí)際效果分析可知：隨著轉(zhuǎn)速增加交錯(cuò)斜槽排肥輪的交錯(cuò)槽結(jié)構(gòu)發(fā)揮作用，交錯(cuò)槽兩邊同時(shí)進(jìn)行排肥工作，避免了單邊排肥的齒槽相間的脈動(dòng)施肥，施肥性能較好，當(dāng)轉(zhuǎn)速較大時(shí)，槽輪出現(xiàn)充肥不足，施肥性能有所下降；隨著交錯(cuò)斜槽排肥輪工作長度的增加，斜槽的結(jié)構(gòu)發(fā)揮作用，肥料顆粒經(jīng)過斜槽連續(xù)性排出，施肥穩(wěn)定，當(dāng)斜槽工作長度較長時(shí)，斜槽內(nèi)的肥料顆粒排出滯后，排肥器施肥性能略有下降；分肥器角度對(duì)施肥裝置兩行施肥量一致性變異系數(shù)影響極顯著，在分肥器角度較小時(shí)，出現(xiàn)渦流現(xiàn)象，肥料在肥口出現(xiàn)堆積，影響分肥效果，分肥器角度較大時(shí)分肥器的角度導(dǎo)流效果變差，分肥效果也變差。

對(duì)臺(tái)架試驗(yàn)結(jié)果分析并采用Design-Expert 8.0軟件得出影響施肥裝置施肥性能的主要因素最佳工作參數(shù)組合：排肥軸轉(zhuǎn)速16 r·min-1，槽輪工作長度22 mm，分肥器角度114°，此時(shí)，總排肥量穩(wěn)定性變異系數(shù)為1%，各行排肥量一致性變異系數(shù)為2.45%，施肥裝置滿足再生稻追肥要求。下一步可考慮結(jié)合變量施肥技術(shù)對(duì)再生稻雙側(cè)施肥裝置施肥量實(shí)時(shí)控制，使其適用于不同田塊的再生稻施肥。