基于光滑粒子流法的振動旋耕減阻效果仿真分析

摘要：為研究振動旋耕的減阻效果，基于光滑粒子流法對多工況下的旋耕切削過程進行仿真分析。將耕深為0.10 m和0.11 m時仿真得到的刀輥受力與通過經驗公式計算得到的力對比，誤差分別為11.87%、12.86%；土壤破碎線與刀尖入土角垂直，隨深耕增大，土壤粒子速度下降，驗證仿真模型的有效性。隨著耕深和轉速增大，土壤動能、內能、刀具動能均增大；當耕深從0.10 m增大到0.14 m時，土壤獲得的動能占比從77.51%增大到78.86%；當轉速從3.6 r/s增大到5.2 r/s時，土壤獲得的動能占比從79.41%下降為72.95%，表明只增大轉速，對改善土壤耕作效果有限。對由刀輥偏心引起的振動旋耕過程進行仿真分析，當偏心距為0.05 m時，轉速為4.0 r/s、4.4 r/s和4.8 r/s時，平均降阻率為58.46%；當耕深為0.10 m、0.12 m、0.14 m時，降阻率分別為50.6%、46.3%、43.1%，降阻率隨耕深增加而下降。隨著偏心距的增大，刀輥受力呈波動式下降。研究結果可以為振動式旋耕設備的研發提供理論依據。

關鍵詞：旋耕刀；振動旋耕；光滑粒子流法；減阻降耗

中圖分類號：S225.31" " " 文獻標識碼：A" " " 文章編號：2095?5553 （2024） 08?0284?07

Simulation analysis of anti?drag effect of rotary vibrating tillage based on"smoothed particle hydrodynamics method

Liu Xiaochan1， Hao Wenchang1， Wang Yuan2， Wang Yunxiao1， Chen Yong1， Zhang Xiuli1

（1. College of Mechanical and Electrical Engineering， Henan Agricultural University， Zhengzhou， 450002， China；2. Beijing Institute of Space Launch Technology， Beijing， 100076， China）

Abstract： In order to study the anti?drag effect of vibratory rotary tillage， this paper conducted a simulation analysis of the rotary tillage cutting process under multiple working conditions by using the Smooth Particle Hydrodynamics （SPH） method. The simulated forces on the cutter roll， obtained for tilling depths of 0.10 m and 0.11 m， were compared to those calculated with empirical formulas， resulting in errors of 11.87% and 12.86%， respectively. The soil fragmentation line was perpendicular to the blade tip. As the tilling depth increased， the soil particle velocities decreased， confirming the effectiveness of the simulation model. As the tilling depth and rotational speed increased， soil kinetic energy， internal energy and tool kinetic energy all increased. When the tilling depth increased from 0.10 m to 0.14 m， the proportion of soil's acquired kinetic energy increased from 77.51% to 78.86%. When the rotational speed increased from 3.6 r/s to 5.2 r/s， the proportion of soil's acquired kinetic energy decreased from 79.41% to 72.95%， indicating that only increasing the rotational speed had limited effects on improving soil tillage effect. A simulation analysis of the vibratory rotary tillage process induced by the eccentricity of the rotary tiller was conducted. When the eccentricity was 0.05 m and the rotational speed was 4.0 r/s， 4.4 r/s， and 4.8 r/s， the average anti?drag rate was 58.46%. When the tilling depth was 0.10 m， 0.12 m， and 0.14 m， the anti?drag rates were 50.6%， 46.3%， and 43.1%， respectively， and the anti?drag rate decreased with increasing tilling depth. With increasing eccentricity， the forces acting on the cutter roll showed fluctuating decreases. The research results can provide a theoretical basis for the development of vibratory rotary tillage equipment.

Keywords： rotary tiller; vibrating rotary tillage; smoothed particle hydrodynamics method; drag reduction and consumption reduction

0 引言

旋耕是刀具和土壤相互作用的過程，旋耕刀是典型的高耗能部件。旋耕減阻降耗理論和技術的發展，是研究熱點之一[1]。研究表明可以通過刀輥結構優化、仿生減阻[2]、振動減阻和農機農藝結合的方法達到旋耕過程減阻降耗的目的。其中，振動不僅能夠有效地減小土壤切削阻力，還有助于減少土壤對刀具的黏附[3]，達到更好的耕作效果。

振動式土壤切削在農業機械領域的應用研究得到了一定的發展。辛麗麗等[4]建立了考慮分段土壤作用力的振動深松機土壤系統的力學模型，分析了深松鏟入土角、振幅和振頻對切削阻力的影響。王東偉等[5]設計了一款超聲振動土壤切削裝置，其試驗研究表明，土壤硬度越大，振動降阻效果越明顯。王家勝等[6]研究了超聲波高頻振動條件下農機觸土部件與土壤的相互作用力學特性，解析了觸土部件碎土與降阻機理。Totten等[7]研究指出，對于在硬黏性土中的深耕鏟，沿水平方向和機組前進的垂直方向添加振動，均能夠有效地降低機具的牽引阻力。Gupta等[8]研究表明當振動激勵信號的頻率接近土壤的固有頻率時，振動切削的效果最好。王文明等[9]對旋轉中耕機的關鍵部件進行了減阻設計，其中旋轉單體的彈簧剛度對碎土效果的影響最大，作業功耗下降了32%。

振動式土壤切削減阻效果研究主要依賴于數值仿真方法，主要包括有限元法、離散元法和光滑粒子流法（SPH法）。其中有限單元法采用黏彈性模型模擬土壤動態切削行為，但是土壤切削過程是大變形問題，采用有限元仿真土壤切削過程時，存在網格畸變，導致仿真誤差增大[10]；作為一種復雜的多相混合物，土壤具有很強的離散特性，離散元法是分析土壤動態行為變化過程的較理想方法，在農業散體物料和土壤的仿真中得到了廣泛的應用[11]，但離散元法仿真的準確性取決于參數標定的準確性[12]，例如武濤[13]、孫景彬[14]等分別對黏性土壤和黃土高原土壤進行了離散元參數標定。光滑粒子流法（SPH法），是一種無網格方法，在處理大變形、大位移及畸變等方面，具有獨到優勢和特點，可以很好模擬土壤粒子動態行為，可用于觸土部件的動力學行為的仿真。宋魯鵬等[15]基于SPH法對圓盤覆土器切削土壤的行為進行了仿真。Hu等[16]基于SPH法分別對黏性土壤和非黏性土壤的刀具切削過程進行了仿真分析。Major等[17]基于SPH仿真模型，得到了在不同牽引速度下旋轉耕作部件受到的水平力。Li等[18]將SPH法應用于分析旋耕刀片結構優化對旋耕能耗減小的作用。蔣建東等[19]應用SPH法分析了外加激勵的振型、頻率及振幅對旋耕過程土壤切削阻力的影響及變化規律。

但是振動的利用，需要新增機械和控制系統，增加設備的能量消耗，導致設備的動力學特性更加復雜，是振動旋耕技術發展的難點。振動旋耕的實現，需要對振動旋耕理論進行更深入的研究。綜上所述，本文將基于光滑粒子流法對旋耕過程進行仿真分析，研究振動旋耕的減阻效果，為振動式旋耕設備的設計及振動特性分析提供理論基礎。

1 旋耕過程仿真模型的建立

1.1 土壤模型的建立

Mohr_Coulomb 屈服面準則計算如式（1）所示。

[F=-Psinφ+KθJ2-ccosφ=0] （1）

式中： P——壓力，Pa；

[φ]——內摩擦角，rad；

[K（θ）]——張量平面角的函數；

[J2]——應力偏張量的第二不變量；

c——內聚力。

LS-DYNA軟件材料庫中的MAT147 MAT_FHWA_SOIL材料本構模型對標準Mohr_Coulomb屈服面準則進行了修正，修正后的屈服面準則公式為

[F=-Psinφ+J2Kθ2+a2sin2φ-ccosφ=0] （2）

式中： a——修正后的屈服面貼合系數。

基于SPH粒子，建立了土壤模型，土壤模型為體積為1 m×0.5 m×0.8 m立方體，共有204 800個土壤粒子。在LS-DYNA的預處理軟件LS-PrePost軟件中，通過關鍵詞定義土壤的體材料參數如表1所示。

1.2 刀輥模型的建立

1WGQ4-70型手扶式微耕機及其刀輥結構如圖1所示，整機主要結構包括單缸柴油發動機、扶手架、變速箱、耕深限制裝置、機架、刀輥。發動機為單缸柴油發動機，其額定功率為4.1 kW，發動機額定轉速為3 600 r/min。柴油發動機輸出轉速通過變速箱將轉速傳遞給刀輥，驅動微耕機進行旋耕作業。旋耕刀輥通過銷軸與輸出軸連接。

以旋耕刀輥刀軸的旋轉中心為原點，旋耕刀的前進方向是X軸方向，垂直方向是Y軸，假定旋耕刀順時針旋轉，并建立坐標系，建立了旋耕刀端點的運動軌跡計算如式（3）所示。

[x=Rcosωt+Vmty=Rsinωt] （3）

式中： R——刀輥回轉半徑；

t——時間；

[Vm]——旋耕刀前進速度；

[ω]——旋耕刀旋轉角速度。

在SOLIDWORKS軟件中建立旋耕刀輥的三維幾何模型如圖2所示，在旋耕作業時，首先刀尖將土壤表面刺破，刀具進入土壤并進行旋耕運動，利用刀體將土壤推擠破碎，并使土壤憑借慣性飛出，進而使土壤進一步破碎。

將旋耕刀輥的三維模型導入Hypermesh軟件，基于六面體網格，對旋耕刀輥進行網格劃分，建立旋耕過程旋耕刀輥的有限元模型。旋耕刀輥仿真模型的參數如表2所示。

將刀輥模型導入LS-DYNA，得到旋耕刀與土壤相互作用的仿真模型如圖3所示，通過Entity Creation模塊中的*Set Data_*Set_Node，定義所有土壤粒子為組1，左、右和底面土壤粒子為組2。組1不進行位移約束，以組 2 為邊界并對 x、y、z 方向進行位移約束，形成SPH粒子的全約束邊界，用以模擬深層堅硬土壤，防止其余土壤粒子飛散。

將刀具的運動分解為直線運動與圓周運動，直線運動速度為0.25 m/s，圓周運動設置為4 r/s。接觸形式設置為點面接觸（SPH方法默認接觸形式），運算時間為1 s，步長為0.002 s。將設置完成的k文件導入ANSYS LS-DYNA971 求解器中求解，求解完成后使用LS-PrePost進行處理。

2 仿真及結果分析

2.1 仿真結果的驗證

旋耕過程仿真結果如圖4所示。土壤粒子在旋耕刀具的作用下的獲得初始速度，開始四散飛出，顏色深淺則體現了土壤所受應力大小，粒子的分布情況則體現了土壤在旋耕作用下的形變情況。土壤破碎呈現出表層土壤與刀具運動切土同步，而深層土壤稍有滯后，土壤擾動呈現出一個傾角，經測量為32°，而刀具入土時與土壤夾角約為60°，兩者基本相互垂直。表明該角度下的土壤破碎為刀具入土后的刀尖擠壓土壤所產生。

為了更直觀地表現出刀具對不同深度土壤的破碎效果，如圖5所示，現提取模型中深度分別為0.07 m、0.21 m和0.35 m的土壤粒子群A、B和C，并提取粒子群速度曲線。粒子群A、B、C速度呈現逐步降低的趨勢，且獲得速度越來越晚。說明旋耕刀的耕作效果隨土壤深度增大而減小，與經驗所得相符。

當轉速為4 r/s時，旋耕刀軸整體垂直方向的受力，如圖6所示，可以看出，刀軸的受力呈現出一定的周期性，這也是導致微耕機旋耕過程振動的原因之一。

刀具與土體之間的相互作用力相當復雜，發動機的轉速和微耕機的前進速度對切削力都有很大的影響，很難直接測試或計算得到精確的土壤切削力。微耕機刀輥的切削力經驗計算如式（4）所示。

[T=9 550PwnQ=T0.9Rθs=arccosR-hR-20°Fs=Qsinθs] （4）

式中： T——微耕機刀輥產生的扭矩，N ? m；

[Pw]——柴油機功率，kW；

n——發動機轉速，r/min；

Q——切割力，N；

R——刀具旋轉半徑，m；

h——工作深度，m；

[Fs]——切割力Q在y方向上的分量，N;

[θs]——切割力與刀輥的中心的合力所形成的垂直平面與直線之間的夾角，rad。

將通過式（4）計算得到的旋耕刀輥受力與仿真得到的刀輥受力進行對比，如表3所示，誤差分別為11.87%和12.86%，驗證了仿真模型的有效性。

2.2 參數對旋耕刀輥受力的影響

為了討論轉速對旋耕刀輥受力的影響，分別對轉速為3.6 r/s、4 r/s、4.4 r/s、4.8 r/s和5.2 r/s時的旋耕過程進行仿真，并計算旋耕刀輥受力的平均值，結果如圖7所受，隨著轉速增加，刀具垂直方向受力也增加，但隨著轉速的增大刀具受力增大的速度減緩，刀具受力隨轉速的增加不是線性的。

對轉速為4 r/s，耕深分別為0.10 m、0.11 m、0.12 m、0.13 m和0.14 m時的旋耕過程進行仿真，刀輥垂直方向所受力如圖8所示，隨著耕深增大，刀具垂直方向受力也同樣增大，但是增大幅度較小，耕深對刀輥阻力的影響有限。

2.3 參數對旋耕過程能耗的影響

在LS-PrePost中提取各個工況下旋耕刀輥切土過程中一個切土周期的總能耗，主要包括刀具的動能、土壤的動能、系統的內能，內能主要是旋耕刀切削土壤過程中的摩擦產熱。不同耕深工況下，旋耕過程的能量分布如表4所示，隨著耕深從0.10 m增大到0.14 m，刀具動能、土壤動能、內能和總能量都增大，同時，土壤動能占總能量的比值從77.51%增大到78.86%，這表明增大耕深可以在一定程度上有效提高耕作效果。

隨著轉速的增大，旋耕過程能量分布如表5所示，可以看出隨著轉速從3.6 r/s增到5.2 r/s，刀具動能、土壤動能與內能均增大，說明提高轉速確實有助于提高土壤耕作效果，但是土壤動能占比從79.41%下降為72.95%，這說明單純提高轉速來改善土壤耕作效果不可行。

3 振動旋耕減阻效果分析

振動式切削是已經被驗證有效減少土壤耕作阻力和能源消耗的方法。對于小型微耕機，通過刀軸偏心的方式，實現刀輥的振動，結構簡單，體積小、質量小、振幅可控。刀輥偏心時，刀尖的運動計算如式（5）所示。

[x=（R+e）cosωt+Vmty=（R+e）sinωt] （5）

式中： e——刀軸的偏心量。

將刀輥偏心情況下的運動數據直接導入LS-DYNA的前處理軟件LS-PROPOST的關鍵詞CURVE中，然后將設置完成的k文件導入ANSYS-LS-DYNA971求解器中進行求解，對振動旋耕過程進行仿真分析。

在刀輥的偏心距為0.05 m時，旋耕刀輥垂直方向受力如圖9所示，與無偏心時對比，刀輥受力整體減小，但是刀輥受力的波動增大。

當存在偏心距時，刀輥受力相對無偏心距時減小的比值，定義為降阻率。分別對偏心距為0.05 m，轉速為4.0 r/s、4.4 r/s和4.8 r/s時的振動旋耕過程進行仿真分析，并提取旋耕刀輥受力的平均值，計算對應的降阻率。如圖10所示，隨著轉速的增大，刀具垂直方向受力先下降后上升，但其整體波動并不大；與無偏心時對比，刀輥垂直方向受力大幅減小，減阻效果明顯，對應的降阻率分別為50.6%、64.9%、59.9%，平均降阻率為58.46%。

當刀輥的偏心距為0.05 m時，分別對耕深為0.10 m、0.12 m和0.14 m時的工況進行仿真分析，提取旋耕刀輥受力的平均值，并計算對應的降阻率。如圖11所示，刀輥垂直方向受力隨耕深的增大而增大，與無偏心相比，對應的降阻率分別為50.6%、46.3%和43.1%，平均降阻率為46.67%，降阻率隨耕深增大而減小，減阻效果明顯。

在水平前進速度為0.25 m/s，旋轉速度4 r/s的情況下，分別對偏心距為0.04 m、0.045 m、0.05 m、0.055 m和0.06 m情況下，刀輥的受力進行仿真分析，刀輥受力如圖12所示。

與無振動情況時相比，刀輥受力均有所下降，說明偏心式振動在旋耕減阻是有效可行的。但是隨著偏心距的增加，旋耕刀具收到的阻力并不是直線下降的，而是呈現出波動式的下降。這種波動式下降的原因，尚待討論，但是這將導致振動旋耕對設備動力學特性的影響更加復雜。

4 結論

為研究振動旋耕的減阻效果，基于光滑粒子流法，對旋耕刀切削土壤過程進行仿真。

1） 耕深為0.10 m和0.11 m時，仿真得到的土壤切削力與根據經驗公式計算得到的切削力相比，誤差分別為11.87%和12.86%；土壤擾動線傾角為32°，垂直于刀尖入土角；隨著耕深增加，土壤粒子速度逐步下降。土壤粒子的運動規律與刀具切削土壤的實際情況相符，驗證了模型的有效性。

2） 隨著轉速和耕深的增加，刀輥受力逐漸增大。當耕深從0.10 m增大到0.14 m時，土壤獲得的動能占比從77.51%增大到78.86%；當轉速從3.6 r/s增大到5.2 r/s時，土壤獲得的動能占比從79.41%下降為72.95%。這表明增大耕深可以提高耕作效果，但是僅僅增加轉速不能提高耕作效果。

3） 通過使旋耕刀軸偏心，實現振動旋耕，建立了振動旋耕過程的仿真模型。與無振動旋耕對比，當偏心距為0.05 m時，轉速為4.0 r/s、4.4 r/s、4.8 r/s時，平均降阻率為58.46%；當耕深為0.10 m、0.12 m、0.14 m時，降阻率分別為50.6%、46.3%和43.1%，降阻率隨耕深增加而下降；隨著偏心距增加，刀輥受力呈波動式下降，使得振動旋耕對設備動力學特性的影響更加復雜。

參 考 文 獻

[ 1 ] 孫志全， 段潔利， 楊洲. 旋耕減阻降耗技術研究進展[J]. 中國農機化學報， 2021， 42（1）： 37-45.

Sun Zhiquan， Duan JieLi， Yang Zhou. Research progress of drag reduction and consumption reduction technology for rotary tillage [J]. Journal of Chinese Agricultural Mechanization， 2021， 42（1）： 37-45.

[ 2 ] Yu H， Han Z， Zhang J， et al. Bionic design of tools in cutting： Reducing adhesion， abrasion or friction [J].Wear， 2021， 482（10）： 203955.

[ 3 ] 劉國陽， 夏俊芳， 鄭侃， 等. 振動橫擋阻隔式旋耕防粘結刀輥設計與試驗[J]. 農業工程學報， 2022， 38（23）： 29-40.

Liu Guoyang， Xia Junfang， Zheng Kan， et al. Design and experiments of the barrier type rotary anti?adhesion blade roller with vibration crosspiece [J]. Transactions of the Chinese Society of Agricultural Engineering， 2022， 38（23）： 29-40.

[ 4 ] 辛麗麗， 李傳琦， 梁繼輝， 等. 考慮分段土壤作用力的振動減阻分析[J]. 農業機械學報， 2014， 45（2）： 136-140.

Xin Lili， Li Chuanqi， Liang Jihui， et al.Vibrating drag reduction considering acting force of piecewise soil [J].Transactions of the Chinese Society for Agricultural Machinery， 2014， 45（2）： 136-140.

[ 5 ] 王東偉， 王家勝.基于超聲振動的土壤切削挖掘裝置設計與試驗[J]. 農業機械學報， 2020， 51（11）： 85-92.

Wang Dongwei， Wang Jiasheng. Design and test of soil cutting and digging device based on ultrasonic vibration [J]. Transactions of the Chinese Society for Agricultural Machinery， 2020， 51（11）： 85-92.

[ 6 ] 王家勝， 王東偉， 趙智豪.超聲波振動下農機觸土部件—土壤相互作用力學特性[J]. 農業工程學報， 2021， 37（23）： 35-41.

Wang Jiasheng， Wang Dongwei， Zhao Zhihao. Mechanical properties of soil?engaging components interacted with ultrasonic vibration in agricultural machinery [J]. Transactions of the Chinese Society of Agricultural Engineering， 2021， 37（23）： 35-41.

[ 7 ] Totten D S， Kaufman L C. Dynamic response of a clay soil to an oscillating subsoiler [J]. Transaction of the ASAE， 1971， 14 （5）： 793-796.

[ 8 ] Gupta C P， Rajuot D S. Effect of amplitude and frequency on soil break?up by an oscillating tool bin experiment [J]. Soil amp; Tillage Research， 1993， 25（4）： 329-338.

[ 9 ] 王文明， 王天元， 郭兵， 等. 基于振動減阻原理的旋轉中耕機關鍵部件設計與試驗[J]. 農業機械學報， 2019， 50（3）： 35-45.

Wang Wenming， Wang Tianyuan， Guo Bing， et al. Design and experiment of key components of rotary cultivator based on vibration reducing resistance principle [J].Transactions of the Chinese Society for Agricultural Machinery， 2019， 50（3）： 35-45.

[10] Liu Xiaohong， Yu Yan， Qiu Lichun. An overview of the finite element method on the tool?soil interacting problem of tillage [J]. Applied Mechanics and Materials， 2015， 707（2）： 397-400.

[11] 曾智偉， 馬旭， 曹秀龍， 等. 離散元法在農業工程研究中的應用現狀和展望[J].農業機械學報， 2021， 52（4）： 1-20.

Zeng Zhiwei， Ma Xu， Cao Xiulong， et al. Critical review of applications of discrete element method in agricultural engineering [J]. Transactions of the Chinese Society for Agricultural Machinery， 2021， 52（4）： 1-20.

[12] 王憲良， 胡紅， 王慶杰， 等. 基于離散元的土壤模型參數標定方法[J]. 農業機械學報. 2017， 48（12）： 78-85.

Wang Xianliang， Hu Hong， Wang Qingjie，et al. Calibration method of soil contact characteristic parameters based on DEM theory [J].Transactions of the Chinese Society for Agricultural Machinery，2017， 48（12）： 78-85.

[13] 武濤， 黃偉鳳， 陳學深， 等.考慮顆粒間黏結力的黏性土壤離散元模型參數標定[J].華南農業大學學報， 2017， 38（3）： 93-98.

Wu Tao， Huang Weifeng， Chen Xueshen， et al. Calibration of discrete element model parameters for cohesive soil considering the cohesion between particles [J]．Journal of South China Agricultural University， 2017， 38（3）： 93-98.

[14] 孫景彬， 劉琪， 楊福增， 等.黃土高原坡地土壤與旋耕部件互作的離散元仿真參數標定[J].農業機械學報， 2022， 53（1）： 63-73.

Sun Jingbin， Liu Qi， Yang Fuzeng，et al. Calibration of discrete element simulation parameters of sloping soil on loess plateau and its interaction with rotary tillage components [J]．Transactions of the Chinese Society for Agricultural Machinery， 2022， 53（1）： 63-73.

[15] 宋魯鵬， 楊先海， 李倩.基于SPH無網格法的精密播種圓盤覆土器仿真研究[J].中國農機化學報， 2016， 37（2）： 29-32.

Song Lupeng， Yang Xianhai， Li Qian. Simulations study of precision seeding disc covering device based on SPH's mesh less method [J]. Journal of Chinese Agricultural Mechanization， 2016， 37（2）： 29-32.

[16] Hu M， Gao T， Dong X， et al. Simulation of soil?tool interaction using smoothed particle hydrodynamics （SPH） [J]. Soil and Tillage Research， 2023， 229（5）： 105671.

[17] Major T， Csanády V. Combined FEM-SPH simulation method for the modeling of the interaction of tillage tools and the soil [J]. Erdészettudományi K?zlemények， 2015， 5（1）： 7-19.

[18] Li Shoutai， Chen Xiaobing， Chen Wei， et al. Soil?cutting simulation and parameter optimization of handheld tiller's rotary blade by smoothed particle hydrodynamics modelling and Taguchi method [J]. Journal of Cleaner Production， 2018， 179（4）： 55-62.

[19] 蔣建東， 高潔， 趙穎娣， 等. 土壤旋切振動減阻的有限元分析[J]. 農業機械學報， 2012， 43（1）： 58-62.

Jiang Jiandong， Gao Jie， Zhao Yingdi， et al. Finite element simulation and analysis on soil rotary tillage with external vibration excitation [J].Transactions of the Chinese Society for Agricultural Machinery， 2012， 43（1）： 58-62.