ABAQUS UMAT DIFFUSION FICK

Hello.

     Could anyone tell me how to  achieve the basic unsteady diffusion equation (Fick's second law) of some solute for a given Abaqus model by UMAT & User-defined field in Abaqus? or some other subroutines?

    Is there any reference good for this problem ?

  Thanks very much!