Abstract: In order to solve a frictional contact problem, the isogeometric analysis combined with B differential equation method (IGA-BDE method) is developed. The geometry of boundaries of two contact bodies are represented accurately by NURBS. The contact conditions are formulated as a B differential equation and satisfied rigorously, and the convergence of algorithm to solve the B differential equation is guaranteed by the theory of mathematical programming. Numerical results validate the effectiveness and accuracy of the proposed IGA-BDE method to solve the elastic frictional contact problems. In contrast with the finite element method by employing an ANSYS contact model, the number of degree of freedom can be saved to a great extent.