祝清震，武廣偉，陳立平，趙春江，孟志軍

?

槽輪結(jié)構(gòu)參數(shù)對直槽輪式排肥器排肥性能的影響

祝清震1,2，武廣偉1※，陳立平3，趙春江1,2，孟志軍3

（1. 北京農(nóng)業(yè)信息技術(shù)研究中心，北京 100097；2. 西北農(nóng)林科技大學(xué)機械與電子工程學(xué)院，楊凌 712100；3. 北京農(nóng)業(yè)智能裝備技術(shù)研究中心，北京 100097）

外槽輪排肥器是實施精準變量施肥過程中的關(guān)鍵控制載體，為了探尋外槽輪的凹槽截面形狀和參數(shù)設(shè)置對排肥器排肥性能的影響規(guī)律，該文運用離散元仿真技術(shù)、3D快速成型技術(shù)以及正交試驗，以直槽輪式排肥器作為研究對象，分析了槽輪半徑、凹槽數(shù)目、有效工作長度以及凹槽截面形狀4個因素對排肥性能的影響。選用排肥均勻度變異系數(shù)作為評價指標，設(shè)置排肥器的前進速度0.5 m/s，工作高度500 mm，槽輪轉(zhuǎn)速20 r/min，對4個影響因素進行交互判斷試驗，試驗結(jié)果表明，槽輪半徑與凹槽數(shù)目之間存在一定的交互作用。根據(jù)交互作用判別結(jié)果選用L18（37）標準正交表進行仿真試驗，試驗結(jié)果表明，影響排肥性能的主次因素依次為凹槽數(shù)目>槽輪半徑>有效工作長度>槽輪半徑×凹槽數(shù)目>凹槽截面形狀，其中凹槽數(shù)目對排肥性能具有非常顯著影響，槽輪半徑對排肥性能具有顯著影響，最佳結(jié)構(gòu)參數(shù)組合是槽輪半徑30 mm，凹槽數(shù)目7個，有效工作長度20 mm和凹槽截面形狀為圓弧形，此水平組合下的排肥均勻度變異系數(shù)為1.75%，為了檢驗結(jié)構(gòu)參數(shù)對直槽輪排肥性能影響的仿真試驗結(jié)果，進行最優(yōu)結(jié)構(gòu)與參數(shù)組合下的直槽輪排肥器臺架試驗，試驗結(jié)果顯示，排肥均勻度變異系數(shù)為5.48%，仿真試驗結(jié)果和實測結(jié)果基本吻合，說明應(yīng)用離散元仿真技術(shù)開展排肥性能的影響因素分析研究是可行的。該研究結(jié)果可為直槽輪排肥器的槽輪結(jié)構(gòu)參數(shù)優(yōu)化及提升其排肥性能提供參考。

農(nóng)業(yè)機械；離散元法；試驗；凹槽截面；結(jié)構(gòu)參數(shù)；排肥性能；排肥器

0 引 言

施用化肥是農(nóng)作物生產(chǎn)管理過程中的一個重要環(huán)節(jié)，已有研究表明，進行化肥的精準施用對提高肥料的利用率、減少化肥的施用量具有重要意義[1-3]。為了提高肥料的精準施用效果，科研人員對精量排肥控制系統(tǒng)進行了大量研究[4-5]，左興健[6-7]等開發(fā)了基于作業(yè)機具實時速度的精準排肥控制系統(tǒng)，孟志軍等設(shè)計研發(fā)了基于處方圖的精準變量施肥控制系統(tǒng)[8-9]，汪小旵等結(jié)合近地光譜技術(shù)，設(shè)計了冬小麥實時變量追肥控制系統(tǒng)[10-11]，苑進等[12-13]對精量排肥的控制算法進行優(yōu)化研究，以上研究都是通過實時改變外槽輪排肥器轉(zhuǎn)速或者外槽輪有效工作長度達到精準變量排肥的目的，作為重要的實施載體，外槽輪排肥器在精量排肥過程中發(fā)揮了重要作用，然而外槽輪排肥器工作過程中存在脈動性大、排肥穩(wěn)定性差的問題，另外一方面中國生產(chǎn)的外槽輪排肥器品類繁多，排肥效果差異性較大，開展與其結(jié)構(gòu)參數(shù)相關(guān)的排肥性能影響分析研究非常必要。

為了提高外槽輪排肥器的作業(yè)效果，科研人員對其結(jié)構(gòu)參數(shù)已經(jīng)進行了大量的研究，施印炎等將排肥舌替代原有的塑料毛刷，以避免排肥過程中出現(xiàn)肥料顆粒泄露現(xiàn)象[14-15]，張濤等應(yīng)用離散元法對外槽輪的有效工作長度進行了仿真研究[16]，楊洲等應(yīng)用離散元仿真技術(shù)和3D打印成型技術(shù)對外槽輪排肥器的凹槽半徑和螺旋升角等結(jié)構(gòu)參數(shù)進行了優(yōu)化研究[17]，汪博濤等綜合應(yīng)用離散元法、二次正交回歸試驗的方法對外槽輪的有效工作長度、排肥舌開度進行了研究[18]，潘世強對外槽輪的槽數(shù)、有效工作長度、排肥舌開度及其制造材料進行了試驗研究[19]?，F(xiàn)階段關(guān)于外槽輪排肥器結(jié)構(gòu)參數(shù)方面的研究大多局限在槽輪工作有效長度、排肥舌開度、槽數(shù)等，影響外槽輪排肥器作業(yè)性能的結(jié)構(gòu)參數(shù)眾多，還包括槽輪的半徑、凹槽的截面形式等，關(guān)于這些關(guān)鍵結(jié)構(gòu)參數(shù)優(yōu)化的研究還比較欠缺，另外一方面，現(xiàn)有研究還沒有考慮到各結(jié)構(gòu)參數(shù)之間的交互作用，因此在現(xiàn)有研究的基礎(chǔ)上，開展外槽輪不同結(jié)構(gòu)參數(shù)對排肥性能的影響規(guī)律研究十分必要。

針對現(xiàn)有外槽輪排肥器品類繁多，槽輪結(jié)構(gòu)參數(shù)的差異將影響排肥器的工作性能，本文擬開展外槽輪凹槽截面形狀和參數(shù)設(shè)置與排肥器排肥性能影響關(guān)系的分析，以直槽輪排肥器作為研究對象，在分析直槽輪排肥器關(guān)鍵結(jié)構(gòu)參數(shù)組成及各因素的交互作用判斷的基礎(chǔ)上，設(shè)計了6因素3水平的正交試驗，應(yīng)用離散元仿真技術(shù)對不同結(jié)構(gòu)參數(shù)下的直槽輪排肥器的排肥效果進行極差和方差分析，并應(yīng)用3D快速成型技術(shù)將最優(yōu)結(jié)構(gòu)參數(shù)組合下的直槽輪排肥器制造出來，進行臺架驗證試驗，以期為直槽輪排肥器的槽輪結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計和提升其排肥性能提供參考。

1 直槽輪式排肥器工作原理及結(jié)構(gòu)參數(shù)分析

1.1 直槽輪式排肥器結(jié)構(gòu)組成

直槽輪式排肥器作為廣泛應(yīng)用的排肥器種類之一，其結(jié)構(gòu)主要包括直槽輪、肥盒、毛刷和排肥舌等，如圖1所示，其中直槽輪是排肥器完成排肥作業(yè)的關(guān)鍵部件。

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

1.2 直槽輪式排肥器工作原理

直槽輪式排肥器工作時，肥料首先依靠自身重力充滿槽輪凹槽，隨著槽輪的旋轉(zhuǎn)，被強制排出，處于槽輪凹槽外側(cè)的肥料顆粒也會在槽輪外圓的撥動和肥料顆粒間的摩擦力作用下被帶動起來，帶動層的線速度從槽輪圓周逐漸向外遞減直至靜止層，由槽輪強制帶出和帶動層帶出的肥料從排肥舌上掉入排肥管，完成整個排肥過程。根據(jù)排肥工作原理，直槽輪式排肥器的槽輪旋轉(zhuǎn)1周的排肥量可以按式（1）進行計算[13]。

式中1表示被槽輪凹槽強制排出的肥料質(zhì)量，g/r；2表示帶動層排出的肥料質(zhì)量，g/r；表示槽輪旋轉(zhuǎn)1周的肥料排出質(zhì)量，g/r；表示肥料顆粒的密度，g/cm3；表示凹槽內(nèi)肥料的充滿系數(shù)；表示凹槽數(shù)目；表示單個凹槽的截面積，mm2；表示槽輪的有效長度，mm；表示槽輪的半徑，mm；表示肥料顆粒的帶動層系數(shù)，mm。

1.3 直槽輪關(guān)鍵結(jié)構(gòu)參數(shù)分析

直槽輪凹槽截面形狀要便于肥料的填充和排出，現(xiàn)有凹槽端面主要有圓弧形凹槽、直角凹槽和錐形圓弧凹槽3種形狀，3種凹槽截面形狀的剖視圖如圖2所示。

注：O表示槽輪的軸心；O1表示圓弧形凹槽的圓心；O2表示錐形圓弧凹槽的圓心；α表示相鄰凹槽之間的節(jié)距角，rad；β表示單個凹槽在槽輪橫截面上的跨度，rad；R表示槽輪的半徑，mm；r1表示圓弧形凹槽的半徑，mm；r2表示錐形圓弧凹槽的根圓半徑，mm；r3表示錐形圓弧凹槽的半徑，mm。

為了便于對單個凹槽截面積進行數(shù)字化分析，這里相鄰凹槽之間的節(jié)距角、單個凹槽在槽輪橫截面上的跨度、凹槽圓弧半徑和單個凹槽的截面積之間應(yīng)滿足

式中表示相鄰凹槽之間的節(jié)距角，rad；表示單個凹槽在槽輪橫截面上的跨度，rad；1表示圓弧形凹槽的半徑，mm；2表示錐形圓弧凹槽的根圓半徑，mm；3表示錐形圓弧凹槽的半徑，mm；1表示單個圓弧形凹槽的截面積，mm2；2表示單個直角凹槽的截面積，mm2；3表示單個錐形圓弧凹槽的截面積，mm2。

對于特定的排施對象，由式（1）可知，影響直槽輪排肥器作業(yè)效果的槽輪結(jié)構(gòu)因素主要有凹槽數(shù)目、槽輪有效工作長度、槽輪的半徑和單個凹槽的截面積等。由式（2）可知，單個凹槽的截面積大小可以用凹槽數(shù)目和槽輪半徑表示；當凹槽的截面積相同時，其截面形狀卻不一定相同，圓弧形凹槽、直角凹槽和錐形圓弧凹槽作為廣泛采用的3種凹槽截面形狀，探索這3種凹槽截面形狀對排肥器工作性能的影響也具有一定意義。綜上分析，直槽輪的關(guān)鍵結(jié)構(gòu)參數(shù)主要有凹槽數(shù)目、槽輪有效長度、槽輪的半徑和凹槽截面形狀。

2 離散元仿真平臺搭建

直槽輪式排肥器的作業(yè)過程中伴隨著大量的顆粒運動，在肥料顆粒被排出的過程中存在重力、摩擦力、槽輪給予顆粒的強制力等多種作用力的綜合作用。離散元法（discrete element method，DEM）作為有效的散粒體運動分析方法，已經(jīng)在農(nóng)業(yè)工程領(lǐng)域被廣泛應(yīng)用[20-22]。采用離散元分析軟件EDEM對排肥器的工作過程進行仿真分析，以確定直槽輪各結(jié)構(gòu)參數(shù)對排肥性能影響的主次，并確定所研究結(jié)構(gòu)參數(shù)的最優(yōu)水平組合，可以降低試驗成本，提高研究工作效率。在進行離散元仿真試驗前，需要搭建直槽輪排肥過程的離散元仿真平臺，鑒于前文已有圖1所示的直槽輪式排肥器結(jié)構(gòu)模型，這里主要對肥料顆粒模型和離散元接觸模型進行介紹。

2.1 肥料顆粒離散元模型

山東肥業(yè)有限公司生產(chǎn)的藍精靈復(fù)合肥被中國黃淮海地區(qū)廣泛使用，因此本文選用該肥料作為化肥顆粒離散元建模的試驗材料，參照化肥顆粒三維建模方法[23-24]，從肥料堆中隨機取樣200粒，使用游標卡尺（桂林廣陸數(shù)字測控股份有限公司，K15G278418數(shù)顯卡尺）對化肥顆粒的三維尺寸（長、寬、高）進行測量，長度、寬度和高度分別為3.09、3.06和2.58 mm，按照式（3）和（4）計算化肥顆粒的等效直徑和球形率，分別為2.90 mm和93.86 %。

式中表示肥料顆粒的等效直徑，mm；表示肥料顆粒的球形率，%；、、分別表示肥料顆粒的長、寬、高，mm。

經(jīng)過測量和計算，化肥顆粒樣本的球形率已經(jīng)大于90%，說明該化肥顆粒樣本的個體具有較高的球形率分布特征，因此選擇直徑為2.90 mm的球作為肥料顆粒的三維離散元模型是合適的。

應(yīng)用OriginPro 2017軟件對樣本肥料顆粒的群體分布規(guī)律進行統(tǒng)計分析，繪制柱狀圖，如圖3所示，并使用式（5）所示的正態(tài)分布函數(shù)對樣本分布規(guī)律數(shù)據(jù)進行擬合。

式中表示因變量，表示自變量，0、、0和表示正態(tài)分布函數(shù)的系數(shù)，擬合后0=1.42，=24.29，0=2.84，=0.28，其中擬合優(yōu)度2=0.98，可以認為化肥顆粒的等效直徑服從正態(tài)分布規(guī)律，所以設(shè)置顆粒工廠的肥料顆粒生成規(guī)律服從正態(tài)分布。

2.2 離散元接觸模型

EDEM軟件內(nèi)置了多種接觸模型，其中Hertz-Mindlin (no slip)模型是EDEM軟件中默認接觸模型，在力的計算方面精確且高效[25]。鑒于肥料的排出過程不涉及肥料顆粒之間的粘結(jié)作用，因此本研究采用Hertz-Mindlin (no slip)模型作為肥料顆粒與肥料顆粒，肥料顆粒與排肥器（肥箱、排肥盒、直槽輪、排肥舌和排肥管）、地面之間的接觸模型，這里排肥器的肥箱、排肥盒、直槽輪、排肥舌和排肥管均采用PLA材料的屬性。本研究采用試驗測定和查閱相關(guān)文獻相結(jié)合的方法，確定肥料顆粒、排肥器和地面相關(guān)的材料和接觸力學(xué)參數(shù)，如表1所示。

圖3 肥料顆粒的等效直徑分布柱狀圖

表1 離散元模型的材料及接觸力學(xué)參數(shù)

3 仿真正交試驗

3.1 試驗方法

將排肥器結(jié)構(gòu)模型導(dǎo)入離散元仿真軟件EDEM 2.7.0，為了觀察肥料在水平地表的分布狀況，在距離排肥器肥盒下端高度為500 mm的地方設(shè)置一個長1200 mm，寬250 mm的平面，用于模擬地面，其中該模擬地面的中軸線與排肥系統(tǒng)的中軸線位于同一平面。研究發(fā)現(xiàn)，直槽輪排肥器在轉(zhuǎn)速15~30 r/min時，具有較好的排肥性能[16]，所以本文選擇直槽輪的工作轉(zhuǎn)速為20 r/min；參考文獻[24]，本文設(shè)置排肥器的前進速度為0.5 m/s，為了降低排肥器運動的復(fù)雜性，這里設(shè)置地面以-0.5 m/s的速度向后方運動，用于模擬排肥裝置的前進速度，并參照表1設(shè)置相關(guān)的材料物理和接觸模型參數(shù)。

在肥箱上端設(shè)置一個顆粒工廠，使其在仿真過程中不斷生產(chǎn)化肥顆粒，顆粒工廠每秒產(chǎn)生3 000個化肥顆粒，生成總量為10 000個，當顆粒數(shù)達到一定程度，分別設(shè)置排肥輪轉(zhuǎn)動和地面運動，進行仿真試驗，仿真試驗過程如圖4所示。

圖4 排肥器工作過程仿真

3.2 排肥性能評價方法

不同的評價指標使得影響排肥器排肥性能的因素顯著性也不相同，為了準確評價在EDEM仿真試驗過程中，不同結(jié)構(gòu)參數(shù)直槽輪排肥器的排肥效果，需設(shè)定其評價標準。本研究參照《JB/T9783-2013，播種機外槽輪排種器》規(guī)定的試驗方法，以排肥均勻度變異系數(shù)作為排肥器排肥性能的評價指標。選取模擬地面中間1 000 mm區(qū)域作為排肥效果的取樣區(qū)域，在模擬地面上設(shè)置Grid Bin Group，將中間1 000 mm區(qū)域橫向均分為5個網(wǎng)格單元（每個單元網(wǎng)格尺寸200 mm×280 mm×50 mm），并進行編號，如圖5所示。

注：1~5表示數(shù)據(jù)采集網(wǎng)格單元的編號。

仿真結(jié)束后，統(tǒng)計每一個網(wǎng)格單元（grid bin group）內(nèi)所有肥料顆粒的總質(zhì)量，設(shè)置第個取樣網(wǎng)格單元內(nèi)肥料顆粒的質(zhì)量為m。利用式（6），求解5個網(wǎng)格單元內(nèi)肥料顆粒的平均質(zhì)量。

根據(jù)式（7），計算取樣區(qū)域內(nèi)所有網(wǎng)格單元內(nèi)總化肥顆粒質(zhì)量的標準差。

式中表示取樣網(wǎng)格單元內(nèi)化肥顆?？傎|(zhì)量的標準差，g。

根據(jù)式（8），計算排肥器的排肥均勻度變異系數(shù)。

式中表示排肥器的排肥均勻度變異系數(shù)，%。

排肥均勻度變異系數(shù)作為評價排肥器排肥穩(wěn)定性和均勻度效果的參量，數(shù)值越小，說明排肥器的排肥穩(wěn)定性和均勻性越好。利用排肥均勻度變異系數(shù)可以用于對比分析不同直槽輪結(jié)構(gòu)參數(shù)形式下的排肥器排肥性能。

3.3 試驗因素水平選取

根據(jù)直槽輪排肥器工作原理，及影響排肥效果的槽輪結(jié)構(gòu)參數(shù)影響因素分析結(jié)果，試驗選取影響排肥器排肥性能的直槽輪結(jié)構(gòu)參數(shù)因素為槽輪半徑、凹槽數(shù)目、槽輪有效長度和凹槽截面形狀。通過對市場調(diào)研和文獻檢索，直槽輪排肥器的凹槽數(shù)目一般在5~8之間[19]，槽輪半徑、有效工作長度的最大值分別不超過40和60 mm[16]，并結(jié)合實際應(yīng)用需求，本文為4個因素分別選取3個水平，如表2所示。

表2 正交試驗因素水平

3.4 交互作用判別

在進行正交試驗設(shè)計過程中，必須考慮到所選因素之間是否存在一定的交互作用，其中2個因素間的交互作用稱為一級交互作用，3個或3個因素以上的因素間交互，稱為高級交互作用（在絕大多數(shù)的實際問題中，高級交互作用可以忽略[29]）。所以有必要對選取的4個因素進行交互作用的判別，參考文獻[30]分別進行槽輪半徑與槽輪有效長度、槽輪半徑與凹槽數(shù)目、槽輪半徑與凹槽截面形狀、凹槽數(shù)目與槽輪有效長度、凹槽數(shù)目與凹槽截面形狀、槽輪有效長度與凹槽截面形狀之間的交互判別試驗。以槽輪半徑和槽輪有效長度2個因素之間的交互作用判別試驗方法為例，保持凹槽數(shù)目、凹槽截面形狀不變，槽輪半徑分別取25，30和35 mm，有效長度分別取20，35和50 mm，槽輪半徑和有效長度2因素之間共計9種不同的試驗組合，每種組合各做1次仿真試驗，以排肥均勻度變異系數(shù)作為評價指標，其余5組判別試驗方法與此類似，4個因素之間的交互判別試驗結(jié)果如圖6所示。

因素之間是否存在交互作用的判別方法參考文獻[30]，由于試驗誤差的存在，如果3條曲線變化趨勢近似，也可以認為2個因素之間無交互作用，或者交互作用可以忽略不計，圖6a、6c、6d、6e和6f中的3條曲線的變化趨勢基本一致，說明槽輪半徑與有效工作長度、槽輪半徑與凹槽截面形狀、凹槽數(shù)目與有效工作長度、凹槽數(shù)目與凹槽截面形狀、槽輪有效工作長度與凹槽截面形狀之間不存在明顯的交互作用；在圖6b中，3條曲線變化趨勢不一致，這是因素之間存在相交作用的一種表現(xiàn)形式，說明槽輪半徑與凹槽數(shù)目之間可能存在一定的交互作用。

3.5 正交試驗表頭設(shè)計

由于所研究影響排肥性能的槽輪結(jié)構(gòu)參數(shù)的因素之間存在一定交互作用，必須在正交表中必須設(shè)置相應(yīng)的列，稱為交互作用列。由于槽輪半徑與凹槽數(shù)目之間可能存在一定的交互作用的狀況下，同時3水平因素之間的交互作用需要占用2列，所以本研究至少應(yīng)該按照6因素3水平的情況來選擇正交表，這里選擇滿足這一條件的最小正交試驗表L18（37），設(shè)置因素槽輪半徑、凹槽數(shù)目、有效工作長度和凹槽截面形狀分別用A、B、C和D進行表示，相關(guān)因素列設(shè)置如表3所示。

3.6 仿真結(jié)果分析與討論

3.6.1 極差分析

根據(jù)設(shè)計的L18（37）表頭進行正交試驗，共計18個試驗組合，得到的試驗結(jié)果如表4所示。由排肥性能評價標準可知，排肥均勻度變異系數(shù)越小，排肥器的排肥效果越好，因此影響排肥器排肥性能的主次因素順序：B，A，C，A×B，D。

圖6 因素之間交互作用仿真試驗結(jié)果

表3 L18（37）表頭設(shè)計

注：第7列為空列，未在表中列出，A表示槽輪半徑，B表示凹槽數(shù)目，C表示槽輪有效工作長度，D表示凹槽截面形狀。

Note: The seventh column is null column, which is not listed in the table. A is radius of flute wheel, B is groove number, C is effective length of flute wheel, D is section of groove.

由表4可知，交互作用A×B對直槽輪排肥器排肥性能的影響程度小于因素A和B，所以在設(shè)置的3個水平下，最優(yōu)的結(jié)構(gòu)參數(shù)的水平組合為A2B2C1D1，即槽輪半徑為30 mm，槽數(shù)為7，有效長度為20 mm和凹槽截面形狀為圓弧形時，直槽輪排肥器的排肥均勻度變異系數(shù)最小，具有較好的排肥工作性能。

3.6.2 方差分析

為了精確估計直槽輪的4個結(jié)構(gòu)參數(shù)對排肥性能影響的重要程度，特別是考慮到因素之間的交互作用，有必要對仿真試驗結(jié)果進行方差分析。方差分析結(jié)果如表5所示，通過查表0.05（2,10）=4.10，0.01（2,10）=7.56，所以因素B對試驗結(jié)果有非常顯著的影響，因素A對試驗結(jié)果有顯著的影響，因素A×B、C和D對試驗結(jié)果沒有顯著影響。

3.6.3 討 論

結(jié)合表4與表5的分析結(jié)果可知，凹槽截面形狀對排肥器的排肥變異系數(shù)沒有顯著影響，且是4個因素中對排肥變異系數(shù)影響最小的因素，出現(xiàn)這種試驗結(jié)果的原因應(yīng)該與前文式（2）的單個凹槽截面積數(shù)字化建模過程有關(guān)，當3種凹槽截面形狀的直槽輪凹槽截面積相等時，由式（1）可知，直槽輪旋轉(zhuǎn)1周的排肥量也是相等的，然而本文建立的排肥器排肥性能評價標準是在式（6）的排肥量基礎(chǔ)上開展，反映在表4中，凹槽截面形狀的極差最小，即相較于其他3個因素，凹槽截面形狀對排肥變異系數(shù)影響最??；反映在表5中，凹槽截面形狀的平方和遠小于校正前殘差的平方和，即對排肥變異系數(shù)的影響不顯著。

表4 正交試驗結(jié)果

注：表示排肥均勻度變異系數(shù)，%。

Note:is coefficient of variation, %.

注：* 表示影響顯著（F0.01>>F0.05）；** 表示影響非常顯著（>F0.01）；表示校正前的殘差；△表示校正后的殘差。

Note: * indicate that this impact is significant (F0.01>>F0.05); ** indicate that this impact is very significant (>F0.01);is residual error before correction;△is residual error after correction.

4 仿真試驗結(jié)果驗證

根據(jù)直槽輪排肥器結(jié)構(gòu)參數(shù)對排肥性能的仿真試驗結(jié)果，構(gòu)建最優(yōu)結(jié)構(gòu)參數(shù)水平組合下的直槽輪排肥器結(jié)構(gòu)模型，即槽輪半徑為30 mm、槽數(shù)為7、有效長度為20 mm和凹槽截面形狀為圓弧形的直槽輪排肥器，進行仿真排肥性能試驗，其排肥均勻度變異系數(shù)為1.75%，相對于18組仿真正交試驗結(jié)果，該結(jié)構(gòu)參數(shù)水平組合的直槽輪排肥器具有最優(yōu)的排肥性能。應(yīng)用3D打印機（上海派恩科技有限公司，N2 Plus FFF 3D Printer）加工該結(jié)構(gòu)參數(shù)組合下的直槽輪及排肥器其他結(jié)構(gòu)部件（PLA材質(zhì)），在北京市昌平區(qū)小湯山國家精準農(nóng)業(yè)研究示范基地，開展直槽輪排肥器最佳結(jié)構(gòu)參數(shù)組合下的臺架驗證試驗，如圖7所示。參照仿真試驗的試驗條件，設(shè)置臺架試驗過程中直槽輪的轉(zhuǎn)速為20 r/min，傳送帶的前進速度為0.5 m/s，排肥管的高度為500 mm，肥料選用山東肥業(yè)有限公司生產(chǎn)的藍精靈復(fù)合肥，試驗結(jié)束后，參照圖5對傳送帶上肥料顆粒的分布規(guī)律進行統(tǒng)計，重復(fù)臺架試驗3次，取其平均值，得到的排肥器的排肥均勻度變異系數(shù)為5.48%。

圖7 臺架驗證試驗

對比仿真正交驗證試驗和臺架驗證試驗的直槽輪排肥器排肥均勻度變異系數(shù)，二者之間存在一定的偏差。通過對偏差出現(xiàn)的原因進行分析，主要是因為在仿真試驗過程，建立的肥料顆粒個體并非標準意義上的球體，顆粒的群體離散程度也并非完全意義上的正態(tài)分布，然而在實際臺架試驗中所使用肥料的均勻度和離散程度相對較差，因此臺架驗證試驗得到的變異系數(shù)略高于仿真試驗的變異系數(shù)，所以這里可以認為仿真試驗結(jié)果和實測結(jié)果基本吻合。通過對仿真正交試驗結(jié)果進行臺架驗證可以得出，借助離散元仿真軟件對直槽輪排肥器的排肥過程進行仿真，以研究其結(jié)構(gòu)參數(shù)對排肥性能影響是可行的。

5 結(jié) 論

1）在直槽輪排肥器工作原理的基礎(chǔ)上，其關(guān)鍵部件直槽輪的結(jié)構(gòu)參數(shù)進行理論分析，建立了單個凹槽端面面積與槽數(shù)、槽輪半徑之間的數(shù)學(xué)關(guān)系，最終將可能影響排肥性能的直槽輪結(jié)構(gòu)參數(shù)轉(zhuǎn)化為凹槽數(shù)目、槽輪有效工作長度、槽輪半徑和凹槽截面形狀。

2）應(yīng)用EDEM軟件搭建直槽輪排肥器工作過程的離散元仿真平臺，進行4個因素之間的交互判定試驗和正交試驗，仿真試驗結(jié)果表明，槽輪半徑與凹槽數(shù)目之間存在一定的交互作用，影響排肥器排肥性能的主次因素依次為凹槽數(shù)目、槽輪半徑、槽輪有效工作長度、槽輪半徑×凹槽數(shù)目和凹槽截面形狀，其中凹槽數(shù)目對排肥均勻度變異系數(shù)影響非常顯著，槽輪半徑對排肥均勻性度系數(shù)影響顯著，最優(yōu)的結(jié)構(gòu)參數(shù)水平組合為槽輪半徑30 mm、凹槽數(shù)目7、有效長度20 mm和凹槽截面形狀為圓弧形，此時直槽輪排肥器的排肥均勻性度變異系數(shù)最小，其值為1.75%。

3）根據(jù)仿真正交試驗得出的最優(yōu)結(jié)構(gòu)參數(shù)水平組合，應(yīng)用3D打印快速成型技術(shù)制造該類型直槽輪排肥器，進行臺架試驗，其排肥均勻度變異系數(shù)為5.48%，對比仿真與臺架試驗結(jié)果，證明借助離散元仿真軟件對直槽輪排肥器的工作過程進行仿真分析，以研究不同結(jié)構(gòu)參數(shù)對排肥性能的影響是可行的。

[1] Zhang Weifeng, Cao Guoxin, Li Xiaolin, et al. Closing yield gaps in China by empowering smallholder farmers[J]. Nature, 2016, 537 (7622): 671.

[2] 羅錫文，廖娟，胡煉，等. 提高農(nóng)業(yè)機械化水平促進農(nóng)業(yè)可持續(xù)發(fā)展[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(1)：1－11.Luo Xiwen, Liao Juan, Hu lian, et al. Improving agricultural mechanization level to promote agricultural sustainable development[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(1): 1－11. (in Chinese with English abstract)

[3] 趙春江，楊信廷，李斌，等. 中國農(nóng)業(yè)信息技術(shù)發(fā)展回顧及展望[J]. 農(nóng)學(xué)學(xué)報，2018，8(1)：172－178. Zhao Chunjiang, Yang Xinting, Li Bin, et al. The retrospect and prospect of agricultural information technology in China[J]. Journal of Agriculture, 2018, 8(1): 172－178. (in Chinese with English abstract)

[4] 陳桂芬，馬麗，陳航. 精準施肥技術(shù)的研究現(xiàn)狀與發(fā)展趨勢[J]. 吉林農(nóng)業(yè)大學(xué)學(xué)報，2013，35(3)：253－259. Chen Guifen, Ma Li, Chen Hang. Research status and development trend of precision fertilization technology[J]. Journal of Jilin Agricultural University, 2013, 35(3): 253－259. (in Chinese with English abstract)

[5] 付衛(wèi)強，孟志軍，黃文倩，等. 基于CAN總線的變量施肥控制系統(tǒng)[J]. 農(nóng)業(yè)工程學(xué)報，2008，24(增刊2)：127－132.Fu Weiqiang, Meng Zhijun, Huang Wenqian, et al. Variable rate fertilizer control system based on CAN bus[J] Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(Supp.2): 127－132. (in Chinese with English abstract)

[6] 左興健，武廣偉，付衛(wèi)強，等. 風(fēng)送式水稻側(cè)深精準施肥裝置的設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(3)：14－21. Zuo Xingjian, Wu Guangwei, Fu Weiqiang, et al. Design and experiment on air-blast rice side deep precision fertilization device[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(3): 14－21. (in Chinese with English abstract)

[7] Tola E, Kataoka T, Burce M, et al. Granular fertiliser application rate control system with integrated output volume measurement[J]. Biosystems Engineering, 2008, 101(4): 411－416.

[8] 孟志軍，趙春江，付衛(wèi)強，等. 變量施肥處方圖識別與位置滯后修正方法[J]. 農(nóng)業(yè)機械學(xué)報，2011，42(7)：204－209. Meng Zhijun, Zhao Chunjiang, Fu Weiqiang, et al. Prescription map identification and position lag calibration method for variable rate application of fertilizer[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(7): 204－209. (in Chinese with English abstract)

[9] 孟志軍，趙春江，劉卉，等. 基于處方圖的變量施肥作業(yè)系統(tǒng)設(shè)計與實現(xiàn)[J]. 江蘇大學(xué)學(xué)報(自然科學(xué)版)，2009，30(4)：338－342. Meng Zhijun, Zhao Chunjiang, Liu Hui, et al. Development and performance assessment of map-based variable rate granule application system[J]. Journal of Jiangsu University (Natural Science Edition), 2009, 30(4): 338－342. (in Chinese with English abstract)

[10] 陳滿，施印炎，汪小旵，等. 基于光譜探測的小麥精準追肥機設(shè)計與試驗[J]. 農(nóng)業(yè)機械學(xué)報，2015，46(5)：26－32. Chen Man, Shi Yinyan, Wang Xiaochan, et al. Design and experiment of variable rate fertilizer applicator based on crop canopy spectral reflectance[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(5): 26－32. (in Chinese with English abstract)

[11] 汪小旵，陳滿，孫國祥，等. 冬小麥變量施肥機控制系統(tǒng)的設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(增刊2)：88－92. Wang Xiaochan, Chen Man, Sun Guoxiang, et al. Design and test of control system on variable fertilizer applicator for winter wheat[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(Supp.2): 88－92. (in Chinese with English abstract)

[12] 苑進，劉成良，古玉雪，等. 基于相關(guān)向量機的雙變量施肥控制序列優(yōu)化[J]. 農(nóng)業(yè)機械學(xué)報，2011，42(增刊1)：184－189，171. Yuan Jin, Liu Chengliang, Gu Yuxue, et al. Bivariate fertilization control sequence optimization based on relevance vector machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(Supp.1): 184－189, 171. (in Chinese with English abstract)

[13] Su Ning, Xu Taosheng, Song Liangtu, et al. Variable rate fertilization system with adjustable active feed-roll length[J]. International Journal of Agricultural & Biological Engineer- ing, 2015, 8(4): 19－26.

[14] 施印炎，陳滿，汪小旵，等. 稻麥精準變量施肥機排肥性能分析與試驗[J]. 農(nóng)業(yè)機械學(xué)報，2017，48(7)：97－103. Shi Yinyan, Chen Man, Wang Xiaochan, et al. Analysis and experiment of fertilizing performance for precision fertilizer applicator in rice and wheat fields[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(7): 97－103. (in Chinese with English abstract)

[15] 施印炎，陳滿，汪小旵，等. 稻麥精準追肥機執(zhí)行機構(gòu)的設(shè)計與試驗[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報，2015，36(6)：119－124. Shi Yinyan, Chen Man, Wang Xiaochan, et al. Design and experiment of precision fertilizer applicator actuator of rice and wheat[J]. Journal of South China Agricultural University, 2015, 36(6): 119－124. (in Chinese with English abstract)

[16] 張濤，劉飛，劉月琴，等. 離散元模擬外槽輪排肥器排量分析[J]. 農(nóng)機化研究，2015，37(9)：198－201. Zhang Tao, Liu Fei, Liu Yueqin, et al. Discrete element simulation of outer groove wheel type fertilizer discharging device capacity analysis[J]. Journal of Agricultural Mechani- zation Research, 2015, 37(9): 198－201. (in Chinese with English abstract)

[17] 楊洲，朱卿創(chuàng)，孫健峰，等. 基于EDEM和3D打印成型的外槽輪排肥器排肥性能研究[J]. 農(nóng)機化研究，2018，40(5)：175－180. Yang Zhou, Zhu Qingchuang, Sun Jianfeng, et al. Study on the performance of fluted roller fertilizer distributor based on EDEM and 3D printing[J]. Journal of Agricultural Mechani- zation Research, 2018, 40(5): 175－180. (in Chinese with English abstract)

[18] 汪博濤，白璐，丁尚鵬，等. 外槽輪排肥器關(guān)鍵工作參數(shù)對排肥量影響的仿真與試驗研究[J]. 中國農(nóng)機化學(xué)報，2017，38(10)：1－6. Wang Botao, Bai Lu, Ding Shangpeng, et al. Simulation and experimental study on impact of fluted-roller fertilizer key parameters on fertilizer amount[J]. Journal of Chinese Agricultural Mechanization, 2017, 38(10): 1-6. (in Chinese with English abstract)

[19] 潘世強，趙亞祥，金亮，等. 2BFJ-6型變量施肥機外槽輪式排肥器的設(shè)計與試驗研究[J]. 中國農(nóng)機化學(xué)報，2016，37(1)：40－42. Pan Shiqiang, Zhao Yaxiang, Jin Liang, et al. Design and experimental research of external grooved wheel fertilizer apparatus of 2BFJ-6 type variable rate fertilizer applicator[J]. Journal of Chinese Agricultural Mechanization, 2016, 37(1): 40－42. (in Chinese with English abstract)

[20] Pamornnak B, Somchai L, Mitchai C. Oil content determi- nation scheme of postharvest oil palm for mobile devices[J]. Biosystems Engineering, 2015, 134: 8－19.

[21] 方會敏，姬長英，張慶怡，等. 基于離散元法的旋耕刀受力分析[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(21)：54－59.Fang Human, Ji Zhangying, Zhang Qingyi, et al. Force analysis of rotary blade based on distinct element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(21): 54－59. (in Chinese with English abstract)

[22] 賀一鳴，吳明亮，向偉，等. 離散元法在農(nóng)業(yè)工程領(lǐng)域的應(yīng)用進展[J]. 中國農(nóng)學(xué)通報，2017，33(20)：133－137. He Yiming, Wu Mingliang, Xiang Wei, et al. Application progress of discrete element method in agricultural engineering[J]. Chinese Agricultural Science Bulletin, 2017, 33(20): 133－137. (in Chinese with English abstract)

[23] 呂昊. 外槽輪排肥器優(yōu)化設(shè)計新方法研究[D]. 長春：吉林大學(xué)，2014.Lü Hao. A new kind of method for optimized design of outer groove-wheel fertilizer apparatuses[D], Changchun: Jilin University, 2014. (in Chinese with English abstract)

[24] 苑進，劉勤華，劉雪美，等. 配比變量施肥中多肥料摻混模擬與摻混腔結(jié)構(gòu)優(yōu)化[J]. 農(nóng)業(yè)機械學(xué)報，2014，45(6)：125－132. Yuan Jin, Liu Qinhua, Liu Xuemei, et al. Simulation of multi-fertilizers blending process and optimization of blending cavity structure in nutrient proportion of variable rate fertilization[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(6): 125－132. (in Chinese with English abstract)

[25] 陳進，周韓，趙湛，等. 基于EDEM的振動種盤中水稻種群運動規(guī)律研究[J]. 農(nóng)業(yè)機械學(xué)報，2011，42(10)：79－83. Chen Jin, Zhou Han, Zhao Zhan, et al. Analysis of rice seeds motion on vibrating plate using EDEM[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(10): 79－83. (in Chinese with English abstract)

[26] 劉彩玲，黎艷妮，宋建農(nóng)，等. 基于EDEM的離心甩盤撒肥器性能分析與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(14)：32－39. Liu Cailing, Li Yanni, Song Jiannong, et al. Performance analysis and experiment on fertilizer spreader with centrifugal swing disk based on EDEM[J]. Transactions of the Chinese Society of Agricultural Engineering (Transac- tions of the CSAE), 2017, 33(14): 32－39. (in Chinese with English abstract)

[27] 陳磊，張繼春，王世博. 新型塑料顆粒3D打印機技術(shù)方案的研究[J]. 機械設(shè)計與制造，2018(1)：83－85. Chen Lei, Zhang Jichun, Wang Shibo. The research of technical solutions of the new plastic particles 3D printer[J]. Machinery Design & Manufacture, 2018(1): 83－85. (in Chinese with English abstract)

[28] 鄭侃，何進，李洪文，等. 基于離散元深松土壤模型的折線破土刃深松鏟研究[J]. 農(nóng)業(yè)機械學(xué)報，2016，47(9)：62－72. Zheng Kan, He Jin, Li Hongwen, et al. Research on polyline soil-breaking blade subsoiler based on sub- soiling soil model using discrete element method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(9): 62－72. (in Chinese with English abstract)

[29] 呂勇，楊堅，梁兆新，等. 單圓盤甘蔗切割器影響破頭率的運動學(xué)仿真[J]. 農(nóng)業(yè)機械學(xué)報，2008，39(4)：50－55. Lü Yong, Yang Jian, Liang Zhaoxin, et al. Simulative kinematics analysis on the affecting factors of rate of broken biennial root of single base cutter of sugarcane harvester[J]. Transactions of the Chinese Society for Agricultural Machinery, 2008, 39(4): 50－55. (in Chinese with English abstract)

[30] 李云雁，胡傳榮. 試驗設(shè)計與數(shù)據(jù)處理[M]. 北京：化學(xué)工業(yè)出版社，2017.

Influences of structure parameters of straight flute wheel on fertilizing performance of fertilizer apparatus

Zhu Qingzhen1,2, Wu Guangwei1※, Chen Liping3, Zhao Chunjiang1,2, Meng Zhijun3

(1.1000972.7121003100097)

Application of chemical fertilizer is an important part in the process of crop production management. The research showed that precise application of chemical fertilizer had great significance to improving the utilization rate of fertilizer and reducing the amount of fertilizer application. The outer fluted roller fertilizer apparatus is the key carrier in the process of precise variable rate fertilization. In order to optimize the fertilizer performance of the existing outer fluted roller fertilizer apparatus comprehensively, the discrete element simulation technology, orthogonal test and 3D rapid prototyping technology were used in this paper. The straight flute wheel structural was taken as research object, and influence regularity of 4 key structural parameters such as groove wheel radius, groove numbers, effective working length and groove section on fertilizer performance was analyzed. A discrete element simulation platform was set up in this paper for outer fluted roller fertilizer apparatus, which mainly consists of fertilization bin, fertilization box, the outer groove roller, fertilization tongue, fertilization pipe, soil ground and fertilizer particles. In the simulation test, the advancing speed of outer fluted roller fertilizer apparatus was 0.5 m/s, the working height was 500 mm, and the rotation speed of outer fluted roller was 20 r/min. The coefficient of variation of the fertilizer was selected as the performance evaluation index of outer fluted roller fertilizer apparatus. The test results of the interaction of 3 factors showed that there was no interaction between groove wheel radius and groove numbers, and between effective working length and groove numbers, there was a certain interaction between groove wheel radius and number of the grooves. The L18(37) standard orthogonal table was selected according to the result of interaction. The orthogonal result showed that factors affecting fertilizer performance, from importance to secondary, mainly include groove numbers, groove wheel radius, effective working length, interaction between groove wheel radius and groove numbers, and groove section. The groove numbers had very significant influence on fertilizer performance. The groove wheel radius had significant influence on fertilizer performance. The effective working length, interaction between groove wheel radius and groove numbers, and groove section had no significant influence on fertilizer performance. The optimum combination of parameters is that the groove numbers is 7, the radius of grooves is 30 mm, the effective working length is 20 mm, and the groove section is circular-arc, and the coefficient of variation of the uniformity of fertilizer for this type of outer fluted roller fertilizer apparatus is 1.75%. In order to test the simulation structure parameters optimization results of outer fluted roller fertilizer apparatus, a bench test was carried out under the optimal combination of structural parameters. The bench test mainly consisted of fluted roller fertilizer apparatus, power motor of conveyor belt, power motor of external groove, control box, conveyor belt and ruler. The experimental results showed that the coefficient of variation of fluted roller fertilizer apparatus is 5.48%. Considering that there may be some errors in the bench test, the results of the bench test were basically consistent with the simulation results, which proved that the discrete element method is feasible to optimize the structural parameters of fluted roller fertilizer apparatus. This research results can provide references for optimizing the structural parameters of fluted roller fertilizer apparatus and improving the performance of the fertilizer.

agricultural machinery;discrete element method; experiments; groove section; structural parameters; fertilizing performance;fertilizer apparatus

10.11975/j.issn.1002-6819.2018.18.002

S224.21

A

1002-6819(2018)-18-0012-09

2018-05-10

2018-06-19

國家重點研發(fā)計劃項目（2016YFD0200600）；國家重點研發(fā)計劃課題（2016YFD0200601）；北京市農(nóng)林科學(xué)院院級科技創(chuàng)新團隊（JNKYT20 1607）

祝清震，男，山東菏澤人，博士生，主要從事小麥基肥精準施用與裝備研究。Email：zhenforyou@163.com。

武廣偉，男，河北趙縣人，博士，副研究員，主要從事精準農(nóng)業(yè)智能裝備研究。Email：wugw@nercita.org.cn。

祝清震，武廣偉，陳立平，趙春江，孟志軍. 槽輪結(jié)構(gòu)參數(shù)對直槽輪式排肥器排肥性能的影響[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(18)：12－20. doi：10.11975/j.issn.1002-6819.2018.18.002 http://www.tcsae.org

Zhu Qingzhen, Wu Guangwei, Chen Liping, Zhao Chunjiang, Meng Zhijun. Influences of structure parameters of straight flute wheel on fertilizing performance of fertilizer apparatus[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 12－20. (in Chinese with English abstract) doi：10.11975/j.issn.1002- 6819.2018.18.002 http://www.tcsae.org