Models proposed to study the multiaxial fatigue damage phenomenon generally lack probabilistic interpretation due to their deterministic form. This implies failure compulsory happening at the plane exhibiting the maximum damage value, whereas the remaining planes are disregarded. Nevertheless, the random orientation of the predominant defect evidences the possibility of failure being initiated as a function of the predominant defect presence without requiring, necessarily, maximum values of the damage parameter, which emphasizes the need of introducing probabilistic concepts into the failure prediction analysis. In this paper, a probabilistic model is presented that enables the failure probability to be found for any selected plane orientation by considering the damage gradient as a parameter for both proportional and non-proportional loading. The applicability of the model is elucidated by means of an example. Assuming the cdf for the local failure of the material to be known, the probability of failure is calculated for a cross shaped specimen in which shift between the principal stresses sigma xx and sigma yy ranges from 0º to 180º.