This paper illustrates how to implement effectively solvers for elasto-plastic problems. We consider the time step problems formulated by nonlinear variational equations in terms of displacements. To treat nonlinearity and nonsmoothnes we use semismooth Newton method. In each Newton iteration we have to solve linear system of algebraic equations and for the numerical solution of the linear systems we use TFETI algorithm. In our benchmark we compute von Misses plasticity with isotropic hardening and use return mapping concept