Influences of adiabatic index γ on the parameters of different complex wave zones in a planar detonation
-
摘要: 爆炸气体产物冲击膨胀过程中会形成多种复合波区,当爆炸气体绝热指数γ不同时其波区衰减特性差异较大。为研究不同γ条件下(γ>3,γ=3,γ<3)复合波区的特性差异,基于特征线法,对一平面爆轰过程中不同复合波区的波系相交特性进行了规律分析,并利用MATLAB对该平面爆轰过程进行流场模拟,验证并分析了不同复合波区流场内的参数变化特性。对比发现,γ不同时复合波区衰减特性的差异主要体现在与质点速度和气体声速相关的u-c平面特性上,其中在两中心稀疏波相交的复合波区,其差异还体现当γ≠3时相交的中心稀疏波不再具有中心发散特性。对爆炸过程中各波区特性的分析可为全面了解各特征参数的衰减规律提供参考。Abstract: A variety of complex wave zones are formed during the impact and expansion of explosive gas products. When the adiabatic index γ of explosive gas is different, the attenuation characteristics of the wave zones are quite different. In order to understand the characteristics of the complex wave zone under the different γ conditions (γ>3, γ=3, γ<3), the intersection characteristics of different complex wave zones in a planar detonation were analyzed based on the characteristic line method. The flow field in the planar detonation was simulated by MATLAB to verify and analyze the parameter change characteristics of the flow fields in the different complex wave zones. Comparisons display that the differences in the attenuation characteristics of the wave zones are mainly reflected in the differences in the u-c plane characteristics related to the particle velocity and the gas sound velocity. Among them, in the complex wave zone where two rarefaction waves intersect, the difference is also reflected in that when γ≠3, the rarefaction wave no longer has the characteristic of central divergence. The analysis result on the characteristics of each wave zone in the explosion process provides a reference for comprehensively understanding the attenuation of each characteristic parameter.
-
恐怖爆炸袭击的对象一般是办公楼、酒店、商场等人员相对密集集中的建筑,这些建筑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. 聚脲加固砖墙抗爆试验
1.1 试验装置及试件制作
爆炸试验过程中,空气冲击波在墙体边缘形成的稀疏波以及墙后绕射会影响墙体的真实抗爆动力响应。为此,本文专门改进设计了如图1所示的防绕流爆炸试验装置,试验前在特制的钢筋混凝土(reinforced concrete, RC)框架中预制砖墙,通过螺栓将养护后的砖墙框架固定在图1所示的爆炸试验装置上。爆炸试验装置通过RC腔室将边缘稀疏波以及墙后绕射波隔绝,可以测得仅正面冲击波作用下墙体的真实响应。通过在填充墙体试件的背爆面喷射聚脲涂层进行抗爆加固,喷涂完成后通过四周角钢把聚脲层边界固定在RC框架上,以防出现边界剪切破坏,加固后的砖墙如图1所示。
框架填充墙墙体试件选用黏土空心砖,单个黏土空心砖的几何尺寸为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 砖 1200 897 0.15 5.5 5.5 14.1 砂浆 2100 913 0.25 3.5 3.5 7.03 表 2 聚脲的力学参数Table 2. Polyurea material parameters密度/(kg·m−3) 杨氏模量/MPa 泊松比 抗拉强度/MPa 屈服强度/MPa 切线模量/MPa 真实失效应变 1150 80 0.17 15 5.5 14.1 1.2 1.2 试验装置与量测
图2为试验装置及测量方案。试验前在RC框架中预制砖墙,通过螺栓将养护后的砖墙框架固定在试验装置上。试验中主要量测墙体上的爆炸荷载超压和墙体位移时程数据,超压和位移测点布置如图2所示。其中,墙体迎爆面空气冲击波压力测点有3个,为了防止墙体变形过程中传感器遭到破坏,压力传感器(PCB 102B04)布置在与墙体试件相同距离、同一水平高度的特制钢板上,测量范围为0~6.8 MPa,对应于墙体的超压测点(P1、P2、P3)位置如图2所示。墙体构件位移测点有3个(D1、D2、D3),使用接触式LVDT位移传感器,传感器垂直安装在防绕流爆炸试验装置内,传感器端头用角钢固定在背爆面上,位移传感器量程为±200 mm,位移传感器布置如图2所示;试验数据采集由DH8302N高速数采完成,采样频率为0.5 MHz。
1.3 试验工况
为了保证爆炸超压测试数据与墙面实际超压一致,降低炸药形状的影响,试验采用铸药的方式制作了圆柱形TNT药柱,药柱轴线方向与墙面平行,通过空中爆炸对砖墙施加冲击波荷载。炸药中心距地面1.5 m,距墙面3.0 m,共计开展了4个炮次的爆炸试验,炸药当量和试验工况见表3。
表 3 试验工况Table 3. Test conditions试验 聚脲厚度/mm 比例爆距/(m·kg−1/3) 装药当量/kg 1 0 1.89 4 2 0 1.39 10 3 6 1.89 4 4 6 1.39 10 1.4 试验结果和分析
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”状双向板弯曲破坏形态。分析结果表明,聚脲加固不但可以增加结构抗力,而且可以使砖墙的破坏形态由局部剪切的脆性破坏转变为延性弯曲破坏。
图4给出了爆炸后试验3、4中聚脲加固砖墙背爆面的聚脲变形情况。可以发现,试件2背爆面聚脲涂层光滑、平整、无裂缝,无明显变化。当荷载增加至10 kg时,背爆面聚脲虽未发生破坏,但是也出现了褶皱现象。对比图3(d),发现图4(b)的褶皱位置与正面裂缝对应,裂缝与褶皱分布位置与图3(b)中为加固墙体剪切破洞位置接近,进一步说明了背爆面喷涂聚脲可以有效限制爆炸作用下砖墙的局部剪切脆性破坏,大幅降低砖墙的倒塌风险。
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中测点P1、P2的压力差值较大。整体来看,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]方法 4 306 400 23.5 10 824 790 4.3 1.4.3 位移响应
根据规范UFC-3-340-02[11] 中砌体墙破坏准则,支座转角被用来作为分类损伤等级的指标,如表5所示。对于双向墙而言,爆炸作用下砖墙支座转角的可修复限值为0.5°,不可修复限值为2°,经过反算,两种限值对应试验中墙体试件的跨中点挠度分别为8.72和34.9 mm。
图6为试验中测得的不同装药设计条件下聚脲加固(polyurea reinforcement, PR)和未加固墙体D1测点的位移时程曲线。如图6(a)所示,试验1中未加固砖墙D1测点的最大正向位移达到21 mm,反向位移为12.5 mm,损伤等级为可修复损伤;而砖墙加固后,在试验3中D1测点的最大位移仅为16 mm,反向位移为18.5 mm。两者损伤等级虽同样为可修复损伤,但聚脲加固后试件在相同荷载作用下最大位移降低了23.8%,砖墙抗爆性能显然得到了大幅提升。
图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%。说明加固墙体的振动过程中,背爆面喷涂的聚脲涂层对加固墙体正反向振动变形刚度的贡献不一样,对于正向变形的刚度增加贡献较大,而对于反向反弹变形的刚度贡献较小。图8给出了墙体受荷振动变形过程中,跨中点D1附近局部砖砌体的受力情况。墙体正向运动过程中墙体发生弯曲变形,背爆面局部砖块间出现如图8(b)所示的拉伸变形,此时聚脲涂层则在砖块间拉伸变形方向上提供反向拉力T,且作用力位置与砖块间拉伸变形位置重合。此时,聚脲拉力T和砂浆与砖块之间粘接力
fvm 可共同为加固砖墙提供抗力。墙体反向运动时,墙体迎爆面局部砖块间出现如图8(a)的拉伸变形,此时背爆面聚脲变形较小,因而提供的拉力T较小,且拉力作用点位置相对于砖块间拉伸变形位置偏移了一个墙厚的距离,对加固砖墙抵抗变形的抗力贡献较小,因此,相同变形情况下,背爆面聚脲加固砖墙正向弯曲抗力较负向弯曲抗力更大。如图3(a)所示,4 kg TNT爆炸作用下未加固砖墙出现明显裂缝,这是因为墙体构件材料局部损伤导致砖墙变形过程中能量耗散,从而造成反向位移峰值的下降。相同荷载作用下加固砖墙并未出现明显裂缝,正向最大位移未达到损伤标准,能量耗散较少,且墙体反向变形的抗力较正向变形抗力低,导致墙体试件反向最大位移增加。此时若聚脲与墙体表面的黏结力不足,有可能导致墙体发生反向反弹破坏。1.5 聚脲加固砖墙的失效破坏机理
出于安全考虑,试验中各工况设计的爆炸比例爆距较为适中,试验工况中聚脲涂层并未出现明显的失效破坏。为进一步研究聚脲加固砖墙的失效破坏机理,选择文献[19]的相关数据进行了比较分析,该文献报告了近距离爆炸作用下聚脲加固砖墙的一系列试验结果,其试验设计如表6所示。
结合文献[19]中报告的爆炸试验结果和本文爆炸试验结果可以发现,随着试验比例爆距和喷涂聚脲厚度d不断降低,砖墙背爆面聚脲涂层出现撕裂、鼓包和剥离。当比例爆距达到0.584 m/kg1/3时背爆面出现裂缝;比例爆距为0.35 m/kg1/3时,背爆面出现局部震塌破坏,该位置附近的聚脲涂层出现明显鼓包和剥离,而聚脲涂层有效的限制了背爆面砖墙碎片的飞散,如图9所示(图中Z为比例爆距)。
可以认为,砖墙背爆面聚脲涂层的失效破坏特征主要与聚脲厚度及比例爆距相关,爆炸作用下聚脲涂层失效主要为局部撕裂或鼓包。因此,如图10所示,可根据比例爆距Z和涂层厚度d将聚脲加固砖墙的失效模式划分为以下3种:
(1)当Z≥1.89 m/kg1/3时,墙体呈现典型的双向板弯曲破坏形态;如图10(b)所示,聚脲涂层会随着局部砖块与砖块间的弯曲开裂而发生拉伸变形,通常极限拉伸应变会超过砖块高度的1.5倍[6];
(2) 当Z≤0.584 m/kg1/3且d≤3 mm时,砖块间的剪切变形超过了聚脲的极限剪切失效应变,聚脲材料本身发生断裂,此时结构会出现局部剪切破坏,如图10(c)所示;
(3) 当Z≤0.35 m/kg1/3且d≥6 mm时,爆炸冲击波的波动效应导致砖墙背爆面震塌,震塌位置的聚脲与墙体发生剥离,出现局部鼓包破坏,如图10(d)所示。
2. 聚脲加固砖墙抗爆的简化计算方法
2.1 ESDOF分析方法
由于试验中聚脲加固砖墙在爆炸荷载作用下主要呈弯曲破坏,因此可以采用等效单自由度模型(ESDOF)来进行简化计算。由达朗贝尔原理建立聚脲加固砖墙体系的等效动力微分方程为:
KMLm¨z(t)+R(z)=F(t) (1) 式中:
KML=KM/KL 为弹性阶段等效质量荷载系数,KM 为质量系数,KL 为荷载系数,z为垂直于墙面方向的动位移,本文中即为挠度;R(z) 为构件抗力函数;F(t) 为瞬时动态作用力函数。根据运动方程和能量方程,可以得到等效体系的抗力系数KR、荷载系数KL和质量系数KM为KR=KL=1A∬Aψ(x,y)dxdy (2) KM=1A∬Aψ2(x,y)dxdy (3) 式中:
ψ(x,y) 为双向墙体的振形函数,A为墙体面积;根据Biggs的建议[10],弹性阶段KL=0.64,KM=0.5;塑性阶段构件KL=0.5,KM=0.33。2.2 聚脲加固双向填充墙抗力模型
图11为在爆炸荷载作用下双向墙体结构弯曲变形示意图,图中,a、b分别为结构长短边长度,S1、S2为挠曲线长度,θa、θb为支座转角,z为挠度,假定加固材料形成的膜与结构之间无粘结[8]。
此时砖墙结构抗力函数可以表示为:
R1(z)={Kezz<zeRuze<z<zm (4) 式中:
R1(z) 为未加固砖墙的抗力函数;Ke 为初始刚度[10],由UFC-3-340-02规范[17]计算得到1.77×107 N/m;ze 、zm 分别为弹、塑性极限挠度,支座的允许转角θa、θb为2°,Ru 为塑性阶段的抗力限值,可由公式Ru=5(Mfa+Msa)ab/z2 计算得到,其中,Mfa 和Msa 为固定边界及跨中处单位长度弯矩抗力。对于聚脲涂层,其抗弯刚度很小,聚脲涂层可简化与墙体平行的为面力,其应力-应变关系满足胡克定律,砖墙弯曲时聚脲沿长边和短边方向的应力为:
σ1=E(S1−b)b (5) σ2=E(S2−a)a (6) 式中:
σ1 和σ2 为长边和短边方向的应力;E 为聚脲材料的弹性模量;S1 和S2 为挠曲线长度,a、b 为长短边长度。得到聚脲的抗力函数Rv 为:Rv=Ta sinθ1+Tb sinθ2=naσa sinθa+nbσb sinθb (7) 式中:
Ta 、Tb 为聚脲涂层短边和长边方向拉应力;n为应变强化指数,n取4[4]。加固砌体墙结构的抗力函数为:
R(z)={R1(z)+Rvz<zmRvzm≤z (8) 需要注意的是,基于1.3节的研究可以发现,爆炸作用下加固砖墙振动过程中等效刚度是变化的,因此抗力模型仅适用于正向弯曲条件下加固砖墙体系变形求解。
3. 求解方法和模型验证
3.1 求解方法
将迎爆面的爆炸荷载等效为均布载荷[21]:
F=−12F1+F2+76F3−23F4 (9) 式中:F为等效均布荷载峰值,
F1 为结构表面中心处反射超压峰值,F2、F3 为长边、短边中点反射超压峰值,F4 为边角处反射超压峰值。试验中测得的4 kg TNT等效均布荷载峰值为477 kPa。荷载作用时间Δt为:
Δt=2I1F (10) 式中:
I1 为墙中心爆炸波正反射冲量。振型函数为:ψ(x,y)=sinπxasinπyb (11) 将爆炸空气冲击波荷载简化为突加线性衰减的冲击波,将此种形式的荷载作为动载代入式(1),本文采用预估-校正形式的显式Newmark方法[22],求解得到等效单自由度体系的位移时程解。
3.2 模型验证
以试验中聚脲加固砖墙为计算对象,采用数值模拟的方法对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方法基于实际工况分别考虑了砖墙及聚脲对墙体抗力的贡献,其计算结果与试验结果相比具有更好的一致性。
10 kg TNT爆炸作用计算过程中发现,砖墙边界出现及墙面出现如图7(d)的剪切滑移。ESDOF模型则是采用固定边界的假定,故ESDOF模型在D1位置理论解比试验低了10.5%。而数值模拟采用了删除单元算法,所以构件质量下降刚度降低,导致模拟中D1位置的最大位移比试验大12.6%。总体而言,ESDOF模型与数值模拟相比可以更加准确地计算出砖墙的正向位移响应过程。
4. 结 论
基于一种改进的大型爆炸试验装置,开展了聚脲加固框架砖填充墙的原型爆炸试验,分析了爆炸荷载作用下加固砖墙的动力响应特征和破坏过程及模式,揭示了其失效破坏机理,建立了聚脲加固双向墙抗爆性能的理论计算模型,得到以下主要结论:
(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 阶段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 1630 0 3 620.500 2.34 ④区 介于③区和②区之间,参数按中心稀疏波规律衰减 2.34 ②区 28.29 138.05 4 730.127 692.510 2.34 ①区 28.29 7.20 4 730.127 2 345.567 1.40 ⓪区 0.101 1.225 0 340.294 1.40 -
[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. -