刚体运动学

我们如何描述一个物体的运动？虽然我们可以轻易地追踪其在空间中的路径，但当物体本身改变形状时，一个更深层次的问题就出现了。我们如何将简单的旋转或平移（即刚体运动）与导致材料拉伸和扭曲的真实变形区分开来？发展一种不被观察视角所迷惑的数学语言，是力学中的一个根本挑战。一个正确的变形描述必须独立于观察者自身的运动，这一原则确保了我们的物理定律具有普适性。

本文探讨了刚体运动运动学的核心概念及其在定义变形中的深远作用。在第一章“原理与机制”中，我们将为这项任务构建数学框架。我们将介绍变形梯度，发现为何某些应变度量优于其他度量，并将关键的客观性概念形式化。在这一理论基础之后，第二章“应用与跨学科联系”将展示这些原理如何不仅仅是抽象概念，而是现代科学和工程中的重要工具，从验证物理定律到构建从桥梁到分子等万物的可靠计算模拟。

在我们理解世界的旅程中，我们通常从简单地观察事物运动开始。但对物理学家而言，描述本身是不够的；我们寻求的是其潜在的原理，即支配混乱运动的优雅规则。我们如何建立一种语言来讨论固体的扭转、拉伸和翻滚？更深刻的是，我们如何能确定我们的描述捕捉到了变形的真实物理现实，而不仅仅是我们自己特定的观察视角？这正是我们故事的真正起点——寻求一种既强大又真实的运动描述。

让我们想象一块黏土。在使其变形之前，我们可以用其位置矢量（我们称之为 $\boldsymbol{X}$）来标记黏土中的每一个微小粒子，这处在​​参考构型​​中。这就像在黏土上绘制一个完美的、未变形的坐标网格。现在，我们挤压并扭转它。原来位于 $\boldsymbol{X}$ 的粒子现在处在一个新位置 $\boldsymbol{x}$，这处在​​当前构型​​中。这个运动是一个映射，$\boldsymbol{x} = \boldsymbol{\varphi}(\boldsymbol{X}, t)$，它告诉我们每个粒子去了哪里。

但这个映射告诉我们的是整个物体的命运。如果我们想知道局部，就在某一个点周围，发生了什么呢？考虑一个无穷小的矢量，就像在我们的参考网格上画的一个箭头，我们称之为 $d\boldsymbol{X}$。变形后，这个小箭头变成了一个新的矢量 $d\boldsymbol{x}$——很可能指向一个新的方向并具有新的长度。原始箭头和新箭头之间的关系是关键。对于一个平滑的运动，这种关系是线性的，并由一个强大而单一的实体所捕捉：​​变形梯度​​ $\boldsymbol{F}$。

$d\boldsymbol{x} = \boldsymbol{F} d\boldsymbol{X}$

变形梯度 $\boldsymbol{F} = \frac{\partial \boldsymbol{x}}{\partial \boldsymbol{X}}$ 是一个张量，在我们故事中扮演主角。它是一个局部的“放大镜”，告诉我们关于材料在某一点上变化的一切信息。它包含了关于局部拉伸和局部旋转的所有信息。如果你想象在参考构型中有一个无穷小的材料立方体，$\boldsymbol{F}$ 就是将其转换为在当前构型中它所变成的无穷小平行六面体的算子。为避免物质相互穿透，这个映射必须是可逆的，这意味着 $\boldsymbol{F}$ 的行列式（称为​​雅可比行列式​​ $J=\det\boldsymbol{F}$）必须为正。这个雅可比行列式告诉我们体积如何变化，因为 $dV = J dV_{0}$。

我们可以通过想象一个由刚性旋转和一些附加位移组合而成的运动来感受这一点。如果运动由 $\boldsymbol{x} = \boldsymbol{R}\boldsymbol{X} + \boldsymbol{u}(\boldsymbol{X})$ 给出，其中 $\boldsymbol{R}$ 是一个常数旋转矩阵，$\boldsymbol{u}(\boldsymbol{X})$ 是一个位移场，那么变形梯度就是 $\boldsymbol{F} = \boldsymbol{R} + \nabla\boldsymbol{u}$。这似乎暗示了旋转和变形之间存在一种清晰的加性分离。但自然界的故事要微妙一些，而且正如我们将看到的，要优雅得多。

我们所说的“真实变形”是什么意思？从物理角度看，变形是导致材料产生应力的原因。如果你拿一根钢梁并弯曲它，它会产生内应力。但如果你只是在空中旋转同一根钢梁而不弯曲它，则不会出现新的应力。刚体运动——纯粹的平移或旋转——本身并不构成变形。因此，任何真实的应变度量都必须完全“无视”刚体运动。

让我们来检验这个想法。考虑一个只经历刚性旋转的物体，其运动由 $\boldsymbol{x} = \boldsymbol{R}\boldsymbol{X}$ 描述，其中 $\boldsymbol{R}$ 是一个旋转张量。对于这个运动，变形梯度就是 $\boldsymbol{F} = \boldsymbol{R}$。我们如何从 $\boldsymbol{F}$ 构建一个忽略这种旋转的量？

让我们不看矢量的变化，而是看其长度的平方的变化——这是一个不依赖于方向的标量。我们无穷小矢量 $d\boldsymbol{X}$ 的长度平方是 $d\boldsymbol{X} \cdot d\boldsymbol{X} = d\boldsymbol{X}^{\mathsf{T}} d\boldsymbol{X}$。变形后，其新的长度平方是 $d\boldsymbol{x} \cdot d\boldsymbol{x} = d\boldsymbol{x}^{\mathsf{T}} d\boldsymbol{x}$。代入我们的基本关系 $d\boldsymbol{x} = \boldsymbol{F} d\boldsymbol{X}$，我们得到：

看看我们发现了什么！长度平方的变化完全由中间的张量 $\boldsymbol{C} = \boldsymbol{F}^{\mathsf{T}}\boldsymbol{F}$ 所控制。这就是著名的​​右 Cauchy-Green 变形张量​​。现在，让我们看看它在纯旋转情况下告诉我们什么。当 $\boldsymbol{F} = \boldsymbol{R}$ 时，我们得到 $\boldsymbol{C} = \boldsymbol{R}^{\mathsf{T}}\boldsymbol{R}$。由于 $\boldsymbol{R}$ 是一个旋转张量，它的转置是它的逆，所以 $\boldsymbol{R}^{\mathsf{T}}\boldsymbol{R} = \boldsymbol{I}$，即单位张量。

这意味着对于纯旋转，$\boldsymbol{C} = \boldsymbol{I}$，新的长度平方是 $d\boldsymbol{X}^{\mathsf{T}}\boldsymbol{I} d\boldsymbol{X} = d\boldsymbol{X}^{\mathsf{T}} d\boldsymbol{X}$。长度根本没有改变！张量 $\boldsymbol{C}$ 成功地滤除了刚性旋转，告诉我们没有发生拉伸。这是一个优美的结果。

由此，我们可以定义​​Green-Lagrange 应变张量​​为 $\boldsymbol{E} = \frac{1}{2}(\boldsymbol{C} - \boldsymbol{I})$。这个张量直接度量长度平方的变化。对于我们的纯旋转，$\boldsymbol{E} = \frac{1}{2}(\boldsymbol{I} - \boldsymbol{I}) = \boldsymbol{0}$。我们已经找到了我们真实的应变度量——一个当且仅当局部运动是纯刚性旋转时才为零的量。

我们的物理描述必须对刚性运动不敏感，这一思想实际上是物理学的一个深刻原理，称为​​物质标架无关性原理​​，或​​客观性​​。这是关于物理现实本质的论断：自然法则不能依赖于观察者的位置或旋转运动。如果我在实验室长凳上研究橡胶的特性，而你在一个旋转的旋转木马上研究完全相同的这块橡胶，我们必须推导出相同的材料定律。你旋转的视角不能神奇地改变橡胶的刚度。

我们可以通过一个思想实验将其形式化。假设我们有一个运动 $\boldsymbol{x}(\boldsymbol{X},t)$。第二个观察者相对于第一个观察者在旋转和平移，他将看到物质点在不同的位置 $\boldsymbol{x}^{\star}(t) = \boldsymbol{c}(t) + \boldsymbol{Q}(t)\boldsymbol{x}(t)$，其中 $\boldsymbol{c}(t)$ 是平移，$\boldsymbol{Q}(t)$ 是一个随时间变化的旋转。那些具有独立于观察者的物理现实的量被称为​​客观的​​。

我们的运动学张量在这种变换下表现如何？

• 一个标量，要成为客观的，必须是不变的：$s^{\star} = s$。
• 一个在未变形体上定义的物质（或参考）张量，也必须是不变的：$\boldsymbol{S}^{\star} = \boldsymbol{S}$。
• 一个在当前的、运动的标架中定义的空间张量，必须简单地随观察者标架旋转：$\boldsymbol{T}^{\star} = \boldsymbol{Q}\boldsymbol{T}\boldsymbol{Q}^{\mathsf{T}}$。

让我们检查一下我们的关键量。 新的变形梯度变为 $\boldsymbol{F}^{\star} = \boldsymbol{Q}\boldsymbol{F}$。它不是不变的，所以 $\boldsymbol{F}$ 不是一个客观的物质张量。它内在地包含了关于观察者视角的信息。 但是右 Cauchy-Green 张量呢？

它是不变的！$\boldsymbol{C}$（以及由此衍生的 $\boldsymbol{E}$）是一个​​客观物质张量​​。这为我们使用它作为应变度量提供了最深刻的理由。它捕捉了材料状态的物理现实，与观察者无关。这就是为什么任何物理定律，比如材料的储存能函数 $\psi$，必须只通过像 $\boldsymbol{C}$ 这样的客观组合来依赖于变形梯度 $\boldsymbol{F}$，即 $\psi(\boldsymbol{F}) = \hat{\psi}(\boldsymbol{C})$。一块材料的能量不能仅仅因为我们决定围绕它旋转而改变！

到目前为止，我们只看了变形的一个快照。但是动态过程，即运动的影片呢？这里我们需要讨论“率”。自然的起点是​​速度梯度​​ $\boldsymbol{L} = \nabla\boldsymbol{v}$，它描述了空间中各点粒子速度的变化。这个张量可以清晰地分解为其对称和反对称部分：

• $\boldsymbol{D} = \frac{1}{2}(\boldsymbol{L} + \boldsymbol{L}^{\mathsf{T}})$，​​变形率张量​​，描述瞬时拉伸率。
• $\boldsymbol{W} = \frac{1}{2}(\boldsymbol{L} - \boldsymbol{L}^{\mathsf{T}})$，​​自旋张量​​，描述瞬时旋转率（涡量）。

现在我们必须问我们那个关键问题：这些率是客观的吗？我们应用我们叠加的刚体运动，发现速度梯度变换为 $\boldsymbol{L}^{\star} = \boldsymbol{Q}\boldsymbol{L}\boldsymbol{Q}^{\mathsf{T}} + \dot{\boldsymbol{Q}}\boldsymbol{Q}^{\mathsf{T}}$。那个额外的项 $\dot{\boldsymbol{Q}}\boldsymbol{Q}^{\mathsf{T}}$ 是观察者参考系的角速度！这意味着 $\boldsymbol{L}$ 不是客观的——它的值取决于观察者旋转的速度。自旋张量 $\boldsymbol{W}$ 也是如此。

但请看变形率 $\boldsymbol{D}$ 发生了什么。因为观察者的自旋项 $\dot{\boldsymbol{Q}}\boldsymbol{Q}^{\mathsf{T}}$ 是反对称的，当我们取对称部分得到 $\boldsymbol{D}^{\star}$ 时，这些来自 $\boldsymbol{L}^{\star}$ 及其转置的额外项完美地抵消了！结果是 $\boldsymbol{D}^{\star} = \boldsymbol{Q}\boldsymbol{D}\boldsymbol{Q}^{\mathsf{T}}$。变形率 $\boldsymbol{D}$ 的变换方式与一个客观空间张量应有的方式完全一致。它是一个真实的、客观的拉伸率度量。

这不仅仅是一个数学上的奇趣；它具有至关重要的实际意义。许多先进的材料模型，特别是用于流体或经历快速过程的金属的模型，被表述为应力变化率与变形率之间的关系：$\stackrel{\circ}{\boldsymbol{\sigma}} = \mathbb{C}:\boldsymbol{D}$。为了使这个物理定律是客观的，如果 $\boldsymbol{D}$ 是客观的，那么应力率 $\stackrel{\circ}{\boldsymbol{\sigma}}$ 也必须是客观的。

这里的重磅消息是：你在微积分中学到的普通时间导数 $\dot{\boldsymbol{\sigma}}$ 并不是客观的。在 constitutive law 中使用它是产生物理谬误的根源。

让我们用一个简单的、毁灭性的例子来看看为什么会这样。想象一根已经承受一定拉力 $\sigma_{0}$ 的杆，我们只是以恒定的角速度旋转它。没有新的拉伸发生，所以变形率 $\boldsymbol{D}$ 处处为零。一个不正确的、非客观的定律，如 $\dot{\boldsymbol{\sigma}} = \mathbb{C}:\boldsymbol{D}$，会预测 $\dot{\boldsymbol{\sigma}} = \boldsymbol{0}$。这意味着在固定的实验室参考系中看到的应力张量保持不变。但这在物理上是荒谬的！杆中的拉力是一个与材料绑定的物理实体；拉力的方向必须随杆一起旋转。错误的定律预测应力保持指向，比如说，水平方向，而杆在其下方旋转。这就产生了一种没有物理基础的“伪应力”。这个误差的量级不小；在旋转角度为 $\theta_f$ 后，应力误差的 Frobenius 范数为 $\sqrt{2}\sigma_{0}|\sin(\theta_{f})|$，这甚至可能比实际应力还要大！

为了解决这个问题，我们必须使用一个​​客观应力率​​。这些是巧妙构造的导数（如​​Jaumann 率​​或​​Green-Naghdi 率​​），它们实质上是在一个与材料共旋的参考系中计算时间导数，从而移除了由纯自旋引起的非客观部分。在我们那个 $\boldsymbol{D}=\boldsymbol{0}$ 的旋转杆例子中，任何有效的客观应力率都会正确地得出零，这表明没有由变形产生新的应力。

在工程软件（如有限元法）中不使用客观率将是灾难性的。当一个部件经受大旋转时，即使实际应变很小，程序也会预测应力从无到有地出现。想象一下，用一个连旋转的基本物理都搞错的程序来设计一个旋转的涡轮叶片或汽车悬架部件。

因此，我们从一个如何描述运动的简单问题开始的旅程，一直被一个深刻的对称性原理——客观性所引导。这个原理迫使我们发现了真实的、与观察者无关的应变（$\boldsymbol{C}$、$\boldsymbol{E}$）和应变率（$\boldsymbol{D}$）的度量，并揭示了以尊重物理定律基本统一性的方式描述变化是何等微妙而又至关重要。

现在我们已经熟悉了刚体运动的语言，您可能会倾向于认为这是一个相当专门的主题，是物理学中一个迷人但有些古老的角落，处理的是旋转的陀螺和行星轨道。事实远非如此。实际上，我们刚刚发展的这些概念，并不仅仅是力学在一种特定类型物体上的应用；它们是我们理解几乎所有力学系统的基石。刚体运动的理想状态充当了一个基本的基准线，一个“无变化”的完美状态，与之相比，更为复杂的拉伸、弯曲和流动的现实才得以衡量。它是一个无声的、不变的背景，赋予了变形那嘈杂而充满活力的交响乐以意义。在本章中，我们将看到这一理想如何在从最抽象的材料科学理论到最实际的工程挑战，再到分子的复杂舞蹈等领域中，扮演着惊人多样且至关重要的角色。

刚体运动学最深刻的作用之一，是作为我们更普适的可变形材料理论的一致性检验。任何声称能描述一个柔性、可变形物体的理论，都必须在其特殊情况下，正确地再现刚体的行为。如果它通不过这个简单的测试，那么这个理论就是有缺陷的。

考虑任何一块材料的角动量积分平衡。一个深刻而优美的结果，我们称之为柯西第二运动定律，是这个平衡定律意味着柯西应力张量的对称性，$\boldsymbol{\sigma} = \boldsymbol{\sigma}^{\mathsf{T}}$。令人惊讶的是，这个对称性的推导不需要对运动本身做任何假设。它是在微观层面不存在“内力偶”的结果。这意味着无论物体是在变形、平移，还是仅仅作为一个整体刚性旋转，应力张量都必须是对称的。刚体情况远非一个例外，它是普适定律必须满足的一个关键情景。

这个思想被形式化为一个宏大的原理，即​​物质标架无关性​​，或称客观性。它是现代力学的一大支柱，指出材料的本构律——即关联应力与应变的规则——必须独立于观察者。观察者意味着什么？观察者只是从他们自己的参考系中观察运动的人，而这个参考系本身可能在平移和旋转。改变观察者在数学上等同于在现有变形上叠加一个刚体运动。

所以，客观性原理要求我们编写物理定律的方式，必须使其预测不会因为我们这些观察者决定旋转而改变。这带来了巨大的影响。当我们为像肌腱这样的软生物组织建模时，我们需要测量它被拉伸了多少。一个幼稚的变形度量可能会被肌腱在空间中的摆动所迷惑。标架无关性迫使我们寻找对这类刚体运动“视而不见”的数学量。其中一个量就是右 Cauchy-Green 张量，$\boldsymbol{C} = \boldsymbol{F}^{\mathsf{T}}\boldsymbol{F}$，其中 $\boldsymbol{F}$ 是变形梯度。在叠加刚性旋转下，$\boldsymbol{F}$ 会改变，但 $\boldsymbol{C}$ 却奇迹般地保持不变。一个只依赖于 $\boldsymbol{C}$ 的应变能函数因此保证是客观的；它测量的是真实的变形，而不是物体作为一个整体的无关紧要的旋转。通过这种方式，刚体运动的概念就像一个强大的过滤器，让我们能够从数学上的无稽之谈中筛选出有物理意义的理论。

客观性的抽象原理在计算力学中找到了最具体的体现，这个领域使我们能够为从摩天大楼到航天器的所有东西构建虚拟原型。

为一个结构（比如一个钢框架）建模的第一步，是决定如何描述其节点的运动。在一个简单的二维模型中，我们从运动学中知道，平面上一个刚性截面的运动由三个数字描述：两个用于平移，一个用于旋转。在三维中，是六个数字：三个用于平移，三个用于旋转。这些不是随意的选择；它们是刚体的基本自由度，并且它们成为我们有限元模型中的基本变量。这是我们书写结构分析语言的字母表。

但刚体运动也像“机器中的幽灵”一样出现。想象一下模拟一个没有固定在地上的结构——比如一个漂浮在太空中的卫星。如果你推它一下，它既会变形又会进行刚体运动。我们的模拟软件必须足够聪明来识别这一点。一个无约束的结构可以自由平移或旋转而不会产生任何内部应变，因此，也不会有任何恢复力。在数学上，这表现为一个“奇异的”全局刚度矩阵 $\mathbf{K}$。一个对应于刚体运动的非零位移矢量 $\mathbf{u}_{rbm}$ 产生零应变，导致零应变能，$\frac{1}{2} \mathbf{u}_{rbm}^{\mathsf{T}} \mathbf{K} \mathbf{u}_{rbm} = 0$。这意味着 $\mathbf{K} \mathbf{u}_{rbm} = \mathbf{0}$，所以矩阵 $\mathbf{K}$ 有一个零空间。这不是一个错误；这是物理在正确地告诉我们物体可以自由移动。理解刚体运动的运动学对于诊断和恰当约束这些零能模态至关重要。

当我们分析结构稳定性时，这个问题变得更加关键。如果我们想知道一根细长柱子将在哪个临界载荷下屈曲，我们需求解一个广义特征值问题。然而，如果我们的虚拟柱子是无约束的，它也具有那些相同的刚体模态。屈曲方程会变得不适定，允许多个解，因为自由漂浮的柱子无法决定是屈曲还是仅仅漂走。在能够提出关于其稳定性的有意义问题之前，我们必须首先通过施加足够的边界条件来消除这些刚体运动。

当然，有时我们希望复杂装配体的某些部分表现得像刚体。在一个汽车发动机的模型中，将发动机缸体视为一个单一刚体，而像皮带和软管等更柔性的组件附着于其上，这可能是合理的。刚体运动的运动学为我们提供了精确的数学工具来强制实现这一点。我们可以编写约束方程，使我们网格中的一组节点“被奴役”，迫使它们作为一个单一的刚性部件一起运动，仅由少数几个主自由度控制。

也许最微妙和最美妙的联系出现在模拟动力学中。想象一下模拟一个正在高速旋转的柔性物体。客观性的基本原理告诉我们，纯粹的刚性旋转不应产生任何内部应变能。然而，许多简单的数值算法却违反了这一点！它们可能错误地将方向的改变解释为拉伸或压缩，导致能量的虚假产生。经过数千个时间步长，这种“数值热”可以完全破坏模拟的物理真实性。这一发现推动了复杂的“共旋”和“几何”积分格式的发展，这些格式经过精心构建，能够免疫于这种人为效应。要正确模拟汽车碰撞，其物理正确性取决于对刚性转动运动学的深刻尊重。

刚体运动学的影响范围远远超出了传统的工程结构，下探到微观甚至纳米世界。

考虑两个表面接触。这种相互作用的最初始阶段——整个接触力学和摩擦学领域的基础——是一个几何问题。我们可以通过想象一个刚性压头向一个表面移动来模拟接触的开始。它们之间的间隙是压头刚体运动和两个表面形状的函数。首次接触的时刻发生在某个点上这个间隙首次关闭为零时。通过分析这种接近的简单运动学，我们可以预测接触将在何处以及何时发生，这是理解摩擦、粘附和磨损的第一步。

让我们再放大一些，到分子的尺度。一个纳米粒子如何在流体中翻滚，或者一个蛋白质如何折叠？将这些复杂的实体建模为原子的刚性簇通常是很有用的。要在计算机中模拟这一点，我们必须回到我们研究过的运动方程。我们使用牛顿第二定律 $\dot{\mathbf{p}} = \mathbf{F}$ 来追踪质心。而对于方向，我们在一个体固坐标系中求解欧拉运动方程，这自动考虑了变化的惯性张量的复杂性。我们使用四元数来表示方向，这种方式既稳健又避免了奇异点。 如果我们想在一个恒定温度下模拟这个系统，我们必须将它耦合到一个“控温器”。例如，一个 Nosé–Hoover 控温器通过控制动能来工作。对于一个刚性纳米粒子，我们必须正确定义转动动能，$K_r = \frac{1}{2} \boldsymbol{\omega}^{\mathsf{T}} \mathbf{I} \boldsymbol{\omega}$，并认识到一个通用的三维物体恰好有三个转动自由度。支配行星运动的同样刚体动力学定律，也支配着单个分子的热运动。

于是，我们看到了一种宏大的统一。刚体运动这个抽象而完美的概念，不仅仅是一个课堂练习。它是我们衡量物理理论的黄金标准，是构建我们世界的实用工具，是一种同样适用于桥梁、骨骼和分子的普适运动语言。它是一个优雅简约的舞台，在其上，变形那无比复杂的戏剧得以展开。