# Simulation and Experimental Tests of Ballistic Impact on Composite Laminate Armor

Authors: Ali Murat Soydan; Bahadır Tunaboylu; Ahmed Galal Elsabagh; Abdul Kadir Sarı; Recep Akdeniz
Journal: Advances in Materials Science and Engineering
 (2018)

---

## Abstract

This paper presents the experimental testing and simulation results of ballistic impact tests on laminated armor samples that consist of three layers of different materials: fiber-cement, Kevlar fabric, and steel. In experimental tests, a 9 mm FMJ bullet was launched towards a 100 cm2 sample of the armor from the fiber-cement side. Ansys Workbench Explicit Dynamics and Ansys AUTODYN 3D were used to model and simulate the ballistic impact. Experimental testing and simulation results were compared to analyze the behavior of composite armor designs, and a good agreement was observed.


---

## Body

## 1. Introduction

The purpose of armor systems is preventing entry of bullets into the protected area. The armor to be used should be durable and easy to produce, use, and repair. The traditional method of armoring is the use of thick steel plates. However, multilayer armor was developed in decades and improved penetration resistance/weight ratio relative to steel. In multilayered armors, the outer ceramic layer deflects the bullets, the artificial fibers hold the bullet particles, and the metal part stops the bullets. This kind of armor systems significantly less weigh than one-layer solid steel armors providing equivalent protection [1]. Enhancing the ballistic resistance of armors by using composite materials is a generally accepted practice in armor design. In the past few decades, laminate armors have been investigated to achieve good durability and lightweight armors [2]. Kevlar and steel are commonly used materials in ballistic armor industry especially in structural armor types [3, 4]. Armor design is typically performed in two stages especially for those that contain multiple stacks for maximum efficacy. First steps of this stage are the armor design and selection of appropriate materials for the application and the threat under consideration. Second stage includes simulation and testing the structure which define the ballistic limits of the materials being used in designing the armor and total efficiency of the design [5].According to EN 1063 Ballistic Standards in Table1, there are different types of projectiles and testing conditions of the ballistic test. Tests are classified into levels listed from BR1 to BR7 based on the geometry, the tip shape, the shooting distance, the mass, the speed, and the material of the projectile.Table 1
Ballistic standards [6].
CEN (Committee European Normalization) BS/EN 1063 Ballistic Standards
Level
Gun Type
Core Type
Caliber
Mass (gr)
Shooting distance (m)
Bullet speed (m/s)
BR1
Riffle
0.22 LR
L/RN
2.6 ± 0.1
10.00 ± 0.5
360 ± 10
BR2
Hand gun
9 mm parabellum
FJ1/RN/SC
8.0 ± 0.1
5.00 ± 0.5
400 ± 10
BR3
Hand gun
0.357 magnum
FJ1/CB/SC
10.2 ± 0.1
5.00 ± 0.5
430 ± 10
BR4
Hand gun
0.44 magnum
FJ2/FN/SC
15.6 ± 0.1
5.00 ± 0.5
440 ± 10
BR4+
Riffle
7.62 × 39 mm
FJ2/PB/HC1
8.0 ± 0.1
10.00 ± 0.5
710 ± 10
BR5
Riffle
5.56 × 45 mm
FJ2/PB/SCP1
4.0 ± 0.1
10.00 ± 0.5
950 ± 10
BR6
Riffle
7.62 × 51 mm
FJ1/PB/SC
9.5 ± 0.1
10.00 ± 0.5
830 ± 10
BR6+
Riffle
5.56 × 45 mm
FJ/PB/SC
3.6 ± 0.1
10.00 ± 0.5
980 ± 10
BR7
Riffle
7.62 × 51 mm
FJ2/PB/HC1
9.8 ± 0.1
10.00 ± 0.5
820 ± 10
L: lead; CB: conic bullet; FJ: full steel jacket; FN: flat nose bullet; RN: rounded nose bullet; HC1: steel hard core with incendiary; SCP1: soft core with steel penetrator; PB: pointed bullet; SC: soft core.There are two major applications of ballistic armors, the body armor and the structure armors. In the battlefield, humans are usually encountered BR2 level of threat, while vehicles and buildings are encountered a higher level of threat. Material type, geometry, application, and working conditions are the key parameters of the armor design. In the simulation process, the modeling conditions of both target and projectile are playing an important role in the results. Meshing tool, mechanical properties, and physical properties of the materials and material modeling characterizing like the type of elasticity and failure conditions are the most important modeling conditions in ballistic armor simulation process [6, 7].Commercial software tools like SolidWorks, CATIA V5, and Ansys design modeler are commonly used to model the geometry of both armor and projectile. While tools like ABAQUS, Ansys with LS-DYNA, and AUTODYN 3D are used for material and test modeling [8, 9].In this study, composite-layered laminate armor structures were designed, and test samples were manufactured. Simulations, as well as experimental tests at BR2 ballistic level, were performed, and results are presented.


## 2. Experimental Conditions and Results

### 2.1. Armor Structure

The composite structures consisted of layers of Kevlar 29 fabric of 200 g/m2, 1 mm, 1.5 mm, or 3 mm thick S235JR metal sheets and a fiber-cement layer. HEXION MGS L285 resin impregnated with MGS H160 hardener was used for bonding the materials assembly. The damage in composite laminates at BR2 level ballistic impact with shear wave velocity in the target of 2 mm/µsec is localized within a zone of 3–10 cm near the impact area; therefore, samples with 10 cm × 10 cm were found to be sufficient to test the ballistic behavior of the laminate structure [10]. The samples were performed by pressing with 10–20 bar pressure at 60°C for 10 minutes with a hydraulic press.The structure has a targeted special design of fiber-cement board bonded to a steel plate via Kevlar, and an optimization study was carried out. The composite laminates were produced for add-on armored protection in prefabricated structures or buildings where the fiber-cement layer plays as a front layer which look alike the structure wanted to be protected, Kevlar layer as an aid layer to slow down the bullet, and steel layer to stop the bullet. Add-on systems will make it easier to repair structures in areas under attack. Only deformed plates will be removed and replaced with new ones. To provide protection with the cheapest and most common materials in the market, S235JR steel sheets, Kevlar 29, cellulose fibers, and fiber-cement were preferred. It was aimed to provide maximum protection at BR2 level with minimum weight, so fiber-cement and Kevlar fabrics were kept same for all samples and minimum metal sheet thickness for maximum protection was tried to be found.


### 2.2. Test Conditions

The produced composite armor samples were tested at the BR2 ballistic level against Parabellum 9 mm round tip, soft core bullet projectile which was shot towards the fiber-cement layer at the ballistic laboratories of Makina Kimya Enstitusu (MKE) in Ankara. The bullets were placed in the machine shown in Figure1, and shooting was made at 5 m distance to the composite armor specimens mounted in the center of 50 × 84 × 3 mm (Figure 2) metal board. The middle of the metal plate was cut to size 10 × 10 cm so as not to interfere with the armor feature. The produced armor is attached to the metal plate with the help of Z-shaped clamps (Figure 2). The left of the clamps is stationary, but the right side can be adjusted at the desired distance on the move and can be fixed in the position specified after the armor is placed. So that when the bullet hits the armor after the ignition, the firmly positioned armor does not move and keeps the impact entirely on itself.Figure 1
Images of test facility and ballistic testing mechanism.
(a)
(b)Figure 2
An image of the metal board (50 × 84 × 3 mm) and Z-shaped clamps.
(a)
(b)


## 2.1. Armor Structure

The composite structures consisted of layers of Kevlar 29 fabric of 200 g/m2, 1 mm, 1.5 mm, or 3 mm thick S235JR metal sheets and a fiber-cement layer. HEXION MGS L285 resin impregnated with MGS H160 hardener was used for bonding the materials assembly. The damage in composite laminates at BR2 level ballistic impact with shear wave velocity in the target of 2 mm/µsec is localized within a zone of 3–10 cm near the impact area; therefore, samples with 10 cm × 10 cm were found to be sufficient to test the ballistic behavior of the laminate structure [10]. The samples were performed by pressing with 10–20 bar pressure at 60°C for 10 minutes with a hydraulic press.The structure has a targeted special design of fiber-cement board bonded to a steel plate via Kevlar, and an optimization study was carried out. The composite laminates were produced for add-on armored protection in prefabricated structures or buildings where the fiber-cement layer plays as a front layer which look alike the structure wanted to be protected, Kevlar layer as an aid layer to slow down the bullet, and steel layer to stop the bullet. Add-on systems will make it easier to repair structures in areas under attack. Only deformed plates will be removed and replaced with new ones. To provide protection with the cheapest and most common materials in the market, S235JR steel sheets, Kevlar 29, cellulose fibers, and fiber-cement were preferred. It was aimed to provide maximum protection at BR2 level with minimum weight, so fiber-cement and Kevlar fabrics were kept same for all samples and minimum metal sheet thickness for maximum protection was tried to be found.


## 2.2. Test Conditions

The produced composite armor samples were tested at the BR2 ballistic level against Parabellum 9 mm round tip, soft core bullet projectile which was shot towards the fiber-cement layer at the ballistic laboratories of Makina Kimya Enstitusu (MKE) in Ankara. The bullets were placed in the machine shown in Figure1, and shooting was made at 5 m distance to the composite armor specimens mounted in the center of 50 × 84 × 3 mm (Figure 2) metal board. The middle of the metal plate was cut to size 10 × 10 cm so as not to interfere with the armor feature. The produced armor is attached to the metal plate with the help of Z-shaped clamps (Figure 2). The left of the clamps is stationary, but the right side can be adjusted at the desired distance on the move and can be fixed in the position specified after the armor is placed. So that when the bullet hits the armor after the ignition, the firmly positioned armor does not move and keeps the impact entirely on itself.Figure 1
Images of test facility and ballistic testing mechanism.
(a)
(b)Figure 2
An image of the metal board (50 × 84 × 3 mm) and Z-shaped clamps.
(a)
(b)


## 3. Modeling and Impact Simulation

### 3.1. Modeling Methodology

The armor has been modeled with Ansys Workbench v18.2 and 19.00 [5]; first, the geometry of both the projectile and the armor has been modeled with Ansys Design Modeler, and then a separate simulation of the projectile against each layer has been done. Finally, a simulation of the whole system was performed.Ansys Explicit Dynamics with AUTODYN has been used to solve the model and generate modeling results. The materials of each part have been modeled into Ansys engineering data, and then the model properties like mesh and contact conditions between parts have been setup. Finally, the solver properties like initial conditions, system statics, and dynamic properties, and desired output have been defined.


### 3.2. Projectile

The projectile geometry was modeled as two parts. The inner part is the bullet core, and outer part is the bullet jacket. According to the European EN 1063 Ballistic Standards in Table1, the assigned material to the core is the lead and the jacket is full metal (FMJ). In the simulation, the lead and the brass in the explicit material database of Ansys were assigned to the core and the jacket, respectively. Figure 3 shows the geometry and the 3D mesh of the bullet created using Ansys Explicit Dynamics mesh modeler with size multizone method and a hex-tetra element for the jacket and hex element for the core all with medium relevance center, coarse span angle center, and default element size. Both the brass and the lead have been modeled as explicit material [6]. The mechanical properties of both of these materials are shown in Tables 2 and 3, respectively.Figure 3
The representation in Ansys of the geometry and 3D mesh of the bullet.
(a)
(b)Table 2
Material modeling of bullet jacket brass [11].
Brass shock EOS linear
Response
Constants
Gruneisen coefficient
Parameter C1 (m·ts−1)
Parameter S1
Parameter quadratic S2 (s·m−1)
Shear modulus (Pa)
Density (kg·m−3)
2.04
3726
1.434
0
4.6e+ 010
8450Table 3
Material modeling of bullet core lead [12].
Lead shock EOS linear
Gruneisen coefficient
Parameter C1 (m·s−1)
Parameter S1
Parameter quadratic S2 (s·m−1)
2.74
2006
1.429
0
Constant and response
Density (kg·m−3)
Specific heat constant pressure J (kg−1·C−1)
Shear modulus (Pa)
11340
124
8.6e+ 009
Steinberg Guinan strength
Initial yield stress Y (Pa)
Maximum yield stress Ymax (Pa)
Hardening constant B
Hardening exponent n
Derivative (dG/dP G’P)
Derivative dG/dT G’T Pa·C−1
8e+ 006
1e+ 008
110
0.52
1
−9.976e+ 006
Derivative (dY/dP Y'P)
Melting temperature Tmelt (C°)
9.304e − 004
486.85


### 3.3. Ballistic Armor

The armor is consisting of three layers of different materials, the fiber-cement, Kevlar woven fabrics, and the steel. Each material had to be modeled separately into Ansys material engineering data. The experimental test has been done for different samples as shown in Section2. It was chosen to model two samples of different thicknesses for each layer to investigate the ballistic behavior of the whole system because the exact measurements of deformation and stress of the armor were difficult to obtain [6, 13]. The geometry of the stack has been modeled into Ansys Design Modeler as solid materials with a different thickness of each layer. Samples had x, y dimensions of 10 cm × 10 cm, and the thickness varied as follows: Sample 1 has the thickness dimensions of 8 mm of fiber-cement, ten layers of Kevlar fabric, and 3 mm thickness of the steel layer. Sample 2 has the thickness dimensions of 8 mm of fiber-cement, ten layers of Kevlar fabric, and 1.5 mm thickness of the steel layer. Sample 3 has thickness dimensions of 8 mm of fiber-cement, ten layers of Kevlar fabric, and 1 mm thickness of steel layer. The Kevlar layers were modeled with macrohomogeneous model that considers the whole layers are homogenous in geometry with orthotropic mechanical properties. This model is commonly used in modeling Kevlar material which gives sufficient result with limited resources [14]. Figure 4 shows sample 2 and the bullet as modeled on Ansys AUTODYN, while Figure 5 shows the geometry of Kevlar 29 layers.Figure 4
Sample 2 and the firing direction as modeled on Ansys AUTODYN.Figure 5
Geometry of 10 layers of Kevlar 29.


### 3.4. Fiber-Cement

After understanding that asbestos is harmful to human health, fiber-cement sheets are made of cement mix of various natural cellulose fibers. Fiber-cement is convenient for all construction systems, easy to use and install.In all samples, the fiber-cement layer had the same thickness of 8 mm. In this study, it was modeled as a composite material with linear isotropic elasticity and bilinear isotropic hardening. The mechanical properties of the material have been calculated by doing 3-point tensile test [15]. Table 4 shows the modeling properties of fiber-cement layer. The mesh for the layer was created using Ansys Explicit dynamic mesh modeler as a body mesh with default element size and high smoothing. This was also applied to the other layers of the armor. Figure 4 shows the 2D mesh of the armor as a stack.Table 4
Material modeling of fiber-cement layer.
Isotropic elasticity
Temperature (°C)
Young’s modulus (Pa)
Poisson’s ratio
Bulk modulus (Pa)
Shear modulus (Pa)
25
7.62e+ 009
0.1
3.175e+ 009
3.4636e+ 009
Constants
Failure response
Density (kg·m−3)
Tensile ultimate strength (Pa)
1390
2.9e+ 009
Bilinear isotropic hardening
Yield strength (Pa)
Tangent modulus (Pa)
6e+ 006
1.2e+ 007


### 3.5. Kevlar Fabric

One of the materials used in the experimental test was Kevlar woven fabric whose mechanical properties are dependent on the geometry of yarn and the thickness of the layer [16]. It is possible to classify the fabrics used in the ballistic applications as knitted and woven fabrics. The woven fabrics are generally preferred in military areas due to their complex production methods and costs [17]. In this study, fabric direction at an angle of ± 45° and 2D plain woven Kevlar fabric were used. Different simulations have been performed to conduct the modeling conditions of Kevlar woven fabric. In the literature, it was found that Kevlar 29, Kevlar 149, and Kevlar 49 were the most relevant materials tested as Kevlar in terms of mechanical and ballistic properties [6, 13, 14]. In stack simulation, Kevlar 29 was modeled as orthotropic elasticity mechanical properties following the MAT_054/055 composite material standard [14] which also includes fiber and matrix failure conditions [6, 14]. Figure 5 shows the geometry of 10 layers of Kevlar 29 used in the study, and the impregnated final stack can have a total thickness of 2.4 mm. Table 5, on the other hand, shows its mechanical properties.Table 5
Mechanical properties of Kevlar 29.XY axes lie in the plane of the layers whereas Z is perpendicular to the layer [14].
Orthotropic elasticity
Temperature (°C)
Young’s modulusX direction (Pa)
Young’s modulusY direction (Pa)
Young’s modulusZ direction (Pa)
Poisson’s ratioXY
Poisson’s ratioYZ
Poisson’s ratioXZ
Shear modulusXY (Pa)
Shear modulusYZ (Pa)
Shear modulusXZ (Pa)
25
1.85e+ 010
1.85e+ 010
6e+ 009
0.25
0.33
0.33
7.7e+ 008
5.43e+ 009
5.43e+ 009
Constant and response
Density (kg·m−3)
Specific heat constant pressure J·kg−1·C−1
1440
1420
Yield strengths
Compressive yield strength (Pa)
Tensile yield strength (Pa)
1.85e+ 008
1.85e+ 009


### 3.6. Steel

Soft structural steel is easy to weld and bend. In the experimental test, S235JR steel was used as the third layer of the stack. In the simulation steel, S235JR was modeled as an explicit material with grade code SAE 1006 from Ansys Explicit material data sources which indicates that it has the same mechanical and ballistic properties of S235JR steel [18]. Table 6 shows the mechanical and ballistic properties of steel 1006. Stack Mesh was illustrated in Figure 6.Table 6
Mechanical and ballistic properties of steel 1006.
Constant and response
Density (kg·m−3)
Specific heat constant pressure (J·kg−1·C−1)
7896
452
Shock EOS linear
Gruneisen coefficient
Parameter C1 (m·s−1)
Parameter S1
Parameter quadratic S2 (s·m−1)
2.17
4569
1.49
0
Shear modulus (Pa)
Initial yield stress Pa
Hardening constant Pa
Hardening exponent
Strain rate constant
Thermal softening exponent
Melting temperature C
Reference strain rate (sec)
3.5e+ 008
2.75e+ 008
0.36
2.2e− 002
1
1537.8
1
8.18e+ 010
Source: Ansys explicit material database.Figure 6
Stack mesh.


### 3.7. Analysis Settings

The model has been solved in Ansys Explicit dynamics in Ansys workbench v18.2 and 19 using Ansys AUTODYN 3D solver [5, 6]. The initial condition of the test has been simulated by assigning an initial component velocity to the bullet parts equals to 400 m/s in the Z direction and fixed support boundary condition at X and Y faces of the stack. The end time of the analysis sets to 7e − 004 s with maximum number of cycles equals to 1e+ 07. Bonded contact was added between the fiber-cement layer and Kevlar layer as well as between the Kevlar layer and the steel layer in the Ansys model simulation settings to simulate the contact condition between the layers of the samples. In the literature, there are several detailed descriptions of the AUTODYN 3D code as used in explicit dynamic problems [19, 20]. However, in brief, the code solves the conservation laws of mass and momentum in Eulerian or Lagrangian-form based on boundary initial conditions. User must define the mechanical properties and failure response of the materials modeled into the code, where the program calculated and presented the stresses in terms of pressure, volume, and internal energy [21]. The code also has the capability of solving different types of stress including equivalent (von Mises) stress and shear stress. To predict the ballistic behavior of the system, the solver output sets to calculate total and Z directional deformation of the stack layers. The mechanical properties and failure conditions have been illustrated in the material modeling section. In addition, von Mises stress was set to show the residual stress of the whole system.


## 3.1. Modeling Methodology

The armor has been modeled with Ansys Workbench v18.2 and 19.00 [5]; first, the geometry of both the projectile and the armor has been modeled with Ansys Design Modeler, and then a separate simulation of the projectile against each layer has been done. Finally, a simulation of the whole system was performed.Ansys Explicit Dynamics with AUTODYN has been used to solve the model and generate modeling results. The materials of each part have been modeled into Ansys engineering data, and then the model properties like mesh and contact conditions between parts have been setup. Finally, the solver properties like initial conditions, system statics, and dynamic properties, and desired output have been defined.


## 3.2. Projectile

The projectile geometry was modeled as two parts. The inner part is the bullet core, and outer part is the bullet jacket. According to the European EN 1063 Ballistic Standards in Table1, the assigned material to the core is the lead and the jacket is full metal (FMJ). In the simulation, the lead and the brass in the explicit material database of Ansys were assigned to the core and the jacket, respectively. Figure 3 shows the geometry and the 3D mesh of the bullet created using Ansys Explicit Dynamics mesh modeler with size multizone method and a hex-tetra element for the jacket and hex element for the core all with medium relevance center, coarse span angle center, and default element size. Both the brass and the lead have been modeled as explicit material [6]. The mechanical properties of both of these materials are shown in Tables 2 and 3, respectively.Figure 3
The representation in Ansys of the geometry and 3D mesh of the bullet.
(a)
(b)Table 2
Material modeling of bullet jacket brass [11].
Brass shock EOS linear
Response
Constants
Gruneisen coefficient
Parameter C1 (m·ts−1)
Parameter S1
Parameter quadratic S2 (s·m−1)
Shear modulus (Pa)
Density (kg·m−3)
2.04
3726
1.434
0
4.6e+ 010
8450Table 3
Material modeling of bullet core lead [12].
Lead shock EOS linear
Gruneisen coefficient
Parameter C1 (m·s−1)
Parameter S1
Parameter quadratic S2 (s·m−1)
2.74
2006
1.429
0
Constant and response
Density (kg·m−3)
Specific heat constant pressure J (kg−1·C−1)
Shear modulus (Pa)
11340
124
8.6e+ 009
Steinberg Guinan strength
Initial yield stress Y (Pa)
Maximum yield stress Ymax (Pa)
Hardening constant B
Hardening exponent n
Derivative (dG/dP G’P)
Derivative dG/dT G’T Pa·C−1
8e+ 006
1e+ 008
110
0.52
1
−9.976e+ 006
Derivative (dY/dP Y'P)
Melting temperature Tmelt (C°)
9.304e − 004
486.85


## 3.3. Ballistic Armor

The armor is consisting of three layers of different materials, the fiber-cement, Kevlar woven fabrics, and the steel. Each material had to be modeled separately into Ansys material engineering data. The experimental test has been done for different samples as shown in Section2. It was chosen to model two samples of different thicknesses for each layer to investigate the ballistic behavior of the whole system because the exact measurements of deformation and stress of the armor were difficult to obtain [6, 13]. The geometry of the stack has been modeled into Ansys Design Modeler as solid materials with a different thickness of each layer. Samples had x, y dimensions of 10 cm × 10 cm, and the thickness varied as follows: Sample 1 has the thickness dimensions of 8 mm of fiber-cement, ten layers of Kevlar fabric, and 3 mm thickness of the steel layer. Sample 2 has the thickness dimensions of 8 mm of fiber-cement, ten layers of Kevlar fabric, and 1.5 mm thickness of the steel layer. Sample 3 has thickness dimensions of 8 mm of fiber-cement, ten layers of Kevlar fabric, and 1 mm thickness of steel layer. The Kevlar layers were modeled with macrohomogeneous model that considers the whole layers are homogenous in geometry with orthotropic mechanical properties. This model is commonly used in modeling Kevlar material which gives sufficient result with limited resources [14]. Figure 4 shows sample 2 and the bullet as modeled on Ansys AUTODYN, while Figure 5 shows the geometry of Kevlar 29 layers.Figure 4
Sample 2 and the firing direction as modeled on Ansys AUTODYN.Figure 5
Geometry of 10 layers of Kevlar 29.


## 3.4. Fiber-Cement

After understanding that asbestos is harmful to human health, fiber-cement sheets are made of cement mix of various natural cellulose fibers. Fiber-cement is convenient for all construction systems, easy to use and install.In all samples, the fiber-cement layer had the same thickness of 8 mm. In this study, it was modeled as a composite material with linear isotropic elasticity and bilinear isotropic hardening. The mechanical properties of the material have been calculated by doing 3-point tensile test [15]. Table 4 shows the modeling properties of fiber-cement layer. The mesh for the layer was created using Ansys Explicit dynamic mesh modeler as a body mesh with default element size and high smoothing. This was also applied to the other layers of the armor. Figure 4 shows the 2D mesh of the armor as a stack.Table 4
Material modeling of fiber-cement layer.
Isotropic elasticity
Temperature (°C)
Young’s modulus (Pa)
Poisson’s ratio
Bulk modulus (Pa)
Shear modulus (Pa)
25
7.62e+ 009
0.1
3.175e+ 009
3.4636e+ 009
Constants
Failure response
Density (kg·m−3)
Tensile ultimate strength (Pa)
1390
2.9e+ 009
Bilinear isotropic hardening
Yield strength (Pa)
Tangent modulus (Pa)
6e+ 006
1.2e+ 007


## 3.5. Kevlar Fabric

One of the materials used in the experimental test was Kevlar woven fabric whose mechanical properties are dependent on the geometry of yarn and the thickness of the layer [16]. It is possible to classify the fabrics used in the ballistic applications as knitted and woven fabrics. The woven fabrics are generally preferred in military areas due to their complex production methods and costs [17]. In this study, fabric direction at an angle of ± 45° and 2D plain woven Kevlar fabric were used. Different simulations have been performed to conduct the modeling conditions of Kevlar woven fabric. In the literature, it was found that Kevlar 29, Kevlar 149, and Kevlar 49 were the most relevant materials tested as Kevlar in terms of mechanical and ballistic properties [6, 13, 14]. In stack simulation, Kevlar 29 was modeled as orthotropic elasticity mechanical properties following the MAT_054/055 composite material standard [14] which also includes fiber and matrix failure conditions [6, 14]. Figure 5 shows the geometry of 10 layers of Kevlar 29 used in the study, and the impregnated final stack can have a total thickness of 2.4 mm. Table 5, on the other hand, shows its mechanical properties.Table 5
Mechanical properties of Kevlar 29.XY axes lie in the plane of the layers whereas Z is perpendicular to the layer [14].
Orthotropic elasticity
Temperature (°C)
Young’s modulusX direction (Pa)
Young’s modulusY direction (Pa)
Young’s modulusZ direction (Pa)
Poisson’s ratioXY
Poisson’s ratioYZ
Poisson’s ratioXZ
Shear modulusXY (Pa)
Shear modulusYZ (Pa)
Shear modulusXZ (Pa)
25
1.85e+ 010
1.85e+ 010
6e+ 009
0.25
0.33
0.33
7.7e+ 008
5.43e+ 009
5.43e+ 009
Constant and response
Density (kg·m−3)
Specific heat constant pressure J·kg−1·C−1
1440
1420
Yield strengths
Compressive yield strength (Pa)
Tensile yield strength (Pa)
1.85e+ 008
1.85e+ 009


## 3.6. Steel

Soft structural steel is easy to weld and bend. In the experimental test, S235JR steel was used as the third layer of the stack. In the simulation steel, S235JR was modeled as an explicit material with grade code SAE 1006 from Ansys Explicit material data sources which indicates that it has the same mechanical and ballistic properties of S235JR steel [18]. Table 6 shows the mechanical and ballistic properties of steel 1006. Stack Mesh was illustrated in Figure 6.Table 6
Mechanical and ballistic properties of steel 1006.
Constant and response
Density (kg·m−3)
Specific heat constant pressure (J·kg−1·C−1)
7896
452
Shock EOS linear
Gruneisen coefficient
Parameter C1 (m·s−1)
Parameter S1
Parameter quadratic S2 (s·m−1)
2.17
4569
1.49
0
Shear modulus (Pa)
Initial yield stress Pa
Hardening constant Pa
Hardening exponent
Strain rate constant
Thermal softening exponent
Melting temperature C
Reference strain rate (sec)
3.5e+ 008
2.75e+ 008
0.36
2.2e− 002
1
1537.8
1
8.18e+ 010
Source: Ansys explicit material database.Figure 6
Stack mesh.


## 3.7. Analysis Settings

The model has been solved in Ansys Explicit dynamics in Ansys workbench v18.2 and 19 using Ansys AUTODYN 3D solver [5, 6]. The initial condition of the test has been simulated by assigning an initial component velocity to the bullet parts equals to 400 m/s in the Z direction and fixed support boundary condition at X and Y faces of the stack. The end time of the analysis sets to 7e − 004 s with maximum number of cycles equals to 1e+ 07. Bonded contact was added between the fiber-cement layer and Kevlar layer as well as between the Kevlar layer and the steel layer in the Ansys model simulation settings to simulate the contact condition between the layers of the samples. In the literature, there are several detailed descriptions of the AUTODYN 3D code as used in explicit dynamic problems [19, 20]. However, in brief, the code solves the conservation laws of mass and momentum in Eulerian or Lagrangian-form based on boundary initial conditions. User must define the mechanical properties and failure response of the materials modeled into the code, where the program calculated and presented the stresses in terms of pressure, volume, and internal energy [21]. The code also has the capability of solving different types of stress including equivalent (von Mises) stress and shear stress. To predict the ballistic behavior of the system, the solver output sets to calculate total and Z directional deformation of the stack layers. The mechanical properties and failure conditions have been illustrated in the material modeling section. In addition, von Mises stress was set to show the residual stress of the whole system.


## 4. Results

There are different approaches to comparing experimental and simulation ballistic test data [2]. Both the model and the experimental test could give a good indicator to the damage conditions of the armor in terms of shape and size. In addition, the simulation model could predict the total and directional deformations and shows the residual stress of the bullet and target as well.


### 4.1. Sample 1

Figure7(a) shows resulting images with a deformation comparison between the experimental (on the left-hand side) and simulation (on the right-hand side) of sample 1 which indicate that the sample was able to stop the 9 mm FMJ bullet. Additionally, Figure 7(b) shows that the bullet goes through the fiber-cement layer and causes some deformation in the middle section of the Kevlar and steel layer of the sample. Table 7 shows the measurement of the after-impact deformation conducted by Ansys Explicit dynamics v19. Finally, the residual stress of both the bullet and stack of sample 1 is shown in Figure 8.Figure 7
(a) Front and the back side of deformation of sample 1. (b) Total deformation of sample 1 as shown in experimental and simulation conditions (3D cross section of sample 1 deformation).
(a)
(b)Table 7
Directional deformation inZ direction of steel layer of sample 1 during the time of simulation.
Time (s)
Minimum (m)
Maximum (m)
Average (m)
1.1755e − 038
0
0
0
3.5025e − 005
−1.0591e − 005
2.2814e − 003
5.7685e − 005
7.0025e − 005
−6.6276e − 005
6.554e − 003
2.5768e − 004
1.0506e − 004
−7.5644e − 005
7.6125e − 003
5.2781e − 004
1.4002e − 004
−7.6373e − 005
8.1679e − 003
8.2409e − 004
1.7503e − 004
−1.2504e − 004
8.5275e − 003
1.1284e − 003
2.1003e − 004
−2.0057e − 004
8.7608e − 003
1.4199e − 003
2.4503e − 004
−1.1072e − 004
8.8981e − 003
1.6691e − 003
2.8003e − 004
0
8.9576e − 003
1.8495e − 003
3.1502e − 004
8.954e − 003
1.9437e − 003
3.5e − 004
8.8967e − 003
1.9528e − 003
3.8505e − 004
8.7837e − 003
1.8952e − 003
4.2003e − 004
8.6263e − 003
1.7947e − 003
4.5502e − 004
8.446e − 003
1.67e − 003
4.9001e − 004
8.2587e − 003
1.5314e − 003
5.2506e − 004
8.0735e − 003
1.3848e − 003
5.6005e − 004
7.8903e − 003
1.234e − 003
5.9505e − 004
7.7064e − 003
1.0826e − 003
6.3005e − 004
−1.4061e − 005
7.5225e − 003
9.3409e − 004
6.6505e − 004
−1.6921e − 004
7.3395e − 003
7.8644e − 004
7.0005e − 004
−2.722e − 004
7.1569e − 003
6.418e − 004Figure 8
3D cross section of sample 1 shows the residual stress in both the sample and the bullet.


### 4.2. Sample 2

Impact test on a sample of 8 mm thick fiber-cement, 10 layers of Kevlar fabric, and 1.5 mm of steel layer is also simulated. Figures9(a) and 9(b) show the total deformation of the sample in both the experimental (on the left-hand side) and simulation (on the right-hand side) conditions, respectively. As shown in Table 8, the deformation of the total laminate appears to be larger than that of sample 1. Figure 10 shows the residual stresses of sample 2 and the bullet.Figure 9
(a) Total deformation of sample 2 as shown in experimental and simulation conditions. (b) 3D cross section of total deformation of sample 2 and the 9 mm FMJ bullet.
(a)
(b)Table 8
Directional deformation inZ direction of steel layer of sample 1 during the time of simulation.
Time (s)
Minimum (m)
Maximum (m)
Average (m)
1.1755e − 038
0
0
0
3.5056e − 005
−1.3921e − 005
2.9322e − 003
7.9885e − 005
7.004e − 005
−7.4402e − 005
8.5653e − 003
3.6095e − 004
1.0502e − 004
−1.0713e − 004
1.0281e − 002
8.3111e − 004
1.4002e − 004
−1.0694e − 004
1.131e − 002
1.3797e − 003
1.7501e − 004
−8.5951e − 005
1.2086e − 002
1.9626e − 003
2.1002e − 004
−3.7253e − 004
1.2694e − 002
2.5367e − 003
2.4505e − 004
−4.4633e − 004
1.3006e − 002
3.0206e − 003
2.8003e − 004
−2.5738e − 004
1.3004e − 002
3.3204e − 003
3.1501e − 004
0
1.2813e − 002
3.4271e − 003
3.5005e − 004
1.2635e − 002
3.3909e − 003
3.8502e − 004
1.2606e − 002
3.3043e − 003
4.2005e − 004
1.2689e − 002
3.2335e − 003
4.5502e − 004
1.2788e − 002
3.1917e − 003
4.9006e − 004
−3.5527e − 018
1.2841e − 002
3.1658e − 003
5.2503e − 004
1.2848e − 002
3.1412e − 003
5.6006e − 004
1.2834e − 002
3.1079e − 003
5.9503e − 004
1.2818e − 002
3.0547e − 003
6.3001e − 004
1.2783e − 002
2.9853e − 003
6.6504e − 004
1.2739e − 002
2.9168e − 003
7.0001e − 004
1.2712e − 002
2.8615e − 003Figure 10
3D cross section of residual stress of sample 2 and the 9 mm FMJ bullet.


### 4.3. Sample 3

Finally, to predict the ballistic limit of the armor, a third sample was simulated. To make the comparison easier, the sample was chosen to have the same geometry of the other two samples in terms of Kevlar and the fiber-cement layer. The only change was the thickness of the steel layer which has been selected to be 1 mm. Figure11 shows the experimental (on the left-hand side) and simulation (on the right-hand side) results of the sample. Figure 12 shows a cross section of the model of the sample which clearly indicates that the sample has failed to stop the bullet and the armor structure has been penetrated by the 9 mm FMJ bullet. It was demonstrated from these analytical and test results of composite structures that a thinner layer of steel was not sufficient for these types of threat levels and a larger thickness, 1.5 mm, or higher was needed for proper protection.Figure 11
Total deformation of sample 3 as shown in experimental and simulation conditions.Figure 12
3D cross section of sample 3.


## 4.1. Sample 1

Figure7(a) shows resulting images with a deformation comparison between the experimental (on the left-hand side) and simulation (on the right-hand side) of sample 1 which indicate that the sample was able to stop the 9 mm FMJ bullet. Additionally, Figure 7(b) shows that the bullet goes through the fiber-cement layer and causes some deformation in the middle section of the Kevlar and steel layer of the sample. Table 7 shows the measurement of the after-impact deformation conducted by Ansys Explicit dynamics v19. Finally, the residual stress of both the bullet and stack of sample 1 is shown in Figure 8.Figure 7
(a) Front and the back side of deformation of sample 1. (b) Total deformation of sample 1 as shown in experimental and simulation conditions (3D cross section of sample 1 deformation).
(a)
(b)Table 7
Directional deformation inZ direction of steel layer of sample 1 during the time of simulation.
Time (s)
Minimum (m)
Maximum (m)
Average (m)
1.1755e − 038
0
0
0
3.5025e − 005
−1.0591e − 005
2.2814e − 003
5.7685e − 005
7.0025e − 005
−6.6276e − 005
6.554e − 003
2.5768e − 004
1.0506e − 004
−7.5644e − 005
7.6125e − 003
5.2781e − 004
1.4002e − 004
−7.6373e − 005
8.1679e − 003
8.2409e − 004
1.7503e − 004
−1.2504e − 004
8.5275e − 003
1.1284e − 003
2.1003e − 004
−2.0057e − 004
8.7608e − 003
1.4199e − 003
2.4503e − 004
−1.1072e − 004
8.8981e − 003
1.6691e − 003
2.8003e − 004
0
8.9576e − 003
1.8495e − 003
3.1502e − 004
8.954e − 003
1.9437e − 003
3.5e − 004
8.8967e − 003
1.9528e − 003
3.8505e − 004
8.7837e − 003
1.8952e − 003
4.2003e − 004
8.6263e − 003
1.7947e − 003
4.5502e − 004
8.446e − 003
1.67e − 003
4.9001e − 004
8.2587e − 003
1.5314e − 003
5.2506e − 004
8.0735e − 003
1.3848e − 003
5.6005e − 004
7.8903e − 003
1.234e − 003
5.9505e − 004
7.7064e − 003
1.0826e − 003
6.3005e − 004
−1.4061e − 005
7.5225e − 003
9.3409e − 004
6.6505e − 004
−1.6921e − 004
7.3395e − 003
7.8644e − 004
7.0005e − 004
−2.722e − 004
7.1569e − 003
6.418e − 004Figure 8
3D cross section of sample 1 shows the residual stress in both the sample and the bullet.


## 4.2. Sample 2

Impact test on a sample of 8 mm thick fiber-cement, 10 layers of Kevlar fabric, and 1.5 mm of steel layer is also simulated. Figures9(a) and 9(b) show the total deformation of the sample in both the experimental (on the left-hand side) and simulation (on the right-hand side) conditions, respectively. As shown in Table 8, the deformation of the total laminate appears to be larger than that of sample 1. Figure 10 shows the residual stresses of sample 2 and the bullet.Figure 9
(a) Total deformation of sample 2 as shown in experimental and simulation conditions. (b) 3D cross section of total deformation of sample 2 and the 9 mm FMJ bullet.
(a)
(b)Table 8
Directional deformation inZ direction of steel layer of sample 1 during the time of simulation.
Time (s)
Minimum (m)
Maximum (m)
Average (m)
1.1755e − 038
0
0
0
3.5056e − 005
−1.3921e − 005
2.9322e − 003
7.9885e − 005
7.004e − 005
−7.4402e − 005
8.5653e − 003
3.6095e − 004
1.0502e − 004
−1.0713e − 004
1.0281e − 002
8.3111e − 004
1.4002e − 004
−1.0694e − 004
1.131e − 002
1.3797e − 003
1.7501e − 004
−8.5951e − 005
1.2086e − 002
1.9626e − 003
2.1002e − 004
−3.7253e − 004
1.2694e − 002
2.5367e − 003
2.4505e − 004
−4.4633e − 004
1.3006e − 002
3.0206e − 003
2.8003e − 004
−2.5738e − 004
1.3004e − 002
3.3204e − 003
3.1501e − 004
0
1.2813e − 002
3.4271e − 003
3.5005e − 004
1.2635e − 002
3.3909e − 003
3.8502e − 004
1.2606e − 002
3.3043e − 003
4.2005e − 004
1.2689e − 002
3.2335e − 003
4.5502e − 004
1.2788e − 002
3.1917e − 003
4.9006e − 004
−3.5527e − 018
1.2841e − 002
3.1658e − 003
5.2503e − 004
1.2848e − 002
3.1412e − 003
5.6006e − 004
1.2834e − 002
3.1079e − 003
5.9503e − 004
1.2818e − 002
3.0547e − 003
6.3001e − 004
1.2783e − 002
2.9853e − 003
6.6504e − 004
1.2739e − 002
2.9168e − 003
7.0001e − 004
1.2712e − 002
2.8615e − 003Figure 10
3D cross section of residual stress of sample 2 and the 9 mm FMJ bullet.


## 4.3. Sample 3

Finally, to predict the ballistic limit of the armor, a third sample was simulated. To make the comparison easier, the sample was chosen to have the same geometry of the other two samples in terms of Kevlar and the fiber-cement layer. The only change was the thickness of the steel layer which has been selected to be 1 mm. Figure11 shows the experimental (on the left-hand side) and simulation (on the right-hand side) results of the sample. Figure 12 shows a cross section of the model of the sample which clearly indicates that the sample has failed to stop the bullet and the armor structure has been penetrated by the 9 mm FMJ bullet. It was demonstrated from these analytical and test results of composite structures that a thinner layer of steel was not sufficient for these types of threat levels and a larger thickness, 1.5 mm, or higher was needed for proper protection.Figure 11
Total deformation of sample 3 as shown in experimental and simulation conditions.Figure 12
3D cross section of sample 3.


## 5. Conclusion

Different materials and geometry designs were developed using composite laminate materials for ballistic structural armors. In this research, a 100 cm2 sample of a laminate plate consists of three different materials. Composite laminate structures with fiber-cement, Kevlar woven fabric, and steel layers were modeled with Ansys simulation program to investigate the technical feasibility of the armor design. Experimental testing was done on samples for validation. The result shows that a fiber-cement layer of 8 mm thickness, Kevlar 29 layer of 2.4 mm total thickness, and steel 1006 plate of 3 mm thickness can stop a 9 mm FMJ bullet with only slight deformation. Using the model, simulation can reduce expenses in the developing process of ballistic armor as it can predict the ballistic behavior and limits of the design and give significant insights through the product development process.



---

*Source: 4696143-2018-08-23.xml*