大连理工大学数学科学学院 - 大连理工大学精品课程.pdf

微分方程数值解法 信息与计算科学专业核心课 大连理工大学数学科学学院 第3章 Galerkin有限元法 Galerkin有限元法 有限元法，实质上就是Ritz-Galerkin法。它和传统的Ritz- Galerkin法的主要区别在于，它应用样条函数方法提供了一种选取 “局部基函数”或“分片多项式空间”的新技巧，从而在很大程度 上克服了古典Ritz-Galerkin法选取基函数时的固有的困难。 有限元法首先成功地用于结构力学和固体力学，以后又用于流 现在，有限元和差分法一样，已 体力学、物理学和其他工程科学。 经成为求解偏微分方程，特别是椭圆型偏微分方程的一种有效的数 值方法。 有限元的基本问题可归纳为一下六个方面： 把问题转化为变分问题 选定单元形状，对求解域作剖分 构造基函数或单元形状函数 有限元的基本问题： 有限元的基本问题： 形成有限元(R-G)方程 提供有限元方程的有效解法 收敛性及误差估计 总之，有限元方法通过剖分插值构有限子空间V 的基函数的 n 方法保持了古典Ritz-Galerkin法的优点,克服了它的不足。 有限元 方法尤其适用于区域的形状复杂、定解条件中含有第二、三类 边界条件、方程的系数有间断点等情况。因此，有限元方法和 限差分法并列成为有效的求解偏微分方程定解问题的数值方法。 例如，两点边值问题 d ⎛ du ⎞ Lu =− p +qu f , ⎜ ⎟ a x b , dx ⎝ dx ⎠ u a 0, u′ b 0 ( ) ( ) 其中p ∈C1 (I),p (x)≥p min ＞0,q ∈C0 (I),f ∈L2 (I),I=[a,b] 下面我们将从Ritz法和Galerkin法两种观点出发,导出解上述 边值问题的有限元方法 1.1 从Ritz法出发 首先，构造广义解空间，即Sobolev空 1 ( ) H I E 的有限维子空间 v 。 h 为此，对求解区间进行网格剖分，节点为 a x x x x b L L 0 1 i n