Application of Digital One-Cycle Control for Current in Electromagnetic Suspension System
-
摘要: 在串级控制的电磁悬浮系统中,电流环的响应速度和精度对整个悬浮控制起着至关重要的作用. 为了加快悬浮系统电流环的响应速度以及减小跟随误差,基于TMS320F28335设计了EMS (electromagnetic suspension system)的数字单周期控制(digital one-cycle control,D-OCC)电流控制器. 以悬浮斩波器为研究对象,建立起D-OCC的数学模型,对额定悬浮工作点处斩波器电流的D-OCC算法进行了详细推导;通过Simulink平台对算法进行仿真验证,并将D-OCC的电流环投入到实际悬浮系统中进行悬浮实验. 实验结果表明:对频率为5 Hz,幅值为3 A的方波信号进行跟随时,传统PID控制在方波上升沿和下降沿均存在一定的超调,且稳定后存在不小于20 mA的跟随误差,D-OCC在调节过程中不存在超调,且稳定后没有跟随误差,说明D-OCC算法能够实现对指令电流快速、准确跟随;采用电流环D-OCC的悬浮系统起浮过程需要约0.4 s的调整时间,并且悬浮稳定后可以克服50%荷载扰动和1.5 mm气隙扰动,说明该方法可以实现系统稳定悬浮,且具有较强的鲁棒性能.Abstract: In an electromagnetic suspension system (EMS) with cascade control, the response speed and accuracy of the current loop play a vital role in the overall suspension control. A current controller was designed for digital one-cycle control (D-OCC) based on TMS320F28335 for the suspension system, in order to improve the response speed of the current loop and reduce the following error. Taking a chopper circuit as the research object, the D-OCC’s mathematical model was established and the control algorithm of rated suspension point was derived in detail. The algorithm was validated through Simulink, and the current loop controlled by D-OCC was incorporated into the actual EMS. Experimental results show that the conventional PID control has a certain overshoot at the wave’s rising and falling edges when a square wave signal with frequency of 5 Hz and amplitude of 3 A is followed. Moreover, there is a following error of not less than 20 mA after the stabilization. However, there is no overshoot in the adjustment process and no following error after stabilization, when following the same wave with the method of D-OCC. This means that the D-OCC algorithm can follow the instruction current quickly and accurately. The floating process of the suspension system using the D-OCC current loop requires an adjustment time of about 0.4 s, and the system can overcome the 50% load disturbance and the 1.5 mm air gap disturbance; this shows that the method can achieve the system stable levitation and has strong robust performance.
-
Key words:
- electromagnetic suspension system /
- cascade control /
- current loop /
- chopper /
- digital one-cycle control
-
磁悬浮技术具有无摩擦、微振动、长寿命及高精度等优点,被广泛地应用于高速机械领域. 磁轴承作为核心部件,一般与高速旋转电机配合用于高速旋转机械,如磁悬浮飞轮、磁悬浮控制力矩陀螺及磁悬浮电机等[1-3]. 对于直线运动,通常采用能够实现悬浮及导向功能的电磁铁与直线电机的组合,如磁浮列车、磁浮电梯等. 目前为止,世界很多国家已经对磁浮列车技术研究了较长时间,而且部分实现了商业运行,如中国、德国、日本、韩国等[4-6]. 磁浮列车由电磁铁与轨道功能件之间产生的电磁力支撑,从而实现无接触运行. 电磁铁为悬浮系统的核心执行部件,其电磁力特性决定了列车的承载能力,并与控制器、传感器相互配合实现列车的稳定悬浮. 尤其对于高速运行的磁浮列车,运行载荷更加复杂、苛刻,对承载能力及稳定性要求更高. 因此,电磁力特性分析作为基础研究,对悬浮系统的设计及优化起至关重要作用[7-9].
目前,电磁力建模分析方法主要包含等效磁路法(EMC)及有限元法(FEM),其中,FEM计算精度较高,但效率较低,很难与控制模型联合用于分析系统实时特性. EMC计算速度较高,能够与控制模型联合用于系统实时动态特性分析,但其计算精度较低. 因此,在传统EMC模型中通常会考虑加入补偿系数,通过调整系数校正电磁力结果,使其与FEM结果接近[10-11]. 然而,传统EMC模型仅考虑线性工作区,导磁材料采用恒定的相对磁导率,忽略磁饱和影响,甚至忽略导磁材料磁阻. 这样会导致EMC模型结果在小电流区间内较准确,而在大电流区间就会出现较大偏差[12-13]. 本文在搭建高速磁浮悬浮电磁铁EMC模型时,考虑了导磁材料的磁阻及其非线性. 通过导磁材料B-H (B为磁感应强度;H为磁场强度)曲线的拟合及引入,求解电磁力的准确性大幅度提高,适用范围增加.
高速磁浮列车的悬浮电磁铁共有12个磁极,极性为NS交替,相邻磁极之间通过磁轭连接,磁场经过长定子铁芯形成回路,磁极与长定子之间的磁场产生电磁吸力,实现悬浮功能. 12个磁极分为左、右两组,分别由两个悬浮控制器单独控制,从而形成两个控制回路,每个回路对应两个间隙传感器. 传感器实时监测电磁铁与长定子之间的间隙,并反馈给悬浮控制器,经过控制策略计算,悬浮控制器输出相应电压给悬浮电磁铁,实现动态稳定悬浮.
1. 电磁铁模型搭建
1.1 等效磁路
悬浮电磁铁与长定子的物理模型如图1所示. 其中电磁铁分为左、右两个回路单独控制,为简化计算模型,仅对半个电磁铁进行建模,并忽略两个回路间磁场的影响. 搭建的等效磁路模型如图2所示,磁路中包含了气隙磁阻、漏磁磁阻及导磁材料磁阻. 图中:Raj为磁极与长定子间的气隙磁阻;Rsi为长定子铁芯磁阻;Rei为磁极铁芯及磁轭磁阻;Rli为相邻磁极之间的漏磁磁阻;ϕaj、ϕsi、ϕei及ϕli分别为磁阻Raj、Rsi、Rei及Rli对应的磁通;ϕpj 为磁极磁通;θj为磁极磁动势nI,n为磁极匝数,I为控制回路电流;i=1,2,…,5,j=1,2,…,6.
根据建模需求,定义磁通、磁动势向量为
{ϕ=(ϕ1,ϕl1,ϕ2,ϕl2,⋯,ϕ5,ϕl5)T,ϕs=(ϕ1,ϕ2,⋯,ϕ5)T,ϕa=(ϕ1,ϕ1+ϕ2,⋯,ϕ4+ϕ5,ϕ5)T,ϕe=(ϕ1+ϕl1,ϕ2+ϕl2,⋯,ϕ5+ϕl5)T,ϕp=(ϕ1+ϕl1,ϕ1+ϕl1+ϕ2+ϕl2,⋯,ϕ4+ϕl4+ϕ5+ϕl5,ϕ5+ϕl5)T,θ=(θ1,θ2,⋯,θ6)T. (1) 磁极磁通向量ϕp与磁通向量ϕ的转换关系为
ϕp =Tϕ,
式中:
T=[11111111111111111111]. 在等效磁路中,根据基尔霍夫电压定律,建立关于磁通的方程组,如式(3).
{Ra(j−1)ϕa(j−1)+Rsiϕsi+Rajϕaj+Reiϕei=θj−1+θj,Rliϕli+Reiϕei=θj−1+θj, (3) 式中:j=2,3,⋯,6.
将式(1)表达为矩阵及向量形式为
Aϕ=TTθ=TTnI, (4) 式中:向量n=(n, n, n, n, n, n)T;A为磁阻矩阵,A∈R10×10,
A=[Ra1+Ra2+Rs1+Re1Re1Ra200⋯Re1Rl1+Re1000⋯Ra20Ra2+Ra3+Rs2+Re2Re2Ra3⋯00Re2Rl2+Re20⋯⋮⋮⋮⋮⋮]. 1.2 磁阻计算
1) 磁极与长定子间气隙磁阻
受磁极直线发电机(linear generator, LIG)槽与长定子齿槽结构的影响,磁极与长定子之间的气隙磁通分布较为复杂,如图3所示. 因此,将气隙磁通等效分为主磁通、槽磁通及LIG磁通,分别对应3种磁阻.
气隙磁阻为3种磁阻并联,即
1Raj=1Ra,mj+1Ra,nj+1Ra,Lj, (5) 式中:Ra,mj为主磁通对应的磁阻;Ra,nj为槽磁通对应的磁阻;Ra,Lj为LIG磁通对应的磁阻.
每种磁阻可由式(6)计算.
Ra,oj=sj+hojμ0Aa,oj,o∈[m,n,L], (6) 式中:sj为磁极与长定子齿之间的间隙;μ0为空气磁导率;Aa,oj为相应气隙面积;hoj为额外气隙长度;m、n、L分别对应主磁通、槽磁通和LIG磁通.
2) 相邻磁极间漏磁磁阻
磁极磁通大部分经过长定子回到相邻磁极,小部分未经过长定子而直接回到相邻磁极,该部分磁通为相邻磁极之间的漏磁,对应的磁阻称为漏磁磁阻,可由式(7)计算.
Rli=hliμ0Ali, (7) 式中:Ali为等效气隙面积;hli为相邻磁极间的等效气隙长度.
3) 导磁材料磁阻
对于长定子铁芯、磁极铁芯及磁轭的磁阻,采用分段方式进行求解,且尽可能保证每段的截面积相同. 具体分段如图4所示,其中,磁极及磁轭共分为5段(1~5),长定子分为3段(6~8). 该部分磁阻计算时考虑导磁材料的非线性.
磁极铁芯及磁轭磁阻为
Rei=5∑k=1lei,kμ0μr(ϕei,k)Aei,k; (8) 长定子铁芯磁阻为
Rsi=8∑k=6lsi,kμ0μr(ϕsi,k)Asi,k, (9) 式(8)、(9)中:lei,k, lsi,k分别为悬浮电磁铁侧、长定子侧每段铁芯的长度;Aei,k, Asi,k分别
为悬浮电磁铁侧、长定子侧每段铁芯的截面积;μr为每段铁芯的相对磁导率;ϕsi,k为第i个磁回路中第k段长定子内的磁通量;ϕei,k为第i个回路中第k段磁极铁芯或磁轭内的磁通量. 为在EMC模型中引入导磁材料的非线性,不再将μr简单地设置为恒定值,而是根据每段铁芯的磁通进行计算. 磁极铁芯、磁轭及长定子铁芯均采用硅钢片,牌号为M530-50A,导磁材料的非线性可通过B-H曲线体现,如图5所示.
根据文献[14],B-H曲线可采用式(10)函数进行拟合.
H=α1B+α2Bα3, (10) 式中:α1、α2及α3为拟合函数自变量B的系数,可通过对图5中B-H曲线的拟合确定.
长定子铁芯的相对磁导率可表示为
μr1=Bμ0H=1μ0(α1+α2Bα3−1). (11) 根据磁密、面积及磁通的关系,式(11)可表达为
μr1=1μ0(α1+α2(ϕsiAsi,k)α3−1). (12) 同理可求解磁极铁芯及磁轭的相对磁导率为
μr2=1μ0(α1+α2(ϕeiAei,k)α3−1). (13) 将式(12)、(13)代入式(8)、(9)求解磁极铁芯、磁轭及长定子铁芯磁阻. 通过式(5)~(9)可知,磁阻矩阵A与间隙s及磁通ϕ相关,因此,将其记为A(s,ϕ).
1.3 电路模型
悬浮电磁铁分为两个控制回路,分别由一个控制器进行供电. 控制器的输出为电压,根据控制回路的负载特性转变为相应的负载电流. 在进行电磁力模型与控制模型联合分析时,需要搭建电磁铁的电路模型,其功能是将控制模型的输入电压转变为负载电流,再结合EMC模型计算电磁力. 与EMC模型类似,仅搭建一个回路的模型,6个磁极串联实际可等效为6个电阻与6个电感的串联,如图6所示. 图中:RMj为单个磁极电阻;LMj为单个磁极电感;ULj为单个磁极电感电压;UM为单个控制回路的输入电压.
根据电路模型,单个控制回路的输入电压为
UM=RMI+6∑j=1ULj=RMI+nT∙ϕp, (14) 式中:RM为单个控制回路的总电阻.
整理式(14),控制回路电流可表达为
I=(UM−nT∙ϕp)/RM. (15) 1.4 磁力模型
式(4)给出了磁通与电流关系,式(15)给出了磁通、电流与电压关系. 将式(15)代入式(4),可得磁通与电压的关系为
A(s,ϕ)ϕ=TTn(UM−nT∙ϕp)/RM. (16) 根据磁通关系ϕp=Tϕ,式(16)整理为
TTnnTT∙ϕ = TTnUM−RMA(s,ϕ)ϕ. (17) 经过计算发现,矩阵TTnnTT为奇异矩阵,常微分式(17)很难进行求解. 因此,为便于求解方程,采用中间变量替换原变量.
根据磁通关系ϕp=Tϕ,将式(16)整理成变量为ϕp的方程,如式(18)
A(s,ϕ)T−1ϕp=TTn(UM−nT∙ϕp)/RM. (18) 进一步整理为
nTϕp=nTTA−1(s,ϕ)TTn(UM−nT∙ϕp)/RM. (19) 将6个磁极的磁通之和β作为中间变量,则β可表示为
β=6∑j=1ϕpj. (20) 将式(20)代入式(19),得到关于β的常微分方程为
nβ=nTTA−1(s,ϕ)TTn(UM−n∙β)/RM, (21) 式中:nTTA−1(s,ϕ)TTn不再是一个矩阵或向量,而是一个关于s、ϕ的变量,记为1/M (s,ϕ).
对式(21)进行整理得
˙β=UM/n−RMM(s,ϕ). (22) 依据式(15)、(22),可得I与β间的关系为
I=nM(s,ϕ)β. (23) 依据式(19)、(22)及磁通关系ϕp=Tϕ,得到ϕ与β间的关系为
ϕ=A−1(s,ϕ)TTnM(s,ϕ)nβ. (24) 通过式(22)~(24)计算出悬浮电磁铁电流以及电磁力求解所需的磁通.
根据式(5),计算各部分气隙的磁通为
ϕa,oj=ϕajRajRa,oj. (25) 基于虚功原理,各部分气隙对应的电磁力为
Fmag,oj=ϕ2a,oj2μ0Aa,oj. (26) 单个磁极电磁力为式(26)中各部分电磁力之和,即
Fmag,j=Fmag,mj+Fmag,nj+Fmag,Lj. 最后,求解半个悬浮电磁铁的电磁力为
Fmag=6∑j=1Fmag,j. 根据电磁力解析过程,EMC模型可简化为图7所示的结构框图,输入为电压及间隙,输出为电流及电磁力. 首先进行磁阻计算,并组建磁阻矩阵,而导磁材料磁阻的计算需要将磁通作为输入. 采用磁阻矩阵A进行电流及磁通计算,磁通需通过常微分方程及代数方程求解,磁通求解结果一方面用于电磁力计算,一方面反馈给导磁材料磁阻计算.
2. 模型分析及试验验证
2.1 与传统EMC模型对比
根据悬浮电磁铁及长定子的尺寸及参数(如表1所示),对本文EMC模型进行量化. 此外,将导磁材料的相对磁导率设为恒定值,搭建基于线性导磁材料的传统EMC模型. 对两个模型电磁力的计算结果进行对比分析,如图8所示.
表 1 悬浮电磁铁及长定子参数Table 1. Parameters of maglev electromagnet and long stator项点 取值 项点 取值 定子极距/mm 258.0 铁芯厚度/mm 170.0 电磁铁极距/mm 266.5 磁极匝数 300 定子齿宽度/mm 43.0 额定磁间隙/mm 12.5 定子槽宽度/mm 43.0 恒定相对磁导率 7 000 磁间隙为12.5 mm,电流在0~35 A内,两个模型的电磁力结果非常接近;当电流超过35 A时,随着电流增加,计算结果偏差增大. 原因是实际工作状态下,随着电流增大,导磁部件的磁密增大;当达到材料饱和磁密时,磁密随电流的增加率大幅降低,电磁力也相应地出现饱和现象,而传统EMC模型并未考虑材料的磁饱和. 因此,传统EMC模型适用于小电流区间,一般应用在工作点附近的线性区间20~30 A. 而大电流区间与实际工作情况不符,例如在故障、起浮、降落等特殊工况时,模型精度大幅降低,无法用于电磁力计算及系统特性分析.
2.2 与有限元模型对比
为验证本文提出的EMC模型准确性,搭建了悬浮电磁铁与长定子FEM模型,如图9所示,两者电磁力的计算结果如图10所示. 磁间隙12.5 mm与16.0 mm,电流0~80 A内,两者电磁力的计算结果具备非常高的一致性,均存在饱和现象;磁间隙12.5 mm,电流50 A时,电磁力偏差最大,EMC计算结果为115 kN,FEM计算结果为110 kN,偏差仅为4.5%,这表明了基于非线性材料的EMC模型具有较高的准确性.
2.3 试验验证
针对高速磁浮电磁铁特性研究,搭建了地面试验平台,对悬浮电磁铁的静态电磁力进行测试,如图11所示. 试验台通过液压系统调整长定子与悬浮电磁铁的间隙;采用两路电源供电,但受最大输出电流限制,仅对0~50 A电流进行测试,步长为5 A;通过力传感器检测电磁铁与长定子间的电磁力. 磁间隙12.5 mm下电磁力的测试结果与EMC及FEM的计算结果对比如图12所示.
额定工况下,悬浮电磁铁的工作点:磁间隙为12.5 mm,电流为25 A,电磁力约为46 kN. 在工作点处,EMC、FEM及试验测试的电磁力结果几乎相同,在其他电流值下,电磁力结果偏差也极小,从而进一步验证了本文EMC模型以及所搭建FEM模型的准确性.
3. 结 论
本文基于非线性材料搭建了高速磁浮悬浮电磁铁的磁路、电路及磁力模型,将计算结果与传统EMC模型进行对比分析,并通过有限元及试验验证,通过对模型研究分析,得到以下结论:
1) 搭建悬浮电磁铁EMC模型时,采用了非线性导磁材料,通过引入B-H曲线的拟合函数,将导磁材料的非线性及饱和特性体现在模型中.
2) 无论小电流区,还是大电流区,本文EMC模型求解的电磁力均与实际情况接近,相比传统EMC模型,结果更加准确,适用范围更广.
3) 本文EMC模型能够快速、准确地求解电磁力,且通过电路模型实现与控制模型的良好匹配,因此,可通过联合仿真对悬浮系统动态特性进行深入分析,为悬浮系统设计及参数优化提供了依据.
-
表 1 单电磁铁模型参数
Table 1. Single electromagnet model parameters
参数 符号 取值 电磁铁等效荷载/kg m 6.5 电磁铁磁极面积/m2 A 0.003 75 电磁铁线圈数/匝 N 500 真空磁导率/(H•m–1) μ 4 π × 10–7 初始悬浮气隙/m z0 0.006 5 表 2 仿真参数
Table 2. Simulation parameters
参数 符号 取值 直流母线电压/V UDC 48 开关频率/kHz f 20 线圈等效电阻/Ω R0 2 线圈等效电感/mH L0 90.62 -
朱开锋,董金文,雷学国,等. 电流最小方差控制在电磁吸力悬浮系统中的应用[J]. 电气开关,2010,48(1): 58-60 doi: 10.3969/j.issn.1004-289X.2010.01.020ZHU Kaifeng, DONG Jinwen, LEI Xueguo, et al. Application of the current minimum variance control in EMS system[J]. Electric Switchgear, 2010, 48(1): 58-60 doi: 10.3969/j.issn.1004-289X.2010.01.020 张德魁,赵雷,赵鸿宾. 电流响应速度及力响应速度对磁轴承系统性能的影响[J]. 清华大学学报(自然科学版),2001,41(6): 23-26 doi: 10.3321/j.issn:1000-0054.2001.06.007ZHANG Dekui, ZHAO Lei, ZHAO Hongbin. Effect of current response rate and force response rate on performance of magnetic bearing systems[J]. Jourrnal of Tsinghua University (Science and Technology), 2001, 41(6): 23-26 doi: 10.3321/j.issn:1000-0054.2001.06.007 蒋启龙,胡基士. 磁浮列车斩波器研究[J]. 电力电子技术,1997(2): 60-62JIANG Qilong, HU Jishi. Study of the maglev chopper[J]. Power Electronics, 1997(2): 60-62 李云钢,常文森. 磁浮列车悬浮系统的串级控制[J]. 自动化学报,1999,25(2): 247-251LI Yungang, CHANG Wensen. Cascade control of an EMS maglev vehicle’s levitation control system[J]. Acta Automatic Sinica, 1999, 25(2): 247-251 周媛. 三电平悬浮斩波器性能研究[J]. 电气传动,2011,41(12): 12-15 doi: 10.3969/j.issn.1001-2095.2011.12.003ZHOU Yuan. Performance analysis of a maglev chopper based on three-level control[J]. Electric Drive, 2011, 41(12): 12-15 doi: 10.3969/j.issn.1001-2095.2011.12.003 李云钢,柯朝雄,程虎. 磁浮列车悬浮控制器的电流环分析与优化设计[J]. 国防科技大学学报,2006,28(1): 94-97 doi: 10.3969/j.issn.1001-2486.2006.01.021LI Yungang, HE Chaoxiong, CHENG Hu. Analyzing and optimizing design of current-loop in the magnetic levitation controller on maglev vehicle[J]. Journal of National University of Defense Technology, 2006, 28(1): 94-97 doi: 10.3969/j.issn.1001-2486.2006.01.021 程金路,张承进. 自适应有限拍感应电机电流控制器设计[J]. 控制与决策,2005,20(10): 1097-1102 doi: 10.3321/j.issn:1001-0920.2005.10.003CHENG Jinlu, ZHANG Chenjin. Adaptive deadbeat current controller for induction machine drives[J]. Control and Decision, 2005, 20(10): 1097-1102 doi: 10.3321/j.issn:1001-0920.2005.10.003 张东升,梅雪松,郝晓红,等. 磁悬浮系统的电流控制方法[J]. 西安交通大学学报,2007,41(9): 1096-1100 doi: 10.3321/j.issn:0253-987x.2007.09.021ZHANG Dongshen, MEI Xuesong, HAO Xiaohong, et al. Current control method for magnetic suspension system[J]. Journal of Xi’an Jiaotong University, 2007, 41(9): 1096-1100 doi: 10.3321/j.issn:0253-987x.2007.09.021 SMEDLEY K M, CUK S. One-cycle control of switching converters[J]. IEEE Transactions on Power Electronics, 1995, 10(6): 625-633 doi: 10.1109/63.471281 万蕴杰,周林,张海,等. 单周控制的发展及其应用[J]. 高电压技术,2007,33(4): 163-169 doi: 10.3969/j.issn.1003-6520.2007.04.040WAN YunJie, ZHOU Lin, ZHANG Hai, et al. Development and application of one-cycle control[J]. High Voltage Engineering, 2007, 33(4): 163-169 doi: 10.3969/j.issn.1003-6520.2007.04.040 韦徵,陈新,陈杰,等. 基于D-OCC的三相PFC整流器输入电流相位滞后及闭环补偿[J]. 中国电机工程学报,2013,33(33): 42-49WEI Zheng, CHEN Xin, CHEN Jie, et al. Input current phase lag and closed loop compensation for three-phase PFC rectifiers based on one-cycle control strategy[J]. Proceedings of the CSEE, 2013, 33(33): 42-49 BARBATI M, CALUISI C, CECATI C. One-cycle controlled active rectifier for full digital implementation[C]//36th Annual Conference on IEEE Industrial Electronics Society. Glendale: [s.n.], 2010: 625-630 LOCK A S, FERNANDES D A, SILVA E R C, et al. A current based OCC technique implemented by DSP for a three-phase OCC rectifier[C]//40th Annual Conference of the IEEE Industrial Electronics Society. Dallas: IEEE, 2014: 5136-5142 凌雁波. D-OCC技术的研究和应用[D]. 南京: 南京航空航天大学, 2009 蒋真. 单相/三相逆变器的D-OCC研究[D]. 南京: 南京航空航天大学, 2011 刘程子. 适用于高速电机的混合型磁悬浮轴承设计及控制策略的研究[D]. 南京: 南京航空航天大学, 2015 PETRONE G, VITELLI M, SPAGNULOL G. Digital implementation of one cycle control in back to back converters[C]//IEEE International Symposium on Power Electronics for Distributed Generation Systems. Aalborg: IEEE, 2012: 344-348 Ul-HAQ A, CECATI C, KHALID H A. FPGA based one-cycle control of multilevel cascaded H-bridge inverter[C]//Fortieth Annual Conference of IEEE Industrial Electronics Society. Dallas: IEEE, 2015: 4628-4634 BAGAWADE S, JAIN P. Digital implementation of one-cycle controller (OCC) for AC-DC converters[C]// IEEE International Conference on Power Electronics, Drives and Energy Systems. Trivandrum: IEEE, 2017: 1-6 刘程子,邓智泉,曹鑫,等. 混合型磁悬浮轴承开关功放的单周期数字控制[J]. 中国电机工程学报,2015,35(22): 5899-5907LIU Chengzi, DENG Zhiquan, CAO Xin, et al. An one-cycle control digital control strategy for switching power amplifiers in hybrid magnetic levitation bearing systems[J]. Proceedings of the CSEE, 2015, 35(22): 5899-5907 HAO A, SHE L, LIU H, et al. Research of force feedback controlling on high speed maglev train guidance system[C]//IEEE International Conference on Control and Automation. Guangzhou: IEEE, 2007: 1381-1384 钟毅. 磁悬浮嵌入式控制系统基础理论和关键技术研究[D]. 武汉: 武汉理工大学, 2007 -