Published: 31 December 2017

An analytical modeling for high-velocity impacts on woven Kevlar composite laminates

Chen Fangyu1
Zhou Li2
Tang Yihao3
Zhou Xiaoming4
1, 2State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing, China
3, 4China Helicopter Research and Development Institute, Jingdezhen, China
Corresponding Author:
Zhou Li
Views 168
Reads 83
Downloads 1773

Abstract

In this paper, an analytical model, which based on energy balance, is built to study the process of high velocity impacts on woven Kevlar composite laminates by a cylindrical projectile. Four different mechanisms, such as laminate crushing, linear momentum transfer and tensile fiber failure, and shear plugging, is absorbed by the laminate while impacting. Then, simplification of the model is done to obtain the residual velocity and ballistic limit. The analytical results are validated with the results of experiment, and the perturbation analysis is done to analyze the reason of error.

1. Introduction

Resistance to ballistic impacts is an important requirement in many industrial applications, such as car, space shuttle and aeronautical industry. Woven fabrics are well-known to be good choice to resist ballistic impacts [1-3]. In the past, several researches have been focused on the impact on composite laminates. Elasticity method [4-7] was widely used to predict the behavior of impact. Three phases, which namely indentation, perforation and exit, are divided in the impact event by Zhu et al. [4]. Sun et al. [5] applied different criteria to obtain damage initiation, progression and plug formation. Wen et al. [6, 7] applied projectiles with different shapes to investigate FRP laminates’ perforation and penetration. Another method to predict the behavior of impact was the energy method. The projectile’s residual velocity and ballistic limit are calculated by energy balance equation. Naik et al. [8] and Xuan Haijun et al. [9] study the perforation of the composite plates with all the damage mechanisms and verified the analytical results with experimental results. Impact test is done by Ganesh Babu et al. [10] using different shape projectiles on unidirectional glass/epoxy composite plates. Sikarwar et al. [11, 12] and Wang Yuanbo et al. [13, 14] studied the analytical and experimental results of the residual velocity and ballistic limit of glass/epoxy and Kevlar/epoxy composite laminates respectively, which are subjected to high velocity impact. Yahaya et al. [15] observed a linear increasing tendency of energy absorption against impact velocity for Kevlar-Kenaf hybrid woven composites. Gu Xingjin et al. [16] and Michele et al. [17] calculate energy absorbed in different mechanisms with flat-end cylindrical projectiles. Chen et al. [18] studied simulation and experimental results of damage and energy absorption behavior for woven composites under ballistic impact. The researches mentioned above spend a lot of attention to high precision of the solutions of different equations, so they are high cost. To make a quick assessment at early stages of the design process, Lopez-Puente et al. neglect elastic deformation contribution in impact event [19] and then proposes a low-cost model and obtains closed-form solution to predict the ball projectile’s residual velocity after impacting on thin CFRPs woven laminated plates [20], but the influence of shear plugging [21-23] is not considered in their research.

In this paper, an analytical model is proposed to build the energy balance equation and predict the residual velocity of a cylindrical projectile after impacting at woven Kevlar composites. Simplification of the model is done to obtain the residual velocity and ballistic limit. Also, three dimensionless variables are built to simplify the process of analysis perturbation, which caused by tensile fiber failure. Analysis is done to explain the error between the results of calculation and experiment [13]. The model is low-cost, and its predictive accuracy is improved. The method of simplification could be useful in preliminary assessment of engineering.

2. Analytical model

The analytical model is built by terms of energy balance. In Lopez-Puente et al.’ [17] opinion, the kinetic energy of the projectile Ek is absorbed by the laminate through three different processes. Ec is the energy of crushing, which leads to laminate breakage. Em is the energy of linear momentum, transferred from the projectile to the detached part of laminate. Ef is the energy of tensile fiber failure breakage. Elastic deformation energy is not considered by them in this analytical model. Naik et al. [19, 20]and Brenda et al. [21] consider that Esp, the energy absorbed by shear plugging, is also an important process. So, the basic balance equation can be presented in the following because of the influence of the shear plugging:

1
-dEk=dEc+Em+dEf+dEsp.

Following each energetic term of Eq. (1) is presented.

The kinetic energy of the projectile is expressed as:

2
dEkx=12mpdvx2,

where x is the distance between the position of the projectile and the impact surface of the composite laminate, mp is the mass of projectile, and v is the velocity varies with the position.

The energy of laminate crushing is expressed as:

3
dEcx=σcAxdx,

where Ax is the frontal projectile area that contacts the non-crushed laminate and σc is the out-of-plane compressive strength of the laminate. Here, Ax=πr2, where r denotes the radius of the projectile cross-section.

The energy of linear momentum transfer is expressed as:

4
dEmx=12Axdxρv2x,

where ρ is the density of laminate.

The energy of tensile fiber failure is expressed as:

5
dEfx=ωfdV,

where ωf is the specific energy and dV is the affected material volume at any dx. The specific energy is estimated as:

6
ωf=212Xtεf,

where Xt is the tensile strength in fiber direction and εf is the ultimate strain. Experimental result [9, 20] shows that the transverse section of laminate is a truncated pyramid, with a semi-angle α and β, and its diagonals aligned in the fiber directions. The upper base semi-diagonal is r and the lower base semi-diagonal is L=r+h-x0tanα with x0 the position of the projectile at the time x0, as Fig. 1 shown. The affected material volume dV is given by:

7
dV=2l×lwdx=2dxr+h-x0tanα2tanβ,

and the energy of tensile fiber failure is given by:

8
dEfx=2Xtεfr+h-x0tanα2tanβdx.

The value of x0 can be obtained from following equation:

9
0x0dxv(x)=2hE3/ρ,

where h is the thickness of laminate and E3 is the elastic modulus through-thickness.

Fig. 1Woven composite with cylindrical projectile

Woven composite with cylindrical projectile

Tensile fiber failure happens while x0<x<h, so a function Crx is defined as follow to multiplies the last equation:

10
Crx=0, 0<x<x0,1, x0<x<h,0, xh.

The energy absorbed of shear plugging is expressed as:

11
dEspx=SspSxxdx,

where Ssp is shear plugging strength and Sx is the cross-section perimeter of the projectile. Here, Sx=2πr.

While xx0, the crushing energy and shear plugging energy have no influence. So, a function Csx is defined as follow to multiply the Eq. (3) and Eq. (4):

12
Csx=1, 0<x<x0,0, xx0.

3. Model simplification

As analysis above, the basic balance equation can be written as:

13
-12mpdvx2=σcAxCsxdx+12Axdxρv2x
+2Xtεfr+h-x0tanα2tanβCrxdx+SspSxxCsxdx.

To solve the equation, some simplifications are made and the initial condition is imported:

14
-12mpd2vdx=12Axρv2+SspSxxCsx+σcAxCsx +2Xtεfr+h-x0tanα2tanβCrx,v0=v0.

Here the velocity and the displacement are converted to dimensionless variables:

15
v*=vv0, x*=xh.

By defining a variable y*=v*2, Eq. (14) leads to:

16
-dy*dx*=πr2hρmpy*+4Sspπrh2mpv02Cshx*x*+2πr2hσcmpv02Cshx* +4Xtεfr2hmpv021+hrx*-x0rtanα2tanβCrhx*,y*0=1.

To simplify the process of finding solution, the relative importance of every energy absorption mechanisms are estimated firstly. Then, inconsequential factors can be found and ignored.

The characteristic values of different variables involved are showed in Table 1.

Some variables are defined as a function of the impact velocity to estimate the relative importance of the terms that appear in the equation:

17
Rm=πr2hρmp-10-2,Rsp=4Sspπrh2mpv02-103v02,Rc=2πr2hσcmpv02-103v02,Rf=4Xtεfr2hmpv02-10v02.

Obviously, the value of Rf is very small. Within the velocity range considered, Rf is always lower than the others. The energy of tensile fiber failure is ignored to simplify the model as following:

18
-dy*dx*=πr2hρmpy*+4Sspπrh2mpv02Cshx*x*+2πr2hσcmpv02Cshx*,y*0=1.

The equation can be numerically solved to find the solution:

19
y*=-4Ssphv02rρx*+4Sspmpv02πr3ρ2-2σcv02ρ+e-πr2hρmpx*-4Sspmpv02πr3ρ2+2σcv02ρ+1, 0<x*<x0*,-4Ssphv02rρx0*+4Sspmpv02πr3ρ2-2σcv02ρ+e-πr2hρmpx0*-4Sspmpv02πr3ρ2+2σcv02ρ+1, x*x0*.

x0* can be obtained from following equation:

20
0x0*v0y*(x*)dx*=2hE3/ρ.

The expression for residual velocity is obtained, which is given by:

21
vr=limx*vx*=-4Ssphrρx0*+4Sspmpπr3ρ2-2σcρ+e-πr2hρmpx0*-4Sspmpπr3ρ2+2σcρ+v02.

For ballistic limit v50, vr= 0, so ballistic limit v50 can be expressed by:

22
v50=4Ssphrρx0*-4Sspmpπr3ρ2+2σcρ-e-πr2hρmpx0*-4Sspmpπr3ρ2+2σcρe-πr2hρmpx0*.

Table 1Characteristic values of different variables

Variable
Characteristic values
h [m]
10-3
r [m]
10-3
mp [kg]
10-3
ρ [kg/m3]
103
Xt [N/m2]
109
σc [N/m2]
108
Ssp [N/m2]
108
εf
10-2

4. Model validation and discussion

To validate the analytical model above, analytical results are calculated to compare with experimental result [13]. Here cylindrical projectiles with a diameter of 7.62 mm and a mass of 5.2 g are used to impact 0°/90° layup Kevlar laminates with thickness of 5 mm and 10 mm respectively. The mechanical properties of Kevlar composite are shown in Table 2.

The projectile residual velocities obtained from different laminate thickness values and projectile initial velocities are shown as Fig. 2. It is clear that the differences between results of calculation and experiment reduce as the initial velocity increases for laminates. For laminate with thickness of 5 mm, the error is about 13 % while the initial velocity below 400 m/s, and the error is less than 9 % while the initial velocity over 400 m/s. It shows that the calculated result is closer to experimental result for laminates with thickness of 10 mm than that of 5 mm. Also, good result was found in the residual velocity estimated by the mentioned method [19], but its error is larger than that found by the proposed model.

The ballistic limit obtained for different laminate thickness values are shown as Fig. 3.

It is clear that the difference between calculated result and experimental result is larger while the laminates with thickness of 5 mm than which of 10 mm. The ballistic limit of calculated result [19] is lower than that of proposed model.

Table 2Mechanical properties of Kevlar composites

Mechanical property
Value
ρ [kg/m3]
1440
Modulus in x direction Ex [GPa]
21
Modulus in y direction Ey [GPa]
21
Modulus in z direction Ez [GPa]
4.6
Shear modulus in xy direction Gxy [GPa]
5.2
Xt [MPa]
800
Ssp [MPa]
183
σc [MPa]
80
εf
0.02

Fig. 2Comparison of residual velocities obtained from analysis and experiment for different projectile initial velocities

Comparison of residual velocities obtained from analysis and experiment for  different projectile initial velocities

Fig. 3Comparison of ballistic limit obtained from analysis and experiment for different laminate thickness

Comparison of ballistic limit obtained from analysis and experiment for  different laminate thickness

The key of explaining the error is analyzing the perturbation caused by tensile fiber failure. It is also important to analyze the influence of the energy absorbed by other three different mechanisms. To simplify the process of analysis and estimate the influence of every mechanisms quickly, three dimensionless variables Rsp/Rf, Rm/Rf and Rc/Rf are built.

Fig. 4 shows that the perturbation caused by tensile fiber failure descends as soon as the initial velocity increases. Shear plugging is always the main absorbing energy mechanism (within the velocity range considered). At initial velocity below 370 m/s, the crushing is more important than the momentum transfer. At initial velocity below 100 m/s, the momentum transfer is almost without influence. As the impact velocity increases, the ratio of Rm/Rf rises quickly. The reason is that the ratio is proportional to square of the projectile velocity. It means that the influence of shear plugging descends as the initial velocity increases, so the solutions between the proposed model and method [19] become closer in Fig. 2.

Fig. 5 shows that the perturbation caused by tensile fiber failure descends as soon as the thickness of laminate increases. Shear plugging is the main absorbing energy mechanism while thickness over 1mm. Crushing is always more important than momentum transfer. As the impact velocity increases, the ratio of Rsp/Rf rises quickly. The reason is that the ratio is proportional to the thickness of laminate. It means that the influence of shear plugging rises as the initial velocity increases, so the solutions between the proposed model and method [19] become larger in Fig. 3.

Fig. 4Ratio of different factor for different projectile initial velocities (thickness = 5 mm)

Ratio of different factor for different projectile initial velocities (thickness = 5 mm)

Fig. 5Ratio of different factor for different laminate thickness (v0= 200 m/s)

Ratio of different factor for different  laminate thickness (v0= 200 m/s)

5. Conclusions

An analytical model to predict residual velocity and ballistic limit after high velocity impact on woven Kevlar laminates is proposed. Four different energy absorption mechanisms for the laminate are considered in the model. Then, the model is simplified to obtain the residual velocity and ballistic limit. The analytical result is calculated to compare with that of experiment. The comparison shows a good correlation between the results obtained both numerically and experimental. The perturbation caused by tensile fiber failure descends as soon as the initial velocity or the thickness of laminate increases. The importance of energy absorbing by shear plugging, crushing and the momentum transfer will change by changing of initial velocity of projectile and thickness of laminate. The method of simplification could be useful in preliminary assessment of engineering.

References

  • Bhatnagar A. Lightweight Ballistic Composites Military and Law-Enforcement Applications. Abington Hall, Abington, Cambridge, 2006, p. 364-397.
  • Grujicic M., Arakere G. A ballistic material model for cross-plied unidirectional ultra-high molecular-weight polyethylene fibre-reinforced armor-grade composites. Materials Science and Engineering, Vol. 498, 2008, p. 231-241.
  • Nilakantan G., Keefe M., Gillespie Jr J., Bogetti T. Novel multi-scale modeling of woven fabric composites for use in impact studies. 10th International LSDYNA users conference, Dearborn, Michigan, USA, 2008.
  • Zhu G., Goldsmith Dharan W. C. H. K. Penetration of laminated Kevlar by projectiles-II, analytical model. International Journal of Solids and Structures, Vol. 29, Issue 4, 1992, p. 421-436.
  • Sun C. T., Potti S. V. A simple model to predict residual velocity of thick composite laminates subjected to high velocity impact. International Journal of Impact Engineering, Vol. 18, Issue 3, 1996, p. 339-353.
  • Wen H. M. Predicting the penetration and perforation of FRP laminates struck normally by projectiles with different nose shapes. Composite Structures, Vol. 49, Issue 3, 2000, p. 321-329.
  • Wen H. M. Penetration and perforation of thick FRP laminates. Composites Science and Technology, Vol. 61, 2001, p. 1163-1172.
  • Naik N. K., Shrirao P. Composite structures under ballistic impact. Composite Structures, Vol. 66, Issues 1-4, 2004, p. 579-590.
  • Xuan Haijun, Zhang Na Impact response and damage evolution of triaxial braided carbon/epoxy composites. Part I: Ballistic impact testing. Textile Research Journal, Vol. 83, Issue 16, 2013, p. 1703-1716.
  • Ganesh Babu M., Velmurugan Gupta R. N. K. Energy absorption and ballistic limit of targets struck by heavy projectile. Latin American Journal of Solids and Structures, Vol. 3, Issue 1, 2006, p. 21-39.
  • Sikarwar R. S., Velmurugan R., Madhu V. Experimental and analytical study of high velocity impact on Kevlar/Epoxy composite plates. Central European Journal of Engineering, Vol. 2, Issue 4, 2012, p. 638-650.
  • Sikarwar R. S., Velmurugan R. Ballistic impact on glass/epoxy composite laminates. Defense Science Journal, Vol. 64, Issue 4, 2014, p. 393-399.
  • Wang Y. B., Wang X. J., Hu X. Z., Wang F. Experimental study of ballistic resistance of Kevlar laminates. Engineering Mechanics, Vol. 22, Issue 3, 2005, p. 76-81.
  • Wang Y. B., Wang X. J., Bian L., Yu Y. M. CDM model and its application to numerical simulation on fiber-reinforced laminate under penetration. Explosion and Shock Wave, Vol. 28, Issue 2, 2008, p. 172-177.
  • Yahaya R., Sapuan S. M., Jawaid M., Leman Z., Zainudin E. S. Investigating ballistic impact properties of woven kenaf-aramid hybrid composites. Fibers and Polymers, Vol. 17, Issue 2, 2016, p. 275-281.
  • Gu Xingjin, Xu Xiwu, Huang Jing High velocity impact damage of thin composite laminates. Journal of Nanjing University of Aeronautics and Astronautics, Vol. 40, Issue 3, 2008, p. 370-375.
  • Michele Pasquali, Consuelo Terra, Paolo Gaudenzi Analytical modeling of high-velocity impacts on thin woven fabric composite targets. Composite Structure, Vol. 131, 2015, p. 951-965.
  • Chen Z., Yang G., Zhang C., Wang W. Experimental and modeling study on ballistic impact behavior of a woven carbon/epoxy composite. American Society for Composites Thirty-First Technical Conference, 2016.
  • Lopez Puente J., Zaera R., Navarro C. An analytical model for high velocity impacts on thin CFRPs woven laminated plates. International Journal of Solids and Structures, Vol. 44, 2007, p. 2837-2851.
  • Lopez Puente J., Zaera R., Navarro C. Analytical model for high velocity impaction thin CFRPs woven laminated targets. International Journal of Solids and Structures, 2007, p. 1-15.
  • Naik N. K., Shrirao P., Reddy B. C. K. Ballistic impact behavior of woven fabric composites: formulation. International Journal of Impact Engineering, Vol. 32, 2006, p. 1521-1552.
  • Naik N. K., Doshi A. V. Ballistic impact behavior of thick composites: parametric studies. Composite Structures, Vol. 82, 2008, p. 447-464.
  • Brenda Buitrago L., Shirley Garcia Castillo K., Enrique Barbreo Influence of shear plugging in the energy absorbed by thin carbon-fibre laminates subjected to high-velocity impacts. Composites Part B: Engineering, Vol. 49, 2013, p. 86-92.

Cited by

Numerical simulation of armor materials and optimization using gray relational analysis
Caio Lopes Martins Marques | Shiv Ranjan Kumar | Chandramani Goswami | Rahul Verma
(2020)

About this article

Received
08 October 2017
Accepted
27 October 2017
Published
31 December 2017
Keywords
high velocity impacts
composite laminates
residual velocity
ballistic limit
Acknowledgements

This research is partially supported by the National Natural Science Foundation of China (Grant No. 51475228), the Research Fund of State Key Laboratory of Mechanics and Control of Mechanical Structures (Nanjing University of Aeronautics and Astronautics) (Grant No. MCMS-0517K01) and the Priority Academic Program Development of Jiangsu Higher Education Institutions. The authors wish to thank Mr. Peng Haifeng, Mr. Wushuo and Mr. Li Jin of Chrdi (China) for their technical and financial support. Special thanks also to Mr. Jin Wanzeng and Mr. Mengxiangji for useful discussions.