于慎波，李 紅，翟鳳軍，鐘雙雙，夏鵬澎

(沈陽工業(yè)大學 機械工程學院，沈陽 110870)

基于nCode Design-Life的永磁同步電機主軸的疲勞特性分析*

于慎波，李 紅，翟鳳軍，鐘雙雙，夏鵬澎

(沈陽工業(yè)大學 機械工程學院，沈陽 110870)

針對某型號永磁同步電機主軸局部出現(xiàn)裂紋或斷裂現(xiàn)象展開了分析研究。首先用解析法和數(shù)值法對電機主軸進行了二維電磁力計算分析，確定了作用在主軸上的載荷譜。然后利用nCode Design-Life中的S-N分析方法，對永磁同步電機的主軸進行疲勞特性分析。得到主軸的疲勞損傷云圖，結(jié)果顯示在主軸輸出端第一變截面處出現(xiàn)了應力集中。在瞬態(tài)電磁力和輸出脈動轉(zhuǎn)矩循環(huán)作用下，存在斷軸隱患。提出了在變截面處采用凹圓角的加工方法，可有效的減少應力集中，提高主軸壽命。為永磁同步電機的可靠性設計提供了參考。

疲勞分析；永磁同步電機；電磁力

0 引言

永磁同步電機作為驅(qū)動元件具有維護方便、可控性強、受環(huán)境影響小、電機效率及功率因數(shù)高等諸多優(yōu)點，被廣泛的應用于工業(yè)及生活中[1-2]。然而，由于電機轉(zhuǎn)子要安裝永磁體，其轉(zhuǎn)軸結(jié)構(gòu)尺寸受到限制。加之轉(zhuǎn)軸在瞬態(tài)電磁力和輸出脈動轉(zhuǎn)矩循環(huán)作用下，常常會產(chǎn)生疲勞損傷破壞，因此，永磁同步電機的可靠性問題倍受人們的關(guān)切。

在疲勞特性分析方面，很多專家學者對影響主軸的疲勞特性的因素進行了深入的研究。文獻[3]對電機的軸進行了靜力學的疲勞分析，提出了在軸肩處增大圓角提高主軸疲勞強度的方法。文獻[4]應用Ansys軟件中的疲勞分析模塊對變速器主軸進行了結(jié)構(gòu)疲勞分析，校核了軸的疲勞壽命滿足設計壽命要求。文獻[5]通過試驗的方法證明了粗糙度對壽命的影響，得出粗糙度的增大，疲勞壽命會相應的減小的結(jié)論。

永磁同步電機內(nèi)部定子繞組產(chǎn)生基波磁場和永磁體產(chǎn)生的諧波磁場相互作用，產(chǎn)生不斷變化的徑向電磁力和輸出轉(zhuǎn)矩。徑向電磁力引起主軸的彎曲振動，輸出轉(zhuǎn)矩引起主軸的扭轉(zhuǎn)振動，這都會影響主軸的壽命。但是在實際模擬仿真電機主軸受力時，瞬態(tài)電磁力一直被忽視。本文首先通過解析法和數(shù)值法，計算主軸受到的瞬態(tài)徑向電磁力和輸出脈動轉(zhuǎn)矩，然后再利用疲勞分析軟件對永磁同步電機的主軸進行了分析。并根據(jù)分析結(jié)果，提出相應的改進措施。

1 永磁同步電機電磁場分析

1.1 理論計算

本文在描述永磁同步電機作用狀態(tài)時，采用了有限元分析軟件進行分析的方法。首先通過數(shù)值方法對電磁場進行分析，再利用麥克斯韋應力法的后處理方式，即可獲取永磁同步電機的徑向電磁力和切向電磁力。

用矢量磁位A描述的電機瞬態(tài)電磁場方程和電路方程可以表示為：

通過對上述方程進行離散求解，就可以獲得不同時刻永磁同步電機二維界面的磁力線分布，為下一步求解電磁力的分布打下基礎。

由麥克斯韋應力關(guān)于磁質(zhì)上的合力方程：

式中，T是一個二階張量；S為包圍磁質(zhì)的任意一個閉合曲面；S通常設置在磁質(zhì)周圍的空氣中。

對于二維問題，這個積分面就簡化為一條閉合曲線，合力F可用如下公式表示[7]：

式中，t，n為積分路徑的切向和法向單位矢量；ft為切向電磁力密度；fr為徑向電磁力密度。其中：

其中，x，y是相應單元的坐標。

1.2 數(shù)值計算

本文以一臺8極36槽的表面式永磁同步電機為例，其額定功率為25kW，額定轉(zhuǎn)速為2000r/min，該電機的橫截面幾何結(jié)構(gòu)見圖1。

1.定子 2.永磁體 3.主軸 4.轉(zhuǎn)子 5.氣隙圖1 永磁同步電機橫截面幾何結(jié)構(gòu)

本文建立了永磁同步電機的模型，并利用時步有限元分析法，對永磁同步電機模型進行仿真分析。通過仿真結(jié)果分析，得到了永磁同步電機的磁力線分布和電磁場內(nèi)各節(jié)點的磁密值。圖2為永磁同步電機二維磁力線趨勢圖，從中可以看到磁力線分布。圖3為永磁同步電機二維磁密分布云圖。

圖2 永磁同步電機二維磁力線趨勢圖

圖3 永磁同步電機二維磁密分布云圖

永磁同步電機的氣隙中由定子繞組產(chǎn)生的基波磁場和永磁體產(chǎn)生的諧波磁場相互作用，產(chǎn)生不斷變化的徑向電磁力密度和輸出轉(zhuǎn)矩。用數(shù)值法對永磁同步電機在正弦波供電與變頻器供電條件下的氣隙磁場進行分析，并對計算結(jié)果進行頻譜分析。

圖4 徑向電磁力密度隨時間變化的曲線及其頻譜圖

圖4為永磁同步電機轉(zhuǎn)子徑向電磁力密度隨時間變化的曲線及其頻譜分析圖。從圖中可以看出，隨著頻率的增加，徑向磁通密度的幅值卻逐漸減小。徑向電磁力密度的峰值對應的頻率為極數(shù)與旋轉(zhuǎn)基波頻率的乘積，f=2kpf0=2k×4×2000/60=266.7Hz，533.3Hz，800Hz…，k=1，2，3…。

圖5為永磁同步電機主軸輸出轉(zhuǎn)矩及其頻譜圖，由圖中可以得知輸出轉(zhuǎn)矩會隨著頻率的增加而逐漸減小，但下降率比較小，相對平穩(wěn)。輸出轉(zhuǎn)矩的頻率峰值出現(xiàn)的次數(shù)較多，其頻率值分別為800Hz，1600Hz，2400Hz等。

圖5 輸出轉(zhuǎn)矩及其頻譜圖

徑向電磁力是永磁同步電機主軸產(chǎn)生彎曲振動的主要原因。輸出轉(zhuǎn)矩是永磁同步電機產(chǎn)生扭轉(zhuǎn)振動的主要原因。因此，在進行電機結(jié)構(gòu)設計時，要盡可能的減少電磁力諧波和轉(zhuǎn)矩脈動對電機主軸造成的影響。

2 永磁同步電機主軸的數(shù)值分析

本文的研究對象為永磁同步電機的主軸，其疲勞壽命分析基本思路為：依據(jù)產(chǎn)品的二維圖紙，在三維軟件Pro/E中建立主軸的實體模型。根據(jù)數(shù)值分析基本理論，先利用Ansys Workbench軟件對模型進行定義材料，網(wǎng)格劃分和受力分析。將靜力分析結(jié)果傳給nCode Design-Life。綜合主軸的實際工況及材料的S-N曲線預估永磁同步電機主軸的疲勞壽命。

把永磁同步電機主軸的三維模型導入有限元軟件Ansys Workbench中，設置永磁同步電機主軸的材料為45#鋼。并添加材料屬性，如表1所示。

表1 45#鋼材料屬性

本文采用四面體自由劃分網(wǎng)格模型對主軸進行劃分網(wǎng)格，劃分完后得到永磁同步電機主軸節(jié)點數(shù)為78437個。單元類型為四面體，單元的總數(shù)為51546個。

網(wǎng)格劃分完成后，根據(jù)主軸在永磁同步電機中的實際情況，在軸承位置施加軸承約束，將上一節(jié)計算的徑向電磁力和電磁轉(zhuǎn)矩導入Ansys Workbench中。

通過Ansys Workbench進行有限元計算，得出主軸靜態(tài)分析應力云圖，如圖6所示。初步確定主軸容易發(fā)生斷裂的位置可能在輸出端的第一變截面上。

圖6 永磁同步電機主軸的應力云圖

主軸材料的屈服極限為310MPa，靜強度安全系數(shù)取1.5，得到材料的許用應力為206.7MPa。從應力云圖中可見永磁同步電機主軸的最大應力值為102.9MPa，低于材料的許用應力。表明該主軸有足夠的剛度和強度儲備，在壽命周期內(nèi)不會發(fā)生疲勞損壞。

3 永磁同步電機主軸的疲勞可靠性分析

3.1 確定疲勞破壞類型

疲勞壽命指的是結(jié)構(gòu)或機械直至破壞所作用的循環(huán)載荷的次數(shù)或時間[8]。近年來，國內(nèi)外研究工作人員通過采用疲勞分析軟件MSC.Fatigue和不同的疲勞設計方法相結(jié)合對零件進行疲勞分析[9-11]，而本文主要采用nCode Design-Life進行分析。疲勞破壞類型通常分為高周疲勞和低周疲勞。根據(jù)Ansys Workbench靜力分析結(jié)果，確定了該永磁同步電機主軸屬于高周疲勞類型。

3.2 確定材料的S-N曲線

高精度的疲勞壽命預測離不開真實的材料參數(shù)輸入。選擇電機主軸的材料為45#鋼，材料屬性參見表1。將材料屬性輸入nCode Design-Life中，得出了主軸材料的估計S-N曲線，如圖7所示。因為實驗用的是標準試棒，而實際分析中零件的形狀和真實形狀差別較大，因此需要對S-N曲線進行修正，在nCode Design-Life具有對S-N曲線自動修正的功能。對45#鋼材料修正后的S-N曲線，如圖8所示。

圖7 45#鋼的估計S-N曲線

圖8 45#鋼修正后的S-N曲線

3.3 疲勞可靠性分析

將Ansys Workbench中生成的靜力分析結(jié)果文件導入到nCode Design-Life疲勞分析軟件中，如圖9所示。設置好各個參數(shù)，并進行求解，得到主軸疲勞損傷云圖，如圖10所示。

圖9 可靠性分析界面

圖10 永磁同步電機主軸疲勞損傷云圖

分析以上圖可知，在nCode Design-Life中計算得到主軸輸出端第一個變截面位為危險截面，在電磁力和輸出轉(zhuǎn)矩循環(huán)作用下，存在斷軸隱患。因此，需要采取措施改進主軸的結(jié)構(gòu)，減少應力集中，提高主軸的壽命。

針對主軸結(jié)構(gòu)的應力集中問題，經(jīng)過綜合分析，提出了采用了凹圓角結(jié)構(gòu)，既解決了電機主軸安裝軸承時軸肩倒角與軸承的干涉問題，又解決了產(chǎn)生應力集中的現(xiàn)象。凹圓角的布置如圖11所示。

圖11 凹圓角示意圖

4 結(jié)論

本文針對永磁同步電機主軸發(fā)生疲勞斷裂的情況進行了電磁場分析和疲勞特性分析，主要得到了以下結(jié)論：

(1)永磁同步電機主軸斷裂主要由瞬態(tài)徑向電磁力引起的彎曲振動和電磁轉(zhuǎn)矩脈動作用下的高周疲勞引起的。

(2)通過nCode Design-Life軟件得到的仿真結(jié)果與實際在主軸輸出端第一變截面處發(fā)生斷裂的位置相符，在該處采用凹圓角的結(jié)構(gòu)，并降低該處的表面粗糙度值，可有效的減少應力集中，提高主軸的壽命，為主軸的結(jié)構(gòu)設計提供了參考。

(3)由于加工裝配引起的定轉(zhuǎn)子不對中，與負載連接的聯(lián)軸器安裝偏心可產(chǎn)生單邊磁拉力；主軸材料特性和熱處理的質(zhì)量也是造成主軸的疲勞破壞的主要原因。

[1] 李靜，程小華. 永磁同步電機的發(fā)展趨勢[J]. 防爆電機，2009，5(44):1-4.

[2] 唐麗嬋，齊亮. 永磁同步電機的應用現(xiàn)狀與發(fā)展趨勢[J]. 裝備機械，2011(1)：7-12.

[3] 曾翠華，廖海平．電機軸優(yōu)化設計[J]．機電設備，2005(2):44-47.

[4] 吳勝軍，徐有良．輕型汽車變速箱第一軸的疲勞分析[J].機械設計與制造，2008(5):111-112.

[5] 魚濱濤，佟巍峰.表面粗糙度對國產(chǎn)316LN鋼低周疲勞性能的影響[J]. 原子能科學技術(shù)，2015,49(9):1660-1665.

[6] 于慎波，姜菲菲，王輝，等.永磁同步電主軸分數(shù)槽電機的徑向電磁力分析[J]. 組合機床與自動化加工技術(shù)，2014(6):15-18.

[7] 張文敏，郝清亮.考慮變頻器供電諧波電流作用下的永磁電機電磁力計算與分析[J]. 大電機技術(shù)2011(4):23-29.

[8] 姚衛(wèi)星，結(jié)構(gòu)疲勞壽命分析[M]. 北京:國防工業(yè)出版社,2003.

[9] J Schijve. Fatigue of Structures and Materials in the 20th Century and the State of The Art[J]. International Journal of Fatigue, 2003(25):679-702.

[10] 李成林，宋莎莎，韓振南. 基于nCode Design-Life的某車架疲勞可靠性分析[J]. 圖學學報，2014,3(1):42-45.

[11] Sigmund K A. Fatigue assessment of aluminum automotive structure [D]. Norwegian:Norwegian University of science and Technology,2002.

(編輯 李秀敏)

Fatigue Characteristic Analysis of the Rotor of Permanent Magnet Synchronous Motor Based on nCode Design-Life

YU Shen-bo，LI Hong，ZHAI Feng-jun，ZHONG Shuang-shuang，XIA Peng-peng

(School of Mechanical Engineering, Shenyang University of Technology, Shenyang 110870,China)

This paper researches the phenomenon of the crack and the break appearing on the rotor of Permanent Magnet Synchronous Motor (PMSM).Firstly, 2D electromagnetic force on the rotor is analyzed and computed by analytical and numerical method. Load spectrum on the rotor is defined. The fatigue characteristics of the rotor are analyzed by using S-N analysis method of nCode Design-Life. The stress concentration appears on the first section on the output side of the rotor from the fatigue damage distributing. The danger of broken shaft is hidden under the cycle action of transient electromagnetic force and output ripple torque. Based on this, the paper introduces the machining method of concave fillet at the variable section, which can effectively reduce the stress concentration and increase the life span of rotor. The achieved results can provide the reference of reliability design of PMSM.

fatigue analysis; permanent magnet synchronous motor; electromagnetic force

1001-2265(2017)01-0050-04

10.13462/j.cnki.mmtamt.2017.01.014

2016-09-21;

2016-11-02

國家自然科學基金項目(51175350)；沈陽市科技計劃項目(F15-199-1-13)

于慎波(1958—)，男，沈陽人，沈陽工業(yè)大學教授，博士，研究方向為電機噪聲與振動抑制技術(shù)，(E-mail)yushenbo@126.com；通訊作者：李紅(1989—)，女，山東淄博人，沈陽工業(yè)大學碩士研究生，研究方向為機械設計制造及其自動化，(E-mail)274709167@qq.com。

TH166;TG506

A