何玉靈 彭勃 萬書亭

摘要： 對氣隙靜偏心和定子短路故障的相對位置對定子振動特性的影響進行解析分析、數(shù)值計算和實驗研究。首先推導得到了包含定子繞組短路位置信息的作用于定子單位面積的磁拉力解析表達式，分析得到了定子短路位置對磁拉力的影響；然后建立了算例發(fā)電機正常運行、氣隙靜偏心故障、氣隙靜偏心與定子匝間短路復合故障下的有限元模型，對不同短路位置下的定子磁拉力進行了仿真計算。并實測了MJF306型故障模擬發(fā)電機定子不同短路位置的定子振動信號，三者結(jié)果相互吻合。結(jié)果表明：氣隙靜偏心與定子匝間短路復合故障下，同等故障程度時，發(fā)電機定子振動幅度與短路位置有關(guān)，短路位置越接近最小氣隙處，各倍頻振動幅度越大，越遠離最小氣隙處，各倍頻振動幅度越小。關(guān)鍵詞： 故障診斷； 發(fā)電機； 定子振動； 氣隙靜偏心； 定子匝間短路

中圖分類號： TH165.3； TM311文獻標志碼：A文章編號：10044523（2017）04067909

DOI：10.16385/j.cnki.issn.10044523.2017.04.020

引言

由于制造和運行等原因，發(fā)電機定轉(zhuǎn)子之間的氣隙或多或少存在不均勻狀況（即氣隙偏心），當偏心程度大于10%時，將會引起發(fā)電機的劇烈振動，嚴重時，將會造成定子鐵芯變形、繞組損壞等事故[1]。文獻[23]研究了氣隙偏心故障下，發(fā)電機三相繞組的電流或電壓值的變化；文獻[45]分析氣隙偏心下發(fā)電機的定子并聯(lián)支路環(huán)流特性，并提出相應識別方法。

定子繞組匝間短路作為運行人員主要關(guān)注的電氣故障之一也時有發(fā)生，是一種危害較為嚴重，修復費用較高的故障[6]。文獻[7]分析了定子短路匝間對發(fā)電機感應電動勢的影響，通過故障前后電勢的差異來識別故障；文獻[8]分析了定子匝間短路故障下定子并聯(lián)支路環(huán)流特性的變化，并討論了短路參數(shù)變化對環(huán)流特性的影響；文獻[9]在發(fā)電機氣隙中安裝探測線圈，通過分析端口電壓頻率特征差異判斷對定子匝間短路故障。

事實上，氣隙靜偏心是普遍存在的，定子匝間短路故障發(fā)生時，發(fā)電機往往處于氣隙靜偏心與定子匝間短路的復合故障狀態(tài)。同時由于發(fā)電機是機電耦合的整體，故障將引起氣隙磁場和作用于定轉(zhuǎn)子的磁拉力變化，從而激勵發(fā)電機定轉(zhuǎn)子振動，并且振動信號測試更加便捷。因此已有文獻研究故障下的振動特性，文獻[10]研究了氣隙靜偏心和定子短路復合故障的定子振動特性及其變化規(guī)律，但沒有考慮故障位置對其特性的影響。文獻[11]研究了考慮氣隙偏心時發(fā)電機定子匝間短路位置對電磁轉(zhuǎn)矩波動特性的影響。

本論文將在文獻[1011]基礎(chǔ)上，首次對氣隙靜偏心和定子短路故障的相對位置對定子振動特性的影響進行解析分析、數(shù)值計算和實驗研究，分析結(jié)果是對發(fā)電機氣隙靜偏心與定子匝間短路復合故障診斷的一個重要補充，具備積極意義。

1定子振動特性分析〖*2〗1.1定子受力分析發(fā)電機的定子由定子鐵芯、繞組、機座、端蓋等部件組成，其中定子鐵芯是磁拉力激勵的主要部件。定子鐵芯單位面積磁拉力表達式如下qαm，t=Bαm，t22μ0=fαm，tΛαm，t22μ0（1）式中B為氣隙磁密，αm為氣隙周向角，Λ為單位面積氣隙磁導，f為氣隙磁勢，μ0為真空磁導率。

文獻[10]在推導氣隙靜偏心與定子匝間短路復合故障下的發(fā)電機氣隙磁勢時，認為短路環(huán)電流產(chǎn)生的脈振磁勢恒定，忽略了定子短路點與靜偏心的相對位置對定子振動的影響。事實上，如圖1 所示，以氣隙最小處作為原點，氣隙靜偏心時，氣隙磁密會隨著氣隙周向角的變化而變化。因此，定子匝間短路位置不同（定義定子繞組短路位置角α′m為定子繞組短路匝中心位置與氣隙最小處的夾角），短路環(huán)電流值也將不同，進一步導致其產(chǎn)生的脈振磁勢幅值有所不同。

圖1發(fā)電機靜態(tài)偏心氣隙

Fig.1Static airgap eccentricity of generator

作為改進，本文假設(shè)在氣隙最小處發(fā)生定子匝間短路（短路位置角α′m為0度），設(shè)短路環(huán)中附加環(huán)流產(chǎn)生的脈振磁勢為Fd，且幅值最大（磁導最大），故短路位置產(chǎn)生的脈振磁勢幅值與最大脈振磁勢幅值的比為ρ=fdFd=1-δs1-δscosα′m（2）式中g(shù)0為平均氣隙，δs為相對靜偏心值。以42槽汽輪發(fā)電機為例，假設(shè)最小氣隙位置位于1槽，則當定子匝間短路匝中心位置在1槽至42槽時，ρ值變化趨勢如圖2所示。

圖2短路匝中心位于不同槽下ρ值變化趨勢

Fig.2Developing trend of ρ when interturn shorted center is in different slots

從式（2）與圖2中可看出，定子匝間短路位置越接近最小氣隙處時產(chǎn)生的脈振磁勢幅值越大，反之，定子匝間短路位置越遠離最小氣隙處時產(chǎn)生的脈振磁勢幅值越小。此外，ρ值大小還與偏心程度有關(guān)系，偏心程度越大則上述規(guī)律越明顯。忽略高次諧波，不同位置的定子匝間短路產(chǎn)生的脈振磁勢可表示為fd（αm，t）=ρFdcosωtcospαm=

ρFd+cos（ωt-pαm）+ρFd-cos（ωt+pαm）（3）式中p為發(fā)電機極對數(shù)，ω為電角頻率。

第4期何玉靈，等： 定子匝間短路位置對發(fā)電機定子振動特性的影響振 動 工 程 學 報第30卷對比文獻[10]中公式可發(fā)現(xiàn)，考慮氣隙靜偏心后，式（3）中的第1項以同步速與轉(zhuǎn)子同向旋轉(zhuǎn)，在轉(zhuǎn)子繞組中不感應附加諧波電勢，第2項以同步速與轉(zhuǎn)子逆向旋轉(zhuǎn)，將在轉(zhuǎn)子繞組中感應2倍頻電流值ρIf2。以基頻轉(zhuǎn)速反向旋轉(zhuǎn)的電樞反應磁勢最大幅值變?yōu)棣袴s1-，以3倍頻轉(zhuǎn)速正向旋轉(zhuǎn)的電樞反應磁勢最大幅值變?yōu)棣袴s3+，最終影響定子磁拉力各次諧波幅值，詳細推導過程可見文獻[10]。根據(jù)上述內(nèi)容對定子磁拉力主要諧波成分幅值做出修改，如表1所示。直流分量（0 Hz）是常值，不會激勵定子產(chǎn)生振動，不列入表中。

分析表1可知，正常運行和氣隙靜偏心單故障下，作用于定子單位面積的磁拉力主要諧波成分為2次諧波（100 Hz），但靜偏心下諧波幅值較正常情況要大。復合故障下，定子磁拉力主要諧波成分為2次諧波（100 Hz）、4次諧波（200 Hz）、6次諧波（300 Hz）。復合故障下定子磁拉力各次諧波幅值比單故障大。此外，復合故障下，當保持氣隙靜偏心程度與定子匝間短路程度不變時，短路位置越接近最小氣隙處，定子磁拉力2次諧波、4次諧波以及6次諧波幅值增大越多，反之越少。各次諧波幅值隨周期變化，會激勵定子產(chǎn)生的相應倍頻振動。

1.2基于有限元法求解單位面積磁拉力

由于電機磁場計算時常忽略端部效應，認為磁場軸向均勻分布，故采用二維磁場進行計算。由上文所述，定子振動的本質(zhì)激勵源為單位面積磁拉力，計算定子所受磁拉力時應選取定子內(nèi)圓上的一點作為求解點，最終求解得到的磁拉力即為單位面積磁表1不同情況下定子磁拉力各次諧波成分幅值

Tab.1Stator electromagnetic force amplitude under different conditions

成分運行狀態(tài)幅值表達式2次

諧波正常I2f0N2Λ20/4μ0靜偏心0.5I2f0N2+0.5F2s-If0NFssinΛ20+Λ2s+2Λ0Λs/2μ0復合故障[（ρIf0If2N2+ρFs1+Fs1-+ρFs1+Fs3++ρIf0NFs1-+If0NFs1++ρIf2NFs1++ρ2If2NFs1-+ρIf0NFs3++

0.5ρ2If2NFs3++ρ2Fs1-Fs3++0.5F2s1++0.5ρ2F2s1-+0.5I2f0N2+0.5ρ2I2f2N2）（Λ20+Λ2s/2）+

0.25（ρ2F2s1-+F2s1++ρ2I2f2N2+I2f0N2+2ρIf0If2N2+2If0NFs1++2ρIf0NFs1-+2ρIf2NFs1++

2ρ2If2NFs1-+2ρIf0NFs3++ρ2If2NFs3++2ρFs1+Fs1-+2ρFs1+Fs3++2ρ2Fs1-Fs3+）Λ2s+

（F2s1++ρ2F2s1-+2ρFs1+Fs1-+2ρ2Fs1-Fs3++2ρFs1+Fs3++I2f0N2+ρ2I2f2N2+2ρIf0If2N2+

2If0NFs1++2ρIf2NFs1++2ρIf0NFs1-+2ρ2If2NFs1-+2ρIf0NFs3++ρ2If2NFs3+）Λ0Λs]/2μ04次

諧波正?！o偏心—復合故障[（ρ2I2f2N2/2+ρIf0If2N2+ρIf2NFs1++ρ2If2NFs1-+ρ2If0NFs3++ρ2If2NFs3++2ρ2Fs1-Fs3++

2ρFs1+Fs3+）（Λ20+Λ2s/2）/2+（2ρIf0If2N2+ρ2I2f2N2+ρ2If2NFs1++2ρ2If2NFs1-+

4ρIf0NFs3++2ρ2If2NFs3++4ρ2Fs1-Fs3++4ρFs1+Fs3+）Λ2s/8+（ρ2I2f2N2/2+ρIf0If2N2+

2ρIf0NFs3++ρIf2NFs1++ρ2If2NFs1-+ρ2If2NFs3++2ρ2Fs1-Fs3++2ρFs1+Fs3+）Λ0Λs]/2μ06次

諧波正?！o偏心—復合故障[ρ2（F2s3++I2f2N2/2+If2NFs3+）Λ2s+ρ2（4F2s3++4If2NFs3++2I2f2N2）Λ0Λs+

ρ2（2F2s3++2If2NFs3++I2f2N2）（Λ20+Λ2s/2）]/8μ0

拉力。本文使用電機有限元軟件Ansoft作為計算工具，仿真運行后，可以直接獲取發(fā)電機氣隙磁密在所選求解點處X軸分量Bx（αm，t）與Y軸分量By（αm，t），再通過Field Calculator對公式（4）編輯并進行計算，即可求得求解點處徑向氣隙磁通密度Bαm，t=Bxαm，tcosθ+Byαm，tsinθ（4）式中θ為求解時使用的柱坐標。

聯(lián)立公式（1），（4），最終得到作用于發(fā)電機定子的單位面積磁拉力的表達式為qαm，t=Bxαm，tcosθ+Byαm，tsinθ22μ0（5）再次通過Field Calculator對公式（5）編輯并進行計算，求得單位面積磁拉力，并分析頻率成分。

2有限元仿真與動模實驗〖*2〗2.1實驗設(shè)備與測試方法本文利用MJF306故障模擬發(fā)電機進行實驗驗證，總體外觀如圖3（a）所示。該型模擬發(fā)電機的定轉(zhuǎn)子均無法移動，不能對氣隙靜偏心故障進行模擬。但考慮到制造精度、裝配誤差以及長時間運行等因素的影響，判斷實驗室發(fā)電機可能會存在一定程度的氣隙靜偏心。利用高精塞尺對定轉(zhuǎn)子間的氣隙進行實際測量，如圖3（b）所示，發(fā)現(xiàn)該模擬發(fā)電機勵端視圖上偏左22°區(qū)域為氣隙最小處，gmin=0.80 mm，下偏右22°區(qū)域為氣隙最大處，gmax=0.90 mm，偏心示意圖如圖3（c）所示。此外，該型模擬發(fā)電機引出多個定子繞組短路抽頭，短路抽頭設(shè)置及其接線盤如圖3（d），（e）所示。其中，定子C1支路短路5%抽頭C12所在繞組位置與氣隙最小處位置機械夾角α′m=111°（短路位置遠離最小氣隙處），B1支路短路5%抽頭B12所在繞組位置與氣隙最小處位置機械夾角α′m=60°，C2支路短路5%抽頭C22所在繞組中心與氣隙最小處機械夾角α′m=2°（短路位置接近最小氣隙處）。分別短接抽頭C12與C22，即可模擬靜偏心程度與短路程度一定而短路位置不同的復合故障。

實驗過程，并網(wǎng)運行，勵磁電流If0=1.3 A，P=10 kW，Q=4 kVar。采用CD21C型振動速度傳感器對發(fā)電機定子鐵芯的豎直方向振動進行測取，其靈敏度為30 mV/mm/s，如圖3（f）所示。將振動速度傳感器的輸出端接入至U60116型采集儀，采樣頻率為5 kHz，實驗測試系統(tǒng)如圖3（g）所示。模擬實驗分4組，具體方案如下：

1） 不短接抽頭，模擬氣隙靜偏心單故障，采集定子振動數(shù)據(jù)。

2） 短接抽頭C12，模擬氣隙靜偏心與定子匝間短路復合故障1，采集定子振動數(shù)據(jù)。

3） 短接抽頭B12，模擬氣隙靜偏心與定子匝間短路復合故障2，采集定子振動數(shù)據(jù)。

圖3MJF306型故障模擬機組

Fig.3MJF306 type fault simulating generator set

4）短接抽頭C22，模擬氣隙靜偏心與定子匝間短路復合故障3，采集定子振動數(shù)據(jù)。

2.2有限元建模與仿真

本文對MJF306型隱極故障模擬發(fā)電機進行有限元建模仿真，其相關(guān)參數(shù)如表2所示。表2MJF306型故障模擬發(fā)電機主要參數(shù)

Tab.2Primary parameters of MJF306 type generator

參數(shù)名稱參數(shù)值額定容量30 kVA額定電壓400 V轉(zhuǎn)速1000 r/min定子槽數(shù)54槽轉(zhuǎn)子槽數(shù)30槽定子每槽匝數(shù)72匝轉(zhuǎn)子每槽匝數(shù)88匝氣隙徑向長度0.85 mm在有限元軟件Ansoft的RMxprt電機模塊環(huán)境下，根據(jù)表2輸入相關(guān)參數(shù)，即可自動生成發(fā)電機有限元模型，且該模型能對發(fā)電機正常運行狀態(tài)進行仿真。氣隙靜偏心故障的仿真方法是使轉(zhuǎn)子軸心與定子軸心產(chǎn)生偏移。將RMxprt模塊自動生成的發(fā)電機有限元模型導入Maxwell 2D模塊，按照圖3（c）所示偏心情況對發(fā)電機定子鐵芯、定子繞組及OuterRegion進行移動，即可得到氣隙靜偏心故障下的發(fā)電機有限元模型，且與實驗機組的偏心情況保持一致。定子短路為電氣故障，不會對定轉(zhuǎn)子的氣隙間距產(chǎn)生影響，故直接對氣隙靜偏心故障下的發(fā)電機有限元模型進行二次修改即可。利用Ansoft Maxwell Circuit Editor軟件繪制定轉(zhuǎn)子外電路，如圖4（a）所示，并將其導入Maxwell 2D模塊，使發(fā)電機有限元模型的繞組與外電路相耦合。以B相繞組為例，圖中LB1，LB2，LB3分別為發(fā)電機B相第一條支路繞組，RB1為B相第一條支路總圖4發(fā)電機數(shù)值仿真分析模型

Fig.4Analysis model of turbogenerator for simulation

電阻，RB為B相負載電阻。LshortB3為短路故障下LB3中被短路的線圈，RB3S為LshortB3電阻，RB3為短路故障下LB3中未被短路的線圈電阻。SB3為壓控開關(guān)，R_SB3為SB3的電阻。V_SB3為脈沖電壓源， -1 V為脈沖電壓源的初始電壓，1 V為脈沖電壓源的峰值電壓，用來提供開關(guān)動作電壓。ShortModel_SB3為SB3的模型元件，壓控開關(guān)模型元件需要與壓控開關(guān)一起使用。DC是勵磁電流源。RA是為保證電流源所提供電流全部通過轉(zhuǎn)子繞組而添加的電阻，與繞組本身結(jié)構(gòu)無關(guān)。LF1，LF2，LF3，LF4，LF5，LF6分別是轉(zhuǎn)子繞組。Rf是轉(zhuǎn)子繞組的總電阻。

至此，建立了氣隙靜偏心與定子短路復合故障下的發(fā)電機有限元模型，二維分析模型如圖4（b）所示，有限元網(wǎng)格剖分情況如圖4（c）所示，氣隙及繞組槽內(nèi)網(wǎng)格最大邊長為0.008 mm，得到10710個網(wǎng)格單元。

仿真過程中，對轉(zhuǎn)子繞組施加1.3 A電流，與實驗過程保持一致。3個開關(guān)SB3，SC3和SC6設(shè)置有相應的控制電路，如圖4（a）所示。SB3，SC3和SC6均不閉合時，仿真氣隙靜偏心單故障；SC3閉合時，短路繞組LshortC3被短接，仿真復合故障1；SB3閉合時，短路繞組LshortB3被短接，仿真復合故障2；SC6閉合時，短路繞組LshortC6被短接，仿真復合故障3。

3結(jié)果分析與討論〖*2〗3.1有限元仿真結(jié)果Ansoft仿真計算得到3種運行狀況下的單位面積磁拉力頻譜，如圖5和表3所示。

圖5仿真得到的不同狀態(tài)下單位面積磁拉力頻譜

Fig.5Magnetic pull per unit area （MPPUA） spectra under different conditions from simulation

表3定子單位面積磁拉力各次諧波幅值

Tab.3Harmonics amplitudes of stator MPPUA

運行狀態(tài)定子磁拉力/（N·m-2）2次/100 Hz4次/200 Hz6次/300 Hz正常運行64.5600氣隙靜偏心單故障68.7500復合故障175.0811.864.05復合故障288.9314.664.42復合故障3108.2616.035.16

從圖5和表3可以看出，復合故障2的磁拉力各次諧波幅值都大于復合故障1，與解析分析結(jié)果一致。

3.2實驗測試結(jié)果

實驗室故障模擬發(fā)電機存在氣隙靜偏心，無法采集到嚴格無偏心情況下的定子振動數(shù)據(jù)。作為對比，本文調(diào)取了2002年不設(shè)置短路故障采集的定子振動數(shù)據(jù)，并對其與2015年采集的數(shù)據(jù)做傅里葉變換，分析發(fā)電機振動特征，原始頻譜如圖6所示。

由文獻[10]可知，發(fā)電機氣隙靜偏心和定子短路復合故障下定子振動的特征頻率成分為100，200與300 Hz，而從圖6可看出，原始定子振動信號中故障特征振動頻率成分完全被環(huán)境圖6原始定子振動頻譜

Fig.6Original stator vibration spectra

噪聲淹沒，無法對發(fā)電機故障進行診斷與識別。因此，本文利用文獻[12]提出的一種基于正弦結(jié)構(gòu)元素的自適應改進TopHat變換算法，旨在有效消除原始信號中環(huán)境噪聲，同時增大故障特征振動頻率的幅值，去噪過程中4組信號采用相同的處理參數(shù)，對應的算法流程如圖7所示。

圖7算法流程圖

Fig.7Algorithm flow chart

該算法主要步驟如下：

（1） 構(gòu)造結(jié)構(gòu)元素為正弦，對待分析原始信號進行自互補TopHat變換。由于形態(tài)學濾波所選結(jié)構(gòu)元素的形狀與待分析信號的形狀越相近，濾波效果越好。發(fā)電機在定轉(zhuǎn)子短路、氣隙偏心故障下定子振動信號頻率為50 Hz（電頻率）整數(shù)倍的正弦波，所以構(gòu)造正弦結(jié)構(gòu)元素進行濾波處理，可獲得最佳的濾波效果。

（2） 正弦結(jié)構(gòu)元素的幅值及長度經(jīng)嵌套遍歷循環(huán)選取不同數(shù)值并與原始信號進行自互補TopHat變換，計算變換后信號中故障振動特征頻率的能量與信號總能量的比值，選擇遍歷過程中能量比最大時所對應的幅值和長度作為正弦結(jié)構(gòu)元素的最佳幅值及長度。

（3） 將待分析原始信號與最佳幅值及長度的正弦結(jié)構(gòu)元素進行自互補TopHat變換，計算增強后信號的頻譜，并與理論分析結(jié)果作對比，對發(fā)電機故障進行診斷與識別。

增強后頻譜如圖8所示。

圖8增強后振動頻譜

Fig.8Enhanced stator vibration spectra

觀察圖8（a）可知，2002年定子振動頻譜中2倍頻譜線突出，幅值為0.5914 mm/s。對比圖8（b）可知，歷經(jīng)13年運行，發(fā)電機出現(xiàn)氣隙靜偏心故障，2倍頻譜線幅值從0.5914 mm/s增加到了0.6985 mm/s，增幅為18.1%。復合故障1下2倍頻振幅為0.8151 mm/s，相對靜偏心增長了16.7%；4倍頻振幅為0.1834 mm/s，6倍頻振幅為0.0831 mm/s。復合故障2下2倍頻振幅為1.0329 mm/s，相對復合故障1增長了26.7%；4倍頻振幅為0.2271 mm/s，相對復合故障1增長了23.8%；6倍頻振幅為0.0982 mm/s，相對復合而故障1增長了18.2%。復合故障3下2倍頻振幅為1.1825 mm/s，相對靜偏心相對復合故障1增長了45.1%；4倍頻振幅為0.2504 mm/s，相對復合故障1增長了36.5%；6倍頻振幅為0.1067 mm/s，相對復合故障1增長了27.2%。

綜上可知， 理論解析推導、有限元仿真結(jié)果以及實驗測試數(shù)據(jù)相互吻合。

4結(jié)論

本文對氣隙靜偏心和定子短路故障的相對位置對定子振動特性的影響進行解析分析、數(shù)值計算和實驗研究，結(jié)果表明：

（1） 正常運行和氣隙靜偏心單故障下，發(fā)電機定子產(chǎn)生2倍頻振動；氣隙靜偏心與定子匝間短路復合故障下，發(fā)電機定子產(chǎn)生2倍頻、4倍頻與6倍頻振動。發(fā)電機分別處于正常、單故障以及復合故障下，定子振動會依次增強。

（2） 氣隙靜偏心與定子匝間短路復合故障下，同等故障程度時，發(fā)電機定子振動幅度與短路位置有關(guān)，短路位置越接近最小氣隙處，各倍頻振動幅度越大，越遠離最小氣隙處，各倍頻振動幅度越小。

論文成果在現(xiàn)有研究基礎(chǔ)上進一步完善了氣隙靜偏心與定子匝間短路故障的診斷判據(jù)。

參考文獻：

[1]譚善軍. 汽輪發(fā)電機偏心故障的交叉特征與診斷方法研究[J].中國電業(yè)（技術(shù)版）， 2012， （12）： 39—43.

TAN Shanjun. Research on steam turbine generator eccentric fault cross characteristics and diagnostic methods[J]. China Electric Power （Technology Edition）， 2012， （12）： 39—43.

[2]Blodt M， Regnier J， Faucher J. Distinguishing load torque oscillations and eccentricity faults in induction motors using stator current Wigner distributions [J]. IEEE Transactions on Industry Applications， 2009， 45（6）： 1991—2000.

[3]趙向陽， 葛文韜. 基于定子電流法監(jiān)測無刷直流電動機轉(zhuǎn)子動態(tài)偏心的故障模型仿真研究[J]. 中國電機工程學報， 2011， 31（36）： 124—130.

ZHAO Xiangyang，GE Wentao. Simulation research of fault model of detecting rotor dynamic eccentricity in brushless DC motor based on motor current signature analysis [J]. Proceedings of the CSEE， 2011， 31（36）： 124—130.

[4]萬書亭， 何玉靈， 唐貴基， 等. 發(fā)電機氣隙偏心時定子并聯(lián)支路的環(huán)流特性分析[J].高電壓技術(shù)， 2010， 36（6）： 1547—1553.

WAN Shuting， HE Yuling， TANG Guiji， et al. Analysis on stator circulating current characteristics under eccentricity faults of turbogenerator[J]. High Voltage Engineering， 2010， 36（6）： 1547—1553.

[5]Jian Li，Yunhyun Cho. Dynamic reduction of unbalanced magnetic force and vibration in switched reluctance motor by the parallel paths in windings [J]. Mathematics & Computers in Simulation， 2010， 81（2）： 407—419.

[6]Prasad B J C， Sanker Ram B V. Interturn fault analysis of synchronous generator using finite element method （FEM） [J]. International Journal of Recent Technology and Engineering， 2013， 2（5）： 150—156.

[7]Sarikhani A， Mohammed O A. Interturn fault detection in PM synchronous machines by physicsbased back electromotive force estimation [J]. IEEE Transactions on Industrial Electronics， 2013， 60（8）： 3472—3484.

[8]何玉靈，萬書亭，唐貴基，等. 定子匝間短路對發(fā)電機并聯(lián)支路環(huán)流特性的影響[J].電機與控制學報， 2013， 17（3）：1—7.

HE Yuling，WAN Shuting，TANG Guiji， et al. Effects of stator interturn short circuit on the circulating current inside the parallel branch loop[J]. Electric Machines and Control， 2013， 17（3）： 1—7

[9]孫宇光，余錫文，魏錕，等. 發(fā)電機繞組匝間故障檢測的新型探測線圈[J]. 中國電機工程學報， 2014， 34（6）： 917—924.

SUN Yuguang， YU Xiwen， WEI Kun， et al. A New type of search coil for detecting interturn faults in synchronous machines [J]. Proceedings of the CSEE， 2014，34（6）： 917—924.

[10]唐貴基， 何玉靈， 萬書亭， 等. 氣隙靜態(tài)偏心與定子短路復合故障對發(fā)電機定子振動特性的影響[J].振動工程學報， 2014， 27（1）： 118—127.

TANG Guiji， HE Yuling， WAN Shuting， et al. Effect of static eccentricity & stator short circuit composite faults on stator vibration characteristics of generator [J]. Journal of Vibration Engineering， 2014， 27（1）： 118—127.

[11]何玉靈，王發(fā)林，唐貴基，等. 考慮氣隙偏心時發(fā)電機定子匝間短路位置對電磁轉(zhuǎn)矩波動特性的影響[J].電工技術(shù)學報，2017，3（7）：11—19.

HE Yuling，WANG Falin，TANG Guiji，et al. Effect of stator interturn short circuit position on electromagnetic torque of generator with consideration of airgap eccentricity[J]. Transactions of China Electrotechnical Society， 2017，3（7）：11—19.

[12]何玉靈，蒙玉超，唐貴基，等.基于正弦結(jié)構(gòu)元素的自適應 TopHat 變換及發(fā)電機特征振動信號增強檢測[J].中國電機工程學報，2016，36（15）：4266—4273.

HE Yuling， MENG Yuchao， TANG Guiji， et al.SelfAdapted TopHat transformation based on sinestructure element and its application in enhanced characteristic vibration signal detection of generator[J].Proceedings of the CSEE，2016，36（15）：4266—4273.

Effect of stator interturn short circuit position on stator vibration of generatorHE Yuling， PENG Bo， WAN Shuting（Department of Mechanical Engineering， North China Electric Power University， Baoding 071003， China）Abstract： This paper analyzes the effect of the relative position of static airgap eccentricity and stator interturn short circuit on the stator vibration characteristic of generator via theoretical analysis， numerical simulation， and experimental study. Firstly， the detailed expressions of the stator magnetic pull per unit area with stator interturn short circuit location information are deduced， and the effect of short circuit position on magnetic pull is obtained. Then， the finite element model of example generator is set up to calculate the stator magnetic pull under normal condition， the static airgap eccentricity fault， and the composite fault （composed by static airgap eccentricity & stator interturn short circuit） with different shorted positions. Finally， the stator vibration signals of MJF306 type fault simulating generator at different stator short circuit positions are measured， which has the consistent result with the theoretical analysis and the numerical simulation. It is shown that， the stator vibration amplitude is associated with short circuit position when the generator under static airgap eccentricity & stator interturn short circuit composite fault. The shorted position is more near the min airgap position， the bigger the vibration amplitudes of each harmonic component. On the contrary， the shorted position is more far away the min airgap position， the smaller the vibration amplitudes of each harmonic component.Key words： fault diagnosis； turbo generator； stator vibration； static airgap eccentricity； stator interturn short circuit作者簡介： 何玉靈（1984—），男，博士，副教授。電話：（0312）7525043；Email：heyuling1@163.com

通訊作者： 萬書亭（1970—），男，教授，博士生導師。電話：（0312）7525026；Email：13582996591@139.com