实验2 MATLAB有限元计算

MATLAB电磁场有限元计算实验目的：1、了解有限元算法的原理，熟练运用MATLAB环境的PDE工具。2、熟练运用PDE工具分析简单的电磁场边值问题。实验内容：一、有限元简介对场域进行剖、，这样变分法的node（节点）在电磁场的计算中，仅对那些具有最简单边界条件和场域几何形状规则的问题才有解析解，多数问题的求解必须用数值计算的方法，其场域分布的数值计算内容是学习难点。本实验将有限元法和Matlab结合起来对电磁场教学中的电位分布问题进行计算。结果表明使用Matlab对有限元分析编程中的矩阵进行处理，对场域进行剖、，这样变分法的node（节点）方法，其基本思想,,element（单angularElement）有限元法是以变分原理和是将场域方程等价为一个条＜散为方程组进行求解。对于-函和电场的储能就联系起来了。对于边界为L方法，其基本思想,,element（单angularElement）在有限元分析中，将所研究的区域S划分成有限的n个三角形网格单元。对应m个节点,ds为单元e的面积。对任意三角形单元e中任一点的电位可以认为由该三角形的三个节点（分别设为i、j、k）上的电位u随该点坐标x、y变化而线性确定。因此，对于单元e构造插值函数：CtZZi+CtZZi+OkZik=/J01!Zihv，A'其中ah称为形状函数。那么有插值函数的一阶偏导数为:从而得到能量函数We:A[(V迪mJ迪+『】jf=。捉*«dxJdx'W顷A[(V迪mJ迪+『】jf=。捉*«dxJdx'W顷为dW^u)_duik11*.■k1n.0::,■••0-jfcnl*1■kntj-ZZ»--QI［澜I“,,表示为矩阵形式有：然后进行总体合成，将各单元的能量函数对同一节点的电位一阶偏导数相加，获得所要求解的线性方程组。由以上分析,可知在该场域内电场有限元数学模型为:式中U为n个节点处的待求电位,K为n阶矩阵。最后进行强加边界条件处理,消去已知电位节点在系数矩阵中所在的行和列，得到简化后的方程，继而可以对电位进行求解。流程框图如下图所示。、静电场仿真VZ）=p<VXH=JIvxE二0、V=0静态场满足上方基本方程，式中D为电位移，A为电荷密度，H为磁场强度，J为电流密度,E为电场强度,B为磁感应强度.对于恒定的电场：E--M式中电位满足泊松（Poisson）方程：对于不存在电荷的空间部分有电荷体密度为零，上式退化为拉普拉斯（Laplace）方程：v2^=a利用上述方程，再加上边界条件，利用Matlab中的偏微分工具箱，即可求解带电体周围空间的电场分布.输入pdetool可进入软件环境。两点电荷的电场：两等值异号点电荷±1单位，两者间距为1，求其电势分布.整个求解域取中心为原点,半径为2的圆,两空间电荷点位置为（-0.5,0）和（0.5,0），作为一种近似，画一个尽量小的圆，取半径为0.05.大圆的边界条件是Dirichlet边界条件，取h=1,r=0,这种做法是模拟远处的电势为零.由于大圆与小圆之间的区域没有电荷，满足Laplace方程，因此在选择方程时选取Elliptic（椭圆）方程，其方程类型为:取系数为c=1,a=0,f=0.在表示点电荷的小圆内，我们认为电荷是均匀分布的,满足Poisson方程，在选择方程时也取Elliptic方程，取系数为c=1,a=0,f=0.2.其两点电荷电势分布上图所示,电力线用箭头表示.三、静电场中的导体问题描述：在电场强度为E的静电场中放置一根无限长的导体,研究截面上的电势分布。首先画一个2*2的矩形2,然后在中心原点画半径为0.3的圆E1.然后将Setformula对话框中的公式改为R1-E1，表示求解区域为二者之差.矩形所有的边界条件是Dirichlet边界条件，取h=1,r=y.而在圆的边界取h=1,r=0.由于求解域没有电荷，因此在选择方程时选取Elliptic(椭圆)方程,系数为c=1,a=0,f=0.其电势分布如下图所示，电力线用箭头表示.四、两根载流长直导线的磁场问题描述：两根载流长直导线,相距为0.8，导线直径为0.2,求电流引起的磁场.从麦克斯韦(Maxwell)方程组出发，其磁场强度B和磁感应强度H的关系为：B=hH磁场势A与B有如下关系：B=VA故可简化为椭圆方程:V2A=-|iJ画出大小为2*2的矩形R1,两导线用直径为0.2、相距0.8的两个圆表示.矩形的边界条件是Dirichlet边界条件，取h=1,r=0。这种做法是模拟远处的磁场势为零.在设置方程类型时，选取应用模式为Mangetostatics.故在选择方程时选取Elliptic(椭圆)方程，对于矩形其它部分系数取R=1、J=0.在表示导线的圆内，取R=1,J=1.两根载流长直导线的磁场势和磁力线如下图所示,磁力线用箭头表示.实验要求：有限元法计算电磁场问题的基本思想是什么？求解静电场两个等量异号带电导线电场位函数及电场强度。求解静电场中的导