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.     

The influence of ply sequence and thermoelastic stress field on asymmetric delamination crack growth behavior of embedded elliptical delaminations in laminated FRP composites

The influence of ply lay up and the interaction of residual thermal stresses and mechanical loading on the interlaminar asymmetric embedded delamination crack growth behavior have been investigated. Two sets of full three-dimensional thermo-elastic finite element analyses have been performed for the interlaminar elliptical delaminations, which may be due to manufacturing defects or other reasons and are located symmetrically with respect to the midplane in a quasi-isotropic FRP composite laminate lay up. Depending upon the through-the-thickness location of the embedded elliptical delaminations, four different laminate configurations have been considered. Strain energy release rate (SERR) procedures have been employed to assess the delamination crack growth characteristics at the interfaces. It is found that the individual fracture modes exhibit asymmetric and non self-similar crack growth behavior along the delamination front depending upon the location of the interfacial delaminations; ply sequence and orientation and thermo-elastic anisotropy of the laminae.     

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.     

Enriched finite element analysis for a delamination crack in a laminated composite strip

Edge delamination cracks in laminated composite strips are analyzed with the aid of the enriched finite element method, wherein the asymptotic singular solution for a delamination crack is incorporated into finite elements. The strip is assumed to be in the state of generalized plane deformations including extension (compression), bending or torsion. Comparison of the numerical results with those from other methods is made to confirm the solution. The crack growth stability is examined for a couple of ply orientations in terms of the energy release rate and mode mixity.This work has been partially supported by the Agency for the Defense Development in Republic of Korea: under the Grant No. ADD-92-5-004.     

A Discrete Model for Simulation of Composites Plate Impact Including Coupled Intra- and Inter-ply Failure

Laminated composites can undergo complex damage mechanisms when subjected to transverse impact. For unidirectional laminates it is well recognized that delamination failure usually initiates via intra-ply shear cracks that run parallel to the fibres. These cracks extend to the interface of adjacent orthogonal plies, where they are either stopped, or propagate further as inter-ply delamination cracks. These mechanisms largely determine impact energy absorption and post-delamination bending stiffness of the laminate. Important load transfer mechanisms will occur that may lead to fibre failure and ultimate rupture of the laminate. In recent years most Finite Element (FE) models to predict delamination usually stack layers of ply elements with interface elements to represent inter-ply stiffness and treat possible delamination. The approach is computationally efficient and does give some estimate of delamination zones and damaged laminate bending stiffness. However, these models do not properly account for coupled intra-ply shear failure and delamination crack growth, and therefore cannot provide accurate results on crack initiation and propagation. An alternative discrete meso-scale FE model is presented that accounts for this coupling, which is validated against common delamination tests and impact delamination from the Compression After Impact (CAI) test. Ongoing research is using damage prediction from the CAI simulation as a basis for residual strength analysis, which will be the published in future work.     

Numerical investigation to prevent crack jumping in Double Cantilever Beam tests of multidirectional composite laminates

During the experimental characterization of the mode I interlaminar fracture toughness of multidirectional composite laminates, the crack tends to migrate from the propagation plane (crack jumping) or to grow asymmetrically, invalidating the tests.The aim of this study is to check the feasibility of defining the stacking sequence of Double Cantilever Beam (DCB) specimens so that these undesired effects do not occur, leading to meaningful onset and propagation data from the tests. Accordingly, a finite element model using cohesive elements for interlaminar delamination and an intralaminar ply failure criterion are exploited here to thoroughly investigate the effect of specimen stiffness and thermal residual stresses on crack jumping and asymmetric crack growth occurring in multidirectional DCB specimens.The results show that the higher the arm bending stiffness, the lower the tendency to crack jumping and the better the crack front symmetry. This analysis raises the prospect of defining a test campaign leading to meaningful fracture toughness results (onset and propagation data) in multidirectional laminates.     

Simulation of delamination–migration and core crushing in a CFRP sandwich structure

Following the onset of damage caused by an impact load on a composite laminate structure, delaminations often form propagating outwards from the point of impact and in some cases can migrate via matrix cracks between plies as they grow. The goal of the present study is to develop an accurate finite element modeling technique for simulation of the delamination–migration phenomena in laminate impact damage processes. An experiment was devised where, under a quasi-static indentation load, an embedded delamination in the facesheet of a laminate sandwich specimen migrates via a transverse matrix crack and then continues to grow on a new ply interface. Using data from this test for validation purposes, several finite element damage simulation methods were investigated. Comparing the experimental results with those of the different models reveals certain modeling features that are important to include in a numerical simulation of delamination–migration and some that may be neglected.     

Energy-based prediction of failure in general symmetric laminates

This paper describes a new homogenisation technique for general symmetric laminates that enables progressive ply crack formation to be predicted in any number of plies having a variety of orientations. The approach involves (i) the analysis of non-uniformly spaced discrete ply cracks having a single orientation, (ii) a novel technique to homogenise the properties of the cracked ply so that discrete ply cracking can be analysed in plies having a different orientation, and (iii) the use of energy based methods to predict the progressive formation of ply cracks in any number of plies during loading. The analysis takes full account of the effects of thermally induced residual stresses. A key feature of the approach is the inclusion of a shear coupling term in the stress-strain relations for homogenised plies that ensures that the homogenised laminate has exactly the same effective properties as the laminate having a ply with discrete cracks in place of one of the homogenised plies. The model is applied to the prediction of the significant dependence of ply thickness and ply lay-up on laminate strength, and results for carbon fibre reinforced plastic laminates are compared favourably with published experimental results.     

Model for prediction of simultaneous time-dependent damage evolution in multiple plies of multidirectional polymer composite laminates and its influence on creep

A model to predict time-dependent evolution of simultaneous transverse cracking developed in multiple plies during creep loading and its effects on creep of multidirectional polymer matrix composite laminates is presented. The stress states in the intact regions of the plies are determined using the lamination theory during an incremental change in time. The stored elastic energy, determined using this stress state, is compared with a critical stored elastic energy value for damage to determine if a ply would fracture after the increment. If fracture is predicted, variational analysis is used to determine the perturbation in ply stresses due to cracking. This procedure is repeated to determine the crack evolution and creep strain. Model predictions compared well with experimental results for a [±θ_m/90_n]_s laminate.     

Fracture surface evolution and apparent delamination toughness in split composite beam specimens subjected to mixed mode I–III loading

The influence of specimen twisting during global anti-plane shear loading in composite split beam specimens is studied. Tests were conducted on specimens with different thicknesses and delamination lengths to produce different amounts of specimen twisting prior to fracture. It is shown that specimen twisting causes mode I stresses to develop, thereby producing mixed mode I–III conditions along the delamination front. This causes near-tip transverse cracks to initiate, prior to delamination advance, at an orientation related to the mode mix. Unlike in homogeneous materials, transverse crack extension is accompanied by planar delamination advance, and transverse crack rotation during extension is restricted by the laminate’s fibers. The overall fracture surface evolution is therefore strongly controlled by specimen geometry. The influence of these findings on the apparent delamination toughness as obtained from composite split beam and other types of mode III tests is discussed.     

Misinterpreting the results: How similitude can improve our understanding of fatigue delamination growth

Use of the strain energy release rate for characterizing delamination growth in composite and bonded structures is now commonplace. Analogous to the use of the stress intensity factor range for fatigue crack growth in metals, the strain energy release rate range, defined as the arithmetic difference between maximum and minimum values, is often used to characterize fatigue delamination growth behaviour. The basis for similitude for these two parameters, however, is different and can lead to misinterpretations of delamination growth behaviour if this difference is not understood. This paper examines the basis of similitude for the strain energy release rate range and how it can be redefined in order to avoid potential misinterpretations of fatigue delamination growth.     

Effect of stitch density and stitch thread thickness on low-velocity impact damage of stitched composites

Effect of stitch density and stitch thread thickness on low-velocity impact damage of stitched composites is investigated experimentally. Physical examination on damage surfaces shows that stitches perform as crack initiators, as well as crack arrestors. Longer matrix cracks are observed in densely-stitched composites, in contrast to isolated matrix cracks found in moderately-stitched composites. Ultrasonic C-scan evidently compares the delamination areas and concludes that specimens with higher stitch density and thread thickness are more capable of impeding delamination growth by effectively bridging delamination cracks. Load–time curves reveal that the onset of delamination is not influenced by stitch density and thread thickness. Energy consumption for the impact event is evaluated and discussed with the conclusion that, although absorbed energy is independent of stitch density and thread thickness, the proportion of energy consumption for damage mechanisms like delamination, matrix cracks and stitch debonding are different for laminated composites stitched with different stitch parameters.     

Crack deflection analyses of different peel stopper designs for sandwich structures

This paper relates to a newly developed peel stopper concept for sandwich structures. The proposed concept is a specially designed core insert, which has the ability to confine face sheet debonding/delamination (peeling) by deflecting a delamination crack front away from the face/core interface into the bulk of the sandwich core, and thereby constraining the debonding/delamination to a limited prescribed area. In this paper various peel stopper designs are analysed for their ability to deflect cracks away from propagating along a face–core interface. The crack deflection ability of the studied peel stopper designs leads to design guidelines, which describes the minimum requirements regarding the relation between the two interface toughnesses. The analysis further reveals that compliant peel stopper wedges are preferred because they lead to the lowest interface toughness ratio requirement. This has been confirmed through an experiment with a sandwich beam subjected to three-point bending loading. The experiment has shown that the ability of a peel stopper to deflect cracks is highly dependent on the stiffness of the wedge.     

A study on the influence of cohesive zone interface element strength parameters on mixed mode behaviour

This paper presents a detailed study of the influence of maximum interfacial stress on interface element analyses for composites delamination. The development of the non-linear cohesive zone ahead of a crack tip is analysed with respect to length, stress distribution and mode ratio. The energy absorbed by interface elements is compared with the crack tip strain energy release rate from fracture mechanics analyses. These studies are performed initially on standard fracture toughness specimens, where mode-ratio is fixed by the applied displacement constraints. Results show close agreement with linear elastic fracture mechanics solutions. A simple ply drop specimen is then modelled, where the mode ratio is not constrained by the boundary conditions, and results are compared with the Virtual Crack Closure Technique. In this case maximum interfacial stress has a far greater influence on the numerical results, due to its significant influence on cohesive zone length, mode ratio and energy absorbed.     

Designation of Mode Mix in Orthotropic Composite Delamination Problems

In interfacial fracture modeling of composite delamination, mode mix is typically specified in terms of energy release rates. Other near-tip quantities can be used to designate mode mix, however. This paper considers the designation of mode mix in terms of energy release rates, stress intensity factors, stresses ahead of the crack tip and crack face displacements and the consequences of using different near-crack-tip quantities to designate mode mix in analyzing composite delamination. The problem addressed is two-dimensional debonding between plies or ply groups modeled as in-plane orthotropic materials; however, the conclusions discussed apply to general composite delamination problems. It is shown that use of different quantities to designate mode mix can give significantly different results in matching composite applications to mixed-mode toughness tests. For cases where measured interfacial toughness increases with increasing mode II deformation, it is demonstrated that use of a mode mix designation based on energy release rates could be non-conservative. Based on these findings, it is suggested that practitioners consider the differences in failure load predictions that would result if different near-tip quantities were used to relate composite applications to measured toughnesses. To this end, methods for converting mode mix designations in terms of energy release rates into designations in terms of other fracture quantities are outlined and applied. This revised version was published online in August 2006 with corrections to the Cover Date.     

Matrix cracking and delamination in laminated composites. Part II: Evolution of crack density and delamination

This paper presents a model to predict the propagation of transverse cracks in polymer matrix composite laminates. Different possibilities for the crack pattern are analyzed and the different stress-strain response are compared. Taking into account that matrix cracking promotes delamination between the plies, the propagation of delamination is also simulated. The model predictions are compared with experimental data obtained in composite laminates that accumulate transverse cracks and delaminations before failing catastrophically. The possibility and limitations of a general constitutive law applied at ply level, as a mesomodel, is analyzed and the bounds of applicability of the model are explained.     

Embedded flaws for crack path control in composite laminates

An experimental investigation was conducted on using small flaws purposefully introduced into composite laminates to control growth of interlaminar cracks and through-thickness crack branching. Mode I crack growth specimens were used to study branching through 0°, 90° and 45° plies. The results showed that crack growth through 0° plies could be promoted by a ply gap, but this was not as controllable as combining a ply gap with a pre-crack to create a “crack branch flaw”. Crack branching through 45° plies could be controlled using crack branch flaws, and also promoted controllably using ply gaps. Crack branching through 90° plies was seen without any flaws, but was better controlled with embedded delaminations. Using these outcomes, crack branching through two quasi-isotropic laminates was demonstrated. The results have application to improved damage tolerance and fracture toughness, by taking advantage of high toughness crack growth mechanisms.     

铁素体管线钢的分层裂纹及其对断裂的影响

通过对针状铁素体管线钢缺口根部三维应力状态的有限元分析和不同形式的断裂实验,研究了管线钢分层裂纹产生的条件及其对断裂性能的影响.结果表明裂纹或缺口根部的三维应力状态是产生分层裂纹的必要条件,材料的强度分布影响分层裂纹的形式和方向.分层裂纹均为主裂纹扩展前材料中的弱界面在垂直该弱界面的拉应力作用下产生的,其数量和方向受裂纹端部三维应力场和材料的强度分布状态控制.分层裂纹面上的应力为零,分层裂纹有一定的间距.在断裂过程中产生的分层裂纹使裂纹或缺口根部的构形发生改变,从而对裂尖的应力状态和材料的断裂性能产生巨大的影响.穿透裂纹体的分层裂纹使其有效厚度减小,表面裂纹体的分层裂纹与裂纹扩展方向垂直.在断裂过程中产生分层裂纹需要消耗更多的能量、降低裂端三维应力约束、有效厚度降低或裂尖钝化.这些因素均使断裂扩展更加困难,而使材料韧性得到提高.     

Characterization and evolution of matrix and interface related damage in [0/90]S laminates under tension. Part I: Numerical predictions

This paper considers damage development mechanisms in cross-ply laminates using an accurate numerical method that assumes a Generalized Plane Strain (GPS) state. A 2D Boundary Element Method (BEM) model is generated to investigate the two types of damage progression in a [0/90]_S laminate: transverse cracks in the 90° lamina and delamination between both laminae. The model permits the contact between the surfaces of the cracks. The study is carried out in terms of the dependence of the Energy Release Rates (ERR) of the two types of crack on their respective lengths. A special emphasis is put on the mechanisms of the joining of the two aforementioned types of crack, including the study of the distribution of the stresses along the interface between the two plies when the transverse crack is approaching this interface.