CN110773573B - 一种板形调控功效系数实测数据处理方法 - Google Patents
一种板形调控功效系数实测数据处理方法 Download PDFInfo
- Publication number
- CN110773573B CN110773573B CN201911077070.6A CN201911077070A CN110773573B CN 110773573 B CN110773573 B CN 110773573B CN 201911077070 A CN201911077070 A CN 201911077070A CN 110773573 B CN110773573 B CN 110773573B
- Authority
- CN
- China
- Prior art keywords
- plate shape
- data
- regulation
- coefficient
- fitting
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 13
- 238000003672 processing method Methods 0.000 title description 3
- 238000004364 calculation method Methods 0.000 claims abstract description 59
- 238000005096 rolling process Methods 0.000 claims abstract description 26
- 230000000875 corresponding Effects 0.000 claims description 40
- 238000005452 bending Methods 0.000 claims description 22
- 230000001105 regulatory Effects 0.000 claims description 21
- 230000001276 controlling effect Effects 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000005315 distribution function Methods 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B37/00—Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
- B21B37/28—Control of flatness or profile during rolling of strip, sheets or plates
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Abstract
本发明公开了一种板形调控功效系数实测数据处理方法，属于冶金轧制技术领域，包括以下步骤：步骤S1、计算板形调控功效系数实测值，步骤S2、对实测板形调控功效系数进行排序，步骤S3、排序后对每组实测板形调控功效系数进行权重分配，计算加权板形调控功效系数，步骤S4、加权板形调控功效系数计算结果数据拟合，步骤S5、判断拟合板形调控功效系数是否满足要求。本发明使用变加权计算方法来对实测数据进行处理，从而获得准确的板形调控功效系数，使板形调控功效系数更接近实际情况，可以提高板形自动控制系统的控制能力，提高出口带材质量。
Description
技术领域
本发明涉及一种板形调控功效系数实测数据处理方法，属于冶金轧制技术领域。
背景技术
冷轧带材板形自动控制系统是保证成品带材质量的关键部分，该系统的是基于板形调控功效系数建立的。因此，只有获得准确的板形调控功效系数，才能够准确计算出用于消除板形偏差所需要的弯辊、横移和倾斜等调控机构的调节量。板形调控功效系数从实测板形应力分布的角度进行相关的分析，可以描述任意形态的板形调节性能。与传统模型相比，能够实现对板形测量信息的全面利用，有利于轧机板形控制能力的充分发挥和板形控制精度的提高。板形调控功效是一种板形控制技术的单位调节量作用下，轧机承载辊缝形状沿带钢宽度上各处的变化量。目前板形调控功效系数求解方法是针对特定的轧机和板形调控机构通过轧辊的弹性变形计算和轧件的塑性变形计算，或者通过有限元方法进行离线计算来完成。这两类方法虽然各有特点，但也都存在不足。采用数学解析方式建立的机理模型虽然可以进行在线计算，但在分析轧制过程时做了大量的假设，其计算结果往往与实际有较大的误差。采用有限元仿真模型建立的板形调控功效系数计算方法精度较高，但其计算效率又无法满足在线实时控制要求，只对某几种工况离线计算的方式无法全面反映轧制过程的工况，造成某些轧制条件下无法对板形缺陷进行有效控制。
而基于实测数据的板形调控功效系数自学习模型的计算虽然满足在线实时控制的要求，但其计算结果受到现场复杂工作环境的影响，大量的实测数据存在噪声、离群、遗漏和偏差等，使得计算结果不够理想，实测板形调控功效系数无法与实际情况相匹配。
发明内容
本发明需要解决的技术问题是提供一种板形调控功效系数实测数据处理方法，使用变加权计算方法来对实测数据进行处理，从而获得准确的板形调控功效系数，使板形调控功效系数更接近实际情况，可以提高板形自动控制系统的控制能力，提高出口带材质量。
为解决上述技术问题，本发明所采用的技术方案是：
一种板形调控功效系数实测数据处理方法，所述方法包括以下步骤：
步骤S1、计算板形调控功效系数实测值，具体包括以下步骤：
S11、沿板形宽度方向上设定m个通道，对应板形宽度上设定的m个测量段数，单位时间间隔需要采集的数据包括与m个通道相对应的m个实测板形变化量和与其对应的调控机构的变量，
S12、根据调控功效系数公式计算相应调控机构的实测板形调控功效系数，调控功效系数公式为：
其中Eff为相应调控机构的实测板形调控功效系数，ΔY为实测板形变化量，ΔU为对应的调控机构改变量，调控机构改变量包括轧辊倾斜、轧制力波动、工作辊弯辊、中间辊弯辊和中间辊横移；
步骤S2、对实测板形调控功效系数进行排序，具体包括以下步骤：
S21、将得到的单个实测板形调控功效系数做成表格，表格的每一列都表示数据的通道数，总计m列通道数据，表格的每一行都是单位时间间隔内采集到的实际数据计算得到一组实测板形调控功效系数数据，共有n组数据，
S22、将对应每个通道计算得到的n个实测板形调控功效系数使用排序算法从大到小进行顺序排列，使m个通道内的n组实测板形调控功效系数偏差较大的数据均集中在每个通道数据集的两端；
步骤S3、排序后对每组实测板形调控功效系数进行权重分配，计算加权板形调控功效系数，具体包括以下步骤：
S31、使用正态分布函数分别为排序后的n组实测板形调控功效系数进行权重分配，为了防止偏差较大的数据对结果产生影响，权重的分配规则符合正态分布排列，且每个通道数据集两端偏差较大的部分所分配的权重为0，权重为0的数据个数占数据集总个数的比例为20％，
S32、计算权重因子Rj，
S33、板形调控功效系数加权结果的计算，计算公式为：
式中：n表示实测数据的组数，Effj表示第j组板形调控功效系数，Rj表示第j组数据所对应的权重因子，Effw表示加权板形调控功效系数；
步骤S4、加权板形调控功效系数计算结果数据拟合，包括以下步骤：S41、采用最小二乘法的多项式对加权板形调控功效系数计算结果进行拟合，最小二乘法的多项式公式为：
式中：m表示沿板形宽度方向设定的通道数目，Effi表示第i个通道加权处理后的加权板形调控功效系数，E(i)表示第i个通道板形调控功效系数的数据拟合计算结果，
S42、拟合板形调控功效系数多项式拟合公式为：
E(i)＝a0+a1·i1+a2·i2+a3·i3+a4·i4
式中：E(i)表示第i个通道板形调控功效系数的数据拟合计算结果，a0…a4分别为拟合后的多项式系数；
步骤S5、判断拟合板形调控功效系数是否满足要求，
对比拟合计算得到的第i个通道板形调控功效系数的数据拟合计算结果E(i)与有限元仿真计算得到的理论板形调控功效系数Efi的相似程度，公式如下所示，
式中：m表示沿板宽方向设定的通道数目,E(i)表示第i个通道板形调控功效系数的数据拟合计算结果；Efi表示第i个通道有限元仿真计算得到的理论板形调控功效系数，ε为设定的数值，用来保证输出数据的精度，
对比结果满足设定的条件时，将计算得到的结果保存并用于板形自动控制系统；如果不满足公式的条件则需要重复步骤S1至S5。
本发明技术方案的进一步改进在于：步骤S1中的单位时间间隔取值为0.5s-1s。
本发明技术方案的进一步改进在于：步骤S32中权重因子的计算步骤为：
S321、将排序后m个通道对应的n组数据进行归一化处理，即1到n组数据编号被等距划分到[-1,1]的区间内，并用xj表示，划分公式为：
式中：xj表示第j组数据划分的等距点对应的数值，xj∈[-1,1]，j∈[1,n]，
S322、将xj带入下式，来确定第j组数据对应的概率因子；
式中：rj表示第j组数据对应的概率因子，
S323、对所有的概率因子rj进行求和，得到rz，并通过以下公式确定第j组数据对应的权重因子Rj，
式中：Rj表示第j组数据所对应的权重因子。
本发明技术方案的进一步改进在于：步骤S42中针对不同的调控机构，要采取不同的拟合策略，采用不同的多项式拟合板形调控功效系数，对于轧辊倾斜的拟合板形调控功效系数采用一次多项式拟合即可，即a2,a3和a4此时都为0，中间辊弯辊的拟合板形调控功效系数采用二次多项式拟合，即a3和a4为0，轧制力波动、工作辊弯辊和中间辊横移的拟合板形调控功效系数采用四次多项式拟合，即a0、a1、a2、a3和a4均不为0。
由于采用了上述技术方案，本发明取得的技术进步是：
本发明的一种板形调控功效系数实测数据处理方法，使用变加权计算方法来对实测数据进行处理，从而获得准确的板形调控功效系数，使板形调控功效系数更接近实际情况，可以提高板形自动控制系统的控制能力，提高出口带材质量。
将排序算法和中心极限定理引入数据处理过程，实现了实测数据的变加权处理，得到了与实际情况相接近的板形调控功效系数，可以进行实时在线计算，并将其用于冷轧板形自动控制系统，可以获得较高质量的冷轧带材。
附图说明
图1是本发明板形调控功效系数理论计算值折线图；
图2是本发明工作辊弯辊实测板形调控功效系数数据折线图；
图3是本发明工作辊弯辊排序后实测板形调控功效系数数据折线图；
图4是本发明理论板形调控功效系数、加权板形调控功效系数和拟合板形调控功效系数对比图；
图5是本发明实际轧制过程中的板形云图；
图6是本发明流程图；
图7是本发明实际轧制过程中板形残余均方差图。
具体实施方式
下面结合实施例对本发明做进一步详细说明：
本实施例公开了板形辊长度为1450mm五机架冷连轧机实测板形调控功效系数数据的处理过程。板形调控机构有轧辊倾斜、轧制力波动、工作辊弯辊、中间辊弯辊和中间辊横移，其中图1为四种板形调控机构和轧制力波动的板形调控功效系数有限元仿真计算得到的理论计算值。五机架冷连轧机的主要控制参数及轧制参数如下表1所示。
表1轧制过程主要参数
本实施例一种板形调控功效系数实测数据处理方法，包括以下步骤：
步骤S1、计算板形调控功效系数实测值，具体包括以下步骤：
S11、沿板形宽度方向上设定20个通道，对应板形宽度上设定的20个测量段数，单位时间间隔需要采集的数据包括与20个通道相对应的20个实测板形变化量和与其对应的调控机构的变量，本实施例中调控机构为工作辊弯辊。在采集数据时要保证板形变化量和调控机构的改变量是同步的，且数据采样的单位时间间隔不宜过小，一般取值为0.5s-1s，本实施例中采集数据的单位时间间隔为0.5s；
S12、根据调控功效系数公式计算工作辊弯辊的实测板形调控功效系数，板形调控功效是一种板形控制技术的单位调节量作用下，轧机承载辊缝形状沿带钢宽度上各处的变化量，在实际的应用中，一般沿板形宽度上设定20个测量段数，即为通道数CH为20，板形调节机构的数目为4个，同时将轧制力波动对板形的影响也纳入板形调控机构考虑，则板形调控功效系数可以表示为一个20×5的矩阵形式：
本实施例中，Eff为工作辊弯辊的实测板形调控功效系数，ΔY为实测板形变化量，ΔU为工作辊弯辊改变量；根据以上调节公式对板形测量值和工作辊弯辊改变量处理得到一系列单个工作弯辊的实测板形调控功效系数。
步骤S2、对实测板形调控功效系数进行排序，具体包括以下步骤：
S21、将得到的单个实测板形调控功效系数做成表格，表格的每一列都表示数据的通道数，从CH1到CH20，总计m列通道数据，表格的每一行都是单位时间间隔内采集到的实际数据计算得到一组实测板形调控功效系数数据，从GR1到GRn，共有n组数据，本实施例中带钢规格为宽度1000mm、厚度0.285mm，连续提取出了500组实测数据计算工作辊弯辊的实测板形调控功效数据，表1列出了部分直接计算得到的实测板形调控功效系数数据，共包括1,3,5-19奇数通道的10组数据，具体数值如表1实测板形调控功效系数数据，同时在计算完成之后的数据中，任意挑选出四组比较理想的数据绘制成折线图，如图2所示，从折线图可以看出直接通过实测数据计算出来的工作辊弯辊板形调控功效系数毫无规律、杂乱无章，不可以直接使用，还需要做进一步的处理。
表1实测板形调控功效系数数据
| GR | CH<sub>1</sub> | CH<sub>3</sub> | CH<sub>5</sub> | CH<sub>7</sub> | CH<sub>9</sub> | CH<sub>11</sub> | CH<sub>13</sub> | CH<sub>15</sub> | CH<sub>17</sub> | CH<sub>19</sub> |
| 1 | 32.40 | 1.97 | 6.04 | -1.81 | 16.02 | -10.26 | 3.70 | -7.82 | -10.21 | -17.62 |
| 2 | 7.40 | 4.97 | 1.29 | 1.54 | -2.44 | 0.91 | -0.78 | -1.58 | -3.78 | -3.53 |
| 3 | 0.10 | -2.85 | -1.73 | -3.13 | 0.03 | -1.96 | -0.13 | 2.61 | 2.35 | 5.28 |
| 4 | -0.14 | -2.35 | -2.77 | -3.02 | -1.76 | -1.09 | 1.09 | 1.92 | 3.21 | 5.28 |
| 5 | -40.07 | -14.99 | -44.67 | -20.48 | -24.59 | 11.45 | 29.82 | -1.82 | 32.21 | 39.01 |
| 6 | 1.76 | -3.28 | -1.24 | -4.37 | -0.27 | 0.19 | 0.72 | 0.17 | 1.36 | 3.84 |
| 7 | 5.41 | -2.94 | 2.27 | -3.98 | 2.28 | 1.24 | -1.07 | -2.36 | -2.11 | 3.64 |
| 8 | 0.53 | -0.72 | -1.70 | -2.89 | -2.36 | -1.87 | -0.66 | 3.76 | 4.77 | 1.03 |
| 9 | 1.31 | -0.06 | -2.35 | -3.58 | -2.35 | -0.72 | -1.73 | 2.74 | 3.83 | 1.93 |
| 10 | 1.96 | 0.39 | 0.17 | 0.98 | -0.87 | -1.30 | -1.26 | 0.01 | -1.65 | 1.09 |
S22、将对应每个通道计算得到的n个实测板形调控功效系数使用排序算法从大到小进行顺序排列，使m个通道内的n组实测板形调控功效系数偏差较大的数据均集中在每个通道数据集的两端，即对每一个通道CH内的多组GR1到GRn数据使用排序算法，使其按照从大到小的顺序排列，排序的目的是为了使偏差较大的数据聚集在数据集的两端，方便剔除或者分配较小的权重，减少其对最终结果的影响，能够使结果更接近实际情况。
将500组数据排序处理后，选出四组比较理想的实测调控功效系数数据绘制成如图3所示的排序后的实测调控功效系数数据折线图，从图3中可以看出，曲线形状虽然得到改善，但是还不够理想，同时表2也列出了表1实测板形调控功效系数数据经过排序后得到的结果，如表2所示。
表2排序后的实测板形调控功效系数数据
| GR | CH<sub>1</sub> | CH<sub>3</sub> | CH<sub>5</sub> | CH<sub>7</sub> | CH<sub>9</sub> | CH<sub>11</sub> | CH<sub>13</sub> | CH<sub>15</sub> | CH<sub>17</sub> | CH<sub>19</sub> |
| 1 | 32.40 | 4.97 | 6.04 | 1.54 | 16.02 | 11.45 | 29.82 | 3.76 | 32.21 | 39.01 |
| 2 | 7.40 | 1.97 | 2.27 | 0.98 | 2.28 | 1.24 | 3.70 | 2.74 | 4.77 | 5.28 |
| 3 | 5.41 | 0.39 | 1.29 | -1.81 | 0.03 | 0.91 | 1.09 | 2.61 | 3.83 | 5.28 |
| 4 | 1.96 | -0.06 | 0.17 | -2.89 | -0.27 | 0.19 | 0.72 | 1.92 | 3.21 | 3.84 |
| 5 | 1.76 | -0.72 | -1.24 | -3.02 | -0.87 | -0.72 | -0.13 | 0.17 | 2.35 | 3.64 |
| 6 | 1.31 | -2.35 | -1.70 | -3.13 | -1.76 | -1.09 | -0.66 | 0.01 | 1.36 | 1.93 |
| 7 | 0.53 | -2.85 | -1.73 | -3.58 | -2.35 | -1.30 | -0.78 | -1.58 | -1.65 | 1.09 |
| 8 | 0.10 | -2.94 | -2.35 | -3.98 | -2.36 | -1.87 | -1.07 | -1.82 | -2.11 | 1.03 |
| 9 | -0.14 | -3.28 | -2.77 | -4.37 | -2.44 | -1.96 | -1.26 | -2.36 | -3.78 | -3.53 |
| 10 | -40.07 | -14.99 | -44.67 | -20.48 | -24.59 | -10.26 | -1.73 | -7.82 | -10.21 | -17.62 |
步骤S3、排序后对每组实测板形调控功效系数进行权重分配，计算加权板形调控功效系数，具体包括以下步骤：
从表2以及图3都可以看出，实测板形调控功效系数存在波动大，数据不稳定等情况，因此，对于实测板形调控功效系数的处理引入中心极限定理对数据进行进一步的处理，中心极限定理是概率论中最重要的一类定理，被广泛应用于数理统计学和误差分析学领域，其指出大量独立同分布随机变量近似服从正态分布。因此，将实测的板形调控功效系数数据看成是独立同分布的随机变量，并使用正态分布函数来对每组数据进行权重分配，来实现变加权计算。
S31、使用正态分布函数分别为排序后的n组实测板形调控功效系数进行权重分配，使权重的分配规则符合正态分布排列，按照排列分布，两端所占权重小，中间所占权重大，较准确的数据分配较大的权重，差异偏大的数据分配较小的权重，为了防止偏差较大的数据对结果产生影响，权重的分配规则符合正态分布排列，且每个通道数据集两端偏差较大的部分所分配的权重为0，权重为0的数据个数占数据集总个数的比例为20％。
S32、计算权重因子Rj，
权重因子的计算步骤为：
S321、将排序后m个通道对应的n组数据进行归一化处理，即1到n组数据编号被等距划分到[-1,1]的区间内，并用xj表示，划分公式为：
式中：xj表示第j组数据划分的等距点对应的数值，xj∈[-1,1]，j∈[1,n]，
S322、将xj带入下式，来确定第j组数据对应的概率因子；
式中：rj表示第j组数据对应的概率因子，
S323、对所有的概率因子rj进行求和，得到rz，并通过以下公式确定第j组数据对应的权重因子Rj，
式中：Rj表示第j组数据所对应的权重因子。
表3中的数据为表2中10组数据分别对应的权重因子，如下表3所示：
表3实测板形调控功效系数权重因子
| n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
| R<sub>i</sub> | 0 | 0.006 | 0.043 | 0.155 | 0.294 | 0.294 | 0.155 | 0.043 | 0.006 | 0 |
S33、在权重因子确定后，板形调控功效系数加权结果的计算，计算公式为：
式中：n表示实测数据的组数，Effj表示第j组板形调控功效系数，Rj表示第j组数据所对应的权重因子，Effw表示加权板形调控功效系数。
步骤S4、加权板形调控功效系数计算结果数据拟合，包括以下步骤：
对于加权得到的调控功效系数曲线，还可能存在奇异点，需要进行曲线拟合获得更理想的板形调控功效系数曲线。
S41、采用最小二乘法的多项式对加权板形调控功效系数计算结果进行拟合，最小二乘法的多项式公式为：
式中：m表示沿板形宽度方向设定的通道数目，Effi表示第i个通道加权处理后的加权板形调控功效系数，E(i)表示第i个通道板形调控功效系数的数据拟合计算结果，
S42、拟合板形调控功效系数多项式拟合公式为：
E(i)＝a0+a1·i1+a2·i2+a3·i3+a4·i4
式中：E(i)表示第i个通道板形调控功效系数的数据拟合计算结果，a0…a4分别为拟合后的多项式系数；针对不同的调控机构，要采取不同的拟合策略，采用不同的多项式拟合板形调控功效系数，对于轧辊倾斜的拟合板形调控功效系数采用一次多项式拟合即可，即a2,a3和a4此时都为0，中间辊弯辊的拟合板形调控功效系数采用二次多项式拟合，即a3和a4为0，轧制力波动、工作辊弯辊和中间辊横移的拟合板形调控功效系数采用四次多项式拟合，即a0、a1、a2、a3和a4均不为0。
将加权计算的结果，拟合后的结果与理论计算的结果绘制在图4中，方块点代表的加权板形调控功效系数点线图已经大致符合板形调控功效系数曲线，为了去掉奇异点，对于工作辊弯辊加权计算后的结果使用多项式进行拟合，即将方块点的点线图经过四次多项式拟合成图4中所示的圆点组成的拟合值板形调控功效系数点线图，图中三角点组成的点线图为有限元仿真计算得到的理论计算的板形调控功效系数曲线，通过对比可以发现，轧机操作侧的板形调控功效系数曲线比较接近，而传动侧的板形调控功效系数则存在一定的差别。
步骤S5、判断拟合板形调控功效系数是否满足要求，
对比拟合计算得到的第i个通道板形调控功效系数的数据拟合计算结果E(i)与有限元仿真计算得到的理论板形调控功效系数Efi的相似程度，公式如下所示，
式中：m表示沿板宽方向设定的通道数目,E(i)表示第i个通道板形调控功效系数的数据拟合计算结果；Efi表示第i个通道有限元仿真计算得到的理论板形调控功效系数，ε为设定的数值，用来保证输出数据的精度，
对比结果满足设定的条件时，将计算得到的结果保存并用于板形自动控制系统，如果不满足公式的条件则需要重复步骤S1至S5。
通过板形自动控制系统绘制的板形云图，如图5所示，可以看到，操作侧的板形偏差都控制在都控制在5IU-10IU内，而传动侧的板形偏差波动较大，数值在-15IU-5IU内波动；可见，传动侧的板形偏差控制效果不理想。造成这种现象的原因是因为目前板形自动控制系统所采用的板形调控功效系数曲线为图4中的三角点画线，即有限元仿真计算得到的理论计算的板形调控功效系数曲线，而轧机实际的板形调控功效系数为图4中圆点所示的曲线，因此对于传动侧的板形偏差控制效果不好的是因为实际的板形调控功效系数与在使用中调控功效不匹配。而为了和轧机的实际工作情况相匹配，将本实施例建立的实测板形调控功效系数模型应用到冷轧带材板形自动控制系统，图7为本实施例建立的实测板形调控功效系数模型投入自动控制系统后的板形残余偏差图，从图中可以明显看到新的板形调控功效系数投入使用后，板形残余偏差明显下降。
Claims (4)
1.一种板形调控功效系数实测数据处理方法，其特征在于：所述方法包括以下步骤：
步骤S1、计算板形调控功效系数实测值，具体包括以下步骤：
S11、沿板形宽度方向上设定m个通道，对应板形宽度上设定的m个测量段数，单位时间间隔需要采集的数据包括与m个通道相对应的m个实测板形变化量和与其对应的调控机构改变量，
S12、根据调控功效系数公式计算相应调控机构的实测板形调控功效系数，调控功效系数公式为：
其中Eff为相应调控机构的实测板形调控功效系数，ΔY为实测板形变化量，ΔU为对应的调控机构改变量，调控机构改变量包括轧辊倾斜、轧制力波动、工作辊弯辊、中间辊弯辊和中间辊横移；
步骤S2、对实测板形调控功效系数进行排序，具体包括以下步骤：
S21、将得到的单个实测板形调控功效系数做成表格，表格的每一列都表示数据的通道数，总计m列通道数据，表格的每一行都是单位时间间隔内采集到的实际数据计算得到一组实测板形调控功效系数数据，共有n组数据，
S22、将对应每个通道计算得到的n个实测板形调控功效系数使用排序算法从大到小进行顺序排列，使m个通道内的n组实测板形调控功效系数偏差较大的数据均集中在每个通道数据集的两端；
步骤S3、排序后对每组实测板形调控功效系数进行权重分配，计算加权板形调控功效系数，具体包括以下步骤：
S31、使用正态分布函数分别为排序后的n组实测板形调控功效系数进行权重分配，为了防止偏差较大的数据对结果产生影响，权重的分配规则符合正态分布排列，且每个通道数据集两端偏差较大的部分所分配的权重为0，权重为0的数据个数占数据集总个数的比例为20％，
S32、计算权重因子Rj，
S33、板形调控功效系数加权结果的计算，计算公式为：
式中：n表示实测数据的组数，Effj表示第j组板形调控功效系数，Rj表示第j组数据所对应的权重因子，Effw表示加权板形调控功效系数；
步骤S4、加权板形调控功效系数计算结果数据拟合，包括以下步骤：
S41、采用最小二乘法的多项式对加权板形调控功效系数计算结果进行拟合，最小二乘法的多项式公式为：
式中：m表示沿板形宽度方向设定的通道数目，Effi表示第i个通道加权处理后的加权板形调控功效系数，E(i)表示第i个通道板形调控功效系数的数据拟合计算结果，
S42、拟合板形调控功效系数多项式拟合公式为：
E(i)＝a0+a1·i1+a2·i2+a3·i3+a4·i4
式中：E(i)表示第i个通道板形调控功效系数的数据拟合计算结果，a0…a4分别为拟合后的多项式系数；
步骤S5、判断拟合板形调控功效系数是否满足要求，
对比拟合计算得到的第i个通道板形调控功效系数的数据拟合计算结果E(i与有限元仿真计算得到的理论板形调控功效系数Efi的相似程度，公式如下所示，
式中：m表示沿板宽方向设定的通道数目,E(i)表示第i个通道板形调控功效系数的数据拟合计算结果；Efi表示第i个通道有限元仿真计算得到的理论板形调控功效系数，ε为设定的数值，用来保证输出数据的精度，
对比结果满足设定的条件时，将计算得到的结果保存并用于板形自动控制系统；如果不满足公式的条件则需要重复步骤S1至S5。
2.根据权利要求1所述的一种板形调控功效系数实测数据处理方法，其特征在于：步骤S1中的单位时间间隔取值范围为0.5s-1s。
3.根据权利要求1所述的一种板形调控功效系数实测数据处理方法，其特征在于：步骤S32中权重因子的计算步骤为：
S321、将排序后m个通道对应的n组数据进行归一化处理，即1到n组数据编号被等距划分到[-1,1]的区间内，并用xj表示，划分公式为：
式中：xj表示第j组数据划分的等距点对应的数值，xj∈[-1,1]，j∈[1,n]，
S322、将xj带入下式，来确定第j组数据对应的概率因子；
式中：rj表示第j组数据对应的概率因子，
S323、对所有的概率因子rj进行求和，得到rz，并通过以下公式确定第j组数据对应的权重因子Rj，
式中：Rj表示第j组数据所对应的权重因子。
4.根据权利要求1所述的一种板形调控功效系数实测数据处理方法，其特征在于：步骤S42中针对不同的调控机构，要采取不同的拟合策略，采用不同的多项式拟合板形调控功效系数，对于轧辊倾斜的拟合板形调控功效系数采用一次多项式拟合即可，即a2,a3和a4此时都为0，中间辊弯辊的拟合板形调控功效系数采用二次多项式拟合，即a3和a4为0，轧制力波动、工作辊弯辊和中间辊横移的拟合板形调控功效系数采用四次多项式拟合，即a0、a1、a2、a3和a4均不为0。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911077070.6A CN110773573B (zh) | 2019-11-06 | 2019-11-06 | 一种板形调控功效系数实测数据处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911077070.6A CN110773573B (zh) | 2019-11-06 | 2019-11-06 | 一种板形调控功效系数实测数据处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110773573A CN110773573A (zh) | 2020-02-11 |
| CN110773573B true CN110773573B (zh) | 2021-03-02 |
Family
ID=69389560
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911077070.6A Active CN110773573B (zh) | 2019-11-06 | 2019-11-06 | 一种板形调控功效系数实测数据处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110773573B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111940516B (zh) * | 2020-08-25 | 2022-03-11 | 中冶赛迪工程技术股份有限公司 | 一种确定轧机板形调控能力的方法 |
| CN112676376B (zh) * | 2020-12-30 | 2021-12-14 | 南通铭泰机床有限公司 | 一种提高数字控制的卷板机的控制精度的方法及系统 |
| CN112845617B (zh) * | 2021-01-05 | 2022-07-08 | 武汉钢铁有限公司 | 一种热轧带钢的板型控制方法及装置 |
| CN112916624A (zh) * | 2021-01-29 | 2021-06-08 | 华北电力大学（保定） | 一种ucm轧机板形执行机构调控功效系数获取方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10291013A (ja) * | 1997-04-16 | 1998-11-04 | Nisshin Steel Co Ltd | 冷間タンデム圧延における形状制御方法 |
| JP2001314909A (ja) * | 2000-05-02 | 2001-11-13 | Nkk Corp | 薄板圧延の板形状制御方法 |
| CN101920269B (zh) * | 2009-06-11 | 2012-10-31 | 鞍钢股份有限公司 | 冷轧机板形控制执行器调控功效系数优化方法 |
| CN102172639B (zh) * | 2010-12-30 | 2013-01-30 | 东北大学 | 一种冷轧机工作辊弯辊超限的动态替代调节方法 |
| CN102366760B (zh) * | 2011-09-13 | 2014-04-09 | 中冶南方工程技术有限公司 | 冷轧带钢板形控制执行器调控功效系数在线解耦方法 |
| CN102581026B (zh) * | 2012-01-18 | 2013-11-20 | 中冶南方工程技术有限公司 | 冷轧带钢板形横向整体优化控制方法 |
-
2019
- 2019-11-06 CN CN201911077070.6A patent/CN110773573B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110773573A (zh) | 2020-02-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110773573B (zh) | 一种板形调控功效系数实测数据处理方法 | |
| CN1311922C (zh) | 确定额定值和中间量的计算机辅助方法及轧机机列 | |
| CN102941232A (zh) | 一种热连轧精轧过程控制方法 | |
| CN104785543B (zh) | 一种基于滑动平均滤波的热轧带钢凸度反馈控制方法 | |
| CN102029294A (zh) | 冷轧带钢横向厚差控制方法 | |
| CN105013832A (zh) | 一种兼顾轧制能耗与良好板形的热轧带钢负荷分配方法 | |
| CN1625721A (zh) | 控制工业过程的方法 | |
| CN103506404A (zh) | 一种带钢精轧过程辊缝的控制方法 | |
| CN102189117A (zh) | 基于横向性能检测的冷轧带钢平直度前馈控制方法 | |
| CN104942019A (zh) | 一种带钢冷轧过程宽度自动控制方法 | |
| CN110929347A (zh) | 一种基于梯度提升树模型的热连轧带钢凸度预测方法 | |
| CN106825068B (zh) | 一种轧制过程带钢表面粗糙度的预测方法 | |
| CN108817103B (zh) | 一种轧钢模型钢族层别分类优化方法 | |
| CN111570533B (zh) | 热轧卷取温度与终轧温度对平整变形抗力影响预测方法 | |
| CN1091008C (zh) | 基于板形板厚协调规律的板带轧制过程互联控制方法 | |
| CN101320031B (zh) | 一种奥氏体不锈钢精密钢带的终轧变形量预测方法 | |
| CN110802114B (zh) | 一种冷轧板带轧制力的计算方法 | |
| CN105855297A (zh) | 一种提高热轧首块无取向硅钢头部厚度精度的控制方法 | |
| CN105750336B (zh) | 一种提高森吉米尔轧机非稳态轧制阶段轧制力预设定精度的方法 | |
| CN113290061A (zh) | 一种基于锆合金的多轧程多道次变形制度设计方法 | |
| CN105855298B (zh) | 一种高精度含硼钢热轧轧制力计算方法 | |
| CN111790762A (zh) | 一种基于随机森林的热轧带钢轧制力设定方法 | |
| CN111411215B (zh) | 一种多钢坯对象的炉温综合决策方法 | |
| CN112387791A (zh) | 冷轧带钢的轧制温度确定方法及确定系统 | |
| CN114733913A (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 |