各向同性压缩过程中圆盘体系的基本性质-计算物理.pdf

第３４卷 第２期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３４，Ｎｏ．２ ２０１７年３月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｍａｒ． ，２０１７ 文章编号：１００１⁃２４６Ｘ（２０１７）０２⁃０２４５⁃０８ 各向同性压缩过程中圆盘体系的基本性质 张兴刚，　 胡　 林 （贵州大学理学院，贵阳　 ５５００２５） 摘　 要：结合理论分析与数值模拟研究仅有线性排斥力的二维无序圆盘体系在准静态各向同性压缩过程中的基本 性质．引入平均约化重叠量、约化压强等概念，基于微观的力学分析，导出约化压强与体系的组成成分、几何统计量 及形变统计量之间的关系．基于数值模拟结果，讨论平均接触数增量、各种接触类型的百分比、平均约化重叠量、约 化压强等物理量在大范围压缩的过程中随体积分数增量的变化，给出一个形式简洁的近似约化压强公式．模拟结 果表明圆盘体系在远离临界阻塞态及接近临界阻塞态时的压缩响应不同． 关键词：颗粒物质；各向同性压缩；阻塞态；体积分数增量 中图分类号：Ｏ４１４２ 文献标志码：Ａ ０　 引言 对于以排斥性相互作用占主导的体系来说，当体系的体积分数增大到一定程度时，其体弹性模量等物理 性质一般会发生不连续变化；也就是说，体系由非阻塞态转变为阻塞态，这样的转变称为阻塞转变．近年来， 胶体、颗粒物质等体系的阻塞行为引起了许多物理学研究者的兴趣．早在 １９９８年，Ｌｉｕ 和Ｎａｇｅｌ就提出有可 ［１］ 能在统一的框架下对颗粒物质、胶体系统、玻璃态物质等不同体系中的阻塞行为进行理解 ．随后，许多人 用零温度下二维或三维球体系统研究了颗粒物质等非热力学体系的阻塞转变及其处于临界阻塞态附近的性 质［２－３］．由于单分散球体系统在接近阻塞态时容易晶格化，因此研究者通常采用双分散或多分散球体系统获 ［４］ 得无序堆积．ＯＨｅｒｎ等人 对无摩擦双分散体系的阻塞转变进行了数值模拟，发现当颗粒数Ｎ→ ∞ 时存在 一个临界体积分数ϕ ；当体积分数ϕ 由小于ϕ 向大于ϕ 变化时，平均接触数Ｚ 由０变为Ｚ ＝ ２ｄ（ｄ是维 Ｃ Ｃ Ｃ Ｃ 数），体系由非阻塞态变为阻塞态；在ϕ 附近，体系的压强Ｐ、平均接触数增量ΔＺ ＝ Ｚ －Ｚ 等参量与体积分 Ｃ Ｃ ［５］ 数增量Δϕ ＝ ϕ －ϕ 成幂律增长的关系．Ｂｅｈｒｉｎｇｅｒ小组 通过光弹性实验对二维双分散体系进行讨论，也得 Ｃ 出了类似的结果．更进一步的研究表明［６－７］ ，临界体积分数ϕ 与产生临界阻塞态的算法有关，并非一个确定 Ｃ 的值，但是ϕ 附近的幂律规律不变．最近，人们对光滑球体系统的阻塞态进行了更深入地探讨［８－９］ ；另外，还 Ｃ 对有摩擦的情况展开了研究，此时体系的阻塞行为更加复杂［１０－１３］．可以看到，前人主要在临界阻塞态附近进 行实验或模拟，理论研究工作较少；另外，研究者一般在Δϕ ＜ ０１的范围内讨论问题，此时两个颗粒间由于 相互挤压而产生的名义应变远小于百分之十．在泡沫、水凝胶等体系的实验中［１４－１５］ ，体积分数增量Δϕ 可以 大于０２，此时在接触点处产生的名义应变可大于百分之十．由于实际当中存在比较‘软’的颗粒，为了更全 面地理解颗粒体系的阻塞行为，有必要结合理论与数值模拟在较大的形变范围内考察球体系统阻塞态的性 质特征．