内容发布更新时间 : 2024/5/12 20:06:06星期一 下面是文章的全部内容请认真阅读。

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 阶矩阵。最后进行强加边界条件处理, 消去已知电位节点在系数矩阵中所在的行和列, 得到简化后的方程,继而可以对电位进行求解。流程框图如下图所示。

二、静电场仿真

静态场满足上方基本方程，式中 D 为电位移, ??为电荷密度, H 为磁场强度, J 为电流密度, E为电场强度, B为磁感应强度.对于恒定的电场:

式中电位满足泊松( Poisson)方程:

对于不存在电荷的空间部分有电荷体密度为零,上式退化为拉普拉斯( Laplace) 方程:

利用上述方程, 再加上边界条件, 利用Matlab 中的偏微分工具箱, 即可求解带电体周围空间的电场分布.输入pdetool可进入软件环境。

两点电荷的电场：两等值异号点电荷

单位, 两者间距为1,求其电势分布. 整个

求解域取中心为原点,半径为2 的圆,两空间电荷点位置为(-0.5,0)和( 0.5,0),作为一种近似,画一个尽量小的圆,取半径为0.05. 大圆的边界条件是Di richlet边界条件,取h= 1, r= 0,这种做法是模拟远处的电势为零. 由于大圆与小圆之间的区域没有电荷, 满足Laplace 方程, 因此在选择方程时选取Elliptic(椭圆)方程,其方程类