章计算电磁学.ppt

* * 第4章 积分方程法 4.1 引言 积分方程法的发展： 1）是由于微分方程法对某些问题的解，需要计算机内存很多，计算精度难于提高，如对三维场的计算，节点剧增，需要寻求更适当方法计算； 2）在此同时，一些从事解析方法和均匀磁化法磁场计算的研究人员从微分方程法的离散概念中吸取有用的东西，将源区进行离散，发展成积分方程法，代表程序有GFUN（可计算二维、三维场）、BIM等程序。 第4章 积分方程法 积分方程法泊松方程的积分形式是由泊松方程的微分形式导出来的，又可由积分形式的Maxwell方程导出，可以证明两者是等价的。 本章只做概念性介绍。 积分方程法的应用： 求解区域内介质数目少，交界面形状不特别复杂的情况。 不需要对全部求解区域划分网格，只需对有源区和感应(或极化)产生的“第二次源”区划分网格，节点数大大减少，系数矩阵阶数比有限元(或差分法)的系数阵阶数低得多，在求解时计算时间短、费用低。 4.2 积分方程法的物理概念和基本公式 积分方程法从宏观的角度来描述场，场区中每点的值仅取决于所有场源对它的影响。场点和源点的联系是通过毕奥－萨伐定律实现。由于离散只在源区进行，加上恒定磁场问题的电流分布是已知的，因此实际上离散只在非线性铁区内进行，这使得数据输入和网格剖分大为简化。 第4章 积分方程法 磁场积分方程法主要采取的方法： （1）磁化值直接解法； （2）标量磁位的边界积分法 —— 边界元法； （3）等值面电荷积分方程法。 第4章 积分方程法 磁化值直接解法： （1）积分方程法的物理概念： 积分方程法解磁场问题的基本思想是认为在空间任一点的场是由电流源所产生的磁场和介质被磁化所产生的磁场迭加而成： 下标为 m 的量表示由磁介质磁化所产生的磁场。c 表示电流源存在所产生的磁场。 介质在空间各点所产生的场强的大小不仅取决于它和场点间的几何距离，而且与介质的磁化强度有关，磁化强度除了与本身材料有关以外，还与周围电流源及介质产生的场有关，相互影响，最后趋于稳定，其关系是一个较复杂的非线性关系。 第4章 积分方程法 （2）基本公式： 积分方程法求场的基本思想和物理概念，关键问题是求出产生磁场的源在所求场点的场强表达式。 a）电流源 c 在空间 r 处所产生的磁场 根据毕奥－萨伐定律，对二维平面对称场和轴对称场，讲义中已经给出了推导。因此导线电流在 P 点产生的场量可认为是已知量。 第4章 积分方程法 毕奥－萨伐定律如何得出？ 泊松方程 的特解为： r 为场点坐标，r’ 为源点坐标；场点和源点之间的矩离 R = r - r′，其方向由源点指向场点。 b）磁介质在场点所产生的磁场 第4章 积分方程法 此公式利用矢量磁位与磁矩概念可以导出，推导过程略。 利用此公式可以求出磁介质在空间所产生的磁场分布，但须先知道磁化强度 M。实际上一般 M 是未知的，因此只有采用先假设的方法，把磁介质剖分为许多单元。 I. 假设磁化体单元 a 的磁化强度为 Ma，并认为是一个常数，同样 b 单元中磁化强度为 Mb，也为常数。 II. 取源点为 b，场点为 a，在上式中，单元内 M 看成常数，可以从积分中提出。剩下的积分参数只与单元几何形状、场点位置有关，在确定位置情况下，积分为常数，用 Cab 表示： 第4章 积分方程法 III. 对所有的场源写出联立方程组： 又 （在各向同性媒质中） 代入上式可得 整理可得 得到一个具有三个方向分量的联立方程组，系数 C 是3×3阶张量。 第4章 积分方程法 若磁化率 和 是已知的，可以写成六个线性方程，有六个未知数：Hax，Hay，Haz，Hbx，Hby，Hbz可求出。但实际上 和 不是常数，要经过反复迭代进行求解。 IV. 计算框图： 否 参数输入，假定 值 形成系数 计算电流场量 是 解方程求 H 查 B ～H 曲线求 计算M，H 输出 第4章 积分方程法 积分方程法用于求解平行平面场、轴向场情况可以自学，加深对积分方程法的理解。 （3）积分方程法的特点： 1) 离散区域：只需对电流源和磁性介质进行离散，并且是分别进行考虑，因此可使电流源和磁性介质对称性充分得到利用； 2) 不需要专门对边界条件进行讨论：在计算中考虑到了空间所有场源的作用，与外界无能量交换，可以不考虑边界条件； 3) 由积分方程法得到的方程组其系数矩阵是满矩阵，离散单元数受到限制，故此法对有限小的无源区的封闭边界问题或饱和情况较复杂的铁区