第六章 有限元法.pdfVIP

第6章 Galerkin有限元法 Ritz-Galerkin方法在实际应用中会遇到许多原则性的困难： （1）基函数的选取。 基函数必须满足本质边值条件，对于代数 或三角多项式为基函数而言，当求解区域的形状不规则时难于满足 本质边值条件，进而难于构造子空间V 。 n （2）Ritz-Galerkin方程的形成。 ( , a)ϕ ϕ 在实际计算时, 和 i j ( , )f ϕj 都是在区域上的积分,计算量大,而且系数矩阵 a(ϕ ,ϕ ) i j n×n 一般是稠密矩阵，无法节省计算量和存储量。例如，若取n=100, 则要计算约5000（对称情形）到10000（不对称情形）。 （3）求解Ritz-Galerkin方程。 系数矩阵 a(ϕ ,ϕ ) 的条件数 i j n×n 大，无论用迭代法或消元法都会出现数值不稳定的问题。 Galerkin有限元法 有限元法，实质上就是Ritz-Galerkin法。它和传统的Ritz- Galerkin法的主要区别在于，它应用样条函数方法提供了一种选取 “局部基函数”或“分片多项式空间”的新技巧，从而在很大程度 上克服了古典Ritz-Galerkin法选取基函数时的固有的困难。 有限元法首先成功地用于结构力学和固体力学，以后又用于流 现在，有限元和差分法一样，已 体力学、物理学和其他工程科学。 经成为求解偏微分方程，特别是椭圆型偏微分方程的一种有效的数 值方法。 有限元的基本问题可归纳为一下六个方面： 把问题转化为变分问题 选定单元形状，对求解域作剖分 构造基函数或单元形状函数 有限元的基本问题： 有限元的基本问题： 形成有限元(R-G)方程 提供有限元方程的有效解法 收敛性及误差估计 总之，有限元方法通过剖分、插值构造有限子空间V 的基函 n 数的方法，保持了古典Ritz-Galerkin法的优点,克服了它的不足。 有限元方法尤其适用于区域的形状复杂、定解条件中含有第 二、三类边界条件、方程的系数有间断点等情况。因此，有限元 方法和有限差分法并列成为有效的求解偏微分方程定解问题的数 值方法。 例如，两点边值问题 d ⎛ du ⎞ Lu =− p +qu f ,