gradient crystal plasticity computational numerical FEM

Hi everyone,

I am recently implementing gradient crystal plasticity using FEM. As I understand, in order to obtain gradient of some plastic parameters (e.g. gradient of slip), additional dof (e.g. slip) should be added to each node and be interpolated with shape funcitons. And unlike conventional plasticity, which slip increment is computed locally at every integration point with stress update algorithm, both diaplacement increment and slip increment is computed globally.