The stability of some numerical time integration schemes used in the quasistatic finite element analysis of hardening elasto-viscoplastic structures is investigated. The discrete problem is reduced to a set of ordinary differential equations to which the general theory is applied. Unconditional stability is proved for an implicit algorithm and explicit stability criteria which can be very easily introduced in a finite element code, are derived for the gradient and Euler's procedures.