Vumat损伤子程序 求指导

gepeng101319

最近编写了一个Vumat损伤子程序，计算考虑材料滞后回环能的损伤疲劳寿命问题。
问题中外载荷为周期性外载荷，滞后回环能计算过程中需要用到循环在和下的最大最小应力应变。
因此，将在和步分解，以是的显示计算中可以计算应力应变峰值，并在单次循环在和结束时计算构件滞后回环能，并计算损伤量。
算上计算公式及计算流程参考文献Multi-scale damage mechanics method for fatigue life prediction of additive manufacture structures of Ti-6Al-4V
问题：编写子程序后程序并没有错误，如果执行单次损伤量计算，可以有计算结果，但实际问题中英考虑为单个单元完全损伤才进行下一步计算
因此需将单次损伤进行循环计算，知道单元完全损伤位置，但加入循环后，时间0.01秒时计算一直不收敛（或者出错），无法完成计算
求大神帮忙看看子程序是否有问题，问题出现什么地方。
跪求大神指导

hfuter211

你好，请问问题解决了吗？求解释！

枯藤老树昏鸦

学习一下