材料框架无关性：连续介质力学中的客观性原理

玻尔百科

核心要点

材料框架无关性，或称客观性，是一项基本原理，它指出材料的本构律必须独立于观察者的刚体运动。
为满足此原理，应变和能量等物理量必须通过客观度量来描述，例如右柯西-格林张量 ( $C$ )，它能滤除与观察者相关的转动。
应力的标准时间导数不具有客观性，这要求在经历大转动或流动的材料模型中，必须使用特殊构造的客观应力率（例如Jaumann率）。
此原理是工程领域精确计算模拟的一项关键要求，它能防止对非物理应力的预测，并确保结果独立于物体的刚体运动。

引言

一种材料如何表现——它如何拉伸、弯曲或流动——是其固有的属性，不能依赖于观察者是谁，或者他们如何移动。这个看似不言自明的想法，构成了物理学和工程学中一个严谨且至关重要的原理的基础，即材料框架无关性（或称客观性）。尽管这个原理很直观，但在数学上强制执行它却是一个深刻的挑战，并塑造了我们对材料建模的整个理解。若没有它，我们的物理定律将会产生荒谬的结果，比如预测一个刚性旋转的物体会自发产生应力，或者弹簧中储存的能量取决于你的观察视角。

本文旨在探讨客观性在建立具有物理意义的材料行为描述中的关键作用。它解决了如何将真实的材料变形与由观察者的运动和朝向引起的表观变化分离开来的问题。在接下来的章节中，我们将首先踏上“原理与机制”的探索之旅，探讨描述变形的数学语言，寻求应变和应力的客观度量，以及客观性与储存能之间的深层联系。随后，在“应用与跨学科联系”中，我们将看到这一抽象原理如何成为一种实用且不可或缺的工具，指导着弹性、塑性和损伤模型的发展，并成为可靠工程模拟的基石。

原理与机制

想象你正在观察一根被拉伸的橡皮筋。你能感觉到其中的张力，看到它伸长了多少。现在，如果你转过头，或者从它旁边走过，会发生什么？橡皮筋里的张力会改变吗？当然不会。橡皮筋的物理状态——其内部的应力和应变——与你自身的运动或朝向无关。你的描述可能会改变——原来的“左”现在变成了“右”——但材料的真实情况保持不变。这个简单、近乎不言自明的想法，是连续介质物理学中一个深刻而强大原理的核心：材料框架无关性，亦称客观性原理。

从本质上讲，这是适用于材料的相对性原理。它规定，本构律——即支配材料行为的法则——对于所有非加速运动的观察者必须是相同的。我们尤其关注那些通过刚体运动（即恒定平移和旋转的组合）相互关联的观察者。材料内部的物理定律不应依赖于我们选择观察它的“框架”。让我们踏上征途，看这一个简单的思想如何塑造我们对材料行为的整个理解，从固体中储存的能量到驱动现代工程模拟的方程。

场景变换：观察的数学描述

为了使我们直观的想法变得精确，我们首先需要一种语言来描述变形。在力学中，我们使用一个称为变形梯度的数学对象，用张量 $F$ 表示。你可以把它想象成一本小小的“说明书”，它告诉原始未变形材料中的每一个微小矢量如何拉伸和旋转，以达到其在变形体中的新朝向。

现在，让我们引入第二个观察者。想象你静静地站在地面上，而你的朋友（我们称之为“星标”观察者）正乘坐一个旋转的木马从旁飞驰而过。你们都看着同一个变形的物体。当你看到物体处在某个位置和朝向时，你的朋友从她移动、旋转的视角看到了它。她所见的变形形状只是你所见形状的一个旋转版本。如果我们用一个数学上的旋转张量 $Q(t)$ 来描述这个旋转，她测量的变形梯度（我们称之为 $F^\star$ ）与你的测量结果有何关系呢？

通过简单的微积分法则，我们可以证明她的测量值通过左乘与你的值相关联：

$F^\star = Q F$

这是基本的变换法则。变形梯度本身不是客观的；它的分量取决于观察者的转动。材料框架无关性原理现在要求，任何描述材料内在状态的真实物理定律或物理量，其表述方式都必须独立于这个任意的转动 $Q$ 。该定律产生的标量值必须完全不变，或者张量值必须以一种一致且可预测的方式变换。

寻找客观量：应变的客观度量

这引出了一项宏大的探寻：我们可以从 $F$ 构作出的哪些数学量是真正客观的，能够滤除观察者的旋转视角，而只捕捉材料的内在变形？

让我们测试几个候选者。变形梯度 $F$ 本身被排除了；我们刚刚看到它被 $Q$ “污染”了。那么一个更简单的度量，比如其对角元素之和，即迹 $\text{tr}(F)$ 呢？让我们来做一个快速的思想实验。如果材料未变形， $F$ 是单位矩阵 $I$ ，其迹为3。如果一个星标观察者只是绕着某个轴旋转180度，她的旋转由一个张量 $Q$ 描述。她观察到的变形梯度变为 $F^\star = QI = Q$ 。这个 $Q$ 的迹是-1。由于 $3 \neq -1$ ， $\text{tr}(F)$ 显然不是客观的。它是一种错觉，是观察者视角的属性，而非材料的真实情况。

我们需要更巧妙一些。我们需要一个 $F$ 的组合，能奇迹般地让观察者的转动 $Q$ 消失。让我们尝试将 $F$ 与其自身的转置 $F^T$ 相乘，创建一个新张量，即右柯西-格林变形张量， $C = F^T F$ 。现在让我们看看我们的星标观察者对这个量测量到什么：

$C^\star = (F^\star)^T F^\star = (QF)^T (QF) = F^T Q^T Q F$

奇迹在此发生！由于 $Q$ 代表纯转动，其转置 $Q^T$ 就是其逆，意味着 $Q^T Q = I$ ，即单位张量。这个方程优美地简化为：

$C^\star = F^T I F = F^T F = C$

我们发现了！右柯西-格林张量 $C$ 对两个观察者来说完全相同。它是一个客观张量。它成功地滤除了观察者的任意转动，捕捉了材料纯粹的、内在的拉伸。因此，任何由 $C$ 构建的标量，例如其迹 $I_1 = \text{tr}(C)$ 或其行列式 $I_3 = \det(C)$ ，都是一个客观标量，是材料状态的真实度量。例如， $C$ 的行列式与材料的体积变化有关。具体来说，当前体积与初始体积之比 $J = \det(F)$ ，可以通过 $C$ 求得，即 $J = \sqrt{\det(C)}$ ，这意味着体积变化也是一个客观量，理应如此。

材料的本质：能量与极分解

这种对客观性的探求带来了真正深远的影响，尤其是当我们考虑储存在变形弹性材料中的能量时。这个应变能密度，我们称之为 $\Psi$ ，是一个真实的物理量。一根被拉伸的弹簧中锁定的能量，绝不可能取决于你是站着不动还是在旋转中观察它。因此，从变形 $F$ 计算此能量的函数 $\Psi$ 必须是客观的。从数学上讲，这意味着：

$\Psi(F) = \Psi(F^\star) = \Psi(QF)$

要理解这个方程中蕴含的美妙物理，我们可以借助一个非凡的数学成果，称为极分解。该定理告诉我们，任何变形 $F$ 都可以唯一地分解为一个纯拉伸 $U$ 后跟一个纯转动 $R$ 。所以， $F = R U$ 。你可以想象成首先沿着主轴拉伸材料（这是 $U$ ），然后将其刚性旋转到最终的朝向（这是 $R$ ）。张量 $U$ 是对称正定的，被称为右拉伸张量。

现在，让我们回到我们的客观性条件 $\Psi(F) = \Psi(QF)$ 。使用极分解，我们有 $\Psi(RU) = \Psi(Q(RU))$ 。这对任何观察者转动 $Q$ 都必须成立。让我们巧妙地选择一个特定的观察者：一个与材料自身转动的逆转动同步旋转的观察者，即 $Q=R^T$ 。代入此式得到：

$\Psi(F) = \Psi(RU) = \Psi(R^T R U) = \Psi(I U) = \Psi(U)$

这是一个惊人的结论。客观材料的储存能只能依赖于变形的拉伸部分 $U$ ；它完全独立于转动部分 $R$ 。材料不会仅仅因为在空间中进行刚性旋转而储存能量。这就是我们最初的直觉得到了数学上的严谨证实！

那么，这个客观拉伸 $U$ 和我们的客观柯西-格林张量 $C$ 之间有什么关系呢？很简单： $C = F^T F = (RU)^T(RU) = U^T R^T R U = U^T U = U^2$ 。由于 $U$ 是 $C$ 唯一的正定平方根，它们之间有着密不可分的联系。说能量只依赖于 $U$ 等同于说它只依赖于 $C$ 。这是力学的一个基石：对于任何客观的超弹性材料，无论是橡胶、钢还是单晶，其能量的本构律必须能表达为右柯西-格林张量的函数，即 $\Psi(F) = \hat{\Psi}(C)$ 。

力的语言：客观应力及其率

到目前为止，我们已经找到了描述应变和能量的客观方法。那么应力，即内力的度量呢？我们熟悉的柯西应力张量 $\sigma$ 是在当前变形构形中定义的。如果我们的星标观察者将她的坐标系旋转 $Q$ ，她测量的应力分量 $\sigma^\star$ 将是我们应力的旋转分量：

$\sigma^\star = Q \sigma Q^T$

这是一个客观空间张量的标准、一致的变换法则。因此，任何从变形计算柯西应力的本构律，比如 $\sigma = \hat{\sigma}(F)$ ，必须遵守这个变换法则才能被认为是物理上有效的。也就是说，它必须满足 $\hat{\sigma}(QF) = Q \hat{\sigma}(F) Q^T$ 。

当我们考虑应力不仅取决于当前变形，还取决于变形率的材料时，事情变得更加有趣。想象一下蜂蜜、湿水泥，或是在高温下锻造的金属。为了对这些材料进行建模，我们需要讨论应力的变化率。

我们的第一直觉是简单地取柯西应力的时间导数 $\dot{\sigma}$ 。让我们测试一下这是否是一个客观量。通过对 $\sigma$ 的变换法则对时间求导，我们发现情况有点乱：

$\dot{\sigma}^\star = \dot{Q} \sigma Q^T + Q \dot{\sigma} Q^T + Q \sigma \dot{Q}^T$

很明显， $\dot{\sigma}^\star$ 不等于 $Q \dot{\sigma} Q^T$ ，因为多了涉及 $\dot{Q}$ （观察者坐标系的角速度）的项。应力的简单物质时间导数不是客观的！它的值被观察者的旋转所污染，使其作为本构变量在物理上没有意义。

让我们看看如果我们使用一个非客观的率，情况会变得多么糟糕。考虑一个在无应力环境中已经处于恒定内应力（比如沿x轴的拉伸）下的块体，我们只是对其进行刚性旋转。物理上，应力状态只是随块体一起旋转。对于这种纯刚性旋转，变形率张量 $D$ 为零。一个像 $\dot{\sigma} = \mathbb{C}:D$ （其中 $\mathbb{C}$ 是弹性张量）这样的朴素定律会预测 $\dot{\sigma}=0$ 。这意味着在我们实验室坐标系中的应力分量将保持不变。但这在物理上是荒谬的！这将意味着随着块体旋转，x方向的拉伸应力会奇迹般地消失，并出现虚假的新应力来保持张量分量恒定。这是模型的灾难性失败，是直接违反客观性原理的结果。

解决方案是发明新类型的时间导数，称为客观应力率。它们被专门构造成客观的。它们通常以非客观的物质导数 $\dot{\sigma}$ 开始，并添加涉及应力 $\sigma$ 和一个自旋张量 $\mathcal{W}$ 的校正项，以精确抵消依赖于观察者的部分。最终得到的率，通常表示为 $\overset{\triangledown}{\sigma}$ ，会正确地变换为 $\overset{\triangledown}{\sigma}^\star = Q \overset{\triangledown}{\sigma} Q^T$ 。这些通常被称为协同转动率，因为它们测量的是在与材料一同旋转的假想观察者看来应力的变化率。

这里是最后的转折：没有唯一的方法来定义“材料的自旋”。我们是使用来自速度梯度的自旋 $W$ 吗？这得到了著名的Jaumann率。我们是基于极分解中的旋转 $R$ 来定义自旋，即 $\Omega = \dot{R}R^T$ 吗？这得到了Green-Naghdi率。还是我们从基本的材料框架量推导这个率？这引出了Truesdell率。所有这些率，以及许多其他的率，都是完全客观的。然而，它们并不相同。对于相同的运动，它们可以给出不同的值，这一事实引发了力学领域数十年的研究和辩论。客观率的选择本身也成为本构模型的一部分。

从抽象原理到具体模拟

从一个直观的想法到一个关于导数的微妙辩论的旅程可能看起来很抽象，但其后果在科学和工程的日常中无处不在。当工程师设计汽车并需要模拟底盘在碰撞中如何褶皱时，他们使用基于有限元方法（FEM）的强大软件。这种软件必须将材料框架无关性原理融入其核心。

这通常通过协同转动格式来实现。在这些方法中，被模拟物体的每个微小单元都被赋予其自身的局部坐标系，该坐标系随那块材料平移和旋转。所有应变和应力的计算都在这个局部的、旋转的框架中进行，自然地使用了独立于单元整体刚体运动的客观度量。这保证了模拟结果——关于一个部件是会弯曲还是会断裂的预测——在物理上是真实的，并且不依赖于模拟的汽车是向北还是向西行驶，或者它在撞击前是否在空中翻滚。

所以，下次你看到一个壮观的碰撞模拟或一个旋转粘性流体的模型时，你可以欣赏到在幕后默默工作的优雅原理。这一切都回归到一个简单、优美且严苛的要求：物理定律不关心你从哪个方向看。

应用与跨学科联系

你可能会认为，像“材料框架无关性”这样抽象的原理是只有理论物理学家才关心的东西——一个适合在尘土飞扬的讲堂里辩论的哲学细节。事实远非如此。这个原理不仅仅是哲学上的讲究；它是一个务实的、实践性的设计规范，适用于任何旨在描述真实世界的理论。它是矗立在可靠物理模型与无稽预测之间的沉默守护者。忽略它就如同在沙上建屋。

让我们一起踏上旅程，看看这个简单的想法——支配材料的物理定律不应该在乎我们是静止站立还是在旋转木马上旋转——是如何塑造了整个科学和工程领域的。我们将看到它如何引导我们，纠正我们，并最终带领我们对材料世界有更深刻、更统一的理解。

拉伸弹簧的迷惑性简洁

让我们从最简单的地方开始：熟悉的线性弹性世界。当你稍微拉伸一根弹簧或弯曲一把尺子时，它的行为由一个极其简单的理论来描述。这个理论的基石是小应变张量 $\boldsymbol{\varepsilon}$ 。现在，这个数学对象有一个相当神奇的特性。在微小变形和微小转动的领域里，它的构造巧妙地使其几乎完全“无视”转动。如果你取一小块材料并将其旋转一个微小的角度，小应变张量几乎记录不到任何变化。

正因如此，如果你将整个弹性理论建立在小应变张量之上——例如，通过将储存能 $\psi$ 定义为其函数——你就能免费获得框架无关性！这个理论“通过构造即为客观的”，甚至无需你刻意为之。就好像宇宙在简单情况下给了你一张通行证。这就是为什么经典的小应变弹性力学如此稳健和成功：它的根基尊重了自然的这一基本对称性，即使我们并不总是有意识地去思考它。

问题变得复杂：大转动情形

但世界并不总是那么简单。当一个物体旋转很大，即使它没有怎么拉伸或弯曲时，会发生什么？想想汽车的传动轴、涡轮叶片或体操运动员的身体。在这里，小应变张量的“魔力”消失了。它被大转动所污染，再也不能被信任来只报告真实的变形。

此时，材料框架无关性原理从一个被动的观察者变成了一个积极的指导者。它告诉我们：“你旧的应变度量已经坏了！你必须发明一个新的，对转动免疫的度量。” 这一要求导致了*有限应变度量*的发展，例如格林-拉格朗日应变张量 $\boldsymbol{E}$ 。该张量由变形梯度 $\boldsymbol{F}$ 定义为 $\boldsymbol{E} = \frac{1}{2}(\boldsymbol{F}^{\mathsf{T}}\boldsymbol{F} - \boldsymbol{I})$ 。如果你推导一下数学，你会发现在刚性转动下， $\boldsymbol{E}$ 保持完全不变。它是一个真正客观的变形度量。

通过将我们的超弹性理论建立在依赖于 $\boldsymbol{E}$ （或其近亲，右柯西-格林张量 $\boldsymbol{C}$ ）的储存能函数上，我们可以为像橡胶这样的材料创建模型，这些材料可以承受巨大的拉伸和扭曲，并且完全有信心我们的理论不会被物体的整体转动所蒙骗。这是我们的第一个重要线索：该原理迫使我们极其谨慎地选择我们的数学工具。

率的危机：塑性与蠕变的挑战

当我们考虑那些会“流动”的材料，比如正在锻造的金属块或沿山坡蠕动的冰川时，情节变得更加有趣。对于这些材料，我们通常不仅关心当前的应力状态，还关心它变化的速度。我们需要一个“率形式”的本构律。

最显而易见的方法是使用我们熟知的柯西应力的物质时间导数 $\dot{\boldsymbol{\sigma}}$ 。但在这里我们遇到了一个可怕的问题。要理解这一点，想象你正站在一个旋转的木马上。你相对于地面的速度是你相对于平台的行走速度和平台旋转速度的组合。简单的时间导数 $\dot{\boldsymbol{\sigma}}$ 就像你相对于地面的速度——它无可救药地将源于材料应变的“真实”应力变化率与纯粹由物体刚性转动引起的“虚假”变化率混合在一起。

一个简单的标度分析表明，对于微小转动，这个虚假的转动部分是一个可以安全忽略的高阶项，这就是为什么我们在小应变世界里能够侥幸过关。但对于大转动，这个虚假的部分和真实的部分一样大！基于 $\dot{\boldsymbol{\sigma}}$ 的本构律将以一种惊人的方式违反框架无关性。它会预测一个旋转的、无应力的物体会自发产生应力，这简直是荒谬。

解决方案？如果你的导数被污染了，那就发明一个更好的！这场危机催生了*客观应力率*的发明。其中最著名的是Jaumann率， $\boldsymbol{\sigma}^{\nabla J} = \dot{\boldsymbol{\sigma}} - \boldsymbol{W}\boldsymbol{\sigma} + \boldsymbol{\sigma}\boldsymbol{W}$ ，其中 $\boldsymbol{W}$ 是自旋张量（转动速率）。这个新的导数巧妙地减去了虚假的转动贡献，只留下由变形引起的“真实”应力变化率。这类似于在旋转木马上计算你真实的行走速度，方法是用你的地面速度减去由平台旋转引起的速度。这个由框架无关性原理“逼迫”我们得出的优美数学创造，现已成为塑性和蠕变计算模型的基石。

对优雅的追求：超弹性与亚弹性

客观应力率是一个绝妙的补丁。但事实证明，它仅仅是一个补丁。定义客观率有许多不同的方法，而且令人不安的是，对于同一个问题，它们可以给出不同的预测。有些，比如Jaumann率，甚至可以在某些情况下预测出非物理行为，比如在简单剪切流中产生振荡应力。这种“率的模糊性”表明，也许我们走在了错误的道路上。

材料框架无关性指向一个更优雅、更深刻的解决方案。与其从一个率方程开始，然后“修复”它以使其客观（一种称为*亚弹性的框架），为什么不从一个从一开始就客观的基础出发呢？这就是基于超弹性的塑性*的哲学。

在这个框架中，响应的弹性部分是从一个标量势能函数 $\psi$ 推导出来的。正如我们所学到的，如果我们将这个能量函数构造成一个客观应变度量的函数，那么整个公式就通过构造而变得客观。不再需要特殊的、非唯一的客观应力率。应力只是当前（客观）应变状态的一个明确定义的函数。这种方法不仅更稳健、物理上更一致，而且在热力学上也是自洽的，因为它保证了弹性变形是能量守恒的。在这里我们看到，该原理不仅仅是错误的纠正者，更是指向更深刻、更优美物理理论的指路牌。

机器中的幽灵：计算中非客观性的危害

那么在现实世界中，如果一个编写价值数百万美元工程软件包的程序员忽略了这个原理，会发生什么？对模拟而言，其后果可能是灾难性的。

想象一下使用一个有限元程序来分析飞行中的直升机叶片。叶片正在快速旋转。如果程序员天真地使用了一个非客观的应变度量——比如说，一个基于在固定坐标系上投影的度量——代码就会被愚弄。它会“看到”叶片在旋转时不断地拉伸和压缩，即使它实际上是完全刚性的。这将产生巨大的、完全虚假的内力。模拟可能会预测叶片即将在巨大的应力下失效，或者它正在无中生有地产生热量，导致一场灾难性的（但幸运的是，只是虚拟的）失效。

这个问题不仅仅是学术上的好奇；它是任何处理大转动和大变形的计算方法的基本合理性检查，从有限元方法（FEM）到物质点法（MPM）。确保底层的本构更新是客观的，是防止模拟被这些虚假的、会产生能量的幽灵所困扰的唯一方法。

材料的通用蓝图

材料框架无关性的影响远不止于简单的力学。它为描述几乎任何可以想象的材料的构成提供了一个通用的蓝图。

损伤材料： 我们如何为一个因微裂纹扩展而弱化的材料建模？在连续介质损伤力学中，我们引入内变量来表示损伤状态。框架无关性原理规定了这些变量必须如何表现。一个简单的标量损伤变量必须对转动保持不变。一个可以描述定向裂纹的二阶损伤张量，必须像应力张量一样变换，以确保系统的总能量是客观的。
各向异性材料： 对于具有內建方向性的材料，如木材、纤维增强复合材料或肌肉组织，又该如何？框架无关性告诉我们如何构建它们的本构律。我们引入一个“结构张量”（如 $\boldsymbol{M} = \boldsymbol{a}_0 \otimes \boldsymbol{a}_0$ ，其中 $\boldsymbol{a}_0$ 是沿纤维方向的矢量），并由该张量和应变张量构成的标量不变量来建立我们的自由能函数。这确保了模型能正确捕捉材料的各向异性，同时保持完全的客观性。

从金属和聚合物到复合材料和生物组织，该原理为描述物质的物理现实提供了一个统一、强大且不可或缺的框架。它是一个美丽的例子，展示了一个简单的对称性原理能够为一个广阔而复杂的世界带来秩序和清晰。