abaqus应力诱导氢扩散问题解决

输氢管道的氢扩散模拟，目前有子程序但是应力分布求解得到的结果为0，也就是没有应力分布，子程序使用的是unatht+umat，并使用disp子程序施加弹性裂纹位移场来施加循环位移载荷，但是在我这个空心管道上似乎无法使用。需要给出解决方案来计算静水应力梯度。