Calculation Model of the Axial Dispersion Distribution of Cylindrical-Warhead Fragments

. The dispersion distribution characteristics of prefabricated fragments from an improvised explosive device (IED) are essential in defense technology development. To improve warhead design, it is essential to predict a fragment’s velocity and dispersion distribution accurately. This paper investigates a fragment’s ﬂying angle from a detonated warhead with a cylindrical charge under a one-end central detonation. A modiﬁed formula for calculating the fragment dispersion distribution at each axial position of the warhead was obtained based on experimental data and model analysis. The impacts of the following parameters on the ﬂying angle were considered in the formula: ratio of charge mass to shell mass, charge diameter, axial sparse eﬀect, and relative axial shell position. The formula was veriﬁed by experimentation to be suitable for wide applications using diﬀerent validation samples. Theoretical calculations with the formula show good accuracy in predicting the ﬂying angle of cylindrical-shell fragments, agreeing with experimental data. The proposed model meets the error requirements of engineering applications. This work provides a theoretical foundation for subsequent research on damage eﬀectiveness and IED damage assessment.


Introduction
A cylindrical charge warhead filled with fragments (e.g., metal balls) is a typical structure for warhead design and is called a fragmentation warhead [1].After detonation, the high-pressure fragments are transported at high velocities and are the primary source of damage [2].Being able to predict the fragment's trajectory accurately will improve warhead design.Previous researchers have proposed a series of empirical formulas to estimate the performance of cylindrical charge warheads.Among these formulas, the Gurney formula is representative in calculating the fragment's velocity as follows [3]: where V Gurney is the initial velocity of the fragment, �� � 2E √ is the Gurney energy, and β is the ratio of charge mass to shell mass.e Gurney formula does not consider the sparse wave effect on the fragment's velocity at the two ends of the warhead [4].After the Gurney formula was proposed, many researchers have improved his theory.For example, Zulkoski [5] and Randers Pehrson [6] independently proposed their methods to modify the Gurney formula.However, their model yields large errors for the terminal velocity.erefore, Huang et al. [4] proposed a more accurate formula to try to estimate the fragment's velocity of a cylindrical shell detonated at one end as follows: where L and d are the length and inner diameter of the cylindrical shell, respectively, and x is the axial distance from the detonation location (0 ≤ x ≤ L, where x � 0 represents the detonation location).is formula has been widely used and is a modification to the Gurney formula by applying the axial sparse wave effect.
To predict the warhead flying angle, Taylor [7] proposed a representative analytical formula as follows: where D is the detonation velocity of the explosive, and θ/2 is the angle between the fragment's flying direction and the normal direction of the projectile body, which is called the "Taylor angle." Figure 1 is adapted from Felix et al. [8], Carlucci and Jacobson [9], Deshpande [10], and Walters and Zukas [11], representing the variables in the Taylor formula.
Here, V N is the fragment's velocity when the fragment's flying direction is 90 °from the fragment, and V A is the fragment's velocity when the fragment's flying direction is 90 °from the shell's axial direction.Many experimental data have shown that the Taylor formula is inaccurate in predicting a fragment's flying angle, especially for the fragments near the two ends of the warhead.
To improve the accuracy of the Taylor formula, many researchers [12][13][14][15][16][17] have modified the calculation approach of the Taylor angle.However, there are still significant calculation errors found when compared with empirical data.Other researchers [18,19] proposed an equation based on a simplified simulation method to estimate the initial flying angle of the fragment by applying the influence of fragment acceleration as follows: where τ e is a characteristic acceleration time, and V ′ is the spatial derivative of the fragment's initial velocity along the axial direction of the warhead.It is not easy to estimate these two parameters.Moreover, they cannot be measured experimentally so can only be obtained by simulation.erefore, this formula has considerable limitations.Felix et al. [8] proposed formulas to calculate the initial flying angle based on shell expansion shape as follows: ( where λ is the expansion coefficient of the shell measured from the detonation end (denote λ � 1 at 20% of the axial distance away from the detonation end), and x is the axial position relative from the detonation end; M � 0.6 assumes that both the shell expansion and the fragment's flying velocity reach their maximum values at the relative axial position of 0.6.e use of this formula (6) requires known inputs of the fragment's flying angle reduction at the endplate and the position of the maximum fragment velocity.It is difficult to obtain these two input data experimentally.erefore, the relative position and reduction accuracy are mostly assumed values and hence inaccurate.us, the above-modified formulas also have significant limitations.
In order to overcome the disadvantages in the models mentioned above, in terms of the impact of the axial sparse effect on the fragment dispersion distribution, this paper proposes a model to predict the dispersion distribution of the cylindrical warhead under one-end detonation by considering the axial sparse effect and the relative position along the shell's axial direction.e model was developed based on the experimental data using different ratios of charge mass to shell mass and different charge diameters.Experimental data of a fragment's initial velocity and flying angle under different sample conditions were obtained.e impacts of the detonation wave reaction and sparse wave propagation on fragment dispersion were analyzed.Parametric function curves derived from the proposed model were established.Finally, the theoretical model was validated by using static detonation tests.

Samples and Experimental Layout
In order to study the axial dispersion characteristics of the fragments driven by an explosion, a static explosion test was carried out on an explosive device, as shown in Figure 2. e spatial fragment distribution was measured by the penetration positions on a witness plate.e sample explosive was placed on flat plank support 1.5 m above the ground.
e witness plate was made of carbon structural steel and had dimensions of 1570 mm in length, 3000 mm in height, and 3 mm in thickness.
e covering angle between the warhead and the witness plate was 18 °.e witness plate was placed 5 m away from the warhead.e detonation method was a one-end detonation from the center.e test sample was placed vertically with the detonation end facing upward.A high-speed camera system, which included a FASTCAM ultima APX high-speed camera produced by the Photron Company, was used to monitor and record the fragment's perforation process through the target witness plate.e camera shooting rate was set at 13,000 frames per second.Shock and Vibration e horizontal distance between the explosion-driven device of the warhead and the high-speed camera was 35 m.After the test, the witness plate was recycled and photographed to obtain perforation conditions, and the perforation positions were measured.
According to the structure shown in Figure 3, six sets of cylindrical specimens with different internal diameters and aspect ratios (i.e., the ratio of length to diameter) were prepared.e cylindrical samples were equipped with the JH-2 explosive, 45 steel sleeves with a 3 mm thickness, side end covers with a 5 mm thickness, and steel balls with a 5 mm diameter.e charge was pressed into the samples under a 300 MPa pressure using 300-ton pressing equipment to obtain a charge density greater than 1.69 g/cm 3 .e internal and external sleeves were fixed using a threaded connection.e 5 mm diameter steel ball fragments were put into the sleeve in sequence.An epoxy resin was injected into the sleeve to fix the steel balls.All samples had the same density of the 8,701 explosives and the same sleeve thickness.
e length-to-diameter ratios of the six samples were 2 : 1, 1.8 : 1, and 1.6 : 1, and the charge diameters were 56 mm, 70 mm, and 89 mm. e charge mass to shell mass ratios (C/M) of the six samples were 0.778, 0.579, and 0.743.Table 1 lists detailed specifications of the samples.

Analysis of High-Speed Photography Results.
e highspeed camera recorded the fragment's perforation process through the target witness plate.According to the time marks on the photo, the fragment's velocity at a given distance was calculated to obtain all the fragment data within the coverage angle of the witness plate.After the test, the image was converted into the grayscale bitmap, and brightness and contrast were adjusted using a data acquisition system connected to the high-speed camera to reveal the process of fragment perforation, as shown in Figure 4.
Figure 4 shows that obvious firelight appeared on both sides of the witness plate after the charge exploded.As time went by, the firelight's brightness first increased and then gradually decreased.e frame before the firelight occurrence was set to 0 in the time scale.It was found that after the charge explosion, the firelight produced by fragment perforation on the witness plate appeared at 30.8 ms and lasted till 38.4 ms. e fragments spread and increased from the center of the witness plate to both sides.
e increased number of fragments on the lower part of the witness plate was greater than that on the upper part.e fragment's flying velocity first increased and then decreased from the detonation end to the explosion center and then to the other end.In order to calculate the fragment's trajectory, the trajectory of the fragment was assumed horizontal, and the impacts of air lift and gravity were ignored.us, only the effect of air resistance was considered.According to the velocity attenuation formula of a spherical fragment, the fragment's motion equation was expressed as follows: where v is the instantaneous velocity of the fragment, c f is the air resistance coefficient of fragment flying (c f � 0.97 for a spherical prefabricated fragment), and A is the average windward area of the fragment.Equation ( 7) was integrated to obtain the following equation that links the fragment flying velocity v f and the flying time t f : According to the above equation, the fragment's initial velocity was calculated based on the measured velocity.
en, these data were used as a basis for establishing a model to predict the dispersion distribution of the cylindrical warhead under one-end detonation.e measured values of Shock and Vibration initial warhead velocity by using the high-speed camera are shown in Table 2.

Witness Plate Analysis.
After the test, fragment perforation on the witness plate was analyzed.e sudden change in the gray value of each edge pixel of the perforated hole was used to segment the image, identify the target hole area, and extract the hole shape.e Image-Pro Plus (IPP) and Image J analysis software packages were used to measure and locate the contour of the hole on the witness plate.e number of holes was counted, as shown in Figure 5.
Because the perforated hole was the same as the fragment in both size and shape, the IPP technique was effectively used to obtain the perforated holes produced by the prefabricated fragments.e perforation conditions of the six warhead samples were collected and analyzed to obtain their fragment dispersion distributions.e samples with the same ratio of length to diameter had consistent dispersion distributions.e test data were reliable with few experimental errors, and therefore, they were used to build a model to predict the fragment dispersion distribution, as shown in Figure 6.
e derived fragment dispersion distribution curves showed that the warheads with different aspect ratios were affected by the sparse waves at the two ends of the shell, exhibiting similar nonlinear distributions.e fragments in the middle part flew steadily, and the fragments at the two ends flew with more deviated flying angles.e decrease or increase in the flying angles at the two ends depended on the end effect of the sparse waves, which mainly affected the dispersion distributions on the detonating and nondetonating ends; the detonation end was more affected by the sparse wave than the nondetonation end.As mentioned above, when the explosion was detonated, the axial sparse wave was rapidly transmitted from the detonation end to the explosion products, inhibiting the load given by the detonation wave near the detonation end.e axial sparse wave propagating from the nondetonation end to the explosion products was generated only when the detonation wave   propagated to the nondetonation end.erefore, the nondetonation end was only affected by the axial sparse wave in a very small time interval, and the corresponding fragment flying angle was less affected.e flying angle presented a distribution trend that it first rapidly increased in an arc pattern, then became flat, and finally slowly increased in an exponential form.e impacts of the axial sparse waves at the two ends on the fragment's flying angle gradually separated.e explosion-driven dispersion in the middle part of the shell gradually completed in an intermediate transition period of the axial sparse wave.erefore, the effects of the axial sparse waves from the two ends on fragments flying off the middle part were gradually weakened.

Model Establishment
is function is an odd function that can be decomposed into a convex function and a concave function that is monotonically increasing.e hyperbolic cosine curve could not construct the empirical model because there were negative values in the fragment dispersion curves obtained experimentally.
Previous formulas could not reflect the actual trajectories of warhead fragments under the impact of the sparse waves at the two ends of the shell.e detonation wave had different shapes at different axial positions, causing different impulses and pressures sustained by the inner wall of the shell and different loads along the axial direction at different positions along the projectile body.erefore, there was an axial expansion in the projectile.Related terms of the sparse waves at the two ends were introduced to the model to describe the propagation of the axial sparse waves towards the two free ends to correct the predicted fragment flying angle.
As shown in Figures 7 and 8, the dispersion distribution was affected by the coupling effect of the axial sparse waves from the two ends.Huang et al. [4] studied the fragment velocity distribution of a warhead with a charge of a large aspect ratio.ey found that the empirical formula of a single variable could not describe the interaction between them.ey proposed a widely used two-variable three-dimensional empirical formula to determine the coupling effect of the sparse waves on the fragment dispersion distribution at different coaxial positions.is paper proposes a modification by using a fragment dispersion distribution model tanθ(x, y) to cover the effect of different aspect ratios as follows:  where the independent variable x � Lx/D is the proportional distance of the fragment to the detonation end, y � Ly/D is the proportional distance of the fragment to the nondetonation end, θ is an ideal fragment dispersion distribution model without the impact of the axial sparse waves.f(x) and f(y) represent the impacts of the axial sparse waves from the detonation and nondetonation ends on the fragment dispersion distribution, respectively.e Taylor formula can explain the impact of the axial sparse waves on the dispersion distribution at the two ends.It can also explain why the impact of the axial sparse waves on the maximum flying angle changes with an increase in L/D as a hyperbolic sine function.Based on the empirical model of fragment dispersion using the hyperbolic sine function, the modified term for the detonation end can be constructed as follows: where a 1 and b 1 are the correction coefficients reflecting the impact of the sparse waves from the detonation end.
Similarly, the correction formula for the nondetonation end can be expressed as where a 2 and b 2 are the correction coefficients reflecting the impact of the sparse waves from the nondetonation end.erefore, the calculation model of fragment dispersion distribution can be further developed as follows by using the empirical model with the hyperbolic sine function and the Taylor formula: Test data were substituted into the above model, and the data were fitted with a nonlinear least square regression equation.e modified distribution function curve of the flying angle is shown in Figure 9. Test data determined the correction coefficients.
e fitted correction result is as follows:

Model Error Analysis.
e R 2 value of parameter fitting was 0.9885, indicating that equation ( 14) fit the test data very well and describe the fragment dispersion distributions under the impact of the sparse waves for different diameters and charge aspect ratios, especially the dispersion distributions at the detonation and nondetonation ends.6 Shock and Vibration accurately predict the fragment dispersion distribution at each axial position of the cylindrical warhead under one-end detonation.erefore, the model has significant improvement in accuracy and application range.
Figure 10 shows the comparison between the calculation result and test data.e degrees of freedom (DF), sum of squares, and mean square values of the residuals are 72, 12.86, and 0.19938, respectively.e residuals between the relative axial position and the fitted Y value are mainly distributed between −1 and 1.5, as shown in Figures 10(a) and 10(b).With the line passing the points (−0.93, 0.66) and (0.93, 99.33) as a reference line, the conventional residuals of the percentile are no more than 0.3, as shown in Figure 10(c).
e model-fitting results are accurate and in good agreement with the test data.e relative errors of 97% of the data at all axial positions are within 3%.e average error of the six warhead samples is within ±2%.Some large relative errors are due to the impact of the sparse waves from the detonation and nondetonation ends, causing uneven forces acting on the fragments, and hence, dispersion distributions are unstable.e proposed model can largely avoid these errors.In general, this empirical model can be used to calculate the fragment dispersion distribution under the effect of the sparse waves with one-end detonation for a cylindrical charge warhead.

Test Verification to Validate the Model
5.1.Samples.In order to verify the universal applicability of the proposed model, several different cylindrical charge warheads were designed.Two samples with the JH-2 explosive were assembled with an aluminum shell, a charge diameter of 30 mm, and a length-to-diameter ratio of 2 : 1.
e C/M ratio of the specimen was 0.295.e sleeve thickness was 2 mm.e prefabricated fragment was a tungsten ball of 3 mm diameter.Other detailed structural parameters are listed in Table 3. e JH-2 explosive used in the validation warhead samples was taken from the same batch of the explosive used in earlier tests.e charge was pressed at a 300 MPa pressure with 300 -ton pressing equipment to obtain a charge density greater than 1.69 g/ cm 3 .e charge length of each specimen was the same as the sleeve length.e test site layout is also shown in Figure 2.

Model Verification Results
. According to the high-speed camera photography results and the fragment velocity attenuation formula, the fragment's initial velocity was calculated based on the measured velocity at the measurement location.
e initial fragment velocities of the two test samples were 1205 m/s and 1213 m/s.e fragment dispersion and target penetration conditions were analyzed by using the witness plate after the test.
e images were processed with the previously mentioned approach using the IPP image analysis software to measure the contours of the perforated holes and determine their positions on the witness plate, as shown in Figure 11.e fragment dispersion distribution was obtained by effectively extracting the perforation information from the witness plate, and the proposed empirical formula was used to calculate the fitted dispersion curves.e comparison between the validation test data and the model prediction is shown in Figure 12.

Model Verification Error Analysis.
e R 2 value of parameter fitting was 0.934, indicating a good fit between the model and the validation test data.Figure 13 shows the relative difference between the model result and test data.
e DF, sum of squares, and mean square values of the residuals are 56, 56.48, and 1.22, respectively.e residuals between the relative axial position and the fitted Y value are mainly distributed between −2 and 0, as shown in Figures 13(a) and 13(b).With the line passing the points (−2.49, 0.62) and (0.82, 99.37) as a reference line, the conventional residuals of the percentile are no more than 0.4, as shown in Figure 13(c).e model-fitting results are accurate and in good agreement with the validation test data.
e relative errors of 97% of the data at all axial positions are within 6%. e average error of the two warhead samples is within ±7%.Shock and Vibration e calculated model data in this paper are compared with flying angle prediction results of Taylor and Flix et al. [7,8] formulas, and the comparison results are shown in Figure 14. e prediction of the fragment's flying angle given by the Taylor formula was close to a straight line parallel to the X/L axis.When such a predicted result was compared with the experimental data, the correlation coefficient could not be estimated because the standard deviation of the predicted result was zero.By making a minor adjustment to the straight line, e.g., reducing the value of a selected point on the straight line by 10%, the correlation coefficient became approximately 0.3, indicating that the Taylor formula has little correlation with experimental data.As to the flying angle prediction model developed by Flix, although its predicted trend of the flying angle distribution was correct,   8 Shock and Vibration its accuracy was poor compared with the actual data, especially at the detonation and nondetonation ends where the correlation coefficient was only approximately 0.64.erefore, on the one hand, compared with the previous models, the proposed model in this study can accurately predict the flying angle distribution.On the other hand, the model still needs to be further improved for better accuracy.For example, for the fragment dispersion distribution at the remaining 3% axial positions near the tail of the nondetonation end, distribution errors need to be reduced and subject to future research.

Results and Discussion
In this study, the charge diameter and the ratio of charge mass to shell mass were varied to build different test samples to obtain the test data of a warhead fragment's initial velocity and flying angle under different conditions.e impact of the sparse waves on fragment dispersion was analyzed.A calculation model was proposed to predict fragment dispersion distribution under the effects of the axial sparse waves and the relative axial shell position under one-end detonation for cylindrical warheads.e model was completed by experimental data fitting.e research results show the following findings: (1) Different warheads of the cylindrical charge with different aspect ratios were affected by the sparse waves from the two ends of the shell, showing similar nonlinear distributions.e middle fragments scattered steadily, and the fragments at the two ends scattered with larger deviations in the flying angle.
e decrease or increase trend in the flying angle on both sides was related to the end effect of the sparse waves, especially in the dispersion distributions at the detonation and nondetonation ends.e impact of the sparse waves on the detonation end is greater than that on the nondetonation end.e flying angle presented a nonlinear distribution with two ends rising rapidly and the middle part increasing steadily.
(2) e proposed model was completed using curve fitting with experimental data based on the hyperbolic sine (sinhx) function curve.e model can accurately predict the fragment distribution under the impact of the sparse waves, especially the detonation and nondetonation ends.With the parameters of x/d and (L − x)/d, the model can predict the fragment dispersion distribution at each axial position of the warhead.e residual differences between the relative axial positions and the fitted Y values were between −1 and 1.5, and the conventional residual in the percentile was less than 0.3.e relative error in 97% of the data at all axial positions was within 3%. is means that the proposed model is more accurate than the previous models.given by the Taylor and the Flix formulas, the

Figure 1 :
Figure 1: Variables used in Taylor equation.

Figure 3 :
Figure 3: Schematic diagram of the sample structure.

Figure 4 :
Figure 4: Target penetration process of prefabricated fragments driven by explosion.

4. 1 .
Empirical Model Establishment of the Flying Angle.It is necessary to establish an empirical model to describe the fragment dispersion to study further the impact of the axial sparse waves under different charge diameters and aspect ratios on the dispersion distribution.Based on the above test results of dispersion distribution, the fragment dispersion distribution curve was similar to a hyperbolic sine curve (sinhx) in their function variation characteristics.e hyperbolic sine curve was translated and processed by coefficient transformation so that the curve had positive and negative values at x � 0 -1 (relative position in the axial direction).erefore, such a new hyperbolic sine function was obtained as follows:

Figure 5 :
Figure 5: Picture of pit treatment on witness plate after test.

Figure 6 :
Figure 6: Test value of combat fragmentation dispersion distribution.

Figure 7 :Figure 8 :
Figure 7: Sparse wave propagation diagram at two free ends of the sleeve.

Figure 12 :
Figure 12: Test value of combat fragmentation dispersion distribution.

( 3 )
A separate set of validation tests was designed to confirm the model's accuracy.e validation results showed that the relative errors in 97% of the data along the axial direction of the warhead were less than 6%, indicating good accuracy for different cylindrical warheads.e model is essential for the risk assessment of improvised explosive devices.(4) Compared with the flying angle prediction results

Figure 13 :
Figure 13: Error analysis for verifying the model.(a) Axial relative position error.(b) Dispersion distribution error.(c) Percentile error.

Table 1 :
Parameter specifications of cylindrical-warhead samples with charge and fragments.

Table 2 :
Measured values of the initial velocity of the warhead fragment.

Table 3 :
Specimen parameters used in the validation test for the proposed model.