A robust numerical methodology for fatigue damage evolution simulation in composites

Composite materials, like metals, are subject to fatigue effects, representing one of the main causes for component collapse in carbon fiber‐reinforced polymers. Indeed, when subject to low stress cyclic loading, carbon fiber‐reinforced polymers exhibit gradual degradation of the mechanical properties. The numerical simulation of this phenomenon, which can strongly reduce time and costs to market, can be extremely expensive in terms of computational effort since a very high number of static analyses need to be run to take into account the real damage propagation due the fatigue effects. In this paper, a novel cycle jump strategy, named Smart Cycle strategy, is introduced in the numerical model to avoid the simulation of every single cycle and save computational resources. This cycle jump strategy can be seen as an enhancement of the empirical model proposed by Shokrieh and Lessard for the evaluation of the fatigue‐induced strength and stiffness degradation. Indeed, the Smart Cycle allows quickly obtaining a preliminary assessment of the fatigue behavior of composite structures. It is based on the hypothesis that the stress redistribution, due to the fatigue‐induced gradual degradation of the material properties, can be neglected until sudden fiber and/or matrix damage is verified at element/lamina level. The numerical procedure has been implemented in the commercial finite element code ANSYS MECHANICAL, by means of Ansys Parametric Design Languages (APDL). Briefly, the Smart Cycle routine is able to predict cycles where fatigue failure criteria are likely to be satisfied and to limit the numerical simulation to these cycles where a consistent damage propagation in terms of fiber and matrix breakage is expected. The proposed numerical strategy was preliminarily validated, in the frame of this research study, on 30° fiber‐oriented unidirectional coupons subjected to tensile– tensile fatigue loading conditions. The numerical results were compared with literature experimental data in terms of number of cycles at failure for different percentage of the static strength. Lastly, in order to assess its potential in terms of computational time saving on more complex structures and different loading conditions, the proposed numerical approach was used to investigate the fatigue behavior of a cross‐ply open‐hole composite panel under tension–tension fatigue loading conditions.