圆筒试验JWL状态方程参数的贝叶斯标定

陈华; 周海兵; 刘国昭; 孙占峰; 张树道

doi:10.11883/1001-1455(2017)04-0585-06

圆筒试验JWL状态方程参数的贝叶斯标定

doi: 10.11883/1001-1455(2017)04-0585-06

1.
北京应用物理与计算数学研究所，北京 100094
2.
中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室，四川绵阳 621999

基金项目:

国家自然科学基金项目 11472060

国家自然科学基金项目 11671413

中国工程物理研究院科学技术发展基金项目 ZYSZ1518-13

北京应用物理与计算数学研究所所长青年基金项目 ZYSZ1518-13

详细信息

作者简介:
陈华(1979-)，女，博士，副研究员，chen_hua@iapcm.ac.cn

中图分类号: O381
计量
- 文章访问数: 4275
- HTML全文浏览量: 1286
- PDF下载量: 404
- 被引次数: 0
出版历程
- 收稿日期: 2015-12-30
- 修回日期: 2016-04-25
- 刊出日期: 2017-07-25

Bayesian calibration for parameters of JWL equation of state in cylinder test

1.
Institute of Applied Physics and Computational Mathematics, Beijing 100094, China
2.
Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621999, Sichuan, China

摘要

摘要: 在研究数值模拟的输入参数引入的不确定性时，通常需要人为给定每个输入参数的概率分布，且输入参数概率分布的选择可能会对分析结果产生直接影响。利用贝叶斯方法标定了圆筒试验JWL状态方程参数，得到了标定参数的估计值和后验分布，并研究了不同统计模型假设对标定参数的估计值和后验分布的影响。贝叶斯后验分布融合了基准试验的试验数据的信息，因此将其作为不确定度量化分析时输入参数的初始概率分布，可以尽量减少分布选择引入的认知不确定性。
- 不确定度量化 /
- 参数标定 /
- 后验概率分布 /
- 贝叶斯方法 /
- JWL状态方程 /
- 圆筒试验 /
- 高斯过程
Abstract: Since probability distributions of input parameters are usually assigned subjectively for uncertainty quantification (UQ) in numerical simulations, the selection of distributions may have significant effect on analysis results of UQ. In this paper, to calibrate more precisely the parameters of JWL equation of state in the cylinder test, we proposed to adopt Bayesian methods that provide estimators and posterior distributions of calibration parameters. Further, we investigated the effects of model assumptions on estimators and posterior distributions of calibration parameters. Our study shows that, owing to the information of cylinder tests they contain, the posterior distributions can be used as initial probability distributions of the input parameters in UQ to reduce the epistemic uncertainty.
- uncertainty quantification /
- calibration /
- posterior distribution /
- Bayesian methods /
- JWL equation of state /
- cylinder test /
- Gaussian process

HTML全文

图 1 圆筒试验的试验装置示意图

Figure 1. Device of cylinder test

下载: 全尺寸图片幻灯片

图 2 圆筒试验径向半径和时间的关系

Figure 2. Radius vs. time in cylinder test

下载: 全尺寸图片幻灯片

图 3 GP代理模型预测的结果和数值计算结果的比较

Figure 3. Comparison of predictions of the GP surrogate model with numerical simulations

下载: 全尺寸图片幻灯片

图 4 两模型的后验概率等高线和最大后验估计

Figure 4. Posterior distribution and MPE with different models

下载: 全尺寸图片幻灯片

图 5 径向半径的预测、试验和模拟结果的比较

Figure 5. Comparison of radius histories obtained by prediction, experiment and simulation

下载: 全尺寸图片幻灯片

参考文献(16)

[1]	Oberkampf W L, Roy C J. Verification and validation in scientific computing[M]. Cambridge: Cambridge University Press, 2010.
[2]	Smith R C. Uncertainty quantification: Theory, implementation, and applications[M]. Philadelphia: Society for Industrial and Applied Mathematics, 2013.
[3]	Kury J W, Honig H C, Lee E L, et al. Metal acceleration by chemical explosives[C]//4th Symposium on Detonation. Maryland: White Oak, 1966: 3-13.
[4]	Lan I F, Hung S C, Chen C Y, et al. An improved simple method of deducing JWL parameters from cylinder expansion test[J]. Propellants, Explosives, Pyrotechnics, 1993, 18(1):18-24. doi: 10.1002/(ISSN)1521-4087
[5]	于川, 刘文翰, 李良忠, 等.钝感炸药圆筒试验与爆轰产物JWL状态方程研究[J].高压物理学报, 1997, 11(3):227-233. http://www.cnki.com.cn/Article/CJFDTOTAL-GYWL703.011.htm Yu Chuan, Liu Wenhan, Li Liangzhong, et al. Studies on cylinder test and JWL equation of state of detonation product for insensitive high explosive[J]. Chinese Journal of High Pressure Physics, 1997, 11(3):227-233.. http://www.cnki.com.cn/Article/CJFDTOTAL-GYWL703.011.htm
[6]	陈朗, 黄毅民, 冯长根.含铝炸药圆筒试验及爆轰产物JWL状态方程研究[J].火炸药学报, 2001, 24(3):13-15. doi: 10.3969/j.issn.1007-7812.2001.03.005 Chen Lang, Huang Yimin, Feng Changgen. The cylinder test and JWL equation of state detontion product of aluminized explosives[J]. Chinese Journal of Explosives and Propellants, 2001, 24(3):13-15. doi: 10.3969/j.issn.1007-7812.2001.03.005
[7]	陈朗, 李泽仁.激光速度干涉仪测量法在炸药圆筒试验中的应用[J].爆炸与冲击, 2001, 21(3):229-232. doi: 10.3321/j.issn:1001-1455.2001.03.012 Chen Lang, Li Zeren. The cylinder tests measured by VISAR interferometer[J]. Explosion and Shock Waves, 2001, 21(3):229-232. doi: 10.3321/j.issn:1001-1455.2001.03.012
[8]	徐辉, 孙占峰, 李庆忠.标准圆筒试验数据处理和不确定度评定方法[J].北京理工大学学报, 2010, 30(5):626-630. http://d.old.wanfangdata.com.cn/Periodical/bjlgdxxb201005027 Xu Hui, Sun Zhanfeng, Li Qingzhong. Study on data processing and uncertainty evaluation of standard cylinder test[J]. Transactions of Beijing Institute of Technology, 2010, 30(5):626-630. http://d.old.wanfangdata.com.cn/Periodical/bjlgdxxb201005027
[9]	Lee E L, Tarver C M. Phenomenological model of shock initiation in heterogeneous explosives[J]. Physics of Fluids, 1980, 23(12):2362-2372. doi: 10.1063/1.862940
[10]	Wescott B L, Stewart D S, Davis W C. Equation of state and reaction rate for condensed-phase explosives[J]. Journal of Applied Physics, 2005, 98(5):053514. doi: 10.1063/1.2035310
[11]	孙承纬, 卫玉章, 周之奎.应用爆轰物理[M].北京:国防工业出版社, 2000.
[12]	Kennedy M, O'Hagan A. Bayesian calibration of computer models (with discussion)[J]. Journal of the Royal Statistical Society Series B: Statistical Methodology, 2001, 68:425-464.
[13]	Williams B, Higdon D, Gattiker J, et al. Combining experimental data and computer simulations, with an application to flyer plate experiments[J]. Bayesian Analysis, 2006, 1(1):765-792. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Open J-Gate000000530154
[14]	Higdon D, Gattiker J, Williams B, et al. Computer model calibration using high-dimensional output[J]. Journal of the American Statistical Association, 2008, 103(482):570-583. doi: 10.1198/016214507000000888
[15]	Bayarri M J, Berger J O, Rui P, et al. A framework for validation of computer models[J]. Technometrics, 2007, 49(2):138-154. doi: 10.1198/004017007000000092
[16]	张树道, 周海兵, 熊俊.相容拉氏流体动力学CHAP程序研制及其应用研究: GF-A0091846G[R]. 2005.