电磁场有限元分析复习课程.ppt

工程电磁场数值分析（4）（电磁场有限元法）;第4章 电磁场有限元法（Finite Element Method, FEM）;加权余量法回顾： 对算子方程 用 作为该方程的近似解（试探解）： 代入方程得余量：;设L为线性算子，代入 ，得;利用有限元法求解一维边值问题： （1）单元剖分 如图5个单元，6个节点 （2）选取基函数 ;（3）方程离散 （计算系数阵 [K] 和右端项 [b]） 基函数 Ni 只是一阶可导 的，不能严格满足微分方 程，称为“弱解”。;（3）方程离散;（3）方程离散;总体方程;（4）求解方程;与有限差分法（FDM）相比，有限差分法是对点的离散，得到一系列离散点上的解；而有限元（FEM）是对区域的离散（单元），尽管所求的是节点上的自由度，但它的解在场域中每一个点上都有定义。 所以，即是有限元节点上的解是精确的，有限元的整个解仍然是近似的。好的数据处理技术可以从该近似解中提取更精确的分析结果。 线性单元中，如果所求的自由度是电位j，单元中的电场 E是场量；节点上的 E 取邻近单元的平均。;计算系数阵是有限元分析的主要工作量。所涉及到的积分，如果不是解析可积的，通常要用到数值积分。其中最常用的数值积分方法是Gauss数值积分。;n=4 x(1)= 0.861136311594053d0 x(2)= 0.339981043584856d0 w(1)= 0.347854845137454d0 w(2)= 0.652145154862546d0 n=5 x(1)= 0.906179845938664d0 x(2)= 0.538469310105683d0 x(3)= 0.0d0 w(1)= 0.236926885056189d0 w(2)= 0.478628670499366d0 w(3)= 0.568888888888889d0 n=6 x(1)= 0.932469514203152d0 x(2)= 0.661209386466265d0 x(3)= 0.238619186083197d0 w(1)= 0.171324492379170d0 w(2)= 0.360761573048139d0 w(3)= 0.467913934572691d0;为提高有限元分析精度，有两种方法： 其一：增加节点，细化网格——称为h方法。 其二：增加有限元的阶数——称为p方法。 ;一些补充说明： 二阶单元;一些补充说明： 三阶单元;h方法和p方法的求解精度;作业：要独立完成，凡雷同者没分！！;以二维静电场泊松方程的求解为例。; 场域离散 二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。;基函数 ;三角形单元内的基函数 设三角形三个顶点处待求函数值分别为u1, u2, u3。如果单元足够小，可以采用线性近似，将单元内任意p点的u(x,y)表示为 ; ;记住我们的任务 —寻找基函数;在积分 中，对于确定的 i，j的有效取值为i 本身以及与节点i相联的周围节点，积分的有效区域为以i、j 为公共节点的所有三角形单元 ，在这些单元中Ni、Nj才有交叠。; 这些积分可以分单元进行。例如对右图所示的局部编码，K01、K00以及b0的计算公式为： ;以下把单元e的贡献记为;单元分析：计算单元内积分对系数阵和右端项元素的贡献。;故 ，; 写成一般形式，若一个三角形三个顶点编号为i, j, m（逆时针顺序），则;再看边界部分：; 由于单元很小，做单元分析时通常可以取 f (e) 为常数值（可以认为等于三个顶点上的平均值）。因此;上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而在实际编程中，更有效率的是以单元为序，逐个计算单元系数阵[K(e)]，然后合成整体系数阵[K]。单元系数阵[K(e)]定义为 设 i, j, m 是节点的整体编号，元素Kij在整体矩阵中的实际位置是第i行、j列；因此 必须合成到整体矩阵的第i行、j列元素上。; 整体矩阵合成：; 通过上述过程，对于一个“正常”的内部节点就建立起了一个代数方程。“非正常”的节点包括：媒质交界面衔接条件和场域边界条件。 ; 对于静电场问题，媒质分界面衔接条件为; 由于有限元方法能够自动满足媒质交界面条件，因此有限元法特别适合于处理多层复杂媒质问题。这是其它方法无可比拟的。;第一类边界条件（强加边界条件）;第二类边界条件（自然边界条件）;有限元方法的推导过程虽然看起来有些