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

椭圆截面截卵形刚性弹体正贯穿加筋板能量耗散分析

王浩 潘鑫 武海军 皮爱国 李金柱 黄风雷

引用本文:
Citation:

椭圆截面截卵形刚性弹体正贯穿加筋板能量耗散分析

    作者简介: 王 浩(1990- ),男,博士研究生,wangh@bit.edu.cn;
    通讯作者: 武海军, wuhj@bit.edu.cn
  • 中图分类号: O385

Energy dissipation analysis of elliptical truncated oval rigid projectilepenetrating stiffened plate

    Corresponding author: WU Haijun, wuhj@bit.edu.cn ;
  • CLC number: O385

  • 摘要: 为获得椭圆截面截卵形刚性弹体正贯穿加筋板的剩余速度,根据椭圆截面弹体贯穿靶板的破坏特征,认为贯穿过程中靶板的能量耗散方式主要为塞块剪切变形功与塞块动能、扩孔塑性变形功、花瓣动力功、花瓣弯曲变形功、靶板整体凹陷变形功、加强筋侧向凹陷变形功。推导了每种能量计算方法,计算中定量考虑了靶板扩孔、花瓣弯曲、凹陷变形的应变率效应。根据能量守恒关系,得到了椭圆截面弹体剩余速度和弹道极限速度预测公式。并通过实验结果对模型进行了验证。结果表明:考虑靶板应变硬化、应变率效应的贯穿模型可以准确预测弹体剩余速度;随着椭圆截面弹体长短轴之比的增大,靶板的弹道极限速度近似线性增大;长短轴之比小于3时,加筋板的主要耗能为花瓣弯曲变形能、整体凹陷变形能。
  • 图 1  着靶点位置

    Figure 1.  Impact locations

    图 2  椭圆截面弹体撞击单加筋板示意图

    Figure 2.  Schematic diagram of single stiffened plate impactedby elliptic section projectile

    图 3  加筋板变形破坏模式 [9]

    Figure 3.  Deformation and damage modes of stiffened plate

    图 4  靶板椭圆形扩孔变形示意图

    Figure 4.  Deformation diagram of elliptical hole enlargementof target plate

    图 5  截卵弹体贯穿薄板剖面

    Figure 5.  Plate perforated by truncated oval-nose projectiles

    图 6  椭圆形扩孔动力功计算原理图

    Figure 6.  Schematic diagram of dynamical workduring elliptical hole enlargement

    图 7  面板椭圆形整体凹陷变形计算原理图

    Figure 7.  Schematic diagram of elliptical dishing energy calculation

    图 8  加强筋侧向凹陷变形

    Figure 8.  Dishing deformation of stiffener

    图 9  无量纲弹道极限速度随${k_{{{ab}}}}$变化

    Figure 9.  Dimensionless ballistic limit velocity vs ${k_{{{ab}}}}$

    图 10  椭圆截面与“截面积等效”弹体各耗散能量比随${k_{{{ab}}}}$的变化

    Figure 10.  Variation of energy dissipation ratio of elliptical sectionand "section equivalent" projectiles with ${k_{{{ab}}}}$

    图 11  平板和加筋板各能量耗散比随${k_{{{ab}}}}$的变化

    Figure 11.  Energy dissipation ratio of flat and stiffened plates vs ${k_{{{ab}}}}$

    表 1  截卵形弹体几何参数

    Table 1.  Geometry parameters oftruncated oval-nose projectiles

    b1/mmb0/mmM/gL/mmCRH
    12.54165.327.51.42
    下载: 导出CSV

    表 2  加筋板结构和力学参数

    Table 2.  Geometry and mechanical parameters of stiffened plate

    hp/mmhs/mmbs/mmρ/(g·cm−3)σ0/MPaσy/MPaσu/MPaεfDP
    5,102557.85410~420510585~6100.21.14 ×1045.8
    下载: 导出CSV

    表 3  贯穿实验剩余速度与模型预测结果对比

    Table 3.  Comparison between the theoretical and experimental residual velocities

    实验板型mn${v_0}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$${v_{\rm{f}}}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$$\Delta r/{\rm{mm}}$${w_0}/{\rm{mm}}$${v_{\rm{f}}}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$η/%${v_{\rm{f}}}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$η/%
    实验参数刚塑性模型本文模型
    1F0024413315180.535.71135.21.65
    5SI00.524715612.510191.322.63151.82.69
    10SI*10388144 5252.375.21140.82.22
    13CS1142724012.520323.734.87258.27.85
    15CS1141715712.518311.898.60243.555.10
    16QS0024712816184.043.75139.99.30
    下载: 导出CSV
  • [1] 杜忠华, 曾国强, 余春祥, 等. 异型侵彻体垂直侵彻半无限靶板试验研究 [J]. 弹道学报, 2008, 20(1): 19–21.
    DU Zhonghua, ZENG Guoqiang, YU Chunxiang, et al. Experimental research of novel penetrator vertically penetrating semi- infinite target [J]. Journal of Ballistics, 2008, 20(1): 19–21.
    [2] 杜忠华, 朱建生, 王贤治, 等. 异型侵彻体垂直侵彻半无限靶板的分析模型 [J]. 兵工学报, 2009, 30(4): 403–407. DOI: 10.3321/j.issn:1000-1093.2009.04.005.
    DU Zhonghua, ZHU Jiansheng, WANG Xianzhi, et al. Analytical model on non-circular penetrator impacting semi-infinite target perpendicularly [J]. Acta Armamentarii, 2009, 30(4): 403–407. DOI: 10.3321/j.issn:1000-1093.2009.04.005.
    [3] 高光发, 李永池, 刘卫国, 等. 长杆弹截面形状对垂直侵彻深度的影响 [J]. 兵器材料科学与工程, 2011, 34(3): 5–8. DOI: 10.3969/j.issn.1004-244X.2011.03.002.
    GAO Guangfa, LI Yongchi, LIU Weiguo, et al. Influence of the cross-section shapes of long rod projectile on the vertical penetration depth [J]. Ordnance Material Science and Engineering, 2011, 34(3): 5–8. DOI: 10.3969/j.issn.1004-244X.2011.03.002.
    [4] BLESS S J, LITTLEFIELD D L, ANDERSON C E, et al. The penetration of non-circular cross-section penetrators [C] // Proceedings of the 15th International Symposium on Ballistics. Jerusalem, Israel: IBS, 1995: 21−24.
    [5] BLESS S J. Penetration mechanics of non-circular rods [C] // AIP Conference Proceedings. AIP, 1996: 1119−1122.
    [6] 王文杰, 张先锋, 邓佳杰, 等. 椭圆截面弹体侵彻砂浆靶规律分析 [J]. 爆炸与冲击, 2018, 38(1): 164–173. DOI: 10.11883/bzycj-2017-0020.
    WANG Wenjie, ZHANG Xianfeng, DENG jiajie, et al. Analysis of projectile penetrating into mortar target with elliptical cross-section [J]. Explosion and Shock Waves, 2018, 38(1): 164–173. DOI: 10.11883/bzycj-2017-0020.
    [7] LANDKOF B, GOLDSMITH W. Petalling of thin, metallic plates during penetration by cylindro-conical projectiles [J]. International Journal of Solids & Structures, 1985, 21(3): 245–266. DOI: 10.1016/0020-7683(85)90021-6.
    [8] 张中国, 黄风雷, 段卓平, 等. 弹体侵彻带加强筋结构靶的实验研究 [J]. 爆炸与冲击, 2004, 24(5): 431–436. DOI: 10.3321/j.issn:1001-1455.2004.05.009.
    ZHANG Zhongguo, HUANG Fenglei, DUAN Zhuoping, et al. The experimental research for projectile penetrating the structural target with rebar [J]. Explosion and Shock Waves, 2004, 24(5): 431–436. DOI: 10.3321/j.issn:1001-1455.2004.05.009.
    [9] CHEN Y, WANG Y, TANG P, et al. Impact characteristics of stiffened plates penetrated by sub-ordnance velocity projectiles [J]. Journal of Constructional Steel Research, 2008, 64(6): 634–643. DOI: 10.1016/j.jcsr.2007.12.006.
    [10] SONG W, NING J, WANG J. Normal impact of truncated oval-nosed projectiles on stiffened plates [J]. International Journal of Impact Engineering, 2008, 35(9): 1022–1034. DOI: 10.1016/j.ijimpeng.2007.05.008.
    [11] 徐双喜, 吴卫国, 李晓彬, 等. 截锥形弹穿甲单加筋板的破坏特性 [J]. 爆炸与冲击, 2011, 31(1): 62–68. DOI: 10.11883/1001-1455(2011)01-0062-07.
    XU Shuangxi, WU Weiguo, LI Xiaobin, et al. Falure characteristics of a conical projectile penetrating single stiffened plate [J]. Explosion and Shock Waves, 2011, 31(1): 62–68. DOI: 10.11883/1001-1455(2011)01-0062-07.
    [12] HE Q, XIE Z, XUAN H, et al. Ballistic testing and theoretical analysis for perforation mechanism of the fan casing and fragmentation of the released blade [J]. International Journal of Impact Engineering, 2016, 91: 80–93. DOI: 10.1016/j.ijimpeng.2016.01.001.
    [13] ZHAN T, LI J, LV S, et al. Residual velocity for the truncated ogival-nose projectile into stiffened plates [J]. Ships and Offshore Structures, 2015, 11(6): 636–644. DOI: 10.1080/17445302.2015.1041441.
    [14] 黄涛, 吴卫国, 李晓彬, 等. 截锥形弹体斜穿甲花瓣型破坏模型 [J]. 振动与冲击, 2010, 29(2): 125–127. DOI: 10.3969/j.issn.1000-3835.2010.02.028.
    HUANG Tao, WU Weiguo, LI Xiaobin, et al. Oblique armor-piercing effect of a truncated cylindro-conical projectile [J]. Journal of Vibration and Shock, 2010, 29(2): 125–127. DOI: 10.3969/j.issn.1000-3835.2010.02.028.
    [15] XU S X, WU W G, LI X B, et al. Petal failure characteristics of a conical projectile penetrating a thin plate at high oblique angle [J]. Journal of Shanghai Jiaotong University (Science), 2010, 15(4): 434–440. DOI: 10.1007/s12204-010-1029-8.
    [16] WU Q G, WEN H M. Petalling of a thin metal plate struck by a conical-nosed projectile [J]. Acta Mechanica Solida Sinica, 2015, 28(5): 568–577. DOI: 10.1016/S0894-9166(15)30050-1.
    [17] LEE Y W, WIERZBICKI T. Fracture prediction of thin plates under localized impulsive loading. Part I: dishing [J]. International Journal of Impact Engineering, 2005, 31(10): 1253–1276. DOI: 10.1016/j.ijimpeng.2004.07.010.
    [18] THOMSON W T. An approximate theory of armor penetration [J]. Journal of Applied Physics, 1955, 26(1): 80–82. DOI: 10.1063/1.1721868.
    [19] CORBETT G G, REID S R, JOHNSON W. Impact loading of plates and shells by free-flying projectiles: a review [J]. International Journal of Impact Engineering, 1996, 18(2): 141–230. DOI: 10.1016/0734-743X(95)00023-4.
    [20] JONES N. Structural impact [M]. 2nd ed. Cambridge: Cambridge University Press, 2011: 360−361. DOI: 10.1017/CBO9780511820625
  • [1] 秦健张振华 . 原型和模型不同材料时加筋板冲击动态响应的相似预报方法. 爆炸与冲击, 2010, 30(5): 511-516. doi: 10.11883/1001-1455(2010)05-0511-06
    [2] 王文杰张先锋邓佳杰郑应民刘闯 . 椭圆截面弹体侵彻砂浆靶规律分析. 爆炸与冲击, 2018, 38(1): 164-173. doi: 10.11883/bzycj-2017-0020
    [3] 徐双喜吴卫国李晓彬孔祥韶黄燕玲 . 截锥形弹穿甲单加筋板的破坏特性. 爆炸与冲击, 2011, 31(1): 62-68. doi: 10.11883/1001-1455(2011)01-0062-07
    [4] 朱锡牟金磊王恒张振华 . 水下爆炸载荷作用下加筋板的毁伤模式. 爆炸与冲击, 2010, 30(3): 225-231. doi: 10.11883/1001-1455(2010)03-0225-07
    [5] 牟金磊朱锡张振华王恒 . 水下爆炸载荷作用下加筋板的毁伤模式. 爆炸与冲击, 2009, 29(5): 457-462. doi: 10.11883/1001-1455(2009)05-0457-06
    [6] 侯海量朱锡古美邦 . 爆炸载荷作用下加筋板的失效模式分析及结构优化设计. 爆炸与冲击, 2007, 27(1): 26-33. doi: 10.11883/1001-1455(2007)01-0026-08
    [7] 严成欧卓成段卓平黄风雷 . 脆性材料动态强度应变率效应. 爆炸与冲击, 2011, 31(4): 423-427. doi: 10.11883/1001-1455(2011)04-0423-05
    [8] 汤铁钢李庆忠孙学林孙占峰金山谷岩 . 45钢柱壳膨胀断裂的应变率效应. 爆炸与冲击, 2006, 26(2): 129-133. doi: 10.11883/1001-1455(2006)02-0129-05
    [9] 席丰张云 . 脉冲载荷作用下钢梁动力响应及反常行为的应变率效应. 爆炸与冲击, 2012, 32(1): 34-42. doi: 10.11883/1001-1455(2012)01-0034-09
    [10] 罗鑫许金余白二雷李为民 . 碱类型对地聚合物混凝土应变率效应 影响的对比研究. 爆炸与冲击, 2014, 34(3): 340-346. doi: 10.11883/1001-1455(2014)03-0340-07
    [11] 惠旭龙白春玉刘小川牟让科王计真 . 宽应变率范围下2A16-T4铝合金动态力学性能. 爆炸与冲击, 2017, 37(5): 871-878. doi: 10.11883/1001-1455(2017)05-0871-08
    [12] 王振张超王银茂王祥索涛 . 飞机风挡无机玻璃在不同应变率下的力学行为. 爆炸与冲击, 2018, 38(2): 295-301. doi: 10.11883/bzycj-2016-0186
    [13] 胡亮亮黄瑞源高光发蒋东李永池 . 混凝土类材料SHPB实验中确定应变率的方法. 爆炸与冲击, 2019, 39(6): 063102-1-063102-9. doi: 10.11883/bzycj-2018-0142
    [14] 王峰王肖钧胡秀章刘文韬 . 卵形杆弹对铝靶的斜侵彻. 爆炸与冲击, 2005, 25(3): 265-270. doi: 10.11883/1001-1455(2005)03-0265-06
    [15] 陈国华胡昆周池楼祁帅 . 尖头碎片撞击小尺寸储罐的模拟实验. 爆炸与冲击, 2018, 38(6): 1295-1302. doi: 10.11883/bzycj-2017-0197
    [16] 文祝邱艳宇紫民赵章泳王明洋 . 钙质砂的准一维应变压缩试验研究. 爆炸与冲击, 2019, 39(3): 033101-1-033101-9. doi: 10.11883/bzycj-2018-0015
    [17] 尚兵盛精王宝珍胡时胜 . 不锈钢材料的动态力学性能及本构模型. 爆炸与冲击, 2008, 28(6): 527-531. doi: 10.11883/1001-1455(2008)06-0527-05
    [18] 吴金荣马芹永 . 聚酯纤维对透水沥青混凝土冲击压缩性能的影响. 爆炸与冲击, 2016, 36(2): 279-284. doi: 10.11883/1001-1455(2016)02-0279-06
    [19] 高光发郭扬波 . 高强混凝土动态压缩试验分析. 爆炸与冲击, 2019, 39(3): 033103-1-033103-10. doi: 10.11883/bzycj-2017-0405
    [20] 罗斌强胡时胜 . SnAgCu焊料的动态力学性能. 爆炸与冲击, 2009, 29(5): 542-545. doi: 10.11883/1001-1455(2009)05-0542-04
  • 加载中
图(11)表(3)
计量
  • 文章访问数:  48
  • HTML全文浏览量:  114
  • PDF下载量:  15
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-14
  • 录用日期:  2018-11-26
  • 网络出版日期:  2019-09-25
  • 刊出日期:  2019-10-01

椭圆截面截卵形刚性弹体正贯穿加筋板能量耗散分析

    作者简介:王 浩(1990- ),男,博士研究生,wangh@bit.edu.cn
    通讯作者: 武海军, wuhj@bit.edu.cn
  • 北京理工大学爆炸与科学国家重点实验室,北京 100081

摘要: 为获得椭圆截面截卵形刚性弹体正贯穿加筋板的剩余速度,根据椭圆截面弹体贯穿靶板的破坏特征,认为贯穿过程中靶板的能量耗散方式主要为塞块剪切变形功与塞块动能、扩孔塑性变形功、花瓣动力功、花瓣弯曲变形功、靶板整体凹陷变形功、加强筋侧向凹陷变形功。推导了每种能量计算方法,计算中定量考虑了靶板扩孔、花瓣弯曲、凹陷变形的应变率效应。根据能量守恒关系,得到了椭圆截面弹体剩余速度和弹道极限速度预测公式。并通过实验结果对模型进行了验证。结果表明:考虑靶板应变硬化、应变率效应的贯穿模型可以准确预测弹体剩余速度;随着椭圆截面弹体长短轴之比的增大,靶板的弹道极限速度近似线性增大;长短轴之比小于3时,加筋板的主要耗能为花瓣弯曲变形能、整体凹陷变形能。

English Abstract

  • 异型武器平台的发展对战斗部结构设计提出了新的需求。学者们探索了非圆截面弹体的侵彻性能。杜忠华等[1-2]、高光发等[3]、Bless等[4-5]研究了截面形状(三角形、矩形、圆形、十字形)对长杆弹侵彻半无限金属靶的影响。王文杰等[6]开展了椭圆截面弹体垂直侵彻半无限砂浆靶实验研究,获得了截面弹体长短轴变化对侵彻性能影响。以上研究主要是针对半无限靶,对于舰船类等薄靶目标,目前的研究仍集中在圆形截面弹体方面,关于非圆截面弹体的研究少有报道。Landkof等[7]在刚塑性模型假设的基础上,基于能量守恒原理建立了锥头弹丸正贯穿薄铝板弹道极限速度计算模型,将靶板的能量耗散分为裂纹传播、花瓣弯曲、靶板整体凹陷变形能。张中国等[8]、Chen等[9]开展了加筋板贯穿实验,根据实验数据修正了De Marre贯穿经验公式。Song等[10]、徐双喜等[11]同样基于刚塑性模型假设和能量守恒原理,分别对截卵形弹体和截锥形弹体贯穿加筋钢板进行了理论分析。He等[12]进行了平头柱形和刀型弹丸贯穿TC4平板和加筋板实验,将靶板的能量吸收分为整体凹陷变形能、冲塞变形能、花瓣撕裂和弯曲变形能、弹丸头部墩粗变形能,推导了弹道极限速度预测公式,计算中靶板的屈服强度使用了某个高应变率下靶板的动态强度。Zhan等[13]基于动量守恒定理,建立了弹体贯穿加筋板理论模型。黄涛等[14]、Xu等[15]推导了锥形弹体斜贯穿薄板花瓣形刚塑性破坏模型。Wu等[16]考虑了整体凹陷变形的应变率效应,建立了锥形弹贯穿薄板理论模型。以上模型对穿甲理论发展做出了重要贡献,但还存在一些不足之处:基于动量守恒的理论分析,忽略了靶板强度,低速下误差较大,即使高速下预测得到的速度也很难取得理想的结果[13]。以能量守恒为基础的理论分析或采用刚塑性模型[10-11],或采用某个高应变率下靶板的动态屈服强度代替静态屈服强度[12],或使用整体凹陷变形的应变率代替扩孔、弯曲变形等应变率[16],均缺乏对靶板不同变形模式下应变率效应的定量分析。

    本文中,对椭圆截面截卵形刚性弹体正贯穿加筋板进行理论分析,认为椭圆截面弹体与圆截面弹体贯穿薄板过程中能量耗散机理相同,其贯穿的能量损失可归为塞块剪切变形功与其动能、扩孔塑性变形功、花瓣动力功、花瓣弯曲变形功、靶板整体凹陷变形功、加强筋侧向凹陷变形功。并结合椭圆截面弹体的破坏特征,对每种能量计算方法进行理论推导。近似估计靶板塑性扩孔、花瓣弯曲的平均应变率,继而采用率相关的本构关系引入应变率效应的影响。根据能量守恒原理,建立椭圆截面弹体贯穿加筋板的弹道极限速度与剩余速度预测模型,且通过与实验数据比较,验证模型的准确性。讨论截面长短轴之比对弹道极限速度、靶板各能量耗散的影响。

    • 假设弹体为刚体,垂直撞击加筋板。加强筋截面为矩形。图1给出了弹体撞击十字加筋板着靶点位置。图2为椭圆截面弹体撞击矩形单加筋板示意图。其中,弹体头部截顶横截面长半轴和短半轴为${a_0}$${b_0}$,弹身横截面长半轴和短半轴为${a_1}$${b_1}$,则长短轴之比${k_{{{ab}}}} = {{{a_1}} / {{b_1}}}={{{{a_0}} / b}_0}$。弹体的质量为M,加筋板的面板厚度为${h_{\rm{p}}}$,加强筋的高度与宽度分别为${h_{\rm{s}}}$${b_{\rm{s}}}$。定义与椭圆截面弹体横截面积相等的圆柱形弹体为“截面积等效”弹体,其弹体半径$\bar r = \sqrt {{a_1}{b_1}} $。弹轴与弹体横截面长轴组成的剖面为弹体长轴剖面。弹体的初始撞击速度为${v_0}$。假设加筋板面板与加强筋材料相同,其密度为$\rho $,静态屈服强度为${\sigma _0}$。本文中,Johnson 破坏数${d_{\rm n}} = {{\rho v_0^2} / {{\sigma _0}}}$范围为${10^{-3}}\text{~}{10^{-1}}$,在此范围内靶板的响应以塑性变形为主。

      图  1  着靶点位置

      Figure 1.  Impact locations

      图  2  椭圆截面弹体撞击单加筋板示意图

      Figure 2.  Schematic diagram of single stiffened plate impactedby elliptic section projectile

      为考虑材料应变硬化的影响,基于塑性流动能量等效[17],修正材料的屈服强度为:

      ${\sigma _{\rm{y}}} = \frac{1}{{{\varepsilon _{\rm{f}}}}}\int_0^{{\varepsilon _{\rm{f}}}} {\sigma (\varepsilon )} {\rm{d}}\varepsilon $

      式中:$\sigma (\varepsilon )$表示材料随应变变化的应力值,${\varepsilon _{\rm{f}}}$为单轴拉伸状态下材料的断裂应变。

      为表征靶板在冲击条件下的应变率效应,材料的动态屈服强度表示为:

      ${\sigma _{\rm{d}}} = K{\sigma _{\rm{y}}}$

      式中:K为应变率增强系数,有Cowper-Symonds关系:

      $K = 1 + {\left( {\frac{{\dot \varepsilon }}{D}} \right)^{{1 / P}}}$

      式中:$\dot \varepsilon $为应变率,DP为相关参数,由K$\dot \varepsilon $关系曲线拟合得到。

      弹体垂直撞击加筋板中心位置时,弹头截锥面沿弹体运动方向剪切面板和加强筋形成塞块,塞块随弹体运动飞出。弹头的锥形部分侵入靶板,引起面板扩孔变形并形成花瓣形破坏(见图3(a)),加强筋绕弹肩位置翻转(见图3(b)),并在面板上形成弹身截面形状的孔洞。实验中,由于各种扰动的影响,出现更多的是弹体偏心撞击加强筋(弹体轴线位置偏离加强筋中心位置)的情形。由于弹体头部的挤压作用,加强筋会沿着弹体径向方向发生凹陷变形(见图3(c)),而不是发生剪切和花瓣弯曲变形。在弹体贯穿过程中,靶板整体同时伴随发生盘形凹陷变形。

      图  3  加筋板变形破坏模式 [9]

      Figure 3.  Deformation and damage modes of stiffened plate

    • 通过对截卵形弹体贯穿加筋板破坏过程和失效模式分析可知,当着靶位置为1或2时,将弹丸的能量损失分为塞块剪切变形功与塞块动能、扩孔塑性变形功、花瓣动力功、花瓣弯曲变形功、靶板整体凹陷变形功。弹体偏心正撞击加筋板时(着靶点位置3),能量损失增加强筋侧向凹陷变形功,且不再计及该加强筋的剪切和弯曲变形功。弹体仅撞击面板时(着靶位置4),忽略计算中关于加强筋的各项能量耗散。

    • 忽略塞块剪切变形中的应变率效应,塞块剪切变形功${W_{\rm{s}}}$由面板塞块剪切变形功${W_{{\rm{sp}}}}$与加强筋塞块剪切变形功${W_{{\rm{ss}}}}$组成:

      ${W_{\rm{s}}} = {W_{{\rm{sp}}}} + {W_{{\rm{ss}}}}$

      ${W_{{\rm{sp}}}} = [{\text{π}} {b_{\rm{0}}} + 2\left( {{a_{\rm{0}}} - {b_{\rm{0}}}} \right)]{\tau _{{m}}}{h_{\rm{p}}}^2$

      ${W_{{\rm{ss}}}} = m{b_{\rm{s}}}{\tau _{{m}}}{h_{\rm{s}}}^2\quad\quad 2{b_0} \text{>} {b_{\rm{s}}}$

      式中:${\tau _{{m}}} = {{{\sigma _{\rm{y}}}} / {\sqrt 3 }}$m为发生剪切和弯曲变形的加强筋数量。

    • 由Thomson理论[18]可知,靶板扩孔变形中满足下假设:(1)扩孔变形区域只有环向应力${\sigma _{{\theta }}}$比较显著,忽略径向应力${\sigma _{{r}}}$和法向应力${\sigma _{\textit{z}}}$;(2)靶板材料是不可压缩的或体积不变,即${\varepsilon _{{r}}} + {\varepsilon _{{\theta }}} + {\varepsilon _{\textit{z}}} = 0$;(3)忽略扩孔变形中径向拉伸变形,即${\varepsilon _{{r}}} = 0$图4为靶板椭圆形扩孔变形功计算原理图。靶板微元从初始位置因弹丸碰撞翻转产生的单位体积应变能为:

      图  4  靶板椭圆形扩孔变形示意图

      Figure 4.  Deformation diagram of elliptical hole enlargementof target plate

      $ {\rm{d}}W = \int {{\sigma _{{\theta }}}} {\rm{d}}{\varepsilon _{{\theta }}} = \int_s^{{R_{\rm{p}}}} {{\sigma _{{\theta }}}\frac{{{\rm{d}}s}}{s}} = {\sigma _{\rm{y}}}\ln \left(\frac{{{R_{\rm{p}}}}}{s}\right) $

      式中:${R_{\rm{p}}} = {{{b_{\rm{1}}}} / {\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}$e为椭圆的离心率。

      则靶板总的扩孔变形功为:

      $\begin{split}{E_{\rm{h}}} =& \int_0^{2{\text{π}} } {{\rm{d}}\alpha \int_{{R_{\rm{0}}}}^{{R_{\rm{p}}}} {({\rm{d}}W)} } {h_{\rm{p}}}s{\rm{d}}s =\\& \frac{{{h_{\rm{p}}}{\sigma _{\rm{y}}}{\text{π}} (b_{\rm{1}}^2 - b_{\rm{0}}^2 - 2b_{\rm{0}}^2\ln ({{{b_{\rm{1}}}} / {{b_{\rm{0}}}}}))}}{{2\sqrt {(1 - {e^2})} }}\end{split}$

      式中:${R_{\rm{0}}} = {{{b_{\rm{0}}}} / {\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}$

      图5为截卵弹体贯穿平板剖面图,由假设(3)可知,花瓣微元变形前后长度相等:

      图  5  截卵弹体贯穿薄板剖面

      Figure 5.  Plate perforated by truncated oval-nose projectiles

      $ {b_{\rm{p}}} - r = R(\theta - \omega ) $

      由几何关系可知:

      $\begin{aligned} \theta =& \arcsin \left( {1 - \frac{{d - 2{b_{\rm{p}}}}}{{2R}}} \right),\\ \omega =& \arcsin \left( {1 - \frac{{d - 2{r_{\rm{1}}}}}{{2R}}} \right)\end{aligned} $

      式中:${b_{\rm{p}}},{r_1} \in \left[ {0,{d / 2}} \right]$。为便于理论计算,在该区间内对反正弦函数进行线性拟合,得到:

      $\begin{aligned} \theta =& k{b_{\rm{p}}} + \arcsin (1 - {d / {2R}}),\\ \omega =& k{r_{\rm{1}}} + \arcsin (1 - {d / {2R}})\end{aligned} $

      式中:$k = \dfrac{1}{d}\left[ {{\text{π}} - 2\arcsin \left(1 - \dfrac{d}{{2R}}\right)} \right]$

      将式(11)代入式(9),由体积守恒得到:

      ${\varepsilon _{\rm{\theta }}} = \ln \left(\frac{{{r_{\rm{1}}}}}{r}\right) = \ln \left[ {\frac{{{b_{\rm{p}}}}}{r}\left(1 - \frac{1}{{Rk}}\right) + \frac{1}{{Rk}}} \right]$

      式(12)对时间求导,结合${b_{\rm{p}}} = \dfrac{d}{2}\sin \left(\dfrac{{\text{π}} }{2}\dfrac{{X'}}{L}\right)$得到:

      ${\dot \varepsilon _{{\theta }}} = \frac{{(Rk - 1)\dot X}}{{\dfrac{d}{2}\sin \left(\dfrac{{\text{π}} }{2}\dfrac{{X'}}{L}\right)(Rk - 1) + r}}\frac{{{\text{π}} d}}{{4L}}{\rm{cos}}\left(\frac{{\text{π}}}{2}\frac{{X'}}{L}\right)$

      取平均应变率为$r = d/4$处最大值的一半,则令$X' = 0,\dot X = {v_0},r = d/4$,可得:

      ${\bar {\dot \varepsilon} _{\rm{\theta }}} = \frac{{{{\dot \varepsilon }_{{{\theta {\rm{max}}}}}}}}{2} = \frac{{{\text{π}} (Rk - 1){v_0}}}{{2L}}$

      对于椭圆截面弹体dR均为$\alpha $的函数,为简化分析,可取长短轴的均值近似:

      $d = {a_{\rm{1}}} + {b_{\rm{1}}}$

      $R = {{({R_{\rm{1}}} + {R_{\rm{2}}})} / 2}$

      式中:${R_{\rm{1}}}$${R_{\rm{2}}}$为弹体长轴和短轴剖面头部弧长半径。

    • 花瓣动力功${W_{\rm{D}}}$是惯性力所做的功,包括面板花瓣扩孔动力功${W_{{\rm{Dp}}}}$和加强筋花瓣扩孔动力功${W_{{\rm{Ds}}}}$

      $ {W_{\rm{D}}} = {W_{{\rm{Dp}}}} + {W_{{\rm{Ds}}}} $

      图6为椭圆形扩孔动力功计算原理图。如图所示,当扩孔的短半轴为b时,面板微元(图中阴影部分)形成的花瓣质量为:

      图  6  椭圆形扩孔动力功计算原理图

      Figure 6.  Schematic diagram of dynamical workduring elliptical hole enlargement

      ${\rm{d}}{M_{\rm{p}}} = {h_{\rm{p}}}\rho \left( {\frac{{b_{\rm{p}}^2 - r_{\rm{0}}^2}}{{\rm{2}}}} \right){\rm{d}}\alpha $

      微元质量随时间的变化为:

      $\frac{{{\rm{d}}{M_{\rm{p}}}}}{{{\rm{d}}t}} = b_{\rm{p}}{h_{\rm{p}}}\rho {\rm{d}}\alpha \frac{{{\rm{d}}b_{\rm{p}}}}{{{\rm{d}}t}}$

      图5中的几何关系可知:

      $b_{\rm{p}} = \sqrt {{{\left\{ {d\left[ {{{\left(\frac{L}{d}\right)}^2} + \frac{1}{4}} \right]} \right\}}^2} - {{(L - X')}^2}} - \left\{ {d\left[ {{{\left(\frac{L}{d}\right)}^2} - \frac{1}{4}} \right]} \right\} \approx \frac{{2{b_{\rm{1}}}}}{{\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}\left(\frac{{X'}}{L} - \frac{{{{X'}^2}}}{{2{L^2}}}\right)$

      式中:$d = {{2{b_{\rm{1}}}} / {\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}$

      认为弹体贯穿薄板的瞬态过程为匀速运动,式(20)对时间求导可得:

      $\frac{{{\rm{d}}b_{\rm{p}}}}{{{\rm{d}}t}} = \frac{{2{b_1}}}{{\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}\left(\frac{{\dot X'}}{L} - \frac{{X'\dot X}}{{{L^2}}}\right)$

      $\begin{array}{*{20}{c}} {\dfrac{{{{\rm{d}}^2}b_{\rm{p}}}}{{{\rm{d}}{t^2}}} = } \end{array}\frac{{2{b_{\rm{1}}}}}{{\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}\left( - \frac{{{{\dot X}^2}}}{{{L^2}}}\right)$

      式中:$X' = X + L - H$

      面板微元的惯性力为:

      ${\rm{d}}{F_{\rm{p}}} = \frac{{\rm{d}}}{{{\rm{d}}t}}\left({\rm{d}}{M_{\rm{p}}}\frac{{{\rm{d}}b_{\rm{p}}}}{{{\rm{d}}t}}\right) = {\rm{d}}{M_{\rm{p}}}\frac{{{{\rm{d}}^2}b_{\rm{p}}}}{{{\rm{d}}{t^2}}} + \frac{{{\rm{d}}{M_{\rm{p}}}}}{{{\rm{d}}t}}\frac{{{\rm{d}}b_{\rm{p}}}}{{{\rm{d}}t}}$

      对式(23)积分,得到面板惯性力做功为:

      ${W_{{\rm{Dp}}}} = \int_0^{2{\text{π}} } {{\rm{d}}\alpha } \int_{\textstyle\frac{{{b_{\rm{0}}}}}{{\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}}^{\textstyle\frac{{{b_{\rm{1}}}}}{{\sqrt {1 - {e^2}{{\cos }^2}\alpha } }}} {{\rm{d}}{F_{\rm{p}}}} {\rm{d}}b_{\rm{p}} = \frac{{(2 - {e^2}){\text{π}} {h_{\rm{p}}}\rho }}{{3{{(1 - {e^2})}^{{3 / 2}}}}}{\left( {\frac{{{v_{\rm{0}}}}}{L}} \right)^2}{b_{\rm{1}}}(2{b_{\rm{0}}} + {b_{\rm{1}}}){({b_{\rm{0}}} - {b_{\rm{1}}})^2}$

      同理,加强筋花瓣扩孔动力功为:

      ${W_{{\rm{Ds}}}} = \int_{{b_{\rm{0}}}}^{{b_{\rm{1}}}} {\frac{{\rm{d}}}{{{\rm{d}}t}}\left( {{M_{\rm{s}}}\frac{{{\rm{d}}b}}{{{\rm{d}}t}}} \right)} {\rm{d}}b = 2{b_{\rm{s}}}{h_{\rm{s}}}\rho {b_{\rm{1}}}{\left( {{b_{\rm{1}}} - {b_{\rm{0}}}} \right)^2}{\left( {\frac{{{v_{\rm{0}}}}}{L}} \right)^2}$

    • 花瓣弯曲变形功${E_{\rm{b}}}$包括面板花瓣弯曲变形功${E_{{\rm{bp}}}}$和加强筋花瓣弯曲变形功${E_{{\rm{bs}}}}$

      ${E_{\rm{b}}} = {E_{{\rm{bp}}}} + {E_{{\rm{bs}}}}$

      忽略弯曲变形中塑性铰的移行过程,面板花瓣弯曲变形功为:

      ${E_{{\rm{bp}}}} = [{\text{π}} {b_{\rm{0}}} + 2\left( {{a_{\rm{1}}} - {b_{\rm{1}}}} \right)]{\text{π}} {M_{\rm{p}}}$

      式中:${M_{\rm{p}}} = {{{\sigma _{\rm{y}}}{h_{\rm{p}}}^2}/ 4}$

      加强筋花瓣弯曲变形功为:

      ${E_{{\rm{bs}}}} = m{\text{π}} {b_{\rm{s}}}{M_{\rm{s}}}$

      式中:${M_{\rm{s}}} = {{{\sigma _{\rm{y}}}{h_{\rm{s}}}^2} / 4}$

      由文献[19]可知,考虑应变率效应时,动态塑性极限弯矩${M'_{\rm{0}}}$表示为:

      $\frac{{{M'_{\rm{0}}}}}{{{M_{\rm{0}}}}} = 1 + \frac{{2q}}{{2q + 1}}{\left(\frac{{h\dot \kappa }}{{2D}}\right)^{{1 / q}}}$

      式中:$\dot \kappa $为曲率变化率,h为板的厚度。近似认为加强筋与面板弯曲变形的曲率变化率相同,其表达式为:

      $\dot \kappa = \frac{{\dot \psi }}{l}$

      式中:$\dot \psi $为弯曲变形角速度,$l$为花瓣的长度。

      图5中的几何关系,可知:

      $\dot \psi = \frac{{{{\dot r}_1}\left| {_{r = 0}} \right.}}{{{d / 2} - {r_{\rm{0}}}}}\sin \omega = \left( {1 - \frac{1}{{Rk}}} \right)\frac{{{\text{π}} d\dot X\sin \omega }}{{4L\left( {{d / 2} - {r_{\rm{0}}}} \right)}}\cos \left( {\frac{{{\text{π}} X'}}{{2L}}} \right)$

      $l = \frac{d}{2} - {r_0}$

      由于$\sin \omega \text{≤} 1$,令$X' = 0$$\dot X = {v_0}$,得到最大角加速度:

      ${\dot \psi _{{\rm{max}}}} \approx \left(1 - \frac{1}{{Rk}}\right)\frac{{{\text{π}} d{v_0}}}{{\left( {{d / {2 - {r_0}}}} \right)4L}}$

      取曲率平均变化率为其最大值的一半,于是可得:

      $\bar {\dot \kappa} = \frac{{{{\dot \kappa }_{{\rm{max}}}}}}{2} = \left(1 - \frac{1}{{Rk}}\right)\frac{{{\text{π}} d{v_0}}}{{8{{({d / 2} - {r_0})}^2}L}}$

      对于椭圆截面弹体花瓣形弯曲变形,其曲率平均变化率可参考式(15)~(16),进行近似估计。

    • 加筋板整体凹陷变形能Eg包括面板凹陷变形能Egp和加强筋凹陷变形能Egs

      $ {E_{\rm{g}}} = {E_{{\rm{gp}}}} + {E_{{\rm{gs}}}} $

      文献[12]给出了任意截面弹体贯穿靶板时凹陷变形能的计算方法。如图7所示,假设面板弹孔边界各处横向位移沿着曲率半径Or方向变化规律相同。

      图  7  面板椭圆形整体凹陷变形计算原理图

      Figure 7.  Schematic diagram of elliptical dishing energy calculation

      面板的横向位移可以表示为[20]

      $ w(r) = {w_{\rm{0}}}{{\rm e}^{ - {c_{\rm{1}}}\left( {r - {r_{\rm{p}}}} \right)}} $

      式中:${w_0}$为弹孔边界处横向位移,c1为相关系数,由靶板凹陷变形实验数据拟合得到。

      面板凹陷变形能为:

      $ {E_{{\rm{gp}}}} = \iint {\left( {{N_{{r}}}{\varepsilon _{{r}}} + {N_{\rm{\theta }}}{\varepsilon _{\rm{\theta }}} + {M_{{r}}}\left| {{\kappa _{{r}}}} \right| + {M_{\rm{\theta }}}\left| {{\kappa _{\rm{\theta }}}} \right|} \right)r{\rm{d}}r{\rm{d}}\theta } $

      式中:

      $\begin{aligned} {\varepsilon _{{r}}} = \frac{1}{2}{\left( {\frac{{{\rm{d}}w}}{{{\rm{d}}r}}} \right)^2},\quad\;{\varepsilon _{\rm{\theta }}} = 0,\\\;{\kappa _{{r}}} = - \frac{{{{\rm{d}}^2}w}}{{{\rm{d}}{r^2}}},\quad\;{\kappa _{\rm{\theta }}} = - \frac{1}{r}\frac{{{\rm{d}}w}}{{{\rm{d}}r}} \end{aligned}$

      $ {N_{{r}}} = {N_{{\theta }}} = {N_{\rm{0}}} = {\sigma _{\rm{y}}}h,\quad\;{M_{{r}}} = {M_{{\theta }}} = {M_{\rm{0}}} = \frac{{{\sigma _{\rm{y}}}{h^2}}}{4} $

      对面板,取$h = {h_{\rm{p}}}$,对加强筋,取$h = {h_{\rm{s}}}$

      图7中靶板微元凹陷变形能为:

      ${E_{{\rm{gp}}}}{\left| {{\rm{d}}A} \right._{{\rm{dishing}}}} = {\rm{d}}\theta \int_{{r_{\rm{p}}}}^{ + \infty } {\left[ {{N_0}\frac{1}{2}{{\left( {\frac{{{\rm{d}}w}}{{{\rm{d}}r}}} \right)}^2} + } \right.} \left. {{M_0}\left(\left| {\frac{{{{\rm{d}}^2}w}}{{{\rm{d}}{r^2}}} + \frac{1}{r}\frac{{{\rm{d}}w}}{{{\rm{d}}r}}} \right|\right)} \right]r{\rm{d}}r$

      对式(40)积分,得:

      ${E_{{\rm{gp}}}} = \int_A {{E_{{\rm{gp}}}}\left| {{\rm{d}}{A_{{\rm{dishing}}}}} \right.} = \left(\frac{1}{{\rm{4}}}{N_0}w_0 + {M_0}\right){w_0}{C_{\rm{p}}}{c_1} + 4{\text{π}} {w_0}\left({M_0} + \frac{1}{{16}}{N_0}w_0\right)$

      式中:${C_{\rm{p}}} = 2{\text{π}} {b_{\rm{1}}} + 4\left( {{a_1} - {b_1}} \right)$

      加强筋整体变形凹陷变形能为:

      ${E_{{\rm{gs}}}} = 2{b_{\rm{s}}}\int_{{b_1}}^{ + \infty } {\left[ {{N_0}\frac{1}{2}{{\left( {\frac{{{\rm{d}}w}}{{{\rm{d}}r}}} \right)}^2} + {M_0}\frac{{{{\rm{d}}^2}w}}{{{\rm{d}}{r^2}}}} \right]{\rm{d}}r} = 2{b_{\rm{s}}}\left(\frac{1}{{\rm{4}}}{N_0}w_0 + {M_0}\right){w_0}{c_1}$

      对于靶板整体弯曲变形,使用动态强度代替静态强度的方法近似考虑应变率效应,取其应变率与薄膜拉伸应变率相同。由文献[19],可知:

      $ {\dot \varepsilon _m} = \frac{{4{w_0}{v_0}}}{{3\sqrt 2 d'R'{{\ln }^2}({{d'}/ {2R'}})}} $

      式中:$R' \approx 5d' \approx 5({a_1} + {b_1})$

    • 加强筋的侧向凹陷变形能${E_{\rm{s}}}$包括侧向剪切变形能${E_{{\rm{sb}}}}$和侧向弯曲与拉伸变形能${E_{{\rm{ss}}}}$

      $ {E_{{\rm{sd}}}} = {E_{{\rm{sb}}}} + {E_{{\rm{ss}}}} $

      图8给出了椭圆截面弹体偏心撞击单加强筋侧向凹陷变形能计算原理图。为便于分析,计算过程中,忽略弹头挤压加强筋过程中头部曲率变化和弹体偏转的影响,由弹体穿靶后加强筋的残余变形求解其能量耗散。该计算方法得到是加强筋能量耗散的下限。如图8所示,建立坐标系$aO'w'$${w_{\rm{m}}}$为加强筋最大横向位移,则加强筋上任意点C的横向位移可表示为:

      图  8  加强筋侧向凹陷变形

      Figure 8.  Dishing deformation of stiffener

      ${w_{\rm{a}}} = {w_{\rm{m}}} - \left({b_1} - \sqrt {{{\left( {{{{b_1}} / {\sqrt {1 - {e^2}{{\cos }^2}\chi } }}} \right)}^2} - {a^2}} \right) \approx {w_{\rm{m}}} - \frac{{(1 - {e^2}){a^2}}}{{2{b_{\rm{1}}}}}$

      加强筋侧向凹陷变形能为:

      ${E_{{\rm{sb}}}} = 2n{h_{\rm{s}}}\int_0^{{L_{\rm{s}}}} {{\rm{(}}{N_{{\rm{0s}}}}{\varepsilon _{{r}}} + \left| {{M_{{\rm{0s}}}}{\kappa _{\rm{a}}}} \right|{\rm{)d}}a} = n{h_{\rm{s}}}{\sigma _{\rm{y}}}{b_{\rm{s}}}(2b_{\rm{1}}^2 + 3{b_{\rm{1}}}{b_{\rm{s}}} - 2\Delta {r^2})\frac{{\sqrt {(1 - {e^2})\left( {b_{\rm{1}}^2 - \Delta {r^2}} \right)} }}{{6b_{\rm{1}}^2}}$

      式中:n为发生侧向凹陷变形的加强筋数量,

      ${\varepsilon _r} = \frac{1}{2}{\left( {\frac{{{\rm{d}}{w_{\rm{a}}}}}{{{\rm{d}}a}}} \right)^2} = \frac{1}{2}{\left[ {\frac{{(1 - {e^2})a}}{{{b_1}}}} \right]^2},\quad\;{\kappa _a} = - \frac{{{{\rm{d}}^2}{w_{\rm{a}}}}}{{{\rm{d}}{a^2}}} = - \frac{{(1 - {e^2})}}{{{b_1}}}$

      ${N_{{\rm{0s}}}} = {\sigma _{\rm{y}}}{b_{\rm{s}}},\quad\;{M_{{\rm{0s}}}} = {{{\sigma _{\rm{y}}}{b_{\rm{s}}}^2} / 4}$

      ${L_{\rm{s}}} = \sqrt {\frac{{b_1^2 - \Delta {r^2}}}{{1 - {e^2}}}} $

      由文献[19]可知,加强筋侧向凹陷拉伸平均应变率近似为:

      ${\dot \varepsilon _{{\rm{s}}r}} \approx \frac{{2{b_1}{w_{\rm{m}}}{V_0}}}{{3\sqrt 2 LL_{\rm{s}}^2}}$

      加强筋的侧向剪切变形能为:

      ${E_{{\rm{ss}}}} = \frac{{n{L_{\rm{s}}}{b_{\rm{s}}}^2{\tau _{\rm{m}}}}}{2} = \frac{{n{b_{\rm{s}}}^2{\tau _{\rm{m}}}}}{2}\sqrt {\frac{{b_1^2 - \Delta {r^2}}}{{1 - {e^2}}}} $

    • 剩余动能${E_{\rm{r}}}$包括弹丸的剩余动能${E_{\rm{k}}}$、面板塞块动能${E_{\rm{p}}}$和加强筋塞块动能${E_{\rm{s}}}$。近似取塞块飞散速度与弹体剩余速度相同,有:

      ${E_{\rm{r}}} = {E_{\rm{k}}} + {E_{\rm{p}}} + {E_{\rm{s}}}$

      ${E_{\rm{k}}} = \frac{1}{2}Mv_{\rm{f}}^2$

      ${E_{\rm{p}}} = \frac{1}{2}{m_{\rm{p}}}v_{\rm{f}}^2 = \frac{1}{2}\rho {\text{π}} {a_{\rm{0}}}{b_{\rm{0}}}{h_{\rm{p}}}v_{\rm{f}}^2$

      ${E_{\rm{s}}} = \frac{1}{2}{m_{\rm{s}}}v_{\rm{f}}^2 = \left\{\begin{array}{*{20}{l}} {\rho {b_{\rm{s}}}{b_0}{h_{\rm{s}}}v_{\rm{f}}^2 }& {\text{着靶位置}}\; 1\\ \rho {b_{\rm{s}}}({a_0} + {b_0} - {{{b_{\rm{s}}}} / 2}){h_{\rm{s}}}v_{\rm{f}}^2 & {\text{着靶位置}}\; 2\\ 0 & {\text{着靶位置}} \;3,4 \end{array} \right. $

      式中:M${m_{\rm{p}}}$${m_{\rm{s}}}$分别为弹体质量、面板塞块质量、加强筋塞块质量,$v_{\rm{f}}$为弹体剩余速度。

    • 忽略贯穿过程中弹体变形能、裂纹传播能量、摩擦力做功,根据能量守恒定律得到:

      $\frac{1}{2}Mv_0^2 - {E_{\rm{r}}} = {W_{\rm{s}}} + {E_{\rm{h}}} + {W_{\rm{D}}} + {E_{\rm{b}}} + {E_{\rm{g}}} + {E_{{\rm{sd}}}}$

      将式(52)代入式(56),可得椭圆截面弹体的剩余速度:

      $v_{\rm{f}} = \sqrt {\frac{{Mv_0^2 - 2({W_{\rm{s}}} + {E_{\rm{h}}} + {W_{\rm{D}}} + {E_{\rm{b}}} + {E_{\rm{g}}} + {E_{{\rm{sd}}}})}}{{M + {m_{\rm{p}}} + {m_{\rm{s}}}}}} $

      $v_{\rm{f}} = 0$,得到弹道极限速度:

      $v_{{\rm{bl}}} = \sqrt {\frac{{2({W_{\rm{s}}} + {E_{\rm{h}}} + {W_{\rm{D}}} + {E_{\rm{b}}} + {E_{\rm{g}}} + {E_{\rm{s}}})}}{M}} $

    • 为校验本文模型的准确性,将模型预测结果与文献[9]中穿甲实验数据对比。实验中弹体材料为D6A,屈服强度为1 069 MPa。表1给出了截卵形弹体的几何参数。靶板材料为16MnR。表2给出了实验中平板和加筋板的结构和力学参数,其中DP由文献中靶板材料在不同应变率下的屈服强度值拟合得到。

      b1/mmb0/mmM/gL/mmCRH
      12.54165.327.51.42

      表 1  截卵形弹体几何参数

      Table 1.  Geometry parameters oftruncated oval-nose projectiles

      hp/mmhs/mmbs/mmρ/(g·cm−3)σ0/MPaσy/MPaσu/MPaεfDP
      5,102557.85410~420510585~6100.21.14 ×1045.8

      表 2  加筋板结构和力学参数

      Table 2.  Geometry and mechanical parameters of stiffened plate

      表3给出了截卵形弹体贯穿平板和加筋板实验剩余速度与理论模型计算结果对比。表中,F表示平板,SI表示单加筋板,CS表示十字加筋板,QS表示双十字加筋板。所有靶板面板(除SI*面板厚度为10 mm)尺寸均为500 mm×500 mm×5 mm,加强筋矩形横截面尺寸为5 mm×25 mm。由表3可以看出,刚塑性模型过高估计弹体剩余速度,本文的定量考虑应变率效应后,剩余速度预测结果与实验结果吻合良好。

      实验板型mn${v_0}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$${v_{\rm{f}}}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$$\Delta r/{\rm{mm}}$${w_0}/{\rm{mm}}$${v_{\rm{f}}}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$η/%${v_{\rm{f}}}/({\rm{m}} \cdot {{\rm{s}}^{ - 1}})$η/%
      实验参数刚塑性模型本文模型
      1F0024413315180.535.71135.21.65
      5SI00.524715612.510191.322.63151.82.69
      10SI*10388144 5252.375.21140.82.22
      13CS1142724012.520323.734.87258.27.85
      15CS1141715712.518311.898.60243.555.10
      16QS0024712816184.043.75139.99.30

      表 3  贯穿实验剩余速度与模型预测结果对比

      Table 3.  Comparison between the theoretical and experimental residual velocities

      计算过程需要做以下几点说明:(1)文献中未给出靶板贯穿后的最大横向位移,模型计算数据由文献中照片测量换算得到。实验表明${c_1}$${10^{-2}}\text{~}{10^{-3}}$量级,由式(41)~(42)可知其值对计算结果影响很小,本文计算中${c_1}$取0.005 mm−1。(2)在第15组实验中模型计算剩余速度误差较大,原因是弹体撞击的中心位置没有位于十字加筋中心,弹体受到非对称力作用,弹道发生偏转,由此导致了其他的能量损失。本文模型中未考虑弹体偏转的影响,得到的剩余速度应为实验的上限值。

    • 假设在本文的${k_{{{ab}}}}$范围内,椭圆截面弹体贯穿靶板的横向最大位移与“截面积等效”弹体相同。在${k_{{{ab}}}}$变化过程中,保持弹体的截面积、质量恒定,仅改变截面尺寸。加强筋不发生侧向凹陷变形。图9给出了弹体贯穿各靶板无量纲弹道极限速度随${k_{{{ab}}}}$变化规律。图中,${\bar v_{{\rm{bl}}}}$为“截面积等效”弹体贯穿平板的弹道极限速度。SI-1、SI-2分别表示弹体长轴剖面、短轴剖面与加强筋布置方向垂直。由图9可知,随着${k_{{{ab}}}}$增大,弹体贯穿3类靶板的弹道极限速度呈线性增大变化趋势。当${k_{{{ab}}}}$增大至5时,贯穿平板的弹道极限速度增加约15%,贯穿十字加筋板的弹道极限速度增加约为40%,贯穿单加筋板的增加量居于两者之间。对于圆截面弹体,单根加强筋增加的能量损失约使弹道极限速度提高30%。

      图  9  无量纲弹道极限速度随${k_{{{ab}}}}$变化

      Figure 9.  Dimensionless ballistic limit velocity vs ${k_{{{ab}}}}$

    • 图10给出了椭圆截面截卵形弹体与“截面积等效”弹体以弹道极限速度贯穿平板各耗散能量比值随${k_{{{ab}}}}$的变化。上划线表示“截面积等效”弹体各项能量耗散。由图10可知,塞块剪切和花瓣弯曲变形功随${k_{{{ab}}}}$近似线性增大。当${k_{{{ab}}}}$为2时,塞块剪切和花瓣弯曲变形功增加16%。扩孔塑性变形功与靶板凹陷变形功与${k_{{{ab}}}}$变化无关。当${k_{{{ab}}}}$从1增大到2时,花瓣动力功随其缓慢增加。当${k_{{{ab}}}}$为2时,增约25%。当${k_{{{ab}}}}$在2~5之间时,其值线性快速增加。图9中平板弹道极限速度随${k_{{{ab}}}}$增大,主要是由于弯曲变形功、花瓣动力功的增加。

      图  10  椭圆截面与“截面积等效”弹体各耗散能量比随${k_{{{ab}}}}$的变化

      Figure 10.  Variation of energy dissipation ratio of elliptical sectionand "section equivalent" projectiles with ${k_{{{ab}}}}$

      图11给出了弹体以弹道极限速度贯穿平板和加筋板时各能量耗散与弹丸初始动能(靶板耗散总能量)的比随${k_{{{ab}}}}$变化规律。图中,${w_{\rm s}}$表示塞块剪切变形功比,${e_{\rm h}}$表示扩孔塑性变形功比,${w_{\rm d}}$表示花瓣动力功比,${e_{\rm b}}$表示弯曲变形功比,${e_g}$表示靶板整体凹陷变形功比。

      图  11  平板和加筋板各能量耗散比随${k_{{{ab}}}}$的变化

      Figure 11.  Energy dissipation ratio of flat and stiffened plates vs ${k_{{{ab}}}}$

      由图11,可以得出以下规律:

      (1)扩孔塑性变形功比、整体凹陷变形功比随着${k_{{{ab}}}}$增大而减小。随着参与变形的加强筋长度的增加,弯曲变形功和塞块剪切变形功比由随着${k_{{{ab}}}}$增大逐渐增加变为减少。

      (2)对于平板,靶板的能量耗散主要为整体凹陷变形功,约占40%。其次为花瓣弯曲变形功,约占耗散总能量的25%。

      (3)对于单加筋板,当${k_{{{ab}}}}$<3时,靶板的能量耗散主要为整体凹陷变形功、花瓣弯曲变形功。当${k_{{{ab}}}}$>3时,碰撞条件为SI-2、CS时,花瓣动力功比例逐渐开始超过超过除弯曲变形功的其他能量。

      (4)对于十字加筋板,与单加筋板相比,塞块剪切变形功超过了整体整体凹陷变形功。

    • 通过对椭圆截面截卵形弹体贯穿加筋板能量耗散分析,根据能量守恒原理,建立了考虑应变硬化、应变率效应的椭圆截面弹体贯穿薄加筋板的弹道极限速度与剩余速度预测模型,模型预测结果与实验数据吻合较好。结合模型,讨论了椭圆截面长短轴变化对靶板弹道极限速度、各能量耗散的影响。得到如下主要结论。

      (1)本文在能量耗散计算中考虑了靶板应变硬化、不同变形模式下应变率效应影响,建立了椭圆截面弹体贯穿加筋板分析模型,可以准确预测弹体剩余速度,较刚塑性模型具有更高的精度。令${k_{{{ab}}}}$=1,可退化为圆截面弹体。

      (2)随着${k_{{{ab}}}}$的增大,靶板的弹道极限速度近似线性增大。

      (3)当${k_{{{ab}}}}$<3时,加筋板的能量耗散主要为整体凹陷变形功、花瓣弯曲变形功、剪切变形功。当${k_{{{ab}}}}$>3时,花瓣动力功比例增加,并逐渐超过整体凹陷变形功。对于十字加筋板,剪切变形功大于整体凹陷变形功。

参考文献 (20)

目录

    /

    返回文章
    返回