Variational approach to determine the properties of dislocations at finite deformation

本文为存在位错的有限变形弹性理论建立了变分基础，证明了在有限变形框架中引入这些缺陷并非平凡过程，并导致作用在位错段上的力偏离了经典的 Peach-Koehler 力。

要点

想象一块金属，比如一根铜线或一根钢梁。在肉眼看来，它显得坚实且光滑。但如果你将其放大一百万倍，你会发现它实际上是一个晶格，一个由原子构成的完美有序网格。当你弯曲或拉伸这块金属时，它并不会像橡皮筋那样直接弹回原状；它会发生永久性的形状改变。这被称为塑性变形。

你提供的论文解释了在微观层面上这是如何发生的，并为描述金属在发生显著弯曲时的数学规则奠定了基础。

以下是使用简单类比对论文思想进行的拆解：

1. 问题所在：太多的舞者

在金属内部，引起形状变化的“舞者”被称为位错（dislocations）。你可以把它们想象成在原子网格中移动的微小、灵活的线或褶皱。

• 挑战： 在一段弯曲的金属中，存在着数以万亿计的这些位错。试图追踪每一个位错（就像在拥挤的人群中追踪每一位舞者）对计算机来说太难了。
• 目标： 科学家们想要一种“连续介质理论”。与其追踪单个舞者，他们希望将整个“人群”描述为一个整体流体。这篇论文是关于为这种流体建立规则手册的，但特别针对的是金属被大幅度弯曲的情况（有限变形），而不仅仅是轻微弯曲。

2. 旧规则手册 vs. 新规则手册

长期以来，科学家一直使用“线性弹性理论”来描述这些材料。

• 旧方法（小变形）： 想象你只是轻微拉伸一根橡皮筋。数学非常简单：如果你施加两倍的力量，它就会拉长两倍。作用在位错（即“舞者”）上的力是众所周知的，且易于计算。这就像是 Peach-Koehler 力，一个大家都在使用的标准公式。
• 新方法（大变形）： 现在，想象你拉伸那根橡皮筋，直到它几乎达到断裂点。规则改变了。材料变得更硬，几何结构发生了扭曲，简单的数学不再适用。
• 论文的发现： 作者 István Groma 表明，当你显著拉伸金属时，推动位错的“力”并不是用于小幅度拉伸时的那个简单公式。它需要一个更复杂的新版本力。

3. “切割与滑动”类比

如何在完美的晶体中创造一个位错？

• 隐喻： 想象一副扑克牌。如果你将牌堆从中切开，然后将上半部分向右移动一张牌，你就创造了一个“阶梯”或“褶皱”在中间。这个褶皱就是位错。
• 数学问题： 在论文中，作者必须用数学来描述这种“切割”。他引入了一个概念，叫做塑性畸变（plastic distortion）。
• 转折点： 当金属发生大幅度弯曲时，计算这种“切割”的逆过程（弄清楚如何回到原始形状）非常棘手，因为数学涉及“尖峰”（狄拉克 $\delta$ 函数），这些尖峰代表了切割的锐利边缘。作者展示了如何通过数学手段将这些尖峰“平滑化”，从而使方程不至于崩溃。

4. “能量景观”法

为了确定金属如何稳定在新的形状，作者使用了一种变分法（Variational Approach）。

• 类比： 想象一个球在丘陵地带滚动。球总是倾向于滚向最低点（山谷），因为那是能量最低的状态。
• 应用： 金属就像那个球。它想要寻找能量最低的形状。作者使用一种数学工具（泛函导数）来询问：“如果我稍微扰动一下原子，能量是上升还是下降？”
• 结果： 通过寻找能量不再发生变化的地方（谷底），他推导出了平衡方程。这些规则告诉我们，在弯曲的金属内部，应力是如何分布的。

5. 核心结论：力发生了变化

这篇论文最重要的发现是关于 Peach-Koehler 力 的。

• 在旧世界里： 推动位错的力就像吹向风帆的一阵简单微风。
• 在新世界里（大变形）： 作者证明了当金属发生剧烈变形时，“风”改变了。这种力取决于一种新型的“有效应力”，它考虑到了材料本身已经被拉伸和旋转的事实。
• 为什么重要： 如果你在处理重度弯曲的金属时仍使用旧的简单公式，你的计算将会出错。你需要这种经过修正的新型力，才能准确预测金属的行为。

总结

这是一篇基础性的数学更新论文。它在说：“对于金属轻微弯曲的理论，我们已经做得很好，但当它们大幅度弯曲时，内部力的旧规则就不再适用了。我们使用了一种新的数学方法，推导出了适用于这些大幅度弯曲的正确规则。”

作者指出，这项工作是一个必要的垫脚石。一旦这些规则确立，它们就可以被用来构建一个更好的、更精确的计算机模型，用以预测在剧烈变形的材料中，复杂的位错网络是如何移动和相互作用的。

技术摘要

技术摘要：确定有限变形下位错性质的变分方法

问题陈述
曲率位错连续统理论的最新进展已成功描述了统计存储位错密度、几何必要位错密度以及曲率的时空演化。然而，现有的这一框架（Groma 等，2021）严格局限于小变形机制。将这些概念应用于有限变形理论，特别是在系统中存在单个位错的情况下，存在显著的空白。挑战在于如何在非线性弹性框架内推导平衡方程，并确定作用在位错段上的力，因为在这种框架下，标准的线性弹性假设（例如 Peach-Koehler 力的直接适用性）可能不再成立。

方法论
本文采用变分形式，利用能量对变形场的泛函导数来推导平衡条件。该方法通过以下步骤进行：

1. 无缺陷非线性弹性：作者首先通过将弹性能量 $E$ 定义为变形场 $\vec{x}(\vec{X})$ 的泛函，重新审视无缺陷情况。通过应用虚功原理和泛函微分，他们推导出了涉及第一和第二 Piola-Kirchhoff 应力张量的平衡方程。这一部分建立了处理必要非局部弹性项的形式体系。
2. 塑性变形分解：将变形梯度 $F_{ij}$ 分解为弹性分量 ($F^e_{ij}$) 和塑性分量 ($F^p_{ij}$) ($F_{ij} = F^e_{im}F^p_{mj}$)。弹性能量仅由弹性变形张量 $\epsilon^e_{ij}$ 定义。作者推导了含有塑性畸变的系统的平衡方程，表明平衡方程中的应力项必须被一个有效转换后的应力张量所取代。
3. 引入单个位错：为了模拟单个位错环，通过一个带有位错（伯格斯矢量 $\vec{b}$）的切断面来定义塑性变形梯度。这引入了一个包含 Dirac $\delta$ 函数的塑性畸变 $\beta^p_{ij}$。
   • 正则化：意识到由于 $\delta$ 函数的存在，塑性变形梯度的逆 ($F^{-p}_{ij}$) 会在数学上变得奇异，作者提出了一个正则化程序。他们用一个宽度在伯格斯矢量数量级的正则函数（例如高斯函数）来近似 $\delta$ 函数。这使得能够一致地定义至多包含 $\vec{b}$ 的线性项的 $F^{-p}_{ij}$，从而避免在能量计算中出现奇异性。
4. 力推导：通过计算位错线虚位移导致的总能量变化来推导作用在位错段上的力。这涉及在保持位移场固定的情况下，计算能量对塑性畸量 ($F^{-p}_{ij}$) 的泛函导数，并随后求解平衡方程。

主要贡献与结果

• 修正的平衡方程：本文推导了包含塑性变形的物体的平衡方程。研究表明，在有限变形机制下，出现在平衡方程中的应力是一个有效转换后的应力 ($\sigma^{2PK*}$)，它不同于通过能量对弹性应变求导得到的标准第二 Piola-Kirchhoff 应力。
• 广义 Peach-Koehler 力：主要结果是推导了有限变形下作用在位错段上的力。作者表明，该力并非简单地源自线性弹性的经典 Peach-Koehler 力。相反，它取决于一个有效应力张量 $\sigma^f_{ij}$，定义为：
  $$\sigma^f_{ij} = F^T_{ip} F^e_{pl} \sigma^{2PK}_{lj}$$
  由此产生的单位长度力为 $\vec{f}_{PK} = (\hat{\sigma}^f \vec{b}) \times \vec{t}$，其中 $\vec{t}$ 是线切向量。在小变形极限下，该表达式恢复了标准形式。
• 耗散性质的运动：作者证明，如果位错速度与该力在滑移面上的投影成正比，则弹性能量的变化率是非正的 ($\dot{E} \leq 0$)。这证实了即使在有限变形框架内，位错运动仍保持耗散特性。
• 处理奇异性：本文通过将 $F^{-p}_{ij}$ 定义为主要物理量并利用正则化，提供了一种处理与单个位错环的塑性畸变相关的数学奇异性的系统方法。

意义与主张
本文声称，将位错引入有限变形框架是一项非平凡的任务，需要系统的变分方法。作者断言，其结果为将曲率位错的连续统理论（此前局限于小变形）推广到有限变形机制提供了必要的理论基础（“关键输入”）。

至关重要的是，本文强调，对于基于位错网络集体演化的通用塑性理论，大变形框架是必不可少的。推导出的结果表明，在有限变形机制下，用于平衡方程的“应力度量”与作用在位错上的力并不相同，这一区别必须在未来的连续统模型中予以考虑。作者指出，这些发现将被应用于后续论文中，以实现曲率位错在有限变形下的位错连续统理论的推广。