饱和水泥试样被爆炸激波损伤破碎的尺度研究

利用水中爆炸冲击波使水泥试样损伤破坏,模拟爆炸采油时激波使岩石损伤开裂的现象。实验获得了适合本实验条件的激波峰压衰减规律pm8.2(W1/3/R)1.46,得知压碎区尺度为集中装药特征尺度的2～5倍、拉伸损伤区尺度为集中装药特征尺度的20～30倍,激波使水泥试样破碎、拉裂的能量占总能量的2%～7%。     

球形重质气体物理爆炸特性

基于大涡模拟方法,结合高阶混合格式,对高压重质的SF₆球形气云在空气中爆炸进行了模拟。数值模拟表明,爆炸产生的激波经过气体分界面时分为透射激波以及反射稀疏波,透射激波导致气体分界面处Richtmyer-Meshkov失稳增强,从而加速了2种气体的混合,而反射的稀疏波经过汇聚,在球心处形成二次激波,在该强激波作用下,流场区域基本呈现湍流形态。     

具有引力场的气体中激波传播的近似理论

本文在三维不定常流体力学特征理论的基础上,运用CCW方法,获得了激波速度的变化和激波阵面曲率、引力的大小和方向之间的解析关系,讨论了引力场对激波的加强和减弱的影响以及超新星爆炸中反弹激波的问题。利用得到的激波面积变化公式,给出了二、三维激波在引力场作用下传播问题的模拟处理方法。作为简单的例子,具体地研究了引力场变密度气体中强激波传播问题。     

柱形装药条件下锥形水中爆炸激波管内的冲击波特性

锥形水中爆炸激波管是进行水中爆炸实验的一种装置,该装置能够通过较小装药量在相同距离处实现自由场水中较大装药量爆炸的冲击波峰值。为了获得柱形装药条件下锥形水中爆炸激波管内的冲击波特性,本文通过数值计算的方式,对不同圆锥角和不同柱形装药质量下锥形激波管内的冲击波传播过程进行了模拟,通过对不同工况下激波管内冲击波特性进行分析,发现其初始冲击波的衰减规律符合自由场水中的指数衰减形式,并拟合得到了与自由场水中爆炸相容的冲击波峰值、比冲量和能流密度经验公式;发现其二次脉动压力周期与炸药质量呈反常规的变化规律,并引入等效静水压深度解释了这一现象;发现其二次脉动压力幅值与初始冲击波幅值之比比自由场水中更大,而二次脉动压力的比冲量与初始冲击波冲量之比与自由场水中相当。     

单向砌体填充墙激波管试验和动力行为分析

为研究远场爆炸荷载作用下单向砌体填充墙的动力行为及其失效破坏机理,首先,基于研发的压缩空气驱动大截面（3 m×3 m）激波管开展了两面不同厚度单向实心砌体填充墙的面外加载试验,获得了作用在墙体表面的反射超压荷载时程、墙体面外挠度时程及墙体的变形失效模式。其次,建立了激波管精细化有限元模型,提出了砌体墙简化微观有限元建模方法,以及扩展砖块的Riedel-Hiermaier-Thoma材料模型和接缝的内聚力接触模型参数取值计算方法,对激波管中的压力传播以及试验墙体面外动态响应和损伤破坏开展了数值模拟。最后,基于爆炸荷载作用下单向砌体填充墙面外抗力方程和等效单自由度模型对试验墙体中心点面外挠度时程进行预测。结果表明：减小墙体高厚比可以增大框架拱推力,从而显著提升墙体的抗爆性能,厚105 mm的墙体在经历1次激波管试验后发生倒塌,而厚235 mm的墙体在经历6次激波管试验后仅有轻微损伤;墙体表面反射超压荷载时程的试验和数值模拟结果均为均布脉冲型荷载,且两者吻合很好,验证了激波管设计和激波管精细化有限元模型的合理性;数值模拟和理论预测的墙体动力行为与试验结果吻合较好,可为砌体填充墙抗爆评估与分析...     

空气自由场中强爆炸冲击波传播二维数值模拟

编写了适用于模拟具有高密度比、高压力比的强激波问题的二维柱对称多介质流体计算程序。利用有限体积方法求解流体的Euler方程组,采用level set方法捕捉爆炸产物与空气的运动界面,并通过求解物质界面两侧Riemann问题的精确解来计算爆炸产物与空气之间的数值通量。研制了三角形网格自适应技术来实现网格的自动加密和粗化,在保证捕捉激波峰值的前提下有效地提高了计算效率。利用计算程序对1 kt TNT当量的空气自由场强爆炸问题进行数值模拟,计算得到的峰值超压、冲击波到达时间等物理参数与点爆炸理论结果基本一致。     

激波在波前有非均匀定常流动的变截面管道中传播及波后流场特性的数值研究

本文采用分裂算子的随机选取法,对激波在波前有非均匀定常流动的一维变截面管道中的传播进行了数值模拟,得到了激波强度和波后流动的变化规律,并对结果作了初步分析。     

环形激波绕射、反射和聚焦流场的间断有限元模拟

采用间断有限元方法对环形激波在圆柱形激波管内绕射、反射和聚焦流场进行了数值模拟。将二维守恒方程的间断有限元方法发展到轴对称Euler方程,并对环形激波绕后台阶流动进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉运动激波在圆柱形激波管内传播的复杂流场结构;在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力,在聚焦点附近不会产生振荡或抹平间断现象。     

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成;(2) U型空泡团头部溃灭;(3) U型空泡团腿部溃灭.在U型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象,抑制了下一周期的空化发展.     

冷激波灭火系统中激波对灭火效果和周边环境的影响

为了评估冷激波灭火弹爆炸后形成的激波对灭火效果和周边环境的影响,建立了一套纹影实验装置。通过纹影实验,观察了小尺寸下不同灭火介质爆炸抛撒后激波的形成和传播,并推算了冷激波灭火弹爆炸后油盆边缘附近的波后质点速度。观察发现,水基灭火介质爆炸抛撒后没有形成激波;而粉基灭火介质尽管爆炸抛撒后形成激波,但在油盆边缘附近激波强度较弱,波后质点速度小。与介质抛撒引起的可燃气体介质运动相比,激波对加快可燃气体介质运动的影响可以忽略。最后,通过高速摄影实验验证了激波对灭火效果的影响是可以忽略的。因此,冷激波系统中,激波对灭火效果和周边环境的影响是可以忽略的。     

爆炸波在非饱和钙质砂中的传播规律

以球形TNT药包作为爆源，在密实的非饱和钙质砂中进行了一系列大尺寸爆炸模型试验。主要研究在不同药包质量、埋深及砂土试样含水率条件下，密实钙质砂中爆炸波的主要基本参数随传播距离增加而变化的规律。试验结果表明爆炸波主要以弹塑性波的形式在密实钙质砂中传播，在干燥和潮湿试样中塑性纵波波速随试样初始密度的增大或含水率的降低而增大，且范围分别为250～282 m/s和302～339 m/s。集团装药情况下，非饱和钙质砂中封闭爆炸的临界比例埋深约为2.25 m/kg1/3。在试验范围内，密实钙质砂中爆炸波的法向应力峰值及法向比冲量的衰减均服从爆炸相似律。封闭爆炸时，干燥钙质砂中爆炸波的应力衰减指数在测点比例爆心距大于或小于0.75 m/kg1/3处分别为2.94或1.37；潮湿钙质砂中爆炸波的应力衰减指数随含水率升高而增大，其范围为1.39～1.79。法向比冲量衰减指数随试样含水率升高而减小，其范围为0.97～1.18。     

低温和低压环境下炸药爆炸冲击波的传播特性

针对高海拔或高空的低温、低压环境对炸药爆炸冲击波传播的影响,利用量纲分析理论和AUTODYN有限元软件,研究了低温、低压及海拔高度对炸药爆炸冲击波参量（峰值超压、比冲量和波阵面运动轨迹）的影响规律,建立了相应的计算公式,并通过数值模拟和实验数据进行了对比验证。结果表明,该计算公式可以有效预测低温和低压环境下炸药爆炸冲击波参量。环境压力降低,爆炸冲击波峰值超压和爆炸远场（比例距离Z>0.2 m/kg1/3）比冲量减小,冲击波传播速度增大。环境温度降低,冲击波比冲量增大,传播速度降低,峰值超压影响不大。海拔高度在0～9 000 m范围内,每升高1 000 m冲击波峰值超压和爆炸远场比冲量分别平均降低约3.9%和3.2%。海拔升高,爆炸近场冲击波传播速度升高,爆炸远场冲击波传播速度则降低。高海拔环境下低压对冲击波峰值超压和比冲量的影响大于低温,爆炸近场冲击波传播速度取决于低压的影响,爆炸远场冲击波传播速度取决于低温的影响。     

Effect of laser supported detonation wave confinement on termination conditions

A laser supported detonation (LSD) wave was driven using line-focusing laser optics, in which an induced blast wave expanded laterally from the LSD region to surrounding air in two-dimensional space. The LSD wave was confined in quasi-1D space using a wedge nozzle to restrict the lateral expansion of a blast wave. The LSD termination threshold and the blast wave energy were deduced from shadowgraphs showing the blast wave expansion. The respective threshold laser intensities for cases with and without confinement were estimated as 17 and 34 GW/m², indicating that the lateral expansion strongly influenced on the LSD termination condition.      

空中强爆炸冲击波地面反射规律数值模拟研究

为准确预测空中强爆炸产生的冲击波载荷分布，基于欧拉坐标系建立了能够模拟具有高密度比、高压力比的强激波问题的二维多介质流体数值方法。结合网格自适应技术，对1 kt TNT当量的空中强爆炸在不同爆炸高度下的冲击波地面反射过程进行了数值模拟，并考虑了真实气体状态方程和空气随高度不均匀分布的影响。计算得到了地面上距爆心投影点大尺度范围内的反射超压和冲量等冲击波载荷分布，并给出了冲击波载荷随爆高的变化规律。     

Numerical investigation of inviscid shock wave dynamics in an expansion tube

Highly complicated shock wave dynamics has been numerically calculated by solving the Euler equations for a circular shock tube suddenly expanded three times of the original tube diameter atx=0. Shock waves of different shock Mach number,M _s =1.5 and 2.0, have produced remarkably distinct blast jet structures. A planar shock wave took its final form after the blast by repeated Mach reflections of the blast wave: the first one at the wall and the second one at the central axis. The central Mach disc overtook and merged with the annular Mach stem before the planar shock wave was formed. In contrast to the blast wave which would propagate spherically in an open space, the present blast wave undergoes complex morphological transformation in the restricted flow passage, resulting in an unstable and oscillatory blast jet structure of highly rotational nature. The slipstream tube emanating from the shock tube exit corner decomposed into a chain of small vortex rings that interacted with the barrel shock of the jet, which caused periodic collapse of the jet structure. The finite volume-FCT formulation equipped with the time-dependenth-refinement adaptive unstructured triangular mesh technique in the present paper has contributed to resolution of the intricate physical discontinuities developing in the blast flow fields.     

Failure mechanisms of concrete slab–soil double-layer structure subjected to underground explosion

The failure mechanism of a concrete slab–soil double-layer structure subjected to an underground explosion was investigated by experimental and numerical methods in this paper. Two underground explosion depths of 150 and 350 mm were tested. The typical failure modes such as the conoid spall of concrete, the bulge of the concrete slab and the cavity in the soil were obtained experimentally. Numerical simulations of the experiments were performed using a hydrodynamic code to analyze the effects of both the stress wave and the expansion of the blast products. Based on the experimental and numerical results, the effects of explosive depth, blast wave front and expansion of the blast products on the failure modes and failure mechanisms were discussed. The underground explosion process at different explosion depths was also analyzed. The results show that attenuation of the stress wave in the soil is significant. The blast wave front and the expansion of the blast products play different roles at different explosion depths. At the explosion depth of 150 mm, the failure mode is mainly caused by a point load induced by the blast wave front, whereas at the depth of 350 mm a sphere-shaped load resulting from the expansion of the blast products is a key factor for failure.     

高电压点火有效能量的测量及相关问题

建立了测量高电压点火产生的爆炸波参数装置,利用PCB压力传感器测定了在空气中高压电火花产生的爆炸波压力及其到达时间。改变点火能量和测试距离,得到了爆炸波的变化趋势和传播规律,并探讨了爆炸波各参数的标度尺寸。结果表明,电火花产生的爆炸波符合Hopkinson-Sachs爆炸相似律。把1/4周期放电能量与数值模拟结果进行比较,得知初始1/4周期放电能量(点火能量)与爆炸波能量基本吻合;因此,1/4周期放电能量可被视为作用于直接起爆引起爆轰的有效能量。进一步探讨了高电压点火各能量的分布;结果表明,储存于电容的总能量1/2CU2中约91%消耗于能量损失和欧姆损耗,初始1/4周期放电能量仅约占2%,其余的点火能量只起到加热混合物的作用。     

冲击波对核反应堆安全壳的动力响应研究

本文用时域边界积分方程提出研究半球形核反应堆安全壳在爆炸冲击波作用下动力响应分析的计算模型。从空中爆炸和地面爆炸两方面探求壳体位移和变形过程,分析冲击波在壳体内的传播规律。计算了在水平冲击波作用下壳体的内力分布,其计算结果可以作为我国核电站设计的参考分析依据。     

Evolution of blast wave profiles in simulated air blasts: experiment and computational modeling

Shock tubes have been extensively used in the study of blast traumatic brain injury due to increased incidence of blast-induced neurotrauma in Iraq and Afghanistan conflicts. One of the important aspects in these studies is how to best replicate the field conditions in the laboratory which relies on reproducing blast wave profiles. Evolution of the blast wave profiles along the length of the compression-driven air shock tube is studied using experiments and numerical simulations with emphasis on the shape and magnitude of pressure time profiles. In order to measure dynamic pressures of the blast, a series of sensors are mounted on a cylindrical specimen normal to the flow direction. Our results indicate that the blast wave loading is significantly different for locations inside and outside of the shock tube. Pressure profiles inside the shock tube follow the Friedlander waveform fairly well. Upon approaching exit of the shock tube, an expansion wave released from the shock tube edges significantly degrades the pressure profiles. For tests outside the shock tube, peak pressure and total impulse reduce drastically as we move away from the exit and majority of loading is in the form of subsonic jet wind. In addition, the planarity of the blast wave degrades as blast wave evolves three dimensionally. Numerical results visually and quantitatively confirm the presence of vortices, jet wind and three-dimensional expansion of the planar blast wave near the exit. Pressure profiles at 90° orientation show flow separation. When cylinder is placed inside, this flow separation is not sustained, but when placed outside the shock tube this flow separation is sustained which causes tensile loading on the sides of the cylinder. Friedlander waves formed due to field explosives in the intermediate-to far-field ranges are replicated in a narrow test region located deep inside the shock tube.