user subroutines

Dear everyone:

  I am trying to use a nonlocal damage criteria to simulate the crack growth, which can reduce the mesh dependency.

In abaqus, VUMAT/UMAT is doing computation at integration point for each increment, but I want to do the post-processing after each increment to get the non-local state variable and return this value back to calculation.    

Such as taking average value of one state variable within a distance about the center of one element, then return it back to another state variable to this element.