徐文仲，孫玉凱，王娜，張華敏，孟軍*

（1.東北農(nóng)業(yè)大學(xué)文理學(xué)院，哈爾濱 150030；2.東北農(nóng)業(yè)大學(xué)工程學(xué)院，哈爾濱 150030）

馬鈴薯中耕作業(yè)可疏松土壤、消滅雜草，在一定程度上增加馬鈴薯產(chǎn)量、提高馬鈴薯品質(zhì)[1-2]，中耕培土作業(yè)可塑造馬鈴薯生長地壟[3]，良好培土壟型降低土壤容重，增加土壤松散度；調(diào)節(jié)土壤三相比，增加氣液項，增大總空隙度，增加土壤透氣性；調(diào)節(jié)地溫，為馬鈴薯生長提供適宜溫度環(huán)境，減少病蟲害發(fā)生[4-5]，顯著提高馬鈴薯品質(zhì)和產(chǎn)量[6]。

國外對馬鈴薯中耕機械研究較早，理論較成熟，已開發(fā)出自動化程度高、性能可靠的馬鈴薯中耕機，但目前國外機械進口手續(xù)繁雜，各種零部件更換、維修不便，不適宜在國內(nèi)大面積推廣。國內(nèi)專家學(xué)者對馬鈴薯中耕機械展開研究，呂金慶等研究犁鏵式馬鈴薯中耕機，通過繪制犁體曲面、分析培土壟型等方法設(shè)計一種應(yīng)用于砂壤土環(huán)境的犁鏵式馬鈴薯中耕機，分析培土過程影響中耕效果因素，優(yōu)化作業(yè)參數(shù)，作業(yè)效果良好[7]；衣淑娟等設(shè)計一種圓盤式馬鈴薯中耕機，采用一對傾斜安裝的培土圓盤完成培土作業(yè)，通過改變調(diào)節(jié)把手與機具安裝位置改變培土圓盤之間夾角，適應(yīng)不同地塊培土作業(yè)[8]。東北地區(qū)土壤黏重，傳統(tǒng)中耕機在作業(yè)過程中存在碎土率低，培土效果不佳問題，對馬鈴薯中耕機理論和作業(yè)參數(shù)仍需深入研究[9]。

針對上述問題，本文對馬鈴薯中耕機培土器展開研究，通過分析培土機理，確定影響培土效果主要因素，優(yōu)化田間試驗相關(guān)參數(shù)，以期為馬鈴薯中耕機關(guān)鍵部件設(shè)計與調(diào)整提供技術(shù)參考。

1 整機結(jié)構(gòu)與工作原理

1.1 馬鈴薯中耕機整機結(jié)構(gòu)

馬鈴薯中耕機主要由肥箱、機架、地輪、深松鏟、碎土刀、培土犁體、培土擋板等組成，培土犁體和培土擋板共同構(gòu)成培土器完成培土作業(yè)，其結(jié)構(gòu)如圖1所示，4把碎土刀兩兩之間呈90°安裝在旋轉(zhuǎn)單體上，通過U形卡安裝到機架上，每個旋轉(zhuǎn)單體前均安裝有深松鏟，培土犁體安裝在碎土刀后方，在碎土刀切碎土壤后開展培土作業(yè)，培土擋板安裝在培土犁體后，配合培土犁體完成培土作業(yè)，各旋轉(zhuǎn)單體和地輪間距可根據(jù)作業(yè)地區(qū)壟距調(diào)整，通用性較好。

圖1 馬鈴薯中耕機整機結(jié)構(gòu)Fig.1 Overall structure of driving-type potato cultivator

1.2 工作原理

中耕機由拖拉機帶動工作，拖拉機動力輸出軸為中耕機提供動力，工作中碎土刀高速旋轉(zhuǎn)切碎土壤并除草，培土犁體的鏵尖保持斜下朝向，拖拉機驅(qū)動碎土刀高速旋轉(zhuǎn)將土壤切碎，掉落溝底土壤隨中耕機前進沿培土犁體和兩側(cè)擋板向后傾斜運動，并被推向兩側(cè)薯壟[10-12]，碎土刀拋起的土壤運動到培土犁體和培土擋板上反彈覆到土壟。

2 馬鈴薯培土農(nóng)藝要求

通過查閱文獻[13]及相關(guān)資料，明確馬鈴薯中耕培土要求，中耕培土后壟形斷面尺寸見圖2，圖中L為壟距（mm）；h為壟高（mm）；φ為壟壁土壤自然休止角（°）；a1為壟底寬（mm）；a0為壟頂寬（mm）。

圖2 壟形斷面Fig.2 Ridge section

馬鈴薯培土農(nóng)藝要求規(guī)定主要參數(shù)有：壟距L700~900 mm，壟底寬a184~90 mm，壟頂寬a0310~400 mm，壟壁土壤自然休止角φ40°~50°。

3 培土機理分析

碎土刀與培土犁體位置關(guān)系如圖3所示，中耕機通過動力輸出軸驅(qū)動碎土刀旋轉(zhuǎn)將土壤切碎，壟溝底破碎土壤顆粒在碎土刀運動到溝底時隨碎土刀旋轉(zhuǎn)，沿碎土刀邊緣切線方向向后運動至培土犁體和培土擋板上，破碎土壤顆粒以一定初速度運動至培土犁體和培土擋板上反彈落到薯壟上，未被碎土刀拋起的土壤顆粒及由于碎土刀罩殼阻擋而掉落回溝底的土壤則在中耕機前進過程中沿培土犁體兩側(cè)向后傾斜運動上升，達到預(yù)定溝邊時，推向壟中心，土壤顆粒按自然休止角構(gòu)成壟形。

圖3 土壤顆粒運動狀態(tài)Fig.3 Movement state of the soil particle

為滿足培土農(nóng)藝要求，被碎土刀拋起的土壤顆粒將分別運動培土犁體和培土擋板上反彈，運動到培土犁體上土壤顆粒經(jīng)反彈落到薯壟側(cè)面，運動到培土擋板上土壤顆粒經(jīng)反彈落到薯壟上方，在培土犁體和培土擋板共同作用下保證良好培土效果，如圖4所示。

圖4 培土過程原理Fig.4 Schematic of soil-cultivating process

因土壤顆粒在培土犁體和培土擋板上反彈過程類似，因此以土壤顆粒在培土犁體上反彈過程為例展開分析。

由于培土犁體為傾斜平面，因此土壤顆粒與培土犁體接觸后運動為空間運動，實際田間作業(yè)過程顯示，土壤顆粒與培土犁體接觸后反彈落到土壟側(cè)面，以培土犁體犁尖為原點建立空間直角坐標系OXYZ，如圖5所示，培土犁體面為該空間直角坐標系第一象限內(nèi)過原點的平面，假設(shè)培土犁體面一般方程為：

圖5 土壤顆粒運動Fig.5 Soil particle movement

其中，A、B、C均不為0。

則培土犁體面與坐標平面XOY夾角為：

培土擋板面與坐標平面YOZ夾角為：

以與土壤顆粒接觸的培土犁體面為O1X1Y1面建立空間坐標系O1X1Y1Z1，如圖6所示。土壤顆粒從空間坐標系第三象限以速度V運動到培土犁體面發(fā)生反彈，忽略土壤顆粒與培土犁體碰撞過程能量損失，土壤顆粒仍以速度V反彈離開培土犁體面落到薯壟側(cè)面，γ為速度方向與O1X1Y1坐標平面夾角。以反彈速度V在O1X1Y1平面上投影及Z1軸為坐標軸建立平面直角坐標系O1X2Z1，便于分析土壤顆粒反彈后運動狀態(tài)，如圖6所示，γ為土壤顆粒速度方向與坐標系O1X1Y1平面夾角，δ為平面坐標系O1X2軸與水平軸X1之間夾角。

圖6 培土過程Fig.6 Soil cultivation process

土壤顆粒從培土犁體面上反彈后以相同速度發(fā)生斜上拋運動落到薯壟側(cè)面，分析土壤顆粒拋體運動的水平拋出距離，在該過程中土壤顆粒作斜上拋運動，過程中受自身重力和空氣阻力作用[14-16]，如圖7所示。

圖7 斜拋運動Fig.7 Oblique motion

由牛頓運動定律，分別建立土壤顆粒在水平和豎直方向動力學(xué)方程：

式中，z-土壤顆粒在豎直方向位移（mm）；x-土壤顆粒在水平方向位移（mm）；vx-土壤顆粒速度在水平方向分速度（m·s-1）；vz-土壤顆粒速度在豎直方向分速度（m·s-1）；m-土壤顆粒的質(zhì)量（g）；k-阻力系數(shù)。

分離變量求解微分方程（4）和（5）可得：

式中，t-土壤顆粒運動時間。

當(dāng)土壤顆粒落到薯壟上，即z=0，式（5）可整理為：

式（7）為超越方程，無法求出土壤顆粒水平拋出距離與拋出角度解析式，因此采用近似方法求解[17]，

最終得到培土擋板中心位置到土壤顆粒斜拋落點位置距離為：

式中，γ-土壤顆粒速度方向與坐標系O1X1Y1平面夾角（°）；δ-平面坐標系O1X2軸與水平軸X1之間夾角（°）。

由上述分析可得出，土壤顆粒運動主要與土壤顆粒速度、土壤顆粒質(zhì)量、碎土刀與培土器橫向距離l1、土壤顆粒在豎直方向位移、土壤顆粒拋射角度、土壤顆粒拋射方向與水平方向夾角等有關(guān)，而土壤顆粒在豎直方向位移、拋射角度等受碎土刀與培土器縱向距離l2影響。

4 田間試驗

4.1 試驗準備

2021年6月在東北農(nóng)業(yè)大學(xué)阿城試驗基地開展相關(guān)試驗研究。使用馬鈴薯中耕機對指定地塊進行中耕作業(yè)，試驗地塊為旱作區(qū)，土壤類型為東北黑黏土，土壤含水率為9.8%，深度150 mm，平均土壤緊實度1 230 kPa。試驗過程如圖8所示，配套動力為約翰迪爾1354拖拉機（99.2 kW）。

圖8 田間試驗Fig.8 Field experiment

4.2 試驗指標

試驗指標參考JB/T7864-1999《旱田中耕追肥機試驗方法》[18]，選擇除草率和培土斷面合格率作為田間試驗指標。

4.2.1 除草率

試驗前隨機選取3個長1 m×寬1 m區(qū)域，統(tǒng)計雜草數(shù)目，試驗后再次統(tǒng)計該區(qū)域內(nèi)雜草株數(shù)，重復(fù)5次試驗之后取平均值可得除草率,其計算公式如下：

式中，C為除草率（%）；QZ為試驗前區(qū)域內(nèi)雜草株數(shù)；HZ為試驗后區(qū)域內(nèi)雜草株數(shù)。

4.2.2 培土斷面合格率

在中耕區(qū)域每隔2 m選取一個位置測定繪制中耕后薯壟橫向斷面圖，從橫向斷面圖上測量溝底寬、上寬、邊坡角、深度及培土高度，溝底座土及溝壁浮土厚度等數(shù)據(jù)，橫向斷面圖數(shù)量占總測定數(shù)量百分比表示培土斷面合格率，計算方法如式（11）：

式中，η為培土斷面合格率（%）；n為符合農(nóng)藝要求的薯壟橫向斷面數(shù)量；N為薯壟橫向斷面總測量數(shù)。

4.3 試驗方案與結(jié)果分析

4.3.1 試驗方案

正交試驗是實施多因素試驗并尋求最優(yōu)參數(shù)組合的一種方法，本研究利用Design Expert軟件安排二次旋轉(zhuǎn)正交組合試驗。根據(jù)機理分析結(jié)果確定試驗因素為刀軸轉(zhuǎn)速、碎土刀與培土器橫向距離和縱向距離，馬鈴薯中耕培土農(nóng)藝要求壟距L為700~900 mm，中耕后壟底寬a1為84~90 mm，根據(jù)實測薯壟橫向斷面數(shù)據(jù)及馬鈴薯農(nóng)藝要求，確定行距L為900 mm，中耕后壟底寬a1為85 mm，結(jié)合土壤顆粒拋出距離表達式及中耕機實際工況，確定土壤顆粒從碎土刀與培土器橫向距離l1范圍為510~550 mm，結(jié)合機理分析及中耕深度要求，確定碎土刀與培土器縱向距離l2為200~260 mm，通過查閱參考文獻[19]得到碎土刀軸較優(yōu)轉(zhuǎn)速為200~300 r·min-1，確定試驗指標為除草率和培土斷面合格率，根據(jù)二次旋轉(zhuǎn)正交組合試驗相關(guān)知識，計算γ值后確定各試驗因素水平，表1為試驗因素水平編碼表，試驗方案及結(jié)果，如表2所示。

表1 水平編碼Table 1 Horizontal coding

表2 試驗方案與結(jié)果Table 2 Test plan and experimental data

4.3.2 結(jié)果分析

對試驗結(jié)果進行方差分析，確定各試驗因素與試驗指標相關(guān)程度，并采用多元回歸擬合建立試驗因素與指標之間數(shù)學(xué)模型[20]，得到除草率C和培土斷面合格率η兩個試驗指標回歸方程。

①除草率C方差分析。

使用Design-Expert軟件對除草率C進行方差分析，得到方差分析結(jié)果如表3所示，其中，各模型對試驗指標除草率C影響顯著性順序如下：x1、x3、x2x3、x2、x22、x32、x1x2、x1x3、x12；其中刀軸轉(zhuǎn)速x1、縱向距離x3對除草率C影響極顯著（P＜0.01）；橫向距離x2、橫向距離和縱向距離交互項x2x3對除草率C影響顯著（0.01＜P＜0.05）；刀軸轉(zhuǎn)速和橫向距離交互項x1x2、橫向距離二次項x22和縱向距離二次項x32影響較顯著（0.05＜P＜0.1）；剩余因素影響不顯著（P＞0.1）。將模型中不顯著項并入殘差項，得到方差分析結(jié)果如表3所示，同時得到除草率C與各顯著因素及其交互作用回歸模型如式（12）所示。

表3 除草率方差分析Table 3 Variance analysis for the weeding rate

由除草率方差分析表可得，失擬項P=0.2936，不顯著（P＞0.1），可知不存在其他因素影響試驗效果，上述分析合理。

②培土斷面合格率η方差分析。

使用Design-Expert軟件對培土斷面合格率η進行方差分析，結(jié)果如表4所示，其中，各模型對于試驗指標培土斷面合格率η影響顯著性順序如下：x3、x1、x2、x22、x32、x2x3、x12、x1x3、x1x2；其 中，刀軸轉(zhuǎn)速x1、縱向距離x3對培土斷面合格率η影響極顯著（P＜0.01）；橫向距離x2、橫向距離的二次項x22和縱向距離二次項x32對培土斷面合格率η影響顯著（0.01＜P＜0.05）；橫向距離和縱向距離交互項x2x3對培土斷面合格率η影響較顯著（0.05＜P＜0.1）；剩余因素模型影響不顯著（P＞0.1）。將模型中不顯著項并入殘差項，得到方差分析結(jié)果如表4所示，同時得到培土斷面合格率η與各顯著因素及其交互作用回歸模型，如式（13）所示。

表4 培土斷面合格率η方差分析Table 4 Variance analysis for the pass rate of ridge section

由培土斷面合格率η方差分析表可得，失擬項P=0.3498，不顯著（P＞0.1），可知不存在其他因素影響試驗效果，上述分析合理。

4.3.3 響應(yīng)曲面分析

通過響應(yīng)曲面分析法分析各試驗因素間顯著交互作用對試驗指標的影響情況，利用Design-Expert軟件得出各試驗因素對試驗指標影響的響應(yīng)曲面，如圖9所示。

對試驗指標除草率C，當(dāng)碎土刀與培土器縱向距離為230 mm時，刀軸轉(zhuǎn)速與橫向距離交互作用響應(yīng)曲面如圖9a所示，橫向距離一定時，除草率隨刀軸轉(zhuǎn)速增加呈逐漸增加趨勢，最優(yōu)刀軸轉(zhuǎn)速范圍為250~279.7 r·min-1；當(dāng)?shù)遁S轉(zhuǎn)速一定時，除草率C隨橫向距離增加呈先升后降趨勢，最優(yōu)橫向距離在524.9~535.1 mm范圍內(nèi)，其中，刀軸轉(zhuǎn)速是影響除草率主要試驗因素。

對試驗指標培土斷面合格率η，當(dāng)?shù)遁S轉(zhuǎn)速為250 r·min-1時，橫向距離與縱向距離間交互作用關(guān)系如圖9b所示，當(dāng)橫向距離一定時，培土斷面合格率η隨縱向距離增加呈先增后減趨勢，縱向距離最優(yōu)范圍為230~247.8 mm；當(dāng)縱向距離一定時，培土斷面合格率η隨橫向距離增加呈先增后減趨勢，橫向距離最優(yōu)范圍為524.9~538.5 mm，其中，培土斷面合格率主要影響因素是縱向距離。

圖9 除草率和培土斷面合格率雙因素響應(yīng)曲面Fig.9 Response surface of double parameters about weeding rate and pass rate of ridge section

4.3.4 參數(shù)優(yōu)化

根據(jù)上述響應(yīng)曲面分析結(jié)果，可得刀軸轉(zhuǎn)速x1較優(yōu)范圍為250~279.7 r·min-1，橫向距離x2較優(yōu)范圍為524.9~535.1 mm，縱向距離x3較優(yōu)范圍為230~247.8 mm，利用Design-Expert軟件優(yōu)化功能求解最優(yōu)參數(shù)組合，根據(jù)上述分析并結(jié)合田間作業(yè)要求確定優(yōu)化約束條件為：

通過優(yōu)化求解，得到刀軸轉(zhuǎn)速為270.6 r·min-1，橫向距離為529.2 mm，縱向距離為240.8 mm時，馬鈴薯中耕機作業(yè)除草率和培土斷面合格率最優(yōu)，除草率為96.7%，培土斷面合格率為95.6%。

4.4 田間驗證試驗

根據(jù)優(yōu)化最優(yōu)參數(shù)組合，檢驗優(yōu)化后馬鈴薯中耕機作業(yè)質(zhì)量，確定各試驗因素參數(shù)如下：刀軸轉(zhuǎn)速為270 r·min-1，碎土刀與培土器橫向距離為530 mm，碎土刀與培土器縱向距離為240 mm。根據(jù)國家標準JB/T7864-1999中規(guī)定檢驗指標，確定試驗指標為除草率、培土斷面合格率、碎土率、傷苗率，試驗結(jié)果如表5所示，其中，碎土率和傷苗率測量和計算方式如下：

在每一個行程隨機取一個測區(qū)，稱量中耕后0.5 m×0.5 m面積、0.2 m耕層內(nèi)直徑＜25 mm土塊質(zhì)量和土壤總質(zhì)量，以直徑＜25 mm土塊質(zhì)量占總質(zhì)量百分比表示碎土率，共測3個行程，取平均值。計算方法如式（15）所示。

式中，ψ為碎土率（%）；WS為測區(qū)內(nèi)土壤總質(zhì)量（g）；WZ為測區(qū)內(nèi)直徑＜25 mm土塊質(zhì)量（g）。

在試驗區(qū)內(nèi)往返行程上各取2點，在機具作業(yè)幅寬內(nèi)，長度為1 m區(qū)域內(nèi)測定該面積內(nèi)總株數(shù)。中耕作業(yè)后在測定區(qū)域內(nèi)統(tǒng)計傷苗株數(shù)，隨機重復(fù)5次，取平均值，計算方法如式（16）所示。

式中，s為傷苗率（%）；MZ為測定區(qū)域內(nèi)總苗數(shù)；MZ為測定區(qū)域內(nèi)傷苗總株數(shù)。

優(yōu)化后馬鈴薯中耕機除草率為96.4%，培土斷面率為95.2%，碎土率為95.3%，傷苗率為1.2%，作業(yè)性能滿足馬鈴薯中耕培土農(nóng)藝要求。

5 結(jié)論

根據(jù)JB/T7864-1999《旱田中耕追肥機試驗方法》確定試驗指標，開展試驗研究，通過方差分析和響應(yīng)曲面分析，建立各試驗因素與試驗指標之間數(shù)學(xué)模型，求出最優(yōu)作業(yè)參數(shù)組合，以優(yōu)化后作業(yè)參數(shù)進行田間試驗，檢驗參數(shù)優(yōu)化后馬鈴薯中耕機作業(yè)效果，結(jié)果表明，當(dāng)馬鈴薯中耕機刀軸轉(zhuǎn)速為270 r·min-1、碎土刀與培土器的橫向距離為530 mm、碎土刀與培土器縱向距離為240 mm時，除草率為96.4%，培土斷面合格率為95.2%，碎土率為95.3%，傷苗率為1.2%，可滿足馬鈴薯中耕培土要求。