python简单的计算方法_用python实现简单的有限元方法（二）

华中师范大学 hahakity

第一部分介绍了加权残差法(Weighted Residual Method)并引出有限元算法中经常出现的伽辽金(Galerkin)法。简单说明了有限元与有限差分法的区别：有限差分是对微分算子做差分近似，有限元是对待求函数做基函数展开近似。这一节继续介绍伽辽金有限元算法。

这里特意强调伽辽金有限元算法，是因为还有另一种基于泛函和变分的 Ritz 有限元算法。跳过 Ritz 变分法，直接讲解伽辽金法入门比较简单。

学习目标Galerkin 有限元法

分片基底函数展开

二阶微分算子的降阶

一维泊松方程：已知电荷密度求电势

预备知识加权残差法 (Weighted Residual Method)

Galerkin 有限元法

假设要求解的微分方程有统一的格式，

其中

是微分算子，

是要求解的场，

是源。

将

做函数基底展开 (使用爱因斯坦求和规则，省略

求和符号)，

注意此时

为基底函数，人工设定，为已知量，

为待定系数。

选取基底函数

作为伽辽金法中的检验函数，则对第 i 个检验函数，加权残差为，

其中

是全域，

,

可以将

写成矩阵乘积形式，此处用

表示列向量，

是

构成的矩阵，

是

构成的列向量，

是

构成的列向量。

Galerkin 有限元法令残差

为0，通过解