形变率张量

玻尔百科

关键要点

形变率张量通过分解速度梯度得到，它清晰地将纯形状改变（应变率）与局部刚体旋转（自旋）分离开来。
它构成了连续介质力学中本构律的基础，将材料的内应力与其形变率联系起来。
该张量是计算由黏性引起的能量耗散的基础，并在流体动力学、固体力学甚至生物学中有着广泛的应用。

引言

要理解任何连续材料的运动——从管道中流动的水到正在锻造的钢铁——我们都需要超越简单的速度概念。材料如何拉伸、压缩、剪切和旋转的真实情况，隐藏在速度从一点到其邻近点的变化方式之中。这种复杂的局部运动带来了一个挑战：我们如何能精确量化纯粹的形状变化，并将其与简单的平移或旋转分离开来？答案在于连续介质力学核心的一个强大数学工具：形变率张量。本文将对这一基本概念进行全面概述。第一部分“原理与机制”将深入探讨运动的数学分解，揭示形变率张量如何将应变从自旋中分离出来，并与不可压缩性和能量耗散等物理性质联系起来。接下来的“应用与跨学科联系”部分将展示该张量在定义从简单流体到复杂固体的材料行为中的重要作用，以及其在材料科学乃至生物学等领域出人意料的相关性。

原理与机制

想象一下，你正在观察一条河流。你看到树叶和嫩枝被水流带着走。有些在小涡流中旋转，有些在进入更快的溪流时被拉长，而另一些则只是平静地漂浮着。如果我们放大观察一个微小的、想象中的水立方体，它会经历什么样的运动？它可能被带到下游（平移），可能像陀螺一样旋转（旋转），其形状也可能被挤压、拉伸或扭曲（形变）。为了真正理解一个形变中材料的物理学——无论是河里的水，勺子上滑落的蜂蜜，还是承载负荷的钢梁——我们需要一个工具，能够精确地分离和量化这种形变行为。这个工具就是形变率张量。

运动的局部地图

要描述一个物体如何形变，仅仅知道单一点的速度 $\mathbf{v}$ 是不够的。我们需要知道当我们移动到邻近点时，速度如何变化。如果我们取两个无限接近的点，它们由一个微小的向量 $d\mathbf{x}$ 分隔，它们的速度差 $d\mathbf{v}$ 由速度梯度张量给出，通常记为 $L$ 。在坐标系中，这写为：

L_{ij} = \frac{\partial v_i}{\partial x_j}

你可以将 $L$ 看作是“运动的局部地图”。它包含了某一点紧邻区域内发生的所有信息。它告诉我们，当我们朝任何方向迈出一小步时，速度矢量如何变化。这一个数学对象混合包含了平移、旋转和形变的所有信息。我们下一步，也是最关键的任务，就是将它们分离开来。

伟大的分解：分离自旋与应变

这里蕴含着一段具有深远物理意义的数学魔法。任何方阵——我们的速度梯度 $L$ 正是如此——都可以唯一地分解为一个对称矩阵和一个斜对称（或反对称）矩阵之和。

L = D + W

在这个分解中：

$D$ 是对称部分，称为形变率张量（或应变率张量）。
$W$ 是斜对称部分，称为自旋张量（或涡量张量）。

它们的计算方法如下：

D = \frac{1}{2}(L + L^T) \quad \text{(对称)}

W = \frac{1}{2}(L - L^T) \quad \text{(斜对称)}

这不仅仅是一个数学技巧；它是两种不同物理行为的清晰分离。张量 $D$ 捕捉了所有纯粹的形状改变——拉伸、压缩和剪切。张量 $W$ 捕捉了纯粹的局部刚体旋转。一个运动是刚体运动当且仅当其形变率为零，即 $D = \mathbf{0}$ 。

让我们看看它的实际应用。

纯拉伸与压缩：对角项

考虑一个流动，其中每个粒子都以与其距离成正比的速度远离或朝向原点运动。一个简单的例子是速度场 $\mathbf{v} = (k_1 x_1, k_2 x_2, k_3 x_3)$ 。让我们计算速度梯度 $L$ ：

L = \begin{pmatrix} \frac{\partial v_1}{\partial x_1} & \frac{\partial v_1}{\partial x_2} & \frac{\partial v_1}{\partial x_3} \\ \frac{\partial v_2}{\partial x_1} & \frac{\partial v_2}{\partial x_2} & \frac{\partial v_2}{\partial x_3} \\ \frac{\partial v_3}{\partial x_1} & \frac{\partial v_3}{\partial x_2} & \frac{\partial v_3}{\partial x_3} \end{pmatrix} = \begin{pmatrix} k_1 & 0 & 0 \\ 0 & k_2 & 0 \\ 0 & 0 & k_3 \end{pmatrix}

这个矩阵已经是对称的了！这意味着自旋张量 $W = \frac{1}{2}(L - L^T)$ 为零。这个运动纯粹是形变，没有局部旋转。形变率张量就是 $D = L$ 。

D = \begin{pmatrix} k_1 & 0 & 0 \\ 0 & k_2 & 0 \\ 0 & 0 & k_3 \end{pmatrix}

对角分量 $D_{11}$ 、 $D_{22}$ 和 $D_{33}$ 直接告诉我们沿 $x_1$ 、 $x_2$ 和 $x_3$ 轴的拉伸速率（如果为正）或压缩速率（如果为负）。

如果我们把这些对角项加起来，会发生一件有趣的事情。这个和被称为张量的迹， $\text{tr}(D)$ 。事实证明， $\text{tr}(D) = \frac{\partial v_1}{\partial x_1} + \frac{\partial v_2}{\partial x_2} + \frac{\partial v_3}{\partial x_3} = \nabla \cdot \mathbf{v}$ 。这个量，即速度的散度，衡量的是体积膨胀的速率。如果 $\text{tr}(D)=0$ ，就像在许多液体流动中一样，这种运动被称为不可压缩的——我们微小流体立方体的体积不会改变，即使其形状可能被扭曲。

剪切与旋转：非对角项

那么像 $D_{12}$ 这样的非对角项呢？它们描述了剪切速率，即两条最初垂直的线之间夹角的变化率。

为了看清这一点，让我们对比一下流体力学中的两种经典流动。

简单剪切流： $\mathbf{v} = (ky, 0, 0)$ 。这描述了一种流体层相互滑动的流动，就像一副牌中的牌一样。速度梯度是 $L = \begin{pmatrix} 0 & k \\ 0 & 0 \end{pmatrix}$ 。让我们分解它：
$D = \frac{1}{2}\left( \begin{pmatrix} 0 & k \\ 0 & 0 \end{pmatrix} + \begin{pmatrix} 0 & 0 \\ k & 0 \end{pmatrix} \right) = \begin{pmatrix} 0 & k/2 \\ k/2 & 0 \end{pmatrix}$ $W = \frac{1}{2}\left( \begin{pmatrix} 0 & k \\ 0 & 0 \end{pmatrix} - \begin{pmatrix} 0 & 0 \\ k & 0 \end{pmatrix} \right) = \begin{pmatrix} 0 & k/2 \\ -k/2 & 0 \end{pmatrix}$
在这里， $D$ 和 $W$ 都非零。 $D$ 中非零的非对角项表示剪切形变，而 $W$ 非零则表示流体微元也在旋转。
平面驻点流： $\mathbf{v} = (kx, -ky, 0)$ 。这种流动描述了流体从 $y$ 方向流入，并沿着 $x$ 方向流出。速度梯度是 $L = \begin{pmatrix} k & 0 \\ 0 & -k \end{pmatrix}$ 。正如我们之前看到的，这是一个对称矩阵，代表纯拉伸。形变张量是 $D=L$ ，自旋张量是 $W=\mathbf{0}$ 。因此，一个流体微元在一个方向被拉伸，在另一个方向被压缩，但它不旋转。这是一种无旋流。

这个对比完美地说明了分解的威力：它将改变形状的部分（ $D$ ）与旋转的部分（ $W$ ）分离开来。自旋张量 $W$ 与流动的涡量 $\boldsymbol{\omega} = \nabla \times \mathbf{v}$ 直接相关，后者是物理学家衡量局部旋转的标准。

为何这很重要：形变的物理学

所以，我们有了一种将形变与旋转清晰分开的方法。为什么这如此重要？

首先，物理定律应该与你的观察角度无关。形变率 $D$ 在某种意义上是一个物理“实在”的量，而自旋 $W$ 不是。如果你在一辆转弯的车里，外面的世界似乎在围绕你旋转。你对自旋的测量取决于你自身的旋转。然而，形变率则不然。它是一个客观张量，意味着它的分量在不同观察者之间的变换方式简单而可预测，而自旋张量 $W$ 不具备这一特性。即使在一个以恒定速度移动的观察者的更简单情况下，形变率也保持不变，这巩固了它作为流动本身固有属性的地位。

其次，或许也是最重要的，大自然本身在涉及力和能量时也遵循这种划分。形变材料内部的力——应力，记为 $\sigma$ ——是为了抵抗形状变化而产生的。对于许多常见材料，如水或空气（所谓的牛顿流体），黏性应力与形变率成正比： $\sigma' = 2\mu D$ ，其中 $\mu$ 是黏度。注意，自旋 $W$ 根本没有出现！

这导致了一个关于能量的深刻结论。为使材料形变而做的机械功（通常会以热量形式耗散）的速率，由单位体积的功率 $p = \sigma : L$ 给出。因为应力张量 $\sigma$ 是对称的，它与斜对称的自旋张量 $W$ 的“乘积”总是零。因此，功率简化为：

p = \sigma : D

事实证明，大自然是一位出色的记账员。它精心地将花费在改变形状上的能量与仅仅是旋转的“能量”分离开来。你可以在转盘上旋转一桶水（刚体旋转），而不会使其升温。但如果你用勺子猛烈地搅拌它（产生剪切形变），你所做的功将因黏性而耗散为热量。

形变的自然坐标轴

张量 $D$ 的分量取决于我们选择的 $x,y,z$ 坐标系。但是形变本身并不关心我们的坐标轴。对于任何形变状态，都必定存在一组特殊的、“自然的”坐标轴。在这些方向上，材料线只被拉伸或压缩，没有剪切。这些被称为主应变率轴。

在数学上，这些方向是形变率张量 $D$ 的特征向量。沿这些方向的拉伸速率是相应的特征值，称为主应变率。找到这些特征值揭示了材料中伸长的最大和最小速率，为我们提供了最纯粹的形变物理图像，而与我们选择的坐标系无关。

一个统一的原则

形变率张量的概念是连续介质力学的基石。它在流体和固体的世界之间架起了一座美丽的桥梁。在流体动力学中， $D$ 描述了移动流体微元的瞬时形变率。在固体力学中，对于经历微小形变的材料，正是这个张量 $D$ 等于无穷小应变张量的物质时间导数，即 $D = \dot{\boldsymbol{\varepsilon}}$ 。

此外，对于大的、复杂的形变，形变率表现为更高级应变度量（如Cauchy-Green张量）在当前瞬间的变化率。从最简单的流动到最复杂的材料失效，形变率张量 $D$ 都处于问题的核心，优雅地捕捉了事物如何改变形状的基本运动学。

应用与跨学科联系

现在我们已经熟悉了形变率张量的数学机制，我们可能会想把它留在张量和索引的抽象世界里。但这将是一个天大的错误！因为这个张量不仅仅是一个数学上的奇物；它是一把钥匙，解锁了广阔的物理现象景观，一种描述事物如何改变形状的通用语言。它将你咖啡中奶油的漩涡与冰川缓慢而无情的爬行联系起来，将钢铁的锻造与我们身体的生长联系起来。让我们踏上一段旅程，看看这把万能钥匙适合用在何处。

流动的语言：本构律

想象你正在探测一种神秘物质。你让它形变，你拉伸和剪切它，并倾听它的响应。你问的问题是：“你正在如何形变？”——这个问题由形变率张量 $\mathbf{D}$ 精确地提出。该物质以内部力或应力的形式给出答案，由应力张量 $\boldsymbol{\sigma}$ 描述。问题 $\mathbf{D}$ 和答案 $\boldsymbol{\sigma}$ 之间的关系被称为本构律。这是材料独特的印记，它的个性。

最简单和最熟悉的个性是牛顿流体，比如水或空气。它的回答直截了当且呈线性：它产生的黏性应力与你使其形变的速率成正比。对于不可压缩流体，这种优雅的对话形式为 $\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu\mathbf{D}$ ，其中 $p$ 是背景压力， $\mathbf{I}$ 是单位张量，比例常数 $\mu$ 是我们熟悉的黏度。这种简单的线性响应是著名的Navier-Stokes方程的基石，它支配着极其广泛的流体动力学现象。

但大自然热爱多样性。许多物质具有更复杂的个性。它们是非牛顿的。想想番茄酱：它在瓶子里很稠，但当你用力摇晃时就很容易流动。它的黏度随形变率而变化。为了描述这类材料，我们需要更丰富的语言。Reiner-Rivlin流体模型正是这样做的，它假设应力不仅可能依赖于 $\mathbf{D}$ ，还可能依赖于它的平方 $\mathbf{D}^2$ 。这个非线性项的加入，为描述剪切稀化和剪切稠化等大量新行为提供了词汇，这对于研究聚合物、油漆和血液等复杂物质流动的流变学领域至关重要。

这种对话可以更加微妙。一些材料根据你推动它们的方向而有不同的响应。木材沿纹理比横穿纹理更容易劈开；这类材料是各向异性的。为了捕捉这一点，简单的标量黏度 $\mu$ 就不再足够了。我们必须将其提升为一个四阶黏度张量 $C_{ijkl}$ ，它像一本复杂的词典，将应变率的每个分量翻译成应力的一个分量： $\tau_{ij} = C_{ijkl} D_{kl}$ 。这种描述水平对于理解工程复合材料、液晶，甚至地球地幔中岩石的流动都至关重要。在所有情况下，形变率张量仍然是我们向材料提出的基本问题。

运动的代价：能量耗散

形变从来都不是免费的。推动流体、搅拌汤，或者仅仅是让它在管道中流动都需要做功。抵抗内部黏性力所做的功并没有丢失；它被转化为热能，使材料变暖。形变率张量准确地告诉我们这个“运动的代价”是多少。单位体积耗散的功率由优美而紧凑的表达式 $\dot{q} = \boldsymbol{\sigma} : \mathbf{D}$ 给出，这是应力张量和应变率张量之间的完全“握手”或缩并。这个原理不仅仅是一个学术上的好奇心；它解释了为什么摩天大楼或桥梁中的黏性阻尼器能够有效地将地震的剧烈能量转化为无害的热量，从而保护结构。

对于我们的老朋友牛顿流体，我们可以将其本构律代入这个功率方程。结果是一个简单而深刻的关系：耗散率是 $\Phi = 2\mu \mathbf{D}:\mathbf{D}$ 。这意味着能量损失的速率与黏度和形变率张量的“平方”成正比。这种二次依赖关系在直觉上是说得通的——以两倍的速度搅拌流体需要四倍的功率来克服黏性。这正是黏性是什么的核心：一种将机械能耗散为热量的机制。

从微观到宏观：隐藏的机制

到目前为止，我们一直将材料视为光滑、均匀的“东西”——一个连续介质。但我们知道这是一种错觉。物质是由原子组成的，我们看到的优美流动是微观尺度上疯狂、隐藏的舞蹈的结果。形变率张量在这两个世界之间提供了一座绝佳的桥梁。

考虑一个金属晶体。它的永久性或塑性形变不是平滑的屈服，而是离散的滑移事件的结果，其中原子平面沿着特定的晶体学方向相互滑动。每个滑移系都贡献了少量的剪切。材料工程师测量的宏观塑性应变率张量 $\dot{\boldsymbol{\varepsilon}}^p$ ，仅仅是所有这些微观贡献的总和。每个由其滑移面法线 $\mathbf{m}$ 和滑移方向 $\mathbf{s}$ 定义的滑移系，都贡献一个与其剪切速率 $\dot{\gamma}$ 和几何因子 $\frac{1}{2}(\mathbf{s} \otimes \mathbf{m} + \mathbf{m} \otimes \mathbf{s})$ 成比例的项。连续介质的张量作为无数原子量子跃遷的统计平均而出现。

在湍流的混乱中也展开着类似的故事。当流体运动得足够快时，其运动会分解成一个由各种尺寸的旋转涡流组成的漩涡。在计算机模拟中，我们只能承担计算大涡流的成本，而让小的、快速的涡流无法被解析。但这些小涡流并非无足轻重；它们非常有效地从大尺度运动中抽取能量并混合流体，其行为奇特地类似于一种增强的黏性。Boussinesq假设明确了这一类比：它将亚格子尺度涡流施加的应力建模为与已解析流动的形变率张量成正比， $\boldsymbol{\tau}_{sgs} = -2\nu_{sgs}\bar{\mathbf{D}}$ 。“涡黏度” $\nu_{sgs}$ 不是真正的材料属性，而是未解析混沌集体耗散效应的表示。形变率张量再次为建模一个复杂的多尺度现象提供了框架。

塑造我们的世界：从锻造钢铁到生长骨骼

有了这种深刻、多尺度的理解，我们不仅可以用形变率张量来描述世界，还可以用它来塑造世界。在制造业中，诸如锻造、轧制和挤压等过程都是为了对材料施加特定的形变，以达到期望的形状和内部结构。在像扭转挤压这样的先进技术中，材料受到巨大的剪切，以细化其晶体结构并增强其强度。为了量化这种过程的剧烈程度，工程师计算一个称为Von Mises等效应变率的标量值， $\dot{\bar{\varepsilon}} = \sqrt{\frac{2}{3} D_{ij} D_{ij}}$ ，它直接从形变率张量的分量导出。这个值就像一个“里程表”，记录了一个材料点所经历的总形变。

让我们不要忘记张量本身隐藏的美丽几何。在一个形变介质中的任何一点，无论流动多么复杂，都存在一组特殊的三个正交方向——主轴——材料沿着这些方向只经历纯拉伸或压缩，没有剪切。这种拉伸的速率是形变率张量的特征值，而这些轴是它的特征向量。找到这些轴使我们能够将任何复杂的形变分解为其最基本的组成部分，为流动的局部运动学提供了深刻的见解。

也许最令人惊讶的应用在于生物学领域。如果不是新物质的创造和生命组织的持续形变，生长又是什么呢？生长的连续介质力学使用形变率张量作为中心工具。它的迹 $\text{tr}(\mathbf{D}) = \nabla \cdot \mathbf{v}$ ，测量了局部的体积膨胀率。一个正在穿透土壤的植物根，或者一个为响应压力而重塑自身的骨骼，都是一个正在主动生成形变场的系统。在这些过程的模型中，形变率张量的迹可以直接与局部质量的生成率相关联。在这里，我们的张量超越了力学，成为生命本身语言的一部分。

从水的简单流动到晶体中原子的复杂舞蹈，从湍流的混沌到生长这一寂静而有序的过程，形变率张量提供了一种单一、统一的语言。它证明了数学捕捉物理现实本质的力量，并揭示了贯穿我们宇宙的深刻联系。