雷小葉，盧南方，齊景榮

（1．貴州工業(yè)職業(yè)技術(shù)學(xué)院 機械與電氣工程學(xué)院，貴州 貴陽 550025；2.貴州電子信息職業(yè)技術(shù)學(xué)院 機電工程系，貴州 凱里 556000）

高速永磁無刷電機作為當(dāng)前電機研究領(lǐng)域的熱點之一，其主要特點表現(xiàn)為體積小、功率密度大和效率高等方面。與普通電機轉(zhuǎn)子設(shè)計相比，高速電機為了降低離心力和提高輸出功率，轉(zhuǎn)子結(jié)構(gòu)一般被設(shè)計為細長型。結(jié)構(gòu)一長，剛度降低，為了降低并避開轉(zhuǎn)子臨界轉(zhuǎn)速的振動頻率，轉(zhuǎn)子結(jié)構(gòu)又需要有足夠的剛度。對于高速永磁電機來說，轉(zhuǎn)子上有永磁體，整個轉(zhuǎn)子系統(tǒng)受到不平衡磁拉力的影響，因此，對電機轉(zhuǎn)子磁拉力的精確計算和動力學(xué)分析是高速永磁無刷電機設(shè)計的關(guān)鍵技術(shù)之一。

以某型90000r/min高速永磁無刷電機轉(zhuǎn)子軸系設(shè)計為例，采用NX軟件建立三維模型，通過ANSYS對轉(zhuǎn)子結(jié)構(gòu)進行磁固耦合分析轉(zhuǎn)子的模態(tài)頻率和振型，分析了轉(zhuǎn)子振動特性，為高速永磁高速電機轉(zhuǎn)子設(shè)計提供參考。

1 轉(zhuǎn)子有限元模型的建立

在有限元分析前，采用NX對高速電機外轉(zhuǎn)子進行三維幾何模型建模設(shè)計，轉(zhuǎn)子結(jié)構(gòu)由主軸、釹鐵硼磁鋼、鋁制隔板、兩塊擋板、套筒以及跟蹤磁環(huán)等組成。電機轉(zhuǎn)子磁鋼采用表貼式，為防止磁鋼在高速下飛散，采用熱套方式加護套保護，其結(jié)構(gòu)如圖所示。

圖1 電機轉(zhuǎn)子結(jié)構(gòu)N X模型

材料參數(shù)的選擇是有限元計算結(jié)果誤差的主要來源。通過查閱機械設(shè)計手冊以及材料供應(yīng)商提供的資料，所選各部件的材料及相關(guān)參數(shù)如表1所示。

表1 各部件材料屬性

將NX軟件中建好的三維模型，存儲為“.x_t”格式，將其導(dǎo)入ANSYS Workbench中進行網(wǎng)格劃分。有限元分析中，網(wǎng)格的質(zhì)量和離散點數(shù)對計算精度影響很大，在盡可能用六面體網(wǎng)格的同時對出現(xiàn)模態(tài)振型的轉(zhuǎn)子軸伸處和軸尾處網(wǎng)格加密，有限元模型如圖2所示。

圖2 轉(zhuǎn)子有限元模型

圖3 電磁仿真

2 電磁力數(shù)值計算

對轉(zhuǎn)子進行自振分析得到的是轉(zhuǎn)子結(jié)構(gòu)無軸承支撐狀態(tài)下的模態(tài)，文獻2-4對永磁電機轉(zhuǎn)子自由模態(tài)均作出了分析。但是在實際工作中電機轉(zhuǎn)子受著軸承預(yù)緊力、電磁力場、引力場、旋轉(zhuǎn)離心力結(jié)構(gòu)場等耦合所施加的預(yù)應(yīng)力，這些預(yù)應(yīng)力會對模態(tài)中各階固有頻率有不同程度的降低影響。

本文通過ANSYS maxwell自帶的虛功法計算電磁力密度。虛功法通過公式（1）將有限元中各網(wǎng)格單元所受到的力等效到節(jié)點上，提取節(jié)點上的集中力，計算力密度。

圖4 磁固耦合分析

圖5 轉(zhuǎn)子固有頻率與振型

圖6 結(jié)構(gòu)優(yōu)化示意圖

式中：W為磁場能，s為節(jié)點虛位移，i為節(jié)點數(shù)。

提取了電機模型在電磁環(huán)境下產(chǎn)生的電磁力場，如圖3所示。

3 模態(tài)計算分析

在ANSYS Workbench平臺下與Static Structural模塊的結(jié)構(gòu)場耦合進行預(yù)應(yīng)力計算，分析在約束下的轉(zhuǎn)子模態(tài)如圖4所示。

通過對電機轉(zhuǎn)子施加約束后進行仿真計算，得到轉(zhuǎn)子實際情況下的模態(tài)；分析轉(zhuǎn)子各階固有頻率下轉(zhuǎn)子結(jié)構(gòu)的變化趨勢，從而為進一步的動力學(xué)分析，如臨界轉(zhuǎn)速和諧響應(yīng)分析提供理論支持。若要改善機構(gòu)的振動狀況，可從振型分析中觀察，對結(jié)構(gòu)振動影響較大的部位，從而通過設(shè)計上的改進或其它約束方法來改善變形較大的部位。

軸承一般為滾動軸承，它的剛度系數(shù)為

,阻尼很小，通常忽略不計。本文在軸肩位置施加的軸承約束取最小剛度系數(shù)模擬電機轉(zhuǎn)子實際支撐情況，設(shè)置軸承的剛度矩陣K11與K22值均為2×107N/m；取軸承跨度為安裝空間最大允許跨度，可得最低情況的振動特性。在模態(tài)的計算中，由于低階模態(tài)振型幅度較大，取前六階，其計算結(jié)果如表2所示。

表2 轉(zhuǎn)子前六階固有頻率

在高速電機轉(zhuǎn)子的設(shè)計中，應(yīng)將轉(zhuǎn)子設(shè)計為剛性轉(zhuǎn)子，即電機轉(zhuǎn)子的實際工作頻率應(yīng)不高于轉(zhuǎn)子最低固有頻率，甚至更低來避免共振的趨勢。電機轉(zhuǎn)速為 nN=90000 rpm ，對應(yīng)轉(zhuǎn)子旋轉(zhuǎn)頻率為：

由表2可知，轉(zhuǎn)子旋轉(zhuǎn)頻率在上升過程中將會與一階、二階固有頻率相交并產(chǎn)生劇烈振動，與轉(zhuǎn)子動力學(xué)設(shè)計不符。其振型如圖5所示，振型成對出現(xiàn)，但不在一個平面內(nèi)，振動位置會發(fā)生在軸伸端。

圖7 靈敏度分析

圖8 響應(yīng)曲面圖

4 轉(zhuǎn)子結(jié)構(gòu)優(yōu)化

對高速永磁無刷電機轉(zhuǎn)子進行的模態(tài)分析是其它動力學(xué)分析的基礎(chǔ)，可知升高模態(tài)固有頻率可有效避免共振。電機轉(zhuǎn)子中有許多部位的結(jié)構(gòu)與其余零部件之間形成裝配關(guān)系，部分尺寸為功能尺寸，如為達到電機所需磁場強度的轉(zhuǎn)子磁鋼部位長度、軸尾跟蹤磁環(huán)位置，以及早已選定的軸承軸肩尺寸。因此在不影響電機轉(zhuǎn)子裝配尺寸及功能尺寸的前提下，著重對圖5振型位置作出優(yōu)化，如圖6所示。

當(dāng)修改轉(zhuǎn)子結(jié)構(gòu)時，即使某些尺寸因為給定的幾何關(guān)系或功能限制，仍然有4個尺寸參數(shù)對結(jié)構(gòu)有影響；對4個參數(shù)都進行逐一修改，計算成本仍較大。利用靈敏度分析找出與優(yōu)化振型變化相關(guān)程度較高的影響因子，這樣既實現(xiàn)了優(yōu)化的目的，又不會對設(shè)計尺寸修改較大，符合工程設(shè)計原則，靈敏度分析如圖7所示。

由圖7可以看出，優(yōu)化1段長度對振型影響最大，優(yōu)化1段直徑對固有頻率影響較大。通過ANSYS驅(qū)動優(yōu)化系統(tǒng)，對兩個優(yōu)化段尺寸在裝配空間允許的條件下設(shè)置尺寸變化參數(shù)作有限元分析，得到三者之間關(guān)系的曲面響應(yīng)，如圖8所示。

圖8表明，優(yōu)化1段直徑的尺寸與固有頻率關(guān)系存在一個最大值；優(yōu)化1段長度與振動頻率呈負線性關(guān)系，即取長度最小值，得固有頻率最大值。通過優(yōu)化后得到的尺寸，作再計算得系統(tǒng)一階固有頻率為1979.3Hz，相比于優(yōu)化前表2中結(jié)構(gòu)的一階固有頻率340.6Hz，提高了638.7Hz，有效避開了轉(zhuǎn)速頻率1500Hz，滿足轉(zhuǎn)子剛性的設(shè)計要求。

5 結(jié)語

高速永磁無刷電機的廣泛應(yīng)用促進了其快速發(fā)展，也對轉(zhuǎn)子振動、噪音問題提出了更高的要求。本文利用ANSYS Workbench對高速永磁無刷電機轉(zhuǎn)子結(jié)構(gòu)進行了模態(tài)分析并作出優(yōu)化，區(qū)別于傳統(tǒng)分析的是在模態(tài)分析前進行了電磁力分析，并通過有限元Workbench平臺進行了磁固耦合。這樣的分析，所施加的約束更接近工程實際。通過計算出的低階固有頻率與轉(zhuǎn)子轉(zhuǎn)速頻率作進一步優(yōu)化，提高了轉(zhuǎn)子結(jié)構(gòu)固有頻率，有效降低了轉(zhuǎn)子振動。整個分析與優(yōu)化過程可為電機產(chǎn)品的設(shè)計與改進提供一個參考。

[1]羅林,宋春華.微特電機發(fā)展綜述[J].橡塑技術(shù)與裝備,2016,(06):38-39.

[2]許貞俊,史忠震.基于A NS Y S的無刷直流電機轉(zhuǎn)軸的模態(tài)分析[J].現(xiàn)代機械,2013, (04):29-31.

[3]劉剛,張衛(wèi),陳強.基于A NS Y S的航空電機轉(zhuǎn)子振動分析[J].微特電機,2014,42(04):47-48+52.

[4]楊桃月,張賢信.基于A NS Y S的高速永磁電機轉(zhuǎn)子模態(tài)分析[J].機電技術(shù),2015, (05):36-38.