Modeling of hardening and thermal recovery in metals is considered within the context of unified elastic-viscoplastic theories. Specifically, the choices of internal variables and hardening measures, and the resulting hardening response obtained by incorporating saturation-type evolution equations into two general forms of the flow law are examined. Based on the analytical considerations, a procedure for delineating directional and isotropic hardening from uniaxial hardening data has been developed for the Bodner-Partom model and applied to a nickel-base superalloy, B1900 + Hf. Predictions based on the directional hardening properties deduced from the monotonic loading data are shown to be a good agreement with results of cyclic tests.