趙忠霞 張萬全

（1．蘭州石化職業(yè)技術(shù)大學(xué)機(jī)械工程學(xué)院，甘肅 蘭州 730000；2．哈爾濱電氣國際工程有限責(zé)任公司，黑龍江 哈爾濱 150080）

0 引言

無刷勵(lì)磁機(jī)是為同步電機(jī)旋轉(zhuǎn)磁場(chǎng)提供直流電的電機(jī)，分為定子、轉(zhuǎn)子、旋轉(zhuǎn)整流盤和外罩4個(gè)部分。從哈爾濱某電機(jī)廠引進(jìn)原型機(jī)實(shí)體，使用密閉循環(huán)通風(fēng)方式進(jìn)行內(nèi)部冷卻，但是在對(duì)原型機(jī)進(jìn)行工程仿制的過程中，存在風(fēng)路不暢、局部溫升過高的情況。該文對(duì)原型無刷勵(lì)磁機(jī)空氣流場(chǎng)特性進(jìn)行研究，并嘗試對(duì)通風(fēng)系統(tǒng)進(jìn)行改進(jìn)，對(duì)工程實(shí)際指導(dǎo)具有重要意義。目前，關(guān)于封閉空間旋轉(zhuǎn)部件流場(chǎng)及溫度場(chǎng)的研究有很多，但大多數(shù)研究對(duì)象都是1個(gè)旋轉(zhuǎn)部件，無法體現(xiàn)無刷勵(lì)磁機(jī)內(nèi)部流場(chǎng)的運(yùn)動(dòng)規(guī)律。

1 建立物理及數(shù)學(xué)模型

1.1 物理模型

該文利用Gambit三維軟件構(gòu)建某無刷勵(lì)磁機(jī)物理屬性各不同的組及其組合件的幾何結(jié)構(gòu)圖，考慮真機(jī)結(jié)構(gòu)的對(duì)稱性及流場(chǎng)周期性，同時(shí)為了簡化運(yùn)算，選取1/2的整機(jī)作為研究對(duì)象，通過布爾運(yùn)算得到關(guān)于計(jì)算區(qū)域內(nèi)空氣區(qū)、固體區(qū)的三維物理模型，其物理模型簡圖如圖1所示。

圖1 計(jì)算域內(nèi)的三維物理模型

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

無刷勵(lì)磁機(jī)內(nèi)部流體為不可壓縮流體，即=С（為流體密度，kg/m；C為常數(shù)）的流體，空氣流動(dòng)處于紊流狀態(tài)也就是雷諾數(shù)＞2 320（注：2 320為下臨界雷諾數(shù)）。在旋轉(zhuǎn)參考坐標(biāo)系中建立三維流動(dòng)以及傳熱穩(wěn)態(tài)控制方程，轉(zhuǎn)子、整流盤部分采用穩(wěn)態(tài)流動(dòng)控制方程組，包括質(zhì)量守恒方程及動(dòng)量守恒方程式，能量守恒方程式如公式（1）所示。

式中：為散度，1/s；為流體的平均密度，kg/m；為溫度，K；為絕對(duì)速度矢量，m/s；為擴(kuò)散系數(shù)；S為單位體積內(nèi)熱源產(chǎn)生的熱量與定壓比熱c的比值。

針對(duì)穩(wěn)態(tài)流動(dòng)，流體質(zhì)量守恒方程如公式（2）所示。

式中：、和分別為流體在、和方向上的速度，m/s；為密度，kg/m。

勵(lì)磁機(jī)除轉(zhuǎn)子、整流盤外的其他部分位于三維固定直角坐標(biāo)系中，描寫紊流問題的控制方程組也包括質(zhì)量、動(dòng)量與能量守恒方程及標(biāo)準(zhǔn)k-ε兩方程湍流模型。

當(dāng)應(yīng)用CFD軟件計(jì)算時(shí)，必須選中定義Model選項(xiàng)中的剪切加熱（viscous heating）；當(dāng)采用多重參考系定義邊界條件時(shí)，靜止的固體壁面與流體間無須設(shè)置相應(yīng)參數(shù)，但是在旋轉(zhuǎn)的壁面與不旋轉(zhuǎn)的流體界面處，固體壁面Wall必須設(shè)置相應(yīng)的轉(zhuǎn)速及旋轉(zhuǎn)方向；否則，貼壁黏性流體邊界層及由此產(chǎn)生的流體內(nèi)部剪切摩擦生熱計(jì)算模塊不能正確地與數(shù)學(xué)模型中的其他控制方程進(jìn)行有效聯(lián)立；所有流體與固體壁面交界處均為默認(rèn)耦合對(duì)流邊界，不需要設(shè)定對(duì)流換熱系數(shù)（由換熱微分方程求得）。換熱微分方程為固體與流體溫度場(chǎng)、流場(chǎng)聯(lián)系的關(guān)鍵方程。

2 2個(gè)方案的流場(chǎng)特性比對(duì)

2.1 計(jì)算條件

為進(jìn)一步了解所設(shè)計(jì)的勵(lì)磁機(jī)定轉(zhuǎn)子冷卻風(fēng)路是否合理，筆者對(duì)方案一紊流流場(chǎng)和溫度場(chǎng)進(jìn)行三維數(shù)值模擬，通過三維流場(chǎng)數(shù)值模擬反演真機(jī)的實(shí)際運(yùn)行狀態(tài)。方案一假設(shè)在定轉(zhuǎn)子上方的通風(fēng)口為入口，整流盤上方通風(fēng)口為出口，計(jì)算區(qū)域選取整機(jī)物理模型的1/2（圖1）。

計(jì)算條件如下：轉(zhuǎn)子轉(zhuǎn)速為=4 800 r/min，入口壓力為0 MPa，出口壓力為0 MPa。中心剖面固體為固壁邊界，中心剖面的流體區(qū)域?yàn)橹芷谛赃吔纭＝?jīng)計(jì)算入口雷諾數(shù)遠(yuǎn)大于2 320，屬于紊流范圍。某電機(jī)研究所通過電磁場(chǎng)計(jì)算得到各種損耗數(shù)值（表1），由表1可求出計(jì)算域內(nèi)各熱源數(shù)值，入口風(fēng)筒均接近主機(jī)風(fēng)扇前端，需要在計(jì)算三維流場(chǎng)后才能確定出、入口，根據(jù)電動(dòng)機(jī)相關(guān)規(guī)范選取空氣入口溫度為40 ℃。當(dāng)計(jì)算求解域內(nèi)的流場(chǎng)及溫度場(chǎng)時(shí)，貼壁處須滿足標(biāo)準(zhǔn)壁面函數(shù)要求，方程組采用分離及隱式求解，壓力速度耦合方程組采用SIMPLE算法，方程均采用二階迎風(fēng)格式離散，經(jīng)多次網(wǎng)格調(diào)整，方程組獲得網(wǎng)格獨(dú)立收斂解。將無刷勵(lì)磁機(jī)計(jì)算域劃分成固體區(qū)域、流體區(qū)域，僅轉(zhuǎn)子和整流盤內(nèi)流體是旋轉(zhuǎn)空氣，會(huì)隨主軸一起旋轉(zhuǎn)，其他均為靜止流體。

表1 勵(lì)磁機(jī)內(nèi)定轉(zhuǎn)子部件損耗值（單位：W）

在確定計(jì)算域后，選取勵(lì)磁機(jī)入口以及出口邊界條件作為壓力的入出口，相對(duì)壓強(qiáng)均為0 MPa；中心剖面的流體區(qū)域?yàn)橹芷谛赃吔?，固體區(qū)域?yàn)楣腆w壁面，近壁處采用標(biāo)準(zhǔn)壁面函數(shù)，以=4 800 r/min的運(yùn)行工況進(jìn)行計(jì)算，方程求解控制參數(shù)相對(duì)誤差均為1×10，誤差在允許范圍內(nèi)，方程組獲得收斂解。

通過單一因素及方案一了解無刷勵(lì)磁機(jī)內(nèi)壓力、速度以及溫度等物理量的變化規(guī)律。在其基礎(chǔ)上，可以有針對(duì)性地改變通風(fēng)結(jié)構(gòu)，達(dá)到降低峰值溫度及優(yōu)化設(shè)計(jì)的目的。由于方案一中冷卻轉(zhuǎn)子鐵心及線棒的空氣流到轉(zhuǎn)子鐵心及線棒需要的時(shí)間太長，同時(shí)進(jìn)入遠(yuǎn)離整流盤側(cè)轉(zhuǎn)子風(fēng)溝進(jìn)入的風(fēng)量很少，因此冷卻效果不好。因?yàn)閯?lì)磁機(jī)中最重要的是要轉(zhuǎn)子及鐵心不超溫，所以方案二設(shè)計(jì)的入口在定轉(zhuǎn)子遠(yuǎn)離整流盤的一側(cè)，出口在定轉(zhuǎn)子上方。在消化吸收引進(jìn)技術(shù)過程中，設(shè)計(jì)新型通風(fēng)風(fēng)路，通過比對(duì)確定最優(yōu)通風(fēng)方案，方案一的出、入口如圖2所示，方案二的出、入口如圖3所示。

2.2 流場(chǎng)特性比對(duì)分析

為了分析2個(gè)方案旋轉(zhuǎn)流區(qū)及固定流區(qū)內(nèi)的流動(dòng)特性，該文匯總了進(jìn)出口體積流量、計(jì)算域內(nèi)壓力范圍以及速度范圍等參數(shù)，見表2。

由表2可知，出、入口位置的選擇對(duì)風(fēng)路有很大的影響。方案一設(shè)定整流盤上方為入口、定子上方為出口，通風(fēng)風(fēng)路如圖2所示。經(jīng)入口進(jìn)入勵(lì)磁機(jī)的空氣分為3個(gè)部分：第一部分在整流盤及轉(zhuǎn)子旋轉(zhuǎn)作用，在入口處相互削弱，形成短路；第二部分受整流盤的負(fù)壓吸風(fēng)作用，冷卻遠(yuǎn)離整流盤側(cè)轉(zhuǎn)子，最終進(jìn)入定子背部，部分空氣直接通過定轉(zhuǎn)子氣隙進(jìn)入整流盤側(cè)渦旋區(qū)做循環(huán)運(yùn)動(dòng)；第三部分空氣直接進(jìn)入整流盤與轉(zhuǎn)子間的空腔，形成大的渦旋區(qū)，部分空氣經(jīng)通風(fēng)槽（冷卻定轉(zhuǎn)子）后進(jìn)入定子背部，小部分空氣直接從出口流出，大部分空氣重新參內(nèi)部循環(huán)。

表2 2個(gè)方案流場(chǎng)特性對(duì)比

圖2 方案一的通風(fēng)風(fēng)路圖

方案二設(shè)定定轉(zhuǎn)子左側(cè)上方為入口、定子上方為出口，通風(fēng)風(fēng)路如圖3所示。經(jīng)入口進(jìn)入的空氣一部分在冷卻左側(cè)端部線棒后進(jìn)入定轉(zhuǎn)子氣隙中，冷卻轉(zhuǎn)子鐵芯及繞組部件；同時(shí)，在整流盤離心力和科氏力作用下，一部分空氣由整流盤右側(cè)口進(jìn)入定子機(jī)座內(nèi)，最終進(jìn)入定轉(zhuǎn)子氣隙中，經(jīng)出風(fēng)口流出。

圖3 方案二的通風(fēng)風(fēng)路圖

2.3 溫度場(chǎng)特性比對(duì)分析

無刷勵(lì)磁機(jī)內(nèi)部絕緣材料工作溫度不超過130 ℃，二極管工作溫度不超過110 ℃。根據(jù)經(jīng)驗(yàn)可知，二極管發(fā)熱的可能性極小，可忽略。2個(gè)方案勵(lì)磁機(jī)內(nèi)最高溫度均未出現(xiàn)超溫情況，見表3。從溫度場(chǎng)分布可知，方案一的風(fēng)路布置方式會(huì)直接導(dǎo)致勵(lì)磁機(jī)內(nèi)部溫度分布的最大差值為10 ℃，雖然該方案冷卻效果能滿足安全要求，但是冷卻介質(zhì)在勵(lì)磁機(jī)內(nèi)部的循環(huán)流動(dòng)性較差。

表3 勵(lì)磁機(jī)內(nèi)部溫度臨界值分布

由圖4～圖5可知，勵(lì)磁機(jī)內(nèi)最高溫度為116 ℃，位于端壓圈附近的端部轉(zhuǎn)子邊段鐵心中；轉(zhuǎn)子線棒最高溫度達(dá)112℃，位于遠(yuǎn)離整流盤側(cè)的轉(zhuǎn)子邊段鐵心位置轉(zhuǎn)子線棒?？傮w來說，計(jì)算域溫度呈現(xiàn)入口側(cè)溫度低、出口側(cè)溫度高的溫度分布特點(diǎn)。方案一入出口風(fēng)路的布置方式導(dǎo)致定子部件溫度分布左右不均勻，遠(yuǎn)離整流盤側(cè)的部件溫升較大。

圖4 極角10°及170°溫度分布云圖

圖5 轉(zhuǎn)子鐵心及線棒溫度分布云圖

由圖6～圖7可知，方案二中整流盤上方的外罩與機(jī)座間溫度較高（88 ℃～91 ℃）。2個(gè)旋轉(zhuǎn)體中間腔體內(nèi)空氣溫升較高，大部分空氣區(qū)域溫度為98 ℃～106 ℃，也進(jìn)一步證明了該方案大部分流體參與了內(nèi)部冷卻，且冷卻效果較好。

圖6 極角10°與84°截面位置溫度分布云圖

圖7 勵(lì)磁機(jī)內(nèi)空氣區(qū)溫度分布云圖

3 結(jié)語

在通風(fēng)效果方面，方案一入口冷卻空氣體積流量q為0.044 00 m/s，方案二入口冷卻空氣體積流量q為0.326 52 m/s，從進(jìn)風(fēng)空氣量來說，方案二的空氣量約為方案一的10倍，方案二的通風(fēng)效果更好。

在溫度場(chǎng)方面，方案一內(nèi)勵(lì)磁機(jī)內(nèi)最高溫度為116 ℃，方案二內(nèi)勵(lì)磁機(jī)內(nèi)最高溫度為112 ℃，均不超溫，但是方案一呈現(xiàn)兩側(cè)溫度分布極其不均勻的特點(diǎn)。方案二的整流盤空氣區(qū)溫度為98 ℃～106 ℃，非常接近許用值（110 ℃），應(yīng)進(jìn)一步改變冷卻空氣入出口位置及個(gè)數(shù)，從而使整流盤區(qū)空氣順利出流，進(jìn)而降低二極管附近空氣的溫度。