周 鵬，徐衍亮，陳國偉，王文凱

(1.山東大學 電氣工程學院，濟南 250061；2.國家電網(wǎng)內(nèi)蒙古東部電力綜合能源服務有限公司，呼和浩特 010010)

0 引 言

永磁同步電機因其功率密度、轉(zhuǎn)矩密度大，功率因數(shù)高，具有良好的調(diào)速性能，且在寬廣的轉(zhuǎn)速范圍內(nèi)具有較高工作效率等優(yōu)勢，在工程機械、家電、交通以及軍工國防等領域都得到了十分廣泛的應用[1]。永磁電機因裝配有永磁體，因此當轉(zhuǎn)子轉(zhuǎn)動時，繞組會產(chǎn)生空載反電動勢?？蛰d反電動勢的大小是永磁電機的一項核心設計指標，對電機的運行性能有重要影響作用。實際應用中，永磁電機在額定轉(zhuǎn)速下的空載反電動勢并不是恒定不變的，隨著電機投入工作年限的增加，由于工作環(huán)境的溫度反復變化、長期的機械應力或者可能存在的各類化學反應，都會導致永磁體的磁性能發(fā)生一定程度的下降[2-5]，進而影響電機的反電動勢系數(shù)。反電動勢系數(shù)的下降會造成電機功率因數(shù)降低、有功電流增加、電機發(fā)熱、效率降低等各種問題。因此，對于長期運行于高負荷或是惡劣工況下的永磁電機，我們一般需要定期對其反電動勢進行測量校驗，以確保電機中永磁體磁性能的穩(wěn)定性。

對于永磁同步電機反電動勢的測量，傳統(tǒng)的方法主要有兩種：一是用同極永磁電機或同步電機拖動待測永磁電機運行于額定轉(zhuǎn)速下，檢測待測試永磁電機繞組的開路電壓，該電壓即為額定空載反電動勢。二是使待測樣機空載運行，調(diào)整其輸入電壓，當電樞電流達到最小值時，此時可近似將測得的輸入電壓視為永磁電動機空載反電動勢。這兩種測量方式都屬于離線測量，具有兩點不足之處：一是永磁電動機空載運行所需要的對拖實驗需要在專門的測試平臺上進行，這勢必會中斷設備生產(chǎn)，耗費大量的人力物力，操作極不方便，同時也不便于經(jīng)常性的對永磁電機磁性能進行檢查；二是離線測試的方法只能夠檢驗永磁體的靜態(tài)磁性能，而無法測試永磁體在不同運行工況下的動態(tài)磁性能，提前甄別某些可能存在的隱患，比如永磁體因轉(zhuǎn)子發(fā)熱所引起的可逆退磁現(xiàn)象。電機在實際運行中，永磁體會因為渦流損耗發(fā)熱、轉(zhuǎn)子不易散熱等原因出現(xiàn)可逆退磁現(xiàn)象，而隨著電機卸去負載，電機的磁性能又會逐漸恢復。永磁體可逆退磁的產(chǎn)生是由于電機電磁結(jié)構(gòu)設計的不合理、變頻器輸入電流的復雜時間諧波，以及轉(zhuǎn)子散熱設計不足等多種原因綜合所導致的[6-8]。處于可逆退磁工況下的永磁電機，相比于理想的工況，其損耗和發(fā)熱會增加[9-11]，效率和輸出轉(zhuǎn)矩會下降，極大的影響電機的工作性能和使用壽命。電機若長期工作于可逆退磁工況下，永磁體的工作點會隨溫升的正反饋而不斷下移，最終將導致永磁體發(fā)生不可逆退磁，危害設備及人員安全。

針對傳統(tǒng)離線測量永磁電機反電動勢的缺陷，本文提出了一種能夠在線測試永磁同步電機反電動勢的方法，通過在運行過程中瞬時通斷電機的驅(qū)動電路或者利用電機斷電后的自由減速過程，對電機接線端子輸出的電壓波形進行捕獲，通過一定的計算方法反推出永磁電機在額定轉(zhuǎn)速下的空載反電動勢。文中對該計算方法做出了理論分析，設計了核心部分的原理電路及測試樣機，并通過與示波器的對比實驗印證了理論和實踐的可行性和精確性。

1 反電動勢的測試原理

永磁電機在運行時，電樞繞組產(chǎn)生的反電動勢可由式(1)得到。

E0=pω0ψ0

(1)

式中,E0為永磁同步電機空載反電動勢；p為磁極對數(shù)；ω0為機械角速度；ψ0為穿過繞組的磁鏈，當電機空載運行時，該磁鏈僅為永磁磁鏈。

由式(1)可知，電機的空載反電動勢大小與轉(zhuǎn)子轉(zhuǎn)速、磁極對數(shù)和磁鏈成正比關系。對于一臺設計好的永磁同步電機，其磁極對數(shù)和永磁磁鏈為均已確定。因此，該永磁電機的空載反電動勢大小正比于轉(zhuǎn)子的轉(zhuǎn)速，如式(2)所示。

(2)

該式表明，一臺永磁電機如果已知其在某一轉(zhuǎn)速下的反電動勢，可根據(jù)轉(zhuǎn)速比來推算它在其他轉(zhuǎn)速下的反電動勢。

當永磁同步電機工作于恒速狀態(tài)下時，在一個電周期內(nèi)，穿過一相繞組的磁鏈正弦變化，則其反電動勢波形正弦變化。而當永磁同步電機工作于變速的狀態(tài)下時，比如電機斷電后轉(zhuǎn)子在阻轉(zhuǎn)矩的作用下發(fā)生自由運動的過程，此時反電動勢波形不再是一個標準的正弦波。任一時刻，波形的幅值和頻率都隨著轉(zhuǎn)子轉(zhuǎn)速的變化而變化。圖1為某電機斷電前后幾個電周期內(nèi)的電壓實測波形。可以看到，斷電后的波形雖形似一個正弦波(文中簡稱該波形為“類正弦波”)，但在任一時刻下，波形的幅值和頻率都是轉(zhuǎn)子轉(zhuǎn)速的函數(shù)，隨著轉(zhuǎn)速的下降，幅值和頻率也在同步的衰減。

圖1 永磁電機斷電瞬間接線端子實測電壓變化

轉(zhuǎn)速的下降規(guī)律因負載類型的不同而不同。在下一小節(jié)，本文將根據(jù)不同的負載性質(zhì)研究電機斷電后反電動勢波形的衰減規(guī)律，并對比各負載類型下波形正弦度的優(yōu)劣，對正弦度最差的一類反電動勢波形進行重點研究。針對該波形，本文將提出一種波形等效的方法，可以將衰減的類正弦波等效成一個標準的正弦波。利用等效正弦波的波形頻率和有效值，通過式(2)進行頻率歸算，進而得到電機在額定轉(zhuǎn)速下的空載反電動勢有效值。

2 斷電減速過程的動力學特點

如上所述，反電動勢波形的衰減規(guī)律是與電機的負載性質(zhì)息息相關的。根據(jù)所拖動的負載性質(zhì)的不同，可將負載分為恒功率負載、恒轉(zhuǎn)矩負載和流體性負載。電機穩(wěn)態(tài)運行時，滿足轉(zhuǎn)矩平衡方程如式(3)所示。

Te=T0+T2

(3)

當變頻器切斷供電時，電磁轉(zhuǎn)矩Te=0。此時，施加在轉(zhuǎn)軸上的轉(zhuǎn)矩僅有阻轉(zhuǎn)矩。對于電機，阻轉(zhuǎn)矩類型包括摩擦阻轉(zhuǎn)矩、流體阻轉(zhuǎn)矩和磁阻轉(zhuǎn)矩。在一個電周期內(nèi)，磁阻轉(zhuǎn)矩合力為0，分析時暫不考慮。在其余兩項中，摩擦阻轉(zhuǎn)矩為恒轉(zhuǎn)矩負載，流體阻轉(zhuǎn)矩為流體類負載。斷電后電機轉(zhuǎn)矩的動態(tài)平衡方程式可由式(4)表示。

(4)

式中，J為電機轉(zhuǎn)軸的轉(zhuǎn)動慣量，α(t)為轉(zhuǎn)軸的角加速度，ω0為轉(zhuǎn)子的初速度，ω(t)為任意時刻下轉(zhuǎn)子的角速度。由于阻力轉(zhuǎn)矩T2類型的差異，求解方程得到轉(zhuǎn)子的運動衰減過程是不同的。設電機分別在摩擦阻轉(zhuǎn)矩和流體阻轉(zhuǎn)矩的作用下，于t時刻內(nèi)由轉(zhuǎn)速n下降至n/2，轉(zhuǎn)子轉(zhuǎn)速的衰減規(guī)律對比如圖2所示。

圖2 不同負載類型下轉(zhuǎn)子的速度衰減特征

依照轉(zhuǎn)速的變化規(guī)律，圖3為同一永磁電機對應兩種阻轉(zhuǎn)矩下，反電動勢的衰減變化規(guī)律。

圖3 反電動勢衰減特征

為了比較不同阻轉(zhuǎn)矩類型下反電動勢波形的正弦程度差異，定義波形的正弦度如下

(5)

式中,Amp1和Amp2分別為反電動勢波形在一個電周期內(nèi)的最大值和最小值的絕對值，而Cyc則表示一個電周期內(nèi)正半周波或負半周波的時間，兩者中數(shù)值大者記為Cyc1，數(shù)值較小者記為Cyc2，參數(shù)標注如圖4所示。

圖4 參數(shù)標注

由的定義可知k≥1，當k=1時為正弦波。該值越接近1，波形的正弦度越好；反之正弦度越差。

在兩種阻轉(zhuǎn)矩下，轉(zhuǎn)子降速過程中計算空載反電動勢波形序列的正弦度變化趨勢如圖5所示，其中第n個波形序列是指電機減速過程中第n個電周期內(nèi)的反電動勢波形?？梢娔Σ令愖柁D(zhuǎn)矩的正弦度為基本保持恒定值，而當阻轉(zhuǎn)矩為流體類負載，且轉(zhuǎn)速越接近斷電自由運動的初始速度時，反電動勢的波形正弦度最差。因此對此條件下的反電動勢波形進行等效正弦波研究最具代表性。

圖5 不同阻轉(zhuǎn)矩下波形正弦度的變化趨勢

3 反電動勢的等效和推算

電機空載運行時，A相繞組的電壓方程為

(6)

A相空載反電動勢與A相磁鏈的波形和相位關系如圖6所示。

圖6 反電動勢與磁鏈

在任意的轉(zhuǎn)子運動狀態(tài)下，當轉(zhuǎn)子轉(zhuǎn)過一對磁極，此時反電動勢可視為一個完整的電周期。對一個電周期內(nèi)反電動勢E0的絕對值計算積分:

(7)

式中，T為一個電周期波形的總時間，0～t1之間為正半周波，t1～t2之間為負半周波。由上式可知，在一個電周期內(nèi)，無論波形是否正弦，均可將積分區(qū)間按E0的正負分為兩部分。每一部分對時間的積分都是穿過A相繞組的永磁磁鏈最大值的兩倍。因此，一個電周期內(nèi)反電動勢E0絕對值的積分是一個常數(shù)，等于穿過A相繞組永磁磁鏈最大值的四倍。

將式(7)等式兩端同時除以電周期T，則得到該反電動勢波形的平均值。若T保持不變，則波形的平均值為一常數(shù)。

(8)

由以上推導可知，雖然波形的形狀隨著轉(zhuǎn)速變化存在任意性，但波形在一個電周期內(nèi)的平均值具有確定性。因此，可以將這個任意的波形等效為一個平均值和周期與其相同的標準正弦波，如圖7所示。

圖7 任意波形與同周期的等效正弦波

基于以上思想，就可以設計相關程序，實現(xiàn)在任意轉(zhuǎn)速和任意運動狀態(tài)下反電動勢的快速測量。只要能夠得到電機在到一個電周期內(nèi)的反電動勢波形，就可將其等效為一個頻率和有效值已知的正弦波，并利用式(2)計算電機在額定轉(zhuǎn)速下的空載反電動勢。

以一臺48極216槽、50 Hz額定頻率下轉(zhuǎn)速125 r/min、額定功率1800 kW的泵用低速直驅(qū)永磁同步電機為例進行仿真以驗證上述分析。建立2D有限元模型，如圖8所示。額定轉(zhuǎn)速下，仿真得到空載線反電動勢有效值為5066.2 V。

由于主要承災體未產(chǎn)生變化，邊坡加固后定量風險評估計算參量災害到達承災體概率PT∶L、承災體時空概率PS∶T、承災體易損性VD∶T、Vprop參照上文計算獲得的數(shù)據(jù)(見表3)。

圖8 2D FEM模型及空載反電動勢

泵站電機的阻轉(zhuǎn)矩為流體型負載，設電機在額定轉(zhuǎn)速125 r/min下突然斷電，1 s內(nèi)轉(zhuǎn)速下降至50%，仿真得到減速過程中的反電動勢波形序列，取序列中的第一個電周期內(nèi)的反電動勢波，對其進行傅里葉級數(shù)擬合，并計算得到擬合波形平均值為4.4154 kV，電周期為0.02066 s。

如上所述，將這個類正弦波波形等效為平均值為4.4154 kV，周期為0.02066 s的標準正弦波。同時得到等效反電動勢波的幅值為6.9357 kV，有效值為4.9043 kV。通過式(2)可求得額定轉(zhuǎn)速下有效值E0’為5.0661 kV，與直接進行額定轉(zhuǎn)速下空載反電動勢仿真得到的有效值5.0662 kV相比，計算誤差僅為0.023‰。

進一步的，將同一電機模型在125 r/min，100 r/min和75 r/min下進行3組對比實驗，計算結(jié)果和誤差如下表所示?？梢娛褂迷撚嬎惴椒ㄔ谌我廪D(zhuǎn)速下反推電機的額定空載反電動勢，都能保持計算結(jié)果的一致和準確。

表1 仿真計算結(jié)果對比

4 硬件與程序設計

根據(jù)上述計算原理，基于STM32F4系列主控單元和AD7606系列模數(shù)轉(zhuǎn)換芯片，設計樣機如下。該設備具有針對電機的在線電壓測量和斷電自動計算反電動勢的功能，計算速度快，實時性好，反電動勢從采集到完成計算所需時間小于70 ms，且能夠顯示出采集的波形。

程序的核心是對特定頻率電動勢波形的捕獲和反推算法的實現(xiàn)，簡要的程序流程如圖10所示。

圖9 設備樣機

圖10 簡要程序流程圖

電機反電動勢測試儀的硬件功能框圖如圖11所示。其中AD模塊負責電壓信號的模數(shù)轉(zhuǎn)換；過零檢測模塊用于頻率計算、以及對電動勢波形的電周期起止點進行判斷。

圖11 設備硬件功能框圖

5 實驗驗證

搭建實驗平臺如圖12所示，共進行兩組實驗。第一組實驗為對拖實驗，被試電機為一臺4極永磁同步電機，額定轉(zhuǎn)速1500 r/min。將被試永磁電機拖動至額定轉(zhuǎn)速，用示波器記錄空載反電動勢的波形，并計算有效值，額定轉(zhuǎn)速下的空載反電動勢波形如圖13所示。第二組實驗為單機實驗，使用變頻器驅(qū)動被試電機運行至額定轉(zhuǎn)速后切斷供電，并使用樣機測試儀計算電機的反電動勢有效值。

圖12 實驗測試平臺

圖13 額定轉(zhuǎn)速下的空載反電動勢波形

一共設置3個不同的初始速度，分別比較設備在不同轉(zhuǎn)速下反推額定轉(zhuǎn)速空載反電動勢的準確性，相關計算結(jié)果如表2所示，可見該測試設備在不同轉(zhuǎn)速下反推出的額定轉(zhuǎn)速空載反電動勢結(jié)果一致。對比仿真結(jié)果，計算誤差略微增大，主要原因來自于電壓采集電路中各元器件的精度和測量電路的測試精度的限制。

表2 實驗數(shù)據(jù)對比

6 結(jié) 論

永磁電機的運行性能與電機中永磁體的磁性能具有密切的關系，因此常需要測量永磁電機的空載反電動勢，以檢驗磁性能的穩(wěn)定性。本文針對永磁電機離線測量反電動勢所存在的固有缺陷，提出一種在線反電動勢測量方法。該方法的理論依據(jù)是：當永磁電機空載運行時，一個電周期內(nèi)，穿過一相繞組的磁鏈守恒?；诖耍梢詫⒁粋€衰減的反電動勢波形等效為一個標準正弦波，并利用等效波形的有效值和頻率計算電機在額定轉(zhuǎn)速下的反電動勢。

因轉(zhuǎn)子在流體阻轉(zhuǎn)矩下反電動勢的波形正弦性最差，因此以一臺48極、額定功率1800 kW的泵用永磁同步電機有限元仿真模型為例，得到電機在斷電減速過程中的衰減反電動勢波形序列。選取該序列中的第一個電周期的波形，利用等效原則得到對應的正弦反電動勢波，并計算得出額定轉(zhuǎn)速下的空載反電動勢，該數(shù)值與直接仿真得到的空載反電動勢數(shù)值基本一致。

以一臺4極、額定轉(zhuǎn)速1500 r/min的小型永磁電機為例，搭建實驗平臺，進行了多組不同轉(zhuǎn)速下的實測對比，采樣樣機的在線計算結(jié)果和采用對拖實驗得到的反電動勢數(shù)值也基本吻合。

綜上，本文所提出的反電動勢測量方法和開發(fā)的樣機設備具有一定程度的理論意義和實踐意義。