A variational multiscale method to incorporate strain gradients in a phenomenological plasticity model