基于有限元方法的凸齒鎮(zhèn)壓器運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)分析

張智泓，張廣凱，佟 金，Stephen Carr，賴慶輝，張兆國(guó)， 方秋容

(1.昆明理工大學(xué)現(xiàn)代農(nóng)業(yè)工程學(xué)院，云南昆明 650500；2.吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院，吉林長(zhǎng)春 130025；3.吉林大學(xué)工程仿生教育部重點(diǎn)實(shí)驗(yàn)室，吉林長(zhǎng)春 130025；4.International Soil and Water Renewables，LLC，Indiana Salem 47167 US)

?

基于有限元方法的凸齒鎮(zhèn)壓器運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)分析

張智泓1，張廣凱1，佟 金2，3，Stephen Carr4，賴慶輝1，張兆國(guó)1， 方秋容1

(1.昆明理工大學(xué)現(xiàn)代農(nóng)業(yè)工程學(xué)院，云南昆明 650500；2.吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院，吉林長(zhǎng)春 130025；3.吉林大學(xué)工程仿生教育部重點(diǎn)實(shí)驗(yàn)室，吉林長(zhǎng)春 130025；4.International Soil and Water Renewables，LLC，Indiana Salem 47167 US)

摘要為探索凸齒鎮(zhèn)壓器在作業(yè)過(guò)程中的運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)特性，基于驗(yàn)證有效性后的有限元模型，求解了凸齒鎮(zhèn)壓器外緣點(diǎn)運(yùn)動(dòng)軌跡和速度以及加速度隨時(shí)間的變化規(guī)律曲線。結(jié)果表明：輪緣點(diǎn)所在凸齒沖擊土壤的行為并未在輪緣點(diǎn)豎直方向速度最大時(shí)完成；凸齒鎮(zhèn)壓器輪緣外一點(diǎn)的絕對(duì)運(yùn)動(dòng)軌跡呈波動(dòng)狀態(tài)的余擺線形式；凸齒鎮(zhèn)壓器輪緣點(diǎn)在接地位置有與前進(jìn)方向相反的相對(duì)滑動(dòng)。

關(guān)鍵詞凸齒鎮(zhèn)壓器；土壤表面微形貌加工；有限元方法；運(yùn)動(dòng)軌跡求解；運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)分析

土壤表面微形貌加工是近年來(lái)在干旱與半干旱耕作區(qū)得到快速發(fā)展的一種整地方式。在地面機(jī)械領(lǐng)域，研究橫軸類旋轉(zhuǎn)工作觸土部件的工作機(jī)理時(shí)，探索工作部件特征點(diǎn)的運(yùn)動(dòng)軌跡曲線是分析機(jī)具正常工作的條件及其量化指標(biāo)的重要研究?jī)?nèi)容[1-3]。目前，求解機(jī)具旋轉(zhuǎn)工作部件運(yùn)動(dòng)軌跡主要依靠經(jīng)典力學(xué)的基礎(chǔ)理論推導(dǎo)，然而使用基礎(chǔ)理論求解凸齒鎮(zhèn)壓器輪緣點(diǎn)運(yùn)動(dòng)軌跡具有如下局限性：分析過(guò)程中需要做大量簡(jiǎn)化，與實(shí)際作業(yè)條件有一定的偏離；與地面相互作用時(shí)，受到垂直載荷、凸齒形狀、牽引速度等參數(shù)的影響；土壤流動(dòng)變形特征具有復(fù)雜本質(zhì)，導(dǎo)致一些工作參數(shù)難以測(cè)量。隨著計(jì)算機(jī)計(jì)算能力的逐步提高，借助于計(jì)算機(jī)應(yīng)用有限元等數(shù)值方法求解機(jī)械對(duì)土壤加工的典型工藝問(wèn)題增多，有限元分析方法成為分析土壤與農(nóng)機(jī)觸土部件相互作用規(guī)律的有效工具。為了揭示凸齒鎮(zhèn)壓器的輪緣點(diǎn)運(yùn)動(dòng)特征和運(yùn)動(dòng)學(xué)特性，筆者基于有限元方法，對(duì)凸齒鎮(zhèn)壓器工作過(guò)程中輪緣點(diǎn)的軌跡及其運(yùn)動(dòng)學(xué)特征進(jìn)行求解和分析，旨在為提高凸齒鎮(zhèn)壓器的工作性能和作業(yè)質(zhì)量提供理論基礎(chǔ)。

1材料與方法

1.1土壤表面微形貌加工原理土壤表面微形貌加工原理如圖1所示，通過(guò)在田間土壤表面加工出若干規(guī)則微坑陣列以改變地表微觀地形地貌，從而攔蓄降水產(chǎn)生的田間地表徑流[4]。通過(guò)該整地措施，地表徑流得以匯集并滲入土壤內(nèi)部，不僅有效地控制過(guò)量降水造成的水土流失，還把雨水轉(zhuǎn)化為土壤水，使水分富集于地表土壤以供作物生長(zhǎng)利用，從而趨利避害，在降低過(guò)量降水對(duì)農(nóng)業(yè)生產(chǎn)的負(fù)面影響的同時(shí)集蓄降水，使其對(duì)作物生長(zhǎng)產(chǎn)生積極作用，最大限度地提高天然降水的利用效率，達(dá)到控制水土流失和緩解干旱缺水的雙重目標(biāo)[5-7]。

凸齒鎮(zhèn)壓器是實(shí)現(xiàn)土壤表面微形貌加工作業(yè)的新型特殊滾動(dòng)觸土部件。該機(jī)具圓柱面上均勻分布凸齒，凸齒形狀為4個(gè)端面為斜角的棱柱，凸齒鎮(zhèn)壓器機(jī)具作業(yè)時(shí)，由牽引機(jī)具以一定速度牽引并在田間土壤表面滾動(dòng)，滾動(dòng)時(shí)凸齒擠壓地表土壤，使地表土壤流動(dòng)和變形，在田間地表加工出幾何形狀規(guī)則的微坑形貌[8-9]。在凸齒鎮(zhèn)壓器加工地表微形貌的過(guò)程中，其非圓形輪廓在滾動(dòng)時(shí)重心呈周期性的交替升高和下降，凸齒鎮(zhèn)壓器的勢(shì)能和動(dòng)能周期性轉(zhuǎn)化為沖擊能作用于地面，凸齒鎮(zhèn)壓器的輪緣點(diǎn)運(yùn)動(dòng)軌跡、速度、加速度等因素影響著工作效率和作業(yè)質(zhì)量。

注：a為土壤表面微形貌加工過(guò)程；b為微坑蓄水示意。Note：a.Process of micro-topographical preparation；b.Diagram of water storage in micro-basin.圖1　土壤表面微形貌加工原理Fig.1　Principles of micro-topographical preparation

而在農(nóng)業(yè)機(jī)械領(lǐng)域，旋轉(zhuǎn)工作觸土部件可歸納為驅(qū)動(dòng)式和牽引式，由于工作原理不同，其工作部件的運(yùn)動(dòng)情況也不同。對(duì)于非驅(qū)動(dòng)牽引式旋轉(zhuǎn)工作部件，如旋耕機(jī)，工作時(shí)刀片一邊繞軸正向旋轉(zhuǎn)，一邊隨機(jī)組做直線運(yùn)動(dòng)，因此，刀片的絕對(duì)運(yùn)動(dòng)軌跡是一條由旋轉(zhuǎn)運(yùn)動(dòng)與直線運(yùn)動(dòng)合成的數(shù)學(xué)擺線。數(shù)學(xué)擺線形式由機(jī)組前進(jìn)速度和機(jī)具旋轉(zhuǎn)角速度兩個(gè)參數(shù)同時(shí)決定，由于二者之間的數(shù)值及其組合不同，其合成后的擺線形狀存在較大差異，擺線形式可分為短擺線、滾擺線、余擺線，并且對(duì)旋耕機(jī)最終的工作結(jié)果產(chǎn)生不同影響。對(duì)于非驅(qū)動(dòng)牽引式旋轉(zhuǎn)工作部件，如傳統(tǒng)剛性輪，在地面做無(wú)滑動(dòng)滾動(dòng)時(shí)其輪緣外一點(diǎn)的運(yùn)動(dòng)軌跡是標(biāo)準(zhǔn)滾擺線形式[10-11]。然而，凸齒鎮(zhèn)壓器與其他牽引式典型圓形工作輪[12-14]有本質(zhì)區(qū)別，凸齒鎮(zhèn)壓器為非圓形工作輪，運(yùn)動(dòng)時(shí)其質(zhì)心呈周期性的升高和下降，故作業(yè)過(guò)程中對(duì)土壤產(chǎn)生沖擊作用，具有沖擊壓實(shí)的效果。另外，與剛性輪在地表做無(wú)滑動(dòng)滾動(dòng)不同的是，凸齒鎮(zhèn)壓器在土壤表面工作時(shí)，工作情況具有復(fù)雜特性，一方面，土壤為彈塑性材料，凸齒鎮(zhèn)壓器在土壤表面工作的過(guò)程中土壤會(huì)發(fā)生變形；另一方面，凸齒鎮(zhèn)壓器與土壤的接觸界面存在滑動(dòng)摩擦。1.2凸齒鎮(zhèn)壓器工況設(shè)置凸齒鎮(zhèn)壓器與土壤相互作用時(shí)的工況如圖2所示。相對(duì)運(yùn)動(dòng)：凸齒鎮(zhèn)壓器外緣的一點(diǎn)M相對(duì)于凸齒鎮(zhèn)壓器的中心RP做旋轉(zhuǎn)運(yùn)；牽連運(yùn)動(dòng)：凸齒鎮(zhèn)壓器中心RP被牽引機(jī)構(gòu)牽引，在水平方向向前運(yùn)動(dòng)；在加工微坑過(guò)程中，凸齒鎮(zhèn)壓器在豎直方向載荷與水平方向牽引力的作用下，不斷嵌入土壤并被抬起，因此凸齒鎮(zhèn)壓器中心RP在豎直方向上做往復(fù)運(yùn)動(dòng)。

圖2　凸齒鎮(zhèn)壓器工況Fig.2　Working condition of toothed wheel

1.3有限元模型尺寸和邊界條件凸齒鎮(zhèn)壓器與土壤作用模型如圖3所示。L為模擬土槽的長(zhǎng)度，尺寸設(shè)置為2 500 mm；W為模擬土槽的寬度，尺寸設(shè)置為200 mm；H為模擬土槽的高度，尺寸設(shè)置為250 mm。凸齒鎮(zhèn)壓器以1 000 mm/s的速度沿X軸向前運(yùn)動(dòng)2 000 mm。為了使凸齒鎮(zhèn)壓器各單元的牽引力集中計(jì)算輸出，設(shè)定凸齒鎮(zhèn)壓器滾動(dòng)輪旋轉(zhuǎn)中心為整個(gè)凸齒鎮(zhèn)壓器的剛性參考點(diǎn)(RP)。將土體模型底面設(shè)定為全局約束，限制底面的全部6個(gè)自由度，將通過(guò)原點(diǎn)的YZ平面設(shè)為對(duì)稱面，土體模型的上表面不做約束。參考點(diǎn)RP沿水平方向以固定速度運(yùn)動(dòng)，時(shí)沿豎直方向自由移動(dòng)，輪緣點(diǎn)M繞參考點(diǎn)RP做圓周運(yùn)動(dòng)，同時(shí)做以上幾種運(yùn)動(dòng)的合成運(yùn)動(dòng)。由于該模型為幾何對(duì)稱型，為提高計(jì)算效率，采用沿對(duì)稱面分割的半個(gè)模型進(jìn)行模擬，但輸出結(jié)果為完整模型的結(jié)果。

圖3　土壤模型Fig.3　Size of soil bin model

1.4凸齒鎮(zhèn)壓器的尺寸凸齒鎮(zhèn)壓器的尺寸(圖4)為目前被廣泛采用的凸齒鎮(zhèn)壓器設(shè)計(jì)尺寸[4]，其中D和d分別為凸齒鎮(zhèn)壓器滾動(dòng)輪的外徑和內(nèi)徑，尺寸分別為340和300 mm，U為凸齒鎮(zhèn)壓器的滾動(dòng)輪寬度，尺寸為115 mm。凸齒鎮(zhèn)壓器上沿圓周均勻分布有6個(gè)凸齒，每個(gè)齒為三角棱形狀，凸齒寬度z為90 mm，凸齒高度y為70 mm，凸齒底部長(zhǎng)度b為150 mm，凸齒頂部長(zhǎng)度t為100 mm。

圖4　凸齒鎮(zhèn)壓器模型Fig.4　Size of toothed wheel model

1.5有限元模型網(wǎng)格在有限元模型的建立過(guò)程中，有限元網(wǎng)格劃分的質(zhì)量直接影響計(jì)算的精度和速度。在進(jìn)行網(wǎng)格劃分時(shí)，首先確定單元類型。模擬土槽模型的網(wǎng)格采用八節(jié)點(diǎn)縮減積分實(shí)體單元C3D8R，該類型的單元常被用于連續(xù)材料的三維應(yīng)力應(yīng)變分析，土壤模型單元尺寸設(shè)置為8 mm，節(jié)點(diǎn)總數(shù)為261 248，單元總數(shù)為242 575。凸齒鎮(zhèn)壓器模型作為離散剛體導(dǎo)入有限元軟件，其網(wǎng)格采用三維四面體單元C3D4，單元為9 mm，節(jié)點(diǎn)總數(shù)為5 781，單元總數(shù)為25 006。

1.6模型運(yùn)行計(jì)算機(jī)環(huán)境有限元模型的計(jì)算求解在吉林大學(xué)工程仿生教育部重點(diǎn)實(shí)驗(yàn)室的工作站進(jìn)行。工作站型號(hào)：戴爾Precision T7500(CPU型號(hào)：E5645，CPU主頻：2.4 G，內(nèi)存容量：8*DDR3)。

1.7有限元模型的有效性驗(yàn)證基于前期研究結(jié)果[8]，為了驗(yàn)證模型的準(zhǔn)確性，在吉林大學(xué)工程仿生教育部重點(diǎn)實(shí)驗(yàn)室土槽實(shí)驗(yàn)室進(jìn)行驗(yàn)證實(shí)驗(yàn)。選擇牽引力和土壤流動(dòng)變形結(jié)果作為驗(yàn)證參數(shù)，基于室內(nèi)土槽測(cè)試系統(tǒng)搭建牽引實(shí)驗(yàn)平臺(tái)(圖5)，在載荷為250 N，牽引機(jī)構(gòu)運(yùn)行速度為1 m/s的作業(yè)條件下，考察實(shí)際實(shí)驗(yàn)與有限元分析的誤差。通過(guò)比較有限元與土槽實(shí)驗(yàn)求解的牽引力與時(shí)間的關(guān)系曲線，發(fā)現(xiàn)2條曲線變化規(guī)律一致；比較牽引力平均值，計(jì)算得到相對(duì)誤差為3.40%，證實(shí)有限元分析結(jié)果達(dá)到所需精度；同時(shí)，將有限元方法求解得到的微坑形貌與土槽實(shí)驗(yàn)結(jié)果比較，結(jié)果表明有限元方法預(yù)測(cè)的微坑外部和內(nèi)部的形貌特征與實(shí)測(cè)實(shí)驗(yàn)一致。因此，筆者所建立的有限元模型可有效地模擬凸齒鎮(zhèn)壓器與土壤間的相互作用。

圖5　牽引實(shí)驗(yàn)平臺(tái)Fig.5　Test platform of draft force

2結(jié)果與分析

2.1場(chǎng)變量輸出結(jié)果模型連續(xù)求解時(shí)間合計(jì)260 h，運(yùn)行完成后，場(chǎng)變量輸出結(jié)果如圖6所示。土壤模型網(wǎng)格規(guī)整，未出現(xiàn)網(wǎng)格單元過(guò)渡扭曲現(xiàn)象。

2.2輪緣點(diǎn)運(yùn)動(dòng)軌跡曲線求解所得點(diǎn)M(圖6)的絕對(duì)運(yùn)動(dòng)軌跡如圖7所示。與剛性輪在地表做純滾動(dòng)的標(biāo)準(zhǔn)旋輪線運(yùn)動(dòng)軌跡不同，凸齒鎮(zhèn)壓器在土壤表面運(yùn)動(dòng)時(shí)，由于凸齒鎮(zhèn)壓器圓周上的凸齒依次壓入土壤和抬離地面，凸齒外緣點(diǎn)M的運(yùn)動(dòng)曲線并非光滑曲線，而是呈規(guī)律的波動(dòng)狀態(tài)。施加豎直向下載荷后，點(diǎn)M在0~40 mm位置有一個(gè)向下的位移，對(duì)點(diǎn)RP施加水平速度后，隨著凸齒鎮(zhèn)壓器的向前運(yùn)動(dòng)，凸齒被抬起，點(diǎn)M向上移動(dòng)，在800 mm位置處運(yùn)動(dòng)至最高點(diǎn)，完成提升過(guò)程；此后開(kāi)始向下移動(dòng)，在1 400 mm位置處，點(diǎn)M的豎直方向位移已為負(fù)值，此時(shí)點(diǎn)M位置的凸齒已開(kāi)始貫入微坑；點(diǎn)M在1 450 mm位置運(yùn)動(dòng)至最低點(diǎn)并到達(dá)微坑的底部，完成沖壓過(guò)程。在圖6的框中放大區(qū)域位置，點(diǎn)M有一定的向后移動(dòng)距離，在此處凸齒壓入土壤，然后隨著凸齒鎮(zhèn)壓器繼續(xù)向前運(yùn)動(dòng)凸齒再次被抬離，進(jìn)入下一個(gè)工作周期。

圖6　場(chǎng)變量輸出結(jié)果Fig.6　 Result of field ouput

圖7　凸齒鎮(zhèn)壓器輪緣點(diǎn)的運(yùn)動(dòng)軌跡Fig.7　Motion trail of toothed wheel outer edge point

為進(jìn)一步揭示凸齒鎮(zhèn)壓器的運(yùn)動(dòng)軌跡規(guī)律，將點(diǎn)M的合位移分解為沿水平和豎直方向的分位移。如圖8所示，凸齒鎮(zhèn)壓器運(yùn)動(dòng)開(kāi)始后，點(diǎn)M在豎直方向向上移動(dòng)的過(guò)程中，水平方向位移也在增加；點(diǎn)M在0.9 s時(shí)到達(dá)最高點(diǎn)，此后迅速向下移動(dòng)；在1.5 s時(shí)，點(diǎn)M豎直方向位移值為負(fù)，此時(shí)點(diǎn)M已開(kāi)始貫入微坑，并在1.6 s時(shí)運(yùn)動(dòng)至最低點(diǎn)；點(diǎn)M運(yùn)動(dòng)到最低點(diǎn)后，水平方向位移曲線趨于平緩并有一定程度的下降，此時(shí)凸齒壓入微坑最底端，說(shuō)明凸齒與微坑底部在水平方向有一定的相對(duì)滑動(dòng)；在1.7 s時(shí)，點(diǎn)M豎直方向位移轉(zhuǎn)換為正值，此時(shí)點(diǎn)M所在凸齒已從微坑中抬離。

2.3輪緣點(diǎn)速度隨時(shí)間變化曲線的求解分別對(duì)點(diǎn)M在水平方向和豎直方向的速度分量隨時(shí)間的變化曲線進(jìn)行求解。從圖9可以看出，點(diǎn)M在0.9 s達(dá)到最高點(diǎn)之前，其所在的凸齒屬于提升過(guò)程，水平方向速度分量不斷增加，而豎直方向速度分量先增加后減小，并在0.9 s位置為0；點(diǎn)M在0.9 s達(dá)到最高點(diǎn)之后，其所在的凸齒進(jìn)入沖擊過(guò)程，水平方向速度分量逐漸減少，豎直方向分量沿相反的方向逐漸增大，即點(diǎn)M在向下做加速運(yùn)動(dòng)，在1.3 s位置點(diǎn)M豎直向下的速度達(dá)最大值。然而，點(diǎn)M所在凸齒沖擊土壤的行為并未在其速度最大時(shí)完成，這是由于其所在凸齒開(kāi)始沖擊土壤的同時(shí)也是前一個(gè)凸齒從土壤中抬起的過(guò)程，該過(guò)程阻礙了點(diǎn)M所在凸齒的向下運(yùn)動(dòng)，從而減小了其豎直向下的速度。而點(diǎn)M在1.5 s和土壤接觸時(shí)仍具有較高的豎直向下速度，因此沖擊作用明顯。沖擊發(fā)生后，點(diǎn)M在豎直方向的速度變化有一定減緩，此時(shí)點(diǎn)M所在凸齒貫入微坑，在水平方向上的速度為最低點(diǎn)，隨后出現(xiàn)正負(fù)交替，同樣說(shuō)明此時(shí)點(diǎn)M在水平方向上與土壤有相對(duì)滑動(dòng)。2.4輪緣點(diǎn)加速度隨時(shí)間變化曲線的求解為更加直觀地觀察水平方向和豎直方向加速度分量隨時(shí)間的變化情況，分別對(duì)點(diǎn)M在水平方向和豎直方向的加速度分量隨時(shí)間的變化曲線進(jìn)行求解。從圖10可以看出，沖擊發(fā)生的1.4 s前的瞬間，點(diǎn)M豎直方向加速度為正，說(shuō)明此前點(diǎn)M豎直向下的速度在逐漸減小；沖擊發(fā)生后，點(diǎn)M豎直方向加速度為0，說(shuō)明此時(shí)點(diǎn)M有一段勻速運(yùn)動(dòng)過(guò)程；沖擊完成后，點(diǎn)M加速度迅速增大，說(shuō)明此時(shí)點(diǎn)M貫入微坑過(guò)程中向下運(yùn)動(dòng)的速度迅速減?。蛔罱K，點(diǎn)M離開(kāi)微坑時(shí)向上做加速運(yùn)動(dòng)。

注：a為水平方向的運(yùn)動(dòng)軌跡；b為豎直方向的運(yùn)動(dòng)軌跡。Note：a.Motion trail at horizontal direction；b.Motion trail at vertical direction.圖8　凸齒鎮(zhèn)壓器輪緣點(diǎn)在不同方向的運(yùn)動(dòng)軌跡Fig.8　Motion trail of toothed wheel outer edge point at different directions

注：a為水平方向的運(yùn)動(dòng)速度分量；b為豎直方向的運(yùn)動(dòng)速度分量。Note：a.Velocity at horizontal direction； b.Velocity at vertical direction.圖9　凸齒鎮(zhèn)壓器輪緣點(diǎn)在不同方向的運(yùn)動(dòng)速度分解Fig.9　Velocity of toothed wheel outer edge point at different directions

注：a為水平方向的運(yùn)動(dòng)速度分量；b為豎直方向的運(yùn)動(dòng)速度分量。Note：a.Velocity at horizontal direction；b.Velocity at vertical direction.圖10　凸齒鎮(zhèn)壓器輪緣點(diǎn)在不同方向的運(yùn)動(dòng)加速度分解Fig.10　Velocity of toothed wheel outer edge point at different directions

3結(jié)論與討論

筆者使用有限元方法求解了凸齒鎮(zhèn)壓器輪緣上一點(diǎn)的運(yùn)動(dòng)軌跡，以及該點(diǎn)速度、加速度不同方向的分量隨時(shí)間的變化情況，結(jié)果表明：①?gòu)慕^對(duì)運(yùn)動(dòng)軌跡隨時(shí)間變化曲線可以看出，凸齒鎮(zhèn)壓器輪緣外一點(diǎn)的絕對(duì)運(yùn)動(dòng)軌跡為非光滑曲線，表現(xiàn)為呈波動(dòng)狀態(tài)的余擺線形式；②從輪緣點(diǎn)位移與速度隨時(shí)間變化曲線可以看出，輪緣點(diǎn)運(yùn)動(dòng)到最低點(diǎn)時(shí)，凸齒壓入微坑最底端，水平方向位移曲線趨于平緩并有一定程度的下降，同時(shí)觀察到輪緣點(diǎn)在水平方向上的速度矢量出現(xiàn)正負(fù)交替，證實(shí)凸齒鎮(zhèn)壓器輪緣點(diǎn)在接地位置有與前進(jìn)方向相反的相對(duì)滑動(dòng)；③從速度與加速度隨時(shí)間變化曲線可以看出，沖擊發(fā)生前的瞬間，輪緣點(diǎn)豎直方向加速度為正、方向向上，說(shuō)明此前輪緣點(diǎn)豎直向下的速度逐漸減小，輪緣點(diǎn)所在凸齒沖擊土壤的行為并未在輪緣點(diǎn)豎直方向速度最大時(shí)完成，這是由于輪緣點(diǎn)所在凸齒開(kāi)始沖擊土壤的同時(shí)也是前一個(gè)凸齒從土壤抬離的過(guò)程，該過(guò)程阻礙了輪緣點(diǎn)所在凸齒的向下運(yùn)動(dòng)，從而減小了輪緣點(diǎn)豎直向下的速度。

該研究揭示的基本信息和科學(xué)規(guī)律為今后凸齒鎮(zhèn)壓器的研究提出如下思路：一方面，充分利用凸齒鎮(zhèn)壓器輪緣和土壤的相對(duì)滑動(dòng)現(xiàn)象改變凸齒刃口結(jié)構(gòu)，如將平滑刃口改為鋸齒結(jié)構(gòu)刃口，從而將滑切方式改為鋸切方式，以減小凸齒貫入土壤的阻力，并增加所加工微坑的深度；另一方面，改變凸齒沿輪緣排列方式或改變凸齒外緣輪廓幾何結(jié)構(gòu)，以使凸齒沖擊土壤的行為在輪緣點(diǎn)速度最大時(shí)完成，從而提高所加工微坑的質(zhì)量和效率。

參考文獻(xiàn)

[1] 何岳平，陳青春，何瑞銀，等.以穴寬為優(yōu)化目標(biāo)的大蒜栽植機(jī)關(guān)鍵機(jī)構(gòu)運(yùn)動(dòng)仿真及優(yōu)化設(shè)計(jì)[J].安徽農(nóng)業(yè)科學(xué)，2009，37(36)：8161-8164，8179.

[2] 張華標(biāo)，朱梅，朱德泉，等.高速水稻插秧機(jī)分插機(jī)構(gòu)振動(dòng)特性研究 [J].安徽農(nóng)業(yè)科學(xué)，2012，40(18)：9906-9907，9919.

[3] 趙志國(guó).吊籃-鴨嘴式移栽機(jī)栽植器運(yùn)動(dòng)軌跡分析[J].安徽農(nóng)業(yè)科學(xué)，2011，39(19)：11932，11935.

[4] PATRICK C，KECHAVARZI C，JAMES I T，et al.Developing reservoir tillage technology for semi-arid environments [J].Soil use and management，2007，23(2)：185-191.

[5] SALEM H M，VALERO C，MUOZ M，et al.Effect of integrated reservoir tillage for in-situ rainwater harvesting and other tillage practices on soil physical properties [J].Soil and tillage research，2015，151：50-60.

[6] SALEM H M，VALERO C，MUOZ M，et al.Short-term effects of four tillage practices on soil physical properties，soil water potential，and maize yield [J].Geoderma，2015，237/238：60-70.

[7] SALEM H M，VALERO C，MUOZ M，et al.Effect of reservoir tillage on rainwater harvesting and soil erosion control under a developed rainfall simulator[J].Catena，2014，113：353-362.

[8] 佟金 張智泓，陳東輝，等.凸齒鎮(zhèn)壓器與土壤相互作用的三維動(dòng)態(tài)有限元分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2014，30(10)：48-58，293.

[9] 張智泓，佟金，陳東輝，等.不同材質(zhì)仿生凸齒鎮(zhèn)壓器滾動(dòng)件的模態(tài)分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28(13)：8-15.

[10] KIM K U，SHIN B S.Modeling motion resistance of rigid wheels[J].Journal of Terramechanics，1985，22(4)：225-236.

[11] 高海波.用旋輪線方程求輪緣上點(diǎn)的速度和加速度[J].鄂州大學(xué)學(xué)報(bào)，1996(3)：36-37，46.

[12] WONG J Y.Behaviour of soil beneath rigid wheels [J].Journal of agricultural engineering research，1967，12(4)：257-269.

[13] NAKASHIMA H，KOBAYASHI T.Effects of gravity on rigid rover wheel sinkage and motion resistance assessed using two-dimensional discrete element method [J].Journal of terramechanics，2014，53：37-45.

Kinematics and Dynamics Analysis of Toothed Wheel Based on Finite Element Method

ZHANG Zhi-hong1, ZHANG Guang-kai1, TONG Jin2,3et al(1. School of Modern Agricultural Engineering, Kunming University of Science and Technology, Kunming, Yunnan 650500; 2. School of Biological and Agricultural Engineering, Jilin University, Changchun, Jilin 130025; 3. The Key Laboratory of Bionic Engineering, Jilin University, Changchun, Jilin 130025)

AbstractIn order to investigate the kinematics and dynamics characteristics of operating toothed wheel, we adopted the validated finite element model, solved the toothed wheel rim point motion trajectory, and obtained the change rule curves of velocity and acceleration velocity with the time. It was found that when toothed wheel rim point hit the soil, the speed at vertical direction did not reach its maximum; motion trajectory of toothed wheel rim point presented to be a fluctuating trochoidal curve. At the bottom of micro-basin, toothed wheel rim point motion had relative sliding in the opposite direction of movement.

Key wordsToothed wheel; Micro-topographical preparation; Finite element method; Motion trajectory solution; Kinematics and dynamics analysis

收稿日期2015-12-21

作者簡(jiǎn)介張智泓(1986- )，男，白族，云南昆明人，講師，博士，從事農(nóng)業(yè)機(jī)械化工程和農(nóng)業(yè)工程仿生研究。

基金項(xiàng)目國(guó)家自然科學(xué)基金面上項(xiàng)目(51075185)；云南省科技計(jì)劃項(xiàng)目青年項(xiàng)目(2015FD011)；云南省教育廳科學(xué)研究基金項(xiàng)目(2015Y079)；昆明理工大學(xué)引進(jìn)人才科研啟動(dòng)基金項(xiàng)目(14118940)；昆明理工大學(xué)學(xué)生課外學(xué)術(shù)科技創(chuàng)新基金(2015BA132)。

中圖分類號(hào)s 222.29

文獻(xiàn)標(biāo)識(shí)碼A

文章編號(hào)0517-6611(2016)02-325-05