李文春+王斌+劉曉麗

摘要：對果園避障旋耕機旋耕刀片特殊工況下（機組前進作業(yè)時機架附加橫向運動）的受力狀況進行了理論分析與計算，獲得了單刀受力的數(shù)學(xué)模型。運用SolidWorks軟件建立了旱地用C型旋耕刀實體模型，基于ANSYS Workbench分別對旋耕左刀和右刀進行了正常工況與避障作業(yè)2種情況下的應(yīng)力、應(yīng)變與變形仿真分析。仿真結(jié)果表明，避障動作改變了旋耕刀的受力狀況，但對旋耕刀片最大應(yīng)力、應(yīng)變、變形的位置沒有產(chǎn)生影響。仿真結(jié)果為分析特殊工況下旋耕刀失效形式、安全性能提供了方法，也為研究設(shè)計新型旋耕刀片提供了理論依據(jù)。

關(guān)鍵詞：果園避障旋耕機；C型旋耕刀；ANSYS Workbench；有限元分析

中圖分類號：S222.3 文獻標(biāo)志碼：A 文章編號：1002-1302（2017）01-0193-04

旋耕機具有很強的切土、碎土能力，耕后地表平整、松碎，又不會破壞深層土壤結(jié)構(gòu)，與目前推行的免耕、少耕的保護性耕作模式相適應(yīng)，發(fā)展迅速，是運用最廣泛的中耕作業(yè)機械。但普通旋耕機果園作業(yè)時，受耕作對象的局限會在株距之間形成一條非耕作區(qū)，需要人工進行株間除草松土等二次作業(yè)[1]。為了使工作部件能夠自動伸入株間進行松土除草作業(yè)，筆者設(shè)計了一種液壓式果園避障旋耕機。該旋耕機機架右前方設(shè)置有一個感應(yīng)觸桿，作業(yè)時觸桿首先碰觸到果樹，進而控制機架上液壓缸動作，液壓缸伸出帶動旋耕機架橫向退回，完成避障動作[2-3]。

旋耕刀輥是果園避障旋耕機的主要工作部件，也是旋耕作業(yè)中承受載荷最大的工作部件。通常情況下，由于旋耕作業(yè)時，機具只進行等幅寬作業(yè)，刀輥沿刀軸方向的受力很小，在分析刀片受力時只側(cè)重考慮垂直方向的受力[4-7]；但果園避障旋耕機旋耕作業(yè)時，遇到果樹能夠自動避障，旋耕刀輥在垂直于機組前進方向有明顯的橫向運動，因此，軸向受力不能忽略。在果園避障旋耕機具有橫向避讓動作的工況下，對旋耕刀片進行有限元分析可以計算刀片的結(jié)構(gòu)剛度、靜應(yīng)力等靜態(tài)參數(shù)，有利于分析旋耕刀的失效形式、失效位置，驗證旋耕刀片的安全性能，改善旋耕刀的設(shè)計參數(shù)，為研究設(shè)計新型旋耕刀提供理論依據(jù)。

1 果園避障旋耕機旋耕刀片受力分析

果園避障旋耕機避障作業(yè)時，刀輥受到的阻力來源于各組刀片阻力的合成，旋耕刀單刀的工作受力情況如圖1所示。旋耕刀切削過程中，刀具和土壤之間的相互作用大致可分為3個部分：切削面ABE產(chǎn)生的切削阻力及克服側(cè)切削刃與待切削土壤之間的摩擦力、切削面AEFD產(chǎn)生的土壤剪切阻力和刀具前刃切削中土壤金屬的摩擦力、刀盤沿軸向移動時刀片與土壤間摩擦力及已切土壤與未切土壤之間的摩擦阻力和黏結(jié)阻力等。

上述阻力中的前兩部分即為普通旋耕刀切削土壤時的切削阻力，第三部分為果園避障旋耕機避障作業(yè)時受到的軸向附加阻力，為了計算方便，把兩部分力分開計算，即總阻力為：

Ff=Fxy+Fz。

式中：Fxy=Fqianyin=Zbnh0k0，Z為刀盤數(shù)，bn為單體刀寬，h0為設(shè)計耕深，k0為土壤旋耕比阻[8]；Fz=（1+a）（μ1m1g+μ2m2g），μ1為旋耕刀與土壤的摩擦系數(shù)，m1為機架總質(zhì)量，μ2為土壤的內(nèi)摩擦系數(shù)，m2為刀輥推動土壤的質(zhì)量，a為考慮已切土壤與未切土壤之間黏聚力賦予的阻力系數(shù)。

設(shè)計中有：Z=7，b=23 cm，h0=16 cm，k0=9.645 N/cm2，μ1=0.6，μ2=0.75，m1=300 kg，m2=100 kg，a取0.5，帶入各數(shù)據(jù)得旋耕刀單刀切土阻力Fxy=1 218 N，F(xiàn)z=182 N。

2 旋耕刀片實體建模

新疆南疆果園土壤多為沙質(zhì)土壤，耕作比阻小，由于C型旋耕刀兼具直角刀剛性好、碎土能力強和彎刀不易纏草的特點，適宜旱地作業(yè)。同時，使用C型旋耕刀耕作時土壤有明顯的分層和流動現(xiàn)象，土壤圍繞旋轉(zhuǎn)軸從刀片凹面的前部和后部翻出形成運動流，兼具埋草效果[9]。因此，適用于新疆南疆果園土壤耕作。

在建模過程中，考慮C型旋耕刀結(jié)構(gòu)復(fù)雜，ANSYS軟件的建模能力相對較弱，但提供了豐富的數(shù)據(jù)接口程序，使得在其他3D軟件中建立模型后導(dǎo)入ANSYS變得簡單方便。因此，本研究選擇在SolidWorks軟件中進行三維建模，再通過SolidWorks和ANSYS的數(shù)據(jù)接口，解決ANSYS建模困難的問題。具體實施方法如下：（1）建立SolidWorks與ANSYS的軟件關(guān)聯(lián)；（2）采用描點法測繪C型旋耕刀結(jié)構(gòu)參數(shù)，運用SolidWorks軟件繪制旋耕刀草圖，并通過拉伸、彎曲、切除等命令建立C型旋耕刀片實體模型；（3）從SolidWorks主面板菜單中打開ANSYS Workbench16.0，則ANSYS Workbench 16.0共享獲得C型旋耕刀模型（圖2）。

3 旋耕刀片有限元分析

有限元分析（Finite Element Analysis，F(xiàn)EA）是指利用數(shù)學(xué)近似的方法模擬真實物理系統(tǒng)中的幾何約束和載荷工況，進而用有限數(shù)量的未知量去逼近無限未知量的真實系統(tǒng)分析方法?，F(xiàn)階段應(yīng)用最廣泛的有限元分析軟件是美國ANSYS公司研制的大型通用有限元分析軟件，能與市場上大多數(shù)計算機輔助設(shè)計（computer-aided design，CAD）軟件接口通信，實現(xiàn)數(shù)據(jù)共享和交換，已廣泛應(yīng)用于結(jié)構(gòu)力學(xué)、流體力學(xué)、熱力學(xué)、聲學(xué)、電磁學(xué)、電路學(xué)、化學(xué)化工反應(yīng)等領(lǐng)域[10]。

3.1 材料設(shè)置與網(wǎng)格劃分

旋耕刀是旋耕機工作時直接與土壤相互作用的部件，通過高速旋轉(zhuǎn)和前進運動，將土壤切碎、拋出，為易損件。旋耕刀的材料要求硬度高、耐磨、韌性好，工業(yè)上常選取60 Si2Mn材料并進行熱處理作為旋耕刀的原料[11]。其參數(shù)見表1。

旋耕刀的單元劃分采用ANSYS Workbench自帶的智能網(wǎng)格劃分，該方法可以根據(jù)模型的幾何關(guān)系，自動將網(wǎng)格劃分得稀疏得當(dāng)。單元尺寸控制在0.005 m，劃分好的旋耕刀網(wǎng)格如圖3所示。

3.2 約束與載荷施加

由前述分析可知，旋耕刀在整機避障時受力情況為：在受正常作業(yè)切削阻力的同時受到軸向阻力，因此在旋耕刀正切刃和側(cè)切刃上施加垂直于刃口方向的1 218 N的力，在旋耕刀面沿刀軸方向施加182 N的力（圖4）。旋耕刀通過安裝孔固定在刀盤上，因此在模型安裝孔的位置施加固定約束。

3.3 求解與分析

添加總變形、應(yīng)力和應(yīng)變分析項目進行求解，求解結(jié)果見圖5。

應(yīng)力云圖結(jié)果如圖5所示：果園避障旋耕機正常作業(yè)時左刀最大等效應(yīng)力為88.60 MPa，右刀最大等效應(yīng)力為 87.75 MPa；避障作業(yè)時，左刀最大等效應(yīng)力為114.14 MPa，右刀最大等效應(yīng)力為108.14 MPa，均遠(yuǎn)遠(yuǎn)小于材料屈服強度值， 說明旋耕刀的結(jié)構(gòu)仍然存在很大的優(yōu)化空間。最大應(yīng)力點所在位置為遠(yuǎn)離切土點的安裝孔位置和刀盤與刀接觸位置，即刀與刀盤連接處會產(chǎn)生較大應(yīng)力集中，與實際工作中刀盤的斷裂形式相吻合，驗證了仿真模型的準(zhǔn)確性。應(yīng)力最小值位于旋耕刀正切部。避障動作對左右刀片應(yīng)力值影響不同是因為避障時2個刀片所受到土壤給與的軸向力受力面不同且方向相反。

應(yīng)變云圖結(jié)果如圖6所示：果園避障旋耕機正常作業(yè)和避障作業(yè)時，旋耕刀最大應(yīng)變和最小應(yīng)變位置均與相同工況下最大最小應(yīng)力位置相同，正常作業(yè)時左刀最大應(yīng)變?yōu)?6.44×10-4 m/m，右刀最大應(yīng)變?yōu)?.26×10-4 m/m；避障作業(yè)時左刀最大應(yīng)變?yōu)?.79×10-4 m/m，右刀最大應(yīng)變?yōu)?.26×10-4 m/m，均在旋耕刀與刀盤連接處，說明此處承受載荷較大。

變形云圖結(jié)果如圖7所示：果園避障旋耕機正常作業(yè)和避障作業(yè)時，旋耕刀最大變形位置均在旋耕刀刀尖處，說明作業(yè)中刀尖處剛度表現(xiàn)最差。正常作業(yè)時左刀最大變形量為0.346 mm，右刀最大變形量為0.508 mm；避障作業(yè)時左刀最大變形量為0.359 mm，右刀最大變形量為0.543 mm，在剛度允許范圍內(nèi)。

4 結(jié)論

以果園避障旋耕機C型旋耕刀為研究對象，分析了旋耕機避障作業(yè)時旋耕刀片的受力狀況，運用SolidWorks軟件建立了旋耕刀的三維模型，基于ANSYS Workbench有限元分析軟件對旋耕刀片在正常工況和避障作業(yè)工況下分別進行了應(yīng)力、應(yīng)變和變形分析。分析結(jié)果表明：果園避障旋耕機正常工況下和避障作業(yè)時，旋耕刀片的最大應(yīng)力和應(yīng)變區(qū)域均在刀片與刀盤連接處，與實際工作過程中旋耕刀斷裂處一致，表明該處為旋耕刀在工作過程中最薄弱的環(huán)節(jié)；旋耕刀片總變形最大處均在旋耕刀的正切部，表明旋耕刀正切部剛度最差，應(yīng)采用工藝手段加強處理；由于避障作業(yè)時旋耕刀輥是沿刀軸方向向左回縮，左刀正面受力，右刀背面受力，因此該工況對左刀片應(yīng)力、應(yīng)變的影響大于右刀片，變形的影響小于右刀片，以此推斷當(dāng)避障作業(yè)結(jié)束后刀輥橫向向右伸出時，工況對右刀片應(yīng)力、應(yīng)變的影響大于左刀片，變形的影響小于左刀片，這在后續(xù)仿真中得到了證實。

參考文獻：

[1]萬 暢，李文春，畢新勝，等. 淺談國內(nèi)果園旋耕機存在的問題和解決思路[J]. 新疆農(nóng)機化，2013（5）：15-17.

[2]羅進軍，王 斌，畢新勝，等. 果園避障旋耕機的研究與設(shè)計[J]. 農(nóng)機化研究，2016，38（1）：138-141.

[3]畢新勝，查鑫宇，高小軍，等. 一種果園旋耕機：CN203537758U[P]. 2014-04-16.

[4]張靈芝，朱為國. 基于Solidworks旋耕刀實體建模與有限元分析[J]. 農(nóng)機化研究，2013，35（11）：26-29.

[5]高玉芝，王君玲，雷曉柱. 基于Solidworks的旋耕刀實體建模與有限元分析[J]. 農(nóng)業(yè)機械，2010（5）：136-137.

[6]葛 云，吳雪飛，王 磊，等. 基于ANSYS微型旋耕機旋耕彎刀的應(yīng)力仿真[J]. 石河子大學(xué)學(xué)報（自然科學(xué)版），2007，25（5）：627-629.

[7]朱留憲，楊 玲，朱 超，等. 基于ANSYS Workbench的微耕機旋耕刀有限元分析[J]. 機械研究與應(yīng)用，2014，27（1）：88-89.

[8]中國農(nóng)業(yè)機械化科學(xué)研究院.農(nóng)業(yè)機械設(shè)計手冊[M]. 北京：中國農(nóng)業(yè)科學(xué)技術(shù)出版社，2007.

[9]斯瑞斯克·切克凱特，坦亞·尼亞瑪帕. 不同旋耕刀的轉(zhuǎn)矩變化與耕作比能[J]. 農(nóng)業(yè)工程，2011，1（1）：129-130，132.

[10]CAE應(yīng)用聯(lián)盟. ANSYS Workbench 15.0有限元分析從入門到精通[M]. 北京：機械工業(yè)出版社，2014.

[11]彈簧鋼：GB/T 1222—2007[S]. 北京：中國標(biāo)準(zhǔn)出版社，2007.顏丙囤，梁守真，王 猛，等. 花生葉綠素含量的高光譜遙感估算模型研究[J]. 江蘇農(nóng)業(yè)科學(xué)，2017，45（1）：197-200.