呂躍剛 張文原 蔡騰潛

(1.華北電力大學(xué)控制與計(jì)算機(jī)工程學(xué)院，北京 102206；2.國網(wǎng)冀北電力有限公司秦皇島供電公司，河北 秦皇島 066000)

鼠籠型異步電機(jī)廣泛應(yīng)用于工業(yè)生產(chǎn)，但該類電機(jī)啟停頻繁、運(yùn)行時(shí)間長、故障多發(fā)。轉(zhuǎn)子斷條是最常見的故障之一，約占全部故障的10%。故障初期電機(jī)效率降低，振動增加，嚴(yán)重時(shí)導(dǎo)條連續(xù)斷裂，甚至有可能損壞定子結(jié)構(gòu)，致使電機(jī)報(bào)廢[1,2]。學(xué)者們針對斷條故障提出了多種診斷方法。其中，振動頻譜技術(shù)憑借其原理簡單及起步較早等優(yōu)勢，在國內(nèi)得到了最為廣泛的應(yīng)用。然而這種技術(shù)硬件成本高、信號干擾大，經(jīng)濟(jì)性和有效性也受到一定的制約。近年來，定子電流分析技術(shù)(MCSA)已經(jīng)引起國內(nèi)的重視，如方瑞明等利用提取定子電流特定頻率成分的方法診斷轉(zhuǎn)子故障[3]。然而實(shí)際情況下，定子電流信號中被廣泛認(rèn)可的故障特征頻率(1±2ks)f1成分與基頻成分相距很近，頻率分辨率低及頻譜泄漏等原因?qū)е缕錁O易被湮沒。而在許多以DSP或單片機(jī)為硬件核心的在線診斷系統(tǒng)中，數(shù)據(jù)存儲容量有限，難以通過增長采樣時(shí)間的方式提高分辨率。這一特點(diǎn)限制了定子電流分析技術(shù)的發(fā)展，同時(shí)也是一個(gè)有待解決的難題。

轉(zhuǎn)速是電機(jī)的重要運(yùn)行參數(shù)，測量技術(shù)成熟，成本低廉且準(zhǔn)確性高、響應(yīng)快，目前已被應(yīng)用于內(nèi)燃機(jī)等旋轉(zhuǎn)機(jī)械故障診斷中[4]。筆者以異步電機(jī)瞬時(shí)轉(zhuǎn)速為分析對象，通過對瞬時(shí)轉(zhuǎn)速的時(shí)域、頻域分析，尋找其與斷條故障之間的聯(lián)系。

筆者依次分析了轉(zhuǎn)子斷條故障發(fā)生后電機(jī)內(nèi)部磁場、轉(zhuǎn)矩及電流等參量的變化情況，最終得出故障在瞬時(shí)轉(zhuǎn)速中的表現(xiàn)。利用Maxwell軟件對故障電機(jī)進(jìn)行有限元建模，分別對正常狀態(tài)下和不同程度故障狀態(tài)下的運(yùn)行過程進(jìn)行仿真，實(shí)驗(yàn)表明轉(zhuǎn)速信號中的2sf1成分能良好地反映故障情況，易于識別。此外，轉(zhuǎn)速波動程度、平均轉(zhuǎn)速值還能有效反映故障的嚴(yán)重程度。與傳統(tǒng)的定子電流分析方法進(jìn)行對比發(fā)現(xiàn)，轉(zhuǎn)速分析方法簡單易行、準(zhǔn)確性高，而且可以降低信號采樣率，有利于有限存儲空間下頻譜分辨率的提高，是一種有實(shí)際應(yīng)用價(jià)值的方法。

1 故障模型分析①

傳統(tǒng)定子電流分析方法中，將定子電流頻譜中的(1±2ks)f1頻率成分作為轉(zhuǎn)子斷條的依據(jù)。實(shí)際情況中，上述成分極易被基頻成分淹沒，故障提取具有一定的難度。事實(shí)上轉(zhuǎn)子斷條故障會通過磁場作用進(jìn)一步影響轉(zhuǎn)速信號，使其在時(shí)域、頻域上出現(xiàn)一系列易于識別的故障信息。

1.1 定子電流中(1-2ks)f1成分的產(chǎn)生

假設(shè)異步電機(jī)定子接入頻率為f1的三相對稱交流電時(shí)，轉(zhuǎn)速為ω。此時(shí)轉(zhuǎn)子中感應(yīng)出轉(zhuǎn)差頻率f2=sf1的電流，并產(chǎn)生相對于轉(zhuǎn)子轉(zhuǎn)速為ω2=2πf2/p的旋轉(zhuǎn)磁場。

斷條故障發(fā)生后，轉(zhuǎn)子出現(xiàn)電量不對稱。根據(jù)旋轉(zhuǎn)磁場理論，轉(zhuǎn)子的不對稱性會使氣隙中產(chǎn)生相對于轉(zhuǎn)子轉(zhuǎn)速±2πf2/p的兩個(gè)大小相等、方向相反的附加旋轉(zhuǎn)磁場。其中正向旋轉(zhuǎn)的磁場在轉(zhuǎn)速為ω的轉(zhuǎn)子的承載下，保持與定子電流旋轉(zhuǎn)磁場同速；而反向附加磁場相對于定子繞組的轉(zhuǎn)速ωb可由下式計(jì)算得到：

(1)

式中f——轉(zhuǎn)子轉(zhuǎn)頻；

f1——電網(wǎng)頻率；

p——電機(jī)的極對數(shù)；

s——電機(jī)轉(zhuǎn)差率。

可見，當(dāng)異步電機(jī)出現(xiàn)轉(zhuǎn)子斷條等不對稱故障時(shí)，氣隙中會產(chǎn)生相對于定子轉(zhuǎn)速為ωb的附加旋轉(zhuǎn)磁場，此磁場會在定子繞組中感應(yīng)產(chǎn)生頻率fb=(1-2s)f1的感應(yīng)電流。許多研究將此分量作為定子電流中轉(zhuǎn)子斷條的故障特征，并通過在定子電流中提取該頻率成分的方式進(jìn)行故障識別。然而，對于異步電動機(jī)，由于轉(zhuǎn)差率s的值很小(1%～5%),將導(dǎo)致轉(zhuǎn)子電量不對稱特征頻率(1±2s)f1與50Hz基頻成分相距很近，這樣一來，故障特征成分就極易被幅值強(qiáng)大的基頻成分覆蓋，影響故障識別的準(zhǔn)確性。因此傳統(tǒng)的定子電流頻譜分析方法難以滿足故障診斷的要求。

1.2 轉(zhuǎn)子故障對轉(zhuǎn)速的影響

故障電機(jī)運(yùn)行過程中，故障信息可以通過氣隙磁場的作用影響轉(zhuǎn)子轉(zhuǎn)矩，進(jìn)而影響瞬時(shí)轉(zhuǎn)速。轉(zhuǎn)速信號是電機(jī)運(yùn)行過程中的重要參數(shù)，其中包含大量電機(jī)運(yùn)行狀態(tài)信息，此外，轉(zhuǎn)速信號可通過光電編碼器等傳感器方便地獲取，信號采集技術(shù)成熟，誤差小、響應(yīng)快，造價(jià)低廉。因此，轉(zhuǎn)速信號適用于對電機(jī)進(jìn)行狀態(tài)監(jiān)測。

假設(shè)上述定子電流中故障附加電流分量瞬時(shí)表達(dá)式為：

ib1(t)=Ib1cos[(1-2s)ω1t-αb1]

(2)

其中，Ib1為此電流分量的幅值，αb1為其初相角。

此電流分量與基波磁通Φ=cos(ω1t-αΦ)作用，使電機(jī)電磁轉(zhuǎn)矩發(fā)生波動，頻率為2sf1，則有：

ΔT(t)=3pΦIb1sin[2sω1t-(αΦ-αb1)]

(3)

式中αΦ——定子基波磁通初相角；

Φ——定子基波磁通幅值。

由于轉(zhuǎn)矩的波動，導(dǎo)致轉(zhuǎn)速隨之發(fā)生相應(yīng)變化，轉(zhuǎn)速變化量可由轉(zhuǎn)矩變化量的積分計(jì)算：

(4)

其中，J為發(fā)電機(jī)轉(zhuǎn)子轉(zhuǎn)動慣量。

可見，轉(zhuǎn)子斷條故障下，轉(zhuǎn)速受到轉(zhuǎn)矩脈動的影響，也出現(xiàn)頻率為2sf1的周期性波動，這表明故障電機(jī)瞬時(shí)轉(zhuǎn)速信號中將會包含相同頻率的特征成分。因此對瞬時(shí)轉(zhuǎn)速信號采樣后進(jìn)行時(shí)域及頻域等角度分析，可以提取異步電機(jī)轉(zhuǎn)子斷條的故障信息。此外，穩(wěn)態(tài)下，轉(zhuǎn)速信號低頻段中僅包含直流成分，頻譜相對簡單，易于故障識別和提取。

實(shí)際上，由于轉(zhuǎn)速出現(xiàn)了頻率為2sf1的周期性波動，靜止的定子繞組中的電流會受到頻率為2sf1的調(diào)制，導(dǎo)致其中生成(1±2ks)f1一系列特征頻率成分。

2 有限元故障建模與仿真

眾所周知，在故障狀態(tài)下電機(jī)三相不對稱，傳統(tǒng)交流電機(jī)數(shù)學(xué)模型不再適用，從數(shù)學(xué)模型的角度進(jìn)行故障建模具有極大的難度，而對實(shí)體電機(jī)進(jìn)行故障實(shí)驗(yàn)成本高昂，因此對于電機(jī)故障研究來說，數(shù)據(jù)獲取方式始終是一個(gè)關(guān)鍵難題。Maxwell是一款針對于電磁場仿真的有限元分析軟件，廣泛應(yīng)用于電機(jī)設(shè)計(jì)等領(lǐng)域，參數(shù)設(shè)置靈活，適用于故障分析，結(jié)果具有一定的可靠性。

圖1所示為鼠籠型異步電機(jī)有限元模型，電機(jī)極對數(shù)為2，額定功率為11kW，額定轉(zhuǎn)速為1 462r/min，定子36槽，轉(zhuǎn)子24槽，轉(zhuǎn)子機(jī)械轉(zhuǎn)動慣量為0.096 82N·m?？梢酝ㄟ^減小指定導(dǎo)條的電導(dǎo)率來模擬轉(zhuǎn)子斷條。

圖1 鼠籠型異步電機(jī)有限元模型

令電機(jī)在額定負(fù)載(80N·m)下運(yùn)行，采樣率為200Hz。正常狀態(tài)下的電機(jī)轉(zhuǎn)速、電磁轉(zhuǎn)矩、電流信號及其頻譜如圖2所示?？梢姡顟B(tài)下，電機(jī)轉(zhuǎn)速和轉(zhuǎn)矩平穩(wěn)無周期性波動，電流信號中僅包含50Hz基頻成分。

圖2 正常狀態(tài)下參數(shù)的仿真結(jié)果

圖3所示為瞬時(shí)轉(zhuǎn)速信號經(jīng)去直流運(yùn)算后FFT頻譜的低頻(0～30Hz)部分，可以看出，轉(zhuǎn)速信號中僅含直流成分，頻譜成分簡單清晰。

圖3 正常電機(jī)瞬時(shí)轉(zhuǎn)速去直流后的FFT頻譜

在相同的運(yùn)行和采樣條件下，設(shè)置電機(jī)一根斷條，所得穩(wěn)態(tài)下的仿真結(jié)果如圖4所示。由于斷條，轉(zhuǎn)速出現(xiàn)了周期性變化(s≈0.025，波動周期Tb=1/2sf1≈1/2.5≈0.4s=400ms)。從電流時(shí)域波形可見，在轉(zhuǎn)子故障的影響下，電流波形被調(diào)制而出現(xiàn)失真，然而由于頻譜泄漏等原因，其FFT譜上的故障特征頻率(1±2ks)f1成分并不明顯，只有在(1+4s)f1≈55Hz處出現(xiàn)了尖峰，與故障相對應(yīng)，但與基頻相比，故障特征頻率幅值極小，難以分辨。

圖4 電機(jī)一根斷條時(shí)額定負(fù)載下仿真結(jié)果

將瞬時(shí)轉(zhuǎn)速信號經(jīng)過去直流后，做FFT運(yùn)算，得到如圖5所示的頻譜低頻部分，可清晰地看見包含的2sf1=2.5Hz成分及其倍頻4sf1=5Hz，證明轉(zhuǎn)速在以此頻率進(jìn)行周期性波動。

圖5 電機(jī)一根斷條時(shí)額定負(fù)載下轉(zhuǎn)速的FFT頻譜

值得一提的是，轉(zhuǎn)速信號中頻率成分簡單，特征頻率附近無干擾成分，非常易于故障的識別。此外，由于異步電動機(jī)轉(zhuǎn)差率很小(1%～5%)，故轉(zhuǎn)速中故障特征頻率一般小于5Hz，所以，在工程應(yīng)用中，存儲空間有限時(shí)，可以大幅減小轉(zhuǎn)速信號的采樣率，進(jìn)而延長采樣時(shí)間，將頻譜分辨率提高幾十倍，提高故障分辨的準(zhǔn)確性。

3 故障程度對故障特征的影響

事實(shí)上，故障后由于相鄰導(dǎo)條電流及所受電磁力增大，通常故障會進(jìn)一步發(fā)展為連續(xù)斷條的情況。現(xiàn)通過實(shí)驗(yàn)，分析斷條根數(shù)對瞬時(shí)轉(zhuǎn)速等參數(shù)的影響程度。

在相同的運(yùn)行和采樣條件下，為電機(jī)配置了1～3根連續(xù)斷條。不同故障程度下定子電流和瞬時(shí)轉(zhuǎn)速的頻譜對比如圖6所示?？梢钥闯?，隨著故障程度的增加，定子電流中的特征頻率成分幅值增大(如(1-2s)f1成分幅值由-6上升至2)，且諧波成分也趨于明顯。轉(zhuǎn)速的頻譜中，2sf1成分幅值同樣隨著斷條根數(shù)的增加而變大，且變化更為明顯(由0.4增大至8.0)。在轉(zhuǎn)子為3根斷條的情況下，由于故障嚴(yán)重導(dǎo)致轉(zhuǎn)速波動較大，導(dǎo)致0～1Hz低頻成分幅值增大，但2sf1成分仍可輕易分辨。

圖6 不同故障程度下定子電流和瞬時(shí)轉(zhuǎn)速頻譜對比

轉(zhuǎn)速原始信號高頻噪聲較大，而本方法關(guān)注的頻率范圍是5Hz以下的低頻帶，將3種情況下的轉(zhuǎn)速信號進(jìn)行低通濾波后，所得曲線如圖7所示。可見，轉(zhuǎn)速的波動幅度隨著斷條根數(shù)的增加而增大，且穩(wěn)態(tài)轉(zhuǎn)速的大小有著隨斷條根數(shù)的增加而下降的趨勢。

圖7 不同故障程度下的瞬時(shí)轉(zhuǎn)速

表1統(tǒng)計(jì)了3種故障狀態(tài)下由電機(jī)轉(zhuǎn)速信號計(jì)算而得的統(tǒng)計(jì)量，可見，隨著斷條根數(shù)的增加，電機(jī)穩(wěn)態(tài)時(shí)的平均轉(zhuǎn)速有所下降，而電磁轉(zhuǎn)矩卻未發(fā)生變化，這種現(xiàn)象的主要原因?yàn)椋簲鄺l故障導(dǎo)致轉(zhuǎn)矩和轉(zhuǎn)速信號波動的增加，消耗了部分定子輸入的能量，導(dǎo)致轉(zhuǎn)速的下降。

表1 不同故障程度下轉(zhuǎn)速和轉(zhuǎn)矩?cái)?shù)據(jù)統(tǒng)計(jì)

4 結(jié)束語

筆者提出了一種基于轉(zhuǎn)速分析的異步電機(jī)轉(zhuǎn)子斷條診斷方法，通過利用Maxwell有限元分析軟件進(jìn)行故障建模并仿真，從時(shí)域、頻域多個(gè)角度分析了轉(zhuǎn)速信號中的故障信息。實(shí)驗(yàn)結(jié)果表明：瞬時(shí)轉(zhuǎn)速中2sf1成分能夠有效反映故障情況，頻率成分明顯，清晰易辨。此外，瞬時(shí)轉(zhuǎn)速波動程度、平均轉(zhuǎn)速值均受連續(xù)斷條根數(shù)的影響，可作為判斷故障嚴(yán)重程度的依據(jù)。與傳統(tǒng)的定子電流分析方法進(jìn)行對比發(fā)現(xiàn)，轉(zhuǎn)速分析方法簡單易行、準(zhǔn)確性高。特別地，這種方法降低了對信號采樣率的要求，有利于有限存儲空間下頻譜分辨率的提高，是一種有實(shí)際應(yīng)用價(jià)值的方法。

[1] 宋軍艦，郝點(diǎn)，曹建強(qiáng).高速轉(zhuǎn)子故障模擬實(shí)驗(yàn)研究[J].化工機(jī)械，2011,38(4):417～419.

[2] 常穎，李聰，張風(fēng)波.離心機(jī)轉(zhuǎn)子力學(xué)性能的有限元分析及測試[J].化工機(jī)械，2012,39(4):478～481,502.

[3] 方瑞明，鄭力新，馬宏忠，等.基于MCSA和SVM的異步電機(jī)轉(zhuǎn)子故障診斷[J].儀器儀表學(xué)報(bào)，2004，24(5)：115～119.

[4] 王維琨，江志農(nóng)，張進(jìn)杰.基于神經(jīng)網(wǎng)絡(luò)和瞬時(shí)轉(zhuǎn)速的發(fā)動機(jī)失火故障研究[J].機(jī)電工程，2013，30(7)：824～827.