鮑曉華，程曉巍，方勇，呂強，單麗

(1.合肥工業(yè)大學 電氣與自動化工程學院，安徽合肥 230009;2.合肥恒大江海泵業(yè)股份有限公司，安徽合肥 231131)

大型異步電機定子端部繞組電磁力的研究

鮑曉華1，程曉巍1，方勇1，呂強1，單麗2

(1.合肥工業(yè)大學 電氣與自動化工程學院，安徽合肥 230009;2.合肥恒大江海泵業(yè)股份有限公司，安徽合肥 231131)

針對大型異步電機定子端部繞組受電磁力作用容易發(fā)生故障的實際工程問題，提出了一種通過改變端部繞組錐度角進而改善定子端部繞組力學性能的方法。利用MAXWELL方程結(jié)合分離變量法推導出了端部繞組磁場的表達式，通過時步有限元方法建立端部繞組模型，進而分析繞組端部磁場的分布情況。最后結(jié)合洛倫茲法則，利用局部坐標變換得到端部繞組斜線段電磁力的計算公式，運用數(shù)值分析方法對其仿真分析，在此基礎(chǔ)上結(jié)合具體算例，對比分析了端部繞組在不同錐度角下電磁力的分布情況。分析結(jié)果表明，該方法能夠改善端部繞組的力學性能，延長電機使用壽命，為干式潛水電機定子端部繞組的優(yōu)化設(shè)計提供理論依據(jù)。

潛水電機;有限元方法;錐度角;端部繞組磁密;電磁力密度

0 引 言

大型異步電機廣泛應(yīng)用于關(guān)系到國民經(jīng)濟的重要場合，種類很多，高壓潛水電機是其中重要的一種，廣泛應(yīng)用于江河湖泊中取水、農(nóng)業(yè)排灌、城市或工廠的污水處理、礦山給排水、海上采油平臺抽取海水等重要場合，在國內(nèi)需求不斷增加的刺激下，潛水電機向大型化、高功率化的方向發(fā)展，因此定子繞組端部承受的電磁力也越來越大。電機在長期運行過程中，定子繞組端部承受強大的、周期性的電磁力作用，在其作用下所引發(fā)的電磁振動是導致繞組疲勞斷裂和引線絕緣磨損的主要原因之一[1-2]，特別是發(fā)生短路故障時，作用在定子繞組端部的電磁力超過額定運行狀態(tài)下的100倍，造成的危害及其嚴重。因此對定子端部繞組電磁力進行理論研究[3-4]和優(yōu)化計算，以便對干式潛水電機端部繞組提供更加合理的設(shè)計方案，提高電機運行的安全穩(wěn)定性，具有非常要的理論意義。

近年來，國內(nèi)外的專家學者對電機定子繞組端部電磁力進行了較為廣泛的研究。文獻[5-6]在忽略端區(qū)結(jié)構(gòu)影響的基礎(chǔ)上，運用有限差分法分別計算了電機在運行穩(wěn)定狀態(tài)和瞬時狀態(tài)下大型發(fā)電機定子端部繞組所受到的電磁力。文獻[7]在文獻[5-6]的基礎(chǔ)上針對大型異步電動機，采用有限元分析方法著重分析了電磁力的正弦分量對端部繞組的影響，忽略了其恒定分量的作用。文獻[8]僅考慮了由轉(zhuǎn)子電流激發(fā)的磁場對鼠籠型電機端部繞組所受電磁力的影響。R.Albanese等人在考慮到定轉(zhuǎn)子電流、鐵磁飽和等因素的影響下，采用非線性積分公式法闡述了繞組端部電磁力的計算思路，以及定子繞組在電磁力的作用下發(fā)生的位移變化[9]，但是該方法產(chǎn)生較大的計算誤差。文獻[10]對定子端部繞組整體結(jié)構(gòu)的動力問題進行了研究，推導了定子端部繞組整體振動控制方程，以及所受電磁力的表達式，并且通過數(shù)值計算得到了端部繞組的固有振動頻率、電磁力及動力響應(yīng)的計算結(jié)果。文獻[11]針對具體實例，考慮到定轉(zhuǎn)子鐵心端部的影響，運用有限元方法分析了磁場在定子繞組端部區(qū)域的分布情況，進而更為精確地計算出端部各區(qū)域中電磁力密度的動態(tài)響應(yīng)。以上各文獻在計算電機定子繞組端部電磁力方面取得了一定的成果，具有一定的實用價值。但是針對電磁力對端部繞組的破壞，并沒有提出相應(yīng)的改進辦法用以減弱電磁力對繞組端部的影響。

本文在前人研究的基礎(chǔ)上針對干式潛水電機（1000 kW/10 kV）建立繞組端部模型，通過有限元方法對定子繞組端部電磁場進行分析，在得到磁感應(yīng)場強度分布的基礎(chǔ)上，運用數(shù)值分析方法結(jié)合具體算例對繞組端部所受電磁力進行計算分析，提出了一種通過改變端部繞組錐度角進而改善繞組端部所受電磁力的方法，最后對不同的角度所對應(yīng)的電磁力進行仿真計算。

1 定子繞組端部磁場的分析

1.1 繞組端部磁場的產(chǎn)生

定子繞組端部區(qū)結(jié)構(gòu)復(fù)雜，如圖1所示，其磁場是由分布復(fù)雜的繞組端部電流所產(chǎn)生，又受到形狀復(fù)雜、材料不同的端部結(jié)構(gòu)件的影響，致使端部磁場的計算極為困難。通常僅考慮一些主要結(jié)構(gòu)件（定子鐵心壓板、端蓋、機座壁和轉(zhuǎn)軸等）的影響，并認為這些結(jié)構(gòu)件的形狀是規(guī)則的圓筒或平面；同時忽略其中的渦流阻尼效應(yīng)，把它們看作μ為無窮大，σ為零的理想磁鏡面。

圖1 端部繞組結(jié)構(gòu)Fig.1 Structure of end winding

繞組端部區(qū)域是由分屬于各相的許多線圈組成構(gòu)成，由于端部磁場主要由端部繞組斜線段電流產(chǎn)生，因此為求端部電磁場，將每層中屬于同一相繞組的斜線段電流密度沿x、y、z軸展成傅里葉級數(shù)。然后將三相疊加，得同一層繞組電流的表達式[12]。

斜線伸出段外層（1）段：

在直角坐標系下，通過給出的電流密度展開式，依據(jù)麥克斯韋方程組，即

式中：A(Ax,Ay,Az)為矢量磁位；J(Jx,Jy,Jz)為電流密度矢量；B(Bx,By,Bz)為磁感應(yīng)強度矢量；?2為拉普拉斯微分算子；?為哈密頓算子。

采用分離變量法對式（5）進行求解，將電流產(chǎn)生的磁場進行疊加，最終得到內(nèi)外層繞組磁感應(yīng)強度的表達式為

圖2 啟動過程中磁感應(yīng)強度在繞組端部的分布Fig.2 End-winding’s magnetic fl ux density distribution of starting state

圖3 穩(wěn)定運行過程中磁感應(yīng)強度在繞組端部的分布Fig.3 End-winding’s magnetic fl ux density distribution of steady state

1.2 端部磁場仿真分析

由于繞組端部結(jié)構(gòu)復(fù)雜，二維有限元分析難以達到實際的工程要求，因此通過軟件建立三維有限元模型可以精確地進行計算和仿真。通過FEM分析軟件計算磁感應(yīng)強度在電機定子繞組端部的分布情況，仿真結(jié)果如圖2和圖3所示。

圖2和圖3分別表示干式潛水電機在起動和穩(wěn)定運行過程中，A相電流達到最大值時磁感應(yīng)強度在繞組端部的分布情況。從圖中可以看出起動狀態(tài)下磁感應(yīng)強度比穩(wěn)定運行時磁感應(yīng)強度大,同時磁感應(yīng)強度在A相所屬繞組端部斜線段的分布比較集中，因此較之繞組端部其他部位，斜線段所受的電磁力比較大。當B相、C相電流達到最大值時，可得到相同的結(jié)論，即電磁力主要集中分布于繞組端部斜線段部位。

2 定子繞組端部電磁力的分析

2.1 繞組端部電磁力的數(shù)值計算

電機繞組端部的振動是在電磁力激發(fā)下進行的，同時繞組的疲勞、斷裂破壞等問題也均與電磁力有關(guān)系。因此，繞組端部電磁力的數(shù)值計算是干式潛水電機運行安全性和可靠性的重要前提。

如圖1所示，為求單根繞組所受電磁力，在繞組線棒的局部坐標系(ξ,η,γ)下，結(jié)合洛倫茲力公式可得

式中：J(Jξ,Jn,Jγ)為繞組傳導的電流密度矢量；B(Bξ,Bη,Bγ)為經(jīng)過坐標變換后的各段磁感應(yīng)強度；i,j,k分別為沿坐標ξ,η,γ,方向的單位矢量，而由傳導電流引起的電磁力密度矢量為F(Fξ,Fη,Fγ)。

電磁力主要集中分布于繞組端部斜線段，由傳導電流引起的電磁力密度分量在內(nèi)外層繞組斜線段的表達式為

內(nèi)層：

式中：Im為繞組電流的幅值；?0=0?,120?,240?分別對應(yīng)于電機三相電流的相角；Q1～Q4為坐標和B(Bξ,Bη,Bγ)的函數(shù)[13]；上標n和w分別表示繞組的內(nèi)層和外層。

2.2 端部繞組電磁力的仿真分析

針對干式潛水電機（1 000 kW/10 kV）的實際運行情況，結(jié)合式(8)～式(13)，通過數(shù)值分析方法，在實際算例的基礎(chǔ)上，給出繞組端部所受電磁力密度的仿真結(jié)果。

圖4～圖7分別給出了繞組端部斜線段內(nèi)、外層繞組所受磁力密度隨時間和坐標變化的三維圖。由圖可見內(nèi)外層電磁力密度隨時間作周期性變化，并且隨坐標ξ變化的過程中，三維曲線具有不同的峰值，對應(yīng)不同的最大正、負值。對斜線段內(nèi)層繞組而言，沿ξ的分布主要集中于0.13～0.19m范圍內(nèi)，靠近繞組鼻端。在ξ=0.16m處取得正向最大值為13.6 N/m，負向最大值為29.7 N/m。對斜線段外層繞組而言，沿ξ的分布主要集中于0.15～0.21m范圍內(nèi)，同樣靠近繞組鼻端，在ξ=0.18m處取得正向最大值為37.8 N/m，負向最大值為17.5 N/m。由此可見，電磁力主要集中分布于繞組端部斜線段上的0.13～0.21m范圍內(nèi)。

圖4 斜線段內(nèi)層電磁力密度隨時間t和坐標ξ的變化圖Fig.4 The electromagnetic forces density of inner slant partchange with time t and coordinate ξ

圖5 斜線段內(nèi)層電磁力密度隨時間t和坐標ξ的變化圖Fig.5 The electromagnetic forces density of inner slant partchange with time t and coordinate ξ

圖6 斜線段外層電磁力密度隨時間t和坐標ξ的變化圖Fig.6 The electromagnetic forces density of outer slantpart change with t time and coordinate ξ

圖7 斜線段外層電磁力密度隨時間t和坐標ξ的變化圖Fig.7 The electromagnetic forces density of outer slant part change with time t and coordinate ξ

3 端部繞組力學性能的改進

由以上結(jié)論可知，電磁力在斜線段上主要集中分布于0.13～0.21 m的范圍內(nèi)，并且在0.16 m和0.18m處取得極值。為減少電磁力對端部繞組的影響，提出了一種方法，通過改變直線段與斜線段夾角α（如圖1所示）進而改善端部繞組的力學性能。因此在數(shù)值算例(α=21.7?)的基礎(chǔ)上，分別給出了一個周期內(nèi)不同角度下對應(yīng)的電磁力密度在斜線段0.16 m和0.18 m處的變化規(guī)律。在數(shù)值算例(α=21.7?)的基礎(chǔ)上，工程要求允許范圍內(nèi)，每隔3?取一個角度，得到不同角度下電磁力密度隨時間的變化規(guī)律。

圖8和圖9表示內(nèi)層繞組電磁力密度在不同的角度下隨時間變化的曲線，由圖可知電磁力曲線近似地呈正弦規(guī)律分布。α=27.7?對應(yīng)的曲線幅值較之其它角度下的曲線幅值大，正向最大值為15.9 N/m，負向最大值為33.7 N/m。α=19.7?時對應(yīng)的曲線幅值較之其他角度下的幅值小，正向最大值為11.7 N/m，負向最大值為27.2 N/m。因此兩個角度下曲線對應(yīng)的幅值分別相差為4.2N/m和6.5N/m。圖10和圖11為外層繞組電磁力密度變化規(guī)律，且電磁力的分布曲線近似呈余弦規(guī)律分布。當α=27.7?時，曲線對應(yīng)的正向幅值較之其他角度下的幅值大，其大小為38.4 N/m。當α=24.7?時，曲線對應(yīng)的負向幅值較之其他角度下的幅值大，大小為14.3 N/m。但當α=19.7?時，曲線對應(yīng)的正向幅值和負向幅值都是最小，分別為28.3 N/m和8.8 N/m，因此正向幅值和負向幅值分別相差為10.1 N/m和5.5 N/m。

圖8 不同角度下內(nèi)層繞組電磁力密度在ξ=0.16m處隨時間的變化圖Fig.8 The forces density at ξ=0.16m of inner slant part change with time t corresponding different angles

圖9 不同角度下內(nèi)層繞組電磁力密度在ξ=0.16m隨時間t的變化圖Fig.9 The forces densityat ξ=0.16m of inner slant part change with time t corresponding different angles

圖10 不同角度下外層繞組電磁力密度在ξ=0.18m隨時間t的變化圖Fig.10 The forces density at ξ=0.18m of outer slant part change with time t corresponding different angles

圖11 不同角度下外層繞組電磁力密度在ξ=0.18m隨時間t的變化圖Fig.11 The forces densityat ξ=0.18m of outer slant part change with time t corresponding different angles

綜上所述，α=19.7?時對應(yīng)的電磁力密度幅值相對于其他角度下的電磁力密度幅值小，整體的變化趨勢相對來說比較平緩，由此可以減弱因電磁力作用而引起的端部電磁振動強度，進而降低端部繞組事故的發(fā)生率，延長電機使用壽命。

4 結(jié)論

本文針對干式潛水電機端部繞組易發(fā)生故障的實際工程問題，提出了一種通過改變繞組端部傾斜角α的方法改善端部繞組的力學性能。通過局部坐標變換法，詳細分析了電磁力密度在端部繞組斜線段上的分布規(guī)律，在此基礎(chǔ)上改變端部繞組錐度角α的大小，得出了不同角度下電磁力密度在斜線段上隨時間的分布圖。綜上所述，可得到以下結(jié)論：

1）電磁力密度在繞組端部斜線段上隨坐標和時間呈周期性變化，主要集中分布于端部繞組斜線段的后半段0.13～0.21 m范圍內(nèi)。

2）在工程允許的范圍內(nèi)，改變斜線段傾斜角對端部繞組力學性能進行優(yōu)化分析，通過對比計算得到當α=19.7?時，電磁力密度的幅值在數(shù)值算例的基礎(chǔ)上降低了9.5 N/m，達到了預(yù)期的效果。

該方法綜合運用了電磁場有限元、數(shù)值分析等計算手段，有效地分析了工程中的實際問題，具有較高的計算精確度，對干式潛水電機端部繞組的優(yōu)化設(shè)計具有較大的參考意義和實用價值。

[1]LIN Ranran,LAIHO Antti Nestori,HAAVISTO Ari,et al.End-winding vibrations caused by steady-state magnetic forces in an induction machine[J].IEEE Transactions on Magnetics,2010,46(7):2665-2674.

[2]MORI D,ISHIKAWA Takeo.Force and vibration analysis of induction motors[J].IEEE Transactions on Magnetic,2005,41(5):1948-1951.

[3]SHALLY D,FARRELL M,SULLIVAN K.Generator end winding vibration monitoring[C]//Proceedings of the Universities Power Engineering Conference,September 1-4,2008,Padora,Italy.2008:4651488.

[4]KIM Ki-Chan,LEE Hyung-Woo,CHUN Yon-Do,et al.Analysis of electromagnetic force distribution on end winding for motor reliance[J].IEEE Transactions on Magnetics,2005,41(10):4072-4074.

[5]SCOTT D J,SALON S J,KUSIK G L.Electromagnetic forces on the armature end windings of large turbine generatorsⅠsteady state conditions[J].IEEE Transactions on Magnetics,1981,PAS-100(11):4597-4603.

[6]SALON S J,SCOTT D J,KUSIK G L.Electromagnetic forces on the armature end windings of large turbine generatorsⅡTransientconditions[J].IEEETransactionson Power Apply,1983,34(7):14-19.

[7]LIN Ranran,ARKKIO Antero.3-D fi nite element analysis of magnetic forces on stator end-windings of an induction machine[J].IEEE Transactions on Magnetics,2008,44(11):4045-4048.

[8]WILLIAMSON S,ELLIS M R E.In fl uence of rotor currents on end-winding forces in cage motor[J].IEE Proceedings B:Electric Power Applications,1988,135(6):371-379.

[9]ALBANESE R,CALVANO F,MUT G D,et al.Coupled three dimensional numerical calculation of forces and stresses on the end windings of large turbo generators via integral formulation[J].IEEE Transactions on Magnetics,2012,48(2):875-878.

[10]胡宇達，邱家俊，卿光輝，等.大型汽輪發(fā)電機定子端部繞組整體結(jié)構(gòu)的電磁振動[J].中國電機工程學報，2003，23(7):93-98.

HU Yuda,QIU Jiajun,QING Guanghui,et al.Electromagnetic vibration of integrity end winding of large turbogenerator stator[J].Proceedings of the CSEE,2003,23(7):93-98.

[11]STANCHEVA Rumena D,IATCHEVA Ilona I.3-D Electromagnetic force distribution in the end region of turbogenerator[J].IEEETransactionsonMagnetics,2009,45(3):1000-1003.

[12]李京平,唐任遠.汽輪發(fā)電機定子端部繞組電動力研究[J].沈陽工業(yè)大學學報,1991,13(2):13-25.

LI Jingping,TANG Renyuan.The study of the electromagnetic forces on the end windings of stator in turbogenerators[J].Journal of Shenyang Polytechnic University,1991,13(2):13-25.

[13]湯蘊璆，梁艷萍.電機電磁場的分析與計算[M].北京：機械工業(yè)出版社，2010.

(編輯：于雙)

Research on electromagnetic force in the stator end winding of large induction motor

BAO Xiao-hua1，CHENG Xiao-wei1，F(xiàn)ANG Yong1，LU¨Qiang1，SHAN Li2
(1.School of Electrical Engineering and Automation,Hefei University of Technology,Hefei 230009，China;2.Hefei Hengda Jianghai Pump Co.,Ltd.,Hefei 231131,China)

According to the engineering problems produced by electromagnetic force of end winding in large induction motor,a method that changes taper angle of the end winding to improve the mechanical properties was proposed.Firstly,the magnetic fi eld expression in the region of end winding was derived by adopting the separation of variables method based on Maxwell equations,and then the magnetic fi eld distribution of the end winding was analyzed by Finite Element Method.Finally,combined with Lorentz rule and the local coordinated transformation,the electromagnetic force calculation formula of the end winding was derived,and the simulation of electromagnetic force was done by numerical recipes method.In combination with speci fi c examples of large dry submersible motor,different electromagnetic force distributions were derived by changing the taper angle of the stator end-winding.The simulation results show that the method improves the end winding’s mechanical properties effectively,prolongs the service life of motor and provides fundamental theory for the optimization design of stator end winding.

submersible motors; fi nite element method;taper angle; fl ux density of end-winding;electromagnetic force density

TM 46

A

1007–449X(2013)10–0027–06

2013–01–14

國家自然科學基金(51377039,51177033);中央高?；究蒲袠I(yè)務(wù)專項基金(2012HGZY0003)

鮑曉華(1972—),男,博士,教授,研究方向為電機設(shè)計理論和技術(shù)、電機電磁場理論分析及計算;程曉巍(1988—),男,碩士研究生,研究方向為電機電磁場分析與計算、電機絕緣系統(tǒng);

方 勇(1989—),男,碩士研究生,研究方向為電機電磁場計算;

呂 強(1988—),男,碩士研究生,研究方向為電機設(shè)計理論及電磁場分析;

單 麗(1964—),女,高級工程師,研究方向為電機制造與控制。

鮑曉華