  This study provides an abstract framework to analyze mixed formulations in
viscoelasticity, in the classic saddle point form. Standard hypothesis for
mixed methods are adapted to the Volterra type equations in order to obtain
stability of the proposed problem. Error estimates are derived for suitable
finite element spaces. We apply the developed theory to a bending moment
formulation for a linear viscoelastic Timoshenko beam and for the Laplace
operator with memory terms. For both problems we report numerical results to
asses the performance of the methods.
