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

刚性弹体斜侵彻贯穿混凝土靶的姿态偏转理论模型

段卓平 李淑睿 马兆芳 欧卓成 黄风雷

引用本文:
Citation:

刚性弹体斜侵彻贯穿混凝土靶的姿态偏转理论模型

    作者简介: 段卓平(1965- ),男,博士,研究员,博导,duanzp@bit.edu.cn;
    通讯作者: 马兆芳, 603328@bit.edu.cn
  • 中图分类号: O385

Analytical model for attitude deflection of rigid projectile during oblique perforation of concrete targets

    Corresponding author: MA Zhaofang, 603328@bit.edu.cn ;
  • CLC number: O385

  • 摘要: 为描述刚性弹体斜侵彻贯穿混凝土靶的弹体姿态变化,针对已有贯穿模型存在的问题,在斜侵彻贯穿过程中考虑了弹体转动惯量对姿态偏转的影响,根据弹体贯穿靶板后的成孔特性重新假设了背靶面崩落块形状,并在弹体贯穿出靶的剪切冲塞阶段引入了弹体姿态二次偏转机制,从而建立了刚性弹体斜侵彻贯穿混凝土靶的姿态偏转理论模型,同时给出了混凝土薄靶、中厚靶和厚靶的分类方法。多种侵彻状态的理论模型计算结果均与实验测量结果吻合较好,表明本文理论模型可有效预估弹体斜侵彻贯穿混凝土靶的弹体出靶姿态。
  • 图 1  弹体斜侵彻贯穿混凝土靶过程中各阶段弹体姿态偏转过程示意图

    Figure 1.  Illustrated attitude deflection of the projectile in each stage during the oblique perforation of concrete targets

    图 2  弹体斜侵彻贯穿不同厚度混凝土靶时弹体姿态的偏转过程

    Figure 2.  Attitude deflection of the projectile in oblique perforation of concrete targets with different thickness.

    表 1  垂直侵彻普通强度(48 MPa)钢筋混凝土靶弹体出靶剩余速度的实验结果与计算结果[15]

    Table 1.  Experimental data and numerical results of the residual velocities in normal perforation of the normal strength (48 MPa) reinforced concrete targets[15]

    弹体质量/
    kg
    靶厚/
    m
    初始速度/
    (m·s−1)
    弹体出靶剩余速度/(m·s−1)相对误差/
    %
    实验结果[15]计算结果
    0.50.127 301 0 00
    360 67 7917.9
    38113615010.3
    43421425016.8
    606449482 7.3
    746605641 6.0
    749615645 4.9
    1 058947973 2.7
    下载: 导出CSV

    表 2  斜侵彻贯穿多层间隔混凝土薄靶实验中弹体出靶剩余速度的实验结果和计算结果

    Table 2.  Experimental data and numerical results of the residual velocities in oblique perforation of the multi-layered thin concrete targets

    实验序号靶厚/
    m
    初始速度/
    (m·s−1)
    弹体出靶剩余速度/(m·s−1)相对误差/
    %
    实验结果[10]计算结果
    1-10.308338208270.9
    1-20.188208008182.3
    1-30.188007717983.5
    2-10.306886846810.4
    2-20.186846696791.5
    下载: 导出CSV

    表 3  斜侵彻贯穿多层间隔混凝土薄靶实验中弹体出靶姿态偏转角的实验结果和计算结果

    Table 3.  Experimental data and numerical results of the attitude deflection angles in oblique perforation of the multi-layered thin concrete targets

    实验序号靶厚/
    m
    初始姿态角/
    (°)
    弹体出靶姿态偏转角/(°)绝对误差/
    (°)
    实验结果[10]计算结果
    1-10.3015.0−0.9−1.000.10
    1-20.1814.1−0.1−0.300.20
    1-30.1814.0−0.2−0.300.10
    2-10.3014.7−1.6−1.400.20
    2-20.1813.1−0.4−0.420.02
    下载: 导出CSV
  • [1] 王礼立. 应力波基础[M]. 2版. 北京: 国防工业出版社, 2005: 39−64.
    [2] 薛建峰, 沈培辉, 王晓鸣. 弹体斜侵彻混凝土靶的实验研究及其数值模拟 [J]. 爆炸与冲击, 2017, 37(3): 536–543. DOI: 10.11883/1001-1455(2017)03-0536-08.
    XUE Jianfeng, SHEN Peihui, WANG Xiaoming. Experimental study and numerical simulation of projectile obliquely penetrating into concrete target [J]. Explosion and Shock Waves, 2017, 37(3): 536–543. DOI: 10.11883/1001-1455(2017)03-0536-08.
    [3] 汪斌, 曹仁义, 谭多望. 大质量高速动能弹侵彻钢筋混凝土的实验研究 [J]. 爆炸与冲击, 2013, 33(1): 98–102. DOI: 10.11883/1001-1455(2013)01-0098-05.
    WANG Bin, CAO Renyi, TAN Duowang. Experimental study on penetration of reinforced concrete by a high-speed penetrator with large mass [J]. Explosion and Shock Waves, 2013, 33(1): 98–102. DOI: 10.11883/1001-1455(2013)01-0098-05.
    [4] 何翔, 徐翔云, 孙桂娟, 等. 弹体高速侵彻混凝土的效应实验 [J]. 爆炸与冲击, 2010, 30(1): 1–6. DOI: 10.11883/1001-1455(2010)01-0001-06.
    HE Xiang, XU Xiangyun, SUN Guijuan, et al. Experimental investigation on projectiles’ high-velocity penetration into concrete targets [J]. Explosion and Shock Waves, 2010, 30(1): 1–6. DOI: 10.11883/1001-1455(2010)01-0001-06.
    [5] 吕中杰, 徐钰巍, 黄风雷. 弹体斜侵彻混凝土过程中的方向偏转 [J]. 兵工学报, 2009, 30(2): 301–304.
    LYU Zhongjie, XU Yuwei, HUANG Fenglei. Transverse deflection of projectile obliquely penetrating into concrete [J]. Acta Armamentarii, 2009, 30(2): 301–304.
    [6] CHEN X W, FAN S C, LI Q M. Oblique and normal perforation of concrete targets by a rigid projectile [J]. International Journal of Impact Engineering, 2004, 30: 617–637. doi: 10.1016/j.ijimpeng.2003.08.003
    [7] CHEN X W, LI Q M. Deep penetration of a non-deformable projectile with different geometrical characteristics [J]. International Journal of Impact Engineering, 2002, 27(6): 619–637. doi: 10.1016/S0734-743X(02)00005-2
    [8] LI Q M, CHEN X W. Dimensionless formulae for penetration depth of concrete target impacted by a non-deformable projectile [J]. International Journal of Impact Engineering, 2003, 28: 93–116. doi: 10.1016/S0734-743X(02)00037-4
    [9] 马兆芳, 段卓平, 欧卓成, 等. 弹体斜侵彻多层间隔混凝土靶实验和数值模拟 [J]. 北京理工大学学报, 2016, 36(10): 1001–1005.
    MA Zhaofang, DUAN Zhuoping, OU Zhuocheng, et al. Experimental and simulative research on projectile oblique penetration into concrete targets with multi-layered space structure [J]. Transactions of Beijing Institute of Technology, 2016, 36(10): 1001–1005.
    [10] 马兆芳. 动能弹斜侵彻有限厚混凝土靶体的弹道规律研究 [D]. 北京: 北京理工大学, 2016: 12−26.
    [11] IPSON T W, RECHT R F. Ballistic penetration resistance and its measurement [J]. Experimental Mechanics, 1975, 15: 249–257. doi: 10.1007/BF02318057
    [12] 薛建锋, 沈培辉, 王晓鸣. 弹体斜侵彻混凝土靶面的开坑阶段分析 [J]. 南京理工大学学报, 2016, 40(1): 72–76.
    XUE Jianfeng, SHEN Peihui, WANG Xiaoming. Analysis on crater-forming of projectile obliquely penetrating into concrete target [J]. Journal of Nanjing University Science Technology, 2016, 40(1): 72–76.
    [13] FORRESTAL M J, ALTMAN B S, CARGILE J D, et al. An empirical equation for penetration depth of ogive-nose projectiles into concrete targets [J]. International Journal of Impact Engineering, 1994, 15(4): 395–405. doi: 10.1016/0734-743X(94)80024-4
    [14] DANCYGIER A N. Rear face damage of normal and high-strength concrete elements caused by hard projectile impact [J]. Aci Structural Journal, 1998, 95(3): 291–304.
    [15] HANCHAK S J, FORRESTAL M J, YOUNG E R, et al. Perforation of concrete slabs with 48 MPa (7 ksi) and 140 MPa (20 ksi) unconfined compressive strengths [J]. International Journal of Impact Engineering, 1992, 12(1): 1–7. doi: 10.1016/0734-743X(92)90282-X
  • [1] 薛建锋沈培辉王晓鸣 . 弹体斜侵彻混凝土靶的实验研究及其数值模拟. 爆炸与冲击, 2017, 37(3): 536-543. doi: 10.11883/1001-1455(2017)03-0536-08
    [2] 张欣欣武海军黄风雷皮爱国 . 斜侵彻混凝土靶的刻槽弹体的结构响应. 爆炸与冲击, 2019, 39(3): 033301-1-033301-6. doi: 10.11883/bzycj-2017-0047
    [3] 纪冲龙源万文乾邵鲁中 . 钢纤维混凝土抗侵彻与贯穿特性的实验研究. 爆炸与冲击, 2008, 28(2): 178-185. doi: 10.11883/1001-1455(2008)02-0178-08
    [4] 宋顺成才鸿年王富耻 . 关于弹体冲击和贯穿混凝土的三维数值模拟. 爆炸与冲击, 2006, 26(1): 1-6. doi: 10.11883/1001-1455(2006)01-0001-06
    [5] 孙传杰卢永刚张方举李会敏 . 新型头形弹体对混凝土的侵彻. 爆炸与冲击, 2010, 30(3): 269-275. doi: 10.11883/1001-1455(2010)03-0269-07
    [6] 岳松林王明洋张宁邱艳宇王德荣 . 混凝土板在接触爆炸作用下的震塌和贯穿临界厚度计算方法. 爆炸与冲击, 2016, 36(4): 472-482. doi: 10.11883/1001-1455(2016)04-0472-11
    [7] 张丁山王浩冯国增刘赟郭锦炎 . 钻地弹侵彻混凝土靶过程中弹体温度的变化. 爆炸与冲击, 2010, 30(3): 314-319. doi: 10.11883/1001-1455(2010)03-0314-06
    [8] 王一楠黄风雷段卓平 . 小攻角条件下动能弹体高速侵彻混凝土靶的弹体弯曲. 爆炸与冲击, 2010, 30(6): 598-606. doi: 10.11883/1001-1455(2010)06-0598-09
    [9] 陈伟王明洋顾雷雨 . 弹体在内摩擦介质中的斜侵彻深度计算. 爆炸与冲击, 2008, 28(6): 521-526. doi: 10.11883/1001-1455(2008)06-0521-06
    [10] 陈小伟张方举徐艾民屈明 . 细长薄壁弹体的屈曲和靶体等效分析. 爆炸与冲击, 2007, 27(4): 296-305. doi: 10.11883/1001-1455(2007)04-0296-10
    [11] 王德荣苏杭程怡豪冯淑芳 . 运用极限荷载法研究钢筋混凝土板低速侵彻效应. 爆炸与冲击, 2017, 37(5): 837-843. doi: 10.11883/1001-1455(2017)05-0837-07
    [12] 段建杨黔龙周刚王可慧张颖田亚军初哲 . 串联随进战斗部侵彻混凝土靶实验研究. 爆炸与冲击, 2007, 27(4): 364-379. doi: 10.11883/1001-1455(2007)04-0364-06
    [13] 王浩陶如意 . 截卵形弹头对混凝土靶侵彻性能的试验研究. 爆炸与冲击, 2005, 25(2): 171-175. doi: 10.11883/1001-1455(2005)02-0171-05
    [14] 宋梅利李文彬王晓鸣冯君刘志林 . 弹体高速侵彻效率的实验和量纲分析. 爆炸与冲击, 2016, 36(6): 752-758. doi: 10.11883/1001-1455(2016)06-0752-07
    [15] 赵军陈小伟金丰年徐迎 . 非对称质量磨蚀导致正侵彻弹体的弯曲屈服. 爆炸与冲击, 2011, 31(2): 119-126. doi: 10.11883/1001-1455(2011)02-0119-08
    [16] 谢文波张伟姜雄文 . 钢球斜侵彻碳纤维复合材料板的实验研究. 爆炸与冲击, 2018, 38(3): 647-653. doi: 10.11883/bzycj-2016-0289
    [17] 邓佳杰张先锋葛贤坤陈东东郭磊 . 基于局部相互作用理论的侵彻弹头部形状优化及仿真. 爆炸与冲击, 2017, 37(4): 611-620. doi: 10.11883/1001-1455(2017)04-0611-10
    [18] 邓勇军陈小伟姚勇杨涛 . 基于细观混凝土模型的刚性弹体正侵彻弹道偏转分析. 爆炸与冲击, 2017, 37(3): 377-386. doi: 10.11883/1001-1455(2017)03-0377-10
    [19] 葛涛刘保荣王明洋 . 弹体侵彻与贯穿有限厚度混凝土靶体的力学特性. 爆炸与冲击, 2010, 30(2): 159-163. doi: 10.11883/1001-1455(2010)02-0159-05
    [20] 皮爱国黄风雷 . 大长细比弹体斜侵彻混凝土靶的动力学响应. 爆炸与冲击, 2007, 27(4): 331-338. doi: 10.11883/1001-1455(2007)04-0331-08
  • 加载中
图(2)表(3)
计量
  • 文章访问数:  622
  • HTML全文浏览量:  334
  • PDF下载量:  35
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-22
  • 录用日期:  2018-12-29
  • 网络出版日期:  2019-05-25
  • 刊出日期:  2019-06-01

刚性弹体斜侵彻贯穿混凝土靶的姿态偏转理论模型

    作者简介:段卓平(1965- ),男,博士,研究员,博导,duanzp@bit.edu.cn
    通讯作者: 马兆芳, 603328@bit.edu.cn
  • 1. 北京理工大学爆炸科学与技术国家重点实验室,北京 100081
  • 2. 北京理工大学珠海学院,广东 珠海 519088

摘要: 为描述刚性弹体斜侵彻贯穿混凝土靶的弹体姿态变化,针对已有贯穿模型存在的问题,在斜侵彻贯穿过程中考虑了弹体转动惯量对姿态偏转的影响,根据弹体贯穿靶板后的成孔特性重新假设了背靶面崩落块形状,并在弹体贯穿出靶的剪切冲塞阶段引入了弹体姿态二次偏转机制,从而建立了刚性弹体斜侵彻贯穿混凝土靶的姿态偏转理论模型,同时给出了混凝土薄靶、中厚靶和厚靶的分类方法。多种侵彻状态的理论模型计算结果均与实验测量结果吻合较好,表明本文理论模型可有效预估弹体斜侵彻贯穿混凝土靶的弹体出靶姿态。

English Abstract

  • 弹体斜侵彻贯穿混凝土靶一直是侵彻领域的热点。实际情况下,由于弹体所处环境和受力情况非常复杂,弹体着靶时会与靶板法线存在一定角度,当混凝土靶板厚度有限时,弹体会贯穿靶板,同时反射形成的拉伸波会在背靶面造成一定厚度的崩落[1]。实验研究[2-5]发现,在斜侵彻过程中弹体姿态会发生一定程度的偏转,偏转程度受弹体着靶初速、着靶角等多种因素影响,为定量描述这些因素对弹体姿态偏转的影响,并预测斜侵彻贯穿混凝土靶时弹体姿态变化,需结合实验研究发展相应的理论模型。

    Chen等[6]建立了刚性弹体斜侵彻贯穿混凝土靶的三阶段理论模型,将刚性弹体斜侵彻贯穿混凝土靶的过程分为了初始开坑、隧道和剪切冲塞三个阶段,同时提出了斜锥形剪切冲塞块假设,并利用无量纲量冲击函数I和弹体形状函数N[7-8]分析了弹体的姿态偏转。然而该模型假定弹体仅在初始开坑阶段发生偏转,而在隧道阶段和剪切冲塞阶段弹体姿态不再改变,因此该模型得到的弹体姿态角(即弹体轴线与靶板表面外法线之间的夹角)的变化趋势是增大的。但是弹体斜侵彻贯穿混凝土薄靶实验所得弹体的出靶姿态角是减小的[9-10],由于开坑阶段弹体有跳弹趋势,弹体姿态角一定增大[11-12],因此在剪切冲塞阶段弹体姿态再次发生了偏转,则该模型中关于剪切冲塞段弹体姿态不变的假设与实际不符。同时,在弹体斜侵彻贯穿多层间隔混凝土靶实验中得到的靶板背面的破坏形状是以弹洞为圆心近似为对称圆形[10],并非Chen等[6]提出的不对称斜锥形弹坑。另外该模型仅根据开坑深度将混凝土靶分为薄靶和厚靶,忽略了这两者之间的中厚靶情况,从而限制了模型的适用范围,因此Chen等[6]的弹体斜侵彻贯穿混凝土靶的模型存在一定问题,需要进一步修正。

    本文中基于Chen等[6]的三阶段理论模型,根据实验现象重新假设冲塞块形状,且在弹体贯穿出靶的剪切冲塞阶段引入弹体姿态二次偏转机制,并考虑弹体转动惯量对姿态偏转的影响,通过分析不同阶段的弹体运动状态,从而建立适用范围更广、更符合实际的弹体斜侵彻贯穿混凝土靶的姿态偏转理论模型,同时提出明确的混凝土靶厚分类方法,通过计算结果与实验数据的对比验证本文模型预测弹体姿态变化的有效性。

    • 理论模型的建立基于以下假设条件:(1)弹体为刚体,在混凝土靶体内部做刚体平面运动;(2)混凝土靶体为均质各向同性材料,忽略其中钢筋和骨料结构的影响;(3)混凝土对弹体的阻力作用在弹头部。

      不考虑弹体攻角对弹体姿态的影响,仅关注侵彻过程中弹体姿态角β(弹体轴线与靶体表面外法线之间的夹角)的变化,并定义弹体姿态偏转角Δβ为侵彻过程中弹体瞬时姿态角与弹体初始姿态角之间的差值,Δβ的符号表示弹体姿态的偏转方向,为正时表示弹体姿态向着姿态角变大的方向偏转,为负时表示弹体姿态向着姿态角变小的方向偏转。

      设弹体斜侵彻贯穿混凝土靶时,弹体着靶初速为v0,初始姿态角为β0。如图1所示,弹体直径为d,弹长为L,弹体质心C距弹尖的距离为lC,弹头长度为h,混凝土靶体厚度为H。弹体的斜侵彻贯穿过程可分为初始开坑、隧道和剪切冲塞三个阶段[6],本文中认为弹体在初始开坑阶段和剪切冲塞阶段均会发生姿态偏转,且偏转方向相反,即在初始开坑阶段弹体姿态发生偏转且弹体姿态角β变大(Δβ值为正),隧道阶段弹体姿态不偏转(Δβ=0),剪切冲塞阶段弹体姿态再次偏转且弹体姿态角β变小(Δβ值为负)。

      图  1  弹体斜侵彻贯穿混凝土靶过程中各阶段弹体姿态偏转过程示意图

      Figure 1.  Illustrated attitude deflection of the projectile in each stage during the oblique perforation of concrete targets

    • 在初始开坑阶段,弹体的运动微分方程为:

      $ m{{{{\rm{d}}^2}{x_1}}}/{{{\rm{d}}{t^2}}} = - {F_{{\rm n}1}} \tag{1a} $

      ${{{{\rm{d}}^2}{y_1}}}/{{{\rm{d}}{t^2}}} = - {F_{{\rm t}1}} \tag{1b} $

      ${J_C}{{{{\rm{d}}^2}{\beta _1}}}/{{{\rm{d}}{t^2}}} = {M_C}\left( {{F_{{\rm t}1}}} \right) \tag{1c} $

      式中:弹体轴向合阻力Fn1、横向合阻力Ft1以及过质心C的转动力矩MC(Ft1)满足下式[6, 13]

      ${F_{{\rm n}1}} = {c_1}{x_1} \tag{2a} $

      ${F_{{\rm t}1}} = {F_{{\rm n}1}}\sin {\beta _0} \tag{2b} $

      ${M_C}\left( {{F_{{\rm t}1}}} \right) = {F_{{\rm t}1}}{l_C} \tag{2c} $

      式中:m为弹体质量,x1β1分别为弹体的瞬时侵彻行程和瞬时姿态角,JC为弹体绕垂直弹道平面且穿过质心C的转轴的转动惯量,c1为初始开坑阶段的阻力常数。开坑阶段开始瞬时t = 0,弹体侵彻行程X01 = x1(0) = 0,速度V01 = v1(0) = v0,角速度Ω01 = ω1(0) = 0,姿态角B01 = β1(0) = β0。联立式(1)~(2),由弹体速度v1 = dx1/dt、弹体角速度ω1 = dβ1/dt,利用前述初始条件,则可解得初始开坑阶段弹体瞬时侵彻行程x1、瞬时速度v1和瞬时姿态角β1的表达式。设初始开坑阶段结束瞬时t = T1,弹体侵彻行程X1 = x1(T1),速度V1 = v1(T1),弹体姿态角B1 = β1(T1),弹体姿态偏转角为Δβ1 = B1B01,则得关系式(阻力常数c1、时间T1需通过联立方程组求解):

      $ {X_1} = \frac{{{V_{01}}}}{{\sqrt {{{{c_1}}}/{m}} }}\sin \left( {\sqrt {{{{c_1}}}/{m}} {T_1}} \right)$

      ${V_1} = {V_{01}}\cos \left( {\sqrt {{{{c_1}}}/{m}} {T_1}} \right)$

      $\Delta {\beta _1} = \frac{{m{l_C}{V_{01}}}}{{{J_C}}}\sqrt {\frac{m}{{{c_1}}}} \sin {B_{01}}\left[ {\arccos \left( {\frac{{{V_1}}}{{{V_{01}}}}} \right) - \sqrt {1 - {{\left( {\frac{{{V_1}}}{{{V_{01}}}}} \right)}^2}} } \right]$

    • 隧道阶段开始瞬时t = T1,弹体侵彻行程X02 = x2(T1) = 0,弹体速度V02 = v2(T1) = V1,弹体姿态角B02 = B1。设隧道阶段结束瞬时t = T2,弹体侵彻行程X2 = x2(T2),弹体速度V2 = v2(T2),由于该阶段弹体横向阻力对称,弹体姿态不偏转,则弹体姿态角B2 = β2(T2) = B1 = β0β1,弹体姿态偏转角Δβ2 = 0。该阶段弹体轴向合阻力Fn2满足关系式[13]

      ${F_{n2}} = \frac{{{\rm{{\text{π}} }}{d^2}}}{4}\left( {S{f_C} + {N^*}\rho v_2^2} \right)$

      式中:S为经验常数,与混凝土无侧限抗压强度fC相关;N*为弹头形状系数,ρ为混凝土靶体密度。利用牛顿第二定律并对式(6)进行积分,可得隧道阶段弹体侵彻行程X2为:

      ${X_2} = \frac{{2 m}}{{{\rm{{\text{π}}}}{d^2}\rho {N^*}}}\ln \frac{{S{f_C} + {N^*}\rho V_1^2}}{{S{f_C} + {N^*}\rho V_2^2}}$

    • 将冲塞块形状近似假设为锥形圆台,并取冲塞块的坡度角α=66.1°[14],如图1所示。利用图1中的几何关系可得锥形圆台剪切冲塞块的侧面积As满足关系式:Ascosα = πdH*[sec(β0β1)+(H*/d)tanα],H*为剪切冲塞块厚度,因此剪切冲塞块形成瞬时弹体所受轴向合阻力Fns的表达式为[6]

      ${F_{\rm{ ns}}} = \frac{1}{{\sqrt 3 }}{f_C}{{{\text{π}}{ d}}}{H^*}\left[ {1 + \frac{{{H^*}}}{d}\tan \alpha \cos \left( {{\beta _0} + \Delta {\beta _1}} \right)} \right]$

      剪切冲塞阶段弹体姿态的偏转方向与初始开坑段相反,则该阶段弹体的运动微分方程为:

      $m\frac{{{{\rm{d}}^2}{x_3}}}{{{\rm{d}}{t^2}}} = - {F_{{\rm n}3}} \tag{9a} $

      $m\frac{{{{\rm{d}}^2}{y_3}}}{{{\rm{d}}{t^2}}} = {F_{{\rm t}3}} \tag{9b} $

      ${J_C}\frac{{{{\rm{d}}^2}{\beta _3}}}{{{\rm{d}}{t^2}}} = {M_C}\left( {{F_{{\rm t}3}}} \right) \tag{9c} $

      其中:

      ${F_{{\rm n}3}} = {F_{{\rm {ns}}}} - {c_2}{x_3} \tag{10a} $

      ${F_{{\rm t}3}} = {F_{{\rm n}3}}\sin \left( {{\beta _0} + \Delta {\beta _1}} \right) \tag{10b} $

      ${M_C}\left( {{F_{{\rm t}3}}} \right) = {F_{{\rm t}3}}\left( {{l_C} - {x_3}} \right) \tag{10c} $

      设剪切冲塞阶段随着弹体运动,弹体横向合阻力Ft3的作用位置由弹尖逐渐向弹身方向移动,轴向合阻力Fn3由该阶段开始瞬时的Fns逐渐线性递减为0,如式(10a)所示,其中c2为剪切冲塞段的阻力常数,且Ft3Fn3始终满足关系式(10b)[6]。设使弹体姿态偏转的力臂l随弹体侵彻行程x3线性变化l=lCx3,则转动力矩MC(Ft3)如式(10c)所示。

      设剪切冲塞阶段开始瞬时t = T2,弹体侵彻行程X03 = x3(T2) = 0,弹体速度V03 = v3(T2) = V2,弹体角速度Ω03 = ω3(T2) = 0,弹体姿态角B03 = β3(T2) = B2 = β0 + Δβ1。联立式(9)~(10),利用前述初始条件,可解得剪切冲塞段弹体瞬时侵彻行程x3、瞬时速度v3、瞬时加速度a3以及瞬时姿态角β3的表达式。

      当弹体运动到某一位置时,弹体与冲塞块分离,弹靶之间再无作用力(不考虑分离再接触),则此时弹体加速度为零。设此瞬时即为剪切冲塞阶段结束瞬时t = T3,且此时弹体侵彻行程X3 = x3(T3),速度V3 = v3(T3),加速度a3=0,角速度Ω3 = ω3(T3),姿态角B3 = β3(T3),弹体姿态偏转角Δβ3 = B3B03,则得关系式(阻力常数c2需通过联立方程组求解):

      $ {X_3} = \left( {{V_{03}}\sqrt {\frac{m}{{4{c_2}}}} - \frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}}} \right){{\rm{e}}^{\left( {{T_3} - {T_2}} \right)\sqrt {{{{c_2}}}/{m}} }} - \left( {\frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}} + {V_{03}}\sqrt {\frac{m}{{4{c_2}}}} } \right){{\rm{e}}^{ - \left( {{T_3} - {T_2}} \right)\sqrt {{{{c_2}}}/{m}} }} + \frac{{{F_{{\rm {ns}}}}}}{{{c_2}}}$

      ${V_3} = \sqrt {\frac{{{c_2}}}{m}} \left( {{V_{03}}\sqrt {\frac{m}{{4{c_2}}}} - \frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}}} \right){{\rm{e}}^{\left( {{T_3} - {T_2}} \right)\sqrt {{{{c_2}}}/{m}} }} + \sqrt {\frac{{{c_2}}}{m}} \left( {\frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}} + {V_{03}}\sqrt {\frac{m}{{4{c_2}}}} } \right){{\rm{e}}^{ - \left( {{T_3} - {T_2}} \right)\sqrt {{{{c_2}}}/{m}} }}$

      $\left( {{V_{03}}\sqrt {\frac{m}{{4{c_2}}}} - \frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}}} \right){{\rm{e}}^{\left( {{T_3} - {T_2}} \right)\sqrt {{{{c_2}}}/{m}} }} - \left( {\frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}} + {V_{03}}\sqrt {\frac{m}{{4{c_2}}}} } \right){{\rm{e}}^{ - \left( {{T_3} - {T_2}} \right)\sqrt {{{{c_2}}}/{m}} }} = 0$

      $\begin{split} \Delta {\beta _3} =& - \frac{{\sin {B_{03}}}}{{{J_C}}}\left[ {\frac{{{{mV_{03}^2}}/{2} - {F_{{\rm {ns}}}}{V_{03}}\sqrt {{m}/{{4{c_2}}}} }}{{{{{c_2}}}/{m}}}{{\rm{e}}^{\sqrt {{{{c_2}}}/{m}} \left( {{T_3} - {T_2}} \right)}} - \frac{{{{mV_{03}^2}}/{2} - {{F_{{\rm {ns}}}^2}}/{{2{c_2}}} + {F_{{\rm {ns}}}}{l_C}}}{{{{{c_2}}}/{m}}}} \right. - \\& \left. { \left( {\displaystyle\frac{{mV_{03}^2}}{2} - \frac{{F_{{\rm {ns}}}^2}}{{2{c_2}}}} \right)\frac{{{{\left( {{T_3} - {T_2}} \right)}^2}}}{2} + m{l_C}{V_{03}}\left( {{T_3} - {T_2}} \right)} \right] \\ \end{split} $

      ${{{\varOmega }}_3} = - \frac{{\sin {B_{03}}}}{{{J_C}}}\left[ {\left( {\frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}} - {l_C}} \right)\left( {m{V_{03}} - {F_{\rm{ns}}}\sqrt {\frac{m}{{{c_2}}}} } \right){{\rm{e}}^{\sqrt {{{{c_2}}}/{m}} \left( {{T_3} - {T_2}} \right)}} - \left( {mV_{03}^2 - \frac{{F_{{\rm {ns}}}^2}}{{{c_2}}}} \right)\frac{{\left( {{T_3} - {T_2}} \right)}}{2} + m{l_C}{V_{03}}} \right]$

    • 弹体斜侵彻贯穿不同厚度混凝土靶的姿态变化与靶体厚度有关,不同靶厚下弹体的姿态偏转过程如图2所示。

      图  2  弹体斜侵彻贯穿不同厚度混凝土靶时弹体姿态的偏转过程

      Figure 2.  Attitude deflection of the projectile in oblique perforation of concrete targets with different thickness.

    • 弹体斜侵彻贯穿混凝土厚靶时弹体姿态的偏转过程如图2(a)所示,包括完全开坑阶段、隧道阶段和剪切冲塞阶段[6]。由于隧道段足够长,能够完全阻滞弹体的横向运动,因此弹体姿态在剪切冲塞阶段不偏转(Δβ3=0),仅在开坑阶段发生一次偏转,则弹体出靶时姿态角为β0β1

      完全开坑时弹体的侵彻行程X1 = kd,其中k = 0.707 + h/d为无量纲侵彻经验深度[8]。由于完全开坑阶段结束瞬时即为隧道阶段开始瞬时,此时弹体轴向合阻力Fn1= Fn2,利用式(2a)与式(6),并联立式(3)~(5),通过求解方程组即可得到完全开坑阶段结束瞬时的T1、弹体速度V1和弹体姿态偏转角Δβ1以及开坑阶段阻力常数c1

      $V_1^2 = \frac{{mv_0^2 - \left( {{\rm{{\text{π}} }}k{d^3}/4} \right)S{f_C}}}{{m + \left( {{\rm{{\text{π}} }}k{d^3}/4} \right){N^*}\rho }}$

      ${T_1} = \frac{{\arccos \left( {{{{V_1}}}/{{{v_0}}}}\right)}}{{\sqrt {{{{c_1}}}/{m}} }}$

      $\Delta {\beta _1} = \frac{{m{l_C}{v_0}}}{{{J_C}}}\sqrt {\frac{m}{{{c_1}}}} \sin {\beta _0}\left[ {\arccos \left( {\frac{{{V_1}}}{{{v_0}}}} \right) - \sqrt {1 - {{\left( {\frac{{{V_1}}}{{{v_0}}}} \right)}^2}} } \right]$

      ${c_1} = \frac{{{\rm{{\text{π}}}}d}}{{4 K}}\left( {\frac{{{N^*}\rho v_0^2 + S{f_C}}}{{1 + \displaystyle\frac{{{\rm{{\text{π}} }}k{d^3}}}{{4 m}}{N^*}\rho }}} \right)$

      图2(a)可得弹体隧道区行程和剪切冲塞段行程分别满足关系式:

      ${X_2} = \left( {H - {H^*}} \right)\sec \left( {{\beta _0} + \Delta {\beta _1}} \right) - {X_1}$

      ${X_3} = {H^*}\sec \left( {{\beta _0} + \Delta {\beta _1}} \right)$

      由于隧道阶段结束瞬时即为剪切冲塞阶段开始瞬时,此时弹体轴向合阻力Fn2=Fn3,利用式(6)、(8),并联立式(7)、(11)~(13)与(20)~(21),将式(7)中的ln函数进行一阶Taylor展开,通过求解方程组即可得到弹体侵彻行程X2X3、弹体速度V2V3、剪切冲塞阶段结束瞬时的T3,以及剪切冲塞段的阻力常数c2和剪切冲塞块厚度H*

    • 弹体斜侵彻贯穿混凝土薄靶时弹体姿态的偏转过程如图2(b)所示,仅包括不完全开坑阶段和剪切冲塞阶段[6]、无隧道阶段,因此弹体姿态会在开坑阶段和剪切冲塞阶段发生两次偏转,且偏转方向相反,则弹体出靶时姿态角为β0β1β3

      图2(b)可得弹体初始开坑段和剪切冲塞段的侵彻行程X1X3分别满足关系式:

      ${X_1} = \left( {H - {H^*}} \right)\sec \left( {{\beta _0} + \Delta {\beta _1}} \right)$

      ${X_3} = {H^*}\sec \left( {{\beta _0} + \Delta {\beta _1} + \Delta {\beta _3}} \right)$

      由式(19)可知,阻力常数c1与开坑阶段的侵彻深度无关,因此式(19)也适用于不完全开坑。由于开坑阶段结束瞬时即为剪切冲塞阶段开始瞬时,此时弹体轴向合阻力Fn1=Fn3,利用式(2a)、(8),并联立式(3)~(5)、(11)~(15)与(22)~(23),求解方程组即可得到弹体不完全开坑阶段结束瞬时T1和剪切冲塞段结束瞬时T3、弹体侵彻行程X1X3、弹体速度V1V3、弹体姿态偏转角Δβ1和Δβ3、弹体角速度Ω3,以及剪切冲塞段阻力常数c2和剪切冲塞块厚度H*

    • 弹体斜侵彻贯穿混凝土中厚靶时弹体姿态的偏转过程如图2(c)所示,包括完全开坑阶段、隧道阶段和剪切冲塞阶段。但由于隧道区长度有限,不能完全阻滞弹体的横向运动,弹体姿态在剪切冲塞阶段仍会发生偏转,因此弹体出靶时姿态角为β0β1β3

      利用隧道阶段结束瞬时的弹体速度V2以及隧道段长度X2,假设该弹体以v0t的初速侵彻一虚拟薄靶,令其完全开坑段结束瞬时的弹体速度V1t=V2,则可利用式(16)反推得弹体初速v0t。引入修正系数δx描述弹体斜侵彻贯穿中厚混凝土靶时隧道长度对剪切冲塞段弹体姿态偏转角的影响,其表达式为:

      ${\delta _x} = \frac{{{v_{0t}}}}{{{v_0}}}\left( {1 - \frac{{{X_2}}}{{{X_{\rm{L}}}}}} \right)$

      式中:XL为临界隧道区长度,即能够完全阻滞弹体横向运动、使得弹体姿态在剪切冲塞阶段不发生偏转(即Δβ3=0)的最短隧道长度。若弹体长径比L/d>4,XL可取经验值XL=2d。利用式(14)~(15),则弹体斜侵彻贯穿混凝土中厚靶时剪切冲塞阶段结束瞬时的弹体姿态偏转角Δβ3和弹体角速度Ω3为:

      $\begin{split} {\rm{\Delta }}{\beta _3} = & - {\delta _x}\frac{{\sin {B_2}}}{{{J_C}}}\left[ {\frac{{{{mV_2^2}}/{2} - {F_{{\rm {ns}}}}{V_2}\sqrt {{m}/{{4{c_2}}}} }}{{{{{c_2}}}/{m}}}{{\rm{e}}^{\sqrt {{{{c_2}}}/{m}} \left( {{T_3} - {T_2}} \right)}} - \frac{{{{mV_2^2}}/{2} - {{F_{{\rm {ns}}}^2}}/{{2{c_2}}} + {F_{{\rm {ns}}}}{l_C}}}{{{{{c_2}}}/{m}}}} \right.-\\ &\left. { \left( {\frac{{mV_2^2}}{2} - \frac{{F_{{\rm {ns}}}^2}}{{2{c_2}}}} \right)\frac{{{{\left( {{T_3} - {T_2}} \right)}^2}}}{2} + m{l_C}{V_2}\left( {{T_3} - {T_2}} \right)} \right] \end{split}$

      ${{{\varOmega }}_3} = - {\delta _x}\frac{{\sin {B_2}}}{{{J_C}}}\left[ {\left( {\frac{{{F_{{\rm {ns}}}}}}{{2{c_2}}} - {l_C}} \right)\left( {m{V_2} - {F_{{\rm {ns}}}}\sqrt {\frac{m}{{{c_2}}}} } \right){{\rm{e}}^{\sqrt {{{{c_2}}}/{m}} \left( {{T_3} - {T_2}} \right)}} - \left( {mV_2^2 - \frac{{F_{{\rm {ns}}}^2}}{{{c_2}}}} \right)\frac{{\left( {{T_3} - {T_2}} \right)}}{2} + m{l_C}{V_2}} \right]$

    • 根据弹体侵彻过程中是否存在隧道段区分薄靶和非薄靶,设HL为临界靶厚,即弹体斜侵彻贯穿靶厚为HL的靶板时仅经历完全开坑和剪切冲塞两个阶段,无隧道阶段。因此当靶体厚度HHL时弹体侵彻过程不含隧道段,则靶体为薄靶;当HHL时弹体侵彻过程存在隧道段,则靶体为非薄靶。若弹体的几何形状、弹体和混凝土靶体的材料参数以及弹靶关系确定,则对应有唯一的临界靶厚HL,可参考上述分析过程通过求解方程组得到。

      由于弹体斜侵彻贯穿混凝土厚靶时隧道区足够长,弹体在剪切冲塞阶段不发生偏转;而斜侵彻贯穿中厚靶时隧道区长度有限,弹体在剪切冲塞段仍会再次偏转。因此本文中利用隧道区长度X2与临界隧道区长度XL区分混凝土中厚靶和厚靶,即当0 <X2XL时,靶体为中厚靶;当X2XL时,靶体为厚靶。

    • Hanchak等[15]利用实验手段测得了不同初速下弹体垂直侵彻48 MPa普通强度钢筋混凝土靶体的出靶剩余速度,其中混凝土靶体强度为48 MPa,混凝土密度为2 440 kg/m3,靶板厚度为0.127 m。实验所用弹体为尖卵形弹头,弹体质量为0.5 kg,弹体直径为0.025 4 m,CRH (caliber-radius-head)=3。如表1所示,本文中模型的计算结果与实验数据之间的相对误差均在18%以内,表明本文模型可以较准确地预估弹体垂直侵彻贯穿混凝土靶的出靶剩余速度。

      弹体质量/
      kg
      靶厚/
      m
      初始速度/
      (m·s−1)
      弹体出靶剩余速度/(m·s−1)相对误差/
      %
      实验结果[15]计算结果
      0.50.127 301 0 00
      360 67 7917.9
      38113615010.3
      43421425016.8
      606449482 7.3
      746605641 6.0
      749615645 4.9
      1 058947973 2.7

      表 1  垂直侵彻普通强度(48 MPa)钢筋混凝土靶弹体出靶剩余速度的实验结果与计算结果[15]

      Table 1.  Experimental data and numerical results of the residual velocities in normal perforation of the normal strength (48 MPa) reinforced concrete targets[15]

      在弹体斜侵彻贯穿多层间隔混凝土薄靶的实验研究中,马兆芳等[10]利用高速运动分析系统记录的弹体弹道偏转过程显示每层靶体的弹体出靶姿态角均小于入靶姿态角,即弹体姿态向着姿态角减小的方向发生了偏转,因此实验测得的弹体姿态偏转角(弹体出靶姿态角与入靶姿态角的差值)均为负值。该实验采用的钢筋混凝土靶的等效抗压强度为46 MPa,混凝土密度约2 500 kg/m3,第1层靶板厚度为0.30 m,其余各层厚度均为0.18 m。实验弹体为截卵形弹头,弹体质量为290 kg,弹体直径为0.25 m,弹体长径比为4.8,质心至弹头长度为0.68 m,转动惯量JC为35.32 kg·m2,CRH=1.56。本文仅对两发实验中无攻角影响的前两到三层靶体的弹体出靶剩余速度和姿态偏转角进行了计算,如表23所示,弹体出靶剩余速度的相对误差均在4%以内,弹体出靶姿态偏转角的绝对误差均在0.2°以内,表明本文模型引入的剪切冲塞阶段弹体姿态二次偏转机制是合理的,且能够准确预测弹体斜侵彻贯穿混凝土薄靶时弹体的出靶剩余速度和姿态偏转角。

      实验序号靶厚/
      m
      初始速度/
      (m·s−1)
      弹体出靶剩余速度/(m·s−1)相对误差/
      %
      实验结果[10]计算结果
      1-10.308338208270.9
      1-20.188208008182.3
      1-30.188007717983.5
      2-10.306886846810.4
      2-20.186846696791.5

      表 2  斜侵彻贯穿多层间隔混凝土薄靶实验中弹体出靶剩余速度的实验结果和计算结果

      Table 2.  Experimental data and numerical results of the residual velocities in oblique perforation of the multi-layered thin concrete targets

      实验序号靶厚/
      m
      初始姿态角/
      (°)
      弹体出靶姿态偏转角/(°)绝对误差/
      (°)
      实验结果[10]计算结果
      1-10.3015.0−0.9−1.000.10
      1-20.1814.1−0.1−0.300.20
      1-30.1814.0−0.2−0.300.10
      2-10.3014.7−1.6−1.400.20
      2-20.1813.1−0.4−0.420.02

      表 3  斜侵彻贯穿多层间隔混凝土薄靶实验中弹体出靶姿态偏转角的实验结果和计算结果

      Table 3.  Experimental data and numerical results of the attitude deflection angles in oblique perforation of the multi-layered thin concrete targets

    • (1)在弹体斜侵彻贯穿混凝土靶过程中,考虑了弹体转动惯量对弹体姿态的影响,并在剪切冲塞阶段引入了弹体姿态二次偏转机制,使弹体姿态在该阶段向着姿态角减小的方向偏转,同时将靶体背面产生的崩落块形状修正为对称锥形圆台,即令崩落块中心轴与靶体法线方向重合。计算结果表明本文模型可以准确预测弹体斜侵彻贯穿混凝土薄靶的出靶姿态偏转角,同时也说明剪切冲塞阶段弹体姿态二次偏转机制的引入和崩落块形状的修正是合理的。

      (2)通过定义临界靶厚HL与临界隧道区长度XL,给出了混凝土靶体厚度的定义区间:当靶厚小于HL时,靶体属于薄靶;当靶厚大于HL且隧道区长度大于或等于XL时,靶体属于厚靶;当靶厚大于HL且隧道区长度在0~XL之间时,靶体属于中厚靶。

      (3)刚性弹体斜侵彻贯穿不同厚度混凝土靶时弹体的姿态偏转由不同机制控制:若靶体为厚靶,则主要由初始开坑阶段弹体姿态一次偏转机制控制;若为薄靶,则由初始开坑阶段弹体姿态一次偏转机制和剪切冲塞阶段弹体姿态二次偏转机制共同控制;若靶体为中厚靶,则由初始开坑阶段弹体姿态一次偏转机制、隧道区的不完整夹持机制以及剪切冲塞阶段弹体姿态二次偏转机制共同控制。

参考文献 (15)

目录

    /

    返回文章
    返回