基于时间演化响应面模型的JWL状态方程参数标定方法

doi:10.6040/j.issn.1671-9352.0.2023.061

摘要/Abstract

摘要：

综合试验数据和建模与模拟结果, 标定爆轰产物JWL(Jones-Wilkins-Lee)状态方程中的唯象参数, 对理解爆轰膨胀驱动过程和评价做功能力至关重要。时间演化代理模型能构建重要参数与系统响应量(system response quantity, SRQ)之间较为简洁的函数关系, 为描述爆轰动力行为和参数标定提供一种快速高精度算法。本文首先利用拉丁超立方抽样(Latin hypercube sampling, LHS)技术对爆轰产物JWL状态方程中的不确定参数进行抽样, 选取有效数据并分成两组: 第一组进行高分辨率仿真, 使用具有自主知识产权爆轰数值软件CFD2D计算得到径向壁位置和径向壁速度, 输出结果用于构建不确定JWL参数和SRQ之间的时间演化响应面模型; 第二组用于评估时间演化代理模型的精确度。其次, 定义代理模型与试验数据之间的残差, 进而搜索最小残差获得最优参数。结果表明: 将最优参数代入高分辨率模型后导出的SRQ的输出结果和试验数据吻合较好, 从而完成了参数标定和模型确认。结果为武器型号设计和装备性能评估提供依据。

关键词: 参数标定, 拉丁超立方抽样, 时间演化响应面模型法, 圆筒试验, JWL状态方程

Abstract:

Calibration of phenomenological parameters of Jones-Wilkins-Lee(JWL) equation of state based on experimental data and modeling and simulation result plays a vital role in comprehension of detonation-driven expansion process and evaluating the work performance of detonation product. Time varying response surface methodology (RSM) can establish a simpler functional relation between prominent parameters and system response quantity (SRQ). It provides a rapid and high-accurate computational method of describing the dynamical behavior and calibrating parameters. Firstly, Latin hypercube sampling (LHS) technique is used to sample the uncertainty in the uncertain parameters of JWL equation of state in the detonation production, and then effective data is selected and divided into two groups. High-resolution simulation is implemented based on group Ⅰ, and the radial position and radial velocity is obtained through full intellectual software—CFD2D. Time varying RSM between uncertain JWL parameters and SRQ is constructed through output of the high fidelity simulation. Group Ⅱ is used to evaluate the accuracy of the time varying RSM. Then Residual error between the simulation results of surrogate model and experiment data is constructed. Furthermore, the minimum residual error is searched among all the residual error and the corresponding parameters are obtained as the optimum parameters. The result shows that output of SRQ coincides with and experimental data perfectly, when these optimum parameters are substituted into the high fidelity computational methods. That means, the parameters are calibrated and model validated. The result will provide guidance for weapon design and munitions assessment.

Key words: parameters calibration, Latin hypercube sampling, time varying response surface methodology, cylinder test, JWL equation of state

中图分类号:

O385

梁霄,王朔,王瑞利,陈江涛. 基于时间演化响应面模型的JWL状态方程参数标定方法[J]. 《山东大学学报(理学版)》, 2024, 59(4): 117-126.

Xiao LIANG,Shuo WANG,Ruili WANG,Jiangtao CHEN. Parameters calibration method of JWL equation of state based on time varying response surface methodology[J]. JOURNAL OF SHANDONG UNIVERSITY(NATURAL SCIENCE), 2024, 59(4): 117-126.

图/表 11

图1

图2

图3

图4

图5

图6

图7

图8

图9

表1

图10

参考文献 21

1	王瑞利, 江松. 多物理耦合非线性偏微分方程与数值解不确定度量化数学方法[J]. 中国科学: 数学, 2015, 45 (6): 723- 738.
	WANG Ruili , JIANG Song . Mathematical methods for uncertainty quantification of the nonlinear partial differential equation and numerical solution[J]. Scientia Sinica Mathematica, 2015, 45 (6): 723- 738.
2	梁霄, 王瑞利. 爆轰敏感性分析与验证CFD模型[J]. 物理学报, 2017, 109, 114- 121.
	LIANG Xiao , WANG Ruili . Sensitivity analysis and validation of detonation CFD model[J]. Acta Physica Sinica, 2017, 109, 114- 121.
3	梁霄, 王瑞利. 基于自适应和投影Wiener混沌的圆筒试验不确定度量化[J]. 爆炸与冲击, 2019, 39 (4): 041408.
	LIANG Xiao , WANG Ruili . Uncertainty quantification of cylindrical test through Wiener chaos with basis adaptation and projection[J]. Explosion and Shock Waves, 2019, 39 (4): 041408.
4	FICKETT W , DAVIS D . Detonation Theory and Experiment[M]. Boston: Dover Publications, 1979.
5	MADER C , FOREST C . Two dimensional homogeneous and heterogeneous detonation wave propagation[J]. Chemical Explosive, 1976, 21, 1- 17.
6	王艳莉, 张树道, 周海兵, 等. 基于不确定度量化分析的炸药驱动飞片参数标定和预测方法[J]. 数值计算与计算机应用, 2017, 38 (2): 143- 154.
	WANG Yanli , ZHANG Shudao , ZHOU Haibing , et al. A uncertainty parameter calibration and prediction method with an application to flyer plate experiments[J]. Journal on Numerical Methods and Computer Applications, 2017, 38 (2): 143- 154.
7	王言金, 张树道, 李华, 等. 炸药爆轰产物Jones-Wilkins-Lee状态方程不确定参数[J]. 物理学报, 2016, 65 (10): 106401.
	WANG Yanjin , ZHANG Shudao , LI Hua , et al. Uncertain parameters of Jones-Wilkin-Lee equation of state for detonation products of explosive[J]. Acta Physica Sinica, 2016, 65 (10): 106401.
8	MORTENSEN C , SOUERS P . Optimizing code calibration of the JWL explosive equation-of-state to the cylinder test[J]. Propellants Explosive Pyrotechnics, 2017, 42 (2): 616- 622.
9	ANDREWS S , FRASER A . Estimating physics models and quantifying their uncertainty using optimization with a Bayesian objective function[J]. ASME Journal of Verification, Validation, and Uncertainty Quantification, 2019, 4 (1): 011002.
10	GARNO J , OUELLET F , BAE S , et al. Calibration of reactive burn and Jones-Wilkins-Lee parameters for simulations of a detonation-driven low experiment with uncertainty quantification[J]. Physical Review Fluids, 2020, 5 (12): 123201.
11	LIANG Xiao , WANG Ruili , GHANEM R . Uncertainty quantification of detonation through adapted polynomial chaos[J]. International Journal for Uncertainty Quantification, 2020, 10 (1): 83- 100.
12	梁霄, 陈江涛, 王瑞利. 高维参数不确定爆轰的不确定度量化[J]. 兵工学报, 2020, 4 (4): 723- 734.
	LIANG Xiao , CHEN Jiangtao , WANG Ruili . Uncertainty quantification of detonation with high-dimensional parameter uncertainty[J]. Acta Armamentarii, 2020, 4 (4): 723- 734.
13	OBERKAMPF W , ROY C . Verification and validation in scientific computing[M]. Boston: Cambridge University Press, 2010.
14	JIANG Xiong , HU Xingzhi , LIU Gang , et al. A generalized active subspace for dimension reduction in mixed aleatory-epistemic uncertainty quantification[J]. Computer Methods in Applied Mechanical Engineering, 2020, 370, 113240.
15	温丽晶, 段卓平, 张震宇, 等. 采用遗传算法确定炸药爆轰产物JWL状态方程参数[J]. 爆炸与冲击, 2013, 33 (1): 130- 134.
	WEN Lijing , DUAN Zhuoping , ZHANG Zhenyu , et al. Determination of JWL-EOS parameters for explosive detonation products using genetic algorithm[J]. Explosion and Shock Waves, 2013, 33 (1): 130- 134.
16	TAGUCHI G . Introduction to quality engineering: white plains[M]. New York: Unipub Press, 1987.
17	BREIG S , LUTI K . Response surface methodology: a review on its applications and challenges in microbial cultures[J]. Material Today: Proceedings, 2021, 42 (5): 2277- 2284.
18	WANG Ruili , LIANG Xiao , LIN Zhong . LAD2D: a computer code for detonation fluid dynamics problems with large distortion and internal slip[J]. Detection and Diffusion Forum, 2017, 379 (3): 58- 63.
19	WANG Ruili , LIANG Xiao , LIU Xuezhe . Research on verification and validation strategy of detonation fluid dynamics code of LAD2D[J]. AIP Conference Proceedings, 2017, 1863 (1): 030007.
20	董海山, 周芬芬. 高能炸药及相关物性能[M]. 北京: 科学出版社, 1989.
	DONG Haishan , ZHOU Fenfen . High explosive and its physical property[M]. Beijing: Science Press, 1989.
21	卢校军, 王蓉, 黄毅民, 等. 两种含铝炸药作功能力与JWL状态方程研究[J]. 含能材料, 2005, 13 (3): 144- 147.
	LU Jiaojun , WANG Rong , HUANG Yimin , et al. Study on work ability and JWL equation of state of two aluminized explosives[J]. Chinese Journal of Energetic Materials, 2005, 13 (3): 144- 147.

多维度评价

Viewed

Full text

Abstract

Cited

Shared

Discussed

系统响应量	E_rf	R₁	R₂	ω
径向壁位置	3.36×10^-5	5.164 77	1.622 09	0.405 80
径向壁速度	1.27×10^-3	10.239 66	2.622 20	0.428 11