The subloading surface model has been formulated and applied to the prediction of cyclic loading behavior of metals. In this article, the return-mapping in the fully implicit algorithm based on the extended subloading surface model is formulated aiming at obtaining accurate solutions in a high efficiency in the whole loading state. It is incorporated into the implicit finite element program Abaqus through the User-subroutine UMAT. Simulation result by the fully implicit algorithm is compared to the forward Euler algorithm about the accuracy and the convergence property under the cyclic loading.