An inverse method to reconstruct crack-tip cohesive zone laws for fatigue by numerical field projection

Cohesive zone failure models are widely used to simulate fatigue crack propagation under cyclic loadings, but the functional form of such phenomenological models does not have a strong micromechanical basis. While numerical techniques have been developed to extract the crack-tip cohesive zone law for monotonic fracture from experimental strain information, elucidating the unloading-reloading hysteresis representing accumulated and irreversible damage under fatigue cycling is still challenging. Here, we introduce a novel field projection method (FPM) to reconstruct the cohesive zone laws for steady-state fatigue crack growth from elastic strain field information at various loading and unloading stages within a single steady-state fatigue cycle. The method is based on the Maxwell-Betti's reciprocal theorem, which considers the material response to be linear elastic during loading and unloading, albeit with a reciprocity gap to account for the inhomogeneous residual elastic strain accumulated at the end of each fatigue cycle. Through numerical experiments, we demonstrate that this FPM is capable of accurately extracting the crack-tip cohesive tractions and separations, along with the elusive unloading-reloading hysteresis, which constitute the cohesive zone law for fatigue crack growth. We discuss the errors and uncertainties associated with this inverse approach.