Generation of near-field blast wave by means of shock tube
-
摘要: 激波管可以在实验室环境下模拟爆炸产生冲击波,具有参数易于控制和测量手段准确多样等优势,在爆炸冲击效应的研究中被广泛应用。但与真实爆炸相比,尤其是近场爆炸,激波管产生的冲击波存在正压作用时间难以缩短、超压峰值难以提升的困难。通过对激波管运行理论和数值模拟分析发现:缩短正压作用时间的关键是让反射稀疏波尽快追上入射激波;提升超压峰值的关键是提高驱动气体的驱动能力。为此,设计了一种驱动段为锥形截面的激波管,使得反射稀疏波更快地追上入射激波,从而有效减小激波管设备长度并缩短正压作用时间;同时,采用正向爆轰驱动技术,利用化学能代替高压空气驱动提高驱动气体声速,在低爆轰初始压力下可以获得高的超压峰值。数值计算结果表明,在入射激波马赫数(MS=2.0)相同条件下,相对于等截面驱动方式,采用锥形截面驱动方式时,激波管长度可以减少近2/3,正压作用时间可以缩短近1/2。激波管实验结果表明,锥形截面驱动激波管产生的超压曲线满足近场爆炸冲击波形要求,并获得了超压峰值为64.7~813.4 kPa、正压作用时间为1.7~4.8 ms的爆炸冲击波波形。该研究可为近场爆炸冲击波致伤及装备防护效应评价实验提供参考。Abstract: Shock tubes can simulate blast waves in laboratory settings, offering advantages such as easily controlled parameters and varied measurement methods. It is widely used in the research of blast wave effects. However, in comparison to real blast, particularly in near-field blast, the blast waves generated by shock tubes has challenges in achieving shorter positive pressure durations and higher overpressure values. Through analysis of shock tube theory and numerical simulations, it has been determined that reducing positive pressure durations hinges on ensuring a swift catch-up by the reflected rarefaction wave with the incident shock wave. Similarly, increasing peak overpressure relies on enhancing the driving capability of the driving gas. Therefore, a conical cross-section driving approach is proposed to reduce the positive pressure durations, which allows the reflected rarefaction wave to catch up with the incident shock wave faster. By employing forward detonation driving technology and utilizing chemical energy to replace high-pressure air to increase the sound speed of the driving gas, high peak overpressure can be achieved at low detonation initial pressure. Numerical simulations show that under the same conditions of the incident shock Mach number (MS=2.0), the positive pressure durations can be reduced by nearly half and the device length can be reduced to nearly one-third by implementing the conical section-driven approach. Experimental results from the shock tube show blast wave characteristics, with peak overpressures ranging from 64.7 kPa to 813.4 kPa and positive pressure durations ranging from 1.7 ms to 4.8 ms. In blast wave simulation experiments, it is important to maintain the peak overpressure within a reasonable range to prevent the interface from reaching the test position. However, when the interface does reach the test position, it is possible to simulate the temperature field of the fireball in near-field blast waves. This research provides the necessary experimental conditions for evaluating the impact of near-field blast waves on injuries and investigating the protective performance of equipment.
-
Key words:
- blast wave /
- shock tube /
- detonation driving /
- near-field blast
-
作为多元合金的典型代表,CoCrFeNiMn多主元合金最早由Cantor等[1]发现,因此也将其称为Cantor合金。CoCrFeNiMn多主元合金是目前研究最多的FCC(face-centered cubic)多主元合金之一,其优异的力学性能使其成为低温和高温应用的潜在材料[2]。目前,关于CoCrFeNiMn多主元合金的研究主要集中在不同温度下的准静态分析,对于其在高温-高应变率耦合的极端服役条件下的研究还相对匮乏。在许多重要工程领域中,高温、高应变率是其结构服役过程中频繁面临的苛刻条件,例如发动机的高涡轮入口温度、叶片的高速旋转、起落架的高速冲击等,这些复杂的工况给飞机的飞行安全带来了严重的威胁[3-4]。新型多主元合金较传统合金表现出一种新的特性,即“鸡尾酒效应”。基于这一特性,可以进一步实现对材料性能的定向调控[5],例如通过添加Mo、Nb、Hf、Ta等难熔金属可制备耐高温合金[6],添加Al、Mg、Li等轻金属元素可制备轻质高强合金[7]。因此,新型多主元合金具有较好的强韧性、耐热性、抗腐蚀性、抗氧化性等,被视为航空工业领域极具潜在价值的合金材料[8-10],目前已在航空工业领域得到了一定的应用,如王秒等[11]将CoFeNiCrCu多主元合金作为钎料实现了高强度航空材料SiC陶瓷的钎焊连接,并获得了高温焊接条件下微观组织的演变规律。
针对涡轮叶片面临的高温-高应变率耦合的苛刻服役条件,众多研究者[12-13]一直致力于研发新的材料、发展新的制备工艺以提高承温能力,然而这些研究常常忽略材料在高应变率-高温耦合条件下的力学响应。例如,动态加载中材料由于声子拖曳效应表现出的强应变硬化能力[14]、高应变率与高温耦合下出现的动态应变时效[15]等,Jiang等[16]还发现高强度CoCrFeNiMn多主元合金在动态强迫剪切变形时会出现晶体的非晶化。因此,提高飞行器的结构安全不能仅局限在新型材料研发,明晰材料在极端环境中的变形机理、建立精度高且适用范围广的本构关系也是至关重要的环节。例如,飞鸟、冰雹等外物撞击时由于结构强度预测能力不足造成的飞行器失事已有报道[17],在类似的失效事故分析中,研究者多采用经典的Johnson-Cook(J-C)本构模型来开展数值仿真分析,或采用实验测试的手段来分析外物撞击飞机蒙皮和叶片的损伤失效规律以预测失效强度[18-19]。然而,这些研究大多数都局限在单方面考虑高速冲击的率效应或准静态的温度效应。虽然描述金属塑性变形的本构关系发展了多年,但率温本构关系应用最为广泛的模型仍然是J-C模型。J-C模型因形式简单、使用方便,被大量用于各种金属的工程实践及数值仿真计算中。为提高航空用38CrMoAl钢的使用效率,陈跃良等[20]对其开展了高速冲击下的力学性能研究,并获取了相关的J-C模型参数。Fu等[21]在高温高应变率条件下开展了Ti3Al合金的力热耦合分析,发现了材料的温度敏感指数不是一个恒定值,而是随着变形温度与塑性应变的增加而变化,证实了J-C模型对高温-高应变率耦合条件下材料应力-应变关系预测能力的不足。
对材料在高温、高应变率等苛刻条件下塑性变形行为的认识不足以及动态本构关系的缺乏,限制了高性能航空材料的探索及应用。本文中,采用新型CoCrFeNiMn多主元合金,开展其在高温-高应变率耦合加载条件下的力学性能测试,探索高温与高应变率对变形机制的影响,建立适用于宽温度域的动态本构关系,以期为材料在极端服役条件下的应用提供相应的分析方法和实验支撑。
1. 实验材料和方法
实验材料为CoCrFeNiMn多主元合金,材料加工流程为:铸锭→1473 K热处理12 h→1323 K锻造。为确定材料的化学组成及晶体结构类型,采用能谱仪(energy dispersive spectrometer, EDS)与X射线衍射仪(X-ray diffractometer, XRD)进行微观组成分析。结果显示,Ni、Fe、Cr、Co和Mn的质量分数分别为25.9%、20.1%、19.0%、17.9%和17.2%,与理论值的偏差较小,可以认为,该多主元合金的各元素是近等原子比。该多主元合金的XRD图谱(图1)显示,在43°、50°、75°附近出现了峰值,晶体结构为典型的面心立方晶系。
为进一步探索材料的初始微观结构,对材料开展电子背散射衍射(electron backscatter diffraction, EBSD)面扫分析。通过反极图(图2(a))可以看出,该多主元合金的初始晶粒取向随机分布,组织类型主要以变形态呈现(图2(b)的红色区域),同时还出现了部分回复组织(黄色)以及少量再结晶组织(蓝色)。材料的局部平均取向差(kernel average misorientation, KAM)较小,仅部分晶界处的KAM达到了2°左右(图2(c)),晶界角以大角晶界为主(图2(d)蓝色线条)。
采用配有高温炉的霍普金森压杆分别在298、673、873、1073和1273 K下进行动态加载实验。入射杆和透射杆的直径为14 mm,子弹长度为300 mm;温度控制装置主要包括调压器与控温器,其中温度传感器选用K型热电偶,高温炉尺寸为120 mm×70 mm×150 mm,炉内填充保温石棉,详细实验装置如图3(a)所示。为保证实验过程中试样温度的均匀性,将试样快速加热至目标温度,加热时间约为3 min,然后利用控温器将试样温度控制在目标温度±2 K的范围内并持续10 min。根据霍普金森压杆的实验标准,为获得较高的加载应变率,将试样的径高比设定为1.2,即直径为6 mm,高度为5 mm。不同温度下变形前后的试样如图3(b)所示,可以看出,试样的表观形貌未发生明显的改变,尤其是变形后试样未产生明显的鼓形,表明试样与入射和透射杆端的摩擦效应并不显著,进而证明了试样的整体变形在高度方向是均匀的。变形后试样的微结构变化将在第2节进行详细讨论。
通过气炮发射的子弹撞击入射杆的一端产生应力波,应力波沿入射杆、试样、透射杆进行传播,采用超动态应变仪采集入射杆和透射杆的应力波信号
εr 和εt 。基于一维应力波与应力均匀假设,可推导t时刻加载到试样的应力σ(t) 、应变ε(t) 、应变率˙ε(t) :{σ(t)=E0SAsεt(t)ε(t)=−2c0l0∫t0εr(t)dt˙ε(t)=−2c0l0εr(t) (1) 式中:E0和S分别为杆的弹性模量和横截面积,As和l0分别为试样的横截面面积和高度,c0为弹性波波速。
2. 实验结果和分析
实验中,所有试样采用同一尺寸,动态加载的应变率主要取决于撞击杆的加载速度,本研究重点分析动态加载中CoCrFeNiMn多主元合金的温度效应,不讨论其应变率效应,因此每次加载均控制撞击杆的加载速度约为20 m/s。冲击加载下的典型入射、透射波形如图3(a)的红线、绿线所示,根据式(1)可计算出不同温度下的真应力-真应变曲线,如图4所示。可以看出,在高速加载条件下,CoCrFeNiMn多主元合金表现出显著的温度效应,随着温度的升高,强度逐渐下降,当温度升高到1273 K后,其屈服强度仍能达到200 MPa,表现出了较好的耐高温性能。另外,在1273 K的环境中,CoCrFeNiMn多主元合金进入塑性变形阶段后依旧具有较好的应变硬化能力。
为了进一步探明高应变率条件下温度对CoCrFeNiMn多主元合金塑性变形机制的影响规律,分别对4种温度(673、873、1073和1273 K)时的试样沿高度方向对中心区域进行EBSD面扫分析,结果如图5~6所示。为避免高速加载后高温炉内热量对变形后试样微观结构的影响,加载完成后,采用图3(a)所示的试样夹持装置将试样快速拆卸,然后进行微结构表征分析。通过反极图(图5(a)~(d))可以看出,随着温度的升高,晶粒尺寸较初始形貌出现了明显的改变,673 K时,塑性变形后材料的晶粒发生了细化,出现了大量尺寸小于5 μm的晶粒(图5(a)),还产生了大量的再结晶组织(图5(e))。文献[22]表明,准静态条件下CoCrFeNiMn多主元合金的再结晶温度约为1073 K,可初步推断高应变率加载可能会降低CoCrFeNiMn多主元合金的再结晶温度。当温度升高到873 K后,局部晶粒开始重新长大,出现了大量的回复组织(图5(f)),然而1073 K时样品发生塑性变形后不再出现5 μm以下的晶粒(图5(c)),回复组织与再结晶组织也较少,主要以变形组织为主(图5(g))。动态加载过程中,受绝热温升的影响,动态再结晶的初始温度降低;但在动态变形过程中,由于响应时间较短以及材料固有的缓慢位错动力学特性和较低的层错能,材料的连续动态再结晶过程被抑制,进而导致高应变率变形过程中的动态再结晶能力较低(相较于准静态)。Chen等[23]在CoCrFeNiMn多主元合金(1073 K)的动态压缩测试中也观察到了类似的现象。当环境温度升高到1273 K时,塑性变形后晶粒尺寸均为10~20 μm(图5(d)),而且出现了大尺寸的再结晶组织(图5(h)),该温度下不连续动态再结晶的出现可能导致再结晶晶粒长大。Khan等[24]也发现AlTiNbZrW合金在1273 K时的热变形将同时出现连续动态再结晶和不连续动态再结晶。
图6为不同环境温度下塑性变形后CoCrFeNiMn多主元合金的KAM分布图(图6(a)~(d))和晶界角分布图(图6(e)~(h)),直接反映了几何必要位错(geometrically necessary dislocation, GND)在试样中的分布。可以看出,673 K时,塑性变形后合金具有较大的KAM,KAM大多均集中在2°左右,部分晶界处的取向差达到5°(图6(a)),同时材料内部还出现了大量的小角晶界(图6(e)的绿色线条),表明该温度下塑性变形后试样内部存储较高水平的位错密度,由于此时温度较低,较高的位错密度主要归因于塑性变形过程中的应变硬化效应。当环境温度为873和1073 K时,塑性变形后材料内部的KAM较673 K时更小(图6(b)~(c)),这是由于当温度达到873 K时,温度导致部分位错湮灭,图5(f)~(g)所示的回复组织也可以证实这一推断,与此同时,材料内部小角晶界的增加量也在逐渐降低。另外,通过KAM图可以发现,变形晶粒边缘区域的取向差高于内部区域,表明变形过程中晶界附近比内部区域具有更强的亚结构孕育能力。值得注意的是,当温度达到1273 K时,试样内部的KAM较873和1073 K时更高,表明由于高温与高应变率的耦合作用,GND增殖,进而导致合金具有较高的应变硬化能力[25],1273 K的宏观应力-应变曲线可以证明这一推断。此外,由图5(h)可以发现,1273 K时动态变形后材料内部出现了较大尺寸的动态再结晶结构,大尺寸再结晶结构的形成过程会湮灭大量的位错,变形进一步集中在晶界中,导致位错分布不均匀,表现为较大的KAM。Soares等[26]对动态压缩中CoCrFeNiMn微结构的分析也得到了类似的结论。
高温-高应变率耦合条件下,CoCrFeNiMn多主元合金的应变硬化能力提升,考虑温度效应的本构模型(例如J-C模型)对材料温敏指数以及应力-应变关系的预测能力下降。
3. 本构建模与参数分析
为加速CoCrFeNiMn多主元合金在极端环境中的应用,在探明其塑性变形机理的基础上还需建立动态本构模型,以便仿真计算。动态本构建模的前提就是精确获取材料变形的相关参数以及准确还原塑性变形过程。材料以高应变率进行变形时,没有足够的时间使变形热从材料中消散,这有别于准静态条件下的等温过程。因此,绝热温升成为高温动态本构建模的重要影响因素,建模时,既需要考虑环境温度效应也需要考虑绝热温升。然而,目前多数的高温动态本构建模研究中仅考虑了环境的温度效应,忽略了绝热温升。通常,绝热条件下塑性变形引起的温升
ΔT 可表示为:ΔT=T−T0=ηρcV∫εpεyσdε (2) 式中:T为瞬时温度,即试样的真实变形温度;T0为初始温度,即环境温度;
η 为功热转化效率,通常取0.9;ρ 和cV分别为材料的密度和比定容热容,对于CoCrFeNiMn多主元合金,通常取ρ = 8042 kg/m3,cV = 0.43 J/(g·K)[27];εy 和εp 分别为初始屈服应变和t时刻的塑性应变。通过式(2),可以绘制出不同环境温度下塑性变形与绝热温升的关系,如图7所示。通过多次数值拟合分析,可以进一步将温度增量统一归纳为动态压缩塑性应变与环境温度的函数:ΔT=(163.42−0.26T0+3.23×10−4T20−1.35×10−7T30)(expεp3.09×10−4T0+0.58−1) (3) 式(3)中的已知数值均由拟合得到。首先,在同一变形温度下获得绝热温升与塑性应变的关系,然后将不同温度下拟合所得的绝热温升参数再次与变形温度进行拟合分析,最终可得塑性应变与环境温度共同影响的绝热温升函数。
构建考虑材料温度效应的动态本构关系时,精确确定材料的温度敏感指数是非常关键的一步,目前研究中关于材料温度敏感指数的确定大多未考虑塑性变形以及绝热温升。基于式(2)~(3)获得不同时刻的变形温度,将温度敏感指数(m)修正为如下形式:
m=ln[σT(εp)/σTr(εp)]ln[T(εp)/Tr(εp)] (4) 式中:
σT 和σTr 分别为变形温度和参考温度下的真应力;Tr(εp) 为参考温度,即室温下考虑了绝热温升后的真实变形温度。根据式(4)可进一步计算不同变形温度和塑性应变时的温度敏感指数,如图8所示。由图8可知,不同塑性应变水平下m并不是恒定不变的,不同变形温度下随着塑性应变的增大,m出现了不同程度的波动。当环境温度处于673~1073 K时,m在−0.45上下波动,而且波动幅值不超过0.1;然而当环境温度达到1273 K时,m出现了明显的下降,塑性应变从0增大到0.20后,m从−0.75增大到−0.50。Park等[28]在准静态范围获得的CoCrFeNiMn多主元合金的m为−0.60。可以看出,使用幂律函数关系构建的本构关系(如J-C模型)描述CoCrFeNiMn多主元合金的温度效应时,仅能较准确地描述1073 K以下的动态应力-应变响应,对1273 K时合金动态应力-应变关系的预测能力显著降低。
考虑到塑性变形过程对温度效应的影响,首先要分析初始屈服强度的温度效应,不同变形温度下CoCrFeNiMn多主元合金的初始屈服强度如图9所示。可以看出,随着变形温度的升高,初始屈服强度不断地降低。将初始屈服应力
σ0 分为温度敏感部分σt 和温度不敏感部分σat ,结合温度效应的指数变化率,初始屈服应力与温度的关系可进一步表示为:σ0=σte−βT+σat (5) 式中:
β 为温度影响系数。通过对实验数据的拟合,可得σt 、β 、σat 分别为733.16 MPa、1.23×10−3 K−1、115.42 MPa。基于式(5)和图9,可以推断CoCrFeNiMn多主元合金不受温度影响的屈服强度约为115.42 MPa,这部分应力不受热激活的影响;另一方面,可预测合金在绝对零度下的动态屈服强度约为848.58 MPa。结合式(5)和动态应力-应变曲线的特征,在初始屈服应力的热敏感部分引入金属材料最常用的幂强化模型,可以初步将考虑温度效应的应力-应变关系表示为:
σ=(σt+Aεnp)e−βT+σat (6) 式中:A和n分别为应变硬化系数和应变硬化指数。采用式(6)对图4中298 K时的动态塑性应力-应变曲线进行拟合,可得到A=567.88 MPa,n=0.69;然后,取n恒等于0.69,拟合可得其余4种温度下的A分别为739.82、959.85、1206.15和1435.03 MPa。参数A与T的关系如图10所示。
通过拟合可将A与T的关系定量表示为:
A=aeβ1T+b (7) 式中:a、
β1 、b均为温度相关的材料参数,分别取226.01 MPa、1.36×10−3 K−1、204.90 MPa。将式(7)代入式(6),可得:σ=[σt+(aeβ1T+b)εnp]e−βT+σat (8) 整理后可得考虑温度效应的CoCrFeNiMn多主元合金的动态塑性应力-应变关系:
σ=σte−βT+[ae(β1−β)T+be−βT]εnp+σat (9) 对于CoCrFeNiMn多主元合金,式(9)中的参数取值列于表1。根据式(9)和表1,绘制不同温度下CoCrFeNiMn多主元合金的塑性应力-应变曲线,其与实验结果的对比如图11所示。可以看出,理论结果与实验结果吻合较好,验证了这一本构关系的有效性。另外,通过图11还可以看出,指数形式的单轴本构关系的精度较高,其对高温下出现的热硬化现象具有非常好的描述能力,这一建模思路也可以进一步推广到其他金属材料。
表 1 模型参数Table 1. Parameters of proposed modelσt/MPa β/K−1 a/MPa β1/K−1 b/MPa n σat/MPa 733.16 1.23×10−3 226.01 1.36×10−3 204.90 0.69 115.42 本研究建立的单轴高温动态本构方程将初始屈服应力和塑性硬化应力的温度效应进行了解耦。对于初始屈服阶段,基于数值拟合确定了初始屈服应力部分中的热不敏感应力。对于塑性硬化阶段,将塑性应力的温度效应分解为两部分进行表征:一部分为常规的热软化项;另一部分为热硬化项,用于解决部分FCC合金在高温动态加载过程中由于动态应变时效、声子拖曳机制、相对论效应等引起的热硬化现象。
4. 结 论
研究了不同高温环境下CoCrFeNiMn多主元合金的动态力学响应分析及塑性变形机理,给出了不同温度下绝热温升的发展规律,并建立了一个唯象的计及温度效应的本构方程,得到的主要结论如下。
(1)在20 m/s的加载条件下,CoCrFeNiMn多主元合金的塑性变形表现出显著的温度敏感性。较低温度下,材料变形后内部将存储较高水平的位错密度;随着温度的升高,内部存储的位错将不断湮灭。变形晶粒边缘区域的局部平均取向差高于内部区域,变形过程中晶界附近比内部区域具有更强的亚结构孕育能力。当变形温度达到1273 K时,晶粒出现了明显的粗化现象,晶粒尺寸较变形前更大。
(2)动态加载过程中,绝热温升随塑性应变的增大呈指数形式变化,通过引入环境温度项定量分析了塑性变形和环境温度对绝热温升的影响。综合考虑绝热温升和环境温度的量化分析方法精确捕获了CoCrFeNiMn多主元合金的温度敏感指数,证实了幂律形式的本构关系无法精确描述该材料在宽温度域内的应力-应变关系。通过对初始屈服应力和塑性流动阶段的解耦分析,并引入热不敏感应力和热硬化项,建立了指数形式的本构关系。该本构方程对CoCrFeNiMn多主元合金在298~1273 K的动态应力-应变关系有较高的预测能力。
-
表 1 在H2和O2充气物质的量的比为3∶1和不同初始压力条件下正向爆轰驱动实验获得的超压峰值和正压作用时间
Table 1. Peak overpressure and positive pressure action time obtained in positive detonation driving experiment at different initial inflation pressures and n(H2)∶n(O2)=3∶1
实验状态 初始压力/MPa 超压峰值/kPa 正压作用时间/ms 1 0.55 490.4 4.7 2 0.60 539.5 4.8 3 0.65 624.2 4.8 4 0.85 813.4 4.4 -
[1] ELDER G A, CRISTIAN A. Blast-related mild traumatic brain injury: mechanisms of injury and impact on clinical care [J]. Mount Sinai Journal of Medicine: A Journal of Translational and Personalized Medicine, 2009, 76(2): 111–118. DOI: 10.1002/msj.20098. [2] TURNER R C, NASER Z J, LOGSDON A F, et al. Modeling clinically relevant blast parameters based on scaling principles produces functional & histological deficits in rats [J]. Experimental Neurology, 2013, 248: 520–529. DOI: 10.1016/j.expneurol.2013.07.008. [3] RISDALL J E, MENON D K. Traumatic brain injury [J]. Philosophical Transactions of the Royal Society B: Biological Sciences, 2011, 366(1562): 241–250. DOI: 10.1098/rstb.2010.0230. [4] HERNANDEZ A, TAN C F, PLATTNER F, et al. Exposure to mild blast forces induces neuropathological effects, neurophysiological deficits and biochemical changes [J]. Molecular Brain, 2018, 11(1): 64. DOI: 10.1186/s13041-018-0408-1. [5] NING Y L, ZHOU Y G. Shock tubes and blast injury modeling [J]. Chinese Journal of Traumatology, 2015, 18(4): 187–193. DOI: 10.1016/j.cjtee.2015.04.005. [6] BAKER W E. Explosions in air [M]. Austin: University of Texas Press, 1973. [7] CLEMEDSON C J, CRIBORN C O. A detonation chamber for physiological blast research [J]. Journal of Aviation Medicine, 1955, 26(5): 373–381. [8] FILLER W S. Propagation of shock waves in a hydrodynamic conical shock tube [J]. Physics of Fluids, 1964, 7(5): 664–667. DOI: 10.1063/1.1711266. [9] STEWART J B, PECORA C. Explosively driven air blast in a conical shock tube [J]. Review of Scientific Instruments, 2015, 86(3): 035108. DOI: 10.1063/1.4914898. [10] COURTNEY A C, ANDRUSIV L P, COURTNEY M W. Oxy-acetylene driven laboratory scale shock tubes for studying blast wave effects [J]. Review of Scientific Instruments, 2012, 83(4): 045111. DOI: 10.1063/1.3702803. [11] COURTNEY M W, COURTNEY A C. Note: a table-top blast driven shock tube [J]. Review of Scientific Instruments, 2010, 81(12): 126103. DOI: 10.1063/1.3518970. [12] CELANDER H, CLEMEDSON C J, ERICSSON U A, et al. The use of a compressed air operated shock tube for physiological blast research [J]. Acta Physiologica Scandinavica, 1955, 33(1): 6–13. DOI: 10.1111/j.1748-1716.1955.tb01188.x. [13] CULBERTSON D W. Description and performance of a conical shock tube nuclear air blast simulator [C]// Proceedings of the Seventh International Shock Tube Symposium. Toronto: University of Toronto Press, 1970: 396–409. DOI: 10.3138/9781487595876-024. [14] OPALKA K O, MARK A. The BRL-Q1D code: a tool for the numerical simulation of flows in shock tubes with variable cross-sectional areas: AD-A139631 [R]. Aberdeen: U. S. Army Ballistic Research Laboratory, 1986. [15] YU H R, GU J H, LI Z F, et al. Generation of blast wave by means of the normal shock tube [C]//Proceedings of the International Symposium on Shock Waves. Sendai, Japan, 1992: 897–900. [16] 王正国, 孙立英, 杨志焕, 等. 系列生物激波管的研制与应用 [J]. 爆炸与冲击, 1993, 13(1): 77–83. DOI: 10.11883/1001-1455(1993)01-0077-7.WANG Z G, SUN L Y, YANG Z H, et al. The design production and application of a series of bio-shock tubes [J]. Explosion and Shock Waves, 1993, 13(1): 77–83. DOI: 10.11883/1001-1455(1993)01-0077-7. [17] KIRK D R, FAURE J M, GUTIERREZ H, et al. Generation and analysis of blast waves from a compressed air-driven shock tube [C]//38th Fluid Dynamics Conference and Exhibit. Seattle: AIAA, 2008: 4777. DOI: 10.2514/6.2008-3847. [18] KLEINSCHMIT N N. A shock tube technique for blast wave simulation and studies of flow structure interactions in shock tube blast experiments [D]. Lincoln: The University of Nebraska, 2011. [19] NGUYEN T T N, WILGEROTH J M, PROUD W G. Controlling blast wave generation in a shock tube for biological applications [J]. Journal of Physics: Conference Series, 2014, 500: 142025. DOI: 10.1088/1742-6596/500/14/142025. [20] ANDREOTTI R, COLOMBO M, GUARDONE A, et al. Performance of a shock tube facility for impact response of structures [J]. International Journal of Non-Linear Mechanics, 2015, 72: 53–66. DOI: 10.1016/j.ijnonlinmec.2015.02.010. [21] LI X D, HU Z M, JIANG Z L. Numerical investigation of the effects of shock tube geometry on the propagation of an ideal blast wave profile [J]. Shock Waves, 2017, 27(5): 771–779. DOI: 10.1007/s00193-017-0716-x. [22] FRIEDLANDER F G. The diffraction of sound pulses Ⅰ: diffraction by a semi-infinite plane [J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1946, 186(1006): 322–344. DOI: 10.1098/rspa.1946.0046. [23] LUO K, WANG Q, LI J W, et al. Numerical modeling of a high-enthalpy shock tunnel driven by gaseous detonation [J]. Aerospace Science and Technology, 2020, 104: 105958. DOI: 10.1016/j.ast.2020.105958. -