CN114417236A - Data evaluation-based quality optimization control method for steel rolled product - Google Patents

Data evaluation-based quality optimization control method for steel rolled product Download PDF

Info

Publication number
CN114417236A
CN114417236A CN202210052368.7A CN202210052368A CN114417236A CN 114417236 A CN114417236 A CN 114417236A CN 202210052368 A CN202210052368 A CN 202210052368A CN 114417236 A CN114417236 A CN 114417236A
Authority
CN
China
Prior art keywords
value
calculated
strip
rolling
data
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.)
Granted
Application number
CN202210052368.7A
Other languages
Chinese (zh)
Other versions
CN114417236B (en
Inventor
彭文
刘烁
陈世译
武文腾
李越朗
韩璧保
张鸿源
孙杰
邸洪双
张殿华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN202210052368.7A priority Critical patent/CN114417236B/en
Publication of CN114417236A publication Critical patent/CN114417236A/en
Application granted granted Critical
Publication of CN114417236B publication Critical patent/CN114417236B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B3/00Rolling materials of special alloys so far as the composition of the alloy requires or permits special rolling methods or sequences ; Rolling of aluminium, copper, zinc or other non-ferrous metals
    • B21B3/02Rolling special iron alloys, e.g. stainless steel
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06395Quality analysis or management
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/04Manufacturing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Computational Mathematics (AREA)
  • Operations Research (AREA)
  • Mathematical Optimization (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Mechanical Engineering (AREA)
  • Manufacturing & Machinery (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Control Of Metal Rolling (AREA)

Abstract

The invention belongs to the technical field of steel rolling, and particularly relates to a quality optimization control method of a steel rolled product based on data evaluation. According to the method, data evaluation is carried out on the samples on the basis through data acquisition and processing of multiple samples, and the samples which are applicable to model correction and have the highest accuracy are screened out and used as source data; the correction of the model core parameters, such as roll gap, rolling force and rolling speed, directly related to the product quality in the steel production process is realized through model recalculation; in the correction process, the smooth coefficient is optimally selected, the correction efficiency is guaranteed, meanwhile, the prediction precision of the model is improved, the prediction result of the model is enabled to be faster and more accurate to be close to an actual measurement value, the control effect of quality indexes such as the thickness of steel products is improved, and the purpose of improving the high-quality control of the products in the same batch is finally achieved.

Description

一种基于数据评估的钢铁轧制产品质量优化控制方法A quality control method for steel rolling product quality optimization based on data evaluation

技术领域technical field

本发明属于钢铁轧制技术领域,具体涉及一种基于数据评估的钢铁轧制产品质量优化控制方法。The invention belongs to the technical field of steel rolling, in particular to a method for optimizing the quality of steel rolling products based on data evaluation.

背景技术Background technique

产品厚度是钢铁生产过程中十分重要的质量指标,高精度的头部厚度控制精度是带钢全长厚度控制的基础,高精度的头部厚度有助于快速进入后续的厚度自动控制(AGC)过程,从而保证带钢全长的厚度控制精度。在实际板带生产过程中,轧线上的测厚仪一般安装在最末机架的出口,由于安装环境、投资成本等的限制,在串行的机架之间一般不安装测厚仪,导致带钢头部厚度不能实时进行测量,当得到测量数据之后,带钢头部已完成生产过程,此时带钢的头部厚度已经确定,不能再进行修正;如果出现头部厚度异常,则会出现较大长度范围的头部质量异议,从而影响带钢全长的头部控制精度。若不进行及时的调整,下一带钢的厚度控制精度也会受到影响,从而影响整个生产批次的产品质量。在实际生产过程中,带钢通过测量仪表时,得到的厚度相关的实际检测数据经由数据采集系统进行采集和存储,基于实测过程数据进行模型修正,对模型中关键参数进行优化调整,能够提升产品的厚度控制精度。但在实际生产过程中,采集数据的准确度易受到穿带稳定性的影响,从而直接影响到模型修正的质量,因此必须对数据进行评估,筛选得到质量最高的生产过程数据才可以进一步供模型修正使用。Product thickness is a very important quality index in the steel production process. High-precision head thickness control accuracy is the basis for full-length strip thickness control. High-precision head thickness helps to quickly enter the subsequent automatic thickness control (AGC). process, so as to ensure the thickness control accuracy of the full length of the strip. In the actual strip production process, the thickness gauge on the rolling line is generally installed at the exit of the last stand. Due to the limitation of the installation environment and investment cost, the thickness gauge is generally not installed between the serial stands. As a result, the thickness of the strip head cannot be measured in real time. When the measurement data is obtained, the strip head has completed the production process. At this time, the strip head thickness has been determined and cannot be corrected. If the head thickness is abnormal, the There will be head quality objections in a larger length range, which will affect the head control accuracy of the full length of the strip. If the adjustment is not made in time, the thickness control accuracy of the next strip will also be affected, thereby affecting the product quality of the entire production batch. In the actual production process, when the strip passes through the measuring instrument, the actual detection data related to the thickness obtained is collected and stored through the data acquisition system, and the model is corrected based on the measured process data, and the key parameters in the model are optimized and adjusted, which can improve the product. thickness control accuracy. However, in the actual production process, the accuracy of the collected data is easily affected by the stability of the belt threading, which directly affects the quality of the model correction. Therefore, the data must be evaluated, and the highest quality production process data can be further used for the model. Correct use.

在板带轧制厚度控制方面,前人做了大量的研究工作;中国专利CN101804420B,提出了一种热轧薄板生产中精轧厚度控制的方法,针对于进入闭环控制的过程,设计了AGC串联双环系统,实现AGC自动辊缝控制功能;中国专利CN104741388A,提出一种热连轧精轧厚度控制方法,将Smith预估补偿引入了监控AGC控制系统,用GM方法来直接对轧机的辊缝进行软测量,显著提高了控制系统的响应轧制速度、稳定性和控制精度。高蕾(热加工工艺,2013,42(11),92-95)依据轧制基本理论,优化修正了厚度控制模型,提高厚度预测精度。王建(中南大学学报(自然科学版),2014,45(10),3398-3407)使用了自学习的方法针对某厂热连轧精轧机组预设定模型的自学习模块进行研究,提升了轧制力预报精度,为厚度控制提供了基础;钱京学(河南冶金,2019,27(153),49-53)分析了厚度控制系统自学习模型采样规则对头部厚度控制精度的影响,通过修改厚度自学习头部采样的时间提升模型自学习的效果。In the aspect of strip rolling thickness control, predecessors have done a lot of research work; Chinese patent CN101804420B proposes a method for finishing rolling thickness control in hot-rolled sheet production. For the process of entering closed-loop control, an AGC series connection is designed. Double-loop system realizes the automatic roll gap control function of AGC; Chinese patent CN104741388A, proposes a thickness control method for hot tandem finishing rolling, introduces Smith prediction compensation into the monitoring AGC control system, and uses the GM method to directly control the roll gap of the rolling mill. The soft measurement significantly improves the response rolling speed, stability and control precision of the control system. Gao Lei (Hot Working Technology, 2013, 42(11), 92-95) based on the basic theory of rolling, optimized and revised the thickness control model to improve the thickness prediction accuracy. Wang Jian (Journal of Central South University (Natural Science Edition), 2014, 45(10), 3398-3407) used the self-learning method to study the self-learning module of the preset model of a hot tandem finishing mill in a factory, and improved the The rolling force prediction accuracy provides a basis for thickness control; Qian Jingxue (Henan Metallurgy, 2019, 27(153), 49-53) analyzed the influence of the sampling rules of the self-learning model of the thickness control system on the thickness control accuracy of the head. The effect of model self-learning is improved by modifying the time of thickness self-learning head sampling.

以上研究过程中,通过AGC系统优化的方式提升厚度控制控制精度,必须在带钢通过测厚仪之后,进入厚度控制闭环后,根据厚度的实测值偏差才能实施,对带钢头部的厚度精度控制效果难以保证;头部厚度的控制必需通过模型修正的方式来实现,但在传统的模型修正过程中,未对生产过程中采集到的数据进行筛选和处理,由于生产过程中的实测数据易受到外界环境的影响,存在较大偏差,难以反映实际的生产状态,若使用误差较大的数据进行模型修正,将会产生相反的效果,由此造成的厚度控制精度不但不会得到提升,反而会下降,最终造成模型不适用,导致厚度控制效果不稳定,不利用产品的质量控制。In the above research process, to improve the control accuracy of thickness control through the optimization of the AGC system, it must be implemented according to the deviation of the actual measured value of the thickness after the strip passes through the thickness gauge and enters the closed loop of thickness control. The control effect is difficult to guarantee; the control of head thickness must be realized by model correction, but in the traditional model correction process, the data collected in the production process is not screened and processed, because the measured data in the production process is easy to Affected by the external environment, there is a large deviation, and it is difficult to reflect the actual production state. If the model is corrected using data with a large error, it will have the opposite effect. The resulting thickness control accuracy will not be improved, but will decrease, eventually causing the model to be unsuitable, resulting in an unstable thickness control effect and not utilizing the quality control of the product.

发明内容SUMMARY OF THE INVENTION

针对现有技术的不足,本发明提出一种基于数据评估的钢铁轧制产品质量优化控制方法,应用于将钢铁中间坯进行精轧,制备成成品带钢的过程,如图1所示,具体方案流程如下:In view of the deficiencies of the prior art, the present invention proposes a method for optimizing the quality of rolled steel products based on data evaluation, which is applied to the process of finishing intermediate steel billets to prepare finished steel strips, as shown in FIG. 1 . The program flow is as follows:

步骤1:确定穿带完成后的轧制过程的样本个数l,并确定每个样本中包含的目标采样点个数N。Step 1: Determine the number l of samples in the rolling process after the strip is finished, and determine the number N of target sampling points included in each sample.

本发明在穿带过程采集1个样本,在穿带完成后的轧制过程中采集l个样本。每个样本包含N个目标采样点,N的数值根据要轧制的钢铁产品的中间坯尺寸和成品带钢的尺寸确定。具体来说,N的数值可由以下方法确定:In the present invention, one sample is collected during the belt threading process, and one sample is collected during the rolling process after the belt threading is completed. Each sample contains N target sampling points, and the value of N is determined according to the intermediate billet size of the steel product to be rolled and the size of the finished strip. Specifically, the value of N can be determined by the following methods:

Figure BDA0003474779710000021
Figure BDA0003474779710000021

穿带完成后的轧制过程的样本个数l可为4-8,在后续步骤中要从这l个样本中筛选一个最优样本,并利用其中的数据测定值进行计算。The number of samples l in the rolling process after the strip is completed can be 4-8. In the subsequent steps, an optimal sample should be selected from the l samples, and the data measurement value therein should be used for calculation.

步骤2:带钢依次通过各机架开始穿带过程,对生产过程中产生的实测数据进行采集和存储:Step 2: The strip starts to pass through each rack in turn, and the measured data generated in the production process is collected and stored:

步骤2.1:穿带过程中带钢头部数据的采集和存储:Step 2.1: Collection and storage of strip head data during belt threading:

设轧制采用的机架数目为n。在精轧开始之前,模型设定系统会根据生产目标值的要求,计算出轧制规程,即各机架的压下率、轧制速度、辊缝、轧制力等参数的设定值。Let the number of stands used for rolling be n. Before finishing rolling starts, the model setting system will calculate the rolling schedule according to the requirements of the production target value, that is, the set values of parameters such as reduction ratio, rolling speed, roll gap and rolling force of each stand.

之后开始进行穿带和轧制过程。带钢头部通过第1机架,穿带过程开始,带钢沿轧制方向依次通过各机架和各机架处设置的测量仪表。Then the threading and rolling process begins. The head of the strip passes through the first frame, the strip threading process begins, and the strip passes through each frame and the measuring instruments set at each frame in turn along the rolling direction.

在穿带过程中,按照固定的数据采样周期,依次对带钢头部通过测量仪表后的采样点进行计数,同时采集记录测量仪表测出的该采样点的实测数据;当通过测量仪表的采样点计数等于目标采样点数目N时,计算采集到的实测数据的平均值,并将实测数据的平均值进行存储作为测定值,作为穿带过程的实测数据样本。In the process of belt threading, according to the fixed data sampling period, the sampling points after the strip head passes through the measuring instrument are counted in turn, and the measured data of the sampling point measured by the measuring instrument are collected and recorded at the same time; When the point count is equal to the number of target sampling points N, the average value of the collected measured data is calculated, and the average value of the measured data is stored as the measured value, which is used as the measured data sample of the belt passing process.

为避免穿带过程中,带钢进入各机架时咬钢冲击导致的测量值波动,对于通过测量仪表的前m个(m可取3-5)采样点不进行采集,而是从第m+1个采样点开始记录并存储,穿带过程中样本共包含N个采样点(即第m+1到第m+N个采样点)。In order to avoid the fluctuation of the measurement value caused by the impact of the steel bite when the strip enters each rack during the belt threading process, the first m sampling points (m can be 3-5) passing through the measuring instrument are not collected, but are collected from the m+th sampling point. 1 sampling point starts to be recorded and stored, and the sample contains N sampling points (ie m+1th to m+Nth sampling points) in the process of tape passing.

依次在各个机架处进行上述采集和存储过程,直至带钢头部通过最末测量仪表,并且在最末测量仪表处也完成上述采集和存储过程。The above collection and storage process is performed at each rack in sequence until the strip head passes through the last measuring instrument, and the above collection and storage process is also completed at the last measuring instrument.

穿带过程中,采集和存储的数据至少包括穿带过程带钢头部通过各个机架处的测量仪表时的轧制速度,还可以包括各个机架处的辊缝、轧制力,以及各个机架之间的活套角度等等。During the belt threading process, the data collected and stored at least include the rolling speed when the strip head passes through the measuring instruments at each stand, and may also include the roll gap at each stand, rolling force, and Looper angles between racks, etc.

步骤2.2:穿带完成后的轧制过程中带钢数据的采集和存储:Step 2.2: Collection and storage of strip data during rolling after strip threading:

带钢头部通过最末测量仪表,完成穿带过程及穿带过程中的数据采集、存储后,开始对穿带完成后的轧制过程中带钢实测数据进行采集和存储:按照固定的数据采样周期,对通过测量仪表的采样点进行计数,同时采集记录测量仪表测出的该采样点的实测数据,当通过测量仪表的采样点计数等于目标采样点数目N时,分别计算当前样本内的所有实测数据的平均值,并将当前样本内实测数据平均值进行存储作为测定值,样本数目记为1;而后,采样点清零,重新开始计数,样本数目+1,以此类推,直至完成l个样本的数据采集和存储过程。The head of the strip passes through the final measuring instrument, and after completing the data collection and storage during the strip threading process and the strip threading process, it starts to collect and store the strip steel measured data during the rolling process after the strip threading is completed: according to the fixed data Sampling period: Count the sampling points passing through the measuring instrument, and collect and record the actual measured data of the sampling point measured by the measuring instrument. When the number of sampling points passing through the measuring instrument is equal to the number of target sampling points N, calculate the current sample respectively. The average value of all measured data, and the average value of the measured data in the current sample is stored as the measured value, and the number of samples is recorded as 1; then, the sampling point is cleared, and the counting is restarted, the number of samples + 1, and so on, until the completion l sample data collection and storage process.

穿带完成后的轧制过程中,采集和存储的数据至少包括各个机架处的轧制速度、辊缝、轧制力和活套角度,以及最末机架出口的带钢温度、带钢厚度和带钢宽度。During the rolling process after strip threading, the data collected and stored at least include the rolling speed, roll gap, rolling force and looper angle at each stand, as well as the strip temperature and strip temperature at the exit of the last stand. Thickness and strip width.

步骤3:对已采集样本中的数据进行有效性评估和筛选Step 3: Validation and screening of data in collected samples

步骤3.1:对穿带过程中和穿带完成后采集到的样本中的数据进行有效性评估Step 3.1: Evaluate the validity of the data from the samples collected during and after the threading process

若样本中的任一数据超出有效性区间范围,则转至步骤7,系统报警,过程结束;若样本中的所有数据均在有效性区间范围内,则转至步骤3.2。If any data in the sample exceeds the validity interval, go to step 7, the system alarms, and the process ends; if all the data in the sample is within the validity interval, go to step 3.2.

步骤3.2:从步骤2.2中得到的l个样本数据中筛选最优样本,作为用于计算的样本。筛选最优样本可采用如下方法:Step 3.2: Screen the optimal sample from the l sample data obtained in step 2.2 as a sample for calculation. The optimal samples can be screened by the following methods:

筛选样本数据波动程度最小的样本为最优样本,样本数据波动程度δaim的计算公式如下:The sample with the smallest fluctuation degree of sample data is selected as the optimal sample. The calculation formula of the fluctuation degree of sample data δ aim is as follows:

Figure BDA0003474779710000031
Figure BDA0003474779710000031

其中,选择p个数据种类作为用于计算样本数据波动程度的指标,优选地,可采用最末机架出口的带钢温度、带钢厚度、带钢宽度、各个机架处的活套角度4个数据种类作为指标;δj为样本中第j项指标的波动程度,式中cact,j,k为第j项指标的第k个测定值,caim,j为第j项指标的目标值;Mj为各指标对应的数据测定值个数,比如上述的最末机架出口的带钢温度、带钢厚度、带钢宽度都各只有1个数据测定值,而由于有n个机架,机架之间的活套数目为n-1个,则各个机架处的活套角度这一指标则有n-1个数据测定值。Among them, select p data types as the index for calculating the fluctuation degree of the sample data, preferably, the strip temperature, strip thickness, strip width, and looper angle at each rack at the exit of the last rack can be used. δj is the fluctuation degree of the jth index in the sample, where c act,j,k is the kth measured value of the jth index, c aim,j is the target of the jth index M j is the number of measured data values corresponding to each index. For example, the strip temperature, strip thickness, and strip width at the exit of the last frame have only one measured data value each, and because there are n machines The number of loopers between the racks is n-1, and the index of the looper angle at each rack has n-1 data measurement values.

步骤4:使用筛选出的最优样本的数据确定工艺参数计算值Step 4: Use the data of the selected optimal sample to determine the calculated value of the process parameters

步骤4.1确定各个机架的出口厚度计算值Step 4.1 Determine the calculated value of the outlet thickness of each rack

根据步骤3.2得到的最优样本中的最末机架出口带钢厚度、各个机架处的轧制速度和活套角度计算各机架出口厚度;Calculate the exit thickness of each stand according to the strip thickness at the exit of the last stand in the optimal sample obtained in step 3.2, the rolling speed at each stand and the looper angle;

对于i为1至n-1的第i机架处的出口厚度计算值h计算,i为:For the calculated value h of the outlet thickness at the i-th rack where i is 1 to n-1 , i is:

Figure BDA0003474779710000041
Figure BDA0003474779710000041

式中,fi为第i机架的前滑值,fn为最末机架的前滑值;前滑值可以由模型设定系统通过现有技术中的计算方法计算获得。In the formula, f i is the forward slip value of the i-th rack, and f n is the forward slip value of the last rack; the forward slip value can be calculated by the model setting system through the calculation method in the prior art.

hn、vn分别为最优样本中,最末机架也就是第n个机架处的出口带钢厚度和轧制速度测定值,vi为最优样本中,第i机架处的轧制速度测定值;lθi为根据最优样本中第i机架处的活套角度θi测定值计算出的第i机架至第i+1机架间的带钢长度,lθs为根据目标活套角度θs计算出的第i机架至第i+1机架间的带钢目标长度;最末机架出口厚度计算值h计算,n直接采用hnh n and v n are the measured values of the thickness and rolling speed of the exit strip at the last stand, that is, the n-th stand in the optimal sample, respectively, and v i is the optimal sample at the i-th stand. Measured value of rolling speed; l θi is the strip length from the ith stand to the i+1th stand calculated according to the measured value of the looper angle θ i at the ith stand in the optimal sample, and l θs is The strip steel target length between the i-th frame and the i+1-th frame is calculated according to the target looper angle θ s ; the calculated value h of the exit thickness of the last frame is calculated, and n directly adopts h n .

上述方法中,根据两机架之间的活套角度θ计算两相邻机架间的带钢长度lθ的方法如下:In the above method, the method for calculating the strip length l θ between two adjacent frames according to the looper angle θ between the two frames is as follows:

如图2所示,两机架间的带钢长度lθ为活套最高点与前一机架出口的距离AB、活套最高点与后一机架入口的距离BC之和,根据图2的示意图不难得到:As shown in Figure 2, the strip length l θ between the two racks is the sum of the distance AB between the highest point of the looper and the exit of the previous rack, and the distance BC between the highest point of the looper and the entrance of the next rack, according to Figure 2 The schematic diagram is not difficult to get:

Figure BDA0003474779710000042
Figure BDA0003474779710000042

其中L1为前一机架到活套支点的水平距离,L2为活套支点至轧制平面的高度,L为两机架之间的水平距离,RL为活套臂长度,r为活套辊半径。Among them, L 1 is the horizontal distance from the previous stand to the fulcrum of the looper, L 2 is the height from the fulcrum of the looper to the rolling plane, L is the horizontal distance between the two stands, RL is the length of the looper arm, and r is the length of the looper arm. Looper roll radius.

步骤4.2确定各个机架的出口温度计算值Step 4.2 Determine the calculated outlet temperature for each rack

将精轧出口温度偏差分配到各机架上,计算带钢通过各机架时的温度值:Distribute the temperature deviation at the exit of the finishing rolling to each stand, and calculate the temperature value of the strip passing through each stand:

第i机架的出口温度计算值T计算,i计算公式为:The calculated value T of the outlet temperature of the i-th rack is calculated, and the calculation formula of i is:

Figure BDA0003474779710000043
Figure BDA0003474779710000043

式中,Ti为步骤2.1的穿带过程中第i机架的轧制温度设定值,T实测为最优样本中最末机架出口的带钢温度测定值,T目标为最末机架出口带钢温度目标值。In the formula, T i is the set value of the rolling temperature of the i-th stand during the strip threading process in step 2.1, T is the measured value of the strip temperature at the exit of the last stand in the optimal sample, and T target is the final machine. Rack outlet strip temperature target value.

步骤4.3各机架出口宽度值确定Step 4.3 Determine the outlet width value of each rack

精轧过程忽略宽展,因此将最优样本中最末机架出口的带钢宽度测定值作为带钢通过各机架时的出口宽度值。The width spread is ignored in the finishing rolling process, so the measured value of the strip width at the exit of the last stand in the optimal sample is taken as the exit width value when the strip passes through each stand.

步骤4.4确定各机架辊缝计算值Step 4.4 Determine the calculated value of each rack roll gap

使用弹跳方程确定各机架处的辊缝计算值,各个机架处的辊缝计算值S计算计算公式如下:Use the bouncing equation to determine the calculated value of the roll gap at each frame, and the calculation formula of the calculated value S of the roll gap at each frame is as follows:

Figure BDA0003474779710000051
Figure BDA0003474779710000051

式中,h计算为所要计算的机架处的出口厚度计算值,由步骤4.1得到;P0为调零轧制力,M为轧机刚度,均为根据轧机的性能已经确定好的数值;P实测为最优样本中所要计算的机架处的轧制力测定值。In the formula, h is calculated as the calculated value of the outlet thickness at the stand to be calculated, which is obtained from step 4.1; P 0 is the zero-adjusted rolling force, M is the stiffness of the rolling mill, which are all values that have been determined according to the performance of the rolling mill; P The actual measurement is the measured value of the rolling force at the stand to be calculated in the optimal sample.

步骤4.5确定各机架轧制力计算值Step 4.5 Determine the calculated value of rolling force of each stand

由步骤4.1至步骤4.3的各个机架的出口厚度计算值、各个机架的出口温度计算值、各机架出口宽度值,以及轧制速度的测定值等条件,计算确定各机架轧制力计算值P计算;现有技术中,根据这些已知条件计算轧制力存在多种方法。From steps 4.1 to 4.3, the calculated value of the outlet thickness of each stand, the calculated value of the outlet temperature of each stand, the value of the exit width of each stand, and the measured value of the rolling speed, etc., calculate and determine the rolling force of each stand. Calculated value P is calculated ; in the prior art, there are various methods for calculating the rolling force based on these known conditions.

优选的,可以根据西姆斯公式,确定各机架轧制力计算值:Preferably, the calculated value of rolling force of each stand can be determined according to the Sims formula:

P计算=1.15σslcQPw/1000P calculation = 1.15σ s l c Q P w/1000

式中:w—所要计算的机架的出口带钢宽度,mm,为步骤4.3确定的各机架出口宽度值;In the formula: w—the exit strip width of the rack to be calculated, mm, is the value of the exit width of each rack determined in step 4.3;

σs—所要计算的机架处的变形抗力,MPa;σ s —Deformation resistance at the frame to be calculated, MPa;

Figure BDA0003474779710000052
Figure BDA0003474779710000052

a1~a6—回归系数,其值取决于钢种;a1~ a6 regression coefficient, its value depends on steel grade;

T—所要计算的机架处的热力学温度,无量纲,

Figure BDA0003474779710000053
T计算为步骤4.2确定的所要计算机架的出口温度计算值;T—the thermodynamic temperature at the rack to be calculated, dimensionless,
Figure BDA0003474779710000053
T is calculated as the calculated value of the outlet temperature of the desired computer rack determined in step 4.2;

Figure BDA0003474779710000054
—变形速率,s-1
Figure BDA0003474779710000055
其中v为最优样本中,所要计算的机架处的轧制速度测定值;
Figure BDA0003474779710000054
— deformation rate, s -1 ,
Figure BDA0003474779710000055
where v is the measured value of rolling speed at the stand to be calculated in the optimal sample;

ε—工程应变,%,

Figure BDA0003474779710000061
Δh—所要计算的机架的压下量,mm,为所计算的机架入口厚度计算值与出口厚度计算值的差值;对于第1机架,入口厚度为带钢初始厚度h0,之后每一个机架的入口厚度计算值为上一个机架的出口厚度计算值,各机架出口厚度计算值按照步骤4.1获得;ε—engineering strain, %,
Figure BDA0003474779710000061
Δh—the reduction of the frame to be calculated, mm, is the difference between the calculated value of the inlet thickness of the frame and the calculated value of the outlet thickness; for the first frame, the inlet thickness is the initial thickness h 0 of the strip, and then The calculated value of the inlet thickness of each rack is the calculated value of the outlet thickness of the previous rack, and the calculated value of the outlet thickness of each rack is obtained according to step 4.1;

σ0—在变形温度1000℃,

Figure BDA0003474779710000062
ε=0.4时的变形抗力,MPa;σ 0 — at the deformation temperature of 1000°C,
Figure BDA0003474779710000062
Deformation resistance at ε=0.4, MPa;

lc—所要计算的机架考虑压扁后的接触弧长,mm,

Figure BDA0003474779710000063
l c —Contact arc length after crushing is considered for the frame to be calculated, mm,
Figure BDA0003474779710000063

R′—所要计算的机架的轧辊压扁半径,mm;

Figure BDA0003474779710000064
R'—roll flattening radius of the stand to be calculated, mm;
Figure BDA0003474779710000064

R—所要计算的机架的轧辊半径,mm,为已知量;R—the roll radius of the stand to be calculated, mm, is a known quantity;

QP—所要计算的机架的应力状态影响系数:Q P —Stress state influence coefficient of the rack to be calculated:

Figure BDA0003474779710000065
Figure BDA0003474779710000065

hm—带钢平均厚度,mm,为所要计算的机架处入口厚度计算值与出口厚度计算值的平均值;h m — the average thickness of the strip, mm, is the average value of the calculated value of the inlet thickness and the calculated value of the outlet thickness at the frame to be calculated;

b0~b4为与机架有关的回归系数。b 0 to b 4 are regression coefficients related to the rack.

需要说明的是,上述步骤4.4、4.5中,为了描述方便,没有加入区分不同机架的下标i,但不难理解,辊缝和轧制力计算值都是对应每个机架,根据各个机架的各个计算值、测定值和已知参数分别进行计算。It should be noted that, in the above steps 4.4 and 4.5, for the convenience of description, the subscript i for distinguishing different stands is not added, but it is not difficult to understand that the calculated values of roll gap and rolling force are corresponding to each stand. The calculated values, measured values and known parameters of the rack are calculated separately.

步骤5:确定各机架辊缝模型、轧制力模型和轧制速度模型修正系数的计算值Step 5: Determine the calculated values of the correction coefficients of the roll gap model, rolling force model and rolling speed model of each stand

步骤5.1:计算辊缝模型修正系数Step 5.1: Calculate the Roll Gap Model Correction Factor

使用步骤4.4计算得到的辊缝计算值S计算,步骤3.2筛选得到的最优样本的辊缝测定值S实测,计算辊缝模型修正系数:Use the calculated value S of the roll gap calculated in step 4.4 to calculate , and the measured value S of the roll gap of the optimal sample screened in step 3.2 to be measured , and calculate the correction coefficient of the roll gap model:

辊缝模型修正系数计算公式:

Figure BDA0003474779710000066
The calculation formula of the correction coefficient of the roll gap model:
Figure BDA0003474779710000066

步骤5.2:计算轧制力模型修正系数Step 5.2: Calculate the rolling force model correction factor

使用步骤4.5得到的轧制力计算值P计算,步骤3.2筛选得到的最优样本的轧制力测定值P实测,计算轧制力模型修正系数:Use the calculated value P of the rolling force obtained in step 4.5 to calculate , and the measured value P of the rolling force of the optimal sample screened in step 3.2 is measured , and the correction coefficient of the rolling force model is calculated:

轧制力模型修正系数计算公式:

Figure BDA0003474779710000071
The calculation formula of the correction coefficient of the rolling force model:
Figure BDA0003474779710000071

步骤5.3:计算轧制速度模型修正系数Step 5.3: Calculate the rolling speed model correction factor

使用步骤2.1得到的穿带过程中轧制速度测定值v实测-穿带过程,以及步骤3.2得到的穿带完成后的轧制过程的最优样本中的轧制速度测定值v,计算轧制速度模型修正系数

Figure BDA0003474779710000072
Use the measured value v of the rolling speed during the strip threading process obtained in step 2.1 - the strip threading process , and the measured value v of the rolling speed in the optimal sample of the rolling process after the strip threading obtained in step 3.2, to calculate the rolling Speed model correction factor
Figure BDA0003474779710000072

轧制速度模型修正系数计算公式:

Figure BDA0003474779710000073
The calculation formula of the correction coefficient of the rolling speed model:
Figure BDA0003474779710000073

步骤6:各机架模型修正系数更新Step 6: Update the correction factor of each rack model

步骤6.1:根据修正系数旧值与修正系数计算值间的偏离程度计算平滑系数α;计算公式如下:Step 6.1: Calculate the smoothing coefficient α according to the degree of deviation between the old value of the correction coefficient and the calculated value of the correction coefficient; the calculation formula is as follows:

Figure BDA0003474779710000074
Figure BDA0003474779710000074

式中:Δold为模型修正系数旧值,包括辊缝、轧制力、轧制速度的模型修正系数旧值

Figure BDA0003474779710000075
Figure BDA0003474779710000076
Figure BDA0003474779710000077
Δ计算为步骤5.1-5.3计算得到的模型修正系数计算值,包括
Figure BDA0003474779710000078
Figure BDA0003474779710000079
Δmax为对应的模型修正系数最大值,包括辊缝、轧制力、轧制速度的模型修正系数最大值
Figure BDA00034747797100000710
Figure BDA00034747797100000711
Figure BDA00034747797100000712
Δmin为对应的模型修正系数最小值,包括辊缝、轧制力、轧制速度的模型修正系数最小值
Figure BDA00034747797100000713
Figure BDA00034747797100000714
In the formula: Δ old is the old value of the model correction coefficient, including the old value of the model correction coefficient of the roll gap, rolling force and rolling speed
Figure BDA0003474779710000075
Figure BDA0003474779710000076
and
Figure BDA0003474779710000077
Δ is calculated as the model correction coefficient calculated in steps 5.1-5.3, including
Figure BDA0003474779710000078
and
Figure BDA0003474779710000079
Δmax is the maximum value of the corresponding model correction coefficient, including the maximum value of the model correction coefficient of roll gap, rolling force and rolling speed
Figure BDA00034747797100000710
Figure BDA00034747797100000711
and
Figure BDA00034747797100000712
Δmin is the minimum value of the corresponding model correction coefficient, including the minimum value of the model correction coefficient of roll gap, rolling force and rolling speed
Figure BDA00034747797100000713
and
Figure BDA00034747797100000714

对于辊缝、轧制力、轧制速度三个修正系数,均有若Δ计算≥Δmax,则Δ计算替代并成为新的Δmax值,若Δ计算≤Δmin,则Δ计算替代并成为新的Δmin值,然后再用于平滑系数α的计算;ω为比例系数,可取0.5-0.9。For the three correction coefficients of roll gap, rolling force and rolling speed, if Δcalculation ≥Δmax , then Δcalculation is replaced and becomes the new Δmax value; if Δcalculation≤Δmin , then Δcalculation is replaced and becomes the new Δmax value . The new Δmin value is then used for the calculation of the smoothing coefficient α; ω is the proportional coefficient, which can be 0.5-0.9.

步骤6.2:模型修正系数新值计算Step 6.2: Calculate the new value of the model correction coefficient

采用平滑的方式对模型系数进行修正,平滑后的模型修正系数新值Δnew按下式计算:The model coefficients are corrected in a smoothing manner, and the new value of the smoothed model correction coefficient Δnew is calculated as follows:

Δnew=Δold+α·(Δ计算old) Δnew = Δold +α·( Δcalculation - Δold )

步骤6.3:将计算得到的模型修正系数新值传递给模型设定系统,当下一块带钢进行轧制规程计算时,将辊缝修正系数与模型设定系统的设定结果相加作为新的辊缝设定值,将轧制力修正系数与模型设定结果相乘作为新的轧制力设定值,将轧制速度模型修正系数与模型设定结果相乘作为新的轧制速度设定值,供带钢生产使用。Step 6.3: Transfer the calculated new value of the model correction coefficient to the model setting system, and add the roll gap correction coefficient and the setting result of the model setting system as the new roll when the next strip is calculated for the rolling schedule. The setting value of the rolling force is multiplied by the rolling force correction factor and the model setting result as the new rolling force setting value, and the rolling speed model correction factor is multiplied by the model setting result as the new rolling speed setting. value for strip production.

步骤6.1-6.3的各机架模型修正系数计算更新过程可参见图4。See Figure 4 for the calculation and update process of the correction coefficients of each rack model in steps 6.1-6.3.

对于每次轧制,上一次的模型修正系数新值就变成这次的旧值,通过上述方法进行一次修正和更新。模型修正系数及模型修正系数最大值、最小值在初始时可以取默认值,比如对于辊缝,初始的模型修正系数值可以取0,作为初始的

Figure BDA0003474779710000081
初始的
Figure BDA0003474779710000082
可以分别取0.05和-0.05;对于轧制力和轧制速度,则可以初始模型修正系数取1,初始的最大值和最小值分别取1.05和0.95,之后在每次轧制过程中,不断对模型修正系数进行更新,对轧制规程持续进行优化。For each rolling, the new value of the model correction coefficient of the previous time becomes the old value of this time, and a correction and update are performed by the above method. The model correction coefficient and the maximum and minimum value of the model correction coefficient can take the default values at the beginning. For example, for the roll gap, the initial model correction coefficient value can take 0 as the initial value.
Figure BDA0003474779710000081
Initially
Figure BDA0003474779710000082
It can be taken as 0.05 and -0.05 respectively; for the rolling force and rolling speed, the initial model correction coefficient can be taken as 1, and the initial maximum and minimum values are taken as 1.05 and 0.95 respectively. Model correction factors are updated to continuously optimize the rolling schedule.

步骤7:结束。Step 7: End.

本发明通过多样本的数据采集和处理,在此基础上对样本进行数据评估,筛选出适用模型修正用的准确度最高的样本,以此为源数据;通过模型再计算,实现对钢铁生产过程中的与产品质量直接相关的模型核心参数辊缝、轧制力、轧制速度的修正;在修正过程中,对平滑系数进行最优选取,保证修正效率的同时,提高模型的预测精度,使模型的预测结果更快更准确的接近于实测值,从而实现提升钢铁产品厚度等质量指标的控制效果,最终达到提高同批次产品高质化控制的目的。The invention collects and processes the data of multiple samples, and on this basis, evaluates the data of the samples, selects the samples with the highest accuracy for model correction, and uses this as the source data; Correction of the core parameters of the model, which are directly related to the product quality, roll gap, rolling force and rolling speed; in the correction process, the smoothing coefficient is optimally selected to ensure the correction efficiency and improve the prediction accuracy of the model. The prediction result of the model is faster and more accurate close to the measured value, so as to achieve the control effect of improving the quality indicators such as the thickness of steel products, and finally achieve the purpose of improving the high-quality control of the same batch of products.

附图说明Description of drawings

图1为本发明中方法的流程示意图。FIG. 1 is a schematic flow chart of the method in the present invention.

图2为两相邻机架间的带钢长度计算涉及的各参数示意图。Figure 2 is a schematic diagram of various parameters involved in the calculation of the strip length between two adjacent racks.

图3为本发明实施例中的带钢轧制过程精轧区设备仪表布置图。FIG. 3 is a layout diagram of equipment and instruments in the finishing rolling area of the strip rolling process in the embodiment of the present invention.

图4为本发明的各机架模型修正系数计算流程图。FIG. 4 is a flow chart of calculating the correction coefficient of each rack model according to the present invention.

图5为本发明实施例中设定参数修正后,对于厚度精度的提升效果图。FIG. 5 is a diagram showing the effect of improving the thickness accuracy after the setting parameters are corrected in the embodiment of the present invention.

具体实施方式Detailed ways

下面结合实施例对本发明的方案进行详细说明:Below in conjunction with embodiment, the scheme of the present invention is described in detail:

本实施例对Q235A钢种的热连轧精轧过程的轧制产品质量进行优化:热连轧精轧区由7个机架组成,第7机架为最末机架,轧辊半径为380mm,机架间距离为5000mm。如图3所示,带钢9由精轧区入口依次穿过7个机架(图示中2为第一机架)后,到达精轧出口,完成穿带过程。精轧区入口、精轧区出口分别布置有测温仪1和测温仪6测量带钢的温度,精轧出口即最末机架出口设置测宽仪7测量带钢的实际宽度,精轧出口还设置测厚仪8测量带钢的实际厚度,每个机架安装有位移传感器3、压力传感器4和轧制速度传感器5,分别测量带钢通过各机架时的辊缝、轧制力和轧辊线速度(即轧制速度),相邻机架之间活套10安装有角度编码器,测量生产过程中的活套角度;出口测厚仪8为最末测量仪表。This embodiment optimizes the quality of rolled products in the hot tandem finishing rolling process of Q235A steel: the hot tandem finishing area consists of 7 stands, the 7th stand is the last stand, and the roll radius is 380mm. The distance between racks is 5000mm. As shown in FIG. 3 , the strip 9 passes through 7 stands (2 is the first stand in the figure) from the entrance of the finishing rolling area in sequence, and then reaches the finishing rolling exit to complete the strip threading process. Thermometer 1 and thermometer 6 are respectively arranged at the entrance of the finishing rolling area and the exit of the finishing rolling area to measure the temperature of the strip, and a width measuring instrument 7 is installed at the exit of the finishing rolling, that is, the exit of the last stand to measure the actual width of the strip. The outlet is also equipped with a thickness gauge 8 to measure the actual thickness of the strip, and each stand is equipped with a displacement sensor 3, a pressure sensor 4 and a rolling speed sensor 5 to measure the roll gap and rolling force when the strip passes through each stand. And the linear speed of the roll (ie rolling speed), an angle encoder is installed on the looper 10 between the adjacent stands to measure the angle of the looper in the production process; the outlet thickness gauge 8 is the final measuring instrument.

当前产品的入口温度为1020℃,精轧中间坯厚度为32.0mm,中间坯宽度为1200mm;活套角度目标值为20°,成品带钢的温度目标值为880℃,厚度目标值为2.00mm,宽度目标值为1200mm;数据采样周期为200ms,穿带速度为10×103mm/s。The inlet temperature of the current product is 1020°C, the thickness of the finishing rolled intermediate billet is 32.0mm, and the width of the intermediate billet is 1200mm; the target value of the looper angle is 20°, the target value of the temperature of the finished strip is 880°C, and the target value of the thickness is 2.00mm , the width target value is 1200mm; the data sampling period is 200ms, and the tape threading speed is 10×10 3 mm/s.

生产过程使用的轧制规程由模型设定系统计算得到,带钢按照此轧制规程开始生产过程,轧制规程如表1所示。The rolling schedule used in the production process is calculated by the model setting system. The strip starts the production process according to this rolling schedule. The rolling schedule is shown in Table 1.

表1Table 1

参数名称parameter name 1机架1 rack 2机架2 racks 3机架3 racks 4机架4 racks 5机架5 racks 6机架6 racks 7机架7 racks 入口厚度/mmInlet thickness/mm 32.0032.00 16.6516.65 9.449.44 5.695.69 3.903.90 3.013.01 2.382.38 出口厚度/mmOutlet thickness/mm 16.6516.65 9.449.44 5.695.69 3.903.90 3.013.01 2.382.38 2.002.00 压下率/%Reduction rate/% 48.048.0 43.343.3 39.739.7 31.531.5 22.822.8 20.920.9 16.016.0 轧制速度/m/sRolling speed/m/s 1.051.05 2.012.01 3.323.32 5.015.01 6.586.58 8.428.42 10.0010.00 温度/℃temperature/℃ 989.0989.0 971.5971.5 954.3954.3 937.4937.4 920.8920.8 904.5904.5 888.5888.5 辊缝/mmRoll gap/mm 19.4319.43 11.8911.89 7.937.93 5.305.30 3.523.52 2.782.78 1.881.88 轧制力/kNRolling force/kN 2469224692 2268022680 2144921449 1639916399 1105311053 1039710397 72987298

步骤1:确定穿带过程中的样本个数为1个,穿带完成后的轧制过程的样本个数为5个。Step 1: It is determined that the number of samples in the belt threading process is 1, and the number of samples in the rolling process after the belt threading is completed is 5.

根据带钢的中间坯和成品带钢尺寸,确定模型修正用样本中包含的目标采样点数目。The number of target sampling points included in the sample for model correction is determined based on the intermediate and finished strip dimensions of the strip.

每个样本中包含的目标采样点数目为:The number of target sampling points included in each sample is:

Figure BDA0003474779710000091
Figure BDA0003474779710000091

为避免穿带过程中,带钢进入各机架时咬钢冲击导致的测量值波动,前5个采样点不进行采集,从第6个采样点开始采集记录并存储,每个样本共包含N=20个采样点(即第6到第25个采样点)。In order to avoid the fluctuation of the measurement value caused by the impact of the steel bite when the strip enters each rack during the belt threading process, the first 5 sampling points are not collected, and the recording and storage are started from the 6th sampling point. Each sample contains a total of N = 20 sample points (ie, the 6th to 25th sample points).

步骤2:带钢依次通过各机架开始穿带过程,对生产过程中产生的实测数据进行采集和存储。Step 2: The strip starts to pass through each rack in turn, and the measured data generated in the production process is collected and stored.

步骤2.1:穿带过程中带钢头部数据的采集和存储Step 2.1: Collection and storage of strip head data during belt threading

带钢头部通过第1机架,穿带过程开始,带钢沿轧制方向依次通过各机架和各机架处设置的测量仪表;按照固定的数据采样周期,依次对带钢头部通过测量仪表后的采样点进行计数,同时采集记录测量仪表测出的该采样点的实测数据;当通过测量仪表的采样点计数等于目标采样点数目20时,计算采集到的实测数据的平均值,并将实测数据的平均值进行存储作为测定值,作为穿带过程的实测数据样本;The strip head passes through the first frame, and the strip threading process begins. The strip passes through each frame and the measuring instruments set at each frame in turn along the rolling direction; according to the fixed data sampling period, the strip head passes through in turn. The sampling points behind the measuring instrument are counted, and the measured data of the sampling point measured by the measuring instrument are collected and recorded at the same time; when the sampling point count through the measuring instrument is equal to the number of target sampling points 20, the average value of the collected measured data is calculated, And store the average value of the measured data as the measured value, as the measured data sample of the belt wearing process;

依次在各个机架处进行上述采集和存储过程,直至带钢头部通过最末测量仪表,并且在最末测量仪表处也完成上述采集和存储过程。The above collection and storage process is performed at each rack in sequence until the strip head passes through the last measuring instrument, and the above collection and storage process is also completed at the last measuring instrument.

采集得到的穿带过程中带钢头部样本数据如表2所示。其中第1机架处的活套角度是指第1和第2机架之间的活套角度,其他以此类推。The sample data of the strip steel head during the belt threading process collected are shown in Table 2. The looper angle at the first rack refers to the looper angle between the first and second racks, and so on.

表2Table 2

参数名称parameter name 1机架1 rack 2机架2 racks 3机架3 racks 4机架4 racks 5机架5 racks 6机架6 racks 7机架7 racks 轧制速度/m/sRolling speed/m/s 1.111.11 2.022.02 3.343.34 5.005.00 6.566.56 8.438.43 1010 辊缝/mmRoll gap/mm 19.4319.43 11.8911.89 7.937.93 5.305.30 3.523.52 2.782.78 1.881.88 轧制力/kNRolling force/kN 2400224002 2212222122 2088920889 1629116291 1108811088 99829982 73027302 活套角度/°Loop angle/° 14.614.6 16.816.8 17.017.0 19.419.4 20.420.4 18.218.2 //

步骤2.2:穿带完成后带钢数据的采集和存储Step 2.2: Collection and storage of strip data after strip-threading

带钢头部通过出口最末测量仪表(测厚仪8)后,开始对穿带完成后的轧制过程中带钢实测数据进行采集和存储:按照固定的数据采样周期,对通过测量仪表的采样点进行计数,同时采集记录测量仪表测出的该采样点的实测数据,当通过测量仪表的采样点计数等于目标采样点数目N时,分别计算当前样本内的所有实测数据的平均值,并将当前样本内实测数据平均值进行存储作为测定值,样本数目记为1;而后,采样点清零,重新开始计数,样本数目+1,直至完成5个样本的数据采集和存储过程。采集得到的穿带完成后的轧制过程中带钢5个样本数据如表3所示。After the strip head passes through the last measuring instrument (thickness gauge 8) at the outlet, the actual measurement data of the strip in the rolling process after the strip is completed starts to be collected and stored: according to the fixed data sampling period, The sampling point is counted, and the measured data of the sampling point measured by the measuring instrument is collected and recorded. When the sampling point count through the measuring instrument is equal to the number of target sampling points N, the average value of all the measured data in the current sample is calculated respectively, and The average value of the measured data in the current sample is stored as the measured value, and the number of samples is recorded as 1; then, the sampling point is cleared, and the counting is restarted, and the number of samples is +1 until the data collection and storage process of 5 samples is completed. The collected data of 5 samples of strip steel in the rolling process after strip threading is completed are shown in Table 3.

表3table 3

Figure BDA0003474779710000101
Figure BDA0003474779710000101

Figure BDA0003474779710000111
Figure BDA0003474779710000111

步骤3:对已采集样本中的数据进行有效性评估和筛选Step 3: Validation and screening of data in collected samples

步骤3.1:对穿带过程中和穿带完成后采集到的样本中的数据进行有效性评估Step 3.1: Evaluate the validity of the data from the samples collected during and after the threading process

若样本中的任一数据超出有效性区间范围,则转至步骤7,系统报警,过程结束;若样本中的所有数据均在有效性区间范围内,则转至步骤3.2。If any data in the sample exceeds the validity interval, go to step 7, the system alarms, and the process ends; if all the data in the sample is within the validity interval, go to step 3.2.

各工艺参数的有效性区间范围如表4所示,可见样本所有数据均在有效性区间范围内。The validity interval range of each process parameter is shown in Table 4. It can be seen that all the data of the sample are within the validity interval range.

表4Table 4

Figure BDA0003474779710000112
Figure BDA0003474779710000112

Figure BDA0003474779710000121
Figure BDA0003474779710000121

步骤3.2:从步骤2.2中得到的5个样本数据中筛选最优样本,作为用于计算的样本;具体方法为筛选样本数据波动程度最小的样本为最优样本,样本数据波动程度δaim的计算公式如下:Step 3.2: Select the optimal sample from the 5 sample data obtained in step 2.2 as the sample for calculation; the specific method is to select the sample with the smallest fluctuation degree of the sample data as the optimal sample, and calculate the fluctuation degree of the sample data δ aim The formula is as follows:

Figure BDA0003474779710000122
Figure BDA0003474779710000122

其中,选择p个数据种类作为用于计算样本数据波动程度的指标,本实施例采用最末机架出口的带钢温度、带钢厚度、带钢宽度、各个机架处的活套角度4个数据种类作为指标;δj为样本中第j项指标的波动程度,式中cact,j,k为第j项指标的第k个测定值,caim,j为第j项指标的目标值;Mj为各指标对应的数据测定值个数,对于最末机架出口的带钢温度、带钢厚度、带钢宽度都各只有1个数据测定值,而由于本实施例有7个机架,对应的有6个处于机架间的活套,则各个机架间的活套角度这一指标则有6个数据测定值。Among them, p data types are selected as indicators for calculating the fluctuation degree of the sample data. In this embodiment, the strip temperature, strip thickness, strip width, and looper angle at each rack outlet are used in this embodiment. The data type is used as the indicator; δj is the fluctuation degree of the jth indicator in the sample, where c act,j,k is the kth measured value of the jth indicator, c aim,j is the target value of the jth indicator ; M j is the number of data measurement values corresponding to each index, and there is only one data measurement value for the strip temperature, strip thickness, and strip width at the exit of the last frame, and because this embodiment has 7 machines. If there are 6 loopers between the racks, there are 6 data measurement values for the index of the looper angle between the racks.

本实施例中5个样本的数据波动程度计算结果如表5所示,样本5对应的数据波动程度为0.018,小于其他样本数据波动程序,选为最优样本。The calculation results of the data fluctuation degree of the five samples in this embodiment are shown in Table 5. The data fluctuation degree corresponding to sample 5 is 0.018, which is smaller than that of other sample data fluctuation programs, and is selected as the optimal sample.

表5table 5

样本号sample number 11 22 33 44 55 样本数据波动程度Fluctuation of sample data 0.1590.159 0.1120.112 0.0800.080 0.0440.044 0.0180.018

步骤4:使用筛选出的最优样本的数据确定工艺参数计算值Step 4: Use the data of the selected optimal sample to determine the calculated value of the process parameters

步骤4.1确定各个机架的出口厚度计算值Step 4.1 Determine the calculated value of the outlet thickness of each rack

根据步骤3.2得到的最优样本中的最末机架出口带钢厚度、各个机架处轧制速度和活套角度计算各机架出口厚度;对于i为1至6的第i机架的出口厚度计算值h计算,i为:Calculate the exit thickness of each stand according to the strip thickness at the exit of the last stand, the rolling speed at each stand and the looper angle in the optimal sample obtained in step 3.2; for the exit of the i-th stand where i is 1 to 6 Thickness calculation value h calculation, i is:

Figure BDA0003474779710000123
Figure BDA0003474779710000123

式中,fi为第i机架的前滑值,f7为最末机架的前滑值,根据模型设定系统计算得到;根据表3,h7=2.00mm,v7=10.00m/s;vi为最优样本中,第i机架处的轧制速度测定值;lθi为根据最优样本中第i机架处的活套角度θi测定值计算出的第i机架至第i+1机架间的带钢长度,l20°为根据目标活套角度20°计算出的第i机架至第i+1机架间的带钢目标长度。最末机架出口厚度计算值h计算,7直接采用h7=2.00mm。In the formula, f i is the forward sliding value of the i-th rack, and f 7 is the forward sliding value of the last rack, which is calculated according to the model setting system; according to Table 3, h 7 =2.00mm, v 7 =10.00m /s; vi is the measured value of the rolling speed at the i -th stand in the optimal sample; l θi is the i-th machine calculated from the measured value of the looper angle θ i at the i-th stand in the optimal sample The strip length from the rack to the i+1th rack, l 20° is the strip steel target length from the ith rack to the i+1th rack calculated according to the target looper angle of 20°. The calculated value h of the final frame outlet thickness is calculated, and h 7 = 2.00mm is directly used for 7.

机架间活套结构和对应的参数如图2所示。相邻机架间的带钢长度计算值如表6所示:The looper structure between racks and the corresponding parameters are shown in Figure 2. The calculated value of strip length between adjacent racks is shown in Table 6:

表6Table 6

活套编号Looper number 11 22 33 44 55 66 带钢长度/mmStrip length/mm 5007.15007.1 5008.15008.1 5007.45007.4 5008.15008.1 5007.85007.8 5007.75007.7

相邻机架间的带钢长度具体计算方法为,当活套角度为θ时,根据几何关系可以计算得到相邻机架间的带钢长度:The specific calculation method of the strip length between adjacent frames is, when the looper angle is θ, the strip length between adjacent frames can be calculated according to the geometric relationship:

Figure BDA0003474779710000131
Figure BDA0003474779710000131

以本实施例的第1机架的出口厚度计算值为例,第1机架和第2机架的活套角度为19.5°时,机架间的带钢长度为:Taking the calculated value of the outlet thickness of the first frame in this embodiment as an example, when the looper angle of the first frame and the second frame is 19.5°, the strip length between the frames is:

Figure BDA0003474779710000132
Figure BDA0003474779710000132

当目标活套角度为20.0°时,机架间的带钢长度为:When the target looper angle is 20.0°, the strip length between the frames is:

Figure BDA0003474779710000133
Figure BDA0003474779710000133

由表3可知,v1=1.15m/s,计算得到的第1机架出口厚度计算值为:It can be seen from Table 3 that v 1 =1.15m/s, and the calculated value of the thickness at the exit of the first frame is:

Figure BDA0003474779710000134
Figure BDA0003474779710000134

步骤4.2确定各个机架的出口温度计算值Step 4.2 Determine the calculated outlet temperature for each rack

将精轧出口温度偏差分配到各机架上,计算带钢通过各机架时的温度值;第i机架的出口温度计算值T计算,i计算公式为:The temperature deviation of the finishing rolling outlet is allocated to each stand, and the temperature value of the strip steel passing through each stand is calculated; the calculated value T of the outlet temperature of the i-th stand is calculated, and the calculation formula of i is:

Figure BDA0003474779710000135
Figure BDA0003474779710000135

式中,Ti为步骤2.1的穿带过程中第i机架的轧制温度设定值,T实测为最优样本(样本5)中最末机架出口的带钢温度测定值,如表3所示为877℃,T目标为最末机架出口带钢温度目标值,在本实施例中为880℃。机架数目n=7。In the formula, T i is the set value of the rolling temperature of the i-th stand during the strip threading process in step 2.1, and T is the measured value of the strip temperature at the exit of the last stand in the optimal sample (sample 5), as shown in the table 3 shows 877°C and T target is the last stand exit strip temperature target value, which in this example is 880°C. The number of racks is n=7.

以第一机架为例,根据表1,T1=989.0℃,计算得到的第一机架带钢出口温度计算值为:Taking the first rack as an example, according to Table 1, T 1 =989.0℃, the calculated value of the strip steel outlet temperature of the first rack is:

Figure BDA0003474779710000141
Figure BDA0003474779710000141

步骤4.3各机架出口宽度值确定Step 4.3 Determine the outlet width value of each rack

精轧过程忽略宽展,将最优样本5中精轧出口宽度测定值作为带钢通过各机架时的出口宽度计算值;因此,所有机架的带钢出口宽度值为:In the finishing rolling process, the width spread is ignored, and the measured value of the finishing exit width in the optimal sample 5 is taken as the calculated value of the exit width when the strip passes through each stand; therefore, the strip exit width value of all stands is:

w计算=w实测=1205mmw calculated = w measured = 1205mm

步骤4.4确定各机架辊缝计算值Step 4.4 Determine the calculated value of each rack roll gap

使用弹跳方程确定各机架处的辊缝计算值,各个机架处的辊缝计算值S计算计算公式如下:Use the bouncing equation to determine the calculated value of the roll gap at each frame, and the calculation formula of the calculated value S of the roll gap at each frame is as follows:

Figure BDA0003474779710000142
Figure BDA0003474779710000142

式中,h计算为所要计算的机架处的出口厚度计算值,由步骤4.1得到;P0为调零轧制力,本实施例中P0=8000kN,M为轧机刚度,本实施例中M=6000kN/mm;P实测为最优样本中所要计算的机架处的轧制力测定值,在本实施例中为表3中样本5的各机架轧制力测定值。In the formula, h is calculated as the calculated value of the outlet thickness at the stand to be calculated, obtained from step 4.1; P 0 is the zero-adjusted rolling force, in this embodiment, P 0 =8000kN, M is the stiffness of the rolling mill, in this embodiment M=6000kN/mm; P is the measured value of rolling force at the stand to be calculated in the optimal sample, and in this embodiment, it is the measured value of rolling force of each stand of Sample 5 in Table 3.

以第一机架为例,计算得到的机架辊缝为:Taking the first stand as an example, the calculated roll gap of the stand is:

Figure BDA0003474779710000143
Figure BDA0003474779710000143

步骤4.5确定各机架轧制力计算值Step 4.5 Determine the calculated value of rolling force of each stand

根据西姆斯公式,确定各机架轧制力计算值,各个机架处的轧制力计算值P计算为:According to Sims formula, the calculated value of rolling force of each stand is determined, and the calculated value of rolling force P at each stand is calculated as:

P计算=1.15σslcQPw/1000P calculation = 1.15σ s l c Q P w/1000

式中:w—所要计算的机架的出口带钢宽度,mm,为步骤4.3确定的各机架出口宽度值;In the formula: w—the exit strip width of the rack to be calculated, mm, is the value of the exit width of each rack determined in step 4.3;

σs—所要计算的机架处的变形抗力,MPa;σ s —Deformation resistance at the frame to be calculated, MPa;

Figure BDA0003474779710000144
Figure BDA0003474779710000144

a1~a6—回归系数,其值取决于钢种;a1~ a6 regression coefficient, its value depends on steel grade;

T—所要计算的机架处的热力学温度,无量纲,

Figure BDA0003474779710000151
T计算为步骤4.2确定的所要计算机架的出口温度计算值;T—the thermodynamic temperature at the rack to be calculated, dimensionless,
Figure BDA0003474779710000151
T is calculated as the calculated value of the outlet temperature of the desired computer rack determined in step 4.2;

Figure BDA0003474779710000152
—变形速率,s-1
Figure BDA0003474779710000153
其中v为最优样本中,所要计算的机架处的轧制速度测定值;
Figure BDA0003474779710000152
— deformation rate, s -1 ,
Figure BDA0003474779710000153
where v is the measured value of rolling speed at the stand to be calculated in the optimal sample;

ε—工程应变,%,

Figure BDA0003474779710000154
Δh—所要计算的机架的压下量,mm,为所计算的机架入口厚度计算值与出口厚度计算值的差值;对于第1机架,入口厚度为带钢初始厚度h0,也就是精轧中间坯的厚度;之后每一个机架的入口厚度计算值为上一个机架的出口厚度计算值,各机架出口厚度计算值按照步骤4.1获得;ε—engineering strain, %,
Figure BDA0003474779710000154
Δh—the reduction amount of the frame to be calculated, mm, is the difference between the calculated value of the inlet thickness of the frame and the calculated value of the outlet thickness; for the first frame, the inlet thickness is the initial thickness h 0 of the strip, also is the thickness of the finishing rolling intermediate billet; after that, the calculated value of the entrance thickness of each stand is the calculated value of the exit thickness of the previous stand, and the calculated value of the exit thickness of each stand is obtained according to step 4.1;

σ0—在变形温度1000℃,

Figure BDA0003474779710000155
ε=0.4时的变形抗力,MPa;对于本实施例的钢种为150.6MPa;σ 0 — at the deformation temperature of 1000°C,
Figure BDA0003474779710000155
Deformation resistance when ε=0.4, MPa; for the steel grade of this embodiment, it is 150.6MPa;

lc—所要计算的机架考虑压扁后的接触弧长,mm,

Figure BDA0003474779710000156
l c —Contact arc length after crushing is considered for the frame to be calculated, mm,
Figure BDA0003474779710000156

R′—所要计算的机架的轧辊压扁半径,mm;

Figure BDA0003474779710000157
R'—roll flattening radius of the stand to be calculated, mm;
Figure BDA0003474779710000157

R—所要计算的机架的轧辊半径,mm;R—roll radius of the stand to be calculated, mm;

QP—所要计算的机架的应力状态影响系数:Q P —Stress state influence coefficient of the rack to be calculated:

Figure BDA0003474779710000158
Figure BDA0003474779710000158

hm—带钢平均厚度,mm,为所要计算的机架处入口厚度计算值与出口厚度计算值的平均值;h m — the average thickness of the strip, mm, is the average value of the calculated value of the inlet thickness and the calculated value of the outlet thickness at the frame to be calculated;

b0~b4为与机架有关的回归系数。b 0 to b 4 are regression coefficients related to the rack.

对于本实施例中的钢种和机架,各个回归系数的取值如下表所示:For the steel grade and frame in this embodiment, the values of each regression coefficient are shown in the following table:

回归系数Regression coefficients a<sub>1</sub>a<sub>1</sub> a<sub>2</sub>a<sub>2</sub> a<sub>3</sub>a<sub>3</sub> a<sub>4</sub>a<sub>4</sub> a<sub>5</sub>a<sub>5</sub> a<sub>6</sub>a<sub>6</sub> 取值value 2.8782.878 3.6653.665 0.18610.1861 -0.1216-0.1216 0.37950.3795 1.4021.402 回归系数Regression coefficients b<sub>0</sub>b<sub>0</sub> b<sub>1</sub>b<sub>1</sub> b<sub>2</sub>b<sub>2</sub> b<sub>3</sub>b<sub>3</sub> b<sub>4</sub>b<sub>4</sub> 取值value 0.80490.8049 -0.3393-0.3393 0.24880.2488 0.03930.0393 0.07320.0732

按照步骤4.1-4.5,各工艺参数的计算值如表7所示:According to steps 4.1-4.5, the calculated values of each process parameter are shown in Table 7:

表7Table 7

参数名称parameter name 1机架1 rack 2机架2 racks 3机架3 racks 4机架4 racks 5机架5 racks 6机架6 racks 7机架7 racks 入口厚度计算值/mmCalculated value of inlet thickness/mm 32.0032.00 16.7316.73 9.349.34 5.705.70 3.853.85 2.992.99 2.352.35 出口厚度计算值/mmCalculated value of outlet thickness/mm 16.7316.73 9.349.34 5.705.70 3.853.85 2.992.99 2.352.35 2.002.00 温度计算值/℃Calculated temperature/℃ 988.6988.6 970.6970.6 953.0953.0 935.7935.7 918.7918.7 901.9901.9 885.5885.5 宽度计算值/mmCalculated width/mm 12051205 12051205 12051205 12051205 12051205 12051205 12051205 辊缝计算值/mmRoll gap calculation value/mm 19.4519.45 11.7511.75 7.877.87 5.295.29 3.523.52 2.692.69 1.921.92 轧制力计算值/kNCalculated value of rolling force/kN 2300223002 2423824238 2007120071 1720117201 1102011020 99829982 74877487

步骤5:确定各机架辊缝模型、轧制力模型和轧制速度模型修正系数的计算值Step 5: Determine the calculated values of the correction coefficients of the roll gap model, rolling force model and rolling speed model of each stand

步骤5.1:计算辊缝模型修正系数Step 5.1: Calculate the Roll Gap Model Correction Factor

使用步骤4.4计算得到的辊缝计算值S计算,步骤3.2筛选得到的最优样本的辊缝测定值S实测,计算辊缝模型修正系数:Use the calculated value S of the roll gap calculated in step 4.4 to calculate , and the measured value S of the roll gap of the optimal sample screened in step 3.2 to be measured , and calculate the correction coefficient of the roll gap model:

辊缝模型修正系数计算公式:

Figure BDA0003474779710000161
The calculation formula of the correction coefficient of the roll gap model:
Figure BDA0003474779710000161

以第1机架为例,辊缝模型修正系数为

Figure BDA0003474779710000162
Taking the first frame as an example, the correction coefficient of the roll gap model is
Figure BDA0003474779710000162

各机架处辊缝模型修正系数计算值如下:The calculated value of the correction coefficient of the roll gap model at each stand is as follows:

Figure BDA0003474779710000163
Figure BDA0003474779710000163

步骤5.2:计算轧制力模型修正系数Step 5.2: Calculate the rolling force model correction factor

使用步骤4.5得到的轧制力计算值P计算,步骤3.2筛选得到的最优样本的轧制力测定值P实测,计算轧制力模型修正系数:Use the calculated value P of the rolling force obtained in step 4.5 to calculate , and the measured value P of the rolling force of the optimal sample screened in step 3.2 is measured , and the correction coefficient of the rolling force model is calculated:

轧制力模型修正系数计算公式:

Figure BDA0003474779710000164
The calculation formula of the correction coefficient of the rolling force model:
Figure BDA0003474779710000164

以第1机架为例,轧制力模型修正系数为

Figure BDA0003474779710000165
Taking the first stand as an example, the correction coefficient of the rolling force model is
Figure BDA0003474779710000165

轧制力模型修正系数计算值如下:The calculation value of the correction coefficient of the rolling force model is as follows:

Figure BDA0003474779710000171
Figure BDA0003474779710000171

步骤5.3:计算轧制速度模型修正系数Step 5.3: Calculate the rolling speed model correction factor

使用步骤2.1得到的穿带过程中轧制速度测定值v实测-穿带过程(见表2),以及步骤3.2得到的穿带完成后的轧制过程的最优样本中的轧制速度测定值v(见表3),计算轧制速度模型修正系数

Figure BDA0003474779710000172
Use the measured value v of the rolling speed during the strip threading process obtained in step 2.1 - the strip threading process (see Table 2), and the measured value of the rolling speed in the optimal sample of the rolling process after the strip threading obtained in step 3.2 v (see Table 3), calculate the rolling speed model correction factor
Figure BDA0003474779710000172

轧制速度模型修正系数计算公式:

Figure BDA0003474779710000173
The calculation formula of the correction coefficient of the rolling speed model:
Figure BDA0003474779710000173

以第1机架为例,轧制速度模型修正系数为

Figure BDA0003474779710000174
Taking the first stand as an example, the correction coefficient of the rolling speed model is
Figure BDA0003474779710000174

轧制速度模型修正系数计算值如下:The calculation value of the correction coefficient of the rolling speed model is as follows:

Figure BDA0003474779710000175
Figure BDA0003474779710000175

步骤6:模型修正系数更新Step 6: Model Correction Coefficient Update

步骤6.1:根据修正系数旧值与修正系数计算值计算平滑系数,计算公式如下:Step 6.1: Calculate the smoothing coefficient according to the old value of the correction coefficient and the calculated value of the correction coefficient. The calculation formula is as follows:

Figure BDA0003474779710000176
Figure BDA0003474779710000176

式中:Δold为模型修正系数旧值,分别代表

Figure BDA0003474779710000177
Figure BDA0003474779710000178
Δ计算为步骤5.1-5.3计算得到的模型修正系数计算值,分别代表
Figure BDA0003474779710000179
Figure BDA00034747797100001710
Δmax为对应的模型修正系数最大值,分别代表
Figure BDA00034747797100001711
Figure BDA00034747797100001712
Δmin为对应的模型修正系数最小值,分别代表
Figure BDA00034747797100001713
Figure BDA00034747797100001714
对于三个修正系数,均有若Δ计算≥Δmax,则Δ计算替代并成为新的Δmax值;若Δ计算≤Δmin,则Δ计算替代并成为新的Δmin值;ω为比例系数;此处取ω=0.70;In the formula: Δ old is the old value of the model correction coefficient, representing
Figure BDA0003474779710000177
and
Figure BDA0003474779710000178
Δ is calculated as the calculated value of the model correction coefficient calculated in steps 5.1-5.3, representing
Figure BDA0003474779710000179
and
Figure BDA00034747797100001710
Δmax is the maximum value of the corresponding model correction coefficient, representing
Figure BDA00034747797100001711
and
Figure BDA00034747797100001712
Δmin is the minimum value of the corresponding model correction coefficient, representing
Figure BDA00034747797100001713
and
Figure BDA00034747797100001714
For the three correction coefficients, if Δcalculation≥Δmax , then Δcalculation is substituted and becomes the new Δmax value ; if Δcalculation≤Δmin , then Δcalculation is substituted and becomes the new Δmin value; ω is the proportional coefficient ; take ω=0.70 here;

本实施例中辊缝模型修正系数旧值

Figure BDA00034747797100001715
The old value of the correction coefficient of the roll gap model in this embodiment
Figure BDA00034747797100001715

本实施例中轧制力模型修正系数旧值

Figure BDA00034747797100001716
The old value of the correction coefficient of the rolling force model in this embodiment
Figure BDA00034747797100001716

本实施例中轧制速度模型修正系数旧值

Figure BDA00034747797100001717
The old value of the correction coefficient of the rolling speed model in this embodiment
Figure BDA00034747797100001717

以第1机架轧制力修正系数的平滑系数为例,

Figure BDA00034747797100001718
计算值为1.06,大于
Figure BDA00034747797100001719
因此计算得到的第1机架轧制力平滑系数为:
Figure BDA0003474779710000181
Taking the smoothing coefficient of the rolling force correction coefficient of the first stand as an example,
Figure BDA00034747797100001718
The calculated value is 1.06, which is greater than
Figure BDA00034747797100001719
Therefore, the calculated smoothing coefficient of the rolling force of the first stand is:
Figure BDA0003474779710000181

各机架辊缝、轧制力和轧制速度的平滑系数计算结果如下:The calculation results of the smoothing coefficients of the roll gap, rolling force and rolling speed of each stand are as follows:

参数名称parameter name 1机架1 rack 2机架2 racks 3机架3 racks 4机架4 racks 5机架5 racks 6机架6 racks 7机架7 racks 辊缝平滑系数Roll Gap Smoothing Coefficient 0.140.14 0.510.51 0.410.41 0.070.07 0.000.00 0.450.45 0.280.28 轧制力平滑系数Rolling force smoothing factor 0.380.38 0.410.41 0.350.35 0.210.21 0.070.07 0.070.07 0.000.00 轧制速度平滑系数Rolling speed smoothing factor 0.210.21 0.140.14 0.000.00 0.000.00 0.000.00 0.000.00 0.000.00

步骤6.2:模型修正系数新值计算Step 6.2: Calculate the new value of the model correction coefficient

采用平滑的方式对模型系数进行修正,平滑后的模型系数新值Δnew按下式计算:The model coefficients are corrected in a smoothing manner, and the new value Δnew of the smoothed model coefficients is calculated as follows:

Δnew=Δold+α·(Δ计算old) Δnew = Δold +α·( Δcalculation - Δold )

以第1机架轧制力修正系数新值计算为例:Take the calculation of the new value of the rolling force correction coefficient of the first stand as an example:

计算得到的修正系数新值为:α=1.00+0.38×(1.06-1.00)=1.02The new value of the calculated correction coefficient is: α=1.00+0.38×(1.06-1.00)=1.02

平滑后的模型系数新值计算结果如下表所示:The calculation results of the new values of the model coefficients after smoothing are shown in the following table:

Figure BDA0003474779710000182
Figure BDA0003474779710000182

步骤6.3:将计算得到的模型修正系数新值传递给模型设定系统,当下一块带钢进行轧制规程计算时,将辊缝修正系数与模型设定结果相加作为新的辊缝设定值,将轧制力修正系数与模型设定结果相乘作为新的轧制力设定值,将轧制速度轧制速度模型修正系数与模型设定结果相乘作为新的轧制速度设定值,供带钢生产使用。Step 6.3: Transfer the calculated new value of the model correction coefficient to the model setting system, and add the roll gap correction coefficient and the model setting result as the new roll gap setting value when the next strip is calculated for the rolling schedule , multiply the rolling force correction coefficient and the model setting result as the new rolling force setting value, and multiply the rolling speed rolling speed model correction coefficient and the model setting result as the new rolling speed setting value , for strip production.

步骤7:结束。Step 7: End.

图5为采用本实施例的方法,对辊缝、轧制力、轧制速度等参数的设定值进行修正前后,最末机架出口厚度与目标值的偏差测定结果。可以看出与修正之前相比,对于前一部分采样点,即对应带钢头部的采样点,修正之后的厚度偏差值大幅下降,高精度的头部厚度有助于快速进入后续的厚度自动控制(AGC)过程,从而保证带钢全长的厚度控制精度。Figure 5 shows the measurement results of the deviation between the thickness of the final stand outlet and the target value before and after the set values of the parameters such as roll gap, rolling force, and rolling speed are corrected using the method of this embodiment. It can be seen that compared with before the correction, for the first part of the sampling points, that is, the sampling points corresponding to the strip head, the thickness deviation value after the correction is greatly reduced, and the high-precision head thickness helps to quickly enter the subsequent automatic thickness control. (AGC) process, thereby ensuring the thickness control accuracy of the full length of the strip.

Claims (10)

1.一种基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,包括以下步骤:1. a steel rolling product quality optimization control method based on data evaluation, is characterized in that, comprises the following steps: 步骤1:确定穿带完成后的轧制过程的样本个数l,并确定每个样本中包含的目标采样点个数N;Step 1: Determine the number of samples l in the rolling process after the belt threading is completed, and determine the number of target sampling points N included in each sample; 步骤2:带钢依次通过各机架开始穿带过程,对生产过程中产生的实测数据进行采集和存储;Step 2: The strip starts to pass through each rack in turn, and the measured data generated in the production process is collected and stored; 步骤3:对已采集样本中的数据进行有效性评估和筛选;Step 3: Validate and screen the data in the collected samples; 步骤4:使用筛选出的最优样本的数据确定工艺参数计算值;Step 4: Use the data of the selected optimal sample to determine the calculated value of the process parameters; 步骤5:确定各机架辊缝模型、轧制力模型和轧制速度模型修正系数的计算值;Step 5: Determine the calculated values of the correction coefficients of the roll gap model, rolling force model and rolling speed model of each stand; 步骤6:对各机架的模型修正系数更新;Step 6: Update the model correction coefficient of each rack; 步骤7:结束。Step 7: End. 2.根据权利要求1所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,包括以下步骤:2. the steel rolling product quality optimization control method based on data evaluation according to claim 1, is characterized in that, comprises the following steps: 步骤1:确定穿带完成后的轧制过程的样本个数l,并确定每个样本中包含的目标采样点个数N;Step 1: Determine the number of samples l in the rolling process after the belt threading is completed, and determine the number of target sampling points N included in each sample; 步骤2:带钢依次通过各机架开始穿带过程,对生产过程中产生的实测数据进行采集和存储:Step 2: The strip starts to pass through each rack in turn, and the measured data generated in the production process is collected and stored: 步骤2.1:穿带过程中带钢头部数据的采集和存储:Step 2.1: Collection and storage of strip head data during belt threading: 轧制采用的机架数目为n;The number of stands used for rolling is n; 带钢头部通过第1机架,穿带过程开始,带钢沿轧制方向依次通过各机架和各机架处设置的测量仪表;The head of the strip passes through the first frame, the strip threading process begins, and the strip passes through each frame and the measuring instruments set at each frame in turn along the rolling direction; 按照固定的数据采样周期,依次对带钢头部通过测量仪表后的采样点进行计数,同时采集记录测量仪表测出的该采样点的实测数据;当通过测量仪表的采样点计数等于目标采样点数目N时,计算采集到的实测数据的平均值,并将实测数据的平均值进行存储作为测定值,作为穿带过程的实测数据样本;According to the fixed data sampling period, the sampling points after the strip head passes through the measuring instrument are counted in turn, and the measured data of the sampling point measured by the measuring instrument are collected and recorded at the same time; when the number of sampling points passing through the measuring instrument is equal to the target sampling point number When the target is N, calculate the average value of the collected measured data, and store the average value of the measured data as the measured value, as the measured data sample of the belt wearing process; 依次在各个机架处进行上述采集和存储过程,直至带钢头部通过最末测量仪表,并且在最末测量仪表处也完成上述采集和存储过程;Carry out the above-mentioned collection and storage process at each rack in turn, until the strip head passes through the last measuring instrument, and also complete the above-mentioned collection and storage process at the last measuring instrument; 穿带过程中,采集和存储的数据包括穿带过程各个机架处的轧制速度;During the belt threading process, the collected and stored data includes the rolling speed at each stand during the belt threading process; 步骤2.2:穿带完成后的轧制过程中带钢数据的采集和存储:Step 2.2: Collection and storage of strip data during rolling after strip threading: 带钢头部通过最末测量仪表后,开始对穿带完成后的轧制过程中带钢实测数据进行采集和存储:按照固定的数据采样周期,对通过测量仪表的采样点进行计数,同时采集记录测量仪表测出的该采样点的实测数据,当通过测量仪表的采样点计数等于目标采样点数目N时,分别计算当前样本内的所有实测数据的平均值,并将当前样本内实测数据平均值进行存储作为测定值,样本数目记为1;而后,采样点清零,重新开始计数,样本数目+1,直至完成l个样本的数据采集和存储过程;After the strip head passes through the last measuring instrument, it starts to collect and store the measured data of the strip during the rolling process after the strip is completed: according to the fixed data sampling period, the sampling points passing through the measuring instrument are counted and collected at the same time. Record the measured data of the sampling point measured by the measuring instrument. When the sampling point count through the measuring instrument is equal to the number of target sampling points N, calculate the average value of all the measured data in the current sample respectively, and average the measured data in the current sample. The value is stored as the measured value, and the number of samples is recorded as 1; then, the sampling point is cleared, and the counting is restarted, and the number of samples is +1 until the data collection and storage process of l samples is completed; 穿带完成后的轧制过程中,采集和存储的数据包括各个机架处的轧制速度、辊缝、轧制力和活套角度,以及最末机架出口的带钢温度、带钢厚度和带钢宽度;During the rolling process after strip threading, the data collected and stored include the rolling speed, roll gap, rolling force and looper angle at each stand, as well as the strip temperature and strip thickness at the exit of the last stand and strip width; 步骤3:对已采集样本中的数据进行有效性评估和筛选Step 3: Validation and screening of data in collected samples 步骤3.1:对穿带过程中和穿带完成后采集到的样本中的数据进行有效性评估Step 3.1: Evaluate the validity of the data from the samples collected during and after the threading process 若样本中的任一数据超出有效性区间范围,则转至步骤7;若样本中的所有数据均在有效性区间范围内,则转至步骤3.2;If any data in the sample is outside the validity interval, go to step 7; if all the data in the sample are within the validity interval, go to step 3.2; 步骤3.2:从步骤2.2中得到的l个样本数据中筛选最优样本,作为用于计算的样本;Step 3.2: Screen the optimal sample from the l sample data obtained in step 2.2 as a sample for calculation; 步骤4:使用筛选出的最优样本的数据确定工艺参数计算值Step 4: Use the data of the selected optimal sample to determine the calculated value of the process parameters 步骤4.1确定各个机架的出口厚度计算值Step 4.1 Determine the calculated value of the outlet thickness of each rack 根据步骤3.2得到的最优样本中的最末机架出口带钢厚度、各个机架处轧制速度和活套角度计算各机架出口厚度;对于i为1至n-1的第i机架的出口厚度计算值h计算,i为:Calculate the exit thickness of each stand according to the strip thickness at the exit of the last stand, the rolling speed at each stand and the looper angle in the optimal sample obtained in step 3.2; for the i-th stand where i is 1 to n-1 The calculated value of the outlet thickness h is calculated, and i is:
Figure FDA0003474779700000021
Figure FDA0003474779700000021
式中,fi为第i机架的前滑值,fn为最末机架的前滑值;hn、vn分别为最优样本中,最末机架出口的带钢厚度和轧制速度测定值,vi为最优样本中,第i机架处的轧制速度测定值;
Figure FDA0003474779700000022
为根据最优样本中第i机架处的活套角度θi测定值计算出的第i机架至第i+1机架间的带钢长度,lθs为根据目标活套角度θs计算出的第i机架至第i+1机架间的带钢目标长度;
In the formula, f i is the forward slip value of the i-th frame, f n is the forward slip value of the last frame; h n and v n are the strip thickness and rolling at the exit of the last frame in the optimal sample, respectively. Measured value of rolling speed, v i is the measured value of rolling speed at the i-th stand in the optimal sample;
Figure FDA0003474779700000022
is the strip length between the i-th frame and the i+1-th frame calculated according to the measured value of the looper angle θ i at the i-th frame in the optimal sample, and l θs is calculated according to the target looper angle θ s The target strip length from the i-th frame to the i+1-th frame;
最末机架出口厚度计算值h计算,n直接采用hnCalculate the calculated value h of the final frame outlet thickness , and h n is directly used for n ; 步骤4.2确定各个机架的出口温度计算值Step 4.2 Determine the calculated outlet temperature for each rack 第i机架的出口温度计算值T计算,i计算公式为:The calculated value T of the outlet temperature of the i-th rack is calculated, and the calculation formula of i is:
Figure FDA0003474779700000023
Figure FDA0003474779700000023
式中,Ti为步骤2.1的穿带过程中第i机架的轧制温度设定值,T实测为最优样本中最末机架出口的带钢温度测定值,T目标为最末机架出口带钢温度目标值;In the formula, T i is the set value of the rolling temperature of the i-th stand during the strip threading process in step 2.1, T is the measured value of the strip temperature at the exit of the last stand in the optimal sample, and T target is the final machine. The target value of the strip temperature at the outlet of the rack; 步骤4.3各机架出口宽度值确定Step 4.3 Determine the outlet width value of each rack 将最优样本中最末机架出口的带钢宽度测定值作为带钢通过各机架时的出口宽度值;Take the measured value of the strip width at the exit of the last rack in the optimal sample as the exit width value when the strip passes through each rack; 步骤4.4确定各机架辊缝计算值Step 4.4 Determine the calculated value of each rack roll gap 使用弹跳方程确定各机架处的辊缝计算值,各个机架处的辊缝计算值S计算计算公式如下:Use the bouncing equation to determine the calculated value of the roll gap at each frame, and the calculation formula of the calculated value S of the roll gap at each frame is as follows:
Figure FDA0003474779700000031
Figure FDA0003474779700000031
式中,h计算为所要计算的机架处的出口厚度计算值,由步骤4.1得到;P0为调零轧制力,M为轧机刚度;P实测为最优样本中所要计算的机架处的轧制力测定值;In the formula, h is calculated as the calculated value of the outlet thickness at the stand to be calculated, obtained from step 4.1; P 0 is the zero-adjusted rolling force, M is the stiffness of the rolling mill; P is measured as the calculated stand in the optimal sample. The measured value of rolling force; 步骤4.5确定各机架轧制力计算值Step 4.5 Determine the calculated value of rolling force of each stand 由步骤4.1至步骤4.3的各个机架的出口厚度计算值、各个机架的出口温度计算值、各机架出口宽度值,计算确定各机架轧制力计算值P计算From the calculated value of the outlet thickness of each stand, the calculated value of the outlet temperature of each stand, and the value of the exit width of each stand in steps 4.1 to 4.3, the calculated value P of the rolling force of each stand is calculated and determined; 步骤5:确定各机架辊缝模型、轧制力模型和轧制速度模型修正系数的计算值Step 5: Determine the calculated values of the correction coefficients of the roll gap model, rolling force model and rolling speed model of each stand 步骤5.1:计算辊缝模型修正系数Step 5.1: Calculate the Roll Gap Model Correction Factor 使用步骤4.4计算得到的辊缝计算值S计算,步骤3.2筛选得到的最优样本的辊缝测定值S实测,计算辊缝模型修正系数:Use the calculated value S of the roll gap calculated in step 4.4 to calculate , and the measured value S of the roll gap of the optimal sample screened in step 3.2 to be measured , and calculate the correction coefficient of the roll gap model: 辊缝模型修正系数计算公式:
Figure FDA0003474779700000032
The calculation formula of the correction coefficient of the roll gap model:
Figure FDA0003474779700000032
步骤5.2:计算轧制力模型修正系数Step 5.2: Calculate the rolling force model correction factor 使用步骤4.5得到的轧制力计算值P计算,步骤3.2筛选得到的最优样本的轧制力测定值P实测,计算轧制力模型修正系数:Use the calculated value P of the rolling force obtained in step 4.5 to calculate , and the measured value P of the rolling force of the optimal sample screened in step 3.2 is measured , and the correction coefficient of the rolling force model is calculated: 轧制力模型修正系数计算公式:
Figure FDA0003474779700000033
The calculation formula of the correction coefficient of the rolling force model:
Figure FDA0003474779700000033
步骤5.3:计算轧制速度模型修正系数Step 5.3: Calculate the rolling speed model correction factor 使用步骤2.1得到的穿带过程中轧制速度测定值v实测-穿带过程,以及步骤3.2得到的穿带完成后的轧制过程的最优样本中的轧制速度测定值v,计算轧制速度模型修正系数
Figure FDA0003474779700000034
Use the measured value v of the rolling speed during the strip threading process obtained in step 2.1 - the strip threading process , and the measured value v of the rolling speed in the optimal sample of the rolling process after the strip threading obtained in step 3.2, to calculate the rolling Speed model correction factor
Figure FDA0003474779700000034
轧制速度模型修正系数计算公式:
Figure FDA0003474779700000035
The calculation formula of the correction coefficient of the rolling speed model:
Figure FDA0003474779700000035
步骤6:各机架模型修正系数更新Step 6: Update the correction factor of each rack model 步骤6.1:根据修正系数旧值与修正系数计算值间的偏离程度计算平滑系数;计算公式如下:Step 6.1: Calculate the smoothing coefficient according to the degree of deviation between the old value of the correction coefficient and the calculated value of the correction coefficient; the calculation formula is as follows:
Figure FDA0003474779700000041
Figure FDA0003474779700000041
式中:Δold为模型修正系数旧值,包括辊缝、轧制力、轧制速度的模型修正系数旧值
Figure FDA0003474779700000042
Figure FDA0003474779700000043
Figure FDA0003474779700000044
Δ计算为步骤5.1-5.3计算得到的模型修正系数计算值,包括
Figure FDA0003474779700000045
Figure FDA0003474779700000046
Δmax为对应的模型修正系数最大值,包括辊缝、轧制力、轧制速度的模型修正系数最大值
Figure FDA0003474779700000047
Figure FDA0003474779700000048
Figure FDA0003474779700000049
Δmin为对应的模型修正系数最小值,包括辊缝、轧制力、轧制速度的模型修正系数最小值
Figure FDA00034747797000000410
Figure FDA00034747797000000411
In the formula: Δ old is the old value of the model correction coefficient, including the old value of the model correction coefficient of the roll gap, rolling force and rolling speed
Figure FDA0003474779700000042
Figure FDA0003474779700000043
and
Figure FDA0003474779700000044
Δ is calculated as the model correction coefficient calculated in steps 5.1-5.3, including
Figure FDA0003474779700000045
and
Figure FDA0003474779700000046
Δmax is the maximum value of the corresponding model correction coefficient, including the maximum value of the model correction coefficient of roll gap, rolling force and rolling speed
Figure FDA0003474779700000047
Figure FDA0003474779700000048
and
Figure FDA0003474779700000049
Δmin is the minimum value of the corresponding model correction coefficient, including the minimum value of the model correction coefficient of roll gap, rolling force and rolling speed
Figure FDA00034747797000000410
and
Figure FDA00034747797000000411
对于三个修正系数,均有:若Δ计算≥Δmax,则Δ计算替代并成为新的Δmax值;若Δ计算≤Δmin,则Δ计算替代并成为新的Δmin值;ω为比例系数;For the three correction coefficients, there are: if Δcalculation ≥Δmax , then Δcalculation replaces and becomes the new Δmax value; if Δcalculation ≤Δmin , then Δcalculation replaces and becomes the new Δmin value; ω is the ratio coefficient; 步骤6.2:模型修正系数新值计算Step 6.2: Calculate the new value of the model correction coefficient 采用平滑的方式对模型系数进行修正,平滑后的模型修正系数新值Δnew按下式计算:The model coefficients are corrected in a smoothing manner, and the new value of the smoothed model correction coefficient Δnew is calculated as follows: Δnew=Δold+α·(Δ计算old) Δnew = Δold +α·( Δcalculation - Δold ) 步骤6.3:将计算得到的模型修正系数新值传递给模型设定系统,当下一块带钢进行轧制规程计算时,将辊缝修正系数与模型设定结果相加作为新的辊缝设定值,将轧制力修正系数与模型设定结果相乘作为新的轧制力设定值,将轧制速度模型修正系数与模型设定结果相乘作为新的轧制速度设定值,供带钢生产使用;Step 6.3: Transfer the calculated new value of the model correction coefficient to the model setting system, and add the roll gap correction coefficient and the model setting result as the new roll gap setting value when the next strip is calculated for the rolling schedule , multiply the rolling force correction coefficient and the model setting result as the new rolling force setting value, multiply the rolling speed model correction coefficient and the model setting result as the new rolling speed setting value, supply the strip Steel production and use; 步骤7:结束。Step 7: End.
3.根据权利要求1所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,所述步骤1中,每个样本中包含的目标采样点个数N的确定方法如下:3. the steel rolling product quality optimization control method based on data assessment according to claim 1, is characterized in that, in described step 1, the determination method of the target sampling point number N that is included in each sample is as follows:
Figure FDA00034747797000000412
Figure FDA00034747797000000412
4.根据权利要求2所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,所述步骤2.1中,前m个采样点不进行采集,从第m+1个采样点开始采集和存储,样本共包含N个采样点,即第m+1到第m+N个采样点。4. The steel rolling product quality optimization control method based on data evaluation according to claim 2, characterized in that, in the step 2.1, the first m sampling points are not collected, and start from the m+1th sampling point For collection and storage, the sample contains N sampling points in total, that is, the m+1th to m+Nth sampling points. 5.根据权利要求2所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,所述步骤3.2中筛选最优样本的方法如下:5. the steel rolling product quality optimization control method based on data evaluation according to claim 2, is characterized in that, the method for screening optimal sample in described step 3.2 is as follows: 筛选样本数据波动程度最小的样本为最优样本,样本数据波动程度δaim的计算公式如下:The sample with the smallest fluctuation degree of sample data is selected as the optimal sample. The calculation formula of the fluctuation degree of sample data δ aim is as follows:
Figure FDA0003474779700000051
Figure FDA0003474779700000051
其中,选择p个数据种类作为用于计算样本数据波动程度的指标;δj为样本中第j项指标的波动程度,式中cact,j,k为第j项指标的第k个测定值,caim,j为第j项指标的目标值;Mj为各指标对应的数据测定值个数。Among them, select p data types as the index used to calculate the fluctuation degree of the sample data; δj is the fluctuation degree of the jth index in the sample, where c act,j,k is the kth measured value of the jth index , c aim,j is the target value of the jth index; M j is the number of data measurement values corresponding to each index.
6.根据权利要求5所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,选择最末机架出口的带钢温度、带钢厚度、带钢宽度、各个机架处的活套角度4个数据种类作为指标,最末机架出口的带钢温度、带钢厚度、带钢宽度对应的数据测定值个数为1个,各个机架处的活套角度对应的数据测定值个数为n-1个。6. The steel rolling product quality optimization control method based on data evaluation according to claim 5 is characterized in that, selecting the strip temperature, strip thickness, strip width and The four data types of the looper angle are used as indicators. The number of data measurement values corresponding to the strip temperature, strip thickness, and strip width at the exit of the last rack is one, and the data measurement corresponding to the looper angle at each rack The number of values is n-1. 7.根据权利要求2所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,所述步骤4.1中,根据两机架之间的活套角度θ计算两机架间的带钢长度lθ的方法如下:7. The steel rolling product quality optimization control method based on data evaluation according to claim 2, characterized in that, in the step 4.1, according to the looper angle θ between the two stands, the belt between the two stands is calculated The method for steel length l θ is as follows:
Figure FDA0003474779700000052
Figure FDA0003474779700000052
其中L1为前一机架到活套支点的水平距离,L2为活套支点至轧制平面的高度,L为两机架之间的水平距离,RL为活套臂长度,r为活套辊半径。Among them, L 1 is the horizontal distance from the previous stand to the fulcrum of the looper, L 2 is the height from the fulcrum of the looper to the rolling plane, L is the horizontal distance between the two stands, RL is the length of the looper arm, and r is the length of the looper arm. Looper roll radius.
8.根据权利要求4所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,所述对前m个采样点不进行采集,其中m为3-5。8 . The method for optimizing the quality of rolled steel products based on data evaluation according to claim 4 , wherein the first m sampling points are not collected, wherein m is 3-5. 9 . 9.根据权利要求2所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,所述穿带完成后的轧制过程的样本个数l为4-8。9 . The method for optimizing the quality of rolled steel products based on data evaluation according to claim 2 , wherein the number of samples 1 in the rolling process after the strip passing is completed is 4-8. 10 . 10.根据权利要求2所述的基于数据评估的钢铁轧制产品质量优化控制方法,其特征在于,所述步骤4.5中,确定各机架轧制力计算值P计算的方法如下:10. The steel rolling product quality optimization control method based on data evaluation according to claim 2, characterized in that, in the step 4.5, the method for determining the calculated value P of the rolling force of each stand is as follows: 根据西姆斯公式确定各机架轧制力计算值:Determine the calculated value of rolling force of each stand according to the Sims formula: P计算=1.15σslcQPw/1000P calculation = 1.15σ s l c Q P w/1000 式中:w—所要计算的机架的出口带钢宽度,mm,为步骤4.3确定的各机架出口宽度值;σs—所要计算的机架处的变形抗力,MPa;Where: w—the exit strip width of the rack to be calculated, mm, is the value of the exit width of each rack determined in step 4.3; σ s —the deformation resistance at the rack to be calculated, MPa;
Figure FDA0003474779700000061
Figure FDA0003474779700000061
a1~a6—回归系数,其值取决于钢种;a1~ a6 regression coefficient, its value depends on steel grade; T—所要计算的机架处的热力学温度,无量纲,
Figure FDA0003474779700000062
T计算为步骤4.2确定的所要计算机架的出口温度计算值;
T—the thermodynamic temperature at the rack to be calculated, dimensionless,
Figure FDA0003474779700000062
T is calculated as the calculated value of the outlet temperature of the desired computer rack determined in step 4.2;
Figure FDA0003474779700000063
—变形速率,s-1
Figure FDA0003474779700000064
其中v为最优样本中,所要计算的机架处的轧制速度测定值;
Figure FDA0003474779700000063
— deformation rate, s -1 ,
Figure FDA0003474779700000064
where v is the measured value of the rolling speed at the stand to be calculated in the optimal sample;
ε—工程应变,%,
Figure FDA0003474779700000065
Δh—所要计算的机架的压下量,mm,为所计算的机架入口厚度计算值与出口厚度计算值的差值;对于第1机架,入口厚度为带钢初始厚度h0,之后每一个机架的入口厚度计算值为上一个机架的出口厚度计算值,各机架出口厚度计算值按照步骤4.1获得;
ε—engineering strain, %,
Figure FDA0003474779700000065
Δh—the reduction amount of the frame to be calculated, mm, is the difference between the calculated value of the inlet thickness of the frame and the calculated value of the outlet thickness; for the first frame, the inlet thickness is the initial thickness h 0 of the strip, and then The calculated value of the inlet thickness of each rack is the calculated value of the outlet thickness of the previous rack, and the calculated value of the outlet thickness of each rack is obtained according to step 4.1;
σ0—在变形温度1000℃,
Figure FDA0003474779700000066
ε=0.4时的变形抗力,MPa;
σ 0 — at the deformation temperature of 1000°C,
Figure FDA0003474779700000066
Deformation resistance at ε=0.4, MPa;
lc—所要计算的机架考虑压扁后的接触弧长,mm,
Figure FDA0003474779700000067
l c —Contact arc length after crushing is considered for the frame to be calculated, mm,
Figure FDA0003474779700000067
R′—所要计算的机架的轧辊压扁半径,mm;
Figure FDA0003474779700000068
R'—roll flattening radius of the stand to be calculated, mm;
Figure FDA0003474779700000068
R—所要计算的机架的轧辊半径,mm;R—roll radius of the stand to be calculated, mm; QP—所要计算的机架的应力状态影响系数:Q P —Stress state influence coefficient of the rack to be calculated:
Figure FDA0003474779700000069
Figure FDA0003474779700000069
hm—带钢平均厚度,mm,为所要计算的机架处入口厚度计算值与出口厚度计算值的平均值;h m —the average thickness of the strip, mm, is the average value of the calculated value of the inlet thickness and the calculated value of the outlet thickness at the frame to be calculated; b0~b4为回归系数。b 0 to b 4 are regression coefficients.
CN202210052368.7A 2022-01-18 2022-01-18 A quality optimization control method for steel rolling products based on data evaluation Active CN114417236B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210052368.7A CN114417236B (en) 2022-01-18 2022-01-18 A quality optimization control method for steel rolling products based on data evaluation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210052368.7A CN114417236B (en) 2022-01-18 2022-01-18 A quality optimization control method for steel rolling products based on data evaluation

Publications (2)

Publication Number Publication Date
CN114417236A true CN114417236A (en) 2022-04-29
CN114417236B CN114417236B (en) 2024-09-27

Family

ID=81273906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210052368.7A Active CN114417236B (en) 2022-01-18 2022-01-18 A quality optimization control method for steel rolling products based on data evaluation

Country Status (1)

Country Link
CN (1) CN114417236B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119500788A (en) * 2024-10-29 2025-02-25 宿迁南钢金鑫轧钢有限公司 A production control method and system for small-scale steel continuous rolling process

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20040044325A (en) * 2002-11-21 2004-05-28 주식회사 포스코 Method for Controlling Thickness of Front Portion of Hot-Rolled Steel Sheet
CN107363105A (en) * 2017-07-20 2017-11-21 首钢集团有限公司 A kind of control method and device of hot rolled coil head thickness
CN112893484A (en) * 2021-01-18 2021-06-04 东北大学 Hot continuous rolling force self-setting method based on grey correlation extraction
CN113102516A (en) * 2021-03-05 2021-07-13 东北大学 Hot continuous rolling strip steel head width prediction method integrating rolling mechanism and deep learning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20040044325A (en) * 2002-11-21 2004-05-28 주식회사 포스코 Method for Controlling Thickness of Front Portion of Hot-Rolled Steel Sheet
CN107363105A (en) * 2017-07-20 2017-11-21 首钢集团有限公司 A kind of control method and device of hot rolled coil head thickness
CN112893484A (en) * 2021-01-18 2021-06-04 东北大学 Hot continuous rolling force self-setting method based on grey correlation extraction
CN113102516A (en) * 2021-03-05 2021-07-13 东北大学 Hot continuous rolling strip steel head width prediction method integrating rolling mechanism and deep learning

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
钱京学;毛一标;郭德福;: "热轧无取向硅钢首块钢头部厚度超差原因及控制对策", 河南冶金, no. 01, 15 February 2019 (2019-02-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119500788A (en) * 2024-10-29 2025-02-25 宿迁南钢金鑫轧钢有限公司 A production control method and system for small-scale steel continuous rolling process

Also Published As

Publication number Publication date
CN114417236B (en) 2024-09-27

Similar Documents

Publication Publication Date Title
CN103286142B (en) Dynamical roll gap compensation method during a kind of cold continuous rolling dynamic variable specification
CN111666653B (en) Online judging method for setting precision of strip steel finish rolling model
WO2023130666A1 (en) Strip steel plate convexity prediction method based on data-driving and mechanism model fusion
CN101869914A (en) Thickness control method of finish roller strip steel and device
CN107008758A (en) Cold-strip steel high accuracy plate shape surface roughness On-Line Control Method and system
CN104070070B (en) Control method for improving rolling force of precisely rolled strip steel and thickness precision through tension compensation
CN110052507B (en) A kind of strip steel coiling thickness control method and device
CN103506404A (en) Control method for roll gap in precision rolling process of strip steel
CN102744265A (en) Strip steel C warping control method
CN109226277B (en) Method for controlling flexible rolling thickness of non-ferrous metal
CN112287550B (en) Strip steel head thickness difference process parameter optimization method based on principal component analysis controller
CN111570533B (en) Prediction method for the influence of hot rolling coiling temperature and finish rolling temperature on surface deformation resistance
CN104801548A (en) Method for automatically improving per-second flow balance of strip steel in hot continuous rolling strip threading process
CN103240279B (en) The control device of continuous hot-rolling mill and the control method of continuous hot-rolling mill
CN102350441B (en) Control method of wall thickness of seamless steel tube
CN118023303A (en) Roll gap setting method, device and equipment in thin strip rolling production process
CN114417236B (en) A quality optimization control method for steel rolling products based on data evaluation
CN114091352B (en) Dynamic prediction method for convexity of hot continuous rolling outlet plate based on Elman neural network
CN106424156A (en) Cold tandem mill thickness control system and thickness control method of cold tandem mill process
CN110227725B (en) Roll gap self-learning deviation rectifying method and device
CN115121626A (en) Hot-rolled strip steel transient hot roll shape forecasting method based on error compensation
CN108817100B (en) A method of judging that certain breast roller causes strip exit thickness to fluctuate
CN112474826B (en) A Calculation Method of Thickness Distribution of Intermediate Billet in Hot Continuous Rolling
CN110614274B (en) Improved second flow thickness control method for single-stand reversible rolling mill
CN100409956C (en) Tube sizing rolling control method and sizing rolling control device

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