基于離散元法的分層施肥靴參數(shù)優(yōu)化與試驗(yàn)

宋少龍 張東超 湯智輝 鄭 炫 陳學(xué)庚,3 楊懷君 孟祥金

(1.石河子大學(xué) 機(jī)械電氣工程學(xué)院，新疆 石河子 832000； 2.新疆農(nóng)墾科學(xué)院 機(jī)械裝備研究所，新疆 石河子 832000； 3.農(nóng)業(yè)農(nóng)村部西北農(nóng)業(yè)裝備重點(diǎn)實(shí)驗(yàn)室，新疆 石河子 832000)

施肥技術(shù)對(duì)提高作物肥料利用率具有重要意義[1-3]。目前我國農(nóng)業(yè)生產(chǎn)主要以作物產(chǎn)量為目標(biāo)，大量施用化肥是實(shí)現(xiàn)高產(chǎn)的主要途徑，但由于施肥技術(shù)仍停留在地表分期撒施和全層施肥階段，從而造成了肥料浪費(fèi)及土壤板結(jié)[4-5]。張余莽等[6]指出我國目前化肥使用量是發(fā)達(dá)國家安全上限的1.93倍，但氮肥利用率比發(fā)達(dá)國家低20%～25%，磷肥鉀肥利用率更低。

分層施肥技術(shù)是指將肥料按一定比例分層施于土壤中供作物在不同生長時(shí)期吸收，提高化肥的利用率[7-9]，分層施肥技術(shù)依靠分層施肥裝置實(shí)現(xiàn)。王云霞等[10]調(diào)整分層施肥器的安裝角和施肥片調(diào)節(jié)裝置的工作長度控制肥料配比和施肥量，通過臺(tái)架試驗(yàn)，確定最優(yōu)參數(shù)安裝角為36°，工作長度13～16 mm；其田間試驗(yàn)表明，各層施肥量比例穩(wěn)定，分層效果明顯。姚萬生等[11]設(shè)計(jì)一種可以將種子、種肥和基肥分別置于距地表5.3，11.5和17.6 cm，并且種肥和基肥的比例為3∶7的組合式下位分層施肥裝置，保證不傷害種子、幼苗的同時(shí)也可以滿足中、后期營養(yǎng)的供給。張小麗等[12]設(shè)計(jì)的雙層施肥旋耕播種機(jī)，采用前后布置兩套開溝器完成分層施肥作業(yè)。李輝等[13]設(shè)計(jì)一種將種子和化肥分別施于距地表5.4和11.0 cm處，種肥間層為5.6 cm的垂直分層種施開溝器。上述研究主要通過試驗(yàn)直接驗(yàn)證分層施肥裝置的分肥效果，缺少對(duì)分層施肥結(jié)構(gòu)的理論設(shè)計(jì)過程?；首鳛橐环N顆粒物料，與分層結(jié)構(gòu)的碰撞接觸會(huì)直接影響分層施肥裝置的分層效果[14]。本研究基于化肥顆粒與分層施肥部件的相互作用，擬采用理論計(jì)算的方法設(shè)計(jì)一種新型分層施肥裝置，用響應(yīng)面分析法獲得最有組合參數(shù)，并進(jìn)行臺(tái)架試驗(yàn)和田間試驗(yàn)驗(yàn)證，以期為分層施肥裝置的設(shè)計(jì)研究提供參考。

1 分層施肥靴結(jié)構(gòu)及原理

分層施肥靴結(jié)構(gòu)主要由開溝器、分肥器和均肥器3部分組成(圖1)。開溝器選用對(duì)土壤擾動(dòng)較小的深松鏟結(jié)構(gòu)，由鏟刃和鏟柄組成。均肥器由圓管和周向均布、軸向等距的圓桿組成[15]。分肥器主要由分肥板、排肥板、擋泥板、排肥管和周圍擋板5部分組成，分肥板頂部與上擋板相連，通過調(diào)整分肥板與均肥管軸線的徑向距離，保證上下兩層肥量比為4∶6，底部與排肥管進(jìn)肥口端面相切，確保上下兩層肥料不干涉。

1.開溝器；2.均肥器；3.上擋板；4.排肥板；5.分肥板；6.排肥管前擋板；8.側(cè)擋板；9.擋泥板1.Furrower；2.Fertilizer uniform distribution device；3.Overhead baffle；4.Fertilizer discharging plate；5.Fertilizer separating plate；6.Fertilizer discharging tube；7.Front baffle；8.Side baffle；9.Fender圖1 分層施肥靴總體結(jié)構(gòu)Fig.1 General structure of layered fertilization boot

機(jī)具行進(jìn)時(shí)，開溝器在土壤中開出30 cm的深溝，同時(shí)肥料沿導(dǎo)肥管以一定的初速度進(jìn)入均肥器，肥料顆粒與圓桿發(fā)生碰撞和彈射后，下落形狀由團(tuán)聚狀變?yōu)榫紶?，然后由分肥板按設(shè)計(jì)比例分成兩份，下層肥料沿排肥板流出落在開溝器形成的溝槽，上層肥料沿排肥管流出落在回流土壤的表面，上下兩層肥料實(shí)現(xiàn)分層。

2 分層施肥靴關(guān)鍵部件的理論設(shè)計(jì)

2.1 均肥器

肥料顆粒在均肥器中的運(yùn)動(dòng)過程見圖2，肥料顆粒由導(dǎo)肥管進(jìn)入均肥器與周向均布、軸向等距的圓桿發(fā)生碰撞、彈射，顆粒形態(tài)由團(tuán)聚狀變?yōu)榫紶睢?/p>

圓桿數(shù)量N和圓桿直徑d分別決定了分肥份數(shù)2N和單層通過面積S[16]。肥料顆粒在均肥器的單層通過面積S為：

(1)

式中：D為均肥器內(nèi)徑，mm；d為圓桿直徑，mm。均肥器中心交疊圓桿在圓管內(nèi)形成的通過面積(即均肥器實(shí)際通過面積)Ssum為：

(2)

由式(2)可知，當(dāng)圓管內(nèi)徑不變時(shí)，圓桿直徑越大，數(shù)量越多，則Ssum越小，肥料顆粒在均肥器中的流通效率越低。

1.肥料顆粒；2.圓桿；3.圓管1.Fertilizer particle；2.Round bar；3.Tube圖2 肥料顆粒運(yùn)動(dòng)軌跡與均肥器內(nèi)部結(jié)構(gòu)示意圖Fig.2 Motion trajectory of fertilizer particles and schematic diagram of internal structure of fertilizer equalizer

根據(jù)我國新疆砂礫土質(zhì)肥力情況，棉花和玉米每hm2最大施肥量為450 kg，作業(yè)速度最大值為7 km/h，排肥器單位時(shí)間排肥量151 g/s。外槽輪式排肥器的排肥舌距均肥器入口的距離為30 cm，設(shè)定肥料顆粒離開排肥舌時(shí)的初速度為0 cm/s，則肥料顆粒進(jìn)入均肥器入口的速度為24.5 cm/s，肥料顆粒的堆積密度為0.877 g/cm3，均肥器不發(fā)生堵塞的最大流通面積為：

(3)

式中：Sthrough為最大流通面積，cm2；mt為單位時(shí)間施肥量，g/s；ρ為肥料顆粒的堆積密度，g/cm3；v為肥料流速，cm/s。

肥料顆粒由于重力作用在均肥器內(nèi)豎直方向呈加速運(yùn)動(dòng)狀態(tài)，當(dāng)下肥量一定時(shí)，流速越大，通過面積越小。因此最大通過面積應(yīng)采用肥料顆粒進(jìn)入均肥器的初速度。根據(jù)肥箱的導(dǎo)肥管外徑，設(shè)計(jì)圓管內(nèi)徑為38 mm，考慮到肥料的流通性和加工條件，選取圓桿直徑最大為5 mm，圓桿數(shù)量最多為8根。由此計(jì)算得均肥器最小通過面積Ssum=20.99 cm2，單層最小通過面積S=9.44 cm2，均大于最大通過面積Sthrough=7.03 cm2，符合設(shè)計(jì)要求。

2.2 分肥器

分肥器肥料運(yùn)動(dòng)軌跡及分肥器結(jié)構(gòu)見圖3，肥料顆粒呈均布狀由均肥器落入分肥器。肥料顆粒被均肥器下方的分肥板分成不同比例的上下兩層，其上層落在排肥板形成加速下滑的顆粒流，與分肥板下端碰撞后，沿排肥管流出。下層肥料沿分肥板滑落到排肥板，在慣性的作用下與擋泥板發(fā)生碰撞，經(jīng)下排肥口流出。最后在土壤回流作用下，肥料顆粒在土壤中形成兩層。

1.分肥器；2.肥料顆粒1.Fertilizer separating device；2.Fertilizer particleL1和L2分別為上下兩層肥料出口的縱向和橫向距離，mm；L3和L4為排肥板底端到下層肥料出口的縱向和橫向距離，mm；α為分肥器傾角，(°)；Ltube為排肥管長度，mm；Lplate為排肥板長度，mm；Lsum為排肥板頂端到擋泥板的對(duì)齊距離，mm；δ為排肥板厚度，mm；A為分肥板中心線與圓管中心線的距離，mm；D為圓管內(nèi)徑，mm。L1 and L2 are the vertical and horizontal distances of the fertilizer outlets in the upper and lower layers,mm;L3 and L4 are the vertical and horizontal distances of the bottom of the fertilizer exhaust plate to the fertilizer outlets in the lower layer,mm;α is the inclination angle of the fertilizer separator,(°);Ltube is the length of the fertilizer tube,mm;Lplate is the length of the fertilizer plate,mm;Lsum is the alignment distance from the top of the fertilizer plate to the fender，mm;δ is the thickness of the fertilizer plate,mm;A is the distance between the center line of the fertilizer plate and the center line of the round tube,mm;D is the inner diameter of the round tube,mm.圖3 分肥器肥料運(yùn)動(dòng)軌跡及分肥器結(jié)構(gòu)示意圖Fig.3 Motion trajectory of fertilizer particles and structure of fertilizer separating device

肥料顆粒在分肥器中運(yùn)動(dòng)時(shí)，單個(gè)顆粒i的位移和運(yùn)動(dòng)方程為[17-20]：

(4)

(5)

式(4)表示肥料顆粒新的位移值，將新的位移值代入式(5)得到新的作用力，反復(fù)循環(huán)即為單個(gè)顆粒i在分肥器中任意時(shí)刻的運(yùn)動(dòng)狀態(tài)[21]。

排肥管和排肥板的長度決定上下兩層肥料的相對(duì)位置。排肥管過長，處于土壤回流區(qū)域，容易發(fā)生堵塞。排肥管和排肥板過短，上下兩層肥料出口的距離較近，分層效果不明顯。排肥板過長，下排肥口較小，肥料顆粒與擋泥板發(fā)生碰撞回彈，影響施肥效率。

上下兩層肥料出口的相對(duì)位置關(guān)系為：

(6)

式中：L1為上下兩層肥料出口的縱向距離，mm；L2為上下兩層肥料出口的橫向距離，mm；L3為排肥板底端到下層肥料出口的縱向距離，mm；L4為排肥板底端到下層肥料出口的橫向距離，mm；α為分肥器傾角，(°)；Ltube為排肥管長度，mm；Lplate為排肥板長度，mm；Lsum為排肥板頂端到擋泥板的對(duì)齊距離，mm；δ為排肥板厚度，mm。

根據(jù)棉花和玉米施肥深度10～20 cm的要求，確保分層效果顯著，選取Ltube=110 mm，Lplate=235 mm，Lsum=360 mm，α=50°，δ=2 mm，通過EDEM仿真檢驗(yàn)肥料分層效果。

肥料顆粒經(jīng)過均肥器后，其形態(tài)由團(tuán)聚狀變?yōu)榫紶?，分肥板中心線與圓管中心線的距離A控制著上下兩層肥料在分肥器入口通過面積，即控制上下層肥量的比例，得到上下兩層肥量比例為：

(7)

式中：Sup和Sdown分別為上下層肥料流入面積，mm2；A為分肥板中心線與圓管中心線的距離，mm；D為圓管內(nèi)徑，mm。

根據(jù)上下兩層肥量比為4∶6的設(shè)計(jì)要求，即Sdown∶Sup=4∶6，得到分肥板中心線與圓管中心線的距離A=5.868 mm，為了減小比例誤差，取分肥板的厚度為1 mm。

3 分層施肥靴參數(shù)優(yōu)化仿真試驗(yàn)

分層施肥靴的結(jié)構(gòu)決定肥料分層效果。對(duì)分層施肥靴進(jìn)行三因素三水平的仿真試驗(yàn)，根據(jù)肥料分層效果優(yōu)化最佳工作組合參數(shù)。

3.1 仿真參數(shù)設(shè)定

3.1.1全局變量參數(shù)設(shè)定

肥料顆粒、土壤和分層施肥靴的兩兩相互作用均采用Hertz-Mindlin模型[22-24]，根據(jù)相關(guān)文獻(xiàn)[25-30]，確定仿真變量參數(shù)如表1所示。

表1 仿真變量參數(shù)設(shè)置Table 1 Simulation variable parameters setting

3.1.2土壤、肥料顆粒建模

土壤顆粒的外形主要包括塊狀、核狀和柱狀3種[31]。通過前期對(duì)試驗(yàn)區(qū)域的土壤進(jìn)行分級(jí)，60%以上的土壤屬于粗砂粒(粒徑1～0.25 mm)，利用EDEM自帶的顆粒單元生成土壤顆粒模型，如圖4所示，建立1 500 mm(長)×200 mm(寬)×300 mm(高)的土槽模型，為了節(jié)省運(yùn)算時(shí)間，將組成土壤顆粒的基體半徑設(shè)置為3 mm，土壤密度設(shè)置為測(cè)量得到的1.187 g/cm3。肥料顆粒選用金屯農(nóng)化磷酸二銨(金屯嘉吉三環(huán)肥料廠)，對(duì)其三維尺寸進(jìn)行測(cè)定，得到磷酸二胺顆粒的球形率為0.965，等效直徑為3.56 mm，將肥料顆粒近似為圓形顆粒物體進(jìn)行建模[30]。

圖4 土壤顆粒EDEN仿真模型Fig.4 Simulation model of edaphic particle

3.1.3仿真模型的建立

根據(jù)設(shè)計(jì)參數(shù)使用SolidWorks 2015對(duì)分層施肥靴按1∶1比例建模，將三維模型另存為.x_t格式導(dǎo)入EDEM中。按照表1設(shè)置仿真參數(shù)，在均肥器上方建立Polygon顆粒工廠，生成速率為515 g/s，下落初速度24.5 cm/s，總仿真時(shí)間為8 s，固定時(shí)間步長為Rayleigh時(shí)間步長的25%，數(shù)據(jù)記錄間隔為0.1 s，分層施肥靴的仿真模型見圖5。

1.上層肥料；2.下層肥料；3.分層施肥靴；4.土壤顆粒1.Upper fertilizer particle；2.Lower fertilizer particle；3.Layered fertilization boot；4.Soil particle圖5 分層施肥靴仿真模型Fig.5 Simulation model of layered fertilization boot

3.2 試驗(yàn)因素與指標(biāo)

1)試驗(yàn)因素。主要影響因素有圓桿直徑和圓桿數(shù)量。分肥器的排肥管和排肥板主要影響上下兩層肥料的間隔距離，根據(jù)式(6)進(jìn)行理論計(jì)算可直接設(shè)定，故分肥器不作為影響肥料分層的因素。分層施肥靴的行進(jìn)速度影響土壤回流的快慢，直接關(guān)系到肥料的分層效果，故本研究以圓桿直徑，圓桿數(shù)量以及分層施肥靴的行進(jìn)速度作為影響因素。

2)試驗(yàn)指標(biāo)。用上下兩層肥料間隔距離的變異系數(shù)Cv作為考察指標(biāo)反映肥料的分層效果。試驗(yàn)中，通過測(cè)量分層施肥靴行進(jìn)過后同一縱截面上下兩層肥料的中心距，計(jì)算上下兩層肥料距離的變異系數(shù)。變異系數(shù)計(jì)算公式為：

(8)

式中：SSD為同一縱截面上下兩層肥料距離的標(biāo)準(zhǔn)差，cm；Smean為同一縱截面上下兩層肥料距離的平均值，cm。

3.3 試驗(yàn)設(shè)計(jì)

根據(jù)前期相關(guān)試驗(yàn)得到各變量的變化范圍，選取圓桿直徑4 mm、圓桿數(shù)量6根和行進(jìn)速度 6 km/h 為零水平，以上下兩層肥料距離的變異系數(shù)為評(píng)價(jià)指標(biāo)(表2)。采用響應(yīng)面法分析3個(gè)因子對(duì)響應(yīng)值(上下兩層肥料距離的變異系數(shù))的影響程度，得到分層施肥靴的最佳組合參數(shù)。

表2 響應(yīng)面因素水平表Table 2 Factors and levels of RSM analysis

3.4 結(jié)果分析

3.4.1回歸模型的建立與檢驗(yàn)

采用Design-Expert.V10.0.3的Box-Behnken模塊進(jìn)行三因素三水平的試驗(yàn)設(shè)計(jì)[33]，試驗(yàn)結(jié)果見表3。根據(jù)顯著性檢驗(yàn)，模型的F值為33.87，P<0.01，表眀該模型極顯著。變異系數(shù)模型的決定系數(shù)R2=0.9775，證明該模型擬合程度良好。上下兩層肥料距離變異系數(shù)的回歸方程為：

表3 仿真試驗(yàn)方案及結(jié)果Table 3 Simulation test plan and results

Y=4.60-0.07X1-0.33X2-0.19X3+
0.27X1X2+0.04X1X3+0.50X2X3+
0.85X12+0.85X22+1.11X32

(9)

式中：X1、X2和X3分別為圓桿直徑、圓桿數(shù)量和行進(jìn)速度的水平編碼；Y為上下兩層肥料距離變異系數(shù)，%。

由回歸方程系數(shù)的絕對(duì)值大小得到各因素影響程度主次順序?yàn)閳A桿數(shù)量>施肥靴行進(jìn)速度>圓桿直徑。對(duì)式(9)所表示的回歸方程進(jìn)行方差分析，分析結(jié)果見表4。

3.4.2兩因素對(duì)變異系數(shù)的影響

由表4可知，圓桿直徑和圓桿數(shù)量的交互作用對(duì)上下兩層肥料距離變異系數(shù)的影響顯著，而圓桿數(shù)量和分層施肥靴行進(jìn)速度的交互作用對(duì)上下兩層肥料距離變異系數(shù)的影響極顯著，因此對(duì)圓桿數(shù)量和分層施肥靴的行進(jìn)速度對(duì)上下兩層肥料距離變異系數(shù)的影響進(jìn)行響應(yīng)曲面分析，結(jié)果見圖6。

表4 上下兩層肥料距離的變異系數(shù)回歸模型方差分析Table 4 Analysis of variance for regression model of coefficient of variation of fertilizer distance between upper and lower layers

圖6 施肥靴行進(jìn)速度和均肥器圓桿數(shù)量對(duì)變異系數(shù)的影響Fig.6 Effects of travel speed of layered fertilization boot and pole number of fertilizer equalizer on coefficient of variation

可見，圓桿數(shù)量和分層施肥靴的行進(jìn)速度兩因素的響應(yīng)面為開口向上的拋物曲面。圓桿數(shù)量在低水平時(shí)，分層施肥靴的行進(jìn)速度對(duì)上下兩層肥料距離的變異系數(shù)的影響極顯著，表現(xiàn)為響應(yīng)面拋物曲線比較陡，說明在圓桿數(shù)量較少時(shí)，增加分層深施機(jī)構(gòu)的行進(jìn)速度能明顯降低上下兩層肥料距離的變異系數(shù)。而在圓桿數(shù)量較多時(shí)，分層施肥靴的行進(jìn)速度對(duì)上下兩層肥料距離的變異系數(shù)的影響不是很明顯，表現(xiàn)為響應(yīng)拋物曲線比較平緩，這表明圓桿數(shù)量和分層施肥靴的行進(jìn)速度存在一定的交互作用，同時(shí)根據(jù)方差分析可知二者的交互作用為極顯著[34]。

3.4.3分層施肥靴參數(shù)優(yōu)化

為確保施肥靴具有更好的分層效果，本研究以上下層肥料距離的變異系數(shù)為分層目標(biāo)，進(jìn)行分層施肥靴工作參數(shù)和結(jié)構(gòu)參數(shù)優(yōu)化，利用Design-Expert.V10.0.3軟件中的Optimization-Numerical模塊進(jìn)行優(yōu)化求解，其目標(biāo)函數(shù)與約束條件為：

(10)

根據(jù)上下兩層肥料距離的變異系數(shù)越小越好的原則，優(yōu)化后得到影響因素最佳區(qū)域，通過Design-Expert.V10.0.3選取最佳參數(shù)組合為：圓桿直徑4 mm，圓桿數(shù)量6根，分層施肥靴的行進(jìn)速度6.04 km/h，此條件下上下兩層肥料距離的變異系數(shù)為4.57%。

4 驗(yàn)證試驗(yàn)

通過理論計(jì)算分肥器中肥料顆粒的流通面積，將上下層肥量比固定為4∶6，但經(jīng)過均肥器的肥料顆粒必須為均布狀才能實(shí)現(xiàn)按比例分層，故對(duì)均肥器進(jìn)行仿真試驗(yàn)，觀察肥料顆粒的分布狀態(tài)，利用臺(tái)架試驗(yàn)實(shí)測(cè)上下層分肥比例，最后由田間試驗(yàn)觀察肥料顆粒在土壤中的分層效果。

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

通過臺(tái)架試驗(yàn)對(duì)上下排肥口的肥量進(jìn)行實(shí)測(cè)驗(yàn)證。試驗(yàn)地點(diǎn)為新疆科神農(nóng)業(yè)裝備科技開發(fā)股份有限公司的裝配工廠；試驗(yàn)材料選用金屯農(nóng)化磷酸二銨；試驗(yàn)裝置見圖7，主要包括直流電機(jī)(1016Z1B型，豫鑫電機(jī)公司)、12 V直流電機(jī)調(diào)速控制器(12-2*500全功能型，豫鑫電機(jī)公司)、分層施肥靴(圓桿數(shù)量為6根，圓桿直徑為4 mm)。

1.肥箱；2.調(diào)肥手輪；3.直流電機(jī)；4.分層施肥靴；5.12V直流電機(jī)調(diào)速控制器；6.上層肥料收納袋；7.下層肥料收納袋1.Fertilizer can；2.Fertilizer handwheel；3.DC motor；4.Layered fertilization boot；5.12 V DC motor speed controller；6.Upper fertilizer storage bag；7.Lower fertilizer storage bag圖7 上下層肥量比測(cè)定臺(tái)架試驗(yàn)裝置Fig.7 A bench test device for the determination of the fertilizer ratio of the upper and lower layers

使用直流電機(jī)進(jìn)行分層施肥靴行進(jìn)速度的模擬；用調(diào)肥手輪進(jìn)行施肥量的調(diào)整；使用接肥袋分別收集上下兩層肥料并稱重計(jì)算出其占比。試驗(yàn)按施肥量(20、25、30 kg)分成 3個(gè)大組，每個(gè)大組又根據(jù)行進(jìn)速度(5、6、7 km/h)各分成3個(gè)小組，每個(gè)小組進(jìn)行10次試驗(yàn)取其平均值，試驗(yàn)結(jié)果見圖8。試驗(yàn)結(jié)果表明：根據(jù)理論計(jì)算獲取的均肥器結(jié)構(gòu)參數(shù)可以保證肥料顆粒順暢流動(dòng)，施肥量和分層施肥靴的行進(jìn)速度對(duì)上下層肥量比的影響甚微，分層施肥靴滿足上下層肥量比4∶6的設(shè)計(jì)要求。

圖8 3種速度不同施肥量下施肥靴上下層肥量比Fig.8 The ratio of the upper and lower layers of the boots with different fertilization rates

4.2 田間試驗(yàn)

為了驗(yàn)證離散元仿真試驗(yàn)所得最優(yōu)組合參數(shù)及分層施肥靴的工作性能，在不同施肥量條件下，比較仿真試驗(yàn)與田間試驗(yàn)上下兩層肥料距離變異系數(shù)的相對(duì)誤差。試驗(yàn)地點(diǎn)為伊寧縣青年農(nóng)場，試驗(yàn)區(qū)域土壤粒級(jí)為粗砂粒，密度1.187 g/cm3，含水率14.726%，堅(jiān)實(shí)度0.57 MPa。試驗(yàn)材料為金屯農(nóng)化磷酸二銨；試驗(yàn)裝置主要包括雷沃歐豹TG125型拖拉機(jī)、安裝分層施肥靴的分層施肥機(jī)等(圖9)。根據(jù)仿真試驗(yàn)結(jié)果，分層施肥靴圓桿數(shù)量選用6根，圓桿直徑4 mm，機(jī)具行進(jìn)速度設(shè)為6.04 km/h。

圖9 分層施肥機(jī)與分層施肥靴Fig.9 Layered fertilizer applicator an layered fertilizer boots

田間試驗(yàn)和仿真試驗(yàn)肥料分層效果比對(duì)(圖10)顯示，2種試驗(yàn)均實(shí)現(xiàn)了肥料顆粒的上下分層且效果明顯，滿足設(shè)計(jì)要求。施肥量20 kg時(shí)肥料顆粒在土壤中的位置分布情況見圖11，仿真試驗(yàn)和田間試驗(yàn)都存在一定的波動(dòng)性，波動(dòng)都不超過預(yù)設(shè)的幅度范圍，上層肥料位于距地表10～15 cm，下層肥料位于距地表15～20 cm，上下兩層肥料間隔在5 cm左右，肥料分層效果顯著。

1.上層肥料；2.土壤顆粒；3.下層肥料1.Upper fertilizer particle；2.Soil particle；3.Lower fertilizer particle圖10 田間試驗(yàn)與仿真試驗(yàn)肥料分層效果對(duì)比Fig.10 Comparison of fertilizer stratification effect between field test and simulation test

圖11 施肥量20 kg肥料顆粒分布位置Fig.11 Fertilizer distribution locations of 20 kg fertilizer granules

測(cè)量不同施肥量條件下田間試驗(yàn)和仿真試驗(yàn)中上下兩層肥料間隔距離的變異系數(shù)見表5。可以看出變異系數(shù)的仿真結(jié)果比試驗(yàn)結(jié)果略大，相對(duì)誤差不超過4%，而且具有相同趨向，施肥量越大，肥料顆粒越多，分層更為明顯，故上下兩層肥料間隔距離的變異系數(shù)越小。證明離散元仿真可以很好地輔助關(guān)鍵部件的設(shè)計(jì)與優(yōu)化。

表5 仿真試驗(yàn)與田間試驗(yàn)上下層肥料距離變異系數(shù)相對(duì)誤差Table 5 Relative error of the coefficient of variation of the distance between the top and bottom layers of the simulation test and the field test

5 結(jié)論及討論

1)本研究設(shè)計(jì)了一種分層施肥靴，以上下兩層肥料間隔距離的變異系數(shù)為分層目標(biāo)，對(duì)EDEM仿真試驗(yàn)數(shù)據(jù)通過Design-Expert.V10.0.3軟件處理及優(yōu)化，分層施肥靴最優(yōu)組合參數(shù)為：圓桿數(shù)量6根，圓桿直徑4 mm，分層施肥靴行進(jìn)速度6.04 km/h，此時(shí)變異系數(shù)為4.57%。

2)選用最優(yōu)組合參數(shù)的分層施肥靴進(jìn)行臺(tái)架試驗(yàn)和田間試驗(yàn)，并與仿真試驗(yàn)進(jìn)行對(duì)比，結(jié)果表明，在不同施肥量條件下，分層施肥靴均有較好的肥料分層效果，仿真值與試驗(yàn)值的相對(duì)誤差不超過4%，分層效果基本吻合，施肥效果穩(wěn)定，說明利用離散元仿真試驗(yàn)對(duì)分層施肥裝置進(jìn)行優(yōu)化具有可行性。

由于化肥種類繁多，物料特性存在一定差異，因此本研究中的分層施肥靴試驗(yàn)裝置的試驗(yàn)對(duì)象為物料特性與磷酸二銨類似的化肥，以實(shí)現(xiàn)分層施肥的目的。機(jī)具在作業(yè)過程中，經(jīng)過土壤條件較差的區(qū)域時(shí)，在下排肥口有土壤粘黏，并且施肥靴行進(jìn)阻力略有增加，分析原因有：施肥靴機(jī)構(gòu)存在不足，需進(jìn)一步優(yōu)化；施肥靴材料方面的原因；土壤黏度較大，前期整地效果不好。后續(xù)將進(jìn)一步對(duì)影響施肥靴土壤粘黏和減阻性能的因素進(jìn)行研究。