何浩猛，胡 斌,2※，潘 峰，羅 昕,2，郭孟宇,2，解彥宇，陳學(xué)庚,2

（1. 石河子大學(xué)機(jī)械電氣工程學(xué)院，石河子 832003；2. 農(nóng)業(yè)農(nóng)村部西北農(nóng)業(yè)裝備重點(diǎn)實(shí)驗(yàn)室，石河子 832003）

0 引言

機(jī)械化回收地膜已成為治理新疆棉田白色污染的重要途徑和有效方法，其工作特點(diǎn)是將地膜連同棉稈、棉鈴、棉殼和土壤等雜質(zhì)一同回收[1-2]。因膜雜混合物中地膜和棉稈（簡稱“地膜與棉稈”）相互纏繞打結(jié)不易分離，且缺少有效的分離方法及設(shè)備等問題，回收后的地膜無法進(jìn)一步實(shí)現(xiàn)回收再利用，多采用棄置田間地頭或焚燒處理，間接對環(huán)境造成了二次污染[3]。

國外地膜回收均以卷收式為主且含雜量較低，可直接回收再利用，尚未有關(guān)此方面研究[4]。針對地膜收而未用等問題，國內(nèi)學(xué)者主要從兩方面開展研究，一類是通過研制聯(lián)合作業(yè)機(jī)在收膜過程實(shí)現(xiàn)地膜和雜質(zhì)相分離，如陳學(xué)庚院士團(tuán)隊(duì)研制的隨動式殘膜回收聯(lián)合作業(yè)機(jī)，該機(jī)型已取得了不錯(cuò)田間試驗(yàn)效果，但目前尚未進(jìn)行市場推廣[5-7]；另一類是對回收集堆后的膜雜混合物先粉碎破結(jié)再進(jìn)行機(jī)械分離，如石鑫等[8]和彭強(qiáng)吉等[9]根據(jù)地膜、棉稈和土塊等風(fēng)力懸浮速度差異，設(shè)計(jì)并研制了滾筒篩式風(fēng)選裝置，機(jī)具試驗(yàn)發(fā)現(xiàn)分離效率低但功耗高。胡斌[10]和羅昕[11]從研究機(jī)收地膜與棉稈混合物物料物理參數(shù)屬性出發(fā)，基于固液兩相流運(yùn)動基本理論，提出了以水介質(zhì)作用，實(shí)現(xiàn)膜雜固固分離和清洗的方法。李俊虹等[12-13]設(shè)計(jì)了一種水射流和葉輪驅(qū)動復(fù)合式機(jī)收膜雜水洗分離裝置，可基本實(shí)現(xiàn)殘膜與雜質(zhì)的分離，結(jié)果表明在研究地膜與棉稈水洗分離時(shí)可把地膜與棉稈物料屬性、外部流場和機(jī)械裝置三者視為一個(gè)統(tǒng)一系統(tǒng)，其相互作用和影響，即地膜與棉稈水洗分離是地膜與棉稈物料屬性、外部流場和機(jī)械裝置的復(fù)合函數(shù)作用的過程，但其相互影響關(guān)系并未得到證明。

膜雜分離研究處于初步階段，通過水介質(zhì)作用實(shí)現(xiàn)地膜的清洗和分離的研究更為鮮見。又因膜雜水洗分離流場域環(huán)境復(fù)雜不可控且影響因素多，現(xiàn)今缺乏相關(guān)基礎(chǔ)性研究成果，故在室內(nèi)可控條件下進(jìn)行重復(fù)性機(jī)理試驗(yàn)是研究膜雜水洗分離運(yùn)動規(guī)律的有效方法。本研究通過借鑒造紙纖維懸浮液動力學(xué)和河流泥沙動力學(xué)領(lǐng)域相對成熟理論分析方法，研究了外部流場和物料屬性對地膜與棉稈分離效果的影響，所得結(jié)果可為后續(xù)的機(jī)械裝置設(shè)計(jì)提供相關(guān)參數(shù)。進(jìn)一步地通過錄像后期處理，將地膜與棉稈在沉降聚集過程中的姿態(tài)、軌跡變換情況及地膜與棉稈分離全過程進(jìn)行系統(tǒng)再現(xiàn)，所得結(jié)論可為地膜與棉稈分離方案的實(shí)際應(yīng)用提供相關(guān)理論基礎(chǔ)。

1 材料與方法

1.1 地膜與棉稈水洗分離理論基礎(chǔ)

當(dāng)?shù)啬づc棉稈投入裝水的立式圓筒后，由表1可知其密度和比重均比水小，兩者皆漂浮在水面。當(dāng)對圓筒中的水介質(zhì)施加機(jī)械擾動時(shí)，由拜爾定律和茶葉悖論知[14-16]，棉稈在水介質(zhì)挾裹下運(yùn)動至圓筒中部中央附近，地膜則運(yùn)動至底部中央附近。當(dāng)機(jī)械擾動停止后，地膜與棉稈在浮力、壓力梯度力、Magnus升力、Saffman升力和離心力等作用下向水面作聚集恢復(fù)運(yùn)動。

表1 地膜和棉稈物理參數(shù)對比 Table 1 Comparison of physical parameters between plastic film and cotton stalk

棉稈屬細(xì)長剛性圓柱體，體積大，具有三維旋轉(zhuǎn)自由度，力學(xué)響應(yīng)周期短，可快速達(dá)到自身運(yùn)動穩(wěn)定狀態(tài)，不易受外界流場影響而改變；地膜屬不規(guī)則柔性體，體積較小但與水接觸面積大，相對水介質(zhì)有較大跟隨性，力學(xué)響應(yīng)周期長，易受流場環(huán)境影響而改變自身運(yùn)動軌跡。故此可知地膜與棉稈的外觀體積和抗彎強(qiáng)度等物理屬性不同，兩者對水介質(zhì)傳遞力作用時(shí)有不同的恢復(fù)響應(yīng)差異。在二維垂向平面時(shí)，這種差異具體表現(xiàn)為響應(yīng)時(shí)間差和位移差。參考于此，可以借助因地膜與棉稈在外部流場中因物料物理參數(shù)差異而產(chǎn)生的時(shí)間差和位移差對地膜與棉稈混合物進(jìn)行分時(shí)或分層分離回收。

1.2 試驗(yàn)裝置

試驗(yàn)裝置主要由葉輪、聯(lián)軸器、葉輪軸、立式亞克力圓筒（直徑為200 mm、高為500 mm、壁厚為5 mm，上?？掠袡C(jī)玻璃制作中心）、57步進(jìn)電機(jī)（J-5718HB2401，溫州普菲德電氣有限公司）及配套驅(qū)動器和直流電源、可編程邏輯控制器（Programmable Logic Controllers，PLC）和EOS M6佳能數(shù)碼相機(jī)（日本）組成（圖1）。立式亞克力圓筒為透明且貼有標(biāo)尺刻度，便于錄像、觀察液面高度讀數(shù)及確定地膜與棉稈物料位置；葉輪軸通過聯(lián)軸器和螺栓分別與步進(jìn)電機(jī)和葉輪實(shí)現(xiàn)可拆卸式連接，試驗(yàn)中通過更換葉輪軸實(shí)現(xiàn)葉輪在不同位置擾動作用；通過三菱觸摸屏畫面設(shè)計(jì)軟件GT Designer繪制操作界面寫入三菱觸摸屏；利用三菱觸摸屏編程軟件GX Works編寫程序?qū)懭隤LC，實(shí)現(xiàn)對步進(jìn)電機(jī)的啟停、擾動時(shí)長和轉(zhuǎn)速大小控制；直流電源為步進(jìn)電機(jī)提供動力。

試驗(yàn)用水提前12 h置入圓筒中，以使氣泡消失并達(dá)到熱平衡。試驗(yàn)時(shí)將地膜與棉稈投入圓筒內(nèi)，啟動步進(jìn)電機(jī)帶動葉輪轉(zhuǎn)動，試驗(yàn)由數(shù)碼相機(jī)（圖像分辨率為1280×720，幀速率為每秒50幀）錄像，拍攝時(shí)間為葉輪開始轉(zhuǎn)動到地膜與棉稈恢復(fù)至相對穩(wěn)定狀態(tài)。每次試驗(yàn)后關(guān)閉電源停機(jī)，待地膜與棉稈恢復(fù)靜止，再進(jìn)行重復(fù)試驗(yàn)。

1.3 試驗(yàn)材料

為便于試驗(yàn)觀測，優(yōu)先選擇埋鋪黑色地膜的棉田進(jìn)行取樣。2019年10月對石河子市145團(tuán)3分場1LMLG-7型立桿摟膜機(jī)作業(yè)后的膜雜混合物隨機(jī)取樣3份，每份1000 g。棉花種植品種為新陸早49號，雜質(zhì)中棉稈直徑分布為4～5 mm，其含水率 ＜ 20%；地膜厚度為0.01 mm，田間使用時(shí)長約180 d。為減小棉稈含水率對地膜與棉稈分離效果的影響，每組材料只進(jìn)行1次試驗(yàn)，不作重復(fù)使用。

1.4 試驗(yàn)評價(jià)指標(biāo)

參照造紙纖維懸浮液動力學(xué)[20-21]、河流泥沙運(yùn)動力學(xué)理論研究[22-24]以及殘膜作業(yè)機(jī)膜雜分離質(zhì)量要求[6-7]，選用地膜與棉稈馳豫時(shí)間差（Δt，s）和最大垂向位移差（Δd，cm）作為評價(jià)指標(biāo)，其計(jì)算如式（1）～式（4）

式中t1和t3分別為地膜下降至最深時(shí)和恢復(fù)至穩(wěn)定狀態(tài)的時(shí)刻，s；t2和t4分別為棉稈下降至最深時(shí)和恢復(fù)至穩(wěn)定狀態(tài)的時(shí)刻，s；d1和d2分別為t1時(shí)刻下地膜和棉稈的在圓筒內(nèi)的刻度值，cm；td為地膜馳豫時(shí)長，s；tm為棉稈馳豫時(shí)長，s。為便于試驗(yàn)數(shù)據(jù)處理與統(tǒng)計(jì)，以地膜與棉稈的質(zhì)心為基準(zhǔn)，通過圓筒刻度估讀其垂向位移，不保留小數(shù)位。

1.5 單因素試驗(yàn)

在圓筒水深為350 mm，葉輪位于圓筒刻度線30 cm處，擾動時(shí)長為10 s，轉(zhuǎn)速為200 r/min，葉輪葉數(shù)為3葉，棉稈直徑約為5 mm、長60 mm[25]，長條形黑色地膜面積為1413 mm2的條件下進(jìn)行試驗(yàn)，通過錄像后處理計(jì)算其地膜與棉稈馳豫時(shí)間差和最大垂向位移差。

采用控制變量法和對比分析方法，以地膜與棉稈馳豫時(shí)間差或地膜與棉稈最大垂向位移差為響應(yīng)指標(biāo)，分別進(jìn)行4組如下的外部流場特性單因素試驗(yàn)：1）葉輪分別位于圓筒刻度線35、30、25、20、15和10 cm處；2）擾動時(shí)長分別為5、10、15、20和25 s；3）葉輪轉(zhuǎn)速分別為100、150、200、250和300 r/min；4）葉數(shù)分別為2、3和4葉；每組試驗(yàn)重復(fù)3次，取其平均值作為試驗(yàn)結(jié)果。3組地膜與棉稈物料屬性單因素試驗(yàn)：1）取面積為1413 mm2的地膜，分別裁為半徑為21.2 mm的圓形，長為45.0 mm、寬為31.4 mm的長條形，兩對角線分別為50.0和40.4 mm的正菱形和邊長為57.0 mm的正三角形；2）保持地膜寬度31.4 mm不變，分別取面積為471、942、1413、1884、2355和2826 mm2長條形地膜；3）取直徑約為5 mm，長分別為20、40、60、80和100 mm的棉稈，其對應(yīng)長徑比[26]分別為4、8、12、16和20；每組試驗(yàn)重復(fù)5次，取其平均值作為試驗(yàn)結(jié)果。采用數(shù)據(jù)處理軟件Design Expert V10.0.4下General Factorial模塊對各單因素試驗(yàn)數(shù)據(jù)進(jìn)行方差分析，得到響應(yīng)指標(biāo)的檢驗(yàn)統(tǒng)計(jì)量（FA）。通過F分布表查得各顯著性水平（α）下的臨界值（Fα），并根據(jù)FA和Fα的相對大小，判斷各因素對試驗(yàn)指標(biāo)影響是否顯著。

1.6 正交試驗(yàn)

1.6.1 試驗(yàn)方案設(shè)計(jì)

根據(jù)單因素試驗(yàn)結(jié)果，以地膜與棉稈馳豫時(shí)間差和最大垂向位移差為響應(yīng)指標(biāo)，選定顯著性影響因素：葉輪轉(zhuǎn)速（A, r/min）、地膜面積（B, mm2）和棉稈長徑比（C）為自變量，進(jìn)行3因素5水平中心組合設(shè)計(jì)（Central Composite Design，CCD）正交試驗(yàn)。根據(jù)響應(yīng)曲面法中的CCD試驗(yàn)設(shè)計(jì)要求[27-28]，共進(jìn)行20組試驗(yàn)，其中析因部分試驗(yàn)次數(shù)為8，星號點(diǎn)次數(shù)為6，中心點(diǎn)重復(fù)次數(shù)為6，每組試驗(yàn)重復(fù)3次，取其平均值作為該組試驗(yàn)結(jié)果。采用Design Expert V10.0.4軟件進(jìn)行數(shù)據(jù)處理，建立試驗(yàn)因素與響應(yīng)指標(biāo)間的擬合回歸方程。并通過軟件對回歸擬合方程的檢驗(yàn)統(tǒng)計(jì)量F值、P值、校正決定系數(shù)（adjusted coefficient of determination，RAdj2）和變異系數(shù)（Coefficient of Variance，CV）進(jìn)行計(jì)算并做出評價(jià)，判斷是否可進(jìn)行后續(xù)優(yōu)化研究。F值和P值表示其相關(guān)系數(shù)的顯著性[28]，F(xiàn)值越大，P值越小，其顯著性越高，當(dāng)P＜ 0.01為極顯著；RAdj2值越大，說明擬合回歸方程的相關(guān)性越好[29]；CV值表示試驗(yàn)的精確度，其值越小，代表試驗(yàn)可靠性越高[29]。

各試驗(yàn)因素的實(shí)際水平值通過式（5）進(jìn)行編碼[27]，其結(jié)果如表2所示。

式中zj為因素水平的編碼值；xj為因素水平的實(shí)際值；xj0、xj1、xj2、xjγ、x-jγ分別為試驗(yàn)因素xj的零水平、上水平、下水平、上星號臂水平和下星號臂水平；Δj為試驗(yàn)因素xj的變化間距；γ為星號臂長，根據(jù)試驗(yàn)因素個(gè)數(shù)和中心點(diǎn)次數(shù)，γ取值為1.682[27]。

表2 中心組合設(shè)計(jì)的試驗(yàn)因素水平及編碼 Table 2 Level and coding of experimental factors for the Central Composite Design (CCD)

1.6.2 參數(shù)優(yōu)化與驗(yàn)證

為獲得地膜與棉稈分離最佳的外部流場和地膜與棉稈物料特性參數(shù)，以地膜與棉稈最大馳豫時(shí)間差和最大垂向位移差為優(yōu)化目標(biāo)，利用Design Expert V10.0.4軟件下Optimization-Numerical模塊，采用期望函數(shù)法進(jìn)行參數(shù)優(yōu)化求解，其約束條件和目標(biāo)函數(shù)如式（6）所示

為檢驗(yàn)參數(shù)優(yōu)化求解后的可靠性，運(yùn)用所得的優(yōu)化參數(shù)組合在地膜與棉稈懸浮水力特性試驗(yàn)臺進(jìn)行試驗(yàn)驗(yàn)證，每組試驗(yàn)重復(fù)5次，取其平均值作為試驗(yàn)結(jié)果，并對比優(yōu)化求解的預(yù)測值與試驗(yàn)平均值間的相對誤差判斷其是否可靠。

1.7 地膜與棉稈馳豫過程分析

通過視頻編輯軟件Adobe Premiere CC 2015對錄像進(jìn)行逐幀分析，依次記錄葉輪啟停時(shí)刻、地膜與棉稈向下沉降和向上聚集時(shí)刻、以及其運(yùn)動軌跡和姿態(tài)變化的情況。研究內(nèi)容為1）葉輪轉(zhuǎn)速對地膜與棉稈分離效果影響的單因素試驗(yàn)時(shí)分別在150和250 r/min試驗(yàn)條件下，對錄像每15幀截取1張圖片，觀測并記錄葉輪轉(zhuǎn)動時(shí)棉稈向下沉降和向上聚集馳豫過程；2）地膜相同形狀下，不同面積對其馳豫時(shí)間影響的單因素試驗(yàn)時(shí)，對錄像每30幀截取1張圖片，觀測并記錄地膜的復(fù)雜運(yùn)動過程；3）在正交試驗(yàn)條件下，對錄像每35幀截取1張圖片，觀測并記錄葉輪停轉(zhuǎn)后地膜的恢復(fù)運(yùn)動姿態(tài)及軌跡變化；以及對錄像每30幀截取1張圖片，觀測并研究地膜與棉稈沉降聚集馳豫分離全過程。

2 結(jié)果與分析

2.1 外部流場特性單因素試驗(yàn)結(jié)果分析

2.1.1 葉輪擾動位置對地膜與棉稈分離的影響

隨著葉輪擾動位置在水中的變化，地膜與棉稈馳豫時(shí)間差和位移差呈現(xiàn)出先升后降趨勢（圖2）。當(dāng)葉輪位于刻度線30 cm（距水面5 cm）左右時(shí)，最優(yōu)地膜與棉稈馳豫時(shí)間差和位移差分別為3.92 s和20 cm。葉輪擾動位置靠近筒底時(shí)，地膜與棉稈只能跟隨表層水形成的渦流運(yùn)動，且膜在桿的上方；擾動位置靠近水面時(shí)，葉輪在水表面旋轉(zhuǎn)濺起水浪，達(dá)不到地膜與棉稈分離效果。

2.1.2 葉輪擾動時(shí)長對地膜與棉稈分離的影響

隨擾動時(shí)長的變化，地膜與棉稈馳豫時(shí)間差先增后平緩，位移差隨擾動時(shí)長先增后減（圖3）。當(dāng)擾動時(shí)長在10 s左右時(shí)，最佳地膜與棉稈馳豫時(shí)間差為4.62 s，位移差為23 cm。擾動時(shí)間過短時(shí)，地膜與棉稈來不及向下沉降但葉輪擾動已結(jié)束；當(dāng)?shù)啬づc棉稈達(dá)到穩(wěn)定狀態(tài)后，馳豫時(shí)間和位移不再隨擾動時(shí)長增長而變化。

2.1.3 葉輪轉(zhuǎn)速對地膜與棉稈分離的影響

隨葉輪轉(zhuǎn)速的增大，地膜與棉稈馳豫時(shí)間差先增大至平緩再減小，位移差先增后減（圖4）。當(dāng)葉輪轉(zhuǎn)速為200 r/min時(shí)，最優(yōu)馳豫時(shí)間差和位移差分別為4.32 s和19 cm。葉輪轉(zhuǎn)速低于200 r/min時(shí)，擾動強(qiáng)度過小不足以形成渦流，此時(shí)地膜無法越過葉輪向下沉降；葉輪轉(zhuǎn)速高于200 r/min時(shí)，擾動強(qiáng)度過大，地膜在葉輪和水面的空間作振蕩運(yùn)動或吸附于葉輪隨其一起運(yùn)動，地膜與棉稈分離效果不明顯。單因素方差分析可得葉輪轉(zhuǎn)速對其馳豫時(shí)間差和位移差的檢驗(yàn)統(tǒng)計(jì)量FA值分別為18.11和12.64。當(dāng)顯著性水平α= 0.01時(shí)，F(xiàn)0.01(2,6) = 10.92，F(xiàn)A＞F0.01，故葉輪轉(zhuǎn)速是影響地膜與棉稈馳豫時(shí)間差和位移差的極顯著因素。

2.1.4 葉輪葉數(shù)對地膜與棉稈分離的影響

當(dāng)葉數(shù)為3葉時(shí)地膜與棉稈馳豫時(shí)間差和垂向位移差的試驗(yàn)效果最優(yōu)，分別為3.94 s和21 cm（表3）。相同條件下2葉葉輪比3葉產(chǎn)生的渦流小，無法滿足地膜與棉稈分離所需渦流強(qiáng)度要求；4葉時(shí)可達(dá)到足夠的渦流強(qiáng)度，但因葉輪截面面積過大，地膜越過葉輪平面時(shí)有更大概率碰撞葉輪并吸附表面隨其運(yùn)動，造成分離效果不明顯。

2.2 地膜與棉稈物料屬性單因素試驗(yàn)結(jié)果分析

2.2.1 相同面積下地膜形狀對其馳豫時(shí)間影響

地膜馳豫時(shí)間隨其形狀變化而波動（圖5），但影響效果不顯著。圓形地膜馳豫時(shí)間略小于長條形、菱形和三角形地膜，三角形地膜馳豫時(shí)間波動較大，這是因?yàn)槠湓谟鞣较虻南鄬﹂L度不同，迎流長度越長越易達(dá)到穩(wěn)定狀態(tài)，馳豫時(shí)間波動越小。長條形和菱形地膜馳豫時(shí)間較長且穩(wěn)定，考慮到試驗(yàn)中長條形更易于量化測試，實(shí)際生產(chǎn)中長條形地膜的破結(jié)粉碎更易于實(shí)現(xiàn)，所以后續(xù)研究選用長條形地膜。

2.2.2 相同形狀下，地膜面積對其馳豫時(shí)間影響

地膜的不同面積是影響其馳豫時(shí)間的重要因素。當(dāng)?shù)啬っ娣e分別為942、1413和1884 mm2時(shí)，其馳豫時(shí)間保持在7.00～7.60 s，波動較?。幻娣e為471 mm2的地膜馳豫時(shí)間波動較大，這是因?yàn)榈啬っ娣e過小時(shí)，更易在筒底和水表面間振蕩，且運(yùn)動軌跡隨機(jī)性更大；面積為2826 mm2的地膜更快達(dá)到穩(wěn)定狀態(tài)，馳豫時(shí)間最短，分析可知地膜面積越大，其與水介質(zhì)間摩擦力越大，越易達(dá)到穩(wěn)定狀態(tài)（圖6）。單因素方差分析可得地膜面積對其馳豫時(shí)間影響的檢驗(yàn)統(tǒng)計(jì)量FA值為10.28。當(dāng)顯著性水平α=0.05時(shí)，F(xiàn)0.05(2,12) = 3.89，F(xiàn)A＞F0.05，故地膜面積是影響地膜與棉稈馳豫時(shí)間差的顯著因素。

表3 不同葉輪葉數(shù)對馳豫時(shí)間差和垂向位移差的影響 Table 3 Effect of different number of impellers on time difference and vertical displacement difference

2.2.3 棉稈長徑比對其馳豫時(shí)間影響

棉稈長徑比對其馳豫時(shí)間影響較大且隨其變化而波動（圖7），當(dāng)棉稈長徑比分別為8、12和16時(shí)，其馳豫時(shí)間均約為3.30 s，其中當(dāng)長徑比為12時(shí)棉稈馳豫時(shí)間較穩(wěn)定。長徑比為4和20的棉稈馳豫時(shí)間均約為2.10 s，但長徑比為4的棉稈馳豫時(shí)間波動較大，長徑比為20的棉稈馳豫時(shí)間最穩(wěn)定，這種現(xiàn)象也是因棉稈不同體積與水介質(zhì)間不同摩擦力造成的。單因素方差分析可得棉稈長徑比對其馳豫時(shí)間的檢驗(yàn)統(tǒng)計(jì)量FA值為6.10。當(dāng)顯著性水平α=0.05時(shí)，查F分布表知F0.05(2,12) = 3.89，F(xiàn)A＞F0.05，故棉稈長徑比是影響地膜與棉稈馳豫時(shí)間的顯著因素。

2.3 正交試驗(yàn)結(jié)果分析

2.3.1 試驗(yàn)結(jié)果及回歸方差分析

CCD正交試驗(yàn)設(shè)計(jì)方案及其結(jié)果如表4所示。通過對表4中試驗(yàn)的響應(yīng)值進(jìn)行回歸方差分析，分別建立葉輪轉(zhuǎn)速、地膜面積和棉稈長徑比與地膜與棉稈馳豫時(shí)間差（Δt，s）和最大位移差（Δd，cm）之間的擬合回歸方程為

由地膜與棉稈馳豫時(shí)間差擬合回歸方程一次項(xiàng)系數(shù)可知，在所選取的因素水平范圍內(nèi)各影響因素的檢驗(yàn)統(tǒng)計(jì)量由大到小依次為FA、FB和FC，故各因素對地膜與棉稈馳豫時(shí)間差的影響由大到小依次為葉輪轉(zhuǎn)速、地膜面積和棉稈長徑比；同理可知影響地膜與棉稈最大垂向位移差的因素由大到小依次為葉輪轉(zhuǎn)速、棉稈長徑比和地膜面積。

在用上述擬合回歸方程預(yù)測各因素與響應(yīng)值之間的變化關(guān)系時(shí)（表5），可知其檢驗(yàn)統(tǒng)計(jì)量F值分別為46.38和22.38，相對應(yīng)的P值均 ＜ 0.0001，說明該擬合方程的相關(guān)系數(shù)是極顯著的?；貧w擬合方程的RAdj2值分別為0.956和0.910，說明該方程可解釋95.6%和91.0%的地膜與棉稈馳豫時(shí)間差和垂向最大位移差的波動變化，擬合性良好。CV值分別為1.37%和12.8%，試驗(yàn)操作可信。

表4 中心組合設(shè)計(jì)試驗(yàn)表及結(jié)果 Table 4 Center composite design matrix and experimental results

表5 擬合回歸方程的方差可信度分析 Table 5 Variance reliability analysis of fitted regression equation

2.3.2 參數(shù)優(yōu)化與驗(yàn)證結(jié)果分析

1）參數(shù)優(yōu)化

尋優(yōu)后得到的因素最優(yōu)組合響應(yīng)曲面如圖8所示。其中，當(dāng)棉稈長徑比為13不變，葉輪轉(zhuǎn)速和地膜面積分別為202 r/min和1271.4 mm2時(shí)，地膜與棉稈馳豫時(shí)間差最優(yōu)為4.20 s（圖8a）；當(dāng)?shù)啬っ娣e為1413 mm2不變，葉輪轉(zhuǎn)速和棉稈長徑比分別為202 r/min和13時(shí)，地膜與棉稈最大垂向位移差最優(yōu)為22.1 cm（圖8b）；最終得到的最佳參數(shù)組合：葉輪轉(zhuǎn)速為202 r/min，地膜面積為1271.4 mm2，棉稈長徑比為13。

2）試驗(yàn)驗(yàn)證

在滿足客觀試驗(yàn)條件和便于試驗(yàn)進(jìn)行的前提下，對優(yōu)化求解后的組合數(shù)據(jù)進(jìn)行圓整取值：葉輪轉(zhuǎn)速為200 r/min，地膜面積為1256 mm2（長為40.0 mm，寬為31.4 mm），棉稈長徑比為13（棉稈直徑約為5.0 mm，長為65.0 mm），其對應(yīng)的地膜與棉稈馳豫時(shí)間差和最大垂向位移差預(yù)測值分別為4.20 s和22.1 cm。試驗(yàn)所測得地膜與棉稈馳豫時(shí)間差和最大垂向位移差的試驗(yàn)平均值分別為4.34 s和21.2 cm（表6），與預(yù)測值相對誤差分別為3.23%和4.25%，均＜5%，表明優(yōu)化條件選擇合理。最優(yōu)參數(shù)組合下，通過試驗(yàn)錄像可觀測到顯著的地膜與棉稈分離效果，故所建立的回歸擬合方程是合適的，所得到的優(yōu)化分離參數(shù)也是可行的。

表6 優(yōu)化參數(shù)下地膜與棉稈馳豫時(shí)間差和最大垂向位移差的預(yù)測值及試驗(yàn)平均值 Table 6 Predicted and experimental means of relaxation time difference and maximum vertical displacement difference under optimal condition

2.4 地膜與棉稈的沉降聚集馳豫過程影像分析

2.4.1 葉輪擾動下棉稈的向下沉降和向上聚集錄像分析

棉稈在沉降聚集過程中，其姿態(tài)是在運(yùn)動和穩(wěn)定取向間相互轉(zhuǎn)換的（圖9）。如圖9a中數(shù)字3，棉稈運(yùn)動取向?yàn)樽陨碇行妮S平行于葉輪軸（豎直）方向，且在達(dá)到平行時(shí)的聚集（沉降）速度最大；如圖9b中數(shù)字9，其在運(yùn)動取向的向上聚集過程中甚至?xí)驊T性加速度而躍遷出水面；如圖9a中數(shù)字6和圖9b中數(shù)字12，穩(wěn)定取向?yàn)樽陨碇行妮S垂直于葉輪軸（水平）方向，這與林建忠[20]在研究無碰撞作用下圓柱狀顆粒的沉降規(guī)律時(shí)得到的結(jié)論是一致的。

葉輪開始轉(zhuǎn)動后，如圖9a中數(shù)字2～4，棉稈姿態(tài)從穩(wěn)定取向的水平方向逐漸轉(zhuǎn)為運(yùn)動取向的豎直方向，并跟隨渦流向下螺旋沉降，馳豫時(shí)長為1.00～1.10 s，垂向下降深度為3～8 cm，并在此水平面內(nèi)保持并作回轉(zhuǎn)運(yùn)動，如圖9a中數(shù)字5，馳豫時(shí)長為1.20～1.30 s，垂向位移在此過程無明顯變化。如圖9b所示，棉稈向上聚集過程中運(yùn)動姿態(tài)變換與向下沉降過程相同，但總用時(shí)略短為0.90～1.00 s，向上運(yùn)動恢復(fù)至穩(wěn)定后在水平面內(nèi)一端貼葉輪軸固定，另一端繞葉輪軸作回轉(zhuǎn)圓周運(yùn)動，運(yùn)動軌跡為一組同心圓線。此外，棉稈在沉降聚集過程中不旋轉(zhuǎn)或以較小角速度繞自身中心軸旋轉(zhuǎn)，對其運(yùn)動過程影響不明顯，在后續(xù)的研究中可忽略不計(jì)。

2.4.2 葉輪停轉(zhuǎn)后地膜的向上聚集和復(fù)雜運(yùn)動錄像分析

地膜的向上聚集運(yùn)動是一個(gè)緩慢且復(fù)雜的內(nèi)螺旋周期上升過程（圖10）。如圖10a中數(shù)字2、7、9，地膜在向上聚集過程中運(yùn)動姿態(tài)是不斷變化的：運(yùn)動取向?yàn)殚L邊平行于葉輪軸方向，地膜一邊隨渦流回轉(zhuǎn)，一邊沿豎直方向向上作聚集恢復(fù)運(yùn)動，馳豫時(shí)長為0.30～0.40 s，垂向上升位移2～4 cm；如圖10a中數(shù)字1、6、10，地膜穩(wěn)定取向?yàn)殚L邊垂直于葉輪軸方向，只隨渦流在水平面內(nèi)作回轉(zhuǎn)運(yùn)動，馳豫時(shí)長為1.60～1.70 s，垂向位移無明顯變化。歸納其向上聚集運(yùn)動的1個(gè)周期：地膜先以運(yùn)動取向姿態(tài)在垂向方向上升2～4 cm，再以穩(wěn)定取向姿態(tài)在水平面內(nèi)保持作回轉(zhuǎn)運(yùn)動。重復(fù)以上周期3～6次，至其垂向位移不再發(fā)生變化或停止作水平離心擴(kuò)散運(yùn)動。此外如圖10a中數(shù)字6和7，地膜在運(yùn)動過程中以葉輪軸為中心在三維空間翻轉(zhuǎn)，運(yùn)動軌跡為一組同心螺旋線，且自身姿態(tài)與水平水面始終保持垂向最小接觸面積，這可能是因葉輪轉(zhuǎn)動產(chǎn)生的渦流內(nèi)不均勻的壓力梯度力造成的。

如圖10b中數(shù)字1～4，地膜在穩(wěn)定取向階段作回轉(zhuǎn)運(yùn)動時(shí)，自身會發(fā)生形變但并不相互吸附折疊；如圖10b中數(shù)字5～10，地膜運(yùn)動姿態(tài)重復(fù)無規(guī)律的彎曲-舒展-彎曲交互過程，可能是因地膜在此階段所受張力遠(yuǎn)大于靜電力，這一現(xiàn)象符合Arlov等[30]和Forgacs等[31]對不同撓性強(qiáng)度的纖維運(yùn)動軌跡分類下柔性纖維S型復(fù)雜運(yùn)動結(jié)論。

2.4.3 葉輪擾動下地膜與棉稈沉降聚集分離過程錄像分析

地膜與棉稈沉降聚集分離過程中其運(yùn)動軌跡及姿態(tài)變化如圖11所示。在葉輪擾動圓筒內(nèi)的水形成渦流的同時(shí)，地膜開始向下沉降，其沉降時(shí)間早于棉稈0.80～1.00 s，如圖11中數(shù)字2所示，這一現(xiàn)象表明棉稈具有遲滯沉降效應(yīng)。圖11中數(shù)字2～4為發(fā)生沉降后的棉稈快速恢復(fù)至穩(wěn)定狀態(tài)（隨渦流在葉輪軸處作回轉(zhuǎn)運(yùn)動），馳豫時(shí)長為2.00～2.40 s，垂向位移變化為3～8 cm，而此時(shí)的渦流才剛達(dá)到穩(wěn)定狀態(tài)。如圖11中數(shù)字5～10，向下沉降后的地膜首先吸附于葉輪做回轉(zhuǎn)運(yùn)動，運(yùn)動時(shí)間同葉輪轉(zhuǎn)動時(shí)間；葉輪停轉(zhuǎn)后渦流減弱消散，地膜脫離葉輪沿圓筒內(nèi)壁作螺旋運(yùn)動并再次向下沉降至最深，后向上聚集恢復(fù)至穩(wěn)定狀態(tài)，馳豫時(shí)長為7.00～7.20 s，垂向位移變化為20～30 cm。

結(jié)合圖9～圖11可知，地膜與棉稈在馳豫過程中相同點(diǎn)為1）沉降聚集過程中地膜與棉稈姿態(tài)是不斷變化的。穩(wěn)定取向均為垂直于葉輪軸方向的水平方向，運(yùn)動取向均為自身平行于葉輪軸的豎直方向；2）聚集過程中均以葉輪軸為中心隨渦流作往復(fù)回轉(zhuǎn)運(yùn)動，運(yùn)動軌跡則是以葉輪軸為中心的一組同心圓螺旋線。不同點(diǎn)為1）相較于棉稈的直接快速沉降聚集，地膜則更為間歇緩慢，間接的體現(xiàn)為其馳豫時(shí)間差和位移差，這一現(xiàn)象也對地膜與棉稈跟隨性和響應(yīng)周期差異性進(jìn)行了試驗(yàn)驗(yàn)證；2）地膜與棉稈在水平往復(fù)回轉(zhuǎn)運(yùn)動中，棉稈基本不發(fā)生旋轉(zhuǎn)或以較小角速度繞自身中心軸旋轉(zhuǎn)；地膜則因離心擴(kuò)散作用貼筒壁作圓周運(yùn)動，且在三維空間翻轉(zhuǎn)現(xiàn)象明顯。

2.4.4 葉輪擾動下地膜與棉稈特殊馳豫現(xiàn)象分析

一般試驗(yàn)條件下，棉稈馳豫時(shí)間和垂向位移的變化范圍均小于地膜，且始終位于其上方。在進(jìn)行葉輪擾動位置對其分離效果影響的單因素試驗(yàn)時(shí)，當(dāng)葉輪處于圓筒刻度線10 cm時(shí)，觀測到如圖12的現(xiàn)象：葉輪擾動作用時(shí)地膜始終位于棉稈上方，且有較明顯的分層現(xiàn)象，如圖12中數(shù)字3～6，但運(yùn)動狀態(tài)和軌跡遵循前述規(guī)律。葉輪停轉(zhuǎn)后地膜與棉稈立即恢復(fù)為水平取向的穩(wěn)定狀態(tài)。地膜與棉稈只在葉輪轉(zhuǎn)動時(shí)有分離現(xiàn)象，且持續(xù)時(shí)間較短，此時(shí)葉輪擾動位置和擾動時(shí)長對地膜與棉稈分離效果影響均變得顯著，可能是因渦流場中壓力梯度力和重力綜合作用造成的。

3 結(jié)論

1）本研究以地膜與棉稈馳豫時(shí)間差和最大垂向位移差為試驗(yàn)評價(jià)指標(biāo)，系統(tǒng)地開展了4組外部流場特性和3組地膜與棉稈物料屬性單因素試驗(yàn)，確定了影響地膜與棉稈分離效果的顯著性因素為葉輪轉(zhuǎn)速、地膜面積和棉稈長徑比。

2）在單因素試驗(yàn)結(jié)果上進(jìn)行正交試驗(yàn)，基于回歸方程分析，確定了影響地膜與棉稈馳豫時(shí)間的因素由大到小為葉輪轉(zhuǎn)速、地膜面積和棉稈長徑比，影響地膜與棉稈最大垂向位移差的因素由大到小為葉輪轉(zhuǎn)速、棉稈長徑比和地膜面積。通過尋優(yōu)獲得最佳參數(shù)組合：葉輪轉(zhuǎn)速為202 r/min，地膜面積為1271.4 mm2，棉稈長徑比為13，相對應(yīng)的地膜與棉稈馳豫時(shí)間差和最大垂向位移差預(yù)測值分別為4.20 s和22.1 cm。經(jīng)試驗(yàn)驗(yàn)證，最優(yōu)組合參數(shù)試驗(yàn)下的平均值分別為4.34 s和21.2 cm，與預(yù)測值的相對誤差分別為3.23%和4.25%，以上參數(shù)組合可滿足地膜與棉稈水洗分離方案的要求。當(dāng)?shù)啬づc棉稈物料在文中研究變化范圍內(nèi)，地膜與棉稈水洗分離效果主要受外部流場影響，在確定的試驗(yàn)場地和作業(yè)條件下，可通過合理設(shè)計(jì)機(jī)械裝置改變外部流場特性以期達(dá)到地膜與棉稈最優(yōu)分離效果。

3）通過對試驗(yàn)錄像后期處理，對地膜與棉稈水洗分離理論分析部分進(jìn)行了驗(yàn)證，并獲得了棉稈向下沉降和向上聚集過程、地膜向上聚集和復(fù)雜運(yùn)動過程、地膜與棉稈向下沉降分離過程以及在特定試驗(yàn)條件下膜在桿上6種試驗(yàn)現(xiàn)象。進(jìn)一步對比分析地膜與棉稈運(yùn)動過程中的異同點(diǎn)：地膜向下沉降時(shí)間早于棉稈0.80～1.00 s，棉稈具有遲滯沉降效應(yīng)，但地膜響應(yīng)周期更長。棉稈馳豫時(shí)長為2.00～2.40 s，垂向位移變化為3～8 cm，地膜馳豫時(shí)長為7.00～7.20 s，垂向位移變化為20～30 cm。

本研究系屬新疆棉田殘膜污染治理下新方案的嘗試與探索。通過借鑒固液多相流領(lǐng)域相關(guān)理論，定量開展流場和物料的可重復(fù)試驗(yàn)研究，并對地膜與棉稈運(yùn)動過程進(jìn)行了定性分析。研究結(jié)果為今后將殘膜回收機(jī)的農(nóng)業(yè)機(jī)械與棉田水利灌溉的農(nóng)藝相結(jié)合，如通過在農(nóng)業(yè)灌溉用水地表引水渠中裝設(shè)機(jī)械裝置，實(shí)現(xiàn)一次作業(yè)完成機(jī)收地膜的清洗、分離和回收奠定了基礎(chǔ)。