1 引言

隨著市場需求嚴(yán)苛程度不斷提高，變壓器容量增大，其運行穩(wěn)定性成為了用戶關(guān)注度極高的問題。變壓器性能包括散熱、噪聲、振動、抗短路能力等眾多因素，變壓器作為電站主要設(shè)備之一，并且是變電站主要噪聲源設(shè)備是研究的重點，因此變壓器的噪聲問題一直是設(shè)計人員關(guān)注的重點。本文中根據(jù)GB/T 1094.10變壓器聲級測定標(biāo)準(zhǔn)，結(jié)合變壓器額定負載運行工況，基于ANSYS Workbench平臺實現(xiàn)了變壓器噪聲分析，從而在噪聲產(chǎn)生機理上進行深入研究，不僅可以在變壓器設(shè)計階段預(yù)估噪聲值，還可以為有效降低變壓器噪聲提供科學(xué)依據(jù)。

2 噪聲分析理論基礎(chǔ)

2.1電磁分析基礎(chǔ)

電磁場理論由麥克斯韋方程組(如下圖所示)來描述。求解方法上，數(shù)值法優(yōu)于解析法，近年來電磁場數(shù)值解法在工程及科學(xué)研究上的應(yīng)用也越來越廣泛和高效。電磁場的數(shù)值分析和計算通常歸結(jié)為求微分方程的解，對于偏微分方程，輔助邊界條件和初始條件即可獲得方程的定解。

ANSYS Maxwell 采用有限元法，將求解區(qū)域離散化為”單元“，采用Maxwell方程進行求解。

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

通過電磁場分析得到鐵芯和繞組所受的電磁力分布，對其進行傅里葉變換，可以得到電磁力各諧波分量的幅值和相位角大小，將其作為簡諧激勵源，進行結(jié)構(gòu)的諧響應(yīng)分析。諧響應(yīng)分析的運動控制方程為：

其中假設(shè)F和u做簡諧變化，則：

2.3 噪聲分析基礎(chǔ)

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

計算變壓器聲場分析需要將結(jié)構(gòu)表面的振動速度導(dǎo)入聲學(xué)分析中作為邊界條件，聲學(xué)有限元系統(tǒng)方程形式為：

2.4 耦合分析流程

本次分析首先在MAXWELL進行電磁場分析，求解完成后，對電磁力進行FFT變換，在workbench平臺利用耦合功能，將其導(dǎo)入Mechanical進行簡諧振動分析，得到質(zhì)點振動速度，再將其導(dǎo)入ANSYS Acoustics聲學(xué)仿真模塊，求解聲壓波動方程，進行聲場分析，得到最后的噪聲計算結(jié)果，并根據(jù)GB/T 1094.10進行評定。

Figure.基于ANSYS Workbench的聲學(xué)仿真耦合流程

3 干式變壓器振動噪聲分析

3.1電磁場分析

將變壓器的電磁模型導(dǎo)入Maxwell,給定鐵芯、繞組的材料，設(shè)定好額定工況的激勵、邊界條件、求解參數(shù)，即可進行求解。設(shè)定好的繞組激勵如下圖所示：

① 設(shè)定鐵芯、繞組材料：

Figure.材料設(shè)定

② 施加激勵、求解計算：

Figure.激勵加載&求解設(shè)置

③ 后處理：

Figure.后處理設(shè)置

Figure. 電磁力密度

3.2 結(jié)構(gòu)分析

在mechanical中進行分析前，首先根據(jù)提供的材料在Engineer Data中輸入材料數(shù)據(jù)，由于諧響應(yīng)分析是線性分析類型，并且變壓器結(jié)構(gòu)在實際工作中也不允許超出屈服強度，因此此處以線彈性材料進行簡化輸入。網(wǎng)格劃分過程中，實體單元以四面體、六面體混合。根據(jù)實際工作，掃頻范圍設(shè)置為0~1000Hz。加載時，根據(jù)變壓器實際安裝位置，將下部的底座框架施加固定約束。具體操作如下：

① 網(wǎng)格劃分：針對模型不同部件，Mesh下插入Body Sizing，指定尺寸，生成網(wǎng)格。

② 邊界條件：根據(jù)實際工作情況，將底部進行全約束。在Harmonic Response處右鍵insert插入fixed support

③ 分析設(shè)置：此處根據(jù)前述分析，將頻率區(qū)間設(shè)置為0~1000Hz

Figure. 分析設(shè)置

④ 導(dǎo)入電磁力：在Import Load處，鼠標(biāo)右鍵Insert，選擇Surface Force Density，選擇需要導(dǎo)入電磁力的部件，Surface Force Density右鍵選擇import Load，即可導(dǎo)入。

⑤ 后處理：

3.3噪聲分析

噪聲分析利用ANSYS專業(yè)噪聲仿真模塊Acoustics。噪聲分析需要輸入聲音在介質(zhì)中的傳播速度及介質(zhì)密度等參數(shù)，此處介質(zhì)為空氣，在Engineer Data中輸入相應(yīng)數(shù)據(jù)即可。噪聲分析由于主要分析聲音在介質(zhì)中傳播現(xiàn)象，因此需要設(shè)置空氣域。由于變壓器與空氣接觸部分幾何復(fù)雜，因此對空氣域采用四面體網(wǎng)格劃分方式?；贏NSYS Workbench耦合平臺，將上一步諧響應(yīng)分析計算得到的質(zhì)點振動速度導(dǎo)入噪聲分析中，作為激勵源。通過計算可以得到不同頻率下的聲壓情況，由于輸入正弦激勵，頻率為50Hz，而一次交流過程中會有兩次信號達到峰值，因此振動分析的基礎(chǔ)頻率為100Hz。因此可以查看100,200,300等倍頻噪聲情況，此分析中僅截止到1000Hz。計算完成后，根據(jù)GB/T 1094.10變壓器聲級測定標(biāo)準(zhǔn)，后處理中提取相關(guān)輪廓線處A計權(quán)聲壓，并計算平均值，得到最終結(jié)果。

① 模型處理：進行聲場分析，首先需要建立空氣域，在Design Modeler中利用Enslosure功能可以插入空氣域，同時指定空氣域大小即可。

② 網(wǎng)格劃分：由于空氣域形狀復(fù)雜，此處以四面體方式進行網(wǎng)格劃分，此類特征的幾何模型適合采用Patch Independent算法進行網(wǎng)格劃分。此處Max Element Size指定為250mm。

③ 邊界條件：右鍵單擊Import Load選擇Insert，插入Velocity，插入諧響應(yīng)分析中計算得到的質(zhì)點振動速度作為聲場分析激勵。

④ 分析設(shè)置：進行聲場分析前，需要選擇聲場區(qū)域。在Harmonic Acoustics處右鍵單擊，選擇Insert，選擇Physics Region,選擇我們繪制的聲場區(qū)域。

⑤ 后處理：求解計算完成后，在Solution右鍵單擊，選擇Insert，選擇Acoustics，選擇我們關(guān)心的結(jié)果即可。

Figure. 不同頻率下聲壓變化曲線(前后面最大聲壓)

通過上述曲線，發(fā)現(xiàn)前后面聲壓最大發(fā)生在400Hz時。

通過上述曲線，發(fā)現(xiàn)側(cè)面聲壓最大時為300Hz。

Figure. 300Hz時側(cè)面聲壓分布

根據(jù)實際試驗要求，提取輪廓線處的聲壓，并取平均值。

表格13 各點聲壓值

結(jié)果說明

①通過噪聲分析，發(fā)現(xiàn)變壓器在工作時，前后面的聲壓分布趨勢基本一致，側(cè)面的聲壓分布趨勢基本一致，最大值略有差異。

②通過噪聲分析，發(fā)現(xiàn)該變壓前后面的最大A計權(quán)聲壓為58dB，側(cè)面最大A計權(quán)聲壓為50dB。

③通過噪聲分析后處理，300Hz平均聲壓為50.4dB，400Hz平均聲壓為得到平均為50.8dB。

4 總結(jié)

本文通過基于ANSYS Workbench平臺的干式變壓器振動噪聲仿真，實現(xiàn)了在產(chǎn)品設(shè)計階段對其噪聲值進行預(yù)估的完整流程，可以幫助企業(yè)在探究變壓器噪聲的機理上，對產(chǎn)品及時做出改進，響應(yīng)市場，提高競爭力。