A Numerical Study of the High-Velocity Impact Response of a Composite Laminate Using LS-DYNA
A Numerical Study of the High-Velocity Impact Response of a Composite Laminate Using LS-DYNA
International Journal of Aeronautical and Space Sciences. 2010. Sep, 11(3): 221-226
Copyright ©2010, The Korean Society for Aeronautical Space Science
  • Published : September 15, 2010
Export by style
Cited by
About the Authors
Jeoung-Hee Ahn
Khanh-Hung Nguyen
Yong-Bin Park
Jin-Hwe Kweon
Jin-Ho Choi

The failure of a Kevlar29/Phenolic composite plate under high-velocity impact from an fragment simulation projectile was investigated using the nonlinear explicit finite element code, LS-DYNA. The composite laminate and the impactor were idealized by solid elements, and the interface between the laminas was modeled as a tiebreak type in LS-DYNA. The interaction between the impactor and laminate was simulated using a surface-to-surface eroding contact algorithm. When the stress level meets the given failure criteria, the layer in the element is eroded. Numerical results were verified through existing test results and showed good agreement.
1. Introduction
As the application of composite materials to aircraft structures has increased in various fields (Niu, 2005; Park et al. 2009), ballistic-tolerant design concepts have become a topic of interest for designers. Particularly in the design of fighters, it is necessary to consider the ballistic resistance of structures that are exposed to the potential danger of ballistic impact. In helicopters, survivability after a ballistic impact is required for the rotor blades, control, and the hub.The application of composite materials to armor already has a long history due to not only the ballistic resistance but also the excellent specific stiffness and strength offered by these materials.
The behavior of composite structures during a highvelocity impact is affected by a variety of factors, including the impact velocity, shape, mass and stiffness of the impactor, and laminate properties. The complicated failure mechanism of a composite laminate is an especially important obstacle in numerical simulations. Accordingly,numerical approaches regarding a high-velocity impact are considered to be a challenging issue. For this reason, mainly tests have been used in this area (Jung et al., 2006; Lee et al.,2001; Liu, 2004).
There are two approaches regarding numerical methods,i.e., microscopic and macroscopic. In the microscopic model, the constituents of composite materials such as the reinforcing fiber and matrix are treated as separate materials(Ji and Kim, 2006). In another microscopic method,composite materials are idealized as a combination of unit cells consisting of the fiber and the matrix (Tabiei and Ivanov, 2002; Unosson and Buzaud, 2000). The macroscopic model assumes a composite laminate to be a structure made of a homogenous and anisotropic material. It is useful and efficient for simulating the structural behavior, as the constituents of the composite materials are not separately handled. However, it is somewhat limited in simulations of failure due to the combined effect of the fiber and the matrix(Ji and Kim, 2006).
While a number of papers (Kong et al., 2007; Lee et al.,2008; Park et al., 2008) are available related to low-velocity impacts on composite structures, a limited number of studies related to high-velocity impacts are available to the public.Even in the available literature, only limited information is provided, and data such as the stacking sequence and material properties are frequently excluded.
Tabiei et al. (2001, 2002) performed an analysis of a highvelocity impact on flexible Kevlar fabric using a homogenized microscopic unit cell approach. They then compared their results with test results. LS-DYNA provides the user defined material (UMAT) function to simulate a failure through impact on composite structures. Using this method, Hoof(1999) and Yen (2002) conducted numerical analyses of laminated composite structures.
The present paper has used a numerical analytic method for a Kevla29/Phenolic fabric laminate under the highvelocity impact of an fragment simulation projectile (FSP). A nonlinear explicit finite element program, LS-DYNA version 971 (Livermore Software Technology Corporation, 2007),was applied using the element eroding contact algorithm.In the algorithm, where the stress level exceeds a certain criterion, the associated element is eroded from the finite element model. Delamination is an important failure mode in impact problems. The tiebreak contact interface was used to simulate delamination between composite layers. The finite element results were compared with the test results in an earlier study (Hoof, 1999).
2. Finite Element Analysis
- 2.1 Problem description
Lager Image
Composite laminate under a ballistic impact.
Lager Image
Configuration of the fragment simulation projectile (mm).
Finite element analysis was conducted for the Kevlar29/Phenolic fabric laminated composite plate shown in Fig. 1 .The plate has dimensions of 152.4 (length) × 101.6 (width)× 9.5 (thickness) (unit: mm). The long edges are free, and the short edges are simply supported. The impact velocities considered are 483, 545, and 586 m/s. The impactor FSP shown in Fig. 2 has a diameter of 5.59 mm and weight of 1.1 g. The FSP impacts upon the center of the plate. For finite element modeling, a symmetry condition is used. The material properties of the composite and the impactor are given in Table. 1 and 2, respectively (Hoof, 1999).
- 2.2 Finite element modeling
Key issues for finite element analysis of a composite laminate under a high-velocity impact include definitions of the material properties, development of a proper contact algorithm, and failure criteria. It is very difficult to define material properties that are valid for a wide range of impact velocities. Accordingly, it is important to find an adjustment factor for material properties to correlate the results of finite element analysis with test results. In the early design phase,many parameters for operating LS-DYNA were adjusted through iterative comparisons between the test and analytic results. To adjust the parameters, we use a trial and error method. First, we set the parameters as they are by default in the program and undertake analysis. We then change these parameters slightly when there is a mismatch between the results of numerical analysis and the test. After several iterations, suitable parameters can be found.
The finite element model is shown in Fig. 3 , where a quarter of the plate is idealized using the symmetry condition. An
Material properties of Kevlar29/Phenolic panels
Lager Image
Material properties of Kevlar29/Phenolic panels
Material properties of the impactor
Lager Image
Material properties of the impactor
eight-node hexagonal solid element was used to simulate the plate and the impactor. Each layer of the plate was modeled using two solid elements through the thickness, and the minimal element size was 0.25 mm × 0.83 mm at the center of the plate, which was the point of impact. The failure criterion for the composite material is a stress-based equation from Schweizerhof et al. (1998).
Lager Image
In Eq. (1),
σ 1 , σ 2 : Stresses in the fiber and transverse directions,respectively.
σ 12 , σ 13 , σ 23 : Shear stresses.
Y t, Yc : Tension and compression strengths, respectively, in the fiber direction.
S 12 , S 12 , S 23 : Shear strengths.
When the index, f , is equal to or greater than 0, the element is assumed to fail. It is interesting that in the equation, only five stress components, viz., σ 1 , σ 2 , σ 12 , σ 13 , and σ 23 , affect the failure criterion, with the exception of σ 3 . In the analysis, an additional constraint is imposed on the strain for the case in which further calculation is not possible due to a large amount of deformation. Although an element does not fail in terms of the stress criterion, it is eroded when the strain exceeds a critical value. An eroding surface-to-surface contact algorithm was adopted for the simulation of contact between the composite plate and the impactor (Hoof, 1999;Livermore Software Technology Corporation, 2007). The tiebreak contact interfaces are defined between adjacent layers of the composite plate to determine the delamination,as shown in Fig. 4 . When the combined value of the inter-
Lager Image
Finite element model of a composite plate and an impactor.
Lager Image
Definition of the interfacial contact.
laminar normaland shear stresses exceeds a critical value,the tiebreak contact interface is deleted and delamination occurs (Hoof, 1999; Ji and Kim, 2006; Livermore Software Technology Corporation, 2007). The failure of the tiebreak contact interface is determined by a combination of the fracture modes, I, II, and III, as shown in Fig. 5 . The failure criterion for the tiebreak is given in Eq. (2)
Lager Image
In Eq. (2)
Sn : Inter-laminar tension strength.
SS : Inter-laminar shear strength.
The steel impactor was modeled with an elastic-plastic material. The von-Mises yield criterion for the impactor is defined as follows (Livermore Software Technology Corporation, 2007):
Lager Image
Lager Image
with the following notation.
J2 : Second order stress invariant
sij : Deviatoric stresses.
The yield stress, σ y , is defined in terms of the effective
Lager Image
Inter-laminar fracture modes.
Lager Image
Perforating behavior of a composite panel (t = 200 μs).
Lager Image
Temporal variation of the fragment simulation projectile velocity(V = 586 m/s).
Lager Image
Temporal variation of the backplane deflection (V = 483 m/s).
plastic strain, ε p eff' , and the plastic hardening stiffness, E p .
Lager Image
In Eq. (5),
σ 0 : Specified threshold stress.
- 3. Results and Discussion
The failure, perforation, backplane deflection, and velocity of the impactor are the important parameters in the analytic results for a high-velocity impact on a composite plate.The behavior of the composite panel 200 μs after the initial contact from a 586 m/s-velocity impactor is shown in Fig.6 . As shown in the figure, the panel was nearly completely perforated by the impactor. Extensive delamination occurred between the bottom plies of the panel, whereas this was not observed in the upper plies. Figure 7 shows the velocity-time history of the impactor. The present results of finite-element
Lager Image
Temporal variation of the backplane deflection (V = 545 m/s).
Lager Image
Temporal variation of the backplane deflection (V = 586 m/s).
analysis are in good agreement with the aforementioned test results (Hoof, 1999).
The results regarding the backplane deflection of a composite plate under various impact velocities are shown in Figs. 8 - 10 in comparison with the test and numerical results based on the UMAT model (Hoof, 1999). As discussed in a relevant study, namely (Hoof, 1999), the initial time of commencement of deflection was not precisely measured during the test. Except for the very early stage of the impact event, the deflection was believed to have been measured correctly. Accordingly, the authors of the aforementioned study, namely (Hoof, 1999), did not shift the test curve to position the starting point at the origin. This is also true for the present paper. The finite element results here show excellent agreement with the test results, including the maximum deflection, apart from the early stage of the event.The springback phenomenon is natural after the maximum deflection under the impact occurs. While the test results do
Lager Image
Failure of a composite plate (t = 25 μs).
Lager Image
Failure of a composite plate (t = 50 μs).
not account for this phenomenon, it is easily observable in the finite element analyses.
Figures 11 - 14 show the progress of damage in the composite plate at an impact velocity of 586 m/s. The entire process from the initial impact to the final penetration of the impactor through the plate was successfully simulated. Figure 11 shows the damage at 25 μs after the initial strike by the impactor. Regardless of the short time for the response,the impactor penetrated into more than half the plate and the backplane deflection reached nearly 3 mm. However,the damage zone was limited to the region close to the impact point. The finite elements where the impactor passed through failed and were therefore eroded. Figure 12 shows the results at 50 μs after the initial impact. The maximum backplane deflection was approximately 6 mm, which was double the value at t = 25 μs. Extensive delamination around the impactor and through the thickness are observed in the figure. Figures 13 and 14 show the damage of a composite plate at t = 150 and 200 μs, respectively. Delamination in the layers close to the backplane occurred over nearly the entire plate, and the impactor began to leave the plate after penetration.
4. Conclusions
A finite element model using LS-DYNA was developed to simulate the high-velocity impact response of a Kevlar29/Phenolic composite plate. The interface between the laminas was modeled by a tiebreak contact algorithm, while the interaction between the impactor and the laminate was simulated using a surface-to-surface eroding contact algorithm. Numerical analyses were conducted at impact velocities of 483, 545, and 586 m/s of a steel impactor. The
Lager Image
Failure of a composite plate (t = 150 μs).
Lager Image
Failure of a composite plate (t = 200 μs).
results of finite element analysis were in good agreement with the test results and in especially excellent agreement regarding the maximum backplane deflection. The present finite element model also successfully simulated the progress in damage from the initial impact to the final penetration of the composite plate.
This work was supported by Priority Research Centers Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Education, Science and Technology (2009-0094015) and the Degree and Research Center for Aerospace Green Technology (DRC) of the Korea Aerospace Research Institute (KARI) funded by the Korea Research Council of Fundamental Science & Technology(KRCF).
Hoof J. V , PhD Thesis 1999 Modeling of Impact Induced Delaminationin Composite Materials. Carleton University PhD Thesis
Ji K. H , Kim S. J 2006 Direct numerical simulationof composite laminates under low velocity impact Journal ofthe Korean Society for Composite Materials 19 1 - 8
Jung W. K , Ahn S. H , Lee W. I , Kim H. J , Kwon J. W 2006 Fracture mechanism of ceramic/glass-fiberreinforced-composite laminate by high velocity impact. Journal of the Korean Society of Precision Engineering 23 170 - 176
Kong C. D , Kim Y. K , Lee S. H 2007 Comparative study on low-velocity impact behavior of graphite/epoxy composite laminate and steel plate Journal of the Korean Society for Composite Materials 20 1 - 6
Lee B. W , Lee S. H , Kim S. G , Yoon B. I , Paik J.G 2008 A study on the low speed impact response andfrictional characteristic of shear thickening fluid impregnatedKevlar fabrics Journal of the Korean Society for CompositeMaterials 21 15 - 24
Lee B. L , Walsh T. F , Won S. T , Patts H. M , Song J. W , Mayer A. H 2001 Penetration failure mechanismsof armor-grade fiber composites under impact
Liu D 2004 Characterization of impact properties anddamage process of glass/epoxy composite laminates. Journalof Composite Materials 38 1425 - 1442    DOI : 10.1177/0021998304042741
2007 Livermore Software Technology Corporation Livermore Software Technology Corporation Livermore CA LSDYNA Keyword User Manual
2006 Livermore Software Technology Corporation Livermore Software Technology Corporation Livermore CA LSDYNA Theory Manual
Niu M Composite airframe structures Hong Kong Conmilit Press Ltd Hong Kong
Park I. J , Jung S. N , Kim D. H , Yun C. Y 2009 General purpose cross-section analysis program forcomposite rotor blades. International Journal of Aeronauticaland Space Sicences 10 77 - 85    DOI : 10.5139/IJASS.2009.10.2.077
Park J. B , Kim D. R , Kim H. G , Hwang T. K 2008 Evaluation of residual strength of CFRP pressure vessel after low velocity impact Journal of the Korean Society for Composite Materials 21 9 - 17
Schweizerhof K , Weimar K , Munz T , Rottner T 1998 Crashworthiness analysis with enhanced composite material in ls-dyna-material and limits Proceedings of the 5th LS-DYNA World Conference Southfield MI.
Tabiei A , Chen Q 2001 Micromechanics basedcomposite material model for crash worthiness explicit finiteelement simulation Journal of Composite Materials 14 264 - 289
Tabiei A , Ivanov I 2002 Computational micromechanicalmodel of flexible woven fabric for finite elementimpact simulation International Journal for NumericalMethods in Engineering 53 1259 - 1276    DOI : 10.1002/nme.321
Unosson M , Buzaud E 2000 Scalar and Vectorized User Defined Material Routine in LS-DYNA. Methodology Report FOA Defence Research Establishment Tumba
Yen C. F 2002 Ballistic impact modeling of composite materials Proceedings of the 7th International LS-DYNA Users Conference Detroit MI.