en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
通讯作者:

刘文洋,E-mail:liuwenyang@hnu.edu.cn

中图分类号:O383

文献标识码:A

文章编号:1672-6553-2022-20(2)-078-10

DOI:10.6052/1672-6553-2021-062

参考文献 1
Xue Z,Hutchinson J W.Preliminary assessment of sandwich plates subjected to blast.International Journal of Mechanical Sciences,2003;45:687~705
参考文献 2
Aune V,Fagerholt E,Langseth M,et al.A shock tube facility to generate blast loading on structures.International Journal of Protective Structures,2016,7(3):340~66
参考文献 3
Krauthammer T,Altenberg A.Negative phase blast effects on glass panels.International Journal of Impact Engineering,2000,24(1):1~17
参考文献 4
Aune V,Fagerholt E,Hauge K O,et al.Experimental study on the response of thin aluminium and steel plates subjected to airblast loading.International Journal of Impact Engineering,2016,90:106~121
参考文献 5
Galiev S U.Experimental observations and discussion of counterintuitive behavior of plates and shallow shells subjected to blast loading.International Journal of Impact Engineering,1996,18(7-8):783~802
参考文献 6
Symonds P S,Yu T X.Counterintuitive behavior in a problem of elastic-plastic beam dynamics.Journal of Applied Mechanics,1985,52(3):517~522
参考文献 7
Yu T X.Chapter 9:Elastic effects in the dynamic plastic response of structures.Proceedings of the Third International Symposium on Structural Crashworthiness,1993,295~332
参考文献 8
Li Q M,Zhao L M,Yang G T.Experimental results on the counter-intuitive behaviour of thin clamped beams subjected to projectile impact.International Journal of Impact Engineering,1991,11(3):341~348
参考文献 9
吴桂英,李鑫,周稳.空中爆炸载荷作用下铝合金圆板的反直观响应研究.振动与冲击,2010,29(12):166~168,246(Wu G V;Li X;Zhou W.The anomalous response of the aluminum alloy circular plate subjected to the air blast load.Journal of Vibration and Shock,2010,29(12):166~168,246(in Chinese))
参考文献 10
吴桂英,秦冬祺,杨桂通.在冲击载荷作用下弹塑性圆板的反直观动力行为数值分析.固体力学学报,2004(4):389~393(Wu G Y;Qin D Q;Yang G T.The anomalous response of elastoplastic circular plates under impulsive load.Acta Mechanica Solida Sinica,2004(4):389~393(in Chinese))
参考文献 11
Kinney G F,Graham K J.Explosive shocks in air.Berlin:Springer,1985
参考文献 12
Friedlander F G.The diffraction of sound pulses I.Diffraction by a semi-infinite plane.Proceedings of the Royal society A:Mathematical,Physical and Engineering Sciences,1946,186(1006):322~44
参考文献 13
Swisdak J M.Simplified kingery airblast calculations.Minutes of the Twenty Sixth DOD Explosives Safety Seminar,DTIC Document,Maryland.
参考文献 14
Nassr A A.experimental and analytical study of the dynamic response of steel beams and columns to blast loading.McMaster-Civil Engineering,2012
参考文献 15
Rigby S E,Tyas A,Bennett T,et al.The negative phase of the blast load.International Journal of Protective Structures,2014,5(1):1~19
参考文献 16
Karlos V,Solomos G,Larcher M.Analysis of the blast wave decay coefficient using the Kingery-Bulmash data.International Journal of Protective Structures,2016,7(3):409~429
参考文献 17
Voce E.The relationship between stress and strain for homogenous deformation.Journal of the Japan Institute of Metals and Materials,1948,74:537~562
目录contents

    摘要

    薄板在爆炸冲击载荷作用下的动力行为复杂,具有大变形、瞬时性、高度非线性等特点.当冲击波冲量较小时,薄板容易发生异常动力响应.然而,爆炸冲击波特征对薄板异常动力响应的影响尚不明晰.为此,本文采用经实验数据验证的自主开发爆炸载荷模型,首先成功地模拟了方形薄板的反直观异常动力响应;然后通过数值仿真研究了冲击波负压描述方法、正压衰减系数等对薄板异常动力响应的影响;最后,发现了在距离当量比一定时,随着爆炸距离的增加薄板依次发生继发直观、反直观、混合、始发直观四种行为模式,并分析了不同行为模式下响应历程与塑性耗散能之间的关系.

    Abstract

    The dynamic behavior of thin plate under impulsive loading is of high complexity, characterized by large deformation, instantaneity, and severe nonlinearity. When a thin plate is subjected to a shock wave of low impulse, abnormal dynamic response is prone to occur. However, effects of loading characteristics on the abnormal dynamic response of thin plate are not clear. An impulsive loading model validated by experimental data is successfully used to simulate the counter-intuitive response of a thin square plate, and then the effects of negative pressure and decaying parameter on the abnormal dynamic response are studied. Finally, when the stand-off distance increases with a fixed scaled distance, four different modes-secondary intuitive, counter-intuitive, mixed, and initial intuitive modes-are found. The relationship between the different dynamic behavior modes and the plastic dissipation energy is examined.

  • 引言

  • 薄板在爆炸冲击载荷作用下的变形预测通常只关注爆炸冲击载荷的正相冲击压力(后简称正压),而忽略负相压力(后简称负压).这在爆炸点距离薄板较近或者当量较大的情况下是合理的,此时薄板的变形失效过程受爆炸冲击正压主导[1, 2].由于爆炸冲击的正压极值远大于负压极值,薄板变形方向与冲击波方向一致.但是,当爆炸点距离薄板较远或者当量较小的情况下,忽略爆炸冲击负压则不合理.在反射超压峰值与负压峰值大小相近时,负压对结构动力响应影响不可忽略[3].V.Aune等人发现爆炸冲击载荷的负压阶段会改变薄板最终构型[4].余同希、李庆明等人发现爆炸冲击载荷下的薄板弹塑性动力响应可分成两类[5-8].若中心的永久偏转方向与冲击波方向相同,则称之为直观变形模式;若中心的永久偏转方向与冲击波方向相反,则称之为异常(反直观)变形模式.

  • 薄板异常动力响应受载荷、边界条件、材料、几何形状的影响 [9, 10].大多数薄板异常动力响应的研究采用简化载荷如矩形或三角形脉冲,与真实爆炸冲击载荷相差较大;若采用商业软件中内嵌的爆炸冲击经验模型,则难以分离载荷特征变量,从而无法分析其对异常动力响应的影响.本文使用经实验数据验证的自主开发载荷模型,从爆炸载荷冲击波的特征入手,分析了影响方形薄板动力行为的因素;通过控制距离当量比,发现了方形薄板变形模式的变化规律,并进一步分析了薄板异常动力行为与塑性耗散能之间的关系.

  • 1 爆炸冲击载荷压强-时间历程

  • 图1 描述了爆炸冲击压力随时间的变化.波阵面在起爆后经过时间ta抵达,冲击波压强瞬时达到入射超压极值p maxso,t+d是冲击波正压持续时间.随着冲击波的扩散,压力指数衰减回大气压P0.由于空气过度膨胀后分子间距减小从而产生负压.t-d为负压持续时间,负压极值相比于正压极值较小,而负压持续时间远长于正压持续时间.

  • Kinney和Graham将大量球形、半球形冲击实验的冲击波压强数据通过距离当量比(Z=R/W 1/3) 和等效TNT质量(W)作为自变量组成的高阶多项式对载荷特征参数进行曲线拟合[11].

  • 图1 爆炸冲击载荷压强随时间变化示意图

  • Fig.1 Typical pressure-time history for a blast wave

  • 1.1 冲击波正压

  • 正压压强-时间历程可由修正的Friedlander[12]方程描述,即

  • p(t)=p00<t<tap0+prα1-t-tatd+exp-bt-tatd+,ta<t<ta+td+
    (1)
  • 其中,当冲击波遇到与其方向不平行的结构时发生斜反射,冲击压力增强,形成斜反射压力p.

  • 在Friedlander方程中,用等效TNT质量W与距离当量比Z作为输入参数,来定义描述正相载荷的各项参数.冲击波到达时间ta、正压持续时间td+、入射超压极值psomax、反射压强p ref的计算依据Swisdak [13]提出的拟合公式获得,正相冲量ir+依据Nassr [14]提出的拟合公式获得.而Friedlander方程中构件所受斜反射压力p,由受冲击构件表面法线与爆炸冲击波之间的入射角α以及入射超压极值psomax和反射超压p ref决定

  • prα=psomax1+cosα-2cos2α+prefcos2α
    (2)
  • 在公式(1) 中唯一的未知量为衰减系数b.通过正压压强对时间td+积分得到关于冲量irα+的隐式方程:

  • irα+=tahd+td prα1-t-tatd+exp-bt-tatd+dt=prαtd+b2(b-1+exp(-b))
    (3)
  • 将正相冲量拟合公式代入式(3)可迭代求解衰减系数b.衰减系数b是关于爆炸距离R与等效TNT质量W的函数.该值在一定的R/W区域内呈现高度非线性.图2直观描述了衰减系数b随爆炸距离R与等效TNT质量W的变化.可以发现当爆炸距离R恒定时,等效TNT质量在10g以上衰减系数b的变化较为平缓;而等效TNT质量在0g-10g之间,衰减系数b随等效TNT质量减小急剧增大,并且在等效TNT质量接近0时又快速减小,在该区间内呈现出“∧”型.而在等效TNT质量恒定的情况下,衰减系数b随爆炸距离R的变化相对平稳.这表明爆炸距离与等效TNT质量对衰减系数b的影响程度不同,等效TNT质量对衰减系数b的影响更加显著.

  • 图2 衰减系数b变化趋势

  • Fig.2 Variation trend of decay coefficient b

  • 1.2 冲击波负压

  • 关于爆炸冲击载荷的负压,其描述方式有两种[15]:(1) 双线性描述法;(2) 三次曲线描述法.对于双线性描述法,爆炸冲击载荷的负压被分成两部分,压强随时间线性变化.双线性描述法压强随时间的分段函数如下.其中α为负压时间比例参数,控制负压压强上升时间在整个负压持续时间的占比.

  • p(t)=p0-prmint-ta+td+αtdlin-, ta+td+<t<ta+td++αtdlin-p0-prmint-ta+td++αtdlin (1-α)tdlin-, ta+td++αtdlin-<t<ta+td++tdlin-p0, t>ta+td++tdlin-
    (4)
  • 三次曲线描述法则采用一条三次曲线拟合爆炸冲击载荷的负压变化过程.其中负压上升时间大致占负压持续时长的三分之一.

  • p(t)=p0-prmin274t-ta+td+tdcub-1-t-ta+td+tdcub-, ta+td+<t<ta+td++tdcub-p0, t>ta+td++tdcub-
    (5)
  • 双线性和三次描述法负压参数取自文献[16],负压极值prmin与负压冲量ir-是根据拟合公式计算.

  • pr,min=100 ,Z<0.80379.433Z-1.05 ,Z0.803
    (6)
  • ir-=480W13 ,Z<0.60810-0.87731log(Z)+2.49145 ,Z0.608
    (7)
  • 将式(4)与式(5)冲击波负压分别对时间积分得到负压冲量,再将负压冲量拟合公式(7) 代入,可以求得双线性负压持续时间tdlin-与三次曲线负压持续时间tdcub-.

  • tdlin-=2ir-prmin
    (8)
  • tdcub-=169ir-prmin
    (9)
  • 2 仿真试验模型

  • 爆炸冲击载荷下薄板异常动力响应对外界条件较为敏感.其中不仅爆炸载荷等效TNT质量W、爆炸距离R以及爆炸冲击载荷模型的各项参数对薄板异常动力响应的仿真结果有着显著影响,薄板的几何形状、材料属性、边界约束条件等因素对仿真结果也会造成较大的影响.为模拟试验工况,建立了较为复杂的几何模型.材料本构为修正Johnson-Cook模型.使用ABAQUS用户子程序实现了自主开发爆炸冲击载荷的加载以及材料本构的赋予.

  • 2.1 几何模型

  • 金属薄板两侧用钢制夹具夹持并以螺栓固定.金属薄板尺寸为400mm×400mm×0.8mm,两侧钢制夹具宽度为50mm,螺栓直径为20mm,金属薄板的实际受载区域为300mm×300mm.爆炸载荷加载在金属薄板正中心,整体结构中心对称,为增加计算效率,采用四分之一建模法.钢制夹具外围施加固支约束.其中薄板为三维连续壳单元(S4R);钢制夹具与圆柱螺栓均为三维实体单元(C3D8R),如图3所示.金属薄板采用1050A-H14铝材,材料属性如表1.钢制夹具材料属性如表2.

  • 为模拟实际工况中螺栓螺帽与螺母为薄板两侧的钢制夹具提供预紧力,在模型两侧的钢制夹具外表面螺栓周围一圈施加均布压力500MPa,螺栓上下表面与夹具绑定约束.为模拟金属薄板与夹具接触面之间滑动摩擦力,设置有切向行为表面接触,静摩擦系数为0.15,动摩擦系数为0.1.爆炸载荷等效TNT质量为41.2g,距离铝板中心625mm,四周无反射面.

  • 图3 几何模型

  • Fig.3 Geometric model

  • 2.2 材料本构

  • 为模拟材料弹塑性大变形,将变形率加法分解为弹性和塑性两部分

  • D=De+Dp
    (10)
  • 式中,De为弹性变形率,Dp为塑性变形率.弹性变形率与柯西应力之间的本构关系为

  • De=1+νEσΔJ-νEtrσΔJI
    (11)
  • 式中,E为杨氏模量,ν为泊松比, σΔJ为柯西应力的Jaumann率,I为单位张量.对于各向同性材料,屈服函数采用

  • fσ,εp,T=σeq-σyεp,T
    (12)
  • 式中,σeq为等效应力,εp为等效塑性应变,T为温度.

  • σeq=32σ':σ'
    (13)
  • σ'=σ-13tr(σ)I
    (14)
  • 式中,σ′为应力偏量,应变增量与应力偏量有关.材料屈服应力函数采用修正Johnson-Cook模型[17]

  • σy=σ0+i=12 Q11-exp-Ciεp1+ε˙pε˙0pc1-T*m
    (15)
  • 式中,σ0为初始屈服应力,Q 1,Q 2,C 1,C 2为材料硬化参数,p 0为参考塑性应变率,T*的表达式为

  • T*=T-TrTm-Tr
    (16)
  • 式中,Tr为环境温度,Tm为材料熔点.假设塑性势函数与屈服函数重合,塑性变形率可通过相关联流动法则表达为

  • DP=ε˙pfσ=32ε˙pσ'σeq
    (17)
  • 金属材料在爆炸载荷下迅速发生剧烈塑性变形产生温升,可假设为绝热过程.温升速率的表达式为

  • T˙=χρCpσ:Dp=χρCpσeqε˙p
    (18)
  • 式中,χ为Taylor-Quinney系数,ρ为密度,Cp为比热.

  • 3 仿真验证与结果讨论

  • 3.1 仿真验证

  • 爆炸冲击载荷各项特征参数对薄板反直观异常动力响应有着诸多影响.本文采用经实验数据验证的自主开发爆炸冲击载荷模型对薄板反直观异常动力响应进行模拟研究.在研究过程中发现爆炸载荷负压采用三次曲线描述法时薄板中心挠度到达负方向最大值时间早于双线性描述法,与实验结果差异较大.因此后续仿真结果除特别说明外爆炸载荷负压均采用双线性描述法.图4对比了实验与仿真的结果.其中两次实验的外界条件一致,但薄板反直观行为的敏感性导致结果略有差异.图4(a)中仿真预测的薄板正向变形历程与实验整体一致,回弹曲线位于两次实验曲线中间,薄板达到正方向最大值的时间点完全一致,达到负方向最大值的时间点吻合,但仿真回弹最大挠度略高于实验结果,这可能是由于负压冲量估值略大于实际值.图4(b)中薄板负方向挠度最大值仿真与实验基本一致.通过对比仿真与试验的薄板变形历程与最终构型,验证了爆炸冲击载荷模型的准确性.

  • 图4 实验与仿真结果对比

  • Fig.4 Comparison of experimental and simulation

  • 薄板在爆炸冲击载荷下的变形过程具有“先边缘,后中心”的特性.图5为等时间间隔薄板变形过程.在1ms时波阵面刚抵达,薄板边缘部分先发生变形,而中间部分仍保持平直,表明薄板的变形是由边缘部分向中心部分逐渐扩散的.在3ms时薄板中心外围区域先回弹至负方向,中心位置最晚到达负方向极值.表明回弹过程变形也是从板边缘向中心逐渐发展.

  • 图5 薄板等时间间隔变形图

  • Fig.5 Equal time interval deformation diagram of thin plate

  • 表1 1050A-H14铝材料属性

  • Table1 Material properties of 1050A-H14aluminum

  • 表2 钢材料属性

  • Table2 Material properties of steel

  • 3.2 负压对薄板动力响应的影响

  • 由于负压压强极值远低于正压压强极值,研究结构受爆炸冲击载荷的动力行为时,冲击波负压常常被忽略.这种简化适用于板厚较厚、刚度较强情况下的动力响应仿真,但对于较为敏感的薄板结构,负压会对仿真结果产生不可忽视的影响.图6(a)、(b)分别为薄板在有无负压情况下中心挠度变化历程的仿真结果.在1.2ms之前,由于爆炸冲击波处于正压阶段,两次仿真的结果完全一致.在1.2ms后,有负压作用下的薄板开始回弹,并发生了反直观异常动力响应;没有加载负压的薄板则发生直观行为.结果表明爆炸冲击负压对薄板异常动力响应的发生至关重要.因此,下文对爆炸冲击载荷负压的影响进行分析.

  • 3.2.1 负压描述法的影响

  • 爆炸载荷负压段使得准确预测薄板反直观行为结果更加复杂.图7(a)对比了在同一工况下分别采用双线性与三次曲线两种负压描述方式的薄板中心挠度变化历程.结果表明,三次曲线加载下薄板回弹速度更快,中心挠度先于双线性负压加载到达负方向极值.造成双线性负压加载下薄板回弹时间延长的原因与负压持续时间有关,当负压冲量和负压极值一致时,双线性加载方式负压持续时间略长于三次曲线加载,两者相差12.5%.爆炸过程的高度非线性放大了两者的差异,两种加载方式薄板中心挠度到达负方向极值时间相差22.2%.在双线性、三次曲线两种不同负压加载模式下薄板反直观行为最终挠度值基本相等,这是由于爆炸载荷负压冲量相等.

  • 图6 负压对薄板动力行为的影响

  • Fig.6 Influence of negative pressure on dynamic behavior

  • 使用双线性负压加载方式时,描述负压过程中到达负压极值所需时间在整个负压时间段所占比例的参数α对薄板反直观行为的数值模拟结果有着较大的影响.图7(b)对比了参数α=0.4与α=0.8两种情况下薄板中心挠度变化历程.当参数α=0.8时,薄板中心挠度达到负方向极值的时间晚于参数α=0.4,而薄板中心挠度负方向极值两次模拟结果一致.在冲击波负压段冲量一定的情况下,薄板中心挠度极值基本不变.负压压强随时间历程的改变引起薄板中心挠度时间历程的变化.

  • 图7 负压描述法对薄板中心挠度变化历程影响

  • Fig.7 Influence of negative pressure profile on variation history of central deflection of thin plate

  • 图8进一步表明了不同的负压时间比例参数α下,薄板中心挠度负方向极值和薄板中心挠度负方向极值到达时间的变化趋势.薄板中心挠度最高最低值之差仅有2.14mm,不超过平均值的7%,由此可知,参数α的变化几乎不会影响薄板回弹至负方向的极值.但是,参数α对薄板回弹时间历程有着显著影响.随着负压极值到达时间的延后,薄板中心挠度到达最低值的时间也在延后,其最高最低值之差达到平均值的21%.

  • 综上,爆炸载荷负压段描述方式仅对薄板变形回弹时间历程产生影响,对负方向挠度极值没有影响.

  • 3.2.2 正压冲量与衰减系数b

  • 不同学者提出了多种冲击波正压冲量计算公式.衰减系数b受正压冲量控制,并且两者呈非线性关系.选取了两种基于试验数据的正压冲量拟合公式[13, 15]分别进行仿真实验.图9对比了两种正压冲量计算方法薄板中心挠度时间历程的结果.两者正压冲量相差36%,薄板中心挠度正方向极值差10mm,均在1.2ms达到正方向极值.在1.2ms后两者开始回弹过程,到达负方向极值时间相差2ms,薄板中心挠度负方向极值差8mm.由此可以发现,在薄板向正方向变形的过程中,变形持续时间并未受到正压冲量及衰减系数变化的影响,而这两者的变化引起了薄板挠度正负方向极值与回弹时间历程的不同.

  • 图8 负压时间比例参数α的影响

  • Fig.8 Influence of negative pressure time parameter α

  • 图9 衰减系数对薄板动力响应的影响

  • Fig.9 Influence of decay coefficient on dynamic response of thin plate

  • 4 薄板动力响应与模式转变

  • 大多数研究采用距离当量比Z衡量冲击强度并通过调整Z值来研究结构在不同爆炸冲击强度下的动力行为.但是,薄板在爆炸冲击载荷下的动力行为复杂,其行为模式不单一取决于正负压极值.仅通过调整距离当量比Z难以探明薄板动力行为与爆炸冲击载荷的关系.因此,需要进一步讨论薄板动力行为与其他冲击波控制量的关系.

  • 4.1 薄板动力响应模式

  • 保持距离当量比Z不变,则冲击波正负压极值即冲击强度恒定.通过改变爆炸距离R,可以实现在冲击波强度不变的情况下改变冲击波正负压持续时间,冲击波冲量也随之变化.在此条件下,冲击波冲量与爆炸距离R呈正相关.

  • Z=1.809m/kg1/3,R由375mm逐渐递增至750mm,不同爆炸距离薄板变形模式总结如表3所示.在冲击强度不变的情况下,薄板动力行为仍可以表现出四种模式:模式一、继发直观模式;模式二、反直观模式;模式三、始发直观模式;模式四、混合模式.其中,继发直观模式表现为薄板到达初次变形极值后向负方向发生反向屈曲,在到达负向极值后向正方向再次发生反向屈曲并最终停留在正方向一侧,如图10 (a)所示.而始发直观模式下,薄板在初次正向变形后不会发生反向屈曲,直接停留在正方向一侧.图10(b)为继发直观模式薄板最终构型.通过与图6(a)对比可知,相比于始发直观模式,两者仅在正向变形程度上有所区别,整体构型基本一致.图11(a)为混合变形模式薄板中心点(A点)与靠近边缘处一点(B点)的挠度变化历程.薄板边缘区域在极值点反向屈曲后停留在负方向,发生反直观行为;而中心区域始终在正方向运动,发生直观行为.图11 (b)为薄板混合模式中线构型,其中心区域向正方向凸出,边缘区域向负方向凹入.混合模式表明薄板的不同区域可以同时发生直观行为与反直观行为.

  • 综上,继发直观模式发生多于一次的极值点反向屈曲;反直观模式发生一次极值点反向屈曲;混合模式薄板边缘区域发生一次极值点反向屈曲;始发直观模式没有发生极值点反向屈曲.

  • 表3 距离当量比Z=1.809m/kg1/3时不同爆炸距离R下薄板变形模式

  • Table3 Deformation mode of thin plate under different explosion distance R when scaled distance Z=1.809m/kg1/3

  • 图10 薄板继发直观变形模式

  • Fig.10 Secondary intuitive deformation mode of thin plate

  • 图11 薄板混合变形模式

  • Fig.11 Mixed deformation mode of thin plate

  • 表4 距离当量比Z=1.6m/kg1/3时不同爆炸距离R下薄板变形模式

  • Table4 Deformation mode of thin plate under different explosion distance R when scaled distance Z=1.6m/kg1/3

  • Z减小至1.6m/kg1/3时,在R=375mm的情况下反直观模式替代了原有的继发直观模式;在R=625mm的情况下始发直观模式替代了原有的反直观模式.根据表3的规律可以预测在Z=1.6m/kg1/3的情况下R=425mm-575mm区间内必然存在一个R的范围使薄板发生混合变形模式.通过验证得到R=500mm时薄板发生混合变形模式.薄板变形模式总结如表4所示.

  • 表3和表4结果表明当距离当量比Z一定时,薄板在爆炸冲击载荷下存在顺序固定的四种变形模式.爆炸距离R较小时,首先发生继发直观模式,随着爆炸距离R的递增,依次发生反直观模式、混合模式、始发直观模式,而且这四种变形模式并不随着爆炸距离R等间隔变化.

  • 4.2 薄板响应模式机理分析

  • 薄板不同变形模式塑性耗散能的演化过程有所不同.继发直观模式下塑性耗散能发生三次抬升,与薄板的初次正向变形、回弹及再次正向变形相对应,如图12(a)所示;反直观与混合模式下塑性耗散能在薄板回弹过程同样存在抬升,如图12(b)(c)所示;直观模式下塑性耗散能在初始阶段快速上升后保持平稳,如图12(d)所示.各图中A阶段塑性耗散能的初始快速抬升阶段对应薄板初次正向变形.B阶段与C阶段均发生极值点反向屈曲,其中B阶段塑性耗散能缓慢上升对应薄板的回弹阶段,C阶段塑性耗散能进一步缓慢上升对应薄板的再次正向变形阶段.A阶段薄板初次正向变形产生的塑性耗散能在整个变形过程中占比最多,后续B阶段与C阶段极值点反向屈曲过程占比较少.

  • 当较低冲量造成的薄板初次正向变形较小,无法将冲击波提供能量完全耗散时,薄板会在负压的激励下发生回弹并发生继发直观模式,回弹与再次正向变形过程塑性耗散进一步增加,导致薄板在经历多次变形后最终停留在正方向.随着冲击波冲量的增加,薄板初次正向变形增加,塑性耗散增加,薄板仅需初次正向变形和回弹两个阶段即停止变形并最终停留在负方向,表现出反直观模式.同理,在冲量进一步增加后,薄板初次正向变形的塑性功耗散足以使薄板停止变形,即发生始发直观模式.而引起混合模式的冲击波冲量介于反直观模式与始发直观模式之间,该冲量范围内只能使薄板边缘区域发生回弹.

  • 综上,爆炸冲击载荷下薄板发生何种变形模式取决于初次正向变形塑性功耗散多少.初次塑性耗散能越少则极值点反向屈曲越多;初次塑性耗散能越多则极值点反向屈曲越少.

  • 图12 薄板动力行为模式与塑性耗散能关系

  • Fig.12 Relationship between dynamic behavior mode of thin plate and plastic dissipated energy

  • 5 结论

  • 本文讨论内容与爆炸载荷负压关系密切,爆炸载荷距离当量比较大时负压对结果有着显著影响,因此相关结论的适用范围也是在距离当量比较大的情况下.

  • (1)尽管爆炸冲击载荷的负压压强相对较小,但其对于薄板反直观异常行为的发生至关重要.

  • (2)爆炸冲击载荷负压段的描述方式对薄板反直观异常动力行为的响应历程有着显著影响.

  • (3)爆炸冲击载荷正压衰减过程影响薄板反直观行为的最大挠度与回弹过程持续时间.

  • (4)在距离当量比一定时,随着爆炸距离的增加,薄板依次发生继发直观、反直观、混合、始发直观四种行为模式.这种现象发生的原因是薄板不同正向变形程度的塑性耗散强弱不同,变形程度越小、塑性耗散越少,则越容易发生极值点反向屈曲.

  • 参考文献

    • [1] Xue Z,Hutchinson J W.Preliminary assessment of sandwich plates subjected to blast.International Journal of Mechanical Sciences,2003;45:687~705

    • [2] Aune V,Fagerholt E,Langseth M,et al.A shock tube facility to generate blast loading on structures.International Journal of Protective Structures,2016,7(3):340~66

    • [3] Krauthammer T,Altenberg A.Negative phase blast effects on glass panels.International Journal of Impact Engineering,2000,24(1):1~17

    • [4] Aune V,Fagerholt E,Hauge K O,et al.Experimental study on the response of thin aluminium and steel plates subjected to airblast loading.International Journal of Impact Engineering,2016,90:106~121

    • [5] Galiev S U.Experimental observations and discussion of counterintuitive behavior of plates and shallow shells subjected to blast loading.International Journal of Impact Engineering,1996,18(7-8):783~802

    • [6] Symonds P S,Yu T X.Counterintuitive behavior in a problem of elastic-plastic beam dynamics.Journal of Applied Mechanics,1985,52(3):517~522

    • [7] Yu T X.Chapter 9:Elastic effects in the dynamic plastic response of structures.Proceedings of the Third International Symposium on Structural Crashworthiness,1993,295~332

    • [8] Li Q M,Zhao L M,Yang G T.Experimental results on the counter-intuitive behaviour of thin clamped beams subjected to projectile impact.International Journal of Impact Engineering,1991,11(3):341~348

    • [9] 吴桂英,李鑫,周稳.空中爆炸载荷作用下铝合金圆板的反直观响应研究.振动与冲击,2010,29(12):166~168,246(Wu G V;Li X;Zhou W.The anomalous response of the aluminum alloy circular plate subjected to the air blast load.Journal of Vibration and Shock,2010,29(12):166~168,246(in Chinese))

    • [10] 吴桂英,秦冬祺,杨桂通.在冲击载荷作用下弹塑性圆板的反直观动力行为数值分析.固体力学学报,2004(4):389~393(Wu G Y;Qin D Q;Yang G T.The anomalous response of elastoplastic circular plates under impulsive load.Acta Mechanica Solida Sinica,2004(4):389~393(in Chinese))

    • [11] Kinney G F,Graham K J.Explosive shocks in air.Berlin:Springer,1985

    • [12] Friedlander F G.The diffraction of sound pulses I.Diffraction by a semi-infinite plane.Proceedings of the Royal society A:Mathematical,Physical and Engineering Sciences,1946,186(1006):322~44

    • [13] Swisdak J M.Simplified kingery airblast calculations.Minutes of the Twenty Sixth DOD Explosives Safety Seminar,DTIC Document,Maryland.

    • [14] Nassr A A.experimental and analytical study of the dynamic response of steel beams and columns to blast loading.McMaster-Civil Engineering,2012

    • [15] Rigby S E,Tyas A,Bennett T,et al.The negative phase of the blast load.International Journal of Protective Structures,2014,5(1):1~19

    • [16] Karlos V,Solomos G,Larcher M.Analysis of the blast wave decay coefficient using the Kingery-Bulmash data.International Journal of Protective Structures,2016,7(3):409~429

    • [17] Voce E.The relationship between stress and strain for homogenous deformation.Journal of the Japan Institute of Metals and Materials,1948,74:537~562

  • 参考文献

    • [1] Xue Z,Hutchinson J W.Preliminary assessment of sandwich plates subjected to blast.International Journal of Mechanical Sciences,2003;45:687~705

    • [2] Aune V,Fagerholt E,Langseth M,et al.A shock tube facility to generate blast loading on structures.International Journal of Protective Structures,2016,7(3):340~66

    • [3] Krauthammer T,Altenberg A.Negative phase blast effects on glass panels.International Journal of Impact Engineering,2000,24(1):1~17

    • [4] Aune V,Fagerholt E,Hauge K O,et al.Experimental study on the response of thin aluminium and steel plates subjected to airblast loading.International Journal of Impact Engineering,2016,90:106~121

    • [5] Galiev S U.Experimental observations and discussion of counterintuitive behavior of plates and shallow shells subjected to blast loading.International Journal of Impact Engineering,1996,18(7-8):783~802

    • [6] Symonds P S,Yu T X.Counterintuitive behavior in a problem of elastic-plastic beam dynamics.Journal of Applied Mechanics,1985,52(3):517~522

    • [7] Yu T X.Chapter 9:Elastic effects in the dynamic plastic response of structures.Proceedings of the Third International Symposium on Structural Crashworthiness,1993,295~332

    • [8] Li Q M,Zhao L M,Yang G T.Experimental results on the counter-intuitive behaviour of thin clamped beams subjected to projectile impact.International Journal of Impact Engineering,1991,11(3):341~348

    • [9] 吴桂英,李鑫,周稳.空中爆炸载荷作用下铝合金圆板的反直观响应研究.振动与冲击,2010,29(12):166~168,246(Wu G V;Li X;Zhou W.The anomalous response of the aluminum alloy circular plate subjected to the air blast load.Journal of Vibration and Shock,2010,29(12):166~168,246(in Chinese))

    • [10] 吴桂英,秦冬祺,杨桂通.在冲击载荷作用下弹塑性圆板的反直观动力行为数值分析.固体力学学报,2004(4):389~393(Wu G Y;Qin D Q;Yang G T.The anomalous response of elastoplastic circular plates under impulsive load.Acta Mechanica Solida Sinica,2004(4):389~393(in Chinese))

    • [11] Kinney G F,Graham K J.Explosive shocks in air.Berlin:Springer,1985

    • [12] Friedlander F G.The diffraction of sound pulses I.Diffraction by a semi-infinite plane.Proceedings of the Royal society A:Mathematical,Physical and Engineering Sciences,1946,186(1006):322~44

    • [13] Swisdak J M.Simplified kingery airblast calculations.Minutes of the Twenty Sixth DOD Explosives Safety Seminar,DTIC Document,Maryland.

    • [14] Nassr A A.experimental and analytical study of the dynamic response of steel beams and columns to blast loading.McMaster-Civil Engineering,2012

    • [15] Rigby S E,Tyas A,Bennett T,et al.The negative phase of the blast load.International Journal of Protective Structures,2014,5(1):1~19

    • [16] Karlos V,Solomos G,Larcher M.Analysis of the blast wave decay coefficient using the Kingery-Bulmash data.International Journal of Protective Structures,2016,7(3):409~429

    • [17] Voce E.The relationship between stress and strain for homogenous deformation.Journal of the Japan Institute of Metals and Materials,1948,74:537~562

  • 微信公众号二维码

    手机版网站二维码