In vibration fatigue, dynamic excitation such as base acceleration or force can produce a multiaxial stress response. In the frequency domain, this stress response is characterised by a cross-power-spectral-density (cross-PSD) matrix that contains the spectral content of each stress component and the phase relations between them. For fatigue-life estimation under stationary excitation, the multiaxial stress cross-PSD is reduced to an equivalent uniaxial stress PSD compatible with established spectral damage models. Several criteria have been proposed for this reduction, each differing in the assumptions made about phase relations, the required material parameters, and computational cost. This paper provides an implementation-oriented review of ten such criteria: three critical-plane variants (maximum normal, maximum shear, and combined normal-and-shear stress), the equivalent von Mises stress and its adaptation for out-ofphase stress components, the Carpinteri–Spagnoli frequency-domain reformulation, the frequency-based multiaxial rainflow, a thermoelasticity-based criterion, the Niesłony hydrostatic–deviatoric combination, and the equivalent Lemaitre stress with multiaxial S–N interpolation. As a result of this research, all ten methods are also implemented in FLife (version 2.2.0), an open-source Python package that extends the uniaxial spectral framework to multiaxial inputs. Previously these implementations were distributed across separate source papers; here they share a single input format and a single workflow, so that any subset can be compared on the same dataset without reimplementation. A comparative overview highlights trade-offs between simplicity, sensitivity to non-proportional loading, material-parameter requirements, and computational cost, supporting method selection in multiaxial vibration-fatigue practice. The criteria are also applied to a common dataset, where the equivalent stress varies by up to about a factor of 1.7 across the ten criteria.



