Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759

www.springerlink.com/content/1738-494x(Print)/1976-3824(Online)
DOI 10.1007/s12206-015-0821-7

Investigation of bullet penetration in ballistic gelatin via

finite element simulation and experiment†
Gil Ho Yoon*, Jun Su Mo, Ki Hyun Kim, Chung Hee Yoon and Nam Hun Lim
School of Mechanical Engineering, Hanyang University, 133-791, Korea

(Manuscript Received October 8, 2014; Revised March 17, 2015; Accepted April 15, 2015)

----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Abstract

To qualitatively understand the deformation processes and damages in the human body caused by high-speed impact, we conducted
experimental and computational investigations for bullet penetration into viscoelastic ballistic gelatin blocks. Because it is difficult to
measure the strain rate-dependent material properties of viscoelastic gelatin blocks during high-speed impact, the material properties that
are indirectly defined by the stress relaxation test were used for the computational simulation. We also conducted some firing experi-
ments and analyzed the deformation processes of the structures. In particular, the passing through times and the shapes of the temporary
and permanent cavities inside the ballistic gelatin blocks were analyzed and compared. This data reveals that the employed material mod-
els, with some modifications for the FE simulation, are sufficient for predicting the high-speed impact behaviors. To investigate the
shapes of the permanent cavities and fragments made by bullets inside the gelatin blocks, two-dimensional sectional images were taken
by an industrial CT scanner and a three-dimensional CAD model was constructed based on these images.
Keywords: Wound ballistics; Ballistic gelatin; Temporary cavity; Permanent cavity; 3-D scanner; Failure model
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

influenced by many factors such as the type of bullet, rotation

1. Introduction
and traveling speeds of the bullet, and the impact angle. When
This paper conducts research into wound ballistics by ex- a bullet penetrates a body, which is a mixture of various mate-
perimental and computational approaches. Terminal ballistics rials, the deceleration of a bullet occurs due to the resistance
is the study of the phenomena that occur when a bullet strikes forces of the body. These resistance forces initiate three-
and penetrates various brittle and ductile objects. When the dimensional chaotic rotations of the bullet. As expected, the
object is a person or an animal, this is referred to as sophisti- swirling, tumbling, and travelling of a bullet inside a human
cated wound ballistics. A firm understanding of the impact of body causes serious damage. Furthermore, when a bullet hits
bullets on a human body is important for both improving the bone or organs, the resistance forces exhibit sudden changes in
survival rate and also for developing efficient ammunition and magnitude and direction; these forces cause further decelera-
bullets. To understand wound ballistics, many theoretical and tions and rotations of the bullet. In order to understand these
experimental studies have been carried out [1-8]. However, it chaotic phenomena, there have been some relevant studies
is still a difficult subject due to the chaotic behavior of high- looking into the experimental analysis of ballistic penetration
speed impacts and the various rate-dependent material behav- for many materials (e.g., gelatin blocks and soap) [1-8]. In Ref.
iors [9-12]. From a physical point of view, the most relevant [13], a ballistic experiment was conducted for double-layered
and dominant parameters may be the velocity and deceleration steel plates composed with different materials. A variety of
of a projectile as the energy of a projectile is dissipated in impact conditions, such as the order of layers, nose shape of
terms of its dynamic energy, frature energy, thermal energy, projectiles, incident angle, and ballistic velocity, were investi-
and sound energy throughout a body caused by a very short gated related to the evolution of the perforation mechanisms
impact time. Thus, this study revisits the research on bullet and failure models of the plates and projectiles. For the nose
penetration into ballistic gelatin blocks in order to understand shape of projectiles, blunt-nosed projectiles and ogival-nosed
the phenomenon experimentally and computationally. projectiles were examined. These experiments were conducted
Some relevant research into wound ballistics has been con- with an in-house experimental device for a wide range of pro-
ducted previously [1-8]. Bullet penetration phenomena are jectile velocities above 400 m/s. In Ref. [14], gelatin blocks
*
Corresponding author. Tel.: +82 2 2220 0451 were shot by spherical steel projectiles with a lower impacting
E-mail address: [email protected], gilho,[email protected] velocity (about 120 m/s). With these low speeds, a projectile
†
Recommended by Associate Editor Jun-Sik Kim
© KSME & Springer 2015
3748 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759

often does not pass through the gelatin block completely. In yawing, and tumbling once it enters the block, which is also
their research, the topic of interest was the penetration depths what happens inside of the human body. Another issue that
versus the velocities normalized by the projectile diameter and should be considered is the fact that human bodies and bullets
the threshold velocity below which passing through does not usually consist of several heterogeneous materials that exhibit
occur. To identify the effects of various conditions on the different strengths and failure modes. The failure modes ob-
penetration depth, some experiments were conducted for dif- tained by a standard elongation test are not applicable for the
ferent gelatin concentrations, different projectile diameters, high strain rate deformation case in this scenario because the
and different projectile velocities. In Ref. [15], ballistics with failure modes can be different. To cope with these difficulties,
10% ballistic gelatin blocks behind soft armor were investi- special experiments must be carried out with some empirical
gated. In their research, the bullet does not penetrate into the assumptions. Published experimental results can be utilized as
gelatin but produces a blunt impact on the gelatin block’s sur- in some countries, i.e., the Republic of Korea, and some spe-
face. Therefore, the attenuation of the shockwave could be cial licenses can be obtained by administrative governments to
monitored using five pressure sensors imbedded at different test the penetration phenomena for research purposes. Because
positions in the gelatin block. In Refs. [16, 17], the expansion gelatin blocks with different concentrations can be used, some
behavior of gelatin blocks was investigated. In these studies, material properties of the viscoelastic material model are
projectiles passed through the gelatin blocks completely. In measured for gelatin blocks with different concentrations us-
Ref. [16], both the increase in diameter and the length reduc- ing a spherical indentation test. In addition, note that by
tion of the temporary expansion cavity in the gelatin block changing the densities of the gelatin blocks (i.e., the ratio be-
were taken into account. For different gelatin concentrations, tween the gelatin powder and water), it is possible to make
lengths of gelatin blocks, and impact velocities, the deforma- gelatin blocks with material properties and fracture mecha-
tion rates were compared using a set of predefined parameters. nisms that are similar to those of a human body. Depending on
In Ref. [17], the passing through of gelatin blocks was investi- the subjects of each study, the densities studied for bullet
gated using both experiments and numerical simulation. In penetration can vary. In this study, 10% gelatin blocks are
their research, steel spheres were used as projectiles, and the made at a low temperature (between 0 ℃ and 10 ℃). Finally,
velocity range of the spherical projectiles was between 728 obtaining physically meaningful numerical simulation for the
m/s and 947 m/s. For the simulation case, projectile velocities penetration phenomena for a very short period of time inside
were increased up to 1630 m/s. In addition, the effect of the of a rate-dependent, nonlinear material is still a challenging
instant impact velocity and sphere diameter on the damage of task. In our study, we simulate ballistic gelatin blocks with
the gelatin blocks was investigated. To validate the effective- embedded bone and nerve fiber and conduct real firing tests in
ness of their computational model, the test firing findings and order to show the validity of our simulations and to under-
the computational results were compared related to the impor- stand the physical behaviors of the gelatin blocks. For the FE
tant features of the impact event such as a permanent tunnel simulation, the Equation of state (EOS), Johnson-Cook failure
formed in the gelatin, the time histories of the pressure at a model, and viscoelastic material model are employed for the
certain point, and the time histories of the temporary cavity simulation of the bullet and gelatin. Some material properties
size. published in the literature are also employed.
There are some difficulties in studying wound ballistics. The layout of this paper is organized as follows: the mate-
One of the challenging issues in understanding wound ballis- rial failure models and the material properties are briefly dis-
tics is that these events happen for a very short period of time cussed in Sec. 2; some finite element simulations are also
(i.e., about a few microseconds) and that they show chaotic conducted. In Sec. 3, gun firing tests, using the same simula-
behavior (which is very difficult to reproduce) inside opaque tion conditions of Sec. 2, are carried out and their results are
materials such as blood, skin, muscle, and bone. In other compared. The conclusion section summarizes our findings
words, it is almost impossible to see exactly what happens and presents some future research topics.
when a person is hit by a bullet. To cope with this issue, trans-
parent ballistic gelatin blocks (a pseudo-elastic material) or
2. Finite element simulation of high-speed impact of
soap (a plastic material), where the internal behavior can be
gelatin block
seen and the material properties are similar to those of a hu-
man body, are used. To capture the images of fast and large Among the four sub-areas of ballistics, i.e., interior ballis-
deformations, the use of high-speed cameras is usually tics, intermediate ballistics, exterior ballistics, and terminal
adopted. To take snapshots of a bullet travelling faster than ballistics, this research is interested in terminal ballistics,
Mach 1, several tens of thousands of shots per a second must which studies the penetration of a target (simulating a body)
be taken by a high-speed camera. The speed of a bullet was by a bullet. Compared with the times of the first three ballistic
measured by analyzing the images of a high-speed camera and sub-areas, terminal ballistics takes the longest time. Therefore,
was determined to be 300 to 400 m/s before impacting the it is challenging to computationally analyze the behavior of a
surface of the gelatin block. Due to the resistance forces inside promptly deformed structure and many numerical and ex-
the ballistic gelatin block, the bullet experiences swirling, perimental approaches have been previously developed [16-
 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759 3749

Set Initial Conditions 12

for Node/Element i:
0.01 /s
X i0,a , vi0,a , e i0,ab , s i0,ab , ... 1000 /s
1250 /s
Time Step Size 1350 /s
9

The compressive stress (MPa)

Dt (max(vin ), min(DX ijn ))
1550 /s
New Positions
X n +1
i ,a =X n
i ,a +un +1
i ,a
Strain and Strain Rate
tensors
æ ¶u n ö
e in,ab = f ç i n,a dt ÷ , 6
ç ¶X ÷
Integrate Displacements è i,b ø
n
tn ¶e i ,ab
u n +1
i ,a
n
= ò n-1 v dt
i ,a
e&i ,ab =
n
t ¶t

Integrate Velocities Mass Conservation 3

tn ¶vin,a
tn rin = ò rin -1 dt
vin,a+1 = ò n-1 a in,a dt t n-1 ¶X in, b
t

Energy Conservation
Momentum Conservation 0
1 tn ¶g n
¶vin,a æ ¶s in,ab ö ein = s nÑ × v n + rin dt 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6
r n -1 òt
0 n
a i ,a = = fç
ç ¶X n
, Fcontact ÷
÷
n-1
¶t
¶t è i ,a ø The compressive strain (mm/mm)
Equation of state
Full Stress Tensor pin = p ( r n , e n ) Fig. 2. Strain rates of 10% ballistic gelatin at high deformation rates
s in,ab = Sin,ab - pindab
(Redrawn from Ref. [23]).
Constitutive Equations
Sin,ab ( e in,ab , e&in,ab )

Fig. 1. Finite element solution procedure (see Ref. [9] for the detail Gruneisen equation of state as the equation of state and the
descriptions). Johnson-Cook failure model as the failure model. As stated,
we employed the most commonly used material properties
18]. Among the various computational codes for wound ballis- that have been published in the literature, with the exception
tics, this research adopts the AUTODYN hydrocode to solve of the rate-dependent shear moduli of gelatin blocks with dif-
the mass conservation equation, the energy conservation equa- ferent concentrations. In this research, these shear moduli are
tion, and the Equation of state (EOS), as shown in Fig. 1 [9]. measured by indentation tests or by stress-relaxation tests.
For the sake of accurate simulation, the material properties for
the ballistic gelatin and bullet should be inserted accurately; 2.1.1 Equation of state (EOS)
however, these show a large variance. Thus, in this section, the The equation of state is the equation used for the hydrody-
FE simulation and our measurements of the involved material namic prediction of a solid material. Without shear stress or
properties are described. force, the pressure inside a gas or fluid is only dependent on
the density and the internal energy. In this scenario, the equa-
tion of state is used to calculate the magnitude of the pressure.
2.1 Material model of ballistic gelatin block
When a solid material is under high-speed deformation, it
The numerical simulation of high-speed impact on the bal- behaves like a gas or fluid, and the internal pressure of the
listic gelatin, whose behavior is similar to that of human tissue solid material can be calculated with the equation of state.
without bone or skin, is a very difficult subject. The ballistic Because ballistic gelatin deforms rapidly, the Mie-Gruneisen
gelatin can be modeled as a viscoelastic material whose shear EOS is used. The Mie-Gruneisen EOS is the Taylor approxi-
modulus is rate-dependent and whose local behaviors under mation of the state surface of the measured Hugoniot curve,
high-speed impact become that of a liquid or gas. and it is known that this EOS cannot predict phase changes
In the literature, the rate-dependent stress-strain relationship [23]. With the Mie-Gruneisen EOS, the internal pressure is
has been tested in many ways [19-22], including the Hopkin- determined with respect to the volume and the internal energy
son bar test, which utilizes the velocity of sound in the mate- in the hydrocode (AUTODYN) [24].
rial, as shown in Fig. 2 [23]. Although this method is a stan-
dard measurement, conducting these experiments for a wide 2.1.2 Johnson-cook failure model
range of velocities at our institute is difficult. To address the rate-dependent failure mode, this research
For more accurate FE simulation, we should adopt a me- adopts the Johnson-Cook failure model as follows:
chanical model of the rate-dependent material properties. Ba-
sically, the material models should be accurate enough to pre- De
D=å (1)
dict the behavior of the ballistic gelatin but simple enough to ef
calculate its behavior in a reasonable amount of time. Until e f = é D1 + D2e D s ù é1 + D4 ln e& * ù [1 + D5T * ]
3
*
(2)
now, many theories and approaches have been developed; in ë ûë û
this research, the Equation of state (EOS), the Johnson-Cook D1 = -0.135490, D2 = 0.601496,
(3)
failure model, and the viscoelastic material model are em- D3 = 0.258923, D4 = 0.030127, D5 = 0
ployed for the rate-dependent bulk modulus and the rate-
dependent shear modulus [24]. In our study, we use the Mie where the effective failure strain is denoted as e f . The nor-
3750 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759

where the instantaneous shear modulus and the long-term

shear modulus are denoted as G0 and G¥ , respectively. The
decay constant is denoted as b . From the literature, it has
been shown that the involved constants are highly dependent
on the manufacturing environments [25-28]. Therefore, to
have consistent viscoelastic material properties, we made 10%
gelatin blocks. Indentation and stress relaxation experiments
were conducted.
(a) Creep experiment
2.2 Indentation test for shear modulus

The spherical indentation test measures the change of the

reaction force with respect to time when a spherical indenter is
lowered into a viscoelastic material, as shown in Fig. 4(a) [25,
29]. The rate-dependent material properties are measured for
different gelatin concentrations. In our research, we adopt the
experimental setup shown in Fig. 4(b).
(b) Stress relaxation test To measure the rate-dependent material properties, the pol-
Fig. 3. Experiments for determining the coefficients of a viscoelastic ished stainless ball probe is first lowered down to a 10% gela-
material: (a) creep response to the initial stress s 0 ; (b) the stress tin block with a certain indentation length and the probe is
relaxation test due to the initial strain e 0 (The deviatoric stress and fixed, as shown in Figs. 4(a) and 4(c). While maintaining this
the deviatoric strain are denoted as s deviatoric and e deviatoric , respec- status, the reaction forces from the gelatin block to the probe
tively.).
are measured for a sufficient period of time. A standard meas-
ured force is given in Fig. 4(d). Because the probe is lowered
malized effective plastic stress and rates are denoted as s * with a constant speed, it can be observed that the reaction
and e& * , respectively. The normalized temperature is T * . The force proportionally increases up to 50 gf (Fig. 4(d)). When
involved constants D1 , D2 , D3 , D4 , and D5 can be ob- the probe is fixed, due to the stress relaxation in the viscoelas-
tained by Hopkinson-bar or Taylor-bar tests. From the litera- tic material, the magnitude of the reaction force decreases, as
ture, these constants are set by ignoring the effect of tempera- shown in Fig. 4(d). The mathematical function to fit the sec-
ture in this research [24]. tion of the reaction force after it has been fixed is constructed
with the Eq. (5). The parameters are inserted into Eqs. (6) and
2.1.3 Rate-dependent shear modulus (7) to obtain the rate-dependent shear modulus [25-28]. In
It is known that many polymers, such as gelatin, rubber and these equations, the maximum indentation length of the probe
plastic, often behave as elastic materials, viscous materials, or is denoted asand the time when the stress relaxation starts is
a combination of the two. Thus, the ballistic gelatin consid- denoted as t R .
ered here is modeled as a viscoelastic material showing altered
t
stress and strain curves under altered loading conditions. From 3 -
ti
P (t ) = ( B0 + å i =1 Bi e ) (5)
a mechanical energy point of view, the viscoelastic material
dissipates the internal energy with respect to time. To model Bi
Ci = 3
, b 0 = 1,
the viscoelastic material properties, various finite element b i hmax
2
(8 R / 3)
models have been proposed. In our numerical simulation, the (6)
t i t /t
Kelvin model and the Maxwell model are combined, as im- bi = (e R i
- 1), i > 0
tR
plemented in AUTODYN. Note that determining the coeffi-
t
cients in the viscoelastic material model is an issue that must 1 3 -
G (t ) = (C0 + å i =1 Ci e t i ) . (7)
be solved. One of the conventional models conducts creep 2
experiments and measures the relaxed strain with respect to
time, as shown in Fig. 3(a). Another method is the stress re- After arriving at Eq. (7), we obtain the parameters used for
laxation method, which measures the stress value with respect the FE hydrocode (AUTODYN) with Eq. (8).
to time to compute the involved constants, as shown in Fig.
G (t ) = G¥ + (G0 - G¥ )e - b t (8)
3(b).
By combining the Kelvin model and the Maxwell model, whereis the instantaneous shear modulus,is the long-term
the shear modulus G becomes time dependent, as follows: shear modulus, andis the viscoelastic decay constant. The
values obtained using the above indentation test (Fig. 4) for a
G0 -G¥ 10% gelatin block are shown in Table 1. These obtained val-
G ( t ) = Gv e - b t , b = , Gv = G0 - G¥ (4)
h ues are used in the viscoelastic material model.
 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759 3751

Table 1. The instantaneous shear modulus, long-term shear modulus,

and decay constant.
Spherical indenter
G0 G¥ b
21427.41 Pa 15818.91 Pa 0.00087

Ballistic gelatin

(a)

Fig. 5. Spherical indentation test according to the aging time (10%

ballistic gelatin block).

2.2.1 Effect of preparation time and temperature increase

To prevent both ice formation and liquidation from ad-
versely affecting the material properties of the ballistic gelatin,
to ensure that its material properties are similar to that of a
human body, its temperature should be maintained between
(b) 0°C and 10°C. However, because the real firing tests are con-
ducted at normal room temperature, degradation of the mate-
rial properties of the ballistic gelatin block should be consid-
ered and controlled if possible [19]. In our preparation, the
gelatin blocks are kept inside ice boxes before transportation
to the experiment site. Experimental setup usually took one or
two hours. Even though the exposure time to the hot light
source for the high-speed camera is short, the temperature
increase caused by the heat from the light source cannot be
ignored. To show this effect on material degradation, Fig. 5
shows the recorded reaction forces at a depth of 2 mm at dif-
ferent time intervals; in Fig. 4, the reaction forces are meas-
(c) ured with 1-mm-deep indentations. The first indentation test
shows 73 gf for the maximum reaction force. After 4 h have
passed, another indentation test shows that the maximum reac-
tion force decreased to 39 gf; degradation of the material
strength can be clearly seen in subsequent indentation tests.
This experiment illustrates that the rising temperature and the
increased aging time inevitably cause degradation of the mate-
rial properties; this is one of the sources of the chaotic experi-
mental results that we obtained.

2.3 Finite element modeling and applied material properties

(d) for impact simulation
Fig. 4. Spherical indentation test: (a) diagram of spherical indentation 2.3.1 Modeling of gelatin block
test on the gelatin sample (P: reaction force on the indenter, R: radius Many hydrocodes that can simulate high-speed impact have
of spherical indenter, and h: displacement of indenter from the surface
of the gelatin sample); (b) test machine used for indentation test; (c)
been developed. Among them, this research uses the
measure of displacement of indenter with time; (d) measure of force on AUTODYN hydrocode with the above finite element theory.
the indenter with time. The finite element model is constructed with the ballistic
3752 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759

Table 2. Material properties of the 10% gelatin block.

Viscoelasticity
Density, r
Instantaneous shear modulus, G0 Shear modulus, G¥ Viscoelastic decay constant, b
1030 kg/m3 0.214 MPa 0.158 MPa 0.00087 s -1
Mie-Gruneisen EOS (Shock EOS linear)
Gruneisen coefficient c1 s1 s2
0.17 1553 m × s -1 1.93 0 s × m -1
Johnson-Cook failure model
D1 D2 D3 D4 D5 Melting temperature Reference strain rate
-0.13549 0.6015 0.25892 0.030127 0 20°C 0.001

Table 3. Material properties of the bullet interior (Lead).

Steinberg Guinan strength

Shear
Density, r Initial Maximum modulus
yield stress yield stress
1134 kg/m3 8 MPa 0.1 GPa 11.13 GPa
Shock EOS linear
Gruneisen
c1 s1 s2
coefficient
-1
2.74 2006 m × s 1.429 0 m × s -1

Fig. 6. Finite element model for ballistic gelatin block and bullet. Table 4. Material properties of the bullet coating (Copper).

30 mm Bilinear isotropic hardening

19 mm Shear
Density, r Tangent
Yield strength modulus
3.0 mm modulus
1.3 mm
15 mm 8900 kg/m3 70 MPa 2 GPa 46.4 GPa
Gruneisen
c1 s1 s2
COPPER coefficient
2 3958 m × s -1 1.497 0 m × s -1
GUNPOWDER LEAD
Ø 10 mm
Ø 9.9 mm

The material properties for the 10% gelatin block used in

Ø 9.7 mm
Ø 9.0 mm

the AUTODYN hydrocode are illustrated in Table 2. As

stated, the material properties of the ballistic gelatin block are
(a) measured by the indentation test and the other material
properties are taken from the Refs. [30-32].

2.3.2 Modeling of bullet

The sketch and finite element mesh for the 9 mm 19 mm
parabellum bullet are given in Fig. 7. The internal part of the
bullet is filled with lead and the outer surface is coated with
copper [33].
The elasto-plasticity material model is used for the bullet in
Tables 3 and 4. Strictly speaking, the bullet should also be
(b)
modeled considering the high-speed impact behavior. Because
Fig. 7. Modeling of 9 mm ´ 19 mm parabellum bullet: (a) bullet profile the rigidity of the bullet is sufficiently high, the elasto-
sketch; (b) finite element mesh. plasticity material model is employed [33, 34].

gelatin block and the bullet is modeled as a 9 mm ´ 19 mm 2.3.3 Modeling of cow bone
parabellum bullet. The size of the gelatin block is 30 cm 15 As stated in the introduction, one of our goals is to
cm 15 cm. The finite element meshes of the bullet and the investigate the mechanical behavior of bone under high-speed
gelatin block are shown in Fig. 6. impact. To our knowledge, there has been some research that
 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759 3753

Table 5. Material properties of cow bone.

Viscoelasticity
Density, Instantaneous Shear Viscoelastic Bulk
r shear modulus, modulus, decay modulus
G0 G¥ constant, b
2000 kg/m3 0.675 MPa 0.5 MPa 0.1 s -1 6 GPa
Principal strain failure
Maximum principal strain, e max Maximum shear strain, g max
(a)
0.5 0.5

(b)

Fig. 9. Experiment setup for gelatin blocks: (a) schematic illustration;

(b) real experiment.

(a)

(a) (b)

(c)
(b) Fig. 10. Gelatin blocks: (a) gelatin block without cow bone; (b) with
Fig. 8. Finite element model of gelatin block with bone: (a) position of cow bone; (c) the high-speed camera and a lamp.
bone; (b) shape of bone.

3. High-speed impact experiments

has considered bone under high-speed impacts [10]. As an
3.1 Experiment
extension, our research conducted some experiments with
gelatin blocks and bone. The bone is modeled, as shown in To study the rapid deformation in wound ballistics, this
Fig. 8, and inserted into the upper portion of the gelatin block. section presents our experimental and computational results;
Some static failure models of bone can be found in the Ref. due to the regulations of the ministry of military of the
[35], but it is rare to find studies that investigate the dynamic Republic of Korea for the ballistic experiment, only Sig Saure
failure of bone. Therefore, we inversely obtained the material P226, which is commercially available, is used for the
properties after comparing the finite element results with the following experiments (Fig. 9). A trained shooter shoots the
experimental results; the experimental results are presented in gelatin blocks from a distance of 5 m, as shown in Fig. 9.
a later portion of this paper. The cow bone is modeled with a In these experiments, 15 cm ´ 15 cm ´ 30 cm gelatin blocks
viscoelastic material model with a high stiffness value, as with and without cow bone and nerve fibers are prepared, as
shown in Table 5. shown in Figs. 10(a) and 10(b) [36]. A high-speed camera,
3754 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759

shown in Fig. 10(c), is installed at the sides of the gelatin 0 sec 0.00038 sec 0.00084 sec 0.0013 sec

blocks.
In order to capture the high-speed rate deformation before
0.0018 sec 0.0022 sec
and after impact, a high-speed camera with more than 10,000 0.0027 sec 0.0031 sec

frames per second is installed. Considering the recoding time

related to the memory and resolution of our high-speed
0.0036 sec 0.0041 sec 0.0045 sec 0.0050 sec
camera, a shutter speed of 13029 frames per second is set.
After the bullet’s impact, most of the kinetic energy of the
bullet is dissipated by deforming and tearing up the ballistic
gelatin block and we can observe both temporary and Fig. 11. High-speed camera images taken during impact (13029
permanent cavities. The evolving shapes of temporary cavities frames/ second, model: Phantom V7.3).
can be captured and investigated by the high-speed camera.
However, it is difficult to precisely examine the remaining
shapes of the permanent cavities. Normally, to examine the
shapes of a permanent cavity, the ballistic gelatin blocks are
cut after ballistic testing. Unfortunately, this inevitably distorts
the shapes of the cavities. To overcome this limitation in our
study, the shapes of permanent cavities are instead scanned by
an industrial CT (Computed tomography) and three-
dimensional CAD models of the permanent cavities are built
with the scanned CT images.

3.2 Ballistic gelatin experiment 1

For the first experimental and computational tests, the

instantaneous images of the pure ballistic gelatin blocks before
and during the ballistic collision are captured; the FE results
are compared with these images. It should be mentioned that
the quantitative analysis between the experiment and the finite
element simulation is almost impossible due to the chaotic
behavior of the bullet. However, the qualitative analysis that is
performed here and the deformation tendencies of the wounds Fig. 12. Comparisons of the deformations from the experiment and the
in the gelatin blocks are similar to each other. finite element simulation.
Investigating the images taken with our high-speed camera
(at 13029 frames per second) reveals that, before collision, the reveals that nearly-identical behavior can be predicted by the
bullet moves toward the block of ballistic gelatin at 21 mm per FE procedure; without the Johnson-Cook failure model, the
frame; this is equivalent to a bullet speed of 313 m/s. The time pulsations of the temporary cavities cannot be obtained.
for the bullet to pass through the block of ballistic gelatin is Furthermore, it is expected that the gelatin undergoes some
about 0.0012 seconds in the experiment, whereas the rotational motion due to the centrifugal force exerted on the
simulation predicted a time of 0.0017 seconds. In our ballistic gelatin by the rotation of the bullet. In this experiment,
experiment, the decelerated speed of the bullet after the it was difficult to observe this rotational motion. After the
passing through is impossible to measure because the bullet is temporary cavity reaches maximum volume, it contracts due
covered with deformed gelatin debris. It is also observed that to the restoring force of the ballistic gelatin. Qualitatively, the
the approximate incident angle between the bullet and the contraction behaviors of the temporary cavities of the
surface of the gelatin block is 3 degrees. experiment and the FE simulation are similar to each other.
Some images taken before, during, and after the impact and Yawing phenomena caused by the resistance force from the
the collision are shown in Fig. 11. Just after impact, there are ballistic gelatin block are also observed. Normally, the
temporary cavities and a permanent cavity inside the ballistic resistance forces from the ballistic gelatin block cause bending
gelatin block. The images taken as the bullet is passing and compressive stresses on the surface of the bullet, which
through the block show that the bullet experiences immediate can compress it or even cause it to fracture. However, inside
deceleration and that the ballistic gelatin block at the impact the pure ballistic gelatin block, the bullet (9 mm 19 mm
surface flows away from the bullet to form a pulsating parabellum bullet) maintains its shape and undergoes only
temporary cavity. In Fig. 12, the evolving diameters of the minor changes in form and loses no mass: these are shape-
temporary cavities in this experiment and the finite element stable bullets. It is observed that the translation and rotation
simulation are compared to each other. This comparison kinetic energies of the bullet cause the dramatic deformations
 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759 3755

(a)
Fig. 15. Permanent cavity diameter with the position in the gelatin
block obtained by CT.

Fig. 13(b).

3.3 Ballistic gelatin experiment 2

For the next study, we installed the cow bone and nerve fiber,
whose material properties are similar to those of human bone
(b) and nerve fiber, inside the 10% ballistic gelatin block. Note that
Fig. 13. (a) Permanent cavity after penetration; (b) a three-dimensional the behavior of a projectile during high velocity impact with a
CAD model based on the industrial CT images. purely viscoelastic medium is largely governed by the density
of the medium. From a physical point of view, the higher the
density of the ballistic gelatin block, the greater the pressure
acting at the nose of the bullet. This causes the swirling and
tumbling motion of the bullet. Additionally, some weaker parts
of the bullet are often squeezed out and shattered. For these
reasons, pure gelatin blocks are less than ideal because they do
not have blood vessels, nerve fibers, and/or bones embedded in
them. To overcome this limitation, we installed cow bones and
nerve fibers into the gelatin blocks and conducted both FE
simulation and experiments. The same experimental
(a)
procedures used for the pure gelatin blocks were followed.
Fig. 16 shows some images taken while the bullet was
passing through the gelatin block. Although an experienced
shooter is firing the gun, some detailed impact conditions,
such as the incident angle, the position, and the impact speed,
are different between each experiment. Furthermore, the
material properties of the bone and nerve fiber are arbitrarily
chosen; as such, the FE result exhibits a large discrepancy
with the experimental result. By analyzing the images of the
bullet, it was found that the bullet speed before impact is about
(b) (c)
313 m/s and that the measured time of the bullet’s passing
Fig. 14. (a) Side CT image; (b) front CT image; (c) top-view CT image through the gelatin block is about 0.0012 seconds. However,
of the penetration scene.
the FE simulation in Fig. 17 shows that the time for the bullet
to pass through the block is 0.0014 seconds; this discrepancy
seen inside of the ballistic gelatin block. is partially caused by the erroneous material properties of the
Fig. 13(a) shows the permanent cavity inside the ballistic bone and nerve fiber, the fracture, and the pulverization of the
gelatin block. In order to examine the shape of the permanent bone debris. Additionally, the first and the second experi-
cavity, blocks are typically dismantled and the remaining ments confirm that the motion of the bullet is influenced
shape is investigated. However, in our study, an industrial CT primarily by the density, viscosity, and flowability of the
machine is utilized to take snapshots of inside the gelatin medium it is passing through. One of the observed differences
block. With the help of automated CAD software, it was between the first and second tests is that a secondary cavity
possible to construct a three-dimensional model, as shown in appears near the bone around 0.0038 s and around 0.0046 s.
3756 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759

0 sec 0.00046 sec 0.00092 sec 0.0014 sec

0.0018 sec 0.0023 sec 0.0028 sec 0.0032 sec

0.0037 sec 0.0041 sec 0.0046 sec 0.0049 sec

Fig. 16. Penetration scenes of gelatin block with cow bone and nerve fiber taken by high-speed camera (FPS: 13029 frames/second; an impact
velocity 313 m/s).

(a)

Fig. 17. High-speed camera image and finite element simulation (b)
results.

This secondary cavity is also observed in the FE simulation.

When the bullet passes through the bone, the size of the
temporary cavity behind the bone is observed to become
smaller. This is partially due to the energy dissipation of the
bullet by the bone. Because the reaction force of the bone is
larger than that of the ballistic gelatin block, the bullet
undergoes a larger rotation, which increases the size of the
permanent cavity. (c)
After the experiment, the three-dimensional CT is also Fig. 18. (a) Permanent cavity inside the gelatin block and the cow bone
applied to construct the three-dimensional CAD model after the impact; (b) penetration scene taken by CT; (c) the fracture of
(without the destruction of the ballistic gelatin block), as the bone.
shown in Figs. 18 and 19; here, it should be noted that there
are the two permanent cavities in Fig. 18. After the first shot, the bone and the nerve due to the differences in their strengths.
the shooter fired a secondary shot. Fig. 18(c) shows the This is also an important aspect that should be considered in
fracture of the bone and nerve fiber. The investigation of the the treatment of wounds and bone infection.
bone and nerve fiber shows the sliding and separation between Fig. 20 shows the diameter variation of the permanent
 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759 3757

penetration into viscoelastic ballistic gelatin blocks. The

computational investigations are based on the explicit finite
element analysis AUTODYN. For the computational models
for the viscoelastic gelatin material, metal, and bone, rate-
dependent nonlinear material models are employed. To obtain
the viscoelastic shear modulus, this study uses an in-house
indentation test facility and measurement software. In order to
conduct qualitative studies, some experiments (following
domestic regulations) were conducted. During and after the
experiments, the ballistic gelatin blocks with and without bone
are scanned with a three-dimensional CT and 3-D CAD
(a)
models are constructed to analyze the shapes of the permanent
cavities. By comparing the computational results and the
experiment results, qualitatively similar results in terms of the
shape of the temporary and permanent cavities and the
shockwave propagations can be obtained. However, because
the experiments are statically stochastic, comparison of the
quantitative results is impossible.
To consider how the presence of bone and nerve fiber affect
the high-speed impact, a cow bone is placed inside the
ballistic gelatin block. It is observed that just after direct
impact on the bone, another shockwave at the surface of the
bone is generated. This shockwave causes another cavity, in
(b) (c)
addition to the primary cavities, along the direction of the
Fig. 19. CT images at the: (a) side; (b) front; (c) top-view of the penetration bullet’s penetration. Furthermore, the retardation becomes
scene.
serious and causes the bullet to swirl and tumble. This
shockwave also cracks and damages the nerve and bone.
Delamination, i.e. a mode of failure between the nerve and
bone, is also observed. Because the instantaneous kinetic
energy of the bullet is decreased significantly after impact
inside the ballistic gelatin block, which has a higher
retardation coefficient, the bullet dramatically swirls and
tumbles inside the ballistic gelatin block. This induces a
bigger permanent cavity that is approximately twice as large
compared to the scenario without bone impact.
This research investigates the high-speed penetration of a
bullet in a 10% ballistic gelatin block with and without bone.
Although this type of research has been conducted for many
Fig. 20. Permanent cavity diameter as a function of the position in the decades, these tests do not typically consider both the
gelatin block. computational and experimental results. Particularly, because
the failure modes of viscoelastic gelatin, brittle bones, and
cavity in the ballistic gelatin block. The first cavity collapses ductile bullets are very different from each other, a full-scale,
after a few milliseconds because some of the bullet’s energy is reliable computational model is still very difficult to build.
converted into the kinetic energy of the ballistic gelatin block. Improving the FE prediction accuracy was not attempted here;
Next, a secondary cavity is created; the size of the secondary this should be researched in the near future. In conclusion, the
cavity is smaller than that of the first cavity because less current study can be used to design new protective devices or
kinetic energy is available. This process creates so-called systems for humans or to improve current devices/systems by
cavity pulsation, which occurs between 0 cm and 10 cm in estimating their abilities to reduce the damage caused by the
Fig. 20. high-speed impacts.

4. Conclusions Acknowledgement
To understand and analyze the damage caused in the human This work was supported by the Research fund of Sur-
body due to high-speed impact, this research conducted vivability Technology Defense Research Center of Agency for
experimental and computational investigations for bullet Defense Development of Korea (No. UD120019OD).
3758 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759

Impact of steel spheres on ballistic gelatin at moderate

References
velocities, Int. J. Impact Eng., 62 (2013) 142-151.
[1] G. Seisson, D. Hebert, L. Hallo, J. M. Chevalier, F. Guillet, L. [18] I. Ionescu, J. E. Guilkey, M. Berzins, R. M. Kirby and J. A.
Berthe and M. Boustie, Penetration and cratering experi-ments Weiss, Simulation of soft tissue failure using the material point
of graphite by 0.5-mm diameter steel spheres at various impact method, J. Biomech. Eng-T Asme, 128 (6) (2006) 917-924.
velocities, Int. J. Impact Eng., 70 (2014) 14-20. [19] D. S. Cronin and C. Falzon, Characterization of 10% Balli-stic
[2] N. Ndompetelo, P. Viot, G. Dyckmans and A. Chabotier, Gelatin to Evaluate Temperature, Aging and Strain Rate Effects,
Numerical and experimental study of the impact of small caliber Exp. Mech., 51 (7) (2011) 1197-1206.
projectiles on ballistic soap, J. Phys. Iv, 134 (2006) 385-390. [20] J. Kwon and G. Subhash, Compressive strain rate sensi-tivity
[3] C. C. Shi, M. Y. Wang, J. Li and M. S. Li, A model of depth of ballistic gelatin, J. Biomech., 43 (3) (2010) 420-425.
calculation for projectile penetration into dry sand and [21] D. S. Cronin and C. Falzon, Dynamic characterization and
comparison with experiments, Int. J. Impact Eng., 73 (2014) simulation of ballistic gelain, Proceedings of the SEM Annual
112-122. Conference, Albuquerque New Mexico, USA (2009).
[4] B. Akers and A. Belmonte, Impact dynamics of a solid sphere [22] D. S. Cronin, C. P. Salisbury and C. R. Horst, High rate
falling into a viscoelastic micellar fluid, J. Non-Newton Fluid, characterization of low impedance materials using a poly-meric
135 (2-3) (2006) 97-108. split Hopkinson pressure bar, SEM Conference & Ex-position
[5] D. J. Frew, M. J. Forrestal and J. D. Cargile, The effect of on Experimental & Applied Mechanics, St. Louis, Missouri
concrete target diameter on projectile deceleration and (2006).
penetration depth, Int. J. Impact Eng., 32 (10) (2006) 1584-1594. [23] C. P. Salisbury and D. S. Cronin, Mechanical Properties of
[6] E. Sevkat, Experimental and numerical approaches for Ballistic Gelatin at High Deformation Rates, Exp. Mech., 49 (6)
estimating ballistic limit velocities of woven composite beams, (2009) 829-840.
Int. J. Impact Eng., 45 (2012) 16-27. [24] ANSYS, ANSYS AUTODYN User Mannual.
[7] M. G. Perdekamp, S. Pollak, A. Thierauf, E. Strassburger, M. [25] B. Qiang, J. Greenleaf, M. Oyen and X. M. Zhang, Estima-
Hunzinger and B. Vennemann, Experimental simulation of ting material elasticity by spherical indentation load-relaxa-tion
reentry shots using a skin-gelatine composite model, Int. J. tests on viscoelastic samples of finite thickness, Ieee T Ultrason
Legal Med., 123 (5) (2009) 419-425. Ferr, 58 (7) (2011) 1418-1429.
[8] A. Dorogoy, D. Rittel and A. Brill, Experimentation and [26] M. L. Oyen, Spherical indentation creep following ramp
modeling of inclined ballistic impact in thick polycarbonate loading, J. Mater. Res., 20 (8) (2005) 2094-2100.
plates, Int. J. Impact Eng., 38 (10) (2011) 804-814. [27] M. L. Oyen, Analytical techniques for indentation of vis-
[9] S. Hiermaier, Structures under crash and impact : continuum coelastic materials, Philos. Mag., 86 (33-35) (2006) 5625-5641.
mechanics, discretization and experimental characterization, [28] J. M. Mattice, A. G. Lau, M. L. Oyen and R. W. Kent,
Springer, New York, USA (2008). Spherical indentation load-relaxation of soft biological tissues,
[10] B. Kneubuehl, Wound ballistics: Basics and applications, J. Mater. Res., 21 (8) (2006) 2003-2010.
Springer (2011). [29] X. M. Zhang, B. Qiang and J. Greenleaf, Comparison of the
[11] S. Yang, A. Fafitis and A. Wiezel, Nonlinear impact model of surface wave method and the indentation method for mea-
a tennis racket and a ball, J. Mech. Sci. Technol., 26 (2) (2012) suring the elasticity of gelatin phantoms of different
315-321. concentrations, Ultrasonics, 51 (2) (2011) 157-164.
[12] Y. Kim, J. Yoo and M. Lee, Optimal design of spaced plates [30] H. W. Kang, J. Kim, Y. Yu and J. Oh, Photoacoustic
under hypervelocity impact, J. Mech. Sci. Technol., 26 (5) Response of Magnetic Nanoparticles to Pulsed Laser Irradia-
(2012) 1567-1575. tion, J. Korean Phys. Soc., 55 (5) (2009) 2224-2228.
[13] Y. F. Deng, W. Zhang, Y. G. Yang, L. Z. Shi and G. Wei, [31] G. R. Johson, Fracture characteristics of three metals sub-
Experimental investigation on the ballistic performance of jected to various strains, strain rates, temperatures and pressures,
double-layered plates subjected to impact by projectile of high Engineering Fracture Mechanics, 21 (1) (1985) 31-48.
strength, Int. J. Impact Eng., 70 (2014) 38-49. [32] M. Aihaiti and R. J. Hemley, Equation of state of ballistic
[14] M. V. Swain, D. C. Kieser, S. Shah and J. A. Kieser, Pro- gelatin (II), 55048-EG.1, Carnegie Institute of Washington, DC
jectile penetration into ballistic gelatin, J. Mech. Behav. (2011).
Biomed., 29 (2014) 385-392. [33] S. P. Mates and R. Rhorer, Dynamic deformation of copper-
[15] L. Liu, Y. R. Fan and W. Li, Viscoelastic shock wave in jacketed lead bullets captured by high speed digital image
ballistic gelatin behind soft body armor, J. Mech. Behav. correlation, Proceedings of the SEM Annual Conference,
Biomed., 34 (2014) 199-207. Indianapolis, Indiana, USA (2010) 431-441.
[16] F. Bresson, J. Ducouret, J. Peyre, C. Marechal, R. Delille, T. [34] D. S. Preece and V. S. Berg, Bullet impact on steel and
Colard and X. Demondion, Experimental study of the expan- Kevlar/steel armor-computer modeling and experimental data,
sion dynamic of 9 mm Parabellum hollow point projectiles in The Proceedings of the ASME Pressure Vessels and Piping
ballistic gelatin, Forensic Sci. Int., 219 (1-3) (2012) 113-118. Conference-Symposium on Structures Under Ex-treme Loading,
[17] Y. K. Wen, C. Xu, H. S. Wang, A. J. Chen and R. C. Batra, San Diego, CA, USA (2004).
 G. H. Yoon et al. / Journal of Mechanical Science and Technology 29 (9) (2015) 3747~3759 3759

[35] V. Astier, L. Thollon, P. J. Arnoux, F. Mouret and C. Brunet, Ki Hyun Kim received his B.S. degree
A finite element model of the shoulder for many applications: in mechanical engineering from Han-
trauma and orthopaedics, European HyperWorks Technol-ogy yang University in 2013.
Conference, Berlin, Germany (2007).
[36] J. Jussila, Preparing ballistic gelatine - review and proposal for
a standard method, Forensic Sci. Int., 141 (2-3) (2004) 91-98.

Gil Ho Yoon received his B.S. degree

in mechanical and aerospace engineer-
ing from Seoul National University in Chung Hee Yoon received his B.S.
1998. And he received his M.S. degree degree in mechanical engineering from
and Ph.D. in mechanical and aerospace Hanyang University in 2013.
engineering from Seoul National Uni-
versity in 2000 and 2004, respectively.
Dr. Yoon is currently an associate pro-
fessor at School of Mechanical Engineering, Hanyang Univer-
sity, Seoul, Republic of Korea.

Jun Su Mo received his B.S. degree in Nam Hun Lim received his B.S. degree
mechanical engineering from Hanyang and M.S. degree in mechanical engi-
University in 2013. neering from Hanyang University in
2012 and 2015, respectively.