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

泡沫金属夹芯梁在重复冲击下的动态响应

朱凌 郭开岭 余同希 李应刚

胡宏伟, 宋浦, 王建灵, 郭炜, 徐洪涛, 金鹏刚, 任松涛. 炸药水中爆炸冲击因子的新型计算方法[J]. 爆炸与冲击, 2014, 34(1): 11-16.
引用本文: 朱凌, 郭开岭, 余同希, 李应刚. 泡沫金属夹芯梁在重复冲击下的动态响应[J]. 爆炸与冲击, 2021, 41(7): 073101. doi: 10.11883/bzycj-2020-0198
Hu Hong-wei, Song Pu, Wang Jian-ling, Guo Wei, Xu Hong-tao, Jin Peng-gang, Ren Song-tao. A new calculation method for shock factor of underwater explosion[J]. Explosion And Shock Waves, 2014, 34(1): 11-16.
Citation: ZHU Ling, GUO Kailing, YU Tongxi, LI Yinggang. Dynamic responses of metal foam sandwich beams to repeated impacts[J]. Explosion And Shock Waves, 2021, 41(7): 073101. doi: 10.11883/bzycj-2020-0198

泡沫金属夹芯梁在重复冲击下的动态响应

doi: 10.11883/bzycj-2020-0198
基金项目: 国家自然科学基金(11972269)
详细信息
    作者简介:

    朱 凌(1962- ),男,博士,教授,ZL79111@hotmail.com

    通讯作者:

    郭开岭(1989- ),男,博士,讲师,guokailing2008@163.com

  • 中图分类号: O347.3

Dynamic responses of metal foam sandwich beams to repeated impacts

  • 摘要: 为了研究重复冲击载荷作用下泡沫金属夹芯梁的动态响应,采用Abaqus数值仿真软件,基于可压碎泡沫模型(crushable foam),建立了泡沫金属夹芯梁遭受楔形质量块冲击的有限元模型。通过将仿真获得的夹芯梁上下面板最终挠度与重复冲击实验结果进行对比,验证仿真方法的准确性。在此基础之上,分析了泡沫金属夹芯梁在楔形质量块重复冲击作用下的变形模式、加卸载过程以及能量耗散特性。结果表明,在重复冲击载荷作用下,夹芯梁的变形不断累积,上面板主要出现局部凹陷和整体弯曲,而芯层则是局部压缩,下面板表现为整体弯曲。在重复加卸载过程中,加卸载刚度随着冲击次数的增加而增大。随着冲击次数的增加,上面板和芯层的能量吸收增量不断减小,而下面板的能量吸收增量不断增加,且最终均趋于稳定。泡沫金属夹芯梁的塑性变形能增量不断减小,而回弹系数随着冲击次数逐渐增加,最后趋于稳定值。
  • 冲击波作为爆破战斗部的主要毁伤方式之一,是评估武器毁伤威力的一项重要指标[1]。毁伤威力评估时,需要对战斗部炸点周围的冲击波场进行分布式测量,同时要求多测点实现同步触发。实战环境下战斗部落点及起爆时刻具有一定的随机性,且战斗部爆炸前具有一定的速度,战斗部的运动速度会改变爆炸冲击波的压力场分布。传统的触发方法如断线触发[2]、光触发[3]、无线触发[1]等均难以实现战斗部实战环境下冲击波超压的可靠触发,因此对战斗部动爆压力场的特性分析主要是通过仿真计算,并结合少量试验数据结果和爆炸相似律获得经验公式[4-6],缺少实战环境下的试验研究。

    本文中,提出了一种基于地震波可靠触发的战斗部空中爆炸冲击波超压测试方法,对着靶速度为0、535和980 m/s的战斗部空中爆炸冲击波进行了测试分析。结果表明,基于地震波信号触发测试方法能可靠获取战斗部动爆冲击波超压峰值。试验成果可为实战复杂环境下基于实测数据研究动爆冲击波特性提供依据。

    为了验证基于地震波实现冲击波超压测试触发的可行性,设计了基于加速度信号触发的冲击波超压测试系统并进行了试验。测试系统主要包括传感器、信号调理电路、电源管理、无线通信和信号采集存储等五部分,测试系统组成如图1所示。进行冲击波超压测试时,信号调理电路对传感器获得的冲击波超压信号进行滤波、降噪,在进入FPGA (field programmable gate array)之前进行信号有无和是否达到触发阈值的判断,A/D控制模块将模拟信号转为数字信号后存储在外部同步动态存储器(synchronous dynamic random-access memory, SDRAM)中,最后通过USB (universal serial bus)或无线通信模块将数据上传到上位机上,在上位机上完成数据最终的显示、分析和处理。

    图  1  测试系统
    Figure  1.  Test system

    基于地震波信号的冲击波超压测试触发方法借助配置在各冲击波超压测试节点上的加速度计与信号调理电路,利用地震波传播速度比冲击波传播速度快的特点,在冲击波到达各测试节点之前,提前感应到的加速度信号,经专用调理电路处理后,触发该测试节点的冲击波超压信号存储。其触发原理如图2所示,冲击波超压信号采集缓存区划分为循环采集和时序采集两部分,通过加速度信号启动冲击波超压信号的循环采集完成第1步触发,通过预设超压阈值的比较进行第2步触发,满足触发条件时,立即固化循环采集区,并开始冲击波超压的时序采集。

    图  2  地震波触发冲击波超压测试原理
    Figure  2.  The principle for measuring shock wave overpressure triggered by seismic wave

    爆炸时,形成以炸点为中心向四周传播的地震波,地震波传播速度最快的是纵波,其使地面发生上下振动,在地壳中的传播速度:u0=5.5~7 km/s。

    (1)根据金尼-格雷厄姆公式,计算冲击波超压峰值(pp,MPa):

    pppair=808[1+(fdR4.5)2]1+(fdR0.048)21+(fdR0.032)21+(fdR1.35)2,fd=3pairp0T0Tair
    (1)

    式中:R为比例距离,是观测点到距爆炸中心的距离r(m)与炸药TNT当量W(kg)的立方根之比,即R=r/3W,本文中0.053m/kg1/3R500m/kg1/3pair为试验现场大气压;p0为标准大气压,p0=101.325kPaTair为试验现场大气温度;T0为标准大气温度,T0=288.16K

    (2)根据Rankine-Hugoiot方程,爆炸点空气冲击波传播速度u与冲击波超压峰值pp之间的关系为:

    u=(ppp0×γ+12γ)+1c0
    (2)

    式中:p0为大气压,γ为空气的比热比,c0为波阵面前的空气声速。

    (3)冲击波和地震波传播到观测点的平均时间差Δt为:

    Δt=r(1u1u0)
    (3)

    联立式(1)~(3)可得不同观测距离处冲击波传播平均速度及其与地震波传播平均时间差的关系,如图3所示。由图3(a)可以看出,冲击波超压值随爆距的增大而减小,当距离r5m时,冲击波的传播平均速度u995m/s,该值远小于地震波的传播速度;由图3(b)可以看出,当距离r5m时,冲击波与地震波传播到观测点的平均时间差Δt4.1ms,而地震波信号触发该测点的冲击波信号只需要几十微秒,在该时间差内能够完成触发该测试节点的冲击波信号的存储。因此,当等效TNT装药量不大于100 kg、爆心距不小于5 m时,可以通过安装在测试终端的加速度计采集地震波信号,作为冲击波超压测试的可靠触发信号。

    图  3  冲击波超压峰值和传播时间差随距离的变化
    Figure  3.  Variation of shock wave overpressure peak and propagation time difference with distance

    图4是弹丸爆炸时在距爆心5 m处获取的加速度和冲击波超压信号,可以看出加速度信号的触发阈值先于冲击波信号到达测试节点,提前时长为12.42 ms,这进一步验证了基于地震波信号的冲击波测试触发方法原理的可行性。

    图  4  试验测试的加速度曲线和超压峰值曲线
    Figure  4.  Test curves of acceleration and shock wave overpressure peak

    为分析不同速度战斗部爆炸时冲击波场的分布规律,建立分布式动爆试验测试系统,采用球形裸装药,装药量为1.2 kg。比例距离R分别为4.71、9.41和14.12 m/kg1/3,战斗部着靶速度v0分别为0、535和980 m/s,火炮火药发射获得着靶速度,距目标炸点50 m处顺序放置4台天幕靶,每台间隔20 m,通过区截法获取弹丸在天幕靶处的飞行速度,再结合炮口靶获取的出炮口速度以及制式弹的外弹道模型,计算出目标炸点处弹丸速度。采用模块装药,速度为535和980 m/s弹丸的装药分别为3×B模块和6×B模块,每个B模块装药2.35 kg,引信采用瞬触发引信。战斗部与地靶平面的水平夹角β为45°,测点分布如图5所示,以爆炸中心为原点,建立三维坐标系Oxyz,分别在地靶平面距爆心5、10、15 m处安装冲击波超压测试装置。在Oxy平面上(即地靶平面),爆心到测试点的连线与x轴的夹角为θ,受试验条件限制,共24个测点。理论上测点越多,对爆炸冲击波的重建越有利。

    图  5  战斗部爆炸测点分布
    Figure  5.  Measuring points of blasting warhead explosion

    图6为比例距离R=4.71 m/kg1/3时不同方向上测得的冲击波超压时域曲线。由图6可知,以静爆条件下(v0=0 m/s)的冲击波超压时域曲线为参照,当θ=0°时,动爆冲击波超压远高于静爆冲击波超压;当θ增大到45°时,动爆冲击波超压有所下降,但是仍然高于静爆冲击波超压;当θ增大到90°时,动爆冲击波超压继续降低,v0对冲击波超压的影响变小,不同速度战斗部爆炸的冲击波超压趋于一致,与静爆冲击波压力相当;随着θ的进一步增大,动爆冲击波压力进一步降低,开始低于静爆冲击波压力;当θ增大到180°时,测点处于与战斗部速度完全相反的方向,为压力最低点,此时的冲击波超压远低于静爆冲击波超压。此外,在动爆冲击波超压高于静爆冲击波超压的方向上,即θ在0°~90°和270°~360°范围内时,v0越大,压力越高,冲击波到达时间越短;相反地,在动爆冲击波超压低于静爆冲击波超压的方向上,即θ在90°~270°范围内时,v0越大,压力越低,冲击波到达时间越长。

    图  6  不同方向的冲击波超压时域曲线
    Figure  6.  Shock wave overpressure-time curves in different directions

    读取冲击波的超压峰值,得到不同速度战斗部的爆炸冲击波超压峰值对比曲线,如图7所示。由图7可知:(1)比例距离相同时,战斗部爆炸的冲击波超压峰值随θ的增大近似呈余弦衰减,当θ=0°时,超压峰值最大,θ=180°时,超压峰值最小,并且战斗部着靶速度v0越大,超压峰值衰减得越快;(2)以静爆冲击波超压峰值为参照,在与战斗部速度方向相同的区域(0°~90°和270°~360°),动爆冲击波存在较大的压力升,超压峰值大于静爆状态下的超压峰值,使得冲击波场呈现出局部高压区,而在与战斗部速度方向相反的区域(90°~270°),则存在较大的压力降。由此,可将运动战斗部的爆炸冲击波场分为压力升和压力降两个区域,分界点大约在θ=90°处。

    图  7  不同速度战斗部的冲击波超压峰值曲线
    Figure  7.  Shock wave overpressure peak curves of the blasting warheads with different velocities

    根据战斗部动爆冲击波超压峰值pp,d与战斗部相对靶平面的速度v0及静爆冲击波超压峰值pp,s的关系[4]

    pp,dpp,s=(1+0.31+Rv0c0cosθ)2
    (4)

    战斗部在靶平面的投影速度v0为:

    v0=v0cosβ
    (5)

    由式(1)、(4)、(5)联合可得战斗部动爆冲击波超压峰值的经验公式计算值,实测战斗部动爆冲击波超压峰值pp,e与经验公式计算结果pp,d对比如表1所示,其中ε=(pp,epp,d)/pp,e。由表1可知,实测战斗部静爆冲击波超压峰值与经验公式计算结果一致性较好,除个别点外,实测战斗部动爆冲击波超压峰值与理论值较为接近,相对误差小于20%,且相对误差随着战斗部速度的增大而增大。此外,越靠近战斗部运动方向轴线(θ=0°和θ=180°)的实测冲击波超压峰值相对误差越大,垂直于战斗部运动方向轴线(θ=90°和θ=270°)的实测冲击波超压峰值相对误差较小。

    表  1  试验结果与理论值对比
    Table  1.  Comparison between experimental and theoretical results
    R/(m·kg−1/3θ/(°)v0=0 m/sv0=535 m/sv0=980 m/s
    pp,e/kPapp,d/kPaε/%pp,e/kPapp,d/kPaε/%pp,e/kPapp,d/kPaε/%
    4.710321.39317.351.26434.86355.5718.23483.16389.0219.48
    45343.49317.357.61384.34344.1510.46392.65367.286.46
    90342.92317.357.45285.10317.35−11.31288.32317.35−10.07
    135335.52317.355.41270.90291.64−7.65258.87271.08−4.72
    180320.94317.351.12228.45281.30−23.14187.89252.98−34.64
    225294.20317.35−7.87261.93291.64−11.34215.23271.08−25.95
    270296.57317.35−7.01285.98317.35−10.97294.97317.35−7.59
    315319.73317.350.74368.50344.156.61402.27367.288.70
    9.410112.04104.866.41127.36111.7012.30149.08117.5421.15
    45110.61104.865.19110.62109.670.86118.86113.764.30
    90100.10104.86−4.7696.52104.86−8.6596.21104.86−9.00
    135103.81104.86−1.0288.96100.16−12.5986.6996.34−11.13
    180106.88104.861.8990.3598.25−8.7575.8892.91−22.44
    225102.12104.86−2.6994.81100.16−5.6488.0696.34−9.40
    27098.28104.86−6.7099.25104.86−5.6698.10104.86−6.90
    315111.00104.865.53118.45109.677.41123.47113.767.86
    14.12063.8763.041.3174.8765.8512.0578.2368.2412.77
    4563.2363.040.3066.3565.022.0169.5666.694.12
    9057.7263.04−9.2159.1663.04−6.5559.0363.04−6.78
    13559.0863.04−6.7057.7561.08−5.7850.6359.48−17.48
    18060.5163.04−4.1751.7260.28−16.5645.7358.04−26.91
    22557.7563.04−9.1653.6961.08−13.7748.0659.48−23.77
    27057.7863.04−9.1059.5263.04−5.9257.5263.04−9.58
    31564.7363.042.6169.4265.026.3473.9866.699.85
    下载: 导出CSV 
    | 显示表格

    表1中的冲击波超压峰值为插值点,利用MATLAB数据处理软件的薄板样条插值方法(thin-plate-spline interpolation)对实测数据进行插值(该插值方法可以使得三维超压曲面弯曲能量最小),得到战斗部速度分别为0、535、980 m/s的爆炸冲击波超压峰值场分布和等压曲线,如图810所示。由图8可知,战斗部静爆冲击波超压峰值在各个方向基本相同;由图910可知,战斗部动爆冲击波超压峰值在战斗部运动速度方向增强,在战斗部运动相反方向减弱,且战斗部速度越大,增强和减弱的程度越大。

    图  8  战斗部速度为0 m/s的冲击波超压峰值场分布
    Figure  8.  Shock wave overpressure field for the blasting warhead with the velocity of 0 m/s
    图  9  战斗部速度为535 m/s的冲击波超压峰值场分布
    Figure  9.  Shock wave overpressure field for the blasting warhead with the velocity of 535 m/s
    图  10  战斗部速度为980 m/s的冲击波超压峰值场分布
    Figure  10.  Shock wave overpressure field for the blasting warhead with the velocity of 980 m/s

    提出了一种基于地震波可靠触发的战斗部空中爆炸冲击波超压测试方法,并对速度为0、535和980 m/s的战斗部空中爆炸冲击波进行了测试研究,通过测试结果和经验公式计算值的对比分析,以及重建的战斗部动爆冲击波超压三维可视化模型,可以得出以下结论:

    (1)本文中提出的测试方法能可靠获取战斗部动爆冲击波超压峰值;

    (2)战斗部动爆冲击波超压峰值在战斗部运动速度方向增强,在战斗部运动相反方向减弱,且战斗部速度越大,增强和减弱的程度越大。

  • 图  1  低碳钢和泡沫铝的塑性应力应变曲线

    Figure  1.  Plastic stress-strain curves of mild steel and aluminum foam

    图  2  夹芯梁有限元模型

    Figure  2.  The finite element model for the sandwich beam

    图  3  实验装置

    Figure  3.  Experimental apparent

    图  4  不同冲击次数时夹芯梁最终变形对比

    Figure  4.  Comparison of permanent deflections in different impact numbers

    图  5  冲击力时程曲线对比

    Figure  5.  Comparison of time histories of impact force

    图  6  最终挠度的数值仿真结果与实验结果的对比

    Figure  6.  Comparison of permanent deflections between numerical simulation and impact test

    图  7  夹芯梁重复冲击变形过程

    Figure  7.  Deformation process of the sandwich beam under repeated impacts

    图  8  泡沫金属夹芯梁受力和变形情况

    Figure  8.  Stress state and deformation of the sandwich beam

    图  9  夹芯梁上、下面板中点挠度时程曲线

    Figure  9.  Time histories of deflections of the front and back faces of the sandwich beam

    图  10  冲击力-位移曲线

    Figure  10.  Force-displacement curves

    图  11  泡沫金属夹芯梁各部分的变形能时程曲线

    Figure  11.  Time histories of energy absorption for different parts of the metal foam sandwich beam

    图  12  不同冲击次数下夹芯梁各部分的能量吸收

    Figure  12.  Energy absorption for different parts of the metal foam sandwich beam at different impact numbers

    图  13  不同冲击次数下的回弹特性

    Figure  13.  Rebound characteristics at different impact numbers

    图  14  不同冲击次数下的能量吸收效率

    Figure  14.  Energy absorption efficiency at different impact numbers

    图  15  不同冲击次数下的回弹因数

    Figure  15.  Rebound coefficient at different impact numbers

  • [1] ZHU L. Dynamic inelastic behaviour of ship plates in collision [D]. Glasgow: University of Glasgow, 1990: 45−56.
    [2] ZHU L. Faulkner D. Damage estimate for plating of ships and platforms under repeated impacts [J]. Marine Structures, 1996, 9(9): 697–720. DOI: 10.1016/0951-8339(95)00018-6.
    [3] ZHU L, CAI W, CHEN M S, et al. Dynamic analysis of ship plates under repeated ice floes impacts based on a simplified ship-ice collision model [C]//Proceedings of the 28th International Ocean and Polar Engineering Conference. Sapporo, Japan: International Society of Offshore and Polar Engineers, 2018: 1718−1723.
    [4] ZHU L, SHI S Y, JONES N. Dynamic response of stiffened plates under repeated impacts [J]. International Journal of Impact Engineering, 2018, 117: 113–122. DOI: 10.1016/j.ijimpeng.2018.03.006.
    [5] 杨宝, 汤立群, 刘逸平, 等. 冲击条件下泡沫铝的细观变形特征分析 [J]. 爆炸与冲击, 2012, 32(4): 399–403. DOI: 10.11883/1001-1455(2012)04-0399-05.

    YANG B, TANG L Q, LIU Y P, et al. Meso deformation characteristics analysis of aluminum foam under impact [J]. Explosion and Shock Waves, 2012, 32(4): 399–403. DOI: 10.11883/1001-1455(2012)04-0399-05.
    [6] 王鹏飞, 徐松林, 李志斌, 等. 高温下轻质泡沫铝动态力学性能实验 [J]. 爆炸与冲击, 2014, 34(4): 433–438. DOI: 10.11883/1001-1455(2014)04-0433-06.

    WANG P F, XU S L, LI Z B, et al. An experimental study on dynamic mechanical property of ultra-light aluminum foam under high temperatures [J]. Explosion and Shock Waves, 2014, 34(4): 433–438. DOI: 10.11883/1001-1455(2014)04-0433-06.
    [7] YU J L, WANG X, WEI Z G, et al. Deformation and failure mechanism of dynamically loaded sandwich beams with aluminum-foam core [J]. International Journal of Impact Engineering, 2003, 28(3): 331–347. DOI: 10.1016/S0734-743X(02)00053-2.
    [8] YU J L, WANG E H, LI J R, et al. Static and low-velocity impact behavior of sandwich beams with closed-cell aluminum-foam core in three-point bending [J]. International Journal of Impact Engineering, 2008, 35(8): 885–894. DOI: 10.1016/j.ijimpeng.2008.01.006.
    [9] TAN Z H, LUO H H, LONG W G, et al. Dynamic response of clamped sandwich beam with aluminum alloy foam core subjected to impact loading [J]. Composites Part B: Engineering, 2013, 46: 39–45. DOI: 10.1016/j.compositesb.2012.10.044.
    [10] 敬霖, 王志华, 赵隆茂, 等. 撞击载荷下泡沫铝夹芯梁的塑性动力响应 [J]. 爆炸与冲击, 2010, 30(6): 561–568. DOI: 10.11883/1001-1455(2010)06-0561-08.

    JING L, WANG Z H, ZHAO L M, et al. Dynamic plastic response of foam sandwich beams subjected to impact loading [J]. Explosion and Shock Waves, 2010, 30(6): 561–568. DOI: 10.11883/1001-1455(2010)06-0561-08.
    [11] JING L, WANG Z H, NING J G, et al. The dynamic response of sandwich beams with open-cell metal foam cores [J]. Composites Part B: Engineering, 2011, 42(1): 1–10. DOI: 10.1016/j.compositesb.2010.09.024.
    [12] DESHPANDE V S, FLECK N A. Isotropic constitutive models for metallic foams [J]. Journal of the Mechanics and Physics of Solids, 2000, 48(6−7): 1253–1283. DOI: 10.1016/S0022-5096(99)00082-4.
    [13] QIU X M, DESHPANDE V S, FLECK N A. Finite element analysis of the dynamic response of clamped sandwich beams subject to shock loading [J]. European Journal of Mechanics: A/Solids, 2003, 22(6): 801–814. DOI: 10.1016/j.euromechsol.2003.09.002.
    [14] TILBROOK M T, DESHPANDE V S, FLECK N A. Underwater blast loading of sandwich beams: regimes of behaviour [J]. International Journal of Solids and Structures, 2009, 46(17): 3209–3221. DOI: 10.1016/j.ijsolstr.2009.04.012.
    [15] JING L, SU X Y, CHEN D, et al. Experimental and numerical study of sandwich beams with layered-gradient foam cores under low-velocity impact [J]. Thin-Walled Structures, 2019, 135: 227–244. DOI: 10.1016/j.tws.2018.11.011.
    [16] QIU X M, DESHPANDE V S, FLECK N A. Impulsive loading of clamped monolithic and sandwich beams over a central patch [J]. Journal of the Mechanics and Physics of Solids, 2005, 53(5): 1015–1046. DOI: 10.1016/j.jmps.2004.12.004.
    [17] QIN Q H, WANG T J. Low-velocity heavy-mass impact response of slender metal foam core sandwich beam [J]. Composite Structures, 2011, 93(6): 1526–1537. DOI: 10.1016/j.compstruct.2010.11.018.
    [18] QIN Q H, WANG T J. Low-velocity impact response of fully clamped metal foam core sandwich beam incorporating local denting effect [J]. Composite Structures, 2013, 96: 346–356. DOI: 10.1016/j.compstruct.2012.09.024.
  • 期刊类型引用(0)

    其他类型引用(3)

  • 加载中
图(15)
计量
  • 文章访问数:  975
  • HTML全文浏览量:  376
  • PDF下载量:  216
  • 被引次数: 3
出版历程
  • 收稿日期:  2020-06-16
  • 修回日期:  2021-03-17
  • 网络出版日期:  2021-07-06
  • 刊出日期:  2021-07-05

目录

/

返回文章
返回