引言
线弹性力学研究小变形固体在卸载后能够完全恢复变形的力学行为，是结构强度分析最直接的理论工具。其核心问题是建立应力与应变之间的线性映射——广义胡克定律。对一般各向异性材料，这一映射包含 21 个独立弹性常数；而对工程中绝大多数金属材料，两个常数（如杨氏模量 $E$ 与泊松比 $\nu$）即可完整描述。本文从应力应变的小变形描述出发，推导广义胡克定律的结构与各向同性退化形式，梳理弹性常数之间的换算关系，并以电阻应变花测量为例展示理论在工程测试中的完整应用链条。
应力与应变的小变形描述
在小变形假设下，位移梯度远小于 1，应变由位移场的对称梯度定义：
$$\varepsilon_{ij} = \frac{1}{2}\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right)$$
对角分量 $\varepsilon_{11}$、$\varepsilon_{22}$、$\varepsilon_{33}$ 表示沿坐标方向的线应变，非对角分量（工程剪应变 $\gamma_{xy} = 2\varepsilon_{12}$）表示直角的角度改变量。应力张量 $\sigma_{ij}$ 为对称二阶张量，共 6 个独立分量。线弹性问题共有 15 个未知场量（3 个位移、6 个应变、6 个应力），由 3 个平衡方程、6 个几何方程与 6 个物理方程（本构关系）封闭求解。
广义胡克定律与各向异性
最一般的线性本构关系把每个应力分量表示为全部应变分量的线性组合：
$$\sigma_{ij} = C_{ijkl},\varepsilon_{kl}$$
四阶弹性张量 $C_{ijkl}$ 表观上有 81 个分量，但应力与应变的对称性使其独立分量减至 36 个；若进一步要求变形过程存在应变能势函数（超弹性假设），则 $C_{ijkl} = C_{klij}$，独立常数最终只有 21 个。工程上常用 Voigt 记号把应力、应变排成 6 维列向量，本构关系写成 $6\times 6$ 矩阵形式。材料的对称性越高，独立常数越少：正交各向异性材料（如纤维增强复合材料层板）有 9 个，横观各向同性材料（如沉积岩）有 5 个，各向同性材料仅 2 个。
各向同性弹性体的本构方程
各向同性要求本构关系在任意坐标旋转下形式不变，这迫使 $C_{ijkl}$ 只能由 Kronecker 符号构造，最终得到：
$$\sigma_{ij} = \lambda,\varepsilon_{kk}\delta_{ij} + 2\mu,\varepsilon_{ij}$$
其中 $\lambda$、$\mu$ 为 Lamé 常数，$\mu$ 即剪切模量 $G$。其逆形式用工程常数表示为：
$$\varepsilon_{ij} = \frac{1+\nu}{E}\sigma_{ij} - \frac{\nu}{E}\sigma_{kk}\delta_{ij}$$
该式还可分解为体积响应与畸变响应两部分：$\sigma_{kk} = 3K\varepsilon_{kk}$ 与 $s_{ij} = 2\mu e_{ij}$（$s_{ij}$、$e_{ij}$ 分别为偏应力与偏应变）。这一分解表明：静水压力只引起体积改变，偏应力只引起形状改变，两种响应分别由体积模量 $K$ 与剪切模量 $\mu$ 独立控制——这正是塑性理论把屈服行为与偏应力（而非静水压力）联系起来的物理背景。
应变能密度为 $W = \frac{1}{2}\sigma_{ij}\varepsilon_{ij}$，各向同性时 $W = \frac{1}{2}\lambda\varepsilon_{kk}^2 + \mu\varepsilon_{ij}\varepsilon_{ij}$，其正定性要求 $E > 0$、$-1 < \nu 0, \mu>0$），这也解释了为何橡胶类材料 $\nu \approx 0.5$ 表现为几乎不可压缩。五个常用弹性常数中只有两个独立，换算关系如下表。
| 常数 | 符号 | 用 $E,\nu$ 表示 | 用 $\lambda,\mu$ 表示 |
|---|---|---|---|
| 剪切模量 | $\mu$（即 $G$） | $\dfrac{E}{2(1+\nu)}$ | $\mu$ |
| Lamé 第一常数 | $\lambda$ | $\dfrac{E\nu}{(1+\nu)(1-2\nu)}$ | $\lambda$ |
| 体积模量 | $K$ | $\dfrac{E}{3(1-2\nu)}$ | $\lambda + \dfrac{2\mu}{3}$ |
| 杨氏模量 | $E$ | $E$ | $\dfrac{\mu(3\lambda+2\mu)}{\lambda+\mu}$ |
| 泊松比 | $\nu$ | $\nu$ | $\dfrac{\lambda}{2(\lambda+\mu)}$ |
需要区分平面应力与平面应变两种退化情形：薄板面内问题取 $\sigma_{33} = 0$（平面应力），厚截面长结构取 $\varepsilon_{33} = 0$（平面应变），两者的等效刚度矩阵相差 $1-\nu^2$ 因子，在有限元建模中选错单元类型会带来系统性的应力误差。
工程案例分析：电阻应变花测量应力
电阻应变片只能测量沿栅丝方向的线应变，无法直接得到剪应变与应力。为确定平面应力场，工程上使用三向应变花：沿 $0^\circ$、$45^\circ$、$90^\circ$ 三个方向布置应变片，利用应变转轴公式
$$\varepsilon_\theta = \frac{\varepsilon_x+\varepsilon_y}{2} + \frac{\varepsilon_x-\varepsilon_y}{2}\cos 2\theta + \frac{\gamma_{xy}}{2}\sin 2\theta$$
由三个测量值反解 $\varepsilon_x$、$\varepsilon_y$、$\gamma_{xy}$，再代入平面应力形式的胡克定律换算应力。某钢制构件表面测得 45° 直角应变花读数如下表，材料弹性模量 $E = 210$ GPa，泊松比 $\nu = 0.30$，许用应力 $[\sigma] = 160$ MPa。
| 测点方向 | 读数（με） | 说明 |
|---|---|---|
| $0^\circ$ | 620 | 即 $\varepsilon_x$ |
| $45^\circ$ | 180 | 用于解剪应变 |
| $90^\circ$ | −140 | 即 $\varepsilon_y$ |
由 $45^\circ$ 片读数得 $\gamma_{xy} = 2\varepsilon_{45} - \varepsilon_x - \varepsilon_y = -120 \mu\varepsilon$。以下程序完成应变到应力、再到主应力与 von Mises 等效应力的全过程换算：
import numpy as np
E, nu = 210e3, 0.30 # MPa, 钢材
G = E / (2 * (1 + nu)) # 剪切模量
# 45° 直角应变花读数（单位应变）
e0, e45, e90 = 620e-6, 180e-6, -140e-6
ex, ey = e0, e90
gxy = 2 * e45 - e0 - e90 # 工程剪应变 γ_xy
# 平面应力假设：由应变求应力
sx = E / (1 - nu**2) * (ex + nu * ey)
sy = E / (1 - nu**2) * (ey + nu * ex)
txy = G * gxy
print(f"σx = {sx:7.2f} MPa, σy = {sy:7.2f} MPa, τxy = {txy:7.2f} MPa")
# 面内主应力（自由表面 σ3=0）与 von Mises 等效应力
avg = (sx + sy) / 2
R = np.sqrt(((sx - sy) / 2)**2 + txy**2)
s1, s2 = avg + R, avg - R
vm = np.sqrt(s1**2 - s1 * s2 + s2**2)
print(f"σ1 = {s1:7.2f} MPa, σ2 = {s2:7.2f} MPa")
print(f"von Mises 等效应力 = {vm:.2f} MPa")计算结果：$\sigma_x = 133.38$ MPa，$\sigma_y = 10.62$ MPa，$\tau_{xy} = -9.69$ MPa；主应力 $\sigma_1 = 134.15$ MPa、$\sigma_2 = 9.85$ MPa，von Mises 等效应力约 129.5 MPa，低于许用值 160 MPa，该测点强度储备系数约 1.23。整个链条——应变花读数、转轴公式、平面应力胡克定律、主应力合成——正是广义胡克定律在结构试验监测中的标准用法；反过来，有限元后处理中由节点应变重构应力的过程与之完全同构。
结论与展望
广义胡克定律以 21 个独立常数刻画最一般的线性弹性材料，而各向同性假设将其简化为 $\lambda$、$\mu$ 两个常数，使应力应变分析具备了工程可操作性。弹性常数换算表与平面应力/平面应变的区分是日常使用中最易出错的两个环节。当变形进入有限范围或材料出现各向异性、温度耦合时，需要回到更一般的本构框架；当应力超过弹性极限后，线弹性关系失效，塑性与断裂理论将接过强度分析的接力棒，这也是后续文章的主题。
发表回复