客观应力率

在连续介质力学领域，最基本的任务之一是描述材料在变形时其内部应力状态如何演变。尽管这看似直接，但当材料在拉伸或压缩之外还经历大转动时，一个深刻的挑战便出现了。我们如何才能建立普适的、并且不依赖于我们观察者自身在空间中如何旋转的应力物理定律？

本文所要解决的核心问题是，简单直观的应力张量时间导数存在根本性缺陷。它将真实的材料响应与由刚体旋转引起的虚假效应混合在一起，从而导致所建立的定律不具有客观性。为解决此问题，人们提出了​​客观应力率​​的概念——这是一种“更智能”的导数，它提供了一种纯净的、无旋转影响的度量，用以衡量应力在材料自身参考系中的真实变化。

本文将引导您深入理解这一基本概念。在第一章​​原理与机制​​中，我们将探讨客观性的理论基础，通过类比来理解为何朴素的应力率会失效，并考察不同客观率（如 Jaumann 率和 Truesdell 率）的构造及其后果。随后，在​​应用与跨学科联系​​一章中，我们将看到这些理论思想如何在计算机模拟、结构工程乃至复杂流体研究中产生至关重要的现实影响。

想象一下，您正在尝试描述风。您有一个风向标，它是一只立于尖塔之上的华丽金鸡。风吹来时，金鸡随之转动。作为一名勤奋的科学家，您记下指针改变方向的速率。但这时，一个淘气的朋友走过来，开始旋转整个尖塔。现在，指针改变方向的速度快多了！如果您天真地记下这个新的、更快的速率，您就会错误地断定风发生了巨大变化。您的测量结果被仪器的旋转所污染。要找出风的真实变化，您必须足够聪明，从您的测量值中减去尖塔的旋转。

这正是我们在连续介质力学中描述材料应力变化时所面临的困境。简单的应力材料时间导数 $\dot{\boldsymbol{\sigma}}$ 就像对旋转金鸡的朴素测量。它既看到了因材料实际被挤压和拉伸（变形）而产生的应力变化，也看到了仅仅因为材料本身在空间中旋转而发生的应力分量的“变化”。然而，一条物理定律不应依赖于我们决定让实验室旋转多快。它必须独立于观察者，这一原则我们称之为​​材料坐标系无关性​​或​​客观性​​。我们朴素的率 $\dot{\boldsymbol{\sigma}}$ 在这项检验中彻底失败了。

让我们来看看这种“污染”为何会发生。假设一个参考系中的观察者看到的应力张量为 $\boldsymbol{\sigma}$。第二个观察者的参考系相对于第一个参考系以旋转 $\boldsymbol{Q}(t)$ 运动，他将看到一个不同的应力张量 $\boldsymbol{\sigma}^* = \boldsymbol{Q}\boldsymbol{\sigma}\boldsymbol{Q}^{\mathsf{T}}$。这只是旋转坐标系时张量分量变化的标准法则。

当我们问应力的变化率如何变化时，问题就出现了。利用微分的乘法法则，第二个观察者看到的变化率为： $\dot{\boldsymbol{\sigma}}^* = \dot{\boldsymbol{Q}}\boldsymbol{\sigma}\boldsymbol{Q}^{\mathsf{T}} + \boldsymbol{Q}\dot{\boldsymbol{\sigma}}\boldsymbol{Q}^{\mathsf{T}} + \boldsymbol{Q}\boldsymbol{\sigma}\dot{\boldsymbol{Q}}^{\mathsf{T}}$

如果我们将观察者参考系的自旋定义为反对称张量 $\boldsymbol{\Omega} = \dot{\boldsymbol{Q}}\boldsymbol{Q}^{\mathsf{T}}$，经过少量代数运算可以得到： $\dot{\boldsymbol{\sigma}}^* = \boldsymbol{Q}\dot{\boldsymbol{\sigma}}\boldsymbol{Q}^{\mathsf{T}} + \boldsymbol{\Omega}\boldsymbol{\sigma}^* - \boldsymbol{\sigma}^*\boldsymbol{\Omega}$

瞧！变换后的变化率 $\dot{\boldsymbol{\sigma}}^*$ 并不仅仅是原始变化率的旋转版本 $\boldsymbol{Q}\dot{\boldsymbol{\sigma}}\boldsymbol{Q}^{\mathsf{T}}$。它多了两个额外的项 $\boldsymbol{\Omega}\boldsymbol{\sigma}^* - \boldsymbol{\sigma}^*\boldsymbol{\Omega}$，这两项依赖于观察者的自旋 $\boldsymbol{\Omega}$。这就是“旋转金鸡”问题的数学标记。像 $\dot{\boldsymbol{\sigma}} = \mathbb{C}:\boldsymbol{D}$（其中 $\mathbb{C}$ 是弹性张量，$\boldsymbol{D}$ 是变形率）这样的本构律不可能是自然界的基本定律，因为它的形式本身就依赖于观察者。

我们该如何修正这个问题？我们必须发明一种“更智能”的应力率，一种清除了这些旋转效应的应力率。我们需要定义一种新的导数，我们称之为 $\stackrel{\circ}{\boldsymbol{\sigma}}$，它必须是真正​​客观的​​，这意味着它应该像一个真正的张量那样进行干净利落的变换：$\stackrel{\circ}{\boldsymbol{\sigma}}^* = \boldsymbol{Q}\stackrel{\circ}{\boldsymbol{\sigma}}\boldsymbol{Q}^{\mathsf{T}}$。

最直观的方法是在一个随材料一同旋转的坐标系中测量应力率。这就是​​同转率​​（corotational rate）的精髓。我们取朴素的应力率 $\dot{\boldsymbol{\sigma}}$，并加上一些修正项，这些修正项被设计用来精确抵消虚假的旋转部分。其一般形式如下： $\stackrel{\circ}{\boldsymbol{\sigma}} = \dot{\boldsymbol{\sigma}} - \boldsymbol{\omega}\boldsymbol{\sigma} + \boldsymbol{\sigma}\boldsymbol{\omega}$ 此处，$\boldsymbol{\omega}$ 是一个精心选择的反对称张量，代表了我们对材料自旋速率的最佳猜测。

一个非常常见且自然的选择是使用自旋张量 $\boldsymbol{W}$，它就是速度梯度 $\boldsymbol{L}$ 的反对称部分。这便催生了著名的 ​​Zaremba-Jaumann 率​​（或简称 ​​Jaumann 率​​）： $\stackrel{\nabla}{\boldsymbol{\sigma}} = \dot{\boldsymbol{\sigma}} - \boldsymbol{W}\boldsymbol{\sigma} + \boldsymbol{\sigma}\boldsymbol{W}$ 根据其构造，这个新的率是客观的。一个检验其有效性的绝佳测试是考虑一个只进行纯刚体旋转而没有任何变形的物体。在这种情况下，没有拉伸或挤压，因此 $\boldsymbol{D} = \mathbf{0}$。从物理上我们预期，“真实”的应力变化率应为零。对于纯刚体旋转，事实证明 $\dot{\boldsymbol{\sigma}} = \boldsymbol{W}\boldsymbol{\sigma} - \boldsymbol{\sigma}\boldsymbol{W}$。将此代入 Jaumann 率的定义，得到 $\stackrel{\nabla}{\boldsymbol{\sigma}} = (\boldsymbol{W}\boldsymbol{\sigma} - \boldsymbol{\sigma}\boldsymbol{W}) - \boldsymbol{W}\boldsymbol{\sigma} + \boldsymbol{\sigma}\boldsymbol{W} = \mathbf{0}$。它完美地工作了！我们更智能的率正确地报告出，在材料自身的同转坐标系中，应力没有发生任何有意义的变化。这使得我们能够建立客观的本构律：$\stackrel{\nabla}{\boldsymbol{\sigma}} = \mathbb{C}:\boldsymbol{D}$。

但等等。自旋张量 $\boldsymbol{W}$ 是 $\boldsymbol{\omega}$ 的唯一可能选择吗？当然不是！我们打开了一个充满可能性的潘多拉魔盒。物理学家和工程师们定义了五花八门的客观率，每一种都基于对材料自旋的不同定义。

例如，​​Green-Naghdi 率​​ 使用了变形梯度极分解（$\boldsymbol{F} = \boldsymbol{R}\boldsymbol{U}$）中旋转张量 $\boldsymbol{R}$ 的自旋，它代表了材料“纤维”本身的旋转。另一个重要选择是 ​​Truesdell 率​​，它不严格属于同转率，但它是通过考虑未变形参考构形中应力的变化率并将其映射到当前构形推导出来的。它的形式更为复杂：$\overset{\mathrm{T}}{\boldsymbol{\sigma}} = \dot{\boldsymbol{\sigma}} - \boldsymbol{L}\boldsymbol{\sigma} - \boldsymbol{\sigma}\boldsymbol{L}^{\mathsf{T}} + (\text{tr}\boldsymbol{D})\boldsymbol{\sigma}$。

所有这些率（以及其他许多率！）都完全是客观的。它们都能正确处理叠加的刚体运动。但它们并不相同。如果我们将材料置于简单剪切之下，自旋张量 $\boldsymbol{W}$ 与 Green-Naghdi 率中使用的材料纤维自旋 $\boldsymbol{\Omega}$ 是不同的。这就引出了一个引人入胜且至关重要的问题：如果我们使用不同的客观率来构建本构模型，我们会得到不同的物理预测吗？

答案是响亮的“是”。客观率的选择不仅仅是数学上的细微差别，它具有深远的物理后果。为了说明这一点，让我们考虑一个著名的思想实验：一个弹性块承受大的、连续的简单剪切，就像一副扑克牌从顶部被推动一样。

我们使用一个简单的基于率的本构模型，称为​​次弹性​​（hypoelastic）模型，该模型表明客观应力率与变形率成正比：$\stackrel{\circ}{\boldsymbol{\sigma}} = 2G\,\boldsymbol{D}$，其中 $G$ 是剪切模量。

如果我们使用 ​​Truesdell 率​​，我们会发现剪应力 $\sigma_{12}$ 随剪切量 $\gamma$ 线性增长：$\sigma_{12}^{T}(\gamma) = G\gamma$。这似乎是合理的：剪切越大，应力越大。

但如果我们使用 ​​Jaumann 率​​，就会发生完全奇异的事情。剪应力并不会无限增长，而是会振荡：$\sigma_{12}^{J}(\gamma) = G\sin(\gamma)$！该模型预测，当您持续剪切材料时，应力会上升，然后下降，接着变为负值，如此循环。这非常反直觉，并且与大多数真实材料在大剪切下的行为不符。两种预测之间的差异 $\Delta \sigma_{12}(\gamma) = G(\gamma - \sin(\gamma))$ 随着剪切量的增大而变得巨大。

这个结果是一个警示信号。它告诉我们，这些简单的次弹性模型虽然是客观的，但并未捕捉到弹性的全貌。使用 Jaumann 率构建的模型表现出一种奇怪的​​路径依赖性​​。如果您将其剪切 $2\pi$ 的量，然后再反向剪切回到起点，最终应力为零，但您所做的功却不为零。您在一条闭合回路上莫名其妙地产生或损失了能量，而一个真正的弹性材料绝不应该出现这种情况。

这把我们引向了最深邃的思想。描述弹性行为的黄金标准是​​超弹性​​。超弹性材料是指其应力状态可由一个​​储能势函数​​ $\psi$ 导出的材料，该函数仅依赖于当前的变形状态。把它想象成一个完美的弹簧，其势能为 $\frac{1}{2}kx^2$。力是势的导数，$F = kx$，它只取决于当前的位置 $x$，而与如何到达该位置无关。从 $x_1$ 移动到 $x_2$ 再回到 $x_1$ 所做的功恒为零。这是一个保守系统。

变形过程所做的功可以被想象成在某个地形上的旅程。对于超弹性材料，这个地形就像真实世界中受重力作用的山脉。从一点爬到另一点所做的功只取决于高度的变化（势能），而与所走的具体路径无关。

然而，使用 Jaumann 率的次弹性模型就像一座“魔法山”。您可以沿一条闭合回路行走并回到起点，却发现自己获得或失去了能量。所做的功是路径依赖的。在数学上，我们称“应力功率一范式非恰当（not exact）”。这意味着对于一般的次弹性模型，不存在一个底层的储能函数 $\psi$。这就是为什么基于 Jaumann 率的模型被称为次弹性（hypoelastic，“次等”弹性），而非超弹性（hyperelastic，真正的弹性）。

是否存在既是次弹性又是超弹性的模型？是的，但仅在非常特殊的情况下。例如，使用一种称为“对数率”的特定率的模型，在数学上等价于一个其能量是对数（或 Hencky）应变的函数的超弹性模型。

这给了我们一个关键的洞见。如果一种材料是真正的弹性体，我们根本就不应该使用基于率的模型。我们应该定义一个储能函数，例如，定义为右 Cauchy-Green 张量 $W(\boldsymbol{C})$ 的函数。这种​​超弹性​​中使用的方法从一开始就将客观性内建其中，因为张量 $\boldsymbol{C} = \boldsymbol{F}^{\mathsf{T}}\boldsymbol{F}$ 自动对叠加的旋转保持不变。由于应力是直接从当前变形计算得出，而不是通过对应力率进行积分，整个客观应力率的问题就变得没有必要了。客观率是用于非弹性世界的工具——例如塑性和粘塑性——在这些领域，材料的状态本质上依赖于所经历的路径。

在讨论了这么多之后，您可能会想：如果朴素的率 $\dot{\boldsymbol{\sigma}}$ 如此错误，为什么每本入门工程教科书都在使用它？答案在于一个优美的量级分析。

回想一下，朴素的率被一个形如 $\boldsymbol{\omega}\boldsymbol{\sigma} - \boldsymbol{\sigma}\boldsymbol{\omega}$ 的虚假旋转项所污染。让我们比较一下这个项的大小与来自物理应变的“真实”应力率的大小。

在​​无穷小应变范围​​内，所有应变和旋转都非常小（比如，量级为 $\epsilon \ll 1$），此时应力也很小（$\|\boldsymbol{\sigma}\| \sim \epsilon$），自旋速率也很小（$\|\boldsymbol{\omega}\| \sim \dot{\epsilon}$）。因此，虚假的旋转率的量级约为 $\sim \dot{\epsilon} \cdot \epsilon$。而由应变产生的物理率的量级约为 $\sim \dot{\epsilon}$。虚假项与物理项之比的量级为 $\epsilon$。由于 $\epsilon$ 非常小，旋转误差是一个可以忽略不计的高阶项！对于小变形，我们朴素的率是真实客观率的一个完美近似。

然而，在​​有限应变范围​​内，即使应变很小，旋转也可能很大。自旋速率 $\|\boldsymbol{\omega}\|$ 的量级可以为 1。现在，虚假的旋转率的量级约为 $\sim 1 \cdot \epsilon = \epsilon$。这与物理应力率的量级相同！误差不再可以忽略；它成了一个主导效应。

这是整个谜题的最后一块拼图。对于一根轻微弯曲的梁或一座承受交通荷载的桥梁，旋转是微小的，我们可以安全地使用简化的、应力率之间区别消失的小应变理论。但对于一个旋转的汽车轮胎、一块被冲压成复杂形状的金属板，或高分子聚合物的湍流，旋转是巨大的。在这些世界里，客观应力率的概念不是学术上的好奇心——它是编写任何有意义的物理定律的必要关键。

既然我们已经探讨了客观性这一抽象原理，您心中可能还留有一个挥之不去的问题：这种数学形式主义真的重要吗？一种客观率与另一种之间是否存在实质性差异，或者这都只是学术上的记账游戏？答案是响亮的“是”。选对率往往是区分一个忠实预测真实世界行为的计算机模拟与一个产生纯粹胡言乱语的模拟的关键。

现在，让我们踏上一段旅程，亲眼见证这一原理的实际作用。我们将深入数字碰撞测试的核心，目睹材料在极端应力下的奇特行为，并揭示其与流动聚合物世界之间惊人的联系。我们将发现，这一个单一而优雅的思想，如同一条线索，贯穿于广阔多样的科学和工程领域，确保我们的虚拟世界与物理现实使用同一种语言。

想象一下，为了设计从现代飞机的轻质框架到家用车防撞性能的方方面面，所动用的巨大计算能力。这些模拟的核心是诸如有限元法（FEM）和物质点法（MPM）之类的数值方法，它们将复杂物体分解成由更简单、相互连接的单元组成的拼图。模拟的任务是计算出整个组件在荷载作用下如何变形、受力以及可能发生的失效。

整体行为由宏大的物理定律主导，如虚功原理，它平衡了内力与外力。但模拟的真正灵魂存在于每一个微小单元内部：即*本构模型*。这是一套规则——材料的“DNA”——决定了应力如何响应变形。而正是在这里，在应力状态的瞬时更新中，我们的客观率得以展现其威力。

如我们所学，我们不能简单地使用 Cauchy 应力 $\dot{\boldsymbol{\sigma}}$ 的时间导数。我们需要一个客观率。但该用哪一个呢？连续介质力学的世界提供了一整套选择。最著名的两个是 Zaremba-Jaumann 率（一个历史悠久的主力模型）和 Green-Naghdi 率（一个更复杂的选择）。根本区别在于它们选择与什么“同转”。Jaumann 率随材料的局部流体状涡量 $\boldsymbol{W}$ 一起旋转，而 Green-Naghdi 率则随从变形本身中提取的真实材料旋转 $\boldsymbol{R}$ 一起旋转。

这个选择重要吗？思考一个警示性的案例：大振幅循环剪切。想象一下，你拿一块金属反复地来回剪切，每次都回到起点。直观地看，如果我们在零点附近对称地循环应变，应力响应也应该稳定成一个对称的回线。然而，如果我们使用 Jaumann 率运行这个模拟，奇怪的事情发生了。随着每个循环，平均应力开始漂移，或“棘轮”式地上升到越来越高的值，就好像材料被逐渐地向一个方向推动。这就像在平地上来回摇晃一辆汽车，却看到它奇迹般地自己向前蠕动。这在物理上是错误的。

Green-Naghdi 率通过将其参考系与材料的真实旋转正确对齐，预测不会出现这种漂移。它显示出稳定、对称的应力-应变回线，与实验结果一致。这个简单模型的戏剧性失败揭示了一个深刻的真理：Jaumann 率虽然是客观的，但它并不能“积分”得到一个真正的弹性势。这意味着对于纯弹性材料，它可能预测能量在闭合变形循环中无中生有——这是另一个物理上的不可能！更先进的公式，例如基于对数应变的公式，从一开始就建立在坚实的超弹性储能势基础上，从而优雅地回避了这些问题。率的选择不仅仅是品味问题，它关乎您模拟的物理完整性。

客观性原理的应用远不止于简单的弹塑性模型。世界上充满了行为远为复杂的材料。

考虑那些表现出“记忆”效应的材料。当你弯曲一个回形针时，它不仅变形，而且变得更难进一步弯曲。如果你把它弯回去，它会记住初始加载的方向。在塑性理论中，这通过一个称为随动硬化的概念来建模，该概念引入了一个内部“背应力”张量（通常表示为 $\boldsymbol{\alpha}$），用于在应力空间中追踪屈服面的中心。就像 Cauchy 应力 $\boldsymbol{\sigma}$ 一样，这个背应力张量是一个存在于变形材料中的物理量。因此，为了在大转动下正确模拟材料的记忆效应，背应力本身必须使用客观率进行更新。如果忘了这样做，就好像试图用一个未考虑船只旋转的罗盘在旋转的船上导航一样。

那么极端事件呢，比如高速车祸或射弹撞击装甲？在这些情况下，材料的响应不仅取决于应变的大小，还取决于应变率和温度。像 Johnson-Cook 定律这样的模型就是为这些场景设计的。它们在一个标量流动应力函数中捕捉了这种复杂性。由于硬化的物理机制被打包进标量不变量（自然独立于旋转的量）中，这些模型天生就与任何有效的客观应力率兼容。客观性的基本要求依然存在，但硬化的具体公式并不偏好某一种客观率。

选择正确（或错误）的率所带来的最戏剧性的后果，或许出现在结构稳定性领域。工程师不仅需要预测桥梁在交通荷载下会如何弯曲，还需要预测它在哪个点上可能会灾难性地屈曲或开始失控地颤振，就像臭名昭著的 Tacoma Narrows 大桥一样。

这些生死攸关的问题通过分析结构的切向刚度矩阵 $K_T$ 来回答，该矩阵代表了结构对变形的瞬时抵抗能力。结构的稳定性被编码在该矩阵的特征值中。$K_T$ 的一个关键属性是其对称性。事实证明，模拟中客观应力率的选择对这种对称性有着直接而深远的影响。

基于真实能量势（如超弹性）的公式总能产生对称的 $K_T$。然而，基于 Jaumann 率的次弹性模型通常会产生一个非对称的 $K_T$。Jaumann 率中的自旋项将非对称分量引入了刚度矩阵的数学结构本身。

为什么非对称矩阵如此令人不安？对称系统只能发生屈曲（一种静态失稳）。而非对称系统不仅可以屈曲，还可能“颤振”——一种振动随时间指数增长的动态失稳。一个假设矩阵对称的标准稳定性分析将完全无法察觉这种可能性。一个使用基于 Jaumann 率的模型而没有配备适当的非对称分析工具的工程师，可能会将一个设计认证为稳定，而实际上，它是一个随时可能发生的颤振灾难。再一次，本构模型中一个看似微小的选择，对于预测大型工程系统的安全性和可靠性具有深远的影响。

到目前为止，我们的旅程一直停留在固体领域。但一个基本原理的美妙之处在于其普适性。客观性的概念是否也出现在其他地方？

让我们进入非牛顿流体力学的世界——研究熔融聚合物、油漆和生物流体等奇特流体的学科。当这些复杂流体流动时，它们会拉伸和旋转，产生内应力，很像固体。一个试图模拟塑料流入模具的流体动力学家，面临着与一个模拟钢梁的固体力学家完全相同的问题：联系应力与变形率的本构方程必须是坐标系无关的。

于是，流体动力学家们独立地发展出了他们自己的一套客观导数！你会遇到诸如​​上、下随钻 Oldroyd 导数​​和 ​​Gordon-Schowalter 导数​​等名称。虽然名称和物理背景不同，但它们是我们所研究的客观率在流体力学中的表亲。事实上，可以证明它们在数学上是相关的，共同构成了一个连续的客观率谱系。

这一发现是物理学统一性的一个绝佳例子。同样深刻的原理——自然法则不依赖于观察者——对模拟截然不同物理系统的科学家们提出了相同的数学挑战。无论您是在模拟汽车碰撞、桥梁屈曲，还是聚合物挤出，对客观率的需求都是一个共同的、统一的线索。它是确保我们的计算模型不仅仅是在解方程，而是在真正使用物理宇宙的语言的基本数学工具。