大尺度锤头冲击下固支方板的动响应分析

甘霖 伍星星 王海坤 刘建湖 赵延杰

甘霖, 伍星星, 王海坤, 刘建湖, 赵延杰. 大尺度锤头冲击下固支方板的动响应分析[J]. 高压物理学报, 2021, 35(6): 064103. doi: 10.11858/gywlxb.20210733
引用本文: 甘霖, 伍星星, 王海坤, 刘建湖, 赵延杰. 大尺度锤头冲击下固支方板的动响应分析[J]. 高压物理学报, 2021, 35(6): 064103. doi: 10.11858/gywlxb.20210733
 GAN  Lin, WU Xingxing, WANG Haikun, LIU Jianhu, ZHAO Yanjie. Dynamical Response Behavior of Clamped Square Plates under the Impact of Large-Scale Hammer[J]. Chinese Journal of High Pressure Physics, 2021, 35(6): 064103. doi: 10.11858/gywlxb.20210733
Citation:  GAN  Lin, WU Xingxing, WANG Haikun, LIU Jianhu, ZHAO Yanjie. Dynamical Response Behavior of Clamped Square Plates under the Impact of Large-Scale Hammer[J]. Chinese Journal of High Pressure Physics, 2021, 35(6): 064103. doi: 10.11858/gywlxb.20210733

大尺度锤头冲击下固支方板的动响应分析

doi: 10.11858/gywlxb.20210733
基金项目: 国防基础科研重点项目(B0820132045);国家重点安全基础研究项目(613279)
详细信息
    作者简介:

    甘 霖(2000-),男,本科,主要从事交通安全研究. E-mail:ganlin1230620@163.com

  • 中图分类号: O346

Dynamical Response Behavior of Clamped Square Plates under the Impact of Large-Scale Hammer

  • 摘要: 为探究大尺度锤头冲击下固支方板的动响应行为,结合落锤工装设计了大尺度冲击锤头,通过对固支方板开展不同冲击强度下的撞击试验,获取了固支方板在大尺度锤头撞击下的典型破坏模式。同时依据试验结果及冲击动力学理论,建立了大尺度锤头冲击下固支方板变形评估方法,并结合数值模拟建立了大尺度锤头冲击下方板边界撕裂判据。研究结果表明:大尺度锤头冲击下固支方板主要依靠边界塑性铰、面板内塑性铰耗能,固支方板的撕裂从边界位置开始。基于塑性铰耗能机制建立了方板塑性变形理论评估方法,理论计算结果与试验结果吻合较好。固支方板撕裂过程中,边界起始撕裂处应力三轴度基本在0.6左右,本试验中建立的临界失效判据Bh值可取1.6倍板厚。所建立的评估方法及失效判据进一步完善了落锤冲击下固支方板的动响应行为研究。

     

  • 随着空间科学研究的发展,为进一步开展空间碎片毁伤效应、防护结构以及超高速碰撞等研究工作,提出了10 km/s以上的超高速发射需求。国内外常用的发射技术(设备)包括:二级/三级轻气炮、磁驱动飞片、聚能射流、等离子体加速器、激光驱动飞片装置等[1-6]。其中二级轻气炮是目前应用最广泛、技术最为成熟的发射设备。然而,受到炮体材料性能的限制,二级轻气炮的稳定发射速度基本都低于8 km/s,如果超过该速度,发射管将发生严重的烧蚀破坏[7]。虽然有些发射技术能够实现8 km/s以上的超高速发射,但是这些发射技术要么难以控制弹丸的质量和形状,要么只能发射很薄的飞片状弹丸或质量很小的弹丸[8],因此将其应用到超高速碰撞实验研究中受到了一定的限制。在20世纪60年代,内爆驱动式超高速发射技术已经实现了2 g弹丸12.2 km/s的发射速度[9-10];进入21世纪后,加拿大麦吉尔大学[1, 11-13]又实现了15 g弹丸7.5 km/s左右的发射速度,这是一项发射质量和速度都超过二级轻气炮的超高速发射技术,能够满足空间碎片撞击的研究需求。

    内爆驱动式超高速发射技术(国外有Explosive hypervelocity launcher[9]、Explosively driven gun[10]、Implosion-driven launcher[11-13]等不同说法,国内也有炸药爆轰驱动高速激波管发射技术[14]、炸药驱动枪[15]等不同名称,本文以Implosion-driven launcher说法为准,简称内爆发射器)是一种利用炸药爆炸线性挤压压缩管内轻质气体,通过被压缩的轻质气体驱动弹丸的一种发射技术[11],作用过程与二级轻气炮中活塞压缩高压轻质气体类似,不同的是其活塞速度为炸药爆速。国内一些文献中对该技术进行了相关介绍,但尚未见相关实验和仿真研究报道[14-15]

    为了掌握内爆发射技术,获得10 km/s以上的发射能力,通过与加拿大麦吉尔大学的Higgins教授合作,以8 mm口径内爆发射装置为研究对象,对内爆发射器的发射能力和研究方法进行了初步探索,开展了内爆发射器内弹道分析和结构设计,并利用AUTODYN软件对设计结构进行了初步仿真分析,最终针对典型状态开展了初步验证实验。

    内爆发射器主要由炸药线性压缩段、高压段和发射管组成,其中炸药线性压缩段包括外壳、装药、压缩管、驱动气体等,发射器的结构如图1所示,作用过程如图2所示。左端炸药起爆后,压缩管向内塌缩,形成锥形结构,当炸药爆轰波传播时,锥形结构向右运动并压缩驱动气体。该过程与二级轻气炮中活塞压缩高压轻质气体类似,因此该锥形结构可看成虚拟活塞,当冲击波到达弹丸底部时,高压驱动气体开始推动弹丸加速运动。

    图  1  内爆式超高速发射器结构示意图[11]
    Figure  1.  Structural diagram of implosion-driven hypervelocity launcher[11]
    图  2  内爆式超高速发射器工作原理[11]
    Figure  2.  Working diagram of implosion launcher after ignition[11]

    设计了内径16 mm、长度500 mm的压缩管和内径8 mm、长度300 mm的发射管组成的内爆发射器,高压段采用锥段连接压缩管和发射管,锥段半锥角度为5°。采用AUTODYN软件对该内爆发射器的工作过程进行了仿真计算。

    采用有限元中的流固耦合方法进行仿真。将模型分为Euler模型和Lagrange模型两个部分:Euler模型包括填充的氦气、压缩管、装药以及周围空气介质;Lagrange模型包括压缩段外壳、高压段、发射管以及弹丸。外壳、高压段、压缩管、发射管均采用4340钢。本构模型采用Johnson Cook模型,状态方程为Linear。装药采用爆速为7 km/s左右的季戊四醇四硝酸酯(PETN),状态方程为JWL;驱动气体为氦气,采用理想气体状态方程;弹丸本构模型为Steinberg Guinan,状态方程为Shock。材料参数均选自AUTUDYN数据库。

    为了使弹丸能够约束初始状态的高压氦气,在弹丸与高压气体交界处增加一层Lagrange网格。该网格与高压段和弹丸均采用共节点的方式固定,并将该部分的材料失效应变设置为0.5,当达到失效应变时自动删除该部分的网格,使其不会对弹丸发射造成过大影响。

    针对弹丸材料和充气压力两个参数,共进行4个状态的仿真计算,相关参数以及各状态下获得的弹丸出炮口速度见表1

    表  1  不同弹丸材料和充气压力下仿真计算参数
    Table  1.  Simulation parameters of launchers with different projectile materials and filling pressure
    No.Material of projectileFilling pressure/MPaDiameter of launch tube/mmDiameter of pump tube/mmMuzzle velocity/(km·s–1
    1Aluminum alloy4816 8.62
    2Aluminum alloy5816 9.17
    3Aluminum alloy6816 9.25
    4Magnesium alloy581610.66
    下载: 导出CSV 
    | 显示表格
    2.3.1   内爆发射器工作过程

    图3给出了内爆式超高速发射装置发射铝弹丸的仿真结果,其中弹丸尺寸为$\varnothing$8 mm × 4 mm,质量为0.55 g。零时刻从装药的左端面开始点火,装药爆炸后开始挤压压缩管形成锥形活塞,并在氦气(初始压力5 MPa)中形成压力达到500 MPa左右的冲击波向右传播,外壳达到失效应变后被删除,如图3(b)所示,压力超过1 GPa时显示为红色;在55 μs左右,冲击波到达弹丸底部反射,形成峰值压力高达5 GPa的气体推动弹丸加速运动;70 μs左右炸药爆轰结束;在108 μs左右弹丸离开发射管,达到最高发射速度。

    图  3  8 mm口径内爆式发射装置工作过程
    Figure  3.  Launch process of 8 mm caliber implosion launcher
    2.3.2   不同初始充气压力下弹丸的发射速度

    图4给出了充气压力分别为4、5和6 MPa状态下铝弹丸的速度-时间历程。从速度-时间曲线上可以看到,弹丸加速过程有较为明显的二次冲击加载。初始冲击加载是由先导冲击波到达弹丸底部产生的;二次冲击加载是弹丸底部反射的冲击波到达左端后再次反射冲击波,该冲击波追赶上弹丸产生的。从表1中的数据可知,4、5和6 MPa充气压力下弹丸的出炮口速度分别为8.62、9.17和9.25 km/s。可以看到,炮口速度随充气压力的增加而增加,但增加幅度缩小。可能是由于当充气压力从5 MPa增加至6 MPa时,压缩管、高压段内壁的膨胀显著增加,削弱了驱动压力增加的增速效果。

    图  4  不同充气压力下铝弹丸的速度-时间历史
    Figure  4.  Velocity-time history of aluminum projectiles with different filling pressure
    2.3.3   不同弹丸材料的发射速度

    根据建立的仿真模型对尺寸为$\varnothing$8 mm × 4 mm的铝合金弹丸和镁合金弹丸进行发射过程的计算,初始填充5 MPa氦气,计算结果如图5所示。0.55 g铝合金弹丸的最终发射速度为9.17 km/s;0.37 g镁合金弹丸的发射速度为10.66 km/s,比铝合金弹丸速度高1.49 km/s。

    图  5  不同材料弹丸速度-时间历史
    Figure  5.  Velocity-time history of projectiles with different materials

    发射器总长度为914 mm,最大部位直径为81.4 mm,主要包括发射管、发射管护套、高压段、外壳、压缩管、弹丸和装药等。其中:发射管内径为8 mm,长度为304 mm;压缩管内径为16 mm,长度为495 mm。装药为柱壳装药,厚度约为5 mm,装药量约为300 g,采用两层以黑索金(RDX)为基底的片状挠性炸药包裹而成,挠性炸药的爆速为7.1 km/s。弹丸直径为8 mm,高度为4 mm。弹丸材料有两种,分别是铝合金和镁合金,其中铝合金弹丸质量约为0.55 g,镁合金弹丸质量约为0.37 g,安装完成后的发射器如图6所示。

    图  6  内爆式超高速发射器实物
    Figure  6.  Image of implosion-driven hypervelocity launcher

    实验装置布局如图7所示。实验装置主要包括发射器、测试系统、真空系统和防护系统等。发射器放置在厚壁防护筒内,防止爆炸产生的碎片向外飞散,损伤实验装置和设备。采用激光测速系统和靶网测速装置等测量弹丸的发射速度,并采用序列激光阴影成像仪拍摄模型的飞行姿态和形貌。实验中发射器压缩管内充5 MPa氦气,从端面起爆装药。

    图  7  实验装置布局示意图
    Figure  7.  Layout of experimental equipment

    共开展了7次实验,其中铝合金弹丸实验3次,镁合金弹丸实验4次,具体实验参数如表2所示。

    表  2  内爆发射器实验相关参数
    Table  2.  Parameters of the implosion-driven launchers in tests
    No.Material of projectileMass of projectile/gFilling pressure/MPaDiameter of launch tube/mmDiameter of pump tube/mmMaximum velocity/(km·s–1
    ILT04Aluminum alloy0.555816 7.26
    ILT07 6.96
    ILT14 7.95
    ILT08Magnesium alloy0.375816 9.73
    ILT0910.28
    ILT11 9.36
    ILT12 9.77
    下载: 导出CSV 
    | 显示表格
    3.3.1   回收发射器

    实验后的发射器压缩管和外壳完全碎裂,高压段外壳明显膨胀变形并出现裂纹,回收的典型发射器装置如图8所示。

    图  8  ILT08实验发射器装置回收
    Figure  8.  Recycle launcher of test ILT08
    3.3.2   铝合金弹丸速度测试结果

    3次铝合金弹丸实验结果如图9所示,其中横坐标为测点位置离炮口的距离,纵坐标为测得的速度,各次实验测得的最大速度见表2。其中,实验中测得的弹丸最高出炮口速度为7.95 km/s,比仿真预测速度9.17 km/s低1.22 km/s,相对偏差为15.3%。实验速度偏低可能是许多因素共同影响的结果,例如装置加工误差、装药松紧程度、发射管与弹丸有摩擦阻力等,而仿真状态则相对理想,没有考虑这些因素的影响。另外,实验中弹丸的速度也有1 km/s左右的偏差,挠性炸药装填、装药爆轰的圆周同步性、发射管与弹丸的配合程度等问题都会对实际发射速度造成一定影响,导致重复实验出现一定的偏差。

    图  9  0.55 g铝合金弹丸速度测试结果
    Figure  9.  Experimental muzzle velocities of 0.55 g aluminum projectile
    3.3.3   镁合金弹丸速度测试结果

    4次镁合金有效实验获得的速度结果如图10所示,实验测得的最大速度见表2。4次实验获得的最高发射速度与仿真预测速度10.66 km/s非常接近,相对偏差为3.7%。镁合金弹丸在实验中的速度偏差为1 km/s左右,产生偏差的原因可能与铝合金弹丸相同。

    图  10  0.37 g镁合金弹丸速度测试结果
    Figure  10.  Experimental muzzle velocities of 0.37 g magnesium projectile
    3.3.4   弹丸成像结果

    通过超高速八序列激光阴影成像仪拍摄弹丸在实验过程中的飞行情况,图11为3次实验拍摄的照片,其中ILT11、ILT12为镁合金弹丸,ILT14为铝合金弹丸。在图11中可以看到,实验过程中的弹丸仍保持为一个整体弹丸,但存在一定程度变形,其中实验ILT12的弹丸变形非常明显。发生变形的原因可能是材料强度不够,也可能是发射器结构存在一定的缺陷,有待进一步深入研究。

    图  11  弹丸序列激光阴影成像结果
    Figure  11.  Projectile photos shot by sequence laser shadowgraph imager

    通过与加拿大麦吉尔大学合作,初步设计了8 mm口径内爆式超高速发射器结构,采用AUTODYN软件对内爆式超高速发射器内弹道进行数值模拟,所获得的发射器作用过程中弹底压力最高可达5 GPa,并计算了铝合金弹丸和镁合金弹丸的发射速度。最终通过开展验证实验,证明了设计的内爆发射技术能够将尺寸为$\varnothing$8 mm × 4 mm、质量为0.37 g的柱状镁合金弹丸发射到10.28 km/s,将尺寸为$\varnothing$8 mm × 4 mm、质量为0.55 g的柱状铝合金弹丸发射到7.95 km/s。实验获得的发射速度与数值仿真结果吻合较好。实验结果也证明了内爆式超高速发射技术的可行性,在空间碎片的撞击与防护方面具有较强的应用前景。利用超高速八序列激光阴影成像仪拍摄弹丸的飞行情况,获得了弹丸发射后的形状,发现发射后的弹丸是一个整体弹丸,但与原始柱形结构相比,存在不同程度的变形。未来将深入探索提高内爆发射器速度的方法,分析弹丸发生变形的原因,并在保持弹丸初始形状的发射器结构设计等方面开展进一步的研究工作。

    感谢中国空气动力研究与发展中心的罗庆、龙耀、宋强、周毅、覃金贵、李俊玲、姜林、邹胜宇、李文光、廖强、刘晓龙、廖富强、丁建文等同事在实验过程中提供的无私帮助。

  • 图  落锤试验装置示意图

    Figure  1.  Schematic diagram of the drop hammer test device

    图  冲击锤头模型示意图(单位:mm)

    Figure  2.  Dimension of impact hammer head model (Unit: mm)

    图  方板的典型破坏模式

    Figure  3.  Typical failure modes of square plates

    图  矩形板的塑性铰线和横向速度场

    Figure  4.  Plastic hinge and transverse velocity distribution of square plates

    图  Q345B钢在不同应变率下的应力-应变曲线

    Figure  5.  Stress-strain curves of Q345B steel under different strain rates

    图  锤头撞击方板计算模型及提取单元示意图

    Figure  6.  Simulation model of the square plate under hammer impacting and the selected element

    图  Q345B钢材料试验结果及提取单元所处应力

    Figure  7.  Test results of Q345B steel material and the stress of the selected element

    图  典型固支方板计算模型示意图

    Figure  8.  Typical numerical simulation models of the clamped square plates

    图  计算结果

    Figure  9.  Calculation results

    表  1  试验工况及试验结果

    Table  1.   Experiment cases and test results

    Case
    h/m
    v/(m·s−1)
    M/kg
    E/kJ
    Result
    D/mmDamage
    10.42.802 024 7.93439.3Plastic deformation
    20.53.132 024 9.91745.9Boundary tearing
    30.63.422 02411.90148.6Boundary tearing
    40.73.702 02413.88455.8Boundary tearing
    下载: 导出CSV
  • [1] ZHU L, FAULKNER D, ATKINS A G. The impact of rectangular plates made from strain-rate sensitive materials [J]. International Journal of Impact Engineering, 1994, 15(3): 245–255. doi: 10.1016/S0734-743X(05)80016-8
    [2] ZHU L, FAULKNER D. Dynamic inelastic behaviour of plates in minor ship collisions [J]. International Journal of Impact Engineering, 1994, 15(2): 165–178. doi: 10.1016/S0734-743X(05)80029-6
    [3] SHEN W Q, WONG P S, LIM H C, et al. An experimental investigation on the failure of rectangular plate under wedge impact [J]. International Journal of Impact Engineering, 2003, 28(3): 315–330. doi: 10.1016/S0734-743X(02)00102-1
    [4] SUN B, HU Z Q, WANG G. An analytical method for predicting the ship side structure response in raked bow collisions [J]. Marine Structures, 2015, 41: 288–311. doi: 10.1016/j.marstruc.2015.02.007
    [5] 陈继恩. 基于应力三轴度的材料失效研究[D]. 武汉: 华中科技大学, 2012.

    CHEN J E. Research of material failure basic on stress triaxiality [D]. Wuhan: Huahong University of Science and Technology, 2012.
    [6] GONG Y F, GAO Y X, XIE D, et al. Deflection and fracture of a clamped plate under lateral indentation by a sphere [J]. Ocean Engineering, 2015, 103: 21–30. doi: 10.1016/j.oceaneng.2015.04.075
    [7] 黄东, 龚榆峰, 刘敬喜. 壳板厚度对船舶加筋板结构耐撞性能的影响 [J]. 中国舰船研究, 2016, 11(2): 58–65. doi: 10.3969/j.issn.1673-3185.2016.02.009

    HUANG D, GONG Y F, LIU J X. The influence of variation of thickness for crashworthiness ability of stiffener plate structure [J]. Chinese Journal of Ship Research, 2016, 11(2): 58–65. doi: 10.3969/j.issn.1673-3185.2016.02.009
    [8] 诺曼•琼斯. 结构冲击[M]. 成都: 四川教育出版社, 1994.
    [9] 伍星星, 汪俊, 刘建湖, 等. 底部近距离爆炸下舱段模型毁伤试验研究 [J]. 爆炸与冲击, 2020, 40(11): 111406. doi: 10.11883/bzycj-2020-0067

    WU X X, WANG J, LIU J H, et al. Damaging characteristics of a cabin model under close-in underwater explosion from bottom attacting [J]. Explosion and Shock Waves, 2020, 40(11): 111406. doi: 10.11883/bzycj-2020-0067
    [10] 伍星星, 刘建湖, 孟利平, 等. 金属材料试件在压缩、扭转、拉伸断裂过程中的应力状态变化及表征 [J]. 高压物理学报, 2020, 34(5): 054204. doi: 10.11858/gywlxb.20200517

    WU X X, LIU J H, MENG L P, et al. Variation of stress distribution in metal fracture process under compressive, torsional, and tensile loading [J]. Chinese Journal of High Pressure Physics, 2020, 34(5): 054204. doi: 10.11858/gywlxb.20200517
    [11] 朱锡, 冯刚, 张振华. 爆炸载荷作用下固支方板的应变场及破坏分析 [J]. 船舶力学, 2005, 9(2): 83–89. doi: 10.3969/j.issn.1007-7294.2005.02.013

    ZHU X, FENG G, ZHANG Z H. Strain field and damage analisis of clamped square plates subjected to explosive loading [J]. Journal of Ship Mechanics, 2005, 9(2): 83–89. doi: 10.3969/j.issn.1007-7294.2005.02.013
  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  801
  • HTML全文浏览量:  379
  • PDF下载量:  43
出版历程
  • 收稿日期:  2021-03-08
  • 修回日期:  2021-03-26

目录

/

返回文章
返回