考慮諧波及集膚效應(yīng)的電工鋼片旋轉(zhuǎn)異常損耗計(jì)算與測(cè)量

李勁松，　楊慶新,2，　李永建，　張長(zhǎng)庚

(1.河北工業(yè)大學(xué) 電磁場(chǎng)與電器可靠性省部共建重點(diǎn)實(shí)驗(yàn)室，天津 300130；2.天津工業(yè)大學(xué) 電工電能新技術(shù)天津市重點(diǎn)實(shí)驗(yàn)室，天津 300387)

?

考慮諧波及集膚效應(yīng)的電工鋼片旋轉(zhuǎn)異常損耗計(jì)算與測(cè)量

李勁松1，楊慶新1,2，李永建1，張長(zhǎng)庚1

(1.河北工業(yè)大學(xué) 電磁場(chǎng)與電器可靠性省部共建重點(diǎn)實(shí)驗(yàn)室，天津 300130；2.天津工業(yè)大學(xué) 電工電能新技術(shù)天津市重點(diǎn)實(shí)驗(yàn)室，天津 300387)

摘要:為了準(zhǔn)確計(jì)算電工鋼片的旋轉(zhuǎn)異常損耗，根據(jù)諧波分析原理對(duì)建立的電工鋼疊片有限元模型進(jìn)行時(shí)步有限元仿真；基于旋轉(zhuǎn)鐵心損耗計(jì)算模型，通過(guò)考慮渦流集膚效應(yīng)對(duì)旋轉(zhuǎn)損耗系數(shù)的影響結(jié)合鋼片在中低頻率下的損耗測(cè)試獲得疊片損耗計(jì)算的關(guān)鍵系數(shù)，間接求得電工鋼片中旋轉(zhuǎn)異常損耗的計(jì)算式。利用構(gòu)建的新型三維磁特性測(cè)試系統(tǒng)對(duì)典型電工鋼疊片樣品進(jìn)行橢圓形旋轉(zhuǎn)與交變勵(lì)磁方式下的寬頻鐵耗實(shí)驗(yàn)測(cè)量，并定量地進(jìn)行了對(duì)比與分析。結(jié)果表明：2種勵(lì)磁方式下疊片損耗的變化規(guī)律相類似，但其橢圓形旋轉(zhuǎn)各損耗都要比交變時(shí)的對(duì)應(yīng)損耗大，必須認(rèn)真考慮諧波、集膚效應(yīng)和旋轉(zhuǎn)勵(lì)磁等對(duì)材料特性的影響；所計(jì)算出的旋轉(zhuǎn)異常損耗也是相對(duì)較小，在1 kHz時(shí)也未占到旋轉(zhuǎn)總鐵耗的5%。從而驗(yàn)證了所推導(dǎo)出的計(jì)算式和測(cè)量手段的正確性與可行性。

關(guān)鍵詞:電工鋼疊片；諧波分析；集膚效應(yīng)；旋轉(zhuǎn)異常損耗；三維磁特性

0引言

電氣設(shè)備中的諧波主要包括由于非正弦電源供電而使其繞組中通入諧波電流所產(chǎn)生的時(shí)間諧波和即使是通入三相正弦波電流，由于設(shè)備氣隙磁勢(shì)在空間的分布非正弦而帶來(lái)的一系列空間高次諧波[1]。處在旋轉(zhuǎn)磁場(chǎng)中的電工鋼疊片內(nèi)存在渦流，使其上部的電密增大，下部減小，從而引起電密的分布不均，這種集中在電工鋼疊片表面的現(xiàn)象就是渦流集膚效應(yīng)[2]。

諧波對(duì)電力系統(tǒng)的影響和危害是十分嚴(yán)重的，尤其是對(duì)運(yùn)行頻率在幾百甚至幾千赫茲的高頻電機(jī)和變壓器而言，由于諧波的存在和顯著的渦流集膚效應(yīng)，特別是在旋轉(zhuǎn)勵(lì)磁方式下，使得其由電工鋼片疊置而成的鐵心損耗急劇增加以致過(guò)熱燒毀，從而降低了電氣設(shè)備的使用效率和壽命[3-4]。目前，對(duì)電機(jī)和變壓器的鐵耗計(jì)算大都是建立在工頻正弦交變磁場(chǎng)的假設(shè)基礎(chǔ)之上，但在非正弦供電時(shí)鐵耗的計(jì)算模型發(fā)生了變化[5]，且對(duì)考慮諧波及集膚效應(yīng)的電工鋼片旋轉(zhuǎn)異常損耗的獨(dú)立研究，限于其建模與計(jì)算較為困難和不易測(cè)量與分離等原因，卻很少有文獻(xiàn)深入涉及。

文獻(xiàn)[6]以諧波分析法為基礎(chǔ)，提出了一種在非正弦供電情況下綜合考慮電源諧波、電機(jī)鐵心飽和以及集膚效應(yīng)時(shí)計(jì)算高速變頻電機(jī)鐵耗的簡(jiǎn)單而有效的方法，卻由于比重太小而忽略了異常損耗且研究是在交變勵(lì)磁方式下。文獻(xiàn)[7]在考慮諧波影響時(shí)，通過(guò)對(duì)模型進(jìn)行仿真，分析和評(píng)價(jià)了變壓器在非線性負(fù)載下的損耗分布，從而驗(yàn)證了由于諧波失真將導(dǎo)致變壓器不必要的繞組損耗和典型的溫度上升，降低電能質(zhì)量。但并沒(méi)有對(duì)變壓器損耗分布進(jìn)行實(shí)驗(yàn)測(cè)量，尤其是疊片的鐵耗，更不用說(shuō)異常損耗了。文獻(xiàn)[8]對(duì)電工鋼片和旋轉(zhuǎn)電機(jī)中的旋轉(zhuǎn)鐵耗進(jìn)行了建模與測(cè)量，提出了在電工鋼片中利用圓形和橢圓形旋轉(zhuǎn)磁密矢量計(jì)算旋轉(zhuǎn)磁滯、渦流和異常損耗的新方法。根據(jù)旋轉(zhuǎn)磁場(chǎng)的分布，將這些公式與有限元法結(jié)合使用來(lái)計(jì)算旋轉(zhuǎn)磁通電機(jī)中的鐵耗，理論與實(shí)驗(yàn)結(jié)果都達(dá)到了很好的匹配。結(jié)果表明，在同樣頻率和磁密的條件下，旋轉(zhuǎn)勵(lì)磁方式下的各損耗都要比交變時(shí)大。然而，即使是在旋轉(zhuǎn)勵(lì)磁頻率相對(duì)較高時(shí)，這些計(jì)算與測(cè)量也均未考慮諧波及集膚效應(yīng)的影響。

誠(chéng)然，上述對(duì)電工鋼片交變和旋轉(zhuǎn)鐵耗的計(jì)算與測(cè)量具有重要的工程意義。考慮到目前對(duì)旋轉(zhuǎn)勵(lì)磁的認(rèn)識(shí)還有待深入，且限于產(chǎn)生異常損耗的物理機(jī)理過(guò)程非常復(fù)雜，其在總鐵耗中所占的比重相對(duì)較小，多數(shù)研究者都是通過(guò)統(tǒng)計(jì)學(xué)方法給出數(shù)學(xué)表達(dá)式等原因，本文基于旋轉(zhuǎn)鐵耗計(jì)算模型，綜合考慮旋轉(zhuǎn)磁場(chǎng)諧波及渦流集膚效應(yīng)的影響，對(duì)電工鋼片中的旋轉(zhuǎn)異常損耗進(jìn)行間接計(jì)算與測(cè)量。

1旋轉(zhuǎn)磁場(chǎng)的諧波影響與分析

當(dāng)施加旋轉(zhuǎn)勵(lì)磁時(shí)，無(wú)論是晶粒取向還是非晶粒取向電工鋼片，其飽和磁化后磁場(chǎng)強(qiáng)度H波形的5次諧波含量與不飽和時(shí)相比，都會(huì)有顯著增加[9]，因此考慮旋轉(zhuǎn)磁場(chǎng)中的諧波是非常有必要的。

根據(jù)諧波分析原理，在圖1所示的電工鋼疊片有限元模型中任意一點(diǎn)的磁密波形都可以分解為一系列的諧波分量，即任意磁場(chǎng)波形勵(lì)磁時(shí)產(chǎn)生的損耗都等于其基波與各次諧波分量產(chǎn)生的損耗之和[10]。可設(shè)

(1)

式中：B1和fr分別為基波的旋轉(zhuǎn)峰值磁密和勵(lì)磁頻率；Bn和φn分別為n次諧波分量的旋轉(zhuǎn)峰值磁密和初始相角。

對(duì)于含有n次諧波的橢圓形旋轉(zhuǎn)磁場(chǎng)，可以將其分解為長(zhǎng)軸磁密為Bnx和短軸磁密為Bny的2個(gè)正交的交變磁場(chǎng)。圖2和3所示分別為經(jīng)過(guò)時(shí)步有限元仿真計(jì)算所得圖1中A點(diǎn)的磁密n次諧波矢量形成的橢圓形磁場(chǎng)及各次諧波分量峰值圖。

圖1　電工鋼疊片有限元模型Fig.1　Finite element model of electrical steel sheets

圖2　　　　電工鋼疊片有限元模型中A點(diǎn)的磁密n次　　　諧波矢量形成的橢圓形磁場(chǎng)Fig.2　　　　Elliptical nthharmonic flux density vector　　　 waveforms at position A in the finite　　　element model of electrical steel sheets

2間接計(jì)算式建立

根據(jù)旋轉(zhuǎn)鐵心損耗計(jì)算模型[11]，電工鋼片在旋轉(zhuǎn)勵(lì)磁方式下的總鐵耗Prot可寫(xiě)為

Prot=Phr+Per+Par，W/kg。

(2)

式中：Phr、Per和Par分別代表旋轉(zhuǎn)磁滯、渦流和異常損耗。

2.1計(jì)算式推導(dǎo)

據(jù)前述可知，在計(jì)算橢圓形旋轉(zhuǎn)磁場(chǎng)勵(lì)磁產(chǎn)生的損耗時(shí)，可將其等效為2個(gè)正交的交變損耗。則

電工鋼片單位質(zhì)量的旋轉(zhuǎn)磁滯損耗為

(3)

式中：khr和α同為旋轉(zhuǎn)磁滯損耗系數(shù)。

單位質(zhì)量的旋轉(zhuǎn)渦流損耗為

(4)

式中：ker為旋轉(zhuǎn)渦流損耗系數(shù)。

圖3　　　　電工鋼疊片有限元模型中A點(diǎn)的磁密　　　各次諧波分量峰值圖Fig.3　　　　Peak harmonic flux densities at position A in the　　　 finite element model of electrical steel sheets

采用時(shí)步有限元法對(duì)施加橢圓形旋轉(zhuǎn)勵(lì)磁的電工鋼疊片進(jìn)行仿真計(jì)算，可以得到疊片中各剖分計(jì)算單元的磁密矢量B的波形，則其旋轉(zhuǎn)總鐵耗等于各剖分計(jì)算單元的損耗之和[12]。即

(5)

式中：Ni為電工鋼疊片有限元模型中的剖分計(jì)算單元數(shù)；pi為第i個(gè)剖分計(jì)算單元的單位質(zhì)量鐵耗。

對(duì)式(5)經(jīng)由分析可知：

1)就本電工鋼疊片有限元模型，應(yīng)用式(5)通過(guò)時(shí)步有限元法可計(jì)算考慮15次及以下諧波分量的疊片旋轉(zhuǎn)總鐵耗。

2)所得到的是有限元仿真數(shù)值解，需要編寫(xiě)計(jì)算程序來(lái)實(shí)現(xiàn)。

為便于后續(xù)的實(shí)驗(yàn)測(cè)量，應(yīng)用Poynting定理[13-14]推導(dǎo)出的電工鋼片單位質(zhì)量的二維平均旋轉(zhuǎn)總鐵耗，將式(5)改寫(xiě)為

(6)

式中：Tr為旋轉(zhuǎn)勵(lì)磁周期；ρ為電工鋼片的質(zhì)量密度；Bx(t)和By(t)分別為鋼片中的法向和切向旋轉(zhuǎn)磁密；Hx(t)和Hy(t)分別為法向和切向磁場(chǎng)強(qiáng)度。

對(duì)式(6)經(jīng)由分析可知：當(dāng)施加橢圓形旋轉(zhuǎn)勵(lì)磁時(shí)，通過(guò)測(cè)量電工鋼疊片表面的磁場(chǎng)強(qiáng)度H和內(nèi)部的磁密B隨時(shí)間的變化數(shù)據(jù)，就可計(jì)算出包括旋轉(zhuǎn)磁滯、渦流和異常損耗在內(nèi)的旋轉(zhuǎn)總鐵耗。

聯(lián)立式(2)～式(4)和式(6)，可得電工鋼片在考慮諧波影響時(shí)的旋轉(zhuǎn)異常損耗間接計(jì)算式如下

Par=Prot-(Phr+Per)=

(7)

2.2影響旋轉(zhuǎn)損耗系數(shù)的因素分析與確定

2.2.1旋轉(zhuǎn)勵(lì)磁頻率的變化影響與分析

利用課題組構(gòu)建的新型三維磁特性測(cè)試系統(tǒng)(由測(cè)試儀、功率放大器及信號(hào)采集和數(shù)字信號(hào)處理系統(tǒng)等組成)，在對(duì)牌號(hào)為50WW800、單片厚度為0.50 mm且?guī)в蠦-H傳感線圈(如圖4所示)的冷軋非晶粒取向電工鋼疊片樣品進(jìn)行橢圓形旋轉(zhuǎn)勵(lì)磁實(shí)驗(yàn)時(shí)，將磁場(chǎng)分別按順時(shí)針和逆時(shí)針?lè)较蛐D(zhuǎn)，由于疊片樣品存在磁各向異性，因此所得H軌跡并不完全相同，如圖5所示。

圖4　帶有B-H傳感線圈的50WW800疊片樣品Fig.4　　　　50WW800 laminated specimen with B-H　　　sensing coil

疊片樣品在中低旋轉(zhuǎn)勵(lì)磁頻率(5～100 Hz)下的橢圓形旋轉(zhuǎn)總鐵耗測(cè)試數(shù)據(jù)[15-16]，如圖6所示。

圖5　50 Hz時(shí)疊片樣品在不同旋轉(zhuǎn)方向下的H軌跡Fig.5　　　　H-loci of laminated specimen under different　　　rotating directions at 50 Hz

圖6　　　　中低旋轉(zhuǎn)勵(lì)磁頻率fr下50WW800疊片　　　樣品的橢圓形旋轉(zhuǎn)總鐵耗ProtFig.6　　　　Total elliptical rotation core losses Protof　　　50WW800 laminated specimen at low to　　　medium frequencies fr

圖7所示為經(jīng)過(guò)數(shù)學(xué)擬合后的旋轉(zhuǎn)損耗系數(shù)隨旋轉(zhuǎn)勵(lì)磁頻率的變化關(guān)系曲線。

由圖7可以看出：旋轉(zhuǎn)磁滯損耗系數(shù)khr和α隨著旋轉(zhuǎn)勵(lì)磁頻率fr的增加變化不是很大，但旋轉(zhuǎn)渦流損耗系數(shù)ker卻隨著旋轉(zhuǎn)勵(lì)磁頻率fr的增加出現(xiàn)了減小的趨勢(shì)。經(jīng)由分析可知：同交變勵(lì)磁損耗特性相一致，當(dāng)橢圓形旋轉(zhuǎn)勵(lì)磁頻率較低時(shí)，采用旋轉(zhuǎn)鐵耗計(jì)算模型所得電工鋼片的損耗誤差不會(huì)很大，但當(dāng)旋轉(zhuǎn)勵(lì)磁頻率較高，尤其是達(dá)到幾百甚至上千赫茲時(shí)，將導(dǎo)致鋼片中的渦電密分布不均，此時(shí)鋼片的旋轉(zhuǎn)渦流損耗系數(shù)與旋轉(zhuǎn)勵(lì)磁頻率也有關(guān)系，必須要考慮渦流集膚效應(yīng)的影響。

圖7　　　　旋轉(zhuǎn)損耗系數(shù)隨旋轉(zhuǎn)勵(lì)磁頻率的變化　　　關(guān)系曲線(數(shù)學(xué)擬合)Fig.7　　　　Rotating loss coefficients vs. rotating excitation　　　 frequency(Mathematical Fitting)

2.2.2集膚效應(yīng)的影響與分析

集膚效應(yīng)將影響電工鋼疊片中的渦電場(chǎng)分布，尤其是當(dāng)旋轉(zhuǎn)勵(lì)磁頻率很高時(shí)，渦電流幾乎集中在疊片的表面流動(dòng)，這種渦電密的數(shù)值由表及里逐漸減小的特性通常用集膚(透入)深度δe來(lái)表示[17]

(8)

式中：ω為旋轉(zhuǎn)勵(lì)磁角頻率；σ和μ分別為電工鋼片的電導(dǎo)率和磁導(dǎo)率。

通過(guò)對(duì)圖1的電工鋼疊片有限元模型進(jìn)行時(shí)步有限元仿真計(jì)算，所得集膚(透入)深度為δe≈1.0×10-3m，這一數(shù)值遠(yuǎn)遠(yuǎn)小于模型中鋼片的寬度Yw，即δe<

圖8　　　　沿y軸在z=0處的磁場(chǎng)強(qiáng)度Hy　　　和磁密By的峰值Fig.8　　　　Peak values of magnetic strength Hyand　　　flux density Byalong y-axis at z=0

與之相反，模型中邊界區(qū)域的單層鋼片的厚度d與集膚(透入)深度δe可比擬，即δe～d。由于渦電場(chǎng)的存在而產(chǎn)生感應(yīng)磁場(chǎng)，從而改變內(nèi)部旋轉(zhuǎn)磁化場(chǎng)的分布，表現(xiàn)出與外施勵(lì)磁場(chǎng)更加密切的相關(guān)性，如圖9所示。

圖9　　　　沿z軸在y=0處的磁場(chǎng)強(qiáng)度Hz　　　和磁密Bz的峰值Fig.9　　　　Peak values of magnetic strength Hzand　　　flux density Bzalong z-axis at y=0

結(jié)合式(8)分析圖8和圖9可知：隨著旋轉(zhuǎn)勵(lì)磁頻率的增加，集膚(透入)深度將變小而趨于表面，導(dǎo)致電工鋼疊片厚度上渦流的分布不均，對(duì)其內(nèi)部旋轉(zhuǎn)磁化場(chǎng)分布的影響也會(huì)加深。所以在確定旋轉(zhuǎn)損耗系數(shù)特別是旋轉(zhuǎn)渦流損耗系數(shù)時(shí)，必須要考慮渦流集膚效應(yīng)的影響，以減小損耗計(jì)算與測(cè)量的誤差。

2.2.3旋轉(zhuǎn)損耗系數(shù)的確定

據(jù)式(7)可知，如果想得到電工鋼片在考慮諧波及集膚效應(yīng)時(shí)的旋轉(zhuǎn)異常損耗，必須要先獲得各項(xiàng)損耗的關(guān)鍵系數(shù)khr、α和ker，具體措施如下[18]

1)對(duì)圖6中50WW800疊片樣品在中低旋轉(zhuǎn)勵(lì)磁頻率下的橢圓形旋轉(zhuǎn)總鐵耗實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合與線性回歸分析，使其滿足式(9)的約束條件且使ξ達(dá)到最小。

(9)

2)在計(jì)算旋轉(zhuǎn)渦流損耗系數(shù)時(shí)，考慮到旋轉(zhuǎn)勵(lì)磁頻率的變化及渦流集膚效應(yīng)的影響，在中頻段(50～100 Hz)選取旋轉(zhuǎn)損耗實(shí)驗(yàn)數(shù)據(jù)時(shí)進(jìn)行了取平均與歸一化處理。所得50WW800電工鋼片的旋轉(zhuǎn)損耗系數(shù)如表1所示。

表1　50WW800電工鋼片旋轉(zhuǎn)損耗系數(shù)

因此式(7)可進(jìn)一步寫(xiě)為

Par=Prot-(Phr+Per)=

(10)

式(10)即為考慮諧波及集膚效應(yīng)的電工鋼片旋轉(zhuǎn)異常損耗間接計(jì)算式，它是經(jīng)過(guò)中低頻率下的鋼片損耗測(cè)試結(jié)合時(shí)步有限元法仿真計(jì)算得到的。

3實(shí)驗(yàn)測(cè)量

3.1測(cè)量研究

為了驗(yàn)證所推導(dǎo)出的旋轉(zhuǎn)異常損耗間接計(jì)算式的正確性，依舊利用構(gòu)建的新型三維磁特性測(cè)試系統(tǒng)，選用圖4的50WW800電工鋼疊片樣品進(jìn)行橢圓形旋轉(zhuǎn)勵(lì)磁方式下的寬頻(100～1 000 Hz)鐵耗實(shí)驗(yàn)測(cè)量，如圖10所示。圖11所示為特定旋轉(zhuǎn)峰值磁密時(shí)的總鐵耗分離曲線。(本文的旋轉(zhuǎn)磁密波形B(t)計(jì)算到7次諧波分量，即n=7)

圖10　　　　較寬旋轉(zhuǎn)勵(lì)磁頻率fr下50WW800疊片　　　樣品的橢圓形旋轉(zhuǎn)總鐵耗ProtFig.10　　　　Total elliptical rotation core losses Protof　　　50WW800 laminated specimen under wide　　　 rotating excitation frequencies fr

圖11　　　　Bmr=1.0 T時(shí)疊片樣品的橢圓形　　　旋轉(zhuǎn)總鐵耗分離曲線Fig.11　　　　Total elliptical rotation core losses separation　　　 curves of laminated specimen at Bmr=1.0 T

由圖10和圖11可以看出，電工鋼疊片在旋轉(zhuǎn)勵(lì)磁方式下的損耗變化規(guī)律為

1)在達(dá)到飽和磁化的過(guò)程中，其橢圓形旋轉(zhuǎn)總鐵耗隨著旋轉(zhuǎn)勵(lì)磁頻率的增加而增大，且頻率越高損耗越大。

2)當(dāng)旋轉(zhuǎn)勵(lì)磁頻率較低時(shí)，旋轉(zhuǎn)磁滯損耗是旋轉(zhuǎn)總鐵耗的主要部分；隨著頻率的增加，旋轉(zhuǎn)渦流損耗在旋轉(zhuǎn)總鐵耗中所占的比重越來(lái)越大，當(dāng)旋轉(zhuǎn)勵(lì)磁頻率高達(dá)1 kHz時(shí)，實(shí)測(cè)疊片樣品的旋轉(zhuǎn)渦流損耗幾乎占據(jù)了旋轉(zhuǎn)總鐵耗的80%。

3)旋轉(zhuǎn)異常損耗也是相對(duì)較小，在1 kHz時(shí)也未占到旋轉(zhuǎn)總鐵耗的5%。

3.2對(duì)比與分析

同樣是在交變峰值磁密Bma=1.0 T時(shí)，再次對(duì)同一電工鋼疊片樣品進(jìn)行交變總鐵耗的寬頻實(shí)驗(yàn)測(cè)量，并定量比較電工鋼片在不同勵(lì)磁方式下的損耗特性，如圖12所示。

圖12　　　　Bm=1.0 T時(shí)疊片樣品橢圓形旋轉(zhuǎn)與　　　交變損耗結(jié)果比較Fig.12　　　　Elliptical rotation vs. alternating losses results　　　of laminated specimen at Bm=1.0 T

由圖12不難發(fā)現(xiàn)：電工鋼疊片樣品在橢圓形旋轉(zhuǎn)勵(lì)磁方式下的總鐵耗、渦流損耗、磁滯損耗和異常損耗都要比其在交變勵(lì)磁時(shí)的對(duì)應(yīng)損耗大。這是因?yàn)闄E圓形旋轉(zhuǎn)磁場(chǎng)終究是由2個(gè)交變磁場(chǎng)共同作用形成的。由于電工鋼片在不同勵(lì)磁方式下產(chǎn)生損耗的物理機(jī)理并不完全一樣，尤其是磁滯損耗，所以也不能簡(jiǎn)單地認(rèn)為旋轉(zhuǎn)損耗就等于交變損耗的2倍，它們之間的數(shù)量關(guān)系因材料、勵(lì)磁頻率和峰值磁密等的不同而不同。

綜上所述，在對(duì)電工鋼片損耗進(jìn)行計(jì)算時(shí)，要認(rèn)真考慮諧波、集膚效應(yīng)和旋轉(zhuǎn)勵(lì)磁等對(duì)材料特性的影響，特別是在高頻高磁密時(shí)必須研究損耗系數(shù)的變化規(guī)律，套用經(jīng)典損耗分離模型的方法必將受到限制。

4結(jié)論

本文基于旋轉(zhuǎn)鐵耗計(jì)算模型，根據(jù)諧波分析原理，在考慮渦流集膚效應(yīng)對(duì)旋轉(zhuǎn)損耗系數(shù)的影響時(shí)，間接求得電工鋼片中的旋轉(zhuǎn)異常損耗計(jì)算式。利用構(gòu)建的新型三維磁特性測(cè)試系統(tǒng)對(duì)典型電工鋼疊片樣品進(jìn)行了橢圓形旋轉(zhuǎn)和交變勵(lì)磁方式下的損耗實(shí)驗(yàn)測(cè)量，并定量地進(jìn)行了對(duì)比與分析，驗(yàn)證了所推導(dǎo)出的旋轉(zhuǎn)異常損耗間接計(jì)算式和測(cè)量手段的正確性與可行性。得到了如下結(jié)論：

1)2種不同勵(lì)磁方式下電工鋼片的損耗變化規(guī)律相類似，但橢圓形旋轉(zhuǎn)各損耗都要比交變勵(lì)磁時(shí)的對(duì)應(yīng)損耗大，必須認(rèn)真考慮諧波、集膚效應(yīng)和旋轉(zhuǎn)勵(lì)磁等對(duì)材料特性的影響。

2)由間接計(jì)算式所得旋轉(zhuǎn)異常損耗也是相對(duì)較小，在旋轉(zhuǎn)勵(lì)磁頻率高達(dá)1k Hz時(shí)也未占到旋轉(zhuǎn)總鐵耗的5%。

本文考慮諧波及集膚效應(yīng)的電工鋼片旋轉(zhuǎn)異常損耗間接計(jì)算式研究有助于解決電機(jī)和變壓器等電氣設(shè)備的可能的危險(xiǎn)的局部過(guò)熱問(wèn)題，延長(zhǎng)使用壽命，同時(shí)可促進(jìn)工程中電磁裝備的優(yōu)化與設(shè)計(jì)，具有一定的理論與實(shí)用價(jià)值。

參 考 文 獻(xiàn):

[1]章瑋, 陳伯建, 張平. 雙三相永磁同步電機(jī)諧波電流抑制技術(shù)[J]. 電機(jī)與控制學(xué)報(bào), 2015, 19(1): 23-28.

ZHANG Wei, CHEN Bojian, ZHANG Ping. Techniques to restrain harmonics of dual three-phase permanent magnet synchronous motor[J]. Electric Machines and Control, 2015, 19(1): 23-28.

[2]ZHANG Y.Experimental, theoretical and numerical investigation of core losses in laminated magnetic materials[D]. Potsdam, NY, USA: Clarkson University, 2008.

[3]INDARACK P, DOUANGSYLA S,JOOCHIM C, et al. A harmonic loss calculation of PWM-fed induction motors using loss factor characteristics[C]//2004 IEEE Region 10 Conference (TENCON 2004), Nov. 2004, 3: 236- 239.

[4]黃東洙, 李偉力, 王耀玉, 等. 磁性槽楔對(duì)永磁電機(jī)轉(zhuǎn)子損耗及溫度場(chǎng)影響[J]. 電機(jī)與控制學(xué)報(bào), 2016, 20(1): 60-66.

HUANG Dongzhu, LI Weili, WANG Yaoyu, et al. Influence of magnetic slot wedge on rotor losses and temperature field of PMSM[J]. Electric Machines and Control, 2016, 20(1): 60-66.

[5]張洪亮, 鄒繼斌. 考慮旋轉(zhuǎn)磁通的PMSM鐵心損耗數(shù)值計(jì)算[J]. 電機(jī)與控制學(xué)報(bào), 2007, 11(4): 340-348.

ZHANG Hongliang, ZOU Jibin. Numerical calculation and analysis of stator iron losses due to rotational flux and magnet pole eddy-current loss in PMSM[J]. Electric Machines and Control, 2007, 11(4): 340-348.

[6]方瑞明, 王榕. 基于諧波分析法的高速變頻電機(jī)鐵耗計(jì)算方法[J]. 電機(jī)與控制學(xué)報(bào), 2004, 8(1): 25-28.

FANG Ruiming, WANG Rong. Calculation of iron loss in high-speed induction motor fed by nonsinusoidal voltage source based on harmonic analysis method [J]. Electric Machines and Control, 2004, 8(1): 25-28.

[7]SAIDD M,NOR K M，MAJID M S. Analysis of distribution transformer losses and life expectancy using measured harmonic data[J]. Journal of Nematology, 2014, 46(2): 1-6.

[8]黃允凱. 鐵磁材料損耗及高速軟磁復(fù)合材料電機(jī)的研究[D]. 南京: 東南大學(xué), 2007.

[9]黃平林. 旋轉(zhuǎn)電機(jī)鐵心損耗的分析與計(jì)算[D]. 南京: 東南大學(xué), 2007.

[10]KINNARES V, JMANCHAI P, SUKSAWAT D, et al. Effect of motor parameter changes on harmonic power [C]//Proceedings of the IEEE 1999 International Conference on Power Electronics and Drive Systems (PEDS’99), July，1999, 2: 1061-1066.

[11]ZHU J G, RAMSDEN V S. Improved formulations for rotational core losses in rotating electrical machines[J]. IEEE Transactions on Magnetics, 1998, 34(4): 2234-2242.

[12]江善林, 鄒繼斌, 徐永向, 等. 考慮旋轉(zhuǎn)磁通和趨膚效應(yīng)的變系數(shù)鐵耗計(jì)算模型[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2011, 31(3): 104-110.

JIANG Shanlin, ZOU Jibin, XU Yongxiang, et al. Variable coefficient iron loss calculating model considering rotational flux and skin effect [J]. Proceedings of the CSEE, 2011, 31(3): 104-110.

[13]ZHU J G. Numerical modelling of magnetic materials for computer aided design of electromagnetic devices[D]. Sydney: University of Technology, Sydney, 1994.

[14]戴興建, 鄧占峰, 劉剛, 等. 大容量先進(jìn)飛輪儲(chǔ)能電源技術(shù)發(fā)展?fàn)顩r[J]. 電工技術(shù)學(xué)報(bào), 2011, 26(7): 133-140.

DAI Xingjian, DENG Zhanfeng，LIU Gang, et al. Review on advanced flywheel energy storage system with large scale[J]. Transactions of China Electrotechnical Society, 2011, 26(7): 133-140.

[15]李勁松, 楊慶新, 李永建, 等. 高頻高磁密時(shí)疊置硅鋼片的鐵芯損耗計(jì)算式改進(jìn)[J]. 高電壓技術(shù), 2016, 2(3): 994-1002.

LI Jingsong, YANG Qingxin, LI Yongjian, et al. Improved calculation formula of core losses for laminated silicon steels at high frequency and high flux density[J]. High Voltage Engineering, 2016, 42(3): 994-1002.

[16]LI J S, YANG Q X, LI Y J, et al. Anomalous loss modeling and validation of magnetic materials in electrical engineering[J]. IEEE Transactions on Applied Superconductivity, 2016, 26(4): 8800105.

[17]MTHOMBENI T L. Improved lamination core loss measurements and calculations[D]. Potsdam, NY, USA: Clarkson University, 2006.

[18]DLALAE．A simplified iron loss model for laminated magnetic cores[J]. IEEE Transactions on Magnetics, 2008, 44 (11): 3169-3172.

(編輯：劉素菊)

Measurement and calculation of rotating anomalous loss inelectrical steel sheets considering harmonics and skin effect

LI Jing-song1,YANG Qing-xin1, 2,LI Yong-jian1,ZHANG Chang-geng1

(1.Province-Ministry Joint Key Laboratory of EFEAR, Hebei University of Technology, Tianjin 300130, China;2.Tianjin Key Laboratory of AEEET, Tianjin Polytechnic University, Tianjin 300387, China)

Abstract:In order to accurately calculate the rotating anomalous loss in electrical steel sheets,sheets finite element model was simulated by adopting time-stepping finite element method based on Fourier analysis. Then, the key loss coefficients of the sheets were obtained by experimental measurement at low to medium frequencies with considering the influence of eddy current skin effect to the loss coefficients on the basis of rotating core losses calculation model. So, the calculation formula of rotating anomalous loss in electrical steel sheets was deducd indirectly. By using a novel 3-D magnetic properties testing system, the total broadband core losses under elliptical rotation and alternating excitations were carried out for a typical electrical steel laminated specimen. Meanwhile, losses characteristics of the two excitation modes were compared and analyzed quantitatively. The results show that both the losses variations are analogous, but varieties of elliptical rotation excitation losses are bigger than the corresponding alternating losses, so it is necessary to consider the influences of harmonic, skin effect and rotating excitation, etc, to the materials properties seriously. The calculation of rotating anomalous loss is relatively small and is not accounted for 5% of the total rotating core losses up to 1k Hz. Thus, it validates the correctness and feasibility of the calculation formula and measurement means.

Keywords:electrical steel sheets; Fourier analysis; skin effect; rotating anomalous loss; 3-D magnetic properties

收稿日期:2014-11-03

基金項(xiàng)目：國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(51237005)；河北省自然科學(xué)基金(青年基金)(E2014202137)

作者簡(jiǎn)介：李勁松(1987—)，男，博士研究生，研究方向?yàn)殡姽ご判圆牧洗艙p耗特性模擬、測(cè)量分析與應(yīng)用； 楊慶新(1961—)，男，博士，教授，博士生導(dǎo)師，研究方向?yàn)楣こ屉姶艌?chǎng)數(shù)值分析與磁技術(shù)相關(guān)工作； 李永建(1978—)，男，博士，教授，碩士生導(dǎo)師，研究方向?yàn)殡姽ご判圆牧洗盘匦詼y(cè)量和磁滯建模； 張長(zhǎng)庚(1985—)，男，博士研究生，研究方向?yàn)殡姽ご判圆牧洗盘匦詼y(cè)量和磁滯建模。

通訊作者:李勁松

DOI:10.15938/j.emc.2016.06.006

中圖分類號(hào):TM 271

文獻(xiàn)標(biāo)志碼:A

文章編號(hào):1007-449X(2016)06-0042-08