The actual service axles are often subjected to mixed-mode loading, and predicting the mixedmode I-II crack propagation behaviour using the mode I effective stress intensity factor (SKI) differs from the real service conditions. To effectively predict the fatigue crack propagation behaviour of actual service structures, the I-II stress intensity factor range (SKP-R) considering two closure effects was adopted to describe the fatigue crack propagation under mixed mode loading. A test database was established based on monitoring data of mode I and mixed-mode I-II (30/45/ 60 degrees) crack propagation tests under different stress ratios. Combining domain knowledge and symbolic regression (SR) methods, an angle factor was proposed for constructing correlation functions between SKI and SKP-R. The results showed that the loading angle (alpha) only affects the initial projection of the load parallel and perpendicular to the fatigue crack growth (FCG) direction. Compared with the geometric correction factor, the correlation function acquired by the angle factor constructed by the SR method has higher accuracy, and the balance parameters (SCORE) obtained by the former are significantly higher than those obtained by the latter under the same function complexity. The SR verification results demonstrated that constructing mode I and I-II correlation functions with angle factors has a good predictive effect.