Finite element analysis of postbuckling and delamination of composite laminates using virtual crack closure technique

The two-dimensional and three-dimensional parametric finite element analysis (FEA) of composite flat laminates with two through-the-width delamination types: 0₄/(±θ)₆//0₄ and 0₄//(±θ)₆//0₄ (θ = 0°, 45°, and “//” denotes the delaminated interface) under compressive load are performed to explore the effects of multiple delaminations on the postbuckling properties. The virtual crack closure technique which is employed to calculate the energy release rate (ERR) for crack propagation is used to deal with the delamination growth. Three typical failure criteria: B-K law, Reeder law and Power law are comparatively studied for predicting the crack propagation. Effects of different mesh sizes and pre-existing crack length on the delamination growth and postbuckling properties of composite laminates are discussed. Interaction between the delamination growth mechanisms for multiple cracks for 0₄//(±θ)₆//0₄ composite laminates is also investigated. Numerical results using FEA are also compared with those by existing models and experiments.     

Effect of off-axis ply orientation on 0°-fibre microbuckling

The aim of the present work is to study both experimentally and theoretically the compression failure mechanisms in multi-directional composite laminates, and especially the effect of the off-axis ply orientation on fibre microbuckling in the 0°-plies. The critical mechanism in the compressive fracture of unidirectional polymer matrix composites is plastic microbuckling/kinking. In multi-directional composites with internal 0°-plies, catastrophic failure also initiates by kinking of 0°-plies at the free-edges or manufacturing defects, followed by delamination. When 0°-plies are located at the outside, or in the case of cross-ply laminates, failure rather tends to occur by out-of-plane buckling of the 0°-plies. T800/924C carbon-fibre–epoxy laminates with a [(±θ/0₂)₂]_s lay-up are used here to study the effect of the supporting ply angle θ on the stress initiation of 0°-fibre microbuckling. Experimental data on the compressive strength of laminates with θ equal to 30, 45, 60 or 75° are compared to theoretical predictions obtained from a fibre kinking model that incorporates interlaminar shear stresses developed at the free edges at (0/θ) interfaces. Initial misalignment of the fibres and non-linear shear behaviour of the matrix are also included in the analysis.     

Effects of boundary conditions on buckling and postbuckling responses of composite laminate with various shaped cutouts

The objective of the present work is to study the effects of flexural boundary conditions on buckling and postbuckling behavior of axially compressed quasi-isotropic laminate, (+45/−45/0/90)_2s with various shaped cutouts (i.e., circular, square, diamond, elliptical–vertical and elliptical–horizontal) of various sizes using the finite element method. The FEM formulation is based on first order shear deformation theory and von Karman’s assumptions are used to incorporate geometric nonlinearity. The 3-D Tsai-Hill criterion is used to predict the failure of a lamina while the onset of delamination is predicted by the interlaminar failure criterion. It is observed that the laminates clamped and simply supported on all edges have the highest and lowest buckling and postbuckling strength, respectively, irrespective of cutouts shape and size. It is also noted that a fully clamped laminate with large size elliptical–vertical cutout can take higher compressive buckling load than the laminate without cutout for same boundary condition.     

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.     

Post-impact compressive strength of curved GFRP laminates

Low velocity impacts to fibre reinforced plastic composites cause a pattern of damage consisting in general of delamination, fibre breakage and matrix cracking. Such damage is accidental and may go unnoticed; therefore composite structures must be designed assuming impact damage exists. Previous work on flat composite laminates has resulted in a reasonable understanding of the mechanisms of compressive strength reduction. There are, however, many instances where curved laminates are used in structures where impact is likely. Furthermore, due to the mechanisms of strength reduction, it may be expected that curvature would have a significant effect on the behaviour of the laminates.The work described here consists of experimental measurement of the post-impact compressive strength of curved GFRP laminates. The laminates were of 8 plies of 0.3 mm thick pre-impregnated glass fibre/epoxy tape in a (0, ±45, 0°)_s lay-up. Each laminate was 200 mm in length by 50 mm wide with the plane of curvature normal to the length. Laminates were impacted on the convex surface of the laminate by dropping a steel mass from 1 m vertically above it.Impacted laminates were loaded in compression and the out-of-plane displacements of the top and bottom surfaces were recorded. Final failure was typically due to fibre breakage occurring through the centre of the impacted area of the laminate. Possible differences in the impact response, and measurable differences in the sizes of the impact damage area, were found to arise from these curvatures, and differences were observed in their post-impact buckling behaviour. However, perhaps unexpectedly, the post-impact compressive strength for a curved laminate was found to be similar to that for a flat laminate. The failure loads for the impact damage laminates are shown to be comparable with those for laminates containing artificial delaminations.     

The stability of delamination growth in compressively loaded composite plates

The stability of growth of internal delaminations in composite plates subjected to compressive loading is investigated. Due to the compressive loading, these structures can undergo buckling of the delaminated layer and subsequently growth of the delamination. The study does not impose any restrictive assumptions regarding the delamination thickness and plate length (as opposed to the usual thin film assumptions). The growth characteristics of the delamination under monotonic compressive loading are obtained on the basis of a combined delamination buckling/postbuckling and fracture mechanics model. The postbuckling solution is derived through a perturbation procedure, which is based on an asymptotic expansion of the load and deformation quantities in terms of the distortion parameter of the delaminated layer, the latter being considered a compressive elastica. The closed form solutions for the energy release rate at the delamination tip versus applied compressive strain during the initial postbuckling phase are used to define the combinations of delamination length and applied strain that lead to unstable growth. This would practically cause either contained jump growth or complete (catastrophic) growth of the delamination. Estimates for the lower and upper bounds of the jump distance (unstable growth) are provided. Moreover, a study of the influence of the mode dependence of interface toughness on the conditions of initiation and extent of delamination growth is performed.     

Stability and failure of composite laminates with various shaped cutouts under combined in-plane loads

The objective of this paper is to study stability and failure of a composite laminate with a centrally placed cutout of various shapes (i.e., circular, square, diamond, elliptical-vertical and elliptical-horizontal) under combined action of uni-axial compression and in-plane shear loads. The FEM formulation based on the first order shear deformation theory and von Karman’s assumptions has been utilized. Newton–Raphson method is used to solve nonlinear algebraic equations. Failure of a lamina is predicted by the 3-D Tsai–Hill criterion whereas the onset of delamination is predicted by the interlaminar failure criterion. The effects of cutout shape, direction of shear load and composite lay-up on buckling and postbuckling responses, failure loads and failure characteristics of the laminate has been discussed. An efficient utilization of material strength is observed in the case of laminate with circular cutout as compared to the laminate with other shaped cutouts. In addition, it is also concluded that although the buckling strength of the (0/90)_4s laminate is lower than that of the (+45/?45/0/90)_2s and (45/?45)_4s laminates, but its strength is increased in the advanced stage of postbuckling deformation.     

Fatigue life estimation of graphite/epoxy laminates under consideration of delamination growth

A closed form approach to the assessment of the fatigue life of graphite/epoxy laminates under cyclic tension–compression loading has been developed. The model is mechanistic and uses cyclic energy release rates for prediction of delamination growth and of critical delamination sizes which induce buckling and the final failure of the laminates. Tests performed with graphite/epoxy specimens of stacking order [0_n, ?_m]_s with severed central plies [?], and of stacking order [0₂, +45, 0₂, ?45, 0, 90]_s with a central unloaded hole, indicate good correlation between estimated values and observed delamination growth, critical buckling strength of separated plies and load cycles to failure.     

Buckling and postbuckling behavior of delaminated sandwich beams

An investigation was performed to study the buckling and postbuckling behavior of sandwich beams containing lengthwise and depthwise through-the-width delaminations. An analytical beam model was developed to predict the buckling load of the beam and to describe its postbuckling response for arbitrarily situated delaminations and various combinations of boundary conditions. Griffith's energy release rate model was employed to predict the stability of delamination propagation under external loading and to determine the direction of delamination growth.

Parametric studies over a wide range of beam geometries, damage sizes and locations, composite facings and beam boundary conditions were carried out to study their effects on the overall behavior of the sandwich structure, as well as its damage tolerance. The results demonstrated that a sandwich construction is very ‘sensitive’ to the presence of delaminations situated at the core-faceplate interface. Premature buckling failure occurs at external loads which are significantly lower than the buckling load for a ‘perfect’ sandwich beam; in ‘imperfect’ beams with composite faceplates, the layup sequence affects significantly the load-carrying capacity of the beam; varying either the boundary conditions in a sandwich beam or the lengthwise location of a delamination has a small effect on the postbuckling behavior of the beam. Delaminations located within composite faceplates have less pronounced influence, and as the defect is moved outwards the limit load may reach the buckling load corresponding to that of the ‘perfect’ beam.

The proposed model is capable of analyzing the postbuckling behavior of both sandwich and composite laminated beams for arbitrary locations of the delamination, and various combinations of boundary conditions.     

Postbuckling strengths of composite laminate with various shaped cutouts under in-plane shear

The aim of present investigation is to study the buckling and postbuckling response and strengths under positive and negative in-plane shear loads of simply-supported composite laminate with various shaped cutouts (i.e., circular, square, diamond, elliptical-vertical and elliptical-horizontal) of various sizes using finite-element method. The FEM formulation is based on the first order shear deformation theory which incorporates geometric nonlinearity using von Karman’s assumptions. The 3-D Tsai-Hill criterion is used to predict the failure of a lamina while the onset of delamination is predicted by the interlaminar failure criterion. The effect of cutout shape, size and direction of shear load on buckling and postbuckling responses, failure loads and failure characteristics of quasi-isotropic [i.e., (+45/−45/0/90)_2s] laminate has been discussed. In addition, the effect of composite lay-up [i.e., (+45/−45/0/90)_2s, (45/−45)_4s and (0/90)_4s] has also been reported. It is observed that the cutout shape has considerable effect on the buckling and postbucking behaviour of the quasi-isotropic laminate with large size cutout. It is also observed that the direction of shear load and composite lay-up have substantial influence on strength and failure characteristics of the laminate.     

Study of critical buckling loads and modes of cross-ply laminated annular plates

Buckling of composite annular plates under uniform internal and external radial edge loads have been investigated using energy method. Trefftez rule is used in the stability equations. The symmetric buckling of symmetric cross-ply laminates is considered. In this paper, buckling behavior for the three laminates (90/0)_2s, (90/0₂/90)_s and (90₂/0₂)_s are studied. Influence of some parameters such as thickness, stacking sequence, type of supports and the ratio of hole to sheet radius on buckling loads and modes are investigated. The results of the energy method are compared with the results of numerical method. Based on the results, in the plates with clamped boundary conditions the symmetric buckling assumption is not accurate, contrary to other boundary conditions.     

Local buckling of stitched composite laminate

Due to relatively low interlaminar strength, delamination is a common failure mode of composite laminates. Through-thickness stitching is shown to improve the delamination resistance of laminated composites. Under in-plane compressive loading, significant strength reduction occurs due to coupling between delamination and delamination buckling. In this paper, an energy-based model was developed to predict the effect of critical stitching parameters on the delamination buckling strength of stitched laminates. Excellent agreement was found between the model results and a corresponding finite element analysis.     

Residual compressive strength of laminated plates with delamination

A study of residual compressive strength in delaminated laminates is presented. A methodology is proposed for simulating the whole compressive failure responses, such as initial buckling, postbuckling, contact of delamination front region, delamination propagation, fiber breakage, and matrix cracking etc. An finite element analysis (FEA) of the residual compressive strength is conducted on the basis of the Von Karman's nonlinearity assumption and the first-order shear deformation plate theory, combined with a stiffness degradation scheme. The numerical analysis models and methods are briefly introduced in this paper and some numerical examples are presented to illustrate it. From numerical results and discussion, it is clear that the compressive failure response involves complex multi-failure modes during compressive process. The method and numerical conclusions provide in this paper should of great value to engineers dealing with composite structures.     

Delamination buckling growth in laminated composites using layerwise-interface element

In this paper, a numerical investigation on the buckling of composite laminates containing delamination, under in-plane compressive loads, is presented. For this purpose, delamination propagation is modeled using the softening behavior of interface elements. The full layerwise plate theory is applied for approximating the displacement field of laminates and the interface elements are considered as a numerical layer between any two adjacent layers where the delamination is expected to propagate. A non-linear computer code was developed to handle the numerical procedure of delamination buckling growth in composite laminates using layerwise-interface elements. The load/displacement behavior and the contours of embedded and through-the-width delamination propagation for composite laminates are presented. It is shown that delamination growth can be well predicted using this layerwise-interface elements with decohesive law.     

含椭圆分层损伤复合材料层合板分层前缘接触效应研究

采用基于Mindlin 一阶剪切理论的四节点板单元, 分析了分层前缘闭合接触效应对含椭圆分层复合材料层合板前后屈曲行为的影响。通过在分层区域内引入虚杆单元来防止上、下子板发生相互重叠, 并假设接触区内没有摩擦。结果表明, 分层前缘的闭合接触效应受分层几何形状的影响很大, 当椭圆分层子板的长轴与压缩载荷方向一致时, 闭合接触效应对含分层层合板前后屈曲行为最大。      

On matrix crack saturation in composite laminates

This paper presents a comparative experimental investigation on matrix crack saturation in the constrained laminae of glass–fibre resin composite laminates. Fatigue tests under load and strain control were conducted on [±θ/90₃]_s and [0/±θ₄/0]_T laminates, respectively. The crack density was periodically measured by using digital image method and optical microscope examination during fatigue tests, to determine crack saturation. It is found that matrix crack saturation occurs in all laminates of different lay-ups under both stress and strain control cyclic loading. The crack density at saturation in the constrained plies is dependent on the cyclic loading type, the off-axis angle of the laminates, as well as the constraining ply. Based on the strain energy release rate, a crack saturation criterion is developed and the experimental findings are qualitatively analysed.     

低速冲击后复合材料层合板的压缩破坏行为

对缝纫层合板和无缝纫层合板进行低速冲击后压缩破坏实验,以研究低速冲击后层合板的压缩破坏机理。采用C扫描、X射线、热揭层等技术对层合板内的损伤进行测量和对比。结果表明,界面不是很强的碳纤维增强复合材料层合板低速冲击后受压时,层合板非冲击面的子层屈曲及其扩展是导致层合板冲击后压缩强度下降的重要因素,而且子层屈曲主要是沿垂直载荷的方向(90°)扩展;对于准各向同性板,屈曲子层中与母层相邻的铺层的方向一般为90°。层合板的剩余压缩强度与板的冲击损伤面积无直接关系。      

Improvement of interlaminar mechanical properties of CFRP laminates using VGCF

In order to improve the interlaminar mechanical properties of CFRP laminates, hybrid CFRP/VGCF laminates have been fabricated by using a newly-developed method, i.e., powder method, where the powder of vapor grown carbon fiber (VGCF) is added at the mid-plane of [0°/0°]₁₄ CFRP laminates. Experimental results of double cantilever beam (DCB) tests indicate the improvement on the interlaminar mechanical properties of Mode-I fracture behavior with much higher critical load P_C and fracture toughness G_IC with VGCF interlayer. Crack propagation and fracture surface have also been observed to interpret this improvement mechanism. Moreover, based on experimental G_IC, numerical simulations using finite element method (FEM) with cohesive elements have been carried out to analyze the delamination propagation. The interlaminar tensile strength of hybrid CFRP/VGCF laminates, which is obtained by matching the numerical load–COD (crack opening displacement) curves to experimental ones, is higher than that of base CFRP laminates.     

Interaction of laminate damage and adhesive disbonding in composite scarf joints subjected to combined in-plane loading and impact

Impact tests were carried out on composite laminates and composite scarf repairs, while both were subjected to in-plane loading with tensile pre-strain levels up to 5000 microstrain. The results show that pre-straining of the composite laminates has no noticeable influence on the size of the delamination area for the given impact energy of 8 J, which represents a typical barely-visible impact on thin-skin composite structures. For composite scarf joints, however, resulting damage has been found to be a combination of adhesive disbonding and matrix cracking (delamination and intraply cracking) in the composite laminate. The size of this mixed type of damage increases significantly with increasing pre-strain levels. A finite element model was developed to investigate the interaction between adhesive disbonding and composite delamination. The computational results reveal that both delamination and adhesive disbonding are dominated by the mode II fracture. Since the critical mode II fracture energy release rate for composite laminates (G_IIC = 1.08 kJ/m²) is much less than that pertinent to the adhesive (G_IIC = 3.73 kJ/m²), delamination tends to occur first in the composite laminates, which then shield the growth of disbonding in the adhesive.