间接边界单元法IBEM在高梯度磁选机磁系设计中的应用.pdf

间接边界单元法(IBEM)在高梯度磁选机 磁系设计中的应用 何桂春何平波 (南方冶金学院环境与建筑工程学院赣州341001) 摘要本文利用了边界单元法在求解三维和开域问题上的优点，将它应用于计算高梯度磁选机磁系分选腔内的磁场 分布，并且提出了计算古有铁磁介质边界的间接边界单元法。通过编程计算，结果表明，计算结果与实测值吻合得相当好，说 明该方法简单易行，精确度高，可供实际应用。 关键词高梯度磁选；磁系；设计；边界单元法 1前言 在高梯度磁选设备研制过程中，磁系设计是关 键，对磁系磁场的了解是十分重要的，它对分析磁 选设备的性能，进而确定适宜的磁系结构参数、磁 极形状和尺寸，以及研究磁性颗粒受力情况起着重 要的作用。同时，为了优化结构、节省材料，设计制 造供不同用途且能满足特定磁场要求的磁系，就必 须对其磁场的大小与分布，作出严密的分析与准确 的计算。因此磁系磁场的理论研究是磁选领域中 的一个重要课题。目前，高梯度磁选设备的磁系设 计计算还没有突破传统的计算范围，漏磁系数的确 定也不能摆脱对经验的依赖。在电子计算机技术 迅猛发展、电磁场数值计算日益深入的今天，完全 应该。也完全可能研究出更为简便、准确，建立在电 磁场理论研究基础上的磁系设计计算的新方法。 利用计算机，采用数值计算方法来解决电磁场 问题，是近三十年来发展起来的一门新学科。在磁 选领域，尽管很少有人利用数值分析法来进行计 算，但在高梯度磁选应用中曾有人利用有限差分法 和有限元法对分选腔内矩形和任意多边形磁介质 的磁场特性进行过研究“o J，中南工业大学矿物工程 系何平波”1也曾利用有限元法进行了磁系设计，效 果十分好。但是，他们都没有对高梯高磁选机磁系 磁极头开7L以后分选腔内磁场的分布进行研究，而 且有限差分法和有限元法从原理上而言，不适合于 处理开域问题。本文就是针对该问题应用间接边 界单元法去研究高梯度磁选设备磁系磁场分布情 况的。 边界单元法是本世纪七十年代发展起来的一 种新的数值计算方法，并且在近几十年得到了迅速 108 的发展．获得了广泛的应用。它是在积分定理的基 础上，结合有限元法的插值方法，对所考虑的区域的 边界进行离散化处理，然后将边界积分方程转化为 一多元代数方程，求解即可得到边界上的值。与有 限元法相比，其主要优点在于降维特性，即只需对区 域的边界进行离散，因而具有输人数据少，计算时间 短，对计算机存储容量要求低等特点，便于微机应 用，特别适合与处理三维问题和开域问题。本文正 是利用了边界单元法在求解三维和开域问题上的优 点，将它应用于计算高梯度磁选机磁系分选腔内的 磁场分布．并且提出了计算含有铁磁质边界(铁铠) 的间接边界单元法，从而解决了开域恒定磁场问题， 为研究高效高质的磁选社提供了一个具有广阔前途 的崭新方法。 2铠装螺线管线圈磁场计算的间接边界积分方程 的导出 当磁场中含有铁磁质(铁铠和磁极头)时，在电 流产生的外磁场作用下而被磁化，从而空间任意点 的磁场强度就是有传导电流和磁化了的铁磁质分别 单独作用时产生的磁场之和。当场中无传导电流分 布而仅有磁化了的铁吃斋单独作用时，有电磁场理 论可知，当传导电流为0时，我们可以采用标量磁位 “来描述磁场强度，因此有 v X巩=0 (1) 而比=一Vp (2) 将B=№(比+M)代人V·B=O，有 V·比=一V‘M (3) 进而有如下形式的准泊松方程： V·(VU)=V·M=一pⅢ (4) 几=一．uoV‘M (5) 式中u为标量磁位，‰为体磁荷密度，M为磁 化强度。如果铁磁质处处均匀(媒质特性不因点的 位置而变)，则M为常矢量，因此pⅢ恒为0，铁磁质 内无磁荷分布，而不需考虑体磁荷密度。由于分界 面的存在，铁磁质对空间磁场的贡献将通过磁荷密 度体现出来，且这时面磁荷密度为常量。设面磁荷 密度为％，则场域中任意点的磁位U“’： u2}葡毒葡曲+c (6) 式中R(P，Q)为源点Q与场点P之间的距离，s 为铁磁质边界面，c为积分常数，取决于磁位参考点 的选择，式(6)即为磁位计算的间接边界积分方程。 在实际的磁场计算中，通常是很难预先给出一 系列等磁位面的，所以式(6)的应用极为有限。但我 们关心的某～空间的磁场分布而不是等磁位分布， 故式(6)可以转化为： 凡=一可u2{葡翩凸 (7) 当磁场中同时存在传导电流和磁介质时，根据 迭加原理，任意点的磁场是由传导电流和磁化了的 磁介质分别单独作用时产生的磁场之和。即合成 磁场： H(P) 珏(P)+也(P)=坟(P)一Vu(P) (8) 式H是在无限大真空中仅由传导电流作用时 产生的磁场强度。 根据磁场的边界条件我们可以推导出由传导 电流和磁介质共同作用时的问接边界积分方程”1： 訾“P)．V{磊禺础hn(P)P‘}Ic (P) (9) 式中‰——铁