一种大单重、大断面特厚钢板辊式淬火过程的控制方法

摘要

一种大单重、大断面特厚钢板辊式淬火过程的控制方法，包括比热模型、热传导系数模型，温度场模型以及修正模型。其实现过程为，首先输入钢板参数，包括厚度、长度、含碳量；工艺规程，包括辊道速度，加速度等；实测参数，包括出炉温度、空冷后温度、返红后温度等。利用温度场模型并调用比热模型、热传导系数模型依次对冷却空冷段、水冷段、返红段进行计算，并分别通过修正模型对温度场进行修正。模拟结果为得到一组不同厚度处冷却曲线及冷速曲线。该控制方法的优点在于结合生产实际情况，得到符合实际的温降、冷速曲线，通过模型计算替代部分实际调试过程，可用于指导工艺规程的调整，节约生产成本，并能获得板形良好的产品。

说明书

技术领域

本发明属于冶金技术领域，涉及钢板辊式淬火过程控制方法。

背景技术

特厚钢板淬火后温度场分布，在很大程度上会影响钢板的物理性能和机械加工性能。传统的浸入式淬火，将钢板浸入淬火池或淬火槽中，受容器体积限制，冷却速度较慢，且分布不均。辊式淬火，钢板通过转动的辊道进入淬火机，高压射流喷射到钢板表面，相比于传统淬火方式，冷却速度大幅度增加，钢板在短时间内即可冷却至室温。但由于钢板厚度较大，传热较慢，导致心部与表面温度相差较大，淬火后性能分布不均。淬火过程计算，能够较好的改善这一问题。通过设置不同的参数，计算得到钢板的温度曲线、冷速曲线，直观的观察到冷却中钢板不同位置出的温度分布情况及冷速变化情况，计算结果可以很好的服务于生产，降低生产成本。

专利CN105445319A公开了一种测定钢板表面水冷换热系数的方法及装置，由红外热像仪采集钢板非水冷面图像信息，得到非水冷面的实际冷却曲线，有限元模拟水冷面到非水冷面的全过程，通过不断调整水冷换热系数，得到非水冷面的模拟冷却曲线。专利为CN102507636A公开了一种测定钢的快速冷却过程界面换热系数的方法，通过将热电偶点焊连接到温度采集模块，得到表面温度变化数据，利用热处理软件得到冷却过程中的界面换热系数，再以此模拟工件冷却过程的温度变化，并与实测结果进行对比，直至变化趋势吻合良好；现有技术中均未涉及与大单重、大断面特厚钢板辊式淬火过程有关的控制方法，更未提及利用修正模型对大单重、大断面特厚钢板辊式淬火过程中的换热系数进行修正的具体方法。

发明内容

本发明的目的是提供一种大单重、大断面特厚钢板辊式淬火过程的控制方法，该方法适用于特厚钢板在淬火过程中的模型控制，以数学模型为基础，结合生产实际情况，以曲线的形式直观的表现出钢板头部、尾部在钢板表面、四分之一厚度处以及心部的温度变化及冷速变化，可用于指导工艺规程的调整，节约生产成本，并能获得板形良好的产品。

本发明的具体技术方案为，一种大单重、大断面特厚钢板辊式淬火过程的控制方法，其特征在于，包括如下步骤：

步骤1，读取钢板信息及工艺参数，包括钢板的厚度、长度、含碳量、辊速、加速度；

步骤2，设置初始换热系数，利用温降计算模型，并调用比热模型、热传导系数模型，根据实测空冷段温降利用换热系数的修正模型对空冷段换热系数进行修正，进而得到空冷段温度场；

步骤3，确定水冷段换热系数，包括高压段水冷换热系数的确定和低压段水冷换热系数的确定；高压段水冷换热系数的确定，淬火机高压段的作用主要是将钢板的表面温度快速冷却到一个较低的温度，由于高压段的作用时间远低于低压段且实际生产中无法直接获得，故采用实验所得经验数据；

低压段水冷换热系数的确定，由于钢板出淬火机时只能直接获得表面温度，而在淬火过程中钢板表面仅仅在淬火开始时温度变化剧烈，此后，至出淬火机，其温度趋于稳定，接近对流介质温度，故不能直接用于计算对流换热系数。而钢板内部的温度无法直接获得，故无法通过水冷段对水冷换热系数直接进行修正计算。钢板在离开淬火机时，由于内部的温度仍高于表面温度，通过热传导作用，会将内部的温度传递到钢板的表面，因此，可采用淬火后空冷时的表面返红温度利用换热系数的修正模型对低压段水冷段的换热系数进行修正；具体方法为以出炉空冷后的温度场为水冷段的初始温度场，给定低压段初始换热系数，进行温度场计算，并以水冷后温度场作为返红段初始温度场计算返红后表面节点的温度值，与测量值进行比较，调用换热系数的修正模型对水冷换热系数进行修正，空冷换热系数不变；并重新利用温降计算模型，并调用比热模型、热传导系数模型进行水冷段、返红段温度场计算，至差值在允许误差范围内；

步骤4，得到符合实际情况的钢板淬火过程中不同位置的温降曲线，冷速曲线。

由于无法直接测得淬火过程中的对流换热系数，因此需要通过比较计算温度与实测温度，利用修正模型对换热系数进行修正，进而对温度场进行修正。上述步骤2、步骤3中所述换热系数的修正模型为：

给定初始换热系数区间范围[0,A]，取上限值A作为初始换热系数进行温度场计算，若计算值高于目标值，则取换热系数范围为[A,1.5A]；每次更改区间以上一区间的上限值作为新区间的下限值，并以新区间的下限值的1.5倍作为新区间的上限值，至换热系数取值为该区间的上限，计算值低于目标值为止；表明此时实际对流换热系数处于该区间内，在该区间内采用黄金分割法，不断缩小所在区间，至实测值与目标值差值在允许误差范围内，此时换热系数值即为实际值。

进一步地，上述步骤2、3中：

1)比热模型的计算；比热系数主要与钢板的含碳量及温度有关；含碳量以设定的定值作为界定范围，当含碳量不为上述值时，先确定其所对应的左右界值，通过插值的方式确定含碳量的权重，再比定温度所在区间，从而确定钢板的比热值；

2)热传导系数模型的计算；首先实验得出不同含碳量钢板在不同温度下的比热值及热传导系数值，然后通过插值的方式确定其他含碳量、其他温度对应的比热及热传导系数值；

3)温降计算模型为：

建立笛卡尔坐标系中的一维非稳态导热微分方程：

其中：

x为划分单元格长度；d为钢板厚度；t为时间；T为温度；a为导温系数，为内热源；λ为淬火板热传导系数；ρ为淬火钢板密度；c为淬火板比热；

计算时将钢板冷却过程中的相变潜热计入平均比热内，因此可以忽略内热源；

初始条件为：

T(x,0)＝T0(0＜x＜d,t＞0)

边界条件为：

为了提高傅里叶数的收敛性和稳定性并使模型更具有更小的误差，采用Crank-Nicolson差分方法；

t为时间；i为节点，0≤i≤I；建立温度场如下:

内部节点：

边界节点：

其中

h_x为对流换热系数；T_f为水温；为时间为t时钢板第i个节点所对应的温度值；F_ox为傅里叶数；B_ix为毕渥数；

稳定性条件为:

当已知初始温度场及换热系数的条件下，通过差分计算出之后任一时刻任一节点的温度分布情况。

进一步地，上述温降计算模型中冷却时间的控制：钢板的模型计算根据钢板所在的辊道位置分为三部分，分别为进入淬火机前的空冷段、通过淬火机时的淬火段和通过淬火机以后的返红阶段；

空冷段的时间确定；由于钢板具有一定的长度，不同位置进入淬火机所需要的时间不同，因此对头尾分别进行计算，钢板头部进入淬火机前为匀速运动，通过头部到淬火机的距离及初速度进行计算；当钢板头部进入淬火机以后辊道开始施加一定的加速度，所以钢板的尾部在空冷段即开始进行加速，此时及距离淬火机的位置为钢板长度，通过初速度、距离以及加速度来对尾部空冷时间进行计算；

淬火段时间的确定；淬火段时间分为通过高压段的时间和通过低压段的时间首先确定高压段长度，根据设定好辊道的初速度和加速度，直接计算出钢板头部的淬火时间；由于在钢板头部进入淬火机时即对钢板进行加速，通过之前空冷加速部分的时间和加速度确定尾部进入淬火机时的速度，再根据此时的速度和加速度计算确定其通过高压段所需要的时间，钢板通过低压段时间根据摆动时间确定；

返红段时间的确定；其具体方法为在钢板出淬火机时，通过秒表开始进行计时，针对钢板的同一位置，测量其在不同时刻的返红温度，返红结束后停止计时，取返红温度最大值作为模拟计算中的目标温度，其对应时间作为返红阶段的时间。

进一步地，上述步骤2、3中各段初始温度场模型的建立：以钢板出加热炉时温度作为空冷段的初始温度场，利用温度场模型计算得到的空冷后模拟温度场与淬火机前测温点测得的温度进行比较修正，最终得到符合实际的空冷后温度场，并以此作为水冷段的初始温度场，水冷段不需要进行修正计算，直接将模型计算结果传递到返红阶段，作为返红阶段的初始温度场。

进一步地，上述步骤4输出结果包括空冷段、水冷高低压段换热系数值，钢板头部、尾部的表面、四分之一、中心部的温度变化曲线及冷速变化曲线。

本发明具有以下优点：

1)初始参数符合实际生产规程，考虑到现场生产条件，所需测量数据在生产中可以直接获得；

2)各阶段计算结果分别依据实测值进行修正，计算结果与钢板实际冷却曲线更贴近；

3)输出结果为不同厚度处冷却曲线、冷速曲线，更加直观的展现心、表温差，心、表冷速差，头尾温差，头尾冷速差，可用于指导调整生产工艺；

4)通过计算代替部分现场调试，减少能源消耗，降低生产成本。

附图说明

图1：计算过程流程图。

图2：计算所得温度变化曲线。

图3：计算所得冷却速度变化曲线。

图4：过滤高压段后冷却速度变化曲线。

具体实施方式

1)参数的输入。包括钢板参数：厚度、长度、含碳量；工艺规程：辊道速度，加速度，淬火机高压段长度，出炉温度，换热系数初值；实测参数：出炉温度、进入淬火机前温度、返红时间、返红后温度；

2)空冷段、淬火段、返红段时间的确定。由于钢板具有一定的长度，不同位置进入淬火机所需要的时间不同，因此对头尾分别进行计算，钢板头部进入淬火机前为匀速运动，进入淬火机高压段后为加速运动，此后在低压段进行摆动。钢板尾部的运动方式为进入淬火机前空冷段首先为匀速运动，当钢板头部进入淬火机后为加速运动，高压段为加速运动，此后在低压段进行摆动。空冷段及高压段时间通过距离、初速度、加速度进行计算，低压段时间通过设定的摆动时间确定，返红段通过计时确定；

3)空冷段温度场的计算，以出炉温度建立初始温度场，每个时间步分别对不同节点的比热、热传导系数进行计算，然后调用温度场模型对空冷温度场进行计算，将计算结果与实测值进行比较，调用修正模型，对空冷换热系数进行修正，得到空冷后温度场。

4)水冷段时间的确定。分别计算钢板头部、中间部位、尾部通过淬火机高压段所需要的时间，根据摆动时间确定钢板通过低压段时间；

5)水冷温度场的计算。以计算的空冷后温度场为初始温度场，调用比热模型、导热系数模型、温度场模型，计算获得水冷后温度场；

6)返红后温度场的计算。以水冷后温度场为的初始温度场，换热系数值采用修正后的空冷换热系数值，调用比热模型、导热系数模型、温度场模型对返红温度场进行计算，将计算结果与实测值进行比较，未在允许误差范围内，调用修正模型对水冷换热系数进行修正，重新计算步骤5、6；

7)计算结果的输出。绘制钢板头尾部表面、四分之一厚度处、心部在空冷段及淬火段的温降曲线、冷速曲线。

实施例：

钢板厚度132mm，长度7250mm，含碳量0.15％，辊道速度0.2m/s，加速度0.00015m/s²，淬火机高压段长度3.2m，换热系数初值：空冷段100W/(m²K)，高压段20000W/(m²K)，低压段8000W/(m²K)，水温22.1℃、出炉温度910℃、进入淬火机前温度830℃、空冷时间45s，低压摆动时间1560s，返红时间142s，返红后温度28℃。取单元格长度1mm，时间步长0.5s，计算流程如附图1所示。

计算结果：空冷段换热系数值124.64W/(m²K)，低压段水冷换热系数2250W/(m²K)。温降曲线如附图2所示，由图可知，钢板经由空冷段进入水冷的高压段，表面温度迅速降低，进入低压段后由于内部的温度向外传递，温度有少许的回升，相较于表面，四分之一厚度处及心部温度冷却较缓慢。冷速曲线如附图3所示，与温降曲线相对应可以发现，在高压段，表面冷速迅速升高，进入低压段后，冷速曲线为负值代表温度有所上升，由于表面短时间内冷速远高于其他时刻，不便于观察其他时刻冷速，因此过滤掉部分后冷速如附图4所示，可以发现，四分之一和心部在心表温差较大时冷速逐渐增加，此后随着温度的降低，冷速逐渐降低。