Abstract

To study the influence of different initiation modes on the forming characteristics of the MEFP warhead, numerical simulations were carried out on three types of initiation modes. The numerical simulation results showed that the number of EFPs was the least by double-column multipoint synchronous initiation, the number of EFPs was the largest by the central single-point (multipoint) initiation, and single-column single-point (multipoint synchronous) detonation forms the number of EFPs between central single-point (multipoint synchronous) detonation and double-column multipoint synchronous detonation. For the MEFP warhead of a small-caliber grenade, whether it is center detonation or eccentric detonation, the EFP velocity of multipoint detonation is higher than that of the single-point detonation, the velocity of double-column multipoint eccentric synchronous detonation is 2%–9% higher than that of the single-column single-point (multipoint eccentric synchronous) detonation, the velocity of double-column multipoint eccentric synchronous detonation is 10%–17% higher than that of the central single-point (multipoint synchronous) detonation, and the velocity of single-column single-point (multipoint eccentric synchronous) detonation is 5%–17% higher than that of the central single-point (multipoint synchronous) detonation. Research results show that although the number of EFPs is reduced during eccentric single-point (multipoint simultaneous) detonation of MEFP warheads, a higher velocity can be obtained.

1. Introduction

As a kind of shaped charge projectile, explosively formed projectile (EFP) began to appear in the 1970s. In order to improve the hit rate and damage probability of the projectile, the multiple explosively formed projectile (MEFP) began to appear in the 1980s [1]. MEFP is a highly effective damage warhead developed on the basis of a single EFP warhead [2]. The emergence of EFP warhead technology is to deal with armored targets. The MEFP warhead technology has become one of the effective means of destroying armed personnel with the increase in the protection capabilities of armed personnel. The main factors affecting the forming characteristics of EFP include the shape of liner, liner material, liner thickness, charge length diameter ratio, and initiation model [3]. The current research on the forming characteristics of MEFP warheads mainly focuses on two aspects. First, the axial MEFP warhead, that is, the liner, is arranged along the charge axis [410]; second, the circumferential MEFP warhead, that is, the liner, is arranged along the circumference of the charge [1116]. In literature [415], the research on MEFP mainly focused on the material of the liner, the parameters of the liner, the parameters of the charge, and so on. In the study of initiation mode on the forming characteristics of EFP [1721], the main focus is on the influence of single EFP forming characteristics. In the study of initiation mode on the forming characteristics of MEFP warheads, literature [7] used LS-DYNA simulation software to conduct numerical simulation studies on three initiation modes, namely, single-point initiation, ring initiation, and plane initiation; the research results show the use of plane initiation. In this way, the effective utilization rate of explosives is the highest, and the aerodynamic performance of the projectile and the penetration effect on the target are the best. In [16], the damage effectiveness of circumferential MEFP warhead with eccentric initiation was studied by combining experiment and numerical simulation; static detonation tests were carried out for center initiation and eccentric initiation, respectively. It was found that the eccentric initiation model can effectively improve the average velocity, distribution density, and penetration ability of EFP damage elements. Through the analysis of [421], it can be seen that there are more researches on axial MEFP warheads by researchers at present and relatively few researches on circumferential MEFP. There are more researches on the influence of the detonation model on the forming characteristics of a single EFP; there are few studies on the influence of MEFP forming characteristics. Only [7] and [16] study the influence of detonation methods on the forming characteristics of MEFP warheads. Among them, [7] studies the effects of the detonation model on the forming characteristics of axial MEFP warheads. Literature [16] studies the influence of the initiation model on the forming characteristics of the circumferential MEFP warhead. The combination of multipoint initiation technology and MEFP warhead can not only improve the forming velocity of local EFP but also realize the directional damage function of MEFP warhead.

In order to study the influence of multipoint initiation on the forming characteristics of circumferential MEFP warhead of small-caliber grenade, a total of 12 initiation modes in three categories were designed. The first category is central initiation, including single-point central initiation, two-point central synchronous initiation in the middle, two-point central synchronous initiation near the charge end face, and four-point central synchronous initiation. The second category is single-column eccentric initiation, including single-point eccentric initiation, two-point eccentric synchronous detonation in the middle, two-point eccentric synchronous detonation near the two end faces of the charge, and four-point eccentric synchronous detonation. The third type is double-column eccentric initiation, including two-point eccentric synchronous initiation, central four-point eccentric synchronous initiation, four-point eccentric synchronous initiation near the two end faces of the charge, and eight-point eccentric synchronous initiation. Using ANSYS/ICEM, HyperMesh, and LS-DYNA cosimulation method, the designed MEFP warhead was simulated according to 12 detonation methods, and the numerical simulations were compared and analyzed.

2. Structural Design

2.1. MEFP Warhead Structure Design

The MEFP warhead charge diameter of the small-caliber grenade designed by this research institute is 32 mm; the charge height is 52 mm, the liner diameter is 9 mm, the inner wall curvature radius of liner is 12 mm, the inner wall curvature radius of the liner is 9 mm, the thickness of the center of liner is 1.15 mm, and the thickness of the shell is 1.5 mm. There are 12 liners in each row along the radial direction of the charge and four in each column along the charge circumferential direction. There are 48 liners in the entire MEFP warhead structure. The MEFP warhead structure is shown in Figure 1.

In Figure 1, for the convenience of follow-up research, one column of the MEFP warhead is marked as 1#EFP, 2#EFP, 3#EFP, and 4#EFP.

2.2. Initiation Model Design of MEFP Warhead

In order to study the influence of the initiation model on the forming characteristics of MEFP warheads, the following initiation modes were designed: center initiation model and eccentric initiation mode. Among them, eccentric initiation mode includes single-column eccentric initiation and double-column eccentric initiation.

The central detonation mode includes single-point detonation, two-point synchronous detonation, and four-point synchronous detonation. The structure diagram of the designed detonation point is shown in Figure 2.

In Figure 2, (a) represents a single-point detonation, which is represented by center initiation mode 1, (b) represents central two-point center synchronous initiation, which is represented by center initiation mode 2, and (c) represents a two-point synchronous detonation near the end surface, which is represented by center initiation mode 3, and (d) represents four-point synchronous initiation, which is represented by center initiation mode 4.

Single-column eccentric initiation includes single-point eccentric initiation, two-point eccentric synchronous initiation, and four-point eccentric synchronous initiation. Two-point eccentric synchronous initiation includes two-point eccentric synchronous initiation in the middle and two-point eccentric synchronous initiation near the end surface. The structure diagram of the designed initiation mode is shown in Figure 3.

In Figure 3, single-column mode 1 represents a single-column single-point eccentric initiation, as shown in Figure 3(a). Single-column mode 2 represents single-column two-point eccentric synchronous initiation in the middle of charge, as shown in Figure 3(b). Single-column mode 3 represents a single-column two-point eccentric synchronous detonations near the end face of charge, as shown in Figure 3(c). Single-column mode 4 represents a single-column four-point eccentric synchronous detonation, as shown in Figure 3(d).

Double-column eccentric initiation includes two-point eccentric synchronous initiation, four-point eccentric synchronous initiation, and eight-point eccentric synchronous initiation. Among them, four-point eccentric synchronous initiation includes four points in the middle of charge (two columns and two points) and four-point eccentric synchronous initiation near the end of charge (two columns and two points). The structure diagram of the designed initiation mode is shown in Figure 3. The angle between the initiation points of the two rows is 45°.

In the central initiation structure diagram shown in Figure 4, double-column mode 1 indicates double-column two-point eccentric synchronous initiation, as shown in Figure 4(a); double-column mode 2 refers to double-column four-point eccentric synchronous initiation in the middle of charge, as shown in Figure 4(b); double-column mode 3 refers to double-column four-point eccentric synchronous initiation near the end of charge, as shown in Figure 4(c); double-column mode 4 refers to double-column eight-point eccentric synchronous initiation, as shown in Figure 4(d).

3. Numerical Calculation Model

3.1. Finite Element Model

In this study, the joint simulation method combines the advantages of ANSYS/ICEM, HyperMesh, and LS-DYNA. First, ANSYS/ICEM is used to divide the designed MEFP warhead structure into a structured mesh, and then the structured mesh is transformed into unstructured meshes and export mesh files. HyperMesh will continue the preprocessing of finite element analysis required by LS-DYNA, and finally, LS-DYNA will perform numerical calculation and postprocessing. For the central detonation, because the MEFP warhead designed by this research is a centrally symmetric structure, in order to reduce the calculation amount, a quarter finite element model of the MEFP warhead is established, the finite element model is shown in Figure 5(a). For eccentric initiation, two eccentric initiation modes are designed, namely, single-column eccentric initiation and double-column eccentric initiation. In order to better simulate the detonation wave superimposition effect of single-column eccentric initiation, the established half of the MEFP warhead model is shown in Figure 5(b); in order to better simulate the detonation wave superposition effect of double-column eccentric initiation, the complete finite element model established is shown in Figure 5(c).

3.2. Material Model

The material model is the same as that in Section 3.2 in [22].

4. Numerical Calculation Results

4.1. Numerical Simulation Results of the Center Detonation of the MEFP Warhead
4.1.1. Forming Results

Figure 6 shows the molding result when the MEFP warhead is detonated in the center.

Figure 6(a) shows the forming result of MEFP warhead with single-point initiation, and the setting of initiation point is shown in Figure 2(a); in the case of single-point center initiation, the position of detonation wave front arriving at the liner is different, which results in a certain divergence angle of EFPs except 1#EFP. Figure 6(b) shows the forming result of two-point center synchronous initiation in the middle of charge, the initiation point is set as shown in Figure 2(b), and the difference in the positions of the explosion wave front reaching the 1#liner and 4# liner causes 1#EFP and 4#EFP to have a certain dispersion angle. The detonation wave is superimposed in the middle of the two initiation points, and the superimposed detonation wave front acts on the center of the 2#liner and 3#liner, resulting in horizontal dispersion of 2#EFP and 3#EFP along the radial direction of the charge. Figure 6(c) shows two-point center synchronous initiation near the end face of charge, and the initiation point is set as shown in Figure 2(c). During two-point center synchronous initiation near the end face of charge, the detonation wave front acts on the center of the 1#liner and the 4#liner, respectively, and the 1#EFP and 4#EFP disperse horizontally along the radial direction of the charge. As the charge detonation progresses, the detonation front first acts on the side of the 2#liner and 3#liner; when the two detonation wave fronts are superimposed, the detonation wave pressure acts on the other side of 2#liner and 3#liner, causing 2#EFP and 3#EFP to disperse horizontally along the radial direction of the charge. Figure 6(d) shows four-point center synchronous initiation, and the setting of initiation points is shown in Figure 2(d). During the four-point center synchronous initiation, the detonation wave fronts generated by the four detonation points will act on the center of 1#liner, 2# liner, 3#liner, and 4# liner, respectively, after being superimposed, resulting in horizontal dispersion of 1#EFP, 2#EFP, 3#EFP, and 4#EFP along the radial direction of the charge.

4.1.2. Forming Velocity

The velocity contrast curves of four different initiation modes are shown in Figure 7.

In Figure 7, 1i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in central initiation mode 1 (e.g., 11 is the velocity of 1#EFP). 2i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in central initiation mode 2 (e.g., 21 is the velocity of 1#EFP). 3i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in central initiation mode 3 (e.g., 31 is the velocity of 1#EFP). 4i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in central initiation mode 4 (e.g., 41 is the velocity of 1#EFP). Figure 7(a) shows the comparison curve of forming velocity of 1#EFP in four initiation modes; Figure 7(b) shows the comparison curve of forming velocity of 2#EFP in four initiation modes; Figure 7(c) shows the comparison curve of forming velocity of 3#EFP in four initiation modes; Figure 7(d) shows the forming velocity comparison curve of 4#EFP in the four initiation modes. As can be seen from Figure 7, 11 = 2036 m/s, 12 = 2051 m/s, 13 = 2096 m/s, and 14 = 2142 m/s; 21 = 2008 m/s, 22 = 2249 m/s, 23 = 2247 m/s, and 24 = 1983 m/s; 31 = 2117 m/s, 32 = 2429 m/s, 33 = 2389 m/s, and 34 = 2075 m/s; 41 = 2170 m/s, 42 = 2375 m/s, 43 = 2376 m/s, and 44 = 2124 m/s. In Figure 7(a), the velocity of 1#EFP in the four initiation modes is 41 > 31 > 11 > 21. In Figure 7(b), the velocity of 2#EFP in the four initiation modes is 32 > 42 > 22 > 12. In Figure 7(c), the velocity of 3#EFP in the four initiation modes is 33 > 43 > 23 > 13. In Figure 7(d), the velocity of 4#EFP in the four initiation modes is 14 > 44 > 34 > 24. Based on the above analysis, it can be seen that when the small-caliber MEFP warhead adopts a four-point center synchronous initiation or a two-point center close to the end face of charge, the EFPs can obtain a greater velocity. When the small-caliber MEFP warhead uses a single-point center initiation or a two-point center synchronous initiation in the middle of the charge, the EFP velocity is lower.

4.2. Numerical Simulation Results of the Single-Column Eccentric Initiation of the MEFP Warhead
4.2.1. Forming Results

The forming results of single-column single-point (multipoints) eccentric synchronous initiation of small-caliber grenade MEFP warhead are shown in Figure 8, and the detonation pressure nephogram of eccentric initiation is shown in Figure 9.

In Figure 8, 1#EFP-4#EFP represents the EFPs on the side far from the initiation point, and 5#EFP-8#EFP represents the EFPs on the side close to the initiation point. Figure 8(a) shows single-column single-point eccentric initiation, and the pressure nephogram is shown in Figure 9(a). Figure 8(a) shows single-column two-point eccentric synchronous initiation in the middle of charge, and the pressure nephogram is shown in Figure 9(b). Figure 8(c) shows single-column two-point eccentric synchronous initiation near the end face of charge, and the pressure nephogram is shown in Figure 9(b). It can be seen from Figure 8 that the liner on the side close to the initiation point does not form a better-shaped EFPs, and the liner on the side far from the initiation point can form EFPs with the head turned over. The setting of initiation position for single-column single-point initiation of MEFP warhead is shown in Figure 3(a); due to the different positions of detonation wave front arriving at the liner, there is a certain divergence angle between 2#EFP and 4#EFP. The setting of initiation position for single-column two-point eccentric synchronous initiation in the middle of charge of MEFP warhead is shown in Figure 3(a); although the initiation points are located on the horizontal line where the centers of the 2#liner and 3#liner are located respectively, due to the superposition of detonation waves, as shown in Figure 9(b), the detonation wave front is in contact near the center of the 2#liner and the 3# liner; as a result, 2#EFP and 3#EFP formed a certain divergence angle. The setting of initiation position for single-column two-point eccentric synchronous initiation near the end face of charge of MEFP warhead is shown in Figure 3(c); in the initial stage of initiation, the detonation wave front acts on the center of the 1#EFP and 4#EFP, respectively, so that the 1#EFP and 4#EFP scatter along the radial horizontal direction and the maximum pressure after detonation wave superposition is in the middle of 2#EFP and 3#EFP, as shown in Figure 9(c), resulting in 2#EFP and 3#EFP flying away along a radial direction at a certain divergence angle. As the detonation progresses, the superimposed detonation wave pressure diffuses from the middle of 2#EFP and 3#EFP to 1#EFP and 4#EFP, respectively, which results in the formation of certain dispersion angle of 1#EFP and 4#EFP along the radial direction of the charge. The setting of the initiation position for a single-column four-point eccentric synchronous initiation is shown in Figure 3(d). Due to the comprehensive action of detonation wave generated by initiation point on the horizontal line of center of 4#liner and detonation wave generated by initiation point on the horizontal line with center of 2#liner and 3#liner, this results in the fact that 3#EFP disperses horizontally along the radial direction of the charge, while the 4#EFP has a certain dispersion angle along the radial direction of the charge. Due to the comprehensive action of detonation wave generated by initiation point on the horizontal line of center of 1#liner and detonation wave generated by initiation point on the horizontal line with center of 2#liner and 3#liner, this results in the fact that 2#EFP disperses horizontally along the radial direction of the charge, while the 1#EFP has a certain dispersion angle along the radial direction of the charge. The pressure nephogram of detonation wave of single-column four-point eccentric synchronous initiation is shown in Figure 9(d).

4.2.2. Comparison Diagrams of Forming Velocity

The forming velocity comparison diagram of small-caliber MEFP warhead in single-column eccentric initiation modes is shown in Figure 9.

In Figure 10, k1i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in single-column mode 1 (e.g., k11 is the velocity of 1#EFP). k2i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in single-column mode 2 (e.g., k21 is the velocity of 1#EFP). k3i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in single-column mode 3 (e.g., k31 is the velocity of 1#EFP). k4i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in single-column mode 4 (e.g., k41 is the velocity of 1#EFP). Figure 10(a) shows the comparison curve of forming velocity of 1#EFP in four initiation modes; Figure 10(b) shows the comparison curve of forming velocity of 2#EFP in four initiation modes; Figure 10(c) shows the comparison curve of forming velocity of 3#EFP in four initiation modes; Figure 10(d) shows the forming velocity comparison curve of 4#EFP in the four initiation modes. As can be seen from Figure 10, k11 = 2192 m/s, k12 = 2338 m/s, k13 = 2298 m/s, and k14 = 2251 m/s; k21 = 2311 m/s, k22 = 2492 m/s, k23 = 2498 m/s, and k24 = 2278 m/s; k31 = 2236 m/s, k32 = 2562 m/s, k33 = 2578 m/s, and k34 = 2225 m/s; k41 = 2362 m/s, k42 = 2606 m/s, k43 = 2602 m/s, and k44 = 2320 m/s. In Figure 10(a), the velocity of 1#EFP in the four initiation modes is k41 > k21 > k31 > k11. In Figure 10(b), the velocity of 2#EFP in the four initiation modes is k42 > k22 > k32 > k12. In Figure 10(c), the velocity of 3#EFP in the four initiation modes is k43 > k33 > k23 > k13. In Figure 10(d), the velocity of 4#EFP in the four initiation modes is k44 > k24 > k34 > k14. Based on the above analysis, it can be seen that, during a single-column single-point (multipoints) eccentric synchronous initiation of MEFP warheads, the four-point eccentric synchronous initiation has the highest velocity, and the single-column single-point eccentric synchronous initiation has the lowest velocity.

4.3. Numerical Simulation Results of the Double-Column Eccentric Initiation of the MEFP Warhead
4.3.1. Forming Results

The forming result of the small-caliber grenade MEFP warhead with double-column single-point (multipoints) eccentric synchronous detonation is shown in Figure 11.

Figure 11(e) shows the top view of Figures 11(a), Figure 11(f) shows the top view of Figure 11(b), Figure 11(g) shows the top view of Figures 11(c), and Figure 11(f) shows the top view of Figure 11(d). It can be seen from Figure 11 that, during the double-column multipoints eccentric synchronously initiation the small-caliber grenade MEFP warhead, 5 of the 12 columns designed liner can form EFP, which is marked as 1, 2, 3, 4, and 5 in Figures 11(e)11(h)). Due to the superposition of detonation waves, the velocity of the EFP in the second or third column is the highest. By comparing the velocity of the EFPs in the second and third columns, the speed of the EFP in the third column is higher, and the EFP in the third column is marked as 1#EFP, 2#EFP, 3#EFP, and 4#EFP. It can be seen from Figure 11(a) that, except for 1#EFP, 2#EFP, 3#EFP, and 4#EFP have certain divergence angles. In Figure 11(b), since the initiation point is set at four points in the middle of charge, the maximum pressure after detonation wave superposition is located in the middle of 2#EFP and 3#EFP; therefore, the scattering directions of 1#EFP and 2#EFP deviate from the charge radial level upward, and the scattering directions of 3#EFP and 4#EFP deviate from the charge radial level downward. In Figure 11(c), since the initiation point is set near the end face of charge, the maximum pressure after detonation wave superposition is located in the middle of 2#EFP and 3#EFP; therefore, the scattering directions of 1#EFP and 2#EFP deviate from the charge radial level upward, and the scattering directions of 3#EFP and 4#EFP deviate from the charge radial level downward. In Figure 11(d), due to the interaction of the detonation wave, 2#EFP and 3#EFP are scattered along the horizontal direction of the charge, 1#EFP deviates from the charge radial level upward, and the scattering directions of 4#EFP deviate from the charge radial level downward.

4.3.2. Forming Velocity

The forming result of the small-caliber grenade MEFP warhead double-column multipoints eccentric synchronous initiation is shown in Figure 12.

In Figure 12, (m)1i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in double-column mode 1 (e.g., m11 is the velocity of 1#EFP). m2i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in double-column mode 2 (e.g., m21 is the velocity of 1#EFP). m3i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in single-column mode 3 (e.g., m31 is the velocity of 1#EFP). m4i (i = 1, 2, 3, 4) is used to represent the forming velocity of EFP in double-column mode 4 (e.g., m41 is the velocity of 1#EFP). Figure 12(a) shows the comparison curve of forming velocity of 1#EFP in four initiation modes; Figure 12(b) shows the comparison curve of forming velocity of 2#EFP in four initiation modes; Figure 12(c) shows the comparison curve of forming velocity of 3#EFP in four initiation modes; Figure 12(d) shows the forming velocity comparison curve of 4#EFP in the four initiation modes. As can be seen from Figure 12, m11 = 2346 m/s, m12 = 2471 m/s, m13 = 2456 m/s, and m14 = 2381 m/s; m21 = 2368 m/s, m22 = 2623 m/s, m23 = 2608 m/s, and m24 = 2392 m/s; m31 = 2418 m/s, m32 = 2789 m/s, m33 = 2771 m/s, and m34 = 2423 m/s; m41 = 2487 m/s, m42 = 2814 m/s, m43 = 2807 m/s, and m44 = 2506 m/s. In Figure 12(a), the velocity of 1#EFP in the four initiation modes is m41 > m31 > m21 > m11. In Figure 10(b), the velocity of 2#EFP in the four initiation modes is m42 > m32 > m22 > m12. In Figure 10(c), the velocity of 3#EFP in the four initiation modes is m43 > m33 > m23 > m13. In Figure 10(d), the velocity of 4#EFP in the four initiation modes is m44 > m34 > m24 > m14. Based on the above analysis, it can be seen that, during a double-column multipoints eccentric synchronous initiation of MEFP warheads, the double-column eight-point eccentric synchronous initiation has the highest velocity, and the double-column two-point eccentric initiation has the lowest velocity.

4.4. Velocity Comparison of Different Initiation Modes

In order to more clearly analyze the velocity difference between the three types of initiation modes, the 12 different initiation modes under the three types of initiation modes are divided into four groups. The first group is single-point central initiation, single-column single-point eccentric initiation, and double-column two-point synchronous initiation. The second group is the two-point center synchronous detonation in the middle of charge, the single-column two-point eccentric synchronous initiation in the middle of charge, and the double-column four-point eccentric synchronous initiation in the middle of charge. The third group is the two-point center synchronous initiation near the end face of charge, the single-column two-point eccentric synchronous initiation near the end face of charge, and the double-column four-point eccentric synchronous initiation near the end face of charge. The fourth group is central four-point synchronous detonation, single-column four-point eccentric synchronous initiation, and double row eight-point eccentric synchronous initiation.

4.4.1. The Velocity Comparison of the First Group

The velocity comparison diagrams of the first group are shown in Figure 13.

In the velocity comparison diagrams of the first group, the velocity of 1#EFP is m11 > k11 > , which is (m11 = 2346 m/s) > (k11 = 2192 m/s) > ( = 2036 m/s). The velocity of the double-column model 1 is 6.6% higher than that of the single-column model 1, the velocity of double-column model 2 is 13.2% higher than that center initiation model 1, and the velocity of the single-column model is 7.1% higher than the velocity of the center initiation model 1. The velocity of 2#EFP is m12 > k12 > , which is (m12 = 2471 m/s) > (k12 = 2338 m/s) > ( = 2051 m/s). The velocity of the double-column model 1 is 5.4% higher than that of the single-column model 1, the velocity of the double-column model 1 is 17% higher than that of the center initiation model 1, and the velocity of the single-column model 1 is 12.3% higher than the velocity of the center initiation model 1. The velocity of 3#EFP is m13 > k13 > , which is (m13 = 2456 m/s) > (k13 = 2298 m/s) > ( = 2096 m/s). The velocity of the double-column model 1 is 6.4% higher than that of the single-column model 1, the velocity of the double-column model 1 is 14.7% higher than that of the center initiation model 1, and the velocity of the single-column model 1 is 8.8% higher than the velocity of the center initiation model 1. The velocity of 4#EFP is m14 > k14 > , which is (m14 = 2381 m/s) > (k14 = 2251 m/s) > ( = 2142 m/s). The velocity of the double-column model 1 is 5.5% higher than that of the single-column model 1, the velocity of the double-column model 1 is 10% higher than that of the center initiation model 1, and the velocity of the single-column model 1 is 4.8% higher than the velocity of the center initiation model 1. From the above data, it can be seen that the velocity of the double-column model 1 is 5%–7% higher than that of the single-column model 1, the velocity of the double-column model 1 is 10%–17% higher than that of the center initiation model 1, and the velocity of the single-column model 1 is 5%–13% higher than that of the center initiation model 1.

4.4.2. The Velocity Comparison of the Second Group

The velocity comparison diagrams of the second group are shown in Figure 14.

In the velocity comparison diagrams of the second group, the velocity of 1#EFP is m21 > k21 > , which is (m21 = 2368 m/s) > (k21 = 2311 m/s) > ( = 2008 m/s). The velocity of the double-column model 2 is 2.4% higher than that of the single-column model 2, the velocity of the double-column model 2 is 15.2% higher than that of the center initiation model 2, and the velocity of the single-column model 2 is 13.1% higher than the velocity of the center initiation model 2. The velocity of 2#EFP is m22 > k22 > , which is (m22 = 2623 m/s) > (k22 = 2492 m/s) > ( = 2249 m/s). The velocity of the double-column model 2 is 5% higher than that of the single-column model 2, the velocity of the double-column model 2 is 14.3% higher than that of the center initiation model 2, and the velocity of the single-column model 2 is 9.8% higher than the velocity of the center initiation model 2. The velocity of 3#EFP is m23 > k23 > 23, which is (m23 = 2608 m/s) > (k23 = 2498 m/s) > (23 = 2247 m/s). The velocity of the double-column model 2 is 4.2% higher than that of the single-column model 2, the velocity of the double-column model 2 is 13.8% higher than that of the center initiation model 2, and the velocity of the single-column model 2 is 10% higher than the velocity of the center initiation model 2. The velocity of 4#EFP is m24 > k24 > 24, which is (m24 = 2392 m/s) > (k24 = 2278 m/s) > (24 = 1983 m/s). The velocity of the double-column model 2 is 2.1% higher than that of the single-column model 2, the velocity of the double-column model 2 is 10% higher than that of the center initiation model 2, and the velocity of the single-column model 2 is 17.1% higher than the velocity of the center initiation model 2. From the above data, it can be seen that the velocity of the double-column model 2 is 2%–5% higher than that of the single-column model 2, the velocity of the double-column model 2 is 10%–15% higher than that of the center initiation model 2, and the velocity of the single-column model 2 is 10%–17% higher than that of the center initiation model 2.

4.4.3. The Velocity Comparison of the Third Group

The velocity comparison diagrams of the third group are shown in Figure 15.

In the velocity comparison diagrams of the third group, the velocity of 1#EFP is m31 > k31 > 21, which is (m31 = 2418 m/s) > (k31 = 2236 m/s) > (31 = 2117 m/s). The velocity of the double-column model 3 is 7.5% higher than that of the single-column model 3, the velocity of the double-column model 3 is 12.4% higher than that of the center initiation model 3, and the velocity of the single-column model 3 is 5.3% higher than the velocity of the center initiation model 3. The velocity of 2#EFP is m32 > k32 > 32, which is (m32 = 2789 m/s) > (k32 = 2562 m/s) > (32 = 2429 m/s). The velocity of the double-column model 3 is 8.1% higher than that of the single-column model 3, the velocity of the double-column model 3 is 12.9% higher than that of the center initiation model 3, and the velocity of the single-column model 3 is 5.2% higher than the velocity of the center initiation model 3. The velocity of 3#EFP is m33 > k33 > 33, which is (m33 = 2771 m/s) > (k33 = 2578 m/s) > (33 = 2389 m/s). The velocity of the double-column model 3 is 7% higher than that of the single-column model 3, the velocity of the double-column model 3 is 13.8% higher than that of the center initiation model 3, and the velocity of the single-column model 3 is 7.3% higher than the velocity of the center initiation model 3. The velocity of 4#EFP is m34 > k34 > 24, which is (m34 = 2423 m/s) > (k34 = 2225 m/s) > (34 = 2075 m/s). The velocity of the double-column model 3 is 8.2% higher than that of the single-column model 3, the velocity of the double-column model 3 is 14.4% higher than that of the center initiation model 3, and the velocity of the single-column model 3 is 6.7% higher than the velocity of the center initiation model 3. From the above data, it can be seen that the velocity of the double-column model 3 is 7%–8.5% higher than that of the single-column model 3, the velocity of the double-column model 3 is 13%–15% higher than that of the center initiation model 3, and the velocity of the single-column model 3 is 5%–7.5% higher than that of the center initiation model 3.

4.4.4. The Velocity Comparison of the Fourth Group

The velocity comparison diagrams of the fourth group are shown in Figure 16.

In the velocity comparison diagrams of the third group, the velocity of 1#EFP is m41 > k41 > 41, which is (m41 = 2487 m/s) > (k41 = 2362 m/s) > (41 = 2170 m/s). The velocity of the double-column model 4 is 5% higher than that of the single-column model 4, the velocity of the double-column model 4 is 12.7% higher than that of the center initiation model 4, and the velocity of the single-column model 4 is 8.1% higher than the velocity of the center initiation model 4. The velocity of 2#EFP is m42 > k42 > 42, which is (m42 = 2814 m/s) > (k42 = 2606 m/s) > (42 = 2375 m/s). The velocity of the double-column model 4 is 7.4% higher than that of the single-column model 4, the velocity of the double-column model 4 is 15.6% higher than that of the center initiation model 4, and the velocity of the single-column model 4 is 8.9% higher than the velocity of the center initiation model 4. The velocity of 3#EFP is m43 > k43 > 43, which is (m43 = 2807 m/s) > (k43 = 2602 m/s) > (43 = 2376 m/s). The velocity of the double-column model 4 is 7.3% higher than that of the single-column model 4, the velocity of the double-column model 4 is 15.4% higher than that of the center initiation model 4, and the velocity of the single-column model 4 is 8.7% higher than the velocity of the center initiation model 4. The velocity of 4#EFP is m44 > k44 > 44, which is (m44 = 2506 m/s) > (k44 = 2320 m/s) > (44 = 2124 m/s). The velocity of the double-column model 4 is 7.4% higher than that of the single-column model 4, the velocity of the double-column model 4 is 15.2% higher than that of the center initiation model 4, and the velocity of the single-column model 4 is 8.4% higher than the velocity of the center initiation model 4. From the above data, it can be seen that the velocity of the double-column model 4 is 5%–7.5% higher than that of the single-column model 4, the velocity of the double-column model 4 is 12%–16% higher than that of the center initiation model 4, and the velocity of the single-column model 4 is 8%–9% higher than that of the center initiation model 4.

5. Conclusion

In this study, the circumferential MEFP warhead technology is combined with the directional warhead technology, and three kinds of initiation modes are designed and simulated. The results of numerical simulation show that when MEFP warhead adopts central single-point (multipoint synchronous) initiation, the number of shaped projectiles is the most, and the velocity of the central multipoint synchronous initiation is higher than that of the central single-point initiation; however, compared with the eccentric initiation MEFP warhead, the EFP velocity is lower. Compared with the central initiation, when the MEFP warhead adopts eccentric detonation, the number of formed projectiles is smaller, but the velocity is higher; among them, the MEFP warhead adopts double-column eccentric initiation. The maximum velocity can reach 2814 m/s. Through the comparative analysis of the forming velocity of three types of initiation modes, it can be seen that the velocity of double-column multipoint eccentric synchronous initiation is 2%–9% higher than that of the single-column single-point (multipoint eccentric synchronous) detonation, the velocity of double-column multipoint eccentric synchronous initiation is 10%–17% higher than that of the central single-point (multipoint synchronous) initiation, and the velocity of single-column single-point (multipoint eccentric synchronous) initiation is 5%–17% higher than that of the central single-point (multipoint synchronous) initiation. Different from the directional fragment warhead, the EFP warhead cannot increase the number of EFPs in the directional area, that is to say, the density gain of the kill element. However, it can improve the velocity of local EFP, that is, the velocity gain of the kill element, so as to improve the penetration ability of EFP to targets and improve the damage efficiency of the MEFP warhead.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by the Science and Technology on Electromechanical Dynamic Control Laboratory, China, no. 6142601200408.