趙德金，郭艷玲，李志鵬

(1.東北林業(yè)大學(xué)機(jī)電工程學(xué)院，哈爾濱150040;2.延邊大學(xué)工學(xué)院，吉林省延吉133002;3.東北林業(yè)大學(xué)交通學(xué)院，哈爾濱150040)

目前，利用挖樹機(jī)實(shí)施城市綠化、是保證城市的生態(tài)，加速城市鄉(xiāng)鎮(zhèn)綠化建設(shè)最有利的工具，挖樹鏟是挖樹機(jī)的重要工作部件，挖樹鏟在含有碎石、樹根、雜草以及其根莖的土壤中作業(yè)，所受的工作阻力較大，對挖樹鏟的剛度和抗磨損能力要求很高，目前常用的挖樹鏟按照結(jié)構(gòu)和入土形式可分為U形鏟、弧形鏟和直鏟，直鏟按形狀又分為V字形鏟和曲面鏟等［1－3］。本文主要對直鏟進(jìn)行設(shè)計(jì)和有限元分析。

挖樹鏟在挖樹過程中所受的阻力計(jì)算很復(fù)雜，采用傳統(tǒng)的計(jì)算方法很難計(jì)算，主要受樹木周圍的土壤的實(shí)際情況影響很大，對于直鏟而言，挖樹阻力主要是挖樹鏟與土壤的切削阻力、挖樹鏟與土壤摩擦阻力、挖樹鏟與樹木根莖以及草根的切斷阻力，再就是對于含有碎石的土壤的碎石對挖掘鏟所產(chǎn)生的阻力等，挖樹鏟受到阻力后，發(fā)生變形以及應(yīng)力分布的情況，采用傳統(tǒng)的方法也很難確定，一般采用分組實(shí)驗(yàn)確定經(jīng)驗(yàn)公式，并且計(jì)算結(jié)果主要適用于實(shí)驗(yàn)所采用的土壤條件，一般可計(jì)算出挖樹鏟的受力情況，對于挖樹鏟的應(yīng)力分布及變形仍然不明確。對于這種情況，本文采用有限元軟件對挖樹鏟進(jìn)行靜力學(xué)分析來解決挖樹鏟在工作過程中應(yīng)力分布及變形情況［4－7］。

1 建立挖樹鏟模型與有限元靜力學(xué)分析

挖樹鏟是樹木移植機(jī)的主要工作部件之一，其結(jié)構(gòu)應(yīng)該簡單，制造方便，應(yīng)經(jīng)可能保留樹木的根莖，從而保證被移植樹木的成活率，直鏟具有結(jié)構(gòu)簡單，容易制造的特點(diǎn)，在國內(nèi)外樹木移植機(jī)械中已經(jīng)得到了廣泛的應(yīng)用。

圖1 曲面鏟和V字形鏟結(jié)構(gòu)圖Fig.1 The structure of surface spade and V-shaped spade

1.1 挖樹鏟的設(shè)計(jì)

利用PTC Creo軟件建立挖樹直鏟三維實(shí)體模型如圖1所示。根據(jù)樹木移植機(jī)的設(shè)計(jì)要求確定鏟的結(jié)構(gòu)尺寸，利用PTC Creo的曲面功能設(shè)計(jì)曲面，然后利用曲面加厚功能生成挖樹鏟實(shí)體模型。然后保存為IGES格式，用于Ansys Workbench進(jìn)行有限元分析。

1.2 挖樹鏟的網(wǎng)格劃分

在進(jìn)行有限元分析之前，首先將挖樹鏟的IGES模型讀入Ansys Workbench系統(tǒng)中，如圖2所示。在進(jìn)行網(wǎng)格劃分之前，首先定義材料屬性。

圖2 Ansys中挖樹鏟模型Fig.2 The IGES model of the tree spade in Ansys Workbench

材料與屬性:采用Q235碳素結(jié)構(gòu)鋼，彈性模量E=2.06×105N/mm2;密度ρ=7.8×10－6kg/mm3;泊松比0.3。

本次設(shè)計(jì)挖掘鏟的厚度尺寸為6 mm，根據(jù)機(jī)械設(shè)計(jì)手冊，碳素結(jié)構(gòu)鋼的鋼材厚度小于等于16 mm時，其屈服極限σS=235 MPa。

采用Ansys Workbench對挖樹鏟進(jìn)行網(wǎng)格劃分，默認(rèn)方式自動劃分網(wǎng)格如圖3所示，曲面挖樹鏟有限元模型產(chǎn)生2477個單元，5342個節(jié)點(diǎn)。V字形挖樹鏟有限元模型產(chǎn)生2159個單元，4705個節(jié)點(diǎn)。

圖3 Ansys Workbench中挖樹鏟模型網(wǎng)格劃分Fig.3 The meshes of the tree spade in Ansys Workbench

1.3 挖樹鏟有限元靜力學(xué)分析

對于V字形鏟和和曲面鏟而言，鏟上表面的任意點(diǎn)的法線都與下鏟阻力互相垂直，這樣就使鏟的下鏟阻力主要是鏟刀刃與土壤以及樹木根莖等的阻力，而挖樹鏟上下鏟面與鏟的土壤切削阻力相比很小，本次研究將載荷加在鏟的切削刃上，方向與挖樹鏟動力方向相反，可以使用PTC Creo軟件的分析功能計(jì)算V字形挖樹鏟的截面的面積如圖4所示。同樣的方法也可以計(jì)算曲面鏟的截面面積。

在施加載荷和求解之前，先假定挖樹鏟刃口各處所受阻力方向和大小一致，施加載荷的大小可以參考下鏟阻力經(jīng)驗(yàn)公式［8］:

圖4 V字形鏟橫截面面積Fig.4 The cross-section area of the V-shaped spade

式中:x為鏟刀下鏟的位移量，m;γ為土壤容重，g/cm3;c為土壤內(nèi)聚力，Pa;β為鏟刀體圍角的半角，°。

表1 入土角為60度是土壤對兩種挖樹鏟刃的壓強(qiáng)Tab.1 The pressures on the tree spade when the digging angle is 60°

通過計(jì)算取下鏟阻力為27，這樣就可以確定挖樹鏟刃口壓強(qiáng)見表1，然后將挖樹鏟驅(qū)動銷軸孔曲面固定，限定挖樹鏟的自由度，然后將計(jì)算好的壓強(qiáng)施加在挖樹鏟刃上，使用Static Structural求解。

2 挖樹鏟的應(yīng)力計(jì)算結(jié)果

經(jīng)過求解得到曲面挖樹鏟應(yīng)力圖如圖5所示，應(yīng)變圖如圖6所示。V字形鏟的應(yīng)力圖如圖7所示，應(yīng)變圖如圖8所示。應(yīng)力應(yīng)變極限值見表2。

圖5 曲面挖樹鏟的應(yīng)力圖Fig.5 The stress distribution of surface tree spade

圖6 曲面挖樹鏟的變形圖Fig.6 The strain of surface tree spade

圖7 V字形挖樹鏟的應(yīng)力圖Fig.7 The stress distribution of V-shaped spade

圖8 V字形挖樹鏟的變形圖Fig.8 The strain of surface V-shaped spade

表2 挖樹鏟應(yīng)力應(yīng)變極限值α(60°)Tab.2 The limit value α of the stress and the strain of the tree spade(digging angle=60°)

3 結(jié)論

(1)從兩個挖樹鏟的應(yīng)力分布情況出，應(yīng)力和最大應(yīng)力主要分布在施加驅(qū)動載荷的銷軸孔附近，曲面挖樹鏟鏟最大應(yīng)力為0.45 MPa，符合設(shè)計(jì)要求，V字形鏟最大應(yīng)力為1.659 MPa，也符合設(shè)計(jì)要求，從數(shù)值上看曲面挖樹鏟性能要優(yōu)于V字形鏟。但從結(jié)構(gòu)上看V字形鏟結(jié)構(gòu)簡單，與曲面鏟相比，容易制造。

(2)最大變形量V字形鏟為9.1431×10－6m，主要在發(fā)生鏟尖局部區(qū)域。曲面挖樹鏟1.283×10－5m，主要發(fā)生在鏟尖和切削刃口上，兩者相對于鏟的長度而言，變形量很小，可以忽略不計(jì)。

(3)總體上看，曲面鏟雖然制造困難一些，但是從保留樹根量和承載能力上看都優(yōu)于V字形鏟。

［1］ Lemar Tree Spades［EB/02］.http://www.lemarweb.com/site/

［2］ Dutchman Truck Spade［EB/02］.http://www.dutchmantreespade.com/

［3］ I Dutchman Truck Spade［EB/02］.http://www.dutchmantreespade.com/

［4］賈晶霞，楊德秋，王昕偉，等.基于Solidworks的各型馬鈴薯挖掘鏟應(yīng)力研究［J］.農(nóng)業(yè)工程，2011，1(1):102－104.

［5］呂金慶，陳春富，李世柱.馬鈴薯挖掘機(jī)挖掘鏟的設(shè)計(jì)［J］.農(nóng)業(yè)機(jī)械化研究，2006，26(3):105 －106.

［6］賈晶霞，劉漢武，郝新明，等.馬鈴薯收獲機(jī)挖掘鏟有限元靜力學(xué)分析［J］.農(nóng)業(yè)機(jī)械學(xué)報，2006，37(9):86 －88.

［7］張 凡，鮑際平.重卡離合器壓盤熱應(yīng)力有限元分析與改進(jìn)［J］.林業(yè)機(jī)械與木工設(shè)備，2011，39(4):23 －26.

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