薄板坯连铸连轧过程轧辊温度场及热凸度的研究_郭忠峰

国家自然科学基金仪器专项项目(50527402) 郭忠峰，博士生，zf_guo@

薄板坯连铸连轧过程轧辊温度场及热凸度的研究

郭忠峰1， 李长生1， 徐建忠1， 刘相华1， 王国栋1， 刘太斗2

（1. 东北大学轧制技术及连轧自动化国家重点实验室，东北大学； 2. 本钢热连轧厂，本溪）

摘 要：薄板坯连铸连轧生产中，轧辊的热变形是影响带钢板形的重要因素之一，因此研究薄板坯连铸连轧过程中的轧辊热变形对提高板形控制系统精度具有重要的意义。本文提出了一种简化分析轧辊温度及热凸度的方法，依据生产过程中轧辊实际的边界条件建立了轧辊温度场及热凸度模型，考虑了轧制过程瞬态热接触及复杂的边界条件，并借助有限元分析软件ANSYS 模拟了轧辊的瞬态温度场及热凸度，得到了轧辊表面各点的温度变化及轧辊热凸度变化规律。轧辊热凸度有限元模拟结果与实验测量结果吻合较好，表明了采用简化有限元法得到的结果是正确的。

关键词：温度场；热凸度；有限元；轧辊

Analysis of Roll Temperature Field and Thermal Crown

in TSCR

GUO Zhong-feng 1, LI Chang-sheng 1, XU Jian-zhong 1, LIU Xiang-hua 1,

WANG Guo-dong 1, LIU Tai-dou 2

(1. The State Key Laboratory of Rolling and Automation, Northeastern University, China;

2. BenXi Hot Strip Mill, BenXi, China)

Abstract: Roll thermal crown variation in thin slab continuous casting and rolling (TSCR) is an important factor which affects strip profile. Prediction of roll thermal crown variation in rolling is a difficulty and weakness for many kinds of heat transfer ways and many physical parameters could not be confirmed well and truly. Setting up roll temperature field model which consider different affection factors is the key to solve the difficult problem. A new simplification method of analyzing roll temperature field is presented. Transient roll temperature field and thermal crown in TSCR are simulated by ANSYS software. Transient thermal contact and complex boundary condition are considered. The prediction results of thermal crown are good agreement with measurement data. Key words: TSCR; roll; temperature field; thermal crown; FEM

1 引言

薄板坯连铸连轧现已在世界范围内取得了很大的发展[1,2]。

准确预报热轧工作辊在轧制过程中的热变形，对提高带材板形质量具有重要的实际意义。热轧生产中工作辊复杂的热交换条件，使得轧辊热凸度在线测量困难。目前对温度场和热凸度变化规律的研究取得了许多重要进展，如国外学者Stevens [3]采用实验方法获得了轧辊温度场的分布，

国内学者[4~7]等采用有限元法分析了轧辊温度场，李俊洪[8]等采用差分法进行了分析。对轧辊温度场进行三维有限元分析可获得任意时刻任意节点的温度分布，但需要耗费大量的时间。而在轧辊温度场的二维分析过程中做了大量简化，结果不甚理想。

采用新的简化方法，借助ANSYS 软件，全面考虑各影响因素对热带轧制过程中轧辊温度场进行了数值模拟，模拟结果与实测结果吻合较好。

487

2 边界条件

边界条件如图1及图2所示。在图1中，轧辊表面依据热交换形式被分成了8个不同的区域。

A ～

B 轧辊和带钢的强烈对流区；B ～

C 带钢入口辐射区；C ～

D 入口水冷区；D ～

E 入口自然对流区；E ～

F 工作辊和支撑辊接触区；F ～

G 出口自然对流区；G ～

H 出口水冷区；H ～A 出口带钢辐射区。

如图2所示，点(1-8)为边界条件特征点。工作辊表面点2-3-7-6围成的区域为带钢与轧辊接触区，带钢热量通过热传导与辐射的形式进入轧辊，又被冷却水与支撑辊带走。由点3-4-8-7及点1-2-6-5组成的区域没有与带钢接触，但与冷却水及空气接触。

1—工作辊；2—支撑辊；3—喷嘴；4—集水管；5—挡水板

图1 边界条件横断面示意图

Fig.1 Boundary conditions in cross section

strip

cooling nozzle

图2 轧辊轴向热边界条件示意图

Fig.2 Sketch of WR thermal boundary conditions in

radial

在分析过程中，可对某些热交换形式进行简化，将其忽略。如轧辊与带钢间的摩擦热只能通过实验的方法获得。工作辊和支撑辊之间的热量传递就可忽略，因支撑辊从工作辊带走一部分热量，同时又将摩擦热传递给工作辊，这两部分热量可近似认为抵消。由工作辊轴承产生的热量因离轧辊中部较远，对轧辊热凸度的影响甚微，故也可忽略不计。

3 基本假设

轧辊温度场是一个3维非稳态系统。在轧制过程中，轧辊轴向、径向和周向的温度都要发生变化，考虑到轧辊的回转周期与热凸度对轧制条件变化的响应时间相比为二阶小，所以忽略轧辊在圆周方向的温度变化，这样就可简化为二维非稳态系统。在分析过程中由于轧辊几何形状的对称性，为节约计算时间，选取工作辊的四分之一作为分析对象。

图3 基本假设简图

Fig.3 Schematic plane of basic hypotheses

在计算过程中，事先做了两个基本假设： ①平面假设

因轧辊几何形状的对称性及边界条件的周期

性，工作辊可被分成相当多的平面，并且每个平面都具有相同的热行为，只是时间不同而已。举个例子来说，第l 个平面的热行为等于第k 个平面旋转过β角度（图3）。

②平均初始条件假设

为节约分析时间，现假设一个平面在相同的边界条件区域内时具有相同的初始条件。

根据上述基本假设，整个轧辊的热行为可看作是一个平面在不同时刻的热行为的累积，于是复杂的3维轧辊温度场问题就简化为了简单的2维平面

问题，这必将节约大量分析时间。

4 有限元数学模型

采用PLANE13单元非均匀划分有限元模型（图4），轧辊表面及与带钢边部接触区域网格划分较密，因为在这些区域温度变化较剧烈。整个模型中共有176个节点和150个单元。

根据热平衡原理，2维瞬态温度场可表示为：

22

2

21()T

T T T

v t

r r r z c k ρ∂∂∂∂∂∂∂∂=++ (1)

热流密度可表示为温度梯度的函数:

(,,)T n q z r t k ∂∂=− (2)