摘 要： 針對綠豆收獲清選過程中清選損失率高、籽粒含雜率高的問題，研發(fā)一種雙斜面風(fēng)篩式綠豆清選裝置。運用流固耦合仿真模擬綠豆清選過程，并開展了正交試驗。同時運用二次回歸方程分析和響應(yīng)面分析等方法對籽粒含雜率、清選損失率的影響進行探究。試驗結(jié)果表明，最優(yōu)參數(shù)組合為入風(fēng)口風(fēng)速5.7 m/s、清選篩振動頻率4.71 Hz、上篩橫向振幅7.99 mm、下篩篩面傾角0.68°，此條件下籽粒含雜率2.12%、清選損失率2.51%，符合性能指標(biāo)要求，實物樣機可用于后期田間試驗。

關(guān)鍵詞：綠豆；清選裝置；風(fēng)篩式振動篩；響應(yīng)面分析；綠豆收獲機

中圖分類號：S225 文獻(xiàn)標(biāo)識碼：A 文章編號：2095-1795（2024）06-0100-06

DOI：10.19998/j.cnki.2095-1795.2024.06.017

0 引言

清選裝置是綠豆聯(lián)合收獲機的關(guān)鍵組成部分之一，其性能優(yōu)劣直接影響整機的工作性能[1]。目前關(guān)于綠豆收獲及其脫粒清選機械的研究較少。侯俊銘等[2] 設(shè)計了一種雙層傾斜振動風(fēng)篩式蓖麻清選裝置，通過分析上篩面篩孔排列型式、篩孔直徑、篩面傾角對篩分效率和損失率的影響，表明雙層傾斜振動篩可以實現(xiàn)篩分效率和籽粒透篩率的有效提高。袁子翔等[3] 采用CFD-DEM 耦合方法，對胡麻脫粒物料在胡麻聯(lián)合收獲機風(fēng)篩式清選裝置中的清選效果進行研究，闡明了胡麻脫粒物料在振動篩上實際的篩分效果。

針對綠豆收獲機清選過程中存在清選損失率高、含雜率高等突出問題，結(jié)合仿真軟件模擬試驗，設(shè)計一種低含雜率、低損失率雙斜面風(fēng)篩式綠豆清選裝置，并測得該裝置的最佳參數(shù)組合。

1 清選裝置結(jié)構(gòu)設(shè)計

1.1 結(jié)構(gòu)組成

綠豆清選裝置結(jié)構(gòu)如圖1 所示，主要由風(fēng)機、入風(fēng)口、入料口、上篩、清選室、下篩、機架、集料箱、下篩圓盤式曲柄?連桿?擺桿傳動機構(gòu)、上篩圓盤式曲柄?連桿?擺桿傳動機構(gòu)、上篩驅(qū)動軸、下篩驅(qū)動軸、一級同步帶傳動機構(gòu)、二級同步帶傳動機構(gòu)和電機等組成。

1.2 工作原理

工作時，電機和風(fēng)機同時啟動。電機通過一級同步帶傳動機構(gòu)帶動上篩驅(qū)動軸兩端的圓盤式曲柄轉(zhuǎn)動，進而通過上篩圓盤式曲柄?連桿?擺桿傳動機構(gòu)帶動上篩往復(fù)運動。電機帶動上篩驅(qū)動軸轉(zhuǎn)動的同時，上篩驅(qū)動軸通過二級同步帶傳動機構(gòu)帶動下篩驅(qū)動軸轉(zhuǎn)動，下篩驅(qū)動軸帶動軸端兩側(cè)的上篩圓盤式曲柄運動，進而帶動下篩圓盤式曲柄?連桿?擺桿傳動機構(gòu)運動，令下篩開始往復(fù)運動。同時通過入風(fēng)口將勻速氣流輸送到清選室，自入料口落下的待清選物料在上下篩和風(fēng)力的綜合作用下得到充分清選，最終綠豆籽粒和雜余收集到集料箱的各自倉室中，通過各自出料口可排出倉外，完成脫出物清選。

2 建立脫出物離散元模型與清選裝置模型

2.1 綠豆脫出物離散元模型建立

綠豆清選脫出物主要成分有綠豆籽粒、短莖稈和綠豆豆莢。為提高耦合仿真綠豆清選過程的準(zhǔn)確性，先通過Solidworks 軟件對綠豆清選脫出物各組分進行三維建模。再將其導(dǎo)入EDEM 軟件中，對脫出物的各個模型用不同半徑的球形顆粒進行自動填充，直至接近于實物。其中， 綠豆籽粒、短莖稈和豆莢模型[ 長×寬×厚（高）] 分別為5.16 mm×4.03 mm×4.03 mm、12.35 mm×3.71 mm×3.71 mm 和64.56 mm×4.53 mm×2.27 mm，如圖2 所示[4-5]。

綠豆脫出物仿真參數(shù)的設(shè)定直接影響仿真分析結(jié)果的準(zhǔn)確性。綠豆脫出物的各項物理特性參數(shù)和材料間的接觸系數(shù)如表1 和表2 所示。

2.2 清選裝置模型建立

采用Solidworks2022 進行三維建模，采用1∶1 等比例簡化清選裝置物理模型，如圖3 所示。清選室模型（長×寬×高）為1 300 mm×460 mm×800 mm；進風(fēng)口和出風(fēng)口（長×寬）均為460 mm×450 mm。

采用EDEM2018 與Fluent19.2 進行耦合仿真，具體參數(shù)設(shè)置如表3 和表4 所示。

在EDEM 對試驗數(shù)據(jù)后處理中，分籽粒集料倉和雜余集料倉兩個統(tǒng)計區(qū)域。其中，籽粒含雜率和清選損失率分別由式（1）和式（2）計算[6]。

式中 Y1——籽粒含雜率

M2、M3——籽粒集料倉豆莢、短莖稈的總質(zhì)量

M0——籽粒集料倉物料總質(zhì)量

Y2——清選損失率

M1——籽粒集料區(qū)綠豆籽粒總質(zhì)量

MZ1——雜余集料倉綠豆籽?？傎|(zhì)量

3 EDEM-Fluent 流固耦合仿真試驗分析

3.1 試驗設(shè)計

為尋求上述裝置清選作業(yè)的最佳組合參數(shù)，依次通過單因素試驗和篩選試驗確定了正交試驗中各試驗因素的取值范圍和各試驗因素中對籽粒含雜率影響最高的4 個試驗因素，并確定余下兩個試驗因素：下篩橫向振幅、上篩篩面傾角分別取6 mm 和2°。用x1、x2、x3 和x4 分別表示選出的4 個試驗因素（入風(fēng)口風(fēng)速、振動頻率、上篩橫向振幅及下篩篩面傾角），用Y1、Y2 表示清選性能的試驗指標(biāo)（籽粒含雜率、清選損失率），開展4 因素3 水平正交試驗，并依據(jù)試驗結(jié)果進行相關(guān)分析，各試驗因素水平選取如表5 所示，正交試驗與結(jié)果如表6 所示[7]。

3.2 回歸模型與方差分析

根據(jù)表6 試驗結(jié)果，應(yīng)用Design-Expert 軟件對籽粒含雜率Y1 與清選損失率Y2 進行回歸模型方差分析，分別剔除不顯著項可得到Y(jié)1 與Y2 的二次回歸模型。

Y1=2.37?0.388x1?0.217x2?0.096x3?0.063x4+0.013x1x2+0.005x1 x3 ?0.003x1 x4 +0.003x2 x3 +0.05x2 x4 ?0.005x3 x4?0.113x12?0.028x22+0.228x32+0.268x42 （3）

Y2=2.17+0.408x1+0.263x2+0.057x3+0.144x4?0.093x1x2+0.003x1x3?0.008x1x4+0.003x2x3+0.03x3x4+0.262x12+0.146x22+0.014x32?0.058x42 （4）

由表7 可知，籽粒含雜率Y1 回歸模型Plt;0.000 1，表明回歸模型顯著性極高，說明二次回歸模型擬合效果好。其中，籽粒含雜率回歸模型一次項x1、x2、x3 和x4 對籽粒含雜率影響均極顯著；交互項x2x4 對籽粒含雜率影響顯著；二次項x12、x32 和x42 對籽粒含雜率影響均極顯著。根據(jù)回歸模型各因素影響程度的大小，得到各因素影響程度主次順序為x1、x2、x3 和x4。

由表8 可知，清選損失率Y2 回歸模型中失擬項Plt;0.000 1，表明回歸模型顯著性極高，說明二次回歸模型擬合效果好。其中，清選損失率回歸模型一次項x1、x2 和x4 對清選損失率影響均極顯著，x3 對清選損失率影響達(dá)很顯著；交互項x1x2 對清選損失率影響顯著；二次項x12、x22 對清選損失率影響極顯著，x42 對清選損失率影響顯著。根據(jù)回歸模型各因素影響程度的大小，得到各因素影響程度主次順序為x1、x2、x4和x3。

3.3 交互因素對清選性能的影響

根據(jù)上述二次回歸模型式（3）和式（4），進行雙因素響應(yīng)面分析，結(jié)果如圖4 所示[8-9]。

由圖4a 可知，當(dāng)振動頻率固定時，隨下篩篩面傾角的增大，籽粒含雜率先減小后增大，其主要原因是，隨著篩面傾角增大，籽粒在篩面上的前進速度加快、松散程度增加，雜余有效透篩面積減小，雜余透篩概率降低，故籽粒含雜率下降。同時，導(dǎo)致綠豆脫出物在篩面上停留的時間縮短，降低了籽粒透篩的概率，使清選損失率增加，當(dāng)清選損失率增大至一定程度后篩下雜余比例增加，籽粒含雜率隨之增高。當(dāng)下篩篩面傾角固定時，隨振動頻率的增大籽粒含雜率不斷減小，其主要原因是隨著振動頻率的增大，篩面上物料可以快速分層，雜余在篩面上停留時間縮短，雜余透篩概率降低，故籽粒含雜率減小。

由圖4b 可知，當(dāng)振動頻率固定時，隨風(fēng)速的增大，清選損失率不斷增大，其主要原因是隨著風(fēng)速的增大，部分輕質(zhì)籽粒隨雜余被吹出清選室，故清選損失率增大。當(dāng)風(fēng)速固定時，隨振動頻率的增大，清選損失率不斷增大，其主要原因是上下篩橫向振幅均大于縱向振動，篩體振動頻率越高，物料向前運動的速度越快，籽粒在篩面停留時間變短，籽粒透篩概率降低，故清選損失率增大。

3.4 參數(shù)優(yōu)化和試驗驗證

為獲得清選篩最佳清選性能參數(shù)，利用Design-Expert 軟件對試驗結(jié)果進行參數(shù)優(yōu)化，約束條件為

優(yōu)化求解得到最優(yōu)參數(shù)組合為入風(fēng)口風(fēng)速5.7 m/s、清選篩振動頻率4.71 Hz、上篩橫向振幅7.99 mm、下篩篩面傾角0.68°，此條件下籽粒含雜率2.12% 及清選損失率2.51%。

為驗證正交試驗結(jié)果的準(zhǔn)確性，對最優(yōu)參數(shù)組合進行虛擬試驗驗證，同組試驗重復(fù)3 次，結(jié)果取平均值。由表9 可知，籽粒含雜率、清選損失率的誤差率均不超過2%，表明預(yù)測模型可靠，試驗最優(yōu)解有效，可用于后期田間試驗[10]。

4 結(jié)束語

（1）針對目前綠豆清選過程籽粒含雜率和清選損失率高的問題，設(shè)計一種雙斜面風(fēng)篩式綠豆清選裝置。

（2）通過二次回歸方程模擬可知，各因素對籽粒含雜率影響程度主次順序為入風(fēng)口風(fēng)速、振動頻率、上篩橫向振幅和下篩篩面傾角；各因素對清選損失率影響程度主次順序為入風(fēng)口風(fēng)速、振動頻率、下篩篩面傾角和上篩橫向振幅。通過響應(yīng)曲面分析可知，振動頻率和下篩篩面傾角的交互效應(yīng)對籽粒含雜率影響顯著；入風(fēng)口風(fēng)速和振動頻率的交互效應(yīng)對清選損失率影響顯著。

（3）利用Design-Expert 軟件，對清選裝置參數(shù)進行優(yōu)化，得到最佳參數(shù)組合。該參數(shù)組合經(jīng)驗證符合性能指標(biāo)要求，可用于后期田間試驗。

參考文獻(xiàn)

［1］李洪昌．風(fēng)篩式清選裝置理論及試驗研究[D]．鎮(zhèn)江：江蘇大學(xué)，2011．

LI Hongchang． Theoretical and experimental study on air-and-screencleaning unit[D]．Zhenjiang：Jiangsu University，2011．

［2］侯俊銘，任兆坦，朱紅杰．雙層傾斜振動風(fēng)篩式蓖麻清選裝置設(shè)計與試驗[J]．農(nóng)業(yè)機械學(xué)報，2022，53（S2）：39-51．

HOU Junming， REN Zhaotan， ZHU Hongjie． Design and test ofdouble-layer inclined vibrating air-screen castor cleaning device[J]． Transactions of the Chinese Society for Agricultural Machinery，2022，53（S2）：39-51．

［3］袁子翔，戴飛，趙武云，等．基于CFD-DEM 的胡麻風(fēng)篩式清選裝置仿真分析與試驗[J/OL]． 干旱地區(qū)農(nóng)業(yè)研究： 1-10． http：//kns.cnki.net/kcms/detail/61.1088.S.20230927.0858.004.html．

YUAN Zixiang，DAI Fei，ZHAO Wuyun，et al．Simulation analysisand experiment of air-screen cleaning devicefor flax based on CFD-DEMcoupling[J/OL]．Agricultural Research in the Arid Areas：1-10．http：//kns.cnki.net/kcms/detail/61.1088.S.20230927.0858.004.html．

［4］張勝偉，張瑞雨，陳天佑，等．綠豆種子離散元仿真參數(shù)標(biāo)定與排種試驗[J]．農(nóng)業(yè)機械學(xué)報，2022，53（3）：71-79．

ZHANG Shengwei，ZHANG Ruiyu，CHEN Tianyou，et al．Calibrationof simulation parameters of mung bean seeds using discrete elementmethod and verification of seed-metering test[J]． Transactions of theChinese Society for Agricultural Machinery，2022，53（3）：71-79．

［5］徐福龍，賀俊林，王月華，等．藍(lán)亞麻蒴果離散元模型仿真參數(shù)標(biāo)定[J]．甘肅農(nóng)業(yè)大學(xué)學(xué)報，2021，56（4）：144-153．

XU Fulong， HE Junlin， WANG Yuehua， et al． Simulation parameterscalibration about capsule discrete element model of blueflax[J]．Journal of Gansu Agricultural University，2021，56（4）：144-153．

［6］李姣姣，賀俊林，徐福龍，等．基于EDEM-Fluent 耦合胡麻清選裝置過程模擬分析[J]．農(nóng)業(yè)工程，2021，11（12）：76-81．

LI Jiaojiao， HE Junlin， XU Fulong， et al． Process simulation and analysis of flax cleaning device based on EDEM-Fluent coupling[J]．Agricultural Engineering，2019，11（12）：76-81．

［7］占才學(xué)，陳巧敏，肖宏儒，等．食用豆聯(lián)合收獲脫粒清選裝置適用性分析[J]．中國農(nóng)機化學(xué)報，2020，41（6）：13-19．

ZHAN Caixue，CHEN Qiaomin，XIAO Hongru，et al．Applicabilityanalyses of combined harvest threshing and cleaning device for ediblebeans[J]． Journal of Chinese Agricultural Mechanization， 2020，41（6）：13-19．

［8］徐鵬慶，戴飛，趙武云，等．胡麻脫粒清選裝置作業(yè)過程仿真與試驗[J/OL]． 農(nóng)業(yè)機械學(xué)報： 1-13． http：//kns.cnki.net/kcms/detail/11.1964.S.20230911.1913.034.html．

XU Pengqing， DAI Fei， ZHAO Wuyun， et al． Simulation analysisand experiment of operation process of flax threshing and cleaningdevice[J/OL]．Transactions of the Chinese Society for Agricultural Ma-chinery：1-13．http：//kns.cnki.net/kcms/detail/11.1964.S.20230911.1913.034.html．

［9］吳楠，賀俊林，劉少華，等．基于EDEM-Fluent 耦合的鈣果風(fēng)篩式清選裝置仿真與試驗[J]．農(nóng)業(yè)工程，2021，11（1）：82-87．

WU Nan， HE Junlin， LIU Shaohua， et al． Simulation and experimentof air screen cleaning device for Cerasus humilis based on EDEMFluentcoupling[J]． Agricultural Engineering， 2021， 11（ 1） ： 82-87．

［10］王欣，徐智．全自動紅薯殺秧機的設(shè)計與試驗：基于Design-Expert8.0.6[J]．農(nóng)機化研究，2021，43（6）：89-94．

WANG Xin， XU Zhi． Design optimization and experiment of fullyautomatic sweet potato seedling killer： based on Design-Expert8.0.6[J]． Journal of Agricultural Mechanization Research， 2021，43（6）：89-94．

基金項目： 山西農(nóng)業(yè)大學(xué)科技創(chuàng)新基金項目（Zdpy201906）