abaqu通过子程序C3D8I积分点赋予

通过坐标赋予每个element状态变量作为该分区的材料属性，用C3D8单元就没有问题如上图，但是用C3D8I单元就会出现下面的情况，请各位解答，我这个SDV92我是直接通过子程序赋予的刚度矩阵系数，由于这汇总问题的出现，出现负数，会导致计算不收敛。