Numerical study of the structural behaviour of impacted composite laminates subjected to compression load

Composite materials allow all the benefits which a high specific strength involves, in a design process their application involves many critical problems. Currently, these problems, such as environmental conditions, notch sensitivity, damaging under low velocity impacts, are taken into account by means of the application of conservative design safety factors regarding the ultimate tensile strength. In order to try to reduce these safety factors, this work aimed to study and to understand the impact damage growing mechanisms due to compression loads. To this purpose, compression tests have been experimentally performed on composite panels, which have been previously subjected to low velocity impact phenomena, considering impact energies of 6 J, 10 J and 13 J respectively. Moreover, numerical model able to simulate Low Velocity Impacts (LVI) and Compression After Impacts (CAI) onto CFRP panels is proposed. A single explicit finite element analysis has been carried out by using the Abaqus^® finite element code; the need to build a numerical model, which allows simulation in only one analysis both LVI and CAI steps, depends on the difficulty to import the impact damage distribution into a separate compression analysis. In fact, in only one analysis the compression step can occur directly onto the impacted plate, which allows to consider the effective impact damage distribution as the starting configuration for quasi static analysis under operating loads.     

A progressive failure model for mesh-size-independent FE analysis of composite laminates subject to low-velocity impact damage

An original, ply-level, computationally efficient, three-dimensional (3D) composite damage model is presented in this paper, which is applicable to predicting the low velocity impact response of unidirectional (UD) PMC laminates. The proposed model is implemented into the Finite Element (FE) code ABAQUS/Explicit for one-integration point solid elements and validated against low velocity impact experimental results.     

Predicting low velocity impact damage and Compression-After-Impact (CAI) behaviour of composite laminates

Low-velocity impact damage can drastically reduce the residual strength of a composite structure even when the damage is barely visible. The ability to computationally predict the extent of damage and compression-after-impact (CAI) strength of a composite structure can potentially lead to the exploration of a larger design space without incurring significant time and cost penalties. A high-fidelity three-dimensional composite damage model, to predict both low-velocity impact damage and CAI strength of composite laminates, has been developed and implemented as a user material subroutine in the commercial finite element package, ABAQUS/Explicit. The intralaminar damage model component accounts for physically-based tensile and compressive failure mechanisms, of the fibres and matrix, when subjected to a three-dimensional stress state. Cohesive behaviour was employed to model the interlaminar failure between plies with a bi-linear traction–separation law for capturing damage onset and subsequent damage evolution. The virtual tests, set up in ABAQUS/Explicit, were executed in three steps, one to capture the impact damage, the second to stabilize the specimen by imposing new boundary conditions required for compression testing, and the third to predict the CAI strength. The observed intralaminar damage features, delamination damage area as well as residual strength are discussed. It is shown that the predicted results for impact damage and CAI strength correlated well with experimental testing without the need of model calibration which is often required with other damage models.     

Damage prediction in composite sandwich panels subjected to low-velocity impact

The paper illustrates the application of a finite element tool for simulating the structural and damage response of foam-based sandwich composites subjected to low-velocity impact. Onset and growth of typical damage modes occurring in the composite skins, such as fibre fracture, matrix cracking and delaminations, were simulated by the use of three-dimensional damage models (for intralaminar damage) and interfacial cohesive laws (for interlaminar damage). The nonlinear behaviour of the foam core was simulated by a crushable foam plasticity model. The FE results were compared with experimental data acquired by impact testing on sandwich panels consisting of carbon/epoxy facesheets bonded to a PVC foam. Good agreement was obtained between predictions and experiments in terms of force histories, force–displacement curves and dissipated energy. The proposed model was also capable of simulating correctly nature and size of impact damage, and of capturing the key features of individual delaminations at different depth locations.     

Analytical and experimental studies on the low-velocity impact response and damage of composite laminates under in-plane loads with structural damping effects

In this paper, low-velocity impact response and damage of composite laminates under in-plane loads are analytically and experimentally investigated. The authors recently proposed a modified displacement field of plate theory, considering the effect of initially loaded in-plane strain, and used a finite element program to analyze the structural behavior of the composite laminate. In this study, the program is upgraded to account for the structural damping effect of the laminate. A pendulum type impact test system and an in-plane loading fixture are constructed for the experimental study. The analytical and experimental impact behaviors are compared at different impact energy levels for cases with an initial in-plane tensile load and a compressive load, as well as cases without the initial in-plane load. The results show good correspondence between the analytical and experimental impact force histories. The effect of the initial in-plane load reduces for higher impact energies. The numerical estimation of the damaged area is in good agreement with the results from C-scanning experiments.     

The effect of hybridization on the ballistic impact behavior of hybrid composite armors

In the present study, effect of hybridization on the hybrid composite armors under ballistic impact is investigated using hydrocode simulations. The hybrid composite armor is constructed using various combinations and stacking sequences of fiber reinforced composites having woven form of fibers specifically high specific-modulus/high specific-strength Kevlar fiber (KF), tough, high strain-to-failure fiber Glass fiber (GF) and high strength/high stiffness Carbon fiber (CF). Different combinations of composite armors studied are KF layer in GF laminate, GF layer in KF laminate, KF layer in CF laminate and CF layer in KF laminate at various positions of hybridized layers for a fixed thickness of the target. In this article the results obtained from the finite element model are validated for the case of KF layer in a GF laminate with experimental predictions reported in the literature in terms of energy absorption and residual velocity and good agreement is observed. Further, the effect of stacking sequence, projectile geometry and target thickness on the ballistic limit velocity, energy absorbed by the target and the residual velocity are presented for different combinations of hybrid composite armors. The simulations show that, at a fixed thickness of the hybrid composite armor, stacking sequence of hybridized layer shows significant effect on the ballistic performance. The results also indicate energy absorption and ballistic limit velocity are sensitive to projectile geometry. Specifically, it is found that arranging the KF layer at the rear side, GF layer in the exterior and CF layer on the front side offers good ballistic impact resistance. The hybrid composite armor consisting of a CF layer in KF laminate acquires maximum impact resistance and is the best choice for the design compared to that of other combinations studied.     

A methodology for hydrocode analysis of ultra-high molecular weight polyethylene composite under ballistic impact

Ballistic performance analysis of ultra-high molecular weight polyethylene (UHMW-PE) is critical for the design of armour systems against ballistic threats. However, no validated modelling strategy has been published in literature for UHMW-PE composite that captures the penetration and damage mechanisms of thick targets impacted between 900 m/s and 2000 m/s. Here we propose a mechanistically-based and extensively validated methodology for the ballistic impact analysis of thick UHMW-PE composite. The methodology uses a non-linear orthotropic continuum model that describes the composite response using a non-linear equation of state (EoS), orthotropic elastic–plastic strength with directional hardening and orthotropic failure criteria. A new sub-laminate discretisation approach is proposed that allows the model to more accurately capture out-of-plane failure. The model is extensively validated using experimental ballistic data for a wide range of UHMW-PE target thicknesses up to 102 mm against 12.7 mm and 20 mm calibre fragment simulating projectiles (FSPs) with impact velocities between 400 m/s and 2000 m/s. Very good overall agreement with experimental results is seen for depth of penetration, ballistic limit and residual velocity, while the penetration mechanisms and target bulge behaviour are accurately predicted. The model can be used to reduce the volume of testing typically required to design and assess thick UHMW-PE composite in ballistic impact applications.     

Response of pultruded composite tubes subjected to dynamic and impulsive axial loading

The energy absorption of circular pultruded composite tubes subjected to axial crush load, transmitted by a small attached mass accelerated by means of an explosive load is presented in this paper. Different masses of explosive are used to provide a range of transmitted impulse and crushed distance of the pultruded composite tubes. The influence of the mass of the explosive on the tube response is investigated with regard to crushed distance, the average crushing force and the specific energy absorption (SEA). The crushing distance increases with increasing transmitted impulse. The results and failure mode are also compared with compression tests carried out on a servo-hydraulic machine (type: MTS-309).     

A unit cell approach of finite element calculation of ballistic impact damage of 3-D orthogonal woven composite

This paper presents ballistic impact damages of 3-D orthogonal woven composite in finite element analysis (FEA) and experimental. A unit-cell model of the 3-D woven composite was developed to define the material behavior and failure evolution. A user-defined subroutine VUAMT was compiled and connected with commercial available FEA code ABAQUS/Explicit to calculate the ballistic penetration. Ballistic impact tests were conducted to investigate impact damage of 3-D kevlar/glass hybrid woven composite. Residual velocities of conically-cylindrical steel projectiles (Type 56 in China Military Standard) and impact damage of the composite targets after ballistic perforation were compared both in theoretical and experimental. The reasonable agreements between FEA results and experimental results prove the validity of the unit-cell model in ballistic limit prediction of the 3-D woven composite. We believe such an effort could be extended to bulletproof armor design with the 3-D woven composite.     

Investigation of the response to low velocity impact and quasi-static indentation loading of particle-toughened carbon-fibre composite materials

This work investigates damage caused by low velocity impact and quasi-static indentation loading in four different particle-toughened composite systems, and one untoughened system. For impact tests, a range of energies were used between 25 and 50 J. For QSI, coupons were interrupted at increasing loading point displacement levels from 2 to 5 mm to allow for monitoring of damage initiation and propagation. In both loading cases, non-destructive inspection techniques were used, consisting of ultrasonic C-scan and X-ray micro-focus computed tomography. These techniques are complemented with instrumentation to capture force–displacement data, whereby load-drops are associated with observed damage modes. Key results from this work highlight particular issues regarding strain-rate sensitivity of delamination development and an earlier onset of fibre fracture associated with particle-toughened systems. These issues, in addition to observations on the role of micro-scale events on damage morphology, are discussed with a focus on material development and material testing practices.     

A numerical study of ply orientation on ballistic impact resistance of multi-ply fabric panels

This paper presents a detailed finite element (FE) analysis aiming to investigate numerically the impact deformation of multi-ply fabric panels with angled plies. The purpose of the investigation described in this paper is to study numerically the way in which the multi-ply panels deform and to identify the energy absorption in different panel constructions. The FE model was created using ABAQUS to simulate the transverse impact of a projectile onto various woven fabric panels. Influencing factors such as the impact velocity, panel construction and the number of plies are taken into account in the FE simulations. The numerical predictions show that the orientation of plies significantly affects the energy-absorbing capacity of the multi-ply fabric panels. The angled panels always increase the energy-absorbing capacity, compared with the aligned panel, by as much as 20%, depending on the number of plies in the panel. In addition, the stacking sequence of oriented plies also plays an important role in absorbing the energy. For the multi-ply fabric panel with large numbers of plies, there is an optimised sequence of plies which can maximise the energy-absorbing capacity of the panel. An important aspect of the work is validation of the numerical technique. It is shown that the FE predictions are highly consistent with the experimental study [1].     

Experimental and numerical study of the ballistic performance of steel fibre-reinforced explosively welded targets impacted by a spherical fragment

Steel fibres were used to reinforce the layered targets with surface-to-surface combination. The two- and three-layer metal targets with a total thickness of 5 mm were fabricated by explosive welding. The damage mechanism and the anti-penetration performance of the targets were studied experimentally and numerically using the LS-DYNA 3D finite element code. The effects of layer number and fibre spacing density on the anti-penetration performance were discussed. The results show that the failure modes of the steel front plate were shearing and plugging, and that the failure mode of the aluminium rear plate was ductile prolonging deformation when the tied interface failed by tension (or shearing and plugging when the interface remained connected) for the two-layer target. For the three-layer target, the failure modes of the steel front plate and the aluminium middle plate were shearing and plugging, while the steel rear plate failed by ductile prolonging deformation. At the same time, the steel-fibres failed by bending and tensile deformation. The anti-penetration performance of the three-layer composite targets was better compared with the performance of the two-layer targets when the areal density and fibre spacing density were equal. The reinforced fibres will improve the anti-penetration performance of the targets, and the ballistic resistance decreased with an increase in the fibre spacing distance.     

Experimental and numerical investigations of low velocity impact on functionally graded circular plates

This paper addresses low-velocity impact behaviour of functionally graded clamped circular plates. An experimental work was carried out to investigate the impact behaviour of FG circular plates which is composed of ceramic (SiC) and metal (Al) phases varying through the plate thickness by using a drop-weight impact test system. The influence of the compositional gradient exponent and impactor velocity on the contact forces and absorbed energies was concentrated on the tests. The explicit finite element method, in which a volume fraction based elastic–plastic model (the TTO model) was implemented for the functionally graded materials, was used to simulate their drop-weight impact tests. Effective material properties at any point inside FGM plates were determined using Mori–Tanaka scheme. The experimental and numerical results indicated that the compositional gradient exponent and impactor velocity more effective on the elasto-plastic response of the FG circular plates to a low-velocity impact loading. The comparison at the theoretical and experimental results showed that the use of the TTO model in modelling the elasto-plastic behaviour of FG circular plates results in increasing deviations between the numerical and experimental contact forces for ceramic-rich compositions whereas it becomes more successful for metal-rich compositions.     

Finite element modeling of ballistic impact on multi-layer Kevlar 49 fabrics

This paper presents a material model suitable for simulating the behavior of dry fabrics subjected to ballistic impact. The developed material model is implemented in a commercial explicit finite element (FE) software LS-DYNA through a user defined material subroutine (UMAT). The constitutive model is developed using data from uniaxial quasi-static and high strain rate tension tests, picture frame tests and friction tests. Different finite element modeling schemes using shell finite elements are used to study efficiency and accuracy issues. First, single FE layer (SL) and multiple FE layers (ML) were used to simulate the ballistic tests conducted at NASA Glenn Research Center (NASA-GRC). Second, in the multiple layer configuration, a new modeling approach called Spiral Modeling Scheme (SMS) was tried and compared to the existing Concentric Modeling Scheme (CMS). Regression analyses were used to fill missing experimental data – the shear properties of the fabric, damping coefficient and the parameters used in Cowper-Symonds (CS) model which account for strain rate effect on material properties, in order to achieve close match between FE simulations and experimental data. The difference in absorbed energy by the fabric after impact, displacement of fabric near point of impact, and extent of damage were used as metrics for evaluating the material model. In addition, the ballistic limits of the multi-layer fabrics for various configurations were also determined.     

Multiaxial fatigue life prediction of composite bolted joint under constant amplitude cycle loading

In this paper, a fatigue model of composite is established to predict multiaxial fatigue life of composite bolted joint under constant amplitude cycle loading. Firstly, finite element model is adopted to investigate stress state of composite bolted joint under constant amplitude cycle loading. Secondly, Tsai–Hill criterion is used to calculate equivalent stress of joint. At last, modified S–N fatigue life curve fitted by unidirectional laminate S–N curve which takes ply angle and stress ratio into consideration is adopted to determine fatigue life of composite. Calculation results of equivalent stress model show excellent agreement with experiments of composite bolted joint.     

Experimental investigation of impact on composite laminates with protective layers

This paper presents an experimental study of low energy impacts on composite plates covered with a protective layer. In service, composite materials are subjected to low energy impacts. Such impacts can generate damage in the material that results in significant reduction in material strength. In order to reduce the damage severity, one solution is to add a mechanical protection on composite structures. The protection layer is made up of a low density energy absorbent material (hollow spheres) of a certain thickness and a thin layer of composite laminate (Kevlar). Energy absorption ability of these protective layers can be deduced from the load/displacement impact curves. First, two configurations of protection are tested on an aluminium plate in order to identify their performance against impact, then the same are tested on composite plates. Test results from force–displacement curves and C-scan control are compared and discussed and finally a comparison of impact on composite plates with and without protection is made for different configurations.     

Sensitivity of structurally loaded sandwich panels to localized ballistic penetration

This paper describes a series of tests focused on the combination of structural loading (bending, shear) and simultaneous penetrating impact on sandwich panels with thin GFRP face-sheets, with emphasis on the specific damage morphologies and developments depending on the type and magnitude of structural loading. The test specimens were sandwich panels, length 250 mm and width 150 mm, with carbon fibre prepreg face-sheets ([0°/90°], thickness t_f ≅ 0.5 mm) bonded to the faces of a foam core (density 80 kg/m³, thickness H = 10 mm). The impact velocity was approximately 420 m/s, using a spherical steel impactor, diameter 10 mm, with a mass of 4.1 g. A high-speed camera was used for registration of panel response. It was demonstrated, that, at preload levels above a specific limit, the impact would cause catastrophic failure, i.e., complete or near-complete loss of structural load carrying capacity. Developments of failure morphology, consistent with the observed evidence, were derived and outlined.     

Detection of multiple low-energy impact damage in composite plates using Lamb wave techniques

This work assesses the suitability of the two zero-order Lamb wave modes to detect multiple barely-visible impact damage in composite material. Three specimens were subjected to damage at three different low-energy levels and one was left as an undamaged reference sample. Ultrasonic Lamb wave modes were selectively generated by surface-bonded piezoceramic wafer transducers in two tuned configurations. By using an algorithm based on the Akaike Information Criterion the time-of-flight of the Lamb modes was determined, allowing their threshold detection capabilities for the studied application to be successfully benchmarked. The results were consistently validated by digital shearography, ultrasonic C-scan and optical microscopy. A study of the effects on structural integrity was completed with an assessment of the damping ratio and residual bending strength proving to be sensitive parameters to the induced damage.     

A numerical investigation into the influence of fabric construction on ballistic performance

The use of high-performance fibres has made it possible to produce lightweight and strong personal body armour. Parallel to the creation and use of new fibres, fabric construction also plays an essential role for performance improvement. In this research, finite element (FE) models were built up and used to predict the response of woven fabrics with different structural parameters, including fabric structure, thread density of the fabric and yarn linear density. The research confirmed that the plain woven fabric exhibits superior energy absorption over other structures in a ballistic event by absorbing 34% more impact energy than the fabric made from 7-end satin weave. This could be explained that the maximum interlacing points in this fabric which help transmit stress to a larger fabric area, enabling more secondary yarns to be involved for energy dissipation. It was found that fabric energy absorption decreases as fabric is made denser, and this phenomenon becomes more pronounced in a multi-ply ballistic system than in a single-ply system. The research results also indicated that the level of yarn crimp in a woven fabric is an effective parameter in influencing the ballistic performance of the fabrics. A low level of yarn crimp would lead to the increase of the fabric tensile modulus and consequently influencing the propagation of the transverse wave. In addition, it was found that for fabrics with the same level of yarn crimp, low yarn linear density and high fabric tightness were desirable for ballistic performance improvement.