二维静电场有限元matlab,实验2 MATLAB有限元计算

MATLAB电磁场有限元计算

实验目的：

1、了解有限元算法的原理，熟练运用MATLAB环境的PDE工具。 2、熟练运用PDE工具分析简单的电磁场边值问题。

实验内容：

一、 有限元简介

在电磁场的计算中, 仅对那些具有最简单边界条件和场域几何形状规则的问题才有解析解, 多数问题的求解必须用数值计算的方法,其场域分布的数值计算内容是学习难点。本实验将有限元法和Matlab 结合起来对电磁场教学中的电位分布问题进行计算。结果表明使用Matlab对有限元分析编程中的矩阵进行处理,程序设计清晰简便,易于理解和实现。

node(节点) 有限元法是以变分原理和剖分插值为基础的一种数值计算方法,其基本思想是将场域方程等价为一个条件变分问题,然后由条件变分问题对场域进行剖分离散为方程组进行求解。对于一个电场来说,其储能总是趋于最小,这样变分法的泛函和电场的储能就联系起来了。对于边界为L 的无源空气介质二维静电场中,一个封闭场域S 内的等价能量泛函可以写为:

element(单

在有限元分析中,将所研究的区域 S划分成有限的n 个三角形网格单元。 对应m个节点, ds 为单元e的面积。对任意三角形单元 e 中任一点的电位可以认 为由该三角形的三个节点(分别设为i、 j、 k) 上的电位u 随该点坐标x、 y 变化而线性确定。 因此, 对于单元e 构造插值函数:

其中ah 称为形状函数。那么有插值函数的一阶偏导数为:

从而得到能量函数We:

则将单元e中的能量函数We 对每一个节点电位ul ( l = i, j , k)求一阶偏导数, 得:

表示为矩阵形式有：

然后进行总体合成, 将各单元的能量函数对同一节点的电位一阶偏导数相加, 获得所要求解的线性方程组。 由以上分析,可知在该场域内电场有限元数学模型为: 式中U 为n 个节点处的待求电位, K 为n 阶矩阵。最后进行强加边界条件处理, 消去已知电位节点在系数矩阵中所在的