目的:采用数值模拟方法研究涉及大变形的土柱坍塌过程。创新点:基于临界状态的高级本构框架,利用光滑粒子流体动力学(SPH)方法模拟土柱坍塌过程,并同时研究颗粒形状对坍塌的影响。方法:采用数值仿真和离散元的方法。结论:本文提出的SIMSAND-SPH方法在模拟涉及大变形问题的土柱坍塌过程方面具有很高的效率和计算精度,可为实际工程问题提供参考依据。