Abstract
Systematic three-dimensional finite element (FE) simulations are carried out to study the ballistic protection performance of double-layer sandwich plates having metallic pyramidal lattice truss cores filled with ceramic prism insertions and void-filling epoxy resin. Both normal and oblique projectile impacts are considered in the FE simulations that are validated against experimental measurements. The ballistic limit velocity, the energy absorbed by key constituting elements and the critical oblique angle corresponding to the transition from ballistic perforation to projectile embedment are calculated. As the oblique angle is increased, the evolution of deformation and failure in the double-layer plates as well as the underlying mechanisms are explored. It is demonstrated that the proposed double-layer sandwich plates outperform both the single-layer sandwich plates and the homogeneous (monolithic) metallic plates having equal total mass, and the top layer (the ceramic insertions in particular) of the double-layer configuration plays a more dominant role in energy absorption.
Keywords
Introduction
The concept of layered plates has been explored in order to enhance the ballistic perforation resistance over a monolithic plate (often made of high-strength steel). 1 The ballistic resistance of layered steel plates (both with and without spacing) struck by 7.62 mm standard bullets was investigated experimentally by Almohandes et al. 2 A single steel plate was found to be more effective than a layered plate of an equal total thickness and the perforation resistance of layered plates increased as the number of plates decreased and the thickness of the back plate increased with the same total thickness. Using a combined experimental and numerical approach, Gupta et al. investigated the behaviour of thin layered aluminum plates of different thicknesses under the impact of flat, ogival and hemispherical steel projectiles. For the double-layered plates, the residual velocity of the projectile was found to be comparable to that of a single plate of equivalent thickness. 3 When the number of layers was increased, however, the single plate offered more resistance against perforation which is consistent with the result of Almohandes et al. 2 However, Teng et al. studied numerically the protection performance of double-layered metal shields against blunt and conical-nose projectiles impact and found the double-layered targets exhibited advantage over monolithic plates of equal weight. 4
In addition to layered plates, to increase the ballistic resistance of monolithic plates, lightweight sandwich constructions having highly porous two-dimensional) or three-dimensional (3D) periodic metallic lattice truss cores were investigated. 5 –9 For example, ceramic prisms were inserted into the interstices of pyramidal lattice trusses to construct pyramidal-ceramic hybrid-cored sandwich plates to increase the performance of empty pyramidal lattice truss structures under projectile impact. 10,11 It was demonstrated both experimentally and numerically that the sandwich plates with pyramidal-ceramic-epoxy hybrid core as illustrated schematically in Figure 1(a) possessed superior ballistic resistance and energy absorption capability over those without ceramic inserted or epoxy infiltrated as well as a monolithic metal plate having an equal total mass. Whilst the projectile was severely eroded by the ceramic insertions, the infiltrated epoxy resin adhered all the sub-structures as an integrated whole, thus improving significantly the energy absorption capacity of the hybrid-cored sandwich structure. 11 It was also shown that, with a certain total mass of the sandwich structure, decreasing the thickness of its front face-sheet whilst increasing that of the back face-sheet could enhance further its ballistic protection capability. 11 Bazle Z (Gama) et al. also presented a combined theoretical-semiempirical penetration model of ballistic penetration of thick section composites. 12 Shanazari H et al. developed an analytical model based on wave propagation and energy balance between the projectile and target to analyze hybrid fabric panels for ballistic protection. 13 Tasdemirci A et al. discussed the effect of interlayer on the damage formation in a ceramic/composite armour. 14,15 The effect of interlayer which can reduce stress wave transmission was demonstrated experimentally and numerically.

Schematic of (a) single-layer and (b) double-layer sandwich plate having metallic pyramidal lattice truss core filled with ceramic prisms and epoxy resin. Hemispherical projectile (bullet) impacting plate center at (c) normal angle and (d) oblique angle.
Existing studies of layered metal plates or hybrid-cored sandwich constructions focussed mainly upon the most critical impact condition, namely, the normal impact case where the projectile impacts the target perpendicularly. However, the oblique angle of a projectile affected not only the penetration depth and energy absorbed but also its deformation/failure modes of the target. 16 –21 In general, the penetration resistance of the target increases with increasing angle of obliquity. In the present study, double-layered sandwich plates having pyramidal-ceramic-epoxy hybrid cores as shown in Figure 1(b) are constructed and their ballistic performance under normal projectile impact is examined with numerical simulations. The aim is to investigate whether the layered sandwiches can enhance further the ballistic performance of a single-layered sandwich. The effect of oblique projectile impact on the ballistic resistance of these double-layered sandwiches is quantified and the critical oblique angle at which the penetration process changes from perforation to embedment is identified. The accuracy of the numerical simulations is verified against experimental measurements under normal ballistic impact.
Problem description
With reference to Figure 1(a) and (b), consider both single- and double-layered sandwich plates having metallic pyramidal lattice truss cores filled with ceramic prisms and epoxy resin. The hybrid-cored sandwich plates are clamped and shot by a hemispherical projectile at the plate center. Both normal and oblique impact conditions are considered, as shown schematically in Figure 1(c) and (d), respectively.
For the single-layered plates,
Three-dimensional finite element model
The commercially available finite element package ANSYS is used to generate 3D meshes for both single- and double-layered sandwich plates, which are then transferred to the explicit integration version of the FE code LS-DYNA 971 R5.0 to carry out numerical simulations. 22
To ensure numerical convergence, each double-layered plate is meshed with 10,38,637 elements, whilst the projectile is meshed with 5632 elements. The minimum element size is 0.56 mm for the projectile, 0.15 mm for pyramidal lattice trusses, 0.2 mm for the central ceramic prism directly under impact, 0.8 mm for the face-sheets, 0.48 mm for air, and 0.48 mm for the epoxy bonding layer, respectively. The region directly under projectile impact is modelled with relatively dense meshes to highlight the penetration details. To emulate the penetration process, ‘surface to surface eroding contact’ is employed to model the contacts between the sub-structures of the sandwich and the projectile. Contacts are updated once the elements on free surfaces of the structure are deleted according to relevant material failure criteria. A scale factor of 0.5 for sliding interface penalties is selected to define all the contacts. 11
In the current study, the face-sheets, the lattice truss core and the projectile are made of AISI 304 stainless steel and modelled by the Johnson and Cook model with the consideration of strain rate effect, which has the form of
where,
where, Δεp is the increment of equivalent plastic strain and
For the ceramic (AD 98 alumina) prism insertions, the Johnson-Holmquist-Ceramics (JH-2) constitutive relation and fracture criterion are adopted, 22 which has been validated for describing the fracture behaviour of ceramics under high velocity penetration. 11 The model assumes that the strength of a ceramic is closely related with the exerting pressure, strain rate and damage. The damage is defined as the ratio of cumulative strain to failure strain and the relationship between the pressure and the specific heat capacity is included in the volume effect of the material. The JH-2 constitutive relation has the following form,
where,
where, Δεp represents an integration cycle of the material plastic strain and
Due to the low strength of the void-filling epoxy compared with either steel or ceramic, it is regarded here as a kind of hydrodynamic material and modelled by applying the Elastic-Plastic-Hydro constitutive relation. 23
The epoxy resin used for bonding has a mass density of 1185 kg/m3, shear modulus of 769 MPa and yield stress of 70 MPa. More details of the present FE model as well as material constitutive models and fracture criteria for AISI 304 stainless steel and AD 98 alumina can be found in Ni et al., 11 and hence omitted here for brevity. The tensile modulus of the steel and ceramic mentioned here is 200 GPa and 386 GPa separately.
For single-layered plates, the FE model as described above has been successfully validated against experimental measurements.
11
In the present study, the FE model is firstly employed to assess how the double-layered sandwich construction favours over its single-layer counterpart in terms of ballistic protection and energy absorption capability under normal projectile impact, with further validation of the FE model performed by comparing the simulation results with experimental measurements of double-layered plates. Subsequently, as the oblique angle
Response of double-layered sandwich under normal ballistic impact
In this section, subjected to the constraint of equal total mass, the ballistic responses of three different kinds of double-layered sandwich plates are numerically simulated and compared with their single-layer counterpart that has been investigated experimentally and numerically by Ni et al.
11
Table 1 listed the key geometrical parameters associated with both the single- and double-layered sandwiches. Depending upon the value of core height ratio of
Specifications and simulation results of single- and double-layered sandwich plates having pyramidal-ceramic-epoxy hybrid cores.
Experimental validation
Figure 2 illustrates schematically the experimental setup, with a ballistic rifle used to launch the impact projectile. A pair of brake screens spaced 300 mm apart and connected to a timing device is used to measure the projectile entry velocity. This is achieved by evaluating the ratio of screen spacing divided by the time difference between the moments separately triggered by the projectile impacting the first and second screens. In view of the drift effect in the process of projectile penetration through the sandwich, X-ray tomography is adopted to measure the exit velocity of the projectile. In a separate scheme, the exit velocity of the projectile is determined by evaluating the ratio of the distance between two projectile shadows in the X-ray native divided by the time difference between the two triggers.

Experimental setup for measuring ballistic response of hybrid-cored sandwich plate.
For single-layer sandwich plates, the FE model has been successfully validated against experimental results in Figure 3(a) by Ni et al.
11
Figure 3(a) also shows that the ceramic single plate possesses similar ballistic performance with pyramidal-ceramic sandwich in the same weight. Pyramidal-ceramic sandwich occupies better resistance to the projectile as the initial velocity increases. To verify the accuracy of the present FE model predictions for double-layer sandwich plates, a series of experimental measurements with sandwich 2 specimens (

Projectile exit velocity plotted as a function of initial impact velocity for single layer sandwich (a) and double-layer sandwich 2 (b) (see Table 1 for definition): comparison between experimental measurements and FE simulations (FE simulation results for ceramic single plate which has the same weight with pyramidal-ceramic sandwich). FE: finite element.
For sandwich 2 specimens (see Figure 1 and Table 1) subjected to normal projectile impact, Figure 3(b) compares the experimentally measured projectile exit velocities with those calculated numerically for selected initial impact velocities. In general, the agreement between the test data and FE simulations is reasonable, thus validating further the FE approach employed in the present study. The difference may be attributed to the complexities involved with the ballistic tests, the wobbly flying altitude of the projectile during the test, and the complicacy associated with the 3D numerical simulations. Consequently, in the rest of this study, the FE model as described in Section 3 is employed to evaluate the ballistic performance of double-layer sandwiches and explore the underlying physical mechanisms. The numerically predicted ballistic limit velocities for the three double-layer specimens are summarized in Table 1; for comparison, the corresponding result for their single-layer counterpart having equal mass is also presented. Note that, at the ballistic limit velocity, the projectile just perforates the whole sandwich structure.
Deformation and failure modes
In terms of ballistic limit velocity, the results of Table 1 show that all the three double-layered sandwich plates are considerably superior to their single-layer counterpart, with sandwich 3 enjoying a 42% increase over that of the latter (1850 m/s vs. 1300 m/s). It is worth noting that, also subject to the constraint of equal mass, Ni et al. demonstrated that the single-layered sandwich plate outperformed its monolithic counterpart.
11
In other words, for sandwich plates having pyramidal-ceramic-epoxy hybrid cores, the proposed concept of layered sandwich construction is beneficial so far as ballistic protection is of concern. Further, for double-layered constructions, the ballistic velocity increased as the value of
For the four different sandwich plates listed in Table 1, Figure 4(a)–(d) presented separately the numerically simulated deformation and failure modes under projectile normal impact at the ballistic limit velocity. For comparison, the experimentally observed failure configurations of sandwich 2 specimen after projectile penetration at the initial impact velocity of 1500 m/s are presented in Figure 4(e) and (f).

Evolution of deformation and failure under normal projectile impact at ballistic limit velocity: (a) single-layer sandwich (1300 m/s); (b) double-layer sandwich 1 (1630 m/s); (c) double-layer sandwich 2 (1780 m/s); and (d) double-layer sandwich 3 (1850 m/s). Experimentally observed failure modes of double-layer sandwich 2 under normal projectile impact at 1500 m/s. (e) Side view; (f) surface morphology of the bottom sandwich layer upon removing the bottom face-sheet. Euler elements for epoxy resin were marked for clarity.
As shown in Figure 4(a), when the projectile is penetrating across the single-layer plate, the ceramic prisms directly taking the impact slows down the projectile whilst the stress wave spreads out from the central impact zone to the whole structure. The spreading is attributed mainly to the epoxy bonding effects between the truss members, the face-sheets and the ceramic prisms. Correspondingly, as a result of projectile penetration, apart from the central ceramic prism, the side ceramic prisms have also fractured. The small ceramic fragments are constrained by the epoxy resin, enabling further resistance to projectile penetration. Further, due to the epoxy bonding effects, the lattice truss members in the central impact zone deflects much less in comparison with those in single-layer sandwiches without void-filling epoxy, and the back face-sheet plays a more leading role in energy absorption than the front one (results not shown here for brevity). Towards the end of the penetration (∼90 μs), the projectile is severely eroded by the central ceramic insertion. These results are consistent with those reported by Ni et al. 11 Compared to single AISI 304 stainless steel plate and ceramic-cored sandwich plate having the same total mass, Ni et al. also demonstrated that the single-layer sandwich with pyramidal-ceramic-epoxy core prolonged the penetration period of the projectile, so that each sub-structure of the hybrid-cored sandwich can deform and absorb the impact energy more sufficiently. 11
The numerically predicted deformation and failure modes of the three double-layered sandwich plates are presented in Figure 4(b)–(d). Relative to their single-layer counterpart of Figure 4(a), the increased complexity of the double-layer configuration leads to increased interfacial complicacy, enabling more efficient energy absorption and stress wave deflecting/spreading. Upon perforating the top face-sheet at high speed, the ballistic projectile rapidly enters the top sandwich layer. The ceramic insertions and the supporting metal lattice trusses absorb the impact energy by fragmentation and plastic deformation whilst the bottom sandwich layer acts as a load support. Simultaneously, the propagation speed of the projectile decelerates, accompanied by large deformation, abrasive erosion and mass loss. With further penetration, the middle face-sheet and the bottom sandwich layer start to absorb the remaining impact energy by plastic deformation and ceramic fragmentation. Due to the epoxy bonding effects, apart from the central ceramic insertions, the neighbouring ceramic insertions in both the top and bottom sandwich layers are also fragmentized. Further, subject to the constraint of equal mass, the results of Figure 4(b), (c) and (d) suggest that, as the value of
With reference to Figure 4(b), the projectile at a speed of 1630 m/s starts to penetrate the top layer of sandwich 1 at about 2 μs. Its impact energy is absorbed by a multitude of mechanisms, including localized plastic deformation and perforation fracture of the top face-sheet, plastic deformation, bending flexure, fracture and node failure of the lattice truss members, and fragmentation fracture of the central ceramic insertion. At the same time, progressive erosion of the projectile occurs due to the presence of the ceramic insertion and its fragments (constrained by the surrounding epoxy bonding layer). At about 8 μs, the middle face-sheet started to absorb the impact energy by plastic deformation even though the projectile tip is yet to reach it. At about 12 μs, the projectile reaches the middle face-sheet, causing further deformation of the latter and initiating deformation/fracture of the lattice trusses and ceramic insertions in the bottom layer. At about 20 μs, upon fully perforating the middle face-sheet, the projectile starts to penetrate the remaining sandwich structure. At about 28 μs, the projectile penetrates into the bottom layer completely, interacting now directly with the central ceramic insertion and the supporting lattice truss members. Due to the continuing bending deformation of the middle face-sheet and the lateral spreading of stress waves, fragmentation fracture occurs in the side ceramic insertions whilst the bottom face-sheet start to deform plastically and deflect outwards. The projectile completes the penetration process across sandwich 1 at about 80 μs, severely eroded, with only a small fraction of its initial mass left.
The predicted penetration process of the projectile at ballistic speed across sandwich 2 and 3 (see Figure 4(c) and (d)) is similar to that of sandwich 1 as described above, which agrees reasonably well with the experimentally observed failure configurations of sandwich 2 shown in Figure 4(e) and (f). One notable exception, however, is that, at about 38 μs, even before the projectile fully perforates the middle face-sheet of sandwich 3, severe ceramic fragmentation occurs in the bottom layer although the projectile is yet to interact directly with the ceramic insertions in the bottom layer. The results of Figure 4(a)–(c) also suggest that, as the value of
Energy absorption by sub-structures
For both single- and double-layer sandwich plates, the predicted energy absorption percentage is plotted in Figure 5 as a function of the initial projectile impact velocity. For clarity, only results for the main constituting elements (i.e. ceramic insertions and face-sheets) are presented, as the amount of energy absorbed by the lattice trusses and void-filling epoxy is relatively small. For the single-layer plate, the results shown in Figure 5(a) indicates that when the initial impact velocity is relatively small, the projectile interacts mainly with the upper portion of the sandwich: the amount of energy absorbed by the top face-sheet via plastic deformation and perforation fracture is the highest, followed by the ceramic insertion directly under impact. As the initial velocity increases, whilst the amount of energy absorbed by the top face-sheet monotonically decreases, the percentage of energy absorbed by both the ceramic insertions and the bottom face-sheet gradually increases, each surpassing that absorbed by the top face-sheet at the ballistic velocity of ∼1300 m/s; see Figure 5(a).

Energy absorption percentage of each sub-structure plotted as a function of initial project impact velocity for: (a) single-layer sandwich, (b) double-layer sandwich 1, (c) double-layer sandwich 2 and (d) double-layer sandwich 3.
As shown in Figure 5(b)–(d), when the initial velocity of the projectile impacting a double-layer plate is small (approximately 400 m/s), the energy absorption percentage of the ceramic insertions in the upper layer is comparable with that of the top face-sheet. However, as the initial velocity is increased, the ceramic insertions in the upper layer play a more leading role in energy absorption. This is attributed to the fact that, when the total mass of the double-layer plate is kept the same as that of its single-layer counterpart, the thickness of its three face-sheets decreases from 1.5 mm to 0.5 mm (Table 1), reducing therefore the amount of energy absorbed by the front face-sheet. As the initial velocity further increases, the percentage of energy absorbed by ceramic insertions in both the top and bottom layers increases, whilst that absorbed by the front face-sheet decreases.
For each double-layer plate considered, the energy absorbed by ceramic insertions in the top layer reached a maximum when the initial velocity is increased to about 1150 m/s (i.e. before the ballistic limit is reached). Below this velocity, the projectile interacts mainly with the front face-sheet and the constituting elements of the upper layer; when it penetrates into the bottom layer, its velocity is considerably decreased so that the energy absorption capability of the bottom layer is not efficiently utilized. However, as the initial velocity is increased beyond 1150 m/s, the residual velocity of the projectile entering the bottom layer is sufficiently high so that its main constituents can absorb enough amount of the remaining impact energy. As a result, the energy absorption percentage of both the ceramic insertions and the bottom face-sheet in the bottom layer increases, whilst that of the front face-sheet and top-layer ceramic insertions decreases; see Figure 5(b)–(d).
For all the three double-layer plates considered, the amount of energy absorbed by the top-layer ceramic insertions is always the highest at or below the ballistic velocity. When the bottom layer has a thickness twice that of the top one (i.e. sandwich 1), the maximum energy absorption percentage of the top-layer ceramic insertions is only about 45%. In contrast, for sandwich 3 having a top-layer thickness twice that of the bottom one, the maximum energy absorption percentage of the top-layer ceramic insertions is about 67%. As a result, the ballistic velocity of sandwich 3 is the highest, about 13% higher than that of sandwich 1 and 4% higher than that of sandwich 2.
It should be mentioned that the main focus of the present investigation is placed upon evaluating the concept of double-layer sandwich constructions with pyramidal-ceramic-epoxy hybrid cores for ballistic protection applications. Consequently, compared with single-layer plates having equal mass, the thickness of the face-sheets selected, 0.5 mm, is relatively small. As a result, as shown in Figure 5(b)–(d), the contribution of the face-sheets to total energy absorption is small, particularly at relatively high initial velocities. Also note that, as the initial velocity is increased, the percentage of impact energy absorbed by both the middle and bottom face-sheets increases, each surpassing that of the top face-sheet at the ballistic velocity.
Response of double-layered sandwich under oblique ballistic impact
Consider next the ballistic responses of double-layer sandwich plates subjected to oblique projectile impact. To search for the critical oblique angle at which the projectile penetration process is changed from perforation to embedment, four oblique angles that is, 15°, 30°, 45°and 60° are selected, as shown schematically in Figure 1(d). For each oblique angle selected, the ballistic trajectory and erosion of the projectile, the deformation and failure modes of the sandwich structure, and the energy absorption percentage of each constituting element are calculated.
Critical oblique angle
Figure 6 presents the predicted residual velocity of the impact projectile for selected oblique angles, with its initial velocity fixed at 2500 m/s. It is seen that the projectile with oblique angle equal to 0° (i.e. normal penetration with

Numerically predicted residual velocity of impact projectile penetrating across double-layer sandwich plate at initial velocity of 2500 m/s for selected oblique angles.

Evolution of deformation and failure in sandwich 1 under projectile impact at initial velocity of 2500 m/s for selected angles of obliquity: (a)
Deformation and failure modes
As a complement of Figure 6, the evolution of deformation and failure in sandwich 1 subjected to projectile impact at the initial velocity of 2500 m/s is presented in Figure 7 for varying angles of obliquity.
When the oblique angle is small (
For the case of
When
As shown in Figure 7(d), with its oblique angle equal to
The penetration process of the projectile having varying oblique angles across either sandwich 2 or sandwich 3 is similar to that described above for sandwich 1 and hence not be repeated here. However, it is established that, rather than
Residual mass and transverse movement trajectory of impact projectile
To explore further the mechanisms underlying the penetration process of the projectile having different oblique angles, Figure 8(a), (c) and (e) plot

Effect of oblique angle on ratio of projectile residual mass to initial mass (
For all the three double-layer sandwich plates considered, it is seen from Figure 8 that, for a given
The results of projectile residual mass and trajectory movement shown in Figure 8(a) and (b) for sandwich 1 are noticeably different from those of Figure 8(c)–(f) for sandwich 2 and 3. This is attributable to the fact that the former has a critical oblique angle of
Energy absorption by constituting elements
To compare the efficiency of impact energy absorption among the double-layer sandwich plates as listed of Table 1, Figure 9 plots the energy absorption percentage of key constituting elements as a function of oblique angle (initial impact velocity fixed at 2500 m/s). The results of Figure 9 indicate that ceramic prism insertions play a more dominant role in absorbing the impact energy than the three face-sheets. For the configuration of sandwich 1 (

Energy absorption percentage of key constituting elements plotted as a function of oblique angle with projectile initial impact velocity fixed at 2500 m/s for double-layer plates: (a) sandwich 1; (b) sandwich 2; and (c) sandwich 3.
As to sandwich 2 (
The energy absorption performance of sandwich 3 (
Concluding remarks
The concept of double-layer sandwich construction with metallic pyramidal lattice-ceramic insertions-epoxy filling hybrid core for ballistic protection is evaluated using 3D finite element simulations. Both normal and oblique projectile impact conditions are considered. The numerical predictions are successfully validated by comparing with results measured experimentally using ballistic rifle. It is demonstrated that, under the present simulation conditions, the ballistic protection performance of the proposed double-layer sandwich plate is superior to its single-layer counterpart of equal mass, with the former achieving a ballistic limit velocity of about 1850 m/s under normal projectile impact whilst that of the latter is only about 1300 m/s. The presence of the ceramic insertions causes serious erosion and mass loss of the projectile, whilst the void-filling epoxy resin contributes to further enhance the ballistic resistance of the hybrid-cored sandwich by adhering its constituting elements as an integrated whole.
Subject to the constraint of equal mass, increasing the height (and hence mass) of the top layer facing the initial strike of the impact projectile enhances the ballistic protection capability of the double-layer sandwich. Increasing the angle of obliquity leads to increased transverse movement of the projectile, which in turn prolongs its interaction with the constituting elements of the hybrid core. As a result, the erosion effects of the ceramic insertions on the penetrating projectile intensifies, causing more serious loss of its mass. In comparison with homogeneous (monolithic) metal plates, ricochet rarely occurs in the present double-layer sandwich plates as their relatively thin top face-sheet enables the projectile to perforate with ease. Also, embedment is more likely to occur in the hybrid-cored sandwiches in comparison with monolithic plates, especially at large oblique angles.
The results presented here would be helpful for designing lightweight layered constructions having superior ballistic protection capabilities. Future research is needed to investigate a number of issues not addressed in the present study, including the influence of face-sheet thickness, the position of projectile impact, and the penetration resistance of the hybrid-cored sandwich to multiple projectile impacts.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was financially supported partially by the National Basic Research Program of China (2011CB6103005), partially by the National 111 Project of China (B06024), partially by National Natural Science Foundation of China (11472208) and partially by the Shaanxi Province 13115 Project. The authors gratefully acknowledge the support by the opening project (project number KFJJ13-11 M) of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology) and the open fund SV2014-KF-17 of State Key Laboratory for Strength and Vibration of Mechanical Structures.
