津西H型钢凝固过程温度场和应力场模拟研究.docVIP

津西H型钢凝固过程温度场和应力场模拟研究 杨景军1) 韦士来2) 1)北京科技大学冶金与生态工程学院，北京 100083 2) 河北津西钢铁股份有限公司，唐山 063006 摘 要 采用数值模拟的方法对津西钢铁公司炼钢厂H型钢连铸过程中的温度场、应力场进行模拟研究，分析H型钢腹板裂纹产生的主要原因，并对二次冷却区的冷却方案进行改进。在腹板表面中心区，矫直点温度正好处在Q235b的塑性最低点，同时该位置的热应变也很大，因此在矫直力作用下该位置容易出现裂纹；为避免腹板表面裂纹产生，可将二冷三段分成两个部分，前半部分维持原有水流密度，自三段后半部分至五段切断腹板表面喷水，让腹板自然冷却。 关键词 数值模拟；裂纹；二次冷却 1 引言 随着连铸技术的快速发展，铸坯质量始终是不可忽视的重要方面。众所周知连铸凝固过程的传热行为对连铸坯质量保证有着极其重要的作用。根据不同的铸机特点，铸坯坯型均要采取特定的冷却制度保证连铸坯的无缺陷化和连铸生产高效顺行[1]。 津西钢铁公司是我国大规格H型钢的主要生产厂家。自从2006年5月投产以来，有时出现铸坯腹板裂纹，影响了铸坯的合格率，考虑到有关这方面的研究较少，为此本文采用数值模拟的方法对河北津西钢铁股份有限公司炼钢厂异型坯凝固过程进行模拟，深入分析了连铸过程中腹板裂纹形成的原因，优化津西H型钢的连铸生产中的相关工艺参数。这方面的研究工作对津西钢厂H型钢质量改进具有较大的现实意义。 2 数学模型的建立 2.1 裂纹的特征 H型钢产生的腹板裂纹如下图所示： 2.2 温度场数学模型 热过程和凝固过程模拟的主要任务是对给定过程参数和设备参数下连铸温度场进行模拟，预测铸坯温度场和凝固坯壳厚度，同时发现连铸过程在 传热中存在的问题。该过程二维非稳态热传输方程为： …..（1） 在结晶器上表面，假设钢水温度与中间包浇注温度（）相等，也即时，。 结晶器内边界条件属于第二类边界条件，表达式为q=a－b,式中a为经验常数，b是结晶器冷却水量、冷却水温差、铸坯结构尺寸等参数的函数。 在二冷区，边界结点热流通量与边界温度成线性关系，即,式中为铸坯表面温度；为冷却水温度；为铸坯与水之间的给热系数，这里取，其中w为水速，该值可由冷却水量、冷却水温等参数计算得到。 在空冷区，以辐射传热为主，辐射传热热流密度采用四次方定律。 上述连铸过程数学模型基于下列假设：铸坯两相区和液芯对流传热靠增大钢液导热系数4倍进行考虑；,,仅随温度变化；忽略结晶器振动、结晶器锥度和凝固偏析的影响；忽略弯月面的影响；视弯月面处熔融金属的温度为浇注温度。 2.3 应力场数学模型 连铸应力过程数学模型可以预测铸坯坯壳内的应力和应变分布，为所生产的连铸坯质量，尤其是铸坯裂纹的产生提供非常重要的数据。本文采用二维瞬态热弹塑性方程进行模拟计算。计算过程中考虑了钢水静压力对凝固坯壳的影响；采用von Mises屈服准则来计算钢的应变强化。总应变方程为： （2） 其中, , 分别代表弹性应变量、塑性应变量、热应变量[2]。 3 模拟结果 根据生产现场连铸过程数据（表1）对H型钢连铸温度场和应力场进行了模拟。 表1 输入参数 单位 数据值 铸坯类型 BB1 金属 Q235b 结晶器有效长度 mm 700 结晶器宽面水量消耗 l/s 30 结晶器窄面水量消耗 l/s 16.7 二冷一段长度 mm 660 二冷一段冷却水量 l/s 4.22 二冷二段长度 mm 1500 二冷二段冷却水量 l/s 3.13 二冷三段长度 mm 2400 二冷三段冷却水量 l/s 1.90 二冷四段长度 mm 2400 二冷四段冷却水量 l/s 0.78 二冷五段长度 mm 2400 二冷五段冷却水量 l/s 0.35 矫直点距结晶器液面距离 m 19.14 拉速 m/min 0.98 浇注温度 ℃ 1549 液相线温度 ℃ 1519 固相线温度 ℃ 1464 3.1 温度场模拟结果及验证 为了验证铸坯传热计算过程的合理性，在模拟过程中先取H型钢的几个特殊点（图1）对计算结果进行验证。图2列出了计算得到特殊点的温度历程图，表2对计算温度和实测温度进行了比较。由表2可以看出热模拟计算温度与实测温度比较吻合。图3-6分别列出了几个位置的温度等值线图。 图1 H型钢铸坯横断面计算网格划分及特殊位置点 图2 铸坯上特殊点位置的温度历程 2 计算温度与实测温度比较 测温位置 铸坯测温点 计算温度(℃) 实测温度(℃) 绝对误差 相对误差(%) 二冷出口处 (距离结晶器液面10.06米,对应图1中616秒) 腹板表面 851 829 21 2.53 翼缘表面