In this paper a selected type of elasto-viscoplastic constitutive equations is considered. The viscoplastic model which was proposed by Marquis is generalized by allowing multiple terms describing the isotropic and the kinematic hardening. Furthermore, the presented model formulation enables one to use an arbitrary equivalent plastic strain function to describe the isotropic hardening behavior. What is more, the backstress evolution equation which was proposed by Marquis was modified so that any equivalent plastic strain function can be used in the recovery term now. The general form of the generalized constitutive model obtained this way was subsequently implemented into the finite element method (FEM). For that purpose, the radial-return mapping algorithm was utilized. The consistent tangent operator was derived for the considered class of models and is presented in the paper. A developed user material subroutine (UMAT) which allows one to use the viscoplastic models under consideration in the FEM program CalculiX is attached in the appendix section. A number of numerical simulations were conducted in order to verify the performance of the developed UMAT code.