唐旭， 王秀和， 徐定旺

(1.山東大學(xué) 電氣工程學(xué)院，山東 濟(jì)南 250061；2.青島理工大學(xué) 自動(dòng)化工程學(xué)院，山東 青島 266520)

異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算

唐旭1,2， 王秀和1， 徐定旺1

(1.山東大學(xué) 電氣工程學(xué)院，山東 濟(jì)南 250061；2.青島理工大學(xué) 自動(dòng)化工程學(xué)院，山東 青島 266520)

針對(duì)采用時(shí)步有限元法研究異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中永磁體工作點(diǎn)的變化時(shí)存在計(jì)算時(shí)間長(zhǎng)、不適合電機(jī)設(shè)計(jì)階段的快速計(jì)算的問(wèn)題，通過(guò)將異步起動(dòng)永磁同步電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型和磁路計(jì)算模型結(jié)合，建立了異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型。利用該解析計(jì)算模型計(jì)算了3臺(tái)樣機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的變化，分別得到了電機(jī)起動(dòng)過(guò)程中最大退磁磁場(chǎng)出現(xiàn)時(shí)的永磁體平均工作點(diǎn)和電機(jī)穩(wěn)定運(yùn)行時(shí)的永磁體平均工作點(diǎn)。通過(guò)與有限元法的計(jì)算結(jié)果作比較，驗(yàn)證了該解析計(jì)算模型的準(zhǔn)確性，可為該類(lèi)電機(jī)的快速設(shè)計(jì)提供參考。

異步起動(dòng)永磁同步電動(dòng)機(jī)；動(dòng)態(tài)數(shù)學(xué)模型；磁路計(jì)算模型；起動(dòng)過(guò)程；平均工作點(diǎn)

0 引 言

異步起動(dòng)永磁同步電動(dòng)機(jī)在結(jié)構(gòu)上可以看作是在感應(yīng)電動(dòng)機(jī)的轉(zhuǎn)子鐵心內(nèi)放置了永磁體，依靠轉(zhuǎn)子鼠籠產(chǎn)生的異步轉(zhuǎn)矩實(shí)現(xiàn)自起動(dòng)[1-2]，具有功率因數(shù)高、效率高、經(jīng)濟(jì)運(yùn)行范圍寬等優(yōu)點(diǎn)，在高效節(jié)能場(chǎng)合具有廣闊的應(yīng)用前景[3]。

與其他類(lèi)型的永磁電機(jī)一樣，異步起動(dòng)永磁同步電動(dòng)機(jī)在實(shí)際使用過(guò)程中，當(dāng)永磁體的工作點(diǎn)低于永磁體退磁曲線(xiàn)拐點(diǎn)時(shí)，永磁體會(huì)發(fā)生不可逆退磁[4-5]，這不僅會(huì)導(dǎo)致電機(jī)的反電動(dòng)勢(shì)減小、起動(dòng)電流增大[2]，還會(huì)導(dǎo)致電機(jī)的效率下降，并可能產(chǎn)生轉(zhuǎn)矩波動(dòng)及振動(dòng)噪聲[6-7]。

在異步起動(dòng)永磁同步電動(dòng)機(jī)的起動(dòng)過(guò)程中，定子繞組和轉(zhuǎn)子鼠籠導(dǎo)條中的電流會(huì)產(chǎn)生強(qiáng)退磁磁場(chǎng)，使永磁體的工作點(diǎn)不斷變化，可能導(dǎo)致永磁體發(fā)生不可逆退磁。異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中的永磁體退磁問(wèn)題得到了廣泛的研究：文獻(xiàn)[8]利用二維有限元法計(jì)算了電機(jī)起動(dòng)過(guò)程中永磁體某點(diǎn)磁密值的變化，并分析了永磁體退磁磁場(chǎng)產(chǎn)生的原因；文獻(xiàn)[9]和文獻(xiàn)[10]利用二維有限元法計(jì)算了電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的變化，并研究了電樞磁動(dòng)勢(shì)對(duì)永磁體的退磁作用以及負(fù)載條件、電源電壓等因素對(duì)永磁體退磁的影響。可以看出，現(xiàn)有文獻(xiàn)基本都采用時(shí)步有限元法研究異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中永磁體工作點(diǎn)的變化，尚未有文獻(xiàn)建立電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型。

在異步起動(dòng)永磁同步電動(dòng)機(jī)的設(shè)計(jì)階段，為了保證電機(jī)在實(shí)際運(yùn)行中永磁體不發(fā)生不可逆退磁，需要校核永磁體的最小平均工作點(diǎn)(永磁體遇到最大退磁磁場(chǎng)時(shí)的平均工作點(diǎn))，并使其不低于某一限定值。傳統(tǒng)的計(jì)算方法所基于的假設(shè)是：電機(jī)起動(dòng)至接近同步速時(shí)，轉(zhuǎn)子鼠籠對(duì)轉(zhuǎn)子內(nèi)永磁體的屏蔽作用減弱，電樞磁動(dòng)勢(shì)對(duì)永磁體的去磁作用最為嚴(yán)重[11]；但是并不能保證在電機(jī)的實(shí)際起動(dòng)過(guò)程中，永磁體的平均工作點(diǎn)不會(huì)低于該計(jì)算值。有限元法能夠計(jì)算永磁體各點(diǎn)的磁密分布，可以方便的計(jì)算電機(jī)起動(dòng)過(guò)程中永磁體最小工作點(diǎn)和平均工作點(diǎn)的變化；但是有限元法建模復(fù)雜、計(jì)算周期長(zhǎng)[12-13]，不適合電機(jī)設(shè)計(jì)階段的快速計(jì)算。

本文通過(guò)將異步起動(dòng)永磁同步電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型和磁路計(jì)算模型結(jié)合，建立了電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型，能夠快速計(jì)算并分析起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的變化，可為該類(lèi)電機(jī)的設(shè)計(jì)提供參考。本文以3臺(tái)異步起動(dòng)永磁同步電動(dòng)機(jī)為例，在Matlab/SIMULINK中建立了相應(yīng)的計(jì)算模型，計(jì)算并分析了電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的變化。利用有限元法驗(yàn)證了所建立解析計(jì)算模型的計(jì)算準(zhǔn)確性，并得到了相應(yīng)的結(jié)論。

1 起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型

1.1 異步起動(dòng)永磁同步電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型

本文建立的異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型主要包括電機(jī)的動(dòng)態(tài)數(shù)學(xué)模型和磁路計(jì)算模型兩部分。

三相靜止ABC坐標(biāo)系下的電機(jī)動(dòng)態(tài)數(shù)學(xué)模型十分復(fù)雜，通過(guò)坐標(biāo)變換(本文采用正交變換)將其變換到兩相旋轉(zhuǎn)dq0坐標(biāo)系下，可以使方程得以簡(jiǎn)化，便于計(jì)算分析。文獻(xiàn)[14-16]建立了異步起動(dòng)永磁同步電動(dòng)機(jī)在dq0坐標(biāo)系下的動(dòng)態(tài)數(shù)學(xué)模型，電壓方程為：

(1)

式中：d軸為轉(zhuǎn)子永磁磁場(chǎng)的軸線(xiàn)方向，q軸在空間超前d軸90°電角度；ud、uq分別為定子直、交軸繞組的電壓；u2d、u2q分別為轉(zhuǎn)子直、交軸繞組的電壓；R1為定子每相繞組的電阻；R2d、R2q分別為轉(zhuǎn)子直、交軸繞組的電阻；id、iq分別為定子直、交軸繞組的電流，i2d、i2q分別為轉(zhuǎn)子直、交軸繞組的電流；ψd、ψq分別為定子直、交軸繞組的磁鏈，ψ2d、ψ2q分別為轉(zhuǎn)子直、交軸繞組的磁鏈；ωr為電機(jī)轉(zhuǎn)子旋轉(zhuǎn)的電角速度。需要注意的是，上述轉(zhuǎn)子各量已歸算到定子側(cè)。

磁鏈方程為：

(2)

電磁轉(zhuǎn)矩方程為

Tem=p(ψdid-ψqiq)。

(3)

式中p為電機(jī)的極對(duì)數(shù)。

不計(jì)鐵耗和附加損耗時(shí)，機(jī)械運(yùn)動(dòng)方程為

(4)

式中：J為電機(jī)轉(zhuǎn)子和負(fù)載的轉(zhuǎn)動(dòng)慣量之和；TL為負(fù)載轉(zhuǎn)矩。

利用異步起動(dòng)永磁同步電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型，可以計(jì)算電機(jī)起動(dòng)過(guò)程中各時(shí)刻的定轉(zhuǎn)子直、交軸繞組電流等數(shù)據(jù)，為利用磁路計(jì)算模型計(jì)算永磁體的平均工作點(diǎn)奠定了基礎(chǔ)。

1.2 異步起動(dòng)永磁同步電動(dòng)機(jī)的磁路計(jì)算模型

永磁電機(jī)負(fù)載時(shí)的等效磁路計(jì)算模型如圖1所示[11]。其中，fc、fm、fd分別為永磁體磁動(dòng)勢(shì)、外磁路磁動(dòng)勢(shì)和等效增(去)磁磁動(dòng)勢(shì)的標(biāo)幺值，對(duì)于異步起動(dòng)永磁同步電動(dòng)機(jī)，該增(去)磁磁動(dòng)勢(shì)是定轉(zhuǎn)子繞組的直軸合成磁動(dòng)勢(shì)；λ0、λδ、λσ分別為永磁體內(nèi)磁導(dǎo)、主磁路磁導(dǎo)和漏磁路磁導(dǎo)的標(biāo)幺值；φm、φδ、φσ分別為永磁體、主磁路和漏磁路磁通的標(biāo)幺值。為了便于與電機(jī)動(dòng)態(tài)數(shù)學(xué)模型結(jié)合，本文規(guī)定，定轉(zhuǎn)子繞組直軸合成磁動(dòng)勢(shì)起增磁作用時(shí)為正。

求解等效磁路模型，可以得到永磁體的平均工作點(diǎn)為

(5)

在異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中，根據(jù)各時(shí)刻的定轉(zhuǎn)子繞組直軸合成磁動(dòng)勢(shì)，利用磁路計(jì)算模型，可以計(jì)算對(duì)應(yīng)時(shí)刻的永磁體平均工作點(diǎn)。

圖1 永磁電機(jī)的等效磁路模型Fig.1 Magnetic circuit of permanent magnet machine

1.3 定轉(zhuǎn)子直軸合成磁動(dòng)勢(shì)的計(jì)算

從上述分析可以看出，利用異步起動(dòng)永磁同步電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型可以計(jì)算出各時(shí)刻的定轉(zhuǎn)子繞組直軸電流，而利用磁路模型計(jì)算永磁體平均工作點(diǎn)需要的是定轉(zhuǎn)子直軸合成磁動(dòng)勢(shì)。因此，將電機(jī)的動(dòng)態(tài)數(shù)學(xué)模型和磁路計(jì)算模型結(jié)合起來(lái)計(jì)算電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的關(guān)鍵是確定定轉(zhuǎn)子繞組直軸電流與定轉(zhuǎn)子繞組直軸合成磁動(dòng)勢(shì)之間的關(guān)系。本文在多相繞組合成磁動(dòng)勢(shì)計(jì)算公式的基礎(chǔ)上，推導(dǎo)了采用正交變換時(shí)，定轉(zhuǎn)子繞組直軸電流與定轉(zhuǎn)子繞組直軸合成磁動(dòng)勢(shì)之間的關(guān)系式。

根據(jù)電機(jī)學(xué)基本原理，m相繞組的合成基波磁動(dòng)勢(shì)幅值為

(6)

式中：N為繞組每相串聯(lián)匝數(shù)；kw1為基波磁動(dòng)勢(shì)的繞組因數(shù)；Iφ為相電流的有效值。對(duì)應(yīng)的直軸磁動(dòng)勢(shì)為

(7)

式中θa為合成基波磁動(dòng)勢(shì)與電機(jī)直軸之間的夾角。

由于電機(jī)的直軸電流可以表示為

Id=Icosθa。

(8)

式中I為電流綜合矢量的幅值。因此，式(7)可以變換為

(9)

(10)

在電機(jī)動(dòng)態(tài)數(shù)學(xué)模型中，轉(zhuǎn)子側(cè)參數(shù)已經(jīng)歸算到定子側(cè)；因此，利用動(dòng)態(tài)數(shù)學(xué)模型得到的定轉(zhuǎn)子繞組的直軸電流可以計(jì)算定轉(zhuǎn)子繞組的直軸合成磁動(dòng)勢(shì)

(11)

作用于永磁體的等效增(去)磁磁動(dòng)勢(shì)為

Fd=Fd1Kad。

(12)

式中Kad為直軸電樞磁動(dòng)勢(shì)折算系數(shù)。

1.4 起動(dòng)過(guò)程中永磁體平均工作點(diǎn)計(jì)算模型的建立

利用上述定轉(zhuǎn)子繞組直軸電流與直軸合成磁動(dòng)勢(shì)的關(guān)系式，通過(guò)將異步起動(dòng)永磁同步電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型和磁路計(jì)算模型結(jié)合，可以建立電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型，計(jì)算流程如圖2所示。首先利用動(dòng)態(tài)數(shù)學(xué)模型計(jì)算的各時(shí)刻定轉(zhuǎn)子繞組直軸電流，由式(11)、式(12)計(jì)算出對(duì)應(yīng)的等效增(去)磁磁動(dòng)勢(shì)，然后利用磁路計(jì)算模型計(jì)算出對(duì)應(yīng)時(shí)刻的永磁體平均工作點(diǎn)。

圖2 計(jì)算流程圖Fig.2 Flowchart of the established model

2 起動(dòng)過(guò)程中永磁體平均工作點(diǎn)計(jì)算模型的仿真分析

2.1 仿真樣機(jī)的主要參數(shù)

本文在Matlab/SIMULINK中建立了異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型，以3臺(tái)異步起動(dòng)永磁同步電動(dòng)機(jī)為樣機(jī)，分別計(jì)算了電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的變化，并利用有限元法驗(yàn)證了所建立模型的計(jì)算準(zhǔn)確性。3臺(tái)樣機(jī)的主要設(shè)計(jì)參數(shù)和額定參數(shù)分別如表1～表6所示，其中，樣機(jī)的定轉(zhuǎn)子漏抗等參數(shù)是利用文獻(xiàn)[3]中的解析方法計(jì)算得到的。

表1 樣機(jī)一的主要設(shè)計(jì)參數(shù)Table 1 Design parameters of the first prototype motor

表2 樣機(jī)二的主要設(shè)計(jì)參數(shù)Table 2 Design parameters of the second prototype motor

表3 樣機(jī)三的主要設(shè)計(jì)參數(shù)Table 3 Design parameters of the third prototype motor

表4 樣機(jī)一的主要額定參數(shù)Table 4 Rated parameters of the first prototype motor

表5 樣機(jī)二的主要額定參數(shù)Table 5 Rated parameters of the second prototype motor

表6 樣機(jī)三的主要額定參數(shù)Table 6 Rated parameters of the third prototype motor

2.2 計(jì)算準(zhǔn)確性驗(yàn)證

2.2.1 永磁體最小平均工作點(diǎn)的計(jì)算

文獻(xiàn)[11]指出，在異步起動(dòng)永磁同步電動(dòng)機(jī)的起動(dòng)過(guò)程中，永磁體的最大退磁磁場(chǎng)出現(xiàn)在電機(jī)轉(zhuǎn)速接近同步速時(shí)。

在Matlab/SIMULINK中，利用本文建立的模型計(jì)算電機(jī)起動(dòng)過(guò)程中的永磁體最小平均工作點(diǎn)時(shí)，采取的方法是逐步增大電機(jī)的轉(zhuǎn)動(dòng)慣量直至電機(jī)不能牽入同步。對(duì)于樣機(jī)一，計(jì)算得到的電機(jī)轉(zhuǎn)速曲線(xiàn)和永磁體平均工作點(diǎn)曲線(xiàn)如圖3所示。可以看出，隨著電機(jī)轉(zhuǎn)速的上升，永磁體平均工作點(diǎn)曲線(xiàn)不斷波動(dòng)，退磁磁場(chǎng)強(qiáng)度不斷變化；由于電機(jī)的轉(zhuǎn)動(dòng)慣量過(guò)大，電機(jī)未牽入同步。最終，電機(jī)轉(zhuǎn)速低于同步速小幅波動(dòng)，永磁體的平均工作點(diǎn)曲線(xiàn)也隨著電機(jī)轉(zhuǎn)速的波動(dòng)而波動(dòng)，且波動(dòng)幅值保持不變，最小值維持在0.13 T，對(duì)應(yīng)的電機(jī)轉(zhuǎn)速為988 r/min，此時(shí)的永磁體退磁磁場(chǎng)強(qiáng)度達(dá)到最大值。在該狀態(tài)下，定轉(zhuǎn)子繞組合成磁動(dòng)勢(shì)與永磁磁場(chǎng)軸線(xiàn)的相對(duì)位置周期性變化，二者之間的夾角可以是0～360°電角度之間的任意角度，因此可以找到永磁體的最小平均工作點(diǎn)。對(duì)于樣機(jī)一，利用本文建立的模型計(jì)算得到的永磁體最小平均工作點(diǎn)為0.13 T。

圖3 轉(zhuǎn)動(dòng)慣量過(guò)大導(dǎo)致樣機(jī)一不能牽入同步時(shí)的轉(zhuǎn)速曲線(xiàn)與永磁體平均工作點(diǎn)曲線(xiàn)Fig.3 Speed and magnets′ average operating point curves during starting failure caused by large moment of inertia

在利用有限元法計(jì)算樣機(jī)起動(dòng)過(guò)程中的永磁體最小平均工作點(diǎn)時(shí)，采取相似的方法，不斷增大電機(jī)的轉(zhuǎn)動(dòng)慣量直至電機(jī)不能牽入同步，可以得到永磁體的最小平均工作點(diǎn)。對(duì)于樣機(jī)一，有限元法計(jì)算得到的永磁體最小平均工作點(diǎn)為0.11 T。

對(duì)于樣機(jī)二和樣機(jī)三，采用上述方法，利用本文建立的解析計(jì)算模型和有限元法可以分別計(jì)算其最小平均工作點(diǎn)。表7比較了利用本文建立的解析計(jì)算模型和有限元法計(jì)算得到的3臺(tái)樣機(jī)的最小平均工作點(diǎn)，可以看出，本文所建立的解析計(jì)算模型可以較準(zhǔn)確地計(jì)算異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中的最小平均工作點(diǎn)。在電機(jī)設(shè)計(jì)階段，本文所建立的解析計(jì)算模型可用于校驗(yàn)永磁體的最小平均工作點(diǎn)。

表7 樣機(jī)的最小平均工作點(diǎn)比較

2.2.2 穩(wěn)態(tài)時(shí)永磁體平均工作點(diǎn)的計(jì)算

利用Matlab/SIMULINK中建立的模型計(jì)算得到了樣機(jī)一帶額定負(fù)載轉(zhuǎn)矩、以轉(zhuǎn)子本身轉(zhuǎn)動(dòng)慣量起動(dòng)時(shí)的轉(zhuǎn)速曲線(xiàn)和永磁體平均工作點(diǎn)曲線(xiàn)，如圖4所示?？梢钥闯?，隨著電機(jī)轉(zhuǎn)速的不斷升高，永磁體平均工作點(diǎn)曲線(xiàn)多次波動(dòng)；在a點(diǎn)(0.067 s)處，永磁體的平均工作點(diǎn)曲線(xiàn)出現(xiàn)最小值(0.169 T)，對(duì)應(yīng)的電機(jī)轉(zhuǎn)速為505 r/min；電機(jī)牽入同步穩(wěn)定運(yùn)行后，永磁體的平均工作點(diǎn)維持在0.742 T。

圖4 樣機(jī)一起動(dòng)時(shí)的轉(zhuǎn)速曲線(xiàn)與永磁體平均工作點(diǎn)曲線(xiàn)Fig.4 Speed and magnets′ average operating point curves during starting process of the first prototype motor

利用本文建立的解析計(jì)算模型分別計(jì)算了3臺(tái)樣機(jī)帶不同負(fù)載轉(zhuǎn)矩起動(dòng)至穩(wěn)定運(yùn)行后的永磁體平均工作點(diǎn)，并與有限元法的計(jì)算結(jié)果作比較，結(jié)果如圖5～圖7所示?？梢钥闯?，解析計(jì)算結(jié)果與有限元計(jì)算結(jié)果吻合較好，本文所建立的解析計(jì)算模型可以準(zhǔn)確地計(jì)算異步起動(dòng)永磁同步電動(dòng)機(jī)穩(wěn)定運(yùn)行時(shí)的永磁體平均工作點(diǎn)。

圖5 樣機(jī)一穩(wěn)定運(yùn)行時(shí)的永磁體平均工作點(diǎn)Fig.5 Magnets′ average operating points when the first prototype motor operates stably

圖6 樣機(jī)二穩(wěn)定運(yùn)行時(shí)的永磁體平均工作點(diǎn)Fig.6 Magnets′ average operating points when the second prototype motor operates stably

圖7 樣機(jī)三穩(wěn)定運(yùn)行時(shí)的永磁體平均工作點(diǎn)Fig.7 Magnets′ average operating points when the third prototype motor operates stably

2 結(jié) 論

本文通過(guò)將異步起動(dòng)永磁同步電動(dòng)機(jī)的動(dòng)態(tài)數(shù)學(xué)模型與磁路計(jì)算模型結(jié)合，建立了電機(jī)起動(dòng)過(guò)程中永磁體平均工作點(diǎn)的解析計(jì)算模型，并通過(guò)與有限元法的計(jì)算結(jié)果作比較驗(yàn)證了所建立模型的計(jì)算準(zhǔn)確性。

結(jié)果表明，本文所建立的計(jì)算模型可以較準(zhǔn)確地計(jì)算異步起動(dòng)永磁同步電動(dòng)機(jī)起動(dòng)過(guò)程中的最小平均工作點(diǎn)和電機(jī)穩(wěn)定運(yùn)行時(shí)的永磁體平均工作點(diǎn)；該模型可為異步起動(dòng)永磁同步電動(dòng)機(jī)設(shè)計(jì)階段的快速計(jì)算提供參考。

[1] KIM T H,HONG J P.A study on the irreversible magnet demagnetization in single-phase line-start permanent magnet motor[J].Journal of Applied Physics,2009,105(7): 07F108-07F108-3.

[2] KANG G H,HUR J,NAM H,et al.Analysis of irreversible magnet demagnetization in line-start motors based on the finite-element method[J].IEEE Transactions on Magnetics,2003,39(3): 1488-1491.

[3] 王秀和,楊玉波,朱常青.異步起動(dòng)永磁同步電動(dòng)機(jī)-理論、設(shè)計(jì)與測(cè)試[M].北京: 機(jī)械工業(yè)出版社,2009: 1-6.

[4] HONG J,HYUN D,LEE S B,et al.Automated monitoring of magnet quality for permanent magnet synchronous motors at standstill[J].IEEE Transactions on Industry Applications,2010,46(4): 1397-1405.

[5] RIBA RUIZ J R,ROSERO J A,ESPINOSA A G,et al.Detection of demagnetization faults in permanent-magnet synchronous motors under non-stationary conditions[J].IEEE Transactions on Magnetics,2009,45(7): 2961-2969.

[6] URRESTY J C,ATASHKHOOEI R,RIBA J,et al.Shaft trajectory analysis in a partially demagnetized permanent-magnet synchronous motor[J].IEEE Transactions on Industrial Electronics,2013,60 (8): 3454-3461.

[7] TORREGROSSA D,KHOOBROO A,FAHIMI B.Prediction of acoustic noise and torque pulsation in PM synchronous machines with static eccentricity and partial demagnetization using field reconstruction method[J].IEEE Transactions on Industrial Electronics,2012,59(2): 934-944.

[8] 盧偉甫,劉明基,羅應(yīng)立,等.自起動(dòng)永磁同步電機(jī)起動(dòng)過(guò)程退磁磁場(chǎng)的計(jì)算與分析[J].中國(guó)電機(jī)工程學(xué)報(bào),2011,31(15): 53-60. LU Weifu,LIU Mingji,LUO Yingli,et al.Demagnetization field analysis and calculation for line-start permanent magnet synchronous motor during start process[J].Proceedings of the CSEE,2011,31(15): 53-60.

[9] 盧偉甫,羅應(yīng)立,趙海森.自起動(dòng)永磁同步電機(jī)起動(dòng)過(guò)程電樞反應(yīng)退磁分析[J].電機(jī)與控制學(xué)報(bào),2012,16(7): 29-33. LU Weifu,LUO Yingli,ZHAO Haisen.Armature reaction demagnetization analysis for line-start permanent magnet synchronous motor during start process[J].Proceedings of Electric Machines and Control,2012,16(7): 29-33.

[10] LU Weifu,LIU Mingji,LUO Yingli,et al.Influencing factors on the demagnetization of line-start permanent magnet synchronous motor during its starting process[C]//International Conference on Electrical Machines and Systems,August 20-23,2011,Beijing,China.2011: 1-4.

[11] 唐任遠(yuǎn).現(xiàn)代永磁電機(jī)[M].北京: 機(jī)械工業(yè)出版社,2011: 185-187.

[12] FAROOQ J,SRAIRI S,DJERDIR A,et al.Use of permeance network method in the demagnetization phenomenon modeling in a permanent magnet motor[J].IEEE Transactions on Magnetics,2006,42(4): 1295-1298.

[13] WANG J,WANG W,ATALLAH K,et al.Demagnetization assessment for three-phase tubular brushless permanent-magnet machines[J].IEEE Transactions on Magnetics,2006,44(9): 2195-2203.

[14] 李國(guó)建,王道涵,朱常青.基于MATLAB/ SIMULINK的三相異步起動(dòng)永磁同步電動(dòng)機(jī)的建模與仿真[J].山東大學(xué)學(xué)報(bào)(工學(xué)版),2006,36(5): 40-43. LI Guojian,WANG Daohan,ZHU Changqing.Modeling and simulation of three-phase line-start permanent magnet synchronous motor[J].Journal of Shandong University (Engineering Science),2006,36(5): 40-43.

[15] 劉仲恕.永磁同步電動(dòng)機(jī)狀態(tài)方程研究和起動(dòng)過(guò)程仿真[J].福建工程學(xué)院學(xué)報(bào),2003,1(2): 43-46. LIU Zhongshu.On the synchronous motor state equation of permanent magnetism and emulation of starting course[J].Journal of Fujian University of Technology,2003,1(2): 43-46.

[16] WANG D H,WANG X H,CHEN H Y,et al.Matlab/Simulink-based simulation of line-start PMSM used in pump jacks[C]// IEEE Conference on Industrial Electronics and Applications,May 23-25,2007,Harbin,China.2007: 1179-1181.

[17] 王秀和,孫雨萍.電機(jī)學(xué)[M].北京: 機(jī)械工業(yè)出版社,2009: 322-334.

Analytical calculation of permanent magnets′ average operating point for line-start permanent magnet synchronous motor during starting process

TANG Xu1,2， WANG Xiu-he1， XU Ding-wang1

(1.School of Electrical Engineering,Shandong University,Jinan 250061,China;2.School of Automation Engineering,Qingdao Technological University,Qingdao 266520,China)

When finite element method (FEM) is adopted to study the variation of magnets′ operating point in line-start permanent magnet synchronous motor (LSPMSM) during the starting process,the computation is time consuming.Thus FEM is not suitable for the fast calculation at the motor design stage.To solve this problem,the analytical model to calculate the variation of magnets′ average operating point during the starting process of LSPMSM was established by combining the dynamic mathematical model and magnetic circuit model of LSPMSM.The variations of magnets′ average operating point during the starting process of three prototype motors were calculated by the above analytical model.And magnets’ average operating point when the severest demagnetization field during the starting process occurs as well as magnets’ average operating point when the motor operates stably were obtained.Accuracy of the analytical calculation model was validated by comparing with the result of FEM.Thus,the established analytical model can benefit fast design of LSPMSM.

line-start permanent magnet synchronous motor; dynamic mathematical model; magnetic circuit model; starting process; average operating point

2014-06-17

國(guó)家自然科學(xué)基金(51577107)

唐 旭(1988—)，男，博士研究生，研究方向?yàn)橛来烹姍C(jī)； 王秀和(1967—)，男，博士，教授，博士生導(dǎo)師，研究方向?yàn)橛来烹姍C(jī)； 徐定旺(1992—)，男，碩士研究生，研究方向?yàn)橛来烹姍C(jī)。

王秀和

10.15938/j.emc.2017.05.002

TM 351

A

1007-449X(2017)05-0008-07