A novel methodology is proposed for accelerating the calculation efficiency of the boundary element modeling of rolling contact. This methodology involves the implementation of an initial solution estimate. The method is to provide the initial estimate value by means of simplifying the method, which is used for the iterative calculation of the boundary element method to solve the normal and tangential contact problems. In the normal contact problem, the normal pressure and contact patch are provided as the initial values for the iterative calculation of the boundary element method. In the tangential contact problem, the initial values for tangential stress and adhesion-slip distribution are provided. A first novel aspect is breaking the initial iteration setting of the traditional BEM, which can significantly reduce the iterations. A second novelty is that a method for determining the potential contact area is proposed to ensure the correctness of boundary element modeling without increasing the computational cost. In the following section, an analysis is conducted to ascertain the impact of the initial solution estimate method on the efficacy of boundary element modeling. The result is that efficiency increased by 69.1% of normal contact calculations, with the initial contact patch exerting the most significant influence. The efficiency of the process under investigation increased by 56.9%, and the most pronounced effect is the distribution of adhesion-slip.