The aim of this paper is to study the theoretical and numerical stability of the Bresse system in one-dimensional bounded domain with viscoelastic Kelvin–Voigt damping. We first showed the well posedness of the system. Then stability is obtained by applying the multiplicative techniques. Later a numerical scheme is proposed and analyzed. Finally a priori error estimate is established.