Motion of the detonation products of a point-ignited charge in a cylindrical shell

A numerical solution is obtained for the two-dimensional nonsteady problem of the motion of detonation products from a cylindrical high-explosive charge enclosed in a shell with the initiation of detonation at a central point in the end of the charge. The detonation products propagate in vacuum. The strength of the shell is not considered. A three-term equation of state is used for the detonation products.     

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.     

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.     

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.     

On the motion of flying plate under explosive attack

The one-dimensional problem of the motion of a rigid flying plate under explosive attack has an analytic solution only when the polytropic index of detonation products equals to three. In general, a numerical analysis is required. In this paper, however, by utilizing the weak shock behavior of the reflection shock in the explosive products, and applying the small parameter purterbation method, an analytic, first-order approximate solution is obtained for the problem of flying plate driven by various high explosives with polytropic indices other than but nearly equal to three. Final velocities of flying plate obtained agree very well with numerical results by computers. Thus an analytic formula with two parameters of high explosive (i.e. detonation velocity and polytropic index) for estimation of the velocity of flying plate is established.     

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

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

流固耦合数值方法在点爆发散爆轰驱动平板飞片中的应用

采用Eulerian/Lagrangian耦合方法对炸药与飞片直接接触的爆轰驱动问题进行了数值模拟.其中对炸药和环境气体应用Euler方法、对飞片应用Lagrange方法分别计算,流体与固体接触面采用Ghost Fluid方法和外插技术进行处理,全程模拟了炸药起爆、爆轰波传播、大板变形和运动等过程,着重分析了大板飞片后...     

Theoretical analysis and numerical simulation for the spill procedure of liquid fuel of fuel air explosive with shell

Addressing the problem of spill of liquid fuel of fuel air explosive (FAE) with shell, this paper deduces variation relationship of fracture radius and fracture velocity based on continuity and motion equations of shell using the Taylor fracture criteria of cylindrical shells, and analyzes the correctness of the approximating process in detail. Furthermore, analytical solutions and numerical simulations for the liquid fuel spill process were obtained using Adomian's decomposition method (ADM) and the Runge-Kutta numerical method. Results show that both the numerical solution and the analytical solution based on ADM are closer to the result from experiments and computations from commercial dynamics numerical software.     

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

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

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     

Deformation of rock in the explosion of a cylindrical charge in a coaxial cavity filled with liquid

A study is made of the effect of a water sheath around a cylindrical explosive charge on the parameters of shock waves in rocks. A coupled problem for the detonation-product-water-rock system is examined within the framework of continuum mechanics. A binomial isentropic equation is used for the detonation products, water is modeled by a Tate-type equation, and the rock is modeled by an elastoplastic medium which undergoes dilatation. A numerical solution is obtained by a finite-difference scheme. The dependence of the wave processes in the rock and the dimensions of the crushing zone on the radius of the water sheath is studied with the use of sandstone and siltstone as examples. Institute of Hydromechanics, National Academy of Sciences of Ukraine, Kiev, Ukraine. Translated from Prikladnaya Mekhanika, Vol. 35, No. 3, pp. 83–88, March, 1999.     

Characterization of the Mixing Layer Resulting from the Detonation of Heterogeneous Explosive Charges

A dense, two-phase numerical methodology is used to study the mixing layer developing behind the detonation of a heterogeneous explosive charge, i.e., a charge comprising of a high explosive with metal particles. The filtered Navier–Stokes equations are solved in addition to a sub-grid kinetic energy equation, along with a recently developed Eulerian–Lagrangian formulation to handle dense flow-fields. The mixing layer resulting from the post-detonation phase of the explosion of a nitromethane charge consisting of inert steel particles is of interest in this study. Significant mixing and turbulence effects are observed in the mixing layer, and the rms of the radial velocity component is found to be about 25% higher than that of the azimuthal and zenith velocity components due to the flow being primarily radial. The mean concentration profiles are self-similar in shape at different times, based on a scaling procedure used in the past for a homogeneous explosive charge. The peak rms of concentration profiles are 23–30% in intensity and decrease in magnitude with time. The behavior of concentration gradients in the mixing layer is investigated, and stretching along the radial direction is observed to decrease the concentration gradients along the azimuth and zenith directions faster than the radial direction. The mixing and turbulence effects in the mixing layer subsequent to the detonation of the heterogeneous explosive charge are superior to that of a homogeneous explosive charge containing the same amount of the high explosive, exemplifying the role played by the particles in perturbing the flow-field. The non-linear growth of the mixing layer width starts early for the heterogeneous explosive charge, and the rate is reduced during the implosion phase in comparison with the homogeneous charge. The turbulence intensities in the mixing layer for the heterogeneous explosive charge are found to be nearly independent of the particle size for two different sizes considered in the initial charge. Overall, this study has provided some useful insights on the mixing layer characteristics subsequent to the detonation of heterogeneous explosives, and has also demonstrated the efficacy of the dense, multiphase formulation for such applications.     

炸药爆轰产物驱动不可压缩刚粘塑性柱壳的运动

从质量和动量守恒关系出发 ,结合炸药爆轰产物压力剖面的分析和假设 ,导出了不可压缩刚粘塑性柱壳在柱面散心和聚心爆轰产物推动作用下的运动方程。应用该方程式计算爆轰驱动柱壳的运动与实验结果及动力有限元程序DYNA2D的数值模拟结果十分接近 ,说明散心爆轰波场合波后产物流场可用柱面Taylor波的近似解合理描述 ,有限初始半径聚心爆轰波后产物流场可用平面Taylor波近似描述。     

用光电技术测定爆轰参数

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

爆炸载荷作用下铀气溶胶形成机理研究

针对铀材料在爆炸载荷作用下形成放射性气溶胶的过程,采用光滑粒子流体动力学方法开展了数值模拟和实验研究。通过将颗粒动力学和SPH方法结合,建立了炸药爆轰作用于铀金属壳的数值模拟模型,以铀材料比内能为气溶胶转化判据,获得了铀材料转化为气溶胶的物理过程,得到了在相同爆炸当量下,不同质量铀材料的气溶胶转化效率,并与实验结果进行了对比分析。结果显示,铀材料在爆炸载荷作用下,当其比内能达到1.9 MJ/kg时,即可认为完全转变为气溶胶,对于本文中的爆炸装置结构形式,当炸药质量为铀材料质量的6倍时,转化率超过90%。实验验证了数值模拟结果,表明该方法能够对铀材料的气溶胶转化过程进行准确描述。     

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.     

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

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

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.     

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

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

在两层不同爆速炸药滑移爆轰作用下飞板运动的数值计算研究

本文研究用低爆速炸药复盖高爆速炸药在定滑移爆轰作用下飞板运动的计算方法,并对惯性复盖与炸药复盖的结果作出相互比较。