数学2 第十章 有限元法.pptVIP

第十章 有限元法 有限元法是求解偏微分方程问题的一种重要数值方法，它的基础分两个方面：一是变分原理，二是剖分插值。从第一方面看，有限元法是Ritz-Galerkin方法的一种变形。它提供了一种选取“局部基函数”的新技巧，从而克服了Ritz-Galerkin方法选取基函数的固有困难。从第二方面看，它是差分方法的一种变形。差分法是点近似，它只考虑在有限个离散点上函数值，而不考虑在点的邻域函数值如何变化；有限元方法考虑的是分段（块）的近似。因此有限元方法是这两类方法相结合，取长补短而进一步发展了的结果。在几何和物理条件比较复杂的问题中，有限元方法比差分方法有更广泛的适应性。 一、有限元方法解题分析 为了说明应用有限元方法的解题步骤，以及每一步骤中的要点，下面我们 以两点边值问题为例进行具体分析。 考虑两点边值问题 其中 我们将从Ritz法和Galerkin法两种观点出发，导出 解边值问题（1.1）、（1.2）的线性有限元方法。 （一）从Ritz法出发建立有限元方程 1、写出Ritz形式的变分问题 由变分原理可知，与边值问题（1.1）（1.2） 等价的变分问题是：求 其中 积分表达式（1.3）是应用有限元法求解（1.1）、（1.2） 式的出发点。 2、区域剖分 剖分原则与差分法相同，即将求解区域剖分成 若干个互相连接，且不重叠的子区域，这些子区域 称为单元。单元的几何形状可以人为选取，一般是 规则的，但形状与大小可以不同。对于一维情形最 为简单： 将求解区域 剖分成若干个子区间，其节点为 每个单元 的长度为 单元在区间中分布的疏密程度或单元尺寸的大小， 可根据问题的物理性质来决定，一般来说，在物理量变 化剧烈的地方，单元尺寸要相对小一些，排列要密一些。 3、确定单元基函数 有限元法与Ritz-Galerkin方法的主要区别之一，就 在于有限元方法中的基函数是在单元中选取的。由于各 个单元具有规则的几何形状，而且可以不比考虑 边界条件的影响，因此在单元中选取基函数可遵 循一定的法则。 设 为 的有限维子空间，它的元素为 要构造 ，只需构造单元基函数 。构造单元基 函数应遵循如下原则： （1）、每个单元中的基函数的个数和单元中的 节点数相同，每个节点分别对应一个基函数，本例 中，单元 有两个节点，因此基函数有两个。 （2）基函数应具有下面的性质： 其中 是单元节点序号为k的节点。 若取 为线性函数，则按上述原则，可将 中的基函数取为 显然， 中任一函数 可以表示为基函数 的 线性组合，即 （1.4） 其中， 是 在节点上的值，即 在单元 上， 表示为 可见，单元中的近似函数由单元基函数线性组合 产生，全区域的近似函数由各个单元的近似函数 叠加而成。 从以上可以看出， 是满足下列条件的所有 函数 的集合： 故 是 的一个n维子空间，称为试探函数空间 称为试探函数。 4、有限元方程的形成 与Ritz法一样，以 代替 ，在 上解泛函数 （1.3）的极小问题。 将（1.5）代入（1.3），得 令 便得到确定 的线性代数方程组 称（1.7）为有限元方程。 显然，只要我们分别算出 及 就可以求解（1.7）。 但在工程计算中，并不是按照上述步骤形成 有限元方程的，而是首先建立单元有限元特 征式（称这一过程为单元分析），然后再将 单元的有限元特征式进行累加，合成为总体 有限元方程（这一过程称为总体合成）。 下面分步分析具体的计算方法。 第一步：单元分析。注意到 我们来计算单元 上的积分。为讨论方便，作变换 并引入记号 则在 上， 可写成 或写成 （1.10） 其中， 而 可表示为 式中， 于是有 这里， 称为单元刚度矩阵，其中 对（1.8）式右端第二项积分，同样有 式中， 称 为单元“荷载”向量。 根据以上分析，便有 这样，我们就得到了单元有限元特征式的一般表示形式： 第二步：总体合成 总体合成就是将单元的有限元特征式进行累加， 合成为总体有限元方程。这一过程实际上是将 单元有限元特征式中的系数矩阵逐个累加，合成 为总体系数矩阵（称为总刚度矩阵