散鋆龍，楊會民，王學(xué)農(nóng)，郭文松，侯書林※

?

振動收獲過程中杏果實(shí)脫落的動態(tài)響應(yīng)

散鋆龍1，楊會民2，王學(xué)農(nóng)2，郭文松1，侯書林1※

（1. 中國農(nóng)業(yè)大學(xué)工學(xué)院，北京 100083；2. 新疆農(nóng)業(yè)科學(xué)院農(nóng)業(yè)機(jī)械化研究所，烏魯木齊 830091）

為研究杏果實(shí)脫落不同階段內(nèi)的瞬時響應(yīng)狀態(tài)與運(yùn)動狀態(tài)，該文建立了杏果實(shí)-果枝雙擺動力學(xué)模型，利用雙機(jī)位高速攝像儀記錄果實(shí)受迫振動并發(fā)生脫落全過程。理論分析與試驗(yàn)表明：杏果實(shí)-果枝雙擺系統(tǒng)在振動過程中主要表現(xiàn)為果枝擺動、果枝扭轉(zhuǎn)和果實(shí)傾擺3種運(yùn)動狀態(tài)，果實(shí)在空間運(yùn)動的軌跡接近橢圓形。使用PCC 3.1高速視頻控制分析軟件獲得杏果實(shí)在不同平面內(nèi)相對靜止原點(diǎn)的速度、脫落前瞬時速度和加速度數(shù)據(jù)。試驗(yàn)中不同平面內(nèi)果實(shí)振動脫落速度變化分為3個階段，通過MATLAB對不同平面內(nèi)果實(shí)相對于靜止原點(diǎn)的速度進(jìn)行傅里葉擬合，得到整體及各階段速度擬合曲線及相應(yīng)函數(shù)，同時推導(dǎo)出具有完整周期的前2個階段的速度函數(shù)周期與振幅間的關(guān)系。通過對比杏果實(shí)-果枝實(shí)際連接力，分析果實(shí)振動脫落特性，為優(yōu)化杏振動收獲機(jī)設(shè)計(jì)參數(shù)、提高果實(shí)采收率提供理論依據(jù)。

振動；采摘；模型；杏；高速攝像；運(yùn)動分析

0 引 言

利用振動方式對林果進(jìn)行機(jī)械化收獲作業(yè)是提高整個林果生產(chǎn)率的有效手段之一。自1980年以來，相對于人工采摘，使用機(jī)械進(jìn)行林果收獲具有明顯的經(jīng)濟(jì)優(yōu)勢，利用機(jī)械化收獲桃子效率可提高3倍[1]，采收咖啡豆效率可提高8倍[2]。國外已擁有多種針對不同林果品種的收獲機(jī)械，并對果樹振動收獲理論進(jìn)行了大量的研究。

目前林果機(jī)械化大面積收獲的研究主要集中在振動參數(shù)對收獲的影響[3-5]、果樹振動響應(yīng)狀態(tài)[6-7]和振動能量傳遞與耗散[8-9]等方面。Láng等[10-13]建立了單自由度果樹模型，對果樹受到不同慣性式振動激勵后的響應(yīng)狀態(tài)進(jìn)行了系統(tǒng)的建模分析。Castro-García 等[14-15]對不同頻率、振動時間和振動次數(shù)等振動參數(shù)對五針?biāo)伤勺訖C(jī)械化收獲的影響進(jìn)行了研究，保證果樹嫩枝和樹皮所受損傷最小的同時，獲得最佳的收獲效率。Du等[16-17]對櫻桃樹施加正弦振動激勵，研究主要分枝不同位置的振動響應(yīng)狀態(tài)、能量傳遞與損耗，分析果樹分枝節(jié)點(diǎn)角度與振動傳遞間的關(guān)系，并對山核桃果樹進(jìn)行3D建模并對振動特性進(jìn)行了分析[18-20]。

此外，對收獲前果實(shí)識別與定位的相關(guān)研究較多[21-25]，但分析果實(shí)振動脫落過程與運(yùn)動狀態(tài)的分析研究鮮有報道。王春耀等[26]利用高速攝像儀記錄二維平面內(nèi)果實(shí)受振動脫落過程，研究果實(shí)瞬時運(yùn)動變化規(guī)律，分析果實(shí)振動脫落不同階段運(yùn)動參數(shù)。杜小強(qiáng)等[27-28]借助扁球形電子果實(shí)，獲得果實(shí)受到三維激振收獲系統(tǒng)的振動響應(yīng)，研究分析不同階段果實(shí)所受沖擊加速度變化，為果實(shí)在振動收獲過程中損傷評估和振動機(jī)械參數(shù)優(yōu)化提供參考。

完善林果機(jī)械化收獲理論，研制高效林果收獲機(jī)械，是加快中國林果產(chǎn)業(yè)化發(fā)展，降低生產(chǎn)成本的重要途徑。本文首次采用雙機(jī)位高速攝像儀拍攝杏振動收獲時果實(shí)脫落全過程。研究三維空間內(nèi)杏果實(shí)脫落不同階段內(nèi)的瞬時響應(yīng)狀態(tài)與運(yùn)動狀態(tài)，測量分析果實(shí)各項(xiàng)運(yùn)動參數(shù)。建立杏果實(shí)-果枝動力學(xué)模型，結(jié)合振動收獲機(jī)振動特性，對比果實(shí)-果枝實(shí)際連接力，分析果實(shí)振動脫落特性，以期為優(yōu)化杏振動收獲機(jī)設(shè)計(jì)參數(shù)、提高果實(shí)收獲率提供理論依據(jù)。

1 杏果實(shí)-果枝動力學(xué)模型及求解

1.1 杏果實(shí)-果柄動力學(xué)模型

杏的實(shí)際生長情況相對特殊，果實(shí)果柄極短，果柄長度約等于花萼深度，且果實(shí)生長方向隨機(jī)，不同品種杏樹果實(shí)生長狀態(tài)不同，主要分為2種形式：果實(shí)直接生長在較粗的枝干上（圖1a），受迫振動時隨枝干運(yùn)動，生長在較粗枝干的果實(shí)分析時果實(shí)-果枝可簡化視為單擺模型；果實(shí)垂吊生長在單獨(dú)的細(xì)長枝干上（圖1b），此時細(xì)長果枝類似于蘋果長果柄，運(yùn)動過程中果枝也會發(fā)生相對擺動，垂吊生長的果實(shí)分析時果實(shí)-果枝可簡化為雙擺模型。根據(jù)庫麥提杏實(shí)際生長情況，本文僅對杏果實(shí)-果枝雙擺模型進(jìn)行分析和理論推導(dǎo)。

圖1 杏果實(shí)生長狀態(tài)

杏果實(shí)-果枝雙擺模型在受迫振動過程中主要運(yùn)動狀態(tài)有3種：果枝擺動、果實(shí)相對果枝間的扭轉(zhuǎn)和果實(shí)傾擺（圖2a）。若忽略果柄和空氣阻力對運(yùn)動的影響，可將果枝視為剛性桿，果實(shí)視為球體質(zhì)心為，以果枝與主枝干連接處為原點(diǎn)，為水平方向，為垂直方向，建立具有3個自由度的雙擺模型，如圖2b所示。

注：1：果枝帶動果實(shí)擺動；2：果枝帶動果實(shí)扭轉(zhuǎn)；3：果實(shí)擺動；θ為果枝與垂直方向夾角，（°）；φ為果實(shí)質(zhì)心與垂直方向夾角，（°）；δ為果實(shí)與果枝間的扭轉(zhuǎn)夾角，（°）；k0為果實(shí)與果枝連接處彈性系數(shù)，常數(shù)；k1為果枝彈性系數(shù)，常數(shù)；k2為果實(shí)與果枝間扭轉(zhuǎn)彈性系數(shù)，常數(shù)；l為果枝長度，mm；m1為果枝質(zhì)量，kg；r為果實(shí)半徑，mm；m0為果實(shí)質(zhì)量，kg。

杏果實(shí)-果枝雙擺系統(tǒng)由杏果實(shí)速率與轉(zhuǎn)動慣量產(chǎn)生的動能為

由重力和彈性撓度所產(chǎn)生的勢能為

式中為重力加速度，取9.8 m/s2；0為果實(shí)與果枝連接處彈性系數(shù)，常數(shù)；1為果枝彈性系數(shù)，常數(shù)；2為果實(shí)與果枝間扭轉(zhuǎn)彈性系數(shù)，常數(shù)。

由式（1）動能和式（2）勢能可得拉格朗日函數(shù)為

系統(tǒng)Rayleigh耗散函數(shù)為

可得保守系統(tǒng)拉格朗日方程

將式（3）帶入（5）可得杏果實(shí)-果枝雙擺系統(tǒng)振動動力學(xué)微分方程組為

1.2 杏果實(shí)-果枝雙擺動力學(xué)模型求解

針對杏果實(shí)-果枝雙擺動力學(xué)模型，杏樹在受到振動激勵后，果實(shí)振動規(guī)律與“果樹-收獲機(jī)”振動收獲系統(tǒng)保持一致，果枝在豎直方向上微幅振動，不考慮果實(shí)轉(zhuǎn)動慣量以及Rayleigh耗散，僅考慮果枝在水平方向的振動。

對杏果實(shí)-果枝雙擺系統(tǒng)非線性方程進(jìn)行分析，根據(jù)Kauderer逼近原理可得夾角逼近方程與、方向位移方程為

式中為果枝振幅，mm；為果實(shí)振幅，mm；為系統(tǒng)振幅，mm；為系統(tǒng)振動角速度，rad/s。

則有

由此可得到系統(tǒng)振動角速度與振型常數(shù)間的關(guān)系，進(jìn)而推導(dǎo)出系統(tǒng)頻率與振幅間的關(guān)系。

1.3 杏果實(shí)-果柄分離條件

隨著果實(shí)的形態(tài)學(xué)變化和“果實(shí)-果枝”系統(tǒng)質(zhì)量的增加，成熟果實(shí)與未成熟果實(shí)的振動響應(yīng)具有差異性，不同生長時期杏果實(shí)與果枝間的連接力不同，果實(shí)與果柄間連接力隨著果實(shí)成熟而變小[32]。當(dāng)機(jī)械施加的振動能量沿著果枝傳遞到果實(shí)處，被放大的加速度與果實(shí)自身重力形成慣性作用，促使果實(shí)產(chǎn)生脫落運(yùn)動。

由圖3二維平面杏果實(shí)受力分析可知：當(dāng)果實(shí)-果柄間連接力沿慣性力方向分力與果實(shí)自身重力的法向力之和小于慣性力F法向分力時，杏果實(shí)發(fā)生脫落，即：

注：F為果實(shí)與果柄連接力，N；F為果實(shí)慣性力，N；F為慣性力法向分力，N；F為慣性力切向分力，N；為果實(shí)加速度，m·s-2；a為法向加速度，m·s-2；a為果實(shí)加速度，m·s-2。

Note:Fis the connection force between fruit and carpopodium, N;Fis the inertia force, N;Fis the normal force of inertia force, N; Fis the tangential force of inertia force, N;is the acceleration, m·s-2;ais the normal acceleration, m·s-2;ais the tangential acceleration, m·s-2.

圖3 杏果實(shí)受力分析圖

Fig.3 Apricot fruit stress analysis

果實(shí)受振動激勵產(chǎn)生加速度與自重產(chǎn)生的慣性力F，可分解為果實(shí)分離的法向切應(yīng)力n和產(chǎn)生力矩的切向力F，果實(shí)法向加速度和切向加速度分別為a和a。

實(shí)際收獲過程中，果實(shí)分離主要依靠果實(shí)慣性力的法向分力F，果實(shí)重力沿連接點(diǎn)方向上的分力忽略不計(jì)，因此振動收獲過程中果實(shí)分離力必須滿足果實(shí)的法向切應(yīng)力大于果實(shí)與果枝沿慣性力法向方向連接力。即

2 杏果實(shí)振動脫落試驗(yàn)

為了驗(yàn)證理論分析的正確性，研究杏果實(shí)-果枝雙擺模型系統(tǒng)在振動采收過程中果實(shí)的運(yùn)動特性以及脫落特性，于2017年6月23日杏成熟期，在新疆巴音郭楞蒙古自治州輪臺縣新疆農(nóng)業(yè)科學(xué)院輪臺國家果樹資源圃（41°46′56″N， 84°13′22″E）中進(jìn)行杏果實(shí)振動脫落試驗(yàn)。試驗(yàn)果園水肥管理良好，無套種作物，果樹種植行距6 m，株距3 m，果樹為10 a樹齡的庫麥提杏樹。

2.1 儀器與設(shè)備

杏基本參數(shù)測量儀器：ATAGO 公司PAL-1型可溶性固形物含量測量儀，測試范圍0~53%；華馳HC10002型電子秤，精度0.01 g；上海昂軒儀器公司AM-50型數(shù)顯式推拉力計(jì)和自制夾頭測量果實(shí)-果柄分離力（如圖4），量程0~50 N，精度0.001 N；AY-4型果實(shí)硬度計(jì)測量果實(shí)表面硬度，量程0~50 kg/cm2，精度0.01 kg/cm2。

圖4 成熟期杏果實(shí)特性參數(shù)測量

視頻圖像設(shè)備：美國Vision Research公司Phantom v9.1型高速攝像儀2臺，最高分辨率為1 600×1 200，幀頻率范圍1 000～16 000 Hz； Nikon公司AF 50 mm f/1.8 D全畫幅單反鏡頭，最近對焦距離0.45 m，最大放大倍率0.151倍。

數(shù)據(jù)分析軟件：Phantom Camera Control (PCC) 3.1 高速攝像儀控制與分析軟件，MATLAB 2016數(shù)據(jù)分析軟件。

振動裝置：選用新疆農(nóng)業(yè)科學(xué)院農(nóng)業(yè)機(jī)械化研究所研制的4GZG-25型牽引式振動收獲機(jī)進(jìn)行果實(shí)振動脫落試驗(yàn)，如圖5所示。收獲機(jī)振動部件采用對稱雙偏心塊設(shè)計(jì)，工作時可穩(wěn)定輸出正弦振動激勵，振動頻率范圍10~15 Hz，振動周期范圍0.06~0.1 s，偏心塊轉(zhuǎn)速范圍500~1 000 r/min，一次作業(yè)可完成振動落果、收獲、除雜和裝箱作業(yè)，適用于大型果園收獲。

圖5 4GZG-25牽引式振動收獲機(jī)

2.2 杏果實(shí)特性參數(shù)測定

在試驗(yàn)果樹上不同位置隨機(jī)選擇5片區(qū)域，每個區(qū)域選擇20個果實(shí)，進(jìn)行杏果實(shí)主要特性參數(shù)測定，所有數(shù)據(jù)測量后取平均值。

測得成熟期庫麥提杏果實(shí)平均質(zhì)量為14.26 g；平均外形尺寸為長28.84 mm、寬27.93 mm、厚27.83 mm；成熟果實(shí)表面平均硬度為9.17 kg/cm2；成熟果實(shí)-果柄平均連接力2.923 N。

2.3 杏果實(shí)振動脫落參數(shù)采集

選擇型號相同的2臺高速攝像儀和鏡頭同時拍攝杏果實(shí)振動脫落過程，采集杏果實(shí)振動脫落過程圖像與數(shù)據(jù)。安裝調(diào)試好的高速攝像儀鏡頭中心線與果實(shí)在同一平面上（如圖6），距地面1.75 m，間距相等為1.5 m。以果樹根部中心為原點(diǎn)，水平地面為平面，垂直地面方向?yàn)榉较?，建立空間坐標(biāo)系。一號位與二號位攝像儀鏡頭分別垂直于平面和平面，果實(shí)位于兩攝像儀視野交點(diǎn)處，果實(shí)在2臺攝像機(jī)中方向數(shù)據(jù)相同。振動機(jī)施加振動激勵后，2臺攝像儀可拍攝不同平面內(nèi)果實(shí)從開始振動到脫落的全過程，聯(lián)合可得杏果實(shí)空間運(yùn)動狀態(tài)。

圖6 高速攝像儀機(jī)位安放示意圖

振動收獲前調(diào)整鏡頭焦距，使杏生長區(qū)域圖像在與高速攝像儀連接的電腦上可清晰顯示（圖7），對2臺攝像機(jī)均可見的杏果實(shí)進(jìn)行標(biāo)記。

選取拍攝區(qū)域無遮擋、振動脫落過程無碰撞、不發(fā)生干涉的單個杏果實(shí)進(jìn)行圖像采集與分析。攝像機(jī)拍攝采樣頻率1 000 fps，時間間隔周期1 000s，曝光度990s，圖像分辨率1 280×720。拍攝時間為果實(shí)靜止到完全脫落，圖像相對比例為0.348 mm/pix。

利用PCC 3.1高速視頻控制分析軟件，得到杏果實(shí)受迫振動發(fā)生脫落過程視頻圖像，并導(dǎo)出果實(shí)脫落不同階段在不同平面內(nèi)的瞬時位移坐標(biāo)、相對于果實(shí)靜止時目標(biāo)點(diǎn)的速度以及瞬時速度數(shù)據(jù)。利用MATLAB數(shù)據(jù)處理軟件對所有數(shù)據(jù)進(jìn)行分析，得到果實(shí)在空間上的運(yùn)動軌跡以及動態(tài)數(shù)據(jù)。

圖7 不同位置高速攝像儀拍攝的圖像

3 結(jié)果與分析

3.1 果實(shí)振動脫落運(yùn)動狀態(tài)及軌跡分析

利用PCC 3.1高速攝像儀視頻分析軟件，得到杏果實(shí)從受迫振動到發(fā)生脫落過程中主要的運(yùn)動狀態(tài)圖像（如圖8所示）。通過對比靜止時（圖8a）與發(fā)生受迫振動過程中（圖8b）果枝形態(tài)以及觀測點(diǎn)位置可知：受破振動過程果實(shí)發(fā)生了明顯的形態(tài)變化，果實(shí)-果枝整體位置發(fā)生明顯位移，即果實(shí)-果枝發(fā)生擺動，果枝發(fā)生了扭轉(zhuǎn)運(yùn)動；觀測點(diǎn)相對位置改變，標(biāo)記果實(shí)發(fā)生了明顯的偏轉(zhuǎn)。對比結(jié)果與理論分析中假設(shè)的果實(shí)運(yùn)動狀態(tài)一致。

圖8 杏果實(shí)-果枝運(yùn)動狀態(tài)

由分析軟件記錄的數(shù)據(jù)可知，目標(biāo)果實(shí)從開始振動到果實(shí)脫落的總時間為0.322 s。以目標(biāo)果實(shí)標(biāo)記中心為原點(diǎn)，每隔0.002 s對所跟蹤的杏果實(shí)在平面與平面上的位置坐標(biāo)進(jìn)行1次采集，得到目標(biāo)果實(shí)從開始振動到脫落過程中在2個平面上的軌跡坐標(biāo)數(shù)據(jù)，整合后繪制杏果實(shí)三維空間運(yùn)動軌跡曲線。

由圖9杏果實(shí)空間運(yùn)動軌跡可知：杏果實(shí)從振動到脫落前，整個空間運(yùn)動軌跡近似橢圓形。結(jié)果與Cooke和Rand[29-31]推導(dǎo)得出的果實(shí)-果柄動力運(yùn)動學(xué)軌跡的結(jié)論一致。

3.2 果實(shí)脫落運(yùn)動參數(shù)分析

以果實(shí)靜止時標(biāo)記中心為目標(biāo)原點(diǎn)，研究果實(shí)振動脫落過程中的運(yùn)動狀態(tài)，分析在不同平面上果實(shí)振動脫落過程中每0.002 s時間間隔內(nèi)相對于原點(diǎn)的速度的變化規(guī)律，即受迫振動后杏果實(shí)相對原點(diǎn)速度（）與時間的關(guān)系。同時，利用PCC 3.1軟件測量果實(shí)脫落時瞬時速度，獲得果實(shí)脫落前加速度，進(jìn)而推導(dǎo)出果實(shí)在不同平面上的受力。

圖9 杏果實(shí)空間運(yùn)動軌跡

3.2.1 杏果實(shí)在平面內(nèi)運(yùn)動狀態(tài)

利用MATLAB對目標(biāo)果實(shí)在平面速度與加速度數(shù)據(jù)進(jìn)行擬合分析可知，目標(biāo)果實(shí)受迫產(chǎn)生振動至脫落過程中相對果實(shí)靜止點(diǎn)速度符合傅里葉變換形式，相對果實(shí)靜止位置的速度擬合曲線如圖10所示。

圖10 杏果實(shí)在xoz平面上速度擬合曲線

杏果實(shí)在平面上速度的傅里葉擬合曲線2為0.970 4，其擬合函數(shù)方程為

式中為時間，s。

通過對高速攝像儀拍攝杏運(yùn)動圖像和數(shù)據(jù)分析，由平面上相對原點(diǎn)的速度與時間的傅里葉函數(shù)擬合曲線圖10可知，果實(shí)發(fā)生脫落時最大速度為0.112 m/s，果實(shí)在平面上從開始振動到脫落速度變化共經(jīng)歷的3個階段：0至0.148 s、0.148至0.248 s和0.248至0.322 s。果實(shí)在第3階段未達(dá)到完整周期就已脫落，為保證分析準(zhǔn)確性，僅對具有完整振動周期的第1階段和第2階段的數(shù)據(jù)進(jìn)行對比分析。

杏果實(shí)在平面上第1階段的速度與時間的傅里葉擬合函數(shù)曲線如圖11a所示，2為0.965 2，周期為0.15 s。擬合方程為

杏果實(shí)在平面上第2階段的速度與時間的傅里葉擬合函數(shù)曲線如圖11b所示，2為0.952 5，周期為0.15 s。擬合方程為

通過分析傅里葉擬合曲線圖11可知：果實(shí)速度在第1和第2階段內(nèi)速度均呈現(xiàn)周期性變化，速度變化周期為0.15 s；前2個階段速度函數(shù)振幅均為0.021 m/s。

圖11 杏果實(shí)在xoz平面上第1和第2階段速度擬合曲線

利用PCC 3.1軟件所得測量數(shù)據(jù)，杏果實(shí)在平面上，果實(shí)脫落前0.002 s與脫落時刻間相對速度為1.53 m/s，則此時觀測點(diǎn)的加速度為765 m/s2。根據(jù)牛頓第二定理，結(jié)合成熟果實(shí)平均質(zhì)量，計(jì)算可知，0.188 s時果實(shí)加速度產(chǎn)生的分離力達(dá)到成熟果實(shí)-果枝連接力；果實(shí)脫落時刻在平面上由加速度產(chǎn)生的分離力為11.22 N。

3.2.2 杏果實(shí)在平面內(nèi)運(yùn)動狀態(tài)

利用相同的方法對目標(biāo)果實(shí)在平面速度數(shù)據(jù)進(jìn)行擬合分析可知，相對果實(shí)靜止原點(diǎn)的速度與時間的傅里葉擬合曲線如圖12所示。

圖12 杏果實(shí)在yoz平面上速度擬合曲線

杏果實(shí)在平面上速度的傅里葉擬合曲線2為0.986，其函數(shù)方程為

由平面上相對原點(diǎn)速度與時間的傅里葉函數(shù)擬合曲線圖12可知，杏果實(shí)在平面上相對原點(diǎn)速度的變化也分為與在平面上相同的3個階段。果實(shí)發(fā)生脫落時最大速度為0.166 m/s，果實(shí)在第3階段未達(dá)到完整周期就已脫落，為保證分析準(zhǔn)確性，僅對第1和第2階段數(shù)據(jù)進(jìn)行對比分析。

對平面上運(yùn)動的前2個時間相對原點(diǎn)速度數(shù)據(jù)段進(jìn)行分析可知：

第1階段0至0.148 s，傅里葉擬合函數(shù)2為0.954 5，周期為0.13 s，振幅為0.027 m/s。擬合曲線如13a所示，擬合函數(shù)方程為

第2階段0.148至0.248 s，傅里葉擬合函數(shù)2為0.981 1，周期為0.13 s，振幅為0.054 m/s。擬合曲線如圖13b所示，擬合曲線方程為

圖13 杏果實(shí)在yoz平面上第1和第2階段速度擬合曲線

分析傅里葉擬合曲線圖13可知：在平面上杏果實(shí)相對原點(diǎn)速度在第1和第2階段內(nèi)均呈現(xiàn)周期性變化，速度變化周期為0.13 s；第1階段速度函數(shù)振幅為0.027 m/s，第2階段速度函數(shù)振幅為0.054 m/s。

利用PCC 3.1軟件測得數(shù)據(jù)測得在平面上，脫落前0.002 s與脫落時刻間瞬時速度為1.359 m/s，則此時觀測點(diǎn)的加速度為679.98 m/s2。由牛頓第二定理與成熟果實(shí)平均質(zhì)量計(jì)算可知：0.154 s時杏果實(shí)加速度產(chǎn)生的分離力達(dá)到成熟果實(shí)-果枝連接力；杏果實(shí)脫落時刻在平面上由加速度產(chǎn)生的分離力為9.7 N。

4 結(jié) 論

1）建立了三自由度杏果實(shí)-果枝雙擺模型，并對動力學(xué)模型進(jìn)行分析和求解，明確振動收獲中杏果實(shí)-果柄系統(tǒng)中各運(yùn)動參數(shù)間關(guān)系；對杏果實(shí)-果枝雙擺系統(tǒng)果實(shí)脫落條件進(jìn)行了理論分析。

2）利用雙機(jī)位高速攝像儀及分析軟件杏果實(shí)振動脫落過程進(jìn)行視頻圖像采集。目標(biāo)杏果實(shí)脫落過程時間為0.322 s，通過果實(shí)振動過程空間坐標(biāo)，繪制杏果實(shí)振動脫落空間運(yùn)動軌跡，分析可知杏果實(shí)-果枝系統(tǒng)在受迫振動后整個空間運(yùn)動軌跡近似橢圓形。

3）利用PCC 3.1分析軟件對與兩個平面內(nèi)同一目標(biāo)果實(shí)進(jìn)行運(yùn)動數(shù)據(jù)采集，可得果實(shí)脫落時相對靜止原點(diǎn)最大速度分量分別為0.112和0.166 m/s；脫落時刻瞬時加速度分別為765和679.98 m/s2；結(jié)合成熟杏果實(shí)重量計(jì)算得到加速度在各平面上產(chǎn)生的脫落力分量分別為11.22和9.7 N，均大于成熟杏果實(shí)與果柄間的平均連接力。

4）利用MATLAB對采集的速度數(shù)據(jù)進(jìn)行擬合分析，得到了目標(biāo)果實(shí)在與平面內(nèi)的傅里葉擬合曲線與擬合函數(shù)，其2值分別為0.970 4和0.986。

5）分析果實(shí)振動脫落速度擬合曲線可知，果實(shí)振動脫落可分為0~0.148 s、0.148~0.248 s和0.248 ~0.322 s三個階段。果實(shí)在第3階段未達(dá)到完整周期就已脫落，為保證分析的準(zhǔn)確性，對具有完整周期的第1與第2階段數(shù)據(jù)進(jìn)行分析可知：平面內(nèi)，前2個階段速度擬合函數(shù)2分別為0.965 2和0.952 5，速度變化周期均為0.15 s，函數(shù)振幅均為0.021 m/s；平面內(nèi)前2個階段速度擬合函數(shù)2分別為0.954 5和0.981 1，速度變化周期均為0.13 s，第2階段函數(shù)振幅是第1階段的2倍，分別為0.054和 0.027 m/s。杏果實(shí)-果枝雙擺系統(tǒng)振動周期接近激振源周期。

[1] Torregrosa A. Martin B, Garcia Brunton J, Bernad J J. Mechanical harvesting of processed peaches [J]. Transactions of the ASAE, 2008, 24(6): 723－729.

[2] Filho G R D S, Vieira H D, Rodrigues W P, et al. Comparison between manual and semi-mechanical harvest of coffee fruit in mountainous areas[J]. African Journal of Agricultural Research, 2015, 10(28): 2724－2730.

[3] Torregrosa A, Chaparro O, Ortiz C, et al. Design, construction and testing of an apricot tractor-trailed harvester[J]. Spanish Journal of Agricultural Research, 2008, 3(3): 333－340.

[4] Aristizábal I D, Oliveros C E, Alvarez F. Mechanical harvest of coffee applying circular and multidirectional vibrations[J]. Transactions of the Asae, 2003, 46(2): 205－210.

[5] Erdo?an D, Güner M, Dursun E, et al. Mechanical harvesting of apricots[J]. Biosystems Engineering, 2003, 85(1): 19－28.

[6] Whitney J D, Churchill D B, Hedden S L, et al. Trunk shakers of citrus harvesting. Part I: Measured trunk shaker and tree trunk motion [J]. Applied Engineering in Agriculture, 1988, 4(2): 93－101

[7] Whitney J D. Smerage G H. Block W A. Dynamic analysis of a trunk shaker-post system [J]. Transactions of the ASAE, 1990, 33(4): 1066－1070

[8] Horvath E, Sitkei G. Damping properties of plum trees shaken at their trunks [J]. Transactions of the ASAE, 2005,48(1): 19－25.

[9] Horvath E, Sitkei G. Energy consumption of selected tree shakers under different operational conditions [J]. Journal of Agricultural Engineering Research, 2001, 80(2): 191－199.

[10] Láng Z. One degree of freedom damped fruit tree model[J]. Agricultural Engineering International Cigr Journal, 2008, 51(3): 823－829.

[11] Láng Z. Dynamic modelling structure of a fruit tree for inertial shaker system design[J]. Biosystems Engineering, 2006, 93(1): 35－44.

[12] Láng Z. A fruit tree stability model for static and dynamic loading[J]. Biosystems Engineering, 2003, 85(4): 461－466.

[13] Láng Z, Csorba L. Modelling and comparing two canopy shapes using FEM[J]. International Journal of Horticultural Science, 2003.

[14] Castro-García S, Blanco-Roldán G L, Gil-Ribes J A. Vibrational and operational parameters in mechanical cone harvesting of stone pine (Pinus pinea L.)[J]. Biosystems Engineering, 2012, 112(4): 352－358.

[15] Castro-García S, Blanco-Roldán G L, Gil-Ribes J A. Vibrational and operational parameters in mechanical cone harvesting of stone pine (Pinus pinea L.)[J]. Biosystems Engineering, 2012, 112(4): 352－358.

[16] Du X,Chen D,Zhang Q, et al. Dynamic responses of sweet cherry trees under vibratory excitations[J]. Biosystems Engineering, 2012, 111(3): 305－314.

[17] Du X, Chen D, Zhang Q, et al Dynamic response of sweet cherry tree to oscillating excitations for mechanical harvest[J]. IFAC Proceedings Volumes, 2010, 43(26): 83-88.

[18] He L,Du X,Qiu G,Wu C.3D Reconstruction of Chinese Hickory Trees for Mechanical Harvest[C]//American Society of Agricultural and Biological Engineers Annual International Meeting-ASABE 2012.

[19] Wu C, He L, Du X, et al. 3D reconstruction of Chinese hickory tree for dynamics analysis[J]. Biosystems Engineering, 2014, 119(1): 69－79.

[20] Du X Q, Wu C Y, He L Y, et al. Dynamic characteristics of dwarf Chinese hickory trees under impact excitations for mechanical fruit harvesting[J]. International Journal of Agricultural & Biological Engineering, 2015, 8(1): 17－25.

[21] Thendral R, Suhasini A, Senthil N. A comparative analysis of edge and color based segmentation for orange fruit recognition[C]. International Conference on Communications and Signal Processing. IEEE, 2014: 463－466.

[22] Li D, Zhang L, Liu Y. A Study on Recognition Method of Fruits Based on Machine Vision[J]. Journal of Software Engineering, 2010, 9(4): 895－902.

[23] Villibor G P, Santos F L, Queiroz D M D, et al. Determination of modal properties, of the coffee fruit-stem system, using high speed digital video and digital image processing[J]. Acta Scientiarum Technology, 2016, 38(1): 41.

[24] 呂繼東，趙德安，姬偉，等采摘機(jī)器人振蕩果實(shí)動態(tài)識別[J].農(nóng)業(yè)機(jī)械學(xué)報，2012，43(05)：173－178+196. Lü Jidong, Zhao De’an, Ji wei, et al. Dynamic recognition of oscillating fruit for harvesting robot[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(5): 173－178+196.

[25] 呂繼東，趙德安，姬偉. 采摘機(jī)器人振蕩果實(shí)匹配動態(tài)識別[J]. 農(nóng)業(yè)工程學(xué)報，2013，29(20)：32－39. Lü Jidong, Zhao De’an, Ji Wei. Research on matching recognition method of oscillating fruit for apple harvesting robot[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE) 2013, 29(20): 32－39. (in Chinese with English abstract)

[26] 蔡菲，王春耀，王學(xué)農(nóng)，等. 基于高速攝像技術(shù)的振動落果慣性力研究[J].西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版)，2013，41(04)：208－212. Cai Fei, Wang Chunyao, Wang Xuenong, et al. Inertia force of fruits abscised by vibration based on high-speed video camera technology[J]. Journal of Northwest A&F University (Nat. Sci. Ed), 2013,41(04): 208－212. (in Chinese with English abstract)

[27] 杜小強(qiáng)，李松濤，賀磊盈，等. 三維激振果品采收機(jī)構(gòu)優(yōu)化設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(16)：48－55. Du Xiaoqiang, Li Songtao, He Leiying, et al. Optimal design and experiment on vibratory fruit harvesting mechanism with three-dimensional excitation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(16): 48－55. (in Chinese with English abstract)

[28] 杜小強(qiáng)，李黨偉，賀磊盈，等. 基于電子果實(shí)技術(shù)的機(jī)械振動采收過程果實(shí)運(yùn)動分析[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(17)：58－64. Du Xiaoqiang, Li Dangwei, He Leiying, et al. Fruit motion analysis in process of mechanical vibration harvesting based on electronic fruit technique[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(17): 58－64. (in Chinese with English abstract)

[29] Crooke J R, Rand R H. Vibratory fruit harvesting: A linear theory of fruit-stem dynamics[J]. Journal of Agricultural Engineering Research, 1969, 14(3): 201－200.

[30] Rand R H, Cooke J R. Vibratory fruit harvesting: A non-linear theory of fruit-stem dynamics [J]. Journal of Agricultural Engineering Research, 1970, 15(4): 357－356.

[31] Parchomchuk P, Cooke J R. Vibratory harvesting: an experimental analysis of fruit-stem dynamics[J]. Theoretical Computer Science, 1972, 40(12): 131－148.

[32] 散鋆龍，劉旋峰，牛長河，等. 杏果實(shí)成熟度特性參數(shù)與果柄分離力的相關(guān)性分析[J]. 農(nóng)業(yè)工程學(xué)報，2013，29(23)：62－68. San Yunlong, Liu Xuanfeng, Niu Changhe, et al. Relativity analysis between characteristic parameters of apricot’s ripeness and its fruit removal force[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(23): 62－68. (in Chinese with English abstract)

Dynamic response analysis of apricot fruit dropping during vibration harvesting

San Yunlong1, Yang Huimin2, Wang Xuenong2, Guo Wensong1, Hou Shulin1※

(1.100083; 2.830091)

In order to study the transient response state and movement state of apricot fruit during different stages of exfoliation, a dynamic model of apricot fruit-fruit branch double pendulum hanging on a single long stem branch according to the actual growth of apricot fruit was established in this paper. The apricot fruit-stem separation conditions were analyzed theoretically, and the physical parameters of Kumai apricot fruit at ripening stage were measured. The whole forced vibration process of the target fruitS was recorded by using 2 high-speed cameras with same specifications perpendicular to theplane and theplane, and the relevant datas such as shedding time, speed, and acceleration were obtained. Theoretical analysis and experiments showed that the apricot fruit-fruit branch double pendulum system mainly shows 3 kinds of motion states: fruit branch drives fruit to swing, fruit branch drives fruit to twist and fruit swing. The average weight of the ripe fruit of Kumari apricot was 14.26 g, the average external dimensions were 28.84 mm in length, 27.93 mm in width, and 27.83 mm in thickness; the average surface hardness of mature fruit was 9.17 kg/cm2; and the average connection force of ripe fruit-stem is 2.923 N. Based on the measured spatial coordinates of the fruit in theandplanes and Phantom Camera Control 3.1 high-speed video control software analysis, the whole spatial motion track of apricot fruit is nearly elliptical from vibration to shedding, the total time from vibration to shedding of target fruit was 0.322 s, and the maximum speed of the relative static origin of apricot fruit was 0.112 and 0.166 m/s when be dropped in spatial planesand. The instantaneous speed of the apricot fruit before shedding was 1.53 and 1.359 m/s, respectively, and the acceleration was 765 and 679.98 m/s2, respectively. Using MATLAB software, the velocity data of fruit in different planesandwere compared with the Fourier fitting curves of velocity and time relative to stationary origin, and the fitting curve of the speed of apricot fruit during vibration shedding and the corresponding function were obtained.2value of the fitting function were 0.970 4 and 0.986. Through analyzing the speed function curve, the process of fruit vibration shedding could be divided into 3 periods: 0 to 0.148 s, 0.148 to 0.248 s, and 0.248 to 0.322 s. In the third stage, the fruit has already fallen off before it reaches the full cycle. To ensure the accuracy of the analysis, the analysis of the velocity function in the first 2 stages with complete vibration period showed that: in theplane, the velocity fitting function2of the first 2 stages are 0.965 2 and 0.952 5, respectively, the speed change period is 0.15 s and the function amplitude is 0.021 m/s; the velocity fitting function2in the first 2 stages of theplane is 0.954 5 and 0.981 1, respectively, and the speed change period is 0.13 s, the amplitude of the second stage function is twice that of the first stage, which is 0.054 and 0.027 m/s, respectively. The vibration cycle of the apricot fruit-fruit branch double pendulum system is close to the period of the excitation cycle. According to the calculation of the fruit weight of mature apricot, the shedding force components produced by acceleration on each plane were 11.22 and 9.7 N respectively, both of which were greater than the average connected force between mature apricot fruit and stalk. Based on the above characteristics of fruit vibratory shedding, this study can provide a theoretical basis for clarifying the shedding mechanism during fruit harvesting, optimizing the design parameters of apricot vibrating harvester, and improving the fruit recovery rate.

vibration; picking; models; apricot; high-speed camera; motion analysis

10.11975/j.issn.1002-6819.2018.18.009

S225.93

A

1002-6819(2018)-18-0068-08

2018-05-22

2018-07-27

國家自然科學(xué)基金項(xiàng)目—新疆杏振動脫落動力學(xué)特性及振動采收機(jī)理研究（51465059）

散鋆龍，男，博士生，從事機(jī)械制造及自動化研究。 Email：sanyunlong@hotmail.com

侯書林，教授，博士生導(dǎo)師，主要從事機(jī)械設(shè)計(jì)及制造工藝研究。Email：hsl010@126.com

散鋆龍，楊會民，王學(xué)農(nóng)，郭文松，侯書林. 振動收獲過程中杏果實(shí)脫落的動態(tài)響應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(18)：68－75. doi：10.11975/j.issn.1002-6819.2018.18.009 http://www.tcsae.org

San Yunlong, Yang Huimin, Wang Xuenong, Guo Wensong, Hou Shulin. Dynamic response analysis of apricot fruit dropping during vibration harvesting [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 68－75. (in Chinese with English abstract) doi: 10.11975/j.issn.1002-6819.2018.18.009 http://www.tcsae.org