羅 玲，李 丹，呂曉威，王 震

(西北工業(yè)大學(xué)，陜西西安710072)

0 引 言

盤(pán)式電機(jī)的氣隙磁場(chǎng)是軸向的，又稱(chēng)軸向磁場(chǎng)電機(jī)，具有軸向尺寸短、結(jié)構(gòu)緊湊等優(yōu)點(diǎn)。盤(pán)式無(wú)鐵心永磁同步發(fā)電機(jī)不存在定轉(zhuǎn)子鐵耗，效率高，且消除了齒槽轉(zhuǎn)矩的影響［1］，高磁能積、高剩磁密度、高矯頑力稀土永磁材料的應(yīng)用解決了無(wú)鐵心結(jié)構(gòu)帶來(lái)的氣隙磁密低的問(wèn)題，因此將其應(yīng)用于風(fēng)力發(fā)電領(lǐng)域能實(shí)現(xiàn)微風(fēng)起動(dòng)，提高風(fēng)能利用率［2－3］。

盤(pán)式電機(jī)結(jié)構(gòu)特殊，磁路長(zhǎng)度、等效氣隙長(zhǎng)度及磁路飽和程度等均隨半徑變化，氣隙磁場(chǎng)沿徑向分布不均勻。目前，盤(pán)式電機(jī)電磁場(chǎng)的分析主要采用傳統(tǒng)的等效磁路計(jì)算方法［4］，其中涉及的諸多參數(shù)需憑經(jīng)驗(yàn)選取，計(jì)算精度不高。為了精確地分析盤(pán)式電機(jī)的電磁場(chǎng)分布，需進(jìn)行三維電磁場(chǎng)仿真。然而，三維電磁場(chǎng)仿真建模復(fù)雜，仿真過(guò)程需要占用大量計(jì)算機(jī)資源和時(shí)間，仿真模型的不合理會(huì)導(dǎo)致計(jì)算機(jī)資源不夠或者仿真精度不高。

本文應(yīng)用電磁場(chǎng)仿真軟件MagNet 建立一臺(tái)盤(pán)式無(wú)鐵心永磁同步發(fā)電機(jī)樣機(jī)的三維仿真模型，利用靜態(tài)求解器分析空載磁場(chǎng)分布規(guī)律，采用動(dòng)態(tài)求解器計(jì)算其空載特性，并進(jìn)行樣機(jī)的空載試驗(yàn)，對(duì)實(shí)驗(yàn)結(jié)果與仿真結(jié)果進(jìn)行對(duì)比分析。

1 仿真模型

1.1 樣機(jī)結(jié)構(gòu)及參數(shù)

本文研究的盤(pán)式無(wú)鐵心永磁同步發(fā)電機(jī)為雙邊轉(zhuǎn)子－中間定子結(jié)構(gòu)，如圖1(a)所示。主磁路從一個(gè)轉(zhuǎn)子的N 極出發(fā)，軸向穿過(guò)氣隙和與之相對(duì)的磁極，進(jìn)入另一轉(zhuǎn)子的背鐵并沿周向到達(dá)相鄰磁極，最后沿相同路徑回到第一個(gè)轉(zhuǎn)子的N 極形成閉合回路。

磁極及電樞繞組的尺寸如圖1(b)所示，圖中僅給出了2 個(gè)線圈。樣機(jī)主要參數(shù)如表1 所示，轉(zhuǎn)子背鐵材料為DW310，矩形磁鋼材料為釹鐵硼NNF38SH，其剩磁感應(yīng)強(qiáng)度為1.22 T，矯頑力為907 kA/m。

圖1 樣機(jī)結(jié)構(gòu)參數(shù)示意圖

表1 樣機(jī)參數(shù)

1.2 仿真建模

采用MagNet 進(jìn)行電磁場(chǎng)仿真時(shí)，為提高求解精度，需要細(xì)化網(wǎng)格剖分，對(duì)內(nèi)存為2 G 的計(jì)算機(jī)，以樣機(jī)的整個(gè)三維模型進(jìn)行計(jì)算，常導(dǎo)致計(jì)算機(jī)內(nèi)存溢出。分析樣機(jī)結(jié)構(gòu)不難發(fā)現(xiàn)，其磁場(chǎng)和繞組分布成周期性對(duì)稱(chēng)，對(duì)稱(chēng)周期為一對(duì)極，亦即整個(gè)電機(jī)模型的1 /3，可以建立樣機(jī)的1 /3 模型進(jìn)行電磁場(chǎng)有限元仿真。

利用MagNet 的幾何建模器、材料編輯器可方便地對(duì)電機(jī)進(jìn)行三維建模［5］。磁鋼和轉(zhuǎn)子背鐵的創(chuàng)建只需在平面繪圖的基礎(chǔ)上進(jìn)行簡(jiǎn)單的直線拉伸。但三維繞組的建模較復(fù)雜，需進(jìn)行多維拉伸，拉伸時(shí)首先根據(jù)樣機(jī)繞組的實(shí)際尺寸(圖2(a))，繪制出繞組的拉伸路徑示意圖(圖2(b))，然后根據(jù)示意圖對(duì)坐標(biāo)系進(jìn)行旋轉(zhuǎn)、平移變換，將坐標(biāo)系的XY 平面變換到繞組的拉伸起始面上，并繪制繞組的橫截面(圖2(c))，最后按照拉伸路徑進(jìn)行拉伸，簡(jiǎn)述如下:

(1)如圖2(b)所示，以點(diǎn)1 作為拉伸起始點(diǎn)，逆時(shí)針旋轉(zhuǎn)一小弧度(本文為2°)，繪制出內(nèi)徑處繞組的一半端部;

(2)直線拉伸至繞組外徑(點(diǎn)2)處，順時(shí)針旋轉(zhuǎn)28°，得到繞組外徑處的端部;

(3)直線拉伸繪制出繞組的另外一條有效邊，并順時(shí)針旋轉(zhuǎn)2°，繪出繞組內(nèi)徑處的另外一半端部，但不閉合。

MagNet 中進(jìn)行多維拉伸時(shí)，需在直線拉伸和圓弧拉伸之間加上Blend 進(jìn)行平滑修飾，拉伸后的單個(gè)繞組模型如圖2(d)所示，在建立好一個(gè)繞組模型的基礎(chǔ)上旋轉(zhuǎn)、復(fù)制就可以得到其他兩個(gè)繞組，完整的繞組模型如圖2(e)所示。

將雙邊轉(zhuǎn)子與電樞繞組組合得到如圖2(f)所示的1 /3 仿真模型，同時(shí)創(chuàng)建空氣包作為模型求解的邊界條件，并將氣隙分為兩層，動(dòng)態(tài)仿真時(shí)一層歸為定子，一層歸為轉(zhuǎn)子。

圖2 仿真建模

2 空載磁場(chǎng)

電樞繞組的磁導(dǎo)率與空氣相當(dāng)，計(jì)算空載磁場(chǎng)時(shí)將其作為氣隙處理，軟件在求解過(guò)程中就不對(duì)電樞進(jìn)行重新細(xì)化剖分，可以節(jié)約計(jì)算機(jī)資源，縮短仿真時(shí)間。由于仿真模型為一對(duì)極下的模型，需要設(shè)置偶對(duì)稱(chēng)邊界條件。偶對(duì)稱(chēng)邊界條件的主面為包圍定、轉(zhuǎn)子的空氣包及氣隙與XZ 平面重合的部分，主面旋轉(zhuǎn)120°得到從面，如圖3 所示。為充分利用現(xiàn)有計(jì)算機(jī)資源，提高計(jì)算精度，將網(wǎng)格剖分精度設(shè)定為不大于5 mm。

圖3 靜態(tài)求解邊界條件

該樣機(jī)的空載磁場(chǎng)計(jì)算結(jié)果如圖4 所示。圖4(a)為磁場(chǎng)矢量圖，反映了電機(jī)內(nèi)磁場(chǎng)的走向和分布情況，不難看出磁力線在氣隙中是軸向的，而在轉(zhuǎn)子背鐵中是周向的;圖4(b)為磁密云圖，該圖顯示轉(zhuǎn)子背鐵中的磁通密度最大，僅為1.27 T，小于DW310 的飽和磁通密度1．8 T;圖4(c)為利用Mag-Net 的場(chǎng)量采集器Field Sampler 獲取的氣隙磁通密度幅值隨半徑變化的曲線，可以看出磁鋼平均半徑68．5 mm 處的磁通密度最大為0．36 T，靠近內(nèi)、外徑處的氣隙磁通密度受邊緣效應(yīng)影響而有所下降;圖4(d)為不同半徑處氣隙磁通密度沿軸向分布曲線，不難看出同一半徑氣隙磁通密度沿圓周方向近似為正弦分布;圖4(e)是用場(chǎng)量采集器導(dǎo)出的數(shù)據(jù)繪制的氣隙磁通密度三維分布圖。

圖4 空載磁場(chǎng)

3 空載電壓

對(duì)圖2 的模型進(jìn)行MagNet 動(dòng)態(tài)仿真時(shí)需設(shè)置運(yùn)動(dòng)邊界條件，使得運(yùn)動(dòng)部件在運(yùn)動(dòng)過(guò)程中被空氣包包圍。選擇定子氣隙與轉(zhuǎn)子氣隙的重合面作為運(yùn)動(dòng)邊界的主面，從面由主面旋轉(zhuǎn)120°得到，如圖5所示。

圖5 運(yùn)動(dòng)邊界設(shè)置

若同時(shí)對(duì)一對(duì)極下的三個(gè)線圈進(jìn)行動(dòng)態(tài)求解，則求解區(qū)域太大，導(dǎo)致計(jì)算機(jī)數(shù)據(jù)存儲(chǔ)溢出，仿真被迫終止。因此，每次只針對(duì)一個(gè)線圈進(jìn)行求解獲得每相單個(gè)線圈的空載電壓eφ0，將其乘以3，得到各相電勢(shì)eφ。額定轉(zhuǎn)速500 r/min 時(shí)，三相空載相電壓如圖6 所示，不同轉(zhuǎn)速下空載線電壓如表2 所示。

圖6 500 r·min － 1仿真空載相電壓

表2 不同轉(zhuǎn)速下空載線電壓仿真結(jié)果

4 空載試驗(yàn)及數(shù)據(jù)分析

4.1 空載試驗(yàn)

用一臺(tái)無(wú)刷直流電動(dòng)機(jī)拖動(dòng)該樣機(jī)進(jìn)行空載試驗(yàn)，試驗(yàn)系統(tǒng)框圖如圖7(a)所示，用WT3000 高精度功率分析儀測(cè)量發(fā)電機(jī)輸出線電壓(發(fā)電機(jī)沒(méi)有中線，無(wú)法測(cè)得相電壓)，扭矩傳感儀測(cè)量發(fā)電機(jī)的轉(zhuǎn)速。通過(guò)控制器調(diào)節(jié)原動(dòng)機(jī)的轉(zhuǎn)速，每隔60 r/min記錄一次試驗(yàn)數(shù)據(jù)。轉(zhuǎn)速為500 r/min 時(shí)的空載線電壓波形如圖7(c)所示，不同轉(zhuǎn)速下空載線電壓如表3 所示。

圖7 盤(pán)式電機(jī)空載試驗(yàn)

表3 試驗(yàn)實(shí)測(cè)電機(jī)空載線電壓

4.2 數(shù)據(jù)分析

根據(jù)表2 與表3 的數(shù)據(jù)，繪制仿真空載特性曲線與試驗(yàn)空載特性曲線，如圖8 所示?？梢钥闯觯抡娼Y(jié)果與試驗(yàn)實(shí)測(cè)結(jié)果相比仿真值偏大，但相差不超過(guò)4．4 V。

圖8 空載電壓仿真值與實(shí)驗(yàn)值對(duì)比

5 結(jié) 語(yǔ)

本文基于一臺(tái)盤(pán)式無(wú)鐵心永磁同步發(fā)電機(jī)樣機(jī)結(jié)構(gòu)的對(duì)稱(chēng)性，建立了1 /3 仿真模型，并設(shè)置合理的邊界條件。靜態(tài)求解器求得空載氣隙磁密分布情況:同一半徑的氣隙磁密沿圓周方向呈正弦分布，同一角度下不同半徑處的氣隙磁密大小不同，永磁體平均半徑處的氣隙磁密最大，距平均半徑處距離越大，受端部效應(yīng)越明顯氣隙磁密越小。動(dòng)態(tài)求解器求得額定轉(zhuǎn)速下的空載電壓為110．8 V，與試驗(yàn)值106．4 V 相差4．4 V。仿真得到的空載電壓波形、空載特性曲線均與試驗(yàn)結(jié)果接近。

［1］ 周俊杰，范承志，葉云岳，等．基于斜磁極的盤(pán)式永磁電機(jī)齒槽轉(zhuǎn)矩削弱方法［J］．浙江大學(xué)學(xué)報(bào)(工學(xué)版)，2010．(08):1548－1552．

［2］ Al－Aawar N，Hijazi T M ，Arkadan A A．Design optimization of Axial－Flux Permanent Magnet generator［J］．2010 14th Biennial IEEE Conference on Electromagnetic Field Computation，2010:1．

［3］ 馮勇利，程鵬，李偉力，等．無(wú)槽盤(pán)式永磁風(fēng)力發(fā)電機(jī)有限元分析［J］．防爆電機(jī)，2009，(03):39－43．

［4］ 黃科元，歐金生，黃守道，等．盤(pán)式永磁同步發(fā)電機(jī)的電磁場(chǎng)分析［J］．湖南大學(xué)學(xué)報(bào)(自然科學(xué)版)，2008．(03):41－45．

［5］ ?；萍迹姶艌?chǎng)分析軟件MagNet 的操作說(shuō)明［M］．2007．