A numerical study on toughening laminated composite materials against delamination by initiating multiple interlaminar cracks is presented. Different configurations of interface toughening and weakening patches, that modify the interface properties at selected locations, are investigated as a way to trigger multiple delaminations. Both interface toughening and weakening patches can be configured to toughen the laminated material by initiating multiple delaminations. The initiation of multiple delaminations and the increase in toughness depend on the interface strengths and toughness of the patches. The main mechanisms behind the initiation of multiple delaminations for both cases are presented. An adaptive refinement method implemented within a Matlab Finite Element Analysis code that models the interfaces of the laminate with cohesive elements is used for the analyses. The adaptive refinement framework allows efficient analysis of multiple delaminations with very fine meshes at the wake of the crack tips. A discussion on the overall performance of the toughening concept, and the main parameters affecting the results, i.e. the length of the interface toughening or weakening patches, the distance of the substrate between the affected interfaces, and the material's mechanical properties, is carried out. The results presented in the paper show that a toughening effect against delamination can be achieved using interface toughening or weakening patches to onset multiple delaminations.

A novel and efficient numerical formulation for the modelling of multiple delaminations growth in laminated composite materials subjected to quasi-static loading is presented. The proposed formulation alleviates the high computational cost associated with models featuring cohesive elements by using a novel Adaptive Refinement Scheme and an Adaptive Floating Node Method Element to refine the model effectively during the analysis without modifying the global finite element connectivity. The formulation has been implemented in a MATLAB finite element code and validated with single and multiple delamination numerical models with varying mode mixities. The new formulation provides accurate results comparable to standard fully refined finite element models while drastically lowering the computational time of the analysis.

A novel efficient numerical formulation for the analysis of multiple fatigue-driven delamination cracks is presented. A cohesive zone model is used in combination with an Adaptive Refinement Scheme (ARS) and an Adaptive Floating Node Method (A-FNM) element that refine the model effectively during the analysis. Novel techniques are proposed to track the positions of multiple crack tips and calculate the mode decomposed energy release rates for the individual crack tips using the J-integral. The method has been implemented in a Matlab finite element code and validated with single and multiple delamination cases with varying mode mixities. Comparisons with theoretically based predictions and available experimental data showcase the high accuracy of the method. The presented method lowers the computational time compared to standard, fully refined finite element models by a factor of 4–5.