CN106311754B - 适用于冷连轧机组的乳化液流量动态综合优化设定方法 - Google Patents
适用于冷连轧机组的乳化液流量动态综合优化设定方法 Download PDFInfo
- Publication number
- CN106311754B CN106311754B CN201610825624.6A CN201610825624A CN106311754B CN 106311754 B CN106311754 B CN 106311754B CN 201610825624 A CN201610825624 A CN 201610825624A CN 106311754 B CN106311754 B CN 106311754B
- Authority
- CN
- China
- Prior art keywords
- rack
- emulsion
- flow quantity
- racks
- transferred
- 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
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
-
- 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
- B21B2037/002—Mass flow control
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Control Of Metal Rolling (AREA)
Abstract
一种适用于冷连轧机组的乳化液流量动态综合优化设定方法,它主要包括以下步骤:(1)参数收集;(2)设定各初始值、寻优步长及中间过程参数;(3)定义2#‑5#机架乳化液流量及寻优步长;(4)计算各机架乳化液冷却能力系数;(5)计算稳态轧制速度下各机架理论最佳流量;(6)计算机架乳化液总流量;(7)计算乳化液流量动态综合优化目标函数;(8)输出最优各机架乳化液流量;(9)完成新的轧制速度下各机架乳化液流量的分配,直至完成全部流量的分配。本发明避免了机架间乳化液流量固定分配造成的轧制压力异常增大以及轧制压力“倒挂”等问题,充分发挥了乳化液冷却能力,使得各机架乳化液冷却能力更加均匀,提高了生产的稳定性。
Description
技术领域
本发明属于冷轧技术领域,特别涉及一种冷连轧机组的乳化液流量设定方法。
背景技术
在冷连轧生产过程中,为了降低轧辊与带材的表面温度、减小变形区接触弧表面上摩擦系数和摩擦力、防止金属粘在轧辊表面的同时减少轧辊磨损,需要向轧辊与带材表面喷洒的乳化液。现阶段对于冷连轧机组而言,在总流量与单机架流量不超过最大值的前提下机架之间乳化液流量是根据生产经验按照固定流量比例分配的。考虑到冷连轧机组各个机架带材的强度、压下率、前后张力等轧制工艺以及轧辊直径、表面粗糙度等设备工艺参数是不一样的,因此不同机架保证润滑与冷却所需要的乳化液供给量是不同的,不合理的乳化液流量的分配容易导致各机架冷却能力不均匀,容易造成轧制压力异常增大以及轧制压力“倒挂”等问题。特别的,由于冷连轧过程中轧制速度与摩擦系数、轧制压力、轧制功率、摩擦热、变形热、打滑因子、滑伤指数等关键轧制特征参数之间的关系是非线性的,因此维持冷却与润滑平衡所需要的乳化液的流量与轧制速度之间的关系也是非线性的。而正是由于这种非线性的存在,使得在不同轧制速度下各个机架之间所需要乳化液流量的比例应该是不断变化的。这就是说,冷连轧过程中机架间乳化液流量的分配系数不应该是一个固定值,而应该随着轧制速度的变化动态调整。现阶段,在冷连轧机组乳化液流量分配的过程中,均是单机架乳化液流量固定分配,这样的乳化液流量分配结果,不能很好的保证各机架冷却能力均匀[1-8]。
(参考文献:[1]王康健,郑涛,陈杰.一种冷连轧机组极薄带钢轧制的乳化液流量综合优化方法:中国,CN201410058091.4[P].2015-08-26.[2]白振华,邢雨,李伯阳等.冷连轧过程以成本综合控制为目标的乳化液流量设定方法:中国,CN201510859447.9[P].2016-04-20.[3]王康健,郑涛,陈杰.一种冷连轧机组的工艺润滑制度优化方法:中国,CN201410371135.9[P].2016-02-10.[4]白振华,孙立壮,郭乾勇.冷连轧机组以拉毛防治为目标的工艺润滑制度优化方法:中国,CN201310562473.6[P].2014-03-05.[5]徐俊,邱格君,白振华.高速冷连轧时工艺润滑的优化技术[J].轧钢,2005,22(5):12-14.[6]刘凯,胡滨,李连胜,等.邯钢新冷轧乳化液系统的优化[J].轧钢,2011,28(6):54-56.[7]顾正秋.轧制工艺润滑理论分析[J].辽宁科技大学学报,1981,(3):97-107.[8]李红翠,张虎.济钢冷轧生产中热划伤现象的研究[J].山东冶金,2007,29(6):52-54.)
发明内容
针对冷连轧机组机架之间乳化液流量分配固定,不能随着轧制速度的变化动态调整等问题,本发明提供一种能够充分发挥各机架乳化液的冷却能力、避免乳化液流量供给不足产生的轧制压力异常增大以及轧制压力“倒挂”等问题、提高机组生产稳定性的适用于冷连轧机组的乳化液流量动态综合优化设定方法。
本发明包括以下由计算机执行的步骤:
(a)参数收集;
(a1)收集带材参数,主要包括带材宽度b,各机架入、出口带材厚度Hi、hi,(i=1,2,...,n),各机架带材入、出口速度v0,i、v1,i,(i=1,2,...,n),各机架带材入、出口温度ts0,i、ts1,i,(i=1,2,...,n);
(a2)收集轧制工艺参数,主要包括各机架轧制压力pi,(i=1,2,...,n),各机架压下率εi,(i=1,2,...,n),各机架前张力σ1,i,(i=1,2,...,n),各机架后张力σ0,i,(i=1,2,...,n),轧制速度vi,(i=1,2,...,n);
(a3)收集工艺润滑参数,主要包括乳化液浓度C,乳化液初始温度t0;
(a4)收集轧辊参数,主要包括各机架工作辊直径Di,(i=1,2,...,n);
(b)设定目标函数初始值F0,定义1#机架乳化液流量并初始化w1,min,w1,max,及寻优步长Δw1,设定1#机架乳化液流量寻优中间过程参数k1,并令k1=0;
(c)令w1=w1,min+k1Δw1;
(d)定义2#机架乳化液流量并初始化w2,min,w2,max,及寻优步长Δw2,设定2#机架乳化液流量寻优中间过程参数k2,并令k2=0;
(e)令w2=w2,min+k2Δw2;
(f)定义3#机架乳化液流量并初始化w3,min,w3,max,及寻优步长Δw3,设定3#机架乳化液流量寻优中间过程参数k3,并令k3=0;
(g)令w3=w3,min+k3Δw3;
(h)定义4#机架乳化液流量并初始化w4,min,w4,max,及寻优步长Δw4,设定4#机架乳化液流量寻优中间过程参数k4,并令k4=0;
(i)令w4=w4,min+k4Δw4;
(j)定义5#机架乳化液流量并初始化w5,min,w5,max,及寻优步长Δw5,设定5#机架乳化液流量寻优中间过程参数k5,并令k5=0;
(k)令w5=w5,min+k5Δw5;
(l)计算轧制速度为vj时各机架乳化液冷却能力系数ηi,j,其表达式为
(l1)计算带材轧制过程中的摩擦热Qf,其表达式为Qf=2blμpmnf(SH(l-xn)+Shxn);
式中,b为带钢宽度,mm;l为变形区长度,mm;μ为摩擦系数;pm为轧制区平均单位压力,kN;ηf为摩擦功转化为热的分配系数,一般取2/3;SH为变形区后滑率,Sh为变形区前滑率,v0、v1分别为轧件入出口速度,ms-1;vr为轧辊圆周速度,ms-1;α为变形区咬入角,rad;xn为中性面到出口的距离,mm。
(l2)计算带材轧制过程中的变形热Qb,其表达式为Qb=pmηmVln(H/h);
式中,ηm为塑性变形功转化为热的分配系数,一般取0.9;V为变形区带材的体积,mm3;H、h分别为带材入、出口厚度,mm。
(l3)计算带材轧制过程中上游机架残留的热量Qc,其表达式为
式中:Qc,i为第i机架所在的上游机架残留的热量;Qf,i-1为第i-1机架带材轧制过程中的摩擦热;Qb,i-1为第i-1机架带材轧制过程中的变形热;Qr,i-1为第i-1机架带材轧制过程中乳化液带走的热量;Qd,i-1为第i-1机架与第i机架之间带钢表面热量的散失量。
(l4)计算两个机架之间带钢表面热量的散失量Qd,其表达式为Qd=ρcmV(ts1-ts0);
式中:ρ为带材密度,对于钢材ρ=7800Kg/m3;cm为带材的比热,对于钢材可以取cm=470J/(Kg.℃);ts1为前部机架带材的出口温度,℃;ts0为后部机架带材的入口温度,℃。
(l5)计算带材轧制过程中乳化液带走的热量Qr,其表达式为Qr=αBA(tm-t0);
式中:αB为换热系数,αB=αB0w0.264T-0.213exp(9.45-0.1918C)×1.163;αB0为喷嘴形状、喷射角度影响系数;C为乳化液浓度,%;w为乳化液流量,L/min;T为带材温度,℃;A为变形区接触面积,mm2;tm为带钢与轧辊的平均温度,℃;t0为乳化液初始温度,℃。
(m)计算冷连轧机组第i机架稳态轧制速度下保证Qri=Qfi+Qbi+Qci时的理论最佳流量wli,(i=1,2,...,n);
(n)判断不等式wi≤wli是否成立?如果成立,转入步骤(o),否则直接转入步骤(s);
(o)计算轧制速度为vj时机架乳化液总流量W,其表达式为
(p)判断不等式W≤Wmax是否成立?如果成立,转入步骤(q),否则直接转入步骤(s);
(q)计算轧制速度为vj时乳化液流量动态综合优化设定的目标函数F(X),其表达式为
式中,α为加权系数;为冷连轧机组轧制速度为vj时整体平均乳化液冷却能力系数,
(r)判断不等式F<F0是否成立?如果成立,则令F0=F,1#机架乳化液流量w1y=w1,2#机架乳化液流量w2y=w2,3#机架乳化液流量w3y=w3,4#机架乳化液流量w4y=w4,5#机架乳化液流量w5y=w5,转入步骤(s),否则,直接转入步骤(s);
(s)判断不等式是否成立?如果不等式成立,则令k5=k5+1,转入步骤(k),否则转入步骤(t);
(t)判断不等式是否成立?如果不等式成立,则令k4=k4+1,转入步骤(i),否则转入步骤(u);
(u)判断不等式是否成立?如果不等式成立,则令k3=k3+1,转入步骤(g),否则转入步骤(v);
(v)判断不等式是否成立?如果不等式成立,则令k2=k2+1,转入步骤(e),否则转入步骤(w);
(w)判断不等式是否成立?如果不等式成立,则令k1=k1+1,转入步骤(c),否则转入步骤(x);
(x)输出轧制速度为vj时最优1#机架乳化液流量w1y,j,最优2#机架乳化液流量w2y,j,最优3#机架乳化液流量w3y,j,最优4#机架乳化液流量w4y,j,最优5#机架乳化液流量w5y,j。
(y)令vj=vj+1,转入步骤(b),完成新的轧制速度vj+1情况下各机架乳化液流量的分配,直至最终完成全部轧制速度情况下各机架乳化液流量的分配。
本发明与现有技术相比具有如下优点:
实现了冷连轧机组乳化液流量随着轧制速度进行动态综合优化设定,避免了机架间乳化液流量固定分配造成的轧制压力异常增大以及轧制压力“倒挂”等问题,充分发挥了乳化液冷却能力,同时使得各机架乳化液冷却能力更加均匀,提高了生产的稳定性,为生产企业带来了效益。
附图说明
图1是本发明总计算流程图;
图2是本发明实施例1优化后各机架乳化液流量随着轧制速度的变化曲线图;
图3是本发明实施例2优化后各机架乳化液流量随着轧制速度的变化曲线图。
图中:v-out—轧制速度曲线、w-1#—1#机架流量曲线、w-2#—2#机架流量曲线、w-3#—3#机架流量曲线、w-4#—4#机架流量曲线、w-5#—5#机架流量曲线。
具体实施方式
实施例1
一种适用于冷连轧机组的乳化液流量动态综合优化设定方法,其计算流程如图1所示:
首先,在步骤(a)中,参数收集:
(a1)收集带材参数,主要包括带材宽度b=850mm,各机架入、出口带材厚度Hi={2.1,1.17,0.65,0.40,0.27}(mm)、hi={1.17,0.65,0.40,0.27,0.22}(mm),各机架带材典型入、出口速度(m/s)
各机架带材典型入、出口温度(℃)
(a2)收集轧制工艺参数,主要包括各机架典型轧制压力(kN)
各机架压下率εi={1#44.3%;2#44.2%;3#38.1%;4#32.5%;5#19.3%},
各机架典型前张力(MPa)
各机架典型后张力(MPa)
典型轧制速度vi={260,751.4,1038.9,1371,1156}(m/s);
(a3)收集工艺润滑参数,主要包括乳化液浓度C=5.6%,乳化液初始温度t0=20℃;
(a4)收集轧辊参数,主要包括各机架工作辊直径
Di={1#435;2#435;3#435;4#435;5#435}(mm);
随后,在步骤(b)中,设定目标函数初始值F0=1.0×1020,定义1#机架乳化液流量并初始化w1,min=3000L/min,w1,max=4000L/min,及寻优步长Δw1=50L/min,设定1#机架乳化液流量寻优中间过程参数k1,并令k1=0;
随后,在步骤(c)中,令w1=w1,min+k1Δw1=3000L/min;
随后,在步骤(d)中,定义2#机架乳化液流量并初始化w2,min=3000L/min,w2,max=4000L/min,及寻优步长Δw2=50L/min,设定2#机架乳化液流量寻优中间过程参数k2,并令k2=0;
随后,在步骤(e)中,令w2=w2,min+k2Δw2=3000L/min;
随后,在步骤(f)中,定义3#机架乳化液流量并初始化w3,min=3000L/min,w3,max=4000L/min,及寻优步长Δw3=50L/min,设定3#机架乳化液流量寻优中间过程参数k3,并令k3=0;
随后,在步骤(g)中,令w3=w3,min+k3Δw3=3000L/min;
随后,在步骤(h)中,定义4#机架乳化液流量并初始化w4,min=3500L/min,w4,max=5000L/min,及寻优步长Δw4=50L/min,设定4#机架乳化液流量寻优中间过程参数k4,并令k4=0;
随后,在步骤(i)中,令w4=w4,min+k4Δw4=3500L/min;
随后,在步骤(j)中,定义5#机架乳化液流量并初始化w5,min=3500L/min,w5,max=5000L/min,及寻优步长Δw5=50L/min,设定5#机架乳化液流量寻优中间过程参数k5,并令k5=0;
随后,在步骤(k)中,令w5=w5,min+k5Δw5=3500L/min;
随后,在步骤(l)中,计算轧制速度vj=260m/s时各机架乳化液冷却能力系数ηi,j={0.32,0.38,0.45,0.42,0.47};
其中,在步骤(l1)中,计算轧制速度vj=260m/s时各机架带材在轧制过程中的摩擦热Qfi,j={59726,52707,47353,43934,41096}(KJ);
其中,在步骤(l2)中,计算轧制速度vj=260m/s时各机架带材在轧制过程中的变形热Qbi,j={2100,1036,444,213,76}(KJ);
其中,在步骤(l3)中,计算轧制速度vj=260m/s时各机架带材在轧制过程中的上游机架残留的热量Qci,j={0,15627,2095,17356,263}(KJ);
其中,在步骤(l4)中,计算轧制速度vj=260m/s时两个机架之间带钢表面热量的散失量Qdi,j={1#-2#:4157,2#-3#:2700,3#-4#:2058,4#-5#:923}(KJ);
其中,在步骤(l5)中,计算轧制速度vj=260m/s时带材轧制过程中乳化液带走的热量Qri,j={42042,48948,28383,42961,22084}(KJ);
随后,在步骤(m)中,计算冷连轧机组第i机架稳态轧制速度下保证Qri=Qfi+Qbi+Qci时的理论最佳流量wli={3615L/min,3897L/min,4212L/min,4751L/min,5208L/min};
随后,在步骤(n)中,判断不等式wi≤wli是否成立?由步骤(m)计算结果可知,不等式结果显然成立,转入步骤(o);
随后,在步骤(o)中,计算轧制速度vj=260m/s时机架乳化液总流量W=16000L/min;
随后,在步骤(p)中,判断不等式W≤Wmax=34000L/min是否成立?由步骤(o)计算结果可知,不等式显然成立,转入步骤(q);
随后,在步骤(q)中,计算轧制速度vj=260m/s时乳化液流量动态综合优化设定的目标函数F(X)=0.218;
随后,在步骤(r)中,判断不等式F<F0=1.0×1020是否成立?由步骤(q)计算结果可知,不等式显然成立,则令F0=0.218,1#机架乳化液流量w1y=3000L/min,2#机架乳化液流量w2y=3000L/min,3#机架乳化液流量w3y=3000L/min,4#机架乳化液流量w4y=3500L/min,5#机架乳化液流量w5y=3500L/min,转入步骤(s);
随后,在步骤(s)中,判断不等式是否成立?由于此时k5=1,不等式显然成立,则令k5=2,转入步骤(k)继续循环计算,直到k5=30时转入步骤(t);
随后,在步骤(t)中,判断不等式是否成立?由于此时k4=1,不等式显然成立,则令k4=2,转入步骤(i)继续循环计算,直到k4=30时转入步骤(u);
随后,在步骤(u)中,判断不等式是否成立?由于此时k3=1,不等式显然成立,则令k3=2,转入步骤(g)继续循环计算,直到k3=20时转入步骤(v);
随后,在步骤(v)中,判断不等式是否成立?由于此时k2=1,不等式显然成立,则令k2=2,转入步骤(e)继续循环计算,直到k2=20时转入步骤(w);
随后,在步骤(w)中,判断不等式是否成立?由于此时k1=1,不等式显然成立,则令k1=2,转入步骤(c)继续循环计算,直到k1=20时转入步骤(x);
随后,在步骤(x)中,输出轧制速度vj=260m/s时最优1#机架乳化液流量w1y=3100L/min,最优2#机架乳化液流量w2y=3200L/min,最优3#机架乳化液流量w3y=3350L/min,最优4#机架乳化液流量w4y=3600L/min,最优5#机架乳化液流量w5y=4400L/min。
随后,在步骤(y)中,令vj=vj+1,转入步骤(b),完成新的轧制速度vj+1情况下各机架乳化液流量的分配,直至最终完成全部轧制速度情况下各机架乳化液流量的分配。优化后各机架乳化液流量随着轧制速度的变化曲线如图2所示。
实施例2
首先,在步骤(a)中,参数收集:
(a1)收集带材参数,主要包括带材宽度b=975mm,各机架入、出口带材厚度Hi={2.1,1.17,0.65,0.40,0.27}(mm)、hi={1.17,0.65,0.40,0.27,0.22}(mm),各机架带材典型入、出口速度(m/s)
各机架带材典型入、出口温度(℃)
(a2)收集轧制工艺参数,主要包括各机架典型轧制压力(kN)
各机架压下率εi={1#44.3%;2#44.2%;3#38.1%;4#32.5%;5#19.3%},
各机架典型前张力(MPa)
各机架典型后张力(MPa)
典型轧制速度vi={210,700.4,987.9,1320,1102}(m/s);
(a3)收集工艺润滑参数,主要包括乳化液浓度C=5.6%,乳化液初始温度t0=20℃;
(a4)收集轧辊参数,主要包括各机架工作辊直径
Di={1#435;2#435;3#435;4#435;5#435}(mm);
随后,在步骤(b)中,设定目标函数初始值F0=1.0×1020,定义1#机架乳化液流量并初始化w1,min=3000L/min,w1,max=4000L/min,及寻优步长Δw1=50L/min,设定1#机架乳化液流量寻优中间过程参数k1,并令k1=0;
随后,在步骤(c)中,令w1=w1,min+k1Δw1=3000L/min;
随后,在步骤(d)中,定义2#机架乳化液流量并初始化w2,min=3000L/min,w2,max=4000L/min,及寻优步长Δw2=50L/min,设定2#机架乳化液流量寻优中间过程参数k2,并令k2=0;
随后,在步骤(e)中,令w2=w2,min+k2Δw2=3000L/min;
随后,在步骤(f)中,定义3#机架乳化液流量并初始化w3,min=3000L/min,w3,max=4000L/min,及寻优步长Δw3=50L/min,设定3#机架乳化液流量寻优中间过程参数k3,并令k3=0;
随后,在步骤(g)中,令w3=w3,min+k3Δw3=3000L/min;
随后,在步骤(h)中,定义4#机架乳化液流量并初始化w4,min=3500L/min,w4,max=5000L/min,及寻优步长Δw4=50L/min,设定4#机架乳化液流量寻优中间过程参数k4,并令k4=0;
随后,在步骤(i)中,令w4=w4,min+k4Δw4=3500L/min;
随后,在步骤(j)中,定义5#机架乳化液流量并初始化w5,min=3500L/min,w5,max=5000L/min,及寻优步长Δw5=50L/min,设定5#机架乳化液流量寻优中间过程参数k5,并令k5=0;
随后,在步骤(k)中,令w5=w5,min+k5Δw5=3500L/min;
随后,在步骤(l)中,计算轧制速度vj=210m/s时各机架乳化液冷却能力系数ηi,j={0.41,0.48,0.51,0.49,0.56};
其中,在步骤(l1)中,计算轧制速度vj=210m/s时各机架带材在轧制过程中的摩擦热Qfi,j={78581,69348,62306,57807,54065}(KJ);
其中,在步骤(l2)中,计算轧制速度vj=210m/s时各机架带材在轧制过程中的变形热Qbi,j={2763,1364,584,280,100}(KJ);
其中,在步骤(l3)中,计算轧制速度vj=210m/s时各机架带材在轧制过程中的上游机架残留的热量Qci,j={0,26794,3830,26202,1120}(KJ);
其中,在步骤(l4)中,计算轧制速度vj=210m/s时两个机架之间带钢表面热量的散失量Qdi,j={1#-2#:6557,2#-3#:3318,3#-4#:2042,4#-5#:1140}(KJ);
其中,在步骤(l5)中,计算轧制速度vj=210m/s时带材轧制过程中乳化液带走的热量Qri,j={47993,63564,34646,55827,24953}(KJ);
随后,在步骤(m)中,计算冷连轧机组第i机架稳态轧制速度下保证Qri=Qfi+Qbi+Qci时的理论最佳流量wli={3782L/min,3991L/min,4365L/min,4923L/min,5518L/min};
随后,在步骤(n)中,判断不等式wi≤wli是否成立?由步骤(m)计算结果可知,不等式结果显然成立,转入步骤(o);
随后,在步骤(o)中,计算轧制速度vj=210m/s时机架乳化液总流量W=16000L/min;
随后,在步骤(p)中,判断不等式W≤Wmax=34000L/min是否成立?由步骤(o)计算结果可知,不等式显然成立,转入步骤(q);
随后,在步骤(q)中,计算轧制速度vj=210m/s时乳化液流量动态综合优化设定的目标函数F(X)=0.365;
随后,在步骤(r)中,判断不等式F<F0=1.0×1020是否成立?由步骤(q)计算结果可知,不等式显然成立,则令F0=0.365,1#机架乳化液流量w1y=3000L/min,2#机架乳化液流量w2y=3000L/min,3#机架乳化液流量w3y=3000L/min,4#机架乳化液流量w4y=3500L/min,5#机架乳化液流量w5y=3500L/min,转入步骤(s);
随后,在步骤(s)中,判断不等式是否成立?由于此时k5=1,不等式显然成立,则令k5=2,转入步骤(k)继续循环计算,直到k5=30时转入步骤(t);
随后,在步骤(t)中,判断不等式是否成立?由于此时k4=1,不等式显然成立,则令k4=2,转入步骤(i)继续循环计算,直到k4=30时转入步骤(u);
随后,在步骤(u)中,判断不等式是否成立?由于此时k3=1,不等式显然成立,则令k3=2,转入步骤(g)继续循环计算,直到k3=20时转入步骤(v);
随后,在步骤(v)中,判断不等式是否成立?由于此时k2=1,不等式显然成立,则令k2=2,转入步骤(e)继续循环计算,直到k2=20时转入步骤(w);
随后,在步骤(w)中,判断不等式是否成立?由于此时k1=1,不等式显然成立,则令k1=2,转入步骤(c)继续循环计算,直到k1=20时转入步骤(x);
随后,在步骤(x)中,输出轧制速度vj=210m/s时最优1#机架乳化液流量w1y=3150L/min,最优2#机架乳化液流量w2y=3200L/min,最优3#机架乳化液流量w3y=3400L/min,最优4#机架乳化液流量w4y=3650L/min,最优5#机架乳化液流量w5y=4400L/min。
随后,在步骤(y)中,令vj=vj+1,转入步骤(b),完成新的轧制速度vj+1情况下各机架乳化液流量的分配,直至最终完成全部轧制速度情况下各机架乳化液流量的分配。优化后各机架乳化液流量随着轧制速度的变化曲线如图3所示。
Claims (1)
1.一种适用于冷连轧机组的乳化液流量动态综合优化设定方法,其特征在于:它包括以下由计算机执行的步骤:
(a)参数收集;
(a1)收集带材参数,主要包括带材宽度b,各机架入、出口带材厚度Hi、hi,(i=1,2,...,n),各机架带材入、出口速度v0,i、v1,i,(i=1,2,...,n),各机架带材入、出口温度ts0,i、ts1,i,(i=1,2,...,n);
(a2)收集轧制工艺参数,主要包括各机架轧制压力pi,(i=1,2,...,n),各机架压下率εi,(i=1,2,...,n),各机架前张力σ1,i,(i=1,2,...,n),各机架后张力σ0,i,(i=1,2,...,n),轧制速度vi,(i=1,2,...,n);
(a3)收集工艺润滑参数,主要包括乳化液浓度C,乳化液初始温度t0;
(a4)收集轧辊参数,主要包括各机架工作辊直径Di,(i=1,2,...,n);
(b)设定目标函数初始值F0,定义1#机架乳化液流量最小值w1,min,1#机架乳化液流量最大值w1,max,及寻优步长Δw1,设定1#机架乳化液流量寻优中间过程参数k1,并令k1=0;
(c)令w1=w1,min+k1Δw1;
(d)定义2#机架乳化液流量最小值w2,min,2#机架乳化液流量最大值w2,max,及寻优步长Δw2,设定2#机架乳化液流量寻优中间过程参数k2,并令k2=0;
(e)令w2=w2,min+k2Δw2;
(f)定义3#机架乳化液流量最小值w3,min,3#机架乳化液流量最大值w3,max,及寻优步长Δw3,设定3#机架乳化液流量寻优中间过程参数k3,并令k3=0;
(g)令w3=w3,min+k3Δw3;
(h)定义4#机架乳化液流量最小值w4,min,4#机架乳化液流量最大值w4,max,及寻优步长Δw4,设定4#机架乳化液流量寻优中间过程参数k4,并令k4=0;
(i)令w4=w4,min+k4Δw4;
(j)定义5#机架乳化液流量最小值w5,min,5#机架乳化液流量最大值w5,max,及寻优步长Δw5,设定5#机架乳化液流量寻优中间过程参数k5,并令k5=0;
(k)令w5=w5,min+k5Δw5;
(l)计算轧制速度为vj时各机架乳化液冷却能力系数ηi,j,其表达式为
(11)计算带材轧制过程中的摩擦热Qf,其表达式为Qf=2blμpmηf(SH(l-xn)+Shxn);
式中,b为带钢宽度,mm;l为变形区长度,mm;μ为摩擦系数;pm为轧制区平均单位压力,kN;ηf为摩擦功转化为热的分配系数,取2/3;SH为变形区后滑率,%;Sh为变形区前滑率,%;v0、v1分别为轧件入出口速度,ms-1;vr为轧辊圆周速度,ms-1;α为变形区咬入角,rad;xn为中性面到出口的距离,mm;
(l2)计算带材轧制过程中的变形热Qb,其表达式为Qb=pmηmVln(H/h);
式中,ηm为塑性变形功转化为热的分配系数,取0.9;V为变形区带材的体积,mm3;H、h分别为带材入、出口厚度,mm;
(l3)计算带材轧制过程中上游机架残留的热量Qc,其表达式为
式中:Qc,i为第i机架所在的上游机架残留的热量;Qf,i-1为第i-1机架带材轧制过程中的摩擦热;Qb,i-1为第i-1机架带材轧制过程中的变形热;Qr,i-1为第i-1机架带材轧制过程中乳化液带走的热量;Qd,i-1为第i-1机架与第i机架之间带钢表面热量的散失量;
(l4)计算两个机架之间带钢表面热量的散失量Qd,其表达式为Qd=ρcmV(ts1-ts0);
式中:ρ为带材密度,对于钢材ρ=7800Kg/m3;cm为带材的比热,对于钢材取cm=470J/(Kg.℃);ts1为前部机架带材的出口温度,℃;ts0为后部机架带材的入口温度,℃;
(l5)计算带材轧制过程中乳化液带走的热量Qr,其表达式为Qr=αBA(tm-t0);
式中:αB为换热系数,αB=αB0w0.264T-0.213exp(9.45-0.1918C)×1.163;αB0为喷嘴形状、喷射角度影响系数;C为乳化液浓度,%;w为乳化液流量,L/min;T为带材温度,℃;A为变形区接触面积,mm2;tm为带钢与轧辊的平均温度,℃;t0为乳化液初始温度,℃;
(m)计算冷连轧机组第i机架稳态轧制速度下保证Qri=Qfi+Qbi+Qci时的理论最佳流量wli,(i=1,2,...,n);
(n)判断不等式wi≤wli是否成立?如果成立,转入步骤(o),否则直接转入步骤(s);
(o)计算轧制速度为vj时机架乳化液总流量W,其表达式为
(p)以Wmax表示冷连轧机组机架乳化液总流量最大值,判断不等式W≤Wmax是否成立?如果成立,转入步骤(q),否则直接转入步骤(s);
(q)计算轧制速度为vj时乳化液流量动态综合优化设定的目标函数F(X),其表达式为
式中,α为加权系数;为冷连轧机组轧制速度为vj时整体平均乳化液冷却能力系数,
(r)判断不等式F<F0是否成立?如果成立,则令F0=F,1#机架乳化液流量w1y=w1,2#机架乳化液流量w2y=w2,3#机架乳化液流量w3y=w3,4#机架乳化液流量w4y=w4,5#机架乳化液流量w5y=w5,转入步骤(s),否则,直接转入步骤(s);
(s)判断不等式是否成立?如果不等式成立,则令k5=k5+1,转入步骤(k),否则转入步骤(t);
(t)判断不等式是否成立?如果不等式成立,则令k4=k4+1,转入步骤(i),否则转入步骤(u);
(u)判断不等式是否成立?如果不等式成立,则令k3=k3+1,转入步骤(g),否则转入步骤(v);
(v)判断不等式是否成立?如果不等式成立,则令k2=k2+1,转入步骤(e),否则转入步骤(w);
(w)判断不等式是否成立?如果不等式成立,则令k1=k1+1,转入步骤(c),否则转入步骤(x);
(x)输出轧制速度为vj时最优1#机架乳化液流量w1y,j,最优2#机架乳化液流量w2y,j,最优3#机架乳化液流量w3y,j,最优4#机架乳化液流量w4y,j,最优5#机架乳化液流量w5y,j;
(y)令vj=vj+1,转入步骤(b),完成新的轧制速度vj+1情况下各机架乳化液流量的分配,直至最终完成全部轧制速度情况下各机架乳化液流量的分配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610825624.6A CN106311754B (zh) | 2016-09-14 | 2016-09-14 | 适用于冷连轧机组的乳化液流量动态综合优化设定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610825624.6A CN106311754B (zh) | 2016-09-14 | 2016-09-14 | 适用于冷连轧机组的乳化液流量动态综合优化设定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106311754A CN106311754A (zh) | 2017-01-11 |
CN106311754B true CN106311754B (zh) | 2018-07-17 |
Family
ID=57787920
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610825624.6A Active CN106311754B (zh) | 2016-09-14 | 2016-09-14 | 适用于冷连轧机组的乳化液流量动态综合优化设定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106311754B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106909723B (zh) * | 2017-02-16 | 2020-03-31 | 燕山大学 | 冷轧过程乳化液流量与轧制速度关系曲线优化设定方法 |
CN108160723B (zh) * | 2017-12-07 | 2019-06-28 | 燕山大学 | 一种冷轧乳化液直喷系统喷嘴喷射高度优化方法 |
CN107971355B (zh) * | 2017-12-07 | 2019-03-22 | 燕山大学 | 一种二次冷轧机组乳化液喷嘴的喷射角度优化设定方法 |
CN108160722B (zh) * | 2017-12-07 | 2019-07-26 | 燕山大学 | 一种二次冷轧机组乳化液喷嘴喷射方向角的优化方法 |
CN108160721B (zh) * | 2017-12-07 | 2019-03-22 | 燕山大学 | 一种乳化液喷淋架上喷嘴间距与喷嘴数量综合优化方法 |
CN108380676B (zh) * | 2018-01-25 | 2019-04-09 | 燕山大学 | 二次冷轧过程带钢表面乳化液流量横向分布预报方法 |
CN108723097B (zh) * | 2018-04-10 | 2019-10-11 | 燕山大学 | Dcr机组大变形下以稳轧为目标的轧制参数优化方法 |
CN108714627B (zh) * | 2018-04-10 | 2019-10-11 | 燕山大学 | Dcr机组大变形下以稳轧为目标的工艺润滑优化方法 |
CN108838211B (zh) * | 2018-06-29 | 2020-09-18 | 首钢京唐钢铁联合有限责任公司 | 一种冷连轧过程的乳化液浓度优化方法 |
CN110842031B (zh) * | 2018-07-24 | 2020-10-27 | 宝山钢铁股份有限公司 | 一种抑制冷连轧机组振动的乳化液流量优化方法 |
CN110961464B (zh) * | 2018-09-29 | 2022-06-28 | 宝山钢铁股份有限公司 | 冷连轧机组以振动抑制为目标的乳化液浓度优化方法 |
CN110116059B (zh) * | 2019-05-13 | 2020-07-14 | 燕山大学 | 一种适用于dos油静电涂油机旋转喷雾器高度优化的方法 |
CN110624957B (zh) * | 2019-09-30 | 2021-03-23 | 江苏九天光电科技有限公司 | 湿平整机组以粗糙度为目标的工艺润滑制度控制方法 |
CN111672915B (zh) * | 2020-05-21 | 2022-03-15 | 武汉定飞科技有限公司 | 一种可逆冷轧机供液泵节能运行的方法 |
CN111872131B (zh) * | 2020-07-27 | 2022-04-29 | 广西钢铁集团有限公司 | 动态调整冷连轧机乳化液流量的方法 |
CN113319137B (zh) * | 2021-06-03 | 2022-04-05 | 宝钢湛江钢铁有限公司 | 六机架冷连轧机组超高强钢工艺润滑制度综合优化方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU899181A1 (ru) * | 1980-05-20 | 1982-01-23 | Киевский институт автоматики им.ХХУ съезда КПСС | Система регулировани профил валков |
JPH0716630A (ja) * | 1993-06-22 | 1995-01-20 | Kawasaki Steel Corp | 圧延制御方法 |
DE20006508U1 (de) * | 2000-04-08 | 2000-08-31 | Achenbach Buschhütten GmbH, 57223 Kreuztal | Walzenkühl- und/oder Schmiervorrichtung für Kaltbandwalzwerke, insbesondere Feinband- und Folienwalzwerke |
CN101683660B (zh) * | 2008-09-28 | 2011-07-20 | 宝山钢铁股份有限公司 | 冷连轧机乳化液分段冷却控制方法 |
CN103586286B (zh) * | 2013-11-12 | 2015-06-10 | 燕山大学 | 冷连轧机组以拉毛防治为目标的轧制规程综合优化方法 |
CN105512804B (zh) * | 2015-12-01 | 2019-03-22 | 燕山大学 | 冷连轧过程以成本综合控制为目标的乳化液流量设定方法 |
-
2016
- 2016-09-14 CN CN201610825624.6A patent/CN106311754B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106311754A (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106311754B (zh) | 适用于冷连轧机组的乳化液流量动态综合优化设定方法 | |
CN104858241B (zh) | 一种冷连轧机组极薄带钢轧制的乳化液流量综合优化方法 | |
RU2352414C1 (ru) | Способ подачи смазочного масла при холодной прокатке | |
CN106909723A (zh) | 冷轧过程乳化液流量与轧制速度关系曲线优化设定方法 | |
WO2020020191A1 (zh) | 一种抑制冷连轧机组振动的乳化液流量优化方法 | |
CN106269888B (zh) | 一种实现esp精轧机组在线换辊的逆流换辊方法 | |
CN101214501B (zh) | 兼顾带钢凸度、边降控制和磨损控制的工作辊及使用方法 | |
CN104785538B (zh) | 一种冷连轧机组极薄带钢轧制的压下规程优化方法 | |
CN104289530B (zh) | 双机架四辊轧机的机架间乳化液流量分配方法 | |
EP1870172A4 (en) | METHOD FOR CONTINUOUS CASTING AND ROLLING OF A MEDIUM SHEET | |
CN108405625B (zh) | 一种实现esp精轧机组在线换辊的顺流换辊方法 | |
CN100518976C (zh) | 一种中间坯加热器优化控制方法 | |
CN103586286A (zh) | 冷连轧机组以拉毛防治为目标的轧制规程综合优化方法 | |
CN101602067A (zh) | 五机架ucm冷连轧机组板形与板凸度在线综合控制方法 | |
CN101934290A (zh) | 不锈钢冷连轧负荷分配调整方法 | |
CN108856305A (zh) | 一种无取向硅钢生产轧机负荷分配方法 | |
CN104289527B (zh) | 双四辊机组汽车板冷轧中乳化液浓度优化设定方法 | |
CN103418619A (zh) | 一种冷轧带钢板形预测控制方法 | |
CN107321797A (zh) | 短流程esp精轧机组在线换辊方法 | |
CN110434172A (zh) | 一种炉卷和精轧机组连轧的负荷分配计算方法 | |
CN113333477B (zh) | 一种esp精轧机组在线换辊与动态变规程时控制辊缝的方法 | |
CN104289525A (zh) | 双机架六辊轧机冷轧中乳化液总流量设定方法 | |
CN108838211A (zh) | 一种冷连轧过程的乳化液浓度优化方法 | |
CN107321800B (zh) | 基于动态变规格的带钢厚度在线动态调整控制方法 | |
CN103433298A (zh) | 薄规格轧制集装箱板尾部控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20181024 Address after: 063600 Lingang Industrial Zone, Leting County, Tangshan City, Hebei Patentee after: Tangshan City Rong Ze iron and steel processing Co., Ltd. Address before: 066004 438 west section of Hebei Avenue, Qinhuangdao, Hebei. Patentee before: Yanshan University |
|
TR01 | Transfer of patent right |