The fast multipole boundary element methods (FM BEMs) are developed by means of the multipole expansions in conjunction with the implicit solution GMRES (generalized minimum residual method) of linear system arising from BEM. The Taylor series multipole boundary element method (TSM BEM) updates the solution mode of conventional BEM, enhances the computational efficiency rapidly, and enlarges the scale of solution. The development and research status of TSM BEM are introduced, the basic ideas and principles are described, and the classify methods are presented. The Taylor series expansions for fundamental solutions and basic procedures of boundary elements integrals are presented for multipole BEM. The numerical method is applied to rolling engineering. The elastic deformation of rollers and the contact and deformation analysis of HC mill roll system are provided. The numerical example demonstrates that the TSM BEM is suitable for the solution of large scale rolling engineering problems.