本文引入了一種分?jǐn)?shù)槽PMSM振動頻譜特性及振動源的分析方法，以一臺6槽4極非均勻氣隙內(nèi)置式PMSM為研究對象，理論分析了樣機(jī)電磁力波特征參數(shù)和諧波來源；并建立了樣機(jī)多物理場磁-固耦合仿真模型，計算分析了樣機(jī)的電磁振動頻譜特性；最后，對樣機(jī)的振動頻譜特性進(jìn)行了實(shí)驗驗證。

1

電機(jī)基本參數(shù)及模型

常規(guī)均勻氣隙的內(nèi)置式PMSM，轉(zhuǎn)子結(jié)構(gòu)的形狀為以點(diǎn)O為圓心的圓形。為了改善電機(jī)氣隙磁場磁密波形，在永磁體一個極距范圍內(nèi)，以O(shè)′為轉(zhuǎn)子外圓的圓心，結(jié)構(gòu)如圖1所示，用偏心距dm=OO′的長度來衡量電機(jī)氣隙非均勻程度。

圖1內(nèi)置式PMSM非均勻氣隙結(jié)構(gòu)圖

樣機(jī)主要參數(shù)如表1所示，根據(jù)樣機(jī)參數(shù)建立的二維電磁計算模型如圖2所示。

表1樣機(jī)主要參數(shù)

圖2樣機(jī)基本模型

2

電機(jī)電磁力波的分析

2.1電磁力波特性分析

在不考慮飽和與鐵心磁阻的影響情況下，內(nèi)置式PMSM氣隙磁場正弦電流供電時，氣隙磁密的解析表達(dá)式：

b(θ,t)=f(θ,t)Λ(θ)=[∑fμ(θ,t)+∑fv(θ,t)]·

(1)

式中：∑fμ(θ,t)為轉(zhuǎn)子永磁體諧波磁動勢；∑fv(θ,t)為定子繞組諧波磁動勢；Λ(θ)為相對磁導(dǎo)函數(shù)；μ為轉(zhuǎn)子磁場諧波次數(shù)；v為定子磁場諧波次數(shù)；ω為基波磁勢角頻率；p為電機(jī)極對數(shù)；z為電機(jī)槽數(shù)。

去除切向磁通密度的影響，由麥克斯韋應(yīng)力張量方程可知，徑向電磁力波的瞬時值在單位面積上解析表達(dá)式[19]：

ppm+ppm-s+ps

(2)

式中：b(θ,t)為徑向氣隙磁密；μ0為真空磁導(dǎo)率，μ0=4π×10-7H/m；ppm為轉(zhuǎn)子磁場之間互相作用引起的電磁力密度；ps為定子磁場之間互相作用引起的電磁力密度；ppm-s為轉(zhuǎn)子磁場與定子磁場互相作用引起的電磁力密度。

PMSM徑向力波的空間力波階數(shù)r與電機(jī)極對數(shù)和槽數(shù)有關(guān)。電磁力波引起的振動與力波幅值呈正比，與力波階數(shù)r的4次方呈反比，因此在計算電機(jī)的振動時，可以忽略高階力波的影響[20]。根據(jù)前文理論分析結(jié)果，可以得到電磁力波在轉(zhuǎn)子磁場之間互相作用、轉(zhuǎn)子磁場與定子磁場互相作用兩種情況下的各頻譜分量和諧波來源，如表2和表3所示。

表26槽4極電機(jī)轉(zhuǎn)子磁場之間互相作用電磁力波來源表

表36槽4極電機(jī)轉(zhuǎn)子磁場與定子磁場互相作用電磁力波來源表

表2、表3中為定子基波頻率，n為電機(jī)轉(zhuǎn)速，p為電機(jī)極對數(shù)，z為電機(jī)槽數(shù)。

根據(jù)表2和表3可知，r=2、4階為6槽4極PMSM的空間徑向電磁力波的主要階數(shù)，其中r=2階電磁力波的力波頻率主要有4f，10f，16f，22f，…；r=4階電磁力波的力波頻率主要有2f，8f，14f，20f，…。

分析表2可得，轉(zhuǎn)子基波磁場和轉(zhuǎn)子磁場3次、5次磁密諧波是產(chǎn)生樣機(jī)轉(zhuǎn)子磁場互相作用r=2階電磁力波的4f與r=4階電磁力波的2f的主要諧波次數(shù)。

分析表3可得，轉(zhuǎn)子磁場3次、5次磁密諧波和定子基波磁場、2次、4次磁密諧波是產(chǎn)生樣機(jī)電樞反應(yīng)r=2階電磁力波的4f的主要諧波次數(shù)；轉(zhuǎn)子基波磁場、3次磁密諧波和定子基波磁場、2次、4次磁密諧波是產(chǎn)生樣機(jī)電樞反應(yīng)r=4階電磁力波的2f的主要諧波次數(shù)。

2.2氣隙磁場有限元分析

利用有限元軟件建立6槽4極非均勻氣隙PMSM二維電磁模型，仿真計算出樣機(jī)在額定轉(zhuǎn)速下的空載和負(fù)載氣隙磁通密度，并在某一時間下對其進(jìn)行傅里葉級數(shù)分解，得出氣隙徑向磁通密度及其傅里葉級數(shù)分解，如圖3所示。

(a) 徑向磁通密度

(b)徑向磁通密度傅里葉級數(shù)分解
圖3某一時刻氣隙徑向磁通密度及其傅里葉級數(shù)分解

由圖3可知，徑向磁通密度包含大量的諧波，次數(shù)較大的有3次、5次。負(fù)載工作點(diǎn)相對于空載工作點(diǎn)，各次諧波幅值均增大，諧波之間相互作用可能產(chǎn)生較大的電磁力波，和空載時相比會引起更大的振動噪聲。

3

多物理場電磁振動仿真計算

圖4為本文PMSM電磁場和結(jié)構(gòu)場耦合電磁振動仿真計算流程圖。首先建立用于電磁力計算的電磁有限元模型，然后建立電機(jī)定子的結(jié)構(gòu)有限元模型，把二維電磁場計算得到的電磁力加載到三維電機(jī)定子結(jié)構(gòu)上,計算定子振動的響應(yīng)頻譜特性，最后測試樣機(jī)的振動特性，驗證有限元仿真結(jié)果。

圖4電機(jī)電磁振動仿真計算流程圖

3.1電磁力波仿真分析

采用瞬態(tài)有限元方法對樣機(jī)進(jìn)行分析計算，得到氣隙磁場與定子鐵心相互作用產(chǎn)生的電磁力密度，如圖5所示。仿真計算出樣機(jī)在額定轉(zhuǎn)速下的空載和負(fù)載徑向電磁力，并對其進(jìn)行空間上的傅里葉級數(shù)分解，結(jié)果如圖6所示。

圖5定子齒徑向力密度

由圖6可知，樣機(jī)的空間徑向電磁力波的階數(shù)主要是r=2階、r=4階、r=6階等，造成電機(jī)電磁振動的主要原因為低階電磁力波的存在，對于6槽4極PMSM，r=2階徑向電磁力波是造成電機(jī)電磁振動的主要原因。負(fù)載時各階電磁力波的幅值都增加，r=2階的幅值增加了22%，因此電機(jī)負(fù)載時的振動噪聲加劇。電機(jī)運(yùn)行在負(fù)載情況下，會發(fā)生電樞反應(yīng)，造成電機(jī)相對于空載運(yùn)行時的電磁力波幅值增大。該結(jié)果與電磁力波解析法分析結(jié)果一致。

(a) 徑向電磁力

(b) 徑向電磁力傅里葉級數(shù)分解
圖6空間上徑向電磁力及其傅里葉級數(shù)分解

圖7為電機(jī)運(yùn)行在額定轉(zhuǎn)速下，電機(jī)內(nèi)部齒中心處的一點(diǎn)在時間變化上徑向電磁力的波形及其傅里葉級數(shù)分解結(jié)果。

(a) 徑向電磁力

(b)徑向電磁力傅里葉級數(shù)分解
圖7時間上徑向電磁力及其傅里葉級數(shù)分解

由圖7可知，徑向電磁力波的時間頻率是2f，4f，6f，…，其中2f分量的幅值遠(yuǎn)遠(yuǎn)大于其它諧波分量。由表2和表3可知，電磁力波的2f分量主要由基波產(chǎn)生，因此幅值比較大；和空載時相比，負(fù)載時各諧波分量的幅值增加。這是因為空載時轉(zhuǎn)子磁場相互作用產(chǎn)生徑向力波諧波分量，負(fù)載時電磁力諧波除了轉(zhuǎn)子磁場相互作用產(chǎn)生外，定轉(zhuǎn)子磁場相互作用也會產(chǎn)生大量的電磁力諧波，有限元結(jié)果驗證了前文的理論分析結(jié)果。

3.2樣機(jī)模態(tài)分析

對電機(jī)進(jìn)行固有模態(tài)特性分析，設(shè)計電機(jī)時盡量使電機(jī)定子的固有頻率避開徑向電磁力波的高頻頻率，避免發(fā)生結(jié)構(gòu)共振，是抑制電機(jī)電磁振動噪聲的關(guān)鍵。對電機(jī)定子鐵心進(jìn)行模態(tài)分析，定子鐵心密度為7 410 kg/m3，彈性模量為206 GPa，泊松比為0.3，定子鐵心三維有限元仿真模型與網(wǎng)格剖分圖如圖8所示。

(a) 三維有限元模型

(b)網(wǎng)格剖分圖
圖8定子鐵心三維有限元仿真模型與網(wǎng)格剖分圖

圖9為樣機(jī)定子鐵心模態(tài)振型與模態(tài)頻率仿真結(jié)果。

(a)n=2f=2 015 Hz

(b)n=3f=4 611 Hz

(c)n=4f=8 930 Hz

(d)n=5f=11 690 Hz
圖9樣機(jī)定子結(jié)構(gòu)有限元振動模態(tài)

3.3多物理場電磁振動頻譜分析

建立電磁場和結(jié)構(gòu)場的耦合模型。將電磁力施加到電機(jī)定子齒內(nèi)表面上，得到電機(jī)的電磁振動響應(yīng)頻譜特性，如圖10所示。

(a)有限元仿真響應(yīng)頻譜

(b) 2 020 Hz處振型

(c)8 940 Hz處振型
圖10樣機(jī)電磁振動有限元響應(yīng)頻譜及形變觀測圖

由圖10可知，定子鐵心在2f=340 Hz，4f=700 Hz，2 020 Hz，8 940 Hz處均發(fā)生了明顯的振動。

由前文分析得到的徑向電磁力成分可知，頻率為2f，4f的徑向電磁力幅值大，在這些頻率下電機(jī)振動明顯。

頻率在2 020 Hz，8 940 Hz處對應(yīng)的振動加速度幅值相對較大，分別與電機(jī)2015Hz(n=2)、8 930 Hz(n=4)的固有頻率接近，振動加速度的幅值分別為3.559 m/s2、0.867 m/s2，其中，二階固有頻率處的振動加速度最大，但由于2 020 Hz=12f、8 940 Hz=52f，該兩頻率下的電磁力幅值相對較小，不會引起較大的共振。

綜上，徑向電磁力幅值較大時所對應(yīng)的頻率，遠(yuǎn)離電機(jī)各階模態(tài)所對應(yīng)的固有頻率，本文樣機(jī)的固有頻率避開了電磁力較大的力波階數(shù)的頻譜，有限元仿真結(jié)果表明，該電機(jī)振動幅值較小，不會輻射較大的振動噪聲。

4

樣機(jī)振動實(shí)驗測試

為了驗證前文的分析結(jié)果，實(shí)驗測試樣機(jī)的振動特性。實(shí)驗測試圖如圖11所示，其中，使用的動態(tài)數(shù)據(jù)采集分析儀型號為AVANT MI-7008，最大采樣頻率為204.8 kHz；ICP型三向加速度傳感器采集信號，靈敏度100 mV/g；負(fù)載為磁粉制動器。

圖11樣機(jī)振動實(shí)測圖

圖12為額定工況下電機(jī)運(yùn)行在5 200 r/min時徑向振動加速度與實(shí)測結(jié)果的對比。從圖12中可看出，仿真計算與實(shí)驗測量在2 000 Hz和8 900 Hz附近有較大的振動，驗證了仿真結(jié)果的正確性。實(shí)測振動響應(yīng)頻譜波形波動較大，這是由于電機(jī)振動實(shí)驗中是整機(jī)，而仿真計算僅研究了樣機(jī)的定子鐵心結(jié)構(gòu)，因此實(shí)驗結(jié)果與計算結(jié)果有誤差。

圖12樣機(jī)仿真計算電磁振動與實(shí)測結(jié)果對比

5

結(jié) 語

本文引入了一種分?jǐn)?shù)槽PMSM振動頻譜特性和振動源的分析方法，理論分析了PMSM電磁力波的產(chǎn)生原理和各階力波的諧波來源，對一臺6槽4極非均勻氣隙PMSM建立了電磁場和結(jié)構(gòu)場多物理場耦合有限元模型，分析了電機(jī)電磁振動特性，并進(jìn)行了實(shí)驗驗證。仿真與實(shí)驗結(jié)果表明，對于本文6槽4極樣機(jī)，樣機(jī)各階固有頻率遠(yuǎn)離引起電機(jī)共振的力波頻率，運(yùn)行過程中不會發(fā)生較大的共振。本文方法也適用于其他極槽配合電機(jī)的振動頻譜特性和振動源的分析。

審核編輯 ：李倩