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.

Efficient modelling and optimisation of hybrid multilayered plates subject to ballistic impact

Among the key challenges present in the modelling and optimisation of composite structures against impact is the computational expense involved in setting up accurate simulations of the impact event and then performing the iterations required to optimise the designs. It is of more interest to find good designs given the limitations of the resources and time available rather than the best possible design. In this paper, low cost but sufficiently accurate finite element (FE) models were generated in LS Dyna for several experimentally characterised materials by semi-automating the modelling process and using existing material models. These models were then used by an optimisation algorithm to generate new hybrid offspring, leading to minimum weight and/or cost designs from a selection of isotropic metals, polymers and orthotropic fibre-reinforced laminates that countered a specified impact threat. Experimental validation of the optimal designs thus identified was then successfully carried out using a single stage gas gun. With sufficient computational hardware, the techniques developed in this pilot study can further utilise fine meshes, equations of state and sophisticated material models, so that optimal hybrid systems can be identified from a wide range of materials, designs and threats.

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 investigation of Mode I and Mode II fracture toughness enhancement using aligned carbon nanotubes forests at the crack interface

A novel approach for introducing aligned multi-walled carbon nanotubes (MWCNTs) in a carbon-fibre composite pre-impregnated (prepreg) laminate, to improve the through-thickness fracture toughness, is presented. Carbon nanotube (CNT) 'forests' were grown on a silicon substrate with a thermal oxide layer, using a chemical vapour deposition (CVD) process. The forests were then transferred to a pre-cured laminate interface, using a combination of pressure and heat, while maintaining through-thickness CNT alignment. Standard Mode I and four-point bend end-notched flexure Mode II tests were undertaken on a set of specimens and compared with pristine specimens. Mode I fracture toughness for T700/M21 laminates was improved by an average of 31% while for T700/SE84LV specimens, an improvement of 61% was observed. Only T700/M21 specimens were tested in Mode II which yielded an average fracture toughness improvement of 161%. Scanning Electron Microscopy (SEM) showed good wetting of the CNT forest as well as evidence of penetration of the forest into the adjacent plies. © 2013 Elsevier Ltd.

Experimental investigation of thermoforming carbon fibre-reinforced polyphenylene sulphide composites

Pre-consolidated carbon fibre-reinforced polyphenylene sulphide (CF/PPS) laminates were
thermoformed into V-shaped parts via designed out of autoclave thermoforming experiments.
The different processing conditions tested in the experiment have resulted in final
part angles whose differences ranged from 2.087 to 3.431 from the original mould angle.
The test results show that processing conditions influenced finished part dimensions as the
final sample angles were found to decrease relative to the tooling dimensions, as mould
temperature increases. Higher mould temperature conditions produce thinner parts due
to the thermal expansion of mould tools. The mould temperature of 170C, which can
produce parts with high degree of crystallinity as well as small size of crystal, has been
established as the optimal thermoforming condition for CF/PPS composites.

Numerical investigation of the effects of drill geometry on drilling induced delamination of carbon fiber reinforced composites

Drilling is a major process in the manufacturing of holes required for the assemblies of composite laminates in aerospace industry. Simulation of drilling process is an effective method in optimizing the drill geometry and process parameters in order to improve hole quality and to reduce the drill wear. In this research we have developed three-dimensional (3D) FE model for drilling CFRP. A 3D progressive intra-laminar failure model based on the Hashin's theory is considered. Also an inter-laminar delamination model which includes the onset and growth of delamination by using cohesive contact zone is developed. The developed model with inclusion of the improved delamination model and real drill geometry is used to make comparison between the step drill of different stage ratio and twist drill. Thrust force, torque and work piece stress distributions are estimated to decrease by the use of step drill with high stage ratio. The model indicates that delamination and other workpiece defects could be controlled by selection of suitable step drill geometry. Hence the 3D model could be used as a design tool for drill geometry for minimization of delamination in CFRP drilling. © 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

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

Improving composite damage modelling through automatic placement of cohesive elements

Numerous experimental studies of damage in composite laminates have shown that intralaminar (in-plane) matrix cracks lead to interlaminar delamination (out-of-plane) at ply interfaces. The smearing of in-plane cracks over a volume, as a consequence of the use of continuum damage mechanics, does not always effectively capture the full extent of the interaction between the two failure mechanisms. A more accurate representation is obtained by adopting a discrete crack approach via the use of cohesive elements, for both in-plane and out-of-plane damage. The difficulty with cohesive elements is that their location must be determined a priori in order to generate the model; while ideally the position of the crack migration, and more generally the propagation path, should be obtained as part of the problem’s solution. With the aim of enhancing current modelling capabilities with truly predictive capabilities, a concept of automatic insertion of interface elements is utilized. The consideration of a simple traction criterion in relation to material strength, evaluated at each node of the model (or of the regions of the model where it is estimated cracks might form), allows for the determination of initial crack location and subsequent propagation by the insertion of cohesive elements during the course of the analysis. Several experimental results are modelled using the commercial package ABAQUS/Standard with an automatic insertion subroutine developed in this work, and the results are presented to demonstrate the capabilities of this technique.

Experimental and numerical study of debonding in composite adhesive joints

The mode I and mode II fracture properties of the FM300-2 adhesive bond between 5HS/RTM6 laminates are determined experimentally by DCB and ELS test. The crack propagation is studied numerically by means of interface elements based on the decohesive zone model. The latter is characterized by material degradation, which is usually assumed to be linear. In the present study it is shown that if a non-linear material degradation is used with an increased magnitude of the interface relative displacement at failure it is possible to model more correctly the experimentally observed significant non-linear behaviour before the start of crack propagation. An adhesive stepped flush joint is studied experimentally and numerically. A mixed mode interaction criterion is used together with the nonlinear material degradation of the interface. Sensitivity studies are performed to study the influence of the parameters defining it.

Measurement of fibre fracture toughness using an alternative specimen geometry

This paper presents an experimental and numerical study focused on the tensile fibre fracture toughness characterisation of hybrid plain weave composite laminates using non-standardized Overheight Compact Tension (OCT) specimens. The position as well as the strain field ahead of the crack tip in the specimens was determined using a digital speckle photogrammetry system. The limitation on the applicability of standard data reduction schemes for the determination of the intralaminar fibre fracture toughness of composites 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 stress intensity factor for alternative composite specimen geometries. A comparison between different methods currently available to compute the intralaminar fracture toughness in composites is also presented and discussed. Good agreement between numerical and experimental results using the proposed methodology was obtained.

Modelling the nonlinear behaviour and fracture process of AS4/PEKK thermoplastic composite under shear loading

The accurate determination of non-linear shear behaviour and fracture toughness of continuous carbon-fibre/polymer composites remains a considerable challenge. These measurements are often necessary to generate material parameters for advanced computational damage models. In particular, there is a dearth of detailed shear fracture toughness characterisation for thermoplastic composites which are increasingly generating renewed interest within the aerospace and automotive sectors. In this work, carbon fibre (AS4)/ thermoplastic Polyetherketoneketone (PEKK) composite V-notched cross-ply specimens were manufactured to investigate their non-linear response under pure shear loading. Both monotonic and cyclic loading were applied to study the shear modulus degradation and progressive failure. For the first time in the reported literature, we use the essential work of fracture approach to measure the shear fracture toughness of continuous fibre reinforced composite laminates. Excellent geometric similarity in the load-displacement curves was observed for ligament-scaled specimens. The laminate fracture toughness was determined by linear regression, of the specific work of fracture values, to zero ligament thickness, and verified with computational models. The matrix intralaminar fracture toughness (ply level fracture toughness), associated with shear loading was determined by the area method. This paper also details the numerical implementation of a new three-dimensional phenomenological model for carbon fibre thermoplastic composites using the measured values, which is able to accurately represent the full non-linear mechanical response and fracture process. The constitutive model includes a new non-linear shear profile, shear modulus degradation and load reversal. It is combined with a smeared crack model for representing ply-level damage initiation and propagation. The model is shown to accurately predict the constitutive response in terms of permanent plastic strain, degraded modulus as well as load reversal. Predictions are also shown to compare favourably with the evolution of damage leading to final fracture.

Painéis fenólicos para aplicação em fachadas exteriores

Dissertação para obtenção do grau de Mestre em Engenharia Civil

Uncertainty propagation in inverse reliability-based design of composite structures

An approach for the analysis of uncertainty propagation in reliability-based design optimization of composite laminate structures is presented. Using the Uniform Design Method (UDM), a set of design points is generated over a domain centered on the mean reference values of the random variables. A methodology based on inverse optimal design of composite structures to achieve a specified reliability level is proposed, and the corresponding maximum load is outlined as a function of ply angle. Using the generated UDM design points as input/output patterns, an Artificial Neural Network (ANN) is developed based on an evolutionary learning process. Then, a Monte Carlo simulation using ANN development is performed to simulate the behavior of the critical Tsai number, structural reliability index, and their relative sensitivities as a function of the ply angle of laminates. The results are generated for uniformly distributed random variables on a domain centered on mean values. The statistical analysis of the results enables the study of the variability of the reliability index and its sensitivity relative to the ply angle. Numerical examples showing the utility of the approach for robust design of angle-ply laminates are presented.

Buckling strength of adhesively-bonded single and double-strap repairs on carbon-epoxy structures

This work reports on an experimental and finite element method (FEM) parametric study of adhesively-bonded single and double-strap repairs on carbon-epoxy structures under buckling unrestrained compression. The influence of the overlap length and patch thickness was evaluated. This loading gains a particular significance from the additional characteristic mechanisms of structures under compression, such as fibres microbuckling, for buckling restrained structures, or global buckling of the assembly, if no transverse restriction exists. The FEM analysis is based on the use of cohesive elements including mixed-mode criteria to simulate a cohesive fracture of the adhesive layer. Trapezoidal laws in pure modes I and II were used to account for the ductility of most structural adhesives. These laws were estimated for the adhesive used from double cantilever beam (DCB) and end-notched flexure (ENF) tests, respectively, using an inverse technique. The pure mode III cohesive law was equalled to the pure mode II one. Compression failure in the laminates was predicted using a stress-based criterion. The accurate FEM predictions open a good prospect for the reduction of the extensive experimentation in the design of carbon-epoxy repairs. Design principles were also established for these repairs under buckling.