[论文解读] Elasticity of sphere packings: pressure and initial state dependence

Ivana Agnolin, Jean-Noël Roux|ArXiv.org|Jan 15, 2009

Material Dynamics and Properties参考文献 11被引用 6

一句话总结

本研究探讨了初始制备工艺（特别是振动式与润滑式组装）对各向同性球体密堆积体弹性模量及内部结构的影响。通过施加受控压力循环的数值模拟，发现弹性模量主要由配位数决定，而不仅由密度决定；在配位数较低的体系中，剪切模量表现出异常的 P^{1/2} 依赖关系，挑战了标准均质化模型。

ABSTRACT

Elastic properties and internal states of isotropic sphere packings are studied by numerical simulations. Several numerical protocols to assemble dense configurations are compared. One, which imitates experiments with lubricated contacts, produces well coordinated states, while another, mimicking the effect of vibrations, results, for the same density, in a much smaller coordination number z, as small as in much looser systems. Upon varying the confining pressure P, simulations show a very nearly reversible variation of density, while z is irreversibly changed in a pressure cycle. Elastic moduli are shown to be mainly related to the coordination number. Their P dependence notably departs from predictions of simple homogenization approaches in the case of the shear moduli of poorly coordinated systems.

研究动机与目标

• 理解不同的数值组装工艺如何影响致密球体堆积体的弹性性质和内部结构。
• 研究配位数 z 和初始状态（如润滑与振动）在不同压力下对弹性模量的影响。
• 评估经典均质化模型在预测颗粒材料弹性模量压力依赖性方面的有效性。
• 评估弹性模量是否可作为内部堆积状态（尤其是与实验数据对比）的可靠指标。
• 确定为何完全润滑的堆积体模拟结果会高估刚度，而与实验结果相比，振动式工艺是否能提供更好的一致性。

提出的方法

• 采用周期性边界条件和各向同性应力控制（Parrinello-Rahman 方法）的分子动力学数值模拟。
• 接触力采用赫兹弹性模型（F_N ∝ h^{3/2}），切向力采用改进的 Cattaneo-Mindlin 理论并引入库仑摩擦（μ = 0.3）。
• 四种不同的制备工艺：A（润滑，无摩擦）、B（轻微润滑，μ₀ = 0.02）、C（振动，膨胀后混合）、D（直接压缩，μ = 0.3）。
• 施加从 10 kPa 到 100 MPa 再返回的压力循环，保持准静态条件（应变率 < 10^{-4} × √(m/aP)）。
• 通过接触网络在平衡状态下的刚度矩阵计算弹性模量，使用 Voigt 和 Reuss 估计法推导上下限。
• 对 5 个独立样本进行统计分析，误差棒以一个均方根偏差表示。

实验结果

研究问题

• RQ1 不同的初始组装工艺（润滑与振动）如何影响相同密度下的配位数 z* 和堆积分数 Φ？
• RQ2 在配位数较低的体系中，剪切模量 G 如何随压力变化？其是否偏离简单均质化模型预测的 P^{1/3} 标度？
• RQ3 为何润滑堆积体（类型 A）的模拟结果会高估刚度，尽管其与 MPa 范围内的实验数据吻合？
• RQ4 弹性模量能否作为颗粒堆积体内部状态（如接触网络质量）的可靠指标？
• RQ5 在压力循环过程中，配位数 z* 的不可逆降低与初始制备方法之间有何关联？

主要发现

• 通过类振动工艺制备的样品（C）达到高密度（Φ = 0.635 ± 0.002），但配位数较低（z* = 4.56 ± 0.03），且存在 13.3% 的游动粒子，而润滑样品则为 z* = 6.07，游动粒子仅占 1.3%。
• 所有工艺下，密度 Φ 随压力可逆变化，但在初始配位数较高的体系（如 A 和 B）中，z* 在压力循环过程中出现不可逆降低。
• 弹性模量主要由配位数 z* 决定；在配位数较低的体系中（C 和 D），剪切模量 G 显示出 P^{1/2} 依赖关系，显著偏离 P^{1/3} 的预测。
• Voigt 上限在低 z 体系中高估了 G（例如，G^Voigt = 1.34 × B^Voigt），表明存在强烈的非仿射位移和应变涨落。
• 体积模量 B 的 Reuss 下限与真实值相差不超过约 1.15 倍，表明 B 在所有体系中均近似遵循 P^{1/3} 的标度关系。
• 类振动结构（C）的超声波速度比润滑结构（A/B）更接近实验观测值，表明其更准确地反映了真实实验室中的颗粒堆积状态。