小天体探测是当今太阳系探测的一大热点,对小行星演化的研究有助于人们了解太阳系的起源.演化研究的一项重要内容是小天体结构的演化,即小天体在多种力学机制作用下自身形状与结构的演化.在小天体碎石堆结构的假设之下,一种比较常见的模拟小天体结构演化的方法是离散元仿真算法(Discrete Element Method, DEM),目前国内外有数个团队开发了相关算法软件.介绍了团队开发的《基于DEM仿真算法的多粒子系统模拟软件》的基础理论、实现方法以及加速算法,并使用二体接触模型、声速传播仿真、小行星内部压强、小行星旋转稳定性仿真算例验证了算法的可靠性.