如何使用lsdyna仿真泥石流？

大佬们，我的毕设需要用lsdyna来仿真排障器撞击泥石流。

泥石流是由泥土和石块组成。

泥土用SPH表示。石块用DEM表示，也可以用有限元表示。排障器用有限元表示。

我需要通过耦合光滑粒子流体动力学(SPH)、离散元方法(DEM)和有限元方法(FEM)来考虑复杂的流体-颗粒-结构相互作用(FPSI)。

我现在sph和fem相互作用会做，可是怎么在泥土中混入石块不会做，求大佬们指教，有偿。