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.

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. 

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.

Fracture mechanics using a 3D composite element

The formulation of a 3D composite element and its use in a mixed-mode fracture mechanics example is presented. This element, like a conventional 3D finite element, has three degrees of freedom per node although, like a plate element, the strains are defined in the local directions of the mid-plane surface. The stress-strain property matrix of this element was modified to decouple the stresses in the local mid-plane and the strains normal to this plane thus preventing the element from being too stiff in bending. A main advantage of this formulation is the ability to model a laminate with a single 3D element. The motivation behind this work was to improve the computational efficiency associated with the calculation of strain energy release rates in laminated structures. A comparison of mixed-mode results using different elements of an in-house finite element package are presented. Good agreement was achieved between the results obtained using the new element and coventional higher-order elements

Buckling mode transition in hat-stiffened composite panels loaded in uniaxial compression

A combined experimental and analytical study of a hat-stiffened carbon-fibre composite panel loaded in uniaxial compression was investigated. A buckling mode transition was observed in the panel's skin bay which was not captured using non-linear finite-element analysis. Good correlation between experimental and numerical strain and displacement results was achieved in the prebuckling and initial postbuckling region of the loading history. A Marguerre-type Rayleigh-Ritz energy method was applied to the skin bay using representative displacement functions of permissible mode shapes to explain the mode transition phenomenon. The central criterion of this method was based on the assumption that a change in mode shape occurred such that the total potential energy of the structure was maintained at a minimum. The ultimate strength of the panel was limited by the column buckling strength of the hat-stiffeners.

Nonlinear numerical modelling of lightning strike effect on composite panels with temperature dependent material properties

This paper presents a physics based modelling procedure to predict the thermal damage of composite material when struck by lightning. The procedure uses the Finite Element Method with non-linear material models to represent the extreme thermal material behaviour of the composite material (carbon/epoxy) and an embedded copper mesh protection system. Simulation predictions are compared against published experimental data, illustrating the potential accuracy and computational cost of virtual lightning strike tests and the requirement for temperature dependent material modelling. The modelling procedure is then used to examine and explain a number of practical solutions to minimize thermal material damage. © 2013 Elsevier Ltd.

NONLINEAR NUMERICAL MODELLING OF LIGHTENING STRIKE EFFECT ON COMPOSITE PANELS WITH TEMPERATURE DEPENDANT MATERIAL PROPERTIES

Lightning strike is one of the challenges that the aerospace industry is facing in an effort to increase the percentage of composite materials used in aircraft structures. Lightning strike damage is due to high orthotropic electric resistivity of the composite panels, which leads to high thermal loads that cause decomposition of the epoxy and delimitations of the laminates. Yet, experimental testing of lightning strike on aircraft panels is expensive due to the large number of design parameters that can control the inflicted damage. A coupled thermal-electrical finite element analysis is used to investigate the design variables space that can affect lightning strike damage on epoxy/graphite composite panels. The contribution of this study is modeling the composite panels’ material properties as temperature dependent, which was excluded by other researchers. A number of practical solutions to minimize the damage effect are proposed. Two set of experimental results are used to verify the numerical ones. One experimental set for plain composite panel, and second one for composite panels with joints

Study of localized damage in composite laminates using micro-macro approach

A robust multiscale scheme referred to as micro–macro method has been developed for the prediction of localized damage in fiber reinforced composites and implemented in a finite element framework. The micro–macro method is based on the idea of partial homogenization of a structure. In this method, the microstructural details are included in a small region of interest in the structure and the rest is modeled as a homogeneous continuum. The solution to the microstructural fields is then obtained on solving the two different domains, simultaneously. This method accurately predicts local stress fields in stress concentration regions and is computationally efficient as compared with the solution of a full scale microstructural model. This scheme has been applied to obtain localized damage at high and low stress zones of a V-notched rail shear specimen. The prominent damage mechanisms under shear loading, namely, matrix cracking and interfacial debonding, have been modeled using Mohr–Coulomb plasticity and traction separation law, respectively. The average stress at the notch has been found to be 44% higher than the average stresses away from the notch for a 90 N shear load. This stress rise is a direct outcome of the geometry of the notch.

Crush responses of composite cylinder under quasi-static and dynamic loading

Despite the abundance of studies investigating the performance of composite structures under crush loading, disagreement remains in the literature regarding the effect of increased strain rate on the crush response. This study reports an experimental investigation of the behaviour of a carbon-epoxy composite energy absorber under static and dynamic loading with a strain rate of up to 100s^-1. Consistent damage modes and measured force responses were obtained in samples tested under the same strain rate. The energy absorption was found to be independent of strain rate as the total energy absorption appeared to be largely associated with fibre-dominated fracture, which is independent of strain rate within the studied range. The results from this study are beneficial for the design of energy absorbing structures.

Predicting the crushing behaviour of composite material using high-fidelity finite element modelling

The capability to numerically model the crushing behaviour of composite structures will enable the efficient design of structures with high specific energy absorption capacity. This is particularly relevant to the aerospace and automotive industries where cabin structures need to be shown to be crashworthy. In this paper, a three-dimensional damage model is presented, which accurately represents the behaviour of composite laminates under crush loading. Both intralaminar and interlaminar failure mechanisms are taken into account. The crush damage model was implemented in ABAQUS/Explicit as a VUMAT subroutine. Numerical predictions are shown to agree well with experimental results, accurately capturing the intralaminar and interlaminar damage for a range of stacking sequences, triggers and composite materials. The use of measured material parameters required by the numerical models, without the need to ‘calibrate’ this input data, demonstrates this computational tool's predictive capabilities

Modelling Damage in Thin-Walled Structures [Keynote speaker]

A high-fidelity composite damage model is presented and applied to predict low-velocity impact damage, compression after impact (CAI) strength and crushing of thin-walled composite structures. The simulated results correlated well with experimental testing in terms of overall force-displacement response, damage morphologies and energy dissipation. The predictive power of this model makes it suitable for use as part of a virtual testing methodology, reducing the reliance on physical testing.  

Element-Free Galerkin Modelling of Crack Migration in Composite Laminates

The development of a virtual testing environment, as a cost-effective industrial design tool in the design and analysis of composite structures, requires the need to create models efficiently, as well as accelerate the analysis by reducing the number of degrees of freedom, while still satisfying the need for accurately tracking the evolution of a debond, delamination or crack front. The eventual aim is to simulate both damage initiation and propagation in components with realistic geometrical features, where crack propagation paths are not trivial. Meshless approaches, and the Element-Free Galerkin (EFG) method, are particularly suitable for problems involving changes in topology and have been successfully applied to simulate damage in homogeneous materials and concrete. In this work, the method is utilized to model initiation and mixed-mode propagation of cracks in composite laminates, and to simulate experimentally-observed crack migration which is difficult to model using standard finite element analysis. N

Design of laminated structures using piezoelectric materials

Composite structures incorporating piezoelectric sensors and actuators are increasingly becoming important due to the offer of potential benefits in a wide range of engineering applications such as vibration and noise supression, shape control and precisition positioning. This paper presents a finit element formulation based on classical laminated plate theory for laminated structures with integrated piezoelectric layers or patches, acting as actuators. The finite element model is a single layer triangular nonconforming plate/shell element with 18 degrees of freedom for the generalized displacements, and one electrical potential degree of freedom for each piezsoelectric elementlayer or patch, witch are surface bonded on the laminate. An optimization of the patches position is performed to maximize the piezoelectric actuators efficiency as well as, the electric potential distribuition is search to reach the specified structure transverse displacement distribuition (shape control). A gradient based algorithm is used for this purpose. The model is applied in the optimization of illustrative laminated plate cases, and the results are presented and discussed.