各位老师好，最近在做单单元的压缩模拟，碰到问题特来请教：

如图，长度为10mm的立方体，底部节点组约束Z方向的位移，顶部节点组匀速向下运动，运动速度0.0001m/s（即应变率为0.0001/0.01=0.01），材料为钢，密度7900kg/m3，弹性模量200GPa，泊松比0.3.求解时间为0.3s，通过关键字输出单元信息共1000个。

速度加载关键字：

但是进行求解时求解1步，单元就出现无限大体积，模型也会无限变形，直接压扁

SingleElement.k

请问这个是什么原因？该如何解决呢？谢谢！

按照@田丰翼前辈的回答，修改之后结果如下：

2#算法（Y,Z视图）：

3#算法（Y,Z视图）：

看到这个结果突然想到可能约束有问题，于是对5#节点也进行了水平位移约束，使用2号算法(K文件已上传)，结果如下，应该是正常了：