Computation of stresses in jammed packings modeled with Tresca friction

技术摘要：基于 Tresca 摩擦模型的阻塞填充（Jammed Packings）应力计算

问题陈述
本文研究了刚性多边形单元（如砌体结构或致密颗粒材料）在阻塞填充状态下的应力场计算。不同于通常依赖弹性变形来惩罚相互穿透的传统离散元法（DEM），本研究关注的是细胞完全刚性且不发生弹性变形的系统。在这些“阻塞”构型中，机械平衡完全通过接触力和摩擦力来维持，而不存在显著的弹性应变。具体的挑战在于确定这些静态超定（hyperstatic）系统的内部应力分布，在这些系统中，力变量的数量超过了平衡约束的数量，导致解空间具有非唯一性。作者使用 Tresca 摩擦定律对接触物理进行建模，该定律设定了一个固定的切向滑动阈值，这与依赖应力的库仑摩擦（Coulomb friction）不同。

方法论
所提出的方法遵循一个三阶段计算框架：

1. 原问题能量最小化（Primal Energy Minimization）：
   作者构建了一个受约束的最小化问题，其中总摩擦能被最小化。能量泛函由所有内部界面上的 Tresca 摩擦功组成，定义为 $s_T |[u]_e \cdot t_e|$，其中 $s_T$ 是摩擦系数，$[u]_e$ 是位移跳跃，$t_e$ 是切向量。该最小化过程受以下约束：

   • 非穿透约束： 位移跳跃的正向分量必须非负（$[u]_e \cdot n_e \geq 0$）。
   • 约束条件： 应用 Neumann 边界条件（$g$）以防止单元发生无限滑动，因为该能量函数缺乏弹性回复力。
   • 平移不变性： 添加全局约束以固定质心。
     单元被建模为仅具有平移自由度的刚体；旋转自由度被明确忽略。这种简化是基于以下观察：在阻塞构型中，抵抗旋转（需要剪胀性）的阻力比抵抗滑动高出数个数量级。

2. 对偶极大化（Dual Maximization）：
   由于原问题寻求的是位移（在稳定的阻塞状态下位移实际上为零），作者推导出了一个对偶极大化问题。对偶变量对应于作用在内部边缘上的力。求解该对偶问题可以得到界面处的法向应力和切向应力。在关于约束项的适当假设下，原能量的极小值存在性以及对偶问题的解得以确立。对偶公式有效地从可能的应力状态的超定流形中选择了一个特定的解。

3. 应力重构（Stress Reconstruction）：
   对偶问题的输出提供了边缘处的离散应力值。为了在每个单元内重构连续应力场 $\sigma$，作者采用了最低阶 Raviart–Thomas ($RT_0$) 有限元。这些单元是 $H(\text{div})$ 相容的，这意味着它们的自由度对应于边缘重心处的法向应力。重构过程涉及在每个单元内求解一个局部混合问题，以找到应力张量 $\sigma_c$，使得：

   • $\text{div}(\sigma_c) = 0$（在不存在体积载荷的情况下满足平衡）。
   • 边界上的法向应力与从对偶解中获得的值相匹配。
   • 应力张量的对称性（$\sigma = \sigma^T$）通过拉格朗日乘子进行弱实现。这是必要的，因为底层的离散模型忽略了力矩平衡方程，这意味着重构的应力可能具有代表残余转矩的反对称分量。

核心贡献

• 刚性摩擦模型： 本文引入了一种基于刚性摩擦（Tresca 准则）的模型，该模型不依赖弹性参数，从而区别于标准的 DEM 方法。
• 处理超定性： 该工作明确解决了阻塞填充的静态不确定性问题。它证明了优化方法通过最小化摩擦能量泛函，从自应力状态空间中选择了一个唯一的解。
• 一致的应力重构： 提出了一种使用 $RT_0$ 元进行应力重构的新技术，以从离散界面力生成一致的应力场，确保每个单元内的局部力平衡。
• 忽略旋转： 作者提出了一个简化的模型，忽略了旋转自由度，并认为这在物理上是合理的，因为在受压填充中，平移滑动占主导地位，同时也承认这会导致非对称应力张量。

数值结果
该方法在各种二维多边形网格（具体为 Voronoi 镶嵌）上进行了测试，其发现如下：

• 验证： 在均匀边界条件下，重构的应力场与预期的均匀应力状态相匹配，证实了该方法的一致性。
• 剪切测试： 在剪切加载场景下，重构的应力张量表现出不对称性（$\sigma_{12} \neq \sigma_{21}$）。作者将其归因于离散模型中缺乏强化的力矩平衡，并指出反对称部分解释了残余转矩。
• 均匀压缩： 尽管存在均匀压缩边界条件，求解器仍产生了非均匀应力场。作者将其解释为系统超定性的结果；优化过程从零空间中选择了一个特定的能量最小化状态，这自然导致了“力链”（高应力稀疏网络）的形成以及低应力或零应力区域的出现。
• 压缩中的拉伸： 结果显示，即使在全局压缩下，某些单元也出现了拉伸（正主应力）。作者认为这是一种物理现象，源于摩擦连锁效应和桥接效应，而非数值伪影。

意义与主张
论文声称，其提出的方法通过利用刚性摩擦和能量最小化的物理特性，提供了一种鲁棒且一致的方法来计算阻塞填充中的应力。作者认为，与基于弹性的 DEM 相比，其方法能更好地捕捉阻塞系统的底层物理，特别是在弹性变形可以忽略不计的机制中。

研究结果的意义在于证明了最小化摩擦能（一种 $L_1$ 型目标）自然会导致稀疏的力网络，这与实验中观察到的颗粒材料中的力链现象相吻合。作者谦虚地承认了其局限性：忽略旋转平衡导致了非对称应力张量，且使用 Tresca 摩擦（固定阈值）简化了问题，相比之下，库仑摩擦需要更复杂的非凸拟变不等式（Quasi-Variational Inequalities）。作者建议，未来的工作可以扩展该模型，以包含旋转自由度和更真实的摩擦定律，从而实现完全的机械一致性。