何聯(lián)格，胡書凡，屈 翔，嚴(yán) 新

(1.重慶理工大學(xué) 汽車零部件先進(jìn)制造技術(shù)教育部重點實驗室， 重慶 400054；2.北京信息科技大學(xué) 現(xiàn)代測控技術(shù)教育部重點實驗室， 北京 100192；3.重慶青山工業(yè)有限責(zé)任公司， 重慶 402761)

0 引言

在當(dāng)今能源枯竭的情況下，新能源汽車將逐漸替代傳統(tǒng)燃油車。由于新能源汽車中電機(jī)結(jié)構(gòu)十分緊湊，導(dǎo)致電機(jī)散熱變得比較困難[1]，而永磁同步電機(jī)在新能源汽車上應(yīng)用十分廣泛。電機(jī)在長久的工作過程中，過高溫升往往會引起電機(jī)一系列的問題，而這些問題通常影響著電機(jī)的使用壽命，甚至?xí)闺姍C(jī)遭到破壞[2]。有研究表明，35%左右的永磁同步電機(jī)失效是因為電機(jī)溫升過高而引起的[3]，因此，未來電機(jī)發(fā)展關(guān)鍵因素就是如何提高電機(jī)散熱效率，增強(qiáng)電機(jī)冷卻效果。

目前水冷和油冷是永磁同步電機(jī)最常用的2種冷卻方式，為了提升電機(jī)的冷卻效果，研究人員對這2種冷卻方式展開了大量的研究，如，沈啟平等[4]對水冷永磁同步電機(jī)進(jìn)行仿真分析，研究冷卻水道的溫度分布情況；文獻(xiàn)[5]設(shè)計出一種新型的軸向Z字型水道，不同于傳統(tǒng)的Z字型水道，其水道截面形狀為三角形，并通過實驗證明相對于傳統(tǒng)機(jī)殼冷卻方式，該冷卻方式冷卻效果更為顯著；文獻(xiàn)[6-8]提出了3種水道結(jié)構(gòu)的水冷方式(工字型、Z字型和螺旋型)，并對3種結(jié)構(gòu)進(jìn)行仿真計算分析，結(jié)果表明螺旋型水道的冷卻效果優(yōu)于其他2種，并對螺旋型水道做了深入的研究分析。

水冷方式雖然被廣泛應(yīng)用，但由于不能與電機(jī)內(nèi)部發(fā)熱部件直接接觸，導(dǎo)致某些部件溫度較高，電機(jī)內(nèi)部溫差較大。因此大量研究人員將目光投向油冷，如，杜愛民等[9]針對混合動力汽車用永磁同步電機(jī)結(jié)構(gòu)緊湊、散熱難的特點，采用電磁場與溫度場耦合的有限元分析方法對油冷冷卻條件下的永磁同步電機(jī)的三維溫度場進(jìn)行了計算；瑞典皇家理工學(xué)院的Nategh等[10]提出了在其定子上開油槽的油冷方式，增強(qiáng)了電機(jī)繞組端部的冷卻效果；Park等[11]提出一個將輪轂電機(jī)冷卻效果提25%的方法，就是在其原有的油冷結(jié)構(gòu)基礎(chǔ)上進(jìn)行優(yōu)化，使其冷卻通道從單通道變?yōu)槎嗤ǖ绹娪?。雖然油冷相對于水冷有著較強(qiáng)的優(yōu)勢，但其成本也相對較高，難以廣泛使用。

本文在前人研究的基礎(chǔ)上，針對繞組端部散熱難、局部溫升高這一問題，以一臺額定功率為38 kW、額定轉(zhuǎn)速為4 000 r/min的永磁同步電機(jī)為研究對象，在保持其螺旋型水道的前提下，設(shè)計出繞組端部噴油管，對其在4種典型工況下電機(jī)的溫升特性進(jìn)行了研究。

1 模型及參數(shù)確定

1.1 移動粒子半隱式法

MPS算法為移動粒子半隱式算法[12]。它是計算流體方法中最具權(quán)威的一個體系——基于拉格朗日體系，近年來Khayyer等[13]提出該算法的高階源項模型以及對其壓力模型的進(jìn)一步改進(jìn)。

MPS算法最大的優(yōu)勢是流體無網(wǎng)格，把流體定義為一顆顆細(xì)小的粒子，能夠更高精度地捕捉流體的流動軌跡，使流體仿真更接近實際，特別適用于減速器齒輪攪油和電機(jī)噴油冷卻。

1.1.1控制方程

連續(xù)方程[14]

(1)

Navier-Stokes 方程

(2)

1.1.2MPS算法

除壓力項的顯示計算

(3)

壓力的隱式計算

(4)

通過壓力梯度進(jìn)行速度和位置校正

(5)

式中：P為壓力；ρ為密度；k為時間步長，n為粒子數(shù)密度；*為階段物理量[14](完成顯式計算)；n0為初始狀態(tài)下的粒子數(shù)密度。

1.1.3粒子數(shù)密度與核函數(shù)

Shakibaeinia等[15]建立的核函數(shù)。

(6)

粒子數(shù)密度

(7)

1.1.4粒子作用模型

梯度矢量模型

(8)

式中：d為空間維數(shù)；φ為要計算的物理參數(shù)[14]；n0為初始狀態(tài)下的粒子數(shù)密度。

拉普拉斯模型

(9)

(10)

式中，λ在初始化起點處計算[14](該區(qū)域中所有粒子的平均值在整個瞬態(tài)計算中做常數(shù))。

散度模型：

(11)

1.1.5邊界條件

粒子數(shù)密度滿足以下條件的粒子被認(rèn)為是在自由表面上的粒子。

ni<βn0

(12)

式中：β為密度閾值[16](可示為一個小于1的常數(shù))；n0為初始狀態(tài)下的粒子數(shù)密度。

1.2 噴油冷卻結(jié)構(gòu)設(shè)計

油冷分為直接冷卻和間接冷卻2種方式。直接冷卻是冷卻油直接和發(fā)熱部件接觸并帶走其熱量的方式；間接冷卻是冷卻油流過冷卻油道帶走傳導(dǎo)到油道上的熱量。本文設(shè)計的油水復(fù)合冷卻方案是采用直接油冷的方式，在繞組端部設(shè)計出噴淋油管，通過淋油的方式將冷卻油噴淋到定子繞組端部上方，在重力作用下流動到繞組端部下方，帶走繞組端部的熱量，實現(xiàn)對繞組端部的冷卻。這種噴淋方式會因噴油管的設(shè)計位置導(dǎo)致冷卻不均勻，在臨近噴油管位置的地方流量大且油溫較低，冷卻效果好，繞組端部下段的位置流量小且冷卻油在流動的過程中被逐漸加熱，冷卻效果較差。

如圖1和圖2所示為水冷和油水復(fù)合冷卻下電機(jī)的有限元模型；表1為噴油管結(jié)構(gòu)參數(shù)。

圖1 水冷結(jié)構(gòu)有限元模型示意圖

圖2 油水復(fù)合冷卻有限元模型示意圖

表1 噴油管結(jié)構(gòu)參數(shù)

1.3 電機(jī)損耗計算

本文采用Maxwell軟件對永磁同步電機(jī)進(jìn)行電磁損耗仿真計算。

表2和表3為電機(jī)進(jìn)行電磁損耗分析所需基本參數(shù)以及材料的相對磁導(dǎo)率。硅鋼片的B-H曲線如圖3。

表2 樣機(jī)基本參數(shù)

表3 電磁材料的相對磁導(dǎo)率

圖3 硅鋼片的B-H曲線

電機(jī)有4種典型工況。額定工況1：額定功率(38 kW)，額定轉(zhuǎn)速(4 000 r/min)；額定工況2：額定功率(38 kW)，峰值轉(zhuǎn)速(12 000 r/min)；峰值工況1：峰值功率(95 kW)，轉(zhuǎn)折轉(zhuǎn)速(4 320 r/min)；峰值工況2：峰值功率(95 kW)，峰值轉(zhuǎn)速(12 000 r/min)。在這4種工況下對永磁同步電機(jī)的溫升特性進(jìn)行計算，其計算結(jié)果更能夠直接反映電機(jī)實際的溫升情況。

本文對電機(jī)電磁場仿真計算分析使用的軟件為Maxwell，表4—7所示分別為樣機(jī)在4種典型工況運行時的熱源分布計算值。

表4 額定工況1各構(gòu)件損耗 W

表5 額定工況2各構(gòu)件損耗 W

表6 峰值工況1各構(gòu)件損耗 W

表7 峰值工況2各構(gòu)件損耗 W

2 電機(jī)溫升特性研究

2.1 邊界條件

本文借助于shonDy和shonTA軟件對電機(jī)溫升特性進(jìn)行仿真計算。該模型中需要對電機(jī)冷卻液入口的2種冷卻物質(zhì)進(jìn)行參數(shù)設(shè)置。電機(jī)入口冷卻油流量為6 L/min(質(zhì)量流量為0.088 4 kg/s)，入口冷卻油溫度為80 ℃(353.15 K)；電機(jī)入口冷卻水流量為8 L/min(質(zhì)量流量為0.133 kg/s)，入口冷卻水溫度為70 ℃(343.15 K)。冷卻物質(zhì)的物理性質(zhì)如表8所示。

表8 冷卻物質(zhì)的物理性質(zhì)

2.2 溫升特性

額定工況1和額定工況2下定子繞組的溫度分布云圖如圖4。從圖4可以看出，對繞組端部加增噴油管，定子繞組整體溫度降低十分明顯，圖4(a)、(c)中可以看出水冷方案下繞組端部溫度最高，且高溫區(qū)與低溫區(qū)溫差較大(16、14 K)；從圖4(b)、(d)中可以看出油水復(fù)合冷卻方案下整體溫度降低，雖然最高溫度仍在繞組端部，但已從整個繞組端部轉(zhuǎn)移到繞組端部的下端，且高低溫區(qū)溫差較小(6、4 K)；電機(jī)最高溫度如表9。

圖4 定子繞組額定工況溫度場分布云圖

表9 額定工況電機(jī)最高溫度 ℃

圖5為永磁同步電機(jī)油水符合冷卻方案下不同油水復(fù)合比例下的溫度云圖，是在保持水的流量不變的情況下，采用不同冷卻油的流量；表10為額定工況下不同油水復(fù)合比例下電機(jī)最高溫度。從圖表中可以看出在雖然增大冷卻油流量會使繞組最高溫度進(jìn)一步降低，但是隨著流量的增加其溫度降低的趨勢會越來越小。

圖5 油水混合不同比例下溫度云圖

表10 油水不同比例下電機(jī)最高溫度 ℃

圖6為水冷永磁同步電機(jī)溫升試驗臺架場景，本次溫升試驗是通過在電機(jī)繞組端部選取4處測溫點埋置檢溫計的方法(見表11)，測出電機(jī)實際的溫升情況，以便于與之后的仿真數(shù)據(jù)進(jìn)行對比。

圖6 溫升臺架試驗場景

表11 4處采樣點溫度 ℃

由圖7(a)可看出，額定工況1下，水冷電機(jī)溫升試驗得出繞組最高溫度穩(wěn)定在141 ℃；水冷仿真繞組最后溫度穩(wěn)定在133.85 ℃；油水復(fù)合冷卻下繞組最高溫度穩(wěn)定在106.85 ℃，相比其水冷卻繞組最高溫度降低了27 ℃。由圖7(b)可看出，額定工況2下，水冷電機(jī)溫升試驗得出繞組最高溫度穩(wěn)定在149 ℃；水冷仿真繞組最高溫度穩(wěn)定在142.85 ℃；油水復(fù)合冷卻下繞組最高溫度穩(wěn)定在115.85 ℃，相比其水冷卻繞組最高溫度降低了27 ℃。在2個額定工況運行下，電機(jī)定子鐵芯和繞組穩(wěn)定狀態(tài)下的溫度分布云圖復(fù)合油水復(fù)合冷卻方案的設(shè)計要求。

圖7 電機(jī)額定工況溫升曲線

圖8為峰值工況下電機(jī)溫升曲線。由圖8(a)可以看出，電機(jī)在水冷情況下達(dá)到150 ℃需要44 s，在油水復(fù)合冷卻情況下達(dá)到150 ℃需要380 s；由圖8(b)可以看出電機(jī)在水冷情況下達(dá)到150 ℃需要60 s，在油水復(fù)合冷卻情況下達(dá)到150 ℃需要520 s；試驗得出的標(biāo)準(zhǔn)是電機(jī)達(dá)到150 ℃所需的時間應(yīng)超過30 s，故其性能復(fù)合要求。

圖8 電機(jī)峰值工況溫升曲線

3 結(jié)論

在電動汽車永磁同步電機(jī)水冷的基礎(chǔ)上設(shè)計出油水復(fù)合冷卻方案，對永磁同步電機(jī)溫升特性研究，并對2種冷卻方案下的永磁同步電機(jī)溫度場進(jìn)行對比和分析，得到如下結(jié)論。

1) 2種額定工況下，油水復(fù)合冷卻方案較水冷方案下電機(jī)最高溫度降低27 ℃，冷卻效果提升20%、19%。

2) 2種峰值工況下，電機(jī)在油水復(fù)合冷卻方案下達(dá)到150 ℃所需時間為380、520 s,其時間遠(yuǎn)大于水冷方案下電機(jī)達(dá)到150 ℃所需時間。

3) 定子鐵芯和繞組在采用油水復(fù)合冷卻方案下溫度低于水冷方案，電機(jī)繞組最高溫度大幅度降低，其高低溫區(qū)溫差較小，起到良好的溫度均衡作用。