氣吸圓盤式微型薯排種器充種性能模擬與試驗

賴慶輝 馬文鵬 劉 素 蘇 微 張智泓

(昆明理工大學(xué)現(xiàn)代農(nóng)業(yè)工程學(xué)院， 昆明 650500)

氣吸圓盤式微型薯排種器充種性能模擬與試驗

賴慶輝 馬文鵬 劉 素 蘇 微 張智泓

(昆明理工大學(xué)現(xiàn)代農(nóng)業(yè)工程學(xué)院， 昆明 650500)

為提高氣吸圓盤式微型薯排種器充種性能，以云南麗薯6號微型薯為播種對象，基于離散元法，以種子平均法向應(yīng)力方差為指標，對振動頻率和振動幅度分別進行數(shù)值模擬，并對上述因素進行單因素試驗，試驗結(jié)果與仿真效果一致。結(jié)果表明：增加振動頻率和振動幅度可以增大種子平均法向力方差，增強對種子的擾動性，從而提高充種性能。為尋求最佳工作參數(shù)組合，采用三因素五水平二次旋轉(zhuǎn)正交組合試驗方法，對排種器進行排種性能試驗，并對試驗結(jié)果進行優(yōu)化與驗證。結(jié)果表明：在作業(yè)速度為2.4 km/h，吸種負壓為6 kPa，種層高度為70 mm，振動頻率為6.5～6.9 Hz，振動幅度為20～21 mm時，合格指數(shù)大于95，重播指數(shù)和漏播指數(shù)小于2.5。

微型薯； 氣吸圓盤式排種器； 充種性能； 數(shù)值模擬； 試驗

引言

目前國內(nèi)外先進的排種器普遍采用了氣力式原理[1-6]，這是因為氣力式排種器具有對種子外形要求不嚴，不傷種，能適應(yīng)較高速度播種作業(yè)等優(yōu)點。對于氣力式排種器，充種環(huán)節(jié)尤為重要，其充種性能直接影響播種質(zhì)量。采用種箱振動的方式提高充種性能是最常見的方法[7-9]，即通過振動激勵使振動種盤內(nèi)種群產(chǎn)生“沸騰”運動，以便氣力吸種部件完成吸種過程。周海波[10]采用振動原理，研制了一種電磁振動種室和氣動振盤相結(jié)合的振動式精密播種裝置,減免傷種并實現(xiàn)均勻穩(wěn)定播種；王朝暉等[11]、楊明金等[12]利用振動理論、高速攝影與振動分析儀，采用二次正交旋轉(zhuǎn)組合試驗方法，得出吸種性能與種層厚度、振動頻率等因素間的回歸模型；龔智強[13]分析了變質(zhì)量系統(tǒng)種層厚度與振動種盤振動參數(shù)的關(guān)系,建立了排種裝置工作參數(shù)與吸種距離之間的數(shù)學(xué)模型；鹿芳媛等[14]基于離散元法，模擬振動盤工作過程，探究振動參數(shù)對種子運動規(guī)律的影響，尋求最佳振動頻率；賴慶輝等[15]采用DEM-CFD耦合方法，研究振動參數(shù)對種間平均法向應(yīng)力和供種高度的影響，得到最優(yōu)工作參數(shù)組合；陳進等[16]研究了氣吸振動播種試驗臺內(nèi)種子的運動機理，通過對試驗臺內(nèi)種子運動過程的計算機動態(tài)模擬，得出了振動臺內(nèi)種子的運動規(guī)律。

微型薯屬于大粒種子，為解決其不易充種的問題，設(shè)計了一種振動和氣吸組合式精密排種器，并通過簡易的正交試驗得到了一組關(guān)于作業(yè)速度、吸種負壓和振動頻率的最優(yōu)工作參數(shù)組合，排種合格指數(shù)可達到90以上[17]。但在振動供種過程中，振動幅度和種層高度也會影響種群運動，進而影響充種性能。因此，需深入研究各項振動參數(shù)對充種性能的影響機理，并尋求最優(yōu)振動參數(shù)組合，以進一步提高排種合格指數(shù)。由于種群與振動種盤的碰撞以及種子之間相互碰撞關(guān)系非常復(fù)雜，目前從理論上還無法建立排種裝置吸種過程種群運動學(xué)模型，因此，本文以自行研制的氣吸圓盤式排種器為研究對象，借助離散元仿真軟件EDEM進行振動頻率和振動幅度對種群擾動影響的仿真分析，得出振動過程中種群運動參數(shù)的變化規(guī)律，并進行臺架驗證試驗。采用三因素五水平二次旋轉(zhuǎn)正交組合試驗方法，選取振動頻率、振動幅度、種層高度為試驗因素進行排種性能試驗，尋求各參數(shù)的最佳區(qū)間組合，提高充種性能的同時增強排種器的適用性。

1 排種器工作原理與振動理論

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

振動供種氣吸圓盤式微型薯排種器結(jié)構(gòu)如圖1所示，主要由振動供種機構(gòu)、氣力排種機構(gòu)和機架等組成。

圖1 振動供種氣吸圓盤式排種器結(jié)構(gòu)示意圖Fig.1 Structure diagram of vibrating feeding pneumatic disc seed-metering device1.種箱 2.機架 3.密封槽 4.氣室 5.微型薯 6.吸孔 7.排種盤 8.振動板

氣室外圍設(shè)有密封槽，排種盤與密封條緊密貼合，機架上安裝有種箱和振動供種機構(gòu)。微型薯從種箱進入振動板，在振動作用下種群處于“沸騰”狀態(tài)，接通風機使氣室內(nèi)部形成負壓區(qū)，播種機地輪帶動排種盤轉(zhuǎn)動，微型薯在吸孔內(nèi)外壓差的作用下被吸附在吸孔上，并隨著排種盤一起運動。當微型薯運動到氣室末端時，隨著吸孔移出氣室，負壓消失，微型薯在自身重力作用下從吸孔脫落掉入種溝，實現(xiàn)排種過程。

圖2 振動供種機構(gòu)結(jié)構(gòu)示意圖Fig.2 Structure diagram of vibration mechanism1.振動板 2.連桿 3.曲柄盤 4.調(diào)節(jié)孔 5.支架 6.直線軸承 7.直線軸

振動激勵裝置按照工作原理可分為電液式、電動式和機械式。其中偏心式機械振動機構(gòu)與其他振動激勵方式相比，結(jié)構(gòu)簡單，振動幅度和振動頻率相互獨立可分別進行調(diào)節(jié)，有利于研究振動頻率、振動幅度對種群運動的影響關(guān)系。圖2為振動供種機構(gòu)結(jié)構(gòu)示意圖，振動板在連桿的帶動和直線軸的限位作用下，實現(xiàn)上下振動。其中，曲柄盤的轉(zhuǎn)速由步進電動機控制，振動板振動頻率可調(diào)節(jié)范圍為0～10 Hz；曲柄盤上設(shè)有調(diào)節(jié)孔，可通過改變曲柄長度來調(diào)節(jié)振動幅度，調(diào)節(jié)范圍為10～50 mm。

1.2 振動板運動分析

將振動供種機構(gòu)簡化為曲柄連桿機構(gòu)，機構(gòu)簡圖如圖3所示，取振動板的上止點為振動板運動坐標原點O′，則由圖中幾何關(guān)系可得[18]

x=R+L-Rcosα-Lcosβ

(1)

式中x——O′點位移R——曲柄長度L——連桿長度α——曲柄與豎直方向夾角β——連桿與豎直方向夾角

記λ=R/L，為曲柄半徑與連桿長度之比，當λ?1時，式(1)可近似簡化為

(2)

其中

ω=2πf

(3)

(4)

式中A——振動幅度ω——振動圓頻率f——振動頻率t——時間

于是可求得振動板加速度表達式為

(5)

若使種群處于“沸騰”狀態(tài)，需滿足

max{2Aπ2f2(cos(2πft)+λcos(2πft))}>g

(6)

式中g(shù)——重力加速度

由式(6)可知，振動幅度和振動頻率將對種子在振動盤上的運動產(chǎn)生較大影響，且參數(shù)選擇范圍較廣，為最佳工作參數(shù)的選取造成較大干擾。因此本文采用離散元法，利用EDEM(Engineering discrete element method)軟件探究振動供種過程中種子的運動規(guī)律及參數(shù)對其運動規(guī)律的影響，為確定最佳振動參數(shù)，提高振動播種裝置的播種性能提供參考。

圖3 曲柄滑塊振動機構(gòu)Fig.3 Slider-crank mechanism

1.3 種群運動分析

定義排種器振動部件振動強度K為振動板豎直方向上加速度的最大值與重力加速度的比值[12]，設(shè)種群總層數(shù)為p，且最下層種子為第p層，最上層種子為第1層。

在種子被拋擲的過程中，種群內(nèi)部運動可以看成種子之間的相互碰撞，則位于φ(1≤φ≤p)層的種子經(jīng)過p-φ次碰撞后的振幅Aφ與振動盤振幅A關(guān)系滿足[19]

(7)

式中r——顆粒間恢復(fù)系數(shù)，r<1

則第1層種子拋擲強度需滿足

(8)

因此排種器振動板振動強度需滿足

(9)

由于r<1，由式(9)可知，種層厚度對振動盤內(nèi)種子的運動具有較大影響。種層厚度越大，即p值越大，所需振動板的振動強度越大。因此，需尋求振動頻率、振動幅度以及種層高度之間的最佳參數(shù)組合，以提高充種性能。

2 數(shù)值模擬

2.1 接觸模型

離散元法描述的是顆粒之間碰撞的過程，也就是接觸的產(chǎn)生和發(fā)生作用的過程。離散元中根據(jù)接觸方式的不同可分為硬顆粒接觸和軟顆粒接觸，軟顆粒接觸方式允許顆粒的接觸點間出現(xiàn)重疊部分，并根據(jù)接觸顆粒的物理屬性和法向重疊量、切向位移計算出接觸力，而顆粒間內(nèi)摩擦力即切向力與法向力有著密切的聯(lián)系，因此本文采用軟顆粒接觸模型。常用的軟顆粒接觸模型包括：摩擦電荷接觸模型、線彈性接觸模型、運動表面接觸模型、線性黏附接觸模型、Hertz-Mindlin黏結(jié)接觸模型以及Hertz-Mindlin無滑動接觸模型。其中Hertz-Mindlin無滑動接觸模型為EDEM軟件中默認接觸模型，又考慮到微型薯顆粒之間粘附力較小，故本文選取Hertz-Mindlin無滑移接觸模型[20-24]。軟顆粒接觸使用的是軟球模型，如圖4所示，軟球模型在顆粒i和顆粒j間設(shè)定了彈簧、阻尼器、滑動器和耦合器，引入了剛度k和阻尼系數(shù)c來量化彈簧、阻尼器、滑動器的作用。

圖4 軟球模型對顆粒間接觸力的簡化模型Fig.4 Simplified model of particles contact force from soft-sphere model

耦合器是用來確定發(fā)生接觸的顆粒配對關(guān)系，不引入任何力。在切向上，如果切向力超過屈服值，兩顆粒在法向力和摩擦力作用下滑動，由滑動器實現(xiàn)這一目的[25-27]。

2.2 接觸力計算

法向力Fnij是彈簧和法向阻尼器作用在顆粒i上的彈性力和阻尼力的合力。對于三維球形顆粒，根據(jù)Hertz接觸理論，F(xiàn)nij表示為

(10)

式中kn——法向剛度cn——法向阻尼系數(shù)a——法向重疊量vij——顆粒i相對于顆粒j的速度，m/sn——從顆粒i球心到顆粒j球心的單位矢量

切向力Ftij表示為

Ftij=-ktδ-ctvct

(11)

式中kt——切向剛度ct——切向阻尼系數(shù)vct——接觸點的滑移速度，m/sδ——接觸點的切向位移

如果

|Ftij|>μs|Fnij|

(12)

則顆粒i發(fā)生滑移，切向力為

(13)

式中μs——靜摩擦因數(shù)

式(13)為庫侖摩擦定律，從式中可知，種子間切向力即為內(nèi)摩擦力，而內(nèi)摩擦力是阻礙種子間相對運動或者有相對運動趨勢的主要因素[28]。從公式中可以看出，影響內(nèi)摩擦力的因素有靜摩擦因數(shù)和法向力，而靜摩擦因數(shù)又由種子本身物理屬性決定，因此本文采用法向力Fnij來衡量種子之間的內(nèi)摩擦力[29-30]。

2.3 振動供種過程仿真

在Pro/E軟件中建立振動板幾何模型并導(dǎo)入EDEM，種子顆粒通過多球面組合填充的方式在EDEM中生成(圖5a)并添加經(jīng)測量所得的力學(xué)特性(表1)，種子由振動板模型內(nèi)部的顆粒工廠生成，生成的種子在自身重力作用下分布在振動板底部，如圖5b所示。對振動板添加振頻和振幅，在此作用下，種群出現(xiàn)“沸騰”狀態(tài)，種間出現(xiàn)相對運動與碰撞，EDEM軟件獲取實時顆粒體的位置與接觸信息，離散元仿真過程輸出每個時間步長的接觸平均法向力、接觸總數(shù)、種子總數(shù)。由于從種群生成完畢到種群運動達到穩(wěn)定狀態(tài)需要1 s左右的時間，因此每組仿真數(shù)據(jù)從2.0 s時開始輸出，6.0 s時結(jié)束。i時刻每個種子受到的平均法向力為

(14)

式中Fni——i時刻接觸法向力的平均值，NNi——i時刻種子總數(shù)Nci——i時刻接觸總數(shù)

圖5 振動供種數(shù)值模擬結(jié)果Fig.5 Simulation results of seed-filling

表1 模擬所需參數(shù)Tab.1 Material parameters used in simulation

2.3.1 振動頻率

圖6 不同頻率下平均法向力變化曲線Fig.6 Changing curves of average normal stress under different frequencies

為了能夠反映振動頻率對充種性能的影響，在前期試驗中選取能夠凸顯這一影響的3個水平進行模擬仿真，根據(jù)振動試驗實際情況，選取振幅為20 mm時，模擬振動頻率分別為4、6、8 Hz，獲取種子隨時間變化的平均法向力，如圖6所示。從圖6中可以看出種子平均法向力隨時間出現(xiàn)無規(guī)律波動，并隨著振頻的提高，波動幅度增大且落差也隨著增大。為了定量描述這種波動，引入方差公式

(15)

其中

式中D——方差N——輸出數(shù)據(jù)總數(shù)，N=100z——種子平均法向力總和，N

方差較大的其波動離散程度大，也就是對種子的擾動性較大。通過計算得出振動頻率為4 Hz時，方差為102.54 N2；振動頻率為6 Hz時，方差為113.48 N2；振動頻率為8 Hz時，方差為130.62 N2。較大的振動頻率可以使得種群波動性增強，從而提高種子充填的概率。因此提高振動頻率可以在一定

程度上提高充種性能，為正交優(yōu)化試驗因素及因素水平范圍的選取提供參考。

2.3.2 振動幅度

為了能夠反映振動幅度對充種性能的影響，選取振動頻率為6 Hz，在振動幅度分別為10、20、30 mm時，獲取種間隨時間變化的平均法向力，如圖7所示。從圖7中可以看出隨著振幅的提高，種子平均法向力也隨著增大。通過計算得出振動幅度為10 mm時，平均法向力方差為86.35 N2；振動幅度為20 mm時，平均法向力方差為113.48 N2；振動幅度為30 mm時，平均法向力方差為140.78 N2。因此適當?shù)卦黾诱駝臃瓤梢蕴岣叱浞N性能，為正交優(yōu)化試驗因素及因素水平范圍的選取提供參考。

圖7 不同振幅下平均法向力變化曲線Fig.7 Changing curves of average normal stress under different amplitudes

3 排種性能試驗

3.1 試驗條件

研制振動供種氣吸圓盤式排種器樣機，選取“麗薯6號”微型薯作為排種對象，在JPS-12型計算機視覺排種器性能檢測試驗臺上進行試驗，如圖8所示。該試驗臺以種床帶為種床，排種器在試驗時保持固定狀態(tài)，種床相對于排種器進行運動，微型薯排落在涂有油層的種床帶上，圖像采集處理系統(tǒng)對種床帶上的微型薯進行實時檢測，檢測完畢后，系統(tǒng)根據(jù)GB/T 6973—2005《單粒(精密)播種機試驗方法》對檢測數(shù)據(jù)進行處理保存，輸出粒距合格指數(shù)、漏播指數(shù)和重播指數(shù)。利用AZ8000型手持非接觸式紅外線轉(zhuǎn)速表測定曲柄轉(zhuǎn)速，轉(zhuǎn)換為振動頻率，并通過控制步進電動機轉(zhuǎn)速對頻率進行調(diào)節(jié)；在曲柄盤徑向上開有數(shù)個調(diào)節(jié)孔，連桿一端通過調(diào)節(jié)孔與曲柄鉸接，不同的調(diào)節(jié)孔對應(yīng)不同曲柄長度，進而調(diào)節(jié)振動板的振動幅度(2倍曲柄長度)；在種箱與振動板接合處，安裝一塊擋板，可通過調(diào)整擋板位置來控制種層高度。

圖8 排種性能試驗臺Fig.8 Seed-metering performance test1.數(shù)據(jù)處理單元 2.控制柜 3.排種試驗臺 4.氣吸圓盤式排種器 5.微型薯

3.2 單因素試驗

為驗證數(shù)值模擬試驗結(jié)果的準確性，進行單因素試驗。分別以振動頻率和振動幅度作為試驗因素，以合格指數(shù)、重播指數(shù)及漏播指數(shù)作為評價指標，連續(xù)記錄吸附200粒種子的排種試驗結(jié)果，每次試驗重復(fù)5次取平均值。

3.2.1 振動頻率對排種性能的影響

設(shè)定振動幅度為20 mm，吸種負壓為6 kPa，作業(yè)速度為2.4 km/h，分別在振動頻率為3.5、4.5、5.5、6.5、7.5、8.5、9.5 Hz的條件下進行排種試驗，得到振動頻率對試驗指標的影響，結(jié)果如圖9所示。

圖9 不同振動頻率時排種性能試驗曲線Fig.9 Seeding performance test curves under different vibration frequencies

由圖9可知，隨著振動頻率的增加，排種合格指數(shù)呈先上升后下降的趨勢，重播指數(shù)逐漸下降，漏播指數(shù)逐漸下降后上升。振動頻率在5.5～7.5 Hz的范圍內(nèi)，排種合格指數(shù)大于80，在振動頻率為6.5 Hz時，排種合格指數(shù)達到最大值94.1。振動頻率過小時，振動對于種群的作用較小，種子之間的內(nèi)摩擦力較大，不利于種群“沸騰”，隨著振動頻率的提高，振動板對種群的擾動性增強，漏播指數(shù)逐漸下降，有效地提高了充種效果；但隨著振動頻率進一步提高，種子出現(xiàn)了較強的波動，使得原本被吸附的種子被波動較強的種子撞擊掉落，導(dǎo)致漏播指數(shù)上升。

3.2.2 振動幅度對排種性能的影響

設(shè)定振動頻率為6.5 Hz，吸種負壓為6 kPa，作業(yè)速度為2.4 km/h，分別在振動幅度為5、10、15、20、25、30、35 mm的條件下進行排種試驗，得到振動幅度對試驗指標的影響，結(jié)果如圖10所示。

圖10 不同振動幅度時排種性能試驗曲線Fig.10 Seeding performance test curves under different amplitudes

由圖10可知，隨著振動幅度的增加，排種合格指數(shù)呈現(xiàn)先上升后下降的趨勢，重播指數(shù)逐漸上升，漏播指數(shù)逐漸下降。在振動幅度大于15 mm時，排種合格指數(shù)均大于80，在振動幅度為20 mm時，排種合格指數(shù)達到最大值95.4。振動幅度過小時，種群會隨著振動板一起上下運動，影響充種效果；振動幅度過大，會縮短種子處于“沸騰”狀態(tài)的時間，致使種群不能規(guī)則地隨振動板一起振動。

3.3 二次回歸旋轉(zhuǎn)正交試驗

3.3.1 試驗方案

為尋求最佳參數(shù)組合，進行三因素二次旋轉(zhuǎn)正交組合試驗。選取振動頻率、振動幅度和種層高度作為試驗因素，合格指數(shù)、重播指數(shù)和漏播指數(shù)作為試驗指標。設(shè)定排種器作業(yè)速度為2.4 km/h，吸種負壓為6 kPa，試驗因素和編碼見表2，試驗方案與結(jié)果見表3(X1、X2、X3為因素編碼值)，每組試驗重復(fù)5次取平均值。

表2 試驗因素和編碼Tab.2 Factors and levels of test

表3 試驗設(shè)計方案與結(jié)果Tab.3 Experimental design and results

3.3.2 回歸數(shù)學(xué)模型的建立與顯著性檢驗

采用Design-Expert 8.0.6軟件對試驗數(shù)據(jù)進行多元回歸擬合，對試驗結(jié)果進行回歸分析，可以得到合格指數(shù)Y1、重播指數(shù)Y2和漏播指數(shù)Y3的回歸方程。

(1)合格指數(shù)Y1回歸模型建立與顯著性檢驗

通過試驗以及對試驗數(shù)據(jù)進行多元回歸擬合，得到各因素對排種合格指數(shù)Y1影響的回歸模型為

(16)

回歸方程的顯著性檢驗如表4所示。根據(jù)表4可知，該模型的擬合度極顯著(P<0.01)。但振動幅度、供種高度、振動頻率和振動幅度的交互項(X1X2)的P值、振動頻率和種層高度的交互項(X1X3)的P值以及振動幅度和種層高度的交互項(X2X3)的P值均大于0.05，說明以上各項對排種合格指數(shù)的影響不顯著，其他各項的F檢驗均極顯著或顯著，說明相關(guān)試驗因素對響應(yīng)值的影響存在二次關(guān)系。對于失擬項P=0.2094，不顯著，說明不存在其他影響指標的主要因素。剔除不顯著因素后的回歸模型為

(17)

通過對式(17)回歸系數(shù)的檢驗得出，影響排種合格指數(shù)的因素主次順序為振動頻率、振動幅度和種層高度。

(2)重播指數(shù)回歸模型的建立與顯著性檢驗

通過試驗以及對試驗數(shù)據(jù)進行多元回歸擬合，得到各因素對重播指數(shù)Y2影響的回歸模型為

(18)

回歸方程的顯著性檢驗如表4所示。根據(jù)表4可知，該模型的擬合度極顯著(P<0.01)。但振動幅度、振動頻率和振動幅度交互項(X1X2)的P值、振動頻率和種層高度交互項(X1X3)的P值以及振動幅度和種層高度交互項(X2X3)的P值均大于0.05，說明以上各項對重播指數(shù)的影響不顯著，其他各項的F檢驗均極顯著或顯著，說明相關(guān)試驗因素對響應(yīng)值的影響存在二次關(guān)系。對于失擬項P=0.097 5，不顯著，說明不存在其他影響指標的主要因素存在。剔除不顯著因素后的回歸模型為

(19)

通過對式(19)回歸系數(shù)的檢驗得出，影響重播指數(shù)的因素主次順序為種層高度、振動頻率和振動幅度。

表4 回歸方程方差分析Tab.4 Variance analysis of regression equation

注：*顯著(P<0.05)； ** 極顯著(P<0.01)。

(3)漏播指數(shù)Y3回歸模型建立與顯著性檢驗

通過試驗以及對試驗數(shù)據(jù)進行多元回歸擬合，得到各因素對漏播指數(shù)Y3影響的回歸模型為

(20)

回歸方程的顯著性檢驗如表4所示。根據(jù)表4可知，該模型的擬合度極顯著(P<0.01)。但振動頻率和振動幅度交互項(X1X2)的P值、振動頻率和種層高度交互項(X1X3)的P值以及振動幅度和種層高度交互項(X2X3)的P值均大于0.05，說明各因素交互項對漏播指數(shù)的影響不顯著，其他各項的F檢驗均極顯著或顯著，說明相關(guān)試驗因素對響應(yīng)值的影響存在二次關(guān)系。對于失擬項P=0.078 4，不顯著，說明不存在其他影響指標的主要因素存在。剔除不顯著因素后的回歸模型為

(21)

通過對式(21)回歸系數(shù)的檢驗得出，影響漏播指數(shù)的因素主次順序為振動幅度、振動頻率和種層高度。

3.4 各因素對排種合格指數(shù)的影響

通過Design-Expert 8.0.6對數(shù)據(jù)進行處理，可得到振動頻率、振動幅度和供種高度對合格指數(shù)的影響，其響應(yīng)曲面如圖11所示。任意固定某個因素的水平，根據(jù)響應(yīng)曲面圖，分析其余2個因素間的交互作用對排種合格指數(shù)的影響。

圖11 交互因素對合格指數(shù)的影響Fig.11 Effects of interactive factors on eligible rate

3.4.1 振動頻率和振動幅度的交互作用

圖11a為種層高度為70 mm時，振動頻率和振動幅度對排種合格指數(shù)交互作用的響應(yīng)曲面圖。由圖11a可知，在振動頻率為6.5～7.0 Hz，振動幅度為20～21 mm時，排種合格指數(shù)較高。振動頻率一定時，隨著振動幅度的增大，排種合格指數(shù)先上升后下降。振動幅度一定時，隨著振動頻率的增大，排種合格指數(shù)同樣先上升后下降。振動頻率為6.8 Hz，振動幅度為20 mm時，排種合格指數(shù)最高。

3.4.2 振動頻率和種層高度的交互作用

圖11b是振動幅度為20 mm時，振動頻率和種層高度對排種合格指數(shù)交互作用的響應(yīng)曲面圖。由圖11b可知，在振動頻率為6.5～7.0 Hz，種層高度為67～73 mm時，排種合格指數(shù)較高。振動頻率一定時，隨著種層高度的增大，排種合格指數(shù)先上升后下降。種層高度一定時，隨著振動頻率的增大，排種合格指數(shù)同樣先上升后下降。振動頻率為6.8 Hz，種層高度為70 mm時，排種合格指數(shù)最高。

3.4.3 振動幅度和種層高度的交互作用

圖11c為振動頻率為6.8 Hz時，振動幅度和種層高度對排種合格指數(shù)交互作用的響應(yīng)曲面圖。由圖11c可知，在振動幅度為20～21 mm，種層高度為67～73 mm時，排種合格指數(shù)較高。振動幅度一定時，隨著種層高度的增大，排種合格指數(shù)先上升后下降。種層高度一定時，隨著振動幅度的增大，排種合格指數(shù)同樣先上升后下降。振動幅度為20 mm，種層高度為70 mm時，排種合格指數(shù)最高。

3.5 最佳參數(shù)優(yōu)化

設(shè)定合格指數(shù)大于95，重播指數(shù)小于2.5，漏播指數(shù)小于2.5，優(yōu)化得最佳參數(shù)范圍如圖12所示(x3=70 mm)。

圖12 參數(shù)優(yōu)化分析圖Fig.12 Diagram of parameters optimization and analysis

由圖12可知，在種層高度為70 mm時，黃色區(qū)域為參數(shù)優(yōu)化區(qū)域，即振動頻率為6.5～6.9 Hz，振動幅度為20～21 mm時，可獲得合格指數(shù)大于95，漏播指數(shù)小于2.5，重播指數(shù)小于2.5。

對優(yōu)化后的理論結(jié)果進行試驗驗證。在相同的試驗條件下選取振動頻率為6.5～6.9 Hz，振動幅度為20～21 mm，種層高度為70 mm，進行3次重復(fù)驗證試驗，得到排種器合格指數(shù)平均值為96.6，且均大于95；漏播指數(shù)平均值為2.1，且均小于2.5；重播指數(shù)平均值為2.3，且均小于2.5；試驗結(jié)果與優(yōu)化結(jié)果基本相符。

將本文最終的優(yōu)化結(jié)果與之前正交試驗的優(yōu)化結(jié)果[17]進行對比，排種合格指數(shù)有所提高，且得到了各工作參數(shù)的最優(yōu)區(qū)間，增強了排種器的實際操作性。

4 結(jié)論

(1)基于離散元法，利用EDEM軟件，分別對振動頻率和振動幅度進行數(shù)值模擬，并以種子平均法向力方差為指標，對結(jié)果進行分析，得出增大振動頻率和振動幅度可以有效提高排種器充種性能的結(jié)論。

(2)以合格指數(shù)、重播指數(shù)和漏播指數(shù)為試驗指標，分別對振動頻率和振動幅度進行單因素試驗，試驗現(xiàn)象及效果與仿真分析結(jié)果一致。驗證了所建仿真模型的正確性以及采用離散元法模擬仿真振動供種的可行性。

(3)為尋找最佳工作參數(shù)組合，以振動頻率、振動幅度和種層高度為試驗因素，以合格指數(shù)、重播指數(shù)和漏播指數(shù)作為試驗指標，進行三因素五水平二次旋轉(zhuǎn)正交組合試驗，并對結(jié)果進行分析和優(yōu)化：影響排種合格指數(shù)的因素主次順序為振動頻率、振動幅度和種層高度。在吸種負壓為6 kPa，作業(yè)速度為2.4 km/h，種層高度為70 mm，振動頻率為6.5～6.9 Hz，振動幅度為20～21 mm時，可獲得合格指數(shù)大于95，重播指數(shù)小于2.5，漏播指數(shù)小于2.5。對優(yōu)化結(jié)果進行驗證試驗，得到的驗證結(jié)果與優(yōu)化結(jié)果基本一致。

1 趙佳樂，賈洪雷，姜鑫銘，等. 大豆播種機偏置雙圓盤氣吸式排種器[J/OL].農(nóng)業(yè)機械學(xué)報，2013，44(8)：78-83. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20130814&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2013.08.014. ZHAO Jiale, JIA Honglei, JIANG Xinming, et al.Suction type offset double disc seed metering device of soybean seeder[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(8): 78-83. (in Chinese)

2 王延耀，李建東，王東偉，等.氣吸式精密排種器正交試驗優(yōu)化[J/OL] .農(nóng)業(yè)機械學(xué)報，2012，43(增刊)：54-58，89. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=2012s11&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2012.S0.011. WANG Yanyao, LI Jiandong, WANG Dongwei, et al. Orthogonal experiment optimization on air-suction precision seed-metering device[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(Supp.): 54-58, 89. (in Chinese)

3 左彥軍，馬旭，玉大略，等.水稻芽種窩眼窄縫式氣吸滾筒排種器流場模擬與試驗[J].農(nóng)業(yè)機械學(xué)報，2011，42(2)：58-62. ZUO Yanjun, MA Xu, YU Dalue, et al. Flow field numerical simulation of suction cylinder-seeder for rice bud seed with socket slot[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(2)：58-62. (in Chinese)

4 劉佳，崔濤，張東興，等.機械氣力組合式玉米精密排種器[J/OL] .農(nóng)業(yè)機械學(xué)報，2012，43(2)：43-47. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20120209&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2012.02.009. LIU Jia, CUI Tao, ZHANG Dongxing, et al.Mechanical-pneumatic combined corn precision seed-metering device[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(2): 43-47. (in Chinese)

5 趙曉順，于華麗，張晉國，等.槽縫氣吸式小麥精量排種器[J/OL] .農(nóng)業(yè)機械學(xué)報，2013，44(2)：48-51，74. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20130210&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2013.02.010. ZHAO Xiaoshun, YU Huali, ZHANG Jinguo, et al. Slot-type pneumatic precise wheat seed-metering device[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(2): 48-51, 74. (in Chinese)

6 GAIKWAD B B, SIROHI N P S. Design of a low-cost pneumatic seeder for nursery plug trays[J]. Biosystems Engineering, 2008, 99(3): 322-329.

7 楊麗，史嵩，崔濤，等.氣吸與機械輔助附種結(jié)合式玉米精量排種器[J/OL] .農(nóng)業(yè)機械學(xué)報，2012，43(增刊)：48-53. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=2012s10&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2012.S0.010. YANG Li, SHI Song, CUI Tao, et al. Air-suction corn precision metering device with mechanical supporting plate to assist carrying seed[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(Supp.): 48-53. (in Chinese)

8 趙湛，李耀明，陳進，等. 氣吸滾筒式排種器吸種過程的動力學(xué)分析[J]. 農(nóng)業(yè)工程學(xué)報，2011，27(7)：112-116. ZHAO Zhan, LI Yaoming, CHEN Jin, et al. Dynamic analysis of seeds pick-up process for vacuum-cylinder seeder[J].Transactions of the CSAE, 2011, 27(7): 112-116. (in Chinese)

9 張石平，夏靜，陳進. 氣吸振動式蔬菜穴盤育苗精密播種裝置的研究[J]. 農(nóng)機化研究，2007，29(8)：80-83. ZHANG Shiping, XIA Jing, CHEN Jin. Study on equipment of vibration air-suction vegetable tray precise seeder[J].Journal of Agricultural Mechanization Research, 2007, 29(8): 80-83. (in Chinese)

10 周海波. 水稻秧盤育秧精密播種機的關(guān)鍵技術(shù)研究與應(yīng)用[D]. 長春：吉林大學(xué)，2009. ZHOU Haibo. Research and application on key technologies of precision seeder for tray nursing seedings of rice[D]. Changchun: Jilin University, 2009. (in Chinese)

11 王朝暉，馬旭，賈瑞昌. 工作參數(shù)對超級稻育秧播種部件吸種性能的影響[J]. 農(nóng)業(yè)工程學(xué)報，2009，25(8)：88-92. WANG Zhaohui, MA Xu, JIA Ruichang. The working parameters of suction effect on the performance of super rice seedling planting parts[J]. Transactions of the CSAE, 2009, 25(8): 88-92. (in Chinese)

12 楊明金，邱兵，楊玲，等. 振動氣吸式精密穴播機的排種性能[J]. 農(nóng)業(yè)工程學(xué)報，2010，26(9)：139-143. YANG Mingjin, QIU Bing, YANG Ling, et al. Metering performance of the vibration vacuum precision tray seeder[J].Transactions of the CSAE, 2010, 26(9): 139-143. (in Chinese)

13 龔智強.氣吸振動盤式精密排種裝置理論與試驗研究[D]. 鎮(zhèn)江：江蘇大學(xué)，2013. GONG Zhiqiang. Theoretical and experimental study on vacuum-vibration tray precision seeding device[D]. Zhenjiang:Jiangsu University, 2013. (in Chinese)

14 鹿芳媛，馬旭，齊龍，等. 基于離散元法的雜交水稻振動勻種裝置參數(shù)優(yōu)化與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(10)：17-24. LU Fangyuan, MA Xu, QI Long, et al. Parameter optimization and experiment of vibration seed-uniforming device for hybrid rice based on discrete element method[J].Transactions of the CSAE, 2016,32(10):17-24. (in Chinese)

15 賴慶輝，高筱鈞，張智泓. 三七氣吸滾筒式排種器充種性能模擬與試驗[J/OL]. 農(nóng)業(yè)機械學(xué)報，2016，47(5)：27-37. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20160505&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2016.05.005. LAI Qinghui, GAO Xiaojun, ZHANG Zhihong. Simulation and experiment of seed-filling pneumatic cylinder seed-metering device forPanaxnotoginseng[J/OL].Transactions of the Chinese Society for Agricultural Machinery, 2016,47(5):27-37. (in Chinese)

16 陳進，李耀明.氣吸振動式播種試驗臺內(nèi)種子運動規(guī)律的研究[J] .農(nóng)業(yè)機械學(xué)報，2002，33(1):47-50. CHEN Jin, LI Yaoming. Study on seeds movement law in sowing test stand with suction and vibration[J]. Transactions of the Chinese Society for Agricultural Machinery, 2002, 33(1): 47-50. (in Chinese)

17 賴慶輝，馬文鵬，蘇微.氣吸圓盤式微型薯排種器設(shè)計與試驗[J/OL].農(nóng)業(yè)機械學(xué)報，2016，47(12)：30-37. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20161205&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2016.12.005. LAI Qinghui, MAWenpeng, SU Wei.Design and experiment of pneumatic disc seed-metering device for mini-tuber[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(12): 30-37. (in Chinese)

18 李耀明，邱白晶.氣吸振動式水稻播種試驗臺的振動分析[J].農(nóng)業(yè)機械學(xué)報，1998，29(3)：43-47. LI Yaoming, QIU Baijing. Vibration analysis of rice seeding test bench[J]. Transactions of the Chinese Society for Agricultural Machinery, 1998, 29(3):43-47.(in Chinese)

19 GUARELLA P, PELLERANO A, PASCUZZI S. Experimental and theoretical performance of a vacuum seeder nozzle for vegetable seeds[J]. Journal of Agricultural Engineering Research, 1996,64(1): 29-36.

20 VAN LIEDEKERKE P, TIJSKEN E, DINTWA E, et al. DEM simulation of the particle flow on a centrifugal fertilizer spreader[J]. Powder Technology, 2009, 190(3):348-360.

21 YOUSHIYUKI S, CUNDALL P A. Three dimensional DEM simulation of bulk handling by screw conveyors[J]. Journal of Engineering Mechanics, 2001, 127(9):864-872.

22 胡建平，郭坤，周春健，等. 磁吸滾筒式排種器種箱振動供種仿真與試驗[J/OL]. 農(nóng)業(yè)機械學(xué)報，2014，45(8)：61-65. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20140810&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2014.08.010. HU Jianping, GUO Kun, ZHOU Chunjian, et al. Experiment of supplying seeds in box of magnetic precision cylinder-seeder[J/OL]. Transactions of the Chinese society for Agricultural Machinery, 2014,45(8):61-65. (in Chinese)

23 李帥，李曄，馬履翱，等. 基于EDEM對圓形導(dǎo)料溜槽的仿真分析[J]. 煤礦機械，2014,33(6)：71-73. LI Shuai, LI Ye, MA Lüao, et al. Simulation analysis of cylindrical chute based on EDEM[J]. Coal Mine Machinery, 2014, 33(6):71-73.(in Chinese)

24 陳進，周韓，趙湛，等. 基于EDEM的振動種盤中水稻種群運動規(guī)律研究[J]. 農(nóng)業(yè)機械學(xué)報，2011,42(10)：79-83，100. 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,100. (in Chinese)

25 于建群，申燕芳，牛序堂，等. 組合內(nèi)窩孔精密排種器清種過程的離散元仿真分析[J]. 農(nóng)業(yè)工程學(xué)報，2008,24(5)：105-109. YU Jianqun, SHEN Yanfang, NIU Xutang, et al. DEM simulation and analysis of the clearing process in precision metering device with combination inner-cell[J]. Transactions of the CSAE, 2008, 24(5): 105-109. (in Chinese)

26 王國強，郝萬軍，王繼新. 離散單元法及其在EDEM上的實踐[M]. 西安：西北工業(yè)大學(xué)出版社，2010.

27 胡國明. 顆粒系統(tǒng)的離散元素法分析仿真[M]. 武漢：武漢理工大學(xué)出版社，2010.

28 吳愛祥，孫業(yè)志，劉湘平. 散體動力學(xué)理論及其應(yīng)用[M]. 北京：冶金工業(yè)出版社，2002.

29 騰云楠，曾慶輝，姚紅良，等. 振動摩擦機理及其非線性動力學(xué)特性[J]. 摩擦學(xué)學(xué)報，2010,30(2)：197-202. TENG Yunnan, ZENG Qinghui, YAO Hongliang, et al. Mechanism of vibration friction and its non-linear dynamic characteristics[J]. Tribology, 2010, 30(2):197-202. (in Chinese)

30 孔維姝，胡林，杜學(xué)能，等. 用探測棒研究顆粒堆中最大靜摩擦力[J]. 物理學(xué)報，2007,56(4)：2318-2322. KONG Weishu, HU Lin, DU Xueneng, et al. A simulation static friction force in granular pile with probing road[J]. Acta Physica Sinica, 2007, 56(4):2318-2322.(in Chinese)

Simulation and Experiment on Seed-filling Performance of Pneumatic Disc Seed-metering Device for Mini-tuber

LAI Qinghui MA Wenpeng LIU Su SU Wei ZHANG Zhihong
(CollegeofModernAgriculturalEngineering,KunmingUniversityofScienceandTechnology,Kunming650500,China)

In order to realize the mechanized planting of mini-tuber and solve the problem of large seed filling, a vibrating feeding pneumatic disc seed-metering device for mini-tuber was designed. The seed-metering device was composed of metering plate, air chamber, frame and vibrating seed feeding device. Seeding-filling performance was the key factor that influenced the properties of metering device. In order to improve the seed-filling performance of the pneumatic disc seed-metering device, Lishu No.6 mini-tuber was adopted as sowing object, the method of UDEC was used as research technique, meanwhile, variance yields of average normal stress was considered as the optimizing targets. Through single factor experiment, the vibration frequency and vibration amplitude were numerically simulated respectively. It can be found that the experimental phenomenon and effect were in consistent with the simulation ones. Results showed that vibration frequency and vibration amplitude increased variance yields of average normal stress of seeds, hence disturbance of seeds was enhanced. Based on the three-factor five-level quadratic orthogonal rotating combination test method, the influence of vibration frequency, vibration amplitude and layer height on sowing performance were explored. The regression equations of the three factors and five levels were fitted through the processing by Design-Expert 8.0.6 (experimental design expert) software. Based on the results of examination, it can be found that the fitting of the equation was good, and the best parameter combination was the operating speed of 2.4 km/h, negative pressure of 6 kPa, vibration frequency of 6.5～6.9 Hz, vibration amplitude of 20～21 mm and layer height of 70 mm. Under the optimal condition, the qualified index was greater than 95, and the multiple index as well as the missing index was less than 2.5. It showed that the vibrating feeding pneumatic disc seed-metering device met the standard and requirements by comparing the results with the national standard. The research approach was suitable for the exploitation of mini-tuber seed-metering, and it provided a theoretical reference for the design of pneumatic disc seed-metering device.

mini-tuber; pneumatic seed-metering device; seed-filling performance; numerical simulation; experiment

2017-02-21

2017-03-21

國家自然科學(xué)基金項目(51305187)、云南省重點新產(chǎn)品開發(fā)計劃項目(2014BC007)、云南省科技計劃青年項目(2015FD011)、昆明理工大學(xué)自然科學(xué)基金項目(KKSY201323067、KKSY201323025、14118940)和云南省教育廳項目(2016ZZX048)

賴慶輝(1980—)，男，副教授，主要從事農(nóng)業(yè)機械裝備與計算機測控研究，E-mail: laiqinghui007@163.com

10.6041/j.issn.1000-1298.2017.05.005

S223.2+5

A

1000-1298(2017)05-0044-10