用Monte Carlo方法模拟砂土的自然堆积过程

用MonteCarlo随机模拟方法模拟砂土颗粒在重力作用下,在圆柱与长方体容器中的自然堆积过程。首先用参考网格法生成一个砂土颗粒的松散结构,松散结构中粒子与粒子、粒子与边界间不存在任何接触;然后启动MonteCarlo随机模拟算法,即给处于松散结构中的每个粒子施加随机位移,得到新构形,如果新构形中粒子间或粒子与边界间发生重叠,则放弃这个构形;如果没有重叠存在,则判别粒子体系的势能变化,运用Metropolis准则来判别这个构形是否被接受,重复这个过程可以得到砂土的密集堆积结构。采用Schinner建议的接触发现算法判别粒子间是否存在重叠,同时详细介绍了粒子与边界间的接触发现算法。模拟结果表明,用MonteCarlo方法模拟砂土的自然堆积结构是非常有效的,可以为砂土的流动、压实等的数值模拟工作提供初始构形。