馬小英，杜伊健，閆玉濤，邵鴻媚

(1.沈陽理工大學(xué) a.機(jī)械工程學(xué)院；b.環(huán)境與化學(xué)工程學(xué)院，沈陽 110159；2.東北大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院，沈陽 110819)

土壤在受到外力作用時(shí)產(chǎn)生的動(dòng)態(tài)行為變化通常與作用力系統(tǒng)、工況環(huán)境條件及土壤性質(zhì)相關(guān)，曾德超[1]以吉爾與范登堡的分類為基礎(chǔ)，給出了土壤行為的定義，對(duì)土壤行為的動(dòng)態(tài)變化過程、機(jī)理及影響因素進(jìn)行了分析；但這種基于連續(xù)介質(zhì)理論的土壤力學(xué)模型難以應(yīng)對(duì)復(fù)雜的土壤變化。

隨著計(jì)算機(jī)技術(shù)和數(shù)值模擬技術(shù)的的發(fā)展，一種基于不連續(xù)數(shù)值模型的離散元理論被提出，離散單元法(Discrete Element Method，DEM)[2]常用來描述復(fù)雜的散體系統(tǒng)變化，為解決動(dòng)態(tài)響應(yīng)問題提供了新的思路。Ucgul M等[3]采用離散單元法對(duì)土壤-板型犁的相互作用進(jìn)行了模擬，驗(yàn)證了離散單元法對(duì)預(yù)測相互作用模型的有效性。Wang X等[4]利用離散單元法開發(fā)模型，研究了不同模型顆粒半徑對(duì)土壤的破裂距離、擾動(dòng)面積、密度變化率等參數(shù)的影響。王鵬等[5]利用離散單元法研究了鏟式成穴器在不同傾角下的工作過程，確定了成穴器傾角對(duì)穴孔長度變化與寬度變化的影響。

推土板是一種與土壤發(fā)生接觸響應(yīng)的觸土部件，常用于平整土地，推集泥土等作業(yè)，可針對(duì)各類散體顆粒進(jìn)行工作。劉國敏等[6]對(duì)不同表面形態(tài)的推土板作用下的土壤變化規(guī)律和影響因素進(jìn)行了離散元模擬，解釋了推土板表面形態(tài)對(duì)土壤動(dòng)態(tài)行為的影響。郭志軍等[7]對(duì)推土板的觸土曲面進(jìn)行了不同的觸土曲面準(zhǔn)線設(shè)計(jì)，確定了不同曲面對(duì)工作阻力的影響程度，為推土板的結(jié)構(gòu)優(yōu)化提供了新參考。

國內(nèi)外針對(duì)推土板的的研究多偏向表面形態(tài)、曲線設(shè)計(jì)等外形結(jié)構(gòu)，針對(duì)推土板折彎角度對(duì)推土板性能影響的研究極少，因此對(duì)觸土推板折彎角度的研究很有必要。本文研究推土板在不同折彎角度下與土壤產(chǎn)生動(dòng)態(tài)接觸行為時(shí)土壤的動(dòng)態(tài)響應(yīng)特性，以離散單元法作為研究手段，對(duì)觸土件與土壤顆粒進(jìn)行接觸分析；利用離散單元法分析不同折彎角度的觸土推板工作壓力的變化趨勢、研究觸土推板與土壤顆粒相互作用時(shí)的變化規(guī)律，確定觸土推板在推土作業(yè)中最佳折彎角度，并為后續(xù)的優(yōu)化分析提供理論支撐。

1 觸土推板模型建立

1.1 模型參數(shù)

觸土推板模型的參數(shù)包括長度、寬度、高度、厚度、質(zhì)量、角度等。觸土件與土壤接觸時(shí)會(huì)對(duì)土壤顆粒產(chǎn)生動(dòng)態(tài)響應(yīng)變化，其中角度參數(shù)對(duì)土壤的動(dòng)態(tài)影響不明朗，因此針對(duì)一種具有折彎角度的平面型觸土推板進(jìn)行分析，模型材料選用45號(hào)碳素結(jié)構(gòu)鋼。觸土推板模型如圖1所示。

圖1 觸土推板模型

觸土推板模型的左側(cè)端點(diǎn)為模型末端，右側(cè)端點(diǎn)為模型始端，模型的設(shè)計(jì)將整板分為A和B兩段，兩段長度分別為l1=62.1mm和l2=80mm，兩段的尺寸比接近3∶4。模型的水平方向尺寸為L1，豎直方向尺寸為L2，兩段之間的折彎角度為α，工作時(shí)的切土角度為β。其中：

L1=l1+l2×cosα

(1)

L2=l2×cosβ

(2)

α+β=90°

(3)

觸土推板的折彎角度α一般為0～90°范圍內(nèi)，當(dāng)折彎角度大于45°時(shí)，B段豎直方向的尺寸分量大于B段水平方向的尺寸分量，在固定質(zhì)量參數(shù)條件下豎直方向上過多的尺寸分量會(huì)影響觸土推板的推土能力，造成材料的浪費(fèi)，因此將最大折彎角度設(shè)定為45°。

觸土推板前端土壤的堆積、波動(dòng)等動(dòng)態(tài)行為變化直接受切土角β大小的影響[8]，為確定觸土推板最優(yōu)折彎角度，將折彎角度確定為5個(gè)等級(jí)，模型整體參數(shù)如表1所示。尺寸參數(shù)變化如圖2所示。

表1 模型尺寸參數(shù)表

圖2 模型尺寸參數(shù)變化趨勢

隨折彎角度的增加，L1方向尺寸減小速度逐漸增大，此時(shí)，觸土推板在水平方向的作用范圍縮小的也更快。

1.2 模型運(yùn)動(dòng)分析

觸土推板的工作狀態(tài)為水平直線運(yùn)動(dòng)，在水平進(jìn)給工作時(shí)將土壤散體顆粒推離原始位置。對(duì)觸土推板運(yùn)動(dòng)過程中兩段板上的兩顆粒進(jìn)行分析，其受力和運(yùn)動(dòng)狀態(tài)如圖3所示。

圖3 觸土推板工作狀態(tài)下顆粒運(yùn)動(dòng)趨勢

在不考慮顆粒間作用力及運(yùn)動(dòng)碰撞等情況時(shí)，A段的顆粒受到Fa方向的力，B段的顆粒受到來自Fb方向的力。若A段位置γ角度為0°，B段位置φ角度為90°，模型擁有理想的工作狀態(tài)。

但在實(shí)際運(yùn)動(dòng)狀態(tài)中，觸土推板由于加工誤差、安裝誤差、振動(dòng)、顆粒形狀等因素的影響，推板并不能以理想的運(yùn)動(dòng)狀態(tài)運(yùn)動(dòng)。

單個(gè)土壤顆粒的運(yùn)動(dòng)方向存在不確定性，在多顆粒同時(shí)運(yùn)動(dòng)時(shí)，需要考慮更多的影響因素，包括土壤顆粒間的接觸力、摩擦力、粘性力、范德華力、液橋力等多種作用力。為在更復(fù)雜的工況環(huán)境下進(jìn)行分析，需要使用離散單元法。

2 土壤力學(xué)離散單元模型

2.1 離散單元法顆粒模型原理

離散單元法是把離散體看做具有一定形狀和質(zhì)量的離散顆粒單元的集合。根據(jù)幾何特征可將離散體單元分為顆粒和塊體，土壤屬于一種散體顆粒，在對(duì)土壤顆粒進(jìn)行分析時(shí)常采用顆粒離散單元法。一般情況下認(rèn)為顆粒的運(yùn)動(dòng)是相互獨(dú)立的，只有在發(fā)生相互接觸時(shí)才會(huì)在接觸點(diǎn)處產(chǎn)生力的相互作用，因此顆粒會(huì)表現(xiàn)出離散性，離散性使顆粒在受載和不受載的情形下產(chǎn)生復(fù)雜且不同的運(yùn)動(dòng)。

離散單元法的顆粒模型是將顆粒與顆粒、顆粒與邊界的接觸采用振動(dòng)運(yùn)動(dòng)方程進(jìn)行模擬，模擬過程如圖4所示[9]。

圖4 接觸模型轉(zhuǎn)換成振動(dòng)模型原理圖

圖4中，P1和P2為兩個(gè)顆粒；R1和R2分別為顆粒的接觸半徑；δ為接觸變形量；μn為法向接觸相對(duì)位移；μ為滑動(dòng)位移；Kn、Ks為接觸模型中的法向和切向彈性系數(shù)；ηn、ηs分別為接觸模型中的法向和切向粘滯阻尼系數(shù)；2a為接觸圓直徑。顆粒接觸模型可表示成振動(dòng)模型，振動(dòng)模型由法向振動(dòng)模型、切向振動(dòng)模型和滑動(dòng)模型組合構(gòu)成。

2.2 土壤模型建立

土壤顆粒形狀不規(guī)則且性質(zhì)迥異，采樣地不同，其性質(zhì)參數(shù)會(huì)有很大不同；不同材質(zhì)和成分的土壤具有不同的物理參數(shù)，且在相互作用中的接觸參數(shù)也會(huì)不同，進(jìn)而影響實(shí)際的動(dòng)態(tài)接觸響應(yīng)效果。土壤顆粒的力學(xué)參數(shù)難以直接測量，精確的模擬結(jié)果取決于土壤顆粒模型參數(shù)的設(shè)置和標(biāo)定，合理的土壤模型能夠準(zhǔn)確地反應(yīng)土壤的力學(xué)性質(zhì)[10]。

在土壤模型中將顆粒半徑設(shè)置為4mm，土壤顆粒的粒徑分布呈正態(tài)分布。采用Hertz-Mindlin無滑動(dòng)接觸模型模擬顆粒與顆粒、顆粒與觸土推板間的接觸行為。離散元分析過程需要確定恢復(fù)系數(shù)、摩擦系數(shù)、密度、剪切模量、泊松比等參數(shù)，離散元軟件EDEM參數(shù)設(shè)置如表2所示，土壤為來自吉林省松原市長嶺縣腰井子牧場的砂壤土，經(jīng)參數(shù)標(biāo)定獲得相關(guān)參數(shù)[11]。

表2 EDEM模擬參數(shù)設(shè)置

3 觸土推板流動(dòng)性分析

觸土推板在推土作業(yè)時(shí)，當(dāng)土壤顆粒與推板接觸時(shí)顆粒發(fā)生移動(dòng)，B段斜板位置處的顆粒會(huì)出現(xiàn)向A段平板方向運(yùn)動(dòng)的趨勢；由于推板平動(dòng)作業(yè)時(shí)顆粒運(yùn)動(dòng)趨勢不明顯，因此建立流動(dòng)顆粒-靜物流動(dòng)性分析模型。通過在EDEM中建立邊界完成顆粒堆積，將觸土推板模型放入計(jì)算域內(nèi)作固定設(shè)置，擋板抽離后，顆粒在重力作用下自然下落并與推板發(fā)生碰撞接觸，經(jīng)由推板的阻擋逸出邊界，完成自由狀態(tài)下不同折彎角度推板的顆粒流動(dòng)性分析。流動(dòng)性分析模型如圖5所示。

圖5 流動(dòng)性分析模型

分析模型的總時(shí)間為3s，階段分析步設(shè)置為15%，網(wǎng)格單元設(shè)置為2.5R(R為顆粒半徑)，土壤顆粒由上向下共生成25000個(gè)；當(dāng)模型分析時(shí)間為1s時(shí)，中間部分的擋板抽離，靜止的顆粒開始做自由落體運(yùn)動(dòng)，在運(yùn)動(dòng)一段時(shí)間后開始進(jìn)行流動(dòng)性分析。五種折彎角度的推板分析過程中的運(yùn)動(dòng)速度變化如圖6所示。

圖6 土壤顆粒速度曲線圖

顆粒在1～3s的時(shí)間內(nèi)開始流動(dòng)，因此只針對(duì)1s后的數(shù)據(jù)進(jìn)行分析，并將時(shí)間軸坐標(biāo)歸零處理。由圖6可以看出，顆粒掉落，并在0.2s左右達(dá)到速度峰值，此時(shí)顆粒與觸土推板接觸，速度開始下降，減速的顆粒與逐漸加速的顆粒碰撞沖擊造成速度波動(dòng)，當(dāng)顆粒通過模型的缺口時(shí)開始逸出，避免造成堆積。隨著時(shí)間的增加，顆粒的速度由波動(dòng)變?yōu)榉€(wěn)定下降，45°的觸土推板速度峰值最高，分析結(jié)束時(shí)的速度最低，說明45°推板的顆粒流動(dòng)性更好，角度越大，流動(dòng)性越好。15°觸土推板在速度波動(dòng)周期內(nèi)的整體數(shù)值較低，且分析結(jié)束時(shí)的速度數(shù)值最大，說明推板在這種角度下對(duì)土壤顆粒的流動(dòng)性控制更穩(wěn)定。顆粒數(shù)目的變化曲線如圖7所示。

圖7 顆粒數(shù)目變化曲線圖

由圖7 可以看出，顆粒的數(shù)量變化曲線與速度變化曲線趨勢相同，顆粒的數(shù)目隨時(shí)間的增加而穩(wěn)定下降，但折彎角度不同，顆粒數(shù)目的下降速率不同，分析時(shí)間周期內(nèi)最終留存的顆粒數(shù)目也不同。

圖7中，5°推板與15°推板的留存顆粒數(shù)目相近，說明兩推板對(duì)顆粒數(shù)目的作用效果相近。綜合兩種曲線分析，15°推板的顆粒數(shù)目變化呈線性趨勢，且留存顆粒數(shù)目更多，對(duì)于實(shí)現(xiàn)推土作業(yè)的觸土推板而言，這種綜合性能十分理想。顆粒的運(yùn)動(dòng)狀態(tài)如圖8所示。

圖8 顆粒運(yùn)動(dòng)狀態(tài)圖

由圖8可以看出，顆粒的整體狀態(tài)可劃分為三個(gè)區(qū)域：上端顆粒流動(dòng)區(qū)、中端顆粒擾動(dòng)區(qū)、下端顆粒微動(dòng)區(qū)，三個(gè)區(qū)域很好地反映了顆粒的運(yùn)動(dòng)變化狀態(tài)。

綜合分析認(rèn)為觸土推板的折彎角度為15°時(shí)，土壤顆粒的流動(dòng)性更好，顆粒數(shù)量留存率更高。

4 土壤顆粒動(dòng)態(tài)變化離散元分析

觸土推板的流動(dòng)性分析僅表示土壤顆粒在與推板接觸時(shí)的運(yùn)動(dòng)狀態(tài)，并不能反映土壤顆粒在與推板發(fā)生動(dòng)態(tài)接觸時(shí)的動(dòng)態(tài)響應(yīng)變化，因此需要采用土槽實(shí)驗(yàn)進(jìn)行分析。利用EDEM軟件生成一個(gè)長300mm、寬300mm、高200mm的土槽放置土壤顆粒。土壤顆粒動(dòng)態(tài)變化分析模型如圖9所示。

圖9 土壤顆粒動(dòng)態(tài)變化分析圖

為確定推板在工作時(shí)不同層級(jí)顆粒的運(yùn)動(dòng)變化，在生成顆粒時(shí)將土壤層按深度自下而上分為5層，每層顆粒生成6750個(gè)，共生成33750個(gè)顆粒。設(shè)定推板的前進(jìn)速度為0.05m/s，分析時(shí)間步為20%，網(wǎng)格尺寸為2.5R，模型分析時(shí)間為5s。僅針對(duì)2～5s時(shí)間范圍內(nèi)觸土推板運(yùn)動(dòng)階段進(jìn)行分析，并將此時(shí)的時(shí)間坐標(biāo)歸零。

不同工作狀態(tài)的最終作用效果不同，為測定滿負(fù)載狀態(tài)下推板模型的工作壓力與顆粒層動(dòng)態(tài)變化，將土壤層高度設(shè)置為高于推板高度，土壤層延伸至無限遠(yuǎn)，使用土槽壁模擬運(yùn)動(dòng)阻礙，測定顆粒層的動(dòng)態(tài)變化，如圖10所示。

圖10 土壤層動(dòng)態(tài)變化圖

由圖10可以看出，觸土推板向左移動(dòng)，逐漸對(duì)土壤層顆粒進(jìn)行擠壓，顆粒層受力升高呈現(xiàn)出層級(jí)弧形空間堆積的效果，左上方顆粒層升高越過推板至推板后方，右下方顆粒層僅在觸土推板的破壞范圍內(nèi)呈現(xiàn)出層級(jí)破壞效果，但原有的土壤層堆積形式依舊保持初始狀態(tài)。為直觀地呈現(xiàn)出土壤層動(dòng)態(tài)變化，采用圖11所示的勢能曲線圖進(jìn)行表示。

圖11 土壤顆粒勢能變化趨勢圖

由圖11可以看出，一～三層顆粒在工作時(shí)，勢能隨觸土推板的擠壓而升高；隨著時(shí)間變化，第四層與第五層的顆粒出現(xiàn)勢能先升高后降低的效果，但第五層的勢能變化波動(dòng)與第四層比出現(xiàn)更早。更早的勢能變化說明顆粒越過觸土推板的時(shí)間更早。

對(duì)具有不同折彎角度的觸土推板，前三層基本呈現(xiàn)出折彎角度越小，顆粒勢能的升高越明顯的現(xiàn)象。第三層土壤顆粒在2.8s出現(xiàn)勢能下降的趨勢，此時(shí)顆粒越過了推板。觸土推板的壓力變化曲線如圖12所示。

圖12 觸土推板壓力變化曲線圖

由圖12可以看出，觸土推板在作業(yè)時(shí)，受到的壓力隨時(shí)間增加而增大，并且折彎角度越大，提升至峰值的速率越慢。推板在穩(wěn)定前進(jìn)時(shí)受到的壓力隨時(shí)間產(chǎn)生波動(dòng)，5°推板的壓力變化十分顯著，具有組別中最高的峰值數(shù)據(jù)；15°推板的壓力變化較穩(wěn)定；25°推板的變化較大，并且在分析周期末端產(chǎn)生了劇烈的變化；其余兩條曲線壓力較小，但變化較劇烈。

綜合分析結(jié)果，當(dāng)觸土推板折彎角度為15°時(shí)，推板壓力變化較穩(wěn)定，且此時(shí)觸土推板的工作范圍更大，工作效率更高。

5 結(jié)論

(1)通過對(duì)觸土推板的流動(dòng)性進(jìn)行分析，認(rèn)為觸土推板在15°時(shí)的流動(dòng)性與顆粒數(shù)目留存率更均衡。

(2)通過土槽實(shí)驗(yàn)，對(duì)土壤顆粒進(jìn)行離散元分析，認(rèn)為15°的觸土推板受力更平穩(wěn)。

(3)綜合觸土推板模型尺寸、流動(dòng)性分析評(píng)價(jià)、土槽實(shí)驗(yàn)離散單元法分析，在設(shè)定工況下，15°折彎角度的觸土推板更適合推土工作。