韓 超，楊有剛，張淑偉

(西北農(nóng)林科技大學(xué) 機械與電子工程學(xué)院，陜西 楊凌 712100)

懸掛式鏵式犁機架有限元分析及基于等強度的優(yōu)化設(shè)計

韓 超，楊有剛，張淑偉

(西北農(nóng)林科技大學(xué) 機械與電子工程學(xué)院，陜西 楊凌 712100)

合理的鏵式犁機架結(jié)構(gòu)可以保證在相同的牽引力作用下拖拉機能夠拖動更大的犁鏵，機架所用材料更少，同等作業(yè)條件下機架使用壽命更長。為此，利用ANSYS有限元分析軟件對經(jīng)過長期生產(chǎn)實際檢驗合格的懸掛式鏵式犁機架進行靜力學(xué)及靈敏度分析，模擬真實環(huán)境下的工作狀態(tài)，調(diào)整設(shè)計并使用基于等強度的設(shè)計方法優(yōu)化結(jié)構(gòu)。結(jié)果表明：改進后的機架結(jié)構(gòu)不但節(jié)省材料、滿足了設(shè)計要求，而且在原基礎(chǔ)上擴大了承載范圍，在相同結(jié)構(gòu)下，改進后的機架能夠承載更大的載荷。

鏵式犁；有限元；機架；結(jié)構(gòu)優(yōu)化；標(biāo)準(zhǔn)化設(shè)計

0 引言

犁耕是大田農(nóng)業(yè)生產(chǎn)中最重要的工作環(huán)節(jié)之一，鏵式犁是農(nóng)業(yè)生產(chǎn)的重要農(nóng)機具之一。鏵式犁工作時，主要依靠由犁壁與犁鏵組成的主犁體曲面對土壤進行入土、切碎、破碎和土垡翻轉(zhuǎn)來實現(xiàn)作業(yè)[1]。良好的機架結(jié)構(gòu)有助于主犁體克服來自土壤的粘滯阻力，提升翻垡效率，提升犁體承載能力，提高耕作效率，降低耕作者的勞動強度。因此，本文利用ANSYS有限元分析軟件，對經(jīng)過長期生產(chǎn)實際檢驗合格的懸掛式鏵式犁機架進行靜力學(xué)及靈敏度分析，并根據(jù)實際工作環(huán)境下的狀態(tài)調(diào)整設(shè)計，并基于等強度的設(shè)計方法進行優(yōu)化結(jié)構(gòu)，以期提供更實用的設(shè)計方案，達到保證工作強度需要、降低生產(chǎn)成本的目的。

1 懸掛式鏵式犁

1.1 基礎(chǔ)結(jié)構(gòu)分析

懸掛式鏵式犁主要是由三角牽引懸掛裝置、耕深調(diào)節(jié)裝置、機架及主犁體(幅寬1 050mm，耕深270mm)組成。其中，機架主要由三角懸掛牽引裝置、犁臂、固定螺栓及制動銷等零件構(gòu)成。機架梁局采用20×60的板材，主體材料為Q235碳鋼，切變模量為G=8×1011N/m，機架梁主要連接方式為焊接。

1)懸掛參數(shù)：懸掛軸至犁體支撐面的距離h=1.09m，上下懸掛點的距離H1=(犁架柱高度)懸掛軸長度L=0.1m，兩下懸掛點間距B=0.51m，兩下懸掛點與犁梁的相對位置L1=0.1m。

2)牽引方式：三鏵犁正牽引，在實際工作中用東方紅X850拖拉機牽引，速度范圍v=0.55～0.92m/s，標(biāo)定功率P=62.5kW ，額定轉(zhuǎn)速n=2 400r/min，輸出軸轉(zhuǎn)速為n=720r/min，動力輸出軸最大功率Pmax=55.5kW，牽引力為F≥23.3kN，計算時取F=23.3kN。

1.2 機架數(shù)據(jù)收集與修正

本文所選懸掛式鏵式犁機架梁數(shù)據(jù)來源于本校機械與電子工程學(xué)院農(nóng)機化實驗室，具體數(shù)據(jù)如圖1所示。

圖1 懸掛式機架梁實體數(shù)據(jù)Fig.1 Solid data of suspended frame beam

實地可測的犁體曲面與前進方向成32°夾角,犁的梁架高度H是指犁架下表面至犁底平面的空間高度。為了保證土壤垡片在犁架下翻轉(zhuǎn)，不產(chǎn)生堆土、堆草現(xiàn)象。一般依據(jù)矩形土垡的厚度(按最大耕深計算)與割茬高度的對角線高度計算，即

其中，b為犁體耕寬；amax為最大耕深；h為割茬高度?？傻脩覓焓界f式犁空間高度為

H=[b2+(amax+h)2]0.5=1.09m

2 簡化機架有限元模型

2.1 模型建立

本文采用Pro/E軟件建立機架部分的三維實體模型[2]，如圖2所示。將其導(dǎo)入到ANSYS有限元分析軟件中形成CAE模型[3]，為了準(zhǔn)確模擬懸掛式鏵式犁在工作過程中機架的受力狀態(tài)，采用ANSYS對其進行有限元分析，包括靜力學(xué)及靈敏度分析。由于在ANSYS 軟件中進行分析時圓角、倒角、小孔及凸臺等會降低網(wǎng)格質(zhì)量，使得計算數(shù)量急劇增加，所以在本模型中將有影響的小特征去除。懸掛式鏵式犁機架采用焊接而成，雖然結(jié)構(gòu)較為簡單，焊縫處又可近似認(rèn)為是剛性連接；但是，為了數(shù)據(jù)準(zhǔn)確性及更好地應(yīng)用等強度設(shè)計方法，在有限元分析時把焊縫單獨定義。

圖2 導(dǎo)入后機架的CAD模型Fig.2 The frame CAD model after import

根據(jù)機架材質(zhì)，設(shè)置材料屬性，定義機架的材料為Q235 碳鋼[4]，材料參數(shù)如表1所示。其中，彈性模量為E=2.11×1011N /m2，泊松比為μ=0.3，材料密度為ρ= 7 850kg/m3。

表1 Q235材料參數(shù)

2.2 劃分網(wǎng)格

網(wǎng)格劃分是線性靜力分析必不可缺少的一個環(huán)節(jié)，網(wǎng)格品質(zhì)的優(yōu)劣對結(jié)果的準(zhǔn)確性有著重要的影響。Ansys軟件中提供了多種網(wǎng)格劃分方法，分別是自動劃分法、四面體劃分法、六面體主導(dǎo)法、掃掠劃分法以及多區(qū)劃分法。為了保證計算精度，在本次分析中采用Automatic(自動劃分法)[4-6]。在Mesh 工具欄下載入Method，然后選擇待劃分的體，在窗口中的Geometry中單擊Apply按鈕，單擊Mesh設(shè)置Sizing尺寸控制網(wǎng)格精度，設(shè)置尺寸為3mm，點擊Generate Mesh生成網(wǎng)格，如圖3所示。

圖3 網(wǎng)格劃分Fig.3 Grid division

2.3 定義邊界條件

在進行邊界條件處理時，一種是不添加任何約束，另一種是根據(jù)結(jié)構(gòu)的實際工作狀況添加相應(yīng)約束[7]。為了準(zhǔn)確模擬鏵式犁在工作過程中的受力邊界條件，保留犁的幾何構(gòu)造，并對其劃分網(wǎng)格，施加最大耕深時的邊界條件。此時最大耕深amax=280mm，工作時拖拉機以0.5m/s 速度勻速前行，牽引力F=23.3kN以接觸壓力的形式加載在犁刀上，3個犁刀上的力P=7.767kN。另外，考慮到鏵式犁在設(shè)計作過程中通過三點懸掛系統(tǒng)與拖拉機相連，因此在計算時對機架與拖拉機三點懸掛連接處施加固定約束。施加邊界條件后如圖4所示，Rxy水平面內(nèi)受力示意圖如圖5所示。

圖4 懸掛式鏵式犁機架約束圖Fig.4 Suspension type plow frame constraint graph

圖5 懸掛式鏵式犁機架Rxy水平面內(nèi)受力示意圖Fig.5 Hanging plow frame Rxy horizontal plane force diagram

圖5中，Rx=F=23.3×103N，Ry=Rxcosθ=14.55×103N。

3 計算結(jié)果與分析

本文主要模擬實際生產(chǎn)情況中所受到的主要載荷，對懸掛式鏵式犁機架進行有限元靜力學(xué)和靈敏度分析。

3.1 靜力學(xué)分析

在Solution 中選擇Deformation>Total,在結(jié)構(gòu)樹中選中Deformation Total 并求解，可得位移等效云圖,如圖6 所示。在樹結(jié)構(gòu)圖中選擇Equivalent Stress ，然后求解可得其應(yīng)變等效云圖如圖7 所示，所得的Von-mises應(yīng)力云圖如圖8所示。

圖6 懸掛式鏵式犁機架位移等效云圖Fig.6 Suspension type plow frame equivalent displacement nephogram

圖7 懸掛式鏵式犁機架應(yīng)變等效云圖Fig.7 Suspension type plow frame equivalent strain nephogram

圖8 懸掛式鏵式犁機架Von-mises應(yīng)力云圖Fig.8 Suspension type plow frame equivalent Von-mises stress nephogram

由圖6可知：懸掛式鏵式犁的最大位移出現(xiàn)在第3個犁刀上，變形量約為14.3mm，主要是由于鋼架扭轉(zhuǎn)造成的。結(jié)合應(yīng)變云圖可以明顯地看出，這部分位移量是剛體位移。

由圖7可以看出：機架上的最大應(yīng)變主要出現(xiàn)在肋板連接位置。

由圖8可以看出：機架結(jié)構(gòu)上的大部分結(jié)構(gòu)的應(yīng)力都較小，集中在55MPa左右，遠(yuǎn)小于屈服極限235MPa。這說明，結(jié)構(gòu)可以進行優(yōu)化的空間很大。另外，最大應(yīng)力500MPa出現(xiàn)在立板的連接處，是因為采用剛性連接所致，不影響最后結(jié)果。

由于應(yīng)力較大的部位分別位于肋板、前面支架與犁的連接處、右前懸掛連接處，因此單獨對相應(yīng)的3部分進行拆分分析。肋板應(yīng)力較大處的云圖如圖9所示，前機架應(yīng)力較大處的云圖如圖10所示，懸掛連接處的應(yīng)力云圖如圖11所示。

圖9 肋板應(yīng)力較大處的等效應(yīng)力云圖Fig.9 The equivalent with strong floor stress nephogram

圖9中，最大應(yīng)力為285.43MPa。前面支架與犁的連接處應(yīng)力也較大，最大應(yīng)力為301MPa，如圖10所示；右前懸掛連接處的應(yīng)力較大，最大應(yīng)力為339.6MPa，如圖11所示。這3處的應(yīng)力都超過了Q235鋼的屈服極限，存在強度不足的問題，需要進行補強處理。

圖10 前機架應(yīng)力較大處等效應(yīng)力云圖Fig.10 The frame stress before a big equivalent stress nephogram

圖11 懸掛連接處的等效應(yīng)力云圖Fig.11 The equivalent stress nephogram of hanging junction

3.2 靈敏度分析

為了對懸掛式鏵式犁機架結(jié)構(gòu)進行優(yōu)化，首先確定設(shè)計變量，所以要進行靈敏度分析。這里選取了12個設(shè)計變量，均為懸掛式鏵式犁機架鋼結(jié)構(gòu)的厚度，如圖12所示。

圖12 懸掛式鏵式犁機架設(shè)計變量示意圖Fig.12 Hanging furrow plough frame design variable map

靈敏度分析所用的狀態(tài)變量為最小體積(減重需要)與圖9～圖11中所示危險區(qū)域的應(yīng)力(強度需要)，最后靈敏度分析如表2所示。

表2 靈敏度分析表

由表2可以看出：變量dv_11、dv_12、dv_8、dv_7對最小體積影響較大，而dv_1、dv_12、dv_6、dv_5對危險區(qū)域的應(yīng)力影響較大。綜合考慮，選取dv_1、dv_7、dv_8 、dv_11、dv_12、作為優(yōu)化時的設(shè)計變量。

4 機架結(jié)構(gòu)的優(yōu)化分析

對懸掛式鏵式犁機架梁基于等強度標(biāo)準(zhǔn)化設(shè)計方法進行結(jié)構(gòu)優(yōu)化，主要考慮靈敏度影響比較大的因素，在保證設(shè)計要求的情況下減少用材。所以，在對懸掛式鏵式犁進行結(jié)構(gòu)優(yōu)化時，設(shè)計變量為dv_1、dv_7、dv_8 、dv_11、dv_12，如圖13所示。

圖13 懸掛式鏵式犁機架優(yōu)化變量示意圖Fig.13 Hanging furrow plough frame diagram optimization variables

約束條件為：圖9～圖11中的應(yīng)力小于屈服極限235MPa。

對懸掛式鏵式犁機架進行基于等強度設(shè)計理論和目標(biāo)變量設(shè)計方法進行優(yōu)化，結(jié)果如表3所示。

表3 懸掛式鏵式犁機架尺寸優(yōu)化結(jié)果

對優(yōu)化后的懸掛式鏵式犁機架結(jié)構(gòu)利用Pro/E進行建模，之后在導(dǎo)入ANSYS中進行靜力學(xué)分析，為保持條件一致，施加的載荷及邊界條件與未優(yōu)化時的相同。經(jīng)靜力分析后，得到的位移云圖如圖14所示，應(yīng)變云圖如圖15所示，Von-mises應(yīng)力云圖如圖16所示。

圖14 改進后懸掛式鏵式犁機架等效位移云圖Fig.14 The improved hanging furrow plough frame equivalent displacement nephogram

圖15 改進后懸掛式鏵式犁機架等效應(yīng)變云圖Fig.15 The improved hanging plow frame strain nephogram

圖16 改進后懸掛式鏵式犁機架等效Von-mises應(yīng)力云圖Fig.16 The improved hanging furrow plough frame equivalent Von-mises stress nephogram

改進后，懸掛鏵式犁機架梁靜載荷作用下，應(yīng)力分布等效云圖中應(yīng)力分布均勻，機架梁整體應(yīng)力變化減小。由圖14與圖15可以看出，優(yōu)化后的位移與應(yīng)變都有較大的降低；由圖16可以看出，優(yōu)化后的結(jié)構(gòu)的最大應(yīng)力為192.1MPa，遠(yuǎn)小于屈服極限235MPa。這說明，優(yōu)化后的結(jié)構(gòu)是可靠的。

5 結(jié)論

通過對懸掛式鏵式犁機架進行有限元靜力學(xué)分析，得到了機架的應(yīng)力、應(yīng)變以及位移云圖。通過分析等效云圖可以看出：機架在承受靜載荷失效時，會從多個點處同時發(fā)生失效，最薄弱環(huán)節(jié)在應(yīng)力較大的部位，分別位于肋板、前面支架與犁的連接處，以及右前懸掛連接處。進行靈敏度分析有助于對現(xiàn)有的懸掛式鏵式犁機架進行等強度設(shè)計方法的結(jié)構(gòu)優(yōu)化。從優(yōu)化結(jié)果可以看出：改進后的機架能夠承載更大的載荷，改進后的位移與應(yīng)變都有較大的降低，應(yīng)力變化均勻。通過對懸掛式鏵式犁機架進行有限元分析，為機架的改進設(shè)計提供了理論依據(jù)，為實際應(yīng)用及標(biāo)準(zhǔn)化設(shè)計提供了參考。

[1] 鎮(zhèn)江農(nóng)業(yè)機械學(xué)院.農(nóng)業(yè)機械學(xué)[M].北京:中國農(nóng)業(yè)機械出版社，1981.

[2] 李耀明，孫朋朋，龐靖，等.聯(lián)合收獲機底盤機架有限元模態(tài)分析與試驗[J].農(nóng)業(yè)工程學(xué)報，2013，29 (3):38-42．

[3] 佚名.ANSYSWorkbench15.0—有限元分析從入門到精通[M].北京:機械工業(yè)出版社，2014.

[4] 呂廷，石秀東，張秋菊，等．基于 ANSYS 的破碎機機架模態(tài)分析[J]．機械設(shè)計與制造，2008(11):99-101．

[5] 凌桂龍，丁金濱，溫正．ANSYS Workbench13.0從入門到精通[M]．北京: 清華大學(xué)出版社，2012:137-138．

[6] 劉選偉, 金亮,王景立.基于ANSYSWorkbench深松機機架模態(tài)分析[J]．農(nóng)機化研究，2015,37(5):29-31,35.

[7] 哈爾濱工業(yè)大學(xué)理論力學(xué)教研室.理論力學(xué)(I)[M].北京:高等教育出版社，2002.

[8] 李超群,楊文彩，楊航，等.基于ANSYS三七播種機機架模態(tài)及諧響應(yīng)分析[J]．制造業(yè)自動化，2015(7)：77-79.

ID:1003-188X(2018)02-0061-EA

Suspension Bottom Plow Frame Finite Element Analysis and Optimization Design Based on Equal Strength

Han Chao, Yang Yougang, Zhang Shuwei

(College of Mechanical and Electronic Engineering, Northwest A&F University, Yangling 712100，China)

Abstract： Reasonably designed plow frame structures can ensure the tractor to drag a larger plow under the same traction effect, racks with less material and longer service life under the same operating conditions. In this study, ANSYS finite element analysis software was used to perform the static analysis and sensitivity analysis for hanging plow frame to be qualified through actual production test, and adjust the design by simulating the real environment of working condition to optimize the structure based on the method of equal strength . The results indicated that improved chassis could carry larger loads in the same structure after extending the load range on the basis of the original.

bottom plow; finite element; frame; structure optimization; standardized design

2016-11-22

陜西省專項資金項目(K332021313)；國家自然科學(xué)基金項目(41101201)

韓 超(1991-)，男，山東章丘人，碩士研究生，(E-mail)hc_runner@163.com。

楊有剛(1961-)，男，西安人，教授，碩士生導(dǎo)師，(E-mail)yangyougang@nwsuaf.edu.cn。

S222.12+1

A

1003-188X(2018)02-0061-05