武大电气工程电磁场仿真 实验报告.docxVIP

武汉大学 工程电磁场及高电压综合实验 一、题目 有一极长的方形金属槽，边宽为1cm，除顶盖电位为100sinπx V外，其他三面的电位均为零，试用差分法求槽内电位的分布。 二、解题原理：均匀媒质中的有限差分法 我们在求解场的分布时，当边界形状比较复杂时，解析分析法不再适合了，我们可以采用数值计算的方法，数值计算法的基本思想，是将整体连续的场域划分为若干个细小区域，一般称之为网格或单元，如图1所示，然后用所求的网格交点（一般称为节点或离散点）的数值解，来代替整个场域的真实解。因而数值解，即是所求场域离散点的解。虽然数值解是一种近似解法，但当划分的网格或单元愈密时，离散点数目也愈多，近似解（数值解）也就愈逼近于真实值。 实解。在此处键入公式。  图1 场域的剖分，网格节点及步长 (一）、场域的剖分、网格节点及步长 由边界Γ所界定的二维平行平面场（见图1），若采用直角坐标系则可令该场处在xoy平面内。 所谓场域的剖分就是场域的离散化，即将场域剖分为若干个网格或单元。最常见最简单的剖分为正方形剖分，这种剖分就是在xy平面上作许多分别与x轴及y轴平行的直线，称为网格线。网格线的交点称为节点或离散点，场域内的节点称为内节点，场域边界上的节点称为边界节点。两相邻网格线间距离称为步长，一般以h表示。若步长相等则整个场域就被剖分为许多正方形网格，这就是正方形剖分。节点（离散点）的布局不一定采用正方形剖分，矩形剖分也常采用，正三角形剖分偶尔也被应用，不过最常见的最简单的仍然是正方形剖分。 （二）、差分与微分 从前面的分析可知，稳恒电、磁场的求解问题，归根到底是求解满足给定边界条件的偏微分方程（泊松方程或拉普拉斯方程）的解的问题所谓差分方法，就是用差商近似代替偏微商，或者说用差分代替微分，从而把偏微分方程转换为差分方程，后者实际上为代数方程。因此这种转化有利于方程的求解。 下面分别对一阶及二阶的差分公式进行推导。首先回顾有关偏导数的定义，有 (1) 因此当充分小时，可近似地用或代替，所谓差分公式，即是基于上述观点推得的。 设图1所示场域中的位函数为A，任取一网格节点0，它在xy平面上的坐标为（），记节点0的矢量磁位为，并把与节点0相邻的其他四个节点1、2、3、4的矢量磁位分别记为、、、，将节点0处函数A的一阶偏微商，用1、0两点函数值的差商近似代替，则有 （2） 式（2）中之差商，称为向前差商。上述一阶偏微商也可用0、3两点函数值的差商近似代替，称为向后差商，得 （3） 同理，对于偏微商也可分别用向前或向后差商近似代替，所得结果为 (4) 或 (5) 式（2）~式（5）就是一阶差分公式。 二阶差分公式可以在一阶差分的基础上进一步推出，0点处的二阶偏微商，如对一阶差商再取差商则得二阶差商为 用上式近似代替二阶偏微商就是要推演的二阶差分公式。 又由于 = （6） 所以 （7） 同理 （8） 式（7）及式（8）即为二阶差分公式。 需要说明的是，用差分近似代替偏微分，必定会产生误差。理论分析表明：其误差与步长成比例，因此若网格剖分得愈小，步长h就愈小，从而引起的误差也愈小。 (三)、均匀媒质中泊松与拉普拉斯方程的差分离散格式 设图2所示的平行平面场，场域每边长均为b,场域内电流密度为δ，媒质磁导率为，边界上的矢量磁位值已知，求域内矢量磁位。所提问题为第一类边值问题，则 （9） 应用差分法解图2所示场域的步骤如下： 图2 边界与网格线重合 第一步，若采用正方形网格剖分，即将场域剖分为具有四个内点（即点1,2,3,4），边界与网格线重合的九个离散单元（九个网格）。这就是把连续的场域进行离散化，从而将求解场域内矢量磁位函数的问题，转化为求1,2,3,4各内点的矢量磁位值的问题。 第二步，根据式（7）、式（8）的二阶差分公式，列出1,2,3,4各内点泊松方程的差分表达式。因为h=b/3，则各内点1,2,3,4的矢量磁位的差分方程为 内点1： （10） 内点2：