冷却壁的安全稳定工作是高炉长寿运行的基础, 尤其在炉腰和炉腹高热负荷区域, 良好的冷却制度和工艺条件可使冷却壁表面形成稳定渣皮, 对冷却壁起保护作用, 延长高炉寿命.近年来国内外学者从传热学角度出发, 对炉腹区铜冷却壁热结构就传热特性、挂渣能力等方面进行了研究[1-4].Zhong等[5]对埋纯铜管式铸铜冷却壁进行了系统的可行性分析.郭光胜等[6]基于有限元法研究了冷却比表面积对冷却壁的传热影响.石琳等[7]对埋纯铜管式冷却壁进行热态试验研究, 得出煤气温度对铜冷却壁本体应力分布影响显著.Huo等[8]对轧制铜冷却壁水管断裂破损机理进行论述研究, 并提出修复方法.Kawaoka等[9]对铜冷却壁纯铜材料进行高温下磨损试验, 得出铜冷却壁材料在超过200 ℃环境中强度明显降低.近年来埋纯铜管式铜冷却壁已在国内外多个大型高炉上成功应用[10].然而, 目前铜冷却壁已出现水管开裂和热面破损等现象, 而对此现象产生的机理缺乏定量化的分析, 外加以往对铜冷却壁的研究只局限于圆形冷却水管或某单一因素的影响, 缺乏各因素间相互的影响, 故本文基于传热学理论对目前高炉炉腹中普遍采用的铜冷却壁进行传热特性和力学性能分析, 在优化冷却水管管型的基础上, 以铜冷却壁本体厚度、冷却比表面积、冷却水水速和冷却水温度为比较序列, 以高炉正常工况下冷却壁本体最大热应力为指标, 定量化计算各参数与指标参数间的相互关联程度.结合响应面模型和NSGA-Ⅱ遗传算法以铜冷却壁本体最大热应力为优化目标, 对炉腹区热结构进行结构参数和长寿技术优化, 并就优化前后的炉腹结构传热模型进行对比分析, 验证了响应面法和NSGA-Ⅱ遗传算法相结合对高炉炉腹区热结构进行优化分析的有效性.
1 边界条件及有限元模型以某高炉炉腹区结构为对象进行分析, 由外到内分别为炉壳(含螺栓、定位销等)、填料层、铜冷却壁(含冷却水管, 本文高炉炉腹冷却壁为圆形水管)和镶砖.相邻冷却壁间及冷却壁与炉壳间缝隙由填料填充.冷却壁为“四进四出”埋纯铜管式铸铜冷却壁, 由于每块冷却壁的工作环境基本相同, 为减少计算量, 取单块冷却壁+侧隙填料+填料层+炉壳的扇区作为计算模型.冷却壁上下端面具有相同位移,而左右侧面环向位移为零, 忽略上下左右间的传热, 则单块冷却壁具有周期对称力学结构特征.工况下水温差造成的冷却壁上下游温度梯度对其耦合特性的影响较小, 冷却水取平均温度, 可减少计算量, 柱坐标系下单块冷却壁扇区结构模型如图 1所示.
炉壳与周围空气间综合对流传热系数hs的表达式为
(1) |
式中, t1为炉壳周围环境温度, ℃.
根据吴俐俊等[11]基于边界条件替代法的研究数据, 确定500~1 248 ℃范围内热炉气与镶砖热面间的综合对流传热系数hx随温度变化规律如式(2)所示.
(2) |
式中, t2为热炉气温度, ℃.
铜冷却壁具有高导热性、挂渣稳定、承受高负荷和长寿命等优点, 埋纯铜管式铸铜冷却壁避免了焊接工序, 彻底消除了气隙层.某高炉炉腹区各层结构尺寸参数见表 1.
为更好地比较各管型对冷却壁传热特性的影响, 以圆管为参考标准, 考虑两种情况下异型管冷却壁传热特性, 即与圆管同横截面积和同周长.
2.1 同横截面积时冷却壁传热特性分析各异型管当量直径DG计算式为
(3) |
式中:A为过流断面面积; S为湿周长度.
各异型管当量直径见表 2, 矩形管与椭圆管短长轴比例b/a均先按0.6进行计算.
在给定流场下流体的惯性力与黏性力对比关系可用雷诺数Re来表征, 表达式为
(4) |
式中:v为流体速度; υ为流体运动黏度; ρ为流体密度; μ为流体动力黏度; l为特征长度.
雷诺数是判断流体流动状态的依据, 对于圆管, 当Re≤2 300时为层流流动; 当Re≥104时为紊流流动; 当2 300 < Re < 104时为过渡流流动.以圆管为例, 对壁面温度不变的光滑或垂直长管且流体和壁面温差不大时, 紊流状态冷却水对流换热系数可用狄特斯和玻尔特特征关联式[12]表示:
(5) |
式中:n为指数, 当流体加热时n=0.4.
联立式(4)、式(5)和努赛尔数Nu=h·DG/k, 可得紊流状态下圆管内冷却水对流换热系数表达式:
(6) |
据此分别求得各管型内冷却水对流换热系数h的表达式见表 3.
冷却水在各管型中流动时泵的功率相同, 可利用强化传热评价准则作为各管型强化传热评价标准, 计算式为
(7) |
式中:Nui和Nu0分别为异型管与圆管的努赛尔数; fi和f0分别为异型管及圆管的阻力系数.圆管阻力系数可由Blasius公式求得.
热炉气温度为1 200 ℃.图 2为横截面积相同, 冷却水流速3 m/s条件下高炉炉腹区各管型铜冷却壁本体最高温度tmax图.
从图中计算结果可看出, 各异型管与同横截面积的圆管相比, tmax均有小幅度降低, 其中矩形管(b/a=0.6)及三圆管冷却壁温度下降最为明显.图 3为矩形管不同b/a时tmax及η的变化规律, 随着b/a的不断增大, tmax不断增大, η不断减小, 当b/a大于0.6时, tmax急剧增大, η降幅趋于平缓.从传热强化评价指标来考虑, 各水管同横截面积时最佳管型为b/a不大于0.6的矩形管, 此时η基本上维持在0.78以上.考虑到加工难度及矩形管宽度过大影响冷却壁径向受力等问题, 矩形管b/a在0.5~0.6范围内最为合适.图 4为矩形管(b/a=0.6)冷却壁炉腹区热结构温度场云图.
图 5为各管型同周长、冷却水流速3m/s条件下tmax计算结果, 由图可知同周长时椭圆管冷却壁温度下降幅度最大.图 6为椭圆管不同b/a下tmax和η分布规律, 从传热强化综合评价指标考虑, 各管型同周长时最佳管型应为0.55≤b/a≤0.6的椭圆管, 此时η维持在0.83以上.
分析图 2和图 5计算结果可得, 同横截面积时矩形管铜冷却壁与同周长情况下椭圆管铜冷却壁相比在传热特性方面差距不大, 在强化传热评价指标方面略低.综合考虑传热、流阻特性及矩形管在工作过程中易出现应力集中点等因素, 高炉炉腹区铜冷却壁冷却水管改为0.55≤b/a≤0.6的椭圆管最为合适, 同时椭圆管同原先圆管冷却壁相比可在获得相同冷却效果的同时节约大量冷却水用量, 提高效益.
3 高炉炉腹区结构长寿技术优化 3.1 响应面法简介响应面法基本思想是通过构造具有明确表达形式的多项式来表达隐式的功能函数.假设功能函数Z和随机变量Q=[Q1, Q2, …, QR]间关系如式(8)所示.对随机变量的N个样本进行随机抽样, 统计分析功能函数值(z1, z2, …, zN), 由最小二乘法理论拟合得到该系统功能函数[13].
(8) |
式中:a0, ai, aij(i=1, …, R; j=i, …, R)为函数方程待定系数, 共1+R+R(R+1)/2个.
采用矩阵法对每个随机变量取三个水平点, 根据Box-Behnken抽样法得出中心及边中点作为样本值点, 当随机变量分布符合某种规律时, 可用式(9)确定变量水平qs:
(9) |
式中:f(q)为随机变量概率密度函数; pn为水平点, 取p1=0.01, p2=0.5, p3=0.99.
对S个随机变量样本值进行数值模拟, 得S个输出点(z1, z2, …, zS), 通过回归分析可得式(10), 可确定功能函数(8)的函数关系式.设计方法为BBD法.
(10) |
热应力是影响高炉冷却壁寿命的关键因素之一, 从冷却壁本体热应力数值计算来分析各参数变化时的影响规律, 对高炉长寿有重要意义.本文以表 1中炉腹区结构为基准, 将圆管改为同周长且短长轴比例为0.6的椭圆管, 以工况下铜冷却壁本体最大热应力Smax(即Von Mises等效应力)为研究对象, 以铜冷却壁本体厚度、冷却比表面积、冷却水流速及冷却水温度为比较序列, 参数见表 4.
对Smax计算结果进行多元回归拟合分析, 表 5为分析结果.该模型决定系数R2=0.946 3和调整决定系数Radj2=0.892 5, 说明该模型计算结果对实验数据拟合较好, 可用于铜冷却壁本体最大热应力Smax随各参数变化的理论预测.各参数对Smax的显著性由F检验来判定, P值越小, F值越大, 则显著程度越高.由表 5可得, x2与x4对Smax影响尤为突出, 交互项中x1和x2的相互作用对Smax的影响比较大.限于篇幅只给出部分交互参数对Smax的影响规律结果如图 7所示.
各参数对Smax的影响程度由大到小排名依次为:x2, x4, x3, x1, 交互项中x1和x2相互作用较大.由响应面理论可知, 影响程度越大, 该因素对目标序列的影响越大.Smax随x1增加呈先增大后减小趋势, x2小时, 应适度减薄冷却壁厚度; Smax随x2的增大先迅速减小后趋于平缓, 冷却比表面积从0.7增加到1.1时, Smax下降非常明显, 可适当增大比表面积来降低冷却壁应力; x3从2 m/s增大到4 m/s时, Smax下降幅度较小, 这是由于冷却水对流换热热阻已不是冷却壁传热的限制环节, 提高冷却水水速, 无疑增加了水流量, 同时成倍提高了水流动阻力, 增大运行成本; 当x4较低时, 可减小水速降低成本, 水温较高时再增大水流速以减小冷却壁应力.实际操作中还应注意水温反复变化引起的本体反复热应力, 以免造成疲劳破坏, 故水温应以稳定为主.
据式(10)可得在高炉正常工作时铜冷却壁本体最大热应力Smax经回归分析后的响应面方程为
(11) |
在炉腹铜冷却壁响应面模型的基础上, 利用遗传算法对其本体最大热应力进行多参数下的高效全局寻优.以铜冷却壁优化后的整体质量mT不高于初始质量的95%及单位时间用水量Ls不高于初始单位时间用水量的95%作为约束条件.
优化设计的数学模型可表示为
(12) |
式中:ρ为铜冷却壁材料密度; m0为椭圆管铜冷却壁初始质量; L0为初始单位时间用水量.
取遗传算法交叉概率为0.6, 变异概率为0.01.优化前后各随机变量及目标函数计算结果见表 6.
为验证优化后结果的准确性, 将优化后的参数值代入有限元模型中, 结果如图 8所示.对比响应面方程计算结果198.42 MPa, 近似值误差为0.611%, 表明通过响应面模型和NSGA-Ⅱ遗传算法对优化高炉炉腹区结构参数及长寿技术具有较高的精度.
1) 从传热强化评价指标、加工难度及冷却壁径向承力等方面来考虑, 最佳管型为短长轴比例为0.55~0.6的椭圆管.炉腹区铜冷却壁冷却水管改为椭圆管后, 能够在保证较好冷却效果的同时节约冷却水用量, 有效提高企业效益.
2) 各参数对Smax的影响程度由大到小排名依次为:冷却比表面积、冷却水温度、冷却水速度、冷却壁本体厚度,交互项中冷却壁本体厚度和冷却比表面积相互作用比较大.设计炉腹区结构时可参考各因素对冷却壁应力场分布的影响, 保证冷却壁安全工作的同时对其结构进行优化.
3) 采用响应面法和NSGA-Ⅱ遗传算法相结合的思想对目标函数进行优化.优化后的冷却壁在传热特性方面得到提高, 本体质量、冷却水用量等方面得到降低, 验证了响应面法和遗传算法相结合的高炉炉腹区热结构参数优化方法的有效性.
[1] |
Guo H W, Fan D Q, Chang J, et al.Numerical simulation for BF hearth and bottom using finite difference method[C]//International Conference on Applied Informatics and Communication 2011.Berlin: Springer-Verlag, 2011: 449-455.
http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZNXX201012028076.htm |
[2] |
Wu L L, Zhou W G, Li P, et al.
Study on the equivalent convection coefficient of the hot surface of blast furnace stave[J]. Heat and Mass Transfer, 2007, 43(12): 1303–1309.
DOI:10.1007/s00231-006-0220-9 |
[3] |
Deng K, Wu L J.
Influence of structural parameters on steady-state highest temperature and thermal stress of BF cooling stave[J]. Iron Steel Vanadium Titanium, 2004, 25(4): 47–53.
|
[4] |
Smith M P, Fletcher J, Harvey R W, et al.
Blast furnace cooling stave design[J]. Steel Times International, 2016, 40(2): 33–37.
|
[5] |
Zhong Y, Lu X Z, Fan X M, et al.
Thermal simulation study on copper cooling stave of blast furnace[J]. Advanced Materials Research, 2015, 1088: 739–743.
DOI:10.4028/www.scientific.net/AMR.1088 |
[6] |
郭光胜, 张建良, 焦克新, 等.
冷却比表面积对高炉炉缸铸铁冷却壁传热的影响研究[J]. 铸造, 2016, 65(6): 542–548.
( Guo Guang-sheng, Zhang Jian-liang, Jiao Ke-xin, et al. Study of the effect of the ratio of cooling surface area on the heat transfer of the cast iron cooling stave in blast furnace hearth[J]. Foundry, 2016, 65(6): 542–548. DOI:10.3969/j.issn.1001-4977.2016.06.011 ) |
[7] |
石琳, 郭永茂, 曹福军.
铸铜冷却壁蠕变变形研究[J]. 内蒙古科技大学学报, 2013, 32(1): 42–45.
( Shi Lin, Guo Yong-mao, Cao Fu-jun. Research on cheep deformation of cast copper staves[J]. Journal of Inner Mongolia University of Science and Technology, 2013, 32(1): 42–45. ) |
[8] |
Huo Y, Chen L Y.Analysis on heat transfer characteristics for bosh copper cooling stave with shaped tube[C]//IOP Conference Series: Materials Science and Engineering.Xi'an: IOP Publishing, 2018: 022003.
|
[9] |
Kawaoka K, Anan K, Tsuda A, et al.
Latest blast furnace relining technology at Nippon Steel[J]. Shinnittetsu Giho, 2006, 384: 127–132.
|
[10] |
Qiu H Q, Chen Z D.
The forty years of vermicular graphite cast iron development in China(part Ⅰ)[J]. China Foundry, 2007, 4(2): 91–98.
|
[11] |
吴俐俊, 周伟国.
高炉冷却壁与炉气之间的换热特性[J]. 钢铁研究学报, 2006, 18(11): 7–10.
( Wu Li-jun, Zhou Wei-guo. Heat transfer characteristics between gas flow and stave[J]. Journal of Iron and Steel Research, 2006, 18(11): 7–10. DOI:10.3321/j.issn:1001-0963.2006.11.002 ) |
[12] |
Shah R K, London A L.
Laminar flow forced convection in ducts:a source book for compact heat exchanger analytical data[M]. Pittsburgh: Academic Press, 2014: 1-55.
|
[13] |
程进.
基于响应面法的几何非线性结构概率响应分析[J]. 同济大学学报(自然科学版), 2006, 34(9): 1147–1151.
( Cheng Jin. Probabilistic response analysis of geometrically nonlinear structures based on systematic response surface method[J]. Journal of Tongji University(Natural Science), 2006, 34(9): 1147–1151. DOI:10.3321/j.issn:0253-374X.2006.09.003 ) |