CC BY
27
NUMERICAL ANALYSIS OF A FRONTAL IMPACT OF A 12.7 mm PROJECTILE ON AN ARMOR PLATE
Milos S. Pesica, Aleksandra B. Zivkovicb, Aleksa D. AniciCc, Lazar J. Blagojevicd, Petko M. Boncheve, Predrag R. Pantovicf
a University of Kragujevac, Institute for Information Technologies -National Institute of the Republic of Serbia, Kragujevac, Republic of Serbia,
e-mail: [email protected], corresponding author, ORCID iD: https://orcid.org/0000-0002-3405-5216 b Military Technical Institute, Belgrade, Republic of Serbia, e-mail: [email protected], ORCID iD: https://orcid.org/0000-0001-7973-8590 c Agency for Testing, Stamping and Marking of Weapons, Devices and Ammunition, Kragujevac, Republic of Serbia, e-mail: [email protected],
ORCID iD: https://orcid.org/0000-0002-9490-833X d University of Kragujevac, Faculty of Engineering, Kragujevac, Republic of Serbia, e-mail: lazarblagojevicl [email protected], ORCID iD: https://orcid.org/0000-0001-6034-6888 e Bulgarian Defense Institute "Professor Tsvetan Lazarov", Sofia, Republic of Bulgaria, e-mail: [email protected], ORCID iD: https://orcid.org/0000-0002-0645-1676 f University of Kragujevac, Faculty of Engineering, Kragujevac, Republic of Serbia, e-mail: [email protected], ORCID iD: https://orcid.org/0000-0001-6811-7238
DOI: 10.5937/vojtehg70-38412; https://doi.org/10.5937/vojtehg70-38412
FIELD: Mechanical engineering, Materials ARTICLE TYPE: Original scientific paper
Abstract:
Introduction/purpose: The paper presents a numerical simulation of an impact of a 12.7 mm projectile on an armored metal plate with a velocity of 500 m/s at a distance of 900 m. Numerical simulations offer the possibility of drastically reducing the time required to obtain results in comparison to the time required for planning, organization and execution of experiments. The numerical simulation is done by variations in the thickness of the armor metal plate, specifically an armor metal plate of a thickness of 10 mm, 17 mm, 18 mm, and 23 mm. The mentioned armored plate thicknesses were chosen based on the results in order to determine the limit thickness of the
<1>
(/) _00
o"
Ô >
CM
of
UJ Cd ZD
O o
_J
<
o
X
o
LU
I—
>-
a: <
i—
CO <
-J
CD >o
X LU I—
o
o >
0
armored plate for the projectile perforation limit, as well as for complete ballistic protection.
Methods: Finite element modeling is used for analyzing stresses and deformations of the armored plates. The mentioned method calculates the impact of the projectile on the obstacle, precisely the collision of the projectile and the armor plate.
Results: For the comparative analysis, the parameters used are the values of the stress and the displacement. For each of the above-mentioned thicknesses of the armored metal plate, the values of stress and displacement during projectile impact were determined. The results of this study show how the thickness of the armor plate affects the interaction of the projectile and the armor plate.
Conclusion: If the physical and chemical characteristics of the armored plate remain unchanged, as the thickness of the armored plate increases, the possibility of projectile penetration decreases, and vice versa. This research is of essential importance because it analyzes the stresses and deformation of armor plates whose basic role is the protection of personnel and equipment from the projectile impact. In this regard, the thickness of the armored plate for semi-penetration of the projectile is determined.
Keywords: armor plate, projectile, impact, finite element modeling.
Introduction
Small-caliber bullet protection is a key concern for both military and civilian facilities, especially at distances up to 100 m. The main task is how to protect infantry from the effects of anti-materiel rifles in calibers of 10 mm to 20 mm. Modern war implies that infantry is transported by combat vehicles such as Infantry Fighting Vehicles (IFVs), also known as Mechanized Infantry Combat Vehicles (MICVs), or Mine-Resistant Ambush-Protected (MRAP) wheeled armored vehicles.
Troops transported by such vehicles are a very easy group target, and because of that, it is very important to protect troops inside vehicles from the effect of projectiles. In order to reduce the penetrability of vehicles, armored steel plates are added. Metallic armor plates are often used to protect moving and stationary platforms from a variety of projectiles. However, it is necessary to be careful, because the addition of armor plates affects the overall weight of the vehicle and reduces the mobility and passability of the vehicle.
Large deformation, erosion, high strain rate, dependent nonlinear material behavior, and fragmentation are all problems associated with high-velocity impact and projectile penetration.
The basic task of this paper is to determine the thickness of the plate that will be resistant to the impact of a projectile of 12.7 mm, thus protecting the infantry, and which will not affect the performance of the vehicle. For this study, only a frontal impact of a projectile into a plate of various thicknesses was considered.
The bullet used in this analysis is 12.7 mm and it is shown in Figure 1.
The core of the bullet is made of an alloy of copper and zinc, and the core of this bullet is the projectile used in the simulation. The ballistic characteristics of the core are presented in Table 1.
Figure 1 - Bullet 12.7 x 108 Рис. 1 - Пуля 12.7х 108 Слика 1 - Метак 12,7 х 108
Table 1 - Ballistic characteristics of the bullet core Таблица 1 - Баллистические характеристики сердечника пули Табела 1 - Балистичке карактеристике зрна метка
Projectile velocity (at a distance of 25 m) V25 805 m/s
Projectile velocity (at a distance of 300 m) V300 720 m/s
Pressure Pmax 304 MPa
Precision Rs300 10 cm
Core weight m 51,3 g
The material characteristics of the core of the bullet for the explicit dynamic analysis are given in Table 2.
Table 2 - Material characteristics of the core of the bullet for the explicit dynamic
analysis
Таблица 2 - Характеристики материала сердечника пули и для явного динамического анализа Taбела 2 - Матери^алне карактеристике зрна метка и параметри потребни за експлицитну динамичку анализу
Parameters Values
Johnson- Yield stress A [MPa] 112
Cook Proportionality coefficient B [MPa] 505
parameters Strain rate Impact parameter C 0.009
Temperature impact parameter m 1.68
Reinforcement exponent n 0.42
Melting temperature Tm [K] 1189
Room temperature Tr [K] 293
Constant ¿0 1
Johnson- Damage parameters Di 0.54
Cook D2 4.89
damage Da 3.03
parameters D4 0.014
D5 1.12
EOS Mie-Gruneisen equations M [m/s] 3667
parameters of state parameters Si 1.507
S2 0.000
Sa 0.000
Г 2.086
a 0.485
General Density р [t/mm3] 8.52E-9
parameters Young's modulus E [MPa] 110
Shear modulus G [GPa] 40
Poisson's ratio V 0.375
Specific heat Cp [J/kgK] 385
Finite element modeling
The penetration, damage, and failure mechanisms when the projectile impacts the armor plate were investigated using a computational model based on finite elements (Jena et al, 2019). Theoretical models are used in simulations with real material properties (projectiles and armor plates) to show how the projectile interacts with the armor plate. LS-DYNA (Livermore Software Technology, 2014), a commercially available finite element software, was used for finite element modeling and analysis (Mahfuz et al, 1999).
A two-dimensional finite element model of the armor plate and the bullet core was developed as shown in Figure 2.
Figure 2 - Finite element model Рис. 2 - Модель конечных элементов Слика 2 - Модел коначних елемената
The armor plate was meshed with four-node continuum hexahedral elements. Two-dimensional finite elements provide better computational performance/cost than fully integrated 3D elements. The element size was smallest in the region where the projectile impacted the armor plate and the element size was increased in regions away from the impact point. The overall finite element model had 8400 2D four-noded hexahedral elements. Contact was defined between the projectile and the armor plate with a hard contact definition for normal contact. The developed finite element model was used to investigate the penetration of the projectile through the base armor plate. The numerical calculation was performed at a distance of 900 m when the projectile velocity was 500 m/s.
Theoretical basis
Penetration is the motion process of a penetrator through an obstacle (armor plate for this study). The term, penetrator, means anything that is intended for penetration, and the obstacle is the environment that is exposed to the action of the penetrator. The study of the penetration process is of great importance both in the field of military technology and in the field of civilian application (Feng et al, 2020). Terminal ballistics is one of the basic disciplines that deals with defining the mechanisms of penetration, which significantly contributes to the optimization of the
<u _00 o"
O >
CM
of
UJ Cd ZD
O o
_J
<
o
X
o
LU
I—
>-
a: <
i—
< -j
CD >o
X LU I—
O
o >
0
design of the projectile, penetrating, and destructive action, as well as for the design of armor protection (Meng et al, 2021).
Depending on the outcome of the penetration process, there are four different cases:
Perforation - means the penetration of the entire penetrator through the obstacle (armor plate for this study), forming a regular, approximately cylindrical hole in the obstacle.
Limit perforation - represents the limit case of penetration because the hole in the obstacle is of irregular shape and a smaller area than the cross-sectional area of the penetrator, unlike the perforation, i.e. only parts of the broken penetrator pass through the hole.
Semi penetration - characterizes the stopping (jamming) of the penetrator in the obstacle or its breaking during penetration.
Ricochet - is the repulsion of the penetrator due to sliding on the surface of the obstacle if it is tilted.
The penetrating power of a penetrator is the ability to break through an obstacle. Increasing the penetrating power of the penetrator can be achieved by increasing the length and density of the penetrator, as well as by reducing its diameter. In opposition to this, the ability to resist penetration is the resistance of an obstacle. Increasing the resistance of the obstacle is achieved by increasing its thickness and density, as well as by improving the mechanical properties of the material. When considering penetration, the impact velocity, output velocity, and velocity of the ballistic limit are of greatest importance.
Impact velocity Vs (or Vo) is the instantaneous value of the penetrator line velocity at the moment of initial contact with the obstacle. It is assumed that the impact velocity vector is collinear with the penetrator axis, i.e. the flight of the penetrator with zero angle of attack is always assumed. The effects of the angular velocity of the penetrator around its own axis, in the case of gyro-stabilized penetrators, are not taken into account (Rajole et al, 2020).
Output (residual) velocity Vr is the velocity of the penetrator at the moment of passing the bottom of the penetrator through the plane determined by the rear surface of the obstacle.
The velocity of the ballistic limit is one of the basic characteristics of the penetrator-obstacle system and can be defined in several ways. Theoretically, this is the minimum value of the impact velocity at which the penetration occurs, or the maximum value of the impact velocity at which the penetration through the obstacle does not occur.
Johnson-Cook material model
The Johnson-Cook plasticity model was used to calculate the strain rate-dependent plastic deformation of the projectile core and armor plate material. Metal high-strain rate deformation has been successfully defined using the Johnson-Cook plasticity model (Wang & Shi, 2013). The effects of strain, strain rate, and adiabatic heating on flow stress are included in the Johnson-Cook plasticity model. The Johnson-Cook plasticity model is represented by Equation 1.
(1)
where A, B, C, n, and m are the material parameters determined from experimental data. The temperature is determined from equation 2.
(T - Tref)
T * =
(T - T )
\ melt ref )
(2)
where Tref is the temperature below which material shows no temperature dependence on flow stress. The strain rate is given by equation 3.
** = - O)
The initiation of damage is determined by equation 4, which gives the equivalent plastic strain at the onset of damage.
spl =
d,
-d2e(-d3v)'
1
d4 ln
v so j
(1 + d5T*)
(4)
where di, d2, d3, d4, and ds are the material damage parameters and s0 is the reference strain rate. The damage in the material is defined by
using a parameter D with a value between 0 and 1 where 0 means no damage and 1 means fully damaged material. Material failure occurs when D reaches a value of 1.
Material characteristics and initial conditions
AISI 4340 steel material characteristics were used for the armor plate and are shown in Table 3, while the projectile material characteristics were the ones from a copper alloy and are shown in Table 2.
Table 3 - Material characteristics of the armor plate for the explicit dynamic analysis Таблица 3 - Характеристики материала бронелиста и характеристики для явного динамического анализа Табела 3 - Матери^алне карактеристике балистичке плоче и параметри потребни за експлицитну динамичку анализу
Parameters Values
Johnson- Yield stress A [MPa] 792
Cook Proportionality coefficient B [MPa] 510
parameters Strain rate Impact parameter C 0.014
Temperature impact parameter m 1.03
Reinforcement exponent n 0.26
Melting temperature Tm [K] 1793
Room temperature Tr [K] 293
Constant ¿0 1
Johnson- Damage parameters Di 0.05
Cook D2 3.44
damage Da -2.12
parameters D4 0.002
D5 0.61
EOS Mie-Gruneisen equations M [m/s] 3850
parameters of state parameters Si 1.354
S2 0.000
Sa 0.000
Г 1.707
a 0.430
General Density р [t/mm3] 7.85E-9
parameters Young's modulus E [MPa] 210
Shear modulus G [GPa] 80
Poisson's ratio V 0.29
Specific heat Cp [J/kgK] 477
Results and discussion
In accordance with theoretical and practical knowledge, it is very easy to conclude that with increasing the thickness of the obstacle (it is important to mention that the same physical and chemical characteristics are maintained) the probability of achieving the effect of penetration decreases.
Within this paper, a numerical simulation of the penetration of a 12.7 mm projectile was performed for four different cases, i.e. for four different obstacle thicknesses.
Model 1
Figures 3-7 show the field of distribution of the von Misses equivalent stress for model 1. For this case, the impact projectile velocity was 500 m/s, the armor plate thickness 10 mm, and the simulation time 0.2 ms.
Figure 3 - Von Misses equivalent stress, time of analysis 0.01 ms - Model 1 Рис. 3 - Von Misses эквивалентное напряжение, время анализа 0.01 ms -
Модель 1
Слика 3 - Von Misses-ов еквивалентни напон, време анализе 0,01 ms - модел
1
Figure 4 - Von Misses equivalent stress, time of analysis 0.03 ms - Model 1 Рис. 4 - Von Misses эквивалентное напряжение, время анализа 0.03 ms -
Модель 1
Слика 4 -Von Misses-ов еквивалентни напон, време анализе 0,03 ms - модел 1
Figure 5 - Von Misses equivalent stress, time of analysis 0.05 ms - Model 1 Рис. 5 - Von Misses эквивалентное напряжение, время анализа 0.05 ms - Модель 1 Слика 5 -Von Misses-ов еквивалентни напон, време анализе 0,05 ms - модел 1
Figure 6 - Von Misses equivalent stress, time of analysis 0.07 ms - Model 1 Рис. 6 - Von Misses эквивалентное напряжение, время анализа 0.07 ms - Модель 1 Слика 6 -Von Misses-ов еквивалентни напон, време анализе 0,07 ms - модел 1
Figure 7 - Von Misses equivalent stress, time of analysis 0.1 ms - Model 1 Рис. 7 - Von Misses эквивалентное напряжение, время анализа 0.1 ms - Модель 1 Слика 7 - Von Misses-ов еквивалентни напон, време анализе 0,1 ms - модел 1
As it can be seen from the previous figures, for the case when the thickness of the armor plate is 10 mm, the penetration of the projectile occurs.
The projectile velocity after the impact and penetration is shown on the diagram in Figure 8.
Figure 8 - Velocity of the projectile in relation to time Рис. 8- Скорость снаряда в зависимости от времени Слика 8 - Брзина про^ектила у зависности од времена
From the diagram in Figure 8, it can be seen that the projectile perforates the armor plate after 0.03 ms. It can be noticed that the projectile velocity decreases between 0.03 ms and 0.1 ms, for the perforation required time. The projectile velocity after 0.1 ms is 350 m/s.
The displacement of the armor plate after the impact and perforation is shown on the diagram in Figure 9. The first displacements occur after 0.03 ms.
Figure 9 - Displacement of the armor plate in relation to time Рис. 9 - Смещение в зависимости от времени Слика 9 - Помераше у зависности од времена
Model 2
Figures 10-14 show the field of distribution of the von Misses equivalent stress for model 2. For this case, the impact projectile velocity was 500 m/s, the armor plate thickness 23 mm, and the simulation time 0.2 ms.
Figure 10 - Von Misses equivalent stress, time of analysis 0.01 ms - Model 2 Рис. 10 - Von Misses эквивалентное напряжение, время анализа 0.01 ms - Модель 2 Слика 10 - Von Misses-ов еквивалентни напон, време анализе 0,01 ms - модел 2
Figure 11 - Von Misses equivalent stress, time of analysis 0.04 ms - Model 2 Рис. 11 - Von Misses эквивалентное напряжение, время анализа 0.04 ms - Модель 2 Слика 11 - Von Misses-ов еквивалентни напон, време анализе 0,04 ms - модел 2
L
Figure 12 - Von Misses equivalent stress, time of analysis 0.06 ms - Model 2 Рис. 12 - Von Misses эквивалентное напряжение, время анализа 0.06 ms - Модель 2 Слика 12 - Von Misses-ов еквивалентни напон, време анализе 0,06 ms - модел 2
<u
о
О >
с\Г см о CM
of
Ш
a:
ZD О О
_J
<
о
X
о ш
i-
^
а: <
Figure 13 - Von Misses equivalent stress, time of analysis 0.1 ms - Model 2 Рис. 13 - Von Misses эквивалентное напряжение, время анализа 0.1 ms - Модель 2 Слика 13 - Von Misses-ов еквивалентни напон, време анализе 0,1 ms - модел 2
(Л С
о >о
X ш I-
о —>
О >
ф
Figure 14 - Von Misses equivalent stress, time of analysis 0.2 ms - Model 2 Рис. 14 - Von Misses эквивалентное напряжение, время анализа 0.2 ms - Модель 2 Слика 14 - Von Misses-ов еквивалентни напон, време анализе 0,2 ms - модел 2
As it can be seen from the previous figures, for the case when the thickness of the armor plate is 23 mm, the perforation of the projectile does not occur.
The projectile velocity after the impact and semi-penetration is shown on the diagram in Figure 15.
Time [ms] (E-03) H
O.
Figure 15 - Velocity of the projectile in relation to time ^
Рис. 15 - Скорость снаряда в зависимости от времени oi
Слика 15 - Брзина проjектила у зависности од времена ^
ч—
О
From the diagram in Figure 15, it can be seen that semi-penetration -s occurs.
The displacement of the armor plate after the impact and semi- = penetration is shown on the diagram in Figure 16. The first displacements occur after 0.03 ms. After 0.13 ms of the analysis, the maximum values of ^ the displacements are achieved. ^
Drenlarempril (Л
Time |ms] (E-03)
Figure 16 - Displacement of the armor plate in relation to time Рис. 16 - Смещение в зависимости от времени Слика 16 - Помераше у зависности од времена
Additional numerical simulations
It is of great importance to determine the maximum value of the plate thickness at which the penetration effect occurs, as well as the minimum value of the plate thickness at which the semi-penetration effect occurs, at the same impact velocity.
After presenting the results obtained by the numerical simulation of the penetration process, it is easy to conclude in which cases the projectile has enough energy to break through obstacles of certain thicknesses. In this case, an armor plate made of AISI 4340 alloy was used as an obstacle and it was determined that the 12.7 mm armor projectile at an impact velocity of 500 m/s achieves the effect of penetration on the armor plate with a thickness of 10 mm, while in the case of an armor plate with a thickness of 23 mm it achieves the effect of semi penetration, i.e. no penetration occurs.
In accordance with the previously defined models, using the same initial and boundary conditions, additional numerical simulations were performed and on that occasion, it was determined that the penetration effect is realized on up to 17 mm thick plates, and then the limit penetration effect occurs.
Model 3
Figures 17-21 show the field of distribution of the von Misses equivalent stress for model 3. For this case, the impact projectile velocity was 500 m/s, the armor plate thickness 17 mm, and the simulation time 0.2 ms.
Figure 17 - Von Misses equivalent stress, time of analysis 0.01 ms - Model 3 Рис. 17 - Von Misses эквивалентное напряжение, время анализа 0.01 ms - Модель 3 Слика 17 - Von Misses-ов еквивалентни напон, време анализе 0,01 ms - модел 3
Penetration - Armor plate AISI 4340. e=17mm, Vo=500m/s Time = 4.0999C-05 Contours ot Ettective Stress (v-m) max IP. value min=12.451S. al elem* 788 max=1634.76r at elem« 17413
L 1
effective Stress (v-mi 1.635e*03 1.473e*03 1.3l0e*03 1.148e*03 9.8586+02 8.236P+02 6 614e+02 4.991e+02 3.369e+02 1.747e*02 1.245e*01
Figure 18 - Von Misses equivalent stress, time of analysis 0.04 ms - Model 3 Рис. 18 - Von Misses эквивалентное напряжение, время анализа 0.04 ms - Модель
3
Слика 18 - Von Misses-ов еквивалентни напон, време анализе 0,04 ms - модел 3
СО
см
СП
а> аз !±
<u -t—' го
Cl
о Е
о
JD
О <13
сГ
Е Е
h-C\i
ГО ч— О -t—'
О
го Е
Figure 19 - Von Misses equivalent stress, time of analysis 0.06 ms - Model 3 Рис. 19 - Von Misses эквивалентное напряжение, время анализа 0.06 ms - Модель 3 Слика 19 - Von Misses-ов еквивалентни напон, време анализе 0,06 ms - модел 3
го
ч—
О
ГО
го о
аз Е
го -t—'
<13
<13 Q_
Figure 20 - Von Misses equivalent stress, time of analysis 0.1 ms - Model 3 Рис. 20 - Von Misses эквивалентное напряжение, время анализа 0.1 ms - Модель 3 Слика 20 - Von Misses-ов еквивалентни напон, време анализе 0,1 ms - модел 3
Figure 21- Von Misses equivalent stress, time of analysis 0.2 ms - Model 3 Рис. 21 - Von Misses эквивалентное напряжение, время анализа 0.2 ms - Модель 3 Слика 21 - Von Misses-ов еквивалентни напон, време анализе 0,2 ms - модел 3
As it can be seen from the previous figures, for the case when the thickness of the armor plate is 17 mm, the penetration of the projectile occurs.
The projectile velocity after the impact and penetration is shown on the diagram in Figure 22.
Figure 22 - Velocity of the projectile in relation to time Рис. 22 - Скорость снаряда в зависимости от времени Слика 22 - Брзина про^ектила у зависности од времена
From the diagram in Figure 22, it can be seen that the limit perforation occurs after 0.03 ms. It can be noticed that the projectile velocity decreases between 0.03 ms and 0.1 ms, for the limit perforation required time. The velocity of projectile fragments after 0.1 ms is 140 m/s.
The displacement of the armor plate after the impact and limit perforation is shown on the diagram in Figure 23. The first displacements occur after 0.03 ms.
CO CM <33
<33
аз !±
<u я a
о E
о
JD
О <13
сГ
E E
h-C\i
Я ч—
О '
О
я Е
я
ч—
О
Я
я о
<13
Е
я
<13
<13 Q_
Figure 23 - Displacement of the armor plate in relation to time Рис. 23 - Смещение в зависимости от времени Слика 23 - Помераше у зависности од времена
Model 4
Figures 24-28 show the field of distribution of the von Misses equivalent stress for model 4. For this case, the impact projectile velocity was 500 m/s, the armor plate thickness 18 mm, and the simulation time 0.2 ms.
Figure 24 - Von Misses equivalent stress, time of analysis 0.01 ms - Model 4 Рис. 24 - Von Misses эквивалентное напряжение, время анализа 0.01 ms - Модель 4 Слика 24 - Von Misses-ов еквивалентни напон, време анализе 0,01 ms - модел 4
Figure 25 - Von Misses equivalent stress, time of analysis 0.04 ms - Model 4 Рис. 25 - Von Misses эквивалентное напряжение, время анализа 0.04 ms - Модель 4 Слика 25 - Von Misses-ов еквивалентни напон, време анализе 0,04 ms - модел 4
Figure 26 - Von Misses equivalent stress, time of analysis 0.06 ms - Model 4 Рис. 26 - Von Misses эквивалентное напряжение, время анализа 0.06 ms - Модель 4 Слика 26 - Von Misses-ов еквивалентни напон, време анализе 0,06 ms - модел 4
Figure 27 - Von Misses equivalent stress, time of analysis 0.1 ms - Model 4 Рис. 27 - Von Misses эквивалентное напряжение, время анализа 0.1 ms - Модель 4 Слика 27 - Von Misses-ов еквивалентни напон, време анализе 0,1 ms - модел 4
ф
о h-
Ô >
с\Г
СМ о СМ
of
Ш
cd
ZD О О
_J
<
о
X
о ш
I—
>-
а: <
Figure 28 - Von Misses equivalent stress, time of analysis 0.2 ms - Model 4 Рис. 28- Von Misses эквивалентное напряжение, время анализа 0.2 ms - Модель 4 Слика 28 - Von Misses-ов еквивалентни напон, време анализе 0,2 ms - модел 4
(Л С
о >о
X ш I—
о —>
О >
ф
As it can be seen from the previous figures, for the case when the thickness of the armor plate is 18 mm, the armor plate is splitting but the penetration of the projectile does not occur.
The projectile velocity after the impact and semi penetration is shown on the diagram in Figure 29.
ьоо : ! !
.....4^1.....(
i
0. 35 0 1 0. 5 0
Time [ms] (E-03)
Figure 29 - Velocity of the projectile in relation to time Рис. 29 - Скорость снаряда в зависимости от времени Слика 29 - Брзина проjектила у зависности од времена
From the diagram in Figure 29, it can be seen that semi-penetration occurs. It can also be noticed that the 18 mm thickness of the armor plate does not provide complete ballistic protection.
The displacement of the armor plate after the impact and semipenetration is shown on the diagram in Figure 30. The first displacements occur after 0.03 ms. After 0.2 ms of the analysis, the maximum values of the displacements are achieved.
Figure 30 - Displacement of the armor plate in relation to time Рис. 30 - Смещение в зависимости от времени Слика 30 - Помераше у зависности од времена
Conclusion
Armored projectiles are intended to destroy armored targets. They penetrate armor plates thanks to enormous kinetic energy they have at the moment of collision with an obstacle and the great endurance of their body. Impact modeling for armor obstacles is very complex, extensive, and demanding, and the formed models in a very successful way approximate the real problem of projectile penetration.
It was determined that dynamic phenomena that occur during the process of ballistic penetration largely depend on deformation, strain rate, temperature, and pressure. In order to describe these phenomena in a correct way, it is necessary to define the models of material behavior. The Johnson-Cook material model and the material damage model proved to be the most suitable models for this study.
In this paper, a numerical simulation of the process of a 12.7 mm projectile penetration into armored plates of different thicknesses made of
CO CM CT>
CT>
oo !±
<D Я Cl
о E
о
о <1J
cT
E E
h-C\i
Я ч—
О '
О
я Е
я
ч—
О
Я
я
о ф
Е
я ф
ф Q_
<u _00 o~
O >
CM
of
UJ Cd ZD
O o
_J
<
o
X
o
LU
I—
>-
cd <
I—
< _J
CD >o
X LU I—
o
o >
0
AISI 4340 alloy was performed. In all 4 models, there is a contact between the bullet and the armor plate after 0.03 ms of the analysis. It is clear that when the thickness of the armor plate is 10 mm, there is perforation, and when the armor plate is 23 mm thick, there is semi-penetration.
In models 1 and 3, the armor plate destruction occurs. The velocity of the bullet after perforation through the armored plate in model 1 is 350 m/s, while in model 3 the velocity of the bullet fragments is 140 m/s.
In models 2 and 4, there is no destruction of the armored plate. In model 2, the semi-penetration of the bullet is after 0.13 ms, and in model 4 after 0.2 ms.
In all 4 models, the first displacements occur after 0.03 ms of the analysis.
However, what was also very important in this paper is to determine the limit values of the thickness of obstacles/armor plates in which penetration occurs.
The semi thickness of the armor plate at which the limit penetration occurs is 18mm. With a thickness of 23 mm, the armor plate deforms but withstands the impact of projectiles without splitting, which provides complete ballistic protection.
References
Feng, J., Sun, W., Wang, L., Chen, L., Xue, S. & Li, W. 2020. Terminal ballistic and static impactive loading on thick concrete target. Construction and Building Materials, 251, art.number:118899. Available at: https://doi.org/10.1016/j.conbuildmat.2020.118899.
Jena, D.P., Jena, D.K. & Kumar, S. 2019. Simulation of Bullet Penetration using Finite Element Method. In: 2019 International Conference on Range Technology (ICORT), pp.1-4. Available at:
https://doi.org/10.1109/IC0RT46471.2019.9069609.
-Livermore Software Technology. 2014. LS-DYNA, method: Explicit dynamic analysis, LS-PrePost-4.5-A New Post Processor for Use LSDYNA. Livermore, CA, USA: Livermore Software Technology.
Mahfuz, H., Zhu, Y., Haque, A., Abutalib, A., Vaidya, U., Jeelani, S., Gama, B., Gillespie, J. & Fink, B. 1999. Investigation of high-velocity impact on integral armorusing finite element method. International Journal of Impact Engineering, 24(2), pp.203-217. Available at: https://doi.org/10.1016/S0734-743X(99)00047-0.
Meng, S., Taddei, L., Lebaal, N. & Roth, S. 2021. Advances in ballistic penetrating impact simulations on thin structures using Smooth Particles Hydrodynamics: A state of the art. Thin-Walled Structures, 159, art.number:107206. Available at: https://doi.org/10.1016/j.tws.2020.107206.
Rajole, S., Ravishankar, K.S. & Kulkarni, S.M. 2020. Performance study of jute-epoxy composites/sandwiches under normal ballistic impact. Defence Technology, 16(4), pp.947-955. Available at:
https://doi.org/10.1016/j.dt.2019.11.011.
Wang X. & Shi, J. 2013. Validation of Johnson-Cook plasticity and damage model using impact experiment. International Journal of Impact Engineering, 60, pp.67-75. Available at; https://doi.org/10.1016/j.ijimpeng.2013.04.010.
ЧИСЛЕННЫМ АНАЛИЗ ЛОБОВОГО УДАРА СНАРЯДА 12,7-ММ ПО БРОНЕЛИСТУ
Милош С. Пешича, корессподент, Александра Б. Живкович6, Алекса Д. Аничичв, Лазар Й. Благоевичг, Петко М. Бончевд, Предраг Р. Пантовичг
а Крагуевацкий университет, Институт информационных технологий -Национальный институт Республики Сербия, г. Крагуевац, Республика Сербия
6 Военно-технический институт, г. Белград, Республика Сербия
в Агентство по испытаниям, маркировке и обозначению оружия, устройств и боеприпасов, г. Крагуевац, Республика Сербия
г Крагуевацкий университет, факультет инженерных наук, г. Крагуевац, Республика Сербия д Болгарский институт обороны «Профессор Цветан Лазаров», г. София, Республика Болгария
РУБРИКА ГРНТИ: 78.25.00 Вооружение и военная техника ВИД СТАТЬИ: оригинальная научная статья
Резюме:
Введение/цель: В данной статье представлено численное моделирование удара снаряда 12,7-мм по бронелисту со скоростью 500 м/с на расстоянии 900 м. Численное моделирование позволяет значительно сократить время, необходимое для получения результатов, по сравнению со временем, необходимым для планирования, организации и проведения экспериментов. Численное моделирование проводилось на пластинах разной толщины, их толщина составляла: 10 мм, 17 мм, 18 мм и 23 мм. Упомянутые толщины бронелиста были выбраны на основании полученных результатов с целью определения предельной толщины бронелиста и бронепробиваемости снаряда, а также для полной баллистической защиты.
Методы: Конечно-элементное моделирование используется для анализа напряжений и деформаций бронированных пластин при
ф (Л _00 о" ö > о2 Ol 0 01 ее" LU Cd ZD О О ^ С О X о ш н >- < 1- ^ пробитии снарядом. Упомянутый метод вычисляет удар снаряда о препятствие, а именно столкновение снаряда с бронелистом. Результаты: Для сравнительного анализа использовались параметры, представляющие значения напряжения и смещения. Для каждой из вышеупомянутых толщин бронированной стальной пластины были определены значения напряжений и смещений при ударе снаряда. Результаты данного исследования показывают, как толщина бронелиста влияет на взаимодействие снаряда и броневой плиты. Выводы: Если физические и химические характеристики бронеплиты остаются неизменными, то по мере увеличения толщины бронеплиты вероятность пробития снарядом уменьшается, и наоборот. Данное исследование имеет особую значимость, поскольку в нем анализируются напряжения и деформации бронелистов, основной ролью которых является защита личного состава и оборудования от проникновения снаряда. В связи с этим определяется толщина бронелиста для предотвращения пробития снарядом. Ключевые слова: бронелист, снаряд, удар, метод конечных элементов.
2: сл с О >о 2: X ш н О —> О > НУМЕРИЧКА АНАЛИЗА ФРОНТАЛНОГ УДАРА nPOJEKT/mA 12,7 mm У ПАНЦИРНУ ПЛОЧУ Милош С. Пеши^3, аутор за преписку, Александра Б. Живкови^, Алекса Д. Аничи^в, Лазар J. Блан^еви^1", Петко М. Бончевд, Предраг Р. Пантови^г a Универзитет у Крагу]евцу, Институт за информационе технологи]е -Национални институт Републике Срби]е, Крагу]евац, Ребулика Ср6и]а 6 Во]нотехнички институт, Београд, Република Срби]а в Агенци]а за испитиване, жигосане и обележаване оруж]а, направа и муници]е, Крагу]евац, Република Срби]а г Универзитет у Крагу]евцу, Факултет инженерских наука, Крагу]евац, Република Срби]а д Бугарски институт одбране „Професор Цветан Лазаров", Софи]а, Република Бугарска ОБЛАСТ: машинско инженерство, матери]али КАТЕГОРША (ТИП) ЧЛАНКА: оригинални научни рад Сажетак: Увод: У овом раду представлена jе нумеричка симулац^а удара проjектила 12,7 mm у панцирну плочу брзином од 500 m/s на растоjа^>у од 900 т. Нумеричке симулац^е нуде могуЬност
драстичног смаъеъа времена потребног за добщаъе резултата у поре^еъу са временом потребним за планираъе, организаци]у и изво^еъе експеримената. Нумеричка симулаци}а jе ура^ена за различите деблине плоча: 10 mm, 17 mm, 18 mm и 23 mm. Поменуте деблине панцирних плоча изабране су на основу резултата, а ради одре^иваъа граничне деблине панцирне плоче за продор про]ектила, као и за потпуну балистичку заштиту.
Методе: Метода коначних елемената применена jе како би се анализирали напони и деформаци}е панцирних плоча приликом удара про]ектила. Помогу наведене методе рачуна се удар про/ектила у препреку, односно колизи}а про]ектила и панцирне плоче.
Резултати: За упоредну анализу коришЯени су параметри: вредност напона и апсолутног померак>а. За сваку од наведених деблина панцирне металне плоче одре^ене су вредности напона и апсолутног помераъа при удару про]ектила. Резултати овог истраживаъа показу}у како деблина плоче утиче на интеракци}у про]ектила и панцирне плоче.
Заклучак: Уколико физичко-хеми}ске карактеристике панцирне плоче остану непромеъене, са повеЪаъем деблине панцирне плоче смаъ^е се могуЬност пробо]а про]ектила, и обрнуто. Ово истраживаъе }е од суштинског знача}а, jер анализира напоне и деформаци}е панцирних плоча, чи}а jе основна намена заштита лудства и опреме од де}ства про]ектила. С тим у вези, одре^ена jе деблина панцирне плоче за задор про]ектила.
Клучне речи: панцирна плоча, про]ектил, удар, метода коначних елемената.
Paper received on / Дата получения работы / Датум приема чланка: 11.06.2022. Manuscript corrections submitted on / Дата получения исправленной версии работы / Датум достав^а^а исправки рукописа: 13.10.2022.
Paper accepted for publishing on / Дата окончательного согласования работы / Датум коначног прихвата^а чланка за об]ав^ива^е: 14.10.2022.
© 2022 The Authors. Published by Vojnotehnicki glasnik / Military Technical Courier (www.vtg.mod.gov.rs, втг.мо.упр.срб). This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/rs/).
© 2022 Авторы. Опубликовано в «Военно-технический вестник / Vojnotehnicki glasnik / Military Technical Courier» (www.vtg.mod.gov.rs, втг.мо.упр.срб). Данная статья в открытом доступе и распространяется в соответствии с лицензией «Creative Commons» (http://creativecommons.org/licenses/by/3.0/rs/).
© 2022 Аутори. Об]авио Во]нотехнички гласник / Vojnotehnicki glasnik / Military Technical Courier (www.vtg.mod.gov.rs, втг.мо.упр.срб). Ово jе чланак отвореног приступа и дистрибуира се у складу са Creative Commons licencom (http://creativecommons.org/licenses/by/3.0/rs/).
