CN111125623A - 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 - Google Patents
一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 Download PDFInfo
- Publication number
- CN111125623A CN111125623A CN201911194809.1A CN201911194809A CN111125623A CN 111125623 A CN111125623 A CN 111125623A CN 201911194809 A CN201911194809 A CN 201911194809A CN 111125623 A CN111125623 A CN 111125623A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- fragment
- launcrd
- target
- velocity
- 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
Classifications
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明提出了一种超低空作战条件下的战斗部破片飞散轨迹计算的方法,通过构建在弹目相对速度坐标系内由于空气阻力的作用的战斗部破片飞散轨迹变模型,有效解决了在超低空作战条件下战斗部破片命中目标的精度问题,并可针对大规模仿真计算进行精度设置,满足不同计算精度和计算规模的平衡。相比传统的战斗部破片飞散轨迹计算的方法,该方法能够准确仿真评估超低空作战条件下的引战配合设计、引战系统设计和武器系统作战效能。
Description
技术领域
本发明涉及战斗部破片飞散轨迹计算技术领域,特别是一种超低空作战条件下的战斗部破片飞散轨迹计算的方法。
背景技术
在导弹超低空作战时,空气密度较大,此过程中可以认为目标速度保持不变,空气阻力引起的破片速度变化不能被被忽略,破片在弹目相对速度坐标系内的飞散轨迹不再是直线,而是曲线。以往方法中,破片在弹目相对速度坐标系内的飞散速度未考虑到不断变小的问题,破片飞散轨迹仍按照匀速直线运动计算,因此以往方法不能精确计算出破片命中目标的位置,影响了引战配合设计和评估的置信度。
发明内容
本发明目的在于提供一种超低空作战条件下的战斗部破片飞散轨迹计算的方法,解决以往引战配合设计和评估置信度的问题。
本发明提出一种超低空作战条件下的战斗部破片飞散轨迹计算的方法,其步骤为:
第一步、根据引信起爆位置,初始化计算参数,第二步、根据命中检测算法,求解每段破片飞行轨迹是否命中目标,并进行下一步迭代计算。
其中,第一步、根据引信起爆位置,初始化计算参数包括:
破片在目标坐标系内的初始位置PFInit_TargCrd(n0)为:
弹目相对速度Vrel可用以下公式计算:
其中[VXM_LaunCrd,VYM_LaunCrd,VZM_LaunCrd]为发射坐标系内的导弹速度,[VXT_LaunCrd,VYT_LaunCrd,VZT_LaunCrd]为发射坐标系内的目标速度;
相对速度坐标系内的破片初速VFInit_VrelCrd(n0)可用以下公式计算:
弹体系内的破片初速VFInit_TargCrd(n0)用以下公式计算:
发射坐标系内的破片初速VFInit_LaunCrd(n0)用以下公式计算:
每条轨迹的起始速度与位置计算公式如下:
VFTargCrd(n0,0)=VFInit_TargCrd(n0)
VFLaunCrd(n0,0)=VFInit_LaunCrd(n0)
PFTargCrd(n0,0)=PFInit_TargCrd(n0)
其中PFInit_TargCrd(n0)为破片在目标坐标系内的初始位置,由前文计算得到。
其中,第二步、根据命中检测算法,求解每段破片飞行轨迹是否命中目标,并进行下一步迭代计算包括:
使用命中检测的算法,求解轨迹的长度、末端的位置与破片存速,如果命中目标,计算出此段轨迹命中目标的时间thtTarg,t(n0,n1)=thtTarg,此线段轨迹终止计算;如果未命中目标,则t(n0,n1)=Δtatn,破片飞行距离R(n0,n1)=t(n0,n1)|VFLaunCrd(n0,n1-1)|,发射坐标系内破片存速:其中为上一计算周期内发射坐标系内破片存速,k为衰减系数;
目标坐标系内的破片存速为:
此破片轨迹继续计算至命中目标或者超过设定的破片最大飞行时间。其中EV2L为相对速度坐标系到发射坐标系的转换矩阵,ET2V为目标坐标系到相对速度坐标系的转转矩阵,[VXT_LaunCrd,VYT_LaunCrd,VZT_LaunCrd]为发射坐标系内的目标速度。
本发明的优点是考虑了由于空气阻力引起的战斗部破片飞散速度的衰减,弹目相对速度坐标内的破片飞散轨迹不是直线而是曲线,该方法能够精确计算出破片命中目标的位置,并可针对大规模仿真计算进行精度设置,满足不同计算精度和计算规模的平衡。相比传统的战斗部破片飞散轨迹计算的方法,该方法能够准确仿真评估超低空作战条件下的引战配合设计、引战系统设计和武器系统效能。
具体实施方式
以下对本发明的具体实施方式做出详细说明。
本发明提出一种超低空作战条件下的战斗部破片飞散轨迹计算的方法的具体步骤为:
第一步、根据引信起爆位置,初始化计算参数。
通过起爆点时刻的弹道参数,初始化以下参数:发射坐标系内的导弹速度VM_LaunCrd=[VXM_LaunCrd,VYM_LaunCrd,VZM_LaunCrd]、发射坐标系内的目标速度VT_LaunCrd=[VXT_LaunCrd,VYT_LaunCrd,VZT_LaunCrd]、起爆时刻的导弹高度HM、相对速度坐标系到发射坐标系的转换矩阵EV2L、相对速度坐标系到弹体系的转换矩阵EV2M、目标坐标系到相对速度坐标系的转转矩阵ET2V、破片初始位置PInit(n0)与速度VInit(n0)计算,相对速度坐标系内的破片初始位置PFInit_VrelCrd(n0)、弹目相对速度Vrel,弹体内破片初速VFInit_MslCrd(n0)、相对速度坐标系内的破片初速VFInit_VrelCrd(n0)、弹体系内的破片初速VFInit_TargCrd(n0)和发射坐标系内的破片初速VFInit_LaunCrd(n0)。通过仿真参数设置破片飞行时间T和破片飞行的衰减时间间隔设置Δtatn。
第二步、根据命中检测算法,求解每段破片飞行轨迹是否命中目标,并进行下一步迭代计算。
采用迹线法求解每段破片飞行轨迹是否命中目标、每一段破片轨迹的长度、轨迹末端的位置与破片存速,如果命中目标,计算出此段轨迹命中目标的时间thtTarg,t(n0,n1)=thtTarg,此线段轨迹终止计算;如果未命中目标,则t(n0,n1)=Δtatn,破片飞行距离为R(n0,n1),发射坐标系内破片存速VFLaunCrd(n0,n1),目标坐标系内的破片存速为VFTargCrd(n0,n1),此破片轨迹继续计算至命中目标或者超过设定的最大破片飞行时间。
具体的,本发明一种超低空作战条件下的战斗部破片飞散轨迹计算的方法的具体步骤为:
第一步、根据引信起爆位置,初始化计算参数:
破片在目标坐标系内的初始位置PFInit_TargCrd(n0)为:
弹目相对速度Vrel可用以下公式计算:
其中[VXM_LaunCrd,VYM_LaunCrd,VZM_LaunCrd]为发射坐标系内的导弹速度,[VXT_LaunCrd,VYT_LaunCrd,VZT_LaunCrd]为发射坐标系内的目标速度;
相对速度坐标系内的破片初速VFInit_VrelCrd(n0)可用以下公式计算:
弹体系内的破片初速VFInit_TargCrd(n0)用以下公式计算:
发射坐标系内的破片初速VFInit_LaunCrd(n0)用以下公式计算:
每条轨迹的起始速度与位置计算公式如下:
VFTargCrd(n0,0)=VFInit_TargCrd(n0)
VFLaunCrd(n0,0)=VFInit_LaunCrd(n0)
PFTargCrd(n0,0)=PFInit_TargCrd(n0)
其中PFInit_TargCrd(n0)为破片在目标坐标系内的初始位置,由前文计算得到;
第二步、根据命中检测算法,求解每段破片飞行轨迹是否命中目标,并进行下一步迭代计算:
使用命中检测的算法,求解轨迹的长度、末端的位置与破片存速,如果命中目标,计算出此段轨迹命中目标的时间thtTarg,t(n0,n1)=thtTarg,此线段轨迹终止计算;如果未命中目标,则t(n0,n1)=Δtatn,破片飞行距离R(n0,n1)=t(n0,n1)|VFLaunCrd(n0,n1-1)|,发射坐标系内破片存速:其中为上一计算周期内发射坐标系内破片存速,k为衰减系数;
目标坐标系内的破片存速为:
此破片轨迹继续计算至命中目标或者超过设定的破片最大飞行时间。其中EV2L为相对速度坐标系到发射坐标系的转换矩阵,ET2V为目标坐标系到相对速度坐标系的转转矩阵,[VXT_LaunCrd,VYT_LaunCrd,VZT_LaunCrd]为发射坐标系内的目标速度。
发明考虑了由于空气阻力引起的战斗部破片飞散速度的衰减,弹目相对速度坐标内的破片飞散轨迹不是直线而是曲线,该方法能够精确计算出破片命中目标的位置,并可针对大规模仿真计算进行精度设置,满足不同计算精度和计算规模的平衡。相比传统的战斗部破片飞散轨迹计算的方法,该方法能够准确仿真评估超低空作战条件下的引战配合设计、引战系统设计和武器系统效能。
Claims (3)
1.一种超低空作战条件下的战斗部破片飞散轨迹计算的方法,其特征在于,其步骤为:
第一步、根据引信起爆位置,初始化计算参数;
第二步、根据命中检测算法,求解每段破片飞行轨迹是否命中目标,并进行下一步迭代计算。
2.根据权利要求1所述的超低空作战条件下的战斗部破片飞散轨迹计算的方法,其特征在于,第一步、根据引信起爆位置,初始化计算参数:
破片在目标坐标系内的初始位置PFInit_TargCrd(n0)为:
弹目相对速度Vrel用以下公式计算:
其中[VXM_LaunCrd,VYM_LaunCrd,VZM_LaunCrd]为发射坐标系内的导弹速度,[VXT_LaunCrd,VYT_LaunCrd,VZT_LaunCrd]为发射坐标系内的目标速度;
相对速度坐标系内的破片初速VFInit_VrelCrd(n0)用以下公式计算:
弹体系内的破片初速VFInit_TargCrd(n0)用以下公式计算:
发射坐标系内的破片初速VFInit_LaunCrd(n0)用以下公式计算:
每条轨迹的起始速度与位置计算公式如下:
VFTargCrd(n0,0)=VFInit_TargCrd(n0)
VFLaunCrd(n0,0)=VFInit_LaunCrd(n0)
PFTargCrd(n0,0)=PFInit_TargCrd(n0)
其中PFInit_TargCrd(n0)为破片在目标坐标系内的初始位置,由前文计算得到。
3.根据权利要求1所述的超低空作战条件下的战斗部破片飞散轨迹计算的方法,其特征在于,
第二步、根据命中检测算法,求解每段破片飞行轨迹是否命中目标,并进行下一步迭代计算:
使用命中检测的算法,求解轨迹的长度、末端的位置与破片存速,如果命中目标,计算出此段轨迹命中目标的时间thtTarg,t(n0,n1)=thtTarg,此线段轨迹终止计算;如果未命中目标,则t(n0,n1)=△tatn,破片飞行距离R(n0,n1)=t(n0,n1)|VFLaunCrd(n0,n1-1)|,发射坐标系内破片存速:其中为上一计算周期内发射坐标系内破片存速,k为衰减系数;
目标坐标系内的破片存速为:
此破片轨迹继续计算至命中目标或者超过设定的破片最大飞行时间,其中EV2L为相对速度坐标系到发射坐标系的转换矩阵,ET2V为目标坐标系到相对速度坐标系的转转矩阵,[VXT_LaunCrd,VYT_LaunCrd,VZT_LaunCrd]为发射坐标系内的目标速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911194809.1A CN111125623B (zh) | 2019-11-28 | 2019-11-28 | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911194809.1A CN111125623B (zh) | 2019-11-28 | 2019-11-28 | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111125623A true CN111125623A (zh) | 2020-05-08 |
CN111125623B CN111125623B (zh) | 2023-09-19 |
Family
ID=70496975
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911194809.1A Active CN111125623B (zh) | 2019-11-28 | 2019-11-28 | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111125623B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1482278A (en) * | 1973-11-07 | 1977-08-10 | Dassault Electronique | Method and installation for firing a succession of projectiles |
WO2013105093A1 (en) * | 2012-01-10 | 2013-07-18 | Israel Aerospace Industrie Ltd. | Anti-rocket system |
CN107798208A (zh) * | 2016-08-28 | 2018-03-13 | 南京理工大学 | 对空目标导弹破片飞散最大毁伤算法 |
CN108563893A (zh) * | 2018-04-24 | 2018-09-21 | 南京理工大学 | 杀爆类战斗部破片流密度计算方法 |
CN108694281A (zh) * | 2018-05-14 | 2018-10-23 | 西安工业大学 | 一种基于多维角度破片散布侵彻的目标毁伤概率计算方法 |
CN108984808A (zh) * | 2017-06-04 | 2018-12-11 | 南京理工大学 | 引战配合仿真交会破片命中点参数计算方法 |
-
2019
- 2019-11-28 CN CN201911194809.1A patent/CN111125623B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1482278A (en) * | 1973-11-07 | 1977-08-10 | Dassault Electronique | Method and installation for firing a succession of projectiles |
WO2013105093A1 (en) * | 2012-01-10 | 2013-07-18 | Israel Aerospace Industrie Ltd. | Anti-rocket system |
CN107798208A (zh) * | 2016-08-28 | 2018-03-13 | 南京理工大学 | 对空目标导弹破片飞散最大毁伤算法 |
CN108984808A (zh) * | 2017-06-04 | 2018-12-11 | 南京理工大学 | 引战配合仿真交会破片命中点参数计算方法 |
CN108563893A (zh) * | 2018-04-24 | 2018-09-21 | 南京理工大学 | 杀爆类战斗部破片流密度计算方法 |
CN108694281A (zh) * | 2018-05-14 | 2018-10-23 | 西安工业大学 | 一种基于多维角度破片散布侵彻的目标毁伤概率计算方法 |
Non-Patent Citations (6)
Title |
---|
YIN GUOXIN,ZHAO HUI,FENG QIAOLING,WANG XIAOYAN,WANG GAO: ""Velocity measurement of warhead fragments of parallel net target"", 《PROCEEDINGS OF 2011 INTERNATIONAL CONFERENCE ON ELECTRONICS AND OPTOELECTRONICS》 * |
YIN GUOXIN,ZHAO HUI,FENG QIAOLING,WANG XIAOYAN,WANG GAO: ""Velocity measurement of warhead fragments of parallel net target"", 《PROCEEDINGS OF 2011 INTERNATIONAL CONFERENCE ON ELECTRONICS AND OPTOELECTRONICS》, 8 September 2011 (2011-09-08) * |
徐新林等: "舰空导弹引战配合可视化仿真分析方法", 《现代防御技术》 * |
徐新林等: "舰空导弹引战配合可视化仿真分析方法", 《现代防御技术》, no. 02, 15 April 2011 (2011-04-15) * |
李向东等: "破片式反导导弹引战配合仿真与效率计算", 《上海航天》 * |
李向东等: "破片式反导导弹引战配合仿真与效率计算", 《上海航天》, no. 03, 25 June 2006 (2006-06-25) * |
Also Published As
Publication number | Publication date |
---|---|
CN111125623B (zh) | 2023-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106407596B (zh) | 防空导弹命中毁伤过程建模仿真方法 | |
CN108920743A (zh) | 一种基于炸点空间位置的目标毁伤计算方法 | |
CN113359813B (zh) | 一种多领弹多从弹群组协同制导方法 | |
US11867484B2 (en) | Apparatus for generating a vitals target | |
CN108201695A (zh) | 高环境仿真射击游戏系统 | |
CN115238450A (zh) | 一种制导火箭弹行进间发射攻击区仿真计算方法 | |
CN112464451B (zh) | 基于作战仿真系统的防空导弹武器命中概率修正方法 | |
CN111125623A (zh) | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 | |
CN104713428A (zh) | 定向引战全系统的地面联合试验装置及其试验方法 | |
CN105841550A (zh) | 一种具有高度约束的高置修正比例导引律方法 | |
CN202141380U (zh) | 一种自动瞄准二次加速矩阵炮 | |
RU2707637C1 (ru) | Способ поражения воздушной цели управляемой ракетой | |
CN106643298A (zh) | 一种基于预置命中点的大气层内反导拦截器中制导方法 | |
CN112861330B (zh) | 一种基于matlab的导弹杀伤效果计算和可视化方法 | |
Lim | Predicting the accuracy of unguided artillery projectiles | |
RU2819146C2 (ru) | Способ наведения неуправляемых метаемых снарядов на цель | |
CN111023907A (zh) | 一种具有伤情模拟的智能靶标设备及伤情判定方法 | |
CN114970201A (zh) | 动对动射击条件下坦克炮首发命中概率的确定方法 | |
CN112444157A (zh) | 一种直瞄面杀伤武器弹药的射击模拟方法 | |
RU2744208C1 (ru) | Способ расчета индивидуальной функции сопротивления воздуха неуправляемого артиллерийского снаряда по результатам табличных стрельб по местности | |
Baranowski et al. | The analysis of the 35 mm artillery projectile's motion model parameters' identification based on the recorded flight trajectory | |
US20220357129A1 (en) | Method, computer program and weapons system for calculating a bursting point of a projectile | |
Li et al. | Achieving Interactive Shooting Training: Mathematical Modeling of Gun-launched Missile Ballistic Trajectories in LVC Architecture | |
CN115823972A (zh) | 一种导弹引信炸点控制方法 | |
Totev et al. | Effects of Contactless Suppressors on Accuracy in Shooting |
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 |