Numerical analysis of intralaminar failure mechanisms in composite structures, Part I: FE implementation

This paper presents a three-dimensional continuum damage mechanics-based material model which was implemented in an implicit finite element code to simulate the progressive intralaminar degradation of fibre reinforced laminates. The damage model is based on ply failure mechanisms and uses seven damage variables assigned to tensile, compressive and shear damage at a ply level. Non-linear behaviour and irreversibility were taken into account and modelled. Some issues on the numerical implementation of the damage model are discussed and solutions proposed. Applications of the methodology are presented in Part II

Numerical analysis of intralaminar failure mechanisms in composite structures, Part II: Applications

A three-dimensional continuum damage mechanics-based material model was implemented in an implicit Finite Element code to simulate the progressive intralaminar degradation of fibre reinforced laminates based on ply failure mechanisms. This paper presents some structural applications of the progressive failure model implemented. The focus is on the non-linear response of the shear failure mode and its interaction with other failure modes. Structural applications of the damage model show that the proposed model is able to reproduce failure loads and patterns observed experimentally.

Predicting low-velocity impact damage on a stiffened composite panel

An intralaminar damage model, based on a continuum damage mechanics approach, is presented to model the damage mechanisms occurring in carbon fibre composite structures incorporating fibre tensile and compressive breakage, matrix tensile and compressive fracture, and shear failure. The damage model, together with interface elements for capturing interlaminar failure, is implemented in a finite element package and used in a detailed finite element model to simulate the response of a stiffened composite panel to low-velocity impact. Contact algorithms and friction between delaminated plies were included, to better simulate the impact event. Analyses were executed on a high performance computer (HPC) cluster to reduce the actual time required for this detailed numerical analysis. Numerical results relating to the various observed interlaminar damage mechanisms, delamination initiation and propagation, as well as the model’s ability to capture post-impact permanent indentation in the panel are discussed. Very good agreement was achieved with experimentally obtained data of energy absorbed and impactor force versus time. The extent of damage predicted around the impact site also corresponded well with the damage detected by non destructive evaluation of the tested panel.

Element-Free Galerkin modelling of composite damage

Delamination and matrix cracking are routine damage mechanisms, observed by post-mortem analysis of laminated structures containing geometrical features such as notches or bolts. Current finite element tools cannot explicitly model an intralaminar matrix microcrack, except if the location of the damage is specified a priori. In this work, a meshless technique, the Element-Free Galerkin (EFG) method, is utilized for the first time to simulate delamination (interlaminar) and intralaminar matrix microcracking in composite laminates.

An Experimental and Numerical Study of the Static and Fatigue Performance of a Composite Adhesive Repair

Experimental static and fatigue tension-tension tests were carried out on 5HS/RTM6 composite intact coupons and coupons incorporating adhesively-bonded (FM300-2) stepped flush joints. The results show that the adhesive joint, which is widely used in repairs, significantly reduces the static strength as well as the fatigue life of the composite. Both, the static and the fatigue failure of the ‘repaired’ coupons occur at the adhesive joint and involve crack initiation and propagation. The latter is modelled using interface finite elements based on the decohezive zone approach. The material degradation in the interface constitutive law is described by a damage variable, which can evolve due to the applied loads as well as the number of fatigue cycles. The fatigue formulation, based on a published model, is adapted to fit the framework of the pseudotransient formulation that is used as a numerical tool to overcome convergence difficulties. The fatigue model requires three material parameters. Numerical tests show that a single set of these parameters can be used to recover, very accurately, the experimental S-N relationship. Sensitivity studies show that the results are not mesh dependent.

Automatic Insertion of Cohesive Elements for Delamination

Composite damage modelling with cohesive elements has initially been limited to the analysis of interface damage or delamination. However, their use is also being extended to the analysis of inplane tensile failure arising from matrix or fibre fracture. These interface elements are typically placed at locations where failure is likely to occur, which infers a certain a priori knowledge of the crack propagation path(s). In the case of a crack jump for example, the location of the jump is usually not obvious, and the simulation would require the placement of cohesive elements at all element faces. A better option, presented here, is to determine the potential location of cohesive elements and insert them during the analysis. The aim of this work is to enable the determination of the crack path, as part of the solution process. A subroutine has been developed and implemented in the commercial finite element package ABAQUS/Standard[1] in order to automatically insert cohesive elements within a pristine model, on the basis of the analysis of the current stress field. Results for the prediction of delamination are presented in this paper.

3-D micromechanical model for predicting the elastic behaviour of woven laminates

This paper presents an analytical model for the prediction of the elastic behaviour of plain-weave fabric composites. The fabric is a hybrid plain-weave with different materials and undulations in the warp and weft directions. The derivation of the effective material properties is based on classical laminate theory (CLT).

The theoretical predictions have been compared with experimental results and predictions using alternative models available in the literature. Composite laminates were manufactured using the resin infusion under flexible tooling (RIFT) process and tested under tension and in-plane shear loading to validate the model. A good correlation between theoretical and experimental results for the prediction of in-plane properties was obtained. The limitations of the existing theoretical models based on classical laminate theory (CLT) for predicting the out-of-plane mechanical properties are presented and discussed. 

Numerical analysis of stiffener runout sections

The recent trend of incorporating more composite material in primary aircraft structures has highlighted the vulnerability of stiffened aerostructures to through-thickness stresses, which may lead to delamination and debonding at the skin-stiffener interface, leading to collapse. Stiffener runout regions are particularly susceptible to this problem and cannot be avoided due to the necessity to terminate stiffeners at rib intersections or at cutouts, interrupting the stiffener load path. In this paper, experimental tests relating to two different stiffener runout specimens are presented and the failure modes of both specimens are discussed in detail. A thinner-skinned specimen showed sudden and unstable crack propagation, while a thicker-skinned specimen showed initially unstable but subsequent stable crack growth. Detailed finite element models of the two specimens are developed, and it is shown how such models can explain and predict the behaviour and failure mode of stiffener runouts. The models contain continuum shell elements to model the skin and stiffener, while cohesive elements using a traction-separation law are placed at the skin-stiffener interface to effectively model the debonding which promotes structural failure.

Intralaminar toughness characterisation of unbalanced hybrid plain weave laminates

A numerical and experimental investigation on the mode-I intralaminar toughness of a hybrid plain weave composite laminate manufactured using resin infusion under flexible tooling (RIFT) process is presented in this paper. The pre-cracked geometries consisted of overheight compact tension (OCT), double edge notch (DEN) and centrally cracked four-point-bending (4PBT) test specimens. The position as well as the strain field ahead of the crack tip during the loading stage was determined using a digital speckle photogrammetry system. The limitation on the applicability of the standard data reduction schemes for the determination of intralaminar toughness of composite materials is presented and discussed. A methodology based on the numerical evaluation of the strain energy release rate using the J-integral method is proposed to derive new geometric correction functions for the determination of the stress intensity factor for composites. The method accounts for material anisotropy and finite specimen dimension effects regardless of the geometry. The approach has been validated for alternative non-standard specimen geometries. A comparison between different methods currently available for computing the intralaminar fracture toughness in composite laminates is presented and a good agreement between numerical and experimental results using the proposed methodology was obtained. 

A correction to the analytical solution of the mixed-mode bending (MMB) problem

This paper addresses the analytical solution of the mixed-mode bending (MMB) problem. The first published solutions used a load separation in pure mode I and mode II and were applied for a crack length less than the beam half-span, a <= L. In later publications, the same mode separation was used in deriving the analytical solution for crack lengths bigger than the beam half-span, a > L. In this paper it is shown that this mode separation is not valid when a > L and in some cases may lead to very erroneous results. The correct mode separation and the corresponding analytical solutions, when a > L, are presented. Results, of force vs. displacement and force vs. crack length graphs, obtained using the existing formulation and the corrected formulation are compared. A finite element solution, which does not use mode separation, is also presented

Structural testing and numerical simulation of a 34 m composite wind turbine blade

A full-scale 34 m composite wind turbine blade was tested to failure under flap-wise loading. Local displacement measurement equipment was developed and displacements were recorded throughout the loading history.

Ovalization of the load carrying box girder was measured in the full-scale test and simulated in non-linear FE-calculations. The nonlinear Brazier effect is characterized by a crushing pressure which causes the ovalization. To capture this effect, non-linear FE-analyses at different scales were employed. A global non-linear FE-model of the entire blade was prepared and the boundaries to a more detailed sub-model were extracted. The FE-model was calibrated based on full-scale test measurements.

Local displacement measurements helped identify the location of failure initiation which lead to catastrophic failure. Comparisons between measurements and FE-simulations showed that delamination of the outer skin was the initial failure mechanism followed by delamnination buckling which then led to collapse. 

An automated hybrid procedure for capturing mode-jumping in postbuckling composite stiffened structures

This paper presents a robust finite element procedure for modelling the behaviour of postbuckling structures undergoing mode-jumping. Current non-linear implicit finite element solution schemes, found in most finite element codes, are discussed and their shortcomings highlighted. A more effective strategy is presented which combines a quasi-static and a pseudo-transient routine for modelling this behaviour. The switching between these two schemes is fully automated and therefore eliminates the need for user intervention during the solution process. The quasi-static response is modelled using the are-length constraint while the pseudo-transient routine uses a modified explicit dynamic routine, which is more computationally efficient than standard implicit and explicit dynamic schemes. The strategies for switching between the quasi-static and pseudo-transient routines are presented

Stiffener debonding mechanisms in post-buckled CFRP aerospace panels

This paper describes the fractographic analysis of five CFRP post-buckled skin/stringer panels that were tested to failure in compression. The detailed damage mechanisms for skin/stiffener detachment in an undamaged panel were characterised and related to the stress conditions during post-buckling; in particular the sites of peak twist (at buckling nodes) and peak bending moments (at buckling anti-nodes). The initial event was intralaminar splitting of the +45 degrees plies adjacent to the skin/stiffener interface, induced by high twist at a nodeline. This was followed by mode II delamination, parallel to +/- 45 degrees plies and then lengthwise (0 degrees) shear along the stiffener centreline. The presence of defects or damage was found to influence this failure process, leading to a reduction in strength. This research provides an insight into the processes that control post-buckled performance of stiffened panels and suggests that 2D models and element tests do not capture the true physics of skin/stiffener detachment: a full 3D approach is required.

Predicting the ultimate load of a CFRP wingbox

The finite element method in conjunction with the Soutis-Fleck model is used to predict the residual strength after impact of a carbon-fibre reinforced plastic wingbox subjected to a cantilever type loading. The maximum stress failure criterion further validates the Soutis-Fleck model predictions. The Soutis-Fleck model predicts that the wingbox fails at a tip load of 99.2 kN, approximately 5.5% less than the experimental observation

Capturing mode-switching in postbuckling composite panels using a modified explicit procedure

A postbuckling blade-stiffened composite panel was loaded in uniaxial compression, until failure. During loading beyond initial buckling, this panel was observed to undergo a secondary instability characterised by a dynamic mode shape change. These abrupt changes cause considerable numerical difficulties using standard path-following quasi-static solution procedures in finite element analysis. Improved methods such as the arc-length-related procedures do better at traversing certain critical points along an equilibrium path but these procedures may also encounter difficulties in highly non-linear problems. This paper presents a robust, modified explicit dynamic analysis for the modelling of postbuckling structures. This method was shown to predict the mode-switch with good accuracy and is more efficient than standard explicit dynamic analysis. (C) 2003 Elsevier Science Ltd. All rights reserved.