颗粒尺寸分散度对颗粒体系力学和几何结构特性影响.pdfVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 18 (2013) 184501 颗粒尺寸分散度对颗粒体系力学和 几何结构特性的影响* 1 1† 2 冯旭 张国华 孙其诚 1) ( 北京科技大学物理系, 北京 100083 ) 2) ( 清华大学, 水沙科学与水利水电工程国家重点实验室, 北京 100084 ) ( 2013年5月4 日收到; 2013年6月18 日收到修改稿) 利用颗粒离散元方法, 研究了由2048 个光滑颗粒组成的体系在各向同性压缩条件下, 颗粒尺寸分散度s 对颗粒 体系力学和几何结构特性的影响. 结果表明: 配位数、剪切模量、静态结构因子以及取向序关联函数等都随分散 度的变化而变化, 而力的累积分布不受分散度的影响. 其中, 单分散体系(亦即s 0) 的静态结构因子在低波数区域 (亦即低k) 遵从幂律标度S k 02k−43 , 各分散度下取向序关联函数峰值符合e 指数变化规律g6 r ae−r6 , 其 中序关联长度6 随分散度s 增大而减小. 关键词: 颗粒物质, 尺寸分散度, 结构因子, 取向序关联 PACS: 45.70.−n, 61.43.Bn, 81.05.Rm DOI: 10.7498/aps.62.184501 线仍然以指数方式衰减. 1 引言 近年来, 关于二维体系结构因子的研究也取得 6 了很大进展. Meyer 等 发现在中间波矢区域二 近年来, 无序材料的力学和几何结构特性受到 维聚合物溶液结构因子随k 呈幂律标度S k ∼ k , 了较多关注. Silbert 和Silbert1 对比分析了三维无 且随着链长度N 的增加, 幂律指数 从 −2 变到 7 摩擦和有摩擦颗粒体系的静态结构因子S k, 发现 −3. Wen 等 发现二维颗粒链的静态结构因子 光滑颗粒体系在双对数坐标下S k 曲线在低k 区 S k ∼ k−2 , 与密集的聚合物溶液结构相似. 的线性行为, 其斜率近似为 1. Xu 和Ching2 研究 在本文中, 我们用颗粒离散元法(DEM) 生成了 了三维双分散光滑颗粒体系中, 发现静态结构因 由2048 个二维圆形颗粒、边壁压强为1000Pa 的系 子强烈的依赖颗粒粒径比, 以前发现的低k 部分的 统, 为了研究颗粒尺寸分散度s 对系统性质的影响, S k ∼ k 仅仅是颗粒的粒径比接近1 时的特殊情况. 每个分散度下随机生成 100 个位形. 通过研究体系 与此同时, 关于二维体系力学和几何结构特 的配位数、剪切模量、静态结构因子、取向序关 3 联、力的累积分布等物理量随尺寸分散度的变化 性的研究也取得了很大进展. Han 等 研究了二 维胶体晶体的取向序关联函数, 发现当处于Hex- 规律, 进一步分析了体系的无序程度对二维颗粒体 atic 相 - 液相共存区时, 系统的取向序关联函数 系力学和几何结构特性造成的影响. g6 r ∼ e−r6 , 且取向序关联长度6 随面密