永磁同步電動機結構的模態研究*

李 娉 馬永龍

(1.湖北三江航天險峰電子信息有限公司 武漢 432000)(2.海軍駐中南地區光電系統軍事代表室 武漢 430074)

永磁同步電動機結構的模態研究*

李 娉1馬永龍2

(1.湖北三江航天險峰電子信息有限公司 武漢 432000)(2.海軍駐中南地區光電系統軍事代表室 武漢 430074)

論文先簡要介紹模態分析理論,再運用機電類比的方法推導定子固有頻率的計算方法,并用ANSYS軟件對電機定子建立的三維有限元模型進行模態分析,計算出它們的固有頻率及其振型,避免和電磁力的頻率及次數發生共振,從而降低永磁同步電動機的電磁振動和噪聲。

永磁同步電動機; 模態分析; ANSYS

Class Number TM359.4

1 引言

永磁同步電動機的電磁振動與噪聲不僅與電機中氣隙磁場所產生的電磁激振力波頻率、幅值和力波階數有關,而且還與電機結構的固有頻率的大小及分布有著密切的關系。通過對電機定子結構的模態分析,可以確定電機電樞結構的固有頻率和振型,就可預測電機在電磁力的作用下實際振動響應。

模態分析技術從20世紀60年代后期發展至今已經有40多年的歷史了。模態分析的經典性定義是:將線性定常系統振動微分方程組中的物理坐標變換為模態坐標,使方程組解藕,成為一組以模態坐標及模態參數描述的獨立方程,以便求出系統的模態參數[1]。振動模態是振動系統的各點以特定的頻率作簡諧振動時,表示波節和波腹的振動形態或與其相應的衰減振動狀態;固有模態是自由振動時系統的振動模態,與系統的各個自由度相對應,常簡稱為模態[2]。模態是機械及結構的固有振動特性。每一階模態具有特定的固有頻率、阻尼比和模態振型。這些模態參數可以由有限元計算或試驗分析取得?；诰€性疊加原理,一個復雜的振動系統可以分解為許多模態的疊加。這樣一個分析過程為模態分析[3]。

永磁同步電動機結構的模態分析是判斷電動機定子結構是否發生共振的重要手段,電動機的電磁噪聲和振動主要是通過定子、殼體向外輻射。本節主要是進行電動機結構的模態分析,由于電動機結構的幾何特性決定著定子沿軸向的剛度是變化的,所以采用三維結構有限元模型來分析電動機定子及機殼結構的模態。

2 電機結構的模態分析理論

永磁同步電動機模態分析的關鍵在于找到模態振型矩陣,以構成模態坐標系統,并求得響應量在這一坐標系統中的坐標,具體過程如下[4]:

1) 尋找各階模態振型,組成坐標變換矩陣{ψr}=[{ψ1} {ψ2} {ψ3} … {ψN}];

2) 進行坐標變換{X}={ψ}{q};

3) 求取模態坐標qr,r=1,2,3,…,N,從而得到{q};

4) 建立響應計算模型—模態模型;運用模態模型,便可計算在實際激勵作用下的結構的運動,包括位移、速度、加速度,乃至應力應變。

2.1 質點運動方程

通過牛頓第二定律可以得到質點的運動方程:

(1)

進行拉普拉斯變換后上式可寫為

s2mX(s)+srX(s)+kX(s)=F(s)

兩邊同除以m可變換為

則質點的位移響應為

質點的速度響應為

質點的加速度響應為

因此可得出由力到加速度的傳遞函數為

永磁同步電動機定子振動中,除了有電磁力引起的強迫振動,還有自由振動。下面分析永磁同步電動機的自由振動特性。令系統的運動方程(1)的右邊為0,則可得系統的自由振動方程為

該方程的通解形式為

式中,C1、C2為由系統初始條件決定的常數。

故運動方程的解可表示為

對于無阻尼運動,ζ=0;對于欠阻尼運動,ζ<1;對于臨界阻尼運動,ζ=1對于過阻尼運動,ζ>1。在實際的振動系統中,一般都是欠阻尼運動,否則就不會形成振動了。則對于欠阻尼情況,ζ2?1,運動方程的解可寫為

2.2 永磁同步電動機振動特性的機電類比分析

圖1 電機定子機電類比雙環型模型

永磁同步電動機一般選用同容量的異步電機機座,定子鐵心與機殼內圓采用過盈配合,無連接筋。設定子鐵心軛的剛度為K1,質量為m1,機殼的剛度為K2,質量為m2,不考慮繞組、定子凸極的影響,散熱筋、底腳、接線盒的質量歸入機殼質量;并設各部分的阻尼系數分別為rm1、rm2(假設各部分阻尼比為粘性,且等于鋼或鑄鐵的阻尼比),則其物理模型、機械網絡圖、等效電路圖如圖1所示[5～6]。

圖中p0為作用在定子鐵心內表面上單位面積力的幅值,pn為等效的集中力,即:

Pn=2πRilaP0

式中,Ri為定子磁軛的內半徑,la為定子鐵心軸向長度。

由機電類比法[7]可得對應角頻率ω的總機械阻抗為

則定子振動位移的幅值為

(2)

式中,鐵心及機殼的剛度和質量為[3]

(3)

(4)

設定子鐵芯、機殼的阻尼比ζ1、ζ2為粘性,則阻尼系數為

令機械阻抗的虛部為零,可求得定子的固有頻率為

(5)

共振時,定子振動位移的振幅為

將式(3)、(4)代入式(2)、(5),可得電機強迫振動幅值和固有頻率為

(6)

(7)

由式(6)～(7)可得永磁電機軸向階數為零的平面振動特性[6],即:

1) 振動幅值與力波幅值成正比。因此減小力波幅值是控制振動的主要措施之一;

2) 隨著振形階數n的增大,對應的固有頻率近似以n2的比例增大,而振幅則近似以n4的比例減小,因此,階數較低的力波引起的振動應是防范的重點;

3) 定子直徑越大,固有頻率越低,亦即當激勵力的振幅、頻率均相等時,幾何尺寸越大的電機,振動幅值越大,因而磁噪聲級亦越大;

4) 振動大小與振動系統的阻抗Zm關系很大,阻抗表明振動系統接近諧振的程度。諧振時振動甚至在小的電磁激振力時都可能很大。因此,使固有頻率避開主要的力波頻率是減小振動的主要措施之一;

5) 共振時激振力波與振動速度同相位,這時振動幅值完全由阻尼限制。

3 永磁同步電動機定子模態分析

目前,理論計算固有頻率的方法有兩大類:一類是解析解計算方法,典型的是機電類比法,該方法可以得到固有頻率的解析表達式,但是計算精度較差;另一類是能量法,它有兩種解法,一種是傅里葉級數解法,另一種是有限元法。一般情況下,兩種解法都不能得到解析解,而只能得到數值解。在定子結構對稱時,傅立葉級數的求解精度可滿足一般工程上的要求。有限元解法可以考慮定子結構的不規則性,其計算精度較高。由于電機結構復雜,用解析法計算雖然簡便,但計算誤差較大,因此本文采用有限元方法仿真計算了電機的模態,研究了電機的定子鐵心的模態。

用戶使用ANSYS的模態分析來決定一個結構或者機器部件的振動頻率(固有頻率和振型)。ANSYS產品家族的模態分析是線性分析。任何非線性特性,如塑性和接觸(間隙)單元,即使定義了也將被忽略。

需要記住以下兩個要點[8]:

1) 模態分析中只有線性行為是有效的,如果指定了非線性單元,它們將被當作是線性的。例如,如果分析中包含了接觸單元,則系統取其初始狀態的剛度值并且不再改變此剛度值。

2) 必須指定彈性模量EX(或某種形式的剛度)和密度DENS(或者某種形式的質量)。材料性質可以是線性的或非線性的、各向同性或正交各向異性的、恒定的或與溫度有關的,非線性特性將被忽略。用戶必須對某些指定的單元(COMBIN7,COMBIN14,COMBIN37)進行實常數的定義。

3.1 基于能量法的電機振動模態的有限元分析模型

作為有限元法的一種,能量法的基礎是拉格朗日方程。在能量法計算中,一般用廣義坐標來描述系統,將系統的動能和勢能等表示為廣義坐標及其導數的函數[9]。系統的拉格朗日運動方程可表示為

其中,L=T-U為拉格朗日函數,T為系統動能,U為系統勢能,Fi為系統非保守廣義力,qi為廣義坐標,i=1,2,3…。

在有限元解法中,單元的動能為

勢能為

單元的剛度矩陣為

依據Hamilton原理和應力—應變、應變—位移的關系[10],可導出電機離散后單元的運動方程為

將各個離散的單元方程組聯合成總的系統方程組,經過整理則可得出系統標準的運動方程為

(8)

([K]-ω2[M]){u}={0}

(9)

根據線性方程組理論有非零解的充分必要條件為

|[K]-ω2[M]|=0

(10)

解出滿足方程(9)、(10)的頻率ω2和對應的非零解向量{ur}(r=1,2,3…),則ωr,{ur}分別為電機結構的固有頻率和固有振型。

有限元法ANSYS軟件提供了子空間法、分塊Lanczos法、Power Dynamics法、縮減法、非對稱法、阻尼法和QR阻尼法七種模態提取方法。大多數分析都選用前四種求解方法,非對稱法和阻尼法只在特殊情況下會用到。本文主要采用分塊Lanczos方法對問題展開求解。此方法應用于提取大模型的多階模態,經常應用在具有實體單元或殼單元的模型中,特別是在模型中包含形狀較差的實體和外殼單元中,它的計算速度很快,但要以消耗內存為代價[11]。

分塊Lanczos法的特征值求解器采用缺省求解器,它采用Lanczos算法,是由一組向量來實現Lanczos遞歸運算。這種方法和子空間法一樣精確,但速度更快。無論EQSLV命令指定過何種求解器進行求解,分塊法都將自動采用稀疏矩陣方程求解器。計算某系統特征值譜所包含一定范圍的固有頻率時,采用這種方法特別有效。計算時,求解從頻率譜中間位置到高頻端范圍內的求解收斂速度和求解低階頻率時基本上一樣快。因此,當采用頻移頻率(FREQB)來提取從FREQB(起始頻率)的模態時,提取大于FREQB的n階模態和提取n階低頻模態的速度基本相同。在本節的計算中均采用了這種方法。

3.2 定子鐵芯的模態分析

永磁同步電動機定子鐵芯的彈性模量通常取2.05～2.15×1011N/m2,密度為7305kg/m2,泊松比是0.3。采用有限元軟件進行模態仿真時,模型選用的單元類型為PLANE82,分析類型選用分塊的蘭索斯(Block Lanczos)法,求解的模態數選擇為40,擴展的模態數選擇為40。圖2是定子鐵芯三維有限元模型,圖3是模態分析后的固有頻率及振型。

圖2 定子鐵芯三維有限元模型

(a) 二階振型

(b) 三階振型

(c) 四階振型

(d) 五階振型圖3 定子鐵芯模態分析的固有頻率及振型表1 永磁同步電動機定子三維模型的振型和固有頻率

振型固有頻率振型固有頻率二階1618Hz 1620Hz2251Hz 2253Hz四階5300Hz 5317Hz6225Hz 6228Hz7536Hz 7603Hz三階3791Hz 3805Hz4602Hz 4612Hz4886Hz 4893Hz五階8317Hz 8431Hz9542Hz 9553Hz1174Hz1183Hz

從表1中可以看出三維定子模型更加接近實際模型,可以得到豐富的模態和固有頻率信息。有些頻率相近是由于電機定子結構是對稱的,會出現振型和頻率相同但相位不同的情況;考慮永磁同步電動機定子結構的固有模態,有利于避免電磁共振,減少振動幅值,從而有利于降低電磁振動和噪聲。

4 結語

本文運用機電類比的方法推導出了定子固有頻率的計算方法。利用有限元分析軟件ANSYS建立了永磁同步電動機樣機定子的三維有限元模型,計算得到了樣機的振型及振動固有頻率的關系,從而為避免共振起到了預見性作用,為高轉速和低噪聲的永磁同步電動機定子設計提供理論參考。

[1] 袁安富,陳俊.ANSYS在模態分析中的應用[J].中國制造業信息化,2007,36(11):42-44.

[2] S. J. 楊著.低噪聲電動機[M].北京:科學出版社出版,1985.

[3] 張培強,馮文琴.現代實驗模態分析的基本原理[J].武漢水運工程學院學報,1984,(4):99-115.

[4] 李德葆.振動模態分析及其應用[M].北京:宇航出版社,1989:57.

[5] 陳永校,諸自強,應善成.電機噪聲的分析與控制[M].杭州:浙江大學出版社,1987:1-200.

[6] 孫劍波.開關磁阻電機的減振降噪和低轉矩[D].武漢:華中科技大學,2005.

[7] 王宏華,王治平,江泉.開關磁阻電動機固有頻率解析計算[J].中國電機工程學報,2005,25(12):133-135.

[8] 王建江,胡仁喜,劉英林,等.ANSYS11.0結構與熱力學有限元分析實例指導教程[M].北京:機械工業出版社,2008,3.

[9] 盛劍霓.工程電磁場數值分析[M].西安:西安交通大學出版社,1991.

[10] 吳建華.基于物理模型開關磁阻電機定子模態和固有頻率的研究[J].中國電機工程學報,2004,24(8):109-114.

[11] 佟寧澤.大中型感應電機電磁力及定子振動分析[D].沈陽:沈陽工業大學,2006.

Model Analysis of Permanent Magnet Synchronous Motor

LI Ping1MA Yonglong2

(1. Hubei Sanjiang Aerospace Xianfeng Electronic Information Co., Ltd., Wuhan 432000) (2. Navy Representative Office in EO System of Central-South Area, Wuhan 430074)

The paper introduces the theory of the model analysis. The electromechanical analogy method is used to derive the calculation method of the natural frequency of the stator, and the ANSYS software is used to set up motor stator three-dimensional finite element model and model analysis. The natural frequency and vibration mode are calculated to avoid resonance of frequencies and force with electromagnetic force then reduce electromagnetic vibration and noise of PMSM.

permanent magnet synchronous motor, model analysis, ANSYS

2013年8月5日,

2013年9月27日

李娉,女,工程師,研究方向:電子信息、技術改造等。馬永龍,男,碩士,工程師,研究方向:艦船電子。

TM359.4

10.3969/j.issn1672-9730.2014.02.041