木塑擠出成型特性分析與影響因素數值仿真研究

劉亞秋　李明　孫垚

摘要 [目的]為木塑型材生產提供合理的工藝控制參數，提升木塑擠出型材的質量。[方法]針對木塑擠出成型工藝控制問題，從木塑熔融體到擠出成型的工藝特性角度進行了分析研究，重點得出了模具域熔融體流場（溫度分布、壓力分布）與剪切率（流速）的關系特性及工藝閥值特征。通過建立單螺桿擠出機計量段和法蘭處三維計算模型，對木塑擠出成型的若干典型工況流場進行了數值仿真。[結果]在一定的工藝閥值條件下，加熱段的溫度擾動對熔融體流場影響不大，而螺桿轉速對流場影響較大，隨著螺桿轉速的增大，計量段壓力逐漸減小，剪切率線性增大，粘度減小。[結論]驗證了該文分析結論的正確性和一致性，根據分析結果得到較理想的木塑擠出工藝控制參數。

關鍵詞 木塑擠出機；熔融體流場分析；數值仿真

中圖分類號 S784 文獻標識碼 A 文章編號 0517-6611（2014）03-00823-05

Abstract [Objective]The paper was in order to achieve better controlling parameters of extruder，so as to improve the guality of wood plastic extrusion profiles.[Method]Focusing on process control of WPC extrusion，this paper studyed from the angle of WPC melt and process characteristics of extrusion molding，to get distribution curve and process threshold value of mould segment in flow field （temperature field，pressure field） of melt and screw speed （velocity）.numerical simulation on flow field of some typical working conditions of WPC extrusion.

by constructing three-dimensioal computational modelling of the metering section and flange of single screw extruder.[Result]Experimental results show that heating temperature little influences flow field while revolving speed greatly influences flow field with screw speed increases，pressure of metering section decreases，shear rate increases and viscosity decreases，under the condition of process threshold.[Conclusion]Result verifies the correctness and consistency of the conclusion，according to results of analysis，better controlling parameters of extruder are achieved.

Key words Extruder of wood and plastic；Flow field analysis of melt；Numerical simulation

擠出成型是指物料通過擠出機料筒和螺桿間的作用，邊受熱塑化，邊被螺桿向前推送，連續經機頭和模具而制成各種截面制品或半制品的一種加工方法[1]。擠出成型過程是由多個具有相關聯系的因素共同作用而成，為解決擠出機生產實驗的穩定控制，滿足市場對高品質木塑型材擠出控制系統的要求，提高木塑擠出型材擠出機實驗效率已成為生產及研究最為關鍵的任務。在擠出機中計量段是聚合物熔體成型的關鍵部分，受溫度和螺桿轉速影響比較大，K.Wliczynsk通過構建數學模型，對聚合物質量流體速率、溫度、壓力、能量轉化進行了分析，并進行了試驗對比[2]。B.BroelPlater通過分析擠出機工作機理，對溫度控制進行了改進[3]。E.Lai和DW.Yu構建的新數學模型實現對單螺桿擠出機聚合物流性進行全面分析[4]，但該模型都沒有將實驗參數及相關性進行充分考慮。

總之，擠出機研究在經歷數十年的發展后已有了基本的理論基礎，各種基礎模型的計算軟件也在某種程度上反映了生產實際情況，當前擺在生產廠家和科研院所的中心課題是如何將試驗參數與理論模型研究緊密結合。該文通過研究單螺桿擠出機木塑型材流體影響因素的分布特點，進一步明確木塑熔體主要工藝控制參數對擠出型材影響特性關系，為木塑型材生產提供合理的擠出機工藝控制參數，減少擠出機實驗次數和提高實驗效率，從而提升木塑擠出型材的質量。

1 木塑成型方法與實驗設計

實驗原料及配比：聚乙烯（PE），顆粒狀，30%；木粉，粉狀，65%；偶聯劑，顆粒狀，4%；潤滑劑，粉狀，1%。

采用SJ-45型號單螺桿擠出機對木塑復合型材擠出進行了實驗研究，為數值模擬結果提供數據對比。實驗系統如圖1所示，通過檢測點1），2），3）溫度傳感器測得溫度，4）壓力傳感器檢測出法蘭處流體壓力，通過控制系統設置螺桿轉速。

設定相同材料，溫度一定時，分別改變螺桿轉速為10，20，30 r/min；同理，螺桿轉速一定時，分別改變溫度為150，160，170 ℃，分別進行實驗觀察得到壓力變化情況。實驗測量了圖2所示A處（法蘭處流道中心）的壓力PA，其結果如表1所示。從表1中可以看出，計量段機筒內壁加熱溫度Tw對PA影響較小，螺桿轉速N對PA影響較大，且PA隨著轉速N的增大而增大。

2 計算模型與測試

由于實驗難以獲得詳細的流場信息，同時為了與實驗結果進行比較，采用計算流體動力學（CFD）軟件ANSYS FLUENT 14.0為研究工具，建立了實驗中單螺桿擠出機計量段和法蘭處的三維計算模型，對表1所示6個不同工況下的流場進行了數值模擬，得到了不同螺桿轉速和計量段加熱溫度下的壓力分布和溫度分布等流場信息，分析了螺桿轉速和計量段加熱溫度對流場的影響，以期獲得加熱溫度、轉速和壓力之間的關系，通過與實驗結果的對比，驗證采用CFD仿真模擬實際生產過程和指導實際生產的可行性。

3 計算結果與分析

3.1 實驗結果的比較及分析

通過數值模擬得到了不同工況下流道A處與流道進口的壓力差ΔP，其結果如表1所示。從表1中可以看出，不同加熱溫度Tw下，相同轉速N下的ΔP相同，這可能是由于相同轉速下木塑熔體質量流量相同和沒有考慮溫度對粘度的影響（如Ⅳ式所示）而引起的。而轉速對ΔP有較大影響，且ΔP隨N的增大而增大，這與實驗結果PA在定性上相符。

3.2 熔體流場分析

由于不同工況下流場信息相似，因此以160-20工況為例進行分析。圖4～8分別給出了x=0平面，圓環面以及如圖2所示LINE-1和LINE-2這2條直線上的流場信息。其中，LINE-1穿過螺紋，位于螺槽中，其空間位置為（x=0 m，y=0.021 5 m）；LINE-2則位于螺紋和機筒內壁之間的縫隙中，其空間信息為（x=0 m，y=0.022 75 m）；圓環面為LINE-2所處的圓環面。從圖4所示的速度分布可以看出，x方向上分速度u呈對稱分布，左右兩側分速度大小互為相反數，這是由于螺桿旋轉而造成的；z方向上分速度w（軸向速度）的最大值出現在螺槽中靠近螺紋處，螺槽中軸向速度沿軸向呈先減小后增大的趨勢，w的最小值分布于螺紋和機筒內壁的縫隙中，說明木塑熔體從縫隙中通過的流量很小，絕大部分都在螺桿的驅動下沿螺旋線流動。圖5給出了流場中壓力的分布，需要說明的是，后文圖中給出的壓力值均為相對于流道出口處壓力的相對壓力，從圖中可以看出，在計量段壓力呈螺旋分布，沿螺旋線逐漸增大，如圖5（b）所示，法蘭處壓力基本保持不變；圖5（c）從定量上更直觀地反映了這一分布趨勢，螺槽中壓力沿軸向降低，相鄰螺槽相同位置處壓力逐漸增大，法蘭處壓力保持不變，此結果與彭炯等[7-9]的數值模擬結果在定性上相符，說明了該研究計算結果的準確性；此外，同一軸向位置上不同徑向位置處的壓力基本相同。值得注意的是，圖5（c）所示的LINE-1的壓力分布存在間斷，這是由于LINE-1穿過螺紋，間斷為螺紋所處位置，后文中亦是如此。從圖6的溫度分布可以看出，木塑熔體的溫度在離進口不遠的地方就達到加熱溫度，說明木塑熔體在進入計量段不久后便充分混合。從圖7可以看出，在計量段剪切率呈周期性分布，剪切率最大值存在于螺紋和機筒內壁之間的縫隙處，這與王建等[8，10]的結果相符，且離螺槽底部越近，剪切率越大。圖8給出了木塑熔體粘度的分布，從圖中可以看出，粘度與剪切率呈現相反的分布趨勢，在計量段同樣具有周期性的分布，但在螺紋和機筒內壁之間的縫隙處為最小值，這與王建等[8]的結果相符，且離螺槽底部越近，粘度越小。剪切率和粘度分布趨勢相反可由（Ⅳ）式來解釋，描述木塑熔體粘度的Cross粘度模型說明其具有剪切稀變的特性，即剪切率越大，粘度越小；而在法蘭處木塑熔體粘度較大，這是由于此處木塑熔體不再受螺桿的直接驅動，因而速度較?。▓D4），剪切率也較小[圖7（a）]，且熔體由法蘭流向模具為突縮流動，熔體積聚在近出口流道角落處，此處熔體速度小[圖4（c）]，剪切率小，因而粘度最大。

（1）加熱溫度對流場影響不大，轉速對流場影響較大，

數值模擬結果與實驗結果在定性上相符，流場信息與已有數值模擬結果相符，說明該研究計算結果的準確性，也表明了可以通過CFD數值模擬技術指導實驗和實際生產，一方面數值模擬可以獲得更為詳盡的流場信息，另一方面能有效減少實驗帶來的人力物力，降低實驗和生產成本。

（2）隨著螺桿轉速的增大，計量段壓力逐漸減小，剪切率線性增大，粘度減小。此外，若進一步研究加熱溫度對流場的影響，可在Cross模型中加入溫度對木塑熔體粘度的影響，以獲得更為精確的結果。

參考文獻

[1]

LIU S J，SU P C.Novel three-dimensional in-cavity transient temperature measurements in injection molding and fluid-assisted injection molding [J].Polymer Testing，2009，28：66-74.

[2] WILCZYNSKI K.A mathematical model of the single-crew extrusion process.Part VIII.Extrusion dies[J].Polimery，1998，43（1）： 44-48.

[3] BROEL PLATER B.Optimization of the control of polymer extrusion processes[J].Polimery，1997，42（6）： 386-397.

[4] LAI E，YU D W.Modeling of the plasticating process in a single-screw extruder：A fast-track approach[J].Polymer Engineering and Science，2000，40（5）：1074-1084.

[5] LIU Y Q，CHANG L，ZHU L K.Modelling of airflow circulation systems in wood drying processes[J].Int.J.Modelling，Identification and Control，2011，12（4）：378-385.

[6] LIU Y Q，CHANG L，ZHU L K et al.Researching of airflow circulation system in wood drying kiln based on array fans[J].Int.J.Modelling，Identification and Control，2011，13：126-133.

[7] 彭炯，陳晉南.同向旋轉雙螺桿擠出機計量段中聚合物擠出的模擬[J].中國塑料，2001，15（7）： 39-42.

[8] 王建，郭迪，陳晉南.數值模擬研究螺筒結構對單螺桿擠出機性能的影響[J].塑料科技，2012，40（10）：74-78.

[9] 王小龍，徐百平，劉躍軍，等.嵌入式行星螺桿擠出機內三維流動數值模擬[J].高分子材料科學與工程，2012，28（7）： 171-174.

[10] 朱向哲，袁惠群.三螺桿擠出機擠出特性的數值模擬[J].高分子材料科學與工程，2008，24（3）： 32-35.