关于Abaqus中隐式蠕变子程序的问题

的隐式的在abaqus中实现蠕变需要dh/d(\delta \epsilon_c)作为牛顿-拉夫森迭代中的雅可比矩阵。导数可以展开，如图所示。用绿色笔画出的术语dh/dq，期望用户通过指定DECRA来提供，而用红色笔画出的术语dq/d(\delta \epsilon)让我感到困惑。这一项似乎是由Abaqus计算出来的，但我不知道它是如何完成的。同样的问题也适用于右边最后一个导数。

非常感谢!