In this paper, we consider the Galerkin method to approximate the solution of Volterra-Urysohn integral equations of second kind with a smooth kernel, using piecewise polynomial bases. We show that the exact solution is approximated with the order of convergence O(hr) for the Galerkin method, whereas the iterated Galerkin solutions converge with the order O(h2r) in uniform norm, where h is the norm of the partition and r is the smoothness of the kernel. For improving the accuracy of the approximate solution of the integral equation, the multiGalerkin method is also discussed here and we prove that the exact solution is approximated with the order of convergence O(h3r) in uniform norm for iterated multi-Galerkin method. Numerical examples are given to illustrate the theoretical results. © Springer Nature Singapore Pte Ltd. 2017.