[论文解读] Hydrodynamic Theory of Granular Solids: Permanent, Transient and Granular Elasticity

Yimin Jiang, Mario Liu|ArXiv.org|Jun 10, 2007

Geotechnical and Geomechanical Engineering参考文献 3被引用 3

一句话总结

本文通过将永久弹性和瞬态弹性统一到基于可逆与不可逆热力学的一致框架中，发展了一套颗粒固体的流体动力学理论。该理论推导出一组描述剪切或敲击下弹性能量弛豫的微分方程，关键成果是推导出一个经过实验数据验证的热力学一致的颗粒弹性能量表达式。

ABSTRACT

Although fully elastic when static, granular media become transiently elastic when being slowly sheared -- during which both the elastic energy and stress relax. Starting from this observation, we cogently derive the framework for granular hydrodynamics, a set of differential equations consistent with general principles of physics, especially reversible and irreversible thermodynamics. In addition, an expression for the granular elastic energy is reviewed and further discussed.

研究动机与目标

• 建立一个统一的颗粒固体流体动力学框架，以同时考虑永久弹性和瞬态弹性。
• 推导一组与热力学原理一致的微分方程，包括可逆与不可逆过程。
• 利用来自颗粒静力学的实验数据，提供一个物理解释清晰的颗粒弹性能量表达式。
• 通过将流体动力学理论与能量及输运系数的具体函数形式关联，实现预测建模。

提出的方法

• 从热力学原理推导颗粒流体动力学，将弹性能量视为应变场 $ u_{ij} $、熵 $ s $、质量密度 $ \rho $ 和动量 $ g_i $ 的函数。
• 引入颗粒温度 $ T_g $ 作为独立变量，以模拟剪切流中内部能量的产生，类似于粘性加热。
• 采用应变分解 $ \varepsilon_{ij} = u_{ij} + p_{ij} $，其中 $ u_{ij} $ 为弹性部分，$ p_{ij $ 为塑性部分，以分离可逆与不可逆贡献。
• 对弹性能量 $ w(u_{ij}) $ 应用凸性条件，以确保热力学稳定性，从而对 $ w $ 的二阶导数施加约束。
• 通过条件 $ \vec{m}_1 \parallel \partial g / \partial \vec{\sigma} $ 推导关联流动法则，确保塑性流动垂直于屈服面。
• 利用实验数据验证能量形式，表明其与观测到的颗粒静力学及屈服行为一致。

实验结果

研究问题

• RQ1 如何从热力学的一般原理出发，一致地推导出颗粒固体的流体动力学理论？
• RQ2 颗粒温度 $ T_g $ 在慢速剪切或敲击过程中实现瞬态弹性时起到什么作用？
• RQ3 弹性能量 $ w(u_{ij}) $ 如何依赖于应变不变量，且何种约束可确保热力学稳定性？
• RQ4 颗粒弹性能量的何种函数形式最符合静态颗粒系统中的实验数据？
• RQ5 所推导的理论与已知模型（如超塑性模型）有何关联，对预测建模有何启示？

主要发现

• 该流体动力学理论通过一组从热力学原理推导出的一致微分方程，统一了永久弹性和瞬态弹性，涵盖应力平衡与能量演化。
• 弹性能量 $ w(u_{ij}) $ 被证明是应变的凸函数，其稳定性由海塞矩阵的正定性及对其二阶导数的特定不等式保证。
• 屈服条件 $ \mathcal{B}/\mathcal{A} > 1 $ 作为热力学稳定性的必要条件被推导出，将材料参数与力学响应联系起来。
• 关联流动法则被推导为 $ \vec{m}_1 \parallel \partial g / \partial \vec{\sigma} $，确保塑性流动垂直于由 $ g=0 $ 定义的屈服面。
• 能量表达式 $ w = \mathcal{A} \Delta^2 + \mathcal{B} u_s^2 $ 经实验数据验证，其中 $ \mathcal{A} > 0 $，$ \mathcal{B} > 0 $，且 $ \mathcal{B}/\mathcal{A} > 1 $，确保稳定性。
• 条件 $ \left(\partial P / \partial \Delta\right)_{\sigma_s} > 0 $ 被识别为关键的热力学稳定性判据，等价于能量在应变变量上的凸性。