CN113408059B - 一种汽车复合材料保险杠防撞梁厚度的多目标优化方法 - Google Patents

一种汽车复合材料保险杠防撞梁厚度的多目标优化方法 Download PDF

Info

Publication number
CN113408059B
CN113408059B CN202110737146.4A CN202110737146A CN113408059B CN 113408059 B CN113408059 B CN 113408059B CN 202110737146 A CN202110737146 A CN 202110737146A CN 113408059 B CN113408059 B CN 113408059B
Authority
CN
China
Prior art keywords
thickness
layer
plate
function
collision beam
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
Application number
CN202110737146.4A
Other languages
English (en)
Other versions
CN113408059A (zh
Inventor
陈静
徐森
崔晓凡
胡金旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN202110737146.4A priority Critical patent/CN113408059B/zh
Publication of CN113408059A publication Critical patent/CN113408059A/zh
Application granted granted Critical
Publication of CN113408059B publication Critical patent/CN113408059B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/12—Computing arrangements based on biological models using genetic models
    • G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • 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
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Geometry (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Software Systems (AREA)
  • Biophysics (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Databases & Information Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Operations Research (AREA)
  • Physiology (AREA)
  • Genetics & Genomics (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种汽车复合材料保险杠防撞梁厚度的多目标优化方法，包括：步骤一、分别确定防撞梁的上板、下板、前板、后板和肋板的铺层厚度取值范围，以及各铺层厚度对应的铺层顺序；步骤二、生成初始样本数据集，并确定防撞梁铺层厚度优化的目标函数f(x)和约束函数g(x)；步骤三、根据所述目标函数和所述约束函数，设置高斯随机过程核函数初始参数，并训练得到GPR模型；步骤四、通过GPR模型得到Pareto前沿解集，并在所述Pareto前沿解集中筛选出新的样本，加入所述初始样本数据集中；循环步骤二至步骤四，直到满足迭代次数，得到优化的Pareto前沿解集，并根据所述优化的Pareto前沿解集得到防撞梁的上板、下板、前板、后板和肋板的优化铺层厚度。

Description

一种汽车复合材料保险杠防撞梁厚度的多目标优化方法
技术领域
本发明属于汽车被动安全性研究技术领域，特别涉及一种汽车复合材料保险杠防撞梁厚度的多目标优化方法。
背景技术
目前，多目标优化大都是依据代理模型的优化方法，优化前需建立具备足够精度的代理模型，通常情况下，代理模型精度只有达到0.9以上，其优化结果才准确可信。但是，对于高维高度非线性的数学模型，比如汽车碰撞等，需要采用几百组甚至上千组的试验设计(DOE)，进行仿真与模拟。这个代价是不可承受的。除此之外，大多数传统代理模型对于模型的不确定性并没有直接反映，导致优化的结果产生较大的偏差，当面对更高维度的优化时，优化无法顺利进行。还有，对于一步优化设计，建模具有一次性，没有近似模型提升的过程，若模型存在较大近似误差，基于该模型的优化必然不准确。目前，尚未有结合序列优化与高斯过程回归模型实现多目标优化的优化方法。
发明内容
本发明设计开发了一种汽车复合材料保险杠防撞梁厚度的多目标优化方法，采用高斯过程回归的第三代非支配排序序列优化方法，只需获取少量抽样点，在逐步优化中更新样本点，加入样本点进行模型重构并进行多目标优化即可获得较为准确地Pareto前沿；从而快速准确的得到满足碰撞性能要求的碳纤维复合材料防撞梁厚度优化结果。
本发明提供的技术方案为：
一种汽车复合材料保险杠防撞梁厚度的多目标优化方法，包括如下步骤：
步骤一、分别测量多个防撞梁的上板、下板、前板、后板和肋板的铺层厚度，得到最大值和最小值从而确定取值范围(以测量结果的最大值作为取值范围上限，最小值作为取值范围下限)，以及确定各铺层厚度对应的铺层顺序；
步骤二、生成初始样本数据集，并确定防撞梁铺层厚度优化的目标函数f(x)和约束函数g(x)；
其中，所述初始样本数据集中的样本包括上板、下板、前板、后板和肋板的铺层厚度以及所述铺层厚度对应的防撞梁的质量、防撞梁的碰撞力、防撞梁的吸能量和防撞梁的入侵量；
f(x)＝(m,F)，g(x)＝(E,D)；
式中，m表示最小化防撞梁的质量，F表示最小化防撞梁碰撞力，E表示入侵量为60mm以下，D表示吸能量为400J以上；
步骤三、根据所述目标函数和所述约束函数，设置高斯随机过程核函数初始参数，并训练得到GPR模型；
步骤四、通过GPR模型得到Pareto前沿解集，并在所述Pareto前沿解集中筛选出新的样本，加入所述初始样本数据集中；
循环步骤二至步骤四，直到满足迭代次数，得到优化的Pareto前沿解集，并根据所述优化的Pareto前沿解集得到防撞梁的上板、下板、前板、后板和肋板的优化铺层厚度。
优选的是，在所述步骤一中，确定防撞梁的上板、下板、前板、后板和肋板的铺层厚度取值范围为：
Figure BDA0003140337470000021
式中，x1表示防撞梁的后板的铺层厚度，x2表示防撞梁的前板的铺层厚度，x3表示防撞梁的肋板的铺层厚度，x4表示防撞梁的下板的铺层厚度，x5表示防撞梁的上板的铺层厚度。
优选的是，其特征在于，x4＝x5。
优选的是，在所述步骤一中，各铺层厚度对应的铺层顺序为：
铺层厚度为3mm时，铺层顺序为：[(0/90/±45)2(±45)2]s；
铺层厚度为2.5mm时，铺层顺序为：[(0/90/±45)2(±45)]s；
铺层厚度为2mm时，铺层顺序为：[(0/90/±45)]2s；
铺层厚度为1.5mm时，铺层顺序为：[0/90/±45/±45]s；
铺层厚度为1mm时，铺层顺序为：[0/90/±45]s。
优选的是，在所述步骤四中，在所述Pareto前沿解集中筛选出新的样本，包括如下步骤：
步骤1、计算Pareto前沿解集中每个个体基于欧拉距离的多目标期望提高准则EIMe，其中：
Figure BDA0003140337470000031
Figure BDA0003140337470000032
Figure BDA0003140337470000033
Figure BDA0003140337470000034
式中，S为Pareto前沿解集，下标表示目标个数，上标表示解的个数；EIMe(x)为基于欧拉距离的期望提高准则，EI表示期望提高矩阵，EI每一行表示对于处于第j个前沿位置所有目标函数的提高值，EI每一列表示对于所有前沿位置第j个目标函数的提高值；fi j(x)为Pareto前沿解集中第j号个体的第i个目标函数值，
Figure BDA0003140337470000035
s-fi(x)为第i个目标函数的GPR模型对应的预测值以及误差函数值；φ(.),Φ(.)为标准正态分布的概率密度函数和累积分布函数；
步骤2、选择EIMe值最大的个体作为新的样本。
优选的是，在所述步骤二中，通过最优拉丁超方试验生成初始样本。
优选的是，通过最优拉丁超方试验生成初始样本数量为50个。
本发明的有益效果是：
本发明提供的汽车复合材料保险杠防撞梁厚度的多目标优化方法，采用高斯过程回归的第三代非支配排序序列优化方法，只需获取少量抽样点，在逐步优化中更新样本点，加入样本点进行模型重构并进行多目标优化即可获得较为准确地Pareto前沿；序列优化过程充分考虑了模型的误差信息，在多次建模的过程中实现优化结果的全局性提高；同时也缩短了仿真的时间和成本，且具备较好的准确率；从而快速准确的得到满足碰撞性能要求的碳纤维复合材料防撞梁厚度优化结果。
附图说明
图1为本发明所述的汽车复合材料保险杠防撞梁的结构示意图。
图2为本发明中计算解空间中参考线与归一化目标之间的距离的示意图。
图3为本发明所述的优化Pareto前沿面示意图。
具体实施方式
下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
本发明提供了一种汽车复合材料保险杠防撞梁厚度的多目标优化方法，具体过程如下：
如图1所示，设定设计变量x1代表防撞梁后板110的厚度，x2代表防撞梁前板120的厚度，x3代表防撞梁肋板130的厚度，x4代表防撞梁下板140的厚度，x5代表防撞梁上板150的厚度，定义上下板厚度相同；设计变量采用离散取值的方式，取值的间距为0.5。
因此铺层厚度设计变量的取值范围如下：
Figure BDA0003140337470000041
各厚度下铺层顺序如表1所示：
表1各铺层厚度对应的铺层顺序
铺层厚度 铺层顺序
3mm [(0/90/±45)<sub>2</sub>(±45)<sub>2</sub>]<sub>s</sub>
2.5mm [(0/90/±45)<sub>2</sub>(±45)]<sub>s</sub>
2mm [(0/90/±45)]<sub>2s</sub>
1.5mm [0/90/±45/±45]<sub>s</sub>
1mm [0/90/±45]<sub>s</sub>
定义优化数学模型为：
Minimize fm(x),fF(x)
Subjectto gE(x)≥400J
gD≤60.0mm
1.0≤x1≤2.5
1.0≤x2≤3
1.0≤x3≤2.5
1.0≤x4≤2。
数学模型构建完成后，依据目标函数以及约束函数进行多目标优化设计，采用高斯过程回归的第三代非支配排序遗传序列优化，主要包括如下步骤：
步骤1：在初始设计空间中对复杂系统进行最优拉丁超方试验(OLHS)生成初始样本点，并获得初始样本点处的目标函数f(x)和约束函数g(x)，如表2所示；
其中，f(x)＝[f1(x),f2(x),...fj(x)...,fn(x)]；
g(x)＝[g1(x),g2(x),...gj(x)...,gt(x)]；
式中，n、t分别为目标函数、约束函数的个数，x为设计变量。
表2初始样本表
Figure BDA0003140337470000051
Figure BDA0003140337470000061
Figure BDA0003140337470000071
依据对铺层结构划分，防撞梁设计变量共有5个，因为上下层合板采用相同的铺层方式，可以将设计变量缩减至4个，变量设计范围已经给出，根据层合板铺设的单层厚度，进行50组最优拉丁超立方抽样设计。
在有限元求解后获得目标函数f(x)和约束函数g(x)：
Figure BDA0003140337470000072
步骤2：根据目标函数和约束函数，设置高斯随机过程核函数初始参数；模型训练方法为极大似然估计，依据贝叶斯推断超参数θ的后验密度函数为：
Figure BDA0003140337470000073
P(Y|X)与参数θ无关，所以该分布的对于一个给定观测值的关于参数的后验分布来说相当于是一个常数，并不会影响到后验分布的均值与方差常数，此时超参数的最大后验估计
Figure BDA00031403374700000811
最大化P(Y|X,θ)P(θ),当先验P(θ)为均匀分布时，超参数的最大后验估计就是寻找使得P(Y|X,θ)最大的θ，即极大似然估计。对边际似然取对数作为似然函数((最大化)对数似然边际LML)，对超参数求偏导：
Figure BDA0003140337470000081
Figure BDA0003140337470000082
依据初始样本进行回归模型超参数调优以及模型训练，其中一个关键的设置为高维映射核函数的选择，单一核方法对于混杂着各种数值噪声的数据往往无法得到较为准确的回归模型。故设置初始化先验均值设置为0；meanfunc＝[]；核函数为复合核，具体包括：
平方指数(SE)函数：
Figure BDA0003140337470000083
马特温协方差(Matern)函数：
Figure BDA0003140337470000084
其中，
Figure BDA0003140337470000085
为平滑系数，l为长度尺度参数，Γ(x)为欧拉第二积分，在实数域上的定义为：
Figure BDA0003140337470000086
Figure BDA00031403374700000810
为第二类贝塞尔函数，是下列微分方程的标准解函数y(x):
Figure BDA0003140337470000087
有理二次(RQ)函数：
Figure BDA0003140337470000088
(l为长度尺度参数，α为比例混合系数)；
方差核函数选择有理二次和平方指数核函数，设置为
Figure BDA0003140337470000089
σ初始值为2，l1初始值为0.25，l2初始值为2，α为10-3。
并采用共轭梯度法对超参数进行优化；
Figure BDA0003140337470000091
式中，gp是回归函数，hyp，Hyper为核函数超参数，包括均值函数meanfunc，协方差函数covfunc。Minimize为共轭梯度优化函数，Sample_x为原设计点，此处即为设计变量[x1,x2,x3,x4]。Sample_f为样本点的目标函数，即数据集中的[m,F],要求最小化防撞梁的质量，最小化防撞梁碰撞力。Sample_g为约束函数；为数据集中的[D,E]；约束侵入量不得大于60mm，吸能量不得小于400J。
步骤3：生成参考点信息并初始化个体信息；进行模拟二进制交叉以及多项式变异生成新个体。
首先生成初始参考点，参考点生成方式如下：其中，M为目标函数的维度，H为分段数；
1)产生
Figure BDA0003140337470000092
个维数为M-1的初始参考点集；此处维数为2。
X＝(1/H)*nchoosek(0:H+M-2,M-1)；
其中，nchoosek为抽样函数，第一个参数表示抽样总体，第二个参数表示抽样的个体数；
2)对于
Figure BDA0003140337470000093
进行如下换:aij＝xij-(j-1)/H,获得变换集A；
3)利用变换集A,参考点依照下式生成:
Figure BDA0003140337470000094
通常，因划分段数H以及维数M增大会使参考点集过大，故采用一种内外层机制：H1为边界分段数，H2为内层分段数；
H1,H2大小应满足
Figure BDA0003140337470000095
N为初始生成个体总数，且内层数应小于等于边界分段数；边界层遵照上述3个过程生成参考点S1，内层首先生成参考点S(3个过程)，取映射系数为0.5,S2由下列公式产生,并合并S1,S2；
Figure BDA0003140337470000101
对于2目标而言，H1＝149，H2＝0。设置种群个体数初始化为150，设置迭代次数200，
用GPR模型获取目标函数和约束函数的预测值以及对应的方差，并构建适应度信息；
基于核方法将贝叶斯线性回归映射成高维空间的非线性回归，对于函数f(x)视为下式所示服从高斯分布的高斯随机过程：
P(f*|X,y,X*,σ2)～N(σ-2φ(X*)TA-1φy,φ(X*)TA-1φ(X*))
Figure BDA0003140337470000102
Figure BDA0003140337470000103
μ*＝k(X*,X)(K+σ2I)-1y
σ*＝k(X*,X*)-k(X*,X)(K+σ2I)-1k(X,X*)
获取2中超参数，结合分布信息估计未知个体位置处的目标函数、约束函数信息以及模型的误差，并利用惩罚机制构建用于非支配比较的适应度函数；
Figure BDA0003140337470000104
式中，a11,a12,…,an1,…ant为惩罚系数，PopDec是初始化个体位置信息，sf(x),sg(x)为目标函数和约束函数的预测值的误差函数，F为适应度函数；依照上述步骤进行种群个体目标函数，约束函数预测值和方差估计：PopDec是经过遗传操作产生的个体位置信息，维度与变量相同。
[f1,sf1]＝gp(hyp1,'infGaussLik',[],covfunc2,@likGauss,Sample_x,Sample_f1,PopDec)；
[f2,sf2]＝gp(hyp2,'infGaussLik',[],covfunc2,@likGauss,Sample_x,Sample_f2,PopDec)；
[g1,sg1]＝gp(hyp3,'infGaussLik',[],covfunc2,@likGauss,Sample_x,Sample_g1,PopDec)；
[g2,sg2]＝gp(hyp4,'infGaussLik',[],covfunc2,@likGauss,Sample_x,Sample_g2,PopDec)；
f12＝[f1 f2]；
f12_interval＝[sf1 sf2]；
g12＝[-g1+400g2-60]；
g12_interval＝[sg1 sg2]；
z＝[f12 g12]；
f_penalty＝100*max(0,-g1+400)+100*max(0,g2-60)；
y＝f12+f_penalty；
y为构造的适应度函数，Samplex是原始训练数据中的设计变量，Sample_f1、Sample_f2、Sample_g1、Sample_g2对应质量、碰撞力、吸能量、侵入量响应。f1为防撞梁质量模型预测值，g1为防撞梁吸能模型预测值，f2为防撞梁碰撞力模型预测值，g2为防撞梁吸能模型预测值。Sf,Sg分别是目标函数和约束函数预测值方差。
根据遗传操作进行个体位置的产生：
ILower＝repmat(Lower,N,1)
IUpper＝repmat(Upper,N,1)
Indiv_Pos＝unifrnd(ILower,IUpper)
其中，Upper、Lower为变量上、下限，
lower＝[1 1 1 1]；upper＝[2.5 3 2.5 2]；
repmat为重复构造函数，N＝150；unifrnd为随机数函数，产生一个NxN的介于Lower与Upper的随机数矩阵，保证样本的随机性、多样性；之后，按照遗传操作流程，首先确定交配池大小，并根据交配池进行筛选，其具体实现方法为：
随机产生行数为2，列数为N的个体序号(可重复，序号介于1～N),根据粒子违反约束的函数值，该值总是小于0的，且值越小对于约束的满足度越高，先随机产生两组个体，依据违约值对两个体取最优：
Tool＝randi(N,M,N)；
对违约值进行升序排序，得到对应原来位置的索引：
[～,index1]＝Sort(Popcon，[],1)；
本例中，Popcon＝0；再对索引进行升序排序，得到index2，此时index2便包含原位置个体对应的大小排序序号:
[～,index2]＝sortrows(index1)
根据index2信息，可知交配池中个体的大小排序，因为序列号越小，满足约束的可能性越大，所以对比两个个体，将对应较小个体排序序号选出；对于子代个体通过模拟二进制交叉以及多项式变异进行；模拟二进制交叉，模拟遗传算法二进制编码中的单点交叉：
主要遵循交叉前后父代与子代的均值相等：
Figure BDA0003140337470000121
定义传播因子(Spread Factorβ)，可得c1与c2:
Figure BDA0003140337470000122
Figure BDA0003140337470000123
Figure BDA0003140337470000124
其中，传播因子β可由下式确定：
Figure BDA0003140337470000125
其中，o为分布系数；同样，对于多项式变异，遵循以下公式：
indivk'＝indivk+δ(Upperk-Lowerk)
Figure BDA0003140337470000131
δ1＝(indivk-Lowerk)/(Upperk-Lowerk)
δ1＝(indivk-Upperk)/(Upperk-Lowerk)
本例中，r为一个随机算子，o,η＝20，若随机算子小于等于0.5，则按照上面式子计算传播因子，变异系数。否则按照第二个式子计算传播因子，变异系数。
步骤4：对步骤3生成的个体(总数为2*N)，依照适应度评估，划分非支配等级，定义储备集St(N(St)≥N),定义Fk为等级为k的非支配解集；对储备解集St进行归一化，在非支配等级划分时采用一种高效非支配排序方法ENS(EffcientNon_dominatedSort)减少计算中时间和空间复杂度，定义MaxFNo为当前最高支配等级(数字越小支配等级越高),初始值设定为0；FrontNo为前沿集合，初始设置为INF，涉及的主要步骤为：
1)依照对所有个体的适应度第一维度进行降序排序，获得排序信息；
2)定义大循环结束条件为FrontNo中不为INF的总个体数超过N，每次循环开始对MaxFNo进行加一操作，记录满足截止条件时最大非支配等级；对于当前个体，将其与之前处于前沿解集的个体进行非支配比较，比较从第二维度开始进行，若个体被支配，即进行下个个体的非支配比较；若个体未被支配，个体位置对应的FrontNo为当前MaxFNo；
3)上述步骤2结束，根据排序信息以及MaxFNo，获取非支配等级小于MaxFNo的所有个体，定义为Pt+1,非支配等级等于MaxFNo的个体定义为Fm,St为Pt+1和Fm的并集；
对获取到的St进行如下所示的归一化操作：
4)理想点位置获取，对于最小值问题，Zmin＝min(f1,f2,…fn)；
5)坐标系转换，以理想点作为原点，变化后函数目标f’＝f-Zmin对于多目标问题，需要获取对应各个目标轴的极值点以及在各个目标轴的截距，其获取机制为:
Figure BDA0003140337470000141
XA＝1
dir为方向向量，大小为nxn,主对角线为1，其余位置均为1e-6；X为n个极值点的集合，大小为nxn，A为目标截距集合，大小为nx1，1是全为1且大小为nx1的列向量。
如果矩阵X的秩小于n，那么这n个极限点就不能构成一个n维的超平面。甚至即使超平面能够建立，也可能在某些方向上得不到截距或某些截距ai不满足
Figure BDA0003140337470000142
在所有上述情形下，对每个i∈{1,2,...,n}，
Figure BDA0003140337470000143
设置为St中的非支配解在目标fi上的最大值。对于St中各个目标函数进行归一化处理：
Figure BDA0003140337470000144
步骤5：如图2所示，以三维为例，计算解空间中参考线(参考点与理想点之间的连线)与归一化目标之间的距离，依照下式生成距离矩阵D：
Figure BDA0003140337470000145
Figure BDA0003140337470000146
其中，q代表St中个体数量，r对应参考点的数量；St每个个体考虑距离每个参考点距离，选择具有最小距离的参考点进行依附；
步骤6：考虑前MaxFNo-1前沿等级的情况下，根据参考点被选择的个数，选择参考点被选择较少的作为最后一前沿种群个体选择的依据，进而实现选择。个体选择依照小生境保持算子，步骤5取得的St所有个体对应参考点信息，定义ρj为Pt+1依附在参考点j的个体总数目，进行如下判断：
1)参考点选择，Jmin＝{j:argminρj}，若Jmin存在多个，随机选择
Figure BDA0003140337470000147
2)进行Fm(最后前沿解集)的选择，存在两种情况：
首先，若Fm中并未发现存在个体与
Figure BDA0003140337470000148
相对应，那么将
Figure BDA0003140337470000149
在当前迭代步不在考虑，从剩余参考点位置按照步骤1筛选新的参考点；若Fm存在个体j与关联，同样存在两种情况：第一种情况为，对应最少数目个体的参考点值为0，表示对于之前非支配前沿中并未存在个体与之对应，为了保持个体的多样性，从Fm中(如果存在多个，依照步骤5计算各个个体距离该参考点对应参考线的距离，取距离最小的个体)选取个体加入Pt+1，ρj的数值随着加入个体而增大；第二种情况，对应最少数目个体的参考点值不为0，表示已经存在个体与参考点与之对应，从Fm中依附于该参考点的众多个体(如果存在多个)进行随机挑选加入；ρj的数值随着加入个体而增大；
上述过程重复执行，直至Pt+1中的个体数达到需求的个体数N；依据设置迭代次数，个体不断演变进化，得到近似Pareto近似前沿解；
步骤7：将所述近似Pareto前沿解集作为真实函数的真实Pareto前沿解集，并获取每个个体基于欧拉距离的多目标期望提高准则：
当前最优解是由多个多维点共同组成的近似Pareto前沿：
Figure BDA0003140337470000151
Figure BDA0003140337470000152
Figure BDA0003140337470000153
Figure BDA0003140337470000154
式中，S为近似Pareto解的情况，下标表示目标个数，上标表示近似解的个数；EIMe(x)为基于欧拉距离的期望提高准则，EI表示期望提高矩阵，EI每一行表示对于处于第j个前沿位置所有目标函数的提高值，EI每一列表示对于所有前沿位置第j个目标函数的提高值；fi j(x)为近似Pareto前沿解集中第j号个体的第i个目标函数值，
Figure BDA0003140337470000155
s-fi(x)为第i个目标函数的高斯过程回归模型对应的预测值以及误差函数值；φ(.),Φ(.)为标准正态分布的概率密度函数和累积分布函数；
选择EIMe值最大的个体作为新的样本点加入原样本点中生成新的样本点集；将所述新的样本点集作为初始样本点，重复步骤1-7，直至满足迭代次数，输出最终的样本点集。
按照上述步骤1～7，设置优化结果如图3所示，图3为优化Pareto前沿面，优化后依据基于熵值法的Toposis准则(Ci越接近1，个体越好)选取前沿Pareto中一个目标解：
熵值法(确定权值)：若前沿pareto解集中含有m个个体，每个个体为n维度目标函数。如下式：
Figure BDA0003140337470000161
对上式进行按列归一化处理:
Figure BDA0003140337470000162
计算n个目标函数各自的熵权值：
Figure BDA0003140337470000163
Toposis方法：对PY每个元素进行判断，若属于极大型指标(越大越好)，不做处理，若属于极小型指标，进行指标正向化：
Figure BDA0003140337470000164
得到同向化处理矩阵:
Figure BDA0003140337470000165
最优方案为同向化处理后每一列的最大值
Figure BDA0003140337470000166
最劣方案为每一列的最小值
Figure BDA0003140337470000167
计算每个个体距离最优方案以及最劣方案的接近程度：
Figure BDA0003140337470000171
Figure BDA0003140337470000172
Figure BDA0003140337470000173
[x1,x2,x3,x4,x5]＝[1,1.5262,1.4610,1.1295,1.1295]
[x1,x2,x3,x4,x5]＝[1.0,1.5,1.5,1.0,1.0]
[E,F,D,M]＝[435.5802,10.6882,59.9988,0.8668]
优化结果满足目标以及约束条件。
尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。

Claims (6)

1.一种汽车复合材料保险杠防撞梁厚度的多目标优化方法，其特征在于，包括如下步骤：
步骤一、分别测量多个防撞梁的上板、下板、前板、后板和肋板的铺层厚度，得到最大值和最小值从而确定取值范围，以及确定各铺层厚度对应的铺层顺序；
步骤二、生成初始样本数据集，并确定防撞梁铺层厚度优化的目标函数f(x)和约束函数g(x)；
其中，所述初始样本数据集中的样本包括上板、下板、前板、后板和肋板的铺层厚度以及所述铺层厚度对应的防撞梁的质量、防撞梁的碰撞力、防撞梁的吸能量和防撞梁的入侵量；
f(x)＝(m,F)，g(x)＝(E,D)；
式中，m表示最小化防撞梁的质量，F表示最小化防撞梁碰撞力，E表示入侵量为60mm以下，D表示吸能量为400J以上；
步骤三、根据所述目标函数和所述约束函数，设置高斯随机过程核函数初始参数，并训练得到GPR模型；
步骤四、通过GPR模型得到Pareto前沿解集，并在所述Pareto前沿解集中筛选出新的样本，加入所述初始样本数据集中；
在所述Pareto前沿解集中筛选出新的样本，包括如下步骤：
步骤1、计算Pareto前沿解集中每个个体基于欧拉距离的多目标期望提高准则EIMe，其中：
Figure FDA0003591384610000011
Figure FDA0003591384610000012
Figure FDA0003591384610000021
Figure FDA0003591384610000022
式中，S为Pareto前沿解集，下标表示目标个数，上标表示解的个数；EIMe(x)为基于欧拉距离的期望提高准则，EI表示期望提高矩阵，EI每一行表示对于处于第j个前沿位置所有目标函数的提高值，EI每一列表示对于所有前沿位置第j个目标函数的提高值；fi j(x)为Pareto前沿解集中第j号个体的第i个目标函数值，
Figure FDA0003591384610000023
s_fi(x)为第i个目标函数的GPR模型对应的预测值以及误差函数值；φ(.),Φ(.)为标准正态分布的概率密度函数和累积分布函数；
步骤2、选择EIMe值最大的个体作为新的样本；
循环步骤二至步骤四，直到满足迭代次数，得到优化的Pareto前沿解集，并根据所述优化的Pareto前沿解集得到防撞梁的上板、下板、前板、后板和肋板的优化铺层厚度。
2.根据权利要求1所述的汽车复合材料保险杠防撞梁厚度的多目标优化方法，其特征在于，在所述步骤一中，确定防撞梁的上板、下板、前板、后板和肋板的铺层厚度取值范围为：
Figure FDA0003591384610000024
式中，x1表示防撞梁的后板的铺层厚度，x2表示防撞梁的前板的铺层厚度，x3表示防撞梁的肋板的铺层厚度，x4表示防撞梁的下板的铺层厚度，x5表示防撞梁的上板的铺层厚度。
3.根据权利要求2所述的汽车复合材料保险杠防撞梁厚度的多目标优化方法，其特征在于，x4＝x5。
4.根据权利要求3所述的汽车复合材料保险杠防撞梁厚度的多目标优化方法，其特征在于，在所述步骤一中，各铺层厚度对应的铺层顺序为：
铺层厚度为3mm时，铺层顺序为：[(0/90/±45)2(±45)2]s；
铺层厚度为2.5mm时，铺层顺序为：[(0/90/±45)2(±45)]s；
铺层厚度为2mm时，铺层顺序为：[(0/90/±45)]2s；
铺层厚度为1.5mm时，铺层顺序为：[0/90/±45/±45]s；
铺层厚度为1mm时，铺层顺序为：[0/90/±45]s。
5.根据权利要求4所述的汽车复合材料保险杠防撞梁厚度的多目标优化方法，其特征在于，在所述步骤二中，通过最优拉丁超方试验生成初始样本。
6.根据权利要求5所述的汽车复合材料保险杠防撞梁厚度的多目标优化方法，其特征在于，通过最优拉丁超方试验生成初始样本数量为50个。
CN202110737146.4A 2021-06-30 2021-06-30 一种汽车复合材料保险杠防撞梁厚度的多目标优化方法 Active CN113408059B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110737146.4A CN113408059B (zh) 2021-06-30 2021-06-30 一种汽车复合材料保险杠防撞梁厚度的多目标优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110737146.4A CN113408059B (zh) 2021-06-30 2021-06-30 一种汽车复合材料保险杠防撞梁厚度的多目标优化方法

Publications (2)

Publication Number Publication Date
CN113408059A CN113408059A (zh) 2021-09-17
CN113408059B true CN113408059B (zh) 2022-06-07

Family

ID=77680560

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110737146.4A Active CN113408059B (zh) 2021-06-30 2021-06-30 一种汽车复合材料保险杠防撞梁厚度的多目标优化方法

Country Status (1)

Country Link
CN (1) CN113408059B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115481493B (zh) * 2022-10-10 2024-04-30 南京林业大学 一种基于耐撞性的变截面长玻纤增强聚丙烯保险杠横梁优化方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109145490A (zh) * 2018-09-10 2019-01-04 北京理工大学 基于数据特征挖掘的行星进入可达集最优子集计算方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102390335A (zh) * 2011-09-30 2012-03-28 力帆实业(集团)股份有限公司 前保险杠防撞梁及其加工方法
CN104229155A (zh) * 2013-06-08 2014-12-24 北京航空航天大学 考虑屈曲约束与铺层顺序的复合材料机翼壁板铺层确定方法
CN104239952B (zh) * 2014-08-26 2016-02-24 西北工业大学 复合材料层合板纤维连续削层一维优化基础模型及方法
CN107729701A (zh) * 2017-11-23 2018-02-23 江苏理工学院 一种汽车防撞梁铺层优化设计方法
FR3077219B1 (fr) * 2018-01-31 2020-06-26 Saint-Gobain Glass France Vitrage feuillete a dispositif electrocommandable et fabrication
CN108595787A (zh) * 2018-04-04 2018-09-28 江苏理工学院 一种frp防撞梁结构/材料协同优化方法
CN109409013B (zh) * 2018-12-10 2023-02-24 国电联合动力技术有限公司 一种低风速风电机组风轮智能优化设计方法
CN110046378B (zh) * 2019-02-28 2022-09-13 昆明理工大学 一种基于进化多目标优化的选择性分层集成高斯过程回归软测量建模方法
CN110110480B (zh) * 2019-05-21 2020-09-11 吉林大学 一种考虑铺层相容性的碳纤维保险杠防撞梁结构优化设计方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109145490A (zh) * 2018-09-10 2019-01-04 北京理工大学 基于数据特征挖掘的行星进入可达集最优子集计算方法

Also Published As

Publication number Publication date
CN113408059A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
CN108985335B (zh) 核反应堆包壳材料辐照肿胀的集成学习预测方法
CN102185735B (zh) 一种网络安全态势预测方法
EP1598751B1 (en) Estimation of distribution algorithm (EDA)
CN112101630A (zh) 一种薄壁塑件注射成型工艺参数多目标优化方法
CN102694800A (zh) 网络安全态势预测的高斯过程回归方法
CN106778846A (zh) 一种基于支持向量机的风电场风速预测方法
CN113408059B (zh) 一种汽车复合材料保险杠防撞梁厚度的多目标优化方法
CN103440493A (zh) 基于相关向量机的高光谱影像模糊分类方法及装置
Hota Diagnosis of breast cancer using intelligent techniques
Fong et al. A novel feature selection by clustering coefficients of variations
Qiao et al. Identification of fuzzy neural networks by forward recursive input-output clustering and accurate similarity analysis
CN116822214A (zh) 一种基于异方差高斯过程回归的多可信度数据融合方法
CN113836823A (zh) 一种基于负荷分解和优化双向长短期记忆网络的负荷组合预测方法
Zadeh et al. An efficient Bi-level hybrid multi-objective reliability-based design optimization of composite structures
Wignall et al. Comparisons of performance metrics and machine learning methods on an entry descent and landing database
Mashrei Neural network and adaptive neuro-fuzzy inference system applied to civil engineering problems
Cheraghi et al. An adaptive artificial neural network for reliability analyses of complex engineering systems
CN110363232B (zh) 基于bp神经网络的毫米波探测器干扰效果评估方法
Shemyakin et al. Online identification of large-scale chaotic system
CN115166811A (zh) 一种模拟工作现场中子谱快速构建方法
CN116663404A (zh) 一种耦合人工智能和贝叶斯理论的洪水预报方法及系统
CN106203520A (zh) 基于深度相关向量机的sar图像分类方法
Roushangar et al. Novel application of robust GWO-KELM model in predicting discharge coefficient of radial gates: a field data-based analysis
Min et al. Machine learning for predicting fatigue properties of additively manufactured materials
Tovias-Alanis et al. Learning smooth dendrite morphological neurons for pattern classification using linkage trees and evolutionary-based hyperparameter tuning

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