龔文軍，　黃劭剛，　楊玉文，　饒利強(qiáng)，　曾　裕

(南昌大學(xué) 信息工程學(xué)院，江西 南昌　330031)

?

ANSYS Maxwell硅鋼片鐵耗系數(shù)提取方法的實(shí)現(xiàn)*

龔文軍，黃劭剛，楊玉文，饒利強(qiáng)，曾裕

(南昌大學(xué) 信息工程學(xué)院，江西 南昌330031)

摘要：詳細(xì)介紹了Maxwell軟件的鐵耗計(jì)算模型以及鐵耗參數(shù)提取計(jì)算方法及步驟。采用最小二乘法對(duì)硅鋼片材料M19_29G的鐵損曲線進(jìn)行擬合分析，得到了可直接輸入Maxwell進(jìn)行計(jì)算的鐵耗系數(shù)。采用Maxwell自帶的鐵耗參數(shù)提取功能驗(yàn)證了計(jì)算結(jié)果的正確性。為Maxwell軟件鐵耗系數(shù)提取的實(shí)現(xiàn)提供了另外一種途徑，使得數(shù)據(jù)處理過(guò)程更加開放、靈活，不再受限于軟件本身。所述Maxwell軟件鐵耗參數(shù)提取過(guò)程對(duì)其他商用有限元軟件的鐵耗系數(shù)提取具有重要的參考價(jià)值。

關(guān)鍵詞：ANSYS Maxwell； 最小二乘法； 鐵心損耗； 鐵耗系數(shù)

0引言

鐵磁材料是電機(jī)的重要組成部分，鐵耗在總損耗中占有較大比例，與效率直接相關(guān)，是影響電機(jī)性能的一個(gè)敏感因素。因此，鐵耗的準(zhǔn)確計(jì)算對(duì)電機(jī)性能分析具有重大意義。傳統(tǒng)的方法是將齒部和軛部?jī)啥未怕贩謩e進(jìn)行計(jì)算[1]，并用適當(dāng)?shù)慕?jīng)驗(yàn)系數(shù)進(jìn)行修正。對(duì)傳統(tǒng)電機(jī)而言，該方法具有一定的準(zhǔn)確性。但隨著新型電機(jī)的不斷發(fā)展，以及對(duì)電機(jī)參數(shù)精確度更高要求的背景下，傳統(tǒng)計(jì)算方法顯得心有余而力不足。有限元法可對(duì)磁密不均勻分布、諧波磁場(chǎng)、旋轉(zhuǎn)磁化、機(jī)械加工等因素進(jìn)行計(jì)算處理[2]，因此，相對(duì)傳統(tǒng)磁路法計(jì)算，有限元計(jì)算具有計(jì)及因素全面、準(zhǔn)確度高等優(yōu)點(diǎn)，已成為電機(jī)研究領(lǐng)域的一個(gè)熱點(diǎn)。利用有限元法對(duì)電機(jī)鐵耗進(jìn)行計(jì)算，鐵耗系數(shù)的準(zhǔn)確提取則顯得尤為重要。

文獻(xiàn)[3]對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行擬合分析，得到了鐵耗分離模型中的鐵耗系數(shù)，但得到的鐵損系數(shù)無(wú)法直接應(yīng)用于商業(yè)有限元軟件鐵耗計(jì)算；文獻(xiàn)[4]介紹了Ansoft Maxwell軟件鐵耗計(jì)算原理和方法，但沒(méi)有給出具體步驟。本文針對(duì)Maxwell軟件所采用的鐵耗計(jì)算模型，對(duì)硅鋼片材料的鐵損曲線進(jìn)行擬合分析，得到了Maxwell軟件鐵耗計(jì)算所需系數(shù)，最后用實(shí)例驗(yàn)證了本方法的可行性。本文僅討論Maxwell軟件的鐵耗系數(shù)提取過(guò)程，且不考慮直流偏磁效應(yīng)。

1硅鋼片材料的建模

鐵耗的基本計(jì)算模型主要有以下兩種：

1.1常系數(shù)兩項(xiàng)式模型

常系數(shù)兩項(xiàng)式模型是電機(jī)鐵耗計(jì)算中的基本鐵耗模型?；捐F耗是由主磁場(chǎng)在鐵心內(nèi)發(fā)生變化時(shí)產(chǎn)生的，主要計(jì)及磁滯損耗和渦流損耗兩項(xiàng)。

PFe=KhfBmα+Kc(fBm)2

(1)

式中：PFe——鐵耗；

Kh——磁滯損耗系數(shù)；

f——頻率；

Bm——磁通密度幅值；

α——Steinmetz系數(shù)，取決于材料本身特性；

Kc——經(jīng)典渦流損耗系數(shù)。

常系數(shù)兩項(xiàng)模型主要用于工程上的粗略估算，較少有限元被采用。

1.2常系數(shù)三項(xiàng)式模型

常系數(shù)三項(xiàng)式模型為意大利學(xué)者Bertotti于1988年總結(jié)提出的。按照Bertottti鐵耗分離計(jì)算模型，不考慮集膚效應(yīng)時(shí)，鐵心損耗可以用式(2)表達(dá)，該損耗為正弦交變磁場(chǎng)下單位重量硅鋼片產(chǎn)生的損耗：

Pir=Ph+Pc+Pe=

khfBmα+kc(fBm)2+ke(fBm)1.5

(2)

式中：Pir——鐵耗；

Ph——磁滯損耗；

Pc——經(jīng)典渦流損耗；

Pe——附加損耗；

f——頻率；

Bm——磁密幅值；

kh、kc、ke、α——與材料有關(guān)的系數(shù)。

常系數(shù)三項(xiàng)式模型對(duì)鐵耗進(jìn)行了分離，包括磁滯、渦流及附加損耗三項(xiàng)。雖然Bertotti鐵耗分離模型在一定程度上能較為準(zhǔn)確地計(jì)算鐵心損耗，但在一些特殊應(yīng)用場(chǎng)合誤差仍較高，因此，在Bertotti鐵耗分離模型的基礎(chǔ)上衍生出了一些改進(jìn)計(jì)算模型。比如，針對(duì)非正弦供電，文獻(xiàn)[6]在Bertotti分立鐵耗計(jì)算模型的基礎(chǔ)上，給出了PWM逆變器供電下硅鋼片損耗的解析計(jì)算模型。文獻(xiàn)[7]針對(duì)變頻電機(jī)鐵耗分析中存在的鐵磁材料建模，以及非正弦激勵(lì)時(shí)鐵耗計(jì)算模型等問(wèn)題，提出了適用于任意磁場(chǎng)波形下的諧波分析法。

2Maxwell軟件鐵耗系數(shù)提取方法[8]

2.1單頻鐵耗系數(shù)提取

ANSYS Maxwell有限元仿真軟件，是將Bertotti提出的磁心鐵耗分離模型經(jīng)過(guò)等效進(jìn)行計(jì)算，且令式(2)中a=2，可得式(3)。

pv=pc+ph+pe=

K1Bm2+K2Bm1.5

(3)

式中渦流損耗pc=kc(fBm)2，磁滯損耗ph=khfBm2，附加損耗pe=ke(fBm)1.5，故可得K1=khf+kcf2，K2=kef1.5。

(1)K1、K2的計(jì)算。以式(4)作為擬合函數(shù)，采用最小二乘法對(duì)實(shí)測(cè)的損耗曲線進(jìn)行擬合分析，使得式(4)達(dá)到最小，則可得到K1、K2。

err(K1,K2)= ∑i[pvi-(K1Bmi2+

K2Bmi1.5)]2=min

(4)

式中：Pvi，Bmi——損耗特性曲線第i點(diǎn)對(duì)應(yīng)的數(shù)據(jù)。

(2)kc的計(jì)算。按式(5)進(jìn)行計(jì)算可得經(jīng)典渦流損耗系數(shù)kc

(5)

式中：σ——硅鋼片電導(dǎo)率；

d——硅鋼片厚度；

σFe——硅鋼片密度。

(3)kh、ke的計(jì)算。將K1、K2、kc和頻率f0，代入式(6)、(7)計(jì)算可得kh，ke。至此，完成了單頻條件下Maxwell軟件鐵耗系數(shù)的提取工作。

(6)

(7)

2.2多頻鐵耗系數(shù)提取

不考慮直流偏磁，在正弦磁通條件下，鐵耗可按修改后的Bertotti鐵耗分離模型進(jìn)行計(jì)算，如式(8)所示。

pv=ph+pc+pe=Khf(Bm)2+

Kc(fBm)2+Ke(fBm)1.5

(8)

式中：Bm——交流磁通分量的幅值；

f——頻率；

Kh——磁滯鐵耗系數(shù)；

Kc——渦流鐵耗系數(shù)；

Ke——附加鐵耗系數(shù)。

利用實(shí)測(cè)的鐵損曲線，采用最小二乘法對(duì)式(8)進(jìn)行擬合分析，得出磁滯鐵耗系數(shù)Kh， 渦流鐵耗系數(shù)Kc以及附加鐵耗系數(shù)Ke。

3實(shí)例計(jì)算與比較

本文以Maxwell自帶牌號(hào)為M19_29G的硅鋼片為例，對(duì)鐵耗系數(shù)進(jìn)行擬合求解，并與Maxwell提取出的鐵耗系數(shù)進(jìn)行對(duì)比分析。

3.1單頻鐵耗系數(shù)提取

(1)K1、K2的計(jì)算。本文采用MATLAB lsqcurvefit或者nlinfit命令對(duì)50Hz條件下?lián)p耗曲線進(jìn)行擬合分析，使得式(4)達(dá)到最小，得到K1=1.0025，K2=0.054301。50Hz條件下，損耗曲線如圖1所示。

圖1　50Hz條件下，鐵心損耗曲線

(2)kc的計(jì)算。將電導(dǎo)率σ=1960000S/m，硅鋼片厚度d=0.35e-3m，硅鋼片密度σFe=7872kg/m3代入式(5)，計(jì)算得經(jīng)典渦流損耗系數(shù)kc。

(3)kh、ke的計(jì)算。將K1、K2、kc和頻率f0，代入式(6)、(7)計(jì)算可得kh、ke，最終鐵耗系數(shù)的如表1所示。

表1　單頻條件下，鐵損系數(shù)

3.2多頻鐵耗系數(shù)提取

根據(jù)硅鋼片廠商提供的數(shù)據(jù)，在50Hz、100Hz、200Hz、400Hz、1000Hz時(shí)的M19_29G損耗曲線如圖2所示。

依據(jù)圖2曲線數(shù)據(jù)，采用最小二乘法對(duì)式(8)進(jìn)行擬合分析，得出磁滯鐵耗系數(shù)Kh， 渦流鐵耗系數(shù)Kc以及附加鐵耗系數(shù)Ke，如表2所示。

表2　多頻條件下，鐵損系數(shù)

從表1、表2中可以看出，本方法計(jì)算結(jié)果與Maxwell軟件計(jì)算結(jié)果完全吻合，證明了本文計(jì)算結(jié)果的正確性?，F(xiàn)有主流商用電磁場(chǎng)有限元仿真軟件是將現(xiàn)有的某些經(jīng)過(guò)修改后的模型應(yīng)用于仿真模型，雖然本文僅計(jì)算了Maxwell軟件所需的鐵耗系數(shù)，但對(duì)于其他沒(méi)有參數(shù)提取功能的商用有限元軟件，可參照本文所述方法，按照相應(yīng)軟件鐵耗計(jì)算模型修改擬合函數(shù)，并對(duì)現(xiàn)有的損耗曲線進(jìn)行擬合，得到所需參數(shù)，最終再進(jìn)行鐵耗計(jì)算，因此，本文所述鐵耗參數(shù)提取方法對(duì)其他商用有限元軟件的鐵耗系數(shù)提取同樣具有重要參考價(jià)值。

4結(jié)語(yǔ)

針對(duì)單頻和多頻兩種情況，本文利用第三方軟件對(duì)Maxwell軟件鐵耗計(jì)算所需的鐵耗系數(shù)進(jìn)行了擬合計(jì)算，得到了硅鋼片M19_29G的鐵損系數(shù)，并與Maxwell計(jì)算出的鐵耗系數(shù)完全吻合，證明了計(jì)算過(guò)程的正確性。本文對(duì)鐵耗系數(shù)求取過(guò)程是開放的，數(shù)據(jù)處理更加靈活，且求取出的鐵耗系數(shù)可直接輸入有限元計(jì)算軟件并對(duì)鐵耗進(jìn)行計(jì)算。同時(shí)，本文所述鐵耗參數(shù)提取方法對(duì)其他商用有限元軟件的鐵耗系數(shù)提取也具有重要參考價(jià)值。

【參 考 文 獻(xiàn)】

[1]陳世坤.電機(jī)設(shè)計(jì)[M].2版.北京： 機(jī)械工業(yè)社，2000.

[2]劉曉芳，趙海森，陳偉華，等.電機(jī)鐵耗的有限元計(jì)算方法研究進(jìn)展及有待解決的問(wèn)題.電機(jī)與控制應(yīng)用，2010，37(12)： 1-6.

[3]崔楊，胡虔生，黃允凱.任意頻率正弦波條件下鐵磁材料損耗的計(jì)算[J].微電機(jī)，2007，40(8)： 1-3.

[4]LIN D， ZHOU P， FU W N， et al. A dynamic core loss model for soft ferromagnetic and power ferrite materials in transient finite element analysis[J]. IEEE Transactions on Magnetics， 2004，40(2)： 1318-1321.

[5]BERTOTTI G. General properties of power losses in soft ferromagnetic material[J]. IEEE Transactions on Magnetics， 1998，24(1)： 621-630.

[6]黃平林，胡虔生，崔楊，等.PWM逆變器供電下電機(jī)鐵心損耗的解析計(jì)算[J].中國(guó)電機(jī)工程學(xué)報(bào)，2007，27(12)： 19-23.

[7]劉萬(wàn)太，彭曉，謝衛(wèi)才，等.變頻電機(jī)鐵心損耗計(jì)算方法的研究[J].電機(jī)與控制應(yīng)用，2010，37(11)： 11-14.

[8]ANSYS Company. Maxwell Online Help[G].2012.

Realization of the Method of Extracting Iron Loss Coefficient of Silicon Steel Sheet in ANSYS Maxwell*

GONGWenjun,HUANGShaogang,YANGYuwen,RAOLiqiang,ZENGYu

(College of Information Engineering, Nanchang University, Nanchang 330031, China)

Abstract：The iron loss calculation model of Maxwell software and method and steps of iron loss coefficient extraction were introduced. Iron loss curve of silicon steel sheet material M19_29G fitting was analyzed using the least square method. The iron loss coefficient which could be directly was applied to Maxwell could be obtained. The accuracy of the results was verified by using the Maxwell’s own extracted iron loss coefficient. Provided a new way for the realization of Maxwell software iron loss coefficient extraction, which made the data processing more open and flexible, and was no longer limited by the software itself, and the process of extracting iron loss coefficient for the Maxwell software was an important reference value for other commercial finite element software.

Key words：ANSYS Maxwell; least square method; core loss; core loss coefficient

收稿日期：2015-08-31

中圖分類號(hào)：TM 304

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

文章編號(hào)：1673-6540(2016)04- 0082- 04

作者簡(jiǎn)介：龔文軍(1991—)，男，碩士研究生，研究方向?yàn)殡姍C(jī)電磁場(chǎng)及電機(jī)控制。黃劭剛(1948—)，男，碩士，教授，研究方向?yàn)殡姍C(jī)電磁場(chǎng)及電機(jī)控制。楊玉文(1990—)，男，碩士研究生，研究方向?yàn)殡姍C(jī)電磁場(chǎng)分析。

*基金項(xiàng)目：江西省研究生創(chuàng)新專項(xiàng)資金資助(YC2014—S068)