基于饱和响应时间的封闭空间内爆炸载荷等效方法研究

针对舱内爆炸载荷形式复杂、作用时间长、缺乏有效的简化描述方法的问题，首先采用显式动力学计算程序开展了内爆载荷作用下钢板动态响应的数值计算，在与试验结果对比验证的基础上分析了金属板的内爆载荷饱和冲量。通过对216种不同爆炸载荷加载时长与金属板响应关系的分析，提出了内爆炸载荷作用下结构最大变形所对应的饱和时间计算经验公式，并给出了饱和时间的无量纲系数建议值。考虑到内爆载荷初始冲击波的影响，结合爆炸载荷饱和作用时间的规律，提出了封闭空间爆炸载荷的矩形载荷等效方法，对比了18组简化载荷与耦合载荷分别作用下钢板的动力响应，验证了等效方法的有效性。     

舱室内爆下舰船结构损伤的一种计算方法

为了评估舱室内爆多载荷耦合作用下舰船结构的损伤范围,设计了大尺度舱段模型,并开展了舱室内爆毁伤试验,试验后测量了舱室结构的破坏范围及破坏模式,分析了舱室内爆多载荷耦合作用下舰船结构的损伤机理,据此建立了舱室内爆下舰船结构损伤的计算方法。结果表明：(1)舱室内爆下形成的强冲击波载荷和准静态压力载荷可对舰船结构造成大范围的损伤,形成多种破坏模式;(2)舱室内爆下准静态压力载荷是舱室结构损伤破坏的主要毁伤元;(3)建立的舱室内爆载荷下结构损伤变形计算方法可同时考虑强冲击载荷和准静态压力载荷对结构的损伤破坏,理论计算结果与试验结果吻合较好。     

坑道内化爆冲击波的传播规律

装药在坑道内部爆炸时将产生沿坑道传播的空气冲击波,由于受到坑道洞壁的限制,空气冲击波在坑道内的传播规律异于在自由大气中的传播。利用三维数值模拟计算程序,对长坑道中的化爆流场进行了数值计算,归纳出空气冲击波沿坑道方向的传播规律。计算结果与试验结果符合较好。     

Whitham公式在球形内爆实验中的适用性研究

球形内爆实验是实现超高压状态方程测量的方法之一,Whitham公式可对球面聚心冲击波的传播特性进行描述。文中给出了确定Whitham公式系数的方法,针对使用两种不同厚度样品球形内爆实验,通过Whitham系数不确定度分析及其求解方程组的性态分析,探讨了Whitham公式在球形内爆实验中的适用范围。结果表明,样品厚度过薄时,Whitham公式不能应用;当0．8〈A1=R1／R0〈0．857且A2=R2／R0〈0．657时,在内爆实验的结果分析和处理中应用Whitham公式可能会取得较好的结果。     

基于爆能等效原理大型模爆器燃气爆炸冲击加载的数值模拟

运用非线性显式动力有限元程序LS-DYNA,基于多物质Euler算法,对TNT炸药和乙炔-空气混合气体两种爆炸源在自由大气场中爆炸产生的冲击波荷载特征参数进行数值模拟,比较两种爆源产生的冲击波压力传播规律。基于爆能等效原理,按超压相等的原则给出了气体爆炸名义比例距离计算公式。结果表明,基于Euler算法可以较好地描述乙炔-空气混合气体爆炸空气冲击波传播规律,爆炸压力随着距爆源距离的增大而迅速衰减,且两种爆源产生的冲击波超压峰值误差随着冲击波传播距离的增大而逐渐减小。采用名义比例距离公式修正后,气体爆炸与炸药爆炸冲击波计算误差可以得到有效控制。当爆炸冲击波超压小于0.5MPa时,可以采用乙炔-空气混合气体代替化学炸药进行模爆器内爆炸实验加载。     

基于LS-DYNA的液电效应冲击波数值模拟

液电效应机理复杂,鲜有成熟的商用数值模拟软件能够描述等离子体通道内部特性,为了将液电效应产生的冲击波运用于已有的数值模拟软件中,以满足工程需要,介绍了两种基于显式动力学软件LS-DYNA间接模拟液电效应产生冲击波的方法:水下爆炸等效（分为爆炸能量等效与冲击波能量等效）和理想气体等效,并进行了比较与改进,分析了不同沉积能量下采用不同等效方法得到的峰压计算结果的差异。结果显示,在沉积能量相同的条件下,基于爆炸能量等效方法得到的冲击波峰值压力最高,基于冲击波能量等效方法得到的冲击波峰值压力次之,基于理想气体等效方法得到的冲击波峰值压力最低,理想气体等效模拟的峰压相较于前两种等效方法小1～2个数量级;爆炸能量等效与冲击波能量等效的冲击波波速相等,且高于理想气体等效的冲击波波速;沉积能量减小会使得3种等效方法模拟的峰压均有不同程度的减小,但大小顺序不发生变化;改进后的等效爆炸方法能够适应沉积能量的变化,与Touya经验公式拟合较好;基于LS-DYNA对液电效应冲击波峰值压力进行准确模拟,除了选取适合的等效方法,还应结合具体的放电条件,建立适当的数值模型,在满足计算要求的条件下实现冲击波峰值压力的快速计算。     

化爆冲击波在T型通道内到时规律的实验研究

对空气冲击波在T型坑道内传播的到时规律及其影响因素进行了实验研究,在量纲一分析基础上建立了简化的工程模型。通过对实验数据拟合建立了可以对高能炸药在T型坑道内爆炸产生的空气冲击波到达时间进行预估的公式,利用该公式可以估算空气冲击波在T型坑道内的传播速度。     

多级柱面炸药内爆磁通量压缩技术研究

柱面内爆磁通量压缩发生器是利用炸药内爆压缩其内部磁通量至轴线附近小体积内从而实现超高磁场,传统的单级装置因受到金属套筒内爆失稳等影响性能指标受限。开展了多级内爆磁压缩技术研究,突破多项关键技术,包括研制特殊结构的密绕螺线管、脉冲功率源及大电流放电开关等,具备在直径135 mm套筒空间内实现20 T以上初始磁场产生能力,并建立了动态磁光测量系统。利用磁流体力学编码SSS-MHD开展多级装置设计,计算显示,设计的多级装置能够将约42%的初始磁通量压缩至轴线附近直径7 mm的空间内。最终研制成功多级内爆磁压缩装置CJ-150,在亚立方厘米以上空间实现轴向峰值磁场强度906 T,数据不确定度5.35%。10余发动态考核实验显示,CJ-150装置工作稳定,能够满足物理实验需要。利用经实验验证的磁流体模型计算显示,CJ-150具备1 000 T以上超强磁场产生能力,能够对大尺寸样品实现500 GPa以上的准等熵加载。     

内爆压缩多层密绕螺线管的数值模拟

利用AUTODYN二次开发接口建立三维多层密绕螺线管数值模型,并实现周期性边界条件,应用光滑粒子动力学方法对内爆压缩多层密绕螺线管过程及其界面不稳定性开展数值模拟。计算结果表明,内爆压缩螺线管结构过程存在扰动快速增长至后期的界面失稳现象,与对应的实验结果较为相符。同时,计算显示螺线管结构参数对界面不稳定性发展具有显著影响,螺旋角度减小,结构压缩后期的界面不稳定性趋于严重;铜线直径减小,结构压缩后期的界面不稳定性趋于减弱。     

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序，对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为，进行了数值模拟研究。结果表明，首次冲击后至约12 μs，界面发展以RM（Richtmyer-Meshkov）不稳定性为主；12 μs后至冲击波聚心反弹加载前，界面聚心运动处于加速减速状态，界面发展由RT (Rayleigh-Taylor)不稳定性主导；冲击波聚心反弹加载后，界面发展又由RM不稳定性主导。另外，还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示:初始振幅较大时振幅增长也较大；初始波长较小（模数较大）时振幅增长较小，而且存在一个截止波长；钢壳厚度会抑制扰动增长，也存在一个截止厚度；几何汇聚效应会使扰动增长速度更快。     

锥形长药柱水下爆炸冲击波参数计算方法

为了计算锥形长药柱水下爆炸冲击波压力,以及研究长脉宽冲击波的传输特性,基于叠加原理建立了冲击波压力-时间曲线的计算方法,通过实验验证了该方法的有效性,在此基础上分析了锥形长药柱不同方位冲击波压力的分布规律。研究结果表明:锥形长药柱产生的冲击波压力具有各向异性,在起爆端一侧形成的是具有厚波头特征的低幅值长脉宽冲击波,在装药径向形成的是接近指数衰减的高幅值冲击波,而在远离起爆端的冲击波压力幅值和脉宽则介于前两者之间。锥形长药柱与球形装药冲击波分布的差异是由于装药形状和起爆方式的改变所导致的,由于装药不同部位起爆的时间差,导致水下爆炸冲击波在不同位置的叠加效果存在明显差异,药柱周围流场中形成的冲击波压力具有方向性。利用提出的计算方法得到的计算结果与实验结果和数值模拟结果吻合较好,研究结果可为锥形长药柱水下爆炸冲击波威力场和毁伤评估提供参考和依据。     

Simulation of converging cylindrical GPa-range shock waves generated by wire array underwater electrical explosions

Results of one-dimensional numerical simulations of the parameters of the converging strong shock wave generated by electrical underwater explosions of a cylindrical wire array with different array radii and different deposited energies are presented. It was shown that for each wire array radius there exists an optimal duration of the energy deposition into the exploding array, which allows one to maximize the shock wave pressure and temperature in the vicinity of the implosion axis. The simulation results agree well with the 130-GPa pressure in the vicinity of the implosion axis that was recently obtained, which strongly indicates the azimuthal symmetry of the converging shock wave at these extreme conditions. Also, simulations showed that using a pulsed power generator with a stored energy of ~200 kJ, the pressure and temperature at the shock wave front reaches ~220 GPa and 1.7 eV at 0.1 mm from the axis of implosion in the case of a 2.5 mm radius wire array explosion. It was found that, in spite of the complicated equation of state of water, the maximum pressure at the shock wave front at radius r can be estimated as P ≈ (P*(r*/r) ^α , where P* is the known value of pressure at the shock wave front at radius r* ≥ r and α is a parameter that equals 0.62±0.02. A rough estimate of the implosion parameters of the hydrogen target after the interaction with the converging strong shock wave is presented as well.     

Super-spherical cumulation: shock collapse which is more intense than the spherical one

A shock wave implosion in an axisymmetric chamber with a convex bounding wall is studied experimentally, analytically and numerically. The converging shock front area in this geometry shrinks quickly as the shock wave approaches the center point. The analytical theory predicts that the corresponding rate of post-shock pressure and density increase in this case exceeds essentially that achieved in the classical cylindrical or spherical shock implosions, hence, the phenomenon is referred to as “super-spherical cumulation”. The experiments confirm higher intensity of the super-spherical implosion compared with the cylindrical one, both driven by identical high-current pulsed electric discharges. The converging shock stability is analyzed in the framework of the CCW theory. The numerical results obtained using a locally-adaptive unstructured grid technique agree well with the theoretical predictions of the converging shock wave intensity. Received 28 January 1998 / Accepted 6 November 1998     

Computational study of shock wave focusing in a log-spiral duct

The mechanism of shock wave focusing in a two-dimensional log-spiral duct has been investigated here numerically using a finite volume method. This approach is based on a MUSCL TVD scheme with flux-vector splitting applied to the Euler equations. The isopycnics determined from the calculations are compared with the experimental results obtained by use of holographic interferometric photography and are found to be in excellent qualitative agreement with the experiments. The computational results clarify the details of the wave interactions very near to the focus. In particular, phenomena such as the formation of secondary shock waves prior to the implosion, their interaction with the reflected shock and the formation of vortices after the implosion have been examined.     

Experimental and numerical study of shock wave propagation in water generated by pulsed arc electrohydraulic discharges

The objective of this study is to simulate the propagation of the shock wave in water due to an explosion. The study is part of a global research program on the development of an alternative stimulation technique to conventional hydraulic fracturing in tight gas reservoirs aimed at inducing a distributed state of microcracking of rocks instead of localized fracture. We consider the possibility of increasing the permeability of rocks with dynamic blasts. The blast is a shock wave generated in water by pulsed arc electrohydraulic discharges. The amplitude of these shock waves is prescribed by the electrohydraulic discharges which generate high pressures of several kilobars within microseconds. A simplified method has been used to simulate the injected electrical energy as augmentation of enthalpy in water locally. The finite element code EUROPLEXUS is used to perform fluid fast dynamic computation. The predicted pressure is consistent with the experimental results. In addition, shock wave propagation characteristics predicted with simulation can be valuable reference for design of underwater structural elements and engineering of underwater explosion.     

一种水下爆炸冲击波压力调控方法

起爆位置和装药形状对水下爆炸冲击波压力有较为显著的影响,这使得利用小当量装药在局部方向形成与大当量装药一定程度等效的冲击波成为可能。为了能够在小当量装药条件下开展舰船结构及设备抗水下爆炸冲击实验,基于细长装药结构和参数优化设计,设计了一种冲击波压力幅值和持续时间可调的装药方法。首先,基于简单波理论给出了水下爆炸冲击波压力调控的原理,以及装药参数优化设计的目标函数和约束条件;然后,采用自主数值模拟软件研究了细长装药的水下爆炸能量输出规律,通过实验验证了数值模拟的置信度,研究发现起爆位置和装药形状对水下爆炸冲击波压力峰值和持续时间的影响是显著的,在炸药爆速一定的情况下,长药柱水下爆炸冲击波压力的持续时间可通过几何近似确定;最后,为了进一步考察该方法的有效性,以1000 kg TNT和100 m爆距的水下爆炸冲击波压力-时间曲线作为原型,设计了2种与该原型冲击波压力等效的装药方案,并通过数值模拟进行了验证。研究结果表明:设计的装药能够在预定的持续时间内,在装药起爆端一侧形成与原型等效的冲击波压力-时间曲线。由于没有考虑对气泡载荷的等效,因此该调控方法仅适用于中远场爆炸冲击问题。     

The shock–vortex interaction patterns affected by vortex flow regime and vortex models

We have used a third-order essentially non-oscillatory method to obtain numerical shadowgraphs for investigation of shock–vortex interaction patterns. To search different interaction patterns, we have tested two vortex models (the composite vortex model and the Taylor vortex model) and as many as 47 parametric data sets. By shock–vortex interaction, the impinging shock is deformed to a S-shape with leading and lagging parts of the shock. The vortex flow is locally accelerated by the leading shock and locally decelerated by the lagging shock, having a severely elongated vortex core with two vertices. When the leading shock escapes the vortex, implosion effect creates a high pressure in the vertex area where the flow had been most expanded. This compressed region spreads in time with two frontal waves, an induced expansion wave and an induced compression wave. They are subsonic waves when the shock–vortex interaction is weak but become supersonic waves for strong interactions. Under a intermediate interaction, however, an induced shock wave is first developed where flow speed is supersonic but is dissipated where the incoming flow is subsonic. We have identified three different interaction patterns that depend on the vortex flow regime characterized by the shock–vortex interaction.      

直坑道内防护门上的爆炸冲击波荷载

为可靠开展抗爆结构设计与评估,基于理论计算与数值分析对直坑道内爆炸冲击波荷载的计算问题进行了研究。定量对比分析了现有不同方法计算结果的差异,基于结构响应对不同方法进行了评价,并结合模型实验对近距离爆炸情况下直坑道内防护门上的设计荷载取值问题进行了研究。结果表明,在进行荷载简化时有必要考虑荷载形式与结构响应的耦合,而现有坑道内爆炸冲击波荷载的简化计算方法普遍没有考虑结构动态特性影响,且相互之间存在很大差异,严重影响设计或评估工作的可靠性;近距离爆炸情况下,取门中线上距门边1/4宽度处的压力或门上平均压力作为坑道内防护门上的设计荷载,在比较宽的结构频率范围内是合理的。研究结果可为坑道内结构的防护设计提供参考。     

Reflection of a ring shock wave from a rigid wall

A new type of an implosion has been observed experimentally and simulated in a numerical calculation: the formation of a quasi-spherical converging shock wave after the reflection of a ring shock wave from a solid wall. The conversion of the ring shock wave into the quasi-spherical converging shock wave intensifies the local implosion properties.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

聚心火焰与激波相互作用的数值研究

基于带化学反应的2维轴对称Euler方程,利用带有monotonized centered(MC)限制器的波传播算法,在两端敞开的圆桶中对惰性介质的聚心激波和氢气-空气混合物的聚心火焰与激波的相互作用进行了数值模拟。数值结果表明,在惰性介质中激波在轴心的每次汇聚均可成长为马赫干,马赫干的追赶使激波得到一定程度的增强,但整体还呈下降趋势。在氢气-空气混合物中,燃烧诱导的激波,由于与火焰的反复作用,使激波在轴心处产生马赫干的频率和强度皆高于惰性介质中的情形。同时,火焰在与激波的相互作用过程中发生失稳变形,使其形状呈扁平头部的蘑菇云。