崔 友，席榮盛，謝培利，何茂興

0 引言

近年來隨著能源枯竭及環(huán)境污染問題日益突出，各國大力發(fā)展新能源電動汽車。但電動車發(fā)展的瓶頸是電池續(xù)航能力不足，為此増程式電動汽車應(yīng)運(yùn)而生[1]。增程式電動車可以有效降低蓄電池的成本和重量，實(shí)現(xiàn)燃油汽車向純電動汽車的平穩(wěn)過渡。在美國市場，作為典型增程式電動車型雪佛蘭Volt在2016年全美銷量排行第二，目前別克推出的VELITE 5增程式混動汽車也是基于Volt平臺開發(fā)的。國內(nèi)部分車企也一直致力于增程式電動車的開發(fā)，奇瑞、吉利，廣汽傳祺都在進(jìn)行研發(fā)工作，收購了美國Karma公司萬向汽車也有計(jì)劃投資增程式電動車。由于受使用空間限制，要求增程式電動車交流永磁同步電機(jī)功率密度相比驅(qū)動電機(jī)高，但同時(shí)損耗密度也增加，使電機(jī)工作時(shí)各部件的溫度較高，嚴(yán)重時(shí)會影響電機(jī)運(yùn)行。因而要求工程師全面了解電機(jī)的整體溫度分布。目前傳統(tǒng)的電機(jī)溫度場分析基本都是從電機(jī)整體出發(fā)考慮，計(jì)算結(jié)果大多為平均溫度值，而暫態(tài)有限元計(jì)算法可以得到電機(jī)每一點(diǎn)的溫度變化值。本文充分研究增程式電動車永磁同步電機(jī)定轉(zhuǎn)子鐵芯、定子繞組及磁鋼渦流損耗分布，建立了電機(jī)三維暫態(tài)溫度場數(shù)值分析模型，利用有限元軟件進(jìn)行計(jì)算。通過對樣機(jī)進(jìn)行溫升測試，與計(jì)算結(jié)果相符，為增程式電動車交流永磁電機(jī)的實(shí)際應(yīng)用提供了參考依據(jù)。

1 增程式電動車交流永磁同步電機(jī)三維暫態(tài)溫度場計(jì)算模型

傳統(tǒng)的熱場分析包括熱傳導(dǎo)、熱對流及熱輻射。對于電機(jī)來講，溫度場計(jì)算以熱傳導(dǎo)方式為主。近代研究電機(jī)的溫度計(jì)算以穩(wěn)態(tài)運(yùn)行為主。穩(wěn)定狀態(tài)下，由傅里葉方程推導(dǎo)的增程式電動車交流永磁電機(jī)三維熱傳導(dǎo)方程為[2]：

式中，T 為電機(jī)的溫度（℃）；Kx、Ky、Kz分別為電機(jī)各介質(zhì) x、y、z方向的導(dǎo)熱系數(shù)（w/mm℃）；q為熱源密度 （w/mm3）；c為熱容（w.s/kg℃）；γ 為密度（kg/mm3）；τ為時(shí)間（s）.Te為周圍介質(zhì)的溫度（時(shí)間的函數(shù)）（℃）；α為散熱系數(shù)（w/mm2℃）；K為熱傳導(dǎo)系數(shù)。

增程式電動車交流永磁同步電機(jī)穩(wěn)態(tài)運(yùn)行時(shí)，由于電機(jī)在結(jié)構(gòu)上軸向?qū)ΨQ，認(rèn)為電機(jī)軸向中間部分為熱平衡狀態(tài)，由于電機(jī)周向?qū)ΨQ結(jié)構(gòu)，電機(jī)溫度場也周向?qū)ΨQ分布。由于端部繞組產(chǎn)生的損耗相比電機(jī)其他部件可忽略，因此在計(jì)算中為簡化模型，提高計(jì)算效率，可忽略繞組端部的影響?；谥芟?qū)ΨQ結(jié)構(gòu)，取轉(zhuǎn)子一極的范圍，軸向取半個(gè)鐵芯長度作為計(jì)算邊界條件[3]，溫度場計(jì)算模型如圖1所示。

圖1 計(jì)算模型

為簡化計(jì)算，對模型求解域做如下基本假設(shè)：

（1）不考慮繞組集膚效應(yīng)；

（2）定轉(zhuǎn)子鐵芯端面及機(jī)殼表面的散熱系數(shù)取平均值；

（3）假設(shè)某一時(shí)刻下定子鐵芯損耗恒定；

（4）不考慮極弧系數(shù)對溫度分布的影響。

2 增程式電動車永磁同步電機(jī)損耗及散熱系數(shù)的確定

2.1 銅耗的確定

在計(jì)算永磁同步電機(jī)瞬態(tài)溫度場中，電樞繞組的相電阻隨溫度變化而變化，銅損可采用下式計(jì)算：

式中，Pcu為銅損；Iφ為電機(jī)相電流；b為繞組線規(guī)直徑；n為線圈總導(dǎo)體數(shù)；bs為定子槽寬；Rs為定子繞組相電阻；fn為電機(jī)運(yùn)行頻率。

2.2 鐵芯損耗的確定

鐵芯損耗包括定子鐵芯損耗及轉(zhuǎn)子鐵芯損耗，轉(zhuǎn)子磁鋼渦流損耗。電機(jī)運(yùn)行過程中，由于軸向?qū)ΨQ，可認(rèn)為磁場沿鐵心軸向無變化，可采用二維電磁場有限元法取得電機(jī)不同時(shí)刻周向磁場分布，結(jié)合定轉(zhuǎn)子鐵芯材料特性，計(jì)算出一極下鐵芯不同時(shí)刻的鐵耗分布。本文由JMAG有限元電磁場仿真軟件計(jì)算[4]。

電機(jī)負(fù)載工作時(shí)，由于定子繞組電樞反應(yīng)，定子磁動勢及齒諧波磁動勢與轉(zhuǎn)子運(yùn)行頻率不一致，會在轉(zhuǎn)子鐵芯及磁鋼中產(chǎn)生渦流，同時(shí)磁鋼會有附加渦流損耗，高溫下磁鋼會有退磁風(fēng)險(xiǎn)。對于磁鋼渦流損耗，在計(jì)算時(shí)，給定磁鋼材料電導(dǎo)率并賦予零電流，計(jì)算一個(gè)電周期下磁鋼渦流損耗[5]。但實(shí)際電機(jī)磁鋼軸向分塊，可減小磁鋼渦流損耗。在進(jìn)行磁鋼渦流有限元分析時(shí)，剖分精度對計(jì)算結(jié)果影響很大，為提高計(jì)算精度，應(yīng)對每個(gè)剖分單元進(jìn)行渦流損耗瞬時(shí)值計(jì)算：

式中：σe為磁鋼剖分單元的電導(dǎo)率；Je（t）為單元中的電流密度瞬時(shí)值；Ve為剖分單元包含的計(jì)算區(qū)域。

則轉(zhuǎn)子旋轉(zhuǎn)一個(gè)電周期時(shí)產(chǎn)生渦流損耗平均值為[7]：

式中：T0為電機(jī)旋轉(zhuǎn)一個(gè)電周期的時(shí)間；E為磁鋼剖分單元總體數(shù)。

建立永磁電機(jī)模型之后，對永磁體設(shè)置電導(dǎo)率的方法來求解渦流損耗。

2.3 電機(jī)導(dǎo)熱系數(shù)的確定

由于增程式永磁同步電機(jī)溫度場以傳熱為主，因此導(dǎo)熱系數(shù)的確定對溫度場分析的準(zhǔn)確程度至關(guān)重要。對于結(jié)構(gòu)復(fù)雜永磁同步電機(jī)，各部件導(dǎo)熱系數(shù)差別較大。其導(dǎo)熱系數(shù)參考表1.導(dǎo)熱系數(shù)是隨材料的溫度變化而變化的。

表1 電機(jī)各部件材料導(dǎo)熱系數(shù)

2.4 電機(jī)表面主要散熱系數(shù)的確定

本文設(shè)計(jì)的增程式永磁同步電機(jī)采用水冷方式，需確定各部件表面的散熱系數(shù)。圖2為定子鐵芯與轉(zhuǎn)子鐵芯氣隙表面散熱系數(shù)，可采用下式計(jì)算[6]：

圖2 額定點(diǎn)定轉(zhuǎn)子磁密分布

式中，Vn為機(jī)殼內(nèi)冷卻介質(zhì)的流速

3 算例與結(jié)果分析

3.1 三維暫態(tài)溫度場計(jì)算

根據(jù)對熱源及散熱條件的分析，本文對一款增程式電動車永磁同步電機(jī)進(jìn)行了計(jì)算。電機(jī)額定及結(jié)構(gòu)數(shù)據(jù)如下：

PN＝ 68 kW；Un=380 V；nN=3200 rpm；fn=320 Hz；定子鐵芯外徑D0=290 mm；定子鐵芯內(nèi)徑D1=198.4 mm；鐵芯軸向長L=90 mm；極數(shù)P=12；槽數(shù)Z=54，氣隙g=0.8 mm，定子繞組6Y型連接，轉(zhuǎn)子采用嵌入式一字型磁鋼，定子機(jī)殼水循環(huán)冷卻，見圖3.

圖3 某一時(shí)刻下電機(jī)內(nèi)部溫度場分布

在進(jìn)行瞬態(tài)溫度場計(jì)算時(shí)，需將各個(gè)時(shí)刻損耗轉(zhuǎn)換為熱源，并賦予對應(yīng)時(shí)刻下電機(jī)相應(yīng)各部件。從初始時(shí)刻開始分析，記錄時(shí)刻初始溫度為，計(jì)算此時(shí)刻電磁場，此時(shí)電機(jī)未啟動。將轉(zhuǎn)子旋轉(zhuǎn)到下一時(shí)刻，計(jì)算該時(shí)刻下電磁場得到電機(jī)各部件損耗分布，轉(zhuǎn)換為熱源后進(jìn)行此時(shí)刻溫度場分布，確定時(shí)刻的溫度分布。再以時(shí)刻溫度為起點(diǎn)，按以上計(jì)算方法繼續(xù)計(jì)算，得到時(shí)刻的溫度分布，不斷迭代計(jì)算步驟，就可得到從起動到電機(jī)穩(wěn)定運(yùn)行過程中的暫態(tài)溫度場分布情況。見圖4.

圖4 實(shí)驗(yàn)計(jì)算溫升對比曲線

3.2 樣機(jī)溫升測試

根據(jù)以上計(jì)算結(jié)果，對試驗(yàn)樣機(jī)進(jìn)行溫升測量實(shí)驗(yàn)。用驅(qū)動器驅(qū)動電機(jī)至額定轉(zhuǎn)速，并施加額定負(fù)載，記錄不同時(shí)刻定子繞組溫度，可得到如圖4的溫升曲線。實(shí)驗(yàn)結(jié)果表明電機(jī)溫度運(yùn)行120 min以后定子繞組溫度變化趨于穩(wěn)態(tài)，溫度溫度為142℃左右。由圖4的計(jì)算與試驗(yàn)溫升對比曲線可知，本文中采用的增程式永磁同步電機(jī)溫度場計(jì)算模型能夠較合理的分析計(jì)算該電機(jī)暫態(tài)溫升過程。

4 結(jié)論

（1）本文根據(jù)傳熱學(xué)建立了三維增程式電動車永磁同步電機(jī)溫度場計(jì)算模型，并與試驗(yàn)樣機(jī)實(shí)測結(jié)果進(jìn)行對比，結(jié)果表明該計(jì)算模型能夠滿足工程上對增程式電機(jī)暫態(tài)溫度場仿真計(jì)算的需求。

（2）電機(jī)轉(zhuǎn)子采用永磁體，無勵(lì)磁損耗，只產(chǎn)生可忽略的渦流損耗，相對于定子繞組溫升來說轉(zhuǎn)子磁極溫升較低，對電機(jī)永磁體性能及退磁的影響很小。由于電機(jī)氣隙的存在，定轉(zhuǎn)子間的熱交換相對獨(dú)立，因此后續(xù)可以將定轉(zhuǎn)子鐵芯分開計(jì)算溫度場。

[1]許家群，于慎波，徐衍亮，等.電梯曳引驅(qū)動系統(tǒng)的現(xiàn)狀及發(fā)展前景[J].微特電機(jī)，2002（3）：5-7．

[2]葉安麗，周渡海，李惠升.對永磁同步電動機(jī)在電梯中應(yīng)用的探討與分析[J].電工技術(shù)雜志，2003（6）：61-63．

[3]李德基，白亞民，曹國宣.發(fā)電機(jī)暫態(tài)三維溫度場的計(jì)算[J].中國電機(jī)工程學(xué)報(bào)，1989，9（5）：56-63．

[4]胡錫棟，高志剛.發(fā)電電動機(jī)轉(zhuǎn)子磁極暫態(tài)溫度場分析[J].杭州電子科技大學(xué)學(xué)報(bào)，2012，32（3）：25-28.

[5]周 封，熊 斌，李偉力，等.大型電機(jī)定子三維流體場計(jì)算及其對溫度場的影響[J].中國電機(jī)工程學(xué)報(bào)，2005，25（24）：128-132．

[6]溫志偉，顧國彪.實(shí)心磁極同步電動機(jī)轉(zhuǎn)子溫度場計(jì)算[J].大電機(jī)技術(shù)，2005（2）：1-5.

[7]楊 通，周理兵，蔣 煒.實(shí)心籠型異步電機(jī)轉(zhuǎn)子渦流分析與計(jì)算[J].電機(jī)與控制應(yīng)用，2006，33（4）：3-6．