The incremental finite element method is applied to find the numerical solution of the plasticity problem with strain-hardening. Following Watwood and Hartz, the stress field is approximated by equilibrium triangular elements with linear functions. The field of the strain-hardening parameter is considered to be piecewise linear. The resulting nonlinear optimization problem with constraints is solved by the Lagrange multipliers method with additional variables. A comparison of the results obtained...