CN103302108B - 一种六辊平整机组轧制力优化设定方法 - Google Patents
一种六辊平整机组轧制力优化设定方法 Download PDFInfo
- Publication number
- CN103302108B CN103302108B CN201310276037.2A CN201310276037A CN103302108B CN 103302108 B CN103302108 B CN 103302108B CN 201310276037 A CN201310276037 A CN 201310276037A CN 103302108 B CN103302108 B CN 103302108B
- Authority
- CN
- China
- Prior art keywords
- planisher
- roll
- unit
- force
- enter step
- 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
Landscapes
- Metal Rolling (AREA)
- Control Of Metal Rolling (AREA)
Abstract
一种六辊平整机组轧制力优化设定方法,属于带钢平整生产工艺和控制技术。经过大量的生产试验与理论研究,充分考虑到六辊平整机平整轧制过程的设备与工艺特点,根据该方法,为每个钢种每个规格确定一组平整机最优轧制力设定值,保证在第一次生产新规格产品时,***可以给定合理的轧制力设定值。本计算在保证打滑因子在临界点以内、计算平整机轧制力冲击响应幅值不超负荷,最终结合实际生产情况,确定出平整机轧制力设定值,减少因为设定不准造成的生产不稳定、轧断等不良现象出现。可以提高平整机板形控制稳定性,提高成材率。
Description
技术领域
本发明属于带钢平整生产工艺和控制技术,特别涉及一种六辊平整机组轧制力优化设定方法。
背景技术
随着社会的工业化程度越来越广泛,社会经济发展水平大幅度提高,带动了冶金领域板带钢(家电板和汽车板等)的市场需求,同时也对板带钢质量提出了更高的要求。平整作为冷轧板带生产中一道后序精加工艺,对产品质量的控制尤为重要。平整工序可以通过控制轧制力、张力和弯辊力等工艺参数,以保证退火后产品的机械性能与外形质量。所以平整机的运行状态最大程度地决定了产品成材率,影响企业的经济效益。
通常工程应用中,平整机控制根据设定的延伸率对轧制力进行闭环控制,控制过程中对于平整机组轧制力的设定仅仅考虑到板形治理与最大轧制压力的限制幅度,给出一个所有规格带钢生产时的平均值。在每次平整机起停车时,平整机控制***先把设定的轧制力作为初始值进行控制。因此,该轧制力初始设定值的准确程度直接影响到控制期间产品延伸率、延伸率调整到目标值的过渡时间,以及过渡期间的板形好坏。
由于当前一些平整机轧制力设定值误差较大,导致闭环***在控制初期会产生较大幅度的调整量,轧制力、延伸率的波动剧烈,板形也受影响。受反馈周期、设备响应周期以及***灵敏性等因素影响,控制***需要较长一段时间才能将轧制力调节到合适值,有时需要近30s,对带钢的影响长度约30m左右。本发明专利针对该情况,通过理论研究和工程数据分析相结合的方法,找到了一种适合对平整机轧制力设定值进行优化的方法,并通过生产应用,显著改善了控制***初期的平整机轧制力跳跃情况。
发明内容
本发明的目的在于提供一种六辊平整机组轧制力优化设定方法,充分考虑到六辊平整机组平整轧制过程的设备与工艺特点,在保证平整机组不出现打滑与冲击响应满足要求的前提下,给出合适的平整机轧制力设定值。提高成品带钢的板形与延伸率达标率,提高产品质量和成材率。
本发明的技术方案包括以下步骤(如附图1所示):
1、针对所分析平整机组,统计整理其设备及主要工艺参数,包括工作辊直径DW,支撑辊、中间辊和工作辊的等效质量Mb、Mm、Mw,平整机牌坊的等效刚度Km,支承辊与中间辊间接触刚度Kb,工作辊与中间辊间接触刚度Kw,轧制压力最大值Pmax,轧制压力最小值Pmin,进入步骤2;
2、根据收集的实际生产数据,以及建立的轧制力计算方程,采用回归拟合的方法,计算出带钢的平整轧制力计算公式;进入步骤3;
式中:P——总轧制力,单位: N;f——单位宽度轧制力,单位:N/m;B——带钢宽度,单位:m;L——轧制变形区中轧辊与带钢接触弧长度,单位:mm;a0——平整钢种影响系数;a1——工况影响系数;D——工作辊直径,单位:mm;ε——带钢延伸率;μ——摩擦系数;h0——入口厚度,单位:mm;σp——当量变形抗力,单位:Pa;e——应变速率;v——轧制速度,单位:m/s;k3——变形抗力影响系数;k2——后张力影响系数;k1——前张力影响系数;σs——带钢屈服强度;单位:Pa;σ0——前张力,单位:Pa;σ1——后张力,单位:Pa;k4、k5、k6——平整机速度修正系数。
3、确定平整机的带钢厚度h0,带钢宽度B,带钢平整延伸率ε,带钢屈服强度σs,产线速度v,平整机前后张应力σ0、σ1,根据步骤2中的公式计算出平整机轧制力设定值P,如果计算结果大于Pmax,则轧制力P=Pmax;当计算结果小于Pmin,则轧制力P=Pmin;进入步骤4;
4、根据公式,计算出平整机咬入角,根据公式:,计算出中性角,进入步骤5;
式中:f——单位压力;B——带钢宽度(m);Q1,Q0——作用在轧件上的前张力、后张力(N); α——咬入角(弧度);
5、计算出打滑因子fγ=γ/α,判断是否满足0.2<fγ<0.8?如果满足,则进入步骤6;如果不满足,则令单位宽度轧制力f=f-sign(γ/α-1/2)*sign(Q1-Q0)*10000N,进入步骤4,如果循环次数大于1000仍不满足0.2<fγ<0.8,取最后一次迭代循环的单位宽度轧制力 f 值。进入步骤6;
6、根据平整机结构,建立3自由度振动模型,力学模型详见附图2,采用MATLAB软件求出振动方程的数值解,计算出施加平整机轧制力后的平整机的纵向振动响应,得到振动时轧制力最大幅值Pa,判断Pa< Pmax是否满足?如果满足,则进入步骤7;如果不满足,则P=P-100000N,如果此时满足P> Pmin,则进入步骤4;如果P> Pmin条件不满足,则P=Pmin,进入步骤7;
7、将完成计算得到的轧制力设定值P写入平整机轧制力设定表,并植入控制***。
本发明经过大量的生产试验与理论研究,充分考虑到六辊平整机平整轧制过程的设备与工艺特点,发明了一种平整机轧制力优化设定方法,根据该方法,为每个钢种每个规格确定一组平整机最优轧制力设定值,保证在第一次生产新规格产品时,***可以给定合理的轧制力设定值。本计算在保证打滑因子在临界点以内、计算平整机轧制力冲击响应幅值不超负荷,最终结合实际生产情况,确定出平整机轧制力设定值,减少因为设定不准造成的生产不稳定、轧断等不良现象出现。可以提高平整机板形控制稳定性,提高成材率。
附图说明
通过附图对本发明较佳实施例的描述,可以进一步理解本发明的目的、特征和优点。
图1是本发明中平整机轧制力计算及优化设定的程序计算流程框图。
图2是平整机纵向振动简化力学模型。
具体实施方式
以下描述本发明的较佳实施例。
实施例1
为了阐述本发明的基本思想,现以某1850六辊平整机组为例,借助于图1来描述规格为1.2m*1.0mm、钢种为SPCC的带钢在特定的平整机组上的张力优化设定过程。
首先,在步骤1中收集平整机组的设备及工艺参数,在步骤1中收集平整机组的设备及工艺参数,包括工作辊直径550mm,支撑辊、中间辊和工作辊的等效质量5.4*104kg、7.02*103kg、4.71*103kg,平整机牌坊的等效刚度7.5*1010N/m,支承辊与中间辊间接触刚度1.0*1010N/m,工作辊与中间辊间接触刚度0.5*1010N/m,轧制压力最大值12*106N,轧制压力最小值6*105N,进入步骤2;
在步骤2中,整理现场生产数据,将平整机和带钢相关参数代入平整机轧制力计算公式
采用matlab多维拟合函数fminsearch对样本进行拟合,拟合后得到:a0=-11.776573, a1=-4.7147225, k3=0.000580898, a=64153.172, k1=0.3, k2=0.7, k4=-178.83734, k5=9159.3527,k6=30439.591。进入步骤3;
随后,在步骤3中,将待平整带钢的特征参数代入步骤2的公式中,包括带钢宽度B=1.2m、厚度h=1mm、延伸率ε=0.01,抗拉强度σb=360*106Pa;平整机前后张力Q1=Q0=45*103N,平整速度为6m/s,摩擦系数为0.11。计算得到平整机轧制力为3.32*106N。进入步骤4;
随后,在步骤4中,根据公式,计算出平整机咬入角为0.005934弧度(0.34°)。然后将咬入角代入到公式:,得到中性角为约为0.00288弧度(0.165°),进入步骤5。
在步骤5中,根据公式打滑因子fγ=γ/α,计算出打滑因子=0.48,满足“<0.8且>0.2”,进入步骤6;。
在步骤6中,将平整机参数代入振动模型,得到平整机轧制力最大幅值4.3*106N,满足条件“<轧制压力最大值12*106N”,进入步骤,7;
在步骤7中,将轧制力设定值P写入平整机轧制力设定表,并植入控制***。
实施例2
为了进一步阐明专利的设计思想,现以某1850六辊平整机组为例,借助于图1来描述规格为1.7m*0.65mm、钢种为Trip780的带钢在特定的平整机组上的张力优化设定过程。
首先,在步骤1中收集平整机组的设备及工艺参数,在步骤1中收集平整机组的设备及工艺参数,包括工作辊直径550mm,支撑辊、中间辊和工作辊的等效质量5.4*104kg、7.02*103kg、4.71*103kg,平整机牌坊的等效刚度7.5*1010N/m,支承辊与中间辊间接触刚度1.0*1010N/m,工作辊与中间辊间接触刚度0.5*1010N/m,轧制压力最大值12*106N,轧制压力最小值0.6*106N,进入步骤2;
在步骤2中,整理现场生产数据,将平整机和带钢相关参数代入平整机轧制力计算公式
采用matlab多维拟合函数fminsearch对样本进行拟合,拟合后得到,a0=-11.776573, a1=-4.7147225, k3=0.000580898, a=64153.172, k1=0.3, k2=0.7, k4=-178.83734, k5=9159.3527,k6=30439.591。进入步骤3;
随后,在步骤3中,将待平整带钢的特征参数代入步骤2的公式中,包括带钢宽度B=1.7m、厚度h=0.65mm、延伸率ε=0.008,抗拉强度σb=780*106Pa;平整机前后张力Q1=Q0=85*103N,平整速度为6m/s,摩擦系数为0.11。计算得到平整机轧制力为15.9*106N,由于计算得到的轧制力满足条件“>Pmax=12*106N”,因此设定平整机轧制力为12*106N。
随后,在步骤4中,根据公式,计算出平整机咬入角为0.004294弧度(0.246°)。然后将咬入角代入到公式:,得到中性角为约为0.002094弧度(0.12°),进入步骤5;
在步骤5中,根据公式打滑因子fγ=γ/α,计算出打滑因子=0.49,满足<0.8且>0.2,进入步骤6;
在步骤6中,将平整机参数代入振动模型,得到平整机轧制力最大幅值17.3*106N,不满足条件“<轧制压力最大值12*106N”,令平整机轧制力P=P-100000N *106N=11.9*106,因为满足P>0.6*106,进入步骤4;
……
步骤4中,将P=11.9*106N代入到计算公式,最终计算出打滑因子fγ=γ/α=0.49,进入步骤5;
在步骤5中,打滑因子=0.49,满足<0.8且>0.2,进入步骤6;
在步骤6中,将平整机参数代入振动模型,得到平整机轧制力最大幅值17.14*106N,不满足条件“<轧制压力最大值12*106N”。令平整机轧制力P= P-100000N =11.8*106,因为满足P>0.6*106,则进入步骤4;
……
步骤4中,将P=8.3*106N代入到计算公式,最终计算出打滑因子fγ=γ/α=0.49,进入步骤5;
在步骤5中,打滑因子=0.49,满足<0.8且>0.2,进入步骤6;
在步骤6中,将参数代入振动模型,得到平整机轧制力最大冲击响应幅值11.9*106N,满足条件平整机最大轧制力幅值<轧制压力最大值12*106N;进入步骤7;
在步骤7中,将轧制力设定值P=8.3*106N,写入平整机轧制力设定表,并植入控制***。
Claims (1)
1.一种六辊平整机组轧制力优化设定方法,其特征在于,包括以下步骤:
(1)针对所分析平整机组,统计整理设备及工艺参数,包括工作辊直径DW,支撑辊、中间辊和工作辊的等效质量Mb、Mm、Mw,平整机牌坊的等效刚度Km,支承辊与中间辊间接触刚度Kb,工作辊与中间辊间接触刚度Kw,轧制压力最大值Pmax,轧制压力最小值Pmin,进入步骤(2);
(2)根据收集的实际生产数据,以及建立的轧制力计算方程,采用回归拟合的方法,计算出带钢的平整轧制力计算公式;进入步骤(3);
P=fB
σp=k3(σs+alg1000e)-(k1σ0+k2σ1)
式中:P——总轧制力,单位:N;f——单位宽度轧制力,单位:N/m;B——带钢宽度,单位:m;L——轧制变形区中轧辊与带钢接触弧长度,单位:mm;a0——平整钢种影响系数;a1——工况影响系数;D——工作辊直径,单位:mm;ε——带钢延伸率;μ——摩擦系数;h0——入口带钢厚度,单位:mm;σp——当量变形抗力,单位:Pa;e——应变速率;a——应变速率系数;v——轧制速度,单位:m/s;k3——变形抗力影响系数;k2——后张力影响系数;k1——前张力影响系数;σs——带钢屈服强度;单位:Pa;σ0——前张应力,单位:Pa;σ1——后张应力,单位:Pa;k4、k5、k6——平整机速度修正系数;
(3)确定平整机的入口带钢厚度h0,带钢宽度B,带钢平整延伸率ε,带钢屈服强度σs,产线速度v,平整机前后张应力σ0、σ1,根据步骤2中的公式计算出平整机轧制力设定值P,如果计算结果大于Pmax,则轧制力P=Pmax;当计算结果小于Pmin,则轧制力P=Pmin;进入步骤4;
(4)根据公式计算出平整机咬入角,根据公式:计算出中性角,进入步骤(5);
式中:f——单位宽度轧制力;B——带钢宽度,单位:m;Q1,Q0——作用在轧件上的前张力、后张力,单位:N;α——咬入角,单位:弧度;
(5)计算出打滑因子fγ=γ/α,判断是否满足0.2<fγ<0.8?如果满足,则进入步骤(6);如果不满足,则令单位宽度轧制力f=f-sign(γ/α-1/2)*sign(Q1-Q0)*10000N,进入步骤(4),如果循环次数大于1000仍不满足0.2<fγ<0.8,取最后一次迭代循环的单位宽度轧制力f值;进入步骤(6);
(6)根据平整机结构,建立3自由度振动模型,力学模型,采用MATLAB软件求出振动方程的数值解,计算出施加平整机轧制力后的平整机的纵向振动响应,得到振动时轧制力最大幅值Pa,判断Pa<Pmax是否满足?如果满足,则进入步骤(7);如果不满足,则P=P-100000N,如果此时满足P>Pmin,则进入步骤(4);如果P>Pmin条件不满足,则P=Pmin,进入步骤(7);
(7)将完成计算得到的轧制力设定值P写入平整机轧制力设定表,并植入控制***。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310276037.2A CN103302108B (zh) | 2013-07-02 | 2013-07-02 | 一种六辊平整机组轧制力优化设定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310276037.2A CN103302108B (zh) | 2013-07-02 | 2013-07-02 | 一种六辊平整机组轧制力优化设定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103302108A CN103302108A (zh) | 2013-09-18 |
CN103302108B true CN103302108B (zh) | 2015-04-22 |
Family
ID=49128077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310276037.2A Active CN103302108B (zh) | 2013-07-02 | 2013-07-02 | 一种六辊平整机组轧制力优化设定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103302108B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104785539B (zh) * | 2014-01-21 | 2017-01-04 | 宝山钢铁股份有限公司 | 一种用于轧制力调节的张力优化补偿的方法 |
CN104001739B (zh) * | 2014-05-14 | 2015-12-09 | 中冶南方工程技术有限公司 | 一种冷轧平整轧制压力的获取方法及装置 |
CN104128376B (zh) * | 2014-06-23 | 2016-04-20 | 中国重型机械研究院股份公司 | 一种极薄带钢高速精整机组振动计算方法 |
CN106707742A (zh) * | 2015-12-01 | 2017-05-24 | 江苏九天光电科技有限公司 | 适用于镀镍连续生产线转向辊前后张力优化方法 |
CN107008766A (zh) * | 2017-06-14 | 2017-08-04 | 中山中粤马口铁工业有限公司 | 一种单机架四辊平整机生产dr‑ba马口铁基板的方法 |
CN107609301B (zh) * | 2017-09-28 | 2021-03-02 | 中冶南方工程技术有限公司 | 一种利用自适应算法计算带材处理线张力设定的方法 |
CN109590338B (zh) * | 2017-09-30 | 2020-07-28 | 宝山钢铁股份有限公司 | 用于减小二次冷轧间的轧制最小变形量的参数优化方法 |
CN108326049B (zh) * | 2017-12-22 | 2019-06-07 | 中冶南方工程技术有限公司 | 一种冷连轧轧制力的自学习方法 |
CN108733901A (zh) * | 2018-05-02 | 2018-11-02 | 燕山大学 | 一种双平整机组以粗糙度控制为目标的工艺参数优化方法 |
CN110756593B (zh) | 2018-07-26 | 2020-10-27 | 宝山钢铁股份有限公司 | 一种抑制冷连轧机组振动的张力制度优化方法 |
CN109840369B (zh) * | 2019-01-17 | 2020-08-11 | 燕山大学 | 冷轧带材板形仪传动方案设计方法 |
CN112207139B (zh) * | 2020-11-04 | 2023-04-18 | 攀枝花学院 | 六辊平整机轧制力的计算方法 |
CN113649420B (zh) * | 2021-07-28 | 2023-10-24 | 首钢京唐钢铁联合有限责任公司 | 一种平整机轧制力获取方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1923392A (zh) * | 2005-08-30 | 2007-03-07 | 宝山钢铁股份有限公司 | 平整机轧制压力设定方法 |
CN102294364A (zh) * | 2010-06-22 | 2011-12-28 | 宝山钢铁股份有限公司 | 极薄板平整机的轧制力预设定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04238616A (ja) * | 1991-01-08 | 1992-08-26 | Kawasaki Steel Corp | 鋼帯の調質圧延における圧下力制御方法 |
-
2013
- 2013-07-02 CN CN201310276037.2A patent/CN103302108B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1923392A (zh) * | 2005-08-30 | 2007-03-07 | 宝山钢铁股份有限公司 | 平整机轧制压力设定方法 |
CN102294364A (zh) * | 2010-06-22 | 2011-12-28 | 宝山钢铁股份有限公司 | 极薄板平整机的轧制力预设定方法 |
Non-Patent Citations (4)
Title |
---|
单机架平整机平整工艺参数的优化研究;张晓明等;《轧钢》;19991231;第209-212页 * |
宝钢1580热轧两辊平整机轧制力计算模型的研究;丁鸿儒;《宝钢技术》;20031231(第3期);第25-30页 * |
热镀锌机组平整机轧制力预设定模型研究;薛栋梁等;《重型机械》;20101231(第6期);第18-21页 * |
适用于双机架薄板平整机稳态轧制负荷的数学模型;唐成龙等;《上海交通大学学报》;20090131;第43卷(第1期);第19-24页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103302108A (zh) | 2013-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103302108B (zh) | 一种六辊平整机组轧制力优化设定方法 | |
CN103567230B (zh) | 一种微张力控制***和方法 | |
CN102672003B (zh) | 不锈钢带钢拉矫机机组的工艺参数设定方法 | |
CN102688898B (zh) | 冷轧带钢双机架平整机板形控制方法 | |
CN102756013B (zh) | 一种冷轧带钢翘曲的改善方法及装置 | |
CN104889175B (zh) | 一种提高平整稳定性和产品表面质量的张力设定方法 | |
CN102266865A (zh) | 热冷轧负荷分配方法 | |
CN102836882B (zh) | 一种热轧带钢穿带时产生翘头的消除方法 | |
CN103962392B (zh) | 一种热连轧机精轧机组动态负荷控制方法 | |
CN105251780B (zh) | 一种平整机延伸率波动的控制方法 | |
CN105251778B (zh) | 单锥度工作辊窜辊轧机边部减薄反馈控制方法 | |
CN102784814B (zh) | 宽厚金属板材矫直机弯辊补偿方法 | |
CN102989787B (zh) | 一种热轧带钢头尾板形控制方法 | |
CN104985005B (zh) | 板带轧机不对称板形的补偿方法 | |
CN110385343A (zh) | 一种热轧平整机翘曲控制方法 | |
CN105057364A (zh) | 一种镁合金板材轧制边裂预判及控制方法 | |
CN103008361B (zh) | 一种控制钢板头部弯曲的方法 | |
CN104070070A (zh) | 通过张力补偿提升精轧带钢轧制力和厚度精度的综合控制方法 | |
CN110170535B (zh) | 一种镀锌光整机过焊缝轧制力自动控制方法 | |
CN104324939B (zh) | 一种带钢轧制方法 | |
CN106607460A (zh) | 森吉米尔20辊轧机边部板形控制方法及装置 | |
CN102581030B (zh) | 冷轧带钢板板形闭环控制周期的确定方法 | |
CN103230943A (zh) | 改善中厚板头部弯曲的方法 | |
CN110893418B (zh) | 平整机过焊缝的控制方法 | |
CN104550258B (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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100041 Shijingshan Road, Shijingshan District, Shijingshan District, Beijing Patentee after: Shougang Group Co. Ltd. Address before: 100041 Shijingshan Road, Shijingshan District, Shijingshan District, Beijing Patentee before: Capital Iron & Steel General Company |
|
CP01 | Change in the name or title of a patent holder |