实验2matlab有限元计算

《实验2matlab有限元计算》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、MATLAB电磁场有限元计实验目的：1、了解有限元算法的原理，熟练运用MATLAB环境的PDE工具。2、熟练运用PDE工具分析简单的电磁场边值问题。实验内容：一、有限元简介在电磁场的计算屮，仅对那些具有最简单边界条件和场域几何形状规则的问题冰有解析解，多数问题的求解必须用数值计算的方法,其场域分布的数值计算内容是学习难点。本实验将有限元法和Matlab结合起来对电磁场教学中的电位分布问题进行计算。结果表明使用Matlab对有限元分析编程中的矩阵进行处有限元法是以变分原理和剖分插值为基础的一种数值计算方法，其基本思想是将场域方程等价为一个条件变分问题，然后由条件变分问题对场域进行剖分离散为方程组

2、进行求解。对于一个电场来说，其储能总是趋于最小，这样变分法的泛函和电场的储能就联系起来丫。对于边界为L的无源空气介质二维静电场屮,一个封闭场域S内的等价能量泛函可以写为：，陶=ieJ/r^2+（p2jds=mm50=乙Z0在有限元分析中，将所研究的区域S划分成有限的n个三角形网格单元。对应m个节点，ds为单元e的面积。对任意三角形单元e中任一点的电位可以认为由该三角形的三个节点（分别设为i、j、k）上的电位u随该点坐标x、y变化而线性确定。因此，对于单元e构造插值函数：Ue=OUzi+QiZj+OkUk=ChiUhh番ij9k其中ah称为形状函数。那么冇插值函数的一阶偏导数为:从而得到能量函数

3、We:dl¥2j?hAstTTA-nvS3Oh丁Hhi，j,kW,(u)+（S心/，o•，人•，则将单元e中的能量函数We对毎一个节点电位ill（l=i，j，k）求一阶偏导数，得:er（）uI1(S+h^,vlh)ds=()表示为矩阵形式有:'Am•••kInu=•••■A*n1•••鬱•••“kmu••••-=0•()Ac“/，kflat(ixrfen然后进行总体合成，将各单元的能量函数对同一节点的电位一阶偏导数相加，获得所要求解的线性方程组。由以上分析，可知在该场域内电场冇限元数学模型为:IKJ-[U]=0式中U为！1个节点处的待求电位，K为n阶矩阵。最后进行强加边界条件处理，消去已

4、知电位节点在系数矩阵中所在的行和列，得到简化后的方程，继而可以对电位进行求解。流程框图如下图所示。份解，结也二、静电场仿真▽•/）=P▽x//=J

5、异号点电荷±1单位，两者间距为1，求其电势分布.整个求解域取中心为原点，半径为2的圆，两空间电荷点位置为（-0.5,0）和（0.5,0），作为一种近似，画一个尽量小的圆，取半径为0.05.大圆的边界条件是Dirichlet边界条件，取h=1，r=0,这种做法是模拟远处的电势为零.由于大圆与小圆之间的区域没有电荷，满足Laplace方程，因此在选择方程时选取Elliptic（椭圆）方程，其方程类型为:—V•(cVI)+au-f.取系数为1，0,0.在表示点电荷的小圆内，我们认为电荷是均匀分布的，满足Poisson方程，在选择方程时也取Elliptic方程，取系数为c=1，a=0,f=0.2.其两

6、点电荷电势分布上图所示，电力线用箭头表示.三、静电场中的导体问题描述：在电场强度为E的静电场中放置一根无限长的导体，研宄截面上的电势分布。首先画一个2*2的矩形R1,然后在中心原点画半径为0.3的圆E1.然后将Setformula对话框屮的公式改为R1-E1，表示求解区域为二者之差.矩形所有的边界条件是Dirichlet边界条件，取h=l，r=y.而在圆的边界取h=l,r=0.由于求解域没有电荷，因此在选择方程时选取Elliptic(椭圆)方程，系数为c=l,a=0,f=0.其电势分布如下图所示，电力线用箭头表示.四、两根载流长直导线的磁场问题描述：两根载流长直导线，相距为0.8,导线直径为0

7、.2,求电流引起的磁场.从麦克斯韦(Maxwell)方程组出发，其磁场强度B和磁感应强度H的关系为：B=fiH磁场势A与B有如下关系：B=VxA故可简化为椭圆方程:▽2A=-"J画出大小为2*2的矩形R1，两导线用直径为0.2、相距0.8的两个圆表示.矩形的边界条件是Dirichlet边界条件，取h=l,r=Oo这种做法是模拟远处的磁场势为零.在设置方程类型时，选取应用模式为Mangetostat