摘 要：

在電機運行過程中，氣隙偏心是一種常見的故障狀態(tài)，會危害到電機的安全可靠運行，而該問題在大容量高速永磁同步電機中更為嚴(yán)重。因此，本文以一臺兆瓦級高速永磁同步電動機為例，研究了氣隙偏心對電機電磁性能以及其引發(fā)的不平衡磁拉力對轉(zhuǎn)子振動特性的影響。首先，推導(dǎo)了靜、動以及混合偏心狀態(tài)下的不平衡磁拉力表達式。然后，建立了各偏心狀態(tài)下的電機有限元模型，定量分析了不同偏心程度所引發(fā)的不平衡磁拉力變化情況。在此基礎(chǔ)上，計算了電機轉(zhuǎn)子系統(tǒng)的動力響應(yīng)，詳細(xì)討論了各種偏心對轉(zhuǎn)子系統(tǒng)軸心軌跡和電磁性能的影響。結(jié)果表明，不同偏心狀態(tài)下不平衡磁拉力存在明顯差異，在混合偏心時動偏心會削弱靜偏心對轉(zhuǎn)子振動的影響，偏心將增加氣隙磁密的偶次諧波且增大轉(zhuǎn)矩脈動。

關(guān)鍵詞：高速永磁電機；氣隙偏心;不平衡磁拉力；振動特性；電磁性能；軸心軌跡

DOI：10.15938/j.emc.2024.09.008

中圖分類號：TM355

文獻標(biāo)志碼：A

文章編號：1007-449X（2024）09-0081-14

收稿日期： 2023-08-21

基金項目：國家自然科學(xué)基金區(qū)域創(chuàng)新發(fā)展聯(lián)合基金重點支持項目（U21A20145）；國家自然科學(xué)基金（51977052）；黑龍江省自然科學(xué)基金重點項目（ZD2022E006）

作者簡介：謝 穎（1974—），女，博士，教授，博士生導(dǎo)師，研究方向為電機內(nèi)綜合物理場計算、新能源汽車用電機設(shè)計及多目標(biāo)優(yōu)化；

楊艷會（1996—），女，碩士研究生，研究方向為高速永磁同步電機的設(shè)計及振動特性分析；

蔡 蔚（1959—），男，博士，教授，博士生導(dǎo)師，研究方向為驅(qū)動電機、功率電子控制器及汽車電動化電驅(qū)動系統(tǒng)等；

李道璐（1998—），女，博士研究生，研究方向為永磁同步電機設(shè)計及振動噪聲計算分析；

王澤兵（1999—），男，碩士研究生，研究方向為高速永磁同步電機的設(shè)計與優(yōu)化；

張 燚（1999—），男，碩士研究生，研究方向為磁場調(diào)制式復(fù)合電機。

通信作者：謝 穎

Analysis of rotor vibration characteristics of megawatt high speed permanent magnet synchronous motor under eccentric state

XIE Ying， YANG Yanhui， CAI Wei， LI Daolu， WANG Zebing， ZHANG Yi

（School of Electrical and Electronic Engineering， Harbin University of Science and Technology， Harbin 150080， China）

Abstract：

During motor operation， air-gap eccentricity is a common fault situation， which will endanger the safe and reliable operation of the motor， and the problem is more serious in the high-power high-speed permanent magnet synchronous motor. Therefore， the effects of air-gap eccentricity on electromagnetic performance of the motor and unbalanced magnetic pull （UMP） due to air-gap eccentricity on vibration characteristics of the rotor were studied. Firstly， the formulas of unbalanced magnetic tension in static， dynamic and mixed eccentricity were deduced. Then the finite element model of the motor under each eccentricity situation were established， and the variations of the UMP caused by different eccentricities were quantitatively analyzed. On this basis， the dynamic response of the motor rotor system was calculated， and the influence of various eccentricities on the axis locus of the rotor system and electromagnetic performance of the motor were discussed in detail. The results show that there are obvious differences in UMP under different eccentricity situations. Dynamic eccentricity will reduce the influence of static eccentricity on rotor vibration during mixed eccentricity. Eccentricity will increase even harmonics of air-gap flux density and torque ripple.

Keywords：high speed permanent magnet motor; air-gap eccentricity; unbalanced magnetic pull; vibration characteristic; electromagnetic performance

0 引 言

永磁同步電機具有功率密度大，響應(yīng)速度快和效率高等顯著優(yōu)點［1］，在諸多領(lǐng)域得到了廣泛的應(yīng)用。但是由于電機制造和安裝誤差、軸承磨損等因素的影響，電機定轉(zhuǎn)子間可能會發(fā)生偏心，使得氣隙磁場發(fā)生畸變產(chǎn)生不平衡磁拉力［2］，進一步加劇電機振動，加快軸承磨損，嚴(yán)重情況下還可能引起永磁電機定子鐵心與轉(zhuǎn)子接觸，危害電機安全可靠運行。此外，偏心會引起磁場分布不均勻、磁阻變化等，進而影響電機的輸出轉(zhuǎn)矩等［3］。因此，研究偏心對高速永磁電機轉(zhuǎn)子振動及電磁性能的影響具有重要意義。

轉(zhuǎn)子振動主要受氣隙偏心所引發(fā)的不平衡磁拉力的影響［4］，因此不平衡磁拉力的研究對轉(zhuǎn)子振動至關(guān)重要。最先針對不平衡磁拉力的研究均是基于線性解析法，文獻［5］中根據(jù)假設(shè)，第一次建立了線性的不平衡磁拉力解析式，但是后來通過工程實踐發(fā)現(xiàn)線性方法誤差較大，隨后大量學(xué)者開始研究不平衡磁拉力的非線性解析式。文獻［4-7］對外轉(zhuǎn)子永磁電機的偏心氣隙磁場進行直接解析求解，得到氣隙中各處的磁密分布，然后對轉(zhuǎn)子表面的Maxwell應(yīng)力進行積分得到磁拉力解析式，但其推導(dǎo)過程及結(jié)果過于復(fù)雜，約束了該方法的應(yīng)用。文獻［8］通過將氣隙磁導(dǎo)展開為級數(shù)形式，從而得到了不同極對數(shù)下三相同步電機空載時的不平衡磁拉力，然后分析了不同極數(shù)下剛性支承轉(zhuǎn)子系統(tǒng)的振動特性。文獻［9］研究了永磁同步電機負(fù)載時的不平衡磁拉力表達，代入建立的電機轉(zhuǎn)子系統(tǒng)彎扭耦合方程，討論了靜偏心和質(zhì)量偏心以及負(fù)載類型對轉(zhuǎn)子系統(tǒng)振動特性的影響，但是沒有對動靜復(fù)合偏心進行分析。文獻［10］對動偏心和動靜復(fù)合偏心情形下不考慮重力時和考慮重力時轉(zhuǎn)子系統(tǒng)的軸心軌跡，討論了靜偏心方向?qū)D(zhuǎn)子系統(tǒng)振動響應(yīng)的影響，但是所涉及的電機是小型低轉(zhuǎn)速電機。偏心時，氣隙長度不再均勻，使得氣隙磁場發(fā)生畸變進而影響電機性能。文獻［11］以一臺低速永磁同步電機為例，建立轉(zhuǎn)子偏心的電機模型，分析了偏心對電機電磁力特性的影響。文獻［12-15］研究了轉(zhuǎn)子偏心對永磁同步電機氣隙磁場的影響，在此基礎(chǔ)上分析了氣隙諧波磁場的變化對電機轉(zhuǎn)子表面渦流電密的影響。綜上所述，目前針對轉(zhuǎn)子偏心對電機振動和電磁性能影響的研究，主要集中在小型電機和中低速電機上，對兆瓦級高速永磁電機的研究較為有限。

本文首先設(shè)計一臺1.2 MW，20 000 r/min的高速永磁同步電機，利用有限元法驗證電磁設(shè)計的合理性。然后建立各種偏心狀態(tài)下的二維電機模型，定量分析不同偏心量以及不同偏心狀態(tài)對不平衡磁拉力的影響。通過建立轉(zhuǎn)子系統(tǒng)動力模型研究其在不平衡磁拉力作用下的動力響應(yīng)，詳細(xì)討論各種偏心狀態(tài)下轉(zhuǎn)子系統(tǒng)的軸心軌跡和振動頻譜的變化。最后對偏心前后電機性能進行對比分析，得出偏心對電磁性能的影響規(guī)律。

1 電機設(shè)計與電磁性能分析

本文設(shè)計了一臺兆瓦級高速永磁同步電機，電機采用4極轉(zhuǎn)子結(jié)構(gòu)，永磁體徑向分為2段。采用碳纖維保護套綁扎轉(zhuǎn)子鐵心，來滿足轉(zhuǎn)子結(jié)構(gòu)的機械強度要求。定子采用27槽矩形槽結(jié)構(gòu)，電機結(jié)構(gòu)如圖1所示，設(shè)計參數(shù)如表1所示。

高速永磁同步電機按照額定轉(zhuǎn)速20 000 r/min旋轉(zhuǎn)的情況下，其轉(zhuǎn)動頻率為

fr=n60。（1）

對本文所設(shè)計的高速永磁同步電機進行電磁特性仿真，額定轉(zhuǎn)速運行時的空載相反電勢曲線如圖2所示?？蛰d相反電勢趨于平頂波，有效值約為1.67 kV，滿足性能要求。

額定工況下的輸出轉(zhuǎn)矩曲線如圖3所示，負(fù)載轉(zhuǎn)矩在563～584 N·m之間波動，平均值約為573 N·m，轉(zhuǎn)矩脈動約為3.6%。其中，轉(zhuǎn)矩脈動計算公式可表示為

Tripple=Tmax－TminTavg×100%。（2）

式中：Tmax為電機轉(zhuǎn)矩最大值；Tmin為電機轉(zhuǎn)矩最小值；Tavg為電機平均轉(zhuǎn)矩。

2 轉(zhuǎn)子系統(tǒng)運動方程的建立

2.1 不平衡磁拉力計算

電機的氣隙偏心主要分為靜態(tài)偏心和動態(tài)偏心。電機偏心情況如圖4所示。靜態(tài)偏心主要由定子鐵心變形、定轉(zhuǎn)子裝配誤差等因素引起，其特點是轉(zhuǎn)子的旋轉(zhuǎn)中心O與幾何中心O重合，但與定子圓心Os偏移一定距離，最小氣隙位置固定不變。動態(tài)偏心主要由轉(zhuǎn)軸受力彎曲、轉(zhuǎn)子質(zhì)量偏心等因素造成，其特點是轉(zhuǎn)子的旋轉(zhuǎn)中心Os與定子圓心Os重合，但與轉(zhuǎn)子幾何中心O偏移一定距離，最小氣隙位置隨轉(zhuǎn)子運動而改變，混合偏心是動、靜偏心同時存在的狀態(tài)。

若式（11）中的ω為0，則為靜偏心下的徑向電磁力波表達式。根據(jù)解析式推導(dǎo)出偏心狀態(tài)下徑向電磁力波的空間階次和頻率特性變化，如表2所示（其中fc為電機額定頻率，p為極對數(shù)，fr為電機額定轉(zhuǎn)速下的轉(zhuǎn)頻頻率，n∈0、1、2、3…）。

對偏心下的徑向電磁力波沿圓周積分并分解到X軸和Y軸方向上，得到兩個方向上的不平衡磁拉力表達式如下：

FX=RL∫2π0prcosθdθ=RL2μ0∫2π0Br（a，t）2cosθdθ；（12）

FY=RL∫2π0prsinθdθ=RL2μ0∫2π0Br（a，t）2sinθdθ。（13）

式中：R為轉(zhuǎn)子半徑；L為氣隙軸向長度。

2.2 轉(zhuǎn)子系統(tǒng)動力學(xué)模型

當(dāng)電機發(fā)生氣隙偏心時，由此產(chǎn)生的不平衡磁拉力將作用在轉(zhuǎn)子系統(tǒng)上，進而加劇轉(zhuǎn)子振動。為研究偏心下電機轉(zhuǎn)子的振動特性，采用集總參數(shù)法建立轉(zhuǎn)子系統(tǒng)模型見圖5所示。根據(jù)結(jié)點總數(shù)N應(yīng)滿足N≥1+5.23r的關(guān)系［18］（r為要計算的固有頻率的最高階數(shù)，因低階固有頻率對系統(tǒng)振動影響較大，故本文只考慮前三階固有頻率，選定r=3，將更高階舍棄），轉(zhuǎn)子系統(tǒng)整體被劃分為21個節(jié)點，彈性軸段為20個BEAM188梁單元，2個軸承為COMBI214彈簧單元。

若轉(zhuǎn)子激勵頻率達到或接近固有頻率時，振動將會加劇產(chǎn)生共振，為分析轉(zhuǎn)子系統(tǒng)在固有頻率處的振動情況，本文首先對轉(zhuǎn)子系統(tǒng)進行模態(tài)分析，得到前三階固有頻率及臨界轉(zhuǎn)速如表3所示。

3 偏心對電機電磁性能的影響

本文中記動偏心為ed，靜偏心為es，動、靜兩種偏心故障下的偏心量均選取0.5、1、1.5 mm分別占?xì)庀堕L度的20%、40%、60%，偏心方向選取圖4所示的X軸正方向。電機實際運行中往往存在動靜混合偏心，依據(jù)動、靜偏心量的設(shè)置，混合偏心共有9種組合，但由于電機氣隙長度為2.5 mm，因此偏心距離大于或等于氣隙長度的組合不符合實際情況，具體混合偏心組合如表4所示。偏心導(dǎo)致電機內(nèi)氣隙長度變化，磁場分布不均勻，進而影響電機輸出轉(zhuǎn)矩，因此本節(jié)主要對偏心前后的電機氣隙磁密和輸出轉(zhuǎn)矩進行分析，其中由于氣隙磁場是隨時間和空間變化的函數(shù)，故對氣隙磁場從時間和空間兩個方面展開分析。

3.1 偏心對氣隙磁場時間諧波的影響

沿偏心方向取電機氣隙中的一點記為特征點，偏心方向記為轉(zhuǎn)子偏向側(cè)，與偏心方向相反的一側(cè)記為轉(zhuǎn)子偏離側(cè)，特征點示意圖如圖6所示。通過有限元法計算得到，一個電周期內(nèi)該特征點在不同偏心狀態(tài)下，氣隙磁密隨時間的變化曲線如圖7所示，圖中0 ms時刻對應(yīng)轉(zhuǎn)子偏向側(cè)，電機轉(zhuǎn)動一個電周期即1.5 ms時轉(zhuǎn)至偏離側(cè)。由圖6可知，靜偏心時，該點的磁密強度隨偏心率的增加而增大。動偏心時，電機在轉(zhuǎn)子偏向側(cè)時該點的磁密強度隨偏心率的增加而增大，隨著電機向偏離側(cè)轉(zhuǎn)動，該點的磁密強度逐漸減小，并會低于正常狀態(tài)的磁密強度，這是因為動偏心時這一點的氣隙長度隨電機轉(zhuǎn)動而逐漸增大?；旌掀臅r由圖6（c）和圖6（d）可知，當(dāng)靜偏心率大于或等于動偏心率，即靜偏心的影響占主要因素時，其磁密強度均大于其他偏心組合的磁密強度。

對上述隨時間變化的氣隙磁密進行傅里葉分解，圖8為靜、動偏心下的氣隙磁密幅頻特性，由圖可知，靜偏心時，氣隙磁密各次諧波kfc（k=1、3、5、7，…）的幅值，隨靜偏心率的增加逐漸增大。動偏心時，主要諧波kfc的幅值變化不大，但引入了kfc±fr的頻率成分，如333 Hz、1 000 Hz等，引入的頻率成分幅值隨動偏心率的增加而增大。

表5為各偏心下氣隙磁密的諧波幅值，整體而言，偏心使得諧波幅值發(fā)生改變，其中靜偏心對主要諧波幅值的影響較大。以基波幅值（fc=667 Hz）為例，靜偏心時，基波幅值逐漸增大，動偏心時，基波幅值基本不發(fā)生改變?；旌掀臅r，當(dāng)靜偏心率保持不變，動偏心率由20%增加到60%，基波幅值變化較小；當(dāng)動偏心率保持不變，靜偏心率由20%增加到60%，基波幅值從0.92 T逐漸增大到1.07 T。

3.2 偏心對氣隙磁場空間諧波的影響

計算得到各偏心下電機氣隙磁場的空間分布，對于某一時刻電機磁場的空間分布，可將動偏心看作不同時間的靜偏心，所以本文僅繪制了靜偏心和保持靜偏心率不變，動偏心率從20%增加到60%的氣隙磁密空間分布曲線如圖9所示。從圖可知，隨著偏心率的增加，轉(zhuǎn)子偏向側(cè)的磁密強度逐漸增大，偏離側(cè)的磁密強度逐漸減小，并出現(xiàn)不同程度的畸變。

為進一步分析偏心對電機氣隙磁密空間分布的影響，對偏心下的氣隙磁密波形進行傅里葉分解。各偏心狀態(tài)下氣隙磁密諧波幅值如表6所示。由表可知，偏心后原有的奇數(shù)次諧波略微減小，偏心使得氣隙磁密增加了偶數(shù)次諧波，偶數(shù)次諧波幅值隨偏心率的增加而增大。

整體而言，偏心后電機的輸出轉(zhuǎn)矩平均值和轉(zhuǎn)矩脈動略微增大，符合式（16）的推導(dǎo)，但由于偏心數(shù)值較小，所以增幅較小。其中，對比三種偏心狀態(tài)對輸出轉(zhuǎn)矩平均值和轉(zhuǎn)矩脈動的影響發(fā)現(xiàn)，存在單一偏心故障時，靜偏心引起的轉(zhuǎn)矩脈動比動偏心引起的轉(zhuǎn)矩脈動要大。若為混合偏心故障，當(dāng)靜偏心的影響占主導(dǎo)因素即靜偏心率大于動偏心率，此時的轉(zhuǎn)矩脈動略大于動偏心占主導(dǎo)因素時的轉(zhuǎn)矩脈動，這是因為靜偏心時，氣隙最小位置固定不變，使得磁導(dǎo)和磁密分布具有傾向性，對氣隙磁場的影響較大。

4 各偏心故障下轉(zhuǎn)子振動特性分析

4.1 靜偏心對轉(zhuǎn)子振動的影響

電機偏心導(dǎo)致電磁力分布不均勻，從而產(chǎn)生不平衡磁拉力，靜偏心下不平衡磁拉力在X軸和Y軸的分量FX和FY如圖10所示。電機未發(fā)生偏心時，氣隙均勻分布，此時不平衡磁拉力為0。隨著偏心率從20%增至60%，F(xiàn)X的最大值從73 N增加到218 N，平均值由65 N增加到209 N，而FY的最大值和平均值變化不大，始終在0上下波動，這是因為偏心方向在X軸，不平衡磁拉力主要集中在X軸方向。

不平衡磁拉力作為轉(zhuǎn)子系統(tǒng)振動的激勵源，對轉(zhuǎn)子振動有著重要影響，對求得的不平衡磁拉力進行頻率特性分析如圖11所示。由于不平衡磁拉力是對徑向電磁力波沿圓周積分所得，所以其頻率特性與徑向電磁力波保持一致（見表2）。靜偏心下不平衡磁拉力頻率成分為2nfc的關(guān)系，各頻率對應(yīng)的幅值隨偏心率的增加而增大。

不平衡磁拉力將加劇轉(zhuǎn)子振動，分析轉(zhuǎn)子系統(tǒng)在不平衡磁拉力作用下的動力響應(yīng)，得到如圖12所示的軸心軌跡。由圖可知，靜偏心時轉(zhuǎn)子軸心軌跡的形態(tài)不規(guī)則，隨著偏心率的增加轉(zhuǎn)子振動位移增大，且軌跡中心偏向X軸正方向，這是由于偏心使得X軸正方向氣隙變小，從而不平衡磁拉力增大，故軌跡中心向X軸正方向偏移。

以最大靜偏心率60%為例，對X、Y方向的振動響應(yīng)和頻譜進行分析如圖13所示。位移頻譜中最大振動位移對應(yīng)零頻成分，其次是一階固有頻率和轉(zhuǎn)頻成分。因為發(fā)生靜偏心時，原有0階力波偏心產(chǎn)生1階力波（見表2），1階力波對應(yīng)頻率為0 Hz，又因低階次力波更易產(chǎn)生振動，故零頻下的振動幅值較大。

4.2 動偏心對轉(zhuǎn)子振動的影響

由解析式（11）可知，動態(tài)偏心對氣隙磁場的影響較為復(fù)雜，因此求得動態(tài)偏心下電機的不平衡磁拉力如圖14所示。

當(dāng)電機發(fā)生動態(tài)偏心時，不平衡磁拉力FX和FY都呈現(xiàn)周期性變化，二者周期相等，幅值相等，但在相位上相差90°，不平衡磁拉力最大值隨偏心率的增加而增大，動偏心率從0%增大到60%，最大值從0增加到190 N。

由上文可知，動偏心率僅影響不平衡磁拉力幅值，故以動偏心率60%為例對其頻率成分進行分析，如圖15所示。動態(tài)偏心下不平衡磁拉力的主要頻率成分為2nfc，在此基礎(chǔ)上產(chǎn)生了±fr和±2fr的頻率成分，符合解析式（11）推導(dǎo)的規(guī)律。相比靜態(tài)偏心，零頻對應(yīng)的諧波幅值明顯降低，而轉(zhuǎn)動頻率fr對應(yīng)的諧波幅值突出，高階頻率的諧波幅值與靜偏心相同且均較小。

各動態(tài)偏心率下的軸心軌跡如圖16所示，隨著偏心率的增大，轉(zhuǎn)子的振動幅值增加，但轉(zhuǎn)子軸心軌跡的中心不隨偏心量改變，保持同心狀態(tài)，原因是中心位置主要由外加作用力的平均值決定，而動偏心下不平衡磁拉力的平均值不發(fā)生變化。圖17為動偏心率60%時的X方向的振動響應(yīng)及頻譜圖，由圖可知，動偏心時轉(zhuǎn)子在一階固有頻率處的振幅最大，其次為轉(zhuǎn)頻。

4.3 混合偏心對轉(zhuǎn)子振動的影響

計算混合偏心狀態(tài)下的不平衡磁拉力如圖18所示，當(dāng)發(fā)生混合偏心后，不平衡磁拉力幅值隨著偏心程度的增加而增大，且均呈現(xiàn)周期性變化。

對混合偏心狀態(tài)下的轉(zhuǎn)子系統(tǒng)進行動力響應(yīng)分析，得到額定轉(zhuǎn)速下的軸心軌跡如圖19所示。圖19（a）、圖19（b）、圖19（c）中為保持靜偏心率20%不變，動偏心率從20%增加到60%的軌跡，對比圖12（a）靜偏心率20%的軌跡可知，動偏心的加入使得軌跡振幅增大，軌跡形態(tài)規(guī)則化。圖19（a）、圖19（d）、圖19（e）中為保持動偏心率20%不變，靜偏心率從20%增加到60%的軌跡，對比圖16（a）動偏心率20%的軌跡可知，隨著靜偏心率的增加，軌跡X方向的振幅逐漸增大，但Y方向的振幅變化不明顯，這是因為靜偏心時Y方向的不平衡磁拉力始終較小。

圖19（f）為動靜偏心率均為40%的軌跡，對比圖12（b）靜偏心率40%、圖16（b）動偏心率40%的軌跡，可知偏心程度的加劇，使轉(zhuǎn)子振動位移增大。

分析混合偏心下轉(zhuǎn)子的振動頻譜，以動靜偏心率相差較大的混合偏心為例，得到圖20所示兩種混合偏心下X方向的振動響應(yīng)及頻譜。從圖中可以看出，零頻對應(yīng)的諧波幅值隨靜偏心的增加而增大，轉(zhuǎn)頻fr的諧波幅值隨動偏心的增加而增大。混合偏心下的主要頻率成分是動、靜偏心的結(jié)合，沒有再引入其他的頻率成分，同時高階頻率幅值非常小，可忽略不計。

5 結(jié) 論

本文首先設(shè)計了一臺1.2 MW，20 000 r/min的高速永磁同步電機，然后針對電機氣隙偏心問題，分別建立了動、靜、混合偏心狀態(tài)下的電機模型，通過有限元法，計算了各偏心下的不平衡磁拉力，研究了不同偏心情況下的電磁性能和轉(zhuǎn)子動力響應(yīng)，得到以下結(jié)論：

1）電機存在偏心故障時，氣隙磁密偶次諧波明顯增加，其各諧波頻率幅值主要受靜偏心的影響，動偏心使得氣隙磁密引入kfc±fr的頻率成分。由于氣隙偏心后氣隙磁導(dǎo)發(fā)生變化，進而影響電機轉(zhuǎn)矩，使得電機轉(zhuǎn)矩脈動增加，轉(zhuǎn)矩平均值略微增大。

2）靜偏心對不平衡磁拉力頻率成分不產(chǎn)生影響，動偏心的不平衡磁拉力附加了2nfc±fr、±2fr的頻率成分，混合偏心時的不平衡磁拉力也附加了2nfc±fr、±2fr的轉(zhuǎn)頻成分，且不平衡磁拉力均呈現(xiàn)周期性變化。

3）動、靜偏心時轉(zhuǎn)子徑向位移振幅隨偏心量的增加而增大，但靜偏心量影響轉(zhuǎn)子軌跡的中心位置，而動偏心量對轉(zhuǎn)子軌跡的中心不產(chǎn)生影響?；旌掀臅r，動偏心對軸心軌跡形態(tài)的影響占主導(dǎo)作用。

參 考 文 獻：

［1］ 唐任遠. 現(xiàn)代永磁電機—理論與設(shè)計［M］. 北京： 機械工業(yè)出版社， 2005.

［2］ 王天煜. 高速永磁電機轉(zhuǎn)子綜合設(shè)計方法及動力學(xué)特性的研究［D］. 沈陽： 沈陽工業(yè)大學(xué)， 2010.

［3］ 孔漢. 永磁同步電機故障對電機綜合物理場的影響機理研究［D］. 西安： 西北工業(yè)大學(xué)， 2016.

［4］ 于占洋， 李巖， 井永騰. 基于混合磁場解析法的磁極偏心型表貼式永磁同步電機空載特性分析［J］. 電工技術(shù)學(xué)報， 2020， 35（18）： 3811.

YU Zhanyang， LI Yan， JING Yongteng. No-load characteristics analysis of magnetic pole eccentric surface mount permanent magnet synchronous motor based on hybrid magnetic field analysis method［J］. Transactions of China Electrotechnical Society， 2020， 35（18）： 3811.

［5］ BEHREND B A. On the mechanical forces in dynamos caused by magnetic attraction ［J］. Transactions of the American Institute of Electrical Engineers， 1900， （XVII）： 613.

［6］ KIM U， LIEU D K. Magnetic field calculation in permanent magnet motors with rotor eccentricity： without slotting effect［J］. IEEE Transactions on Magnetics， 1998， 34（4）： 2243.

［7］ KIM T J， HANG S M， PARK N G. Analysis of vibration for permanent magnet motors considering mechanical and magnetic couple effects［J］.IEEE Transactions on Magnetics，2000，36（4）：1346.

［8］ 郭丹， 何永勇， 褚福磊. 不平衡磁拉力及對偏心轉(zhuǎn)子系統(tǒng)振動的影響［J］. 工程力學(xué)， 2003， 20（2）： 116.

GUO Dan， HE Yongyong， CHU Fulei. The calculation of unbalanced magnetic pull and its effect on vibration of an eccentric rotor［J］. Engineering Mechanics， 2003， 20（2）： 116.

［9］ 岳二團， 甘春標(biāo)， 楊世錫. 氣隙偏心下永磁電機轉(zhuǎn)子系統(tǒng)的振動分析［J］. 振動與沖擊， 2014， 33（8）： 29.

YUE Ertuan， GAN Chunbiao， YANG Shixi. Vibration characteristics analysis of a rotor for a permanent magnet motor with air-gap eccentricity［J］.Journal of Vibration and Shock，2014，33（8）：29.

［10］ 徐學(xué)平， 韓勤鍇， 褚福磊， 靜載荷作用下偏心轉(zhuǎn)子電磁振動特性［J］. 清華大學(xué)學(xué)報（自然科學(xué)版）， 2016， 56（2）： 176.

XU Xueping， HAN Qinkai， CHU Fulei. Electromagnetic vibration characteristics of an eccentric rotor with a static load［J］. Journal of Tsinghua University（Science and Technology）， 2016， 56（2）： 176.

［11］ 馮立帥. 轉(zhuǎn)子偏心對低速永磁同步電機運行特性影響的研究［D］. 沈陽： 沈陽工業(yè)大學(xué)， 2022.

［12］ 蔣山. 偏心狀態(tài)下永磁同步電機電磁性能分析［D］. 哈爾濱： 哈爾濱理工大學(xué)， 2017.

［13］ 周生通， 郭維年， 肖乾. 偏心故障下動車組牽引電機轉(zhuǎn)子振動特性分析［J］. 工程力學(xué)， 2021， 38（7）： 216.

ZHOU Shengtong， GUO Weinian， XIAO Qian. Vibration characteristics analysis of EMU traction motor rotor under eccentric faults［J］. Engineering Mechanics， 2021， 38（7）： 216.

［14］ 黨文娟， 不平衡磁拉力作用下偏心轉(zhuǎn)子-軸承系統(tǒng)非線性振動特性分析［D］. 西安： 西安理工大學(xué)， 2019.

［15］ 張霄霆， 張炳義， 曹永鵬. 偏心狀態(tài)下的軸向磁通永磁電機受力分析［J］. 電工技術(shù)學(xué)報， 2020， 35（1）： 110.

ZHANG Xiaoting， ZHANG Bingyi， CAO Yongpeng， Force analyze of axial magnetic flux permanent magnet motor under eccentricity condition［J］. Transactions of China Electrotechnical Society， 2020， 35（1）： 110.

［16］ 唐貴基， 鄧瑋琪， 何玉靈. 不同種類氣隙偏心故障對汽車發(fā)電機轉(zhuǎn)子不平衡磁拉力的影響［J］. 振動與沖擊， 2017， 36（15）： 1.

TANG Guiji， DENG Weiqi， HE Yuling. Effect of different kinds of air-gap eccentricity faults on rotor UMP of a turbo-generator［J］. Journal of Vibration and Shock， 2017， 36（15）： 1.

［17］ 孫晶晶. 考慮轉(zhuǎn)子偏心的永磁同步電機電磁力特性分析［D］. 杭州： 浙江大學(xué)， 2016.

［18］ 鐘一諤， 何衍宗. 轉(zhuǎn)子動力學(xué)［M］. 北京： 清華大學(xué)出版社， 1987.

［19］ FENG Lishuai， YU Siyang， ZHANG Fengge. Study on performance of low-speed high-torque permanent magnet synchronous motor with dynamic eccentricity rotor ［J］. SCI Energy Reports， 2022， 8（5）： 1421.

［20］ 田擁勝， 孫巖樺， 虞烈. 高速永磁電機電磁軸承轉(zhuǎn)子系統(tǒng)的動力學(xué)及實驗研究［J］. 中國電機工程學(xué)報， 2012， 32（9）： 116.

TIAN Yongsheng， SUN Yanhua， YU Lie. Dynamical and experimental researches of active magnetic bearing rotor systems for high speed PM machines［J］. Proceedings of the CSEE， 2012， 32（9）： 116.

［21］ 喬芳， 盤軸-松動偏心耦合故障轉(zhuǎn)子系統(tǒng)動力學(xué)建模與振動特性研究［D］. 南昌： 南昌航空大學(xué)， 2021.

［22］ SHEN Quan， ZHOU Jin， MA Yanchao， et al. Dynamic modeling and analysis of magnetic suspension rotor system under foundation excitation［J］.Journal of Vibration and Shock，2022，41（17）：35.

［23］ MA Conggan， CUI Hongwei， ZHENG Ping. Influence of static eccentricity on unbalanced magnetic force of external rotor permanent magnet brushless direct current motor used as in-wheel motor［J］. IET Electric Power Applications， 2019， 13（4）： 538.

［24］ 周曉燕， 李琛， 仇志堅， 等. 基于攝動法的交替級永磁電機偏心磁場解析計算［J］. 電工技術(shù)學(xué)報， 2013， 28（9）： 321.

ZHOU Xiaoyan， LI Chen， QIU Zhijian， et al. Magnetic field calculation for consequent-pole motor with rotor eccentricity based on perturbation method［J］. Transactions of China Electrotechnical Society， 2013， 28（9）： 321.

［25］ 張巖巖， 耿海鵬， 周健， 等. 永磁同步電機轉(zhuǎn)子偏心磁場解析計算［J］. 西安交通大學(xué)學(xué)報， 2019， 53（2）： 93.

ZHANG Yanyan， GENG Haipeng， ZHOU Jian， et al. Analytical calculation for magnetic field of permanent magnet synchronous motors with rotor eccentricity［J］. Journal of Xi’an Jiaotong University， 2019， 53（2）：93.

［26］ LI Yanxin， LU Qinfen， ZHU Z Q. Unbalanced magnetic force prediction in permanent magnet machines with rotor eccentricity by improved superposition method ［J］. IET Electric Power Applications， 2017， 11（6）： 1095.

［27］ 戈寶軍， 劉海濤， 王立坤. LNG泵用低溫高速永磁電機偏心故障動態(tài)電磁力［J］. 電機與控制學(xué)報， 2022， 26（5）： 55.

GE Baojun， LIU Haitao， WANG Likun. Dynamic electromagnetic force of cryogenic high-speed permanent magnet motors with eccentric fault［J］.Electric Machines and Control，2022，26（5）：55.

［28］ 劉蓉暉， 劉錦坤， 章君達. 基于雙曲余切變換的Halbach陣列表貼式永磁電機轉(zhuǎn)子偏心氣隙磁場解析模型［J］. 電工技術(shù)學(xué)報， 2023， 38（6）： 1433.

LIU Ronghui， LIU Jinkun， ZHANG Junda. Analytical model for air-gap magnetic field in Halbach arrays surface-mounted permanent magnet motor with rotor eccentricity based on hyperbolic cotangent transformation［J］. Transactions of China Electrotechnical Society， 2023， 38（6）： 1433.

［29］ 韓雪巖， 王勇， 高俊. 低速永磁同步電機轉(zhuǎn)子偏心的抑制措施［J］. 電機與控制學(xué)報， 2023， 27（11）： 58.

HAN Xueyan， WANG Yong， GAO Jun. Suppression of rotor eccentricity of low-speed permanent magnet synchronous motor［J］. Electric Machines and Control， 2023， 27（11）： 58.

（編輯：劉素菊）