[發明專利]一種超低空作戰條件下的戰斗部破片飛散軌跡計算的方法

有效

• 申請號: 201911194809.1
• 申請日: 2019-11-28
• 公開（公告）號: CN111125623B
• 公開（公告）日: 2023-09-19
• 發明（設計）人: 趙琦;苗立功;孔凡勛;袁晉;李文平
• 申請（專利權）人: 北京電子工程總體研究所
• 主分類號: G06F17/16
• 分類號: G06F17/16
• 代理公司: 中國航天科工集團公司專利中心 11024
• 代理人: 葛鵬

1.一種超低空作戰條件下的戰斗部破片飛散軌跡計算的方法，其特征在于，其步驟為：

第一步、根據相對速度坐標系內的破片初始位置破片在彈體系內的初始位置和導彈在相對速度坐標系內的起爆位置為得到：

破片在目標坐標系內的初始位置PF_{Init_TargCrd}(n₀)為：

彈目相對速度V_rel用以下公式計算：

其中[VX_{M_LaunCrd},VY_{M_LaunCrd},VZ_{M_LaunCrd}]為發射坐標系內的導彈速度，[VX_{T_LaunCrd},VY_{T_LaunCrd},VZ_{T_LaunCrd}]為發射坐標系內的目標速度；

相對速度坐標系內的破片初速VF_{Init_VrelCrd}(n₀)用以下公式計算：

其中VF_{Init_MslCrd}(n₀)為彈體內破片初速，[VX_rel,VY_rel,VZ_rel]為彈目相對速度，E_V2M為相對速度坐標系到彈體系的轉換矩陣，為彈體內破片初速；

彈體系內的破片初速VF_{Init_TargCrd}(n₀)用以下公式計算：

其中E_T2V為目標坐標系到相對速度坐標系的轉轉矩陣，為相對速度坐標系內的破片初速；

發射坐標系內的破片初速VF_{Init_LaunCrd}(n₀)用以下公式計算：

其中E_V2L為相對速度坐標系到發射坐標系的轉換矩陣，為發射坐標系內的導彈速度；

破片飛行時間時間常數a∈{1.0,2.0}，每一枚破片所需計算的循環次數N₁＝ceil(T/Δt_atn)，Δt_atn為設定的時間步長，PM_{Init_VrelCrd}為導彈在相對速度坐標系內的起爆位置；

每條軌跡的起始速度與位置計算公式如下：

VF_TargCrd(n₀,0)＝VF_{Init_TargCrd}(n₀)

VF_LaunCrd(n₀,0)＝VF_{Init_LaunCrd}(n₀)

PF_TargCrd(n₀,0)＝PF_{Init_TargCrd}(n₀)

其中PF_{Init_TargCrd}(n₀)為破片在目標坐標系內的初始位置，由前文計算得到；

第二步、根據命中檢測算法，求解每段破片飛行軌跡是否命中目標，并進行下一步迭代計算。

2.根據權利要求1所述的超低空作戰條件下的戰斗部破片飛散軌跡計算的方法，其特征在于，

第二步、根據命中檢測算法，求解每段破片飛行軌跡是否命中目標，并進行下一步迭代計算：

使用跡線法，求解軌跡的長度、末端的位置與破片存速，如果命中目標，計算出此段軌跡命中目標的時間t_htTarg，t(n₀,n₁)＝t_htTarg，此線段軌跡終止計算；如果未命中目標，則t(n₀,n₁)＝Δt_atn，破片飛行距離R(n₀,n₁)＝t(n₀,n₁)|VF_LaunCrd(n₀,n₁-1)|，發射坐標系內破片存速：其中為上一計算周期內發射坐標系內破片存速，k為衰減系數；

目標坐標系內的破片存速為：

此破片軌跡繼續計算至命中目標或者超過設定的破片最大飛行時間，其中E_V2L為相對速度坐標系到發射坐標系的轉換矩陣，E_T2V為目標坐標系到相對速度坐標系的轉轉矩陣，[VX_{T_LaunCrd},VY_{T_LaunCrd},VZ_{T_LaunCrd}]為發射坐標系內的目標速度。