薛敬業(yè),陳志輝,茆林儀

(南京航空航天大學(xué) 自動化學(xué)院,南京 211100)

0 引 言

無刷同步發(fā)電機(jī)在航空發(fā)電機(jī)領(lǐng)域得到了長期應(yīng)用,目前多數(shù)飛機(jī)的航空發(fā)電系統(tǒng)均采用了無刷同步發(fā)電機(jī)[1]。無刷同步發(fā)電機(jī)具有功率密度高的特點,在運行過程中會帶來較多的熱問題,如果超溫運行,會引起電機(jī)繞組絕緣材料的損壞,直接影響電機(jī)使用壽命。溫升是電機(jī)設(shè)計優(yōu)化需要關(guān)注的關(guān)鍵問題,而溫升最高點通常出現(xiàn)在電樞繞組處,所以研究無刷同步發(fā)電機(jī)中電樞繞組的溫升情況具有重要意義。

隨著計算機(jī)技術(shù)的發(fā)展,計算流體力學(xué)(以下簡稱CFD)方法成為了研究多物理場耦合問題的重要手段。文獻(xiàn)[2-7]采用CFD方法計算了高鐵異步牽引電機(jī)、發(fā)電機(jī)組、軸向磁通永磁電機(jī)、永磁發(fā)電機(jī)、雙定子開關(guān)磁阻電機(jī)和無刷直流電機(jī)的溫度場特性。文獻(xiàn)[8]以風(fēng)冷永磁電機(jī)為對象,研究了電機(jī)整體流體域的速度分布,并以此對電機(jī)溫升進(jìn)行估計。文獻(xiàn)[9]采用流固共軛傳熱數(shù)值計算方法對永磁電機(jī)的散熱性能進(jìn)行了仿真計算,并通過溫升測試驗證了仿真結(jié)果。文獻(xiàn)[10]提出了一種基于磁-熱-流體的多物理場模擬方法,對高速永磁同步電機(jī)進(jìn)行了熱分析。文獻(xiàn)[11]采用對稱性方法建立凸極風(fēng)力發(fā)電機(jī)模型,并進(jìn)行了溫度場分析,優(yōu)化了冷卻結(jié)構(gòu)。文獻(xiàn)[12]仿真計算了各種工況下永磁同步電機(jī)的轉(zhuǎn)子溫升,并采用溫度傳感器的測量實驗方案驗證了計算結(jié)果,轉(zhuǎn)子的仿真計算結(jié)果和溫升實驗結(jié)果存在13%的誤差。文獻(xiàn)[13]利用CFD方法,重點分析了轉(zhuǎn)速和轉(zhuǎn)子表面粗糙度對空氣摩擦損耗的影響規(guī)律。

由于現(xiàn)有計算方法存在較大的溫升誤差,為了更準(zhǔn)確得到電樞繞組的溫升,本文采用CFD方法,借助熱仿真工具研究30 kVA無刷同步發(fā)電機(jī)的電樞繞組溫升,得到了額定負(fù)載情況下電樞繞組溫度分布,并采用整體繞組分析法和繞組分段分析法這兩種磁熱耦合的迭代方法計算電樞繞組溫升。根據(jù)溫升測試結(jié)果,采用繞組分段方法可以提高電樞繞組溫升分析的精度,同時也為準(zhǔn)確研究電機(jī)的散熱情況提供了方法。

1 無刷同步發(fā)電機(jī)參數(shù)及模型

本文以采用自通風(fēng)冷卻方式的無刷同步發(fā)電機(jī)為研究對象,發(fā)電機(jī)絕緣等級為C級,主要參數(shù)如表1所示。

無刷同步發(fā)電機(jī)本身是一個強(qiáng)耦合的模型,所以考慮流固耦合的因素來建立等效的流固實體區(qū)域。由于無法考慮電樞繞組的股線絕緣、槽絕緣和層間絕緣,因而對定子槽內(nèi)及繞組端部的計算區(qū)域采用等效絕緣體的方法進(jìn)行處理,并對定子槽作如下假設(shè):

a) 電機(jī)浸漆狀態(tài)良好,浸漬漆填充均勻;

b) 銅線的絕緣漆均勻分布;

c) 忽略股線間絕緣漆膜存在所造成的溫差;

d) 認(rèn)為定子槽內(nèi)繞組的發(fā)熱情況相同,忽略繞組的集膚效應(yīng)。

根據(jù)以上假設(shè)將股線絕緣、層間絕緣以及槽絕緣等效為一個絕緣實體。其示意圖如圖1所示。

圖1 定子槽截面

圖1中,用槽面積減去裸導(dǎo)線面積得到等效絕緣體的總面積,再除以等效絕緣體的槽周長即可得到它的厚度d。電樞繞組采用一個整體緊密排列的銅線進(jìn)行等效。

電機(jī)熱穩(wěn)定的過程是周期對稱的,為提高網(wǎng)格剖分的質(zhì)量,減少計算資源的浪費,取發(fā)電機(jī)圓周的1/4區(qū)域作為CFD數(shù)值計算域,建立1/4的固體模型,如圖2所示。

圖2 1/4固體區(qū)域模型

由圖2可以看出,發(fā)電機(jī)主要由定子、轉(zhuǎn)子、電樞繞組、勵磁繞組、壓塊以及其他固定部件組成。

電機(jī)采用自通風(fēng)的冷卻方式,流體主要通過定子通風(fēng)道和轉(zhuǎn)子通風(fēng)道與發(fā)電機(jī)進(jìn)行對流換熱,從而對電機(jī)進(jìn)行散熱。發(fā)電機(jī)的1/4流體區(qū)域示意圖如圖3所示,主要包括定子通風(fēng)道流體、轉(zhuǎn)子通風(fēng)道流體、氣隙流體和外部流體4個部分。

圖3 1/4流體區(qū)域模型

2 熱仿真分析

2.1 數(shù)學(xué)模型

本文對無刷同步發(fā)電機(jī)額定負(fù)載狀態(tài)下的穩(wěn)態(tài)溫度場進(jìn)行數(shù)值計算,熱量傳遞主要以熱傳導(dǎo)和熱對流的形式進(jìn)行,不考慮熱輻射的影響。由于電機(jī)內(nèi)部流體雷諾數(shù)大,流動速度遠(yuǎn)小于聲速,所以將流動狀態(tài)假設(shè)為湍流流動,湍流模型選擇k-omega SST。

依據(jù)模型的簡化和求解模型,仿真過程中涉及到的k-omega SST控制方程如下:

(1)

(2)

式中:Uj是j方向的速度分量;xj是j方向上的坐標(biāo)分量;μ是流體黏度系數(shù);μt是渦流運動黏度系數(shù);k是湍流動能;τij是雷諾應(yīng)力;ω是湍流動能耗散率;F1是混合函數(shù);σk,β*,σω,γ,vT,β,σω2為模型系數(shù)。

以上控制方程是x,y平面下k方程和omega方程的對應(yīng)形式,x,z和y,z平面下的方程只要在該方程的基礎(chǔ)上更改成對應(yīng)的坐標(biāo)分量和速度分量。最終的空間變量由三個方程的分量合成得到。

2.2 網(wǎng)格剖分

本文采用多面體和六面體混合的網(wǎng)格剖分,邊界層網(wǎng)格采用多面體網(wǎng)格,其余部分采用六面體網(wǎng)格。與四面體網(wǎng)格相比,同樣條件下網(wǎng)格數(shù)量少,網(wǎng)格節(jié)點之間的過渡性好,計算精度高。網(wǎng)格剖分結(jié)果如圖4所示。

圖4 網(wǎng)格剖分結(jié)果

2.3 邊界條件

1) 固體

定子和電樞繞組采用靜止坐標(biāo)系,轉(zhuǎn)子、勵磁繞組及轉(zhuǎn)子上的部件采用固體運動,轉(zhuǎn)速為12 000 r/min。固體部分的材料屬性如表2所示。

表2 材料屬性

2) 流體

流體模型為理想氣體,采用旋轉(zhuǎn)參考坐標(biāo)系,轉(zhuǎn)速12 000 r/min,進(jìn)風(fēng)口為速度入口,通過實驗測得入口速度6.7 m/s,出風(fēng)口為壓力出口,出口表壓為1個大氣壓,流固接觸壁面為耦合壁面,環(huán)境溫度設(shè)為20 ℃。

2.4 熱源分布

電機(jī)各部分的損耗是計算溫度場的熱源,損耗計算完成之后,在熱仿真工具中轉(zhuǎn)換為單位體積的熱生成率,即可進(jìn)行熱計算。無刷同步發(fā)電機(jī)的總損耗主要由電樞銅損、勵磁銅損、阻尼銅損和定轉(zhuǎn)子鐵損組成。電樞銅損和勵磁銅損由式(5)計算,其余部分的損耗由有限元仿真計算,20 ℃下各部分損耗如表3所示。

表3 電機(jī)損耗

2.5 電樞繞組溫度結(jié)果

利用分離求解器SIMPLE對流固模型進(jìn)行計算,得到電樞繞組的溫度分布,如圖5所示。為了能夠?qū)﹄姌欣@組的溫度進(jìn)行實時測量,樣機(jī)加工過程中在電樞繞組上埋設(shè)了兩處測溫電阻,測溫點位置如圖5中所示。

圖5 電樞繞組溫度結(jié)果

發(fā)電機(jī)進(jìn)風(fēng)口側(cè)的溫度較低,測溫點1處的溫度58.49 ℃,出風(fēng)口側(cè)溫度較高,測溫點2處的溫度73.14 ℃。

3 磁熱耦合法

采用磁熱耦合法對電樞繞組溫升分析時,進(jìn)行如下假設(shè):

a) 電樞電流有效值不變,勵磁電流不變,忽略定轉(zhuǎn)子鐵損的變化;

b) 勵磁繞組與電樞繞組不直接換熱,勵磁繞組的溫升對電樞繞組溫升影響較小,忽略勵磁繞組銅損的變化;

c) 忽略空氣摩擦損耗。

根據(jù)以上假設(shè),只考慮電樞繞組電阻隨溫度的變化對電樞繞組溫升的影響。由于電樞繞組電阻隨溫度的升高而增大,電樞銅損增加,電樞繞組的溫升也會變化,所以初始條件下計算的電樞繞組溫升會與實際運行過程中有較大的差異。本文提出磁熱耦合方法,通過迭代對電樞繞組溫升進(jìn)行計算。磁熱耦合法流程圖如圖6所示。

圖6 磁熱耦合法流程圖

首先確定20 ℃下電樞繞組的電阻和銅損,進(jìn)行仿真,根據(jù)溫度結(jié)果,由積分求得電樞繞組平均溫度,再計算此溫度下電樞繞組的電阻和銅損,迭代求解,直到電樞銅損和溫升之間達(dá)到平衡,輸出溫升結(jié)果。

電樞繞組銅線電阻率和溫度的關(guān)系式:

ρT=ρT0[α(T-T0)+1]

(3)

式中:ρT為溫度T時的電阻率;ρT0為溫度T0時的電阻率。當(dāng)T0=20 ℃時,ρ20=1.72×10-8Ω·m;α=4.27×10-3K-1為溫度系數(shù)。

由式(3)可推出電樞繞組電阻計算公式:

R1=R2[α(T1-T2)+1]

(4)

式中:R1為溫度為T1時的電阻;R2為溫度為T2時的電阻。

電樞銅損的計算公式:

pCu=mIR2

(5)

式中:pCu是電樞銅損;m是發(fā)電機(jī)相數(shù);I為相電流有效值;R是相繞組電阻。

20 ℃下電樞繞組的相電阻為0.019 Ω,相電流有效值87 A,圖5是以此條件的仿真結(jié)果。隨著電樞繞組溫度的升高,電阻增大,電樞銅損增加,又會影響電樞繞組溫升。采用迭代計算的方法,由式(4)、式(5)求得平均溫度下電樞繞組的電阻和銅損,將每次的結(jié)果作為熱源條件再次計算,最終得到電樞繞組溫升。

本文重點對整體繞組分析法和繞組分段分析法進(jìn)行研究。

3.1 整體繞組分析法

根據(jù)20 ℃環(huán)境溫度下熱仿真數(shù)值計算結(jié)果,求得整體電樞繞組上的平均溫度,以此求出該溫度下電樞繞組的電阻值。

根據(jù)積分求得整體繞組上平均溫度65.60 ℃,電樞繞組相電阻變?yōu)?.023 6 Ω,電樞銅損變?yōu)?35.89 W。以此銅損作為熱源條件,迭代計算,得到電樞繞組的溫度分布如圖7所示。

圖7 整體繞組分析法第一次迭代溫度結(jié)果

測溫點1處的溫度64.90 ℃,測溫點2處的溫度80.37 ℃,積分求得電樞繞組平均溫度為72.15 ℃,相電阻為0.024 4 Ω,電樞銅損為544.05 W,進(jìn)行第二次迭代計算。

第二次迭代得到測溫點1處的溫度66.59 ℃,測溫點2處的溫度82.44 ℃。平均溫度73.68 ℃,此時電樞銅損和平均溫度72.15 ℃時相差不大,達(dá)到磁熱平衡的狀態(tài)。三次計算的溫升結(jié)果如表4所示。

表4 整體繞組分析法溫升結(jié)果

若采用熱仿真工具進(jìn)行迭代計算,最終可以得到穩(wěn)態(tài)時電樞繞組溫升,但這種方法的計算周期太長。采用擬合的方法,對得到的測溫點1、2處的溫升結(jié)果與計算次數(shù)的關(guān)系進(jìn)行擬合,擬合出的函數(shù)關(guān)系如下:

ΔT1=47.20-23.79·e-1.23n

(6)

ΔT2=63.27-35.38·e-1.25n

(7)

式中:ΔT1,ΔT2為測溫點1、2處的溫升;n為計算次數(shù)。擬合出的關(guān)系曲線如圖8所示。

圖8 溫升和計算次數(shù)關(guān)系曲線

當(dāng)計算次數(shù)為4時,由上式得測溫點1處的溫升47.03 ℃,測溫點2處的溫升63.03 ℃。將第二次迭代計算的結(jié)果繼續(xù)進(jìn)行仿真,得到測溫點1處的溫升46.86 ℃,測溫點2處的溫升62.77 ℃。仿真計算結(jié)果和函數(shù)計算結(jié)果基本一致,驗證了擬合函數(shù)關(guān)系的正確性。即采用上式得到穩(wěn)態(tài)時測溫點1處的溫升47.20 ℃,測溫點2處的溫升63.27 ℃。

3.2 繞組分段分析法

電樞繞組端部的溫度與平均溫度的差值較大,采用整體繞組分析法求得的平均溫度來代替繞組各部分的溫度,進(jìn)而求解的電樞銅損會有較大的偏差。為了更準(zhǔn)確計算電樞繞組溫升,采用繞組分段分析法進(jìn)行計算,將電樞繞組分兩段,對每段繞組分別求平均溫度,以此溫度來計算每段繞組的電阻和電樞銅損,再進(jìn)行迭代計算。

繞組分兩段的示意圖如圖9所示,在定子軸向中線的位置進(jìn)行分段。在截面位置只對電樞繞組進(jìn)行分段,其他部位并未進(jìn)行分段,對分段后的模型采用相同的方法進(jìn)行網(wǎng)格剖分,并在分段位置添加兩層邊界層網(wǎng)格。仿真計算時,將分段位置電樞繞組截面設(shè)置為耦合傳熱壁面,其余求解域采用的邊界條件和整體繞組分析法相同。

圖9 繞組分兩段示意圖

如圖9,根據(jù)a段和b段的長度比,求得20 ℃下a段的相電阻為0.009 67 Ω,b段的相電阻為0.009 33 Ω,a段電樞銅損219.58 W,b段電樞銅損211.86 W,其余損耗不變,采用相同的邊界條件,得到a段和b段繞組溫度分布如圖10所示。

圖10 繞組分兩段溫度結(jié)果

測溫點1處的溫度59.13 ℃,測溫點2處的溫度71.50 ℃。求得兩段平均溫度,即得到a段電樞銅損變?yōu)?66.13 W;b段電樞銅損變?yōu)?68.40 W。進(jìn)行再次計算,得到a、b段繞組的溫度分布情況。

迭代計算后,測溫點1處的溫度65.24 ℃,測溫點2處的溫度79.21 ℃。根據(jù)求得每段的平均溫度,a段電樞銅損變?yōu)?73.85 W,b段電樞銅損變?yōu)?78.39 W,此時電樞銅損和溫升之間還未完全達(dá)到平衡,需要再次計算。

第二次迭代計算后,測溫點1處的溫度67.57 ℃,測溫點2處的溫度80.52 ℃,和第一次迭代的計算結(jié)果相比,a段和b段的平均溫度基本相同,即此時發(fā)電機(jī)磁熱已經(jīng)平衡。

三次計算的溫升結(jié)果如表5所示。

表5 繞組分兩段的溫升結(jié)果

和整體繞組分析法相同,對測溫點1、2處的溫升結(jié)果和計算次數(shù)進(jìn)行擬合,擬合得到的函數(shù)關(guān)系式如下:

ΔT1=49.00-25.90·e-0.97n

(8)

ΔT2=60.79-54.67·e-1.77n

(9)

由上式可得,采用繞組分兩段的分析方法,得到穩(wěn)態(tài)時測溫點1處的溫升49.00 ℃,測溫點2處的溫升60.79 ℃。

若電樞繞組的分段越細(xì),每段溫度就越接近求得的平均溫度,即計算的電阻值和電樞銅損誤差也會越小。為了研究電樞繞組的分段數(shù)對溫升的影響,對電樞繞組進(jìn)行分四段處理,分四段示意圖如圖11所示。和繞組分兩段的方法相同,在對分段模型進(jìn)行網(wǎng)格剖分時,在分段位置分別添加兩層邊界層網(wǎng)格,求解時分段截面設(shè)置為耦合傳熱壁面,其余邊界條件保持不變。

圖11 繞組分四段示意圖

計算每段的電阻和電樞銅損,添加熱源,以相同的條件進(jìn)行計算,得到每段的平均溫度,進(jìn)行迭代。同樣經(jīng)過3次計算后,每段的平均溫度和電樞銅損基本不變化,已經(jīng)達(dá)到磁熱平衡的狀態(tài)。電樞繞組分四段的溫升結(jié)果如表6所示。

表6 繞組分四段的溫升結(jié)果

采用相同的方法進(jìn)行擬合,擬合的測溫點1、2處的溫升和計算次數(shù)的函數(shù)關(guān)系如下:

ΔT1=49.25-21.58·e-0.96n

(10)

ΔT2=60.48-60.88·e-1.84n

(11)

由上式可得,繞組分四段,穩(wěn)態(tài)時測溫點1處的溫升49.25 ℃,測溫點2處的溫升60.48 ℃。電樞繞組分兩段時,測溫點1處的溫升49.00 ℃,測溫點2處的溫升60.79 ℃。依據(jù)實驗結(jié)果,電樞繞組分四段得到的溫升,更接近實驗測得的溫升。因為分段越細(xì),求得的平均溫度來代替每段的溫度產(chǎn)生的誤差也就越小,計算的電樞銅損也更準(zhǔn)確,從而溫升也更接近實際運行中的溫升。

同時,電樞繞組分兩段的結(jié)果和分四段的結(jié)果比較相近,溫升誤差沒有大幅降低。而且增加電樞繞組的分段數(shù),模型的復(fù)雜度會提高,計算所消耗的時間和資源也會更多。因此在溫升誤差允許的范圍內(nèi),對電樞繞組分兩段處理,經(jīng)過3次計算,通過擬合的函數(shù)關(guān)系式即可得到穩(wěn)態(tài)時電樞繞組溫升。

4 實驗驗證

實驗平臺如圖12所示,主要由原動機(jī)、齒輪箱、無刷同步發(fā)電機(jī)和負(fù)載箱組成。由原動機(jī)經(jīng)過齒輪箱,拖動無刷同步發(fā)電機(jī)運行至額定轉(zhuǎn)速,完成額定阻性負(fù)載狀態(tài)下電樞繞組溫升測試實驗。

圖12 實驗平臺

為避免溫度過高對電機(jī)部件造成損壞,無刷同步發(fā)電機(jī)電樞繞組的進(jìn)風(fēng)口端部和出風(fēng)口端部分別埋設(shè)了鉑電阻(Pt100)作為傳感器測量電樞繞組溫度,埋設(shè)位置圖如圖13所示。溫度傳感器的引線通過位于機(jī)殼內(nèi)的鉆孔引出到電機(jī)外部,然后將引出線接至溫度顯示儀,可以實時顯示電樞繞組的溫度。

無刷同步發(fā)電機(jī)進(jìn)行溫升測試時,發(fā)電機(jī)外接額定阻性負(fù)載,由原動機(jī)拖動至額定轉(zhuǎn)速,運行至熱穩(wěn)定的狀態(tài),同時記錄下測溫點1和測溫點2處的溫度。溫升實驗總共進(jìn)行30 min,由于后期6分鐘內(nèi)溫升小于0.1 ℃,即認(rèn)為電機(jī)已經(jīng)達(dá)到該工況下的穩(wěn)態(tài)溫度。其中1組實驗結(jié)果如圖14所示。

圖13 測溫位置圖

圖14 實驗結(jié)果

進(jìn)行3組測量后取平均值,實驗數(shù)據(jù)如表7所示。

表7 溫升實驗結(jié)果

由表7看出,實測的測溫點1和2每次的溫度基本一致,實驗環(huán)境溫度20 ℃,實驗得到測溫點1處溫升51.33 ℃,測溫點2處溫升59.17 ℃。額定負(fù)載情況下電樞繞組溫升最高59.17 ℃,電機(jī)的絕緣等級為C級,有較強(qiáng)的過載能力,可以在惡劣的環(huán)境下應(yīng)用。

采用三種方法的溫升誤差分析結(jié)果如表8所示。

表8 溫升誤差分析

溫升誤差是通過擬合預(yù)測的溫升和實驗測量溫升之間的誤差。由表8可得,采用整體繞組分析法對電樞繞組溫升計算后,測溫點1處的溫升誤差8.05%,測溫點2處的溫升誤差6.93%,較為準(zhǔn)確度地計算了電樞繞組溫升。

采用繞組分段分析法計算后,對電樞繞組分兩段,得到測溫點1處的溫升誤差4.54%,測溫點2處的溫升誤差2.74%。比較整體繞組分析法,兩個測溫點的溫升誤差大幅減小,準(zhǔn)確度相對更高,仿真計算結(jié)果也更接近同步發(fā)電機(jī)運行過程中的溫升。對電樞繞組分四段進(jìn)行計算,測溫點1處的溫升誤差4.05%,測溫點2處的溫升誤差2.21%,經(jīng)過比較,溫升誤差結(jié)果和分兩段的結(jié)果基本一致。

5 結(jié) 語

無刷同步發(fā)電機(jī)在運行過程中溫升最高點出現(xiàn)在電樞繞組處,本文重點分析電樞繞組溫升情況。電樞繞組電阻和電樞銅損隨溫度升高而變化,很難準(zhǔn)確計算電樞繞組溫升。采用分段-磁熱耦合法預(yù)測電樞繞組溫升有重要意義。

1) 采用繞組分段分析法得到的電樞繞組的溫升結(jié)果比整體繞組分析法的溫升結(jié)果更加準(zhǔn)確,溫升誤差也更小,其中對電樞繞組分兩段分析,擬合得到的測溫點2處的溫升誤差減小了4.19%,為預(yù)計繞組溫升提供了方法。

2) 繞組分段越細(xì),計算的溫升結(jié)果更準(zhǔn)確。但為了節(jié)省仿真計算時間,此類電機(jī)采用分兩段的方法進(jìn)行溫升計算即可得到較高的溫升計算精度。