弹道极限试验值异常性检验方法研究

弹道极限试验值异常性检验方法研究赵晓旭1,2，王树山1，徐豫新1*（1.北京理工大学爆炸科学与技术国家重点实验室，北京100081；2.首都师范大学信息工程学院，北京100048）摘要：针对基于量纲分析获得弹道极限计算模型函数基础上通过若干试验值拟合获得计算模型的方法，研究弹体侵彻靶体试验所获弹道极限值异常性的检验方法。根据统计学原理，提出了特定置信水平下标准化残差值在标准正态分布置信区间落入度的弹道极限试验值检验方法；同时，进行了3.0g、4.5g、6.0g、7.5g、9.0g和10.0g共6种质量35CrMnSi钢FSP弹体对4mm和5mm两种厚度典型高强度低合金（HSLA）钢板的侵彻试验获得了12个弹道极限试验值，基于量纲分析及试验数据获得了相应的弹道极限计算模型；并以此模型为实例根据95%置信水平下标准化残差值在（-1.96，1.96）区间落入度进行了弹道极限计算模型的检验和修正，表明所提方法可用于弹道极限试验值异常性的检验。关键词：终点效应；数理统计；弹道极限；残差中图分类号：TJ012.4文献标识码：AResearchonValidationMethodofBallisticLimitExperimentDataAbnormalityZHAOXiao-xu1,2,WANGShu-shan1,XUYu-xin1*(1.StateKeyLaboratoryofExplosionScienceandTechnology,BeijingInstituteofTechnology,Beijing100081,China;2.CollegeofInformationEngineering,CapitalNormalUniversity,Beijing100048,China)Abstract:Basedonthemethodofballisticlimitcalculationmodelobtainedfromthedimensionalanalysis,andcalculationmodelfittedfromexperimentdata,thevalidationmethodofprojectileimpactingtargetplateexperimentdataabnormalityisstudiedinthispaper.Accordingtostatisticstheory,theballisticlimitexperimentdatavalidationmethodaccordingtheconfidenceintervalofstandardizedresidualfallsintostandardnormaldistributionatspecificconfidencelevelisproposed;inaddition,withtheexperimentof3.0g,4.5g,6.0g,7.5g,9.0gand10.0gtotal6kindsofmass35CrMnSisteelFSPprojectileimpacting4mmand5mmthicknesstypicalhighstrengthlowalloy(HSLA)steelplates,total12ballisticlimitexperimentdataisobtained,withdimensionalanalysisandexperimentdata,thecorrespondingballisticlimitcalculationmodelisobtained;andbasedonthestandardizedresidualfallinginto(-1.96,1.96)confidenceintervalat95%confidencelevel,theballisticlimitcalculationmodelisvalidatedandrevised,whichindicatesthattheproposedmethodcanbeusedtovalidatetheabnormalityoftheballisticlimitexperimentdata.Keywords:Ballisticlimit;experimentdata;residual;validationmethod;基金项目：国家自然科学基金项目（11402027）北京理工大学基础科研基金资助项目（20130242005）作者简介：赵晓旭（1977—），女，讲师，博士生，E-mail:[e-mail]通讯作者：徐豫新（1982—），男，讲师，硕士生导师，E-mail:[e-mail]弹道极限速度（V50）是指弹体在特定着角条件下50%概率贯穿一定厚度靶体的入射速度，表现为统计学意义上的弹体贯穿速度，作为材料和结构抗弹体侵彻性能的重要表征参量，在军用和警用装甲防护中得到广泛应用[1,2]。长期以来，关于不同弹靶材料之间弹道极限速度的试验和理论模型研究一直都是研究的热点[3-8]。这些研究中，弹道极限计算模型获取的最常见办法是基于量纲分析后得到相应的计算模型函数，并基于试验获得若干弹道极限值后拟合获得计算模型[9]。拟合获得计算模型计算精度的检验标准通常是用若干个试验点处计算值与试验值相对统计误差的极差值（极差又称全距，是用来表示统计资料中的变异量数，即最大值减最小值后所得之数据）规定的。这在实际应用中，就存在两个问题：①同一极差值下，随着弹道极限的提高，由极差值产生的允许区间将增加，如：规定极差值不超过15%，1000m/s弹道极限速度的允许区间为150m/s，100m/s弹道极限速度的允许区间仅为15m/s，因基数的不同造成在弹道极限较低时过于严格，在弹道极限较高时过于宽松；②在实际的试验中，由于装药量、弹体着角等客观因素的影响，试验结果会存在随机性，且弹道极限本身就是50%穿透概率下的特征撞击速度，试验点可能存在异常，但异常试验点尚无法发现，这就会使拟合获得计算模型的适用度难以明确。本工作基于数理统计的方法提出特定置信水平下标准化残差值在标准正态分布置信区间落入度的弹道极限试验值异常试验点检验方法。为了检验该方法的合理性，通过试验获得了3.0g、4.5g、6.0g、7.5g、9.0g和10.0g共6种质量35CrMnSi钢FSP弹体对4mm和5mm两种厚度典型高强度低合金（HSLA）钢板的弹道极限速度，基于量纲分析和试验数据获得了弹道极限计算模型，采用所提方法检验出了试验异常点，并进行了计算模型的修正。1检验方法在数理统计中，残差是指实际观察值与估计值（拟合值）之间的差值。在回归分析中，测定值与按回归方程预测的值之差，即为残差，残差遵从正态分布。显然，有多少对数据，就有多少个残差，即可根据多个残差获取残差的数字特征。定义标准化残差为残差与残差均值之差与残差标准差的比值，即式(1)所示。Siiˆ*(1)式中，*i为第i个标准化残差值；i为第i个残差；ˆ为N个残差的均值；S为N个残差的标准差。“残差”蕴含了有关模型基本假设的重要信息，对于通过回归分析获得的数学模型，实质是用连续曲线近似地刻画或比拟平面上离散点组，以表示坐标之间函数关系的一种数据处理方法。残差分析就是利用残差表征数据点与它在回归直线上相应位置的差异信息，来考察模型假设的合理性及数据的可靠性。因此，可以根据试验数据和基于试验数据拟合获得的模型，通过计算获得残差值和标准化残差值。因残差服从正态分布，因此标准化残差值服从标准正态分布，可根据标准正态分布表获得标准化残差值落在（-A，A）区间以外的概率为p≤a。那么，根据统计学原理，若某一变量的标准化残差值落在（-A，A）区间以外，则可在(1-a)100%置信度将其判为异常试验点，不参与回归拟合。若A=1.645时，a=0.1；A=1.96时，a=0.05；A=2.24时，a=0.025。2方法应用实例为了进行上述方法的适用性验证，设计了长径比在1.5~1.68之间的3.0g、4.5g、6.0g、7.5g、9.0g和10.0g共6种质量35CrMnSi钢质FSP弹体对4mm和5mm两种厚度典型高强度低合金（HSLA）钢板（屈服强度：498.75MPa，拉伸强度：585MPa）的垂直（即着角为0度）侵彻试验，按常用6射弹法获得了12个弹道极限的试验值，列于表1中。表1弹体侵彻钢板试验数据Tab.1Experimentdataofprojectilepenetratingsteelplate弹体质量/g靶板厚度/mm3发贯穿最低速度/(m/s)3发嵌入最高速度/(m/s)平均速度/(m/s)标准差/(m/s)速度1速度2速度3速度4速度5速度63.04757.3733.5775.5727.6735.5730.4743.319.05826.0820.2815.1780.9792.3783.6803.019.84.54672.7637.0665.3631.4629.1625.8643.620.25723.5716.4732.0702.5710.3691.8712.814.56.04582.3592.6586.9554.9550.6558.4571.018.35706.8709.0699.2675.8679.7689.6693.413.97.54561.6557.8553.7522.8520.2534.5541.818.35671.6666.7653.1654.3628.0637.0651.816.89.04520.1516.3513.0481.1485.0480.2499.319.05637.3638.0640.0630.0612.2603.8626.915.210.04494.6486.6473.8461.8456.4452.5471.017.05617.6591.5578.3586.8585.4596.5592.713.6假设FSP直径为pd，长度为pL，弹速为pv，靶体运动速度为tv，靶体厚度为tH，弹体的着角为φ，上述符号的下标p是弹的标志，符号的下标t是靶体的标志；此外，假设弹体和靶体材料的惯性和强度由密度、弹性常数和屈服极限代表，分别为p，pE，pY和t，tE，tY。因为弹体和靶体均为钢质，弹体速度远低于材料声速，可认为材料的可压缩性并不重要[10-14]，根据量纲分析可得：pttptttpttptpppptPYvLHYEYYYELdfYv，，，；，；；(2)根据上述具体试验问题，弹、靶材料确定，靶体速度为0m/s，弹体着角为0度，则式(2)可简化为：ptppLHLdfv；50(3)由弹体尺寸可知，弹体长径比相差不大，则长径比可认为常数，根据试验结果可以获得弹道极限随弹靶相对厚度的变化关系。那么，根据试验获得的弹道极限速度采用最小二乘法进行拟合，获得弹道极限速度与弹靶相对厚度的变化关系曲线，如图1所示。拟合得到的弹道极限速度计算模型为式(4)。4.221507.20.8069tptpHLHLv，　50(4)根据试验值和计算模型式(4)，采用上述方法对弹道极限试验值异常性进行检验，具体如下：将试验值、式(4)计算所得的弹道极限、相对统计误差以及标准化残差值列于表2中。2.02.22.42.62.83.03.23.43.63.84.04.24.4450500550600650700750800850弹道极限速度/(m/s)弹靶相对厚度实验点拟合曲线图1弹道极限随弹靶相对厚度的变化关系Fig.1Ballisticlimitvelocityvarieswithprojectile-targetrelativethickness表2式(4)计算值与试验值对比Tab.2Comparisonofexperimentdataandcalculationdatafromequation(4)序号弹体质量/g靶体厚度/mm