应变张量

我们如何精确地描述一个物体的拉伸、挤压和扭转？虽然我们能直观地理解变形，但要用数学方式捕捉它，就需要一种能够区分真实形状变化与简单移动或旋转的工具。这正是应变张量理论所要解决的根本挑战，它是连续介质力学的基石，为量化材料如何变形提供了语言。没有这个工具，预测从橡皮筋到构造板块等一切物体的行为都将是不可能的。

本文深入探讨了应变张量的精妙世界，探索了其理论基础和广泛的实际应用。第一章“原理与机制”将引导您了解应变的概念发展，从简单的位移概念开始，逐步建立起能够处理大型复杂变形的稳健的有限应变张量。我们将揭示如何从数学上将纯应变与刚体运动分离开来。在此之后，“应用与跨学科联系”一章将展示这一个数学概念如何作为工程学、材料科学、地球物理学乃至固态物理学等领域的统一原理，使我们能够设计、预测和理解物理世界。

想象你正在看一块柔软的黏土。你可以挤压它、拉伸它、扭转它。你将如何用数学精度来描述黏土发生了什么？你不能只说“它移动了”，因为如果你只是将整块黏土从桌子的一边滑到另一边，它的形状根本没有改变。真正的变化，即变形，是关于黏土内部的点相对于彼此如何移动。

让我们说得更正式一些。我们可以将未变形的黏土看作点的集合。我们可以用每个黏土粒子的起始位置来标记它，我们称之为 $\mathbf{X}$。在我们挤压它之后，同一个粒子移动到了一个新位置 $\mathbf{x}$。我们可以写下的最简单的东西是​​位移矢量​​，$\mathbf{u} = \mathbf{x} - \mathbf{X}$。这个矢量告诉我们每个粒子移动了多远。

但正如我们从滑动黏土的例子中看到的，位移 $\mathbf{u}$ 本身并不是变形的度量。黏土块中的每个粒子可能都有相同的大位移，但如果它们都一起移动，就没有应变。应变不是关于绝对运动，而是关于相对运动。它是关于你在材料中的邻居相对于你如何移动。

为了捕捉这种相对运动，我们需要放大来看。让我们不仅考虑一个点 $\mathbf{X}$，还考虑一个非常近的邻居，位于 $\mathbf{X} + \mathrm{d}\mathbf{X}$。连接它们的微小矢量是 $\mathrm{d}\mathbf{X}$。变形后，这两个点移动到 $\mathbf{x}$ 和 $\mathbf{x} + \mathrm{d}\mathbf{x}$。连接它们的新微小矢量是 $\mathrm{d}\mathbf{x}$。局部变形的全部信息都编码在每个可能的微小矢量 $\mathrm{d}\mathbf{X}$ 如何转换为新矢量 $\mathrm{d}\mathbf{x}$ 的过程中。

事实证明，对于任何平滑的变形，这种转换在局部是线性的。这意味着有一个矩阵——或者更普遍地说，一个张量——来执行这种映射。我们称之为​​变形梯度张量​​，$\mathbf{F}$。它就像一个同时也会扭转和旋转的局部放大镜，精确地告诉我们一个点周围的邻域是如何变换的：

这个张量的分量就是最终坐标相对于初始坐标的所有可能的偏导数，$F_{ij} = \frac{\partial x_i}{\partial X_j}$。这个张量 $\mathbf{F}$ 功能非常强大；它包含了关于拉伸、剪切以及——重要的是——旋转的所有局部信息。

这里我们遇到了一个微妙而优美的问题。$\mathbf{F}$ 本身是我们正在寻找的应变度量吗？让我们用一个简单的思想实验来检验它。拿一把刚性尺子，只旋转它。它没有被拉伸、压缩或以任何方式变形。没有应变。

然而，尺子上的每个小矢量 $\mathrm{d}\mathbf{X}$ 都被旋转成了一个新的矢量 $\mathrm{d}\mathbf{x}$。这意味着变形梯度 $\mathbf{F}$ 是一个旋转矩阵，而不是单位矩阵。如果我们用 $\mathbf{F}$ 作为应变度量，我们会错误地得出结论，认为旋转后的尺子处于应变状态。

这给我们上了一堂深刻的课：一个真实的应变度量必须完全“无视”刚体运动。如果一个物体只进行平移或旋转，我们的应变度量必须恰好为零。因此，我们的任务是找到一种方法，从 $\mathbf{F}$ 中精准地移除旋转部分，只留下变形的纯粹本质。

我们如何将拉伸与旋转分开？这里有一个绝妙的物理洞见：旋转保持长度不变，而拉伸改变长度。所以，让我们不要看矢量本身，而是看它们的长度——或者更好的是，它们长度的平方，这样可以避免处理平方根。

我们原始微小矢量的长度平方是 $(\mathrm{d}S)^2 = \mathrm{d}\mathbf{X} \cdot \mathrm{d}\mathbf{X}$。

新矢量的长度平方是 $(\mathrm{d}s)^2 = \mathrm{d}\mathbf{x} \cdot \mathrm{d}\mathbf{x}$。

现在，我们使用我们的定义 $\mathrm{d}\mathbf{x} = \mathbf{F} \mathrm{d}\mathbf{X}$ 并将其代入第二个方程。使用矩阵乘法规则，这变成 $(\mathrm{d}s)^2 = (\mathbf{F} \mathrm{d}\mathbf{X})^T (\mathbf{F} \mathrm{d}\mathbf{X}) = \mathrm{d}\mathbf{X}^T \mathbf{F}^T \mathbf{F} \mathrm{d}\mathbf{X}$。

看看我们发现了什么！新的长度平方通过张量组合 $\mathbf{C} = \mathbf{F}^T \mathbf{F}$ 与旧矢量 $\mathrm{d}\mathbf{X}$ 相关。这就是​​右 Cauchy-Green 变形张量​​。如果我们只有一个旋转，$\mathbf{C}$ 会发生什么？在这种情况下，$\mathbf{F}$ 是一个旋转张量 $\mathbf{Q}$，其性质为 $\mathbf{Q}^T \mathbf{Q} = \mathbf{I}$（单位张量）。所以，对于纯旋转，$\mathbf{C} = \mathbf{I}$。通过对长度进行平方，我们巧妙地使旋转消失了！张量 $\mathbf{C}$ 是一个类似度量的对象，它只关心材料的拉伸和剪切，而不关心其在空间中的整体方向。

张量 $\mathbf{C}$ 量化了拉伸状态，但“应变”应该量化与原始未拉伸状态相比的变化。如果根本没有变形，那么 $\mathbf{F} = \mathbf{I}$，因此 $\mathbf{C} = \mathbf{F}^T \mathbf{F} = \mathbf{I}$。所以，“应变量”必须与 $\mathbf{C}$ 与单位张量 $\mathbf{I}$ 的差异有关。

这引出了​​Green-Lagrange 应变张量​​的定义，这是有限变形最基本的度量之一：

$\frac{1}{2}$ 这个因子是一个方便的约定，我们稍后会体会到它的好处。根据其构造，$\mathbf{E}$ 对于任何刚体运动都为零，这使其成为真实应变度量的理想候选者。它的分量直接告诉我们长度平方的变化。像 $E_{11}$ 这样的正对角分量意味着在第一个方向上的拉伸，而负值则意味着压缩。像 $E_{12}$ 这样的非对角分量量化了最初垂直的线之间角度的变化——这就是剪切应变。

我们也可以用位移矢量梯度 $\nabla_{\mathbf{X}}\mathbf{u}$ 直接表示 $\mathbf{E}$。经过一些代数运算后，我们得到了一个优美的结果：

这个表达式非常有启发性。它包含一个线性部分 $\frac{1}{2}(\nabla_{\mathbf{X}}\mathbf{u} + (\nabla_{\mathbf{X}}\mathbf{u})^{T})$ 和一个非线性的二次部分 $\frac{1}{2}(\nabla_{\mathbf{X}}\mathbf{u})^{T}(\nabla_{\mathbf{X}}\mathbf{u})$。这个非线性项是准确描述大变形的秘诀。

如果变形非常非常小，比如钢桥中的微小振动，会怎么样？在这种情况下，位移梯度是微小的数字，而 $\mathbf{E}$ 方程中的二次项——两个小数的乘积——变得可以忽略不计。如果我们忽略它，我们会得到一个简单得多的表达式：

这就是著名的​​无穷小应变张量​​，通常称为 Cauchy 应变。它就是位移梯度的对称部分。几个世纪以来，这都是弹性理论的基石，而且理由充分。它满足三个优雅的特性，使其成为小变形世界的完美工具：

1. 对于任何无穷小刚体运动，它都为零。
2. 它是“真实”有限 Green-Lagrange 应变的直接一阶近似。
3. 它与应力张量在能量上“共轭”，这是一个深刻的原则，确保我们对力和变形的描述在热力学上是一致的。

然而，一旦变形变大——想想拉伸一根橡皮筋或软生物组织的运动——无穷小理论就失效了。$\mathbf{E}$ 中的非线性项就不能再被忽略。例如，如果你对一个块体施加一个大的简单剪切，线性理论 ($\mathbf{\epsilon}$) 预测的是一个纯剪切状态。但完整的 Green-Lagrange 理论 ($\mathbf{E}$) 正确地预测了一个沿对角线之一的额外拉伸效应，这可以通过在厚橡皮筋上画一个正方形并对其进行剪切来观察到。这个非线性项不是一个数学上的复杂问题；它是对真实物理现象的描述。

还有一个更深层次的微妙和优美之处有待揭示。到目前为止，我们所有的描述都来自我们所谓的​​拉格朗日​​视角。我们一直“跟随着”材料粒子，通过参考它们在未变形体中的原始位置 $\mathbf{X}$ 来描述变形。Green-Lagrange 张量 $\mathbf{E}$ 是一个拉格朗日度量，因为它是在参考构形上定义的。

但如果我们想描述一条河流的流动呢？追踪每一个水分子的源头是荒谬的。更明智的做法是站在岸边一个固定的点 $\mathbf{x}$ 上，描述恰好此时此刻流过该点的水的速度和性质。这就是​​欧拉​​视角。

我们能从这个空间性的欧拉视角定义一个应变度量吗？是的，我们可以。我们不再问曾经在 $\mathbf{X}$ 处的材料如何变形，而是问：当前在 $\mathbf{x}$ 处的材料变形了多少？这涉及到将当前几何形状与参考几何形状进行比较，但用当前构形来表示一切。这条路径导向了​​Euler-Almansi 应变张量​​，$\mathbf{e}$。它被定义为 $\mathbf{e} = \frac{1}{2}(\mathbf{I} - \mathbf{b}^{-1})$，其中 $\mathbf{b} = \mathbf{F}\mathbf{F}^T$ 是左 Cauchy-Green 张量，是 $\mathbf{C}$ 的空间对应物。

对于同一个物理变形，拉格朗日应变和欧拉应变的数值是不同的。考虑一个简单的单轴拉伸，拉伸因子为 $\lambda$。拉格朗日轴向应变为 $E_{11} = \frac{1}{2}(\lambda^2 - 1)$，而欧拉轴向应变为 $e_{11} = \frac{1}{2}(1 - \lambda^{-2})$。它们不一样！但它们并不矛盾。它们是用不同的“尺子”来描述同一个现实。拉格朗日应变度量的是相对于初始长度的长度变化，而欧拉应变度量的是相对于最终长度的长度变化。这就像描述一个人的身高猛增：你可以说他们相对于初始身高长高了20%（拉格朗日），或者说他们的初始身高比最终身高矮了16.7%（欧拉）。两种说法都是正确的；它们只是使用了不同的参考点。

这一系列张量——$\mathbf{F}$、$\mathbf{C}$、$\mathbf{E}$、$\mathbf{b}$、$\mathbf{e}$、$\mathbf{\epsilon}$——并不是一堆随意定义的、令人困惑的集合。它是一个联系紧密的工具家族，是一组数学对象的交响曲，使我们能够以极高的精度和从不同的角度描述变形力学。它们都通过前推和后拉操作相互关联，使我们能够随心所欲地在拉格朗日世界和欧拉世界之间进行转换。

这些张量优雅的数学背后有着深刻的物理意义。材料在变形过程中体积不变（一种​​等容​​变形）的条件，可以表述为一个涉及 Green-Lagrange 应变张量不变量的简单代数方程。

在最前沿的应用中，从设计新材料到创建逼真的虚拟手术模拟器，选择正确的应变度量至关重要。一个变形物体中储存的应变能——比如说，一块被手术工具戳动的虚拟组织——必须使用像 $\mathbf{C}$ 或 $\mathbf{E}$ 这样的客观度量来计算。这确保了外科医生通过触觉反馈设备感受到的力是真实的弹性力，而不是仅仅因为在空间中移动或旋转组织而产生的虚假产物。

因此，从“物体是如何拉伸的？”这个简单、直观的问题出发，我们进入了一个丰富的张量世界。这些数学结构优雅地将拉伸与旋转分离开来，提供了连贯的拉格朗日和欧拉视角，并构成了现代连续介质力学的基础，使我们能够理解和预测从行星地壳到我们身体细胞的一切行为。

在掌握了应变张量的数学工具后，我们可能感觉有点像一个花了几周时间学习磨凿子和调锯子的学徒。这是一项基本技能，但真正的乐趣在于当我们最终走进作坊，开始创造出奇妙的东西时。现在就是那个时候。让我们带着我们新获得的工具，去科学世界巡游一番，看看它们能做什么。我们会惊讶地发现，这个单一而优雅的概念——一个精确描述变形的张量——是一把万能钥匙，解开了在初看起来彼此毫无关联的领域中的秘密。

让我们从工程学的有形世界开始，在这里，材料被弯曲、拉伸和扭转，以创造我们日常生活中的结构和设备。我们如何确保一座桥不会屈曲，一架飞机的机翼不会折断，或者一根简单的橡皮筋会恢复到原来的形状？答案始于量化应变。

考虑一个最有趣和最简单的例子：给一个玩具气球充气。当你向气球吹气时，它的半径 $R$ 从初始尺寸 $R_0$ 增长。橡胶本身发生了什么？对于表面上的任何一小块区域，材料在表面内所有方向上都受到均等的拉伸。我们的应变张量完美地捕捉了这种各向同性。代表拉伸的对角分量是相同的，并且取决于比率 $(R/R_0)^2$，而非对角剪切分量为零。该张量用精确的术语告诉我们，我们的直觉已经怀疑的事情：变形是一种纯粹、均匀的二维扩张。

现在，让我们尝试一些更复杂的东西。想象一下，拿一块薄而平的矩形板，把它弯成一个完美的圆弧，就像把它包裹在一个巨大的隐形圆柱体上一样。这就是弯曲梁中发生的事情的本质。应变张量在这里揭示了什么？它讲述了一个美丽的故事。应变不是均匀的！曲线外侧的材料 ($X_2 > 0$) 被拉伸 ($E_{11} > 0$)，而曲线内侧的材料 ($X_2 0$) 被压缩 ($E_{11} 0$)。而在正中间，沿着一条“中性线”($X_2 = 0$)，根本没有拉伸。Green-Lagrange 张量的这个单一分量 $E_{11}$，随离中心线的距离而变化，包含了弯曲的全部秘密。正是这种应变变化产生了内部应力，使得梁能够抵抗弯曲。

现实世界的过程通常比单一、干净的动作更复杂。工厂里的一块金属可能会被拉伸，然后剪切，然后轧制。连续介质力学的强大之处在于我们可以将这些操作链接在一起。通过将每个连续步骤的变形梯度张量相乘，我们可以得到总变形。例如，如果我们取一个物体，将其均匀拉伸，然后施加一个简单的剪切，得到的 Green-Lagrange 应变张量并不仅仅是拉伸应变和剪切应变之和。它包含了揭示大变形耦合、非线性性质的交叉项。最终状态取决于其整个变形历史。

钢和混凝土的经典工程学常常依赖于应变非常小的简化假设。但在现代世界，我们越来越多地设计和使用那些注定要经历巨大变形的材料——想想软体机器人抓手里的柔软、灵活的致动器，或生物医学设备中的超弹性聚合物。对于这些材料，小应变理论是无用的，而像 Green-Lagrange 张量这样的有限应变张量的全部威力是不可或缺的。

想象一块柔软的、类似橡胶的材料块受到一种奇特的非线性剪切，其中水平位移取决于垂直位置的平方。无穷小应变分析会给出一个非常误导的图像。然而，Green-Lagrange 张量揭示了真实情况。它不仅正确预测了预期的剪切应变，还预测了一个非零的*正应变*——一种纯粹由大变形的几何形状引起的拉伸或压缩。软体机器人工程师依靠这种精确的理解来设计能够以复杂方式弯曲、扭转和抓握的致动器。一个更引人注目的例子是一个厚球壳的完全外翻——即内外翻转，这是某些类型的软致动器的模型。在这里，用球坐标计算的应变张量使我们能够分析材料在经历这种拓扑变形时发生的巨大拉伸和压缩。

除了仅仅描述变形，应变张量还是材料“本构”——其定义性行为法则——的基本输入。应变能密度函数 $W$ 告诉我们，在给定的应变量下，材料中储存了多少能量。对于一个简单的各向同性材料，$W$ 只取决于应变张量的不变量（即不随参考坐标系旋转而改变的量）。但对于先进的复合材料，比如用坚固的平行纤维增强的塑料呢？这种材料是各向异性的；它在纤维方向上要硬得多。我们如何捕捉这一点？我们建立一个应变能函数，它不仅包括应变张量 $\mathbf{E}$，还包括一个结构张量 $\mathbf{M} = \mathbf{a} \otimes \mathbf{a}$，该张量代表了优选的纤维方向 $\mathbf{a}$。现在的能量取决于像 $\mathbf{a}\cdot\mathbf{Ea}$ 这样的项，它度量了沿纤维方向的应变。这个框架使得材料科学家能够为从航空航天部件到生物医学植入物的应用定制设计具有特定方向性质的材料。

应变理论的影响范围远远超出了工程材料，为基础物理学的不同分支提供了一种通用语言。

在地球物理学中，构造板块的运动、地幔的流动以及岩石在压力下的压实都是变形过程。一个计算模拟可以用一个简单的变形梯度来描述岩石样本的变形，例如，沿一个轴拉伸它，同时沿另一个轴压缩它。由此，我们既可以计算 Green-Lagrange 张量 $\mathbf{E}$（它度量相对于岩石初始状态的应变），也可以计算 Euler-Almansi 张量 $\mathbf{e}$（它度量相对于最终变形状态的应变）。此外，我们可以计算雅可比行列式 $J = \det(\mathbf{F})$。这个简单的标量告诉我们体积的变化：如果 $J 1$，岩石正在被压实，其孔隙度被挤出；如果 $J > 1$，它正在膨胀。这个单一的数字对于理解从油藏动力学到地球地壳深处过程的现象至关重要。

也许最深刻和最美丽的联系是与固态物理学的微观世界。我们倾向于认为晶体是原子完美有序、刚性的排列。但当晶体受力时，这个晶格会变形。我们如何描述这一点？我们可以用其原胞基矢 $\mathbf{a}_i$ 来定义晶体结构。晶格的几何形状被编码在度规张量 $g_{ij} = \mathbf{a}_i \cdot \mathbf{a}_j$ 中，它包含了基矢之间的所有长度和角度。当晶体变形时，晶格矢量发生变化，度规张量也随之变化。事实证明，度规张量的变化 $\Delta g_{ij}$ 与拉格朗日应变张量成正比。这是一个惊人的发现！我们用来描述橡胶片整体拉伸的抽象应变张量，在晶体基本几何蓝图的扭曲中有一个直接的物理对应物。例如，一个剪切应变分量直接对应于晶体单胞角度的变化。

最后，应变张量提供了变形（运动学）与引起变形的力及能量（动力学和热力学）之间的关键联系。当我们使一个物体变形时，我们对它做功。这些能量去哪儿了？它可以作为弹性势能储存起来，也可以作为热量耗散掉。我们对单位初始体积的材料做功的速率 $P_0$ 有一个极其简单而强大的表达式：$P_0 = \mathbf{S}:\dot{\mathbf{E}}$。

这就是功共轭关系。它表明，应力功率是第二 Piola-Kirchhoff 应力张量 $\mathbf{S}$ 与其“共轭”运动学变量——Green-Lagrange 应变张量 $\mathbf{E}$ 的变化率——的双点积。这不仅仅是一个简洁的公式；它是一个深刻的物理一致性的陈述。它确保了我们对应力和应变的数学描述在能量上是相容的。它构成了本构建模的基础，使我们能够通过简单地对超弹性材料的应变能函数求关于应变张量的导数来推导出应力。这种关系是驱动我们能够创建从车祸到心脏瓣膜功能的预测性模拟的引擎。

从给气球充气这个简单的动作，到各向异性材料的复杂理论，再到晶体弹性的量子力学基础，应变张量是贯穿其中的共同线索。它证明了数学在多样化的物理世界结构中寻找统一性的力量，无论我们是在观察山脉、复合材料飞机机翼，还是单个金属晶体，它都揭示了相同的基本原理在起作用。我们确实创造了一些奇妙的东西。