基于LS—DYNA的深松鏟阻力仿真

李霞+張思遠+王維新+湯明軍+郭金龍

摘要： 為降低耕作阻力，揭示深松鏟與土壤之間的關系特性，根據深松鏟切削土壤的工作特點，利用LS-DYNA模擬分析深松鏟切削土壤的過程，獲得切削土壤的應力變化規(guī)律。結果表明，當深松鏟以0.49 m/s的初速度切削、深松深度為250 mm時，單個深松鏟的最大切削阻力為2 230 N，土壤在1.75 s發(fā)生崩裂，達到深松效果。通過試驗測試，單個深松鏟受力約為2 332.5 N，與仿真誤差為4.5%，驗證了仿真的合理性，仿真數據結果可以為實際深松過程提供參考。

關鍵詞： 土壤；切削；有限元法；深松鏟；LS-DYNA模擬分析；阻力仿真

中圖分類號： S222.19 文獻標志碼： A

文章編號：1002-1302（2017）22-0257-03

深松鏟是深松的核心部件，通過深松可以改善土壤三相比，從而促進植物的生長[1]。然而，深松過程中有特別大的阻力。因此，須要探明在深松土壤切削過程中阻力的變化規(guī)律，達到減少阻力的目的。

在以往關于土壤切削的研究中，主要采用試驗方法，但試驗具有受外界環(huán)境因素影響大且經費多、周期長等缺點[2]。近50年來，隨著數值分析模擬技術的發(fā)展，為研究土壤切削提供了一種新的方法。數學模擬切削土壤方法主要包括傳統(tǒng)分析方法、有限單元法（finite element method，簡稱FEM）、離散單元法（distinct element method，簡稱DEM）[3]。近年來，多利用有限元法研究土壤切削問題，通過模擬仿真，可以方便清楚地觀察到土壤和工作部件的情況，例如土壤的運移狀況、破壞位置和工件部件各部分的受力狀態(tài)。這些模擬仿真數據為耕作部件的設計優(yōu)化提供了理論依據。

針對土壤非線性不均勻的特性，本研究采用動力分析方法探討深松鏟對土壤的切削過程，通過顯式非線性動力分析軟件ANSYS/LS_DYNA來進行深松鏟切削土壤的數值模擬，建立深松鏟深松工作過程的有限元仿真模型，對深松過程進行三維數值仿真研究。

1 有限元建模

1.1 土壤模型

ANSYS/LS-DYNA 12.1程序材料庫中193號材料（MAT_DRUCKER_PRAGER）模型是一種彈塑性本構模型。由于這種材料本構模型能準確地體現土壤特性，簡單實用而且能夠得到準確的結果，所以被廣泛地用于模擬巖土材料中，具體參數取值見表1[4]。Mohr-Coulomb屈服條件通常被使用在土壤切削仿真中，但其有不可忽略的缺點：它的屈服面在平面上表現為特殊的等邊六邊形（具有棱角奇異性的不等角），這大大增加了數值計算的困難。與Mohr-Coulomb屈服條件相比，Drucker-Prager屈服準則不僅計算簡單，而且在選擇恰當材料常數的前提下，就可以尋找到與之匹配的Coulomb模型。基于以上比較，本研究采用改進擴展后的Drucker-Prager準則[5]，即：

式中：C為偏應力第3應力不變量；I為單軸拉伸屈服應力與單軸壓縮屈服應力比值；F1為凝聚力；為材料摩擦角；F2為Mises等效應力；F3為平均壓應力；σ1、σ2、σ3為3個方向的主應力；t為應力。

選定土壤的實體模型為800×800×800的立方體，由于模型比較簡單，所以直接在ANSYS中生成。土壤選用solid164單元，這種單元常用于三維顯示結構實體，由8個節(jié)點構成而且支持所有許可的非線性特性，對土壤采用智能網格劃分，選擇土壤模型的12條邊，設置網格數為20，單元形狀采用六面體8個節(jié)點，劃分為9 261個節(jié)點，8 000個單元[6]。

1.2 深松鏟模型

ANSYS在構建復雜實體模型方面比較薄弱，所以本研究利用SolidWorks構建深松鏟模型并通過接口導入到ANSYS中[7]。深松鏟采用 65Mn材料，其密度為7.8×103 kg/m3，彈性模量為2.10萬MPa，泊松比為0.3。由于深松鏟是一個整體，所以采用計算精度較高的四面體單元solid45單元類型，采用智能自由網格劃分方法對深松鏟模型進行網格劃分，設定尺寸級別為5級，為減小應力集中的影響，在深松鏟鏟刃部分進行網格細化，深松鏟劃分為12 724個節(jié)點， 19 046 個單元。在建立土壤-深松鏟數值模擬模型時，對模型進行網格劃分后如圖1所示。

1.3 創(chuàng)建單元

為了給后面定義接觸類型和施加邊界條件作基準，要先創(chuàng)建單元（簡稱PART）。PARTS是1個單元集合體，它是由相同單元類型、材料號和實常數組成的。若要生成PARTS，必須要求選擇所有實體。本案例生成的PARTS結果如圖2所示，其中1代表深松鏟模型，2代表土壤模型。網格劃分數量相當是能夠計算的必要條件，從圖2可以看出，深松鏟模型和土壤模型網格數量相當，這為后續(xù)計算提供了可能。

1.4 定義接觸

ANSYS/LS_DYNA12.1接觸類型多種多樣，而且與接觸類型對應的還有許多可選擇接觸的分析方式，所以在描述深松的過程中，要選擇恰當的類型和參數。深松過程是利用深松鏟來使堅硬的土壤破碎，在計算機模擬系統(tǒng)中，這種工作狀況就是土壤材料單元的失效。在破碎土壤后，土壤材料的外部單元常常無法與深松鏟接觸，為了避免這一情況的發(fā)生，將深松鏟與土壤之間的接觸類型定義為面與面接觸中的侵蝕接觸（eroding surface-to-surface contact，簡稱ESTS）。深松鏟與土壤之間靜摩擦系數為0.60，動摩擦系數為0.32。

1.5 施加邊界條件和載荷

由于實際土壤立方體材料尺寸足夠大，所以為保證土壤結構在整個侵徹過程中不發(fā)生移動，定下表面為固定表面，施加全約束，土壤頂部為自由邊界，其他面均設想為無限大區(qū)域，采用無反射邊界約束[8]。為了節(jié)約計算時間，將土壤的左右側面也固定起來，為土壤前后端面施加無反射邊界條件。由于深松鏟采用的是剛體材料，只須施加在其質心上。對深松施加Z方向的位移約束，限制深松鏟的左右晃動，并約束深松鏟3個方向的轉動自由度。endprint

2 結果與分析

2.1 求解

上述操作完成后，進入求解分析階段。首先設置求解時間，若設置值較大，則會延長計算時間，而且一般物理過程的計算終止時間多為毫秒級，因此，設置求解時間為0.4 s，設置輸出步數為20步。設置完成后，選擇輸出文件類型，輸出關鍵字文件遞交給LS-DYNA求解，并通過后處理程序打開結果文件，進行分析。

2.2 結果分析

載荷大小影響土壤狀態(tài)。隨著外力的變化，土壤會依次經過線彈性—塑性—屈服—硬化這幾個狀態(tài)，在土壤硬化后繼續(xù)加載，則土壤開始失效，直至最終破裂[9]。圖3描述了深松深度為250 mm的土壤深松過程。在深松工作過程中，首先是深松鏟鏟尖接觸土壤，由于受到深松鏟鏟尖的擠壓和剪切作用，土壤發(fā)生變形和破壞。隨著深松工作的繼續(xù)，深松鏟前進并且產生一定振動，此時深松鏟鏟柄也和土壤相互接觸，導致土壤被剪切擠壓的越發(fā)明顯，受到的力也顯著加大，被破壞的面積和進給方向的深度也增加，后面還未與深松鏟接觸的土壤由于受到前面被深松鏟擠壓破壞后土壤的擠壓也逐漸開始變形，最終實現整個深松土壤的變形破壞，達到深松的目的[10]。觀察1個周期深松鏟對土壤的深松過程，整個過程結果輸出20步。在1.25 s時，深松鏟開始和土壤接觸，并以049 m/s的初始速度開始深松土壤，由于受到深松鏟的剪切力作用，土壤狀態(tài)開始發(fā)生改變；在1.5 s時，因受深松鏟的剪切力和擠壓力作用，土壤在沿深松鏟前刀面方向上發(fā)生變形失效，此時切削阻力最大，達到2 230 N。在1.75 s時，土壤發(fā)生崩裂。深松鏟繼續(xù)前進，后面的土壤也被深松。由此可見，通過仿真模擬可以較為直觀地觀察出深松鏟深松工作過程中土壤的狀態(tài)變化。

3 試驗驗證

3.1 試驗準備

為驗證仿真結果的準確性，在新疆石河子市石河子大學教學試驗場二連進行了試驗。在試驗田內選擇1塊地勢比較平坦的未耕作土地，并在試驗前按照振動深松機組試驗方案的要求對試驗區(qū)進行規(guī)劃。確定長60 m（兩側各自留出 10 m 長的調整區(qū)，用于調整拖拉機的工作狀態(tài)）、寬50 m的區(qū)域為振動深松機樣機的試驗場地，并在試驗場地的兩側插上標桿，保證拖拉機在規(guī)劃的區(qū)域內行駛。試驗設備主要有振動深松機樣機1臺，雷沃M800-D型拖拉機1臺，黑龍江省農業(yè)機械工程科學研究院研制的田間機械動力學參數遙測儀1套。

3.2 試驗內容

將振動深松機樣機與拖拉機3點懸掛連接，拖拉機的后置動力輸出軸通過萬向聯軸器與樣機的動力輸入軸相連。拖拉機以慢二擋前進（秒表測得行駛40 m耗時82 s），當拖拉機進入測定區(qū)域時，記錄遙測儀數據，并且重復試驗3次。

3.3 試驗結果

當拖拉機以慢二擋前進（秒表測得行駛40 m耗時82 s）時，重復3次試驗，采集的數據如圖4所示。

由圖4可見，在深松作業(yè)過程中，3次試驗的平均阻力分別為9.35、9.74、8.9 kN，則試驗深松機阻力平均值為9.33 kN；在試驗過程中，深松機共有4鏟同時工作，且工況相同，所以單鏟的力約為F=9.33/4=2.33 kN，試驗力與仿真得到的力的誤差為4.5%，因此利用仿真的方法來研究深松鏟的土壤切削是合理可行的，同時也驗證了該切削模型的正確性。

4 結論

通過建立深松鏟-土壤的切削有限元模型，利用LS-DYNA進行數值分析，結果顯示，當深松鏟以0.56 m/s的初速度切削、深松深度250 mm時，單個深松鏟最大切削阻力為2 230 N，土壤在1.75 s開始發(fā)生崩裂，達到深松效果，并通過試驗驗證了其準確性。仿真結果為優(yōu)化深松鏟結構、提高拖拉機功率的利用率提供了理論依據。

參考文獻：

[1] 李洪文，陳君達，李問盈. 保護性耕作條件下深松技術研究[J]. 農業(yè)機械學報，2000，31（6）：42-45.

[2]王 敏，李萍萍，王紀章，等. 基于LS-DYNA茶園立式旋耕刀具土壤切削過程模擬[J]. 江蘇農業(yè)科學，2015，43（9）：422-425.

[3]Wegner J L，Yao M M，Zhang X.Dynamic wave-soil-structure interaction analysis in the time domain[J]. Computers&Structures，2005，83（7）：2206-2214.

[4]姚仰平，張丙印，俊 高. 土的基本特性、本構關系及數值模擬研究綜述[J]. 土木工程學報，2012，45（3）：127-150.

[5]蘇翠俠，王燕群，蔡宗熙，等. 盾構刀盤掘進載荷的數值模擬[J]. 天津大學學報，2011，44（6）：522-528.

[6] 張 建，王 穎，王淑紅. 基于LS-dyna的油菜靴鞋式開溝器工作阻力動態(tài)仿真[J]. 甘肅農業(yè)大學學報，2015，50（4）：156-159.

[7]蔣紅旗，王繁生. 起重機吊臂結構有限元模態(tài)分析[J]. 農業(yè)機械學報，2006，37（3）：20-22.

[8]張冬梅，張 亞. 基于ANSY/SLS-DYNA頭部形狀對鉆地斜侵徹土壤影響的數值分析[J]. 現代機械，2010（5）：21-22.

[9]劉 輝. 基于LS-DYNA的旱地旋耕刀工作機理研究[D]. 重慶：西南大學，2012.

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