This work evaluates the ability of cohesive zone modeling-based approaches to predict delamination in composite materials that develop large process zones under complex loading conditions. The R-curve effects subjected to static and fatigue loading under multiple loading modes, considering the loading history, are analyzed. To this end, the delamination predictions of a state-of-the-art CZM-based simulation strategy are evaluated by blind simulation of a validation benchmark test. The validation test promotes a non-self-similar delamination scenario, including a process zone that evolves under different loading mode conditions with a non-straight leading delamination front. Good delamination prediction accuracy is achieved. In addition, insights into the relationship between the features of the simulation strategy and the physics of the delamination process are discussed. With regard to the limitations of the simulation strategy, particular attention should be paid to modeling the contribution of an evolving process zone based on the loading mode history.

The main objective of this paper is to present a delamination benchmark test concept for composite materials that develop non-self-similar delamination in characterization specimens. The non-self-similar delamination is induced by rotating the loading blocks. The simplicity of the test allows for analyzing the loading mode history by concatenating different loading conditions, such as static and fatigue loading, under multiple loading modes. The methodology introduced in this paper can be particularized for any given composite material set and any sequence of loading conditions. To demonstrate the capabilities of the benchmark test, a case study is presented using AS4D/PEKK-FC thermoplastic composite material, which exhibits strong R-curve behavior. A sequence of opening and shear failure modes was applied under static and fatigue loading, providing an experimental data set that is ready to be used as a part of the validation of numerical predictive delamination models. The delamination process was monitored by X-ray radiography, and the final fracture surfaces were analyzed with scanning electron microscopy (SEM), giving a physical insight into the contribution of the fracture mechanisms to the delamination process.

The complex failure mechanisms involved in failure of interfaces requires the use of an accurate description of the cohesive law. In recent years, there have been many developments to determine the full shape of the cohesive law. However, most of the existing cohesive zone models assume a simplified shape, such as bilinear, trapezoidal or exponential, which are usually simple to model. Their accuracy is found to be rather limited, especially in the presence of a large fracture process zone due to either plastic deformation or fibre bridging. In this work, a new cohesive element description is proposed to formulate a general cohesive zone model to overcome these limitations. The benefit of the new approach is that it allows for convenient implementation of any arbitrary shape of the cohesive law obtained experimentally. The authors present a new procedure based on the superposition of n-bilinear cohesive zones to obtain an equivalent multilinear cohesive law that fits any experimental measurement. The new element formulation has been implemented in the commercial finite element software ABAQUS, using user element subroutine. Verification of the methodology is performed at the single element level and the approach is validated for different material systems (adhesives and composites) using the double cantilever beam, end-notched flexure and mixed-mode bending tests. Excellent correlation between all numerical predictions and experimental results is obtained, demonstrating the robustness of the proposed methodology.

Through the application of innovative production processes, thermoplastic composites might help the aviation industry become more sustainable. However, there is currently not much experimental understanding on the fatigue behaviour, and validated analysis methodologies on thermoplastic composites are rather limited. In this work, the fatigue onset and propagation interlaminar properties of AS4D/PEKK-FC thermoplastic composite were characterized under mode I loading. Testing costs were reduced by using a multi-fatigue testing rig that allows loading 6 specimens at the same time. The offset in the fatigue crack growth rate curves tested at different severities was explained by the fatigue R-curve effects. The offset was modelled by means of superposed fatigue cohesive laws. Different approaches on how to consider fatigue damage were compared with experimental data, giving an insight on how fatigue damage evolves, and developing a novel modelling strategy based on a robust method for fatigue model parameter identification for fatigue delamination when several failure mechanisms interact.

Thermoplastic composites can enable the development of new manufacturing techniques to make the aviation industry more sustainable while at the same time greatly benefit cost-efficient and high-volume production. One of the thermoplastic composite materials that can enable this transition is AS4D/PEKK-FC. In this work, the interlaminar properties of AS4D/PEKK-FC thermoplastic composite are characterized and analyzed by means of Mode I, II and Mixed Mode I/II at 50:50 tests, while considering fiber bridging and R-curve effects. In order to achieve stable crack propagation the test configurations are adjusted to account for the large fracture process zone ahead of the crack tip and an appropriate data reduction method is selected. The experimental data is reduced using an inverse methodology to extract cohesive laws based on only the load–displacement curves. Additionally, the use of this methodology provides new insights into the validity of two different mode II tests and the influence of fiber bridging on the mixed-mode interlaminar behavior. The interlaminar damage mechanisms are investigated by means of scanning electron microscopy. The resulting cohesive laws are implemented in commercial finite element software in tabular form, without the need for user-subroutines. All experimental test configurations are analyzed using a single material card and it is shown that fiber bridging and R-curve effects are well captured.