CN110826172B - 一种高炉炉缸横截面侵蚀边界二维快速计算方法 - Google Patents
一种高炉炉缸横截面侵蚀边界二维快速计算方法 Download PDFInfo
- Publication number
- CN110826172B CN110826172B CN201910897604.3A CN201910897604A CN110826172B CN 110826172 B CN110826172 B CN 110826172B CN 201910897604 A CN201910897604 A CN 201910897604A CN 110826172 B CN110826172 B CN 110826172B
- Authority
- CN
- China
- Prior art keywords
- erosion
- calculating
- calculation
- blast furnace
- boundary
- 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
- 230000003628 erosive effect Effects 0.000 title claims abstract description 59
- 238000004364 calculation method Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 20
- 238000007711 solidification Methods 0.000 claims abstract description 9
- 230000008023 solidification Effects 0.000 claims abstract description 9
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 8
- 238000013213 extrapolation Methods 0.000 claims description 4
- 229910052742 iron Inorganic materials 0.000 claims description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000012937 correction Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005034 decoration Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009851 ferrous metallurgy Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
- C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
- C22B5/00—General methods of reducing to metals
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Manufacturing & Machinery (AREA)
- Materials Engineering (AREA)
- Mechanical Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Waste-Gas Treatment And Other Accessory Devices For Furnaces (AREA)
Abstract
本发明公开了一种高炉炉缸横截面侵蚀边界二维快速计算方法，包括以下步骤：高炉炉缸横截面为圆环，采用固定计算网格的方法绘制计算域结构化网格，得到初始凝固线位置；通过一系列假设的侵蚀线的位置得到所有靠近圆心圈热电偶的计算温度；利用线性插值方法初步计算出初始的内壁侵蚀线，并根据内壁侵蚀线计算出所有热电偶所在位置对应的温度；将所有热电偶所在位置对应的温度与实际热电偶温度值多次比较并修正，直到误差范围满足要求。本发明有助于缓解现有技术存在的侵蚀检测迭代方法耗时过长、计算量大、不利于工程上快速预估获得结果且普通的计算硬件无法满足计算要求的技术问题。
Description
技术领域
本发明涉及高炉炉缸侵蚀检测技术领域，尤其是涉及一种高炉炉缸横截面侵蚀边界二维快速计算方法。
背景技术
高炉是钢铁冶金生产中关键设备之一，是炼铁工艺过程的主要设备载体。高炉的使用寿命主要由高炉炉缸部位限制。在影响高炉安全的因素中，内衬侵蚀是最重要的一个因素。由于炉缸内衬侵蚀部位被熔融铁水覆盖，导致炉缸内衬侵蚀状态不能直接观测，这是侵蚀预测的一大难题。为了防止高炉炉缸被烧穿，建立描述高炉炉缸热过程的数学模型，进行炉缸内衬的侵蚀分析计算或诊断技术是非常必要的。在高炉服役中后期，需要判断炉缸内衬侵蚀形貌和厚度，科学安排护炉和大修期，确保炉缸结构的安全。同时也要求尽可能充分地利用现有的热工测量条件参数，作深入综合的数据挖掘，最大限度地真实反映炉缸内衬工作状态。
近年来，在传热学的基础上，根据温度、热流量等热工测量参数所建立的多种数学模型在不断发展中。一般来说，针对高炉炉缸侵蚀问题的研究分为正解模型和逆解模型。正解模型多根据侵蚀机理去预测侵蚀位置，但是由于炉缸内衬侵蚀机理非常复杂，涉及到氧化作用、热应力、机械冲刷等等许多方面，产量和介质也非常复杂，所以从正向问题去预测侵蚀线的位置非常困难。而逆解模型则是直接利用服役高炉炉缸的温度和热流量等现有可测得的热工参数来确定内衬侵蚀状态。即在已知部分几何边界、导热系数、温度、热流密度等参数情况下，求解未知部分几何边界。因此，逆解模型比正解模型更具有实用性，是目前分析高炉炉缸侵蚀问题的主要方法。
依据傅立叶传热原理建立逆解为预测侵蚀边界的可靠途径之一。传统经典的一维计算模型易于掌握，但是精度不够高；三维模型则尚未实用。因此，二维计算模型相对一维能大大提高分析精度，又有很高的技术价值，是目前最实用的方法。其难点在于不能经过一次计算获得可信的侵蚀边界，需要迭代、搜索、逼近多次计算方能获得。常用的二维迭代方法耗时过长，每次迭代需要重新进行网格划分，极大的耗费资源，不利于工程上快速预估获得结果，且普通的计算硬件无法满足计算要求。
发明内容
本发明的目的是提供一种高炉炉缸横截面侵蚀边界二维快速计算方法，以解决现有技术存在的侵蚀检测迭代方法耗时过长、计算量大、不利于工程上快速预估获得结果且普通的计算硬件无法满足计算要求的技术问题。
为实现上述目的，本发明采用的技术方案为：
一种高炉炉缸横截面侵蚀边界二维快速计算方法，包括以下步骤：
步骤一，高炉炉缸横截面为圆环，采用固定计算网格的方法绘制计算域结构化网格，得到初始凝固线位置；
步骤二，通过一系列假设的侵蚀线的位置得到所有靠近圆心圈热电偶的计算温度；
步骤三，利用线性插值方法初步计算出初始的内壁侵蚀线，并根据内壁侵蚀线计算出所有热电偶所在位置对应的温度；
步骤四，将所有热电偶所在位置对应的温度与实际热电偶温度值多次比较并修正，直到误差范围满足要求。
所述步骤一中，计算域结构化网格边界外层由外圈热电偶Pi,out所在位置确定，计算域结构化网格内边界根据一维傅里叶定律反推，由所有的热电偶组Pi依次计算，得到不同的初始凝固线位置；取计算结果中侵蚀厚度的最大值记为δ0，再取△δ作为凝铁层厚度，令δmax＝δ0+△δ，此δmax即为固定网格所描述的计算域的厚度值，由此确定计算域结构化网格的内边界位置。
所述步骤二中，将高炉炉缸横截面的圆心到内壁面的半径值记为r0，由r0开始取间隔△r，获得一系列r1，r2，r3，……，根据所取间隔△r不停将侵蚀线的位置向外层进行假设推进；此时外层热电偶Pi,out的温度值和坐标为已知量；针对不同半径值ri，其坐标和温度值也为已知量，温度为铁水凝固温度1150℃；由上述所取一系列半径值ri和已知的边界条件，分别对每一半径值ri所对应的热电偶Pi,in所在位置对应的温度进行计算，记为Tical,in，直到所有的Tical,in>Ti,in，其中，Ti,in为该热电偶Pi,in所在位置实际温度值。
所述步骤三中，将步骤二得到的一系列半径值ri和计算值Tical,in进行整理，通过线性插值的方法初步预估推测，得到初始的内壁侵蚀线；由初始内壁侵蚀线计算得到一组新的Tical,in，将与实际的Ti,in计算误差范围，当温度数值误差在1.5％以内时，认为此处壁厚即为实际值；对误差范围大于1.5％的进行修正计算。
所述步骤四中，当根据侵蚀线所计算出的所有热电偶Pi,in所在位置对应的温度Tical,in与实际值Ti,in的误差范围都在1.5％以内时，认为此条线即为所求侵蚀线。
有益效果：本发明提供的在固定计算域网格下的高炉炉缸侵蚀边界横截面二维快速计算方法，有助于缓解现有技术存在的侵蚀检测迭代方法耗时过长、计算量大、不利于工程上快速预估获得结果且普通的计算硬件无法满足计算要求的技术问题。
附图说明
图1为本发明的计算域网格示意图；
图2为本发明侵蚀线计算修正示意图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
以下结合附图对本发明的实施例进行详细说明，但是本发明可以根据权利要求限定和覆盖的多种不同方式实施。
本发明的一种高炉炉缸横截面侵蚀边界二维快速计算方法包括以下步骤：
第一步：绘制计算域结构化网格。为了方便展示截取四分之一圆环进行说明。本发明二维方法针对高炉炉缸横截面，其横截面为圆环。为了防止生成凹形边界，导致负网格的产生，同时减小计算量，节约计算时间，本发明采取固定计算网格的方法。
计算域网格边界外层由外圈热电偶Pi,out所在位置确定，热电偶位置如图1所示，每组热电偶都由靠近内壁和外壁的两个热电偶组成。计算域内边界根据一维傅里叶定律反推，由所有的热电偶组Pi依次计算，得到不同的初始凝固线位置。取计算结果中侵蚀厚度的最大值记为δ0，再取△δ作为凝铁层厚度，令δmax＝δ0+△δ，此δmax即为固定网格所描述的计算域的厚度值，在图1中也标注出了δmax所在位置。由此就完成了确定计算域内边界位置。
第二步：如图2所示，将圆心到内壁面的半径值记为r0，由r0开始取间隔△r，获得一系列r1，r2，r3，……，即根据所取间隔△r不停将侵蚀线的位置向外层进行假设推进。此时外层热电偶Pi,out的温度值和坐标为已知量。针对不同半径值ri，其坐标和温度值也为已知量，温度即为铁水凝固温度1150℃。
由上述所取一系列半径值ri和已知的边界条件，分别对每一半径值ri所对应的热电偶Pi,in所在位置对应的温度进行计算，记为Tical,in，直到所有的Tical,in>Ti,in，其中，Ti,in为该热电偶Pi,in所在位置实际温度值。
第三步：将上述一系列半径值ri和计算值Tical,in进行整理，通过线性插值的方法初步预估推测，可以得到初始的内壁侵蚀线。
由初始内壁侵蚀线计算得到一组新的Tical,in，将与实际的Ti,in计算误差范围，当温度数值误差在1.5％以内时，认为此处壁厚即为实际值；对误差范围大于1.5％的进行修正计算。在图2中，浅色区域部分为最新计算获得侵蚀线和外边界所包含的区域，深色线则为上一次计算的得到的侵蚀线。深色区域为两次侵蚀线和外边界包含区域的改变区域。
第四步：当根据侵蚀线所计算出的所有热电偶Pi,in所在位置对应的温度Tical,in与实际值Ti,in的误差范围都在1.5％以内时，认为此条线即为所求侵蚀线。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种高炉炉缸横截面侵蚀边界二维快速计算方法，其特征在于：包括以下步骤：
步骤一，高炉炉缸横截面为圆环，采用固定计算网格的方法绘制计算域结构化网格，得到初始凝固线位置；
步骤二，通过一系列假设的侵蚀线的位置得到所有靠近圆心圈热电偶的计算温度；
其中，将高炉炉缸横截面的圆心到内壁面的半径值记为r 0 ，由r 0 开始取间隔Δr，获得一系列r 1 ，r 2 ，r 3 ，……，根据所取间隔Δr不停将侵蚀线的位置向外层进行假设推进；此时外层热电偶P i,out 的温度值和坐标为已知量；针对不同半径值r i ，其坐标和温度值也为已知量，温度为铁水凝固温度1150℃；由上述所取一系列半径值r i 和已知的边界条件，分别对每一半径值r i 所对应的热电偶P i,in 所在位置对应的温度进行计算，记为T ical,in ，直到所有的T ical,in > T i,in ，其中，T i,in 为该热电偶P i,in 所在位置实际温度值；
步骤三，利用线性插值方法初步计算出初始的内壁侵蚀线，并根据内壁侵蚀线计算出所有热电偶所在位置对应的温度；
其中，将步骤二得到的一系列半径值r i 和计算值T ical,in 进行整理，通过线性插值的方法初步预估推测，得到初始的内壁侵蚀线；由初始内壁侵蚀线计算得到一组新的T ical,in ，将与实际的T i,in 计算误差范围，当温度数值误差在1.5%以内时，认为此处壁厚即为实际值；对误差范围大于1.5%的进行修正计算；
步骤四，将所有热电偶所在位置对应的温度与实际热电偶温度值多次比较并修正，直到误差范围满足要求。
2.根据权利要求1所述的高炉炉缸横截面侵蚀边界二维快速计算方法，其特征在于：所述步骤一中，计算域结构化网格边界外层由外圈热电偶P i,out 所在位置确定，计算域结构化网格内边界根据一维傅里叶定律反推，由所有的热电偶组P i 依次计算，得到不同的初始凝固线位置；取计算结果中侵蚀厚度的最大值记为δ 0 ，再取Δδ作为凝铁层厚度，令δ max =δ 0 +Δδ，此δ max 即为固定网格所描述的计算域的厚度值，由此确定计算域结构化网格的内边界位置。
3.根据权利要求1所述的高炉炉缸横截面侵蚀边界二维快速计算方法，其特征在于：所述步骤四中，当根据侵蚀线所计算出的所有热电偶P i,in 所在位置对应的温度T ical,in 与实际值T i,in 的误差范围都在1.5%以内时，认为此条线即为所求侵蚀线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910897604.3A CN110826172B (zh) | 2019-09-23 | 2019-09-23 | 一种高炉炉缸横截面侵蚀边界二维快速计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910897604.3A CN110826172B (zh) | 2019-09-23 | 2019-09-23 | 一种高炉炉缸横截面侵蚀边界二维快速计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110826172A CN110826172A (zh) | 2020-02-21 |
| CN110826172B true CN110826172B (zh) | 2021-09-17 |
Family
ID=69548205
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910897604.3A Active CN110826172B (zh) | 2019-09-23 | 2019-09-23 | 一种高炉炉缸横截面侵蚀边界二维快速计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110826172B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101275829A (zh) * | 2007-03-29 | 2008-10-01 | 上海梅山钢铁股份有限公司 | 一种测量高炉炉衬侵蚀厚度方法 |
| WO2011089621A2 (en) * | 2010-01-19 | 2011-07-28 | Aditya Birla Science & Technology Co. Ltd. | A soft sensor based on-line decision system for metallurgical processes |
| CN105005632A (zh) * | 2014-04-17 | 2015-10-28 | 宝山钢铁股份有限公司 | 多层耐火砖炉墙结构的高炉炉缸侵蚀预测方法 |
| CN108220513A (zh) * | 2017-12-26 | 2018-06-29 | 柳州钢铁股份有限公司 | 高炉节能护炉方法 |
| CN108517384A (zh) * | 2018-05-04 | 2018-09-11 | 中南大学 | 一种高炉炉缸侧壁的监测方法 |
| CN108707712A (zh) * | 2018-05-24 | 2018-10-26 | 山西太钢不锈钢股份有限公司 | 一种高炉炉缸侧壁残厚的判断方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102876825A (zh) * | 2012-09-13 | 2013-01-16 | 江苏永钢集团有限公司 | 一种高炉炉缸侵蚀线位置的判定方法 |
| CN108090293B (zh) * | 2017-12-26 | 2021-04-30 | 国家超级计算天津中心 | 一种高炉炉缸炉底侵蚀包络面确定方法 |
-
2019
- 2019-09-23 CN CN201910897604.3A patent/CN110826172B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101275829A (zh) * | 2007-03-29 | 2008-10-01 | 上海梅山钢铁股份有限公司 | 一种测量高炉炉衬侵蚀厚度方法 |
| WO2011089621A2 (en) * | 2010-01-19 | 2011-07-28 | Aditya Birla Science & Technology Co. Ltd. | A soft sensor based on-line decision system for metallurgical processes |
| CN105005632A (zh) * | 2014-04-17 | 2015-10-28 | 宝山钢铁股份有限公司 | 多层耐火砖炉墙结构的高炉炉缸侵蚀预测方法 |
| CN108220513A (zh) * | 2017-12-26 | 2018-06-29 | 柳州钢铁股份有限公司 | 高炉节能护炉方法 |
| CN108517384A (zh) * | 2018-05-04 | 2018-09-11 | 中南大学 | 一种高炉炉缸侧壁的监测方法 |
| CN108707712A (zh) * | 2018-05-24 | 2018-10-26 | 山西太钢不锈钢股份有限公司 | 一种高炉炉缸侧壁残厚的判断方法 |
Non-Patent Citations (3)
| Title |
|---|
| Model for Fast Computation of Blast Furnace Hearth Erosion and Buildup Profiles;Johnny Brannbacka,et al.;《industrial & Engineering Chemistry Research》;20080910;第40卷(第20期);7793–7801 * |
| The thermophysical bases of monotoring of the fireproff lining wear in the blast furnace hearth;A.N. Dmitriev,et al.;《Defect & Diffusion Forum》;20170120;第370卷;113-119 * |
| 二维法在高炉炉缸内衬侵蚀模型中的应用研究;陈良玉等;《钢铁》;20080831;第43卷(第8期);12-15 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110826172A (zh) | 2020-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110765623A (zh) | 一种高炉炉缸纵截面侵蚀边界二维快速计算方法 | |
| CN103614498B (zh) | 一种高炉炉壁三维温度场重建方法及计算机监控系统 | |
| Yuan et al. | Development of an improved CBR model for predicting steel temperature in ladle furnace refining | |
| CN102095516A (zh) | 一种废杂铜冶炼炉温测量方法 | |
| Zagaria et al. | Monitoring erosion and skull profile in blast furnace hearth | |
| CN110826172B (zh) | 一种高炉炉缸横截面侵蚀边界二维快速计算方法 | |
| CN108090293B (zh) | 一种高炉炉缸炉底侵蚀包络面确定方法 | |
| CN103966377A (zh) | 高炉炉温在线检测采集系统及在线检测方法 | |
| CN110781566B (zh) | 一种炉缸凝铁层计算方法、系统、存储介质及电子终端 | |
| CN112800612A (zh) | 一种基于元胞自动机的镍基合金腐蚀层动态演变分析方法 | |
| WO2023070694A1 (zh) | 一种判定高炉最优气流分布的方法、电子设备和存储介质 | |
| Botnikov et al. | Development of the metal temperature prediction model for steel-pouring and tundish ladles used at the casting and rolling complex | |
| CN113111549B (zh) | 高炉炉缸浇注修复后的侵蚀模型建模方法及建模系统 | |
| CN103728333A (zh) | 双样杯球铁或蠕铁炉前快速分析方法及装置 | |
| JP2008063593A (ja) | 容器壁の厚み推定方法、装置、コンピュータプログラム | |
| CN116167669B (zh) | 一种基于电力消费回归的碳排放评估方法 | |
| Longwell | A graphical method for solution of freezing problems | |
| WO2024001047A1 (zh) | 温度分布图库的建立方法和晶圆表面温度的获取方法 | |
| Ding et al. | Spatial Differences in Water Quality and Spatial Autocorrelation Analysis of Eutrophication in Songhua Lake | |
| CN114996997A (zh) | 一种炉缸纵截面侵蚀边界计算的有限元方法 | |
| CN112930465A (zh) | 用于监测高炉的耐火炉衬的磨损的方法 | |
| CN111961776B (zh) | 一种高炉炉缸内衬角落区域热电偶位置映射方法 | |
| Prasolov et al. | Development of a simulation model of the heat transfer process in the hot-blast stove checkerwork | |
| JP4559708B2 (ja) | 反応容器内の状態を評価する評価方法、操業管理方法、評価装置、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 | |
| CN114418464B (zh) | 一种基于案例推理的二步法冶金流程能源评价方法及系统 |
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 |