A relaxation modulus-based matrix splitting iteration method is proposed for solving a class of nonlinear complementarity problems. The convergence theory is established when the system matrix is H+and the choice of relaxation parameters is given. Numerical examples show that the proposed methods are efficient and can accelerate the convergence performance of the modulus-based matrix splitting method with less iteration steps and CPU time.