专利名称：板带热轧过程中s型变步长法预测瞬态温度场方法
技术领域：
本发明属于轧制技术领域，特别涉及一种板带热轧过程s型变步长法预测瞬态温度场方法。
背景技术：
在板带热轧制过程中，从板坯出炉到层流冷却结束，轧件瞬态温度场分布计算对轧制过 程轧件组织演变、力学性能和板形等有重要影响。过去经常采用现场测试和数学模型进行确 定，由于受实际条件和模型限制影响，降低了温度预测精度和温度分布信息。有限元法作为 一种有效的数值分析方法已成为预测板板坯空冷过程瞬态温度场的重要手段。目前在有限元 法预测板带热轧过程的瞬态温度场中，往往采用较多单元划分或者较小的定时间步长来提高 计算精度，获得更为详尽的温度分布信息。然而单元划分过多和时间步长过小增加了计算时 间，降低了计算效率。
发明内容
为了克服上述方法计算速度慢，计算时间和计算精度不能同时满足的缺点，本发明提出一
种s型变步长法预测瞬态温度场方法，其目的是保证板带温度预测精度的情况下，縮短计算
时间，提高计算效率。
实现本发明目的技术解决方案如下
1、 采集轧制过程数据，包括轧制参数，材料热物性参数，单元划分信息。 轧制参数初始时间，轧制时间，轧件宽度，轧件厚度，初始温度，轧件周围介质温度 材料热物性参数热传导系数，黑度，比热，密度
单元划分信息宽度单元数和厚度单元数
2、 根据单元划分数据、轧件宽度和厚度尺寸建立有限元分析模型(见图1)，然后进行单 元节点编号、确定换热边界和计算节点坐标。
3、 根据轧制过程实际条件和轧制阶段确定换热系数A和内热源强度々。热
轧过程包括空冷阶段，除磷阶段，轧制阶段。整个轧制过程不同阶段的换热系数和内热源强 度计算如下
(1)热轧板带在空冷过程中，^值为零；其表面换热方式主要为辐射和自然对流，换热 系数A通过式(1)和式(2)计算
朋巧.e.(r《)(r2+r。,2) (i)
(2)
其中/汉为辐射系数，o"为Stefan-Boltzman常数，o" = 5.67xl0—8『/(附2 .《4); s为黑度系数，s与温度的关系式为s = 0.125(771000)2 - 0.38(771000) +1.1 。
热轧板带在空冷过程的对流方式为自然对流，其表达式为 、
^^=i.ix{(r-r。;,)/6^ 上表面 i/q = 0.5 x《r - r。， 下表面
其中r(《)为板带表面温度；7^为环境温度；6表示板带宽度。
(2) 热轧板带在高压水除磷过程中，^值为零；表面换热方式为强迫对流和侧面辐射，换 热系数通过式(1)和式(3)计算
<formula>formula see original document page 5</formula> (3)
其中w (丄/min'm2)为水流密度；T(《)板带表面温度。
(3) 在轧制过程中，板带与轧辊之间接触换热是主要热损失方式，忽略塑性变形和摩擦 做功，4 = 0。接触换热系数与轧制压力有关。换热系数A通过式(4)计算
<formula>formula see original document page 5</formula> (4)
式中A(M尸a)-轧制压力
4、 利用有限元基本原理，计算四边形等参单元的形函数iV、 B矩阵和雅克 比矩阵J。
5、 以二维热传导基本方程为基础，利用欧拉方程建立等效泛函，确定温度场求解的系统 方程。
以热力学第一定律为依据建立无内热源强度的二维热传导的微分方程为
<formula>formula see original document page 5</formula>
(5)
其中:
r-瞬时温度(/0
p-材料密度(&g/m3) c-材料比热 卜时间(s)
A:-热传导系数(『/(m'i0) 利用欧拉方程将二维热传导问题方程(5)变为等效泛涵
<formula>formula see original document page 5</formula>
(6)
根据热传导问题的变分原理，对泛函式(6)求一阶偏导数并置零，得到温度求解的系统方 <formula>formula see original document page 6</formula>
式中[尺r]-温度刚度矩阵，fc]=ffcW]+[At)]); [A]-变温矩阵，[《3]=tkW]; W-
常数项列式，fp卜l;fy"; {"-温度列式；五-单元总数；上标e表示每个单元。 'av, 3乂 抓，抓/
<formula>formula see original document page 6</formula>
对每个单元来说，刚度矩阵、变温矩阵和常数项可以通过下式求解
抓，抓7 + ，，抓乂.
<formula>formula see original document page 6</formula>
必
其中:
"热传导系数(『/(m'x));
/ -材料密度ag/m3); c-材料比热(々(&.尺))； //-换热系数； 7V-形函数；
/， 7'节点编号。
6、利用二点向后差分格式，将系统方程转化为瞬态温度场求解的线性方程组c 将系统方程(7)中的温度对时间偏导数表示为二点向后差分格式<formula>formula see original document page 6</formula>
将时间向后差分格式(8)带入系统方程式(7)得到温度场求解的线性方程组
(<formula>formula see original document page 6</formula>
7、利用S型变步长法设定线性方程组(9)中的时间步长Af ， S型变时间步长的基本原理
和数学模型。得到最终的温度场求解的线性方程组-
〔fl-6xexp(-cx,)) 3 , (a-6xexp(-cxf")
8、 采用一维变带宽存储法求解线性方程组(10)可以获得板带轧制过程瞬态温度场。
9、 计算时间在初始时间基础上加上时间步长，根据计算时间和设定轧制时间判断程序计
算是否结束，如果计算时间大于轧制时间退出，否则继续返回计算。根据轧制时间判断求解
过程是否结束，并将S型变步长法应用到瞬态温度场的有限元求解过程。 其中S型变步长原理与模型如下
S型变步长基本思想是在温度场计算初期采用较小的时间步长，然后时间步长虽时间增
加逐渐增加，最后成为一恒定值，变时间步长的S型曲线见图3所示。由于热分析初期，边
界条件和物体本身温差较大，温度变化较为剧烈，因此采用较小的迭代时间步长，得到瞬时
温度变化信息，提高计算精度。随着热传导过程的进行，整个系统逐渐趋于热平衡，采用小
的时间步长会导致计算时间急剧增加，因此增大时间步长，以縮短计算时间。时间步长不能
一味的增加，时间步长过大，不仅会降低温度求解精度，而且一定程度上不能满足瞬时温度
场求解的要求。S型变步长不仅能够保证计算精度，而且能够极大的縮短计算时间。
为了更方便的实现S型变步长在温度场有限元计算中的作用，对S型变步长曲线进行了 回归，得出如下模型-
= a-6xexp(-ex/1" (11)
其中t为时间，a和b主要控制迭代时间步长的初始值和终值的大小，A0 = a-6， Wc^a， c和d主要控制S型曲线的变化率。
本发明在传统算法的基础上通过采用S型变步长法预测瞬态温度场，在不影响计算效 率的情况下有效地克服了振荡现象，保证了计算的稳定性。
图l本发明方法有限元分析模型图， 图2本发明方法具体操作软件流程图，
图3本发明方法S型变步长曲线图，
图4实施例中第一组方案采用本发明方法与采用定步长计算结果结果比较图， 图5实施例中第二组方案采用本发明方法与采用定步长计算结果结果比较图， 图6实施例中第三组方案采用本发明方法与采用定步长计算结果结果比较图， 图7实施例中第一组方案采用本发明方法与采用定步长计算结果结果比较图。 图中i为单元编号，j为节点编号，H为厚度，W为宽度，1为换热边界，a为主要 控制迭代时间步长的初始值的大小，b为主要控制迭代时间步长的终值的大小，A:0-a-&，Hoc" ， A-定步长平均温度，B-定步长表面温度，C-定步长中心温度，D-变步长平均温度， E-变步长表面温度，F-变步长中心温度。
具体实施例方式
选择一热轧板坯从出炉到除磷前的空冷过程温度变化进行预测分析，比较定步长法和 变步长法对板坯温度计算结果和计算时间的影响。
例计算条件见表l
表l
热物性参数 数值
截面尺寸(宽x厚)1800mm X 220mm
热传导率A: = 30『/(w.iC)
比热c = 620
密度
黑度0.8
波尔滋曼常数5.67 XI0-8
初始温度1200°C
环境温度30。C
冷却时间300s
计算分析了四组方案，每组方案的单元划分、定时间步长和变步长分别为:
方案1单元划分20X10;定时间步长At=ls;变步长设定a=10.0， b = 9.0， C二O.OOOl，
d=2.5
方案2单元划分30X10，定时间步长At-ls;变步长设定a=10.0, b=9.0， c=0.001， d=2.5
方案3单元划分为40X20，定时间步长At二0.5s;变步长设定a=15.0，b二 14.5，c二0.0001 ， d=2.5
方案4单元划分为50X20，定时间步长At二0.1s;变步长设定a=15.0，b = 14.9，c=0.0001 ， d=2.5
图4， 5， 6， 7所示为各个方案下定步长和变步长温度变化比较，可以看出在整个冷却时 间内，定步长温度计算结果和变步长温度计算结果基本相同，只有心部温度有一定差别，但 是相对误差不超过0.2%，所以变步长和定步长温度计算精度相同。通过表2可以看出，变时 间步长的计算时间远远小于定时间步长，这种计算速度的优越性随着单元数目的增加更加明
显，变时间步长在保证计算精度的情况下，极大的縮短了计算时间，提高了计算效率。
表2
方案一方案二 ，方案二方案四
tl(ms) t2(ms)tl(ms) t2(ms)tl(ms) t2(ms)tl(ms) t2(ms)
287 40564 716350 4卯482卯 1475
本发明方法针对板坯出炉空冷过程温度变化进行了求解，变步长计算结果和定步长温度
计算基本相同，相对误差不超过0.2%，而变步长法縮短了计算时间，提高了计算速度。
权利要求
1、板带热轧过程中S型变步长法预测瞬态温度场方法，其特征包括以下步骤(1)采集轧制过程数据，包括轧制参数，材料热物性参数，单元划分信息轧制参数初始时间，轧制时间，轧件宽度，轧件厚度，初始温度，轧件周围介质温度，材料热物性参数热传导系数，黑度，比热，密度，单元划分信息宽度单元数和厚度单元数；(2)根据单元划分数据、轧件宽度和厚度尺寸建立有限元分析模型，进行单元节点编号、确定换热边界和计算节点坐标；(3)根据不同轧制过程，确定边界换热系数热轧板带在空冷过程中，其表面换热方式主要为辐射和自然对流，辐射系数表述为HR＝σ·ε·(T+Tair)(T2+Tair2)式中HR为辐射系数σ为Stefan-Boltzman常数σ＝5.67×10-8W/(m2·K4)ε为黑度系数，ε与温度的关系式为ε＝0.125(T/1000)2-0.38(T/1000)+1.1热轧板带在高压水除磷过程中，主要换热方式为强迫对流和侧面辐射，辐射计算方法同上，对流系数表达式为HCW＝124.7×w0.663×10-0.00147(T-273.16)式中w为水流密度T板带表面温度在轧制过程中，板带与轧辊之间接触换热是主要热损失方式，接触换热系数表达式为IHTC＝695pm-34400(W/m2K)式中pm-轧制压力；(4)利用有限元基本原理，计算四边形等参单元的形函数N、B矩阵和雅克比矩阵J；(5)利用欧拉方程建立等效泛函，确定温度场求解的系统方程以热力学第一定律为依据建立无内热源强度的二维热传导的微分方程为
2、按照权利要求1所述的板带热轧过程中S型变步长法预测瞬态温度场方法，其特征在 于所述的S型变步长法数学模型如下<formula>formula see original document page 3</formula>式中a为控制迭代时间步的初始值大小b为控制迭代时间步的终值大小c和d主要控制S型曲线的变化率t为时间。
全文摘要
板带热轧过程中S型变步长法预测瞬态温度场方法，S型变步长的时间步数学模型为Δt＝a-b×exp(-c×t<sup>d</sup>)，利用S型变步长方法进行瞬态温度场有限元求解不仅提高了计算精度，满足了初始温度场设定要求，而且大大缩短了计算时间。这种方法在有限元法求解瞬态温度场时解决了计算精度低、速度慢的问题，提高了计算效率。
文档编号G06F17/50GK101178747SQ20071015898
公开日2008年5月14日 申请日期2007年12月18日 优先权日2007年12月18日
发明者刘相华, 李长生, 梅瑞斌 申请人:东北大学