MW級永磁同步風(fēng)力發(fā)電機電磁場仿真

岳菁鵬，張新燕

(新疆大學(xué)，新疆烏魯木齊 830047)

0 引 言

近些年來由于風(fēng)力發(fā)電機技術(shù)的迅猛發(fā)展，風(fēng)力發(fā)電機裝機容量在整個電力系統(tǒng)中占據(jù)的比例也越來越重，且風(fēng)能作為一種清潔能源日益被世界各國廣泛應(yīng)用。發(fā)電機是風(fēng)力發(fā)電機組將風(fēng)能轉(zhuǎn)化為電能的重要裝置，大型風(fēng)力發(fā)電機的實際運行是一個非常復(fù)雜的物理過程，為了描述這一物理過程，對電磁場和溫度場的分析以及多場相互耦合的分析是非常必要的［2－6］。傳統(tǒng)的基于“場路”法、等效網(wǎng)絡(luò)法等的分析方法，由于需要假定發(fā)電機的電磁參數(shù)，進(jìn)而建立相應(yīng)的數(shù)學(xué)模型，從而得到發(fā)電機仿真后的各種信息，這樣的結(jié)果與實際情況會存在偏差，不能準(zhǔn)確反應(yīng)發(fā)電機的運行工況［9－10］。文獻(xiàn)［3－4］均利用建立數(shù)學(xué)模型、解析法、圖解法等傳統(tǒng)方法對發(fā)電機在正常運行和故障時電磁場和溫度場進(jìn)行計算，這些方法中有時需用到經(jīng)驗公式，其精度和準(zhǔn)確度不高，計算結(jié)果可能與實際偏差較大。文獻(xiàn)［5］基于有限元方法對同步風(fēng)力發(fā)電機的靜態(tài)電磁場進(jìn)行仿真時采用簡化處理，但只計算了靜態(tài)電磁場。文獻(xiàn)［6］利用有限元方法對風(fēng)力發(fā)電機進(jìn)行了三種方案的設(shè)計，得出當(dāng)電機的極槽匹配關(guān)系不合理時電機磁場和溫度場分布不均勻，且電機損耗較大，熱負(fù)荷較高。文獻(xiàn)［7］得到了雙饋風(fēng)力發(fā)電機的靜、瞬態(tài)磁力線分布圖和反電動勢曲線，其中反電動勢諧波含量較大。

本文通過Ansoft軟件的RMxprt模塊建立永磁同步風(fēng)力發(fā)電機的精確物理模型，并且利用Maxwell 2D有限元分析軟件的電磁場分析計算功能與后處理能力，有效準(zhǔn)確地反映風(fēng)力發(fā)電機在正常工況及故障狀態(tài)下的風(fēng)力發(fā)電機內(nèi)部電磁場變化情況，進(jìn)而為風(fēng)力發(fā)電機的設(shè)計和故障診斷提供科學(xué)的依據(jù)。

1 建立永磁同步風(fēng)力發(fā)電機模型

1.1 Ansoft仿真步驟介紹

利用Ansoft建立電機的有限元模型有兩種方法:一種是在Maxwell 2D模塊里根據(jù)風(fēng)力發(fā)電機結(jié)構(gòu)和尺寸數(shù)據(jù)直接完成物理建模;二是先利用An-soft RMxprt模塊生成風(fēng)力發(fā)電機二維模型，然后將模型導(dǎo)入Maxwell 2D，再利于求解器(2D Transient Solver)進(jìn)行靜態(tài)和瞬態(tài)磁場的求解。本文采用第二種方法，具體流程如圖1所示。

圖1 Ansoft仿真流程圖

1.2 風(fēng)力發(fā)電機基本參數(shù)

本文對1.5 MW三相(Y型連接)永磁風(fēng)力發(fā)電機(以下簡稱PMSG)進(jìn)行建模及仿真分析，主要參數(shù)如表1所示。

表1 PMSG主要參數(shù)

1.3 建立物理模型及設(shè)定材料

根據(jù)風(fēng)力發(fā)電機物理參數(shù)在RMxprt模塊下選擇永磁電機模塊建立物理模型，指定氣隙(Air－gap)屬性為空氣;繞組(Coil)材料屬性為銅;定子鐵心(Stator)及轉(zhuǎn)子軛(Yoke)材料屬性為 DW465－50，該材料滿足指定B－H曲線;永磁體矯頑力Hc=－947000 A/m，剩磁磁感應(yīng)強度 Br=1.25 T，轉(zhuǎn)子磁極采用斜極結(jié)構(gòu)。RMxprt模塊中模型如圖2所示。

圖2 RMxprt模塊中PMSG模型

1.4 邊界條件定義及加載

Maxwell 2D默認(rèn)為自然邊界條件，即不同媒質(zhì)交界面場量的切向和法向邊界條件屬于自然邊界條件，邊界(Boundaries)參數(shù)設(shè)置為0，即假設(shè)邊界無磁場通過。

2 PMSG瞬態(tài)磁場的計算

2.1 磁場求解器基本原理

在靜磁場中，磁場強度滿足安培環(huán)路定律:

根據(jù)麥克斯韋方程:

將式(3)代入式(2)，得:

根據(jù)亥姆霍茲定理在恒定磁場中可得矢量函數(shù)A和磁通密度B的關(guān)系:

將式(5)代入式(4)可得關(guān)于矢量磁位A滿足的方程:

式中:Jz(x，y)為電流流動截面的電流密度;μr為求解域中材料的相對磁導(dǎo)率;μ0為真空中的磁導(dǎo)率;Az(x，y)為矢量磁位z軸分量。

上式理解靜磁場求解器可以在給定的激勵源Jz(x，y)情況下計算出求解域內(nèi)各個點的矢量磁位。

2.2 仿真結(jié)果及分析

2.2.1 氣隙磁密分析

氣隙磁密波形如圖3所示，發(fā)電機在運行時氣隙磁密呈周期性變化，最大為0.7 T。但本文的例子中轉(zhuǎn)子按照廠家提供的數(shù)據(jù)轉(zhuǎn)子磁極采用斜極，可有效減小齒諧波磁場，不容易導(dǎo)致發(fā)電機振動，基本滿足電機優(yōu)化條件，如若出現(xiàn)諧波含量大，可以調(diào)節(jié)轉(zhuǎn)子磁極的形狀和尺寸，來滿足氣隙磁場的要求。

圖3 氣隙磁密波形

2.2.2 磁力線分布及其磁場密度分布云圖

在額定運行狀態(tài)下磁力線分布如圖4所示，磁密分布云圖如圖5所示。由圖可知，在定子槽、氣隙、定子沖片、轉(zhuǎn)子軛之間磁力線較均勻分布，定子齒在靠近磁極的磁力線較密，邊界基本無漏磁。由于在定子齒頂處的磁云最密，在風(fēng)力發(fā)電機運行時要注意該部分溫度過高而引發(fā)故障。

圖4 額定狀態(tài)下的發(fā)電機內(nèi)部磁力線分布

圖5 磁密分布云圖

3 結(jié) 語

通過上述對永磁直驅(qū)同步風(fēng)力機建模和求解過程的介紹，可以看出利用有限元Ansoft軟件可以較直觀全面地完成對大型風(fēng)力發(fā)電機內(nèi)部主要部件的建模。

MW級直驅(qū)永磁風(fēng)力發(fā)電機由于轉(zhuǎn)子極數(shù)多，相應(yīng)定子槽數(shù)多，相應(yīng)氣隙的齒諧波磁場影響大，本方案在設(shè)計電機時轉(zhuǎn)子磁極采用斜極結(jié)構(gòu)可以減少氣隙磁場密度受齒諧波磁場的影響，有效預(yù)防發(fā)電機故障。

利用Maxwell 2D求解器對風(fēng)力發(fā)電機的瞬態(tài)磁場的求解，能反映不同時刻發(fā)電機的磁場分布和磁場密度云圖，為風(fēng)力發(fā)電機的優(yōu)化設(shè)計和故障診斷奠定基礎(chǔ)。

［1］雷亞洲，Lightbody G.國外風(fēng)力發(fā)電導(dǎo)則及動態(tài)模型簡介［J］.電網(wǎng)技術(shù)，2005，29(12):27 －32.

［2］趙博，張洪亮.ANSOFT12在工程電磁場中的應(yīng)用［M］.北京:中國水利水電出版社，2010.

［3］孫宇光，王祥衍，桂林，等.場路偶合法計算同步發(fā)電機定子繞組內(nèi)部故障暫態(tài)過程［J］.中國電機工程學(xué)報，2004，24(1):136－141.

［4］霍菲陽，李偉力，王冬梅.大型水輪發(fā)電機阻尼條數(shù)對電磁參數(shù)和附加損耗的影響［J］.電機與控制學(xué)報，2011(5):89－94.

［5］段曉田，張新燕，張俊.基于ANSYS的直驅(qū)永磁發(fā)電機電磁場仿真［J］.華東電力，2009(12):1997 －2000.

［6］李偉力，程鵬，張美巍，等.1.5MW永磁風(fēng)力發(fā)電機電磁場與溫度場計算與分析［J］.電機與控制學(xué)報，2010(10):52－57.

［7］梅柏衫，劉海華，張金萍.兆瓦級雙饋風(fēng)力發(fā)電機電磁場有限元分析［J］.微電機，2010(10):26－29.

［8］何山，王維慶，張新燕，等.基于有限元方法的大型直驅(qū)同步風(fēng)力發(fā)電機電磁場計算［J］.電網(wǎng)技術(shù)2010(3):157－160.

［9］李俊卿，馬少麗，李和明.汽輪發(fā)電機不同負(fù)荷時電磁場的分析與計算［J］.華北電力大學(xué)學(xué)報，2009(7):1 －5.

［10］李偉清.汽輪發(fā)電機故障檢測及預(yù)防［M］.北京:中國電力出版社，2002.

［11］劉國清.基于ANSYS永磁同步發(fā)電機的電磁分析［J］.機電工程技術(shù)，2009，39(5):68－70.