CN112214896A - 一种活性破片战斗部毁伤威力场的获取方法 - Google Patents
一种活性破片战斗部毁伤威力场的获取方法 Download PDFInfo
- Publication number
- CN112214896A CN112214896A CN202011093810.8A CN202011093810A CN112214896A CN 112214896 A CN112214896 A CN 112214896A CN 202011093810 A CN202011093810 A CN 202011093810A CN 112214896 A CN112214896 A CN 112214896A
- Authority
- CN
- China
- Prior art keywords
- fragment
- fragments
- warhead
- active
- distribution
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种活性破片战斗部毁伤威力场的获取方法,基于惰性破片战斗部毁伤威力场建立相关理论,结合单枚活性破片碰撞靶板毁伤性能分析,提出了一种活性破片战斗部毁伤威力场的研究方法。本发明的方法可以分析战斗部结构参数和终点弹道对毁伤威力场的影响,为战斗部设计和终点弹道优化提供指导和参考。
Description
技术领域
本发明属于高效毁伤和活性材料应用技术领域,具体涉及一种活性破片战斗部毁伤威力场的获取方法。
背景技术
活性破片一般可由氟聚物基体和金属颗粒经混合、模压成型和烧结硬化等工艺制备而成,在碰撞和侵彻目标过程中强动载作用下,活性破片将会被激活发生化学反应,特别是在贯穿目标防护层后,可实现对靶后目标的高效毁伤。
正是由于活性破片这种独特的热力耦合响应,使得其在高效毁伤领域的研究越来越广泛,受到世界各国的大力研究,特别是在撞击起爆活性材料配方制备、力学性能、能量释放特性、撞击起爆敏感性以及毁伤效应等方面取得了巨大进展。目前,受活性破片强度和爆炸驱动安定性的影响,活性破片在战斗部中的应用技术还未成熟。但是,毫无疑问的是,利用活性破片替代现有战斗部中的惰性金属破片,是实现战斗部毁伤威力提升的有效途径之一,特别是战斗部的引燃和引爆能力。另外,受经费和实验条件限制,采用实验测试方法研究活性破片战斗部毁伤威力往往难以获得规律性数据。因此,如何通过现有理论和相关基础实验,科学计算活性破片战斗部毁伤威力场,具有重要的意义。
发明内容
发明目的:针对现有技术中由于各方面条件的限制导致活性破片战斗部毁伤威力的研究较为困难的问题,本发明提供一种活性破片战斗部毁伤威力场的获取方法。
技术方案:一种活性破片战斗部毁伤威力场的获取方法,包括以下步骤:
(1)判断目标是否为易燃/易爆物,根据判断结果确定单枚活性破片的最小杀伤动能;
(2)计算战斗部活性破片分布的概率密度曲线,由活性破片分布的概率密度曲线计算破片的动态分布曲线,根据破片的动态分布曲线及战斗部的起爆方式计算单枚活性破片的速度,得到战斗部活性破片的飞散特性;
(3)基于破片速度衰减规律,并结合步骤(1)得到的单枚活性破片的最小杀伤动能,计算破片最大杀伤距离;建立直角坐标系,并根据步骤(2)得到的活性破片的飞散特性及破片最大杀伤距离绘制活性破片打击轨迹线,得到活性破片战斗部毁伤威力场。
进一步地,所述步骤(1)中,单枚活性破片的最小杀伤动能的确定方法为:
设活性破片的最小杀伤动能为Emin,
若目标为非易燃/易爆物,则Emin为相同尺寸和质量惰性金属破片对目标的最小杀伤动能;
若目标为易燃/易爆物,则Emin为贯穿目标防护层靶板的最小能量,即:
其中,vs为弹道极限速度;mp为活性破片质量;h为靶板厚度;A为活性破片平均着靶面积;k,α,β为实验系数。
进一步地,所述步骤(2)中,计算破片分布的概率密度曲线的方法为:
进一步地,所述步骤(2)中,计算破片的动态分布曲线的具体方法为:
利用误差函数计算破片分布的积分函数为:
设vc为战斗部的终点速度,vx为战斗部壳体上某一微元的速度,则有:
进一步地,所述步骤(2)中,根据起爆方式计算单枚活性破片的速度,具体方法为:
若起爆方式为两端起爆,设药柱的长度为2l,则单枚破片的速度为:
若起爆方式为中心起爆,设药柱的长度为2l,则单枚破片的速度为:
若起爆方式为一端起爆,设药柱的长度为l,则单枚破片的速度为:
其中,D为炸药爆速;α=x/l,lx为壳体各微元离引爆点的距离,r0为药柱的半径,r为壳体破裂时半径,β为装填系数。
进一步地,所述步骤(3)中,计算破片最大杀伤距离的方法为:
活性破片战斗部的破片总数为:
式中,N为破片总数;M为战斗部总质量;β为装填系数;m1为外壳和缓冲层质量;mp为破片质量;
单枚破片速度衰减规律为
对于战斗部中所有活性破片均做公式(14)的计算,即nx取值范围为1≤nx≤N,则可得到破片战斗部中所有破片的最大杀伤距离。
进一步地,所述步骤(3)中,绘制活性破片打击轨迹线的方法为:
在空间中建立直角坐标系,假设战斗部在空间爆炸时的坐标为(x0,y0,z0),通过步骤(2)破片飞散特性分析,得到第nx枚破片的飞散方向,假设其速度矢量在该坐标系中可表示为(Vx,Vy,Vz),绘制射击线时,起点为战斗部位置(x0,y0,z0),设终点坐标为(xR,yR,zR),坐标值可表示为
对于战斗部中所有活性破片均做公式(15)的计算,则可得到战斗部中所有活性破片的终点位置,根据战斗部所有活性破片的终点位置绘制整个活性破片战斗部的毁伤威力场。
有益效果:相比较现有技术,本发明提供一种活性破片战斗部毁伤威力场的获取方法,可以分析战斗部结构参数和终点弹道对毁伤威力场的影响,为战斗部设计和终点弹道优化提供指导和参考。通过该方法可计算活性破片战斗部毁伤威力场,从而提高计算效率,而且,可以获得规律性的数据,优化设计战斗部结构和终点弹道参数,不受经费和实验条件限制。
附图说明
图1为活性破片战斗部毁伤威力场的获取方法的流程图;
图2为活性破片战斗部的结构示意图;
图3(a)为单枚破片飞散示意图;
图3(b)为战斗部破片飞散区域示意图;
图4(a)为落角为90°的活性破片战斗部毁伤威力场的绘制图;
图4(b)为落角为60°的活性破片战斗部毁伤威力场的绘制图;
图4(c)为落角为45°的活性破片战斗部毁伤威力场的绘制图;
具体实施方式
下面根据具体实施方式,进一步阐述本发明。
如图1所示,一种活性破片战斗部毁伤威力场的获取方法,包括以下三部分内容:
(1)单枚活性破片毁伤能力研究:
设活性破片的最小杀伤动能为Emin,判断目标是否为易燃/易爆物,根据判断结果确定单枚活性破片作用不同类型目标的最小杀伤动能;当活性破片的能量大于等于最小杀伤动能时,可以对目标造成毁伤。
活性破片在作用目标过程中是通过动能碰撞和化学能释放的共同作用实现对目标的联合毁伤,若目标为如人员、雷达等非易燃/易爆物时,可以认为在杀伤非易燃/易爆性目标时,不考虑爆燃化学反应的影响,则单枚活性破片的杀伤准则与惰性金属破片相当,即Emin为相应(同尺寸和质量)惰性金属破片对目标的最小杀伤动能;
然而,当活性破片作用易燃/易爆类目标时,由于活性破片在碰撞过程中会被激活发生爆燃/爆炸化学反应,可有效提高引燃、引爆能力。因此,近似认为当活性破片贯穿防护层靶板后,爆燃/爆炸化学反应即可完成引燃、引爆,实现对目标的毁伤。也就是说,活性破片毁伤易燃/易爆类目标所需的最小动能要比惰性金属破片小得多。若目标为易燃/易爆物,则Emin为能够贯穿目标防护层靶板的最小能量,即:
而活性破片贯穿等效铝靶的弹道极限速度为:
其中,vs为弹道极限速度(m/s);mp为活性破片质量(g);h为靶板厚度(cm);A为活性破片平均着靶面积(cm2);k,α,β为实验系数。本实施例的活性破片采用活性破片侵彻铝板制作,k,α,β取值分别为1855.7,0.4143,-0.5549。
(2)破片空间飞散特性研究:
活性破片战斗部示意如图2所示,对称轴5往外依次为炸药4、缓冲层3、活性破片2及壳体1。需要注意的是,缓冲层主要用于防止活性破片在驱动过程中发生破碎或反应。设起爆点在战斗部的轴线上,破片在战斗部径向上分布均匀,且静态时破片在轴向的半平面上服从正态分布,正态分布的参数(期望和方差)由破片静态飞散角以及静态方向角确定,则静态破片分布密度函数为:
得到破片分布的积分曲线利用求出破片的动态分布曲线实际战斗部的破片初速根据起爆方式的不同沿着战斗部的壳体分布不同,这必然对战斗部破片的动态分布带来影响。若考虑破片初速沿壳体的变化,设vc为战斗部的终点速度,vx为战斗部壳体上某一微元的速度,则有:
战斗部的起爆方式不同,vx的分布关系也有所差别,下面根据破片的动态分布曲线及战斗部的起爆方式计算单枚活性破片的速度,得到活性破片的飞散特性:
若起爆方式为两端起爆,设药柱的长度为2l,则单枚破片的速度为:
若起爆方式为中心起爆,设药柱的长度为2l,则单枚破片的速度为:
若起爆方式为一端起爆,设药柱的长度为l,则单枚破片的速度为:
其中,D为炸药爆速;α=x/l,lx为壳体各微元离引爆点的距离,r0为药柱的半径,r为壳体破裂时半径,β为装填系数。
(3)战斗部毁伤威力场绘制:
活性破片战斗部的破片总数为:
式中,N为破片总数;M为战斗部总质量;β为装填系数;m1为外壳和缓冲层质量;mp为破片质量;
单枚破片速度衰减规律为
按照上述方法,通过循环计算的方法,即nx取值范围为1≤nx≤N,则可得到破片战斗部中所有破片的最大杀伤距离。
在空间中建立直角坐标系,假设战斗部在空间爆炸时的坐标为(x0,y0,z0),通过步骤(2)破片飞散特性分析,得到第nx枚破片的飞散方向,假设其速度矢量在该坐标系中可表示为(Vx,Vy,Vz),绘制射击线时,起点为战斗部位置(x0,y0,z0),设终点坐标为(xR,yR,zR),坐标值可表示为
按照上述方法,通过循环计算的方法,则可得到战斗部中所有活性破片的终点位置,从而绘制出整个破片战斗部的毁伤威力场。
本实施例选取战斗部质量为100kg,破片质量为10g的活性破片战斗部进行了实验,破片形状为球形,爆高为20m,战斗部终点速度为300m/s,起爆方式为中心起爆,壳体和缓冲层质量为20kg,最小杀伤动能为80J,以落角为90°、60°、45°三种条件绘制出活性破片毁伤威力场,如图4(a)、图4(b)、图4(c)所示,该毁伤场可较为清晰的表示破片战斗部的打击范围。
Claims (7)
1.一种活性破片战斗部毁伤威力场的获取方法,其特征在于,包括以下步骤:
(1)判断目标是否为易燃/易爆物,根据判断结果确定单枚活性破片的最小杀伤动能;
(2)计算战斗部活性破片分布的概率密度曲线,由活性破片分布的概率密度曲线计算破片的动态分布曲线,根据破片的动态分布曲线及战斗部的起爆方式计算单枚活性破片的速度,得到战斗部活性破片的飞散特性;
(3)基于破片速度衰减规律,并结合步骤(1)得到的单枚活性破片的最小杀伤动能,计算破片最大杀伤距离;建立直角坐标系,并根据步骤(2)得到的活性破片的飞散特性及破片最大杀伤距离绘制活性破片打击轨迹线,得到活性破片战斗部毁伤威力场。
6.根据权利要求5所述的活性破片战斗部毁伤威力场的获取方法,其特征在于,所述步骤(3)中,计算破片最大杀伤距离的方法为:
活性破片战斗部的破片总数为:
式中,N为破片总数;M为战斗部总质量;β为装填系数;m1为外壳和缓冲层质量;mp为破片质量;
单枚破片速度衰减规律为
对于战斗部中所有活性破片均做公式(14)的计算,即nx取值范围为1≤nx≤N,则可得到破片战斗部中所有破片的最大杀伤距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011093810.8A CN112214896B (zh) | 2020-10-14 | 2020-10-14 | 一种活性破片战斗部毁伤威力场的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011093810.8A CN112214896B (zh) | 2020-10-14 | 2020-10-14 | 一种活性破片战斗部毁伤威力场的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112214896A true CN112214896A (zh) | 2021-01-12 |
CN112214896B CN112214896B (zh) | 2022-09-16 |
Family
ID=74053340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011093810.8A Active CN112214896B (zh) | 2020-10-14 | 2020-10-14 | 一种活性破片战斗部毁伤威力场的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112214896B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114018109A (zh) * | 2021-09-26 | 2022-02-08 | 中国人民解放军63921部队 | 一种杀伤破甲子母弹靶场飞行试验落区危险区确定方法 |
CN114882354A (zh) * | 2022-04-15 | 2022-08-09 | 北京理工大学 | 一种活性毁伤元爆裂毁伤快速精确图像识别方法 |
CN115060123A (zh) * | 2022-05-07 | 2022-09-16 | 中国人民解放军国防科技大学 | 一种活性破片综合威力测试装置及测试方法 |
CN116822378A (zh) * | 2023-08-25 | 2023-09-29 | 中国人民解放军陆军装甲兵学院 | 基于回归分析的活性芯体弹靶碰撞效应预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177191A (zh) * | 2013-04-11 | 2013-06-26 | 北京理工大学 | 一种战斗部毁伤效能评估方法 |
RU2520584C1 (ru) * | 2013-02-06 | 2014-06-27 | Открытое акционерное общество "Атомэнергопроект" | Устройство для защиты внутреннего объема помещения от осколков поврежденной стены ударом извне |
CN104143028A (zh) * | 2014-08-12 | 2014-11-12 | 中国人民解放军装甲兵工程学院 | 一种基于仿真实验的装甲装备毁伤规律分析方法 |
CN108563893A (zh) * | 2018-04-24 | 2018-09-21 | 南京理工大学 | 杀爆类战斗部破片流密度计算方法 |
CN109063331A (zh) * | 2018-08-02 | 2018-12-21 | 贵州理工学院 | 一种小脱靶量破片战斗部对再入弹头毁伤效应的评估方法 |
-
2020
- 2020-10-14 CN CN202011093810.8A patent/CN112214896B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2520584C1 (ru) * | 2013-02-06 | 2014-06-27 | Открытое акционерное общество "Атомэнергопроект" | Устройство для защиты внутреннего объема помещения от осколков поврежденной стены ударом извне |
CN103177191A (zh) * | 2013-04-11 | 2013-06-26 | 北京理工大学 | 一种战斗部毁伤效能评估方法 |
CN104143028A (zh) * | 2014-08-12 | 2014-11-12 | 中国人民解放军装甲兵工程学院 | 一种基于仿真实验的装甲装备毁伤规律分析方法 |
CN108563893A (zh) * | 2018-04-24 | 2018-09-21 | 南京理工大学 | 杀爆类战斗部破片流密度计算方法 |
CN109063331A (zh) * | 2018-08-02 | 2018-12-21 | 贵州理工学院 | 一种小脱靶量破片战斗部对再入弹头毁伤效应的评估方法 |
Non-Patent Citations (9)
Title |
---|
何正文等: "基于战斗部数值仿真系统的杀伤概率计算", 《航空兵器》 * |
初善勇: "杀伤爆破弹毁伤威力等效评估研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
史志鑫等: "预制破片的形状对破片飞散性能影响的数值模拟研究", 《兵器装备工程学报》 * |
段丰安等: "爆破子母弹战斗部毁伤威力场研究", 《战术导弹技术》 * |
罗振华等: "某型活性破片毁伤后效试验研究与分析", 《科学技术与工程》 * |
胡杰等: "超高速破片发射系统内弹道参数设计及验证", 《西安工业大学学报》 * |
赵传等: "反舰导弹战斗部近距空爆的毁伤威胁性分析", 《中国舰船研究》 * |
邱从礼等: "压制火炮对暴露有生力量毁伤效能评估方法", 《火力与指挥控制》 * |
马艳丽: "杀爆战斗部对典型雷达毁伤评估研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114018109A (zh) * | 2021-09-26 | 2022-02-08 | 中国人民解放军63921部队 | 一种杀伤破甲子母弹靶场飞行试验落区危险区确定方法 |
CN114018109B (zh) * | 2021-09-26 | 2023-08-11 | 中国人民解放军63921部队 | 一种杀伤破甲子母弹靶场飞行试验落区危险区确定方法 |
CN114882354A (zh) * | 2022-04-15 | 2022-08-09 | 北京理工大学 | 一种活性毁伤元爆裂毁伤快速精确图像识别方法 |
CN114882354B (zh) * | 2022-04-15 | 2024-05-31 | 北京理工大学 | 一种活性毁伤元爆裂毁伤快速精确图像识别方法 |
CN115060123A (zh) * | 2022-05-07 | 2022-09-16 | 中国人民解放军国防科技大学 | 一种活性破片综合威力测试装置及测试方法 |
CN115060123B (zh) * | 2022-05-07 | 2024-01-30 | 中国人民解放军国防科技大学 | 一种活性破片综合威力测试装置及测试方法 |
CN116822378A (zh) * | 2023-08-25 | 2023-09-29 | 中国人民解放军陆军装甲兵学院 | 基于回归分析的活性芯体弹靶碰撞效应预测方法 |
CN116822378B (zh) * | 2023-08-25 | 2023-12-08 | 中国人民解放军陆军装甲兵学院 | 基于回归分析的活性芯体弹靶碰撞效应预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112214896B (zh) | 2022-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112214896B (zh) | 一种活性破片战斗部毁伤威力场的获取方法 | |
CN107798208B (zh) | 对空目标导弹破片飞散最大毁伤算法 | |
CN114492170A (zh) | 一种典型目标易损性毁伤评估方法及系统 | |
CN108563893B (zh) | 杀爆类战斗部破片流密度计算方法 | |
TR201816245T4 (tr) | Gelişmiş parça tesirli delici mühimmat. | |
Van Der Voort et al. | A statistical description of explosion produced debris dispersion | |
US7387060B1 (en) | Rocket exhaust defense system and method | |
Deniz | Ballistic penetration of hardened steel plates | |
CN113468752A (zh) | 考虑作战需求的多轴特种车辆轮胎系统抗毁伤能力评估方法 | |
CN114508968B (zh) | 复合结构含能微弹丸毁伤云的制备方法与评估方法 | |
Zecevic et al. | Lethal influence factors of natural and preformed fragmentation projectiles | |
Ji et al. | Dynamic behaviors of multi-layered steel targets with air gaps subjected to the impact of EFP simulants | |
CN113609574B (zh) | 基于等毁线的多轴特种车辆轮胎系统抗毁伤能力表征方法 | |
TIAN et al. | Attitude deflection mechanism of projectiles with variable elliptical cross-sections obliquely perforating thin targets | |
Chen et al. | Shock-induced detonation of high explosives by high velocity impact | |
Bai et al. | Analysis of criteria for assessing safety distance for focused warhead fragments based on CDEM | |
CN109740266B (zh) | 一种基于复杂效应场计算的地下结构毁伤评估方法 | |
CN110081783B (zh) | 一种低损伤炮射试验弹及其回收方法 | |
Choudha et al. | Parametric study of single confined fragment launch explosive device | |
JP6025525B2 (ja) | 弾頭 | |
Petkov et al. | Optimization of the Design of Fragmentation Warheads Forming an Axial Flow of Preformed Fragments | |
Shi et al. | The Research on Damage Effect of Target under the Action of Fragmentation Warhead | |
CN114508970A (zh) | 高强度含能微弹丸毁伤云的制备方法与评估方法 | |
Chen et al. | Simulation analysis of equipment battle damage test layout design | |
CN214892907U (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 |