張德偉，王興宇，汪傳生，邊慧光，李紹明

（1.青島科技大學 山東省高分子材料先進制造技術重點實驗室, 山東 青島 266061; 2.山西省大同市廣靈縣宜興中心學校, 山西 大同 037500）

基于ADINA對六棱同步轉(zhuǎn)子混煉流場的有限元分析

張德偉1，王興宇2，汪傳生1，邊慧光1，李紹明1

（1.青島科技大學 山東省高分子材料先進制造技術重點實驗室, 山東 青島 266061; 2.山西省大同市廣靈縣宜興中心學校, 山西 大同 037500）

專用于短纖維補強橡膠復合材料混煉的新型六棱同步轉(zhuǎn)子的構型特殊，故混煉過程中膠料的流動狀態(tài)也十分復雜。利用有限元分析軟件ADINA對短纖維補強橡膠復合材料混煉過程中六棱同步轉(zhuǎn)子的流場，進行了模擬分析，內(nèi)容主要包括壓力場、黏度場，剪切應力場和速度場。分析了模擬流場的流動特性，以進一步判斷和分析膠料的混煉狀態(tài)。

ADINA有限元分析軟件；混煉；流場；同步轉(zhuǎn)子；有限元分析

0 前 言

短纖維橡膠復合材料是一種新型的高分子材料，人們對它進行了深入的研究，獲得了廣泛的應用[1-8]?；鞜捠侵苽涠汤w維橡膠復合材料的第一步。短纖維橡膠復合材料的混煉有特殊要求[8]：短纖維在橡膠基體中要均勻分散；在靜態(tài)或動態(tài)條件下，短纖維與彈性體界面間應有良好的粘接性能；依力學性能、熱性能、耐化學腐蝕等不同要求，合理選擇短纖維和彈性體基體的品種；保證短纖維的長徑比（L/D）大于臨界值，只有這樣才能發(fā)揮短纖維的補強特性。因此，我們研究設計了新型六棱同步轉(zhuǎn)子，以滿足短纖維橡膠復合材料的特殊混煉要求。該轉(zhuǎn)子結構復雜，轉(zhuǎn)子棱頂與密煉室內(nèi)壁的間隙各不相同，但不同于VCMT轉(zhuǎn)子，這就使得短纖維橡膠復合材料在新型六棱同步轉(zhuǎn)子混煉區(qū)域內(nèi)的流動機理十分獨特。應用ADINA有限元分析軟件對流場進行了模擬分析，可得到壓力場、黏度場、剪切應力場、速度場等場量數(shù)據(jù)，有助于分析短纖維橡膠復合材料的混煉狀態(tài)。ADINA最早出現(xiàn)于1975年，是Automatic Dynamic Incremental Nonlinear Analysis的縮寫。ADINA具有強大的計算分析功能，而且它的準確性、高效性已獲得廣泛的認同[9-11]。

1 建立模型

1.1 物理模型和有限元模型

用ADINA有限元分析軟件對膠料流場進行分析,首先需要建立相應的流場物理模型和有限元模型。

所謂新型六棱同步轉(zhuǎn)子流場的物理模型，是在密煉室被膠料完全充滿狀態(tài)下，減去轉(zhuǎn)子的體積所形成的，或者說，是由密煉室內(nèi)壁和兩轉(zhuǎn)子外表面之間的空隙所組成的，如圖1所示。建立有限元模型則是對相應的物理模型進行有限元網(wǎng)格劃分。

圖1 物理模型

對新型六棱同步轉(zhuǎn)子流場的物理模型進行網(wǎng)格劃分所得到的有限元模型如圖2所示。所建立的有限元模型的基本單元和節(jié)點信息如圖3所示，平均劃分了66 404個節(jié)點，且流場模型共包含了328 395個四面體。

圖2 有限元模型

圖3 網(wǎng)格劃分結果

轉(zhuǎn)子的流場選用笛卡兒坐標系，坐標原點設在整個流場模型的正中心位置。為了保證計算結果的合理和準確性，采用了國際單位制，即壓力單位為Pa，長度單位為m，時間單位為s，質(zhì)量單位為kg。

1.2 數(shù)學模型

考慮到因轉(zhuǎn)子流場的幾何形狀、物料性質(zhì)、流動狀態(tài)、混煉條件等諸因素造成的流場的復雜性，為便于研究作以下假設：

1）流體為冪律流體；

2）流場為穩(wěn)定流場，即各加載物理量不隨時間變化；

3）雷諾數(shù)較小，可認為流體流動為層流；

4）慣性力、重力等遠小于黏滯力，忽略不計；

5）流體為不可壓縮流體；

6）流體在流場壁面上無滑移；

7）流體在流場內(nèi)是完全充滿的；

8）流場為等溫流場，即流場中各點溫度相同，各性能參數(shù)不隨溫度變化而變化，故不考慮能量轉(zhuǎn)換。

根據(jù)冪律本構方程，三維流體物料運動的基本方程[2]以及Carreau流體本構方程[3]即可進行流場有限元模擬分析。

冪律本構方程為：

式中：τ為剪切應力，Pa；D為剪切速度，s-1；K為冪律流體稠度系數(shù)，N·sn/m；n為冪律指數(shù)。對于牛頓流體，n=1；對于漲塑性流體，n>1，對于擬塑性流體n<1。

三維流體物料運動的基本方程為：

Carreau流體本構方程為：

式中：η0為初始黏度，也就是零剪切速率時的黏度；λ為黏彈性的特征時間；η∞為無窮剪切黏度。η為冪律指數(shù)。

2 確定邊界條件

2.1 速度邊界條件

由于是主要研究在密閉環(huán)境下，完全充滿膠料的流場模型，根據(jù)假設，流場為等溫度場，即各點溫度相同，故不考慮能量轉(zhuǎn)換。對于速度邊界條件，采用轉(zhuǎn)子旋轉(zhuǎn)、密煉室壁靜止的真實速度邊界條件。

2.2 壓力邊界條件

首先設定約束邊界條件，設置壓力為零點，也就是約束了壓力自由度，這樣，整個壓力場就有了基準，可進行對比，相應的求解模型就會變得穩(wěn)定下來。

綜上所述，對新型六棱同步轉(zhuǎn)子流場進行了有限元分析，所設置的邊界條件和物理性能參數(shù)如表1所示。

表1 轉(zhuǎn)子流場分析中的物理性能參數(shù)及邊界條件

3 模擬結果及討論

3.1 壓力場

圖4為新型六棱同步轉(zhuǎn)子流場的壓力分布圖。從圖中可以看出，轉(zhuǎn)子棱頂處的壓力非常大，這是由于膠料從轉(zhuǎn)子棱頂與密煉室壁間的狹小間隙流過時，變形劇烈，產(chǎn)生較大的壓強所造成的，而且轉(zhuǎn)子的剪切面一側(cè)的壓力值也非常大，另外一側(cè)的壓力急劇減小。這是因為位于剪切面的楔形區(qū)域的空間由小變大，膠料在小的空間里產(chǎn)生了很大的壓強。當膠料越過棱頂間隙以后，由于空間突然變大，壓強隨之急劇減小，這樣，膠料在受過強烈的剪切擠壓后，得到很好的翻滾，易于和短纖維等物料混合，其余部分的壓力分布則相對平緩。

圖4 轉(zhuǎn)子流場壓力分布圖

從圖4可以明顯地觀察到，轉(zhuǎn)子棱頂?shù)那安繅毫χ底罡?，間隙稍大一點的區(qū)域的壓力則有較大程度的減弱，而間隙最大的棱頂壓力梯度不明顯。間隙較小的棱頂區(qū)域產(chǎn)生的高壓區(qū)，有利于將短纖維擠入橡膠基體中；間隙稍大的棱頂區(qū)域可以起到分散短纖維的作用。如果膠料中的短纖維連續(xù)三次經(jīng)過較小的高剪切區(qū)域的話，就有可能使短纖維被剪斷，從而影響短纖維的補強性能；而大間隙則恰恰減少了短纖維連續(xù)三次經(jīng)過小間隙的可能性，從而避免了短纖維被剪斷的危險。

3.2 黏度場

黏度場可以反映物料在流場的不同區(qū)域內(nèi)承受剪切的程度。圖5為新型六棱同步轉(zhuǎn)子流場的黏度分布圖。

圖5 轉(zhuǎn)子流場黏度分布圖

從圖5可以看出，轉(zhuǎn)子棱頂部膠料的黏度較低，這是由于這些區(qū)域內(nèi)存在著較高的剪切應力和壓力梯度，使膠料產(chǎn)生了剪切變形行為造成的。此外，在轉(zhuǎn)子體的兩側(cè)以及轉(zhuǎn)子的中間部位存在著高黏度區(qū)域，這是因為這些區(qū)域的剪切應力較小的緣故。新型六棱同步轉(zhuǎn)子的流場位于轉(zhuǎn)子的剪切區(qū)域，剪切面前側(cè)黏度值很低，膠料混煉的后側(cè)黏度較高，在整個流場中沒有出現(xiàn)大面積的高黏度值區(qū)域，這也有利于膠料的均勻混煉。

另外，從圖5還可以發(fā)現(xiàn)，轉(zhuǎn)子頂隙部分的黏度值最低，因為此部位為高剪切區(qū)，短纖維可以很好地混入橡膠基體中。

3.3 剪切應力場

圖6所示為六棱剪切型轉(zhuǎn)子流場的剪切應力分布圖。從圖6可以得知，轉(zhuǎn)子棱頂處的剪切應力較大，短纖維和其它填料的分散也主要發(fā)生在這一區(qū)域。另外，轉(zhuǎn)子圓盤處的剪切應力也較大，其余部分的剪切應力較小，沒有明顯的剪切應力存在。這是因為這些區(qū)域的間隙較大，因而剪切應力較小，膠料在這些區(qū)域不容易受到剪切作用。

圖6 轉(zhuǎn)子流場的剪切應力分布圖

加之，黏度是剪切速率的函數(shù)，剪切速率越高，黏度愈低，所以也可以從剪切應力的分布云圖上分析出轉(zhuǎn)子流場的黏度高低。在轉(zhuǎn)子的棱頂區(qū)域，存在著較高的剪切應力，它可使膠料產(chǎn)生剪切變形行為，因而該區(qū)域的黏度較低，而在轉(zhuǎn)子體的兩側(cè)以及轉(zhuǎn)子的中間部位剪切應力較小，因而這些區(qū)域的黏度就較高。剪切應力大的區(qū)域，膠料黏度就低；反之，剪切應力小的區(qū)域，膠料黏度就高。剪切應力分布圖和圖5的黏度場分布趨勢是完全相反的。

3.4 速度場

在混煉過程中，密煉室內(nèi)膠料的流動行為對膠料所產(chǎn)生的剪切和拉伸作用，是影響分布混合和分散混合效果的重要因素，尤其對短纖維的分散而言。因此，研究膠料在轉(zhuǎn)子流場內(nèi)的流動方式和流動特點是非常必要的。

圖7是新型六棱同步轉(zhuǎn)子流場的速度分布云圖。從圖7可以看出，由于轉(zhuǎn)子表面各點的回轉(zhuǎn)半徑不同，其速度值也各不相同，靠近突棱頂部的流動速度較高，由轉(zhuǎn)子體本身形成的膠料流動速度較低。

速度本身是矢量，具有高低和方向性，圖7只表示了速度的高低，對于膠料的真正流向還缺乏直觀的表示，乃需要進一步分析一系列速度的矢量分布，以研究膠料流動的方向性。圖8 所示是新型六棱同步轉(zhuǎn)子流場上的速度矢量分布圖。從圖8可以看出，轉(zhuǎn)子表面的速度隨著各點半徑的增大而增高，半徑越大，間隙越窄，徑向速度梯度就越大。因此，轉(zhuǎn)子突棱的推力面采用較小的楔入角，便可以增大楔形區(qū)域內(nèi)的徑向速度梯度，從而更好地對膠料實施剪切作用。

圖7 轉(zhuǎn)子流場速度分布云圖

圖8 轉(zhuǎn)子流場上的速度矢量分布圖

分析了新型六棱同步轉(zhuǎn)子流場的速度矢量分布圖，并研究了膠料周向運動的狀態(tài)后，為了再度分析膠料的軸向運動情況，給出了流場軸向的速度分布情況。如圖9所示，軸向速度場分布是指沿Z軸方向的速度分量。從圖9中可以直觀地觀察到軸向速度的分布。新型六棱同步轉(zhuǎn)子流場基本上呈左右軸對稱分布結構，轉(zhuǎn)子流場就好像分為兩個小循環(huán)，左側(cè)一進一出，右側(cè)一進一出。

圖9 轉(zhuǎn)子流場截面上的軸向速度分布圖

4 結 語

ADINA有限元分析軟件能夠?qū)π滦土馔睫D(zhuǎn)子的流場進行較為準確的模擬計算。通過模擬計算的結果，可以對膠料的混煉狀態(tài)進行分析。新型六棱同步轉(zhuǎn)子使短纖維—橡膠復合材料的混煉達到了混煉要求，使混煉膠的力學性能有較大幅度的提高。因此，新型六棱同步轉(zhuǎn)子適用于短纖維—橡膠復合材料的混煉工藝。另外，在有限元分析領域ADINA的應用也越來越廣泛。

[1] Wang C S, Liu C J, Bian H G. Radial Orientation Mechanism and Experimental Research of Short Fiber in Tread Compound[J]. Journal of Donghua University (Eng. Ed.), 2009, 26(6):666-672.

[2] 袁小紅.短纖維補強發(fā)泡橡膠復合材料的性能研究[J].橡膠工業(yè)，2011，58(7)：410-414.

[3] 汪傳生，張德偉，邊慧光.混煉參數(shù)對短纖維-橡膠復合材料性能的影響[J].功能材料，2011，42（8）：1448-1452.

[4] Pramila D S. Nair B, Jabin T, et al. Mechanical, thermal, and microwave properties of conducting composites of polypyrrole/polypyrrole-coated short nylon fibers with acrylonitrile butadiene rubber[J]. Journal of Applied Polymer Science, 2012, 126(6): 1965-1976.

[5] Shirazi M, Talma A G, Noordermee J M W. Viscoelastic properties of short aramid fibers-reinforced rubbers[J]. Journal of Applied Polymer Science, 2012, 128(4): 2255-2261.

[6] Supri A G, Tan S J, Ismail H, et al. Enhancing interfacial adhesion performance by using poly(vinyl alcohol) in (low density polyethylene)/(natural rubber)/(water hyacinth fiber) composites[J]. Journal of Vinyl and Additive Technology, 2013, 19(1): 47-57.

[7] Wang C S, Zhang D W, Bian H G, et al. Experimental Study of Screw Speed Effects on Quality of Short Fiber-Rubber Composite Material[J]. Advanced Materials Research, 2011, 221: 369-372.

[8] 朱玉俊.彈性體的力學改性-填充補強及共混[M].北京:北京科學技術出版社, 1991.

[9] 邊慧光，楊文超，汪傳生. 密煉機的三維流場模擬[J]. 橡膠工業(yè)，2007,52(12): 750-753.

[10] 張軍,鄭捷慶.擬塑性冪律流體的摻氣管輸減阻的數(shù)值模擬[J]. 管道技術與設備，2006(6):1-3.

[11] 邊慧光.密煉機混煉過程中內(nèi)部流場和溫度場的模擬研究[D]. 青島：青島科技大學，2006.

[責任編輯：張啟躍]

FEA of Flow Fluid for Six-Wing Synchronous Based on ADINA

Zhang Dewei1, Wang Xingyu2, Wang Chuansheng1, Bian Huiguang1, Li Shaoming1
(1.Shandong Provincial Key Laboratory of Polymer Material Advanced Manufacturings Technology, College of Electromechanical Engineering, Qingdao University of Science & Technology, Qingdao 266061, China;2.Guangling Central School of Datong City, Datong 037500, China)

Short fi ber-rubber composite material has special mixing requirements. So, the new type six-wing synchronous rotors which have been designed for the mixing of short fi ber-rubber composite material have special structure, and considering the characteristic of the short fi ber-rubber composite material, the flow situation of the mixed rubber is very complex. However, the flow field of the sixwing synchronous varying clearance rotors can be analyzed by ADINA. The fl ow situation and mixing situation of the mixed rubber can be analyzed and judged by the FEA results such as the pressure fi eld, viscosity fi eld, shearing force fi eld and velocity fi eld.

Automatic Dynamic Incremental Nonliner Analysis(ADINA); Mixing; Flow Field; Synchronous Rotor; Finite Element Analysis(FEA)

TQ 330.6+3

B

1671-8232(2015)04-0050-05

2014-06-21

國家自然科學基金項目（51175274），山東省科技發(fā)展計劃項目（2013TD16006），山東省自然科學基金項目（ZR2013EMM002）。

張德偉（1983— ），男，山東煙臺人，青島科技大學在站博士后，主要從事高分子材料加工機械成套設備及機電一體化研究。