The paper presents the method of determination of two-dimensional probability distribution P_f of crack initiation versus fatigue life N and the fatigue damage parameter : P_f − N − ϭ. The proposed distribution P_f 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.