Abstract. We study the numerical approximation of a static innitesimal plasticity model of kinematic hardening with a nonlocal extension involving the curl of the plastic variable. Here, the free energy to be minimized is a combination of the elastic energy and an additional term depending on the curl of the plastic variable. In a rst step, we introduce the stress as dual variable and provide an equivalent primal-dual formulation resulting in a local ow rule. To obtain optimal a priori estimates, the nite element spaces have to satisfy a uniform inf-sup condition. Finally, we show that the associated nonlinear mixed formulation can be solved iteratively by a classical radial return algorithm. Numerical results illustrate the convergence of the applied discretization and the solver.Key words. gradient plasticity, nite element estimates AMS subject classications. 65M12, 65M60, 65N22, 74H15, 74S051. Introduction. The abstract setting for variational inequalities provides a powerful framework for the analysis of innitesimal plasticity and its nite element discretization, see [11] and the references therein. In the static case, the elastoplastic solution is determined by minimizing a (primal) functional for the displacement and the plastic variable. Unfortunately for perfect plasticity this functional is not uniformly convex, and the minimizer exists only in a weak sense. We refer to [28] where suitable Banach spaces for the displacement and the plastic variable are discussed. On the other hand, the idealistic model of perfect plasticity does not include hardening nor size eects reecting internal length scales. There are dierent possibilities for including such eects, and in most cases this leads to a more regular model with an associated uniformly convex primal functional, and thus the theory of standard Sobolev spaces can be applied.Here, we consider a class of nonlocal models which can be obtained by an extension of the classical plasticity model with kinematic hardening. More precisely, we study a subclass of gradient plasticity, where the corresponding elasto-plastic energy includes only the curl of the plastic strain, see, e.g., [7,8,9,13,14,26,27] for the analytical, mechanical, and physical properties of such models. In [17,18] it is shown that a representative nonlocal model of such type can also be transformed into a variational inequality.The nite element analysis for variational inequalities in plasticity with kinematic hardening yields a priori estimates involving the best approximation error and additional nonconformity terms, see [1,3]. Here, we extend these results to our non-local formulation and work with curl-conforming nite elements for the plastic variable. A dierent approach for a discontinuous Galerkin formulation of a gradient plasticity model with full gradient terms is presented in [5,6,20], and in [19] a discontinuous Galerkin approximation of a model containing the curl of the plastic strain is considered.Following [11], we consider plasticity models which are completely determined by