煙稈拔稈破碎機振動特性與刀輥動力學(xué)分析

摘要：拔稈刀輥作為一體式煙稈拔稈破碎機的關(guān)鍵部件，直接決定煙稈的抓取性能。由于作業(yè)環(huán)境、土壤特性及煙稈特征參數(shù)的差異性，致使樣機容易發(fā)生共振且難以對拔稈刀輥動力性能和功耗進(jìn)行準(zhǔn)確計算。為解決上述問題，提出利用EDEM數(shù)值模擬方法計算拔稈刀輥的作業(yè)功耗。首先構(gòu)建拔稈刀輥旋耕刀、拔稈橫刀工作阻力的數(shù)學(xué)模型來分析刀輥的運動軌跡和功耗，再利用EDEM仿真計算刀輥切削煙壟的阻力、扭矩和功耗，最后開展整機振動測試研究，驗證樣機的振動特性及可靠性。結(jié)果表明，EDEM仿真的刀輥切削煙壟平均功耗為3.02 kW，與理論計算值比較其誤差為5.63%；整機的振動幅值在83.894 Hz或其倍頻處達(dá)到峰值，避開激勵頻率和一階固有頻率區(qū)間，樣機不會發(fā)生共振，具有良好的可靠性。

關(guān)鍵詞：煙稈拔稈破碎機；拔稈刀輥；動力學(xué)；振動特性

中圖分類號：S225.3

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

文章編號：2095-5553 （2024） 12-0009-07收稿日期：2024年6月6日

修回日期：2024年8月6日

*基金項目：國家重點研發(fā)計劃項目（2021YFD1300300）；貴州省高層次創(chuàng)新型人才項目（黔科合平臺人才—GCC［2022］007—1）；貴州開放大學(xué)（貴州職業(yè)技術(shù)學(xué)院）科學(xué)研究項目（2023ZD01）；茅臺學(xué)院高層次人才科研啟動經(jīng)費項目（mygccrc ［2022］110）

第一作者：楊桃月，女，1982年生，貴州天柱人，高級實驗師；研究方向為機械裝備設(shè)計。E-mail：44719998@qq.com

通訊作者：曹陽，男，1981年生，湖南益陽人，博士，副教授，碩導(dǎo)；研究方向為農(nóng)業(yè)智能裝備。E-mail：yangc377@163.com

Vibration characteristics and dynamic analysis of knife roll of tobacco stalk pulling and crushing machine

Yang Taoyue¹， Zhang Dabin²， Shu Chengsong²， Yu Chaojing³， Cao Yang²

（1. School of Aviation， Guizhou Open University （Guizhou Vocational Technology Institute）， Guiyang， 550023， China；2. School of Mechanical Engineering， Guizhou University， Guiyang， 550025， China; 3. Maotai College， Zunyi， 564507， China）

Abstract： As the key component of the integrated tobacco stalk pulling crusher， the stalk pulling knife roller directly determines the performance of tobacco straw grabbing. As the operating environment， soil characteristics and the differences of characteristic parameters of tobacco stalks， the prototype is easily resonated， and it is difficult to accurately calculate the dynamic performance and power consumption of straw pulling knife roller. In order to solve these problems， this paper proposed to calculate the working power consumption of straw pulling knife roller by using EDEM numerical simulation method. Firstly， the mathematical models of rotary knife and straw pulling horizontal knife working resistance were constructed to analyze the movement trajectory and power consumption of the knife rollers， and then the resistance， torque and power consumption of the knife rollers cutting tobacco ridge were calculated by using EDEM simulation. Finally， the whole machine vibration test research was carried out to verify the vibration characteristics and reliability of the prototype. The results show that the average power consumption of EDEM simulation of knife roller cutting tobacco ridge is 3.02 kW， and its error is 5.63% compared with the theoretical calculated value. The vibration amplitude of the whole machine peaks at 83.894 Hz or its multiple frequency， which avoids the excitation frequency and the first-order intrinsic frequency interval， and the prototype will not resonate， and has good reliability.

Keywords： tobacco stalk pulling and crushing machine; stalk pulling knife roller; dynamics; vibration characteristics

0 引言

拔稈刀輥作為一體式煙稈拔稈破碎機拔稈的關(guān)鍵部件，直接決定煙稈拔稈的漏稈率及運輸破碎效果。由于刀輥的作業(yè)過程為一個較為復(fù)雜的運動系統(tǒng)，受土壤特性、刀輥結(jié)構(gòu)及功耗等多種因素的耦合作用影響，通過常規(guī)理論公式難以準(zhǔn)確計算拔稈刀輥的作業(yè)功耗^［1］。通過有限元的方法進(jìn)行功耗模擬是一種有效的手段，但目前的仿真方法準(zhǔn)確性較差。此外，煙稈拔稈破碎機振動對樣機的振動疲勞及作業(yè)性能有重要影響，設(shè)計時應(yīng)保證各運動部件的激勵頻率避開樣機的固有頻率，防止發(fā)生共振^［2］。針對上述問題，通過構(gòu)建準(zhǔn)確的刀輥運動學(xué)模型和仿真模型來進(jìn)行刀輥功耗仿真計算，開展煙稈拔稈破碎機整機振動特性分析，對提高樣機作業(yè)穩(wěn)定性和降低功耗具有重要意義。

目前國內(nèi)外學(xué)者利用離散元法（distinct elementmethod， DEM）對土壤—機具、秸稈—土壤—機具的相互作用進(jìn)行了相關(guān)研究^［3-5］。如方長敏等^［3］利用EDEM 對秸稈—土壤—旋耕刀系統(tǒng)中的秸稈位移進(jìn)行仿真分析，得到初始位置不同的秸稈顆粒具有不同的運動。廖慶喜^［6］、王金武^［7］等利用EDEM研究了排種器排種性能的數(shù)值模擬與試驗。由于煙稈獨特的力學(xué)特性和煙蔸根系發(fā)達(dá)性，目前尚未發(fā)現(xiàn)針對拔稈刀輥切削煙壟功耗分析的研究。針對農(nóng)業(yè)機械振動特性國內(nèi)外研究者做了大量的研究，主要集中在人機舒適性、動力執(zhí)行元件振動特性^{［8， 9］}等方面，如龐靖等^［10］基于變分模態(tài)分解方法對谷物清選篩進(jìn)行振動分析與結(jié)構(gòu)優(yōu)化以減小系統(tǒng)的不平衡振動。在現(xiàn)有公開報道中尚未見到針對一體式煙稈拔稈破碎機整機振動特性方面的系統(tǒng)研究。因此，在分析拔稈刀輥工作機理的基礎(chǔ)上，研究作業(yè)過程中刀輥切削煙壟的功耗及運動規(guī)律對拔稈機構(gòu)作業(yè)性能的影響，并分析整機的振動特性具有較好的研究價值。

為進(jìn)行拔稈刀輥功耗數(shù)值模擬與整機的振動特性分析，本文以自行研制的二代一體式煙稈拔稈破碎機為研究載體，建立旋耕刀和拔稈橫刀工作阻力的數(shù)學(xué)模型，分析刀輥橫刀運動軌跡和刀輥功耗。在此基礎(chǔ)上運用EDEM離散元軟件建立旋耕刀切土仿真模型進(jìn)行拔稈刀輥作業(yè)功耗的數(shù)值模擬，并開展一體式煙稈拔稈破碎機整機振動特性試驗，分析其振動特性，驗證整機的可靠性。

1 一體式煙稈拔稈破碎機結(jié)構(gòu)與工作原理

一體式煙稈拔稈破碎機由拖拉機通過三點懸掛機構(gòu)牽引作業(yè)，其拔稈深度可通過拖拉機懸掛高度進(jìn)行調(diào)節(jié)。在拔稈作業(yè)時拖拉機沿著煙壟方向前進(jìn)將煙稈壓倒，并利用萬向節(jié)與分動箱將拖拉機的輸出動力傳送至拔稈刀輥帶動其轉(zhuǎn)動。拔稈作業(yè)時在刀輥左右旋耕刀的作用下，實現(xiàn)煙壟土壤切削及疏松并將煙蔸掘出地面，并通過拔稈橫刀與旋耕刀的相互作用將整顆煙稈拋送至鏈?zhǔn)秸駝虞斔蜋C構(gòu)實現(xiàn)煙稈傳輸及泥土、碎石等去除。最后在上下兩對夾持喂入對輥的作用下將煙稈投入至破碎機進(jìn)行煙稈粉碎，并將煙稈顆粒收集至收集箱，實現(xiàn)煙稈的連續(xù)拔稈破碎作業(yè)。設(shè)計的一體式煙稈拔稈破碎機對煙壟及煙稈具有較好的適應(yīng)性，可實現(xiàn)煙壟耕整及煙稈的拔稈破碎功能，整體上具有拔稈效率高、整機功耗低、適合山地農(nóng)業(yè)環(huán)境作業(yè)等優(yōu)勢，其原理如圖1所示。

2 拔稈刀輥動力學(xué)分析

2.1 旋耕刀阻力分析

拔稈刀輥主要通過旋耕刀及拔稈橫刀切削煙壟和掘起煙蔸，并將其輸送至傳輸機構(gòu)。對整個拔稈刀輥進(jìn)行力學(xué)性能分析，對于單片的旋耕刀，其作業(yè)時主要受到切削土壤的阻力F_c、擠壓土壤的阻力F_r、旋耕刀與土壤的摩擦力F_f和旋耕刀刨土的阻力F_i的相互作用，相關(guān)阻力的大小與旋耕刀的設(shè)計參數(shù)、土壤力學(xué)特性及旋耕切削參數(shù)密切相關(guān)^［11］，建立的旋耕刀作業(yè)力學(xué)模型如圖2所示。

通過圖2的力學(xué)性能分析可知拔稈刀輥旋耕刀所受工作阻力F可表示為

{F_x=-F_ccosη+F_rcos（2π-θ+γ₀）-

F_fsin（2π-θ+γ₀）+

F_icosπ/2-η+γ/2

F_y=-F_csinη-F_rsin（2π-θ+γ₀）+

F_fcos（2π-θ+γ₀）-

F_isinπ/2-η+γ/2

F=F_x²+F_y² （1）

式中：F_x——旋耕刀在x方向受到的阻力，N；

F_y——旋耕刀在y方向受到的阻力，N；

F——旋耕刀受到的阻力合力，N；

θ——刀輥上正反兩把相鄰旋耕刀之間的夾角，（°）；

γ₀——旋耕刀安裝角，（°）；

γ——旋耕刀正切面的切土角，（°）；

η——旋耕刀背切面的切土角，（°）。

2.2 橫刀阻力分析

拔稈橫刀作用功能為將旋耕刀切削的煙壟進(jìn)行二次切削并實現(xiàn)拋土，其阻力分析如圖3所示。在拖拉機牽引力作用下沿著煙壟方向切削土壤時受到的水平阻力F₁和豎直方向的阻力F₂、橫刀前刀面法向載荷F_N和橫刀與土壤之間的摩擦阻力f^［12］。

由圖3的受力分析可知，煙壟方向切削土壤時受到的水平阻力F₁和豎直方向的阻力F₂可表示為

F₁=F_p-F_N²+f²sinα+arccosF_N/F_N²+f²） （2）

F₂=G－F_N²+f²sinα－arcsinF_N/F_N²+f²） （3）

式中：F_p——拔稈刀輥所受的拖拉機牽引力，N；

G——拔稈刀輥的重力，N；

α——拔稈橫刀入土切削角，（°）。

而橫刀與土壤之間的摩擦阻力f與橫刀前刀面的法向載荷F_N之間的關(guān)系可表示為f=μF_N，其中μ為煙壟土壤與橫刀之間的摩擦系數(shù)，因此可得到橫刀在水平和豎直方向所受的阻力最大值為

F₁=F_p+F_N1+μ² （4）

F₂=G+F_N1+μ² （5）

將式（4）和式（5）進(jìn)行矢量合成，可計算出拔稈橫刀作業(yè)時受到的阻力為

F_p²+G²－2F_N（F_p+G）1+μ²+2F_N²（1+μ²）lt;F_合lt;F_p²+G²+2F_N（F_p+G）1+μ²+2F_N²（1+μ²） （6）

2.3 拔稈橫刀運動學(xué)分析

一體式煙稈拔稈破碎機作業(yè)拔稈多作業(yè)時，刀輥橫刀的絕對運動主要為拖拉機牽引的沿?zé)焿欧较虻乃竭\動和刀輥自身旋轉(zhuǎn)運動的合成。要完成煙稈的有效掘起并傳輸至傳輸機構(gòu)，減少漏稈率，要求橫刀端面的線速度大于機組行進(jìn)速度。以機組前進(jìn)方向為x方向建立拔稈橫刀的運動學(xué)模型，如圖4所示，可看出拔稈橫刀的運動軌跡為余擺線，在任意t時刻，拔稈橫刀刀面端點的位移平衡方程為

{x=V_mt+R′cosωt

y=R′sinωt （7）

式中：V_m——拔稈破碎機前進(jìn)速度，m/s；

R′ ——拔稈橫刀旋轉(zhuǎn)半徑，mm；

ω——刀輥角速度，rad/s。

由式（7）可獲得橫刀端面的速度為

{v_x=V_m-R′ωsinωt

v_y=R′ωcosωt （8）

為實現(xiàn)拔稈刀輥的連續(xù)穩(wěn)定作業(yè)，將煙稈拔起并拋送至后續(xù)的鏈?zhǔn)秸駝觽鬏敊C構(gòu)，按上述分析可知，必須保證 R′ωsinωtgt;V_m。

3 拔稈刀輥功耗分析

一體式煙稈拔稈破碎機的拔稈刀輥作業(yè)原理為通過刀輥反轉(zhuǎn)從煙壟底部開始切削土壤，并將煙蔸從煙壟中掘出，在拔稈橫刀與齒梳的共同作用下將煙稈拋送至鏈?zhǔn)秸駝虞斔蜋C構(gòu)完成拔稈作業(yè)。通過旋耕刀、拔稈橫的阻力分析及運動學(xué)分析可知，拔稈刀輥的功耗主要為旋耕刀及橫刀切削煙壟克服土壤阻力及拋送煙稈所做的功，即切土功耗和拋土功耗，可表達(dá)為

{P=1/η_dBd/100V_mp₀+v_d²/2δ

v_d=πRn/30 （9）

式中：P——刀輥功耗，kW；

η_d——傳動效率，取值0.8～0.9；

v_d——刀輥線速度，m/s；

R——刀輥半徑，m；

n——刀輥轉(zhuǎn)速，r/min；

p₀——旋耕切土比阻，N/m²；

B——煙壟耕幅，m；

d——煙壟耕深，cm。

根據(jù)作業(yè)參數(shù)及設(shè)計參數(shù)，煙稈拔稈破碎機的前進(jìn)速度、半徑和轉(zhuǎn)速取值為0.37 m/s，0.25 m，110 r/min，結(jié)合貴州煙草種植條件及土壤環(huán)境，經(jīng)實測得出旋耕切土比阻、耕幅和耕深值分別為4.1×10⁴ N/m²、0.7 m、20 cm，由式（9）可計算出拔稈刀輥的功率消耗為2.85 kW。

4 拔稈刀輥切削煙壟數(shù)值模擬

4.1 仿真模型的建立及參數(shù)設(shè)置

在有限元仿真軟件中建立拔稈刀輥切削煙壟的仿真模型，刀輥材料為45#鋼。模型構(gòu)建時土壤顆粒平均直徑D為10 mm，土壤顆粒按四面體和立方體排列，其顆粒數(shù)量N由顆粒直徑和構(gòu)成的立方體或正四面體邊長決定，假設(shè)立方體（數(shù)量為N₁）和正四面體（數(shù)量為N₂）的土壤顆粒之間呈均等排列，則實際土壤球體顆粒數(shù)可表示為^［13］

N=（N₁+N₂）/2≈1.2V/D³ （10）

式中：V——土壤顆粒在土槽中所占的體積，mm³。

通過計算得出切削的煙壟土壤顆粒為290 000顆，將仿真中的拔稈刀輥前進(jìn)速度、旋轉(zhuǎn)速度和切削時間分別設(shè)置為0.37 m/s、110 r/min和1～3.5 s。土壤顆粒與刀輥和土槽之間的接觸類型設(shè)置為Hert-Mindin wit Bonding接觸，仿真步長為3.738 42×10^-5，構(gòu)建的拔稈刀輥—煙壟的離散元模型如圖5所示。

4.2 數(shù)值模擬結(jié)果分析

選取0～3 s內(nèi)任一時刻拔稈刀輥反轉(zhuǎn)切削煙壟的模擬過程如圖6所示。由圖6可知，刀輥沿著切削煙壟方向能夠?qū)崿F(xiàn)土壤的掘起和拋送。

拔稈刀輥的功耗主要為克服旋耕刀和拔稈橫刀旋轉(zhuǎn)的動能及土壤掘起所做的功，主要由拔稈刀輥刀軸的扭矩T和水平阻力F_z決定，其中F_z為旋耕刀和拔稈橫刀在水平方向受到的阻力F_x和F₁的合力。因此，利用EDEM后處理模塊提取拔稈刀輥切削煙壟時的阻力和曲線，如圖7所示。

由圖7可知，拔稈刀輥在切削煙壟及煙稈拔稈過程中的水平阻力呈周期性交替變化，這主要是由于刀輥在旋轉(zhuǎn)過程中切削煙壟和拋土所需的功耗不同所致。在切削過程中，阻力逐漸增大，在完全切削入煙壟后阻力值穩(wěn)定在1 189 N附近，且交替變化的阻力值分布趨于穩(wěn)定，這主要是由于沿著機組的前進(jìn)方向刀輥上均勻布置的旋根刀依次切削土壤，并在拔稈橫刀切削煙壟至最大耕深時達(dá)到峰值。而針對刀輥的扭矩，從旋耕刀及拔稈橫刀開始切削煙壟開始，逐漸由0 N·m增大，直至完全切入煙壟后趨于穩(wěn)定，呈現(xiàn)與阻力相同的周期性變化規(guī)律，平均扭矩值為43 N·m，檔拔稈橫刀切入至煙壟底部時扭矩達(dá)到峰值。

利用獲取的扭矩T和水平分力F_z可計算出拔稈刀輥切削煙壟的功率曲線如圖8所示。由圖8可知，刀輥的功耗在切削煙壟的過程中呈周期性變化，與扭矩變化趨勢保持一致，在2～3.5 s時趨于穩(wěn)定，此時刀輥已完全切削入煙壟，平均功耗數(shù)值為3.02 kW，與理論計算的功耗2.85 kW相比，誤差為5.63%。誤差的主要來源主要是仿真模型對煙壟土壤及條件做了一定的簡化，且實際作業(yè)環(huán)境中土壤的粘度、土質(zhì)環(huán)境與仿真有差異所致，總體而言，利用離散元數(shù)值模擬方法來研究拔稈刀輥切削煙壟的功耗具有一定的可行性，可為煙稈拔稈破碎機的功耗設(shè)計提供參考依據(jù)。

5 整機振動特性測試與分析

5.1 測點布置及測試方案

田間試驗在貴州省平壩縣煙葉站煙田開展，拔稈破碎作業(yè)時拖拉機主軸輸出軸最大轉(zhuǎn)速720 r/min，帶動破碎機的汽油機轉(zhuǎn)速為1 600 r/min，利用DH5925動態(tài)信號分析儀和DH3111E三相加速度傳感器（X、Y、Z三個方向，對應(yīng)拔稈破碎機的橫向、前進(jìn)方向、豎直方向）對一體式煙稈拔稈破碎機振動信號進(jìn)行采集，采集工況為田間作業(yè)實際工況，在樣機上布置9個測點，測點布置位置如圖9所示。

各測點的振動強度用加速度均方根值來表示^［14］，可表達(dá)為

a=a_x²+a_y²+a_z²/3 （11）

式中：a——測點振動強度總量，m/s²；

a_x——X方向的加速度均方根值，m/s²；

a_y——Y方向的加速度均方根值，m/s²；

a_z——Z方向的加速度均方根值，m/s²。

5.2 結(jié)果及分析

圖10為煙稈拔稈破碎機上布置的9個測點的振動總量加速度均方根值，可以看出，測點9的振動強度最激烈，達(dá)到13.8 m/s²，測點1、測點3、測點4、測點5的振動強度相近。其中測點1為拔稈刀輥旋轉(zhuǎn)部件及三點懸掛機構(gòu)連接處，振動強度僅次于測點9，且在X、Y、Z三個方向的振動加速度均方根值相接近，如圖11（a）所示。

測點1的振動主要是拖拉機發(fā)動機振動、刀輥切削煙壟和煙稈拔稈作業(yè)、壟溝不平衡產(chǎn)生的激勵頻率的振動的耦合作用，即拔稈作業(yè)時拖拉機發(fā)動機的振動通過三點懸掛機構(gòu)傳遞至測點1，刀輥低速旋轉(zhuǎn)切削煙壟和煙稈掘蔸刨土產(chǎn)生的低頻振動與樣機的低階固有頻率相接近，故而加劇了振動強度。測點3、測點4、測點5為破碎機附近位置測點，其振動強度相對較大的主要原因為破碎機刀盤高速旋轉(zhuǎn)切割煙稈時產(chǎn)生的不平衡慣性力產(chǎn)生的激勵振動。測點9布置在汽油機支撐架上，汽油機的高速旋轉(zhuǎn)會產(chǎn)生激烈振動，因此需在汽油機支撐架與機架橫梁的底座安裝處增加減震器、阻尼塊等減振裝置來實現(xiàn)汽油機振動的隔振。

將采集的測點1～測點9的振動加速度均方根值進(jìn)行傅里葉變換，可獲得各測點的頻域信號，見圖11（b）及圖12。由圖12可以發(fā)現(xiàn)，測點9的振動幅值最大，說明破碎機汽油機的振動是整機的主要振源，由圖11（b）可知汽油機的振動主要表現(xiàn)為豎直方向振動，其振幅為9.66 m/s²，是X和Y方向的3倍左右。由此可見，汽油機豎直方向振動是整機的主要振源，為減少劇烈振動而引起整機振動疲勞，在汽油機與其支架之間增加磁力隔振墊，并在汽油機支架與整機機架接觸處增加阻尼塊及相應(yīng)的隔振結(jié)構(gòu)，可減少汽油機振動源的振動傳遞。此外，測點3、測點5、測點6在X或Y方向產(chǎn)生的振動幅值也高于其它測點，主要是由于輸送機構(gòu)、破碎機構(gòu)等動力執(zhí)行部件相對于機架存在一定的偏心，后續(xù)可將動力執(zhí)行部件居中安裝布置消除偏心而減少振動。通過頻域分析可知，一體式煙稈拔稈破碎機上各測點均在振動頻率為83.894 Hz及其倍頻（35.156 Hz、128.909 Hz、339.851 Hz）處振動的振幅峰值在某一方向上達(dá)到最大。針對整個樣機而言，各回轉(zhuǎn)部件產(chǎn)生的振動頻率與其轉(zhuǎn)速呈正比，轉(zhuǎn)速最大的部件為破碎機粉碎煙稈，其轉(zhuǎn)速高達(dá)1 600 r/min，此時產(chǎn)生的振動頻率為 26.67 Hz，在ANSYS中對樣機進(jìn)行模態(tài)分析，獲得機架的一階固有頻率為39.653 Hz，激振頻率遠(yuǎn)低于樣機一階固有頻率，因此，樣機不會發(fā)生共振，驗證了設(shè)計的一體式煙稈拔稈破碎機具有良好的可靠性。

6 結(jié)論

1） 研制一體式煙稈拔稈破碎機，建立旋耕刀、拔稈橫刀工作阻力及刀輥運動軌跡的運動學(xué)模型，通過理論計算得出拔稈刀輥的功率消耗為2.85 kW。

2） 在離散元軟件 EDEM 中對刀輥切削煙壟的阻力、扭矩和功耗進(jìn)行數(shù)值計算，獲得刀輥穩(wěn)定切削煙壟拔稈作業(yè)的平均功耗為3.02 kW，與理論計算的誤差為5.63%。

3） 對一體式煙稈拔稈破碎機進(jìn)行田間試驗振動測試，得出整機上各測點均在振動頻率為83.894 Hz及其倍頻處振動的振幅峰值達(dá)到最大，設(shè)計的樣機避開激勵頻率和一階固有頻率區(qū)間，樣機不會發(fā)生共振，具有良好的可靠性。

參 考 文 獻(xiàn)

［1］ 方會敏， 姬長英， 張慶怡， 等. 基于離散元法的旋耕刀受力分析［J］. 農(nóng)業(yè)工程學(xué)報， 2016， 32（21）： 54-59.

Fang Huimin， Ji Changying， Zhang Qingyi， et al. Force analysis of rotary blade based on distinct element method ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2016， 32（21）： 54-59.

［2］ 楊桃月， 顧苑婷， 舒成松， 等. 煙稈拔稈破碎機機架輕量化設(shè)計與隨機振動疲勞分析［J］. 中國農(nóng)機化學(xué)報， 2021， 42（5）： 100-106.

Yang Taoyue， Gu Yuanting， Shu Chengsong， et al. Lightweight design and random vibration fatigue analysis for the frame of tobacco stalk pulling and crushing machine ［J］. Journal of Chinese Agricultural Mechanization， 2021， 42（5）： 100-106.

［3］ 方會敏， 姬長英， Ahmed Ali Tagar， 等. 秸稈—土壤—旋耕刀系統(tǒng)中秸稈位移仿真分析［J］. 農(nóng)業(yè)機械學(xué)報， 2016， 47（1）： 60-67.

Fang Huimin， Ji Changying， Ahmed Ali Tagar， et al. Simulation analysis of straw movementin straw-soil-rotary blade system ［J］. Transactions of the Chinese Society for Agricultural Machinery， 2016， 47（1）： 60-67.

［4］ Li J， Webb C， Pandiella S S， Campbell G M. Discrete particle motion on sieves-a numerical study using the DEM simulation ［J］. Powder Technology， 2003， 133： 190-202.

［5］ Cleary P W， Sawley M L. DEM modeling of industrial granular flows： 3D case studies and the effect of particle shape on hopper discharge ［J］. Applied Mathematical Modeling， 2002， 26（6）： 89-111.

［6］ 廖慶喜， 張朋玲， 廖宜濤， 等. 基于EDEM的離心式排種器排種性能數(shù)值模擬［J］. 農(nóng)業(yè)機械學(xué)報， 2014， 45（2）： 109-114.

Liao Qingxi， Zhang Pengling， Liao Yitao， et al. Numerical simulation on seeding performance of centrifugal rape-seed metering device based on EDEM ［J］. Transactions of the Chinese Society for Agricultural Machinery， 2014， 45（2）： 109-114.

［7］ 王金武， 唐漢， 王奇， 等. 基于EDEM軟件的指夾式精量排種器排種性能數(shù)值模擬與試驗［J］. 農(nóng)業(yè)工程學(xué)報， 2015， 31（21）： 43-50.

Wang Jinwu， Tang Han， Wang Qi， et al. Numerical simulation and experiment on seeding performance of pickup finger precision seed-metering device based on EDEM ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2015， 31（21）： 43-50.

［8］ 張紅梅， 張晨明， 朱晨輝， 等. 玉米免耕播種機振動特性及對排種器排種性能的影響［J］. 中國農(nóng)機化學(xué)報， 2024， 45（5）： 1-8， 337.

Zhang Hongmei， Zhang Chenming， Zhu Chenhui， et al. Analysis of vibration characteristics of corn no-till seeder and its influence on seed metering performance ［J］. Journal of Chinese Agricultural Mechanization， 2024， 45（5）： 1-8， 337.

［9］ 王順， 馬心坦， 賀明佳. 拖拉機駕駛室內(nèi)低頻噪聲控制研究［J］. 中國農(nóng)機化學(xué)報， 2024， 45（5）： 134-139.

Wang Shun， Ma Xintan， He Mingjia. Research on control of low frequency noise in tractor cab ［J］. Journal of Chinese Agricultural Mechanization， 2024， 45（5）： 134-139.

［10］ 龐靖， 林毅， 王升升， 等. 基于VMD的谷物清選篩振動分析與結(jié)構(gòu)優(yōu)化［J］. 農(nóng)業(yè)工程學(xué)報， 2023， 39（12）： 1-9.

Pang Jing， Lin Yi， Wang Shengsheng， et al. Vibration analysis and structure optimization of grain cleaning screen based on VMD ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2023， 39（12）： 1-9.

［11］ 宋建農(nóng)， 李自華. 反轉(zhuǎn)旋耕理論分析［J］. 北京農(nóng)業(yè)工程大學(xué)學(xué)報， 1990， 10（3）： 15-22.

［12］ 余朝靜. 移動式煙稈拔稈破碎機拔稈機構(gòu)的研制及試驗［D］. 貴陽： 貴州大學(xué)， 2018.

Yu Chaojing. Development and experiment of straw pulling mechanism for moving tobacco straw pulling and crushing machine ［D］. Guiyang： Guizhou University， 2018.

［13］ 周百冬， 牛琪， 喻煒， 等. 雙螺旋振動式馬鈴薯供種裝置離散元仿真及優(yōu)化［J］. 中國農(nóng)機化學(xué)報， 2024， 45（4）： 23-31.

Zhou Baidong， Niu Qi， Yu Wei， et al. Discrete element simulation research and optimization of double channel vibration potato seed feeder ［J］. Journal of Chinese Agricultural Mechanization， 2024， 45（4）： 23-31.

［14］ 舒成松. 煙稈拔稈破碎機關(guān)鍵部件結(jié)構(gòu)優(yōu)化及試驗［D］. 貴陽： 貴州大學(xué)， 2019.

Shu Chengsong. Structural optimization and experiment of key components of tobacco straw pulling crusher machine ［D］. Guiyang： Guizhou University， 2019.