CN116522747A - 一种两阶段优化的挤压铸造工艺参数优化设计方法 - Google Patents
一种两阶段优化的挤压铸造工艺参数优化设计方法 Download PDFInfo
- Publication number
- CN116522747A CN116522747A CN202310180272.3A CN202310180272A CN116522747A CN 116522747 A CN116522747 A CN 116522747A CN 202310180272 A CN202310180272 A CN 202310180272A CN 116522747 A CN116522747 A CN 116522747A
- Authority
- CN
- China
- Prior art keywords
- neural network
- extrusion casting
- casting process
- parameter
- optimization
- 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
- 238000005266 casting Methods 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 93
- 238000005457 optimization Methods 0.000 title claims abstract description 77
- 238000001125 extrusion Methods 0.000 title claims abstract description 65
- 238000013461 design Methods 0.000 title claims abstract description 31
- 241000287127 Passeridae Species 0.000 claims abstract description 69
- 230000008569 process Effects 0.000 claims abstract description 47
- 238000013528 artificial neural network Methods 0.000 claims abstract description 38
- 238000009716 squeeze casting Methods 0.000 claims abstract description 9
- 238000003062 neural network model Methods 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 14
- 230000002431 foraging effect Effects 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 230000000739 chaotic effect Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 210000002569 neuron Anatomy 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 244000062645 predators Species 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 235000003642 hunger Nutrition 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000006872 improvement Effects 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 230000037351 starvation Effects 0.000 claims 1
- 238000010845 search algorithm Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Classifications
-
- 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
-
- 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/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- 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)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Geometry (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种两阶段优化的挤压铸造工艺参数优化设计方法,属于挤压铸造工艺参数优化设计领域,以BP神经网络建立挤压铸造工艺参数与性能(质量)指标的关系模型,引入信息熵权法,确定各性能指标对挤压铸造综合质量的影响权重。考虑到传统BP神经网络的低收敛速度、局部最小化以及现有的智能优化算法的局限性,首先使用改进的麻雀算法(LCSSA)优化BP神经网络的初始权重和阈值,以建立挤压铸造工艺参数的设计模型。然后,以挤压铸造的性能(质量)指标参数为优化目标,利用LCSSA进一步优化了设计模型来获得最优工艺参数。该方法无需过多操作,只需要提供实验数据就可得到对应挤压铸件的最优工艺参数,且可以进行多参数的并行设计,适应不同需求。
Description
技术领域
本发明涉及挤压铸造工艺参数优化设计领域,尤其涉及一种两阶段优化的挤压铸造工艺参数优化设计方法。
背景技术
当前,获取挤压铸造工艺参数主要通过物理实验“试错”方法,大量研究围绕特定材料或特定形状的挤压铸件,选择一定范围或者特定工艺参数进行实验,通过分析实验结果从测试的参数水平中选取最优的工艺参数。这种方式其工艺参数受测试水平或者范围限制,并不一定能得到全局最优解;其次,挤压铸造工艺参数较多,受实验成本限制,常常需要固定一定的实验参数,无法实现所有工艺参数的并行设计。而现有挤压铸造研究只涉及利用优化算法优化神经网络模型预测质量和性能数据,或利用未经优化的网络模型寻找工艺参数组合,预测质量、性能数据和寻找工艺参数组合的效果都不好的问题。因此,如何提高传统神经网络模型的预测能力,利用智能优化算法的两个环节来寻找工艺参数,是实现挤压铸造工艺参数智能优化设计的关键。
发明内容
本发明的目的在于提供一种两阶段优化的挤压铸造工艺参数优化设计方法,解决现有挤压铸造研究只涉及利用优化算法优化神经网络模型预测质量和性能数据,或利用未经优化的网络模型寻找工艺参数组合,预测质量、性能数据和寻找工艺参数组合的效果都不好的技术问题。用于挤压铸造工艺参数智能优化设计的现有优化算法具有收敛速度低、优化时间长以及预测不稳定性的优化神经网络模型等局限性。
为了实现上述目的,本发明采用的技术方案如下:
一种两阶段优化的挤压铸造工艺参数优化设计方法,所述方法包括如下步骤:
步骤1:由于挤压铸造需要铸件在特定温度范围内才能成形,因此相关工艺参数并非在全局范围内搜索,而是有一个初始范围。构建挤压铸造工艺参数优化模型;
步骤2:对于第一阶段,结合改进麻雀算法LCSSA与BP神经网络构建的优化神经网络模型,提高挤压铸造工艺参数神经网络模型的预测性能及鲁棒性;
步骤3:对于第二阶段,在步骤2基础上再次应用改进麻雀算法与信息熵权值法来求解最优挤压铸造工艺参数组合。
进一步地,步骤1中,挤压铸造工艺参数优化模型为:
式中n和σ分别代表挤压铸造工艺参数个数和目标参数个数;Yi代表第i个目标参数,1≤i≤σ;Xj表示第j个挤压铸造工艺参数,1≤j≤t。
进一步地,步骤2中,将神经网络设为输入变量,以铸件性能(或质量)目标参数作为输出变量,以实验数据作为学习样本,通过BP神经网络学习构建工艺参数与目标函数的智能预测模型,隐藏层神经元节点数确定的具体算式为::其中,b1为输入层节点数,b2为输出层节点数,d为隐藏层个数,e为[1,10]之间的常数。
进一步地,步骤2中利用改进麻雀算法对神经网络模型的具体优化过程为:
步骤2.1:为避免改进麻雀算法在后期迭代周期内降低种群多样性,提高初始解的质量,增加种群多样式,结合挤压铸造工艺参数优化模型设置算法初始种群的上下限和优化维度,如下式:
L={[lb1,ub1],…,[lbj′,ubj′]…,[lb|OV|,ub|OV|]}, (3)
式中lbj′和ubj′分别是变量的下限和上限。
而为使算法的震荡性降低,引入Logistic混沌映射增加种群的多样性,其表达式为:
μ=aμ0(1-μ0), (4)
式中μ0∈[0,1],a∈[0,4],当a的值越接近4时,μ的取值范围越能平均分布到[0,1]区域;
将式(3)和式(4)结合起来得到初始种群中个体的数值p*,j′为:
p*,j′=μ(ubj′-lbj′)+lbj′, (5)
步骤2.2:LCSSA优化的主要规则如下:
根据式(5)得到的个体位置生成l个初始麻雀,即t=0。每只麻雀代表一个解决方案。
式中0≤i′≤l;spai′代表群体中的第i只麻雀。
所有麻雀在优化模型中的适应度函数矩阵表示为:
式中每一行的值代表个体的适应度值;
因此,为了对适应度值进行排序,我们有当前最优(fmax(t))和最差(fmin(t))适应度值(或个体)。
式中f(i′,t)表示第i只麻雀在第t次迭代时的当前适应度值。
模拟麻雀算法中发现者、加入者、警戒者的觅食状态,得到最优解。发现者根据式(4-5)的初始种群生成新种群。其位置在每次迭代中按如下方式更新。
式中α∈(0,1]是均匀随机数;R∈[0,1]是预警值;ST∈[0.5,1]是安全阈值;Q是服从正态分布的随机变量。当R<ST时,觅食环境是安全的,发现者进入广域搜索模式以实现更大的适应性,而当R>ST时,觅食环境是危险的,所有麻雀都需要迅速飞到其他安全区域。
麻雀的位置更新可以表示为:
式中是迭代t次维度j′中当前最差的全局位置;/>是第(t+1)次迭代时在维度j′中的全局位置;L是维度为1×l的矩阵,其中所有元素均为1;A是维度为1×l的矩阵,其中每个元素随机为1或-1,并且A+=AT(AAT)。当i′>0.5l时,第i个适应度值较低(即饥饿)的麻雀飞到其他地方以获得足够的能量。
警戒者约占执行预警功能的麻雀的10%-20%。它们的位置更新公式如下:
式中β是方差为1,均值为0正态分布随机数;K∈[-1,1]为随机数;ε是无穷小常数。当f(i′,t)>fmax(t)时,麻雀处于种群边缘,容易受到捕食者的攻击;当f(i′,t)≤fmax(t)发生时,处于种群中间的麻雀意识到危险,需要靠近种群中的其他麻雀,以降低被捕食的概率。
步骤2.3:运用改进麻雀算法LCSSA优化BP神经网络过程如下:
在训练过程中使用LCSSA优化BP神经网络,将BP模型的结构参数(即权重和阈值)设置为LCSSA优化的变量。因此,一组BP模型的这些结构参数对应一只麻雀;一组结构参数的每个值对应麻雀的每个位置。根据式(2)我们有:
|OV|=b1·d+d·b2+d+b2, (12)
由于BP神经网络的权值和阈值都分布在[-1,1],因此|OV|权值阈值的上下限设置为[-1,1]。
同时,设BP神经网络的训练集和测试集的预测值与真实值之差的均方差之和为LCSSA的适应度函数。差值的均方误差和越小,误差越小。
g(k,i′,t)=mse1(k,i′,t)+mse2(k,i′,t)
f(i′,t)=g(1,i′,t)+g(2,i′,t)+…+g(σ,i′,t), (14)
式中r表示有r个训练集,其余为测试集,mse1(k,i′,t)和mse2(k,i′,t)分别为第k个目标参数的第neti′个神经网络在第t次迭代时的均方误差,表示第neti′个神经网络在迭代第t次输出的第k个目标参数。
根据上述设置,通过迭代得到具有最佳适应度值的BP网络模型的最佳权值和阈值,即构建最优神经网络模型。
进一步地,步骤3中挤压铸造工艺参数智能优化具体过程为:
步骤3.1:使用加权集成方式将多个神经网络模型的输出目标整合为单个目标:
f=w1*Y1+…+wσ*Yσ, (15)
式中Yσ代表第σ个目标参数值,wσ代表第σ个权重值。
步骤3.2:利用信息熵权值法来确定挤压铸造工艺参数各目标参数的权值系数。
进一步地,步骤3.2的具体过程为:
设有铸件质量和性能的σ项铸件参数指标,每个铸件目标参数有m个样本,初始铸件目标参数矩阵为:
式中ym,σ代表第m行第σ列目标参数数据;
为消除因各质量指标量纲、数据级不同对结果的影响,对各铸件目标参数进行标准化处理:
式中y′i,k为标准化值,为第k个目标参数的平均值,Rk为第k个目标参数的标准差。
计算第k个目标参数下第i个值的比重zi,k如下:
式中0≤zi,k≤1;使用式(18)建立数据的标准化比重矩阵Z={zi,k}mσ。
第k个目标参数的信息熵权重计算如下:
其中K是一个常数,K=1/lnσ;e为信息熵权值。然后我们有:
dk=l-ek, (20)
第k个目标参数的权重可由下式得到:
通过计算铸件质量和性能指标的权重值(w),可以有效减少因指标重要性的大小对挤压铸造工艺参数智能优化模型构建的影响。
为了使用LCSSA实现最佳工艺参数组合,将每个工艺参数设置为LCSSA的每个变量,并且工艺参数的任意组合在LCSSA中成为麻雀。因此,根据式(3)设置每个挤压铸造工艺参数组合的上限和下限。
然后,式(16)基于要求1中优化的神经网络,用作LCSSA的适应度函数。
基于反向传播(BP)神经网络模型的挤压铸造工艺参数智能优化设计方法与基于混沌映射的改进麻雀搜索算法(LCSSA)相结合的新方法。引入信息熵权法,确定各质量指标对挤压铸造综合质量的影响权重。其次,考虑到传统BP神经网络的低收敛速度、局部最小化以及现有的智能优化算法的局限性,例如优化时间长和容易陷入局部最优,改进麻雀搜索算法用于优化BP神经网络的初始权重和阈值,并建立挤压铸造工艺参数的设计模型。最后,以挤压铸造的性能参数为优化目标,利用改进麻雀搜索算法进一步优化了设计模型。与其他智能优化算法结合BP神经网络建立的优化设计模型相比,通过两阶段优化证明了该方法的有效性。此外,对工艺参数的单目标和多目标优化进行了相关性分析,从而为挤压铸造工艺参数的智能优化设计提供了一种实用的方法。
本发明由于采用了上述技术方案,具有以下有益效果:
(1)本发明基于BP神经网络模型,采用改进的麻雀搜索算法(LCSSA)结合信息熵权法进行两阶段优化过程。在第一阶段,LCSSA算法和BP神经网络结合建立一个优化神经网络模型,预测挤压铸件质量及性能参数指标。在第二阶段,将建立的网络模型与信息熵权法和LCSSA算法相结合,获取了挤压铸造最优工艺参数组合;
(2)设计了两阶段的优化模型,实现对挤压铸造工艺参数的智能优化设计,通过两阶段的优化证明所提方法的有效性,为挤压铸造工艺参数的智能优化设计提供一种切实可行的途径。
附图说明
图1是本发明的LCSSA-BP神经网络流程图。
图2是本发明的挤压铸造工艺参数智能优化设计流程图。
图3是本发明优化算法重复6次完成时间对比。
图4是本发明的不同算法优化的工艺参数组合。
图5是本发明的单目标和多目标优化模型的工艺参数差异。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举出优选实施例,对本发明进一步详细说明。然而,需要说明的是,说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解,即便没有这些特定的细节也可以实现本发明的这些方面。
如图1-2所示,一种两阶段优化的挤压铸造工艺参数优化设计方法,所述方法包括如下步骤:
步骤1:为了提高挤压铸件的性能,减少铸件缺陷和保证生产率,选择铸件质量及性能参数为优化目标,通过构建模型寻找铸件优化目标最小的情况下的工艺参数组合,以此来获得最佳的铸件质量及力学性能。由于挤压铸造需要铸件在特定温度等范围内才能成形,如浇注温度必须在材料的固液相温度以上,因此相关工艺参数并非在全局范围内搜索,而是有一个初始范围。
这样可建立如式(1)中的挤压铸造工艺参数优化模型。
式中n和σ分别代表挤压铸造工艺参数个数和目标参数个数;Yi代表第i个目标参数,1≤i≤σ;Xj表示第j个挤压铸造工艺参数,1≤j≤t。
步骤2:对于第一阶段,结合改进麻雀算法与BP神经网络构建的优化设计的神经网络模型,提高挤压铸造工艺参数神经网络模型的预测性能及鲁棒性;
神经元连接权重、阈值及隐藏层神经元节点的确定是BP神经网络构建的关键问题。隐藏层神经元节点太少难以建立复杂的映射关系、造成误差过大,节点太多,容易造成过拟合。设定隐藏层神经元节点数可由式(2)来决定。
其中,b1为输入层节点数,b2为输出层节点数,d为隐藏层个数,e为[1,10]之间的常数。
传统的麻雀算法采用简单的随机生成方法初始化种群,容易导致最优解周围的个体较少,容易落入局部最优。为避免算法在后续迭代周期内种群多样性减少,提高初始解的质量和增加种群多样性,结合式(1)设置算法初始种群的上下限和优化维度,如式(3)所示。
L={[lb1,ub1],…,[lbj′,ubj′],…,[lb|OV|,ub|OV|]}, (3)
式中lbj′和ubj′分别是变量的下限和上限。
然后,引入Logistic混沌映射增加种群的多样性,以此来使算法的求解趋于稳定。其Logistic混沌映射的数学表达式为:
μ=aμ0(1-μ0), (4)
式中μ0∈[0,1],a∈[0,4],当a的值越接近4时,μ的取值范围越能平均分布到[0,1]区域,当a取为4时,系统处于完全混沌状态,映射分布均匀性达到极值。
将式(3)和式(4)结合起来得到初始种群中个体的数值p*,j′为:
p*,j′=μ(ubj′-lbj′)+lbj′, (5)
步骤2.2:在觅食过程中,通过不断更新三者位置,完成资源的获取和躲避捕食者的攻击。LCSSA优化的主要规则如下:
根据式(5)得到的个体位置生成l个初始麻雀,即t=0。每只麻雀代表一个解决方案。
式中0≤i′≤l;spai′代表群体中的第i只麻雀。
所有麻雀在优化模型中的适应度函数矩阵表示为:
式中每一行的值代表个体的适应度值;
因此,为了对适应度值进行排序,我们有当前最优(fmax(t))和最差(fmin(t))适应度值(或个体)。
式中f(i′,t)表示第i只麻雀在第t次迭代时的当前适应度值。
模拟麻雀算法中发现者、加入者、警戒者的觅食状态,得到最优解。发现者根据式(4-5)的初始种群生成新种群。其位置在每次迭代中按如下方式更新。
式中α∈(0,1]是均匀随机数;R∈[0,1]是预警值;ST∈[0.5,1]是安全阈值;Q是服从正态分布的随机变量。当R<ST时,觅食环境是安全的,发现者进入广域搜索模式以实现更大的适应性,而当R>ST时,觅食环境是危险的,所有麻雀都需要迅速飞到其他安全区域。
麻雀的位置更新可以表示为:
式中是迭代t次维度j′中当前最差的全局位置;/>是第(t+1)次迭代时在维度j′中的全局位置;L是维度为1×l的矩阵,其中所有元素均为1;A是维度为1×l的矩阵,其中每个元素随机为1或-1,并且A+=AT(AAT)。当i′>0.5l时,第i个适应度值较低(即饥饿)的麻雀飞到其他地方以获得足够的能量。
警戒者约占执行预警功能的麻雀的10%-20%。它们的位置更新公式如下:
式中β是方差为1,均值为0正态分布随机数;K∈[-1,1]为随机数;ε是无穷小常数。当f(i′,t)>fmax(t)时,麻雀处于种群边缘,容易受到捕食者的攻击;当f(i′,t)≤fmax(t)发生时,处于种群中间的麻雀意识到危险,需要靠近种群中的其他麻雀,以降低被捕食的概率。
步骤2.3:根据上述规则重复更新群体和适应度值,直到找到最优适应度值。
运用改进麻雀算法优化BP神经网络过程如下:
步骤2.3:运用改进麻雀算法LCSSA优化BP神经网络过程如下:
在训练过程中使用LCSSA优化BP神经网络,将BP模型的结构参数(即权重和阈值)设置为LCSSA优化的变量。因此,一组BP模型的这些结构参数对应一只麻雀;一组结构参数的每个值对应麻雀的每个位置。根据式(2)我们有:
|OV|=b1·d+d·b2+d+b2, (12)
由于BP神经网络的权值和阈值都分布在[-1,1],因此|OV|权值阈值的上下限设置为[-1,1]。
同时,设BP神经网络的训练集和测试集的预测值与真实值之差的均方差之和为LCSSA的适应度函数。差值的均方误差和越小,误差越小。
g(k,i′,t)=mae1(k,i′,t)+mse2(k,i′,t)
f(i′,t)=g(1,i′,t)+g(2,i′,t)+·…g(σ,i′,t), (14)
式中r表示有r个基于DB的训练集,其余为测试集,mse1(k,i′,t)和mse2(k,i′,t)分别为第k个目标参数的第neti′个神经网络在第t次迭代时的均方误差,表示第neti′个神经网络在迭代第t次输出的第k个目标参数。
根据上述设置,通过迭代得到具有最佳适应度值的BP网络模型的最佳权值和阈值,即构建最优神经网络模型。
通过上述步骤完成基于改进麻雀算法优化的BP神经网络模型,以此来预测挤压铸造质量及性能参数目标。
步骤3:对于第二阶段,在步骤2基础上以挤压铸件性能参数为优化目标,结合信息熵权值法与改进麻雀算法对一阶段构建的设计模型进行求解,以此来获得最优挤压铸造工艺参数组合。
步骤3.1:基于第一阶段中的神经网络模型已具备设计特定材料铸件工艺参数的能力和预测特定工艺参数下性能,但确定最优工艺参数的效率较低。将神经网络模型的输出作为适应度值,继续使用改进的麻雀算法对其工艺参数进行智能优化。当神经网络模型的输出(铸件质量和性能目标)超过一个时,需要将多个目标整合为单个目标以适应改进麻雀算法求解。利用信息熵权值法来确定挤压铸造工艺参数各参数指标的权值系数
f=w1*Y1+…+wσ*Yσ, (15)
式中Yσ代表第σ个目标参数值,wσ代表第σ个权重值
步骤3.2:铸件质量及性能参数系数的确定直接影响挤压铸造工艺参数智能优化模型的可靠性和有效性。由于已有工艺参数-铸件质量、性能数据,此处利用信息熵权值法来确定挤压铸造工艺参数各性能(或质量)目标参数的权值系数,步骤如下:
设有铸件质量和性能等σ项铸件目标参数指标,每个铸件目标参数有m个样本,初始铸件目标参数矩阵为:
式中ym,σ代表第m行第σ列目标参数数据;
为消除因各质量指标量纲、数据级不同对结果的影响,对各铸件参数指标进行标准化处理:
式中y′i,k为标准化值,为第k个目标参数的平均值,Rk为第k个目标参数的标准差。
计算第k个目标参数下第i个值的比重zi,k如下:
式中0≤zi,k≤1;使用式(17)建立数据的标准化比重矩阵Z={zi,k}mσ。
第k个目标参数的信息熵权重计算如下:
其中K是一个常数,K=1/lnσ;e为信息熵权值。然后我们有:
dk=1-ek, (19)
第k个目标参数的权重可由下式得到:
通过计算铸件质量和性能指标的权重值(w),可以有效减少因指标重要性的大小对挤压铸造工艺参数智能优化模型构建的影响。
为了使用LCSSA实现最佳工艺参数组合,将每个工艺参数设置为LCSSA的每个变量,并且工艺参数的任意组合在LCSSA中成为麻雀。因此,根据式(3)设置每个挤压铸造工艺参数组合的上限和下限。
然后,式(15)基于要求1中优化的神经网络,用作LCSSA的适应度函数。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种两阶段优化的挤压铸造工艺参数优化设计方法,其特征在于,所述方法包括如下步骤:
步骤1:构建挤压铸造工艺参数优化模型;
步骤2:使用改进麻雀算法与BP神经网络进行构建的优化神经网络模型,提高挤压铸造工艺参数神经网络模型的预测性能及鲁棒性;
步骤3:再次应用改进麻雀算法与信息熵权值法来求解最优挤压铸造工艺参数组合。
2.根据权利要求1所述的一种两阶段优化的挤压铸造工艺参数优化设计方法,其特征在于:步骤1中,挤压铸造工艺参数优化模型为:
式中n和σ分别代表挤压铸造工艺参数个数和目标参数个数;Yi代表第i个铸件性能(或质量)目标参数,1≤i≤σ;Xj表示第j个挤压铸造工艺参数,1≤j≤t。
3.根据权利要求1所述的一种两阶段优化的挤压铸造工艺参数优化设计方法,其特征在于:步骤2中,将神经网络设为输入变量,以铸件性能或者质量目标参数作为输出变量,以实验数据作为学习样本,通过BP神经网络学习构建工艺参数与目标函数的智能预测模型,隐藏层神经元节点数确定的具体算式为:
其中,b1为输入层节点数,b2为输出层节点数,d为隐藏层个数,e为[1,10]之间的常数。
4.根据权利要求1所述的一种两阶段优化的挤压铸造工艺参数优化设计方法,其特征在于:步骤2中利用改进麻雀算法LCSSA对神经网络模型的具体优化过程为:
步骤2.1:为避免改进麻雀算法在后期迭代周期内降低种群多样性,提高初始解的质量,增加种群多样式,结合挤压铸造工艺参数优化模型设置算法初始种群的上下限和优化维度,如下式:
L={[lb1,ub1],…,[lbj′,ubj′],…,[lb|OV|,ub|OV|]},(3)式中lbj′和ubj′分别是变量的下限和上限;
而为使算法的震荡性降低,引入Logistic混沌映射增加种群的多样性,其表达式为:
μ=aμ0(1-μ0),(4)
式中μ0∈[0,1],a∈[0,4],当a的值越接近4时,μ的取值范围越能平均分布到[0,1]区域;
将式(3)和式(4)结合起来得到初始种群中个体的数值p*,j′为:
p*,j′=μ(ubj′-lbj′)+lbj′,(5)
步骤2.2:LCSSA优化的主要规则如下:
根据式(5)得到的个体位置生成l个初始麻雀,即t=0,每只麻雀代表一个解决方案:
式中0≤i′≤l;spai′代表群体中的第i只麻雀;
所有麻雀在优化模型中的适应度函数矩阵表示为:
式中每一行的值代表个体的适应度值;
因此,为了对适应度值进行排序,我们有当前最优适应度值fmax(t)和最差适应度值fmin(t);
式中f(i′,t)表示第i只麻雀在第t次迭代时的当前适应度值;
模拟麻雀算法中发现者、加入者、警戒者的觅食状态,得到最优解,发现者根据式(6)的初始种群生成新种群。其位置在每次迭代中按如下方式更新;
式中α∈(0,1]是均匀随机数;R∈[0,1]是预警值;ST∈[0.5,1]是安全阈值;Q是服从正态分布的随机变量,当R<ST时,觅食环境是安全的,发现者进入广域搜索模式以实现更大的适应性,而当R>ST时,觅食环境是危险的,所有麻雀都需要迅速飞到其他安全区域;
麻雀的位置更新可以表示为:
式中是迭代t次维度j′中当前最差的全局位置;/>是第(t+1)次迭代时在维度j′中的全局位置;L是维度为1×l的矩阵,其中所有元素均为1;A是维度为1×l的矩阵,其中每个元素随机为1或-1,并且A+=AT(AAT),当i′>5.0l时,第i个适应度值较低(即饥饿)的麻雀飞到其他地方以获得足够的能量;
警戒者约占执行预警功能的麻雀的10%-20%。位置更新公式如下:
式中β是方差为1,均值为0正态分布随机数;K∈[-1,1]为随机数;ε是无穷小常数,当f(i′,t)>fmax(t)时,麻雀处于种群边缘,容易受到捕食者的攻击;当f(i′,t)≤fmax(t)发生时,处于种群中间的麻雀意识到危险,需要靠近种群中的其他麻雀,以降低被捕食的概率;
重复更新群体和适应度值,直到找到最优适应度值;
步骤2.3:运用改进麻雀算法LCSSA优化BP神经网络过程如下:
在训练过程中使用LCSSA优化BP神经网络,将BP模型的结构参数设置为LCSSA优化的变量,结构参数即为权重和阈值,一组BP模型的这些结构参数对应一只麻雀,一组结构参数的每个值对应麻雀的每个位置,根据式(2)得到:
|OV|=b1·d+d·b2+d+b2, (12)
由于BP神经网络的权值和阈值都分布在[-1,1],因此|OV|权值阈值的上下限设置为[-1,1]。
同时,设BP神经网络的训练集和测试集的预测值与真实值之差的均方差之和为LCSSA的适应度函数,差值的均方误差和越小,误差越小;
g(k,i′,t)=mse1(k,i′,t)+mse2(k,i′,t)
f(i′,t)=g(1,i′,t)+g(2,i′,t)+…+g(σ,i′,t),(14)
式中r表示有r个基于DB的训练集,其余为测试集,mse1(k,i′,t)和mse2(k,i′,t)分别为第k个目标参数的第neti′个神经网络在第t次迭代时的均方误差,表示第neti′个神经网络在迭代第t次输出的第k个目标参数;
根据上述设置,通过迭代得到具有最佳适应度值的BP网络模型的最佳权值和阈值,即构建最优神经网络模型。
5.根据权利要求1所述的一种两阶段优化的挤压铸造工艺参数优化设计方法,其特征在于:步骤3中挤压铸造工艺参数智能优化具体过程为:
步骤3.1:使用加权集成方式将多个神经网络模型的输出目标整合为单个目标:
f=w1*Y1+…+wσ*Yσ, (15)
式中Yσ代表第σ个目标参数值,wσ代表第σ个权重值
步骤3.2:利用信息熵权值法来确定挤压铸造工艺参数各目标参数的权值系数。
6.根据权利要求5所述的一种两阶段优化的挤压铸造工艺参数优化设计方法,其特征在于:步骤3.2的具体过程为:
设有铸件质量和性能的σ项铸件参数指标,每个铸件目标参数有m个样本,初始铸件目标参数矩阵为:
式中ym,σ代表第m行第σ列目标参数数据;
为消除各质量指标不同维度和数据水平对结果的影响,应将各目标参数标准化如下:
式中y′i,k为标准化值,为第k个目标参数的平均值,Rk为第k个目标参数的标准差;
计算第k个目标参数下第i个值的比重zi,k如下:
式中0≤zi,k≤1;使用式(18)建立数据的标准化比重矩阵Z={zi,k}mσ。
第k个目标参数的信息熵权重计算如下:
其中K是一个常数,K=1/lnσ;e为信息熵权值。然后我们有:
dk=1-ek,(20)
第k个目标参数的权重可由下式得到:
通过计算铸件质量和性能目标参数的权重值(w),可以有效减少因目标参数重要性的大小对挤压铸造工艺参数智能优化模型构建的影响;
为了使用LCSSA实现最佳工艺参数组合,将每个工艺参数设置为LCSSA的每个变量,并且工艺参数的任意组合在LCSSA中成为麻雀;因此,根据式(3)设置每个挤压铸造工艺参数组合的上限和下限;
然后,式(15)基于要求1中优化的神经网络,用作LCSSA的适应度函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310180272.3A CN116522747A (zh) | 2023-02-28 | 2023-02-28 | 一种两阶段优化的挤压铸造工艺参数优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310180272.3A CN116522747A (zh) | 2023-02-28 | 2023-02-28 | 一种两阶段优化的挤压铸造工艺参数优化设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116522747A true CN116522747A (zh) | 2023-08-01 |
Family
ID=87398314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310180272.3A Pending CN116522747A (zh) | 2023-02-28 | 2023-02-28 | 一种两阶段优化的挤压铸造工艺参数优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116522747A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117371341A (zh) * | 2023-12-07 | 2024-01-09 | 小米汽车科技有限公司 | 工艺参数确定方法、目标件制造方法、装置及介质 |
CN117574721A (zh) * | 2023-11-20 | 2024-02-20 | 华中科技大学 | 一种工艺参数概率模型优化方法 |
-
2023
- 2023-02-28 CN CN202310180272.3A patent/CN116522747A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117574721A (zh) * | 2023-11-20 | 2024-02-20 | 华中科技大学 | 一种工艺参数概率模型优化方法 |
CN117574721B (zh) * | 2023-11-20 | 2024-09-06 | 华中科技大学 | 一种工艺参数概率模型优化方法 |
CN117371341A (zh) * | 2023-12-07 | 2024-01-09 | 小米汽车科技有限公司 | 工艺参数确定方法、目标件制造方法、装置及介质 |
CN117371341B (zh) * | 2023-12-07 | 2024-04-09 | 小米汽车科技有限公司 | 工艺参数确定方法、目标件制造方法、装置及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111262858B (zh) | 基于sa_soa_bp神经网络的网络安全态势预测方法 | |
CN116522747A (zh) | 一种两阶段优化的挤压铸造工艺参数优化设计方法 | |
CN111680786B (zh) | 一种基于改进权重门控单元的时序预测方法 | |
CN112734002B (zh) | 一种基于数据层和模型层联合迁移学习的寿命预测方法 | |
CN113722980B (zh) | 海洋浪高预测方法、系统、计算机设备、存储介质、终端 | |
CN112215446A (zh) | 一种基于神经网络的单位动态火灾风险评估方法 | |
CN111861013A (zh) | 一种电力负荷预测方法及装置 | |
CN113947182B (zh) | 基于双阶段堆叠图卷积网络的交通流预测模型构建方法 | |
CN116542382A (zh) | 基于混合优化算法的污水处理溶解氧浓度预测方法 | |
CN112733997A (zh) | 基于woa-lstm-mc的水文时间序列预测优化方法 | |
CN112766603A (zh) | 一种交通流量预测方法、系统、计算机设备及存储介质 | |
CN115018193A (zh) | 基于lstm-ga模型的时间序列风能数据预测方法 | |
CN112765902A (zh) | 基于TentFWA-GD的RBF神经网络软测量建模方法及其应用 | |
CN118095570A (zh) | 台区智能负荷预测方法、系统、电子设备、介质及芯片 | |
CN117910823A (zh) | 基于注意力机制的优化lstm短期负荷预测方法 | |
CN113887570B (zh) | 一种基于神经网络的太阳耀斑二分类预测方法 | |
CN115996135A (zh) | 一种基于特征组合优化的工业互联网恶意行为实时检测方法 | |
CN113884936A (zh) | 一种基于issa耦合delm的锂离子电池健康状态预测方法 | |
Huang et al. | Wine quality evaluation model based on artificial bee colony and BP neural network | |
CN118364342B (zh) | 一种火控计算机电源模块故障诊断方法 | |
CN111709140B (zh) | 一种基于内在可塑性回声状态网络的船舶运动预报方法 | |
CN111753470B (zh) | 一种基于粗糙集和深度置信神经网络的建筑能耗预测方法 | |
Hou et al. | Deep Learning-Based Comparative Modeling of Carbon Emissions Projections | |
CN117828305A (zh) | 一种矿井涌水量预测方法、系统、设备及介质 | |
CN118098439A (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 |