機械式高速混合機內(nèi)部多粒徑顆?；旌咸匦苑抡?/h1>
2018-04-25 01:24:08謝紅笑劉雪東尹傳忠溫傳美
中國粉體技術(shù) 2018年2期
謝紅笑，劉雪東，尹傳忠，馬 乾，溫傳美
（常州大學a.機械工程學院；b.江蘇省綠色過程裝備重點實驗室，江蘇常州 213000）
干燥粉體物料的混合廣泛應(yīng)用于化工、飼料、制藥、食品、材料（粉體冶金）等領(lǐng)域。2種或2種以上不同特性的物料，在機械或是其他外力的作用下相互進入到對方物料所在區(qū)域，物料顆粒的運動速度、方向、位置不斷發(fā)生變化，使各特性的物料顆粒在總體區(qū)域及各個子空間內(nèi)分布均勻，這樣的一個過程就是混合[1]。機械式混合具有混合快、混合均勻、機器操作方便、布置靈活、對環(huán)境無污染、維修簡單等突出優(yōu)勢[2-3]。攪拌槳與高速飛刀是高速混合機進行顆粒混合的核心構(gòu)件，研究顆粒的運動行為對于提高高速混合機混合質(zhì)量和混合效率具有重要意義[4-5]。
Alian等[6]開展了一系列顆?；旌蠈嶒灒ㄟ^改變初始裝填方式、物料裝填系數(shù)、轉(zhuǎn)速等條件，探究固體顆粒攪拌速率的變化規(guī)律和混合均勻度的變化規(guī)律。Veerakiet等[7]通過改變攪拌槳的槳葉數(shù)研究混合擴散系數(shù)的變化規(guī)律，進而提高混合機的混合性能。Halidan等[8]采用離散單元法（DEM）研究立式高速混合機內(nèi)顆粒的大小與密度對顆粒運動行為的影響，進而提高高速混合機的混合性能。國內(nèi)外學者的實驗或數(shù)值仿真重點多集中在高速混合機的混合效率以及混合均勻度，對高速混合機內(nèi)的顆粒運動規(guī)律探究較少，對多粒徑顆粒的混合研究更少。鑒于此，采用離散單元法對高速混合機內(nèi)部多粒徑顆粒的混合過程進行數(shù)值模擬，對探索顆粒運動特性的研究以及優(yōu)化高速混合機結(jié)構(gòu)設(shè)計具有一定的參考意義。
1 計算域模型及設(shè)置
1.1 模型建立
機械式高速混合機結(jié)構(gòu)示意圖如圖1所示。釜體呈圓筒形，為多層結(jié)構(gòu)，分內(nèi)外層、夾套和隔熱層；攪拌槳由驅(qū)動軸帶動旋轉(zhuǎn)，其驅(qū)動裝置由調(diào)速電機、傳動軸、皮帶輪等組成。
圖1 高速混合機結(jié)構(gòu)示意圖Fig.1 Schematic of high speed mixer
采用離散單元軟件EDEM[9-10]建立高速混合機仿真計算模型，如圖2所示。模型主要由攪拌槳、高速飛刀、釜體、擋板四部分組成。其中釜體直徑Φ=180 mm，高度H=180 mm；攪拌槳直徑db=172 mm；高速飛刀直徑dc=37 mm；擋板長度L=50 mm。
圖2 高速混合機仿真計算模型Fig.2 Simulated computational model of high speed mixer
1.2 控制方程
離散單元法是根據(jù)離散單元間的相互作用和牛頓第二定律，通過不斷循環(huán)計算獲得顆粒在每個時間步的位置信息和速度信息，同時能夠獲得顆粒間發(fā)生接觸碰撞時的接觸力，得到顆粒運動的詳細微觀數(shù)據(jù)信息，進而建立宏觀運動與微觀行為的聯(lián)系，更好地掌握顆粒運動規(guī)律[11]。圖3為離散單元法接觸模型，顆粒與顆粒間、顆粒與幾何體之間均采用振動方程進行模擬，接觸模型采用Hertz-Mindlin接觸模型，法向部分來自Hertz接觸理論，切向部分則來自Mindlin接觸理論[12]。
圖3 離散元方法接觸模型Fig.3 Contact model in DEM
顆粒間法向力和法向阻尼力為
式中：E*為等效彈性模量，MPa；R*為等效粒子半徑，m；Un為顆粒的法向重疊量，m；β為阻尼系數(shù)；kn為法向剛度N/m；m*為等效質(zhì)量，kg；為相對速度的法向分量值，m/s。
式中：E 為顆粒的彈性模量，MPa；μ 泊松比；Ri、Rj為球體接觸半徑，m；mi、mj為顆粒 i，j的等效質(zhì)量，kg。
顆粒間的切向力和切向阻尼力為
其中：
式中：Ut為切向重疊量，m；kt為切向剛度，N/m；G 為等效剪切模量，MPa為切向相對速度，m/s。
1.3 模擬參數(shù)的設(shè)置
DEM設(shè)置。采用模擬軟件EDEM 2.7.1，EDEM的相關(guān)參數(shù)設(shè)置[13]如表1所示?；旌蠈ο鬄橘|(zhì)量比為1∶1∶1∶1 的 4 種均勻球形亞克力顆粒組成的顆粒群,每種顆粒質(zhì)量為0.4 kg，顆粒直徑分別為2、3、4、5 mm（顆粒數(shù)量分別為 79 577、23 578、9 947、5 092），并且以小粒徑顆粒在下大粒徑顆粒在上的原則依次填充，混合顆粒群的粒徑分布如圖4所示?；旌线^程中，攪拌槳轉(zhuǎn)速為1 400 r/min；高速飛刀轉(zhuǎn)速為3 000 r/min，轉(zhuǎn)動方向與攪拌槳相同。
圖4 混合顆粒群的粒徑分布圖Fig.4 Particle size distribution of mixing particle group
表1 EDEM參數(shù)設(shè)置Tab.1 Parameters setting of EDEM
1.4 混合質(zhì)量評價
對整個高速混合機進行12×12×12的空間網(wǎng)格劃分,細分網(wǎng)格如圖5所示。以0.2 s為時間間隔將0～5 s不同時刻的各個網(wǎng)格中的顆粒總數(shù)以及每個網(wǎng)格中各種顆粒數(shù)量進行導(dǎo)出計算[14]。
顆粒靜止后的環(huán)形網(wǎng)格模型劃分如圖6所示。對每個網(wǎng)格區(qū)域內(nèi)的顆粒數(shù)目進行統(tǒng)計，得到各區(qū)域內(nèi)4種不同粒徑顆粒的數(shù)目，并計算其絕對平均粒徑 d[15]，計算公式為
式中：di為顆粒 i的直徑，m；mi為顆粒 i的質(zhì)量，kg。為了表征顆粒靜止后各個環(huán)形網(wǎng)格區(qū)域內(nèi)顆粒粒度大小的差異，引入相對平均粒度，它是各區(qū)域內(nèi)顆粒的絕對平均粒徑與整個高速混合機內(nèi)顆粒的絕對平均粒徑的比值。
圖5 細分網(wǎng)格Fig.5 The subdivided mesh
圖6 顆粒靜止后的環(huán)形網(wǎng)格模型Fig.6 Annular mesh model after particle mixing process
為量化高速混合機內(nèi)部顆粒的偏析分布，統(tǒng)計如圖6b所劃分的環(huán)形網(wǎng)格區(qū)域內(nèi)大顆粒的體積分數(shù)。使用不同網(wǎng)格區(qū)域內(nèi)大顆粒所占體積分數(shù)之間的差值，表征高速混合機內(nèi)各網(wǎng)格區(qū)域顆粒是否形成偏析結(jié)構(gòu)。針對偏析結(jié)構(gòu)，定義偏析指數(shù)SI：
即用大顆粒在第n個網(wǎng)格中的體積分數(shù)φn，與大顆粒在整個高速混合機內(nèi)顆粒群的體積分數(shù)φ0之間差值的絕對值，來表示偏析程度。SI值越低，表明偏析程度越小。
2 結(jié)果與討論
2.1 顆?；旌线^程
2.1.1 顆?；旌暇鶆蚨?/h4>
圖7為離散系數(shù)隨時間的變化規(guī)律。由圖可知，離散系數(shù)隨時間不斷減小，混合均勻度不斷提高，混合時間為3 s時，離散系數(shù)達到0.15，3 s后離散系數(shù)基本保持不變，顆粒混合均勻。
圖7 離散系數(shù)隨時間的變化規(guī)律Fig.7 Relative standard deviation of mixing with time
2.1.2 顆粒運動軌跡
圖8為高速混合機內(nèi)顆粒在0、0.5、1、3 s 4個不同時刻的混合狀況，0 s時，4種顆粒幾乎是完全分離的狀態(tài)，隨著攪拌槳與高速飛刀的轉(zhuǎn)動，4種顆粒逐漸混合，最終獲得一個均勻度高的混合物。
圖8 不同時刻顆粒混合狀況Fig.8 Movement condition of particles
圖9為直徑為4 mm和5 mm顆粒的速度矢量圖。由圖可知，高速旋轉(zhuǎn)的攪拌槳借助表面與物料的磨擦力及側(cè)面對物料的推力使物料沿攪拌槳切向運動，由于離心力的作用，又被拋向釜體內(nèi)壁，并且沿壁面上升。一部分顆粒運動至高速飛刀位置被分散，另外一部分顆粒升高到一定高度時，由于重力作用又落回到攪拌槳中心，接著又被拋起，這種上升運動與切向運動結(jié)合，使物料處于連續(xù)的螺旋狀上下運動狀態(tài)。并且整個過程中伴隨著碰撞、折射、反彈等劇烈運動。
圖9 直徑為4 mm和5 mm顆粒速度矢量圖Fig.9 Velocity vector diagram of 4 and 5 millimeter particles
2.1.3 顆粒粒度分布及偏析行為
顆粒在不同區(qū)域內(nèi)的相對平均粒徑如圖10所示。x＝0處為高速混合機對稱中心，將各數(shù)值與1進行比較，其中，高速飛刀一側(cè)自邊緣往中心呈現(xiàn)出由小變大再變小的粒度分布，擋板一側(cè)呈現(xiàn)由小變大再變小的粒度分布，顆粒相對平均粒徑波動范圍在-0.05～0.05之間，在顆粒靜止過程中混合機邊緣出現(xiàn)小顆粒增多現(xiàn)象。
圖10 顆粒在不同區(qū)域的相對平均粒徑Fig.10 Relative to the average particle size in different regions
圖11為顆粒在環(huán)形網(wǎng)格內(nèi)的偏析指數(shù)分布，其值在0～0.07之間波動，相對于高速飛刀一側(cè)，擋板一側(cè)偏析指數(shù)波動稍大，并且位于高速混合機邊緣部分偏析指數(shù)SI增大，出現(xiàn)少量小顆粒偏析現(xiàn)象。靜止過程中，高速混合機擋板位置的特殊構(gòu)造以及大小顆粒各異的流動行為差異導(dǎo)致了顆粒在擋板一側(cè)和高速飛刀一側(cè)之間存在不同的粒度偏析結(jié)構(gòu)。
圖11 顆粒在環(huán)形網(wǎng)格內(nèi)的偏析指數(shù)分布Fig.11 The segregation index of particles in annular mesh
2.2 顆粒受力與速度分析
利用EDEM的后處理工具，將模型進行分層網(wǎng)格劃分，如圖12所示。研究不同高度和不同粒徑顆粒的運動和受力情況，將數(shù)據(jù)進行分析整理得到顆粒的合力和速度的變化規(guī)律。
圖12 分層網(wǎng)格Fig.12 Layered mesh
研究高速混合機內(nèi)部顆粒的運動學特性分析對提高高速混合機性能具有重要意義。圖13為顆粒平均速度與受力隨高度的變化規(guī)律。顆粒運動的平均速度與受力的變化規(guī)律大致相同，顆粒的平均速度與受力在釜體底部最大，隨著高度的升高，顆粒的速度與受力逐漸減小，并且，越高層顆粒速度與受力降低越慢。
圖13 顆粒平均速度和受力隨高度的變化規(guī)律Fig.13 Average velocity and force of particles with height
圖14為顆粒受力與速度隨粒徑的變化規(guī)律。由圖可知，顆粒在運動過程中受力（重力除外）和速度的變化規(guī)律大致相同。隨著顆粒尺寸的增大，顆粒的運動速度與受力逐漸增大。原因是顆粒的直徑越大，體積越大，同一時刻可以與更多的顆粒發(fā)生碰撞運動，導(dǎo)致大顆粒受到的力較大。高速混合機內(nèi)，顆粒速度隨著粒徑的增大，顆粒速度變大。這主要是由兩方面原因造成：一方面，在封閉的高速混合機內(nèi)單個大尺寸顆粒，所占空間較大，與更多顆粒發(fā)生碰撞的概率也大，周圍的小尺寸顆粒在運動過程中受到大顆粒的阻礙，導(dǎo)致速度降低；另一方面，在混合過程中，攪拌槳與高速飛刀的機械能轉(zhuǎn)化為顆粒的動能，如圖15所示，大尺寸顆粒獲得更高的能量，使得大尺寸顆粒運動速度較快。
圖14 顆粒受力與速度隨粒徑的變化規(guī)律Fig.14 The force and velocity of particles with diameter
圖15 法向力、切向力和總能量隨粒徑的變化規(guī)律Fig.15 Normal force,Tangential force and total energy of particles with diameter
圖15為切向力、法向力和總能量隨粒徑的變化規(guī)律。隨著顆粒直徑的變大，顆粒法向力與切向力逐漸增大。顆粒法向力大于切向力，顆粒尺寸越大，法向力與切向力之間差距越大。原因是高速轉(zhuǎn)動的攪拌槳使得顆粒受到離心力作用，使得顆粒被拋向釜體內(nèi)壁，顆粒之間相互碰撞、擠壓，顆粒尺寸越大，碰撞、擠壓越明顯。攪拌槳側(cè)面對顆粒的推力使得顆粒沿切線方向快速轉(zhuǎn)動，轉(zhuǎn)動過程中阻力較小。
3 結(jié)論
采用離散單元法對機械式高速混合機內(nèi)部4種不同粒徑球形亞克力顆粒組成的顆粒群的混合過程進行數(shù)值模擬，得到如下結(jié)論：
1）顆粒在攪拌槳與高速飛刀共同作用下，呈現(xiàn)連續(xù)的螺旋上下運動狀態(tài)，中心形成漩渦；
2）高速混合機的特殊構(gòu)造以及大小顆粒各異的流動行為導(dǎo)致顆粒在擋板一側(cè)和高速飛刀一側(cè)之間存在不同的粒度偏析結(jié)構(gòu)，靜止過程中邊緣出現(xiàn)少量小顆粒偏析現(xiàn)象；
3）顆?；旌现饕性诟咚倩旌蠙C中下部，顆粒受力與速度隨著高度的增加不斷減小，并且越高層顆粒受力與速度降低越慢；
4）顆粒在混合過程中，顆粒受力與顆粒速度隨著顆粒尺寸的不斷增大而增大；
5）高速旋轉(zhuǎn)的攪拌槳使顆粒受較大離心力作用，顆粒法向力大于切向力，并且顆粒尺寸越大，法向力與切向力差距越大。
參考文獻（References）：
[1]葉濤.多組分粉體混合過程的理論分析與實驗研究[D].武漢：武漢理工大學，2009.
[2]ZHOU Y C，YU A B，STEWART R L，et al.Microdynamic analysis of the particle flow in a cylindrical bladed mixer[J].Chemical Engineering Science，2004，59（6）：1343-1364.
[3]陳俊英，李紅偉.化工行業(yè)常用攪拌器研究熱點[J].現(xiàn)代化工，2011，31（10）：12-16.
[4]王港，黃銳，陳曉媛，等.高速混合機的應(yīng)用及研究進展[J].中國塑料，2001，15（7）：11-14.
[5]CHANDRATILLEKE G R，YU A B，STEWART R L，et al.Effects of blade rake angle and gap on particle mixing in a cylindrical mixer[J].Powder Technology，2009，193（3）：303-311.
[6]ALIAN M，F(xiàn)ARHAD E M，SIMANT R，et al.Using discrete element method to analyze the mixing of the solid particles in a slant cone mixer[J].Chemical Engineering Research and Design，2015，93（1）：318-329.
[7]BOONKANOKWONG V，REMY B，KHINAST J G，et al.The effect of the number of impeller blades on granular flow in a bladed mixer[J].Powder Technology，2016，302（11）：333-349.
[8]HALIDAN M，CHANDRATILLEKE G R，CHAN S L I，et al.Prediction of the mixing behaviour of binary mixtures of particles in a bladed mixer[J].Chemical Engineering Science，2014，120（12）：37-48.
[9]周志勇.基于DEM的強力混合機混勻過程模擬及結(jié)構(gòu)優(yōu)化[D].湘潭：湘潭大學，2014.
[10]胡國明.顆粒系統(tǒng)的離散元素法分析仿真[M].武漢：武漢理工大學出版社，2010：1-3.
[11]SAKAI M，SHIGETO Y，BASINSKAS G，et al.Discrete element simulation for the evaluation of solid mixing in an industrial blender[J].Chemical Engineering Journal，2015，279（11）：821-839.
[12]CLEARY P W，SINNOTT M D.Assessing mixing characteristics of particle-mixing and granulation devices[J].Particuology，2008，6（6）：419-444.
[13]LIU W M，XU J，LIU X D.Numerical study on collision characteristics for non-spherical particles in venturi powder ejector[J].Vacuum，2016，131（5）：285-292.
[14]MUZZIO F J，LLUSA M，GOODRIDGE C L，et al.Evaluating the performance of a ribbon blender[J].Powder Technology，2008，186（3）：247-254.
[15]朱明華，杜屏，劉建波.高爐爐頂料罐裝料過程計算及分析[C]//全國大高爐煉鐵學術(shù)年會，2013：78-83.
猜你喜歡
犁刀式混合機飛刀軸結(jié)構(gòu)設(shè)計及臨界轉(zhuǎn)速計算*機械研究與應(yīng)用(2022年5期)2022-11-28 06:26:42 李飛刀的飛刀（外一篇）小小說月刊·下半月(2021年12期)2021-12-25 12:00:37 藝高人膽大西江月(2018年6期)2018-11-15 02:28:20 與鳥相撞飛機受力幾何中學生數(shù)理化·八年級物理人教版(2018年5期)2018-06-21 08:01:38 改善60Si2MnA彈簧鋼小方坯中心碳偏析的研究上海金屬(2016年4期)2016-11-23 05:39:12 受力分析的三個“囑托”中學生數(shù)理化·高一版(2016年10期)2016-03-01 03:45:54 原鋁旋轉(zhuǎn)偏析法凈化工藝大連工業(yè)大學學報(2015年4期)2015-12-11 04:06:50 改善82B 小方坯碳偏析的連鑄工藝研究新疆鋼鐵(2015年1期)2015-11-07 03:11:03 SWRH82B小方坯中心偏析的改善上海金屬(2014年1期)2014-12-18 06:52:12 底排藥受力載荷及其分布規(guī)律火炸藥學報(2014年3期)2014-03-20 13:17:44
謝紅笑，劉雪東，尹傳忠，馬 乾，溫傳美
（常州大學a.機械工程學院；b.江蘇省綠色過程裝備重點實驗室，江蘇常州 213000）
干燥粉體物料的混合廣泛應(yīng)用于化工、飼料、制藥、食品、材料（粉體冶金）等領(lǐng)域。2種或2種以上不同特性的物料，在機械或是其他外力的作用下相互進入到對方物料所在區(qū)域，物料顆粒的運動速度、方向、位置不斷發(fā)生變化，使各特性的物料顆粒在總體區(qū)域及各個子空間內(nèi)分布均勻，這樣的一個過程就是混合[1]。機械式混合具有混合快、混合均勻、機器操作方便、布置靈活、對環(huán)境無污染、維修簡單等突出優(yōu)勢[2-3]。攪拌槳與高速飛刀是高速混合機進行顆粒混合的核心構(gòu)件，研究顆粒的運動行為對于提高高速混合機混合質(zhì)量和混合效率具有重要意義[4-5]。
Alian等[6]開展了一系列顆?；旌蠈嶒灒ㄟ^改變初始裝填方式、物料裝填系數(shù)、轉(zhuǎn)速等條件，探究固體顆粒攪拌速率的變化規(guī)律和混合均勻度的變化規(guī)律。Veerakiet等[7]通過改變攪拌槳的槳葉數(shù)研究混合擴散系數(shù)的變化規(guī)律，進而提高混合機的混合性能。Halidan等[8]采用離散單元法（DEM）研究立式高速混合機內(nèi)顆粒的大小與密度對顆粒運動行為的影響，進而提高高速混合機的混合性能。國內(nèi)外學者的實驗或數(shù)值仿真重點多集中在高速混合機的混合效率以及混合均勻度，對高速混合機內(nèi)的顆粒運動規(guī)律探究較少，對多粒徑顆粒的混合研究更少。鑒于此，采用離散單元法對高速混合機內(nèi)部多粒徑顆粒的混合過程進行數(shù)值模擬，對探索顆粒運動特性的研究以及優(yōu)化高速混合機結(jié)構(gòu)設(shè)計具有一定的參考意義。
1 計算域模型及設(shè)置
1.1 模型建立
機械式高速混合機結(jié)構(gòu)示意圖如圖1所示。釜體呈圓筒形，為多層結(jié)構(gòu)，分內(nèi)外層、夾套和隔熱層；攪拌槳由驅(qū)動軸帶動旋轉(zhuǎn)，其驅(qū)動裝置由調(diào)速電機、傳動軸、皮帶輪等組成。
圖1 高速混合機結(jié)構(gòu)示意圖Fig.1 Schematic of high speed mixer
采用離散單元軟件EDEM[9-10]建立高速混合機仿真計算模型，如圖2所示。模型主要由攪拌槳、高速飛刀、釜體、擋板四部分組成。其中釜體直徑Φ=180 mm，高度H=180 mm；攪拌槳直徑db=172 mm；高速飛刀直徑dc=37 mm；擋板長度L=50 mm。
圖2 高速混合機仿真計算模型Fig.2 Simulated computational model of high speed mixer
1.2 控制方程
離散單元法是根據(jù)離散單元間的相互作用和牛頓第二定律，通過不斷循環(huán)計算獲得顆粒在每個時間步的位置信息和速度信息，同時能夠獲得顆粒間發(fā)生接觸碰撞時的接觸力，得到顆粒運動的詳細微觀數(shù)據(jù)信息，進而建立宏觀運動與微觀行為的聯(lián)系，更好地掌握顆粒運動規(guī)律[11]。圖3為離散單元法接觸模型，顆粒與顆粒間、顆粒與幾何體之間均采用振動方程進行模擬，接觸模型采用Hertz-Mindlin接觸模型，法向部分來自Hertz接觸理論，切向部分則來自Mindlin接觸理論[12]。
圖3 離散元方法接觸模型Fig.3 Contact model in DEM
顆粒間法向力和法向阻尼力為
式中：E*為等效彈性模量，MPa；R*為等效粒子半徑，m；Un為顆粒的法向重疊量，m；β為阻尼系數(shù)；kn為法向剛度N/m；m*為等效質(zhì)量，kg；為相對速度的法向分量值，m/s。
式中：E 為顆粒的彈性模量，MPa；μ 泊松比；Ri、Rj為球體接觸半徑，m；mi、mj為顆粒 i，j的等效質(zhì)量，kg。
顆粒間的切向力和切向阻尼力為
其中：
式中：Ut為切向重疊量，m；kt為切向剛度，N/m；G 為等效剪切模量，MPa為切向相對速度，m/s。
1.3 模擬參數(shù)的設(shè)置
DEM設(shè)置。采用模擬軟件EDEM 2.7.1，EDEM的相關(guān)參數(shù)設(shè)置[13]如表1所示?；旌蠈ο鬄橘|(zhì)量比為1∶1∶1∶1 的 4 種均勻球形亞克力顆粒組成的顆粒群,每種顆粒質(zhì)量為0.4 kg，顆粒直徑分別為2、3、4、5 mm（顆粒數(shù)量分別為 79 577、23 578、9 947、5 092），并且以小粒徑顆粒在下大粒徑顆粒在上的原則依次填充，混合顆粒群的粒徑分布如圖4所示?；旌线^程中，攪拌槳轉(zhuǎn)速為1 400 r/min；高速飛刀轉(zhuǎn)速為3 000 r/min，轉(zhuǎn)動方向與攪拌槳相同。
圖4 混合顆粒群的粒徑分布圖Fig.4 Particle size distribution of mixing particle group
表1 EDEM參數(shù)設(shè)置Tab.1 Parameters setting of EDEM
1.4 混合質(zhì)量評價
對整個高速混合機進行12×12×12的空間網(wǎng)格劃分,細分網(wǎng)格如圖5所示。以0.2 s為時間間隔將0～5 s不同時刻的各個網(wǎng)格中的顆粒總數(shù)以及每個網(wǎng)格中各種顆粒數(shù)量進行導(dǎo)出計算[14]。
顆粒靜止后的環(huán)形網(wǎng)格模型劃分如圖6所示。對每個網(wǎng)格區(qū)域內(nèi)的顆粒數(shù)目進行統(tǒng)計，得到各區(qū)域內(nèi)4種不同粒徑顆粒的數(shù)目，并計算其絕對平均粒徑 d[15]，計算公式為
式中：di為顆粒 i的直徑，m；mi為顆粒 i的質(zhì)量，kg。為了表征顆粒靜止后各個環(huán)形網(wǎng)格區(qū)域內(nèi)顆粒粒度大小的差異，引入相對平均粒度，它是各區(qū)域內(nèi)顆粒的絕對平均粒徑與整個高速混合機內(nèi)顆粒的絕對平均粒徑的比值。
圖5 細分網(wǎng)格Fig.5 The subdivided mesh
圖6 顆粒靜止后的環(huán)形網(wǎng)格模型Fig.6 Annular mesh model after particle mixing process
為量化高速混合機內(nèi)部顆粒的偏析分布，統(tǒng)計如圖6b所劃分的環(huán)形網(wǎng)格區(qū)域內(nèi)大顆粒的體積分數(shù)。使用不同網(wǎng)格區(qū)域內(nèi)大顆粒所占體積分數(shù)之間的差值，表征高速混合機內(nèi)各網(wǎng)格區(qū)域顆粒是否形成偏析結(jié)構(gòu)。針對偏析結(jié)構(gòu)，定義偏析指數(shù)SI：
即用大顆粒在第n個網(wǎng)格中的體積分數(shù)φn，與大顆粒在整個高速混合機內(nèi)顆粒群的體積分數(shù)φ0之間差值的絕對值，來表示偏析程度。SI值越低，表明偏析程度越小。
2 結(jié)果與討論
2.1 顆?；旌线^程
2.1.1 顆?；旌暇鶆蚨?/h4>
圖7為離散系數(shù)隨時間的變化規(guī)律。由圖可知，離散系數(shù)隨時間不斷減小，混合均勻度不斷提高，混合時間為3 s時，離散系數(shù)達到0.15，3 s后離散系數(shù)基本保持不變，顆粒混合均勻。
圖7 離散系數(shù)隨時間的變化規(guī)律Fig.7 Relative standard deviation of mixing with time
2.1.2 顆粒運動軌跡
圖8為高速混合機內(nèi)顆粒在0、0.5、1、3 s 4個不同時刻的混合狀況，0 s時，4種顆粒幾乎是完全分離的狀態(tài)，隨著攪拌槳與高速飛刀的轉(zhuǎn)動，4種顆粒逐漸混合，最終獲得一個均勻度高的混合物。
圖8 不同時刻顆粒混合狀況Fig.8 Movement condition of particles
圖9為直徑為4 mm和5 mm顆粒的速度矢量圖。由圖可知，高速旋轉(zhuǎn)的攪拌槳借助表面與物料的磨擦力及側(cè)面對物料的推力使物料沿攪拌槳切向運動，由于離心力的作用，又被拋向釜體內(nèi)壁，并且沿壁面上升。一部分顆粒運動至高速飛刀位置被分散，另外一部分顆粒升高到一定高度時，由于重力作用又落回到攪拌槳中心，接著又被拋起，這種上升運動與切向運動結(jié)合，使物料處于連續(xù)的螺旋狀上下運動狀態(tài)。并且整個過程中伴隨著碰撞、折射、反彈等劇烈運動。
圖9 直徑為4 mm和5 mm顆粒速度矢量圖Fig.9 Velocity vector diagram of 4 and 5 millimeter particles
2.1.3 顆粒粒度分布及偏析行為
顆粒在不同區(qū)域內(nèi)的相對平均粒徑如圖10所示。x＝0處為高速混合機對稱中心，將各數(shù)值與1進行比較，其中，高速飛刀一側(cè)自邊緣往中心呈現(xiàn)出由小變大再變小的粒度分布，擋板一側(cè)呈現(xiàn)由小變大再變小的粒度分布，顆粒相對平均粒徑波動范圍在-0.05～0.05之間，在顆粒靜止過程中混合機邊緣出現(xiàn)小顆粒增多現(xiàn)象。
圖10 顆粒在不同區(qū)域的相對平均粒徑Fig.10 Relative to the average particle size in different regions
圖11為顆粒在環(huán)形網(wǎng)格內(nèi)的偏析指數(shù)分布，其值在0～0.07之間波動，相對于高速飛刀一側(cè)，擋板一側(cè)偏析指數(shù)波動稍大，并且位于高速混合機邊緣部分偏析指數(shù)SI增大，出現(xiàn)少量小顆粒偏析現(xiàn)象。靜止過程中，高速混合機擋板位置的特殊構(gòu)造以及大小顆粒各異的流動行為差異導(dǎo)致了顆粒在擋板一側(cè)和高速飛刀一側(cè)之間存在不同的粒度偏析結(jié)構(gòu)。
圖11 顆粒在環(huán)形網(wǎng)格內(nèi)的偏析指數(shù)分布Fig.11 The segregation index of particles in annular mesh
2.2 顆粒受力與速度分析
利用EDEM的后處理工具，將模型進行分層網(wǎng)格劃分，如圖12所示。研究不同高度和不同粒徑顆粒的運動和受力情況，將數(shù)據(jù)進行分析整理得到顆粒的合力和速度的變化規(guī)律。
圖12 分層網(wǎng)格Fig.12 Layered mesh
研究高速混合機內(nèi)部顆粒的運動學特性分析對提高高速混合機性能具有重要意義。圖13為顆粒平均速度與受力隨高度的變化規(guī)律。顆粒運動的平均速度與受力的變化規(guī)律大致相同，顆粒的平均速度與受力在釜體底部最大，隨著高度的升高，顆粒的速度與受力逐漸減小，并且，越高層顆粒速度與受力降低越慢。
圖13 顆粒平均速度和受力隨高度的變化規(guī)律Fig.13 Average velocity and force of particles with height
圖14為顆粒受力與速度隨粒徑的變化規(guī)律。由圖可知，顆粒在運動過程中受力（重力除外）和速度的變化規(guī)律大致相同。隨著顆粒尺寸的增大，顆粒的運動速度與受力逐漸增大。原因是顆粒的直徑越大，體積越大，同一時刻可以與更多的顆粒發(fā)生碰撞運動，導(dǎo)致大顆粒受到的力較大。高速混合機內(nèi)，顆粒速度隨著粒徑的增大，顆粒速度變大。這主要是由兩方面原因造成：一方面，在封閉的高速混合機內(nèi)單個大尺寸顆粒，所占空間較大，與更多顆粒發(fā)生碰撞的概率也大，周圍的小尺寸顆粒在運動過程中受到大顆粒的阻礙，導(dǎo)致速度降低；另一方面，在混合過程中，攪拌槳與高速飛刀的機械能轉(zhuǎn)化為顆粒的動能，如圖15所示，大尺寸顆粒獲得更高的能量，使得大尺寸顆粒運動速度較快。
圖14 顆粒受力與速度隨粒徑的變化規(guī)律Fig.14 The force and velocity of particles with diameter
圖15 法向力、切向力和總能量隨粒徑的變化規(guī)律Fig.15 Normal force,Tangential force and total energy of particles with diameter
圖15為切向力、法向力和總能量隨粒徑的變化規(guī)律。隨著顆粒直徑的變大，顆粒法向力與切向力逐漸增大。顆粒法向力大于切向力，顆粒尺寸越大，法向力與切向力之間差距越大。原因是高速轉(zhuǎn)動的攪拌槳使得顆粒受到離心力作用，使得顆粒被拋向釜體內(nèi)壁，顆粒之間相互碰撞、擠壓，顆粒尺寸越大，碰撞、擠壓越明顯。攪拌槳側(cè)面對顆粒的推力使得顆粒沿切線方向快速轉(zhuǎn)動，轉(zhuǎn)動過程中阻力較小。
3 結(jié)論
采用離散單元法對機械式高速混合機內(nèi)部4種不同粒徑球形亞克力顆粒組成的顆粒群的混合過程進行數(shù)值模擬，得到如下結(jié)論：
1）顆粒在攪拌槳與高速飛刀共同作用下，呈現(xiàn)連續(xù)的螺旋上下運動狀態(tài)，中心形成漩渦；
2）高速混合機的特殊構(gòu)造以及大小顆粒各異的流動行為導(dǎo)致顆粒在擋板一側(cè)和高速飛刀一側(cè)之間存在不同的粒度偏析結(jié)構(gòu)，靜止過程中邊緣出現(xiàn)少量小顆粒偏析現(xiàn)象；
3）顆?；旌现饕性诟咚倩旌蠙C中下部，顆粒受力與速度隨著高度的增加不斷減小，并且越高層顆粒受力與速度降低越慢；
4）顆粒在混合過程中，顆粒受力與顆粒速度隨著顆粒尺寸的不斷增大而增大；
5）高速旋轉(zhuǎn)的攪拌槳使顆粒受較大離心力作用，顆粒法向力大于切向力，并且顆粒尺寸越大，法向力與切向力差距越大。
參考文獻（References）：
[1]葉濤.多組分粉體混合過程的理論分析與實驗研究[D].武漢：武漢理工大學，2009.
[2]ZHOU Y C，YU A B，STEWART R L，et al.Microdynamic analysis of the particle flow in a cylindrical bladed mixer[J].Chemical Engineering Science，2004，59（6）：1343-1364.
[3]陳俊英，李紅偉.化工行業(yè)常用攪拌器研究熱點[J].現(xiàn)代化工，2011，31（10）：12-16.
[4]王港，黃銳，陳曉媛，等.高速混合機的應(yīng)用及研究進展[J].中國塑料，2001，15（7）：11-14.
[5]CHANDRATILLEKE G R，YU A B，STEWART R L，et al.Effects of blade rake angle and gap on particle mixing in a cylindrical mixer[J].Powder Technology，2009，193（3）：303-311.
[6]ALIAN M，F(xiàn)ARHAD E M，SIMANT R，et al.Using discrete element method to analyze the mixing of the solid particles in a slant cone mixer[J].Chemical Engineering Research and Design，2015，93（1）：318-329.
[7]BOONKANOKWONG V，REMY B，KHINAST J G，et al.The effect of the number of impeller blades on granular flow in a bladed mixer[J].Powder Technology，2016，302（11）：333-349.
[8]HALIDAN M，CHANDRATILLEKE G R，CHAN S L I，et al.Prediction of the mixing behaviour of binary mixtures of particles in a bladed mixer[J].Chemical Engineering Science，2014，120（12）：37-48.
[9]周志勇.基于DEM的強力混合機混勻過程模擬及結(jié)構(gòu)優(yōu)化[D].湘潭：湘潭大學，2014.
[10]胡國明.顆粒系統(tǒng)的離散元素法分析仿真[M].武漢：武漢理工大學出版社，2010：1-3.
[11]SAKAI M，SHIGETO Y，BASINSKAS G，et al.Discrete element simulation for the evaluation of solid mixing in an industrial blender[J].Chemical Engineering Journal，2015，279（11）：821-839.
[12]CLEARY P W，SINNOTT M D.Assessing mixing characteristics of particle-mixing and granulation devices[J].Particuology，2008，6（6）：419-444.
[13]LIU W M，XU J，LIU X D.Numerical study on collision characteristics for non-spherical particles in venturi powder ejector[J].Vacuum，2016，131（5）：285-292.
[14]MUZZIO F J，LLUSA M，GOODRIDGE C L，et al.Evaluating the performance of a ribbon blender[J].Powder Technology，2008，186（3）：247-254.
[15]朱明華，杜屏，劉建波.高爐爐頂料罐裝料過程計算及分析[C]//全國大高爐煉鐵學術(shù)年會，2013：78-83.
猜你喜歡
機械研究與應(yīng)用(2022年5期)2022-11-28 06:26:42
小小說月刊·下半月(2021年12期)2021-12-25 12:00:37
西江月(2018年6期)2018-11-15 02:28:20
中學生數(shù)理化·八年級物理人教版(2018年5期)2018-06-21 08:01:38
上海金屬(2016年4期)2016-11-23 05:39:12
中學生數(shù)理化·高一版(2016年10期)2016-03-01 03:45:54
大連工業(yè)大學學報(2015年4期)2015-12-11 04:06:50
新疆鋼鐵(2015年1期)2015-11-07 03:11:03
上海金屬(2014年1期)2014-12-18 06:52:12
火炸藥學報(2014年3期)2014-03-20 13:17:44
