2021年3月电工技术学报Vol.36 No. 6 第36卷第6期TRANSACTIONS OF CHINA ELECTROTECHNICAL SOCIETY Mar. 2021
DOI: 10.19595/s.200029
基于精确磁路的新型混合型轴向-径向
磁悬浮轴承研究
禹春敏1邓智泉1梅磊2庞古才3
（1. 南京航空航天大学自动化学院南京 211106
2. 南京工业大学电气工程与控制科学学院南京 211816
3. 南京磁之汇电机有限公司南京 211106）
摘要该文基于一种实心同极性轴向-径向磁悬浮轴承结构，其具有磁场分布均匀、轴向空间较小的优点。针对该新型混合型轴向-径向磁悬浮轴承结构紧凑、磁场分布复杂等特点，结合有限元二维仿真结果，建立考虑漏磁的精确磁路模型。根据磁路模型，计算该结构的气隙磁通密度、刚度和承载力，并与有限元
仿真进行对比，仿真与计算结果基本一致。实验结果表明转子能够在较高转速下实现稳定悬浮，验证了该结构的可行性。
关键词：混合型轴向-径向磁悬浮轴承等效磁路漏磁有限元
中图分类号：TH133.3
Rearch of New Hybrid Axial-Radial Magnetic Bearing Bad on
Accurate Magnetic Circuit
Yu Chunmin1 Deng Zhiquan1Mei Lei2 Pang Gucai3
（1. College of Automation Engineering Nanjing University of Aeronautics and Astronautics
Nanjing 211106 China
2. College of Electrical Engineering and Control Science Nanjing Tech University
Nanjing 211816 China
3. Nangjing Inomag Motor Co. Ltd Nanjing 211106 China）
Abstract This paper propos a structure of solid homo-polar axial-radial magnetic bearing, which has the advantages of uniform magnetic field distribution and small axial space. Considering the characteristics of this new hybrid axial-radial magnetic bearing, such as compact structure and complex magnetic field distribution, an accurate magnetic circuit model considering magnetic flux leakage was established combined with the two-dimensional finite element simulation results. According to the magnetic circuit model, the magnetic density in air, stiffness and bearing capacity of the structure were calculated. The simulation results are basically consistent with the calculation results. It is shown that the rotor can achieve stable suspension at high speed, which verifies the feasibility of the structure.
Keywords：Hybrid, axial-radial magnetic bearing, equivalent magnetic circuit, leakage magnetic field, finite element
咪咪影
国家自然科学基金资助项目（51577087）。
收稿日期 2020-01-07 改稿日期 2020-02-13
1220 电工技术学报 2021年3月
0引言
主动型磁悬浮轴承借助电磁铁、反馈回路、传感器、功率放大器等单元实现对转子位移的主动控制，被动型磁悬浮轴承利用永磁体与永磁体之间、永磁体与软磁材料之间的斥力或吸力来产生不可控的悬浮力[1]。混合型磁悬浮轴承结合了两者的特点，利用永磁体提供偏置磁场，电磁线圈产生控制磁场，显著减小了磁轴承的功耗，成为近年来磁悬浮轴承研究热点之一[2-7]。
混合型磁悬浮轴承又分为同极性结构和异极性结构[8-13]，同极性结构具有较低的磁滞损耗，被广泛用于飞轮储能、高速涡轮机机械等场合[14-15]。用于支承高速旋转体时，传统的磁悬浮电机采用两个径向磁悬浮轴承和一个轴向磁悬浮轴承[16-17]，会带来体积增大、临界转速降低及控制的复杂性提高等弊端。为了提高磁悬浮轴承的集成度，相关研究人员提出了轴向-径向磁悬浮轴承结构，永磁体同时提供轴向、径向偏置磁场，可以同时控制轴向单自由度和径向两自由度。文献[18]中并列式的轴向径向磁悬浮轴承，轴向气隙与永磁体不对称，存在偏置磁场不均匀、轴向空间大等问题。文献[19]提出一种径向充磁的永磁偏置轴向径向磁悬浮轴承，当需要提供较大的偏置磁通时，需要增加永磁体的轴向长度，从而导致磁轴承的轴向长度增加。
针对目前轴向-径向混合型磁悬浮轴承存在的问题，本文选用一种新型同极性轴向-径向混合磁悬浮轴承，永磁体采用轴向充磁，使磁悬浮轴承具有较小的轴向尺寸，轴向气隙与永磁体呈对称关系，不存在偏置磁场不均匀的问题。本文介绍了新型磁悬浮轴承的结构及工作原理。通过对磁场分布的分析，建立了考虑漏磁的精确等效磁路模型，基于磁路模型，分析了该模型的气隙磁通密度、刚度和承载力，
并将有限元仿真结果与磁路法计算结果进行了对比分析。基于设计参数制作了实验样机，实现了24 000r/min的稳定悬浮。
1轴向-径向磁悬浮轴承结构与工作原理
1.1轴向-径向磁悬浮轴承基本结构
混合型轴向-径向永磁偏置磁轴承的结构示意图如图1所示，定、转子铁心均由实心软磁材料制成，主要由两个定子圆盘、套筒、永磁体、转子、径向磁极、径向磁轭、轴向控制线圈、径向控制线圈组成。两个对称定子圆盘和套筒的内腔内设置轴向控制线圈，轴向定子圆盘在轴向上与转子留有轴向气隙；径向定子的内腔壁向轴心外延伸出圆周方向均匀分布的四个径向磁极，四个相同的径向磁极上分别绕有相同的径向控制线圈且在径向上与转子的外壁之间均留有径向气隙；径向磁轭设置在径向磁极的外壁；在径向磁轭的两侧分别设置轴向充磁的环形永磁体环，永磁体环的N极均朝向定子圆盘，S极均朝向径向磁轭。
图1 轴向-径向永磁偏置磁轴承结构示意图
Fig.1 Structural schematic diagram of
axial-radial permanent magnet biad magnetic bearings 为了便于分析，本文使用的结构参数如图2所示。
图2 轴向-径向永磁偏置磁轴承结构参数Fig.2 Structure parameters of axial-radial permanent
magnet biad magnetic bearings
1.2轴向-径向磁悬浮轴承的工作原理
收费英语
偏置磁场由两个轴向充磁的环形永磁体建立，通过永磁体、定子圆盘、轴向气隙、转子、径向气隙、径向磁极、径向磁轭形成磁回路，如图1中的实线部分所示。
以垂直方向为例，当转子受到一个向下的冲击力使其偏离中心位置向下运动，导致上半部分的径向气隙增大，偏置磁场磁通密度减小；下半部分的径向气隙减小，偏置磁场磁通密度增大。此时转子受到的合力方向也向下，使得转子偏离径向中心位置。为了使转子回到中心位置，需要在径向线圈中通入一定的电流，产生的径向控制磁场如图1经过径向磁极的虚线部分所示，在上半部分的气隙中，
第36卷第6期禹春敏等基于精确磁路的新型混合型轴向-径向磁悬浮轴承研究 1221
偏置磁场与径向控制磁场的方向一致，合成磁通密度增加；在下半部分的气隙中，偏置磁场与径向控制磁场的方向相反，合成磁通密度减小，转子受到的合力向上，将其拉向中心位置。因此，通过调节径向控制线圈中电流的大小和方向，能够产生所需的磁场力，将转子稳定悬浮在径向中心位置。
轴向与径向水平方向类似，当转子受到一个向右的冲击力使其偏离中心位置向右运动，左半部分的轴向气隙增大，偏置磁场磁通密度减小；右半部分的轴向气隙减小，偏置磁场磁通密度增大。此时转子受到的合力方向也向右，使得转子偏离轴向中心位置。为了使转子回到平衡位置，需要在轴向线圈中通入一定的电流，产生的轴向控制磁场如图1经过定子圆盘的虚线部分，在左半部分的气隙中，偏置磁场与轴向控制磁场的方向一致，合成磁通密度增加；在右半部分的气隙中，偏置磁场与轴向控制磁
场的方向相反，合成磁通密度减小，转子受到的合力向左，将其拉向中心位置。因此，通过调节轴向控制线圈中电流的大小和方向，能够产生所需的磁场力，将转子稳定悬浮在轴向中心位置。
综上所述，通过调节径向控制电流和轴向控制电流，使转子产生径向和轴向上的力，使转子稳定悬浮在中心位置。
2轴向-径向磁悬浮轴承等效磁路模型
2.1 考虑漏磁后的磁路分布
等效磁路法可以清晰地考虑磁轴承的磁场分布，故在磁轴承本体设计中被广泛使用。根据磁路分布，混合型轴向-径向永磁偏置磁悬浮轴承可以从三个部分考虑磁通分布：①偏置磁通；②轴向控制磁通；③径向控制磁通。
根据仿真结果，偏置磁通经永磁体、定子圆盘、轴向气隙、转子、径向气隙、径向磁极、径向磁轭建立偏置磁场，如图3a所示。轴向控制磁通经套筒、定子圆盘、轴向气隙、转子、轴向气隙、定子圆盘建立轴向控制磁场，如图3b所示。径向控制磁通经径向磁轭、径向磁极、径向气隙、转子建立径向控制磁场，如图3c所示。图3a中，偏置磁场漏磁路径LF pm1～LF pm22对应的漏磁磁阻为R pm1～R pm22；图3b中，轴向控制磁场漏磁路径LF z c1～LF z c19对应的漏磁磁阻为R l z c1～R l z c19；图3c中，径向控制磁场漏磁路径LF r c对应的漏磁磁阻为R l r c。
磁场路径几何分布如图4所示。根据图1偏置磁场磁路的分布，将模型铁心材料分成9个部分。
（a）偏置磁通（b）轴向控制磁通
（c）径向控制磁通
图3 考虑漏磁后的磁通分布
Fig.3 Flux distribution considering magnetic leakage
如图4a所示，5代表径向磁极，3、7代表定子圆盘，4、6代表转子，1、9代表径向磁轭，2、8代表永磁体和定子圆盘之间的气隙。
（a）偏置磁场路径
（b）轴向控制磁场路径（c）径向控制磁场路径
图4 磁场路径几何分布
Fig.4 Division of the geometry of magnetic field path 根据图1轴向控制磁场磁路的分布，将模型分成6个部分。如图4b所示，1、3代表轴向气隙，2代表转子，4、6代表定子圆盘，5代表套筒。
根据图1径向控制磁场磁路的分布，将模型分成3个部分。如图4c所示，1代表径向磁轭，2代表径向磁极，3代表转子。
2.2 等效磁路模型
根据以上分析及图4的磁路几何分布，建立考
1222
电 工 技 术 学 报 2021年3月
虑漏磁的精确磁路模型，如图5所示。
（a ）偏置磁场磁路模型
（b ）轴向控制磁场磁路模型
（c ）径向控制磁场磁路模型
图5 考虑漏磁影响的磁路模型
Fig.5 Magnetic circuit model considering leakage effects
图5a 中，F pm 为永磁体的矫顽磁动势，R pm 为永磁体的磁阻，有
pm c pm F H L = （1）
pm
pm m m
L R S μ=
（2）
式中，L pm 为永磁体轴向厚度；H c 为永磁体的矫顽力；μm 为永磁体磁导率；S m 为永磁体中性面面积。
R r g x+、R r g x −、R r g y+、R r g y −分别为径向四个气隙的磁阻，有
g 0r r x r
g x
R A μ+
+= （3） g 0r r x r
g x
R A μ−−=
（4） g +0r r y r
g y
R A μ+=
嗓子疼吃什么药（5） g 0r r y r
g y
R A μ−−=
（6） 式中，g r 为径向定转子间的气隙长度；x 、y 分别为
转子径向x 、y 方向的偏移量；μ0为直空磁导率；A r 为径向气隙的面积；R r g x +、R r g y +分别为x 、y 方向径向气隙增加的一侧气隙磁阻；R r g x −、R r g y −分别为x 、y 方向径向气隙减小的一侧气隙磁阻。
R z g +、R z g −分别为轴向两个气隙的磁阻，有
g 0z z z
g z
R A μ++=
（7） g 0z z z
g z
R A μ−−=
（8） 式中，g z 为轴向定转子间的气隙长度；z 为转子轴向的偏移量；A z 为轴向气隙的面积。
为了计算偏置磁通，定义
()
1lpm1lpm2lpm3lpm4lpm5lpm6=//////R R R R R R R ++（9） ()
2lpm22lpm19lpm20lpm21lpm18lpm17=//////R R R R R R R ++
（10）
3lpm7lpm8lpm9lpm10=//////R R R R R （11） 4lpm13lpm14lpm15lpm16=//////R R R R R （12）
5lpm11lpm12=//R R R （13） g g g g g =//////r r x r x r y r y R R R R R +−+− （14）
式中，R i1、R i3、R i4、R i5、R i6、R i7、R i9为铁心材料的磁阻；R i2、R i8为永磁体与圆盘之间的气隙磁阻；F 1、F 2、F 3为所在点的磁动势，根据节点电压法，计算出偏置磁场。
当转子悬浮在平衡位置时，轴向偏置磁通和径向偏置磁通分别为
()23
pm 3g 3g i3i43
g z z z z F R R R
R R R R R R
φ−
−−
=
⎛⎞
+++⎜⎟⎜⎟+⎝
⎠
（15）
(
)(
)
15g g g pm 5i55g i5g
圆通老板
g g g g g g g g g g g g r x r y r y r r r r x r x r y r x r x r y r x r y r y r x r y r y F R R R R R R R R R R R R R R R R R R R R R R φ−+−
+−++−−+−+−+−
−=
+++++ （16）
第36卷第6期
禹春敏等 基于精确磁路的新型混合型轴向-径向磁悬浮轴承研究 1223
图5b 中，N z 为轴向线圈匝数，i z 为轴向控制电流。R z e2、R z e4、R z e5、R z e6为轴向控制磁路经过的铁心磁阻，R z e1、R z e3为轴向气隙磁阻，在平衡位置，R z g+与R z g −两者相等，有
e1g z z R R −= （17）
e3g z z R R += （18）
6l c1l c2l c3////z z z R R R R = （19）
7l c14l c15l c16////z z z R R R R = （20）
()()8l c1967l c4l c18l c13l c17////z z z z z R R R R R R R R =++++
（21）
9l c5l c6l c7l c8//////z z z z R R R R R = （22） 10l c9l c10l c11l c12//////z z z z R R R R R = （23）
为了计算整个回路的总磁阻，引入系数o z c1、o z c2、o z c3，有
8
c18t1
z z R o R R =+ （24）
9
c29e1
z z R o R R =
+ （25）
10
c3
10e3
z z R o R R =+ （26） t1e2e4e6e1c2e3c3z z z z z z z z R R R R R o R o =++++ （27）
整个轴向控制磁场回路的磁阻为
t e5t1c1z z z z R R R o =+ （28）
轴向控制磁场磁通为
c1c2i t
z z z z z z N i o o
R φ= （29）
图5c 中，N r 为径向线圈匝数，i x 、i y 分别为径向x 、y 方向上的控制电流。R r e1、R r e2、R r e3分别为径向磁轭、径向磁极、转子的磁阻，有
pmo pmi
e10s p
ln 2πr r r r R T μμ=
（30） pmi pi
p e20s p
ln
22ππ
r r r r r r R T αμμ=
⋅ （31）
ro
e30r 2r r r
r R A μμ=
（32）
式中，r pmo 、r pmi 分别为永磁体的外径和内径，数值上也等于径向磁轭的外径和内径；r r pi 为径向磁极的内径；T r p 为径向磁极的厚度；r r ro 为转子的外径；
μs 为定子铁心的相对磁导率；μr 为转子铁心的相对
磁导率；αr p 为径向磁极的极弧度。
由于径向x 和y 是对称结构，以径向y 方向为例，仅有y 方向上的线圈通电时，该径向控制磁通回路的总磁阻为
e1l c
t e2g g e3e1l c
2r r r r r y r y r r r R R R R R R R R R +−=++++
+（33）
径向控制磁场磁通为
i t 2r y
r r N i R φ= （34）
轴向承载力和径向承载力分别为
pm i 02z z z z
F A φφμ= （35）
pm i
brooding02r r r r
F A φφμ=
（36）
径向电流刚度和径向位移刚度分别为
()
15
i 0t 5i55g i5g +r r r r r r F N R k A R R R R R R R μ−=
+ （37） ()朗朗乾坤是什么意思
222
105g
s 2
5i55g i5g 05i5g
222r r r r r r r r r r r A F R R k g R R g R R g R R g A R R R μμ=
+++
（38）
轴向电流刚度和轴向位移刚度分别为
(
)
23c1c2
i 3g 0t
3g i3i43g 2z z z z z
z z z z F N R o o k R R A R R R R R R R μ=
⎛⎞+++
⎜⎟⎜⎟+⎝
⎠
（39）
2222203
s 2
33i3i403033234
203
3
43i3i403032222
304
2
34i6i703043234
304
4
i6i7042()2()2()2()z z z z z z z z z z z z z z z z z z z z z z A F R k R g R R g A R g A R A F R R g R R g A R g A R A F R R g R R g A R g A R A F R R R R g A R μμμμμμμμμμμ=
+
⎛⎞
++
+⎜⎟+⎝⎠
+
⎛⎞
++
+⎜⎟+⎝⎠
+
⎛⎞
++
+⎜⎟+⎝⎠
++
+3
40042)z z z z g A g A R μμ⎛⎞
⎜⎟+⎝⎠
（
（40）
其中
g g g 0
我的偶像
r z z r z z R R R −
+=== （41）100000用英语怎么说
0g g g g g 0
x r y r r x r x r y y R R R R R =−
+−+===== （42）
