基于ANSYS對振動果樹枝干“Y”型響應的研究

呂夢璐，王春耀，羅建清，魏庭鵬，劉子龍

(新疆大學 機械工程學院，烏魯木齊 830047)

?

基于ANSYS對振動果樹枝干“Y”型響應的研究

呂夢璐，王春耀，羅建清，魏庭鵬，劉子龍

(新疆大學 機械工程學院，烏魯木齊 830047)

以機械振動采收果樹為研究對象，利用ANSYS Workbench軟件對機械振動采收果樹枝干進行動力學研究。通過Pro/E 5.0參數化方法建立果樹枝干三維實體模型,將模型簡化處理后導入ANSYS Workbench中進行仿真分析[1]，對果樹枝干進行10階模態分析，得到果樹枝干的固有頻率和振型；對其進行諧響應分析，計算出果樹枝干產生共振的頻率、共振幅值和共振幅度最大的位置。進行室內沖擊試驗，找到共振頻率，與有限元分析結果進行比較，表明有限元分析法可以為撞擊式收獲機的優化和設計提供理論依據。

果樹模型；Pro/E；ANSYS Workbench；有限元分析；振動

0 引言

隨著我國農業林業的繁榮發展、果品產量和品質要求的提高，機械振動器采收果實得到廣泛應用。機械振動器的采收率取決于振動器的工作參數和果樹的生物特性。隨著計算機輔助技術在農業應用領域的不斷發展，在控制理論、相對理論、信息處理技術和計算機技術基礎上實現面向果樹的自適應機械振動收獲，成為當前應用最廣泛的技術之一[2]。鄭甲紅等[3]應用Pro/E建立了果樹模型，在激振力不變的情況下，調整夾持位置，通過對模型進行模態分析和諧響應分析，得到最優夾持位置，但有關果樹材料屬性賦予方面的問題沒有提及。

本文從果樹的生物特性角度出發，以果樹枝干為研究對象，首先，確定果樹材料屬性，如含水率、密度、彈性模量等；然后，使用Pro/E建立果樹的三維模型；最后，利用ANSYS Workbench從模態、諧響應分析兩方面分析果樹的動力學特性[4]。在課題組試驗的基礎上，設法構造既能反映樹枝干特征又能符合試驗要求的三維模型，以達到了解或設計系統的目的。

1 果樹枝干材料屬性

1.1 果樹枝干幾何特征

圖1顯示用游標卡尺測量果樹枝干不同部位得到的長軸和短軸及其之間的關系。其中，長軸與短軸分別取自一棵果樹從樹根到樹枝分叉處的14個不同橫截面。

圖1 果樹枝干橫截面長短軸相關性圖

從圖1中可以看出：果樹樹枝干截面的形狀并非理想的圓形，它的長軸與短軸呈現顯著的正比例關系，比例值約為1.028 5，R2=0.861 0[5]。

相關系數R是衡量兩個隨機變量之間線性相關程度的指標。通常|R|大于0.8時，認為兩個變量有很強的線性相關性。為了方便果樹枝干的三維建模和有限元分析，假設樹枝干的橫截面形狀為理想圓形，且半徑從根部到頂端呈線性下降趨勢[6]。

1.2 果樹枝干含水率及密度

測量果樹枝干含水率采用的是烘干法。取6段沒有夾皮、節疤、腐蝕、蟲蛀等缺陷的木材，烘烤前用電子秤稱其質量，分別為13.70、17.77、20.13、23.91、29.25、33.73g；經測試、計算得6組含水率分別為43.39%、43.21%、43.07%、41.86%、40.11%、38.92%。即平均含水率W=41.76%[7]，平均密度ρ=936.86kg/m3。

1.3 果樹枝干靜態彈性模量

本研究的實驗設備采用微機控制電子萬能實驗機，設備型號為 CMT5305。利用四點彎曲試驗法測得木材順紋抗彎彈性模量[8]，軟件界面獲取壓力—位移曲線和數值。其彈性模量計算公式為

(1)

其中，E為靜彈性模量(MPa)；F為載荷—變形圖中直線段載荷的增加量(kN)；wE為在E點橫截面形心在垂直于x軸方向的線位移(mm)；a為上壓頭與下壓頭之間的長度。

用四點彎曲試驗[9]測得負荷-變形曲線線性部分如圖2所示。此階段果樹枝干屬于彈性材料，符合廣義胡克定律，由公式(1)求得果樹枝干彈性模量如圖3所示，并可得果樹枝干線性階段彈性模量的平均值為E=1.330GPa。

圖2 四點彎曲果樹枝干負荷-變形

圖3 四點彎曲測果樹枝干彈性模量

2 建模與有限元分析

利用Pro/E建立果樹振動動力學的簡化模型，并導入到ANSYS Workbench中，進行理論研究，分析果樹振動特性[3]。

2.1 三維建模

將果樹按兩級建模即分為主干和樹枝，分枝模式屬于合軸分枝單歧聚傘花序，且每個分叉點上最多只連接2個子分枝[10]。樹干和樹枝的生長形態可近似于歐拉—伯努利的彈性梁，樹干和樹枝的幾何橫截面可以簡化為圓形，將樹簡化為Y字型梁[11]，其直徑和梁長呈線性變化。使用Pro/E掃描混合和倒圓角功能建立果樹模型。其中，分枝角度為35°。

對于樹干而言，由于土壤影響，樹根部錨固情況不同，可以簡化為剛性固定端；對樹枝而言，連接情況復雜，則認為是彈性固定端。設定材料屬性，密度ρ=936.86kg/m3，彈性模量E=1.330E+09Pa，泊松比ν=0.33。可算得剪切模量G=5E+08Pa，體積模量K=1.3039E+09Pa。自動劃分網格，單元尺寸為5mm[12]。

對果樹枝干三維模型進行以下兩方面的動力學分析[13]：模態分析、諧響應分析。

2.2 模態分析

樹根部位施加固定約束面，由ANSYS Workbench進行模態分析，得果樹枝干振動特性的固有頻率和振型參數。計算前10階變形分析云圖、應力分析云圖，前6階變形云圖如圖4所示，各階模態頻率如表1所示。

圖4 前6階變形云圖

模態頻率/Hz11．01121．01337．8048．90510．20610．70724．80824．90937．40

結果分析：第1階變形到第6階變形主要表現為以彎曲與扭轉為主的變形，最大振幅出現在樹枝分叉的兩端部。第1階(頻率為1.011Hz)變形和第2階(頻率為1.013Hz)變形云圖基本一致，第1階最大振幅為61.467mm，第2階最大振幅為37.532mm，第2階最大振幅比第1階最大振幅大0.077mm，波動范圍較小；第3階(頻率為7.752 4Hz)變形和第4階(頻率為8.899 8Hz)變形云圖基本一致，第3階最大振幅為61.959mm，第2階最大振幅為65.305mm，第4階最大振幅比第3階最大振幅大3.346mm，且第3階左側枝以彎曲變形為主，第4階左側枝以扭轉變形為主；第5階振幅相對第4階振幅有所減小；第5階(頻率為10.177Hz)變形和第6階(頻率為10.677Hz)變形云圖基本一致，第5階最大振幅為37.455mm，第6階最大振幅為73.806mm，第6階最大振幅比第5階最大振幅大12.339mm，波動很大，且第5階右側枝以彎曲變形為主，第6階右側枝以扭轉變形為主。當外界激勵和系統本身的固有頻率相等時，系統共振[14]。

2.3 諧響應分析

對模型進行諧響應分析即頻率響應分析，確定果樹枝干在已知頻率和幅值的正弦荷載作用下的穩定響應，采用模態疊加法求解。在距樹根固定端10cm截面處施加740 N的集中力(此力由試驗確定)，由本課題組試驗[15]得到頻率為5.371Hz，阻尼比為0.050 84。設置諧響應最大頻率應為43.924/1.5=29.283Hz[11]，ANSYS Workbench中輸入30Hz，即只考慮0～30Hz頻率范圍內的結構響應情況。設置阻尼參數，頻率輸5.371,阻尼比輸0.050 84。

通過計算可得果樹枝干三維模型在諧響應分析過程中變形分析云圖、變化曲線及應力分布云圖，結果如圖5所示。

圖5 三維模型的諧響應分析

結果分析：頻率在5～10Hz時，振幅快速遞增，從0.144 2mm到0.877 75mm；頻率在10～15Hz時，振幅快速遞減，從0.877 75mm到0.208 64mm；15Hz到30Hz時，振幅緩慢變化，從0.208 64mm到0.123 87mm。在不同激振頻率的荷載作用下，位移響應最大峰值出現在頻率為10Hz的情況，它對果樹枝干的動態性能影響最大，即在樹枝端部產生最大振幅。觀察加速度-頻率響應曲線，共振頻率出現在9、11、26Hz。為此，撞擊式收獲機的激振頻率應設計為10Hz左右[15]。

3 沖擊實驗

進行室內擺錘沖擊實驗，沖擊部位在距根部10cm處，在距果樹根部80cm位置放置IEPE式電壓加速度傳感器，DHDAS5922N動態信號測試分析系統測量加速度-時間曲線[16]，應用MatLab軟件導出加速度-頻率曲線，如圖6所示。

圖6 主枝最高點位置響應曲線

實驗得到加速度響應與有限元得到加速度響應趨勢基本吻合，果樹Y型結構的共振頻率在10Hz左右。

4 結論

1)通過建立果樹枝干三維模型，進行模態分析，獲得樹干各階振型和固有頻率，預測果樹達到共振時的頻率。

2)在模態分析基礎上進行諧響應，獲得果樹枝干在頻率0～30Hz、幅值740N的激振力作用下的最大位移響應、加速度響應曲線。

3)通過對果樹枝干動力學特性的分析研究，快速有效地找到果樹固有頻率及加速度響應曲線。與實驗進行對比，驗證模型有效性，證明用有限元模擬實驗是可行的方法，可為有效地設計、調整采收機械提供必要的參數[17]。

[1] 王行仁.建模與仿真技術的發展和應用[J].機械制造與自動化,2010(1):1-6,45.

[2] 陳騰浩,馬振利,韓建榮.系統建模與仿真技術探討[J].中國儲運，2012(3):106-107.

[3] 鄭甲紅,毛俊超,韓冰冰.振動式采摘機振動夾持位置的仿真研究[J].陜西科技大學學報:自然科學版,2014(1):142-147.

[4] Zoltán Láng. A Fruit Tree Stability Model for Static and Dynamic Loading[J].Biosystems Engineering,2003,85(4):461-466.

[5] 周概容.概率論與數理統計[M].北京：高等教育出版社,2006.

[6] 地質部長春計算站方法組.相似系數與相關系數的關系[J].長春地質學院學報,1980(2):70-92.

[7] 曾其蘊,鮑賢镕.木材含水率測定方法的比較研究[J].木材工業,1990,4(2):20-23.

[8] 王維新.彎曲彈性模量計算公式的推導[J].人造板通訊,1996,5(5):14-16.

[9] 劉家泳,汪堯,戈曉,等.四點彎曲試驗研究[J].物理測試,2010(2):44-47.

[10] 賀磊盈.面向振動采收的果樹枝干三維重建方法及動力學特性研究[D].杭州：浙江理工大學,2014.

[11] 賀磊盈,武傳宇,杜小強.基于雙輪廓同步跟蹤的果樹枝干提取及三維重建[J].農業工程學報,2014,30(7):182-189.

[12] 許進峰.ANSYS Workbench 15.0[M].北京：電子工業出版社,2014.

[13] Zoltán Láng.A simple model for fruit tree shaking harvest[J].International Journal of Horticultural Science,2007,13(1):33-36.

[14] 路娜.基于ANSYS的煤炭采樣機采樣臂動力學分析[D].西安：長安大學,2014.

[15] 邵景范,李文斌.基于ANSYS的電鍍CBN砂輪模態及諧響應分析[J].機械工程與自動化,2011(2):67-69.

[16] 張敏敏.果樹粘彈性振動結構模型的研究[D].烏魯木齊：新疆大學,2015.

[17] 陶嗣巍.樹木風振特性試驗研究與有限元分析[D].北京:北京林業大學,2013.

The Finite Element Analysis of the Vibration of Fruit Trees Based on ANSYS Workbench

Lv Menglu,Wang Chunyao, Luo Jianqing,Wei Tingpeng，Liu Zilong

Based on mechanical vibration of fruit trees as the object of study, we do kinetics research on mechanical vibration of fruit trees by using ANSYS Workbench. The 3D model of fruit trees was constructed by Pro/E 5.0 parameter method, simplified model was transmitted to ANSYS Workbench.The modal analysis is carried out to obtain the natural frequencies and vibration

fruit tree model; Pro/E; ANSYS Workbench; finite element analysis; vibration

2016-03-20

國家自然科學基金項目(51465054)

呂夢璐(1990- )，女，沈陽人，碩士研究生，(E-mail)lvmenglu@126.com。

王春耀(1956- )，男， 四川萬源人，教授，碩士生導師，(E-mail)wangchun_yao@126.com。

S225.93

A

1003-188X(2017)02-0037-05

-mode vectors of the first ten orders; its harmonic response analysis is carried out to calculate the resonance frequencies and amplitudes of the fruit trees, and the position of maximum resonance amplitude.Doing impact test in the lab, we can find resonance frequency, and the analysis results of finite element were verified by comparing the method of finite element analysis,so it can provide a theoretical basis for the design and optimization of mechanical vibrators.