摘要
由于双凸极结构及磁饱和特性,开关磁阻电机较它型电机有较大的径向振动.通过多物理场仿真方法,开展定子壳体径向电磁振动耦合仿真.先通过二维电磁场有限元,分别在三相和六相工作模式下,计算额定转速下电机各定子齿所受径向电磁力密度的时程.再利用定子振动特征方程,进行了定子的模态分析.最后,将电磁力密度施加到相应定子齿作为外部激励,通过瞬态结构场分析获得了定子壳体的整体变形,以及关键节点加速度在时域和频域的振动特性.结果表明,在三相模式下,正方变形模态是壳体的典型振型;在六相模式下,椭圆变形是壳体振动的典型振型.正方变形时的加速度幅值要小于椭圆变形.该工作能从振动方面指导开关磁阻电机工作模式的选择和优化.
开关磁阻电机(SRM)以其结构简单、能量利用效率高、成本低等诸多优点,在许多领域得到了广泛的应用.然而,它却呈现出较高水平的噪声和振动,限制了其在噪声敏感领域的应
众所周知,当激励基频或谐波频率与SRM定子壳体的固有频率重合时,会产生较大的振动和噪声.不同结构参数对定子固有频率的影响已进行了广泛研究,包括定子形状、定子轭厚、外壁厚度、散热筋板、安装底座
第一类是建立径向振动与运行电气参数的关系模型.朱
另一类是基于多物理场的耦合振动仿真.文献[
遵循多物理场仿真分析方法,通过将仿真的电磁力密度耦合到定子齿,研究开关磁阻电机定子壳体的径向电磁振动特性.主要分析步骤为:基于12/8极开关磁阻电机,首先,利用麦克斯韦应力方法,通过电磁仿真获得SRM内部磁密和定子齿的径向电磁力密度的时空分布;然后,进行定子壳体的模态分析和谐响应分析,作为定子壳体避振的基础;最后,施加仿真的电磁力密度到对应的定子齿作为外激励,通过瞬态仿真研究定子变形和关键点的振动特性,从振动角度为SRM工作模式的选择提供依据.
SRM定子径向耦合振动多物理场有限单元仿真流程,如

图1 SRM定子径向振动多物理场分析流程
Fig.1 Radial vibration analysis process of the stator in SRM by multiphysical simulation method
我们使用一款12/8极SRM,电机基本结构和单个定子齿所受径向力如

图2 SRM电机结构和单个定子齿所受径向力
Fig.2 SRM structure and radial force of single tooth
SRM的技术指标见
磁场强度H和磁感应密度B是进行电磁场分析的基本参量.在电磁有限元中,引入矢量磁位A求解H和B,进而获取磁场力等其它参数.
瞬态磁场中,有限元模型中各点的矢量磁位可通过运动方程求
(1) |
式中,为材料磁导率;为源电流密度;为导体电导率;为运动物体的速度;为铁磁材料矫顽力.
求得矢量磁位A后,可得到磁感应密度B和磁场强度H分别为:
(2) |
(3) |
电机中的磁力线会穿过定转子铁芯和气隙,由于气隙磁阻较大,气隙的磁压降占据磁路中的绝大部分.气隙中的磁场,既有径向分量,也有切向分量,但以径向分量为主,我们通常所说的气隙磁场、气隙磁密都默认为径向分
(4) |
其中,为真空磁导率(H/m),表示径向力方向,、分别为径向和切向气隙磁密(T).
电磁分析采用Maxwell 16.0软

(a) 控制驱动电路
(a) Control drive circuit

(b) 三相绕组的电流
(b) Winding currents
图3 控制驱动电路和绕组电流(三相运行时)
Fig.3 Control drive circuit and winding currents (three phase run)
之后,将三相绕组电流添加到电机本体模型的线圈激励,电机本体及网格划分见

图4 电机本体网格
Fig.4 Finite element method (FEM) mesh of SRM body
转子取额定转速,通过电磁场仿真可得转子运转过程中,某一时刻电机内部磁密及径向气隙磁密沿圆周的分布,见

(a) 磁密分布云图
(a) Distribution of magnetic density

(b) 气隙磁密沿圆周的分布
(b) Airgap magnetic density along the circumference
图5 三相运行时,电机内部磁密及径向气隙磁密沿圆周的分布
Fig.5 Distribution of magnetic density and radial airgap magnetic density at particular moment when SRM runs by three phase

图6 三相运行时,单齿径向电磁力密度沿圆周的分布
Fig.6 Distribution of single tooth radial magnetic force density along the circumference when SRM runs by three phase
在额定转速下,提取

(a) 定子单齿平均径向电磁力密度时程
(a) Time history of single tooth mean radial magnetic force density

(b) 平均径向电磁力密度的幅值谱
(b) Amplitude spectrum of single tooth radial magnetic force density
图7 三相运行时,定子单齿平均径向电磁力密度时程
Fig.7 Time history and amplitude spectrum of single tooth mean radial magnetic force density when SRM runs by three phase
SRM六相模式运行时分析过程类似,其获得的平均径向电磁力密度直接加载到定子相应齿.
当径向电磁力密度幅值谱中的某次谐波频率,与定子壳体固有频率一致时,定子壳体会发生共振,导致较大的振动.因此,开展模态分析以确定定子的固有频率,便成为壳体瞬态振动分析的前提.
在模态分析中,首先将连续结构离散化,依据Hamilton能量原理,每个离散后弹性结构单元的振动方程
(5) |
组合离散的单元方程组,生成总体振动方程组,整理之后,标准的强迫振动方程为
(6) |
式中,M为质量矩阵;C为阻尼矩阵;K为刚度矩阵;{u(t)}为位移矢量;{F(t)}为外力矢量.计算固有频率时,系统为无阻尼自由振动.系统特征方程为
(7) |
求解该方程,可得到固有角频率以及特征向量.该特征向量即为固有角频率(或以固有频率表示)对应的定子振型.
结构分析采用Ansys 18.0软件.将实际SRM定子壳体进行适当简化,其组成和结构网格划分如

图8 定子壳体组成(定子铁心、绕组和外壁)和网格
Fig.8 Stator composition (stator core, windings and outer wall)
对定子壳体进行模态计算时,定子采用自由边界条件.定子铁芯与外壁之间为摩擦frictional约束.定子铁芯与12个线圈为接触bonded约束.模态分析后,获得定子壳体的前四阶周向模态振型,见




图9 定子的前四阶模态振型
Fig.9 The first four order modes of the stator
考虑到电磁力密度的对称分布,当SRM以四个正交线圈为一相绕组,三相模式运行时,壳体以正方变形为显著模态.当以相对的两个线圈为一相绕组,六相模式运行时,壳体以椭圆变形为显著模态.选择电机运行转速时,应尽量避免电磁力密度的某次谐波接近定子的固有频率.
谐响应分析是确定在已知频率正弦载荷作用下,结构响应的技术.输入是已知大小和频率的谐波载荷(如:力、力密度和强迫位移),输出是关注位置点的位移、速度或加速度等.目的是探测共振响应频率,避免共振发生.在0~8,000Hz的范围内,分别在两种激励模式下,对SRM电机的一相定子齿,施加幅值恒定力密度谐波载荷,提取定子外壁中间节点的加速度幅值,绘制谐响应曲线如


(a) 三相激励及六相模式下的外部载荷
(a) External force load by three phase and six-phase excitation

(b) 三相激励下,定子外壁中间节点的谐响应曲线
(b) Harmonic response of outerwall middle node by three phase excitation

(c) 六相激励下定子外壁中间节点的谐响应曲线
(c) Harmonic response of outerwall middle node by six phase excitation
图10 定子外壁中间节点的谐响应曲线
Fig.10 Harmonic response of outerwall middle node
可知,在三相激励下,由于四个定子齿的力密度载荷呈π/2对称分布,定子壳体易发生正方变形,即在4阶固有频率处,加速度幅值出现极大值.而在六相模式激励下,两个定子齿力密度载荷呈π分布,定子壳体容易发生椭圆变形,即在2阶固有频率处,加速度幅值出现极大值.
瞬态振动分析也叫时程响应分析,是用来确定结构在随时间变化载荷作用下振动响应的方法.为使瞬态分析准确,结构瞬态响应的时间步长需小于电磁力密度的时间步长.这里,利用电磁有限元仿真出的沿着轴向均匀分布的径向电磁力密度,作为定子壳体的外激励.
以均布载荷的方式,分别在A、B、C三相定子齿(每相4个定子齿,呈π/2的角度分布),依次施加仿真的径向电磁力密度.三相运行时,各相径向电磁力密度的加载位置及时程曲线,分别见


图11 径向电磁力密度的加载位置与时程(三相激励)
Fig.11 Loading position and time history of radial magnetic force density (three phase excitation)
在径向电磁力密度的激励作用下,得到定子的整体变形、A相定子齿正对外壁中心点的加速度时程及对应的幅值谱,如

(a) 某一时刻定子的整体变形云图(三相激励)
(a) Whole deformation of the stator at particular moment (three phase)

(b) 加速度的时间历程
(b) Time history of radial acceleration

(c) 加速度的幅值谱
(c) Amplitude spectrum of radial acceleration
图12 整体变形、外壁中间节点径向振动的时程和幅值谱
Fig.12 Whole deformation, time history and amplitude spectrum of radial acceleration in outerwall intermediate node (three phase excitation)
由于SRM的四个对称分布的线圈同时受到电流激励,相应电磁力密度同时作用在对称的四个定子齿,定子易于发生正方变形,因此,电机以三相模式运行时,4阶正方变形是定子的典型振动模态.
以均布载荷的方式,分别在A、B、C、E、F、G六相定子齿(每相由正对的2个定子齿组成),施加径向电磁力密度.电机以六相运行时,径向电磁力密度的加载位置及时程曲线过程,见


图13 径向电磁力密度的加载位置与时程(六相激励)
Fig.13 Loading position and time history of radial magnetic force density (six phase excitation)
同样,在径向电磁力密度的激励作用下,得到定子的整体变形、A相定子齿正对外壁中心点的加速度时程及幅值谱,见如

(a) 某一时刻定子的整体变形云图(六相激励)
(a) Whole deformation of the stator at particular moment (six phase)

(b) 加速度的时间历程
(b) Time history of radial acceleration

(c) 加速度的幅值谱
(c) Amplitude spectrum of radial acceleration
图14 整体变形、中间节点径向振动的时程和幅值谱(六相激励)
Fig.14 Whole deformation, time history and amplitude spectrum of radial acceleration in outerwall intermediate node (six phase excitation)
由于两个对称分布的线圈受到同时电流激励,相应电磁力密度作用在相对的两个定子齿,定子壳体易于发生椭圆变形,因此,电机以六相模式运行时,2阶椭圆形模态是定子的典型模态.
经由以上分析,电机运转时电磁力密度具有对称分布特性,通过定子结构的振动传递函数之后,会激起相应定子壳体的典型模态.相较于六相激励的椭圆变形振动,三相激励的正方变形振动加速度峰值降低.因此,12/8极电机通常以三相模式运行,这有利于显著降低电机的振动及噪声.
在文献[
然而,壳体的振动会引起气隙磁密的重新分布,进而改变定子齿受到的电磁力,反过来影响壳体的径向振动,因此,考虑电磁和结构双向耦合的壳体径向振动研究,仍有待深入研究.
本文主要通过多物理场仿真方法,研究了开关磁阻电机定子壳体的径向振动特征.先建立了电磁分析模型,仿真了径向电磁力密度,并进行了分析;然后,对定子结构进行了模态分析和谐响应分析,得到典型振动模态;最后在绕组采用三相或六相激励下,分别获得了定子壳体的整体变形、外壁中心节点振动加速度的时频域特性.
外壁中心点振动信号的频谱分布,主要是转速频率与转子齿数乘积为基频的倍频成分.SRM以三相模式运行时,四阶正方变形为典型模态.六相模式运行时,二阶椭圆变形为典型模态.
该多物理场仿真方法,可快速定性获得定子齿上的径向电磁力密度分布特性,定子的固有频率及模态.通过将其施加到相应的定子齿,能快速获得定子壳体的整体变形、外壁中间节点加速度的时频域特征,奠定了壳体辐射噪声的研究基础.该方法适合在设计阶段预估开关磁阻电机的径向振动特性,为降噪设计提供技术支撑.
参考文献
Cameron D E, Lang J H, Umans S D. The origin and reduction of acoustic noise in doubly salient variable-reluctance motors. IEEE Transactions on Industry Applications, 1992, 28(6):1250~1255 [百度学术]
Wu C Y, Pollock C. Analysis and reduction of vibration and acoustic noise in the switched reluctance drive. IEEE Transactions on Industry Applications,1995,(31):91~98 [百度学术]
Takiguchi M, Sugimoto H, Kurihara N, et al. Acoustic noise and vibration reduction of SRM by elimination of third harmonic component in sum of radial forces. IEEE Transactions on Energy Conversion, 2015, 30(3):883~891 [百度学术]
Anwar M N, Husain I. Radial force calculation and acoustic noise prediction in switched reluctance machines. IEEE Transactions on Industry Applications, 1999, 36(6):1589~1597 [百度学术]
Wang H, Bao J, Xue B, et al. Control of suspending force in novel Permanent-magnet-biased bearingless switched reluctance motor. IEEE Transactions on Industrial Electronics, 2015, 62(7):4298~4306 [百度学术]
Guo X Q , Zhong R , Ding D S. Origin of resonance noise and analysis of randomising turn-on angle method in switched reluctance motor. IET Electric Power Applications, 2017, 11(7):1324~1332 [百度学术]
Pillay P, Cai W. An investigation into vibration in switched reluctance motors. IEEE Transactions on Industry Applications, 1999, 35(3):589~596 [百度学术]
Cai W, Pillay P. Resonant frequencies and mode shapes of switched reluctance motors. IEEE Transactions on Energy Conversion, 2001, 16(1):43~48 [百度学术]
Castano S M, Bilgin B, Fairall E, et al. Acoustic noise analysis of a high-speed high-power switched reluctance machine: frame effects. IEEE Transactions on Energy Conversion, 2016, 31(1):69~77 [百度学术]
Li J, Song X G,Cho Y. Comparison of 12/8 and 6 /4 switched reluctance motor: noise and vibration aspects. IEEE Transactions on Magnetics, 2008, 44(11):4131~4134 [百度学术]
Kurihara N, Bayless J, Sugimoto H, et al.Noise reduction of switched reluctance motor with high number of poles by novel simplified current waveform at low speed and low torque region. IEEE Transactions on Industry Applications, 2016, 52(4):3013~3021 [百度学术]
Liang X, Li G, Ojeda J, et al. Comparative study of classical and mutually coupled switched reluctance motors using multiphysics finite-element modeling. IEEE Transactions on Industrial Electronics, 2014, 61(9):5066~5074 [百度学术]
Ahn J W, Park S J, Lee D H. Hybrid excitation of SRM for reduction of vibration and acoustic noise. IEEE Transactions on Industrial Electronics, 2004, 51(2):374~380 [百度学术]
Dong J N, Jiang J W, Howey B, et al. Hybrid acoustic noise analysis approach of conventional and mutually coupled switched reluctance motors. IEEE Transactions on Energy Conversion, 2017,32(3):1041~1051 [百度学术]
Yang H Y, Lim Y C, Kim H C. Acoustic noise/vibration reduction of a single-phase srm using skewed stator and rotor. IEEE Transactions on Industrial Electronics, 2013, 60(10):4292~4300 [百度学术]
Gan C, Wu J, Shen M, et al. Investigation of skewing effects on the vibration reduction of three-phase switched reluctance motors. IEEE Transactions on Magnetics, 2015, 51(9): 1~9 [百度学术]
Isfahani A H, Fahimi B.Comparison of mechanical vibration between a double-stator switched reluctance machine and a conventional switched reluctance machine. IEEE Transactions on Magnetics, 2014, 50(2):293~296 [百度学术]
Zhu Z Q, Liu X, Pan Z.Analytical model for predicting maximum reduction levels of vibration and noise in switched reluctance machine by active vibration cancellation. IEEE Transactions on Energy Conversion, 2011, 26(1): 36~45 [百度学术]
Guo X, Zhong R, Zhao L, et al.Method for radial vibration modelling in switched reluctance moto.IET Electric Power Applications, 2016, 10(9): 834~842 [百度学术]
Lin C, Fahimi B. Prediction of radial vibration in switched reluctance machines. IEEE Transactions on Energy Conversion, 2013, 28(4): 1072~1081 [百度学术]
Guo X, Zhong R, Zhang M, et al.Fast computation of radial vibration in switched reluctance motors. IEEE Transactions on Industrial Electronics, 2018,65(6): 4588~4597 [百度学术]
Tang Z, Pillay P, Omekanda A M. Vibration prediction in switched reluctance motors with transfer function identification from shaker and force hammer tests. IEEE Transactions on Industry Applications, 2003,39(4): 978~985 [百度学术]
Torregrossa D,Fahimi B, Peyraut F,et al. Fast computation of electromagnetic vibrations in electrical machines via field reconstruction method and knowledge of mechanical impulse response. IEEE Transactions on Industrial Electronics, 2012, 59(2): 839~847 [百度学术]
Boesing M, Hofmann A, De Doncker R.Universal acoustic modelling framework for electrical drives. IET Power Electronics, 2015, 8(5): 693~699 [百度学术]
Li J, Choi D, Cho Y.Analysis of rotor eccentricity in switched reluctance motor with parallel winding using FEM. IEEE Transactions on Magnetics, 2009, 45(6):2851~2854 [百度学术]
张鑫, 王秀和, 杨玉波. 基于转子齿两侧开槽的开关磁阻电机振动抑制方法研究. 中国电机工程学报, 2015, 35(6):1508~1515 [百度学术]
Zhang X,Wang X H,Yang Y B.Vibration reduction of a switched reluctance motor using new rotor tooth with slot on each side. Proceedings of the CSEE, 2015, 35(6):1508~1515 (in Chinese) [百度学术]
刘慧娟. Ansoft Maxwell 13电机电磁场实例分析. 北京:国防工业出版社, 2014 (Liu H J. Ansoft Maxwell 13 motor electromagnetic field example analysis. BeiJing:National Defense Industry Press, 2014 (in Chinese)) [百度学术]