应变张量与几何变形
在连续介质力学的宏伟蓝图中，理解物体如何在受力作用下发生形状和尺寸的改变是核心任务之一。无论是设计承受高载荷的航空发动机叶片，还是模拟复杂的地质构造运动，精确描述物质的几何变形都是建立物理方程的第一步。本文将深入探讨连续介质力学中描述变形的核心数学工具——应变张量，并剖析其背后的几何意义与计算方法。
在连续介质力学中，我们假设物质是连续分布的，忽略其微观的原子间隙。当物体受到外力或温度场作用时，其内部的质点会从初始位置运动到新的空间位置。
设物质在变形前的初始构型（Reference Configuration）中的位置坐标为 $X$，变形后的当前构型（Current Configuration）中的位置坐标为 $x$。质点的运动可以描述为一个映射：
$$x = \chi(X, t)$$
在这个过程中，物体的变形通常包含两部分：刚体运动（整体平移和整体旋转）以及纯变形（拉伸、压缩和剪切）。刚体运动不会引起物体内部的相对位移，因而不产生内力和应力。我们所关心的应变，本质上是定量描述“纯变形”的物理量。
质点位移与位移梯度
为了量化变形，我们引入位移矢量 $u$，它定义为质点在变形前后的位置差：
$$u(X, t) = x(X, t) - X$$
通过对位移场求导，可以得到位移梯度张量 $\nabla u$。其分量形式在笛卡尔坐标系下表示为：
$$u_{i,j} = \frac{\partial u_i}{\partial X_j}$$
位移梯度包含了局部旋转和局部变形的信息。为了分离出纯变形，我们需要构造能够自动滤除刚体旋转的数学量——应变张量。
核心数学工具：格林应变与阿曼西-阿尔曼应变
在有限变形（大变形）理论中，通常使用拉格朗日描述（基于初始构型）或欧拉描述（基于当前构型）。
有限变形格林应变张量（Green-Lagrangian Strain Tensor, $E$）：
该张量基于初始构型，能够客观地描述大变形。其定义与变形梯度张量 $F = \frac{\partial x}{\partial X}$ 相关：
$$E = \frac{1}{2} (F^T F - I)$$
其中 $I$ 为单位张量。格林应变张量在刚体旋转时保持为零，完美地排除了刚体运动的干扰。阿曼西-阿尔曼应变张量（Almansi Strain Tensor）：
与格林应变相对，该张量基于当前构型，常用于流体力学或大变形欧拉分析中。
工程应用中的简化：无穷小应变张量
在许多工程实际问题中（例如常规金属材料的弹性变形），物体受力产生的位移和变形是非常微小的，即位移梯度的模远小于 1（$|\nabla u| \ll 1$）。在这种情况下，非线性项可以忽略不iphatic，从而进入小变形理论（Infinitesimal Strain Theory）的范畴。
此时，非线性的格林应变张量退化为线性的无穷小应变张量（或称工程应变张量）$\varepsilon$：
$$\varepsilon_{ij} = \frac{1}{2} \left( \frac{\partial u_i}{\partial X_j} + \frac{\partial u_j}{\partial X_i} \right)$$
在分量形式下：
- 线应变（法向应变）：对角线分量（如 $\varepsilon_{11}$）表示沿坐标轴方向的长度相对变化率。
- 剪应变：非对角线分量（如 $\varepsilon_{12}$）的一半表示原本相互垂直的两个方向之间的夹角变化。
应变张量的物理意义与应用全景
应变张量作为一个二阶对称张量，具有以下重要性质：
- 坐标无关性：虽然应变的分量会随着坐标系的旋转而改变，但其主应变（Principal Strains）和不变量（Invariants）是客观存在的物理属性，与坐标系的选择无关。
- 主应变与主方向：在任何变形点，总可以找到三个互相垂直的方向，在这三个方向上只有伸缩而没有剪切变形，这三个方向称为主方向，对应的拉伸率称为主应变。
在现代科学与工程中，应变张量的计算是连接宏观受力与材料本构关系（如广义胡克定律）的桥梁。无论是通过有限元分析（FEA）预测汽车碰撞时的结构失效，还是在地球物理学中模拟板块挤压引起的岩石形变，准确把握应变张量与几何变形的数学本质，都是进行精确力学建模不可或缺的基础。