不连续散粒体的离散单元法

建立了球体元三维离散单元法的单元接触模型。该模型把两个球体空间接触时的相互作用简化为四个弹簧阻尼器和三个滑动摩擦器。根据赫兹接触理论和非完全弹性碰撞理论，合理地确定了各接触参数。文中详细讨论了接触球体间相互作用力的分析计算方法和球体间产生接触的判断及检索方法，给出了球体元三维离散单元法的原理、算法和数据结构，并编制了相应的程序。利用该程序对球体自重下落和颗粒阻尼减振器分别进行了仿真计算，理论结果和实验结果有较好的一致性，验证了本文方法和程序的正确性。