Imprecise Probabilistic Model Updating Using A Wasserstein Distance-based Uncertainty Quantification Metric
-
摘要: 复杂物理系统的数学代理模型往往包含多类不确定性因素。在实际工程问题如机械系统可靠性优化设计中,可结合系统响应的部分实测数据,校准模型关键参数取值,修正模型结构,提高代理模型的保真性。但对于具有混合不确定性的非精确概率模型,传统基于欧式距离的模型修正方法并不适用。针对这一问题,提出一种基于Wasserstein距离测度的模型修正方法,该方法基于Wasserstein距离测度构建核函数,利用p维参数空间中Wasserstein距离的几何性质以量化不同概率分布之间的差异,相较于现有模型修正方法,可校准模型的高阶超参数,显著降低模型结构及参数不确定性。针对工程实际需求,进一步采用近似贝叶斯推理与切片分割技术以降低计算成本。通过受迫振动钢板本构参数校核问题与NASA Langley多学科不确定性量化问题验证了本方法在静力学与动力学等实际工程问题中的有效性。
-
关键词:
- Wasserstein距离 /
- 贝叶斯方法 /
- 非精确概率 /
- 不确定性量化 /
- 近似推理
Abstract: Uncertainty factors are usually contained in the mathematical proxy model of complex physical system. In practical engineering problems such as mechanical system reliability optimization design, the key parameters of the model can be calibrated and the model structure can be modified to improve the fidelity of the proxy model. However, for imprecise probabilistic models with mixed uncertainties, the traditional model updating method based on the Euclidean distance is not applicable. To solve this problem, a new model updating method based on the Wasserstein distance measure is proposed, which builds the kernel function based on the Wasserstein distance measure, and uses the geometric properties of Wasserstein distance in P-dimensional parameter space to quantify the differences between different probability distributions. Compared with the existing model updating methods, high-order hyper-parameters of the model can be calibrated to significantly reduce the uncertainty of model structure and parameters. In order to reduce the calculation cost, the approximate Bayesian inference and sliced segmentation technology is further adopted to meet the engineering requirements. The validity of this method for practical engineering problems, such as statics and dynamics, is verified by the constitutive parameter checking problem of forced vibration steel plate and the multidisciplinary uncertainty quantification problem of NASA Langley.-
Key words:
- Wasserstein distance /
- Bayesian methods /
- imprecise probability /
- model updating /
- approximate reasoning
-
表 1 因素水平表
尺寸参数 长/mm 宽/mm 高/mm 数值 600 120 3 本构参数 剪切模量/GPa 弹性模量/GPa 质量密度/(kg·m−3) 数值 83 210 7 860 表 2 认知不确定性建模
参数 不确定性模型 认知不确定性建模 E 弹性模量$ E\sim N\left({\mu }_{E}, {\sigma }_{E}^{2}\right) $ $ {{{\theta }} }_{1}={\mu }_{E}\sim U\left[180, 230\right] $ 均值${\mu _E} \in \left[ {180, 230} \right]$ $ {{{\theta }} }_{2}={\sigma }_{E}^{2}\sim U\left[0, 10\right] $ 方差$\sigma _E^2 \in \left[ {0, 10} \right]$ G 剪切模量$ G\sim N\left({\mu }_{G}, {\sigma }_{E}^{2}\right) $ $ {{{\theta }} }_{3}={\mu }_{G}\sim U\left[60, 100\right] $ 均值${\mu _G} \in \left[ {60, 100} \right]$ $ {{{\theta }} }_{4}={\sigma }_{G}^{2}\sim U\left[0, 10\right] $ 方差$\sigma _G^2 \in \left[ {0, 10} \right]$ 表 3 设定认知不确定性参数的名义值
不确定性参数 名义值 $ E\sim N\left({\mu }_{E}, {\sigma }_{E}^{2}\right) $ ${\mu _{_E}} = 200.5, \sigma _E^2 = 2.0$ $ G\sim N\left({\mu }_{G}, {\sigma }_{G}^{2}\right) $ ${\mu _G} = 81.4, \sigma _G^2 = 1.5$ 表 4 后验结果与真值比较数据
参数 真值/GPa W距离极大似然估计/GPa W距离修正95%置信区间 E距离极大似然估计/GPa E距离修正95%置信区间 ${\mu _E}$ 200.5 200.261 1 [200.257 5, 200.264 7] 200.302 2 [200.298 3, 200.306 0] $\sigma _E^2$ 2.0 2.078 1 [2.064 5, 2.091 8] — — ${\mu _G}$ 81.4 81.567 7 [81.565 2, 81.570 1] 81.553 6 [81.550 5, 81.556 8] $\sigma _G^2$ 1.5 1.408 8 [1.402 8, 1.414 9] — — 表 5 认知不确定性建模
参数 类别 不确定性模型 认知不确定性建模 ${p_1}$ Ⅲ 单峰Beta分布,${\mu _1} \in \left[ {0.6, 0.8} \right], \sigma _1^2 \in \left[ {0.02, 0.04} \right]$ ${{{\theta }} _1} = {\mu _1}, {{{\theta }} _2} = \sigma _1^2$ ${p_2}$ Ⅱ 常数,${p_2} \in \left[ {0.0, 1.0} \right]$ ${{{\theta }} _3} = {p_2}$ ${p_3}$ Ⅰ 均匀分布,$ {p}_{3}\sim U\left(0.0, 1.0\right) $ - ${p_4}, {p_5}$ Ⅲ 二元正态分布,${\mu _i} \in \left[ { - 5.0, 5.0} \right], \sigma _i^2 \in \left[ {0.002\;5, 4.0} \right], \rho \in \left[ { - 1.0, 1.0} \right], i = 4, 5$ ${{{\theta }} _4} = {\mu _4}, {{{\theta }} _5} = {\mu _5}, {{{\theta }} _6} = \sigma _4^2, {{{\theta }} _7} = \sigma _5^2, {{{\theta }} _8} = \rho $ 表 6 不同近似贝叶斯推理结果比较
修正参数 95%置信区间 极大似然估计 真值 分位数 E距离 W距离 分位数 E距离 W距离 ${{{\theta }} _1}$ [0.702 2, 0.709 1] [0.700 6, 0.702 9] [0.683 4, 0.684 4] 0.705 7 0.701 8 0.683 9 0.636 4 ${{{\theta }} _2}$ [0.027 5, 0.028 2] [0.029 9, 0.030 6] [0.031 9, 0.032 3] 0.027 9 0.030 2 0.032 1 0.035 6 ${{{\theta }} _3}$ [0.431 5, 0.474 4] [0.451 1, 0.485 2] [0.297 0, 0.330 8] 0.452 9 0.468 1 0.313 9 1 ${{{\theta }} _4}$ [0.601 5, 0.953 4] [−0.177 8, 0.189 2] [−3.703 6, −3.547 8] 0.777 5 0.005 7 −3.625 7 4 ${{{\theta }} _5}$ [−0.760 7, −0.457 2] [1.685 1, 1.986 3] [−3.771 4, −3.621 4] −0.609 0 1.835 7 −3.696 4 −1.5 ${{{\theta }} _6}$ [1.779 7, 1.923 7] [1.772 8, 1.914 4] [1.819 6, 1.946 5] 1.851 7 1.843 6 1.883 0 0.04 ${{{\theta }} _7}$ [1.959 8, 2.115 1] [1.777 4, 1.921 3] [1.110 9, 1.253 1] 2.037 5 1.849 3 1.182 0 0.36 ${{{\theta }} _8}$ [−0.095 8, −0.026 5] [−0.009 1, 0.059 8] [0.128 6, 0.195 6] −0.061 1 0.025 4 0.162 1 0.5 -
[1] 王龙华. 数学模型校准问题研究[D]. 长沙: 国防科学技术大学, 2012.WANG Longhua. Research on the calibration of mathematical model[D]. Changsha: National University of Defense Technology, 2012. [2] 朱战霞, 张红文, 罗秋月. 自由漂浮空间机器人神经网络自适应控制律设计[J]. 西北工业大学学报, 2017, 35(1): 59-65. https://www.cnki.com.cn/Article/CJFDTOTAL-XBGD201701010.htmZHU Zhanxia, ZHANG Hongwen, LUO Qiuyue. Design of neural network adaptive control law for free-floating space robot[J]. Journal of Northwestern Polytechnical University, 2017, 35(1): 59-65. https://www.cnki.com.cn/Article/CJFDTOTAL-XBGD201701010.htm [3] 喻天翔, 赵庆岩, 尚柏林, 等. 考虑间隙不确定性的花键概率疲劳寿命预测方法[J/OL]. 机械工程学报, 1-12[2022-04-05]. http://kns.cnki.net/kcms/detail/11.2187.TH.20210514.1138.002.htmlYU Tianxiang, ZHAO Qingyan, SHANG Bolin, et al. Probabilistic fatigue life prediction method of spline considering clearance uncertainty[J]. Journal of Mechanical Engineering, 1-12[2022-04-05]. http://kns.cnki.net/kcms/detail/11.2187.TH.20210514.1138.002.html [4] 张小强. 随机与认知不确定性下机械系统可靠性分析与优化设计方法研究[D]. 成都: 电子科技大学, 2018.ZHANG Xiaoqiang. Research on reliability analysis and optimization design method of mechanical system under stochastic and cognitive uncertainty[D]. Chengdu: University of Electronic Science and Technology of China, 2018. [5] 李彦锋, 黄洪钟, 黄意贤. 太阳翼驱动机构的故障模式影响分析与时变可靠性研究[J]. 机械工程学报, 2020, 56(5): 108-115. doi: 10.3901/JME.2020.05.108LI Yanfeng, HUANG Hongzhong, HUANG Yixian. Analysis of failure mode influence and time-varying reliability of solar wing driving mechanism[J]. Journal of Mechanical Engineering, 2020, 56(5): 108-115. doi: 10.3901/JME.2020.05.108 [6] KHODAPARAST H H, MOTTERSHEAD J E, BADCOCK K J. Interval model updating with irreducible uncertainty using the Kriging predictor[J]. Mechanical Systems and Signal Processing, 2011, 25(4): 1204-1226. doi: 10.1016/j.ymssp.2010.10.009 [7] ABU HUSAIN N, KHODAPARAST H H, OUYANG H. Parameter selection and stochastic model updating using perturbation methods with parameter weighting matrix assignment[J]. Mechanical Systems and Signal Processing, 2012, 32: 135-152. doi: 10.1016/j.ymssp.2012.04.001 [8] 彭卫文, 黄洪钟, 李彦锋, 等. 基于数据融合的加工中心功能铣头贝叶斯可靠性评估[J]. 机械工程学报, 2014, 50(6): 185-191. http://www.cjmenet.com.cn/CN/Y2014/V50/I6/185PENG Weiwen, HUANG Hongzhong, LI Yanfeng, et al. Bayesian reliability evaluation of machining center functional milling head based on data fusion[J]. Journal of Mechanical Engineering, 2014, 50(6): 185-191. http://www.cjmenet.com.cn/CN/Y2014/V50/I6/185 [9] BI S, PRABHU S, COGAN S, et al. Uncertainty quantification Metrics with varying statistical information in model calibration and validation[J]. AIAA Journal, 2017, 55(10): 3570-3583. doi: 10.2514/1.J055733 [10] YANG L C, BI S F, FAES M G R, et al. Bayesian inversion for imprecise probabilistic models using a novel entropy-based uncertainty quantification metric[J]. Mech. Syst. Signal Proc., 2022, 162: 21. [11] BI S, BROGGI M, BEER M. The role of the Bhattacharyya distance in stochastic model updating[J]. Mechanical Systems and Signal Processing, 2019, 117: 437-452. doi: 10.1016/j.ymssp.2018.08.017 [12] 张建军, 乔松珊. 正态总体下参数的优化极大似然估计方法[J]. 统计与决策, 2012(2): 16-18. https://www.cnki.com.cn/Article/CJFDTOTAL-TJJC201202005.htmZHANG Jianjun, QIAO Songshan. Optimized maximum likelihood estimation method for parameters in normal population[J]. Statistics and Decision, 2012(2): 16-18. https://www.cnki.com.cn/Article/CJFDTOTAL-TJJC201202005.htm [13] MOTTERSHEAD J E, FRISWELL M I. Model updating in structural dynamics: A survey[J]. Journal of Sound and Vibration, 1993, 167(2): 347-375. doi: 10.1006/jsvi.1993.1340 [14] ROY C J, OBERKAMPF W L. A comprehensive framework for verification, validation, and uncertainty quantification in scientific computing[J]. Computer Methods in Applied Mechanics and Engineering, 2011, 200(25-28): 2131-2144. doi: 10.1016/j.cma.2011.03.016 [15] RENS V D S, SARAH D, RUTH K, et al. Bayesian statistics and modelling[J]. Nature Reviews Methods Primers, 2021, 1(1): 3. https://doi.org/10.1038/ s43586-020-00003-0 doi: 10.1038/s43586-020-00003-0 [16] CRANMER K, BREHMER J, LOUPPE G. The frontier of simulation-based inference[J]. Proceedings of the National Academy of Sciences of the United States of America, 2020, 117(48): 30055-30062. doi: 10.1073/pnas.1912789117 [17] BARAGATTI M, PUDLO P. An overview on approximate Bayesian computation[J]. ESAIM: Proceedings, 2014, 44: 291-299. doi: 10.1051/proc/201444018 [18] CHING J, CHEN Y C. Transitional Markov Chain Monte Carlo method for Bayesian model updating, model class selection, and model averaging[J]. Journal of Engineering Mechanics, 2007, 133: 816-832. [19] RUBNER Y, TOMASI C, GUIBAS L J. The Earth mover's distance as a metric for image retrieval[J]. International Journal of Computer Vision, 2000, 40(2): 99-121. [20] PANARETOS V M, ZEMEL Y. Statistical aspects of wasserstein distances[J]. Annual Review of Statistics and Its Application, 2019, 6: 405-431. [21] NADJAHI K, DE BORTOLI V, DURMUS A, et al. Approximate Bayesian computation with the sliced-wasserstein distance[C]// ICASSP 2020-2020 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), May 1, 2020, Barcelona. SPAIN: IEEE, 2020: 5470-5474. [22] FANG S E, ZHANG Q H, REN W X. An interval model updating strategy using interval response surface models[J]. Mechanical Systems and Signal Processing, 2015, 60-61: 909-927. [23] CRESPO L G, KENNY S P, GIESY D P. The NASA langley multidisciplinary uncertainty quantification challenge[C]//AIAA Non-Deterministic Approaches conference, January 13-17, 2014, Maryland. National Harbor: AIAA, 2014: 13-47. [24] SAFTA C, SARGSYAN K, NAJM H N, et al. Probabilistic methods for sensitivity analysis and calibration in the NASA challenge problem[J]. Journal of Aerospace Information Systems, 2015, 12(1): 219-234. [25] SRIVASTAVA A, SUBRAMANIYAN A K, WANG L. Hybrid Bayesian solution to NASA langley research center multidisciplinary uncertainty quantification challenge[J]. Journal of Aerospace Information Systems, 2015, 12(1): 114-139.