|
1 引言 隨著市場需求嚴苛程度不斷提高，變壓器容量增大，其運行穩(wěn)定性成為了用戶關注度極高的問題。變壓器性能包括散熱、噪聲、振動、抗短路能力等眾多因素，變壓器作為電站主要設備之一，并且是變電站主要噪聲源設備是研究的重點，因此變壓器的噪聲問題一直是設計人員關注的重點。本文中根據(jù)GB/T 1094.10變壓器聲級測定標準，結(jié)合變壓器額定負載運行工況，基于ANSYS Workbench平臺實現(xiàn)了變壓器噪聲分析，從而在噪聲產(chǎn)生機理上進行深入研究，不僅可以在變壓器設計階段預估噪聲值，還可以為有效降低變壓器噪聲提供科學依據(jù)。 2 噪聲分析理論基礎 2.1電磁分析基礎 電磁場理論由麥克斯韋方程組（如下圖所示）來描述。求解方法上，數(shù)值法優(yōu)于解析法，近年來電磁場數(shù)值解法在工程及科學研究上的應用也越來越廣泛和高效。電磁場的數(shù)值分析和計算通常歸結(jié)為求微分方程的解，對于偏微分方程，輔助邊界條件和初始條件即可獲得方程的定解。 ANSYS Maxwell 采用有限元法，將求解區(qū)域離散化為”單元“，采用Maxwell方程進行求解。 2.2 結(jié)構(gòu)分析基礎 通過電磁場分析得到鐵芯和繞組所受的電磁力分布，對其進行傅里葉變換，可以得到電磁力各諧波分量的幅值和相位角大小，將其作為簡諧激勵源，進行結(jié)構(gòu)的諧響應分析。諧響應分析的運動控制方程為： 其中假設F和u做簡諧變化，則： 采用聲學有限元法求解聲學Helmholtz方程來計算聲場。通過聲波的連續(xù)方程、運動方程、物態(tài)方程可以推導得到Helmholtz波動方程，進一步通過傅里葉變換可以得到均勻流體中傳播的基本聲學方程頻域形式為： 計算變壓器聲場分析需要將結(jié)構(gòu)表面的振動速度導入聲學分析中作為邊界條件，聲學有限元系統(tǒng)方程形式為： 本次分析首先在MAXWELL進行電磁場分析，求解完成后，對電磁力進行FFT變換，在workbench平臺利用耦合功能，將其導入Mechanical進行簡諧振動分析，得到質(zhì)點振動速度，再將其導入ANSYS Acoustics聲學仿真模塊，求解聲壓波動方程，進行聲場分析，得到最后的噪聲計算結(jié)果，并根據(jù)GB/T 1094.10進行評定。 Figure.基于ANSYS Workbench的聲學仿真耦合流程 3 干式變壓器振動噪聲分析Figure.變壓器三維模型圖 Figure.噪聲分析耦合流程圖 3.1電磁場分析 將變壓器的電磁模型導入Maxwell,給定鐵芯、繞組的材料，設定好額定工況的激勵、邊界條件、求解參數(shù)，即可進行求解。設定好的繞組激勵如下圖所示： ① 設定鐵芯、繞組材料： Figure.材料設定 ② 施加激勵、求解計算： Figure.激勵加載&求解設置 ③ 后處理： Figure.后處理設置
Figure. 電磁力密度 3.2 結(jié)構(gòu)分析 在mechanical中進行分析前，首先根據(jù)提供的材料在Engineer Data中輸入材料數(shù)據(jù)，由于諧響應分析是線性分析類型，并且變壓器結(jié)構(gòu)在實際工作中也不允許超出屈服強度，因此此處以線彈性材料進行簡化輸入。網(wǎng)格劃分過程中，實體單元以四面體、六面體混合。根據(jù)實際工作，掃頻范圍設置為0~1000Hz。加載時，根據(jù)變壓器實際安裝位置，將下部的底座框架施加固定約束。具體操作如下： ① 網(wǎng)格劃分：針對模型不同部件，Mesh下插入Body Sizing，指定尺寸，生成網(wǎng)格。 Figure.網(wǎng)格劃分操作設置 Figure. 變壓器網(wǎng)格劃分 ② 邊界條件：根據(jù)實際工作情況，將底部進行全約束。在Harmonic Response處右鍵insert插入fixed support Figure.插入邊界條件 Figure. 變壓器邊界條件加載 ③ 分析設置：此處根據(jù)前述分析，將頻率區(qū)間設置為0~1000Hz Figure. 分析設置 ④ 導入電磁力：在Import Load處，鼠標右鍵Insert，選擇Surface Force Density，選擇需要導入電磁力的部件，Surface Force Density右鍵選擇import Load，即可導入。 Figure.導入電磁力設置 Figure. Import電磁力 ⑤ 后處理： Figure.后處理插入速度設置 Figure. 質(zhì)點振動速度云圖 3.3噪聲分析 噪聲分析利用ANSYS專業(yè)噪聲仿真模塊Acoustics。噪聲分析需要輸入聲音在介質(zhì)中的傳播速度及介質(zhì)密度等參數(shù)，此處介質(zhì)為空氣，在Engineer Data中輸入相應數(shù)據(jù)即可。噪聲分析由于主要分析聲音在介質(zhì)中傳播現(xiàn)象，因此需要設置空氣域。由于變壓器與空氣接觸部分幾何復雜，因此對空氣域采用四面體網(wǎng)格劃分方式。基于ANSYS Workbench耦合平臺，將上一步諧響應分析計算得到的質(zhì)點振動速度導入噪聲分析中，作為激勵源。通過計算可以得到不同頻率下的聲壓情況，由于輸入正弦激勵，頻率為50Hz，而一次交流過程中會有兩次信號達到峰值，因此振動分析的基礎頻率為100Hz。因此可以查看100,200,300等倍頻噪聲情況，此分析中僅截止到1000Hz。計算完成后，根據(jù)GB/T 1094.10變壓器聲級測定標準，后處理中提取相關輪廓線處A計權聲壓，并計算平均值，得到最終結(jié)果。 ① 模型處理：進行聲場分析，首先需要建立空氣域，在Design Modeler中利用Enslosure功能可以插入空氣域，同時指定空氣域大小即可。 Figure.插入空氣域 Figure. 空氣域的建立 ② 網(wǎng)格劃分：由于空氣域形狀復雜，此處以四面體方式進行網(wǎng)格劃分，此類特征的幾何模型適合采用Patch Independent算法進行網(wǎng)格劃分。此處Max Element Size指定為250mm。 Figure.網(wǎng)格劃分設置 Figure. 空氣域網(wǎng)格劃分
③ 邊界條件：右鍵單擊Import Load選擇Insert，插入Velocity，插入諧響應分析中計算得到的質(zhì)點振動速度作為聲場分析激勵。 Figure.Import Load設置 Figure. 導入所有頻段的質(zhì)點振動速度 ④ 分析設置：進行聲場分析前，需要選擇聲場區(qū)域。在Harmonic Acoustics處右鍵單擊，選擇Insert，選擇Physics Region,選擇我們繪制的聲場區(qū)域。 Figure.插入聲場區(qū)域設置
Figure.生成的聲場區(qū)域
⑤ 后處理：求解計算完成后，在Solution右鍵單擊，選擇Insert，選擇Acoustics，選擇我們關心的結(jié)果即可。 Figure.后處理設置 Figure. 100Hz聲壓分布（前后面） Figure. 100Hz聲壓分布（左右面） Figure. 不同頻率下聲壓變化曲線（前后面最大聲壓）
通過上述曲線，發(fā)現(xiàn)前后面聲壓最大發(fā)生在400Hz時。 Figure. 400Hz時前后面聲壓分布 Figure. 不同頻率下聲壓變化曲線（側(cè)面最大聲壓） 通過上述曲線，發(fā)現(xiàn)側(cè)面聲壓最大時為300Hz。 Figure. 300Hz時側(cè)面聲壓分布 根據(jù)實際試驗要求，提取輪廓線處的聲壓，并取平均值。 表格13 各點聲壓值 ①通過噪聲分析，發(fā)現(xiàn)變壓器在工作時，前后面的聲壓分布趨勢基本一致，側(cè)面的聲壓分布趨勢基本一致，最大值略有差異。結(jié)果說明 ②通過噪聲分析，發(fā)現(xiàn)該變壓前后面的最大A計權聲壓為58dB，側(cè)面最大A計權聲壓為50dB。 ③通過噪聲分析后處理，300Hz平均聲壓為50.4dB，400Hz平均聲壓為得到平均為50.8dB。 4 總結(jié) 本文通過基于ANSYS Workbench平臺的干式變壓器振動噪聲仿真，實現(xiàn)了在產(chǎn)品設計階段對其噪聲值進行預估的完整流程，可以幫助企業(yè)在探究變壓器噪聲的機理上，對產(chǎn)品及時做出改進，響應市場，提高競爭力。 |
