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

绝热指数γ对平面爆轰过程中不同复合波区的参数特性影响分析

张娅 李晓彬 彭帅 施锐

张娅, 李晓彬, 彭帅, 施锐. 绝热指数γ对平面爆轰过程中不同复合波区的参数特性影响分析[J]. 爆炸与冲击, 2021, 41(10): 102102. doi: 10.11883/bzycj-2020-0458
引用本文: 张娅, 李晓彬, 彭帅, 施锐. 绝热指数γ对平面爆轰过程中不同复合波区的参数特性影响分析[J]. 爆炸与冲击, 2021, 41(10): 102102. doi: 10.11883/bzycj-2020-0458
ZHANG Xin-hua, TANG Zhi-ping, XU Wei-wei, TANG Xiao-jun, ZHENG Hang. Experimental study on characteristics of shock-induced phase transition and spallation in FeMnNi alloy[J]. Explosion And Shock Waves, 2007, 27(2): 103-108. doi: 10.11883/1001-1455(2007)02-0103-06
Citation: ZHANG Ya, LI Xiaobin, PENG Shuai, SHI Rui. Influences of adiabatic index γ on the parameters of different complex wave zones in a planar detonation[J]. Explosion And Shock Waves, 2021, 41(10): 102102. doi: 10.11883/bzycj-2020-0458

绝热指数γ对平面爆轰过程中不同复合波区的参数特性影响分析

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

    张 娅(1989- ),女,博士研究生,clairzya@126.com

    通讯作者:

    李晓彬(1971- ),男,博士,教授,lxbmark@163.com

  • 中图分类号: O381

Influences of adiabatic index γ on the parameters of different complex wave zones in a planar detonation

  • 摘要: 爆炸气体产物冲击膨胀过程中会形成多种复合波区,当爆炸气体绝热指数γ不同时其波区衰减特性差异较大。为研究不同γ条件下(γ>3,γ=3,γ<3)复合波区的特性差异,基于特征线法,对一平面爆轰过程中不同复合波区的波系相交特性进行了规律分析,并利用MATLAB对该平面爆轰过程进行流场模拟,验证并分析了不同复合波区流场内的参数变化特性。对比发现,γ不同时复合波区衰减特性的差异主要体现在与质点速度和气体声速相关的u-c平面特性上,其中在两中心稀疏波相交的复合波区,其差异还体现当γ≠3时相交的中心稀疏波不再具有中心发散特性。对爆炸过程中各波区特性的分析可为全面了解各特征参数的衰减规律提供参考。
  • 恐怖爆炸袭击的对象一般是办公楼、酒店、商场等人员相对密集集中的建筑,这些建筑70%由砌体结构构成,且多未考虑抗爆设计。当爆炸发生时,往往会造成砌体构件破坏或结构倒塌,并伴随大量高能碎片[1-2],从而对人员产生二次杀伤。因此,如何提高砌体结构的抗爆性能一直是研究中的热点。

    聚脲在低应变率下表现为高延性和超弹性,在高应变率下具有率敏感性,理论上具有较好的抗爆性能,目前已在防护领域中得到应用[3-5]。Davidion等[6]通过多组全尺寸爆炸试验,发现聚脲加固可以显著提高混凝土砌块墙的抗爆能力,同时延性较高的聚脲比其他高刚度复合材料可以更有效地阻止破片飞散和墙体倒塌。蔡桂杰[7]通过砖墙抗爆试验发现,背爆面喷涂聚脲的砖墙抗爆性能优于在迎爆面喷涂,指出喷涂聚脲厚度在3~9 mm范围内时,砖墙抗爆性能随喷涂厚度增加而增强。Lqbal等[8]发现聚脲喷涂厚度为6 mm可以发挥最大抗爆效果,增加聚脲厚度可能导致聚脲与结构剥离失效,但剥离后的聚脲涂层仍可限制高速碎片的产生[9]。然而,上述试验研究中多存在冲击波在墙体试件边缘发生绕射的干扰,造成爆炸荷载作用机制不清。且已有研究主要通过定性的分析破坏形态,总结强动载作用下加固结构的位移极限,缺少针对砖墙抗爆加固设计的定量计算方法。

    砖墙构件的抗爆设计计算一般是基于等效单自由度法(equivalent single degree of freedom, ESDOF)。Biggs[10]、UFC 3-340-2中[11]和Urgessa等[12]分别对单向砌体墙和双向混凝土板的系统单自由度等效过程、转换系数和抗力函数给出了建议取值。基于上述方法计算的砖墙等构件的超压-冲量曲线(pressure-impulse, P-I),也常用于工程抗爆设计和评估[13-15]。但是,在基于ESDOF方法研究聚脲加固砖墙的抗爆性能时,通常会将聚脲和砖墙两者耦合简化为单一均质材料,并忽视了聚脲涂层面的内张力。Irshidat等[16]将聚脲与砖墙的抗力分别计算,提出了聚脲加固单向砖墙的抗力函数,但并不适用于双向砌体墙的抗爆加固计算。

    本文基于一种经过防绕流改进的大型爆炸试验,分别针对聚脲喷涂加固及未加固的原型黏土砖框架填充墙试件开展爆炸加载试验,量测冲击波超压、冲量和砖墙挠度等时程数据,分析不同当量TNT爆炸作用下两种墙体试件的变形响应特征,归纳聚脲加固砖墙的爆炸失效模式和破坏特征;分别考虑砖墙及聚脲对墙体抗力的贡献,构建聚脲加固双向砌体墙的ESDOF模型和求解方法。

    爆炸试验过程中,空气冲击波在墙体边缘形成的稀疏波以及墙后绕射会影响墙体的真实抗爆动力响应。为此,本文专门改进设计了如图1所示的防绕流爆炸试验装置,试验前在特制的钢筋混凝土(reinforced concrete, RC)框架中预制砖墙,通过螺栓将养护后的砖墙框架固定在图1所示的爆炸试验装置上。爆炸试验装置通过RC腔室将边缘稀疏波以及墙后绕射波隔绝,可以测得仅正面冲击波作用下墙体的真实响应。通过在填充墙体试件的背爆面喷射聚脲涂层进行抗爆加固,喷涂完成后通过四周角钢把聚脲层边界固定在RC框架上,以防出现边界剪切破坏,加固后的砖墙如图1所示。

    图  1  试验装置
    Figure  1.  Test device

    框架填充墙墙体试件选用黏土空心砖,单个黏土空心砖的几何尺寸为24 cm×10 cm×4.5 cm,砖块间砌筑的砂浆厚度为1 cm,具体力学参数见表1,砌筑方式为全顺。试件墙面净跨为2 m×3 m,共制作了4面试验墙体,其中2面墙为未加固墙体试件,2面墙为背爆面喷涂聚脲加固试件。试验加固用聚脲抗爆涂层材料的聚脲断裂应变为2.5~3.4,真实失效应变约为1.2。加固施工过程中采用高压喷涂技术施工,在砖墙背爆面一次喷涂成型,不产生流挂现象,涂层厚度为0.6 cm。聚脲抗爆涂层材料的基本力学参数见表2

    表  1  砖和砂浆的力学参数
    Table  1.  Material parameters of brick and mortar
    材料密度/(kg·m−3)杨氏模量/MPa泊松比抗拉强度/MPa剪切强度/MPa屈服强度/MPa
    12008970.155.55.514.1
    砂浆21009130.253.53.57.03
    下载: 导出CSV 
    | 显示表格
    表  2  聚脲的力学参数
    Table  2.  Polyurea material parameters
    密度/(kg·m−3)杨氏模量/MPa泊松比抗拉强度/MPa屈服强度/MPa切线模量/MPa真实失效应变
    1150800.17155.514.11.2
    下载: 导出CSV 
    | 显示表格

    图2为试验装置及测量方案。试验前在RC框架中预制砖墙,通过螺栓将养护后的砖墙框架固定在试验装置上。试验中主要量测墙体上的爆炸荷载超压和墙体位移时程数据,超压和位移测点布置如图2所示。其中,墙体迎爆面空气冲击波压力测点有3个,为了防止墙体变形过程中传感器遭到破坏,压力传感器(PCB 102B04)布置在与墙体试件相同距离、同一水平高度的特制钢板上,测量范围为0~6.8 MPa,对应于墙体的超压测点(P1P2P3)位置如图2所示。墙体构件位移测点有3个(D1D2D3),使用接触式LVDT位移传感器,传感器垂直安装在防绕流爆炸试验装置内,传感器端头用角钢固定在背爆面上,位移传感器量程为±200 mm,位移传感器布置如图2所示;试验数据采集由DH8302N高速数采完成,采样频率为0.5 MHz。

    图  2  测试方案
    Figure  2.  Test scheme

    为了保证爆炸超压测试数据与墙面实际超压一致,降低炸药形状的影响,试验采用铸药的方式制作了圆柱形TNT药柱,药柱轴线方向与墙面平行,通过空中爆炸对砖墙施加冲击波荷载。炸药中心距地面1.5 m,距墙面3.0 m,共计开展了4个炮次的爆炸试验,炸药当量和试验工况见表3

    表  3  试验工况
    Table  3.  Test conditions
    试验聚脲厚度/mm比例爆距/(m·kg−1/3)装药当量/kg
    101.89 4
    201.3910
    361.89 4
    461.3910
    下载: 导出CSV 
    | 显示表格
    1.4.1   失效破坏特征

    图3为不同爆炸工况下墙体试件的损伤破坏情况。其中,4 kg TNT爆炸作用下试验1墙体中部出现轻微发育裂隙,墙体发生轻微的X状裂缝;相同荷载条件下,加固墙体(试验3)并未产生裂隙。随着TNT装药增加达到10 kg,砖墙受到的冲击波荷载明显增大,试验2墙体发生局部剪切破坏和坍塌,如图3(b)所示,试件2中部出现了尺寸约为1.0 m×1.5 m的局部坍塌;10 kg TNT爆炸荷载作用下加固砖墙破坏如图3(d)所示,可以发现,虽然在墙体边界和中部均出现了多条较大的裂缝,但并未出现倒塌,墙体整体呈现典型的“X”状双向板弯曲破坏形态。分析结果表明,聚脲加固不但可以增加结构抗力,而且可以使砖墙的破坏形态由局部剪切的脆性破坏转变为延性弯曲破坏。

    图  3  墙体试件破坏模式
    Figure  3.  Failure modes of wall specimens

    图4给出了爆炸后试验3、4中聚脲加固砖墙背爆面的聚脲变形情况。可以发现,试件2背爆面聚脲涂层光滑、平整、无裂缝,无明显变化。当荷载增加至10 kg时,背爆面聚脲虽未发生破坏,但是也出现了褶皱现象。对比图3(d),发现图4(b)的褶皱位置与正面裂缝对应,裂缝与褶皱分布位置与图3(b)中为加固墙体剪切破洞位置接近,进一步说明了背爆面喷涂聚脲可以有效限制爆炸作用下砖墙的局部剪切脆性破坏,大幅降低砖墙的倒塌风险。

    图  4  加固砖墙爆炸后背爆面涂层情况
    Figure  4.  Coatings of reinforcement brick wall after blasting
    1.4.2   爆炸荷载

    图5(a)为试验测得的爆炸工况3、4下墙面各测点的典型反射冲击波超压时程曲线,由此可计算出图5(b)所示的冲量曲线。对比4 kg TNT爆炸作用下测点P1的反射超压峰值为1.01 MPa,较测点P2的反射超压峰值增加了10%。10 kg TNT爆炸作用下P1测得的反射超压峰值为4.61 MPa,较测点P2的反射超压峰值2.14 MPa增加约35%。与4 kg TNT爆炸试验相比,10 kg TNT爆炸试验中测点P1的超压峰值提升了356%,测点P2的超压峰值提升了137%。由于T4试验药柱长度更长,造成试验4中测点P1P2的压力差值较大。整体来看,P1点的冲量计算结果与相同比例距离下依据TM-5-855-1规范[17]得到的计算结果对比吻合度较高,如表4所示。10 kg TNT爆炸作用下测点P1冲量为824 Pa·s,相比4 kg TNT作用下P1受到的冲量提升了269%;测点P2冲量为335 Pa·s,相比4 kg TNT作用下P2受到的冲量提升了123%。可以发现,随着装药增加砖墙受到爆炸超压和冲量都显著增大,试件更易发生剪切破坏。

    表  4  爆炸荷载验证
    Table  4.  Explosion load verification
    装药/kg冲量/(Pa·s)误差/%
    测点P1文献[17]方法
    430640023.5
    10824790 4.3
    下载: 导出CSV 
    | 显示表格
    图  5  试验荷载时程曲线
    Figure  5.  Curves of test load
    1.4.3   位移响应

    根据规范UFC-3-340-02[11] 中砌体墙破坏准则,支座转角被用来作为分类损伤等级的指标,如表5所示。对于双向墙而言,爆炸作用下砖墙支座转角的可修复限值为0.5°,不可修复限值为2°,经过反算,两种限值对应试验中墙体试件的跨中点挠度分别为8.72和34.9 mm。

    表  5  砌体墙的破坏准则[11]
    Table  5.  Failure criteria of masonry walls [11]
    破坏等级边界条件支座转角/(°)跨中允许挠度/mm
    可修复单向0.58.72
    双向0.58.72
    不可修复单向117.45
    双向234.90
    下载: 导出CSV 
    | 显示表格

    图6为试验中测得的不同装药设计条件下聚脲加固(polyurea reinforcement, PR)和未加固墙体D1测点的位移时程曲线。如图6(a)所示,试验1中未加固砖墙D1测点的最大正向位移达到21 mm,反向位移为12.5 mm,损伤等级为可修复损伤;而砖墙加固后,在试验3中D1测点的最大位移仅为16 mm,反向位移为18.5 mm。两者损伤等级虽同样为可修复损伤,但聚脲加固后试件在相同荷载作用下最大位移降低了23.8%,砖墙抗爆性能显然得到了大幅提升。

    图  6  不同荷载作用下两种墙体的位移响应
    Figure  6.  Displacement responses of two walls under different loads

    图6(b)为10 kg TNT爆炸荷载作用下加固和未加固墙体的最大位移均明显超过规范的允许值。其中,试验2中,未加固墙体试件的跨中点D1位移达到150 mm后,墙体中心部位发生局部剪切破坏而倒塌。而试验4中,加固后墙体试件D1测点最大位移达到了193 mm,填充墙体与框架边界及墙面均出现较大裂缝,但墙体并未发生倒塌破坏。聚脲加固后墙体最大延性至少可提升28.7%。对比可以发现,喷涂聚脲可显著增加砖墙加固墙体抗爆能力,大幅增大构件延性,并阻止墙体出现剪切局部破坏。

    1.4.4   单面聚脲加固作用机理

    图6(a)中4 kg TNT爆炸作用下未加固框架填充墙的跨中点最大正向位移比反向最大位移大68%,而加固砖墙的正向最大位移反而比反向最大位移小15.6%。图7(a)给出了进一步提取的4 kg TNT爆炸作用下加固砖墙D1测点的位移时程曲线。抓取测点的平衡位置到振动峰值位置的运动时间,得到4个时间段内砖墙的平均自振频率(ω)。计算发现加固砖墙在4个时间段的平均振动频率并非定值,在一个振动周期内平均振动频率逐渐下降,并呈周期性变化,如图7(b)所示。假定构件质量不发生变化,由ω=K/M可以反推得到加固砖墙刚度K[18]也呈周期性变化,加固砖墙正向振动过程中加固砖墙刚度较反向振动增加了133%。说明加固墙体的振动过程中,背爆面喷涂的聚脲涂层对加固墙体正反向振动变形刚度的贡献不一样,对于正向变形的刚度增加贡献较大,而对于反向反弹变形的刚度贡献较小。

    图  7  背爆面加固作用机理
    Figure  7.  Mechanism of back burst surface strengthening

    图8给出了墙体受荷振动变形过程中,跨中点D1附近局部砖砌体的受力情况。墙体正向运动过程中墙体发生弯曲变形,背爆面局部砖块间出现如图8(b)所示的拉伸变形,此时聚脲涂层则在砖块间拉伸变形方向上提供反向拉力T,且作用力位置与砖块间拉伸变形位置重合。此时,聚脲拉力T和砂浆与砖块之间粘接力fvm可共同为加固砖墙提供抗力。墙体反向运动时,墙体迎爆面局部砖块间出现如图8(a)的拉伸变形,此时背爆面聚脲变形较小,因而提供的拉力T较小,且拉力作用点位置相对于砖块间拉伸变形位置偏移了一个墙厚的距离,对加固砖墙抵抗变形的抗力贡献较小,因此,相同变形情况下,背爆面聚脲加固砖墙正向弯曲抗力较负向弯曲抗力更大。如图3(a)所示,4 kg TNT爆炸作用下未加固砖墙出现明显裂缝,这是因为墙体构件材料局部损伤导致砖墙变形过程中能量耗散,从而造成反向位移峰值的下降。相同荷载作用下加固砖墙并未出现明显裂缝,正向最大位移未达到损伤标准,能量耗散较少,且墙体反向变形的抗力较正向变形抗力低,导致墙体试件反向最大位移增加。此时若聚脲与墙体表面的黏结力不足,有可能导致墙体发生反向反弹破坏。

    图  8  不同振动位置砌块间受力
    Figure  8.  Force between blocks at different vibration positions

    出于安全考虑,试验中各工况设计的爆炸比例爆距较为适中,试验工况中聚脲涂层并未出现明显的失效破坏。为进一步研究聚脲加固砖墙的失效破坏机理,选择文献[19]的相关数据进行了比较分析,该文献报告了近距离爆炸作用下聚脲加固砖墙的一系列试验结果,其试验设计如表6所示。

    表  6  本文与文献试验工况对比
    Table  6.  Comparison of test conditions between this paper and the literature
    工况聚脲厚度/mm装药/kg比例爆距/(m·kg−1/3)
    1[19]0, 350.584
    2[19]6, 650.35
    本文试验30, 641.89
    本文试验40, 610 1.39
    下载: 导出CSV 
    | 显示表格

    结合文献[19]中报告的爆炸试验结果和本文爆炸试验结果可以发现,随着试验比例爆距和喷涂聚脲厚度d不断降低,砖墙背爆面聚脲涂层出现撕裂、鼓包和剥离。当比例爆距达到0.584 m/kg1/3时背爆面出现裂缝;比例爆距为0.35 m/kg1/3时,背爆面出现局部震塌破坏,该位置附近的聚脲涂层出现明显鼓包和剥离,而聚脲涂层有效的限制了背爆面砖墙碎片的飞散,如图9所示(图中Z为比例爆距)。

    图  9  试验现象[19]
    Figure  9.  Test phenomena [19]

    可以认为,砖墙背爆面聚脲涂层的失效破坏特征主要与聚脲厚度及比例爆距相关,爆炸作用下聚脲涂层失效主要为局部撕裂或鼓包。因此,如图10所示,可根据比例爆距Z和涂层厚度d将聚脲加固砖墙的失效模式划分为以下3种:

    图  10  动载响应下聚脲涂层失效机理
    Figure  10.  Failure mechanism of polyurea coating under dynamic load response

    (1)当Z≥1.89 m/kg1/3时,墙体呈现典型的双向板弯曲破坏形态;如图10(b)所示,聚脲涂层会随着局部砖块与砖块间的弯曲开裂而发生拉伸变形,通常极限拉伸应变会超过砖块高度的1.5倍[6]

    (2) 当Z≤0.584 m/kg1/3d≤3 mm时,砖块间的剪切变形超过了聚脲的极限剪切失效应变,聚脲材料本身发生断裂,此时结构会出现局部剪切破坏,如图10(c)所示;

    (3) 当Z≤0.35 m/kg1/3d≥6 mm时,爆炸冲击波的波动效应导致砖墙背爆面震塌,震塌位置的聚脲与墙体发生剥离,出现局部鼓包破坏,如图10(d)所示。

    由于试验中聚脲加固砖墙在爆炸荷载作用下主要呈弯曲破坏,因此可以采用等效单自由度模型(ESDOF)来进行简化计算。由达朗贝尔原理建立聚脲加固砖墙体系的等效动力微分方程为:

    KMLm¨z(t)+R(z)=F(t) (1)

    式中:KML=KM/KL为弹性阶段等效质量荷载系数,KM为质量系数,KL为荷载系数,z为垂直于墙面方向的动位移,本文中即为挠度;R(z)为构件抗力函数;F(t)为瞬时动态作用力函数。根据运动方程和能量方程,可以得到等效体系的抗力系数KR、荷载系数KL和质量系数KM

    KR=KL=1AAψ(x,y)dxdy (2)
    KM=1AAψ2(x,y)dxdy (3)

    式中:ψ(x,y)为双向墙体的振形函数,A为墙体面积;根据Biggs的建议[10],弹性阶段KL=0.64,KM=0.5;塑性阶段构件KL=0.5,KM=0.33。

    图11为在爆炸荷载作用下双向墙体结构弯曲变形示意图,图中,ab分别为结构长短边长度,S1S2为挠曲线长度,θaθb为支座转角,z为挠度,假定加固材料形成的膜与结构之间无粘结[8]

    图  11  墙体结构弯曲变形示意
    Figure  11.  Deformation schematic of wall structure

    此时砖墙结构抗力函数可以表示为:

    R1(z)={KezzzeRuzezzm (4)

    式中:R1(z)为未加固砖墙的抗力函数;Ke为初始刚度[10],由UFC-3-340-02规范[17]计算得到1.77×107 N/m;zezm分别为弹、塑性极限挠度,支座的允许转角θaθb为2°,Ru为塑性阶段的抗力限值,可由公式Ru=5(Mfa+Msa)ab/z2计算得到,其中,MfaMsa为固定边界及跨中处单位长度弯矩抗力。

    对于聚脲涂层,其抗弯刚度很小,聚脲涂层可简化与墙体平行的为面力,其应力-应变关系满足胡克定律,砖墙弯曲时聚脲沿长边和短边方向的应力为:

    σ1=E(S1b)b (5)
    σ2=E(S2a)a (6)

    式中:σ1σ2为长边和短边方向的应力;E为聚脲材料的弹性模量;S1S2为挠曲线长度,ab为长短边长度。得到聚脲的抗力函数Rv为:

    Rv=Ta sinθ1+Tb sinθ2=naσa sinθa+nbσb sinθb (7)

    式中:TaTb为聚脲涂层短边和长边方向拉应力;n为应变强化指数,n取4[4]

    加固砌体墙结构的抗力函数为:

    R(z)={R1(z)+RvzzmRvzmz (8)

    需要注意的是,基于1.3节的研究可以发现,爆炸作用下加固砖墙振动过程中等效刚度是变化的,因此抗力模型仅适用于正向弯曲条件下加固砖墙体系变形求解。

    将迎爆面的爆炸荷载等效为均布载荷[21]

    F=12F1+F2+76F323F4 (9)

    式中:F为等效均布荷载峰值,F1为结构表面中心处反射超压峰值,F2F3为长边、短边中点反射超压峰值,F4为边角处反射超压峰值。试验中测得的4 kg TNT等效均布荷载峰值为477 kPa。

    荷载作用时间Δt为:

    Δt=2I1F (10)

    式中:I1为墙中心爆炸波正反射冲量。振型函数为:

    ψ(x,y)=sinπxasinπyb (11)

    将爆炸空气冲击波荷载简化为突加线性衰减的冲击波,将此种形式的荷载作为动载代入式(1),本文采用预估-校正形式的显式Newmark方法[22],求解得到等效单自由度体系的位移时程解。

    以试验中聚脲加固砖墙为计算对象,采用数值模拟的方法对ESDOF计算结果进行对比验证。有限元计算模型中,砌块与砂浆接触方式采用共节点,砌体墙与框架、聚脲的接触为固连失效接触。砖和砂浆的压缩和拉伸动力放大因子为1.2[23-24],砌块的失效应变为0.02。文献[25]已经验证了上述数值计算模型的可靠性。

    图12为ESDOF计算结果与数值模拟和试验结果对比。其中,4 kg TNT爆炸作用下加固墙体位移响应的三种计算结果如图12(a)所示,ESDOF模型计算得到的D1正向最大位移为16.7 mm,数值模拟中D1最大位移为13.2 mm,与试验相比两者分别增加5%和−17.5%。说明,有限元分析结果与试验结果相比,其接触算法高估了砌块与砂浆的接触力,体系具有更大的抗力;ESDOF方法基于实际工况分别考虑了砖墙及聚脲对墙体抗力的贡献,其计算结果与试验结果相比具有更好的一致性。

    图  12  计算结果对比
    Figure  12.  Comparison of calculation results

    10 kg TNT爆炸作用计算过程中发现,砖墙边界出现及墙面出现如图7(d)的剪切滑移。ESDOF模型则是采用固定边界的假定,故ESDOF模型在D1位置理论解比试验低了10.5%。而数值模拟采用了删除单元算法,所以构件质量下降刚度降低,导致模拟中D1位置的最大位移比试验大12.6%。总体而言,ESDOF模型与数值模拟相比可以更加准确地计算出砖墙的正向位移响应过程。

    基于一种改进的大型爆炸试验装置,开展了聚脲加固框架砖填充墙的原型爆炸试验,分析了爆炸荷载作用下加固砖墙的动力响应特征和破坏过程及模式,揭示了其失效破坏机理,建立了聚脲加固双向墙抗爆性能的理论计算模型,得到以下主要结论:

    (1) 喷涂聚脲抗爆涂层可以显著增加框架填充墙的抗爆性能,增加墙体的延性,且可将墙体的局部剪切脆性破坏转变为典型的弯曲延性破坏;

    (2) 墙体受爆炸荷载作用来回振动过程中,背爆面喷涂加固产生的聚脲拉力对砖块间的拉伸变形与墙体振动变形抗力的贡献机制不同,导致墙体抵抗正向变形抗力高,抵抗负向变形抗力低,不同时段墙体的最大平均变形刚度较最小平均变形刚度增加了133%;

    (3) 根据比例爆距不同确定了聚脲加固砖墙的三种不同失效破坏模式;当爆炸比例爆距大于等于1.89 m/kg1/3时,墙体呈现典型的双向板弯曲破坏形态;当爆炸比例爆距不超过0.584 m/kg1/3且聚脲厚度不超过3 mm时,墙体呈现局部剪切破坏;当比例爆距不超过0.35 m/kg1/3且聚脲厚度大于等于6 mm时,聚脲与墙体发生剥离,出现局部鼓包破坏;

    (4) 提出的改进ESDOF方法可以较为准确地预测爆炸作用下背爆面加固双向砖墙的正向位移响应过程,为相似工程设计提供理论计算参考和依据。

  • 图  1  平面爆轰过程中复合波区的形成及相对位置(+右行,−左行)

    Figure  1.  Formation and relative positions of complex wave zones in plane detonation (+ showing a going-right wave, − showing a going-left wave)

    图  2  两稀疏波在原点相遇时流场内的特征参数分布

    Figure  2.  Distribution of characteristic parameters when two rarefaction waves meet at the origin

    图  3  两中心稀疏波相交的x-t平面

    Figure  3.  x-t plane of intersection of two central rarefaction waves

    图  4  应用c值均布特性后系统流场特性变化

    Figure  4.  Changes of flow field characteristics after the application of c-value uniform distribution

    图  5  中心稀疏波的发散特性

    Figure  5.  Divergence characteristics of central rarefaction waves

    图  6  两同等强度稀疏波相交的平面特性

    Figure  6.  Plane characteristics of the intersection of two rarefaction waves with equal intensity

    图  7  两同等强度压缩波相交的平面特性

    Figure  7.  Plane characteristics of the intersection of two compression waves with equal intensity

    图  8  不同强度压缩波和稀疏波相交的平面特性

    Figure  8.  Plane characteristics of the intersection of a compression wave and a rarefaction wave with different intensities

    图  9  γ≠3时两同等强度稀疏波相交的平面特性示意图

    Figure  9.  Plane characteristics of intersection of two rarefaction waves with equal intensity at γ≠3

    图  10  对复合波区的迭代思路

    Figure  10.  The iterative method for the complex wave zone

    图  11  MATLAB计算生成的特征线流场

    Figure  11.  Characteristic line flow field generated by MATLAB calculation

    图  12  T2T3Tn时刻的流场参数分布

    Figure  12.  Distribution of flow field parameters at times T2, T3 and Tn

    图  13  爆轰流场中复合波⑤区的u-c特性

    Figure  13.  u-c characteristics of complex wave zone ⑤ in the post-detonation flow field

    图  14  爆轰流场中复合波⑦区的u-c特性

    Figure  14.  u-c characteristics of complex wave zone ⑦ in the post-detonation flow field

    表  1  阶段1各区域的状态参量

    Table  1.   Characteristic parameters of each region in stage 1

    区域压力p/
    MPa
    密度ρ/
    (kg·m−3
    质点速度u/
    (m·s−1
    气体声速c/
    (m·s−1
    绝热指数γ
    ③区9130.8 163003 620.5002.34
    ④区介于③区和②区之间,参数按中心稀疏波规律衰减2.34
    ②区 28.29 138.054 730.127 692.5102.34
    ①区 28.29 7.204 730.1272 345.5671.40
    ⓪区 0.101 1.2250 340.2941.40
    下载: 导出CSV
  • [1] 杨科之, 刘盛. 空气冲击波传播和衰减研究进展 [J]. 防护工程, 2020, 42(3): 1–10. DOI: 10.3969/j.issn.1674-1854.2020.03.001.

    YANG K Z, LIU S. Progress of research on propagation and attenuation of air blast [J]. Protective Engineering, 2020, 42(3): 1–10. DOI: 10.3969/j.issn.1674-1854.2020.03.001.
    [2] KAHALI S, TOWNSEND M, NGUYEN M M, et al. The evolution of secondary flow phenomena and their effect on primary shock conditions in shock tubes: Experimentation and numerical model [J]. PLoS ONE, 2020, 15(1): e0227125. DOI: 10.1371/journal.pone.0227125.
    [3] 徐春光, 白晓征, 刘瑜, 等. 爆炸激波管管口稀疏波对试验段的影响 [J]. 国防科技大学学报, 2011, 33(4): 1–5. DOI: 10.3969/j.issn.1001-2486.2011.04.001.

    XU C G, BAI X Z, LIU Y, et al. Research on the influence of rarefaction wave to the experimental section in blast shock tube [J]. Journal of National University of Defense Technology, 2011, 33(4): 1–5. DOI: 10.3969/j.issn.1001-2486.2011.04.001.
    [4] 北京工业学院八系《爆炸及其作用》编写组. 爆炸及其作用: 上册: 气体动力学基础和爆轰理论[M]. 北京: 国防工业出版社, 1979: 189−193.
    [5] 张守中. 爆炸与冲击动力学[M]. 北京: 兵器工业出版社, 1993: 120−126.
    [6] 赵铮, 陶钢, 杜长星. 爆轰产物JWL状态方程应用研究 [J]. 高压物理学报, 2009, 23(4): 277–282. DOI: 10.11858/gywlxb.2009.04.007.

    ZHAO Z, TAO G, DU Z X. Application research on JWL equation of state of detonation products [J]. Chinese Journal of High Pressure Physics, 2009, 23(4): 277–282. DOI: 10.11858/gywlxb.2009.04.007.
    [7] 王成, 徐文龙, 郭宇飞. 基于基因遗传算法和γ律状态方程的JWL状态方程参数计算 [J]. 兵工学报, 2017, 38(S1): 167–173.

    WANG C, XU W L, GUO X F. Calculation of JWL equation of state parameters based on genetic algorithm and γ equation of state [J]. Acta Armamentarii, 2017, 38(S1): 167–173.
    [8] TRZCIŃSKI W A, CUDZIO S. Characteristics of high explosives obtained from cylinder test data [J]. Chinese Journal of Energetic Materials, 2006, 14(1): 1–7.
    [9] 刘文祥, 李捷, 钟方平. 组合γ状态方程在数值计算上的应用及分析 [J]. 爆炸与冲击, 2009, 29(2): 209–212. DOI: 10.11883/1001-1455(2009)02-0209-04.

    LIU W X, LI J, ZHONG F P. Application and analysis of a combinational γ equation of state for detonation products [J]. Explosion and Shock Waves, 2009, 29(2): 209–212. DOI: 10.11883/1001-1455(2009)02-0209-04.
    [10] 徐维铮, 吴卫国. JWL状态方程及其等效多方状态方程在内爆炸计算中的应用分析 [J]. 中国舰船研究, 2019, 14(3): 83–91. DOI: 10.19693/j.issn.1673-3185.01328.

    XU W Z, WU W G. Application analysis of JWL EOS and the equivalent polytropic EOS in internal explosion calculation [J]. Chinese Journal of Ship Research, 2019, 14(3): 83–91. DOI: 10.19693/j.issn.1673-3185.01328.
    [11] 恽寿榕, 赵衡阳. 爆炸力学[M]. 北京: 国防工业出版社, 2005: 31−34.
    [12] 张连玉. 爆炸气体动力学基础[M]. 北京: 北京工业学院出版社, 1987: 220−222.
  • 加载中
图(14) / 表(1)
计量
  • 文章访问数:  547
  • HTML全文浏览量:  203
  • PDF下载量:  63
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-07
  • 修回日期:  2021-03-01
  • 网络出版日期:  2021-09-14
  • 刊出日期:  2021-10-13

目录

/

返回文章
返回