In this paper, a simple efficient method for the numerical solution of a class of nonlinear Volterra integral equations (VIEs) is presented. The approach starts by expanding the existing functions in terms of Bernoulli polynomials. Subsequently, using the new introduced Bernoulli operational matrices of integration and the product along with the so-called collocation method, the considered problem is reduced into a set of nonlinear algebraic equations with unknown Bernoulli coefficients. The error analysis and rate of convergence are also given. Finally, some tests of other authors are included and a comparison has been done between the results.