輪勺式大蒜單粒取種裝置設計與試驗

李玉華 張智龍 李天華 吳彥強 牛子孺 侯加林

(1.山東農業大學機械與電子工程學院，泰安 271018；2.山東省農業裝備智能化工程實驗室，泰安 271018)

0 引言

大蒜單粒取種[1]是實現大蒜機械化播種的關鍵技術之一。大蒜顆粒大、形狀不規則及表面粗糙，極易產生顆粒堆積[2-3]、結拱堵塞等現象[4]，從而導致漏充率高。因此，相對比較成熟的玉米取種器[5-6]、小麥取種器[7]和大豆取種器[8]等，無法通過簡單改進而直接應用于大蒜單粒取種。

現有大蒜單粒取種方式主要有勺鏈式、勺帶式、夾持式及振動式等[9-10]。大蒜屬于中大籽粒種子，其形狀尺寸變異系數大，大蒜取種器漏充及重充現象難以避免，勺鏈及勺帶式取種器結構簡單，但中高速時振動明顯，易引起掉種。近年來有學者對大蒜取種進行了研究，耿愛軍等[11]采用鏈勺式結構進行大蒜單粒取種，對取種勺及鏈板進行了優化；此外有學者對大籽粒種子取種器[12-15]進行了研究，為大蒜單粒取種裝置設計提供了技術參考。

針對大蒜單粒取種的漏播及重播問題，闡述輪勺式大蒜單粒取種裝置的工作原理，利用數理統計[16-17]方法對取種勺的結構尺寸區間進行分析。通過離散元仿真軟件EDEM[18-20]對取種過程及取種勺結構參數進行仿真分析，利用Box-Benhnken中心組合設計理論設計試驗方案[21-23]，并進行臺架驗證試驗，以期為大蒜單粒取種裝置的設計提供參考。

1 取種器結構與工作原理

1.1 取種器結構

輪勺式大蒜取種器主要由種箱、驅動電機、取種輪、取種勺組成，結構如圖1所示。

圖1 輪勺式大蒜取種器結構示意圖

取種勺與取種輪通過圓弧面配合，可方便拆卸及更換。驅動電機通過傳動軸驅動取種輪運動，取種輪帶動取種勺旋轉，取種勺與大蒜通過攪動擠壓使其進入取種勺空腔實現取種。

1.2 單粒取種原理

大蒜單粒取種過程分為4個階段，如圖2所示，分別為：取種區、輸種區、排種區和過渡區。

圖2 輪勺式大蒜取種器工作示意圖

圖3 取種過程大蒜受力分析

取種過程分析：取種區中進入取種勺的待取大蒜受到自身重力、周邊大蒜及取種勺的作用力，如圖3a，受力公式為

(1)

式中m——大蒜質量，kg

Fi(i=1,2,…,n)——周圍大蒜對待取大蒜的作用力，N

N——取種勺對大蒜的支持力，N

a——待取大蒜的加速度，m/s2

Mo——大蒜所受力對支撐點的力矩，N·m

按取種勺內有大蒜和無大蒜兩種情況對大蒜運動狀態進行分析：取種勺內無大蒜時，待取大蒜將在取種勺周圍大蒜的擠壓及勺壁的支撐力作用下進入取種勺的空腔內；取種勺內已有大蒜時，勺內大蒜受到周圍大蒜及取種勺的作用力，此時取種勺周圍大蒜處于競爭狀態，勺內大蒜將趨于更加穩定，否則將被周邊大蒜替代，因此經過取種區后勺內大蒜處于穩定狀態，不易在輸種區滑落而引起漏種。

輸種過程分析：輸種區中取種勺的大蒜將在自身重力和取種勺的支撐力作用下做勻速圓周運動，如圖3b，受到重力及取種勺的支撐力，徑向方向受力公式為

(2)

式中ω——取種輪角速度，rad/s

r——取種輪半徑，m

θ——取種勺、取種輪中心連線與水平方向夾角，rad

Fn——大蒜勻速圓周運動所需向心力，N

當重力產生的分力小于向心力時取種勺中的大蒜將沿切向方向滑出而產生漏種，取種勺提供的向心力與ω、r、θ及取種勺輪廓密切相關，需滿足ω2r≤gsinθ。

排種過程分析：排種區的大蒜將在自身重力和取種勺的支撐力作用下做勻速圓周運動，如圖3c，隨θ的減小重力分量無法提供足夠的向心力，使大蒜沿切向方向飛出完成排種。

2 關鍵結構設計與仿真

大蒜具有籽粒大、形狀極不規則和個體差異大等特點，因此易產生漏播和重播現象，本文以概率統計的大蒜外形尺寸為依據對大蒜單粒取種器進行了設計，分析了取種勺的形狀、尺寸及取種輪轉速對單粒取種性能的影響。

2.1 大蒜三維尺寸

大蒜外形為非規則的類扇形體，其長度Ld、寬度Wd、厚度Td對取種勺尺寸影響最大，本文以金鄉大蒜為研究對象，隨機選取100粒飽滿的大蒜進行測量，結果如表1所示。

表1 大蒜三維尺寸

對所測數據的概率分布進行分析，概率密度曲線如圖4所示，各尺寸基本符合正態分布，其中Ld～N(32.38,2.222),Wd～N(16.32,3.162),Td～N(18.98,2.912)，取種勺尺寸參數設計以此大蒜三維尺寸分布為依據。

圖4 大蒜三維尺寸概率分布

2.2 離散元仿真模型

為縮短試驗時間，提高試驗針對性，本文采用離散元分析軟件EDEM對取種器取種性能進行仿真分析。以飽滿的金鄉大蒜為仿真分析對象，大蒜顆粒的尺寸根據體積按正態分布進行設置。顆粒間及顆粒與設備間的接觸模型均采用Hertz-Mindlin(no slip)模型，仿真參數如表2所示[24]。

表2 仿真參數

2.3 取種勺設計與仿真分析

取種勺是大蒜單粒取種器的核心部件，其結構和尺寸均對漏充率和重充率具有直接的影響。

2.3.1取種勺結構的確定

通過離散元分析軟件對多種不同結構的取種勺進行仿真分析，對比分析發現圖5所示的4種取種勺結構的取種性能相對較好，其中D型取種勺漏充率和重充率在相同條件下比其他3種結構低，本文選用D型取種勺進行結構參數優化。

A型取種勺漏充率及重充率均相對較高，因該取種勺弧高恒定，取種勺容積大，容易一次取兩粒大蒜，且取種勺兩端開口在輸種階段大蒜會從兩端掉落導致較高的漏充率；B型取種勺弧高恒定且兩端封閉，雖然漏充率較低但重充率較高；C型取種勺的中間大兩端小，重充率低，但兩端開口導致漏充率較高；D型取種勺弧高中間大兩端小，且兩端的邊緣凸起能避免大蒜在輸種過程中滑落，克服了B、C兩種結構的缺點而保留了其優點，因此D型取種勺具有相對較低的重充率和漏充率。

2.3.2取種勺尺寸的確定

取種勺的關鍵結構參數為半徑R、長度L和弧高H。為保證單粒取種成功率，取種勺尺寸需要在合適的范圍內。

為保證在取種勺長度方向不能同時容納1粒以上蒜種且空間盡量大以提高取種率，取種勺長度L需滿足

(3)

由多個相互獨立隨機變量函數分布的性質知：2Ld～N(64.76，4.442),Ld+Wd～N(48.7，3.862),Ld+Td～N(51.36，3.662)。式(3)各約束條件對應的分布函數曲線如圖6所示，隨L取值增大重充率升高而漏充率降低，為保證單粒取種合格率，計算分析得L取值范圍為35.00～45.00 mm時能兩者兼顧，當L為45.00 mm時重充率為14.50%。

圖6 取種勺長度L取值區間

為保證在取種勺深度方向取種后不滑落且不能同時容納1粒以上蒜種，取種勺的弧高H需滿足

(4)

由多個相互獨立隨機變量函數分布的性質知：1.5Wd～N(24.48，4.742),1.5Td～N(28.47，4.372),0.5Wd～N(8.16，1.582),0.5Td～N(9.49，1.452)。式(4)各約束條件對應的分布函數曲線如圖7所示，隨H取值增大重充率升高而漏充率降低，為保證單粒取種合格率，計算分析得H取值范圍為13.00～17.00 mm。仿真和單因素試驗發現弧高H在此區間范圍內對單粒取種性能影響不大。

圖7 取種勺弧高H取值區間

為保證在取種勺寬度方向單粒取種成功率，取種勺的弦長S需滿足

(5)

由多個相互獨立隨機變量函數分布的性質知：2Wd～N(32.64，6.322),2Td～N(37.96，5.822),Wd～N(16.32，3.162),Td～N(18.98，2.912)。式(5)各約束條件對應的分布函數曲線如圖8所示，隨S取值增大重充率升高而漏充率降低，為保證單粒取種合格率，同時考慮弧高的限制，并通過先期試驗驗證，S取值范圍為26.5～34.0 mm時取種效果較好。

圖8 取種勺弦長S取值區間

由三角幾何關系知取種勺半徑R滿足

(6)

把弧高H及弦長S的取值范圍代入式(6),計算得取種勺半徑R取值范圍為13.30～17.00 mm。

2.3.3取種輪設計與仿真分析

取種輪外圓弧面與取種勺配合形成取種空腔，取種輪外形對取種合格率有顯著影響。通過仿真及試驗發現：取種輪外圓面采用由多個凸起圓柱面組成的異型面結構時(圖9a)，重種現象較少，落種位置一致，因為在圓柱面的法向力作用下大蒜無法在取種輪的表面堆積，且凸起的圓柱面對將進入取種勺的大蒜進行攪動避免產生堆積和結拱堵塞現象的發生；取種輪為規則圓面時(圖9b),大蒜會在取種勺與取種輪的表面產生堆積，重種現象顯著，落種時刻提前且蒜種下落時易與上一個取種勺發生碰撞導致落種軌跡不一致。

圖9 取種輪外形對取種效果的影響

3 試驗與結果分析

3.1 試驗材料與設備

試驗設備為搭建的大蒜單粒取種試驗臺，結構如圖10所示，主要包括支架、種箱、直流驅動電機、調速器、取種輪、取種勺等。驅動電機轉速0～40 r/min連續可調。支架及種箱采用透明PVC板加工而成，便于取種過程觀察。

圖10 取種器試驗裝置

3.2 試驗方法

試驗方法參考GB/T 6973—2005《單粒(精密)播種機試驗方法》進行。由于漏充率、重充率和合格率之和為100%，因此選漏充率y1和合格率y2為評價指標，漏充率和合格率數據通過試驗過程錄像統計得到。

為確定最優單粒取種裝置工作參數，選取對取種性能影響較大的取種勺半徑R、長度L和取種輪轉速為試驗因素，漏充率和合格率為響應指標。根據前期理論計算、仿真分析及單因素試驗確定各因素取值范圍：半徑R為14.0～18.0 m，取種勺長度L為35～45 mm，取種輪直徑為130.00 mm時，取其轉速為5.0～20.0 r/min。依據Box-Benhnken中心組合設計理論進行試驗設計，試驗因素編碼如表3所示。

表3 試驗因素編碼

3.3 結果與分析

3.3.1試驗結果

試驗設計方案與結果如表4所示，其中x1、x2、x3表示各因素編碼值。

表4 試驗設計方案及結果

3.3.2回歸模型建立與顯著性檢驗

對表4中的試驗結果，運用Design-Expert 11.1.2.0數據分析軟件進行多元回歸擬合分析，得到漏充率y1和合格率y2的回歸方程

y1=6.18-4.42x1-2.02x2+1.6x3+

0.18x1x2-0.58x1x3-0.18x2x3+

(7)

y2=91.3-2.38x1-4.59x2-1.09x3-

0.85x1x2+1.05x1x3+3.13x2x3-

(8)

表5 回歸模型方差分析

注：*表示影響差異顯著(P<0.05)；** 表示影響差異極顯著(P<0.01)。

3.3.3影響因素分析

(1)各因素對漏充率的影響

固定某個因素水平為零，分析另外2個因素的相互作用對漏充率的影響，其響應曲面如圖11所示。

圖11a為取種輪轉速處于中心水平時，取種勺半徑及其長度交互作用對漏充率的影響。隨取種勺半徑及長度的增加，漏充率逐漸減小。圖11b為取種勺長度處于中心水平時，取種勺半徑及取種輪轉速交互作用對漏充率的影響。取種勺半徑一定時，隨取種輪轉速升高，漏充率逐漸增大。當取種勺半徑增大時，會減弱因取種輪轉速升高引起的漏充率升高的趨勢。圖11c為取種勺半徑處于中心水平時，取種勺長度及取種輪轉速交互作用對漏充率的影響。取種勺長度一定時，隨取種輪轉速升高，漏充率逐漸增大。當取種輪轉速一定時，隨取種勺長度增加漏充率降低。

圖11 各因素對漏充率的影響

圖12 各因素對合格率的影響

(2)各因素對合格率的影響

固定一個因素水平為零，分析另外2個因素的相互作用對合格率的影響，其響應曲面如圖12所示。

圖12a為取種輪轉速處于中心水平時，取種勺半徑及其長度交互作用對合格率的影響。取種勺長度一定時，隨著取種勺半徑增大，合格率先升后降。取種勺半徑一定時，隨著取種勺長度增大，合格率先升后降。圖12b為取種勺長度處于中心水平時，取種勺半徑及取種輪轉速交互作用對合格率的影響。取種勺半徑一定時，隨著取種輪轉速增大，合格率先升后降。取種輪轉速一定時，隨著取種勺半徑增大，合格率先升后降。圖12c為取種勺半徑處于中心水平時，取種勺長度及取種輪轉速交互作用對合格率的影響。取種勺長度一定時，隨著取種輪轉速增大，合格率先升后降。取種輪轉速一定時，隨著取種勺長度增大，合格率先升后降。

3.3.4參數優化與試驗驗證

以漏充率y1最低、合格率y2最高為目標函數，以取種勺的半徑x1、長度x2和取種輪轉速x3為約束條件，建立大蒜單粒取種裝置工作參數規劃數學模型

(9)

利用軟件Design-Expert 11.1.2.0 Optimization Numerical模塊進行優化，得最佳參數組合為：取種勺的半徑16.32 mm、長度38.46 mm，取種輪轉速9.9 r/min。此時模型預測的漏充率為5.49%，合格率為91.59%。

為驗證優化結果的準確性，設定工作參數為：取種勺半徑16.30 mm、長度38.50 mm，取種輪轉速10.0 r/min。在大蒜單粒取種試驗臺進行3次重復試驗，取平均值為試驗值，試驗結果為漏充率5.50%，合格率91.10%。結果對比表明理論優化值和試驗驗證值非常接近，回歸模型可靠。

3.4 討論

本文的研究對象是未經過篩選分級的大蒜，采用基于數理統計的理論進行取種勺的結構參數設計。由于大蒜外形尺寸變異系數較大，為滿足所有大蒜的成功取種率采用了折中的方案，因此漏充率及重播率略高，如果大蒜播種前進行篩選分級，并針對不同等級的大蒜采用不同結構參數的取種勺，取種合格率將會顯著提高。

4 結論

(1)設計了一種輪勺式大蒜單粒取種裝置，利用取種勺與取種輪形成的空腔進行單粒取種，取種勺與周邊大蒜相互擠壓作用力下大蒜處于相互競爭狀態，并最終使某粒大蒜在空腔內穩定，實現單粒取種。

(2)對取種區、輸種區和排種區的大蒜進行了受力分析，闡述輪勺式大蒜單粒取種的原理。通過離散元仿真分析對不同取種勺及取種輪的結構形狀進行了對比優化，確定了取種勺及取種輪最優結構。采用數理統計的方法確定了取種勺的尺寸區間，為取種勺設計提供了數據支持。

(3)采用Box-Benhnken中心組合方法進行了三因素三水平回歸正交試驗，分別建立漏充率和合格率多元回歸模型，以漏充率和合格率為目標，對各影響因素參數進行了優化。最優參數組合為：取種勺半徑16.30 mm、取種勺長度38.50 mm、取種輪轉速10.0 r/min。在最優參數組合下進行了臺架試驗，此時漏充率為5.50%，合格率為91.10%，與模型預測結果吻合。