第4章-热连轧机轧辊温度场及热凸度研究.doc

精品文档，仅供学习与交流，如有侵权请联系网站删除第4章 热连轧机轧辊温度场及热凸度研究在热轧带钢生产中，实时变化的轧辊热凸度是影响带钢板形的重要因素在带钢生产中，轧辊热交换十分复杂，包括带钢向轧辊传递热量，带钢与轧辊相对运动产生的摩擦热，轧辊与空气、集管冷却水以及与轧辊轴承的热交换等因此，研究和开发高精度的轧辊温度场及热凸度模型具有十分重要的意义[148]4.1 传热学基本定律传热的基本方式有三种：热传导、对流和辐射在计算轧辊温度场及热凸度时需同时考虑上述三种传热方式[149,150]1) 热传导的富立叶简化导热定律热传导即物质内部或物质之间的热传递，在这里为轧辊层或段间节点之间的热交换富立叶简化导热定律为： (4.1)式中，为物质间时间内传递的热量，J；为物质的导热系数，W/(mm·℃)；为垂直于热流的横截面积，mm2；为热流方向上的路程，mm；、为两端介质的温度，℃；为传热时间，s通常以热流密度来表示富立叶传导定律即： (4.2)上式是由典型的单一介质两端传热得到的，但其仍具有普遍意义，只不过一项将要有所改变如图4.1所示图4.1 不同介质间的热传导Fig. 4.1 Heat transforming between different mediator对于和两种不同的介质，厚度分别为和，导热系数分别为和，两种介质间的传热量为： (4.3)由式(4.1)可得： (4.4)(2) 对流传热的牛顿定律对流传热是固体表面与其相邻的运动流体之间的换热方式。

在这里为轧辊与其周围气体及冷却水之间的热量交换对流传热用牛顿定律描述为： (4.5)式中，为交换的热量，；为表面换热系数，W/(mm2·℃)；为流体与固体之间界面面积，； 为热量交换时间，；为流体与固体温差，℃根据牛顿冷却定律，可得冷却水与工作辊之间的热交换公式： (4.6)式中，为i机架工作辊j单元由t到t+Δt时间内冷却水吸收的热量，J；为冷却水与工作辊表面换热系数，W/(mm2·℃)；为工作辊与冷却水接触部分占整个轧辊圆周的比例，一般取0.6；为t时刻i机架冷却水的温度，℃；为t时刻i机架工作辊的j单元温度，℃；为单元长度，mm；为轧辊与冷却水热量交换时间，s同理，可得空气与工作辊之间的热交换公式： (4.7)式中，为i机架工作辊j单元由t到t+Δt时间内空气吸收的热量，J；为空气与工作辊表面换热系数，W/(mm2·℃)，实际应用时取=10 W/(m2·℃)；为工作辊与空气接触部分占整个轧辊圆周的比例，一般取0.4；为t时刻i机架空气的温度，℃；为t时刻i机架工作辊的j单元温度，℃；为单元长度，mm；为轧辊与空气热量交换时间，s3) 能量守恒定律现设某个体系的质量不变，那么可借助能量守恒定律来描述该体系的能量变化及其与周围介质的联系。

此时，能量守恒定律可表示为： (4.8)式中，为进入体系的所有形式的热量；为体系本身产生的热量，即内热源产生的热量；为流出体系的所有形式的热量；为体系内储能量的变化对于板带轧制过程，进入辊系的热量主要有金属向轧辊的热传导，假设轧辊不存在内热源项，即；项主要体现在轧辊温度的变化；项主要体现为轧辊与周围介质，例如与冷却水或空气的对流传热，故能量守恒定律可表示为： (4.9)4.2 轧辊温度场及热凸度计算模型有限差分格式的建立主要有两种方法：(1)从能量守恒观点出发建立差分格式的解法，这种方法被盐崎所采用；(2)从热传导方程出发建立差分格式的解法，这种方法被有村所采用第二种方法虽然数学概念清晰，但存在边界节点温度方程的截断误差与内节点不一致的问题，而且当采用非均匀网格时所得到的节点温度方程较复杂而第一种方法，物理概念清晰，较易解决上述问题，特别是在处理热交换边界条件时存在极大的灵活性4.2.1 工作辊模型单元划分轧辊温度场是一个三维非稳态系统随着轧制过程的进行，轧辊轴向、径向和周向的温度都要发生变化，考虑到轧辊的回转周期与热凸度对轧制条件变化的响应时间相比为二阶小，可忽略轧辊在圆周方向的温度变化，这样就将复杂的三维温度场问题简化为轴对称问题。

同时为简化计算还忽略了轧辊与带钢之间摩擦热和带钢的变形热图4.2为轧辊四分之一有限差分模型图4.2 工作辊单元划分Fig. 4.2 Division of work roll elements认为轧辊轴承处绝热，即轧辊与轴承不发生热传导在半个轧辊辊身划分10段，在径向上划分4层因为轧辊外层单元和高温带钢接触，是轧辊温度和热凸度变化敏感区，故采用非均匀单元划分法，轧辊由表及里，各层厚度逐渐增加4.2.2 轧制过程传热数学模型4.2.2.1 无钢时的传热计算模型(a) 轧辊水冷时传热计算模型(1) 无钢水冷时轧辊的散热量 (4.10)式中，为水冷时轧辊与水的散热传输参数，W/K；为水温，℃；为轧辊外层的温度，℃；为轴向段数2) 无钢水冷时轧辊的吸热量 (4.11)(b) 轧辊空冷时传热计算模型(1) 无钢空冷时轧辊的散热量 (4.12)式中，htres为待轧空冷时轧辊与空气的散热传输参数，W/K；为空气的温度，℃；为轧辊外层的温度，℃；为轴向段数2) 无钢空冷时轧辊的吸热量同(4.11)式4.2.2.2 轧制过程轧辊水冷时传热计算模型带钢轧制过程中，带钢与轧辊接触形式如图4.3所示。

图4.3 轧辊与带钢接触形式Fig. 4.3 Contact fashion of roll and strip轧辊各段与带钢接触状态可分为完全接触段、部分接触段和非接触段3种形式由于轧辊各段与带钢接触状态不同，轧辊各段传热计算模型也不同a) 轧辊与带钢完全接触段传热计算模型(1) 轧辊的散热量 (4.13)(2) 轧辊的吸热量 (4.14)式中，为带钢设定温度，℃；为轴向段数；为冷却水的温度，℃；为冷却水温度的修正值，℃；为轧辊外层的温度，℃；为轧辊与带钢的热传输参数，W/K；为有钢时轧辊表层各段与冷却水之间热传输参数，W/Kb) 轧辊与带钢部分接触段传热计算模型(1) 轧辊的散热量 (4.15)(2) 轧辊的吸热量 (4.16)式中，为带钢设定温度，℃；为轴向段数；为轧辊轴向段的长度，mm；为冷却水温度，℃；为冷却水温度修正值，℃；为轧辊外层温度，℃；为轧辊与带钢热传输系数，W/K；为有钢时轧辊表层各段与冷却水之间热传输参数，W/K；为无钢时轧辊各段与水的热传输系数，W/K；为轧辊轴向段与带钢接触轧辊长度，mmc) 轧辊与带钢非接触段传热计算模型(1) 空冷时非接触段传热计算模型 (4.17)(2) 水冷时非接触段传热计算模型 (4.18)式中，为轧辊各段与空气的热传输参数，W/K；为轧辊各段之间及轧辊与水的热传输参数，W/K；为冷却水的温度，℃；为空气的温度，℃；为轧辊外层k段的温度，℃；为轴向段数。

4.2.2.3 轧辊径向单元之间的热传输模型 (4.19)式中，为单位时间内轧辊散发的热量，W；为单位时间内轧辊吸热量，W；为轧辊第i层第k段交换的热量，W；为轧辊径向单元间的热传输系数，W/K；为轧辊外层的温度，℃；为轧辊第i+1层第k段的温度，℃；i为径向层数；为轴向段数4.2.2.4 轧辊轴向单元之间的热传输模型 (4.20)式中，为轧辊第i层第k段轴向散热量，W；为轧辊轴向单元间的热传输系数，W/K；为轧辊外层的温度，℃；为轧辊第i+1层第k段的温度，℃；i为径向层数；为轴向段数4.2.3 轧辊轴对称温度场计算模型轧辊温度场是求解轧辊热凸度的前提，即先从能量守恒观点出发求解轧辊温度场体系增加的热量使其内能发生变化，温度升高根据式(4.9)求出轧辊温度变化设温度变化率为，则在时间内，体积为的物体储能变化与温度变化之间的关系为： (4.21)式中，为体系内储能量的变化；为物质密度，；c为物质比热，；为物体体积，；为时间间隔，；为温度变化率，由于对工作辊划分了单元并且温度变化的时间很短，可把上式写为： (4.22)将式(4.9)带入式(4.22)可得： (4.23)下面将以式(4.23)为基础，计算轧辊各层各段的温度变化。

轧辊单元热流如图4.4所示，流入单元热量为正，流出单元热量为负图4.4 单元热流图Fig. 4.4 Heat flow graph of elements(1) 0层0段轧辊温度变化量计算数学模型 (4.24)式中，ftew[0]为第0层的单位热容，J/K；deltat为时间间隔，deltat=1.0s2) 0段1~3层轧辊温度计算数学模型 (4.25)式中，ftew[i]为第i层的单位热容，J/K3) 0层1~9段轧辊温度变化量计算数学模型 (4.26)(4) 1~3层其余段轧辊温度变化量计算数学模型 (4.27)式中，deltat为时间间隔，s；tstep为温度变化步长，3℃；dtmax为轧辊外面两层各段温差的最大值，℃；ts1为上一次迭代计算轧辊温度变化累计时间，s；ts为轧辊温度变化的累计时间，s在轧辊温度场模拟计算中，用式(4.23)先求得轧辊在一定时间(给定初始值为1s)内的温度变化，再以3.0℃为温度变化步长计算整个轧制时间内轧辊温度变化并判断轧辊温度变化累计时间ts是否超过整个轧制时间timex若，则进行下一个温度步长计算，计算该步长内单元实际温度，直到，其中，时间步长为最外两层各段单元温度变化最大时所需时间；当剩余时间不足一个温度步长时，则计算在剩余时间结束后单元温差，直到轧辊温度计算结束。

4.2.4 相关热传递参数的确定轧制时冷却液与工作辊间的热传递系数，对工作辊温度场及热凸度影响很大图4.5为国内某1250mm热连轧厂精轧机组工作辊冷却系统图轧制过程中的轧辊热传递主要包括与空气的自然对流换热和与冷却水的强制对流换热空气自然对流换热系数为=5~15W/(m2·K)图4.5 轧辊冷却水系统分布图Fig. 4.5 Diagram of work roll cooling system distribution因为轧辊中部经常与高温轧件接触，为降低轧辊热凸度，通过增大冷却水流量的方法强化轧辊中部冷却图4.6为轧辊冷却水喷嘴的布置图，图中各喷嘴倾斜角度皆为15°a）（b）（c）（d）图4.6 轧辊冷却水喷嘴分布图Fig. 4.6 Cooling water nozzles distribution(a) 入口上喷水集管 (b) 入口下喷水集管 (c) 出口上喷水集管 (d) 出口下喷水集管轧辊温度场和热凸度计算精确与否，主要依赖于边界条件的确定研究表明：轧辊与冷却水的对流换热系数与水的压力、喷嘴与轧辊表面的距离、喷射角度以及喷射区内水量密度密切相关4.2.4.1 工作辊表层各段与冷却水热传递系数的确定工作辊表层单元的轴向热传递系数并不是恒定的，在有钢轧制和无钢空过的状态下其热传递系数不同。

1) 轧辊不与带钢接触或无钢时与冷却水间的对流换热系数轧辊中部 (4.28)式中，aa0为水冷状态及单元位置对换热系数的影响系数；aa1为水冷状态及单元位置对换热系数的影响增益系数；rwc为轧辊与水的对流换热系数，W/(m2·K)；其值见表3.1；twd为轧辊直径，m；aew为轧辊单元长度，mm；fht为修正系数，取值为0.6；alp为带钢咬入角，弧度；dh为带钢压下量，mm；he为入口带钢厚度，mm；ha为。