CN113536452A - 一种适用于多材料多厚度结构的优化设计方法 - Google Patents
一种适用于多材料多厚度结构的优化设计方法 Download PDFInfo
- Publication number
- CN113536452A CN113536452A CN202110437799.0A CN202110437799A CN113536452A CN 113536452 A CN113536452 A CN 113536452A CN 202110437799 A CN202110437799 A CN 202110437799A CN 113536452 A CN113536452 A CN 113536452A
- Authority
- CN
- China
- Prior art keywords
- response
- sequence
- design
- calculating
- correlation coefficient
- 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
- 238000013461 design Methods 0.000 title claims abstract description 49
- 239000000463 material Substances 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000004044 response Effects 0.000 claims abstract description 55
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000005457 optimization Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000000694 effects Effects 0.000 claims abstract description 9
- 238000004088 simulation Methods 0.000 claims abstract description 7
- 238000002474 experimental method Methods 0.000 claims abstract description 4
- 238000010606 normalization Methods 0.000 claims abstract description 3
- 238000012545 processing Methods 0.000 claims abstract 2
- 238000010521 absorption reaction Methods 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000000513 principal component analysis Methods 0.000 claims description 8
- 239000013598 vector Substances 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 2
- 230000006872 improvement Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 238000003466 welding Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004026 adhesive bonding Methods 0.000 description 1
- 238000000418 atomic force spectrum Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
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
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种适用于多材料多厚度结构的优化设计方法,包括以下步骤:1)根据待优化的多材料多厚度结构,确定设计变量及其水平值和响应的优化目标;2)设计正交实验进行仿真计算并得到相应的响应数据矩阵;3)计算响应的信噪比;4)将各指标序列按照归一化处理;5)获得比较序列并计算灰色关联系数矩阵;6)计算得到灰色关联度序列;7)对灰色关联度序列进行方差分析,获得显著影响变量及量化每个变量对响应的贡献率;8)根据所确定的因子水平效应表确定最优参数组合,从而实现多目标优化设计。本发明提出一种适用于多材料多厚度结构的优化设计方法,找出最优材料厚度匹配方案,使结构在优化前后的安全性和轻量化得到了显著提升。
Description
技术领域
本发明涉及汽车车身轻量化及结构设计领域,尤其涉及一种适用于多材料多厚度结构的优化设计方法。
背景技术
随着汽车行业的快速发展以及全球能源紧缺,人们对于汽车的燃油消耗和安全方面提出了更高的要求,而汽车轻量化是实现车辆燃油经济性,节能减排的重要措施之一。
对于车身同一部件的不同位置有着不同的性能要求,传统的设计是根据某一位置最大性能要求进行的保守结构设计,这种方法造成材料的浪费,给轻量化设计留下设计空间。随着对轻量化结构的探索,通过采用焊接、熔接、胶粘、螺接和铆接等方式,将不同材料或厚度的结构进行连接成同一部件,从而实现同一部件不同位置的变强度设计。因此,如何合理匹配结构的材料和厚度等变量,对于车身结构的轻量化具有重要意义。
连接结构的各变量的匹配是一个多变量多水平的离散优化设计问题。对于这类结构的设计,传统的做法是一种不断试错的方法,即对某一部件进行不断替换,这类做法造成巨大的实验成本,极大依赖于工程师的经验。同时,传统的离散优化算法是基于离散变量的连续化进行的,这类方法造成很大的计算成本,迭代次数多,效率比较低,同时由于近似模型的引入造成优化结果的精度比较低。针对以上问题,本文提出一种适用于多材料多厚度结构的优化设计方法,找出最优材料厚度匹配方案,使结构在优化前后的安全性和轻量化得到了显著提升。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷,提供一种适用于多材料多厚度结构的优化设计方法。
本发明解决其技术问题所采用的技术方案是:一种适用于多材料多厚度结构的优化设计方法,包括以下步骤:
1)根据待优化的汽车车身的多材料多厚度结构,确定设计变量及其水平值和响应的优化目标;所述设计变量为部件的材料等级和厚度,设计变量的水平值为部件的材料等级和厚度的取值集合;所述响应目标包括峰值加速度最小、总吸能最大和质量最小;
2)根据设计变量及其水平数,设计正交实验进行仿真计算并得到相应的响应数据矩阵;
3)对响应目标峰值加速度、总吸能和总质量进行Taguchi分析,计算响应的信噪比;
4)不同的响应指标的量纲和单位不同,为了保证数据序列间的可比性,将各指标序列按照归一化处理;
5)根据所获得归一化后的信噪比矩阵及其本身的望大质量特性,获得比较序列并计算灰色关联系数矩阵;
6)对响应的灰色关联系数矩阵进行主成分分析,获得各响应的权重比,计算得到灰色关联度序列;
7)对灰色关联度序列进行方差分析,获得显著影响变量及量化每个变量对响应的贡献率,并更新因子水平效应表;
8)根据所确定的因子水平效应表确定最优参数组合,对设计参数组合进行试验或者仿真计算,根据优化结果的改善程度或者迭代次数要求,判断是否需要进一步迭代更新正交表,若不需要则输出最优参数组合,从而实现多目标优化设计。
按上述方案,所述步骤3)中计算响应的信噪比具体如下:
根据不同的目标响应各自的质量特性,响应的信噪比矩阵计算如下:
望大特性(响应值越大越好)信噪比计算公式:
望小特性(响应值越小越好)信噪比计算公式:
按上述方案,所述步骤5)中对于灰色关联系数矩阵的计算应首先确定参考序列,由于这里的灰色计算序列规范到0到1之间的信噪比,属于望大特性,因此相应的理想参考序列可设置为1,从而得出各响应指标序列与理想序列1的灰色关联系数:
其中,xi(0)为理想参考序列,ρ为分辨系数,一般取值为0.5。
按上述方案,所述步骤6)对响应的灰色关联系数矩阵进行主成分分析,具体如下:
将主成分分析运用到响应权重比的测量中,首先计算各响应灰色关联系数序列的相关系数矩阵:
相关系数阵列的特征值和特征向量由下式确定:
(R-λkIn)·Vik=0
因此,可以根据对应于第一特征值的特征向量的平方来计算得到每个响应的权重比。
按上述方案,所述步骤6)根据所计算出的响应的灰色关联系数矩阵及其权重比,可得灰色关联度序列:
按上述方案,所述步骤7)具体如下:
(7.1)为得出各变量的贡献率,对灰色关联度序列进行方差分析,计算如下:
其中,SST是所有因子总的离均差平方和,SSa是某个变量的离均差平方和, ravg是响应系统总的平均灰色关联度,MS是均方,F是显著值,C是贡献率,D 是自由度。
(7.2)建立设计变量因子水平效应表,分析各变量及其水平与灰色关联度的变化规律,确定最优的参数组合,实现多目标问题的设计参数优化。
本发明产生的有益效果是:本发明提出一种适用于多材料多厚度结构的优化设计方法,找出最优材料厚度匹配方案,使结构在优化前后的安全性和轻量化得到了显著提升。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的方法流程图。
图2是本发明实施例的拼焊板台车正面碰撞有限元模型示意图;
图3是本发明实施例的拼焊板结构示意图;
图4是本发明实施例的优化前后碰撞力曲线对比图;
图5是本发明实施例的优化前后比吸能曲线对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明所提出一种适用于多材料多厚度结构的优化设计方法,流程如图1所示,主要包括Taguchi分析,灰色关联分析,主成分分析和方差分析,具体步骤如下:
(1)建立拼焊板的台车正碰数值模型,如图2所示,将拼焊板焊接到台车前端。其中所要优化的目标拼焊板如图3所示,部件1与部件2,部件3与部件4 之间是通过焊缝连接在一起的;部件1与部件3,部件2与部件4是通过焊点连接在一起,并且通过试验验证了模型的有效性。
对于此具体案例研究的具体问题是拼焊板的耐撞性和轻量化,因此,这里选择四个部件的材料等级和厚度作为设计变量,以峰值加速度、总吸能和质量作为响应的优化目标。这里以一次迭代为案例,其相应的设计变量水平取值如下表1所示:
表1:
根据设计变量及其水平数,选择L27(38)正交表作为方案进行仿真计算,得到的仿真结果,响应数据序列Y如表2所示:
表2:
(2)对响应目标峰值加速度、总吸能和总质量进行Taguchi分析,由于峰值加速度和总质量属于望小特性,故采用下面公式(1)计算响应的信噪比:
总吸能属于望大特性,故采用下面公式(2)计算响应的信噪比:
计算各响应数据序列得到信噪比序列,并用公式(3)对三个信噪比响应序列进行归一化。峰值加速度、总吸能和总质量的信噪比和归一化后的序列如表3 所示:
表3:
(3)根据所获得归一化后的信噪比序列及其本身的望大质量特性,获得比较序列为1,并且通过下面公式4计算相应的灰色关联系数如表4所示:
表4:
(4)对三个目标响应对应的灰色关联系数进行主成分分析,以确定权重比。
首先利用公式(5)计算各灰色关联系数的关联系数矩阵:
然后利用公式(6)可确定相关系数阵列的特征值和特征向量。
(R-λkIn)·Vik=0
6)
因此,可以根据对应于第一特征值的特征向量的平方值即可获得每个响应的权重比。特征向量及权重比如表5所示:
表5:
根据公式(7)计算出的各响应的灰色关联系数及其权重比,可得灰色关联度:
(5)利用公式(8)-(10),对上述所求的灰色关联度进行方差分析。
其中SST是所有因子总的离均差平方和,SSa是某个变量的离均差平方和, ravg是响应系统总的平均灰色关联度,MS是均方,F是显著值,C是贡献率,D 是自由度。
获得显著影响变量及量化每个变量对响应的贡献率,如表6所示:
表6:
可以看出影响结构耐撞性最为显著的变量是部件2的厚度,其次分别是部件1的厚度、部件4的材料、部件2的材料、部件3的材料、部件1的材料和部件4的厚度。
结合灰色关联度以及正交表,可获得每个变量水平的效应表,如表7所示:
表7:
(6)根据所确定的因子水平效应表获得最优参数组合为G1E3F2D1C1B1A3H1,对设计参数组合进行试验或者仿真计算,最终的响应优化结果与原始的对比如图4所示,可以看出优化后的结构能大幅度降低峰值加速度同时提高比吸能。同时,由于经过一次迭代,拼焊板结构已经有了很大的改善,故并没有更新正交表。若对于其它具体实例需要反复迭代正交表进行进一步优化设计时,可以将上一步迭代所求得的最优参数组合作为设计的中心,即各变量的水平2,变量水平1和3则需要另外选取,以此进行更新正交表,实现迭代的优化,从而实现更为复杂问题的多目标优化设计。从本实例可以看出,采用本发明的设计方法可以使拼焊板的性能得到显著提升,同时只进行了一次迭代,提高了优化效率。
所描述的拼焊板部件优化实施例仅仅是本发明一个实施例,
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (6)
1.一种适用于多材料多厚度结构的优化设计方法,其特征在于,包括以下步骤:
1)根据待优化的多材料多厚度结构,确定设计变量及其水平值和响应的优化目标;所述设计变量为部件的材料等级和厚度,设计变量的水平值为部件的材料等级和厚度的取值集合;所述响应目标包括峰值加速度最小、总吸能最大和质量最小;
2)根据设计变量及其水平数,设计正交实验进行仿真计算并得到相应的响应数据矩阵;
3)对响应目标峰值加速度、总吸能和总质量进行Taguchi分析,计算响应的信噪比;
4)将各指标序列按照归一化处理;
5)根据所获得归一化后的信噪比矩阵及其本身的望大质量特性,获得比较序列并计算灰色关联系数矩阵;
6)对响应的灰色关联系数矩阵进行主成分分析,获得各响应的权重比,计算得到灰色关联度序列;
7)对灰色关联度序列进行方差分析,获得显著影响变量及量化每个变量对响应的贡献率,并更新因子水平效应表;
8)根据所确定的因子水平效应表确定最优参数组合,对设计参数组合进行试验或者仿真计算,根据优化结果的改善程度或者迭代次数要求,判断是否需要进一步迭代更新正交表,若不需要则输出最优参数组合,从而实现多目标优化设计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110437799.0A CN113536452A (zh) | 2021-04-22 | 2021-04-22 | 一种适用于多材料多厚度结构的优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110437799.0A CN113536452A (zh) | 2021-04-22 | 2021-04-22 | 一种适用于多材料多厚度结构的优化设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113536452A true CN113536452A (zh) | 2021-10-22 |
Family
ID=78095321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110437799.0A Pending CN113536452A (zh) | 2021-04-22 | 2021-04-22 | 一种适用于多材料多厚度结构的优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113536452A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114329788A (zh) * | 2021-12-31 | 2022-04-12 | 武汉理工大学 | 一种基于田口法和熵权灰色关联分析的车门优化设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104809306A (zh) * | 2015-05-13 | 2015-07-29 | 梧州学院 | 基于多质量指标的注塑成型工艺参数优化方法 |
CN110489903A (zh) * | 2019-08-26 | 2019-11-22 | 四川大学 | 基于可拓学-灰色关联理想解法的床身结构优化设计方法 |
-
2021
- 2021-04-22 CN CN202110437799.0A patent/CN113536452A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104809306A (zh) * | 2015-05-13 | 2015-07-29 | 梧州学院 | 基于多质量指标的注塑成型工艺参数优化方法 |
CN110489903A (zh) * | 2019-08-26 | 2019-11-22 | 四川大学 | 基于可拓学-灰色关联理想解法的床身结构优化设计方法 |
Non-Patent Citations (1)
Title |
---|
SUO ZHANG等: "A two-stage hybrid optimization for honeycomb-type cellular structures under out-of-plane dynamic impact", 《APPLIED MATHEMATICAL MODELLING》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114329788A (zh) * | 2021-12-31 | 2022-04-12 | 武汉理工大学 | 一种基于田口法和熵权灰色关联分析的车门优化设计方法 |
CN114329788B (zh) * | 2021-12-31 | 2024-04-12 | 武汉理工大学 | 一种基于田口法和熵权灰色关联分析的车门优化设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xiong et al. | Lightweight optimization of the side structure of automobile body using combined grey relational and principal component analysis | |
CN107391891B (zh) | 一种基于模型融合方法的大展弦比机翼优化设计方法 | |
CN106650016A (zh) | 基于粒子群算法的车身侧围结构多工况协同优化实现方法 | |
CN113190917A (zh) | 一种多胞吸能结构组合优化设计方法 | |
CN111651929A (zh) | 一种基于Dynaform与智能算法融合的多目标优化方法 | |
CN113536452A (zh) | 一种适用于多材料多厚度结构的优化设计方法 | |
CN114329788B (zh) | 一种基于田口法和熵权灰色关联分析的车门优化设计方法 | |
CN112231839A (zh) | 针对高耗时约束的飞行器追峰采样智能探索方法 | |
CN114235330A (zh) | 关联风洞试验和计算数据的多源气动载荷模型构建方法 | |
CN113836823A (zh) | 一种基于负荷分解和优化双向长短期记忆网络的负荷组合预测方法 | |
Booth et al. | Multi-objective optimization of a multi-cellular aluminum extruded crush rail subjected to dynamic axial and oblique impact loading conditions | |
CN110020466B (zh) | 基于代理模型的负泊松比结构吸能盒协同优化设计方法 | |
CN114595510A (zh) | 一种空天飞行器融合结构外形参数的综合建模方法 | |
Honda et al. | Construction of design guidelines for optimal automotive frame shape based on statistical approach and mechanical analysis | |
Zhang et al. | Multiobjective optimization design for lightweight and crash safety of body-in-white based on entropy weighted grey relational analysis and MNSGA-II | |
CN114072807B (zh) | 基于小样本迁移学习的铣削机器人多模态频响预测方法 | |
CN113158514A (zh) | 汽车车身材料结构匹配轻量化设计方法、系统及存储介质 | |
Izanloo et al. | Optimal design of automotive body B-pillar using simplified finite element model of body-in-prime combined with an optimization procedure | |
Mohamed et al. | Generalized Gaussian mixture model | |
Aulig et al. | Preference-based topology optimization of body-in-white structures for crash and static loads | |
Korta et al. | Multi-objective optimization of a car body structure | |
CN113312827A (zh) | 一种汽车骨架多目标优化方法 | |
CN115099093A (zh) | 一种基于熵权topsis的白车身结构多目标优化设计方法 | |
CN115130303A (zh) | 一种机理与数据相结合的胎噪不确定性优化方法 | |
El-Beltagy et al. | Gaussian processes for model fusion |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211022 |