皮奥拉-基尔霍夫应力

在分析物体如何响应力时，应力——单位面积上的力——这一概念是基础。然而，对于经历大变形的材料，如橡胶或生物组织，一个重大的挑战出现了：力作用的形状和面积在不断变化。直观上的“真实”应力，即柯西应力，是定义在这种变化的、已变形的形状上的，这使得计算变得复杂，并模糊了材料的内在属性。因此，需要一个更稳定的分析框架。

本文通过引入皮奥拉-基尔霍夫应力张量来弥合这一差距，该张量通过将变形状态下的力与物体的初始、未变形形状联系起来，提供了一种强大的替代方法。在接下来的章节中，您将对这一基本概念有全面的理解。“原理与机制”一节将揭示其理论基础，区分柯西应力、第一皮奥拉-基尔霍夫（PK1）应力和第二皮奥拉-基尔霍夫（PK2）应力，并解释它们的数学和能量关系。之后，“应用与跨学科联系”一节将展示这些理论工具如何在计算力学、材料科学到生物力学等领域中不可或缺，从而能够准确模拟和理解我们复杂的物理世界。

要真正理解可变形物体的世界——从被拉伸的橡皮筋到生物组织的复杂折叠——我们必须适应一种“双重人格”。我们需要同时思考一个物体“现在”的样子和它“过去”的样子。这种思维上的二元性是连续介质力学的核心，迫使我们在两个截然不同的世界之间穿梭：​​当前构形​​和​​参考构形​​。

当前（或空间）构形是物体此时此地、处于变形状态下的形状。这是我们能看到和触摸到的世界。参考（或物质）构形是物体在某个初始、未变形状态下的快照。这是“之前”的世界。我们用大写字母，如 $\mathbf{X}$，表示参考世界中的点，用小写字母，如 $\mathbf{x}$，表示当前世界中的点。物体的整个运动过程由一个函数 $\mathbf{x} = \boldsymbol{\varphi}(\mathbf{X}, t)$ 捕捉，该函数告诉我们每个起始于 $\mathbf{X}$ 的质点在时间 $t$ 移动到了哪里。材料的局部拉伸和旋转由一个极其重要的量来描述，即​​变形梯度​​ $\mathbf{F} = \partial\mathbf{x}/\partial\mathbf{X}$。

这种双重世界的观点不仅仅是一种哲学选择，更是一种实践上的必需。想象一下，要对一场车祸进行计算机模拟。金属部件发生剧烈变形。如果对应力的描述与当前扭曲的形状绑定，那么计算网格就必须不断地扭曲和更新——这是一个艰巨的挑战。如果能在汽车原始、完好的形状上进行所有计算，岂不是更方便？为此，我们需要能够在两个世界之间架起桥梁的工具。这就是不同类型应力发挥作用的地方。

当我们初次学习应力时，我们认为它是单位面积上的力。这个简单、直观的图像就是​​柯西应力​​张量的本质，用 $\boldsymbol{\sigma}$ 表示。它是“真实”应力，因为它是在此时此地测量的：力是当前的力，面积是变形体实际的、当前的面积。

Augustin-Louis Cauchy 的天才之处在于，他意识到在材料内部，任何假想切割面上的面力 $\mathbf{t}$ 都与该面的方向（由其法向量 $\mathbf{n}$ 给出）存在简单的线性关系。这种关系定义了应力张量：

这个方程中的每一项——单位面积上的力 $\mathbf{t}$、应力张量 $\boldsymbol{\sigma}$ 和面法线 $\mathbf{n}$——都是在当前变形构形中测量的。经典材料中角动量守恒的一个深刻推论是，柯西应力张量总是对称的（$\boldsymbol{\sigma} = \boldsymbol{\sigma}^T$），这意味着它本身没有内在的“扭转”。这是我们最熟悉的应力，是你在加载的桥梁上用应变计测量到的应力。

虽然柯西应力在物理上很直观，但其定义在不断变化的当前构形上，使得它在许多计算中显得 cumbersome。这促使我们发明一种新的应力度量，以便于我们从“参考世界”的视角进行分析。让我们定义一种应力，它将我们“现在”看到的力与“最初”的面积联系起来。

这就是​​第一皮奥拉-基尔霍夫（PK1）应力​​张量 $\mathbf{P}$ 的任务。它回答了这样一个问题：“作用在一个在未变形状态下单位面积的面上的当前力是多少？” 我们定义一个​​名义面力​​ $\mathbf{T}_R$，即单位“参考”面积上的当前力。PK1 应力将这个面力与“参考”法线 $\mathbf{N}$ 联系起来：

注意 $\mathbf{P}$ 的奇特性质。它接收一个来自参考构形的向量（$\mathbf{N}$），并将其映射到存在于当前构形中的一个力向量（$\mathbf{T}_R$）。因为它连接了两个不同的世界，所以它通常被称为​​两点张量​​。

为了将这个新应力与我们熟悉的柯西应力联系起来，我们使用一个简单而强大的物理原理：作用在一小块材料上的实际力是一个物理现实，与我们的数学描述无关。作用在一小块当前面积 $da$ 上的力必须等于作用在其对应的参考面积 $dA$ 上的力。这给了我们 $\mathbf{t}\,da = \mathbf{T}_R\,dA$。将此与应力的定义结合起来，我们得到 $\boldsymbol{\sigma}(\mathbf{n}\,da) = \mathbf{P}(\mathbf{N}\,dA)$。这个谜题的最后一块是知道有向面积元如何变换。这个几何关系就是​​南森公式​​，它表明 $\mathbf{n}\,da = J\mathbf{F}^{-T}\mathbf{N}\,dA$，其中 $J = \det(\mathbf{F})$ 是局部体积变化。代入这个关系，我们得到基本的转换公式：

这个优美的公式是我们的桥梁，允许我们将应力从当前世界转换到更方便的混合描述中。对于不可压缩材料，体积保持不变，$J=1$，关系简化为 $\mathbf{P} = \boldsymbol{\sigma}\mathbf{F}^{-T}$。

从这种关系中浮现出 PK1 应力一个有趣的特性。即使柯西应力 $\boldsymbol{\sigma}$ 是完全对称的，PK1 应力 $\mathbf{P}$ 通常是​​不对称​​的。这可能看起来很奇怪，但它并不违反任何物理定律。这是变形混合了各个分量的数学结果。例如，一个简单的剪切变形可以从一个简单的单轴柯西应力中产生一个不对称的 $\mathbf{P}$。

PK1 应力是一个重大的进步，因为它允许我们使用来自固定参考构形的法线 $\mathbf{N}$。但它仍然是一个两点张量，在当前世界中产生力。我们能否创建一个“完全”定义在参考世界中的应力度量？

答案是肯定的，其结果就是​​第二皮奥拉-基尔霍夫（PK2）应力​​张量 $\mathbf{S}$。我们通过变形梯度将 PK1 应力在数学上“拉回”来定义它：

这个定义可能看起来像一个纯粹的形式技巧，但它有一个奇妙的结果。如果我们代入 $\mathbf{P}$ 的表达式，我们发现 $\mathbf{S} = J\mathbf{F}^{-1}\boldsymbol{\sigma}\mathbf{F}^{-T}$。快速检查表明，如果 $\boldsymbol{\sigma}$ 是对称的，那么 $\mathbf{S}$ 也是对称的！。我们成功地构建了一个应力张量，它既完全存在于方便、不变的参考构形中，“又”与真实的柯西应力共享舒适的对称性。因此，$\mathbf{S}$ 是理论家们的最爱，并且对于构建描述材料内在行为的本构律至关重要。

不同的应力度量不仅仅是任意的数学构造。它们与能量和功的物理学有着深刻而优美的联系。在力学中，我们说一个力和一个速度是​​能量共轭​​的，因为它们的乘积得到功率。在连续介质力学中，特定的应力和应变率度量对是共轭的。无论我们使用哪种描述，单位体积材料所做的功率率必须相同，这决定了这些配对。

这个功率一致性原则揭示了一个完美对应的度量交响乐：

• ​​柯西应力​​ $\boldsymbol{\sigma}$ 与​​变形率​​ $\mathbf{d}$ 共轭。单位“当前”体积的功率是 $\boldsymbol{\sigma}:\mathbf{d}$。
• ​​第一皮奥拉-基尔霍夫应力​​ $\mathbf{P}$ 与​​变形梯度​​的变化率 $\dot{\mathbf{F}}$ 共轭。单位“参考”体积的功率是 $\mathbf{P}:\dot{\mathbf{F}}$。
• ​​第二皮奥拉-基尔霍夫应力​​ $\mathbf{S}$ 与​​格林-拉格朗日应变​​张量的变化率 $\dot{\mathbf{E}}$ 共轭。单位“参考”体积的功率是 $\mathbf{S}:\dot{\mathbf{E}}$。

这些配对不是可选的。任何用于新材料（无论是聚合物、金属还是活组织）的热力学一致模型都“必须”尊重这些能量伙伴关系。一个动手计算可以证实，对于任何给定的变形，以 $\mathbf{P}:\dot{\mathbf{F}}$ 计算的功率与以 $\mathbf{S}:\dot{\mathbf{E}}$ 计算的功率得到的数值完全相同，这表明了该框架的完美一致性。

在了解了这些不同的世界和应力度量之后，人们可能会想：为什么这没有在入门物理学中教授？答案在于当变形非常小时会发生什么。在一座下垂几毫米的钢桥的情况下，变形是微小的。变形梯度 $\mathbf{F}$ 几乎是单位矩阵（$\mathbf{F} \approx \mathbf{I}$），体积变化 $J$ 基本上为 1。

如果我们将这些近似值代入我们的转换公式，一个小小的奇迹发生了：

在小变形的极限下，所有三种应力度量——柯西、PK1 和 PK2——都收敛到相同的值！不同构形的复杂性消失了，我们只剩下一种单一、明确的“应力”。这就是为什么这种区别对于软体机器人、生物力学和橡胶弹性等领域如此关键（在这些领域，大变形是常态），但对于许多传统的土木和机械工程问题来说，却是不必要的复杂化。

最后，这些应力如何融入牛顿第二定律 $\mathbf{F}=m\mathbf{a}$？在连续体中，这个定律用应力散度来表示，它代表了由应力变化产生的净力。在当前构形中，运动方程是 $\rho_0 \mathbf{a} = \operatorname{Div}_X(\mathbf{P}) + \rho_0 \mathbf{b}_0$，其中所有量都在固定的参考坐标系中定义。为了将其与更熟悉的柯西应力联系起来，我们需要一个能在两个世界之间转换散度算子的“罗塞塔石碑”。这就是​​皮奥拉恒等式​​：

这个恒等式使我们能够以最方便的构形来书写基本的运动定律，并确信我们的物理预测将是相同的。它是统一两个世界的最后环节，使我们能够以优雅、一致和强大的方式描述复杂的变形之舞。

在了解了皮奥拉-基尔霍夫应力的原理和机制之后，我们可能会问自己：“所有这些数学工具是用来做什么的？” 这是一个合理的问题。既然我们已经有了直观的柯西应力，即存在于变形物体此时此地的“真实”应力，为什么还要发明新的方式来谈论应力呢？答案，正如物理学中常有的情况一样，是通过退后一步并改变我们的视角，我们可以发现对世界更深刻、更简单、更强大的理解。皮奥拉-基尔霍夫应力不仅仅是一种替代方案；它们是一把钥匙，解锁了描述和预测在各种科学和工程学科中经历复杂大变形的材料行为的能力。

想象一下，你正站在一张正在拉伸和扭曲的橡胶薄片上，试图描述它的地理状况。你的地标会不断移动，你的距离会改变，你画的地图在你画完的那一刻就已经过时了。这就是柯西应力的世界。现在，想象一下从一个固定的有利位置描述同一张拉伸的薄片，使用它原始的、平坦的、未变形的形状作为一张永久的地图。这就是皮奥拉-基尔霍夫应力的世界。它提供了一个固定的，或称“参考”构形，来描述一个变化的现实。

这种视角的转变是深刻的。一种材料的内在属性——它的刚度、它的弹性、它的本质——并不会仅仅因为被变形而改变。这些属性属于材料的参考状态。皮奥拉-基尔霍夫应力允许我们在这个不变的参考框架中阐述材料行为的规律，将材料的内在物理特性与其变形的纯几何后果分离开来。

这一思想在材料科学领域，特别是在为橡胶、聚合物和生物软组织等材料开发“本构模型”时，得到了最优雅的体现。这些被称为超弹性材料。它们的行为由一种储存的应变能来控制，这是一个势函数 $\Psi$，取决于材料偏离其自然状态的变形程度。

测量这种变形的自然方法是使用右柯西-格林张量 $\mathbf{C} = \mathbf{F}^T \mathbf{F}$，它完全是相对于参考构形定义的。这种方法的美妙之处在于应力自然地从能量中导出。具体来说，​​第二皮奥拉-基尔霍夫应力张量 $\mathbf{S}$​​ 作为能量对应变的导数而出现：$\mathbf{S} = 2 \frac{\partial \Psi}{\partial \mathbf{C}}$。这个应力张量是对称的，并且完全存在于参考框架中。它是“材料应力”，即材料自身所感受到的应力。例如，对于用于描述橡胶的复杂模型，如 Mooney-Rivlin 材料，其应力张量 $\mathbf{S}$ 可以直接从其应变能函数中导出，从而提供了对材料响应的完整且热力学一致的描述。

这种“以材料为中心”的观点可以揭示出惊人的物理现象。考虑一个承受简单剪切变形的材料块。如果我们在参考框架中使用对称张量 $\mathbf{S}$ 来描述这个简单的剪切状态，然后将这个状态“前推”到变形世界中看它是什么样子，结果得到的柯西应力 $\boldsymbol{\sigma}$ 惊人地复杂。它不仅包含预期的剪应力，还包含正应力——材料在一个没有被直接剪切的方向上向外推！这种现象，被称为 Poynting 效应，是非线性弹性的一个标志。它是大变形几何的直接后果，但其根本原因从 $\mathbf{S}$ 的参考视角来看最容易理解。简单的材料响应产生了复杂的空间应力。

同样，当一个各向同性材料沿其主轴拉伸时，“真实”柯西应力的主值 $\sigma_i$ 与“材料”第二皮奥拉-基尔霍夫应力的主值 $S_i$ 之间的关系异常简单：$\sigma_i = \frac{\lambda_i^2}{J} S_i$，其中 $\lambda_i$ 是主伸长率，$J$ 是体积变化。这个方程优雅地告诉我们，我们在变形世界中观察到的应力是内在的材料响应 $S_i$，被拉伸的几何形状 $\lambda_i^2$ 所放大。

拉格朗日（参考框架）视角的威力在计算科学中真正展现出来。当工程师模拟车祸、跳动的心脏或地震中土壤的行为时，他们通常使用有限元法（FEM）。对于大变形问题，最稳健的方法通常是​​全拉格朗日（TL）列式​​。在这种方法中，计算机在物体的初始、未变形状态下建立一个模型，并相对于这个固定的网格进行“所有”计算。

在这个计算世界中，皮奥拉-基尔霍夫应力不仅有用，而且是不可或缺的。虚拟物体内的内力是使用第二皮奥拉-基尔霍夫应力 $\mathbf{S}$ 和格林-拉格朗日应变 $\mathbf{E}$ 这对功共轭对来计算的。使用这对张量是因为它们都是“客观的”——如果物体仅仅在空间中刚性旋转，它们的值不会改变。这确保了模拟的能量平衡在物理上是正确的。

那么​​第一皮奥拉-基尔霍夫应力 $\mathbf{P}$​​ 在哪里发挥作用呢？这个张量由于其通常缺乏对称性而显得格格不入，但它扮演着两个框架之间信使或“运输者”的关键角色。它将参考框架中的一个面与作用在当前框架中该面上的力联系起来。这使得它成为在 TL 模拟中施加边界条件（如压力或面力）的完美工具。工程师可以在原始形状上指定一个力，而 $\mathbf{P}$ 确保该力在整个变形过程中被正确计算。

此外，为了求解复杂的非线性运动方程，计算机需要知道内力如何随着变形的变化而变化。材料的这种“切线刚度”由材料弹性张量 $\mathbb{C}_{iJkL}$ 描述，它被定义为第一皮奥拉-基尔霍夫应力对变形梯度的导数，$\mathbb{C} = \partial \mathbf{P} / \partial \mathbf{F}$。这使得 $\mathbf{P}$ 不仅在描述应力状态中至关重要，而且在求解算法本身中也处于中心地位。

这些概念的应用横跨了广泛的尺度和学科。在​​地质力学​​中，模拟冰川的流动或构造板块的折叠涉及长时间尺度上的巨大变形，这使得拉格朗日视角至关重要。

也许最令人兴奋的前沿是​​生物力学​​。人体是一个超弹性材料的宇宙。肌肉、皮肤、血管和韧带都在其正常功能中经历显著而复杂的变形。理解它们的力学对于设计医疗设备、规划手术和诊断疾病至关重要。皮奥拉-基尔霍夫应力和应变能函数的语言是描述这些软组织的原生语言。

但这个框架也告诉我们其自身的局限性以及与更简单理论的联系。考虑分析咀嚼负荷下人类下颌骨皮质骨等硬组织。在这里，应变通常非常小，大约在 $0.001$ 的量级。在这种无穷小变形的极限下，参考构形和当前构形之间的区别消失了。变形梯度 $\mathbf{F}$ 变得几乎是单位矩阵，体积变化 $J$ 基本上为一。结果，所有三种应力度量都趋于一致：$\boldsymbol{\sigma} \approx \mathbf{P} \approx \mathbf{S}$。我们回到了熟悉的线性弹性世界，在那里我们不需要区分不同的应力定义。这是一个美丽的教训：有限应变力学的复杂框架并没有抛弃更简单的理论，而是将它们作为一个特例包含在内。它向我们精确地展示了更简单的模型在何时以及为何是足够的，并在它们不足时提供了完整的理论。

归根结底，皮奥拉-基尔霍夫应力张量远不止是一个数学上的奇物。它们是一个统一的透镜。通过让我们采纳材料自身的视角，它们简化了其行为的基本规律，使强大的计算工具能够模拟我们的世界，并连接了我们对跨学科和跨尺度力学的理解。它们在看似混乱的变形之舞中揭示了潜在的简单性和秩序，这是物理学深刻而常常令人惊讶的美的证明。