楊會民，散鋆龍，陳毅飛，王學(xué)農(nóng)，牛長河，侯書林

?

不同振動特性參數(shù)對杏樹振動響應(yīng)的影響

楊會民1，散鋆龍2，陳毅飛1，王學(xué)農(nóng)1※，牛長河1，侯書林2

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

為研究振動時間、振動頻率和振動激勵點(diǎn)振幅等不同振動特性參數(shù)對杏樹振動的影響，該文利用ANSYS軟件對杏樹進(jìn)行了有限元建模分析；通過三因素三水平試驗(yàn)分析不同振動特性參數(shù)對杏樹振動檢測點(diǎn)的影響，利用Design Expert軟件進(jìn)行優(yōu)化分析，并進(jìn)行實(shí)驗(yàn)室驗(yàn)證試驗(yàn)。杏樹自由模態(tài)振動響應(yīng)分析表明最佳杏樹振動采收響應(yīng)頻率范圍為0～20 Hz；諧振動響應(yīng)分析可知在最佳頻率范圍內(nèi)，杏樹振動激勵點(diǎn)振幅為5、10和15 mm時，同一頻率下，隨著激勵振幅的增大，相同位置加速度增大，但振動曲線整體變化規(guī)律和趨勢一致。試驗(yàn)分析可知，各因素影響檢測點(diǎn)1和2加速度的強(qiáng)弱順序一致：激勵點(diǎn)振幅>振動頻率>振動時間；各因素影響檢測點(diǎn)3加速度的強(qiáng)弱順序?yàn)椋赫駝訒r間>激勵點(diǎn)振幅>振動頻率；建立3個檢測點(diǎn)的響應(yīng)方程，由下至上3個檢測點(diǎn)的回歸方程決定系數(shù)分別為0.906 7、0.879 3和0.973 3；多目標(biāo)參數(shù)優(yōu)化結(jié)果為：振動時間7.207 s，振動頻率15 Hz，激振點(diǎn)振幅10 mm，通過驗(yàn)證試驗(yàn)可知由下至上各檢測點(diǎn)加速度為10.4、10.2和9.3，與優(yōu)化值相近。該研究可為杏振動采收機(jī)械參數(shù)設(shè)計提供參考。

振動；收獲；機(jī)械化；杏；振動響應(yīng)試驗(yàn)

0 引 言

杏（Lam.）屬薔薇科（）杏屬（Mill.）落葉喬木，是世界重要的經(jīng)濟(jì)果樹之一[1]。中國杏樹栽培歷史悠久，根據(jù)遺傳學(xué)研究表明，新疆伊犁河北岸野生普通杏是世界栽培杏的起源種群[2]。截止2016年底，新疆維吾爾自治區(qū)杏樹種植面積12.487萬hm2，產(chǎn)量115.42萬t[3]。近年來新疆杏產(chǎn)業(yè)化發(fā)展迅速，但其收獲方式仍停留在利用長梯或長桿輔助的傳統(tǒng)人工采摘，作業(yè)風(fēng)險較高，同時增大了收獲成本，收獲環(huán)節(jié)成為制約新疆林果產(chǎn)業(yè)發(fā)展的首要問題。振動采收是目前林果采收的主要手段之一，利用機(jī)械化采收方式可有效降低勞動成本，提高林果采收效率。

國外已形成了一套完整的林果振動采收體系，林果振動采收理論完善，機(jī)具成熟，其中不同振動采收參數(shù)對林果收獲的影響是主要的研究方向之一[4-11]。由于果樹自身生長特性與材料特性，不同品種果樹的振動響應(yīng)不同，采收機(jī)械振動頻率、振幅、振動時間以及振動激勵點(diǎn)位置等振動特性參數(shù)的選擇和設(shè)定將直接影響林果采收效率與質(zhì)量[11-12]。早在1957年美國就已展開往復(fù)式振動采收機(jī)工作頻率和振幅對李子采收、果實(shí)損傷以及果樹搖動所需動力的研究[13-14]。Aristizábal等[15]研究使用圓周和多方向振動采收咖啡過程中，振動頻率、振幅、振動時間以及振動激勵位置對成熟果實(shí)采收率的影響。Polat等[16]為獲得巴丹木最佳振動采收參數(shù)，利用慣性式懸臂振動機(jī)進(jìn)行振動采收，研究相同振動時間下，不同振動頻率和振幅參數(shù)組合對果實(shí)脫落率的影響。

國內(nèi)開展了大量林果振動采收機(jī)構(gòu)設(shè)計、振動傳遞特性以及果樹響應(yīng)狀態(tài)等方面的基礎(chǔ)研究，但對振動特性參數(shù)對林果采收影響研究較淺[17-26]。范雷剛等[27]研究了海棠果樹受不同振動特性參數(shù)組合時，不同位置加速度變化。杜小強(qiáng)等[28]對所設(shè)計的三維激振林果采收機(jī)構(gòu)進(jìn)行優(yōu)化設(shè)計，并對系統(tǒng)振動響應(yīng)進(jìn)行了分析。王東等[29]利用有限元進(jìn)行果樹響應(yīng)分析，研究不同樹形果樹振動響應(yīng)狀態(tài)。

為研究不同振動采收參數(shù)對林果機(jī)械采收的影響，指導(dǎo)杏林果振動采收機(jī)械設(shè)計，優(yōu)化杏振動采收機(jī)械性能，本文利用ANSYS建立杏樹有限元模型，進(jìn)行杏樹自由模態(tài)與諧響應(yīng)仿真分析，確定杏振動采收參數(shù)范圍。在實(shí)驗(yàn)室內(nèi)進(jìn)行振動試驗(yàn)，分析振動試驗(yàn)結(jié)果，研究不同振動特性參數(shù)組合對振動采收的影響，得到杏最佳振動參數(shù)組合，并進(jìn)行實(shí)驗(yàn)室驗(yàn)證試驗(yàn)，以期為杏振動采收機(jī)械參數(shù)設(shè)計提供設(shè)計及理論依據(jù)。

1 杏振動采收機(jī)理

根據(jù)4GZG-25型牽引式偏心式林果采收樣機(jī)工作方式與振動特性，建立偏心塊對稱放置的振動采收機(jī)-果樹模型，如圖1所示[30-31]。在振動采收機(jī)與杏樹加持位置水平面上，以杏樹夾持位置中心為坐標(biāo)系原點(diǎn)，方向?yàn)橄到y(tǒng)水平運(yùn)動方向，方向?yàn)榕c水平垂直運(yùn)動方向。杏樹動力學(xué)特性由加持位置杏樹樹干等效彈性系數(shù)（常數(shù)）與等效阻尼系數(shù)（常數(shù)）表示。

注：m為偏心塊質(zhì)量，kg；ω為角速度，rad·s-1；t為時間，s；L1為振動箱與杏樹間距離，m； L2為偏心塊旋轉(zhuǎn)軸中心與x軸間距離，m；r為偏心塊轉(zhuǎn)動中心與質(zhì)心的距離，mm；k為杏樹樹干等效彈性系數(shù)，常數(shù)；c為杏樹樹干等效阻尼系數(shù)，常數(shù)。

對稱放置的2個偏心塊質(zhì)量相等，轉(zhuǎn)動方向相反，結(jié)合杏樹-振動機(jī)動力學(xué)模型，由牛頓第三定律可知，由于偏心塊做圓周運(yùn)動，采收機(jī)振動過程中偏心機(jī)構(gòu)整體在方向產(chǎn)生的力相互疊加，在方向上的力相互抵消，即方向受力為零。由此可的偏心塊產(chǎn)生的慣性力為

式中F為偏心機(jī)構(gòu)方向上產(chǎn)生的力，N；F為偏心機(jī)構(gòu)方向上產(chǎn)生的力，N；為偏心機(jī)構(gòu)在方向上的位移，mm。

對系統(tǒng)整體進(jìn)行受力分析，并簡化可得杏樹-振動采收機(jī)振動微分方程為

2 杏樹振動響應(yīng)試驗(yàn)材料與方法

由于所設(shè)計的4GZG-25型牽引式偏心式林果采收樣機(jī)正弦振動激勵輸出特性，為優(yōu)化樣機(jī)工作參數(shù)提供參考，研究杏樹受迫振動后的響應(yīng)狀態(tài)，獲得最佳的振動響應(yīng)參數(shù)；通過有限元建模與實(shí)驗(yàn)室振動響應(yīng)試驗(yàn)相結(jié)合對杏樹振動響應(yīng)進(jìn)行研究，對優(yōu)化后的參數(shù)進(jìn)行實(shí)驗(yàn)室驗(yàn)證試驗(yàn)。杏樹振動響應(yīng)試驗(yàn)路線如圖2所示。

圖2 杏樹振動響應(yīng)試驗(yàn)路線

根據(jù)牛頓第二定律，振動過程中杏果實(shí)在質(zhì)量一定的情況下，加速度越大，果實(shí)由于慣性產(chǎn)生的脫落力越大。因此本文以杏樹檢測點(diǎn)處加速度為主要評價指標(biāo)進(jìn)行試驗(yàn)，加速度單位為（為重力加速度，取9.8 m/s2）。

2.1 試驗(yàn)材料

以新疆主要杏品種庫麥提杏樹為主要研究對象，樹形結(jié)構(gòu)為自然開心型[29]。2018年6月底杏成熟期，將新疆農(nóng)業(yè)科學(xué)院輪臺國家果樹資源圃內(nèi)樹齡為3 a的果樹，移栽到新疆農(nóng)業(yè)科學(xué)院農(nóng)業(yè)機(jī)械化研究所物料特性實(shí)驗(yàn)室內(nèi)，進(jìn)行測繪與基本參數(shù)測量。

2.2 試驗(yàn)儀器與軟件

試驗(yàn)儀器：蘇州試驗(yàn)儀器設(shè)備公司DC-300振動試驗(yàn)臺，正弦推力2.94 kN，工作頻率5～5 000 Hz，最大位移40 mm；SA-3振動平臺功率放大器，額定輸出功率15 kVA；YMC壓電式IEPE型加速度傳感器，靈敏度50 mV/ms2，測試范圍±10，頻率范圍0～4 000 Hz。美國NI公司NI DAQ 9 178四通道信號采集器，NI 9 234振動信號采集卡。上海一恒科學(xué)儀器有限公司BPG-9 140A型精密鼓風(fēng)干燥箱，功率2 050 W，溫度范圍RT+10～25 ℃。美特斯工業(yè)系統(tǒng)（中國）有限公司CMT 6 103系列電子萬能試驗(yàn)機(jī)，最大應(yīng)力1 kN。

試驗(yàn)軟件：ANSYS Workbench 17.0有限元分析軟件；SolidWorks 2017三維建模軟件；Labview振動測試虛擬軟件；Design expert 10.0.3試驗(yàn)分析軟件。

2.3 試驗(yàn)方法

利用SolidWorks 2017軟件對已測繪的杏樹進(jìn)行三維建模，使用ANSYS Workbench 17.0軟件對杏樹進(jìn)行有限元分析，有限元模型如圖3所示。分析杏樹不附加任何振動激勵狀態(tài)下的自由模態(tài)，結(jié)合模態(tài)云圖仿真結(jié)果中杏樹響應(yīng)的一致性，獲得最佳的振動頻率范圍；研究杏樹不同振幅下受正弦振動激勵時沿樹干3個檢測位置的諧響應(yīng)狀態(tài)，每個檢測位置間距離為400 mm，以加速度大小表示檢測位置的響應(yīng)狀態(tài)。

圖3 杏樹有限元模型

根據(jù)有限元分析相應(yīng)結(jié)果，利用DC-300振動實(shí)驗(yàn)平臺對杏樹施加正弦振動激勵，研究振動時間、振動頻率和加持點(diǎn)振幅等振動采收主要因素對杏樹不同位置加速度的影響。杏樹從果園中砍伐后立即移入實(shí)驗(yàn)室，如圖4所示。去掉細(xì)小樹枝與樹葉，利用自制的連接器將杏樹與振動平臺連接，杏樹根部固定。沿杏樹樹干安裝加速度傳感器，傳感器安裝位置與有限元分析時果樹3個檢測點(diǎn)位置相同，通過Labview振動檢測軟件得到NI DAQ 9 178振動信號采集器所采集的加速度數(shù)據(jù)。采用Design Expert 10.0.3 軟件進(jìn)行試驗(yàn)方案設(shè)計，分析杏樹不同位置在不同振動激勵參數(shù)下的響應(yīng)狀態(tài)，得到最佳的振動激勵參數(shù)，并進(jìn)行實(shí)驗(yàn)室驗(yàn)證試驗(yàn)。

圖4 杏樹振動檢測平臺

3 杏樹有限元建模與分析

根據(jù)國家標(biāo)準(zhǔn)GB/T 1 931-2009木材含水率測定方法和GB/T 1 936.2-2009木材抗彎彈性模量測定方法測得庫麥提杏樹木材含水率為41.01%，彈性模量為425.35 MPa，密度為0.856 g/cm3，泊松比取0.330，由此可得杏樹木材剪切模量為159.9 MPa。將所測得的杏樹木材基本參數(shù)代入ANSYS Workbench建立的杏樹模型中，對杏樹進(jìn)行有限元分析。

3.1 杏樹自由振動模態(tài)分析

利用ANSYS Workbench中Modal模塊對杏樹進(jìn)行1～50階模態(tài)分析，得到反應(yīng)杏樹模型自由模態(tài)下階次、響應(yīng)頻率和杏樹整體形變關(guān)系的3D曲面圖，如圖5所示。

圖5 杏樹自由模態(tài)響應(yīng)分析

由圖5可知，杏樹在15、26、30和40階時具有典型自由模態(tài)響應(yīng)，對應(yīng)的頻率分別為5.5、10.5、15.1和29.3 Hz，果樹最大形變分別為286.5、261.9、267.2和273.5 mm。

分析杏樹典型自由模態(tài)云圖6可知，第15和26階典型自由模態(tài)中杏樹枝干末端響應(yīng)性較高；第30階果樹整體響應(yīng)一致性高，杏樹末端形變響應(yīng)最大；第40階29.3 Hz時果樹部分枝條變形嚴(yán)重，樹形整體結(jié)構(gòu)易被破壞。由此可知：杏樹振動采收過程中最佳頻率范圍可控制在5～15 Hz之間。

圖6 杏樹典型模態(tài)響應(yīng)云圖

3.2 杏樹諧振動響應(yīng)分析

通過對杏樹模型進(jìn)行諧響應(yīng)分析，確定杏樹枝干在固定幅值的正弦位移作用下的穩(wěn)定響應(yīng)，采用完全法進(jìn)行求解。根據(jù)果樹振動采收激振點(diǎn)振幅特性[31]，夾持點(diǎn)振幅設(shè)定為5、10和15 mm進(jìn)行杏樹諧振動響應(yīng)分析。由自由模態(tài)仿真結(jié)果設(shè)定頻率范圍為0～20 Hz，進(jìn)行諧振動分析可得到檢測點(diǎn)1、2、3位置加速度值變化曲線，如圖7。

圖7 不同振幅激勵下檢測點(diǎn)X軸方向加速度變化曲線

由圖7可知，激勵點(diǎn)振幅為5、10和15 mm時，在0～20 Hz振動頻率區(qū)間內(nèi)，檢測點(diǎn)1在4.67、7.33、15.33和18.66 Hz時具有典型加速度峰值；檢測點(diǎn)2在6.34和17.33 Hz時具有典型加速度峰值；檢測點(diǎn)3在7.01和17.3 Hz時具有典型加速度峰值。檢測點(diǎn)位置由下至上，振幅為5 mm時，最大加速度分別為7.52、6.57和4.93；振幅為10 mm時，最大加速度分別為14.89、12.96和9.87；振幅為15 mm時，最大加速度分別為22.36、19.71和14.8。

由此可知，在0～20 Hz內(nèi)，振幅為5、10和15 mm時，同一頻率下，隨著激勵振幅的增大，相同位置加速度增大，但振動曲線整體變化規(guī)律和趨勢一致。

4 杏樹實(shí)驗(yàn)室振動試驗(yàn)

為了驗(yàn)證有限元分析結(jié)果，獲得最優(yōu)振動采收參數(shù)組合，指導(dǎo)振動采收樣機(jī)參數(shù)設(shè)計，根據(jù)樣機(jī)工作原理與ANSYS有限元仿真結(jié)果，在2018年6月底在新疆農(nóng)業(yè)科學(xué)院農(nóng)業(yè)機(jī)械化研究所物料特性實(shí)驗(yàn)室中，按照試驗(yàn)路線圖2，利用振動試驗(yàn)平臺對剛移栽入實(shí)驗(yàn)室的成熟期杏樹施加振動激勵，進(jìn)行振動試驗(yàn)，試驗(yàn)數(shù)據(jù)由加速度傳感器與NI 9 234振動信號采集卡采集。

4.1 評價指標(biāo)

選擇振動時間（s）、激勵點(diǎn)振幅（mm）和振動頻率（Hz）3個試驗(yàn)因素，在實(shí)驗(yàn)室進(jìn)行振動試驗(yàn)，每組試驗(yàn)重復(fù)3次。樹枝加速度是果實(shí)振動脫落的主要評價指標(biāo)，本文選擇3個主要檢測點(diǎn)所測得的加速度作為評價指標(biāo)。

4.2 振動試驗(yàn)設(shè)計

進(jìn)行三因素三水平組合試驗(yàn)，根據(jù)杏樹有限元自由模態(tài)響應(yīng)與諧振動響應(yīng)分析結(jié)果，試驗(yàn)因素與水平編碼如表1所示。

利用Design Expert 10.0.3 軟件進(jìn)行試驗(yàn)方案設(shè)計，共進(jìn)行17組試驗(yàn)，每組重復(fù)3次，試驗(yàn)結(jié)果取平均值，試驗(yàn)數(shù)據(jù)如表2。

表1 試驗(yàn)因素和水平編碼

表2 試驗(yàn)設(shè)計方案及響應(yīng)值結(jié)果

注：1、2、3為圖3中3個檢測點(diǎn)的加速度。

Note:1,2, and3are the accelerations of the three detection points in Fig.3.

4.3 試驗(yàn)結(jié)果分析與回歸模型建立

通過Design Expert 10.0.3軟件對試驗(yàn)結(jié)果進(jìn)行回歸分析，并進(jìn)行多元回歸擬合，得到3個檢測點(diǎn)加速度1、2和3，3個試驗(yàn)指標(biāo)的回歸方程，并驗(yàn)證其顯著性，試驗(yàn)分析結(jié)果如表3所示。

4.3.1 檢測點(diǎn)1加速度的顯著性分析

由表3可知，整體模型極其顯著（<0.01），各因素影響檢測點(diǎn)1加速度的強(qiáng)弱順序?yàn)?；其中模型回歸項(xiàng)和顯著，因素不顯著（>0.05）可忽略，擬合后得到回歸方程決定系數(shù)2為0.906 7，擬合程度高，方程為

式中為激勵點(diǎn)振幅，mm；為振動頻率，Hz。

4.3.2 檢測點(diǎn)2加速度的顯著性分析

由表3可知，整體模型極其顯著（<0.01），各因素影響檢測點(diǎn)2加速度的強(qiáng)弱順序?yàn)?>，與檢測點(diǎn)1一致，其中模型回歸項(xiàng)和顯著，因素不顯著（>0.05）可忽略，擬合后得到回歸方程決定系數(shù)2為0.879 3，方程為

4.3.3檢測點(diǎn)3加速度的顯著性分析

由表3可知，整體模型極其顯著（<0.01），整體模型極其顯著，各因素影響檢測點(diǎn)3加速度的強(qiáng)弱順序?yàn)?>。其中除模型回歸項(xiàng)、與2的檢測點(diǎn)3加速度值大于0.05不顯著可忽略，其他回歸項(xiàng)均為顯著。由于檢測點(diǎn)3位于整個果樹末端，檢測點(diǎn)位置、振動過程中位移、果樹樹枝粗細(xì)以及果果枝形狀等原因會造成因素對不同檢測點(diǎn)的影響強(qiáng)弱順序不同?；貧w方程決定系數(shù)2為0.973 3，擬合程度極高，方程為

式中為振動時間，s。

表3 檢測點(diǎn)加速度回歸模型方差F檢驗(yàn)與分析

注：<0.01，為極其顯著；0.01<<0.05為顯著。

Note:<0.01, extremely significant; 0.01<<0.05, significant.

4.4 各項(xiàng)振動特性參數(shù)優(yōu)化

振動過程果實(shí)脫落的響應(yīng)加速度越大，果實(shí)所受脫落力越大，果實(shí)更易脫落，即3個檢測點(diǎn)加速度響應(yīng)均為最大值時滿足優(yōu)化條件。根據(jù)上述有限元分析結(jié)果、振動試驗(yàn)與項(xiàng)目組前期的試驗(yàn)結(jié)果[30-31]，為避免振動過程造成果樹損傷，3個檢測點(diǎn)加速度優(yōu)化范圍為0～20；根據(jù)檢測點(diǎn)由下至上所處的位置，設(shè)定3個檢測點(diǎn)值的權(quán)重，檢測點(diǎn)1、2設(shè)置權(quán)重為1，檢測點(diǎn)3接近果實(shí)生長位置權(quán)重為2。

根據(jù)上述參數(shù)設(shè)定，由Design Expert10.0.3軟件中Optimization模塊（響應(yīng)優(yōu)化）對檢測點(diǎn)進(jìn)行多目標(biāo)優(yōu)化。

目標(biāo)編碼變量優(yōu)化函數(shù)為

對目標(biāo)函數(shù)進(jìn)行優(yōu)化得到最優(yōu)振動采收參數(shù)組合為：振動時間7.207 s，振動頻率15 Hz，激振點(diǎn)振幅10 mm。

4.5 實(shí)驗(yàn)室驗(yàn)證試驗(yàn)

為驗(yàn)證最終優(yōu)化參數(shù)，按優(yōu)化后得到參數(shù)調(diào)整振動實(shí)驗(yàn)平臺設(shè)置參數(shù)，在實(shí)驗(yàn)室中進(jìn)行振動試驗(yàn)，試驗(yàn)重復(fù)3次取平均值，優(yōu)化值與實(shí)際測量值之間分析驗(yàn)證結(jié)果如表4所示，實(shí)驗(yàn)值與優(yōu)化參數(shù)相近，偏差較小。

表4 優(yōu)化值與實(shí)際值對比結(jié)果

5 結(jié) 論

1）對建立的杏樹模型進(jìn)行ANSYS有限元分析，通過杏樹自由模態(tài)振動響應(yīng)分析可知最佳杏樹振動采收響應(yīng)頻率范圍為0～20 Hz；通過諧振動響應(yīng)分析可知在0～20 Hz內(nèi)，杏樹振動激勵點(diǎn)振幅為5、10和15 mm時，同一頻率下，隨著激勵振幅的增大，相同位置加速度增大，但振動曲線整體變化規(guī)律和趨勢一致。

2）研究杏樹振動特性參數(shù)（振動時間、激勵點(diǎn)振幅和振動頻率）對不同檢測點(diǎn)振動加速度的影響規(guī)律，進(jìn)行三因素三水平試驗(yàn)分析，各參數(shù)對檢測點(diǎn)1和2的影響強(qiáng)弱順序一致為>>；對檢測點(diǎn)3的影響強(qiáng)弱順序?yàn)?>。

3）建立3個檢測點(diǎn)響應(yīng)方程，由下至上3個檢測點(diǎn)的回歸方程決定系數(shù)分別為0.906 7、0.879 3和0.973 3，為杏樹振動特性分析提供參考。

4）利用Design Expert10.0.3軟件響應(yīng)優(yōu)化模塊對檢測點(diǎn)響應(yīng)方程進(jìn)行優(yōu)化，獲得最優(yōu)振動采收參數(shù)組合為：振動時間7.207 s，振動頻率15 Hz，激振點(diǎn)振幅10 mm，經(jīng)試驗(yàn)驗(yàn)證可知由下至上3個檢測點(diǎn)加速度分別為10.4、10.2和9.3，試驗(yàn)值與優(yōu)化參數(shù)相符，該研究可為杏振動采收機(jī)械參數(shù)設(shè)計提供設(shè)計及理論依據(jù)。

[1] 包文泉，烏云塔娜，王淋，等. 野生杏和栽培杏的遺傳多樣性和遺傳結(jié)構(gòu)分析[J]. 植物遺傳資源學(xué)報，2017，18(2)：201－209.

Bao Wenquan, Wuyun Tana, Wang Lin, et al. Genetic diversity and population structure of the wild apricot and cultivation apricot[J]. Journal of Plant Genetic Resources, 2017, 18(2): 201－209. (in Chinese with English abstract)

[2] 廖康. 新疆野生果樹資源研究[M]. 烏魯木齊：新疆人民出版社，2013.

[3] 中華人民共和國國家統(tǒng)計局. 中國統(tǒng)計年鑒[M]. 北京：中國統(tǒng)計出版社，2017.

[4] BlancoRoldán G L, GilRibes J A, Kouraba K, et al. Effects of trunk shaker duration and repetitions on removal efficiency for the harvesting of oil olives[J]. Applied Engineering in Agriculture, 2009, 25(3): 329－334.

[5] Pezzi F, Caprara C. Mechanical grape harvesting: Investigation of the transmission of vibrations[J]. Biosystems Engineering, 2009, 103(3): 281－286.

[6] Torregrosa A, Ortí E, Martín B, et al. Mechanical harvesting of oranges and mandarins in Spain[J]. Biosystems Engineering, 2009, 104(1): 18－24.

[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] Láng Z. A one degree of freedom damped fruit tree model[J]. Agricultural Engineering International Cigr Journal, 2008, 51(3): 823－829.

[10] 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.

[11] Torregrosa A, Martin B, Garcia Brunton J, et al. Mechanical harvesting of processed peaches[J]. Applied Engineering in Agriculture, 2008, 24(6): 723－729.

[12] 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.

[13] Castro-García S, Blanco-Roldán G L, Gil-Ribes J A, et al. Dynamic analysis of olive trees in intensive orchards under forced vibration[J]. Trees, 2008, 22(6): 795－802.

[14] Adrian P A, Fridley R B. Mechanical fruit tree shaking[J]. California Agriculture, 1958, 12: 3－15.

[15] 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.

[16] Polat R, Guner M, Dursun E, et al. Mechanical harvesting of almond with an inertia type limb shaker[J]. Asian Journal of Plant Sciences, 2007, 6(3): 528－532.

[17] 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.

[18] 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.

[19] 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.

[20] 杜小強(qiáng)，倪柯楠，潘珂，等. 可調(diào)振幅單向拽振式林果采收機(jī)構(gòu)參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(16)：25－32.

Du Xiaoqiang, Ni Ke’nan, Pan Ke, et al. Parameter optimization of stroke-adjustable and monodirectional pulling fruit harvester[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE) , 2014, 30(16): 25－32. (in Chinese with English abstract)

[21] 杜小強(qiáng)，倪柯楠，武傳宇. 基于外旋輪線軌跡的果品振動采收機(jī)構(gòu)研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報，2016，47(3)：59－66.

Du Xiaoqiang, Ni Ke’nan, Wu Chuanyu. Vibratory harvesting mechanism for tree fruit based on epitrochoid[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3): 59－66. (in Chinese with English abstract).

[22] 杜小強(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)

[23] 喬園園，牛長河，孟詳金，等. 牽引式林果振動采收機(jī)的設(shè)計與田間試驗(yàn)[J]. 新疆農(nóng)業(yè)科學(xué)，2015，52(3)：528－534.

[24] 王長勤，許林云，周宏平，等. 偏心式林果振動采收機(jī)的研制與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(16)：10－16.

Wang Changqin, Xu Linyun, Zhou Hongping, et al. Development and experiment of eccentric-type vibratory harvester for forest-fruits[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(16): 10－16. (in Chinese with English abstract)

[25] 蔡菲，王春耀，王學(xué)農(nóng)，等. 基于高速攝像技術(shù)的振動落果慣性力研究[J]. 西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版，2013，41(4)：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(4): 208－212. (in Chinese with English abstract)

[26] 付威，何榮，曲金麗，等. 自走式矮化密植紅棗收獲機(jī)的設(shè)計[J]. 農(nóng)機(jī)化研究，2014，36(4)：106－109.

[27] 范雷剛，王春耀，劉夢霞，等. 振動參數(shù)對果樹采收影響的試驗(yàn)研究[J]. 農(nóng)機(jī)化研究，2016，38(10)：165－168.

[28] 杜小強(qiáng)，李松濤，賀磊盈，等. 三維激振果品采收機(jī)構(gòu)優(yōu)化設(shè)計與試驗(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)

[29] 王冬，陳度，王書茂，等. 基于有限元方法的整形果樹振動收獲機(jī)理分析[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(增刊1)：56－62.

Wang Dong, Chen Du, Wang Shumao, et al. Analysis on vibratory harvesting mechanism for trained fruit tree based on finite element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(Supp.1): 56－62. (in Chinese with English abstract)

[30] 散鋆龍，楊會民，王學(xué)農(nóng)，等. 振動方式和頻率對杏樹振動采收響應(yīng)的影響[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(8)：10－17.

San Yunlong, Yang Huimin, Wang Xuenong, et al. Effects of vibration mode and frequency on vibration harvesting of apricot trees[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(8): 10－17. (in Chinese with English abstract)

[31] 散鋆龍，楊會民，王學(xué)農(nóng)，等. 振動收獲過程中杏果實(shí)脫落的動態(tài)響應(yīng)分析[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(18)：68－75.

San Yunlong,Yang Huimin,Wang Xuenong, et al. 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)

Influence of different vibration characteristic parameters on vibration response of apricot trees

Yang Huimin1, San Yunlong2, Chen Yifei1, Wang Xuenong1※, Niu Changhe1, Hou Shulin2

(1.830091,; 2.100083)

In this paper, the combination of finite element modeling and experiment was used to study the effects of different vibration parameters such as vibration time (), amplitude () of the vibration excitation point on apricot tree vibration and vibration frequency (), and optimize the performance of forest fruit vibration harvesting prototype. The effects of different vibration parameters on the vibration detection points of apricot trees were analyzed by vibration response test. The response analysis of apricot tree free modal vibration using ANSYS Workbench software showed that apricot trees had typical free modal responses at 15, 26, 30 and 40 orders, the corresponding frequencies were 5.5, 10.5, 15.1, and 29.3 Hz, respectively; and the maximum displacement deformation of fruit trees was 286.5, 261.9, 267.2, and 273.5 mm, respectively. The typical free modal cloud image analysis showed that the typical free modalities of the 15th and 26th orders were higher in the terminal end responsiveness of the apricot branches. The 30th order fruit tree had a high overall response consistency, and the apricot tree had the largest deformation response at the end; at the 40th order of 29.3 Hz, some branches of the fruit trees were severely deformed, and the overall structure of the tree was easily destroyed. The optimal response frequency of apricot tree vibration harvesting ranged from 0 to 20 Hz. The harmonic response analysis showed that in the optimal frequency range, the acceleration of the same position increased with the increase of the excitation amplitude at the same frequency, but the overall variation of the vibration curve was consistent with the trend; When the amplitude of the vibration excitation point of the apricot tree was 5, 10 and 15 mm, at the same frequency, as the excitation amplitude increased, the acceleration at the same position increased, but the overall variation of the vibration curve was consistent with the trend. The three-factor and three-level vibration response tests were conducted to study the effects of vibration time, vibration frequency and excitation point amplitude on the acceleration of three different detection points. The multivariate regression analysis of variance showed that the accelerations1and2of detection points 1 and 2 were less than 0.000 1, the overall model was highly significant (<0.01), and the model regression termsandwere significant. The acceleration3of the detection point 3 was 0.000 1, and the overall model was highly significant (<0.01), except that the acceleration3value of the detection point 3 of the model regression terms,and2was not significant (>0.05), and other regression terms were significant (<0.05). The factors which affected the acceleration of the detection points 1 and 2 were the same as the amplitude of the excitation point, vibration frequency and vibration time. The order of the magnitude of the acceleration affecting the detection point 3 was the vibration time, the amplitude of the excitation point, and the vibration frequency. Through the response equations of the three detection points, the coefficient of determination2of the regression equations from the bottom to the top three detection points were 0.906 7, 0.879 3 and 0.973 3, respectively. Using the Design-Expert 10.0.3 software to optimize the detection point response equation, the optimal vibration recovery parameter combination was that the vibration time was 7.207 s, the vibration frequency was 15 Hz, and the amplitude of the excitation point was 10 mm. The verification test showed that the acceleration from the bottom to the top of each detection point was 10.4, 10.2 and 9.3, which was similar to the optimized value. These conclusions can provide design and theoretical basis for the design of mechanical parameters of apricot vibration harvesting.

vibrations; harvesting; mechanization; apricot; vibration response test

10.11975/j.issn.1002-6819.2019.02.002

S225.93

A

1002-6819(2019)-02-0010-07

2018-09-10

2019-01-09

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

楊會民，助理研究員，主要從事農(nóng)業(yè)工程研究。Email：yhm_shz@163.com

王學(xué)農(nóng)，研究員，主要從事農(nóng)業(yè)工程與裝備研究。Email：xjwxn2010@ sina.com

楊會民，散鋆龍，陳毅飛，王學(xué)農(nóng)，牛長河，侯書林. 不同振動特性參數(shù)對杏樹振動響應(yīng)的影響[J]. 農(nóng)業(yè)工程學(xué)報，2019，35(2)：10－16. doi：10.11975/j.issn.1002-6819.2019.02.002 http://www.tcsae.org

Yang Huimin, San Yunlong, Chen Yifei, Wang Xuenong, Niu Changhe, Hou Shulin. Influence of different vibration characteristic parameters on vibration response of apricot trees[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(2): 10－16. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.02.002 http://www.tcsae.org