第二章 有限元基础理论回顾：从弹性力学到离散化
各位同学，欢迎来到第二章。说实话，很多搞光机结构的人，一听到“弹性力学”、“虚功原理”就开始头疼。我当年也一样，觉得这些理论离实际工程太远。
但做了十几年项目后，我越来越发现：不懂这些基础，你连有限元软件报错都看不懂。有一次我算一个反射镜支撑结构，结果应力集中位置完全算错了。后来回头一查，是边界条件没理解透——说白了，就是弹性力学里的平衡方程没吃透。
这一章，咱们就把这些“硬骨头”啃下来。我会尽量用工程语言讲，不搞纯数学推导。
2.1 弹性力学三大基本方程
弹性力学研究的是：物体在外力作用下，内部应力、应变、位移之间的关系。这三大方程，就是描述这三者关系的“宪法”。
2.1.1 平衡方程（力的平衡）
平衡方程描述的是：物体内部任意一点，各个方向上的力必须平衡。用数学语言说，就是应力梯度与体积力相平衡。
在直角坐标系下，三维问题的平衡方程长这样：
∂σx/∂x + ∂τxy/∂y + ∂τxz/∂z + fx = 0
∂τyx/∂x + ∂σy/∂y + ∂τyz/∂z + fy = 0
∂τzx/∂x + ∂τzy/∂y + ∂σz/∂z + fz = 0其中 σ 是正应力，τ 是剪应力，f 是体积力（比如重力）。
2.1.2 几何方程（变形与位移）
几何方程描述的是：应变与位移之间的关系。说白了，就是物体怎么“变形”的。
对于小变形问题，几何方程是线性的：
εx = ∂u/∂x
εy = ∂v/∂y
εz = ∂w/∂z
γxy = ∂u/∂y + ∂v/∂x
γyz = ∂v/∂z + ∂w/∂y
γzx = ∂w/∂x + ∂u/∂z这里 u、v、w 分别是 x、y、z 方向的位移，ε 是正应变，γ 是剪应变。
2.1.3 物理方程（应力-应变关系）
物理方程，也叫本构方程。它描述的是材料本身的力学行为。对于各向同性的线弹性材料，就是大家熟悉的胡克定律：
εx = (1/E)[σx - ν(σy + σz)]
εy = (1/E)[σy - ν(σx + σz)]
εz = (1/E)[σz - ν(σx + σy)]
γxy = τxy/G
γyz = τyz/G
γzx = τzx/G其中 E 是弹性模量，ν 是泊松比，G 是剪切模量，三者关系为 G = E/[2(1+ν)]。
| 材料 | E (GPa) | ν | 典型应用 |
|---|---|---|---|
| 铝合金 6061 | 68.9 | 0.33 | 镜筒、结构件 |
| 钛合金 TC4 | 110 | 0.34 | 高刚度轻量化件 |
| 殷钢 | 141 | 0.28 | 低热膨胀镜架 |
| 光学玻璃 BK7 | 81 | 0.21 | 透镜、棱镜 |
2.2 虚功原理与最小势能原理
这两个原理，是有限元法的数学基础。你想想看，直接求解那15个偏微分方程有多难？尤其是复杂几何形状。所以我们需要一个“等效”的方法。
2.2.1 虚功原理
虚功原理说的是：如果一个物体处于平衡状态，那么对于任意微小的虚位移，外力做的虚功等于内力做的虚功。
数学表达式：
∫V σ·δε dV = ∫V f·δu dV + ∫S t·δu dS左边是内力虚功（应力在虚应变上做的功），右边是外力虚功（体积力和面力在虚位移上做的功）。
这个原理的好处是：把微分方程变成了积分方程。积分方程对函数的光滑性要求更低，更容易用数值方法求解。
2.2.2 最小势能原理
最小势能原理是虚功原理的推论。它说：在所有可能的位移场中，真实位移场使系统的总势能取最小值。
总势能 Π = 应变能 U - 外力势能 W
用变分法表示：δΠ = 0
这个原理直接给出了有限元法的求解框架：
- 假设一个位移场（用形函数表示）
- 计算总势能 Π
- 对位移求变分，令 δΠ = 0
- 得到刚度矩阵方程 [K]{u} = {F}
2.3 有限元离散化思想——以杆单元为例
好了，理论讲完了。现在咱们看看有限元到底是怎么工作的。我选最简单的杆单元（1D 线性单元）来演示。
2.3.1 为什么从杆单元开始？
杆单元只有两个节点，每个节点一个自由度（轴向位移）。它的刚度矩阵是 2×2 的，手算都能算。你想想看，如果连杆单元都搞不明白，怎么去理解复杂的四面体单元？
2.3.2 杆单元的推导过程
假设一个杆单元，长度 L，截面积 A，弹性模量 E，两个节点编号 1 和 2。
第一步：假设位移场
用形函数 N 表示单元内任意点的位移：
u(x) = N1(x)·u1 + N2(x)·u2
其中：
N1(x) = 1 - x/L
N2(x) = x/L第二步：计算应变
由几何方程：ε = du/dx = (u2 - u1)/L
第三步：计算应力
由物理方程：σ = E·ε = E·(u2 - u1)/L
第四步：应用最小势能原理
单元应变能：U = (1/2)∫V σ·ε dV = (1/2)·(EA/L)·(u2 - u1)²
对 u1、u2 求偏导，得到单元刚度矩阵：
[k] = (EA/L) · [ 1 -1
-1 1 ]第五步：组装全局刚度矩阵
把每个单元的 [k] 按照节点编号“对号入座”，叠加到全局矩阵 [K] 中。
2.3.3 从杆单元到一般单元
杆单元的推导思路，可以推广到任何单元类型：
| 步骤 | 杆单元 | 一般单元 |
|---|---|---|
| 1. 位移假设 | 线性形函数 | 多项式/等参形函数 |
| 2. 应变-位移 | B 矩阵（1×2） | B 矩阵（维度取决于单元） |
| 3. 应力-应变 | D 矩阵（标量 E） | D 矩阵（6×6 或更复杂） |
| 4. 单元刚度 | [k] = ∫ Bᵀ·D·B dV | 相同公式，数值积分 |
| 5. 全局组装 | 直接叠加 | 相同原理 |
你看，核心公式 [k] = ∫ Bᵀ·D·B dV 是通用的。区别只在于 B 矩阵和 D 矩阵的维度，以及积分方式（解析积分 vs 数值积分）。
2.4 本章知识体系图
下面这张图，把本章的核心逻辑串起来了。建议你保存下来，后面学完整个课程再回来看，会有更深的理解。
2.5 本章小结
这一章的内容，是有限元分析的理论基石。我建议你：
- 理解三大方程：不用背公式，但要明白它们描述的是什么物理关系
- 掌握变分原理：知道 [K]{u}={F} 是怎么来的
- 吃透杆单元推导：这是理解所有单元类型的起点