碎石堆小行星自旋演化过程的离散元数值模拟

近年的探测研究表明，由引力和范德华粘附力聚合成碎石堆结构可能是小行星的一种重要形成方式。为探究这类小行星的演化机制，文章基于离散元软件MatDEM进行二次开发，实现了单元的万有引力计算。建立直径为1.2 km“碎石堆”小行星数值模型，通过数值模拟分析了不同粘附强度和体积密度下，球形聚集体的旋转重塑和破坏过程。结果表明：球形聚集体受到其角速度初始增量的影响，以启动聚集体的重塑过程，导致其发生变形或者破坏；粘附强度越小，模型越容易被破坏，破坏前小行星会通过向扁圆状椭球体变形来保持自身构型的稳定性；对于粘附强度一定的小行星模型，体积密度的增大会有效增加其结构的稳定性，可以延后小行星破坏行为发生的时间，且颗粒脱落的数量随着密度的增大而减少；采用离散元分析可以较好地模拟小行星自旋演化过程。进一步与真实小行星数据相结合，有利于深入、全面地探索小行星自旋演化破坏机理。