衛(wèi) 韋，張富貴，黃化剛，符德龍，吳雪梅

(1.貴州省煙草公司畢節(jié)市公司,貴州 畢節(jié) 551700；2.貴州大學(xué) a.工程訓(xùn)練中心；b.機(jī)械工程學(xué)院，貴陽 550025)

0 引言

為解決南方山區(qū)多年使用微耕機(jī)進(jìn)行土壤耕整，造成土地耕層淺，淺耕層下形成板結(jié)的犁底層，影響作物根系下扎，水、肥入滲困難等問題[1-2]，研制出了適合南方山區(qū)復(fù)雜地形條件使用而且耕深可以打破犁底層的深耕機(jī)—小型履帶式鏈?zhǔn)缴罡麢C(jī)，其耕深可以達(dá)到25cm。ANSYS/LS-DYNA軟件具有同時利用ANSYS軟件強(qiáng)大的前后處理功能和LS-DYNA軟件非線性求解功能的特點，能夠有效地進(jìn)行切削過程的有限元模擬分析[3-6]。近3年來，越來越多的專家學(xué)者采用這個軟件對各種刀具切削土壤進(jìn)行仿真模擬。例如，2017年，中交疏浚技術(shù)裝備國家工程研究中心有限公司的楊晨等人采用LS-DYNA軟件對絞刀切削土壤進(jìn)行仿真模擬[7]；2016年，西南大學(xué)的韓煜杰采用LS-DYNA軟件對旋耕刀切削土壤進(jìn)行仿真模擬[8]。2015年，華南農(nóng)業(yè)大學(xué)的齊龍等人采用該軟件對松土刀切削土壤進(jìn)行有限元仿真[9]。本文基于ANSYS/LS-DYNA軟件，采用LSTC公司開發(fā)的LS-PREPOST專用后處理器直接讀取LS-DYNA計算生成結(jié)果，進(jìn)行數(shù)據(jù)匯整及二次運(yùn)算，進(jìn)行了所研制深耕機(jī)刀具切削土壤的顯式動力學(xué)分析[10-11]。經(jīng)過對切削過程的研究，對單個刀具切削土壤過程進(jìn)行分析，旨在揭示單個刀具-土壤的工作機(jī)理。同時，對其切削力以及切削深度進(jìn)行了分析，驗證了所研制出的深耕機(jī)可以達(dá)到25cm的設(shè)計耕深要求，同時也為分析切削力及擬為深耕機(jī)刀具與土壤相互作用的數(shù)值模擬方面的研究提供參考。

1 材料與方法

1.1 土壤本構(gòu)理論

通常情況下，在給土體施加外載荷后其狀態(tài)與所施加載荷的數(shù)值有著密不可分的關(guān)系。當(dāng)所施加載荷小到足以忽略時，土體就會呈現(xiàn)出線彈性。如果外載荷不間斷的持續(xù)加大時，土體所受應(yīng)力就會高于彈性極限，此時土體的應(yīng)力應(yīng)變關(guān)系就會發(fā)生變化，不再是之前所呈現(xiàn)出的理想彈性狀態(tài)，土體的一些部分開始向塑性狀態(tài)轉(zhuǎn)變。判斷是否轉(zhuǎn)變?yōu)樗苄誀顟B(tài)的條件稱之為屈服條件。但是，有些土壤類型，即使達(dá)到屈服條件以后也并不發(fā)生破壞，而且繼續(xù)對于加載依然能夠承受，此時應(yīng)力和應(yīng)變之間的關(guān)系就會發(fā)生改變，不再是之前所呈現(xiàn)出的線彈性關(guān)系，而是開始強(qiáng)化，應(yīng)力和應(yīng)變之間的關(guān)系變?yōu)榉蔷€性關(guān)系[12]。本文的土壤模型所采用的屈服準(zhǔn)則是Mohr-Coulomb條件(莫爾-庫倫屈服準(zhǔn)則)和Drucker-Prager條件。

1)Mohr-Coulomb條件(見圖1)。Coulomb給出了土體每一受力平面上所受的切應(yīng)力大小達(dá)到極限值時，土體就逐步踏入屈服極限狀態(tài)，其表達(dá)式為

τ=c+σntanφ

(1)

式中φ—土壤內(nèi)摩擦角；

σn—每一受力平面上的正應(yīng)力；

n—受力面外法線方向；

c—粘聚力。

圖1 Mohr-Coulomb條件Fig.1 Mohr-Coulombconditions

后經(jīng)Mohr.優(yōu)化，推廣到一般情況。Mohr-Coulomb條件也可以寫成以下形式

(2)

在三維情況下，該條件可表示為

(3)

式中I1—應(yīng)力張量的第1個不變量；

Mohr-Coulomb條件在主應(yīng)力空間為一棱椎體。

2)Drucker-Prager條件。Drucker和Prager提出可對Mohr-Coulomb條件的改進(jìn)。其表達(dá)式為

(4)

(5)

(6)

實際上，I1=σ1+σ2+σ3是發(fā)生體積變化時的靜水壓力狀態(tài)，如果比較大的靜水壓力作用于金屬材料上時，由于其體積發(fā)生了非常小的改變，故不考慮在內(nèi)，式(4)中的首項就刪掉。對于土體來講，多數(shù)時候下的體積改變是顯而易見的。

1.2 數(shù)值模擬中所采用的土壤本構(gòu)模型

為使刀具切削土壤數(shù)值仿真模擬更接近現(xiàn)實，所以土壤本構(gòu)關(guān)系的選用尤為重要。本文選擇的土壤模型材料為MAT147(MAT_FHWA_SOIL)，其應(yīng)用的是修正的Drucker-Prager塑性模型，如圖2所示。Mohr-Coulomb與Drucker-Prager屈服表面的結(jié)合使模擬結(jié)果更接近真實破壞平面(同時避免了Mohr-Coulomb屈服準(zhǔn)則的奇異性)[13]。

圖2 修正的Drucker-Prager塑性模型Fig.2 Modified Drucker-Prager plastic model

Abbo和Sloan給出了修正的Drucker-Prager塑性模型。當(dāng)剪切力微乎其微時，屈服面就成為平滑的面，同時壓力軸成直角關(guān)系。改善后的Mohr-Coulom屈服準(zhǔn)則式為

(7)

式中P—壓力；

φ—內(nèi)摩擦角；

J2—應(yīng)力偏張量的第2個不變量；

K(θ)—張量平面角的函數(shù)；

c—粘聚力。

當(dāng)Ahyp=0時，式(7)為標(biāo)準(zhǔn)Mohr-Coulom屈服面。針對這種仿真來說，Ahyp取值就要低于ccosφ，參照Lewis文獻(xiàn)，通常依據(jù)下式選取，則

(8)

2 深耕刀具切削土壤模型確定

2.1 具切削土壤初始位置模型

選擇土壤MAT147的本構(gòu)模型建立了土壤模型，通過SolidWorks搭建了深耕刀具的三維模型,將其另存為.x-t格式的文件,再導(dǎo)進(jìn)ANSYS—LS-DYNA中。圖3為刀具與土壤初始位置模型。

圖3 刀具切削土壤初始位置模型Fig.3 Tool cutting soil initial position model

2.2 ANSYS—LS-DYNA中模型參數(shù)設(shè)定

現(xiàn)將刀具切削土壤模型的模擬過程中相關(guān)參數(shù)的設(shè)定如下：

1)設(shè)定ANSYS—LS-DYNA中的SOLID164作為土壤和刀具的單元類型，為防止在網(wǎng)格劃分過程中出現(xiàn)沙漏能會影響結(jié)果的準(zhǔn)確性，本文采用的是全積分算法。將刀具用手動來劃分有限元網(wǎng)格， 劃分單元數(shù)為3 880個；將土壤用手動來劃分有限元網(wǎng)格，劃分單元數(shù)為55 000個。圖4為刀具切削土壤模型的網(wǎng)格劃分。

圖4 刀片切削土壤模型的網(wǎng)格劃分Fig.4 Mesh generation of blade cutting soil model

2)通過實際情況的分析，選取了MAT-147塑性模型和MAT-GIRID剛性體材料模型分別作為本文仿真模擬研究的土壤模型和刀具模型。表1和表2是土壤模型及刀具模型的材料參數(shù)。

表1 土壤材料參數(shù)Table 1 Soil material parameters

表2 刀具材料特性Table 2 Tool material characteristics

3)將刀具和土壤分別設(shè)置為接觸Part和目標(biāo)Part，把刀具與土壤之間的接觸類型設(shè)置為面面侵蝕接觸，此接觸選項能夠做到如果在出現(xiàn)表面單元失效的時候，程序能夠自動定義新的接觸面。

4)設(shè)置刀具沿X方向平動切削土壤，大小為1.1m/s。

5)刀具切削土壤模型的約束設(shè)置。在本文中，需要對土壤及其刀具加上相關(guān)的約束條件；將土壤底部設(shè)置為全約束，目的是使得土壤模型在施加切削力后可以真實的體現(xiàn)土壤的變形狀態(tài)。由于要滿足刀具沿x軸方向平動切削的要求，所以對刀具z軸方向進(jìn)行了全約束，圖5為刀具切削土壤模型施加約束。

圖5 刀具切削土壤模型施加約束Fig.5 Tool cutting soil model imposing constraints

2.3 切削模型k文件輸出及其關(guān)鍵字修改

將在SolidWorks中建好的刀具模型導(dǎo)入ANSYS-LS-DYNA前處理器中，并在其中建立好土壤模型，通過選定單元類型、定義刀具和土壤的材料模型、劃分網(wǎng)格、指定接觸類型、對邊界條件進(jìn)行設(shè)置、定義刀具初始速度以及設(shè)置求解控制設(shè)置等一系列操作后生成K文件；然后，根據(jù)模擬研究的要求對生成的K文件進(jìn)行關(guān)鍵字的更改；將修改好的K文件保存并導(dǎo)進(jìn)ANSYS-LS-DYNA求解器中進(jìn)行運(yùn)行求解。

2.4 ANSYS-LS-DYNA求解器求解

將K文件進(jìn)行修改或替換后調(diào)入ANSYS-LS-DYNA求解器中進(jìn)行運(yùn)行計算，整個求解過程所用的時間大約7h。

2.5 結(jié)果后處理

用ANSYS-LS-DYNA自帶的后處理器LS-PREPOST導(dǎo)入求解器運(yùn)行后生成的d3plot結(jié)果文件，進(jìn)行可視化處理觀察動態(tài)仿真過程并進(jìn)行分析。

3 結(jié)果與討論

3.1 刀具切削土壤過程分析

在仿真過程中，刀具沿直線切削土壤。刀具以速度1.1m/s沿前進(jìn)方向(即仿真中X方向)對土壤進(jìn)行切削。圖6為1個周期的作業(yè)過程。觀察1個切削周期，發(fā)現(xiàn)切削過程中刀具刀尖首先接觸土壤，該處土壤受刀具切削力作用發(fā)生剪切破壞從而凹陷變形；隨著刀片繼續(xù)前進(jìn)，刀具對土壤的接觸面也隨之增大，擠壓變形也更加劇烈，刀具不間斷地進(jìn)行土壤切削，對土壤的破壞效應(yīng)也急劇加大，然后還未與刀具接觸的土壤受前面土壤的擠壓也依次發(fā)生變形效應(yīng)。

(a) 0.03s

(b) 0.06s

(c) 0.09s圖6 切削過程組圖Fig.6 Cutting process diagram

3.2 切削力分析

刀具切土壤的過程是非常復(fù)雜的，因而過程是不易控制的。專家學(xué)者經(jīng)大量研究發(fā)現(xiàn)：土壤的破壞一般情況下是剪切破壞。因為土壤本身特性，自身有強(qiáng)度的存在，在刀具剛開始與土壤接觸發(fā)生相互作用時，刀具上的切削力小于對土壤破壞的強(qiáng)度，此時土壤只會產(chǎn)生擠壓變形，隨著切削過程的進(jìn)行切削力也會隨之增強(qiáng)。當(dāng)切削力不斷增大突破土壤的破壞強(qiáng)度時，土壤開始遭到瓦解，刀具所受切削力此時不再繼續(xù)增加；當(dāng)所有土壤都被刀具切削過后，刀具切削力立即減下來并漸漸趨于水平狀態(tài)。

如圖7是刀具切削土壤的一個周期的切削力變化情況，選擇了4個時間段的Von Mises Stress圖像。其等效應(yīng)力大多數(shù)情況下是聚集在土壤變形最大處，出現(xiàn)這種情況的原因是由于刀具隨著切削土壤的深度增大，產(chǎn)生的切削力自然也隨之增強(qiáng)導(dǎo)致的。土壤變形最大處成為應(yīng)力集中區(qū)域。當(dāng)?shù)毒叻謩e切削到0.03、0.06、0.09、0.12s時，對應(yīng)的最大等效應(yīng)力分別為1.831e-06GPa，1.258e-06GPa，6.955e-07GPa，6.535e-07GPa。

(a)

(b)

(c)

(d)圖7 1個周期內(nèi)刀具切削土壤的Von Mises Stress等效應(yīng)力云圖Fig.7 The Von mises stress equivalent stress clouds of cutting soil in a cycle soil in a cycle

根據(jù)不同時刻刀具切削土壤的等效應(yīng)力變化規(guī)律用曲線畫出，如圖8所示。觀察發(fā)現(xiàn)：在刀具剛進(jìn)入土壤的一段時間，等效應(yīng)力隨時間的延長而持續(xù)加大，隨之刀具一步步將土壤結(jié)構(gòu)損壞，然后等效應(yīng)力就下降。在切削到0.09s后，刀具切削土壤的切削力漸漸處于水平狀態(tài)，變化幅度不大，等效應(yīng)力也慢慢保持不變。這時，刀具所受等效應(yīng)力主要是土壤被破壞后對刀具的阻力、粘滯力以及土壤重力。

圖8 等效應(yīng)力隨時間變化曲線Fig.8 Equivalent stress versus time curve

依照等效應(yīng)力隨時間變化曲線圖可以分析出刀具上該單元或者周邊單元切削還未切的土壤時的應(yīng)力數(shù)值。分析刀具切削土壤的受力情況，為刀具參數(shù)的選取以及優(yōu)化提供了參考依據(jù)。

在土壤模型上選取A(H31384)、B(31385)、C(30013)單元如圖9所示。在整個作業(yè)過程中，A、B、C3個單元Y方向(即深耕的深度方向)位移由小變大的變化情況如圖10所示。

在土壤模型中取點A、B、C單元，由圖10觀察到它們的Y向位移。在約0.09s時，A單元土壤的Y向(即深耕的深度方向)位移約為0.25m，說明本文設(shè)計的深耕機(jī)可以達(dá)到耕深25cm的深度。

圖9 ABC 3個單元Fig.9 Three-units of the ABC

圖10 ABC 3單元Y向位移變化Fig.10 The Y-direction displacement of the three-units of the ABC

4 結(jié)論

應(yīng)用有限元分析理論，應(yīng)用SolidWorks建立刀具模型，導(dǎo)進(jìn)ANSYS/LS-DYNA軟件中，再創(chuàng)建土壤的幾何模型；然后運(yùn)用軟件對單個刀具切削土壤過程進(jìn)行分析，旨在揭示單個刀具-土壤的工作機(jī)理。同時，對其切削力以及切削深度進(jìn)行了分析，驗證了所研制出的深耕機(jī)可以達(dá)到25cm的設(shè)計耕深要求，同時也在分析切削力及擬為深耕機(jī)刀具與土壤相互作用的數(shù)值模擬方面的研究提供了參考。