不同地应力条件下切缝药包爆破的数值模拟

针对切缝药包定向爆破的特点,考虑岩石介质非均匀性的基础上,把岩石爆破视为爆炸应力波动态作用和爆生气体压力准静态作用的过程,基于损伤力学理论建立岩石爆破的力学模型,并对不同地应力条件下切缝药包爆破的裂纹演化规律进行数值模拟,分析不同地应力条件对切缝药包爆破效果的影响。模拟结果表明:采用切缝药包爆破时,裂纹主要萌生于切缝周边,沿切缝方向扩展,切缝对定向裂纹的控制作用明显;当考虑地应力作用,且最大地应力方向与切缝方向垂直时,不利于定向裂纹的扩展;最大地应力方向与切缝方向平行时,有利于定向裂纹的扩展。裂纹的扩展方向受控于切缝角度和最大地应力方向这2个条件,裂纹扩展规模则受到地应力的限制。     

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

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

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

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

爆炸冲击波对墙体绕射效应的研究

采用试验与数值模拟相结合的方法,对爆炸冲击波作用于防爆墙的载荷与冲击波绕过 防爆墙的规律进行了研究. 采用压力传感器测压技术获得了防爆墙前后不同距离的压力曲线. 数值模拟基于ALE方法与JWL状态方程,试验结果与数值模拟基本相符,给出了冲击波作用 于防爆墙迎爆面上的载荷并研究了冲击波绕过防爆墙详细过程. 通过综合分析,得到了爆炸 冲击波绕流的内在机理为防爆墙的设计与加固提供依据和参考.     

水下成组药包爆炸作用下冰盖动态响应数值模拟

利用大型有限元软件ANSYS/LS-DYNA建立了成组药包水下爆炸冰盖动态响应模型,通过数值计算,得到了冲击波峰值压力变化规律,峰值压力与理论计算结果基本吻合。本文分析了冰盖在成组药包水下爆炸冲击载荷下的应力分布及垂直位移响应特征。结果表明,冰盖迎爆面为压缩破坏,背爆面为拉伸破坏,两者数值均十分接近,其中迎爆面最大压力值可达18.12MPa,冰盖近爆炸点最大垂直位移为1.211cm,两药包连线中点最小位移为0.15cm,达到脆性冰盖形成贯通裂隙的基本条件,从而确定了冰盖在水下成组药包大间距布设条件下,其动态破坏形式是以产生裂隙为主要特征。     

微型药柱有限水域爆炸径向压力衰减规律研究

为了研究药柱水下爆炸冲击波压力的衰减规律, 通过实验研究了不同装药量微型药柱 有限水域爆炸固定距离的压力曲线, 得到装药量与峰值压力间的关系, 通过数值模拟方式验 证了实验的正确性. 在同一装药量情况下研究了药柱水下爆炸径向压力的衰减规律, 通过曲 线拟合, 得到径向峰值压力随爆心距离的衰减曲线和方程.     

切槽爆破中V形切槽产生的力学效应研究

本文详细研究了切槽爆破中V形切槽在爆炸冲击波的动态压力和爆生气体的准静态压力作用下产生的力学效应。在此基础上,进一步阐明了V形切槽对岩石定向破裂的控制机理,提出了切槽爆破中裂缝扩展的力学判据。最后,对所得结论用实验结果进行了验证。     

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

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

水下爆炸冲击波的数值模拟研究

应用商业有限元程序MSC.DYTRAN数值模拟了球形药包在无限水域中爆炸产生的冲击波。通过和经验公式计算结果的比较,证明采用合理的计算参数和有限元模型能够较好地模拟水下爆炸冲击波的传播过程。通过调整水的状态方程参数,达到了提升冲击波峰值应力的效果,从而有效地降低了单元数量和计算时间。     

离心机水下爆炸缩比实验原理及数值研究

对离心机水下爆炸缩比实验方法进行探究,对离心机水下爆炸缩比实验的相似理论进行了推导,通过数值计算分析,探究了原模型、离心机缩比实验及常规缩比实验的冲击波载荷、气泡载荷以及气泡动力学行为。结果表明:常规缩比实验不能对气泡行为及垂直方向的近场载荷进行准确的预报,若要保证远场气泡脉冲峰值误差小于10%,则爆距需大于9.5倍气泡最大半径。而离心机缩比实验能够对原模型进行准确的预报,以小当量装药模拟大当量装药水下爆炸整个物理过程,且冲击波和气泡两个阶段完全相似。同时,水深也可以进行几何缩比,克服了常规缩比方法的缺陷。     

聚奥-9C装药的传爆管殉爆

为了研究冲击波作用下引信传爆装置的响应规律,进行了以主发炸药为RDX-8701、被发装置为聚奥-9C（JO-9C）装药的传爆管（含导爆药柱）的殉爆实验。通过观测残留传爆药、壳体和见证块变形,判断传爆管的爆炸程度,分析了殉爆过程中JO-9C爆轰波的成长历程及传播规律。采用AUTODYN软件建立了殉爆实验有限元模型,计算模型中主要考虑了主发炸药产生的爆炸冲击波对传爆管的冲击响应。基于流固耦合方法,通过调整距离模拟计算得到了传爆管的临界殉爆距离和殉爆安全距离。结果表明,传爆管上端的侧角受到爆炸冲击后产生的爆轰波沿斜下方传播,使传爆药柱完全爆轰,并引起导爆药柱发生殉爆;数值模拟结果显示,JO-9C装药的传爆管临界殉爆距离为5.7 mm,殉爆安全距离为8.8 mm。     

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

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

小不耦合系数装药爆破孔壁压力峰值计算方法

基于轮廓爆破孔壁压力峰值计算方法的相关研究,充分考虑空气冲击波的传播与爆轰产物膨胀的过程,理论分析了小不耦合系数装药爆破过程中空气冲击波与炮孔壁的相互作用,建立了三维空气介质径向不耦合装药单孔爆破有限元模型,研究了工程爆破中常用的多种小不耦合系数装药组合工况下,炸药单点起爆后的炮孔壁压力峰值,并获得了相应工况下的孔壁压力峰值较爆生气体准静态等熵膨胀压力的压力增大倍数。结果表明:小不耦合系数装药爆破过程中,爆轰产物参数会对空气冲击波波后物质参数产生显著影响,揭示了小不耦合系数装药爆破与轮廓爆破在孔壁压力峰值计算方法上的本质差异;柱状装药结构爆轰波沿轴向传播使得空气冲击波撞击炮孔壁时存在叠加效应,孔壁压力峰值也相应增大,通过统计分析不同炸药类型、不同岩石类型工况下压力增大倍数与不耦合系数的关系,发现压力增大倍数随不耦合系数的增大近似呈线性增长;基于理论推导结果及常用爆破孔壁压力峰值计算形式,综合考虑炸药性能、孔壁岩石介质条件、不耦合装药系数对空气冲击波撞击炮孔壁后压力增大倍数的影响,提出了不耦合系数较小时爆破孔壁压力峰值计算方法。     

Investigations into water mitigation using a meshless particle method

It is very difficult for traditional numerical methods to simulate the problems of water mitigation which has been increasingly used to reduce blast effects. This paper studies water mitigation problems by using smoothed particle hydrodynamics (SPH), which is a meshless, Lagrangian method appealing in treating large deformation explosion events with significant inhomogeneities. Numerical verifications considering high explosive detonation and underwater explosion shock waves have demonstrated the effectiveness of the SPH method, the solution procedure and the code. Contact and non-contact water mitigation simulations have been carried out and are compared with the case without mitigation. For either contact or non-contact water shield, the peak shock pressure and the equilibrium gas pressure are reduced to different levels according to the relevant geometry of the system setup. An optimum contact water shield thickness is found to produce the best mitigation effect for a given high explosive charge, while the non-contact water shield, if properly designed, can result in further reduction of the peak shock pressure and the equilibrium gas pressure.     

Breakdown of a discontinuity,taking account of the real equation of state of the detonation products from a condensed explosive

The problem is considered concerning the breakdown of an arbitrary discontinuity, due to the emergence of a detonation wave at the boundary of a condensed explosive charge. The real equation of state of the detonation products of Hexogen was used in the numerical calculations. A u vs p diagram is constructed, which allows graphical calculations to be carried out of the discontinuity breakdown for different media. A comparison is carried out of the calculated values of the initial shock-wave velocities with the experimental data obtained at a certain distance from the explosive charge. It is shown that an increase of the pressure of the gas in which dispersion of the explosion products occurs leads to a reduction of the initial shock wave velocity and to an enhancement of its attenuation during its further motion in the shock tube.     

Blast wave propagation in air from a gaseous charge

In the point explosion problem it is assumed that an instantaneous release of finite energy causing shock wave propagation in the ambient gas occurs at a space point. The results of the solution of the problem of such blasts are contained in [1–4]. This point model is applied for the determination of shock wave parameters when the initial pressure in a sphere of finite radius exceeds the ambient air pressure by 2–3 orders of magnitude. The possibility of such a flow simulation at a certain distance from the charge is shown in papers [4, 5] as applied to the blast of a charge of condensed explosive and in [6, 8] as applied to the expansion of a finite volume of strongly compressed hot gas. In certain practical problems the initial pressure in a volume of finite dimensions exceeds atmospheric pressure by a factor 10–15 only. Such cases arise, for example, in the detonation of gaseous fuel-air mixtures. The present paper considers the problem of shock wave propagation in air, caused by explosion of gaseous charge of spherical or cylindrical shape. A numerical solution is obtained in a range of values of the specific energy of the charge characteristic for fuel-air detonation mixtures by means of the method of characteristics without secondary shock wave separation. The influence of the initial conditions of the gas charge explosion (specific energy, nature of initiation, and others) is investigated and compared with the point case with respect to the pressure difference across the shock wave and the positive overpressure pulse.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 110–118, May–June, 1986.     

PBX-01炸药水中爆轰产物状态方程研究

提出通过水中实验确定炸药的水中爆轰产物JWL状态方程参数的方法;选择PBX-01高能炸药进行水中实验,利用ANSYS/LS-DYNA程序建立炸药的水中实验模型,将实验结果与数值计算结果进行对比,确定PBX-01炸药水中爆轰产物的JWL状态方程参数。研究结果显示,圆筒实验确定的JWL参数在反映炸药水中爆轰产物的膨胀状态时有所不足,水中实验确定的JWL状态方程参数能够更准确地描述PBX-01炸药水中爆轰产物的膨胀过程,因此对水中爆炸的研究需要通过水中爆炸实验建一套状态方程参数。     

用光电技术测定爆轰参数

本文介绍了一种用光电技术同时测定爆速和爆温、爆压和爆温的方法。用光电比色技术测定炸药爆轰温度的同时,用两个爆轰面上爆轰的时间差计算爆轰波速度;或用透明介质中冲击波速度来反算炸药中的爆轰压力。方法原理可靠,技术简便。     

Investigation on charge parameters of underwater contact explosion based on axisymmetric SPH method

This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.