滾筒篩式膜雜分離裝置數值模擬與參數優化

摘要：滾筒篩式膜雜分離裝置存在膜雜混合物經篩分后膜中含雜率高等問題，利用計算流體力學（CFD）與離散元（DEM）耦合方法，對滾筒篩式膜雜分離裝置清選室氣流場動態分布以及裝置內物料運動情況進行數值模擬，確定工作參數范圍。選取氣流角度、進氣口風速和滾筒轉速為試驗因素，進行氣固兩相流耦合仿真試驗，分析各試驗因素對清選室氣流場以及對裝置內物料運動情況的影響，以產出物含膜率、雜中含膜率及殘膜產出量為評價指標，開展三因素三水平正交試驗，通過多目標優化的方法得到最佳工作參數組合：進氣口風速為8.5 m/s、氣流角度為8°、滾筒轉速為26 r/min，在此條件下進行驗證試驗，結果表明：產出物含膜率為89.31%、雜中含膜率為1.99%、殘膜產出量為0.032 6 kg/s，試驗結果與理論優化值相對誤差小于7%，各指標分別優于工作參數優化前0.87個、3.82個、6.57個百分點。該研究驗證基于DEM-CFD氣固兩相流耦合仿真對滾筒篩式膜雜分離裝置流場分析的可行性，為滾筒篩式膜雜分離裝置作業參數優化提供理論依據。

關鍵詞：DEM-CFD耦合；膜雜分離；數值模擬；正交試驗；響應面分析

中圖分類號：S223.5

文獻標識碼：A

文章編號：2095-5553 （2024） 04-0001-11

收稿日期：2022年4月6日" 修回日期：2022年8月11日

基金項目：國家自然科學基金項目（51805305，52175238）；山東省自然科學基金重點項目（ZR2020KE045）；山東省薯類產業技術體系（SDAIT—16—10）

第一作者：解臣碩，男，1998年生，山東日照人，碩士；研究方向為農田殘膜污染技術及裝備。E-mail： 18805485693@163.com

通訊作者：侯加林，男，1963年生，山東泰安人，博士，教授；研究方向為智能農業裝備。E-mail： jlhou@sdau.edu.cn

Numerical simulation and parameter optimization of trommel screen type membrane

and impurity separation device

Xie Chenshuo1， 2， Kang Jianming2， 3， Peng Qiangji2， 3， Lin Xiangyang4， Hou Jialin1， 2

（1. College of Mechanical and Electronic Engineering， Shandong Agricultural University， Tai’an， 271018， China；

2. Shandong Academy of Agricultural Machinery Sciences， Jinan， 250100， China; 3. Huang Huai Hai Key

Laboratory of Modern Agricultural Equipment， Ministry of Agriculture and Rural Affairs， Jinan， 250100， China；

4. Bureau of Agriculture and Rural Development， Qixia City， Shandong Province， Yantai， 265300， China）

Abstract：

Aiming at the problems of high impurity content in the membrane after sieving the membrane impurity mixture caused by the lack of theoretical support of the trommel screen type membrane and impurity separation device and the unreasonable setting of working parameters. The dynamic distribution of the air flow field in the sorting chamber of the trommel screen type membrane and impurity separation device and the movement of materials in the device are numerically simulated to determine the best combination of working parameters， by using computational fluid dynamics and discrete element coupling method. The airflow angle， air inlet wind speed and trommel screen speed were selected as test factors to conduct gas-solid two-phase flow coupling simulation test， to analyze the influence of various test factors on the air flow field in the sorting chamber and the movement of materials in the device. Taking the film content rate of the output， the film content rate of impurities and the output of residual film as evaluation indicators， a three-factor， three-level orthogonal test was carried out. The test results were subjected to multiple regression analysis， and the best combination of working parameters was obtained through the multi-objective optimization method： the air inlet wind speed was 8.5 m/s， the airflow angle was 8°， and the trommel speed was 26 r/min. The verification tests were carried out under these conditions， the results showed that the film content rate of the output was 89.31%， the film content rate of impurities was 1.99%， and the output of residual film was 0.032 6 kg/s. The relative error between the experimental results and the theoretical optimization value was less than 7%， and the indicators were 0.87， 3.82 and 6.57 percentage points better than that before the optimization of working parameters， respectively. The study verified the feasibility of analyzing the flow field of the trommel screen type membrane and impurity separation device based on DEM-CFD gas-solid two-phase flow coupling simulation， and provided a theoretical basis for the optimization of the operating parameters of the trommel screen type membrane and impurity separation device.

Keywords：

DEM-CFD coupling; membrane impurity separation; numerical simulation; orthogonal experiment；response surface analysis

0 引言

地膜覆蓋栽培技術在我國農業生產中具有重要的地位，對利用有限的水資源發展旱地農業、保障我國糧食安全起到了重要的支撐作用［1］。我國每年使用數百萬噸塑料地膜，由于地膜未及時回收且很難自然降解，對農田造成了嚴重的污染。同時，地膜屬于聚乙烯材料，是一種寶貴的可再生資源，可用來加工塑料顆粒，形成地膜“使用—回收—加工—再利用”的良性循環機制［2］。當前，國內院校及相關企業研制了多種農田殘膜回收清選機械，但回收后的殘膜中混雜有大量碎土、棉稈等雜質，難于重復利用，只能隨意堆放、掩埋、焚燒，不僅浪費資源，而且對環境造成了二次污染。

機械回收后的膜雜混合物需清雜后才能二次利用，滾筒篩式膜雜分離裝置已經廣泛用于膜雜混合物篩分工作［3， 4］，近年來，學者們對膜雜分離裝置進行了大量研究，康建明等［5］通過對膜雜混合物在滾筒篩內的運動規律進行分析，闡明了滾筒轉速與直徑的相互關系；彭強吉等［6］針對氣力式圓筒篩篩分性能波動大的問題，增設運移裝置、重置圓筒篩篩孔排布與大小，改善了雜質通過性能；石鑫等［7］試制了膜雜分離裝置整機，并對關鍵部件進行分析說明。以上研究均未對膜雜分離裝置篩分機理進行分析，工作參數以及結構參數均依靠經驗選擇，但依靠經驗設計的滾筒篩式膜雜分離裝置缺乏理論支持、工作參數設置不合理，因此迫切需要研究殘膜—土壤—棉稈混合物篩分機理，研制殘膜與雜質篩分裝備，優化篩分裝備工作參數。

計算流體力學和離散元法已廣泛應用于農業領域［8］。李洪昌等［9］采用流固耦合的方法模擬物料在風篩式清選裝置中的運動情況并與試驗結果進行對比。史嵩等［10］通過DEM-CFD氣固兩相流耦合方法對排種器進行分析，從宏觀尺度、微觀尺度對充種過程各階段臨界點進行劃分，并確定了各階段充種性能評價指標以及排種盤導種槽最佳參數組合。韓丹丹等［11］運用DEM-CFD耦合的方法，對內充氣吹式玉米排種器進行仿真試驗，確定最佳安裝結構參數，優化工作參數組合。DEM-CFD耦合常用于顆粒—流體系統的仿真，但滾筒篩式膜雜分離裝置流固耦合數值模擬目前鮮有報道。

本文以滾筒篩式膜雜分離裝置為對象，借助DEM-CFD耦合數值計算方法，分析殘膜與棉稈、土壤等雜質分離過程以及影響分離性能的主要因素，選取氣流角度、進氣口風速和滾筒轉速為試驗因素，開展三因素三水平正交試驗研究，優化滾筒篩式膜雜分離裝置工作參數組合，并通過試驗驗證最優工作參數組合下設備工作性能，以期通過耦合仿真的方法為滾筒篩式膜雜分離裝置最優工作參數選擇提供理論參考。

1 滾筒篩式膜雜分離裝置結構與工作原理

1.1 滾筒篩式膜雜分離裝置結構

滾筒篩式膜雜分離裝置結構如圖1所示，主要由風機、進料口、滾筒篩、密封罩、機架、電機、集膜箱等組成，滾筒篩由激光打孔鋼板彎折成形，直徑為1 000 mm，長為2 750 mm，帶有篩孔一段長度為2 500 mm，篩孔的直徑為30 mm，滾筒篩內部帶有螺旋葉片，集膜箱由絲網圍成，安裝在分離裝置后端。

1.2 工作原理

作業時，進氣口風速、氣流角度以及滾筒轉速通過控制器進行調節。將利用CMJD-1500型殘膜撿拾打包作業機收取、剪碎的膜雜混合物通過進料口喂入分離裝置，在滾筒篩轉動與氣流吹動的共同作用下，膜雜混合物在裝置中被打散，實現物料離散化，螺旋葉片在旋轉過程中將膜雜混合物向后運移，并進一步增大膜雜混合物的分散度，提高物料與氣流的接觸時間，充分發揮氣流的清選作用。密封罩可以減少塵土飛揚，同時保證內部空間的流場穩定性。膜雜混合物中各成分懸浮速度各不相同，密度大、懸浮速度高的土塊、棉稈等雜質，在重力場、氣流場的共同作用下，從滾筒篩前段由篩孔排出分離裝置；密度較小、懸浮速度低的地膜，在流場的作用下從滾筒篩另一側排出分離裝置進入集膜箱。

2 膜雜分離裝置數值模擬及分析

2.1 膜雜混合物運動分析及分離裝置重要工作參數

膜雜混合物喂入分離裝置后，在下落過程中受到氣流阻力Ft、重力mg的作用。建立坐標，對裝置內的物料進行運動分析，如圖2所示。

為滿足分離條件，取Sy為滾筒直徑為1 m，則殘膜水平運動距離Sx1應大于2.75 m，為減少集膜箱內的雜質含量，雜質的水平距離Sx2應小于2.75 m。綜上，將各結果代入式（4）中可得，滾筒篩式膜雜分離裝置的分離效果與氣流速度v2、氣流角度α、氣流相對物料速度與水平方向夾角γ、物料密度ρs以及物料迎風面積A有關，確定氣流速度、氣流角度為影響分離裝置分離效果的重要工作參數。

2.2 數學模型

湍流和層流狀態利用雷諾數Re進行判斷。

Re=ρuLμ

（7）

式中：

ρ——流體密度，kg/m3；

u——空氣流速，m/s；

L——特征長度，m；

μ——動力黏度，Pa·s。

選取空氣密度ρ=1.25 kg/m3，空氣流速u=3～9 m/s，取u最小值為3 m/s，特征長度L=1 m，動力黏度μ=1.85×10-5 Pa·s，代入式（1）求得Re=202 702.7，對于內部流動，通常認為雷諾數大于2 300為湍流，低于2 300為層流，則該裝置內部流場為湍流［12］，對數學模型進行數值模擬分析。

2.3 簡化模型及劃分網格

利用SolidWorks軟件對分離裝置進行結構簡化，簡化后的結構如圖3（a）所示，分別對裝置中的運動區域和靜止區域進行網格劃分，整體采用六面體網格進行劃分，滾筒篩以及滾筒篩與密封罩之間的流體區域網格邊長為10 mm，對于篩孔等物理量變化劇烈的位置，進行網格加密處理，網格邊長為2 mm，篩孔處通過interface進行計算域流通。

將繪制好的動網格區域和靜網格區域進行合并，檢查合并后的網格數量與質量，合并后的網格如圖3（b）所示，網格總數為1 939 410個，網格質量Skewness為0.36，網格數量與質量符合要求，導出.msh文件。

2.4 流體相數值模擬

2.4.1 前處理設置

利用Fluent軟件對滾筒篩式膜雜分離裝置內部流場區域進行數值模擬。將劃分好的網格導入Fluent軟件中，選擇Standard k-epsilon湍流模型，設置各邊界條件，在相互接觸的篩孔邊界處通過interface對計算結果數據進行插值傳遞，對Fluent軟件初始化后進行仿真迭代，當標度殘差曲線圖中的連續性曲線降到合適的收斂位置時，仿真完成。

2.4.2 流場分析

當分離裝置進氣口風速為9 m/s、滾筒轉速為26 r/min、氣流角度為0°時滾筒內部流場的速度矢量圖、速度矢量局部放大圖和壓力云圖如圖4和圖5所示，通過手持式熱敏風速儀，對滾筒內部多個位置的氣流速度進行測量，與仿真結果進行對比。結果表明，仿真的流場環境與實際流場環境相符，證明流場仿真結果具有較高的準確性。

如圖4（a）所示，建立直角坐標系，x軸、y軸正方向如圖所示；從中心處開始，將徑向半徑為200 mm的圓柱劃分為區域Ⅰ，從區域Ⅰ邊緣至滾筒壁劃分為區域Ⅱ，滾筒篩壁到密封罩壁之間劃分為區域Ⅲ；將進氣口一端視為滾筒前端，出氣口一端視為滾筒后端。

從圖4（a）速度矢量圖可知，沿x軸方向，區域Ⅰ的氣流速度呈現衰減的趨勢，由8.25～9.71 m/s逐漸降低到4.86～6.31 m/s；區域Ⅱ處由于氣流擴展，使得滾筒后端的氣流速度高于滾筒前端氣流速度，且滾筒前端氣流呈現不規則的運動趨勢，后端氣流運動方向朝向x軸正方向；區域Ⅲ的空氣在壓力差的作用下，滾筒前端的空氣通過篩孔由區域Ⅲ流入區域Ⅱ，滾筒后端的空氣通過篩孔由區域Ⅱ流入區域Ⅲ。沿y軸方向，由于氣流角度為0°，且滾筒可近似看成薄壁對氣流影響較小，速度矢量沿x軸兩側呈對稱分布。

如圖4（a）、圖5所示，滾筒前端區域Ⅰ處氣流速度要大于區域Ⅱ以及區域Ⅲ的氣流速度，由伯努利方程可得，滾筒前端區域Ⅰ處壓強要小于區域Ⅲ處壓強，由局部放大圖M可知，滾筒前端篩孔處空氣在壓力差的作用下由滾筒外側流入滾筒內側，氣流速度為1.46～2.92 m/s；滾筒后端區域Ⅰ處氣流速度要大于區域Ⅱ以及區域Ⅲ的氣流速度，且滾筒后端氣流與滾筒前端氣流有較大的速度差，故滾筒后端區域Ⅰ的壓強要大于區域Ⅲ的壓強，由局部放大圖N可知，滾筒后端篩孔處空氣在壓力差的作用下由滾筒內側流入滾筒外側，氣流速度為0.491～1.46 m/s。

對滾筒式膜雜分離裝置內部流場分析可知，由于裝置內部流場復雜，各篩孔處氣流方向、大小均不相同，故進氣口風速在滿足殘膜清選條件下，應盡可能取較大值，避免殘膜運動到篩孔處受到流場提供的吸附力，造成雜中含膜率高的問題。故根據懸浮試驗測量以及裝置內部流場分析，裝置工作過程進氣口風速應取3～9 m/s。

2.5 DEM-CFD耦合仿真分析

DEM-CFD耦合方法基本思路：通過Fluent軟件求解流場，使用EDEM軟件計算顆粒系統的運動受力情況，二者通過耦合接口進行質量、動量和能量等的傳遞，實現耦合［11］。文中利用Fluent軟件進行流體仿真模擬，對分離裝置清選室氣流場動態分布情況進行觀察，導出流場數據文件，通過EDEM軟件的API二次開發功能，將得到的流場數據加載到EDEM軟件中，得到分離裝置內部流場分布，如圖6所示。觀察顆粒受力及運動情況，并監測滾筒下端和集膜箱內物料質量。

測量殘膜、棉稈等物料的外形參數，在EDEM軟件中采用球顆粒聚合方法建立物料離散模型，其中殘膜顆粒直徑為1 mm，殘膜模型尺寸（長×寬）為60 mm×30 mm；土壤顆粒直徑為1 mm，土壤模型由兩個相同的球形顆粒聚合而成；棉稈顆粒直徑為9 mm，棉稈模型長度為60 mm，物料離散模型如圖7所示。

通過排水法測量殘膜以及棉稈的密度；借助萬能試驗機、質構儀分別對棉稈、殘膜進行拉伸試驗，通過伸長量與截面收縮量之間的比值，確定物料的泊松比；自制斜面儀，通過測量兩物料碰撞后彈跳的距離，計算不同物料之間的接觸參數，其他物料的物理特性參數通過查閱文獻獲得［13-15］，各類參數如表1所示，在軟件中進行相應設置。殘膜在裝置內部體積分數較小，為減小計算量提高仿真效率、忽略殘膜的厚度對分離裝置內部流場以及殘膜運動情況的影響，適當增加殘膜顆粒直徑的同時減小殘膜顆粒密度，以保持殘膜質量不變。

將裝置模型導入到EDEM軟件中，設置滾筒篩、螺旋葉片等工作部件的運動。增加顆粒工廠、膜雜混合物中各物料質量比例為土壤、碎土塊等43%，棉稈等21%，廢舊農膜36%，按照測定的膜雜混合物中各成分質量比例設置物料的投入量。為保證裝置內部具有充足的物料且喂料口不堵塞，根據前期試驗，設置喂入量為175 kg/h，滾筒篩式膜雜分離裝置工作過程模擬如圖8所示。

2.6 單因素試驗

理想狀態下，膜雜混合物投入分離裝置內部，殘膜在氣流帶動下飄向裝置后端的集膜箱，只有少量殘膜滯留在滾筒內部，故裝置內部受力主要由棉稈與土壤顆粒提供；對分離裝置滾筒篩受力變化進行分析，通過受力變化的大小衡量裝置內部物料運動的劇烈程度，若裝置內部受力異常增大，則表明大量的殘膜滯留在裝置內部流場中，紊亂裝置內部流場穩定，故該工作參數組合不利于膜雜混合物的分離。

完成仿真，導出滾筒篩在各時刻受力的最大值，在其他工作參數相同的情況下改變氣流角度［16-18］，如圖9（a）、圖9（c）所示，在氣流角度為0°和10°的情況下，滾筒受到的碰撞力較大且波動劇烈，當氣流角度為0°時最大受力為10.28 N，氣流角度為10°時最大受力為7.02 N；如圖9（b）所示，在氣流角度為5°的情況下，滾筒受到的碰撞力較小且波動較為規律，最大受力為4.15 N。綜上，氣流角度的變化對分離裝置分離過程具有一定的影響，確定氣流角度的最優值位于0°～10°之間。

在其他工作參數相同的情況下改變滾筒轉速，如圖10（a）、圖10（c）所示，在滾筒轉速為16 r/min和36 r/min的情況下，滾筒受到的碰撞力較大且波動劇烈。當滾筒轉速為16 r/min時滾筒最大受力為8.35 N，滾筒轉速為36 r/min時滾筒最大受力為7.32 N；如圖10（b）所示，在滾筒轉速為26 r/min的情況下，滾筒受到的碰撞力較小且波動較為規律，最大受力為4.68 N。

當滾筒轉速較低時，裝置內的物料無法及時運出，滯留在裝置內部的物料對裝置反復進行碰撞，裝置受到的碰撞力較大且波動劇烈；當滾筒轉速較高時，拋灑的物料與高速旋轉的滾筒發生碰撞，彈起的棉稈、土壤等物料混入吹起的殘膜中，膜中含雜率升高。綜上，滾筒轉速的變化對分離裝置分離過程具有一定的影響，確定滾筒轉速最優值位于16～36 r/min之間。

通過單因素試驗，滾筒轉速、氣流角度的改變影響分離裝置內部流場穩定以及分離效率，確定為影響分離裝置工作的重要因素。

3 參數優化試驗

3.1 試驗參數與方法

進行實機試驗，將不同工作參數組合下測定的產出物含膜率Y1、雜中含膜率Y2、殘膜產出量Y3作為滾筒篩式膜雜分離裝置工作性能的評價指標。選取氣流角度、進氣口風速、滾筒轉速為試驗因素，根據單因素試驗以及懸浮試驗測得物料懸浮速度，確定三因素參數范圍：進氣口風速為3～9 m/s、氣流角度為0°～10°、滾筒轉速為16～36 r/min。進氣口風速過大，殘膜產出量提高，但較大的風力將棉稈等吹到集膜箱內導致產出物含膜率較高，進氣口風速過小殘膜不能被順利吹出分離裝置，導致殘膜產出量下降以及雜中含膜率上升；氣流角度過小殘膜漂浮距離近，雜中含膜率高，氣流角度過大漂浮的殘膜容易與滾筒上側發生碰撞，于螺旋葉片纏繞，使得殘膜產出量下降；滾筒轉速過大部分未及時分散的膜雜混合物被直接輸送到后方的集膜箱內，導致產出物含膜率上升，滾筒轉速過小膜雜混合物打散效果差。為了獲得三個因素的最佳參數組合，借助Design-Expert軟件，選擇Box-Behnken試驗設計原理，對氣流角度X1、進氣口風速X2、滾筒轉速X3進行響應面試驗。試驗因素與水平如表2所示。

試驗結束，測量分離裝置滾筒后方集膜箱處以及滾筒篩下方的土壤、殘膜、棉稈質量。滾筒篩式膜雜分離裝置各評價指標計算公式如式（8）、式（9）所示。

Y1=m2m1×100%

（8）

式中：

Y1——產出物含膜率，%；

m1——集膜箱內物料的總質量，kg；

m2——集膜箱內殘膜的質量，kg。

Y2=m4m3×100%

（9）

式中：

Y2——雜中含膜率，%；

m3——滾筒篩下方物料的總質量，kg；

m4——滾筒篩下方殘膜的質量，kg。

在保證雜中含膜率最低的條件下，為保證較高的殘膜產出量，定義Y3為仿真單位時間內集膜箱中收集的殘膜總質量，單位為kg/s，用來衡量地膜產出量的高低。

3.2 正交試驗與結果分析

3.2.1 試驗結果

試驗設計方案及響應值如表3所示，其中X1、X2、X3為各因素編碼值。

3.2.2 回歸模型建立與顯著性檢驗

運用Design-Expert軟件對表3數據進行多元回歸擬合分析［19］，得到各因素與產出物含膜率Y1、雜中含膜率Y2、殘膜產出量Y3的回歸方程，如式（10）～式（12）所示，回歸方程的方差分析結果如表4所示。

從表4可以看出，響應面模型中產出物含膜率Y1、雜中含膜率Y2、殘膜產出量Y3的響應面模型的P值均小于0.05，說明回歸模型顯著；其決定系數R2值分別為0.893 4、0.967 5、0.981 9，表明這3個模型可以解釋89%以上的評價指標［20］。因此該模型可以優化滾筒篩式膜雜分離裝置的工作參數。

產出物含膜率Y1模型中X1、X12兩個回歸項對模型影響極顯著（P＜0.01），X1X3對模型影響顯著（P＜0.05）；雜中含膜率Y2模型中X1、X3、兩個回歸項對模型影響極顯著（P＜0.01），X12對模型影響顯著（P＜0.05）；殘膜產出量Y3模型中X1、X2、X3、X1X3、X2X3、X12、X22、X32八個回歸項對模型影響極顯著（P＜0.01）。

各因素對滾筒篩式膜雜分離裝置分離效果的影響可通過分析貢獻率K值得出，各因素對產出物含膜率影響順序為：進氣口風速X1＞氣流角度X2＞滾筒轉速X3；各因素對雜中含膜率影響順序為：滾筒轉速X3＞進氣口風速X1＞氣流角度X2；各因素對殘膜產出量影響順序為：進氣口風速X1＞滾筒轉速X3＞氣流角度X2［21］。

3.2.3 交互因素對性能影響規律分析

根據上述回歸方程分析結果，利用Design-Expert軟件繪制響應面圖，根據響應面圖分析進氣口風速、氣流角度、滾筒轉速交互因素對響應值的影響［22］。

1） 交互因素對產出物含膜率的影響規律分析。進氣口風速、氣流角度、滾筒轉速交互因素對產出物含膜率影響的響應面圖如圖11所示。圖11（a）為滾筒轉速位于中心水平（26 r/min）時，進氣口風速與氣流角度對產出物含膜率交互作用的響應面圖。從圖11（a）可以看出，在同一氣流角度下產出物含膜率隨著進氣口風速的增大先增大后減小；當進氣口風速較小時，產出物含膜率隨著氣流角度的增加而增大，當進氣口風速較大時，產出物含膜率隨著氣流角度的增大先減小后增大。

圖11（b）為氣流角位于中心水平（5°）時，進氣口風速與滾筒轉速對產出物含膜率交互作用的響應面圖。從圖11（b）可以看出，當滾筒轉速較高時，產出物含膜率隨著進氣口風速的增大先增大后減小，當滾筒轉速較低時，產出物含膜率隨著進氣口風速的增加而增加；當進氣口風速較低時，產出物含膜率隨著滾筒轉速的增大而增大，當進氣口風速較高時，產出物含膜率隨著滾筒轉速的增大而減小。從圖11（c）可以看出，氣流角度與滾筒轉速交叉因素對指標影響不顯著。總體影響趨勢為進氣口風速、氣流角度越大、滾筒轉速適中，則產出物含膜率越高。

2）" 交互因素對雜中含膜率的影響規律分析。各因素對雜中含膜率影響的響應面圖如圖12所示。

圖12（a）為滾筒轉速位于中心水平（26 r/min）時，進氣口風速與氣流角度對產出物含膜率交互作用的響應面圖。從圖12（a）可以看出，在同一氣流角度下雜中含膜率隨著進氣口風速的增大先增大后減小；當進氣口風速較小時，產出物含膜率隨著氣流角度的增加而增大，當進氣口風速較大時，雜中含膜率隨著氣流角度的增大而減小。圖12（b）為氣流角位于中心水平（5°）時，進氣口風速與滾筒轉速對雜中含膜率交互作用的響應面圖。從圖12（b）可以看出，在同一滾筒轉速下進氣口風速對雜中含膜率的影響不顯著；在同一進氣口風速下雜中含膜率隨著滾筒轉速的增加而減小。圖12（c）為進氣口風速位于中心水平（6 m/s）時，氣流角度與滾筒轉速對雜中含膜率交互作用的相應面圖。從圖12（c）可以看出，在同一滾筒轉速下，氣流角度對雜中含膜率的影響不顯著；在同一氣流角度下雜中含膜率隨著滾筒轉速的增大而減小。總體影響趨勢為進氣口風速、氣流角度、滾筒轉速越大，則雜中含膜越低。

3）" 交互因素對殘膜產出量的影響規律分析。各因素對殘膜產出量影響的響應面圖如圖13所示。圖13（a）為滾筒轉速位于中心水平（26 r/min）時，進氣口風速與氣流角度對殘膜產出量交互作用的響應面圖。從圖13（a）可以看出，在同一氣流角度下殘膜產出量隨著進氣口風速的增大先增大后減小；在同一進氣口風速下殘膜產出量隨著氣流角度的增大先增大后減小。圖13（b）為氣流角位于中心水平（5°）時，進氣口風速與滾筒轉速對殘膜產出量交互作用的響應面圖。從圖13（b）可以看出，在同一滾筒轉速下殘膜產出量隨著進氣口風速的增大先增大后減小；在同一進氣口風速下殘膜產出量隨著滾筒轉速的增大先增大后減小。圖13（c）為進氣口風速位于中心水平（6 m/s）時，氣流角度與滾筒轉速對殘膜產出量交互作用的相應面圖。從圖13（c）可以看出，在同一滾筒轉速下殘膜產出量隨著氣流角度的增大先增大后減小，在同一氣流角度下殘膜產出量隨著滾筒轉速的增大先增大后減小。總體影響趨勢為進氣口風速、氣流角度、滾筒轉速越大，則殘膜產出量越高。

4 參數優化與驗證試驗

4.1 參數優化

為提升滾筒篩式膜雜分離裝置的作業性能，需滿足產出物含膜率高、雜中含膜率低、殘膜產出量高，綜合考慮各因素對響應值的影響關系，通過多目標優化，獲取最佳的工作參數組合。

本研究針對滾筒篩式膜雜分離裝置工作參數優化，為滿足產出物含膜率高、雜中含膜率低、殘膜產出量高的作業要求，根據滾筒篩式膜雜分離裝置的實際工作條件、作業性能要求，選擇優化約束條件為

maxY1（X1，X2，X3）

minY2（X1，X2，X3）

maxY3（X1，X2，X3）

Yi＞0，i=1，2，3

-1≤Xj≤1，j=1，2，3

（13）

利用Design-Expert軟件對各參數進行優化求解。當進氣口風速為8.54 m/s、氣流角度為8.01°、滾筒轉速為26.36 r/min時，此時產出物含膜率為87.67%、雜中含膜率為2.1%、殘膜產出量為0.034 8 kg/s。

4.2 試驗驗證

為了驗證模型預測的準確性，采用上述參數在山東省農業機械科學研究院試制工廠進行3次重復試驗。

對理論最優工作參數進行圓整，將進氣口風速設置為8.5 m/s、氣流角度為8°、滾筒轉速為26 r/min，在此優化方案下進行整機試驗，利用手持熱敏式風速儀（測量范圍：0～30 m/s，誤差：±1%）確定進氣口風速，三量數顯傾斜角儀（測量范圍：0°～90°）確定進風口管道角度，UT372型高精度非接觸式轉速儀（測量范圍：0～99 999 r/min）確定滾筒轉速，分離物利用電子秤確定其中殘膜及雜物的質量，計算評價指標值，結果如表5所示。

從表5可知，試驗值與理論優化值接近，相對誤差均小于7%，因此參數優化模型可靠，驗證了基于DEM-CFD氣固兩相流耦合仿真對滾筒篩式膜雜分離裝置流場分析的可行性。滾筒篩式膜雜分離裝置優化工作參數前，產出物含膜量為88.54%、雜中含膜率為2.13%、殘膜產出量為0.031 4 kg/s；采用優化工作參數組合，即進氣口風速為8.5 m/s、氣流角度為8°、滾筒轉速為26 r/min，此時產出物含膜率為89.31%、雜中含膜率為1.99%、殘膜產出量為0.032 6 kg/s，與工作參數優化前各指標相比，產出物含膜量、殘膜產出量分別提高0.87個、3.82個百分點，雜中含膜率降低6.57個百分點，為滾筒篩式膜雜分離裝置工作參數優化提供參考。

5 結論

1）" 以現有滾筒篩式膜雜分離裝置為基礎，借助Fluent軟件對裝置內部流場進行分析，利用DEM-CFD耦合仿真的方法，對分離裝置工作過程進行仿真分析，根據仿真得到的結果，通過正交試驗優化分離裝置工作參數，驗證了DEM-CFD耦合仿真對滾筒篩式膜雜分離裝置分析的可行性。

2）" 通過響應曲面分析，各因素對產出物含膜率影響順序為：進氣口風速＞氣流角度＞滾筒轉速；各因素對雜中含膜率影響順序為：滾筒轉速＞進氣口風速＞氣流角度；各因素對殘膜產出量影響順序為：進氣口風速＞滾筒轉速＞氣流角度。

3）" 滾筒篩式膜雜分離裝置最優工作參數組合進氣口風速為8.5 m/s、氣流角度為8°、滾筒轉速為26 r/min，此時產出物含膜率為89.31%、雜中含膜率為1.99%、殘膜產出量為0.032 6 kg/s，與工作參數優化前各指標相比，產出物含膜量、殘膜產出量分別提高0.87個、3.82個百分點，雜中含膜率降低6.57個百分點，為滾筒篩式膜雜分離裝置工作參數優化提供參考。

參 考 文 獻

［1］

張德奇， 廖允成， 賈志寬. 旱區地膜覆蓋技術的研究進展及發展前景［J］. 干旱地區農業研究， 2005， 23（1）： 208-213.

Zhang Deqi， Liao Yuncheng， Jia Zhikuan. Research advances and prospects of film mulching in arid and semi-arid areas ［J］. Agricultural Research in Arid Areas， 2005， 23（1）： 208-213.

［2］

張恒， 康建明， 張國海， 等. 黃淮海地區農膜污染現狀及對策分析［J］. 中國農機化學報， 2019， 40（1）： 156-161.

Zhang Heng， Kang Jianming， Zhang Guohai， et al. Analysis of current situation and countermeasures of agricultural film pollution in Huang-Huai-Hai area ［J］. Journal of Chinese Agricultural Mechanization， 2019， 40（1）： 156-161.

［3］

唐紅俠， 趙由才. 滾筒篩篩分生活垃圾的理論研究［J］. 環境工程學報， 2007， 1（12）： 124-127.

Tang Hongxia， Zhao Youcai. Research on theories of the trommel screen separating municipal solid waste ［J］. Chinese Journal of Environmental Engineering， 2007， 1（12）： 124-127.

［4］

石鑫， 牛長河， 喬園園， 等. 塑料垃圾分選技術在廢舊地膜與雜質分離中的應用研究進展［J］. 農業工程學報， 2016， 32（S2）： 22-31.

Shi Xin， Niu Changhe， Qiao Yuanyuan， et al. Application of plastic trash sorting technology in separating waste plastic mulch films from impurities ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2016， 32（S2）： 22-31.

［5］

康建明， 張恒， 張國海， 等. 殘膜物料空氣動力學特性與膜雜分離裝置試驗［J］. 中國農機化學報， 2020， 41（1）： 167-172.

Kang Jianming， Zhang Heng， Zhang Guohai， et al. Aerodynamic characteristics of residual film materials and test of membrane separation device ［J］. Journal of Chinese Agricultural Mechanization， 2020， 41（1）： 167-172.

［6］

彭強吉， 李成松， 康建明， 等. 氣力式圓筒篩膜雜分離機改進設計與試驗［J］. 農業機械學報， 2020， 51（8）： 126-135.

Peng Qiangji， Li Chengsong， Kang Jianming， et al. Improved design and test on pneumatic cylinder sieve film hybrid separator ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2020， 51（8）： 126-135.

［7］

石鑫， 牛長河， 王學農， 等. 滾筒篩式廢舊地膜與雜質風選裝置設計［J］. 農業工程學報， 2017， 33（18）： 19-26.

Shi Xin， Niu Changhe， Wang Xuenong， et al. Design of roller sieve waste plastic film and trash winnowing machine ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2017， 2017， 33（18）： 19-26.

［8］

蔣恩臣， 孫占峰， 潘志洋， 等. 基于CFD-DEM的收獲機分離室內谷物運動模擬與試驗［J］. 農業機械學報， 2014， 45（4）： 117-122.

Jiang Enchen， Sun Zhanfeng， Pan Zhiyang， et al. Numerical simulation based on CFD-DEM and experiment of grain moving laws in inertia separation chamber ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2014， 45（4）： 117-122.

［9］

李洪昌， 李耀明， 唐忠， 等. 風篩式清選裝置振動篩上物料運動CFD-DEM數值模擬［J］. 農業機械學報， 2012， 43（2）： 79-84.

Li Hongchang， Li Yaoming， Tang Zhong， et al. Numerical simulation of material motion on vibrating screen of air-and-screen cleaning device based on CFD-DEM ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2012， 43（2）： 79-84.

［10］

史嵩， 劉虎， 位國建， 等. 基于DEM-CFD的驅導輔助充種氣吸式排種器優化與試驗［J］. 農業機械學報， 2020， 51（5）： 54-66.

Shi Song， Liu Hu， Wei Guojian， et al. Optimization and experiment of pneumatic seed metering device with guided assistant filling based on EDEM-CFD ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2020， 51（5）： 54-66.

［11］

韓丹丹， 張東興， 楊麗， 等. 內充氣吹式玉米排種器工作性能EDEM-CFD模擬與試驗［J］. 農業工程學報， 2017， 33（13）： 23-31.

Han Dandan， Zhang Dongxing， Yang Li， et al. EDEM-CFD simulation and experiment of working performance of inside-filling air-blowing seed metering device in maize ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2017， 33（13）： 23-31.

［12］

胡坤， 胡婷婷， 馬海峰， 等. ANSYS CFD入門指南——計算流體力學基礎及應用［M］. 北京： 機械工業出版社， 2018.

［13］

Martina C L， Bouvarda D， Shimab S. Study of particle rearrangement during powder compaction by the discrete element method ［J］. Journal of the Mechanics and Physics of Solids， 2003， 51（4）： 667-693.

［14］

Cleary P W. Predicting charge motion， power draw， segregation and wear in ball mills using discrete element methods ［J］. Minerals Engineering， 1998， 11（11）： 1061-1080.

［15］

龔明. 材料物料特性參數標定［R］. 北京： EDEM用戶大會報告， 2013.

［16］

王磊， 廖慶喜， 廖宜濤， 等. 氣送式排肥系統分配裝置結構型式對排肥性能的影響［J］. 農業工程學報， 2021， 37（7）： 24-34.

Wang Lei， Liao Qingxi， Liao Yitao， et al. Effects of distributor structure types on fertilizing performance in an air-assisted applicator ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2021， 37（7）： 24-34.

［17］

Landry H， Thirion F， Lague C， et al. Numerical modeling of the flow of organic fertilizers in land application equipment ［J］. Computers amp; Electronics in Agriculture， 2006， 51（1/2）： 35-53.

［18］

Sugirbay A M， Zhao J， Nukeshev S O， et al. Determination of pin-roller parameters and evaluation of the uniformity of granular fertilizer application metering devices in precision farming ［J］. Computers amp; Electronics in Agriculture， 2020， 179（12）： 1-11.

［19］

康建明， 陳學庚， 溫浩軍， 等. 基于響應面法的梳齒式采棉機采收臺優化設計［J］. 農業機械學報， 2013， 44（S2）： 57-61.

Kang Jianming， Chen Xuegeng， Wen Haojun， et al. Optimization of comb-type cotton picker device based on response surface methodology ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2013， 44（S2）： 57-61.

［20］

吳騰， 胡良龍， 王公仆， 等. 步行式甘薯碎蔓還田機的設計與試驗［J］. 農業工程學報， 2017， 33（16）： 8-17.

Wu Teng， Hu Lianglong， Wang Gongpu， et al. Design and test of walking sweet potato（Ipomoea batatas） vines crushing and returning machine ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2017， 33（16）： 8-17.

［21］

嚴偉， 胡志超， 吳努， 等. 鏟篩式殘膜回收機輸膜機構參數優化與試驗［J］. 農業工程學報， 2017， 33（1）： 17-24.

Yan Wei， Hu Zhichao， Wu Nu， et al. Parameter optimization and experiment for plastic film transport mechanism of shovel screen type plastic film residue collector ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2017， 33（1）： 17-24.

［22］

孫岳， 簡建明， 田玉泰， 等. 殘膜回收機旋轉式起膜裝置起膜機理分析與試驗［J］. 農業機械學報， 2018， 49（S1）： 304-310.

Sun Yue， Jian Jianming， Tian Yutai， et al. Analysis and experiment of filming mechanism of rotary film-lifting device of residual film recycling machine ［J］. Transaction of the Chinese Society for Agricultural Machinery， 2018， 49（S1）： 304-310.