高速永磁電機轉子系統(tǒng)臨界轉速仿真研究

關 棟，楊小輝，劉 更，佟瑞庭，馬尚君

(西北工業(yè)大學，陜西西安710072)

0 引 言

隨著航空、航天等領域對永磁電機功率重量比和體積等要求的不斷提高，日益通過提高轉速的方式達到縮小體積、省去傳動裝置的目的［1-2］。電機轉速提高不可避免加劇轉子系統(tǒng)的振動，甚至可能接近或達到其臨界轉速而降低可靠性。因此，對高速永磁電機轉子系統(tǒng)的臨界轉速分析就成為設計時的重要內(nèi)容［3-8］。準確計算出高速電機轉子系統(tǒng)的臨界轉速，給出提高其臨界轉速的工程措施就顯得尤為重要。

鐘一諤［9］、聞邦椿［10］、虞烈［11］等學者在軸承 -轉子系統(tǒng)、臨界轉速、響應計算等方面做了大量的研究工作，對旋轉機械的工程設計提供了豐富的理論基礎。但是這些研究主要集中于壓縮機、水輪機、汽輪發(fā)電機等轉子系統(tǒng)，而對中小功率的電機轉子研究較少。文獻［3］～［5］、［8］、［12］～［17］采用有限元法分別研究了開關磁阻電機和永磁電機轉子的動力學特性。劉闖［3，13］和王鳳翔［8，15-16］分別研究了陶瓷滾動軸承和磁力軸承支撐結構對電機轉子臨界轉速和軸心軌跡等的影響。其中王鳳翔［15-16］等人進一步分析了轉子結構對電機臨界轉速的影響，并采用樣機振動實驗驗證其計算方法的正確性。文獻［17］建立3D有限元模型計算永磁無刷電機轉子系統(tǒng)的固有頻率，并考慮陶瓷軸承剛度對高速電機固有頻率的影響。

由于工作原理的差別，文獻［3］、［13］所研究的開關磁阻電機轉子工作原理和結構形式上與本文探討的永磁電機轉子結構形式有所不同，故其動力學特性也有所差別。而文獻［15］～［17］所研究磁力軸承和陶瓷滾動軸承的支撐結構形式在航空航天領域較少采用。本文在分析滾動軸承支撐特性的基礎上，采用有限元法對電機轉子進行動力學建模，模型中考慮轉子高速旋轉時的陀螺效應。針對某航天用永磁電機轉子進行動力學仿真分析，并采用數(shù)值計算的方法對仿真結果進行驗證和分析。最后研究了采用滾動軸承支撐時不同支撐剛度對永磁電機臨界轉速的影響，給出提高電機臨界轉速的工程措施。

1 滾動軸承－轉子動力學模型

1.1 模型示意圖

永磁電機磁鋼所采用的燒結釹鐵硼材料具有較好的抗壓能力。如圖1所示，一般在磁鋼表面安裝非導磁材料的護環(huán)，護環(huán)與磁鋼之間采用過盈配合，這樣可以保證永磁轉子在高速下的安全運行。

圖1 轉子系統(tǒng)動力學模型

由圖1可知，轉軸兩端由2個滾動軸承支撐并考慮軸承的彈性支撐及阻尼效應。同時考慮轉子上的質量和磁拉力產(chǎn)生的不平衡效應。電機轉子高速旋轉時，由于上述兩項不平衡力的激勵作用將促使轉子產(chǎn)生橫向振動。圖中，m為轉子不平衡質量;e為偏心距;C1、K1、C2、K2分別表示轉子與前后兩個軸承連接處的阻尼和支撐剛度;Fx1、Fy1為前端軸承的支撐反力;Fx2、Fy2為后端軸承的支撐反力;o1、o2為轉子的幾何中心和轉子質心。

1.2 動力學方程

航空航天用高速永磁電機的主要特點是高速旋轉的轉子圓周速度可達200 m/s以上，是定子中繞組電流和鐵心磁通具有很高的頻率。由于電機轉子在高速旋轉時具有強烈的自穩(wěn)定能力即陀螺效應，而且高速旋轉時的陀螺效應對其固有頻率、臨界轉速等特性影響很大。為了更好地分析此轉子動力學特性，在考慮陀螺效應的情況下其動力學方程如下［15，18］:

式中:M、C、K、X、Y分別是轉子的質量、阻尼、剛度和節(jié)點在x、y兩個方向的位移矩陣(向量);Fc和Fb分別表示旋轉時不平衡質量產(chǎn)生的離心力和電機工作時轉子所受的不平衡磁拉力;G是反對稱陣，為考慮慣性力的陀螺矩陣。轉子所受離心力，陀螺矩陣和轉速矩陣如下:

結合圖1，上述方程中m為轉子不平衡質量，e表示偏心距，ρ表示密度，N為單元形函數(shù)矩陣。式(1)的系數(shù)矩陣是由單元矩陣組集合而成，由于永磁直流電機在高速下旋轉，公式考慮了電機工作時轉子慣性力的作用，能夠更加精確地模擬其真實的工作情況。

2 永磁電機滾動軸承－轉子臨界轉速計算

2.1 模型建立及邊界條件的設定

2.1.1 轉子相關參數(shù)

本文以某電機轉子系統(tǒng)為例，其結構形式如圖1所示。兩端采用普通深溝球軸承支撐，轉子相關參數(shù)如表1所示。

表1 電機參數(shù)

2.1 .2模型轉換及邊界條件設置

電機轉子在滾動軸承連接處的模型如圖2所示。

軸頸與軸承配合處采用旋轉副，轉子軸頸在旋轉副的約束下僅有一個旋轉自由度，如圖2中速度 v所示。與軸頸配合的旋轉副四周采用彈簧-阻尼單元模擬滾動軸承支撐。旋轉副外側與彈簧連接點添加軸向約束，彈簧另一端全約束。因此設定軸承剛度時，可以僅考慮軸承的徑向剛度，忽略其軸向剛度。由表1完成相關仿真參數(shù)的設定。

圖2 轉子支撐簡化模型

2.2 臨界轉速計算

轉子的固有頻率除了與轉子結構和支撐參數(shù)有關外，還隨轉子渦動轉速和自轉轉速的變化而變化［19］。因為不平衡力的存在，轉子通常作正同步渦動。當其渦動轉速和固有頻率相等時即出現(xiàn)共振，與之對應的轉速即為臨界轉速。

本文采用有限元法對轉子的臨界轉速進行計算時出現(xiàn)正進動和反進動兩種狀態(tài)。考慮到陀螺效應的影響，正進動的固有頻率隨著轉速的升高而增大，反進動的固有頻率將隨著轉速的升高而下降。在對永磁電機轉子臨界轉速進行分析時應只分析其正進動的固有頻率。

圖3、圖4分別為電機轉子的坎貝爾圖和其1階固有頻率隨轉速的變化圖。圖中橫坐標表示轉子的轉速，縱坐標表示轉子的固有頻率。圖3給出了轉子前四階固有頻率隨轉速的變化趨勢。其中第2階為反進動的狀態(tài)，可以不予考慮。第1、3、4階為轉子的正進動，結合圖4可以看出，正進動時轉子的固有頻率隨著其轉速的升高而增大。其中第1階固有頻率的增幅較小，而3、4階增速變化明顯。與轉子1階固有頻率相交于521.48Hz處，由nc=60f可知，其臨界轉速約為31 289 r/min。為保證永磁電機良好的運轉性能，避免轉子發(fā)生共振，要求電機的實際轉速滿足如下條件:

圖3 坎貝爾圖

圖4 轉子1階固有頻率隨轉速的變化圖

此處，0.75Nc1為 23 466.6 r/min。由表1 知，電機的額定轉速為16 000 r/min，計算得到的1階臨界轉速遠大于其額定轉速，故電機轉子為剛性轉子，而且在額定轉速工作時不會因轉子本身運轉而發(fā)生共振現(xiàn)象。

2.3 結果驗證

為驗證仿真結果的可信性，本文對電機轉子的1階臨界轉速進行理論計算，相關公式如下［20］:

式中:E為轉子軸的楊氏模量;Gi表示軸上不同截面的自重及其所受載荷;l為支撐跨度;式(6)中J的表達式:

式中:di表示軸上的不同截面直徑。將轉子結構尺寸參數(shù)代入式(6)、式(7)可計算其1階臨界轉速。仿真值與理論值對比如表2所示。

表2 仿真和計算值對比

由表2可以看出，在不考慮陀螺效應的情況下仿真值與理論值比較接近，誤差僅為1.32%。但是考慮陀螺效應后，采用仿真方法計算出的臨界轉速明顯大于理論計算結果，誤差達到10.39%。出現(xiàn)這樣的結果可以用轉子高速旋轉時產(chǎn)生的陀螺效應解釋。由于式(6)在計算臨界轉速時未考慮陀螺效應的作用，所以其計算的臨界轉速無論是否考慮陀螺效應均為同一值。如圖4所示，隨著電機轉子轉速升高，其1階固有頻率逐步上升，在采用考慮陀螺效應的有限元仿真方法計算時，其值明顯大于傳統(tǒng)理論計算結果。因此在計算高速永磁電機轉子系統(tǒng)臨界轉速時必須將陀螺效應考慮在內(nèi)。

3 軸承剛度對臨界轉速的影響

通常情況下高速永磁電機可采用滾動軸承、磁力軸承等方式支撐其軸系。但是綜合考慮可靠性等方面的因素，目前航空航天用永磁電機支撐結構大多仍采用滾動軸承。

軸承剛度是表征軸承動態(tài)特性的重要參數(shù)，其大小及變化特性將直接影響軸承-轉子系統(tǒng)的臨界轉速。因此，研究軸承剛度的支撐特性對永磁電機轉子系統(tǒng)來說是非常必要的。本文主要考慮滾動軸承支撐對電機轉子臨界轉速的影響，結合前面所建模型主要考慮滾動軸承徑向剛度改變時電機前兩階臨界轉速的變化趨勢。

采用前文所述永磁電機轉子作為算例，取初始軸承剛度為7.5×107N/m，忽略阻尼對臨界轉速的影響，計算出其1、2階臨界轉速分別為31 289 r/min和55 462 r/min。在一定范圍內(nèi)逐漸提高滾動軸承的支撐剛度，所得到的電機轉子1、2階臨界轉速如圖5所示。從圖中可知，隨著支撐剛度的提高臨界轉速明顯升高，而且轉子的2階臨界轉速增速明顯大于其1階臨界轉速。

圖5 軸承剛度對臨界轉速的影響

為進一步探究這種變化趨勢，繼續(xù)加大軸承支撐剛度，如表3所示。由表3知，當滾動軸承的支撐剛度在7.5×107～5×109N/m范圍內(nèi)變化時，其前兩階固有頻率增速較為明顯。而當軸承剛度大于5×109N/m以后，固有頻率增速放緩。

表3 軸承剛度對固有頻率的影響

由上述理論分析可知，提高永磁電機臨界轉速可通過增加軸承支撐剛度實現(xiàn)。目前航空航天領域采用普通滾動軸承支撐的電機臨界轉速具有較大的提升空間。對于同型號的滾動軸承，可通過適當增加軸承預緊力的方式提高轉子系統(tǒng)的臨界轉速。同時還可在設計階段采用較大型號或寬系列的軸承提高轉子系統(tǒng)的支撐剛度，以達到提高系統(tǒng)臨界轉速的目的。但是加大軸承尺寸或寬度同時又將降低其極限轉速，與提高永磁電機的轉速產(chǎn)生矛盾，因此需要電機設計人員在設計階段根據(jù)需要進行取舍。當普通軸承剛度不能滿足轉子系統(tǒng)對臨界轉速的要求時可采用陶瓷滾動軸承取代，這樣可以大幅度提高轉子系統(tǒng)的臨界轉速，使電機轉子在更高轉速下能夠穩(wěn)定運轉。

4 結 語

通過對考慮陀螺效應的高速永磁電機滾動軸承-轉子系統(tǒng)進行動力學建模和有限元仿真分析，可得到如下結論:

(1)由于航空航天用高速永磁電機滾動軸承-轉子系統(tǒng)結構固有的特性，在對其進行臨界轉速計算時必須考慮陀螺效應的影響，尤其是高階臨界轉速更不能忽略。

(2)提高滾動軸承剛度可大幅度增加永磁電機轉子系統(tǒng)的臨界轉速。所以目前航空航天領域在用永磁電機可以通過增加軸承剛度的方式達到提高臨界轉速的目的。

(3)對于航空航天用永磁電機的設計可以通過增加軸承預緊力、增大軸承尺寸、采用寬系列軸承或高可靠性的陶瓷滾動軸承取代普通滾動軸承等方式增加其臨界轉速。但是增大軸承尺寸或寬度又影響到軸承極限轉速的提高，不利于提高永磁電機轉子轉速。因此需要設計人員根據(jù)高速永磁電機的應用場合在增加軸承剛度和提高臨界轉速間折中取舍。

(4)隨著軸承支撐剛度的增加，永磁電機轉子前兩階固有頻率之間的工作區(qū)域逐步變大。這為航空航天用永磁電機繼續(xù)提高工作轉速而在前兩階臨界轉速間工作提供了可能性。

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

［2］ 唐任遠，趙清，周挺.稀土永磁電機正進入大發(fā)展的新時期［J］.沈陽工業(yè)大學學報，2011，33(1):1-8.

［3］ 劉闖，朱旭勇，卿湘文.開關磁阻電機轉子動力學建模與分析［J］.中國電機工程學報，2008，28(3):83-89.

［4］ Bae D，Chung J.Design sensitivity analysis of dynamic responses for a BLDC motor with mechanical and electromagnetic interactions［J］.nal of Sound and Vibration，2012，331(9):2070-2079.

［5］ Srinivas K N，Arumugam R.Static and dynamic vibration analyses of switched reluctance motors including bearings，housing，rotor dynamics，and applied Loads［J］.IEEE Transactions on Magnetics，2004，40(4):1911-1919.

［6］ Werner U.Theoretical rotor dynamic analysis of two-pole induction motors regarding excitation due to static rotor eccentricity［J］.Archive of Applied Mechanics，2011，81(2):241-262.

［7］ 王鳳翔，高速電機的設計特點及相關技術研究［J］.沈陽工業(yè)大學學報，2006，28(3)，258-263.

［8］ 王繼強，王鳳翔，宗鳴.高速電機磁力軸承-轉子系統(tǒng)臨界轉速的計算［J］.中國電機工程學報，2007，27(27):94-98.

［9］ 鐘一鄂，何衍宗，王正，等.轉子動力學［M］.北京:清華大學出版社，1987.

［10］ 聞邦椿，顧家柳，夏松波，等.高等轉子動力學［M］.北京:機械工業(yè)出版社，2000.

［11］ 虞烈，劉恒.軸承-轉子動力學［M］.西安:西安交通大學出版社，2001.

［12］ Kyung H H，Jung P H，Gyu H T，et al.Orbital analysis of rotor due to electromagnetic force for switched reluctance motor［J］.IEEE Transactions on Magnetic，2000，36(4):1400-1407.

［13］ 周強，劉闖，朱學忠，等.超高速開關磁阻電動機設計［J］.中國電機工程學報，2009，29(9):87-92.

［14］ Kim T J，Hwang S M，Park N G.Analysis of vibration for permanent magnet brushless machines motors considering mechanical and magnetic coupling effects［J］.IEEE Transactions On Magnetic，2000，36(4):1346-1350.

［15］ 王天煜，王鳳翔，方程.高速永磁電機機組軸系振動研究［J］.振動與沖擊，2011，30(9):111-115.

［16］ 王天煜，王鳳翔，方程，等.高速電機轉子臨界轉速計算與振動模態(tài)分析［J］.遼寧工程技術大學學報(自然科學版)，2009，28(5):804-808.

［17］ Ede J D，Zhu Z Q，Howe D.Rotor resonances high-speed permanent magnet brushless machines［J］.IEEE Transactions on Industry Applications，2002，38(6):1542-1548.

［18］ Greenhill L M，Bickford W B，Nelson H D.A Conical beam finite elememt for rotor dynamics analysis［J］.Journal of Vibration，Stress，and Reliability in Design，1985，107(10):421-430.

［19］ 顧家柳，夏松波，張文.轉子動力學研究的現(xiàn)狀及展望［J］.振動工程學報，1988，1(2):63-70.

［20］ 關醒凡.現(xiàn)代泵技術手冊［M］.北京:宇航出版社，1995.