王恩民，李立君，湯剛車(chē)，歐陽(yáng)益斌

(中南林業(yè)科技大學(xué) 機(jī)電工程學(xué)院，長(zhǎng)沙 410004)

0 引言

林業(yè)機(jī)械的發(fā)展是機(jī)械化水平提高的重要保證，尤以深松機(jī)為重。目前，保護(hù)性耕作受到越來(lái)越多的重視[1-3]，土壤深松技術(shù)作為保護(hù)性耕作關(guān)鍵技術(shù)之一更是一種新型的土壤耕作方法。土壤深松技術(shù)是指在土壤不被翻轉(zhuǎn)的情況下進(jìn)行土壤疏松，對(duì)土壤中固、液、氣3種成分比例進(jìn)行調(diào)節(jié)，考慮到對(duì)林地作物情況影響，土壤耕深在200～300mm之間有利于根系生長(zhǎng)，促進(jìn)作物增產(chǎn)增收，對(duì)于保水保墑也發(fā)揮著重要作用[4]。

目前，對(duì)深松機(jī)研究多集中在旋耕刀或深松鏟的研究，以及基于這兩種刀具的受力分析及振動(dòng)減阻原理研究，研究取得一系列成果[5-7]；但是，耕深不穩(wěn)定、阻力大、功耗大、深松后土壤硬度及平整度達(dá)不到要求等一系列問(wèn)題仍沒(méi)有得到很好的解決。針對(duì)這一現(xiàn)狀，設(shè)計(jì)了一種基于旋轉(zhuǎn)切削原理的立式深松刀具，確定了刀具工作參數(shù)，并采用SPH法[8-11]建立刀具切削土壤模型，開(kāi)展試驗(yàn)驗(yàn)證仿真可信度，為后續(xù)優(yōu)化提供技術(shù)依據(jù)。

1 深松工作參數(shù)確定

1.1 深松刀結(jié)構(gòu)設(shè)計(jì)

深松刀由法蘭，加強(qiáng)肋、深松刀軸及深松刀片組成，結(jié)構(gòu)如圖1所示。

1.法蘭 2.加強(qiáng)肋 3.深松刀軸 4.深松刀片圖1 深松刀結(jié)構(gòu)圖Fig.1 Deep loosening tool structure

深松刀作業(yè)時(shí)，依靠刀軸帶動(dòng)刀片，刀片在刀軸上呈等螺距分布，工作頭升降使刀片可以對(duì)不同深度的土壤進(jìn)行深松，且土壤不被翻轉(zhuǎn)，滿(mǎn)足保護(hù)性耕作的要求。

切削過(guò)程中，深松刀切削阻力與切削角正切成正比[12]，但一定的切削角才能使土壤微粒向上運(yùn)輸，不至于堵塞造成阻力變大，同時(shí)保證切土效果。切削角取值范圍為15°～25°。

1.2 深松刀運(yùn)動(dòng)分析

深松過(guò)程中，假設(shè)底盤(pán)以固定速度前進(jìn)，深松刀軸以固定速度旋轉(zhuǎn)，深松刀為勻速直線(xiàn)運(yùn)動(dòng)與勻速旋轉(zhuǎn)運(yùn)動(dòng)的復(fù)合。其運(yùn)動(dòng)軌跡可用余擺線(xiàn)方程表示，即

式中R—深松刀旋轉(zhuǎn)半徑(m)；

ω—深松刀旋轉(zhuǎn)速度(rad/s)；

v—深松刀前進(jìn)速度(m/s)。

1.3 深松過(guò)程最低旋轉(zhuǎn)速度確定

深松刀工作過(guò)程，為了使土壤微粒能隨著深松刀旋轉(zhuǎn)向上運(yùn)輸，深松刀旋轉(zhuǎn)速度必須大于一個(gè)臨界值。

1.3.1 土壤微粒速度分析

土壤微粒速度分析如圖2所示。

(a) xoz平面內(nèi)土壤微粒速度分析圖

(b) xoy平面內(nèi)土壤微粒速度分析圖圖2 土壤微粒速度分析圖Fig.2 Analysis of soil particle velocity

設(shè)深松刀角速度為ω，前進(jìn)速度為v0，對(duì)土壤微粒在刀邊緣進(jìn)行速度分布。

深松刀在旋轉(zhuǎn)時(shí)，土壤微粒也在以ve的速度隨刀旋轉(zhuǎn)的同時(shí)，以vr的速度向上運(yùn)動(dòng)，va為vr與ve的合速度，與水平面夾角為β，α為刀片的傾角,則

(1)

水平方向上

vt=ve-vrx

(2)

vt=va·cosβ

(3)

ve=ω·R

(4)

vrx=vz·cotα

(5)

vz=va·sinβ

(6)

因?yàn)関t=ωa·R=va·cosβ，則

(7)

將式(3)～式(7)代入式(2)中得

(8)

式中ω—深松刀旋轉(zhuǎn)速度；

α—深松刀片傾角；

β—土壤微粒絕對(duì)速度與水平面的夾角；

ωa—土壤微粒絕對(duì)角速度。

1.3.2 土壤微粒受力分析

土壤微粒受力分析如圖3所示。

由式(3)得

dF1=N1·f1

(9)

dF2=N2·f2

(10)

式中dF1—深松刀片與土壤微粒間摩擦力，與vr方向相反；

dF2—溝壁與土壤微粒間摩擦力；

N1—刀片對(duì)土壤微粒支撐力；

N2—溝壁對(duì)土壤微粒支撐力；

f1—土壤與螺旋面摩擦因數(shù)；

f2—土壤內(nèi)摩擦因數(shù)。

據(jù)圖3可知：在N1的方向上有

N1=dF2·cos(90°-α-β)+dm·(g+az)cosα

(11)

式中dm—土壤微粒的質(zhì)量；

az—土壤微粒的徑向加速度。

在垂直FN1方向上有

(12)

(13)

式中FN1—dF1與N1的合力；

φ1—FN1與N1的夾角；

Fe—土壤微粒所受離心力；

N2—溝壁對(duì)土壤微粒的支持力。

(a) xoz平面內(nèi)土壤微粒受力分析圖

(b) xoy平面內(nèi)土壤微粒受力分析圖圖3 土壤微粒受力分析圖Fig.3 Analysis of soil particles by force

將式(10)、式(13)代入式(12)中得

代入式(8)得

令β=0得到深松刀需滿(mǎn)足的最低角速度為

式中ω0—滿(mǎn)足土粒升運(yùn)最低角速度；

g—重力加速度；

R—深松刀片外徑；

α—深松刀片傾角；

φ1—土壤質(zhì)點(diǎn)與螺旋面摩擦角。

計(jì)算得：ω0=33.9rad/s。

1.4 深松前進(jìn)速度確定

配套深松機(jī)前進(jìn)速度如表1所示。

表1 前進(jìn)速度表Table 1 Forward speedTable m/s

根據(jù)深松作業(yè)效率要求得深松機(jī)前進(jìn)速度需大于0.5m/s。

2 刀具切削土壤過(guò)程模擬

2.1 SPH算法基本原理

刀具切削土壤屬于高速切削，采用有限元網(wǎng)格法進(jìn)行仿真易造成網(wǎng)格畸變導(dǎo)致結(jié)果不收斂。SPH算法是一種無(wú)網(wǎng)格算法，將土壤離散成一系列具有相互作用且?guī)в匈|(zhì)量和速度等物理量的點(diǎn)集合，通過(guò)粒子模擬流體運(yùn)動(dòng)規(guī)律[13]。其核心是插值理論，要求解各個(gè)質(zhì)點(diǎn)在任意時(shí)刻能量和速度，需引入核函數(shù)；場(chǎng)函數(shù)經(jīng)過(guò)核函數(shù)“光滑化”，再在整個(gè)求解域上積分，便得到了表征各個(gè)質(zhì)點(diǎn)的動(dòng)力學(xué)特征函數(shù)[14-16]，即

其中，W為核函數(shù)，使用輔助函數(shù)θ定義為

其中，v為空間維數(shù)；h為光滑長(zhǎng)度。

2.2 刀具切削土壤模型建立

模型的前處理及數(shù)據(jù)的后處理均在Ls-Prepost中進(jìn)行，仿真流程如圖4所示。

土壤SPH模型材料采用MAT147，并采用了Drucke_Prager塑性模型，以Mohr_Cloulomb屈服準(zhǔn)則進(jìn)行修正[17]，考慮了粘塑性和切削應(yīng)變失效等方面的因素，使土壤模型更加符合實(shí)際，使仿真結(jié)果更加可靠。土壤模型參數(shù)如下：

土壤密度/kg·m-3： 2.35

體積模量/Pa： 4.6×107

剪切模量/Pa：
1.86×107

土壤水分含量/%： 3.4

粘塑性參數(shù)：
1.1

土壤比重： 2.79

在Ls-Prepost前處理中，在土壤底部和兩側(cè)面添加固定約束，耕深為230mm，接觸靜摩擦因數(shù)為0.2，動(dòng)摩擦因數(shù)為0.18，深松刀具與土壤接觸方式為點(diǎn)面接觸，最后通過(guò)定義曲線(xiàn)的方式來(lái)定義深松刀具初始速度。

圖4 仿真流程圖Fig.4 Simulation process

根據(jù)參數(shù)分析，切削角取值范圍為15°～25°，刀具最低旋轉(zhuǎn)速度為33.9rad/s，最低前進(jìn)速度為0.5m/s，以牽引力為評(píng)價(jià)指標(biāo)，因素水平表如表2所示。

表2 正交試驗(yàn)因素水平表Table 2 Simulation test factor level

2.3 仿真結(jié)果分析

前處理K文件在Ls-dyna中求解完畢后，深松過(guò)程牽引力隨時(shí)間變化，如圖5所示。求得其平均牽引力為1.092 3kN。

按表2進(jìn)行9組仿真，試驗(yàn)安排及結(jié)果如表3所示。

由表3中極差分析可得：對(duì)牽引力影響顯著性排序?yàn)榍邢鹘?、旋轉(zhuǎn)速度、前進(jìn)速度，同時(shí)得出最優(yōu)組合為A1B3C3，即切削角為15°、旋轉(zhuǎn)速度為58rad/s、前進(jìn)速度為0.7m/s時(shí)，深松牽引力最小。

圖5 牽引力變化圖Fig.5 Traction tendency chart表3 試驗(yàn)安排及結(jié)果Table 3 Experiment arrangement and result

試驗(yàn)序號(hào)A切削角/(°)B旋轉(zhuǎn)速度/rad·s-1C前進(jìn)速度/m·s-1牽引力/kN11111.092321221.176131330.855642121.295752231.164162311.095473111.356483231.264593321.1654K13.1243.7443.544K23.5553.6053.637K33.7863.1163.284k11.0411.2481.181k21.1851.2021.212k31.2621.0391.095R0.2210.2090.118最優(yōu)組合A1B3C3

3 深松刀牽引力試驗(yàn)

3.1 試驗(yàn)?zāi)康募霸囼?yàn)地概況

為驗(yàn)證仿真結(jié)果可信度，為后續(xù)優(yōu)化提供依據(jù)，開(kāi)展深松刀林地作業(yè)試驗(yàn)，驗(yàn)證仿真結(jié)果的正確性。

深松試驗(yàn)于2017年夏季在湖南省某林地進(jìn)行，該地土壤常年未進(jìn)行深松作業(yè)，含水量較低，地面硬度比較大。為了使試驗(yàn)場(chǎng)地土壤的含水量等情況基本一致，選擇同一坡面及等高線(xiàn)平面進(jìn)行試驗(yàn)。為了完成牽引力的測(cè)試，使用1臺(tái)履帶式除草機(jī)和深松機(jī)配合作業(yè)。

試驗(yàn)儀器設(shè)備為F0024艾德堡推拉力計(jì)，量程2 000kg，示值誤差±0.5%；筆記本電腦；皮尺；土壤硬度計(jì)；秒表；標(biāo)桿。

3.2 試驗(yàn)設(shè)計(jì)

3.2.1 性能指標(biāo)

開(kāi)展以刀具切削角、刀具旋轉(zhuǎn)速度及刀具前進(jìn)速度三水平三因素的正交試驗(yàn)，以深松過(guò)程牽引力大小作為評(píng)價(jià)指標(biāo)，試驗(yàn)因素水平表如表2所示。

3.2.2 試驗(yàn)步驟

試驗(yàn)過(guò)程中，采用牽引法測(cè)試牽引力大小，在履帶式除草機(jī)和深松機(jī)之間安裝拉力傳感器，試驗(yàn)時(shí)，牽引車(chē)輛從靜止到勻速的過(guò)程盡可能緩慢，以減少?zèng)_擊載荷影響；記錄下在相同作業(yè)長(zhǎng)度內(nèi)深松機(jī)前進(jìn)阻力及安裝深松刀具后深松牽引力，二者相減即為深松刀具的牽引力。試驗(yàn)過(guò)程如圖6所示。

(a) 試驗(yàn)前硬度測(cè)量 (b) 傳感器安裝調(diào)試圖6 試驗(yàn)過(guò)程Fig.6 Experiment process

3.3 試驗(yàn)結(jié)果及分析

試驗(yàn)以牽引力為試驗(yàn)指標(biāo)，傳感器軟件實(shí)時(shí)記錄牽引力數(shù)值，如圖7所示。

求得深松牽引力平均值為2.63kN，以及深松機(jī)前進(jìn)阻力1.48kN，二者相減即為深松刀牽引力1.15kN。正交試驗(yàn)的安排及結(jié)果如表4所示。

3.4 數(shù)據(jù)分析及對(duì)比

1)得到與仿真結(jié)果相同的最優(yōu)組合A1B3C3，對(duì)牽引力影響顯著性排序?yàn)榍邢鹘?、旋轉(zhuǎn)速度及前進(jìn)速度。

2)仿真誤差按數(shù)值表示方法分為絕對(duì)誤差和相對(duì)誤差，即

絕對(duì)誤差=仿真值-實(shí)測(cè)值

試驗(yàn)與仿真結(jié)果對(duì)比如表5所示。

圖7 牽引力變化圖Fig.7 Traction tendency chart表4 試驗(yàn)安排及結(jié)果Table 4 Experimental arrangement and results

試驗(yàn)序號(hào)A切削角/(°)B旋轉(zhuǎn)速度/rad·s-1C前進(jìn)速度/m·s-1牽引力/kN11111.1521221.2731330.8142121.3852231.3162311.0873111.4183231.3493321.29K13.233.943.64K23.773.923.94K34.043.183.46k11.081.311.21k21.261.311.31k31.351.061.15R0.270.250.16最優(yōu)組合A1B3C3

表5 試驗(yàn)與仿真結(jié)果對(duì)比Table 5 Data analysis and contrastTable

續(xù)表5

4 結(jié)論

1)設(shè)計(jì)了一種立式深松刀具，對(duì)深松過(guò)程工作參數(shù)進(jìn)行分析計(jì)算，得出了切削角范圍為15°～25°，最低旋轉(zhuǎn)速度為33.9rad/s，最低前進(jìn)速度為0.5m/s。

2)分析切削土壤特點(diǎn)，基于SPH法對(duì)刀具切削土壤進(jìn)行了仿真，采用3因素3水平正交試驗(yàn)表模擬仿真過(guò)程，得到切削角對(duì)阻力數(shù)值影響最大，且當(dāng)切削角為15°、旋轉(zhuǎn)速度為58rad/s、刀具前進(jìn)速度為0.7m/s時(shí)，深松阻力最小。

3)開(kāi)展了實(shí)地試驗(yàn)，采用牽引法收集牽引力數(shù)值，與仿真數(shù)據(jù)進(jìn)行對(duì)比，誤差在可接受的范圍內(nèi)，為下一步刀具優(yōu)化提供技了術(shù)依據(jù)。