在非線性磁場(chǎng)的求解和分析中，可以通過(guò)引入不同的位和位組，簡(jiǎn)化計(jì)算或提高效率?；诖攀噶课籄的有限元法在二維、三維非線性磁場(chǎng)問(wèn)題的求解中得到了廣泛的應(yīng)用。然而在三維場(chǎng)的計(jì)算中，由于每一個(gè)節(jié)點(diǎn)上A有三個(gè)方向分量，使得總未知量相對(duì)較多，進(jìn)而占用較多內(nèi)存；相對(duì)而言，基于磁標(biāo)量位Ф的計(jì)算則有未知量總數(shù)少的優(yōu)點(diǎn)。

另一方面，磁場(chǎng)問(wèn)題中作為激勵(lì)的外加電流的處理十分關(guān)鍵，尤其是在考慮場(chǎng)路耦合的條件下。一般直接在激勵(lì)源網(wǎng)格中施加電流密度或者通過(guò)積分的方法計(jì)算電流對(duì)簡(jiǎn)化位的貢獻(xiàn)，前者難以保證激勵(lì)源與實(shí)體網(wǎng)格的準(zhǔn)確吻合，后者則需要額外的積分運(yùn)算，且不易處理場(chǎng)路耦合關(guān)系。

已有學(xué)者針對(duì)基于磁標(biāo)量位的有限元法進(jìn)行了研究。

• 有學(xué)者提出了對(duì)靜態(tài)磁場(chǎng)進(jìn)行求解的單標(biāo)量磁位法，并對(duì)相關(guān)模型進(jìn)行了求解。
• 有學(xué)者對(duì)簡(jiǎn)化標(biāo)量位法、雙標(biāo)量位法和差場(chǎng)標(biāo)量位法三種方法進(jìn)行了分析，可以對(duì)一些問(wèn)題進(jìn)行比較準(zhǔn)確的求解。
• 有學(xué)者使用雙標(biāo)量磁位法對(duì)開(kāi)域靜磁場(chǎng)進(jìn)行了研究。
• 有學(xué)者給出了有源區(qū)域的處理方法并應(yīng)用到三維電磁場(chǎng)的求解中。

但上述研究大多數(shù)針對(duì)無(wú)源區(qū)域的磁場(chǎng)進(jìn)行分析，在有源區(qū)域的磁場(chǎng)需要通過(guò)積分的方法得到，在處理問(wèn)題的過(guò)程中具有比較大的局限性。國(guó)外專家學(xué)者對(duì)磁標(biāo)量位問(wèn)題的研究相對(duì)較多。L. F. Yann等基于磁標(biāo)量位對(duì)三維磁場(chǎng)的計(jì)算方法進(jìn)行了深入研究。O.Bíró等在時(shí)域下采用磁標(biāo)量位對(duì)渦流問(wèn)題進(jìn)行了計(jì)算和分析。M. Gérard等基于磁標(biāo)量位對(duì)塊狀導(dǎo)體中的場(chǎng)路耦合關(guān)系和多連通域的處理方法進(jìn)行了研究。

太陽(yáng)磁暴的周期性變化產(chǎn)生的準(zhǔn)直流、高壓直流輸電中單極大地回路方式或雙極以不平衡方式運(yùn)行和電力系統(tǒng)中電力電子器件特性的影響均會(huì)導(dǎo)致直流偏磁現(xiàn)象的發(fā)生。直流偏磁的發(fā)生會(huì)使變壓器工作在半飽和狀態(tài)，引起勵(lì)磁電流發(fā)生畸變、損耗增加、溫升增加和振動(dòng)加劇等后果，嚴(yán)重影響變壓器正常運(yùn)行以及電力系統(tǒng)的安全穩(wěn)定。

直流偏磁條件下變壓器鐵心內(nèi)的磁場(chǎng)求解屬于典型的強(qiáng)非線性問(wèn)題。有相關(guān)研究工作采用A在時(shí)域和頻域?qū)χ绷髌艈?wèn)題進(jìn)行計(jì)算和分析。當(dāng)系統(tǒng)時(shí)間常數(shù)較大時(shí)，時(shí)步法需要較長(zhǎng)的過(guò)渡時(shí)間，而頻域法則可以克服這方面的缺陷，同時(shí)還能夠與定點(diǎn)技術(shù)相結(jié)合，在頻域內(nèi)對(duì)各次諧波進(jìn)行并行計(jì)算。

本文采用基于復(fù)指數(shù)級(jí)數(shù)的諧波平衡法建立基于磁標(biāo)量位的頻域非線性磁場(chǎng)方程，引入電矢量位對(duì)激勵(lì)源進(jìn)行處理從而有效地簡(jiǎn)化計(jì)算。將定點(diǎn)技術(shù)與頻域算法相結(jié)合，考慮場(chǎng)路耦合關(guān)系，在頻域內(nèi)對(duì)直流偏磁磁場(chǎng)和勵(lì)磁電流進(jìn)行求解。利用疊片鐵心模型進(jìn)行直流偏磁實(shí)驗(yàn)，將基于磁標(biāo)量位的計(jì)算結(jié)果同實(shí)驗(yàn)結(jié)果進(jìn)行比較，驗(yàn)證了算法的有效性。

圖1 疊片鐵心實(shí)驗(yàn)?zāi)Ｐ?/p>

總結(jié)

• 1）提出了基于磁標(biāo)量位的定點(diǎn)頻域分解算法，對(duì)疊片鐵心模型的直流偏磁問(wèn)題進(jìn)行求解和分析。通過(guò)與實(shí)驗(yàn)測(cè)量結(jié)果和其他數(shù)值方法計(jì)算結(jié)果的對(duì)比分析，驗(yàn)證了本文算法的有效性和準(zhǔn)確性。
• 2）通過(guò)引入定點(diǎn)磁導(dǎo)率實(shí)現(xiàn)非線性系數(shù)矩陣的線性化處理，對(duì)各次諧波進(jìn)行解耦，降低了計(jì)算中系統(tǒng)占用的內(nèi)存，且適用于并行和分布式計(jì)算。在強(qiáng)加電流區(qū)域，引入與單位電流激勵(lì)對(duì)應(yīng)的電矢量位t0，利用其線性特征，基于棱單元進(jìn)行插值計(jì)算。對(duì)于形狀規(guī)則的繞組，該方法能夠有效地實(shí)現(xiàn)場(chǎng)路耦合建模。
• 3）基于磁標(biāo)量位的頻域分解算法可以對(duì)直流偏磁條件下變壓器內(nèi)部的磁場(chǎng)進(jìn)行較為準(zhǔn)確和快速的計(jì)算，既避免了時(shí)步法中的過(guò)渡過(guò)程，又可以為研究由直流偏磁引起的變壓器鐵心損耗、局部溫升、振動(dòng)加劇及相關(guān)多物理場(chǎng)耦合問(wèn)題提供重要支撐。
• 4）本文中的計(jì)算模型為二維模型，磁標(biāo)量位占用內(nèi)存少的優(yōu)勢(shì)難以體現(xiàn)，但其在三維電磁場(chǎng)數(shù)值計(jì)算中的優(yōu)勢(shì)將十分明顯。后續(xù)研究工作需要進(jìn)一步實(shí)現(xiàn)針對(duì)三維模型的計(jì)算程序，并實(shí)現(xiàn)本文方法在解決非線性渦流問(wèn)題時(shí)的有效性分析及其驗(yàn)證。