CN104289524B - 一种适合于四辊平整机组的复杂浪形控制方法 - Google Patents
一种适合于四辊平整机组的复杂浪形控制方法 Download PDFInfo
- Publication number
- CN104289524B CN104289524B CN201410467812.7A CN201410467812A CN104289524B CN 104289524 B CN104289524 B CN 104289524B CN 201410467812 A CN201410467812 A CN 201410467812A CN 104289524 B CN104289524 B CN 104289524B
- Authority
- CN
- China
- Prior art keywords
- amount
- roller
- inequality
- roll shifting
- shifting amount
- 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.)
- Expired - Fee Related
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B37/00—Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
- B21B37/28—Control of flatness or profile during rolling of strip, sheets or plates
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Control Of Metal Rolling (AREA)
Abstract
一种适合于四辊平整机组的复杂浪形控制方法,它主要包括以下由计算机执行的步骤:(1)收集设备参数;(2)收集产品的轧制工艺参数;(3)采用非对称手段对复杂浪形进行控制:①给定上下窜辊量过程变量、倾辊量的初始值以及弯辊力优化步长、窜辊量优化步长、倾辊量优化步长;②计算出上下窜辊量、倾辊量的最佳左、右弯辊力,③计算出左右弯辊力、倾辊量的最佳上、下窜辊量,④计算出左右弯辊力、上下窜辊量的最佳倾辊量;⑤输出最佳左弯辊力、右弯辊力、上窜辊量、下窜辊量、倾辊量,对平整机进行复杂浪形控制。本发明能够有效地控制非对称的局部浪形、复合浪等复杂浪形、提高成品板形质量,具有进一步的推广和使用价值。
Description
技术领域
本发明属于轧钢技术领域,特别涉及一种带钢板形的控制方法。
背景技术
近年来,随着大部分板带用户由低端转向高端,对带钢的板形质量提出了越来越高的要求,平整作为板带生产中最接近成品的一道工序,对提高成品带钢的板形质量起着举足轻重的影响。与冷轧及热轧工序不一样,平整过程压下量很小(一般都在0.5%至2%之间),因此轧制过程中变形热与摩擦热很少,轧辊热凸度几乎可以忽略不计,所以无法采用分段冷却的方法对平整机组板形进行精细控制。目前,对于四辊平整机组而言,其常规的板形控制手段主要包括倾辊、对称弯辊以及对称窜辊等三种,所能控制的板形缺陷主要有单边浪、双边浪以及中浪,而对非对称的局部浪形、复合浪等复杂浪形则无能为力。这样,如何利用平整机组对复杂浪形进行精细控制就成为现场技术攻关的重点。
参考文献:
[1]连家创,刘宏民.板厚板形控制[M].兵器工业出版社.1995
[2]白振华,韩林芳,李经洲,等.四辊轧机非常态轧制时板形模型的研究.机械工程学报,2012,48
发明内容
本发明的目的在于提供一种能够控制非对称的局部浪形、复合浪等复杂浪形、提高成品板形质量的适合于四辊平整机组的复杂浪形控制方法。本发明主要是利用轧机的基本数据采集***与板形机理模型,提出了一套采用非对称式结构来替代对称性结构,将板形控制手段从传统的对称弯辊(对称窜辊)、倾辊等二维增加到左弯辊、右弯辊、上窜辊、下窜辊以及倾辊等五维的板形控制思想,以实际板形与目标板形的差值最小为目标建立了一套适合于四辊平整机组的复杂浪形控制方法,
本发明主要包括以下由计算机执行的步骤:
(a)基本设备参数的收集,主要包括工作辊的辊径Dw、支撑辊的辊径Db、工作辊的辊身长度Lw、支撑辊的辊身长度Lb、工作辊弯辊缸距离lw、支撑辊压下螺丝中心距lb、轧机所允许的最大正弯辊力轧机所允许的最大负弯辊力轧机所允许的最大窜辊量δmax、轧机所允许的最大倾辊量ηmax;
(b)收集典型规格产品的基本轧制工艺参数,主要包括带材厚度H,宽度B,变形抗力k,轧制力P,轧制速度v,延伸率入口平均张力T0,出口平均张力T1,目标板形系数a0、a1、a2、a3、a4;
(c)四辊平整机组采用非对称手段控制复杂浪形的实现,主要包括以下由计算机执行的步骤:
c1)给定上窜辊量δs、下窜辊量δx、以及倾辊量η的初始值以及弯辊力优化步长ΔS、窜辊量优化步长Δδ、倾辊量优化步长Δη;
c2)计算出上窜辊量为δs、下窜辊量为δx、倾辊量为η时的最佳的左、右
弯辊力Szbest、Sybest,主要包括以下由计算机执行的步骤:
c2-1)定义最佳弯辊计算过程变量iS、jS,目标函数初始值F0;
c2-2)令iS=0,jS=0,取F0=1010;
c2-3)令左弯辊力为
c2-4)令右弯辊力为
c2-5)计算出上窜辊量为δs、下窜辊量为δx、倾辊量为η、左弯辊力为Sz、右弯辊力为Sy时带材的前张力分布值σ1i;
c2-6)计算出以I-Unit为单位板形分布值式中E为带材的弹性模量、v为泊松比;
c2-7)计算出板形控制过程函数g1(X)、g2(X)的值,其中g2(X)=max(|αi-βi|),αi为目标板形,i为带材在横向的条元号,n为带材在横向总的条元数,xi为带材在横向第i条元的坐标;
c2-8)计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值,其中A为加权系数;
c2-9)判断不等式F(X)<F0是否成立?如果不等式成立则令F0=F(X)、Szbest=Sz、Sybest=Sy,转入步骤c2-10);如果不等式不成立,则直接转入步骤c2-10);
c2-10)判断不等式是否成立?如果不等式成立则令js=js+1,转入步骤c2-4);如果不等式不成立则转入步骤c2-11);
c2-11)判断不等式是否成立?如果不等式成立则令is=is+1,转入步骤c2-3);如果不等式不成立则转入步骤c2-12);
c2-12)得到最佳的左、右弯辊力Szbest、Sybest;
c3)计算出左弯辊力为Szbest、右弯辊力为Sybest、倾辊量为η时的最佳上、下窜辊量δsbest、δxbest,主要包括以下可由计算机执行的步骤:
c3-1)定义最佳窜辊计算过程变量iδ、jδ,δ′s、δ′x;
c3-2)令iδ=0,jδ=0;
c3-3)令上窜辊量为δ′s=-δmax+iδΔδ;
c3-4)令下窜辊量为δ′x=-δmax+jδΔδ;
c3-5)计算出上窜辊量为δ′s、下窜辊量为δ′x、倾辊量为η、左弯辊力为Szbest、右弯辊力为Sybest时带材的前张力分布值σ1i;
c3-6)计算出以I-Unit为单位板形分布值板形控制过程函数g2(X)=max(|αi-βi|)的值;
c3-7)计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值;
c3-8)判断不等式F(X)<F0是否成立?如果不等式成立则令F0=F(X),δsbest=δ′s,δxbest=δ′x,转入步骤c3-9);如果不等式不成立则直接转入步骤c3-9);
c3-9)判断不等式是否成立?如果不等式成立则令jδ=jδ+1,转入步骤c3-4);如果不等式不成立则转入步骤c3-10);
c3-10)判断不等式是否成立?如果不等式成立则令iδ=iδ+1,转入步骤c3-3);如果不等式不成立则转入步骤c3-11);
c3-11)得到最佳的上、下窜辊量δsbest、δxbest;
c4)计算出左弯辊力为Szbest、右弯辊力为Sybest、上窜辊量为δsbest、下窜辊量为δxbest时的最佳倾辊量ηbest;
c4-1)定义最佳倾辊计算过程变量iη,η';
c4-2)令iη=0;
c4-3)令倾辊量为η'=-ηmax+iηΔη;
c4-4)计算出上窜辊量为δsbest、下窜辊量为δxbest、倾辊量为η'、左弯辊力为Szbest、右弯辊力为Sybest时带材的前张力分布值σ1i;
c4-5)计算出以I-Unit为单位板形分布值板形控制过程函数g2(X)=max(|αi-βi|)的值;
c4-6)计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值;
c4-7)判断不等式F(X)<F0是否成立?如果不等式成立则令F0=F(X),ηbest=η',转入步骤c4-8);如果不等式不成立则直接转入步骤c4-8);
c4-8)判断不等式是否成立?如果不等式成立则令iη=iη+1,然后转入步骤c4-3);如果不等式不成立则转入步骤c4-9);
c4-9)得到最佳的倾辊量ηbest;
c5)判断不等式 是否成立?如果不等式不成立则给令δs=δsbest、δx=δxbest、η=ηbest,转入步骤c2)处;如果不等式成立转入步骤c6);
c6)输出最佳左弯辊力Szbest、右弯辊力Sybest、上窜辊量δsbest、下窜辊量δxbest、倾辊量ηbest,对平整机进行复杂浪形控制。
本发明与现有技术相比具有如下优点和效果:
通过左弯辊、右弯辊、上窜辊、下窜辊以及倾辊等板形控制参数的优化设定,能够有效地控制非对称的局部浪形、复合浪等复杂浪形、提高成品板形质量,弥补了平整机组因无法采用精细冷却而带来的板形控制先天性不足问题,具有进一步的推广和使用价值。
附图说明
图1是本发明总计算流程图;
图2是本发明非对称手段控制复杂浪形的实现流程图;
图3是本发明窜辊力优化计算流程图;
图4是本发明窜辊量优化计算流程图;
图5是本发明倾辊量优化计算流程图;
图6是本发明实例1优化过程中的前张应力分布图;
图7是本发明实例1优化过程中的板形值分布图;
图8是本发明实例1优化后的板形值对比图;
图9是本发明实例2优化过程中的前张应力分布图;
图10是本发明实例2优化过程中的板形值分布图;
图11是本发明实例2优化后的板形值对比图。
具体实施方式
实施例1
一种适合于四辊平整机组的复杂浪形控制方法,如图1所示,首先,在步骤(a)中,收集基本设备参数,主要包括工作辊的辊径Dw=500mm、支撑辊的辊径Db=1100mm、工作辊的辊身长度Lw=1450mm、支撑辊的辊身长度Lb=1450mm、工作辊弯辊缸距离lw=2300mm、支撑辊压下螺丝中心距lb=2300mm、轧机所允许的最大正弯辊力轧机所允许的最大负弯辊力轧机所允许的最大窜辊量δmax=200mm、轧机所允许的最大倾辊量ηmax=150μm。
随后,在步骤(b)中,收集厚度为0.22mm、宽度为1220mm、钢种为SPCC的典型规格产品的基本轧制工艺参数,主要变形抗力k=300MPa,轧制力P=4200KN,轧制速度v=600m/min,延伸率入口平均张力T0=80MPa,出口平均张力T1=120MPa,目标板形系数a0=0、a1=0、a2=0、a3=0、a4=0。
随后,在步骤(c)中,四辊平整机组采用非对称手段对复杂浪形进行控制,计算过程如图2所示:
首先,在步骤(c1)中,给定上窜辊量过程变量δs=0mm、下窜辊量过程变量δx=0mm、以及倾辊量η=0μm做为初始值以及弯辊力优化步长ΔS=1KN、窜辊量优化步长Δδ=5mm、倾辊量优化步长Δη=2μm;
随后,在步骤(c2)中,计算出上窜辊量为δs=0mm、下窜辊量为δx=0mm、倾辊量为η=0μm时的最佳的左、右弯辊力Szbest、Sybest,主要包括以下可由计算机执行的步骤,如图3所示:
首先,在步骤(c2-1)中,定义最佳弯辊计算过程变量iS、jS,目标函数初始值F0;
随后,在步骤(c2-2)中,令iS=0,jS=0,取F0=1010;
随后,在步骤(c2-3)中,令左弯辊力为
随后,在步骤(c2-4)中,令右弯辊力为
随后,在步骤(c2-5)中,计算出上窜辊量为δs=0mm、下窜辊量为δx=0mm、倾辊量为η=0μm、左弯辊力为Sz=-500KN、右弯辊力为Sy=-500KN时带材的前张力分布值σ1i,其前张力分布值如附图6中曲线a所示;
随后,在步骤(c2-6)中,计算出以I-Unit为单位板形分布值式中E=2.1×105MPa为带材的弹性模量、v=0.3为泊松比,其板形值如图7中曲线a所示;
随后,在步骤(c2-7)中,计算出板形控制过程函数g1(X)、g2(X)的值,其中 、g2(X)=max(|αi-βi|)=max(|βi|)=8.90I,αi为目标板形,i为带材在横向的条元号,28为带材在横向总的条元数,xi为带材在横向第i条元的坐标
随后,在步骤(c2-8)中,由公式F(X)=A·g1(X)+(1-A)·g2(X)=0.5×5.406+0.5×8.908=7.16I,其中A=0.5为加权系数;
随后,在步骤(c2-9)中,判断条件F(X)=7.16<F0=1010成立,记录下最佳的目标函数F0=F(X)=7.16I和最佳的左、右弯辊力Szbest=Sz=-500KN,Sybest=Sy=-500KN然后转入步骤(c2-10);
随后,在步骤(c2-10)中,判断条件成立,则js=js+1=0+1=1,然后转入步骤(c2-4);
随后,在步骤(c2-11)中,判断条件成立,则is=is+1=0+1=1,然后转入步骤(c2-3);
随后,在步骤(c2-12)中,得到最佳的左、右弯辊力Szbest=-459KN、Sybest=-459KN;
随后,在步骤(c3)中,计算出左弯辊力为Szbest=-459KN、右弯辊力为Sybest-459KN、倾辊量为η=0μm时的最佳上、下窜辊量δsbest、δxbest,主要包括以下可由计算机执行的步骤,如图4所示;
随后,在步骤(c3-1)中,定义最佳窜辊计算过程变量iδ、jδ,δ′s、δ′x;
随后,在步骤(c3-2)中,令iδ=0,jδ=0;
随后,在步骤(c3-3)中,令上窜辊量为δ′s=-δmax+iδΔδ=-200+0×5=-200mm;
随后,在步骤(c3-4)中,令下窜辊量为δ′x=-δmax+jδΔδ=-200+0×5=-200mm;
随后,在步骤(c3-5)中,计算出上窜辊量为δs=-200mm、下窜辊量为δx=-200mm、倾辊量为η=0μm、左弯辊力为Sz=-459KN、右弯辊力为Sy=-459KN时带材的前张力分布值σ1i,其前张力分布值如附图6中曲线b所示;
随后,在步骤(c3-6)中,计算出以I-Unit为单位板形分布值其板形值如图7中曲线b所示,板形控制过程函数 g2(X)=max(|αi-βi|)=max(|βi|)=41.72I;
随后,在步骤(c3-7)中,计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值,F(X)=0.5×16.219+0.5×41.726=28.97I;
随后,在步骤(c3-8)中,判断条件F(X)=28.97<F0=6.28不成立,直接转入步骤(c3-9);
随后,在步骤(c3-9)中,判断条件成立,则jδ=jδ+1=0+1=1,然后转入步骤(c3-4);
随后,在步骤(c3-10)中,判断条件成立,则iδ=iδ+1=0+1=1,然后转入步骤(c3-3);
随后,在步骤(c3-11)中,得到最佳的上、下窜辊量δsbest=90mm、δxbest=-95mm;
随后,在步骤(c4)中,计算出左弯辊力为Szbest=-459KN、右弯辊力为Sybest=-459KN、上窜辊量为δsbest=90mm、下窜辊量为δxbest=-95mm时的最佳倾辊量ηbest,如图5所示;
随后,在步骤(c4-1)中,定义最佳倾辊计算过程变量iη,η'代;
随后,在步骤(c4-2)中,令iη=0;
随后,在步骤(c4-3)中,令倾辊量为η'=-ηmax+iηΔη=-150+0×2=-150μm;
随后,在步骤(c4-4)中,计算出上窜辊量为δsbest=90mm、下窜辊量为δxbest=-95mm、倾辊量为η'=-150μm、左弯辊力为Szbest=-459KN、右弯辊力为Sybest=-459KN时带材的前张力分布值σ1i,其前张力分布值如图6中曲线c所示;
随后,在步骤(c4-5)中,计算出以I-Unit为单位板形分布值其板形值如图7中曲线c所示,板形控制过程函数 g2(X)=max(|αi-βi|)=max(|βi|)=24.94I;
随后,在步骤(c4-6)中,计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值F(X)=0.5×10.887+0.5×24.948=17.92I;
随后,在步骤(c4-7)中,判断条件F(X)=17.92<F0=5.76不成立,直接转入步骤(c4-8);
随后,在步骤(c4-8)中,判断条件成立,则iη=iη+1=0+1=1,然后转入步骤(c4-3);
随后,在步骤(c4-9)中,得到当前最佳的倾辊量ηbest=0μm;
随后,在步骤(c5)中,判断条件(δ′s-δsbest)2+(δ′x-δxbest)+(η'-ηbest)2=(0-90)2+(0-(-95))2+(0-0)2=17125<ε0=0.01不成立则给定δs=δsbest=90mm、δx=δxbest=-95mm、η=ηbest=0μm作为初始值转入步骤(c2)处。
最后,在步骤(c6)中,输出最佳左弯辊力为Szbest=236KN、最佳右弯辊力为Sybest=-500KN、最佳上窜辊量为δsbest=125mm、最佳下窜辊量为δxbest=-40mm、最佳倾辊量为ηbest=112μm对平整机进行复杂浪形控制,如图8所示,图中曲线c为最佳的板形目标函数值F(X)=3.69I;图中曲线b为常规手段控制的最佳的板形目标函数值F(X)=5.84I,与常规控制手段相比,本实施例的目标函数减小36.8%;图中曲线a为不加控制手段的板形分布,板形最差。
实施例2
首先,在步骤(a)中收集基本设备参数,主要包括工作辊的辊径Dw=500mm、支撑辊的辊径Db=1100mm、工作辊的辊身长度Lw=1450mm、支撑辊的辊身长度Lb=1450mm、工作辊弯辊缸距离lw=2300mm、支撑辊压下螺丝中心距lb=2300mm、轧机所允许的最大正弯辊力轧机所允许的最大负弯辊力轧机所允许的最大窜辊量δmax=200mm、轧机所允许的最大倾辊量ηmax=150μm。
随后,在步骤(b)中收集厚度为0.22mm、宽度为1220mm、钢种为SPCD的典型规格产品的基本轧制工艺参数,主要包括带材变形抗力k=300MPa,轧制力P=4200KN,轧制速度v=600m/min,延伸率入口平均张力T0=80MPa,出口平均张力T1=120MPa,目标板形系数a0=0、a1=0、a2=0、a3=0、a4=0。
随后,在步骤(c)中,四辊平整机组采用非对称手段对复杂浪形进行控制:
首先,在步骤(c1)中,给定上窜辊量过程变量δs=0mm、下窜辊量过程变量δx=0mm、以及倾辊量η=0μm做为初始值以及弯辊力优化步长ΔS=1KN、窜辊量优化步长Δδ=5mm、倾辊量优化步长Δη=2μm;
随后,在步骤(c2)中,计算出上窜辊量为δs=0mm、下窜辊量为δx=0mm、倾辊量为η=0μm时的最佳的左、右弯辊力Szbest、Sybest,主要包括以下可由计算机执行的步骤:
首先,在步骤(c2-1)中,定义最佳弯辊计算过程变量iS、jS,目标函数初始值F0;
随后,在步骤(c2-2)中,令iS=0,jS=0,取F0=1010;
随后,在步骤(c2-3)中,令左弯辊力为
随后,在步骤(c2-4)中,令右弯辊力为
随后,在步骤(c2-5)中,计算出上窜辊量为δs=0mm、下窜辊量为δx=0mm、倾辊量为η=0μm、左弯辊力为Sz=-500KN、右弯辊力为Sy=-500KN时带材的前张力分布值σ1i,其前张应力分布如图9中曲线a所示;
随后,在步骤(c2-6)中,计算出以I-Unit为单位板形分布值式中E=2.1×105MPa为带材的弹性模量、v=0.3为泊松比,其板形分布值如图10中曲线a所示;
随后,在步骤(c2-7)中,计算出板形控制过程函数g1(X)、g2(X)的值,其中 g2(X)=max(|αi-βi|)=max(|βi|)=12.68I,αi为目标板形,i为带材在横向的条元号,28为带材在横向总的条元数,xi为带材在横向第i条元的坐标
随后,在步骤(c2-8)中,由公式F(X)=A·g1(X)+(1-A)·g2(X)=0.5×5.406+0.5×8.908=9.20I,其中A=0.5为加权系数;
随后,在步骤(c2-9)中,判断条件F(X)=9.20<F0=1010成立,记录下最佳的目标函数F0=F(X)=9.20I和最佳的左、右弯辊力Szbest=Sz=-500KN,Sybest=Sy=-500KN然后转入步骤(c2-10);
随后,在步骤(c2-10)中,判断条件 成立,则js=js+1=0+1=1,然后转入步骤(c2-4);
随后,在步骤(c2-11)中,判断条件 成立,则is=is+1=0+1=1,然后转入步骤(c2-3);
随后,在步骤(c2-12)中,得到最佳的左、右弯辊力Szbest=-100KN、Sybest=95KN;
随后,在步骤(c3)中,计算出左弯辊力为Szbest=-100KN、右弯辊力为Sybest=95KN、倾辊量为η=0μm时的最佳上、下窜辊量δsbest、δxbest,主要包括以下可由计算机执行的步骤;
随后,在步骤(c3-1)中,定义最佳窜辊计算过程变量iδ、jδ,δ′s、δ′x;
随后,在步骤(c3-2)中,令iδ=0,jδ=0;
随后,在步骤(c3-3)中,令上窜辊量为δ′s=-δmax+iδΔδ=-200+0×5=-200mm;
随后,在步骤(c3-4)中,令下窜辊量为δ′x=-δmax+jδΔδ=-200+0×5=-200mm;
随后,在步骤(c3-5)中,计算出上窜辊量为δs=-200mm、下窜辊量为δx=-200mm、倾辊量为η=0μm、左弯辊力为Sz=-100KN、右弯辊力为Sy=95KN时带材的前张力分布值σ1i,其前张力分布如图9中曲线b所示;
随后,在步骤(c3-6)中,计算出以I-Unit为单位板形分布值其板形分布值如图10中曲线b所示,板形控制过程函数 g2(X)=max(|αi-βi|)=max(|βi|)=47.31I;
随后,在步骤(c3-7)中,计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值,F(X)=0.5×24.895+0.5×47.31=36.10I;
随后,在步骤(c3-8)中,判断条件F(X)=28.97<F0=8.18不成立,直接转入步骤(c3-9);
随后,在步骤(c3-9)中,判断条件成立,则jδ=jδ+1=0+1=1,然后转入步骤(c3-4);
随后,在步骤(c3-10)中,判断条件成立,则iδ=iδ+1=0+1=1,然后转入步骤(c3-3);
随后,在步骤(c3-11)中,得到最佳的上、下窜辊量δsbest=30mm、δxbest=-60mm;
随后,在步骤(c4)中,计算出左弯辊力为Szbest=-100KN、右弯辊力为Sybest=95KN、上窜辊量为δsbest=30mm、下窜辊量为δxbest=-60mm时的最佳倾辊量ηbest;
随后,在步骤(c4-1)中,定义最佳倾辊计算过程变量iη,η'代;
随后,在步骤(c4-2)中,令iη=0;
随后,在步骤(c4-3)中,令倾辊量为η'=-ηmax+iηΔη=-150+0×2=-150μm;
随后,在步骤(c4-4)中,计算出上窜辊量为δsbest=30mm、下窜辊量为δxbest=-60mm、倾辊量为η'=-150μm、左弯辊力为Szbest=-100KN、右弯辊力为Sybest=95KN时带材的前张力分布值σ1i,其前张力分布值如图9中c所示;
随后,在步骤(c4-5)中,计算出以I-Unit为单位板形分布值其板形分布值如图10中c所示,板形控制过程函数 g2(X)=max(|αi-βi|)=max(|βi|)=47.31I;
随后,在步骤(c4-6)中,计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值F(X)=0.5×24.895+0.5×47.31=36.10I;
随后,在步骤(c4-7)中,判断条件F(X)=36.10<F0=7.66不成立,直接转入步骤(c4-8);
随后,在步骤(c4-8)中,判断条件成立,则iη=iη+1=0+1=1,然后转入步骤(c4-3);
随后,在步骤(4-9)中,得到当前最佳的倾辊量ηbest=10μm;
随后,在步骤(c5)中,判断条件(δ′s-δsbest)2+(δ′x-δxbest)+(η'-ηbest)2=(0-30)2+(0-(-60))2+(0-10)2=4600<ε0=0.01不成立则给定δs=δsbest=30mm、δx=δxbest=-60mm、η=ηbest=10μm作为初始值转入步骤(c2)处。
最后,在步骤(c6)中,输出最佳左弯辊力为Szbest=196KN、最佳右弯辊力为Sybest=136KN、最佳上窜辊量为δsbest=50mm、最佳下窜辊量为δxbest=10mm、最佳倾辊量为ηbest=42μm,对平整机进行复杂浪形控制,如附图11所示,图中曲线c为最佳的板形目标函数值F(X)=1.84I;图中曲线b为常规手段控制的非对称控制后最佳的板形目标函数值F(X)=2.66I,与常规控制手段相比,目标函数减小30.8%;图中曲线a为不加控制手段的板形分布图,板形最差。
Claims (1)
1.一种适合于四辊平整机组的复杂浪形控制方法,其特征在于:它主要包括以下由计算机执行的步骤:
(a)基本设备参数的收集,主要包括工作辊的辊径Dw、支撑辊的辊径Db、工作辊的辊身长度Lw、支撑辊的辊身长度Lb、工作辊弯辊缸距离lw、支撑辊压下螺丝中心距lb、轧机所允许的最大正弯辊力轧机所允许的最大负弯辊力轧机所允许的最大窜辊量δmax、轧机所允许的最大倾辊量ηmax;
(b)收集典型规格产品的基本轧制工艺参数,主要包括带材厚度H,宽度B,变形抗力k,轧制力P,轧制速度v,延伸率入口平均张力T0,出口平均张力T1,目标板形系数a0、a1、a2、a3、a4;
(c)四辊平整机组采用非对称手段控制复杂浪形的实现,主要包括以下由计算机执行的步骤:
c1)给定上窜辊量δs、下窜辊量δx、以及倾辊量η的初始值以及弯辊力优化步长ΔS、窜辊量优化步长Δδ、倾辊量优化步长Δη;
c2)计算出上窜辊量为δs、下窜辊量为δx、倾辊量为η时的最佳的左、右弯辊力Szbest、Sybest,主要包括以下由计算机执行的步骤:
c2-1)定义最佳弯辊计算过程变量iS、jS,目标函数初始值F0;
c2-2)令iS=0,jS=0,取F0=1010;
c2-3)令左弯辊力为
c2-4)令右弯辊力为
c2-5)计算出上窜辊量为δs、下窜辊量为δx、倾辊量为η、左弯辊力为Sz、右弯辊力为Sy时带材的前张力分布值σ1i;
c2-6)计算出以I-Unit为单位板形分布值式中E为带材的弹性模量、v为泊松比;
c2-7)计算出板形控制过程函数g1(X)、g2(X)的值,其中g2(X)=max(|αi-βi|),αi为目标板形,i为带材在横向的条元号,n为带材在横向总的条元数,xi为带材在横向第i条元的坐标;
c2-8)计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值,其中A为加权系数;
c2-9)判断不等式F(X)<F0是否成立?如果不等式成立则令F0=F(X)、Szbest=Sz、Sybest=Sy,转入步骤c2-10);如果不等式不成立,则直接转入步骤c2-10);
c2-10)判断不等式是否成立?如果不等式成立则令js=js+1,转入步骤c2-4);如果不等式不成立则转入步骤c2-11);
c2-11)判断不等式是否成立?如果不等式成立则令is=is+1,转入步骤c2-3);如果不等式不成立则转入步骤c2-12);
c2-12)得到最佳的左、右弯辊力Szbest、Sybest;
c3)计算出左弯辊力为Szbest、右弯辊力为Sybest、倾辊量为η时的最佳上、下窜辊量δsbest、δxbest,主要包括以下可由计算机执行的步骤:
c3-1)定义最佳窜辊计算过程变量iδ、jδ,δ′s、δ'x;
c3-2)令iδ=0,jδ=0;
c3-3)令上窜辊量为δ′s=-δmax+iδΔδ;
c3-4)令下窜辊量为δ'x=-δmax+jδΔδ;
c3-5)计算出上窜辊量为δ′s、下窜辊量为δ'x、倾辊量为η、左弯辊力为Szbest、右弯辊力为Sybest时带材的前张力分布值σ1i;
c3-6)计算出以I-Unit为单位板形分布值板形控制过程函数g2(X)=max(|αi-βi|)的值;
c3-7)计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值;
c3-8)判断不等式F(X)<F0是否成立?如果不等式成立则令F0=F(X),δsbest=δ′s,δxbest=δ′x,转入步骤c3-9);如果不等式不成立则直接转入步骤c3-9);
c3-9)判断不等式是否成立?如果不等式成立则令jδ=jδ+1,转入步骤c3-4);如果不等式不成立则转入步骤c3-10);
c3-10)判断不等式是否成立?如果不等式成立则令iδ=iδ+1,转入步骤c3-3);如果不等式不成立则转入步骤c3-11);
c3-11)得到最佳的上、下窜辊量δsbest、δxbest;
c4)计算出左弯辊力为Szbest、右弯辊力为Sybest、上窜辊量为δsbest、下窜辊量为δxbest时的最佳倾辊量ηbest;
c4-1)定义最佳倾辊计算过程变量iη,η';
c4-2)令iη=0;
c4-3)令倾辊量为η'=-ηmax+iηΔη;
c4-4)计算出上窜辊量为δsbest、下窜辊量为δxbest、倾辊量为η'、左弯辊力为Szbest、右弯辊力为Sybest时带材的前张力分布值σ1i;
c4-5)计算出以I-Unit为单位板形分布值板形控制过程函数g2(X)=max(|αi-βi|)的值;
c4-6)计算出目标函数F(X)=A·g1(X)+(1-A)·g2(X)的值;
c4-7)判断不等式F(X)<F0是否成立?如果不等式成立则令F0=F(X),ηbest=η',转入步骤c4-8);如果不等式不成立则直接转入步骤c4-8);
c4-8)判断不等式是否成立?如果不等式成立则令iη=iη+1,然后转入步骤c4-3);如果不等式不成立则转入步骤c4-9);
c4-9)得到最佳的倾辊量ηbest;
c5)判断不等式 是否成立?如果不等式不成立则给令δs=δsbest、δx=δxbest、η=ηbest,转入步骤c2)处;如果不等式成立转入步骤c6);
c6)输出最佳左弯辊力Szbest、右弯辊力Sybest、上窜辊量δsbest、下窜辊量δxbest、倾辊量ηbest对平整机进行复杂浪形控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410467812.7A CN104289524B (zh) | 2014-09-15 | 2014-09-15 | 一种适合于四辊平整机组的复杂浪形控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410467812.7A CN104289524B (zh) | 2014-09-15 | 2014-09-15 | 一种适合于四辊平整机组的复杂浪形控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104289524A CN104289524A (zh) | 2015-01-21 |
CN104289524B true CN104289524B (zh) | 2016-03-02 |
Family
ID=52309603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410467812.7A Expired - Fee Related CN104289524B (zh) | 2014-09-15 | 2014-09-15 | 一种适合于四辊平整机组的复杂浪形控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104289524B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105127211B (zh) * | 2015-09-24 | 2017-03-22 | 胡小东 | 一种金属板带横向仿形矫形方法及装置 |
CN105436211B (zh) * | 2015-12-17 | 2017-07-14 | 太原重工股份有限公司 | 单机架六辊可逆轧机组板形闭环反馈控制方法 |
CN105598183B (zh) * | 2016-01-14 | 2017-05-17 | 北京科技大学 | 兼顾浪形和断面的热轧高次曲线工作辊窜辊策略控制 |
CN105975771B (zh) * | 2016-04-29 | 2018-12-07 | 燕山大学 | 一种四辊轧机工作辊与支撑辊凸度转换关系的计算方法 |
CN107774719B (zh) * | 2016-08-29 | 2020-06-23 | 宝山钢铁股份有限公司 | 20辊轧机中间辊倒角位置动态模型控制方法 |
CN107685075A (zh) * | 2017-07-10 | 2018-02-13 | 首钢京唐钢铁联合有限责任公司 | 应用于平整机预防浪形发生的方法 |
CN111215454B (zh) * | 2018-11-26 | 2021-07-16 | 宝山钢铁股份有限公司 | 一种热轧带钢平整板形控制方法 |
CN112974545B (zh) * | 2019-12-12 | 2022-09-06 | 上海梅山钢铁股份有限公司 | 一种极薄规格t5料s形中浪缺陷的防控方法 |
CN113399472A (zh) * | 2021-08-04 | 2021-09-17 | 宝钢湛江钢铁有限公司 | 一种五机架六辊冷连轧机组板形模糊优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3934438A (en) * | 1973-05-09 | 1976-01-27 | Nippon Kokan Kabushiki Kaisha | Method of long-edge shape control for tandem rolling mill |
CN101301659A (zh) * | 2008-03-15 | 2008-11-12 | 燕山大学 | 双ucm平整机组基于机理模型的板形参数在线设定方法 |
CN102407239A (zh) * | 2010-09-21 | 2012-04-11 | 上海宝钢工业检测公司 | 一种过平整机轧辊边浪产生的判断方法 |
CN103990651A (zh) * | 2013-02-19 | 2014-08-20 | 宝山钢铁股份有限公司 | 一种针对薄宽带钢高次边浪的控制装置及其控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07115055B2 (ja) * | 1990-06-27 | 1995-12-13 | 株式会社日立製作所 | 多段圧延機の形状制御方法及び多段圧延機及び多段圧延機の形状制御装置及び多段圧延機の圧延方法 |
JP3004456B2 (ja) * | 1992-03-30 | 2000-01-31 | 三菱重工業株式会社 | クロスロール圧延機の板形状制御装置 |
-
2014
- 2014-09-15 CN CN201410467812.7A patent/CN104289524B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3934438A (en) * | 1973-05-09 | 1976-01-27 | Nippon Kokan Kabushiki Kaisha | Method of long-edge shape control for tandem rolling mill |
CN101301659A (zh) * | 2008-03-15 | 2008-11-12 | 燕山大学 | 双ucm平整机组基于机理模型的板形参数在线设定方法 |
CN102407239A (zh) * | 2010-09-21 | 2012-04-11 | 上海宝钢工业检测公司 | 一种过平整机轧辊边浪产生的判断方法 |
CN103990651A (zh) * | 2013-02-19 | 2014-08-20 | 宝山钢铁股份有限公司 | 一种针对薄宽带钢高次边浪的控制装置及其控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104289524A (zh) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104289524B (zh) | 一种适合于四辊平整机组的复杂浪形控制方法 | |
CN104275352B (zh) | 一种带材冷轧机跑偏与板形自动控制方法 | |
CN101513647B (zh) | 二次冷轧机组生产带材的平整方法 | |
CN106391725B (zh) | 一种适用于冷轧过程轧制压力随轧制速度变化预报方法 | |
CN103567229B (zh) | 一种针对六辊轧机的弯辊力组合板形控制方法 | |
CN104324948B (zh) | 一种热连轧粗轧过程轧件宽度控制方法 | |
CN101507978B (zh) | 冷连轧机基于机理与工况相结合的工作辊磨损预报方法 | |
CN113373295B (zh) | 基于多层感知器的退火炉数字孪生智慧报警***及方法 | |
CN103394520B (zh) | 一种冷轧带钢板形模糊控制方法 | |
CN105512804A (zh) | 冷连轧过程以成本综合控制为目标的乳化液流量设定方法 | |
CN104492822B (zh) | 一种适合于平整轧制的as-ucm机型轧机的板形控制方法 | |
CN102553945A (zh) | 一种适合于四辊轧机的非常态板形预报方法 | |
CN105234188B (zh) | 冷连轧过程中以效益控制为目标的轧制速度优化方法 | |
CN105057364A (zh) | 一种镁合金板材轧制边裂预判及控制方法 | |
CN105013835B (zh) | 冷连轧机组极薄带轧制中基于热凸度的原始辊缝设定方法 | |
CN102671959B (zh) | 六辊平整机组用虚拟板形仪进行板形闭环控制的方法 | |
CN103203371A (zh) | 冷轧机辊缝位置压力双闭环控制方法 | |
CN101722198B (zh) | 一种限动芯棒连轧机组中芯棒限动速度的处理方法 | |
CN104537136A (zh) | 一种六辊轧机辊径大小头缺陷补偿方法 | |
Bu et al. | Rolling schedule multi-objective optimizationbased on influence function for thin gauge steel strip in tandem cold rolling | |
Zhang et al. | Simulation of warping deformation in thin steel strips | |
CN107214199B (zh) | 一种适合于八辊冷轧机组的板形预报方法 | |
CN107526854A (zh) | 一种三辊行星轧机轧制力的计算方法 | |
CN110947771B (zh) | 一种热连轧轧制能耗计算方法 | |
张晓峰 et al. | Model for the generation of surface topography in steel strip temper rolling |
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 | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160511 Address after: 252500, 389, Zhenxing East Road, Liaocheng, Shandong, Guanxian Patentee after: Shandong Guanzhou Co., Ltd. Address before: Hebei Street West Harbor area, 066004 Hebei city of Qinhuangdao province No. 438 Patentee before: Yanshan University |
|
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160302 Termination date: 20180915 |
|
CF01 | Termination of patent right due to non-payment of annual fee |