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

刚性平头弹正侵彻钢筋混凝土靶的阻力模型

刘永佑 杨华伟 张杰 王志勇 王志华

引用本文:
Citation:

刚性平头弹正侵彻钢筋混凝土靶的阻力模型

    作者简介: 刘永佑(1992- ),男,硕士,liuyongyou_cy@163.com;
    通讯作者: 王志华, wangzh@tyut.edu.cn
  • 中图分类号: O385

A resistance model for a rigid flat projectile penetrating a reinforced concrete target

    Corresponding author: WANG Zhihua, wangzh@tyut.edu.cn
  • CLC number: O385

  • 摘要: 本文基于素混凝土侵彻理论,将钢筋混凝土中钢筋的失效模式简化为弯曲剪切失效后,建立了刚性平头弹侵彻钢筋混凝土靶的阻力模型,侵彻深度的计算结果与Young公式吻合良好,表明本文提出的理论模型可较为合理地预测侵彻深度。进一步分析了不同着靶点的位置对弹体侵彻的影响,结果表明:当弹体直径与钢筋网眼尺寸的比值小于1时,弹体撞击到网眼中心处侵彻深度最大;当弹体直径与网眼尺寸的比值大于1时,最不利着靶点位置视其比值而定。最后,基于防护角度的最不利工况,建立了侵彻深度的工程计算公式。
  • 图 1  单根钢筋受力分析图

    Figure 1.  Stress analysis diagram of single steel bar

    图 2  弹体侵彻半无限厚靶过程示意图

    Figure 2.  Schematic diagram of a projectile penetrating a semi-infinite thick target

    图 3  弹体侵彻半无限厚靶算法示意图

    Figure 3.  Schematic diagram of a projectile penetrating a semi-infinite thick target

    图 4  不同着靶点位置

    Figure 4.  Different target positions

    图 5  侵彻深度随时间变化图(弹体直径为100 mm)

    Figure 5.  Penetration depth varing time at different target positions for the projectile with the diameter of 100 mm

    图 6  侵彻深度随时间变化图(弹体直径为140 mm)

    Figure 6.  Penetration depth varing time at different target positions for the projectile with the diameter of 140 mm

    表 1  计算侵彻深度

    Table 1.  Calculating penetration depth

    工况弹体参数弹体与钢筋的相对位置侵彻深度/mm误差/%
    直径/mm长度/mm质量/kg撞击速度/(m·s−1)
    18053213600钢筋交汇处 6490.6
    钢筋网格一边中点 6490.6
    钢筋网眼中心 6490.6
    Young公式 653
    29059918600钢筋交汇处 7122.47
    钢筋网格一边中点 7112.6
    钢筋网眼中心 7112.6
    Young公式 730
    310066625600钢筋交汇处 8013.84
    钢筋网格一边中点 8003.96
    钢筋网眼中心 8003.96
    Young公式 833
    410066625410钢筋交汇处 5137.57
    钢筋网格一边中点 5137.57
    钢筋网眼中心 5127.75
    Young公式 555
    510066625780钢筋交汇处1 0464.56
    钢筋网格一边中点1 0454.65
    钢筋网眼中心1 0454.65
    Young公式1 096
    下载: 导出CSV
  • [1] YOUNG C W. Depth prediction for earth-penetrating projectiles [J]. Journal of the Soil Mechanics and Foundations Division, 1969, 95: 803–818.
    [2] YOUNG C W. Penetration equations [J]. Office of Scientific and Technical Information Technical Reports, 1997, 33(1−12): 837–846.
    [3] LUK V K, FORRESTAL M J. Penetration into semi-infinite reinforced-concrete targets with spherical and ogival nose projectiles [J]. International Journal of Impact Engineering, 1987, 6(4): 291–301. DOI: 10.1016/0734-743X(87)90096-0.
    [4] 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.
    [5] 黄民荣, 顾晓辉, 高永宏. 刚性弹丸侵彻钢筋混凝土的实验和简化分析模型 [J]. 实验力学, 2009, 24(4): 283–290.
    HUANG M R, GU X H, GAO Y H. Experiment and simplified analytical model for penetration of rigid projectile in a reinforced concrete target [J]. Journal of Experimental Mechanics, 2009, 24(4): 283–290.
    [6] 欧阳春, 赵国志, 李文彬, 等. 弹丸垂直侵彻钢筋混凝土介质的工程模型 [J]. 弹箭与制导学报, 2004, 24(S7): 140–142.
    OUYANG C, ZHAO G Z, LI W B, et al. An engineering model for penetration normally into semi-infinite reinforced concrete targets [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2004, 24(S7): 140–142.
    [7] 周宁, 任辉启, 沈兆武, 等. 弹丸侵彻钢筋混凝土的工程解析模型 [J]. 爆炸与冲击, 2007, 27(6): 529–534. DOI: 10.11883/1001-1455(2007)06-0529-06.
    ZHOU N, REN H Q, SHEN Z W, et al. An engineering analytical model for projectiles to penetrate into semi-infinite reinforced concrete targets [J]. Explosion and Shock Waves, 2007, 27(6): 529–534. DOI: 10.11883/1001-1455(2007)06-0529-06.
    [8] 穆朝民, 任辉启. 弹丸对混凝土中钢筋结构侵彻效应研究 [J]. 兵器材料科学与工程, 2011, 34(6): 1–5.
    MU C M, REN H Q. Penetration effect of projectile into reinforced concrete targets [J]. Ordnance Material Science and Engineering, 2011, 34(6): 1–5.
    [9] 刘志林, 孙巍巍, 王晓鸣, 等. 卵形弹丸垂直侵彻钢筋混凝土靶的工程解析模型 [J]. 弹道学报, 2015(3): 84–90. DOI: 10.3969/j.issn.1004-499X.2015.03.016.
    LIU Z L, SUN W W, WANG X M, et al. Engineering analytical model of ogive-nose steel projectiles vertically penetrating reinforced concrete target [J]. Journal of Ballistics, 2015(3): 84–90. DOI: 10.3969/j.issn.1004-499X.2015.03.016.
    [10] 杨华伟.尖卵形长杆弹侵彻半无限混凝土靶的动力学行为研究[D].太原: 太原理工大学, 2018: 1−2.
    [11] 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(1): 93–116. DOI: 10.1016/S0734-743X(02)00037-4.
    [12] 欧阳春, 赵国志, 杜中华, 等. 弹丸垂直侵彻钢筋混凝土介质的工程解析模型 [J]. 爆炸与冲击, 2004, 24(3): 273–277.
    OUYANG C, ZHAO G Z, DU Z H, et al. An engineering analyticaI modeI for projectiles to penetrate normally into semi-infinite reinforced concrete targets [J]. Explosion and Shock Waves, 2004, 24(3): 273–277.
    [13] YU T X, CHEN F L. A further study of plastic shear failure of impulsively loaded clamped beams [J]. International Journal of Impact Engineering, 2000, 24(6): 613–629.
  • [1] 周宁任辉启沈兆武何祥刘瑞朝吴飚 . 弹丸侵彻钢筋混凝土的工程解析模型. 爆炸与冲击, 2007, 27(6): 529-534. doi: 10.11883/1001-1455(2007)06-0529-06
    [2] 楼建锋王政朱建士张凤国洪滔 . 含筋率和弹着点对钢筋混凝土抗侵彻性能的影响. 爆炸与冲击, 2010, 30(2): 178-182. doi: 10.11883/1001-1455(2010)02-0178-05
    [3] 邓勇军宋文杰陈小伟姚勇 . 钢筋混凝土靶侵彻的可压缩弹-塑性动态空腔膨胀阻力模型. 爆炸与冲击, 2018, 38(5): 1023-1030. doi: 10.11883/bzycj-2017-0043
    [4] 马天宝武珺宁建国 . 弹体高速侵彻钢筋混凝土的实验与数值模拟研究. 爆炸与冲击, 2019, 39(10): 103301-1-103301-11. doi: 10.11883/bzycj-2018-0275
    [5] 邓勇军陈小伟钟卫洲何丽灵 . 弹体正侵彻钢筋混凝土靶的试验及数值模拟研究. 爆炸与冲击, 2020, 40(2): 023101-1-023101-11. doi: 10.11883/bzycj-2019-0001
    [6] 温丽晶张春明郭超段璞张连生段卓平 . 飞机模型高速撞击钢筋混凝土载荷特性实验研究. 爆炸与冲击, 2018, 38(4): 811-819. doi: 10.11883/bzycj-2016-0337
    [7] 高飞王明洋张先锋何勇李猛深 . 强动载荷下钢筋混凝土结构计算模型简评. 爆炸与冲击, 2017, 37(2): 365-376. doi: 10.11883/1001-1455(2017)02-0365-12
    [8] 张德海朱浮声邢纪波 . 梁-颗粒模型在混凝土侵彻问题中的应用. 爆炸与冲击, 2005, 25(1): 85-89. doi: 10.11883/1001-1455(2005)01-0085-05
    [9] 吴成沈晓军王晓鸣姚文进 . 细观混凝土靶抗侵彻数值模拟及侵彻深度模型. 爆炸与冲击, 2018, 38(6): 1364-1371. doi: 10.11883/bzycj-2017-0123
    [10] 李德聪陈力丁雁生 . 装药弹体侵彻混凝土厚靶中的炸药摩擦起爆模型. 爆炸与冲击, 2009, 29(1): 13-17. doi: 10.11883/1001-1455(2009)01-0013-05
    [11] 林华令丁育青汤文辉 . 混凝土侵彻数值模拟的影响因素. 爆炸与冲击, 2013, 33(4): 425-429. doi: 10.11883/1001-1455(2013)04-0425-05
    [12] 屈明陈小伟 . 钢筋混凝土穿甲的数值模拟. 爆炸与冲击, 2008, 28(4): 341-349. doi: 10.11883/1001-1455(2008)04-0341-09
    [13] 何翔徐翔云孙桂娟沈俊杨建超金栋梁 . 弹体高速侵彻混凝土的效应实验. 爆炸与冲击, 2010, 30(1): 1-6. doi: 10.11883/1001-1455(2010)01-0001-06
    [14] 徐伟芳张方举陈裕泽毛勇建任时成陶俊林 . 混凝土薄板的侵彻破坏响应实验. 爆炸与冲击, 2013, 33(2): 169-174. doi: 10.11883/1001-1455(2013)02-0169-06
    [15] 王元博王肖钧卞梁余育苗 . CDM模型及在纤维增强层合材料侵彻数值模拟中的应用. 爆炸与冲击, 2008, 28(2): 172-177. doi: 10.11883/1001-1455(2008)02-0172-06
    [16] 王明洋郑大亮白晓燕 . 弹体对钢筋混凝土板-钢板的贯穿计算问题. 爆炸与冲击, 2005, 25(4): 289-295. doi: 10.11883/1001-1455(2005)04-0289-07
    [17] 张凤国刘军梁龙河楼建锋王政 . 数值建模时骨料对混凝土侵彻及毁伤问题的影响. 爆炸与冲击, 2013, 33(2): 217-221. doi: 10.11883/1001-1455(2013)02-0217-04
    [18] 赵小龙马铁华徐鹏范锦彪 . 弹丸侵彻混凝土加速度信号测试及分析. 爆炸与冲击, 2014, 34(3): 347-353. doi: 10.11883/1001-1455(2014)03-0347-07
    [19] 黄风雷张雷雷段卓平 . 大锥角药型罩聚能装药侵彻混凝土实验研究. 爆炸与冲击, 2008, 28(1): 17-22. doi: 10.11883/1001-1455(2008)01-0017-06
    [20] 皮爱国黄风雷 . 大长细比弹体斜侵彻混凝土靶的动力学响应. 爆炸与冲击, 2007, 27(4): 331-338. doi: 10.11883/1001-1455(2007)04-0331-08
  • 加载中
图(6)表(1)
计量
  • 文章访问数:  9
  • HTML全文浏览量:  31
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-12
  • 录用日期:  2019-05-16
  • 网络出版日期:  2019-12-25
  • 刊出日期:  2020-03-01

刚性平头弹正侵彻钢筋混凝土靶的阻力模型

    作者简介:刘永佑(1992- ),男,硕士,liuyongyou_cy@163.com
    通讯作者: 王志华, wangzh@tyut.edu.cn
  • 1. 太原理工大学应用力学研究所,山西 太原 030024
  • 2. 材料强度与结构冲击山西省重点实验室,山西 太原 030024
  • 3. 青海大学土木工程学院,青海 西宁 810016

摘要: 本文基于素混凝土侵彻理论,将钢筋混凝土中钢筋的失效模式简化为弯曲剪切失效后,建立了刚性平头弹侵彻钢筋混凝土靶的阻力模型,侵彻深度的计算结果与Young公式吻合良好,表明本文提出的理论模型可较为合理地预测侵彻深度。进一步分析了不同着靶点的位置对弹体侵彻的影响,结果表明:当弹体直径与钢筋网眼尺寸的比值小于1时,弹体撞击到网眼中心处侵彻深度最大;当弹体直径与网眼尺寸的比值大于1时,最不利着靶点位置视其比值而定。最后,基于防护角度的最不利工况,建立了侵彻深度的工程计算公式。

English Abstract

  • 钢筋混凝土因其优异的性能在军事和民用建筑领域得到了广泛应用。在现代战争中,为了更好地摧毁目标或发挥防护功能,对于钢筋混凝土在侵彻打击下动态响应的研究具有十分重要的意义。

    目前,国内外学者对于弹体侵彻素混凝土的研究已相对成熟。钢筋混凝土方面,Young[1-2]在经历了数千次实验后,提出了弹体侵彻岩石、土、混凝土、钢筋混凝土的统一经验公式,可较好地预测侵彻深度,但受限于经验模型本身的特点,无法给出侵彻阻力等中间变量和运动参量随时间变化的关系。Luk等[3]、Forrestal等[4]认为侵彻过程中钢筋只起到阻碍混凝土裂纹开展的作用,忽略了其提供的直接阻力,并基于此建立了侵彻模型,但其计算结果与实验数据相差20%以上。黄民荣等[5]通过实验发现钢筋的破坏模式有两种:弯曲拉伸破坏和弯曲剪切破坏,并据此建立了两种失效模式下弹体侵彻钢筋混凝土的阻力模型,然而只适用于弹体低速撞击靶体,导致模型应用受限。欧阳春等[6]、周宁等[7]、穆朝民等[8]在空腔膨胀理论的基础上,考虑了钢筋对弹体侵彻提供的直接阻力,认为钢筋发生弯曲变形破坏,建立了阻力模型;刘志林等[9]认为钢筋发生弯曲拉伸耦合断裂失效,分析其对弹体的阻力作用。此外,他们注意到侵彻后弹体接触到的钢筋发生断裂,未接触到的钢筋则完好无损。

    大量研究表明,根据不同速度条件下弹体响应的机制,可将侵彻划分为4个阶段[10],即:(1)刚性弹侵彻阶段(侵彻速度小于800 m/s);(2)变形弹侵彻阶段(侵彻速度介于800~1 500 m/s);(3)半流体侵彻阶段(侵彻速度介于1 500 ~3 000 m/s);(4)流体侵彻阶段(侵彻速度大于3 000 m/s)。本文中将研究撞击速度为800 m/s以下的刚性弹侵彻问题,并从钢筋弯曲剪切破坏方面分析,且基于以下假设:(1)钢筋与混凝土的粘结为理想粘结;(2)忽略钢筋与弹体之间、钢筋与混凝土之间的摩擦力;(3)弹体与钢筋接触时,只在接触点发生相互作用;(4)弹体侵彻钢筋混凝土靶体时,受力均匀且弹道不发生偏转。以钢筋网眼中心、钢筋交汇处、钢筋网格一边中点等3个特殊位置进行研究,建立弹体侵彻钢筋混凝土靶体的阻力模型,将计算结果与Young公式进行对比,结果吻合良好,表明本文中提出的理论模型能够较为合理地预测侵彻深度及其他侵彻过程量并给出弹体各运动参量随时间变化关系,进一步得到最不利的撞击位置及在此条件下的侵彻深度计算公式。

    • 弹体侵彻钢筋混凝土厚靶的弹道轨迹包含了两部分,即在入靶面形成的锥形前坑区域以及与弹体等直径的隧道区域。对于平头弹,Li等[11]利用滑移线理论得出前坑深度为Hrc=0.707d0,其中d0为弹体直径。

    • 单根钢筋的受力分析如图1所示(取一半为研究对象),AH为塑性铰位置,$\lambda $为塑性铰AH间的水平距离,q为混凝土提供给钢筋的压力,近似简化为均布荷载,MpQp为钢筋塑性极限弯矩、剪力,Fz为弹体轴向侵彻阻力。

      图  1  单根钢筋受力分析图

      Figure 1.  Stress analysis diagram of single steel bar

      设钢筋具有有限的剪切强度Qp,并忽略在极限状态时弯矩M和剪力Q之间的相互作用,则在塑性铰A处有MA=MpQA=Qp

      在小变形阶段,对AH段列竖直方向的动量定理:

      $\frac{\rm d}{{{\rm d}{t_2}}}\left(\frac{1}{2}\rho {\lambda _i}{v_{i0}}\right) = {Q_{\rm{p}}} - q{\lambda _i}$

      式中:$\rho $为钢筋线密度;λi为弹体撞击到第i层钢筋时钢筋上塑性铰AH间的水平距离;$v_{i{\rm{0}}}^{}$为弹体撞击到第i层钢筋到钢筋失效间的钢筋截面A的任意时刻速度;${Q_{\rm{p}}} = \dfrac{{{\text{π}} d_1^2{\sigma _{\rm{s}}}}}{{4\sqrt 3 }}$为钢筋极限剪力,$q = {d_1}{Y_{\rm{c}}}\left( {2 +{\text{π}} } \right)$[12]为混凝土对钢筋的阻力,d1为钢筋直径,σs为钢筋屈服应力,Yc为混凝土抗剪强度;t2为以弹体撞击到第i层钢筋为起始点的时间。

      弹体和与其直接接触部分钢筋的总质量$m' = {m_0} + {d_0}\rho $m0为弹体质量),由动量定理得:

      $\frac{\rm d}{{{\rm d}{t_2}}}(m'{v_{i1}}) = - 2{Q_{\rm{p}}} - q{d_0} - {F_{\textit{zi}}}$

      式中:${v_{i{\rm{1}}}}$为弹体撞击到第i层钢筋到钢筋失效间的弹体任意时刻速度;${F_{{{\textit{z}}i}}}$为弹体撞击到第i层钢筋时弹体受到的轴向侵彻阻力。

      将式(1)、(2)对时间积分得:

      $\frac{1}{2}\rho {\lambda _i}{v_{i0}} = {Q_{\rm{p}}}{t_2} - q{\lambda _i}{t_2}$

      $ 2{Q_{\rm{p}}}{t_2} + q{d_0}{t_2} + {F_{\textit{zi}}}{t_2} = m'{v_i} - m'{v_{i1}} $

      式中:${v_i}$为弹体第i次碰到钢筋时的速度。

      AH段绕H点作转动,列动量矩定理:

      $\frac{\rm d}{{{\rm d}{t_2}}}\left(\frac{1}{6}\rho {\lambda _i^2}{v_{i0}}\right) = 2{M_{\rm{p}}} - q\frac{{{\lambda _i^2}}}{2}$

      式中:为钢筋塑性极限弯矩。

      式(5)对时间积分得:

      $\frac{1}{6}\rho {\lambda _i^2}{v_{i{\rm{0}}}} = 2{M_{\rm{p}}}{t_2} - \frac{{q{\lambda _i^2}}}{2}{t_2}$

      联立式(3)和(6)得:

      $q{\lambda _i^2} + 2{Q_{\rm{p}}}{\lambda _i} - 12{M_{\rm{p}}} = 0$

      则:

      ${\lambda _i} = - \frac{{{Q_{\rm{p}}}}}{q} + \sqrt {{{\left( {\frac{{{Q_{\rm{p}}}}}{q}} \right)}^2} + 12\frac{{{M_{\rm{p}}}}}{q}} $

      由于塑性极限弯矩Mp和剪切强度Qp、混凝土对钢筋的阻力q都是常数,式(8)表明,在此阶段塑性铰H是一个驻定铰。

      由式(3)得塑性铰A所在的钢筋端部的横向速度:

      $\mathop v\nolimits_{i{\rm{0}}} = \frac{{2{Q_{\rm{p}}}{t_2} - 2q{\lambda _i}{t_2}}}{{\rho {\lambda _i}}}$

      弹体的速度与塑性铰A所在的钢筋端部的横向速度并不相同,由式(4)可知:

      ${v_{i{\rm{1}}}} = {v_i} - \frac{{2{Q_{\rm{p}}}{t_2} + q{d_0}{t_2} + {F_{{{\textit{z}}i}}}{t_2}}}{{m'}}$

      加速度为:

      ${\tilde a_i} = - \frac{{2{Q_{\rm{p}}} + q{d_0} + {F_{{{\textit{z}}i}}}}}{{m'}}$

      联立式(3)、(4)、(9),截面A两侧的相对速度为:

      ${v_{i1}} - {v_{i0}} = {v_i} - \left( {1 + \frac{{\rho {\lambda _i}}}{{m'}}} \right)\left( {\frac{{2{Q_{\rm{p}}}{t_2} - 2q{\lambda _i}{t_2}}}{{\rho {\lambda _i}}}} \right) - \frac{{2q{\lambda _i}{t_2} + q{d_0}{t_2} + {F_{{{\textit{z}}i}}}{t_2}}}{{m'}}$

      将相对速度对时间积分可得A截面两侧的相对位移:

      ${S_{1i}} = \int\limits_0^{{t_2}} {\left( {\mathop v\nolimits_{i{\rm{1}}} - \mathop v\nolimits_{i0} } \right){\rm d}{t_2}} = {v_i}{t_2} - \left(1 + \frac{{\rho {\lambda _i}}}{{m'}}\right)\left(\frac{{{Q_{\rm{p}}}t_2^2 - q{\lambda _i}t_2^2}}{{\rho {\lambda _i}}}\right) - \frac{{2q{\lambda _i}t_2^2 + q{d_0}t_2^2 + {F_{{{\textit{z}}i}}}t_2^2}}{{2m'}}$

      根据余同希等[13]的研究成果,S1i的极限值取为0.3d1且当钢筋剪切失效时的弹体侵彻位移:

      ${S_i} = {v_i}{t_2} - \frac{{2{Q_{\rm{p}}}t_2^2 + q{d_0}t_2^2 + {F_{{{\textit{z}}i}}}t_2^2}}{{2m'}}$

    • 由1.1.1节知,当弹体撞击到一根钢筋时,塑性铰H是驻定铰,则当弹体撞击到n根钢筋时,每根钢筋的塑性铰位置都为${\lambda _i}$

      弹体和与其直接接触部分的n根钢筋的总质量为$m = {m_0} + \displaystyle\sum\limits_{j = 1}^n {{m_j}} $(mj为与弹体接触的第j根钢筋的质量),列动量定理:

      $\frac{\rm d}{{{\rm d}{t_2}}}(m{v_{i{\rm{1}}}}) = - 2n{Q_{\rm{p}}} - q\sum\limits_{j = {\rm{1}}}^n {{b_j}} - {F_{{{\textit{z}}i}}}$

      式中:bj为第j根钢筋与弹头部分重合的长度。

      对其积分得:

      $2n{Q_{\rm{p}}}{t_2} + q{t_2}\sum\limits_{j = 1}^n {{b_j}} + {F_{{{\textit{z}}i}}}{t_2} = m{v_i} - m{v_{i1}}$

      可得,弹体的速度:

      ${v_{i1}} = {v_i} - \frac{{2n{Q_{\rm{p}}}{t_2} + q{t_2}\displaystyle\sum\limits_{j = 1}^n {{b_j}} + {F_{{{\textit{z}}i}}}{t_2}}}{m}$

      加速度:

      ${\tilde a_i} = - \frac{{2n{Q_{\rm{p}}} + q\displaystyle\sum\limits_{j = 1}^n {{b_j}} + {F_{{{\textit{z}}i}}}}}{m}$

      由式(3)、(9)、(19),每根钢筋截面A两侧的相对位移为:

      ${S_{1i}} = {v_i}{t_2} - \left(\frac{{n\rho {\lambda _i}}}{m} + 1\right)\left(\frac{{{Q_{\rm{p}}}t_2^2 - q{\lambda _i}t_2^2}}{{\rho {\lambda _i}}}\right) - \frac{{qt_2^2\displaystyle\sum\limits_{j = 1}^n {{b_j}} + {F_{{{\textit{z}}i}}}t_2^2 + 2nq{\lambda _i}t_2^2}}{{2m}}$

      弹体头部接触钢筋,且钢筋剪切失效时的弹体侵彻位移:

      ${S_i} = {v_i}{t_2} - \frac{{2n{Q_{\rm{p}}}t_2^2 + qt_2^2\displaystyle\sum\limits_{j = 1}^n {{b_j}} + {F_{{{\textit{z}}i}}}t_2^2}}{{2m}}$

    • Forrestal等[4]利用空腔膨胀理论得出弹体在隧道区受到混凝土的阻力为:

      ${F_{{{\textit{z}i}}}} = {\text{π}} \frac{{d_0^2}}{4}\left( {S{{f'_{\rm{c}}}} + N{\rho _0}{v_{i}(t)}^2} \right)$

      式中:$S = 82.6{{f'}_{\rm{\!\!\!\!c}}^{ - 0.544}}$${f'_{\rm{c}}}$为混凝土无侧限抗压强度;$N = 1$为平头弹弹头形状系数,ρ0为混凝土密度,vi(t)为弹体任意时刻t的速度。

      任意时刻的加速度、速度、位移表达式为:

      ${{a_i}{(t + \Delta t)}} = - \frac{{{F_{{{\textit{zi}}}}}}}{{{m_0}}}$

      ${{v_i}{(t + \Delta t)}} = {v_{i}(t)} - \frac{{{F_{{{\textit{zi}}}}}}}{{{m_0}}}\Delta t$

      ${{s_i}{(t + \Delta t)}} = {v_{i}(t)}\Delta t - \frac{{{F_{{{\textit{zi}}}}}}}{{2{m_0}}}\Delta {t^2}$

    • 图2为弹体侵彻多层钢筋的示意图,其中,n为钢筋层数,h1为首层钢筋的埋置深度,hn为钢筋网间距。在隧道阶段,依据弹体是否与钢筋直接撞击,分为两种工况:(1)弹体撞击到钢筋,侵彻阻力包括钢筋与混凝土两部分的作用;(2)弹体未撞击到钢筋或钢筋失效,此时弹体仅受到混凝土提供的侵彻阻力。

      图  2  弹体侵彻半无限厚靶过程示意图

      Figure 2.  Schematic diagram of a projectile penetrating a semi-infinite thick target

      图3为弹体侵彻半无限厚靶的算法流程图。其中,首先输入的参数包括(1)弹体参数:弹体质量、直径、长度、侵彻速度;(2)钢筋参数:钢筋密度、屈服强度、直径、首层钢筋埋置深度、钢筋网间距;(3)混凝土参数:混凝土的无侧限抗压强度、密度等。根据首层钢筋埋置深度与开坑深度的大小关系,判定弹体在进入隧道阶段时,是否撞击到钢筋,之后不断循环钢筋与混凝土阶段,直到弹体速度为零,其中在每一个阶段,取适当时间步长△t计算每个阶段的参量,不断循环,直到达到其边界条件,跳出此阶段进入下一个阶段。

      图  3  弹体侵彻半无限厚靶算法示意图

      Figure 3.  Schematic diagram of a projectile penetrating a semi-infinite thick target

    • Young[1-2]通过开展大量实验研究,得出了侵彻土、岩石、钢筋混凝土的统一公式:

      $H \!\!=\!\! \left\{ \begin{aligned} & 0.000\;8\;KRN{(M/A)^{0.7}}{\rm{ln}}(1\!\! +\!\! 2.15v_{\rm{s}}^2 \!\times\! {10^{ - 4}}),\\&{v_{\rm{s}}} {\text{<}} 61\;{\rm{m/s}} \\ & 0.000\;018\;KRN{(M/A)^{0.7}}({v_{\rm{s}}} - 30.5),\\&{v_{\rm{s}}} {\text{≥}} 61\;{\rm{m/s}} \\ \end{aligned} \right.$

      式中:M为弹体质量;A为弹体的截面面积;vs为弹体侵彻靶体的初始撞击速度;K为缩尺效应系数,当M<182 kg时,K=0.46M0.15;当M≥182 kg时,K=1。N为弹头性能系数,对于平头弹N=0.56;R为阻力系数,对混凝土R =0.085$ {K_{\rm{c}}}(11 - P){({t_{\rm{c}}}{h_{\rm{c}}})^{ - 0.06}}{({{35} / {{\sigma _{\rm{c}}}}})^{0.3}}$,其中σc为实验时混凝土的无侧限抗压强度;${K_{\rm{c}}} = {\left( {{F / {{W_1}}}} \right)^{0.3}}$W1为靶体宽度与弹体直径的比值,F=20,如果W1FKc=1;P为混凝土中按体积计算的含钢百分率;tc为混凝土的凝固时间;hc为混凝土目标的厚度。在没有足够的数据无法计算R值时,采用R=0.9。

      由于缺乏大口径平头弹侵彻钢筋混凝土厚靶的公开实验数据,本文中将不同弹体直径和质量(80、90、100 kg作为工况1~3)、不同撞击速度(600、410、780 m/s作为工况3~5)等5种工况的计算结果与Young公式的计算结果对比来验证模型的有效性。选取的靶体的参数:靶体尺寸为2 m×2 m×2 m,无侧限抗压强度fc'=35 MPa,钢筋密度为7 850 kg/m3、钢筋屈服强度为400 MPa,钢筋直径d1=8 mm,网眼尺寸60 mm,首层钢筋埋置深度为200 mm,体积配筋率2%,钢筋层间距为160 mm;弹体参数见表1

      工况弹体参数弹体与钢筋的相对位置侵彻深度/mm误差/%
      直径/mm长度/mm质量/kg撞击速度/(m·s−1)
      18053213600钢筋交汇处 6490.6
      钢筋网格一边中点 6490.6
      钢筋网眼中心 6490.6
      Young公式 653
      29059918600钢筋交汇处 7122.47
      钢筋网格一边中点 7112.6
      钢筋网眼中心 7112.6
      Young公式 730
      310066625600钢筋交汇处 8013.84
      钢筋网格一边中点 8003.96
      钢筋网眼中心 8003.96
      Young公式 833
      410066625410钢筋交汇处 5137.57
      钢筋网格一边中点 5137.57
      钢筋网眼中心 5127.75
      Young公式 555
      510066625780钢筋交汇处1 0464.56
      钢筋网格一边中点1 0454.65
      钢筋网眼中心1 0454.65
      Young公式1 096

      表 1  计算侵彻深度

      Table 1.  Calculating penetration depth

      此外,由于弹体撞击到钢筋的位置具有随机性,本文中,如图4所示,以3种典型位置的计算结果与Young公式进行比较,比较结果见表1表1中计算结果与Young公式吻合良好,进一步验证了模型的可靠性。

      图  4  不同着靶点位置

      Figure 4.  Different target positions

    • 2.1节的分析表明:当钢筋直径较小时,不同着靶位置对于侵彻深度的影响较小。本节中,将弹靶参数改为:弹体长度为666 mm,质量m0=25 kg,弹体撞击速度为304 m/s;靶体尺寸为2 m×2 m×2 m,无侧限抗压强度$f'_{\rm c} $=35 MPa,钢筋直径d1=30 mm,网眼尺寸60 mm,首层钢筋埋置深度100 mm,体积配筋率2%,钢筋层间距为50 mm。并分别对弹体直径为100 mm和140 mm进行分析,着重研究不同着靶点对于侵彻深度的影响,计算结果见图56。从图5中可以看出,弹体撞击到钢筋交汇处侵彻深度最大;从图6中可以看出,弹体撞击到钢筋网眼中心侵彻深度最大。

      图  5  侵彻深度随时间变化图(弹体直径为100 mm)

      Figure 5.  Penetration depth varing time at different target positions for the projectile with the diameter of 100 mm

      图  6  侵彻深度随时间变化图(弹体直径为140 mm)

      Figure 6.  Penetration depth varing time at different target positions for the projectile with the diameter of 140 mm

    • 根据防护角度的最不利位置,平头弹侵彻钢筋混凝土厚靶的工程计算公式推导如下:

      由2.2知,弹体直径与网眼尺寸的比值不同,最大侵彻深度的着靶点不同。由于无法预估撞击到的钢筋层数,且钢筋对于弹体的阻滞作用较小,因此,通过计算弹体侵彻素混凝土靶体的侵彻深度,可根据钢筋层间距估算弹体撞击的钢筋层数和本文提出的解析模型,得到最终的侵彻深度。

      以在侵彻阶段首次撞击到钢筋时为计算起始点,此时速度为v1

      (1)根据Forrestal等[4]计算侵彻素混凝土靶体的侵彻深度s

      $s= \frac{{2m{}_0}}{{{\text{π}} d_0^2\rho _0^2}}{\rm{ln}}\left(1 + \frac{{{\rho _0}}}{{S{{f'_{\rm{c}}}}}}v^2\right) + 0.707{d_0}$

      式中:$v^{\rm{2}} = {{\left( {2{m_0}v_0^2 - {\text{π}} d_0^3S{{f'_{\rm{c}}}}} \right)} \big/ {\left( {2{m_0} + {\text{π}} d_0^3{\rho _0}} \right)}}$

      (2)估算弹体撞击的整钢筋层间距的层数x$ x = \left[\dfrac{s}{{{h_i}}}\right]$,其中[]为取整函数,hi为钢筋层间距;

      (3)得到弹体撞击到x+1层钢筋时的速度vx+1

      ①由混凝土阻力引起速度的减小,从弹体首次撞击到钢筋到第x层钢筋,依据动能定理得:

      $ - \displaystyle\int\limits_0^{x{h_i}} {{F_{\textit z}}{\rm d}x} = \displaystyle\int\limits_{{v_1}}^{{v_{x+1}}} {\dfrac{1}{2}{m_0}{\rm d}{v^2}},\quad\quad\quad v_{x+1}^2 = \frac{{{{\rm e}^{ - \frac{2}{{{m_0}}}x{h_i}B}}(A + Bv_1^2) - A}}{B}$

      式中:${F_{\textit{z}}} = {\text{π}} \dfrac{{d_0^2}}{4}\left( {S{{f'_{\rm{c}}}} + N{\rho _0}{v^2}} \right) = A + B{v^2}$$A = {\text{π}} \dfrac{{d_0^2}}{4}S{f'_{\rm{c}}}$$B = {\text{π}} \dfrac{{d_0^2}}{4}N{\rho _0}$

      ②由钢筋阻力引起的速度减小量:

      $\Delta {v^2} = \dfrac{{2\left(2n{Q_{\rm{p}}} + q\displaystyle\sum\limits_{j = 1}^n {{b_j}} \right)}}{{{m_0}}} \times 0.3{d_1} x $,则:

      x+1层钢筋时弹体速度为:

      $v_{x + 1}^2 = \frac{{{{\rm e}^{ - \frac{2}{{{m_0}}}x{h_i}B}}(A + Bv_1^2) - A}}{B} - \frac{{2\left(2n{Q_{\rm{p}}} + q\displaystyle\sum\limits_{j = 1}^n {{b_j}} \right)}}{{{m_0}}} \times 0.3{d_1} x$

      (4)隧道阶段的侵彻深度${S'_1}$为:

      $\left\{ \begin{split}&{S'_1} = x{h_{\rm{i}}} - \frac{{v_{x + 1}^2}}{{2{{\tilde a}_{x + 1}}}},\;{\text{若}}\;v_{x + 1}^2 {\text{>}} 0\;\;{\text{且}}\;\;{S_{x + 1}} = - \frac{{v_{x + 1}^2}}{{2{{\tilde a}_{x + 1}}}} {\text{≤}}0.3{d_1}\\ & {S'_1} = x{h_{\rm{i}}} + 0.3{d_1} - \frac{{{m_0}}}{{2B}}l{\rm{n}}\frac{A}{{A + B(2{{\tilde a}_{x + 1}} \times 0.3{d_1} + v_{x + 1}^2)}},\;\;\;\;{\text{若}}\;\;v_{x + 1}^2 {\text{>}} 0\;\;{\text{且}}\;{S_{x + 1}} = - \frac{{v_{x + 1}^2}}{{2{{\tilde a}_{x + 1}}}} {\text{>}} 0.3{d_1} \end{split}\right. $

      $v_{x + 1}^2 {\text{<}} 0$时,此时弹体撞击到第x层,因此将式(2)中的xx=x−1带入式(3)、(4)中计算。

      当弹体直径小于钢筋网眼尺寸时,弹体撞击到钢筋网眼中心时只受到混凝土提供的阻力,相比其他情况,可知此时受到的阻力越小,侵彻深度越大,越为不利。

      综上所述,总侵彻深度L的公式为:

      $\left\{ \begin{split}&L = - \frac{{{m_0}}}{{2B}}{\rm{ln}}\frac{A}{{A + Bv_0^2}} + 0.707{d_0},\;\;{\text{若}}\;\;k {\text{<}}1\\ & L = x{h_i} - \frac{{v_{x + 1}^2}}{{2{{\tilde a}_{x + 1}}}} + {h_1}=xh_i+S_{x+1}+h_1\;,\;{\text{若}}\;\;k{\text{>}}1\;\;{\text{且}}\;\; {S_{x + 1}} = - \frac{{v_{x + 1}^2}}{{2{{\tilde a}_{x + 1}}}} {\text{≤}} 0.3{d_1}\\ & L = x{h_i} + 0.3{d_1} - \frac{{{m_0}}}{{2B}}\ln\frac{A}{{A + BS_{x+1}(2{{\tilde a}_{x + 1}} \times 0.3{d_1} /S_{x+1}-2{{\tilde a}_{x + 1}})}} + {h_1}\;\;{\text{若}}\;k {\text{>}} 1,\;\;{\text{且}}\;\;{S_{x + 1}} = - \frac{{v_{x + 1}^2}}{{2{{\tilde a}_{x + 1}}}} {\text{>}} 0.3{d_1}\end{split}\right. $

      式中:k为弹体直径与钢筋网眼尺寸的比值,h1为首层钢筋埋置深度。

      $v_{x + 1}^2 = \frac{{{{\rm e}^{ - \frac{2}{{{m_0}}}x{h_i}B}}(A + Bv_1^2) - A}}{B} - \frac{{2\left(2n{Q_{\rm{p}}} + q\displaystyle\sum\limits_{j = {\rm{1}}}^n {{b_j}} \right)}}{{{m_0}}} \times 0.3{d_1} \times x;$

      $ v_1^2 = \frac{{{{\rm e}^{ - \frac{2}{{{m_0}}}({h_1} - 0.707{d_0})B}}(A + Bv_0^2) - A}}{B}, x = \left[ {\frac{s}{{{h_i}}}} \right]{\text{,}}{\text{若}}\; v_{x + 1}^2 {\text{<}} 0{\text{,}}{\text{则}}\;x = \left[ {\frac{s}{{{h_i}}}} \right] - 1; $

      ${\tilde a_{x + 1}} = - \frac{{2n{Q_{\rm{p}}} + q\displaystyle\sum\limits_{j = 1}^n {{b_j}} + {F_{{{\rm{z}}_{x + 1}}}}}}{{{m_0}}};$

      $\displaystyle\sum\limits_{j = {\rm{1}}}^n {{b_j}} = \min \left\{ {\displaystyle\sum {{b_{{n_1}}}} ,\displaystyle\sum {{b_{{n_2}}}} ,\displaystyle\sum {{b_{{n_3}}}} } \right\}$n取与$\displaystyle\sum\limits_{j = 1}^n {{b_j}} $对应的ni值,对应如下:

      钢筋交汇处:${n_1} = 4\left[\dfrac{k}{2}\right] + 2$$\displaystyle\sum {{b_{{n_1}}}} = \left(8\displaystyle\sum\limits_{i = 1}^{[\frac{k}{2}]} {\sqrt {\dfrac{1}{4} - \frac{{{i^2}}}{{{k^2}}}} } + 2\right){d_0}$,

      钢筋网眼中心:${n_2} = 4\left[\dfrac{k}{2} - \dfrac{1}{2}\right] + 4$$\displaystyle\sum {{b_{{n_2}}}} = 8\displaystyle\sum\limits_{i = {\rm{1}}}^{\frac{{{n_2}}}{4}} {\sqrt {\dfrac{1}{4} - {{\left(i - \dfrac{1}{2}\right)}^2}\dfrac{1}{{{k^2}}}} } {d_0}$

      一根钢筋中点:

      ${n_3} = 2\left[\dfrac{k}{2}\right] + 2\left[\dfrac{k}{2} - \dfrac{1}{2}\right] + 3$$\displaystyle\sum {{b_{{n_3}}}} = \left(4\displaystyle\sum\limits_{i = 1}^{[\frac{k}{2}]} {\sqrt {\dfrac{1}{4} - \dfrac{{{i^2}}}{{{k^2}}}} } + 4\displaystyle\sum\limits_{i = 1}^{[\frac{k}{2} - \frac{1}{2}] + 1} {\sqrt {\dfrac{1}{4} - {{\left(i - \dfrac{1}{2}\right)}^2}\dfrac{1}{{{k^2}}}}} + 1\right){d_0}$

    • 本文基于弹体侵彻素混凝土的经典理论,将钢筋的失效模式考虑为弯曲剪切破坏后,建立平头弹侵彻钢筋混凝土的阻力模型。将计算所得的结果与经典Young公式进行对比,结果吻合良好,表明了模型的可靠性,并能够较为合理地预测侵彻深度及其它侵彻过程量并给出弹体各运动参量随时间变化关系,通过分析不同着靶点位置对于侵彻深度的影响,发现当弹体直径小于钢筋网眼尺寸时,弹体撞击到钢筋网眼中心时,弹体受到的阻力较小,侵彻深度越大;当弹体直径大于钢筋网眼尺寸时,最不利着靶点位置视弹体直径与网眼尺寸的比值而定。最后,根据防护角度的最不利位置提出了平头弹侵彻钢筋混凝土厚靶的工程计算公式。

参考文献 (13)

目录

    /

    返回文章
    返回