In this paper, a numerical simulation method of mixed-mode fatigue crack propagation was explored using the extended finite element method (XFEM) and the Virtual Crack Closure Technique (VCCT). Both 2D and 3D numerical models were selected to simulate the fatigue crack propagation of steel specimens. Two coefficients were proposed to calculate the equivalent energy release rate (Geq) for a better simulation of the mixed-mode fatigue crack propagation of S355 steel. The Walker equation and the calculation formula of Geq were realized by a user-defined subroutine. A set of optimal correction coefficients of mode II energy release rate (GII) and mode III energy release rate (GIII) were quantitatively comparing the simulation results and test data. The results will contribute to fatigue crack propagation prediction of steel structures in the civil engineering field.@en