华中师范大学 hahakity
有限元算法(Finite Element Method,简称 FEM)是一种非常流行的求解偏微分方程的数值算法。有限元被广泛应用于结构受力分析、复杂边界的麦克斯韦方程求解以及热传导等问题。这一节介绍有限元方法的基本原理,以及如何用 Python 从头实现一个有限元算法,数值求解麦克斯韦方程。
学习内容
- 筑基:加权残差法 (Weighted Residual Method)
- 心法:有限元与有限差分算法的区别
加权残差法(Weighted Residual Method )
理解了加权残差法,有限元算法就理解了一半。
这里参考文献【1】中例子,介绍加权残差法,为有限元算法筑基。
举例:求解
, 使其满足如下微分方程,
初始条件
,
求解区域
。
这个方程的解析解是
,
在计算机上展开到无穷阶,就能给出方程的精确解。现实中显然只能考虑有限项,后面忽略的项称作截断误差。假设截断到二阶,
代入原方程产生残差
,
越接近0,表示近似解越接近精确解。
明知道精确解
不在
这三个函数支撑起的函数空间中,我们仍可做出改进。
使用待定系数展开,