張志弘 韓勤鍇 徐學平 顏益峰 彭 超

基于保角變換與等效磁路法的永磁直驅(qū)發(fā)電機氣隙磁場計算

張志弘1韓勤鍇2徐學平3顏益峰1彭 超1

（1. 新疆金風科技股份有限公司 烏魯木齊 830026 2. 清華大學機械工程系 北京 100084 3. 北京航空航天大學前沿科學技術(shù)創(chuàng)新研究院 北京 100191）

該文開展了負載工況下永磁直驅(qū)發(fā)電機氣隙磁場研究?；诰_域求解和保角變換，計算考慮齒槽效應(yīng)的氣隙磁場標量磁勢；采用離散線電流元，定量模擬發(fā)電機定轉(zhuǎn)子局部出現(xiàn)的磁飽和現(xiàn)象；基于考慮幾何尺寸的定轉(zhuǎn)子非線性磁阻計算方法，建立包含發(fā)電機定轉(zhuǎn)子結(jié)構(gòu)以及氣隙長度的等效磁路，給出相應(yīng)的迭代求解流程。以實際在役兆瓦級永磁直驅(qū)風力發(fā)電機為對象，計算考慮齒槽效應(yīng)的氣隙磁場分布，分析了空載和負載工況下磁飽和效應(yīng)對氣隙磁場分布曲線波形和諧波成分的影響，并與有限元計算結(jié)果對比，驗證了所提出的分析模型和求解方法的有效性。

永磁直驅(qū)發(fā)電機 氣隙磁場 等效磁路 負載工況 磁飽和 齒槽效應(yīng) 非線性磁阻

0 引言

由于低速性能好、發(fā)電效率高、維護成本低、運行可靠等突出優(yōu)點，兆瓦級永磁直驅(qū)風力發(fā)電機是大型陸上及海上風力發(fā)電的理想選擇[1-3]。在永磁直驅(qū)風力發(fā)電機性能分析與優(yōu)化設(shè)計中，準確高效地獲得氣隙磁場分布是重要前提。

盡管有限元法普適性強、計算精度高，但建模和計算時間較長[4]；而計算量小、速度快且物理概念清晰的解析法仍然受到學者的廣泛關(guān)注。由于定子開槽導致氣隙磁通密度發(fā)生畸變，諧波成分增加，所以基于解析法準確計算氣隙磁通密度變得困難。通過經(jīng)驗系數(shù)的修正，采用傳統(tǒng)磁路法可獲得考慮齒槽效應(yīng)的氣隙磁通密度最大值[5]。因定子槽形、定轉(zhuǎn)子軛部尺寸和形狀等參數(shù)的影響，修正系數(shù)難以準確得到，致使考慮齒槽效應(yīng)的傳統(tǒng)磁路法計算精度不高[6]。子域法[7-9]將電機域劃分為氣隙、永磁體和齒槽三個子域，通過交界條件建立各個子域的拉普拉斯方程或泊松方程的聯(lián)系，進而采用解析法求解?；诒＝亲儞Q，將齒槽等不規(guī)則域變換為規(guī)則求解域，利用子域法的解析結(jié)果可獲得氣隙磁場分布[10-11]。與傳統(tǒng)磁路法相比，基于保角變換的解析法計算精度較高。但上述方法均難以考慮因鐵磁材料非線性引起的磁飽和效應(yīng)，因而具有局限性[6]。

磁飽和導致定子齒端部的相對磁導率降低。F.R.Alam和K. Abbaszadeh[12-13]提出一種改進的保角變換，通過增加定子齒槽開口寬度模擬磁飽和效應(yīng)。因經(jīng)驗飽和因子的引入，使得該方法的精度不易保證?；诖怕坊鶢柣舴蚨傻牡刃Т怕贩╗14]，是一種用于電機定量分析的方法，可考慮磁飽和效應(yīng)。但是與有限元法相比，其網(wǎng)格劃分較粗略，雖然計算效率較高，但精度較差。經(jīng)過S.D.Sudhoff等[15]的不斷改進，等效磁路法的計算精度已有很大提升，且可計及電機精細幾何尺寸、三維空間磁場等復雜因素的影響。目前，等效磁路法已廣泛用于多種形式的永磁電機穩(wěn)態(tài)和瞬態(tài)磁場分析中，取得很好的效果。很多學者結(jié)合保角變換和等效磁路法，研究了考慮飽和與否對表面永磁電機氣隙磁場分布的影響[6,16-21]。但這些研究主要針對空載工況，針對負載工況的研究則較少[22]。與空載相比，負載工況下的氣隙磁場變化更為明顯，諧波成分更為復雜，磁飽和效應(yīng)更為顯著。

因此，本文開展負載工況下永磁直驅(qū)發(fā)電機氣隙磁場研究?；诰_域求解和保角變換，計算考慮齒槽效應(yīng)的氣隙磁場標量磁勢；采用離散線電流元，定量模擬發(fā)電機定轉(zhuǎn)子局部出現(xiàn)的磁飽和現(xiàn)象；基于考慮幾何尺寸的定轉(zhuǎn)子非線性磁阻計算方法，建立包含發(fā)電機定轉(zhuǎn)子結(jié)構(gòu)以及氣隙長度的等效磁路，給出相應(yīng)的迭代求解框架。以實際在役兆瓦級永磁直驅(qū)風力發(fā)電機為對象，計算考慮齒槽效應(yīng)的氣隙磁場分布，并與有限元計算結(jié)果對比，驗證所提出的分析模型和求解方法。在此基礎(chǔ)上，分析了空載和負載工況下磁飽和效應(yīng)對氣隙磁場分布曲線波形和諧波成分的影響。

1 氣隙磁場計算模型

1.1 保角變換

典型兆瓦級永磁直驅(qū)風力發(fā)電機結(jié)構(gòu)如圖1所示，為外轉(zhuǎn)子結(jié)構(gòu)。徑向磁化的永磁體，粘貼于轉(zhuǎn)子內(nèi)徑。磁極對數(shù)為，齒槽數(shù)為s。2倍磁極對數(shù)與齒槽數(shù)的最大公約數(shù)為g=gcd(2,s)?？紤]結(jié)構(gòu)周期性并忽略局部細節(jié)尺寸（如定/轉(zhuǎn)子磁軛的散熱結(jié)構(gòu)），建立永磁發(fā)電機對稱分析模型，如圖2所示。其中所考慮的永磁體個數(shù)為2g，齒槽數(shù)為s/g。

圖1 兆瓦級永磁直驅(qū)風力發(fā)電機結(jié)構(gòu)

圖2 永磁發(fā)電機氣隙磁場分析

暫不計入永磁體和負載繞組對標量磁勢的影響。為了提高計算效率，僅分析定子齒鄰近半個齒槽與氣隙的簡化結(jié)構(gòu)，如圖2b所示。其余氣隙和齒槽區(qū)域的解，通過“鏡像對稱”的方式獲得。對于圖2b所示的氣隙結(jié)構(gòu)，在物理坐標系（）下，需要6個頂點來定義。各個頂點的坐標(x,y)為

式中，od、id分別為定子外徑和內(nèi)徑；tfw、sw分別為定子齒寬度和齒槽寬度；pm為永磁體厚度；為氣隙長度，如圖2a所示。

基于保角變換，可將原物理平面上復雜幾何域映射成像平面上的單位矩形平面、半無限平面等。多數(shù)情況下，難以獲得保角變換的解析解，需借助數(shù)值方法。本文利用Matlab中的Schwarz–Christoffel保角變換工具包[11]，將包含齒槽的氣隙區(qū)域變換為規(guī)則的矩形區(qū)域，如圖2c所示。矩形的上下兩邊（1-4和2-3）分別為轉(zhuǎn)子和定子鐵心（磁導率可近似為無窮大），而1-2和3-4則為氣隙邊界。若該矩形區(qū)域存在一個離散電流元，坐標為(v,v)。任意一點(,)的標量磁勢可表示為

1.2 離散線電流元

對于徑向磁化的永磁體，其磁場可通過磁體兩側(cè)的電流片模擬。電流片可用mms個線電流元離散。在這種情況下，單個電流元的幅值為

式中，Θpm為永磁體的磁動勢；Hc為永磁體的矯頑力。將所分析的永磁體分別采用離散電流元進行等效，離散線電流元模擬永磁體、繞組以及定轉(zhuǎn)子飽和效應(yīng)如圖3所示。類似地，采用離散線電流元模擬負載繞組的作用，表示為Iwi（i=1,2,, Qs/Ng）。

1.3 等效磁路

圖4 永磁直驅(qū)發(fā)電機等效磁路模型

式中，fe為電機軸向長度。對于第1個定子齒，注意到邊界條件為第s/g+1個齒的右半部分的磁勢大小相等，但與第1個齒的左半部分的磁勢符號相反，則有

用于模擬磁阻非線性的線電流與穿過鐵心磁阻的磁動勢相同，可從等效磁路模型中獲得。在已知磁通量的前提下，定子齒線電流可以計算為

式中，Ti為定子齒的非線性磁阻，如圖4所示。根據(jù)磁阻的定義，本研究所關(guān)注的定子磁阻可表示為

其中

式中，為定子鐵心磁導率，與磁感應(yīng)強度ΦTi /Sav有關(guān)；幾何參數(shù)如圖5中所示。

定子鐵心的線電流可表示為

其中非線性定子軛部磁阻為

式中，sy和s分別為定子軛部高度和齒槽節(jié)圓跨度。顯然，定子軛部磁阻與磁感應(yīng)強度也是非線性關(guān)系。類似地，轉(zhuǎn)子軛部的線電流為

式中，ryi為轉(zhuǎn)子軛部的非線性磁阻，其表達式為

式中，ry和p分別為轉(zhuǎn)子軛部高度和節(jié)圓跨度。

式中，pm為弧形永磁體的跨度；o,pm、i,pm分別為永磁鐵的內(nèi)徑和外徑；pmi為流經(jīng)永磁體的磁通量，有

因此，在新的迭代計算步中，需將永磁體模擬線電流更新為

1.4 求解流程

基于保角變換與等效磁路法的永磁直驅(qū)發(fā)電機氣隙磁場計算流程如圖6所示，步驟如下。

圖6 永磁直驅(qū)發(fā)電機氣隙磁場求解流程

（1）在給定的初始條件下（包括結(jié)構(gòu)參數(shù)、永磁體參數(shù)、鐵心磁化曲線等），開展Schwarz-Christoffel保角變換，將包含齒槽的氣隙區(qū)域變換至規(guī)則的矩形區(qū)域。

（2）采用離散線電流元模擬永磁體和負載繞組作用，包括mms、wi，以及定子齒端部、定子鐵心和轉(zhuǎn)子鐵心的磁飽和效應(yīng)，包括ti、syi、ryi。

（3）計算規(guī)則矩形區(qū)域內(nèi)的標量磁勢，包括A、pmi、ryi等；構(gòu)建等效磁路模型，計算磁通量源（Ti、syi、ryj等），確定轉(zhuǎn)/定子非線性磁阻（Ti、syi、ryi），得到更新后的用于模擬飽和效應(yīng)的離散線電流元（ti、syi、ryi），以及用于修正永磁體磁動勢的線電流Δpmi。

若更新后的離散線電流元與前次結(jié)果相比，相對變化量大于給定的閾值（比如1×10-4），重復步驟（3）；否則，停止迭代，輸出氣隙磁場分布結(jié)果。需要說明的是，空載時離散線電流元的初值可設(shè)置為零，即ti=syi=ryi= 0；負載時，由于磁飽和效應(yīng)較為突出，可將空載時的離散線電流值作為初值進行迭代，以提高收斂速度。

2 模型驗證

2.1 發(fā)電機參數(shù)

以金風1.5MW永磁直驅(qū)發(fā)電機為研究對象。極對數(shù)=44，齒槽數(shù)s=576，主要結(jié)構(gòu)尺寸可參考文獻[4]。永磁體的矯頑力c=1 018kA/m，相對反沖磁導率為1.031 8。通常僅可獲得定/轉(zhuǎn)子硅鋼片相對磁導率隨磁感應(yīng)強度變化的離散點。為了便于分析計算，需對其進行擬合。根據(jù)C. Heida和P. Regtien[23]的研究，可選擇擬合函數(shù)為

式中，為磁感應(yīng)強度；為擬合系數(shù)。依據(jù)相對磁導率數(shù)據(jù)，采用式（19）進行擬合，得到硅鋼片的相對磁導率結(jié)果如圖7所示。擬合系數(shù)的具體值見表1，擬合置信度為95%。

表1 硅鋼片相對磁導率擬合系數(shù)結(jié)果

Tab.1 Fitting coefficient results of relative permeability of silicon steel sheet

2.2 與有限元結(jié)果對比

構(gòu)建永磁直驅(qū)發(fā)電機的電磁有限元模型，如圖8所示。以1/2氣隙長度處為例，空載條件下本文結(jié)果與有限元結(jié)果的對比如圖9所示。以空載時氣隙最大磁感應(yīng)強度為尺度，對圖中結(jié)果做了無量綱化處理。由于2倍磁極對數(shù)與齒槽數(shù)的最大公約數(shù)為g=8，計算時截取72槽與11塊永磁體所對應(yīng)的扇區(qū)（45°）進行計算。從波形上看，本文結(jié)果與有限元結(jié)果基本保持一致。最大磁感應(yīng)強度出現(xiàn)在永磁體（非齒槽開口處）。由于齒槽的存在，使得氣隙磁感應(yīng)強度出現(xiàn)波動。從幅值上看，本文所得到的最大氣隙磁感應(yīng)強度值稍低于有限元結(jié)果。上述對比驗證了本文模型的正確性。

圖8 有限元模型

圖9 空載條件下本文結(jié)果與有限元結(jié)果的對比（1/2氣隙長度處）

3 分析與討論

3.1 空載工況

考慮定/轉(zhuǎn)子鐵心磁飽和效應(yīng)、永磁體磁導率大于1等因素，重新計算了空載下氣隙磁感應(yīng)強度結(jié)果，如圖10所示。圖中結(jié)果也做了無量綱化處理，尺度為空載氣隙最大磁感應(yīng)強度?？紤]飽和效應(yīng)與否，二者差別很小（不超過0.5%），表明空載條件下定/轉(zhuǎn)子鐵心的磁飽和效應(yīng)較弱，可忽略不計。氣隙磁感應(yīng)強度的諧波次數(shù)和幅值如圖11所示?？梢钥闯觯l率及其奇數(shù)倍頻（c, 3c, 7c, 9c）的幅值較為明顯。此外，齒槽頻率與基波頻率的組合成分也可觀察到，即s±c。

圖10 空載條件下考慮飽和和不考慮飽和時氣隙磁感應(yīng)強度分布

圖11 空載氣隙磁感應(yīng)強度諧波幅值

3.2 負載工況

對于本研究的永磁直驅(qū)發(fā)電機，共有兩套三相繞組。第1套，標記為1U1, 1V1, 1W1，相位差為120°；第2套，標記為2U1, 2V1, 2W1，相位差120°。兩套繞組的相位差為27.5°，聯(lián)結(jié)方式為2Y。已知額定電流的有效值為642A，并聯(lián)支路數(shù)為8。每個齒槽含有上、下兩層繞組，匝數(shù)為10。根據(jù)繞組聯(lián)結(jié)方式，可確定特定時刻（例如=0），上下兩層繞組的電流值分布，如圖12所示。以電流最大值（642A×1.414×8/10=1134.7A）為尺度，對圖12中結(jié)果進行了無量綱處理。

圖12 特定時刻上下兩層繞組的電流值分布

負載工況下考慮磁飽和與否氣隙磁感應(yīng)強度分布結(jié)果如圖13所示。圖中結(jié)果也做了無量綱化處理，尺度為空載氣隙最大磁感應(yīng)強度。與空載時的結(jié)果相比，氣隙磁感應(yīng)強度的變化幅值明顯增大，最大相對磁感應(yīng)強度已超過1.4。磁飽和效應(yīng)的影響也顯著增強，考慮飽和后氣隙磁感應(yīng)強度與不考慮飽和時的結(jié)果相差可達15%以上。由于負載繞組的引入，使得氣隙磁場的變化變得劇烈，在定子齒端部、定子鐵心和轉(zhuǎn)子鐵心出現(xiàn)了不同程度的磁飽和現(xiàn)象。負載氣隙磁感應(yīng)強度的諧波次數(shù)和幅值如圖14所示。諧波次數(shù)與空載時的情況基本一致，包括基波頻率及其奇數(shù)倍頻，以及齒槽頻率與基波頻率的組合成分。幅值要高于空載時的結(jié)果，尤其負載時基波頻率c的幅值已接近0.6。此外，在基波頻率附近，考慮繞組負載后，出現(xiàn)了幅值相對較低的諧波階次，如圖14所示。

圖13 負載條件下考慮飽和和不考慮飽和時氣隙磁感應(yīng)強度分布

圖14 負載氣隙磁感應(yīng)強度諧波幅值

4 結(jié)論

基于保角變換和等效磁路模型，開展了空載和負載工況下永磁直驅(qū)發(fā)電機氣隙磁場研究。以實際在役兆瓦級永磁直驅(qū)發(fā)電機為對象，獲得了考慮齒槽效應(yīng)的氣隙磁場分布，并與有限元計算結(jié)果對比，驗證了所提出分析模型和求解方法的有效性。主要結(jié)論如下：

1）空載時考慮飽和效應(yīng)與否氣隙磁感應(yīng)強度基本不變，表明空載條件下定/轉(zhuǎn)子鐵心的磁飽和效應(yīng)較弱，可忽略不計。諧波成分包括基波頻率及其奇數(shù)倍頻、齒槽頻率與基波頻率的組合成分。

2）負載時考慮飽和后氣隙磁感應(yīng)強度與不考慮飽和時的結(jié)果相差可達15%以上，表明磁飽和效應(yīng)的影響顯著增強。與空載時的結(jié)果相比，氣隙磁感應(yīng)強度的變化幅值明顯增大。諧波成分與空載時的情況基本一致，但幅值要高于空載時的結(jié)果，尤其是基波頻率的幅值增長最為明顯。

本文所提出的分析模型和求解方法，為兆瓦級永磁直驅(qū)風力發(fā)電機磁固耦合分析和性能優(yōu)化設(shè)計提供重要的前提條件。

[1] 王鳳翔. 永磁電機在風力發(fā)電系統(tǒng)中的應(yīng)用及其發(fā)展趨向[J]. 電工技術(shù)學報, 2012, 27(3): 12-24.

Wang Fengxiang. Application and development tendency of PM machines in wind power generation system[J]. Transactions of China Electrotechnical Society, 2012, 27(3): 12-24.

[2] 高劍, 黃守道, 張文娟, 等. 基于變流器控制策略的直驅(qū)永磁風力發(fā)電機優(yōu)化設(shè)計[J]. 電工技術(shù)學報, 2013, 28(7): 103-109.

Gao Jian, Huang Shoudao, Zhang Wenjuan, et al. Optimal design for permanent magnet wind power generators based on converter controlling algorithm[J]. Transactions of China Electrotechnical Society, 2013, 28(7): 103-109.

[3] 許水清, 陶松兵, 何怡剛, 等. 基于相電流瞬時頻率估計的永磁直驅(qū)風電變流器開路故障診斷[J]. 電工技術(shù)學報, 2022, 37(2): 433-444.

Xu Shuiqing, Tao Songbing, He Yigang, et al. Open-circuit fault diagnosis for back-to-back converter of PMSG wind generation system based on estimated instantaneous frequency of phase current[J]. Transactions of China Electrotechnical Society, 2022, 37(2): 433-444.

[4] 何山, 王維慶, 張新燕, 等. 基于有限元方法的大型永磁直驅(qū)同步風力發(fā)電機電磁場計算[J]. 電網(wǎng)技術(shù), 2010, 34(3): 157-161.

He Shan, Wang Weiqing, Zhang Xinyan, et al. Electromagnetic field calculation of high capacity direct-driven permanent magnet synchronous wind power generator based on finite element method[J]. Power System Technology, 2010, 34(3): 157-161.

[5] 唐任遠. 現(xiàn)代永磁電機: 理論與設(shè)計[M]. 北京: 機械工業(yè)出版社, 2016.

[6] 龐古才, 鄧智泉, 張忠明. 基于改進廣義磁路法的表貼式永磁電機空載氣隙磁場解析計算[J]. 電工技術(shù)學報, 2019, 34(22): 4623-4633.

Pang Gucai, Deng Zhiquan, Zhang Zhongming. Analytical calculation of no-load air gap magnetic field in surface-mounted permanent magnet motor based on improved generalized magnetic circuit method[J]. Transactions of China Electrotechnical Society, 2019, 34(22): 4623-4633.

[7] 范堅堅, 吳建華. 極間隔斷Halbach型磁鋼的永磁同步電機氣隙磁場解析計算及參數(shù)分析[J]. 電工技術(shù)學報, 2010, 25(12): 40-47.

Fan Jianjian, Wu Jianhua. Analytical calculation and parameter analysis of air-gap magnetic field distribution in permanent magnet synchronous motor with partition-between-poles Halbach magnet[J]. Transactions of China Electrotechnical Society, 2010, 25(12): 40-47.

[8] 章躍進, 李春江, 屠關(guān)鎮(zhèn), 等. 面貼式永磁力矩電機氣隙主磁場解析數(shù)值分析法[J]. 電工技術(shù)學報, 2011, 26(9): 13-17.

Zhang Yuejin, Li Chunjiang, Tu Guanzhen, et al. Analytical method for air-gap main magnetic field computation of surface mounted permanent magnet torque motors[J]. Transactions of China Electrotechnical Society, 2011, 26(9): 13-17.

[9] 李節(jié)寶, 井立兵, 周曉燕, 等. 表貼式永磁無刷電機直接解析計算方法[J]. 電工技術(shù)學報, 2012, 27(11): 83-88.

Li Jiebao, Jing Libing, Zhou Xiaoyan, et al. Exact analytical method for surface-mounted permanent-magnet brushless motors[J]. Transactions of China Electrotechnical Society, 2012, 27(11): 83-88.

[10] Zarko D, Ban D, Lipo T A. Analytical calculation of magnetic field distribution in the slotted air gap of a surface permanent-magnet motor using complex relative air-gap permeance[J]. IEEE Transactions on Magnetics, 2006, 42(7): 1828-1837.

[11] O'Connell T C, Krein P T. A schwarz–christoffel-based analytical method for electric machine field analysis[J]. IEEE Transactions on Energy Conversion, 2009, 24(3): 565-577.

[12] Alam F R, Abbaszadeh K. Magnetic field analysis in eccentric surface-mounted permanent-magnet motors using an improved conformal mapping method[J]. IEEETransactions on Energy Conversion, 2016, 31(1): 333-344.

[13] Abbaszadeh K, Alam F R. On-load field component separation in surface-mounted permanent-magnet motors using an improved conformal mapping method[J]. IEEE Transactions on Magnetics, 2016, 52(2): 1-12.

[14] Ostovi? V. Dynamics of Saturated Electric Machines[M]. New York: Springer, 1989.

[15] Sudhoff S D, Kuhn B T, Corzine K A, et al. Magnetic equivalent circuit modeling of induction motors[J]. IEEE Transactions on Energy Conversion, 2007, 22(2): 259-270.

[16] Lim D K, Yi K P, Woo D K, et al. Analysis and design of a multi-layered and multi-segmented interior permanent magnet motor by using an analytic method[J]. IEEE Transactions on Magnetics, 2014, 50(6): 1-8.

[17] 程明, 文宏輝, 曾煜, 等. 電機氣隙磁場調(diào)制行為及其轉(zhuǎn)矩分析[J]. 電工技術(shù)學報, 2020, 35(5): 921-930.

Cheng Ming, Wen Honghui, Zeng Yu, et al. Analysis of airgap field modulation behavior and torque component in electric machines[J]. Transactions of China Electrotechnical Society, 2020, 35(5): 921-930.

[18] 戈寶軍, 毛博, 林鵬, 等. 無刷雙饋電機轉(zhuǎn)子偏心對氣隙磁場的影響[J]. 電工技術(shù)學報, 2020, 35(3): 502-508.

Ge Baojun, Mao Bo, Lin Peng, et al. Effect of rotor eccentricity fault on air gap magnetic field in brushless doubly-fed machine[J]. Transactions of China Electrotechnical Society, 2020, 35(3): 502-508.

[19] 堯磊, 秦雪飛, 蔡順, 等. 內(nèi)置式永磁電機轉(zhuǎn)子硅鋼片疊裝錯位對氣隙磁場的影響[J]. 電工技術(shù)學報, 2021, 36(5): 1096-1100.

Yao Lei, Qin Xuefei, Cai Shun, et al. Influence of rotor lamination stacking misalignment on airgap field in interior permanent magnet motors[J]. Transactions of China Electrotechnical Society, 2021, 36(5): 1096-1100.

[20] 高鋒陽, 齊曉東, 李曉峰, 等. 不等寬不等厚Halbach部分分段永磁同步電機電磁性能解析計算與優(yōu)化分析[J]. 電工技術(shù)學報, 2022, 37(6): 1398-1414.

Gao Fengyang, Qi Xiaodong, Li Xiaofeng, et al. Analytical calculation and optimization analysis of electromagnetic performance of Halbach partially-segmented permanent magnet synchronous motors with unequal width and thickness[J]. Transactions of China Electrotechnical Society, 2022, 37(6): 1398-1414.

[21] Hanic A, Zarko D, Hanic Z. A novel method for no-load magnetic field analysis of saturated surface permanent-magnet machines using conformal mappingand magnetic equivalent circuits[J]. IEEE Transactions on Energy Conversion, 2016, 31(2): 740-749.

[22] Hanic A, Zarko D, Kuhinek D, et al. On-load analysis of saturated surface permanent magnet machines using conformal mapping and magnetic equivalent circuits[J]. IEEE Transactions on Energy Conversion, 2018, 33(3): 915-924.

[23] Heida C, Regtien P. Estimation of permeability for magnetic flux leakage modeling[C]//10th IMEKO TC7 International Symposium on Advances of Measurement Science, Saint-Petersburg, 2004: 321-326.

Air Gap Magnetic Field Calculation of Permanent Magnet Direct Drive Generator Based on Conformal Mapping and Magnetic Equivalent Circuit Method

Zhang Zhihong1Han Qinkai2Xu Xueping3Yan Yifeng1Peng Chao1

(1. Xinjiang Goldwind Science Technology Co. Ltd Urumqi 830026 China 2. Department of Mechanical Engineering Tsinghua University Beijing 100084 China 3. Research Institute for Frontier Science & Technology Beihang University Beijing 100191 China)

Owing to its outstanding advantages such as good low-speed performance, high power generation efficiency, low maintenance cost and reliable operation, megawatt permanent magnet direct drive wind turbine is an ideal choice for large-scale onshore and offshore wind power generation. In the performance analysis and optimal design of permanent magnet direct drive wind turbines, accurate and efficient air gap magnetic field distribution is an important premise.

In this paper, the air gap magnetic field of permanent magnet direct drive generator under load condition is studied. Based on the exact domain solution and conformal mapping, the scalar magnetic potential of the air gap magnetic field considering the cogging effect is calculated. Discrete line current element is used to quantitatively simulate the local magnetic saturation of generator stator and rotor. Based on the calculation method of the nonlinear reluctance of the stator and rotor considering the geometric size, the magnetic equivalent circuit model including the structure of the stator and rotor and the length of the air gap is established, and the corresponding iterative solution process is presented. Taking the actual megawatt permanent magnet direct drive generator in service as an example, the air gap magnetic field distribution considering the cogging effect is calculated, and the influence of magnetic saturation effect on the air gap magnetic field distribution curve waveform and harmonic components under no-load and load conditions is analyzed. The proposed analysis model and solution method are verified by comparing with the finite element results.

Under the no-load condition, the air gap magnetic field intensity is basically unchanged whether the saturation effect is considered or not. The harmonic component includes the fundamental frequency and its odd frequency multiples, and the combination of the tooth slot frequency and the fundamental frequency. Under the load condition, the relative difference between the magnetic field intensity with and without saturation effects can exceed 15%, indicating that the influence of magnetic saturation effect is significantly enhanced and cannot be ignored in the load condition. Compared with the result in the no-load condition, the variation amplitude of the air gap magnetic field intensity increases obviously. The harmonic component is basically the same as that in no-load condition, but the amplitude is higher than that in no-load condition, especially for the increase of the amplitude of fundamental wave frequency.

The analysis model and solution method proposed in this paper provide important preconditions for the magneto-solid coupling analysis and performance optimization design of megawatt permanent magnet direct drive wind turbines.

Permanent magnet direct drive generator, air gap magnetic field, magnetic equivalent circuit, load condition, magnetic saturation, cogging effect, nonlinear reluctance

10.19595/j.cnki.1000-6753.tces.210769

TM351

國家自然科學基金資助項目（11872222）。

2021-05-27

2021-10-22

張志弘 男，1982年生，工程師，研究方向為永磁直驅(qū)風電機組動力學、發(fā)電機磁固耦合分析。E-mail：zhangzhihong@goldwind.com.cn

韓勤鍇 男，1982年生，博士，副教授，研究方向為旋轉(zhuǎn)機械動力學、永磁直驅(qū)電機磁固耦合分析。E-mail：hanqinkai@mail.tsinghua.edu.cn（通信作者）

（編輯 赫蕾）