Effect of clamping force on the delamination onset and growth in bolted composite laminates

Clamping force is a key element that alters the mechanism and sequence of failure in bolted joints of composite laminates. The mode of failure in bolted joints can be controlled by geometrical parameters and the preferred fail safe mode of failure is ‘bearing’ which generally consists of matrix cracks, delamination and fibre microbuckling. Three-dimensional (3-D) pinned (without clamping force) and bolted (1 kN clamping force) joint models were developed in [0/90]_s carbon fibre reinforced plastic (CFRP) laminates to show the clamping force effect on the onset and growth of delamination. It is shown that delamination was resulted from the shear stress components (Mode II & III) at the interface and the contribution of the out-of-plane component (Mode I - opening), so the clamping force, was negligible without modelling the in-plane failure modes and their coupling with delamination, which will be considered in future work.     

XFEM simulation of delamination in composite laminates

In this paper, the extended finite element method (XFEM) is extended to simulate delamination problems in composite laminates. A crack-leading model is proposed and implemented in the ABAQUS® to discriminate different delamination morphologies, i.e., the 0°/0° interface in unidirectional laminates and the 0°/90° interface in multidirectional laminates, which accounts for both interlaminar and intralaminar crack propagation. Three typical delamination problems were simulated and verified. The results of single delamination in unidirectional laminates under pure mode I, mode II, and mixed mode I/II correspond well with the analytical solutions. The results of multiple delaminations in unidirectional laminates are in good agreement with experimental data. Finally, using a recently proposed test that characterizes the interaction of delamination and matrix cracks in cross-ply laminates, the present numerical results of the delamination migration caused by the coupled failure mechanisms are consistent with experimental observations.     

Experimental study on delamination migration in composite laminates

The transition of delamination growth between different ply interfaces in composite tape laminates, known as migration, was investigated experimentally. The test method used promotes delamination growth initially along a 0/θ ply interface, which eventually migrates to a neighbouring θ/0 ply interface. Specimens with θ = 60° and 75° were tested. Migration occurs in two main stages: (1) the initial 0/θ interface delamination turns, transforming into intraply cracks that grow through the θ plies; this process occurs at multiple locations across the width of a specimen, (2) one or more of these cracks growing through the θ plies reaches and turns into the θ/0 ply interface, where it continues to grow as a delamination. A correlation was established between these experimental observations and the shear stress sign at the delamination front, obtained by finite element analyses.Overall, the experiments provide insight into the key mechanisms that govern delamination growth and migration.     

Finite element simulation of delamination growth in composite materials using LS-DYNA

In this paper, a modified adaptive cohesive element is presented. The new elements are developed and implemented in LS-DYNA, as a user defined material subroutine (UMAT), to stabilize the finite element simulations of delamination propagation in composite laminates under transverse loads. In this model, a pre-softening zone is proposed ahead of the existing softening zone. In this pre-softening zone, the initial stiffness and the interface strength are gradually decreased. The onset displacement corresponding to the onset damage is not changed in the proposed model. In addition, the critical energy release rate of the materials is kept constant. Moreover, the constitutive equation of the new cohesive model is developed to be dependent on the opening velocity of the displacement jump. The traction based model includes a cohesive zone viscosity parameter (η) to vary the degree of rate dependence and to adjust the maximum traction. The numerical simulation results of DCB in Mode-I is presented to illustrate the validity of the new model. It is shown that the proposed model brings stable simulations, overcoming the numerical instability and can be widely used in quasi-static, dynamic and impact problems.     

Predicting mode-II delamination suppression in z-pinned laminates

A finite element model for predicting delamination resistance of z-pin reinforced laminates under the mode-II load condition is presented. End notched flexure specimen is simulated using a cohesive zone model. The main difference of this approach to previously published cohesive zone models is that the individual bridging force exerted by z-pin is governed by a specific traction-separation law derived from a unit-cell model of single pin failure process, which is independent of the fracture toughness of the unreinforced laminate. Therefore, two separate traction-separation laws are employed; one represents unreinforced laminate properties and the other for the enhanced delamination toughness owing to the pin bridging action. This approach can account for the so-called large scale bridging effect and avoid using concentrated pin forces in numerical models, thus removing the mesh-size dependency and permitting more accurate and reliable computational solutions.     

Cryogenic delamination growth in woven glass/epoxy composite laminates under mixed-mode I/II fatigue loading

This paper investigates the fatigue delamination growth behavior in woven glass fiber reinforced polymer (GFRP) composite laminates under mixed-mode I/II conditions at cryogenic temperatures. Fatigue delamination tests were performed with the mixed-mode bending (MMB) test apparatus at room temperature, liquid nitrogen temperature (77 K) and liquid helium temperature (4 K), in order to obtain the delamination growth rate as a function of the range of the energy release rate, and the dependence of the delamination growth behavior on the temperature and the mixed-mode ratio of mode I and mode II was examined. The energy release rate was evaluated using three-dimensional finite element analysis. The fractographic examinations by scanning electron microscopy (SEM) were also carried out to assess the mixed-mode fatigue delamination growth mechanisms in the woven GFRP laminates at cryogenic temperatures.     

A cohesive zone model for predicting delamination suppression in z-pinned laminates

This paper presents a cohesive zone model based finite element analysis of delamination resistance of z-pin reinforced double cantilever beam (DCB). The main difference between this and existing cohesive zone models is that each z-pin bridging force is governed by a traction-separation law derived from a meso-mechanical model of the pin pullout process, which is independent of the fracture toughness of unreinforced laminate. Therefore, two different traction-separation laws are used: one representing the toughness of unreinforced laminate and the other the enhanced delamination toughness owing to the pin bridging action. This approach can account for the large scale bridging effect and avoid using concentrated pin forces, thus removing the mesh dependency and permitting more accurate analysis solution. Computations were performed using a simplified unit strip model. Predicted delamination growth and load vs. displacement relation are in excellent agreement with the prediction by a complete model, and both models are in good agreement with test measured load vs. displacement relation. For a pinned DCB specimen, the unit strip model can reduce the computing time by 85%.     

An interface element for the simulation of delamination in unidirectional fiber-reinforced composite laminates

Unidirectional fiber-reinforced composite laminates are widely used in aerospace industry for a great variety of structural parts. In order to enhance the exploitation of material reserves, there is a need for the integration of progressive damage scenarios in the design phase. Due to their hazardous effects on the load-carrying capacity of composite structures, this work focusses on the simulation of delaminations. A finite element based on a cohesive zone approach is developed. Two constitutive laws are proposed. One is characterized by linear degradation after delamination onset, the other is governed by exponential softening response. The damage process is history-dependent leading to an irreversible stiffness degradation in damaged zones. The practicability of the proposed model and the assets and drawbacks of the two material laws are shown by some numerical examples.     

Environmental effects on mode II fatigue delamination growth in an aerospace grade carbon/epoxy composite

An investigation of the effects of water, hydraulic fluid and deicing fluid exposure on mode II delamination propagation in an aerospace grade composite is presented. All exposed specimens suffered a loss in delamination toughness and an increase in fatigue delamination growth rate, which was particularly significant for deicing fluid exposure. The number of cycles for delamination onset was also reduced by these exposures, although scanning electron micrographs showed no significant differences between the fracture pattern of dry and exposed specimens. It was also shown that environmental effects can be simply accounted for in a cohesive zone based finite element model.     

Modelling damage evolution in composite laminates subjected to low velocity impact

In this paper, the impact damage of composite laminates in the form of intra- and inter-laminar cracking was modelled using stress-based criteria for damage initiation, and fracture mechanics techniques to capture its evolution. The nonlinear shear behaviour of the composite was described by the Soutis shear stress–strain semi-empirical formula. The finite element (FE) method was employed to simulate the behaviour of the composite under low velocity impact. Interface cohesive elements were inserted between plies with appropriate mixed-mode damage laws to model delamination. The damage model was implemented in the FE code (Abaqus/Explicit) by a user-defined material subroutine (VUMAT). Numerical results in general gave a good agreement when compared to experimentally obtained curves of impact force and absorbed energy versus time. The various damage mechanisms introduced during the impact event were observed by non-destructive technique (NDT) X-ray radiography and were successfully captured numerically by the proposed damage evolution model.     

Modeling delamination migration in cross-ply tape laminates

A novel modeling approach is proposed that combines the Floating Node Method (FNM) with the Virtual Crack Closure Technique (VCCT) to capture delamination migration in cross-ply tape laminates. Delamination migration is the damage process by which a delamination propagating at an interface relocates to a different interface via one or multiple matrix cracks. In the approach proposed, delamination, matrix cracking, and their interaction, are represented in a single element. The kinematics of both delamination and matrix cracks are represented explicitly. Migration onset location, and subsequent path, are determined as part of the solution, in a mesh-independent fashion. Delamination growth, matrix cracking, and migration onset, are all modeled using fracture mechanics based failure and migration criteria. The proposed approach is applied to the modeling of the Delamination Migration (DM) test, showing good qualitative and quantitative agreement with experiments.     

Intersonic delamination in curved thick composite laminates under quasi-static loading

Dynamic delamination in curved composite laminates is investigated experimentally and numerically. The laminate is 12-ply graphite/epoxy woven fabric L-shaped laminate subject to quasi-static loading perpendicular to one arm. Delamination initiation and propagation are observed using high speed camera and load–displacement data is recorded. The quasi-static shear loading initiates delamination at the curved region which propagates faster than the shear wave speed of the material, leading to intersonic delamination in the arms. In the numerical part, the experiments are simulated with finite element analysis and a bilinear cohesive zone model. Cohesive interface elements are used between all plies with the interface properties obtained from tests. The simulations predict a single delamination initiating at the corner under pure mode-I stress field propagating to the arms under pure mode-II stress field. The crack tip speeds transition from sub-Rayleigh to intersonic in conjunction with mode change. In addition to intersonic mode-II delamination, shear Mach waves emanating from the crack tips in the arms are observed. The simulations and experiments are found to be in good agreement at the macro-scale, in terms of load-displacement behavior and failure load, and at the meso-scale, in terms of delamination initiation location and crack propagation speeds. Finally, a mode dependent crack tip definition is proposed and observation of vibrations during delamination is presented. This paper presents the first conclusive evidence of intersonic delamination in composite laminates triggered under quasi-static loading.     

Cohesive zone modelling of delamination response of a composite laminate with interleaved nylon 6,6 nanofibres

This work simulates numerically Double Cantilever Beam and End Notched Flexure experiments on Carbon Fibre Epoxy Resin specimens that have been performed by some of the authors in a previous work. Specimens have been nanomodified by interleaving plies with a layer of electrospun nanofibres in the delaminated interface. Eight different configurations of nanofibres have been used as interleave, for a total of 9 configurations (8 nanomodified plus the virgin one) to be simulated for both kind of tests to identify the cohesive zone parameters corresponding to the effect of nanofibre diameter, nanolayer thickness and nanofibre orientation on the delamination behaviour of the composite.Results showed that a bilinear damage law is necessary for almost all nanomodified configurations, and presented a clear relationship between nanomat layer parameters and the cohesive energy of the interface.     

Finite element analysis of postbuckling and delamination of composite laminates using virtual crack closure technique

The two-dimensional and three-dimensional parametric finite element analysis (FEA) of composite flat laminates with two through-the-width delamination types: 0₄/(±θ)₆//0₄ and 0₄//(±θ)₆//0₄ (θ = 0°, 45°, and “//” denotes the delaminated interface) under compressive load are performed to explore the effects of multiple delaminations on the postbuckling properties. The virtual crack closure technique which is employed to calculate the energy release rate (ERR) for crack propagation is used to deal with the delamination growth. Three typical failure criteria: B-K law, Reeder law and Power law are comparatively studied for predicting the crack propagation. Effects of different mesh sizes and pre-existing crack length on the delamination growth and postbuckling properties of composite laminates are discussed. Interaction between the delamination growth mechanisms for multiple cracks for 0₄//(±θ)₆//0₄ composite laminates is also investigated. Numerical results using FEA are also compared with those by existing models and experiments.     

Cohesive zone length in numerical simulations of composite delamination

Accurate analysis of composite delamination using interface elements relies on having sufficient elements within a softening region known as the cohesive zone ahead of a crack tip. The present study highlights the limitations of existing formulae used to predict numerical cohesive zone length and demonstrates modifications necessary for improved accuracy. Clarification is also provided regarding the minimum number of interface elements within the cohesive zone. Finally, appropriate values of numerical interfacial strength are examined. The results presented will aid the application of mesh design techniques that both preserve numerical accuracy, whilst minimising computational expense.     

Effect of damage levels and curing stresses on delamination growth behaviour emanating from circular holes in laminated FRP composites

This paper deals with the study of the effect of drilling induced delamination damage levels and residual thermal stresses (developed during manufacturing process of cooling the laminate form curing temperature to room temperature) on delamination growth behaviour emanating form circular holes in graphite/epoxy laminated FRP composites. Two sets of full three dimensional finite element analyses (one with the residual thermal stresses developed while curing the laminate and the other without residual thermal stresses i.e. with mechanical loading only) have been performed to calculate the displacements and interlaminar stresses along the delaminated interfaces responsible for the delamination onset and propagation. Modified crack closure integral (MCCI) techniques based on the concepts of linear elastic fracture mechanics (LEFM) have been used to calculate the distribution of individual modes of strain energy release rates (SERR) to investigate the interlaminar delamination initiation and propagation characteristics. Asymmetric variations of SERR obtained along the delamination front are caused by the overlapping stress fields due to the coupling effect of thermal and mechanical loadings. It is found that parameters such as ply orientation, drilling induced damage levels and material heterogeneity at the delaminated interface dictate the interlaminar fracture behaviour of laminated FRP composites.     

Tensile failure of laminates containing an embedded wrinkle; numerical and experimental study

The failure of a quasi-isotropic composite laminate containing an embedded out-of-plane fibre wrinkle defect was investigated under tension loading. Laboratory test specimens with controlled severity of fibre waviness were manufactured. Along with recording load–displacement data, high resolution camera images were taken at regular intervals which monitored the initiation and interaction of different damage mechanisms during test. Three-dimensional FE models were built following the geometry of actual test specimens. The information obtained from the tests was used to develop user material subroutines, implemented in Abaqus/Explicit as continuum damage and cohesive zone models for intra- and inter-ply failure respectively. The results of the simulations showed very good correlation with test observations in terms of failure load, location of damage initiation and interaction between different damage mechanisms for a range of waviness cases tested.     

A numerical investigation into the size effect in the transverse crack tension test for mode II delamination

In experimental studies, a size effect has been measured for the fracture energy in the transverse crack tension test. This paper presents a numerical investigation into the cause of this size effect. A finite element model has been developed that includes delamination, friction and shear nonlinearity. After calibration of the model, the size effect was reproduced well. It is shown that shear nonlinearity and friction separately contribute to the measured size effect and that significant amplification of the size effect takes place because of their interaction. As a consequence of their interaction, the unstable crack growth that was observed for the thicker specimens in the experiments is also present in the model results.     

Studies of mode I delamination in monotonic and fatigue loading using FBG wavelength multiplexing and numerical analysis

Bridging by intact fibers in composite materials is one of the most important toughening mechanisms. However, a direct experimental assessment of its contribution is not easy to achieve. In this work a semi-experimental method is proposed to quantify its contribution to fracture of unidirectional carbon fiber/epoxy double cantilever beam (DCB) specimens in mode I delamination under monotonic and 1 Hz fatigue loads. In each specimen, an embedded optical fiber with an array of eight wavelength-multiplexed fiber Bragg gratings is used to measure local strains close to the crack plane. The measured strain distribution serves in an inverse identification procedure to determine the tractions in the bridging zone in monotonic and fatigue loads. These tractions are used to calculate the energy release rate (ERR) associated with bridging fibers. The results indicate that the ERR due to bridging is about 40% higher in fatigue. The bridging tractions are further included in a cohesive element model which allows to predict precisely the complete load displacement curve of monotonic DCB tests. Using the principle of superposition and the identified tractions, the total stress intensity factor (SIF) is calculated. The results show that the SIF, at initiation, is very close to the one calculated at crack propagation and bridging by intact fibers is responsible for the entire increase in toughness seen in the DCB specimens used herein.     

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.