全拉格朗日

描述一个物体在力作用下如何变形，是物理学和工程学中的一个核心挑战。当这些变形很大，涉及形状和方向的显著变化时，选择一个一致的参考系就变得至关重要。全拉格朗日（TL）列式提供了一个强大而优雅的解决方案，它采用了一个单一、固定的视角：物体的初始、未变形状态。所有后续的运动、应变和应力都相对于这个原始“蓝图”来描述，从而创建了一个关于变形过程的统一历史记录。这种方法与那些不断更新其参考系的方法形成对比，它为复杂的非线性问题揭示了深刻的见解并带来了计算上的好处。

本文将对全拉格朗日列式进行全面探索。在第一章 ​​“原理与机制”​​ 中，我们将剖析该框架的数学语言，从基本的变形梯度到构成其稳健性基础的格林-拉格朗日应变和皮奥拉-基尔霍夫应力张量等关键概念。随后，​​“应用与跨学科联系”​​ 一章将展示该理论的实践威力，揭示其与超弹性材料建模的天然契合，其预测结构屈曲的能力，以及其在地质力学和现代计算方法等不同领域中的相关性。

想象你是一位负责记录一座城市生命历程的历史学家。你可以站在街角，描述人流和车流的实时变化——这只是时间的快照。或者，你可以拿出城市的原始蓝图，即它的创始地图，并描述数个世纪以来，每一栋建筑和街道相对于那个初始规划是如何拉伸、收缩或移动的。第一种方法是短暂的，其视角总在更新。第二种方法则是锚定的，通过一个单一、不变的参考系的镜头来观察整个复杂变形历史。

在物理学和工程学领域，当我们研究物体如何弯曲、扭转和变形时，我们面临同样的选择。​​全拉格朗日（TL）​​ 列式是第二种哲学的宏伟体现。它是一个强大而优雅的框架，选择物体初始的、未变形的状态——其“创始地图”——作为所有时间的唯一参考。所有的运动、所有的力、所有的应变都从这个固定的、历史的视角来观察。这个选择并非随意的；正如我们将看到的，它为我们揭示了对变形力学深刻而统一的理解。它的对应方法，即​​更新拉格朗日（UL）​​ 列式，则更像是街角的观察者，不断将其参考系更新到最新已知的构型。通过坚守原始蓝图，全拉格朗日方法获得了独特的清晰度，并且在许多情况下，获得了显著的计算效率。

为了使我们的历史记录精确，我们需要一种数学语言。假设我们物体中的一个粒子在原始蓝图上（我们称之为​​参考构型​​，$\Omega_0$）的初始位置为 $\mathbf{X}$。在受到一些力作用后，这个粒子移动到我们当前所见世界中的一个新位置 $\mathbf{x}$（即​​当前构型​​，$\Omega_t$）。每个粒子的运动轨迹都由一个称为​​运动​​的映射或函数记录下来：$\mathbf{x} = \boldsymbol{\varphi}(\mathbf{X})$。

这个故事中最重要的角色是​​变形梯度​​，记为 $\mathbf{F}$。它被定义为运动相对于原始坐标的梯度：

$\mathbf{F} = \frac{\partial \mathbf{x}}{\partial \mathbf{X}}$

不要被这个符号吓到。你可以把 $\mathbf{F}$ 想象成一台小机器。如果你给它输入一个来自原始蓝图的微小箭头（一个向量），它会告诉你这个箭头在变形后的物体中是如何被拉伸和旋转成一个新箭头的。

让我们把它具体化。想象一个简单的例子，其中变形是均匀的，由一个线性变换 $\mathbf{x} = \mathbf{A}\mathbf{X}$ 描述，其中 $\mathbf{A}$ 是某个矩阵。在这种情况下，变形梯度 $\mathbf{F}$ 就是常数矩阵 $\mathbf{A}$ 本身。假设我们已知一小块材料的该矩阵：

接近1的对角项告诉我们沿轴线的拉伸情况，而非对角项则告诉我们剪切或角度变化的情况。这单个张量 $\mathbf{F}$ 包含了整个局部变形的信息。

这台机器还有一个秘密功能。如果我们计算它的行列式，$J = \det(\mathbf{F})$，我们会得到一个单一而强大的数字。对于上面的矩阵，计算结果为 $J \approx 1.04$。这是什么意思呢？这是局部的体积变化！$J=1.04$ 的值告诉我们，原始物体中一个微小的体积在变形后的物体中膨胀了约4%。如果 $J$ 小于1，表示压缩；如果 $J$ 恰好为1，则变形是保持体积的。雅可比行列式 $J$ 是当前微元体积与其原始体积之比，$J = \frac{dv}{dV}$。

现在我们面临一个微妙但关键的问题。如果我们拿一把刚性钢尺简单地旋转它，它的粒子已经移动了。变形梯度 $\mathbf{F}$ 将不是单位矩阵。然而，这把尺子根本没有被拉伸或产生应力。我们对“应变”的度量必须足够智能，能够区分真实的形状变化和纯粹的刚体转动。它必须是“客观的”，或者说是标架无关的。

在这里，数学提供了一个优美的解决方案：​​极分解​​。它告诉我们，任何变形 $\mathbf{F}$ 都可以唯一地分解为两部分：一个纯转动 $\mathbf{R}$ 和一个纯拉伸 $\mathbf{U}$，使得 $\mathbf{F} = \mathbf{R}\mathbf{U}$。拉伸张量 $\mathbf{U}$ 描述形状的变化，而转动张量 $\mathbf{R}$ 描述方向的变化。我们的目标是找到一个只依赖于 $\mathbf{U}$ 并且完全不受 $\mathbf{R}$ 影响的应变度量。

诀窍是计算一个称为​​右柯西-格林张量​​的量 $\mathbf{C}$：

$\mathbf{C} = \mathbf{F}^{\mathsf{T}}\mathbf{F} = (\mathbf{R}\mathbf{U})^{\mathsf{T}}(\mathbf{R}\mathbf{U}) = \mathbf{U}^{\mathsf{T}}\mathbf{R}^{\mathsf{T}}\mathbf{R}\mathbf{U} = \mathbf{U}^{\mathsf{T}}\mathbf{U} = \mathbf{U}^2$

看发生了什么！因为对于一个转动矩阵，$\mathbf{R}^{\mathsf{T}}\mathbf{R}=\mathbf{I}$（单位矩阵），所以转动部分 $\mathbf{R}$ 从表达式中奇迹般地消失了。张量 $\mathbf{C}$ 只依赖于拉伸 $\mathbf{U}$。它成功地将纯变形从刚体转动中分离了出来。

由此，我们定义了全拉格朗日列式的基石应变度量：​​格林-拉格朗日应变张量​​ $\mathbf{E}$：

$\mathbf{E} = \frac{1}{2}(\mathbf{C} - \mathbf{I}) = \frac{1}{2}(\mathbf{F}^{\mathsf{T}}\mathbf{F} - \mathbf{I})$

如果没有变形（$\mathbf{F}=\mathbf{I}$），那么 $\mathbf{E}=0$。如果只有刚体转动，我们刚看到 $\mathbf{C}=\mathbf{I}$，所以 $\mathbf{E}$ 也为零。只有当物体确实被拉伸或剪切时，格林-拉格朗日应变 $\mathbf{E}$ 才不为零。这个优雅的特性使得TL列式在处理涉及大转动但实际应变很小的问题时非常稳健，例如飞机机翼或风力涡轮机叶片的弯曲。它能正确地报告纯转动的应变为零，从而避免了虚假的、非物理的应力。

我们现在有了一个合适的应变度量。但应力呢？事实证明，在大变形的世界里，“应力”不是一个单一的概念，而是一个家族，每个成员都有不同的用途，就像不同国家的货币一样。

最直观的应力是我们在入门物理学中学到的：​​柯西应力​​ $\boldsymbol{\sigma}$。它是作用在当前变形状态下“真实”面积上的“真实”力。这是如果你被嵌入材料中会物理上感受到的应力。

支配平衡的最高法则是​​虚功原理​​。这是一个关于能量平衡的深刻陈述：对于任何无限小的、运动学上可能的（“虚”）运动，内应力所做的功必须等于外部施加的力所做的功。在当前构型中，这被写成一个包含柯西应力 $\boldsymbol{\sigma}$ 的积分。

但全拉格朗日列式的哲学要求我们将整个法则转换回参考构型 $\Omega_0$。这个“拉回”操作就像将一家全球公司的所有金融交易都转换为单一的本国货币。当我们进行这种数学转换时，一种新型的应力诞生了。内功项从 $\Omega_t$ 拉回到 $\Omega_0$ 时，自然产生了​​第一皮奥拉-基尔霍夫（PK1）应力​​ $\mathbf{P}$。在TL列式中，平衡方程的弱形式变为：

$\int_{\Omega_0} \mathbf{P} : \nabla_0 \delta \mathbf{u}\,\, \mathrm{d}V_0 = \text{External Virtual Work}$

PK1应力 $\mathbf{P}$ 是一个奇特的混合体。它表示作用在原始参考构型面积上的当前构型中的力。它不像柯西应力那样物理直观，但它恰恰是在参考系中陈述虚功法则所需要的“货币”。我们说 $\mathbf{P}$ 与变形梯度 $\mathbf{F}$（或更精确地说，其变化率 $\dot{\mathbf{F}}$）是​​能量共轭​​的。它们是在拉格朗日描述中表达机械功率的自然配对。

但故事并未就此结束。对于许多材料，如橡胶，其变形抗力来自于储存的弹性能，这一概念被称为​​超弹性​​。这个储存的能量 $\Psi$ 不能依赖于物体在空间中的方向（这是一个称为​​客观性​​的原理）。因此，$\Psi$ 不能是包含转动的完整变形梯度 $\mathbf{F}$ 的函数。它必须是一个纯粹的、客观的变形度量（如格林-拉格朗日应变 $\mathbf{E}$）的函数。所以，我们写作 $\Psi = \Psi(\mathbf{E})$。

现在，谜题的最后一块拼图就位了。通过对储存能量 $\Psi$ 关于应变 $\mathbf{E}$ 求导而产生的应力，必须是与 $\mathbf{E}$ 能量共轭的应力。这一发现过程将我们引向​​第二皮奥拉-基尔霍夫（PK2）应力​​ $\mathbf{S}$。它的定义使得 $\mathbf{S} = \frac{\partial \Psi}{\partial \mathbf{E}}$，并且事实表明，内功率密度可以等效地表示为 $\mathbf{S}:\dot{\mathbf{E}}$。

这揭示了理论内部一个深刻而美丽的二元性：

• ​​第一皮奥拉-基尔霍夫应力（$\mathbf{P}$）​​ 是动力学的应力。当虚功原理被转换到参考系时，它自然产生。
• ​​第二皮奥拉-基尔霍夫应力（$\mathbf{S}$）​​ 是材料本构的应力。它自然地来自于储存能量的物理学和客观性原理。

这两种应力“货币”并非相互独立；它们通过变形本身联系在一起，通过简单的变换关系 $\mathbf{P} = \mathbf{F}\mathbf{S}$。这种相互联系提供了一个统一的结构，将抽象的运动定律与材料的实际行为联系起来。弱形式可以用任一种应力度量来书写，因为 $\mathbf{P} : \nabla_0 \delta \mathbf{u}$ 等价于 $\mathbf{S} : \delta\mathbf{E}$。

我们为什么要接受这个看似复杂的多重应力和应变世界？回报体现在计算领域，当我们使用有限元法（FEM）来解决现实世界的工程问题时。在这里，全拉格朗日列式的优雅转化为了具体的优势。

首先，​​它解决了网格畸变的问题​​。在模拟中，所有积分都是在有限元网格的每个单元内的特定“高斯点”上进行数值计算的。在TL列式中，这些计算总是在初始的、未变形的网格上执行。在当前构型中，网格可能会严重弯曲、拉伸和扭曲，但我们的计算仍然在原始蓝图的那个原始、有序的网格上进行。这极大地提高了数值的准确性和稳定性，尤其是在涉及大转动的问题中。

其次，​​它在计算上是高效的​​。因为参考构型永不改变，模拟所需的许多量——例如单元形函数的梯度——可以在最开始计算一次并储存起来。相比之下，更新拉格朗日列式必须在每个求解迭代中都在变形的网格上重新计算这些量。TL中的这种预计算节省了大量工作，使得模拟的每一步都更快。

最后，​​它简化了材料本构律的实现​​。通过使用客观对偶（$\mathbf{S}, \mathbf{E}$），超弹性材料模型的实现变得清晰而直接。人们避免了在UL列式中为了正确处理转动而需要的复杂且有时棘手的“客观应力率”，这简化了算法并减少了潜在的错误源。

归根结底，全拉格朗日列式远不止是一种数学技巧。它是一种深刻的物理选择，也是一种哲学上一致的选择。通过将我们的整个视角锚定在不变的过去，我们获得了一个清晰、强大且计算高效的框架，来理解事物变形的丰富而复杂的故事。

在遍历了全拉格朗日列式的原理和机制之后，人们可能会倾向于将其视为一种相当形式化的，甚至有些抽象的数学工具。但这样做无异于见树不见林。这种方法的真正美妙之处，正如物理学中的许多事物一样，不在于其形式主义，而在于其描述我们周围世界的非凡能力。通过在不变的参考构型中站稳脚跟，我们获得了一个固定的有利位置，从中观察变形那壮丽而往往复杂的舞蹈如何展开。本章就是对这种力量的探索——一次巡礼，看看在哪些领域，全拉格朗日的视角不仅提供了答案，更带来了深刻的见解。

全拉格朗日列式的天然归宿是超弹性的世界——橡胶等材料的柔软、可拉伸的领域。为什么这种结合如此完美？想象一下，试图描述一根被拉伸的橡皮筋中储存的能量。最自然的方式似乎是将该能量与其相对于初始、未拉伸状态的拉伸程度联系起来。超弹性材料正是如此：其应变能 $W$ 是变形本身的函数。全拉格朗日列式接纳了这种直觉。通过将应变能定义为格林-拉格朗日应变张量 $\mathbf{E}$ 的函数，或等效地定义为右柯西-格林张量 $\mathbf{C}$ 的函数，我们自动满足了一个深刻的物理原理：材料标架无关性。储存的能量应该只取决于拉伸和剪切，而不取决于物体是否也在空间中进行刚性转动。由于 $\mathbf{E}$ 对这种刚性转动是“盲目”的，将 $W$ 定义为 $W(\mathbf{E})$ 是一种将此原理构建到我们模型根基中的优雅方式。

这个优雅的选择带来了美好的结果。因为应力和应变度量（$\mathbf{S}$ 和 $\mathbf{E}$）是从单一能量势函数导出的功共轭对，所以得到的系统本质上是能量守恒的。对于数值模拟来说，这不仅仅是一种美学上的胜利。它导出了对称的切向刚度矩阵，这是用于求解非线性方程的牛顿-拉夫逊法的核心引擎。对称矩阵意味着更快、更稳定的算法。这意味着对于涉及大弹性变形的问题，全拉格朗日方法通常提供了一条通往解决方案的路径，该方案不仅准确，而且具有二次方收敛性，这是数值分析中的黄金标准。

但全拉格朗日列式的影响范围远不止于简单的弹性体。考虑一下结构工程的广阔世界，及其桥梁、飞机和建筑物。在这里，我们经常遇到应变很小（钢材本身几乎不变形），但位移和转动可能巨大的情况。一根在荷载下屈曲的细长柱子就是一个典型例子。

结构分析中一个有趣的现象是“P-Delta”（$P$-$\Delta$）效应，即一个柱子抵抗侧向力的能力会因为存在压缩轴向荷载 $P$ 而降低。这似乎是一个复杂的二阶相互作用，但在全拉格朗日框架内，这种效应以惊人的自然方式出现。当我们一致地线性化控制方程时，切向刚度矩阵自动分裂为两部分：一个依赖于材料杨氏模量的材料刚度，和一个直接依赖于结构中预先存在的应力的几何刚度 $K_G$。正是这个几何刚度项精确地捕捉了 $P$-$\Delta$ 效应。该列式不仅仅是让我们加入这种效应；它从第一性原理出发，预测了其存在和形式。同样稳健的框架可以扩展到分析复杂的壳和板结构，其中结构表面的大转动通过跟踪一个导向场来处理，所有这些都从固定的参考构型进行一致描述。

几何刚度 $K_G$ 不仅仅是一个修正项；它是理解力学中最具戏剧性的事件之一——结构失稳——的关键。屈曲不是材料的失效，而是平衡路径的失效。它是一个分岔点，一个结构突然发现一个新的、截然不同的形状在能量上更有利的点。

我们如何预测这个临界点？全拉格朗日列式给了我们一个明确的答案。结构的总刚度是材料刚度（通常是正的，试图恢复原始形状）和几何刚度（在受压时可以是负的，并促使变形）的组合。随着压缩荷载的增加，几何刚度变得越来越负，从而“软化”了结构。屈曲发生在总切向刚度矩阵变为奇异的确切时刻——即它失去抵抗某种变形模态的能力的时刻。忽略几何刚度 $K_G$ 将会错过整个现象；它是问题的核心。

此外，全拉格朗日列式，特别是对于超弹性材料，非常适合追踪结构在*屈曲后*的行为。这种“后屈曲”分析对于理解结构的安全性和剩余强度至关重要。该列式的能量守恒特性确保了预测的路径是物理上现实的，没有可能困扰其他方法的人为能量耗散，从而即使在存在非常大的转动时也能实现稳健而准确的模拟 [@problem_d:2673060]。

在这一点上，您可能会想，这个美妙的框架是否仅限于清洁、可逆的弹性世界。对于像金属那样会屈服和永久变形的、具有复杂历史相关行为的材料呢？或者对于土壤和岩石的复杂响应呢？

值得注意的是，全拉格朗日运动学也为这些问题提供了稳健的支架。对于弹塑性，材料的内部状态（其屈服程度）成为我们必须跟踪的额外变量集。将应力与应变联系起来的本构律不再是一个简单的函数，而是一个增量法则。然而，我们仍然可以定义一个一致的材料切线模量 $C^{ep}$，它描述了材料应力如何响应一个小的附加应变。TL列式的整体结构保持不变：我们有一个材料刚度部分和一个几何刚度部分。关键在于，在与其他列式进行比较时，要正确地将材料切线模量“推前”到空间构型中，确保无论采用何种描述，物理学都保持一致。这种模块化——一个能够容纳不同本构模型的运动学框架——是一个强大物理理论的标志。

在地质力学中，如滑坡或土体固结等会发生巨大变形的场景，描述方式的选择成了一个深刻的哲学问题。全拉格朗日方法追踪单个土壤“包裹”的命运，追问“这个始于位置$\boldsymbol{X}$的物质最终会到哪里去？”这与流体力学中常见的纯欧拉描述形成鲜明对比，后者将其视点固定在空间位置$\boldsymbol{x}$上，并提问“现在有什么物质流过我这里？”两种观点都有其用武之地，但拉格朗日观点对于追踪固体材料变形时的演化属性、损伤或历史通常更为自然。

我们所讨论的原理并非过去的遗物，仅限于经典的有限元教科书。它们在计算力学的前沿领域依然充满活力。在现代无网格方法中，如在模拟流固耦合和颗粒流中广受欢迎的光滑粒子流体动力学（SPH），连续体由一团相互作用的粒子表示。即使在这个看似不同的世界里，全拉格朗日列式的概念也找到了其天然的归宿。梯度是相对于初始粒子位置计算的，整个变形历史通过参照那个原始的、未变形的状态来捕捉。

这种统一的力量也延伸到处理现实世界边界条件的复杂性。如果一个力，如流体压力，作用于一个本身正在变形的表面的法向，会发生什么？这是一个“跟随荷载”，它提出了一个挑战，因为荷载的方向取决于解本身。全拉格朗日列式通过一个巧妙的变换来处理这个问题。利用一个名为Nanson公式的数学工具，我们可以“拉回”这个空间定义的荷载，并将其表示为作用在固定的参考表面上的等效力。这使我们能够将荷载整合到我们的方程中。一个有趣的转折随之出现：这个跟随荷载的一致线性化会在我们原本对称的刚度矩阵中引入一个非对称项。这是一个美妙的例子，理论告诉我们一些深刻的东西：为了准确捕捉非保守力的物理特性，我们有时必须牺牲对称性的数学优雅。当然，一旦力和位移在参考构型上定义好，我们就可以使用各种标准的数值技术来施加它们，无论是通过划分系统直接施加，还是通过罚函数法或拉格朗日乘子进行近似施加。

在这次巡礼中，我们经常将全拉格朗日（TL）列式与其姊妹方法——更新拉格朗日（UL）列式进行对比，后者在每一步都重新定义参考状态。必须记住，这些不是相互竞争的物理理论；它们是针对相同底层物理的不同记账系统。对于任何给定的最终变形形状，如果两种方法都正确实现，它们必须预测出相同的最终物理状态——相同的应力、相同的应变。在此过程中，我们张量中的数字可能看起来不同（$\mathbf{E}$ vs. $\mathbf{e}$），但它们最终描述的是同一个现实。

当我们观察小变形的极限时，这种统一性最为明显。随着转动和应变变得无穷小，当前构型和参考构型之间的区别消失了。变形梯度 $\mathbf{F}$ 趋近于单位张量 $\mathbf{I}$，柯西应力 $\boldsymbol{\sigma}$ 与第二皮奥拉-基尔霍夫应力 $\mathbf{S}$ 变得无法区分，并且TL和UL列式都优雅地简化为我们都最先学习的那个熟悉的线性理论。非线性力学的复杂机器在其应有的极限下，与更简单的线性世界展现出一种美丽的对应关系。

最终，全拉格朗日列式不仅仅是一种巧妙的计算技巧。它是一种哲学选择。它选择相信，通过采用一个固定的、不变的视角——材料的诞生状态——人们可以解开最复杂的变形历史。它有力地证明了一个观点：在力学中，如同在许多事物中一样，了解你从哪里来，是知道你身在何处的关键。