Simulation and Optimization in Materials TechnologyView this Special Issue
Modelling of Ballistic Impact over a Ceramic-Metal Protection System
The behavior of ceramic-metal protection against a projectile impact is modeled. The model takes into account the mass and velocity for each stage of the phenomenon. A former model was modified considering more realistic parameters such as geometries and deformation profile. To analyze the model, simulations on different parameters have been run. The impact results of different ballistic projectiles were simulated, and the movement was plotted. In addition, a deterministic simulation on the mechanical properties of the back metal plate properties was done.
Gonçalves et al.  have developed a mathematical approach to the behavior of a two-layer system protection of ceramic-metal system against high speed impact. The scheme of the protection is presented in Figure 1. The ceramic layer erodes the projectile, and the back metal plate is responsible for absorbing the remaining energy by plastic deformation. The mathematical model separates the penetration into three different stages and proposes a deterministic equation that explains the deflection of the metal plate. The movement equations were manipulated to find a single problem for velocity and mass separately for each stage. The data and initial parameters provided by the original work  were used for the differential equations found. The metal deformity was reevaluated by deducing the deterministic equation. A more realistic shape parameter was considered for the postimpact curvature, and the volume approach was improved. Finally, the solution was plotted and compared to former results .
The model was better understood once the solutions were used to simulate different impacts over the system. Initially, properties of some ballistic projectiles were used to prevent the shock absorption. In addition, the mechanical properties of the metallic plate were modified intending to observe the effect on the final deformation.
A previous work  has focused on the study of the movement evolution of the projectile. This work has continued to analyze the penetration of the projectile searching for the understanding of the final deformation.
The set of data generated in the simulations together with the newly developed equations for impact and shock absorption were, in conclusion, observed as an advance for the understanding of high speed impact phenomena.
2. Materials and Methods
The software Maplesoft Maple V12 was used to solve the system. Gonçalves et al. work  provided all materials, constant values, and initial parameters necessary. Moreover, the projectile specifications used were from rifle bullets.
The mathematical model of the penetration process is divided into three stages. The first represents the initial impact and erosion of the head of the projectile. In this situation there is no penetration into the protection. The force against the projectile is and the variation in the geometry is given by where is the dynamic yielding of the projectile and is the effective contact area of the projectile. Also, is the density of the projectile’s material. According to Tabor , the dynamic hardness of a metal is the pressure with which it resists local indentation by a rapidly moving indenter. The first approximation for the value is a constant and not necessarily the same as the static pressure necessary to cause plastic flow. The actual value of the dynamic yield pressure strongly depends on the velocity of the impact and the projectile geometry. Tabor  also cites different methods to estimate this value. For simplification, following the original work , this value will be considered a fraction of the static hardness of the projectile material.
The impact generates a shock wave that travels through the material and reflects back cracking the ceramic plate. This shock wave is extremely fast and runs through the protection in a small fraction of a second. Theoretically, the sum of the incoming with the returning wave generates a region of high stress. This region is easily observed as a cone in the solids, which fits with the cracking region. Figure 2 demonstrates a single scheme for the initial impact and the propagating shock waves.
The shape and size of the ceramic tile can also affect the protection performance. The crack cone formation allows the necessary scattering to provide enough ceramic particles for the projectile erosion. With smaller geometry, the crack cone would not be properly formed once the constructive interference of the shock waves would not only occur longitudinally but also transversally. Those considerations were not taken into account in this model.
In the second stage the penetration starts pushing an interface projectile-ceramic with velocity . This interface is a small portion of ceramic in contact with the projectile which moves together but with different initial velocity. Thus, the force against the projectile is the same, but the difference of the velocities and gives the erosion ratio:
It is important to note that the erosion ceases when is equal to . In this moment, the interface reaches the same velocity as the projectile. With this, if there is no difference between both velocities, there is no stress applied in the projectile or interface material.
The modified hydrodynamic theory given by Tate  introduced the dynamic yielding () and the resistance against penetration () adopted for the ceramic in this case. Tate’s law is given by where is the density of the ceramic material. Tate’s equation is used as an important auxiliary equation to solve the model. This will be present in the next section of this work.
The second stage ends when . After this, the final fraction of the movement is given by where is the remaining mass of the projectile. The metal plate absorbs the final kinetic energy by plastic deformation. Simultaneously to the penetration of the projectile, the metallic base will move and will be deformed elastically. However, the energy produced by this elastic energy will be low and will be neglected . Once the first and second stage end, the main part of the projectile kinetic energy is absorbed by ceramic scattering and projectile erosion and can be computed using and the velocity obtained in the beginning of the third stage. Due to the deceleration of the bullet caused by the ceramic layer, the metallic plate may not be perforated. In this case, the secondary layer will suffer plastic deformation, and the plastic energy consumed by the plate can be expressed in terms of effective stress and strain as
It is considered a stress-strain curve given by the power law and the deflection profile of the plate given by where is the maximum deflection of the dimple shape, is the radius of the projectile, (dimensionless) is a deflection profile constant that can be determined experimentally, and is the radial distance of the impacting projectile center.
Assuming that the material will be bulged in an axis symmetric mode, the effective strain can be written in terms of the radial strain as
Considering the physical nature of most constants, it is possible to suppress the negative sign generated by the derivation demonstrated in (11). Also, the radial variation of the deflection profile is mirrored in the x-axis, which demonstrates that this ratio can be used as its own module. Now it gives
The solution for the expression is given by
After the initial impact, the movement of the interface starts deforming the metallic plate. It can be argued that, initially, the metal layer is compressed due to high pressure generated by the impact. In addition, the plate does not move significantly because of the low interface velocity. In this way, it is possible to affirm that the deflection of stainless steel plate starts at the end of the second stage of penetration.
Considering that the plastic energy absorbed by the metal plate is equal to the kinetic energy of the projectile in the end of the second stage, it is possible to write the final deflection of the plate as where is the velocity of the projectile at the end of the second stage.
3. Results and Discussion
The movement solutions were obtained divided by the stages and by mass and velocity. The equations needed to be manipulated together to find a unique solvable problem for mass and velocity for each stage. First it was necessary to find a differential equation with a unique function. For the first stage it was for the velocity and for the mass. The auxiliary Tate’s equation was used in the manipulation to formulate a second degree equation with the function as a result. Naturally, one of the results was not considered due to its particularity for being a complex solution, which does not have a physical interpretation. The written as a function of could be used to formulate the second stage equations. These sentences are given by for the velocity and for the mass. For the deflection of the metal plate, the new deflection profile was adopted considering the null derivative in : where is the maximum deflection of the dimple shape, is the radius of the projectile, (m−1) is a modified deflection profile constant, which is determined experimentally, and is the radial distance of the impacting projectile center.
Also, the dimple shape after the impact, shown by Figure 3, was then approximated to the volume of a small paraboloid given by
The newly deduced metal final deformity equation is now given by where is the modified profile constant and is the mathematical gamma function, which is an extension of the factorial function, with its argument shifted down by 1, to real and complex numbers.
In addition, for a solid two-layered system, it is necessary to consider the time needed for the generated wave to reflect in the back part of the protection and return to the initial point and, furthermore, start to crack the ceramic layer. The time needed for the wave to travel and reflect back is given by where and are, respectively, the thickness and the longitudinal velocity of sound for the ceramic material and and are the same constants for the metallic material. It is known that the shock wave travels with the longitudinal velocity of sound in the solid, and this velocity can be calculated by  where is the elastic modulus of the material.
The first results were obtained using the projectile and system data presented in the former work . Table 1 presents the projectile data. In the former work, basically two different compositions were used for the ceramic plate during the experimental tests. The composition and reference properties of the ceramic plates are shown in Tables 2 and 3, respectively.
In Figures 4(a) and 4(b), the upper curve, , corresponds to the projectile velocity and the lower curve, , shows the equivalent results obtained in a previous work [1, 2]. Figure 4(a) presents the results obtained in this work, and Figure 4(b) shows the former work  results.
In both graphs it is possible to observe the different stages of penetration in the graph. During the first stage, there is no interface velocity. In the beginning of the second stage, the interface begins to move tending to reach the same velocity as the projectile. The third and final stage begins with the projectile and interface velocities equalized. During this final fraction of the movement, both projectile and interface decelerate together. With the solutions it is possible to estimate the fraction of mass, velocity, and energy lost in each stage. These values are presented in Table 4.
There were more two ceramic compositions presented in the material selection. However, both of these structures were not stable enough to be molded in the system.
The values of the constants were given by the former work  together with the initial velocities and mass. The different ceramic plate thicknesses tested by the authors of the former work are also declared in the table. In addition, the first stage duration, calculated using (24), is also presented due to the different thickness and composition. Most of the constants and mechanical properties [10, 11] used in the model are presented in Table 5.
These data (21) together with the solved movement equations can give the results for the maximum deflection of the back metal plate. The comparisons are presented in Table 6. It is possible to note that the proportional errors both in the original and in the modified theory (OT and MT) are of the same magnitude. This fact is related to the simplicity of the model and the amount of neglected variable in the energy absorption, such as impact angle, ceramic porosity, and metal grain size, among others.
However, the results of the maximum deflection showed that the new theories can generate good results not only for the high speed impacts, but also for lower velocities, differently from the former theory. Simulations for different calibers  were performed using the new developed theory. The tested calibers with their specification are shown in Table 7. The movement evolutions for the different calibers are presented in Figure 5. Finally, the results are presented in Table 8. In addition, Table 8 also presents the penetration into the ceramic layer by each simulated caliber. This can be calculated by numeric integration of second stage movement equations, considering that the main penetration occurs in this stage.
In Figure 5 it is possible to observe the marked region. This region in the graph corresponds to the probable area of velocity and time where the erosion does not occur. In addition, the shifts in the curvatures of the velocities represent the beginning of the third stage. This helps to investigate the effect of the ceramic’s characteristics in the impact absorption and in the erosion. This graph and this model of simulation can be used to search for specific properties in the ceramic plate and/or to improve its shock absorption property.
The not highlighted region in Figure 5 represents a penetration stage where projectile erosion and ceramic scattering occur. If it is desired not to have a large deformation in the metallic plate, it is necessary to improve ceramic impact absorption. These improvements will affect , and , which will also modify the graph.
Considering the back metal plate as another important energy absorber and integrity maintainer for the system, a deterministic simulation was run using the new deformation law together with the deflection profile function. Then it gives
The variable is the radial distance of the projectile impact. The maximum deflection as a function of the hardening exponent and the radial distance is presented in Figure 6. In addition, the effect of the thickness of the plate together with the radial distance of the impact center is presented in Figure 7.
Based on Figures 6 and 7, it is possible to affirm that the mechanical property and the thickness of the metal layer have an important effect on the maximum deflection. The constant demonstrates the ductile behavior of the material. However, if a lower deflection is required, it is important to note that a harder material could shear more easily depending on the energy of the projectile. Moreover, a change in the thickness of the metal plate will certainly affect the structure weight, which could be an aggravating factor in the carrier’s mobility.
The values obtained from the model demonstrate how important the ceramic plate is in the protection, responsible for absorbing approximately 85% of the total energy. Moreover, the value for the deflection of the metal plate was close to the experimental data, which validates the model. Together with the optimized deflection equation the model could also prevent the effects of different projectile impacts.
In addition, the simulation demonstrated different aspects of the model and some predictable effects in collision phenomena. The control of the properties and its effects were analyzed with the developed computational method. Some of the effects predicted by the program could not be studied experimentally. However, future studies can use the presented results to validate the theory and, then, analyze some internal phenomena in a deeper way.
The present work is a theoretical analysis of the impact phenomenon in a specified type of protection. In the future studies, the current model and its solutions and simulations can be used to perform experimental test to evaluate the reliability of the model. Also, if the model can be considered valid for the impact, the considerations and formulations can be kept or improved to a more advanced modeling as finite element modeling (FEM). In addition, the experimental confirmation of the presented modeling also permits deeper investigation of the impact phenomenon, such as the effect of the interfacial friction between the projectile and sheared surfaces.
The authors wish to thank CAPES and CNPq for supporting this work.
D. Tabor, The Hardness of Metals, Oxford University Press, 1951.
A. Tate, “A theory for the deceleration of long rods after impact,” Journal of the Mechanics and Physics of Solids, vol. 15, no. 6, pp. 387–399, 1967.View at: Google Scholar
R. Hill, The Mathematical Theory of Plasticity, Oxford University Press, 1950.
H. A. Al-Qureshi and J. D. Bressan, “Investigation of the degree of biaxiality on the limit strains in sheet metal stretching,” in Proceedings of the 9th North American Manufacturing Research Conference (NAMRC '81), pp. 538–541, 1981.View at: Google Scholar
D. Ishikura and H. A. Al-Qureshi, “An investigation of perforation of metallic and composite plates by projectiles,” in Proceedings of the 5th Pan-American Congress of Applied Mechanics (PACAM '97), pp. 194–197, 1997.View at: Google Scholar
M. L. Wilkins, “Mechanics of penetration and perforation,” International Journal of Engineering Science, vol. 16, no. 11, pp. 793–807, 1978.View at: Google Scholar
D. Halliday, R. Resnick, and J. Walker, Fundamentals of Physics, vol. 2, John Wiley & Sons, New York, NY, USA, 4th edition, 1996.
U.S. National Institute of Justice, “Ballistic resistance of body armor,” NIJ Standard–0101.06, 2008.View at: Google Scholar