工程电磁场-第7章-二维泊松方程的有限元法课件

11223344556677889910101111121213131414151516167.3二维泊松方程的有限元法177.3二维泊松方程的有限元法17181819192020212122222323242425252626272728282929以二维静电场泊松方程的求解为例。

目标：依据加权余量法，利用分域基，建立离散的代数方程组，即确定系数{Kij}

和{bi}。30以二维静电场泊松方程的求解为例。 目标：依据加权余量法，利用场域离散 二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。

单元：互不重叠，覆盖全部场域；每个单元内介质是 单一、均匀的。

节点：网格的交点，待求变量的设置点。 该步骤需要记录的信息： 节点编号、节点坐标 节点属性(激励源、是否边界等) 单元编号 单元节点编号 单元介质 31场域离散 单元：互不重叠，覆盖全部场域；每个单元内介质是 三角形单元内的基函数 设三角形三个顶点处待求函数值分别为u1,u2,u3。如果单元足够小，可以采用线性近似，将单元内任意p点的u(x,y)表示为 代入三个顶点的坐标和函数值，可以解出a、b、c。得到32三角形单元内的基函数 代入三个顶点的坐标和函数值，可以解出

单元节点的编号按逆时针方向排列！其中，33 单元节点的编号按逆时针方向排列！其中，33记住我们的任务—寻找基函数对比可得基函数Ni常被称为插值函数或者形状函数，具有以下性质：（1）是插值的；（2）（3）在相邻单元的公共边界上，Ni是连续的，从而通过Ni构造的逼近函数也是连续的。34记住我们的任务对比可得基函数Ni常被称为插值函数或者形状函数在积分中，对于确定的i，j的有效取值为i本身以及与节点i相联的周围节点，积分的有效区域为以i、j为公共节点的所有三角形单元，在这些单元中Ni、Nj才有交叠。计算系数阵

35在积分

这些积分可以分单元进行。例如对右图所示的局部编码，K01、K00以及b0的计算公式为：

计算系数阵

36 这些积分可以分单元进行。例如对右图所示的局部编码，K01、以下把单元e的贡献记为这样，就有

每个或的计算都在具体的单元内单独考虑（称为单元分析）。37以下把单元e的贡献记为这样，就有 每个

由于单元很小，做单元分析时通常可以取f(e)为常数值（可以认为等于三个顶点上的平均值）。因此

右端项元素：38 由于单元很小，做单元分析时通常可以取f(e)为常数值上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而在实际编程中，更有效率的是以单元为序，逐个计算单元系数阵[K(e)]，然后合成整体系数阵[K]。单元系数阵[K(e)]定义为

设i,j,m是节点的整体编号，元素Kij在整体矩阵中的实际位置是第i行、j列；因此必须合成到整体矩阵的第i行、j列元素上。

单元矩阵：39上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而

整体矩阵合成：40 整体矩阵合成：40第一类边界条件（强加边界条件）第一类边界节点是指边界上函数值已知。因此处理方法是，合成整体系数阵之后，将该节点所在行的主元素置1，其它元素均置零，同时将右端项中对应元素设为已知函数值。要保持对称性；有更简便的做法41第一类边界条件（强加边界条件）第一类边界节点是指边界上函数值42143244345446547648749850951105211531254135514561557167.3二维泊松方程的有限元法587.3二维泊松方程的有限元法17591860196120622163226423652466256726682769287029以二维静电场泊松方程的求解为例。

目标：依据加权余量法，利用分域基，建立离散的代数方程组，即确定系数{Kij}

和{bi}。71以二维静电场泊松方程的求解为例。 目标：依据加权余量法，利用场域离散 二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。

单元：互不重叠，覆盖全部场域；每个单元内介质是 单一、均匀的。

节点：网格的交点，待求变量的设置点。 该步骤需要记录的信息： 节点编号、节点坐标 节点属性(激励源、是否边界等) 单元编号 单元节点编号 单元介质 72场域离散 单元：互不重叠，覆盖全部场域；每个单元内介质是 三角形单元内的基函数 设三角形三个顶点处待求函数值分别为u1,u2,u3。如果单元足够小，可以采用线性近似，将单元内任意p点的u(x,y)表示为 代入三个顶点的坐标和函数值，可以解出a、b、c。得到73三角形单元内的基函数 代入三个顶点的坐标和函数值，可以解出

单元节点的编号按逆时针方向排列！其中，74 单元节点的编号按逆时针方向排列！其中，33记住我们的任务—寻找基函数对比可得基函数Ni常被称为插值函数或者形状函数，具有以下性质：（1）是插值的；（2）（3）在相邻单元的公共边界上，Ni是连续的，从而通过Ni构造的逼近函数也是连续的。75记住我们的任务对比可得基函数Ni常被称为插值函数或者形状函数在积分中，对于确定的i，j的有效取值为i本身以及与节点i相联的周围节点，积分的有效区域为以i、j为公共节点的所有三角形单元，在这些单元中Ni、Nj才有交叠。计算系数阵

76在积分

这些积分可以分单元进行。例如对右图所示的局部编码，K01、K00以及b0的计算公式为：

计算系数阵

77 这些积分可以分单元进行。例如对右图所示的局部编码，K01、以下把单元e的贡献记为这样，就有

每个或的计算都在具体的单元内单独考虑（称为单元分析）。78以下把单元e的贡献记为这样，就有 每个

由于单元很小，做单元分析时通常可以取f(e)为常数值（可以认为等于三个顶点上的平均值）。因此

右端项元素：79 由于单元很小，做单元分析时通常可以取f(e)为常数值上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而在实际编程中，更有效率的是以单元为序，逐个计算单元系数阵[K(e)]，然后合成整体系数阵[K]。单元系数阵[K(e)]定义为

设i,j,m是节点的整体编号，元素Kij在整体矩阵中的实际位置是第i行、j列；因此必须合成到整体矩阵的第i行、j列元素上。

单元矩阵：80上述以节点为序的分析过程对于有限元原理的说明是易于理解的