第 1 章 · 从弹簧拼接开始：离散系统的有限元

有限元法最好的入口不是抽象的数学，而是一根弹簧。整根弹簧的力学只有一条规律：$F = ku$（力 = 刚度 × 伸长）。把结构看成大量弹簧的拼接，有限元的一切——刚度矩阵、组装、边界条件——都能从这里长出来。

单根弹簧的刚度矩阵

一根弹簧两端节点 1、2，位移 $u_1, u_2$，轴向力平衡给出：

这个 $2\times2$ 的 $\mathbf{k}$ 就是单元刚度矩阵。注意它的两个特征：对称（$k_{12}=k_{21}$，功的互等）与奇异（两行相加为零——刚体平移不产生力，所以行列式为零）。

拼接：组装的真面目

两根弹簧串联共用中间节点，把两个 $2\times2$ 按"节点归属"叠加进总刚度矩阵的对应行列，就得到总阵 $\mathbf{K}$。组装 = 按公共节点自由度对号入座求和，没有任何更神秘的东西。得到的方程组：

施加边界条件（固定端 $u=0$）删行删列后，$\mathbf{K}$ 不再奇异，方程可解。

为什么这条路线通向一切

弹性体可以看成无限多根"弹簧"的连续集合，热传导问题里温度对应位移、热流对应力——只要物理能用"势 + 梯度通量"描述，有限元结构完全一样。这就是为什么同一个 FEM 程序框架能算结构、传热、渗流、电磁。

老师的提醒

• 动手算一个两弹簧串联问题（$k_1, k_2$、一端固定、另一端受力），手算出 $3\times3$ 总阵并解出位移——这是理解组装最快的方法。
• 奇异阵报错时，先找"哪里没约束住"，而不是怀疑求解器。

自检

• 单元刚度矩阵为什么必然对称、奇异？
• 组装的数学操作是什么？
• 哪些物理场可以套用同一套有限元框架？

更深一层：从"弹簧"到"连续体"只差一个极限

第 1 节用弹簧搭起有限元的地基，值得把这个类比的极限说透。把弹性杆想象成 $n$ 根短弹簧串联：总刚度 $k/n$（串联变软），每根弹簧的伸长 $\Delta/n$。令 $n \to \infty$，离散位移序列变成连续函数 $u(x)$，节点力变成分布载荷——有限维线性系统连续化为微分方程。有限元法做的事恰好相反：把微分方程退回有限维线性系统，用矩阵代数这个人类最擅长的工具处理它。这条"离散-连续-再离散"的往返不是循环论证，而是认识结构：连续理论告诉我们解长什么样（存在性、正则性），离散方法告诉我们怎么把它算出来。

刚度矩阵的三条代数性质都有物理孪生：对称性来自功的互等（或最小势能原理的二次型）；奇异性来自刚体运动不被约束；稀疏性来自局部作用（每个节点只与力学意义上的邻居耦合）——稀疏性是有限元能算百万自由度问题的根本原因，满阵时代一个 $10^4$ 自由度问题就需要 GB 级内存，稀疏存储把它压到 MB 量级。

推导细读：最小势能原理——另一条更"力学"的路

第 1 章用虚功原理引入，其实还有一条平行的推导路线：弹性系统在平衡时总势能取最小值。总势能 $\Pi(u) = \frac12\int EA(u')^2dx - \int bu\,dx - P\,u(L)$ 是位移函数的泛函；有限元把 $u$ 限制在形函数张成的有限维子空间里，于是求泛函极小变成求多元函数极小：

——与伽辽金弱形式殊途同归（对称算子时两者完全等价）。这条路线的额外红利是变分不等式的直觉：极小值处的二阶变分正定 ⟹ $\mathbf{K}$ 正定 ⟹ 解唯一且共轭梯度法适用。学过优化的读者会立刻认出：有限元离散 = 在函数空间里做 Ritz-Galerkin 投影，求解 = 解一个凸二次问题——两门学科的词汇表在这里合流。

在真实世界里：一次真实结构分析的完整决策记录

以某支架的静力校核为例记录真实工作流：几何简化（去掉工艺倒角、保留应力考核圆角——圣维南允许远端简化）；单元选择（实体六面体二次单元，减缩积分，因为板厚方向只有三层单元怕锁死）；网格策略（圆角半径内至少 3 层单元、渐变比 1.3、其余粗网格）；约束模拟（螺栓孔用弹簧单元模拟预紧刚度而非刚性固定——过刚约束会虚增孔边应力 15% 以上）；载荷（面压力而非集中力，避免加载点伪应力）；结果验收（网格加密 30% 后考核点应力变化 3.2% < 5% 内控标准，通过）。这十条决策里只有一条涉及"解方程"，其余全是建模判断——有限元工程师的核心技能分布由此可见。

常见误区与辨析

• "自由度越多越准"：在奇异点（点载荷、尖角）无限加密只会无限追高伪应力；收敛性研究的对象应是光滑区的代表性量。
• "对称结构取一半建模没问题"：约束对称面上法向位移成立，但若后续要算扭转或不对称载荷，对称模型直接失效——建模前先列全部工况。
• "软件算出来的云图就是应力分布"：单元平均、节点磨平各有假设；考核应力必须用未磨平的积分点/单元值并检查该区域的离散质量。

专题深入：矩阵的"零能模式"与约束设计

刚度矩阵的奇异性（刚体模态）在有限元里有个更精细的概念谱系：刚体模态（平移与旋转——无应变的整体运动）与零能模式（沙漏——单元内部的"假变形"）。两者的区别在于：刚体模态必须被约束（不约束则刚度阵奇异、方程无解）、沙漏模式必须被控制（不控制则解被污染）。奇异矩阵报警的排障清单按出现频率：刚体模态未约束（约束不足）→ 部件未连接（共节点或接触）→ 零刚度单元（$\det J \approx 0$ 的压扁单元或材料参数漏填）→ 约束冲突（乘子法的矛盾方程）→ 接触初始穿透/分离抖动。90% 的奇异报警落在前两条——排障的第一步永远是"哪里没约束住"，不是怀疑求解器。

数学细节展开：从虚功到刚度的"一步一步"

从虚功原理到刚度矩阵的推导值得逐步走完，因为它是有限元的"出生证明"。第一步：虚功方程 $\int \delta\boldsymbol\epsilon^T\boldsymbol\sigma\,dV = \int\delta\mathbf{u}^T\mathbf{f}\,dV$（内力虚功 = 外力虚功）。第二步：离散化——$\mathbf{u} = \mathbf{N}\mathbf{d}$（形函数插值位移）、$\boldsymbol\epsilon = \mathbf{B}\mathbf{d}$（应变-位移矩阵）。第三步：本构——$\boldsymbol\sigma = \mathbf{D}\boldsymbol\epsilon$。第四步：代入并约去任意虚位移 $\delta\mathbf{d}$：

每一步都是一个"物理翻译成代数"的动作——虚功是平衡、离散是逼近、本构是材料、约去任意性是"对一切虚位移成立"。这个推导链在传热、渗流、电磁的所有位移型有限元里一字不改。

工程案例细节：一个支架静力校核的完整清单

真实支架校核的完整清单（每条都是血泪教训）：几何简化（去工艺倒角、留考核圆角——圣维南允许远端简化）；单元选择（实体六面体二次单元，减缩积分，薄壁方向三层单元防锁死）；网格策略（圆角半径内至少 3 层单元、渐变比 1.3）；约束模拟（螺栓孔用弹簧单元模拟预紧刚度而非刚性固定——过刚约束虚增孔边应力 15% 以上）；载荷（面压力而非集中力，避免加载点伪应力）；结果验收（网格加密 30% 后考核点应力变化 3.2% < 5% 内控标准，通过）。十条决策里只有一条涉及"解方程"，其余全是建模判断——有限元工程师的核心技能分布由此可见。

常见误区与辨析（补充）

• "刚度阵对称来自材料对称"：不——来自功的互等（或最小势能的二次型），与各向异性材料无关（非对称摩擦/陀螺项才是对称性的例外）。
• "奇异 = 无解"：奇异 + 相容右端 = 无穷多解（差一个刚体运动）；纯奇异在无约束时才是"无解报警"。