请教各位大佬一个问题

我在用LS-DYNA做一个土中爆炸液化的案例时，总是刚开始算第一步时就报错说出现负体积，请问这是什么原因呢？
图1为底部固定约束，顶面为自由面，其余各面为法向约束，考虑孔隙水效应关键字时的网格变形情况；

图2为底部固定约束，顶面为自由面，其余各面为法向约束，不考虑孔隙水效应关键字时的网格变形情况；

图3为报错界面，两种情况报错内容一样，均为负体积。

我在建模时采用均匀的六面体网格建模，土体采用拉格朗日单元，FHWA材料模型，加了control_pore_fluid和BOUNDARY_PORE_FLUID关键字用于输出孔隙水压力；炸药和空气采用1 point ALE单元，流固耦合算法。