1. 课程导入与有限元弱形式概览
有限元方法（Finite Element Method, FEM）是工程数值分析中最主流的求解框架。它的核心思想可以概括为一句话：把连续的求解域离散成有限个单元，在每个单元内用简单的形函数近似未知场，再通过积分形式的等效方程组装成代数方程组求解。理解这句话的关键在于两个层面——离散化（几何与场的离散）和等效方程（从强形式到弱形式）。本课程后续所有推导与代码，都围绕这两个层面展开。
为什么需要弱形式？考虑一维杆的轴向受力问题。设杆长 L、截面积 A、弹性模量 E、分布轴向载荷 f(x)，位移场 u(x) 满足强形式方程：
d/dx(EA du/dx) + f = 0, x ∈ (0, L)
边界条件：u(0) = 0（固定端），EA du/dx |_{x=L} = P（自由端受集中力）强形式要求 u 二阶可导，这在解析上很苛刻；而工程中载荷、截面常分段变化，甚至存在集中力，严格满足二阶可导的解往往不存在。弱形式通过引入试探函数 v（也称权函数、虚位移），将方程两边乘以 v 并在全域积分，再用分部积分把二阶导数转移到 v 上，从而把对 u 的光滑性要求从二阶降到一阶。推导过程如下：
∫₀ᴸ [d/dx(EA du/dx)] v dx + ∫₀ᴸ f v dx = 0
分部积分：∫₀ᴸ EA (du/dx)(dv/dx) dx = [EA du/dx · v]₀ᴸ + ∫₀ᴸ f v dx
代入边界：v(0)=0，EA du/dx|_{L}=P
得弱形式：∫₀ᴸ EA (du/dx)(dv/dx) dx = ∫₀ᴸ f v dx + P·v(L)这个式子就是有限元离散的出发点。左端是双线性型 a(u,v)，右端是线性型 l(v)，最终离散后得到 K u = F。可以看到，弱形式天然容纳了自然边界条件（集中力 P 直接进入载荷向量），而本质边界条件（u(0)=0）需要在试探空间和刚度矩阵中强制施加。这是强形式与弱形式最本质的区别之一。
上图把本节的主线串起来：强形式经“乘权函数—积分—分部积分”得到弱形式，弱形式再经“形函数插值—单元积分—组装”得到代数系统。同时下方的一维杆示意给出了课程项目的具体对象：3 个线性单元、4 个节点、左端固定、右端受集中力、杆身受分布载荷。后续章节将围绕这个对象，逐步完成从数学推导到程序实现的完整链条。
弱形式推导的一般流程可以归纳为四步，但每一步都有需要留意的细节：
- 写出强形式与边界条件：明确控制方程、本质边界条件（位移约束）和自然边界条件（力/通量）。一维杆中，u(0)=0 是本质边界条件，EA u'(L)=P 是自然边界条件。
- 选取试探函数与权函数空间：试探函数 u 需满足本质边界条件，权函数 v 在本质边界处取零。这一步决定了离散后哪些自由度被约束。
- 乘权函数并分部积分：把高阶导数转移到权函数上，同时把自然边界条件以边界积分项的形式带入。分部积分的方向和符号必须与坐标系、外法向约定一致。
- 离散化与组装：用形函数展开 u 和 v，得到单元刚度矩阵与单元载荷向量，再按节点编号组装成全局方程，最后施加本质边界条件求解。
课程项目采用一维杆问题作为贯穿案例，原因在于它足够简单，可以手算验证；又足够完整，涵盖了有限元实现的全部关键环节：单元刚度矩阵推导、数值积分、全局组装、边界条件处理、求解与后处理。建议在开始后续章节前，先手算一个两单元杆的刚度矩阵，例如取 L=1 m、A=1×10⁻⁴ m²、E=200 GPa、P=1000 N、f=0，得到节点位移后与解析解 u(x)=Px/(EA) 对比。这个手算练习能帮助你在写代码时快速定位符号错误或组装错误。
需要强调的是，弱形式不是数学上的“近似妥协”，而是有限元方法的严格数学基础。它把求解偏微分方程转化为在有限维子空间中寻找满足积分等效关系的近似解，Galerkin 方法正是取权函数与形函数同空间的特例。理解了这一点，后续遇到二维、三维、热传导、结构动力学等问题时，推导思路是一致的，只是积分维度和本构关系不同。
本节建立了课程的整体框架：有限元的核心是离散化加等效积分；强形式到弱形式的关键操作是乘权函数与分部积分；弱形式到程序实现的关键操作是形函数插值与单元组装。接下来的章节将沿着这条主线，先完成一维杆弱形式的严格推导，再逐步进入形函数构造、数值积分、程序结构与结果验证。