基于永磁體結構優化的動態渦流抑制方法研究

耿惠 李永建 栗浩森 楊占力

摘要 針對永磁電機的時間諧波和空間諧波引起電機中永磁體渦流損耗增加的問題，本文搭建了永磁體綜合磁特性測試裝置測量永磁體諧波激勵下的動態渦流損耗，從磁特性角度解釋了諧波產生大量渦流損耗的原因;又分別從理論分析、數值計算和磁特性測量實驗驗證的角度研究了抑制永磁體渦流損耗的方法-永磁體分割法。首先，分析了永磁體分割法抑制渦流損耗的原理。然后，建立了不同分割方式的永磁體三維有限元仿真模型，分析永磁體電流密度和渦流損耗隨分塊數的變化規律。最后，應用永磁體磁特性測量系統測量了不同分割方式永磁材料釹鐵硼（NdFeB）的動態磁特性。結果證明在一定分割塊數范圍內，可以有效降低永磁體渦流損耗的80%。實際工況下，應根據永磁體有效尺寸及永磁材料利用率合理選擇分塊數。研究結果對永磁電機永磁體的結構優化設計有重要意義。

關 鍵 詞 永磁電機;諧波;永磁材料;渦流損耗;磁特性測試;永磁體分割

中圖分類號 TM351? ? ?文獻標志碼 A

文章編號：1007-2373（2021）06-0030-08

Abstract To solve the problem of increased eddy current loss in permanent magnet machines caused by time and space harmonics， this article builds a permanent magnet comprehensive magnetic properties test device to measure the permanent magnet dynamic eddy current loss under harmonic excitation， and explains the reason for the large amount of eddy current loss caused by harmonics from the perspective of magnetic properties. The article also analyzes the method to reduce eddy current loss of permanent magnet—permanent magnet segmentation from the perspective of theoretical analysis， numerical calculation and experimental verification. First， the principle of the permanent magnet segmentation method to suppress eddy current loss is analyzed. Then， 3-D finite element models of different segmentation methods of NdFeB are established to study the change of the current density and eddy current loss with segmentation numbers of permanent magnet. Finally， the dynamic magnetic properties of NdFeB permanent magnet materials adopting different segmentation methods are measured using the properties measurement system of permanent magnet. The results show that the eddy current loss of the permanent magnet can be effectively reduced by 80% under specific segmentation methods. Under actual working conditions， the number of blocks should be selected reasonably according to the effective size and the utilization rate of the permanent magnet. The results are of great significance in the optimization design of permanent magnet machines.

Key words permanent magnet machines; harmonic; permanent magnetic material; eddy current loss; the magnetic properties measurement; permanent magnet segmentation

0 引言

稀土永磁電機以其效率高、功率大、結構簡單、節能效果顯著等優點在工業生產和日常生活中得到廣泛應用[1]。近年來，高耐熱性、高磁能積釹鐵硼永磁體的成功開發，使新型稀土永磁電機的研發取得突破性進展[2-3]。

由于永磁電機采用電力電子變換器驅動，會產生大量的時間諧波，電機分布繞組會在氣隙中產生空間諧波，定子開槽引起的氣隙磁導變化會在電機氣隙中產生齒諧波，這些諧波會在永磁體中引起大量渦流損耗[4-6]。如內轉子大功率永磁同步電機工作在額定功率450 kW下，渦流損耗高達2.4 kW，在總電機損耗中的占比可達20%[7]。受永磁電機轉子體積和散熱的限制，渦流損耗引起的溫升會導致永磁體退磁或失磁，降低了電機運行的安全性與可靠性[8-9]。因此，研究諧波是如何對永磁體動態渦流損耗產生影響以及降低永磁體中的渦流損耗具有重要意義。

受永磁電機結構的限制，從各類諧波的抑制或消除的角度研究相對困難，因此，可以從永磁體尺寸與渦流損耗的定量關系考慮如何降低渦流損耗。由于最大渦流損耗的位置在永磁體的邊線上，文獻[10]提出了將永磁同步電機中永磁體進行部分分割的結構優化方法，渦流損耗降低且電機的運行性能顯著提升。文獻[11]通過對永磁體渦流損耗函數求偏導數得出其最小值點，證明永磁體對稱分割時，渦流損耗的抑制效果最好。文獻[12-13]建立了電機二維和三維有限元模型預測軸向分割結構優化后永磁體渦流損耗的變化，結果表明軸向分割永磁體是降低渦流損耗的有效手段。文獻[14]建立了永磁同步電機軸向和周向上對永磁體分割的三維有限元模型，證明了軸向和周向分割永磁體都可以有效降低渦流損耗。上述文獻研究大多是從理論或者仿真的角度分析永磁體分割法對其渦流損耗的影響，部分采用永磁樣機進行驗證。而從永磁材料磁特性測量的角度系統分析分割法對降低永磁體渦流損耗的影響卻鮮有報導。

本文重點從永磁材料磁特性測量[15-16]的角度研究諧波激勵對永磁體渦流損耗的影響及定量分析了永磁體分塊數與渦流損耗的關系。首先，利用永磁體磁特性測試裝置測量不同諧波激勵下永磁體動態磁特性的變化情況，分析諧波產生大量渦流損耗的原因。然后，對永磁體結構進行優化抑制其渦流損耗，建立了實驗樣品鍍鋅釹鐵硼的渦流損耗計算模型，計算釹鐵硼永磁體渦流損耗與樣品體積的關系;建立三維仿真模型對釹鐵硼永磁體實現不同分割方式的仿真，分析渦流損耗與永磁體分塊數的關系。最后對釹鐵硼樣品進行結構優化，應用永磁測試系統，測試優化后永磁體動態磁特性。通過分析實驗數據，定量討論了相同體積下永磁體分塊數量對降低渦流損耗的影響，并給出了不同優化方式下渦流損耗的降低比例。

1 永磁體綜合磁特性測試

1.1 磁特性測試系統

永磁體磁特性測試系統包括激磁裝置、樣品、磁場強度傳感線圈（H線圈）、磁通密度傳感線圈（B線圈）、信號發生器、功率放大器、水冷電阻和匹配電容箱。其包含的設備型號如表1所示。

測試系統如圖1所示，其中激磁裝置主要是由磁極、繞置在磁極上的激磁繞組和樣品組成，裝置中磁極是由高取向超薄硅鋼片疊置而成，可以產生較強的磁場，激磁裝置實現對樣品激磁和感應信號輸出的功能。選取30 mm×10 mm×10 mm的長方體鍍鋅釹鐵硼樣品作為測試對象。樣品被分割之后，其接觸面用耐高溫絕緣紙隔開，防止相鄰樣品間發生短路。H傳感線圈緊貼在樣品表面，實現樣品表面磁場強度的測量。B傳感線圈采用線徑為0.2 mm的漆包線纏繞在樣品四周，實現對樣品內部的平均磁密的測量。功率放大器用于驅動激磁線圈，以保證待測樣品在給定頻率下可控磁化。匹配電容箱用來消除激磁繞組感抗的影響，使電容和電感串聯諧振，整個電路呈阻性，便于激磁。

在交變電流激勵條件下，永磁材料產生磁滯回線。對于磁性材料的損耗通常是計算整個磁滯回線的面積，它是由樣品表面的磁場強度H和樣品內的平均磁通密Bav形成的[17]。公式如下：

式中：P為總損耗;V為樣品的體積;f 為激磁頻率。測量的總損耗包含渦流損耗和磁滯損耗。對于釹鐵硼永磁材料，它的磁導率接近于空氣的磁導率[μ0]，一般認為磁滯損耗幾乎可以忽略不計。由于永磁材料較高的電導率，故此永磁材料測試系統測得的結果近似為永磁材料的渦流損耗。

1.2 磁特性測試系統

1.2.1 諧波磁場條件下的激勵模型

式中：[B1]為基波幅值;i為諧波次數;[Bi]和[θi]為各次諧波幅值和相角。若[Bi]為零，得到的波形即為正弦波。通過控制公式中的變量得到不同的諧波激勵模型，進而對諧波激勵對磁特性的影響進行分析。

1.2.2 諧波對永磁體動態磁特性的影響

設定諧波的含量為50%不變，改變諧波次數i（i=1，3，5，7），激磁模型如式（3）所示，觀察諧波次數對磁特性的影響，相應的磁滯回線如圖2所示。

為了更加直觀的分析諧波含量改變對于磁心損耗的影響，繪制不同諧波次數條件下損耗曲線，如圖3所示。整體來看永磁體的動態渦流損耗隨基波疊加諧波次數的增加而增加。這是因為隨著諧波次數增加，導致永磁體磁滯回線波形中出現了個數更多的小磁通翻轉，這些翻轉使磁滯回線的等效面積變大導致了更多的損耗。諧波次數越大，小磁通回環個數越多，渦流損耗越大。

2 永磁體結構優化設計

2.1 永磁體分割法理論分析

正處于交變磁場中的導體內會感生渦流，交變磁場的頻率越高，渦流越大。由渦流引起的渦流損耗是永磁體主要的發熱源。研究永磁體渦流損耗的分布規律可以有助于降低永磁體發熱。本文首先研究永磁體的電流密度分布規律并基于此建立如圖4所示的渦流損耗計算模型，在該模型中，永磁體的長為a，高為b，寬為c。圖4中的體積微元代表（x，y）坐標處的某回路，永磁體內的磁場為z方向且隨時間正弦變化。本文所采用的頻率最高為2 kHz，在此頻率下釹鐵硼材料的集膚深度仍然遠大于實驗采用永磁體的尺寸，故集膚效應的影響可以忽略。

對于永磁體渦流損耗，可以應用渦流密度分布 J直接計算得到，計算公式為

式中：J為電流密度;E為磁場強度;ρ為磁性材料的電阻率;V為永磁體的體積。

根據法拉第電磁感應定律得：

式中：Ψ為磁通量;B為磁感應強度。

當渦流路徑右上頂點坐標為（x，y）時，渦流按圖1中微元路徑流通：

式中：[Jx]為x方向上的電流密度;[Jy]為y方向上的電流密度。

由電流連續性定理得：

由以上公式可知，塊狀永磁體電流密度的分布與其所處的位置（x， y）、磁通密度的變化率dB/dt及永磁體電阻率ρ有關。在永磁體電阻率及磁通密度的變化恒定時，改變永磁體的尺寸，優化永磁體的結構，即可間接改變其所處的位置（x， y），可以改變永磁體的電流密度分布。結合圖4b）及式（9）、式（10）可知，當把永磁體分塊時，Jx、Jy會變小，電流密度J也將減小。在相同體積下，根據渦流損耗求解式（1），電流密度J減小，其渦流損耗也減小。從原理上講，將永磁體分塊之后，原有的渦流流通路徑被阻斷，渦流被限制在多個分塊的狹小回路中流通，等效回路電阻變大。實驗過程中，穿過永磁體的磁通密度不變，故其感應電動勢不變，總體渦流損耗減小。

2.2 永磁體渦流損耗有限元計算

如圖5所示，利用仿真軟件Ansys進行了永磁體結構優化對降低渦流損耗的定量分析，對激磁磁路建立了1∶1三維模型。未分割的永磁體樣品的截面設置為10 mm×10 mm，其高度設置為30 mm，并在高度上將整塊永磁體樣品分別均勻分割為多個小塊。為了模擬絕緣，塊間設置0.05 mm氣隙。圖5a）為測試裝置的仿真模型，由一對雙C型結構的磁路，激磁繞組以及樣品組成。激磁繞組繞在磁路上，通過磁路在磁軛中間產生均勻磁場，樣品安裝在磁路中間。圖5b）為不同測試樣品的仿真模型。不同分割方式樣品施加相同激勵。

取不同優化方式下同一時刻同一截面上的電流密度J進行比較，仿真結果如圖6所示。可以看出，永磁體的電流密度分布不均勻，電流密度從永磁體中心到永磁體邊界呈現增大趨勢，與理論分析一致。仿真結果表明，隨著永磁體分塊數增多，永磁體的電流密度減小。

根據裝置仿真模型，計算不同分割方式時永磁體的渦流損耗。計算結果如表2所示。由結果可知，永磁體的渦流損耗隨著分塊數增大而減小。

釹鐵硼在不同分塊數量下的渦流損耗相對于不分割情況渦流損耗下降百分比δ計算公式如下：

式中：[P1]為不分割時對應的渦流損耗;[Pn]為相同激磁條件下;分塊數為n時對應的渦流損耗，n>1。

相鄰分塊數下釹鐵硼的渦流損耗的下降幅度變化Δδ計算公式如下：

根據式（11）、式（12）對仿真的渦流損耗數據進行處理，結果如圖7所示。

根據不同分塊數下釹鐵硼的渦流損耗相對于不分塊時渦流損耗的下降百分比δ線可得，將同體積釹鐵硼樣品均勻分割為n （n = 2，3，4，5，6，7，8，9，10）塊時渦流損耗相對于不分塊樣品依次下降了27.94%，46.88%，60.47%，72.00%，78.60%，82.90%，86.24%，89.22%，91.00%。根據相鄰分塊數下釹鐵硼的渦流損耗的下降幅度變化Δδ線可得，將同體積釹鐵硼樣品均勻分割為n （n =2，3，4，5，6，7，8，9，10）塊時渦流損耗相對于n-1塊方式依次下降了27.94%，18.94%，13.59%，11.53%，6.60%，4.30%，3.34%，2.58%，1.78%。由此可以得出兩個結論：一是釹鐵硼材料的渦流損耗隨著分塊數量的增加而不斷下降。二是隨著分塊數量的增加，釹鐵硼材料的渦流損耗的下降幅度Δδ逐漸變小。由于永磁體有效尺寸及利用率等因素的存在，永磁體分塊數不能無限增大。從以上數據可以看出，永磁體分塊達到6時，其渦流損耗已下降78.6%;永磁體分塊大于6時，δ隨分塊數變化已不明顯，Δδ已下降到5%以下。

2.3 永磁體結構優化對降低渦流損耗影響實驗驗證

2.3.1 不同優化方式下永磁體動態B-H曲線

如圖8所示，為不同磁場強度下，f = 2 kHz時，未分割鍍鋅釹鐵硼的動態磁滯回線。從圖中可以看出，隨著磁場強度的增大，磁滯回線的面積明顯增大，即樣品的渦流損耗隨激磁磁場的增大而增大。

2.3.2 不同優化方式下永磁體動態磁特性測量

圖9給出了f=2 kHz，激磁磁場H為10 kA/m時，樣品不分塊以及分3塊情況下的磁滯回線。從圖中可以看出，樣品不分塊時的磁滯回線的面積大于樣品分3塊時磁滯回線的面積，即樣品不分塊時的渦流損耗大于樣品分3塊時的渦流損耗。

2.3.3 不同優化方式下永磁體渦流損耗對比

永磁電機運轉時，實際上使永磁體產生渦流損耗的激磁磁場不一定是正對著永磁體的充磁方向，因為實際的永磁體形狀各異，引起渦流的磁場可以是各個方向。因此本節分別從充磁方向和非充磁方向研究永磁體結構優化對永磁體渦流損耗的影響。在頻率f分別為100 Hz、500 Hz、1 kHz、2 kHz時，比較不同分塊數樣品的渦流損耗隨激磁磁場H的變化曲線。圖10 a）～d）展示了不同頻率時充磁方向上分塊數為n （n=2，3，4，5，6，10）時，釹鐵硼渦流損耗隨磁密的變化趨勢。結果表明，永磁體的渦流損耗隨著分塊數的增加而減小，但是從圖中可以看出隨著分塊數的增加渦流損耗減小的幅度變小。而且隨頻率的不同，其優化的效果基本相同。

圖11 a）～d）展示了不同頻率時非充磁方向上分塊數為n （n=2，3，4，5，6，10）時，釹鐵硼渦流損耗隨磁密的變化趨勢。從圖中可以看出永磁體的渦流損耗隨著分塊數的增加而減小，但是同樣隨著分塊數的增加減小的幅度變小，與仿真分析一致。

由圖10、圖11可以得到隨著分塊數的增加渦流損耗減小的幅度變小，由于在分割后每塊之間是絕緣的，在分界面處的渦流路徑增加，因此采用分割方式的永磁材料的渦流回路電阻增加，故渦流損耗減小。當分塊數量不斷增加時，每增加一個分塊，在分界面處的渦流路徑只增加兩條，因此分割方式下的渦流電阻增加的變化率不斷減小，導致渦流損耗的降低量并不是線性減小的，而是非線性減小，其減小幅度不斷降低。在分塊數量比較大時，分塊數量的小幅改變不會引起渦流損耗的明顯降低。

圖12給出了充磁方向上在頻率為2 kHz時，鍍鋅釹鐵硼樣品在不同分塊數下相對于不分塊情況下渦流損耗下降的百分比。從圖中可以看出，相比于不分塊樣品，當釹鐵硼分塊數為n （n=2，3，4，5，6，10）塊時渦流損耗下降百分比依次在30%，45%，60%，70%，77%，90%左右穩定波動。與表1中的仿真結果具有很好的一致性。

圖13給出了非充磁方向上在頻率為2 kHz時，鍍鋅釹鐵硼樣品在不同分塊數下相對于不分塊情況下渦流損耗下降的百分比。相比于不分塊樣品，當釹鐵硼分塊數為n （n=2，3，4，5，6，10）塊時渦流損耗下降百分比依次在30%，45%，60%，70%，80%，90%左右穩定波動。與表2中的仿真結果同樣具有很好的一致性。

3 結論

本文主要對永磁體渦流損耗形成的主要因素之一諧波進行了測量分析，又從理論計算，模型仿真，磁特性測量實驗驗證的角度定量研究了永磁體分塊數與損耗抑制的關系，進行了永磁體結構優化設計。得出以下結論：

1）通過對永磁材料諧波激勵下的交變磁特性進行測量，說明了注入諧波對磁特性的影響，諧波次數會影響小磁滯回線出現的個數，進而增大永磁體渦流損耗。

2）理論推導，仿真及測量結果均表明，永磁體樣品的電流密度隨分塊數的增加而減小，進而導致渦流損耗也隨分塊數的增加而減小。但隨著分塊數的增多，渦流損耗減小的幅度減緩。即在一定分割塊數范圍內，永磁體分割法可以有效抑制永磁體的渦流損耗。

3）分別從永磁體充磁方向和非充磁方向上，對永磁體結構優化設計進行分析，結果表明兩個方向的優化對渦流損耗的抑制效果基本相同。且在同一充磁方向上，優化效果基本不受頻率和激磁磁場的影響。

在實際電機設計中，永磁體分塊數不能無限增大，永磁體分塊數的具體優化應根據永磁體的工作點的選擇、損耗允許范圍及永磁體工藝實現情況綜合考慮來抑制永磁體渦流損耗。其結果對于永磁電機永磁體結構優化設計的和改進具有重要應用意義。

參考文獻：

[1]? ? 王要強，馮玉濤，馬小勇，等. 永磁同步電機轉子位置復合檢測及起動運行策略[J]. 儀器儀表學報，2019，40（8）：230-238.

[2]? ? 陳吉. 繞組短路故障對表貼式永磁同步電機失磁影響的研究[D]. 重慶：重慶大學，2011.

[3]? ? 楊慶新，李永建. 先進電工磁性材料特性與應用發展研究綜述[J]. 電工技術學報，2016，31（20）：1-29.

[4]? ? YANG C X，ZHANG Y，QIU H B. Influence of output voltage harmonic of inverter on loss and temperature field of permanent magnet synchronous motor[J]. IEEE Transactions on Magnetics，2019，55（6）：1-5.

[5]? ? JUMAYEV S，MERDZAN M，BOYNOV K O，et al. The effect of PWM on rotor eddy-current losses in high-speed permanent magnet machines[J]. IEEE Transactions on Magnetics，2015，51（11）：1-4.

[6]? ? ZHOU P，LIN D，XIAO Y，et al. Temperature-dependent demagnetization model of permanent magnets for finite element analysis[J]. IEEE Transactions on Magnetics，2012，48（2）：1031-1034.

[7]? ? 朱衛光. 電動車輛永磁同步電機轉子永磁體渦流損耗及溫度場研究[D]. 北京：北京理工大學，2014.

[8]? ? 黃東洙，李偉力，王耀玉，等. 磁性槽楔對永磁電機轉子損耗及溫度場影響[J]. 電機與控制學報，2016，20（1）：60-66.

[9]? ? 司紀凱，張露鋒，封海潮，等. 異步起動表面-內置式永磁轉子同步電機特性分析及優化[J]. 電工技術學報，2017，32（9）：63-69.

[10]? WILLS D A，KAMPER M J. Reducing PM eddy current rotor losses by partial magnet and rotor yoke segmentation[C]//The XIX International Conference on Electrical Machines-ICEM 2010. September 6-8，2010，Rome，Italy. IEEE，2010：1-6.

[11]? 羅建華. 容錯式磁通切換永磁電機永磁體渦流損耗研究[D]. 鎮江：江蘇大學，2018.

[12]? EDE J D，ATALLAH K，JEWELL G W，et al. Effect of axial segmentation of permanent magnets on rotor loss in modular permanent-magnet brushless machines[J]. IEEE Transactions on Industry Applications，2007，43（5）：1207-1213.

[13]? CHEN C，WANG Y L，WEI J L，et al. Research on method for reducing eddy current loss of magnet in high-speed permanent magnet synchronous motor[C]//2019 22nd International Conference on Electrical Machines and Systems （ICEMS）. August 11-14，2019，Harbin，China. IEEE，2019：1-4.

[14]? SERGEANT P，VAN DEN BOSSCHE A. Segmentation of magnets to reduce losses in permanent-magnet synchronous machines[J]. IEEE Transactions on Magnetics，2008，44（11）：4409-4412.

[15]? 趙志剛，魏樂，溫濤，等. 變壓器鐵心模型諧波磁損耗的計算與分析[J]. 儀器儀表學報，2018，39（11）：19-27.

[16]? 吳斌，李強光，劉秀成，等. 基于磁特性曲線的桿件拉力測量及影響因素研究[J]. 儀器儀表學報，2015，36（3）：560-567.

[17]? LI Y J，JIANG B L，ZHANG C G，et al. Eddy-current loss measurement of permanent magnetic material at different frequency[J]. International Journal of Applied Electromagnetics and Mechanics，2019，61：S13-S22.