CN110348174B - 一种加热炉的钢坯温度计算方法 - Google Patents
一种加热炉的钢坯温度计算方法 Download PDFInfo
- Publication number
- CN110348174B CN110348174B CN201910726081.6A CN201910726081A CN110348174B CN 110348174 B CN110348174 B CN 110348174B CN 201910726081 A CN201910726081 A CN 201910726081A CN 110348174 B CN110348174 B CN 110348174B
- Authority
- CN
- China
- Prior art keywords
- billet
- furnace
- radiation
- exchange area
- unit
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control Of Heat Treatment Processes (AREA)
- Waste-Gas Treatment And Other Accessory Devices For Furnaces (AREA)
Abstract
本发明涉及一种加热炉的钢坯温度计算方法,属于加热炉优化控制技术领域。该方法通过辐射交换面积计算钢坯表面热流,包括以下步骤:S1:将加热炉沿炉长方向分为若干个模型段;S2:首先计算各段炉墙、钢坯、炉气之间的辐射交换面积,再通过辐射交换面积计算导来辐射系数,进而获得钢坯表面热流;S3:采用差分的方法计算钢坯截面温度分布。本发明采用辐射交换面积的方法计算钢坯表面热流,把黑体材料的影响引入到了辐射热流计算中,更适用于有黑体材料的加热炉的模型计算。
Description
技术领域
本发明属于加热炉优化控制技术领域,涉及一种加热炉的钢坯温度计算方法。
背景技术
加热炉是轧钢生产的重要设备,其加热钢坯的质量直接制约着成品质量;同时,加热炉也是轧钢生产的主要耗能设备,其能耗约占整个轧钢工序的70%。为了能降低加热炉能耗,人们研究开发了很多节能材料和节能技术。在加热炉上采用黑体材料,就是一种有效可行的降低加热炉能耗的措施,近年来应用越来越多。
在加热炉炉墙上安装黑体原件,可以扩大炉膛面积,提高炉墙黑度,从而提高加热效率,达到节能降耗的目的。
在传统的加热炉模型中,钢坯表面热流的计算大都采用导来辐射系数法,在导来辐射系数的计算中,并没有涉及到炉墙黑度这一参数,也就是对于普通加热炉和黑体加热炉,传统加热炉模型中对导来辐射系数的计算是相同的,因此,有必要针对黑体加热炉,把炉墙黑度引入到模型计算中,从而体现黑体材料对钢坯温度的影响。
目前,国内就黑体材料提出了一些专利申请。专利CN201310220035.1提出一种黑体材料制备成分及由黑体材料制成的节能辐射杯;专利201811011279.8介绍了一种黑体辐射节能涂料及其制备方法;专利201720574059.0发明了一种高温辐射黑体元件,重点介绍了黑体原件的结构。在加热炉模型专利方面,大都是适用于普通加热炉的计算方法,尚无针对黑体加热炉模型方面的专利。
发明内容
有鉴于此,本发明的目的在于提供一种加热炉的钢坯温度计算方法,该方法利用辐射交换面积计算导来辐射系数,进而计算钢坯表面热流,最后采用一维差分获得钢坯截面温度分布。该方法把炉墙黑度引入到了炉膛辐射换热计算中,不但适用于普通加热炉,更适用于采用黑体材料的节能加热炉。
为达到上述目的,本发明提供如下技术方案:
一种加热炉的钢坯温度计算方法,该方法包括以下步骤:
S1:根据加热炉尺寸和沿炉长方向,将加热炉分为若干个模型段;
S2:在每个模型段内,计算钢坯、炉墙和炉气的辐射直接交换面积;根据直接交换面积计算钢坯、炉墙和炉气三者之间的辐射全交换面积,再进一步求得导来辐射系数,计算钢坯表面热流;
S3:求解一维差分方程,计算钢坯内部温度,获得钢坯厚度方向上的温度分布。
进一步,所述步骤S2具体为:
S21:在每个模型段内,根据段内钢坯、炉墙、炉气的尺寸和位置关系,计算钢坯、炉墙、炉气三者之间的辐射直接交换面积:
式中,为炉气对钢坯的辐射直接交换面积,单位为m2;/>为炉气对炉墙的辐射直接交换面积,单位为m2;/>为炉墙对钢坯的辐射直接交换面积,单位为m2;/>为炉墙对炉墙的辐射直接交换面积,单位为m2;Fs为模型段内钢坯的表面积,单位为m2;εg为炉气黑度;为炉墙对钢坯的角系数;
S22:根据辐射直接交换面积计算钢坯、炉墙、炉气之间的全交换面积:
式中,为炉气对炉墙的辐射全交换面积单位为,m2;/>为炉气对钢坯的辐射全交换面积,单位为m2;/>为炉墙对钢坯的辐射全交换面积,单位为m2;Fs为模型段内钢坯的表面积,单位为m2;Fw为模型段内炉墙的表面积,单位为m2;εg为炉气黑度;εw为炉墙黑度;εs为钢坯黑度;ρw为炉墙的反射率;ρs为钢坯的反射率;
S23:根据辐射全交换面积计算导来辐射系数:
εgws=εgs+εg-w-s
式中,εgws为导来辐射系数;
S24:根据导来辐射系数计算钢坯表面热流:
式中,Qgws为钢坯表面热流,单位为W;σ为斯蒂芬-波尔兹曼常数,为5.67×10-8W/(m2K4);Tg为炉气温度,单位为K;Ts为钢坯表面温度,单位为K。
进一步,所述步骤S3具体为:
划分网格,求解钢坯一维厚度方向的差分方程,钢坯一维隐式差分方程为:
式中,qu为钢坯上表面热流密度,单位为W/m2;qb为钢坯下表面热流密度,单位为W/m2;Fo为傅立叶数;λ为钢坯导热系数,单位为W/(mk);△τ为时间步长,单位为s;△y为空间步长,单位为m;
式中,c为钢坯比热,单位为J/(kg.K)。
本发明的有益效果在于:本发明提出了一种加热炉的钢坯温度计算方法,把炉墙黑度引入到了炉膛辐射换热中,有利于分析采用黑体材料后,加热炉节能降耗的效果,不仅适用于普通加热炉,更适用于采用黑体材料的节能加热炉。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为本发明的流程图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
其中,附图仅用于示例性说明,表示的仅是示意图,而非实物图,不能理解为对本发明的限制;为了更好地说明本发明的实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
本发明实施例的附图中相同或相似的标号对应相同或相似的部件;在本发明的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本发明的限制,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
本实例中选取的加热炉有效长度为25.5m,炉膛内宽10.7m,共设置预热段、加热段、均热段三个段,在两侧炉墙和炉顶安装了黑体原件。
如图1所示,本发明的方法包括以下步骤:
S1:根据加热炉尺寸,沿炉长方向,将加热炉分为25个模型段,前24段每1m一段,最后一段1.5m。
S2:在每个模型段内,利用钢坯、炉墙、炉气之间的辐射交换面积计算钢坯表面热流,具体分为以下步骤:
S21:在每个模型段内,根据段内炉墙、钢坯、炉气的尺寸和位置关系,计算三者之间的辐射直接交换面积。
式中,为炉气对钢坯的辐射直接交换面积,m2;/>为炉气对炉墙的辐射直接交换面积,m2;/>为炉墙对钢坯的辐射直接交换面积,m2;/>为炉墙对炉墙的辐射直接交换面积,m2;Fs为模型段内钢坯的表面积,m2;εg为炉气黑度;/>为炉墙对钢坯的角系数。
S22:根据辐射直接交换面积计算炉墙、钢坯、炉气之间的全交换面积。
式中,为炉气对炉墙的辐射全交换面积,m2;/>为炉气对钢坯的辐射全交换面积,m2;/>为炉墙对钢坯的辐射全交换面积,m2;Fs为模型段内钢坯的表面积,m2;Fw为模型段内炉墙的表面积,m2;εg为炉气黑度;εw为炉墙黑度;εs为钢坯黑度;ρw为炉墙的反射率;ρs为钢坯的反射率。
S23:根据辐射全交换面积计算导来辐射系数。
εgws=εgs+εg-w-s
式中,εgws为导来辐射系数。
钢坯信息如下表所示。
表1钢坯信息表
钢种 | 钢种 | 钢坯厚度/mm | 钢坯长度/mm | 钢坯宽度/mm | 目标温度/℃ |
Q345 | Q345 | 220 | 10000 | 1200 | 1250 |
根据以上步骤,计算得到增加黑体原件前后导来辐射系数分别为0.5023和0.5917,说明黑体材料增加了炉墙对钢坯的辐射传热。
表2计算结果
参数 | 普通加热炉 | 黑体加热炉 |
εgws | 0.5023 | 0.5917 |
S24:根据导来辐射系数计算钢坯表面热流。
式中,Qgws为钢坯表面热流,W;σ为斯蒂芬-波尔兹曼常数,5.67×10-8W/(m2K4);Tg为炉气温度,K;Ts为钢坯表面温度,K。
S3:求解一维差分方程,计算钢坯内部温度,获得钢坯厚度方向上的温度分布,钢坯一维隐式差分方程为:
式中,qu为钢坯上表面热流密度,W/m2;qb为钢坯下表面热流密度,W/m2;Fo为傅立叶数;λ为钢坯导热系数,W/(mk);△τ为时间步长,s;△y为空间步长,m。
式中,c为钢坯比热,J/(kg.K)。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种加热炉的钢坯温度计算方法,其特征在于:该方法包括以下步骤:
S1:根据加热炉尺寸和沿炉长方向,将加热炉分为若干个模型段;
S2:在每个模型段内,计算钢坯、炉墙和炉气的辐射直接交换面积;根据直接交换面积计算钢坯、炉墙和炉气三者之间的辐射全交换面积,再进一步求得导来辐射系数,计算钢坯表面热流;
S3:求解一维差分方程,计算钢坯内部温度,获得钢坯厚度方向上的温度分布;
所述步骤S2具体为:
S21:在每个模型段内,根据段内钢坯、炉墙、炉气的尺寸和位置关系,计算钢坯、炉墙、炉气三者之间的辐射直接交换面积:
式中,为炉气对钢坯的辐射直接交换面积,单位为m2;/>为炉气对炉墙的辐射直接交换面积,单位为m2;/>为炉墙对钢坯的辐射直接交换面积,单位为m2;/>为炉墙对炉墙的辐射直接交换面积,单位为m2;Fs为模型段内钢坯的表面积,单位为m2;εg为炉气黑度;/>为炉墙对钢坯的角系数;
S22:根据辐射直接交换面积计算钢坯、炉墙、炉气之间的全交换面积:
式中,为炉气对炉墙的辐射全交换面积,单位为m2;/>为炉气对钢坯的辐射全交换面积,单位为m2;/>为炉墙对钢坯的辐射全交换面积,单位为m2;Fs为模型段内钢坯的表面积,单位为m2;Fw为模型段内炉墙的表面积,单位为m2;εg为炉气黑度;εw为炉墙黑度;εs为钢坯黑度;ρw为炉墙的反射率;ρs为钢坯的反射率;
S23:根据辐射全交换面积计算导来辐射系数:
εgws=εgs+εg-w-s
式中,εgws为导来辐射系数;
S24:根据导来辐射系数计算钢坯表面热流:
式中,Qgws为钢坯表面热流,单位为W;σ为斯蒂芬-波尔兹曼常数,为5.67×10-8W/(m2K4);Tg为炉气温度,单位为K;Ts为钢坯表面温度,单位为K;
所述步骤S3具体为:
划分网格,求解钢坯一维厚度方向的差分方程,钢坯一维隐式差分方程为:
式中,qu为钢坯上表面热流密度,单位为W/m2;qb为钢坯下表面热流密度,单位为W/m2;Fo为傅立叶数;λ为钢坯导热系数,单位为W/(mk);△τ为时间步长,单位为s;△y为空间步长,单位为m;
式中,c为钢坯比热,单位为J/(kg.K)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910726081.6A CN110348174B (zh) | 2019-08-07 | 2019-08-07 | 一种加热炉的钢坯温度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910726081.6A CN110348174B (zh) | 2019-08-07 | 2019-08-07 | 一种加热炉的钢坯温度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110348174A CN110348174A (zh) | 2019-10-18 |
CN110348174B true CN110348174B (zh) | 2023-06-16 |
Family
ID=68184199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910726081.6A Active CN110348174B (zh) | 2019-08-07 | 2019-08-07 | 一种加热炉的钢坯温度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110348174B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110918655A (zh) * | 2019-11-30 | 2020-03-27 | 宝钢特钢韶关有限公司 | 一种精细化加热控制方法 |
CN111763819A (zh) * | 2020-08-03 | 2020-10-13 | 重庆赛迪热工环保工程技术有限公司 | 一种低温辊底式热处理炉钢板加热的控制方法 |
CN112697290B (zh) * | 2020-11-27 | 2022-03-22 | 攀钢集团攀枝花钢铁研究院有限公司 | 钢坯入炉温度的检测方法 |
CN113850004A (zh) * | 2021-10-08 | 2021-12-28 | 重庆赛迪热工环保工程技术有限公司 | 一种钢板叠装时的温度计算方法 |
CN115065710B (zh) * | 2022-04-29 | 2023-07-25 | 燕山大学 | 一种加热炉智慧温控pc端及移动端远程云测控系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102982243A (zh) * | 2012-11-30 | 2013-03-20 | 中冶南方(武汉)威仕工业炉有限公司 | 一种加热炉总括热吸收率计算方法 |
CN102994730A (zh) * | 2012-11-30 | 2013-03-27 | 中冶南方(武汉)威仕工业炉有限公司 | 一种加热炉内钢坯温度跟踪的方法 |
CN106906352A (zh) * | 2017-03-30 | 2017-06-30 | 重庆赛迪热工环保工程技术有限公司 | 一种轧钢加热炉钢坯混装时的加热方法 |
CN109926675A (zh) * | 2019-03-29 | 2019-06-25 | 安徽双桦热交换系统有限公司 | 一种汽车散热器钎焊工艺的nb连续式钎焊芯体加热方法 |
CN110057712A (zh) * | 2019-04-24 | 2019-07-26 | 首钢集团有限公司 | 一种氧化烧损率测定系统及方法 |
CN110081711A (zh) * | 2019-04-22 | 2019-08-02 | 东北大学 | 一种双圆台式黑体辐射节能元件及其制备、使用方法 |
-
2019
- 2019-08-07 CN CN201910726081.6A patent/CN110348174B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102982243A (zh) * | 2012-11-30 | 2013-03-20 | 中冶南方(武汉)威仕工业炉有限公司 | 一种加热炉总括热吸收率计算方法 |
CN102994730A (zh) * | 2012-11-30 | 2013-03-27 | 中冶南方(武汉)威仕工业炉有限公司 | 一种加热炉内钢坯温度跟踪的方法 |
CN106906352A (zh) * | 2017-03-30 | 2017-06-30 | 重庆赛迪热工环保工程技术有限公司 | 一种轧钢加热炉钢坯混装时的加热方法 |
CN109926675A (zh) * | 2019-03-29 | 2019-06-25 | 安徽双桦热交换系统有限公司 | 一种汽车散热器钎焊工艺的nb连续式钎焊芯体加热方法 |
CN110081711A (zh) * | 2019-04-22 | 2019-08-02 | 东北大学 | 一种双圆台式黑体辐射节能元件及其制备、使用方法 |
CN110057712A (zh) * | 2019-04-24 | 2019-07-26 | 首钢集团有限公司 | 一种氧化烧损率测定系统及方法 |
Non-Patent Citations (1)
Title |
---|
"加热炉过程建模与控制仿真研究";王利涛;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20150515;正文11-34页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110348174A (zh) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110348174B (zh) | 一种加热炉的钢坯温度计算方法 | |
CN106906352B (zh) | 一种轧钢加热炉钢坯混装时的加热方法 | |
CN104827245B (zh) | 一种轴承套圈的锻造工艺 | |
CN107604283B (zh) | 一种生产大规格Ti55531合金饼坯的方法 | |
CN106969632B (zh) | 一种多层箱式混合加热炉 | |
CN105239001B (zh) | 油汀用冷轧钢板及其制备方法 | |
CN108251608A (zh) | 一种高温轴承钢制轴承套圈整形加工方法 | |
CN111763819A (zh) | 一种低温辊底式热处理炉钢板加热的控制方法 | |
CN105002346A (zh) | 回火炉热风循环系统 | |
WO2008094093A1 (en) | Method for the heat treatment of extended steel products | |
CN201924061U (zh) | 一种退火炉中的喷管 | |
CN104220618B (zh) | 具有优异材料均匀度的高碳热轧钢板及其制造方法 | |
CN104384359A (zh) | 一种热冲压成形坯料的立体快速加热装置及加热方法 | |
CN103639248B (zh) | 一种大型环类零件热处理后椭圆变形的校正方法 | |
US10619925B2 (en) | Heating device for hot stamping | |
CN110193536A (zh) | 一种柔性热辊弯成形设备及方法 | |
CN106399651B (zh) | 一种薄壁大直径模具整体真空淬火工艺 | |
CN109518122A (zh) | 薄壁大尺寸非对称回转类钛合金零件离子氮化控制方法 | |
CN204325475U (zh) | 一种卧式多用离子氮化炉 | |
CN103343288B (zh) | 一种缠绕气瓶用无缝钢管及其制造方法 | |
CN106854686A (zh) | 一种超高强度钢薄壁壳体淬火变形的整形方法 | |
CN105369029A (zh) | 台车式钢球回火炉 | |
CN106399664B (zh) | 一种环形加热炉加热工艺优化方法 | |
CN204874656U (zh) | 回火炉热风循环系统 | |
Fang et al. | Batch annealing model for cold rolled coils and its application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |