积分点的应力梯度怎么算？

xsxcn

请教一下大家，怎么算积分点的应力/应变梯度（就是应力/应变对空间的偏导数）。
我用的是六面体八节点单元。

bbssbb

1。高阶位移插值
2。看看还能不能分步积分了
3。如果上面方法都不行，就引入新的未知量，用拉格朗日乘子弱约束一下。