In this work, an initial boundary value problem for a system of viscoelastic wave equations with nonlinear boundary source term of the form (ui)tt - Delta(ui)(tt) + integral(t)(0) gi(t - s)Delta ui(s)ds - Delta(u(i))t = 0, in Omega x (0,T), ui(x, 0) = phi(i)(x), (u(i))t(x, 0) = psi(i)(x), in Omega, u(i) (x, t) = 0, on Gamma(0) x (0, T), partial derivative(v)(u(i))(tt) + partial derivative(v)u(i) -integral(t)(0) (t - s)partial derivative(v)u(i)(s)ds +partial derivative(v)(u(i))(t) + f(i)(u) = 0, on Gamma(1) x (0,T), where i = 1,..., l (1 >= 2) is considered in a bounded domain Omega in R-N (N >= 1). By the Faedo Galerkin approximation method we obtain existence and uniqueness of weak solutions. Under appropriate assumptions on initial data and the relaxation functions, we establish general decay and blow up results associated to solution energy. Estimates for lifespan of solutions are also given. (C) 2017 Elsevier Inc. All rights reserved.