賈洪雷，鄭 健，趙佳樂※，郭明卓，莊 健，王增輝

（1.吉林大學(xué)工程仿生教育部重點實驗室，長春 130025；2.吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院，長春 130025；3.東北師范大學(xué)人文學(xué)院，長春 130117）

0 引 言

東北地區(qū)大豆種植面積和產(chǎn)量約占中國的40%以上[1]，但東北地區(qū)地處寒溫帶，氣候干燥易旱，冬季土壤封凍期較長，春季耕種時一直存在土壤耕種層積溫不夠和含水率不足，嚴(yán)重影響了大豆根系生長，延緩了大豆的生長發(fā)育。因此，相比于水稻和玉米等高產(chǎn)作物，大豆產(chǎn)量較低，同時受到進口轉(zhuǎn)基因大豆的沖擊，造成大豆種植面積逐年遞減，2012-2016年已減少13.4%以上[2]。因此，通過提升土壤耕種層溫度和含水率提高中國大豆產(chǎn)量，是保障中國大豆生產(chǎn)良性發(fā)展的有效途徑。

土壤溫度和含水率作為土壤主要物理參量之一，影響著土壤中水鹽的運移，微生物的數(shù)量與活性，農(nóng)田溫室氣體排放，土壤碳氮等物質(zhì)的化學(xué)和生物學(xué)過程及平衡[3-6]，進而影響種子的萌發(fā)和作物的生長。土壤溫度和含水率除主要受太陽輻射、大氣循環(huán)和降水等外部因素影響以外，還受到土壤自身孔隙結(jié)構(gòu)的影響[7-8]。土壤孔隙主要分為大孔隙和毛細(xì)孔隙[9-10]，Holden等[11]研究得出增加土壤大孔隙的數(shù)量，可加快土壤與外界的物質(zhì)能量交換，可有效提高地溫，但也會加速表層土壤中的水分蒸發(fā)；Skopp等[12]研究得出增加土壤毛細(xì)孔隙的數(shù)量，利用其毛細(xì)作用吸收耕種層下方土壤中的蓄水，可提高耕種層土壤含水率；因此，通過構(gòu)建出合理的土壤孔隙結(jié)構(gòu)，可提升東北地區(qū)春播時的耕種層土壤溫度和含水率。

農(nóng)業(yè)機械中的耕作部件可改善土壤孔隙結(jié)構(gòu)，Haukka等[13]研究得出對土壤進行疏松作業(yè)可有效提高土壤大孔隙數(shù)量；Asare等[14-15]研究得出對土壤進行鎮(zhèn)壓作業(yè)可增加土壤毛細(xì)孔數(shù)量，有效提高土壤蓄水能力。因此，本文設(shè)計一種對表層土壤進行疏松作業(yè)，對下層土壤進行鎮(zhèn)壓作業(yè)，集松土、溝底鎮(zhèn)壓和開溝功能為一體的耕整機構(gòu)，使耕種層土壤達(dá)到表層松軟，下層緊實的狀態(tài)，從而加速耕種層土壤對外部熱量和下層蓄水的吸收，達(dá)到迅速提升耕種層土壤溫度和含水率的目的。

近年來，仿生學(xué)在農(nóng)業(yè)機械設(shè)計中的應(yīng)用愈發(fā)廣泛，如趙淑紅等[16]根據(jù)旗魚頭部的流線型曲線，設(shè)計了仿旗魚頭部曲線型開溝器，優(yōu)化了作業(yè)效果；馬云海等[17]依據(jù)土壤動物與生俱來的減黏脫土特性，設(shè)計了仿生波紋形開溝器，具有較強的減黏降阻性能；王文君[18]研制的仿形仿生鎮(zhèn)壓裝置，可有效提高土壤含水率。本文本著師法自然的思想，在自然界尋找可以改善土壤孔隙結(jié)構(gòu)的生物。蚯蚓是土壤中最常見的生物之一，Lavelle等[19]的研究表明其特殊的運動方式可有效疏松土壤，改善土壤孔隙結(jié)構(gòu)的作用。

基于上述原因，本文的目的是設(shè)計出一種設(shè)有仿蚯蚓動態(tài)松土開溝機構(gòu)和鎮(zhèn)壓力調(diào)節(jié)機構(gòu)的仿蚯蚓運動多功能開溝器，其可模擬出蚯蚓疏松土壤的生物學(xué)行為（主要包括結(jié)構(gòu)和運動方式），并通過理論分析和試驗優(yōu)化等方法，明確仿蚯蚓運動多功能開溝器的結(jié)構(gòu)參數(shù)對中國東北地區(qū)土壤溫度、含水率和大豆生長發(fā)育的影響。

1 仿蚯蚓運動多功能開溝器設(shè)計及其作業(yè)原理

1.1 生物原型分析

如圖1所示，蚯蚓的體壁具有環(huán)肌和縱肌2種肌肉層[20-21]，其主要依靠體壁肌肉層和剛毛的密切配合在土壤中運動[22-25]。蚯蚓運動時，首先將前端的口前葉伸入土壤縫隙并膨脹，像楔子一樣插入土中，后部體節(jié)的剛毛伸出并插入土壤，起到支撐和固定的作用，而后再收縮前端的環(huán)肌產(chǎn)生向前的推力，使前部體節(jié)向前推進，呈收縮狀態(tài)（舒張環(huán)肌層，收縮縱肌層）。環(huán)節(jié)的收縮沿著體節(jié)一節(jié)一節(jié)地向后傳輸，使蚯蚓身體后端向前拉，呈舒張狀態(tài)（舒張縱肌層，收縮環(huán)肌層），環(huán)肌與縱肌協(xié)調(diào)地收縮使蚯蚓產(chǎn)生波浪式蠕動[26-28]。

圖1 蚯蚓的運動過程Fig.1 Motion process of earthworm

1.2 仿蚯蚓運動多功能開溝器設(shè)計與作業(yè)原理

圖2為仿蚯蚓運動多功能開溝器結(jié)構(gòu)示意圖。如圖2a所示，仿蚯蚓運動多功能開溝器由仿蚯蚓動態(tài)松土開溝機構(gòu)、鎮(zhèn)壓力調(diào)節(jié)機構(gòu)、波紋圓盤、軸、機架等組成。如圖2b所示，密封擋板1、波紋圓盤1、軸、鎮(zhèn)壓力調(diào)節(jié)機構(gòu)、波紋圓盤2和橡膠擋板2依次焊合。仿蚯蚓動態(tài)松土開溝機構(gòu)的內(nèi)圓盤和外圓盤的焊合面則與軸上波紋圓盤的安裝位置中心相對應(yīng)，并與滾筒焊合，軸的兩端與機架通過螺母固連。如圖2c所示，機具作業(yè)時，在土壤阻力的作用下，仿蚯蚓動態(tài)松土開溝機構(gòu)做逆時針轉(zhuǎn)動，而波紋圓盤和鎮(zhèn)壓力調(diào)節(jié)機構(gòu)做平移運動。

1.2.1 鎮(zhèn)壓力調(diào)節(jié)機構(gòu)設(shè)計

圖3為鎮(zhèn)壓力調(diào)節(jié)機構(gòu)示意圖。如圖3所示，4個鎮(zhèn)壓力調(diào)節(jié)機構(gòu)的彈簧一端套接于軸上，另一端設(shè)有滾輪，滾輪與套筒相接觸，4個鎮(zhèn)壓力調(diào)節(jié)機構(gòu)于軸上呈90°交錯對稱排列。鎮(zhèn)壓力調(diào)節(jié)機構(gòu)的滾輪安裝在設(shè)置于仿蚯蚓動態(tài)松土開溝機構(gòu)套筒內(nèi)的導(dǎo)軌槽中，通過調(diào)節(jié)軸上的螺母位置，可調(diào)節(jié)仿蚯蚓運動多功能開溝器對土壤的鎮(zhèn)壓力，其鎮(zhèn)壓力由式（1）得出，通過仿蚯蚓運動多功能開溝器對種溝底部的鎮(zhèn)壓作用，提高耕種層下方土壤的堅實度，進而提高耕種層下方土壤內(nèi)的毛細(xì)作用，提高耕種層土壤的含水率。Celik等[29]研究發(fā)現(xiàn)選用60 kPa的壓實力，可有效提高土壤含水量，因此本文選用的溝底鎮(zhèn)壓強度設(shè)計為60 kPa。

圖2 仿蚯蚓運動多功能開溝器示意圖Fig.2 Schematic diagram of earthworm-like multi-function opener

圖3 鎮(zhèn)壓力調(diào)節(jié)機構(gòu)示意圖Fig.3 Schematic diagram of pressure regulating mechanism

式中m外為外圓盤質(zhì)量，kg；m內(nèi)為內(nèi)圓盤質(zhì)量，kg；m套為套筒質(zhì)量，kg；k為彈簧彈性系數(shù)，N/m；x為彈簧原長，m；x1、x2為彈簧受力壓縮后的長度，m。

1.2.2 仿蚯蚓動態(tài)松土開溝機構(gòu)設(shè)計

圖4為仿蚯蚓動態(tài)松土開溝機構(gòu)結(jié)構(gòu)示意圖。如圖4所示，仿蚯蚓動態(tài)松土開溝機構(gòu)由外圓盤、內(nèi)圓盤、套筒、內(nèi)圓盤、外圓盤焊合而成，外圓盤上設(shè)有橡膠擋板。每個外圓盤和內(nèi)圓盤表面均設(shè)有12個仿蚯蚓柔性凸表面機構(gòu)，每個仿蚯蚓柔性凸表面機構(gòu)均設(shè)有1個驅(qū)動機構(gòu)，驅(qū)動機構(gòu)一端與仿蚯蚓柔性凸表面機構(gòu)相套接，另一端設(shè)有滾輪，滾輪與波紋圓盤表面相接觸。

圖4 仿蚯蚓動態(tài)松土開溝機構(gòu)結(jié)構(gòu)示意圖Fig.4 Schematic diagram of earthworm-like dynamic scarification furrow mechanism

圖5為仿蚯蚓柔性凸表面機構(gòu)運動示意圖。如圖5所示，機具作業(yè)時驅(qū)動機構(gòu)與波紋圓盤發(fā)生相對轉(zhuǎn)動，當(dāng)驅(qū)動機構(gòu)由波紋圓盤的波谷向波峰運動時，仿蚯蚓柔性凸表面機構(gòu)可模擬出蚯蚓環(huán)肌擴張，同時縱肌收縮的運動狀態(tài)；當(dāng)驅(qū)動機構(gòu)由波紋圓盤的波峰向波谷運動時，仿蚯蚓柔性凸表面機構(gòu)可模擬出蚯蚓環(huán)肌收縮，同時縱肌擴張的運動狀態(tài)；隨著驅(qū)動機構(gòu)在波紋圓盤表面上循環(huán)往復(fù)的逆時針轉(zhuǎn)動，即可驅(qū)使仿蚯蚓動態(tài)松土開溝機構(gòu)在完成開溝作業(yè)的同時，模擬出蚯蚓疏松土壤的生物學(xué)行為，增強了對土壤的疏松作用，加速了耕種層土壤與外部環(huán)境之間的物質(zhì)能量交換，進而更快速地提高地溫。

圖5 仿蚯蚓柔性凸表面機構(gòu)運動示意圖Fig.5 Motion schematic diagram of earthworm-like flexible convex surface mechanism

1.2.3 仿蚯蚓柔性凸表面機構(gòu)布置結(jié)構(gòu)設(shè)計

為降低增設(shè)仿蚯蚓柔性凸表面機構(gòu)對仿蚯蚓運動多功能開溝器作業(yè)阻力的影響，每個仿蚯蚓動態(tài)松土開溝機構(gòu)表面上的12個仿蚯蚓柔性凸表面機構(gòu)呈阿基米德螺旋線分布，從而在最大程度上減小仿蚯蚓運動多功能開溝器的工作阻力。阿基米德螺旋線標(biāo)準(zhǔn)方程為

式中b為阿基米德螺旋線系數(shù)，表示每旋轉(zhuǎn)1°時極徑的增加(或減小)量；θ為極角，表示阿基米德螺旋線轉(zhuǎn)過的總度數(shù)；a為當(dāng)θ=0時的極徑。

圖6為仿蚯蚓柔性凸表面機構(gòu)布置示意圖。如圖6a所示，仿蚯蚓動態(tài)松土開溝機構(gòu)的最大半徑值選取為最常見的單元盤開溝器結(jié)構(gòu)參數(shù)值135 mm[30]，大豆的播種深度為50 mm。為保證作業(yè)時12組仿蚯蚓柔性凸表面機構(gòu)能夠全部入土作業(yè)，仿蚯蚓柔性凸表面機構(gòu)布置方案如圖6b所示。仿蚯蚓柔性凸表面機構(gòu)于阿基米德螺旋線上呈30°間隔依次排列，當(dāng)阿基米德螺旋線轉(zhuǎn)過630o（即A處）時放置第一組仿蚯蚓柔性凸表面機構(gòu)，其距外圓盤軸心100 mm；當(dāng)阿基米德螺旋線轉(zhuǎn)過990o（即B處）時，此處距圓盤軸心120 mm。

圖6 仿蚯蚓柔性凸表面機構(gòu)布置示意圖Fig.6 Layout diagram of earthworm-like flexible convex surface mechanism

又由式（2）可得參數(shù)

又知阿基米德螺旋線平面直角坐標(biāo)方程為

因此由式（2）、（3）、（4）可得仿蚯蚓柔性凸表面機構(gòu)在仿蚯蚓動態(tài)松土開溝機構(gòu)表面上的布置結(jié)構(gòu)面直角坐標(biāo)曲線方程為

1.2.4 波紋曲面結(jié)構(gòu)設(shè)計

仿蚯蚓運動多功能開溝器作業(yè)時，仿蚯蚓柔性凸表面機構(gòu)的驅(qū)動機構(gòu)在波紋圓盤的軸向方向上進行往復(fù)運動，為使驅(qū)動機構(gòu)的往復(fù)運動既無剛性沖擊亦無柔性沖擊，應(yīng)使驅(qū)動機構(gòu)運動至波紋圓盤波峰或波谷位置時的速度和加速度值均為0，因此驅(qū)動機構(gòu)的往復(fù)運動方式應(yīng)為正弦加速度運動，如圖7所示，由此得出波紋圓盤凸起輪廓曲線方程應(yīng)為

圖7 波紋圓盤的輪廓曲線示意圖Fig.7 Outline curve diagram of corrugated disk

式中H為波紋圓盤最大凸起高度（波峰與波谷之間的高度差）。

圖8為波紋圓盤結(jié)構(gòu)示意圖。如圖8所示，本文所設(shè)計波紋圓盤直徑應(yīng)小于外圓盤直徑，因此波紋圓盤直徑設(shè)置為250 mm。又可知

式中R1為波紋圓盤半徑，mm，機具中R1=125 mm；N為波紋圓盤波峰數(shù)量。

圖8 波紋圓盤結(jié)構(gòu)示意圖Fig.8 Schematic diagram of corrugated disk

因此由式（6）、（7）可知波紋圓盤的輪廓曲線方程為

仿蚯蚓柔性凸表面機構(gòu)作業(yè)時凸起頻率可由式（9）、（10）與（11）計算。

式中T為仿蚯蚓柔性凸表面機構(gòu)作業(yè)時凸起周期；R2為外圓盤半徑，mm，機具中R2=135 mm；V為機具作業(yè)速度，m/s；N為紋圓盤波峰數(shù)量；ω為仿蚯蚓動態(tài)松土開溝機構(gòu)角速度，rad/s。

因此由式（9）、（10）和（11）可知仿蚯蚓柔性凸表面機構(gòu)作業(yè)時凸頻率f為

由式（8）可知，H值越大，仿蚯蚓柔性凸表面機構(gòu)凸起高度越大；由式（12）可知，在機具作業(yè)速度V一定時，N決定了仿蚯蚓柔性凸表面機構(gòu)作業(yè)時的凸起頻率f，N越大，頻率越高。顯而易見H值越大，N值越大，仿蚯蚓運動多功能開溝器對土壤的疏松程度越高，對耕種層土壤的疏松程度更強，耕種層土壤的大孔隙度越高，土壤溫度提升越快，耕種層土壤水分蒸發(fā)越嚴(yán)重。因此，通過設(shè)置不同的H、N值，可得到不同結(jié)構(gòu)參數(shù)的波紋圓盤，形成不同的耕種層土壤結(jié)構(gòu)，進而產(chǎn)生不同的土壤物理性狀(溫度和含水率)，因此本文選取最大凸起高度H和波紋圓盤波峰數(shù)量N為試驗指標(biāo)進行參數(shù)優(yōu)化試驗，進而完成對波紋圓盤曲面結(jié)構(gòu)的設(shè)計。

2 材料與試驗

2.1 試驗條件與儀器設(shè)備

試驗時間為2017年5月1日-7月10日，試驗地點為吉林大學(xué)教學(xué)科研基地（43°55￠0.262N，125°16￠8.362E）。試驗田為大豆玉米輪作地，2016年作物為玉米，秋季收獲后進行了秸稈覆蓋還田作業(yè)，2017年春季播種前地表覆蓋有玉米秸稈，土壤緊實度為40.2 kPa，土壤容積密度為1.237 g/cm3，土壤含水率為21.2%，土壤溫度為13.2℃，全氮質(zhì)量分?jǐn)?shù)為0.13%，速效鉀和速效磷質(zhì)量分?jǐn)?shù)分別為173.2和16.5 mg/mg。2017年春季選用的大豆種子品種為禾豐50，種植模式為壟上雙行交錯播種，隨機選定3個試驗小區(qū)，每個試驗小區(qū)的寬度為機具作業(yè)幅寬，長度為20 m，離地頭5 m以上。

試驗主要儀器設(shè)備如下：1臺CASE210拖拉機（動力為154 kW，最大行進速度為40 km/h）、1臺2BGD-6大豆耕播機（黑龍江省勃農(nóng)興達(dá)機械有限公司）、卷尺、耕深尺、環(huán)刀組件（100 cm3）、MP-5002型電子天平(上海精密儀器儀表有限公司)、MS-350型水分測定儀(青島拓科儀器有限公司)、SC-900型土壤緊實度儀(美國Spectrum技術(shù)公司)、11 000型土壤溫度計(美國Deltatrak技術(shù)公司)等。

2.2 試驗方法

本文試驗分為2個部分，第一部分為研究仿蚯蚓運動多功能開溝器結(jié)構(gòu)參數(shù)對平均土壤溫度（0～100 mm）、平均土壤含水率（0～100 mm）和大豆平均出苗時間的影響，由本文的分析可知仿蚯蚓柔性凸表面機構(gòu)作業(yè)時的凸起高度H和波紋圓盤波峰數(shù)量N決定了仿蚯蚓運動多功能開溝器對土壤的疏松能力，因此本部分試驗選取的3個試驗因素為作業(yè)速度、凸表面機構(gòu)作業(yè)時的凸起高度H和波紋圓盤波峰數(shù)量N。試驗指標(biāo)為平均土壤溫度（0～100 mm）、平均土壤含水率（0～100 mm）和大豆平均出苗時間，試驗因素編碼表如表1所示。運用Design-Expert軟件中Box-Behnken Design響應(yīng)曲面設(shè)計法，共進行17組試驗，其中12組為析因點，5組為零點以估計誤差，并對試驗結(jié)果進行分析，檢驗各因素對試驗指標(biāo)的顯著性影響和各因素之間是否具有交互作用，并得出響應(yīng)曲面和回歸方程，優(yōu)化出各因素的最佳參數(shù)組合；第二部分為仿蚯蚓運動多功能開溝器結(jié)構(gòu)參數(shù)的驗證對比試驗。根據(jù)第一部分試驗所得出最優(yōu)參數(shù)組合加工出新的仿蚯蚓運動多功能開溝器，并選取傳統(tǒng)的單元盤開溝器和雙圓盤開溝器，并對3個開溝裝置進行田間對比試驗，對最優(yōu)結(jié)構(gòu)參數(shù)條件下的仿蚯蚓運動多功能開溝器相較于傳統(tǒng)開溝器的性能優(yōu)越性進行驗證。

表1 因素水平編碼Table 1 Coding list of factors and levels

2.3 試驗指標(biāo)的測試方法

根據(jù)JB/T10295-2001的相關(guān)規(guī)定，本試驗以土壤平均含水率（0～100 mm）、平均土壤溫度（0～100 mm）、平均出苗時間作為響應(yīng)指標(biāo)。

1）平均土壤含水率（0～100 mm）測定方法。播種15天后，于試驗小區(qū)的對角線上，隨機取樣5點，于每個測定點處使用MS-350水分測定儀測定0～100 mm深度的平均土壤含水率，重復(fù)3次，取其平均值作為此點的0～100 mm平均土壤含水率。取5點的平均值作為試驗小區(qū)0～100 mm平均土壤含水率。每個試驗小區(qū)重復(fù)一次。

2）平均土壤溫度（0～100 mm）測定方法。播種15天后，于0～100 mm平均土壤含水率測定點上測定其0～100 mm平均土壤溫度。使用11 000型土壤溫度計測定0～100 mm深度的平均土壤溫度，重復(fù)3次，取其平均值作為此點的0～100 mm平均土壤溫度。同樣取5點的平均值作為試驗小區(qū)0～100 mm平均土壤含水率。每個試驗小區(qū)重復(fù)一次。

3）大豆平均出苗時間測定方法。從出苗開始到出苗結(jié)束期間，需要每2 d測量一次5 m內(nèi)苗數(shù)，重復(fù)3次，每個試驗小區(qū)重復(fù)一次。平均出苗時間MET由式（13）計算[31]。

式中Nn為從前一個時間點算起的出苗數(shù)，n=1,...,n；Tn為播種后的天數(shù)。

3 試驗結(jié)果與討論

3.1 參數(shù)優(yōu)化試驗結(jié)果

參數(shù)優(yōu)化試驗現(xiàn)場如圖9所示。以各影響因素水平編碼值為自變量，以平均土壤含水率（0～100 mm）、平均土壤溫度（0～100 mm）和平均出苗時間為響應(yīng)指標(biāo)的試驗結(jié)果，如表2所示。

圖9 田間參數(shù)優(yōu)化試驗Fig.9 Parameter optimization of field experiment

表2 響應(yīng)曲面試驗設(shè)計與結(jié)果Table 2 Experiment design and result of response

對表2的數(shù)據(jù)經(jīng)Design-Expert軟件進行二次多元回歸擬合，得到平均土壤含水率（0～100 mm）Y1、平均土壤溫度（0～100 mm）Y2和平均出苗時間Y3對編碼自變量的二次多元回歸方程分別為

對于平均土壤含水率（0～100 mm）對編碼自變量的二次多元回歸方程，回歸診斷顯示，因素對對平均土壤含水率（0～100 mm）的影響均為顯著項(P<0.05)。失擬項P值為0.192，顯然P>0.05，回歸方程不失擬。決定系數(shù)與校正決定系數(shù)均接近于1，表明回歸方程擬合度很高，精密度值為12.528>4，表明該回歸方程在設(shè)計域內(nèi)預(yù)測性能良好。各因素對土壤平均含水率（0～100 mm）顯著性的影響順序依次為：波峰數(shù)量、凸起高度、作業(yè)速度。

對于平均土壤溫度（0～100 mm）對編碼自變量的二次多元回歸方程，回歸診斷顯示，因素對平均土壤溫度（0～100 mm）的影響均為顯著項(P<0.05)。失擬項P值為0.137 6，顯然P>0.05，回歸方程不失擬。決定系數(shù)與校正決定系數(shù)均接近于1，表明回歸方程擬合度很高，精密度值為10.108>4，表明該回歸方程在設(shè)計域內(nèi)預(yù)測性能良好。各因素對土壤（0～100 mm）平均含水率顯著性的影響順序依次為：波峰數(shù)量、作業(yè)速度、凸起高度。

對于平均出苗時間對編碼自變量的二次多元回歸方程，回歸診斷顯示，因素對平均出苗時間的影響均為顯著項(P<0.05)。失擬項P值為1.06，顯然P>0.05，回歸方程不失擬。決定系數(shù)與校正決定系數(shù)均接近于1，表明回歸方程擬合度很高，精密度值為15.526>4，表明該回歸方程在設(shè)計域內(nèi)預(yù)測性能良好。各因素對平均出苗時間顯著性的影響順序依次為：凸起高度、波峰數(shù)量、作業(yè)速度。

3.2 參數(shù)優(yōu)化試驗結(jié)果分析與討論

固定3因素中的1個因素為零水平，應(yīng)用響應(yīng)曲面法分析其他2個因素對試驗指標(biāo)的影響和交互作用。運用Matlab軟件對Design-Expert軟件求出的回歸方程進行尋優(yōu)，得出最佳參數(shù)組合。

3.2.1 試驗各因素對平均土壤含水率（0～100 mm）的影響規(guī)律

運用Design-Expert 6.0.10軟件得到作業(yè)速度V、凸起高度H和波峰數(shù)量N對平均土壤含水率（0～100 mm）影響的響應(yīng)曲面，如圖10所示。

由圖10可知，三因素均對平均土壤含水率（0～100mm）具有顯著性影響（P<0.01），作業(yè)速度分別與凸起高度和波峰數(shù)量具有交互作用，凸起高度和波峰數(shù)量之間無交互作用。平均土壤含水率（0～100 mm）隨著前進速度的加快而下降，且下降趨勢逐步減緩；隨著凸起高度的提升而下降，且下降趨勢逐步減緩；隨著波峰數(shù)量的增加而下降，且下降趨勢逐步減緩。

3.2.2 試驗各因素對平均土壤溫度（0～100 mm）的影響規(guī)律

運用Design-Expert 6.0.10軟件得到作業(yè)速度V、凸起高度H和波峰數(shù)量N對平均土壤溫度（0～100 mm）影響的響應(yīng)曲面，如圖11所示。

圖10 各因素對平均土壤含水率（0～100 mm）影響的響應(yīng)曲面Fig.10 Response surfaces influenced by various factors on average soil moisture content(0-100 mm)

由圖11可知，三因素均對平均土壤溫度（0～100mm）具有顯著性影響（P<0.01），作業(yè)速度分別與凸起高度和波峰數(shù)量具有交互作用，凸起高度和波峰數(shù)量之間無交互作用。平均土壤溫度（0～100 mm）隨著前進速度的加快而上升，且上升趨勢逐步減緩；隨著凸起高度的提升而上升，且上升趨勢逐步減緩；隨著波峰數(shù)量的增加而上升，且上升趨勢逐步減緩。

3.2.3 試驗各因素對大豆平均出苗時間的影響規(guī)律及參數(shù)優(yōu)化

運用 Design-Expert 6.0.10軟件得到作業(yè)速度V、凸起高度H和波峰數(shù)量N對大豆平均出苗時間影響的響應(yīng)曲面，如圖12所示。

由圖12可知，三因素均對大豆平均出苗時間具有顯著性影響（P<0.01），作業(yè)速度分別與凸起高度和波峰數(shù)量具有交互作用，凸起高度和波峰數(shù)量之間無交互作用。隨著前進速度的加快，大豆平均出苗時間先下降而后上升；隨著凸起高度的提升，大豆平均出苗時間先下降而后上升；隨著波峰數(shù)量的增加，大豆平均出苗時間先下降而后上升。

依據(jù)試驗結(jié)果分析以及模型擬合，運用Matlab軟件對式（16）進行尋優(yōu)求解，以平均出苗時間最早為主要尋優(yōu)指標(biāo)，獲得最優(yōu)參數(shù)取值方案：作業(yè)速度為8.55 km/h，凸起高度為6.23 mm，波峰數(shù)量為14.42，此時平均出苗時間最短（9.62 d），平均土壤含水率（0～100 mm）為18.07%，平均土壤溫度（0～100 mm）為16.5℃。由于作業(yè)速度易受田間復(fù)雜地況的影響而較難保證其精準(zhǔn)性，且規(guī)定波峰數(shù)量N只能為自然數(shù)，并綜合考慮加工精度等問題，最終選取的最優(yōu)參數(shù)組合為：作業(yè)速度為8 km/h，凸起高度為6 mm，波峰數(shù)量為15。

3.3 對比驗證試驗

依照最優(yōu)方案參數(shù)組合，加工仿蚯蚓運動多功能開溝器，并與傳統(tǒng)的單元盤開溝器和雙圓盤開溝器進行對比試驗。試驗前，將試驗田平均分為9塊（A1、A2、A3…A9），每塊試驗田有6條壟，壟高300 mm、壟寬300 mm、每壟長為200 m，試驗時選用2BGD-6大豆耕播機先后安裝單元盤開溝器、雙圓盤開溝器和仿蚯蚓運動多功能開溝器后，單元盤開溝器在A1、A2和A3地塊進行播種作業(yè)，雙圓盤開溝器在A4、A5和A6地塊進行播種作業(yè)，仿蚯蚓運動多功能開溝器在A7、A8和A9地塊進行播種作業(yè)，并分別測定3種開溝器作業(yè)區(qū)域的平均土壤含水率、平均土壤溫度和平均出苗時間，試驗結(jié)果如圖13所示。

圖11 各因素對平均土壤溫度（0～100 mm）影響的響應(yīng)曲面Fig.11 Response surfaces influenced by various factors on average soil temperature(0-100 mm)

圖12 各因素對平均出苗時間影響的響應(yīng)曲面Fig.12 Response surfaces influenced by various factors on average emergence time

圖13 3種不同類型開溝器對大豆出苗狀況的影響Fig.13 Effect on emergence of soybean of three different kinds of openers

如圖13可知，仿蚯蚓運動多功能開溝器實際田間結(jié)果與計算機仿真結(jié)果基本一致，且相較于其他2種傳統(tǒng)開溝裝置，可顯著提高0～100 mm平均土壤含水率（P<0.01）和0～100 mm平均土壤溫度（P<0.01），并提前大豆平均出苗時間（P<0.01）。相較于2種傳統(tǒng)開溝裝置分別提高0～100 mm土壤平均含水率5.92%和4.86%，分別提高0～100 mm土壤平均溫度0.5和0.4℃，分別提前大豆平均出苗時間0.5和0.92 d。以上試驗結(jié)果表明：本文所設(shè)計仿蚯蚓運動多功能開溝器可有效疏松耕種層土壤，增加耕種層土壤孔隙度，進而有效提升地溫，雖然同時加速了表層土壤中的水分蒸發(fā)，但通過鎮(zhèn)壓力調(diào)節(jié)機構(gòu)，增大了對耕種層下方土壤的毛細(xì)現(xiàn)象，使下層土壤中蓄存的水資源得到有效利用，從而在總體上提升了耕種層土壤的含水率。

4 結(jié)論

1）本文所設(shè)計仿蚯蚓運動多功能開溝器設(shè)有仿蚯蚓動態(tài)松土開溝機構(gòu)和鎮(zhèn)壓力調(diào)節(jié)機構(gòu)，其可通過模擬蚯蚓的運動方式疏松土壤，通過溝底鎮(zhèn)壓作業(yè)壓實下層土壤，其關(guān)鍵結(jié)構(gòu)參數(shù)和作業(yè)速度均對土壤平均含水率（0～100 mm）、土壤平均溫度（0～100 mm）和大豆平均出苗時間具有顯著性影響，其中大豆平均出苗時間隨著凸起高度、波峰數(shù)量和作業(yè)速度的增大而先下降后上升，運用Design-Expert和Matlab軟件對試驗結(jié)果進行響應(yīng)曲面分析和回歸分析，并考慮到作業(yè)速度易受田間復(fù)雜地況的影響而較難保證其精準(zhǔn)性以及加工精度等問題，且規(guī)定波峰數(shù)量N只能為自然數(shù)，最終得出最優(yōu)參數(shù)組合為：作業(yè)速度8 km/h、松土機構(gòu)凸起高度6 mm和波紋圓盤波峰數(shù)量15。

2）本文所設(shè)計仿蚯蚓運動多功能開溝器相較于兩種傳統(tǒng)開溝裝置，可分別提高土壤平均含水率（0～100 mm）5.92%和4.86%，分別提高土壤平均溫度（0～100 mm）0.5和0.4℃，分別提前大豆出苗時間0.5和0.92 d，表明采用模擬蚯蚓生物學(xué)行為的松土機構(gòu)，可有效疏松耕種層土壤，增加耕種層土壤孔隙度，進而有效提升地溫，同時通過增設(shè)鎮(zhèn)壓力調(diào)節(jié)機構(gòu)對種溝下層土壤進行鎮(zhèn)壓作業(yè)，可增大對耕種層下方土壤中水分的吸收，進而提升耕種層土壤含水率，最終通過對耕種層土壤物理性狀的改善，達(dá)到促進大豆發(fā)育的目的。

[1]薛慶喜.中國及東北三省30年大豆種植面積、總產(chǎn)、單產(chǎn)變化分析[J].中國農(nóng)學(xué)通報，2013，29(35)：102－106.Xue Qingxi.Analysis on the change of 30 year's soybean areas,production and yield in China and Northeast China[J].Chinese Agricultural Science Bulletin,2013,29(35):102－106.(in Chinese with English abstract)

[2]董非非，劉愛民，封志明，等.大豆傳統(tǒng)產(chǎn)區(qū)種植結(jié)構(gòu)變化及影響因素的定量化評價:以黑龍江省嫩江縣為例[J].自然資源學(xué)報，2017，32(1)：40－49.Dong Feifei,Liu Aimin,Feng Zhiming,et al.Changes of planting structure and quantitative evaluation of influencing factors in traditionalsoybean producing areas:Taking nenjiang county in heilongjiang province as an example[J].Journal of Natural Resources,2017,32(1):40－49.(in Chinese with English abstract)

[3]孫仕軍，許志浩，張旭東，等.地膜覆蓋對玉米田間土壤含水率和地溫變化的影響[J].玉米科學(xué)，2015，23(3)：91－98.Sun Shijun,Xu Zhihao,Zhang Xudong,et al.Effect of film mulching on soil moisture and soil temperature in eastern inner mongolia rain-fed black soil area[J].Journal of Maize Sciences,2015,23(3):91－98.(in Chinese with English abstract)

[4]Xie Z K,Wang Y J,Li F M.Effect of plastic mulching on soil water use and spring wheat yield in arid region of northwest China[J].Agriculture Water Management,2005,75(1):71－83.

[5]Sarkar S,Paramanick M,Goswami S B.Soil temperature,water use and yield of yellow sarson(Brassica napusL.var.glauca)in relation to tillage intensity and mulch management under rained lowland ecosystem in eastern India[J].Soil＆Tillage Research,2007,93(1):94－101.

[6]員學(xué)峰，吳普特，汪有科.地膜覆蓋保墑灌溉的土壤水、熱以及作物效應(yīng)研究[J].灌溉排水學(xué)報，2006，25(1)：25－29.Yuan Xuefeng,Wu Pute,Wang Youke.Study on the effect of irrigation under plastic preservation of soil moisture on soil and crop[J].Journal of Irrigation and Drainage,2006,25(1):25－29.(in Chinese with English abstract)

[7]劉斐，陳軍，慕軍營，等.土壤溫度檢測及其與含水率關(guān)系研究[J].干旱地區(qū)農(nóng)業(yè)研究，2013，31(3)：95－99，117.Liu Fei,Chen Jun,Mu Junying,et al.Detection of soil temperature and its relationship with moisture content[J].Agricultural Research in the Arid Areas,2013,31(3):95－99,117.(in Chinese with English abstract)

[8]Eric H,Hamed P.Evaluation of several dielectric mixing models for estimating soi l moisture content in sand,loam and clay soils[J].ASAE,2003(7):26－30.

[9]伍海兵，李愛平，方海蘭，等.綠地土壤孔隙度檢測方法及其對土壤肥力評價的重要性[J].浙江農(nóng)林大學(xué)學(xué)報，2015，32(1)：98－103.Wu Haibing,Li Aiping,Fang Hailan,et al.Green-belt soil testing methods for porosity and the importance of porosity on soil fertility evaluation[J].Journal of Zhejiang A＆F University,2015,32(1):98－103.(in Chinese with English abstract)

[10]張學(xué)權(quán).不同植被恢復(fù)土壤容重和孔隙度特征分析[J].成都大學(xué)學(xué)報：自然科學(xué)版，2017，36(3)：325－327.Zhang Xuequan.Analysis of soil bulk density and porosity in differentvegetation restoration[J].JournalofChengdu University:Natural Science Edition,2017,36(3):325－327.(in Chinese with English abstract)

[11]Holen N M,War S M.Quantification of water storage in fingers associated with preferential flow in milled peat stockpiles[J].Soil Science Society of America Journal,1999,63(3):480－486.

[12]Skopp J,Gardner W R,Tyler E J.Solute movement in structured soils:Two-region model with small interaction[J].Soil Science Society of America Journal,1981,45(5):837－842.

[13]Haukka J.Spatial distribution and formation of earthworm burrows[J].Pedobiologia,1991,37:175－278.

[14]Asare S N,Rudra R P,Dickinson W T,et al.Soil macroporosity distribution and trends in a no-till plot using a volume computer tomography scanner[J]. Agricultural Engineering Research,2001,78(4).437－447.

[15]Taser O F,Kara O.Silage maize(Zea mays L.)Seedlings emergence as influenced by soil compaction treatments and contact pressures[J].Plant,Soil and Environment,2005,51(7):289－295.

[16]趙淑紅，劉宏俊，譚賀文，等.仿旗魚頭部曲線型開溝器設(shè)計與性能試驗[J].農(nóng)業(yè)工程學(xué)報，2017，33(5)：32－39.Zhao Shuhong,Liu Hongjun,Cao Hewen,et al.Design and performance experiment of opener based on bionic sailfish head curve[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2017,33(5):32－39.(in Chinese with English abstract)

[17]馬云海，馬圣勝，賈洪雷，等.仿生波紋形開溝器減黏降阻性能測試與分析[J].農(nóng)業(yè)工程學(xué)報，2014，30(5)：36－41.Ma Yunhai,Ma Shengsheng,Jia Honglei,et al.Measurement and analysis on reducing adhesion and resistance of bionic ripple opener[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2014,30(5):36－41.(in Chinese with English abstract)

[18]王文君.大豆精密播種機仿形仿生鎮(zhèn)壓裝置[D].長春：吉林大學(xué)，2016.Wang Wenjun.Bionic Press Device with Profiling Mechanism for Soybean Precision Planter[D].Changchun:Jilin University,2016.(in Chinese with English abstract)

[19]Lavelle P.Diversity of soil fauna and ecosystem function[J].Biol Int,1996,33:3－16.

[20]劉國敏.蚯蚓體表減粘降阻功能耦合仿生研究[D].長春：吉林大學(xué)，2009.Liu Guomin.Coupling Bionic Research on the Adhesion and Resistance Reduction of the Earthworm Surface[D].Changchun:Jilin University,2009.(in Chinese with English abstract)

[21]劉國敏，鄒猛，李建橋.蚯蚓體表自潤滑仿生應(yīng)用試驗研究[J].江西農(nóng)業(yè)大學(xué)學(xué)報，2010，32(2)：378－382.Liu Guomin,Zou Meng,LiJianqiao.Self-lubricating biomimetic surface earthworm experimental studies[J].Acta Agriculturae Universitatis Jiangxiensis,2010,32(2):378－382.(in Chinese with English abstract)

[22]南春容，張海智.蚯蚓運動的觀察[J].生物學(xué)通報，1998(11)：14.Nan Chunrong,Zhang Haizhi.Observation of earthworm movement[J].Bulletin of Biology,1998(11):14.(in Chinese with English abstract)

[23]安瑞永，黃薇.蚯蚓運動行為的觀察與研究[J].生物學(xué)通報，1995，30(6)：42.An Ruiyong,Huang Wei.Observation and study of earthworm movement behavior[J].Bulletin of Biology,1995,30(6):42.(in Chinese with English abstract)

[24]劉國敏，鄒猛，李建橋.蚯蚓體表與土壤接觸界面動態(tài)行為仿真[J].吉林大學(xué)學(xué)報：工學(xué)版，2010，40(6)：1609－1613.Liu Guomin,Zou Meng,Li Jianqiao.Interfacial dynamics simulation between soil and earthworm surface[J].Journal of Jilin University:Engineering and Technology Edition,2010,40(6):1609－1613.(in Chinese with English abstract)

[25]許亞婷.蚯蚓體表形貌及黏附特性研究[D].長春：吉林大學(xué)，2006.Xu Yating.Research on the Surface Feature of Earthworm and its Adhesion Characteristic[D].Changchun: Jilin University,2006.(in Chinese with English abstract)

[26]李安琪，任露泉，陳秉聰.蚯蚓體表液的組成及其減黏脫土機理分析[J].農(nóng)業(yè)工程學(xué)報，1990，6(3)：8－14.Li Anqi,Ren Luquan,Chen Bingcong.Earthworm surface liquid composition and mechanism analysis visbreaking bentonite[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),1990,6(3):8－14.(in Chinese with English abstract)

[27]許亞婷，李建橋，孫久榮.蚯蚓體表形貌及潤濕性的研究[J].華中農(nóng)業(yè)大學(xué)學(xué)報，2005，6(1)：40－43.Xu Yating,Li Jianqiao,Sun Jiurong.Earthworm surface morphology and wettability study[J].Journal of Huazhong Agricultural University,2005，6(1)：40－43.(in Chinese with English abstract)

[28]任露泉，王云鵬，李建橋，等.典型生物柔性非光滑體表的防粘研究[J].農(nóng)業(yè)工程學(xué)報，1996，12(4)：31－36.Ren Luquan,Wang Yunpeng,Li Jianqiao,et al.Bionic research on flexible nonsmooth surface of typical animals[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),1996,12(4):31－36.(in Chinese with English abstract)

[29]Celik A,Ahmet S,Altikat T R,et al.Strip tillage width effects on sunflower seed emergence and yield[J].Soil&Tillage Research,2013,131(7):20－27.

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

[31]孫桓，陳作模，葛文杰，等.機械原理[M].北京：高等教育出版社，2006：238－239.