How to include Hill Stress/ Anisotropic stress using SUBROUTINE CREEP?

Hello,

I am trying to do anisotropic creep simulation using subroutine creep. For validation purpose I am simulating hyperbolic sine creep law with hill ratio. I can do simulation with default Abaqus material library. But when I try to modify subroutine creep then it gives me an error. Also I am able to use subroutine creep for sine hyperbolic creep.

When I go through this reference: https://classes.engineering.wustl.edu/2009/spring/mase5513/abaqus/docs/v6.6/books/sub/default.htm?startat=ch01s01asb01.html

I find that QTILD can take both VM and Hill stress. So when I input potentials in creep library it also gives me error.

Is there a way to add hill ratios in subroutine creep code and compute QTILD in the code?

what are other options I can try using SUBROUTINE CREEP to do anisotropic creep simulation.

Please let me know.

Thanks