[en] In this paper, we present a stress update scheme for non-deviatoric elastoviscoplastic models, which can be considered as an extension of the radial return method for classical plasticity. We use the unconditionally stable backward Euler scheme to obtain the viscoplastic solution at each increment. To solve the large sparse non symmetric iteration matrix, we use the GMRES method. This fast iterative solver, in conjunction with an appropriate preconditioner, can save a lot a memory and CPU time when 3D problem are considered. A comparison between direct and iterative solvers will be discussed on two applications of cold compaction of powders.