劉思幸，柴巖，柳天虹，李爽，繆宏，韓非

(1. 揚(yáng)州大學(xué)機(jī)械工程學(xué)院，江蘇揚(yáng)州，225127； 2. 揚(yáng)州大學(xué)信息工程學(xué)院，江蘇揚(yáng)州，225127)

0 引言

中國(guó)是農(nóng)業(yè)大國(guó)，蔬菜產(chǎn)業(yè)是我國(guó)除糧食作物外栽培面積最廣、經(jīng)濟(jì)地位最重要的作物。而設(shè)施蔬菜種植打破了季節(jié)、氣候等條件的限制，保證了蔬菜的產(chǎn)量供應(yīng)，提高了蔬菜產(chǎn)品質(zhì)量，有利于蔬菜生產(chǎn)的產(chǎn)業(yè)化[1-2]。

設(shè)施蔬菜產(chǎn)業(yè)的快速發(fā)展，對(duì)其配套的播種裝備提出了新的要求。對(duì)于不同的設(shè)施蔬菜種植，其對(duì)應(yīng)的株距各不相同[3]。另外，為了有效利用土地，使蔬菜增產(chǎn)增效，設(shè)施蔬菜往往采用套種模式，不同蔬菜的株距、密度不一，需要每一行都有獨(dú)立的排種器根據(jù)不同株距調(diào)整轉(zhuǎn)速[4-6]。在播種機(jī)上設(shè)置合理的排種裝置，對(duì)解決設(shè)施蔬菜播種過(guò)程中的株距不一致問(wèn)題，具有極其重大的意義[7-9]。傳統(tǒng)排種驅(qū)動(dòng)方式容易由于地輪打滑造成漏播現(xiàn)象，降低株距合格率，采用步進(jìn)電機(jī)直接驅(qū)動(dòng)排種軸為該問(wèn)題提供了解決辦法，但是電機(jī)控制容易受到噪聲干擾，使靈敏度和精度降低，無(wú)法實(shí)現(xiàn)精量播種。對(duì)此，在控制系統(tǒng)中加入濾波算法，減少噪聲影響，成為播種機(jī)精量播種的重要研究方向[10-15]。隨著我國(guó)農(nóng)業(yè)機(jī)械化水平的不斷提升，對(duì)于小粒徑的設(shè)施蔬菜精量播種機(jī)的研究逐漸增多。解彬彬[16]、張靜[17]依據(jù)測(cè)速傳感器和角編碼器采集數(shù)據(jù)，通過(guò)模糊算法自動(dòng)調(diào)整排種器轉(zhuǎn)速，實(shí)現(xiàn)精量排種。王冰媛[18]提出了小粒徑種子流體噴灑播種方法，設(shè)計(jì)了噴播式小粒種子流體播種機(jī)，滿足蔬菜種子播種要求。孫靜鑫等[19-20]結(jié)合了環(huán)槽式和窩眼輪式兩種排種器優(yōu)點(diǎn)，利用排種器的護(hù)種板輔助完成排種工作。李金鳳等[21-23]采用負(fù)壓內(nèi)側(cè)吸種，正壓吹氣投種的方式，解決播種機(jī)傳動(dòng)系統(tǒng)和氣流管路復(fù)雜、布局困難的問(wèn)題。國(guó)外發(fā)達(dá)國(guó)家的設(shè)施蔬菜精量播種機(jī)主要采用氣力式排種器進(jìn)行播種，技術(shù)成熟且裝備完善，具有較高的智能化、精量化和通用化水平[24-25]。綜上，國(guó)外上述播種機(jī)體型龐大，價(jià)格高昂，配件供應(yīng)不及時(shí)，且土壤種植條件、氣候、降水等均與國(guó)內(nèi)存在明顯差異，并不適應(yīng)中國(guó)的設(shè)施蔬菜種植模式；而國(guó)內(nèi)學(xué)者對(duì)精量播種排種器結(jié)構(gòu)設(shè)計(jì)進(jìn)行了深入研究，但對(duì)排種裝置控制系統(tǒng)的無(wú)噪聲干擾方面的研究較少。

本文針對(duì)我國(guó)設(shè)施機(jī)械的現(xiàn)狀和設(shè)施蔬菜種子精量播種過(guò)程中漏種、株距不一致、株距合格率低等問(wèn)題，設(shè)計(jì)了一種基于卡爾曼濾波PID控制方法的精量排種器裝置。并進(jìn)行旋轉(zhuǎn)正交臺(tái)架試驗(yàn)，得出該排種裝置精量排種作業(yè)的最優(yōu)參數(shù)組合，驗(yàn)證所設(shè)計(jì)排種裝置的性能是否滿足設(shè)施蔬菜精量種植要求。

1 結(jié)構(gòu)和工作原理

1.1 排種器結(jié)構(gòu)

如圖1、圖2所示，播種機(jī)排種裝置主要由排種步進(jìn)電機(jī)、排種器、種箱、測(cè)速地輪和編碼器5部分組成。播種機(jī)在機(jī)架側(cè)面設(shè)有步進(jìn)電機(jī)安裝架，安裝架通過(guò)螺栓固定在機(jī)架上，步進(jìn)電機(jī)通過(guò)聯(lián)軸器與排種軸連接，步進(jìn)電機(jī)外側(cè)裝有編碼器，測(cè)速地輪置于機(jī)架前方，用于檢測(cè)牽引機(jī)前進(jìn)速度，輪架上設(shè)有拉伸彈簧，使地輪保持貼地，并具有減震功能。

(a) 排種器關(guān)鍵部件 (b) 地輪附件

圖2 排種部件結(jié)構(gòu)圖

1.2 工作原理

播種機(jī)工作時(shí)，首先由整地輪對(duì)播種之前的作業(yè)平面進(jìn)行平整，隨后由監(jiān)測(cè)地輪上的脈沖采集裝置采集到地輪轉(zhuǎn)動(dòng)生成的脈沖信號(hào)，結(jié)合地輪直徑計(jì)算出播種機(jī)前進(jìn)速度，之后控制器通過(guò)計(jì)算采集到的播種機(jī)行進(jìn)的速度信號(hào)，針對(duì)不同設(shè)施蔬菜株距要求，驅(qū)動(dòng)排種步進(jìn)電機(jī)工作，進(jìn)而驅(qū)動(dòng)排種輪的轉(zhuǎn)動(dòng)，并加以反饋控制，從而完成整個(gè)排種作業(yè)。

2 控制系統(tǒng)設(shè)計(jì)

2.1 總體設(shè)計(jì)

排種控制系統(tǒng)由主控制器、電平轉(zhuǎn)換單元、執(zhí)行單片機(jī)、步進(jìn)電機(jī)驅(qū)動(dòng)器、步進(jìn)電機(jī)、編碼器和排種軸組成。主控制器和執(zhí)行器通過(guò)RS-485信號(hào)進(jìn)行排種指令的發(fā)送和執(zhí)行的反饋，主控制器發(fā)出排種數(shù)字信號(hào)，執(zhí)行器將接收到的信號(hào)轉(zhuǎn)為脈沖信號(hào)下發(fā)給驅(qū)動(dòng)器，進(jìn)而驅(qū)動(dòng)步進(jìn)電機(jī)帶動(dòng)排種軸轉(zhuǎn)動(dòng)相應(yīng)角度，設(shè)施種子沿排種盤下落，完成一次排種作業(yè)。其控制系統(tǒng)總體框圖如圖3所示，CAN總線通訊和步進(jìn)電機(jī)驅(qū)動(dòng)電路設(shè)計(jì)圖如圖4所示。

圖3 控制系統(tǒng)總體框圖

(a) CAN總線通訊電路

(b) 步進(jìn)電機(jī)驅(qū)動(dòng)電路

2.2 控制原理

設(shè)施蔬菜精量播種需要對(duì)排種器排種進(jìn)行精準(zhǔn)控制，而在實(shí)際播種過(guò)程中，地塊不平整、打滑、磁場(chǎng)等因素嚴(yán)重影響排種器排種精度，引入一種卡爾曼濾波PID控制方法，對(duì)步進(jìn)電機(jī)進(jìn)行實(shí)時(shí)反饋控制。播種過(guò)程中，控制器根據(jù)用戶設(shè)定株距值并結(jié)合測(cè)速地輪采集的實(shí)時(shí)速度數(shù)據(jù)，計(jì)算出當(dāng)前電機(jī)所需脈沖數(shù)，與編碼器采集的實(shí)時(shí)電機(jī)轉(zhuǎn)速進(jìn)行分析和判斷，對(duì)于產(chǎn)生的數(shù)據(jù)偏差通過(guò)控制方法及時(shí)修正，從而產(chǎn)生當(dāng)前應(yīng)該輸出的控制信號(hào)，經(jīng)過(guò)執(zhí)行機(jī)構(gòu)施加到步進(jìn)電機(jī)上，實(shí)現(xiàn)精量排種的目的。

輸入脈沖數(shù)與理論株距計(jì)算關(guān)系如下：株距調(diào)節(jié)依靠步進(jìn)電機(jī)獨(dú)立驅(qū)動(dòng)實(shí)現(xiàn)，設(shè)理論設(shè)定株距為s(cm)，排種器轉(zhuǎn)速即步進(jìn)電機(jī)轉(zhuǎn)速為n1(r/min)，排種器每轉(zhuǎn)排出種子數(shù)為k，則在時(shí)間t(s)內(nèi)播種機(jī)前進(jìn)的距離

L=n1tks

(1)

因此要滿足在時(shí)間t內(nèi)排出種子的株距要求，牽引機(jī)需要前進(jìn)同樣的距離L(m)。已知牽引機(jī)前進(jìn)速度由測(cè)速地輪獲得，具體為由安裝在地輪上的編碼器采集到時(shí)間t內(nèi)的脈沖數(shù)計(jì)算而得，設(shè)編碼器的分辨率為P(即為編碼器旋轉(zhuǎn)一周輸出的脈沖數(shù))，時(shí)間t內(nèi)編碼器輸出的脈沖數(shù)為ΔP，則編碼器測(cè)得地輪的轉(zhuǎn)速

(2)

設(shè)測(cè)速地輪的直徑為d(mm)，可得時(shí)間t內(nèi)地輪走過(guò)的距離LZ(m)，有

LZ=πdnt

(3)

可知LZ即牽引機(jī)前進(jìn)的距離，即LZ=L，將式(1)、式(3)聯(lián)立，可得

πdn=n1ks

(4)

進(jìn)而可得牽引機(jī)前進(jìn)速度V(m/s)，有

(5)

步進(jìn)電機(jī)轉(zhuǎn)動(dòng)靠輸入脈沖信號(hào)實(shí)現(xiàn)，每輸入一個(gè)脈沖信號(hào)給電機(jī)，轉(zhuǎn)子轉(zhuǎn)過(guò)一個(gè)角度(步距角)，設(shè)步距角為β，輸入脈沖數(shù)為kc，則步進(jìn)電機(jī)轉(zhuǎn)速

(6)

將式(5)、式(6)聯(lián)立，可知要滿足理論設(shè)定株距的要求，需要向步進(jìn)電機(jī)輸入的脈沖數(shù)

(7)

2.3 控制算法實(shí)現(xiàn)

傳統(tǒng)PID算法存在參數(shù)整定不良、超調(diào)大、適應(yīng)性差、控制性能不佳、干擾后調(diào)整時(shí)間長(zhǎng)、穩(wěn)態(tài)誤差大等問(wèn)題，不適用于精量播種系統(tǒng)調(diào)節(jié)，并且由于噪聲的影響，控制器的參數(shù)難以自動(dòng)調(diào)節(jié)，不能達(dá)到理想的控制效果。在傳統(tǒng)PID控制系統(tǒng)中添加極小極快無(wú)量綱高斯白噪聲后，并選取其調(diào)節(jié)參數(shù)分別為10、5、1，其MATLAB中Simulink控制模型搭建和仿真結(jié)果如圖5所示，可以看出，此時(shí)由于微分項(xiàng)的影響，在系統(tǒng)達(dá)到穩(wěn)定值后處于高頻波動(dòng)狀態(tài)，不利于電機(jī)反饋控制，難以實(shí)現(xiàn)精量排種。

(a) Simulink控制模型

卡爾曼濾波是一種高效率的遞歸濾波器，它能夠從一系列的不完全及包含噪聲的測(cè)量中，估計(jì)動(dòng)態(tài)系統(tǒng)的狀態(tài)。其可以在貝葉斯濾波的基礎(chǔ)上做出如下假設(shè)，其中，式(8)為預(yù)測(cè)方程，式(9)為觀測(cè)方程。

XK=FXK-1+QK

(8)

YK=hXK+RK

(9)

式中：QK——預(yù)測(cè)噪聲，QK服從高斯分布；

RK——觀測(cè)噪聲，RK服從高斯分布。

根據(jù)貝葉斯濾波的預(yù)測(cè)步和更新步方程可計(jì)算出卡爾曼濾波的預(yù)測(cè)步和更新步方程為

(10)

K——卡爾曼增益。

對(duì)傳統(tǒng)PID控制器進(jìn)行卡爾曼濾波融合后，在MATLAB中編寫控制程序并在Simulink中搭建控制模型完成仿真測(cè)試，其控制仿真模型如圖6所示，系統(tǒng)響應(yīng)結(jié)果如圖7所示，由結(jié)果可以得出，采用卡爾曼濾波后，對(duì)于相同量級(jí)的高斯白噪聲干擾，其信號(hào)的諧波幅值明顯降低，穩(wěn)定性提高，趨近于理想信號(hào)。顯然，卡爾曼濾波的引入，對(duì)干擾起到了抑制與降低，表明該系統(tǒng)有著較好的噪聲抑制能力，這對(duì)于電機(jī)精確控制及排種器精量播種的實(shí)現(xiàn)起到了關(guān)鍵作用。

圖6 控制仿真模型

(a) 融合卡爾曼濾波前

3 排種試驗(yàn)

3.1 評(píng)價(jià)指標(biāo)

試驗(yàn)選取設(shè)施蔬菜種子中具有代表性的丸?；S瓜種子，千粒質(zhì)量為28 g，含水率為7%。根據(jù)GB/T 6973—2005《單粒(精密)播種機(jī)試驗(yàn)方法》規(guī)定的試驗(yàn)方法和指標(biāo)，選取株距合格率、重播率、漏播率和株距變異系數(shù)作為排種性能評(píng)價(jià)指標(biāo)，進(jìn)行設(shè)施蔬菜智能精量播種機(jī)排種試驗(yàn)評(píng)價(jià)。

1) 株距合格率

(11)

式中：n1——株距合格數(shù)；

N′——區(qū)間數(shù)。

2) 重播率

(12)

式中：n2——重播數(shù)。

3) 漏播率

(13)

式中：n0——漏播數(shù)。

4) 株距變異系數(shù)

(14)

(15)

C=σ×100%

(16)

式中：C——株距變異系數(shù)，%；

σ——標(biāo)準(zhǔn)差。

3.2 試驗(yàn)方案與結(jié)果分析

3.2.1 單因素試驗(yàn)

通過(guò)實(shí)地調(diào)試，確定排種盤轉(zhuǎn)速和播種機(jī)行走速度是影響設(shè)施蔬菜種子株距水平的主要因素。搭建排種試驗(yàn)臺(tái)架，如圖8所示。

圖8 臺(tái)架試驗(yàn)圖

為了得到理想穩(wěn)定的株距水平，在排種試驗(yàn)臺(tái)上對(duì)排種盤轉(zhuǎn)速和作業(yè)速度進(jìn)行單因素試驗(yàn)，確定排種盤轉(zhuǎn)速20 r/min、行走速度2 km/h 為排種器單因素試驗(yàn)初步優(yōu)選值。

3.2.2 正交旋轉(zhuǎn)多因素試驗(yàn)

為尋求排種盤轉(zhuǎn)速和作業(yè)速度兩個(gè)因素對(duì)于株距水平影響的最優(yōu)參數(shù)組合，以株距合格率、重播率、漏播率和株距變異系數(shù)作為評(píng)價(jià)指標(biāo)，進(jìn)行兩因素五水平二次回歸正交旋轉(zhuǎn)組合試驗(yàn)。根據(jù)單因素試驗(yàn)結(jié)果確定試驗(yàn)因素水平設(shè)置，試驗(yàn)因素水平編碼表如表1所示，試驗(yàn)方案及試驗(yàn)結(jié)果如表2所示。

表1 試驗(yàn)因素水平Tab. 1 Experimental factors and levels

運(yùn)用Design Expert 8.0.6軟件對(duì)試驗(yàn)結(jié)果進(jìn)行方差分析，并進(jìn)行多元回歸擬合，得到株距合格率Y1、重播率Y2、漏播率Y3和株距變異系數(shù)Y4各試驗(yàn)指標(biāo)的回歸方程。其方差分析結(jié)果如表3～表6所示。

Y1=94.91-0.57X1-1.59X2+0.13X1X2-

0.31X12+0.38X22

(17)

Y2=2.47-0.20X1-1.20X2+0.019X1X2-

0.12X12-0.30X22

(18)

Y3=2.88+0.21X1+1.40X2-0.15X1X2-

0.20X12-0.077X22

(19)

Y4=14.43+0.77X1+2.10X2-0.59X1X2-

0.42X12-0.55X22

(20)

由表3～表6的回歸模型方差分析結(jié)果可知，四個(gè)回歸模型的P值均小于0.000 1，表明回歸模型高度顯著；在試驗(yàn)因素中，由排種盤轉(zhuǎn)速和行走速度的P值可判斷行走速度相比于排種盤轉(zhuǎn)速對(duì)4個(gè)試驗(yàn)指標(biāo)的影響更為顯著，株距合格率、重播率、漏播率和株距變異系數(shù)回歸模型中的P值均小于0.05，表明4個(gè)回歸項(xiàng)在回歸模型中交互影響顯著。對(duì)于重播率Y2，回歸項(xiàng)X22也影響極顯著；對(duì)于株距變異系數(shù)Y4，回歸項(xiàng)X1X2也影響顯著。四模型的決定系數(shù)R2均接近1，變異系數(shù)分別為0.83%、0.96%、0.95%、0.87%，說(shuō)明該株距合格率、重播率、漏播率和株距變異系數(shù)擬合回歸模型具有較高的可靠性。

表3 株距合格率方差分析Tab. 3 Variance analysis of qualified rate of plant spacing

表4 重播率方差分析Tab. 4 Variance analysis of multiple rate

表5 漏播率方差分析Tab. 5 Variance analysis of missing rate

表6 株距變異系數(shù)方差分析Tab. 6 Variance analysis of coefficient of

經(jīng)過(guò)數(shù)據(jù)處理和分析，通過(guò)Design Expert 8.0.6軟件得到各因素對(duì)于試驗(yàn)指標(biāo)的響應(yīng)曲面。排種盤轉(zhuǎn)速和行走速度對(duì)株距合格率、重播率、漏播率和株距變異系數(shù)的響應(yīng)曲面分別如圖9～圖12所示，行走速度由1 km/h提高到3 km/h過(guò)程中，排種盤轉(zhuǎn)速越快，株距合格率和重播率越低；當(dāng)行走速度由快變慢時(shí)，株距合格率和重播率隨著排種盤轉(zhuǎn)速的減小而提高。當(dāng)排種盤轉(zhuǎn)速由10 r/min增加到30 r/min過(guò)程中，行走速度越高，漏播率和株距變異系數(shù)越高；當(dāng)排種盤轉(zhuǎn)速由高變低時(shí)，漏播率和株距變異系數(shù)隨著行走速度的減小而減小。

圖9 株距合格率指數(shù)的雙因素響應(yīng)曲面

圖10 重播率指數(shù)的雙因素響應(yīng)曲面

圖11 漏播率指數(shù)的雙因素響應(yīng)曲面

圖12 株距變異系數(shù)指數(shù)的雙因素響應(yīng)曲面

隨著行走速度的提高，株距合格率和重播率呈下降趨勢(shì)，而漏播率和株距變異系數(shù)呈上升趨勢(shì)；在同一設(shè)定速度下，隨著排種盤轉(zhuǎn)速的提高，株距合格率和重播率呈下降趨勢(shì)，而漏播率和株距變異系數(shù)呈上升趨勢(shì)。排種盤轉(zhuǎn)速一致時(shí)，行走速度減小，重播率的增幅較大；行進(jìn)速度增大，漏播率的增幅相對(duì)較大。基于交互效應(yīng)分析可知，排種盤轉(zhuǎn)速和行走速度同時(shí)減小時(shí)，株距合格率指數(shù)、漏播率指數(shù)和株距變異系數(shù)指數(shù)顯優(yōu)；排種盤轉(zhuǎn)速和行走速度同時(shí)增大時(shí)，重播率指數(shù)顯優(yōu)。

3.3 參數(shù)優(yōu)化

為了得到各試驗(yàn)因素的最佳水平組合，運(yùn)用Design Expert 8.0.6軟件中的Optimization功能對(duì)回歸模型進(jìn)行求解，以株距合格率最大、重播率最小、漏播率最小、株距變異系數(shù)最小為條件，求解模型得到的排種盤轉(zhuǎn)速為10 r/min，行走速度為1.6 km/h時(shí)，滿足上述條件，其株距合格率為95.9%，重播率為2.9%，漏播率為1.9%，株距變異系數(shù)為12.1%，均符合標(biāo)準(zhǔn)規(guī)定。

4 結(jié)論

1) 針對(duì)設(shè)施蔬菜排種器存在的漏播、重播、精度低等問(wèn)題，對(duì)排種器裝置進(jìn)行了結(jié)構(gòu)和控制系統(tǒng)優(yōu)化設(shè)計(jì)，利用步進(jìn)電機(jī)直接驅(qū)動(dòng)排種器進(jìn)行排種作業(yè)，通過(guò)控制系統(tǒng)實(shí)現(xiàn)精量播種，實(shí)驗(yàn)得出影響設(shè)施蔬菜精量播種的主要因素為排種盤的轉(zhuǎn)速和機(jī)具的行走速度。

2) 對(duì)傳統(tǒng)PID控制器進(jìn)行優(yōu)化設(shè)計(jì)，融合卡爾曼濾波，仿真結(jié)果表明，其可以有效減少噪聲對(duì)電機(jī)的干擾，對(duì)提升系統(tǒng)控制過(guò)程的平穩(wěn)性起到重要作用。

3) 通過(guò)分析排種盤轉(zhuǎn)速和行走輪行走速度對(duì)排種質(zhì)量的單因素試驗(yàn)，得到理想排種質(zhì)量的各參數(shù)合理變化范圍：排種盤轉(zhuǎn)速為10～30 r/min，行走輪行走速度為1～3 km/h。

4) 采用兩因素五水平正交試驗(yàn)，建立各個(gè)因素與評(píng)價(jià)指標(biāo)間的回歸模型，得到影響株距合格率、重播率、漏播率和株距變異系數(shù)的因素為：行走速度>排種盤轉(zhuǎn)速。確定最佳參數(shù)組合：當(dāng)排種盤轉(zhuǎn)速為10 r/min，行進(jìn)速度為1.6 km/h時(shí)，株距合格率為95.9%，重播率為2.9%，漏播率為1.9%，株距變異系數(shù)為12.1%，滿足設(shè)施蔬菜的精量播種要求。