[论文解读] The nature of quasistatic deformation in granular materials

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

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

一句话总结

该论文通过单分散球体颗粒系的离散元模拟，揭示了准静态颗粒变形中存在两种截然不同的流变行为：I型（通过接触刚度实现的弹性变形）和II型（接触网络的塑性重排）。关键贡献在于表明：在小应变和载荷反向时，I型行为占主导，此时应变与接触刚度成反比；而在大应变下，II型行为出现，表现出显著的力场不确定性与结构非均匀性。

ABSTRACT

Strain in granular materials in quasistatic conditions under varying stress originate in (I) contact deformation and (II) rearrangements of the contact network. Depending on sample history and applied load, either mechanism might dominate. One may thus define rheological regimes I and II accordingly. Their properties are presented and illustrated here with discrete numerical simulation results on sphere packings. Understanding the microscopic physical origin of strain enables one to clarify such issues as the existence of macroscopic elasticity, the approach to stress-strain relations in the large system limit and the sensitivity to noise.

研究动机与目标

• 阐明准静态颗粒变形中应变的微观起源，区分接触变形与网络重排的贡献。
• 定义并表征两种截然不同的流变行为：I型（以弹性接触变形为主导）与II型（以接触网络重排为主导）。
• 研究宏观力学行为（包括弹性与塑性）如何依赖于接触网络结构与加载历史。
• 探讨力场不确定性与配位数在调控两种行为之间转变中的作用。
• 评估载荷反向与扰动对应变不可逆性及行为主导性的影响力。

提出的方法

• 采用Hertz-Mindlin接触定律对单分散球体颗粒系在三轴压缩条件下的离散元方法（DEM）模拟。
• 使用无量纲惯性参数 $ I = \dot{\epsilon}_{a} \sqrt{m/aP} $ 衡量偏离准静态条件的程度。
• 引入法向黏性阻尼，其中阻尼系数 $ \zeta $ 以临界阻尼的分数形式表示，以保证数值稳定性。
• 定义无量纲接触刚度参数 $ \kappa = (E/P)^{2/3} $，用于量化弹性变形量。
• 将完整模拟与固定接触列表（无新接触形成）的模拟进行对比，以隔离I型行为。
• 分析不同样品制备工艺（A与C）下的应力-应变响应、弹性模量及力网络的不确定性。

实验结果

研究问题

• RQ1 是什么因素决定了颗粒材料中的应变主要来源于弹性接触变形，还是接触网络的重排？
• RQ2 I型或II型行为的主导性如何依赖于样品制备方式、配位数及加载路径？
• RQ3 在I型行为中，应变与接触刚度的标度关系如何？该关系在II型行为中是否有所不同？
• RQ4 载荷反向如何影响颗粒材料中不可逆（塑性）应变的大小与性质？
• RQ5 微观特征（如力场不确定性与配位数）在多大程度上可预测两种行为之间的转变？

主要发现

• I型行为以弹性接触变形为主导，其特征为应变与接触刚度 $ E_{m} $ 成反比，该行为在初始加载及载荷反向阶段均可观察到。
• 在I型行为中，应力-应变曲线的初始斜率与宏观杨氏模量 $ E_m $ 一致，且 $ E_m $ 随 $ \kappa $ 变化，应变与 $ 1/E_m $ 呈线性比例关系。
• 对于载荷反向过程，初始卸载斜率等于 $ E_m $，且不可逆应变出现在约 $ \sim 0.2P $ 的区间内，证实I型应变并非完全弹性。
• 在II型行为中，应变幅值对接触刚度不敏感，而依赖于颗粒间间隙；该行为在大应变下占主导，表现为网络的持续性破坏与重构。
• 行为转变由力场不确定性决定：I型行为发生在力网络高度超稳定（力解空间极大）时，而II型行为在摩擦滑移使该解空间缩小至自由度的约10%时出现。
• 采用固定接触列表的模拟（无新接触形成）在I型行为区域与完整模拟结果无法区分，证实该区域接触网络变化可忽略。