求高手帮我把下面的公式转化为FORTRAN子程序,我在ABAQUS模拟中要调用,主要用该公式来描述物体的变形.温度变化范围20-1200℃,10℃/min.谢啦。。。