基于EDEM的傾斜圓盤勺式大豆排種器投種性能優(yōu)化研究

許健 甘義權(quán) 蔡宗壽

摘要 為提高傾斜圓盤勺式大豆排種器投種穩(wěn)定性和均勻性，建立投種過程的運動學(xué)模型，分析得出排種盤轉(zhuǎn)速和投種角是影響投種性能的主要因素。利用離散元軟件EDEM進行仿真研究，設(shè)計二因素二次回歸正交旋轉(zhuǎn)組合試驗，運用SPSS軟件對試驗數(shù)據(jù)進行處理，以合格指數(shù)和變異系數(shù)作為排種器性能評價指標(biāo)，分別建立其與排種盤轉(zhuǎn)速、投種角的回歸方程，利用Matlab繪制三維等值線圖，確定了影響合格指數(shù)和變異系數(shù)的主要因素。對試驗因素進行優(yōu)化計算，得出最優(yōu)參數(shù)組合為排種盤轉(zhuǎn)速41.7 r/min，投種角19.2°，此時合格指數(shù)為94.88%，變異系數(shù)為12.17%。

關(guān)鍵詞 大豆排種器；投種機理；離散元法；旋轉(zhuǎn)正交試驗

中圖分類號 S223.2 文獻標(biāo)識碼 A 文章編號 0517-6611（2018）21-0193-04

Abstract In order to improve the seeding stability and uniformity of the inclined disc scoop soybean seeding device， the kinematics model of the seeding process was established. It was analyzed that the rotational speed and the seeding angle of the seeding disc were the main factors affecting the seeding performance. Discrete element software EDEM was used for simulation study. A twofactor quadratic regression orthogonal rotation combination experiment was designed. SPSS software was used to process the experimental data. Seed index and coefficient of variation were used as evaluation indexes of seedmeters to establish the regression equation of the rotation speed and the seeding angle of the seed disk， the threedimensional contour map was drawn using Matlab， and the main factors affecting each performance index were determined. The optimization factors of the experimental factors were calculated. The optimal parameter combination was 41.7 r/min seeding disc rotation and 19.2° seeding angle. The passing index was 94.88% and the coefficient of variation was 12.17%.

Key words Soybean seed metering device；Seeding mechanism；Discrete element method；Orthogonal rotational test

精密播種是利用精密播種機將種子合理均勻地播入預(yù)定位置的先進技術(shù)，是現(xiàn)代播種作業(yè)的發(fā)展趨勢。排種器是實現(xiàn)精密播種技術(shù)的核心部件，其工作性能和穩(wěn)定性將直接影響播種機的工作質(zhì)量[1-2]。目前國內(nèi)外學(xué)者對傾斜圓盤勺式大豆排種器的研究主要集中在充種、取種相關(guān)部件的理論研究和改進設(shè)計，對投種環(huán)節(jié)的研究甚少[3]。因此，研究排種器投種過程運動學(xué)模型，分析影響排種器投種性能的因素，對提高排種器排種性能具有深遠意義。

離散元法具有計算大量顆粒在給定條件下如何運動的能力，被廣泛應(yīng)用于排種器性能分析[4-5]。胡建平等[6]運用顆粒離散元法研究影響磁吸板式排種器充種性能的主要因素。張濤等[7]運用離散元法分析振動對排種器內(nèi)玉米種群運動的影響規(guī)律。劉月琴等[8]用離散元法得出氣吸式排種器最佳工作參數(shù)并驗證了基于離散元法仿真優(yōu)化的可靠性。筆者以傾斜圓盤勺式大豆排種器為研究載體，對其投種機理進行研究分析并建立投種過程運動學(xué)模型，在此基礎(chǔ)上，運用離散元軟件EDEM作二次正交旋轉(zhuǎn)組合試驗，分析排種盤轉(zhuǎn)速與投種角對排種均勻性和穩(wěn)定性的影響，旨在利用一定的計算方法獲得排種器的最佳工作參數(shù)組合。

1 排種器投種機理分析

1.1 排種器結(jié)構(gòu)和工作原理

傾斜圓盤勺式大豆排種器結(jié)構(gòu)如圖1a所示，主要由取種勺、排種盤、隔板、導(dǎo)種輪、殼體等部件組成。其工作原理主要分為充種、清種、遞種、護種和投種5個串聯(lián)階段，如圖1b所示。排種器正常工作時，種子在排種盤的旋轉(zhuǎn)攪動下進行分種，種子在自身重力、種子間的相互作用力等力共同作用下，側(cè)向充填進取種勺內(nèi)，完成充種過程；排種盤繼續(xù)向上轉(zhuǎn)動，穩(wěn)定存在于取種勺內(nèi)的多粒種子，在力系發(fā)生轉(zhuǎn)變的情況下，逐漸回落至充種區(qū)，直至取種勺內(nèi)僅存一粒種子，完成清種過程；單粒種子由隔板開口處進入導(dǎo)種輪，完成遞種環(huán)節(jié)；導(dǎo)種輪與排種盤同步轉(zhuǎn)動直至來到投種口投放，完成護種、投種過程。

1.2 投種過程運動學(xué)分析

分析大豆顆粒在投種過程的三維運動軌跡是研究投種機理的重要參考，當(dāng)大豆顆粒在投種口投出時，大豆顆粒將作空間拋物線運動。對投送過程中的大豆種子進行運動學(xué)分析，以排種盤轉(zhuǎn)動中心為坐標(biāo)原點，建立空間直角坐標(biāo)系，如圖2所示。大豆顆粒的空間軌跡位移如下：

式中，ω為排種盤角速度（rad/s）；R為排種盤半徑（mm）；α為投種角（°）；β為排種器垂直傾角（°）；t為種子投送運動時間（s）。

通過對投種過程運動學(xué)分析可知，在排種器整體結(jié)構(gòu)參數(shù)確定之后，排種器投種性能的均勻性和穩(wěn)定性主要與排種盤轉(zhuǎn)速、排種器投種角有關(guān)。因此，在后續(xù)基于離散元的仿真試驗階段，將采用二因素二次回歸正交旋轉(zhuǎn)組合的設(shè)計方法分析研究[9]。

2 EDEM虛擬仿真

2.1 參數(shù)設(shè)置

仿真模型通過Pro/E創(chuàng)建并通過igs.格式導(dǎo)入到EDEM中，由于仿真模擬只需將與種子接觸的幾何部件導(dǎo)入，所以將排種器模型進行簡化，如圖3所示。

仿真研究將大豆創(chuàng)建為四面體構(gòu)型，并以中黃39大豆為研究對象，隨機取樣1 200粒，分別對其長、寬、厚的三維尺寸進行統(tǒng)計，按其尺寸頻率分布均值建立顆粒模型。統(tǒng)計得出，中黃39大豆的球形率高達97.3%，因此在EDEM中設(shè)置顆粒為球體。由于大豆種子表面無黏附作用，因此選擇Hertz-MindLin（no slip）built-in為仿真接觸模型。仿真參數(shù)設(shè)置如表1所示[10]。

為了保證仿真時排種器工作有足夠的顆粒和仿真時間，設(shè)置生成1 200顆大豆種子模型，半徑平均值為3.452 mm，采用正態(tài)分布的方式生成顆粒尺寸，標(biāo)準(zhǔn)差為0.086，設(shè)置固定時間步長為Rayleigth時間步長的16%，仿真時間總長20 s，其中前2 s為充種過程，仿真過程如圖4所示。

2.2 仿真方案

以二因素二次回歸正交旋轉(zhuǎn)組合的設(shè)計方法安排仿真試驗，選取合格指數(shù)和變異系數(shù)為試驗評價指標(biāo)：

式中，n1為單粒排種數(shù)；N為理論排種數(shù)；σ為排種標(biāo)準(zhǔn)差；A為合格指數(shù)；V為變異系數(shù)。

根據(jù)前述理論分析及實際排種器作業(yè)要求，合理的控制試驗因素變化范圍，每個試驗重復(fù)3次，取其平均值作為試驗結(jié)果，因素水平編碼如表2所示。

2.3 結(jié)果與分析

2.3.1 回歸分析。

仿真試驗方案與結(jié)果如表3所示。試驗結(jié)果采用SPSS數(shù)據(jù)分析軟件進行回歸分析，以確定各指標(biāo)的擬合優(yōu)度及其回歸方程。

合格指數(shù)回歸方程檢驗表明，F(xiàn)=13 403.39，對應(yīng)的P<0.000 1，因此回歸方程高度顯著，其回歸方程模型如下：

A=58.256+1.429X1+0.631X2-0.016X12-0.011X22-0.004X1X2（4）

同理可得變異系數(shù)回歸方程模型如下：

V=30.721-0.536X1-0.801X2+0.007X12+0.019X22+0.001X1X2（5）

2.3.2 排種性能分析。

為了更加直觀地得到各個試驗因素與性能評價指標(biāo)的關(guān)系，利用Matlab畫出三維等值線圖，如圖5～6所示。

通過各性能評價指標(biāo)的回歸方程和三維等值線圖可知，影響其性能的排種盤轉(zhuǎn)速和取種角存在交互作用。分析圖5表明，當(dāng)投種角位于零水平時，合格指數(shù)隨著排種盤轉(zhuǎn)速的增加呈現(xiàn)先上升后下降的趨勢；當(dāng)排種盤轉(zhuǎn)速位于零水平時，合格指數(shù)隨著取種角的增加呈現(xiàn)先上升后下降的趨勢；由于當(dāng)排種盤轉(zhuǎn)速變化時，相應(yīng)的合格指數(shù)變化區(qū)域更大，因此主要影響合格指數(shù)的因素是排種盤轉(zhuǎn)速。分析圖6表明，當(dāng)投種角位于零水平時，變異系數(shù)隨著排種盤轉(zhuǎn)速的增加呈現(xiàn)上升的趨勢；當(dāng)排種盤轉(zhuǎn)速位于零水平時，變異系數(shù)隨著投種角的增加呈現(xiàn)先下降后上升趨勢；由于當(dāng)取種角變化時，相應(yīng)的變異系數(shù)變化區(qū)域更大，因此主要影響變異系數(shù)的因素是取種角。

2.3.3 排種器性能優(yōu)化。

利用Matlab中非線性優(yōu)化fmincon函數(shù)，優(yōu)化處理結(jié)果如下：排種盤轉(zhuǎn)速41.70 r/min，取種角19.2°。在此條件下進行仿真驗證試驗，試驗表明排種器合格指數(shù)為94.88%，變異系數(shù)為12.17%。

3 結(jié)論

（1）以傾斜圓盤勺式大豆排種器為研究載體，對其投種機理進行研究分析并建立投種過程運動學(xué)模型，分析得出了排種盤轉(zhuǎn)速和投種角是影響投種性能的主要因素。

（2）運用EDEM作正交排種虛擬試驗，采用二次回歸正交旋轉(zhuǎn)組合試驗設(shè)計方法建立了合格指數(shù)和變異系數(shù)的回歸方程。運用Matlab繪制2個考察指標(biāo)的等值線圖，得到了試驗因素對各指標(biāo)的影響趨勢并確定了影響合格指數(shù)和變異系數(shù)的主要因素。

（3）確定了該型排種器的最佳工作參數(shù)組合，當(dāng)排種盤轉(zhuǎn)速為41.70 r/min時，投種角為19.2°，排種穩(wěn)定性和均勻性最優(yōu)，此時排種器合格指數(shù)為94.88%，變異系數(shù)為12.17%。

參考文獻

[1] 王延耀，李建東，王東偉，等.氣吸式精密排種器正交試驗優(yōu)化[J].農(nóng)業(yè)機械學(xué)報，2012，43（S1）：54-58，89.

[2] 王傳鵬，何瑞銀.基于單片機的精密排種器性能檢測裝置的研究[J].科學(xué)技術(shù)與工程，2011，11（33）：8299-8302.

[3]王金武，唐漢，王金峰，等.指夾式玉米精量排種器導(dǎo)種投送運移機理分析與試驗[J].農(nóng)業(yè)機械學(xué)報，2017，48（1）：29-37，46.

[4] 史嵩，張東興，楊麗，等.基于EDEM軟件的氣壓組合孔式排種器充種性能模擬與驗證[J].農(nóng)業(yè)工程學(xué)報，2015，31（3）：62-69.

[5] 廖慶喜，張朋玲，廖宜濤，等.基于EDEM的離心式排種器排種性能數(shù)值模擬[J].農(nóng)業(yè)機械學(xué)報，2014，45（2）：109-114.

[6] 胡建平，周春健，侯沖，等.磁吸板式排種器充種性能離散元仿真[J].農(nóng)業(yè)機械學(xué)報，2014，45（2）：94-98.

[7] 張濤，劉飛，趙滿全，等.基于離散元的排種器排種室內(nèi)玉米種群運動規(guī)律[J].農(nóng)業(yè)工程學(xué)報，2016，32（22）：27-35.

[8] 劉月琴，趙滿全，劉飛，等.基于離散元的氣吸式排種器工作參數(shù)仿真優(yōu)化[J].農(nóng)業(yè)機械學(xué)報，2016，47（7）：65-72.

[9] GAO G H，F(xiàn)ENG T X，YANG H，et al.Development and optimization of endeffector for extraction of potted anthurium seedlings during transplanting[J].Applied engineering in agriculture，2016，32（1）：37-46.

[10] TAVAREZ F A，PLESHA M E.Discrete element method for modelling solid and particulate materials[J].International journal for numerical methods in engineering，2007，70（4）：379-404.

[11] 耿端陽，張明源，何珂，等.傾斜雙圓環(huán)型孔圓盤式玉米排種器設(shè)計與試驗[J].農(nóng)業(yè)機械學(xué)報，2018，49（1）：68-76.