CN113408046B - 一种基于耐撞性的铆接点变间距布置的多学科优化方法 - Google Patents
一种基于耐撞性的铆接点变间距布置的多学科优化方法 Download PDFInfo
- Publication number
- CN113408046B CN113408046B CN202110549245.XA CN202110549245A CN113408046B CN 113408046 B CN113408046 B CN 113408046B CN 202110549245 A CN202110549245 A CN 202110549245A CN 113408046 B CN113408046 B CN 113408046B
- Authority
- CN
- China
- Prior art keywords
- riveting
- points
- optimization
- design
- model
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000013461 design Methods 0.000 claims abstract description 86
- 238000012549 training Methods 0.000 claims abstract description 40
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 238000012795 verification Methods 0.000 claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 12
- 238000005070 sampling Methods 0.000 claims abstract description 11
- 230000004044 response Effects 0.000 claims description 62
- 238000004088 simulation Methods 0.000 claims description 23
- 239000003795 chemical substances by application Substances 0.000 claims description 22
- 239000013598 vector Substances 0.000 claims description 16
- 238000010521 absorption reaction Methods 0.000 claims description 13
- 238000005452 bending Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000002790 cross-validation Methods 0.000 claims description 6
- 230000002068 genetic effect Effects 0.000 claims description 6
- 238000010200 validation analysis Methods 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 abstract description 8
- 238000003466 welding Methods 0.000 description 7
- 239000011295 pitch Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 101100421300 Schizosaccharomyces pombe (strain 972 / ATCC 24843) set8 gene Proteins 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 229910000679 solder Inorganic materials 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/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Automation & Control Theory (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Biophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于耐撞性的铆接点变间距布置的多学科优化方法，属于运载装备结构设计技术领域。基于铆接点变间距布置优化问题的定义，建立以一系列铆接点间距为设计变量的参数化有限元模型，以混料试验设计方法采集训练样本，以随机取样方法选取验证样本，使用留一法交叉验证策略选择代理模型，在完成代理模型精度的验证后，采用启发式多目标优化算法获得优化方案并验证，进而完成铆接点变间距布置的优化设计。本发明的技术方案解决了现有技术中的铆接点优化方法无法实现多学科载荷条件下铆接点变间距布置设计优化的问题，可在保证结构性能的前提下，同时优化铆接点数目和铆接点布置位置。
Description
技术领域
本发明属于运载装备结构设计技术领域，涉及一种铆接点布置多学科优化方法，特别涉及一种基于耐撞性的铆接点变间距布置的多学科优化方法。
背景技术
现今，各种牌号的铝合金、高强钢以及复合材料在汽车、航空、航天等领域已得到广泛应用，其应用方式也从采用某种单一材料发展到采用多种材料的综合制造。材料的改变带来了连接方式的改变，早期的焊点连接已不适合连接异质材料，铆接成为连接异质材料的最主要方式。航空航天和汽车领域产品因其结构尺寸大、性能要求高导致其一种结构上就具有几千个铆接点，而铆接点的多少和位置稀疏直接影响到结构的安全性和轻量化。基于此，本发明开展了一种基于耐撞性的铆接点变间距布置的多学科优化方法研究。
传统的铆接点布置方案多是基于工程经验给出，存在铆接点不足导致的性能缺陷和铆接点过多导致的成本增加等问题。拓扑优化技术是最常用且有效的铆接点布置设计方法之一，但是由于其理论基础的限制，不适用于考虑碰撞等非线性工况的多学科优化问题。
基于代理模型技术的尺寸优化是目前用于铆接点布置多学科优化设计的最先进方法。该方法以铆接点布置几何参数为设计变量，以运载装备结构性能要求和铆接点间距为约束条件，以铆接点数目最少或结构刚强度最大化为优化目标。在定义铆接点布置优化问题后，通过合适的试验设计方法获取训练和验证样本点，创建合理的代理模型并验证其精度，采用启发式算法寻找最优设计方案并验证结果的准确性。因此，铆接点布置优化问题的定义，试验设计方法和代理模型的准确性直接影响铆接点布置设计的优化效果。
铆接点与焊点均采用离散布置的形式连接部件，且均可直接影响结构性能与制造成本，其设计方法相似。中国专利号：201410214771.0，专利名称为：车身B柱焊点布置优化方法，发明人崔岸等提出了一种基于有限元模型与连续二次规划寻优算法的焊点布局设计方法，该方法通过将现有的焊点布置尺寸优化采用分区等间距布置参数为设计变量，即每条焊线内等间距布置，不同焊线采用不同的焊点间距变量表示，并采用直接寻优的方式进行焊点布置设计。
这类专利中的设计方案一般基于等间距布置或分区等间距布置展开，并采用给定的代理模型，甚至直接寻优方法进行布置设计，基于此，现有专利存在以下不足：一是未考虑单一铆接线内铆接点变间距布置对结构性能的提升；二是代理模型基于工程经验给定，难以推广到多学科优化设计应用场景；三是直接寻优迭代次数多，不适用于碰撞等仿真计算成本高的非线性工况。
发明内容
根据上述提出的现有铆接点布置尺寸优化方法无法实现同一铆接线内变间距布置，以及现有代理模型选择方法和直接寻优方法不适用于铆接点布置多学科优化设计的技术问题，而提供一种基于耐撞性的铆接点变间距布置的多学科优化方法。本发明主要利用基于铆接点变间距布置优化问题的定义，建立以一系列铆接点间距为设计变量的参数化有限元模型，以混料试验设计方法采集训练样本，以随机取样方法选取验证样本，使用留一法交叉验证策略选择代理模型，在完成代理模型精度的验证后，采用启发式多目标优化算法获得优化方案并验证，进而完成铆接点变间距布置的优化设计，从而起到可在保证结构多学科性能的前提下，同时优化铆接点数目和铆接点布置位置的优点。
为了达到上述目的，本发明采用的技术方案如下：
一种基于耐撞性的铆接点变间距布置的多学科优化方法，包括以下步骤：
S1、定义铆接点变间距布置的多学科优化问题，具体步骤如下：
S11、设计变量的定义，对于单条铆接线，给定铆接点数目N，以铆接点间距作为设计变量，考虑边缘铆接点与两端的距离以及相邻铆接点间距，设计变量数为n＝N+1，且变量之和等于铆接线长度L；多条铆接线应分别定义；若考虑对称性，可合并相互关联的设计变量；
S12、约束条件的定义，采用有限元仿真方法计算原有铆接点布置方案对应的结构多学科性能响应列向量Y，其中Y＝(y1,y2,y3,y4,y5)T，y1,y2,y3,y4,y5分别对应结构吸能量EA、最大碰撞力峰值Fmax、弯曲刚度Kb、扭转刚度Kt以及一阶模态响应f1的值，并将其分别放松5％作为设计基线；
S13、优化目标的确定，采用结构碰撞安全性能评价指标：结构吸能量最大化和碰撞最大峰值力最小化，优化列式如公式(1)所示：
其中，X表示铆接点布置方案，其包含n个设计变量x1,x2,…,xn；xi是第i个铆接点设计变量，i＝1,n时分别对应的是铆接线上第一个铆接点和最后一个铆接点到结构两端的距离，i＝2,3,…,n-1时对应的是相邻铆接点间距；xmin和xmax分别为铆接点端距的最小值和最大值，x'min和x'max为相邻铆接点间距的最小值和最大值，其与等式约束一起决定了铆接点布置设计域；目标函数-EA(X)为结构吸能量的相反数，Fmax(X)为最大碰撞力峰值；Kb0、Kt0、f10、EA0、Fmax0分别是弯曲刚度下限值、扭转刚度下限值、模态响应下限值、吸能量下限值和最大峰值力上限值。
S2、基于步骤S12中的有限元仿真模型和步骤S13中的优化问题定义，建立参数化有限元模型，具体步骤如下：
S21、基于步骤S12中的有限元仿真模型，删除其中的连接单元后获得只包含部件的有限元模型；基于步骤S13中的优化问题定义，建立铆接点间距和数目可变的参数化有限元模型；
S3、采用混料试验设计方法选取t个铆接点布置方案TL＝(X1,X2,…,Xt)，基于步骤S2建立的参数化有限元模型分别计算各铆接点布置方案Xa(a＝1,2,…,t)对应的性能响应列向量Ya(a＝1,2,…,t)，获得包含t个样本点的训练集采用随机取样方法选取v个铆接点布置方案VL＝(Xt+1,Xt+2,…,Xt+v)，基于步骤S2建立的参数化有限元模型分别计算各铆接点布置方案Xb(b＝t+1,t+2,…,t+v)对应的性能响应列向量Yb(b＝t+1,t+2,…,t+v)，获得包含v个样本点的验证集/>具体包括以下步骤：
S31、基于混料试验设计方法的训练集采样，具体包括：
S311、在设计域中随机确定t个的初始铆接点布置方案；
S312、添加可增加协方差矩阵的秩的额外样本点，直至得到满秩的协方差矩阵；
S313、采用D-最优算法使样本点集对应的协方差矩阵行列式最小，计算各个点对最优判据的贡献值；
S314、以步骤S313获得的贡献值作为排序依据，从高到低排列，从最后一个点开始，将其沿若干方向移动一定的距离并重新评价其对最优判据的贡献；
S315、如果移动后其对最优判据的贡献增加，保留移动后的点并更新；否则，该点保持不变；
S316、重新执行更新过程，直至遍历样本点，如果样本点均保持不变则得到一个局部最优的点集；
S317、重复步骤S311-S316数次，获取t个铆接点布置方案TL＝(X1,X2,…,Xt)；
S318、基于步骤S2建立的参数化有限元模型，分别计算各铆接点布置方案Xa(a＝1,2,…,t)对应的性能响应列向量Ya(a＝1,2,…,t)，获得训练集
S32、基于随机取样方法的验证集采样，具体包括：
S321、在步骤S13定义的设计域中随机生成r个前n-1个设计变量组成的向量其中c＝1,2,…,r；
S322、基于公式(1)中计算向量Z(c)对应的/>
S323、选取v个满足公式(1)中设计变量约束条件的铆接点布置方案VL＝(Xt+1,Xt+2,…,Xt+v)，其中
S324、基于步骤S2建立的参数化有限元模型分别计算各铆接点布置方案Xb(b＝t+1,t+2,…,t+v)对应的性能响应列向量Yb(b＝t+1,t+2,…,t+v)，获得验证集
S4、基于步骤S31获得的训练集TS，采用留一法交叉验证策略选择代理模型，具体步骤如下：
S41、基于步骤S31获得的训练集从中任取一个样本点作为测试点，剩余的t-1个样本点作为训练点建立代理模型，采用公式(2)计算该测试点处各种性能响应yk(k＝1,2,…,5)的预测误差Ejk；
其中，yjk为该测试点对应的第k种性能响应的仿真结果，为采用除样本点TSj外的t-1个样本点建立的代理模型在测试点TSj处对应的第k种性能响应的预测结果。
S42、重复步骤S41，遍历整个训练集，使t个训练样本点依次单独作为测试点，并以公式(3)计算各种性能响应yk(k＝1,2,…,5)的泛化平方差GMSEk；
S43、更改代理模型类型，重复步骤S41与S42；
S44、针对各性能响应yk(k＝1,2,…,5)，分别独立选取可使其泛化平方差GMSEk最小的代理模型类型用于后续优化流程；
S5、基于步骤S4中为各性能响应yk(k＝1,2,…,5)分别独立选取的代理模型类型，验证代理模型拟合精度和预测精度，具体步骤如下：
S51、基于步骤S3获得的t个训练样本点，计算代理模型的拟合精度，如计算公式(4)所示，
其中，表示第k种性能响应对应的代理模型的拟合精度；/>为训练集中第k种性能响应的平均值；/>是采用t个训练样本点建立的代理模型在训练样本点TSj处第k种性能响应的预测结果；yjk表示在训练样本点TSj处第k种性能响应的仿真结果。若拟合精度大于0.95则认为满足要求；若拟合精度小于0.95则应返回步骤S3添加额外的训练样本点后重新计算。
S52、基于步骤S3获得的v个验证样本点，计算代理模型的预测误差，分别为公式(5)和公式(6)所示：
其中，MREk为第k种性能响应对应的代理模型的最大相对误差，RMSEk为第k种性能响应对应的代理模型的均方根误差；ymk表示在验证样本点处第k种性能响应的仿真结果；/>是采用t个训练样本点建立的代理模型在验证样本点VSm处第k种性能响应的预测结果。
如最大相对误差MREk(k＝1,2,…,5)小于10％，且均方根误差RMSEk(k＝1,2,…,5)显著小于响应值，则认为代理模型预测精度满足要求，代理模型可用于后续的铆接点布置优化求解；否则，精度不满足要求，应返回步骤S3，添加训练样本点后重新计算；
S6、基于步骤S13建立的优化列式，步骤S2建立的参数化有限元模型以及步骤S5验证过的代理模型，采用遗传算法获得优化方案，并对优优化方案进行验证分析，具体步骤如下：
S61、基于步骤S13中的优化列式以及步骤S5验证过的代理模型，采用遗传算法获得优化解集；
S62、从优化解集中随机选u个最优设计方案PL＝(Xt+v+1,Xt+v+2,…,Xt+v+u)；
S63、基于步骤S2所建立的参数化有限元模型，计算各个最优设计方案Xp(p＝t+v+1,t+v+2,…,t+v+u)对应的结构多学科性能响应Yp(p＝t+v+1,t+v+2,…,t+v+u)；
S64、基于步骤S63得到的有限元仿真结果与步骤S61得到的代理模型预测结果，计算其相对误差，若误差在5％以内，则认为优化解集精度满足要求；否则，优化解集不满足精度要求，应返回步骤S3，添加训练样本点后重新计算；
S65、若优化解集精度满足要求，且结构性能响应满足公式(1)中的约束条件，则优化解集可用于指导运载装备铆接点布置设计；若优化解集精度满足要求，但结构性能响应不满足公式(1)中的约束条件，则铆接点数目应取N＝N+1，并重新从步骤S1开始迭代。
相较于现有技术，本发明具有以下有益效果：
1)本发明中提出的基于耐撞性的铆接点变间距布置的多学科优化方法中关于铆接点布置设计变量的定义，使单一铆接线内铆接点变间距布置对结构性能的提升得到考虑，可进一步提升结构性能。
2)本发明中提出的基于耐撞性的铆接点变间距布置的多学科优化方法的使用，使设计变量的等式约束条件以及铆接工艺不等式约束条件同时得到考虑，实现单一铆接线内铆接点变间距布置试验设计。
3)本发明中提出的基于耐撞性的铆接点变间距布置的多学科优化方法中留一法交叉验证策略的使用，可避免代理模型选择阶段所需的工程经验或试错过程，且不受多学科工况变化的局限，提升了设计效率和结果的准确性。
基于上述理由本发明可在运载装备结构铆接点布置设计等领域广泛推广。
附图说明
图1为本发明流程图。
图2为本发明铆接点变间距布置设计变量示意图。
图3为本发明多学科工况示意图。
图4为本发明优化方案示意图。
图中：1轴向冲击工况；2弯曲工况；3扭转工况；4模态分析；图4(a)铆接点布置设计方案1；图4(b)铆接点布置设计方案2。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
本发明基于铆接点变间距布置优化问题的定义，建立以一系列铆接点间距为设计变量的参数化有限元模型；以混料试验设计方法采集训练样本；以随机取样方法选取验证样本；考虑的多学科工况包括弯曲、扭转、模态分析以及冲击载荷；以原有铆接点布置方案对应的结构多学科性能响应作为优化问题中约束条件设置的参考，保证优化后的结构性能；使用留一法交叉验证策略选择代理模型，并基于验证样本计算代理模型精度；基于满足精度要求的代理模型，采用启发式多目标优化算法获得优化方案；最后，基于优化后的铆接点布置方案，通过有限元分析验证优化方案是否可行。
结合附图和技术方案本实施例详细说明了本发明提供的一种基于耐撞性的铆接点变间距布置的多学科优化方法的具体实施方式，方法流程如图1所示。
本发明针对目前铆接点尺寸优化中未考虑单一铆接线内铆接点变间距布置对结构性能的提升，采用混料试验设计方法选取满足变间距布置约束条件的训练样本点，使得同一铆接线内的铆接点间距可随机变化。本发明采用铆接点变间距多学科优化方法，基于代理模型进行铆接点数目和布置形式变化后结构性能的预测，通过使用交叉验证的方法选择最合适的代理模型。在满足结构多学科性能的前提下，以结构耐撞性作为优化目标，最终获得满足性能要求的铆接点布置方案。本发明适用于运载装备结构铆接点布置优化设计，下面以双帽形梁为例具体说明本铆接点优化设计方法：
S1、定义铆接点变间距布置的多学科优化问题，具体步骤包括：
S11、设计变量的定义，针对双帽形梁的单条铆接线，给定铆接点数目迭代初始值N＝8；以铆接点间距作为设计变量，考虑到帽形梁的对称性，合并相互关联的设计变量；考虑边缘铆接点与两端的距离以及相邻铆接点间距，设计变量数为n＝9，且变量之和等于铆接线长度L＝350mm，如图2所示；
S12、约束条件的定义，双帽型梁翻边宽度为30mm，长方形梁截面的尺寸为1000mm*80mm，多学科工况如图3所示，包括结构碰撞冲击工况1、弯曲工况2、扭转工况3以及模态分析4，其中结构碰撞冲击工况1为刚性墙以15m/s的速度匀速下压240mm，弯曲与扭转工况加载单位载荷；选用的材料为DP590高强钢。采用有限元仿真技术计算两侧各12个铆接点以30mm等间距布置的方案对应的结构多学科性能响应列向量Y，得到Y＝(32.62,345.70,1611.72,276.70,419.74)T，分别为结构吸能量EA，最大碰撞力峰值Fmax，弯曲刚度Kb、扭转刚度Kt以及一阶模态响应f1的值，并将各性能响应分别放松5％作为设计基线，结果如公式(7)所示；
S13、优化目标的确定，采用结构碰撞安全性能评价指标：结构吸能量相反数和最大碰撞峰值力最小化，优化列式如下，
其中，X表示铆接点布置方案，其包含9个设计变量x1,x2,…,x9；xi是第i个铆接点设计变量，i＝1,9时分别对应的是铆接线上第一个铆接点和最后一个铆接点到结构两端的距离，i＝2,3,…,8时对应的是相邻铆接点间距；铆接点端距的最小值和最大值分别为10mm和35mm，相邻铆接点间距的最小值和最大值分别为10mm和47.8mm，其与等式约束一起决定了铆接点布置设计域；目标函数-EA(X)为结构吸能量的相反数，Fmax(X)为最大碰撞力峰值；弯曲刚度下限值为1563.36N/mm、扭转刚度下限值为268.40N·m/deg、模态响应下限值为407.15Hz、吸能量下限值为31.64kJ和最大峰值力上限值为356.07kN。
S2、基于步骤S12中的有限元仿真模型和步骤S13中的优化问题定义，建立参数化有限元模型，具体步骤包括：
S21、基于步骤S12中的有限元仿真模型，删除其中的连接单元后获得只包含部件的有限元模型；基于步骤S13中的优化问题定义，建立铆接点间距和数目可变的参数化有限元模型；
S3、采用混料试验设计方法选取选取78个铆接点布置方案TL＝(X1,X2,…,X78)，基于步骤S2建立的参数化有限元模型分别计算各铆接点布置方案Xa(a＝1,2,…,78)对应的性能响应列向量Ya(a＝1,2,…,78)，获得包含78个样本点的训练集采用随机取样方法选取8个铆接点布置方案VL＝(X79,X80,…,X86)，基于步骤S2建立的参数化有限元模型分别计算各铆接点布置方案Xb(b＝79,80,…,86)对应的性能响应列向量Yb(b＝79,80,…,86)，获得包含8个样本点的验证集
S4、基于步骤S3获得的训练集TS，使用留一法交叉验证策略选择代理模型，本案例中采用3种代理模型，包括多项式响应面PRS，径向基神经网络RBNN，以及克里金KRG代理模型，用于预测5种结构多学科性能响应。
基于步骤S3获得的78个训练样本点，执行步骤S41-S43，计算得3种代理模型预测5种不同性能响应的泛化平方误差，结果如表1所示，
表1留一法交叉验证结果(铆接点数目N＝8)
基于表中数据，执行步骤S44，针对各性能响应yk(k＝1,2,…,5)，分别独立选取可使其泛化平方差GMSEk最小的代理模型类型。其中，吸能量EA选用KRG，最大峰值力Fmax、弯曲刚度Kb、扭转刚度Kt、以及模态响应f1选用PRS，这些代理模型类型被确定用于后续的优化流程；
S5、基于步骤S4中为各性能响应yk(k＝1,2,…,5)分别独立选取的代理模型类型，验证代理模型拟合精度和预测精度，执行步骤S51和步骤S52计算5种性能响应的拟合精度和预测精度，结果如表2所示
表2代理模型拟合精度与预测精度
基于表中数据，最大相对误差MRE均小于5％，且均方根误差RMSE显著小于响应值的数量级，可认为代理模型预测精度满足要求，代理模型可用于后续的铆接点布置优化求解；
S6、基于步骤S13建立的优化列式，步骤S2建立的参数化有限元模型以及步骤S5验证过的代理模型，采用遗传算法获得优化方案，并对优优化方案进行验证分析，具体步骤如下：
S61、基于步骤S13中的优化列式以及步骤S5验证过的代理模型，采用遗传算法获得优化解集；
S62、从优化解集中随机选2个最优设计方案，后续以设计方案1、设计方案2表示，如图4所示，图4(a)为设计方案1，铆接点间距分别为x1＝10.04mm,x2＝36.95mm,x3＝46.51mm,x4＝41.67mm,x5＝47.80mm,x6＝36.45mm,x7＝47.80mm,x8＝47.80mm,x9＝34.98mm；图4(b)为设计方案2，铆接点间距分别为x1＝10.00mm,x2＝45.18mm,x3＝46.14mm,x4＝38.28mm,x5＝47.80mm,x6＝36.69mm,x7＝43.63mm,x8＝47.80mm,x9＝34.48mm；
S63、基于步骤S2所建立的参数化有限元模型，分别计算设计方案1和设计方案2对应的结构多学科性能响应；
S64、基于步骤S61得到的代理模型预测结果与步骤S63得到的有限元仿真结果，计算其相对误差，如表3所示；
表3部分优化方案对应的预测结果与仿真结果
基于表3中数据，相对误差在5％以内，认为优化解集精度满足要求；
S65、优化解集精度满足要求，且结构性能响应满足公式(1)中的约束条件，则优化解集可用于指导运载装备铆接点布置设计；铆接点优化流程终止迭代。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (4)
1.一种基于耐撞性的铆接点变间距布置的多学科优化方法，其特征在于，包括以下步骤：
S1、定义铆接点变间距布置的多学科优化问题
S11、设计变量的定义；
S12、约束条件的定义，采用有限元仿真方法计算原有铆接点布置方案对应的结构多学科性能响应列向量Y，其中Y＝(y1,y2,y3,y4,y5)T，y1,y2,y3,y4,y5分别对应结构吸能量EA、最大碰撞力峰值Fmax、弯曲刚度Kb、扭转刚度Kt以及一阶模态响应f1的值，并将其分别放松5％作为设计基线；
S13、优化目标的确定，采用结构碰撞安全性能评价指标：结构吸能量最大化和碰撞最大峰值力最小化，优化列式如公式(1)所示：
其中，X表示铆接点布置方案，其包含n个设计变量x1,x2,…,xn；xi是第i个铆接点设计变量，i＝1,n时分别对应的是铆接线上第一个铆接点和最后一个铆接点到结构两端的距离，i＝2,3,…,n-1时对应的是相邻铆接点间距；xmin和xmax分别为铆接点端距的最小值和最大值，x'min和x'max为相邻铆接点间距的最小值和最大值；目标函数-EA(X)为结构吸能量的相反数，Fmax(X)为最大碰撞力峰值；Kb0、Kt0、f10、EA0、Fmax0分别是弯曲刚度下限值、扭转刚度下限值、模态响应下限值、吸能量下限值和最大峰值力上限值；
S2、基于步骤S12中的有限元仿真模型和步骤S13中的优化问题定义，建立参数化有限元模型：
S21、基于步骤S12中的有限元仿真模型，删除其中的连接单元后获得只包含部件的有限元模型；基于步骤S13中的优化问题定义，建立铆接点间距和数目可变的参数化有限元模型；
S3、采用混料试验设计方法选取t个铆接点布置方案TL＝(X1,X2,…,Xt)，基于参数化有限元模型分别计算各铆接点布置方案Xa(a＝1,2,…,t)对应的性能响应列向量Ya(a＝1,2,…,t)，获得包含t个样本点的训练集采用随机取样方法选取v个铆接点布置方案VL＝(Xt+1,Xt+2,…,Xt+v)，基于步骤S2建立的参数化有限元模型分别计算各铆接点布置方案Xb(b＝t+1,t+2,…,t+v)对应的性能响应列向量Yb(b＝t+1,t+2,…,t+v)，获得包含v个样本点的验证集/>
S4、基于训练集TS，采用留一法交叉验证策略选择代理模型，具体如下：
S41、基于训练集从中任取一个样本点/>(j＝1,2,…,t)作为测试点，剩余的t-1个样本点作为训练点建立代理模型，采用公式(2)计算该测试点处各种性能响应yk(k＝1,2,…,5)的预测误差Ejk；
其中，yjk为该测试点对应的第k种性能响应的仿真结果，为采用除样本点TSj外的t-1个样本点建立的代理模型在测试点TSj处对应的第k种性能响应的预测结果；
S42、重复步骤S41，遍历整个训练集，使t个训练样本点依次单独作为测试点，并以公式(3)计算各种性能响应yk(k＝1,2,…,5)的泛化平方差GMSEk；
S43、更改代理模型类型，重复步骤S41与S42；
S44、针对各性能响应yk(k＝1,2,…,5)，分别独立选取可使其泛化平方差GMSEk最小的代理模型类型用于后续优化流程；
S5、基于步骤S4中为各性能响应yk(k＝1,2,…,5)分别独立选取的代理模型类型，验证代理模型拟合精度和预测精度，具体如下：
S51、基于步骤S3获得的t个训练样本点，计算代理模型的拟合精度，如计算公式(4)所示，
其中，表示第k种性能响应对应的代理模型的拟合精度；/>为训练集中第k种性能响应的平均值；/>是采用t个训练样本点建立的代理模型在训练样本点TSj处第k种性能响应的预测结果；yjk表示在训练样本点TSj处第k种性能响应的仿真结果；若拟合精度满足要求；若拟合精度不满足要求，则应返回步骤S3添加额外的训练样本点后重新计算；
S52、基于v个验证样本点，计算代理模型的预测误差，分别为公式(5)和公式(6)所示：
其中，MREk为第k种性能响应对应的代理模型的最大相对误差，RMSEk为第k种性能响应对应的代理模型的均方根误差；ymk表示在验证样本点(m＝t+1,t+2,…,t+v)处第k种性能响应的仿真结果；/>是采用t个训练样本点建立的代理模型在验证样本点VSm处第k种性能响应的预测结果；
如最大相对误差MREk(k＝1,2,…,5)小于10％，且均方根误差RMSEk(k＝1,2,…,5)显著小于响应值，则认为代理模型预测精度满足要求，代理模型可用于后续的铆接点布置优化求解；否则，精度不满足要求，应返回步骤S3，添加训练样本点后重新计算；
S6、基于步骤S13建立的优化列式，步骤S2建立的参数化有限元模型以及步骤S5验证过的代理模型，采用遗传算法获得优化方案，并对优优化方案进行验证分析，具体如下：
S61、基于步骤S13中的优化列式以及步骤S5验证过的代理模型，采用遗传算法获得优化解集；
S62、从优化解集中随机选u个最优设计方案PL＝(Xt+v+1,Xt+v+2,…,Xt+v+u)；
S63、基于步骤S2所建立的参数化有限元模型，计算各个最优设计方案Xp(p＝t+v+1,t+v+2,…,t+v+u)对应的结构多学科性能响应Yp(p＝t+v+1,t+v+2,…,t+v+u)；
S64、基于步骤S63得到的有限元仿真结果与步骤S61得到的代理模型预测结果，计算其相对误差，若误差在5％以内，则认为优化解集精度满足要求；否则，优化解集不满足精度要求，应返回步骤S3，添加训练样本点后重新计算；
S65、若优化解集精度满足要求，且结构性能响应满足公式(1)中的约束条件，则优化解集可用于指导运载装备铆接点布置设计；若优化解集精度满足要求，但结构性能响应不满足公式(1)中的约束条件，则铆接点数目应取N＝N+1，并重新从步骤S1开始迭代。
2.根据权利要求1所述的一种基于耐撞性的铆接点变间距布置的多学科优化方法，其特征在于，所述的步骤S11中，对于单条铆接线，给定铆接点数目N，以铆接点间距作为设计变量，考虑边缘铆接点与两端的距离以及相邻铆接点间距，设计变量数为n＝N+1，且变量之和等于铆接线长度L；多条铆接线应分别定义；若考虑对称性，可合并相互关联的设计变量。
3.根据权利要求1所述的一种基于耐撞性的铆接点变间距布置的多学科优化方法，其特征在于，所述的步骤3具体如下：
S31、基于混料试验设计方法的训练集采样；
S32、基于随机取样方法的验证集采样，具体包括：
S321、在步骤S13定义的设计域中随机生成r个前n-1个设计变量组成的向量其中c＝1,2,…,r；
S322、基于公式(1)中计算向量Z(c)对应的/>
S323、选取v个满足公式(1)中设计变量约束条件的铆接点布置方案VL＝(Xt+1,Xt+2,…,Xt+v)，其中b＝t+1,t+2,…,t+v；
S324、基于步骤S2建立的参数化有限元模型分别计算各铆接点布置方案Xb(b＝t+1,t+2,…,t+v)对应的性能响应列向量Yb(b＝t+1,t+2,…,t+v)，获得验证集
4.根据权利要求1所述的一种基于耐撞性的铆接点变间距布置的多学科优化方法，其特征在于，所述的步骤S51中，若拟合精度大于0.95则认为满足要求，否则认为不满足要求。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110549245.XA CN113408046B (zh) | 2021-05-20 | 2021-05-20 | 一种基于耐撞性的铆接点变间距布置的多学科优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110549245.XA CN113408046B (zh) | 2021-05-20 | 2021-05-20 | 一种基于耐撞性的铆接点变间距布置的多学科优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113408046A CN113408046A (zh) | 2021-09-17 |
| CN113408046B true CN113408046B (zh) | 2023-09-22 |
Family
ID=77679185
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110549245.XA Active CN113408046B (zh) | 2021-05-20 | 2021-05-20 | 一种基于耐撞性的铆接点变间距布置的多学科优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113408046B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114117763B (zh) * | 2021-11-18 | 2024-05-03 | 浙江理工大学 | 一种基于移一交叉验证法的模型筛选方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107256289A (zh) * | 2017-05-18 | 2017-10-17 | 吉林大学 | 汽车碰撞简化参数化有限元模型的建立方法 |
| CN109063357A (zh) * | 2018-08-15 | 2018-12-21 | 大连理工大学 | 基于拓扑优化的自适应多步变域的汽车构件焊点布局优化方法 |
| CN109933836A (zh) * | 2019-01-03 | 2019-06-25 | 重庆长安汽车股份有限公司 | 一种基于车身性能约束的白车身焊点优化布置方法 |
| CN111898291A (zh) * | 2020-06-30 | 2020-11-06 | 西北工业大学 | 一种基于子结构的大壁板连续铆接变形预测方法 |
| WO2020244325A1 (zh) * | 2019-06-04 | 2020-12-10 | 南京依维柯汽车有限公司 | 一种基于相对灵敏度的非承载式白车身轻量化方法 |
-
2021
- 2021-05-20 CN CN202110549245.XA patent/CN113408046B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107256289A (zh) * | 2017-05-18 | 2017-10-17 | 吉林大学 | 汽车碰撞简化参数化有限元模型的建立方法 |
| CN109063357A (zh) * | 2018-08-15 | 2018-12-21 | 大连理工大学 | 基于拓扑优化的自适应多步变域的汽车构件焊点布局优化方法 |
| CN109933836A (zh) * | 2019-01-03 | 2019-06-25 | 重庆长安汽车股份有限公司 | 一种基于车身性能约束的白车身焊点优化布置方法 |
| WO2020244325A1 (zh) * | 2019-06-04 | 2020-12-10 | 南京依维柯汽车有限公司 | 一种基于相对灵敏度的非承载式白车身轻量化方法 |
| CN111898291A (zh) * | 2020-06-30 | 2020-11-06 | 西北工业大学 | 一种基于子结构的大壁板连续铆接变形预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 真空高压铸造铝合金车身后纵梁轻量化设计;林佳武;李玄霜;陈宗明;陈东;李永祥;耿富荣;;汽车工程(第03期);108-114 * |
| 考虑气动阻力和横风稳定的汽车车身多目标优化设计;亓昌;韩元吉;杨姝;吕振华;;汽车安全与节能学报(第01期);57-64 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113408046A (zh) | 2021-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107909179B (zh) | 一种插电式混合动力车辆行驶工况的预测模型构建方法及车辆能量管理方法 | |
| CN102867101B (zh) | 一种确定桁架结构参数的方法 | |
| CN104866692A (zh) | 一种基于自适应代理模型的飞行器多目标优化方法 | |
| CN105260532B (zh) | 基于序列近似优化的薄板拉伸变压边力不确定性设计方法 | |
| CN111709082A (zh) | 一种高效的汽车侧面碰撞安全可靠性设计优化方法 | |
| CN110210155B (zh) | 一种飞机壁板铆接路径数值规划方法 | |
| CN106055791A (zh) | 基于预估校正算法的飞行器全局气动优化方法 | |
| CN113408046B (zh) | 一种基于耐撞性的铆接点变间距布置的多学科优化方法 | |
| CN112182950B (zh) | 一种气体介质绝缘能力评估方法 | |
| CN108319799A (zh) | 一种自主水下航行器的外形多保真度优化设计方法 | |
| CN112036062A (zh) | 一种金属材料弯曲成形回弹角预测方法 | |
| CN116306156B (zh) | 车身优化方法、装置、存储介质及电子设备 | |
| CN111745653B (zh) | 基于双机械臂的船体外板曲面成形协同加工的规划方法 | |
| CN112417666A (zh) | 一种带筋壁板预应力喷丸成形的数值模拟方法 | |
| CN110928261B (zh) | 分布式异构流水车间的分布估计调度方法及系统 | |
| CN112817312A (zh) | 一种基于双重搜索优化算法的路径规划方法 | |
| Tang et al. | Novel solution framework for inverse problem considering interval uncertainty | |
| CN114564787A (zh) | 用于目标相关翼型设计的贝叶斯优化方法、装置及存储介质 | |
| CN111369039B (zh) | 基于多目标优化辅助主动学习的合金热处理工艺优化方法 | |
| CN110210072B (zh) | 基于近似模型及差分进化算法求解高维优化问题的方法 | |
| CN106934186A (zh) | 一种结构可靠性优化求解的融合方法 | |
| CN109326319B (zh) | 一种基于二级结构知识的蛋白质构象空间优化方法 | |
| Wang et al. | An improved system for efficient shape optimization of vehicle aerodynamics with “noisy” computations | |
| CN112016243B (zh) | 一种基于响应曲面的交通流量预测模型参数标定方法 | |
| CN117689001B (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 |