位移梯度与应变的关系
在固体力学中，描述物体变形状态最直观的方式是通过位移场。设物体在参考构型（未变形状态）下的位置矢量为 $\mathbf{X}$，在变形后的当前构型下的位置矢量为 $\mathbf{x}$。位移矢量 $\mathbf{u}$ 定义为两者之差，即 $\mathbf{u}(\mathbf{X}) = \mathbf{x}(\mathbf{X}) - \mathbf{X}$。
为了量化局部变形，我们需要考察位移在空间中的变化率。位移梯度张量 $\mathbf{H}$（通常记为 $\nabla \mathbf{u}$）正是描述了位移矢量在参考构型中的空间变化率。其分量表示为 $H_{ij} = \frac{\partial u_i}{\partial X_j}$。位移梯度张量包含了物体变形的全部信息，但它混合了刚体旋转和纯变形两种效应。为了分离出纯粹的变形部分，我们需要引入应变张量。
小变形假设下的线性应变
在工程实践中，大多数结构分析基于小变形假设，即位移及其梯度都非常小（$|\nabla \mathbf{u}| \ll 1$）。在此假设下，位移梯度张量 $\mathbf{H}$ 可以分解为对称部分和反对称部分：
$$ \mathbf{H} = \frac{1}{2}(\mathbf{H} + \mathbf{H}^T) + \frac{1}{2}(\mathbf{H} - \mathbf{H}^T) $$
其中：
- 对称部分 $\mathbf{\epsilon} = \frac{1}{2}(\mathbf{H} + \mathbf{H}^T)$ 称为小应变张量（或柯西应变张量）。
- 反对称部分 $\mathbf{\omega} = \frac{1}{2}(\mathbf{H} - \mathbf{H}^T)$ 称为旋转张量，代表刚体旋转。
小应变张量 $\mathbf{\epsilon}$ 直接反映了物体局部的拉伸、压缩和剪切变形。其分量具有明确的物理意义：
- 对角分量 $\epsilon_{11}, \epsilon_{22}, \epsilon_{33}$ 分别表示沿坐标轴方向的正应变（线应变）。
- 非对角分量 $\epsilon_{12}, \epsilon_{13}, \epsilon_{23}$ 表示工程剪应变的一半，即 $\epsilon_{ij} = \frac{1}{2}(\gamma_{ij})$。
例如，在二维平面问题中，若位移场为 $u_x = ax + by$，$u_y = cx + dy$，则位移梯度为：
$$ \mathbf{H} = \begin{bmatrix} a & b \ c & d \end{bmatrix} $$
对应的线性应变张量为：
$$ \mathbf{\epsilon} = \begin{bmatrix} a & \frac{b+c}{2} \ \frac{b+c}{2} & d \end{bmatrix} $$
这里，$a$ 和 $d$ 是正应变，$\frac{b+c}{2}$ 是剪应变分量。
有限变形下的应变度量
当变形较大时，线性应变张量不再适用，因为几何非线性效应显著。此时，必须使用基于变形梯度张量 $\mathbf{F} = \mathbf{I} + \mathbf{H}$ 的非线性应变度量。最常用的有两种：
格林-拉格朗日应变张量（Green-Lagrange Strain Tensor, $\mathbf{E}$）：
定义在参考构型上，表达式为：
$$ \mathbf{E} = \frac{1}{2}(\mathbf{F}^T\mathbf{F} - \mathbf{I}) $$
该张量在材料本构关系（如超弹性材料）中非常常用，因为它基于初始构型，便于处理大变形问题。阿尔曼西应变张量（Almansi Strain Tensor, $\mathbf{e}$）：
定义在当前构型上，表达式为：
$$ \mathbf{e} = \frac{1}{2}(\mathbf{I} - \mathbf{F}^{-1}\mathbf{F}^{-T}) $$
该张量适用于应力分析，因为应力通常定义在当前构型上。
应变不变量与物理意义
无论使用哪种应变张量，其本质都是度量线元长度和夹角的变化。对于任意两点间的微元线段，其变形前后的长度变化可以通过应变张量精确计算。
- 正应变：反映纤维长度的相对变化。例如，$\epsilon_{11} > 0$ 表示沿 $X_1$ 方向的纤维被拉伸。
- 剪应变：反映直角的改变。例如，$\epsilon_{12} \neq 0$ 表示原本垂直的 $X_1$ 和 $X_2$ 方向纤维之间的夹角发生了变化。
在实际计算中，从位移场获取应变是有限元分析（FEA）的核心步骤之一。通过节点位移插值函数，可以计算单元内的位移梯度，进而得到应变场。随后，结合本构方程（如胡克定律 $\boldsymbol{\sigma} = \mathbf{D}:\boldsymbol{\epsilon}$），即可求解应力分布。
总结
位移梯度是连接运动学与力学的桥梁。在小变形情况下，其对称部分直接给出线性应变；在大变形情况下，需通过变形梯度构造格林-拉格朗日或阿尔曼西应变张量。正确选择应变度量对于建立准确的力学模型至关重要，尤其是在涉及橡胶、生物组织等超弹性材料或金属成型等塑性大变形问题时。理解位移梯度与应变的数学关系，是掌握固体力学基础理论的关键一步。