It is demonstrated that the eXtended Finite Element Method (XFEM) is of remarkable efficiency in simulating crack evolution by eliminating the need for remeshing and refinement. In this paper, it is shown how to enhance the solution efficiency through a comprehensive mathematical investigation of the solution process using XFEM. A typical example is presented to illustrate the disparities in nodal displacements along the two symmetric faces of the crack resulting from the approximation of XFEM. By analysing the structure and components of the global stiffness matrix, the underlying causes of these discrepancies are identified. Building upon these findings, two improvements of the solution are proposed to gain an acceptable accuracy in computing the nodal displacements. The first improvement consists of the subdivision of the enriched elements depending on the characteristic of the distribution of Gauss points. The second improvement is set by determining the optimal number of Gauss points in each sub-element near the crack tip. To calculate the stress intensity factor of the crack under surface pressure, such improvements are applied in conjunction with the interaction integral method, which significantly reduces computational time and eliminates the influence of surface tractions. The numerical solution is validated by comparing it with the analytical solution and the standard XFEM solution. The proposed improvements can enhance both the accuracy of the solution and the computational efficiency of XFEM.