非晶薄帶銅輥冷卻流道的數值模擬分析

程舟濟 張斐然 雷威

DOI：10.19850/j.cnki.2096-4706.2021.09.006

摘? 要：對非晶薄帶制取過程中的冷卻銅輥進行了研究，根據冷卻銅輥內的實際參數，建立換熱過程計算模型。采用流-固耦合傳熱方法對該過程進行數值模擬分析，得到冷卻水體的流場、壓力場以及溫度場分布，并探討三者之間的相互關系。通過改變銅輥水環厚度，分析冷卻水體壓力場與溫度場分布，得到了不同水環厚度對銅輥冷卻效果的影響。

關鍵詞：冷卻銅輥;流場;壓力場;溫度場;數值模擬

中圖分類號：TG139;TP391.9? ? ? ? ? ? ? ?文獻標識碼：A 文章編號：2096-4706（2021）09-0019-04

Numerical Simulation Analysis of Cooling Channel of?Amorphous Thin Strip Copper Roller

CHENG Zhouji，ZHANG Feiran，LEI Wei

（Wuhan Second Ship Design and Research Institute，Wuhan? 430064，China）

Abstract：In this paper，the cooling copper roller in the production process of the amorphous thin strip is studied. According to the actual parameters in the cooling copper roller，the calculation model of the heat exchange process is established. The fluid-structure interaction heat transfer method is used to make numerical simulation analysis of the process，and the distribution of flow field，pressure field and temperature field of the cooling water body are obtained，and the relationship among the three are discussed. By changing the thickness of copper roller water ring，analyzing the distribution of pressure field and temperature field of the cooling water body，the influence of different water ring thickness on the cooling effect of copper roller is obtained.

Keywords：cooling copper roller;flow field;pressure field;temperature field;numerical simulation

0? 引? 言

非晶態合金是由超急冷凝固，合金凝固時原子來不及有序排列結晶，得到短程有序、長程無序的組織結構，同時不具備傳統金屬材料固定的晶粒和晶界^[1]。非晶薄帶的制取采用快速凝固工藝與技術，平面流鑄（PFC）是常用的快速凝固的方法之一^[2]，具有設備簡單、操作方便、材料致密和成本低廉等特點^[3]。

冷卻銅輥是非晶制帶設備中的核心部件，其對熔潭的冷卻效果直接影響非晶晶粒結構及電磁性能的變化^[4，5]。冷卻銅輥通過內部循環的冷卻水將熱量交換并帶走，冷卻銅輥內部流道的形式和結構是直接影響薄帶質量的關鍵因素。國內外學者針對非晶薄帶制取過程中的銅輥進行了研究，LI Yong等^[6]針對銅輥的結構，建立了銅輥的冷卻過程的傳熱模型;張偉堂^[7]計算了超導熱急冷設備和銅急冷設備下非晶合金快速凝固冷卻速率，得到在銅輥表面涂鎳和不銹鋼對合金快速凝固的影響;李福森^[8]等利用Fluent軟件對雙向回流式螺旋流道的流延輥的換熱過程進行了數值模擬，得到了流延輥內部流道結構、水流量對換熱強度的影響。

本文基于銅輥的冷卻水流道，建立內流道計算模型并進行數值模擬，得到了冷卻水體的流場、壓力場以及溫度場分布。

1? 冷卻流道建模

1.1? 冷卻水流道結構設計

非晶薄帶制取時，熔融態的合金液從坩堝底部的狹縫噴嘴中流出，在噴嘴和銅輥表面狹窄的縫隙之間迅速形成熔潭，當熔體接觸到冷卻銅輥表面時發生熱量交換迅速凝固，高速旋轉的銅輥將凝固成型的薄帶從熔潭中迅速甩出。

冷卻銅輥及內流道結構主要包括密封法蘭、支撐內芯、換熱銅輥以及主軸等組成，如圖1所示，在銅輥的正上方表面接觸到熔潭，為熱量交換區域。冷卻水從主軸的進水口流入，并經主軸上小孔進入冷卻銅輥左側內流道，沿內芯徑向擴散至熱交換區域。融化的合金液噴射到冷卻銅輥上方，形成熔潭，熱量通過熱傳導方式進入銅輥，并與內芯冷卻水發生熱交換，經過熱量交換的冷卻水經過銅輥右側內流道以及主軸小孔匯集到主軸的右側出水口流出。

1.2? 冷卻水流道模型建立

為了提高計算速度、簡化計算模型，對該數值模擬過程進行適當簡化：

（1）冷卻水對換熱銅輥的冷卻過程主要集中在熔潭與銅輥表面的換熱區域，而薄帶在被拉出熔潭后溫度快速下降，因此忽略薄帶與銅輥表面的熱量交換。

（2）熔潭的尺寸相對于冷卻銅輥很小，在計算過程中簡化熔潭為長142.00 mm、寬1.15 mm、高0.40 mm的長方體塊，整體為合金液的溫度。

（3）根據格努塞爾數計算得到冷卻水與密封法蘭、支撐內芯的換熱系數，在計算模型中省略密封法蘭和支撐內芯。

（4）在其他較小的換熱區域內，忽略熱量傳遞，均看作為絕熱壁面。

（5）忽略該過程熱輻射對于換熱效果造成的影響。

基于以上簡化，得到冷卻流道的換熱計算模型如圖2所示，由熔潭、換熱銅輥以及流道冷卻水3個部分組成。

銅輥流道的設計參數如表1所示。

1.3? 邊界條件及初始條件

1.3.1? 邊界條件設置

針對建立的換熱模型，對邊界條件進行設置。

其中各個參數的初始設置如表2所示。

1.3.2? 初始條件設置

主軸高速旋轉，設定水體計算域和換熱銅輥計算域的旋轉速度與主軸的旋轉速度保持一致，即繞Z軸的圓周運動轉速為300 rev/min;設置環境壓強為標準大氣壓，初始溫度為298 K;設置整個計算域的初始溫度為298 K;設置初始湍流強度為中等強度，有助于計算結果快速收斂。

2? 結果分析

2.1? 流場分析

根據以上設置的條件對冷卻流道換熱模型進行數值模擬計算，得到換熱銅輥冷卻水體的流線圖如圖3所示。從圖中可以看出，主軸高速旋轉，帶動水體速度產生較高的旋轉分量;冷卻水由入口水環進入入水流道時，由于內壁的阻擋，流體切向速度逐漸減小為零，法向速度逐漸增大;冷卻水流出入水流道進入水環速度突增并產生分流，一部分迅速流入出水流道，另一部分碰撞到密封法蘭，速度反向，疊加流體域旋轉產生的切向速度，冷卻水在水環產生渦旋;進入出水流道的冷卻水速度先減小后增大，流入出口水環。

根據流場得到速度矢量圖如圖4所示，從圖中可以看出，速度達到最高值出現在入水流道右壁面位置、水環與入水流道交界處、出水流道右壁面位置，速度的峰值約為8.8 m/s。水道壁面隨主軸旋轉產生角速度，冷卻水受旋轉影響，不斷接觸到水道右壁面，根據壁面無滑移條件，壁面會對冷卻水提供法向力，從而使冷卻水流動加快;由于水環厚度小于流道深度，根據流體連續性方程，流體進入水環速度會增加，并且速度方向呈現與旋轉方向相反的輻射狀分布，如圖4（a）中放大圖。

圖4（b）為水環中的速度矢量分布，根據矢量箭頭的疏密反應流量的大小，長短則代表流速的大小，可以看出在相鄰的流道之間，流量變小并且根據箭頭方向得到在水環中形成了周期性的回流。

2.2? 壓力場分析

如圖5所示為數值模擬得到的流體域壓力場分布，可看出由于離心力的作用，壓力隨距主軸軸線的距離變大而增大，并在水環處達到最大。在內流道與水環的接觸位置，壓力會出現一定的回落，是因為速度在該處方向發生了變化，壓力出現一定的損失，從而導致壓力降的產生。同時在水環內環，相鄰流道中間由于出現周期性的回流出現壓力升。

根據壓力場分布可以計算得到流道出口與入口的壓降為ΔP=7.1×10⁴Pa，壓降在一定程度上反映了流體在流動時能量的損失，因此在生產中，應適當減小出水口相對入水口的壓力降。

2.3? 溫度場分析

如圖6所示得到銅輥最高溫度出現在與熔潭接觸的位置，即熱量集中交換區域，該處最高溫度為747.5 K。換熱銅輥的溫度場近似一組同心圓，分別代表不同的等溫線;由于內流道出口流體方向的改變，會在一定程度上造成壓力損失，在銅輥內壁與水環接觸位置，溫度場出現周期性的溫度峰值，根據圖6（a）在水環上出現周期性回流區域，不利于熱量的傳遞，熱量的積聚提高了溫度。根據圖6（b），銅輥截面的溫度分布，等溫線梯度出現沿換熱銅輥旋轉方向分量，通過計算得到銅輥內表面平均溫度為337.3 K。

冷卻水的溫度分布如圖7，換熱主要集中在水環外表面，等溫線呈層狀分布，構成一組近似的同心圓。冷卻水形成的回流區域，不利于溫度的耗散從而形成環狀高溫區。在該條件下，測得水環外部最高溫度為354.5 K，平均溫度值為330.3 K，出口的平均溫度值為304.5 K。同時根據水道截面圖可以得到與銅輥內側發生換熱集中在水環的外表面，而在水環內側，水溫升高幅度很小。

根據對冷卻銅輥內流道數值模擬得到，內流道的設計采用輻射直流道，方便加工、降低成本的同時，可以有效實現對熔潭的換熱，以實現非晶薄帶的連續生產。

3? 結? 論

依據設備參數簡化設計出銅輥冷卻流道結構，建立物理模型、劃分網格并進行數值模擬，得到生產參數下冷卻水的流場、壓力場以及溫度場;并分析在不同流量下水環厚度的溫度場以及壓力場變化情況，得到了冷卻銅輥傳熱過程及溫度分布規律，得出了以下結論：

（1）冷卻水由入水流道流出，速度方向發生改變并受主軸旋轉作用，會在水環處形成周期性回流。

（2）受在水環中形成的回流影響，不利于熱量的傳遞，會形成周期性溫度集中均布在水環上。

參考文獻：

[1] 屈冬冬.ZrCuNiAl非晶合金的原子尺度結構演變和變形行為研究 [D].哈爾濱：哈爾濱工業大學，2013.

[2] 巴發海，沈寧福.平面流鑄快速凝固過程的數值模擬研究進展 [J].材料科學與工程，2001（4）：97-104+96.

[3] 宋言明，楊洋.非晶態薄帶冷卻銅輥用內反饋靜壓軸承數值模擬 [J].航空學報，2014，35（4）：1157-1164.

[4] 李宇明，余紅雅，劉仲武，等.冷卻速度對快淬Fe-Si-B-Cu合金結構和磁性能的影響 [J].熱加工工藝，2010，39（6）：135-138.

[5] KANE S N，SARABHAI S，GUPTA A，et al. Effect of quenching rate on crystallization in Fe73.5Si13.5B9Cu1Nb3 alloy [J].Journal of Magnetism and Magnetic Materials，2000，215：372-374.

[6] 李勇，李家棟，劉宇佳，等.Analysis and optimization of heat loss for water-cooled furnace roller [J].Journal of Central South University，2013，20（8）：2158-2164.

[7] 張偉堂，白敏麗，楊洪武.單輥法制備非晶合金中冷卻速率的數值計算 [J].金屬功能材料，2002（1）：12-18.

[8] 李福森，王栓虎，李欣興，等.塑料薄膜與流延輥換熱過程的模擬分析 [J].機電工程，2010，27（8）：56-59.

作者簡介：程舟濟（1989—），男，漢族，湖北黃岡人，工程師，碩士，研究方向：船舶系統。

收稿日期：2021-04-06