使用Abaqus的VUMAT子程序做内聚力单元的本构，abaqus会传入一个strainInc(nblock，ndir + nshr )应变增量，我输出了ndir=1，nshr=2，刚好第一个是法向应变增量，后面两个是切向应变增量，然后切向的是真实剪应变，需要乘2倍后换成工程剪应变再使用。

对于法向应变增量计算方法就是内聚力单元上下表面的法向分离量除以厚度（默认为1）；我还查看了帮助手册，里面说两个切向的名义应变是积分点的切向位移除以厚度（默认为1），但我用这种方法计算出来的和abaqus传入的并不相等。

所以请各位大佬解惑，这个VUMAT传入的剪应变是怎样计算出来的呢？主要针对COH3D6零厚度内聚力单元。