[论文解读] Entropic origin of stress correlations in granular materials

Gregg Lois, Junlong Zhang+5 more|arXiv (Cornell University)|May 30, 2009

Granular flow and fluidized beds参考文献 2被引用 3

一句话总结

本文提出了一种基于机械稳定性和接触力正性而非能量最小化的场论性、熵基框架，用于描述近jamming转变处颗粒材料中的应力涨落。该理论预测，在剪切作用下由于正性约束，应力相关性呈现各向异性，理论、模拟与实验结果定量吻合，解释了力链网络中的剪切诱导各向异性。

ABSTRACT

We study the response of granular materials to external stress using experiment, simulation, and theory. We derive an entropic, Ginzburg-Landau functional that enforces mechanical stability and positivity of contact forces. In this framework, the elastic moduli depend only on the applied stress. A combination of this feature and the positivity constraint leads to stress correlations whose shape and magnitude are extremely sensitive to the applied stress. The predictions from the theory describe the stress correlations for both simulations and experiments semiquantitatively.

研究动机与目标

• 开发一种非平衡统计力学框架，用于描述无热颗粒系统中的应力涨落，其中基于能量的方法失效。
• 解释在剪切作用下实验与模拟中观测到的各向异性应力相关性（在各向压缩下不存在）的起源。
• 通过机械稳定构型的熵，建立受 jammed 颗粒堆积制备历史与其宏观力学响应之间的联系。
• 通过摩擦less 与摩擦性圆盘的模拟以及受控压缩与剪切下光弹性颗粒的实验，验证该理论。

提出的方法

• 基于力与力矩平衡，构建类似 Ginzburg-Landau 的熵泛函，使用 Airy 应力函数 Ψ 表示粗粒化的应力场。
• 将应力矩张量 Σ 定义为广延的、拓扑不变量，并推导出一个场论，其中涨落 ψ 由仅依赖于 ψ 的二阶导数的拉格朗日量 LΣ[ψ] 所支配。
• 从 ψ 上的路径积分构造配分函数 Z(Σ) = exp[S(Σ)]，通过泛函平均生成相关函数。
• 在傅里叶空间中推导压力相关函数 S(q)，显示其依赖于应力比 τ/Γ 与压缩方向的各向异性、偶极子形式。
• 利用接触力的正性约束以确保物理可实现性，并解释沿压缩方向涨落的抑制。
• 将理论预测与无摩擦及有摩擦圆盘的模拟结果以及 2D 光弹性颗粒在各向压缩与纯剪切下的实验数据进行比较。

实验结果

研究问题

• RQ1 应力相关性在颗粒材料中如何依赖于制备工艺，如各向压缩与纯剪切？
• RQ2 在剪切颗粒堆积中观测到的各向异性应力相关性模式的起源是什么？
• RQ3 能否基于机械稳定性与力的正性，通过场论性、熵基方法解释观测到的应力涨落，而无需依赖能量最小化？
• RQ4 为何应力相关性在 q 空间中表现出方向依赖性，特别是在小波矢量极限下？

主要发现

• 该理论预测，在纯剪切下，应力相关性在 q 空间中呈现偶极子、各向异性的形式，沿更高压缩方向的关联更强。
• 相关函数 S(q) 的低 q 极限依赖于接近 q=0 的方向，满足 S(qy=0) > S(qx=0)，与模拟和实验结果定量一致。
• 当 τ/Γ 较小时，比值 S(qy=0)/S(qx=0) 趋近于 1 + 4(τ/Γ)，在 τ/Γ ≈ 0.51 的实验数据中，该比值接近 3，与理论预测相符。
• 应力相关性的各向异性源于接触力的正性约束，该约束抑制了沿压缩方向的涨落。
• 理论预测的压力相关函数与模拟结果（有摩擦与无摩擦）以及 2D 光弹性圆盘的实验测量结果均表现出良好的半定量一致。
• 该框架成功解释了接近 jamming 时颗粒材料的历史依赖性力学响应，将制备工艺与宏观应力相关性结构联系起来。