在球体离散元数值模拟中,颗粒的初始排列状态是影响计算效率和计算结果的重要环节.本文采用前进面几何构造算法,提出了一种基于网格搜索的球形颗粒随机排列高效算法.通过求解空间三边方程,满足了粒径设置的任意大小的颗粒依次置入前进面的外侧,并与构成前进面的三个颗粒相互接触.为获得高体积分数的颗粒簇,该算法允许颗粒改变其粒径大小.采用颗粒网格化方法可以简化前进面的搜索,并由此提高排列效率.通过计算平均配位数、体积分数和二阶结构张量的特征值,对不同粒径比下得到的立方体试样进行了分析,得到试样配位数及体积分数均随着粒径比的增大而增大,且得到的试样为各向同性.此外,空间网格的大小和初始颗粒的生成点对随机排列的效率均会产生显著的影响.最后,对非规则铁路道砟进行了精细构造及压碎模拟,发现DEM模拟得到的应力-应变曲线与试验结果基本吻合,验证了该算法得到的颗粒试样在模拟道砟裂纹起裂、扩展等过程的有效性.
球形颗粒排列 前进面法 空间网格法 结构对称性 非规则试样 离散元方法