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.     

Mesh‐independent matrix cracking and delamination modeling in laminated composites

The initiation and evolution of transverse matrix cracks and delaminations are predicted within a mesh‐independent cracking (MIC) framework. MIC is a regularized extended finite element method (x‐FEM) that allows the insertion of cracks in directions that are independent of the mesh orientation. The Heaviside step function that is typically used to introduce a displacement discontinuity across a crack surface is replaced by a continuous function approximated by using the original displacement shape functions. Such regularization allows the preservation of the Gaussian integration schema regardless of the enrichment required to model cracking in an arbitrary direction. The interaction between plies is anchored on the integration point distribution, which remains constant through the entire simulation. Initiation and propagation of delaminations between plies as well as intra‐ply MIC opening is implemented by using a mixed‐mode cohesive formulation in a fully three‐dimensional model that includes residual thermal stresses. The validity of the proposed methodology was tested against a variety of problems ranging from simple evolution of delamination from existing transverse cracks to strength predictions of complex laminates withouttextita priori knowledge of damage location or initiation. Good agreement with conventional numerical solutions and/or experimental data was observed in all the problems considered. Published 2011. This article is a US Government work and is in the public domain in the USA.     

Matrix crack-induced delamination in composite laminates under transverse loading

Fibre-reinforced multidirectional composite laminates are observed in experiments under transverse static or low-velocity impact loading to suffer considerable delamination damage. The intensity of this damage depends on the difference in the ply angles above and below the interface. In this paper a fracture mechanics model is presented for investigating the role of matrix cracks in triggering delaminations and the influence of ply angles in adjacent plies on delamination cracking. The fracture mechanics analysis shows that for a graphite fibre-reinforced composite laminate containing a transverse intraply crack, the crack-induced largest interfacial principal tensile stress is a maximum when the difference between the ply angles across the interface is 90 °, and it attains a minimum when the difference is 40 °. When the crack tips touch the interfaces, the minimum mode II stress singularity, which is weaker than the usual square-root type, appears when the difference between the ply angles is about 45 ° for one glass fibre-reinforced laminate and three graphite fibre-reinforced laminates. These results are in agreement with the experimental observation that the largest delaminations appear at the interface across which the difference between the ply angles is the largest i.e. 90 °.     

A study of the propagation of an embedded crack in a composite laminate of finite thickness

In this paper, the mechanism for the propagation of a crack embedded in the central layer of a composite laminate is studied based on energy criteria. As the crack may propagate in the transverse direction or in the tunneling direction, the critical stresses for these two propagations are presented. The theoretical analysis is applied to a fibre-reinforced composite laminate. Critical stresses for the initiation of crack propagation are calculated taking into account the influence of the adjacent ply angle and the crack size. It is found that for the considered typical fibre-reinforced composite laminate, the critical stress for the transverse propagation is less than that for the tunneling propagation. Therefore, transverse cracking is a more possible fracture mode in composite laminates with initial crack-like defects. This finding supports the calculation of the in situ transverse strength based upon the transverse crack propagation.     

Fracture toughness of transverse cracks in graphite/epoxy laminates at cryogenic conditions

Stress singularity of a transverse crack normal to ply-interface in a composite laminate is investigated using analytical and finite element methods. Four-point bending tests were performed on single-notch bend specimens of graphite/epoxy laminates containing a transverse crack perpendicular to the ply-interface. The experimentally determined fracture loads were applied to the finite element model to estimate the fracture toughness. The procedures were repeated for specimens under cryogenic conditions. Although the fracture loads varied with specimen thickness, the critical stress intensity factor was constant for all the specimens indicating that the measured fracture toughness can be used to predict delamination initiation from transverse cracks. For a given crack length and laminate configuration, the fracture load at cryogenic temperature was significantly lower. The results indicate that fracture toughness does not change significantly at cryogenic temperatures, but the thermal stresses play a major role in fracture and initiation of delaminations from transverse cracks.     

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.     

Integrated XFEM-CE analysis of delamination migration in multi-directional composite laminates

Progressive damage and failure in composites are generally complex and involve multiple interacting failure modes. Depending on factors such as lay-up sequence, loading and specimen configurations, failure may be dominated by extensive matrix crack-delamination interactions, which are very difficult to model accurately. The present study further develops an integrated extended finite element method (XFEM) and cohesive element (CE) method for three-dimensional (3D) delamination migration in multi-directional composite laminates, and validates the results with experiment performed on a double-cantilever beam (DCB). The plies are modeled by using XFEM brick elements, while the interfaces are modeled using CEs. The interaction between matrix crack and delamination is achieved by enriching the nodes of cohesive element. The mechanisms of matrix fracture and delamination migration are explained and discussed. Matrix crack initiation and propagation can be predicted and delamination migration is also observed in the results. The algorithm provides for the prediction of matrix crack angles through the ply thickness. The proposed method provides a platform for the realistic simulation of progressive failure of composite laminates.     

Energy-based delamination theory for biaxial loading in the presence of thermal stresses

This paper addresses the issue of using energy balance methods and crack closure concepts to predict the growth of delaminations associated with ply cracks during the progressive loading of cross-ply laminates subject to a combination of in-plane biaxial stresses and thermal residual stresses. When the effective applied stresses and the temperature are held fixed during delamination growth, and there is negligible interaction of the delamination tips with the ply cracks, very simple analytical formulae for the energy release rate can be derived for unconstrained and generalised plane strain conditions, which are exact when the ply crack separation tends to infinity.     

Modelling multiple cohesive crack propagation using a finite element–scaled boundary finite element coupled method

This paper presents an extension of the recently-developed finite element–scaled boundary finite element (FEM–SBFEM) coupled method to model multiple crack propagation in concrete. The concrete bulk and fracture process zones are modelled using SBFEM and nonlinear cohesive interface finite elements (CIEs), respectively. The CIEs are automatically inserted into the SBFEM mesh as the cracks propagate. The algorithm previously devised for single crack propagation is augmented to model problems with multiple cracks and to allow cracks to initiate in an un-cracked SBFEM mesh. It also addresses crack propagation from one subdomain into another, as a result of partitioning a coarse SBFEM mesh, required for some mixed–mode problems. Each crack in the SBFEM mesh propagates when the sign of the Mode-I stress intensity factor at the crack tip turns positive from negative. Its propagation angle is determined using linear elastic fracture mechanics criteria. Three concrete beams involving multiple crack propagation are modelled. The predicted crack propagation patterns and load–displacement curves are in good agreement with data reported in literature.     

Delamination analysis of composites by new orthotropic bimaterial extended finite element method

The extended finite element method (XFEM) is further improved for fracture analysis of composite laminates containing interlaminar delaminations. New set of bimaterial orthotropic enrichment functions are developed and utilized in XFEM analysis of linear‐elastic fracture mechanics of layered composites. Interlaminar crack‐tip enrichment functions are derived from analytical asymptotic displacement fields around a traction‐free interfacial crack. Also, heaviside and weak discontinuity enrichment functions are utilized in modeling discontinuous fields across interface cracks and bimaterial weak discontinuities, respectively. In this procedure, elements containing a crack‐tip or strong/weak discontinuities are not required to conform to those geometries. In addition, the same mesh can be used to analyze different interlaminar cracks or delamination propagation. The domain interaction integral approach is also adopted in order to numerically evaluate the mixed‐mode stress intensity factors. A number of benchmark tests are simulated to assess the performance of the proposed approach and the results are compared with available reference results. Copyright © 2011 John Wiley & Sons, Ltd.     

Cohesive models for damage evolution in laminated composites

A trend in the last decade towards models in which nonlinear crack tip processes are represented explicitly, rather than being assigned to a point process at the crack tip (as in linear elastic fracture mechanics), is reviewed by a survey of the literature. A good compromise between computational efficiency and physical reality seems to be the cohesive zone formulation, which collapses the effect of the nonlinear crack process zone onto a surface of displacement discontinuity (generalized crack). Damage mechanisms that can be represented by cohesive models include delamination of plies, large splitting (shear) cracks within plies, multiple matrix cracking within plies, fiber rupture or microbuckling (kink band formation), friction acting between delaminated plies, process zones at crack tips representing crazing or other nonlinearity, and large scale bridging by through-thickness reinforcement or oblique crack-bridging fibers. The power of the technique is illustrated here for delamination and splitting cracks in laminates. A cohesive element is presented for simulating three-dimensional, mode-dependent process zones. An essential feature of the formulation is that the delamination crack shape can follow its natural evolution, according to the evolving mode conditions calculated within the simulation. But in numerical work, care must be taken that element sizes are defined consistently with the characteristic lengths of cohesive zones that are implied by the chosen cohesive laws. Qualitatively successful applications are reported to some practical problems in composite engineering, which cannot be adequately analyzed by conventional tools such as linear elastic fracture mechanics and the virtual crack closure technique. The simulations successfully reproduce experimentally measured crack shapes that have been reported in the literature over a decade ago, but have not been reproduced by prior models.     

Composite toughness enhancement with interlaminar reinforcement

An experimental investigation of a newly proposed through-thickness reinforcement approach aimed to increase interlaminar toughness of laminated composites is presented. The approach alters conventional methods of creating three-dimensional fiber-reinforced polymer composites in that the reinforcing element is embedded into the host laminate after it has been cured. The resulting composite is shown to possess the benefits of a uniform surface quality and consolidation of the original unreinforced laminate. This technique was found to be highly effective in suppressing the damage propagation in delamination double-cantilever beam (DCB) test samples under mode I loading conditions. Pullout testing of a single reinforcing element was carried out to understand the bridging mechanics responsible for the improved interlaminar strength of reinforced laminate and stabilization and/or arrest of delamination crack propagation. The mode I interlaminar fracture of reinforced DCB samples was modeled using two-dimensional cohesive finite-element scheme to support interpretation of the experiments.     

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.     

High-cycle fatigue characteristics of quasi-isotropic CFRP laminates over 10 cycles (Initiation and propagation of delamination considering interaction with transverse cracks)

High-cycle fatigue features of over 10⁸ cycles, particularly the initiation and propagation of edge delamination considering the effects of transverse cracks, were investigated using quasi-isotropic carbon-fiber-reinforced plastic (CFRP) laminates with a stacking sequence of [45/0/−45/90]_s in this study. In the relationship between a transverse crack density and initiation and growth of edge delamination, it was found that fatigue damage growth behavior varied depending on applied stress. It was observed that edge delamination initiated and grew at parts where transverse cracks were dense at ordinary applied stress, whereas it was observed that edge delamination grew before or simultaneously with transverse crack propagation at a low applied stress and high-cycle loading. In addition, the critical transverse crack density where delamination begins growing was calculated to evaluate the interaction between transverse crack and edge delamination growth.     

含面芯界面缺陷的蜂窝夹芯板侧向压缩破坏模式

为了对含面芯层间脱胶缺陷的蜂窝夹芯板在侧向压缩载荷作用下的典型破坏模式进行数值预报, 建立了基于蔡-希尔破坏准则和粘结模型的计算模型。该计算模型是建立在对蜂窝夹芯板的双悬臂梁(DCB)和单臂梁(SLB) 试验中所发现的一种新的破坏模式的分析基础之上的。对蜂窝夹芯板的侧向压缩破坏行为的数值预报中, 发现一种新的破坏模式: 位于脱胶区域的面板首先发生局部屈曲失稳, 随后面板内部靠近芯子的45°/0°层间出现分层, 与此同时最靠近芯子的45°铺层发生断裂, 伴随着45°/0°层间分层的扩展, 面板发展成为对称性整体屈曲失稳。与侧向压缩试验测试结果对比发现, 计算模型模拟中所预报的破坏模式在实验测试中也得到了很好的验证。      

Numerical analysis of size effects on open-hole tensile composite laminates

The tensile strength of open-hole fibre reinforced composite laminates depends on in-plane, thickness and ply lay-up scaling. Translaminar (fibre direction) mode I fracture toughness has recently been experimentally determined to be thickness dependent. This paper presents a computational study of the tensile strength prediction of open-hole laminates using a cohesive zone model. To the authors’ knowledge, it is for the first time in the literature that the thickness-dependence of translaminar fracture toughness is accounted for in the numerical modelling of composites. The thickness size effect in the strength of open-hole composite laminates failed by pull-out is accurately predicted for the first time by a deterministic model. It is found that neglecting delamination in the numerical models will lead to mesh-dependency and over-estimation on the predicted strength. Smeared crack model with cohesive elements to model delamination is able to predict the correct failure mode; but it is found not suitable for accurate strength predictions for laminates failed by delamination.     

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.     

Stiffness and fracture analysis of laminated composites with off-axis ply matrix cracking

Matrix cracks parallel to the fibres in the off-axis plies is the first intralaminar damage mode observed in laminated composites subjected to static or fatigue in-plane tensile loading. They reduce laminate stiffness and strength and trigger development of other damage modes, such as delaminations. This paper is concerned with theoretical modelling of unbalanced symmetric laminates with off-axis ply cracks. Closed-form analytical expressions are derived for Mode I, Mode II and the total strain energy release rates associated with off-axis ply cracking in [0/θ]_s laminates. Stiffness reduction due to matrix cracking is also predicted analytically using the Equivalent Constraint Model (ECM) of the damaged laminate. Dependence of the degraded stiffness properties and strain energy release rates on the crack density and ply orientation angle is examined for glass/epoxy and carbon/epoxy laminates. Suitability of a mixed mode fracture criterion to predict the cracking onset strain is also discussed.     

Thermal fatigue failure induced by delamination in thermal barrier coating

The paper presents the experimental and theoretical investigation on the thermal fatigue failure induced by delamination in thermal barrier coating system. Laser heating method was used to simulate the operating state of TBC (thermal barrier coating) system. The non-destructive evaluation such as acoustic emission (AE) detect was used to study the evolution of TBC system damage. Micro-observation and AE detect both revealed that fatigue crack was in two forms: surface crack and interface delamination. It was found that interface delamination took place in the period of cooling or heating. Heating or cooling rate and temperature gradient had an important effect on interface delamination cracking propagation. A theoretical model on interface delamination cracking in TBC system at operating state is proposed. In the model, a membrane stress P and a bending moment M are designated the thermal loads of the thermal stress and temperature gradient in TBC system. In this case, the coupled effect of plastic deformation, creep of ceramic coating as well as thermal growth oxidation (TGO) and temperature gradient in TBC system was considered in the model. The thermal stress intensity factors (TSIFs) in non-FGM (functional gradient material) thermal barrier coating system is analytical obtained. The numerical results of TSIFs reveal some same results as obtained in experimental test. The model is based on fracture mechanics theory about heterogeneous materials and it gives a rigorous explanation of delaminations in TBC system loaded by thermal fatigue. Both theoretical analysis and experimental observation reveal an important fact: delaminations are fatigue cracks which grow during thermal shocks due to compressive stresses in the loading, this loads the delaminations cracks in mixed I and II mode.