留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

同步扫描周视脉冲激光引信多探测点最佳起爆建模及仿真

张合 李红霞 丁立波 查冰婷

张合, 李红霞, 丁立波, 查冰婷. 同步扫描周视脉冲激光引信多探测点最佳起爆建模及仿真[J]. 红外与激光工程, 2020, 49(4): 0403001-0403001-7. doi: 10.3788.IRLA202049.0403001
引用本文: 张合, 李红霞, 丁立波, 查冰婷. 同步扫描周视脉冲激光引信多探测点最佳起爆建模及仿真[J]. 红外与激光工程, 2020, 49(4): 0403001-0403001-7. doi: 10.3788.IRLA202049.0403001
Zhang He, Li Hongxia, Ding Libo, Zha Bingting. Modeling and simulation of multi-detection point optimal initiation of synchronous scanning panoramic pulse laser fuze[J]. Infrared and Laser Engineering, 2020, 49(4): 0403001-0403001-7. doi: 10.3788.IRLA202049.0403001
Citation: Zhang He, Li Hongxia, Ding Libo, Zha Bingting. Modeling and simulation of multi-detection point optimal initiation of synchronous scanning panoramic pulse laser fuze[J]. Infrared and Laser Engineering, 2020, 49(4): 0403001-0403001-7. doi: 10.3788.IRLA202049.0403001

同步扫描周视脉冲激光引信多探测点最佳起爆建模及仿真

doi: 10.3788.IRLA202049.0403001
基金项目: 国家自然科学基金(51709147)
详细信息
    作者简介:

    张合(1957-),男,教授,博士生导师,主要研究方向为目标探测与识别。Email:hezhangz@njust.edu.cn

    通讯作者: 丁立波(1977-),男,副教授,硕士生导师,主要研究方向为近程探测与控制。Email: ding_libo@hotmail.com
  • 中图分类号: TJ43

Modeling and simulation of multi-detection point optimal initiation of synchronous scanning panoramic pulse laser fuze

  • 摘要:

    周视激光引信通常仅根据单次测量的目标方位和距离解算最佳起爆时间和起爆方位角,精度不高,另一方面,引信无法判别目标真伪,易发生“虚警”或误炸。基于同步扫描周视脉冲激光引信工作原理,分析了弹目交会过程中多个探测点的方位、距离及其位置坐标,建立了基于多探测点的目标速度模型,对目标速度进行计算以辨别目标真伪。解算了目标特征几何中心位置,建立了最佳起爆时间和方位角计算模型,使战斗部获得最大毁伤效能,并分析了最佳起爆时间和方位角的影响因素。研究可为周视激光引信目标真伪辨别和定向起爆控制提供参考。

  • 图  1  脉冲激光同步扫描探测方案

    Figure  1.  Pulse laser synchronous scanning detection scheme

    图  2  目标探测点${F_{{k}}}$在弹体坐标系${o_m}{x_m}{y_m}{z_m}$中的位置示意图

    Figure  2.  Location diagram of target detection point ${F_{{k}}}$ in missile coordinate system ${o_m}{x_m}{y_m}{z_m}$

    图  3  大地坐标系${o_0}{x_0}{y_0}{z_0}$中系统探测目标示意图

    Figure  3.  System detecting target in geodetic coordinate system ${o_0}{x_0}{y_0}{z_0}$

    图  4  目标表面起始探测点与终止探测点示意图

    Figure  4.  Diagram of start detection point and end detection point on target surface

    图  5  破片与目标交会模型。(a)大地坐标系中破片与目标交会过程;(b)弹体坐标系中起爆方位角

    Figure  5.  Fragment and target intersection model. (a) Intersection process of fragments and targets in the geodetic coordinate system; (b) Initiation azimuth in the projectile coordinate system

    图  6  最佳起爆时间和方位角与破片飞散速度变化关系。(a)最佳起爆时间与破片飞散速度关系图;(b)最佳起爆方位角与破片飞散速度关系图

    Figure  6.  Relationship between the optimal initiation time, initiation azimuth and the fragment speed. (a)Optimal initiation time vs fragment speed; (b)Optimal initiation azimuth vs fragment speed

    图  7  最佳起爆时间和方位角与弹丸速度变化关系。(a)最佳起爆时间与弹丸速度关系图;(b)最佳起爆方位角与弹丸速度关系图

    Figure  7.  Relationship between the optimal initiation time, initiation azimuth and the projectile speed. (a)Optimal initiation time vs projectile speed; (b)Optimal initiation azimuth vs projectile speed

    图  8  最佳起爆时间和方位角与弹目距离变化关系。(a)最佳起爆时间与弹目距离关系图;(b)最佳起爆方位角与弹目距离关系图

    Figure  8.  Relationship between the optimal initiation time, initiation azimuth and the distance between the projectile with the target.(a)Optimal initiation time vs the distance between the projectile and target;(b)Optimal initiation azimuth vs the distance between projectile and target

    图  9  最佳起爆时间和方位角与目标特征长度变化关系。(a)最佳起爆时间与目标特征长度关系图;(b)最佳起爆方位角与目标特征长度关系图

    Figure  9.  Relationship between the optimal initiation time, initiation azimuth and the target characteristic length.(a)Optimal initiation time vs the target characteristic length;(b)Optimal initiation azimuth vs the target characteristic length

    图  10  最佳起爆时间和方位角与目标速度变化关系。(a)最佳起爆时间与目标速度关系图;(b)最佳起爆方位角与目标速度关系图

    Figure  10.  Relationship between the optimal initiation time, initiation azimuth and the target speed. (a)Optimal initiation time vs target speed; (b)Optimal initiation azimuth vs target speed

  • [1] 查冰婷, 张合. 单光束脉冲激光引信扫描频率和脉冲频率研究[J]. 红外与激光工程, 2014, 43(7): 2081−2086. doi:  10.3969/j.issn.1007-2276.2014.07.007

    Zha Bingting, Zhang He. Scanning frequency and pulse frequency of single-beam pulsed laser fuze [J]. Infrared and Laser Engineering, 2014, 43(7): 2081−2086. (in Chinese) doi:  10.3969/j.issn.1007-2276.2014.07.007
    [2] 孔德浩, 路明, 苏益德. 多光束周视激光引信与定向战斗部的配合[J]. 太赫兹科学与电子信息学报, 2019, 17(3): 402−407.

    Kong Dehao, Lu Ming, Su Yide. Muli-bean circum-directional laser fuze and aimed warhead [J]. Journal of Terahertz Science and Electronic Information Technology, 2019, 17(3): 402−407. (in Chinese)
    [3] 郭泽荣, 李洛, 温玉全, 等. 六光束脉冲激光探测定向战斗部最佳起爆方位角研究[J]. 北京理工大学学报, 2015, 35(4): 347−350, 356.

    Guo Zerong, Li Luo, Wen Yuquan, et al. Study on the optimum detonating azimuth of six laser range-finder for aimed wared [J]. Journal of Beijing University of Science and Technology, 2015, 35(4): 347−350, 356. (in Chinese)
    [4] 石志彬, 高敏, 杨锁昌, 等. 瞄准式战斗部威力半径与起爆延时研究[J]. 现代防御技术, 2014, 42(3): 1−5, 30. doi:  10.3969/j.issn.1009-086x.2014.03.001

    Shi Zhibin, Gao Min, Yang Suochang, et al. Lethal radius and initiation delay of gimbaled warhead [J]. Modern Defense Technology, 2014, 42(3): 1−5, 30. (in Chinese) doi:  10.3969/j.issn.1009-086x.2014.03.001
    [5] 王马法, 卢芳云, 李翔宇. 战斗部最佳起爆延时的计算模型及其应用研究[J]. 弹道学报, 2011, 23(4): 95−99.

    Wang Mafa, Lu Fangyun, Li Xiangyu. Calculation model of optimum delay time of detonation of warhead and its application [J]. Journal of Ballistics, 2011, 23(4): 95−99. (in Chinese)
    [6] Steinvall O. Effects of target shape and reflection on laser radar cross sections [J]. Applied Optics, 2000, 39(24): 4381−4391. doi:  10.1364/AO.39.004381
    [7] Nyongesa H O, Kent S. Genetic programming for anti-air missile proximity fuze delay-time algorithms [J]. IEEE Aerospace and Electronic Systems Magazine, 2001, 16(1): 41−45. doi:  10.1109/62.894177
    [8] 刘博虎, 宋承天, 段亚博. 周视激光成像引信设计及仿真[J]. 红外与激光工程, 2018, 47(4): 0406003.

    Liu Bohu, Song Chengtian, Duan Yabo. Design and simulation of a panoramic laser imaging fuze [J]. Infrared and Laser Engineering, 2018, 47(4): 0406003. (in Chinese)
    [9] 查冰婷, 袁海璐, 马少杰, 等. 单光束扩束扫描激光周视探测系统参数对探测能力的影响[J]. 物理学报, 2019, 68(7): 79−88.

    Zha Bingting, Yuan Hailu, Ma Shaojie, et al. Influence of single-beam expanding scanning laser circumferential detection system parameters on detection capability [J]. Journal of Physics, 2019, 68(7): 79−88. (in Chinese)
    [10] 秦帅, 杨锁昌, 李乐, 等. 瞄准式战斗部最佳起爆延时控制研究[J]. 弹箭与制导学报, 2012, 32(4): 92−94. doi:  10.3969/j.issn.1673-9728.2012.04.023

    Qin Shuai, Yang Suochang, Li Le, et al. The research on optimal detonation delay control of aimed warhead [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2012, 32(4): 92−94. (in Chinese) doi:  10.3969/j.issn.1673-9728.2012.04.023
  • [1] 任维贺, 张月, 苏云, 张学敏, 邓红艳, 柳祎.  环境扰动下空中动目标探测技术综述 . 红外与激光工程, 2022, 51(9): 20210843-1-20210843-18. doi: 10.3788/IRLA20210843
    [2] 洪汉玉, 吴世康, 时愈, 吴锦梦, 孙春生.  非合作水雷目标图像非均匀强噪声去除方法 . 红外与激光工程, 2021, 50(3): 20200344-1-20200344-10. doi: 10.3788/IRLA20200344
    [3] 苏宏, 侯亚丽, 王刚, 李岗.  激光引信12分元探测定向起爆控制建模及仿真 . 红外与激光工程, 2020, 49(4): 0403003-0403003-7. doi: 10.3788/IRLA202049.0403003
    [4] 张俊, 张洪健, 孙大开, 王立, 武延鹏, 李春艳, 钟红军, 卢欣.  高灵敏度空间目标自主探测技术研究 . 红外与激光工程, 2020, 49(5): 20201008-20201008-7. doi: 10.3788/IRLA20201008
    [5] 梁天全, 张晓云, 段朋, 于会山, 张保华, 汤庆新.  改进暗通道方法的水下强散射介质目标探测 . 红外与激光工程, 2020, 49(2): 0203012-0203012. doi: 10.3788/IRLA202049.0203012
    [6] 张伟, 李玉钊, 汪哲思, 刘岩.  周视激光引信超低空探测地面回波特性研究 . 红外与激光工程, 2020, 49(4): 0403008-0403008-8. doi: 10.3788/IRLA202049.0413008
    [7] 李佳昕, 白廷柱, 崔志刚, 宋翠芬, 成泽明.  像素偏振成像系统的耦合实现 . 红外与激光工程, 2019, 48(9): 925001-0925001(7). doi: 10.3788/IRLA201948.0925001
    [8] 彭波, 钟昆, 赵慧, 李中云.  水下目标激光周向扫描探测模型与仿真分析 . 红外与激光工程, 2019, 48(12): 1205002-1205002(7). doi: 10.3788/IRLA201948.1205002
    [9] 谢绍禹, 叶茂, 赵毅强, 吴环宝, 贾晓东.  偏振探测激光引信目标判别技术 . 红外与激光工程, 2018, 47(S1): 5-9. doi: 10.3788/IRLA201847.S106002
    [10] 刘博虎, 宋承天, 段亚博.  周视激光成像引信设计及仿真 . 红外与激光工程, 2018, 47(4): 406003-0406003(8). doi: 10.3788/IRLA201847.0406003
    [11] 于玉亭, 舒敬荣, 丁伯圣.  激光引信与预制破片战斗部引战配合设计 . 红外与激光工程, 2018, 47(3): 303003-0303003(6). doi: 10.3788/IRLA201847.0303003
    [12] 徐超, 何利民, 王霞, 金伟其.  红外偏振成像系统高速处理模块设计 . 红外与激光工程, 2017, 46(2): 204002-0204002(8). doi: 10.3788/IRLA201746.0204002
    [13] 王霞, 梁建安, 龙华宝, 姚锦华, 夏润秋, 贺思, 金伟其.  典型背景和目标的长波红外偏振成像实验研究 . 红外与激光工程, 2016, 45(7): 704002-0704002(7). doi: 10.3788/IRLA201645.0704002
    [14] 程纯, 王克逸, 曹兆楼, 马孟超.  曲面复眼中对数型锥透镜成像特性研究 . 红外与激光工程, 2015, 44(4): 1284-1289.
    [15] 黄富瑜, 沈学举, 何永强, 周冰.  超大视场成像系统对空间目标的探测能力分析 . 红外与激光工程, 2015, 44(10): 3134-3140.
    [16] 方义强, 陈卫, 孙晓军, 马东辉, 程正东, 张发强.  云对空中目标红外探测的影响 . 红外与激光工程, 2015, 44(7): 2034-2038.
    [17] 许俊峰, 姜春兰, 毛亮, 王在成, 李明.  测距成像一体化引信与可瞄准战斗部配合技术 . 红外与激光工程, 2014, 43(6): 1794-1800.
    [18] 王霞, 夏润秋, 金伟其, 刘敬, 梁建安.  红外偏振成像探测技术进展 . 红外与激光工程, 2014, 43(10): 3175-3182.
    [19] 田昌会, 杨百愚, 蔡明, 范琦, 王斌科, 王伟宇, 屈绍波.  大气背景对红外目标探测的影响 . 红外与激光工程, 2014, 43(2): 438-441.
    [20] 甘霖, 张合, 张祥金, 冯颖.  激光近炸引信单光束脉冲周向探测技术 . 红外与激光工程, 2013, 42(1): 84-89.
  • 加载中
图(10)
计量
  • 文章访问数:  1214
  • HTML全文浏览量:  548
  • PDF下载量:  68
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-10
  • 修回日期:  2020-03-26
  • 刊出日期:  2020-04-01

同步扫描周视脉冲激光引信多探测点最佳起爆建模及仿真

doi: 10.3788.IRLA202049.0403001
    作者简介:

    张合(1957-),男,教授,博士生导师,主要研究方向为目标探测与识别。Email:hezhangz@njust.edu.cn

    通讯作者: 丁立波(1977-),男,副教授,硕士生导师,主要研究方向为近程探测与控制。Email: ding_libo@hotmail.com
基金项目:  国家自然科学基金(51709147)
  • 中图分类号: TJ43

摘要: 

周视激光引信通常仅根据单次测量的目标方位和距离解算最佳起爆时间和起爆方位角,精度不高,另一方面,引信无法判别目标真伪,易发生“虚警”或误炸。基于同步扫描周视脉冲激光引信工作原理,分析了弹目交会过程中多个探测点的方位、距离及其位置坐标,建立了基于多探测点的目标速度模型,对目标速度进行计算以辨别目标真伪。解算了目标特征几何中心位置,建立了最佳起爆时间和方位角计算模型,使战斗部获得最大毁伤效能,并分析了最佳起爆时间和方位角的影响因素。研究可为周视激光引信目标真伪辨别和定向起爆控制提供参考。

English Abstract

    • 激光探测由于精度高、方向性好、抗电磁干扰能力强,常用于近炸引信实现高精度定距和方位识别[1-2]。周视激光引信能探测360°范围内的目标,测量目标距离、识别目标与引信的相对方位。若将周视激光引信配装于定向战斗部,可在最佳时机向最佳方位起爆战斗部,实现目标的高效毁伤,其中起爆时机和方位的解算结果将直接影响战斗部的毁伤效果[3-8]。与此同时,由于激光探测本身不具备识别目标特性的能力,激光引信探测过程中易受干扰目标的影响,发生“虚警”或误炸。

      文中针对脉冲激光同步扫描体制的周视激光引信,基于在探测过程中获取的多个目标探测点的距离和方位信息,建立目标速度模型,通过对目标速度进行估算,以辨别目标真伪。同时综合目标上多探测点信息,预估目标特征几何中心位置,解算引信最佳起爆时间和方位角,使战斗部获得最大毁伤效能。

    • 脉冲激光同步扫描系统通过单根脉冲激光束360°旋转扫描探测周向目标,随扫描装置的旋转,激光束将形成如图1所示的圆锥状光束场。同时激光接收系统随扫描装置与发射光束同步旋转,接收来自目标反射的回波信号[9]。在探测目标的过程中,系统将逐点扫描目标表面的多个位置,形成探测点,根据各探测点反射的光脉冲信号,可获得各点的距离和方位角。

      图  1  脉冲激光同步扫描探测方案

      Figure 1.  Pulse laser synchronous scanning detection scheme

    • 以弹丸发射点${o_0}$为原点,建立大地坐标系${o_0}{x_0}{y_0}{z_0}$${x_0}$轴在水平面内沿弹丸发射方向,${y_0}$轴铅直向上,${z_0}$轴与${x_0}$轴、${y_0}$轴构成笛卡尔坐标系。以激光探测系统的发射中心${o_m}$为原点,建立弹体坐标系${o_m}{x_m}{y_m}{z_m}$${x_m}$轴沿弹丸轴线指向其前进方向,${y_m}$轴沿弹轴垂直方向指向上,${z_m}$轴与${x_m}$轴、${y_m}$轴构成笛卡尔坐标系。以目标的几何特征中心${o_t}$为原点,建立目标坐标系${o_t}{x_t}{y_t}{z_t}$${x_t}$轴沿目标轴线并指向其头部,${y_t}$轴沿目标轴线垂直方向,${z_t}$轴与${x_t}$轴、${y_t}$轴构成笛卡尔坐标系。

      若弹丸和目标的俯仰角为${\beta _m}$${\beta _t}$,偏航角为${\psi _m}$${\psi _t}$,所有角度以逆时针方向为正,可得三个坐标系之间的转换关系。弹体坐标系到大地坐标系的转换关系为:

      $$\left[ {\begin{aligned} {{x_0}} \\ {{y_0}} \\ {{z_0}} \end{aligned}} \right] = C_0^m\left[ {\begin{aligned} {{x_m}} \\ {{y_m}} \\ {{z_m}} \end{aligned}} \right]$$ (1)

      式中:$C_0^m$表示弹体坐标系到大地坐标系的方向余弦矩阵。

      $$C_0^m = \left[ {\begin{array}{*{20}{c}} {\cos {\beta _m}\cos {\psi _m}}&{ - \sin {\beta _m}\cos {\psi _m}}&{\sin {\psi _m}} \\ {\sin {\beta _m}}&{\cos {\beta _m}}&0 \\ { - \cos {\beta _m}\sin {\psi _m}}&{\sin {\beta _m}\sin {\psi _m}}&{\cos {\psi _m}} \end{array}} \right]$$ (2)

      大地坐标系与目标坐标系的转换关系为:

      $$\left[ {\begin{aligned} {{x_t}} \\ {{y_t}} \\ {{z_t}} \end{aligned}} \right] = C_t^0\left[ {\begin{aligned} {{x_0}} \\ {{y_0}} \\ {{z_0}} \end{aligned}} \right]$$ (3)

      式中:$C_t^0$表示从大地坐标系到目标坐标系的方向余弦矩阵。

      $$C_t^0 = \left[ {\begin{array}{*{20}{c}} {\cos {\beta _t}\cos {\psi _t}}&{\sin {\beta _t}}&{ - \cos {\beta _t}\sin {\psi _t}} \\ { - \sin {\beta _t}\cos {\psi _t}}&{\cos {\beta _t}}&{\sin {\beta _t}\sin {\psi _t}} \\ {\sin {\psi _t}}&0&{\cos {\psi _t}} \end{array}} \right]$$ (4)

      则弹体坐标系与目标坐标系之间的转换关系为:

      $$\left[ {\begin{aligned} {{x_t}} \\ {{y_t}} \\ {{z_t}} \end{aligned}} \right] = C_t^m\left[ {\begin{aligned} {{x_m}} \\ {{y_m}} \\ {{z_m}} \end{aligned}} \right] = C_t^0C_0^m\left[ {\begin{aligned} {{x_m}} \\ {{y_m}} \\ {{z_m}} \end{aligned}} \right]$$ (5)

      假设在周视激光引信探测过程中次序接收到目标表面${F_1}$${F_2}$、…、${F_k}$点反射的光脉冲信号,并解算出上述各探测点在弹体坐标系中的距离和方位角,从而能够获得各点在弹体坐标系中的三维坐标。图2表示探测点${F_k}$在弹体坐标系${o_m}{x_m}{y_m}{z_m}$中的位置。

      图  2  目标探测点${F_{{k}}}$在弹体坐标系${o_m}{x_m}{y_m}{z_m}$中的位置示意图

      Figure 2.  Location diagram of target detection point ${F_{{k}}}$ in missile coordinate system ${o_m}{x_m}{y_m}{z_m}$

      图2中,光束倾斜角${\alpha _t}$表示激光束与弹体坐标系${x_m}$的夹角,探测点${F_k}$到原点${o_m}$的距离为${R_k}$,方位角φ为原点${o_m}$和探测点${F_k}$的连线在${y_m}{o_m}{z_m}$平面上的投影与${y_m}$轴的夹角。根据投影定理和三角形法则,可知探测点${F_k}$在弹体坐标系中的坐标$ ({}^{{m}}{{{x}}_k},{}^{{m}}{{{y}}_k},{}^{{m}}{{{z}}_k})$为:

      $$\begin{split} ({}^{{m}}{{{x}}_k},{}^{{m}}{{{y}}_k},{}^{{m}}{{{z}}_k}) =\; & ({R_k}\cos {\alpha _t},{R_k}\sin {\alpha _t} \cos {\psi _k},\\ & {R_k}\sin {\alpha _t}\sin {\psi _k}) \end{split} $$ (6)
    • 在大地坐标系${o_0}{x_0}{y_0}{z_0}$中激光探测系统探测目标如图3所示。

      图  3  大地坐标系${o_0}{x_0}{y_0}{z_0}$中系统探测目标示意图

      Figure 3.  System detecting target in geodetic coordinate system ${o_0}{x_0}{y_0}{z_0}$

      图3中,定义激光坐标系${o_l}{x_l}{y_l}{z_l}$,坐标系原点${o_l}$${o_m}$重合,激光出射方向为${x_l}$轴。激光坐标系${o_l}{x_l}{y_l}{z_l}$可以看作先将弹体坐标系${o_m}{x_m}{y_m}{z_m}$${z_m}$轴逆时针旋转${\alpha _t}$,再绕${y_m}$轴旋转$\varphi $形成的,令$C_m^l$为激光坐标系${o_l}{x_l}{y_l}{z_l}$和弹体坐标系${o_m}{x_m}{y_m}{z_m}$之间的方向余弦矩阵,则有:

      $$C_{{m}}^l = \left[ {\begin{array}{*{20}{c}} {\cos \varphi \cos {\alpha _t}}&{ - \sin {\alpha _t}}&{ - \sin \varphi \cos {\alpha _t}} \\ {\cos \varphi \sin {\alpha _t}}&{\cos {\alpha _t}}&{ - \sin \varphi \sin {\alpha _t}} \\ {\sin \varphi }&0&{\cos \varphi } \end{array}} \right]$$ (7)

      故激光坐标系与弹体坐标系之间的转换为:

      $$\left[ {\begin{aligned} {{x_m}} \\ {{y_m}} \\ {{z_m}} \end{aligned}} \right] = C_{{m}}^{{l}}\left[ {\begin{aligned} {{x_l}} \\ {{y_l}} \\ {{z_l}} \end{aligned}} \right]$$ (8)

      联立公式(1),得激光坐标系与大地坐标系之间的转换关系为:

      $$\left[ {\begin{aligned} {{x_0}} \\ {{y_0}} \\ {{z_0}} \end{aligned}} \right] = C_0^{{m}}C_{{m}}^{{l}}\left[ {\begin{aligned} {{x_l}} \\ {{y_l}} \\ {{z_l}} \end{aligned}} \right]$$ (9)

      图3中,${{{r}}_t}$为探测点在大地坐标系上的位矢,${{{r}}_{{m}}}$为弹体坐标系和激光坐标系原点的位矢,${{{r}}_{mt}} = {{{r}}_t} - {{{r}}_m}$为探测点在激光坐标系上的位矢。根据速度的定义,可知在大地坐标系下目标的速度矢量${{{v}}_t}$、弹丸的速度矢量${{{v}}_{{m}}}$和目标与弹丸的相对速度矢量${{{v}}_{{{mt}}}}$分别为:

      $$ {{{v}}_{{t}}} = \frac{{{{\rm{d}}{{{r}}}_t}}}{{{\rm{d}}t}},{{{{v}}}_{{m}}} = \frac{{{\rm{d}}{{{{r}}}_m}}}{{{\rm{d}}t}},{{{{v}}}_{{{mt}}}} = \frac{{{\rm{d}}{{{{r}}}_{mt}}}}{{{\rm{d}}t}} $$ (10)

      ${{{{v}}}_{{{mt}}}}$在激光坐标系上投影,可得

      $${}^{{l}}{{{{v}}}_{{{mt}}}} = {(C_0^mC_m^l)^{\rm{T}}}({}^0{{{{v}}}_t} - {}^0{{{{v}}}_m})$$ (11)

      式中:左上标表示其所处的坐标系,其中

      $${}^{{l}}{{{v}}_{{{mt}}}} = \left\{ {\begin{aligned} {{}^{{l}}{v_{{{mtx}}}}} \\ {{}^{{l}}{v_{{{mty}}}}} \\ {{}^{{l}}{v_{{{mtz}}}}} \end{aligned}} \right\},{}^{{0}}{{{v}}_{{t}}} = \left\{ {\begin{aligned} {{}^{{0}}{v_{{{tx}}}}} \\ {{}^{{0}}{v_{{{ty}}}}} \\ {{}^{{0}}{v_{{{tz}}}}} \end{aligned}} \right\},{}^{{0}}{{{v}}_{{m}}} = \left\{ {\begin{aligned} {{}^{{0}}{v_{{{mx}}}}} \\ {{}^{{0}}{v_{{{my}}}}} \\ {{}^{{0}}{v_{{{mz}}}}} \end{aligned}} \right\}$$ (12)

      根据探测点${F_k}$的回波信号,可得点${F_k}$到弹体坐标系原点${{{o}}_{{m}}}$的距离${R_{{k}}}$。由多个探测点的距离${R_{{k}}}$,求解出${}^l{{{v}}_{{{mt}}}}$中的第一个元素${}^l{{{v}}_{{{mtx}}}}$,即在激光坐标系下弹丸与目标之间的相对速度${{{v}}_{mt}}$${x_l}$轴上的投影。由于采样频率较高,相邻探测点探测间隔时间短,可以计算相邻探测点的平均速度代替目标瞬时速度,已知$\left| {{{{r}}_{{\rm{mt}}}}} \right| = {R_k}$,可得:

      $${}^{{l}}{v_{{{mtx}}}} = \frac{{{R_{{{k}} + {{1}}}} - {R_k}}}{{\Delta t}}$$ (13)

      式中:$\Delta {{t}}$表示探测相邻点${F_{{k}}}$${F_{{{k}} + {{1}}}}$的时间间隔,一般情况下即为探测系统的脉冲周期${T_{{m}}}$。为提高目标运动速度的计算精度,通过探测目标上多个相邻点来计算${}^l{{{v}}_{{{mtx}}}}$的平均值。

      $$L = {(C_0^mC_m^l)^{\bf{T}}} = \left[ {\begin{aligned} {{l_{11}}}\quad {{l_{12}}} \quad {{l_{13}}} \\ {{l_{21}}}\quad {{l_{22}}} \quad {{l_{23}}} \\ {{l_{31}}}\quad {{l_{32}}} \quad {{l_{33}}} \end{aligned}} \right]$$ (14)

      联立公式(12)~(14),得:

      $$\left\{ {\begin{aligned} {{}^{{l}}{v_{{{mtx}}}}} \\ {{}^{{l}}{v_{{{mty}}}}} \\ {{}^{{l}}{v_{{{mtz}}}}} \end{aligned}} \right\} = \left[ {\begin{aligned} {{l_{11}}}\quad {{l_{12}}}\quad {{l_{13}}} \\ {{l_{21}}}\quad {{l_{22}}}\quad {{l_{23}}} \\ {{l_{31}}}\quad {{l_{32}}}\quad {{l_{33}}} \end{aligned}} \right]\left\{ {\begin{aligned} {{}^0{v_{tx}} - {}^0{v_{mx}}} \\ {{}^0{v_{ty}} - {}^0{v_{my}}} \\ {{}^0{v_{tz}} - {}^0{v_{mz}}} \end{aligned}} \right\}$$ (15)

      提取公式(15)中的第一个方程:

      $${l_{11}}{}^0{v_{tx}} + {l_{12}}{}^0{v_{ty}} + {l_{13}}{}^0{v_{tz}} = {}^lv{}_{mtx} + {l_{11}}{}^0{v_{mx}} + {l_{12}}{}^0{v_{my}} + {l_{13}}{}^0{v_{mz}}$$ (16)

      仅考虑目标做平面运动,其绝对速度矢量为$^0{{{v}}_t}$,而在公式(16)中包含$^0{v_{tx}}$$^0{v_{t{{y}}}}$$^0{v_{t{{z}}}}$三个未知量,可联立三个探测点信息的计算方程构成方程组:

      $$\left\{ {\begin{aligned} & {l_{11}^1{}^0{v_{tx}} + l_{12}^1{}^0{v_{ty}} + l_{13}^1{}^0{v_{tz}} = {}^lv_{{}_{mtx}}^1 + l_{11}^1{}^0{v_{mx}} + l_{12}^1{}^0{v_{my}} + l_{13}^1{}^0{v_{mz}}} \\ & {l_{11}^2{}^0{v_{tx}} + l_{12}^2{}^0{v_{ty}} + l_{13}^2{}^0{v_{tz}} = {}^lv_{{}_{mtx}}^2 + l_{11}^2{}^0{v_{mx}} + l_{12}^2{}^0{v_{my}} + l_{13}^2{}^0{v_{mz}}} \\ & {l_{11}^3{}^0{v_{tx}} + l_{12}^3{}^0{v_{ty}} + l_{13}^{31}{}^0{v_{tz}} = {}^lv_{{}_{mtx}}^3 + l_{11}^3{}^0{v_{mx}} + l_{12}^3{}^0{v_{my}} + l_{13}^3{}^0{v_{mz}}} \end{aligned}} \right.$$ (17)

      根据公式(17),通过三个探测点的回波信息即可求解目标运动速度的各个分量,从而计算出目标速度矢量$^0{{{v}}_t}$,采用多次测量求平均的方法可进一步提高测量精度。当弹药被用于攻击某一特定目标时,可预先在周视激光引信中装定攻击目标的速度范围,在激光探测目标过程中,通过目标上多探测点信息估算目标速度,并与既定目标速度范围进行比较,实现真伪目标的辨别,避免引信发生“虚警”或误炸。

    • 对于爆破战斗部,为实现最大毁伤效能,应使爆炸最大密度破片恰好攻击目标的易损部位[10]。通常简化认为目标几何中心为其易损部位,则要求引信能够在探测过程中尽可能准确地估计出目标的特征几何中心位置。

      在大地坐标系中,系统对目标的探测过程如图4所示。

      图  4  目标表面起始探测点与终止探测点示意图

      Figure 4.  Diagram of start detection point and end detection point on target surface

      图4中的${F_1}$${F_2}$分别表示目标的起始探测点和终止探测点,对应在大地坐标系中的位矢分别为${{{r}}_{t1}}$${{{r}}_{t{{2}}}}$,两点间的距离为$\left| {{{{r}}_{t{{2}}}}{{ - }}{{{r}}_{t{{1}}}}} \right|$,探测两点的时间间隔为t,则有:

      $$\frac{{\left| {{{{{r}}}_{t{{2}}}}{{ - }}{{{{r}}}_{t{{1}}}}} \right| + L}}{{\left| {{}^{{0}}{{{{v}}}_{{t}}}} \right|}} = t$$ (18)

      公式(18)中L表示目标体的特征长度,由于$\left| {{}^0{{{v}}_t}} \right|$t$\left| {{{{r}}_{t{\rm{2}}}}{\rm{ - }}{{{r}}_{t{\rm{1}}}}} \right|$均已知,则根据公式(18)即可求解出L,从而计算探测结束时刻目标几何中心点${F_c}$的位矢${{{r}}_c}$

      $${{{r}}_c} = {{{r}}_{{{t}}2}} + \frac{L}{2}\frac{{{}^0{{{v}}_t}}}{{\left| {{}^0{{{v}}_t}} \right|}}$$ (19)

      假设弹丸静止时,爆炸后的破片以速度v沿垂直弹轴方向飞散,作用半径为R,且在飞散过程中速度不衰减。探测结束ΔT时间后起爆弹丸,破片与目标交会模型如图5所示。

      图  5  破片与目标交会模型。(a)大地坐标系中破片与目标交会过程;(b)弹体坐标系中起爆方位角

      Figure 5.  Fragment and target intersection model. (a) Intersection process of fragments and targets in the geodetic coordinate system; (b) Initiation azimuth in the projectile coordinate system

      图5(a)中点${M_2}$${M_3}$分别为弹丸探测结束和最终起爆位置,${{{r}}_{m2}}$${{{r}}_{m3}}$分别表示点${M_2}$${M_3}$的位矢,${{{v}}}$${{{v}}_{{m}}}$分别表示破片静态飞散速度矢量和弹丸运动速度矢量。考虑弹丸速度,战斗部爆轰后破片作用范围如图5(a)中阴影部分所示,则破片的绝对速度${{{{v}}}_{{h}}}$可表示为:

      $${{{{v}}}_h} = {{{v}}} + {{{{v}}}_m}$$ (20)

      图5(b)所示,起爆方位角${\alpha _{{m}}}$为弹体坐标系中破片静态飞散速度矢量v${y_m}$轴之间的夹角,则有:

      $${{{{v}}}_{{h}}} = \left[ {\begin{array}{*{20}{c}} {\left| {{{{{v}}}_m}} \right|} \\ {\left| {{{v}}} \right|\cos {\alpha _m}} \\ {\left| {{{v}}} \right|\sin {\alpha _m}} \end{array}} \right]$$ (21)

      最佳起爆条件应使弹丸起爆$\Delta {T_h}$时间后破片恰好击中目标的易损部位即其特征几何中心点,此时距探测结束已历经了ΔtTTh的时间,目标的特征几何中心点的位置为:

      $${{{{r}}}_{t3}} = {{{{r}}}_{{c}}} + (\Delta T + \Delta {T_h}){{{{v}}}_t}$$ (22)

      最佳起爆条件满足:

      $${{{{r}}}_{m3}} + \Delta {T_{{h}}}{}^0{{{{v}}}_h} = {{{{r}}}_{{c}}} + (\Delta T + \Delta {T_h}){}^0{{{{v}}}_t}$$ (23)

      其中,

      $${{{{r}}}_{m3}} = {{{{r}}}_{m2}} + \Delta T{{{{v}}}_m}$$ (24)

      联立公式(19)、(23)、(24),得:

      $$({{{{v}}}_{{m}}} - {{{{v}}}_t})\Delta T + (C_0^m\left[ {\begin{aligned} {\left| {{{{{v}}}_{{m}}}} \right|}\;\;\;\\ {\left| {{{v}}} \right|{{\cos}}{\alpha _{{m}}}} \\ {\left| {{{v}}} \right|\sin {\alpha _{{m}}}} \end{aligned}} \right]{{ - }}{{{{v}}}_t})\Delta {T_{{h}}} = \frac{L}{2}\frac{{{{{{v}}}_t}}}{{\left| {{{{{v}}}_t}} \right|}} + {{{r}}}_{mt2}^{}$$ (25)

      式中:${{{{r}}}_{m{{t}}2}} = {{{{r}}}_{t2}} - {{{{r}}}_{m2}}$,表示探测结束时刻的弹目距离,其值可以通过探测点回波脉冲信号直接得到。在大地坐标系下,弹丸速度${{{{v}}}_m}$为:

      $${{{{v}}}_{{m}}} = C_0^m{\left[ {\begin{array}{*{20}{c}} {\left| {{{{{v}}}_m}} \right|}&0&0 \end{array}} \right]^{\bf{T}}}$$ (26)

      公式(25)中,${{{{v}}}_m}$${{{{v}}}_t}$$\left| {{{v}}} \right|$${{{r}}}_{{{mt2}}}^{}$L可以通过多探测点数据,基于目标速度模型求解,代入即可求解同步扫描周视激光脉冲探测系统的最佳起爆时间$\Delta T$、最佳起爆方位角${\alpha _{{m}}}$和起爆后破片击中目标时间$\Delta {T_h}$

    • 根据上述建立的最佳起爆模型可知,最佳起爆时间$\Delta T$和最佳起爆方位角${\alpha _{{m}}}$与弹丸自身俯仰角${\beta _{{m}}}$、偏航角${\psi _{{m}}}$、速度大小$\left| {{{{{v}}}_{{m}}}} \right|$、破片静态飞散速度大小$\left| {{{v}}} \right|$、目标速度${{{{v}}}_{{t}}}$、弹目距离${{{r}}}_{mt2}^{}$以及目标特征长度L有关。

    • 假设已知${\psi _{{m}}} = - \dfrac{\pi }{4}$$\left| {{{{{v}}}_{{m}}}} \right|$=500 m/s,${{{r}}}_{mt{\rm{2}}}^{}$=$[{\rm{3}}\quad{\rm{5}}\quad{\rm{6}}]^{\rm{T}}$ m,R=10 m,L=2.5 m,${{{{v}}}_{{t}}}$=700$\left[{ - \dfrac{{\sqrt 3 }}{6}}\quad { - \dfrac{{\sqrt 3 }}{2}}\quad { - \dfrac{1}{3}}\right]^{\bf{T}}$ m/s,代入公式(25),可得最佳起爆时间和方位角与破片飞散速度$\left| {{{v}}} \right|$的变化关系如图6所示。

      图  6  最佳起爆时间和方位角与破片飞散速度变化关系。(a)最佳起爆时间与破片飞散速度关系图;(b)最佳起爆方位角与破片飞散速度关系图

      Figure 6.  Relationship between the optimal initiation time, initiation azimuth and the fragment speed. (a)Optimal initiation time vs fragment speed; (b)Optimal initiation azimuth vs fragment speed

      图6可知在既定破片作用范围内,弹丸俯仰角${\beta _m}$$\pi /6$$\pi /4$时,破片飞散速度越大,最佳起爆时间越长,${\beta _m}$$\pi /3$时,破片飞散速度越大,最佳起爆时间越短;弹丸俯仰角确定时,最佳起爆方位角为定值,不随破片飞散速度发生变化;同一破片飞散速度情况下,${\beta _m}$$\pi /3$时,最佳起爆方位角最小,${\beta _m}$$\pi /4$时,最佳起爆方位角最大。

    • 假设已知${\psi _{{m}}} = - \dfrac{\pi }{4}$$\left| {{{v}}} \right|$=2 000 m/s,${{{r}}}_{mt{\rm{2}}}^{}$=$[{\rm{3}} \quad{\rm{5}}\quad{\rm{6}}]^{\rm{T}}$ m,R=10 m,L=2.5 m,${{{{v}}}_{{t}}}$=700$\left[ { - \frac{{\sqrt 3 }}{6}}\quad{ - \frac{{\sqrt 3 }}{2}}\quad{ - \frac{1}{3}}\right]^{\bf{T}}$ m/s,代入公式(25),可得最佳起爆时间和方位角与弹丸速度$\left| {{{{{v}}}_{\rm{m}}}} \right|$的变化关系如图7所示。

      图  7  最佳起爆时间和方位角与弹丸速度变化关系。(a)最佳起爆时间与弹丸速度关系图;(b)最佳起爆方位角与弹丸速度关系图

      Figure 7.  Relationship between the optimal initiation time, initiation azimuth and the projectile speed. (a)Optimal initiation time vs projectile speed; (b)Optimal initiation azimuth vs projectile speed

      由图7可知在既定破片作用范围内,$ {\beta _m}$$ \pi /6$$ \pi /4$$ \pi /3$时,弹丸速度越大,最佳起爆时间越短,最佳起爆方位角越小;弹丸自身俯仰角越大,最佳起爆方位角变化越缓慢。相同目标特征长度情况下,$ {\beta _m}$分别为$ \pi /6$$ \pi /4$$ \pi /3$时,弹丸俯仰角越大,最佳起爆时间越短,最佳起爆方位角越大。

    • 假设已知${\psi _{{m}}} = - \dfrac{\pi }{4}$$\left| {{{v}}} \right|$=2000 m/s,$\left| {{{{{v}}}_{{m}}}} \right|$=500 m/s,${{{r}}}_{mt{{2}}}^{}$=$\left| {{{{{r}}}_{{{mt2}}}}} \right|$$\left[ {\dfrac{3}{{\sqrt {70} }}}\quad{\dfrac{5}{{\sqrt {70} }}}\quad{\dfrac{6}{{\sqrt {70} }}} \right]^{\bf{T}}$ m,R=10 m,${{{{v}}}_{{t}}}$=700 $\left[ { - \dfrac{{\sqrt 3 }}{6}}\quad{ - \dfrac{{\sqrt 3 }}{2}}\quad{ - \dfrac{1}{3}}\right]^{\bf{T}}$ m/s,代入公式(25),可得最佳起爆时间和方位角与弹目距离$\left| {{{{{r}}}_{{{mt2}}}}} \right|$的变化关系如图8所示。

      图  8  最佳起爆时间和方位角与弹目距离变化关系。(a)最佳起爆时间与弹目距离关系图;(b)最佳起爆方位角与弹目距离关系图

      Figure 8.  Relationship between the optimal initiation time, initiation azimuth and the distance between the projectile with the target.(a)Optimal initiation time vs the distance between the projectile and target;(b)Optimal initiation azimuth vs the distance between projectile and target

      图8可知在既定破片作用范围内,弹丸俯仰角${\beta _m}$$\pi /6$$\pi /4$$\pi /3$时,令最佳起爆时间为0时弹目距离为临界弹目距离,弹目距离小于临界弹目距离时,最佳起爆时间为负值,表示探测结束时已错过最佳起爆时机,此时应立即起爆战斗部;弹目距离大于临界弹目距离时,最佳起爆时间越短,最佳起爆方位角越大,最佳起爆时间随弹目距离变化呈线性关系;弹丸自身俯仰角越大,最佳起爆方位角变化越缓慢。相同目标特征长度情况下,${\beta _m}$分别为$\pi /6$$\pi /4$$\pi /3$时,弹丸俯仰角越大,最佳起爆时间越短,最佳起爆方位角越大。

    • 假设已知${\psi _{{m}}} = - \dfrac{\pi }{4}$$\left| {{{v}}} \right|$=2 000 m/s,$\left| {{{{{v}}}_{{m}}}} \right|$=500 m/s,${{{r}}}_{mt{\rm{2}}}^{}$=$[3\;\; 5\;\; 6 ]^{\rm{T}}$ m,R=10 m,${{{{v}}}_{{t}}}$=700$\left[{ - \dfrac{{\sqrt 3 }}{6}}\;\;{ - \dfrac{{\sqrt 3 }}{2}}\;\;{ - \dfrac{1}{3}} \right]^{\bf{T}}$ m/s,代入公式(25),可得最佳起爆时间和方位角与目标特征长度L的变化关系如图9所示。

      图  9  最佳起爆时间和方位角与目标特征长度变化关系。(a)最佳起爆时间与目标特征长度关系图;(b)最佳起爆方位角与目标特征长度关系图

      Figure 9.  Relationship between the optimal initiation time, initiation azimuth and the target characteristic length.(a)Optimal initiation time vs the target characteristic length;(b)Optimal initiation azimuth vs the target characteristic length

      图9可知在既定破片作用范围内,弹丸俯仰角${\beta _m}$$\pi /6$$\pi /4$$\pi /3$时,目标特征长度越大,最佳起爆时间越短,最佳起爆方位角越大,最佳起爆时间随目标特征长度呈线性关系;弹丸自身俯仰角越大,最佳起爆方位角变化越缓慢。相同目标特征长度情况下,${\beta _m}$分别为$\pi /6$$\pi /4$$\pi /3$时,弹丸俯仰角越大,最佳起爆时间越短,最佳起爆方位角越大。

    • 若已知${\psi _{{m}}} = - \dfrac{\pi }{4}$$\left| {{{v}}} \right|$=2 000 m/s,$\left| {{{{{v}}}_{{m}}}} \right|$=500 m/s,${{{r}}}_{mt{\rm{2}}}^{}$=$\left[3\;5\;6\right]^{\rm{T}}$ m,R=10 m,L=2.5 m,${{{{v}}}_{{t}}}$=$\left| {{{{{v}}}_{{t}}}} \right|$$\left[{ - \dfrac{{\sqrt 3 }}{6}}\;{ - \dfrac{{\sqrt 3 }}{2}}\;{ - \dfrac{1}{3}}\right]^{\bf{T}}$ m/s。代入公式(25)可得最佳起爆时间和方位角与目标速度的变化关系如图10所示。

      图  10  最佳起爆时间和方位角与目标速度变化关系。(a)最佳起爆时间与目标速度关系图;(b)最佳起爆方位角与目标速度关系图

      Figure 10.  Relationship between the optimal initiation time, initiation azimuth and the target speed. (a)Optimal initiation time vs target speed; (b)Optimal initiation azimuth vs target speed

      图10可知在既定破片作用范围内,弹丸俯仰角${\beta _m}$$\pi /6$$\pi /4$$\pi /3$时,目标速度越大,最佳起爆时间越短,最佳起爆方位角越大;弹丸自身俯仰角越大,最佳起爆方位角变化越缓慢。目标速度相同情况下,${\beta _{\rm{m}}}$$\pi /6$$\pi /4$$\pi /3$时,弹丸俯仰角越大,最佳起爆时间越短,最佳起爆方位角越大。

    • 文中阐述了同步扫描周视脉冲激光引信的工作原理,分析了目标探测过程中可探测的多个目标位置点,计算了各探测点在弹体坐标系中的位置坐标。根据各探测点的距离和方位信息,建立了基于多探测点的目标速度模型,通过目标速度模型对目标速度进行解算,与攻击目标正常速度范围进行比对,来辨别目标真伪。在目标速度模型的基础上,以破片攻击目标特征几何中心为目的,求解出引信最佳起爆时间和方位角,并分析了最佳起爆模型的影响因素及其对最佳起爆时间和方位角的影响。

参考文献 (10)

目录

    /

    返回文章
    返回