A quasi‐static boundary value problem in multi‐surface elastoplasticity: part 2—numerical solution