CN104036100B - 不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法 - Google Patents
不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法 Download PDFInfo
- Publication number
- CN104036100B CN104036100B CN201410310063.7A CN201410310063A CN104036100B CN 104036100 B CN104036100 B CN 104036100B CN 201410310063 A CN201410310063 A CN 201410310063A CN 104036100 B CN104036100 B CN 104036100B
- Authority
- CN
- China
- Prior art keywords
- design
- reliability
- bayesian
- optimization
- deviation
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012937 correction Methods 0.000 title claims abstract description 20
- 238000013461 design Methods 0.000 claims abstract description 62
- 238000005457 optimization Methods 0.000 claims abstract description 42
- 230000004044 response Effects 0.000 claims abstract description 32
- 238000004088 simulation Methods 0.000 claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 13
- 238000002474 experimental method Methods 0.000 claims abstract description 5
- 230000008569 process Effects 0.000 claims description 8
- 238000013401 experimental design Methods 0.000 claims description 6
- 238000013213 extrapolation Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 238000013400 design of experiment Methods 0.000 abstract description 14
- 238000005516 engineering process Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 238000010993 response surface methodology Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，属于汽车可靠性设计优化技术领域。该方法包括以下步骤：步骤一：定义基于可靠性设计优化(RBDO)问题；步骤二：为贝叶斯推理偏差模型以及初始响应面模型构建试验设计(DOE)矩阵；步骤三：使用步骤二中所述的偏差模型修正初始响应面模型并量化来自于重复试验和CAE仿真的不确定性；步骤四：运行RBDO优化程序寻最优、最可靠解；步骤五：进行蒙特卡洛仿真(MCS)验证所得解的可靠性。本方法考量了各种不确定性因素对产品性能参数的影响，能有效保证产品整体性能，从而在设计中既保证经济效益又保证实际运行中的安全可靠。
Description
技术领域
本发明属于汽车可靠性设计优化技术领域，涉及一种不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法。
背景技术
基于有限元技术仿真的工程设计与优化是一种解决汽车被动安全设计问题的有效的工具。模型验证即是通过比较CAE模型输出及所测试验结果以评估CAE模型在拟用途上的有效性和准确性的过程。成功的模型验证能明显减少在汽车设计中样机建立和试验的投入。
在实际仿真设计过程中，不考虑设计变量的不确定性将直接影响车身结构设计的性能指标，易导致车身结构指标不能满足实际使用要求，从而使得车身结构丧失在规定的使用条件下、规定的时间内完成规定功能的能力，因此，考虑这种固有不确定性的可靠性设计技术至关重要。
可靠性优化设计(RBDO)方法是结构参数化设计过程中考虑设计变量不确定因素及其对目标和约束函数影响的一种优化设计技术。可靠性优化设计把产品的总体可靠度作为性能约束的优化，产生与合理安全性相协调的平衡设计。考量各种不确定性因素对机械结构零部件性能参数的影响，采用可靠性优化设计方法对车身轻量化参数设计，能有效保证车身结构实际工程应用中的结构整体性能，从而在车身设计中既保证产品的经济效益又保证实际运行中的安全可靠。
响应面或替代模型(Response Surface、Surrogate Model)技术是用来开发、改进、优化的统计和数学方法。用一个完整有限元模型进行直接的设计优化或稳健设计时，仿真仍然会花费大量时间。完成一个有限元模型有时会花费数小时甚至数天。为了减少开发时间，响应面模型(RSM)被广泛用作有限元模型的代理。
无论使用哪种响应面方法，在实际FE模型和RSM之间始终会有偏差，于是需要对随机模型内插和外推方法进行研究，以修正RSM和FE之间的偏差以及提供修正偏差后RSM的预测区间。模型更新技术不但可以量化数据不确定性大小，而且可以进一步提高模型的预测能力。本发明所提供的贝叶斯推理技术可达到响应面模型修正的目的并展现出良好的效果。
发明内容
有鉴于此，本发明的目的在于提供一种不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，该方法考虑已有CAE模型的可控和不可控变量的不确定性，以得到与传统确定性优化解相比具有更高可靠性和稳健性的设计。
为达到上述目的，本发明提供如下技术方案：
一种不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，包括以下步骤：步骤一：定义基于可靠性设计优化(RBDO)问题；步骤二：为贝叶斯推理偏差模型以及初始响应面模型构建试验设计(DOE)矩阵；步骤三：使用步骤二中所述的偏差模型修正初始响应面模型并量化来自于重复试验和CAE仿真的不确定性；步骤四：运行RBDO优化程序寻最优、最可靠解；步骤五：进行蒙特卡洛仿真(MCS)验证所得解的可靠性。
进一步，步骤一中所述的可靠性设计优化(RBDO)是指结构参数化设计过程中考虑设计变量不确定因素及其对目标和约束函数影响的一种优化设计；可靠性优化设计把产品的总体可靠度作为性能约束的优化，产生与合理安全性相协调的平衡设计，考量各种不确定性因素对机械结构零部件性能参数的影响，从而在设计中既保证产品的经济效益又保证实际运行中的安全可靠。
进一步，在步骤二中为贝叶斯推理偏差模型构建试验设计DOE矩阵，通过对DOE矩阵下的变量组合进行试验仿真，所得结果用于响应面模型建立。
进一步，所述步骤三中贝叶斯推理偏差修正响应面模型基于贝叶斯理论，将初始响应面模型预测与CAE仿真之间的差值作为偏差先验信息，通过贝叶斯推理得到后验偏差并建立偏差响应面进行验证域的偏差外推预测，从而达到初始响应面修正的目的。
进一步，所述步骤四中的优化程序区别于传统确定性优化将优化解逼近约束边界，RBDO考虑来自变量的不确定性，所得优化解在一定变动范围内仍能够保证稳健性和可靠性。
进一步，所述步骤五具体包括：通过蒙特卡洛仿真(MCS)所得的FE结果与随机偏差修正后的RSM预测一致，则此过程完成，否则修改DOE矩阵直到达到满意结果为止。
本发明的有益效果在于：本发明所述的汽车可靠性设计优化方法考量了各种不确定性因素对产品性能参数的影响，能有效保证产品整体性能，从而在设计中既保证经济效益又保证实际运行中的安全可靠。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：
图1为本发明所述方法的流程示意图。
具体实施方式
下面将结合附图，对本发明的优选实施例进行详细的描述。
图1为本发明所述方法的流程示意图，如图所示，本发明所述的不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，包括以下步骤：步骤一：定义基于可靠性设计优化(RBDO)问题；步骤二：为贝叶斯推理偏差模型以及初始响应面模型构建试验设计(DOE)矩阵；步骤三：使用步骤二中所述的偏差模型修正初始响应面模型(RSM)并量化来自于重复试验和CAE仿真的不确定性；步骤四：运行RBDO优化程序寻最优、最可靠解；步骤五：进行蒙特卡洛仿真(MCS)验证所得解的可靠性。
具体来说：
在步骤二中为贝叶斯推理偏差模型构建试验设计DOE矩阵，通过对DOE矩阵下的变量组合进行试验仿真，所得结果用于响应面模型建立。
步骤三中贝叶斯推理偏差修正响应面模型基于贝叶斯理论，将初始响应面模型预测与CAE仿真之间的差值作为偏差先验信息，通过贝叶斯推理得到后验偏差并建立偏差响应面进行验证域的偏差外推预测，从而达到初始响应面修正的目的。
步骤四中的优化程序区别于传统确定性优化将优化解逼近约束边界，RBDO考虑来自变量的不确定性，所得优化解在一定变动范围内仍能够保证稳健性和可靠性。
步骤五具体包括：通过蒙特卡洛仿真(MCS)所得的FE结果与随机偏差修正后的RSM预测一致，则此过程完成，否则修改DOE矩阵直到达到满意结果为止。
下面利用本发明所述的不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化流程对某汽车安全系统进行可靠性优化设计。
下面通过具体实施例来说明本方法的具体实施：
在本实施例中，碰撞条件为100％的正面碰撞，汽车以56.6km/h的速度撞上刚性墙面。设计优化目标为在满足汽车安全性能的条件下得到车身轻量化的目的。
本例的RBDO问题公式化如下：
Find：
Minimize：μweight
Subject to：P{CG≤CGT arget}≥99％,CGT arget＝65
P{CD≤CDT arget}≥99％,CDT arget＝750
Weight＝6.012x1+3.166x2+2.078x3+1.237x4+1.463x5+4.369x6+3.547x7+2.306x8
CG＝84.699-7.7668x6+0.7635x7-13.133x1-0.999x2x5+4.40899x1x6-0.3187x4x8+0.2922x1x5
CD＝922.51-2.5605x6x7+0.6625x4-88.269x1+13.929x1 2-1.2664x3x6+0.4711x4x5-8.2049x2x6-4.6859x4x8
其中，Lxi和Uxi是设计变量的上下界，CG和CD作为汽车碰撞安全响应量分别代表假人胸部加速度和车身压缩距离。
为了达到对比的目的，在本实施例中使用了三种方法进行设计优化：传统确定性优化、使用原始低保真度CG和CD响应面模型的传统RBDO以及使用贝叶斯推理修正偏差的模型的RBDO，不仅考虑设计变量的变化，还考虑模型的不确定性。
其中基于贝叶斯推理的响应面修正过程始于在设计空间中用于验证的试验设计(DOE)矩阵，随后反复进行CAE仿真和试验。计算出二者间的差作为偏差分布超参数的贝叶斯推理的依据。得到先验分布后，预测偏差超参数的后验分布得以计算。
试验输出Yt(x)和模型输出Ym(x)之间的关系常概括为：
Yt(x)＝Ym(x)+δ(x)+ε(x)
因此输出的实值Yture，等于试验观察值减去试验误差，同时也等于CAE模型预测加上模型偏差:
Yture(x)＝Yt(x)-ε(x)＝Ym(x)+δ(x)
假设偏差τ(x)＝δ(x)+ε(x)服从一个已知方差的正态分布：
设方差已知，而均值τμ需根据已知数据得出。在这里τμ的先验分布可设为：
τμ0～N(μ0,σ0 2)
在贝叶斯理论中，后验分布等于似然函数与先验分布p(θ)的乘积除以p(x)的概率：
则τμ后验分布可表示为：
τμ1～N(μ1,σ1 2)
其中
表1所示为三种优化方法所得结果：
为了验证所得优化解是否达到了可靠性的目的，根据传统RBDO结果的平均值生成了10个MC样本，并运行高保真度的有限元仿真，得到CG和CD相应的值。这是一个失败的解：10个设计中有1个超出了CG的约束范围内，而10个设计的CD都在约束范围之外。此外，原始多项式的RSM的预测能力有限。由MCS得到CG和CD的概率分布与10个有限元仿真结果偏离过大。其中6个CG的有限元仿真以及10个CD的仿真超出了预测分布，两约束的可靠度都未达到99％。
为了验证随机偏差修正的RBDO的解，另运行10个FE仿真，得到CG和CD相应的值。10个MC有限元仿真所得结果与随机偏差修正的RSM相吻合，且所有的设计都在CG和CD的正态分布中。两约束的可靠度均达到设计目标。由此可见，本发明所述的汽车可靠性设计优化方法能有效保证产品整体性能。
最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。
Claims (4)
1.一种不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，其特征在于：包括以
下步骤：
步骤一：定义基于可靠性设计优化问题；
步骤二：为贝叶斯推理偏差模型以及初始响应面模型构建试验设计矩阵；
步骤三：使用步骤二中所述的偏差模型修正初始响应面模型并量化来自于重复试验和CAE 仿真的不确定性；
步骤四：运行可靠性设计优化程序寻最优、最可靠解；
步骤五：进行蒙特卡洛仿真验证所得解的可靠性；
步骤一中所述的可靠性设计优化是指结构参数化设计过程中考虑设计变量不确定因素及其对目标和约束函数影响的一种优化设计；可靠性优化设计把产品的总体可靠度作为性能约束的优化，产生与合理安全性相协调的平衡设计，考量各种不确定性因素对机械结构零部件性能参数的影响，从而在设计中既保证产品的经济效益又保证实际运行中的安全可靠；
在步骤二中为贝叶斯推理偏差模型构建试验设计矩阵，通过对试验设计矩阵下的变量组合进行试验仿真，所得结果用于响应面模型建立。
2.根据权利要求1 所述的不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，其特征在于：所述步骤三中贝叶斯推理偏差修正初始响应面模型基于贝叶斯理论，将初始响应面模型预测与CAE 仿真之间的差值作为偏差先验信息，通过贝叶斯推理得到后验偏差并建立偏差响应面进行验证域的偏差外推预测，从而达到初始响应面模型修正的目的。
3. 根据权利要求1 所述的不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，其特征在于：所述步骤四中的优化程序区别于传统确定性优化将优化解逼近约束边界，可靠性设计优化考虑来自变量的不确定性，所得优化解在一定变动范围内仍能够保证稳健性和可靠性。
4.根据权利要求1 所述的不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法，其特征在于：所述步骤五具体包括：通过蒙特卡洛仿真所得的FE 结果与随机偏差修正后的响应面模型预测一致，则此过程完成，否则修改试验设计矩阵直到达到满意结果为止。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410310063.7A CN104036100B (zh) | 2014-07-01 | 2014-07-01 | 不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410310063.7A CN104036100B (zh) | 2014-07-01 | 2014-07-01 | 不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104036100A CN104036100A (zh) | 2014-09-10 |
| CN104036100B true CN104036100B (zh) | 2017-05-10 |
Family
ID=51466871
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410310063.7A Expired - Fee Related CN104036100B (zh) | 2014-07-01 | 2014-07-01 | 不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104036100B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10449958B2 (en) | 2017-02-15 | 2019-10-22 | Ford Global Technologies, Llc | Feedback-based control model generation for an autonomous vehicle |
| CN111709082B (zh) * | 2020-04-28 | 2022-06-07 | 湖南大学 | 一种高效的汽车侧面碰撞安全可靠性设计优化方法 |
| CN111753463A (zh) * | 2020-05-22 | 2020-10-09 | 重庆长安汽车股份有限公司 | 车辆行驶跑偏主动控制方法 |
| CN114896693B (zh) * | 2022-05-09 | 2025-05-27 | 浙江吉利控股集团有限公司 | 汽车可靠性设计方法、装置、设备及可读存储介质 |
| CN116186934B (zh) * | 2023-02-28 | 2025-09-30 | 重庆长安汽车股份有限公司 | 一种螺栓结构优化设计方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102663176A (zh) * | 2012-03-28 | 2012-09-12 | 北京航空航天大学 | 针对高可靠机械产品的主动可靠性分析评价方法 |
| CN103294857A (zh) * | 2013-05-24 | 2013-09-11 | 重庆大学 | 以频响位移为优化目标的约束阻尼板的拓扑优化方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7493574B2 (en) * | 2006-02-23 | 2009-02-17 | Cadence Designs Systems, Inc. | Method and system for improving yield of an integrated circuit |
-
2014
- 2014-07-01 CN CN201410310063.7A patent/CN104036100B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102663176A (zh) * | 2012-03-28 | 2012-09-12 | 北京航空航天大学 | 针对高可靠机械产品的主动可靠性分析评价方法 |
| CN103294857A (zh) * | 2013-05-24 | 2013-09-11 | 重庆大学 | 以频响位移为优化目标的约束阻尼板的拓扑优化方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于稳健与可靠性优化设计的轿车车身轻量化研究;张宇;《中国博士学位论文全文数据库<工程科技II辑>》;20130315(第3期);第1-162页 * |
| 结合定性知识的定量贝叶斯模型外推方法;郑凯等;《上海交通大学学报》;20120630;第46卷(第6期);第994-998页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104036100A (zh) | 2014-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104036100B (zh) | 不确定性下基于贝叶斯偏差修正的汽车可靠性设计优化方法 | |
| Guo et al. | A reinforcement learning decision model for online process parameters optimization from offline data in injection molding | |
| Soize | A comprehensive overview of a non-parametric probabilistic approach of model uncertainties for predictive models in structural dynamics | |
| CN103136428B (zh) | 基于两种不确定的轿车车身结构稳健设计方法 | |
| CN103268082B (zh) | 一种基于灰色线性回归的热误差建模方法 | |
| CN103092078A (zh) | 多阶段间歇生产过程的全程优化方法 | |
| CN104573281B (zh) | 一种考虑回弹补偿的复杂空间曲面薄板成型模面设计方法 | |
| CN102270302B (zh) | 一种基于灰色支持向量机的多应力加速寿命试验预测方法 | |
| RU2016132376A (ru) | Моделирование деформации топливного элемента | |
| Liu et al. | Virtual modelling integrated phase field method for dynamic fracture analysis | |
| CN103544335A (zh) | 在制造钣金件的计算机模拟中确定结构破坏的方法和系统 | |
| CN111310318A (zh) | 基于数字孪生的工艺余量处理方法、系统及机械制造总成 | |
| CN104021432A (zh) | 基于改进灰色预测模型的电力负荷中长期预测方法 | |
| CN108109206A (zh) | 锂电池安全阀开启压力有限元建模方法 | |
| CN114297779A (zh) | 一种飞行器气动力系数快速预测方法、系统及设备 | |
| CN103116698A (zh) | 一种基于三次样条的gm(1,1)模型预测方法 | |
| CN118036416B (zh) | 接触疲劳损伤分析方法、装置、设备及存储介质 | |
| Tang et al. | Online updating belief-rule-base using Bayesian estimation | |
| CN109902389B (zh) | 基于改进通用似然估计的不确定性有限元模型修正方法 | |
| CN116415651A (zh) | 安全强化学习和安全控制方法和装置、智能体和存储介质 | |
| CN103530435A (zh) | 一种基于敏感度的船体型线设计方法 | |
| CN111915079B (zh) | 一种混合knn风电功率预测方法及系统 | |
| Zeng et al. | Vehicle crashworthiness performance prediction through fusion of multiple data sources | |
| CN106874561B (zh) | 基于牛顿迭代的多学科不确定性传播分析方法 | |
| CN103678738A (zh) | 基于时域多步积分的互连线模型降阶方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170510 Termination date: 20180701 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |