徐照耀，周勇，胡夢(mèng)杰，柯?tīng)Z彬，張夢(mèng)月，呂文

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

我國(guó)棉花主要種植方式分為直播與育苗移栽2種方式[1]，而棉花機(jī)械化直播可省本節(jié)工、效率高，是實(shí)現(xiàn)棉花生產(chǎn)可持續(xù)發(fā)展的重要途徑之一[2].排種器作為機(jī)械化直播的核心工作部件，其排種性能直接決定播種質(zhì)量[3].精量排種器從原理上分為氣力式排種器和機(jī)械式排種器，氣力式排種器對(duì)種子的適應(yīng)性較強(qiáng)，對(duì)種子的損傷較低，但結(jié)構(gòu)較為復(fù)雜，技術(shù)和成本較高.機(jī)械式排種器雖然對(duì)種子尺寸要求嚴(yán)格，存在高速作業(yè)適應(yīng)性差及種子破損率高等問(wèn)題，但有生產(chǎn)成本低、結(jié)構(gòu)簡(jiǎn)單、利于推廣等優(yōu)點(diǎn)[4-6]，其中內(nèi)充種式精量排種器作業(yè)質(zhì)量好、結(jié)構(gòu)簡(jiǎn)單、使用維護(hù)方便，適于作物精播作業(yè)[7-8].

國(guó)內(nèi)對(duì)排種器工作原理的研究主要為設(shè)計(jì)-制造-試驗(yàn)-數(shù)據(jù)分析-改進(jìn)，而后又重復(fù)試驗(yàn)、改進(jìn)、制造成品這一傳統(tǒng)方法，設(shè)計(jì)過(guò)程成本高、周期長(zhǎng),往往影響排種器的開(kāi)發(fā)和推廣.因此，探索新的低成本、短周期的排種器性能研究方法就顯得非常必要[9-10].

離散單元法是專門用于解決不連續(xù)介質(zhì)的一種數(shù)值模擬方法，在排種器性能設(shè)計(jì)及檢驗(yàn)上得到了廣泛的應(yīng)用.劉亞夫基于EDEM軟件建立了雙排型孔輪式油菜排種器的離散元模型，不同排種軸轉(zhuǎn)速下種群擾動(dòng)強(qiáng)度進(jìn)行分析，研究排種軸轉(zhuǎn)速和不同排型孔對(duì)排種器排種量和排種均勻性變異系數(shù)的影響[11].張朋玲建立了離心式排種器的離散元仿真模型，對(duì)油菜籽在排種器中的運(yùn)動(dòng)過(guò)程進(jìn)行了仿真分析，確定了油菜籽的運(yùn)動(dòng)規(guī)律并研究了內(nèi)錐筒中種量與臨界轉(zhuǎn)速的關(guān)系、排種器轉(zhuǎn)速與總排量的關(guān)系，并且仿真結(jié)果與臺(tái)架試驗(yàn)結(jié)果一致，證明采用離散元法分析離心式排種器性能是可行的[12].石林榕基于離散單元法對(duì)水平圓盤式精量排種器進(jìn)行了排種數(shù)值模擬，得到玉米籽粒模型排種、重播、漏播的形成過(guò)程，并通過(guò)臺(tái)架試驗(yàn)驗(yàn)證仿真試驗(yàn).證明基于離散單元法的排種器仿真試驗(yàn)為排種器性能參數(shù)的確定提供參考是可行的[13].

本文以棉花種子為研究生對(duì)象，開(kāi)展內(nèi)充式棉花排種器性能設(shè)計(jì)工作，采用離散元法分析排種器的工作過(guò)程及性能，并將離散元法分析結(jié)果與排種器的臺(tái)架試驗(yàn)結(jié)果進(jìn)行對(duì)比，證明采用離散元法分析內(nèi)充式排種器工作過(guò)程和性能的可行性和有效性，為內(nèi)充式棉花排種器的研究和優(yōu)化設(shè)計(jì)建立一種新方法.

1 內(nèi)充種式棉花排種器結(jié)構(gòu)和工作原理

1.1 結(jié)構(gòu)

如圖1所示內(nèi)充種式棉花排種器主要由排種內(nèi)盤、排種外盤、排種軸、凸耳、毛刷、進(jìn)種管、前殼體、后殼體及入種口高度調(diào)節(jié)板等部分組成.其中排種盤是棉花排種器的核心，其設(shè)計(jì)參數(shù)直接影響棉種的排種質(zhì)量.排種盤由內(nèi)外盤組合而成，內(nèi)外盤貼合在一起形成組合型孔，型孔參數(shù)如圖1-C所示，型孔參數(shù)由矩形型孔高度b(大于最大棉種寬度)、矩形型孔長(zhǎng)度a(大于棉種最大長(zhǎng)度)；梯形型孔下底面長(zhǎng)度g(大于最大棉種寬度小于兩倍棉種最大厚度)、梯形型孔上底面長(zhǎng)度e組成.內(nèi)外盤相對(duì)轉(zhuǎn)動(dòng)可以改變內(nèi)窩孔深度s，控制播量，內(nèi)窩孔寬度為h.

1：毛刷；2：后殼體；3：凸耳；4：排種內(nèi)盤；5：前殼體；6：進(jìn)種管；7：入種口高度調(diào)節(jié)板；8：排種軸；9：排種外盤.1：Brush；2：Back shell；3：Lug；4：Inner seeding plate；5：Front shell；6：Intake pipe；7：Altitude adjusting plate for seed entrance；8：Seeding shaft；9：Outer seeding plate.圖1 排種器結(jié)構(gòu)示意圖Figure 1 Structural diagram of seed-metering device

1.2 工作原理

排種器工作時(shí)，棉種經(jīng)進(jìn)種管，流入充種腔形成種子群；排種盤在排種軸的帶動(dòng)下做周期性逆時(shí)針轉(zhuǎn)動(dòng)；棉種在重力和離心力作用下充入排種盤的型孔；排種盤型孔里多余的種子被毛刷刮掉，并落回充種腔；排種盤型孔中的種子經(jīng)護(hù)種區(qū)，轉(zhuǎn)至投種口，在重力作用下排出.

2 EDEM仿真模型建立

2.1 棉種模型的建立

2.1.1 棉種外部尺寸 棉種的外形尺寸參數(shù)為內(nèi)充種式排種器的尺寸設(shè)計(jì)提供重要參考依據(jù).農(nóng)業(yè)物料的外形尺寸測(cè)定方法較多，其中軸向尺寸法較適用于測(cè)量種子等小物料[14]，故本文采用軸向尺寸法測(cè)定棉種三軸尺寸.如圖2所示定義棉種的長(zhǎng)度(L)、寬度(W)和厚度(T)方向所在的直線分別與三維坐標(biāo)系相互垂直的三軸(x、y、z軸)重合，其中，長(zhǎng)度是指棉種在平面投影中的最大尺寸，寬度是指垂直于長(zhǎng)度方向所在平面棉種邊緣的最大直線距離，厚度是指同時(shí)垂直于長(zhǎng)度與寬度方向棉種邊緣的直線尺寸.為保證離散元仿真的準(zhǔn)確性，對(duì)‘鄂抗棉-10’棉花種子進(jìn)行外形尺寸統(tǒng)計(jì).選取500粒棉花種子用數(shù)顯游標(biāo)卡尺對(duì)其長(zhǎng)度、寬度和厚度直接進(jìn)行測(cè)量，統(tǒng)計(jì)分析所得數(shù)據(jù)，得出‘鄂抗棉-10號(hào)’棉種三軸尺寸平均值.測(cè)試結(jié)果見(jiàn)表1.

圖2 棉花種子三軸尺寸示意圖Figure 2 Triaxial size diagram of cottonseed

表1 棉種三軸尺寸平均值及標(biāo)準(zhǔn)差

2.1.2 棉花種子模型建立 由于棉種顆粒外形尺寸不規(guī)則，其近似為長(zhǎng)橢球體.為使仿真更加真實(shí)地反映出排種器工作過(guò)程，以所測(cè)棉種的平均三軸尺寸為依據(jù)，在EDEM軟件前處理模塊中采用多球面組合構(gòu)型方式創(chuàng)建一個(gè)橢球體組合顆粒模型[15]，棉種顆粒模型由12個(gè)接觸半徑不同的球面粘合而成，并通過(guò)EDEM軟件計(jì)算棉花種子模型的質(zhì)心和質(zhì)量.各球面參數(shù)設(shè)置參數(shù)及棉花種子模型如表2和圖3所示.

表2 棉種顆粒參數(shù)設(shè)置表

圖3 棉種顆粒模型Figure 3 The model of cottonseed particle

2.2 棉花排種器模型建立

將EDEM軟件中顆粒所接觸的幾何體仿真模型可在軟件中直接定義，也可從外界CAD系統(tǒng)中導(dǎo)入，本文采用后者.為提高排種器仿真求解運(yùn)算速度，將排種外盤簡(jiǎn)化為一個(gè)整體，通過(guò)Pro/E軟件繪制出排種器簡(jiǎn)化后的三維模型，并將已裝配完成的排種器三維模型另存為“.stp”格式，在前處理模塊中導(dǎo)入EDEM軟件，其實(shí)體模型如圖4所示.模型透明度設(shè)置為0.3，以便于觀察排種器內(nèi)部棉種運(yùn)動(dòng)規(guī)律.

圖4 排種器模型Figure 4 The model of seed-metering device

3 EDEM仿真試驗(yàn)

3.1 仿真設(shè)置

在材料與物理屬性設(shè)置中，排種盤采用3D打印快速成型技術(shù)構(gòu)造而成，其材料選取為ABS樹(shù)脂，毛刷材料簡(jiǎn)化為橡膠屬性，其余零部件材料均設(shè)置成鋼(steel).參考前人研究成果及前期棉花種子物料特性試驗(yàn)測(cè)定結(jié)果[8]，各材料特性參數(shù)設(shè)置如表3～4所示.

表3 仿真模型材料參數(shù)

Table 3 Material parameter of simulation model

材料Material泊松比 Poisson ratio剪切模量/Pa Shear modulus密度/(kg·m-3)Density棉花Cotton0.141.9×106641.39鋼Steel0.287×1047 850ABS樹(shù)脂ABS Resins0.393.19×1081.05×103

表4 仿真模型材料接觸參數(shù)

排種盤起始旋轉(zhuǎn)時(shí)間1 s，故排種器前1 s為充種階段，排種盤旋轉(zhuǎn)結(jié)束時(shí)間為120 s.棉種總顆粒數(shù)設(shè)定為2 500粒，顆粒生成的速率為1 500粒/s.仿真求解模塊中，時(shí)間步長(zhǎng)設(shè)定為Rayleigh時(shí)間步長(zhǎng)的20%，設(shè)定仿真時(shí)間為120 s，時(shí)間步為1×10-6s.棉種生成平面網(wǎng)格劃分的單元尺寸一般設(shè)置為顆粒最小球面半徑的2倍[16].

3.2 排種器工作過(guò)程仿真

為探究棉種在排種器中的運(yùn)動(dòng)軌跡，選擇棉花種子群中進(jìn)入型孔的籽粒作為觀察對(duì)象.為了較清楚觀察籽粒的運(yùn)動(dòng)過(guò)程，排種器模型以Filled形式顯示opacity設(shè)置為0.1.圖5為內(nèi)充種組合型孔式棉花排種器正常排種仿真過(guò)程.仿真時(shí)間為0 s時(shí)、種子生成平面開(kāi)始生成種子模型，當(dāng)仿真時(shí)間達(dá)到1 s時(shí)排種盤開(kāi)始轉(zhuǎn)動(dòng)，帶動(dòng)棉種群充入排種器組合型孔中，完成充種過(guò)程.當(dāng)仿真時(shí)間為 3.25 s時(shí)內(nèi)外排種盤組合型孔轉(zhuǎn)動(dòng)至種刷區(qū)，多余棉種在毛刷和自身重力的作用下掉落至充種區(qū)，完成清種，而型孔內(nèi)棉種隨排種盤轉(zhuǎn)動(dòng)進(jìn)入護(hù)種區(qū)；當(dāng)仿真時(shí)間為5.58s時(shí)棉種隨型孔到達(dá)排種口，棉種在自身重力和離心力的作用下在排種口處排出進(jìn)入排種管，此時(shí)完成一次充種、清種、護(hù)種、排種的完整過(guò)程.

圖5 排種器排種仿真過(guò)程Figure 5 Simulation process of seed-metering device

由于棉種排種過(guò)程中的隨機(jī)性會(huì)出現(xiàn)重播和漏播兩種情況，通過(guò)EDEM仿真可以發(fā)現(xiàn)問(wèn)題原因.如圖6所示重播是由于種子在型孔內(nèi)排布的隨機(jī)性和排種盤轉(zhuǎn)速過(guò)大的綜合原因下導(dǎo)致4顆棉花種子在型孔中未能及時(shí)清種.如圖7所示漏播主要是由于型孔轉(zhuǎn)入清種區(qū)時(shí)仍存在有卡滯的棉種無(wú)法通過(guò)自重落回充種腔，排種盤繼續(xù)旋轉(zhuǎn)至型孔再次進(jìn)入充種區(qū)后，充入充填孔內(nèi)的棉種因受阻無(wú)法進(jìn)入內(nèi)窩孔，從而導(dǎo)致漏播，且在重力與離心力的共同作用下，卡滯的棉種進(jìn)一步壓實(shí)，如此循環(huán)，此型孔處一直無(wú)棉種排出.

圖6 排種器重播Figure 6 Seed-metering device reseeding

圖7 排種器漏播Figure 7 Seed-metering device miss seeding

通過(guò)仿真分析可以發(fā)現(xiàn)排種器重播、漏播主要由充填孔結(jié)構(gòu)形狀造成的.因此，排種器充填孔進(jìn)行仿真研究得到合理形狀是十分必要的.

3.3 排種器充填孔對(duì)排種器性能影響的單因素試驗(yàn)

3.3.1 試驗(yàn)指標(biāo) 按照農(nóng)藝要求，棉花精量穴播要求每穴(2±1)粒.參照行業(yè)標(biāo)準(zhǔn)NY/T 987-2006《鋪膜穴播機(jī)作業(yè)質(zhì)量》與國(guó)家標(biāo)準(zhǔn)GB/T 6973-2005單粒(精密)播種機(jī)試驗(yàn)方法[17-18]，以排種合格率、漏播率、重播率為試驗(yàn)指標(biāo).其計(jì)算公式如下：

式中，n0指穴粒數(shù)為0粒的總穴數(shù)；n1指穴粒數(shù)為(2±1)粒的總穴數(shù)；n2指穴粒數(shù)為4粒及以上的總穴數(shù)；N為理論排種穴數(shù).

3.3.2 試驗(yàn)因素 為研究充填孔結(jié)構(gòu)對(duì)內(nèi)充種式棉花精量排種器充種性能的影響，以排種器排種合格率、漏播率、重播率為試驗(yàn)指標(biāo)，開(kāi)展充填孔形狀單因素仿真試驗(yàn).充填孔結(jié)構(gòu)分別設(shè)定為凸型、直線型與凹型三個(gè)水平，其所對(duì)應(yīng)的排種盤型孔二維圖如圖8所示.試驗(yàn)中，排種盤內(nèi)窩孔深度設(shè)置為8 mm，入種口高度為60 mm，排種盤轉(zhuǎn)速為25 r/min.

3.3.3 仿真試驗(yàn)結(jié)果 在內(nèi)窩孔深度為8 mm，入種口高度為60 mm，排種盤轉(zhuǎn)速為25 r/min的條件下，凸型、直線型與凹型充填孔結(jié)構(gòu)對(duì)試驗(yàn)指標(biāo)的影響結(jié)果如表5所示.

由表5可知當(dāng)排種器充填孔形狀為凸型時(shí)，排種合格率較小，漏播嚴(yán)重，充填孔形狀采用直線型時(shí)，充填孔容積相應(yīng)增大，排種漏播率減小，合格率增大，當(dāng)充填孔形狀采用凹型時(shí)，排種器達(dá)到較優(yōu)狀態(tài)，此時(shí)，排種合格率為94.6%，漏播率為1.8%，滿足農(nóng)藝要求.

圖8 排種盤型孔二維圖Figure 8 Two-dimensional diagram of seeding plate hole

表5 充填孔形狀對(duì)排種性能的影響

3.4 排種性能仿真試驗(yàn)

3.4.1 試驗(yàn)設(shè)計(jì) 參考研究團(tuán)隊(duì)前期研究成果，取內(nèi)窩孔深度A、入種口高度B和排種盤轉(zhuǎn)速C為主要影響因素，它們的取值范圍為：內(nèi)窩孔深度 6.0～10.0 mm；入種口高度 40.0～80.0 mm；排種盤轉(zhuǎn)速10.0～40.0 r/min，試驗(yàn)采用Design-Expert軟件中的 Box-Behnken方案進(jìn)行三因素三水平響應(yīng)面試驗(yàn).因素水平編碼表如表6所示[19-20].

表6 因素水平編碼表

Table 6 Level coding table mm

3.4.2 試驗(yàn)結(jié)果與分析 根據(jù) Box-Behnken方案共進(jìn)行17組試驗(yàn)，試驗(yàn)結(jié)果見(jiàn)表7，合格率的方差分析見(jiàn)表8，漏播率的方差分析見(jiàn)表9.

表7 響應(yīng)面試驗(yàn)結(jié)果

表8 合格率方差分析結(jié)果

根據(jù)表8和表9的分析結(jié)果可知：對(duì)于合格率，A影響極顯著，C和B2影響顯著，即內(nèi)窩孔深度對(duì)穴粒數(shù)合格率影響極顯著，排種盤轉(zhuǎn)速對(duì)穴粒數(shù)合格率影響顯著，入種口高度對(duì)穴粒合格率影響不顯著.各因素對(duì)棉花排種器的影響顯著性由大到小分別為內(nèi)窩孔深度、排種盤轉(zhuǎn)速、入種口高度；對(duì)于漏播率，C、C2影響極顯著，B2影響顯著，即排種盤轉(zhuǎn)速對(duì)漏播率影響極顯著，各因素對(duì)棉花排種器的影響顯著性由大到小分別為排種盤轉(zhuǎn)速、內(nèi)窩孔深度、入種口高度；各因素交互作用對(duì)排種合格率和漏播率的影響較大.

表9 漏播率的方差分析

對(duì)排種合格率進(jìn)行響應(yīng)面分析，如圖9所示.由圖9-A、圖9-C 可知，排種合格率隨著入種口高度的增大均呈先增大后減小的趨勢(shì)，合格率在入種口高度達(dá)到58 mm時(shí)最大；由圖9-B 和圖9-C 可知，排種合格率隨排種盤轉(zhuǎn)速的增加呈減小趨勢(shì)；由圖9-A 和圖9-B可知，排種合格率隨內(nèi)窩孔增加而減小的趨勢(shì).

圖9 交互因素的響應(yīng)面Figure 9 The response surface of interaction factors

3.4.3 參數(shù)優(yōu)化 當(dāng)排種器內(nèi)窩孔深度為6.0～10.0 mm，入種口高度為40.0～80.0 mm，排種盤轉(zhuǎn)速為10.0～40.0 r/min、采用多重響應(yīng)方法中的主目標(biāo)函數(shù)法對(duì)排種器影響因素進(jìn)行優(yōu)化，分別以合格指數(shù)、重播指數(shù)和漏播指數(shù)作為性能指數(shù)函數(shù)，進(jìn)行優(yōu)化求解，其目標(biāo)函數(shù)和約束條件為：

應(yīng)用軟件Design-Expert中的 Optimization 功能，以合格率最高首要目的，其次為漏播率最低，最后低重播率，進(jìn)行優(yōu)化分析，得最優(yōu)參數(shù)組合為：內(nèi)窩孔深度6.12 mm、投種口高度58.23 mm、轉(zhuǎn)速12.11 r/min.預(yù)測(cè)效果為：合格率98.58%、漏播率0.08%、重播率為1.34%.

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

內(nèi)充種式棉花排種器安裝在JPS-12排種器性能檢測(cè)試驗(yàn)臺(tái)上，采用優(yōu)化參數(shù)組合進(jìn)行臺(tái)架驗(yàn)證試驗(yàn)，測(cè)定臺(tái)架試驗(yàn)中的排種合格率和漏播率，如圖10所示.臺(tái)架試驗(yàn)共進(jìn)行了 10次，對(duì)獲取的排種穴粒數(shù)合格率、漏播率(空穴率)、穴距合格率、穴距變異系數(shù)數(shù)據(jù)各取平均值.臺(tái)架試驗(yàn)結(jié)果如表10所示，內(nèi)充種式棉花排種器的穴粒數(shù)合格率95.31%、漏播率1.66%、穴距合格率93.44%、變異系數(shù)24.12.

圖10 排種器性能試驗(yàn)臺(tái)Figure 10 The structural diagram of the seeder test

表10 驗(yàn)證試驗(yàn)結(jié)果

由于仿真試驗(yàn)中棉花籽粒模型尺寸統(tǒng)一，而臺(tái)架試驗(yàn)中棉花籽粒形狀大小不一，導(dǎo)致臺(tái)架試驗(yàn)的合格率低于仿真試驗(yàn)，重播指數(shù)和漏播指數(shù)較高.而穴粒數(shù)合格率與仿真結(jié)果誤差為3.27%<5%.因此，內(nèi)充種式棉花排種器真試驗(yàn)優(yōu)化排種性能參數(shù)具有可行性.

5 結(jié)論

1) 內(nèi)充種式棉花排種器的仿真試驗(yàn)結(jié)果分析可知：內(nèi)窩孔深度對(duì)重點(diǎn)指標(biāo)穴粒合格率影響極顯著，排種盤轉(zhuǎn)速對(duì)穴粒合格率影響顯著，入種口高度對(duì)穴粒合格率影響不顯著.各因素對(duì)棉花排種器的影響顯著性由大到小分別為內(nèi)窩孔深度、排種盤轉(zhuǎn)速、入種口高度；對(duì)于漏播率，排種盤轉(zhuǎn)速對(duì)漏播率影響極顯著.

2) 以穴粒數(shù)合格率最高首要目的，其次為漏播率最低，最后低重播率，進(jìn)行優(yōu)化分析，得最優(yōu)參數(shù)組合為:內(nèi)窩孔深度6.12mm、投種口高度58.23 mm、轉(zhuǎn)速12.11 r/min.此時(shí)棉花排種器排種合格率為98.58%、漏播率為0.08%、重播率為1.34%.

3) 通過(guò)臺(tái)架試驗(yàn)，在較優(yōu)組合條件下測(cè)得排種合格率為95.31%、漏播率為1.66%、重播率D為3.03%，重點(diǎn)指標(biāo)合格率誤差3.27%(<5%).因此，基于EDEM的排種器離散元仿真試驗(yàn)為排種器性能參數(shù)的確定是可行的.