CN112270140A - 一种用于爆炸冲击波超压预测的远场反演近场分布的方法 - Google Patents
一种用于爆炸冲击波超压预测的远场反演近场分布的方法 Download PDFInfo
- Publication number
- CN112270140A CN112270140A CN202011047335.0A CN202011047335A CN112270140A CN 112270140 A CN112270140 A CN 112270140A CN 202011047335 A CN202011047335 A CN 202011047335A CN 112270140 A CN112270140 A CN 112270140A
- Authority
- CN
- China
- Prior art keywords
- field
- far
- overpressure
- shock wave
- explosion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000035939 shock Effects 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000009826 distribution Methods 0.000 title claims abstract description 50
- 239000002360 explosive Substances 0.000 title abstract description 15
- 238000004880 explosion Methods 0.000 claims abstract description 75
- 238000006243 chemical reaction Methods 0.000 claims abstract description 22
- 230000009466 transformation Effects 0.000 claims description 33
- 238000005474 detonation Methods 0.000 claims description 17
- 238000010606 normalization Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 9
- 230000001131 transforming effect Effects 0.000 claims description 5
- 238000012417 linear regression Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 10
- 238000004088 simulation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000001154 acute effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009863 impact test Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种用于爆炸冲击波超压预测的远场反演近场分布的方法，将爆炸冲击波远场数据对处理后得到变换后数据对，将冲击波峰值超压与比例距离的三次多项式用等效转换因子进行归一化，将所得到空间归一化三次多项式，与构造的冲击波近场空间等效转换因子方程相乘，得到冲击波近场不同空间位置峰值超压反演模型。本发明是为了解决冲击波峰值超压大尺度远近场分布模型构建的问题，提供一种用于爆炸冲击波峰值超压预测的远场数据反演近场分布的方法，将远近场变换方程和近场空间等效转换模型综合起来，在已知比例爆高时建立起近场峰值超压不同空间位置的关系，从而通过测量远场峰值超压反演出近场超压分布，满足常规和复杂工况毁伤评估需求。
Description
技术领域
本发明涉及冲击测试技术领域，具体涉及一种用于爆炸冲击波超压预测的远场反演近场分布的方法。
背景技术
爆炸冲击波峰值超压是高能战斗部杀伤目标的重要毁伤元之一。定量测试实际战场上战斗部冲击波超压数据，既可指导弹药研制又有助于分析毁伤效果。
学者们对空中爆炸冲击波峰值超压传播规律进行了深入研究。早在1915年，Hopkinson B.就根据试验现象和数据分析得出爆炸具有相似性的结论。随后，研究人员进行了大量的TNT爆炸试验，分析收集到的峰值超压数据，得到一系列不同适应范围的峰值超压与比例距离的关系式。包括Mills公式、Henrych公式、Sadovskyi公式、Brode公式、叶晓华公式及国防工程设计规范等在内的众多峰值超压一维分布规律的经验公式，对工程实际应用起到很好的指导作用。上述公式中峰值超压仅是比例距离的多项式关系，并基于相同比例距离时冲击波峰值超压处处一致的假设，未考虑不同空间分布的测点处峰值超压不同。这种冲击波峰值超压的空间非均匀性在满足无限空中爆炸条件下，比例爆高(即爆心离地面高度与当量三次根之比，单位m·kg-1/3)越小越明显。
构建上述经验公式的数据主要来自于自由场测试结果，要求炸药比例爆高不小于0.35m·kg-1/3，且以地面和近地自由场布设测点为主。爆炸冲击入射波会在抵达地面后发生发射，并与入射波叠加形成马赫波。三种波叠加为三波点导致波阵面的分布不均匀。随着毁伤评估需求的发展，典型战斗部爆炸冲击波的测试方法逐渐向适应特殊工况环境方向发展。常规地面和近地场布设测点的冲击波压力测量的方式已不再适应新形势下的战斗部压力毁伤元测试。通过空中测点实现冲击波压力测量的方法可满足上述测试需求，但对冲击波峰值超压大尺度远近场分布模型的构建提出迫切需求。
发明内容
本发明是为了解决冲击波峰值超压大尺度远近场分布模型构建的问题，提供一种用于爆炸冲击波峰值超压预测的远场数据反演近场分布的方法，将远近场变换方程和近场空间等效转换模型综合起来，在已知比例爆高时建立起近场峰值超压不同空间位置的关系，从而通过测量远场峰值超压反演出近场超压分布数据，满足常规和复杂工况毁伤评估需求。
本发明提供一种用于爆炸冲击波超压预测的远场反演近场分布的方法，包括以下步骤：
S1、获取爆炸冲击波远场不同比例距离的多个测点的数据对，利用远近场变换方程对数据对进行变换，得到变换后数据对；
S2、基于爆炸相似律和多项式回归模型，构造冲击波峰值超压与比例距离的三次多项式；
将变换后数据对代入三次多项式，利用最小二乘法求解，得到三次多项式系数；
S3、计算每个测点与爆心连线后相对于与地面平行的水平面的夹角，并计算得到夹角的平均值α；
计算比例爆高，代入近场空间等效转换模型，得到七种角度对应的等效转换因子，通过线性插值法得到以角度为自变量的等效转换因子方程；
将α带入等效转换因子方程，得到冲击波近场空间等效转换因子；
S4、使用冲击波近场空间等效转换因子修正三次多项式，得到空间归一化三次多项式；
S5、将空间归一化三次多项式与等效转换因子方程相乘，得到冲击波近场不同空间位置峰值超压反演模型。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S1中，爆炸为近地空中爆炸，满足无限空中爆炸条件，爆炸条件为爆炸高度与当量三次根之比不小于0.35m·kg-1/3。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S1中，测点至少为四组；测点位置分布在与爆心连线的同一条直线上，测点距爆心的比例距离与当量三次根之比不小于15m·kg-1/3，比例距离是指测点和爆心的直线距离与当量三次根之比；测点峰值超压结果为正值。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S1中，近场是比例距离1-15m·kg-1/3；
数据对由比例距离及其对应峰值超压构成；其中，
表示远场比例距离，单位为m·kg-1/3，Δpd表示对应峰值超压，单位为kPa，第i组远场数据对为
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S1中远近场变换方程如下：
Δpdi′＝1.295·Δpdi
远近场变换方程包括比例修正和坐标变换；比例修正因子为1.295，只对峰值超压起作用；坐标变换通过使比例修正后的数据对
逆时针旋转15°实现。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S2中三次多项式为：
其中，
表示比例距离，单位为m·kg-1/3，Δpm表示峰值超压，单位为kPa,用A1、A2和A3表示无量纲的三次多项式系数。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S3中夹角为0°至90°。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S3中空间等效转换模型是基于线性回归模型构造的比例爆高与空间各角度处峰值超压归一化值的关系模型，用线性方程表示；
空间等效转换模型包含7组线性方程，每组线性方程对应一个空间角度，分别为0°、15°、30°、45°、60°、75°和90°。每个角度对应的峰值超压归一化值的方程如下所示，
其中，x表示比例爆高，即爆心离地面高度与当量三次根之比，单位为m·kg-1/3，角标数字代表空间角度。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S3中的峰值超压归一化值是各角度处同一比例距离的峰值超压与90°对应峰值超压之比。
本发明所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，作为优选方式，步骤S4中冲击波近场空间等效转换因子修正三次多项式的方式是：将三次多项式除以等效转换因子。
本方法适用于近地空中爆炸远场超压毁伤元测试结果的应用，也适用于近地空中爆炸近场超压毁伤元空间分布的求解，对于特殊原因无法提前入场布设超压测点而只能通过远场测试获取超压数据的场合尤其适用。
本发明具有以下优点：
本发明与现有技术相比的优点在于：
(1)本发明针对近地空中爆炸冲击波峰值超压远近场关系首次提出了远近场变换方程，利用该方程将远场峰值超压变换为符合近场分布规律的数据。
(2)本发明针对近地空中爆炸冲击波峰值超压近场空间分布规律首次提出了近场空间等效转换模型，利用该模型，在已知比例爆高时建立起近场峰值超压不同空间位置的关系。
(3)本发明针对近地空中爆炸冲击波峰值超压远场数据反演近场峰值超压分布提出了将远近场变换方程和近场空间等效转换模型综合起来，利用该综合模型，在已知比例爆高和远场不同比例距离同一直线分布的至少四处正峰值超压数据时，可反演比例距离1-19m·kg-1/3内的任意空间位置峰值超压；在2-15m·kg-1/3内且在45°内的空间位置预测结果较准确，其余位置可供参考。
附图说明
图1为一种用于爆炸冲击波超压预测的远场反演近场分布的方法流程图；
图2为一种用于爆炸冲击波超压预测的远场反演近场分布的方法——利用AUTODYN对100kg TNT爆炸高度3m的典型工况数值仿真至98ms时的截图；
图3为一种用于爆炸冲击波超压预测的远场反演近场分布的方法——测点设置及夹角示意图；
图4为一种用于爆炸冲击波超压预测的远场反演近场分布的方法——在典型工况500kg TNT爆高4.5m且同一角度时采用远近场变换方程对远场数据变换前后结果，以变换后远场数据为基础构建远场模型以及以近场数据为基础构建模型的对比图；
图5为一种用于爆炸冲击波超压预测的远场反演近场分布的方法——在比例爆高0.496m·kg-1/3时采用近场空间等效转换模型对近场空间各角度峰值超压转换为90°处效果图；
图6为一种用于爆炸冲击波超压预测的远场反演近场分布的方法——对典型工况500kg TNT爆高4.5m的远场数据采用综合模型处理后相对于原始仿真结果的误差图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
实施例1
如图1所示，一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：包括以下步骤：
S1、获取爆炸冲击波远场不同比例距离的多个测点的数据对，利用远近场变换方程对数据对进行变换，得到变换后数据对；
S2、基于爆炸相似律和多项式回归模型，构造冲击波峰值超压与比例距离的三次多项式；
将变换后数据对代入三次多项式，利用最小二乘法求解，得到三次多项式系数；
S3、计算每个测点与爆心连线后相对于与地面平行的水平面的夹角，得到平均值α；
计算比例爆高，代入近场空间等效转换模型，得到七种角度对应的等效转换因子，通过线性插值法得到以角度为自变量的等效转换因子方程；
将α带入等效转换因子方程，得到冲击波近场空间等效转换因子；
S4、使用冲击波近场空间等效转换因子修正三次多项式，得到空间归一化三次多项式；
S5、将空间归一化三次多项式与冲击波近场空间等效转换方程相乘，得到冲击波近场不同空间位置峰值超压反演模型。
实施例2
如图1所示，一种用于爆炸冲击波超压预测的远场反演近场分布的方法，包括以下步骤：
S1、获取爆炸冲击波远场不同比例距离的至少四组测点的数据对，测点位置分布在与爆心连线的同一条直线上，测点距爆心的比例距离与当量三次根之比不小于15m·kg-1/3，比例距离是指测点和爆心的直线距离与当量三次根之比；测点峰值超压结果为正值；
利用远近场变换方程对数据对进行变换，得到变换后(比例距离1-15m·kg-1/3)数据对；
爆炸为近地空中爆炸，满足无限空中爆炸条件，爆炸条件为爆炸高度与当量三次根之比不小于0.35m·kg-1/3；
数据对由比例距离及其对应峰值超压构成；其中，
表示远场比例距离，单位为m·kg-1/3，Δpd表示对应峰值超压，单位为kPa，第i组远场数据对为
远近场变换方程如下：
Δpdi′＝1.295·Δpdi
远近场变换方程包括比例修正和坐标变换；比例修正因子为1.295，只对峰值超压起作用；坐标变换通过使比例修正后的数据对
逆时针旋转15°实现；
S2、基于爆炸相似律和多项式回归模型，构造冲击波峰值超压与比例距离的三次多项式；
三次多项式为：
其中，
表示比例距离，单位为m·kg-1/3，Δpm表示峰值超压，单位为kPa,用A1、A2和A3表示无量纲的三次多项式系数；
将变换后数据对代入三次多项式，利用最小二乘法求解，得到三次多项式系数；
S3、计算每个测点与爆心连线后相对于与地面平行的水平面的夹角(0°至90°)，并计算得到夹角的平均值α；
计算比例爆高，代入近场空间等效转换模型，得到七种角度对应的等效转换因子，通过线性插值法得到以角度为自变量的等效转换因子方程；
空间等效转换模型是基于线性回归模型构造的比例爆高与空间各角度处峰值超压归一化值的关系模型，用线性方程表示；
峰值超压归一化值是各角度处同一比例距离的峰值超压与90°对应峰值超压之比。
空间等效转换模型包含7组线性方程，每组线性方程对应一个空间角度，分别为0°、15°、30°、45°、60°、75°和90°。每个角度对应的峰值超压归一化值的方程如下所示，
其中，x表示比例爆高，即爆心离地面高度与当量三次根之比，单位为m·kg-1/3，角标数字代表空间角度；
将α带入等效转换因子方程，得到冲击波近场空间等效转换因子；
S4、使用冲击波近场空间等效转换因子修正三次多项式，得到空间归一化三次多项式；
冲击波近场空间等效转换因子修正三次多项式的方式是：将三次多项式除以等效转换因子。
S5、将空间归一化三次多项式与等效转换因子方程相乘，得到冲击波近场不同空间位置峰值超压反演模型。
实施例3
如图1-6所示，本发明的实现步骤如下：
S1、获取爆炸冲击波远场不同比例距离的至少四组数据对；
如图2所示，利用AUTODYN仿真软件建立比例爆高大于0.35m·kg-1/3的近地等效TNT空中爆炸数值模型，设置好与爆心连线在同一直线上的比例距离不同的四个高斯点，上述四个高斯点比例距离均大于15m·kg-1/3。数值仿真时间开销较大，设置较多测点可在一次仿真结束后获得较多测试数据。
利用远近场变换方程对远场数据对进行变换，得到变换后数据对；
近场为比例距离1-15m·kg-1/3的范围，而远场为比例距离15m·kg-1/3以外的范围，数据对由比例距离及其对应峰值超压组成。变换操作的对象为远场数据对，操作方法为将远场数据对代入远近场变换方程中运算，操作结果是得到符合近场分布规律的变换后数据对。
远近场变换方程由两部分组成，一是比例修正，二是坐标变换。比例修正因子为1.295，坐标变换通过逆时针旋转15°实现。
用
表示远场比例距离，单位为m·kg-1/3，Δpd表示对应峰值超压，
比例修正后的远场中心数据对。第i组远场数据对的完整修正方程如下所示：
Δpdi′＝1.295·Δpdi
远场中心数据对按如下方法选取：当远场数据对数为偶数时，取数据量的一半对应位置的数据对为中心数据对；当远场数据对数量为奇数时，取数量加1后一半的位置的数据对为中心数据对。
以500kg TNT爆高4.5m的工况为例，其对应四组远场数据对为(15.886，3.093)、(17.147，2.319)、(18.408，1.673)和(19.290，1.28)。远场中心数据对为(17.147，2.319)，变换后数据对为(15.659，3.620)、(17.147，3.004)、(18.590，2.547)和(19.578，2.301)。
S2、基于爆炸相似律和多项式回归模型，构造冲击波峰值超压与比例距离的三次多项式，系数待定；
用r表示比例距离，单位为m·kg-1/3，Δpm表示峰值超压，单位为kPa,用A1、A2和A3表示无量纲的系数，则构造的三次多项式形式为：
将步骤S1中得到的变换后数据对带入步骤S2上述的三次多项式中，利用最小二乘法求解系数；
步骤S1中将产生至少四对数据对，是求解三次多项式系数的最低要求。
以
和
表示四组远近场变换后的变换后数据对，采用变量替换方法，对步骤S2中的三次多项式进行改造，然后将四组数据对代入多项式中，最后将因变量、自变量和待定系数分别表示，即：
Δpdr′＝[Δpdr1′ Δpdr2′ Δpdr3′ Δpdr4′]T
A＝[A1 A2 A3]T
待求解问题化为求下列方程的极小解：
即求下列法方程的解即可得到系数矩阵：
求解500kg TNT爆高4.5m的工况，得到三个系数分别-7.647，1007，116.7，拟合最大误差为3.04％。
S3、计算每个测点与爆心连线后相对于水平面(与地面平行)的夹角，并求均值得到α；
如图3所示，夹角为测点与爆心连线与水平面所成锐角，范围为0°至90°，含两个边界角度值。首先计算测点与爆心连线的比例距离，然后计算测点垂直比例距离，以垂直比例距离除以连线比例距离后取反正弦运算，得到锐角夹角。最后对求出的每个锐角夹角求均值即可。
以500kg TNT爆高4.5m的工况为例，远场四个测点夹角均为0°，均值为0°。
计算比例爆高，代入近场空间等效转换模型，得到七种角度对应的等效转换因子，通过线性插值法得到以角度为自变量的等效转换因子方程；
比例爆高通过将爆心离地面高度除以炸药当量的三次根得到，500kg TNT爆高4.5m工况的比例爆高为0.567m·kg-1/3。与角度0°、15°、30°、45°、60°、75°和90°对应的等效转换因子分别为1.55794、1.37600、1.15828、1.01213、0.99916和1.01352。基于线性插值法得到的等效转换因子方程如下所示，因变量为等效转换因子，自变量为角度，由第一个方程至第六个方程，自变量范围依次为：0≤x＜15°、15°≤x＜30°、30°≤x＜45°、45°≤x＜60°、60°≤x＜75°及75°≤x≤90°：
将角度均值α代入等效转换因子方程，得到冲击波近场空间等效转换因子；
角度均值α为0°，代入等效转换因子方程得到等效转换因子为1.5794。
S4、根据步骤S3中得到的等效转换因子，对步骤S2中确定了系数的三次多项式进行修正，得到空间归一化三次多项式；
将步骤S2中确定了系数的三次多项式除以步骤S3得到的等效转换因子，从而将0°测点处的峰值超压与比例距离的关系“归一化”为90°处关系，即得到空间归一化三次多项式为：
S5、将步骤S4中空间归一化三次多项式与近场空间等效转换因子方程相乘，得到冲击波近场不同空间位置峰值超压反演模型。
基于冲击波远场测试数据的近场不同空间位置峰值超压反演模型如下所示，六个方程中角度x的范围与步骤S3中的定义一致：
典型工况500kg TNT爆高4.5m时，比例爆高为0.567m·kg-1/3，采用远近场变换方程对远场数据进行处理；基于变换后的远场数据，构建峰值超压分布模型；利用近场数据建立近场峰值超压分布规律；将转换前后数据、远场峰值超压模型和近场峰值超压模型共同展示于图4中。转换后的远场峰值超压模型对比例距离2-15m·kg-1/3内的近场峰值超压进行预测时，相对于近场数值仿真结果，最大预测误差仅为9.7％，平均绝对误差仅为4.5％；对比例距离1-2m·kg-1/3内的近场峰值超压进行预测时，最大预测误差为74.8％，平均绝对误差为63.1％，误差较大，仅供参考。
典型工况1kg TNT爆高0.5m时，比例爆高0.496m·kg-1/3，采用近场空间等效转换模型对近场空间各角度峰值超压统一转换为90°对应比例距离处数值，以90°处原数值的±25％作为上下误差线，将变换前后数据、90°峰值超压连线及误差线共同绘制在图5中，结果显示，变换前不同角度处峰值超压相对于90°处数值偏差大于+25％，变换后均处于±25％误差内，效果提升明显。
典型工况500kg TNT爆高4.5m时，比例爆高0.567m·kg-1/3，远场数据先经过变换得到变换后数据对，然后建立近场预测模型，然后预测近场不同空间位置的峰值超压，最后计算相对于原始仿真结果的误差，如图6所示。图中根据误差分布规律分为五个区域，分别标注为A-E。A区域为比例距离1-2m·kg-1/3以内，45°角以内，其误差较大，平均绝对误差约为62.7％；B区域为比例距离2-15.5m·kg-1/3以内，测点角度在45°以内，其误差分布较小且集中，平均绝对误差约为10.9％；C区域为比例距离15.5-20m·kg-1/3以内，测点角度在45°以内，其误差分布不规律，平均绝对误差约为27.9％；D区域为比例距离1-8.5m·kg-1/3以内，测点角度在45°以外，其误差不大，平均绝对误差约为16.1％；E区域为比例距离8.5-20m·kg-1/3以内，测点角度在45°以外，其误差大，平均绝对误差约为47.0％。
综合看，预测效果按照B、D、C、E和A依次递减；B、D和C区域预测结果可直接采纳，E和A区域预测误差较大，仅供参考。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：包括以下步骤：
S1、获取爆炸冲击波远场不同比例距离的多个测点的数据对，利用远近场变换方程对所述数据对进行变换，得到变换后数据对；
S2、基于爆炸相似律和多项式回归模型，构造冲击波峰值超压与比例距离的三次多项式；
将所述变换后数据对代入所述三次多项式，利用最小二乘法求解，得到所述三次多项式系数；
S3、计算每个所述测点与爆心连线后相对于与地面平行的水平面的夹角，并计算得到所述夹角的平均值α；
计算比例爆高，代入近场空间等效转换模型，得到七种角度对应的等效转换因子，通过线性插值法得到以角度为自变量的等效转换因子方程；
将α带入所述等效转换因子方程，得到冲击波近场空间等效转换因子；
S4、使用所述冲击波近场空间等效转换因子修正所述三次多项式，得到空间归一化三次多项式；
S5、将所述空间归一化三次多项式与所述等效转换因子方程相乘，得到冲击波近场不同空间位置峰值超压反演模型。
2.根据权利要求1所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：步骤S1中，所述爆炸为近地空中爆炸，满足无限空中爆炸条件，所述爆炸条件为爆炸高度与当量三次根之比不小于0.35m·kg-1/3。
3.根据权利要求1所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：步骤S1中，所述测点至少为四组；所述测点位置分布在与所述爆心连线的同一条直线上，所述测点距所述爆心的比例距离不小于15m·kg-1/3，所述比例距离是指所述测点和所述爆心的直线距离与所述当量三次根之比；所述测点峰值超压为正值。
4.根据权利要求1所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：步骤S1中，所述近场是所述比例距离1-15m·kg-1/3；
所述数据对由所述比例距离及其对应所述峰值超压构成；其中，
表示远场比例距离，单位为m·kg-1/3，Δpd表示对应峰值超压，单位为kPa，第i组所述远场数据对为
5.根据权利要求4所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：步骤S1中所述远近场变换方程如下：
Δpdi′＝1.295·Δpdi
所述远近场变换方程包括比例修正和坐标变换；所述比例修正因子为1.295，只对峰值超压起作用；所述坐标变换通过使比例修正后的数据对
逆时针旋转15°实现。
6.根据权利要求1所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：步骤S2中所述三次多项式为：
其中，
表示所述比例距离，单位为m·kg-1/3，Δpm表示所述峰值超压，单位为kPa,用A1、A2和A3表示无量纲的所述三次多项式系数。
7.根据权利要求1所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：步骤S3中所述夹角为0°至90°。
8.根据权利要求1所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：
步骤S3中所述空间等效转换模型是基于线性回归模型构造的所述比例爆高与空间各角度处所述峰值超压归一化值的关系模型，用线性方程表示；
所述空间等效转换模型包含7组线性方程，每组所述线性方程对应一个空间角度，分别为0°、15°、30°、45°、60°、75°和90°。每个角度对应的所述峰值超压归一化值的方程如下所示，
其中，x表示比例爆高，即所述爆心离地面高度与当量三次根之比，单位为m·kg-1/3，角标数字代表所述空间角度。
9.根据权利要求8所述的一种用于爆炸冲击波超压预测的远场反演近场分布的方法，其特征在于：所述峰值超压归一化值是各角度处同一所述比例距离的所述峰值超压与90°对应所述峰值超压之比。
10.根据权利要求9所述的一种用于爆炸冲击波峰值超压预测的远场数据反演近场分布的方法，其特征在于：所述步骤S4中所述冲击波近场空间等效转换因子修正所述三次多项式的方式是：将所述三次多项式除以所述等效转换因子。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011047335.0A CN112270140B (zh) | 2020-09-29 | 2020-09-29 | 一种用于爆炸冲击波超压预测的远场反演近场分布的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011047335.0A CN112270140B (zh) | 2020-09-29 | 2020-09-29 | 一种用于爆炸冲击波超压预测的远场反演近场分布的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112270140A true CN112270140A (zh) | 2021-01-26 |
| CN112270140B CN112270140B (zh) | 2023-10-17 |
Family
ID=74349329
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011047335.0A Active CN112270140B (zh) | 2020-09-29 | 2020-09-29 | 一种用于爆炸冲击波超压预测的远场反演近场分布的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112270140B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113267287A (zh) * | 2021-06-29 | 2021-08-17 | 中北大学 | 冲击波超压三维时空场重建方法 |
| CN114166399A (zh) * | 2021-11-05 | 2022-03-11 | 武汉大学 | 一种用于反演坝面水下爆炸压力的方法 |
| CN114741863A (zh) * | 2022-03-31 | 2022-07-12 | 西北核技术研究所 | 一种城市街道中偶然爆炸冲击波超压峰值的快速计算方法 |
| CN114894363A (zh) * | 2022-04-21 | 2022-08-12 | 中国人民解放军国防科技大学 | 多方位爆炸冲击波超压等效测量装置及方法 |
| CN117574730A (zh) * | 2023-11-28 | 2024-02-20 | 中国航空研究院 | 一种超声速民机机动飞行聚焦声爆的数值预测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080208545A1 (en) * | 2007-02-27 | 2008-08-28 | National Taiwan University Of Science And Technology | Non-destructive method for inverse-calculating fiber probe aperture size and prediction method of fabrication profile of near field photolithography |
| CN104820218A (zh) * | 2015-05-26 | 2015-08-05 | 中国科学院声学研究所东海研究站 | 一种基于频域自相关的浅海海底单参数反演方法 |
| CN105975713A (zh) * | 2016-05-20 | 2016-09-28 | 北京航空航天大学 | 一种采用源重建方法预测电子电路电磁辐射的方法 |
| CN108828008A (zh) * | 2018-04-19 | 2018-11-16 | 西安近代化学研究所 | 一种基于冲击波超压和准静压的炸药内爆炸威力评价方法 |
-
2020
- 2020-09-29 CN CN202011047335.0A patent/CN112270140B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080208545A1 (en) * | 2007-02-27 | 2008-08-28 | National Taiwan University Of Science And Technology | Non-destructive method for inverse-calculating fiber probe aperture size and prediction method of fabrication profile of near field photolithography |
| CN104820218A (zh) * | 2015-05-26 | 2015-08-05 | 中国科学院声学研究所东海研究站 | 一种基于频域自相关的浅海海底单参数反演方法 |
| CN105975713A (zh) * | 2016-05-20 | 2016-09-28 | 北京航空航天大学 | 一种采用源重建方法预测电子电路电磁辐射的方法 |
| CN108828008A (zh) * | 2018-04-19 | 2018-11-16 | 西安近代化学研究所 | 一种基于冲击波超压和准静压的炸药内爆炸威力评价方法 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113267287A (zh) * | 2021-06-29 | 2021-08-17 | 中北大学 | 冲击波超压三维时空场重建方法 |
| CN113267287B (zh) * | 2021-06-29 | 2022-09-23 | 中北大学 | 冲击波超压三维时空场重建方法 |
| CN114166399A (zh) * | 2021-11-05 | 2022-03-11 | 武汉大学 | 一种用于反演坝面水下爆炸压力的方法 |
| CN114166399B (zh) * | 2021-11-05 | 2023-11-10 | 武汉大学 | 一种用于反演坝面水下爆炸压力的方法 |
| CN114741863A (zh) * | 2022-03-31 | 2022-07-12 | 西北核技术研究所 | 一种城市街道中偶然爆炸冲击波超压峰值的快速计算方法 |
| CN114894363A (zh) * | 2022-04-21 | 2022-08-12 | 中国人民解放军国防科技大学 | 多方位爆炸冲击波超压等效测量装置及方法 |
| CN114894363B (zh) * | 2022-04-21 | 2022-12-16 | 中国人民解放军国防科技大学 | 多方位爆炸冲击波超压等效测量装置及方法 |
| CN117574730A (zh) * | 2023-11-28 | 2024-02-20 | 中国航空研究院 | 一种超声速民机机动飞行聚焦声爆的数值预测方法 |
| CN117574730B (zh) * | 2023-11-28 | 2024-05-10 | 中国航空研究院 | 一种超声速民机机动飞行聚焦声爆的数值预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112270140B (zh) | 2023-10-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112270140B (zh) | 一种用于爆炸冲击波超压预测的远场反演近场分布的方法 | |
| CN104992000B (zh) | 一种基于l型阵列天线的波束形成及波束图优化方法 | |
| CN102607552B (zh) | 基于神经网络的工业机器人空间网格精度补偿方法 | |
| CN112613118B (zh) | 火箭发动机内部不可测装配质量数字孪生建模与追溯方法 | |
| CN111506968B (zh) | 一种基于bp神经网络算法的船型优化方法 | |
| CN104298808A (zh) | 输电铁塔非线性柔性构件的应力计算方法 | |
| CN109141266B (zh) | 一种钢结构测量方法及系统 | |
| CN110702364B (zh) | 针对桨尖马赫数影响的高空螺旋桨风洞试验数据修正方法 | |
| CN115575104A (zh) | 基于逆有限元重构图像识别的风机叶片损伤快速检测方法 | |
| CN112258062A (zh) | 一种体系贡献度确定方法及装置 | |
| CN105785320A (zh) | 实心球阵列三维声源识别的函数型延迟求和方法 | |
| CN115994465A (zh) | 一种意外爆炸空气冲击波对桥梁毁伤效果kNN评估方法 | |
| CN113553739A (zh) | 一种混合炸药爆炸输出特性计算方法 | |
| CN109344522B (zh) | 一种用于轴流压缩机的静叶磨削量的计算方法及系统 | |
| CN110879919B (zh) | 一种爆炸作用下毒物扩散的分段式模拟方法 | |
| CN111639419A (zh) | 一种船舱内爆小变形塑性毁伤模式问题的判断方法 | |
| CN111125845A (zh) | 基于质心连线法的空间直线度误差评定方法 | |
| CN116187141A (zh) | 内爆增压环境下爆炸物威力场的仿真验证方法和装置 | |
| CN115795874A (zh) | 一种基于目标结构响应的动态炸点计算方法 | |
| CN109815577B (zh) | 落弹坐标声波定位方法及系统、计算机可读取的存储介质 | |
| CN114861315A (zh) | 一种基于机器学习的叶轮机二维叶型优化方法 | |
| CN108627629B (zh) | 圆柱体炸药轴向径向冲击波超压峰值确定方法 | |
| CN110346105B (zh) | 直接地冲击作用下地下结构响应模型试验方法 | |
| CN106250601A (zh) | 双曲线钢结构冷却塔风压极值的确定方法及装置 | |
| CN108196087B (zh) | 数据处理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |