摘要：為解決傳統(tǒng)深松機(jī)具觸土部件破土困難、耕作阻力大等問(wèn)題，基于離散元法分析以鯊魚(yú)背鰭設(shè)計(jì)的深松鏟耕作行為，以期減小深松鏟作業(yè)阻力和能耗。設(shè)計(jì)仿鯊魚(yú)背鰭鏟尖試樣，并與無(wú)鰭深松鏟試樣進(jìn)行性能對(duì)比。建立離散元模型，探究垂直貫入耕作方式下不同結(jié)構(gòu)的土壤阻力；分析深松過(guò)程中土壤擾動(dòng)、土壤運(yùn)動(dòng)以及帶鰭深松鏟的受力情況。結(jié)果表明，在帶鰭深松鏟的工作過(guò)程中，各層的擾動(dòng)程度從大到小依次為耕作層、犁底層和心土層；仿真試驗(yàn)分析相同深松深度條件下，不同工作速度對(duì)帶鰭深松鏟的土壤顆?；旌锨闆r、土層擾動(dòng)效應(yīng)以及耕地阻力的影響，得出3.60～4.33 km/h的最佳工作速度。結(jié)合田間試驗(yàn)可知，當(dāng)深松深度相同時(shí)，帶鰭深松鏟的真實(shí)耕作阻力（1 904 N）與離散元模擬值（1 941 N）的偏差為1.9%。所提出的帶鰭式深松鏟機(jī)能夠有效地破開(kāi)犁底，為我國(guó)東北地區(qū)耕地科學(xué)研究提供借鑒。

關(guān)鍵詞：深松鏟；離散元；土壤擾動(dòng)；耕作阻力；仿生法

中圖分類號(hào)：S222

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

文章編號(hào)：2095-5553 （2024） 12-0015-06收稿日期：2023年11月23日

修回日期：2024年7月3日

*基金項(xiàng)目：國(guó)家重點(diǎn)研發(fā)計(jì)劃課題（2022YFD1500605—03）

第一作者：陳菁，男，1987年生，遼寧錦州人，碩士，高級(jí)農(nóng)藝師；研究方向?yàn)檗r(nóng)業(yè)機(jī)械化工程。E-mail：169333333@qq.com

通訊作者：張旭，男，1979年生，沈陽(yáng)人，碩士，正高級(jí)工程師；研究方向?yàn)檗r(nóng)業(yè)機(jī)械化工程。E-mail：lnami@126.com

Simulation analysis of cultivation behavior of shark-fin-like deep pine shovel based on discrete elements

Chen Jing^{1， 2}， Shan Hongdao^{1， 2}， Wei Chuansheng^{1， 2}， Zhang Xu^{1， 2}

（1. Liaoning Provincial Institute of Agricultural Mechanization， Shenyang， 110161， China； 2. Key Laboratory of Remanufacture and Innovation of Agricultural Machinery and Equipment in Liaoning Province， Shenyang， 110161， China）

Abstract： In order to solve the problems of difficulty in breaking the soil and high tillage resistance of the touching part of the traditional deep loading equipment， the tillage behavior of the deep loading shovel designed with shark dorsal fin was analyzed based on the discrete element method， so as to reduce the operational resistance and energy consumption of the deep loading shovel. In this paper， a shark-like dorsal fin shovel tip specimen was designed， and the performance was compared with that of an undesigned deep pine shovel without fins. A discrete element model was established to solve the vertical penetration soil resistance of different structures， soil disturbance， soil movement and the force of the deep-pine shovel with fins were analyzed during the deep-pine process. The results show that during the working process of the deep-pine shovel with fins， the disturbance degree of each layer is the tillage layer， the plow subsoil layer and the heart soil layer in descending order， based on different working speeds， the mixing of soil particles， the soil disturbance effect and the tillage resistance of the deep-pine shovel with fins are investigated in the same deep-pine depth， and the optimal working speeds of 3.60-4.33 km/h are obtained. Combined with the field test， it can be seen that the deviation between the real tillage resistance （1 904 N） and the discrete element simulation value （1 941 N） of the band-fin deep-pine shovel is 1.9% when the deep-pine thickness is the same. The experimental results proves that the deep-pine shovel with fins proposed in this paper can effectively break the plow bottom and provide a reference for scientific research on cultivated land in the northeast region of China.

Keywords： subsoiler; discrete element; soil disturbance; tillage resistance; biomimetic method

0 引言

深松是一種能有效緩解耕地緊實(shí)度、提高產(chǎn)量的重要保護(hù)性耕作技術(shù)，在國(guó)內(nèi)外已被廣泛采用。然而，當(dāng)前農(nóng)田作業(yè)能耗較高，約占農(nóng)田作業(yè)總量的40%～60%^［1］。深松工作阻力大、能量損失高及作業(yè)效率低，這是在深松領(lǐng)域中迫切需要解決的問(wèn)題。當(dāng)前，深松減阻的方法主要有：電子滲透減阻，分層深松減阻，振動(dòng)深松減阻以及結(jié)構(gòu)優(yōu)化減阻^［2］。目前的脫層與振動(dòng)深松仍然面臨著能耗高、結(jié)構(gòu)復(fù)雜的難題。

大量的研究表明，從結(jié)構(gòu)上對(duì)深松鏟進(jìn)行優(yōu)化，通過(guò)改變破土刃的弧度，以及鏟柄和鏟尖的幾何形狀，可以降低耕作的阻力，減少能量的損耗^［3-5］。近幾年來(lái)，在農(nóng)業(yè)機(jī)械設(shè)計(jì)中，越來(lái)越多地運(yùn)用了仿生學(xué)的方法。田鼠、野豬及狗獾等土壤挖掘動(dòng)物的幾何特征被用于對(duì)深松鏟進(jìn)行優(yōu)化設(shè)計(jì)，并且對(duì)其減阻有較好的效果^［6-8］。

農(nóng)機(jī)具耕作過(guò)程的模擬仿真是一項(xiàng)復(fù)雜的工作，而離散元法從微觀的角度對(duì)仿真系統(tǒng)中的每個(gè)散體顆粒進(jìn)行分析，能夠很大程度提高模擬準(zhǔn)確性，使仿真結(jié)果更加接近于實(shí)際情況。有學(xué)者對(duì)鯊魚(yú)魚(yú)鰭的減阻效果進(jìn)行了CFD仿真分析，試驗(yàn)證明了鯊魚(yú)背鰭在鯊魚(yú)運(yùn)動(dòng)過(guò)程中受到的阻力小于鯊魚(yú)胸鰭^［9］。張智泓等^［10］以砂魚(yú)蜥蜴的頭顱作為仿生學(xué)模型，通過(guò)逆向工程方法，獲取了其特有的形狀，并將其用于深松鏟頭的設(shè)計(jì)，從而有效地解決了常規(guī)深松機(jī)具的觸土部件難以破土和耕作阻力較大的難題。Zhang等^［11］對(duì)粘性土的力學(xué)特性進(jìn)行了數(shù)值模擬，模擬結(jié)果能夠反映粘性土在壓力下的變形破壞，但由于是2D模型，其適用范圍尚不明確。于建群等^［12］采用2D CAD技術(shù)，分別構(gòu)建了帶芯犁頭的犁頭犁、播種肥料犁頭與土粒的2D離散元分析模型，并采用2D方法，對(duì)犁頭在不同工況下的工作阻力進(jìn)行了數(shù)值模擬，并對(duì)其進(jìn)行了試驗(yàn)驗(yàn)證，結(jié)果表明，模擬結(jié)果與實(shí)測(cè)結(jié)果之間的相對(duì)誤差為10%～20%。因此可利用SolidWorks進(jìn)行深松鏟建模，采用離散元法對(duì)觸土部件工作性能進(jìn)行評(píng)價(jià)和分析逐漸成為研究熱點(diǎn)，是一種用于模擬并分析散體介質(zhì)系統(tǒng)動(dòng)力學(xué)行為的數(shù)值方法。

本文以帶鰭和不帶鰭式兩種深松鏟為研究對(duì)象，利用離散元分析軟件EDEM對(duì)松鏟的耕作過(guò)程進(jìn)行數(shù)值模擬，從土壤擾動(dòng)、土壤顆粒運(yùn)動(dòng)、工作阻力工多個(gè)維度對(duì)深松鏟的耕作行為進(jìn)行分析。將模擬與實(shí)測(cè)結(jié)果進(jìn)行對(duì)比，以驗(yàn)證離散元方法的有效性；利用離散元法模擬土壤顆粒分別與仿生深松鏟和無(wú)背鰭鑿式深松鏟相互作用后的力學(xué)行為及顆粒運(yùn)動(dòng)情況，為深松作業(yè)的節(jié)能減阻提供理論基礎(chǔ)和設(shè)計(jì)依據(jù)。

1 仿真模型建立

1.1 深松鏟模型

根據(jù)JB/T 9788—2020《深松鏟和深松鏟柄》，基于SolidWorks對(duì)深松鏟結(jié)構(gòu)進(jìn)行建模，結(jié)構(gòu)由深松鏟尖和深松鏟柄組成，其二維平面結(jié)構(gòu)如圖1所示。其中，鏟柄長(zhǎng)度為152 mm，鏟寬度為60 mm，鏟柄高度為686 mm，鏟柄工作刃高度為320 mm，鏟柄外圓弧半徑為321 mm，鏟鰭邊長(zhǎng)為70 mm，鏟鰭內(nèi)角為60°，鏟柄刃角為23°。

1.2 土壤顆粒模型

由于長(zhǎng)期耕作下的農(nóng)田會(huì)形成分層結(jié)構(gòu)，該分層結(jié)構(gòu)分別為耕作層、犁底層和心土層。根據(jù)土壤顆粒一般形狀（核狀、條狀、片狀和團(tuán)聚體等結(jié)構(gòu)）基于EDEM離散元軟件對(duì)土壤分層結(jié)構(gòu)中的土壤顆粒進(jìn)行建模，圖2是土壤顆粒仿真模型，模型顆粒大小為6 mm^［13］。

1.3 土壤離散元接觸模型參數(shù)確定

在DEM計(jì)算中，合適的接觸模式是保證計(jì)算結(jié)果精度的關(guān)鍵。研究對(duì)象為農(nóng)田土壤顆粒，考慮到農(nóng)田土壤顆粒之間的黏聚作用，采用Hertz-Mindling with JKR 模型作為耕作層和心土層的土壤顆粒間接觸模型^［14］。由于犁底土顆粒比耕層、心土顆粒更硬，并具有很強(qiáng)的“液橋”黏著力，故選擇一種基于Hertz-Mindling and Bond模型來(lái)描述犁底層土壤顆粒間的接觸行為。表1為模型的特定參數(shù)^［15］。

1.4 虛擬土槽仿真模型建立

利用SolidWorks開(kāi)發(fā)仿鯊魚(yú)背鰭深松鏟的三維構(gòu)造模型，將其導(dǎo)入EDEM離散元仿真軟件。模擬模型固定時(shí)間步長(zhǎng)為8.27×10^-6 s，Rayleigth時(shí)間步長(zhǎng)為4.13×10^-5 s。為分析深松鏟對(duì)不同深度土層土壤的運(yùn)動(dòng)及擾動(dòng)的差異，設(shè)定耕作層、犁底層、心土層3層土壤模型。在EDEM模型設(shè)置中，結(jié)合實(shí)際耕作情況將土層自上而下設(shè)置耕作層厚度為150 mm，犁底層厚度為140 mm，心土層厚度為140 mm。根據(jù)文獻(xiàn)［15］可知深松的最大工作深度為200～350 mm，本文研究的最大工作深度為300 mm。圖3為EDEM中建立的虛擬土槽仿真模型。

2 離散元仿真

2.1 深松鏟耕作仿真分析

為探索仿鯊魚(yú)背鰭深松鏟的最佳工作速率，對(duì)不同轉(zhuǎn)速下深松鏟的工作效果進(jìn)行分析。深松仿真時(shí)間為1.5～4.4 s，深松深度為300 mm，選取作業(yè)速度2.88 km/h（仿真時(shí)間為1.5～2.0 s）、3.6 km/h（仿真時(shí)間為2.0～2.5 s）、4.33 km/h（仿真時(shí)間為2.5～3.0 s）、5.66 km/h（仿真時(shí)間為3.0～3.5 s）。通過(guò)對(duì)深松鏟土體的混合狀況、土層的擾動(dòng)效應(yīng)、耕地阻力等因素的分析，確定最優(yōu)工作速率。由圖4可知，深松鏟對(duì)土壤土層間的顆?；旌狭侩S作業(yè)速度的增加而逐漸增多。當(dāng)作業(yè)速度為3.6 km/h時(shí)，各層土壤顆粒間開(kāi)始出現(xiàn)混合現(xiàn)象；當(dāng)作業(yè)速度為2.88 km/h時(shí)，各層土壤顆粒間混合不明顯；當(dāng)作業(yè)速度為4.33 km/h時(shí)，耕作層與犁底層之間的土壤顆粒混合現(xiàn)象較為明顯；在深松鏟作業(yè)速度為5.66 km/h時(shí)，土壤顆?；旌犀F(xiàn)象大量出現(xiàn)。

根據(jù)仿真結(jié)果可知，深松鏟作業(yè)速度達(dá)到4.33 km/h后，由于高速碰撞使土體在碰撞過(guò)程中產(chǎn)生了大量的松散顆粒，并伴隨著深松鏟運(yùn)動(dòng)向其他土層遷移，使各土層間的土粒混雜程度大大提高。當(dāng)深松鏟作業(yè)速度為4.33～5.66 km/h，可獲得最優(yōu)的混土效果，且能保證工作效率。

從圖5中可以看出，隨著工作速率的提高，深松鏟對(duì)土壤的擾動(dòng)程度也隨之提高。在鏟尖、鏟翼、鏟柄、鏟刃等工具的剪切、擠壓下，深松鏟鏟身側(cè)的土體會(huì)隨著鏟頭的推進(jìn)而向上隆起，從而擴(kuò)大土層的擾動(dòng)范圍。在2.88 km/h的速度下，深松鏟不會(huì)對(duì)土壤造成顯著的破壞；在3.6 km/h的速度下，對(duì)地層的擾動(dòng)較大；當(dāng)速度為4.33～5.66 km/h時(shí)，深松鏟對(duì)犁底、耕層的干擾逐漸加強(qiáng)，變化逐步趨于平穩(wěn)。

通過(guò)EDEM仿真計(jì)算，研究帶鰭深松鏟在相同耕作深度不同耕作速度下的耕作阻力。圖6為不同耕作速度下的耕作阻力對(duì)比。當(dāng)耕作速度為2.88 km/h時(shí)，深松鏟耕作阻力均值為781 kN；當(dāng)耕作速度為3.6 km/h時(shí)，深松鏟耕作阻力均值為1 351 kN；當(dāng)耕作速度為4.33 km/h時(shí)，深松鏟耕作阻力均值為1 534 kN；當(dāng)耕作速度為5.66 km/h時(shí)，深松鏟耕作阻力均值為1 941 kN。根據(jù)仿真計(jì)算可知，隨著工作速度的增大，土壤對(duì)鏟頭、鏟柄和鏟鰭的壓力逐漸增大，從而產(chǎn)生更大的耕作阻力。

根據(jù)仿真分析結(jié)果可知，深松鏟耕作速度越大，耕作阻力越大。深松鏟耕作速度的變化使土壤深松效果不同，結(jié)合圖5和圖6可知，當(dāng)深松鏟耕作速度在3.6～4.33 km/h，各土層之間的土粒摻雜很少，且對(duì)土層的擾動(dòng)作用顯著，試驗(yàn)所得深松鏟耕作阻力處于深松鏟最優(yōu)作業(yè)速度范圍。

2.2 深松效果對(duì)比

深松鏟耕作時(shí)應(yīng)在達(dá)到更大松土面積、更小溝槽寬度和地表平整度的前提下，減少土壤擾動(dòng)。

2.1.1 松土面積

為進(jìn)一步分析兩種深松鏟的松土面積，通過(guò)仿真分析仿形繪制深松鏟在作業(yè)過(guò)程中土壤的擾動(dòng)輪廓，記錄關(guān)鍵點(diǎn)的坐標(biāo)數(shù)據(jù)即可得到兩種深松鏟作業(yè)時(shí)的松土效果圖。圖7為兩種深松鏟松土效果對(duì)比圖，松土面積S可表示為

S=∫^b_af（x）+∫^d_cg（x）+L·H （1）

式中： a、b——深松鏟左側(cè)松土區(qū)域邊界點(diǎn)；

f（x）——深松鏟左側(cè)松土區(qū)域邊界近似函數(shù)；

c、d——深松鏟右側(cè)松土區(qū)域邊界點(diǎn)；

g（x）——深松鏟左側(cè)松土區(qū)域邊界近似函數(shù)；

x——深松鏟松土區(qū)域任意點(diǎn)；

L——深松鏟炳寬度，mm;

H——深松深度，mm;

S——松土面積，mm²。

由圖7可知，在相同作業(yè)深度條件下，無(wú)背鰭鑿式深松鏟最大松土寬度W₁小于仿鯊魚(yú)背鰭深松鏟松土寬度W₂，通過(guò)計(jì)算可以得知仿鯊魚(yú)背鰭深松鏟的松土面積大于無(wú)背鰭鑿式深松鏟的松土面積，從而表明仿鯊魚(yú)背鰭深松鏟的松土量更大，更有利于土壤深松。

2.2.2 土壤擾動(dòng)效果

對(duì)土壤的干擾效應(yīng)是對(duì)深松操作進(jìn)行評(píng)估的一個(gè)重要指標(biāo)。在兩種深松鏟作業(yè)時(shí)，土壤顆粒的受力情況如圖8所示。其中，上層土壤顆粒顏色較深，表示土壤顆粒受力較大；下層土壤顆粒顏色最淺，表示受力較??；中層土壤顆粒顏色較淺，表示受力大小介于上層和下層二者之間?？梢钥闯觯瑤挼纳钏社P與不帶鰭的深松鏟相比，對(duì)每一土層中的土粒的撞擊更少，對(duì)土層的損害也更小，帶鰭式深松鏟對(duì)耕作層和心土層土壤顆粒的沖擊比犁底層要大。圖8（a）中土壤中部的犁底層顆粒深色部分較少，也就是帶鰭式深松鏟對(duì)犁底層土壤起到了滑切作用，不僅可以促進(jìn)土壤中的雨水下滲，還可以降低深松鏟對(duì)土層破壞的程度。與此同時(shí)，從圖8（b）可以看出，與無(wú)鰭式深松鏟相比，帶鰭式深松鏟可以有效地降低對(duì)各個(gè)土壤層的土壤的混合量，達(dá)到不亂土層的效果，減少了較低養(yǎng)分含量的芯土層的上移，因此可以獲得更好的深松效果。

2.2.3 地表平整度和溝槽寬度

為測(cè)定深松鏟對(duì)地表平整度與溝槽寬度的影響，利用離散元軟件中的Clipping模塊對(duì)仿真土槽進(jìn)行區(qū)域劃分。提取該區(qū)域穩(wěn)定后的坐標(biāo)值，進(jìn)行土壤輪廓標(biāo)定，繪制出深松前后的地表線，并在深松后的地表線上過(guò)最高點(diǎn)作水平線為基準(zhǔn)線。如圖9所示，在整個(gè)深松鏟深松寬度L_Y范圍以d₀=30 mm為等分標(biāo)記測(cè)點(diǎn)，在兩種深松鏟作業(yè)前后測(cè)定土壤輪廓曲線至水平基準(zhǔn)線的垂直距離^［16］。第i次測(cè)量的地表平整度D_i計(jì)算如式（2）所示。

D_i=∑n/j=1（y_ij－y_i）²/n－1 （2）

式中： y_i——第i次測(cè)量深松前后地表至水平基準(zhǔn)線垂直距離，mm;

y_ij——第i次測(cè)量中第j個(gè)等分點(diǎn)深松前后地表至水平基準(zhǔn)線垂直距離，mm;

n——測(cè)試點(diǎn)個(gè)數(shù)；

D_i——第i次測(cè)量的地表平整度。

選取深松深度為280 mm，作業(yè)速度為2.88 km/h的作業(yè)條件進(jìn)行深松平整度和溝槽寬度分析。深松鏟深松后溝槽呈V形，以標(biāo)準(zhǔn)差表示地表平整度，分別取兩側(cè)最高點(diǎn)。試驗(yàn)分別測(cè)量5次求平均值。仿真結(jié)果表明，仿鯊魚(yú)背鰭深松鏟地表平整度提高了20.39%。

3 田間試驗(yàn)

3.1 試驗(yàn)條件

試驗(yàn)設(shè)備包括仿鯊魚(yú)背鰭深松鏟和無(wú)背鰭鑿式深松鏟、深松機(jī)拖拉機(jī)2臺(tái)、數(shù)顯架、S型拉力傳感器（7T）、米尺。試驗(yàn)在錦州義縣試驗(yàn)田進(jìn)行，試驗(yàn)田長(zhǎng)220 m，寬160 m，土壤含水率35%，耕作層、犁底層和心土層的土壤堅(jiān)實(shí)度分別為1 482 kPa、2 173 kPa、2 151 kPa。

3.2 試驗(yàn)方法

在深松機(jī)架上設(shè)置5個(gè)深松鏟，并在每臺(tái)牽引車和每臺(tái)牽引車之間的懸吊裝置間裝有張力傳感器，用來(lái)測(cè)定在工作中的深松阻力，先測(cè)定在無(wú)載荷情況下的耕地阻力，然后測(cè)定在有載荷情況下的深松阻力，二者的差值即為耕作阻力，耕地阻力可由高精度測(cè)力儀直接讀取。

測(cè)試區(qū)根據(jù)需要進(jìn)行了劃分，每個(gè)劃分區(qū)域大小為60 m×30 m。在測(cè)試之前，在深松鏟的位置挖出一個(gè)深的土坑，將深松鏟降到試驗(yàn)設(shè)置的深度，把機(jī)架調(diào)到水平的位置，打開(kāi)數(shù)據(jù)采集設(shè)備，并根據(jù)測(cè)試情況，維持工作速度不變完成測(cè)試工作。

為了與離散元模擬相一致，測(cè)試過(guò)程中，每臺(tái)機(jī)器工作速度為2.88～5.66 km/h。在規(guī)定范圍內(nèi)，每一臺(tái)測(cè)試機(jī)車行駛穩(wěn)定后，在20 m的距離上進(jìn)行測(cè)試，每一組測(cè)試重復(fù)5次，取平均值。

3.3 試驗(yàn)結(jié)果分析

根據(jù)JB/T 10295—2014的有關(guān)規(guī)定，將耕地阻力作為評(píng)價(jià)指標(biāo)，對(duì)帶鰭深松鏟的性能進(jìn)行測(cè)試。機(jī)器空載時(shí)耕作阻力平均值為570 N。結(jié)果表明，當(dāng)深松厚度為300 mm時(shí)，帶鰭深松鏟的真實(shí)耕作阻力（1 904 N）與離散元模擬值（1 941 N）的偏差為1.9%。田間試驗(yàn)結(jié)果表明，離散元仿真得出的分析結(jié)果有效。相對(duì)于錦州義縣常用的直柄鑿式深松鏟，仿鯊魚(yú)背鰭深松鏟作業(yè)時(shí)，耕作阻力較小，土層擾動(dòng)量少，深松面積大，適用于錦州義縣地區(qū)的土壤深松。

4 結(jié)論

1） 針對(duì)錦州義縣農(nóng)田土層厚、質(zhì)地黏稠、含水量高等特點(diǎn)，提出一種新型的仿鯊魚(yú)背鰭深松鏟，該鏟斗可以有效地減小深松工作阻力，有助于打破犁底層，降低土壤堅(jiān)實(shí)度。

2） 采用EDEM2022對(duì)耕層、犁底和心土層三個(gè)層次的土體進(jìn)行模擬，并對(duì)其模擬結(jié)果進(jìn)行分析。研究在相同的深松深度、不同的工作速度下，帶鰭深松鏟的土壤顆?；旌锨闆r、土層擾動(dòng)效應(yīng)以及耕地阻力，獲得3.60～4.33 km/h的最佳工作速度。

3） 將仿真與田間試驗(yàn)進(jìn)行對(duì)比驗(yàn)證分析，當(dāng)深松深度相同時(shí)，帶鰭深松鏟的真實(shí)耕作阻力（1 904 N）與離散元模擬值（1 941 N）的偏差為1.9%。這說(shuō)明該土壤模型與錦州義縣地區(qū)土壤的特點(diǎn)基本一致。

4） 試驗(yàn)表明，提出的帶鰭式深松鏟機(jī)可以有效地破開(kāi)犁底，降低耕地壓力，對(duì)耕地表層結(jié)構(gòu)和表面形貌的影響不大，可為我國(guó)東北地區(qū)耕地科學(xué)研究提供借鑒。

參 考 文 獻(xiàn)

［1］ 胡樹(shù)平， 包海柱， 孟天天， 等. 深松對(duì)土壤物理性質(zhì)及油用向日葵產(chǎn)量性狀的影響［J］.內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2020， 41（3）： 4-9.

Hu Shuping， Bao Haizhu， Meng Tiantian， et al. Effect of sub-soiling on soil characteristics and yield characters of oil sunflower ［J］. Journal of Inner Mongolia Agricultural University （Natural Science Editon）， 2020， 41（3）： 4-9.

［2］ 劉俊安， 王曉燕， 李洪文， 等. 基于土壤擾動(dòng)與牽引阻力的深松鏟結(jié)構(gòu)參數(shù)優(yōu)化［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2017， 48（2）： 60-67.

Liu Junan， Wang Xiaoyan， Li Hongwen， et al. Optimization of structural parameters of subsoiler based on soil disturbance and tractionresistance ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2017， 48（2）： 60-67.

［3］ 張金波， 佟金， 馬云海. 仿生減阻深松鏟設(shè)計(jì)與試驗(yàn)［J］.農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2014， 45（4）： 141-145.

Zhang Jinbo， Tong Jin， Ma Yunhai. Design and experiment of bionic anti-drag subsoiler ［J］. Transactionof the Chinese Society for Agricultural Machinery， 2014， 45（4）： 141-145.

［4］ 周華， 張文良， 楊全軍， 等. 滑切型自激振動(dòng)減阻深松裝置設(shè)計(jì)與試驗(yàn)［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2019， 50（5）： 71-78.

Zhou Hua， Zhang Wenliang， Yang Quanjun， et al. Design and experiment of sliding cutting self-excited vibration dragreduction subsoiling device ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2019， 50（5）： 71-78.

［5］ 張強(qiáng)， 張璐， 于海業(yè)， 等. 復(fù)合形態(tài)深松鏟耕作阻力有限元分析與試驗(yàn)［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2012， 43（8）： 61-65.

Zhang Qiang， Zhang Lu， Yu Haiye， et al. Finite element analysis and experiment of soil resistance of multiplex-modality subsoiler ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2012， 43（8）： 61-65.

［6］ 郭志軍， 杜干， 李忠利， 等. 觸土曲面形式對(duì)推土板減阻性能影響的正交試驗(yàn)［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2015， 46（7）： 372-378.

Guo Zhijun， Du Gan， Li Zhongli， et al. Orthogonal experiment on resistance reduction by soil-engaging surfaces of bulldozer blade ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2015， 46（7）： 372-378.

［7］ Sun J， Wang Y， Ma Y， et al. DEM simulation of bionic subsoilers （tillage depthgt; 40 cm） with drag reduction and lower soil disturbance characteristics ［J］. Advances in Engineering Software， 2018， 119： 30-37.

［8］ 王金武， 李響， 高鵬翔， 等. 胡蘿卜聯(lián)合收獲機(jī)高效減阻松土鏟設(shè)計(jì)與試驗(yàn)［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2020， 51（6）： 93-103.

Wang Jinwu， Li Xiang， Gao Pengxiang， et al. Design and experiment of high efficiency drag reducing shovel for carrot combine harvester ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2020， 51（6）： 93-103.

［9］ 劉爽， 呂超， 饒勇， 等. 基于仿真的鯊魚(yú)鰭結(jié)構(gòu)性能分析與研究［J］. 系統(tǒng)仿真學(xué)報(bào)， 2018， 30（6）： 2398-2404.

Liu Shuang， Lü Chao， Rao Yong， et al. Performance analysis and research of shark fin structure based on simulation ［J］. Journal of System Simulation， 2018， 30（6）： 2398-2404.

［10］ 張智泓， 甘帥匯， 左國(guó)標(biāo)， 等. 以砂魚(yú)蜥頭部為原型的仿生深松鏟尖設(shè)計(jì)與離散元仿真［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2021，52（9）： 33-42.

Zhang Zhihong， Gan Shuaihui， Zuo Guobiao， et al. Bionic design and performance experiment of sandfish head inspired subsoiler tine ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2021， 52（9）： 33-42.

［11］ Zhang R， Li J. Simulation on mechanical behavior of cohesive soil by Distinct Element Method ［J］. Journal of Terramechanics， 2005， 43（3）： 303-316.

［12］ 于建群， 錢立彬， 于文靜， 等.開(kāi)溝器工作阻力的離散元法仿真分析［J］.農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2009， 40（6）： 53-57.

Yu Jianqun， Qian Libin， Yu Wenjing， et al. Dem analysis of the resistances applied on furrow openers ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2009， 40（6）： 53-57.

［13］ 李吉成， 孫坤鵬， 張永華， 等. 基于離散元法的雙翼深松鏟耕作行為的仿真分析［J］. 內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2023， 44（2）： 67-75.

Li Jicheng， Sun Kunpeng， Zhang Yonghua， et al. Simulation of tillage behavior of double winged subsoiler based on discrete element method ［J］. Journal of Inner Mongolia Agricultural University （Natural Science Editon）， 2023， 44（2）： 67-75.

［14］ 張喜瑞， 曾望強(qiáng)， 劉俊孝， 等. 基于離散元法的磚紅壤斜柄折翼式深松鏟設(shè)計(jì)與試驗(yàn)［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2022， 53（3）： 40-49.

Zhang Xirui， Zeng Wangqiang， Liu Junxiao， et al. Design and experiment of iateritic soil inclined handle folding wing subsoiling shovel based on discrete element method ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2022， 53（3）： 40-49.

［15］ 高連龍. 深松過(guò)程中玉米秸稈運(yùn)動(dòng)研究及深松鏟的設(shè)計(jì)［D］. 哈爾濱：東北農(nóng)業(yè)大學(xué)， 2022.

Gao Lianlong. Study on movement of corn straw in subsoiling process and design of subsoiler ［D］. Harbin： Northeast Agricultural University， 2022.

［16］ 丁啟朔， 李楊， BELAL Eisa Adam， 等. 基于田間攝像的多參數(shù)水稻土深松擾動(dòng)行為與效應(yīng)研究［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2019， 50（10）： 44-55.

Ding Qishuo， Li Yang， BELAL Eisa Adam， et al. Subsoiler-induced paddy soil disturbance and effects based on video-assisted multi-index quantification ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2019， 50（10）： 44-55.