圆柱形装药驱动轴向预制破片飞散特性

为了增强柱形战斗部轴向威力,在无壳柱形战斗部底面布置单层离散的预制破片。开展圆柱形TNT装药驱动轴向预制破片飞散试验,获得预制破片的最大初速、飞散角等特征参数;运用LS-DYNA软件对装药驱动预制破片过程进行数值模拟,阐述预制破片群飞散过程;对装药驱动整体平板理论计算公式进行改进,获得预制破片的最大初速。结果表明:破片初速理论计算结果、数值计算结果和试验结果吻合良好;随着与装药底部中心距离的加大,破片初速、径向飞散角分别近似呈“抛物线”减小、增大;试验实测、理论计算得到的破片最大初速值超过2500 m/s,试验实测的径向飞散角最大约为22°,而周向飞散角则普遍较小,均值在5°以内。     

双球破片冲击柱壳装药的临界起爆条件

为了研究导弹战斗部(柱壳装药)在破片场中的累积毁伤问题,在现有单破片起爆平板装药的Jacobs-Roslund经验准则的基础上,分别建立了考虑破片尺寸、破片撞击角度、柱壳装药的装药半径和壳体厚度的单球破片、双球破片冲击柱壳装药临界起爆条件的工程分析模型.该模型计算结果与数值模拟结果和现有试验结果相吻,证明利用该模型能较...     

Parametric study of single confined fragment launch explosive device

In this article, parametric study of single confined fragment launch device was carried out. The configuration proposed was further studied to derive the empirical relationship for effect of fragment size, charge size, confinement thickness on fragment velocity. The simulations were carried out using ANSYS-AUTODYNE explicit solver. Fragment velocities were estimated as a function of different parametric combinations of explosive quantities, charge length to diameter ratio, fragment height to diameter ratio, confinement thickness, fragment material and fragment mass. The data was further converted to charge to metal ratio under fragment and confinement. It was observed that, increase in confinement thickness, charge quantity and decrease in fragment height increases the fragment velocity. It is also noted that, charge to metal mass ratio under fragment significantly affects the fragment velocity. At the end, an empirical relationship for fragment velocity interms of all these parameters was established. Using these relations, two velocities 1831.92 m/s and 2523.9 m/s required for NATO STANAG 4496 IM test were estimated. The design parameters for these velocities are presented. Also, the results estimated using the empirical relationship has been compared with published experimental data. Error in the predicted velocities is within the acceptable range. The empirical relationship proposed will be useful for finalization of design of the fragment launch device.     

Simulation of natural fragmentation of rings cut from warheads

Natural fragmentation of warheads that detonates causes the casing of the warhead to split into various sized fragments through shear or radial fractures depending on the toughness,density,and grain size of the material.The best known formula for the prediction of the size distribution is the Mott formulae,which is further examined by Grady and Kipp by investigating more carefully the statistical most random way of portioning a given area into a number of entities.We examine the fragmentation behavior of radially expanding steel rings cut from a 25 mm warhead by using an in house smooth particle hydrodynamic(SPH) simulation code called REGULUS.Experimental results were compared with numerical results applying varying particle size and stochastic fracture strain.The numerically obtained number of fragments was consistent with experimental results.Increasing expansion velocity of the rings increases the number of fragments.Statistical variation of the material parameters influences the fragment characteristics,especially for low expansion velocities.A least square regression fit to the cumulative number of fragments by applying a generalized Mott distribution shows that the shape parameter is around 4 for the rings,which is in contrast to the Mott distribution with a shape parameter of 1/2.For initially polar distributed particles,we see signs of a bimodal cumulative fragment distribution.Adding statistical variation in material parameters of the fracture model causes the velocity numerical solutions to become less sensitive to changes in resolution for Cartesian distributed particles.     

Fragment spatial distribution of prismatic casing under internal explosive loading

Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control. The fragment spatial distribution of prismatic casings is more complex than that of traditional cylindrical casings. In this study, numerical and experimental investigations into the fragment spatial distribution of a prismatic casing were conducted. A new numerical method, which adds the Lagrangian marker points to the Eulerian grid, was proposed to track the multi-material interfaces and material dynamic fractures. Physical quantity mappings between the Lagrangian marker points and Eulerian grid were achieved by their topological relationship. Thereafter, the fragment spatial distributions of the prismatic casing with different fragment sizes, fragment shapes, and casing geometries were obtained using the numerical method. Moreover, fragment spatial distribution experiments were conducted on the prismatic casing with different fragment sizes and shapes, and the experimental data were compared with the numerical results. The effects of the fragment and casing geometry on the fragment spatial distributions were determined by analyzing the numerical results and experimental data. Finally, a formula including the casing geometry parameters was fitted to predict the fragment spatial distribution of the prismatic casing under internal explosive loading.     

Deformation,fragmentation and acceleration of a controlled fragmentation charge casing

Two different finite element software, LS-DYNA and Impetus, have been evaluated to test their ability to predict the deformation, fragmentation and acceleration of a controlled fragmentation charge casing. The general-purpose program LS-DYNA was used with a multi-material ALE formulation and a mass-preserving erosion criterion coupled to a Johnson-Cook fracture criterion. In the Impetus simulations, a third order Lagrangian element formulation was used for the casing and a node-splitting element erosion treatment coupled to a Cockcroft-Latham failure criterion was used to describe casing fracture. The high-explosive gases were described by a discrete particle formalism.In order to acquire data to validate our computational tools and constitutive models, a series of experiments have been performed using a laboratory charge with an internal grooved casing. In the test series, the charge geometry was fixed except that the groove depth were varied from very shallow to very deep resulting in different deformation patterns, fracture modes and terminal velocities. Various diagnostic tools captured the different stages of the expansion and fragmentation of the casing. A high-speed framing camera depicted the deformation pattern before fragmentation and was used to determine the moment when the casing failed. Three different complementary techniques were used to follow the acceleration of the fragments; a Photon Doppler velocimetry to determine the initial acceleration of the casing, double exposed radiographs to estimate the fragment velocity after break-up and a high-speed video to determine the terminal velocity of the fragment after leaving the fireball. In addition, the fragments were soft recovered in a set of sawdust pit tests and their final shape and weight were measured. A SEM was used to characterise the fracture surfaces and to determine the modus of fracture (tensile or shear failure).Comparisons to experiments show that both software can predict the change in deformation behaviour when the groove depth increases, from tangential necking for shallow grooves to radial punching for deep groves. Both software could also reasonable well predict the acceleration of the fragments, though both overestimates the terminal velocity for the charge with the deepest grooves.     

Experimental studies of explosion energy output with different igniter mass

For study the energy output law of cylindrical charge with shell induced by different input energies, four different black powder masses were selected to ignite the main charge. Fragmentation degree of the shell was qualitatively analyzed by the area of holes on the witness plate and the recovered fragments mass. Through theoretical analysis, established the functional relationship between the average mass of fragments and the relative energy output of warhead, obtained how the relative energy output of charge changed with different initial energy input. The results showed that the change of input energy could lead to obvious variation in fragment characteristics, and could also control the output of charge. When the igniter mass increases from 1.55 g to 5.00 g, the relative energy output of the charge increases by 26.28%. Excessive initial input energy will destroy the shell confine in advance, resulting in a decrease in the relative energy output of charge.     

Numerical study on the case effect of a bomb air explosion

When considering the bomb explosion damage effect, the air shock wave and high-speed fragments of the bomb case are two major threats. In experiments, the air shock wave was studied by the bare ex-plosives superseding the real cased bomb; in contrast, the bomb case influence was ignored to reduce risk. The air explosion simulations of the MK84 warhead with and without the case were conducted. The numerical simulation results showed that the bomb case significantly influenced the shock wave generated by the bomb: the spatial distribution of shock wave in the near field changed, and the peak value of shock wave was reduced. Breakage of the case and kinetic energy of the fragmentation consumed 3 and 38% of the explosion energy, respectively. The increasing factors of the peak over-pressure induced by the bare explosive on the ground and in the air were 1.43-3.04 and 1.37-1.57, respectively. Four typical stages of case breakage were defined. The mass distribution of the fragments follows the Mott distribution. The initial velocity distribution of the fragments agreed well with the Gurney equation.     

Numerical simulation study on penetration performance of depleted Uranium (DU) alloy fragments

Due to its high strength, high density, high hardness and good penetration capabilities, Depleted ura-nium alloys have already shined in armor-piercing projectiles. There should also be a lot of room for improvement in the application of fragment killing elements. Therefore, regarding the performance of the depleted uranium alloy to penetrate the target plate, further investigation is needed to analyze its advantages and disadvantages compared to tungsten alloy. To study the difference in penetration per-formance between depleted uranium alloy and tungsten alloy fragments,firstly, a theoretical analysis of the adiabatic shear sensitivity of DU and tungsten alloys was given from the perspective of material constitutive model. Then, taking the cylindrical fragment penetration target as the research object, the penetration process and velocity characteristics of the steel target plates penetrated by DU alloy frag-ment and tungsten alloy fragment were compared and analyzed, by using finite element software ANSYS/LS-DYNA and Lagrange algorithm. Lastly, the influence of different postures when impacting target and different fragment shapes on the penetration results is carried out in the research. The results show that in the penetration process of the DU and tungsten alloy fragments, the self-sharpening properties of the DU alloy can make the fragment head sharper and the penetrating ability enhance. Under the same conditions, the penetration capability of cylindrical fragment impacting target in vertical posture is better than that in horizontal posture, and the penetration capability of the spherical fragment is slightly better than that of cylindrical fragment.     

多点偏心起爆对破片速度增益的影响

多点偏心起爆后,爆轰波相互作用在定向区域内产生了马赫波超压,使得该区域破片速度产生增益.利用爆轰波反射理论,建立多点偏心起爆后爆轰波相互作用特性的理论模型,并确定定向区内破片速度的增益特性.3种不同起爆方式的仿真计算和静爆试验结果对比表明,定向区内不同起爆方式得到的破片增益程度较为吻合.     

蒙特卡洛剖分投影法的爆震弹破片平均比动能评估

为了评估爆震弹封装壳体破片致死、致伤半径,并削弱其平均比动能,基于蒙特卡洛剖分投影法建立破片平均比动能计算模型。通过LS-DYNA及自编程联合仿真方法对爆震弹自然破片、半预制破片比动能计算模型进行仿真求解,获得了全破片全时域的质量分布、初始速度、垂直靶分布、平均比动能阈值及安全半径等指标。结果表明:在相同装药参数下,半预制破片相比于自然破片,质量分布、初始速度阈值更低且更集中;在小于等于2.5 m范围内半预制破片平均比动能阈值比自然破片更低,大于2.5 m后半预制破片平均比动能阈值比自然破片更高。相比于自然破片,采用半预制破片可显著减小爆震弹破片的致死半径,但并不能明显增大其安全半径。     

Field test and probabilistic analysis of irregular steel debris casualty risks from a person-borne improvised explosive device

Person-borne improvised explosive devices (PBIEDs) are often used in terrorist attacks in Western countries. This study aims to predict the trajectories of PBIED fragments and the subsequent safety risks for people exposed to this hazard. An explosive field test with a typical PBIED composed of a plastic explosive charge and steel nut enhancements was performed to record initial fragment behaviour, including positions, velocity, and trajectory angles. These data were used to predict the full trajectory of PBIED fragments using a probabilistic analysis. In the probabilistic analyses a probability of fatality or serious injury was computed. Based on the results presented, many practical conclusions can be drawn, for instance, regarding safe evacuation distances if a person were exposed to a suspected PBIED.     

破片侵彻多层间隔铝靶的试验及数值模拟研究

对破片侵彻多层间隔铝靶进行了试验和数值模拟研究,获得了破片在不同速度、不同侵彻角时破片的侵彻深度,在此基础上利用拟合的方式建立了计算公式,并将数值模拟结果和计算公式的计算结果进行了对比,结果表明两者吻合较好,所建立的计算公式可靠。     

雷诺数对沟槽减阻特性影响的数值分析

采用雷诺平均N-S方程和RNGk-ε湍流模型计算V型沟槽面的湍流边界层流动和黏性阻力,通过改变来流速度大小和沟槽面布置位置,研究了雷诺数对沟槽减阻特性的影响规律。计算结果表明,来流速度对沟槽减阻率的影响很大,对于一种尺度的V型沟槽,存在着一个具有较好减阻效果的来流速度范围,最大减阻率可达8.6%;沟槽面在沿来流方向上的布置位置对其减阻效果的影响则非常小。     

Real-time calculation of fragment velocity for cylindrical warheads

To simulate explosion fragments, it is necessary to predict many variables such as fragment velocity, size distribution and projection angle. For active protection systems these predictions need to be made very quickly, before the weapon hits the target. Fast predictions also need to be made in real time simulations when the impact of many different computer models need to be assessed. The research presented in this paper focuses on creating a fast and accurate estimate of one of these variables - the initial fragment velocity. The Gurney equation was the first equation to calculate initial fragment velocity. This equation, sometimes with modifications, is still used today where finite element analysis or complex mathematical approaches are considered too computationally expensive. This paper enhances and improves Breech’s two-dimensional Gurney equation using available empirical data and the principals of conservation of momentum and energy. The results are computationally quick, providing improved accuracy for estimating initial fragment velocity. This will allow the developed model to be available for real-time simulation and fast computation, with improved accuracy when compared to existing approaches.     

圆柱壳受不对称径向冲击时的塑性动力屈曲分析

本文发展了Lindberg、Duffey等人的思想,考虑不对称冲击下圆柱壳的塑性动力屈曲问题,旨在研究内凹变形与屈曲的关系。用能量法导出了不对称基本变形的微分方程,近似地将之与对称情形下屈曲扰动控制方程联合,构成一个初值问题,然后进行数值求解。考虑了壳厚变化,放弃了屈曲过程中速度保持不变的假定(对称情形中采用)。本文的研究结果很好地解释了高速冲击时倒塌模式上无屈曲波形,屈曲波形遍布加载范围内且并不明显减弱,局部冲击临界速度高于对称情况形下临界速度等许多实验与已有的分析结果,从而明确了屈曲在塑性动响应中的地位。     

Study of simple plane wave generator with an air-metal barrier

Plane wave generators (PWGs) are used to accelerate flyer plates to high velocities with their generated plane waves, which are widely used in the test of dynamic properties of materials. The traditional PWG is composed of two explosives with different detonation velocities. It is difficult to implement the related fabrication processes and control the generated waves due to its complicated structures. A simple plane wave generator is presented in this paper, which is composed of two identical cylindrical high explosive (HE) charges and an air-metal barrier. A theoretical model was established based on two different paths of the propagation of detonation waves, based on which the size of air-metal barrier was calculated for a given charge. The corresponding numerical simulations were also carried out by AUTODYN-2D^® based on the calculated results, which were used to compare with the theoretical calculations. A detonation wave with a flatness of 0.039 μs within the range of 70-percent diameter of the main charge was obtained through the simulations.     

具有在线测量功能的三工位紧塞具压入机设计研究

分析目前使用的药筒紧塞具压入设备的缺点和不足,研究在线测量药筒实际装药高度的方法,提出在线测量药筒实际装药高度,同时根据装药高度控制紧塞具实际压入深度的设计思想,设计了一种具有在线测量功能的采用并行工作方式的三工位药筒紧塞具压入设备。该设备具有自动化程度高和能够精确控制紧塞具压入深度等特点。     

高速破片撞击充液容器形成液压水锤的试验研究

为研究高速破片撞击充液容器形成的液压水锤效应，设计一套试验装置和测试方法，试验破片撞击充液容器产生空腔的变化过程、液体中不同位置处的压力变化以及前后面板的变形情况。试验发现：破片撞击容器后面板时会出现一圈空化气泡，气泡在后面板内表面从撞击点位置迅速沿径向扩展；液压水锤初始冲击阶段，距离撞击点较近区域在初始冲击波压力脉冲过后会出现一个较大负压，而距离撞击点较远区域不受负压影响；破片撞击速度对容器前面板最大变形影响较小，变形范围随着撞击速度的增大沿撞击点向四周扩展。后面板的最大变形及变形范围都随着破片撞击速度的增大而变大。     

Research on the intermediate phase of 40CrMnSiB steel shell under different heat treatments

In this study, 40CrMnSiB steel cylindrical shells were tempered at 350, 500 and 600 ℃ to study the effect of tempering temperature on the dynamic process of expansion and fracture of the metal shell. A mid-explosion recovery experiment for the metal cylinder under internal explosive loading was designed, and the wreckage of the casings at the intermediate phase was obtained. The effects of different tempering temperatures on the macroscopic and microscopic fracture characteristics of 40CrMnSiB steel were studied. The influence of tempering temperatures on the fracture characteristic parameters of the recovered wreckage were measured and analyzed, including the circumferential divide size, the thick-ness and the number of the circumferential divisions. The results show that as the tempering temper-ature was increased from 350 to 600 ℃, at first, the degree of fragmentation and the fracture characteristic parameters of the recovered wreckage changed significantly and then became essentially consistent. Scanning electron microscopy analysis revealed flow-like structure characteristics caused by adiabatic shear on different fracture surfaces. At the detonation initiation end of the casing, fracturing was formed by tearing along the crack, which existed a distance from the initiation end and propagated along the axis direction. In contrast, the fracturing near the middle position consists of a plurality of radial shear fracture units. The amount of alloy carbide that was precipitated during the tempering process increased continuously with tempering temperature, leading to an increasing number of spherical carbide particles scattered around the fracture surface.