50Abaqus中如何计算零厚度内聚力单元COH3D6的应变？特别是VUMAT中？

使用Abaqus的VUMAT子程序做内聚力单元的本构，abaqus会传入一个strainInc(nblock，ndir + nshr )应变增量，我输出了ndir=1，nshr=2，刚好第一个是法向应变增量，后面两个是切向应变增量，然后切向的是真实剪应变，需要乘2倍后换成工程剪应变再使用。对于法向应变增量计算方法就是内聚力单元上下表面的法向分离量除以厚度（默认为1）；我还查看了帮助手册，里面说两个切向的名义应变是积分点的切向位移除以厚度...