松弛迭代法求解非线性磁场的收敛性质的研究.pdfVIP

松弛迭代法求解非线性磁场的收敛性质研究’ 吴迅雷丁耀根 中国科学院电子学研究所 100080北京市中关村 引言 目前市场上具有磁场计算功能的大型软件已经可以见到十几种，例如，M瓶a’ 特别专业化，只作磁场计算的小软件又很难买到。对于很多需要作磁场计算的工程技术 人员来说，自己编个程序进行计算仍然是一个不错的选择。实际编程中可以借鉴各种书 刊中关于磁场计算的小程序或程序框图∽】。也可以从国际互联网络上下载商业软件的 演示版作参考。 不管采用有限差分法还是积分方程法，求解函数是矢势还是标势，最终都会得到～ 个线性方程组。用于求解线性方程组的方法，目前比较流行的有超松弛迭代法、共轭梯 度法等。当磁场中包含有非线性的媒质，如永磁材料AlNiCo，SmCo，软磁材料电铁4等 时，磁感应强度随磁场强度非线性地变化，这样的磁场叫做非线性的磁场。在非线性磁 场中，当磁饱和现象发生时，要非常小心 地保持计算的收敛性。 中科院电子所高功率微波器件研究室 开发有一个约2，700行名为CTA的Fortran 程序用于2D轴对称磁场的计算。这个程序 采用有限差分法和超松弛迭代法(SOR)。 它的计算中只能假设铁块中的磁导率是固 定不变的，如果顾及磁导率分布的变化， 则铁块中的磁感应强度会变得异常的高， 并且迭代将最终发散(见图1，图中纵坐标 的变量是柱坐标对称轴上的B，)。 图1迭代收敛时与发散时的磁感应强度 为改善CTA的收敛性，作者对这个程序做了一些修改，这些改进取得了良好的效果。 作者愿意把程序修改中获得的经验公布出来，希望那些想为非线性磁场问题找到简单解 决办法的科技工作者们能从此受益。 求解办法 ’本文工作包含于合同号臣家自然科学基金资助项目中． 如果用矩阵方程AX=b代表一个线性方程组，则CTA程序对它的求解过程可以图 示为 椭 gf zy’———◆．{ 咖 、l●卜●_、 础)二{b3一妒”， 图2磁场计算中线性方程组的求解 其中Ⅳ代表Z的第Ⅳ个迭代循环，每个循环包括一百次或一百次以上的迭代。 CTA中的工为磁通函数o．矩阵4与Ⅳ的分布有关，常数项b则由边界条件决定。 关于已知系数矩阵和常数项线性方程组的求解已有成熟的理论∞-91。磁场问题却与此不 尽相同。在以上求解过程中，Ⅳ和t，的实际分布最初是无法预知的，因此在每个迭代 循环前对彳和b必须给出恰当的初值，并且在下一个迭代循环前进行修正。CTA程序采 用超松弛迭代法求取Z的值，用矩阵表示为 Z㈣)=GX(‘)+d (1)． 其中G是迭代矩阵，k表示迭代轮次。SOR方法中G的具体形式为 (z)， G=O—coD-t￡)_1[(1一国p+国D‘1u】 其中』是单位矩阵， 厶U，D分别是A的上三角，下三角和对角阵，CO是松弛因子。 我们的工作主要是找出迭代发散的原因，确保迭代过程有很好的收敛性质。图2中可以 看到，磁场计算的收敛过程除了包括每个迭代循环中X的迭代，还包括其后对A和b的 修正。 1．超松弛法 对于一个系数矩阵爿和常数项b固定的线 性方程组，总有一个最佳松弛因子国。使得迭代 Ⅵ 过程最快地收敛。理论上说，它的数值接近2。 ： L 通常在实际操作中∞。显得略微大了点。CTA中 (-o取值为1．250)^一0．5．这对于线性磁场的计 算已足够，因为铁块中没有磁饱和现象时，铁块 唧八 中的“值可以近似取为固定值而不作修正。当