局部系数矩阵元素的计算.ppt

1. 有限元法 1. 有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 2. 一维有限元法 3. 二维有限元法 3. 二维有限元法 3. 二维有限元法 3. 二维有限元法 3. 二维有限元法 3. 二维有限元法 3.二维有限元法 3. 二维有限元法 3. 二维有限元法 3. 二维有限元法 3. 二维有限元法 3.二维有限元法 该问题用变分法从一般化的偏微分方程出发，也可得到与加权余量法一样的局部系数矩阵表达式 3. 二维有限元法 用变分法的一般化目标泛函得问题的泛函 该问题的局部系数矩阵元素值： 3. 二维有限元法 结合剖分情况，根据上式计算各单元局部系数矩阵： 3. 二维有限元法 * * 上一讲，利用加权余数法和变分法将偏微分方程转化为代数方程组求解 该矩阵方程包括系数矩阵、激励源矩阵和边界矩阵，而计算这些矩阵的元素时，常常用到分部积分法。如果为了计算精度而选取很多个尝试函数，那么计算这些为数众多的分部积分既十分复杂又很费时间，并且很难用计算机进行数值计算。 因此，我们需要寻找一个改进的方法来简化计算，并设法利用计算机进行处理，有限元法就是其中的一种。 通过尝试函数的选取，近似解满足1类边界条件， 基本思想： 在有限元方法中，场域被分割成许多很小的子区域，通常称为“单元”或“有限元”。 对所有子区域进行独立的处理和运算，便对一个整体问题进行局部化处理。 通过选取恰当的尝试函数，使每个单元的计算都变得非常简单，经过对每个单元重复而简单的计算，再将其结果总和起来，便可以得到用整体矩阵表达的整个区域的解， 这一整体矩阵又常常是稀疏短阵，可以更进一步简化和加快求解过程。由于计算机非常适合于重复性的计算和处理过程，所以整体矩阵的形成过程很容易使用计算机来实现， 以满足帕松方程的有源静电场为例： 由对称性，将上述问题化为一维问题，用有限元方法求解 场域单元化 （剖分） 场域(0，1)分割成四个“单元”，e1、e2、e3、e4。 单元大小和数量与精度和计算量有关。 处理复杂场域和激励源带来极大方便。 作为一种数值计算方法，有限元法并非用来寻求问题的解析解。实际上，很多工程问题目前都无法找到解析解。有限元的作用就在于求解分布场的势函数在每个节点上的近似值，而势函数在单元的其它位置的值可用插值来表示：如果采用线性捅值的方法来表示分布势函数，则称为一阶有限元法，如果采用高阶插值法表示分布势函数，则称为高阶有限元法，通常一阶插值即可。 本例，场域分割成4个单元，5个节点，求场域内电势分布，转化为求5个节点的电位即可。 场域内其它点（各单元内）的电位，由5个节点电位来插值表示。(一阶插值、高阶插值) 对一 维场域来说，单元就是一个线段； 对二维场域，有限元单元形状可为二角形、矩形等，单元形状对有限元的简化有影响，通常为三角形 单元化后，整体场域积分转化为局部单元求和 每个单元可独立处理 把两个要求解的量联系起来，有限元中令待定系数就是节点电位，当然尝试函数要重新确定 总体系数矩阵化为了局部系数矩阵的求和 这个场域的求解变为求解剖分后节点上的电位 重新确定尝试函数的思想 尝试函数本质上是在场域内对真解的近似。 以前的尝试函数是针对整个场域来选取的（随意性大，有技巧经验因素） 现在整个场域被剖分为N个单元，每个单元内的真解由该单元尝试函数线性组合来近似，（即尝试函数也局部化，） 当单元足够多（剖分足够细）时每单元内的尝试函数可以非常简单（线段） 尝试函数：一维一阶有限元尝试函数为直线段；二维一阶为平面 (形函数) 任意单元内的形函数只与该单元节点坐标有关，与其它单元无关 某单元内尝试函数确定 某单元内尝试函数的表达（局部坐标） 任意单元内尝试函数的形式相同，单元内才不为0，构成的局部系数矩阵对整体矩阵有贡献 在考查的单元之外，尝试函数为0，即i, j不属于同一单元时，构成的局部系数矩阵为0，只有当i, j＝i+1属于同一单元时，构成的局部系数矩阵才不为0，才对整体矩阵有贡献。这样每个单元的计算分开、独立了，且每个单元的计算有是相同的、重复的，便于计算机求解。 局部系数矩阵的计算 局部系数矩阵元素的计算 q是问题的激励 局部系数矩阵的计算 任意单元e的局部系数矩阵，只与该单元的长度（坐标）有关，与其它单元无关，该矩阵的i,i+1的行列上有值，其它位置为0。 其它所有单元的系数矩阵结构形式完全相同，只是不同单元的长度le不同而已。 这样极便于计算机作重复