馬 帥，徐麗明，袁全春，牛 叢，王爍爍，袁訓騰

（中國農(nóng)業(yè)大學工學院，北京 100083）

0 引 言

中國葡萄種植區(qū)域廣泛，是歐洲以外種植面積最大的國家[1-2]，其中大部分優(yōu)質(zhì)葡萄產(chǎn)區(qū)位于北方，為露天籬架式栽植，由于北方地區(qū)冬季寒冷干燥，故每當冬季來臨之前，需將剪枝后的葡萄藤下架并埋入土中，待第二年春季天氣回暖時，再將其出土上架[3]。由于葡萄藤出土作業(yè)季節(jié)性強、工作量大、勞動力成本高[4-5]，故采用機械化作業(yè)是必然發(fā)展趨勢。

近年來，國內(nèi)相關(guān)領(lǐng)域?qū)＜遗c學者對葡萄藤防寒土清土機的設(shè)計與應(yīng)用進行了研究，如曾保寧等[6-7]研究了以刮板為主要清土部件的防寒土清土機，該機結(jié)構(gòu)簡單、制造成本低、適用范圍廣、作業(yè)效率高，但清土不徹底，需要人工進行二次作業(yè)；劉芳建等[8-9]研究了以旋轉(zhuǎn)攪龍為主要清土部件的防寒土清土機，作業(yè)時，機器通過攪龍旋轉(zhuǎn)傳送土壤，出土量大，且能將清除的葡萄藤防寒土均勻地平鋪到行間，但攪龍部件為剛性體，在作業(yè)過程中容易損傷土壤中的葡萄藤；王文斌[10]研制了一種葡萄藤防寒土清除開溝機，通過觸桿識別水泥柱的位置，采用螺旋旋耕的方式進行破土和側(cè)向輸送土壤，具有一定的清土效果，但葡萄藤附近的防寒土不能清除，并且整機采用四連桿機構(gòu)，移動部件總質(zhì)量較大，慣性大，作業(yè)操作難度大；謝冬等[11-12]研制的葡萄藤防寒土清土機也是以旋轉(zhuǎn)攪龍為清土部件，但僅應(yīng)用于葡萄藤防寒土覆蓋之前先鋪設(shè)彩條布的情況，通過回收彩條布進而將防寒土清除干凈；馬帥等[13]研制的葉輪旋轉(zhuǎn)式葡萄藤防寒土清土機，采用刮土板、葉輪和旋轉(zhuǎn)攪龍組合的方式對葡萄藤防寒土進行清除，但僅適用于防寒土埋在水泥柱單側(cè)的情況，若兩側(cè)埋土，水泥柱位于土壟中間，則水泥柱中間的土壤不能被有效清除。國外葡萄產(chǎn)區(qū)主要集中在溫帶及亞熱帶地區(qū)，冬季氣候溫和濕潤，葡萄生產(chǎn)作業(yè)沒有冬季埋土防寒和春季清土上架環(huán)節(jié)，故未見有關(guān)葡萄藤防寒土清土機方面的研究[3,5]。

目前關(guān)于葡萄藤防寒土清土機的研究大多以剛性體清土部件為主，作業(yè)時不僅會損傷葡萄藤，而且清土不徹底，此外關(guān)于清土機的整體研究理論較少、研究方法單一。近些年來，隨著計算機技術(shù)的發(fā)展，離散元法（Discrete Element Method，DEM）及其仿真軟件EDEM 在農(nóng)業(yè)機械設(shè)計中得到了廣泛使用，越來越多的學者將其應(yīng)用于土壤離散元仿真參數(shù)的標定[14-16]和觸土作業(yè)機具的設(shè)計及優(yōu)化[17-19]，其中觸土部件的研究目前主要集中在開溝施肥機、深松鏟和排種器等方面。葡萄藤防寒土清土機作為一種典型的觸土作業(yè)機具，很少有人借助于該方法進行相關(guān)研究，此外，為清除水泥柱間的防寒土，清土機在作業(yè)時不僅與土壤間相互接觸作用，而且需要避開水泥柱，為綜合模擬整機作業(yè)效果，研究機具避障機構(gòu)運動的同時直觀分析觸土部件的清土效果，需利用多體動力學仿真軟件與離散元仿真軟件耦合，目前關(guān)于此耦合仿真試驗方法的應(yīng)用也鮮有研究。

綜上所述，為解決葡萄藤防寒土清土機自動化程度低、清土不徹底和易傷藤等問題，本文基于多體動力學軟件RecurDyn 和離散元軟件EDEM 耦合仿真的方法，針對中國北方砂壤土地區(qū)葡萄園內(nèi)的葡萄藤防寒土清除，提出采用柔性刷子與橡膠組合的清土部件，設(shè)計了一種自動避障式葡萄藤防寒土清土機，并對其工作原理進行分析，確定其關(guān)鍵部件的結(jié)構(gòu)及參數(shù)；在耦合仿真試驗的基礎(chǔ)上，開展四因素二次回歸正交旋轉(zhuǎn)中心組合模擬試驗，優(yōu)化作業(yè)參數(shù)組合并通過田間試驗驗證，以期在解決葡萄藤防寒土清除問題的同時，為觸土機具避障機構(gòu)的整體設(shè)計及優(yōu)化提供參考。

1 自動避障式葡萄藤防寒土清土機設(shè)計

1.1 整機結(jié)構(gòu)

自動避障式葡萄藤防寒土清土機主要由機架、避障擺動機構(gòu)、避障信號采集機構(gòu)、控制器、清土部件、擋土板、傳動部件和限深輪等組成，整機結(jié)構(gòu)如圖1 所示。避障擺動機構(gòu)主要包括擺動支架和避障油缸，擺動支架通過中間旋轉(zhuǎn)軸以鉸接方式安裝在機架的右側(cè)，并將清土部件安裝在擺動支架的最右端，清土部件采用柔性刷子和橡膠組合的結(jié)構(gòu)形式；避障信號采集機構(gòu)包括觸桿、避障位移傳感器、氣彈簧和同步位移傳感器等，同步位移傳感器安裝在避障油缸的底部并與其保持平行，以保證同步位移傳感器和避障油缸的伸縮運動一致，觸桿安裝在擺動支架的前端，用于感知前方障礙物（水泥柱），當觸桿轉(zhuǎn)動的角度達到控制器中預先設(shè)定的閾值時，控制器通過控制電磁換向閥使避障油缸縮回，進而帶動清土部件避開水泥柱；擋土板安裝在機架底部，可以在機架底部的前后橫梁上左右移動，用來防止被清除的土壤過度向清土部件在機架上安裝位置的另一側(cè)拋灑；傳動部件采用2 級帶傳動；2 個限深輪分別安裝在機架的后端兩側(cè)，通過調(diào)節(jié)限深輪，改變機具在作業(yè)過程中的離地高度，增強機具作業(yè)時的穩(wěn)定性。

1.2 自動避障系統(tǒng)

葡萄藤防寒土清土機自動避障系統(tǒng)主要由電源部分、避障感知部分、控制部分和液壓執(zhí)行部分組成，如圖2 所示。避障信號采集機構(gòu)中的觸桿為自動避障系統(tǒng)的感知部分，用于識別障礙物的存在；STM32 單片機作為控制部分中的主控制器，用于采集避障位移傳感器和同步位移傳感器的信號，然后輸出信號控制電磁換向閥，進而控制避障油缸的伸縮；液壓部分主要由液壓閥組和避障油缸組成，其中液壓閥組包括溢流閥、調(diào)速閥和電磁換向閥，溢流閥用于穩(wěn)定油路壓力，調(diào)速閥用于調(diào)節(jié)避障油缸的伸縮速度，電磁換向閥采用三位四通且中位機能為M 型結(jié)構(gòu)，用于控制避障油缸的伸縮，避障油缸作為執(zhí)行部件，用于帶動清土部件避開障礙物。工作時，液壓油由拖拉機后側(cè)液壓輸出口提供，當觸桿碰到水泥柱時，觸桿繞自身旋轉(zhuǎn)軸轉(zhuǎn)動，當轉(zhuǎn)動角度達到STM32單片機程序中預先設(shè)定的角度閾值時，STM32 單片機輸出信號給繼電器以控制電磁換向閥閥芯的移動，進而控制避障油缸縮回，從而帶動清土部件避開障礙物。避開障礙物后，觸桿在氣彈簧的作用下逐漸恢復到初始位置，當觸桿轉(zhuǎn)動的角度恢復到小于STM32 單片機程序中設(shè)定的角度閾值時，STM32 單片機再次輸出信號給繼電器以控制電磁換向閥閥芯的移動，進而控制避障油缸伸出，從而帶動清土部件伸入水泥柱間繼續(xù)進行清土作業(yè)，依此反復。在避障油缸伸縮的同時，與其平行的同步位移傳感器也隨之伸縮，當STM32 單片機通過同步位移傳感器檢測到避障油缸完全縮回或完全伸出時，輸出控制信號使電磁換向閥的閥芯位于中位，由于采用的是M 型中位機能的電磁換向閥，故由液壓源輸出的液壓油直接經(jīng)相通的進、回油口直接回到油箱，從而使液壓油達到卸荷的目的。

圖1 自動避障式葡萄藤防寒土清土機結(jié)構(gòu)圖 Fig.1 Schematic diagram of automatic obstacle-avoiding grapevine cold-proof soil cleaners

1.3 工作原理

自動避障式葡萄藤防寒土清土機作業(yè)時，由拖拉機牽引前進，同時通過動力輸出裝置（Power Take-Off，PTO）將動力傳給變速箱，由變速箱輸出軸將動力經(jīng)二級帶傳動傳給清土部件，從而實現(xiàn)其旋轉(zhuǎn)清土的功能。在機具前進過程中遇到水泥柱時，觸桿與其碰觸，當系統(tǒng)檢測到避障信號后，避障油缸縮回進而帶動清土部件避開水泥柱，在避障油缸縮回的同時，同步位移傳感器也隨之縮回，當檢測到避障油缸完全縮回時，電磁換向閥不通電，液壓油卸荷；當避開水泥柱后，觸桿和避障位移傳感器在氣彈簧的作用下恢復到初始位置，避障油缸帶動清土部件伸出，同步位移傳感器也隨之伸出，清土部件重新進入水泥柱間進行清土作業(yè)，當檢測到避障油缸完全伸出時，電磁換向閥再次斷電，液壓油卸荷，以此往復，完成葡萄藤防寒土的避障清土作業(yè)過程。

圖2 自動避障系統(tǒng) Fig.2 Automatic obstacle-avoiding system

1.4 整機主要技術(shù)參數(shù)

結(jié)合中國北方地區(qū)葡萄種植的行距和葡萄藤春季清土情況，確定自動避障式葡萄藤防寒土清土機的主要技術(shù)參數(shù)，如表1 所示。

2 關(guān)鍵部件及其主要技術(shù)參數(shù)

2.1 避障信號采集機構(gòu)

避障信號采集機構(gòu)主要由安裝支架、觸桿、避障位移傳感器、氣彈簧、限位裝置和同步位移傳感器等組成，如圖3 所示（圖中未顯示同步位移傳感器）。設(shè)觸桿未碰到水泥柱時的位置為P1，碰到水泥柱后，觸桿繞O 點旋轉(zhuǎn)，旋轉(zhuǎn)到P2位置時，此時觸桿末端位于清土部件外側(cè)邊緣位置，同時觸桿另一端由a1點運動到a2點，并帶動避障位移傳感器和氣彈簧伸出。當觸桿轉(zhuǎn)動的角度大于STM32 單片機程序中預先設(shè)定的角度閾值時，避障油缸縮回進而帶動清土部件向后擺動避開水泥柱。為了防止觸桿在作業(yè)過程中由于旋轉(zhuǎn)角度過大而損壞傳感器，故在安裝支架的右側(cè)設(shè)置了限位裝置，主要包括前限位擋塊和后限位擋塊，以保證觸桿在一定范圍內(nèi)轉(zhuǎn)動。

觸桿由一段直線部分和一段圓弧部分組成，由于經(jīng)常與水泥柱碰觸，為防止其變形，采用直徑d 為8 mm 的彈簧鋼制成，其中圓弧部分的半徑R 與清土部件的回轉(zhuǎn)半徑相同，為275 mm，圓弧部分的弧度角α0為65°，為了保證清土部件作業(yè)時不碰到水泥柱，觸桿轉(zhuǎn)動前，其末端在水平橫向上的位置應(yīng)到達清土部件外側(cè)邊緣所在的位置，結(jié)合避障信號采集機構(gòu)的結(jié)構(gòu)參數(shù)和在擺動支架上的安裝位置，取觸桿的直線部分長度l2為505 mm，當觸桿繞O 點由P1轉(zhuǎn)動到P2位置時，此時觸桿的轉(zhuǎn)動角度θ0為16.5°，達到了最大轉(zhuǎn)動角度，由此可計算出觸桿直線部分末端由b1點運動到b2點的位移l4為145 mm。觸桿的另一端從a1點運動到a2點時，避障位移傳感器的伸長量為l3，由于STM32 單片機采集得到的避障位移傳感器的輸出信號為電壓值，并且該值隨避障位移傳感器伸長量的增加而增大，考慮到STM32 單片機采集避障位移傳感器輸出電壓值的準確性與快速性[20]，在此取避障位移傳感器的伸長量為30 mm，即l3=30 mm，由圖3 中ΔOO1a2∽ΔOO2b2可得

式中l(wèi)1為避障位移傳感器前端安裝支點到旋轉(zhuǎn)中心O 點的距離，mm；由l2=505 mm，l3=30 mm，l4=145 mm，求得l1=104 mm，即避障位移傳感器前端安裝支點到觸桿旋轉(zhuǎn)中心的距離為104 mm。

圖3 避障信號采集機構(gòu)示意圖 Fig.3 Schematic diagram of obstacle-avoiding signal acquisition mechanism

2.2 避障擺動機構(gòu)

避障擺動機構(gòu)主要由擺動支架和避障油缸組成，其中避障油缸是擺動機構(gòu)的執(zhí)行部件，其缸徑和行程是2 個關(guān)鍵的結(jié)構(gòu)參數(shù)，因此需要通過對避障擺動機構(gòu)進行動力學和運動學分析以確定其參數(shù)。

2.2.1 動力學分析

分析清土機的作業(yè)過程可知，當碰到水泥柱時，避障油缸將帶動擺動支架向清土機后方轉(zhuǎn)動以避開水泥柱，由于清土部件是逆時針旋轉(zhuǎn)，隨著清土機的前進，其后方的土壤被完全清除，故擺動機構(gòu)向清土機后方轉(zhuǎn)動避障時，由于沒有土壤阻力的作用，可認為避障油缸伸縮所需要克服的阻力主要與中間旋轉(zhuǎn)軸在軸承中轉(zhuǎn)動的摩擦力有關(guān)，由于此力較小，可忽略不計，活塞桿僅受拉壓作用，本文選取缸徑40 mm、活塞桿直徑25 mm、推力20 100 N 和拉力12 270 N 的HSG 型工程液壓缸[21-22]。

2.2.2 運動學分析

為確定避障擺動機構(gòu)中避障油缸的行程，需要對避障擺動機構(gòu)的運動狀態(tài)進行分析。為了將水泥柱行間的防寒土清除干凈，清土部件外側(cè)邊緣伸入水泥柱行間的最大距離應(yīng)大于等于水泥柱本身橫截面的寬度，同時，為保證清土部件避開水泥柱，當避障油缸完全縮回時，清土部件外側(cè)邊緣到水泥柱外邊緣的距離至少為0 mm，根據(jù)葡萄園常用的水泥柱尺寸規(guī)格，初定清土部件水平橫向運動的幅度為80 mm。

以O(shè)3為坐標原點（即中間旋轉(zhuǎn)軸的位置），將避障擺動機構(gòu)的運動簡化為二維平面運動模型，如圖4 所示。當避障油缸完全伸出時，擺動支架處于水平位置，與x軸正方向重合，此時清土部件位于E1位置，當避障油缸完全縮回時，清土部件避開水泥柱，此時清土部件位于E2位置。

分析可知，避障擺動機構(gòu)的運動滿足公式（2）～（6）：

式中l(wèi)5、l6、l7、γ 和θ 為避障擺動機構(gòu)的設(shè)計尺寸，l5=502 mm，l6=240 mm，l7=367 mm，γ=78°，θ=29°，其中清土部件橫向運動距離k=80 mm，由此求得L=135 mm?？紤]到葡萄園內(nèi)水泥柱行的直線度誤差，同時，拖拉機手操作機器作業(yè)時也存在一定偏差，故需增大避障油缸的行程，本文取L=200 mm。

圖4 避障擺動機構(gòu)運動學分析 Fig.4 Kinematics analysis of swing mechanism for obstacle-avoiding

2.3 清土部件

清土部件結(jié)構(gòu)如圖5 所示，主要由旋轉(zhuǎn)滾筒、刷子和橡膠片等組成，刷子由多層柔性塑料刷絲從旋轉(zhuǎn)滾筒底部依次向上套壓而成，共15 層，底部用壓板固定；在壓板的底面均勻地焊接固定板，通過動板與固定板的連接固定安裝橡膠片；固定部分均用螺栓連接，便于更換。由前期田間試驗可知，清土部件可有效清除葡萄藤防寒土，由于土壟下部的緊實度大于上部，故清土部件下部的磨損比上部嚴重，因此下部材料采用耐磨的夾線橡膠片。根據(jù)葡萄藤下架后壓倒在地面上的尺寸范圍和埋土后形成土壟的外形尺寸[13]，本文設(shè)計清土部件的直徑D1為550 mm，有效作業(yè)高度H 為350 mm，其中刷子部分的高度為270 mm，橡膠片的高度a 為80 mm，厚度D 為8 mm，長度b 為228 mm。底部壓板上橡膠片的數(shù)量對土壟下部的清掃范圍具有一定的影響，隨著橡膠片數(shù)量的增加，漏掃區(qū)域減小，清掃重疊區(qū)域增大[23]，為保證清土質(zhì)量，本文確定橡膠片數(shù)量為5 片。

圖5 清土部件 Fig.5 Soil cleaning parts

3 清土作業(yè)仿真試驗

由清土機的作業(yè)過程可知，機器前進速度、避障油缸速度、清土部件轉(zhuǎn)速、STM32 單片機程序中預設(shè)的觸桿轉(zhuǎn)動角度閾值（當觸桿轉(zhuǎn)動角度大于等于該值時，避障擺動機構(gòu)開始工作）和田間環(huán)境狀況都可能會對清土作業(yè)效果產(chǎn)生一定的影響。為了分析各因素對清土作業(yè)效果的影響，進而為機具的優(yōu)化及試驗參數(shù)的選取提供參考依據(jù)，本文利用多體動力學軟件RecurDyn 和離散元軟件EDEM 耦合仿真試驗，前者用來設(shè)置機具的運動參數(shù)，后者用來生成土壤顆粒模型并分析清土作業(yè)效果。

3.1 仿真模型建立

建模時若按照1∶1 的比例模擬葡萄園內(nèi)作業(yè)環(huán)境，需生成大量顆粒，由于顆粒數(shù)量對仿真時間影響巨大[24]，故為了減少仿真計算量，提高仿真效率，本文進行縮比建模試驗。首先將SolidWorks 中的機具三維模型進行簡化，簡化后的模型主要包括機架、觸桿、擺動支架、避障油缸和清土部件等，然后將簡化后的所有結(jié)構(gòu)尺寸縮小為原來的一半，保存為.igs 格式并導入RecurDyn 軟件中，同時添加相應(yīng)的材料屬性參數(shù)，根據(jù)機具的實際運動情況，為其各個部件之間添加固定副、移動副、轉(zhuǎn)動副等約束條件，并在RecurDyn中建立直徑為40 mm的圓柱體代替水泥柱。為了模擬觸桿與水泥柱間的碰撞，在觸桿與水泥柱之間添加接觸力；為了使觸桿避開水泥柱后自動復位，在觸桿與擺動支架之間創(chuàng)建線性彈簧阻尼器，同時通過設(shè)定IF 函數(shù)模擬觸桿碰水泥柱的過程，當觸桿的轉(zhuǎn)動角度大于STM32單片機程序中預設(shè)的觸桿轉(zhuǎn)動角度閾值后，避障油缸帶動清土部件避開水泥柱；在機架與地面間的移動副和清土部件與擺動支架間的旋轉(zhuǎn)副上分別添加驅(qū)動來模擬機具的前進速度v1和清土部件轉(zhuǎn)速n。

在RecurDyn 軟件中完成機構(gòu)建模和運動參數(shù)的設(shè)置后，將其中與土壤直接接觸作業(yè)的清土部件導入EDEM中[25-26]，2 個軟件中相同部件的坐標位置一致，根據(jù)清土部件在EDEM 中的位置，建立長×寬×高為3 000 mm× 2 000 mm×220 mm 的土槽作為仿真作業(yè)的環(huán)境范圍，建立土壤離散元模型時，相關(guān)模型參數(shù)參考文獻[27]，采用2 個長×寬為2 500 mm×150 mm 的四邊形作為顆粒工廠，以8 000 個/s 的速率生成初速度為2 m/s 的土壤顆粒模型，每個顆粒工廠生成40 000 個，共生成80 000 個，生成總時間為6 s，以保證顆粒的完全生成與下落沉積[28-29]，最終生成上底寬150 mm、下底寬400 mm、高125 mm、有效作業(yè)長度2 500 mm 的葡萄藤防寒土土壟模型，如圖6 所示。

圖6 仿真作業(yè)模型 Fig.6 Simulation operation model

3.2 仿真作業(yè)過程分析

仿真試驗時，設(shè)置機具位于水泥柱的一側(cè)進行作業(yè)，根據(jù)實際作業(yè)狀態(tài)，先以某一常規(guī)工況進行仿真作業(yè)：在RecurDyn 軟件中設(shè)置機具前進速度為0.3 m/s，避障油缸速度為80 mm/s，清土部件轉(zhuǎn)速為450 r/min，觸桿轉(zhuǎn)動角度閾值為8°，同時為保證仿真的連續(xù)性，設(shè)置仿真的時間步為500，總時間為10 s，由仿真試驗可知，在機具前進過程中，不僅能夠使清土部件順利避開水泥柱，而且可將水泥柱側(cè)邊和水泥柱行間的土壤清除，圖7a 為機具避開水泥柱時的狀態(tài)，圖7b 為機具避開水泥柱后的正常清土作業(yè)狀態(tài)。

圖7 仿真清土作業(yè)過程 Fig.7 Simulation of soil cleaning operation process

為分析影響清土機作業(yè)效果的主要因素，進而得到理想的作業(yè)參數(shù)組合，在仿真模型基礎(chǔ)上，進行多因素、多水平試驗研究。為提高試驗的可控性與準確性，選取機具前進速度、避障油缸速度、清土部件轉(zhuǎn)速和觸桿轉(zhuǎn)動角度閾值為試驗因素，由于清土機中液壓部分的液壓油來自拖拉機后側(cè)液壓輸出口，在液壓油輸出流量的限制下，避障油缸實際速度不超過100 mm/s，根據(jù)仿真預試驗，在避障油缸速度設(shè)置為最大值的情況下，當機具前進速度大于0.6 m/s 時，清土部件來不及避開水泥柱，當避障油缸速度較小、觸桿轉(zhuǎn)動角度閾值較大時，也會出現(xiàn)同樣的問題，考慮到機具的作業(yè)效率和避障效果，選取機具前進速度范圍為0.3～0.6 m/s、避障油缸速度范圍為60～100 mm/s，觸桿轉(zhuǎn)動角度閾值范圍為5～10°。此外，當清土部件轉(zhuǎn)速較小時，被清除的部分土壤會旋回土壟，達不到有效清土的目的，但當清土部件轉(zhuǎn)速較高時，又會加大機器的振動，故為了保證清土機的作業(yè)效果和運行的平穩(wěn)性，確定清土部件的轉(zhuǎn)速范圍為350～550 r/min，仿真作業(yè)結(jié)束后，以土壤清除率作為評價指標，土壤清除率的計算公式為

式中C 為土壤清除率，%；N 為清土前生成的土壤顆?？倲?shù)（本研究中N=80 000）；N1為清土后土壟上剩余的土壤顆?？倲?shù)。

3.3 仿真優(yōu)化試驗

以機具前進速度v1、避障油缸速度v2、清土部件轉(zhuǎn)速n 和觸桿轉(zhuǎn)動角度閾值θ 為試驗因素，以土壤清除率C為試驗指標，結(jié)合仿真預試驗確定的各試驗因素范圍，采用4 因素二次回歸正交旋轉(zhuǎn)中心組合模擬試驗，由試驗設(shè)計原則，確定編碼系數(shù)r 為2.0[30-31]，試驗因素編碼結(jié)果如表2 所示。

表2 清土作業(yè)仿真試驗因素編碼表 Table 2 Code table for simulation test factors of soil cleaning operation

共進行31 組試驗（0 水平試驗重復7 次），試驗方案與結(jié)果如表3 所示，采用Design-expert 8.0.6 軟件對試驗數(shù)據(jù)進行分析，采用Quadratic 模型建立土壤清除率C的回歸模型，回歸模型的顯著性檢驗結(jié)果如表4 所示。

根據(jù)土壤清除率C 的分析結(jié)果可知，回歸模型的P值小于0.001，而失擬項的P 值大于0.05，說明回歸模型極其顯著且失擬不顯著，回歸有效。在4 個試驗因素中，X1、X2、X3、X4項的系數(shù)均達到顯著性水平；交互項因素中只有X3X4項對土壤清除率影響顯著，其余因素的交互影響不顯著；二次方項因素中，X22對土壤清除率影響顯著，其余影響不顯著。將不顯著項刪除后得到土壤清除率C 的回歸方程為

表3 清土作業(yè)仿真試驗方案與結(jié)果 Table 3 Simulation test scheme and results of soil cleaning operation

表4 回歸模型顯著性分析 Table 4 Significant analysis of regression model

為得到理想的作業(yè)參數(shù)組合，利用Design-expert 8.0.6軟件Optimization-Numerical 模塊對回歸模型進行優(yōu)化，以土壤清除率最大為尋優(yōu)目標，得到4 個試驗因素的優(yōu)化結(jié)果為機具前進速度0.43 m/s、避障油缸速度60 mm/s、清土部件轉(zhuǎn)速550 r/min 和觸桿轉(zhuǎn)動角度閾值9.96°，此時最佳土壤清除率理論值為54.33%。為便于后期在田間試驗中準確測定并驗證觸桿轉(zhuǎn)動角度閾值大小，將其圓整為10°，根據(jù)優(yōu)化得到的作業(yè)參數(shù)組合再次進行仿真試驗，得到土壤清除率為54.65%，與理論值基本一致。

4 田間驗證試驗

田間試驗于2019 年10 月底在山東高密益豐機械有限公司試驗地進行，土壤質(zhì)地為砂壤土，含水率6%～8%，類似于葡萄園內(nèi)的葡萄藤防寒土，試驗前，采用直徑為40 mm的鐵管代替水泥柱，按行距3 m、前后間距3 m 插在地上，共2 行，每行5 根，模擬葡萄園內(nèi)的水泥柱行，其中一行用來試驗機具的避障清土作業(yè)過程，另一行用來模擬葡萄行，以驗證擋土板能否將清土部件清除的土壤全部擋在行間。為了使土壟具有的一定堅實度（葡萄園內(nèi)土壟從上部到底部堅實度為69～135 kPa，并逐漸增大），利用葡萄藤埋土機從兩邊取土覆蓋到其中一行的圓管鐵柱附近，并用鐵鍬將土壟周圍進行修整與壓實，最終形成類似于葡萄園內(nèi)的葡萄藤埋土土壟，土壟橫截面近似于等腰梯形，上底寬約為300 mm，下底寬約為800 mm，高約為250 mm，如圖8 所示。試驗設(shè)備與儀器包括五征NS604C 拖拉機、自動避障式葡萄藤防寒土清土機、電腦、數(shù)據(jù)采集卡NI USB-6008（美國國家儀器有限公司，模擬輸入采樣速率為 10 kHz，采集電壓范圍為-10～10 V，精度為84.8 mV）、智能數(shù)字轉(zhuǎn)速表（分辨率0.1 r/min，非接觸式激光測量，激光射程有效距離50～500 mm）、SC-900 土壤緊實度儀（測量壓力范圍為0～7 000 kPa，測量深度為0～450 mm）、TDR 150 土壤含水率測定儀（精度±3%，分辨率0.1%）、數(shù)顯角度尺（量程0～360°，分辨率0.05 度）、卷尺（浙江奧奔工具有限公司，0～5 m，精度為1 mm）、鋼尺（量程50 cm，精度1 mm）、手機（用作計時秒表）等。

圖8 田間清土試驗 Fig.8 Field soil cleaning test

試驗前，在拖拉機靜止狀態(tài)下，在電腦在程序中設(shè)定觸桿轉(zhuǎn)動角度閾值為10°，并通過USB線下載到STM32單片機中，采用數(shù)顯角度尺測量觸桿的實際轉(zhuǎn)動角度。利用電腦處理數(shù)據(jù)采集卡測得的同步位移傳感器的輸出信號計算避障油缸速度，采用智能數(shù)字轉(zhuǎn)速表測量清土部件的轉(zhuǎn)速，通過控制拖拉機油門的大小使避障油缸速度和清土部件轉(zhuǎn)速分別穩(wěn)定在60 mm/s和550 r/min左右，并在作業(yè)時保持。試驗時，利用秒表測定拖拉機勻速行駛一定距離的時間，以計算其前進速度，使其穩(wěn)定在0.43 m/s 左右，使機具的田間作業(yè)參數(shù)與仿真優(yōu)化試驗保持一致。田間試驗的土壤清除率采用作業(yè)前和作業(yè)后土壟橫截面面積的變化量來衡量（作業(yè)前，土壟橫截面形狀近似為等腰梯形，計算到其面積為137 500 mm2），計算公式如下：

式中C＇為田間試驗土壤清除率，%；S 為清土作業(yè)后的土壟橫截面面積，mm2.

試驗采用多點測量法，在每一次土壟清除試驗完成后，隨機選擇10 個測量點，包括圓管鐵柱之間和圓管鐵柱附近的位置，通過卷尺和鋼尺分別測量出這10 個測量點位置處的土壟橫截面尺寸，并計算其面積大小，共重復3 次清土試驗，得到30 組數(shù)據(jù)，計算其平均值。試驗得到清土機單邊清土作業(yè)時的土壤清除率C＇為59.73%，與仿真優(yōu)化試驗結(jié)果的相對誤差為8.50%，產(chǎn)生誤差的原因除了與土壤條件差異、土壤中的雜物、地面平整度等有關(guān)，還受到拖拉機行駛穩(wěn)定性的影響。清土作業(yè)效果如圖9 所示，從圖中剩余土壟的邊緣可明顯觀察到清土部件的避障運動軌跡，同時可看出水泥柱側(cè)邊和水泥柱行間的土壤被有效清除，此外，清土部件清除的土壤都被擋土板擋在了行間。總體上，機具在田間以仿真優(yōu)化的工作參數(shù)組合進行穩(wěn)定作業(yè)能夠較好地完成葡萄藤防寒土的清除。

圖9 田間試驗清土作業(yè)效果 Fig.9 Soil cleaning operation effects of field test

5 結(jié)論與討論

1）設(shè)計了一種主要由機架、避障擺動機構(gòu)、避障信號采集機構(gòu)、控制器、清土部件、擋土板、傳動部件和限深輪等組成的自動避障式葡萄藤防寒土清土機，該機器中的清土部件采用柔性刷子和橡膠組合結(jié)構(gòu)，通過觸桿感知自動避開水泥柱，可將水泥柱側(cè)邊和水泥柱行間的土壤全部清除干凈，確定了避障信號采集機構(gòu)、避障擺動機構(gòu)和清土部件等關(guān)鍵部件的結(jié)構(gòu)及參數(shù)，其中避障油缸行程為200 mm、缸徑為40 mm，清土部件的回轉(zhuǎn)直徑為550 mm。

2）基于多體動力學軟件RecurDyn 和離散元軟件EDEM 耦合仿真，以機具前進速度、避障油缸速度、清土部件轉(zhuǎn)速和觸桿轉(zhuǎn)動角度閾值為試驗因素，以土壤清除率為評價指標，設(shè)計了四因素二次回歸正交旋轉(zhuǎn)中心組合模擬試驗，利用Design Expert 8.0.6 軟件對仿真試驗數(shù)據(jù)進行回歸模型顯著性分析，機具前進速度、避障油缸速度、清土部件轉(zhuǎn)速和觸桿轉(zhuǎn)動角度閾值均對仿真試驗結(jié)果影響顯著，以土壤清除率最大化為優(yōu)化目標值，采用響應(yīng)面法獲得了機具的最優(yōu)工作參數(shù)組合：機具前進速度0.43 m/s、避障油缸速度60 mm/s、清土部件轉(zhuǎn)速550 r/min、觸桿轉(zhuǎn)動角度閾值10 °，在優(yōu)化后的作業(yè)參數(shù)組合下再次進行仿真試驗，得到土壤清除率理論值為54.65%。

3）加工物理樣機并進行田間試驗，采用優(yōu)化后的工作參數(shù)組合進行田間作業(yè)性能測試，清土作業(yè)完成后，依據(jù)作業(yè)前后土壟橫截面面積的變化，通過多點測量，計算得出田間試驗土壤清除率的平均值為59.73%，與仿真優(yōu)化結(jié)果的相對誤差為8.50%。

該文試驗是在土壟中沒有埋入葡萄藤的情況下進行的，主要分析了自動避障式葡萄藤防寒土清土機的避障效果和土壤清除率，而在實際葡萄園中的生產(chǎn)作業(yè)中，入冬前先將修剪后的葡萄藤下架壓在地面，然后采用葡萄藤埋土機進行覆土，壓的葡萄藤的土壟高度和寬度各異，有的位于水泥柱行間，有的則超出水泥柱截面范圍，由于拖拉機的行走誤差，在葡萄園內(nèi)進行清土作業(yè)時，清土部件不可避免地會與埋在土中的葡萄藤接觸，故清土部件需在保證清除土壤的同時不能損傷葡萄藤，此外，材料的耐磨性也必須可靠，下一步將深入分析清土部件材料與土壤和葡萄藤之間的相互接觸作用。