聶華偉，史忠震，許貞俊

(1.貴州交通職業(yè)技術(shù)學(xué)院,貴陽 550008；2.貴州省機(jī)械工業(yè)學(xué)校,貴陽 550001)

?

電機(jī)定子繞組熱模型等效方法的有限元分析研究

聶華偉1，史忠震1，許貞俊2

(1.貴州交通職業(yè)技術(shù)學(xué)院,貴陽 550008；2.貴州省機(jī)械工業(yè)學(xué)校,貴陽 550001)

通過對實(shí)際電機(jī)定子繞組結(jié)構(gòu)的分析研究，得出了繞組熱模型的建立方法，并針對電機(jī)熱場研究建立了三種簡化模型，運(yùn)用ANSYS熱仿真系統(tǒng)分析和實(shí)驗(yàn)研究對比得到最為合理的繞組熱模型。得出的簡化模型方法對電機(jī)熱場研究精確度提高有十分重要意義，同時(shí)對電機(jī)相關(guān)參數(shù)性能有限元分析等效模型的建立有一定的參考意義。

定子繞組；等效模型；溫升研究；有限元

0 引 言

電機(jī)作為動(dòng)力裝置，在各種機(jī)械設(shè)備上得到了廣泛應(yīng)用，也是其不可缺少的一部分，它的性能直接影響著設(shè)備的正常運(yùn)行，而熱是直接影響電機(jī)性能的參數(shù)之一，如果溫度超過本身所承受范圍將使電機(jī)壽命以及使用性能受到影響，因此在電機(jī)設(shè)計(jì)過程中應(yīng)對溫度重點(diǎn)考慮。

在電機(jī)熱場計(jì)算過程中電機(jī)模型的等效、熱源以及各部件溫度參數(shù)的計(jì)算精度直接影響有限元法的溫度場計(jì)算精度，其中熱源以及各部件溫度參數(shù)的計(jì)算國內(nèi)外學(xué)者已經(jīng)做了大量工作，并取得了豐碩的成果，但對于電機(jī)有限元簡化模型建立研究不多[12]。在電機(jī)模型等效過程中，繞組模型的合理性建立是最困難的，若根據(jù)電機(jī)繞組實(shí)際模型去建立，是根本不可能達(dá)到的，即便可以完整的建立，在分析過程中也會(huì)由于網(wǎng)格的劃分問題，不能得到電機(jī)正確的溫度場分布，因此繞組模型的等效合理性已成為有限元法研究電機(jī)溫升的關(guān)鍵。

本文通過運(yùn)用ANSYS仿真軟件分別對三種簡化模型分析，根據(jù)溫升實(shí)驗(yàn)研究對比得出3種等效模型中最為合理的繞組熱模型。得出的簡化模型方法對電機(jī)熱場研究精確度提高有十分重要意義，同時(shí)對電機(jī)相關(guān)參數(shù)性能有限元分析等效模型的建立有一定的參考意義。

1 電機(jī)繞組的等效模型

該實(shí)驗(yàn)研究電機(jī)定子的具體參數(shù)為：定子槽數(shù)為12，每槽導(dǎo)體數(shù)44個(gè)，6線并繞，漆包線線徑為0.5 mm，銅線線徑為0.45 mm，槽滿率為40%，實(shí)際形狀如圖1所示。

圖1 定子實(shí)際結(jié)構(gòu)圖

建立模型前，為便于更好的進(jìn)行有限元分析計(jì)算，對此進(jìn)行相應(yīng)假設(shè)：

(1)每個(gè)定子槽嵌入的導(dǎo)線均勻排列，且均溫；

(2)定子疊片內(nèi)的絕緣材料和繞組銅線絕緣材料為完全均質(zhì)；

(3)定子疊片內(nèi)部完全填充絕緣漆并且均勻[10]；

(4)定子疊片層層無間隙是一個(gè)整體。

其中建立的繞組模型應(yīng)與實(shí)際模型相接近，其中包括繞組模型面積、重量、槽滿率。

通過以上假設(shè)以及必須遵循的原則建立了以下等效模型。

模型一：1) 定子疊片簡化成整體；2) 不考慮槽內(nèi)絕緣層其中的浸漆層；3) 將定子槽內(nèi)的繞組模型建成一個(gè)整體，外部為混合絕緣層；4) 外部空間由空氣包裹。具體模型結(jié)構(gòu)圖如下所示。

圖2 定子模型簡化方案一

模型二：1) 定子疊片簡化成整體；2) 考慮槽內(nèi)絕緣層其中的浸漆層的存在；3) 將定子槽內(nèi)的繞組模型建成一個(gè)整體，外部為混合絕緣層；4) 外部空間由空氣包裹。具體模型結(jié)構(gòu)圖如下所示。

圖3 定子模型簡化方案二

模型三：根據(jù)定子的實(shí)際結(jié)構(gòu)特點(diǎn)建立定子模型。1) 不忽略任何絕緣介質(zhì)的存在；2) 槽內(nèi)繞組等效成層狀，每兩個(gè)繞組層之間為絕緣層。具體模型結(jié)構(gòu)如圖4所示。

圖4 定子模型簡化方案三

2 電機(jī)定子的熱分析

運(yùn)用仿真軟件ANSYS分別對三種模型進(jìn)行熱場分析研究。

2.1 熱參數(shù)的計(jì)算

分析可知熱參數(shù)包括熱源、散熱系數(shù)以及導(dǎo)熱系數(shù)[10]。

(1)產(chǎn)熱源

該有限元研究電機(jī)運(yùn)行狀態(tài)為靜止情況，因此不存在機(jī)械損耗，鐵心損耗相對于繞組銅耗來說忽略不計(jì)。根據(jù)實(shí)驗(yàn)所得繞組在起始通電流的瞬時(shí)銅耗為21.5 W，當(dāng)通電一定時(shí)間繞組溫度達(dá)到穩(wěn)定狀態(tài)，這時(shí)銅耗為26.8 W。由于在溫度場分析中熱源加載以生熱率進(jìn)行加載，因此計(jì)算出起始瞬時(shí)和穩(wěn)定狀態(tài)時(shí)生熱率分別為1.13×106W/m3和1.41×106W/m3。

(2)熱傳導(dǎo)系數(shù)

該電機(jī)定子零部件的有關(guān)熱傳導(dǎo)系數(shù)如下表所示[10]。

表1 相關(guān)熱傳導(dǎo)系數(shù)

根據(jù)經(jīng)驗(yàn)公式，多種混合體材料的熱傳導(dǎo)系數(shù)[1]：

(1)

式中：λeq為等效導(dǎo)熱系數(shù)；δi(i=1,2,3,…,n)為每種材料的厚度；λi(i=1,2,3,…,n)為每種材料的熱傳導(dǎo)系數(shù)。

通過以上經(jīng)驗(yàn)公式得出3種簡化模型槽內(nèi)絕緣體的等效導(dǎo)熱系數(shù)分別為：模型一λeq=0.29 W/(m·K)；模型二和模型三λeq=0.17 W/(m·K)。

(3)散熱系數(shù)

由于該電機(jī)定子溫度場研究是在室內(nèi)，相對電機(jī)來說可視為無限空間，因此散熱系數(shù)根據(jù)以下經(jīng)驗(yàn)公式計(jì)算[2]：

(2)

當(dāng)104

(3)

當(dāng)109

(4)

式中：Nu，Gr，θ，L，Pr分別為努塞爾系數(shù)、格拉曉夫系數(shù)、電機(jī)外表面與環(huán)境溫差、電機(jī)與空氣接觸面長度和普朗特系數(shù)；β，μ，k，cp分別為空氣的相關(guān)參數(shù)：體積膨脹系數(shù)、動(dòng)力粘度、熱傳導(dǎo)系數(shù)和比熱容。

由于所測電機(jī)定子周圍空氣溫度15 ℃，根據(jù)式(2)、式(3)、式(4)得出與周圍空氣接觸的散熱系數(shù)h=7.1 W/(m2·K)。

2.2 計(jì)算結(jié)果云圖

將以上求得相關(guān)參數(shù)加載分析，分析結(jié)果云圖如圖5所示。

(a)(b)

(c)

3 實(shí)驗(yàn)研究

通過實(shí)驗(yàn)研究分析來找出以上3種簡化模型最合理的簡化方法。

在定子繞線時(shí)將3個(gè)同型號的溫度傳感器繞于內(nèi)部，但繞放的位置不同，分別放在繞組中心、繞組中和繞組最外層，具體放置位置如圖6所示。

圖6 傳感器放置位置示意圖

3.1 搭建實(shí)驗(yàn)平臺(tái)

由于該實(shí)驗(yàn)研究針對溫度場，因此實(shí)驗(yàn)所用裝置主要包括電源(直流穩(wěn)壓)、電壓和電流傳感器、熱巡檢儀以及熱電偶傳感器等組成，分布示意圖如圖7所示。

圖7 分布圖

研究實(shí)驗(yàn)步驟:

(1)根據(jù)定子溫度場實(shí)驗(yàn)研究示意圖將實(shí)驗(yàn)裝置連接。

(2)通入穩(wěn)定直流電流6.26 A，記錄起始電壓值，計(jì)算出剛通電流瞬時(shí)繞組阻值。

(3)當(dāng)溫度平穩(wěn)時(shí)即達(dá)到熱穩(wěn)定狀態(tài)，這時(shí)記錄電壓值并計(jì)算出穩(wěn)定狀態(tài)下繞組的阻值。

(4)記錄穩(wěn)態(tài)下溫度值。

連接出實(shí)驗(yàn)研究裝置如下。

圖8 實(shí)驗(yàn)研究裝置

3.2 實(shí)驗(yàn)研究數(shù)據(jù)分析

將儀器中所測熱數(shù)據(jù)導(dǎo)出，并作出統(tǒng)計(jì)曲線如下圖所示。

圖9 定子繞組溫度測試曲線圖

根據(jù)圖9分析可知，定子繞組在通入直流電流大約6 h后達(dá)到穩(wěn)態(tài)，此時(shí)TMAX=124.6，在繞組中心處，其它兩處溫度依次為T1=122.7和T2=121.3，符合實(shí)際溫度分布情況，計(jì)算得平均溫度122.9℃。

3.3 對比分析

仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)對比如表2所示。

表2 對比表

通過上表得出結(jié)論：(1)模型一與實(shí)驗(yàn)對比誤差在 -20%以上，遠(yuǎn)高于所要求的±10%工程誤差范圍，誤差較大。分析可知該誤差是由于等效模型將定子槽內(nèi)的浸漆層忽略，而浸漆層材料的導(dǎo)熱系數(shù)很小，使得定子繞組產(chǎn)生的熱量能夠相對較好的傳遞出去，從而導(dǎo)致仿真出的繞組溫度相對較低。

(2)模型二的繞組各溫度誤差均接近工程要求10%誤差，誤差相模型一較小，但仍不能滿足工程精確度要求高的場合。

(3)模型三的繞組最高溫度、最低溫度和平均溫度與實(shí)際實(shí)驗(yàn)數(shù)據(jù)相接近，誤差較小處于工程要求10%范圍內(nèi)，精度較高。

前兩個(gè)簡化模型相比，模型二的誤差較小并且誤差接近工程要求，因此在精度要求范圍低的情況下，這種簡化方法可以采用。對比以上三類等效方法，模型三繞組溫度誤差均處在10%以內(nèi)，精度相對較高，該等效方法相對較為合理，滿足要求。

4 結(jié) 語

1) 通過對三種模型分析數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)分析對比發(fā)現(xiàn)定子槽內(nèi)絕緣體不能忽略，它對溫度場有較大的影響，這主要?dú)w結(jié)于絕緣材料傳到系數(shù)很小，形成隔層，熱量不容易散出去，因此建立等效模型時(shí)絕緣層的簡化是十分重要的。

2) 對比三種等效模型可知：通過等效絕緣材料熱傳導(dǎo)系數(shù)和分層簡化方法可以使仿真誤差降低。

3) 根據(jù)實(shí)驗(yàn)研究可得：最接近實(shí)際溫度的為第三種簡化模型的熱場分析，誤差均處在10%以內(nèi)，精度相對較高，該等效方法與實(shí)際接近，簡化方法最為合理。模型二的誤差雖在要求范圍之外，但工程要求范圍相接近，因此在精度要求范圍低的情況下，這種簡化方法可以采用。

4) 通過對永磁電機(jī)定子溫升實(shí)驗(yàn)研究和有限元分析研究對比得出電機(jī)定子熱模型的等效方法，該方法的得出對電機(jī)溫度場有限元研究精度的提高有重要意義，同時(shí)對電機(jī)相關(guān)參數(shù)性能有限元分析等效模型的建立有一定的參考意義。

[1] 魏永田,孟大偉.電機(jī)內(nèi)熱交換[M].北京:機(jī)械工業(yè)出版社,1998.

[2] Motor Design Ltd.Motor-CAD software manual[Z].Motor Design Ltd,2006.

[3] 鮑里先科 A N,丹科 B T.電機(jī)中的空氣動(dòng)力學(xué)與傳熱學(xué)[M].北京:機(jī)械工業(yè)出版社,1985.

[4] 才家剛.電機(jī)試驗(yàn)手冊[M].北京：中國電力出社,1997.

[5] 馬少麗.汽輪發(fā)電機(jī)定子三維溫度場有限元計(jì)算分析[D].保定：華北電力大學(xué)(保定)，2007.

[6] TRIGEOL J F,BERTINY,LAGONOTTE P.Thermal modeling of an induction machine through the associationof two numerical approaches[J].IEEE Transactions on Energy Conversion,2006,21(2):314-323.

[7] 史忠震,楊立,李青青,等.永磁直流無刷電動(dòng)機(jī)三維溫度場分析[J].微特電機(jī),2013,41(6):21-23.

[8] 史忠震,楊立,許貞俊,等.基于ANSYS的直流無刷電機(jī)模態(tài)分析[J].制造業(yè)自動(dòng)化,2013,35(17):119-121.

[9] 趙振奎，李勇，徐永向，等.低速永磁力矩電動(dòng)機(jī)發(fā)熱的仿真分析與實(shí)驗(yàn)測試[J].微特電機(jī),2009,37(11):15-17.

[10] 史忠震.基于ANSYS電機(jī)定子繞組熱模型等效方法的分析研究[J].現(xiàn)代機(jī)械,2015(4):54-57.

[11] 張穎博,王風(fēng)翔,邢軍強(qiáng),等.高速永磁電機(jī)鐵耗計(jì)算方法研究[C]//中國電工技術(shù)學(xué)會(huì)大電機(jī)專業(yè)委員會(huì)2009年學(xué)術(shù)交流會(huì).北京：中國電工技術(shù)學(xué)會(huì)大電機(jī)專委會(huì)，2009.

[12] 任欽.高功率密度永磁電機(jī)磁場與溫度場計(jì)算[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2009.

Finite Element Analysis and Research of Motor Stator Winding Thermal Model Equivalent Method

NIE Hua-wei1，SHI Zhong-zhen1，XU Zhen-jun2

(1.Guizhou Polytechnic College of Communications,Guiyang 550008,China;2.School of Machinery Industry of Guizhou Province,Guiyang 550001,China)

According to the actual distribution of windings,three kinds of motor winding equivalent model was established.The finite element analysis software ANSYS was used to study on temperature field analysis of three kinds of equivalent models.Through temperature rise experiment research in the three kinds of equivalent model the most reasonable winding thermal model. The winding equivalent thermal model method to the improvement of motor temperature field finite element research precision has a certain guiding value,and finite element analysis for the motor model rationality has certain reference significance.

stator winding; equivalent model; study of temperature; finite element

2015-11-13

TM301.4+1

A

1004-7018(2016)10-0030-04

聶華偉(1984-)，男，博士，研究方向?yàn)槟Σ?、潤滑與表面工程。