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.     

Reinforcement delamination of metallic beams strengthened by FRP strips: Fracture mechanics based approach

The delamination failure of metallic beams reinforced by externally bonded fibres reinforced polymers (FRP) is addressed in this paper and a simplified fracture mechanics based approach for the edge delamination of the reinforcement strips is illustrated. The criterion is based on the evaluation of the energy release rate (ERR) using both analytical and numerical models. The analytical models consist of a simplified version of a “two parameters elastic foundation” and “transformed section” model while the numerical analyses refer to the modified virtual crack closure technique (MVCCT). The main aim of the paper is to establish a fracture mechanics failure criterion based on the ERR and the specific fracture energy of the bonded strips. The criterion is an alternative approach to the well known stress based method to asses the load carrying capacity of the adhesive joint. The accuracy of the simplified approaches is shown through a numerical example which refers to a steel beam strengthened by carbon fibres reinforced polymers (CFRP).     

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%.     

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.     

Modelling delamination onset and growth in pin loaded composite laminates

A three-dimensional (3D) finite element (FE) model is created with cohesive zone elements (CZE) to simulate a mechanically fastened [0°/90°]_s pin-loaded joint in a composite laminate. The model incorporates fully integrated solid elements in the pin-loaded area to accurately capture the high stress gradients. Contact based cohesive elements with a bilinear traction–separation law are inserted between the layers to capture the onset and growth of delamination. The stress distribution around the pin-loaded hole was verified with the widely used cosine stress distribution model. Results from the FE model show that delamination damage initiated at the point of maximum average shear stress at the 0°/90° interface. The delaminated area develops an elliptical shape which grows in a non-self similar manner with increasing pin displacement. It is concluded that a progressive damage model should be included to provide a full understanding of the failure sequence, work that the authors are currently engaged with.     

On the robustness of finite element procedures based on Virtual Crack Closure Technique and fail release approach for delamination growth phenomena. Definition and assessment of a novel methodology

Numerical procedures based on the combined use of the Virtual Crack Closure Technique and of a fail release approach have been widely used to simulate delamination growth phenomena of composite material structures. This paper starts explaining why this kind of methodologies might not be robust due to mesh and load step size dependency and introduces a novel approach able to cope with the problems identified. Finally the effectiveness and robustness of the proposed procedure, implemented into a commercial finite element software by means of user subroutines, are assessed by comparing the obtained numerical results for a delamination growth phenomenon against literature experimental data on a stiffened panel with a circular embedded delamination under compressive load.     

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.     

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.     

A critical plane-based fracture criterion for mixed-mode delamination in composite materials

A new critical plane-based mixed-mode delamination failure criterion is proposed in this study. First, many existing models are reviewed and their capability to handle the mixed-mode fracture of general anisotropic materials are discussed. Following this, a previously developed critical plane approach is extended to analyze the interfacial fracture of composite materials by considering the anisotropic fracture resistance under mixed-mode loadings. Next, comparison with extensive experimental data available in the literature is performed to demonstrate the validity of the proposed criterion. A general good agreement is observed between the model's predictions and experimental observations. Finally, some conclusions and future work are drawn based on the proposed study.     

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.     

Debonding failure modes of flexural FRP-strengthened RC beams

In this paper, different types of debonding failure modes are described. Then, experimental results of four-point bending tests on FRP strengthened RC beams are presented and debonding failure mechanisms of strengthened beams are investigated using analytical and finite element solutions. Reasonable results could be obtained for modelling of debonding failure load of tested beams.

Existing international codes and guidelines from organizations such as ACI, fib, ISIS, JSCE, SIA, TR55, etc. are presented and compared with the results from the experiments and calculations. A discrepancy of up to 250% was seen between different codes and guidelines for predicting the debonding load. Furthermore, a new recommendation for debonding control is given.     

Composite delamination modelling using a multi-layered solid element

This paper focuses on the latest development of a solid hexahedron element for composite delamination analysis. The 8-node solid is derived from a 20-node hexahedron. It is transformed into two physical independent 4-node shell elements according to the propagation of delamination process within the element.     

Formulation and assessment of an enhanced finite element procedure for the analysis of delamination growth phenomena in composite structures

An existing procedure based on the combined use of the Virtual Crack Closure Technique and of a fail release approach for the analysis of delamination growth phenomena in composite structures has been enhanced with a front-tracing algorithm and suitable expressions for the evaluation of the Strain Energy Release Rate when dealing with non-smoothed delamination fronts. The enhanced procedure has been implemented into a commercial finite element software by means of user subroutines and applied to the analysis of a composite stiffened panel with an embedded delamination under compressive load. The effectiveness and robustness of the enhanced procedure have been assessed by comparing literature experimental data and numerical results obtained by using different mesh densities in the damaged area (global/local approach).     

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.     

Interlaminar fracture analysis in the GII-GIII plane using prestressed transparent composite beams

This work presents the mixed-mode II/III prestressed split-cantilever beam specimen for the fracture testing of composite materials. In accordance with the concept of prestressed composite beams one of the two fracture modes is provided by the prestressed state of the specimen, and the other one is increased up to fracture initiation by using a testing machine. The novel beam-like specimen is able to provide any combinations of the mode-II and mode-III ERRs. Data reduction is made by using the virtual crack-closure technique. The applicability and the limitations of the novel fracture mechanical test are demonstrated using unidirectional glass/polyester composite specimens. If only crack propagation onset is involved then the mixed-mode beam specimen can be used to obtain the fracture criterion of transparent composite materials in the G_II-G_III plane in a relatively simple way.     

Multi-objective design optimization of variable stiffness composite cylinders

The bending-induced buckling improvement in a variable stiffness (VS) composite cylinder (made by fiber steering) is studied. For such a cylinder, the effect of the variation of the direction of the load on its buckling performance of the cylinder is also examined. Compromise programming, as a multi-objective optimization method, is used to design for buckling of the VS cylinder subjected to bending load in either of the two opposite directions. Different combinations of weight factors for the structural performance in the two opposite directions were also applied to obtain the Pareto frontier as the main decision making tool for the designers in a multi-objective design problem.     

FEM modeling of multilayered textile composites based on shell elements

The paper presents a shell element based unit cell approach for numerical homogenization of fiber reinforced textile laminates. The modeling strategy is set up within the framework of the Finite Element Method. Multilayer laminates comprising equal weaves are considered and the constituents, i.e. the tows as well as the unreinforced matrix pockets are discretized by shell elements only which are coupled appropriately. A study on the effective extensional laminate-shell stiffnesses is presented, the results are discussed, and are compared to approaches found in the literature. Additionally, geometrically nonlinear simulations are conducted and the results are compared with experimental tests from literature.     

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.