請教各位大佬一個問題

我在用LS-DYNA做一個土中爆炸液化的案例時，總是剛開始算第一步時就報錯說出現負體積，請問這是什么原因呢？
圖1為底部固定約束，頂面為自由面，其余各面為法向約束，考慮孔隙水效應關鍵字時的網格變形情況；

圖2為底部固定約束，頂面為自由面，其余各面為法向約束，不考慮孔隙水效應關鍵字時的網格變形情況；

圖3為報錯界面，兩種情況報錯內容一樣，均為負體積。

我在建模時采用均勻的六面體網格建模，土體采用拉格朗日單元，FHWA材料模型，加了control_pore_fluid和BOUNDARY_PORE_FLUID關鍵字用于輸出孔隙水壓力；炸藥和空氣采用1 point ALE單元，流固耦合算法。