CN103191919B - 一种带钢轧制在线控制摩擦系数模型优化方法 - Google Patents
一种带钢轧制在线控制摩擦系数模型优化方法 Download PDFInfo
- Publication number
- CN103191919B CN103191919B CN201210002100.9A CN201210002100A CN103191919B CN 103191919 B CN103191919 B CN 103191919B CN 201210002100 A CN201210002100 A CN 201210002100A CN 103191919 B CN103191919 B CN 103191919B
- Authority
- CN
- China
- Prior art keywords
- coefficient
- centerdot
- friction
- value
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005096 rolling process Methods 0.000 title claims abstract description 55
- 229910000831 Steel Inorganic materials 0.000 title claims abstract description 42
- 239000010959 steel Substances 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005457 optimization Methods 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000003044 adaptive effect Effects 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 17
- 230000006978 adaptation Effects 0.000 claims description 16
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 7
- 238000009785 tube rolling Methods 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 5
- 238000012821 model calculation Methods 0.000 claims description 3
- 238000005097 cold rolling Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000008186 active pharmaceutical agent Substances 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 206010017389 Frotteurism Diseases 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004566 building material Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Landscapes
- Control Of Metal Rolling (AREA)
Abstract
本发明涉及一种带钢轧制在线控制摩擦系数模型优化方法,将带钢轧制过程的实测轧制力代入轧制力理论计算模型中。由于轧制力是摩擦系数的函数,因此通过变式处理得到含有以摩擦系数为自变量的非线性方程。除了将实测轧制力值代入轧制力理论计算模型外,还要将轧制过程中带钢的实际入口厚度、实际出口厚度、实际张力等数据同时代入轧制力计算模型。通过求解该非线性方程获得摩擦系数的实际计算值,以实际摩擦系数为基础,对摩擦系数理论计算模型进行自适应优化。通过该方法可以提高摩擦系数在线设定精度。
Description
技术领域
本发明涉及金属材料加工工程技术领域,特别涉及一种带钢轧制工艺与控制技术。
背景技术
冷轧带钢具有优良的力学性能、表面质量和高精度的厚度和板形几何精度,广泛应用于汽车、家电和建材等国民经济领域的制造部门。带钢冷轧生产过程中,需要进行高精度的控制,因此对轧制过程控制数学模型的精度要求非常高。
在冷轧轧制过程数学模型系统中,轧制力模型起到重要作用。冷轧带钢轧制过程中,摩擦系数是影响轧制力最重要的因素,直接决定了轧制力的大小。同时,摩擦系数模型和变形抗力模型的计算精度决定了轧制力模型的计算精度。因此,可以说摩擦系数对带钢冷轧过程的厚度和板形精度有较大的影响。为此,有必要对摩擦系数模型进行深入的研究。
带钢冷轧过程中需要适度的摩擦条件。如果摩擦系数太大则轧制负荷增加,造成轧制过程无法正常进行。如果摩擦系数太小,容易出现轧制打滑,出现轧机振动。因此,准确地计算出冷轧带钢摩擦系数对于产品质量提高具有重要作用。
带钢冷轧在线控制过程中,摩擦系数采用模型进行设定计算。但由于实际生产过程中,影响轧辊与带钢摩擦状态的因素很多,在线控制模型不可能一一考虑,并且无法通过在线仪表直接测量出摩擦系数。
发明内容
本发明的目的是提供一种带钢轧制在线控制摩擦系数模型优化方法,通过实际测量轧制力间接计算出摩擦系数的方法,根据计算出的实际摩擦系数对摩擦系数模型设定值进行自适应优化,该方法可有效提高轧制过程摩擦系数模型的计算精度。
为实现上述目的,本发明通过以下技术方案实现:
一种带钢轧制在线控制摩擦系数模型优化方法,该方法包括以下步骤:
1)在线检测带钢轧制过程的实测轧制力、实测带钢轧制入口厚度及出口厚度、实测带钢轧制过程的前后张力,将上述工艺参数代入轧制力计算模型中,构成以摩擦系数为未知量的非线性方程,求解该方程得到轧制过程摩擦系数的后计算值,将该摩擦系数的后计算值作为实际摩擦系数值;
所述的轧制力计算模型为:
式中:式中:μi back摩擦系数后计算值实际值;ri back压下率后计算实际值;Ri 'mdl不同机架号的轧辊压扁半径模型计算值;hi back带钢厚度后计算实际值;Fi measure轧制力实测值;kmi mdl带钢变形抗力模型计算值;ti back带钢张力后计算实际值;Δhi back带钢压下量后计算实际值;i机架号。
2)将上述的摩擦系数后计算值作为实际摩擦系数值,与摩擦系数模型计算值进行比较,计算出用于补偿摩擦系数模型计算值与实际值偏差的模型自适应学习系数,通过该模型自适应学习系数可以提高轧制过程摩擦系数模型的计算精度;
所述的摩擦系数模型为:
μ=μL·μv+μmin (5)
(1)轧制长度影响项
(2)速度影响项
μ 摩擦系数
L 工作辊轧制长度
v 出口带钢速度
μmin 摩擦系数基准值
aL~eL 轧制长度影响系数
av~fv 速度影响系数
v0 标准速度
C1 轧制长度影响项的自适应学习系数
C2 速度影响项的自适应学习系数
e模型常数值
通过公式(13)得出的摩擦系数的后计算值即实际摩擦系数值,再结合摩擦系数模型公式(5)、公式(6)、公式(7),可以得出轧制长度影响项的自适应学习系数C1及速度影响项的自适应学习系数C2;C1i、C2i为第i对机架的自适应学习系数;
极限检查:C1i lower≤C1i current≤C1i upper (15)
学习系数更新C1i next=(1-αC1i)·C1i current+αC1i·C1i before (16)
式中
极限检查:C2i lower≤C2i current≤C2i upper (18)
学习系数更新:C2i next=(1-αC2i)·C2i current+αC2i·C2i before (19)
式中C1i new,C2i new—下一次计算摩擦系数模型自适应学习系数
C1i current,C2i current—当前卷摩擦系数模型自适应学习系数计算值
C1i before,C2i before—上一次计算摩擦系数模型自适应学习系数
αC1i,αC2i—摩擦系数模型自适应学习系数平滑因子
i—机架号(1、2、3、4、5……)
3)将计算出的轧制长度影响项的自适应学习系数C1i及速度影响项的自适应学习系数C2i,代入摩擦系数模型公式(6)、公式(7),根据公式(5)计算出摩擦系数计算值。
与现有技术相比,本发明的有益效果是:
通过在线检测带钢轧制过程的实测轧制力、实测带钢轧制前后厚度、实测带钢轧制过程的前后张力等工艺参数,将上述工艺参数代入轧制力计算模型中,构成以摩擦系数为未知量的非线性方程,求解该方程得到轧制过程摩擦系数的后计算值,将该摩擦系数的后计算值作为实际摩擦系数值;根据计算出的实际摩擦系数对摩擦系数模型设定值进行自适应优化,该方法可有效提高轧制过程摩擦系数模型的计算精度。
附图说明
图1是五机架冷连轧机组的示意图;
图2是采用本发明提出的摩擦系数模型优化方法前后的计算值曲线图;
图3是采用本发明提出的摩擦系数模型优化方法前后的轧制力设定值、计算值、实测值曲线图。
具体实施方式
下面通过具体实施方式对本发明带钢轧制在线控制摩擦系数模型优化方法作详细的描述。
一种带钢轧制在线控制摩擦系数模型优化方法,通过在线测量带钢轧制过程的实测轧制力、实测带钢轧制入口厚度及出口厚度、实测带钢轧制过程的前后张力,将上述工艺参数代入轧制力计算模型中,构成以摩擦系数为未知量的非线性方程,求解该方程得到轧制过程摩擦系数的后计算值,将该摩擦系数的后计算值作为实际摩擦系数值;
工艺参数,间接计算出实际摩擦系数,根据计算出的实际摩擦系数对摩擦系数模型设定值进行自适应优化,该方法可有效提高轧制过程摩擦系数模型的计算精度。
一种带钢轧制在线控制摩擦系数模型优化方法,具体内容如下:
1、实际工艺参数处理
通过仪表实测出带钢轧制线上的相关工艺参数的实际值,i为机架号。
轧制力:实测轧制力Fi measure为操作侧测量值FWS和驱动侧测量值FDS的和,并去除弯辊力Pw的影响。
i为机架编号
总张力:机架间总张力T为工作侧张力TWS和驱动侧张力TDS的和。
单位张力t是在单位宽度、单位厚度的带钢上作用的张应力。
测厚仪测得的带钢厚度h:各测量点处的带钢厚度为厚度基准值href和偏差值hdev之和。W为带钢的宽度。
2、摩擦系数模型
μ=μL·μv+μmin (5)
(1)轧制长度影响项的摩擦系数μL
(2)速度影响项的摩擦系数μv
μ 摩擦系数
L 工作辊轧制长度
v 出口带钢速度
μmin 摩擦系数基准值
aL、bL、cL、dL、eL 轧制长度影响系数
av、bv、cv、dv、ev、fv 速度影响系数
v0 标准速度
C1 轧制长度影响项的摩擦系数模型学习系数
C2 速度影响项的摩擦系数模型学习系数
e 模型常数值
3、轧制力模型
F=Fp+Fe (8)
(1)塑性区轧制力
ξ=α·tin+β·tout
(2)弹性区轧制力
F 轧制力
Fp 塑性区轧制力
Fe 弹性区轧制力
Fe in 弹性压缩区轧制力
Fe out 弹性回复区轧制力
hin 入口厚度
hout 出口厚度
W 带钢宽度
tin 入口单位张力
tout 出口单位张力
km 平均变形抗力
μ 摩擦系数
R’ 轧辊压扁半径
QF 轧制力外摩擦影响系数
r 压下率
ν 泊松比
E 杨氏模量
α 入口张力影响系数
β 出口张力影响系数
4、变式处理
将轧制力计算模型进行变形处理,由公式(8)、(9)、(10)可得出:
求解上式可以得到摩擦系数,如下式:
5、摩擦系数后计算
利用实测轧制力Fi measure和其它实际工艺参数代入式(12)中,可求出摩擦系数的后计算实际值。
式中:μi back摩擦系数后计算值实际值;ri back压下率后计算实际值;R′mdl轧辊压扁半径模型计算值;hi back带钢厚度后计算实际值;Fi measure轧制力实测值;kmi mdl带钢变形抗力模型计算值;ti back带钢张力后计算实际值;Δhi back带钢压下量后计算实际值;i机架号。
6、摩擦系数模型的自适应优化
将摩擦系数后计算实际值μi back代入公式(5)、(6),可求得摩擦系数模型的自适应学习系数计算值,极限检查后用指数平滑法求得下次模型计算所需的自适应学习系数C1:
极限检查:C1i lower≤C1i current≤C1i upper (15)
学习系数更新C1i next=(1-αC1i)·C1i current+αC1i·C1i before (16)
将摩擦系数后计算实际值μi back代入公式(5)、(7),可求得摩擦系数模型的自适应学习系数计算值,极限检查后用指数平滑法求得下次模型计算所需的自适应学习系数C2:
式中
极限检查:C2i lower≤C2i current≤C2i upper (18)
学习系数更新:C2i next=(1-αC2i)·C2i current+αC2i·C2i before (19)
式中C1i new,C2i new—下一次计算摩擦系数模型自适应学习系数
C1i current,C2i current—当前卷摩擦系数模型自适应学习系数计算值
C1i before,C2i before—上一次计算摩擦系数模型自适应学习系数
αC1i,αC2i—摩擦系数模型自适应学习系数平滑因子
i—机架号(1~5)
其它符号的含义同前面所述。
实施例:以5机架冷连轧生产机组为例如图1,采用本专利提出的模型优化方法前后摩擦系数模型计算值如图2,采用摩擦系数优化计算前后值所计算出的轧制力设定计算值与实测值对比如图3。通过对比分析可知,采用优化后的摩擦系数模型设定值使轧制力设定精度提高。
Claims (1)
1.一种带钢轧制在线控制摩擦系数模型优化方法,其特征在于,该方法包括以下步骤:
1)在线检测带钢轧制过程的实测轧制力、实测带钢轧制入口厚度及出口厚度、实测带钢轧制过程的前后张力,将上述工艺参数代入轧制力计算模型中,构成以摩擦系数为未知量的非线性方程,求解该方程得到轧制过程摩擦系数的后计算值,将该摩擦系数的后计算值作为实际摩擦系数值;
所述的轧制力计算模型为:
式中:μi back摩擦系数后计算值实际值;ri back压下率后计算实际值;Ri 'mdl不同机架号的轧辊压扁半径模型计算值;hi back带钢厚度后计算实际值;Fi measure轧制力实测值;kmi mdl带钢变形抗力模型计算值;ti back带钢张力后计算实际值;Δhi back带钢压下量后计算实际值;i机架号;
2)将上述的摩擦系数后计算值作为实际摩擦系数值,与摩擦系数模型计算值进行比较,计算出用于补偿摩擦系数模型计算值与实际值偏差的模型自适应学习系数,通过该模型自适应学习系数可以提高轧制过程摩擦系数模型的计算精度;
所述的摩擦系数模型为:
μ=μL·μv+μmin (5)
(1)轧制长度影响项
(2)速度影响项
μ 摩擦系数
L 工作辊轧制长度
v 出口带钢速度
μmin 摩擦系数基准值
aL~eL 轧制长度影响系数
av~fv 速度影响系数
v0 标准速度
C1 轧制长度影响项的自适应学习系数
C2 速度影响项的自适应学习系数
e模型常数值
通过公式(13)得出的摩擦系数的后计算值即实际摩擦系数值,再结合摩擦系数模型公式(5)、公式(6)、公式(7),可以得出轧制长度影响项的自适应学习系数C1及速度影响项的自适应学习系数C2;C1i、C2i为第i对机架的自适应学习系数;
极限检查:C1i lower≤C1i current≤C1i upper (15)
学习系数更新C1i next=(1-αC1i)·C1i current+αC1i·C1i before (16)
式中
极限检查:C2i lower≤C2i current≤C2i upper (18)
学习系数更新:C2i next=(1-αC2i)·C2i current+αC2i·C2i before (19)
式中C1i new,C2i new—下一次计算摩擦系数模型自适应学习系数
C1i current,C2i current—当前卷摩擦系数模型自适应学习系数计算值
C1i before,C2i before—上一次计算摩擦系数模型自适应学习系数
αC1i,αC2i—摩擦系数模型自适应学习系数平滑因子
i—机架号(1、2、3、4、5……)
3)将计算出的轧制长度影响项的自适应学习系数C1i及速度影响项的自适应学习系数C2i,代入摩擦系数模型公式(6)、公式(7),根据公式(5)计算出摩擦系数计算值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210002100.9A CN103191919B (zh) | 2012-01-05 | 2012-01-05 | 一种带钢轧制在线控制摩擦系数模型优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210002100.9A CN103191919B (zh) | 2012-01-05 | 2012-01-05 | 一种带钢轧制在线控制摩擦系数模型优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103191919A CN103191919A (zh) | 2013-07-10 |
CN103191919B true CN103191919B (zh) | 2015-05-06 |
Family
ID=48714838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210002100.9A Active CN103191919B (zh) | 2012-01-05 | 2012-01-05 | 一种带钢轧制在线控制摩擦系数模型优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103191919B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526854B (zh) * | 2016-06-22 | 2020-08-21 | 中国科学院金属研究所 | 一种三辊行星轧机轧制力的计算方法 |
CN106734250B (zh) * | 2017-01-04 | 2019-02-19 | 河钢股份有限公司邯郸分公司 | 一种在线检测冷轧板表面摩擦系数的系统及其使用方法 |
CN112912185B (zh) * | 2018-11-13 | 2023-08-01 | 松下知识产权经营株式会社 | 辊压装置、以及控制装置 |
CN113434968B (zh) * | 2020-03-23 | 2023-11-07 | 上海梅山钢铁股份有限公司 | 一种基于钢板化学成分预测冷轧变形抗力的方法 |
CN113434994B (zh) * | 2020-03-23 | 2023-10-10 | 上海梅山钢铁股份有限公司 | 一种基于热轧原板热轧工艺参数预测冷轧变形抗力的方法 |
CN111570532B (zh) * | 2020-05-19 | 2021-05-25 | 燕山大学 | 热轧卷取温度与终轧温度对平整摩擦系数影响预测方法 |
CN111570534B (zh) * | 2020-06-18 | 2021-08-17 | 燕山大学 | 酸轧机组热轧特性对摩擦系数影响预测方法 |
CN116984385B (zh) * | 2023-09-26 | 2024-01-05 | 太原理工大学 | 冷连轧带钢由薄到厚规格轧制中轧制力确定方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1589986A (zh) * | 2003-08-29 | 2005-03-09 | 东北大学 | 金属板带轧制自动控制工艺参数的优化方法 |
CN1640573A (zh) * | 2004-01-16 | 2005-07-20 | 宝山钢铁股份有限公司 | 冷带钢连轧机轧辊工艺参数的优化控制方法 |
WO2007041179A2 (en) * | 2005-09-29 | 2007-04-12 | University Of Pittsburgh - Of The Commonwealth System Of Higher Education | System for controlling a rolling mill and method of controlling a rolling mill |
CN101025767A (zh) * | 2007-03-21 | 2007-08-29 | 燕山大学 | 冷连轧高速轧制过程中摩擦系数预报与设定方法 |
-
2012
- 2012-01-05 CN CN201210002100.9A patent/CN103191919B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1589986A (zh) * | 2003-08-29 | 2005-03-09 | 东北大学 | 金属板带轧制自动控制工艺参数的优化方法 |
CN1640573A (zh) * | 2004-01-16 | 2005-07-20 | 宝山钢铁股份有限公司 | 冷带钢连轧机轧辊工艺参数的优化控制方法 |
WO2007041179A2 (en) * | 2005-09-29 | 2007-04-12 | University Of Pittsburgh - Of The Commonwealth System Of Higher Education | System for controlling a rolling mill and method of controlling a rolling mill |
CN101025767A (zh) * | 2007-03-21 | 2007-08-29 | 燕山大学 | 冷连轧高速轧制过程中摩擦系数预报与设定方法 |
Non-Patent Citations (5)
Title |
---|
王伟,等.采用混合摩擦模型预报冷轧薄板轧制力.《钢铁研究学报》.CNKI,2000,第12卷(第1期),第10-13页. * |
王军生,等.冷连轧过程控制变形抗力模型的自适应学习.《东北大学学报(自然科学版)》.CNKI,2004,第25卷(第10期),第973-976页. * |
白金兰,等.提高冷轧过程控制轧制力模型的设定精度.《钢铁研究学报》.CNKI,2006,第18卷(第3期),第21-25页. * |
马庆龙,等.基于神经网络和自适应预报模型参数的平整轧制力模型.《塑性工程学报》.CNKI,2008,第15卷(第3期),第191-194页. * |
魏立新,等.基于自适应遗传算法的冷连轧轧制力模型自学习.《轧钢》.CNKI,2010,第27卷(第3期),第7-10页. * |
Also Published As
Publication number | Publication date |
---|---|
CN103191919A (zh) | 2013-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103191919B (zh) | 一种带钢轧制在线控制摩擦系数模型优化方法 | |
CN102189117B (zh) | 基于横向性能检测的冷轧带钢平直度前馈控制方法 | |
CN103286142B (zh) | 一种冷连轧动态变规格时的动态辊缝补偿方法 | |
CN102553923B (zh) | 平整机延伸率自适应补偿控制方法 | |
CN104841700B (zh) | 一种热轧带钢轧制力的优化设定方法 | |
CN102672003B (zh) | 不锈钢带钢拉矫机机组的工艺参数设定方法 | |
CN103418618B (zh) | 冷连轧机性能前馈厚度控制方法 | |
CN104338757B (zh) | 一种控制轧机启动轧制阶段轧制力的方法 | |
CN101716607B (zh) | Hc轧机非对称弯辊非对称横移板形控制方法 | |
CN103341503B (zh) | 一种自适应凸度变化热轧板形控制模型 | |
CN101134207A (zh) | 冷轧板形设定计算中热轧来料凸度的处理方法 | |
CN102744265B (zh) | 带钢c翘控制方法 | |
CN103962392B (zh) | 一种热连轧机精轧机组动态负荷控制方法 | |
CN101927267A (zh) | 一种精轧带钢辊缝的控制方法及装置 | |
CN103128107B (zh) | 一种热连轧粗轧短行程曲线参数的在线计算方法 | |
CN104324948A (zh) | 一种热连轧粗轧过程轧件宽度控制方法 | |
CN103567228B (zh) | 一种六辊轧机极薄带非常态轧制时板形与压靠预报方法 | |
CN104451118B (zh) | 一种适合于连续退火过程的带钢板形逐段演变预报方法 | |
CN102363159A (zh) | 一种单张板精密冷轧板厚测量系统的厚度控制方法 | |
CN106984651A (zh) | 一种提高轧件厚度控制精度的在线控制系统 | |
CN105251778A (zh) | 单锥度工作辊窜辊轧机边部减薄反馈控制方法 | |
CN111570534A (zh) | 酸轧机组基于大数据理论的热轧特性对摩擦系数影响预测方法 | |
CN106345817A (zh) | 工作辊长径比大且直径绝对值小的六辊轧机板形预报方法 | |
CN101927269A (zh) | 一种三辊轧机辊缝控制器调节方法 | |
CN101927268B (zh) | 冷连轧带钢厚度控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |