Summary: In the context of infinitesimal strain plasticity with hardening, we derive a stochastic homogenization result. We assume that the coefficients of the equation are random functions: elasticity tensor, hardening parameter and flow-rule function are given through a dynamical system on a probability space. A parameter \(\varepsilon > 0\) denotes the typical length scale of oscillations. We derive effective equations that describe the behavior of solutions in the limit \(\varepsilon\to 0\). The homogenization procedure is based on the fact that stochastic coefficients “allow averaging”: For one representative volume element, a strain evolution \([0,T]\ni t \mapsto \xi(t)\in \mathbb{R}^{d\times d}_s\) induces a stress evolution \([0,T]\ni t \mapsto \varSigma(\xi)(t)\in \mathbb{R}^{d\times d}_s\). Once the hysteretic evolution law \(\varSigma\) is justified for averages, we obtain that the macroscopic limit equation is given by \(-\nabla \cdot \varSigma(\nabla^s u) = f\).