• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

基于机械振动理论的垂直侵彻弹靶作用模型

程祥利 赵慧 李林川 叶海福

引用本文:
Citation:

基于机械振动理论的垂直侵彻弹靶作用模型

    作者简介: 程祥利(1984- ),男,博士,助理研究员,chengxiangli126@126.com;
    通讯作者: 叶海福, yehaifu@126.com
  • 中图分类号: O385

Projectile target response model for normal penetration process based on mechanical vibration theory

    Corresponding author: YE Haifu, yehaifu@126.com
  • CLC number: O385

  • 摘要: 为了给侵彻引信抗高过载优化设计提供准确的力学输入,将机械振动理论引入侵彻过程建模领域,提出了一种侵彻战斗部刚体运动与一阶轴向振动相结合的垂直侵彻弹靶作用模型。在垂直侵彻过程受力分析的基础上,基于牛顿第二定律建立了战斗部刚体运动模型,基于单自由度弹簧-质量-阻尼系统建立了战斗部一阶轴向振动模型,并采用数值积分的方法获得了垂直侵彻过程中各物理量的变化规律。和火炮试验实测加速度信号的对比分析结果表明:考虑战斗部一阶轴向振动后的垂直侵彻弹靶作用模型能更准确地描述侵彻过程,能更有效地指导侵彻引信的抗高过载优化设计。
  • 图 1  垂直侵彻过程受力分析图

    Figure 1.  Force diagram in normal penetration process

    图 2  单自由度弹簧-质量-阻尼系统

    Figure 2.  Single DOF spring-mass-damper system

    图 3  程序流程图

    Figure 3.  Flow diagram for programming

    图 4  垂直侵彻过程中的刚体过载变化曲线

    Figure 4.  Rigid acceleration curve in normal penetration process

    图 5  垂直侵彻过程中的弹性过载变化曲线

    Figure 5.  Elastic acceleration curve in normal penetration process

    图 6  弹性过载的频谱分析结果

    Figure 6.  Frequency spectrum of elastic acceleration signal

    图 7  战斗部一阶轴向振动的频率响应特性

    Figure 7.  Frequency response characteristics of the first order axial vibration for projectiles

    图 8  垂直侵彻过程中的总过载变化曲线

    Figure 8.  Total acceleration curve in normal penetration process

    图 9  战斗部总过载的频谱分析结果

    Figure 9.  Frequency spectrum of total acceleration signal

    图 10  靶标布置

    Figure 10.  Schematic diagram of the target

    图 11  试验弹

    Figure 11.  Schematic diagram of the projectile

    图 12  实测加速度信号

    Figure 12.  Acceleration data collected

    图 13  实测数据的频谱

    Figure 13.  Frequency spectrum of data collected

    表 1  每层靶的碰靶时刻

    Table 1.  Penetration moment of every layer target

    层数计算结果/ms高速摄像判读结果/ms绝对误差/ms相对误差/%
    第 1 层 0.00 0.00.00
    第 2 层 4.20 4.00.205.00
    第 3 层 8.30 8.00.303.75
    第 4 层12.5012.00.504.17
    第 5 层16.8515.90.955.97
    第 6 层21.3020.31.004.92
    第 7 层25.9025.00.903.60
    第 8 层30.6529.80.852.85
    下载: 导出CSV
  • [1] 李晓峰. 侵彻弹药引信技术 [M]. 北京: 国防工业出版社, 2016: 1−3.
    [2] 王伟力, 黄雪峰, 杨雨潼. 半穿甲战斗部侵彻过程中装药安定性研究 [J]. 海军航空工程学院学报, 2010, 25(1): 79–82. DOI: 10.3969/j.issn.1673-1522.2010.01.019.
    WANG Weili, HUANG Xuefeng, YANG Yutong. Research on the grain safety during the penetration process of semi-armor-piercing warhead [J]. Journal of Naval Aeronautical and Astronautical University, 2010, 25(1): 79–82. DOI: 10.3969/j.issn.1673-1522.2010.01.019.
    [3] 张建新. 侵彻引信炸点控制理论及试验研究 [D]. 南京: 南京理工大学, 2012: 1−2.
    [4] 赵生伟, 初哲, 李明. 抗侵彻过载战斗部装药安定性实验研究 [J]. 兵工学报, 2010, 31(S1): 284–287.
    ZHAO Shengwei, CHU Zhe, LI Ming. Experiment investigation on stability of explosive in anti-overload warhead [J]. Acta Armamentarii, 2010, 31(S1): 284–287.
    [5] 赵南, 王可慧, 李明, 等. 薄壁弹体高速侵彻钢筋混凝土实验研究 [J]. 实验力学, 2017, 32(4): 573–579. DOI: 10.7520/1001-4888-16-179.
    ZHAO Nan, WANG Kehui, LI Ming, et al. Experimental study of high speed penetration of thin-wall projectile in steel reinforced concreete [J]. Journal of Experimental Mechanics, 2017, 32(4): 573–579. DOI: 10.7520/1001-4888-16-179.
    [6] FORRESTAL M J, TZOU D Y. A spherical cavity-expansion penetration model for concrete targets [J]. International Journal of Solids and Structures, 1997, 34(31-32): 4127–4146. DOI: 10.1016/S0020-7683(97)00017-6.
    [7] BENDOR G, DUBINSKY A, ELPERIN T. Analytical solution for penetration by rigid conical impactors using cavity expansion models [J]. Mechanics Research Communications, 2000, 27(2): 185–189. DOI: 10.1016/S0093-6413(00)00080-X.
    [8] FORRESTAL M J, FREW D J, HICKERSON J P, et al. Penetration of concrete targets with deceleration-time measurements [J]. International Journal of Impact Engineering, 2003, 28(5): 479–497. DOI: 10.1016/S0734-743X(02)00108-2.
    [9] 虞青俊, 李玉龙, 金连宝, 等. 侵彻多层混凝土目标弹丸过载特性研究 [J]. 探测与控制学报, 2007, 29(1): 13–17. DOI: 10.3969/j.issn.1008-1194.2007.01.004.
    YU Qingjun, LI Yulong, JIN Lianbao, et al. Research of deceleration-time curves during penetration of multi-plate concrete targets [J]. Journal of Detection & Control, 2007, 29(1): 13–17. DOI: 10.3969/j.issn.1008-1194.2007.01.004.
    [10] 周栋, 吴俊斌. 动能战斗部侵彻混凝土力学响应研究 [J]. 战术导弹技术, 2012(4): 16–19.
    ZHOU Dong, WU Junbin. Research on penetrating concrete effective of kinetic warhead [J]. Tactical Missile Technology, 2012(4): 16–19.
    [11] 徐文亮, 何春, 李朝君. 侵彻爆破型战斗部侵彻性能总体评估系统研究 [J]. 战术导弹技术, 2013(1): 93–100.
    XU Wenliang, HE Chun, LI Chaojun. General evaluation study of penetration warhead’s penetrate capability [J]. Tactical Missile Technology, 2013(1): 93–100.
    [12] 皮爱国, 黄风雷. 大长细比结构弹体侵彻2024-O铝靶的弹塑性动力响应 [J]. 爆炸与冲击, 2008, 28(3): 252–260. DOI: 10.11883/1001-1455(2008)03-0252-09.
    PI Aiguo, HUANG Fenglei. Elastic-plastic dynamic response of slender projectiles penetrating into 2024-O aluminum targets [J]. Explosion and Shock Waves, 2008, 28(3): 252–260. DOI: 10.11883/1001-1455(2008)03-0252-09.
    [13] 王琳, 王富耻, 王鲁, 等. 空心弹体垂直侵彻混凝土靶板的应变测试研究 [J]. 北京理工大学学报, 2002, 22(4): 453–456. DOI: 10.3969/j.issn.1001-0645.2002.04.014.
    WANG Lin, WANG Fuchi, WANG Lu, et al. Strain measurement in hollow projectiles impacting concrete targets [J]. Journal of Beijing Institute of Technology, 2002, 22(4): 453–456. DOI: 10.3969/j.issn.1001-0645.2002.04.014.
    [14] 程兴旺, 王富耻, 王鲁, 等. 钨合金壳体侵彻混凝土靶板过程壳体应变的实验测试 [J]. 兵工学报, 2004, 25(1): 102–105. DOI: 10.3321/j.issn:1000-1093.2004.01.026.
    CHENG Xingwang, WANG Fuchi, WANG Lu, et al. Experimental study on the strain history of critical section during a normal penetration of tungsten alloy shell into a concrete target [J]. Acta Armamentarii, 2004, 25(1): 102–105. DOI: 10.3321/j.issn:1000-1093.2004.01.026.
    [15] 韩学平, 芮筱亭, 王国平, 等. 基于小波的弹性弹丸膛内引信过载研究 [J]. 系统仿真学报, 2008, 20(13): 3496–3499.
    HAN Xueping, RUI Xiaoting, WANG Guoping, et al. Research of fuze overload in bore of flexibility pills based on wavelet method [J]. Journal of System Simulation, 2008, 20(13): 3496–3499.
    [16] 陈学强, 闫明明, 徐晓辉, 等. 微加速度计在高冲击下的断裂失效分析 [J]. 仪表技术与传感器, 2014(2): 16–19. DOI: 10.3969/j.issn.1002-1841.2014.02.006.
    CHEN Xueqiang, YAN Mingming, XU Xiaohui, et al. Fracture failure analysis of micro-accelerometer under high impact [J]. Instrument Technique and Sensor, 2014(2): 16–19. DOI: 10.3969/j.issn.1002-1841.2014.02.006.
    [17] 刘燕芳, 郭海波, 潘启智, 等. 多层陶瓷电容器的失效分析 [J]. 电子元件与材料, 2010, 29(11): 72–74. DOI: 10.3969/j.issn.1001-2028.2010.11.021.
    LIU Yanfang, GUO Haibo, PAN Qizhi, et al. Failure analysis of multi-layer ceramic capacitor [J]. Electronic Components and Materials, 2010, 29(11): 72–74. DOI: 10.3969/j.issn.1001-2028.2010.11.021.
    [18] 何涛, 文鹤鸣. 靶体响应力函数的确定方法及其在侵彻力学中的应用 [J]. 中国科学技术大学学报, 2007, 37(10): 1249–1261. DOI: 10.3969/j.issn.0253-2778.2007.10.017.
    HE Tao, WEN Heming. Determination of the analytical forcing function of target response and its applications in penetration mechanics [J]. Journal of University of Science and Technology of China, 2007, 37(10): 1249–1261. DOI: 10.3969/j.issn.0253-2778.2007.10.017.
    [19] 刘波, 杨黎明, 李东杰, 等. 侵彻弹体结构纵向振动频率特性分析 [J]. 爆炸与冲击, 2018, 38(3): 677–682. DOI: 10.11883/bzycj-2016-0282.
    LIU Bo, YANG Liming, LI Dongjie, et al. Analysis of axial vibration frequency for projectile structure in penetration [J]. Explosion and Shock Waves, 2018, 38(3): 677–682. DOI: 10.11883/bzycj-2016-0282.
    [20] 罗梦翔, 刘涛, 蔡国平. 导弹振动的动力学建模和频率分析 [J]. 中国科技论文, 2015, 10(16): 1924–1927. DOI: 10.3969/j.issn.2095-2783.2015.16.012.
    LUO Mengxiang, LIU Tao, CAI Guoping. Dynamics modeling and frequency analysis of missile vibration [J]. China Science Paper, 2015, 10(16): 1924–1927. DOI: 10.3969/j.issn.2095-2783.2015.16.012.
    [21] 季文美. 机械振动 [M].北京: 科学出版社, 2016: 68−69.
  • [1] 李睿黄正祥祖旭东肖强强贾鑫 . 靶板在爆炸成型弹丸垂直侵彻下的层裂. 爆炸与冲击, 2018, 38(5): 1039-1044. doi: 10.11883/bzycj-2017-0055
    [2] 王文杰张先锋邓佳杰郑应民刘闯 . 椭圆截面弹体侵彻砂浆靶规律分析. 爆炸与冲击, 2018, 38(1): 164-173. doi: 10.11883/bzycj-2017-0020
    [3] 邢柏阳刘荣忠张东江陈亮侯云辉郭锐 . 变截面爆炸成型弹丸垂直侵彻装甲钢板靶后破片质量模型. 爆炸与冲击, 2019, 39(7): 074202-1-074202-11. doi: 10.11883/bzycj-2018-0187
    [4] 朱建生赵国志杜忠华王贤治 . PELE垂直侵彻薄靶的机理分析. 爆炸与冲击, 2009, 29(3): 281-288. doi: 10.11883/1001-1455(2009)03-0281-08
    [5] 尹放林王明洋钱七虎 . 弹体垂直侵彻深度工程计算模型. 爆炸与冲击, 1997, 17(4): 333-339.
    [6] 吴群彪沈培辉方海峰范纪华蔡李花 . 前置组合杆体垂直侵彻钢靶简化模型. 爆炸与冲击, 2019, 39(1): 013302-1-013302-12. doi: 10.11883/bzycj-2017-0287
    [7] 欧阳春赵国志杜中华李文彬 . 弹丸垂直侵彻钢筋混凝土介质的工程解析模型. 爆炸与冲击, 2004, 24(3): 273-277.
    [8] 王政倪玉山曹菊珍王元书张文 . 卵形头部刚性弹侵彻厚靶的半解析模型. 爆炸与冲击, 2004, 24(3): 212-218.
    [9] 邓勇军宋文杰陈小伟姚勇 . 钢筋混凝土靶侵彻的可压缩弹-塑性动态空腔膨胀阻力模型. 爆炸与冲击, 2018, 38(5): 1023-1030. doi: 10.11883/bzycj-2017-0043
    [10] 汪斌曹仁义谭多望 . 大质量高速动能弹侵彻钢筋混凝土的实验研究. 爆炸与冲击, 2013, 33(1): 98-102. doi: 10.11883/1001-1455(2013)01-0098-05
    [11] 何丽灵陈小伟范瑛 . 先进钻地弹高速侵彻实验中质量磨蚀金相分析. 爆炸与冲击, 2012, 32(5): 515-522. doi: 10.11883/1001-1455(2012)05-0515-08
    [12] 王峰王肖钧胡秀章刘文韬 . 卵形杆弹对铝靶的斜侵彻. 爆炸与冲击, 2005, 25(3): 265-270. doi: 10.11883/1001-1455(2005)03-0265-06
    [13] 陈克黄德武S.G.Psakhie . 模拟长杆弹侵彻混凝土靶的MCA方法. 爆炸与冲击, 2004, 24(2): 127-132.
    [14] 李继承陈小伟 . 刚性弹侵彻不同靶材的侵彻深度比较. 爆炸与冲击, 2009, 29(3): 225-230. doi: 10.11883/1001-1455(2009)03-0225-06
    [15] 张欣欣武海军黄风雷段卓平皮爱国 . 刻槽弹侵彻混凝土受力模型研究. 爆炸与冲击, 2016, 36(1): 75-80. doi: 10.11883/1001-1455(2016)01-0075-06
    [16] 吴成沈晓军王晓鸣姚文进 . 细观混凝土靶抗侵彻数值模拟及侵彻深度模型. 爆炸与冲击, 2018, 38(6): 1364-1371. doi: 10.11883/bzycj-2017-0123
    [17] 陈晓周宏 . 叠层靶板弹击实验及弹道侵彻机理的数值模拟研究. 爆炸与冲击, 2003, 23(6): 509-515.
    [18] 许沭华王肖钧张刚明刘文韬 . 预扭转钨合金杆弹侵彻钢靶的数值模拟. 爆炸与冲击, 2002, 22(3): 210-215.
    [19] 王金涛余文力王涛罗永锋王少龙 . SPH 算法在长杆弹侵彻多层间隔靶中的应用. 爆炸与冲击, 2011, 31(5): 533-569. doi: 10.11883/1001-1455(2011)05-0533-07
    [20] 蒋东李永池于少娟邓世春 . 钨合金长杆弹侵彻约束AD95陶瓷复合靶. 爆炸与冲击, 2010, 30(1): 91-95. doi: 10.11883/1001-1455(2010)01-0091-05
  • 加载中
图(13)表(1)
计量
  • 文章访问数:  69
  • HTML全文浏览量:  222
  • PDF下载量:  18
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-07-03
  • 录用日期:  2018-09-08
  • 网络出版日期:  2019-08-25
  • 刊出日期:  2019-09-01

基于机械振动理论的垂直侵彻弹靶作用模型

    作者简介:程祥利(1984- ),男,博士,助理研究员,chengxiangli126@126.com
    通讯作者: 叶海福, yehaifu@126.com
  • 中国工程物理研究院电子工程研究所,四川 绵阳 621999

摘要: 为了给侵彻引信抗高过载优化设计提供准确的力学输入,将机械振动理论引入侵彻过程建模领域,提出了一种侵彻战斗部刚体运动与一阶轴向振动相结合的垂直侵彻弹靶作用模型。在垂直侵彻过程受力分析的基础上,基于牛顿第二定律建立了战斗部刚体运动模型,基于单自由度弹簧-质量-阻尼系统建立了战斗部一阶轴向振动模型,并采用数值积分的方法获得了垂直侵彻过程中各物理量的变化规律。和火炮试验实测加速度信号的对比分析结果表明:考虑战斗部一阶轴向振动后的垂直侵彻弹靶作用模型能更准确地描述侵彻过程,能更有效地指导侵彻引信的抗高过载优化设计。

English Abstract

  • 现代战争中,越来越依靠超高速侵彻战斗部(弹)对防御工事、水面舰艇、机场跑道、地面建筑、飞机掩蔽库等高强度硬目标(靶)进行毁伤[1]。侵彻过程中,弹靶作用机理异常复杂,给战斗部、装药、引信的抗高过载优化设计带来很大困难[2-3],已成为具有重大军事需求的力学难题。

    然而,现有研究主要关注战斗部的结构强度、侵彻能力以及装药安定性[4-5],较少涉及到引信,导致引信抗高过载优化设计的力学输入不明确。针对上述难点,一种行之有效的解决方法是通过建立合理、有效、简化的弹靶作用模型获得侵彻过程中的加速度信息并等效为引信设计输入。

    传统的弹靶作用模型建模方法主要考虑战斗部刚体运动[6-11],即侵彻过程中战斗部不发生弹性变形。首先,在侵彻过程受力分析的基础上采用空腔膨胀理论和微分面元法获得侵彻过程中的弹靶作用力,随后,基于牛顿第二定律建立刚体运动模型,即侵彻过程中的运动微分方程,最后,采用数值积分的方法求解刚体运动模型,获得侵彻过程中各物理量(含位移、速度、加速度)的变化规律,并分析各种因素的影响。

    但是,打靶试验结果显示战斗部在侵彻过程中都会发生不同程度的变形[12-14],在理论上都应视为弹性体[15]。另外,部分电子元器件在正弦或交流载荷作用下更易损坏[16-17],说明侵彻过程中的战斗部弹性变形会影响内部组件的生存性及工作可靠性,特别是含多种电子元器件的引信。

    为了揭示战斗部弹性变形对侵彻引信受力环境的影响规律,本文将机械振动理论引入侵彻过程建模领域,提出一种战斗部刚体运动与一阶轴向振动相结合的垂直侵彻弹靶作用模型。首先,将战斗部等效为具有阻尼特性的线性弹簧,基于单自由度弹簧-质量-阻尼系统建立一阶轴向振动模型,并与传统的刚体运动模型相叠加作为修正后的弹靶作用模型;其次,在垂直侵彻过程受力分析的基础上,采用空腔膨胀理论和微分面元法获得弹靶作用力,并作为弹靶作用模型的输入;最后,采用数值积分的方法获得垂直侵彻过程中各物理量的变化规律,并和实测数据进行对比以验证模型的准确性。

    • 作为求解侵彻问题的一种近似方法,空腔膨胀理论根据质量守恒定律和动量守恒定律建立法向应力与空腔膨胀速度的关系,并通过对弹头表面积分得到弹靶作用力[18]

      战斗部表面微元处的法向应力σn与空腔膨胀速度Vn、靶标材料特征参数AB密切相关,如式(1)所示。参数AB的表达式分别如式(2)和式(3)所示[10]

      ${\sigma _n} = A + BV_n^2$

      $A = \frac{{2{\tau _0}}}{3}\left( {1 - \ln {\eta ^*}} \right) = \frac{2}{3}\left( {1 - \ln {\eta ^*}} \right)\left( {1 - \frac{\lambda }{3}} \right){\sigma _0}$

      $B = {\rho _0}\left[ {\frac{{3{\tau _0}}}{E} + {\eta ^*}{{\left( {1 - \frac{{3{\tau _0}}}{{2E}}} \right)}^2} + \frac{{3{{ {{\eta ^*}} }^{\frac{2}{3}}} - {\eta ^*}\left( {4 - {\eta ^*}} \right)}}{{2\left( {1 - {\eta ^*}} \right)}}} \right]{\left[ {{{\left( {1 + \frac{{{\tau _0}}}{{2E}}} \right)}^2} - \left( {1 - {\eta ^*}} \right)} \right]^{{\rm{ }} - \frac{2}{3}}}$

      式中:η*为锁变体积常数,ρ0为初始密度,λ为强化模量参数,E为弹性模量,τ0为侧向抗剪强度,σ0为抗压强度。

      战斗部表面微元处的切向应力στ如下式所示。

      ${\sigma _\tau } = \mu {\sigma _n}$

      式中:μ为动摩擦因数。

    • 假设战斗部垂直侵彻过程中攻角始终为零,即战斗部不发生姿态偏转。此时,战斗部主要承受靶标施加的侵彻阻力作用,即弹靶作用力Fx,如图1所示。图中,H为靶标厚度,x为侵彻深度,V为侵彻速度,L为卵形弹头的长度。

      图  1  垂直侵彻过程受力分析图

      Figure 1.  Force diagram in normal penetration process

      弹靶作用力可根据微分面元法求解。对于卵形弹垂直侵彻,在垂直于弹轴方向的圆截面上取微面元,记为dS,如式(5)所示。

      ${\rm d}S = 2{\text{π}} {\chi^2}{D^2}\left( {\sin \theta - \frac{{\chi - 0.5}}{{\chi}}} \right){\rm d}\theta $

      式中:$\chi $为弹头系数,D为弹径,θ为卵形母线的圆心角。

      微面元dS所受法向力dFn和切向力dFτ如式(6)所示。

      $\left\{ \begin{array}{l} {\rm d}{F_n} = {\sigma _n}{\rm d}S \\ {\rm d}{F_\tau } = {\sigma _\tau }{\rm d}S \end{array} \right.$

      将微面元dS的法向力和切向力向弹轴方向投影,可得微面元dS的轴向力dFx,如式(7)所示。

      ${\rm d}{F_x} = {\rm d}{F_n}\cos \theta + {\rm d}{F_\tau }\sin \theta = 2{\text{π}} {\chi^2}{D^2}{\sigma _n}\left( {\sin \theta - \frac{{\chi - 0.5}}{{\chi}}} \right)\left( {\cos \theta + \mu \sin \theta } \right)d\theta $

      对微面元dS的轴向力dFx沿战斗部表面进行积分,可得整个战斗部的轴向力,即弹靶作用力,如式(8)所示。

      ${F_x} = \int_{{\theta _1}}^{{\theta _2}} {{\rm d}{F_x}} = 2{\text{π}} {\chi^2}{D^2}\int_{{\theta _1}}^{{\theta _2}} {{\sigma _n}\left( {\sin \theta - \frac{{\chi - 0.5}}{{\chi}}} \right)\left( {\cos \theta + \mu \sin \theta } \right){\rm d}\theta } $

    • 侵彻过程中,战斗部既有刚体运动,又有弹性振动,可分别建立刚体运动模型和弹性振动模型。

    • 垂直侵彻时,弹靶作用力使战斗部做减速运动,可简化为质心的刚体运动。假设战斗部头部的位移为x1,则战斗部各个位置的位移均为x1。根据牛顿第二定律可得运动微分方程,即刚体运动方程为:

      $\frac{{{{\rm d}^2}{x_1}}}{{{\rm d}{t^2}}} = \frac{{{\rm d}V}}{{{\rm d}t}} = - \frac{{{F_x}}}{M}$

      式中:t为侵彻时间,M为战斗部质量。

    • 垂直侵彻时,弹靶作用力会激发起战斗部的多种振动模态,使得战斗部剧烈振动。由于轴向载荷起主要作用,可近似用一阶轴向振动模型描述弹性振动。

      假设战斗部为两端自由的一维圆杆,可根据应力波传播理论计算出战斗部一阶轴向振动的固有频率f1,理论值如下式所示,且与模态仿真结果吻合较好[19]

      ${f_1} = \frac{1}{{2L}}\sqrt {\frac{{{E_0}}}{\rho }} $

      式中:L为战斗部长度,E0为战斗部壳体材料的弹性模量,ρ为战斗部壳体材料的密度。

      由于真实战斗部与理想的一维圆杆有差异,固有频率理论值会比实际值偏大。另外,飞行速度、侵彻阻力均会降低一阶轴向振动频率[20]。因此,取一修正系数k对固有频率理论值进行修正。结合大量实测数据的频谱分析结果,k一般取0.7~0.9,此时一阶轴向振动的固有频率可近似为:

      ${f_1} \approx \frac{k}{{2L}}\sqrt {\frac{{{E_0}}}{\rho }} $

      战斗部的一阶轴向振动频率已知后,便可借鉴机械振动理论进行战斗部弹性振动模型建模。将战斗部等效为具有阻尼特性的线性弹簧,质量集中于战斗部尾部,则一阶轴向振动模型可简化为单自由度弹簧-质量-阻尼系统,如图2所示,与基础运动引起的强迫振动类似[21]。图中,战斗部头部位移为x1,与刚体运动一致,战斗部尾部(即引信安装位置)位移为x2,战斗部等效刚度K和等效阻尼C分别如下式所示:

      图  2  单自由度弹簧-质量-阻尼系统

      Figure 2.  Single DOF spring-mass-damper system

      $K = M\omega _n^2 = M{\left( {2{\text{π}} {f_1}} \right)^2} = 4M{{\text{π}} ^2}f_1^2$

      $C = 2\xi {\omega _n}M = 2\xi \sqrt {MK} $

      式中,ωnf1对应的圆周固有频率,单位为rad/s,ξ为无量纲阻尼比。对于金属材料而言,阻尼比远小于1。

      根据图2可知,弹靶作用力经弹簧传递后再作用于战斗部,此时战斗部受两部分作用力,一是阻尼力,其方向与振动方向相反,大小与振动速度成正比,二是弹簧力,其方向与振动方向相反,大小与振动位移成正比。由于战斗部尾部对于头部有相对位移x2x1,因而作用于战斗部的弹簧力为K(x1x2),阻尼力为C(dx1−dx2)/dt

      以战斗部为研究对象可得到运动微分方程,如下式所示:

      $M\frac{{{{\rm d}^2}{x_2}}}{{{\rm d}{t^2}}} = C\left( {\frac{{{\rm d}{x_1}}}{{{\rm d}t}} - \frac{{{\rm d}{x_2}}}{{{\rm d}t}}} \right) + K\left( {{x_1} - {x_2}} \right)$

      两边同时减去Md2x1/dt2并移项、化简可得:

      $M\left( {\frac{{{{\rm d}^2}{x_2}}}{{{\rm d}{t^2}}} - \frac{{{{\rm d}^2}{x_1}}}{{{\rm d}{t^2}}}} \right) + C\left( {\frac{{{\rm d}{x_2}}}{{{\rm d}t}} - \frac{{{\rm d}{x_1}}}{{{\rm d}t}}} \right) + K\left( {{x_2} - {x_1}} \right) = - M\frac{{{{\rm d}^2}{x_1}}}{{{\rm d}{t^2}}} = {F_x}$

      记战斗部一阶轴向振动引起的长度变化量为ΔL,则ΔL=x2x1。当ΔL>0时,战斗部被压缩,反之被拉伸。此时,式(15)可简化为式(16),即一阶轴向振动模型,描述了战斗部轴向振动变形与弹靶作用力之间的关系。

      $M\frac{{{{\rm d}^2}\Delta L}}{{{\rm d}{t^2}}} + C\frac{{{\rm d}\Delta L}}{{{\rm d}t}} + K\Delta L = {F_x}$

    • 将式(9)和式(16)所示的微分方程相叠加,便可得到战斗部刚体运动与一阶轴向振动相结合的弹靶作用模型,如式(17)所示。可以看出,弹靶作用力决定了战斗部的刚体运动特性与一阶轴向振动特性。其中,刚体运动用战斗部头部位移x1及其一阶导数、二阶导数描述,一阶轴向振动用战斗部长度变化量ΔL及其一阶导数、二阶导数描述。

      $\left\{ \begin{aligned} &M\frac{{{{\rm d}^2}{x_1}}}{{{\rm d}{t^2}}} = - {F_x} \\ &M\frac{{{{\rm d}^2}\Delta L}}{{{\rm d}{t^2}}} + C\frac{{{\rm d}\Delta L}}{{{\rm d}t}} + K\Delta L = {F_x} \end{aligned} \right.$

      对于刚体运动,可直接根据积分求解出侵彻深度x1、速度dx1/dt、加速度d2x1/dt2,如下式所示:

      $\left\{ \begin{aligned} & \frac{{{{\rm d}^2}{x_1}}}{{{\rm d}{t^2}}} = - \frac{{{F_x}}}{M} \\ & \frac{{{\rm d}{x_1}}}{{{\rm d}t}} = {V_0} - \int_0^t {\frac{{{F_x}}}{M}{\rm d}t} \\ & {x_1} = {V_0}t - \int {\int_0^t {\frac{{{F_x}}}{M}{\rm d}{t^2}} } \end{aligned} \right.$

      式中:V0为初始速度。

      假设弹靶作用力Fx的加载、卸载非常快,可近似为峰值为Fmax的脉冲函数,则一阶轴向振动微分方程的解等效为二阶欠阻尼系统的脉冲响应。此时,战斗部长度变化量ΔL为:

      $\Delta L = \frac{{{F_{\max }}}}{{M{\omega _n}\sqrt {1 - {\xi ^2}} }}{{\rm e}^{ - \xi {\omega _n}t}}\sin \left( {{\omega _n}\sqrt {1 - {\xi ^2}} t} \right)$

      对其进行求导可得出变形速度、变形加速度。可以看出,ΔL是振荡衰减的,频率近似为ωn

      为了分析方便,记d2x1/dt2为战斗部刚体过载,该值仅与战斗部刚体运动有关;记d2ΔL/dt2为战斗部弹性过载,该值仅与战斗部的一阶轴向振动有关;记d2x2/dt2为战斗部总过载,并等效为引信设计输入,该值既与刚体运动有关,又与一阶轴向振动密切相关。此时,战斗部总过载与刚体过载、弹性过载之间的关系为:

      $\frac{{{{\rm d}^2}{x_2}}}{{{\rm d}{t^2}}} = \frac{{{{\rm d}^2}{x_1}}}{{{\rm d}{t^2}}} + \frac{{{{\rm d}^2}\Delta L}}{{{\rm d}{t^2}}}$

    • 由于侵彻过程中的弹靶作用力计算非常复杂,而且存在加载与卸载,无法等效为脉冲函数,因此直接求解式(17)的解析解是非常困难的。为此,本文采用数值积分的方法预测侵彻过程中的战斗部刚体运动特性和一阶轴向振动特性。

      首先,根据已知的战斗部参数、靶标参数、速度参数,采用空腔膨胀理论和微分面元法求解弹靶作用力;随后,以弹靶作用力为输入求解刚体运动微分方程,获得侵彻深度、侵彻速度、刚体过载等物理量,以弹靶作用力为输入求解一阶轴向振动微分方程,获得长度变化量、变形速度、弹性过载等物理量,并判断是否满足程序结束条件。如果不满足,则计算下一时刻的值直至满足程序结束条件。

      详细的程序流程如图3所示。

      图  3  程序流程图

      Figure 3.  Flow diagram for programming

    • 本文以某火炮试验弹侵彻钢筋混凝土靶标为例进行计算,详细的参数如下:

      试验弹弹径0.155 m,弹长0.89 m,弹头系数3.0,弹重约90 kg,无量纲阻尼比0.02;弹体材料为高强度合金钢,弹性模量210 GPa,密度7 800 kg/m3,一阶轴向振动的固有频率修正系数0.8。根据以上参数计算出试验弹的一阶轴向振动的固有频率约2.3 kHz。

      靶标分强度靶(single-layer concrete target)和多层靶(multi-layer concrete target)两种。其中,强度靶为3.5 m厚C40钢筋混凝土;多层靶为8层C40钢筋混凝土,首靶的厚度0.3 m,其余靶的厚度0.14 m,靶间距3 m。C40钢筋混凝土的抗压强度40 MPa,弹性模量32.5 GPa,密度2 400 kg/m3。采用空腔膨胀理论计算弹靶作用力时,动摩擦系数0.02,锁变体积常数0.1,强化模量参数0.65。

      数值积分算法选用四阶龙格库塔算法。计算时,时间步长1 μs,采样频率20 kHz,即每隔50 μs输出一次计算结果。

      取试验弹以830 m/s速度垂直侵彻强度靶和多层靶为典型工况,对每种工况分别采用传统方法(traditional method)和本文的修正方法(modificatory method)进行计算。其中,传统方法仅考虑战斗部刚体运动,修正方法综合考虑战斗部刚体运动与一阶轴向振动。

    • 侵彻强度靶时,持续时间(以弹尾完全出靶为准)约9 ms,剩余速度约355 m/s,侵彻多层靶时,持续时间(以弹尾完全出最后一层靶为准)约33 ms,剩余速度约645 m/s。

      垂直侵彻过程中,战斗部刚体过载变化曲线如图4所示。可以看出,无论是强度靶还是多层靶,刚体过载均是单向的,本质上与弹靶作用力的变化趋势一致。

      图  4  垂直侵彻过程中的刚体过载变化曲线

      Figure 4.  Rigid acceleration curve in normal penetration process

      垂直侵彻过程中,战斗部弹性过载变化曲线如图5所示。可以看出,无论是强度靶还是多层靶,弹性过载的变化规律是一致的,均表现出明显的周期性,而且由于阻尼耗散作用,幅值是衰减的,直到下一次碰靶。

      图  5  垂直侵彻过程中的弹性过载变化曲线

      Figure 5.  Elastic acceleration curve in normal penetration process

      对弹性过载进行频谱分析,结果如图6所示。可以看出,无论是强度靶还是多层靶,信号能量主要集中于2.3 kHz,与一阶轴向振动频率一致,这种现象是由战斗部的频率响应特性决定的。将式(16)进行拉普拉斯变换,可得到以刚体过载为输入、以弹性过载为输出的传递函数,如式(21)所示,对应的频率响应特性如图7所示。可以看出,侵彻过程中的战斗部实际上是一个自动控制系统,该系统响应外界输入时,信号的能量主要集中于固有频率f1附近,与图6结果一致。

      图  6  弹性过载的频谱分析结果

      Figure 6.  Frequency spectrum of elastic acceleration signal

      图  7  战斗部一阶轴向振动的频率响应特性

      Figure 7.  Frequency response characteristics of the first order axial vibration for projectiles

      $G\left( s \right) = \frac{{M{s^2}\Delta L\left( s \right)}}{{ - {F_x}\left( s \right)}} = \frac{{M{s^2}}}{{M{s^2} + Cs + K}} = \frac{{{s^2}}}{{{s^2} + 2\xi {\omega _n}s + \omega _n^2}}$

      垂直侵彻过程中,同时考虑战斗部刚体运动和一阶轴向振动的战斗部总过载变化曲线如图8所示,对应的频谱分析结果如图9所示。

      图  8  垂直侵彻过程中的总过载变化曲线

      Figure 8.  Total acceleration curve in normal penetration process

      图  9  战斗部总过载的频谱分析结果

      Figure 9.  Frequency spectrum of total acceleration signal

      可以看出,两种方法的差异主要体现在计算结果的频率成分不同。采用传统方法时,加速度信号是单向的,未叠加高频信号,而采用修正方法时,加速度信号是振荡的,甚至是正负交替变化的,含有明显的高频成分,且高频成分的能量主要集中于一阶轴向振动频率附近。因此,直接将传统方法计算出的加速度值作为引信受力是不合理的。

    • 为了验证计算结果的准确性,利用155火炮平台开展了侵彻8层钢筋混凝土靶试验,靶标布置和试验弹分别如图10图11所示。其中,靶标是倾斜80°安装的,与垂直侵彻略有区别。

      图  10  靶标布置

      Figure 10.  Schematic diagram of the target

      图  11  试验弹

      Figure 11.  Schematic diagram of the projectile

      试验时,通过高速摄像记录完整的穿靶图像,并判读出侵彻速度约833 m/s,通过弹载记录装置采集侵彻过程中的加速度信号。本文从两个方面对计算结果和试验结果进行了对比,一是每层靶的碰靶时刻(以碰首靶为零点),目的是验证刚体运动模型的准确性,二是加速度信号,目的是验证一阶轴向振动模型的准确性。

      每层靶的碰靶时刻对比如表1所示。表中,绝对误差等于计算结果减去高速摄像判读结果,相对误差等于绝对误差除以高速摄像判读结果。可以看出:(1)计算出的碰靶时刻要比高速摄像判读结果略有延迟,说明采用空腔膨胀理论和微分面元法计算出的侵彻阻力偏大;(2)绝对误差最大值仅为1.00 ms,发生在侵彻第6层靶时,相对误差最大值仅为5.97%,发生在侵彻第5层靶时,说明侵彻第5层靶时开始出现较为明显的弹体姿态偏转并持续至碰靶结束,与图10中弹道向上偏转的规律相似,导致高速摄像判读出的后三层的碰靶时刻延后,即后面三层的绝对误差和相对误差均减小。以上分析说明战斗部刚体运动模型的计算结果与试验结果一致,能满足工程估算的需要。

      层数计算结果/ms高速摄像判读结果/ms绝对误差/ms相对误差/%
      第 1 层 0.00 0.00.00
      第 2 层 4.20 4.00.205.00
      第 3 层 8.30 8.00.303.75
      第 4 层12.5012.00.504.17
      第 5 层16.8515.90.955.97
      第 6 层21.3020.31.004.92
      第 7 层25.9025.00.903.60
      第 8 层30.6529.80.852.85

      表 1  每层靶的碰靶时刻

      Table 1.  Penetration moment of every layer target

      弹载记录装置采集到的加速度信号如图12所示,对应的频谱分析结果如图13所示。从时域角度看,实测加速度信号是振荡的,变化规律与修正方法计算结果基本一致,从频域角度看,实测加速度信号含有明显的高频成分,且信号的能量主要集中于约2.5 kHz的频段,与试验弹的一阶轴向振动频率基本吻合。

      图  12  实测加速度信号

      Figure 12.  Acceleration data collected

      图  13  实测数据的频谱

      Figure 13.  Frequency spectrum of data collected

      综合对比分析时域特性和频域特性可以看出,修正方法的计算结果更接近于实测数据,说明战斗部刚体运动与一阶轴向振动相结合的弹靶作用模型能更合理地描述侵彻过程,能为侵彻引信抗高过载优化设计提供准确的力学输入。

    • (1)将侵彻战斗部等效为具有阻尼特性的线性弹簧,基于单自由度弹簧-质量-阻尼系统建立了战斗部一阶轴向振动模型,对传统的仅考虑战斗部刚体运动的垂直侵彻弹靶作用模型进行了修正;

      (2)采用数值积分的方法求解战斗部刚体运动模型和一阶轴向振动模型,获得了垂直侵彻过程中各物理量的变化规律,重点分析了战斗部过载的变化规律;

      (3)利用155火炮平台开展了侵彻8层钢筋混凝土靶试验,通过高速摄像获得了每层靶的碰靶时刻,通过弹载记录装置获得了侵彻过程中的加速度信号;

      (4)计算结果与实测数据的对比分析结果表明:考虑战斗部一阶轴向振动的弹靶作用模型能更准确地预测侵彻过程中的过载变化规律,能为侵彻引信抗高过载优化设计提供更明确的力学输入。

参考文献 (21)

目录

    /

    返回文章
    返回