导读:Zienkiewicz 有限元体系的起点不是单元,而是「弱形式」。理解了从强形式到弱形式的两步操作,有限元的一切推导都变成机械劳动。
第一步:乘权函数积分
把 Poisson 方程 $-\nabla^2 u = f$ 两边乘上权函数 $v$ 并在域内积分,用分部积分(格林公式)把二阶导数降为一阶:$\int \nabla u \cdot \nabla v\,dV = \int f v\,dV + $ 边界项。
这一步的收获有二:对解的连续性要求从 $C^1$ 降到 $C^0$;边界条件中的 Neumann 项自然出现——自然边界条件「免费」进入弱形式。
第二步: trial 空间近似
把 $u$ 用基函数展开 $u \approx \sum N_j a_j$,权函数取与基函数相同(伽辽金选择),得到线性方程组 $K a = F$,其中 $K_{ij} = \int \nabla N_i \cdot \nabla N_j\,dV$ 就是刚度阵。
弱形式的本质:把「逐点满足方程」放松为「在加权平均意义下满足」。牺牲逐点严格性,换来可解性与对复杂几何的适应力。
为什么这个顺序不能反
先设计单元再凑理论的做法,会在遇到混合单元、非线性问题时全面崩溃。Zienkiewicz 的路线是:弱形式 → 插值 → 单元矩阵,每一步都可检查。
常见误区
- 忘记对本质边界(Dirichlet)做强制处理,导致刚度阵奇异。
- 以为弱形式是「近似」,实际上它是与强形式等价的数学陈述(解足够光滑时)。
吃透自检
- 能独立完成从强形式到弱形式的两步推导
- 能说明自然边界条件在弱形式中出现的位置
- 能写出刚度阵元素的定义并计算一维两个线性单元的例子