基于離散元法的臥式雙軸攪拌機攪拌優(yōu)化分析

肖舫

摘 要：提高顆粒物料混合的效率作為研究重點長期被人們關(guān)注。為了研究臥式雙軸攪拌機結(jié)構(gòu)參數(shù)和工藝參數(shù)對混合性能的影響，以兩種顆粒物料混合作為實驗對象，針對臥式雙軸攪拌機的結(jié)構(gòu)參數(shù)和攪拌工藝參數(shù)設(shè)計了多因素多水平的正交模擬實驗。以EDEM軟件為平臺，基于離散元法對顆粒模型進行了攪拌混合數(shù)值模擬實驗，繼而進行了極差和方差分析，得出結(jié)論：軸間距1000mm，葉片間距610mm，雙軸轉(zhuǎn)速為80rpm的參數(shù)水平下混合度參數(shù)結(jié)果最佳;攪拌結(jié)構(gòu)參數(shù)和工藝參數(shù)中，軸間距因素對混合度影響最大，雙軸轉(zhuǎn)速次之，攪拌軸葉片間距影響最小。表面合適的軸間距和恰當?shù)臄嚢杷俣葘嚢铏C混合性能影響顯著。混合攪拌機的參數(shù)優(yōu)化對眾多行業(yè)都有著極為重要意義。

0 引言

粉料工廠的攪拌機生產(chǎn)效率決定了工廠的產(chǎn)量和企業(yè)規(guī)模，而粉料產(chǎn)品的質(zhì)量往往也取決于顆?；旌显O(shè)備的混合性能，攪拌機是粉料業(yè)最為關(guān)鍵的加工設(shè)備[1]。國內(nèi)的攪拌機的研制和發(fā)展都朝著快速、高混合度、低耗、低殘留量以及系列化等目標[2]。其中，臥式雙軸槳葉攪拌機是目前最常見的機型，它采用雙槳葉攪拌機構(gòu)，可以縮短物料的混合時間，能夠?qū)⑽锪嫌行У幕旌希旌系木鶆蚨缺葐螛~的好。設(shè)計新型物料攪拌機、優(yōu)化攪拌機的結(jié)構(gòu)參數(shù)和工藝參數(shù)、提高攪拌機的工作效率一直是臥式雙軸槳葉攪拌機的研究重點[3]。

近幾年，采用數(shù)值仿真的方法開展臥式攪拌機關(guān)鍵技術(shù)的研究也越來越廣泛。涂文靜[4]通過使用軟件SolidWorks簡化雙軸攪拌機模型，并校核傳動軸強度及剛度，進而對攪拌機結(jié)構(gòu)進行優(yōu)化。詹民民[5]通過使用軟件Fluent分析不同的攪拌葉片安裝間距和不同的攪拌臂長度攪拌槳對臥式粉料攪拌機攪拌槽內(nèi)的流場分布以及固體相的體積分數(shù)分布情況的影響。丁智勇[6]通過使用軟件EDEM對瀝青混合料攪拌過程進行仿真，研究混合均勻性在不同轉(zhuǎn)速條件下的變化情況，進而確定最佳的攪拌速度。攪拌機的混合效率是粉料加工企業(yè)長期以來的研究內(nèi)容，提高攪拌混合機的效率也是為提高生產(chǎn)效率，實現(xiàn)收益增加的技術(shù)重點。因而，研究提高臥式雙軸攪拌機的工作性能有著重要的實際意義。

1攪拌機三維模型及顆粒模型

1.1攪拌機模型

研究開展所用的攪拌機模型通過UG繪制部件裝配，設(shè)計了包括箱體，葉片，檢修門，側(cè)擋板，攪拌軸，葉片體在內(nèi)的主要部件等。為了提高解算效率，在進行EDEM分析時，將攪拌軸、攪拌槳葉及攪拌外筒的三維模型在保證仿真真實性的情況下進行了部件的簡化，通過中轉(zhuǎn)文件完成后導(dǎo)入Geometry模塊。對EDEM中的雙軸臥式攪拌機模型進行了簡化改動，在不影響攪拌性能的情況下，攪拌保留了主箱體和雙攪拌軸主要部件，如圖1所示。

1.2顆粒模型

建立兩種材料顆粒模型如圖2所示，顏色分別定為白色和黑色，顆粒直徑分別為15mm和20mm。為了便于分析物料混合的均勻性受到各種條件的影響情況，需要采用數(shù)值化的方法評定攪拌腔內(nèi)散體物料混合程度，這里采用混合度作為評定物料分散均勻性和攪拌性能的指標參數(shù)[7]。

2 實驗設(shè)計

2.1正交實驗參數(shù)

在攪拌顆?；旌蠈嶒炛?，采用了多因素多水平的正交實驗具體分析了攪拌軸中心距、單軸上的葉片間距、雙軸轉(zhuǎn)速結(jié)構(gòu)參數(shù)對工藝性能的影響[8-9]?；陔x散元法對預(yù)先準備的物料進行混合攪拌正交仿真分析，相比于單因素水平實驗，正交試驗不僅能夠更準確地分析混合均勻性受各條件的影響情況，同時還能夠提高實驗效率。表1為實驗采用的因素水平。

本實驗因素為A：攪拌軸間距、B：軸上葉片的間距、C：攪拌軸轉(zhuǎn)速，是L9（33）三因素三水平正交實驗，建立的正交實驗方案如表2所示。

2.2仿真實驗設(shè)計

攪拌開始前，首先將建立的兩種物料顆粒模型按照一定的比例投放， 分別落人轉(zhuǎn)筒內(nèi)的混合腔?；旌蠑嚢枰欢ǖ臅r間后， 停機取樣。在EDEM中的攪拌實驗過程如圖3所示，為了便于觀察顆粒分布情況，將攪拌機箱體隱藏。

圖3? 顆?；旌螮DEM仿真

Figure 3 EDEM simulation of particle mixing

通過觀察混合腔的截面圖可以初步判斷物料的混合情況。以第一組因素水平實驗為例，經(jīng)過不同攪拌時間后混合狀態(tài)如圖4所示，可以看出，經(jīng)過10s的攪拌混合后，兩種顆粒的混合效果有了明顯的提高。為了得到攪拌完成后兩種物料顆粒的混合度準確值，需要將攪拌實驗的區(qū)域劃分為多個尺寸合理等份的子區(qū)域，借助EDEM的result獲取各組間兩種物料的組分比，以顆粒數(shù)50為篩選指標進行挑選，對有統(tǒng)計意義的有效組bin數(shù)據(jù)結(jié)果進行處理，進而獲取相關(guān)參數(shù)并計算得出混合度。

3? 實驗結(jié)果分析

圖5為各組實驗攪拌50s后顆粒混合情況，可以看出第1組、第3組及第9組中依然存在大面積未與白色顆粒相混合的黑色顆粒堆積，混合效果明顯差于其它幾組。經(jīng)過計算，第1組至第9組的最終混合度計算結(jié)果分別為：2.982、3.980、3.402、5.300、4.988、6.480、4.419、5.832、2.904。

3.1極差分析

差分析法簡稱R法，在正交實驗的數(shù)據(jù)處理中，極差分析法可以獲取影響結(jié)果因素的重要性將因素進行排列。表3為對混合度的極差進行分析結(jié)果。

從極差分析結(jié)果可以看出：對于因素A軸間距 kj2 > kj3 > kj1，表明軸間距因素在1000mm水平下，50s內(nèi)的混合度更高。對于因素B雙軸轉(zhuǎn)速 kj2 > kj3 > kj1，說明在50s內(nèi)葉片間距參數(shù)水平在610mm下的混合效果較優(yōu)。對于因素C雙軸轉(zhuǎn)速 kj2 > kj3 > kj1，說明在50s測試時間內(nèi)，雙軸轉(zhuǎn)速參數(shù)在80rpm水平下的混合效果較優(yōu)。

對ABC三因素的極差分析有： RA > RC > RB，說明在三種因素實驗所取水平梯度中，軸間距的顯著性最大，轉(zhuǎn)速的顯著性次之，葉片間距對目標參數(shù)的影響最小。

綜上所述：實驗因素在軸間距為1000mm，葉片間距610mm，雙軸轉(zhuǎn)速為80rpm的參數(shù)水平下混合度參數(shù)結(jié)果最佳，相應(yīng)的實驗組因素組合為A2B2C2實驗組。其次，軸間距對混合度影響是最大的，其次是雙軸轉(zhuǎn)速，對試驗時間內(nèi)的物料混合程度的影響最小因素是攪拌片間距。

3.2方差分析

通過采用方差分析方法對試驗結(jié)果進行分析可以對影響試驗結(jié)果的各因素的重要程度給以準確的定量估計，并對極差分析法的實驗數(shù)據(jù)波動進行分析。表4為對混合度的方差分析結(jié)果。

根據(jù)表4的數(shù)據(jù)，采用方差分析的計算方法進行進一步分析，計算結(jié)果如表5所示：

結(jié)果表明：因素A軸間距的檢驗F 值為 7.328， F0.1>7.328> F0.25，所以因素A軸間距在α = 0.25水平上對混合度的影響較為顯著;因素C轉(zhuǎn)軸轉(zhuǎn)速的檢驗F值為3.281，F(xiàn)0.1>3.281> F0.25，所以因素B 雙轉(zhuǎn)軸轉(zhuǎn)速在α = 0.25 水平上對混合時間的影響較顯著。因素B葉片間距的因素檢驗F 值為1< F0.25，所以因素B葉片間距對混合時間的影響不顯著。

4? 結(jié)論

本文采用了離散元的分析方法，針對臥式攪拌機的結(jié)構(gòu)及工藝參數(shù)對攪拌效率的影響進行了研究。通過分析軟件EDEM模擬顆?；旌舷到y(tǒng)，精確模擬散體顆?；旌戏植紶顟B(tài)，并對攪拌結(jié)構(gòu)參數(shù)與工藝參數(shù)匹配，以多因素正交實驗和極差方差分析為研究手段，評定了攪拌軸中心距，雙臂排布參數(shù)和葉片轉(zhuǎn)速因素在多水平下的混合度指標，并得到了以下結(jié)論：在三因素三水平L9 （ 33）正交實驗中，軸間距為1000mm，葉片間距610mm，雙軸轉(zhuǎn)速為80rpm的參數(shù)水平下混合度參數(shù)結(jié)果最佳;在攪拌結(jié)構(gòu)參數(shù)和工藝參數(shù)中，軸間距因素對混合度影響最大，雙軸轉(zhuǎn)速次之，攪拌軸葉片間距影響最小。表面合適的軸間距和恰當?shù)臄嚢杷俣葘嚢铏C混合性能影響明顯，合適的攪拌軸葉片間距也在一定程度上影響攪拌效率。