科学仿真

从强形式到弱形式:伽辽金方法的两步推导

导读: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)做强制处理,导致刚度阵奇异。
  • 以为弱形式是「近似」,实际上它是与强形式等价的数学陈述(解足够光滑时)。

吃透自检

  • 能独立完成从强形式到弱形式的两步推导
  • 能说明自然边界条件在弱形式中出现的位置
  • 能写出刚度阵元素的定义并计算一维两个线性单元的例子

← 科学仿真 书目精读列表 经典研读