Prediction of Mode I interlaminar fracture toughness of stitched flax fiber composites

This paper aims to propose a simulation procedure to predict the interlaminar fracture toughness of stitched flax fiber composites through a virtual double cantilever beam test. The proposed procedure is constituted of two steps. First, the interlaminar failure of unstitched flax fiber laminate, as the parent laminate, is modeled using cohesive elements with a nonlinear softening law in order to model the large-scale fiber bridging occurred during delamination. The experimental results are used to calibrate the parameters of the cohesive law. Second, two-node beam elements are superposed onto the cohesive interface of the parent laminate at a prescribed stitch density and distribution to model the bridging stitches present in the validation samples. The stitch material behavior and properties are obtained from the tensile test of impregnated stitch fibers. The out-of-plane flax yarn stitching was found to generate a twofold increase in the delamination resistance of the composite laminate at a medium stitch density. The FE analysis results agreed well with the experimental results, where a good fit between the predicted and experimental R-curves was achieved.     

Effect of stitch density and stitch thread thickness on compression after impact strength and response of stitched composites

In this paper, the damage failure and behaviour of stitched composites under compression after impact (CAI) loading are experimentally investigated. This study focuses on the effect of stitch density and stitch thread thickness on the CAI strength and response of laminated composites reinforced by through-thickness stitching. Experimental findings show that stitched composites have higher CAI failure load and displacement, which corresponds to higher energy absorption during CAI damage, mainly attributed to greater energy consumption by stitch fibre rupture. The coupling relationships between CAI strength, impact energy, stitch density and stitch thread thickness are also revealed. It is understood that the effectiveness of stitching has high dependency on the applied impact energy. At low impact energy range, CAI strength is found to be solely dependent on stitch density, showing no influence of stitch thread thickness. It is however observed that stitch fibre bridging is rendered ineffective in moderately stitched laminates during compressive failure, as local buckling occurs between stitch threads, resulting in unstitched and moderately stitched laminates have similar CAI strength. The CAI strength of densely stitched laminates is much higher due to effective stitch fibre bridging and numerous stitch thread breakages. At high impact energy level, CAI strength is discovered to be intimately related to both stitch density and stitch thread thickness. Since CAI failure initiates from impact-induced delamination area, stitch fibre bridging is considerable for all specimens due to the relatively large delamination area present. Stitch threads effectively bridge the delaminated area, inhibit local buckling and suppress delamination propagation, thus leading to increased CAI strength for laminates stitched with higher stitch density and larger stitch thread thickness. Fracture mechanisms and crack bridging phenomenon, elucidated by X-ray radiography are also presented and discussed. This study reveals novel understanding on the effectiveness of stitch parameters for improving impact tolerance of stitched composites.     

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 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.     

含分层损伤缝合复合材料层板的剩余压缩强度

基于渐进损伤方法,研究了含单脱层缝合复合材料层板在压缩载荷下的剩余强度。通过商用软件ABAQUS建立了含单脱层缝合复合材料层板剩余压缩强度计算模型,考虑了子层屈曲和分层扩展对剩余强度的影响。通过UMAT子程序实现了层板失效、层间失效和缝线失效的模拟。通过嵌入式杆单元结构模拟了缝线桥联作用及失效。采用Hashin准则及刚度折减法对纤维拉压、基体拉压失效进行了模拟。通过渐进损伤分析,揭示了缝合情况下含单脱层复合材料层板的失效机理,讨论了缝合参数对剩余压缩强度的影响。所预测的破坏模式和剩余强度结果与实验能较好地吻合。分析表明缝合可以明显提高含分层损伤复合材料层板的子层屈曲载荷,抑制分层扩展,并提高层板的剩余压缩强度。     

纤维增强复合材料层板高速冲击损伤数值模拟

推导了复合材料应变率相关三维本构关系, 并将其用于复合材料层板高速冲击损伤的数值模拟。该模型在复合材料层间引入界面单元模拟层间分层, 结合三维Hashin失效准则进行单层板面内损伤识别, 引入材料刚度退化, 采用非线性有限元方法, 研究了复合材料层板高速冲击的破坏过程及层板的损伤特性。数值分析结果表明, 剩余速度预报结果与实验结果吻合较好, 层板的主要损伤形式是层间分层、 基体微裂纹和纤维断裂, 减小弹体直径、 增大铺层角度和层板厚度能够有效降低层板损伤面积。      

缝合复合材料层板低速冲击损伤数值模拟

建立了缝合复合材料层板在低速冲击载荷下的渐进损伤分析模型。模型中采用空间杆单元模拟缝线的作用;采用三维实体单元模拟缝合层板,通过基于应变描述的Hashin准则,结合相应的材料性能退化方案模拟层板的损伤和演化;采用界面单元模拟层间界面,结合传统的应力失效判据和断裂力学中的应变能释放率准则判断分层的起始和扩展规律。通过对碳800环氧树脂复合材料（T800/5228）层板的数值仿真结果和试验结果相比较,验证了模型的正确性,同时讨论了不同冲击能量下缝合层板的损伤规律。研究结果表明：缝线能够有效地抑制层板的分层损伤扩展;相同冲击能量下缝合与未缝合层板的基体损伤和纤维损伤在厚度分布上相似,缝合层板的损伤都要小于未缝合层板。     

Impact and compression after impact experimental study of a composite laminate with a cork thermal shield

The aim of this paper is to present an experimental study of impact and compression after impact (CAI) tests performed on composite laminate covered with a cork thermal shield (TS) intended for launchers fairing. Drop weight impact tests have been performed on composite laminate sheets with and without TS in order to study its effect on the impact damage. The results show the TS is a good mechanical protection towards impact as well as a good impact revealing material. Nevertheless, totally different damage morphology is obtained during the impact test with or without TS, and in particular at high impact energy, the delaminated area is larger with TS. Afterwards, CAI tests have been performed in order to evaluate the TS effect on the residual strength. The TS appears to increase the residual strength for a same impact energy, but at the same time, it presents a decrease in residual strength before observing delamination. In fact, during the impact tests with TS, invisible fibres’ breakages appear before delamination damage contrary to the impacts on the unshielded sheets.     

Micro-mechanical finite element analysis of Z-pins under mixed-mode loading

This paper presents a three-dimensional micro-mechanical finite element (FE) modelling strategy for predicting the mixed-mode response of single Z-pins inserted in a composite laminate. The modelling approach is based upon a versatile ply-level mesh, which takes into account the significant micro-mechanical features of Z-pinned laminates. The effect of post-cure cool down is also considered in the approach. The Z-pin/laminate interface is modelled by cohesive elements and frictional contact. The progressive failure of the Z-pin is simulated considering shear-driven internal splitting, accounted for using cohesive elements, and tensile fibre failure, modelled using the Weibull’s criterion. The simulation strategy is calibrated and validated via experimental tests performed on single carbon/BMI Z-pins inserted in quasi-isotropic laminate. The effects of the bonding and friction at the Z-pin/laminate interface and the internal Z-pin splitting are discussed. The primary aim is to develop a robust numerical tool and guidelines for designing Z-pins with optimal bridging behaviour.     

复合材料缝编肋片的分析

为了提高复合材料机翼的整体强度和刚度, 降低制造成本, 机翼的主要结构件采用缝编技术。肋片由肋腹板、缘条组成, 它是由整块缝编层合板剪裁加工制成。肋水平缘条和垂直缘条分别缝纫于蒙皮及桁腹板上, 桁缘条也缝纫于蒙皮上。采RTM 工艺成型。当在肋腹板方向施加拉力后, 肋水平缘条和蒙皮之间存在着剪应力, 引起层合板的分层, 肋垂直缘条和桁腹板之间存在着拉应力, 引起它们之间的撕裂。采用了二维和三维有限元混合模型进行分析, 获得了肋水平缘条与蒙皮、以及肋垂直缘条与桁腹板交界面上的应力分布, 并由试验研究验证了分析计算的正确性。      

Modelling delamination onset and growth in pin loaded composite laminates

A three-dimensional (3D) finite element (FE) model is created with cohesive zone elements (CZE) to simulate a mechanically fastened [0°/90°]_s pin-loaded joint in a composite laminate. The model incorporates fully integrated solid elements in the pin-loaded area to accurately capture the high stress gradients. Contact based cohesive elements with a bilinear traction–separation law are inserted between the layers to capture the onset and growth of delamination. The stress distribution around the pin-loaded hole was verified with the widely used cosine stress distribution model. Results from the FE model show that delamination damage initiated at the point of maximum average shear stress at the 0°/90° interface. The delaminated area develops an elliptical shape which grows in a non-self similar manner with increasing pin displacement. It is concluded that a progressive damage model should be included to provide a full understanding of the failure sequence, work that the authors are currently engaged with.     

Fastener pull-through in a carbon fibre epoxy composite joint

This paper considers the capability of finite element (FE) modelling to accurately predict fastener pull-through failure of composite laminates. Such failures are dominated by inter-ply delamination and through-thickness shear failure of the laminate and the common modelling approach is to use computationally expensive, detailed three-dimensional models that include delamination for every ply interface, fastener contact and prestress. This paper considers a simplified FE modelling strategy achieved through judicious use of symmetry boundary conditions, hybrid shell/solid modelling and reduced numbers of interfaces for delamination. The LS-DYNA FE software was used for this study using the available composite material and cohesive failure models. The conclusion drawn from this work is that the use of simplified FE models does have merit in modelling fastener pull-through provided the material is quasi-isotropic and the boundary conditions are uniform around a circular perimeter. Additional work is however required to determine suitable cohesive properties and progressive shear failure parameters.     

Delamination Effect on Bending Behaviour in Carbon–Epoxy Composites

Abstract: This paper presents the results of experimental and numerical studies on the bending behaviour of carbon fibre‐reinforced epoxy composite containing delamination located at different positions along the laminate thickness. Experimental tests were conducted in three‐point bending using specimens with and without delamination to evaluate their bending behaviour. Numerical simulations were also performed in order to evaluate the maximum load as a function of the defect position and its size. The numerical model includes two‐dimensional solid elements of the ABAQUS software and a cohesive mixed mode damage model to simulate delamination propagation. The numerical and experimental results concerning the maximum load were found to be concordant. It was concluded that delaminations affect the bending behaviour of laminates mainly due to alterations in shear stress profiles.     

Motorcyclist helmet composite outer shell characterisation and modelling

In this study a new finite element model of composite outer shell of motorcyclist helmet is proposed, by modelling each layer of the composite material that builds the laminated structure of the outer shell of the helmet. Elastic and rupture properties of the laminate are taken into account for developing the finite element (FE) model and are extracted experimentally. A coupled experimental–numerical method combined with experimental modal analysis on beam samples is used to obtain the elastic characteristics of each layer of the outer shell. The rupture properties for each layer are extracted by experimental impact tests. The FE model of the outer shell is then validated with experimental data for elastic and rupture behaviour.     

Prediction of delamination in braided composite T-piece specimens

This paper presents an approach to predict the delamination of braided composite T-piece specimen using cohesive models. As part of an investigation on simulation of delamination in T-piece specimens, cohesive elements from ABAQUS were employed in forming a cohesive model to study the progressive delamination. Predictions given by the model of single delamination together with experimental results are presented. These results suggest that prediction of progressive delamination using cohesive models is feasible. Finally this paper proposes future work for precise prediction of delamination of braided composite T-piece specimens.     

Unification of strength scaling between unidirectional,quasi-isotropic,and notched carbon/epoxy laminates

An investigation into size effects and notch sensitivity in quasi-isotropic carbon/epoxy laminates was carried out. The purpose is to draw a complete picture of the strength scaling in unidirectional, quasi-isotropic, and notched carbon/epoxy laminates. A link was established between the strength scaling of the unidirectional and quasi-isotropic laminates. Efforts were made to understand the relationship between unnotched and open-hole strengths. For very small holes, the notched strengths approach the unnotched strength limit. A scaling law based on Weibull statistics was used to predict the unnotched laminate strengths. For very large holes, the same scaling law in conjunction with a detailed 3D ply-by-ply FE analysis with matrix cracks in the 90° plies and delamination cohesive interface elements was used to predict the large notched strengths. A good agreement between the modelling and experimental results was achieved. The effects of 90° matrix cracks on unnotched and notched strengths were also studied.     

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.     

Modelling complex progressive failure in notched composite laminates with varying sizes and stacking sequences

The emergence of advanced computational methods and theoretical models for damage progression in composites has heralded the promise of virtual testing of composite structures with orthotropic lay-ups, complex geometries and multiple material systems. Recent studies have revealed that specimen size and material orthotropy has a major effect on the open hole tension (OHT) strength of composite laminates. The aim of this investigation is develop a progressive failure model for orthotropic composite laminates, employing stepwise discretization of the traction–separation relationship, to predict the effect of specimen size and laminate orthotropy on the OHT strength. The results show that a significant interaction exists between delamination and in-plane damage, so that models without considering delamination would over-predict strength. Furthermore, it is found that the increase in fracture toughness of blocked plies must be incorporated in the model to achieve good correlation with experimental results.     

On the Distribution of Delamination in Composite Structures and Compressive Strength Prediction for Laminates with Embedded Delaminations

In this study, large numbers of aircraft composite structures were inspected, and the distribution of delamination sizes and though thickness positions in the composite laminates are investigated. An experiment is conducted to probe into the influence of delamination sizes and through thickness positions on the compressive strengths of laminates with single embedded circular delamination, with the most dangerous delamination sizes and positions defined from the distribution. Furthermore, a shell model is established for compressive strength prediction, with delamination propagation assessed using a mixed mode criterion. The finite element (FE) prediction comes out to be in good agreement with the experimental measurements, for the predicted compressive strengths stand within 10% error of experimental results. It was observed that the compressive strength was highly influenced by the delamination size, while the through thickness position of delamination did not have significant effect on the compressive strength.     

Compressive behaviour of thin-skin stiffened composite panels with a stress raiser

The in-plane compressive behaviour of thin-skin stiffened composite panels with a stress concentrator in the form of an open hole or low velocity impact damage is examined analytically. Drop weight impact in laminated polymer composites causes matrix cracking, delaminations and fibre breakage, which together can seriously degrade the laminate compressive strength. Experimental studies, using ultrasonic C-scan images and X-ray shadow radiography, indicated that the overall damage resembles a hole. Under uniaxial compression loading, 0° fibre microbuckling surrounded by delamination grows laterally (like a crack) from the impact site as the applied load is increased. These local buckled regions continued to propagate, first in discrete increments and then rapidly at failure load. The damage pattern is very similar to that observed in laminated plates with open holes loaded in compression. Because of this resemblance, a fracture mechanics model, developed initially to predict notched compressive strength, was applied to estimate the compression-after-impact (CAI) strength of a stiffened panel; in the analysis the impact damage is replaced with an equivalent open hole. Also, the maximum stress failure criterion is employed to estimate the residual compressive strength of the panel. The unnotched compressive strength of the composite laminate required in the analysis is obtained from a three-dimensional stability theory of deformable bodies. The influence of the stiffener on the compressive strength of the thin-skin panel is examined and included in the analysis. A good agreement between experimental measurements and predicted values for the critical failure load is obtained.