張磊,溫旭輝
(1.中國石油大學(xué)(華東)信息與控制工程學(xué)院,山東 青島266580;2.中國科學(xué)院電工研究所,北京100080)
振動和噪聲已成為衡量電機性能的主要技術(shù)指標之一。電機的噪聲來源主要有三種:電磁噪聲、機械噪聲和空氣動力學(xué)噪聲[1-3]。其中,徑向電磁噪聲是電機噪聲的主要來源,而電機定子是電磁噪聲的主要產(chǎn)生源,因此綜合考慮這兩個因素是多數(shù)文獻的研究重點[1-7]。
國外對電機的電磁噪聲的理論研究起步較早,許多學(xué)者對電磁噪聲的產(chǎn)生機理和影響電磁噪聲的各種因素進行了比較深入的研究。文獻[8]較早地指出,將感應(yīng)電機定子作為一個單環(huán)來研究電磁振動和噪聲,而忽略電機外殼影響,與實際存在較大差距,需要考慮附加部件的影響。文獻[9]指出電機內(nèi)部的磁場是產(chǎn)生機械振動和噪聲的主要原因之一,目前用數(shù)值法求解電機磁場比較容易,然而對定子施加機械力激勵較為困難,主要是計算定子內(nèi)表面所受力時存在較大的計算誤差。而文獻[10-11]針對輸入諧波對感應(yīng)電機噪聲的影響進行了較為深入的研究,指出諧波電流和諧波電壓所產(chǎn)生的噪聲級會很高。并且當諧波頻率接近定子的自然頻率,足夠產(chǎn)生明顯的噪聲。
我國對電機噪聲的研究雖然起步較晚,但在理論分析和工程實踐中取得了一定的研究成果。前浙江大學(xué)諸自強等人[1]利用電機定子表面的振動速度和表面附近空氣中的聲壓計算電機輻射的電磁聲功率。合肥工業(yè)大學(xué)黃禮文對感應(yīng)電機電磁噪聲的成因及消除進行了相關(guān)的研究[12-13],并成功應(yīng)用了二維振動理論分析噪聲。唐任遠院士及其團隊[14]借助三維有限元軟件進行了變頻器供電的盤式永磁電機振動研究和實驗對比。文獻指出:與采用正弦供電的不同,采用變頻器供電的永磁電機其振動噪聲主要分布在開關(guān)頻率及其倍數(shù)附近,文獻具有較高的參考價值。文獻[15]雖然也嘗試建立永磁電機的解析解模型,但卻忽略了定子電流的影響;而實際上,負載輕重、電流諧波與電磁振動和噪聲大小直接相關(guān),需要進行研究。
綜上所述,國內(nèi)外關(guān)于振動模型的研究,主要還是傳統(tǒng)的交流感應(yīng)電機,較少涉及到目前應(yīng)用廣泛的永磁同步電機。同時,車用永磁同步電機相比普通電機具有更加緊湊的結(jié)構(gòu),并且廣泛采用分數(shù)槽集中繞組式結(jié)構(gòu),相比而言存在更復(fù)雜的空間、時間諧波以及更復(fù)雜的振動頻譜。針對此問題,本文建立了永磁同步電機沿空間分布的徑向電磁力波解析式,從電磁力波角度研究分析永磁同步電機電磁振動成因;并分析其成分和分布規(guī)律,以及與電磁振動相關(guān)的影響因素。
當永磁電機繞組施加三相對稱電源后,由于空間和時間諧波影響,會在定子內(nèi)表面形成空間分布的交變電磁力波,而與振動相關(guān)的主要是徑向力波分量。當徑向電磁力波的頻率與定子的固有頻率較為接近時,就可能發(fā)生共振。因此研究電機徑向電磁力波模型是研究電機振動特性的一項關(guān)鍵內(nèi)容。
首先,令:p為電機極對數(shù);θ為電機機械角度;Qs為定子槽數(shù);υ為定子電樞反應(yīng)磁動勢極對數(shù)(υ=1表示磁動勢波形的最低諧波周期是一個機械周期)。
對于內(nèi)轉(zhuǎn)子表面式結(jié)構(gòu)永磁電機,若考慮因脈寬調(diào)制等因素所產(chǎn)生的時間諧波分量,則定子磁動勢表達式為

式中:fA(n,k)為定子磁動勢諧波幅值;ωe電源基波角頻率;βwn是定子磁動勢諧波初始相位。
轉(zhuǎn)子磁鋼的等效磁動勢為

式中:fRm為轉(zhuǎn)子磁鋼等效磁動勢諧波幅值;αwm是轉(zhuǎn)子磁鋼等效磁動勢諧波初始相位。
對于具有矩形槽的氣隙磁導(dǎo)為

式中:g0為氣隙磁導(dǎo)的恒定分量;gk為諧波分量。
氣隙磁密的表達式為

若忽略氣隙磁密的切向分量,則根據(jù)麥克斯韋張量公式,位置為θ處的力波為

因此,可簡化為

從上式可以看出,在定子內(nèi)表面分布的合成力波有3部分組成:一是pA(θ,t),它是由定子電流獨立產(chǎn)生的部分,與定子繞組結(jié)構(gòu)和電流大小有關(guān),而與永磁磁場無關(guān);第二部分是pAR(θ,t),它是定子電流產(chǎn)生的磁場與永磁磁場相互作用產(chǎn)生,這里稱之為耦合項,與定轉(zhuǎn)子結(jié)構(gòu)均相關(guān);第三部分是由永磁磁場的獨立作用項pR(θ,t)。
a)定子電流獨立項 pA(θ,t)產(chǎn)生的徑向電磁力波

其中:

b)耦合項pAR(θ,t)產(chǎn)生的徑向電磁力波

其中:

c)永磁磁場獨立項 pR(θ,t)產(chǎn)生的徑向電磁力波

其中:

當空間分布的徑向電磁力波周期性地作用到定子圓周時,會激勵定子產(chǎn)生周期性徑向變形,并以聲音的形式向外傳播;其振動幅值不僅受激勵大小影響,也直接與電機固有振動特性相關(guān)。因此,研究電機振動模型是研究電機振動特性的另一個重要內(nèi)容。
辣椒去籽,粉碎過篩,稱取 5.0 g,加入40 mL蒸餾水和0.2% (W/W)的纖維素酶,50 ℃反應(yīng)2 h后進行離心分離,取上清液,殘渣用40 mL丙酮浸提劑,超聲提取30 min,過濾,得到包含辣椒紅色素的濾液。
這里主要以定子電流項所產(chǎn)生徑向電磁力波為對象進行相關(guān)的研究。若考慮電流的諧波,以及磁路結(jié)構(gòu)的不對稱和繞組分布因素所造成的,則相應(yīng)的徑向電磁力波可表示為

其中:i,j是定子電流諧波的次數(shù);n1,n2是空間諧波次數(shù),并滿足 i≠j或者 n1≠n2。
在式(10)中,與時間無關(guān)而僅與空間位置有關(guān)的直流項,只會產(chǎn)生恒定力并不會引起電機振動,可以將其忽略。因此,將式(10)中與振動無關(guān)項消除,則可簡化為

式中pm是各力波的幅值。
這樣,可以得到如表1所示的電磁力波成分,并根據(jù)力波與氣隙磁導(dǎo)諧波項的關(guān)系分為3種。

表1 永磁電機負載下電樞反應(yīng)力波組成分量Table 1 Components of armature reaction force wave of PMSMs at load
可以看出,電磁力波最小階數(shù)直接取決于電樞反應(yīng)磁動勢的極對數(shù)ν(序號2,8)。因此對于ν值小(如:分數(shù)槽集中繞組8極9槽,ν=1)的電機結(jié)構(gòu),受電樞反應(yīng)所產(chǎn)生諧波影響最大。對于整數(shù)槽分布繞組結(jié)構(gòu)的電機(ν=p),力波的最小階數(shù)較大。由于鐵心動態(tài)變形的振幅與力波階數(shù)r4成反比,因此相對于分數(shù)槽電機,整數(shù)槽結(jié)構(gòu)的電磁振動問題相對不明顯。
同理,定轉(zhuǎn)子磁場相互耦合所產(chǎn)生的徑向電磁力波表達式為

這樣,可以得到如表2所示的電磁力波成分,并根據(jù)力波與氣隙磁導(dǎo)諧波項的關(guān)系也分為3種。

表2 永磁電機負載下耦合項力波組成分量Table 2 Components of coupling force wave of PMSMs at load
從表2可以看出,由于電樞反應(yīng)與永磁磁場相互作用后,除產(chǎn)生與各自磁場相關(guān)的力波階數(shù)外,也會產(chǎn)生相互影響的新階數(shù):其中,力波最低次階數(shù)是2nν次(序號2),其頻率對應(yīng)為電源頻率的偶數(shù)次2jf。永磁磁場所產(chǎn)生力波階數(shù)主要有:序號1、3、4、6、8、9、11、12,均與電機極對數(shù) p 相關(guān)。當電機極對數(shù)p≥3時,永磁磁場所產(chǎn)生力波階數(shù)較大,對于振動的作用弱于電樞反應(yīng)所產(chǎn)生的力波。
結(jié)合前面的分析,下面對兩臺電動汽車用高功率密度永磁電機進行了相關(guān)的仿真研究工作,如表3所示。兩臺電機采用不同的電磁設(shè)計方案:I型電機采用分數(shù)槽集中繞組結(jié)構(gòu),II型電機則采用整數(shù)槽分布繞組結(jié)構(gòu)。

表3 兩種高功率密度實驗樣機參數(shù)Table 3 Parameters and items of two types of prototype high power-density motors
1)I型電機有限元仿真模型
I型電機仿真模型如圖1所示。
根據(jù)分析,得出了I型電機的一些結(jié)構(gòu)模態(tài)的基本振型和共振頻率數(shù)據(jù),如表4表6所示。從中可以看出,I型電機定子的一些主要結(jié)構(gòu)的基本模態(tài)共振頻率較低,會產(chǎn)生較大的噪聲。同時結(jié)合分析,I型電機完整定子在4~8 kHz頻率段,存在較多共振點,而這一頻率段是開關(guān)頻率諧波的范圍,因此會產(chǎn)生較為明顯的噪聲。

圖1 I型電機有限元仿真模型Fig.1 FEA simulation model of I-type prototype motor

表4 I型電機一些典型結(jié)構(gòu)處的共振頻率Table 4 Typical fundamental resonance frequency of I-type prototype motor
2)II型電機有限元仿真模型
II型電機仿真模型如圖2所示。

圖2 II型電機有限元仿真模型Fig.2 FEA simulation model of II-type prototype motor
從表5中可以看出,對于定子鐵心而言,繞組相當于附加質(zhì)量,對于其各模態(tài)頻率具有較大的影響,使得各振型的頻率均明顯增加,尤其是2階振型以上的模態(tài)增加更為明顯。
對比鐵心繞組與水套機殼的模態(tài)和振型,可以看出II型電機的機殼相同模態(tài)的共振頻率更低,在現(xiàn)有的驅(qū)動條件下更易產(chǎn)生共振。但由于采用整數(shù)槽分布繞組結(jié)構(gòu),最小力波階數(shù)較高,反而不如I型電機的振動明顯。結(jié)合表6可以發(fā)現(xiàn),在4~7 kHz這一段頻率段上,不存在明顯的共振點,電機在此頻率范圍內(nèi)振動和噪聲較小。

表5 II型電機一些典型結(jié)構(gòu)處的共振頻率Table 5 Typical fundamental resonance frequency of II-type prototype motor

表6 I和II型電機完整定子的共振頻率Table 6 Full stator resonance frequency of I and II-type motor
為了進行有效的對比研究,對I型和II型電機進行了噪聲和振動的測量分析工作。傳聲器布放于電機旁一米處,加速度計布放于電機徑向、軸向和電機支架上。通過加速度計測量電機運行時的振動情況,以及聲強計測試電機所輻射的噪聲,對I型和II型電機進行了相關(guān)的測試工作。
由于I型電機的分數(shù)槽結(jié)構(gòu),所帶來的磁密諧波豐富的原因,最小振型較低,造成其在低頻和中頻段相對II型電機振動明顯。
從圖3中可以看出,I型電機振動頻譜非常豐富,證明了前面所頻譜分析結(jié)論。主要可得出如下幾點:
1)I型電機由于采用分數(shù)槽結(jié)構(gòu),在整個頻率段內(nèi),均存在大量與電頻率相關(guān)的諧波分量:如表1和表2所示,在低頻0~2 kHz主要是基波電流及其諧波所產(chǎn)生的力波所引起。階數(shù)較低的奇數(shù)次力波為(n1±n2)ν,頻率為(i±j)f;階數(shù)較低的偶數(shù)次力波為 2nν、2nν ±2kQs,其頻率為 2jf。
2)從4~9 kHz均存在明顯的共振,最高振動幅值達到130 dB以上,這一點與模態(tài)分析結(jié)論(表5:2~8次共振頻率)在4~9 kHz較為集中的共振點也是一致;
3)在載波頻率(7.5 kHz)附近,對比模態(tài)仿真數(shù)據(jù),第6,7次共振頻率在其附近。在載波相關(guān)力波的激勵下,出現(xiàn)明顯的振動(130 dB);
4)I型電機在低頻段(5 kHz以下),軸向振動(100~120 dB)要比徑向振動(80 dB)明顯。而在5 kHz以上,徑向振動則更加突出(100~130 dB)。

圖3 8 000 r/min下I型電機徑向和軸向振動測試頻譜Fig.3 Radial and axial vibration test spectrum of I-type prototype motor at 8 000 r/min
同時,圖4展示了I型電機在不同轉(zhuǎn)速下振動加速度的頻譜分布。
從圖4(a)圖可以看出:由基波頻率2倍頻所引起的振動最明顯,兩個轉(zhuǎn)速下2倍頻振動加速度均達到130 dB,其次是4倍頻處振動(120 dB)。除此之外電源頻率奇數(shù)次。與前面的理論分析結(jié)論一致。
而從圖4(b)圖中可以看出,由于脈寬調(diào)制所引起的電流諧波,與電機共振頻率接近,因此也會引起較為明顯的振動。在開關(guān)頻率附近,存在著明顯的高頻脈寬調(diào)制諧波與低頻調(diào)制波諧波作用的電磁力波(在圖中指示),明顯提高了開關(guān)頻率附近的振動強度。
在相同的臺架上對II型電機進行了類似的實驗測試工作。從圖5可以看出,II型電機在基波頻率二倍頻處(1 037 Hz)出現(xiàn)較為明顯的振動。同時,也在3.5 kHz(軸向振動)和7 kHz(徑向和軸向)附近出現(xiàn)共振情況,這一點與模態(tài)分析的共振頻率是符合的。正如前所述,在4~7 kHz振動和噪聲明顯減少。除此之外,在載波頻率附近也出現(xiàn)了較為明顯的振動,這是由于電樞反應(yīng)磁場與轉(zhuǎn)子磁場共同作用所產(chǎn)生。振動點與有限元分析模態(tài)頻率點可以較好地吻合,說明了有限元模態(tài)分析有效性。

圖5 II型電機10 370 r/min下振動和聲壓級測量數(shù)據(jù)Fig.5 Vibration and sound pressure level spectrum of II-type prototype motor at 10 370 r/min
因此,電動汽車用高功率密度永磁同步電機相比普通永磁電機具有更加緊湊的結(jié)構(gòu),其相對于傳統(tǒng)電機存在更嚴重的時間、空間諧波以及更復(fù)雜的電磁振動成分。因此,對電機系統(tǒng)振動噪聲的要求更高,減振降噪要求更加迫切。本文針對車用高功率密度永磁電機的振動特性,進行了較為深入的研究。建立了徑向電磁力波解析式,分析其由于激勵源不同所產(chǎn)生的各種諧波電磁力密和空間分布規(guī)律,并研究與電磁振動相關(guān)的影響因素。同時,對實驗樣機高速下的振動和噪聲特性,進行了實驗測量和數(shù)值分析,驗證了徑向電磁力波解析式以及模態(tài)分析結(jié)論與實驗測試數(shù)據(jù)頻譜特征的一致性。當然,研究是在忽略磁路凸極效應(yīng)的前提下開展的,與實際情況還存在一定差異,后續(xù)工作可以在此方面進行更深入研究。
[1]陳永校,諸自強,應(yīng)善成.電機噪聲的分析和控制[M].杭州:浙江大學(xué)出版社,1987:135-151.
[2]VAN DER GIET M,F(xiàn)RANCK D,ROTHE R.Fast-and-easy acoustic optimization of PMSM by means of hybrid modeling and FEM-to-measurement transfer functions[C]//XIX International Conference on Electrical Machines,September 6 - 8,2010,Rome,Italy.2010:1 -6.
[3]VIJAYRAGHAVAN P,KRISHNAN R.Noise in electric machines:a review[J].IEEE Transactions on Industry Applications,1999,35(5):1007 -1013.
[4]GIRGIS R S,VERMA S P.Resonant frequencies and vibration behavior of stators of electrical machines as affected by teeth,windings,frame and laminations[J].IEEE Transactions on Power Apparatus and Systems,1979,98(4):1446 -1455.
[5]COLBY R S,MOTTIER F M,MILLER T J E.Vibration modes and acoustic noise in a four-phase switched reluctance motor [J].IEEE Transactions on Industry Applications,1996,32(6):1357-1364.
[6]DELAERE K,HEYLEN W,BELMANS R.Comparison of induction machine stator vibration spectra induced by reluctance forces and magnetostriction[J].IEEE Transactions on Magnetics,2002,38(2):969-972.
[7]楊艷,鄧智泉,曹鑫.無軸承開關(guān)磁阻電機徑向電磁力模型[J].電機與控制學(xué)報,2009,13(3):377-382.YANG Yan,DENG Zhiquan,CAO Xin.Magnetic radial force model of bearingless switched reluctance motors[J].Electric Machines and Control,2009,13(3):377 -382.
[8]ERDELYI E.Predetermination of sound pressure levels of magnetic noise of polyphase induction motors[J].Power Apparatus and Systems,Part III Transactions of the American Institute of Electrical Engineers,1955,74(3):1269 -1280.
[10]GIRGIS R S.Method for accurate determination of resonant frequencies and vibration behaviour of stators of electrical machines[J].IEE Proceedings of Electric Power Applications,1981,128(1):1-11.
[10]LO W C,CHAN C C,ZHU Z Q.Acoustic noise radiated by PWM-controllel induction machine drives[J].IEEE Transactions on Industrial Electronics,2000,47(4):880 -889.
[11]BELMANS R J M,D'HONDT L.Analysis of the audible noise of three-phase squirrel-cage induction motors supplied by inverters[J].IEEE Transactions on Industry Applications 1987,IA -23(5):842-847.
[12]黃禮文,盛宏玉,張偉星.準確控制異步電機電磁噪聲[J].中小型電機,1994,21(6):10-14.HANG Liwen,SHENG Hongyu,ZHANG Weixing.Accuracy control of electromagnetic noise of induction motors[J].S&M E-lectric Machines,1994,21(6):10-14.
[13]黃禮文,王宗培.電動機噪聲理論和控制技術(shù)的進展[J].電工技術(shù)學(xué)報,2000,15(5):34-38.HANG Liwen,WANG Zongpei.Development of motor noise theory and noise control technique[J].Transactions of China Electrotechnical Society,2000,15(5):34 - 38.
[14]唐任遠,宋志環(huán),于慎波.變頻器供電對永磁電機振動噪聲源的影響研究[J].電機與控制學(xué)報,2011,14(3):12-17.TANG Renyuan,SONG Zhihuan,YU Shenbo.Study on source of vibration and acoustic noise of permanent magnet machines by inverter.[J].Electric Machines and Control,2011,14(3):12-17.
[15]張冉,王秀和,喬東偉.基于輔助槽的永磁電機激振力波削弱方法[J].中國電機工程學(xué)報,2010,30(18):103-108.ZHANG Ran,WANG Xiuhe,QIAO Dongwei.Reduction of exciting force wave for permanent magnet motors by teeth notching[J].Proceedings of the CSEE,2010,30(18):103 -108.