林赫超 何延?xùn)| 朱向哲

摘 ?????要： 建立了三維行星式攪拌器的有限元模型，應(yīng)用Polyflow軟件對流場進(jìn)行模擬，對不同轉(zhuǎn)速下的攪拌器內(nèi)糖漿的混合流動進(jìn)行了對比分析，同時對中心區(qū)和銷附近區(qū)域的混合效果進(jìn)行對比。結(jié)果表明，攪拌器不存在流動混合的死區(qū)，且銷附近區(qū)域的混合能力比中心區(qū)的混合能力更強(qiáng)，流場內(nèi)的粒子始終受到周期性的折疊和拉伸作用;轉(zhuǎn)速越大，攪拌器內(nèi)的流體混合效果越好。

關(guān) ?鍵 ?詞：行星式攪拌器;混合效率;剪切速率;拉伸指數(shù)

中圖分類號：TQ051.72 ?????文獻(xiàn)標(biāo)識碼： A ??????文章編號： 1671-0460（2019）01-0196-06

Abstract： The three-dimensional finite element model of planetary mixer was established and the flow field was simulated by using Polyflow software， the mixing flow of syrup in the mixer at different speeds was compared and analyzed. Besides， the effect of mixing in the central area and the region near the pin was compared. The results showed that there was not a dead zone of flow mixing in the mixer and the mixing capacity in the region near the pin was stronger than that in the center. The particles in the flow field were always subjected to periodic folding and stretching. The?faster?the?pin?speed was，?the better the mixing effect in the mixer was.

Key words： Planetary mixer; Mixing efficiency; Shear rate; Stretch logarithm

行星攪拌器廣泛應(yīng)用于流體加工，常用于食品、醫(yī)藥、冶金、化工等加工過程[1]。內(nèi)部的攪拌釜不僅能繞攪拌器中心軸進(jìn)行公轉(zhuǎn)，同時它自身還進(jìn)行自轉(zhuǎn)，促使混合物在極短時間內(nèi)達(dá)到一定的混合效果[2]。運(yùn)行過程中，還可以根據(jù)加工工藝和流體黏度的不同來更換傳動比，以此得到理想的主軸轉(zhuǎn)速[3]。自轉(zhuǎn)較公轉(zhuǎn)的快慢，影響了混合的均勻性與激烈程度[4]。為設(shè)計(jì)出性能更加優(yōu)良的攪拌器，眾多學(xué)者從不同角度對行星攪拌器進(jìn)行了研究。T.Jongen[5]對行星式面團(tuán)攪拌器內(nèi)部流動進(jìn)行數(shù)值分析，主要描述面團(tuán)在攪拌器中受到的各種擠壓變形。RK Connelly[6，7]對雙相葉片混合器及三維行星攪拌器中的牛頓流體進(jìn)行了混合分析。朱霖龍[8]等對長幅內(nèi)旋輪線行星攪拌器進(jìn)行研究，對比分析了新型攪拌器與外嚙合行星攪拌器的葉片攪拌軌跡和線速度梯度。楊伶[9]等對行星攪拌釜內(nèi)的流場進(jìn)行了三維數(shù)值模擬，發(fā)現(xiàn)最高剪切速率主要位于兩個攪拌槳之間。王超明[10]等對比了行星輪式和普通式攪拌器的流場速度，發(fā)現(xiàn)行星輪式攪拌器比普通式攪拌器的攪拌效果更優(yōu)。

本文采用有限元法、利用Polyflow軟件，對行星式混合器內(nèi)的流場進(jìn)行分析，計(jì)算了不同轉(zhuǎn)速比下混合器中流體速度場、混合指數(shù)、剪切速率、拉伸指數(shù)、瞬時混合效率等多種參數(shù)。同時對行星攪拌器中心區(qū)和銷附近區(qū)域的動態(tài)混合特性也進(jìn)行了深入的研究，為后續(xù)研究不同運(yùn)動參數(shù)對攪拌器混合能力的影響奠定了基礎(chǔ)。

1 ?有限元模型

1.1 ?數(shù)學(xué)模型

有限元模型采用如下假設(shè)：（1）牛頓黏度為常數(shù);（2）流體等溫、不可壓縮;（3）慣性力和重力忽略不計(jì); （4）流道壁面無滑移[11]。

三維行星銷釘攪拌器內(nèi)流體為糖漿，49 ℃下黏度為5 400 cP，密度為1 409 kg/m3。攪拌器內(nèi)牛頓流體的連續(xù)性方程、動量方程和本構(gòu)方程分別為[12]：

1.2 ?物理模型

三維行星銷釘攪拌器幾何模型如圖1所示，定位銷釘為黑色，轉(zhuǎn)動銷釘為灰色。尺寸如下：攪拌器半徑為15 mm，高度為30 mm，轉(zhuǎn)動銷釘直徑為1.7 mm，固定銷釘直徑為1.5 mm。為模擬流場瞬時變化，在Polyflow中采用網(wǎng)格疊加技術(shù)（MST），流域軸向分成60等份，節(jié)點(diǎn)數(shù)為98 460個，流域有限元模型如圖2所示。為了簡化銷釘?shù)倪\(yùn)動形式，在Polyflow模擬時采用旋轉(zhuǎn)參考框架。墻壁和固定銷釘相對“靜止”，都設(shè)置成順時針方向，旋轉(zhuǎn)銷釘繞著旋轉(zhuǎn)中心逆時針旋轉(zhuǎn)，通過這種方法既能實(shí)現(xiàn)外旋輪線形式的運(yùn)動軌跡又簡化了模型，使計(jì)算簡便。在有限元計(jì)算時，瞬態(tài)流場計(jì)算采用的積分方法為Crank-Nicolson，坐標(biāo)采用線性插值，速度采用mini-element插值方法，壓力為線性插值，收斂精度為1×10-4。動銷釘?shù)耐庑喚€如圖3所示，外旋輪線的參數(shù)方程如下：

2 ?流場混合強(qiáng)度評價參數(shù)

為了評估分散混合效率，采用混合指數(shù)，定義如下：

3 ?計(jì)算結(jié)果

3.1 ?速度分布

在攪拌器中選取3個監(jiān)測點(diǎn)——中心點(diǎn)p1、靜銷附近的點(diǎn)p2、動銷附近的點(diǎn)p3，如圖4所示。圖5是當(dāng)靜銷釘與動銷釘?shù)霓D(zhuǎn)速比為4：3、8：6和12：9時3個點(diǎn)的速度分布。由圖可知，在同一轉(zhuǎn)速下， p1點(diǎn)速度變化小且穩(wěn)定在一定的區(qū)間內(nèi)，但p1點(diǎn)的速度始終不為0，說明中心區(qū)域不存在糖漿的停滯現(xiàn)象，它是始終流動和混合的。p2、p3點(diǎn)的速度分布變化較大，這是由于銷在轉(zhuǎn)動的過程中其附近區(qū)域的流體受到反復(fù)剪切和擠壓的作用最強(qiáng)，使流體產(chǎn)生很大的無規(guī)則流動。攪拌器中心區(qū)域附近流體的速度之所以較低，是因?yàn)檫@些位置大多不是銷在轉(zhuǎn)動過程中能夠直接掃到的區(qū)域。

3.2 ?分散性混合

本文利用銷轉(zhuǎn)動一周的分布指數(shù)來分析三維行星銷釘攪拌器的分散混合能力。糖漿的分散性混合主要由剪切和拉伸共同作用實(shí)現(xiàn)。而拉伸流動相比于剪切流動對高粘度糖漿混合的促進(jìn)作用更強(qiáng)， 所以攪拌器混合效果的好壞主要取決于拉伸流動區(qū)域面積哪個更大，停留時間更長。

由圖可知，3個點(diǎn)位置的混合指數(shù)都呈周期性變化。p1點(diǎn)的混合指數(shù)波動范圍在0.5到0.6之間，變化平穩(wěn)，屬于簡單的拉伸流動;p2和p3點(diǎn)數(shù)值波動范圍較大，混合指數(shù)峰值也較大，主要表現(xiàn)為拉伸運(yùn)動。通過對比分析可知，銷附近區(qū)域的混合效果比中心區(qū)域的混合效果更好。

由圖中可以看出，三種轉(zhuǎn)速流場內(nèi)銷附近位置的混合指數(shù)較高，數(shù)值大多在0.5到1之間，拉伸效果明顯，其他區(qū)域混合指數(shù)較低，數(shù)值處于0到0.5之間，這顯示了其他區(qū)域沒有銷附近區(qū)域混合的充分，但是整個混合器內(nèi)流體多表現(xiàn)為剪切流動。中心區(qū)不為0，表明中心區(qū)流體處于流動狀態(tài)。

圖8為轉(zhuǎn)速不同時整個流域的平均混合指數(shù)。由圖可見，當(dāng)轉(zhuǎn)速比為12∶9時，混合指數(shù)值最大且趨于平穩(wěn)，在0.41到0.43之間波動，糖漿整體流域多表現(xiàn)為剪切流動。這說明，在相同時間內(nèi)，銷轉(zhuǎn)速越大越有利于混合的進(jìn)行。

由圖9可以看出，無論哪個點(diǎn)，在低轉(zhuǎn)速的條件下，剪切速率都是相對較低且平穩(wěn)的，最大數(shù)值均不超過1.5。隨著轉(zhuǎn)速增加，剪切速率的波峰增高，最大值相應(yīng)變大，剪切作用增強(qiáng)。

這是由于轉(zhuǎn)速增加，流場會受到相當(dāng)強(qiáng)烈的推擠和擠壓作用，使其流場的速度增加，這樣流場就會獲得較大的剪切速率。同時可以發(fā)現(xiàn)轉(zhuǎn)速越大，峰值差越大，而且同一轉(zhuǎn)速下p2、p3點(diǎn)剪切速率的峰值較p1點(diǎn)大。由于剪切速率是衡量物料分散混合性能的一個重要參數(shù)，因此銷附近區(qū)域較中心區(qū)域具有更好的分散混合特性。

攪拌器中糖漿至少會經(jīng)歷一次最大剪切速率，而最大剪切速率所占的百分比，是評價物料混合效果好壞的重要指標(biāo)之一。圖10為3種轉(zhuǎn)速下最大剪切速率的對比，由圖可知，最大剪切速率依次向右偏移，也就是說，三種轉(zhuǎn)速下最大剪切速率所占的比例依次增大，剪切作用依次增強(qiáng)，效果明顯。當(dāng)最大剪切速率為10 s-1時，各模型大于10 s-1的概率分別為，模型1：2% ，模型2：48%，模型3： 82%。這說明轉(zhuǎn)速越大分散混合能力越強(qiáng)。

3.3 ?分布混合

分布混合是混合的一種重要方式，可以通過分離尺度、拉伸指數(shù)、瞬時混合效率和平均時間混合效率進(jìn)行有效評價。

圖11為在3種轉(zhuǎn)速下的分離尺度對比曲線。由于后期波動范圍較為接近，同時采用了局部放大曲線圖。在糖漿粒子剛進(jìn)入流道時，三種轉(zhuǎn)速下的模型都有一段急劇下降，這意味著在此階段分布混合能力都較強(qiáng)，隨后分離尺度曲線逐漸平穩(wěn)，各自在一個小范圍內(nèi)波動。由局部放大圖看出，轉(zhuǎn)速越大分離指數(shù)越小，因?yàn)殡S著銷轉(zhuǎn)速差值的增加，流場中糖漿的交換速度會有所不同，交換的越頻繁，越有利于粒子的分布混合，差速越大分布混合能力越強(qiáng)分離尺度越小，說明適當(dāng)?shù)脑黾愚D(zhuǎn)速有助于提高攪拌器的混合效率。

由于糖漿在流場內(nèi)流動時粒子受到復(fù)雜的拉伸和壓縮，所以任意時刻的混合指數(shù)都不完全相同，圖12為3種轉(zhuǎn)速下的瞬時平均混合效率曲線圖。由圖可知，3種轉(zhuǎn)速下的瞬時混合效率均在0.1上下做周期性波動，無論轉(zhuǎn)速大小如何，瞬時混合效率的參數(shù)總是大于0，因此可以確定流場內(nèi)的粒子受到周期性的折疊與拉伸交替的作用，這對流場內(nèi)的分布混合十分有利。同時可以看出轉(zhuǎn)速為12∶9的瞬時混合效率更優(yōu)。

4 ?結(jié) 論

通過對不同轉(zhuǎn)速下三維行星式攪拌器內(nèi)流場的數(shù)值模擬，得出以下結(jié)論。

雖然三維行星攪拌器的中心區(qū)存在一個不被銷直接掃到的區(qū)域，但是中心區(qū)的速度、混合指數(shù)、剪切速率等值始終大于零，攪拌器內(nèi)并不存在流動混合的死區(qū)。通過混合指數(shù)、剪切速率等數(shù)值的對比分析可知，銷附近區(qū)域的混合效果比中心區(qū)域的混合效果更好。適當(dāng)提高銷的轉(zhuǎn)速，會使物料混合過程加快，在相同時間內(nèi)轉(zhuǎn)速越大，物料混合能力越強(qiáng)、效果越好。

參考文獻(xiàn)

[1]宋吉昌， 李慶領(lǐng)， 楊洪順. 行星輪式攪拌器及其性能分析[J]. 化工機(jī)械， 2009， 36 （4）： 347-350.

[2]任春華， 林朝陽， 周國成. 行星攪拌機(jī)的性能特點(diǎn)及其發(fā)展[J]. 瀘天化科技， 2011 （1）： 54-62.

[3]陳明義. 新型行星式攪拌器開發(fā)及性能研究[J]. 化學(xué)工程， 2010， 38 （7）： 47-50.

[4]陳勇棠， 蒲明輝， CHENYong-tang， 等. 基于UGNX的外嚙合行星攪拌機(jī)運(yùn)動特性研究[J]. 柳州職業(yè)技術(shù)學(xué)院學(xué)報， 2017 （4）： 114-120.

[5]Jongen T. Characterization of batch mixers using numerical flow simulations[J]. Aiche Journal， 2000， 46 （11）： 2140-2150.

[6]Connelly R K， Valenti-Jordan J. Mixing analysis of a Newtonian fluid in a 3D planetary pin mixer[J]. Chemical Engineering Research & Design， 2008， 86 （12）： 1434-1440.

[7]Connelly R K， Kokini J L. Mixing simulation of a viscous Newtonian liquid in a twin sigma blade mixer[J]. Aiche Journal， 2006， 52 （10）： 3383–3393.

[8]朱霖龍， 劉雅文， 趙利軍. 長幅內(nèi)旋輪線行星攪拌機(jī)研究[J]. 工程機(jī)械， 2013， 44 （4）： 13-18.

[9]楊伶， 彭炯， 王曉瑾. 行星式攪拌釜內(nèi)三維流場的數(shù)值模擬[J]. 化工進(jìn)展， 2011 （s1）： 42-48.

[10]王超明， 鄭茂盛， 田忠， 等. 行星輪式攪拌器攪拌特性分析[J]. 化工機(jī)械， 2016， 43 （6）： 791-797.

[11]余彬彬， 楊建勇， 劉強(qiáng). 不同射流速度的油罐旋轉(zhuǎn)噴射攪拌器數(shù)值模擬[J]. 當(dāng)代化工， 2015 （6）： 1412-1414.

[12]Zhu X Z， Wang T S， Wang G. Evaluations of Flow and Mixing Efficiency in the Kneading Disks of a Novel Tri-Screw Extruder[J]. Journal of Applied Fluid Mechanics， 2016， 9 （1）： 51-60.