1. 课程导入与单元矩阵组装概览

有限元分析的本质，是把一个连续体离散成有限个单元，在每个单元内部用形函数近似位移场，再由能量原理或虚功原理得到单元层面的力学方程，最后把所有单元方程"拼"成一个整体方程组求解。这条主线可以概括为：结构离散 → 单元分析 → 整体组装 → 引入边界条件 → 求解 → 后处理。其中"单元分析"和"整体组装"是承上启下的两个关键环节，而本课程聚焦的正是后者——如何把一个个孤立的单元刚度矩阵，按照节点编号的对应关系，准确地叠加进整体刚度矩阵。

很多初学者能推导出单根杆或单个梁单元的刚度矩阵，却在面对多单元结构时卡住：明明每个单元矩阵都算对了，组装出来的整体方程却解出错误的位移。问题往往不在单元推导，而在组装阶段的"对号入座"——哪个单元的哪个自由度，应该累加到整体矩阵的哪一行哪一列。这正是本课程要系统解决的问题。

组装在整体刚度方程中的位置

整体刚度方程的标准形式为 K·U = F，其中 K 是整体刚度矩阵，U 是整体节点位移向量，F 是整体节点载荷向量。组装环节的输入是各单元的刚度矩阵 k^e 和单元节点自由度与整体自由度的映射关系，输出是组装完成的 K 和 F。用数学语言描述，组装就是按映射关系做累加：

K[LM(i), LM(j)] += k_e[i, j]
F[LM(i)]          += f_e[i]

这里的 LM（Location Matrix，定位矩阵）记录了单元第 i 个局部自由度对应的整体自由度编号。理解并正确构造 LM，是组装工作的核心。组装完成后，K 通常是对称、稀疏、带状分布的；在引入边界条件之前，它还是奇异的（存在刚体位移），这也是后续章节要处理的内容。

一维问题的典型场景

一维问题是学习组装的最佳载体，因为自由度少、映射关系直观，可以把注意力集中在"叠加规则"本身。课程将覆盖三类典型对象：

• 弹簧单元：每个节点仅 1 个轴向自由度，单元刚度矩阵为 2×2，形式为 k·[[1,−1],[−1,1]]。它是最纯粹的组装练习，适合建立 LM 映射的第一直觉。
• 杆单元（桁架杆）：同样每节点 1 个轴向自由度，但刚度系数为 EA/L，涉及材料与几何参数。多杆串联、并联、带支座约束的组装都在此展开。
• 梁单元：每个节点有挠度和转角 2 个自由度，单元矩阵为 4×4，组装时自由度编号从"每节点 1 个"变为"每节点 2 个"，是理解自由度扩展与 LM 构造的关键一步。

以一根由 3 个等长杆单元串联、总长 3L、弹性模量 E、截面积 A 的受拉杆为例（示例参数：E=200 GPa，A=100 mm²，L=500 mm，端部受集中力 P=10 kN）。单元①连接节点 1、2，单元②连接节点 2、3，单元③连接节点 3、4。每个单元的 LM 分别为 [1,2]、[2,3]、[3,4]。组装时，单元①的 k22 与单元②的 k11 都落在整体自由度 2 上，必须相加；这正是图中橙色与绿色块所强调的叠加关系。若漏加或错位，整体矩阵将不再对称，求解结果也会失真。

学习路线图

本课程按"概念—规则—实现—验证"的递进组织：先建立组装在整体方程中的位置感，再掌握 LM 定位矩阵的构造方法，随后分别对弹簧、杆、梁三类单元完成手工组装与程序化组装，最后通过对称性检查、刚体位移检验和与解析解对比来验证组装正确性。每一节都会给出可复算的小规模算例，建议读者在纸上或代码中同步演算，把"看会"变成"做对"。

需要强调的是，组装不是机械的"填表"，而是对结构拓扑关系的数值表达。节点编号顺序会直接影响 K 的带宽，进而影响求解效率；对于大规模问题，合理的节点编号（如减小带宽的 Cuthill-McKee 排序）本身就是工程实践的一部分。本课程先保证"组装正确"，再讨论"组装高效"，为后续二维、三维问题的组装打下可迁移的基础。