High fatigue load, which exists widely in steel building structures, likely leads to brittle failure at the joints, supports, and so on. This can lead to the partial or total damage of the structure and even to cause the collapse of the whole structure. This article aims to provide a method to simulate high-cycle crack propagation in tubular joints, which is one of the most common types occurring in steel structures. Firstly, sixteen T-shaped tubular joint models under different load conditions and initial crack dimensions were built through the coordinate mapping method. Secondly, based on the extended finite element method (XFEM), an algorithm was developed by combining the secondary development in Abaqus and a quasistatic simulation method to simulate high-cycle crack growth in tubular joints under a constant amplitude. The results of the simulations were compared with experimental data. The study found that the surface stress calculated from the tubular joint models using the coordinate mapping method was close to the experimental data. Through the comparison of the crack propagation rate and the crack growth process between the simulation and experiment results, the simulation method was validated. When a crack penetrated the tube wall, the difference in the load cycles between the simulations and the experiment was 9.5%. The initial crack dimension had an impact on the crack propagation, with the decrease in the a/c and KII generally becoming the dominant factor with respect to the crack growth, while the fatigue life of the joints tended to increase.