Penetration of ballistic gelatin by explosion-driven inert metal particles

In urban combat, sub-millimeter inert metal powder is used to replace fragments in explosive bombs, which can effectively reduce collateral damage. In order to investigate the damage effect of high-speed metal particles on the human body, a particle ring filled with a charge was designed to create an explosion-driven particle cloud for the penetration ballistic gelatin. The particle ring was made of polytetrafluoroethylene (PTFE) and sub-millimeter tungsten powder. The dispersion properties of the particle cloud driven by the explosion were studied with high-speed photography and ballistic gelatin. Furthermore, the numerical simulation models of particle-penetrating gelatin based on the experimental results were established with the finite element method. The influences of the particle size and velocity on the penetration depth and cavity diameter was obtained. The relationship between the critical interference distance of the cavity and the distance between particles was determined. This paper provides references and theoretical support for the design of low collateral damage ammunition based on inert metal powder.


INTRODUCTION
In recent decades, local armed conflicts and terrorism have become more frequent.Statistics from the Global Terrorism Database (GTD) show there were over 100,000 terrorist attacks between 2005 and 2019.Most terrorist attacks occur in densely populated urban environments and civilians are often close to the terrorists during an engagement.Therefore, avoiding civilian casualties when eliminating terrorists must be considered in anti-terrorism operations (Walsh. 2014).Fragmentation bombs that can produce high-speed metal fragments are very effective in eliminating enemy personnel who are within the killing radius.However, those bombs are not suitable for anti-terrorism operations because the fragments still have a large range of lethality when flying outside of the killing radius, which is the so-called high collateral lethality.Therefore, extensive research has been conducted to develop low collateral lethality ammunitions.
Usually, the killing range of fragment grenades can be effectively reduced and the harm to innocent people can be lowered by reducing the amount of charge, using non-metallic shells and other measures (Yang et al. 2018, Wang et al. 2020, Yu et al. 2023, Liang et al. 2017).However, these measures also lead to a smaller killing radius and a lower power of ammunition.Many studies have shown that the use of the sub-millimeter size of inert metal powder instead of fragments as the damage element in fragment grenades is an effective way to achieve low collateral damage.One method is to evenly mix the metal powder into explosives, forming heterogeneous explosives.When this kind of explosive explodes, the metal powder can produce a uniform metal powder flow with a high initial speed during detonation (Zhang et al. 2012, Frost et al. 2007).Many researchers have conducted in-depth studies on its gas-solid two-phase flow mechanism (Sugiyama et al. 2020, Xue et al. 2014, Zhang et al. 2001), sensitive conditions for particles through shock waves (Balakrishnan and Menon 2008), and detonation propagation and quenching rules (Ju and Law 2002).However, some studies have indicated that heterogeneous explosives containing metal particles are prone to producing hot spots when subjected to impacts, leading to the issue of munition safety.An alternative method, that is used in this research, is wrapping the metal particles around the explosives to form a dense powder ring that acts as a fragile shell.The warheads based on this structure are much safer because they can substantially reduce the possibility of producing local hot spots.
Many researchers have analyzed the killing effect of the metal particle cloud of this structure through experiments.Yao et al. (2010) compared the pressure curves after the explosions of a bare charge, non-metallic shell charge, and metal powder shell charge, and they found that the metal powder shell charge had a larger peak pressure and a significantly longer action time after an explosion, proving that the metal powder shell charge had a smaller killing area and a stronger killing effect within a smaller killing area.Chen et al. (2011) used the multi-layer paper target method to analyze the kinetic energy characteristics such as the number density, velocity and distribution of metal particle groups driven by the explosive.In previous experiments (Chen et al. 2017), it was found that metal particles were prone to sintering at the high temperature and pressure of detonation, resulting in uneven particle spatial distribution.Bai et al. (2010) evaluated the metal powder dispersion situation of different loading processes within an effective distance base on aspects such as the particle spatial number density, the particle penetration ability for the target plate, the process and safety.Furthermore, the addition of a dispersant is beneficial for reducing the possibility of metal powder sintering and agglomeration.Thus, the possibility of PTFE powder as an anti-sintering agent was explored based on this in this research.
The acceleration and dispersing process of metallic particles by the detonation of an explosive has been well studied.However, the penetration effects of high-speed metallic particles on human tissues have rarely been investigated.Gelatin is often used for a design that simulates human tissues (Liu et al. 2017, Moxnes et al. 2016, Liu et al. (2021b), Han et al. 2020) because its mechanical properties are close to those of human muscle.Thus, the ballistic impact into gelatin has become a conventional method for evaluating the ballistic wound of a human (Grimal et al. 2004).Gilson et al. (2020) conducted experiments analysis on the non-penetration ballistic response of a composite plate shielding gelatin block and evaluated the front surface deformation, pressure wave amplitude value, and duration of the gelatin block.Liu et al. (2021a) studied the translation and rotation of rifle bullets in ballistic gelatin.Wen et al. (2012) numerically studied the penetration of steel balls into gelatin.It was shown that the Lagrange solver had the highest calculation accuracy and efficiency.Ye et al. (2022) conducted dimensionless analysis and a finite element simulation on the penetration of medium-low speed millimeter-level steel balls into ballistic gelatin, which reproduced the experimental observation results well.Shuangshuang et al.(2020) developed an original numerical model based on SPH to study the dynamical phenomena during microscale impact of spheres into ballistic gelatin.Mitchell et al. (2003) described the results of an investigation into the impact of model micro-particles to ex vivo buccal mucosa (the cheek) of pigs and beagle dogs.Veysset et al. (2018) investigated the high-velocity impact response of gelatin and synthetic hydrogel samples using a laser-based microballistic platform for launching and imaging supersonic micro-particles.The above studies focused on the movement of millimeter-level projectiles in ballistic gelatin and the response of ballistic gelatin, but there have been relatively few studies on the damage effect of high number density micro-particles with smaller sizes on ballistic gelatin.
Latin American Journal of Solids and Structures, 2024, 21(3), e535 3/18 In this investigation, the penetration effects of the high number density metallic particles into ballistic gelatin were studied experimentally and numerically.Explosion experiments were first performed on a gelatin target using a tungsten particle ring filled with different-mass high-explosive charges.The numerical models of the tungsten particles penetrating the ballistic gelatin were established and validated by the experimental results of the penetration depths.Parametric studies were then carried out using this numerical model to determine the influence of the particle size and the spacing between two particles on the penetration depth and maximum cavity.

Preparation of particle ring and charge
The materials for preparing dense particle rings include tungsten powder (average particle size: 300 μm, density: 19.35 g/cm 3 ) and PTFE powder (average particle size: 40 μm, density: 2.15 g/cm 3 ).Table 1 shows the physical properties of the two powders.The particle ring (ω PTFE =13.3 %, ω W =86.7 %) preparation process is shown in Figure 1.PTFE and tungsten powder are uniformly mixed for 15 minutes with a small drum mixer.Pour the mixed powder into a tool, and press with a pressure of 320 MPa for 3 minutes.The density of the particle ring was 7.99 g/cm 3 .Figure 1 Particle ring preparation process: (a) mold for compressing the particle ring and (b) process of compressing the particle ring.
The particle ring was sheathed on the outer side of the JH-2 cylindrical explosive.The charge density was about 1.68 g/cm 3 .Nylon was used for the top and bottom covers.Between the explosive and the particle ring, the outside of the particle ring is a CFRP shell for enhancing the overall strength.The specific structure is shown in Figure 2 (a).The charges were named D40, D80, and D93, as shown in Figure 2 (b).The number represents its outside diameter.The parameters of the particle ring and charge are shown in Table 2.

Experimental set-up
The particle ring filled with charges was hung up vertically at the selected position, with the center 1.5 m from the ground.The three 10%wt gelatin blocks (side length was 300mm) were placed on three rigid support, as shown in Figure 3 (a) and Figure 3 (b).The supports were placed at a distance R i from the blast center (where i =1, 2, 3).Each set of experimental detonation distances is determined by considering R 2 , which is 40 times the maximum diameter of this structure.R 1 is less than 40 times the maximum diameter, while R 3 is greater than 40 times the maximum diameter, as indicated in Table 3.The field of the explosion experiment was arranged as shown in Figure 3 (c).The processes of explosive detonation and particle dispersion were recorded by a high-speed video recorder (6000 fps).
Table 3 The distances between the explosive centers of the three types of charges and their corresponding gelatin.Latin American Journal of Solids and Structures, 2024, 21(3), e535  5/18

Particle cloud boundary initial velocity and decay
Figure 4 shows the typical explosion moments from the perspective of high-speed photography.The explosive volume expanded rapidly to form a high-temperature and high-pressure detonation product.The particle ring was compressed, and then expanded outward through momentum exchange with the detonation product.The particle ring broke into fragments, and then gradually detached from the detonation product.Due to the interface disturbance of the detonation products and the influence of the velocity gradient of the particle fragments, the particle fragments continued to break up and disperse to form particle flow.The arrival position of the powder boundary at various time intervals is delineated on the high-speed image, with reference to the baseline direction in Figure 4.The average radius from the explosion center to the cloud boundary was measured at different time points, as depicted in Figure 5.The modified Gurney equation accounting for the initial velocity of particle cloud which is primarily a function of the mass ratio, the ratio between the particles and the charge.
where λ=m/m p , 2E is a test constant that represents the performance of explosives, k is the energy loss coefficient considering the outflow of detonation gas.It was measured that the maximum throwing velocity of particles for the three cases were basically the same, at 1178.4 m/s, 1103.4 m/s, and 1138.3 m/s.The average speed of the three working conditions was 1140.0 m/s.It is finally determined that k is 0.89 which is more reasonable.A simplified deceleration model was employed, assuming one-dimensional particle motion and disregarding inter-particle interactions.When the particles were imparted with a specific initial velocity to travel through still air, the primary influence on the particles was aerodynamic resistance.In accordance with aerodynamics, the equation governing particle motion could be derived as follows: where m p0 is the mass of the particle, v is the velocity of the particle in the air, t is the time, C x is the resistance coefficient, ρ air is the air density, and A is the current windward area of the particle.The resistance coefficient C x of a compressible viscous fluid flowing around a particle can be expressed as (Roache. 1976): where d p is the particle size, and μ is the aerodynamic viscosity coefficient.At room temperature of 25 °C, μ is 0.01834 mPa•s.Then the relationship between the particle's velocity and time could be expressed as where is the characteristic time of the velocity relaxation process and φ is the particle shape correction coefficient.As it approaches 1.00, the particle shape is closer to spherical.The initial velocity of the particle is v 0 , and when t =0, v(0) =v 0 , and the particle cloud radius R c could be expressed as It is shown that when the initial velocity is constant, the change in the particle cloud radius with time depends on the particle size and density.The particles used in the experiments were identical.A series of results captured by high-speed photography were fitted to the particle shape correction coefficient φ by Eq. ( 5).The correction coefficient of the particle shape is φ=0.95,where R 2 =0.995.It proves that the metal particles used in the experiment are relatively close to spherical.The curves of the particle cloud radius and particle velocity with time are shown in Figure 5.The ballistic gelatin after the experiment was recovered, as shown in Figure 6.It could be clearly observed from the experiment results that a large number of small blind holes were present.These holes were evenly and densely distributed This phenomenon indirectly proved that the addition of PTFE powder to metal powder could effectively reduce the sintering and agglomeration of particles.
The embedding of the particles in the gelatin could be seen more clearly by shining bright light on the gelatin.Due to the influence of uncertain factors such as particle speed, particle size, shape and angle of incidence, the penetration depth of particles was different.The gelatin image was binarized, as shown in Figure 7.The black dots in the picture are particles stuck in the gelatin.The particle cloud boundary velocity and the corresponding gelatin penetration depth at different detonation distances are shown in Table 4.The particle cloud number density was obtained by recovering particles in gelatin.First, the gelatins after the experiments were placed in a water bath and heated to melt, and the liquid was allowed to stand for precipitation.The particles were precipitated as much as possible.After the removal of most of the clear solution in the upper layer, the remaining liquid was filtered with filter paper.The collected particles were washed repeatedly and dried to remove excess water.Finally, the mass of the tungsten particles in gelatin was obtained.The number density of particles is represented by the number of particles per unit area, and the results are shown in Figure 8.The factors affecting the number density of particles were mainly the mass of the particle ring and the detonation distance.The geometric relationship of the particle cloud impacting the gelatin is shown in Figure 9.It was assumed that the distribution of the particle cloud in the circumferential direction was uniform.As shown in Figure 9 (a), the value of α is very small and can be approximated as a/R i .a is the gelatin side length, m W is the metal particle mass in the particle ring, and R i is the detonation distance (i= 1, 2, 3).Then the particle mass in the projection m α of angle α is It is also assumed that the dispersion of the particle cloud in the axial direction is uniform and symmetrical along the horizontal plane of the center of the charge.As shown in Figure 9 (b), the area between two rays AA' and BB' is the dispersion range of the particle cloud.β is the angle between AA' and the horizontal direction, which represents the concentration of particle dispersion in the axial direction.The distance between A' and B' is L. β and L can be obtained from the geometric relationship that Latin American Journal of Solids and Structures, 2024, 21(3), e535 9/18 Then the ratio of the particle mass m G invading the gelatin to the particle mass m α in the α angle is the ratio of the length of the A'B' segment to the size of the gelatin: Substituting Eqs. ( 6), (7), and (8) into Eq.( 9) The above equation shows the relationship between the particle mass m G invading the gelatin and the total tungsten particle mass m W and the distance R i .When the particle recovery mass at different distances is known by the experiment, β can be obtained by Eq (9), as shown in Figure 10.The dispersion angle β is 13-14°.It could be concluded that the particle dispersion driven by the explosion of this structure is quite concentrated  10) (Fitting the experimental results to get β, D40 is 14.36°, D80 is 13.33°, D93 is 13.01°.)

Finite element model
The numerical simulation model of the particle penetration process was established using Ls-DYNA commercial finite element software.For submillimeter particles and gelatin blocks with large differences in size, it was unnecessary and a waste of computing resources to establish gelatin blocks of the same size as the experiment.It was only necessary to ensure that the gelatin established by the numerical simulation could reflect the temporary cavity expansion and elastic-plastic deformation region.In order to save computing resources, using a cylinder instead of a square model could effectively reduce the number of elements and increase the calculation speed.The geometric dimensions of the gelatin model were a cylinder of φ40mm, and the axial length was appropriately adjusted to ensure that the gelatin was not perforated.The Lagrange solver was used, and a three-dimensional 1/2 symmetric model was established.
One such discretization is exhibited in Figure 11 (a) and Figure 11 (b).: (I) The direct impact region within a diameter of 4 mm contained fine elements.(II) The elastic-plastic deformation region within a diameter of 12 mm contained medium elements.(III) The remaining region used coarse elements.Transition meshes were used between the three regional grids.This method could greatly reduce the number of meshes, ensure the regularity of the mesh shape, and improve computational accuracy.The tungsten particles were divided into six elements along the radial direction.
The CONTACT_ERODING_SURFACE_TO_SURFACE erosion contact definition was used to simulate the interaction between particles and gelatin.The viscous hourglass control algorithm with hourglass coefficient = 0.1 was employed.Symmetry constraints were set on the symmetry plane, and all the binding surfaces of gelatin and particles were considered non-traction except for mutual contact.

Material model
The ballistic gelatin parameters of wt10% were measured by Wen et al. (2013).The density was 1.03 g/cm 3 .The yield strength can be expressed as: where σ 0 is the initial yield strength, ε p is the effective plastic strain, E is Young's modulus, and E t is the tangent modulus.The abovementioned constitutive relationship was supplemented as a polynomial of pressure P versus specific volume or density: where μ is a dimensionless parameter that defines the ratio of the current density to the initial density, and C 0 , C 1 , C 2 , and C 3 are constants.
Gelatin is described by the MAT_ELASTIC_PLASTIC_HYDRO strength model and the EOS_LINEAR_POLYNOMIAL equation of state.The specific material parameters are shown in Table 5.
1.03 850 10 220 0 2.38 7.14 11.9 It was noted that the strength and stiffness of the tungsten material were much greater than that of gelatin.The particles were assumed to be rigid.The particle was described by MAT_RIGID.The material parameters are shown in Table 6.

Grid independence analysis
The elements within region (I) were directly affected by the particle impact, and the size of these elements played a crucial role in determining the penetration results.A grid-independent analysis of region (I) was conducted separately in both the radial and axial directions.The initial input conditions included a particle diameter of 0.3mm and an initial velocity of 700m/s.
Different numbers of meshes were employed in the radial direction: 20, 28, 36, 44, and 52, as illustrated in Figure 12 (a).Correspondingly, the radial element sizes were 0.100mm, 0.071 mm, 0.056 mm, 0.045 mm, and 0.038 mm.For the axial element size of 0.08mm, the time histories of numerically calculated penetration depths for these five finite element meshes are compared in Figure 13 (a).It is noteworthy that the difference between the calculated penetration depths decreases as the minimum element size in the gelatin finite element meshes diminishes.At t = 250 μs, the penetration depths of meshes with 36, 44, and 52 elements differ by less than 4%.Considering the trade-off between calculation speed and efficiency, the discretization with 36 elements (0.056mm) yields the optimal results.
In the axial direction, various grid sizes were used: 0.32 mm, 0.24 mm, 0.16 mm, 0.08 mm, and 0.04 mm, as depicted in Figure 12 (b).With a radial element size of 0.056mm, the time histories of numerically calculated penetration depths for these five finite element meshes are compared in Figure 13 (b).Remarkably, at t = 250 μs, the penetration depths of the 0.08 mm and 0.04 mm grids differ by less than 3%.In summary, the optimal grid discretization conditions entail a radial direction element size of 0.056mm and an axial direction element size of 0.08mm.

Figure 13
Grid independence result: (a) In the radial direction, the discretization with 36 elements (0.056mm) yields the optimal results.(b) In the axial direction, the discretization with 0.08mm yields the optimal results.

Particle movement and cavity expansion
Simulation experiments were conducted to assess the accuracy of the numerical simulations in modeling single particle penetration.These simulations considered a range of particle velocities and sizes.Specifically, particle penetration scenarios were simulated for particles with diameters spanning from 0.2mm to 0.4mm, corresponding to characteristic velocities within the range of 103 m/s to 907 m/s (as indicated in Table 4).Figure 14 illustrates typical images of particle penetration, with particle sizes of 0.3mm and velocities of 500m/s, showcasing the formation of a temporary cavity characterized by the maximum diameter (d max ) and the maximum penetration depth (x max ).The viscous fluid resistance model (Swain et al. 2014) of projectiles penetrating gelatin proposed is used to describe the maximum penetration depth of particles.It is suitable for the problem of a projectile penetrating gelatin in 100-1500 m/s, with consideration of the inertia resistance and viscous resistance of gelatin.The force of submillimeter particles in gelatin can be described as where m p0 is the particle mass, v is the particle velocity, ρ is the gelatin density, S is the cross-sectional area, C 1 is the inertial resistance coefficient, and C v is the viscous resistance coefficient.μ is the velocity coefficient, b is the boundary layer thickness, and μ and b are unknown parameters related to the gelatin material.A=-ρC 1 S/m p0 , B=-μC v S/bm p0 .The relationship between the maximum penetration depth x max and the initial velocity v 0 of particles penetrating gelatin is obtained as The maximum temporary cavity can be considered directly proportional to the particle diameter (d p ) and initial velocity (v 0 ): where K is a cavity expansion coefficient.
Figure 15 (a) uses the viscous fluid resistance model to fit the numerical simulation results of the maximum penetration depth of particles with three different particle sizes from 0.2mm to 0.4mm into gelatin.The fitting values of A are -73.050, -48.479, and -36.525 respectively. The fitting values of B are -28863.380, -9322.929, and -4232.273 respectively.Latin American Journal of Solids and Structures, 2024, 21(3), e535 14/18 The particle mass m p0 and cross-sectional area S are known.C 1 is the inertial resistance coefficient, and the C 1 values of the three particle sizes are 0.183, 0.182, and 0.183 respectively.The inertial drag coefficient has nothing to do with particle size.The value of μC v /b is related to gelatin viscosity and boundary layer thickness, and μC v /b decreases with the increase of particle size.Since it was not possible to directly observe the size of the temporary cavity in the gelatin during the simulations, Figure 15 (b) provides an analysis of the numerical simulation results for the temporary cavity diameter.According to the simulation results, the curve between the initial velocity v 0 and the maximum temporary cavity diameter d max is obtained by least squares fitting, as shown in Figure 15(b).Among them, the fitted K values are 26.35,27.8 and 31.9 respectively.In Figure 16, a comparison is made between experimental results and numerical simulation results for the maximum penetration depth (x max ) under the same velocity and particle diameter conditions.The experimental results represent the depths at which particles with diameters of 0.2 and 0.4 mm penetrate into the gelatin at speeds ranging from 103 m/s to 907 m/s.Upon comparing the numerical simulation with the experimental results, significant differences are observed at low velocities, but the experimental and numerical results converge at higher speeds.It is also noticeable that there are substantial differences between the experimental and numerical results for larger particle diameters, while they are relatively close for smaller diameters.When comparing the depth range of particle distribution between experimental and numerical simulation results, the average error is 30%.Possible sources of errors mainly include: human errors caused by inconspicuous image differences between small-sized carbon fiber fragments and particles, random errors in the binary identification process, and accidental errors in particle penetration during the intrusion process.The former resulted in errors in the measurement of the penetration depth of particles, while the latter led to an increase in the penetration depth of particles at higher velocities.

DOUBLE-PARTICLE PENETRATION AND CAVITY INTERFERENCE
The presence of temporary cavities formed by the particle cloud during gelatin penetration could introduce interference.In order to explore the potential coupling effect during the penetration of dense particles into gelatin, simulations were conducted involving two particles with equal diameters and velocities of 300 m/s, 500 m/s, 700 m/s, and 900 m/s. Figure 17 shows typical moments of double-particle penetration into gelatin, with particle sizes of 0.3mm and velocities of 500m/s.In the simulations, it was observed that the gelatin in the central region of the cavity was subjected to fractures and breakages, which were attributed to the tension and collision of stress waves.These effects were manifested as element failures and deletions within the numerical simulation.It could be known that when the initial spacing between the two particles was large, the interference behavior would not occur.Then there had to be a case in which the cavity generated by the two particles did not interfere.This space could be named the critical cavity interference distance Δy.Ideally, the critical cavity interference distance should be exactly the size of the largest temporary cavity diameter of a single particle, that is Δy=d max .Due to the transmission of a stress wave in gelatin, the expansion environment of adjacent cavities changes, which may make the critical cavity interference distance different from the maximum temporary cavity diameter.Therefore, numerical simulations were carried out for different particle spacing situations.
By comparing the penetration results of a single particle and double particles, it was found that the penetration depth of the double particle was basically consistent with that of the single particle and had nothing to do with the particle spacing, as shown in Figure 18.Submillimeter particle cloud has no coupling effect on the penetration depth of gelatin.
It could be seen that as the distance between particles increased, the process of a temporary cavity from interference to non-interference is shown in Figure 19.The numerical simulation results showed that the critical cavity interference distance was smaller than the maximum temporary cavity diameter, Δy < d max , as shown in Table 7.A possible explanation for the phenomenon described above can be offered.During particle penetration, an extremely high pressure is generated in the region surrounding the particle contact surface.This pressure propagates both radially and axially, as illustrated in Figure 20 (a).The shock pressure wave on one side exerts a force that deforms the gelatin within the cavity on the opposite side, subsequently causing the gelatin adjacent to the cavity to be displaced into the cavity.This expansion and compression between the two particles result in a narrowing of the cavity, leading to the observed effect, as depicted in Figure 20 (b).As indicated in the illustration, the temporary cavity exhibits asymmetry, which results in a critical cavity interference distance smaller than the maximum diameter of the temporary cavity during single-particle penetration.

Figure 2
Figure 2 Three types of charged particle rings: (a) schematic diagram and (b) images.

Figure 3
Figure 3 Explosion experiment set-up: (a) top view, (b) front view, and (c) actual explosion experiment field.

Figure 5
Figure 5The change in particle cloud radius and boundary velocity over time, where the dispersed data points represent measurements from high-speed recordings.The black curve illustrates the variation of the particle cloud radius over time, while the red curve depicts the change in boundary velocity over time.

Figure 6 Figure 7
Figure 6 Results on the gelatin surface after the impact of the particle cloud.(Black circles represent burns caused by fragments, while red circles denote carbon fiber fragments.The concentrated black dots within the white box indicate the points where particles have penetrated the gelatin.)

Figure 8
Figure 8 Number density by particle recovery at different distances (the average particle mass was approximately 1.13×10-4 g).

Figure 9 18 Figure 10
Figure 9 Geometric relationship of particle cloud dispersion: (a) Top view and (b) Front view

Figure 11
Figure 11 Gelatin and particle finite element model: (a) ballistic gelatin zone division and (b) discretization and transition meshes

Figure 14
Figure 14 Typical moments of particle penetration into gelatin (particle size is 0.3mm and velocitiy is 500m/s)

Figure 15
Figure 15 The fitting curve of the simulation results: (a) Maximum penetration depth and (b) Maximum temporary cavity diameter

Figure 16
Figure 16 Comparison of numerical simulation and experimental penetration results

Figure 17
Figure 17Typical moments of double-particle penetration into gelatin (particle sizes are 0.3mm and velocities are 500m/s)

Figure 18 18 Figure 19
Figure 18 Comparison of double-particle penetration depth and single-particle penetration depth

Figure 20
Figure 20 Gelatin in double-particle penetration: (a) Pressure propagation in contact surface and (b) Temporary cavity compression

Table 1
Properties of PTFE and tungsten powder.

Table 2
Sizes and mass of particle ring filled with charges.

Table 4
Boundary velocity of the particle cloud at different detonation distances and the corresponding gelatin penetration depths.

Table 6
Material parameters of tungsten particles.

Table 7
Maximum temporary cavity diameter and critical cavity interference distance