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

实验2matlab有限元计算.doc

MATLAB电磁场有限元计算实验目的1、了解有限元算法的原理，熟练运用MATLAB环境的PDE工具。2、熟练运用PDE工具分析简单的电磁场边值问题。实验内容一、有限元简介在电磁场的计算中,仅对那些具有最简单边界条件和场域几何形状规则的问题才有解析解,多数问题的求解必须用数值计算的方法,其场域分布的数值计算内容是学习难点。本实验将有限元法和MATLAB结合起来对电磁场教学中的电位分布问题进行计算。结果表明使用MATLAB对有限元分析编程中的矩阵进行处理,程序设计清晰简便,易于理解和实现。NODE节点ELMENT单元节点单元有限元法是以变分原理和剖分插值为基础的一种数值计算方法,其基本思想是将场域方程等价为一个条件变分问题,然后由条件变分问题对场域进行剖分离散为方程组进行求解。对于一个电场来说,其储能总是趋于最小,这样变分法的泛函和电场的储能就联系起来了。对于边界为L的无源空气介质二维静电场中,一个封闭场域S内的等价能量泛函可以写为在有限元分析中,将所研究的区域S划分成有限的N个三角形网格单元。对应M个节点,DS为单元E的面积。对任意三角形单元E中任一点的电位可以认为由该三角形的三个节点分别设为I、J、K上的电位U随该点坐标X、Y变化而线性确定。因此,对于单元E构造插值函数其中AH称为形状函数。那么有插值函数的一阶偏导数为从而得到能量函数WE则将单元E中的能量函数WE对每一个节点电位ULLI,J,K求一阶偏导数,得表示为矩阵形式有然后进行总体合成,将各单元的能量函数对同一节点的电位一阶偏导数相加,获得所要求解的线性方程组。由以上分析,可知在该场域内电场有限元数学模型为式中U为N个节点处的待求电位,K为N阶矩阵。最后进行强加边界条件处理,消去已知电位节点在系数矩阵中所在的行和列,得到简化后的方程,继而可以对电位进行求解。流程框图如下图所示。二、静电场仿真静态场满足上方基本方程，式中D为电位移,为电荷密度,H为磁场强度,J为电流密度,E为电场强度,B为磁感应强度对于恒定的电场式中电位满足泊松POISSON方程对于不存在电荷的空间部分有电荷体密度为零,上式退化为拉普拉斯LAPLACE方程利用上述方程,再加上边界条件,利用MATLAB中的偏微分工具箱,即可求解带电体周围空间的电场分布输入PDETOOL可进入软件环境。两点电荷的电场两等值异号点电荷单位,两者间距为1,求其电势分布整个求解域取中心为原点,半径为2的圆,两空间电荷点位置为05,0和05,0,作为一种近似,画一个尽量小的圆,取半径为005大圆的边界条件是DIRICHLET边界条件,取H1,R0,这种做法是模拟远处的电势为零由于大圆与小圆之间的区域没有电荷,满足LAPLACE方程,因此在选择方程时选取ELLIPTIC椭圆方程,其方程类型为取系数为C1,A0,F0在表示点电荷的小圆内,我们认为电荷是均匀分布的,满足POISSON方程,在选择方程时也取ELLIPTIC方程,取系数为C1,A0,F02其两点电荷电势分布上图所示,电力线用箭头表示三、静电场中的导体问题描述在电场强度为E的静电场中放置一根无限长的导体,研究截面上的电势分布。首先画一个22的矩形R1,然后在中心原点画半径为03的圆E1然后将SETULA对话框中的公式改为R1E1,表示求解区域为二者之差矩形所有的边界条件是DIRICHLET边界条件,取H1,RY而在圆的边界取H1,R0由于求解域没有电荷,因此在选择方程时选取ELLIPTIC椭圆方程,系数为C1,A0,F0其电势分布如下图所示,电力线用箭头表示四、两根载流长直导线的磁场问题描述两根载流长直导线,相距为08,导线直径为02,求电流引起的磁场从麦克斯韦MAXWELL方程组出发,其磁场强度B和磁感应强度H的关系为BH磁场势A与B有如下关系A故可简化为椭圆方程2AJ画出大小为22的矩形R1,两导线用直径为02、相距08的两个圆表示矩形的边界条件是DIRICHLET边界条件,取H1,R0。这种做法是模拟远处的磁场势为零在设置方程类型时,选取应用模式为MANGETOSTATICS故在选择方程时选取ELLIPTIC椭圆方程,对于矩形其它部分系数取1、J0在表示导线的圆内,取1,J1两根载流长直导线的磁场势和磁力线如下图所示,磁力线用箭头表示实验要求1有限元法计算电磁场问题的基本思想是什么2求解静电场两个等量异号带电导线电场位函数及电场强度。3求解静电场中的导体等等势面，写出导体与电介质分界面上的衔接条件，并与实验结果对比。4写出方程的推导过程。求解两根载流长直导线的磁场磁密分布。2AJ5以上所有内容在做完实验后将截图及资料保存在WORD文档里，打印后上交、张献