郭嘉成，王 森，趙陽陽，閆鴻魁，呂宗樞

(沈陽工程學(xué)院，沈陽 110136)

電動飛機(jī)主推進(jìn)永磁同步電動機(jī)振動特性研究

郭嘉成，王 森，趙陽陽，閆鴻魁，呂宗樞

(沈陽工程學(xué)院，沈陽 110136)

分析了電動飛機(jī)主推進(jìn)永磁同步電動機(jī)振動產(chǎn)生的原因，給出了不同階數(shù)力波作用下固有頻率的表達(dá)式，利用有限元仿真軟件分別對永磁同步電動機(jī)的定子與機(jī)殼和整機(jī)進(jìn)行振動特性的仿真分析。通過分析可知，優(yōu)化電機(jī)結(jié)構(gòu)能夠抑制電機(jī)振動，從而保證電動飛機(jī)動力系統(tǒng)的安全性。

振動；徑向電磁力；永磁同步電動機(jī)；有限元

0 引 言

飛機(jī)要求在空中飛行安全、可靠，主推進(jìn)電機(jī)的性能直接影響其穩(wěn)定性。傳統(tǒng)飛機(jī)采用內(nèi)燃機(jī)作為主推進(jìn)系統(tǒng)，內(nèi)燃機(jī)燃燒化學(xué)燃料，不僅工作效率低，而且污染嚴(yán)重，不符合當(dāng)前國家大力提倡節(jié)能減排的要求。近些年迅速發(fā)展起來的永磁電機(jī)有著體積小、運行穩(wěn)定、效率高、節(jié)能環(huán)保等眾多優(yōu)點，現(xiàn)已在航空航天、醫(yī)療器械、精密電子儀器、能源電力等眾多領(lǐng)域廣泛地應(yīng)用。永磁同步電動機(jī)用在飛機(jī)的主推進(jìn)系統(tǒng)上，能夠使飛機(jī)在飛行中更加高效、可靠。然而，由于電動飛機(jī)主推進(jìn)永磁同步電動機(jī)在工作時會產(chǎn)生振動，對飛機(jī)的穩(wěn)定性會造成一定影響，所以必須對電機(jī)的振動特性作研究分析。

隨著科技的發(fā)展，近些年來，很多專家學(xué)者針對電機(jī)振動問題展開了一系列研究，并且取得了一定的研究成果。文獻(xiàn)[1]對具有分?jǐn)?shù)槽集中繞組結(jié)構(gòu)的永磁電機(jī)的振動特性進(jìn)行了探討研究，通過對氣隙中徑向電磁力的仿真分析，得出相對于整數(shù)槽電機(jī)，具有分?jǐn)?shù)槽集中繞組結(jié)構(gòu)的永磁電機(jī)內(nèi)部產(chǎn)生的徑向電磁力波，出現(xiàn)了使電機(jī)振動加強(qiáng)的低階空間諧波。文獻(xiàn)[2]則分析探討了無刷直流電動機(jī)的振動受著齒槽定位力矩的影響，解析了消除定位轉(zhuǎn)矩的機(jī)理并借助簡單的模型分析由于轉(zhuǎn)子永磁體不對稱引起的定位轉(zhuǎn)矩的影響和由于永磁體磁動勢幅值和相位誤差造成的影響，最后進(jìn)行實驗驗證。文獻(xiàn)[3]根據(jù)Fourier變換分析方法計算了籠型異步電動機(jī)的電磁力波，得出了導(dǎo)致電機(jī)振動的主要力波組成成分。文獻(xiàn)[4]研究了無刷直流電機(jī)的振動與噪聲的問題，得到電機(jī)相電動勢矢量圖和繞組構(gòu)成，找到了電機(jī)徑向不平衡電磁力產(chǎn)生的原因。文獻(xiàn)[5]則利用解析算法詳細(xì)分析計算出了永磁同步電動機(jī)的電磁場。文獻(xiàn)[6]是以磁-固耦合的方法為依據(jù)對異步電動機(jī)的振動進(jìn)行了分析研究，利用有限元分析軟件對異步電機(jī)的變電磁力做了頻譜分析和校正，通過優(yōu)化電機(jī)定、轉(zhuǎn)子齒槽數(shù)可降低電機(jī)的電磁振動。

本文從原理上定性討論了電動飛機(jī)主推進(jìn)永磁電機(jī)振動問題，考慮到飛機(jī)在實際飛行過程中的狀況，利用有限元仿真軟件對永磁同步電動機(jī)進(jìn)行建模并根據(jù)仿真結(jié)果分析了永磁同步電動機(jī)的振動特性。

1 振動特性分析

1.1振動原因分析

(1) 永磁電機(jī)定、轉(zhuǎn)子氣隙磁場的徑向分量相互作用所產(chǎn)生的電磁力波是引起電機(jī)電磁振動的主要原因。一些學(xué)者所作大量的研究表明：雖然電磁力是直接作用在齒和磁極上，但定子與機(jī)殼主要振動的根源并不是齒和磁極的振動，而是定子軛的振動。此外，定子上的繞組因為受到電磁力的作用會產(chǎn)生繞組系統(tǒng)的工頻或者倍頻振動。

在利用解析算法分析計算電機(jī)氣隙磁場時，一般不考慮鐵心磁阻和飽和的影響，所以得到永磁同步電動機(jī)徑向氣隙磁密的解析表達(dá)式[7]：

br(θ,t)=f(θ,t)·λ(θ,t)

(1)

式中：f(θ,t)為氣隙磁動勢；λ(θ,t)為氣隙比磁導(dǎo)。

電機(jī)氣隙中單位面積電磁力的瞬時值可表示：

(2)

式中：br(θ,t)為電機(jī)徑向氣隙磁密；bt(θ,t)為電機(jī)切向氣隙磁密；θ為機(jī)械角位移；μ0為空氣磁導(dǎo)率；t為時間。

(2) 由于永磁同步電動機(jī)在轉(zhuǎn)子上開有很多槽孔，這樣就難免會使轉(zhuǎn)子出現(xiàn)軸不對稱情況。同時，在加工和裝配時的誤差、磨損等諸多不利因素加上電機(jī)運行過程中引起的熱不對稱，均會導(dǎo)致轉(zhuǎn)子相對于定子偏心，從而產(chǎn)生不平衡磁拉力，使電動機(jī)產(chǎn)生振動。

由轉(zhuǎn)子相對于定子偏心產(chǎn)生的不平衡磁拉力：

F0=k0e0

(3)

式中：k0為不平衡磁拉力剛度；e0為初始偏心。

(3) 永磁電機(jī)機(jī)座產(chǎn)生振動與軸承的磨損、老化會造成電機(jī)的振動。在電機(jī)實際運行當(dāng)中，機(jī)座產(chǎn)生振動的原因可總結(jié)為以下兩方面：一方面是轉(zhuǎn)子出現(xiàn)振動產(chǎn)生的激振力作用到電機(jī)機(jī)座上產(chǎn)生的振動，另一方面則是電機(jī)鐵心所產(chǎn)生的電磁振動經(jīng)由鐵心與機(jī)座的連接部位傳遞到機(jī)座上而引起電機(jī)機(jī)座的振動[8]。

(4) 此外，固有頻率是物體本身具有的一種屬性，當(dāng)電機(jī)內(nèi)部徑向電磁力的頻率等于電機(jī)本身的固有頻率時就會發(fā)生共振現(xiàn)象，從而產(chǎn)生劇烈的振動，對物體結(jié)構(gòu)有嚴(yán)重的危害。本文計算了各階力波對固有頻率的影響。

1.2振動模態(tài)的數(shù)學(xué)模型

振動模態(tài)的分析，即利用系統(tǒng)本身模態(tài)的正交性，用變換矩陣替代模態(tài)向量所組成的模態(tài)矩陣，基于線性變換法對系統(tǒng)的物理坐標(biāo)進(jìn)行轉(zhuǎn)換，使得振動系統(tǒng)原本互相耦合的運動方程組變換成一組互相獨立的振動模態(tài)方程。

依據(jù)Hamilton原理和應(yīng)力、應(yīng)變、位移之間的關(guān)系，并忽略阻尼的影響，可推導(dǎo)得出電機(jī)的運動方程[9]：

(4)

式中：K為剛度矩陣；R為阻尼矩陣；M為質(zhì)量矩陣；u為節(jié)點的位移矢量。若時間導(dǎo)數(shù)用jω代替(其中ω為角頻率)，那么電機(jī)的無阻尼振動模態(tài)計算可轉(zhuǎn)變?yōu)榍蠼夥匠探M特征值的問題。有：

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

(5)

依據(jù)線性代數(shù)中的理論，若使一組齊次線性方程組具有非零解，則充分必要條件：

(6)

1.3固有頻率計算

物體的固有頻率僅由物體自身性質(zhì)所決定，與外部條件無關(guān)，它由一組模態(tài)參數(shù)(固有頻率、模態(tài)振型、模態(tài)阻尼比等)定量描述，實際中可通過無阻尼自由振動方程計算得出。

計算時假設(shè)定子軛為圓柱形的殼體，一系列在時間上呈周期性變化的r階力波沿圓周均勻的分布在定子軛上。利用有限元軟件建立定子軛的簡單模型，并進(jìn)行三維模態(tài)仿真，得到r≤4階模態(tài)和振型，如圖1所示。

圖1 定子軛振動的空間型式

本文分以下3種情況加以討論：

(1) 力波階數(shù)r=0固有振動頻率：

(7)

(2) 力波階數(shù)r=1固有振動頻率：

(8)

(3) 力波階數(shù)r≥2固有振動頻率：

(9)

式中：E為彈性模量；m為定子軛上平均半徑處的單位表面質(zhì)量；hf為定子軛的徑向厚度；Rf1為定子軛平均半徑。

由上述分析可知，在其他條件不變的下，定子軛的半徑越大，其固有振動頻率越低。

2 有限元仿真分析

本文分別對電機(jī)定子與機(jī)殼系統(tǒng)和整機(jī)進(jìn)行振動仿真分析。利用有限元仿真軟件，建立仿真模型，得到了仿真結(jié)果并分析電機(jī)振動特性。

2.1永磁電機(jī)定子與機(jī)殼的有限元仿真

通過對電機(jī)定子與機(jī)殼系統(tǒng)進(jìn)行有限元仿真，分別得到頻率為897 Hz下定子與機(jī)殼的振動仿真結(jié)果(如圖2所示)和永磁同步電動機(jī)定子與機(jī)殼的振動質(zhì)量參與系數(shù)表(如表1所示)。

圖2 897 Hz下定子與機(jī)殼的振動仿真結(jié)果圖

模數(shù)頻率f/HzX方向Y方向Z方向16760.171080.311646.9136×10-526800.311580.172215.5819×10-637983.8759×10-55.2094×10-60.53131048973.2264×10-61.2037×10-50.016255

由仿真結(jié)果可知，定子與機(jī)殼組成的系統(tǒng)整體振動幅度不大。機(jī)殼后端蓋振動幅度非常小，而機(jī)殼前端蓋靠近轉(zhuǎn)軸處的螺栓連接部分則振動的比較劇烈。在低階的情況下，定子與機(jī)殼軸向振動比較劇烈，而徑向振動幅度不大。在高階情況下，定子與機(jī)殼的徑向振動幅度變大，而軸向振動幅度變小。

此外，電機(jī)的機(jī)殼結(jié)構(gòu)、螺孔固定的位置、繞組排列的方式等因素對電機(jī)的振動模態(tài)及共振頻率都有一定影響。而氣隙大小、定子齒厚、定子齒和軛寬度等因素對固有頻率影響則很小。通過此部分的研究，選擇高強(qiáng)度材料，加大機(jī)殼厚度，加大螺栓連接處的厚度，使連接更加緊密，能夠有效解決定子與機(jī)殼振動對電動飛機(jī)主驅(qū)動永磁同步電動機(jī)穩(wěn)定性影響的問題。

2.2永磁電機(jī)整機(jī)的有限元仿真

通過對永磁電機(jī)的整機(jī)進(jìn)行有限元仿真，分別得到頻率為181 Hz整機(jī)振動仿真結(jié)果圖3和永磁同步電動機(jī)整機(jī)的振動質(zhì)量參與系數(shù)表2。

圖3 181 Hz整機(jī)振動仿真結(jié)果圖

模數(shù)頻率f/HzX方向Y方向Z方向11208.6461×10-52.4975×10-57.7882×10-521300.197680.134180.1611632400.00402650.014790.012431044200.00286170.112390.0071478

由仿真結(jié)果可知，相比于定子與機(jī)殼的振動，永磁同步電動機(jī)整機(jī)的振動幅度較小，且質(zhì)量參與系數(shù)較小，表明將定子與機(jī)殼和轉(zhuǎn)子結(jié)合到一起,起到了抑制振動的作用。由于飛機(jī)啟動過程中，其自身的振動頻率是逐漸上升的，通過結(jié)構(gòu)的優(yōu)化設(shè)計，使電機(jī)整機(jī)的共振頻率降低，這樣就能夠減小共振范圍，避免與飛機(jī)發(fā)生共振。

3 結(jié) 語

本文研究了電動飛機(jī)主推進(jìn)永磁同步電動機(jī)的振動問題。給出了引起電機(jī)振動的原因，利用有限元軟件仿真得到永磁同步電動機(jī)定子與機(jī)殼系統(tǒng)和整機(jī)的振動結(jié)果，得出以下結(jié)論：

(1) 電機(jī)定子與轉(zhuǎn)子間氣隙中產(chǎn)生的電磁力波是引起永磁同步電動機(jī)電磁振動的主要原因，定子與機(jī)殼的主要振動是由電磁力作用于定子軛產(chǎn)生的，轉(zhuǎn)子的主要振動來源是由于轉(zhuǎn)子偏心產(chǎn)生的不平衡磁拉力。此外，機(jī)座和軸承的振動對電機(jī)振動產(chǎn)生一定的影響。

(2) 定子與機(jī)殼整體振動幅度不大，選擇高強(qiáng)度材料，加大機(jī)殼厚度，加大螺栓連接處的厚度，使連接更加緊密，可有效抑制定子與機(jī)殼的振動。

(3) 整機(jī)振動相對于定子與機(jī)殼和轉(zhuǎn)子單獨振動幅度要小，表明兩個系統(tǒng)結(jié)合到一起起到了相互抑制振動的效果，通過結(jié)構(gòu)的優(yōu)化設(shè)計，能夠降低電動飛機(jī)主驅(qū)動永磁電機(jī)的振動強(qiáng)度，提高飛機(jī)飛行的穩(wěn)定性。

[1] 呂長朋.分?jǐn)?shù)槽集中繞組永磁電機(jī)振動特性研究[J].船電技術(shù)，2015,35(5)：52-55.

[2] 程樹康,葛新,高宏偉,等.分?jǐn)?shù)槽無刷直流電動機(jī)齒槽定位力矩的研究[J].中國電機(jī)工程學(xué)報，2008,28(21)：107-111.

[3] 王荀,邱阿瑞.籠型異步電動機(jī)徑向電磁力波的有限元計算[J].電工技術(shù)學(xué)報,2012,27(7):109-117.

[4] 陳磊，高宏偉，柴鳳，等.小型無刷直流電動機(jī)振動與噪聲的研究[J].中國電機(jī)工程學(xué)報，2006,26(24)：148-152.

[5] 趙鏡紅，張俊洪，方芳，等.徑向充磁圓筒永磁直線同步電機(jī)磁場和推力解析計算[J].電工技術(shù)學(xué)報，2011.26(7)：154-160.

[6] 王玎．基于有限元的異步電機(jī)電磁振動分析[J].振動與沖擊，2012,31(2)：140-144.

[7] 陳世坤．電機(jī)設(shè)計[M]．2版.北京：機(jī)械工業(yè)出版社，2004．

[8] 郝飛,賈非.電機(jī)振動的產(chǎn)生及控制管窺[J].科技向?qū)В?015, 21(4)：132-133.

[9] VEDRINES M,GASSMANN V,KNITTEL D.Moving web-tensi on determination by out-of-plane vibration measurements using a laser[J].IEEE Transactions on Instrumentation and Measurement,2009,58(1):207-213.

ResearchonVibrationCharacteristicsofPermanentMagnetSynchronousMotorofElectricAircraftMainDrive

GUOJia-cheng，WANGSen，ZHAOYang-yang，YANHong-kui，LüZong-shu

(Shenyang Institute of Engineering, Shenyang 110136, China)

The reason for vibration of permanent magnet synchronous motor was analyzed and the expression of natural frequency under different orders force wave was calculated. The finite element analysis software was used to simulate and analyze the vibration characteristics for the stator and casing and the whole motor. Based on the analytical results, the structure of the motor is optimized to reduce the influence of vibration, which ensures the security of electric aircraft power system.

vibration; radial electromagnetic force; permanent magnet synchronous motor; finite element

2016-02-02

TM351；TM341

：A

：1004-7018(2016)11-0037-03

郭嘉成(1988-)，男，碩士研究生，研究方向為永磁電機(jī)設(shè)計。