基于活跃子空间的爆压不确定度传递分析.pdfVIP

第39卷第2期高压物理学报Vol.39,No.2

2025年2月CHINESEJOURNALOFHIGHPRESSUREPHYSICSFeb.,2025

DOI:10.11858/gywlxb

基于活跃子空间的爆压不确定度传递分析

1123

梁霄，范孟君，王言金，王瑞利

（1.山东科技大学数学与系统科学学院,山东青岛266590；

2.北京应用物理与计算数学研究所,北京100094；

3.河南省科学院,河南郑州450052）

摘要：爆压间接标定中的不确定度无法消除，不确定度量化能提高模型的可信度和预测能

力。然而，爆压间接标定函数具有复杂非线性结构耦合多输入变量等特征，使得爆压不确定度

传播研究遇到“维数灾难”等问题。活跃子空间是处理爆压不确定度量化的有效工具。首先，导

出系统响应量（systemresponsequantity,SRQ）的梯度协方差矩阵；然后，基于MonteCarlo

方法，寻找活跃变量，即SRQ变化最快的方向；接着，将高维输入不确定度转化成一维空间处理，

避免了“维数灾难”；最后，建立基于一维活跃变量的四阶多项式响应面模型。结果表明，活跃子

空间方法成功刻画了输入不确定度对SRQ的影响，且试验结果落在代理模型预测值的置信区间

内，确认了爆压模型的预测能力。研究还发现，爆压的离散程度较大，与孙承纬的结论吻合。此

外，建立了一种新的爆压模型。该模型是仿射变换与多项式函数的复合运算，具有形式简洁、光

滑性好、鲁棒能力强、运算速度快的特点，且系统输入量是随机变量而非固定值，多项式拟合系

数不因输入不确定度的变化而改变。该研究方法具备体系性，可以推广到其他类型的炸药爆压

预测。

关键词：活跃子空间；不确定度量化；爆压；响应面模型；降维

中图分类号：O241;O521.3文献标志码：A

爆压是刻画炸药爆轰性能的重要指标，关系到炸药的工业应用能力和爆轰数值模拟的精度和可信

度。然而，爆轰是极度剧烈的瞬态（约100ns）物理化学反应，空间尺度极小（约100μm），爆压甚至高达

100GPa[1–5]。爆轰的极端特征给爆压标定带来巨大挑战。直接标定爆压需要复杂的流程、昂贵的设

备、精密的仪器仪表。与之相比，炸药密度的测量成本较低，爆速标定试验也相对简单。因此，先标定

易于测量的物理量，再使用工程经验和物理定律间接导出爆压，成为爆压标定的常用方法。特别是研

[6–9]

制新配方炸药时，间接标定爆压几乎成为必用的手段。但是，间接标定使用的数学模型具有模型结

构复杂、多系统输入量（multi-input）等特征，且在一定程度上直接或间接使用最小二乘法，从而导致系

统响应量（systemresponsequantity,SRQ）对系统输入量的赋值敏感。爆压标定值伴随着系统输入量的

变化而剧烈抖动。事实上，爆轰领域的专家学者已普遍达成共识，即受炸药组分的异质性和测试仪器

的精度限制，系统输入量的标定值必然会在一定范围内波动[1–2,10]。爆轰的复杂性导致爆压间接测试模

型中引入了数学物理假设、模型简化和近似，这是爆压标定中不确定度的又一来源。

如上所述，爆压间接标定中的输入不确定度是不可避免的[3–4,9,11]。不确定度的存在又在一定程度

上限制了间接标定的精度、效率和可靠性。遗憾的是，目前，绝大多数的爆压标定公式都是建立在确