CN109711113A - 一种大方坯凸形辊连铸重压下偏析行为预测方法 - Google Patents

Abstract

本发明属于钢铁冶金工业中的大方坯连铸生产领域，具体涉及一种大方坯凸形辊连铸重压下偏析行为预测方法,先准确描述大方坯连铸过程坯壳形貌，然后准确描述大方坯连铸过程坯壳表面速度，最后进行多相凝固模型耦合计算。本发明是高速度、高效率、高精度的实现连铸凝固末端凸形辊重压下实施过程中不同浇铸、压下工艺条件下的大方坯内部宏观偏析行为定量计算，阐明凸形辊重压下过程中压下变形对大方坯凝固偏析的改善机理，为实现高均质度铸坯的稳定生产提供可靠的理论支撑。

Description

一种大方坯凸形辊连铸重压下偏析行为预测方法

技术领域

本发明属于钢铁冶金工业中的大方坯连铸生产领域，具体涉及一种大方坯凸形辊连铸重压下偏析行为预测方法。

背景技术

连铸大方坯内部宏观偏析缺陷作为典型的内部质量问题，且该缺陷不能通过轧制和后处理等手段有效去除，严重的影响最终轧材的机械性能。目前，连铸大方坯凝固末端凸形辊重压下技术是替代传统轻压下工艺的最为有效的技术手段之一。相比较传统的轻压下技术，大方坯凝固末端重压下技术能够通过较大的机械压下变形，有效的促进凝固前沿富集溶质的反向流动与二次分布，达到改善铸坯宏观偏析行为的工艺效果。同时，相比平辊重压下工艺，大方坯凝固末端凸辊重压下技术能够避开大方坯两侧已凝坯壳，显著降低辊坯接触抗力，在液压系统未整体升级条件下提升单辊压下量2.5倍以上，极大的提高压下力向铸坯中心的渗透效果，且能够进一步的促进凝固前沿富集溶质的反向传输，更好的改善宏观偏析行为。

鉴于大方坯凝固末端凸辊重压下技术在改善铸坯宏观偏析方面存在明显技术优势，该技术已得到许多研究者的重视，并开展了相关的科研工作。

中国专利CN201710948669.7公开了一种用于低圧缩比轧制的连铸坯生产方法。其在凝固末端采用凸辊重压下技术轧制，通过调整铸机的浇铸工艺与重压下设备的压下量及压下分配，实现对铸坯中心质量的有效控制。但是，未考虑实施凝固末端凸辊重压下工艺后的偏析改善效果。

中国专利CN201611029498.X公开了一种连铸坯中心偏析行为研究的数值模拟方法。其对连铸坯流场进行计算区域划分，并通过界面数据转换实现多维度模型数据之间的链接，利用ANSYS与Fluent软件实现了连铸坯中心偏析行为的计算。但是，该模型并不能考虑凝固末端机械压下作用下的坯壳变形对铸坯偏析行为的影响。

韩国浦项公司通过在2#大方坯连铸机凝固末端施加凸形辊(ISIJInternational，Vol.52(2012)，No.7，pp.1266-1272)，使技术实施后的大方坯中心偏析评级从1.6降低为0.2，相比平辊压下显著的改善了大方坯内部的宏观偏析行为。但是该技术尚处于工业试验研究，并没有介绍凸形辊压下过程的大方坯偏析行为的定量预测。

目前，已有众多关于压下工艺及装备设计的相关技术专利，但凝固末端凸形辊压下过程铸坯宏观偏析行为的定量预测方法还鲜有报道。

发明内容

针对连铸大方坯凝固末端凸形辊压下过程宏观偏析行为定量预测方法上的研究不足，本发明公开了一种大方坯凸形辊连铸重压下偏析行为预测方法，目的是高速度、高效率、高精度的实现连铸凝固末端凸形辊重压下实施过程中不同浇铸、压下工艺条件下的大方坯内部宏观偏析行为定量计算，阐明凸形辊重压下过程中压下变形对大方坯凝固偏析的改善机理，为实现高均质度铸坯的稳定生产提供可靠的理论支撑。

具体技术方案如下：

一种大方坯凸形辊连铸重压下偏析行为预测方法，包括如下步骤：

步骤1：准确描述大方坯连铸过程坯壳形貌

基于大方坯连铸过程中凝固末端轻压下与凸形辊重压下工艺，忽略铸坯坯壳的z向展宽，对铸坯的坯壳形貌进行准确描述将连铸铸流划分为4个分区，分别为Zone 1结晶器与二冷段、Zone 2轻压下段、Zone 3凸形辊重压下段、Zone 4水平段，将铸流沿y方向同时被划分为5个子区域；

各个分区内的y方向高度计算公式如下：

y₁(x,z)＝T,x＝x₀～x₁,z＝z₀～z₅ (1)

y₄(x,z)＝y₃(x₄,z),x＝x₃～x₄,z＝z₀～z₅ (4)

式中，y₁(x,z)、y₂(x,z)、y₃(x,z)、y₄(x,z)为铸流各个分区内的y方向高度，mm；x₀、x₁、x₂、x₃、x₄分别为铸流各段表面形貌函数对应的x轴坐标，mm；z₀、z₁、z₂、z₃、z₄、z₅分别为铸流内各子区域表面形貌函数对应的z轴坐标，mm；T为大方坯的初始厚度，mm；r_a ^S和r_a ^H分别为凝固末端轻压下工艺压下量和凸形辊重压下工艺压下量，mm；

对于Zone 4水平段，在凸形辊重压下工艺实施后，该区域内内沿y方向5个子区域内表面形貌的具体计算公式如下：

式中，a₀、a₁、a₂、a₃、a₄分别为Zone 3凸形辊压下段内z₁～z₂侧弧线过渡区的高次多项式函数的辊型系数；b₀、b₁、b₂、b₃、b₄为Zone 3凸形辊压下段内z₄～z₅侧弧线过渡区的高次多项式函数的辊型系数；n为Zone 3凸形辊压下段内中心凸形辊凸起段对应的压下厚度，mm；

步骤2：准确描述大方坯连铸过程坯壳表面速度

由于铸坯坯壳z向展宽被忽略，所以Zone 1～4内z轴方向的速度全部为0；根据式(1)～式(5)所示的坯壳表面形貌公式，可以推出Zone 1～4内沿x、y方向铸坯内弧侧的坯壳表面速度公式如下：

式中，为Zone 1～4内沿x方向铸坯内弧侧的坯壳表面速度，m/min；为Zone 1～4内沿y方向铸坯内弧侧的坯壳表面速度，m/min；为了精准的表征大方坯凸辊重压下实施后铸坯内部的速度分布，其Zone 1～4内沿x、y方向的铸坯内部速度公式如下：

式中，为Zone 1～4内沿x方向的铸坯内部速度，m/min；为Zone 1～4内沿y方向的铸坯内部速度，m/min；y₂ ⁱⁿ(x,z)和y₃ ⁱⁿ(x,z)分别为Zone 2和Zone 3区域铸坯内任一节点的y方向数值，mm；

步骤3：多相凝固模型耦合计算过程

在准确描述铸坯坯壳形貌及速度分布的基础上，建立体积平均多相凝固耦合计算模型，对凝固过程中液相、柱状晶相、等轴晶相的溶质传输方程、质量传输方程、动量传输方程、热量传输方程进行耦合计算；上述方程的具体表述形式如下：

式中，i为各传输方程中的相参数，分别代表液相l、柱状晶相c和等轴晶相e；f_i为各相的体积分数，％；ρ_i为各相的密度，kg/m³；c_i为各相的溶质浓度，wt％；为各相速度，m/min；H_i为各相的焓，J/mol；C_s、M_s、D_s、H_s分别为溶质传输方程、质量传输方程、动量传输方程、热量传输方程的源项；

通过“相耦合Simple”算法，实现体积平均多相凝固模型的耦合计算；采用精度和收敛性高的耦合隐式算法，并控制收敛残差在10^-5以下，得到连铸凝固末端凸形辊重压下实施过程中不同浇铸、压下工艺条件对大方坯内部宏观偏析缺陷的改善作用规律。与现有技术相比，本发明具有如下有益技术效果：

本发明针对连铸大断面方坯，基于凝固末端凸形辊重压下后的坯壳形貌，提出了适用于连铸凝固末端凸形辊重压下工艺实施过程中，不同浇铸、压下工艺条件下的大方坯内部宏观偏析行为定量计算方法，有助于高速度、高效率、高精度的定量评估不同生产工艺对大方坯内部宏观偏析的改善效果，从而为浇铸、压下工艺及相关设备开发提供定量化的数据支撑。通过模型的计算，发现相比传统轻压下工艺，采用大方坯凝固末端凸形辊重压下工艺可以使大方坯中心偏析度由1.18降低到1.10。同时，通过实际检测铸坯断面的元素分布规律，验证了模拟计算结果的准确性，进一步定量说明了该模型有助于高速度、高效率、高精度的定量评估大方坯凝固末端的宏观偏析缺陷。此外，本发明也填补了凸形辊压下过程大方坯内部宏观偏析计算方法上的研究空白，丰富了连铸大方坯凝固末端重压下工艺理论体系。

附图说明

图1为连铸大方坯铸流计算区域划分示意图；

图2为连铸大方坯横向子区域划分示意图；

图3为连铸大方坯坯壳厚度变化趋势；

图4为连铸大方坯宏观偏析计算结果；

图5为平辊轻压下工艺生产的铸坯纵断面低倍照片；

图6为凸形辊重压下工艺生产的铸坯纵断面低倍照片；

图7为平辊轻压下与凸型辊重压下工艺下的偏析行为实测对比。

具体实施方式

下面结合附图和具体实施例对本发明进行详细说明，但本发明的保护范围不受实施例所限。

实施例1：大方坯压下过程宏观偏析行为预测：

结合现场的大方坯连铸机实际生产工艺，对凸形辊压下过程的宏观偏析进行模拟预测，主要分为以下步骤：

步骤1：大方坯连铸过程坯壳形貌的准确描述

图1为针对本例建立的大方坯凸形辊重压下过程三维几何模型。本例中，大方坯的厚度为280mm，宽度为380mm，几何模型中的大方坯采用6面体网格划分，网格尺寸为5mm²，网格数量为3,400,000。大方坯凝固末端重压下凸形辊安装在6#拉矫机，位于距结晶器27.95m的位置上。结合现场生产实际，附图1中模型对应的纵向区域划分参数如表1所示，附图2中模型对应的横向子区域划分参数如表2所示。

表1计算模型纵向区域划分

表2计算模型横向子区域划分

凝固末端轻压下工艺压下量r_a ^S为2～13mm，凸形辊重压下工艺压下量r_a ^H为9～15mm。在凸形辊横向辊身轮廓中，两侧过渡区内高次多项式函数的辊型系数如表3所示，n为30mm。

表3高次多项式函数的辊型系数

步骤2：大方坯连铸过程坯壳表面速度的准确描述

在步骤1得到的大方坯凸形辊重压下过程的准确坯壳表面形貌的基础上，通过公式(6)～(14)实现对大方坯连铸过程坯壳表面速度的准确描述。其中拉坯速度为0.72m/min。

步骤3：多相凝固模型的耦合计算

采用“相耦合Simple”算法及精度和收敛性高的耦合隐式算法，收敛残差为10^-6。

步骤4：大方坯凸形辊重压下宏观偏析预测结果

采用该模型进行连铸大方坯凝固末端凸形辊重压下偏析行为的计算，可以得出图3中0.72m/min拉速下对应的坯壳厚度变化趋势。从坯壳厚度计算结果可以得出，凝固终点位于距结晶器弯月面28.32m处，处于6#拉矫机凸形辊位置(27.95m)之后，表明该位置处的凸形辊重压下可以有效的促进大方坯内部中心富集溶质熔体发生反向流动，进而改善偏析。

此外，根据模型计算的大方坯宏观偏析模拟结果，如图4所示，能够得到沿y方向的偏析度分布趋势。模型预测结果可以看出，相比传统轻压下工艺，采用大方坯凝固末端凸形辊重压下工艺可以使大方坯中心偏析度由1.18降低到1.10。

将大方坯凝固末端渐变曲率凸型辊投入使用，根据现场试验的铸坯低倍结果对比，图5中平辊轻压下工艺生产的铸坯纵断面中心位置存在着明显的中心宏观偏析行为，而图6凸形辊重压下工艺生产的铸坯纵断面中心位置的宏观偏析行为则较轻。由此可见，实施凸型辊重压下后，铸坯中心偏析改善效果明显。

为了定量表征实际生产中凸形辊重压下工艺实施后的宏观偏析改善程度，对图5和图6中的铸坯进行了横向的钻屑实验，并利用碳硫分析仪对钻屑进行了成分检测，结果如图7所示。可以看出，大方坯凝固末端凸形辊重压下工艺实施后，相比传统的平辊轻压下工艺，铸坯的中心偏析度可以由1.17降低到1.11。将图4中的模拟预测结果与图7的实际生产结果进行对比，验证了模拟计算结果的准确性，且进一步定量说明了大方坯凝固末端凸形辊重压下工艺可以有效的改善铸坯内部的宏观偏析缺陷。

以上所述的实施例仅为本发明较佳的具体实施方式，本发明的保护范围不限于此，任何熟悉本领域的技术人员在本发明披露的技术范围内，可显而易见地得到该技术方案的简单变化或等效替换，并应用于连铸板坯凸形辊凝固末端压下工艺中，均属于本发明的保护范围之内。

Claims (1)

1.一种大方坯凸形辊连铸重压下偏析行为预测方法，其特征在于，包括如下步骤：

步骤1：准确描述大方坯连铸过程坯壳形貌

基于大方坯连铸过程中凝固末端轻压下与凸形辊重压下工艺，忽略铸坯坯壳的z向展宽，对铸坯的坯壳形貌进行准确描述将连铸铸流划分为4个分区，分别为Zone 1结晶器与二冷段、Zone 2轻压下段、Zone 3凸形辊重压下段、Zone 4水平段，将铸流沿y方向同时被划分为5个子区域；

各个分区内的y方向高度计算公式如下：

y₁(x,z)＝T,x＝x₀～x₁,z＝z₀～z₅ (1)

y₄(x,z)＝y₃(x₄,z),x＝x₃～x₄,z＝z₀～z₅ (4)

式中，y₁(x,z)、y₂(x,z)、y₃(x,z)、y₄(x,z)为铸流各个分区内的y方向高度，mm；x₀、x₁、x₂、x₃、x₄分别为铸流各段表面形貌函数对应的x轴坐标，mm；z₀、z₁、z₂、z₃、z₄、z₅分别为铸流内各子区域表面形貌函数对应的z轴坐标，mm；T为大方坯的初始厚度，mm；r_a ^S和r_a ^H分别为凝固末端轻压下工艺压下量和凸形辊重压下工艺压下量，mm；

对于Zone 4水平段，在凸形辊重压下工艺实施后，该区域内内沿y方向5个子区域内表面形貌的具体计算公式如下：

式中，a₀、a₁、a₂、a₃、a₄分别为Zone 3凸形辊压下段内z₁～z₂侧弧线过渡区的高次多项式函数的辊型系数；b₀、b₁、b₂、b₃、b₄为Zone 3凸形辊压下段内z₄～z₅侧弧线过渡区的高次多项式函数的辊型系数；n为Zone 3凸形辊压下段内中心凸形辊凸起段对应的压下厚度，mm；

步骤2：准确描述大方坯连铸过程坯壳表面速度

由于铸坯坯壳z向展宽被忽略，所以Zone 1～4内z轴方向的速度全部为0；根据式(1)～式(5)所示的坯壳表面形貌公式，可以推出Zone 1～4内沿x、y方向铸坯内弧侧的坯壳表面速度公式如下：

式中，为Zone 1～4内沿x方向铸坯内弧侧的坯壳表面速度，m/min；为Zone 1～4内沿y方向铸坯内弧侧的坯壳表面速度，m/min；为了精准的表征大方坯凸辊重压下实施后铸坯内部的速度分布，其Zone 1～4内沿x、y方向的铸坯内部速度公式如下：

式中，为Zone 1～4内沿x方向的铸坯内部速度，m/min；为Zone 1～4内沿y方向的铸坯内部速度，m/min；y₂ ⁱⁿ(x,z)和y₃ ⁱⁿ(x,z)分别为Zone 2和Zone 3区域铸坯内任一节点的y方向数值，mm；

步骤3：多相凝固模型耦合计算过程

在准确描述铸坯坯壳形貌及速度分布的基础上，建立体积平均多相凝固耦合计算模型，对凝固过程中液相、柱状晶相、等轴晶相的溶质传输方程、质量传输方程、动量传输方程、热量传输方程进行耦合计算；上述方程的具体表述形式如下：

式中，i为各传输方程中的相参数，分别代表液相l、柱状晶相c和等轴晶相e；f_i为各相的体积分数，％；ρ_i为各相的密度，kg/m³；c_i为各相的溶质浓度，wt％；为各相速度，m/min；H_i为各相的焓，J/mol；C_s、M_s、D_s、H_s分别为溶质传输方程、质量传输方程、动量传输方程、热量传输方程的源项；

通过“相耦合Simple”算法，实现体积平均多相凝固模型的耦合计算；采用精度和收敛性高的耦合隐式算法，并控制收敛残差在10^-5以下，得到连铸凝固末端凸形辊重压下实施过程中不同浇铸、压下工艺条件对大方坯内部宏观偏析缺陷的改善作用规律。