張培建，邢鴻雁，衛(wèi)靜怡

基于熱固耦合的雙螺桿擠壓機機筒有限元分析

張培建1,2，邢鴻雁1,2，衛(wèi)靜怡1,2

（1.天津市輕工與食品工程機械裝備集成設計與在線監(jiān)控重點實驗室，天津 300222；2.天津科技大學機械工程學院，天津 300222）

研究雙螺桿擠壓機在食品加工過程中出現(xiàn)卡頓、抱死的原因，為雙螺桿擠壓機設計及應用提供理論指導。應用ABAQUS軟件構建雙螺桿擠壓機仿真模型，通過傳感器反饋的溫度值設定溫度載荷，同時以雙螺桿擠壓機運轉時的低壓工況、正常工況、極限工況分別設定壓力載荷，運用熱固耦合理論對機筒的溫度分布、應力和變形進行仿真分析。在最大溫度工況條件下，機筒軸向熱變形量為2.356 mm，總熱變形量為2.358 mm，方向上的熱變形量為0.1324 mm，方向上的熱變形量為0.1592 mm；在最大溫度工況和極限壓力耦合作用下，機筒總變形量為2.088 mm。溫度引起的熱膨脹是機筒變形的主要原因，機筒的軸向熱變形量與總熱變形量相當，并且遠大于其他兩方向的熱變形量；與常溫環(huán)境相比，在加熱溫度工況條件下機筒的變形量隨著壓力的增大而減小。要充分考慮溫度對擠壓機性能的影響，應用過程中要合理的設置溫度參數(shù)。

雙螺桿擠壓機；機筒；熱固耦合；有限元

隨著生活水平的提高，人們對擠壓膨化食品的營養(yǎng)和口感有更高的要求；膨化食品的生產通常是由雙螺桿擠壓機完成的。由于不同食品原材料加工所需的工藝參數(shù)存在差異，在加工時的高溫高壓工作環(huán)境中，有時會使雙螺桿擠壓機出現(xiàn)卡頓、抱死等現(xiàn)象，嚴重影響食品的加工效率和質量。

雙螺桿擠壓機中機筒的作用是配合螺桿完成對食品物料的剪切、擠壓工序；機筒徑向變形過大是引起雙螺桿擠壓機卡頓、抱死的主要原因之一[1-3]。研究機筒的變形，就需要了解機筒的溫度場分布和應力分布情況。通過調研發(fā)現(xiàn)，國內外對雙螺桿擠壓機的研究主要在螺桿構型設計[4-5]和內部流場的壓力和速度分布上[6-10]。機筒的研究方面，張宇昕[11]利用有限元分析方法對雙螺桿擠出機中的雙通道機筒進行了不同壓力工況的熱固耦合分析，驗證了機筒設計的合理性；萬田瑞等[12-13]對擠出機機筒進行了設計及優(yōu)化；楊迎[14]在Fluent仿真環(huán)境中對擠壓造粒機的機筒預熱、冷卻的熱傳遞過程進行了分析，確定了機筒溫度具有非線性、大滯后性的特性。Alabadan等[15]對食品用的螺桿擠壓機冷卻系統(tǒng)進行了設計。

目前對食品加工領域所用的擠壓機機筒整體進行熱固耦合分析未見報道。目前的研究主要是針對塑料行業(yè)使用的大型擠出機，對內部含有流道的機筒進行了熱固耦合、結構優(yōu)化和傳熱方面的研究，而食品加工生產中所用的擠壓機機筒通常是沒有流道的。該研究在熱固耦合分析的理論基礎上，借助ABAQUS有限元分析軟件，以雙螺桿擠壓機運轉時的低壓工況、正常工況和極限工況為壓力工況，同時以3種食品物料加工過程中溫度傳感器反饋的數(shù)值設定溫度工況；對雙螺桿擠壓機中的機筒進行熱固耦合分析，研究機筒在不同工況作用下的應力和變形。

1 熱固耦合分析理論

由于在物料加工過程中機筒存在溫差，自由膨脹受到約束，同時機筒內部以及端面受到壓力的作用，機筒會產生熱應力和機械應力。由線性理論可知，當溫度載荷和機械載荷同時作用在機筒上時，應力由溫度變化產生的應力和壓力產生的應力疊加而成[16]，其中由溫度變化產生的應力與溫度變化成正比，應力在所有方向上相同，因此，通過溫度場理論公式，建立物理模型，對機筒的溫度場進行分析，將得到的計算結果作為載荷施加到壓力場作用時的機筒表面，計算結果即為熱固耦合的分析結果；其結果滿足如下方程。

1.1 溫度場理論公式

穩(wěn)態(tài)熱分析的能量平衡方程為[17]：

熱傳導為完全接觸的2個物體之間或1個物體不同部分之間由于溫度梯度而引起的內能交換。由于電熱片與機筒緊密接觸進行熱傳導，由傅里葉定律可得：

1.2 物理模型

機筒在溫度和壓力共同作用下的應力公式為：

式中：、、e、分別為單元的彈性系數(shù)矩陣、幾何矩陣、節(jié)點位移向量和應力矩陣。

文中采用ABAQUS軟件對機筒進行熱固耦合分析，通過熱分析模擬機筒的穩(wěn)態(tài)溫度場，將溫度場作為載荷施加到模型中，可以得到機筒的熱應力和熱變形。將溫度場作為載荷施加到單獨壓力作用時的機筒表面，得到熱固耦合的分析結果[19]，具體計算分析流程見圖1。

圖1 計算分析流程

2 機筒仿真模型

2.1 幾何模型的建立

文中以螺桿直徑為65 mm的雙螺桿擠壓機為參考，機筒外圓直徑為150 mm，螺桿與機筒內壁間隙為0.4 mm，去除不必要的倒角；給定機筒結構見圖2。其中圖2a為三維模型，2b為有限元模型。

2.2 幾何模型材料

為有效提高機筒壽命，通常選屈服強度較高和耐磨性良好的38CrMoAlA作為機筒的制作材料，38CrMoAlA的材料屬性見表1。

2.3 邊界條件

為避免物料糊化，通常將溫度沿出料口遞增來設置；文中將溫度分為4段施加到機筒外壁，同時按照加工物料的不同要求，將溫度設置分為3種工況，3種溫度工況設置情況見表2。按照螺桿的輸送段、壓縮段、定量段在機筒內對應的位置將壓力分為3段施加到機筒內壁，由于出口壓力是由傳感器測得的固定值，按照輸送段/壓縮段、壓縮段/定量段為1/1.3的比例反推出各段壓力的大小，具體壓力設置見表3。文中采用的約束類型為一端固定，另一端懸臂方式，分析時不考慮機筒自重的影響，邊界條件見圖3。

圖2 機筒結構

表1 38CrMoAlA的材料屬性

Tab.1 Material properties of 38CrMoAlA

表2 溫度工況

Tab.2 Temperature condition

表3 壓力工況

Tab.3 Pressure condition

圖3 邊界條件

3 計算結果及分析

3.1 不同壓力工況作用下機筒的應力和變形規(guī)律

根據(jù)螺桿各功能段對應機筒的位置，將機筒沿出料方向分為3段，并按照每段壓力遞增的方式施加到機筒內壁及端面，同時施加工作時的約束，得到不同壓力工況作用下機筒的應力和變形。

極限工況下的分析結果見圖4，其中圖4a為 應力云圖，圖4b為變形云圖。從圖4可以看出機筒的最大應力在出料口內壁較薄的位置，數(shù)值為333.3 MPa，最大變形量在出料口末端頂部圓柱面，數(shù)值為0.502 3 mm，由于機筒頂部開設進料口，機筒徑向變形為非對稱分布。

3種不同壓力工況下的應力和變形分析結果見表4，從表4中可看出隨著壓力的增大，機筒的最大等效應力和各方向的變形量也隨之增大。

3.2 不同溫度工況作用下機筒的熱應力和變形規(guī)律

將每種溫度工況下對應的溫度大小作為載荷，對應施加到機筒外壁的各溫度段，得到機筒的穩(wěn)態(tài)熱分布。給出了物料3溫度工況下機筒的穩(wěn)態(tài)熱截面圖見圖5。從圖5中可以看出，機筒內部的溫度分布與給定的外壁溫度載荷相當，由于不同溫度段交界處溫差較大，熱傳導比較明顯。

圖4 極限工況下分析結果

表4 3種壓力工況下機筒的最大等效應力和變形量

Tab.4 Maximum equivalent stress and deformation of cylinder under three pressure conditions

將溫度場的計算結果作為載荷施加到機筒表面，再將機筒施加工作約束條件，分析不同物料溫度工況下機筒的熱應力和熱變形。物料3溫度工況下機筒熱應力云圖見圖6，機筒的最大等效應力在法蘭的約束端，數(shù)值為426.5 MPa；每段溫差越大，交界處的熱應力越明顯。在物料3溫度工況下機筒的變形結果見圖7，其中圖7a、b為機筒沿方向的變形量，遠小于圖8c中機筒軸向變形量2.356 mm，而軸向變形量與圖7d總變形量2.358 mm大小相當；同時可以看出機筒整體變形比較均勻，變形量主要沿軸向并逐漸增加。

對應表2給出的3種溫度工況，3種物料溫度工況下機筒的熱應力和熱變形結果見表5。從表5中可以看出由于物料1的溫度工況到物料2的溫度工況跨度較大，熱應力急劇增大；而物料3的溫度工況是在物料2溫度工況基礎上每段增加5 ℃，熱應力增加了近40 MPa；3種物料溫度工況下機筒沿方向的變形量均小于機筒沿軸向的變形量，并且機筒的軸向變形量與總變形量相當。

圖5 機筒穩(wěn)態(tài)熱截面圖

圖6 機筒熱應力云圖

圖7 機筒變形云圖

Fig.8 Cylinder deformation results

表5 3種溫度工況下機筒的最大等效應力和變形量

Tab.5 Maximum equivalent stress and deformation of cylinder under three temperature conditions

3.3 不同工況下機筒熱固耦合變形規(guī)律

把溫度場的計算結果作為載荷施加到壓力工況作用時的機筒表面，對機筒進行熱固耦合分析，得到機筒在溫度和壓力共同作用下的應力和變形結果。

為更直觀地比較4種溫度工況（常溫、物料1、物料2、物料3）和4種壓力工況（無壓、低壓、正常、極限）的分析結果，不同工況分析結果見表6—7，其中表6為應力結果，表7為變形結果。從表6—7中可以看出，在常溫條件下，機筒的應力和變形都是隨著壓力的增大而增大。從表6中可以看出，加工溫度相對較高的物料（如物料2和物料3），其機筒應力反而隨著機筒內部壓力的增大而減小。從表7中可以看出，溫度是引起機筒變形的主要原因，加工溫度相對較高的物料，其機筒的總變形量隨著內部壓力的增大而減小，這是由于壓力阻礙了機筒的熱變形。

表6 不同工況下機筒應力分析結果

表7 不同工況下機筒變形分析結果

Tab.7 Analysis results of barrel deformation under different working conditions mm

4 結語

通過對雙螺桿擠壓機機筒進行有限元分析，得到了機筒分別在壓力工況、溫度工況，以及熱固耦合作用下的應力和變形結果。結果表明：機筒的最大等效應力小于材料的許用應力，滿足強度要求。溫度引起的熱膨脹是機筒變形的主要原因，機筒的軸向熱變形量與總熱變形量相當，并且遠大于其他兩方向的熱變形量；與常溫環(huán)境相比，在加熱溫度工況作用下機筒的變形量隨著壓力的增大而減小，因此，在應用雙螺桿擠壓機時，要首先考慮溫度對機筒變形的影響，合理地控制溫度避免因機筒徑向變形過大導致螺桿抱死等問題的發(fā)生。

[1] 劉曉強. 擠壓機螺桿與機筒磨損的分析[J]. 機械工程與自動化, 2007(5): 179-181.

LIU Xiao-qiang. Analysis of Wear and Tear on the Screw and Barrel of Screw Extruder[J]. Mechanical Engineering & Automation, 2007(5): 179-181.

[2] ABDULLAH D, ISMAIL T, POLYCHRONOPOULOS N D, et al. The Role of Calender Gap in Barrel and Screw Wear in Counterrotating Twin Screw Extruders[J]. Polymers, 2021, 13(7): 990.

[3] 李錦偉, 唐躍, 孫召成. 雙螺桿擠出機筒結構參數(shù)優(yōu)化設計與分析[J]. 塑料, 2015, 44(3): 56-60.

LI Jin-wei, TANG Yue, SUN Zhao-cheng. Parameters Optimization Design and Analysis of Twin Screw Extruder Barrel Structure[J]. Plastics, 2015, 44(3): 56-60.

[4] 郎珊珊. 雙螺桿構型及擠壓參數(shù)的建模與試驗研究[D]. 蘭州: 蘭州理工大學, 2014: 46-50.

LANG Shan-shan. Modelling and Experimental Study on Twin-Screw Configuration and Extrusion Parameter[D]. Lanzhou: Lanzhou University of Technology, 2014: 46-50.

[5] 韓進. 組合式雙螺桿擠壓膨化機設計及有限元分析[D]. 沈陽: 沈陽化工大學, 2018: 5-25.

HAN Jin. Numerical Simulation and Design a New Type of Combined Twin Screw Extruder[D]. Shenyang: Shenyang University of Chemical Technology, 2018: 5-25.

[6] DONG Tian-wen, JIANG Shun-liang, WU Jian-chun, et al. Simulation of Flow and Mixing for Highly Viscous Fluid in a Twin Screw Extruder with a Conveying Element Using Parallelized Smoothed Particle Hydrodynamics[J]. Chemical Engineering Science, 2020, 212: 115311.

[7] CHEN Ya-zhou, HE Qiao-lin, SHI Xiao-ding, et al. Unsteady Flow of Shear-Thinning Non-Newtonian Incompressible Fluid through a Twin-Screw Extruder[J]. Acta Mathematicae Applicatae Sinica, English Series, 2016, 32(2): 423-436.

[8] BRUNO V. Average Shear Rates in the Screw Elements of a Corotating Twin-Screw Extruder[J]. Polymers, 2021, 13(2): 304.

[9] TANG Hao, ZONG Yuan, ZHAO Ling. Numerical Simulation of Micromixing Effect on the Reactive Flow in a Co-Rotating Twin Screw Extruder[J]. Chinese Journal of Chemical Engineering, 2016, 24(9): 1135-1146.

[10] 武萌, 談慶鵬, 王君, 等. 新型交叉軸雙螺桿膨脹機轉子的設計與受力性能研究[J]. 流體機械, 2021, 49(5): 21-25.

WU Meng, TAN Qing-peng, WANG Jun, et al. Design and Mechanical Performance Research of Novel Cross-Axis Twin Screw Rotors for Expanders[J]. Fluid Machinery, 2021, 49(5): 21-25.

[11] 張宇昕. 雙螺桿擠出機機筒結構分析及機筒內聚合物流動與傳熱研究[D]. 北京: 北京化工大學, 2010: 9-19.

ZHANG Yu-xin. Structural Analysis of Twin-Screw Extruder Barrel and Numerical Simulation of Flow and Heat Transfer of Polymer[D]. Beijing: Beijing University of Chemical Technology, 2010: 9-19.

[12] 萬田瑞, 顧海, 李彬, 等. 單螺桿FDM線材擠出機機筒結構分析與優(yōu)化設計[J]. 塑料科技, 2020, 48(3): 71-74.

WAN Tian-rui, GU Hai, LI Bin, et al. Barrel Structure Analysis and Optimization Design of Single Screw FDM Wire Extruder[J]. Plastics Science and Technology, 2020, 48(3): 71-74.

[13] 楊宥人, 張煒, 殷秋娟, 等. 擠出機雙流道機筒加熱冷卻裝置: CN101288997A[P]. 2008-10-22.

YANG You-ren, ZHANG Wei, YIN Qiu-juan, et al. Extruder Double Passage Barrel Heating Cooling Device: CN101288997A[P]. 2008-10-22.

[14] 楊迎. 擠壓造粒機機筒傳熱及溫度控制研究[D]. 大連: 大連理工大學, 2018: 8-13.

YANG Ying. Research on Heat Transfer and Temperature Control of Machine Barrel of Extruding Granulator[D]. Dalian: Dalian University of Technology, 2018: 8-13.

[15] ALABADAN B, OGUNJIMI S, BABALOLA T, et al. Design of a Cooling System for a Food Extruder[J]. Journal of Emerging Trends in Engineering and Applied Sciences, 2018, 9(4): 211-217.

[16] 王飛. 航空活塞發(fā)動機典型零件的熱-固耦合強度有限元分析[D]. 南京: 南京航空航天大學, 2013: 9-12.

WANG Fei. Thermal-Mechanical Coupling Finite Element Strength Analysis of Typical Parts of Aviation Engine[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2013: 9-12.

[17] 王智博, 李志峰, 晁瑞, 等. 雙螺桿真空泵螺桿轉子熱力耦合分析[J]. 陜西理工大學學報(自然科學版), 2020, 36(4): 7-12.

WANG Zhi-bo, LI Zhi-feng, CHAO Rui, et al. Thermal Mechanical Coupling Analysis of Screw Vacuum Pump Rotor[J]. Journal of Shaanxi University of Technology (Natural Science Edition), 2020, 36(4): 7-12.

[18] 齊威. ABAQUS 6.14超級學習手冊[M]. 北京: 人民郵電出版社, 2016: 326-328.

QI Wei. ABAQUS 6.14 Super Learning Manual[M]. Beijing: People's Posts and Telecommunications Press, 2016: 326-328

[19] 龔金科, 何偉, 鐘超, 等. 基于熱固耦合的柴油機氣缸蓋有限元分析[J]. 湖南大學學報(自然科學版), 2017, 44(2): 34-39.

GONG Jin-ke, HE Wei, ZHONG Chao, et al. Finite Element Analysis of Diesel Engine Cylinder Head Based on Thermosetting Coupling[J]. Journal of Hunan University (Natural Sciences), 2017, 44(2): 34-39.

Finite Element Analysis of Twin-screw Extruder Cylinder Based on Thermo-mechanical Coupling

ZHANG Pei-jian1,2, XING Hong-yan1,2, WEI Jing-yi1,2

(1.Tianjin Key Laboratory of Integrated Design and Online Monitoring of Light Industry and Food Engineering Machinery and Equipment, Tianjin 300222, China; 2.School of Mechanical Engineering, Tianjin University of Science and Technology, Tianjin 300222, China)

The work aims to study the causes of clatter and locking in the food extrusion process of twin-screw extruder, and to provide theoretical guidance for the design and application of twin-screw extruder. ABAQUS software was used to construct the simulation model of the twin-screw extruder. The temperature load was set according to the temperature feedback from the sensor. At the same time, the pressure load was set under the low pressure condition, normal condition and limit condition of the twin-screw extruder. The temperature distribution, stress and deformation of the cylinder were simulated and analyzed by the thermo-mechanical coupling theory. Under the maximum temperature condition, the axial thermal deformation of the cylinder was 2.356 mm, the total thermal deformation was 2.358 mm, the thermal deformation in thedirection was 0.1324 mm, and the thermal deformation in thedirection was 0.1592 mm. Under the coupling effect of maximum temperature condition and ultimate pressure, the total deformation of the barrel was 2.088 mm. The thermal expansion caused by temperature was the main reason for the deformation of the cylinder, and the axial thermal deformation of the cylinder was equivalent to the total thermal deformation, and was much larger than that of the other two directions. Compared with the normal temperature environment, the deformation of the cylinder decreased with the increase of pressure under the condition of heating temperature. The effects of temperature on the performance of the extruder should be considered and reasonable temperature parameters should be set in application.

twin screw extruder; cylinder; thermo-solid coupling; finite element

TB486；TH114+.7

A

1001-3563(2022)07-0218-07

10.19554/j.cnki.1001-3563.2022.07.028

2021-06-02

寧夏回族自治區(qū)重點研發(fā)計劃（22019BFF02004）

張培建（1993—），男，天津科技大學碩士生，主攻多場耦合及優(yōu)化設計。

邢鴻雁（1969—），女，碩士，天津科技大學副教授，主要研究方向為CAD/CAE。

責任編輯：曾鈺嬋