Analytic Nonlinear Theory of Shear Banding in Amorphous Solids

本文提出了一种针对非热剪切非晶固体中剪切带形成的解析非线性理论，通过推导计及塑性诱导偶极子屏蔽效应的方程，阐明了失稳机制，预测了剪切带宽度，并确定了失效的临界应力阈值。

要点

想象一块玻璃或一堆沙子。在物理学世界中，这些被称为“非晶态固体”。与原子排列成完美行列的晶体（如钻石）不同，这些材料中的原子是随机混乱堆积的，就像一场没有固定座位的音乐会上的观众人群。

长期以来，科学家们试图用与完美晶体相同的规则来预测这些材料如何断裂或变形。但这些规则失败了。当你推挤玻璃或沙子时，它不会仅仅发生弯曲；它会突然断裂，或形成一条狭窄、尖锐的损伤线，称为剪切带。这就像挡风玻璃上出现裂纹，但并非单一线条，而是一个材料相对于自身发生滑移的区域。

Avanish Kumar 和 Itamar Procaccia 的这篇论文提出了一种新的数学“配方”，用于精确预测这些剪切带如何以及为何形成，以及它们的外观。以下是用通俗语言进行的分解：

1. 问题：“隐藏”的混乱

当你推挤完美晶体时，它会平滑地拉伸。但当你推挤非晶态固体时，内部会发生微小、混乱的重排。作者将这些称为**“塑性事件”**。

• 类比：想象一个拥挤的房间。如果你推挤人群，人们不会仅沿直线移动；他们会相互碰撞、向侧面挪动，并形成微小的漩涡状运动。在论文中，这些漩涡被称为**“四极子”**（四种方向的运动形状）。
• 旧理论：先前的理论将这些漩涡视为均匀分布，就像糖溶解在茶水中一样。这在微小推挤时有效，但无法解释剪切带的突然、剧烈形成。
• 新见解：作者意识到，当材料受到应力时，这些漩涡不再均匀分布。它们开始聚集，形成**“偶极子”（两点力），其作用类似于屏蔽电荷**。
  • 隐喻：将这些偶极子想象成一群撑着伞的人。如果它们均匀分布，雨水（应力）会均匀地打在每个人身上。但如果它们聚集在一起，就会形成一种“盾牌”或“屏障”，在某些地方阻挡雨水，而在其他地方让雨水倾泻而下。这种屏蔽效应产生了一个特定的“长度尺度”——即损伤区域的自然宽度。

2. 重大突破：非线性数学

该论文认为，要理解剪切带，不能使用简单的直线数学（线性方程）。你需要非线性数学。

• 类比：想象驾驶汽车。在低速时，如果你稍微转动方向盘，车子就会稍微转向（线性）。但在高速时，方向盘的微小转动可能导致车辆失控旋转（非线性）。
• 作者推导出一组新方程，以解释材料的这种“高速”行为。他们包含了两个主要的非线性效应：
  1. 材料变形时形状的变化（应变 - 位移关系）。
  2. 当“漩涡”变得拥挤时，它们彼此之间的相互作用（偶极子相互作用）。

3. 结果：预测“裂纹”

通过求解这些复杂方程，作者找到了一种预测剪切带轮廓的方法。

• “延性”（较软）情况：在稍微更具灵活性的材料中，剪切带宽阔且平滑。
  • 隐喻：就像缓慢、平缓的斜坡。材料在广阔区域上逐渐滑移。数学预测这种形状看起来像双曲正切（tanh）曲线——一个平滑的 S 形。
• “脆性”（较硬）情况：在非常刚硬的材料中，剪切带极其尖锐且狭窄。
  • 隐喻：就像悬崖边缘。材料在一侧保持静止，而在另一侧瞬间滑移。数学表明，在这种情况下，带的“核心”行为与边缘不同，从而产生非常尖锐的过渡。

4. “不稳定性”开关

该论文还解释了何时会发生这种情况。

• 类比：想象将铅笔笔尖朝上平衡。只要风很小，它就能站立。但在特定的临界风速下，它会变得不稳定并倒下。
• 作者计算出了材料失去稳定性的确切“临界应力”（即风速）。他们发现，当某个特定数学值（“海森矩阵”的特征值，这只是一个 fancy 的说法，指稳定性计算器）降至零时，就会发生这种情况。
• 一旦达到这一点，材料就无法再均匀地保持其形状，剪切带便“突然”形成。

5. 为什么这很重要（根据论文）

先前的理论可以说“剪切带会形成”，但如果没有猜测或使用计算机模拟，它们无法告诉你它会有多宽或形状会是什么样。

• 这篇论文提供了一种解析理论，意味着它给出了一个直接公式。
• 它解释了剪切带的宽度是由材料的刚度与内部漩涡的“屏蔽”效应之间的竞争决定的。
• 它基于这些方程的数学，区分了脆性（尖锐、突然的断裂）和延性（缓慢、宽阔的滑移）材料。

总结

简而言之，作者建立了一个新的数学模型，将非晶态固体（如玻璃或沙子）视为复杂的移动粒子群，而非简单的弹簧。通过考虑这些粒子如何相互“屏蔽”彼此的运动，以及它们在应力下如何表现出非线性行为，他们推导出了一个公式，能够精确预测材料何时会断裂，以及 resulting 的“裂纹”（剪切带）会是什么样子——从平滑的滑移到尖锐的断裂。

技术摘要

技术摘要：非晶固体剪切带形成的解析非线性理论

问题陈述
非晶固体（如玻璃、颗粒介质）在无热准静态剪切下的力学响应与晶体固体存在根本差异。晶体在屈服前仅表现出纯弹性响应，而非晶固体在任何应变水平下都会发生塑性变形，通常表现为四极子形式的“埃谢尔比（Eshelby）夹杂”。这些材料中的一个关键不稳定性是剪切带，即变形局域化形成狭窄的带状区域。这种现象涵盖了从跨越平面的位移发生尖锐、类脆性跳跃，到更平滑、类延性的局域化。

现有理论面临显著局限：

1. 线性理论：先前的解析方法依赖于应变与位移之间的线性近似以及拉格朗日量的二次展开。虽然这些方法成功描述了屏蔽效应和重整化弹性模量，但它们无法预测剪切带的具体轮廓或失稳发生的临界应力。
2. 弹塑性模型：标准的介观尺度模型能够捕捉雪崩统计和应力重分布，但通常在空间上是局域的。它们缺乏内禀的连续介质长度尺度，这意味着剪切带的宽度由数值正则化或块体尺寸决定，而非物理材料属性。因此，它们无法唯一确定内部应变轮廓。

作者旨在开发一种非线性解析理论，推导位移场方程，预测失稳所需的临界累积应力，提供剪切带轮廓的解析解，并根据材料属性区分延性和脆性特征。

方法论
作者构建了一种基于拉格朗日表述的变分理论，其中纳入了两个此前被忽视的关键非线性项：

1. 非线性应变 - 位移关系：超越线性关系 $u_{\alpha\beta} = \frac{1}{2}(\partial_\alpha d_\beta + \partial_\beta d_\alpha)$，该理论采用了完整的几何定义：
   $$u_{\alpha\beta} = \frac{1}{2} [\partial_\alpha d_\beta + \partial_\beta d_\alpha] + \frac{1}{2}\partial_\beta d_\gamma \partial_\alpha d_\gamma$$
2. 高阶偶极子展开：虽然塑性事件诱导四极子，但四极子密度的梯度会产生有效的偶极子。作者将拉格朗日量展开以包含偶极子场（$P^\alpha$）的四次项，这对于捕捉失稳的非线性饱和是必要的。

推导过程遵循以下步骤：

• 非线性方程的推导：从包含背景应力（$\Sigma_{\alpha\beta}$）、弹性能以及四极子/偶极子相互作用的能量泛函出发，作者针对位移场（$d$）和四极子场（$Q$）最小化能量。这导出了位移场的非线性偏微分方程（公式 32）。
• 投影至标量振幅方程：为求解复杂的矢量方程，作者采用基于系统海森矩阵（Hessian）软模的策略。他们论证道，在失稳附近，位移场可以投影到沿带法线 $\xi = \mathbf{m} \cdot \mathbf{x}$ 变化的标量函数 $f(\xi)$ 上。这将问题简化为一维“振幅方程”（公式 61）。
• 能量泛函与海森矩阵分析：构建一个能量泛函，其欧拉 - 拉格朗日方程与推导出的振幅方程相匹配。分析该泛函的海森矩阵以识别最低本征值消失的临界点，从而标志失稳的开始。
• 解析解：在两个极限下求解非线性振幅方程：
  • 延性极限：斜率较小，忽略高阶梯度项。
  • 脆性极限：带芯由四次梯度项主导，需要匹配渐近展开（内解和外解）。

主要贡献与结果

1. 非线性振幅方程的推导：
   核心结果是公式 (61)，即剪切带轮廓 $f(\xi)$ 的非线性微分方程：
   $$(\mu + \Sigma(m)) f''(\xi) + \frac{3}{2}(\lambda + 2\mu) (f'(\xi))^2 f''(\xi) = -A f(\xi) + B f^3(\xi)$$
   其中，$\Sigma(m)$ 代表背景应力的投影，$A$ 和 $B$ 是由耦合张量和屏蔽参数导出的参数。

2. 临界应力的预测：
   通过分析能量泛函的线性化海森矩阵，作者推导出了失稳开始的判据。临界条件发生在海森矩阵的最低本征值消失时：
   $$(\mu + \Sigma(m)) \left(\frac{2\pi}{L}\right)^2 = A$$
   这提供了触发剪切带所需的临界累积应力的预测，将其直接与屏蔽参数 $A$（与 $\tilde{\kappa}_e^2$ 相关）联系起来。

3. 延性与脆性带的解析轮廓：

   • 延性解：在小斜率极限下，解为双曲正切（tanh）轮廓：
     $$f(\xi) = f_0 \tanh\left(\frac{\xi - \xi_0}{\ell}\right)$$
     宽度 $\ell$ 由屏蔽参数和弹性模量控制：$\ell \sim \sqrt{(\mu + \Sigma)/A}$。
   • 脆性解：在大斜率极限下，带芯由四次项主导，导致正弦波状带芯与双曲正切尾部相匹配。带芯宽度 $L_c$ 由三次偶极子项的系数（$G$）决定，不同于屏蔽长度。这解释了脆性剪切带更尖锐、更局域化的特性。

4. 内禀长度尺度：
   该理论证明，剪切带形成需要由拉格朗日量中的耦合张量生成的内禀长度尺度。该尺度源于破坏性应力项与正则化梯度项（屏蔽）之间的竞争。与宽度任意的弹塑性模型不同，这里的宽度是由材料参数（$\mu, \lambda, A, B, G$）决定的物理属性。

5. 非线性海森矩阵的零模：
   作者证明，剪切带解的导数 $f'(\xi)$ 是非线性海森矩阵的零模（本征值为零的本征函数）。这将扭结（kink）解的平移不变性与稳定性分析直接联系起来。

意义与主张
本文声称提供了首个能够同时预测剪切带临界应力和带详细空间轮廓的解析理论，并能区分延性与脆性行为。

• 解决先前的失败：作者认为，先前的理论（如 Rudnicki 和 Rice 的理论）未能提供解，是因为它们缺乏使失稳饱和所需的非线性项。振幅方程中的三次项（$B f^3$）抑制了线性不稳定性，而非线性梯度项控制了带芯形状。
• 局域化的物理选择：该理论提出，剪切局域化是一个“真正的屏蔽连续介质对象”。塑性偶极子对长程弹性相互作用（埃谢尔比场）的屏蔽至关重要；如果没有这种屏蔽，变形将全局扩散。该理论自然地生成了有限宽度的带，而无需人为的正则化。
• 与弹塑性模型的区别：本文强调，标准的弹塑性模型在预测剪切带轮廓方面是不完整的，因为它们缺乏内禀长度尺度。相比之下，该理论通过耦合参数将长度尺度直接构建到拉格朗日量中，使带宽度成为一个可预测的、依赖于材料的量。
• 对参数的谦逊态度：作者承认，虽然解的函数形式是解析推导的，但特定参数的数值（如决定脆性极限内长度尺度的 $B$）尚未被测量。他们提出，未来需要对经典玻璃形成模型进行数值模拟，以确定这些参数并充分验证该理论。

总之，本文提供了一个严格的变分框架，统一了塑性屏蔽、非线性弹性和不稳定性理论的概念，以解释非晶固体中剪切带的产生、临界性和形态。