章啟祥 劉發(fā)英 楊振宇 李繼偉 韓夢杰 李學強

摘要：在馬鈴薯輸送與分選過程中，針對機械損傷引起表皮破損導致儲存時極易腐爛變質的問題。通過對馬鈴薯的碰撞分析確定馬鈴薯損傷的影響因素，獲取分選裝置運行參數(shù)的取值范圍，并以傷薯率為評估指標，以輸送帶速度、仿形滾輪移動速度、馬鈴薯滑落落差為試驗因素進行離散元仿真分析的二次正交旋轉回歸試驗，分析各因素對評估指標的影響規(guī)律。參數(shù)優(yōu)化結果表明，輸送帶速度為0.32 m/s，馬鈴薯滑落落差為103.5 mm，仿形滾輪移動速度為0.49 m/s時，最低傷薯率為1.98%。經試驗驗證，滿足馬鈴薯收獲分選時的國家標準和行業(yè)標準需求，可為同類型的馬鈴薯機械的設計和優(yōu)化提供技術參考。

關鍵詞：馬鈴薯；分選裝置；機械損傷；離散元；參數(shù)優(yōu)化

中圖分類號：S233.73

文獻標識碼：A

文章編號：2095-5553 （2023） 04-0096-08

Abstract： In the process of potato transportation and sorting， aiming at the problem that the skin damage caused by mechanical damage is very easy to rot and deteriorate during storage. Through the collision analysis of potatoes， the influencing factors of potato damage are determined， and the value range of the operating parameters of the sorting device is obtained. Taking the potato damage rate as the evaluation index， the conveyor belt speed， the movement speed of profiling roller and the potato slip drop as the test factors， the quadratic orthogonal rotation regression test of discrete element simulation analysis is carried out to analyze the influence law of various factors on the evaluation index. The results of parameter optimization showed that when the speed of conveyor belt was 0.32 m/s， the falling drop of potato was 103.5 mm， and the moving speed of profiling roller was 0.49 m/s， the minimum potato injury rate was 1.98%. The experimental results show that it meets the requirements of national standards and industrial standards during potato harvesting and sorting， and can provide technical reference for the design and optimization of the same type of potato machinery.

Keywords：? potato; separator; mechanical damage; discrete element; parameter optimization

0 引言

2015年，原農業(yè)部提出馬鈴薯主糧化策略，馬鈴薯成為第四大主糧，經濟地位日益重要，馬鈴薯機械化生產對馬鈴薯的增產增收、提高經濟效益有著重大作用［1-3］。國外對馬鈴薯損傷進行了一些研究［4-7］，國外學者的主要研究內容是：采用掃描電鏡（SEM）分析、沖擊鐘擺實驗等新技術，通過試驗研究馬鈴薯損傷機理并確定各個馬鈴薯機械的參數(shù)；國內學者近年來也做了部分研究，桑永英等［8］對馬鈴薯碰撞傷害的原因進行了分析，并建立了馬鈴薯碰撞模型；郭世魯?shù)龋?］使用運動學對馬鈴薯從馬鈴薯收獲機分離的過程進行分析，得到轉速﹑臨界機械碰撞的關系；洪翔等［10］介紹了一種測量馬鈴薯受到碰撞損傷程度的方法；兌瀚［11］對馬鈴薯收獲機升運過程中的機械損傷進行了分析與實驗；魏忠彩等［12］研究了馬鈴薯機械化輸送過程中的傷薯問題；馮斌［13］對馬鈴薯碰撞損傷機理進行了研究。目前，對馬鈴薯損傷的研究主要是馬鈴薯收獲過程中機械損傷的損傷機理研究［14-16］。在馬鈴薯分選過程中造成的機械損傷研究較少，然而馬鈴薯分選是其在入庫或出庫過程中的一個重要環(huán)節(jié)，由于機械損傷在此過程中造成的馬鈴薯表面破損將導致其在儲存時極易腐爛變質，因此需要確定馬鈴薯分選過程分選裝置的最佳運行參數(shù)，降低馬鈴薯表面破損率。

由于馬鈴薯在分選入庫過程中，產生機械損傷主要與馬鈴薯分選裝置的運行參數(shù)（馬鈴薯滑落落差、輸送帶速度和仿形滾輪移動速度）有關。通過分析馬鈴薯分選過程中的碰撞損傷過程，確定影響馬鈴薯碰撞損傷的主要因素，根據(jù)馬鈴薯分選裝置的運行參數(shù)的取值范圍，利用EDEM軟件進行離散元仿真分析，獲取最佳的分選裝置運行參數(shù)，最后在馬鈴薯分選裝置試驗臺上驗證離散元仿真分析的結果。

1 馬鈴薯分選過程的碰撞分析

1.1 馬鈴薯分選系統(tǒng)組成與工作原理

馬鈴薯分選系統(tǒng)由輸送裝置、馬鈴薯分選裝置和收集框組成，如圖1所示。其中，輸送裝置由支架、輸送帶和電機等部分組成；馬鈴薯分選裝置由底座支架、輸送翻轉裝置、圖像采集裝置、剔除裝置和電動滾筒等部分組成；輸送翻轉裝置由電動滾筒驅動鏈傳動，鏈傳動在帶動仿形滾輪移動的同時使得仿形滾輪滾動，實現(xiàn)其支撐的馬鈴薯在輸送過程中翻滾，圖像采集檢測裝置包括工業(yè)相機、控制器、工控機、觸摸屏和傳感器等部分組成，可實現(xiàn)對馬鈴薯表面全面域圖像的獲取、顯示和分析處理，并將處理信號通過控制器傳輸給剔除裝置；剔除裝置由驅動氣缸和剔除杠桿組成，通過氣缸驅動剔除杠桿實現(xiàn)不合格馬鈴薯的剔除；兩個收集框位于馬鈴薯分選裝置的右面和后面，可分別收集不合格薯和合格薯。

馬鈴薯出入庫過程中，馬鈴薯通過輸送裝置不斷地輸送到馬鈴薯分選裝置的仿形滾輪上，在仿形滾輪自轉和平移的過程中，利用仿形滾輪與馬鈴薯表面產生的摩擦力使馬鈴薯翻轉。通過工業(yè)相機獲取馬鈴薯全面域特征圖像，采用機器視覺技術和缺陷檢測算法對馬鈴薯外部缺陷進行在線識別和檢測分選，發(fā)出控制指令。剔除裝置接收控制系統(tǒng)的指令，可實時對被檢測的不合格馬鈴薯進行剔除操作。在整個馬鈴薯分選系統(tǒng)中，馬鈴薯分選裝置是核心裝置，其主要技術參數(shù)如表1所示。

馬鈴薯損傷等級按以下方式進行判斷：（1）特級馬鈴薯：大小均勻，外觀新鮮，硬實，清潔、無泥土、無雜物，成熟度好，薯形好，基本無表皮破損、無機械損傷，無內部缺陷及外部缺陷造成的損傷，單薯的質量不低于150 g。（2）一級馬鈴薯：大小較均勻，外觀新鮮，硬實，清潔、無泥土、無雜物，成熟度較好，薯形較好，輕度表皮破損及機械損傷，內部缺陷及外部缺陷造成的輕度損傷，單薯的質量不低于100 g。（3）二級馬鈴薯：大小較均勻，外觀較新鮮，較清潔、允許有少量泥土和雜物，中度表皮破損，無嚴重畸形，無內部缺陷及外部缺陷造成的嚴重損傷，單薯的質量不低于50 g。

1.2 馬鈴薯碰撞分析

馬鈴薯通過輸送帶滑落到分選裝置的仿形滾輪上的過程如圖2所示，馬鈴薯由輸送帶提升輸送，帶有一定的平拋初速度滑落到仿形滾輪上。在馬鈴薯從輸送帶滑落到仿形滾輪上時，馬鈴薯與仿形滾輪的碰撞過程中，馬鈴薯的機械損傷程度大，馬鈴薯碰撞過程中接觸區(qū)的應力分布可用準靜態(tài)Hertz理論進行分析。

按照Hertz理論，根據(jù)馬鈴薯分選的環(huán)境狀況，作出如下假設。

1） 馬鈴薯的表面是不間斷的，而且是非協(xié)調的。

2） 馬鈴薯與仿形滾輪的接觸區(qū)的有效尺寸遠小于相對曲率半徑。

3） 馬鈴薯與仿形滾輪都可看作半彈性空間，接觸區(qū)的有效尺寸遠小于馬鈴薯與仿形滾輪的有效半徑。

4） 馬鈴薯與仿形滾輪接觸時，無面內摩擦，切向力為0。

馬鈴薯與仿形滾輪碰撞過程的分析原理如圖3所示，馬鈴薯與仿形滾輪碰撞的形變狀態(tài)如圖4所示。

從式（5）中可以看出，馬鈴薯與仿形滾輪間的受力與馬鈴薯的質量、輸送帶速度、仿形滾輪移動速度、馬鈴薯滑落落差有關。

2 仿真試驗

2.1 基于離散元的模型建立

2.1.1 馬鈴薯離散元模型的建立

結合商品馬鈴薯外觀形狀特征，設置離散元仿真分析的馬鈴薯顆粒模型，如圖5所示。為了盡量還原馬鈴薯真實的形狀并且提高仿真效率及精度，馬鈴薯顆粒模型采用多球聚合模型。為了簡化模型，考慮到不需要對馬鈴薯表面的整個表面進行受力分析，而只需對馬鈴薯整體受力進行仿真分析即可，因此采用4個圓形顆粒構建馬鈴薯顆粒模型，進行受力仿真分析，設置其質量范圍為160～320 g。

2.1.2 輸送帶與分選裝置離散元模型的建立

為了能夠合理有效的仿真模擬與計算，結合馬鈴薯分選裝置的參數(shù)，對馬鈴薯分選裝置進行簡化，去掉機架，保留輸送帶和輸送仿形滾輪。

由于馬鈴薯是從輸送裝置后段的水平輸送帶滑落到仿形滾輪上，因此輸送帶模型只使用水平輸送帶的這一部分，進行離散元模擬仿真分析，簡化后的輸送帶與分選裝置的離散元模擬仿真分析模型如圖6所示。

2.1.3 仿真模型參數(shù)確定

由于馬鈴薯與輸送裝置、馬鈴薯與分選裝置之間不存在黏附力，在EDEM軟件內進行接觸模型設置，選取Hertz-Mindlin（no slip）［17］模型作為馬鈴薯與馬鈴薯、馬鈴薯與輸送裝置、馬鈴薯與分選裝置之間的接觸模型。

2.1.4 仿真設置

考慮到馬鈴薯的機械損傷與碰撞材料、下落高度、塊莖質量、含水率、跌落方向和馬鈴薯品種等因素有關［18-20］。

結合本文研究內容的實際情況，按照山東滕州地區(qū)種植的馬鈴薯“希森3號”的特征參數(shù)進行仿真參數(shù)設置；輸送帶采用橡膠材料（TPU）進行仿真參數(shù)設置；仿形滾輪采用塑料材料（PVC）進行仿真參數(shù)設置；輸送機兩側護板材料采用鋼板（Q235A）進行仿真參數(shù)設置。馬鈴薯與輸送裝置、分選裝置的離散元仿真參數(shù)如表2所示。

為了獲得馬鈴薯分選裝置的輸送帶速度v0、仿形滾輪移動速度v1和馬鈴薯滑落落差h的最佳運行參數(shù)，根據(jù)Box-Behnken試驗設計原理，設置試驗因素和水平如表3所示。

2.2.2 EDEM參數(shù)設置

馬鈴薯顆粒模型生成在平傳送帶上，顆粒的生成速度由輸送帶速度決定，重力方向設置為Z軸向下。為保證馬鈴薯在輸送帶上流暢運動，按照Reyleigh時間步長的16%來設置固定時間步長，其值為0.37 ms左右，網格單元尺寸為最小顆粒半徑的2.5倍；為了使仿真結果更精確，設置仿真次數(shù)為30次。

3 仿真結果分析與參數(shù)優(yōu)化

仿真試驗結果如表4所示，表中v0、h、v1表示各因素編碼值。

3.2 響應面分析

依據(jù)建立的馬鈴薯傷薯率回歸模型，依次將其中一個仿真分析的因素置于零水平，考慮其他兩個因素對仿真分析的指標影響，繪制其響應面圖。如圖7（a）所示，當仿形滾輪移動速度處于中心水平，輸送帶速度試驗因素水平一定時，隨著馬鈴薯滑落落差增加，傷薯率隨之增高；馬鈴薯滑落落差試驗因素水平一定時，輸送帶速度增大，傷薯率隨之增高。

如圖7（b）所示，當馬鈴薯滑落落差處于中心水平，輸送帶速度試驗因素水平一定時，隨著仿形滾輪移動速度的增快，傷薯率隨之增高；仿形滾輪移動速度試驗因素水平一定時，輸送帶速度增大，傷薯率隨之增高；如圖7（c）所示，當輸送帶速度處于中心水平，仿形滾輪移動速度試驗因素水平一定時，隨著馬鈴薯滑落落差的增大，傷薯率隨之增高；馬鈴薯滑落落差試驗因素水平一定時，仿形滾輪移動速度加快，傷薯率隨之增高。

通過優(yōu)化參數(shù)，得到100組優(yōu)化結果，將這100組優(yōu)化結果從小到大排列，得到輸送帶速度為0.32 m/s，馬鈴薯滑落落差為103.5 mm，仿形滾輪移動速度為0.49 m/s時，傷薯率最低為1.98%。從優(yōu)化后的結果來看，馬鈴薯滑落落差接近約束條件的偏小值，輸送帶速度和仿形滾輪移動速度均是其約束條件的中間值。實際生產過程中，馬鈴薯滑落落差越小，其機械損傷的概率越低，但受輸送帶驅動滾輪尺寸的限制，馬鈴薯滑落落差不可能無限制小。另外，仿形滾輪移動速度高于輸送帶速度有利于馬鈴薯在仿形滾輪上均隔分布。因此，以上優(yōu)化后的輸送帶速度、馬鈴薯滑落落差和仿形滾輪移動速度的參數(shù)可用。

3.4 試驗驗證

使用優(yōu)化后的輸送帶速度、馬鈴薯滑落落差和仿形滾輪移動速度的運行參數(shù)在馬鈴薯分選裝置進行試驗驗證。輸送帶速度和仿形滾輪移動速度使用VC6236P（Victor）測速儀進行檢測和調定。

試驗采用山東滕州地區(qū)種植的馬鈴薯“希森3號”作為分選對象，其平均含水率70%左右，平均質量230 g 左右。參照原農業(yè)部馬鈴薯收獲機械行業(yè)標準（NY/T 1130—2006）和馬鈴薯商品薯分級與檢驗規(guī)程國家標準（GB/T 31784—2015），進行傷薯率試驗。從100組優(yōu)化結果中選擇具有代表性的5組優(yōu)化結果，這5組優(yōu)化結果分別是：第1、25、50、75、100組的優(yōu)化結果，各組的輸送帶速度、馬鈴薯滑落落差和仿形滾輪移動速度的優(yōu)化參數(shù)和傷薯率優(yōu)化結果，如表6所示。按照優(yōu)化參數(shù)調定馬鈴薯分選裝置的運行參數(shù)，每組分別對60個馬鈴薯進行試驗，并依次進行機械損傷標記和統(tǒng)計。傷薯率的試驗結果，如表6所示。從傷薯率的優(yōu)化結果和試驗結果來看，試驗結果與優(yōu)化結果基本保持一致，沒有產生較大差異。隨著優(yōu)化結果的傷薯率不斷增大，試驗結果傷薯率的變化幅度更大，但在允差范圍之內。馬鈴薯分選裝置運行參數(shù)優(yōu)化后的最佳傷薯率2.11%的試驗結果符合國家標準和行業(yè)標準對商品薯進行分級與檢驗的要求。

4 結論

1） 通過對馬鈴薯分選過程中，其與輸送裝置和分選裝置之間的受力分析，獲得馬鈴薯機械損傷的主要影響因素是輸送帶速度、馬鈴薯滑落落差和仿形滾輪移動速度。

2） 建立基于離散元的分選裝置馬鈴薯運動過程模型，以輸送帶速度、馬鈴薯滑落落差、仿形滾輪移動速度為試驗指標，利用EDEM軟件對馬鈴薯在輸送帶與仿形滾輪上的運動過程進行離散元模擬仿真。使用Box-Behnken試驗設計原理進行試驗，確定了影響傷薯率因素的順序的從高到低排列順序依次為馬鈴薯滑落落差、輸送帶速度、仿形滾輪移動速度。

3） 利用Design Expert進行輸送帶速度、馬鈴薯滑落落差、仿形滾輪移動速度的參數(shù)優(yōu)化，得出輸送帶速度為0.32 m/s，馬鈴薯滑落落差為103.5 mm，仿形滾輪移動速度為0.49 m/s時，最低傷薯率為1.98%；采用三個優(yōu)化后的運行參數(shù)在樣機上進行試驗，傷薯率最低為2.11%，可滿足馬鈴薯收獲時分選作業(yè)的生產要求。

參 考 文 獻

［1］ 李娜， 周進， 崔中凱， 張華， 等. 山東省馬鈴薯生產全程機械化現(xiàn)狀與對策建議［J］. 中國農機化學報， 2019， 40（1）： 198-204.

Li Na， Zhou Jin， Cui Zhongkai， Zhang Hua， et al. Present situation and countermeasures of potato production mechanization in Shandong Province ［J］. Journal of Chinese Agricultural Mechanization， 2019， 40（1）： 198-204.

［2］ 竇青青， 孫永佳， 孫宜田， 等. 國內外馬鈴薯收獲機械現(xiàn)狀與發(fā)展［J］. 中國農機化學報， 2019， 40（09）： 206-210.

Dou Qingqing， Sun Yongjia， Sun Yitian， et al. Current situation and development of potato harvesting machinery at home and abroad ［J］. Journal of Chinese Agricultural Mechanization， 2019， 40（9）： 206-210.

［3］ 柳國光， 王濤， 姚愛萍， 等. 浙江省馬鈴薯機械化收獲關鍵技術探討［J］. 中國農機化學報， 2019， 40（8）： 48-52.

Liu Guoguang， Wang Tao， Yao Aiping， et al. Discussion on key technologies of potato mechanized technologies of potato mechanized harvesting in Zhejiang Province［J］. Journal of Chinese Agricultural Mechanization， 2019， 40（8）： 48-52.

［4］ Nikara S， Ahmadi E， Nia A A. Finite element simulation of the micromechanical changes of the tissue and cells of potato response to impact test during storage by scanning electron microscopy ［J］. Postharvest Biology and Technology， 2020， 164（1）： 111153.

［5］ Rady A M， Soliman S N. Evaluation of mechanical damage of Lady Rosetta potato tubers using different methods ［J］. International Journal of Postharvest Technology and Innovation， 2015， 5（2）： 125-148.

［6］ Azam M M，? Eissa A H A. Comprehensive evaluation of dynamic impact as a measure of potato quality ［J］. International Journal of Food Engineering and Technology， 2015， 1（1）： 1-10.

［7］ Finney E E， Hall C W， Thompson N R. Influence of variety and time upon the resistance of potatoes to mechanical damage ［J］. American Potato Journal， 1964， 41（6）： 178-186.

［8］ 桑永英， 張東興， 張梅梅. 馬鈴薯碰撞損傷試驗研究及有限元分析［J］. 中國農業(yè)大學學報， 2008， 13（1）： 81-84.

Sang Yongying， Zhang Dongxing， Zhang Meimei. Study on bruising damage experiment of potato and finite element analysis［J］. Journal of China Agricultural University， 2008， 13（1）： 81-84.

［9］ 郭世魯， 王衛(wèi)兵， 陳紹杰， 等. 小型馬鈴薯收獲機薯機碰撞有限元分析［J］. 中國農機化學報， 2016， 37（5）： 14-17.

Guo Shilu， Wang Weibing， Chen Shaojie， et al. Finite element analysis for crash safety of the potato and mechanism of small potato harvester working ［J］. Journal of Chinese Agricultural Mechanization， 2016， 37（5）： 14-17.

［10］ 洪翔， 萬陳， 王軍. 一種馬鈴薯臨界損傷跌落高度的測定方法［J］. 包裝學報， 2012， 4（3）： 30-33.

Hong Xiang， Wan Chen， Wang Jun. A method to determine the critical damage dropping height of potato ［J］. Packaging Journal， 2012， 4（3）： 30-33.

［11］ 兌瀚. 馬鈴薯收獲機升運過程機械損傷分析與試驗［D］. 哈爾濱： 東北農業(yè)大學， 2019.

Dui Han. Mechanical damage analysis and experimental research on potato harvester elevator ［D］. Harbin：Northeast Agricultural University， 2019.

［12］ 魏忠彩， 李學強， 孫傳祝， 等. 馬鈴薯收獲與清選分級機械化傷薯因素分析［J］. 中國農業(yè)科技導報， 2017， 19（8）： 63-70.

Wei Zhongcai， Li Xueqiang， Sun Chuanzhu， et al. Analysis of potato mechanical damage in harvesting and cleaning and sorting storage ［J］. Journal of Agricultural Science and Technology， 2017， 19（8）： 63-70.

［13］ 馮斌. 收獲期馬鈴薯塊莖物理特性及損傷機理研究［D］. 蘭州： 甘肅農業(yè)大學， 2018.

Feng Bin. Study on physical characteristics and damage of potato tubers at harvesting stage ［D］. Lanzhou： Gansu Agricultural University， 2018.

［14］ 盧琦. 馬鈴薯損傷機理試驗研究及聯(lián)合收獲機設計［D］. 楊凌： 西北農林科技大學， 2016.

Lu Qi. Experimental research on damage mechanism of potato and combine harvester designed ［D］. Yangling： Northwest A ＆ F University， 2016.

［15］ 胡奔. 馬鈴薯跌落損傷機理與防損傷裝置研究［D］. 成都： 西華大學， 2018.

Hu Ben. Study on potato dropping damage mechanism and damage-proof device ［D］. Chengdu： Xihua University， 2018.

［16］ 呂金慶， 楊曉涵， 呂伊寧， 等. 馬鈴薯挖掘機升運分離過程塊莖損傷機理分析與試驗［J］. 農業(yè)機械學報， 2020， 51（1）： 103-113.

Lü Jinqing， Yang Xiaohan， Lü Yining， et al. Analysis and experiment of potato damage in process of lifting and separating potato excavator［J］. Transactions of the Chinese Society for Agricultural Machinery， 2020， 51（1）： 103-113.

［17］ 蒙建國， 王春光， 謝勝仕， 等. 馬鈴薯恢復系數(shù)測定試驗分析［J］. 中國農業(yè)大學學報， 2017， 22（9）： 93-100.

Meng Jianguo， Wang Chunguang， Xie Shengshi， et al. Measurement test analysis of potato restitution coefficient ［J］. Journal of China Agricultural University， 2017， 22（9）： 93-100.

［18］ 馮斌， 孫偉， 石林榕， 等. 收獲期馬鈴薯塊莖碰撞恢復系數(shù)測定與影響因素分析［J］. 農業(yè)工程學報， 2017， 33（13）： 50-57.

Feng Bin， Sun Wei， Shi Linrong， et al. Determination of restitution coefficient of potato tubers collision in harvest and analysis of its influence factors ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2017， 33（13）： 50-57.

［19］ 石林榕， 孫偉， 趙武云， 等. 馬鈴薯種薯機械排種離散元仿真模型參數(shù)確定及驗證［J］. 農業(yè)工程學報， 2018， 34（6）： 35-42.

Shi Linrong， Sun Wei， Zhao Wuyun， et al. Parameter determination and validation of discrete element model of seed potato mechanical seeding ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2018， 34（6）： 35-42.

［20］ 王萬虎. 基于離散元的馬鈴薯撿拾裝置試驗研究［D］. 楊凌： 西北農林科技大學， 2017.

Wang Wanhu. Potato collecting device experimental study based on discrete element ［D］. Yangling： Northwest A & F University， 2017.