For a large class of elastoplastic models, the present paper proposes an integration scheme which updates stress points on yield surfaces automatically. Associated plastic flow models with the rules of kinematic, isotropic, and distortional hardening/softening are contained in the model class. The underlying structure of the class of elastoplastic models is exposed and utilized to develop the return-free integration scheme that keeps the stress points on the yield surfaces without any extra enforcement. The return-free capability of the integration scheme is underpinned by the theory of Lie group and Lie algebra. Numerical demonstrations show that the stress points are updated on the yield surfaces automatically and exactly within the machinery round-off error.