基于有限元法的變壓器電磁振動噪聲分析

胡靜竹　劉滌塵　廖清芬　晏　陽　梁姍姍

(武漢大學(xué)電氣工程學(xué)院　武漢　430072)

?

基于有限元法的變壓器電磁振動噪聲分析

胡靜竹劉滌塵廖清芬晏陽梁姍姍

(武漢大學(xué)電氣工程學(xué)院武漢430072)

變壓器噪聲是變電站主要噪聲來源之一，研究變壓器噪聲對合理評估變電站噪聲水平及對變壓器降噪分析具有重要意義。從變壓器噪聲產(chǎn)生機理入手，分析變壓器鐵心和繞組的電磁振動噪聲，建立了其電磁-結(jié)構(gòu)-聲場有限元模型。通過瞬態(tài)電磁場分析，并基于虛位移法得到鐵心和繞組所受電磁力的時域波形，采用FFT變換對結(jié)果進(jìn)行后處理獲取電磁力的主要諧波分量大小，將其作為結(jié)構(gòu)諧響應(yīng)分析的激勵源，通過頻域內(nèi)的振動分析得到鐵心和繞組表面各節(jié)點的振動位移，并將其作為變壓器聲場分析的邊界條件，進(jìn)一步求解變壓器聲場模型，分析得到變壓器周圍空間場點在噪聲集中頻率100 Hz和200 Hz上的聲壓級，并與實測值進(jìn)行對比分析，驗證了該模型的可行性，可為變壓器噪聲預(yù)測提供理論依據(jù)和計算方法。

變壓器噪聲電磁-結(jié)構(gòu)-聲場模型有限元法FFT電磁力諧波

0　引言

變電站噪聲環(huán)境問題隨著電力系統(tǒng)的發(fā)展日益嚴(yán)重，而變壓器作為變電站主噪聲源設(shè)備是研究的重點。變壓器的噪聲主要來源于電磁振動，而從噪聲產(chǎn)生機理上進(jìn)行研究不僅可在變壓器初步設(shè)計階段預(yù)估其噪聲值，還可為有效降低變壓器噪聲提供依據(jù)。

對于變壓器電磁振動的研究國外最早開始于20世紀(jì)20年代，主要是一些大型電力變壓器制造公司和相關(guān)研究機構(gòu)[1]。文獻(xiàn)[2]研究了變壓器振動鐵心在磁致伸縮下的強迫振動以及在諧波頻率下的共振，會導(dǎo)致變壓器噪聲值的變大。該文還對變壓器鐵心進(jìn)行了大量模擬，分析出鐵心的模態(tài)共振頻率，從而使變壓器避開共振，減少變壓器的噪聲[3]。文獻(xiàn)[4]采用有限元數(shù)值分析方法，將大型電力變壓器的電磁場與結(jié)構(gòu)力場進(jìn)行耦合計算，得出變壓器鐵心的振動情況。國內(nèi)對大型電力變壓器噪聲的研究主要集中在噪聲測量分析[5-7]、噪聲控制技術(shù)研究[8，9]、采用經(jīng)驗公式[10-12]計算變壓器輻射噪聲等。文獻(xiàn)[13-15]將電磁場理論與結(jié)構(gòu)力學(xué)理論相結(jié)合，建立了電力變壓器鐵心電磁振動的數(shù)學(xué)模型，就目前來說，振動數(shù)學(xué)模型將是研究振動計算的一種發(fā)展趨勢。

本文建立了變壓器鐵心和繞組的電磁-機械-聲場有限元模型，在Ansys workbench中建立了變壓器繞組、鐵心和油箱的有限元模型，基于麥克斯韋電磁理論與虛位移法，對變壓器鐵心和繞組進(jìn)行瞬態(tài)電磁場分析得到其磁場分布與電磁振動受力情況。并在此基礎(chǔ)上，對鐵心和繞組進(jìn)行結(jié)構(gòu)諧響應(yīng)分析得到結(jié)構(gòu)表面振動位移，進(jìn)一步在聲學(xué)軟件LMS中進(jìn)行聲場分析得到變壓器周圍噪聲分布情況。

1　電磁-結(jié)構(gòu)-聲場分析模型基礎(chǔ)

1.1電磁分析基礎(chǔ)

基于Maxwell方程且考慮硅鋼片磁致伸縮應(yīng)力影響，忽略鐵心的渦流效應(yīng)，可得變壓器鐵心區(qū)域瞬態(tài)電磁場方程為

(1)

繞組區(qū)域求解方程為

(2)

(3)

磁場儲能通過磁場的勢函數(shù)計算獲得，在有限元計算中的表達(dá)式為

(4)

求解得到磁場密度和儲能可進(jìn)一步計算鐵心和繞組所受的電磁力，在計算變壓器鐵心和繞組所受電磁力時采用虛功法，瞬時電磁力為

(5)

式中，W(s，i)為系統(tǒng)的磁場儲能，J；i為建立磁場的電流，A。

1.2結(jié)構(gòu)分析基礎(chǔ)

通過瞬態(tài)電磁場分析得到鐵心和繞組所受的電磁力為隨時間變化的瞬態(tài)值，對其進(jìn)行離散傅里葉變換，可得到電磁力的各諧波分量幅值和相角大小，將其作為簡諧激勵源，然后進(jìn)行穩(wěn)態(tài)結(jié)構(gòu)諧響應(yīng)振動分析。

諧響應(yīng)分析的運動方程為

(-ω2M+iωC+K)(u1+iu2)=(F1+iF2)

(6)

式中，ω為簡諧激勵的角頻率，Hz；M為質(zhì)量矩陣，kg；K為剛度矩陣，N/m；C為阻尼矩陣，N/(m/s)；u1、u2分別為振動位移的實部與虛部，m；F1、F2分別為結(jié)構(gòu)受力的實部與虛部，N。

1.3聲場分析基礎(chǔ)

采用聲學(xué)有限元法求解聲學(xué)Helmholtz方程來計算聲場。通過聲波的連續(xù)方程、運動方程和物態(tài)方程可推導(dǎo)得到Helmholtz波動方程[16]，進(jìn)一步通過傅里葉變換可得均勻流體中傳播的基本聲學(xué)方程頻域形式為

(7)

式中，k為波數(shù)，k=ω/c=2πf/c，m-1，其中c為聲波在流體中的聲速(m/s)；ρ0為流體密度，kg/m3；q0為外部作用于流體的質(zhì)量源。

計算變壓器輻射聲場屬于邊界封閉的外聲場求解問題，邊界條件包括兩部分封閉的邊界和無限遠(yuǎn)處的邊界。在封閉的邊界上應(yīng)滿足速度邊界條件，即將結(jié)構(gòu)表面振動響應(yīng)結(jié)果導(dǎo)入Virtual.Lab Acoustics作為聲學(xué)分析的邊界條件，則聲學(xué)有限元系統(tǒng)方程形式為

(K+jωC-ω2M)·pi=Vni

(8)

式中，M為聲學(xué)質(zhì)量矩陣，kg；K為聲學(xué)剛度矩陣，N/m；C為聲學(xué)阻尼矩陣，N/(m/s)；pi為節(jié)點聲壓，Pa；Vni為節(jié)點邊界條件，m/s。

1.4電磁-結(jié)構(gòu)-聲場聯(lián)合分析基礎(chǔ)

瞬態(tài)電磁分析與機械振動計算均在ANSYS中完成，聲場計算在LMS Virtual.Lab Acoustics中進(jìn)行，電磁-機械-聲場分析步驟流程如圖1所示。

圖1　電磁振動噪聲計算步驟流程Fig.1　Process of electromagnetic vibration noise calculation steps

電磁瞬態(tài)分析是振動分析的基礎(chǔ)，在ANSYS MAXWELL中進(jìn)行瞬態(tài)仿真計算時，每隔0.5 ms記錄一組結(jié)果數(shù)據(jù)，求解完后將所有電磁力計算結(jié)果導(dǎo)出為.txt格式，然后導(dǎo)入Matlab進(jìn)行FFT變換，得到頻域數(shù)據(jù)作為結(jié)構(gòu)分析基礎(chǔ)。結(jié)構(gòu)分析時電磁力數(shù)據(jù)通過命令流在ANSYS中直接讀入，并將求解得到的結(jié)構(gòu)表面振動位移數(shù)據(jù)寫入結(jié)果文件中。ANSYS計算結(jié)果保存為*.rst文件，模型網(wǎng)格保存為.cbd文件，可直接導(dǎo)入聲學(xué)計算軟件LMS中進(jìn)行聲場分析。

2　變壓器電磁振動噪聲有限元分析

對110 kV油浸式電力三柱式雙繞組變壓器進(jìn)行仿真建模分析，其基本參數(shù)為：高壓繞組額定電流為165.3 A，低壓繞組額定電流為1 818.7 A，高壓繞組每相匝數(shù)為737匝，低壓繞組每相匝數(shù)為116匝，高壓線圈直徑為948～1 111 mm，低壓線圈直徑為598～710 mm，線圈高1 216 mm，鐵心直徑為565 mm，鐵心軛高為540 mm，油箱體積為5 m×2 m×3 m。建立變壓器的三維實體模型，包括變壓器鐵心、繞組和變壓器油，繞組建模呈圓柱筒形狀，計算時高低壓繞組之間的撐條通過采用施加約束條件模擬，如圖2a所示。有限元法在求解模型時，需將模型進(jìn)行網(wǎng)格劃分，本文采用solid45結(jié)構(gòu)單元對變壓器鐵心和繞組進(jìn)行劃分，采用fluid30流體單元對變壓器油進(jìn)行劃分，網(wǎng)格模型如圖2b所示，模型劃分后共有66 989個節(jié)點。

圖2　變壓器有限元計算模型Fig.2　Finite element calculation model of transformer

2.1電磁場分析

將變壓器鐵心設(shè)置為硅鋼片材料屬性，疊積系數(shù)取為0.95，分別定義軋制方向和垂直于軋制方向的非線性B-H曲線數(shù)據(jù)，為考慮磁致伸縮效應(yīng)，根據(jù)相對磁導(dǎo)率與應(yīng)力的電磁關(guān)系[17](見式(9))以及應(yīng)力與磁感應(yīng)強度的關(guān)系(見式(10))采用插值迭代法對其進(jìn)行修正。

(9)

Bσ=(μ+Δμ)H+λσ

(10)

式中，μ為磁導(dǎo)率；λ為磁致伸縮系數(shù)；λm為磁飽和情況下的磁致伸縮系數(shù)；Bm為飽和磁感應(yīng)強度，T；σ為應(yīng)力，N/m2；Bσ為應(yīng)力作用下的磁感應(yīng)強度，T。

在堆疊方向取相對磁導(dǎo)率為1，忽略鐵心的渦流效應(yīng)，將相對電導(dǎo)率取為1，變壓器繞組和變壓器油均取相對磁導(dǎo)率為1。給繞組賦予銅的材料屬性，在三相低壓繞組截面上分別施加相位相差120°的正弦電流，電流平均分布在截面上。分別設(shè)置高低壓繞組的匝數(shù)。由于大型變壓器油箱壁內(nèi)側(cè)都有磁屏蔽，可將油箱磁屏蔽體看成邊界面，對整個模型外側(cè)施加狄里克萊邊界條件，即磁力線平行于整個模型外表面。

變壓器繞組和鐵心及變壓器油參數(shù)見表1。

表1　變壓器結(jié)構(gòu)材料參數(shù)Tab.1　Material properties of transformer

對變壓器進(jìn)行瞬態(tài)電磁場分析，對每相繞組和鐵心設(shè)置力參數(shù)，時間設(shè)置為0～150 ms，每隔0.5 ms計算一個值，可得到每相繞組和鐵心的整體結(jié)構(gòu)受力分別在X、Y、Z方向的分力。

加正弦電流激勵源進(jìn)行瞬態(tài)分析時，經(jīng)過50 ms波動后達(dá)到穩(wěn)定狀態(tài)，鐵心和繞組(以A相低壓繞組為例)受電磁力波形如圖3和圖4所示。

圖3　鐵心瞬態(tài)受力波形Fig.3　The transient stress wave of transformer core

圖4　A相低壓繞組瞬態(tài)受力波形Fig.4　Transient stress wave of low voltage winding of phase A

圖3和圖4表明電磁力瞬態(tài)波形均為周期波，周期為正弦電流源的一半，T=10 ms。各瞬態(tài)電磁力波形為非正弦，含有不同階次的諧波。

2.2結(jié)構(gòu)振動諧響應(yīng)分析

由于聲場分析在頻域內(nèi)進(jìn)行，需要變壓器鐵心和繞組在頻域內(nèi)的振動位移數(shù)據(jù)，因此對變壓器鐵心和繞組進(jìn)行諧響應(yīng)分析，即在頻域內(nèi)分析變壓器鐵心和繞組振動情況。而電磁分析得到的電磁力為時域值，因此對繞組和鐵心所受電磁力進(jìn)行FFT變換，得到各主要諧波的幅值與相角，由于電磁力在不同方向上分力相差較大，在進(jìn)行振動分析時，忽略幅值較小的力。占主導(dǎo)作用的分力的FFT變換結(jié)果如圖5所示(繞組受力FFT分析結(jié)果以三相低壓繞組為例)，進(jìn)行FFT變換時將結(jié)果中幅值較小的諧波濾去。

圖5　瞬態(tài)電磁力FFT變換Fig.5　FFT transformation of transient electromagnetic force

可見變壓器鐵心和繞組所受電磁力諧波分量在不同頻率上的大小和相位均不同，諧波主要集中在基頻的倍頻帶100～500 Hz上，變壓器噪聲能量集中在低頻段。而且100 Hz和200 Hz最為突出，因此本文重點分析了變壓器輻射噪聲中100 Hz和200 Hz噪聲分量的大小。

進(jìn)行諧響應(yīng)分析，分別將X、Y、Z方向正弦電磁力諧波平均施加到繞組和鐵心的每個節(jié)點上，施加的力的參數(shù)為100 Hz和200 Hz的幅值和相角，并對每個高低壓繞組12個撐條截面的位置和鐵心底面施加全約束，即約束節(jié)點X、Y、Z三個方向的位移，使UX、UY和UZ均為零。分別求解100 Hz和200 Hz的振動情況。依據(jù)計算聲學(xué)基本原理[18]，鐵心和繞組上每個與空氣介質(zhì)接觸的面都獨立地輻射噪聲，因此提取繞組和鐵心所有表面節(jié)點在不同頻率下的振動位移作為后續(xù)聲場分析的基礎(chǔ)。

鐵心和繞組的振動位移如圖6和圖7所示。

圖6　變壓器繞組振動位移Fig.6　Displacements of transformer windings

圖7　變壓器鐵心振動位移Fig.7　Displacements of transformer windings

2.3聲場分析

采用有限元法和PML聲吸收邊界條件進(jìn)行聲場計算，即在聲學(xué)輻射邊界增加幾層網(wǎng)格吸收聲學(xué)量，在增加的匹配層外圍邊界上設(shè)置聲學(xué)量為零。結(jié)構(gòu)網(wǎng)格和聲學(xué)網(wǎng)格均在ANSYS中建立并導(dǎo)入LMS，將結(jié)構(gòu)分析計算得到的變壓器鐵心和繞組表面節(jié)點100 Hz和200 Hz的振動位移數(shù)據(jù)分別導(dǎo)入LMS Virtual.Lab Acoustics中，并將結(jié)構(gòu)網(wǎng)格上的振動數(shù)據(jù)按照加權(quán)平均的關(guān)系映射到聲學(xué)網(wǎng)格上作為聲場分析的邊界條件，給聲場網(wǎng)格賦予變壓器油的材料屬性，分別計算變壓器前后和左右兩個面的輻射噪聲，將場點布置在距離變壓器油箱外1 m的平面上，位于變壓器正面的場點平面位置示意圖如圖8所示。

圖8　場點位置示意圖Fig.8　Position of field point

通過對聲場有限元的求解和聲場分布的計算得到變壓器外1 m處平面上場點的聲壓云圖如圖9所示，由于模型結(jié)構(gòu)的對稱性，變壓器前后兩個面和左右兩個側(cè)面輻射的噪聲分布是一致的。

圖9　變壓器周圍1 m處輻射聲壓大小Fig.9　Noise level at 1 m away from transformer

可見110 kV變壓器正面輻射噪聲比側(cè)面大，100 Hz噪聲在正面主要分布在62～67 dB，側(cè)面主要分布在57～63 dB左右，200 Hz噪聲在正面主要分布在73～77 dB，側(cè)面主要分布在60～65 dB。

2.4計算結(jié)果與實測值對比分析

對110 kV變電站變壓器噪聲頻譜進(jìn)行了測量，測量時變壓器在額定功率下正常運行，風(fēng)扇關(guān)閉。在變壓器箱體四周分別布置測點(1、2、3、4)，見圖10。

圖10　變壓器噪聲測量布點Fig.10　Layout of transformer noise measuring points

聲強探頭與變壓器基準(zhǔn)發(fā)射面垂直距離為1 m，測量高度為1.3 m，采用聲強法測量，最大程度減小防火墻對噪聲的影響，準(zhǔn)確測得變壓器本體噪聲大小[19]。得到變壓器噪聲1/3倍頻帶中心頻率上的噪聲值，4個方位上各點在100 Hz和200 Hz的噪聲值見表2。

表2　110 kV變壓器噪聲測量值Tab.2　Measured noise level of 110 kV transformer

通過對比分析可知仿真計算值與實測值較為接近，誤差在3～5 dB，驗證了此模型在計算變壓器輻射聲場時的準(zhǔn)確性。在變壓器設(shè)計初始階段可為預(yù)測變壓器的噪聲分布情況提供有效的分析方法，另外可為變壓器電磁噪聲降噪措施的降噪性能評估提供理論依據(jù)和計算方法。

3　結(jié)論

本文對110 kV油浸式變壓器進(jìn)行了磁場-機械-聲場的三維有限元分析，分析了變壓器鐵心和繞組的受力情況與變壓器周圍電磁振動噪聲的分布，并與實測結(jié)果進(jìn)行對比分析，得到以下結(jié)論：

1)建立了變壓器鐵心、繞組和變壓器絕緣油的有限元模型，對其進(jìn)行瞬態(tài)電磁場分析得到變壓器鐵心和繞組在交變電流作用下所受的交變電磁力，并進(jìn)行了FFT后處理變換得到電磁力的主要諧波分量，其分布在頻率100～500 Hz上，且100 Hz和200 Hz最為突出。

2)分別分析了變壓器鐵心和繞組在100 Hz和200 Hz的諧波電磁力激勵下的振動情況，并在此基礎(chǔ)上對變壓器進(jìn)行聲場分析，得到變壓器周圍空間噪聲分布情況，并與實測值進(jìn)行對比，仿真計算值與測量值相符，證實了該分析模型的正確可靠性。此模型可在設(shè)計初期分析變壓器噪聲分布情況，可為變壓器降噪與設(shè)計提供依據(jù)。

[1]顧曉安，沈榮瀛，徐基泰.國外變壓器噪聲研究動向[J].變壓器，2002，39(6)：33-38.

Gu Xiaoan，Shen Rongying，Xu Jitai.Trend of transformer noise research abroad[J].Transformer，2002，39(6)：33-38.

[2]Fahnoe H.A study of sound levels of transformers[J].Electrical Engineers，1941，60(6)：748-749.

[3]Masti R S，Heylen W，Leuven K U，et al.On the influence of core laminations upon power transformer noise[C]//Proceedings of the International Conference on Noise and Vibration Engineering，Leuven，Belgium，2004：3851-3862.

[4]Foster S L，Reiplinger E.Characteristics and control of transformer sound[J].IEEE Transactionson Power Systems，1981，PAS-100(3)：1072-1077.

[5]韓輝，吳桂芳，瞿雪弟，等.我國±500 kV換流站設(shè)備可聽噪聲的測量分析及降噪措施[J].電網(wǎng)技術(shù)，2008，32(2)：38- 41.Han Hui，Wu Guifang，Zhai Xuedi，et al.Measurement and analysis of audible noise bringing about by equipments in domestic ±500 kV converter stations and noise reduction measures[J].Power System Technology，2008，32(2)：38- 41.

[6]周兵，裴春明，倪園，等.特高壓交流變電站噪聲測量與分析[J].高電壓技術(shù)，2013，39(6)：1447-1453.Zhou Bing，Pei Chunming，Ni Yuan，et al.Measurement and analysis of noise in UHV AC substation[J].High Voltage Engineering，2013，39(6)：1447-1453.

[7]朱華，劉佳，朱子斌，等.500 kV 輸變電工程可聽噪聲測量及分析研究[J].電氣應(yīng)用，2011，30(13)：42- 45.

Zhu Hua，Liu Jia，Zhu Zibin，et al.Analysis of the audible noise of 500 kV power transmission project[J].Electrotechnical Application，2011，30(13)：42- 45.

[8]周建國，李莉華，杜茵，等.變電站、換流站和輸電線路噪聲及其治理技術(shù)[J].中國電力，2009，42(3)：75-78.

Zhou Jianguo，Li Lihua，Du Yin，et al.Review on methods of noise reducing for equipment in substations converter stations and transmission lines[J].Electric Power，2009，42(3)：75-78.

[9]楊一鳴，章旭雯.特高壓直流換流站設(shè)備的降噪措施[J].高電壓技術(shù)，2006，32(9)：149-152.

Yang Yiming，Zhang Xuwen.Study on noise reduce for equipment at UHVDC converter station[J].High Voltage Engineering，2006，32(9)：149-152.

[10]陳玉紅，王立忠，洛君婷，等.變壓器噪聲的有限元輔助算法[J].變壓器，2013，50(6)：5-9.

Chen Yuhong，Wang Lizhong，Luo Junting，et al.Finite element aided algorithms of transformer noise[J].Transformer，2013，50(6)：5-9.

[11]徐祿文，劉小玲.變電站環(huán)境噪聲三維空間衰減模型及算法研究[J].中國電機工程學(xué)報，2012，32(13)：175-180.

Xu Luwen，Liu Xiaoling.Study on the three dimension attenuated model and the algorithm of environmental noise in substations[J].Proceedings of the CSEE，2012，32(13)：175-180.

[12]吳高強，程勝高，黃磊，等.戶外220kV變電站噪聲環(huán)境影響預(yù)測研究[J].噪聲與振動控制，2007，6(3)：135-137.

Wu Gaoqiang，Cheng Shenggao，Huang Lei，et al.Prediction on noise of 220 kV outdoor substation to environmental infection[J].Noise and Vibration Control，2007，6(3)：135-137.

[13]祝麗花，楊慶新，閆榮格，等.考慮磁致伸縮效應(yīng)電力變壓器振動噪聲的研究[J].電工技術(shù)學(xué)報，2013，28(4)：1-6，19.

Zhu Lihua，Yang Qingxin，Yan Rongge，et al.Research on vibration and noise of power transformer cores including magnetostriction effects[J].Transactions of China Electrotechnical Society，2013，28(4)：1-6，19.

[14]王志敏，顧文業(yè)，顧曉安，等.大型電力變壓器鐵心電磁振動數(shù)學(xué)模型[J].變壓器，2004，41(6)：1-6.

Wang Zhimin，Gu Yewen，Gu Xiaoan，et al.Research on electromagnetic vibration mathematical model of large power transformer core[J].Transformer，2004，41(6)：1-6.

[15]白保東，劉闖，王佳音.直流偏磁下變壓器振動的仿真及實驗研究[J].電工技術(shù)學(xué)報，2013，28(增2)：427- 433.

Bai Baodong，Liu Chuang，Wang Jiayin.Simulation and Experimental research of transformer vibration under DC magnetic bias[J].Transactions of China Electrotechnical Society，2013，28(S2)：427- 433.

[16]李增剛，詹福良.Virtual.Lab Acoustics聲學(xué)仿真計算高級應(yīng)用實例[M].北京：國防工業(yè)出版社，2010.

[17]王社良，王威，蘇三慶，等.鐵磁材料相對磁導(dǎo)率變化與應(yīng)力關(guān)系的磁力學(xué)模型[J].西安科技大學(xué)學(xué)報，2005，25(3)：288-291.

Wang Sheliang，Wang Wei，Su Sanqing，et al.A magneto-mechanical model on differential permeability and stress of ferromagnetic material[J].Journal of Xi’an University of Science and Technology，2005，25(3)：288-291.

[18]Kubiak W，Witczak P.Vibration analysis of small power transformer[J].The International Journal for Computation and Mathematics in Electrical and Electronic Engineering，2010，29(4)：1116-1124.

[19]GBT 1094.10—2003，電力變壓器-第10部分：聲級測定[S].

Analysis of Transformer Electromagnetic Vibration Noise Based on Finite Element Method

Hu JingzhuLiu DichenLiao QingfenYan YangLiang Shanshan

(School of Electrical EngineeringWuhan UniversityWuhan430072China)

The noise of transformer is one of the main noise sources in substations.The analysis of the transformer noise is of vital importance for noise prediction and control.In this paper，a magneto-mechanical-acoustic coupled finite element model is presented through studying the transformer noise generating mechanism and analyzing the electromagnetic vibration noise of the transformer core and winding.First，the flux density is calculated by transient electromagnetic field analysis and the electromagnetic forces of the core and winding are calculated using the virtual displacement method.Then，F(xiàn)FT analysis is carried out to obtain the main harmonics of the electromagnetic forces which are set as the excitation sources for structure vibration harmonic response analysis.By vibration analysis in the frequency domain，the vibration displacements of each point on the surface of transformer core and winding are obtained and the results are set as the boundary conditions for acoustic field analysis.Finally，the sound pressure levels of noise of 100 Hz and 200 Hz around the transformer space are calculated by solving the acoustic model and compared to the measured data to confirm the validity of the proposed method.Results show that the model is available for transformer noise prediction.

Transformer noise，magneto-mechanical-acoustic coupled model，finite element method，F(xiàn)FT，electromagnetic force harmonic

2015-04-12改稿日期2015-08-04

胡靜竹女，1990年生，碩士研究生，研究方向為電力系統(tǒng)運行控制與電磁兼容。

E-mail：jzhu_5@163.com(通信作者)

劉滌塵男，1953年生，教授，研究方向為電力系統(tǒng)分析與運行控制。

E-mail：Dcliu@whu.edu.cn

國家電網(wǎng)公司大電網(wǎng)重大專項資助項目(GYW17201300115)。