CN104741387A - 一种基于动态规划法和粒子群算法的轧制规程优化方法 - Google Patents
一种基于动态规划法和粒子群算法的轧制规程优化方法 Download PDFInfo
- Publication number
- CN104741387A CN104741387A CN201510091232.7A CN201510091232A CN104741387A CN 104741387 A CN104741387 A CN 104741387A CN 201510091232 A CN201510091232 A CN 201510091232A CN 104741387 A CN104741387 A CN 104741387A
- Authority
- CN
- China
- Prior art keywords
- rolling
- met
- dynamic programming
- schedule
- calculated
- 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.)
- Pending
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/16—Control of thickness, width, diameter or other transverse dimensions
- B21B37/18—Automatic gauge control
- B21B37/20—Automatic gauge control in tandem mills
-
- 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)将各道次的压下量在的范围内随机离散;(4)粒子群空间寻优;(5)当程序满足终止条件时输出最优解;(6)判断是否满足约束条件,如满足则继续下一步,如不满足则返回到第(3)步,直至满足约束条件;(7)运用动态规划法求各机架压下量、决策变量,由状态转移方程求状态变量;(8)求目标函数值并进行比较,找出最小目标函数并计算对应的最优压下规程及优化后的板凸度与平直度;(9)输出并存储最优轧制规程。
Description
技术领域
本发明属于板带材热连轧技术领域,具体涉及一种基于动态规划法和粒子群算法的轧制规程优化方法。
背景技术
轧制规程是轧制生产过程的基础。一个合理的轧制规程是实现高效、优质生产的基础。轧制规程的优化也一直是业内研究的热点问题。
传统的连轧生产轧制规程的优化主要是指对各个机架压下量与轧制力(功率)的合理分配。优化的目的是在轧制设备的力能、速度等参数满足一定约束条件的情况下,使轧制过程与产品质量处于最佳状态。然而,在板带材热连轧实际生产中轧制规程的优化技术也存在一些问题。其中最主要的就是由于控制参数众多,在多个优化目标与约束条件下,计算量大、耗时多、在线适时应用困难。而且由于板形理论的不完善和控制手段的耦合与干扰,使得板形板厚协调控制的精度不高,产品质量有待进一步提高。
我国大型板带热连轧机轧制规程的数学模型大多数是从国外引进的,总计在100套以上。例如武钢1700热连轧机是从日本引进的,轧制规程的设定采用的是能耗分配系数法;鞍钢1700热连轧机和宝钢2050热连轧机使用的是西门子的压下率分配系数法;本钢1700热连轧机则采用AEG提供的单调递减函数法;攀钢1450热连轧机最初是使用意大利ANSALDO公司的能耗法,厚度分布计算采用迭代计算法。攀钢在后期技术改造中,对精轧过程数学模型进行了全面升级改造,负荷分配采用了标准压下率负荷分配法。上述从国外引进的数学模型在生产中基本能够正常使用,但由于需要进行复杂的控制计算,无法满足板形板厚在线协调控制的要求。
国内外在轧制规程优化方面做了大量研究工作,主要有最小平方和目标函数的非线性规化和动态规化方法。但这些方法计算量大,而且只能对分配系数做一些改进,很难实现板形板厚在线适时控制的目的。
发明内容
本发明提出一种基于动态规划法和粒子群算法相结合的轧制规程优化方法,能够有效提高优化运算速度,实现对板带热连轧机轧制规程按照板形板厚协调控制的要求进行在线优化的方法。
本发明的特征在于其操作步骤如下:
(1)采集现有轧制规程以及相关的轧机参数、轧制力数据并分类整理保存在Excel数据文档中;
(2)调取现有轧制规程数据,计算该规程对应的板凸度与平直度;
(3)离散压下量:在保持第一机架入口厚度及最末机架出口厚度不变的前提下,将各道次的压下量在±10%的范围内随机离散;
(4)粒子群空间寻优:在粒子的每一次运动中计算出对应的满足二次型目标函数要求的板凸度、平直度与轧制力参数,并不断更新自己的位置寻优;
(5)当程序满足终止条件时输出最优解,终止条件有二:一是寻找到满足目标条件的最优解,二是达到程序设定的最大迭代次数;
(6)判断是否满足约束条件,如满足则继续下一步,如不满足则返回到第(3)步,直至满足约束条件,约束条件包括咬入条件、压下率约束条件、力能参数约束条件、等比例凸度条件;
(7)运用动态规划法求各机架压下量、决策变量,由状态转移方程求状态变量;
(8)求目标函数值并进行比较,找出最小目标函数并计算对应的最优压下规程及优化后的板凸度与平直度;
(9)输出并存储最优轧制规程。
具体实施方式
为了阐述本发明专利的具体实施方法,以国内某钢铁公司7机架板带热连轧生产机组为例,介绍应用本方法完成轧制规程优化的具体实施步骤。
(1)采集现有轧制规程以及相关的轧机参数、轧制力等数据并分类整理保存在Excel数据文档中。表1为某两个钢卷对应的现有轧制规程及有关数据;
表1
(2)调取现有轧制规程数据,计算该规程对应的板凸度与平直度,表2为初始板凸度与平直度;
表2
(3)离散压下量:在保持第一机架入口厚度及最末机架出口厚度不变的前提下,将各道次的压下量在±10%的范围内随机离散,表3为压下量离散后各道次压下量和出口厚度取值范围;
表3
(4)粒子群空间寻优:在粒子的每一次运动中计算出对应的满足二次型目标函数要求的板凸度、平直度与轧制力等参数,并不断更新自己的位置寻优;
(5)当程序满足终止条件时输出最优解,终止条件有二:一是寻找到满足目标条件的最优解,二是达到程序设定的最大迭代次数;
(6)判断是否满足约束条件,如满足则继续下一步,如不满足则返回到第(2)步,直至满足约束条件;
(7)运用动态规划方法求出各机架压下量、决策变量,由状态转移方程求状态变量;
(8)求目标函数值并进行比较,找出在最小目标函数并计算对应的最优压下规程以及优化后的板凸度与平直度;
(9)输出并存储最优轧制规程,表4为优化后的轧制规程。
表4
对比表1、表2与表4的相关数据可以看出,轧制规程优化后,轧制力的最大变化量为6.82%,不会对设备带来安全隐患。在满足出口厚度精度的前提下,实例中两卷钢卷成品机架的出口板凸度分别由63.73μm、57.57μm降至42.82μm、38.92μm,出口平直度分别由31.72I、27.82I降至10.81I、14.06I,优化效果明显,表5为优化效果。
表5
钢卷号 | R3063504100 | R3049705600 |
目标厚度/mm | 1.83 | 2.00 |
出口厚度/mm | 1.83 | 2.00 |
优化前板凸度/μm | 63.73 | 57.57 |
优化后板凸度/μm | 42.82 | 38.92 |
优化前平直度/I | 31.72 | 27.82 |
优化后平直度/I | 10.81 | 14.06 |
板凸度优化率/% | 32.8 | 32.4 |
平直度优化率/% | 65.9 | 49.5 |
。
Claims (1)
1.一种基于动态规划法和粒子群算法的轧制规程优化方法,其特征在于操作步骤如下:
(1)采集现有轧制规程以及相关的轧机参数、轧制力数据并分类整理保存在Excel数据文档中;
(2)调取现有轧制规程数据,计算该规程对应的板凸度与平直度;
(3)离散压下量:在保持第一机架入口厚度及最末机架出口厚度不变的前提下,将各道次的压下量在的范围内随机离散;
(4)粒子群空间寻优:在粒子的每一次运动中计算出对应的满足二次型目标函数要求的板凸度、平直度与轧制力参数,并不断更新自己的位置寻优;
(5)当程序满足终止条件时输出最优解,终止条件有二:一是寻找到满足目标条件的最优解,二是达到程序设定的最大迭代次数;
(6)判断是否满足约束条件,如满足则继续下一步,如不满足则返回到第(3)步,直至满足约束条件,约束条件包括咬入条件、压下率约束条件、力能参数约束条件、等比例凸度条件;
(7)运用动态规划方法求各机架压下量、决策变量,由状态转移方程求状态变量;
(8)求目标函数值并进行比较,找出最小目标函数并计算对应的最优压下规程及优化后的板凸度与平直度;
(9)输出并存储最优轧制规程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510091232.7A CN104741387A (zh) | 2015-02-28 | 2015-02-28 | 一种基于动态规划法和粒子群算法的轧制规程优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510091232.7A CN104741387A (zh) | 2015-02-28 | 2015-02-28 | 一种基于动态规划法和粒子群算法的轧制规程优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104741387A true CN104741387A (zh) | 2015-07-01 |
Family
ID=53581967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510091232.7A Pending CN104741387A (zh) | 2015-02-28 | 2015-02-28 | 一种基于动态规划法和粒子群算法的轧制规程优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104741387A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110722007A (zh) * | 2019-10-17 | 2020-01-24 | 北京交通大学 | 获取精轧机出口厚度的方法 |
US11938528B2 (en) | 2018-07-19 | 2024-03-26 | Sms Group Gmbh | Method for ascertaining control variables for active profile and flatness control elements for a rolling stand and profile and average flatness values for hot-rolled metal strip |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH062288B2 (ja) * | 1982-10-20 | 1994-01-12 | 新日本製鐵株式会社 | 圧延機の設定方法 |
JPH11104718A (ja) * | 1997-09-29 | 1999-04-20 | Sumitomo Metal Ind Ltd | 可逆式圧延機の圧延方法 |
CN101332473A (zh) * | 2007-06-28 | 2008-12-31 | 株式会社日立制作所 | 热轧轧机的凸度控制装置以及方法 |
CN101890433A (zh) * | 2010-01-30 | 2010-11-24 | 太原科技大学 | 一种板带热连轧机组板形板厚在线协调控制方法 |
-
2015
- 2015-02-28 CN CN201510091232.7A patent/CN104741387A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH062288B2 (ja) * | 1982-10-20 | 1994-01-12 | 新日本製鐵株式会社 | 圧延機の設定方法 |
JPH11104718A (ja) * | 1997-09-29 | 1999-04-20 | Sumitomo Metal Ind Ltd | 可逆式圧延機の圧延方法 |
CN101332473A (zh) * | 2007-06-28 | 2008-12-31 | 株式会社日立制作所 | 热轧轧机的凸度控制装置以及方法 |
CN101890433A (zh) * | 2010-01-30 | 2010-11-24 | 太原科技大学 | 一种板带热连轧机组板形板厚在线协调控制方法 |
Non-Patent Citations (1)
Title |
---|
王强: "基于动态规划与粒子群算法的轧制规程优化方法的研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11938528B2 (en) | 2018-07-19 | 2024-03-26 | Sms Group Gmbh | Method for ascertaining control variables for active profile and flatness control elements for a rolling stand and profile and average flatness values for hot-rolled metal strip |
CN110722007A (zh) * | 2019-10-17 | 2020-01-24 | 北京交通大学 | 获取精轧机出口厚度的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021189739A1 (zh) | 一种基于生产计划的钢铁企业氧气负荷预测方法 | |
CN105013832A (zh) | 一种兼顾轧制能耗与良好板形的热轧带钢负荷分配方法 | |
CN103028615A (zh) | 一种预测带钢热连轧过程温度演变的方法 | |
Peng et al. | Strip shape modeling and its setup strategy in hot strip mill process | |
CN103962390A (zh) | Vc辊平整机湿平整过程中张力与轧制压力综合设定方法 | |
CN105204333B (zh) | 一种提高钢铁企业能源利用率的能耗预测方法 | |
CN103372575B (zh) | 一种酸轧轧机负荷分配方法 | |
CN104375478A (zh) | 一种轧钢生产过程产品质量在线预测与优化的方法及装置 | |
CN105268748B (zh) | 按产品分类的热轧负荷分配方法 | |
CN103439926B (zh) | 一种钢铁企业煤气优化调度装置 | |
CN103418618A (zh) | 冷连轧机性能前馈厚度控制方法 | |
CN103586286A (zh) | 冷连轧机组以拉毛防治为目标的轧制规程综合优化方法 | |
CN102266865A (zh) | 热冷轧负荷分配方法 | |
CN103544340A (zh) | 五机架冷连轧机组极薄带轧制中乳化液浓度的设定方法 | |
CN1091008C (zh) | 基于板形板厚协调规律的板带轧制过程互联控制方法 | |
CN104741387A (zh) | 一种基于动态规划法和粒子群算法的轧制规程优化方法 | |
CN107716560B (zh) | 一种板带热连轧负荷分配方法 | |
CN102489524A (zh) | 一种降低热轧带钢轧制过程能耗的机架负荷分配方法 | |
CN104338753A (zh) | 一种冷连轧机的动态变规格控制方法 | |
CN106777594B (zh) | 一种基于连续曲面的轧制过程数学模型的自学习方法 | |
Bu et al. | Rolling schedule multi-objective optimizationbased on influence function for thin gauge steel strip in tandem cold rolling | |
CN104214772A (zh) | 一种循环流化床机组agc负荷指令响应的控制方法 | |
Huang et al. | Development of a geometric modelling strategy for roll pass optimal design | |
CN104209345B (zh) | 连轧机孔型设计方法 | |
CN105512823A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150701 |