abaqus 蠕变本构方程子程序二次开发？

求助这样的一个蠕变本构模型怎么用CREEP编写，大致看懂了帮助文档的实例程序，求大佬点拨下。

ε_c = Asinh{B(σ-σ₀)(1-H)}

 H=[1-H/H*]ε_c

ε_{c为蠕变速率，H为蠕变第一阶段应变速率变化，σ为应力，其他为常数}