abaqus vumat子程序计算损伤变量出现inf极大值问题？ 200

复合材料层合板准静态压缩，子程序调试了许多次，限制了损伤变量范围在0-1，然而实操下来结果大概在e30次方的范围，找不到子程序的问题在哪里？

邀请回答我来回答