左拉伸张量

当一种材料发生变形时，它会经历拉伸、压缩和旋转的复杂组合。虽然这些现象易于观察，但为了科学和工程目的精确描述它们却是一项重大挑战。变形梯度张量提供了完整的数学描述，但它内在地混合了纯拉伸和刚体旋转的效应。为了建立准确的物理模型——例如，理解金属何时会屈服或橡胶密封圈如何储存能量——我们必须解开这两种截然不同的现象。这正是左拉伸张量概念所要解决的核心问题。

本文为理解这一连续介质力学中的关键工具提供了全面的指南。我们将把变形过程分解为其基本组成部分，重点关注左拉伸张量的意义和功用。整个过程分为两个关键部分。

• 首先，在​​原理与机制​​部分，我们将探讨极分解的理论基础，区分左、右拉伸张量，并详细介绍它们的计算方法。
• 接着，在​​应用与跨学科联系​​部分，我们将看到这个抽象概念如何成为材料科学、计算工程等领域的强大工具，支撑着应变度量、本构法则和稳健模拟技术的建立。读完本文，左拉伸张量将不再仅仅是一个数学上的奇特概念，而是描述物质世界的一块基石。

想象一下，你拿一块橡皮泥，对它进行挤压、扭转和拉伸。最终的形状是一个复杂过程的结果。但如果你观察橡皮泥内部一个微小的尘埃点，它从起点到终点的运动可以归结为两个基本动作：​​拉伸​​和​​旋转​​。物理学及其数学描述的魔力就在于，它能将复杂的现实分解为简单而优美的思想。本章的使命正是为材料变形做到这一点。

关于一个点附近微小邻域如何变形的完整信息，被一个称为​​变形梯度​​的数学对象所捕捉，我们用符号 $\mathbf{F}$ 表示。它就像是变形的局部“配方”。如果你在原始未变形的物体（​​参考构型​​）中取一个无穷小的箭头 $d\mathbf{X}$，$\mathbf{F}$ 就会告诉你这个箭头在变形后的物体（​​当前构型​​）中变成了什么：$d\mathbf{x} = \mathbf{F} d\mathbf{X}$。

现在，你可能会想，“这很简单嘛。”但这个方程式看似简单的外表下隐藏着一个有趣的精妙之处。考虑一种被称为​​简单剪切​​的变形，即一块材料的变形如同其顶部相对于底部发生了侧向滑动。其变形梯度相当简单：

其中 $\gamma$ 表示剪切量。看到这个矩阵，你可能会倾向于认为没有涉及旋转——只是滑动而已。你可能还会说，由于对角线元素为1，所以在水平或垂直方向上没有拉伸。然而，我们即将看到，我们的直觉有时会出人意料地错误。连续介质力学为我们提供了一个强大的工具——​​极分解​​，来剖析 $\mathbf{F}$，揭示其中隐藏的真实拉伸和旋转。我们关于​​左拉伸张量​​的故事就从这里开始。

​​极分解定理​​是力学的一块基石。它告诉我们，任何变形 $\mathbf{F}$（只要其在物理上可能，即 $\det(\mathbf{F}) > 0$）都可以唯一地分解为一次纯拉伸后跟一次纯旋转，或者一次纯旋转后跟一次纯拉伸。这为我们讲述变形提供了两种“故事”。

​​故事1：先拉伸，后旋转​​ 在这个版本中，我们写成 $\mathbf{F} = \mathbf{R}\mathbf{U}$。运算从右向左进行：

1. 一个材料元首先被​​右拉伸张量​​ $\mathbf{U}$ 拉伸。
2. 然后，被拉伸的材料元通过​​旋转张量​​ $\mathbf{R}$ 进行刚性旋转，达到其最终的朝向。

​​故事2：先旋转，后拉伸​​ 或者，我们也可以写成 $\mathbf{F} = \mathbf{V}\mathbf{R}$。其顺序为：

1. 一个材料元首先通过同一个旋转张量 $\mathbf{R}$ 进行刚性旋转。
2. 然后，被旋转的材料元通过​​左拉伸张量​​ $\mathbf{V}$ 拉伸成其最终形状。

你可以看到这里的美妙之处：最终的变形 $\mathbf{F}$ 是相同的，两个故事中的刚性旋转 $\mathbf{R}$ 也是相同的。唯一的区别在于两个拉伸张量 $\mathbf{U}$ 和 $\mathbf{V}$，以及它们作用的顺序。这自然引出了一个问题：它们之间到底有什么区别？

右拉伸张量 $\mathbf{U}$ 和左拉伸张量 $\mathbf{V}$ 之间的区别在于视角的不同。这取决于你从哪个“世界”进行观察。

​​右拉伸张量 $\mathbf{U}$ 存在于参考构型中​​。它从物体的未变形状态出发，描述拉伸过程。它有一组特殊的三个正交方向（其特征向量），这些方向上的材料纤维只被纯粹拉伸，而没有任何剪切。在变形前，它们沿着这些特殊方向；经过 $\mathbf{U}$ 的拉伸后，它们变长或变短，但仍然指向相同的特殊方向。这些是在参考构型中的​​主拉伸方向​​。

​​左拉伸张量 $\mathbf{V}$ 存在于当前构型中​​。它描述了在变形后的物体中看到的拉伸结果。它也有一组特殊的正交主方向。在最终变形的物体中，沿着这些方向之一的线段之所以特殊，是因为它来自原始物体中一个只被拉伸而未被剪切的线段。这些是在当前构型中的​​主拉伸方向​​。

那么它们之间有什么联系呢？这种联系惊人地简单。左拉伸张量 $\mathbf{V}$ 的主方向，正是右拉伸张量 $\mathbf{U}$ 主方向经过旋转后的版本。如果 $\mathbf{n}$ 是 $\mathbf{U}$ 的一个主方向，那么 $\mathbf{V}$ 对应的那个主方向就是 $\mathbf{R}\mathbf{n}$。想象一下，在未变形的橡皮泥上画了一组三个特殊的正交坐标轴；经过完全变形后，这些轴被拉伸并旋转。最终形状中的这组新轴就是 $\mathbf{V}$ 的主方向。

那么沿这些方向的拉伸量又如何呢？​​主拉伸​​（$\mathbf{U}$ 和 $\mathbf{V}$ 的特征值）对于这两个张量是完全相同的。这在物理上完全合理：一根纤维被拉伸的量不应该取决于你是在刚性旋转之前还是之后来描述这个过程。将这一切联系起来的关系是一个简单的张量旋转：$\mathbf{V} = \mathbf{R}\mathbf{U}\mathbf{R}^{\mathsf{T}}$。这意味着 $\mathbf{V}$ 与 $\mathbf{U}$ 是相同的拉伸，只是在旋转后的坐标系中观察而已。一个优美的例子 表明，对于一个由沿轴向拉伸后跟一个角度为 $\theta$ 的旋转组成的简单二维变形，$\mathbf{U}$ 的主方向是原始坐标轴，而 $\mathbf{V}$ 的主方向是这些坐标轴旋转 $\theta$ 角后的方向。

这一切听起来非常抽象，但我们如何从一个给定的变形 $\mathbf{F}$ 中实际找到 $\mathbf{U}$ 和 $\mathbf{V}$ 呢？我们不必靠猜测。我们用一个巧妙的技巧来消除旋转 $\mathbf{R}$。想想看，你可以通过平方来求一个数的绝对值，从而消除其符号。我们在这里做类似的事情。

我们定义两个辅助张量：

• ​​右柯西-格林张量​​：$\mathbf{C} = \mathbf{F}^{\mathsf{T}}\mathbf{F}$。如果我们代入 $\mathbf{F}=\mathbf{R}\mathbf{U}$，我们得到 $\mathbf{C} = (\mathbf{R}\mathbf{U})^{\mathsf{T}}(\mathbf{R}\mathbf{U})=\mathbf{U}^{\mathsf{T}}\mathbf{R}^{\mathsf{T}}\mathbf{R}\mathbf{U}$。因为 $\mathbf{R}$ 是一个旋转，所以 $\mathbf{R}^{\mathsf{T}}\mathbf{R}=\mathbf{I}$（单位矩阵），又因为 $\mathbf{U}$ 是对称的，所以 $\mathbf{U}^{\mathsf{T}}=\mathbf{U}$。结果就是简单的 $\mathbf{C} = \mathbf{U}^2$。
• ​​左柯西-格林张量​​：$\mathbf{B} = \mathbf{F}\mathbf{F}^{\mathsf{T}}$。如果我们代入 $\mathbf{F}=\mathbf{V}\mathbf{R}$，我们得到 $\mathbf{B} = (\mathbf{V}\mathbf{R})(\mathbf{V}\mathbf{R})^{\mathsf{T}}=\mathbf{V}\mathbf{R}\mathbf{R}^{\mathsf{T}}\mathbf{V}^{\mathsf{T}}$。这可以简化为 $\mathbf{B} = \mathbf{V}^2$。

这就是关键！要找到 $\mathbf{U}$，我们只需计算 $\mathbf{C} = \mathbf{F}^{\mathsf{T}}\mathbf{F}$ 并求其唯一的​​对称正定​​ (SPD) 平方根。要找到 $\mathbf{V}$，我们计算 $\mathbf{B} = \mathbf{F}\mathbf{F}^{\mathsf{T}}$ 并求其 SPD 平方根。SPD 属性至关重要；它保证了这些张量代表了真实的物理拉伸（正特征值），而不是某种数学上的奇特产物。

让我们通过一个具体案例来看看它的实际应用。考虑一个由 $\mathbf{F} = \begin{pmatrix} \sqrt{3} & -3/4 & 0 \\ 1 & 3\sqrt{3}/4 & 0 \\ 0 & 0 & 1 \end{pmatrix}$ 给出的变形。计算 $\mathbf{C} = \mathbf{F}^{\mathsf{T}}\mathbf{F}$ 会得到一个惊人地简单的对角矩阵：

求平方根现在变得非常简单——我们只需对对角线元素取平方根：

这告诉我们，材料在 $X_1$ 方向上被拉伸了2倍，在 $X_2$ 方向上被拉伸了1.5倍，而在 $X_3$ 方向上完全没有拉伸。有了 $\mathbf{U}$，我们就可以找到旋转 $\mathbf{R} = \mathbf{F}\mathbf{U}^{-1}$，结果是绕 $X_3$ 轴旋转了 $\pi/6$ 弧度（30度）。然后，我们可以找到左拉伸张量 $\mathbf{V} = \mathbf{RUR}^{\mathsf{T}}$。抽象的概念通过计算变得生动起来！

对于那些喜欢在数学中看到统一性的人来说，这里还有更深的一层。这种物理上的极分解与线性代数中的一个基本工具——​​奇异值分解 (SVD)​​——有着内在的联系。SVD 表明，任何矩阵 $\mathbf{F}$ 都可以写成 $\mathbf{F} = \mathbf{W}\mathbf{\Sigma}\mathbf{V}_{\text{svd}}^{\mathsf{T}}$，其中 $\mathbf{W}$ 和 $\mathbf{V}_{\text{svd}}$ 是旋转矩阵，而 $\mathbf{\Sigma}$ 是一个由正数组成的对角矩阵，这些正数被称为奇异值。事实证明，我们极分解的各个组成部分，只是 SVD 组元的优雅组合：奇异值就是主拉伸，而张量 $\mathbf{U}$、$\mathbf{V}$ 和 $\mathbf{R}$ 可以直接由 $\mathbf{W}$、$\mathbf{\Sigma}$ 和 $\mathbf{V}_{\text{svd}}$ 构造出来。这是一个美丽的印证，表明物理直觉和数学结构是同一枚硬币的两面。

那么，我们为什么要费这么大劲呢？这仅仅是数学体操吗？完全不是。这种分解对于建立物理定律至关重要。基本的​​材料坐标系无关性原则​​（或称客观性原则）指出，材料中储存的能量应仅取决于它被拉伸的程度，而不是它作为整体如何被刚性旋转。你的汽车轮胎不会仅仅因为汽车在转弯而储存更多的能量。

这一原则带来了一个深远的后果。如果我们将材料的储存能表示为完整变形的函数，即 $\Psi(\mathbf{F})$，那么要确保遵循这一原则会非常困难。然而，如果我们将能量写成右拉伸张量的函数 $\Psi(\mathbf{U})$，或右柯西-格林张量的函数 $\Psi(\mathbf{C})$，客观性就自动满足了！这是因为，如果我们对系统施加一个刚性旋转 $\mathbf{Q}$，$\mathbf{F}$ 会变为 $\mathbf{Q}\mathbf{F}$，但 $\mathbf{U}$ 和 $\mathbf{C}$ 却完全保持不变。它们是衡量纯变形的“客观”度量。

有趣的是，左拉伸张量 $\mathbf{V}$ 不具备此特性；它会随系统一起旋转。一个能量函数 $\Psi(\mathbf{V})$ 只有在材料是​​各向同性​​（在所有方向上表现相同）时才是客观的。这一微妙的区别引导着物理学家和工程师为从钢梁到橡胶轮胎再到活体组织等各种事物建立精确的模型。通过将变形分解为其最基本的部分——拉伸和旋转——我们解锁了对支配物质世界法则的更深层次的理解。左拉伸张量 $\mathbf{V}$ 不仅仅是 $\mathbf{U}$ 的一个替代品；它在这个宏大的故事中扮演着关键角色，从此时此地的视角提供了对拉伸状态的观察。

现在我们已经掌握了连续介质变形的原理和机制，你可能会好奇，“所有这些复杂的理论是做什么用的？”我们已经拆解了变形梯度张量 $\mathbf{F}$，利用极分解将其分离为纯拉伸和纯旋转。特别地，我们认识了左拉伸张量 $\mathbf{V}$，它从最终变形构型的视角描述了拉伸状态。

这仅仅是一场数学游戏吗？一种线性代数的巧妙技巧？远非如此。拉伸与旋转的分离是力学中最深刻、最有用的思想之一，其影响深远，触及材料科学、工程学、计算机模拟乃至纯粹几何学。它使我们能够提出并回答一些非常合理的问题。如果我们弯曲一根金属棒，其中有多少是简单的旋转，又有多少是可能导致其失效的材料实际拉伸？左拉伸张量 $\mathbf{V}$ 正是此时此地这些信息的保管者。让我们来探索它所编织的美丽画卷。

第一个也是最基本的应用是为“应变”或“某物被拉伸了多少”这一概念赋予精确的含义。想象一张被拉伸的橡胶薄片。不同的部分以不同的量和不同的方向伸展。左拉伸张量 $\mathbf{V}$ 捕捉了各处的局部拉伸状态。它的特征值，即主拉伸 $\lambda_1, \lambda_2, \lambda_3$，告诉我们沿三个相互垂直方向的拉伸因子——这些方向上，无穷小的球体变形为了椭球体。

但物理学家或工程师想要一个在没有变形时为零的数字。这正是应变张量的作用。定义应变的方式不止一种；它取决于你的视角。

如果你是一个在最终变形构型中的观察者，看着扭曲的物体，你会很自然地定义一个与你当前所见几何相关的应变。这就是欧拉-阿尔曼西应变张量 $\mathbf{e}$ 的精神所在。这里就有了第一个美妙的联系：这个空间应变度量与左拉伸张量 $\mathbf{V}$ 通过一个直接而优雅的关系联系在一起：

你可以看到其中的逻辑。如果没有拉伸，$\mathbf{V}=\mathbf{I}$，应变 $\mathbf{e}$ 为零。如果材料被拉伸，$\mathbf{V}$ 的主值为 $\lambda_i > 1$，使得应变的主值 $e_i = \frac{1}{2}(1 - \lambda_i^{-2})$ 为正。如果被压缩，$\lambda_i 1$，应变为负。它的表现完全符合我们的直觉。这表明 $\mathbf{V}$ 不仅仅是矩阵分解中的一个抽象因子；它是我们所见世界中测量变形的物理基础。

值得注意的是，如果我们采取不同的视角，即从初始未变形状态向前看（拉格朗日视角），我们会定义一个不同的应变度量，即格林-拉格朗日应变 $\mathbf{E}$。你可能已经猜到，这个度量并不直接与 $\mathbf{V}$ 相关，而是与它的对应物——右拉伸张量 $\mathbf{U}$ 相关。自然界提供了一种美妙的对称性：两种视角，两种拉伸张量，每一种都有其自然的应变度量。

知道如何测量拉伸只是故事的一半。另一半，对材料科学家来说也许更激动人心的一半，是理解材料如何响应拉伸。一根钢梁的响应与一根橡皮筋或一块橡皮泥截然不同。材料的这种“特性”被我们称为本构模型，或应力-应变关系。

考虑一种超弹性材料，如橡胶。当你拉伸它时，它储存能量；当你松开它时，它释放能量，迅速恢复原状。对于一个简单的各向同性材料（其性质在所有方向上都相同），储存的能量应仅取决于拉伸的量，而不应取决于它所经历的任何刚体旋转。毕竟，一个橡胶球不关心你是否旋转它；它的内能不会改变。

这正是极分解力量的闪光之处。储存能 $W$ 不能依赖于完整的变形 $\mathbf{F}$，因为 $\mathbf{F}$ 包含旋转。它必须只依赖于拉伸。主拉伸 $\lambda_i$（$\mathbf{V}$ 的特征值）是完美的候选者。对于各向同性材料，应变能密度是这些拉伸的对称函数：

任何复杂的应变不变量张量函数都可归结为这个简单、直观的思想。材料的整个复杂响应都编码在它如何根据三个主拉伸值储存能量的方式中。

由此，应力——材料为抵抗变形而施加的内力——可以直接求得。Kirchhoff 应力的主值 $\tau_i$（一个物理上重要的应力度量）与能量之间通过一个异常简单的表达式相关联：

这为工程师提供了一个直接的“配方”：如果你能写出一种材料的能量函数（如用于橡胶的 Ogden 模型），你就可以立即计算出任何给定变形下的应力。这是计算工程学的核心，让我们能够模拟和设计从垫圈、密封件到轮胎和生物医学设备的一切，而所有这些都建立在由 $\mathbf{V}$ 捕捉的拉伸基本概念之上。

但是那些无法弹回的材料呢？当你弯曲一个回形针时，它会保持弯曲。这就是塑性的世界，是永久、不可恢复变形的领域。我们简洁的弹性理论似乎在这里失效了。但奇妙的是，它并没有。这些概念只是提升到了一个更高的抽象层次。

现代有限塑性理论使用了一个绝妙的思想：变形的乘法分解。它认为总变形 $\mathbf{F}$ 可以被看作一个两步过程：首先是永久的塑性变形 $\mathbf{F}_p$，它重新排列了材料的内部结构；然后是从这个新状态出发的可恢复的弹性变形 $\mathbf{F}_e$。

其神奇之处在于，我们所有关于弹性的推理现在都适用于弹性部分 $\mathbf{F}_e$。我们可以对 $\mathbf{F}_e$ 进行极分解，找到一个弹性左拉伸张量 $\mathbf{V}_e$。正是这个张量决定了材料的应力响应。在某种意义上，材料的记忆很短；它当前的应力状态只取决于它从最近的塑性变形构型被弹性拉伸了多少。

从 $\mathbf{V}_e$ 出发，我们可以定义更复杂的量，如对数弹性应变 $\mathbf{h}_e = \ln(\mathbf{V}_e)$。这个度量有一个方便的特性：对于小的弹性变形（即使总塑性变形很大），这些应变表现出可加性，这对计算机模拟是一个巨大的简化。弹性左拉伸张量 $\mathbf{V}_e$ 充当了一个门户，使我们能够将弹性的清晰逻辑应用于永久变形的更复杂世界，这是金属成形、岩土工程和碰撞分析的基石。

我们最后的联系将我们带入动态的运动世界。想象一下，试图描述在一锅浓稠的蜂蜜中搅拌的勺子的状态。勺子在旋转，同时蜂蜜也在变形。如果你想谈论蜂蜜中应力累积的速率，你就会遇到一个问题：你看到的改变有多少是由于蜂蜜被拉伸造成的，又有多少仅仅是因为你正在观察的那块蜂蜜在旋转？

物理学必须是客观的——与观察者的旋转参考系无关。这意味着应力张量的简单时间导数 $\dot{\boldsymbol{\sigma}}$ 不是一个具有物理意义的量。我们需要创建“客观率”，以智能地减去刚体旋转的影响。左拉伸张量 $\mathbf{V}$ 及其时间演化是定义这些率的关键。

这个故事中的两个主要“角色”是变形率张量 $\mathbf{D}$（拉伸的速率）和涡量张量 $\mathbf{W}$（旋转的速率）。对于没有旋转的特殊运动，如纯拉伸，涡量为零。在这种简单情况下，许多不同的客观率定义是一致的。但对于一般运动，用于校正的不同涡量张量的选择会导致不同的客观率，如 Zaremba-Jaumann 率或对数率。例如，对数率使用的涡量张量与拉伸张量 $\mathbf{V}$ 主轴的旋转密切相关。

这个看似深奥的主题是计算固体和流体力学的核心。有限元法 (FEM) 软件中使用的公式，即所谓的共旋公式，正是建立在这个思想之上。对于模拟结构中的每一个小块（单元），程序会从极分解中计算其整体旋转 $\mathbf{R}$。然后，它会虚拟地“反旋转”该单元，在这个简单、未旋转的、拉伸较小的框架中计算应力和应变，然后将结果旋转回全局坐标系。这使得工程师能够准确模拟经历复杂大旋转的结构，如飞机的扑动翅膀或柔性桥梁的屈曲，而不会在变形与涡旋的令人眼花缭乱的舞蹈中迷失方向。

从其在矩阵分解中的卑微起源 开始，左拉伸张量 $\mathbf{V}$ 展现为一个具有深刻统一性的概念。它是测量应变的关键，是定义材料特性的关键，是分离弹性和塑性的关键，也是在不断运动的世界中保持客观视角的关键。这是一个美丽的例子，说明一个优雅的数学思想如何提供一个晶莹剔透的镜头，来观察和解释物理世界。