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

?

振動收獲過程中杏果實脫落的動態(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è)機械化研究所，烏魯木齊 830091）

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

振動；采摘；模型；杏；高速攝像；運動分析

0 引 言

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

則有

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

1.3 杏果實-果柄分離條件

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

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

注：F為果實與果柄連接力，N；F為果實慣性力，N；F為慣性力法向分力，N；F為慣性力切向分力，N；為果實加速度，m·s-2；a為法向加速度，m·s-2；a為果實加速度，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 杏果實受力分析圖

Fig.3 Apricot fruit stress analysis

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

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

2 杏果實振動脫落試驗

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

2.1 儀器與設(shè)備

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

圖4 成熟期杏果實特性參數(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è)機械化研究所研制的4GZG-25型牽引式振動收獲機進行果實振動脫落試驗，如圖5所示。收獲機振動部件采用對稱雙偏心塊設(shè)計，工作時可穩(wěn)定輸出正弦振動激勵，振動頻率范圍10~15 Hz，振動周期范圍0.06~0.1 s，偏心塊轉(zhuǎn)速范圍500~1 000 r/min，一次作業(yè)可完成振動落果、收獲、除雜和裝箱作業(yè)，適用于大型果園收獲。

圖5 4GZG-25牽引式振動收獲機

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

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

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

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

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

圖6 高速攝像儀機位安放示意圖

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

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

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

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

3 結(jié)果與分析

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

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

圖8 杏果實-果枝運動狀態(tài)

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

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

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

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

圖9 杏果實空間運動軌跡

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

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

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

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

式中為時間，s。

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

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

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

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

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

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

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

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

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

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

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

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

第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 杏果實在yoz平面上第1和第2階段速度擬合曲線

分析傅里葉擬合曲線圖13可知：在平面上杏果實相對原點速度在第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，則此時觀測點的加速度為679.98 m/s2。由牛頓第二定理與成熟果實平均質(zhì)量計算可知：0.154 s時杏果實加速度產(chǎn)生的分離力達到成熟果實-果枝連接力；杏果實脫落時刻在平面上由加速度產(chǎn)生的分離力為9.7 N。

4 結(jié) 論

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

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

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

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

5）分析果實振動脫落速度擬合曲線可知，果實振動脫落可分為0~0.148 s、0.148~0.248 s和0.248 ~0.322 s三個階段。果實在第3階段未達到完整周期就已脫落，為保證分析的準(zhǔn)確性，對具有完整周期的第1與第2階段數(shù)據(jù)進行分析可知：平面內(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。杏果實-果枝雙擺系統(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] 呂繼東，趙德安，姬偉，等采摘機器人振蕩果實動態(tài)識別[J].農(nóng)業(yè)機械學(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] 呂繼東，趙德安，姬偉. 采摘機器人振蕩果實匹配動態(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] 杜小強，李松濤，賀磊盈，等. 三維激振果品采收機構(gòu)優(yōu)化設(shè)計與試驗[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] 杜小強，李黨偉，賀磊盈，等. 基于電子果實技術(shù)的機械振動采收過程果實運動分析[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ù)與果柄分離力的相關(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é)基金項目—新疆杏振動脫落動力學(xué)特性及振動采收機理研究（51465059）

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

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

散鋆龍，楊會民，王學(xué)農(nóng)，郭文松，侯書林. 振動收獲過程中杏果實脫落的動態(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