變壓器鐵心電磁振動仿真及影響因素研究

聲學仿真初學者

2023年6月30日 16:07

瀏覽：4967 評論：2 收藏：4

進入 21 世紀以來，伴隨著中國社會經濟的發展，電力行業蓬勃發展，電網規模持續擴大，用戶對供電質量及穩定性的要求也越來越高。電力變壓器不僅在輸配電方面起到重要作用，也是電力系統中消耗能源的主要設備之一。目前城市建設的供電需求很大，電力部門一般采用具有更高負載能力的變壓器，變壓器的負載等級越高，鐵心的振動噪聲越大，會在不同程度上干擾附近居民的生活和學習。因此，變壓器運行過程中的振動噪聲問題一直是變壓器生產制造企業面臨的嚴峻挑戰。

自 20 世紀 70 年代以來，中外針對電力變壓器鐵心的振動特性開展了大量實驗與仿真研究，且研究規模不斷擴大^［1］。趙莉華等^［2］通過實驗研究分析了變壓器鐵心的振動，得到了不同工況下鐵心的振動頻譜特性。韓芳旭等^［3］基于磁致伸縮力-熱應力比擬的數值計算方法建立電磁場數值模型，求解鐵心每個節點不同時刻的磁密值，加載試驗測得的硅鋼片磁致伸縮特性曲線，仿真得到鐵心每個時間步各個節點的磁致伸縮力，導入到結構場計算模型中求得鐵心本體的振動位移。在鐵心振動模型的研究方面，朱葉葉等^［4］、張黎等^［5］建立了鐵心材料磁致伸縮的本質模型，利用彈性力學原理描述硅鋼片材料的本構關系，將不同磁感應強度下的磁致伸縮應變轉化為應力，采用弱耦合的形式對鐵心應力場進行仿真分析。祝麗花^［6］采用方圈法測試硅鋼片磁化特性以及材料的磁致伸縮數據，建立電磁-結構耦合模型，仿真獲得了鐵心磁場與振動位移。王佳音^［7］詳細測量了多種取樣方向硅鋼片的磁化曲線與磁致伸縮曲線，獲得了比較詳細的材料各向異性數據，便于模擬各種情形下的仿真條件。張哲^［8］建立了考慮材料磁致伸縮特性的磁-機械耦合模型，相比于硅鋼片電機，非晶合金電機鐵心振動量更大，且磁致伸縮受應力影響程度更加明顯。張鵬寧等^［9］從直流偏磁機理和振動噪聲基本原理著手，將電磁場、結構力場和聲場進行耦合計算完成直流偏磁下鐵心振動和噪聲問題的研究，分析了偏磁狀態下鐵心本體的振動情況，得到了一般性結論。祝麗花等^［10］利用能量變分原理以強耦合方式描述鐵心磁場與應力場之間的關系，建立了電力變壓器鐵心磁致伸縮三維仿真模型，計算了鐵心在空載運行時的磁感應強度分布和振動分布。魏亞軍等^［11］、莫娟等^［12］通過建立電磁場與應力場的偏微分方程，研究了鐵心和繞組在不同負載條件下的磁通分布和應力應變。劉宏亮^［13］提出在變壓器鐵心的接縫處，硅鋼片的磁感應強度分布較為復雜，在垂直硅鋼片材料軋制方向上，磁致伸縮增大了數十倍，比較合理地解釋了鐵心在接縫處振動較大的問題。張鵬寧等^［14］對比不同磁致伸縮數學模型，利用 COMSOL多場仿真軟件建立了一臺高壓電抗器鐵心本體結構的振動計算模型。Tong 等^［15］針對飽和電抗器進行減振方法研究，認為鐵心中的應力始終存在并影響硅鋼片材料屬性，硅鋼片在直流激勵條件下表現出不同的磁化特性和磁致伸縮特性，在不同拉伸應力下的磁化強度和磁致伸縮曲線也不相同。Tanzer等^［16］研究了兩種材料的磁致伸縮特性的差異，將兩種材料制作成鐵心樣機進行振動測試，比較了兩者之間的振動趨勢。Liu 等^［17］研究了一種全耦合預測變壓器鐵心變形的磁-機械耦合方法，考慮用簡化的多尺度本構模型來描述材料磁化特性和磁致伸縮各向異性，將材料磁致伸縮作為結構場的系統輸入以獲得鐵心變形和位移場。

綜合來看，對變壓器鐵心振動特性進行數值分析研究時，不得不考慮硅鋼片材料的非線性特征及其復雜性。外加應力不僅會改變材料的磁化特性，影響電磁場分析的精度，同時也會改變材料的磁致伸縮特性，進而影響結構場分析的準確度。因此在進行數值分析時需要充分了解材料特性的各部分影響因素，以期得到科學準確的仿真結果。在磁-結構耦合方面，大多數學者基于廣義胡克定律對鐵心的磁致伸縮力進行直接等效，沒有考慮到外界裝配約束條件對鐵心振動的影響，導致鐵心各方向振動求解結果與實際情況存在很大差異，具有一定的局限性。為了研究裝配工藝參數條件的差異導致鐵心振動發生較大變化的問題，分析鐵心電磁振動的影響因素，通過改變材料的磁致伸縮增量以及疊片方向上的彈性模量，利用應變能原理對磁致伸縮力進行等效，對鐵心振動量進行參數研究，為企業的產品設計與生產工藝改進提供參考依據。

1 電磁-結構場耦合模型

1. 1 電場和磁場耦合模型

鐵心在空載運行時的標準勵磁電壓為 400 V，通過每匝線圈的電壓幅值計算得到所需線圈匝數，對線圈進行參數設置。在電路中對 3 個線圈施加三相交流電，相位差依次相差 2π/3，線圈匝數設定為24，并設置線圈電阻，電路中激勵電壓的方程為

在鐵心的電磁場分析中，硅鋼片材料的 B-H( B為磁感應強度，H 為磁場強度) 磁化特性是非常重要的輸入參數，直接影響電磁場分析結果的準確度。硅鋼片材料 B-H 特性曲線是非線性的，工程上直接采用硅鋼卷出廠時的材料特性，沒有考慮剪切導致的材料特性變化。由于剪切加工后的材料 B-H特性與原料的差別較大，磁導率下降，因此將實驗測試獲得的鐵心磁化性能參數引入鐵心電磁模型。利用樣條插值擬合曲線的方法設置硅鋼片材料的磁化特性，如圖 1 所示。

通過電路模塊在線圈中施加電壓，可以計算出線圈中總電流密度 J ，將其作為激勵施加在磁場中，電場和磁場之間的關系通常使用麥克斯韋方程組表示，即

1. 2 磁場和結構場耦合模型

鐵心在外加電壓的激勵下會產生磁場，即在鐵心的閉合回路中產生了磁感應強度。硅鋼片所受磁致伸縮效應與鐵心電磁場分布緊密相關，當鐵心閉合回路中磁感應強度發生變化時，鐵心等效磁致伸縮力會發生變化，鐵心應力場中的激勵源發生變化，進而導致鐵心表面振動發生改變。求解磁場和振動的方程組為

由于鐵心磁致伸縮效應和磁感應強度直接相關，等效磁致伸縮力導致的鐵心變形對鐵心的電磁場分布影響很小，因此可以采用弱耦合的方式來描述鐵心中磁場和應力場的關系。以應變為耦合變量，建立鐵心磁-結構耦合模型為

建立鐵心的振動方程前需要明確鐵心所受到的外界激勵狀態，利用彈性力學原理可以將鐵心材料應變轉化為激勵力。在確定鐵心各個方向上的彈性模量后，利用應變能原理將硅鋼片磁致伸縮效應等效為磁致伸縮力，將其作為鐵心振動的輸入激勵信號，鐵心的振動方程為

2 鐵心電磁振動仿真分析

仿真模擬對象為某型號干式變壓器鐵心，鐵心結構形式為三相三柱式，每一相的初級繞組線圈和次級繞組線圈均安裝在同一根心柱上。首先建立鐵心的三維幾何模型，如果在仿真中考慮到每一片硅鋼片則導致計算量過大，因此忽略每一級的疊片形式，保留鐵心各級的尺寸，將鐵心的幾何模型按照實際情況劃分為 5 個部分，包括 2 個鐵心鐵軛、1 個鐵心心柱以及 2 個鐵心邊柱。仿真中只添加了初級線圈繞組，忽略次級線圈及其他夾件、底腳、絕緣材料等細節性的因素。在 COMSOL 軟件中采用四面體網格的方式對鐵心幾何模型進行劃分，為了保證較為精確的計算結果，將最大網格的邊長設置為 50 mm。鐵心的有限元網格如圖 2 所示。

2. 1 電磁場分析

對電磁場進行瞬態求解，求解器步長設置為0. 000 5 s，總求解時間為 0. 05 s，計算得到鐵心的磁感應強度分布。在 t = 0. 01 s 時鐵心表面磁感應強度與線圈電流密度分布云圖如圖 3 所示。

從磁感應強度分布圖中可以看出，鐵心磁感應強度在 0. 01 s 時間間隔內不斷變化，磁感應強度分布呈周期性的分布規律，鐵心的閉合回路中各部位磁感應強度幅值不同，變化規律也不相同，如圖 4 所示。其中，在鐵心窗口內側的磁感應強度幅值較大，在鐵心窗口外側的磁感應強度幅值較小，鐵心心柱以及接縫處磁感應強度較大。

在時間為 0. 01 s 時，鐵心在各個方向上的磁感應強度分布如圖 5 所示。其中，X 方向表示鐵心的窗寬方向，Y 方向表示在鐵心的窗高方向，Z 方向表示鐵心的疊片方向。磁感應強度主要分布在 X 方向與 Y 方向上，Z 方向存在較小的磁感應強度。在其他區域，鐵心電磁場存在以下分布規律: X 方向磁感應強度主要分布在鐵心的上軛和下軛處，Y 方向磁感應強度主要分布在鐵心的 3 個心柱上; Z 方向磁感應強度非常弱，幅值接近于 0。

2. 2 磁-結構耦合

鐵心的實際運行振動情況較為復雜，鐵心 Z 方向即疊片方向振動明顯，而直接利用彈性力學原理對磁致伸縮力進行等效會導致 Z 方向的體積力很小，因此采用應變能原理對鐵心磁-結構耦合關系進行修正。

由于鐵心表面的振動為周期性振動，在鐵心運行過程中外界夾緊應力幾乎不變，機械應力對鐵心單元體的做功很小。因此忽略外界機械應力的影響，僅考慮等效磁致伸縮力做功，認為鐵心磁致伸縮效應的應變能全部轉化為鐵心的振動能量，基于彈性體的應變能公式建立鐵心的能量公式。彈性體的應變能的定義式為

在企業對鐵心產品的各方向振動情況進行統計，發現鐵心振動速度分布具有一般性結論，即鐵心心柱和鐵軛上的振動速度在各個方向上的分布存在比值關系。在鐵心 X 方向上的振動速度近似為 Y 方向的 1. 2 倍，對于鐵心心柱而言，Z 方向上的振動速度近似為 Y 方向的 3 倍，對于鐵軛而言，Z 方向上的振動速度近似為 Y 方向的 8 倍。

根據鐵心振動速度在各個方向上的差異，將鐵心心柱及鐵軛處不同方向的振動速度比值關系代入式( 9) 中，推導得到鐵心心柱單元體及鐵軛單元體的能量公式為

將體積力作為結構場輸入變量引入到鐵心振動模型中，鐵心在受到磁致伸縮作用時的體積力分布如圖 6 所示。從圖 6 中可以看到，鐵心在各個方向上均存在等效體積力分布，且與鐵心的磁感應強度分布存在一定的關系，在 Z 方向上鐵心的體積力最大。與直接采用彈性力學原理等效的方法相比，在總振動能量不變的情況下，該方法使得鐵心的一部分受力從其他兩個方向轉移到了 Z 方向上。

2. 3 等效磁致伸縮力作用仿真結果

為了得到磁致伸縮力對鐵心振動的影響程度，利用應變能原理將等效得到的磁致伸縮力作為輸入變量引入到鐵心結構場中，仿真獲得鐵心的瞬態振動分布，圖 7 表示鐵心 3 個時刻的振動速度分布。

從圖 7 中可以看出，鐵心上軛及心柱上半部分的振動速度較大，符合立式鐵心的振動分布規律，仿真最大值不超過8×10^－4m / s^－1，與大量研究中測試的鐵心振動分布趨勢相一致。由于鐵心結構的外在影響因素很多，即便同一型號的多臺鐵心振動值也存在較大差異，除異常振動測點外，統計了某型號干式變壓器鐵心表面振動最大值為 ( 6. 2～9. 7)×10^－4m / s，仿真結果與實驗可以相互印證。在仿真過程中，通過考慮鐵心各方向振動分布的差異性，改變磁場與結構場的耦合方式，能夠獲得較為準確的仿真計算結果。

3 鐵心振動影響因素

大量實驗研究表明，加載勵磁電壓后鐵心振動的大小主要取決于硅鋼片的材料特性，在硅鋼片材料方面主要表現在磁致伸縮增量與等效彈性模量。材料在外加磁場狀態下會發生磁致伸縮應變，且這種應變隨著外界應力載荷變化而變化，直接關系到鐵心的振動大小。同時鐵心在疊片方向上的等效彈性模量也不容忽視，不同的外界載荷條件會使得硅鋼片層間結合面特性發生變化，其疊片方向上的等效彈性模量發生變化，進而間接影響鐵心的振動大小。鐵心的裝配工藝主要改變了鐵心層間夾緊力大小，不同的裝配載荷會改變硅鋼片鐵心的磁致伸縮以及彈性模量，下面分別對這兩種影響因素進行分析。

3. 1 硅鋼片磁致伸縮增量

對鐵心進行夾件裝配夾緊時，硅鋼片主要受到壓應力的影響。以往的研究表明［7］，硅鋼片磁致伸縮幅值隨壓應力變化很大，拉應力對其影響程度很小，壓應力變大，磁致伸縮增長迅速，不同應力條件下的硅鋼片磁致伸縮增量如圖 8 所示。

隨著預緊壓應力增大，磁致伸縮增量比值增大。利用磁致伸縮增量比值修正鐵心材料的磁致伸縮單值曲線，并導入到鐵心結構場中進行分析，計算得到不同磁致伸縮增量條件下鐵心的振動量。

圖 9 所示為鐵心不同磁致伸縮增量條件下鐵心各方向上的振動速度均值。

從圖 9 中可以看出，隨著壓應力增大，鐵心材料的磁致伸縮增量比值逐漸增大，鐵心在各方向上的振動均明顯增大且增速明顯。

3. 2 疊片方向彈性模量

由于鐵心為疊片結構，疊片方向彈性模量遠小于硅鋼片材料屬性。在鐵心進行夾件裝配夾緊時，硅鋼片層間的結合面參數會發生變化，進而整個鐵心疊片方向的等效彈性模量發生變化，設置不同彈性模量參數的鐵心具有不同的動力學性能。仿真得到鐵心在不同該參數下各方向上的振動速度均值，如圖 10 所示。

從圖 10 中可以看出，隨著疊片方向彈性模量從50 GPa 增加至 200 GPa，鐵心各方向上的振動速度逐漸降低。通過對鐵心振動的兩種影響因素進行仿真分析，可以定性地得到以下結論: 在鐵心進行夾緊裝配時，夾緊壓強低于 0. 1 MPa，隨著夾緊壓強增大，鐵心剛度提高，疊片方向彈性模量為主要影響因素，磁致伸縮增量比值相對較小，導致鐵心振動降低; 夾緊壓強高于 0. 15 MPa，隨著夾緊壓強增大，鐵心磁致伸縮增量為主要影響因素，磁致伸縮效應加劇，導致鐵心振動增大，因此鐵心夾緊壓強理論上存在最優的夾緊區間。