[论文解读] Internal state of granular assemblies near random close packing

Jean-Noël Roux|arXiv (Cornell University)|May 17, 2004

Granular flow and fluidized beds参考文献 4被引用 5

一句话总结

本研究通过分子动力学模拟表明，接近随机密堆积（RCP）的颗粒堆积体在密度、配位数和取向分布上表现出独立变化，挑战了RCP对应唯一内部状态的传统观念。尽管在Φ ≈ 0.637附近密度相同，润滑压缩与振动压缩两种路径产生的结构截然不同：润滑样品具有较高的配位数（z* ≈ 5.79）和较高的弹性模量，而振动样品则具有较低的配位数（z* ≈ 4.56）且更符合实验测得的声速，揭示了弹性性质可编码隐藏的内部状态变量。

ABSTRACT

The structure of random sphere packings in mechanical equilibrium in prescribed stress states, as studied by molecular dynamics simulations, strongly depends on the assembling procedure. Frictionless packings in the limit of low pressure are devoid of dilatancy, and consequently share the same random close packing density, but exhibit fabric anisotropy related to stress anisotropy. Efficient compaction methods can be viewed as routes to circumvent the influence of friction. Simulations designed to resemble two such procedures, lubrication and vibration (or ``tapping'') show that the resulting granular structures differ, the less dense one having, remarkably, the larger coordination number. Density, coordination number and fabric can thus vary independently. Calculations of elastic moduli and comparisons with experimental results suggest that measurable elastic properties provide information on those important internal state variables.

研究动机与目标

• 挑战颗粒材料中随机密堆积（RCP）对应唯一内部状态的假设。
• 研究不同压缩工艺——润滑与振动——对接近RCP的无摩擦颗粒堆积体内部结构的影响。
• 确定弹性模量是否可作为探测配位数和取向各向异性等隐藏内部状态变量的工具。
• 将模拟得到的弹性性质与玻璃微珠堆积体在压力下的实验声速测量结果进行比较。

提出的方法

• 采用赫兹法向力和简化的Cattaneo-Mindlin-Deresiewicz关系模拟切向力的分子动力学模拟。
• 通过约束势能W = -P∑εα的最小化控制各向同性压力P，其中应变参数εα相对于参考箱体尺寸定义。
• 采用周期性边界条件，并对5个独立的4000球体样本进行平均，以确保统计稳健性。
• 模拟两种压缩路径：(A) 理想润滑（无摩擦、低压压缩）与(C) 振动敲击（敲击诱导松弛）。
• 利用标准关系式计算弹性模量：VP = √[(B + ⁴⁄₃G)/ρ]，VS = √[G/ρ]，其中B为体变模量，G为剪切模量。
• 将模拟得到的声速（VP, VS）与Domenico（1977）和Jia & Mills（2001）在10 kPa至100 MPa压力范围内的实验数据进行比较。

实验结果

研究问题

• RQ1 若不同压缩历史在相同密度下产生截然不同的内部结构，随机密堆积能否被唯一定义？
• RQ2 在接近RCP时，无摩擦润滑与振动压缩路径下，配位数与取向各向异性如何变化？
• RQ3 在致密颗粒堆积体中，弹性模量（VP, VS）在多大程度上与配位数相关，而非固体量分数？
• RQ4 能否利用压力下测得的实验声速测量值推断内部状态变量（如配位数与取向）？

主要发现

• 润滑压缩（A）得到的RCP具有Φ = 0.637 ± 0.002，z* = 5.79 ± 0.007，f₀ = 1.7%的游动颗粒，与理想RCP一致。
• 振动压缩（C）产生更致密的堆积（Φ = 0.635 ± 0.001），但配位数显著降低（z* = 4.56 ± 0.03），游动颗粒分数更高（f₀ = 13.3% ± 0.6%）。
• 尽管密度相近，润滑与振动样品在配位数和取向分布上存在显著差异，证明密度、配位数与取向可独立变化。
• 振动样品（C）的声速（VP与VS）与Domenico（1977）和Jia & Mills（2001）的实验数据更吻合，尤其在log(V)与log(P)关系的斜率上：C样品中VS斜率接近0.25，VP斜率接近0.22。
• 润滑样品的声速压力依赖性遵循幂律，指数约为1/6，与有效介质理论一致；而振动样品表现出更陡峭、非普适的标度行为。
• 润滑样品中更高的弹性模量归因于其更高的配位数，而非密度，表明配位数是弹性响应的关键决定因素。