電工軟磁材料旋轉(zhuǎn)磁滯損耗測量及建模

電工軟磁材料旋轉(zhuǎn)磁滯損耗測量及建模張長庚;楊慶新;李永建【摘要】電力變壓器和電機(jī)中存在旋轉(zhuǎn)鐵心損耗是其損耗預(yù)測不準(zhǔn)的主要原因之一.針對該問題，提出電工軟磁材料旋轉(zhuǎn)磁滯損耗測量及建模方法.首先，矢量磁滯損耗分解為切向損耗和法向損耗兩部分，分別根據(jù)圓形旋轉(zhuǎn)損耗和交變損耗建模.其次，對軟磁復(fù)合材料和無取向電工鋼片進(jìn)行模型參數(shù)辨識，并比較兩種材料的損耗特性和模型參數(shù).最后，利用三維磁特性測量裝置進(jìn)行多種勵磁模式下的損耗測量，并對比實(shí)驗(yàn)結(jié)果與模型預(yù)測值結(jié)果表明，在旋轉(zhuǎn)復(fù)雜激勵下，所提出的模型比傳統(tǒng)的Steinmetz模型有更高的精度.％Therotationalfluxinthecoreofpowertransformerandmotorcausestheinaccuracyefficiencycomputationofmagneticapparatus.Inordertoaddresstheproblem,thisworkfocusesonmagneticmeasurementofrotationalmagnetizationanditshysteresislossmodeling.Firstly,thevectorhysteresislossisdecomposedintothetransverseandradialcomponentswhicharedefinedbasedonrotationallossandalternatingloss,respectively.Secondly,themodelparametersofsoftmagneticcompositesandnon-orientedsiliconsteelareidentified.Theirmagneticpropertiesandmodelparametersarecomparedandanalyzed.Finally,aseriesofcomplexexcitationmagnetizationexperimentsonnon-orientedsiliconsteelarecarriedout.Themeasurementresultsarecomparedwiththepredictedvaluesbythemodel.ItshowsthattheproposedmodelgivesmoreexactpredictedresultsthanthetraditionalSteinmetzmethod.【期刊名稱】《電工技術(shù)學(xué)報》【年(卷)，期】2017(032)011【總頁數(shù)】9頁(P208-216)【關(guān)鍵詞】電機(jī)鐵耗;磁滯損耗;電磁場數(shù)值模擬;旋轉(zhuǎn)損耗【作者】張長庚;楊慶新;李永建【作者單位】河北工業(yè)大學(xué)電磁場與電氣可靠性省部共建重點(diǎn)實(shí)驗(yàn)室天津300130;河北工業(yè)大學(xué)電磁場與電氣可靠性省部共建重點(diǎn)實(shí)驗(yàn)室天津300130沃津工業(yè)大學(xué)電工電能新技術(shù)天津市重點(diǎn)實(shí)驗(yàn)室天津300387;河北工業(yè)大學(xué)電磁場與電氣可靠性省部共建重點(diǎn)實(shí)驗(yàn)室天津300130【正文語種】中文【中圖分類】TM15高效率的電機(jī)和電力變壓器的使用能夠節(jié)約社會總能源消耗，是響應(yīng)國家節(jié)能減排發(fā)展戰(zhàn)略的有效措施。設(shè)計高效率的電機(jī)和電力變壓器，一方面改進(jìn)或提出新的結(jié)構(gòu)、工藝和控制策略，例如三維磁通電機(jī)；另一方面減少磁性材料本身的鐵耗，例如采用新型非晶納米晶材料［1］。兩方面都依賴準(zhǔn)確的鐵心材料的磁特性測量和建模早在電機(jī)等電磁裝置用于工業(yè)生產(chǎn)之初，1892年S.Steinmetz提出了著名的磁滯損耗計算Steinmetz公式［2］。該公式適用于正弦穩(wěn)態(tài)下的損耗預(yù)測，隨著HVDC和電力電子設(shè)備的廣泛應(yīng)用，大量諧波引入電機(jī)和變壓器設(shè)備，為解決復(fù)雜磁通密度波形的激勵，提出了多種修正的Steinmetz公式。文獻(xiàn)［3］引入等價頻率的概念，考慮非正弦磁通密度勵磁下的損耗改進(jìn)Steinmetz公式。文獻(xiàn)［4］提出廣義的Steinmetz公式，進(jìn)一步考慮小磁滯回環(huán)對損耗的貢獻(xiàn)。文獻(xiàn)［5］考慮頻率和磁通密度對Steinmetz公式系數(shù)的影響，通過擬合變系數(shù)的公式，提高損耗計算精度。文獻(xiàn)［6］考慮電力變壓器的雜散損耗和屏蔽鋼板的鐵耗。為了獲得更高精度磁滯損耗的估計，各種經(jīng)典標(biāo)量磁滯模型，如Preisach模型、Play模型和Jiles-Atherton(JA)模型等相繼被提出。在電力變壓器和電機(jī)實(shí)際運(yùn)行中發(fā)現(xiàn)旋轉(zhuǎn)磁場激勵下?lián)p耗大于單方向的交變損耗。對旋轉(zhuǎn)損耗的研究分為兩類：①基于Steinmetz公式方法和矢量磁滯模型方法，文獻(xiàn)［7］在頻域空間把二維空間磁通密度B軌跡分解為圓形旋轉(zhuǎn)磁通密度和交變磁通密度兩部分，根據(jù)交變損耗系數(shù)和旋轉(zhuǎn)損耗系數(shù)分別利用Steinmetz公式進(jìn)行損耗估計；②建立矢量的磁滯模型，如新型混合矢量磁滯模型。文獻(xiàn)［8］利用Stoner-Wohlfarth模型求解磁滯算子矢量磁化強(qiáng)度方向，再利用Preisach模型矢量加和所有單個磁滯算子的磁化強(qiáng)度。因而，準(zhǔn)確的磁滯模型和損耗模型有助于從鐵心磁化本質(zhì)上理解變壓器的磁化和損耗特性，解決相應(yīng)的工程問題如勵磁涌流問題［9,10］。基于磁滯模型的損耗計算方法，通過有限元方法耦合標(biāo)量或矢量磁滯模型，解法器迭代次數(shù)多，算法復(fù)雜計算量大［11,12］；基于Steinmetz公式類交變損耗和旋轉(zhuǎn)損耗方法，計算過程簡單但精度較低［13］，廣泛用于電機(jī)本體設(shè)計和效率分析［14］。傳統(tǒng)Steinmetz公式是穩(wěn)態(tài)計算方法，為了避免電磁有限元計算存儲每一步的計算結(jié)果，文獻(xiàn)［15］提出瞬態(tài)的Steinmetz計算方法，有效減少了計算量。無論是材料的精確磁滯模型還是Steinmetz公式的經(jīng)驗(yàn)?zāi)Ｐ停Ｐ蛻?yīng)用的前提是參數(shù)的辨識，首先獲得模型所需要的材料磁特性測量數(shù)據(jù)?？煽康牟牧夏Ｐ蛿?shù)據(jù)獲得依賴準(zhǔn)確的可重復(fù)的磁特性實(shí)驗(yàn)測量。電工鋼片一維國家測量標(biāo)準(zhǔn)是愛潑斯坦方圈方法；國際上一些科研團(tuán)隊研制了不同結(jié)構(gòu)形式的二維單片磁特性測量系統(tǒng)，根據(jù)勵磁結(jié)構(gòu)和磁化樣品的不同，可分為水平形式方形樣片旋轉(zhuǎn)測量儀［16］、水平圓形樣片旋轉(zhuǎn)測量儀［17］、垂直型旋轉(zhuǎn)測量儀等［18］。悉尼科技大學(xué)針對軟磁復(fù)合材料研發(fā)了首個三維軟磁材料磁特性測量儀［19］。國內(nèi)沈陽工業(yè)大學(xué)對硅鋼片軟磁材料進(jìn)行了二維磁特性測量［20］。河北工業(yè)大學(xué)研制了寬頻(20~1000Hz)大功率軟磁材料的三維磁特性測量裝置［21］。本文提出一種基于旋轉(zhuǎn)磁特性測量的瞬態(tài)矢量磁滯損耗模型。該模型統(tǒng)一描述了交變磁滯損耗和旋轉(zhuǎn)磁滯損耗，能夠預(yù)測對空間任意磁通密度軌跡的瞬時磁滯鐵耗。其基本原理是將空間磁通密度軌跡的運(yùn)動分解為切向運(yùn)動和法向運(yùn)動，其中法向磁通密度軌跡運(yùn)動采用交變損耗計算，切向磁通密度運(yùn)動采用旋轉(zhuǎn)損耗計算方法。利用三維磁特性測量裝置對無取向硅鋼片進(jìn)行交變損耗、旋轉(zhuǎn)損耗和變軸比的橢圓激勵的磁滯損耗測量。根據(jù)本文提出的矢量磁滯損耗模型，預(yù)測了變軸比的復(fù)合型橢圓激勵下的磁滯損耗，驗(yàn)證了該模型的有效性。1.1改進(jìn)矢量磁滯模型的定義本文基于傳統(tǒng)的矢量Play磁滯模型［12］提出考慮旋轉(zhuǎn)磁化的改進(jìn)型矢量磁滯模型。傳統(tǒng)的矢量Play模型依賴于圓形磁滯算子Hk(H)的疊加運(yùn)算，即式中，B為磁通密度；H為磁場強(qiáng)度；n為磁滯算子的數(shù)量；fk為非線性的映射函數(shù)從磁場強(qiáng)度空間映射到磁通密度矢量空間；Hk(H)為圓形的Play磁滯算子。磁場強(qiáng)度超前的空間距離即為磁滯算子半徑。圖1為軟磁材料橢圓磁化軌跡和改進(jìn)磁滯模型，圖中橢圓表示磁滯算子，橢圓形磁滯算子的邊界到算子中心的空間距離表示磁場強(qiáng)度的空間超前矢量，磁滯算子隨磁化軌跡的運(yùn)動而運(yùn)動。磁滯算子橢圓長軸和短軸在不同空間位置的長度表征磁性材料磁特性，根據(jù)二維或三維矢量磁特性測量建模和辨識。橢圓圓心表示無磁滯磁場強(qiáng)度Hnon，橢圓形算子的長軸或短軸之一位于橢圓圓心與坐標(biāo)原點(diǎn)的連接線上。算子的圓心位于磁化軌跡的圓弧上，表示無磁滯的磁場強(qiáng)度軌跡。磁場強(qiáng)度軌跡逆時針旋轉(zhuǎn)，相應(yīng)磁通密度軌跡也逆時針旋轉(zhuǎn)，兩者關(guān)系由磁場飽和、磁場各向異性和磁滯等因素決定。dB/dt的射線代表磁通密度矢量的空間變化方向，一維dB/dt表示磁通密度的增加或減少，從而決定磁場強(qiáng)度超前于磁通密度的方向。宏觀磁場強(qiáng)度Hhys分解為兩部分，即式中，Hcr為矯頑磁場強(qiáng)度；Hnon僅為磁通密度B的函數(shù)。Hcr是在橢圓中心點(diǎn)處作斜率為dB/dt的射線與磁滯算子的相交線段。當(dāng)交變勵磁時，dB/dt的射線沿軌跡的軸線方向；當(dāng)圓形勵磁時，射線沿圓形軌跡的切向方向；當(dāng)橢圓形勵磁時，射線位于兩者之間。Hnon計算式為對于各向同性材料，矢量Hnon方向與B方向相同，函數(shù)g的幅值可選為正切函數(shù)或Langevin函數(shù)；對于各向異性軟磁材料函數(shù)，g是空間二維矢量函數(shù)。綜上，矢量磁滯模型的計算流程如圖2所示。本文提出的改進(jìn)矢量磁滯模型與傳統(tǒng)的矢量Play模型區(qū)別如下：1） 矢量Play模型以磁場強(qiáng)度H為變量，而改進(jìn)模型以B為變量，因此后者也被稱為磁滯逆模型。當(dāng)進(jìn)行基于磁矢位A的有限元分析時，逆磁滯模型避免迭代運(yùn)算，提高了計算效率。2） 傳統(tǒng)矢量Play模型磁滯算子為圓形，改進(jìn)模型磁滯算子為橢圓形且考慮了旋轉(zhuǎn)磁化過程。3） 兩種模型中的非線性映射函數(shù)是一對互逆的函數(shù)。1.2矢量磁滯損耗分解軟磁材料的比損耗，即單位時間內(nèi)的磁場做功，計算式為Hnon磁場分量瞬時吸收或釋放能量，但一個周期內(nèi)做功的積分為零。磁通密度增加時儲存靜磁能，磁通密度減小時釋放靜磁能，本質(zhì)上不消耗能量，從而使磁性器件外特性為感性。磁滯能量的耗散僅來自于矯頑磁場Hcr部分，如圖1所示，矯頑磁場部分依賴橢圓形磁滯算子進(jìn)行插值計算，Hcr總是超前于無磁滯的磁場強(qiáng)度Hnon。商用軟件如Ansys-Maxwell中，矯頑磁場強(qiáng)度沿直角坐標(biāo)系正交分解為兩部分，即式中，Hxcr、Hycr分別是x、y方向的矯頑力。則總磁滯損耗相應(yīng)分解為兩部分之和，即直角坐標(biāo)系下的正交分解不能有效表征旋轉(zhuǎn)磁滯損耗的物理過程。旋轉(zhuǎn)磁化的鐵磁學(xué)機(jī)理是磁疇壁的旋轉(zhuǎn)，而在低場強(qiáng)的交變磁化時本質(zhì)是磁疇壁的移動。實(shí)驗(yàn)發(fā)現(xiàn)接近飽和時的旋轉(zhuǎn)磁滯損耗接近于零，因此理論上矯頑磁場強(qiáng)度也應(yīng)該為零。式(5)中，Hxcr、Hycr獨(dú)立采用交變矯頑力建模方法，故式(6)計算的兩個方向獨(dú)立的交變損耗之和大于旋轉(zhuǎn)磁化損耗，這與實(shí)驗(yàn)結(jié)果不符。為了解決該問題，本文將矯頑磁場強(qiáng)度分解為法向和切向兩部分。橢圓的兩個軸分別是切向軸和法向軸，法向軸與坐標(biāo)原點(diǎn)和橢圓中心連線重合，切向軸與連線垂直。沿法向軸分解的磁損耗受交變磁化影響，沿切向軸分解的磁損耗受旋轉(zhuǎn)磁化影響，則式中，ir、i。分別為橢圓軌跡的法向、切向單位方向；分別為法向、切向磁場矯頑力。下文分別利用幅值交變損耗和旋轉(zhuǎn)損耗計算。1.3幅值交變磁滯損耗建模幅值矯頑磁場強(qiáng)度的建模，根據(jù)等價的橢圓模型，假設(shè)勵磁波形是標(biāo)準(zhǔn)正弦，矯頑磁場強(qiáng)度符合余弦規(guī)律變化，則式中，Bm為磁化歷史記錄的最大磁通密度；Hm為最大的矯頑力；3為時變磁場的角頻率。聯(lián)立式(4)和式(8)，幅值交表磁滯損耗簡化為式中，f為勵磁頻率。根據(jù)Berotti損耗分離定律交變磁滯損耗假設(shè)得式中，kh為磁滯損耗系數(shù)。式⑺與式(6)聯(lián)立解得故瞬時交變磁滯損耗計算式為式中，Br為矢量B在橢圓軌跡法向方向的分量，當(dāng)交變磁化時Br為矢量B的幅值。1.4旋轉(zhuǎn)磁滯損耗建模圖3為軟磁材料圓形旋轉(zhuǎn)磁化過程磁場強(qiáng)度，磁通密度幅值Br不變，B滯后于H空間角度。，且。是圓形旋轉(zhuǎn)磁通密度幅值Br的函數(shù)。當(dāng)磁性材料的磁化軌跡為圓形時，瞬時功率式(4)簡化為根據(jù)圖3幾何關(guān)系的定義可知，是經(jīng)過磁滯算子中心從橢圓圓心做dB/dt方向的射線與磁滯算子相交所截取的線段，即在有限元計算中，Hnon可以通過磁通密度非線性關(guān)系直接計算，H通過磁滯模型計算，故切向計算式更實(shí)用。由于磁性材料SMC和無取向鋼片都存在一定的空間各向異性，當(dāng)磁通密度軌跡運(yùn)動時，空間角度差。不恒定。為了進(jìn)一步簡化損耗模型，取空間角度差的積分平均值作為損耗模型的特征參數(shù)，即式中，a為磁通密度矢量的空間角度。直接根據(jù)平均值計算損耗會出現(xiàn)一定的誤差，下面分析誤差來源并矯正該誤差。磁性材料旋轉(zhuǎn)圓形磁化的磁通密度軌跡為式中，B1為基波幅值。由于磁性材料的非線性，磁場強(qiáng)度H產(chǎn)生高次諧波的畸變?yōu)槭街校琀i為各奇次諧波幅值；&為各次諧波相位。單位時間平均功率的積分表達(dá)式為式中，。1為Hx的基波分量與Bx基波的相位差。式(18)說明對于圓形磁通密度的勵磁模型，旋轉(zhuǎn)損耗僅決定于基頻下的相角差01,Hx、Hy的高次諧波部分對磁滯損耗無貢獻(xiàn)。此時H1僅是H基波的幅值，故在式(12)增加修正系數(shù)，即當(dāng)樣品未達(dá)到飽和時，系數(shù)31;當(dāng)樣品過飽和時，系數(shù)k略小于1。進(jìn)一步給出損耗模型旋轉(zhuǎn)部分的計算式為綜上，通過式(12)和式(20)分別計算幅值交變損耗和切向瞬時旋轉(zhuǎn)損耗。文獻(xiàn)［7］中旋轉(zhuǎn)損耗的分解算法是在頻域下，需要有限元后處理存儲每一時刻每一網(wǎng)格位置的磁場。本文提出的算法在時域進(jìn)行磁通密度軌跡的分解，顯著節(jié)約內(nèi)存，降低計算根據(jù)本文提出的旋轉(zhuǎn)損耗模型對軟磁復(fù)合材料SMC-Somaloy500和無取向硅鋼片35WW270進(jìn)行模型參數(shù)的辨識和數(shù)值實(shí)現(xiàn)。2.1交變損耗模型參數(shù)辨識無論采用哪種動態(tài)磁特性測量方法，在交變或旋轉(zhuǎn)激勵下，磁性材料測量總鐵耗不僅包含磁滯損耗Ph而且含有經(jīng)典渦流損耗Pc和異常損耗Pe。隨著勵磁頻率的升高，渦流損耗和異常損耗在總損耗的占比增加起主導(dǎo)作用。僅在極低測量頻率(1Hz)時，可忽略剩余兩部分，認(rèn)為測量損耗僅是磁滯損耗。常用的實(shí)驗(yàn)損耗分離和處理方法是根據(jù)經(jīng)典的Steinmetz公式理論，將磁性材料總損耗Pv進(jìn)行理論分解，即式中，kh、kc、ke分別為交變磁滯損耗系數(shù)、渦流損耗系數(shù)和異常損耗系數(shù)。通過測量多個頻率的Bm-P損耗曲線，擬合模型損耗參數(shù)，求解最優(yōu)化問題，即優(yōu)化問題采用基因算法求解。對于無取向硅鋼片的優(yōu)化擬合，采用20Hz、50Hz、100Hz和200Hz的損耗測量數(shù)據(jù)；對于SMC材料，采用20Hz、100Hz和200Hz的損耗測量數(shù)據(jù)進(jìn)行辨識。圖4是無取向鋼片和SMC兩種材料20~200Hz的Bm-P損耗曲線，通過Steinmetz公式得到的計算值的擬合結(jié)果和測量結(jié)果吻合較好。SMC鐵耗在低磁通密度時遠(yuǎn)高于無取向鋼片，隨著頻率升高，鐵耗中渦流損耗占比不斷增大，兩者的總鐵耗差距縮小。當(dāng)交變勵磁頻率為200Hz時，SMC鐵耗高于35WW270—倍。表1為無取向鋼片和SMC材料的Steinmetz損耗系數(shù)，是通過在寬頻(20~200Hz)下的交變損耗曲線擬合得到。該模型認(rèn)為擬合的損耗系數(shù)不隨頻率的變化而變化，但是總損耗與頻率的關(guān)系是P倍頻、2倍頻和1.5倍頻三部分的組合關(guān)系。SMC材料的磁滯損耗系數(shù)遠(yuǎn)大于無取向硅鋼片，SMC異常損耗系數(shù)小于無取向鋼片。根據(jù)實(shí)驗(yàn)數(shù)據(jù)，兩者渦流損耗系數(shù)相差不大。SMC材料的渦流損耗機(jī)制復(fù)雜，不僅依賴于微小鐵基顆粒的尺寸而且依賴絕緣層的絕緣性能；既存在鐵基小顆粒內(nèi)部的微觀渦流，也存在穿過絕緣層的宏觀渦流［22］。對于硅鋼片材料，當(dāng)交變勵磁或旋轉(zhuǎn)勵磁時，硅鋼片內(nèi)部渦流Jx和Jy僅在沿疊片方向有顯著變化，在平行疊片的平面內(nèi)變化不顯著。任意復(fù)雜勵磁下，硅鋼疊片的三維渦流場問題可分解為x、y方向的獨(dú)立交變磁化過程，分別用一維渦流方程描述［23］，因此渦流損耗也可相應(yīng)地分解為x、y方向交變渦流損耗。類似地，假設(shè)異常損耗也可分解為x、y方向獨(dú)立的交變異常損耗，使用交變損耗辨識的kc、ke系數(shù)。故瞬時渦流損耗功率和異常損耗功率分別表示為根據(jù)式(23)和式(24)，從測量的總鐵耗中去除渦流損耗和異常損耗部分，剩余部分即不隨頻率變化的磁滯損耗，采用本文算法進(jìn)行模型驗(yàn)證。2.2旋轉(zhuǎn)損耗模型參數(shù)辨識軟磁材料的旋轉(zhuǎn)磁特性測量方法沒有統(tǒng)一的國際標(biāo)準(zhǔn)，國內(nèi)夕卜實(shí)驗(yàn)室廣泛采用控制圓形或橢圓形B軌跡的勵磁模型為式中，9為橢圓軌跡軸比。標(biāo)準(zhǔn)的圓形磁軌跡(9=1)不包含諧波，易于損耗建模和特征參數(shù)的提取，同時反映軟磁材料磁各向異性和旋轉(zhuǎn)磁損耗特性。當(dāng)9=0時，勵磁模型退化為單軸的交變勵磁，反映了軟磁材料的單軸交變損耗特性。在電力變壓器的T型接口處，磁通密度的軌跡接近菱形，可采用橢圓形磁通密度軌跡近似模擬。圖5對比了無取向硅鋼片35WW270和SMC材料在50Hz時圓形磁通密度激勵下的磁場強(qiáng)度軌跡。隨著磁通密度幅值的增大，磁場強(qiáng)度軌跡在45°左右凸出，相比0°和90°軸磁性能更弱，材料對外表現(xiàn)為較強(qiáng)的磁各向異性，而且無取向鋼片的磁各向異性強(qiáng)于SMC材料。因此，H矢量與B矢量的夾角0并不恒定，根據(jù)式(15)和磁特性實(shí)驗(yàn)繪制關(guān)系曲線，如圖6所示。隨著旋轉(zhuǎn)磁化幅值的增加，兩種材料的B-H空間矢量角都單調(diào)減小，其中無取向鋼片角度下降速度更快。在損耗建模中，磁通密度峰值與角度關(guān)系曲線采用分段插值的多項式函數(shù)光滑擬合。3.1三維磁特性測量三維磁測量裝置的主磁路由三對正交的C型鐵心和六個多層柔性勵磁線圈組成，如圖7所示［21］。勵磁裝置的測量樣品是SMC樣品或硅鋼疊片。樣品切割或疊置為邊長為22mm的立方體，放置于B-H復(fù)合立方體傳感箱內(nèi)，B-H線圈的尺寸和形狀如圖8所示，樣品四周相距1mm，是勻場保護(hù)層，保護(hù)層的材料和疊置方向與樣品相同，勻場保護(hù)層作為磁測量樣品的延伸，Ansys-Maxwell的磁場仿真結(jié)果表明保護(hù)層磁場分布不均勻，但在樣品表面獲得了均勻的磁場，從而避免磁場的邊緣效應(yīng)。三維磁測量系統(tǒng)的任意橢圓寬頻勵磁模型為式中，Ix、Iy分別為x、y軸勵磁電流；￡y為橢圓軌跡軸比。三維磁特性測量系統(tǒng)采用自適應(yīng)的數(shù)字控制算法［24］。以工控機(jī)為控制核心，通過多功能采集卡發(fā)出三路模擬控制信號，經(jīng)三路功率放大器和阻抗匹配電路放大獲得大勵磁電流，在測量樣品的內(nèi)部(SMC或硅鋼疊片)生成交變或旋轉(zhuǎn)的磁場，從而測量不同勵磁模式下軟磁材料的磁各向異性和損耗特性。3.2磁損耗的測量和模型驗(yàn)證本文提出的磁滯損耗計算模型同時兼顧交變磁滯損耗和旋轉(zhuǎn)磁滯損耗。采用三維磁特性測量裝置對無取向鋼片進(jìn)行多種勵磁模式的損耗測量。首先，通過一組交變磁化實(shí)驗(yàn)和一組圓形磁化實(shí)驗(yàn)辨識損耗模型的參數(shù)；其次，利用不同軸比橢圓軌跡下的損耗特性驗(yàn)證本文損耗模型的有效性。圖9對比了四種磁特性測量實(shí)驗(yàn)下的鐵心總損耗，磁滯損耗部分采用本文提出的損耗模型，經(jīng)典渦流損耗和異常損耗分別采用式(23)和式(24)進(jìn)行建模分析。根據(jù)勵磁模型式(25)，四組實(shí)驗(yàn)磁通密度軌跡分別是9=0、1、1/2、1/4,軌跡的形狀通過磁場反饋控制算法保證。圖9a是x軸交變損耗實(shí)驗(yàn)結(jié)果及兩種模型的預(yù)測結(jié)果，因?yàn)閮煞N模型的參數(shù)都依賴交變勵磁實(shí)驗(yàn)結(jié)果，故兩種模型都能夠很好地逼近測量結(jié)果。圖9b是圓形旋轉(zhuǎn)磁化的實(shí)驗(yàn)結(jié)果，本文模型的參數(shù)使用該組測量數(shù)據(jù)進(jìn)行辨識，故兩種模型與實(shí)驗(yàn)結(jié)果吻合很好。而基于Steinmetz的方法認(rèn)為旋轉(zhuǎn)磁化是兩個方向獨(dú)立交變磁化的疊加。鐵磁學(xué)的物理磁化過程是當(dāng)軟磁材料接近飽和時，旋轉(zhuǎn)磁化的磁性材料磁疇逐漸減少至消失，因此磁滯損耗部分接近于零，故旋轉(zhuǎn)磁特性實(shí)驗(yàn)損耗低于兩個方向交變的磁特性損耗之和，Steinmetz方法計算的預(yù)測旋轉(zhuǎn)損耗高于測量旋轉(zhuǎn)損耗，即為模型的主要缺陷。本文模型能夠直接從旋轉(zhuǎn)磁化角度建模分析，與實(shí)驗(yàn)結(jié)果吻合較好。圖9c和圖9d是不同軸比的橢圓激勵的測量結(jié)果，混合激勵下的磁疇壁運(yùn)動行為復(fù)雜，既有磁疇壁的平動也包含磁疇壁的轉(zhuǎn)動，基于Steinmetz的模型僅假設(shè)兩個方向獨(dú)立的磁疇壁移動，其預(yù)測損耗結(jié)果低于實(shí)際測量結(jié)果和本文模型預(yù)測結(jié)果。綜上，對于四種不同的勵磁模型，本文提出的損耗預(yù)測方法與實(shí)驗(yàn)結(jié)果誤差較小，尤其在旋轉(zhuǎn)勵磁下，明顯優(yōu)于傳統(tǒng)的Steinmetz方法。本文提出一種旋轉(zhuǎn)磁滯損耗模型，并對SMC材料和無取向硅鋼片進(jìn)行實(shí)驗(yàn)驗(yàn)證，詳細(xì)闡述了磁滯損耗模型原理及其實(shí)現(xiàn)過程。該模型有如下優(yōu)點(diǎn)：1） 統(tǒng)一了交變損耗和旋轉(zhuǎn)損耗模型，實(shí)現(xiàn)了任意軌跡二維磁滯損耗的計算。2） 損耗模型是基于瞬態(tài)實(shí)現(xiàn)的，無需存儲全過程的磁化歷史并計算處理，故占用內(nèi)存少，而且易于集成到電磁有限元分析中。3） 能夠計算過度飽和狀態(tài)下旋轉(zhuǎn)磁滯損耗，此時旋轉(zhuǎn)磁滯損耗顯著下降，鐵耗主要包括渦流損耗和剩余損耗兩部分，該方法比Steinmetz方法有更高的精度。張長庚男，1985年生，博士研究生，研究方向?yàn)楣こ屉姶艌雠c磁技術(shù)、三維磁特性測量與建模。E-mail:李永建男，1978年生，教授，研究方向?yàn)楣こ屉姶艌雠c磁技術(shù)、三維磁特性測量與建模。E-mail:.cn（ 通信作者）【相關(guān)文獻(xiàn)】佟文明，朱曉鋒,朱龍飛，等.不同供電方式對非晶合金永磁同步電機(jī)鐵耗的影響[J].電工技術(shù)學(xué)報,2015,30(10):10-18.TongWenming,ZhuXiaofeng,ZhuLongfei,etal.Theimpactofdifferentsupplymodesoncorelossesofamorphousalloypermanentmagnetsynchronousmotor[J].TransactionsofChinaElectrotechnicalSociety,2015,30(10):10-18.BerotttiG.Hysteresisinmagnetismforphysicists,materialsscientists,andengineers[M].NewYork:AcademicPress,1998.ReinertJ,BrockmeyerA,DonckerRD.Caculationoflossesinferro-andferromagneticmaterialsbasedonthemodifiedsteinmetzequation[J].IEEETransactionsonIndustryApplications,2001,37(4):1055-1061.AbdallahJ,SullivanC.Improvedcalculationofcorelosswithnonsinusoidalwaveforms[C]//Thirty-SixthIASAnnualMeeting,2001,4:2203-2210.IonelD,DellingerSJ,MillerTJE,etal.Onthevariationwithfluxandfrequencyofthecorelosscoefficientsinelectricalmachines[J].IEEETransactionsonIndustryApplications,2006,42(3):658-667.趙志剛,程志光,劉福貴，等.基于漏磁通補(bǔ)償?shù)膶?dǎo)磁鋼板直流偏磁雜散損耗特性模擬[J].電工技術(shù)學(xué)報,2015,30(8):38-44.ZhaoZhigang,ChengZhiguang,LiuFugui,etal.ModelingoftheDC-biasedstray-fieldlossofmagneticsteelplatebasedoncompensatorofleakageflux[J].TransactionsofChinaElectrotechnicalSociety,2015,30(8):38-44.ZhuJianguo,RamsdenVS.Improvedformulationsforrotationalcorelossesinrotatingelectricalmachines[J].IEEETransactionsonMagnetics,1998,34(4):2234-2242.李丹丹,劉福貴,李永建，等.一種新的混合矢量磁滯模型磁滯算子定義方法[J].電工技術(shù)學(xué)報，2014,29(11):1-7.LiDandan,LiuFugui,LiYongjian,elal.Anewdefinitionmethodofthehysteroninhybridvectorhysteresismodel[J].TransactionsofChinaElectrotech-nicalSociety,2014,29(11):1-7.傅偉，趙莉華,梁勇,等多臺變壓器空載合閘勵磁涌流及抑制方案的研究[J].電力系統(tǒng)保護(hù)與控制,2015,43(1):28-33.FuWei,ZhaoLihua,LiangYong,etal.Studyonno-loadclosinginrushcurrentoftransformergroupandsuppressionmeasures[J].PowerSystemProtectionandControl,2015,43(1):28-33.姚東曉，鄧茂軍,倪傳坤，等.變壓器多側(cè)勵磁涌流產(chǎn)生機(jī)理及對差動快速動作區(qū)影響研究[J].電力系統(tǒng)保護(hù)與控制，2016,44(5):36-41.YaoDongxiao,DengMaojun,NiChuankun,etal.Transformer'smulti-sideinrushcurrentgenerationmechanismanditsinfluenceonthedifferentialprotection'sfastactionzone[J].PowerSystemProtectionandControl,2016,44(5):36-41.EnokizonoM,SodaN.DirectmagneticlossanalysisbyFEMconsideringvectormagneticproperties[J].IEEETransactionsonMagnetics,1998,34(5):3008-3011.LinDingsheng,ZhouPing,BergqvistA.Improvedvectorplaymodelandparameteridentificationformagnetichysteresismaterials[J].IEEETransactionsonMagnetics,2014,50(2):357-60.趙小軍，崔燦,李琳,等基于定點(diǎn)諧波平衡法的鐵心磁滯與損耗特性分析[J].電工技術(shù)學(xué)報，2014,29(7):10-18.ZhaoXiaojun,CuiCan,LiLin,etal.Analysisofthehysteresisandlosscharacteristicsinthelaminatedcorebyfixedpointharmonic-balancedmethod[J].TransactionsofChinaElectrotechnicalSociety,2014,29(7):10-18.趙海森，羅應(yīng)立，劉曉芳,等.異步電機(jī)空載鐵耗分布的時步有限元分析[J].中國電機(jī)工程學(xué)報，2010,30(30):99-106.ZhaoHaisen,LuoYingli,LiuXiaofang,etal.Analysisonno-loadironlossesdistributionofasynchronousmotorswithtime-steppingfiniteelementmethod[J].ProceedingsoftheCSEE,2010,30(30):99-106.ZhouPing,LinDingsheng,LuChuan,etal.Anewalgorithmtoconsidertheeffectsofcorelosseson3-Dtransientmagneticfields[J].IEEETransactionsonMagnetic,2014,50(2):365-368.PfutznerH,MulasalihovicE,YamaguchiH,etal.Rotationalmagnetizationintransformercores-areview[J].IEEETransactionsonMagnetics,2011,47(11):4523-4533.SievertJ.Recentadvancesintheone-andtwo-dimensionalmagneticmeasurementtechniqueforelectricalsheetsteel[J].IEEETransactionsonMagnetics,1990,26(3):2553-2558.ZurekS