2024年5月5日发(作者:奕雪晴)
第15卷第7期
电 机 与 控 制 学 报
Vol_15 No.7
2011年7月
ELECTRI C MACHINES AND C0NTR0L
July 2011
永磁电动式磁悬浮装置的研究
秦伟, 范瑜, 朱熙, 吕刚, 李硕, 张晓东
(北京交通大学电气工程学院,北京100044)
摘要:针对目前电磁式磁悬浮系统悬浮气隙小且需要复杂的闭环控制等问题,永磁电动式磁悬浮
系统具有悬浮气隙大和本征稳定的特性不需要复杂的闭环控制即可实现稳定的悬浮等优点。从电
磁场的角度出发,建立永磁电动式磁悬浮系统的数学模型,研究其动态特性,得到悬浮力和制动力
的解析表达式;建立Ansofl有限元模型,分析得出了系统的磁场分布,电磁力等的分析结果;搭建实
验平台对永磁电动式系统的基本特性进行研究,主要是悬浮力和制动力的测试;比较解析解、数值
解和试验结果,表明本文建立的数学模型是合理的、有效的。
关键词:电动式磁悬浮;永磁体;悬浮力;制动力;有限元
中图分类号:TM 302 文献标志码:A 文章编号:1007—449X(2011)07—0077—05
Analysis of permanent-magnet electro-dynamic maglev device
QIN Wei,FAN Yu,ZHU Xi,埔Gang,LI Shuo,ZHANG Xiao—dong
(School of Electrical Engineering,Beijing Jiaotong University,Beijing 100044,China)
Abstract:As an electro—magnetic levitation device,permanent—magnet electro—dynamic maglev system
(EMS)has many advantages.Relative to EMS,it is not necessary for the magnetic levitation to possess
a complex closed—loop control and have a wide air—gap.The electromagnetic part of a permanent magnet
electro—dynamic system analytical equation was established and the solution was given.The results of the
distribution of magnetic field,solid loss and the electromagnetic levitation force were shown.Comparisons
of simulation solution,analytical solution and measurement reveal that two methods are consistent,and
the analytical model is effective and efifcient.
Key words:electro—dynamic maglev;permanent magnet;litf force;drag force;finite—element
0 引 言
相对于其他形式的磁悬浮列车更有独到之处:结构
简单、悬浮气隙大,造价低,可以同时应用于市内和
磁悬浮列车是一种新型地面无接触高速交通运
城际交通,同时它并不需要复杂的闭环控制,因此将
输工具,其具有无污染、噪音小、能耗低、易维护、安
是磁悬浮发展的重要方向之一¨j。
全性高等特点。永磁电动式悬浮装置是通过车载永
本论文从研究导体板涡流分布出发,在总结前
磁体在导体轨道中感应出涡流实现悬浮和导向的,
人相关研究成果的基础上,应用Maxwell电磁场方
收稿日期:2010—11—13
基金项目:国家自然科学基金(51077003,50807004);教育部博士点基金(2009009110025)
作者简介:秦伟(1986一),男,博士研究生,研究方向为特种电机控制和电磁场分析;
范瑜(1954一),男,教授,博士生导师,研究方向为磁悬浮技术和特种电机,电力系统控制与建模;
朱 熙(1983一),男,博士研究生,研究方向为特种电机和电磁场计算;
吕 刚(1976一),男,博士,讲师,研究方向为直线电机分析与控制、机车电力传动;
李硕(1984一),男,博士研究生,研究方向为特种电机和电磁场计算;
张晓东(1983一),男,硕士研究生,研究方向为特种电机和电磁场计算。
78 电机与控制学报 第15卷
程,通过数学推导和分析,得出磁轮电磁部分的磁场
表达式,进而得出考量悬浮质量的悬浮力的表达式,
并用数值方法验证解析方法的正确性。针对永磁体
和导体板磁悬浮系统的特点,提出一些假设条件,将
复杂的磁悬浮系统简化为二维系统,得到悬浮力、磁
阻力的解析表达式。把计算结果和用有限元计算的
数值计算结果进行了对比,最后设计制造试验装置
来验证磁悬浮系统动态特性的分析结果。
1 电动式磁悬浮装置的结构分析
本文研究的磁悬浮机构的模型图如图1所示,
图中的矩形导体板,其上方不连续的矩形永磁体,导
体板静止,永磁体沿 轴运动,上述装置的工作原理
就是非铁磁性导体盘(本文选择的电工铝)和永磁
体产生相对运动,在导体板中感应出涡流磁场,涡流
磁场和永磁体磁场相互作用实现悬浮。
图1永磁电动式磁悬浮系统2D模型
Fig.1 Permanent-magnet electro-dynamic maglev 2-model
2数学模型和解析解
首先是对模型作如下假设 ]:
1)气隙磁场强度沿z轴是恒定的值;
2)永磁体沿 轴运动,导体板中的感应电流只
有 轴分量;
3)各层的物理参数是均匀的,各向同性的。
基于以上假设建立坐标系如图1所示的坐标
系,坐标原点固定在厚度为D,在 方向为无限长的
导体板上表面;永磁体是位于导体板上表面h处并
沿 轴以速度 匀速运动,即h为气隙长度。
在导体区域内,根据Maxwell方程 有
V×B=IX0.,,、
V×E一
.,: E。 J
联立上述方程组可得
7 2B=ix
oo"iaB (1)
是永磁体与导轨之间的相对线速度,则有
OB
= =
以
。
(2)
将式(2)代人式(1)中可得
02B
O 2
+
+ 2
BB
o0"19
035
:A、1
i,
,
(3)j,
d
}.E
_
=Ix
Cl a
式中:Al=IXoo'dv;IX0=41T×10 为空气磁导率,or=
3.8×10 为铝的电导率。
在磁体运动方向上,电磁场的分布为一个行
波,由于相对运动速度远小于光速,所以认为电磁场
是似稳电磁场,也就是说沿 方向运动,变量带有
e
““七 因子,其中 =2 ,k=2盯 ,其中A。是
永磁体有效长度。在导体板中的电磁场是时谐场,
变换成复向量形式有
B =b (Y)e‘‘ 一 。
代入式(3)中可得
:
( 。一i从 )6 。
o’
解上式可得
B =(A1e一 +A2e融)ei‘ 卜h’。 (4)
同理可得
B =(A3e一 +A4e )e ‘“ 。 (5)
又有
堡 =0
d d o
展开上式可得
[(一i 1一肛A3)e一 +(一i 2+ A4)e ]e ‘ 一h’=0。
分解上式可得
ikAl+ 3=0, (6)
一
ikA2+ 4=0。 (7)
在气隙中建立Maxwell方程,可得
02B
+ :0
a ay
可得
B =(A5e一 +A6e ’ )e’ h , (8)
B =(A7e一 + 8e )ei‘ ’。 (9)
同理 + :0,可得
。 oy
ikA5+JB1A7=0, (10)
一
ikA6+JBlA8=0。 (11)
永磁体区域,以均匀速度 运行中的磁悬浮列
车的车载永磁体,可以视为工作在一个静态系统中,
这里不考虑列车在起浮和加速过程中,由于系统悬
浮力变化造成悬浮气隙中磁场强度的变化。永磁场
强度分布如图2所示。
第7期 秦伟等:永磁电动式磁悬浮装置的研究 79
(A・mm )
厂 ] 厂_]
U
图2永磁场强度分布
Fig.2 Distribution of magnetic strength
为了便于计算我首先将间断的永磁体等效成一
个连续的永磁体面,对于每一块永磁体可知 ]
B=/-to(日+ )= 0( + ),
式中: 为真空中的磁导率; 为永磁体的相对磁
导率;Mo为剩余磁化强度,为常数。其方向为导体
板相垂直的方向。
。:
旦一 :
,
O 0
式中 是剩磁感应强度。
则由傅里叶变化可得出沿Y轴方向的磁化强度为
一
,y=∑ cos[n(wt一 )],
式中 是各个谐波分量的幅值,可以表示为
: sin ,n:1,3,5,…,
其中a为磁体占空比。为求解简单方便取复向量为
: e “一
。
在导体板表面和磁体表面,切线方向的称连续
,
B 相等,视导体板与空气的磁导率相同,都等于
真空磁导率,结合公式可得
A1+A2=A5+A6, (12)
A3+A4=A7+A8, (13)
(A5e一培+ 6e姆)= 。 (14)
其中g=h。
假设涡流分布厚度为导体板厚度D,将系统在
导体板下表面“封闭”,这个假设有其合理性的。因
为,在一个实际磁悬浮列车运营系统中,如果透过导
体板的漏磁通比较大,就会对环境和乘客造成不利
影响,而如果导体板比涡流的分布厚度厚很多,说明
导体板尺寸优化设计不好,造成了材料上的浪费,所
以最优的情况就是涡流分布基本等于导体板厚度,
也就是说导体板大致将磁体磁场完全屏蔽,即在导
体板的下表面B =0,即可以得出_8]
(A3e一 +A4e )=0。 (15)
综合方程(6),(7),(10),(11),(12),(13),
(14),(15)可以得到求解得出:A ,A ,A ,A ,A ,A6,
A,,A 。即可以求出气隙磁场,通过代人化解可得到
导体板表面级Y=0的气隙磁场为
-
iA
l
k
 ̄
=
。sh培。 c_h),(16)
一
iA
c
-
l/t:s
e 一h
。
(17)
可以求得在导体板表面的电流密度为
dvBr- - iAlk ̄
ei( )
。
(18)
式中:AI= 0trvd;k =ikMp ;d导体板厚度。导体板
中的涡流和磁场相互作用产生悬浮和阻力力,根据
洛伦兹定律 有
F=J×B 。
则悬浮力为
,
2
。
( s)2_
Sln h
Kg+ U
_d
v cosh。
n
tr
l ‘Kg
(19)
悬浮阻力为
Fx tzola ̄oo'dk。2 _
。
(20)
sl
n
n +l t‘n
o'ctv
J
cosh
g
由上述公式可知悬浮力随着磁体和导体板的相
对速度的增加而增加,随着速度的逐渐增大变化量
逐渐减小;而导体板悬浮的电磁阻力随相对速度
的变化是非线性的,在低速时电磁阻力随相对速度
的升高而升高,但是随着相对速度的继续增加,金属
板中磁场变化频率升高,由于涡流的集肤效应,电磁
阻力反而下降,最后趋于一个渐近值。
3有限元计算
利用有限元计算软件Ansofl进行数值计算,为
了更好的反映悬浮系统的动态特性选择Ansofl瞬态
求解器,建立参数如下的二维模型,具体参数为:永
磁体长为18 mm;高为5 mm;宽为18 mm;极对数为
4;极距为36 mm;次级铝板厚为8 mm;电导率为3.8×
10 s/m;气隙长为10 mm。然后经过指定材料属性,
指定边界条件和激励源,设定求解参数,求解分析等
步骤后将得到在不同转速下的磁场分布图如图3和
图4所示。
通过建立永磁电动式磁悬浮系统的简化二维模
型,计算了在不同气隙下,不同永磁体相对速度下的
悬浮力和磁阻力。比较图3(e)、图3(d)可以看出
随着气隙的增加,与导体板铰链的磁力线减少,从而
使得悬浮力减小;比较图3(a)、图3(b)、图3(C)可
看出随着速度在增加进入导体板中的磁力线增加,
从而为悬浮系统提供更大的悬浮力,这与解析表达
80 电机与控制学报 第l5卷
式的分析结果是一致的。比较图4和图5,发现在
永磁体在导体板表面感应出的磁场方向自身相同,
根据同极相斥实现悬浮,这正是电动式磁悬浮的悬
4试验验证
为了验证解析计算和有限元计算的正确性,设
浮原理。
£_
(a)v=0,g=0
鬈—・x
(h)"=l0 rn/s.g=10 mm
—一—一—一… —毫
:二二: .;
((1)v=20m/s.g=20mm
图3磁通线分布图
Fig.3 Flux distribution
1.5
1.0
0.5
0.0
0.5
1.0
1.5
图4 t,=20 m/s。永磁体表面磁密
Fig.4 tJ=20m/s。the B distribution of the magnet surface
0.25
0.20
0.15
0.10
-r O.O5
0.00
0.O5
0.10
一
O.15
0.20
O
距离/mm
图5"=20 m/s。导体板表面磁密
Fig.5 l,=20 m/s,the B distribution of hte aluminum plate
计制造了永磁电动式磁悬浮的试验装置,本装置优
点是磁体下气隙均匀,可以做出任意厚度,特别是可
以非常方便地进行不同层数金属片叠加安装,便于
比较实验,非常理想。离心力径向造成变形很小,相
对转动速度选择范围大。缺点是沿径向速度分布不
再均匀。
由于次级导体铝板用感应出的涡流场与永磁体
相互作用使得次级导体板有相对运动,因此导体板
受到了沿导体板切线方向的旋转的力和使其悬浮的
垂直力。悬浮力是用拉压传感器测量。本装置都采
用非磁性材料,硬铝或不锈钢,避免因被磁体磁化带
来受力干扰。为了能方便地换装内部的金属片,壳
体设计成可以拆分的几部分,以螺栓或螺钉
联接 。
首先要把测得的转速换算成磁体中心线的线速
度,传感器测到的转矩压力换算成磁体中心线上的
等效阻力。试验系统的数据记录是非连续的,为一
系列的散点,由于装置的直径和电机转速的限制,同
时还有拉压传感器的测量范围限制,本文只测试了
速度在3~20 m/s之间的悬浮力和阻力,如图6所
示,其中F 、F 为解析计算结果;Fx 、F 为试验测量
结果;F F以为有限元计算结果,由图6看出实验得
到的数值要比解析解和有限元小,这是 为解析和
有限元数值分析都是基于简化的二维模型,二维模
型中认为磁场在 方向是均匀的,导致计算误差,但
是可以看出测试和解析结果和有限元分析的结果基
本趋势是一致的。
60
50
4O
蚕3o
20
1O
O
0 20 40 60 8tJ lIll】
/(m・s。 )
图6悬浮力和阻力的计算对比图
Fig.6 Lift and drag-force
5 结 论
1)本文提出的永磁磁场电动式磁悬浮系统在
第7期 秦伟等:永磁电动式磁悬浮装置的研究 81
一
定的相对速度下能够产生足够的悬浮力,实现稳
定的电动式磁悬浮。使用解析方法和有限元法得到
了导体板和气隙的磁场分布,推导出了悬浮法向力
和悬浮阻力的计算表达式,为永磁电动式磁悬浮系
统的特性分析和设计提供了参考。
2)永磁电动式磁悬浮装置能够使非铁磁性导
体次级静止稳定悬浮的特点使其能够应用到金属悬
浮熔炼等其他场合,弥补普遍采用EMS悬浮熔炼技
术只能用于铁磁性金属的空白。
3)该装置产生稳定悬浮力的同时会产生较大
的推进阻力,故最佳的悬浮力和阻力之比是装置设
计和运行的重点之一。该装置的悬浮性能仍存在较
大的提升空间,如采用永磁Halbach阵列,提高气隙
磁密,增大悬浮力,减小永磁体另外一侧磁密降低对
外部的辐射。
参考文献:
[1] YAN Luguang.Development and application of the magnet tech-
nology in China[C]//Proceedings of Fitfeenth International Con—
ference on Magnet Technology。October 20—24,1997,Beijing,
China.1997:30—35.
[2]LEQUESNE B,LIU Buyun,NEHL T W.Eddy current machines
with permanent magnets and solid rotors[J].IEEE Transactions on
Industry Applications,1997,33(5):1289—1294.
[3]NEHLTW,LEQUESNE B,GANGLA V,et a1.Nonlineartwo—
dimensional finite element modeling of permanent magnet eddy
current couplings and brakes[J].IEEE Transactions on Magnet—
its,1994,3o(5):3000—30o3.
[4] 吕刚,范瑜,马云双,等.直线感应电机推力和法向力的解析计
算与分析[J].电机与控制学报,2010,14(3):77—82.
L0 Gang,FAN Yu,MA Yunshuang,et a1.Characteristic analysis
of linear induction traction motor for urban rail transit f J].Electri.
cal Machine and Contorl,2010,14(3):77—82.
[5] 唐永春.超导磁浮直线电机及其控制和涡流制动系统研究
[D].浙江:浙江大学电气工程学院,2006:35—39.
[6] 李春生,杜玉梅,夏平畴,等.直线型Halbach磁体和导体板构
成的电动式磁悬浮系统的分析及实验[j].电工技术学报,
2009,24(1):18—22.
LI Chunsheng,DU Yumei,XIA Pingchou et a1.Analysis and ex—
perimental testing of eds maglev with linear Halbach and conduc—
ting sheet[J].Transactions of China Eleetortechnical Society,
2009.24(1):18—22.
[7]
EDWARD J D,JAYAWANT B V,DAWSON W R C,et a1.Per—
manent magnet linear eddy--current brake with a nonmagnetic reac--
tion plate[J].IEE Proceedings—Electfc Power Applications 1999,
146(6):627—631.
[8] 王厚生,永磁电动式导体板磁悬浮列车轨道结构及柑关研究
[D].北京:中国科学院电工研究所,2004:49—55.
[9] 吕刚,孙守光,马云双.城轨交通中直线感应牵引电机的效率
最优控制[J].电机与控制学报,2009,13(4):490—495,
L12 Gang,SUN Shouguang,MA Yunshuang.Efifciency optimal
control of linear induction motor for urban rail transit『J].Electri.
cal Machine and Contorl,2009,13(4):490—495.
(编辑:张诗阁)
2024年5月5日发(作者:奕雪晴)
第15卷第7期
电 机 与 控 制 学 报
Vol_15 No.7
2011年7月
ELECTRI C MACHINES AND C0NTR0L
July 2011
永磁电动式磁悬浮装置的研究
秦伟, 范瑜, 朱熙, 吕刚, 李硕, 张晓东
(北京交通大学电气工程学院,北京100044)
摘要:针对目前电磁式磁悬浮系统悬浮气隙小且需要复杂的闭环控制等问题,永磁电动式磁悬浮
系统具有悬浮气隙大和本征稳定的特性不需要复杂的闭环控制即可实现稳定的悬浮等优点。从电
磁场的角度出发,建立永磁电动式磁悬浮系统的数学模型,研究其动态特性,得到悬浮力和制动力
的解析表达式;建立Ansofl有限元模型,分析得出了系统的磁场分布,电磁力等的分析结果;搭建实
验平台对永磁电动式系统的基本特性进行研究,主要是悬浮力和制动力的测试;比较解析解、数值
解和试验结果,表明本文建立的数学模型是合理的、有效的。
关键词:电动式磁悬浮;永磁体;悬浮力;制动力;有限元
中图分类号:TM 302 文献标志码:A 文章编号:1007—449X(2011)07—0077—05
Analysis of permanent-magnet electro-dynamic maglev device
QIN Wei,FAN Yu,ZHU Xi,埔Gang,LI Shuo,ZHANG Xiao—dong
(School of Electrical Engineering,Beijing Jiaotong University,Beijing 100044,China)
Abstract:As an electro—magnetic levitation device,permanent—magnet electro—dynamic maglev system
(EMS)has many advantages.Relative to EMS,it is not necessary for the magnetic levitation to possess
a complex closed—loop control and have a wide air—gap.The electromagnetic part of a permanent magnet
electro—dynamic system analytical equation was established and the solution was given.The results of the
distribution of magnetic field,solid loss and the electromagnetic levitation force were shown.Comparisons
of simulation solution,analytical solution and measurement reveal that two methods are consistent,and
the analytical model is effective and efifcient.
Key words:electro—dynamic maglev;permanent magnet;litf force;drag force;finite—element
0 引 言
相对于其他形式的磁悬浮列车更有独到之处:结构
简单、悬浮气隙大,造价低,可以同时应用于市内和
磁悬浮列车是一种新型地面无接触高速交通运
城际交通,同时它并不需要复杂的闭环控制,因此将
输工具,其具有无污染、噪音小、能耗低、易维护、安
是磁悬浮发展的重要方向之一¨j。
全性高等特点。永磁电动式悬浮装置是通过车载永
本论文从研究导体板涡流分布出发,在总结前
磁体在导体轨道中感应出涡流实现悬浮和导向的,
人相关研究成果的基础上,应用Maxwell电磁场方
收稿日期:2010—11—13
基金项目:国家自然科学基金(51077003,50807004);教育部博士点基金(2009009110025)
作者简介:秦伟(1986一),男,博士研究生,研究方向为特种电机控制和电磁场分析;
范瑜(1954一),男,教授,博士生导师,研究方向为磁悬浮技术和特种电机,电力系统控制与建模;
朱 熙(1983一),男,博士研究生,研究方向为特种电机和电磁场计算;
吕 刚(1976一),男,博士,讲师,研究方向为直线电机分析与控制、机车电力传动;
李硕(1984一),男,博士研究生,研究方向为特种电机和电磁场计算;
张晓东(1983一),男,硕士研究生,研究方向为特种电机和电磁场计算。
78 电机与控制学报 第15卷
程,通过数学推导和分析,得出磁轮电磁部分的磁场
表达式,进而得出考量悬浮质量的悬浮力的表达式,
并用数值方法验证解析方法的正确性。针对永磁体
和导体板磁悬浮系统的特点,提出一些假设条件,将
复杂的磁悬浮系统简化为二维系统,得到悬浮力、磁
阻力的解析表达式。把计算结果和用有限元计算的
数值计算结果进行了对比,最后设计制造试验装置
来验证磁悬浮系统动态特性的分析结果。
1 电动式磁悬浮装置的结构分析
本文研究的磁悬浮机构的模型图如图1所示,
图中的矩形导体板,其上方不连续的矩形永磁体,导
体板静止,永磁体沿 轴运动,上述装置的工作原理
就是非铁磁性导体盘(本文选择的电工铝)和永磁
体产生相对运动,在导体板中感应出涡流磁场,涡流
磁场和永磁体磁场相互作用实现悬浮。
图1永磁电动式磁悬浮系统2D模型
Fig.1 Permanent-magnet electro-dynamic maglev 2-model
2数学模型和解析解
首先是对模型作如下假设 ]:
1)气隙磁场强度沿z轴是恒定的值;
2)永磁体沿 轴运动,导体板中的感应电流只
有 轴分量;
3)各层的物理参数是均匀的,各向同性的。
基于以上假设建立坐标系如图1所示的坐标
系,坐标原点固定在厚度为D,在 方向为无限长的
导体板上表面;永磁体是位于导体板上表面h处并
沿 轴以速度 匀速运动,即h为气隙长度。
在导体区域内,根据Maxwell方程 有
V×B=IX0.,,、
V×E一
.,: E。 J
联立上述方程组可得
7 2B=ix
oo"iaB (1)
是永磁体与导轨之间的相对线速度,则有
OB
= =
以
。
(2)
将式(2)代人式(1)中可得
02B
O 2
+
+ 2
BB
o0"19
035
:A、1
i,
,
(3)j,
d
}.E
_
=Ix
Cl a
式中:Al=IXoo'dv;IX0=41T×10 为空气磁导率,or=
3.8×10 为铝的电导率。
在磁体运动方向上,电磁场的分布为一个行
波,由于相对运动速度远小于光速,所以认为电磁场
是似稳电磁场,也就是说沿 方向运动,变量带有
e
““七 因子,其中 =2 ,k=2盯 ,其中A。是
永磁体有效长度。在导体板中的电磁场是时谐场,
变换成复向量形式有
B =b (Y)e‘‘ 一 。
代入式(3)中可得
:
( 。一i从 )6 。
o’
解上式可得
B =(A1e一 +A2e融)ei‘ 卜h’。 (4)
同理可得
B =(A3e一 +A4e )e ‘“ 。 (5)
又有
堡 =0
d d o
展开上式可得
[(一i 1一肛A3)e一 +(一i 2+ A4)e ]e ‘ 一h’=0。
分解上式可得
ikAl+ 3=0, (6)
一
ikA2+ 4=0。 (7)
在气隙中建立Maxwell方程,可得
02B
+ :0
a ay
可得
B =(A5e一 +A6e ’ )e’ h , (8)
B =(A7e一 + 8e )ei‘ ’。 (9)
同理 + :0,可得
。 oy
ikA5+JB1A7=0, (10)
一
ikA6+JBlA8=0。 (11)
永磁体区域,以均匀速度 运行中的磁悬浮列
车的车载永磁体,可以视为工作在一个静态系统中,
这里不考虑列车在起浮和加速过程中,由于系统悬
浮力变化造成悬浮气隙中磁场强度的变化。永磁场
强度分布如图2所示。
第7期 秦伟等:永磁电动式磁悬浮装置的研究 79
(A・mm )
厂 ] 厂_]
U
图2永磁场强度分布
Fig.2 Distribution of magnetic strength
为了便于计算我首先将间断的永磁体等效成一
个连续的永磁体面,对于每一块永磁体可知 ]
B=/-to(日+ )= 0( + ),
式中: 为真空中的磁导率; 为永磁体的相对磁
导率;Mo为剩余磁化强度,为常数。其方向为导体
板相垂直的方向。
。:
旦一 :
,
O 0
式中 是剩磁感应强度。
则由傅里叶变化可得出沿Y轴方向的磁化强度为
一
,y=∑ cos[n(wt一 )],
式中 是各个谐波分量的幅值,可以表示为
: sin ,n:1,3,5,…,
其中a为磁体占空比。为求解简单方便取复向量为
: e “一
。
在导体板表面和磁体表面,切线方向的称连续
,
B 相等,视导体板与空气的磁导率相同,都等于
真空磁导率,结合公式可得
A1+A2=A5+A6, (12)
A3+A4=A7+A8, (13)
(A5e一培+ 6e姆)= 。 (14)
其中g=h。
假设涡流分布厚度为导体板厚度D,将系统在
导体板下表面“封闭”,这个假设有其合理性的。因
为,在一个实际磁悬浮列车运营系统中,如果透过导
体板的漏磁通比较大,就会对环境和乘客造成不利
影响,而如果导体板比涡流的分布厚度厚很多,说明
导体板尺寸优化设计不好,造成了材料上的浪费,所
以最优的情况就是涡流分布基本等于导体板厚度,
也就是说导体板大致将磁体磁场完全屏蔽,即在导
体板的下表面B =0,即可以得出_8]
(A3e一 +A4e )=0。 (15)
综合方程(6),(7),(10),(11),(12),(13),
(14),(15)可以得到求解得出:A ,A ,A ,A ,A ,A6,
A,,A 。即可以求出气隙磁场,通过代人化解可得到
导体板表面级Y=0的气隙磁场为
-
iA
l
k
 ̄
=
。sh培。 c_h),(16)
一
iA
c
-
l/t:s
e 一h
。
(17)
可以求得在导体板表面的电流密度为
dvBr- - iAlk ̄
ei( )
。
(18)
式中:AI= 0trvd;k =ikMp ;d导体板厚度。导体板
中的涡流和磁场相互作用产生悬浮和阻力力,根据
洛伦兹定律 有
F=J×B 。
则悬浮力为
,
2
。
( s)2_
Sln h
Kg+ U
_d
v cosh。
n
tr
l ‘Kg
(19)
悬浮阻力为
Fx tzola ̄oo'dk。2 _
。
(20)
sl
n
n +l t‘n
o'ctv
J
cosh
g
由上述公式可知悬浮力随着磁体和导体板的相
对速度的增加而增加,随着速度的逐渐增大变化量
逐渐减小;而导体板悬浮的电磁阻力随相对速度
的变化是非线性的,在低速时电磁阻力随相对速度
的升高而升高,但是随着相对速度的继续增加,金属
板中磁场变化频率升高,由于涡流的集肤效应,电磁
阻力反而下降,最后趋于一个渐近值。
3有限元计算
利用有限元计算软件Ansofl进行数值计算,为
了更好的反映悬浮系统的动态特性选择Ansofl瞬态
求解器,建立参数如下的二维模型,具体参数为:永
磁体长为18 mm;高为5 mm;宽为18 mm;极对数为
4;极距为36 mm;次级铝板厚为8 mm;电导率为3.8×
10 s/m;气隙长为10 mm。然后经过指定材料属性,
指定边界条件和激励源,设定求解参数,求解分析等
步骤后将得到在不同转速下的磁场分布图如图3和
图4所示。
通过建立永磁电动式磁悬浮系统的简化二维模
型,计算了在不同气隙下,不同永磁体相对速度下的
悬浮力和磁阻力。比较图3(e)、图3(d)可以看出
随着气隙的增加,与导体板铰链的磁力线减少,从而
使得悬浮力减小;比较图3(a)、图3(b)、图3(C)可
看出随着速度在增加进入导体板中的磁力线增加,
从而为悬浮系统提供更大的悬浮力,这与解析表达
80 电机与控制学报 第l5卷
式的分析结果是一致的。比较图4和图5,发现在
永磁体在导体板表面感应出的磁场方向自身相同,
根据同极相斥实现悬浮,这正是电动式磁悬浮的悬
4试验验证
为了验证解析计算和有限元计算的正确性,设
浮原理。
£_
(a)v=0,g=0
鬈—・x
(h)"=l0 rn/s.g=10 mm
—一—一—一… —毫
:二二: .;
((1)v=20m/s.g=20mm
图3磁通线分布图
Fig.3 Flux distribution
1.5
1.0
0.5
0.0
0.5
1.0
1.5
图4 t,=20 m/s。永磁体表面磁密
Fig.4 tJ=20m/s。the B distribution of the magnet surface
0.25
0.20
0.15
0.10
-r O.O5
0.00
0.O5
0.10
一
O.15
0.20
O
距离/mm
图5"=20 m/s。导体板表面磁密
Fig.5 l,=20 m/s,the B distribution of hte aluminum plate
计制造了永磁电动式磁悬浮的试验装置,本装置优
点是磁体下气隙均匀,可以做出任意厚度,特别是可
以非常方便地进行不同层数金属片叠加安装,便于
比较实验,非常理想。离心力径向造成变形很小,相
对转动速度选择范围大。缺点是沿径向速度分布不
再均匀。
由于次级导体铝板用感应出的涡流场与永磁体
相互作用使得次级导体板有相对运动,因此导体板
受到了沿导体板切线方向的旋转的力和使其悬浮的
垂直力。悬浮力是用拉压传感器测量。本装置都采
用非磁性材料,硬铝或不锈钢,避免因被磁体磁化带
来受力干扰。为了能方便地换装内部的金属片,壳
体设计成可以拆分的几部分,以螺栓或螺钉
联接 。
首先要把测得的转速换算成磁体中心线的线速
度,传感器测到的转矩压力换算成磁体中心线上的
等效阻力。试验系统的数据记录是非连续的,为一
系列的散点,由于装置的直径和电机转速的限制,同
时还有拉压传感器的测量范围限制,本文只测试了
速度在3~20 m/s之间的悬浮力和阻力,如图6所
示,其中F 、F 为解析计算结果;Fx 、F 为试验测量
结果;F F以为有限元计算结果,由图6看出实验得
到的数值要比解析解和有限元小,这是 为解析和
有限元数值分析都是基于简化的二维模型,二维模
型中认为磁场在 方向是均匀的,导致计算误差,但
是可以看出测试和解析结果和有限元分析的结果基
本趋势是一致的。
60
50
4O
蚕3o
20
1O
O
0 20 40 60 8tJ lIll】
/(m・s。 )
图6悬浮力和阻力的计算对比图
Fig.6 Lift and drag-force
5 结 论
1)本文提出的永磁磁场电动式磁悬浮系统在
第7期 秦伟等:永磁电动式磁悬浮装置的研究 81
一
定的相对速度下能够产生足够的悬浮力,实现稳
定的电动式磁悬浮。使用解析方法和有限元法得到
了导体板和气隙的磁场分布,推导出了悬浮法向力
和悬浮阻力的计算表达式,为永磁电动式磁悬浮系
统的特性分析和设计提供了参考。
2)永磁电动式磁悬浮装置能够使非铁磁性导
体次级静止稳定悬浮的特点使其能够应用到金属悬
浮熔炼等其他场合,弥补普遍采用EMS悬浮熔炼技
术只能用于铁磁性金属的空白。
3)该装置产生稳定悬浮力的同时会产生较大
的推进阻力,故最佳的悬浮力和阻力之比是装置设
计和运行的重点之一。该装置的悬浮性能仍存在较
大的提升空间,如采用永磁Halbach阵列,提高气隙
磁密,增大悬浮力,减小永磁体另外一侧磁密降低对
外部的辐射。
参考文献:
[1] YAN Luguang.Development and application of the magnet tech-
nology in China[C]//Proceedings of Fitfeenth International Con—
ference on Magnet Technology。October 20—24,1997,Beijing,
China.1997:30—35.
[2]LEQUESNE B,LIU Buyun,NEHL T W.Eddy current machines
with permanent magnets and solid rotors[J].IEEE Transactions on
Industry Applications,1997,33(5):1289—1294.
[3]NEHLTW,LEQUESNE B,GANGLA V,et a1.Nonlineartwo—
dimensional finite element modeling of permanent magnet eddy
current couplings and brakes[J].IEEE Transactions on Magnet—
its,1994,3o(5):3000—30o3.
[4] 吕刚,范瑜,马云双,等.直线感应电机推力和法向力的解析计
算与分析[J].电机与控制学报,2010,14(3):77—82.
L0 Gang,FAN Yu,MA Yunshuang,et a1.Characteristic analysis
of linear induction traction motor for urban rail transit f J].Electri.
cal Machine and Contorl,2010,14(3):77—82.
[5] 唐永春.超导磁浮直线电机及其控制和涡流制动系统研究
[D].浙江:浙江大学电气工程学院,2006:35—39.
[6] 李春生,杜玉梅,夏平畴,等.直线型Halbach磁体和导体板构
成的电动式磁悬浮系统的分析及实验[j].电工技术学报,
2009,24(1):18—22.
LI Chunsheng,DU Yumei,XIA Pingchou et a1.Analysis and ex—
perimental testing of eds maglev with linear Halbach and conduc—
ting sheet[J].Transactions of China Eleetortechnical Society,
2009.24(1):18—22.
[7]
EDWARD J D,JAYAWANT B V,DAWSON W R C,et a1.Per—
manent magnet linear eddy--current brake with a nonmagnetic reac--
tion plate[J].IEE Proceedings—Electfc Power Applications 1999,
146(6):627—631.
[8] 王厚生,永磁电动式导体板磁悬浮列车轨道结构及柑关研究
[D].北京:中国科学院电工研究所,2004:49—55.
[9] 吕刚,孙守光,马云双.城轨交通中直线感应牵引电机的效率
最优控制[J].电机与控制学报,2009,13(4):490—495,
L12 Gang,SUN Shouguang,MA Yunshuang.Efifciency optimal
control of linear induction motor for urban rail transit『J].Electri.
cal Machine and Contorl,2009,13(4):490—495.
(编辑:张诗阁)