CN115062431B - 一种基于CS-Elman神经网络模型的热轧板凸度预测方法 - Google Patents
一种基于CS-Elman神经网络模型的热轧板凸度预测方法 Download PDFInfo
- Publication number
- CN115062431B CN115062431B CN202210734614.7A CN202210734614A CN115062431B CN 115062431 B CN115062431 B CN 115062431B CN 202210734614 A CN202210734614 A CN 202210734614A CN 115062431 B CN115062431 B CN 115062431B
- Authority
- CN
- China
- Prior art keywords
- neural network
- hot rolled
- convexity
- rolled plate
- data
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000003062 neural network model Methods 0.000 title claims abstract description 16
- 238000005096 rolling process Methods 0.000 claims abstract description 36
- 238000004519 manufacturing process Methods 0.000 claims abstract description 21
- 238000013528 artificial neural network Methods 0.000 claims description 85
- 238000012546 transfer Methods 0.000 claims description 18
- 238000010606 normalization Methods 0.000 claims description 16
- 230000002159 abnormal effect Effects 0.000 claims description 15
- 210000002569 neuron Anatomy 0.000 claims description 15
- 241000544061 Cuculus canorus Species 0.000 claims description 14
- 235000013601 eggs Nutrition 0.000 claims description 12
- 238000009499 grossing Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 235000005770 birds nest Nutrition 0.000 claims description 6
- 238000012806 monitoring device Methods 0.000 claims description 6
- 235000005765 wild carrot Nutrition 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000005452 bending Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 4
- 229910000831 Steel Inorganic materials 0.000 description 6
- 239000010959 steel Substances 0.000 description 6
- 230000007547 defect Effects 0.000 description 4
- 238000010845 search algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 241000272177 Cuculiformes Species 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005098 hot rolling Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000011031 large-scale manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013000 roll bending Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B38/00—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product
- B21B38/04—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product for measuring thickness, width, diameter or other transverse dimensions of the product
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B38/00—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product
- B21B2038/004—Measuring scale thickness
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Computer Hardware Design (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mechanical Engineering (AREA)
- Feedback Control In General (AREA)
- Control Of Metal Rolling (AREA)
Abstract
本发明公开了一种基于CS‑Elman神经网络模型的热轧板凸度预测方法,涉及轧制过程自动化生产技术领域。在板带热连轧过程中,板凸度、板厚度及板宽度是衡量板带产品质量和尺寸精度是否合格的重要标准,而三者之间往往互相影响,要实现板凸度的精确控制,需要在板厚度及宽度精确控制的基础上进行。本发明使用Elman算法实现对板厚度和宽度的预报,并采用CS算法优化Elman的各层间权值和阈值,然后将板厚度和宽度也作为输入量,进而实现板凸度的在线预报,具有参数少、操作简单、易实现、寻优能力强、可收敛于全局最优等优点,为提高热轧板凸度的预测精度提供了新方法。
Description
技术领域
本发明公开了一种基于CS-Elman神经网络模型的热轧板凸度预测方法,涉及轧制过程自动化生产技术领域。
背景技术
板形是衡量带钢产品质量的指标之一,通常衡量板形好坏的指标主要是板凸度,有效控制带钢板凸度不但可以防止如楔形等缺陷的出现,而且能够保证带钢平直度。随着用户对于产品质量的要求日益提高,一般的板形控制方法,已经很难满足其要求,因此,根据生产过程数据对带钢出口凸度提前进行准确预测,对于实现带钢板凸度的精准控制具有非常重要的意义。
目前板凸度计算模型通过复杂的数学关系建立,模型极为复杂,而且模型系数往往根据经验并通过反复地调参确定,使得调试人员的工作量大,且效率低下,而且随着轧制设备的更新,其影响因素也在增加,板凸度数学模型的精度已经远远不能满足实际需求,这样的生产状况成为企业急需解决的问题。而国内热连轧生产线大都具有完备的数据采集和存储系统,但现场数据大部分未能进行有效的开采、挖掘、加工,没有产生应有的价值,在数据中挖掘板凸度的规律,能够极大降低时间成本与建模复杂度,并且提升预测精度。
板凸度、板厚度及板宽度作为衡量板带尺寸精度的重要标准,三者之间互相影响较大。对轧制过程的下一卷带钢进行板凸度预测时,由于此时板厚度和宽度未知,传统的预测手段仅将轧制过程特征参数作为输入值,忽略板厚度和宽度这两个重要因素的影响,往往不能实现准确的板凸度预报。
发明内容
针对上述现有技术存在的不足,本发明提供了一种基于CS-Elman神经网络模型的轧制过程热轧板凸度预测方法,该网络能够存储和利用过去时刻输出信息,具有参数少、操作简单、易实现、寻优能力强、可收敛于全局最优等优点。
本发明的方法包括以下步骤:
S1、利用生产现场现有的监测装置采集数据,提取热连轧生产现场监测装置采集的历史数据中,与热轧板凸度、厚度及宽度有关的特征参数,特征参数主要包括各个机架的轧制力、轧制速度、弯辊力、辊缝值、窜辊量,以及轧制的入口温度、出口温度等,还可以包括其他影响热轧板最终板形的特征参数。另外还要提取采用相应特征参数所生产的热轧板的实际凸度、厚度和宽度值。
S2、对采集的数据进行归一化处理。
归一化处理可采用以下公式:
对于某一种数据,比如轧制的入口温度数据,min为采集的所有样本数据中的最小值,max为采集的所有样本数据中的最大值,xi为要进行归一化的第i个样本数据,是第i个数据进行归一化后的映射数值。归一化的效果是原始数据被规整到[0,1]范围内,消除量纲的影响。
在进行数据归一化之前,可以按照以下步骤2.1~2.4中的方法,针对每一种特征参数数据或实际凸度、厚度、宽度值数据,进行以下预处理:
S2.1计算所有样本中某种数据(例如入口温度数据)的平均值与标准差:
所述为该种数据的平均值,Sx为该种数据的标准差,n为该种数据的个数,xi为第i个数据;
S2.2计算每个样本数据与总体平均值之差的绝对值
S2.3将上述绝对值中超过3倍标准差的数据点,即满足的数据点作为异常值予以剔除;
S2.4将剔出异常值后的数据进行平滑处理。
之后再将平滑处理后的数据进行归一化处理。其中S2.4中的平滑处理优选采用5点3次平滑法。
S3、建立热轧板厚度-宽度预报神经网络:
以特征参数为输入层的输入量,以热轧板的厚度和宽度值为输出层,构建Elman神经网络。通过CS算法(布谷鸟搜索算法)搜索最佳的Elman各层的权值和阈值赋给Elman神经网络,完成热轧板厚度-宽度预报神经网络的建立。
S4、建立热轧板凸度预报神经网络:
以特征参数连同对应的热轧板的厚度和宽度值一起作为输入层,以热轧板的凸度为输出层,构建Elman神经网络。通过CS算法搜索最佳的Elman神经网络各层的权值和阈值赋给Elman神经网络,完成热轧板凸度预报神经网络的建立。
S5、利用建立好的热轧板厚度-宽度预报神经网络和热轧板凸度预报神经网络,进行实际热连轧生产过程中热轧板凸度的预测:
S5.1、以热连轧生产中特征参数的实际值为热轧板厚度-宽度预报神经网络的输入量,通过热轧板厚度-宽度预报神经网络获得热轧板厚度和宽度的预测值;
S5.2、以热连轧生产中特征参数的实际值连同S5.1获得的热轧板厚度和宽度的预测值共同作为热轧板凸度预报神经网络的输入量,通过热轧板凸度预报神经网络获得热轧板凸度的预测值。
步骤S3和S4的Elman神经网络的构建,主要需要确定隐含层神经元节点数、隐含层和输出层的传递函数,以及各层中的权值和阈值。
隐含层神经元节点数的可选范围hiddennum可以采用以下公式确定:
hiddennum=sqrt(m+n)+a
其中,m为输入层节点个数,n为输出层节点个数,a取为1-10之间的整数。对于热轧板厚度-宽度预报神经网络,m值为与热轧板凸度、厚度及宽度有关的特征参数的种类数,n值为2,分别对应热轧板厚度和宽度;对于热轧板凸度预报神经网络,m值为特征参数的实际值连同厚度和宽度的预测值,n值为1,对应热轧板凸度。
通过以上公式可以得出隐含层神经元节点数的可选范围为hiddennum个。在确定可选范围后,可以以mae平均绝对值误差、R2决定系数为误差指标,遍历隐含层神经元节点数的可选范围,找到最佳隐含神经元的数量用于Elman神经网络。
关于步骤S3和S4的Elman神经网络中隐含层和输出层的传递函数,可以从purelin、tansig两种传递函数中选择,进而隐含层-输出层传递函数可以有4种不同组合,可通过遍历所有传递函数组合,找到最适合Elman神经网络模型的神经网络隐含层和输出层的传递函数,用于Elman神经网络。
如前文所述,通过CS算法(布谷鸟搜索算法)搜索最佳的Elman各层的权值和阈值赋给Elman神经网络,具体方法包括以下步骤:
S3.1、初始化CS算法的参数;
1)对每个布谷鸟数目初始化参数,在设定范围内遍历每个布谷鸟数目,找出当前最佳布谷鸟数目;
2)对察觉异常鸟蛋的概率初始化参数,在设定范围内遍历察觉异常鸟蛋的概率,找出当前最佳察觉异常鸟蛋的概率;
3)对最大进化代数初始化参数,在设定范围内遍历每个最大进化代数,找出当前最佳最大进化代数;
S3.2、随机选取一组鸟窝对应Elman神经网络的各层权值和阈值;
S3.3、将Elman神经网络预测值与数据集中实际数据结果的均方误差设置为适应度函数;
显然所述实际数据结果对于步骤S3中的热轧板厚度-宽度预报神经网络而言是热轧板实际厚度和宽度,对于步骤S4中的热轧板凸度预报神经网络而言是热轧板实际凸度;
S3.4、计算当前的最优适应度值,得到当前最优的鸟窝位置;
S3.5、开始进化,通过Lexy飞行产生新解,更新适应度较优的鸟窝,根据发现概率发现鸟蛋,舍弃糟糕的鸟窝;
S3.6、列出当前最佳鸟窝,得到最佳适应度;
S3.7、将最佳适应度下的各层权值和阈值赋给Elman神经网络。
本发明具有的优点是:
本发明利用CS算法(布谷鸟搜索算法)结合Elman神经网络预测板凸度,具有参数少、操作简单、易实现、寻优能力强、可收敛于全局最优等优点。CS算法的寻优过程通过莱维飞行大步长小步长交替使用的方法,其中少量的大步长可以扩大搜索范围,有利于跳出局部最优;而小步长有利于局部寻优,提高算法精度,解决了Elman网络容易陷入局部最优的缺点,进而克服了热轧生产过程中参数检测困难精度差的缺陷,不仅精度高,而且运算速度快,利用现场采集的大量生产数据通过直接在计算机上编程,即可实现投入使用,成本较低廉,推广能力较强。
本发明采用热轧板厚度-宽度预报和热轧板凸度预报的两个神经网络,先通过轧制过程中的轧制规程和实时测定值中的各种特征参数预测热轧板的厚度和宽度值,然后将预测板凸度时不能实时测定的厚度和宽度两个参数也作为预测凸度的条件,增强了预测的精确性。
附图说明
图1为本发明实施例中对部分数据进行异常值剔除的对比图,左图为剔除之前,右图为剔除之后。
图2为本发明实施例中预测值与实际测量值的散点图。
图3为本发明实施例中预测误差柱状分布图。
具体实施方式
为了使本发明的目的、技术方案及优势更加清晰,下面结合附图和具体实施例对本发明做进一步详细说明。此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本实施例中的热连轧过程具有7个道次,相应有7个机架,对该热连轧过程的热轧板凸度,采用以下方法进行预测:
S1、利用生产现场现有的监测装置采集数据,提取热连轧生产现场监测装置采集的历史数据中,与热轧板凸度、厚度及宽度有关的特征参数,充分考虑热连轧的轧制机理,提取出49个相关的特征参数,并通过相关性分析对比后,筛选得到其中37个作为特征参数,也即神经网络的输入特征,包括:F1~F7机架的轧制力、F1~F7机架轧制速度、F1~F7机架弯辊力、F1~F7机架辊缝值、F1~F7机架的窜辊量、轧制过程的入口温度、轧制过程的出口温度。
另外还要提取采用相应特征参数所生产的热轧板的实际凸度、厚度和宽度值。
S2、对采集的数据(数据集)进行归一化处理。
首先针对每一种特征参数数据或实际凸度、厚度、宽度值数据,进行以下预处理:
S2.1计算所有样本中某一种数据的平均值与标准差:
所述为该种数据的平均值,Sx为该种数据的标准差,n为该种数据的个数,xi为第i个数据;
S2.2计算每个样本数据与总体平均值之差的绝对值
S2.3将上述绝对值中超过3倍标准差的数据点,即满足的数据点作为异常值予以剔除;异常值剔除前后的对比图如图1所示;
S2.4将剔出异常值后的数据进行平滑处理,具体采用5点3次平滑法。
之后再将平滑处理后的数据进行归一化处理,归一化处理可采用以下公式:
对于某一种数据,比如轧制的入口温度数据,min为采集的所有样本数据中的最小值,max为采集的所有样本数据中的最大值,xi为要进行归一化的第i个样本数据,是第i个数据进行归一化后的映射数值。归一化后,原始数据被规整到[0,1]范围内,消除了量纲的影响。
S3、建立热轧板厚度-宽度预报神经网络:
以特征参数为输入层的输入量,以热轧板的厚度和宽度值为输出层,构建Elman神经网络。
热轧板厚度-宽度预报神经网络的隐含层神经元节点数的可选范围hiddennum采用以下公式确定:
hiddennum=sqrt(m+n)+a
其中,m为输入层节点个数,本实施例中对于热轧板厚度-宽度预报神经网络,m值为与热轧板凸度、厚度及宽度有关的特征参数的种类数,根据上文可知共37个,n值为2,分别对应热轧板厚度和宽度;可以得出热轧板厚度-宽度预报神经网络的隐含层神经元节点数可取8~16个,以MAE平均绝对值误差、R2决定系数为误差指标,遍历隐含层神经元节点数的可选范围,找到最佳隐含神经元的数量为8个,用于Elman神经网络。
热轧板厚度-宽度预报神经网络中隐含层和输出层的传递函数,可以从purelin、tansig两种传递函数中选择,可以有purelin-purelin、purelin-tansig、tansig-purelin、tansig-tansig四种组合,遍历该4种组合,找到最适合本模型的神经网络隐含层和输出层的传递函数,找到最优的传递函数组合为tansig-purelin(隐含层-输出层)。
通过CS算法(布谷鸟搜索算法)搜索最佳的Elman各层的权值和阈值赋给Elman神经网络,完成热轧板厚度-宽度预报神经网络的建立,具体方法为:
S3.1、初始化CS算法的参数;
1)对每个布谷鸟数目初始化参数,设定布谷鸟数目范围为5到30个,间隔5个为一组,遍历每组布谷鸟数目,找出当前最佳布谷鸟数目;
2)对察觉异常鸟蛋的概率初始化参数,设察觉异常鸟蛋的概率为0.1到0.9,间隔0.2为一组,遍历察觉异常鸟蛋的概率,在设定范围内遍历察觉异常鸟蛋的概率,找出当前最佳察觉异常鸟蛋的概率;
3)对最大进化代数初始化参数,设最大进化代数30到100个,间隔10个为一组,在设定范围内遍历每个最大进化代数,找出当前最佳最大进化代数;
S3.2、随机选取一组鸟窝对应Elman神经网络的各层权值和阈值;
S3.3、将Elman神经网络预测值与数据集中实际数据结果的均方误差设置为适应度函数;
S3.4、计算当前的最优适应度值,得到当前最优的鸟窝位置;
S3.5、开始进化,通过Lexy飞行产生新解,更新适应度较优的鸟窝,根据发现概率发现鸟蛋,舍弃糟糕的鸟窝;
S3.6、列出当前最佳鸟窝,得到最佳适应度;
S3.7、将最佳适应度下的各层权值和阈值赋给Elman神经网络。
S4、建立热轧板凸度预报神经网络:
以特征参数连同热轧板的厚度和宽度值一起作为输入层,以热轧板的凸度为输出层,构建Elman神经网络。
热轧板凸度预报神经网络的隐含层神经元节点数的可选范围hiddennum的确定方法与热轧板厚度-宽度预报神经网络类似,采用以下公式确定:
hiddennum=sqrt(m+n)+a
其中,m为输入层节点个数,本实施例中,对于热轧板凸度预报神经网络,m值为39,分别对应37个特征参数,以及热轧板厚度和宽度,n值为1,对应热轧板凸度。可以得出热轧板凸度预报神经网络的隐含层神经元节点数可取8~16个,以MAE平均绝对值误差、R2决定系数为误差指标,遍历隐含层神经元节点数的可选范围,找到最佳隐含神经元的数量为8个,用于Elman神经网络。
热轧板厚度-宽度预报神经网络中隐含层和输出层的传递函数的确定方法与步骤S3中类似,隐含层-输出层传递函数组合为tansig-purelin。
通过CS算法搜索最佳的Elman神经网络各层的权值和阈值赋给Elman神经网络,完成热轧板凸度预报神经网络的建立。具体步骤与步骤S3中类似。
S5、利用建立好的热轧板厚度-宽度预报神经网络和热轧板凸度预报神经网络,进行实际热连轧生产过程中热轧板凸度的预测:
S5.1、以热连轧生产中特征参数的实际值为热轧板厚度-宽度预报神经网络的输入量,通过热轧板厚度-宽度预报神经网络获得热轧板厚度和宽度的预测值;
S5.2、以热连轧生产中特征参数的实际值连同S5.1获得的热轧板厚度和宽度的预测值共同作为热轧板凸度预报神经网络的输入量,通过热轧板凸度预报神经网络获得热轧板凸度的预测值。
多次预测获得的预测值与实际测量值的差距如图2所示,预测绝对误差分布如图3所示。结果显示98.66%的预测凸度绝对误差小于5μm,根据工业生产的要求,实际凸度值在设定值±10μm内视为产品合格,因此本发明提出的预测模型符合要求。
上面仅对本发明的较佳实施例作了详细说明,但是本发明并不限于上述实施例,对于其他实施例,其依然可以对本发明所展示的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换,因而这些修改、替换均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,包括下列步骤:
S1、利用生产现场现有的监测装置采集数据,提取热连轧生产现场监测装置采集的历史数据中,与热轧板凸度、厚度及宽度有关的特征参数,以及相应热轧板的轧后实际凸度、厚度和宽度值;
S2、对采集的数据进行归一化处理;
S3、建立热轧板厚度-宽度预报神经网络:
以特征参数为输入层的输入量,以热轧板的厚度和宽度值为输出层,构建Elman神经网络;通过CS算法搜索最佳的Elman各层的权值和阈值赋给Elman神经网络,完成热轧板厚度-宽度预报神经网络的建立;
S4、建立热轧板凸度预报神经网络:
以特征参数连同热轧板的厚度和宽度值一起作为输入层,以热轧板的凸度为输出层,构建Elman神经网络;
通过CS算法搜索最佳的Elman各层的权值和阈值赋给Elman神经网络,完成热轧板凸度预报神经网络的建立;
S5、利用建立好的热轧板厚度-宽度预报神经网络和热轧板凸度预报神经网络,进行实际热连轧生产过程中热轧板凸度的预测:
S5.1、以热连轧生产中特征参数的实际值为热轧板厚度-宽度预报神经网络的输入量,通过热轧板厚度-宽度预报神经网络获得热轧板厚度和宽度的预测值;
S5.2、以热连轧生产中特征参数的实际值连同S5.1获得的热轧板厚度和宽度的预测值共同作为热轧板凸度预报神经网络的输入量,通过热轧板凸度预报神经网络获得热轧板凸度的预测值;
所述与热轧板凸度、厚度及宽度有关的特征参数,包括各个机架的轧制力、轧制速度、弯辊力、辊缝值、窜辊量,以及轧制的入口温度、出口温度。
2.根据权利要求1所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,所述步骤S2中,在进行归一化之前,针对每一种特征参数数据或实际凸度、厚度、宽度值数据,进行以下步骤S2.1~2.4的预处理:
S2.1、计算所有样本中该种数据的平均值与标准差:
所述为该种数据的平均值,Sx为该种数据的标准差,n为该种数据的个数,xi为第i个数据;
S2.2、计算每个样本数据与总体平均值之差的绝对值
S2.3、将上述绝对值中超过3倍标准差的数据点,即满足的数据点作为异常值予以剔除;
S2.4、将剔出异常值后的数据进行平滑处理;
之后再将平滑处理后的数据进行归一化处理。
3.根据权利要求2所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,所述S2.4中的平滑处理采用5点3次平滑法。
4.根据权利要求1所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,所述S2中归一化处理采用以下公式:
min为数据的最小值,max为数据的最大值,xi为要进行归一化的第i个数据,是第i个数据进行归一化后的映射数值。
5.根据权利要求1所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,所述步骤S3和S4的Elman神经网络中,隐含层神经元节点数的可选范围hiddennum采用以下公式确定:
hiddennum=sqrt(m+n)+a
其中,m为输入层节点个数,n为输出层节点个数,a取为1-10之间的整数。
6.根据权利要求5所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,隐含层神经元节点数的可选范围为hiddennum个,以mae平均绝对值误差、R2决定系数为误差指标,遍历隐含层神经元节点数的可选范围,找到最佳隐含神经元的数量用于Elman神经网络。
7.根据权利要求1所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,所述步骤S3和S4的Elman神经网络中,隐含层和输出层的传递函数从purelin、tansig两种传递函数中选择。
8.根据权利要求7所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,遍历所有传递函数组合,找到最适合Elman神经网络模型的神经网络隐含层和输出层的传递函数,用于Elman神经网络。
9.根据权利要求1所述的基于CS-Elman神经网络模型的热轧板凸度预测方法,其特征在于,所述步骤S3和S4中,通过CS算法搜索最佳的Elman神经网络各层的权值和阈值赋给Elman神经网络,其采用的方法包括以下步骤:
S3.1、初始化CS算法的参数;
1)对每个布谷鸟数目初始化参数,在设定范围内遍历每个布谷鸟数目,找出当前最佳布谷鸟数目;
2)对察觉异常鸟蛋的概率初始化参数,在设定范围内遍历察觉异常鸟蛋的概率,找出当前最佳察觉异常鸟蛋的概率;
3)对最大进化代数初始化参数,在设定范围内遍历每个最大进化代数,找出当前最佳最大进化代数;
S3.2、随机选取一组鸟窝对应Elman神经网络的各层权值和阈值;
S3.3、将Elman神经网络预测值与数据集中实际数据结果的均方误差设置为适应度函数;S3.4、计算当前的最优适应度值,得到当前最优的鸟窝位置;
S3.5、开始进化,通过Lexy飞行产生新解,更新适应度较优的鸟窝,根据发现概率发现鸟蛋,舍弃糟糕的鸟窝;
S3.6、列出当前最佳鸟窝,得到最佳适应度;
S3.7、将最佳适应度下的各层权值和阈值赋给Elman神经网络。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210734614.7A CN115062431B (zh) | 2022-06-27 | 2022-06-27 | 一种基于CS-Elman神经网络模型的热轧板凸度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210734614.7A CN115062431B (zh) | 2022-06-27 | 2022-06-27 | 一种基于CS-Elman神经网络模型的热轧板凸度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115062431A CN115062431A (zh) | 2022-09-16 |
CN115062431B true CN115062431B (zh) | 2024-05-31 |
Family
ID=83201533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210734614.7A Active CN115062431B (zh) | 2022-06-27 | 2022-06-27 | 一种基于CS-Elman神经网络模型的热轧板凸度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115062431B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20040043911A (ko) * | 2002-11-20 | 2004-05-27 | 주식회사 포스코 | 신경 회로망을 이용한 압연 하중 예측 장치 및 방법 |
CN108637020A (zh) * | 2018-05-09 | 2018-10-12 | 北京科技大学 | 一种自适应变异pso-bp神经网络带钢凸度预测方法 |
CN109558677A (zh) * | 2018-11-29 | 2019-04-02 | 东北大学 | 一种基于数据驱动的热轧板凸度预测方法 |
CN110276483A (zh) * | 2019-06-04 | 2019-09-24 | 湖北大学 | 基于神经网络的制糖原材料预测方法 |
CN110929347A (zh) * | 2019-10-25 | 2020-03-27 | 东北大学 | 一种基于梯度提升树模型的热连轧带钢凸度预测方法 |
CN111209967A (zh) * | 2020-01-08 | 2020-05-29 | 太原科技大学 | 一种基于支持向量机的轧制过程板凸度预测方法 |
CN114021290A (zh) * | 2022-01-04 | 2022-02-08 | 东北大学 | 一种基于数据驱动和机理模型融合的板带钢凸度预测方法 |
CN114091352A (zh) * | 2021-12-01 | 2022-02-25 | 东北大学 | 基于Elman神经网络的热连轧出口板凸度动态预测方法 |
-
2022
- 2022-06-27 CN CN202210734614.7A patent/CN115062431B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20040043911A (ko) * | 2002-11-20 | 2004-05-27 | 주식회사 포스코 | 신경 회로망을 이용한 압연 하중 예측 장치 및 방법 |
CN108637020A (zh) * | 2018-05-09 | 2018-10-12 | 北京科技大学 | 一种自适应变异pso-bp神经网络带钢凸度预测方法 |
CN109558677A (zh) * | 2018-11-29 | 2019-04-02 | 东北大学 | 一种基于数据驱动的热轧板凸度预测方法 |
CN110276483A (zh) * | 2019-06-04 | 2019-09-24 | 湖北大学 | 基于神经网络的制糖原材料预测方法 |
CN110929347A (zh) * | 2019-10-25 | 2020-03-27 | 东北大学 | 一种基于梯度提升树模型的热连轧带钢凸度预测方法 |
CN111209967A (zh) * | 2020-01-08 | 2020-05-29 | 太原科技大学 | 一种基于支持向量机的轧制过程板凸度预测方法 |
CN114091352A (zh) * | 2021-12-01 | 2022-02-25 | 东北大学 | 基于Elman神经网络的热连轧出口板凸度动态预测方法 |
CN114021290A (zh) * | 2022-01-04 | 2022-02-08 | 东北大学 | 一种基于数据驱动和机理模型融合的板带钢凸度预测方法 |
Non-Patent Citations (1)
Title |
---|
遗传算法优化神经网络的热轧带钢弯辊力预报模型;王振华;龚殿尧;李广焘;张殿华;;东北大学学报(自然科学版);20181215(12);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115062431A (zh) | 2022-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111241657B (zh) | 一种基于深度学习的冷轧带钢板形预测方法 | |
CN109992921B (zh) | 一种燃煤电厂锅炉热效率的在线软测量方法及系统 | |
CN108335487B (zh) | 一种基于交通状态时序性的道路交通状态预测系统 | |
CN109558677A (zh) | 一种基于数据驱动的热轧板凸度预测方法 | |
CN109919356B (zh) | 一种基于bp神经网络区间需水预测方法 | |
CN110674858B (zh) | 一种基于时空关联与大数据挖掘的交通舆情检测方法 | |
CN111639823A (zh) | 一种基于特征集构建的建筑冷热负荷预测方法 | |
CN108961460B (zh) | 基于稀疏esgp与多目标优化的故障预测方法及装置 | |
CN115856204B (zh) | 基于三维回声状态网络的掘进工作面瓦斯浓度预测方法 | |
CN115841004B (zh) | 基于多维数据的带钢热轧过程力学性能软测量方法及装置 | |
CN116088453A (zh) | 生产质量预测模型训练方法及装置、生产质量监控方法 | |
CN115438897A (zh) | 一种基于blstm神经网络的工业过程产品质量预测方法 | |
CN115062431B (zh) | 一种基于CS-Elman神经网络模型的热轧板凸度预测方法 | |
CN114548494A (zh) | 一种可视化造价数据预测智能分析系统 | |
CN117851908A (zh) | 一种改进的在线低压台区电能表失准监测方法及装置 | |
CN115374858B (zh) | 基于混合集成模型的流程工业生产品质的智能诊断方法 | |
CN115985411A (zh) | 基于高斯过程回归模型的高分子熔融指数软测量方法 | |
CN115017818A (zh) | 基于注意力机制和多层lstm的电厂烟气含氧量智能预测方法 | |
CN114943328A (zh) | 基于bp神经网络非线性组合的sarima-gru时序预测模型 | |
CN111258996B (zh) | 一种面向工业大数据的产品质量多源深度融合预报方法 | |
CN113312587A (zh) | 基于arima预测和回归预测的传感器采集数据缺失值处理方法 | |
CN113112166A (zh) | 基于灰色模糊层次分析的设备状态变量选择方法及设备 | |
CN112816211A (zh) | 一种基于pso-bp算法的带式输送机故障诊断的方法 | |
CN114798763A (zh) | 一种加热炉出钢板坯粗轧出口温度的预测方法及系统 | |
CN111650894A (zh) | 一种基于隐变量的贝叶斯网络复杂工业过程软测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |