Numerical solution of the two-dimensional unsteady-state problem of the motion of a shell under the action of the products of an axial detonation

A numerical solution is obtained to the problem of the motion of an incompressible cylindrical shell with a charge of explosive, with excitation of the detonation simultaneously along the whole axis of the charge. The strength of the shell is not taken into consideration. A three-term equation of state is adopted for the products of the detonation. In [1] a numerical solution is obtained to the problem of the one-dimensional motion of a shell with the axial detonation of a charge of explosive.     

Numerical solution of the two-dimensional nonstationary problem of the motion of a shell under the action of detonation products

The problem of the motion of an incompressible cylindrical shell with an explosive charge is solved numerically for the propagation of a plane detonation wave from the end of the charge. The strength of the shell is not taken into account. A three-term equation of state [1] is assumed for the detonation products. A comparison is made with the one-dimensional case.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 76–79, July–August, 1972.The authors thank G. S. Roslyakov and V. M. Paskonov for assistance in the work and for helpful advice.     

Hurling of shells by hollow charges

Results of the numerical solution of the problem of one-dimensional hurling of shells by hollow explosive charges are elucidated. The results of the numerical solution are compared with asymptotic formulas. Numerous domestic and foreign papers have been devoted to the question of hurling shells by explosive charges. A numerical solution of the problem of convergence of a ring to the center under the effect of detonation products is presented in [1–3]. The problem of hurling a shell by a hollow explosive charge with an internal lining is considered in [4]; the solution of the problem of hurling a shell by a hollow explosive charge without the cavity lining is presented in [5] on the basis of the energy-balance equations; however, the complete picture of the processes occurring in the detonation products is not considered.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 161–166, May–June, 1976.     

Explosive projection of liquid from a thick-walled cylindrical container

The process of projection of a layer of an ideal liquid enclosed into a cylindrical elastoplastic shell by products of instantaneous detonation of a high explosive charge is studied numerically in a two-dimensional plane formulation. The processes of shell fracture and liquid exhaustion through the resultant slots are considered. Numerical results are analyzed, and analytical relations for angular distributions of radial velocity and mass of the liquid escaping through the slots are derived. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 56–63, May–June, 2009.     

杆式弹对厚壁壳体装药冲击起爆机制模拟分析

为研究高速杆式弹冲击厚壁壳体装药的起爆机制,运用冲击物理显式欧拉型动力学SPEED软件,开展了不同弹径和弹长的钨合金杆式弹与厚壁壳体Comp-B装药相互作用过程的数值模拟,采用升降法获得弹体起爆装药临界着速及装药起爆位置变化。研究结果表明:弹体起爆装药临界着速随弹径增大而显著降低,随弹长增大呈先降低后平缓变化的规律;弹体以临界着速起爆装药时,存在2种装药起爆机制,即弹体贯穿壳体后的宏观剪切起爆和未贯穿壳体的低速冲击起爆,且其机制随弹体着速在临界着速以上继续提高会发生转变,最终均会转变为高速冲击起爆机制;装药起爆位置均发生在炸药壳体交界面后一定距离处,相同机制下此距离随弹体着速提高而减小。     

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.     

高-低爆速圆板炸药串联爆轰引起平面爆轰波的变凸现象

为分析高-低爆速圆板炸药于空中串联爆轰时平面爆轰波在传播过程中演变成凸面波的现象,分别对平面爆轰波阵面后爆轰产物状态参数的分布、有效药量与柱面装药几何尺寸的关系、爆轰产物的状态方程及强爆轰关系式进行了讨论。并以100 mm50 mm的TNT与RC炸药串联爆轰为例,描述了平面爆轰波演变为凸面爆轰波的过程,预估了爆轰波的平面范围和波形差。预估结果与实验结果基本符合。     

Meshfree particle simulation of the detonation process for high explosives in shaped charge unlined cavity configurations

The numerical simulation of the detonation of a high explosive (HE) is generally not an easy task for traditional grid based methods. Smoothed particle hydrodynamics (SPH) method, as a meshfree, Lagrangian and particle method, provides a very attractive approach in dealing with large deformations and large inhomogeneities in the extremely transient high explosive detonation and later expansion process. This paper presents the application of SPH to simulate and analyze the detonation process of high explosive in shaped charge. A three-dimensional SPH code is developed and applied to simulate the shaped charge detonation process in different scenarios. It is observed that for high explosive in a shaped charge, the detonation produced gaseous products experience strong convergence that forms an extremely high-pressure gas jet. Factors such as different charge cavity shapes and different detonation models lead to quite different behavior of the gas jet convergence and later divergence. Further analyses reveal that a critical value for the charge head length exists. Beyond this critical value, increasing the charge head length will not result in improvement on the gas jet convergence performance.Received: 11 March 2002, Accepted: 9 December 2002, Published online: 28 April 2003     

Transient dynamic behavior of a closed gas-filled shell subjected to internal blast loading

The paper concerns the numerical solution of a problem on the forced vibration of a shell system under a blast load. The mathematical model of the transient dynamic process in the given multicomponent system is based on the use of a system of hydrodynamic equations to describe the motion of the detonation products and air and equations from shell theory to determine the stress-strain state of the shell system. The hydrodynamic equations are integrated in moving grids by the method developed by S. K. Godunov. The shell-theory equations are integrated by the integro-interpolation method. The distribution of the loads on the shells, the distribution of the maximum strains along the generatrix of the shell system, and the transformations of the detonation products at characteristic points are presented. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev, Ukraine. Translated from Prikladnaya Mekhanika, Vol. 34, No. 3, pp. 76–82, March, 1999.     

Gasdynamics of the scattering of the detonation products of a cylindrical explosive charge with inert metal particles

A mathematical model of the scattering of the detonation products of a condensed explosive with inert metal particles when the wave travels along the axis of a cylindrical charge is proposed and numerical calculations are carried out. Detonation product scattering is simulated by a two-phase nonequilibrium axisymmetric jet flow in a supersonic external airstream. A technique for calculating the product gas-suspension behind the detonation wave is developed. The optimal values of the difference scheme parameters ensuring the required calculation accuracy are found. The gasdynamic properties of the process associated with the particles and the lateral spread are investigated.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 101–111, September–October, 1995.     

Effect of the casing on the initial parameters of the underwater explosion of a cylindrical charge of explosive

An experimental determination was made of the initial parameters of shock waves in water with the explosion of cylindrical charges of TNT in casings. It is shown that these parameters depend mainly on the dynamic rigidity of the material and the relative weight of the shell. It is established that during the process of the expansion of the casing of the charge with an explosion in limited volumes of water there is formed a region of extremely rapid expansion, whose boundary can be identified with the boundary between the detonation products and the water after destruction of the casing, coinciding in time with the arrival of the cavitation front.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 165–168, January–February, 1973.     

考虑约束爆炸后产物发生化学反应的约束空间内准静态温度计算

考虑爆炸产物发生化学反应产生的影响,对约束爆炸后约束空间内准静态温度的计算进行了研究。以能量守恒定律为基础,考虑爆炸产物的化学反应动力学过程,推导得到约束空间准静态温度的计算公式和方法,使用C++语言编写了计算程序,并对TNT炸药约束爆炸的情况进行了计算。计算结果表明,对于约束爆炸,爆炸产物发生的化学反应对约束空间内温度的变化有明显影响,且不同的药量体积比条件下,准静态温度的变化趋势不同。研究结果可为更准确的计算约束爆炸后的准静态温度及其他爆炸参数提供有效的方法。     

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

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

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

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

壳装炸药殉爆实验和数值模拟

进行了壳装固黑铝炸药殉爆实验,通过观察残留炸药、壳体和见证板变形,判断被发炸药的爆炸情况,得到了炸药临界殉爆距离。建立了壳装炸药殉爆实验计算模型,采用非线性有限元计算方法,对壳装固黑铝炸药殉爆实验进行了数值模拟。计算中采用预设壳体单元破片方法描述主发炸药壳体破片的形成和破片对被发炸药的撞击起爆,炸药临界殉爆距离的计算结果与实验结果基本一致。主要是主发炸药中部的壳体破片撞击到被发炸药,被发炸药起爆位置也在装药中部。炸药壳体厚度主要影响破片速度和质量、被发炸药的防护性能,从而影响炸药临界殉爆距离。 更多还原     

爆轰波碰撞聚能无网格MPM法数值模拟

爆轰聚能经常被应用于工程爆破切割,为了探索新的聚能爆破形式,本文应用无网格MPM法对圆柱筒装药和导爆索多点起爆形式的爆轰波碰撞聚能问题进行三维数值模拟,采用显式积分算法完成爆轰波形成和汇聚过程以及爆轰压力场分布的数值计算,并将数值模拟结果与实验结果做了对比分析。由此验证了无网格MPM法的准确性,也为炸药聚能爆破提供了新的思路和高效的操作方法。     

爆轰波推动金属球壳的初始态分析

爆轰波推动金属介质是经常涉及到的工程问题。本文应用极曲线方法,分析了爆轰波作用下金属球壳初始运动状态,依赖于爆轰波速度和爆轰波与金属界面夹角。并与Freund,H.U.等所做的实验结果相比较。同时讨论了单层介质和复合层介质初始态的差异。     

水下接触爆炸下沉箱码头毁伤效应

为研究水下接触爆炸下沉箱码头毁伤效应和毁伤机理，通过LS-DYNA有限元软件建立沉箱码头水下接触爆炸模型，进行数值模拟研究，并通过试验验证模型准确性。结果表明:运用有限元方法能够较好地模拟水下接触爆炸作用下沉箱码头的毁伤效应，沉箱码头的破坏过程可分为两个阶段:冲击波阶段，沉箱外墙产生初始破口和环状裂缝；气泡膨胀阶段，爆轰产物从破口涌入仓格加速了仓格的变形和毁伤，仓格顶部变形严重导致码头面板破坏，气泡由于冲出水面提前溃灭，码头毁伤在0.14倍的气泡第一次脉动周期基本停止。对比不同爆炸深度，水域中部接触爆炸下沉箱毁伤最为严重，近水面接触爆炸对码头面板的毁伤作用更强。     

Explosive acceleration of liquid from a cylindrical container using the implosive technique

Tthe acceleration of a liquid from a cylindrical container by the products of instantaneous detonation of a superimposed explosive charge (implosion method) is studies numerically in a two-dimensional formulation. Analytical formulas for the asymptotic velocities of liquid expansion and acceleration of the container body were obtained using numerical calculations in a one-dimensional approximation. The effect of the fracture of the body on the velocity of liquid expansion was studied in a two-dimensional formulation.