兩自由度往復式滾筒篩茶葉篩分模擬研究

王成軍，章天雨，鄭艷，陳金燕，胡松

1. 安徽理工大學 礦業工程博士后流動站，安徽 淮南 232001；2. 安徽理工大學 機械工程學院，安徽 淮南 232001

兩自由度往復式滾筒篩茶葉篩分模擬研究

王成軍1,2，章天雨2，鄭艷2，陳金燕2，胡松2

1. 安徽理工大學 礦業工程博士后流動站，安徽 淮南 232001；2. 安徽理工大學 機械工程學院，安徽 淮南 232001

為研究兩自由度往復式滾筒篩的運動和結構參數對茶葉篩分效率的影響規律，以涌溪火青茶為對象，基于離散元法運用EDEM軟件在兩自由度往復式滾筒篩中進行茶葉顆粒篩分模擬。研究了滾筒在往復式回轉運動時滾筒傾角、往復角和往復頻率對篩分過程的影響規律，得出滾筒篩的理想參數取值范圍為：傾角為4～8°，往復頻率為2～2.5 Hz，往復角為60～75°。進一步對滾筒施加X向振動激勵，采用正交試驗的方法研究激振頻率、振幅、往復角、對篩分過程的影響。仿真結果表明，激振頻率和往復角對篩分效率影響顯著；兩自由度篩分模式下，滾筒篩的理想參數為：激振頻率為8 Hz，振幅為8 mm，往復角為70°。

涌溪火青茶；往復式滾筒篩；兩自由度；篩分；離散元法

目前，國內外學者根據茶鮮葉及葉梗的物理特性，依托風選、篩選、重力選、色選、靜電及光電原理，并融合機械、電子、圖像識別等技術，研制出茶葉分選設備，提高了工作效率，但茶葉分選的質量和精度不夠理想[1]。平面圓篩機和滾筒式分級機，因結構簡單、成本低，在茶葉精制過程中應用廣泛。Jeffrey等[2]給出了滾筒直徑、長度的計算公式；楊茂成[3]經研究得出當滾筒直徑與長度比為1∶4～1∶6時篩分效果較為理想，滾筒的最佳轉速為11/R1/2～15/R1/2r·min-1；唐萌[4]、張蘭蘭[5]的研究表明，將滾筒篩機與風選機組合篩選，有利于提高茶鮮葉分級效率；胡永光等[6]設計了一種參數可調的等直徑滾筒式篩分機，得出滾筒傾角、進料率和滾筒轉速的最優組合。謝俊等[7]、張曉慶[8]、王南等[9]將并聯機構引入茶葉篩分機構，表明多自由度運動有利于茶葉的分層、分散和透篩。此外，EDEM在農業物料的仿真模擬中得到了廣泛應用[10-11]。

現有茶葉滾筒篩的研究都是基于滾筒作單向整周回轉運動，物料的運動和篩網利用率受限。本文結合EDEM，以足干后的涌溪火青茶為研究對象，分析滾筒篩作往復式回轉運動和兩自由度運動時的篩分過程，得到工作參數的優化組合，以揭示在新的運動形式下滾筒篩內茶葉顆粒的運動及透篩規律。

1 滾筒篩試驗條件設置

1.1顆粒形狀

本文以足干后的涌溪火青為對象，對大量此類成品茶葉（圖1）進行觀察和測量，發現茶葉顆粒外形大多類似由多個球體組合的類橢球顆粒[12]。在EDEM中，設置茶葉顆粒為多球體組合元（圖2），且顆粒的平均長徑比L/d=1.5，改變球體直徑、球體數目、相鄰球體位置，即可產生不同形狀的茶葉顆粒模型。

1.2顆粒接觸模型

足干后的涌溪火青茶的篩分是茶葉機制工藝的最后一步[13]，篩前茶葉含水量僅為鮮葉的26%。因此，模擬實驗中可采用軟干球接觸模型進行茶葉顆粒間的碰撞模擬。在離散元模擬中，每個顆粒的平移和轉動通過牛頓和歐拉方程進行追蹤。

圖3為兩個多球顆粒i和j的碰撞簡圖[15]，其中球面l和k接觸。

基于線性彈簧-阻尼器模型的接觸力的法向分量可表示為：

接觸力的切向分量可用受庫侖條件約束的線性彈簧模型表示為：

1.3滾筒篩模型

考慮到計算機的內存和計算量，對滾筒篩進行簡化和等比例縮小（圖4）。初步選取滾筒篩模型直徑D=0.2 m，滾筒傾角α=5°，參照Jeffrey[2]的計算公式，進行滾筒尺寸設計，當滾筒轉速取臨界轉速的50%時，可得滾筒長度為：

圖1　成品茶葉顆粒Fig. 1 Tea particles of final product

圖2　茶葉顆粒模型Fig. 2 The tea particle model

圖3　兩個多球顆粒的碰撞簡圖Fig. 3 The collision diagram of two multi-sphere

圖4　滾筒篩仿真模型Fig. 4 simulation model of roller

2 數值模擬及分析

滾筒為D=0.2 m，L=0.51 m的圓柱體，篩孔尺寸7 mm，開孔率為75%。入料茶葉由圖2所示的5 250個顆粒組成，粒徑為4、5、6、7 mm的顆粒分別占總體顆粒數量的10%、60%、20%、10%。為簡化模擬顆粒設置，取每種粒徑的顆粒對應的長徑比及其所占此種顆粒總數的百分比分別為：長徑比1.2，30%；長徑比1.5，50%；長徑比1.7，20%。入料時間為15 s，仿真時間為20 s。茶葉顆粒和篩面的材料特性及接觸參數見表1。

2.1單自由度運動時篩分過程影響因素研究

式（6）中，η為t時刻可篩物料的篩凈率，At為0～t時刻篩下物料總質量（kg）；St為0～t時刻，可篩物料總入料質量（kg）。

雜余顆粒在篩面長時間滯留，大量堆積必然影響篩分效率。因此，研究雜余排出率對物料篩分過程的影響十分必要，雜余排出率可定義為：

式（7）中，Z為t時刻雜余排出率，Bt為0～t時刻雜余排出總質量（kg）；Mt為0～t時刻，雜余總入料質量（kg）。

相對誤排率可定義為：

在式（8）中，Wi為第i種工作狀態下的可篩物料相對誤排率，Ki為第i種工作狀態下排料口中可篩物料質量（kg）；Kmax為在一組試驗中對應的排料口中可篩物料的最大質量（kg）。

2.1.1滾筒傾角對篩分過程的影響

表1　茶葉顆粒和篩面的材料特性和接觸參數Table 1 Material properties and contact characteristic parameters of the tea particles and screening surface

圖5　滾筒傾角對篩分過程影響規律Fig. 5 The effect of roller inclination on screening process

2.1.2往復頻率對篩分過程的影響

2.1.3往復角對篩分過程的影響

圖6　往復頻率對篩分過程影響規律Fig. 6 The effect of reciprocating frequency on screening process

圖7　往復角對篩分過程影響規律Fig. 7 The effect of reciprocating angle on screening process

作為正弦旋轉運動的另一重要參數，往復角φ決定滾筒在1/4周期內的擺動角度。取往復頻率f1=2 Hz，傾角α=6°，分別對往復角φ=45～95°的滾筒篩進行模擬，仿真結果如圖7所示。當φ=45～60°時，篩凈率從46.56%上升到95.75%；當φ=60～75°時，篩凈率略有提高，最大為99.62%；當φ＞75°時，篩凈率隨往復角的增大而出現小幅下降。當φ=45～70°時，相對誤排率基本為0；當φ＞70°時，相對誤排率增幅明顯。當φ=45～55°時，雜余排出率基本為0；當φ＞60°時，雜余排出率以較大增幅增至φ=95°時的80.19%。其原因是，當往復角φ=45～60°時，茶葉運動幅度小，分散效果差，堆積現象明顯，茶葉不能及時透篩與排出；當往復角φ＞60°時，茶葉在篩網上滑動的面積增加，拋擲能力增強，料層厚度均勻，透篩率明顯增加，茶葉流向排料端的速度和數量增加，篩凈率和雜余排出率提高。但過大的往復角，會導致物料運動速度過大，被排出的可篩茶葉顆粒增多，相對誤排率增大，同時降低了篩凈率。

在本試驗條件下，經過以上對傾角、往復頻率、往復角的單因素試驗，綜合考慮篩分效率和生產率，滾筒篩往復式回轉運動的理想參數取值范圍為：傾角α=4～8°，往復頻率f1=2～2.5 Hz，往復角φ=60～75°。

2.2兩自由度運動對篩分過程的影響

2.2.1因子水平的確定

為研究兩自由度運動對篩分過程的影響，對滾筒施加X向振動激勵。由2.1.2分析得到，滾筒往復頻率f1=1.5 Hz時，篩分效果最差，在這種工作狀態下，對滾筒施加X向振動激勵，有利于研究X向激振參數對物料篩分過程的影響。因此取滾筒往復頻率f1=1.5 Hz，傾角α=6°，并對滾筒施加X向激振。采用正交試驗的方法，根據仿真模擬研究，選取X向激振頻率f2、振幅A和往復角φ，作為試驗因素，確定各因素水平，如表2所示。采用正交表L16(45)安排試驗方案，如表3所示。

2.2.2正交試驗結果及分析

以模擬進行到15 s時的滾筒篩篩凈率作為表3中試驗結果。根據極差分析得到，影響篩凈率的因子主次順序為激振頻率f2(因素B）、往復角φ（因素D）、振幅A（因素C）。對于單因素，按主次排列，因素B對篩凈率的影響順序為B3、B1、B2、B4；因素C對篩凈率影響順序為C4、C2、C3、C1；因素D對篩凈率影響順序為D1、D3、D2、D4。

表2　因子水平表Table 2 Table of factor levels

結合表4的方差分析結果可知：對于給定顯著性水平α=0.05，因素B（激振頻率）和因素D（往復角）對篩凈率有顯著影響。

綜合以上分析，當滾筒往復頻率f1=1.5 Hz，傾角α=6°時，兩自由度往復式滾筒篩的最佳工作參數組合為B3C4D1，即激振頻率f2=8 Hz，振幅A=8 mm，往復角φ=70°。

2.2.3優選方案的數值模擬與對比分析

結合正交試驗結果，在滾筒傾角α=6°，往復頻率f1=1.5 Hz,往復角φ=70°時，取X向激振頻率f2=8 Hz，振幅A=8 mm，在EDEM中，分別模擬滾筒篩作單自由度運動和兩自由度運動時的篩分過程，得到如圖8所示的兩種篩分模式的終篩狀態對比圖。其中，圖8-a為單自由度篩分模式下，即滾筒篩僅作往復式回轉運動時，滾筒篩的終篩狀態；圖8-b為兩自由度篩分模式下，即在施加X向激振后，滾筒篩的終篩狀態。

仿真結果表明，單自由度篩分模式下，茶葉顆粒在篩面大量堆積，t=20 s時，茶葉尚無排出，雜余排出率為0，篩凈率僅40.21%；施加X向激振后，雜余排出率為88%，篩凈率為99.26%。因此，X向激振有利于茶葉在往復式滾筒篩中的分散和透篩，這與文獻[7]和文獻[10]中的研究結果類似。

表3　正交試驗方案〔L16(45)〕及結果Table 3 Design [L16(45)] and results of orthogonal test

表4　方差分析表Table 4 Variance analysis

圖8　兩種篩分模式的終篩狀態對比圖Fig. 8 The contrast diagram of the final state under two screening modes

3 結論

（1）提出以往復式回轉運動替代傳統滾筒篩的單向整周轉動，研究了滾筒篩傾角、往復頻率和往復角，對篩分過程的影響規律，得到滾筒篩往復式回轉運動的理想參數取值范圍為：傾角α=4～8°，往復頻率f2=2～2.5 Hz，往復角φ=60～75°。

（2）對滾筒篩施加X向激振，正交試驗結果表明，在激振頻率、振幅和往復角，3個因素中，激振頻率和往復角對篩凈率影響顯著。得到兩自由度往復式滾筒篩的最佳工作參數組合為：f2=8 Hz，A=8 mm，φ=70°。通過對比兩自由度茶葉滾筒篩和單自由度茶葉滾筒篩的篩分結果，得到在兩自由度篩分模式下，篩凈率提高幅度達59.05%，雜余排出率提高幅度達88%，證明了X向激振有利于往復式滾筒篩中茶葉的分散和透篩。

（3）運用三維離散元法，結合EDEM數值模擬，為滾筒篩在兩自由度往復式運動模式下的篩分過程研究及篩分設備的研制提供參考依據。

[1] 張皓臻, 肖宏儒, 梅松, 等. 茶葉分選技術的現狀和未來[J]. 農業工程, 2014, 4(5): 59-63.

[2] Jeffrey W S, Rodfers M H, James F S. The place of the trommel in resource recovery [C]//The American Society of Mechanical Engineers. Proceedings of 1992 National Waste Processing Conference: Solid waste processing into the 21st century. New York: American Society of Mechanical Engineers, 1992: 1-8.

[3] 楊茂成. 滾筒圓篩機工藝設計計算[J]. 貴州茶葉, 1994, 1(77): 34-36.

[4] 唐萌. 茶園名優茶機械化采摘集成技術研究[D]. 杭州:浙江大學, 2007.

[5] 張蘭蘭. 機采名優茶鮮葉分級技術及分級研究[D]. 杭州:浙江大學, 2012.

[6] 胡永光, 李建鋼, 陸海燕, 等. 等直徑滾筒式茶鮮葉分級機設計與實驗[J]. 農業機械學報，2015, 46(增刊): 116-121.

[7] 謝俊, 張曉慶, 王成軍, 等. 三平移并聯茶葉篩分機構設計及運動仿真[J]. 工程設計學報學報, 2012, 19(3): 208-212.

[8] 張曉慶. 基于混聯機構的茶葉篩分機的設計與研究[D].鎮江: 江蘇大學, 2013.

[9] 王南, 周海棟, 崔國華, 等. 基于三轉動并聯機構的茶葉篩分裝置運動學性能研究[J]. 工程設計學報, 2014, 21(3): 266-272.

[10] 王成軍, 劉瓊, 馬履中, 等. 棉籽顆粒在三自由度混聯振動篩面上的運動規律[J]. 農業工程學報, 2015, 31(6): 49-55.

[11] 賈富國, 姚麗娜, 韓燕龍, 等. 基于離散元法的糙米勻料盤仿真優化設計[J]. 農業工程學報, 2016, 32(4): 235-241.

[12] 趙歷. 組合橢球顆粒模型的三維離散單元方法的研究[D].長春: 吉林大學, 2010.

[13] 董永泓. 涌溪火青茶機制工藝標準[J]. 茶葉機械雜志, 2001(2): 20-21.

[14] Munjiza A, Latham J P, John N W M. 3D dynamics of discrete element systems comprising irregular discrete elements-integration solution for finite rotations in 3D[J]. International Journal for Numerical Methods in Engineering, 2002, 56(17): 33-55.

[15] Frederik Elskamp, Harald Kruggel-Emden, Manuel Hennig, et al. Benchmarking of process models for continuous screening based on discrete element simulations [J]. Minerals Engineering, 2015(83): 78-96.

[16] Kruggel-Emden H, Wirtz S, Scherer V. A study on tangential force laws applicable to the discrete element method (DEM) for materials with viscoelastic or plastic behavior [J]. Chemical Engineering Science, 2008, 63(6): 1523-1541.

Simulation Study on the Screening of the Reciprocating Roller Screen with Two Degrees of Freedom for Tea

WANG Chengjun1,2, ZHANG Tianyu2, ZHENG Yan2, CHEN Jinyan2, HU Song2
1. Post-doctoral Mobile Station of Mining Engineering, Anhui University of Science and Technology, Huainan 232001, China; 2. School of Mechanical Engineering, Anhui University of Science and Technology, Huainan 232001, China

Taking Yongxi Huoqing tea as research object, the EDEM based on the discrete element method was used for simulation study the influence of the kinematic and structure parameters of reciprocating roller screen with two degrees of freedom on screening efficiency for tea leaves. The effects of roller inclination, reciprocating angle and reciprocating frequency on screening was analyzed when the roller was on reciprocating rotational movement. It was found that the ideal parameters of the roller were: 4-6° (inclination), 2-2.5 Hz (reciprocating frequency) and 60-75° (reciprocating angle). Furthermore, an X-direction excitation was added to the roller screen to study the effects of excitation frequency, amplitude and reciprocating angle on the screening process by orthogonal test. The simulation results showed that the ideal parameters of the roller screen were: 8 Hz (excitation frequency), 8 mm (amplitude) and 70° (reciprocating angle) under screening mode with two degrees of Freedom.

Yongxi Huoqing tea, reciprocating roller screen, two degrees of freedom, screen, discrete element method

TS272.3

A

1000-369X（2016）06-613-08

2016-07-11

2016-08-12

中國博士后科學基金資助項目（2013M540507）

王成軍，男，副教授，碩士生導師，工學博士，主要從事并聯機構及多維振動篩分技術研究。E-mail: cumt1279@163.com