張會(huì)燕 王昱皓 張俊杰 司良英 祝麗花

摘要 并聯(lián)電抗器是電站設(shè)備中噪音的主要來(lái)源，設(shè)計(jì)和研究低噪聲電抗器具有重要意義。依據(jù)間隙鐵心式電抗器鐵心的結(jié)構(gòu)特點(diǎn)，需要螺栓的緊固作用，在工作狀態(tài)下電抗器的間隙式鐵心部分受壓應(yīng)力作用。為了準(zhǔn)確分析間隙鐵心式電抗器的工作磁特性和電磁振動(dòng)，開(kāi)展了不同應(yīng)力下多種鐵心硅鋼的磁化和磁致伸縮特性測(cè)量，并基于測(cè)量數(shù)據(jù)對(duì)電抗器樣機(jī)電磁振動(dòng)進(jìn)行了分析計(jì)算。文中首先改進(jìn)磁致伸縮測(cè)量系統(tǒng)，測(cè)量獲取首鋼3種典型鐵心用硅鋼在不同應(yīng)力，特別是壓應(yīng)力下的磁化和磁致伸縮特性;考慮鐵心磁致伸縮效應(yīng)，對(duì)不同硅鋼牌號(hào)制造的間隙電抗器鐵心，分別對(duì)其不同應(yīng)力下的磁場(chǎng)和機(jī)械場(chǎng)進(jìn)行分析計(jì)算，獲得了鐵心磁通密度和應(yīng)力的分布，以及不同位置的振動(dòng)情況。本文的研究成果可為低噪聲間隙鐵心式電抗器的設(shè)計(jì)提供不同硅鋼磁特性數(shù)據(jù)和理論分析方法。

關(guān) 鍵 詞 壓應(yīng)力;工作磁特性;間隙鐵心式電抗器;磁致伸縮;電磁振動(dòng)

Abstract Shunt reactors are the main sources of noise in power plant equipment. It is of great significance to design and study low noise reactors. According to the structural characteristics of the gap-core reactor core， the fastening effect of bolts is required， and the partial compressive stress of the gap-core of the reactor acts under the working state. In order to accurately analyze the working magnetic characteristics and electromagnetic vibration of the gap core reactor， the magnetization and magnetostriction characteristics of various core silicon steels under different stresses are measured， and the electromagnetic vibration of the reactor prototype is analyzed and calculated based on the measured data. Firstly， the magnetostrictive measurement system is improved to measure and obtain the magnetization and magnetostrictive characteristics of three typical silicon steels for iron cores in Shougang under different stresses， especially compressive stresses. Considering the magnetostrictive effect of iron core， the magnetic field and mechanical field of gap reactor iron cores made of different silicon steel grades under different stresses are analyzed and calculated respectively， and the distribution of magnetic flux density and stress as well as the vibration at different positions of the iron core are obtained. The research results can provide different silicon steel magnetic characteristic data and theoretical analysis methods for the design of low noise gap core reactor.

Key words compressive stress; working magnetic characteristics; gapped-core reactor; magnetostriction; electromagnetic vibration

0 引言

自2018年1月起，環(huán)保稅全面推進(jìn)，特高壓變電站被劃分為2類聲環(huán)境功能區(qū)，GB12348—2008中規(guī)定2類聲環(huán)境功能區(qū)要求廠界外噪聲限值為晝間60 dB（A），夜間50 dB（A）。相對(duì)于電站其他設(shè)備，并聯(lián)電抗器噪聲水平較高。并聯(lián)電抗器采用隔聲罩降噪單臺(tái)造價(jià)高，最嚴(yán)重的問(wèn)題是在失火情況下，BOX-in的密閉結(jié)構(gòu)將阻礙火勢(shì)初始段的最佳消防時(shí)機(jī)，延誤滅火;常規(guī)BOX-in中的玻璃棉也有助燃絕緣油火勢(shì)的作用。如果設(shè)計(jì)和研究一種低噪聲并聯(lián)電抗器本體，將噪音從源頭上降低下來(lái)，將更具有意義。

電工硅鋼片是廣泛應(yīng)用在電氣產(chǎn)品中的一種磁性材料，其磁特性關(guān)系到電氣產(chǎn)品的節(jié)能降耗和優(yōu)化設(shè)計(jì)[1]，準(zhǔn)確地掌握硅鋼片在實(shí)際應(yīng)用狀態(tài)下的磁特性，并根據(jù)產(chǎn)品應(yīng)用要求研制相應(yīng)的硅鋼片，是設(shè)計(jì)出性能優(yōu)良的電氣產(chǎn)品的必要條件。由于拉螺桿的緊固作用，間隙鐵心式并聯(lián)電抗器的鐵心餅在自然和工作狀態(tài)下都是在壓應(yīng)力作用下，然而壓應(yīng)力會(huì)嚴(yán)重影響硅鋼片的磁特性[2]，因此，為了研究應(yīng)力對(duì)硅鋼片磁特性的影響，需要對(duì)電工鋼片的磁特性進(jìn)行仔細(xì)地測(cè)量。安德森等測(cè)量了硅鋼片在磁通密度為1.5 T時(shí)不同應(yīng)力下的磁致伸縮量[3];賁彤[4]測(cè)量了硅鋼片在自由狀態(tài)及拉應(yīng)力下的磁特性;劉松[5]利用無(wú)取向硅鋼立方疊片三維磁特性并在此基礎(chǔ)上對(duì)飽和電抗器電磁振動(dòng)進(jìn)行了計(jì)算分析;Gao等[6-7]討論了采用高硬度絕緣材料填充電抗器鐵心間隙的降噪措施并且對(duì)電抗器鐵心拐角處進(jìn)行結(jié)構(gòu)優(yōu)化以降低損耗;Rossi等[8]分析了一臺(tái)帶有氣隙的U型電抗器，并通過(guò)磁致伸縮力與麥克斯韋力相互抵消來(lái)減小電抗器垂直方向的振動(dòng);張棟等[9]對(duì)大型并聯(lián)電抗器的漏磁場(chǎng)特性和電磁結(jié)構(gòu)特點(diǎn)進(jìn)行了系統(tǒng)討論，運(yùn)用節(jié)點(diǎn)力的方法計(jì)算了電抗器鐵心的受力情況;車福來(lái)等[10]設(shè)計(jì)了1臺(tái)并聯(lián)電抗器等效模型，此模型在各方面性能等效于電抗器產(chǎn)品，對(duì)產(chǎn)品結(jié)構(gòu)設(shè)計(jì)和工藝優(yōu)化有促進(jìn)作用，為產(chǎn)品級(jí)本體降噪新結(jié)構(gòu)應(yīng)用研究提供了基礎(chǔ)參數(shù)和測(cè)試實(shí)驗(yàn)平臺(tái);常晨等[11]基于等效模型對(duì)并聯(lián)電抗器的振動(dòng)和噪聲特性展開(kāi)研究，在不同預(yù)壓力和電壓下進(jìn)行振動(dòng)和噪聲多次實(shí)驗(yàn)，得出相應(yīng)分析結(jié)果，有助于并聯(lián)電抗器減振降噪、優(yōu)化設(shè)計(jì)及建模仿真等方面的工程應(yīng)用。

目前，針對(duì)間隙鐵心式電抗器鐵心工作磁特性的測(cè)量及電磁振動(dòng)分析相關(guān)研究還較少。本文首先改進(jìn)磁特性測(cè)量系統(tǒng)，測(cè)量了首鋼三種不同牌號(hào)硅鋼片在不同應(yīng)力下的磁化和磁致伸縮特性;基于測(cè)量數(shù)據(jù)，對(duì)比分析了間隙鐵心式電抗器鐵心使用不同牌號(hào)硅鋼片、不同壓應(yīng)力下的磁密及分布和電磁振動(dòng)，給出了鐵心的振動(dòng)信息和對(duì)比結(jié)果，為產(chǎn)品使用廠商提供基礎(chǔ)實(shí)用的工作特性，為設(shè)計(jì)更好的電抗器產(chǎn)品提供設(shè)計(jì)分析依據(jù)。

1 壓力下硅鋼磁特性測(cè)量系統(tǒng)及測(cè)量方法

德國(guó)BROCKHAUS的標(biāo)準(zhǔn)磁致伸縮測(cè)量系統(tǒng)，是基于單片磁性測(cè)量原理實(shí)現(xiàn)的，如圖1所示。該磁致伸縮測(cè)量系統(tǒng)包括主機(jī)、電源、單片磁測(cè)量裝置、氣動(dòng)隔振平臺(tái)、激光干涉儀、氣動(dòng)預(yù)應(yīng)力系統(tǒng)、空氣壓縮機(jī)。主機(jī)是用于設(shè)置硅鋼片參數(shù)、線圈參數(shù)、激勵(lì)形式、激勵(lì)頻率、去磁頻率、激勵(lì)大小等參數(shù)。單片磁測(cè)量裝置主要由線圈和上下鐵軛組成，被試電工硅鋼片放入線圈中。激光干涉儀是用來(lái)測(cè)量電工硅鋼片在磁場(chǎng)中的長(zhǎng)度變化量，分辨率為10 nm/m。

被試樣片的規(guī)格為600 mm × 100 mm，如圖2所示。定義矩形硅鋼樣片的長(zhǎng)度方向?yàn)檐堉品较颍≧olling Direction，RD），寬度方向?yàn)榇怪避堉品较颍═ransverse Direction，TD）。

單片硅鋼片所受的不同預(yù)應(yīng)力大小由預(yù)應(yīng)力控制器顯示和調(diào)節(jié)，預(yù)應(yīng)力控制器如圖3所示。通過(guò)控制樣片軋制方向的外加應(yīng)力大小，實(shí)現(xiàn)不同預(yù)應(yīng)力下樣片的磁化和磁致伸縮特性測(cè)量。

滿足標(biāo)準(zhǔn)的磁致伸縮測(cè)量裝置通常是測(cè)量硅鋼片單端自由狀態(tài)下的磁致伸縮，本臺(tái)磁致伸縮測(cè)量裝置也僅能夠施加少量壓應(yīng)力。當(dāng)在單片硅鋼片上施加較大壓應(yīng)力時(shí)，鋼片會(huì)彎曲變形，嚴(yán)重影響鋼片上的磁場(chǎng)和磁特性。本文對(duì)德國(guó)BROCKHAUS的標(biāo)準(zhǔn)磁致伸縮測(cè)量裝置進(jìn)行了部分改造，使其能夠施加10 MPa范圍內(nèi)的壓應(yīng)力。具體方法是：在樣片上放入適當(dāng)厚度的環(huán)氧樹(shù)脂板，并用聚酯薄膜調(diào)節(jié)厚度，以填充樣片和線圈之間的空隙;在硅鋼片無(wú)支撐的部分放入支撐墊塊，并用上鐵軛壓住，以確保測(cè)量過(guò)程中測(cè)量樣片不發(fā)生彎曲變形，墊塊與樣片接觸的表面涂鉛筆芯以減小摩擦力。拉應(yīng)力測(cè)量時(shí)要保證樣片緊貼下鐵軛，上鐵軛懸起，系統(tǒng)如圖4所示，壓應(yīng)力測(cè)量時(shí)系統(tǒng)如圖5所示。

2 不同壓應(yīng)力下硅鋼片磁特性的測(cè)量分析

本文利用符合國(guó)家和國(guó)際標(biāo)準(zhǔn)的Epstein 700測(cè)量工頻50 Hz下硅鋼片的磁化特性，所用樣片為首鋼生產(chǎn)的3種典型硅鋼牌號(hào)，其樣片牌號(hào)分別為20SQGD070、27SQGD085、30SQG120。牌號(hào)為20SQGD070的硅鋼片質(zhì)量為88.5 g，密度為7650 g/dm3，截面面積為19.28 mm2;牌號(hào)為27SQGD085的硅鋼片質(zhì)量為117.6 g，密度為7650 g/dm3，截面面積為25.62 mm2;牌號(hào)為30SQG120的硅鋼片質(zhì)量為132 g，密度為7 650 g/dm3，截面面積為28.76 mm2。

通過(guò)磁致伸縮測(cè)量裝置對(duì)上述3種硅鋼片施加不同拉力和壓力，測(cè)量了磁化特性和磁致伸縮特性，測(cè)量結(jié)果如圖6～圖8所示。圖6測(cè)量結(jié)果表明，拉應(yīng)力下20SQGD070牌號(hào)硅鋼片磁化特性變化較小，拉應(yīng)力對(duì)磁致伸縮特性的影響比較小，而壓應(yīng)力對(duì)磁致伸縮特性影響較大。圖6b）蝴蝶曲線表明最大壓應(yīng)力下，磁致伸縮是負(fù)磁致伸縮。由圖6c）可知樣品的磁致伸縮峰峰值隨壓應(yīng)力增大而增大，特別是壓應(yīng)力超過(guò)7 MPa后磁致伸縮增加更為明顯，遠(yuǎn)大于無(wú)應(yīng)力下的磁致伸縮值。該牌號(hào)硅鋼片，磁致伸縮性微應(yīng)變量值不高，性能優(yōu)異，特別在并聯(lián)電抗器工作點(diǎn)1 T左右的低磁密區(qū)域，壓應(yīng)力下磁致伸縮微應(yīng)變量值也不高。

27SQGD085牌號(hào)硅鋼不同應(yīng)力下磁致伸縮測(cè)量結(jié)果如圖7所示，測(cè)量結(jié)果顯示樣片所受應(yīng)力超過(guò)6 MPa后磁致伸縮峰峰值變化較大，且數(shù)值較20SQGD070牌號(hào)大。30SQG120牌號(hào)硅鋼的測(cè)量結(jié)果如圖8所示，測(cè)量結(jié)果同樣表明硅鋼樣品磁化特性隨拉應(yīng)力變化較小，幾乎可以忽略，磁致伸縮受拉應(yīng)力影響可以忽略，受壓應(yīng)力影響較大;該樣品從壓應(yīng)力超過(guò)5 MPa后磁致伸縮峰峰值顯著增大，數(shù)值比27SQGD085牌號(hào)更大。不同于20SQGD070牌號(hào)，27SQGD085、30SQG120牌號(hào)的蝴蝶曲線表明，隨著壓力的增大，磁致伸縮由負(fù)磁致伸縮逐漸變?yōu)檎胖律炜s。

由圖6～圖8可看出不同牌號(hào)的硅鋼片，拉應(yīng)力對(duì)磁化特性影響不大;對(duì)于磁致伸縮特性，拉應(yīng)力越大，磁致伸縮量略減小，壓應(yīng)力越大，磁致伸縮量越大，超過(guò)一定壓力后，磁致伸縮值遠(yuǎn)大于無(wú)應(yīng)力下的磁致伸縮值。因此，電抗器用的硅鋼片不能用常規(guī)無(wú)應(yīng)力下的磁致伸縮性能來(lái)表征和篩選，某些特制的低磁致伸縮硅鋼片可能不適用于并聯(lián)電抗器。上述結(jié)果也表明，壓力下的磁致伸縮量隨著硅鋼片的厚度增加而增加。這是因?yàn)樵诖艌?chǎng)的作用下，磁性材料中的疇壁將會(huì)移動(dòng)，磁疇將發(fā)生旋轉(zhuǎn)。對(duì)于正磁致伸縮材料，張力方向是易磁化方向，即硅鋼片在拉應(yīng)力方向易磁化，拉伸應(yīng)力將使磁化沿拉伸方向上升，故磁通密度也隨拉應(yīng)力的增加而增加，而壓應(yīng)力反之。同時(shí)，磁致伸縮的方向主要沿磁化方向，為了使磁疇發(fā)生旋轉(zhuǎn)，拉應(yīng)力需要克服磁各向異性，所以磁致伸縮隨拉應(yīng)力的增大而減小，壓應(yīng)力則反之，這與其他研究者測(cè)量和理論分析一致[2，4]。

3 間隙鐵心式電抗器電磁振動(dòng)分析

4 間隙鐵心式電抗器模型計(jì)算與分析

4.1 鐵心磁通密度及應(yīng)力分布

本文以一臺(tái)10 kV、1 250 kVar的間隙鐵心式電抗器作為研究對(duì)象，其相關(guān)參數(shù)值如下表1所示。該電抗器本體結(jié)構(gòu)由鐵軛、鐵心餅、氣隙墊塊、繞組組成，為了提高求解精度和求解速度，對(duì)間隙鐵心部分進(jìn)行了細(xì)化剖分，利用有限元方法對(duì)其求解，得到其磁通密度分布如圖9所示，該電抗器設(shè)計(jì)工作磁密為1.2 T，由圖9的計(jì)算結(jié)果表明計(jì)算磁通密度分布正確。在磁場(chǎng)準(zhǔn)確計(jì)算的基礎(chǔ)上，以應(yīng)力下測(cè)量數(shù)據(jù)為基礎(chǔ)對(duì)間隙電抗器主體部件進(jìn)行固體力學(xué)分析，在進(jìn)行力學(xué)分析時(shí)，以邊柱底端為固定約束，鐵心在4 MPa壓應(yīng)力下的應(yīng)力分布如圖10所示。由于約束和間隙部分的夾緊應(yīng)力作用，圖10所示底柱鐵心的應(yīng)力分布較大。

4.2 數(shù)值計(jì)算結(jié)果

為更準(zhǔn)備表達(dá)鐵心的電磁振動(dòng)信息，選取鐵心中不同點(diǎn)，如圖9所示。選取牌號(hào)為27SQGD085的硅鋼片作為間隙電抗器鐵心材料，應(yīng)用壓應(yīng)力下測(cè)量的硅鋼樣片磁致伸縮特性數(shù)據(jù)，選取如圖7所示特性曲線中4 MPa應(yīng)力時(shí)作為計(jì)算的特性數(shù)據(jù)，計(jì)算得到鐵心A、B、C? 3處的振動(dòng)情況，3點(diǎn)的振動(dòng)加速度如圖11所示，其中A點(diǎn)處振動(dòng)最大，其次是C點(diǎn)，振動(dòng)最小的為B點(diǎn)。結(jié)合圖9磁通密度分布圖可看出，鐵心磁通密度分布與鐵心振動(dòng)情況密切相關(guān)[12]，磁通密度越集中的地方，其振動(dòng)加速度越大，鐵心應(yīng)力也越大。

由于螺栓的緊固作用，電抗器鐵心部分自然工作狀態(tài)是在壓應(yīng)力下，為了研究壓應(yīng)力對(duì)間隙鐵心式電抗器鐵心在工作狀態(tài)下的振動(dòng)影響，先選取某一牌號(hào)硅鋼片在不同應(yīng)力下的磁特性數(shù)據(jù)作為鐵心部分?jǐn)?shù)據(jù)，從而對(duì)鐵心部分振動(dòng)情況進(jìn)行分析。本文選取牌號(hào)為27SQGD085的硅鋼片的磁特性作為裂心鐵心部分的磁特性。將測(cè)得的0、-4.22 MPa、-6.25 MPa、? -7.88 MPa、-10 MPa這幾種壓應(yīng)力下的磁特性數(shù)據(jù)代入有限元分析軟件中進(jìn)行計(jì)算，得到了電抗器鐵心中點(diǎn)A在不同壓應(yīng)力下的振動(dòng)情況，其振動(dòng)加速度和位移的計(jì)算結(jié)果如圖12所示。

由圖12可得，間隙鐵心式電抗器裂心部分在不同壓應(yīng)力下具有不同的振動(dòng)加速度，且振動(dòng)加速度值隨壓應(yīng)力的增加而增大，由圖6～圖8也可看出磁致伸縮峰峰值大小也與壓應(yīng)力值有關(guān)，即隨著壓應(yīng)力的增加，硅鋼片的磁致伸縮峰峰值的增長(zhǎng)基本保持一致，所以硅鋼片的磁致伸縮峰峰值越大，鐵心的振動(dòng)加速度值越大，因此在電抗器的電磁振動(dòng)分析設(shè)計(jì)時(shí)過(guò)大的夾緊力會(huì)增加電抗器的振動(dòng)。

由于不同牌號(hào)的硅鋼片具有不同的磁特性，為了研究不同牌號(hào)下鐵心在相關(guān)工作狀態(tài)下的振動(dòng)情況對(duì)比，將文中3組牌號(hào)的硅鋼片在壓應(yīng)力為4 MPa下的磁特性數(shù)據(jù)作為鐵心磁特性數(shù)據(jù)基礎(chǔ)，分別實(shí)現(xiàn)對(duì)其的電磁振動(dòng)分析，3種牌號(hào)硅鋼鐵心在相同位置（圖9中A點(diǎn)）的振動(dòng)情況如圖13所示。計(jì)算結(jié)果表明，不同牌號(hào)硅鋼片的鐵心具有不同的振動(dòng)加速度，其振動(dòng)加速度大小與硅鋼片的磁致伸縮峰峰值有關(guān)，其磁致伸縮峰峰值越大，鐵心的振動(dòng)加速度和振動(dòng)位移越大，因此，在電抗器設(shè)計(jì)及制造中，選用磁致伸縮較小的硅鋼片可以適當(dāng)降低鐵心的振動(dòng)情況。

5 結(jié)論

本文首先通過(guò)改進(jìn)磁致伸縮測(cè)量裝置實(shí)現(xiàn)了3種首鋼典型硅鋼片在不同應(yīng)力下的磁特性測(cè)量分析，測(cè)量結(jié)果顯示應(yīng)力對(duì)硅鋼片的磁特性有影響：拉應(yīng)力越大，硅鋼片的導(dǎo)磁性越好，磁致伸縮量越小;壓應(yīng)力越大，硅鋼片的磁致伸縮量越大?；陔娍蛊麒F心的三維磁-機(jī)械耦合模型，將測(cè)得的磁特性數(shù)據(jù)代入并進(jìn)行計(jì)算，計(jì)算結(jié)果表明，壓應(yīng)力對(duì)鐵心的振動(dòng)具有一定的影響，同一牌號(hào)的硅鋼片，壓應(yīng)力值越大，鐵心的振動(dòng)越大;不同牌號(hào)的硅鋼片，其磁致伸縮峰峰值越大，鐵心的振動(dòng)越大。因此，在進(jìn)行電抗器的設(shè)計(jì)時(shí)，盡量選擇磁致伸縮峰峰值小的硅鋼片牌號(hào)，同時(shí)，要避免給鐵心施加過(guò)大的夾緊力，但是對(duì)于電抗器本身的結(jié)構(gòu)來(lái)說(shuō)，夾緊力又是必需的，所以應(yīng)采取適當(dāng)?shù)膴A緊力。

參考文獻(xiàn)：

[1]? ? 謝德馨，白保東. 計(jì)算電磁學(xué)中電工鋼片磁特性模型研究的新進(jìn)展[J]. 沈陽(yáng)工業(yè)大學(xué)學(xué)報(bào)，2007，29（3）：289-294.

[2]? ? YUICHIRO K，YUJI T，TAKASHI T，et al. Influence of stress on vector magnetic property under alternating magnetic flux conditions. IEEE Transactions on Magnetics，2011，47（10）：4344-4347.

[3]? ? PHILIP I A，ANTHONY J M，HUGH J S. Assessment of the stress sensitivity of magnetostriction in grain-oriented. IEEE Transactions on Magnetics，2007，43（8）：3467-3476.

[4]? ? BEN T，YANG Q X，YAN R G，et al. Magnetically controlled saturable reactor core vibration under practical working conditions. IEEE Transactions on Magnetics，2017，53：1-4.

[5]? ? 劉松，祝麗花，王斌，等. 基于三維磁特性電抗器電磁振動(dòng)分析[J]. 電氣應(yīng)用，2019，38（3）：35-40.

[6]? ? GAO Y H，NAGATA M，MURAMATSU K，et al. Noise reduction of a three-phase reactor by optimization of gaps between cores considering electromagnetism and magnetostriction[J]. IEEE Transactions on Magnetics，2011，47（10）：2772-2775.

[7]? ? GAO Y H，TAURA D，NAGATA M，et al. Loss reduction of reactor with grain-oriented silicon steel plates[J]. IEEE Transactions on Magnetics，2013，49（5）：1973-1976.

[8]? ? ROSSI M，LE BESNERAIS J. Vibration reduction of inductors under magnetostrictive and maxwell forces excitation[J]. IEEE Transactions on Magnetics，2015，51（12）：1-6.

[9]? ? 張棟，張俊杰，張喜樂(lè)，等. 大型并聯(lián)電抗器的電磁場(chǎng)綜述及受力的仿真[J]. 變壓器，2015，52（10）：19-25.

[10]? 車福來(lái)，張俊杰，劉蘭榮，等. 并聯(lián)電抗器本體降噪技術(shù)的等效模型研究[J]. 黑龍江電力，2018，40（6）：485-491.

[11]? 常晨，劉蘭榮，盧美林，等. 并聯(lián)電抗器等效模型的振動(dòng)噪聲特性試驗(yàn)研究[J]. 電氣技術(shù)，2019，20（2）：37-41，47.

[12]? 趙小軍，崔燦，李琳，等. 基于定點(diǎn)諧波平衡法的鐵心磁滯與損耗特性分析[J]. 電工技術(shù)學(xué)報(bào)，2014，29（7）：10-18.