Numerical simulation of jet formation and penetration characteristics in multi-point initiation mode

Xin Wang1 , Changxiao Zhao2 , Chong Ji3 , Xinghua Li4 , Mingshou Zhong5 , Huayuan Ma6 , Yun Gu7

1, 2, 3, 4, 5, 6College of Field Engineering, Army Engineering University of PLA, Nanjing, China

7Nuclear Industry Nanjing Construction Group Co., Ltd, Nanjing, China

3Corresponding author

Vibroengineering PROCEDIA, Vol. 37, 2021, p. 86-92. https://doi.org/10.21595/vp.2021.21983
Received 2 April 2021; received in revised form 11 April 2021; accepted 18 April 2021; published 7 May 2021

Copyright © 2021 Xin Wang, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF References
Cite this article
Views 11
Reads 2
Downloads 104
CrossRef Citations 0
Abstract.

In order to explore the application feasibility and effective gain of shaped charge jet under multi-point initiation system, the jet forming characteristics of 40 mm diameter shaped charge under multi-point initiation mode were analyzed, and the influence rules of the number of initiation points and initiation radius on jet parameters were obtained. The simulation results of jet penetrating cylindrical shell covered charge show that the multi-point initiation system can effectively improve the jet tip velocity and the length of the jet, which can effectively improve the impact initiation ability of the shell charge.

Numerical simulation of jet formation and penetration characteristics in multi-point initiation mode

Highlights
  • The jet forming characteristics of 40mm diameter shaped charge under multi-point initiation mode were analyzed.
  • The influence rules of the number of initiation points and initiation radius on jet parameters were obtained.
  • The capability of jet penetrating cylindrical shell covered charge under multi-point initiation mode was evaluated.

Keywords: shaped charge jet, multi-point initiation system, formation, penetration.

1. Introduction

Shaped charge jet (SCJ) is a kind of condensed high-speed penetrator, which has high speed and is widely used in penetrating fortifications, armor targets, destroying unexploded ordnance and so on [1]. In practical application, the diversity of targets and the lightweight of warhead require the shaped charge to reduce its weight and size as much as possible when used to destroy targets, which poses a severe challenge to the design of warhead. By changing the number of initiation points, the multi-point initiation system can increase the detonation pressure in the explosive and also change the collapse angle, thus affecting the stress state of the liner and its forming [2], which becomes a potential choice to improve the characteristics of SCJ. At present, multi-point initiation system is widely used in rod penetrator and explosively formed projectile, but the application is limited due to the sensitivity of shaped charge jet to initiation synchronization [3].

In this paper, according to a typical shaped charge structure of 40 mm caliber, the SCJ forming characteristics under different initiation modes were analyzed by numerical simulation. The influence of the number of initiation points and the initiation radius on the characteristic parameters of jet was obtained, and a better initiation mode is optimized. The numerical simulation of penetrating cylindrical shell covered charge was carried out, so as to provide reference for the application of multi-point initiation system in small cone angle shaped charge, and also a basis for the application of multi-point initiation system in the destruction of unexploded ordnance.

2. Numerical simulation model and material parameters

2.1. Numerical simulation model

The shaped charge used in this paper was composed of main charge, liner and shell. The main charge is 40 mm in diameter, 60 mm in height and 40° in cone angle. The main charge was 8701. The liner was made of copper with a wall thickness of 1.6 mm. The shell was 2 mm thick and made of aluminum alloy. The 8-node polyhedral solid element in LS-DYNA software (Solid164) was used in the calculation model. The air domain, explosive and shaped charge structure adopt Euler element, and the shell and charge in cylindrical shell charge used Lagrange mesh. Thus, the model was calculated by fluid structure coupling method. At the same time, a 1/4 model was established to simplify the calculation. Fig. 1 shows the 1/4 finite element model of shaped charge and its components.

Fig. 1. 1/4 finite element model of shaped charge and its components

 1/4 finite element model of shaped charge and its components

a) Shell

 1/4 finite element model of shaped charge and its components

b) Charge

 1/4 finite element model of shaped charge and its components

c) Liner

 1/4 finite element model of shaped charge and its components

d) Shaped charge

2.2. Material model and parameters

In order to better simulate the phenomena of liner collapsing, forming and perforating when interact with the shell charge, Johnson-cook material model with failure fracture criterion [4] was selected for liner and shell of charge, which takes into account the strain hardening, strain rate and temperature effect of material.

For von Mises yield stress model, the material yield stress is expressed as follows:

(1)
σ y = A + B ε - p n 1 + C l n ε ˙ * 1 - T * m ,

where: ε-p is equivalent plastic strain, ε˙* is relative equivalent plastic strain rate, T* is relative temperature, A is yield stress, B is strain hardening, n is strain hardening index, C is strain rate correlation factor, m is temperature correlation factor.

The expression of the fracture strain is expressed as follows:

(2)
ε f = D 1 + D 2 e x p D 3 σ * 1 + D 4 l n ε ˙ * 1 + D 5 T * ,           D = Δ ε y ε f ,

where D is the damage to a material element, Δεyis the increment of accumulated plastic strain, and εfis the accumulated plastic strain to failure under the current conditions of stress triaxiality, strain rate and temperature. Failure occurs when D= 1. Material parameters selected from reference [5].

The explosive of shaped charge was 8701, HIGH_EXPLOSIVE_BURN material model and JWL equation of state were selected, see in reference [3] for specific parameters. The explosive in the shell charge was TNT and ELASTIC_PLASTIC_HYDRO material model and IGNITION_AND_GROWTH_OF_REACTION_IN_HE were selected for the study on shock initiation characteristics. For air, Null material model and LINEAR_POLYNOMIAL equation of state were used. The material parameters were obtained from reference [5].

3. Analysis of numerical simulation results

3.1. Analysis of jet forming process

The typical forming progress of SCJ was shown in Fig. 2. The initiation mode was center point initiation. From the numerical simulation results, it can be seen that at about 10.4 μs after initiation, most of the liner materials have completed the movement to the symmetrical plane. At this time, the jet after collision redistributes the energy. The jet gradually forms at the head, and the pestle body forms at the tail. Due to the velocity difference between the head and tail, the jet was stretched. The final shape of jet formed at about 34.6 μs after initiation, and the velocity of jet head was about 5500 m/s.

Fig. 2. Shaped charge jet forming process with center point initiation

 Shaped charge jet forming process with center point initiation

a)t= 0 µs

 Shaped charge jet forming process with center point initiation

b)t=  12.6 µs

 Shaped charge jet forming process with center point initiation

c)t= 71 µs

 Shaped charge jet forming process with center point initiation

a)t= 21.2 µs

 Shaped charge jet forming process with center point initiation

d)t= 34.6 µs

Fig. 3 shows the propagation process of detonation front in the case of center point initiation. When t= 0.2 μs, the spherical detonation front grows steadily and propagates to the unexploded explosive. When t= 2.2 μs, the spherical detonation wave reaches the top of the liner and starts to impact and squeeze the liner; the element pressure on the top of the liner suddenly jumps to 20.7 GPa. When t= 3.2 μs, the detonation products attached to the detonation wave continue to squeeze the liner, and the detonation wave begins to impact and squeeze the circumferential shell. When t= 7.6 μs, the peak pressure of the detonation front is 31.6 GPa, and the converged detonation wave moves towards the liner direction. After a series of complex interactions, the new wave system and detonation products will have a new round of impact and extrusion on the liner, and the maximum front pressure is 38.9 GPa.

Fig. 3. The propagation process of detonation front of center point initiation

 The propagation process of detonation front of center point initiation

a)t= 0.2 µs

 The propagation process of detonation front of center point initiation

b)t= 2.2 µs

 The propagation process of detonation front of center point initiation

c)t= 3.2 µs

 The propagation process of detonation front of center point initiation

d)t= 7.6 µs

3.2. Influence of initiation point

Based on the above model, numerical simulation under four initiation modes with 2, 4, 8 and 12 initiation points (N) was carried out respectively. The initiation points were distributed in a ring shape and the initiation radius was 20 mm. Fig. 4 shows the propagation state of detonation wave in charge and the pressure distribution on the upper surface of the liner.

When multi-point initiation was adopted, the detonation wave would converge and collide in charge, and the peak pressure on the upper surface of the liner was 23.68 GPa, 26.30 GPa, 26.74 GPa and 25.51 GPa respectively, which were significantly larger than that of single point initiation of 17.13 GPa. The distribution of detonation wave initiated by two points was asymmetric (Fig. 4(d)). With the increase of the number of initiation points, the pressure distribution on the upper surface of the liner was more uniform, which was closer to the effect of ring initiation.

Fig. 4. The propagation process of detonation front of multi-point initiation

 The propagation process of detonation front of multi-point initiation

a)N= 12

 The propagation process of detonation front of multi-point initiation

b)N= 8

 The propagation process of detonation front of multi-point initiation

c)N= 4

 The propagation process of detonation front of multi-point initiation

d)N= 2

Fig. 5 is the variation of jet tip velocity with initiation point. The initial velocity of jet tip increases with the number of initiation points. Compared with single point initiation, the increment of initial velocity of jet head can reach 980 m/s, and the maximum increase amplitude of jet tip velocity can reach 15.63 %.

Fig. 5. The variation of jet tip velocity with initiation point

 The variation of jet tip velocity with initiation point

3.3. Influence of initiation radius

In order to explore the influence of distribution radius of multiple initiation points on jet performance, the circular distribution radius r of eight initiation points was selected, as shown in Fig. 6. The distribution radius r is 4 mm, 8 mm, 12 mm, 16 mm and 20 mm respectively.

By extracting the pressure and impulse data of element (H5601) on the top of liner and drawing the curve in Fig. 7, it shows that the maximum pressure and impulse of element H 5601 are increasing with the increase of initiation radius. We can get enlightenment from the pressure propagation characteristics in Fig. 8.

As shown in the Fig. 8, the detonation wave generated at each initiation point was superimposed for the first time at 0.8 μs, and then gradually converged to the bottom of the charge and the unexploded area under the reflection of the shell; the detonation wave was superimposed at 2.6 μs to form a high pressure area, with an instantaneous high pressure of 119 GPa; as the high pressure superposition area of the initiation ring propagates downward along the central axis and acts on the charge. The top of the liner was then crushed and shaped.

Fig. 6. Schematic diagram of initiation point position and pressure nephogram of liner (t= 6 μs)

Schematic diagram of initiation point position and pressure nephogram of liner (t= 6 μs)
Schematic diagram of initiation point position and pressure nephogram of liner (t= 6 μs)

Fig. 7. The variation of maximum pressure and impulse at observation point with r (t= 6 μs)

The variation of maximum pressure and impulse at observation point with r (t= 6 μs)
The variation of maximum pressure and impulse at observation point with r (t= 6 μs)

Fig. 8. Peak pressure variation of explosive detonation front under r= 20 mm condition

 Peak pressure variation of explosive detonation front under r= 20 mm condition

Fig. 9. The change of jet tip velocity and jet length

 The change of jet tip velocity  and jet length

As the radius of the initiation ring becomes smaller, that is, the initiation point is closer to the central axis, and the detonation wave does not grow sufficiently at the time of the first superposition, the instantaneous high pressure will decrease, and the overall impulse acting on the liner will also decrease. The change of jet tip velocity and jet length recorded in Fig. 9 can well reflect this rule. Theoretically, the bigger the radius of initiation ring is, the better jet state can be obtained.

4. Impact initiation of jet on cylindrical shell covered charge

In order to study the impact initiation ability of SCJ to the charge with shell, the numerical simulation of SCJ detonating charge with shell under single initiation was carried out. The numerical model was shown in Fig. 10. It is necessary to keep the shape intact and the overall characteristics in a good state before the SCJ impinging into the cylindrical shell. Through the numerical analysis in the last section, standoff distance h=3d was selected, namely 120 mm.

Fig. 10. Schematic diagram of shaped charge impacting into cylindrical shell

Schematic diagram of shaped charge impacting into cylindrical shell

Fig. 11 shows the initiation process of 10 mm thick cylindrical shell covered charge by SCJ. Six observation points (#319005, #319083, #319161, #319239, #319317, #319315) are selected along the axis of the charged charge, shown in the Fig. 12, to record the pressure. The pressure of the charge has exceeded the critical initiation pressure of TNT by 10.4 GPa, it is considered that the charge has been successfully initiated. The shell has not been broken down, which indicates that the strong shock wave produced by the jet impinging on the shell is the shock initiation mechanism.

Fig. 11. Detonation growth process of covered charge

Detonation growth process of covered charge

a)t= 47.6 µs

Detonation growth process of covered charge

b)t= 48.8 µs

Detonation growth process of covered charge

c)t= 51.0 µs

Detonation growth process of covered charge

d)t= 53.8 µs

Fig. 12. Observation point diagram and pressure curve

Observation point diagram and pressure curve
Observation point diagram and pressure curve

5. Conclusions

In this paper, the analysis of SCJ forming and penetrating into shell covered charge under multi-point initiation mode was carried out, the detailed process of jet forming was revealed, and the reliable gain of multi-point initiation mode on penetration of SCJ was obtained. The conclusions are as follows:

1) The tip velocity of jet increases with the number of initiation points. Compared with single point initiation, the maximum increase amplitude of jet tip velocity can reach 15.63 %.

2) The increase of initiation radius is helpful to the full growth of detonation wave before collision, which can effectively improve the jet velocity. Under the eight point initiation mode, the jet tip velocity of full-scale initiation is 14.2 % higher than that of single point initiation.

3) The SCJ formed by single point initiation could initiate 10 mm thick shell covered charge through shock initiation mechanism.

References

  1. Hansenberg Divad Consequences of Coaxial Jet Penetration Performance and Shaped Charge Design Criteraia. Naval Postgraduate School, Monterey, California, 2010. [Search CrossRef]
  2. Li W., Wang X., Li W. The effect of annular multi-point initiation on the formation and penetration of an explosively formed penetrator. International Journal of Impact Engineering, Vol. 37, Issue 4, 2010, p. 414-424. [Publisher]
  3. Liu Jian-Ging, et al. Formation of explosively formed penetrator with fins and its flight characteristics. Defense Technology, Vol. 10, Issue 2, 2014, p. 119-123. [Publisher]
  4. Johnson Gordon R., Cook William H. Fracture characteristics of three metals subjected to various strains, strain rates, temperatures, and pressures. Engineering Fracture Mechanics, Vol. 21, 1985, p. 31-48. [Publisher]
  5. LS-DYNAR Keyword User’s Manual V971. Livermore Software Technology Corporation (LSTC), Livermore, 2012. [Search CrossRef]