卜祥利，廖慶喜，孫文成，魏國粱，張青松，王鵬程

華中農(nóng)業(yè)大學(xué)工學(xué)院/農(nóng)業(yè)農(nóng)村部長江中下游農(nóng)業(yè)裝備重點實驗室，武漢 430070

長江中下游是我國冬油菜主產(chǎn)區(qū)，該區(qū)域降雨量充足[1]。油菜根系忌水，種植油菜時需要開溝排水，開好“三溝”(廂溝、腰溝、畦溝)且要求三溝相通是避免油菜漬害的保證[2]。目前開溝方式主要分為人工開溝和機(jī)械開溝2種。傳統(tǒng)油菜種植，人工開溝勞動強(qiáng)度大，開溝效率低，因此在大田種植時，多用機(jī)械開溝作業(yè)。機(jī)械開溝多采用被動式犁體或主動式開溝機(jī)具。主動式開畦溝裝置作業(yè)效果好，但功耗高，而被動式犁體結(jié)構(gòu)簡單，能適應(yīng)較高的作業(yè)速度，在農(nóng)田開溝領(lǐng)域得到了廣泛應(yīng)用[3-5]。

國內(nèi)外學(xué)者對農(nóng)田開溝技術(shù)進(jìn)行了研究[6-10]，但主要針對大田或旱田的開溝機(jī)具，適用于長江中下游油菜種植的開溝器研究相對較少。目前，華中農(nóng)業(yè)大學(xué)工學(xué)院研制的2BFQ-6型油菜聯(lián)合直播機(jī)采用前后犁組合式開畦溝，前犁破土、后犁整形，可形成完整畦溝，在油菜種植時被廣泛應(yīng)用[11]。張青松等[12-13]利用EDEM仿真對犁體曲面進(jìn)行了研究，提出曲面優(yōu)化方法實現(xiàn)犁體降附減阻。但鏵式犁在土壤黏重板結(jié)、含水率高、稻茬留量大的工況下作業(yè)時，開溝后犁易產(chǎn)生粘土、壅土、掛草或堵草現(xiàn)象，導(dǎo)致作業(yè)后溝邊土壤堆積、溝壁垮塌等問題，造成溝型不穩(wěn)定。為提高開溝穩(wěn)定性，包攀峰等[14]研制了雙翼開溝犁與后置分土螺旋器組合完成開溝作業(yè)的分土開溝裝置，為主動與被動開溝裝置相結(jié)合的研究提供了參考。劉曉鵬等[15]設(shè)計了一種配合鏵式前犁完成開畦溝功能的組合式船型開溝器，提高了油菜聯(lián)合直播機(jī)開畦溝的穩(wěn)定性，并通過微分幾何分析與EDEM仿真對犁體曲面進(jìn)行了優(yōu)化，但需與破土前犁配套使用。

本研究在2BFQ-8油菜聯(lián)合直播機(jī)用組合式船型開溝器的研究基礎(chǔ)上，針對長江中下游土壤黏重板結(jié)、含水率波動大、前茬水稻秸稈留存量大的復(fù)雜工況和油菜種植需開畦溝的生產(chǎn)要求，考慮到銳角開溝犁具有開溝阻力小、土壤擾動小的特點，研制了一種與驅(qū)動型犁旋聯(lián)合耕整機(jī)配套使用的仿靴形銳角開溝犁，旨在為長江中下游油菜開畦溝裝置的研究提供參考。

1 材料與方法

1.1 聯(lián)合耕整機(jī)結(jié)構(gòu)與工作原理

配備仿靴形銳角開溝犁的驅(qū)動型犁旋聯(lián)合耕整機(jī)主要由主機(jī)架、驅(qū)動圓盤犁組、旋耕裝置、人字形齒輪傳動箱、平土托板、仿靴形銳角開溝犁等組成，整機(jī)結(jié)構(gòu)如圖1A所示。仿靴形銳角開溝犁安裝于機(jī)組中部，對土壤進(jìn)行側(cè)向擠壓和溝底托平，形成完整溝型。通過犁耕、旋耕相匹配的聯(lián)合作業(yè)，一次作業(yè)完成犁耕、旋耕、秸稈切碎還田、埋茬、開畦溝等多道工序，形成適宜油菜生長的種床。仿靴形銳角開溝犁如圖1B所示，主要由犁尖、滑切刃、犁柱、支撐梁、擋土板、仿靴整形板、擋土板、底板、起土板等組成。犁尖位于開溝犁最前端，犁尖與刃口曲線隨牽引裝置前進(jìn)時破土，隨后由分土板將土壤向兩側(cè)分離，引導(dǎo)土壤至仿靴整形板和擋土板處，并將土壤提升至地面，仿靴整形板對側(cè)邊土壤進(jìn)行擠壓成型，擋土板在對土壤擠壓的同時防止土壤回流，從而形成完整的溝側(cè)壁，仿靴整形板與擋土板成一定角度，保證梯形斜邊的形成。支撐梁起支撐作用，防止擋土板在工作過程中受土壤作用擠壓變形，影響開溝質(zhì)量。底部安裝有底板，底板在拖動作用下對溝底土壤作用形成平整溝底。犁柱上有銷孔，便于安裝和調(diào)節(jié)開溝深度。

A:整機(jī)結(jié)構(gòu)圖 Structure diagram of machine; B:仿靴形銳角開溝犁示意圖 Shoe-shaped acute angle groove plough; 1.三點懸掛 Three-point suspension; 2.機(jī)架 Main frame; 3.驅(qū)動圓盤犁組 Drive disc plows; 4.平土托板 Flat soil pallet； 5.旋耕裝置 Rotary cultivation mechanism; 6.仿靴形銳角開溝犁 Shoe-shaped acute angle opener; 7.人字形傳動齒輪箱 Herringbone transmission gearbox; 8.犁尖 Plow tip; 9.滑切刃Sliding cutting edge; 10.犁柱 Plow column; 11.支撐梁 Support beam;12.擋土板 Retaining plate; 13.仿靴整形板 Boot-shaping plate; 14.底板 Base plate; 15.起土板 Lifting plate.圖1 驅(qū)動型犁旋聯(lián)合耕整機(jī)結(jié)構(gòu)示意圖Fig.1 Structure of the driven plough-rotating combine machine

1.2 仿靴形銳角開溝犁主要參數(shù)設(shè)計

1)滑切刃口曲線設(shè)計。仿靴形銳角開溝犁采用對稱式結(jié)構(gòu)設(shè)計，由前端兩破土曲面相交形成滑切刃口曲線，對仿靴形銳角開溝犁破土阻力影響較大。開溝犁入土后隨機(jī)組前進(jìn)，土壤受刃口曲線作用發(fā)生剪切破壞，參考文獻(xiàn)[6,16]，根據(jù)滑切原理，摩擦角應(yīng)小于滑切角。本研究選用拋物線函數(shù)曲線作為滑切刃口曲線，滑切刃口曲線的形狀如圖2所示。

圖2 滑切刃口曲線示意圖Fig.2 Schematic diagram of smooth cutting edge curve

根據(jù)拋物線方程，設(shè)滑切刃口曲線方程為：

y=ax2+b

(1)

式(1)中，b為常數(shù)，函數(shù)位置隨b的大小而變化，且不會發(fā)生線形變化，故令常數(shù)b為0。設(shè)滑切刃口曲線起始點A(x1，y1)、終止點B(x2，y2)，則A、B兩點的垂直距離為Hq，水平距離為Lq，A、B兩點的斜率為：

(2)

式(2)中，γ1為曲線起始點A的切線與x軸夾角，(°)；γ2為曲線終止點B的切線與x軸夾角，(°)；聯(lián)立式(1)、(2)可得常數(shù)a的關(guān)系式：

(3)

由于滑切角α1、α2與γ1、γ2互為余角，可得最終滑切刃口曲線方程為:

(4)

式(4)中，α1為起始滑切角，(°)；α2為終止滑切角，(°)。如式(4)所示，滑切刃口曲線的形狀受滑切角和A、B兩點垂直距離大小的影響。由圖2可知，滑切角成遞增趨勢，即α2>α1。根據(jù)文獻(xiàn)[17]，取摩擦角為23°，即α1>23°。α1與α2互為余角，所以α2<67°，因此終止滑切角的取值范圍為23°<α2<67°。當(dāng)曲線高度Hq不變時，改變滑切角大小將會得到不同的滑切刃口曲線。

2)犁體入土角及犁尖張角的設(shè)計分析。犁體能否較好入土和分土與入土角和犁尖張角大小密不可分。犁尖在入土的同時，將土壤從溝底提升到地面并將土壤導(dǎo)向兩側(cè)，使犁體曲面能夠發(fā)揮作用。參考犁體入土角的設(shè)計方法，本研究選用三面楔模型對犁尖進(jìn)行分析(圖3)。

圖3 三角楔模型分析Fig.3 Analysis of triangle wedge model

以機(jī)具前進(jìn)方向的反方向為X軸正方向，溝寬方向為Y軸，垂直方向為Z軸，建立空間坐標(biāo)系。其中，面OAC為犁尖左半部分，面ABC為犁尖作業(yè)橫截面，γ為犁尖張角的1/2，δ為入土角，λ為翻土角。機(jī)具以一定速度前進(jìn)時，假設(shè)土壤顆粒以最快速度穿過犁尖進(jìn)入仿靴整形板區(qū)域，則其運動軌跡OE應(yīng)垂直于AC，過F點作與AB成η角的直線EF，則η角為土壤堆積角。由圖3可知：

(5)

由式(5)可得，tanλ=tanδ/tanγ=45°，所以δ=γ。根據(jù)農(nóng)機(jī)設(shè)計手冊，當(dāng)入土角過大時，土層抬高，工作阻力增加；入土角過小時，會增加犁體高度及長度造成浪費。犁尖張角偏大時，土壤易被推向兩側(cè)外翻而亂，同時容易纏繞秸稈；犁尖過小時，切斷草根的能力減弱[18]。綜合考慮設(shè)定入土角為30°、犁尖張角為60°，即δ=γ=30°，即入土角和犁尖張角均為銳角。

3)擋土板與仿靴整形板的設(shè)計分析。擋土板與仿靴整形板采用一體式設(shè)計，兩者共同組成的犁體曲面如圖4所示。其中，面BCDE為擋土曲面，面ABEF為仿靴整形曲面，面OKCA為起土曲面，起土曲面參考農(nóng)機(jī)設(shè)計手冊的水平直元線法形成，根據(jù)滑切刃口曲線高度及犁體高度要求，起土曲面高度為300 mm，曲面寬度為200 mm。

圖4 犁體曲面示意圖Fig.4 Schematic diagram of plough surface

仿靴形銳角開溝犁對土壤作用后，土壤顆粒在梯形溝溝壁上保持穩(wěn)定。假設(shè)開溝深度一定，且機(jī)組勻速前進(jìn)，則擋土板與仿靴形整形板跟土壤接觸的每個點上沿法向方向擠壓土粒，以斜面上單個土壤顆粒為研究對象，將其看作質(zhì)點M，土壤間的相互作用忽略不計，質(zhì)點M的受力情況如圖5所示，土壤顆粒受重力G、溝壁對其產(chǎn)生的支持力FN、摩擦力Ff。土壤回流時，土壤顆粒沿溝壁邊下滑，則沿Y方向的加速度ay為0。

圖5 土壤顆粒受力分析Fig.5 Force analysis of soil slope

則土壤顆粒的平衡方程為：

(6)

式(6)中，m為土壤顆粒的質(zhì)量，g為重力加速度，μ為土壤顆粒間動摩擦系數(shù)，θ為仿靴整形板與擋土板的傾角(簡稱翼板傾角)。由式(6)可知，若μ=tanθ，即土壤的動摩擦系數(shù)等于土壤溝壁堆積角正切值時，土壤顆粒不流入溝底；若μ

ax=gsinθ-μgcosθ

(7)

根據(jù)土壤顆粒在X軸方向的運動方程，可知土壤顆粒沿X軸方向的位移Sx，即土壤回流位移為：

(8)

(9)

式(9)中，SH為土壤回流位移，mm。由于開溝犁以均勻速度和穩(wěn)定深度作業(yè)，則任意一段區(qū)域內(nèi)土壤回流的時間t0不變，視為常數(shù)，則根據(jù)式(9)可知，土壤回流量與擋土板和仿靴整形曲面間的夾角θ有關(guān)。

設(shè)左右仿靴整形板的底端距離和擋土板的上端距離分別為l1、l2，考慮土壤流動性，為保證形成梯形溝，仿靴形整形板的下端距離設(shè)計為200 mm。若擋土板上寬度過大，導(dǎo)致傾角θ過大，會將過多擾動土推出溝外。考慮土壤的擾動和流動特性，取l2=350 mm。為保證開溝器作業(yè)后畦溝無土壤回流，取擋土板長度l為1.5(h+h′)，其中h為開畦溝深度，h′為擾動土壟起高度，取20 mm。根據(jù)《合理耕層構(gòu)建的水旱輪作耕整作業(yè)技術(shù)規(guī)程要求》，開溝深度為150～250 mm，所以擋土板的長度為255 mm

(10)

擋土板與仿靴整形曲面的夾角可由式(10)計算得出：0<θ<79.2°。根據(jù)試驗可得長江中下游地區(qū)土壤含水率在22.86%～32.46%時的土壤堆積角為33.94°～39.76°，為防止溝壁崩塌、溝內(nèi)壅土等現(xiàn)象出現(xiàn)，設(shè)計傾角θ應(yīng)大于40°，所以40°<θ<79.2°。

1.3 田間開溝質(zhì)量試驗

1)試驗器材與條件。試驗器材包括：卷尺(5.0 m)、直尺(500 mm)、土壤堅實度測量儀、取土鉆、烘干箱、500 mm×500 mm自制鐵框等。試驗田塊常年稻-油輪作，土壤類型屬于粉質(zhì)黏壤土。經(jīng)測量，試驗環(huán)境特性參數(shù)如下：土壤深度0～250 mm的平均堅實度為1 173.08 kPa，平均土壤含水率為23.9%，土壤容重為1.12 g/cm3，秸稈覆蓋量為1 256.4 g/m2，秸稈平均留茬高度為330 mm，秸稈平均含水率為38.6%。

2)試驗方法。以仿靴形銳角開溝犁為試驗對象，仿靴形銳角開溝犁與驅(qū)動型犁旋耕整機(jī)配套使用，配套動力均為東風(fēng)井關(guān)T954拖拉機(jī)，參考NY/T 740―2003《田間開溝機(jī)械作業(yè)質(zhì)量》，通過調(diào)節(jié)拖拉機(jī)擋位和手油門保證機(jī)組前進(jìn)速度為2.6 km/h，機(jī)組沿直線行駛50 m，取中間行程30 m為穩(wěn)定測量區(qū)域，在測量區(qū)域等距取10個測量點，測繪每個測量點處的廂面斷面形狀，每組試驗重復(fù)3次。試驗指標(biāo)為溝寬、溝寬穩(wěn)定性、溝深、溝深穩(wěn)定性、土壤黏附量等。

①溝深及其穩(wěn)定性。以繪制的溝型斷面圖為基礎(chǔ)，以耕作平面為基準(zhǔn)，測量溝底與耕作平面的差值，取10個測量點的數(shù)據(jù)平均值作為作業(yè)溝深。根據(jù)計算得出每組溝型的溝深值，溝深穩(wěn)定性系數(shù)按式(11)～(12)計算。

(11)

(12)

式(12)中，a為溝深，mm；ai表示所測行程中第i個測試點的溝深，mm；n表示所測行程的測試點數(shù)；Ua為溝深穩(wěn)定性系數(shù)，%。

②溝寬及其穩(wěn)定性。以繪制的溝型斷面圖為基礎(chǔ)，任意選取溝型斷面一側(cè)的測量點坐標(biāo)為計算基準(zhǔn)，以選取點的縱坐標(biāo)值為基準(zhǔn)在溝型斷面圖中繪制溝寬基準(zhǔn)直線，基準(zhǔn)線與左右溝壁相交所得點的橫坐標(biāo)差值(精確到1 mm)作為畦溝的寬度。依次計算10個測量點的數(shù)據(jù)取平均值即為作業(yè)溝寬。溝寬穩(wěn)定性的計算方法與溝深穩(wěn)定性的計算方法相同。

③土壤黏附量。機(jī)組每行程作業(yè)完成后，將黏附在犁體的土壤秸稈混合物刮下裝入自封袋中稱量，每組試驗重復(fù)3次取平均值，即為該犁體的土壤黏附量。

1.4 田間開溝效果正交試驗

影響開溝犁開溝效果的因素主要有開溝深度、開翼板傾角、作業(yè)速度、土壤工況等[19]。在不考慮土壤工況的前提下，即試驗在同一地塊進(jìn)行，試驗選擇翼板傾角X1、開溝深度X2、作業(yè)速度X3為主要考察因素。本研究以回土量作為評價指標(biāo)。回土量的多少直接反映了畦溝的塑形好壞以及溝壁土壤是否均勻緊實。為更加直觀準(zhǔn)確地測量評價指標(biāo)，使用回土體積與理論溝型體積比作為回土量的計算方法[6]。

在田間試驗過程中，通過調(diào)節(jié)拖拉機(jī)懸掛桿控制銳角開溝犁開溝深度保持一定，調(diào)節(jié)拖拉機(jī)擋位和手油門控制機(jī)組以勻速前進(jìn)，忽略由于地表不平或拖拉機(jī)輪胎滑移造成耕深和前進(jìn)速度波動的影響。沿畦溝中段取30 m的范圍，隨機(jī)取10個橫斷面，實際畦溝截面示意圖如圖6所示。測量實際開溝深度h1、實際上溝寬度d1、下溝寬度測量值d2、去除回土后開溝深度H1、去除回土后實際下溝寬度，約看作為d3=d2，則回土量可通過式(13)計算得出。

H1為實際開溝深度,mm;h1為開溝深度測量值,mm;d1為上溝寬度,mm;d2為下溝寬度測量值，mm;d3為實際下溝寬度,mm。H1 is the actual ditch depth,mm; h1 is the measured ditch depth,mm; d1 is the upper ditch width,mm; d2 is the lower ditch width measurement,mm;d3 is the actual ditch width,mm.圖6 實際畦溝截面示意圖Fig.6 Schematic diagram of actual furrow section

(13)

式(13)中，φ為回土量，%；V′為理論溝形體積，mm3；V為回土體積，mm3。

2 結(jié)果與分析

在湖北省荊州市監(jiān)利縣水稻-油菜(輪作)全程機(jī)械化生產(chǎn)示范基地開展田間試驗，田間試驗分為開溝質(zhì)量試驗和正交試驗。

2.1 田間開溝質(zhì)量

依據(jù)試驗方法，仿靴形銳角開溝犁開出畦溝效果見圖7。試驗結(jié)果表明，仿靴形銳角開溝犁開畦溝溝寬為200～400 mm，溝寬穩(wěn)定性為90.6%，溝深為175～250 mm，溝深穩(wěn)定性為92.4%，畦溝斷面形狀為梯形，作業(yè)穩(wěn)定性高；且機(jī)組作業(yè)完成后，犁體土壤黏附量平均為238 g，犁體上基本無掛草壅土，說明機(jī)組通過性好。各項試驗指標(biāo)均滿足油菜種植農(nóng)藝要求。

圖7 開畦溝作業(yè)效果Fig.7 Effect diagram of trenching experiment

2.2 田間開溝效果

經(jīng)過前期計算設(shè)計，翼板傾角為40°～79.2°，開溝深度為170～250 mm，驅(qū)動型犁旋耕整機(jī)的作業(yè)速度為0.58～1.0 m/s，即仿靴形銳角開溝犁作業(yè)速度為0.58～1.0 m/s。在試驗可控范圍內(nèi)，對翼板傾角X1、開溝深度X2、作業(yè)速度X3分別選擇3個水平，采用三因素三水平Box-Behnken響應(yīng)曲面試驗法進(jìn)行優(yōu)化試驗，采用多元二次回歸方程擬合因素與響應(yīng)值間的函數(shù)關(guān)系，最終對回歸方程分析以尋求最優(yōu)參數(shù)。共17組試驗，因素水平如表1所示。

表1 試驗因素及水平Table 1 Factors and levels

利用Design-Expert軟件對數(shù)據(jù)進(jìn)行處理，得出試驗方案與結(jié)果如表2所示，回土量的方差分析如表3所示。

表2 試驗方案與結(jié)果Table 2 Test scheme and result of experiment

表3 回土量方差分析Table 3 Variance analysis of the quantity of backfill

對表3中的數(shù)據(jù)進(jìn)行多元回歸擬合，可得回土量φ的二次回歸方程為：

由表3分析可知，建立二次回歸模型的P<0.01，表明回歸模型極顯著；失擬項P>0.05，即失擬不顯著，說明模型所擬合的二次回歸方程與實際相符，能正確反映回土量φ與X1、X2、X3的關(guān)系，回歸模型可較好地對試驗中各試驗結(jié)果進(jìn)行預(yù)測。各影響因子對回土量的顯著性順序從小到大依次為：翼板傾角<開溝深度<作業(yè)速度。

2.3 參數(shù)優(yōu)化與驗證

根據(jù)油菜種植農(nóng)業(yè)要求，以開溝回土量最小為優(yōu)化目標(biāo)，在翼板傾角為45°～65°、開溝深度為170～250 mm、作業(yè)速度為0.6～1.0 m/s約束條件下進(jìn)行優(yōu)化，用Design-Expert軟件計算得出最佳參數(shù)組合為：翼板傾角為46.22°、開溝深度為183.06 mm、作業(yè)速度為0.608 m/s時，回土量為2.34%。為驗證優(yōu)化結(jié)果的可行性，按照實際翼板傾角為46°、開溝深度為183 mm、作業(yè)速度為0.6 m/s的參數(shù)組合在田間進(jìn)行3次重復(fù)試驗，驗證結(jié)果顯示開溝回土量為3.28%，回歸方程計算預(yù)測值為2.34%，試驗結(jié)果相差0.94個百分點，驗證了模型的可靠性。

3 討 論

本研究設(shè)計了油菜種床整備開畦溝仿靴形銳角開溝犁。根據(jù)犁體設(shè)計原則，對開溝犁的滑切刃口曲線、犁體入土角及犁尖張角進(jìn)行了設(shè)計分析，選用拋物線作為滑切刃口曲線，通過分析確定了切削刃口的滑切角范圍為23°～67°，并求得刃口曲線方程；通過建立三面楔模型對犁尖進(jìn)行了分析，確定了犁體入土角為30°、犁尖張角為60°，即犁體入土角與犁尖張角均為銳角；建立了仿靴形銳角開溝犁與土壤相互作用的動力學(xué)模型，確定了擋土板的長度為400 mm，擋土板與仿靴整形板的夾角為40°～79.2°。

本研究開展了仿靴形銳角開溝犁的開溝質(zhì)量試驗與正交試驗。開溝質(zhì)量試驗表明：仿靴形銳角開溝器可開出溝寬為200～400 mm、溝深為175～250 mm的梯形溝，溝寬及溝深穩(wěn)定性均達(dá)90%以上。正交試驗表明：翼板傾角、開溝深度、作業(yè)速度均對開溝質(zhì)量有影響，且翼板角度為46°、開溝深度為183 mm、作業(yè)速度為0.6 m/s時，開溝作業(yè)效果滿足油菜種植的農(nóng)藝要求。

本研究僅開展了仿靴形銳角開溝犁在同一工況下的作業(yè)質(zhì)量試驗，仿靴形銳角開溝犁的功耗及不同區(qū)域土壤環(huán)境的適用性還有待于進(jìn)一步研究。