玻尔百科在连续介质力学中，描述经历大形状变化的材料（即有限变形过程）内部的力是一个重大挑战。虽然直观的“真实”应力（柯西应力）是物理上可测量的，但它依赖于变形体不断变化的几何形状，这使得用它来建立预测性材料模型变得困难。这在物理观测和解析公式之间造成了知识鸿沟。本文通过探讨为有限变形分析设计的应力张量族来应对这一挑战。读者将首先浏览“原理与机制”部分，其中定义了柯西、皮奥拉-基尔霍夫和基尔霍夫应力张量，并通过能量和功的原理理解它们之间的基本关系。随后，“应用与跨学科联系”部分将展示这些理论工具，特别是第二类皮奥拉-基尔霍夫应力，在现代材料建模和计算模拟中如何发挥关键作用。
想象一下，你正在揉一个面包团，并试图描述其内部的力。面团在拉伸、折叠和流动。你刚才按下的地方现在已经移动到别处，原本是立方体的材料现在变成了一张扁平的薄片。这就是有限变形的复杂世界，它给物理学家和工程师们提出了一个根本性问题：当面积本身的大小和方向都在不断变化时，你该如何测量应力——这个单位面积上的力的简单概念？
如果你试图测量力，并将其除以你当前看到的面积，你所测量的就是柯西应力 (Cauchy stress)，用希腊字母 表示。这是此时此刻材料内部“真实”的、物理上可感知的应力。它就像一个嵌入面团中的微型假想压力计所读取的数值。然而，仅使用柯西应力来描述材料的行为，就像写一个每句话中人物和场景都在变化的故事一样。这在描述上是一场噩梦。为了建立一门预测性科学，我们需要一个固定的参考系。我们需要将复杂的变形状态与一个更简单的原始状态——未变形的面包团——联系起来。
为了弥合原始、未变形的世界（参考构形）与杂乱、被压扁的世界（当前构形）之间的差距，我们需要一本数学词典。这本词典就是变形梯度 (deformation gradient)，一个用 表示的张量。对于原始面包团中的每一个微小向量， 会告诉你这个向量在揉捏后的面团中变成了什么。它优雅地捕捉了材料所经历的所有拉伸、剪切和旋转。
有了这本词典，我们就可以发明一些新的应力类型，它们虽然在物理上不如柯西应力直接，但在分析上却方便得多。这些就是皮奥拉-基尔霍夫应力。
首先是第一类皮奥拉-基尔霍夫应力 (first Piola-Kirchhoff stress) ()。可以把它看作一个混合体。它测量作用在当前变形形状上的真实力，但将其除以变形发生之前的表面面积。这很有用，但它导致了一个奇怪的产物。 是所谓的“两点张量”(two-point tensor)；它的一个方向指向当前世界，而另一个方向指向参考世界。这就像试图用巴黎的街道名和东京的门牌号来指路。由于这种混合特性， 通常是非对称的。
这种不对称性的根本原因在于一个深刻的物理原理。柯西应力 必须是对称的。这不仅仅是数学上的便利，而是角动量平衡的直接结果。如果 不是对称的，微小的材料立方体会在没有任何外部扭转力的情况下自发旋转，这违反了物理学最基本的定律之一。对称的 和第一类皮奥拉-基尔霍夫应力 之间的数学关系是 ，其中 是体积变化量 ()。由于 通常不是对称的，从对称的 进行这种变换会得到一个非对称的 。
的这种非对称性令人不安。物理学家更喜欢对称张量的优雅。这就引出了第二个兄弟：第二类皮奥拉-基尔霍夫应力 (second Piola-Kirchhoff stress) ()。这种应力度量是完全“材料”的。它是通过将力向量本身从当前构形数学上拉回到参考构形来构建的。现在，力和面积都是从原始、未变形物体的角度来描述的。其变换定义为 ，当与之前的关系式结合时，我们得到回到柯西应力的联系：。
这个定义的美妙之处在于，如果你从对称的柯西应力 出发，得到的第二类皮奥拉-基尔霍夫应力 也是完全对称的。我们找到了一种从固定的参考视角描述应力的方法，同时保留了物理学所要求的优雅对称性。
所以，我们在运动的世界里有了“物理”应力 ，在固定的世界里有了“数学”应力 。那么，基尔霍夫应力 (Kirchhoff stress) 处于什么位置呢？它原来是连接这两种视角的一个绝妙而实用的桥梁。
基尔霍夫应力的定义非常简单，即柯西应力乘以体积变化量：
乍一看，这似乎有些随意。为什么是这个特定的缩放？答案，如同物理学中常见的情况一样，在于能量和功。在力学中，做功的速率——功率——是一个基本量。应力和应变通过功率联系在一起，这种关系被称为功共轭 (work conjugacy)。正如力与位移“共轭”产生功一样，应力张量与应变率张量共轭产生功率密度（单位体积的功率）。
连续介质力学框架的巧妙之处在于它提供了一组一致的配对 [@problem_id:2893483, @problem_id:3564564]：
基尔霍夫应力 在这里找到了它的用武之地。它与柯西应力一样，同变形率张量 功共轭，但它们的乘积给出的是单位参考体积的功率。
这太完美了！因子 正是当前体积与参考体积之比，因此将当前构形中的功率密度乘以 就得到了参考构形中的功率密度。基尔霍夫应力是一种“空间定位”的应力（如 ），它经过预先缩放，以便在能量上与“材料”参考系兼容。这使其成为计算模拟中，尤其是在有限元法 中一个极其宝贵的工具。它允许使用空间量进行计算，同时可以轻松地在固定的参考域中追踪能量。
这个优雅的关系网络可以通过一个强大的等价关系来概括：无论你如何测量，材料块内部所做的内功速率都是相同的：
这种一致性是有限变形理论的基石。基尔霍夫应力与材料应力 之间还有一个非常简单的“前推”(push-forward) 关系，由 给出 [@problem_id:3594867, @problem_id:1549777]。这个方程是现代计算力学的主力。
这些不同的应力度量不仅有不同的公式，它们还有不同的“个性”。让我们看看，如果我们把变形后的物体仅仅进行旋转，而不施加任何额外的拉伸或挤压，它们会如何表现。
柯西应力 和基尔霍夫应力 是空间张量。它们存在于当前构形中，并随物体一起旋转。它们在固定坐标系中的分量会发生变化。
但是第二类皮奥拉-基尔霍夫应力 呢？刚性旋转不会对材料本身增加任何新的变形。它只是改变了材料在空间中的朝向。因此，一个真正代表材料应变状态的应力度量不应受此旋转的影响。而美妙的是，它确实不受影响。第二类皮奥拉-基尔霍夫应力 是客观的 (objective)，或称材料坐标系无关的 (materially frame-indifferent)。如果在刚性旋转前后计算 ，你会得到完全相同的张量。
这个特性使得 成为定义材料内禀行为——其本构律 (constitutive law)——的理想选择。材料的刚度或强度应该取决于它被拉伸了多少，而不是它指向哪个方向。因此，材料模型几乎总是被表述为第二类皮奥拉-基尔霍夫应力 和一种材料应变度量（如格林-拉格朗日应变 ）之间的关系。其他的应力度量，如 和 ，则可以通过变换规则得到。例如，在只有纯拉伸的“无旋”变形中，、 和变形张量 之间的联系变得尤为清晰。
即使在一个看似简单的情况下，比如一个浸没在均匀压力 的流体中的物体，其柯西应力仅为 ，材料应力 也会呈现出更复杂的形式：。这个表达式完美地展示了材料应力不仅取决于外部压力，还取决于变形本身的几何形状，这由体积变化 和变形张量 捕捉。
从一个简单的问题——如何在一个变化的物体中测量应力？——我们揭示了一个相互关联的张量家族，每个成员都有其独特的角色和个性。基尔霍夫应力 作为一个关键的计算和概念纽带，将柯西应力的物理现实与皮奥拉-基尔霍夫应力的解析优雅统一起来，揭示了材料物理学中固有的深刻统一与美感。
在我们完成了对变形数学图景的探索之后，你可能会留下一个挥之不去的问题。我们从应力是作用在面积上的力这个直观概念开始——即传感器在一个被压扁的橡胶块中测得的“真实”或柯西应力 。然后，我们引入了一系列全新的角色，其中最主要的是第二类皮奥拉-基尔霍夫应力张量 。为什么要这么复杂？为什么要去发明一种只存在于未变形参考形状的“幽灵”世界中的新应力类型？这仅仅是一个数学技巧吗？
答案或许令人惊讶，这个数学技巧是整个力学领域最深刻和最有用的思想之一。第二类皮奥拉-基尔霍夫应力不仅仅是为了计算上的方便；它是一条金线，将能量的抽象世界、材料行为的物理现实以及现代工程模拟的数字领域紧密地联系在一起。它是无形的建筑师，让我们能够预测、设计和理解变形事物的世界。要领会它的力量，我们必须看它如何工作。
想象一下拉伸一根橡皮筋。你对它做功，它将这些功以势能的形式储存起来。当你放手时，它释放能量并弹回原状。对于一大类被称为“超弹性”(hyperelastic) 固体的材料——不仅包括橡胶，还包括软生物组织和许多合成聚合物——这种行为由一个主方程控制：应变能函数，通常写作 。这个函数就像材料的 DNA；它编码了任何可能变形的能量代价。
美妙之处在于：第二类皮奥拉-基尔霍夫应力 直接从这个能量函数中诞生。简单地说，它就是能量函数对应变的导数。这个基本关系非常优雅：
其中 是右柯西-格林变形张量，它度量了材料的平方拉伸。这个方程极其强大。它告诉我们，如果我们能写出一个材料储存能量的公式，我们就能立即推导出一个完整的“本构律”，预测其在任何变形下的应力。皮奥拉-基尔霍夫应力不仅仅是真实应力的重新映射；它衡量的是材料的内禀结构在拉伸和剪切时，其储存的能量如何变化。
例如，一个简单的橡胶模型，即新胡克模型 (neo-Hookean model)，有一个非常简单的能量函数。利用它，我们可以推导出在任何情况下 的精确表达式，例如橡胶带的大拉伸。更复杂的模型，如穆尼-里夫林模型 (Mooney-Rivlin model)，使用稍微复杂一些的能量函数来更好地匹配真实材料的行为，但原理是相同的：写下能量函数，求导，然后就得到了第二类皮奥拉-基尔霍夫应力。能量和应力之间的这种直接联系，使得 成为描述材料内禀属性的自然语言，这些属性独立于其旋转或当前形状。
当我们进入更复杂材料的世界时，这种基于能量的方法的威力才真正显现出来。想想肌肉纤维、树枝，或飞机机翼中的碳纤维复合材料。这些材料在所有方向上都不相同；它们是各向异性的 (anisotropic)。肌肉沿其长度方向比横向要强得多。我们如何才能在方程中捕捉到这一点？使用皮奥拉-基尔霍夫应力，这出奇地简单。
我们只需让我们的应变能函数 对优先方向敏感。我们可以引入一个新变量，比如 ，来表示沿纤维方向的拉伸量。通过将这个变量包含在我们的能量函数中，我们就在教导模型关于材料的内部结构。
然后，当我们求导以找到 时，我们的应力方程中会神奇地出现一个新项——这一项明确地解释了沿纤维方向拉伸的额外阻力。科学家和工程师就是这样为生物组织建立预测模型的，从而设计出更好的人工心脏瓣膜或理解肌腱损伤。他们也是这样为航空航天和汽车应用设计先进复合材料的。这个框架非常灵活，它允许我们将极其复杂的材料行为编码到一个单一、优雅的势函数中，而 则是其忠实的信使。
所以，我们有了一种预测应力的方法。但是我们如何用它来解决一个现实世界的问题，比如计算碰撞中汽车框架内的力？几何形状过于复杂，无法使用简单的公式。答案是有限元法 (FEM)，现代工程分析的主力。
当变形很大时，相对于物体的原始、未变形形状进行所有计算要方便得多。这种方法被称为全拉格朗日 (Total Lagrangian) 列式。那么，哪种应力度量存在于这个未变形的参考世界中呢？正是我们的朋友，第二类皮奥拉-基尔霍夫应力 。
每个有限元程序必须求解的核心方程是虚功原理的表达式，它平衡了内力和外力。在全拉格朗日列式中，内虚功——即内应力在微小虚位移期间所做的功——使用 来表达最为自然和优雅。它正是与格林-拉格朗日应变（同样定义在参考构形中的应变度量）功共轭的应力度量。由于这种完美的伙伴关系， 位于那些模拟从气球充气到卫星天线展开等一切事物的软件的计算引擎的核心。
假设我们的有限元模拟已经完成。计算机已经在我们的数字模型的数千个点上计算了变形梯度 和第二类皮奥拉-基尔霍夫应力 。现在怎么办？工程师不关心 ；她想知道最终变形部件中的真实应力 ，以判断它是否会断裂。
这就是我们闭合循环的地方。我们使用之前学到的变换，将计算出的应力 从参考构形“前推”回现实世界的空间构形，以找到柯西应力 ：
其中 是体积变化量。这是工程分析中的一个关键步骤。根据计算出的柯西应力，我们随后可以计算像 von Mises 等效应力这样的失效准则，它告诉我们材料最有可能在哪里屈服或断裂。这个过程——在参考系中使用 进行计算，然后前推得到 以进行解释——是计算固体力学的标准工作流程。
这个框架也与实验世界完美地连接起来。材料科学家可以拿一张橡胶片，在机器中拉伸它，并测量施加的力和产生的变形。从这些原始数据中，可以计算出第二类皮奥拉-基尔霍夫应力张量的分量。然后，可以将这些 的实验值与理论模型（如新胡克模型）的预测进行比较，以检验模型的有效性，或确定材料的特定参数。检查测得的应力分量是否具有预期的对称性，甚至可以验证材料各向同性等假设。这就完成了一个强大的发现循环：实验为理论提供信息，理论为计算提供动力，而计算则预测现实世界的行为。
第二类皮奥拉-基尔霍夫应力，起初可能看起来像一个抽象的麻烦，结果却是连接所有这些领域不可或缺的关键。它是材料原子键中储存的能量与我们在实验室中可测量的力之间的桥梁，也是驱动我们最强大预测模拟的引擎。它证明了这样一个事实：有时，最实用的工具是一种美丽的抽象。