孟兆新，范恒博，王亮亮

(東北林業(yè)大學(xué) 機(jī)電工程學(xué)院，哈爾濱 150040)

樹木移植機(jī)按照配用動力可分為牽引式、自行式和懸掛式，按照鏟刀的形狀可分為U型鏟刀式、彎鏟式和鍬鏟式，按照鏟刀的數(shù)量可以分為三鏟式和四鏟式。目前廣泛應(yīng)用的是四片鍬鏟式樹木移植機(jī)。工作時，鏟刀被固于樹木的根部，在自身液壓系統(tǒng)或者外界設(shè)備系統(tǒng)提供的液壓力作用下開始入土剪切土壤和樹根，合攏的鏟刀使得樹根部的土壤形成一個半圓錐形的土球，土球和樹根部分被鏟刀提升到地面，完成樹木移植工作。

鏟刀結(jié)構(gòu)是樹木移植機(jī)的關(guān)鍵結(jié)構(gòu)，但是沒有計算鏟刀在入土過程中的阻力計算模型，所以在設(shè)計過程無法獲得在設(shè)計過程中的相關(guān)數(shù)據(jù)的準(zhǔn)確值。由于鏟刀結(jié)構(gòu)在入土過程中的受力情況比較復(fù)雜，運(yùn)用ANSYS軟件分析鏟刀在入土過程中與土壤的接觸過程，從ANSYS的后處理模塊中可以直觀地看到鏟刀結(jié)構(gòu)各個部分的受力情況，從而知道鏟刀結(jié)構(gòu)的最大受力點(diǎn)和最小受力點(diǎn)，為以后的鏟刀的改進(jìn)設(shè)計提供數(shù)據(jù)和理論依據(jù)，從而提高鏟刀的壽命和工作效率，并且預(yù)測鏟刀可能出現(xiàn)的問題。

1 鏟刀入土過程力的理論分析

本文所研究的是四鏟式樹木移植機(jī)的鏟刀結(jié)構(gòu)，其結(jié)構(gòu)如圖1所示，四鏟式挖樹完成之后的土球?yàn)榘雸A錐形土球。

圖1 鏟刀的模型

在研究過程中，將鏟刀結(jié)構(gòu)的入土過程簡化，簡化的結(jié)果是：鏟刀在一定的載荷作用下，在入土角的存在下與水平面成一定的角度勻速下鏟，對土壤進(jìn)行剪切。鏟刀結(jié)構(gòu)與土壤間的相互作用力F為：

F=Fu+Fv+Fa

(1)

式中：Fu為鏟刀本身位移產(chǎn)生的相互作用力；Fv為鏟刀本身速度產(chǎn)生的相互作用力；Fa為鏟刀本身加速度產(chǎn)生的相互作用力。

由于在前面已經(jīng)簡化，鏟刀結(jié)構(gòu)是勻速入土，所以忽略由于加速度產(chǎn)生的慣性力Fa。由于此模型適合于理論研究，但是具體的數(shù)值難以計算，所以本文采用顧正平等人得出的直鏟式下鏟阻力的經(jīng)驗(yàn)公式[2-3]：

F=0.130 7x1.138 4γ0.755 3c0.244 7(13.927 5+0.563 3×100.024 3β)

(2)

式中：x為鏟刀結(jié)構(gòu)的下鏟位移量，取鏟刀的鏟高(1 m)；γ為土壤容重(1.34 g/cm3)；c為土壤內(nèi)聚力(44.126 445 kPa)；β為鏟刀體圍角的半角(60°)。

在本文中，鏟刀結(jié)構(gòu)的入土角取45°，通過上式計算得出鏟刀結(jié)構(gòu)在入土過程中的阻力值為：67.2kN。

2 鏟刀有限元模型的分析

對模型進(jìn)行分析的步驟如下：

(1)對模型進(jìn)行網(wǎng)格劃分，得到網(wǎng)格劃分如圖2所示(劃分網(wǎng)格數(shù)9 018，結(jié)點(diǎn)數(shù)為18 579)。

圖2 網(wǎng)格劃分

(2)對模型進(jìn)行約束和加載，在鏟刀的入土過程中，鏟刀面的內(nèi)外側(cè)表面都始終與土壤接觸，并且有接觸力，鏟刀內(nèi)側(cè)的力來自于半圓錐形土球的重力，而且受到了F的力，鏟刀外側(cè)的力、重力和F力平衡。將各個力加載到模型中并且進(jìn)行約束[4-9]。

(3)對模型進(jìn)行優(yōu)化。通過ANSYS的求解，得到的應(yīng)力效果如圖3所示。

從圖中可以看出：受力最大的部分是鏟刀入土側(cè)的兩個刀尖，最大的力是2.75×103MPa；其次是施加力的部位、施加力部位與鏟刀入土側(cè)之間組成的矩形區(qū)域，受到的力在30.5～91.6MPa；受力最小的是鏟刀施加力部位以上的部分和矩形以外的部分，受到的力9 057Pa。

圖3 模型的應(yīng)力圖

從應(yīng)力圖中得到結(jié)論：鏟刀入土側(cè)的兩個刀尖應(yīng)該做成圓角，避免應(yīng)力的集中；矩形區(qū)域是受力交大的部分，應(yīng)該加鏟槽以增加剛度；施加力部位以上的部分的鏟刀應(yīng)該減少厚度，以減小鏟刀在入土過程中收到的阻力。

通過ANSYS的求解，得到效果圖如圖4所示。從中找到不同步優(yōu)化過程中的效果圖和最后優(yōu)化的效果圖如圖5所示。

圖4 第2步的優(yōu)化效果

圖5 第6步的優(yōu)化效果

拓?fù)鋬?yōu)化的最終效果如圖6所示。

圖6 拓?fù)鋬?yōu)化的最終效果

(1)圖6所示的是拓?fù)鋬?yōu)化的最終效果圖為鏟刀的改進(jìn)設(shè)計提供了數(shù)據(jù)和理論依據(jù)，解決了設(shè)計人員沒有準(zhǔn)確數(shù)據(jù)的難題。

(2)在工業(yè)上可以根據(jù)此圖對鏟刀各部位的厚度進(jìn)行合理分布或者改進(jìn)原有的鏟刀，為提高鏟刀的效率和壽命的設(shè)計提供了參考。改進(jìn)措施有：在虛密度大的部分予以加厚或者加以刀槽以加強(qiáng)其強(qiáng)度，對于虛密度較小的部分使其厚度減小以減小鏟刀的入土阻力。

(3)從應(yīng)力云圖中可以看出，鏟刀的中心偏上部位是受到最大應(yīng)力的地方，以應(yīng)力最大的部分為核心逐漸向外擴(kuò)散，應(yīng)力逐漸變小，而鏟刀的刀尖部位沒有明顯的大應(yīng)力。所以可以適當(dāng)?shù)脑黾隅P刀刀尖的鋒利程度，以提高鏟刀剪切土壤和樹木根部的能力。

(4)合理的分布厚度節(jié)約了材料，提高了鏟刀的剛度，提高鏟刀的壽命；減小了刀尖的厚度使得入土阻力減小，提高了工作效率；并且根據(jù)應(yīng)力圖可知鏟刀的受力情況，預(yù)測鏟刀可能出現(xiàn)的問題。

3 結(jié)束語

通過對鏟刀結(jié)構(gòu)的拓?fù)鋬?yōu)化分析，得到鏟刀結(jié)構(gòu)的偽密度云圖，為鏟刀的改進(jìn)設(shè)計提供數(shù)據(jù)和理論依據(jù)，解決設(shè)計人員沒有準(zhǔn)確數(shù)據(jù)的難題；同時對鏟刀質(zhì)量的合理分布提供了模型，保證了材料的最大的利用率，避免鏟刀改進(jìn)設(shè)計的盲目性。

【參 考 文 獻(xiàn)】

[1]宋海兵，王春宏.樹木移植機(jī)下鏟過程土壤與鏟刀接觸受力分析[J].林業(yè)機(jī)械與木工設(shè)備，2013，41(2)：42-43

[2]張英彥，顧正平.直鏟式樹木移植機(jī)下鏟阻力的探討[J].北京林業(yè)大學(xué)學(xué)報，1988，10(2)：38-45

[3]顧正平，張英彥.模型鏟預(yù)測樹木移植機(jī)下鏟阻力的研究[J].林業(yè)科學(xué)，1988，24(2)：185-194.

[4]周 寧.ANSYS機(jī)械工程應(yīng)用實(shí)例[M].北京：中國水利水電出版社，2006.

[5]孫 斌.小型振動式樹木移植機(jī)設(shè)計與仿真[D].西安：西安理工大學(xué)，2011.

[6]薛子萱，呂新民，唐衛(wèi)衛(wèi).螺旋刀具土壤切削過程模擬分析-基于ANSYS/LS-DYNA971[J].農(nóng)機(jī)化研究.2011，4(4)：13-16.

[7]姜 濤.基于SPH法的板結(jié)土壤反旋鑿切有限元仿真研究[D].杭州：浙江工業(yè)大學(xué)，2009.

[8]李曉祥，林嘉棟，侯業(yè)華.應(yīng)用ANSYS對齒輥破碎機(jī)齒輥結(jié)構(gòu)優(yōu)化及應(yīng)力分析[J].林業(yè)機(jī)械與木工設(shè)備，2013，41(3)：31-34.

[9]李芝茹，吳曉峰，羅菊英，等.基于Ansys Workbench 的BYJ-800型油茶嫁接機(jī)穗木切削機(jī)構(gòu)的仿真[J].林業(yè)機(jī)械與木工設(shè)備，2013，41(10)：17-20+25.