CN112414605A - 一种动爆场地表掠过压测量装置及威力场构建方法 - Google Patents

Abstract

本发明公开了一种动爆场地表掠过压测量装置，包括基座和沿基座表面依次排布的膜片，沿基座表面依次排布有M个膜片(M≥2)，膜片通过二阶通孔贴附设置于基座上，并与基座表面保持平齐。膜片通过二阶通孔贴附设置于基座上，并与基座表面保持平齐。本发明还公开了一种动爆场地表掠过压测量装置的标定方法和一种利用动爆场地表掠过压测量装置的动爆威力场构建方法。本发明的装置结构简单，质量轻、体积小，加工及运输成本低，可用于大范围内大量布设；标定方法，能够有效获取压力峰值与膜片变形的关系；动爆威力场构建方法，可获取动爆冲击波压力的空间分布特征，进一步通过数学插值获取弹药动爆下不同特征压力值对应的等压线簇。

Description

一种动爆场地表掠过压测量装置及威力场构建方法

技术领域

本发明属毁伤评估技术领域，主要涉及一种冲击波测量装置及威力场构建方法，特别是一种适用于动爆场的地表掠过压测量装置及动爆条件下的威力场构建方法。

背景技术

弹药动爆条件下，炸点具有随机性，若采用电测法测量，一方面传感器及测试线缆损坏率较高，且测量可靠性差；另外，在野外靶场环境，受到地形、环境等因素的限制，电测法的实施难度大；除此以外，为了尽可能测量到动爆条件下弹药的爆炸冲击波压力，需要在较大范围内布设大量冲击波压力传感器，从而导致测试成本过高。因此，要将电测法应用于动爆场景测量冲击波压力有很大的缺陷与不足。

压力效应靶是一种基于金属膜片机械变形量来表征、测量冲击波压力峰值的工具。根据结构毁伤P-I准则，当膜片的一阶固有振动周期小于冲击波正压脉宽的1/10时，可近似认为膜片的变形只与冲击波压力峰值相关，那么就可以通过理论分析及试验建立冲击波压力峰值和膜片残余变形量的量效关系。在实际使用时，只需将效应靶布设于预定的测点位置，试验后通过工具测量效应靶的变形量，就可以获取该测点的冲击波压力峰值。大量试验结果表明，采用效应靶测量弹药爆炸后中、远场的压力峰值有较好的准确性及重复性。

与常用的压力传感器电测方法相比，效应靶虽然无法获取完整的冲击波压力时域波形，但有以下优点：

①操作效率高。在使用时测试人员仅需将靶体固定安装在测点位置即可，无需搭建复杂电测系统，极大缩短测点布设时间。

②测试成本低。目前效应靶靶体加工成本低，每使用一次仅需更换压力敏感膜片即可重复使用。

③适用范围广。能够适应多种复杂地形、气候条件。

基于上述优点，效应靶作为压力电测法的有益补充，在一些压力电测法难以有效发挥作用的场合中，效应靶能起到不可替代的作用。如：

①密集空投试验爆炸冲击波测量。军事演训中弹药常通过搭载武器系统投放，落点随机性较大且投弹次数较为密集。为了获取弹药冲击波压力，需要在弹药预期落点附近布设大量压力测点。若采用压力电测法，一方面会面临测点通道数超负荷且测试成本过高的问题；另外在每波次投弹结束后都需要较长时间更换传感器，调试测试系统，这在军事演训活动中几乎难以实施。而由于效应靶成本较低，能够大量阵列式布设，且检靶换靶效率高，能够满足这种使用工况。

②恶劣条件下爆炸冲击波测量。在一些人员、物资、电力等保障条件不成熟的靶场中，或地形地貌复杂的环境中，有可能不具备电测条件。效应靶无需电源、引线和数采系统，能够用于山地、高原、戈壁以及极端低温环境。

压力效应靶目前已经被大量用于测量弹药静、动爆条件下的冲击波压力，如发明专利“测量可移动靶标冲击波压力峰值的效应靶结构及测试方法”、“一种适用于动爆场靶标冲击波压力测量的效应靶装置”、“一种适用于动爆场冲击波压力测量的模块化的效应靶装置”等所述的效应靶。

然而目前的效应靶装置仍存在一些不足和缺陷为：

①效应靶感压膜片通过压板固定于基座上，感压膜片凹陷，且与靶面不平齐，从而仅适用于测量斜入射压或正入射压，测量掠过压时，冲击波从靶体表面传播至膜片上时，会受到凹陷下坡度的影响而导致冲击波压力测量不准；

②以上效应靶结构复杂，质量较重，从而导致运输成本高、使用人力成本高，在加工成本上仍有很大的降幅空间，在结构设计上也能进一步优化；

③目前对于测量爆炸反射压力的效应靶可以通过激波管来标定(爆炸冲击波参数薄膜测试法研究，姬建荣等)。然而，由于地表掠过压测量装置在测量时冲击波在膜片表面不同位置的作用时序有差异，不能采用激波管来标定。

另外，由于冲击波威力场能够表征弹药威力，并能用于确定弹药的毁伤范围与区域，等压线是描述冲击波威力场的重要方法。然而，目前由于缺乏高效、经济的手段获取弹药动爆条件下大范围内的地表掠过压，依靠目前的冲击波压力测量手段难以在试验中获取大量冲击波压力数据并建立弹药动爆条件的冲击波威力场。

发明内容

本发明要解决的技术问题是，设计一种冲击波压力测量结构及使用该测量结构构建冲击波威力场的方法，该压力测量结构使用成本低，可大量布设，并且能用于测量动爆条件下地表掠过压，该威力场构建方法能通过试验测量建立弹药动爆条件下的冲击波威力场。

为了解决上述技术问题，本发明采用如下技术方案：

一种动爆场地表掠过压测量装置，包括基座和沿基座表面依次排布的膜片，沿基座表面依次排布有M个膜片(M≥2)，膜片通过二阶通孔贴附设置于基座上，并与基座表面保持平齐。

膜片通过二阶通孔贴附设置于基座上，并与基座表面保持平齐。

进一步地，所述的基座为至少有两个相对对称平行面的空心立体结构，且两个相对对称平行面通过堵头实现可拆卸封堵，沿基座轴向上表面依次设置有与膜片配合的二阶通孔，沿二阶通孔排布方向还设置有与二阶通孔间隔设置的测量定位孔；与二阶通孔所在相邻的侧面还设置有引压孔。

进一步地，基座通过固定支架与地表面固定。

本发明还公开了一种动爆场地表掠过压测量装置的标定方法，包括如下步骤：

步骤一：在静爆靶场划定直径不小于20m，且地表空旷、平整、夯实区域为试验区域，以该区域中心位置为爆心，在爆心处上方一定高度位置放置TNT炸药；

步骤二：以爆心为中心划一条压力测线，在测线上阵列布设N个冲击波压力传感器，传感器爆心距为D_i(i＝1,2,3…N)；在距爆心相同位置处，布设N个地表掠过压测量装置；布设过程中，确保冲击波压力传感器、动爆场地表掠过压测量装置上表面与地表平齐；

步骤三：起爆炸药，然后对冲击波压力传感器波形并进行数据处理，依次记录N个压力峰值P_i(i＝1,2,3…N)；

步骤四：依次记录N个动爆场地表掠过压测量装置上各个膜片的变形量，形成变形量矩阵Y_N×M，矩阵第i(i＝1,2,3,…,N)行第j(j＝1,2,3,…,M)列的元素表示第i个地表掠过压测量装置上第j个膜片的变形；对矩阵Y_N×M的每一行求均值，得到N个膜片变形均值y_i(i＝1,2,3,…,N)；

步骤五：根据压力峰值P_i与变形量y_i的数据序列，采用线性拟合方法，拟合出压力峰值关于膜片变形的函数P＝f(y)。

本发明进一步公开了一种利用动爆场地表掠过压测量装置的动爆威力场构建方法，包括如下步骤：

步骤一：预设动爆场地表掠过压测量装置置布设位置

以弹药瞄准点为中心，采用笛卡尔二维建立坐标系，横坐标为x，纵坐标为y；选定L×L平面范围作为冲击波压力测试区域，该区域在x、y方向的坐标范围均为[-L/2,L/2]；动爆场地表掠过压测量装置在该区域内N行N列网状排布，相邻动爆场地表掠过压测量装置的间距为：

△L＝L/(N-1)

建立测压装置布设坐标矩阵C，该矩阵第i行、第j列(i,j＝1,2,3,…,N)的元素表示第i行、第j爆场地表掠过压测量装置的坐标，即布置位置：

C_i,j＝[x_i,y_j]＝[-L/2+△L(i-1),-L/2+△L(j-1)]

步骤二：动爆场地表掠过压测量装置结构安装

将N×N个动爆场地表掠过压测量装置逐一安装在预设布设位置，安装后确保测压装置上表面与测试区域地表平齐；通过细钢钎穿过安装定位孔将各个动爆场地表掠过压测量装置结构固定在地表；

步骤三：试验实施

开展动爆毁伤试验，记录炸点坐标[x₀,y₀]；

步骤四：动爆场地表掠过压测量装置数据测量及处理

逐一测量N×N个测点位置处测压装置膜片q(q＝1,2,3…N)的特征挠度值，并根据q个膜片的平均特征挠度计算该测压装置的压力值，并形成N行N列的压力矩阵P，矩阵P第i行、第j列(i,j＝1,2,3,…,N)的元素表示第i行、第j列动爆场地表掠过压测量装置测量的掠过压：

步骤五：威力场构建

给定威力场特征压力序列P_c＝[P₁,P₂,…,P_n]，通过插值算法计算等压线簇：

Φ(x,y)＝[f₁(x,y)＝C₁,f₂(x,y)＝C₁,…,f_n(x,y)＝C_n]；

其中f_i(x,y)＝C_i上的压力为P_i。

本发明的有益效果体现在以下几个方面：

(1)本发明提出的冲击波压力测量装置结构简单，加工成本、运输成本及操作人力成本低廉，方便在大范围内大量布设；

(2)本发明提出的冲击波压力测量装置改进了传统的凹陷式膜片结构，采用了平齐式膜片结构，克服了该结构测量掠过压时的不足，能用于测量弹药动爆条件下的地表掠过压；

(3)本发明提出的一种地表掠过压测量装置标定方法，能通过试验快速获取地表掠过压峰值与膜片变形的量效关系；

(4)本发明提出了一种动爆场冲击波威力场构建方法，能够基于上述地表掠过压测量装置获取动爆条件下大范围内的冲击波压力值，进而采用插值算法获取不同压力值对应的等压线，进而确定动爆条件下弹药的冲击波毁伤范围。

附图说明

图1是本发明的动爆场地表掠过压测量装置及其测量方法示意图；

图2是膜片凹陷与平齐两种情况下膜片表面冲击波压力曲线仿真对比；

图3为本发明实施例的动爆场地表掠过压测量装置结构示意图；

图中标号为：1-基座；2-膜片安装孔-1；3-膜片安装孔-2；4-膜片-1；5-膜片

-2；6-测量定位孔；7-堵头；8-引压孔；9-L形支架；10-安装定位孔；

图4是本发明实施例1中地表掠过压测量装置标定结果；

图5是本发明实施例2中的地表掠过压测量装置布设位置；

图6是本发明实施例2中的等压线计算结果；

表1是本发明实施例1中的压力峰值与膜片变形测量结果。

下面结合附图及优选实施例对本发明作进一步的详述。

具体实施方式

参见图1一种动爆场地表掠过压测量装置，包括基座和沿基座表面依次排布的膜片，膜片数量至少两个，膜片通过二阶通孔贴附设置于基座上，并与基座表面保持平齐。膜片为软铝材料，直径与二阶通孔的上台阶孔直径相等。

在现有技术中，测量可移动靶标冲击波压力峰值的方法不同，该装置粘贴于车辆等可移动靶标侧面，在弹药触地爆炸后，冲击波通常会近似垂直入射结构表面，从而该装置中膜片表面冲击波同时到达，且冲击波波阵面传播方向与膜片垂直。

由于地表掠过压的传播方向水平，波阵面通常与地表垂直，这样会导致现有技术压力效应靶测量地表掠过压存在一定的误差，误差有两方面原因：其一为该装置膜片相对于结构表面凹陷，冲击波从膜片表面掠过时，冲击波流动方向会发生改变，这会严重影响冲击波压力峰值测量结果。图2中采用AUTODYN仿真对比了炸药爆炸冲击波掠过地表时，膜片有凹陷与无凹陷两种情况下膜片表面的冲击波压力曲线。从图中可见，膜片无凹陷时，压力峰值为3.49MPa，膜片凹陷时，压力峰值为4.03MPa，在相同工况下，表面凹陷压力效应靶的测量结果比表面平齐效应靶的测量结果高15.5％左右。其二，测量地表掠过压时，膜片表面受力有时序性。相对比膜片直径，冲击波波阵面无限薄，因此，冲击波压力先作用于距离爆心更近的地方，再作用于距离爆心较远的地方，这就导致采用激波管来标定膜片变形与压力峰值的方法在这种场景下不能使用。

本发明的动爆场地表掠过压测量装置，基座为至少有两个相对对称平行面的空心立体结构。该基座的主要作用为约束膜片边界，并给膜片变形留出足够深度。基座的两个相对对称平行面通过堵头实现可拆卸封堵，该堵头能有效防止爆炸冲击波从基座端面进入膜片背面，从而影响膜片变形。

沿基座轴向上表面依次设置有与膜片配合的二阶通孔，沿二阶通孔排布方向还设置有与二阶通孔间隔设置的测量定位孔，为沉孔，用于定位膜片变形测量装置；与二阶通孔所在相邻的侧面还设置有引压孔，该引压孔能有效防止膜片变形导致基座内外产生压力差，从而干扰膜片变形。

基座通过固定支架与地表面固定。固定支架为L形、或工字型等结构，其为不锈钢材料，固定支架上还留有安装孔，用于通过螺钉与基座连接。

为了进一步验证本发明的技术效果，本发明将动爆场地表掠过压测量装置选取一个较优实施例结构，具体参数如下：基座长132mm，宽60mm，高30mm，边缘厚度均为4mm，膜片安装孔-1、膜片安装孔-2均为台阶孔，上台阶沉孔直径50mm，深度0.3mm，下台阶通孔直径30mm，膜片-1、膜片-2均为经过热处理的软铝材料，直径50mm，厚度0.2mm；测量定位孔共3个，为直径3mm、深度3mm的沉孔，用于定位效应靶挠度测量装置；堵头为聚乙烯材料；引压孔为直径2mm的通孔，位于靶体的侧面；将膜片-1、膜片-2通过专用胶水分别粘贴于膜片安装孔-1、膜片安装孔-2上，用堵头堵住基座两端，将L形支架通过螺钉安装在基座上，组成地表掠过压测量装置，如图3所示。

下面通过以下几步标定该地表掠过压测量装置的变形与压力的关系：

步骤一：在静爆靶场划定直径不小于20m，且地表空旷、平整、夯实区域为试验区域，以该区域中心位置为爆心，在爆心处上方一定高度位置放置TNT炸药；

步骤二：以爆心为中心划一条压力测线，在测线上阵列布设N个冲击波压力传感器，传感器爆心距为D_i(i＝1,2,3…N)；在距爆心相同位置处，布设N个地表掠过压测量装置；布设过程中，确保冲击波压力传感器、动爆场地表掠过压测量装置上表面与地表平齐；

步骤三：起爆炸药，然后对冲击波压力传感器波形并进行数据处理，依次记录N个压力峰值P_i(i＝1,2,3…N)；

步骤四：依次记录N个动爆场地表掠过压测量装置上各个膜片的变形量，形成变形量矩阵Y_N×M，矩阵第i(i＝1,2,3,…,N)行第j(j＝1,2,3,…,M)列的元素表示第i个地表掠过压测量装置上第j个膜片的变形；对矩阵Y_N×M的每一行求均值，得到N个膜片变形均值y_i(i＝1,2,3,…,N)；

步骤五：根据压力峰值P_i与变形量y_i的数据序列，采用线性拟合方法，拟合出压力峰值关于膜片变形的函数P＝f(y)。

本发明的还公开了利用动爆场地表掠过压测量装置的动爆威力场构建方法，包括如下步骤：

步骤一：预设动爆场地表掠过压测量装置置布设位置

以弹药瞄准点为中心，采用笛卡尔二维建立坐标系，横坐标为x，纵坐标为y；选定L×L平面范围作为冲击波压力测试区域，该区域在x、y方向的坐标范围均为[-L/2,L/2]；动爆场地表掠过压测量装置在该区域内N行N列网状排布，相邻动爆场地表掠过压测量装置的间距为：

△L＝L/(N-1)

建立测压装置布设坐标矩阵C，该矩阵第i行、第j列(i,j＝1,2,3,…,N)的元素表示第i行、第j爆场地表掠过压测量装置的坐标，即布置位置：

C_i,j＝[x_i,y_j]＝[-L/2+△L(i-1),-L/2+△L(j-1)]

步骤二：动爆场地表掠过压测量装置结构安装

将N×N个动爆场地表掠过压测量装置逐一安装在预设布设位置，安装后确保测压装置上表面与测试区域地表平齐；通过细钢钎穿过安装定位孔将各个动爆场地表掠过压测量装置结构固定在地表；

步骤三：试验实施

开展动爆毁伤试验，记录炸点坐标[x₀,y₀]；

步骤四：动爆场地表掠过压测量装置数据测量及处理

逐一测量N×N个测点位置处测压装置上M个膜片的特征挠度值，并根据M个膜片的平均特征挠度计算该测压装置的压力值，并形成N行N列的压力矩阵P，矩阵P第i行、第j列(i,j＝1,2,3,…,N)的元素表示第i行、第j列动爆场地表掠过压测量装置测量的掠过压；

步骤五：威力场构建

给定威力场特征压力序列P_c＝[P₁,P₂,…,P_n]，通过插值算法计算等压线簇：

Φ(x,y)＝[f₁(x,y)＝C₁,f₂(x,y)＝C₁,…,f_n(x,y)＝C_n]；

其中f_i(x,y)＝C_i上的压力为P_i。

实施例1

针对图3所示的地表掠过压测量装置，通过以下方法获取了膜片变形与压力的关系：

在静爆靶场划定直径不小于20m，且地表空旷、平整、夯实区域为试验区域，以该区域中心位置为爆心，开展了2发试验，第一发试验中，在爆心处上方0.5m处放置9kgTNT炸药，距离爆心3.46m、4.97m处放置了压力传感器及地表掠过压测量装置；第二发试验中，在爆心正上方0.5m放置了6kgTNT炸药，距离爆心5.39m、6.42m、3.46m、4.97m和6.01m分别放置了压力传感器和地表掠过压测量装置；在布设过程中，确保压力传感器及地表掠过压测量装置与地表平齐；

经过两发试验，统计数据，见表1。

表1试验数据结果

测点编号	1	2	3	4	5	6	7
								压力峰值(传感器)/MPa	0.1569	0.167	0.2355	0.2359	0.284	0.562	0.863
变形挠度(膜片-1)/mm	1.81	1.84	2.28	2.38	2.83	4.18	6.24
								变形挠度(膜片-2)/mm	1.72	1.92	2.39	2.38	2.84	4.3	6.84
平均变形挠度/mm	1.765	1.88	2.335	2.38	2.835	4.24	6.54

对压力峰值和平均变形挠度进行数学拟合，可得：

P＝0.1521y-0.1198

其中y为平均变形挠度，P为压力峰值，标定结果如图4所示，R平方为0.99711，说明试验结果线性度很好，采用一阶线性拟合满足精度要求。

实施例2

由于真实试验数据涉及保密事项，下面举一个虚拟实例来阐述本专利的实施方式。假设在某次动爆试验中，采用上述地表掠过压测量装置测量大范围的冲击波压力峰值，并构建冲击波威力场，实施过程包括如下步骤：

步骤一：预设测压装置布设位置

以弹药瞄准点为中心笛卡尔二维建立坐标系，横坐标为x，纵坐标为y；选定40m×10m平面范围作为冲击波压力测试区域，该区域在x、y方向的坐标范围均为[-20,20]；效应靶在该区域内11行11列网状排布，相邻效应靶的间距为

△L＝40/(11-1)＝4m

建立效应靶布设位置矩阵C，该矩阵第i行、第j列(i,j＝1,2,3,…,11)的元素表示第i行、第j列效应靶的坐标：

C_i,j＝[x_i,y_j]＝[-20+4(i-1),-20+4(j-1)]

布设位置如图5所示。

步骤二：效应靶结构安装

将121个效应靶逐一安装在预设布设位置，安装后确保效应靶上表面与测试区域地表平齐；通过细钢钎穿过安装定位孔将结构固定在地表；

步骤三：试验实施

开展动爆毁伤试验，记录炸点坐标[5,-7]；

步骤四：动爆场地表掠过压测量装置数据测量及处理

逐一测量N×N个测点位置处测压装置上M个膜片的特征挠度值，并根据M个膜片的平均特征挠度计算该测压装置的压力值，并形成N行N列的压力矩阵P，矩阵P第i行、第j列(i,j＝1,2,3,…,N)的元素表示第i行、第j列动爆场地表掠过压测量装置测量的掠过压；

步骤五：威力场构建

给定威力场特征压力序列P_c＝[0.1,0.5,1](单位：MPa)，通过插值算法计算等压线簇如图6所示。

Claims (6)

1.一种动爆场地表掠过压测量装置，包括基座和沿基座表面依次排布的膜片，其特征在于：沿基座表面依次排布有M个膜片(M≥2)，膜片贴附设置于基座上，并与基座表面保持平齐。

2.如权利要求1所述的动爆场地表掠过压测量装置，其特征在于：所述的膜片通过二阶通孔贴附设置于基座上，并与基座表面保持平齐。

3.如权利要求1或2所述的动爆场地表掠过压测量装置，其特征在于：所述的基座为至少有两个相对对称平行面的空心立体结构，且两个相对对称平行面通过堵头实现可拆卸封堵，沿基座轴向上表面依次设置有与膜片配合的二阶通孔，沿二阶通孔排布方向还设置有与二阶通孔间隔设置的测量定位孔；与二阶通孔所在相邻的侧面还设置有引压孔。

4.如权利要求3所述的动爆场地表掠过压测量装置，其特征在于：基座通过固定支架与地表面固定。

5.一种权利要求1-4任意一项所述的动爆场地表掠过压测量装置的标定方法，其特征在于：包括如下步骤：

步骤一：在静爆靶场划定直径不小于20m，且地表空旷、平整、夯实区域为试验区域，以该区域中心位置为爆心，在爆心处上方一定高度位置放置TNT炸药；

步骤二：以爆心为中心划一条压力测线，在测线上阵列布设N个冲击波压力传感器，传感器爆心距为D_i(i＝1,2,3…N)；在距爆心相同位置处，布设N个地表掠过压测量装置；布设过程中，确保冲击波压力传感器、动爆场地表掠过压测量装置上表面与地表平齐；

步骤三：起爆炸药，然后对冲击波压力传感器波形并进行数据处理，依次记录N个压力峰值P_i(i＝1,2,3…N)；

步骤四：依次记录N个动爆场地表掠过压测量装置上各个膜片的变形量，形成变形量矩阵Y_N×M，矩阵第i(i＝1,2,3,…,N)行第j(j＝1,2,3,…,M)列的元素表示第i个地表掠过压测量装置上第j个膜片的变形；对矩阵Y_N×M的每一行求均值，得到N个膜片变形均值y_i(i＝1,2,3,…,N)；

步骤五：根据压力峰值P_i与变形量y_i的数据序列，采用线性拟合方法，拟合出压力峰值关于膜片变形的函数P＝f(y)。

6.一种利用权利要求1-4任意一项所述的动爆场地表掠过压测量装置的动爆威力场构建方法，其特征在于：包括如下步骤：

步骤一：预设动爆场地表掠过压测量装置置布设位置

以弹药瞄准点为中心，采用笛卡尔二维建立坐标系，横坐标为x，纵坐标为y；选定L×L平面范围作为冲击波压力测试区域，该区域在x、y方向的坐标范围均为[-L/2,L/2]；动爆场地表掠过压测量装置在该区域内N行N列网状排布，相邻动爆场地表掠过压测量装置的间距为：

△L＝L/(N-1)

建立测压装置布设坐标矩阵C，该矩阵第i行、第j列(i,j＝1,2,3,…,N)的元素表示第i行、第j爆场地表掠过压测量装置的坐标，即布置位置：

C_i,j＝[x_i,y_j]＝[-L/2+△L(i-1),-L/2+△L(j-1)]

步骤二：动爆场地表掠过压测量装置结构安装

将N×N个动爆场地表掠过压测量装置逐一安装在预设布设位置，安装后确保测压装置上表面与测试区域地表平齐；通过细钢钎穿过安装定位孔将各个动爆场地表掠过压测量装置结构固定在地表；

步骤三：试验实施

开展动爆毁伤试验，记录炸点坐标[x₀,y₀]；

步骤四：动爆场地表掠过压测量装置数据测量及处理

逐一测量N×N个测点位置处测压装置上M个膜片的特征挠度值，并根据M个膜片的平均特征挠度计算该测压装置的压力值，并形成N行N列的压力矩阵P，矩阵P第i行、第j列(i,j＝1,2,3,…,N)的元素表示第i行、第j列动爆场地表掠过压测量装置测量的掠过压；

步骤五：威力场构建

给定威力场特征压力序列P_c＝[P₁,P₂,…,P_n]，通过插值算法计算等压线簇：

Φ(x,y)＝[f₁(x,y)＝C₁,f₂(x,y)＝C₁,…,f_n(x,y)＝C_n]；

其中f_i(x,y)＝C_i上的压力为P_i。

Images