我编写了一个C3D20单元的UEL，调用也可以运行成功，但是和abaqus自带的C3D20单元的结果对不上

由于UEL只能输出结点位移，因此我对比了两个C3D20单元的结点位移图，下图为abaqus自带的单元，上图为编写的单元，是一个网格数2×2×2的模型，约束条件为固定一个面z方向位移为0，另一个面的z方位位移为-5。按理来说除了约束住位移的节点之外，其他节点的z方向位移应该都指向z轴反方向且由大变小，abaqus自带单元算的结果正是如此；但是我编写的UEL子程序，除了约束住位移的节点之外，其他节点的z方向位移指向z轴反方向和正方向都存在，且节点位移也不是由大到小，很乱没有规律。我的问题就是以上内容，望有大佬能指点一二。万分感谢！！！！