周 俊 孟一猛 張 娜 嚴(yán)正紅

(南京農(nóng)業(yè)大學(xué)江蘇省智能化農(nóng)業(yè)裝備重點(diǎn)實(shí)驗(yàn)室， 南京 210031)

機(jī)械手不同抓取控制方式對(duì)番茄機(jī)械損傷的影響分析

周 俊 孟一猛 張 娜 嚴(yán)正紅

(南京農(nóng)業(yè)大學(xué)江蘇省智能化農(nóng)業(yè)裝備重點(diǎn)實(shí)驗(yàn)室， 南京 210031)

為了分析機(jī)械手不同抓取控制方式對(duì)番茄造成的機(jī)械損傷。通過Burgers粘彈性模型描述番茄的變形特征，利用蠕變實(shí)驗(yàn)獲得番茄樣本的粘彈性參數(shù)。選擇勻減速、先勻速后減速以及變減速作為3種典型的抓取控制方式，并利用一次函數(shù)、巴特沃斯幅度平方函數(shù)和指數(shù)函數(shù)3種函數(shù)曲線描述末端執(zhí)行器實(shí)際抓取過程中不同抓取控制方式的速度變化。在此基礎(chǔ)上，對(duì)番茄塑性變形進(jìn)行求解，分析在抓取時(shí)間一定時(shí)，不同抓取控制方式下番茄的塑性變形規(guī)律，并以最小塑性變形為依據(jù)，確定最優(yōu)的抓取策略。結(jié)果表明：抓取控制方式3為最優(yōu)抓取控制方式，當(dāng)抓取速度v0=1 mm/s，抓取時(shí)間t0=1 s時(shí)，產(chǎn)生的最小塑性變形最小，紅熟前期和紅熟中期番茄的最小塑性變形分別為0.002 6 mm和0.009 8 mm。以抓取控制方式3進(jìn)行抓取控制實(shí)驗(yàn)，實(shí)驗(yàn)數(shù)據(jù)和理論數(shù)據(jù)的吻合度較高，決定系數(shù)R2為0.99，驗(yàn)證了抓取控制方式3作為最優(yōu)抓取控制方式的合理性和可行性，為優(yōu)化機(jī)械手抓取過程提供了參考。

番茄； 采摘； 機(jī)械手； 抓取控制方式； 機(jī)械損傷

引言

應(yīng)用農(nóng)業(yè)機(jī)器人完成果蔬采摘、分揀等作業(yè)的研究越來越多[1-2]，末端執(zhí)行器抓取果蔬對(duì)象是其中一個(gè)重要的作業(yè)環(huán)節(jié)。末端執(zhí)行器抓取果實(shí)方式不同，對(duì)果實(shí)產(chǎn)生的機(jī)械損傷程度也不盡相同。在抓取過程中，果實(shí)因受到擠壓、沖擊、振動(dòng)等多種載荷作用而產(chǎn)生變形，主要形成以塑性變形為主的現(xiàn)時(shí)損傷和以粘彈性變形為主的延遲損傷[3]，而其中塑性變形對(duì)果品品質(zhì)的影響最為關(guān)鍵[4]。大多數(shù)的果實(shí)為粘彈性體，具有固體和流體的雙重特性[5]，利用粘彈性模型及參數(shù)對(duì)果實(shí)的變形進(jìn)行分析，可以為加工與采收機(jī)械的設(shè)計(jì)提供理論基礎(chǔ)，也可以減少與控制作業(yè)過程中的機(jī)械損傷，為提高產(chǎn)品質(zhì)量提供技術(shù)支持。

目前國(guó)內(nèi)外學(xué)者已通過多種方式對(duì)抓取過程中果實(shí)的變形損傷進(jìn)行研究。LI等[6]對(duì)同一番茄的不同部位進(jìn)行抓取試驗(yàn)，對(duì)比不同抓取位置產(chǎn)生的峰值力和最大變形，為減少機(jī)器人抓取的機(jī)械損傷提供依據(jù)。AHMADI等[7]建立了蘋果的非線性粘彈性模型，通過仿真獲得不同抓取速度下蘋果與剛性平板接觸過程中的變形隨時(shí)間的變化規(guī)律。HIGASHIMORI等[8]基于獲得流變物體的四單元粘彈性模型參數(shù)分析施加的抓取力與變形之間的關(guān)系，在避免對(duì)象過度變形的條件下減少夾持作業(yè)時(shí)間。文獻(xiàn)[9-10]利用有限元的方法分別對(duì)蘋果和柑橘與不同指面類型的夾持手指的接觸進(jìn)行了研究，對(duì)比自動(dòng)抓取過程中不同指面類型作用下接觸應(yīng)力和果實(shí)變形的差異。劉繼展等[11]以果實(shí)的粘彈性模型為基礎(chǔ)，研究了抓取速度與碰撞峰值力和變形之間的關(guān)系，為實(shí)現(xiàn)果實(shí)快速夾持作業(yè)提供借鑒。馬曉麗等[12]對(duì)楊梅果開展壓縮試驗(yàn)，結(jié)合力-變形關(guān)系曲線對(duì)其力學(xué)特性進(jìn)行分析，探討機(jī)械采摘條件下楊梅果的機(jī)械損傷機(jī)理。

但現(xiàn)有研究中依據(jù)塑性變形這一影響果品品質(zhì)的關(guān)鍵因素來分析機(jī)械手抓取果蔬所造成的機(jī)械損傷的研究較為少見。因此，本文以番茄為作業(yè)對(duì)象，利用Burgers粘彈性模型表征其壓縮變形特征。通過蠕變實(shí)驗(yàn)獲得番茄樣本的粘彈性參數(shù)，對(duì)3種典型抓取控制方式下果實(shí)的塑性變形進(jìn)行理論分析，并依據(jù)最小塑性變形原則確定最優(yōu)抓取控制方式，以期減少機(jī)械手在抓取過程中對(duì)果實(shí)產(chǎn)生的機(jī)械損傷。

1 番茄抓取塑性變形

對(duì)于固體農(nóng)業(yè)物料的果蔬而言，可將其受力時(shí)的形變分為彈性變形和塑性變形[8,13-14]。機(jī)械手起始抓取階段，果蔬組織產(chǎn)生的初始小變形為彈性變形，但變形能夠完全恢復(fù)；塑性變形指在施加的外力撤銷后，不能恢復(fù)而保留下來的那部分變形，是物體無(wú)法恢復(fù)原狀的一種物理現(xiàn)狀。物體最終的形狀由塑性變形直接決定，同時(shí)塑性變形也是影響果品品質(zhì)的關(guān)鍵因素。因此，本文將塑性變形作為番茄機(jī)械損傷的衡量指標(biāo)，通過分析接觸產(chǎn)生的塑性變形來評(píng)判當(dāng)前抓取控制方式的優(yōu)劣，從而選擇合適的抓取控制方式，減少被抓物體的機(jī)械損傷。

1.1 番茄粘彈性數(shù)學(xué)模型

包括番茄在內(nèi)的農(nóng)業(yè)物料大多為粘彈性體，具有流變特性。果蔬的流變特性主要研究的是果蔬在外力作用下產(chǎn)生的流動(dòng)與變形，利用粘彈性模型描述其流變特性，可以更好地分析果蔬在受載條件下的變形特征。粘彈性模型最基本的有Maxwell模型和Kelvin模型，但對(duì)于果蔬等農(nóng)業(yè)物料，利用簡(jiǎn)單的兩種模型不足以較為完整地表征果蔬復(fù)雜的力學(xué)特性。將兩種簡(jiǎn)單的組合方式再串聯(lián)起來就形成了伯格斯模型(Burgers)，如圖1所示。

圖1 Burgers模型Fig.1 Burgers model

Burgers模型蠕變加載過程的微分方程為

(1)

式中D(t)——變形量，mm

t——時(shí)間，s

F0——恒定載荷，N

k1——瞬時(shí)彈性系數(shù)，N/mm

k2——延遲彈性系數(shù)，N/mm

c1——串聯(lián)粘性系數(shù)，N·s/mm

c2——并聯(lián)粘性系數(shù)，N·s/mm

Burgers模型可以表示瞬時(shí)彈性變形、延遲彈性變形以及永久塑性變形等多種變形，且根據(jù)前人對(duì)果品特性的研究，Burgers模型能夠較好地描述果品的流變特性[15-16]。因此，本文決定采用伯格斯模型表征番茄果實(shí)的變形特征。

1.2 抓取塑性變形定義

當(dāng)外力f施加在物體上時(shí)，伯格斯模型表征的番茄變形如圖2所示。x、xi及fi分別表示物體的總變形、每一層變形以及彈性單元層作用力(i=1)，粘彈性平行單元層作用力(i=2)和粘性單元作用層作用力(i=3)，kh為機(jī)械手指剛度，xh為從接觸到穩(wěn)定抓取過程中機(jī)械手的閉合位移，則番茄受力和變形的表達(dá)式為[17]

圖2 番茄伯格斯模型受力變形表征Fig.2 Force-deformation characterization of Burgers model

(2)

x(t)=x1(t)

(3)

其中

f1(t)=k1(x1(t)-x2(t))

(4)

(5)

(6)

又因

x(t)=xh(t)-f(t)/kh

(7)

根據(jù)式(2)～(7)，移除fi和xi，可得抓取過程中的力-位移關(guān)系模型為

(8)

其中

B2=c1c2/k2

(9)

B1=c1

(10)

A2=c1c2(k1+kh)/(k1k2kh)

(11)

A1=[c1k1k2+kh(k1c2+k1c1+k2c1)]/(k1k2kh)

(12)

伯格斯模型表征的粘彈性物體的塑性變形主要由粘性系數(shù)c1決定[18]，則番茄的塑性變形為

(13)

式中T——末端執(zhí)行器和果實(shí)的接觸時(shí)間

根據(jù)式(8)和式(13)可知，對(duì)番茄的塑性變形進(jìn)行理論分析需要獲得其粘彈性參數(shù)。因此，本文利用式(1)的蠕變方程，對(duì)番茄開展恒定載荷蠕變實(shí)驗(yàn)獲得其粘彈性參數(shù)[19-21]。

2 番茄粘彈性參數(shù)獲取

2.1 蠕變實(shí)驗(yàn)

選擇番茄品種為皖粉2號(hào)的紅熟前期(著紅面積10%～40%)、紅熟中期(著紅面積40%～60%)番茄各10個(gè)，番茄的赤道直徑在62～80 mm之間，質(zhì)量在166～252 g之間。實(shí)驗(yàn)儀器包括美國(guó)FTC公司的TMS-PRO質(zhì)構(gòu)儀、P75型探頭(直徑75 mm)、工控機(jī)、游標(biāo)卡尺、數(shù)顯電子秤。

對(duì)各成熟期的番茄進(jìn)行編號(hào)，將測(cè)試樣品的果柄水平放置在質(zhì)構(gòu)儀上下探頭間，并通過底座固定，逐一進(jìn)行實(shí)驗(yàn)。設(shè)定探頭加載速度為60 mm/min，數(shù)據(jù)采集頻率為10 Hz，起始力為0.6 N，加載時(shí)間30 s。每個(gè)樣品在果腰處隨機(jī)選取3個(gè)樣本點(diǎn)，并且對(duì)這3個(gè)樣本點(diǎn)分別設(shè)定恒定壓力3、5、7 N進(jìn)行測(cè)試，共進(jìn)行60次蠕變實(shí)驗(yàn)。

2.2 模型參數(shù)獲取

蠕變實(shí)驗(yàn)可以獲得番茄在恒載作用下變形隨時(shí)間變化關(guān)系的蠕變曲線，利用Origin 8.0數(shù)據(jù)分析軟件中的非線性擬合工具箱，通過式(1)對(duì)番茄蠕變數(shù)據(jù)進(jìn)行曲線擬合，得到蠕變模型的粘彈性參數(shù)及相關(guān)系數(shù)。表1為從20個(gè)樣本中隨機(jī)選擇2個(gè)番茄樣本的粘彈性參數(shù)與相關(guān)系數(shù)，并計(jì)算其平均值。決定系數(shù)R2均高于0.992，其余樣本的擬合情況與之類似，驗(yàn)證了所選模型的有效性。本文利用番茄樣本粘彈性參數(shù)的平均值進(jìn)行后續(xù)抓取控制方式的理論和仿真分析。

表1 粘彈性參數(shù)及決定系數(shù)Tab.1 Viscoelastic parameters and determination coefficient

3 不同抓取控制方式對(duì)塑性變形的影響

3.1 3種典型抓取控制方式

對(duì)于采摘機(jī)器人而言，在末端執(zhí)行器從接觸果實(shí)到停止抓取的整個(gè)階段，末端執(zhí)行器閉合速度均為從某一抓取速度最終減小為零。雖然該過程中速度變化形式存在無(wú)數(shù)種可能，但勻減速、先勻速后減速以及變減速3種方式為典型速度變化形式，即：①末端執(zhí)行器以設(shè)定的抓取速度勻速閉合，接觸物體后開始勻減速，直至停止。②末端執(zhí)行器接觸果實(shí)后，仍保持設(shè)定的抓取速度勻速閉合，一段時(shí)間后減速最終停止。③末端執(zhí)行器接觸果實(shí)后以加速度不斷變化的方式減速運(yùn)行，直至停止閉合。因此，本文選擇勻減速、先勻速后減速以及變減速3種具有代表性的抓取控制方式進(jìn)行分析。通常情況下人工完成單次穩(wěn)定抓取動(dòng)作的時(shí)間為1 s，目前，日本采摘機(jī)械手完成穩(wěn)定抓取的時(shí)間大約為4 s[22]，為了分析在一定的采摘時(shí)間內(nèi)，不同的抓取控制方式對(duì)果實(shí)機(jī)械損傷的影響，本節(jié)對(duì)3種典型抓取控制方式下產(chǎn)生的塑性變形規(guī)律進(jìn)行分析，對(duì)比3種典型抓取控制方式對(duì)番茄機(jī)械損傷的影響。

結(jié)合現(xiàn)有實(shí)驗(yàn)平臺(tái)(WSG 50機(jī)械手)的配置參數(shù)，選擇一次函數(shù)、巴特沃斯幅度平方函數(shù)和指數(shù)函數(shù)3種函數(shù)曲線分別描述勻減速、先勻速后減速和變減速3種典型的抓取控制方式，3種函數(shù)表征的抓取速度隨時(shí)間變化的關(guān)系如圖3所示。

圖3 3種抓取控制方式的速度隨時(shí)間變化曲線Fig.3 Diagram of velocity change over time of different grasping control modes

根據(jù)實(shí)驗(yàn)平臺(tái)的實(shí)際參數(shù)配置，確定3種典型抓取控制方式對(duì)應(yīng)的速度隨時(shí)間變化的函數(shù)方程為

(14)

式中v0——設(shè)定的抓取速度

v1(t)——一次函數(shù)下的速度

v2(t)——巴特沃斯幅度平方函數(shù)下的速度

v3(t)——指數(shù)函數(shù)下的速度

t0——穩(wěn)定抓取時(shí)間

3.2 塑性變形理論計(jì)算

由于塑性變形是影響果品品質(zhì)的關(guān)鍵因素，因此決定利用果實(shí)抓取過程中的塑性變形作為評(píng)判抓取控制方式優(yōu)劣的依據(jù)。設(shè)定機(jī)械手抓取速度為v0，依據(jù)果實(shí)質(zhì)量m，可得到當(dāng)前作業(yè)對(duì)象的停止抓取力F0=mg/(2μ)(μ為摩擦因數(shù)，實(shí)驗(yàn)測(cè)得μ=0.22)。依據(jù)2.2節(jié)隨機(jī)選取的番茄樣本及其粘彈性參數(shù)平均值進(jìn)行理論計(jì)算，利用式(8)、(14)進(jìn)行微分方程求解獲得不同抓取控制方式下番茄的受力隨時(shí)間變化關(guān)系。將其代入式(13)可以獲得在一定抓取時(shí)間下，不同抓取控制方式產(chǎn)生的塑性變形量。理論分析的數(shù)學(xué)模型求解參數(shù)如表2所示。

表2 模型求解參數(shù)Tab.2 Solving parameters of model

為了進(jìn)一步說明求解過程，本文以抓取控制方式1在抓取速度v0=1 mm/s，t0=1 s時(shí)的抓取過程為例，對(duì)紅熟前期番茄樣本塑性變形的理論分析過程進(jìn)行詳細(xì)介紹，求解步驟如下：

(15)

(2)式(15)為二階常系數(shù)非齊次線性微分方程，屬于f(t)=eλtPm(t)的類型。由函數(shù)解的性質(zhì)知，求該式的通解可以歸結(jié)為求對(duì)應(yīng)齊次方程

(16)

的通解和非齊次方程式(15)的一個(gè)非零特解y*。

(3)式(16)對(duì)應(yīng)的特征方程為

1 648.17r2+355.47r+1=0

(17)

f*(t)=c1e-0.002 9t+c2e-0.212 8t

(18)

由于λ=0不是特征方程式(17)的根，所以應(yīng)設(shè)特解為

y*=b0t+b1

(19)

將其代入式(15)，求得其1個(gè)特解為

y*=-1 075.41t+376 577.78

(20)

(4)根據(jù)機(jī)械手的穩(wěn)定抓取過程參數(shù)，確定方程(15)的初始條件為：f(0)=0，f(t0)=F0，將初始條件代入式(15)求得其通解為

f(t)=-12 630e-0.002 9t-11.462 6e-0.212 8t-
1 075.41t+376 577.78

(21)

則該抓取過程對(duì)應(yīng)的塑性變形表達(dá)式為

(22)

式(21)即為當(dāng)番茄的粘彈性參數(shù)為k1=4.11 N/mm，k2=12.29 N/mm，c1=1 075.41 N·s/mm，c2=77.41 N·s/mm，穩(wěn)定抓取時(shí)間t0=1 s，抓取速度v0=1 mm/s時(shí)，抓取控制方式1的整個(gè)穩(wěn)定抓取過程中抓取力隨時(shí)間的變化關(guān)系。通過Matlab仿真獲得式(21)對(duì)應(yīng)的力隨時(shí)間變化曲線如圖4所示，將式(21)代入塑性變形表達(dá)式(13)可得抓取控制方式1在該過程中的塑性變形表達(dá)式，求出對(duì)應(yīng)的塑性變形Xp=0.002 7 mm。

同理，利用上述微分方程求解方法對(duì)抓取控制方式2和抓取控制方式3進(jìn)行計(jì)算，得到對(duì)應(yīng)的t0=1，v0=1 mm/s時(shí)的抓取力隨時(shí)間變化關(guān)系如圖5所示，求得對(duì)應(yīng)抓取控制方式2和抓取控制方式3的塑性變形Xp分別為0.006 2、0.002 6 mm。根據(jù)以上求解過程，結(jié)合Matlab工程計(jì)算軟件對(duì)不同成熟度番茄樣本在通常的機(jī)械手抓取時(shí)間t0為1～4 s，抓取速度v0為1～20 mm/s情況下的塑性變形進(jìn)行計(jì)算，番茄紅熟前期和紅熟中期的求解結(jié)果如圖6和圖7所示。

圖6、圖7為采取不同抓取控制方式，在不同抓取時(shí)間t0下，不同成熟度番茄產(chǎn)生的塑性變形隨抓取速度的變化規(guī)律。從圖6、圖7中可以看出，對(duì)同一成熟度的番茄樣本而言，在同一抓取時(shí)間下，隨著抓取速度的增大，番茄的塑性變形增大；在不同抓取時(shí)間下，隨著抓取時(shí)間的增大，機(jī)械手與番茄的接觸時(shí)間延長(zhǎng)，蠕變現(xiàn)象愈明顯，產(chǎn)生的塑性變形越大。對(duì)于不同成熟度的番茄樣本而言，其塑性變形隨抓取速度和抓取時(shí)間的變化規(guī)律一致，且成熟度越高，相同抓取時(shí)間和抓取速度條件下產(chǎn)生塑性變形越大。

圖4 抓取控制方式1的力隨時(shí)間變化曲線Fig.4 Changing curve of force over time of the first grasping control mode

圖5 抓取控制方式2和抓取控制方式3的力隨時(shí)間變化曲線Fig.5 Changing curves of force over time of the second and third grasping control modes

圖6 3種抓取控制方式下紅熟前期番茄的塑性變形規(guī)律Fig.6 Plastic deformations of early maturity tomato with three grasping control modes

圖7 3種抓取控制方式下紅熟中期番茄的塑性變形規(guī)律Fig.7 Plastic deformations of middle maturity tomato with three grasping control modes

通過對(duì)比3種典型抓取控制方式，可以得出：在相同抓取速度下，抓取控制方式2產(chǎn)生的塑性變形最大，抓取控制方式3產(chǎn)生的塑性變形最??；當(dāng)抓取速度v0=20 mm/s，t0=4 s時(shí)，抓取控制方式2產(chǎn)生的塑性變形Xp最大，紅熟前期和紅熟中期番茄的最大塑性變形分別為0.105 mm和0.119 mm；抓取速度v0=1 mm/s，t0=1 s時(shí)，抓取控制方式3產(chǎn)生的塑性變形最小，紅熟前期和紅熟中期番茄的最小塑性變形分別為0.002 6 mm和0.009 8 mm。因此，為了減少作業(yè)過程產(chǎn)生的塑性變形，降低番茄的機(jī)械損傷，采取合適的抓取控制方式具有重要意義。綜上分析，應(yīng)選擇抓取控制方式3為番茄采摘機(jī)械手的最優(yōu)抓取控制方式，即變減速抓取控制策略。

3.3 抓取實(shí)驗(yàn)

為驗(yàn)證抓取過程中番茄塑性變形理論分析的正確性，以及選擇抓取控制方式3作為造成機(jī)械損傷較小的抓取控制方式的合理性和可行性。利用式(8)表征的番茄力-位移關(guān)系模型，結(jié)合當(dāng)前實(shí)驗(yàn)平臺(tái)，在抓取速度v0=1 mm/s，抓取時(shí)間t0=1 s的抓取條件下，對(duì)番茄樣本開展抓取實(shí)驗(yàn)與理論分析。

實(shí)驗(yàn)平臺(tái)采用德國(guó)SCHUNK公司生產(chǎn)的WSG 50型電動(dòng)二指平行夾指，如圖8所示。該末端執(zhí)行器整合了控制器、編碼器和SD存儲(chǔ)卡，與PC機(jī)的通訊協(xié)議為TCP/IP協(xié)議。其中一個(gè)力反饋手指可以實(shí)時(shí)測(cè)量番茄的接觸力，手指的位移由編碼器測(cè)量。設(shè)定機(jī)械手的抓取速度v0=1 mm/s，抓取時(shí)間t0=1 s，利用該實(shí)驗(yàn)平臺(tái)對(duì)上述選擇的紅熟前期的番茄以抓取控制方式3進(jìn)行抓取實(shí)驗(yàn)，分析其可行性。

圖8 抓取實(shí)驗(yàn)平臺(tái)Fig.8 Grasping experimental platform1.信息采集系統(tǒng) 2.WSG 50型機(jī)械手 3.力反饋手指 4.穩(wěn)壓電源

圖9為以抓取控制方式3完成穩(wěn)定抓取過程時(shí)，番茄的接觸力隨時(shí)間變化的實(shí)驗(yàn)數(shù)據(jù)與理論數(shù)據(jù)曲線。從圖9中可以看出，理論數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)的吻合度較高，通過Origin 8.0統(tǒng)計(jì)分析，確定實(shí)驗(yàn)數(shù)據(jù)和理論數(shù)據(jù)的決定系數(shù)R2=0.99，該模型能較準(zhǔn)確地描述抓取過程中番茄的受力特征。由式(13)可知，塑性變形是由抓取過程中接觸力變化最終決定，表明了以此為基礎(chǔ)進(jìn)行的番茄塑性變形理論分析的正確性。因此，可選擇抓取控制方式3為造成機(jī)械損傷較小的抓取控制方式，為優(yōu)化采摘機(jī)械手抓取過程提供依據(jù)。

圖9 抓取控制方式3的理論與實(shí)驗(yàn)曲線Fig.9 Theory and test curves of the third grasping control mode

4 結(jié)論

(1)以Burgers模型描述番茄的變形特征，通過Origin 8.0數(shù)據(jù)分析軟件對(duì)紅熟前期和紅熟中期的番茄蠕變曲線進(jìn)行擬合，得到粘彈性參數(shù)k1、k2、c1、c2，且擬合決定系數(shù)均在0.992以上。

(2)利用一次函數(shù)、巴特沃斯幅度平方函數(shù)以及指數(shù)函數(shù)，分別表征采摘機(jī)械手勻減速、先勻速后減速和變減速3種典型的抓取控制方式。以塑性變形作為衡量番茄機(jī)械損傷的關(guān)鍵因素，根據(jù)理論計(jì)算確定最小塑性變形量，選擇抓取控制方式3作為最優(yōu)抓取控制策略。當(dāng)抓取速度v0=1 mm/s，抓取時(shí)間t0=1 s時(shí)，產(chǎn)生的塑性變形Xp最小，紅熟前期和紅熟中期番茄的最小塑性變形分別為0.002 6、0.009 8 mm。

(3)實(shí)驗(yàn)與理論數(shù)據(jù)對(duì)比結(jié)果表明：實(shí)驗(yàn)和理論模型的吻合度較高，決定系數(shù)R2為0.99，驗(yàn)證了抓取控制方式3作為抓取策略的可行性和有效性。

1 BELFORTE G, DEBOLI R, GAY P, et al. Robot design and testing for greenhouse applications[J]. Biosystems Engineering, 2006, 95(3): 309-321.

2 方建軍. 移動(dòng)式采摘機(jī)器人研究現(xiàn)狀與進(jìn)展[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2004, 20(2): 273-278.

FANG Jianjun. Present situation and development of mobile harvesting robot[J]. Transactions of the CSAE, 2004, 20(2): 273-278.(in Chinese)

3 曹海燕, 盧立新. 梨果實(shí)跌落沖擊特性與流變模型[J]. 食品與生物技術(shù)學(xué)報(bào), 2007, 26(4):11-15.

CAO Haiyan, LU Lixin. The dropping impact behavior and rheological model of pear[J]. Journal of Food Science and Biotechnology, 2007, 26(4):11-15.(in Chinese)

4 盧立新, 王志偉. 基于準(zhǔn)靜態(tài)壓縮的果實(shí)黏彈塑性模型[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2005, 21(12):30-33.

LU Linxin, WANG Zhiwei. Nonlinear viscoelastic plastic model of the fruits under quasi-static compressive loading[J].Transactions of the CSAE, 2005, 21(12):30-33. (in Chinese)

5 TSAI C H D, KAO I, SHIBATA A, et al. Experimental study of creep response of viscoelastic contact interface under force control[C]∥International Conference on Intelligent Robots and Systems, IEEE, 2010:4275-4280.

6 LI Z G, LIU J Z, LI P P, et al. Study on the collision-mechanical properties of tomatoes gripped by harvesting robot fingers[J]. African Journal of Biotechnology, 2009, 8(24):7000-7007.

7 AHMADI E, BARIKLOO H, KASHFI M. Viscoelastic finite element analysis of the dynamic behavior of apple under impact loading with regard to its different layers[J]. Computers and Electronics in Agriculture, 2016, 121: 1-11.

8 HIGASHIMORI M, YOSHIMOTO K, KANEKO M. Active shaping of an unknown rheological object based on deformation decomposition into elasticity and plasticity[C]∥Robotics and Automation (ICRA), IEEE International Conference on Robotics and Automation, IEEE, 2010: 5120-5126.

9 姬偉, 李俊樂, 楊俊, 等. 機(jī)器手采摘蘋果抓取損傷機(jī)理有限元分析及驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2015,31(5): 17-22.

JI Wei, LI Junle, YANG Jun, et al. Analysis and validation for mechanical damage of apple by gripper in harvesting robot based on finite element method[J]. Transactions of the CSAE, 2015,31(5): 17-22. (in Chinese)

10 陳燕, 談建豪, 蔣志林,等. 柑橘機(jī)器人夾持損傷有限元預(yù)測(cè)及試驗(yàn)驗(yàn)證[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào),2016,37(5):98-102.

CHEN Yan, TAN Jianhao, JIANG Zhilin, et al. Finite element prediction and experimental verification for damage on citrus fruit from robot clamping[J]. Journal of South China Agricultural University, 2016, 37(5):98-102. (in Chinese)

11 劉繼展，白欣欣，李萍萍，等. 果實(shí)快速夾持復(fù)合碰撞模型研究[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2014, 45(4):49-54. http:∥www.j-csam.org/ch/reader/view_abstract.aspx?file_no=20140408amp;flag=1amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2014.04.008.

LIU Jizhan, BAI Xinxin, LI Pingping, et al.Complex collision model in High-speed gripping of fruit[J/OL]. Transactions of the Chinese Society for Agricultural Machinery,2014, 45(4):49-54. (in Chinese)

12 馬曉麗, 陳曉英, 閆語(yǔ)絲,等. 楊梅果的機(jī)械損傷試驗(yàn)和生物力學(xué)特性[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012, 28(16):282-287.

MA Xiaoli, CHEN Xiaoying, YAN Yusi, et al. Mechanical damage test and biomechanical characteristics of red bayberry fruit[J].Transactions of the CSAE, 2012, 28(16): 282-287. (in Chinese)

13 熊威，甘忠，許旭東，等. 基于粘彈塑性模型的時(shí)效成形仿真[J]. 塑性工程學(xué)報(bào), 2012, 19(2): 33-37.

XIONG Wei, GAN Zhong, XU Xudong, et al. Simulation for age forming based on viscoelastic plastic model[J]. Journal of Plasticity Engineering, 2012, 19(2): 33-37. (in Chinese)

14 張亞新，鄭曉，林國(guó)祥. 葵花籽壓榨過程中的塑性模型[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2009, 40(7): 138-142.

ZHANG Yaxin, ZHENG Xiao, LIN Guoxiang. Plasticity model for pressing of sunflower seeds[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(7): 138-142. (in Chinese)

15 TUDOR C S N, DANIELA I N, MARIAN V, et al. An experimental method for determining the time of relaxation at static compression of apples[J]. Journal of Engineering Studies amp; Research, 2011, 17(2):18-27.

16 LIU J Z, BAI X X, LI P P. Modified Burger’s model for describing creep behavior of tomato fruits[J]. Transactions of the Chinese Society of Agricultural Engineering, 2013, 29(9):249-255.

17 SAKAMOTO N, HIGASHIMORI M, TSUJI T, et al. An optimum design of robotic food handling by using Burger model[J]. Intelligent Service Robotics, 2008, 2(1):53-60.

18 李瀚如, 潘君拯. 農(nóng)業(yè)流變學(xué)導(dǎo)論[M]. 北京：農(nóng)業(yè)出版社, 1990.

19 楊玲, 張彩霞, 康國(guó)棟,等. 華紅蘋果果肉的流變特性及其主成分分析[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2015, 48(12):2417-2427.

YANG Ling, ZHANG Caixia, KANG Guodong, et al. Rheologic properties of Huahong apple pulp and their principal component analysis[J]. Scientia Agricultura Sinica, 2015, 48(12):2417-2427. (in Chinese)

20 李小昱, 朱俊平, 王為,等. 蘋果蠕變特性與靜載損傷機(jī)理的研究[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)：自然科學(xué)版, 1997，25(6):64-68.

LI Xiaoyu, ZHU Junping, WANG Wei, et al. A study on creep properties and static loading damage mechanisms of apple[J]. Journal of Northwest A amp; F University: Natural Science Edition, 1997,25(6):64-68. (in Chinese)

21 楊曉清, 王春光. 河套蜜瓜靜載蠕變特性的試驗(yàn)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2007, 23(3):202-207.

YANG Xiaoqing, WANG Chunguang. Creep properties of Hetao muskmelon under static loading[J]. Transactions of the CSAE, 2007, 23(3):202-207. (in Chinese)

22 李智國(guó). 番茄采摘機(jī)器人夾持系統(tǒng)的加減速過程研究[D]. 鎮(zhèn)江：江蘇大學(xué), 2009.

LI Zhiguo. Acc/Dec process of the grip system on tomato harvesting robot[D].Zhenjiang: Jiangsu University,2009. (in Chinese)

ImpactofRobotGraspingControlModesonMechanicalDamageofTomato

ZHOU Jun MENG Yimeng ZHANG Na YAN Zhenghong

(KeyLaboratoryofIntelligentAgriculturalEquipmentinJiangsuProvince,NanjingAgriculturalUniversity,Nanjing210031,China)

In order to analyze the impact of robot grasping control modes on mechanical damage of tomato, the Burgers model was used to express the deformation characteristics of tomato, and the viscoelastic parameters of tomato were obtained by the creep experiment. Three kinds of variations of speed, i.e., motion with uniform deceleration, deceleration after uniform motion and motion with variable deceleration were selected to indicate different grasping control modes. The linear function, Butterworth amplitude square function and exponential function were used to represent different grasping control modes during the robot grasping process. On this basis, the model of plastic deformation of tomato was solved, and the rules of tomato plastic deformations in different grasping control modes were analyzed when the grasping time was constant. Then, the optimal grasping strategy was selected based on the minimum plastic deformation. Experimental results indicated that the third grasping control mode was selected as the best grasping strategy, when the grasping speed was 1 mm/s, the grasping time was 1s, the plastic deformation of tomato was the smallest, the early maturity and middle maturity of plastic deformations were 0.002 6 mm and 0.009 8 mm, respectively. Using the third grasping control mode to carry out the grasping control experiment, the simulation results were highly identical with measured results, the determination coefficientR2was 0.99, which proved the rationality and feasibility of the third mode as the optimal grasping control mode. This research can offer a reference for robots grasping process optimization.

tomato; picking; robot; grasping control mode; mechanical damage

10.6041/j.issn.1000-1298.2017.11.003

TP24

A

1000-1298(2017)11-0021-07

2017-03-03

2017-03-23

國(guó)家自然科學(xué)基金項(xiàng)目(31471419)和教育部博士點(diǎn)基金博導(dǎo)類項(xiàng)目(20130097110043)

周俊(1974—)，男，教授，博士生導(dǎo)師，主要從事農(nóng)業(yè)機(jī)器人和機(jī)器視覺與模式識(shí)別研究，E-mail： zhoujun@njau.edu.cn