頓國強 吳星澎 紀欣鑫 紀文義 馬洪巖

(1.哈爾濱劍橋?qū)W院智能農(nóng)機裝備工程實驗室, 哈爾濱 150069; 2.東北林業(yè)大學(xué)機電工程學(xué)院, 哈爾濱 150040;3.東北農(nóng)業(yè)大學(xué)工程學(xué)院, 哈爾濱 150030; 4.黑龍江睿龍創(chuàng)新科技有限責任公司, 哈爾濱 150050)

0 引言

農(nóng)作物種植過程中施肥與追肥是極其重要的一環(huán)[1-4],合理施肥可以提高農(nóng)作物單位面積產(chǎn)量[5]。現(xiàn)有研究表明,不合理的施肥會對環(huán)境造成污染且浪費資源。因此,肥料的精準施用對于肥料的減施增效、降低環(huán)境污染等具有重要意義[6]。

螺旋排肥器具有結(jié)構(gòu)簡單、安裝方便、輸送量可調(diào)、單圈排肥量穩(wěn)定等優(yōu)點[7],在農(nóng)業(yè)生產(chǎn)中應(yīng)用廣泛。工作時可通過螺旋葉片的旋轉(zhuǎn)強制對顆粒狀或粉狀物料進行輸送[8]。現(xiàn)有螺旋排肥器因特殊排肥原理而在排肥螺旋轉(zhuǎn)速較低時,排肥流量波動較大,無法精量排肥[9]。近年來國內(nèi)外相關(guān)學(xué)者為提高肥料的精準施用效果,針對螺旋排肥器進行了大量研究[10]。

肖文立等[11]設(shè)計了一種四頭螺旋雙行排肥器,通過離散元仿真分析得出,排肥螺旋為4頭、螺距為24 mm時排肥性能最佳,并通過田間試驗證明了仿真優(yōu)化的正確性。薛忠等[12]通過仿真的方式分析出影響螺旋排肥器排肥均勻性及排肥量的主次因素依次為排肥螺旋轉(zhuǎn)速、螺距、螺旋葉片直徑。楊文武等[13]針對排肥器排肥流量隨螺旋轉(zhuǎn)動變化而造成排肥均勻性降低的問題,設(shè)計了一種側(cè)面出肥式螺旋排肥器,并通過離散元仿真與臺架試驗相結(jié)合的方法,優(yōu)化得出最佳結(jié)構(gòu)。上述研究證明因螺旋排肥器結(jié)構(gòu)缺陷,其轉(zhuǎn)動時排肥均勻性受到排肥螺旋轉(zhuǎn)速及排肥口結(jié)構(gòu)形狀的影響較大。因此,本研究通過分析瞬時排肥特性,確定螺旋排肥器排肥不均原因,將傳統(tǒng)圓形出肥口改進成斜口出肥口,優(yōu)化排肥器結(jié)構(gòu),以期提高排肥流量均勻性。

本文采用改進螺旋排肥器排肥口的方式提升排肥流量均勻性,以實現(xiàn)精控施肥的目的。并通過基于響應(yīng)曲面法的離散元仿真試驗優(yōu)化斜口排肥口結(jié)構(gòu)參數(shù)。運用3D快速成型技術(shù)制造斜口螺旋排肥器,進行臺架試驗,與仿真結(jié)果及傳統(tǒng)螺旋排肥器的排肥均勻性進行對比分析。并設(shè)計一種通過控制轉(zhuǎn)速來實現(xiàn)精控施肥的排肥控制器,以期為螺旋排肥器的優(yōu)化設(shè)計提供參考。

1 排肥器理論計算與結(jié)構(gòu)優(yōu)化

1.1 螺旋排肥器結(jié)構(gòu)組成

如圖1所示,螺旋排肥器主要由排肥器外殼以及排肥螺旋組成,入肥口及排肥口分別裝配于外殼前后兩側(cè)上下端[14]。依據(jù)市購螺旋排肥器,確定排肥螺旋參數(shù):排肥螺旋外徑D=72 mm,螺旋內(nèi)徑d=30 mm,螺距S=30 mm,厚度b=2 mm,排肥口半徑r1=30 mm。

圖1 螺旋排肥器結(jié)構(gòu)示意圖

1.2 螺旋排肥器關(guān)鍵參數(shù)確定

排肥器作業(yè)時肥箱內(nèi)的肥料通過重力作用落入殼體內(nèi),通過排肥螺旋旋轉(zhuǎn)作業(yè)將入肥口處的肥料輸送至排肥口排出。肥料輸送過程中受到排肥螺旋上螺旋葉片的法向推力和切向摩擦力共同作用產(chǎn)生的軸向和周向移動,顆粒軸向位移過程中伴隨著轉(zhuǎn)動[15]。

1.2.1螺旋排肥器轉(zhuǎn)速

螺旋排肥器排肥關(guān)鍵部件為排肥螺旋,其轉(zhuǎn)速增加時,排肥流量均勻性也隨之提高。但排肥螺旋轉(zhuǎn)速大于極限轉(zhuǎn)速nmax時,排肥螺旋于型腔內(nèi)主要對肥料起攪拌作用,對肥料軸向推進位移較小,無法正常排肥[16]。此時,肥料顆粒所受離心力的最大臨界值與其自身重力之間的關(guān)系為

(1)

(2)

當肥料物理特性不同時,其關(guān)系式為

(3)

(4)

(5)

式中r——排肥螺旋半徑,m

m——肥料顆粒質(zhì)量,kg

g——重力加速度,m/s2

ωmax——排肥螺旋最大角速度,rad/s

K——物料綜合系數(shù)

A——物料綜合特性系數(shù),取50[13]

根據(jù)式(5)計算得出螺旋排肥器極限轉(zhuǎn)速nmax為186.339 r/min。

1.2.2理論排肥量

本文選取江蘇晉煤恒盛化工股份有限公司生產(chǎn)的尿素肥料,抽取其中100粒作為試驗材料,測量三維尺寸L、B、H及密度ρ為0.948 g/cm3。由

(6)

(7)

式中L——顆粒長度,mm

B——顆粒寬度,mm

H——顆粒高度,mm

計算顆粒等效直徑D1為2.05 mm;球形率φ為96.12%。根據(jù)螺旋輸送機構(gòu)工作原理[17-18],有螺旋排肥器螺牙平均長度和單圈排肥量為

(8)

(9)

式中μ——排肥器型腔充填系數(shù),取0.45[19]

S——螺距,mm

計算得出Lp=163 mm,q=40.14 g/r。

根據(jù)農(nóng)藝要求,每秒所需極限排肥量為

(10)

螺旋排肥器極限轉(zhuǎn)速下,每秒排肥量為

(11)

式中Q——每秒所需極限排肥量,g/s

Z——每公頃施肥量,kg/hm2

E——排肥器作業(yè)壟寬,m

vb——螺旋排肥器前進速度,m/s

δ——螺旋排肥器極限排肥量,g/s

查閱相關(guān)文獻,確定排肥器前進速度vb為2.5 m/s,作業(yè)壟寬E為0.6 m,每公頃施肥量Z為750 kg/hm2的條件下[20-21],通過式(10)得出所需極限排肥量Q為112.5 g/s,此時螺旋排肥器所需極限轉(zhuǎn)速為168.16 r/min,小于螺旋排肥器極限轉(zhuǎn)速,并通過式(11)計算得出螺旋排肥器極限排肥量為124.66 g/s。根據(jù)農(nóng)藝要求可知,螺旋排肥器所需排肥量Q應(yīng)小于排肥量,即螺旋排肥器可根據(jù)作物需肥量Q的變化,在極限轉(zhuǎn)速nmax范圍內(nèi)調(diào)節(jié)螺旋排肥器轉(zhuǎn)速n。

結(jié)合式(8)～(10)可知,螺旋排肥器實際作業(yè)時,在螺旋排肥器結(jié)構(gòu)尺寸D、d、b、S,所用肥料顆粒ρ、μ,作業(yè)壟寬E及機具前進速度vb一定的情況下,其排肥量主要與轉(zhuǎn)速n相關(guān),轉(zhuǎn)速變化時其排肥均勻性也隨之變化。因此,對排肥均勻性與轉(zhuǎn)速變化之間的關(guān)系進行分析。

1.3 排肥流量不均勻特性分析

圖2為傳統(tǒng)螺旋排肥器離散元仿真一個轉(zhuǎn)動周期內(nèi),排肥螺旋與排肥口接觸處肥料狀態(tài)時序圖。由圖分析可知,隨著排肥螺旋轉(zhuǎn)動,排肥螺旋與排肥口之間相互位置也隨之變化。3.1 s時,肥料送達排肥口,無肥料排出;3.5～3.8 s時,排肥口處肥料量快速增加,瞬時排出量增大。而排肥螺旋與排肥口之間相互位置周期性變化,且每個循環(huán)只有單個螺距參與排肥,使得排肥口處排肥流量總體呈現(xiàn)出上下波動的狀態(tài)[22-24]。隨著排肥螺旋轉(zhuǎn)速增加,單位時間內(nèi)可同時排出多個螺距內(nèi)肥料,因此排肥流量均勻性受到排肥口結(jié)構(gòu)及排肥螺旋轉(zhuǎn)速的影響較大,均勻性隨著轉(zhuǎn)速增加而變好。

圖2 螺旋排肥器瞬態(tài)排肥效果

相關(guān)研究發(fā)現(xiàn)[13],排肥螺旋轉(zhuǎn)速范圍在30～90 r/min區(qū)間時排肥性能較好。本文采用固定排肥螺旋轉(zhuǎn)速為60 r/min,對排肥口結(jié)構(gòu)進行改進的方式以提高螺旋排肥器排肥均勻性。

1.4 排肥口結(jié)構(gòu)改進

為提高排肥流量均勻性,如圖3所示,將螺旋排肥器排肥口結(jié)構(gòu)改進設(shè)計成斜口式。分析圖3可知,影響排肥均勻性的主要參數(shù)為斜口長度x1、斜口角度x2、開口寬度x3。為使肥料顆粒順利從斜口夾角處落入排肥口,并便于斜口排肥口與下方排肥管相連接,將斜口起始端夾角與終止端出口加工為圓角,起始端及終止端圓角半徑尺寸分別為2 mm及16 mm。且為保證強度,將兩圓角之間過渡區(qū)域加工為變半徑圓角。因排肥口沿螺旋軸軸線方向包覆有多個螺距,肥料隨排肥螺旋轉(zhuǎn)動從排肥口斜口開口緩慢排出,未排出肥料隨排肥螺旋轉(zhuǎn)動繼續(xù)向前移動排肥。使排肥過程中多個螺距同時排肥,提高了排肥流量均勻性。

圖3 斜口螺旋排肥器排肥口改進示意圖

2 斜口螺旋排肥器排肥性能仿真優(yōu)化

2.1 仿真模型建立與參數(shù)設(shè)置

使用EDEM構(gòu)建斜口螺旋排肥器仿真模型[24-25],于肥箱內(nèi)建立肥料顆粒工廠,設(shè)定總顆粒數(shù)為10 000,排肥螺旋轉(zhuǎn)速為60 r/min。肥料質(zhì)量監(jiān)測區(qū)設(shè)置于集肥箱處,用以監(jiān)測排肥器3 s內(nèi)排出肥料顆?？傎|(zhì)量,總仿真過程為18 s,仿真模型如圖4所示。

圖4 仿真模型

尿素顆粒球形率φ為96.12%,半徑標準差為0.199 mm,因此將尿素顆粒看作直徑為2.05 mm的球體。選取EDEM內(nèi)置接觸模型Hertz-Mindlin(no slip)作為尿素顆粒與尿素顆粒、尿素顆粒與排肥裝置之間的接觸模型。參考文獻[10],確定所選用的尿素肥料相關(guān)參數(shù)以及尿素肥料與排肥裝置(PLA塑料)之間的接觸力學(xué)參數(shù)如表1所示。

表1 離散元仿真參數(shù)

2.2 排肥性能評價方法

為分析斜口排肥口結(jié)構(gòu)參數(shù)對排肥流量均勻性的影響,采用網(wǎng)格法對斜口螺旋排肥器的排肥流量變化進行數(shù)據(jù)統(tǒng)計。任意截取排肥器穩(wěn)定排肥后排肥口下方監(jiān)測區(qū)3 s內(nèi)肥料堆積質(zhì)量,以0.5 s時間間隔分成6個時間段,分別統(tǒng)計7個時間點對應(yīng)的監(jiān)測區(qū)內(nèi)肥料堆積質(zhì)量,記為A0、A1、…、A6,計算其0.5 s內(nèi)排肥質(zhì)量為

mi=Ai-Ai-1(i=1,2,…,6)

(12)

式中mi——監(jiān)測區(qū)間第i個時間段內(nèi)的排肥質(zhì)量,g

Ai——第i個時間點監(jiān)測區(qū)內(nèi)肥料堆積質(zhì)量,g

分別計算各時間段內(nèi)肥料平均排肥質(zhì)量、標準差以及排肥流量變異系數(shù)σ。

排肥器的排肥流量均勻性與排肥流量變異系數(shù)σ緊密相關(guān),排肥流量變異系數(shù)σ越低證明排肥流量均勻性越佳。因此將排肥流量變異系數(shù)σ作為試驗評價標準,分析轉(zhuǎn)速固定時不同斜口結(jié)構(gòu)參數(shù)下的排肥流量均勻性變化趨勢。

2.3 試驗因素

采用二次通用旋轉(zhuǎn)組合試驗設(shè)計方法,對斜口長度x1、斜口角度x2、開口寬度x3進行優(yōu)化,其因素水平編碼如表2所示。

表2 因素編碼

肥料顆粒被排肥螺旋推送至斜口排肥口時,一部分肥料顆粒與斜口角度x2相接觸,若斜口開口過小,則此部分肥料不易及時排出而繼續(xù)向前推進。因此,為保證排肥穩(wěn)定性,斜口角度x2及開口寬度x3不宜過小,選擇斜口角度x2為30°～80°,開口寬度x3為30～55 mm。斜口對側(cè)長度L1隨斜口角度x2及斜口長度x1的增加而增加。為確保排肥口出口形狀,當斜口角度x2為30°時,斜口長度x1不應(yīng)小于80 mm,因此,確定斜口長度x1范圍為80～130 mm。

2.4 試驗方案與結(jié)果

表3為斜口螺旋排肥器試驗方案及EDEM仿真試驗結(jié)果。通過Design-Expert 8.0軟件對表3進行回歸分析,分析試驗因素變化時,排肥流量變異系數(shù)的變化規(guī)律,其方差分析結(jié)果見表4。分析得知,排肥流量變異系數(shù)回歸模型為極顯著(P=0.000 2<0.01)?；貧w模型失擬項不顯著(P=0.735 1>0.05),表明其回歸模型擬合較好。對于排肥流量變異系數(shù)σ來說,x1、x2、x3均影響顯著(P<0.05),其中x2、x3影響極顯著(P<0.01),具體影響趨勢為x3>x2>x1,各因素與排肥流量變異系數(shù)σ的回歸方程為

表3 試驗方案與結(jié)果

表4 方差分析

σ=-327.595 8+5.039 9x1+4.075 0x2+0.426 0x3-

0.020 1x1x2+0.017 8x1x3+0.023 3x2x3-

(13)

2.5 響應(yīng)曲面分析

采用響應(yīng)曲面圖的方式,分析試驗因素與試驗指標之間的關(guān)系。

如圖5a所示,隨著斜口角度的增加,排肥流量變異系數(shù)于斜口長度較低處呈增大趨勢,于斜口長度較高處呈先增大后減小趨勢。排肥流量變異系數(shù)隨著斜口長度增加,于斜口角度較低處呈先增大后減小趨勢,于斜口角度較高處呈逐漸下降趨勢。

圖5 試驗因素對排肥流量變異系數(shù)影響的響應(yīng)曲面

分析圖5b可知,隨著斜口角度的增加,排肥流量變異系數(shù)呈增大趨勢,隨著開口寬度的增加,排肥流量變異系數(shù)呈下降趨勢。

2.6 參數(shù)優(yōu)化

為確定最佳斜口排肥口結(jié)構(gòu)參數(shù),對斜口結(jié)構(gòu)參數(shù)進行優(yōu)化,參考文獻[12],將排肥流量變異系數(shù)σ設(shè)定為小于15%。分析表4可知,x1對排肥流量變異系數(shù)σ影響較小,為減小加工成本,確定斜口長度x1為105 mm,對x2、x3進行優(yōu)化,其優(yōu)化結(jié)果為圖6。圖6中黃色陰影區(qū)域為x2及x3最佳組合區(qū)域,分析可知,x2為30°～44°,x3為45.05～55.00 mm時,其排肥流量變異系數(shù)σ小于15%。

圖6 參數(shù)優(yōu)化分析結(jié)果

進行仿真驗證試驗,為便于樣機加工,其斜口排肥口參數(shù)于優(yōu)化區(qū)間內(nèi)選擇x1=105 mm,x2=30°,x3=50 mm。仿真試驗證明,其優(yōu)化后排肥流量變異系數(shù)σ為14.31%,符合理論優(yōu)化最佳值范圍。

3 臺架試驗

如圖7所示,為驗證排肥口參數(shù)優(yōu)化的正確性,通過臺架試驗測量排肥均勻性并進行對比分析,試驗地點為東北林業(yè)大學(xué)智能農(nóng)機裝備實驗室,肥料選取江蘇晉煤恒盛化工股份有限公司生產(chǎn)的尿素肥料。試驗裝置主要包括傳統(tǒng)螺旋排肥器及斜口螺旋排肥器。確定斜口參數(shù)為x1=105 mm,x2=30°,x3=50 mm。試驗轉(zhuǎn)速分別設(shè)置為30、60、90、120、150、180 r/min,傳送帶前進速度為0.5 m/s。以排肥流量變異系數(shù)及排肥流量作為試驗指標,于穩(wěn)定排肥后記錄數(shù)據(jù),進行3次重復(fù)試驗。排肥流量變異系數(shù)對比分析結(jié)果如圖8所示,斜口螺旋排肥器排肥流量變化趨勢如圖9所示。

圖7 臺架驗證試驗

圖8 臺架試驗驗證結(jié)果

圖9 排肥流量與轉(zhuǎn)速變化擬合曲線

對比仿真驗證結(jié)果與圖8臺架試驗結(jié)果可知,斜口螺旋排肥器轉(zhuǎn)速為60 r/min時,臺架試驗結(jié)果為13.59%,與仿真優(yōu)化結(jié)果偏差0.72個百分點,證明了優(yōu)化的正確性。且臺架試驗結(jié)果表明,斜口螺旋排肥器在不同轉(zhuǎn)速情況下,排肥流量變異系數(shù)σ均低于傳統(tǒng)型螺旋排肥器,表明斜口螺旋排肥器排肥效果更佳。

圖9為排肥流量試驗結(jié)果擬合曲線,其函數(shù)方程為y=0.7x+0.47,決定系數(shù)R2=0.99,排肥流量隨著轉(zhuǎn)速的增加而增加,且兩者之間呈線性正相關(guān),因此,可實現(xiàn)通過調(diào)節(jié)排肥器轉(zhuǎn)速來精控排肥量的設(shè)計目的。

4 電控排肥器試驗

如圖10所示,為實現(xiàn)精控配肥,與黑龍江睿龍創(chuàng)新科技有限責任公司合作開發(fā)一種電控排肥控制器。其控制系統(tǒng)基于方程y=0.7x+0.47進行設(shè)計,施肥時,可通過北斗GPS衛(wèi)星導(dǎo)航系統(tǒng)確定施肥機械的作業(yè)速度。為驗證控制器精控排肥性能,于30～180 r/min任選6個轉(zhuǎn)速進行3次重復(fù)驗證試驗,取其均值作為試驗結(jié)果,并與預(yù)設(shè)值進行對比,結(jié)果如圖11所示。分析可知,基于電控系統(tǒng)控制的排肥器排肥流量與預(yù)設(shè)值平均偏差為3.15個百分點,證明了電控系統(tǒng)可通過控制轉(zhuǎn)速來實現(xiàn)精控施肥,滿足設(shè)計要求。

圖10 電控排肥裝置臺架測試試驗

圖11 排肥流量對比曲線

5 結(jié)論

(1)改進傳統(tǒng)螺旋排肥器出肥口結(jié)構(gòu),并通過理論分析其排肥量以及最高轉(zhuǎn)速,確定影響排肥均勻性的主要因素為斜口角度及開口寬度。

(2)通過二次回歸正交組合試驗設(shè)計方法分析影響排肥流量變異系數(shù)的主要因素,確定斜口排肥口最優(yōu)參數(shù)組合為:斜口長度為105 mm、斜口角度范圍為30°～44°、開口寬度范圍為45.05～55.00 mm時排肥流量變異系數(shù)小于15%。

(3)通過臺架試驗分別測試斜口螺旋排肥器與傳統(tǒng)螺旋排肥器排肥均勻性,其結(jié)果表明,轉(zhuǎn)速 60 r/min時,斜口排肥器排肥流量變異系數(shù)為13.59%,與理論優(yōu)化最佳值范圍相吻合,證明了優(yōu)化結(jié)果的準確性,且其排肥均勻性均顯著優(yōu)于傳統(tǒng)排肥器。設(shè)計一種電控排肥系統(tǒng),通過改變轉(zhuǎn)速調(diào)節(jié)斜口排肥流量滿足玉米田間實際需求。