CN101178747B - 板带热轧过程中s型变步长法预测瞬态温度场方法 - Google Patents
板带热轧过程中s型变步长法预测瞬态温度场方法 Download PDFInfo
- Publication number
- CN101178747B CN101178747B CN2007101589827A CN200710158982A CN101178747B CN 101178747 B CN101178747 B CN 101178747B CN 2007101589827 A CN2007101589827 A CN 2007101589827A CN 200710158982 A CN200710158982 A CN 200710158982A CN 101178747 B CN101178747 B CN 101178747B
- Authority
- CN
- China
- Prior art keywords
- partiald
- integral
- temperature
- heat
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Control Of Metal Rolling (AREA)
Abstract
板带热轧过程中S型变步长法预测瞬态温度场方法,S型变步长的时间步数学模型为:Δt=a-b×exp(-c×td),利用S型变步长方法进行瞬态温度场有限元求解不仅提高了计算精度,满足了初始温度场设定要求,而且大大缩短了计算时间。这种方法在有限元法求解瞬态温度场时解决了计算精度低、速度慢的问题,提高了计算效率。
Description
技术领域
本发明属于轧制技术领域,特别涉及一种板带热轧过程S型变步长法预测瞬态温度场方法。
背景技术
在板带热轧制过程中,从板坯出炉到层流冷却结束,轧件瞬态温度场分布计算对轧制过程轧件组织演变、力学性能和板形等有重要影响。过去经常采用现场测试和数学模型进行确定,由于受实际条件和模型限制影响,降低了温度预测精度和温度分布信息。有限元法作为一种有效的数值分析方法已成为预测板板坯空冷过程瞬态温度场的重要手段。目前在有限元法预测板带热轧过程的瞬态温度场中,往往采用较多单元划分或者较小的定时间步长来提高计算精度,获得更为详尽的温度分布信息。然而单元划分过多和时间步长过小增加了计算时间,降低了计算效率。
发明内容
为了克服上述方法计算速度慢,计算时间和计算精度不能同时满足的缺点,本发明提出一种S型变步长法预测瞬态温度场方法,其目的是保证板带温度预测精度的情况下,缩短计算时间,提高计算效率。
实现本发明目的技术解决方案如下:
1、采集轧制过程数据,包括:轧制参数,材料热物性参数,单元划分信息。
轧制参数:初始时间,轧制时间,轧件宽度,轧件厚度,初始温度,轧件周围介质温度材料热物性参数:热传导系数,黑度,比热,密度
单元划分信息:宽度单元数和厚度单元数
2、根据单元划分数据、轧件宽度和厚度尺寸建立有限元分析模型(见图1),然后进行单元节点编号、确定换热边界和计算节点坐标。
HR=σ·ε·(T+Tair)(T2+Tair 2) (1)
其中:HR为辐射系数,σ为Stefan-Boltzman常数,σ=5.67×10-8W/(m2·K4);ε为黑度系数,ε与温度的关系式为ε=0.125(T/1000)2-0.38(T/1000)+1.1。
热轧板带在空冷过程的对流方式为自然对流,其表达式为:
其中:T(K)为板带表面温度;Tair为环境温度;b表示板带宽度。
(2)热轧板带在高压水除磷过程中,值为零;表面换热方式为强迫对流和侧面辐射,换热系数h通过式(1)和式(3)计算:
HCW=124.7×w0.663×10-0.00147(T-7316) (3)
其中w(L/min·m2)为水流密度;T(K)板带表面温度。
IHTC=695pm-34400(W/m2K) (4)
式中:pm(MPa)-轧制压力
4、利用有限元基本原理,计算四边形等参单元的形函数N、B矩阵和雅克比矩阵J。
5、以二维热传导基本方程为基础,利用欧拉方程建立等效泛函,确定温度场求解的系统方程。
以热力学第一定律为依据建立无内热源强度的二维热传导的微分方程为:
其中:
T-瞬时温度(K)
ρ-材料密度(kg/m3)
c-材料比热(J/(kg·K))
t-时间(s)
k-热传导系数(W/(m·K))
利用欧拉方程将二维热传导问题方程(5)变为等效泛涵:
根据热传导问题的变分原理,对泛函式(6)求一阶偏导数并置零,得到温度求解的系统方程:
K3ij=∫∫SρcNiNjdS
pi=∫LhT∞NtdL
对每个单元来说,刚度矩阵、变温矩阵和常数项可以通过下式求解:
其中:
k-热传导系数(W/(m·K));
ρ-材料密度(kg/m3);
c-材料比热(J/(kg·K));
h-换热系数;
N-形函数;
i,j节点编号。
6、利用二点向后差分格式,将系统方程转化为瞬态温度场求解的线性方程组。
将系统方程(7)中的温度对时间偏导数表示为二点向后差分格式:
将时间向后差分格式(8)带入系统方程式(7)得到温度场求解的线性方程组:
7、利用S型变步长法设定线性方程组(9)中的时间步长Δt,S型变时间步长的基本原理和数学模型。得到最终的温度场求解的线性方程组:
8、采用一维变带宽存储法求解线性方程组(10)可以获得板带轧制过程瞬态温度场。
9、计算时间在初始时间基础上加上时间步长,根据计算时间和设定轧制时间判断程序计算是否结束,如果计算时间大于轧制时间退出,否则继续返回计算。根据轧制时间判断求解过程是否结束,并将S型变步长法应用到瞬态温度场的有限元求解过程。
其中S型变步长原理与模型如下:
S型变步长基本思想是在温度场计算初期采用较小的时间步长,然后时间步长虽时间增加逐渐增加,最后成为一恒定值,变时间步长的S型曲线见图3所示。由于热分析初期,边界条件和物体本身温差较大,温度变化较为剧烈,因此采用较小的迭代时间步长,得到瞬时温度变化信息,提高计算精度。随着热传导过程的进行,整个系统逐渐趋于热平衡,采用小的时间步长会导致计算时间急剧增加,因此增大时间步长,以缩短计算时间。时间步长不能一味的增加,时间步长过大,不仅会降低温度求解精度,而且一定程度上不能满足瞬时温度场求解的要求。S型变步长不仅能够保证计算精度,而且能够极大的缩短计算时间。
为了更方便的实现S型变步长在温度场有限元计算中的作用,对S型变步长曲线进行了回归,得出如下模型:
Δt=a-b×exp(-c×td) (11)
其中:t为时间,a和b主要控制迭代时间步长的初始值和终值的大小,k0=a-b,k1∝a,c和d主要控制S型曲线的变化率。
本发明在传统算法的基础上通过采用S型变步长法预测瞬态温度场,在不影响计算效率的情况下有效地克服了振荡现象,保证了计算的稳定性。
附图说明
图1本发明方法有限元分析模型图,
图2本发明方法具体操作软件流程图,
图3本发明方法S型变步长曲线图,
图4实施例中第一组方案采用本发明方法与采用定步长计算结果结果比较图,
图5实施例中第二组方案采用本发明方法与采用定步长计算结果结果比较图,
图6实施例中第三组方案采用本发明方法与采用定步长计算结果结果比较图,
图7实施例中第一组方案采用本发明方法与采用定步长计算结果结果比较图。
图中:i为单元编号,j为节点编号,H为厚度,W为宽度,1为换热边界,a为主要控制迭代时间步长的初始值的大小,b为主要控制迭代时间步长的终值的大小,k0=a-b,k1∝a,A-定步长平均温度,B-定步长表面温度,C-定步长中心温度,D-变步长平均温度,E-变步长表面温度,F-变步长中心温度。
具体实施方式
选择一热轧板坯从出炉到除磷前的空冷过程温度变化进行预测分析,比较定步长法和变步长法对板坯温度计算结果和计算时间的影响。
例:计算条件见表1
表1
计算分析了四组方案,每组方案的单元划分、定时间步长和变步长分别为:
方案1单元划分:20×10;定时间步长:Δt=1s;变步长设定:a=10.0,b=9.0,c=0.0001,
d=2.5
方案2单元划分:30×10,定时间步长Δt=1s;变步长设定:a=10.0,b=9.0,c=0.001,
d=2.5
方案3单元划分为40×20,定时间步长Δt=0.5s;变步长设定:a=15.0,b=14.5,c=0.0001,
d=2.5
方案4单元划分为50×20,定时间步长Δt=0.1s;变步长设定:a=15.0,b=14.9,c=0.0001,
d=2.5
图4,5,6,7所示为各个方案下定步长和变步长温度变化比较,可以看出在整个冷却时间内,定步长温度计算结果和变步长温度计算结果基本相同,只有心部温度有一定差别,但是相对误差不超过0.2%,所以变步长和定步长温度计算精度相同。通过表2可以看出,变时间步长的计算时间远远小于定时间步长,这种计算速度的优越性随着单元数目的增加更加明显,变时间步长在保证计算精度的情况下,极大的缩短了计算时间,提高了计算效率。
表2
本发明方法针对板坯出炉空冷过程温度变化进行了求解,变步长计算结果和定步长温度计算基本相同,相对误差不超过0.2%,而变步长法缩短了计算时间,提高了计算速度。
Claims (2)
1.板带热轧过程中S型变步长法预测瞬态温度场方法,其特征包括以下步骤:
(1)采集轧制过程数据,包括:轧制参数,材料热物性参数,单元划分信息轧制参数:初始时间,轧制时间,轧件宽度,轧件厚度,初始温度,轧件周围介质温度,
材料热物性参数:热传导系数,黑度,比热,密度,
单元划分信息:宽度单元数和厚度单元数;
(2)根据单元划分数据、轧件宽度和厚度尺寸建立有限元分析模型,进行单元节点编号、确定换热边界和计算节点坐标;
(3)根据不同轧制过程,确定边界换热系数
热轧板带在空冷过程中,其表面换热方式主要为辐射和自然对流,辐射系数表述为:
HR=σ·ε·(T+Tair)(T2+Tair 2)
式中: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)以二维热传导基本方程为基础,利用欧拉方程建立等效泛函,确定温度场求解的系统方程;
以热力学第一定律为依据建立无内热源强度的二维热传导的微分方程为:
其中:
T-瞬时温度(K)
ρ-材料密度(kg/m3)
c-材料比热(J/(kg·K))
t-时间(s)
k-热传导系数(W/(m·K))
利用欧拉方程将二维热传导问题方程(5)变为等效泛涵:
根据热传导问题的变分原理,对泛函式(6)求一阶偏导数并置零,得到温度求解的系统方程:
对每个单元来说,刚度矩阵、变温矩阵和常数项可以通过下式求解:
其中:
k-热传导系数(W/(m·K));
ρ-材料密度(kg/m3);
c-材料比热(J/(kg·K));
h-换热系数;
N-形函数;
i,j节点编号;
(6)利用二点向后差分格式,将系统方程转化为瞬态温度场的线性方程组
(7)利用S型变步长法设定线性方程组中的时间步长Δt,得到最终的温度场求解的线性方程组
式中:a为控制迭代时间步的初始值大小
b为控制迭代时间步的终值大小
c和d主要控制S型曲线的变化率
t为时间
(8)采用一维变带宽存储法求解线性方程组板带轧制过程瞬态温度场。
2.按照权利要求1所述的板带热轧过程中S型变步长法预测瞬态温度场方法,其特征在于所述的S型变步长法数学模型如下:
Δt=a-b×exp(-c×td)
式中:a为控制迭代时间步的初始值大小
b为控制迭代时间步的终值大小
c和d主要控制S型曲线的变化率
t为时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101589827A CN101178747B (zh) | 2007-12-18 | 2007-12-18 | 板带热轧过程中s型变步长法预测瞬态温度场方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101589827A CN101178747B (zh) | 2007-12-18 | 2007-12-18 | 板带热轧过程中s型变步长法预测瞬态温度场方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101178747A CN101178747A (zh) | 2008-05-14 |
CN101178747B true CN101178747B (zh) | 2010-06-09 |
Family
ID=39405001
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101589827A Expired - Fee Related CN101178747B (zh) | 2007-12-18 | 2007-12-18 | 板带热轧过程中s型变步长法预测瞬态温度场方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101178747B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8935945B2 (en) * | 2008-11-19 | 2015-01-20 | Toshiba Mitsubishi-Electic Industrial Systems Corporation | Control system |
CN101559511B (zh) * | 2009-05-22 | 2011-12-28 | 清华大学 | 一种以温度为控制变量的焊接数值模拟计算方法 |
CN103028615B (zh) * | 2012-11-29 | 2014-12-10 | 一重集团大连设计研究院有限公司 | 一种预测带钢热连轧过程温度演变的方法 |
CN103852094B (zh) * | 2012-12-07 | 2016-05-18 | 中国核动力研究设计院 | 一种瞬态高温对仪表影响的判断方法 |
CN104462735A (zh) * | 2015-01-16 | 2015-03-25 | 济南璘康光电子信息技术有限公司 | 一种同轴封装tosa温度分布的仿真方法 |
CN105550463B (zh) * | 2015-03-13 | 2018-10-12 | 东北大学 | 钢板电磁感应加热过程温度场的预测方法 |
CN105414205B (zh) * | 2015-12-17 | 2017-04-12 | 东北大学 | 一种基于plc的钢板温度在线预测方法 |
CN106295053A (zh) * | 2016-08-20 | 2017-01-04 | 三峡大学 | 一种基于自适应时间步长的瞬变电磁‑温度场耦合计算方法 |
CN107066737B (zh) * | 2017-04-14 | 2019-05-17 | 北京科技大学 | 一种预测热轧过程板带温度场的二维交替差分方法 |
CN112808781B (zh) * | 2019-11-15 | 2022-08-12 | 中冶华天工程技术有限公司 | 螺纹钢棒轧件轧制过程温度计算的方法 |
CN114861134B (zh) * | 2022-07-08 | 2022-09-06 | 四川大学 | 一种计算水滴运动轨迹的步长确定方法及存储介质 |
-
2007
- 2007-12-18 CN CN2007101589827A patent/CN101178747B/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
包仲南等.带钢热连轧机工作辊瞬态温度场的有限元仿真.北京科技大学学报21 1.1999,21(1),60-63. * |
孔祥伟等.轧辊温度场及轴向热凸度有限元计算.钢铁研究学报12.2000,1251-54. * |
杨利坡等.热连轧工作辊三维瞬态温度场数值模拟.燕山大学学报28 5.2004,28(5),380-383. * |
杨利坡等.热连轧轧辊瞬态温度场研究.钢铁40 10.2005,40(10),38-41. * |
Also Published As
Publication number | Publication date |
---|---|
CN101178747A (zh) | 2008-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101178747B (zh) | 板带热轧过程中s型变步长法预测瞬态温度场方法 | |
CN100495411C (zh) | 一种预测热轧过程板带温度场的有限元方法 | |
Dialameh et al. | Natural convection from an array of horizontal rectangular thick fins with short length | |
Giangaspero et al. | Application of the entropy generation minimization method to a solar heat exchanger: A pseudo-optimization design process based on the analysis of the local entropy generation maps | |
Amanowicz et al. | Validation of CFD model for simulation of multi-pipe earth-to-air heat exchangers (EAHEs) flow performance | |
CN100545849C (zh) | 一种有限元求解轧制过程温度场的集中热容矩阵方法 | |
Korzeń et al. | Modeling of transient response of a plate fin and tube heat exchanger | |
He et al. | Performance prediction of an air-cooled steam condenser using UDF method | |
Lin et al. | A two-dimensional fin efficiency analysis of combined heat and mass transfer in elliptic fins | |
Dubovsky et al. | Analytical model of a PCM-air heat exchanger | |
CN101221416B (zh) | 热轧过程在线计算板带温度的有限元方法 | |
Kundu et al. | An analytical prediction for performance and optimization of an annular fin assembly of trapezoidal profile under dehumidifying conditions | |
CN103761370B (zh) | 一种板带热轧过程表面换热系数的预测方法 | |
Adhikari et al. | Optimizing rectangular fins for natural convection cooling using CFD | |
Zhang et al. | Simulation and experimental investigation of the wavy fin-and-tube intercooler | |
Huang et al. | An impingement heat sink module design problem in determining simultaneously the optimal non-uniform fin widths and heights | |
Kong et al. | Impacts of geometric structures on thermo-flow performances of plate fin-tube bundles | |
Zhao et al. | Numerical study on airside thermal-hydraulic performance of rectangular finned elliptical tube heat exchanger with large row number in turbulent flow regime | |
Chen et al. | Study of heat-transfer characteristics on the fin of two-row plate finned-tube heat exchangers | |
Wu et al. | Modeling the performance of the indirect dry cooling system in a thermal power generating unit under variable ambient conditions | |
Peng et al. | Analysis of heat transfer and flow characteristics over serrated fins with different flow directions | |
Xu et al. | Experimental study on heat transfer performance improvement of wavy finned flat tube | |
Zhang et al. | Numerical analysis of thermal-hydraulic characteristics on serrated fins with different attack angles and wavelength to fin length ratio | |
Campo et al. | Turbulent natural convection in an air-filled isosceles triangular enclosure | |
Zhang et al. | A favorable face velocity distribution and a V-frame cell for power plant air-cooled condensers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100609 Termination date: 20121218 |