基于改进混合密度网络的毁伤效应预测方法.pdf

第 56 卷第 1 期郑 州 大 学 学 报(理 学 版)Vol.56 No.12024 年 1 月J.Zhengzhou Univ.(Nat.Sci.Ed.)Jan.2024收稿日期:2022-10-11基金项目:河南省重点研发与推广专项(212102310039,202102310554);河南省高等学校重点科研项目(20A520035)。第一作者:佘维(1977),男,教授,主要从事复杂系统建模与仿真、机器学习、区块链与数据智能研究,E-mail:wshe 。通信作者:孔德锋(1989),男,工程师,主要从事工程防护与评估研究,E-mail:kdf 。基于改进混合密度网络的毁伤效应预测方法佘维1,2,张人中1,2,田钊1,2,刘炜1,2,孔德锋3(1.郑州大学 网络空间安全学院河南 郑州 450002;2.郑州市区块链与数据智能重点实验室河南 郑州 450002;3.军事科学院国防工程研究院工程防护研究所河南 洛阳 471023)摘要:提出一种基于改进混合密度神经网络的毁伤效应预测方法,解决了现有智能毁伤效应预测方法中仅能输出点预测结果,但难以量化毁伤效应预测结果的不确定性问题。采用鲁棒性更好的 t 分布作为混合分量,利用混合密度网络生成概率密度函数,以反映毁伤效应预测中的不确定性,并根据给定置信水平获得区间预测结果。仿真实验表明,获得的概率密度函数可以较为准确地拟合蒙特卡洛仿真模拟结果,与现有的毁伤效应预测方法相比,可以更好地指导作战筹划。关键词:混合密度网络;毁伤效应预测;t Location-Scale 分布;区间预测中图分类号:TP399文献标志码:A文章编号:1671-6841(2024)01-0009-07DOI:10.13705/j.issn.1671-6841.2022299Prediction Method of Damage Effects Based on Improved Mixture Density NetworkSHE Wei1,2,ZHANG Renzhong1,2,TIAN Zhao1,2,LIU Wei1,2,KONG Defeng3(1.School of Cyber Science and Engineering,Zhengzhou University,Zhengzhou 450002,China;2.Zhengzhou Key Laboratory of Blockchain and Data Intelligence,Zhengzhou 450002,China;3.Institute of Engineering Protection,Institute of Defense Engineering,Academy of Military Sciences,Luoyang 471023,China)Abstract:A damage effect prediction method based on an improved mixture density neural network was proposed.The method aimed to solve the problem that current intelligent damage effect prediction meth-ods could only output point prediction results,but fail to quantify the uncertainty of damage effect predic-tion results.A more robust t distribution as a mixture component and a mixture density network was used to generate a probability density function that reflects the uncertainty in the prediction of damage effects.Interval prediction results could be obtained according to a given confidence level.Simulation experi-ments demonstrated that the probability density function obtained by the proposed method could more ac-curately fit Monte Carlo simulation results and better guide combat planning than existing damage effect prediction methodsKey words:mixture density network;damage effect prediction;t Location-Scale distribution;interval prediction0引言武器毁伤效应是现代军事和作战研究的基础性问题。常用的毁伤效应预测方法包括理论计算法、计算机仿真法、试验法1-3。近年来,支持向量机4、神经网络5等机器学习方法由于计算效率较高、环境适应性较强,已经在郑 州 大 学 学 报(理 学 版)第 56 卷毁伤效应预测领域中得到初步应用6-9。袁辉等提出了基于最小二乘支持向量机的坑道工程动荷段的毁伤仿真实验训练样本约简模型,其利用粒子群算法选取较优的参数,模型具有较好的可行性和分类精度 10。李建光等针对弹体对混凝土材料侵彻深度问题,通过径向基函数神经网络建立了弹体侵彻深度与网络输入量之间的非线性映射关系11。张磊等针对毁伤实验数据少、不均匀、不连续、范围窄等带来的计算精度不高的问题,运用 K-近邻算法、BP 神经网络建立基于数据融合的“三阶段”毁伤效应计算模型进行毁伤效应计算12。上述基于机器学习与数据挖掘的毁伤效应预测方法存在以下两个问题。1)缺乏不确定性量化能力,不能描述预测值可能的上下限以应对使用过程中存在的决策风险。2)使用单值目标函数难以适应毁伤效应预测结果存在的多峰分布情况。针对上述两个问题,本文提出一种基于改进混合密度网络的毁伤效应预测方法。本文的主要贡献如下。1)通过改进混合密度网络生成的基于混合分布的概率密度函数可以反映毁伤效应预测结果,并能很好地适应毁伤效应预测中存在的多峰分布情况。2)改进的混合密度网络可以通过所得概率密度函数量化预测结果的不确定性。经过处理后,既可以得到点预测结果,也可以根据给定置信水平得到相应置信区间。3)本文提出的改进混合密度网络采用鲁棒性更好的 t 分布作为混合分量,降低了离群点和异常点对模型性能的影响。1相关知识1.1混合分布混合分布是多个不同统计特性的分布函数(混合分量)的凸组合,从 而达到拟合 复杂分布的 效果13-14。混合分布概率密度函数表示为f(x)=nk=1kfk(xk),(1)其中:n 为混合分量个数;k为第 k 个混合分量权重,且满足nk=1k=1,k 0;fk(xk)为第 k 个混合分量 的 概 率 密 度 函 数,k为 第 k 个 混 合 分 量参数。1.2混合密度神经网络混合密度神经网络是混合分布模型与前馈神经网络的结合,其中前馈神经网络的输出被用来确定混合模型的参数。对于给定输入样本 x,输出目标值为 y 的条件概率密度函数表示 y 分布上的不确定性,记为 f(yx)15-16。在 混 合 密 度 网 络(mixture density network,MDN)中,通常采用高斯分布作为混合分量,输出的混合高斯分布表示为f(yx)=nk=1k(x)Nk(yk(x),2k(x),(2)其中:Nk(yk(x),2k(x)为 第 k 个 混 合 分 量,k(x)、2k(x)分 别 为 第 k 个 高 斯 分 布 的 均 值 和方差。然而对于厚尾分布或存在噪声的数据集,基于高斯混合分布的混合密度网络存在鲁棒性较差的情况。原始混合密度网络输入是高维的向量,包含两个隐含层,输出是具有多个分量的混合高斯分布。对于输出层的不同类型参数须采用不同的处理方式,即k(x)=exp(hk)ni=1exp(hi),(3)k(x)=hk,(4)2k(x)=exp(hk)。(5)混合密度神经网络使用最大似然法构造损失函数,MDN 的损失函数定义为L(yx)=-log(nk=1k(x)Nk(yk(x),2k(x)。(6)1.3T Location-Scale 分布T Location-Scale 分布是含有尺度参数和位置参数的 t 分布,t Location-Scale 分布概率密度函数表示为f(x)=(v+1)/2)v(v/2)1+1v(x-)/)2-v+12,(7)其中:、2分别为位置参数和尺度参数所对应一般高斯分布的均值和方差;v 为 t 分布的自由度;为gamma 函数。如图 1 所示,给出了方差为 1,均值为 0,不同自由度下 t Location-Scale 分布的不同形状。可以看到t 分布与高斯分布相似,随着自由度 v 的增大,分布形态逐渐接近高斯分布,其极限分布为高斯分布,当v 120 时,可近似为高斯分布处理。鉴于混合高斯分布和混合 t 分布在处理异常值时表现出的鲁棒性差异,混合 t 分布常常作为混合01第 1 期佘维,等:基于改进混合密度网络的毁伤效应预测方法图 1不同自由度下 t 分布形态Figure 1The shape of t distribution under different degree of freedom高斯分布的替代选型而被使用17。2基于 TDMDN 的毁伤效应预测方法本节提出一种改进的 t 分布混合密度网络(t distribution mixture density network,TDMDN),并基于 TDMDN 提出一种毁伤效应预测方法,首先对效应数据库中的数据进行预处理,再将处理后的数据输入 TDMDN 得到混合分布参数,最后根据混合分布参数生成概率分布函数,得到点预测和区间预测结果。其框架如图 2 所示。图 2毁伤效应预测方法框架Figure 2Damage effects prediction framework2.1TDMDN 网络TDMDN 在混合分量类型和网络结构两个方面对 MDN 进行了改进,并采用极大似然法构造损失函数。在 混 合 分 量 类 型 的 选 择 上,TDMDN 使 用t Location-Scale 分布作为混合分量来解决高斯混合密度网络中存在的鲁棒性问题。对于有 n 个输出混合分量的 TDMDN,有 4n 个输出节点,可以分为、2、v 四个部分,每个部分有 n 个节点,其中 为各个 混 合 分 量 的 权 重,、2、v 分 别 对 应 各 个t Location-Scale 分布的分布参数。对于、2的处理方式与原始 MDN 相同,t Location-Scale 分布中的自由度 v 采用vk(x)=exp(hvk)。(8)在网络结构方面,TDMDN 是由隐含层块组成的深度结构,每个隐含层块由多个全连接层、批归一化层、激活层组成。其中:全连接层是神经网络中最常见的层;批处理归一化层是用来减少初始化的影响,加速网络训练;激活层是为模型提供非线性计算,提高网络非线性拟合能力。TDMDN 的结构细节如图 3 所示。图 3TDMDN 网络结构Figure 3TDMDN network structureTDMDN 使用最大似然法构造损失函数,损失函数定义为L(yx)=-log(nk=1k(x)Tk(yk(x),2k(x),vk)。(9)2.2基于 TDMDN 的毁伤效应预测基于 TDMDN 的毁伤效应预测由以下几个步骤组成。1)数据预处理。首先对效应数据库中的数据进行异常数据处理和归一化处理,得到标准数据。2)TDMDN 训练。将处理得到的标准数据输入TDMDN 进行训练,采用正向传播和反向传播得到训练好的 TDMDN 网络。3)毁伤效应预测。将处理得到的标准数据输入训练好的 TDMDN,得到毁伤效应分布所对应混合分布参数。根据混合分布参数生成毁伤效应分布函数,经过处理得到点预测和区间预测结果。2.2.1数据预处理数据预处理首先是对异常数据处理,接着进行数据归一化处理。训练数据集中如果存在异常数据,不仅会降低11郑 州 大 学 学 报(理 学 版)第 56 卷预测成功率,甚至会导致预测结果与真实值的偏差较大。为了保证在输入模型前数据的规范性和正确性,采用拉伊达准则(3 准则)来对效应数据库中的样本集进行异常值检验和剔除。对训练数据进行异常值检验后,在构建毁伤预测模型之前,对数据集进行归一化处理,以加快模型的收敛,本文采用线性函数归一化方法对数据进行归一化处理。所用公式为y=x-xminxmax-xmin,(10)其中:x 为原始数据;xmin和 xmax分别为原始数据中的最小值和最大值;y 为归一化后的数值。2.2.2TDMDN 模型训练TDMDN 模型的训练过程由正向传播与反向传播两部分组成:正向传播时,样本以向量形式从输入层传入,经过隐含层、批归一化层、激活层的计算,从输出层输出;在反向传播时,通过公式(9)计算误差梯度,然后沿着梯度下降的方向逐层返回,修改每一层神经元的权值和偏置值。当误差减少到可以接受的程度或训练达到指定的次数时,训练停止。由于在模型训练过程中,常常出现梯度爆炸问题,本文针对梯度爆炸对损失函数做了一定改进。在文献18中,总结了 MDN 出现梯度爆炸问题的常见原因,并指出梯度爆炸问题主要来自两个方面。1)最大似然估计值接近 0 时,损失计算在取对数时会出现非常接近于零的值,导致梯度爆炸问题。2)自由度过大时,损失计算会有一个较大的指数,导致梯度爆炸问题。当梯度爆炸损失出现时,整个训练过程就会失败。针对上述原因,在 TDMDN 的训练过程中,提出以下解决方案:采用一个大于 0 的小浮点数与对数相加解决情况 1);采用截断法解决情况 2),当自由度大于 120 时,将自由度截断,令其等于 120,避免损失计算中出现极大指数导致梯度爆炸问题。改进后的损失函数为L(yx)=-log(nk=1k(x)fk(x)+),(11)其中:为大于 0 的小浮点数;fk(x)为混合分量的似然,其定义为fk(x)=Tk(k(x),2k(x),min120,vk)。(12)2.2.3毁伤效应预测在 TDMDN 模型训练完成后,采用训练好的 TDMDN 模型进行毁伤效应预测。给定弹药参数、目标参数以及弹目交汇参数,使用TDMDN 模型得到毁伤效应结果的概率分布情况。根据概率分布情况得到点预测结果和区间预测结果。点预测结果可由混合分量期望值给出,y=nk=1k(x)k(x)。(13)由于混合密度网络所得概率分布是一种不规则分布,很难直接根据不同置信水平得到区间预测结果。为了获得更加可靠有效的区间,我们按如下方式构造预测区间。1)考虑到混合权重过小的混合分量的生成区间不具参考性,故在生成置信区间时选取混合权重k(x)(1/n)的混合分量,n 为混合分量个数。2)对每个混合分量按照置信水平取双侧置信区间。3)将由 2)得到的置信区间取并集得到预测区间。3仿真实验及分析3.1实验设置本文利用文献19-20 提出的工程毁伤算法,仿真实验模拟“某型号动能穿甲弹”,采用五点瞄准法对机枪堡目标的破坏过程获得仿真数据集,数据如表 1 所示。本文只展示部分实验数据,剩余数据用省略号表示。该数据集包含工程长度、工程宽度、工程高度、钢板厚度、墙体厚度、覆土厚度、顶盖厚度、弹药数量和入射速度 9 个输入特征,1 个输出特征为震塌比例。表 1毁伤仿真实验训练样本Table 1Damage simulation experiment training samples高度/m 宽度/m长度/m钢板厚度/cm墙体厚度/cm覆土厚度/cm顶盖厚度/cm弹药数量/个入射速度m/s震塌比例/(%)6.819.4712.570.061.420.490.17821 46310.754.156.876.850.050.800.680.59791 25633.899.219.2414.670.040.870.490.19131 1331.977.8111.9411.510.011.920.400.28511 2242.8721第 1 期佘维,等:基于改进混合密度网络的毁伤效应预测方法如图 4 所示,为一组确定参数,利用仿真模型进行多次仿真,得到目标震塌的概率分布直方图。可以看到由于目标各部位材质、厚度等物理性质的不同,对于多瞄准点的打击,其结果呈现多峰分布的情况。图 4效应结果多峰分布Figure 4Multi-peak distribution of effect resultsTDMDN 由 1 个输入层、3 个隐含层、4 个批处理归一化层和 1 个输出层组成,输入层节点数为 9,隐藏层节点数为 84,使用了自适应学习率优化算法AdamW,初始学习率为 0.001。3.2评价指标对于点预测结果,我们采用平均绝对百分比误差(MAPE)、均方根误差(RMSE)、平均绝对误 差(MAE)进行评价。定义为MAPE=1nni=1y i-yiyi,RMSE=1nni=1y i-yi()2,MAE=1nni=1y i-yi,其中:n 为样本个数;yi、y i分别为第 i 个样本的真实值和预测值。对于概率预测结果,需要评估获得的预测区间。首先对区间可靠性进行评估,预测区间覆盖概率(PICP)被广泛应用于评价区间可靠性,定义为PICP=1nni=1ci,其中:当真实值落入预测区间,即 yi Li,Ui 时,ci=1,否则 ci=0;Li、Ui分别为预测区间的上界和下界。狭窄的预测区间往往比宽大的预测区间更有价值,我们采用归一化平均宽度指标(PINAW)评价区间质量,PINAW 越小,表示预测区间越窄,预测性能越好。定义为PINAW=1nni=1Ui-LiD,其中:D 是基础目标范围上、下界之间的差值,在本文中为震塌比例的上、下界之间的差值。为了综合考虑预测区间的覆盖率和区间宽度,采用平均区间锐度(average interval sharpness,AIS)来评价区间整体质量,区间锐度越大,生成的区间质量越好。第 i 个预测区间的区间锐度 S(xi)定义为S(xi)=-2a(Ui-Li)-4(Li-pi),xi Li,-2a(Ui-Li),Li xi Ui,-2a(Ui-Li)-4(xi-Ui),xi Ui,其中:a 是置信度;AIS 的定义为AIS=1nni=1S(xi)。3.3点预测实验结果分析该部分实验选择决策树、随机森林、支持向量回归、BP 神经网络模型及 MDN 模型与本文 TDMDN模型进行点预测对比实验。实验分别对各模型进行了 MAPE、RMSE、MAE 评估,结果如表 2 所示,黑体数据为最优结果。由表 2 可知 TDMDN 采用鲁棒性更好的 t 分布作为混合分量,降低了离群点和异常点对模型性能的影响,故性能优于 MDN 模型。相较于点预测模型中表现最好 BP 神经网络模型,TD-MDN 的 MAPE 降低了 0.67%,RMSE 降低了 0.02,MAE 降低了 0.01,这是由于 TDMDN 使用多个隐含层块组成的深度结构,具有比 BP 神经网络、决策树、随机森林更好的函数逼近和密度估计能力。表 2点预测结果Table 2Point prediction result方法MAPERMSEMAE决策树12.802.451.43随机森林10.822.111.21支持向量回归27.754.743.14BP 神经网络9.581.811.03MDN9.681.831.05TDMDN8.911.791.023.4区间预测结果分析该部分实验选择基于随机森林的分位数回归模型(QRF)、均值方差估计模型(MVE)、MDN 模型与所提 TDMDN 模型在统一置信度下进行对比实验,置信度设置为 95%。实验结果如表 3 所示,黑体数据为最优结果。QRF 采用分位数得到区间预测结果,而不是根据具体的分布情况,故其效果不如基于混合密度网络的方法。而 MVE 由于仅采用单个高斯模型,不符合真实分布情况,在训练过程中会试图通过增大方差来拟合多峰分布,故 MVE 的 PICP 仅31郑 州 大 学 学 报(理 学 版)第 56 卷比 TDMDN 提高了 0.003,但其 PINAW 和 AIS 远不如其他方法。由于异常值干扰,TDMDN 在各个指标上都好于采用高斯核的 MDN 网络。表 3区间预测结果Table 3Interval prediction result方法PICPPINAWAISQRF0.9640.061-0.121MVE0.9750.096-0.187MDN0.9680.061-0.118TDMDN0.9720.055-0.1063.5蒙特卡洛模拟结果对比我们设定一组固定参数,利用仿真模型进行蒙特卡洛模拟,得到目标震塌的概率分布直方图,即输出震塌比例的真实分布情况。对比 TDMDN、MDN、MVE 三个概率模型的输出概率分布对仿真模型真实分布的拟合情况。如图 5 所示,直方图为仿真模型进行 500 次模拟得到的结果,MVE 由于仅采用单个高斯模型,模型很难拟合真实分布情况,效果最差,由于异常值干扰,与 TDMDN 相比采用高斯核的MDN 网络效果较差。TDMDN 对于异常值有很好的鲁棒性,由于采用混合分布模型,所以能较好拟合真实分布情况,效果最好。图 5分布拟合情况Figure 5Distribution fitting results4结论本文针对传统毁伤效应预测模型缺乏量化不确定性及高斯混合模型存在的鲁棒性较差问题,提出了基于改进混合密度神经网络的毁伤效应预测方法,实验表明相对于传统数据挖掘方法,所提出方法更符合毁伤评估的实际需要,其生成概率密度曲线基本拟合仿真模型多次蒙特卡洛模拟结果。与传统点模型的单值输出相比,概率预测的结果具有一定的容错率,可以更好地指导作战筹划。因此,未来有可能将概率预测有效地应用于决策问题,将各种决策活动转化为概率决策。参考文献:1武健,李亚雄,刘新学,等.目标毁伤效果预测研究J.火力与指挥控制,2017,42(9):84-87,92.WU J,LI Y X,LIU X X,et al.Research summary on technology of damage effect predictionJ.Fire control&command control,2017,42(9):84-87,92.2王力超,乔勇军,李永胜,等.目标毁伤评估方法研究J.舰船电子工程,2020,40(5):116-120.WANG L C,QIAO Y J,LI Y S,et al.Research of bat-tle damage assessment methodJ.Ship electronic engi-neering,2020,40(5):116-120.3李峰,石全,张芳,等.破片和冲击波复合作用下装甲板毁 伤 效 应 预 测 J.火 力 与 指 挥 控 制,2020,45(11):101-105.LI F,SHI Q,ZHANG F,et al.Prediction of damage effect of armor plate under combined action of fragments and shock wave J.Fire control&command control,2020,45(11):101-105.4CHERKASSKY V,MA Y Q.Practical selection of SVM parameters and noise estimation for SVM regressionJ.Neural networks,2004,17(1):113-126.5张宗腾,张琳,谢春燕,等.基于改进 GA-BP 神经网络的目标毁伤效果评估J.火力与指挥控制,2021,46(11):43-48.ZHANG Z T,ZHANG L,XIE C Y,et al.Battle damage effect assessment based on improved GA-BP neural net-work J.Fire control&command control,2021,46(11):43-48.6张悦,孙惠香,冯盛辉,等.基于 IGWO-LSSVM 的弹体侵彻地下洞室毁伤效应预测J.空军工程大学学报(自然科学版),2017,18(6):95-100.ZHANG Y,SUN H X,FENG S H,et al.Effect forecast of projectile penetrating underground arched structure based on IGWO-LSSVM J.Journal of air force engi-neering university(natural science edition),2017,18(6):95-100.7王海宽,李文生,熊飞.基于神经网络的长杆弹侵彻能力预测模型J.计算机仿真,2015,32(2):1-5,66.WANG H K,LI W S,XIONG F.Prediction model on long-rod projectile penetrating into steel target based on neural-networkJ.Computer simulation,2015,32(2):1-5,66.8RYAN S,THALER S.Artificial neural networks for characterizing Whipple shield performanceJ.Procedia 41第 1 期佘维,等:基于改进混合密度网络的毁伤效应预测方法engineering,2013,58:31-38.9RYAN S,THALER S,KANDANAARACHCHI S.Ma-chine learning methods for predicting the outcome of hy-pervelocity impact eventsJ.Expert systems with appli-cations,2016,45:23-39.10 袁辉,王凤山,许继恒,等.基于改进 SVM 的坑道毁伤仿真训练样本约简模型J.解放军理工大学学报(自然科学版),2014,15(2):152-157.YUAN H,WANG F S,XU J H,et al.Training sample reduction model of damage simulation of protective engi-neering based on improved SVMJ.Journal of PLA uni-versity of science and technology(natural science edi-tion),2014,15(2):152-157.11 李建光,李永池,王玉岚.人工神经网络在弹体侵彻混凝土深度中的应用 J.中国工程科学,2007,9(8):77-81.LI J G,LI Y C,WANG Y L.Penetration depth of pro-jectiles into concrete using artificial neural networkJ.Engineering sciences,2007,9(8):77-81.12 张磊,吴昊,赵强,等.基于数据挖掘技术的地下工程目标毁伤效应计算方法J.爆炸与冲击,2021,41(3):4-13.ZHANG L,WU H,ZHAO Q,et al.Calculation method of damage effects of underground engineering objectives based on data mining technology J.Explosion and shock waves,2021,41(3):4-13.13 杨茂,董骏城.基于混合分布模型的风电功率波动特性研究J.中国电机工程学报,2016,36(S1):69-78.YANG M,DONG J C.Study on characteristics of wind power fluctuation based on mixed distribution modelJ.Proceedings of the csee,2016,36(S1):69-78.14 AKPINAR S,AKPINAR E K.Estimation of wind energy potential using finite mixture distribution modelsJ.En-ergy conversion and management,2009,50(4):877-884.15 ZHANG H,LIU Y Q,YAN J,et al.Improved deep mix-ture density network for regional wind power probabilistic forecasting J.IEEE transactions on power systems,2020,35(4):2549-2560.16 MEN Z X,YEE E,LIEN F S,et al.Short-term wind speed and power forecasting using an ensemble of mixture density neural networks J.Renewable energy,2016,87:203-211.17 MCLACHLAN G J,LEE S X,RATHNAYAKE S I.Fi-nite mixture modelsJ.Annual review of statistics and its application,2019,6:355-378.18 GUILLAUMES A B.Mixture density networks for distri-bution and uncertainty estimation D.Barcelona:Uni-versitat Politcnica de Catalunya,2017.19 刘光昆,刘瑞朝,汪维,等.地下钢筋混凝土拱形结构在顶爆条件下的抗爆试验J.含能材料,2021,29(2):157-165.LIU G K,LIU R C,WANG W,et al.Blast resistance experiment of underground reinforced concrete arch struc-ture under top explosionJ.Chinese journal of energetic materials,2021,29(2):157-165.20 汪维,刘光昆,赵强,等.近爆作用下方形板表面爆炸载荷分布函数研究J.中国科学:物理学 力学 天文学,2020,50(2):144-152.WANG W,LIU G K,ZHAO Q,et al.Study on load dis-tributing function of square slab surface under close-in blast loadingJ.Scientia sinica(physica,mechanica&astronomica),2020,50(2):144-152.51