1. 引言
永磁体的空间磁场分布是其工程应用 [1] [2] [3] [4] [5] 的基础。目前,针对永磁体的磁场分布多采用有限元分析法和数值分析法,后者包括等效面电流法和等效磁荷法。二者均可辅助实现机械结构的设计和优化,但后者可实现机械结构的动力学建模仿真,以降低研制成本。故有必要对多边形永磁体的空间磁场进行数值分析。
矩形永磁体的使用较为广泛 [3] [4] [5],文 [6] 应用等效磁荷法解析了其磁场分布,文 [7] [8] 应用等效面电流法获得了其磁场分布的解析表达式。文 [9] 基于等效磁荷法研究了梯形永磁体的磁场分布,但其过于简化的推导过程不便于后续研究借鉴,解析式缺乏拓展性。文 [10] 指出,三角形、六边形、八边形等形状的永磁体可能有助于磁机器性能的提升,并基于等效磁荷法,解析了三角形永磁体的空间磁场,但表达式中的参数含义不明,缺乏拓展性。在三维微型数字电磁执行器 [11] [12] [13] [14] [15] 中,六边形永磁体被用作动子,实现了离散位置的控制。文 [12] 将六边形永磁体等效为矩形永磁体来解析其磁场,但因六边形和矩形永磁体结构差异较大,致使二者的磁场分布范围和数值均差异较大,该等效并不合理。磁制冷冰箱中使用八边形阵列的永磁体,提高了其冷却能力 [16];文 [17] 指出,八边形是磁粉成像技术中最合适的磁铁形状。
但目前还没有六边形永磁体 [15] 和八边形永磁体磁场的解析表达式,不便于实现其工程应用中的动力学建模。故本文以六边形永磁体和八边形永磁体为例,解析了多边形永磁体的磁场分布,并使用有限元分析结果对数值计算结果进行验证。基于解析表达式,分析了八边形永磁体较矩形永磁体的磁场分布优势。
2. 多边形永磁体磁场的解析
2.1. 解析原理
等效面电流法以安培分子环流假说和毕奥–萨伐尔定律为基础。永磁体的面磁化电流J和磁化矢量M之间的关系为
(1)
其中n为磁介质表面单位外法线方向,二者在数值上相等。永磁体的磁化强度M、面磁化电流大小J与其剩余磁通量密度Br的关系为
(2)
其中,
为真空磁导率,取值为
,单位为H/m。则薄层电流环dz0的电流强度I为
(3)
由毕奥–萨伐尔定律 [18] 的矢量表达式得到线电流元Idl在空间场点P(x, y, z)处产生的磁通量密度为
(4)
其中,r表示坐标原点到场点P(x, y, z)的矢径;
表示坐标原点到源点(x0, y0, z0)的矢径;
表示源点到场点的矢径;且
(5)
(6)
(4)式在闭合回路及永磁体厚度h上的积分,即为整个永磁体在场点P(x, y, z)处产生的磁通量密度。
2.2. 六边形永磁体
![](//html.hanspub.org/file/6-1251769x18_hanspub.png?20230119082715506)
Figure 1. 3D molecular circulation model of hexagonal permanent magnet
图1. 六边形永磁体的三维分子环流模型
![](//html.hanspub.org/file/6-1251769x19_hanspub.png?20230119082715506)
Figure 2. 2D analytical model of hexagonal permanent magnet
图2. 六边形永磁体的二维解析模型
六边形永磁体的三维分子环流模型和二维解析模型,分别如图1和图2所示。设图2中斜线的斜率为k,截距为b。x1、x2、x3及y1、y2、y3、y4分别为各顶点的坐标值。永磁体的边长为a mm,高为h mm,其磁化电流密度为J。根据2.1节介绍的永磁体磁场分布的解析原理,以图2中标号为3的载流导体l3在场点P(x, y, z)处产生的磁通量密度的法向分量
的计算为例,展开多边形永磁体的磁场分布的解析研究。载流导体3满足(7)式,则由(4)式在闭合回路l及永磁体厚度h上的积分,得到图1所示整个六边形永磁体在场点P(x, y, z)处产生的磁通量密度的法向分量
为(8)式,进一步化简得到(9)式。
(7)
(8)
(9)
式较复杂,按照常规方法难以计算,需先将其化简成文 [19] 中的标准积分形式再计算。在此,定义均以场点坐标分量x、y、z为自变量的函数符号F、F1、F2,具体为
(10)
(11)
(12)
其中,
表示整个
表达式中
分别取值h和0时的差值。
同理,可计算其余载流导体在场点P(x, y, z)处产生的磁通量密度。则由磁场叠加原理整理得到图1所示六边形永磁体的磁通量密度的解析表达式为
(13)
(14)
(15)
其中,
(16)
(17)
(18)
(19)
2.3. 八边形永磁体
图3、图4分别为八边形永磁体的三维分子环流模型和二维解析模型。永磁体的高度为h1 mm,边长为a1 mm,其磁化电流密度为J1。设图4中斜线的斜率为k1,截距为b1,xi (i = 4, 5, 6, 7)和yj (j = 5, 6, 7, 8)为其各坐标顶点。遵循上述解析过程,且定义如(21)~(23)式所示的系数,则图3所示的八边形永磁体的磁通量密度的由(24)~(26)式描述。
![](//html.hanspub.org/file/6-1251769x39_hanspub.png?20230119082715506)
Figure 3. 3D molecular circulation model of octagonal permanent magnet
图3. 八边形永磁体的三维分子环流模型
![](//html.hanspub.org/file/6-1251769x40_hanspub.png?20230119082715506)
Figure 4. 2D analytical model of octagonal permanent magnet
图4. 八边形永磁体的二维解析模型
(20)
(21)
(22)
(23)
(24)
(25)
(26)
其中,函数符号F、F1、F2同上,仅需将
中h变换为h1即可。
3. 有限元仿真验证
在MATLAB 2021a中建立了(15)和(26)式的数值分析模型,得到六边形永磁体和八边形永磁体外部空间磁场分布的数值分析结果。在有限元仿真分析软件ANSYS 2020R2中,分别建立了六边形永磁体和八边形永磁体的三维仿真模型,得到二者外部空间磁场分布的有限元分析结果。仿真参数设置见表1。通过对比有限元分析和数值分析的结果来验证(15)和(26)式的正确性。
在MATLAB 2021a中建立(15)和(26)式的数值分析模型时,首先需定义用于解析Bz分布的(11)和(12)式以及(15)和(26)式中所用参数。其次,每次循环都需计算每一条载流导体在场点P(x, y, z)处产生的Bz值,并进行加和。最后,经过多次循环,得到永磁体在一条线处所产生的Bz分布。两种方法的对比结果见图5、图6。
图5 表明,由(15)式得到的磁通量密度的法向分量Bz能很好地吻合有限元分析的结果。即本文所推导的六边形永磁体的磁场分布的解析表达式是正确的。图6表明,由(26)式得到的磁通量密度的法向分量Bz分布趋势与有限元分析的结果整体上是一致的。但在图形的顶部,有限元分析的结果存在毛刺,这是由于仿真过程中的数值噪声引起的。此外,由于有限元仿真中边界条件的选取和网格的划分等,导致曲线两端的幅值存在较大的误差,两端的波形并未完全重合。但这些偏差都是有限元分析中合理的存在,故本文所推导的六边形、八边形永磁体磁场的数值表达式是正确的。
![](//html.hanspub.org/file/6-1251769x51_hanspub.png?20230119082715506)
Figure 5. Magnetic field distribution of hexagonal permanent magnet
图5. 六边形永磁体的磁场分布
![](//html.hanspub.org/file/6-1251769x52_hanspub.png?20230119082715506)
Figure 6. Magnetic field distribution of octagonal permanent magnet
图6. 八边形永磁体的磁场分布
4. 八边形永磁体的三维磁场分布特性
为研究八边形永磁体的分布特性,对比分析了八边形永磁体和矩形永磁体的磁场分布特性。参数设置见表2。图7是二者在平面R2上的三维磁场对比结果,图8是八边形永磁体在线x = 2 mm,y = [−2.4142, 2.4142] mm上的磁场分布。
![](Images/Table_Tmp.jpg)
Table 2. Parameter settings of octagonal and rectangular permanent magnet
表2. 八边形和矩形永磁体的参数设置
图7中,八边形永磁体和矩形永磁体Bx、Bz的分布趋势一致;八边形永磁体By的幅值较小,且其波形与矩形永磁体的相差半个周期;二者Bz的峰值均是三个分量中最大的。图8也表明,在相同位置,八边形永磁体和矩形永磁体Bx的分布波形相似,八边形永磁体的分布波形具有更好的正弦性,而矩形永磁体的幅值较大;二者By的分布相差半个周期;二者Bz间的差值稳定。综上,虽然八边形永磁体的磁通量密度的幅值较小,但其波形的正弦性更好,且其Bz的幅值是三个分量中最大的,与Bx、By的差值较大,其磁场的综合分布特性较矩形永磁体的好。
![](//html.hanspub.org/file/6-1251769x55_hanspub.png?20230119082715506)
Figure 7. 3D magnetic field distribution of octagonal and rectangular permanent magnets
图7. 八边形和矩形永磁体的三维磁场分布
![](//html.hanspub.org/file/6-1251769x56_hanspub.png?20230119082715506)
Figure 8. Magnetic field distribution of octagonal and rectangular permanent magnets
图8. 八边形和矩形永磁体的磁场分布
5. 结论
有限元分析结果验证了本文全参数化的数值表达式的有效性,有助于多边形永磁体结构的动力学建模和仿真。磁场分布的对比结果表明,八边形永磁体较矩形永磁体具有更好的磁场分布特性。结合坐标变换,本文的研究结果可解析任意磁化方向和任意放置方式的棱柱形永磁体的磁场。此外,也可为今后的研究提供更多的永磁体形状选择和理论依据。
NOTES
*通讯作者。