趙淑紅 王加一 陳君執(zhí) 楊悅乾 譚賀文

(東北農(nóng)業(yè)大學工程學院， 哈爾濱 150030)

0 引言

為了打破堅硬的土壤犁底層，同時降低耕作阻力和表層土壤擾動量，選擇在免耕播種地區(qū)進行深松。目前深松鏟類型多為折線式、圓弧式，折線式深松鏟更多用于深松滅茬機上，耕作阻力較大，對拖拉機的動力要求較高，不適合機械大范圍作業(yè)；圓弧式深松鏟多用于深松機上，其土壤表層擾動量較大，土壤水肥損失嚴重，不能滿足保護性耕作的要求[1-2]。

近些年來，學者對深松鏟的研究取得了大量成果[3-9]，其主要從振動和仿生原理兩方面對深松鏟進行設計與研究。

隨著軟件的發(fā)展，仿真分析以減少樣機的加工及試驗時間的優(yōu)點，逐漸成為模擬試驗的一種方法。離散元法可用來模擬顆粒材料和研究材料間的微觀及宏觀變化。近幾年學者通過離散元軟件仿真模擬觸土部件工作狀態(tài)，并與試驗進行對比，驗證了離散元軟件模擬土壤與試驗誤差較小[10-11]，因此離散元軟件可用于模擬觸土部件與土壤模型。在觸土部件與土壤研究中，觸土部件作用于土壤后，土壤往往在運動行為上反作用于觸土部件。對于深松鏟的研究，鏟柄使土壤沿鏟柄側(cè)向遠離鏟柄，土壤在運動行為上對于鏟柄的影響較小，而鏟尖作用于土壤后，土壤在運動行為上沿鏟柄運動方向較明顯，對鏟柄的影響較大，且土壤軌跡沿鏟柄運動可以增加滑移，減少摩擦[12-15]。因此本文依據(jù)鏟尖對土壤的作用設計鏟柄，建立鏟尖與土壤的仿真模型，獲得鏟尖上方土壤顆粒運動軌跡的擬合曲線，設計擬合曲線型弧狀深松鏟，運用離散元法對鏟柄、鏟尖的互作效應進行仿真試驗，確定設計的合理性，并進行土槽對比試驗與田間對比試驗，驗證設計的優(yōu)良性。

1 離散元仿真分析

深松鏟的主要作用原理是深松鏟鏟尖在前端的深入開土以及后端深松鏟鏟柄的破土。為設計鏟柄形狀，選取鏟尖對土壤的影響進行分析?？紤]到實際中鏟尖不可能單獨運動，設置連接段為任何位置，都會對試驗結(jié)果造成影響，因此選取離散元仿真軟件對其分析。

1.1 仿真模型與鏟尖模型的建立

本試驗選取常用的箭形深松鏟鏟尖進行試驗，為保證鏟尖作業(yè)范圍，設置長1 000 mm、寬400 mm、高500 mm的土壤仿真模型。

依據(jù)文獻[16]規(guī)定的深松鏟鏟尖基本尺寸，建立鏟尖模型，如圖1。保存為igs格式導入EDEM軟件中。材料為Q235鋼，密度為7 850 kg/m3，剪切模量為7.7×1010Pa，泊松比為0.3。

圖1 深松鏟鏟尖結(jié)構(gòu)圖Fig.1 Structure diagrams of subsoiler tip

1.2 土壤模型的建立

1.2.1土壤接觸模型

根據(jù)牛頓第二定律，每個顆粒的運動方程為

(1)

(2)

式中vi——顆粒i的移動速度

Ii——顆粒i的轉(zhuǎn)動慣量

ni——顆粒個數(shù)

ωi——顆粒i的角速度

F——顆粒i的粘聚力

試驗地區(qū)為東北壤土，具有粒度均勻、透氣性好、透水性好、強烈脹縮和擾動特性的特點，土壤有一定粘附力和彈塑性[19]。因此土壤顆粒與深松鏟鏟尖之間設定為Hertz-Mindlin(no slip)接觸模型，土壤顆粒與土壤顆粒之間除設置Hertz-Mindlin(no slip)接觸模型外，考慮到壤土之間內(nèi)聚特性較強，土壤間相互粘結(jié)，因此添加Hertz-Mindlin with bonding接觸模型。依據(jù)文獻[19]對于試驗區(qū)相同類型土壤顆粒微觀參數(shù)的測定和文獻[20]適應大多數(shù)土壤模型所用的微觀參數(shù)，如表1所示。

1.2.2土壤顆粒模型

為研究每一土壤顆粒的運動軌跡，設置土壤顆粒模型為單一顆粒。并修正實際中的土壤顆粒粒徑，設置土壤顆粒半徑為4 mm。

表1 離散元法仿真的微觀參數(shù)Tab.1 Microscopic parameters of DEM simulation

1.3 仿真分析

1.3.1仿真幾何模型的設置

在所設定的仿真模型中生成土壤顆粒，固定時間步長為8.2×10-5s，Rayleigth時間步長為5.47×10-4s?？偼寥李w粒生成時間為8 s，土壤顆粒自然沉降1 s，鏟尖入土深度為300 mm，運動時間為1 s，運動速度選取與深松鏟正常作業(yè)相近的0.5 m/s，沿X軸正向，如圖2所示。

圖2 仿真幾何模型Fig.2 Simulation geometry model

圖5 區(qū)域土壤顆粒的運動軌跡圖Fig.5 Soil particle trajectory maps of various areas

1.3.2結(jié)果與分析

仿真完成后，觀察鏟尖作業(yè)后土壤的運動效果。由于鏟尖作用正上方土壤顆粒在y軸方向上位移不明顯，并參考文獻[13]，截取鏟尖作用的鉛垂面為剖視圖，選擇鏟尖作業(yè)正上方的土壤顆粒，如圖3，并分析其運動軌跡，如圖4所示。

圖3 土壤顆粒的選擇Fig.3 Selection of soil particle

圖4 土壤顆粒軌跡圖Fig.4 Trajectory map of soil particle

由圖4可以看出，每個土壤顆粒在X-Z平面內(nèi)運動位移較大，且運動軌跡可分為3個過程，首先土壤顆粒被鏟尖擾動上升，其次在鏟尖通過后對前方的土粒再次擾動升起的過程中使其后旋，最后在鏟尖完全通過此區(qū)域后回落。根據(jù)顆粒在不同深度的運動規(guī)律的不同，可將其沿Z軸負向平均分為5個區(qū)域(即深度方向)，深度分別為60 mm，分析其運動軌跡，局部圖如圖5所示。

圖6 區(qū)域單一土壤顆粒的軌跡圖Fig.6 Single soil particle trajectory maps of various areas

圖7 土壤顆粒的軌跡擬合曲線Fig.7 Trajectory fitting curves of soil particles

為設計鏟柄，選取回落軌跡段進行分析。對每一區(qū)域土壤顆粒運動軌跡進行分析。觀察土壤顆粒坐標變化，其在Z軸方向最小位移為6 mm，最大位移為25 mm，且土壤顆粒直徑為8 mm。為選取時在Z軸方向土壤軌跡不遺漏，在每一區(qū)域內(nèi)按深度方向每間隔6 mm選取一粒土壤顆粒(選取5粒)，輸出每一顆粒的數(shù)據(jù)結(jié)果分別為每一時刻對應的X、Z坐標，時間步長為0.01 s，其Z坐標值顯著變小時為土壤顆粒的回落軌跡段，得到在同一區(qū)域內(nèi)，5顆土壤顆粒的回落軌跡段對于同一線形擬合方程的擬合度都較高。因此選擇單一顆粒運動軌跡代表單一區(qū)域土壤顆粒運動軌跡，軌跡圖如圖6所示。選取其中單一顆粒，回落軌跡段內(nèi)的坐標值在Matlab中做顆粒軌跡擬合曲線，結(jié)果如圖7所示，并得到第1區(qū)域土壤軌跡擬合曲線回歸方程為

(3)

第2區(qū)域土壤軌跡擬合曲線回歸方程為

y=0.112 2x5+238.9x4+2.034×105x3+

8.658×1010x2+1.843×1010x+1.569×1012

(4)

第3區(qū)域土壤軌跡擬合曲線回歸方程為

y=1.382×10-7x9-3.012×10-5x8+
0.002 905x7-0.162 8x6+5.843x5-
139.2x4+2 202x3-2.229×104x2+
1.311×105x-3.415 2×105

(5)

第4區(qū)域土壤軌跡擬合曲線回歸方程為

y=-9.68×10-5x4+9.908×10-4x3+
0.444 1x2-0.111 3x-433.9

(6)

第5區(qū)域土壤軌跡擬合曲線回歸方程為

y=1.702×10-7x7+1.718×10-5x6+
7.024×10-5x5+0.014 74x4+0.165 3x3+
0.944x2+2.342x-402.6

(7)

各個回歸系數(shù)分別為0.975、0.997、0.997、0.990、0.925。每一區(qū)域數(shù)學函數(shù)曲線的擬合度都較高，可靠性較強，因此可以基于此線形設計深松鏟鏟柄的外輪廓曲線。

2 深松鏟設計

傳統(tǒng)深松鏟一般包括鏟柄與機架連接段、鏟柄和鏟尖，如圖8所示。鏟柄與機架連接段主要與機架連接，要求連接穩(wěn)定性較強，并且可以調(diào)節(jié)深松鏟入土深淺；主要作業(yè)部分為鏟柄和鏟尖。

圖8 傳統(tǒng)深松鏟示意圖Fig.8 Diagram of traditional subsoiler1.鏟柄與機架連接段 2.鏟柄 3.鏟尖

2.1 鏟柄

為與大部分深松機滅茬機械、深松機、中耕施肥機械相配合，要求鏟柄與機架連接段的設計與標準深松鏟尺寸相符合[21]。

2.1.1鏟柄線形的理論設計

為降低深松鏟與土壤的碰撞，減小其擾動面積，同時降低耕作阻力。根據(jù)鏟尖作業(yè)后土壤顆粒的運動軌跡，由表2建立數(shù)學模型平面圖。經(jīng)測量，東北壟作玉米播種深度為(50±10) mm，常年耕作形成的耕作層厚度為120 mm左右，犁底層厚度約為120 mm，心土層厚度為240 mm左右，在AutoCAD中建立各段曲線與土壤層對應關(guān)系，如圖9所示。

2.1.2鏟柄線形的優(yōu)化

由圖9可看出，鏟尖作用后土壤顆?；芈滠壽E1、2段曲線與水平方向夾角為銳角，可以利用滑移原理降低碰撞量和耕作阻力[22-24]。故將1、2段曲線的整體點坐標導入Matlab中擬合為第1段，作業(yè)時對應耕作層土壤；3、4段曲線呈現(xiàn)與水平方向的夾角由銳角向鈍角過渡，因此將3、4段曲線的點坐標導入Matlab中擬合為第2段，作業(yè)時對應犁底層土壤；并將5段曲線的點坐標導入Matlab中擬合為第3段，5段曲線與水平方向夾角為鈍角，可以增加對上方土壤的翻垡[1]，故對應部分在犁底層之下，為心土層土壤，對犁底層的打破起到主要作用。本次擬合各段曲線在各方向上有平移，但曲線線形不變，各段線形擬合曲線如圖10，擬合方程如表2所示。

表2 優(yōu)化后土壤軌跡曲線方程Tab.2 Optimized soil trajectory curvilinear equation

圖9 數(shù)學函數(shù)曲線圖Fig.9 Graph of mathematical function

圖10 優(yōu)化后土壤軌跡曲線Fig.10 Optimized soil trajectory fitting curves

保證各段擬合曲線線形不變的前提條件下在AutoCAD中構(gòu)建深松鏟基本模型，按各層土壤厚度比相應擴大3段曲線。依據(jù)文獻[25-26]，采用線元設計法，將線形與連接不合理處優(yōu)化為由多段圓弧和緩和曲線段(或直線段)順序組合而成，從而減少連接處的應力集中問題，使鏟柄受力更均勻。

在鏟柄后端曲線的設計上，為適應土壤回落軌跡，增加回土量，降低土壤擾動，同時保證深松鏟作業(yè)時受力更均勻，避免應力集中現(xiàn)象，鏟柄后端曲線設計為優(yōu)化后的3段擬合曲線，結(jié)果如圖11。

圖11 擬合曲線型深松鏟示意圖Fig.11 Diagrams of fitting curve subsoiler

2.1.3刃口的設計

現(xiàn)深松鏟刃部多為楔形刃或無刃口，楔形刃錐角設計較小時能夠有效降低牽引阻力，但穩(wěn)定性較差，容易折刃[27]，不開刃深松鏟雖然穩(wěn)定性較強，但其牽引阻力較大[28]。部分增加刃口的深松鏟刃口角度一般取60°，為楔形刃口，由文獻[2]及觀察深松鏟作業(yè)時的實際狀況，東北地區(qū)深松作業(yè)時，會使鏟柄前段堆積為土核，如圖12所示。

圖12 土壤堆積示意圖Fig.12 Diagrams of soil accumulation

現(xiàn)有的楔形刃深松鏟和不開刃深松鏟，在作業(yè)過程中刃口端堆聚的土核較多，增大了作業(yè)阻力，從而使作業(yè)性能受到影響。而國內(nèi)文獻中，針對此討論的文獻較少，單純對深松鏟模型進行受力分析，與實際工況下深松鏟的形態(tài)有一定的差異性。且土核的堆聚多為外凸曲線狀。因此依據(jù)文獻[2]，在刃口的線形設計上采取介于楔形刃和無刃口之間的外凸形輪廓曲線，降低耕作阻力的同時保證刃口的穩(wěn)定性，且在作業(yè)中堆積的土核與輪廓曲線形狀基本相似，可以減少土壤粘附，土核的堆聚可以忽略不計。得到深松鏟刃口結(jié)構(gòu)及受力如圖13、14所示。

圖13 刃口曲線示意圖Fig.13 Diagram of blade curve

圖14 刃口受力示意圖Fig.14 Force diagram of edge

根據(jù)文獻[29]，刃口受到的阻力為

(8)

對式(8)進行推導得到

(9)

其中，將ABC視為直角三角形，A處為直角，則

(10)

式中β——刃口切削刃角在水平方向的投影，(°)

F′N——土壤對深松鏟的法向作用力，N

F′T——土壤對深松鏟刃的作用力，N

μ——深松鏟與土壤間的動摩擦因數(shù)

θ——刃口切削刃角，(°)

α——刃口切削刃角平面與其在水平面間的夾角，(°)

在深松鏟作業(yè)速度一定時，F(xiàn)′N、F′T視為理想常數(shù)值，μ視為定值，則刃口所受到的阻力與刃口切削刃角在水平方向的投影有關(guān)，β越小，受到的切削阻力越小。β又與α、θ有關(guān)，α越大，θ越小，β越小。深松鏟柄切削刃角一般取值為30°～60°，其中壤土的最小切削刃角范圍為40°～45°[30]，外凸形曲線刃口穩(wěn)定性較強，為減小阻力，取刃口切削角為40°。

2.2 深松鏟作業(yè)深度、入土角、刃口高度的確定

犁底層是多年耕作形成的，打破犁底層可以使深松作業(yè)不用每年都進行，因此打破犁底層是設計深松鏟作業(yè)深度時首先要考慮的因素，且犁底層的深度為120～240 mm，本設計中取深松鏟作業(yè)深度為300 mm。深松鏟的入土角是影響深松阻力的重要因素之一。入土角在21°[31]時耕作阻力最小，因此本設計入土角選為21°。

東北地區(qū)的壤土具有強烈脹縮和擾動特性，彈塑性極限較高，所以要對其破壞，刃口高度設置不宜過低，但刃口高度設置過高，穩(wěn)定性會有所降低，容易折刃[27]，因此取刃口高度為15 mm，保證切削刃角的同時，外凸形樣條曲線連接。材料選為65Mn鋼制造，刃部應進行熱處理，淬火區(qū)為整個刃口部位，硬度為HRC48～56。深松鏟如圖15所示。

圖15 擬合曲線型深松鏟Fig.15 Fitting curve subsoiler

2.3 仿真驗證試驗

為證明設計思路的合理性，同時驗證此設計的可行性，選擇在EDEM軟件中對所設計深松鏟進行仿真試驗，進一步驗證鏟尖作用時土壤顆粒的運動與所設計鏟柄的關(guān)系。

2.3.1模型的建立

選擇上述鏟尖仿真試驗的幾何模型和土壤模型，按所設計深松鏟尺寸在CATIA中建立擬合曲線型深松鏟三維圖，如圖16所示，與鏟尖通過螺栓連接裝配，導入EDEM軟件中。材料為65Mn鋼，密度為7 865 kg/m3，剪切模量為7.9×1010Pa，泊松比為0.3。

圖16 擬合曲線型深松鏟三維圖Fig.16 Three-dimensional diagram of fitting curve subsoiler

2.3.2仿真幾何模型的設置

建立深松鏟與土壤的裝配模型。固定時間步長為8.2×10-5s，Rayleigth時間步長為5.47×10-4s。深松鏟入土深度為300 mm，運動時間為1 s，運動速度為0.5 m/s，沿X軸正向，如圖17所示。

圖17 深松鏟仿真幾何模型Fig.17 Simulation geometry model of subsoiler

2.3.3結(jié)果與分析

仿真完成后，觀察深松鏟作業(yè)時土壤的運動軌跡。截取鏟尖作用的鉛垂面為剖視圖，選擇鏟尖作業(yè)正上方的土壤顆粒，并分析其運動軌跡，如圖18所示。

圖18 深松鏟作業(yè)時土壤顆粒軌跡圖Fig.18 Trajectory map of soil particle by subsoiler

由仿真結(jié)果可得鏟尖與鏟柄在共同作業(yè)時，前方土壤顆粒在鏟尖和鏟柄前端及刃口第3段擬合曲線作用下向上翻垡，與第1段擬合曲線產(chǎn)生滑移。在深松鏟作業(yè)后，土壤整體回落趨勢軌跡與鏟柄后端曲線基本一致。

因此仿真試驗驗證了鏟柄前端及刃口設計曲線對土壤的作用形式；驗證了鏟柄后端曲線與土壤整體回落曲線一致；同時驗證了本設計的可行性，以及本文總體設計思路的合理性。

3 土槽對比試驗

3.1 試驗目的

對上述所設計深松鏟，取耕作阻力與土壤擾動量為試驗指標，并探究深松鏟對犁底層深度土壤容重的影響，對擬合曲線型深松鏟與折線式、圓弧式深松鏟進行對比試驗，總結(jié)分析其試驗效果，觀察優(yōu)劣性。

3.2 試驗條件

力的對比試驗對試驗環(huán)境的穩(wěn)定性要求較高，同時為了在對比試驗效果時控制其它無關(guān)因素的一致性，因此對比試驗采取土槽試驗。在黑龍江省農(nóng)業(yè)機械工程科學研究院的室內(nèi)土槽中進行試驗，試驗設備如圖19所示，試驗選取土槽長50 m，寬3 m，土壤厚度1 m，其土壤為典型東北壤土，試驗區(qū)所測土壤深度范圍0～300 mm。人工澆水改變土壤含水率，并利用干燥法測量土壤含水率和容重，利用硬度儀測量土壤硬度。所有數(shù)據(jù)測量3點，取平均值。土壤參數(shù)如表3所示。

圖19 土槽試驗Fig.19 Soil bin test1.土槽車 2.測力架 3.限深輪 4.土槽 5.深松鏟

參數(shù)深度/mm數(shù)值含水率/%0～30014±120±126±1土壤容重/(g·cm-3)0～3001.421.381.300～125515421366硬度/kPa125～250200115671123250～3001018725676

試驗所用器材有全液壓四輪驅(qū)動車(土槽車)、黑龍江省農(nóng)業(yè)機械工程科學研究院智能化研究所研制的六分力三向側(cè)力裝置。側(cè)力裝置分為上懸掛架和下懸掛架，兩個架子間以6個拉力傳感器(3個測量垂直力、2個測量水平力、1個測量側(cè)向力)鉸接相連。擬合曲線型深松鏟與折線式、圓弧式深松鏟分別與懸掛架通過頂絲固定。

3.3 試驗方案

分別取擬合曲線型深松鏟與折線式深松鏟、圓弧式深松鏟單獨在土槽車工作速度為1、2、3 km/h的牽引帶動下對土槽進行深松作業(yè)，深松鏟入土深度都為300 mm。通過拉力傳感器分別測出每種深松鏟作業(yè)時的耕作阻力，利用溝形輪廓儀[32]對3種深松鏟的溝形進行測量，并描繪在坐標紙上，示意圖如圖20，B為擾動寬度，D為堆頂寬度，d為開溝寬度，H為開溝深度，h為回土深度。同時隨機取犁底層深度土壤樣本，計算其土壤容重，為防止取樣偶然性誤差，每種鏟型重復3次試驗，計算其平均值，比較結(jié)果。

圖20 溝形示意圖Fig.20 Diagram of soil groove

圖21 擾動面積與耕作阻力的性能對比Fig.21 Comparisons of performance between disturbed area and tillage resistance

3.4 試驗結(jié)果與分析

犁底層深度土壤容重如表4所示。

表4 土壤容重統(tǒng)計Tab.4 Soil bulk density statistics g/cm3

記錄并計算出試驗所得數(shù)據(jù)，建立3種深松鏟的擾動面積與速度、耕作阻力與速度在不同含水率下的效果對比圖。

由圖21得，擾動面積隨著速度的增加而增大。同一工況下圓弧式深松鏟的土壤擾動面積大于折線式深松鏟，折線式深松鏟土壤擾動面積大于擬合曲線型深松鏟；在牽引阻力方面，折線式深松鏟大于圓弧式深松鏟，圓弧式深松鏟大于擬合曲線型深松鏟。圓弧式、折線式、擬合曲線型深松鏟平均擾土面積為8 749.6、6 377.4、2 959.13 mm2，平均阻力為1 976、1 779.3、1 260.1 N。擬合曲線型深松鏟土壤擾動量比折線式和圓弧式減少53.6%、66.18%，耕作阻力減少36.23%、29.18%。因此在深松鏟的作用效果方面，擬合曲線型深松鏟比其它深松鏟性能較優(yōu)。并由表4得，深松鏟的作用能顯著減小土壤容重，達到對犁底層的深松效果。

4 田間對比試驗

4.1 試驗目的與條件

本次試驗田間土壤平均含水率為(20±1)%，未涉及到粘土量，且土壤板結(jié)較嚴重。因此選擇土壤擾動量為指標，對擬合曲線型、圓弧式、折線式深松鏟進行對比試驗，并探究深松部件對犁底層土壤容重的影響，2017年5月在黑龍江省東北農(nóng)業(yè)大學試驗田進行田間試驗，如圖22所示。

圖22 田間試驗Fig.22 Field tests1.深松施肥機 2.擬合曲線型深松鏟 3.圓弧式深松鏟 4.約翰迪爾484型拖拉機 5.折線式深松鏟 6.限深輪

試驗區(qū)長50 m、寬6 m、土壤厚度1 m，其土壤為典型東北壤土。平均土壤容重為1.50 g/cm3，土壤硬度分別為：1 526 kPa(0～125 mm)、3 649 kPa(125～250 mm)、2 433 kPa(250～300 mm)，調(diào)節(jié)限深輪與鏟尖的高度差保證所設計3段工作刃與試驗區(qū)相應土層相對應。試驗所用器材有約翰迪爾484型拖拉機、深松施肥機(僅深松部件工作)、3種深松鏟。

4.2 試驗方案

拖拉機以正常田間作業(yè)2 km/h的速度進行深松作業(yè)，深松深度為300 mm。借鑒文獻[33-34]對觸土部件工作性能的評價方法，選取回土面積作為深松鏟擾動量的評價指標。利用溝形輪廓儀[32]先對3種深松鏟自然回土的橫截面積進行測量，在坐標紙上繪制溝形曲線；而后將溝內(nèi)回落土壤扒出，使其深度達到開溝深度，利用同一張坐標紙再次繪制溝形曲線。其回土面積示意圖如圖23所示，每種深松鏟所開溝形隨機測量5次，測量出回土面積，求其平均值。并隨機選取各深松鏟作業(yè)后的犁底層土壤樣本5次，計算其土壤容重。

圖23 田間試驗溝形示意圖Fig.23 Diagram of soil groove in field tests

4.3 試驗結(jié)果

作業(yè)效果如圖24所示。

圖24 作業(yè)效果Fig.24 Operation effect

依據(jù)文獻[33-36]，土壤的擾動量可以由其回土量的多少來反映，相同深松深度條件下，回土量越大，其土壤擾動量越小，深松區(qū)域土壤越松散，效果越好；反之，回土量越小，其擾動量越大，深松區(qū)域土壤越少，水肥散失越嚴重。土壤容重則對作物苗期生長有顯著性影響[37-38]，3種鏟型的回土量及犁底層土壤容重如表5所示。

表5 驗證試驗結(jié)果Tab.5 Verification results of test

由表5各種深松鏟的指標平均值可知，擬合曲線型深松鏟的回土面積大于折線式，折線式深松鏟的回土面積大于圓弧式；擬合曲線型深松鏟的土壤擾動量小于折線式，折線式深松鏟的土壤擾動量小于圓弧式，因此擬合曲線型深松鏟深松效果最優(yōu)。3種深松鏟作業(yè)能夠顯著減小犁底層土壤容重，達到打破犁底層的要求。

5 結(jié)論

(1)依據(jù)深松鏟鏟尖運動的離散元(DEM)仿真試驗分析，設計了擬合曲線型深松鏟，并通過鏟柄與鏟尖互作效應仿真試驗驗證了設計思路的可行性與合理性。所設計擬合曲線型深松鏟可以有效減小耕作阻力和表層土壤的擾動量，同時打破犁底層，降低犁底層的土壤容重。

(2)對所設計的擬合曲線型深松鏟與其他類型深松鏟進行土槽對比試驗，測出擬合曲線型深松鏟比折線式深松鏟土壤擾動量減少了53.6%，耕作阻力減少了36.23%；擬合曲線型深松鏟比圓弧式深松鏟土壤擾動量減少了66.18%，耕作阻力減少了29.18%。驗證擬合曲線型深松鏟牽引阻力最小，土壤擾動量最小，指標較優(yōu)。并且減小犁底層土壤容重，達到打破犁底層的要求。

(3)對擬合曲線型深松鏟與其他類型深松鏟進行田間對比試驗，進一步證明擬合曲線型深松鏟能夠減小犁底層土壤容重，達到打破犁底層的要求。測出擬合曲線型、折線式、圓弧式深松鏟的回土面積分別為7 745.4、4 278.6、3 136.4 mm2，分別增加81.03%和146.95%，驗證了擬合曲線型深松鏟比其他類型深松鏟在回土面積指標上較優(yōu)。田間試驗與土槽試驗均說明，擬合曲線型深松鏟深松作業(yè)時土壤擾動量較小，更符合保護性耕作的要求。

1 王燕. 基于離散元法的深松鏟結(jié)構(gòu)與深松效果研究[D]. 長春：吉林農(nóng)業(yè)大學，2014.

WANG Yan. Simulation analysis of structure and effect of the subsoiler based on DEM[D]. Changchun: Jinlin Agricultural University, 2014.(in Chinese)

2 龔皓暉. 基于土壤力學模型的深松鏟有限元分析與結(jié)構(gòu)優(yōu)化[D]. 成都：西華大學，2013.

GONG Haohui.Finite element analysis and structure optimization for subsoiling shovel based on soil mechanics model[D]. Chengdu:Xihua University, 2013.(in Chinese)

3 邱立春，李寶筏. 自激振動深松機減阻試驗研究[J].農(nóng)業(yè)工程學報，2000，16(6)：72-76.

QIU Lichun, LI Baofa. Experimental study on the self-excited vibration subsoiler for reducing draft force[J].Transactions of the CSAE, 2000, 16(6):72-76.(in Chinese)

4 SHAHGOLIA G, FIELKE J, SAUNDERS C, et al. Simulation of the dynamic behaviour of a tractor-oscillating subsoiler system[J]. Biosystems Engineering, 2010, 106(2):147-155.

5 李霞，付俊峰，張東興，等. 基于振動減阻原理的深松機牽引阻力試驗[J]. 農(nóng)業(yè)工程學報，2012，28(1)：32-36.

LI Xia, FU Junfeng, ZHANG Dongxing, et al. Experiment analysis on traction resistance of vibration subsoiler[J].Transactions of the CSAE,2012,28(1):32-36. (in Chinese)

6 張強，張璐，劉憲軍，等. 基于有限元法的仿生鉤形深松鏟耕作阻力[J]. 吉林大學學報：工學版，2012，42(增刊1)：117-121.

ZHANG Qiang, ZHANG Lu, LIU Xianjun, et al. Soil resistance of the bionic hook-shape subsoiler based on the finite element method[J].Journal of Jilin University:Engineering and Technology Edition, 2012, 42(Supp.1):117-121.(in Chinese)

7 龔皓暉，王霜. 仿生深松鏟結(jié)構(gòu)設計與有限元分析[J].農(nóng)機化研究，2013，35(7)：53-57.

GONG Haohui, WANG Shuang. Bionic subsoiler sructural design with the finite element analysis[J].Journal of Agricultural Mechanization Research, 2013, 35(7):53-57.(in Chinese)

8 張金波，佟金，馬云海. 仿生減阻深松鏟設計與試驗[J/OL]. 農(nóng)業(yè)機械學報，2014，45(4)：141-145.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20140422&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2014.04.022.

ZHANG Jinbo,TONG Jin, MA Yunhai. Design and experiment of bionic anti-drag subsoiler[J/OL].Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(4):141-145.(in Chinese)

9 朱風武，佟金. 土壤深松技術(shù)及高效節(jié)能仿生研究的發(fā)展[J].吉林大學學報：工學版，2003，33(2)：95-99.

ZHU Fengwu,TONG Jin. Development of high efficient and energy saving bionic bubsoiling techniques[J].Journal of Jilin University:Engineering and Technology Edition, 2003, 33(2):95-99.(in Chinese)

10 黃玉祥，杭程光，苑夢嬋，等. 深松土壤擾動行為的離散元仿真與試驗[J/OL]. 農(nóng)業(yè)機械學報，2016，47(7)：80-88.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20160712&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2016.07.012.

HUANG Yuxiang,HANG Chengguang, YUAN Mengchan, et al. Discrete element simulation and experiment on disturbance behavior of subsoiling[J/OL].Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(7):80-88.(in Chinese)

11 于建群，錢立彬，于文靜，等. 開溝器工作阻力的離散元法仿真分析[J].農(nóng)業(yè)機械學報，2009，40(6)：53-57.

YU Jianqun, QIAN Libin, YU Wenjing, et al. Discrete element method of the resistances applied on furrow openers[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(6):53-57.(in Chinese)

12 丁啟朔. 耕作力學研究的土壤結(jié)構(gòu)及其評價方法[D]. 南京：南京農(nóng)業(yè)大學，2006.

DING Qishuo. Soil structure and its assessment for soil tillage research[D]. Nanjing:Nanjing Agricultural University, 2006.(in Chinese)

13 丁啟朔，任駿，BELAL Eisa Adam，等. 濕粘水稻土深松過程離散元分析[J/OL].農(nóng)業(yè)機械學報，2017，48(3)：38-48.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20170305&journal_id=jcsam.DOI:10.6041/jissn. 1000-1298.2017.03.005.

DING Qishuo, REN Jun,BELAL Eisa Adam, et al. DEM analysis of subsoiling process in wet clayey paddy soil[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(3):38-48.(in Chinese)

14 孟鳳英. 粘性土壤破碎機理及分形規(guī)律研究[D]. 南京：南京農(nóng)業(yè)大學，2009.

MENG Fengying. A study on the fragmentation mechanism of cohesive soil and its fractal characteristics [D]. Nanjing:Nanjing Agricultural University, 2009.(in Chinese)

15 周祖鍔. 農(nóng)業(yè)物料學[M].北京：中國農(nóng)業(yè)出版社，1994．

16 JB/T 9788—1999 深松鏟和鏟柄[S].1999．

17 胡建平，周春健，侯沖，等. 磁吸板式排種器充種性能離散元仿真[J/OL].農(nóng)業(yè)機械學報，2014，45(2)：94-98.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20140216&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2014.02.016.

HU Jianping, ZHOU Chunjian, HOU Chong, et al. Simulation analysis of seed-filling performance of magnetic plate seed-metering device by discrete element method[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(2):94-98.(in Chinese)

18 韓燕龍，賈富國，唐玉榮，等. 顆粒滾動摩擦系數(shù)對堆積特性的影響[J]. 物理學報，2014，63(17)：165-171.

HAN Yanlong,JIA Fuguo, TANG Yurong, et al. Influence of granular coefficient of rolling friction on accumulation characteristics[J].Acta Physica Sinica, 2014, 63(17):165-171.(in Chinese)

19 曹秀振. 基于離散元法的免耕播種開溝器的設計與試驗[D]. 哈爾濱：東北農(nóng)業(yè)大學，2016.

CAO Xiuzhen. Design and test of no-till seeding opener on discrete element method[D]. Harbin:Northeast Agricultural University, 2016.(in Chinese)

20 MUSTAFA U, JOHN M F, CHRIS S.Three dimensional discrete element modeling DEM of tillage accounting for soil cohesion and adhesion[J]. Biosystems Engineering, 2015, 129:298-306.

21 中國農(nóng)業(yè)機械化科學研究院. 農(nóng)業(yè)機械設計手冊[M].北京：中國農(nóng)業(yè)科學技術(shù)出版社，2007．

22 吉爾 W R,范德伯奇 G E. 耕作和牽引土壤動力學[M]. 耕作和牽引土壤動力學翻譯組,譯. 北京:中國農(nóng)業(yè)機械出版社,1983.

23 DAMORA D, PANDEY K P. Evaluation of performance of furrow openers of combined seed and fertilizer drills[J]. Soil and Tillage Research, 1995, 34(2): 127-139.

24 TESSIER S, SAXTON K E, PAPENDICK R I,et al. Zero-tillage furrow opener effects on seed environment and wheat emergence[J]. Soil and Tillage Research,1991,21(3-4):347-360.

25 趙喜安. 公路線形的線元設計法[J].公路，1995(3)：24-27.

26 張耀，曹小平，王春芬，等. 材料力學[M].北京：清華大學出版社，2015．

27 GILL W R, VANDEN B, GLEN E. Soil dynamics in tillage and traction[M].Agriculture Hanbook No.316.United States Department of Agriculutre,1967.

28 余泳昌，劉文藝，趙迎芳，等. 立柱式深松鏟受力數(shù)學模型及試驗分析[J]. 農(nóng)業(yè)工程學報，2007，23(6)：109-113.

YU Yongchang, LIU Wenyi, ZHAO Yingfang, et al. Force mathematical model and examination analysis of the column subsoiler[J].Transactions of the CSAE,2007,23(6):109-113. ( in Chinese)

29 許春林，李連豪，趙大勇. 北方大型聯(lián)合整地機設計與試驗[M]. 北京：中國農(nóng)業(yè)大學出版社，2014：45-46.

30 鄭侃，何進，李洪文，等.基于離散元深松土壤模型的折線破土刃深松鏟研究[J/OL]. 農(nóng)業(yè)機械學報，2016，47(9):62-72.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20160910&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2016.09.010.

ZHENG Kan, HE Jin, LI Hongwen, et al. Research on polyline soil-breaking blade subsoiler based on subsoiling soil model using discrete element method[J/OL].Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(9):62-72.(in Chinese)

31 劉俊安，王曉燕，李洪文，等.基于土壤擾動與牽引阻力的深松鏟結(jié)構(gòu)參數(shù)優(yōu)化[J/OL]. 農(nóng)業(yè)機械學報，2017，48(2)：60-67.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20170208&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2017.02.008.

LIU Junan, WANG Xiaoyan, LI Hongwen,et al. Optimization of structural parameters of subsoiler based on soil disturbance and traction resistance[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(2):60-67.(in Chinese)

32 AILI Hasimu, CHEN Ying. Soil disturbance and draft force of selected seed openers[J]. Soil & Tillage Research, 2014, 140:48-54.

33 趙淑紅，劉宏俊，張先民，等. 滑推式開溝器設計與作業(yè)性能優(yōu)化試驗[J]. 農(nóng)業(yè)工程學報，2016，32(19)：26-34.

ZHAO Shuhong, LIU Hongjun, ZHANG Xianmin, et al. Design and optimization experiment of working performance of sliding push opener[J].Transactions of the CSAE,2016,32(19):26-34. (in Chinese)

34 陳玉侖. 稻麥聯(lián)合收獲開溝埋草多功能一體機的設計[D]. 南京：南京農(nóng)業(yè)大學，2009.

CHEN Yulun. Design of harvest ditch and stank-disposing machine[D]. Nanjing:Nanjing Agricultural University, 2011.(in Chinese)

35 賈洪雷，鄭嘉鑫，袁洪方，等. 仿形滑刀式開溝器設計試驗[J]. 農(nóng)業(yè)工程學報，2017，33(4)：16-24.

JIA Honglei, ZHENG Jiaxin, YUAN Hongfang, et al. Design and experiment of profiling sliding-knife opener[J].Transactions of the CSAE,2017,33(4):16-24.(in Chinese)

36 李霞，張東興，王維新，等. 受迫振動深松機性能參數(shù)優(yōu)化與試驗[J]. 農(nóng)業(yè)工程學報，2015，31(21)：17-24.

LI Xia, ZHANG Dongxing, WANG Weixin, et al. Performance parameter optimization and experiment of forecd-vibration subsoiler[J].Transactions of the CSAE,2015,31(21):17-24.(in Chinese)

37 李潮海，周順利. 土壤容重對玉米苗期生長的影響[J].華北農(nóng)學報，1994，9(2)：49-54.

LI Chaohai, ZHOU Shunli. Influence of volume weight of soil on maize seedings growth[J]. Acta Agriculturae Boreali-Sinica, 1994, 9(2)：49-54.(in Chinese)

38 李潮海，梅沛沛，王群，等. 下層土壤容重對玉米植株養(yǎng)分吸收和分配的影響[J].中國農(nóng)業(yè)科學，2007，40(7)：1371-1378.

LI Chaohai, MEI Peipei, WANG Qun,et al. Influences of soil bulk density in deep soil layers on absorption and distribution of nitrogen phosophrous and potassium in maize[J]. Scientia Agricultura Sinica, 2007, 40(7):1371-1378.(in Chinese)