Can anyone help me to make the user subroutine UMATHT used in ABAQUS simulation for lithiation process in lithium ion battery (LIB)? I am working on anode part of LIB and need to run an FEM simulation. However, I am not aware of the programming part of this subroutine. I would really appreciate comments from people who know about it.