🔬 condensed matter

Analytic theory of shear localization in amorphous solids confined by Couette geometry

本文提出了一种解析理论，用以解释二维非晶态固体在准静态库埃特（Couette）旋转下的剪切局部化现象，其中应力积累触发塑性雪崩，导致相邻区域向相反方向旋转，从而实现了对所得位移场的先验预测。

原作者： Yang Fu, Yuliang Jin, Itamar Procaccia

发布于 2026-06-26

📖 1 分钟阅读☕ 轻松阅读

原作者： Yang Fu, Yuliang Jin, Itamar Procaccia

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象你有一块玻璃。虽然它看起来像固体，但它的内部结构并不像普通晶体那样整齐排列。科学家们称这种材料为“非晶态固体”。这意味着构成它的微粒并没有按照规则的晶格排列，而是杂乱无章地堆积在一起。常见的例子包括窗户玻璃、像沙子这样的颗粒材料，以及金属玻璃。

现在，想象这种材料被两个巨大的、看不见的圆环困住了：中间是一个小圆环，外面是一个大圆环。

这篇论文中的科学家们决定慢慢旋转内圈，就像转动罐子的盖子一样，同时保持外圈不动。他们想看看里面的材料会如何反应。

以下是他们发现的故事，用简单的语言进行了解释：

问题所在：“意外”的扭转

如果你挤压一根正常的、具有完美弹性的橡皮筋，它会平滑且可预测地拉伸。如果你旋转这个黏稠材料的内圈，你可能会预期整个东西都会平滑地旋转，就像一个螺旋楼梯。

但事实并非如此。相反，这种材料表现得就像神经系统突然发生了“痉挛”。随着应力的积累，材料会突然断裂。当它断裂时，发生了一件怪事：靠近内圈的部分开始向一个方向旋转，而紧挨着它的一个区域却突然开始向相反的方向旋转。

这就像是你正在转动方向盘，突然间仪表盘的一个部分开始自行向后旋转。这被称为剪切局部化（shear localization）。这种“扭转”（或剪切）被卡在了一个非常特定、狭窄的区域内，而不是均匀地扩散开来。

旧理论 vs. 新理论

在此之前，科学家们试图使用标准的“弹性理论”（用于解释弹簧和橡胶的数学模型）来解释这一现象。

旧理论： 预测了一个平滑、温柔的曲线，即材料从中心向外逐渐扭转。
现实情况： 材料的表现却是混乱且锯齿状的。旧的数学模型失败得很彻底，就像试图用计算晴天情况的计算器来预测雷暴天气。

新的解决方案：“拓扑荷”

论文作者构建了一个新的数学模型来解释这种混乱。他们意识到，当材料断裂时，这不仅仅是一个简单的破碎。这就像材料是由微小的、看不见的磁铁（他们称之为“拓扑荷”或“四极矩”）组成的。

可以将这些电荷想象成粘稠液体内部微小的、看不见的漩涡。

设置： 当你旋转内圈时，你正在积聚压力，就像在绕紧一个弹簧。
断裂： 当压力过高时，这些看不见的漩涡会突然激活。它们会推挤并拉扯它们的邻居。
结果： 这些漩涡产生了一种“屏蔽”效应。它们在某些区域抵消了平滑的扭转，而在另一些区域则放大了扭转，从而创造了一个扭转方向发生翻转的清晰边界。

“魔数”

作者发现的一件最酷的事情是，材料的断裂并不是随机发生的。它似乎对特定的“魔数”（他们称为 $\kappa_e$ 的数学值）有某种偏好。

想象你正在调频收音机。你不能随便选一个频率；你必须捕捉到那些信号清晰的特定频道。

如果材料“调频”到其中一个特定数值，断裂就会发生在紧挨着内圈的地方。
如果它调频到另一个数值，断裂就会发生在远离圆环的粘稠物质中间。

论文表明，材料会根据圆环的大小自然地“选择”这些特定的设置，就像吉他弦会根据其长度在特定的音符上振动一样。

这为什么重要（根据论文所述）

作者创建了一个公式，可以精确预测这种“痉挛”会在哪里发生，以及粒子的运动轨迹看起来是什么样的，甚至可以精确到最小的细节。

他们不只是在猜测： 他们将新的数学公式与计算机模拟（虚拟实验）进行了对比，发现两者完美匹配。
他们解释了“为什么”： 他们展示了这些“漩涡”（塑性事件）才是导致材料呈现这种行为的原因，而不仅仅是简单的拉伸。
局限性： 他们指出，该理论适用于缓慢、谨慎的扭转。如果你转得太快或太用力（比如打破玻璃），规则就会改变，这是未来研究的一个课题。

总结

这篇论文就像是找到了混乱舞蹈的规则手册。材料（舞者）通常动作平滑，但当音乐（应力）变得过于响亮时，它们会突然陷入一种特定的、局部的旋转状态。作者弄清楚了告诉我们这种旋转何时发生以及为什么舞者会选择向相反方向运动的数学步骤。他们用一个新的“漩涡”模型取代了失效的“弹簧”模型，这个新模型完美地预测了这种混乱。

技术摘要：非晶态固体剪切局部化的解析理论

问题陈述
剪切局部化和剪切带是无定形固体在剪切应变下普遍存在的现象，常导致如金属玻璃等材料发生灾难性失效。虽然目前已确定这些现象源于塑性事件（雪崩）而非纯弹性，但现有的理论框架缺乏预测或描述在应力降过程中相关的详细位移场的能力。具体而言，对于在库埃特（Couette）几何构型约束下的准静态加载无定形固体，目前尚无能够对塑性雪崩期间的位移场（粒子运动）进行先验预测或描述的现有理论。此外，虽然剪切带（屈服时的系统跨越事件）是截然不同的现象，但控制屈服前后较小尺度剪切局部化的机制仍需要专门的理论方法。

研究方法
作者结合数值模拟与新推导的解析理论来解决这一空白。

数值模拟： 研究利用了一个由 $N=3,534$ 个双分散排斥圆盘（半径分别为 1 和 1.4）组成的二维无定形固体模型，这些圆盘被限制在两个同心环（ $R_{in}=28$ , $R_{out}=80.8$ ）之间。系统在面积分率 $\phi=0.91$ 下制备并达到机械平衡。应用了一种准静态加载协议，其中内环以微小步长（ $\delta\theta = 0.024^\circ$ ）进行逆时针旋转，并在每一步后通过阻尼分子动力学进行系统再平衡。外边界是固定的。研究重点关注与应力降（塑性事件）相关的位移场，而非弹性响应。
解析理论： 作者扩展了将塑性事件视为位移场中拓扑四极子和偶极子的现有理论框架。
- 他们构建了一个拉格朗日量 ( $U$ )，其中包括四极子的代价、四极子-四极子相互作用，以及应变场与四极子场之间的耦合。
- 至关重要的是，该理论考虑了在塑性事件发生前积累的大量背景应力 ( $\Sigma_{r\theta}$ )，这是以往单步加载研究未充分解决的因素。
- 通过对拉格朗日量关于四极子场 ( $Q$ ) 和位移场 ( $d$ ) 求极小值，作者推导出了一个重整化弹性张量。
- 引入偶极子贡献（源自四极子场的梯度）导出了位移场的控制微分方程： $\Delta d + (1 + \tilde{\lambda})\nabla(\nabla \cdot d) + \Gamma d = 0$ 。
- 使用傅里叶级数展开法求解该方程，得到径向 ( $d_r$ ) 和切向 ( $d_\theta$ ) 分量，并受限于库埃特几何的具体边界条件。

主要贡献与结果

位移场的解析预测： 主要贡献是推导出了角平均切向位移场 $d_\theta(r)$ 的解析解（公式 14）。该解通过贝塞尔函数 ( $J_1$ 和 $Y_1$ ) 以及筛选参数 $\kappa_e$ 来表示，成功重现了模拟中观察到的复杂位移剖面，包括经典弹性理论无法预测的粒子旋转方向反转（顺时针 vs 逆时针）。
剪切局部化机制： 理论表明，剪切局部化由离散的筛选参数 $\kappa_e$ $κ_{e}$ 值决定。这些值仅由系统的几何参数（ $R_{in}$ $R_{in}$ 和 $R_{out}$ $R_{o u t}$ ）决定，而与材料的具体状态无关。
- 当 $\kappa_e \approx 0.067$ 时，剪切局部化发生在靠近内边界处。
- 当 $\kappa_e \approx 0.123$ 时，剪切局部化发生在材料本体中。
- 系统优先选择使解析解分母趋于零（奇异点）的 $\kappa_e$ 值，这对应于最大的系统响应。
验证： 解析预测与数值模拟表现出极佳的一致性。该理论正确预测了在分析的 200 个最大应力降中，大多数应力降的局部化位置和形状。论文指出，大的应力降倾向于在边界附近局部化，而较小的应力降则在本体中局部化，这与不同 $\kappa_e$ 模式的选择一致。
与剪切带的区别： 作者澄清，其线性理论（在拉格朗日量中展开至二次项）适用于背景应力为无散形式的准静态剪切步骤。文中明确指出，该理论并不描述与屈服相关的剪切带，后者需要超越本文范围的非线性扩展（例如添加四次偶极项）。

意义
该论文声称提供了首个能够预测剪切无定形固体中塑性降相关的角平均位移场的解析理论。通过引入拓扑电荷（四极子和偶 Dipole）并考虑背景应力，该理论弥合了纯弹性（无法描述这些现象）与复杂数值模拟之间的鸿沟。根据几何约束先验预测剪切局部化发生位置的能力，为无定形材料的延性响应提供了基础性的理解。作者强调，虽然该理论解释了延性固体的局部化，但脆性剪切带的机制仍有待未来的非线性扩展研究。