100 abaqus蠕变子程序帮助文档中的案例公式推导 可有偿？

图片中是ABAQUS帮助文档中的一个双曲正弦蠕变本构方程（见图1上端：红色1处；红色2 处为其偏导数）。在子程序中，DECRA(1)为等效蠕变应变增量（即红色1等号左端形式），DECRA(5)为等效蠕变应变增量偏导增量（即红色2等号左端形式），q为Mises等效应力，A, SIG0, AN为常数，QTILD即为q等效应力，DTIME为时间增量，LEXIMP为显式和隐式蠕变积分的标识符，我也仔细推敲了下红色1处公式（详细推导见图2），但是图片...