The paper presents the method of determination of two-dimensional probability distribution Pf of crack initiation versus fatigue life N and the fatigue damage parameter : Pf − N − ϭ. The proposed distribution Pf uses parameters of the standard fatigue characteristics and allows calculating fatigue life of elements with heterogeneous stress fields at any probability level. The model was successfully verified on experimental test results.