Multiaxial high cycle fatigue modeling of materials is an issue that concerns many industrial domains (automotive, aerospace, nuclear, etc.) and in which many progress still remains to be achieved. Several approaches exist in the literature: invariants, energy, integral and critical plane approaches all of them having their advantages and drawbacks. These different formulations are usually based on mechanical quantities at the micro or mesoscales using localization schemes and strong assumptions to propose simple analytical forms. This study aims to revisit these formulations using a numerical approach based on crystal plasticity modeling coupled with explicit description of microstructure (morphology and texture) and proposes a statistical procedure for the analyses of numerical results in the HCF context. This work has three steps: First, 2.5D periodic digital microstructures based on a random grain sizes distribution are generated. Second, multiaxial cyclic loading conditions corresponding to the fatigue strength at 106 cycles are applied to these microstructures. Third, the mesoscopic Fatigue Indicator Parameters (FIPs), formulated from the different criteria existing in the literature, are identified using the finite element calculations of the mechanical fields. These mesoscopic FIP show the limits of the original criteria when it comes to applying them at the grain scale. A statistical method based on extreme value probability is used to redefine the thresholds of these criteria. These new thresholds contain the sensitivity of the HCF behavior to microstructure attributes. Finally, the biaxiality and phase shift effects are discussed at the grain scale and the loading paths of some critical grains are analyzed.Read less <