劉細(xì)平，曾 劍，徐 晨

(江西理工大學(xué)，贛州 341000)

電動(dòng)汽車(chē)用開(kāi)關(guān)磁阻電機(jī)溫度場(chǎng)分析研究

劉細(xì)平，曾 劍，徐 晨

(江西理工大學(xué)，贛州 341000)

對(duì)一臺(tái)電動(dòng)汽車(chē)用110 kW水冷式12/8極開(kāi)關(guān)磁阻電機(jī)工作時(shí)的溫度場(chǎng)進(jìn)行數(shù)值分析，為其在實(shí)際應(yīng)用中提供可靠理論依據(jù)。利用SolidWorks軟件建立了開(kāi)關(guān)磁阻電機(jī)三維實(shí)體模型；采用ANSYS Maxwell軟件計(jì)算了電機(jī)在額定狀態(tài)下的銅耗和鐵耗，確定了損耗熱源；借助有限元方法對(duì)電機(jī)三維磁熱耦合仿真模型進(jìn)行了數(shù)值分析計(jì)算，深入分析了額定運(yùn)行狀態(tài)時(shí)的電機(jī)溫度場(chǎng)分布規(guī)律。分析結(jié)果表明，在額定工況下，電機(jī)內(nèi)部溫度符合電機(jī)的熱設(shè)計(jì)要求，為該類(lèi)電機(jī)應(yīng)用在電動(dòng)汽車(chē)上提供了可靠的理論依據(jù)。

Maxwell；開(kāi)關(guān)磁阻電機(jī)；溫度場(chǎng)；耦合仿真

0 引 言

開(kāi)關(guān)磁阻電機(jī)(Switched Reluctance Motor，SRM)是一類(lèi)結(jié)構(gòu)簡(jiǎn)單的磁阻式電動(dòng)機(jī)，其定子極和轉(zhuǎn)子極均為凸極，轉(zhuǎn)子上無(wú)任何形式的繞組和永磁體。由于SRM具有轉(zhuǎn)子結(jié)構(gòu)簡(jiǎn)單可靠、耐高溫、成本低、調(diào)速范圍寬、啟動(dòng)轉(zhuǎn)矩大、響應(yīng)速度快、可控參數(shù)多、和效率高等優(yōu)點(diǎn)，一直以來(lái)受到國(guó)內(nèi)外專(zhuān)家學(xué)者的高度重視，為一種具有較好發(fā)展前景的電動(dòng)機(jī)，也是目前研究的熱點(diǎn)[1-2,9]。

近年來(lái)，有關(guān)電機(jī)多物理場(chǎng)計(jì)算的研究屢見(jiàn)不鮮，研究方法也層出不窮，主要的計(jì)算方法有：簡(jiǎn)化公式法，等效熱路法，有限元法[3,5]。其中有限元法以其計(jì)算精確和方便直觀的優(yōu)點(diǎn)，得到了國(guó)內(nèi)外學(xué)者的青睞。本文借助ANSYS有限元分析軟件對(duì)開(kāi)關(guān)磁阻電機(jī)進(jìn)行三維溫度場(chǎng)分析，計(jì)算電機(jī)各部件以及電機(jī)內(nèi)部各個(gè)點(diǎn)的溫度，分析研究電機(jī)溫度分布情況，以保證電機(jī)安全可靠運(yùn)行。

本文針對(duì)一臺(tái)110 kW 電動(dòng)汽車(chē)用3相12/8極開(kāi)關(guān)磁阻電機(jī)，利用SolidWorks軟件建立SRM電機(jī)三維實(shí)體分析模型，然后將三維實(shí)體模型導(dǎo)入至ANSYS Maxwell中分析電機(jī)額定工況下的損耗，再利用數(shù)值分析方法計(jì)算溫度分布規(guī)律；獲取了在自然水冷的情況下電機(jī)定子鐵心、轉(zhuǎn)子鐵心、定子繞組、機(jī)殼和冷卻水管的溫升與溫度分布情況。

1 電機(jī)模型與基本方程

1.1 電機(jī)結(jié)構(gòu)模型和基本假設(shè)

利用SolidWorks建立的12/8極開(kāi)關(guān)磁阻電機(jī)三維模型如圖1所示。其結(jié)構(gòu)簡(jiǎn)單，主要由轉(zhuǎn)子鐵心、定子鐵心、繞組以及轉(zhuǎn)軸組成。轉(zhuǎn)子上既無(wú)永磁體，也無(wú)繞組，可實(shí)現(xiàn)高速運(yùn)行；其中電機(jī)繞組采用集中式繞組，4根雙邊絕緣厚度為0.06 mm、直徑為1.12 mm的銅線繞組并繞在定子齒上；定子鐵心上貫穿12根水管供電機(jī)進(jìn)行水冷循環(huán)。

圖1 開(kāi)關(guān)磁阻電機(jī)結(jié)構(gòu)示意圖

1.2 溫度場(chǎng)求解基本方程及邊界條件

根據(jù)熱學(xué)的熱傳導(dǎo)方程式可知[3]，當(dāng)處于直角坐標(biāo)系時(shí)，電機(jī)在穩(wěn)態(tài)時(shí)的溫度場(chǎng)結(jié)果求解分析可近似歸納為如下：

(1)

式中：λx，λy，λz分別為x,y,z方向上的導(dǎo)熱系數(shù)；qv為熱源密度；α為散熱系數(shù)；Tf為流體溫度。

通過(guò)對(duì)式(1)進(jìn)行變分，可將式可簡(jiǎn)化為如下表達(dá)式：

(2)

最后對(duì)式(2)作離散化處理，既得到如式(3)所示結(jié)果：

KT=F

(3)

式中：T為溫度矩陣，是在已知求解域內(nèi)所有節(jié)點(diǎn)溫度的矩陣；K,F分別為總體系數(shù)矩陣和總體右端矢量。

1.3 導(dǎo)熱系數(shù)的確定

在電機(jī)裝配的過(guò)程中，由于裝配工藝精度的限制，大型電機(jī)定子鐵心和機(jī)殼之間的裝配間隙不能足夠小，且有些間隙極不勻稱(chēng)，故在電機(jī)溫度場(chǎng)計(jì)算中氣隙是不允許忽略的。為方便建立模型，同時(shí)也為了不破壞計(jì)算精度，故本文將這一部分的導(dǎo)熱系數(shù)適當(dāng)放大。開(kāi)關(guān)磁阻電機(jī)主要材料導(dǎo)熱系數(shù)見(jiàn)表1。

表1 主要材料導(dǎo)熱系數(shù)

2 熱源的求取

電機(jī)作為一類(lèi)機(jī)電能量轉(zhuǎn)換的裝置，根據(jù)能量守恒定律可知，在其能量轉(zhuǎn)換的過(guò)程中產(chǎn)生損耗是無(wú)法消除的，而產(chǎn)生的損耗主要也將以熱能的形式表達(dá)出來(lái)，從而使電機(jī)定、轉(zhuǎn)子鐵心，繞組與機(jī)殼的溫度發(fā)生不同程度的變化。本文研究電機(jī)損耗主要為鐵心損耗和繞組銅耗。鐵心損耗主要是定子和轉(zhuǎn)子開(kāi)槽部位所引起的氣隙磁導(dǎo)諧波在對(duì)方鐵心中產(chǎn)生的能耗；以及當(dāng)電機(jī)帶適量負(fù)載運(yùn)行后所產(chǎn)生。繞組銅耗主要是當(dāng)繞組中有電流通過(guò)時(shí)在繞組中產(chǎn)生的能耗；以及工作電流產(chǎn)生的漏磁場(chǎng)和諧波磁場(chǎng)在繞組中產(chǎn)生的能耗。

2.1 發(fā)熱熱源的求取

SRM電機(jī)內(nèi)的發(fā)熱熱源主要來(lái)自于電機(jī)運(yùn)行時(shí)產(chǎn)生的損耗：(1)鐵心損耗，包括渦流損耗和磁滯損耗；(2)繞組損耗，主要是銅耗。

表2 電機(jī)主要熱源參數(shù)

由于開(kāi)關(guān)磁阻電機(jī)的雙凸極結(jié)構(gòu)特點(diǎn)，電壓和電流均為非正弦波，且磁通的波形不僅與電機(jī)實(shí)際工作轉(zhuǎn)速有關(guān)，還與電機(jī)在不同轉(zhuǎn)速時(shí)的時(shí)控方法密切相聯(lián)。即便當(dāng)電機(jī)工作在恒定轉(zhuǎn)速的狀態(tài)下，電機(jī)鐵心等部分磁通波形形狀及其變化規(guī)律也無(wú)法輕易分辨，故計(jì)算SRM電機(jī)損耗時(shí)需與傳統(tǒng)電機(jī)區(qū)別[1-2,9]。

2.2 鐵耗的計(jì)算

由于開(kāi)關(guān)磁阻電機(jī)的磁場(chǎng)所呈現(xiàn)的是非正弦規(guī)律，故可利用有限元的計(jì)算方法求得電機(jī)鐵心各部分的磁密波形。然后根據(jù)所得結(jié)果，對(duì)電機(jī)鐵耗進(jìn)行估算。

(4)

Ke和Kh分別為渦流修正因子和局部磁滯回線因子。這兩個(gè)參數(shù)的表達(dá)式為：

(5)

式中：B為磁密峰值；i為諧波次數(shù)；n為諧波次數(shù)最大值；k為介于0.6和0.7之間的常數(shù)；N為脈動(dòng)次數(shù)；ΔBl為磁密波形脈動(dòng)峰值。

2.3 繞組銅耗計(jì)算

定子銅耗PCu可根據(jù)下式計(jì)算：

(6)

式中：P為定子極數(shù)；I為繞組電流有效值；R(θ)為單相繞組電阻。

式中：Rk75為單相繞組75℃時(shí)的電阻值；Nr為轉(zhuǎn)子級(jí)數(shù)；θp為電流降至零時(shí)的轉(zhuǎn)子位移角。

3 溫度場(chǎng)仿真及結(jié)果分析

根據(jù)前文分析得出的導(dǎo)熱系數(shù)、主要熱源參數(shù)，在ANSYS中對(duì)已建立的開(kāi)關(guān)磁阻電機(jī)模型設(shè)置邊界條件、添加材料屬性、域、散熱面、熱源以及網(wǎng)格剖分，進(jìn)行三維溫度場(chǎng)仿真分析。

3.1 溫度場(chǎng)仿真

首先，在已建立的電機(jī)模型中，對(duì)電機(jī)各部分的材料參數(shù)進(jìn)行定義。

其次，得到有限元分析模型并進(jìn)行網(wǎng)格剖分。綜合考慮到計(jì)算精度及耗時(shí)問(wèn)題，在不影響電機(jī)計(jì)算性能的前提下，對(duì)機(jī)殼、繞組定子鐵心和轉(zhuǎn)子鐵心等區(qū)域采用不一樣的網(wǎng)格剖分。由于繞組模型的銅耗較大，定子模型的鐵耗較大，且溫升變化較明顯，故需對(duì)其網(wǎng)格進(jìn)行較細(xì)剖分。而其它溫升變化較小的模型部分可采用較低的精度進(jìn)行網(wǎng)格剖分。電機(jī)整體網(wǎng)格剖分示意如下圖所示。

圖2 開(kāi)關(guān)磁阻電機(jī)網(wǎng)格剖分示意圖

最后，加載邊界條件。溫度是電機(jī)溫度場(chǎng)分析時(shí)的自由度約束，邊界條件是整個(gè)模型的初始溫度。本文研究電機(jī)取初始溫度為25℃，即298.15 K。

3.2 結(jié)果分析

圖3給出了電機(jī)機(jī)殼、定子鐵心、定子繞組、冷卻水管和轉(zhuǎn)子鐵心在額定狀態(tài)下的溫度分布云圖。

(a)電機(jī)機(jī)殼溫度分布云圖(b)定子鐵心與水管溫度分布云圖

(c)定子繞組溫度分布云圖(d)冷卻水管溫度分布云圖

(e)轉(zhuǎn)子鐵心溫度分布云圖(f)電機(jī)整體溫度分布云圖

圖3 電機(jī)各部分溫度分布

由于電機(jī)損耗主要為繞組銅耗和鐵心損耗，此外，電機(jī)在高速旋轉(zhuǎn)時(shí)，轉(zhuǎn)子與氣隙摩擦?xí)r也會(huì)產(chǎn)生損耗。因此，電機(jī)的定子繞組與轉(zhuǎn)子鐵心溫度一般較高。

從圖3(c)和圖3(e)中的分析結(jié)果可知，電機(jī)定轉(zhuǎn)子鐵心的中部和定子繞組溫升較大。其中，轉(zhuǎn)子鐵心部分的最高溫度升至約95℃，而定子鐵心與繞組的所產(chǎn)生的最高溫度最終穩(wěn)定在約75℃。從圖3(d)中可知，電機(jī)冷卻水管的入水口與出水口溫度相差不大，但分布在定子鐵心中的水管水溫有著不同程度的升高。由于冷卻水管的存在，電機(jī)整體機(jī)殼溫度并未出現(xiàn)大幅上升，見(jiàn)圖3(a)和圖3(f)，其基本維持在30℃左右，屬于合理范圍之內(nèi)。

4 結(jié) 語(yǔ)

本文通過(guò)利用SolidWorks建模功能與ANSYS有限元分析軟件的熱分析功能對(duì)110 kW電動(dòng)汽車(chē)用開(kāi)關(guān)磁阻電機(jī)進(jìn)行三維溫度場(chǎng)仿真分析，得出如下結(jié)論：

(1)在自然水冷的情況下，電機(jī)穩(wěn)定運(yùn)行時(shí)，轉(zhuǎn)子溫度最高，定子繞組的溫度要略高于定子鐵心的溫度，但均符合電機(jī)的熱設(shè)計(jì)要求。

(2)此仿真溫度結(jié)果較為合理，所采取的熱源、導(dǎo)熱系數(shù)可靠準(zhǔn)確，為自然水冷式的開(kāi)關(guān)磁阻電機(jī)電動(dòng)機(jī)溫度場(chǎng)分析提供了參考。

[1] 詹瓊?cè)A.開(kāi)關(guān)磁阻電動(dòng)機(jī)[M].武漢：華中理工大學(xué)出版社，1992.

[2] 吳建華.開(kāi)關(guān)磁阻電機(jī)設(shè)計(jì)與應(yīng)用[M].北京：機(jī)械工業(yè)出版社，2000.

[3] 胡田，唐任遠(yuǎn)，李巖，等．永磁風(fēng)力發(fā)電機(jī)三維溫度場(chǎng)計(jì)算及分析[J].電工技術(shù)學(xué)報(bào)，2013，28(3):122-126.

[4] 陳鵬,朱憲然,魚(yú)振民.開(kāi)關(guān)磁阻電動(dòng)機(jī)流場(chǎng)分析及風(fēng)摩損耗計(jì)算[J].微特電機(jī).2008,36(2):20-23.

[5] 謝穎，辜承林.籠型感應(yīng)電動(dòng)機(jī)三維全域溫度場(chǎng)計(jì)算[J].中國(guó)電機(jī)工程學(xué)報(bào)，2012，32(36):96-101.

[6] 孔小光，王鳳翔，徐云龍．高速永磁電機(jī)鐵耗的分析和計(jì)算[J].電機(jī)與控制學(xué)報(bào)，2010，14(9):26-30．

[7] NAKAMURAK，F(xiàn)UJIMOTOH，F(xiàn)UJITSUNAM．Ironlossandmag-neticfieldsanalysisofpermanentmagnetmotorsbyimprovedfiniteelementmethodwithe&smodelelectronics.conference[J].IEEETransactionsonMagnetics,2001，37(5):3526-3529．

[8] 丁文，周會(huì)軍，魚(yú)振民.基于ANSYS的開(kāi)關(guān)磁阻電機(jī)溫度場(chǎng)分析[J].微電機(jī)，2005，38(5)：13-15.

[9] 王述成.開(kāi)關(guān)磁阻電動(dòng)機(jī)溫升研究[D].武漢：華中科技大學(xué),2004.

[10]RAULINV，RADUNA，HUSAINI.Modelingoflossesinswitchedreluctancemachines[J].IEEETransactionsonIndustryApplications,2004,40(6):1560-1569.

Analysis of the Temperature Field of Switched Reluctance Motor Using for Electric Vehicle

LIUXi-ping,ZENGJian,XUChen

(Jiangxi University of Science and Technology,Ganzhou 341000,China)

A reliable theoretical basis for the practical application for electric vehicle by assigning a numerical value to the working temperature field of a 110 kW water-cooling switched reluctance motor (12-slot 8-pole). The 3D model of the switched reluctance motor was established based on SolidWorks software. By using ANSYS Maxwell software, this paper calculated the copper loss and iron loss of the motor in rated condition to ensure the loss of heat; By using the finite element method, numerical analysis and calculation had been done for 3D motor Magnetic-thermal coupling simulation model, motor working in rated condition temperature field distribution had been deeply analyzed. Analysis results show that: under the rated conditions, internal temperature of the motor conforms to the requirements of the thermal design of the motor. The results provide a reliable theoretical basis for this class motor applied in electric vehicle.

Maxwell; switched reluctance motor; temperature field; coupling simulation

2015-12-23

國(guó)家自然科學(xué)基金項(xiàng)目(51267006)；江西省科技廳工業(yè)支撐項(xiàng)目(20151BBE50109)；江西省青年科學(xué)家培養(yǎng)對(duì)象項(xiàng)目(20151442040049)

TM352

A

1004-7018(2017)01-0023-03