李凱琪， 郭 健， 劉 軍， 張安紅， 俞少佳

(1. 南京航空航天大學(xué)， 江蘇 南京 211106； 2. 錢江電氣集團股份有限公司， 浙江 杭州 311243)

一種高效的電力變壓器繞組軸向固有振動頻率計算模型

李凱琪1， 郭 健1， 劉 軍2， 張安紅2， 俞少佳2

(1. 南京航空航天大學(xué)， 江蘇 南京 211106； 2. 錢江電氣集團股份有限公司， 浙江 杭州 311243)

在電力變壓器設(shè)計過程中，繞組的固有振動頻率是一個非常重要的性能參數(shù)。本文基于彈性動力學(xué)理論，根據(jù)繞組周向墊塊彈性支撐等效原則，提出了一種高效的計算電力變壓器繞組軸向固有振動頻率的二維有限元模型，其固有振動頻率的計算具有與三維模型一致的結(jié)果，但其計算效率得到大大的提高。同時研究表明，增加墊塊的寬度和數(shù)量，以及增大預(yù)緊力都可以提高繞組軸向的固有振動頻率，是避免繞組共振的一種有效應(yīng)對方法。最后通過對固有振動頻率計算的實例校核，驗證了本文所提二維模型方法的有效性。

電力變壓器； 繞組； 高效有限元模型； 軸向固有振動頻率

1 引言

在電力變壓器設(shè)計過程中，繞組的固有振動頻率是一個非常重要的性能參數(shù)[1]，當(dāng)固有振動頻率與電動力的頻率相接近時，將產(chǎn)生諧振現(xiàn)象[2]。這種振動可能會導(dǎo)致比預(yù)緊力大得多的動態(tài)力，從而使油道墊塊與線餅之間分開，進(jìn)而引起繞組的松散、扭轉(zhuǎn)，甚至倒塌，最終使得變壓器繞組失穩(wěn)而無法正常運行[3,4]。對于大型電力變壓器，由于繞組輻向固有振動頻率要遠(yuǎn)遠(yuǎn)大于其運行過程中承受的電動力頻率，故而輻向固有振動頻率在變壓器設(shè)計過程中一般不考慮。而繞組軸向由于采用導(dǎo)線-墊塊-導(dǎo)線的結(jié)構(gòu)形式，其軸向固有振動頻率值大大降低，結(jié)構(gòu)設(shè)計稍不合理，就會導(dǎo)致線圈發(fā)生共振，這也成為變壓器軸向失穩(wěn)或運行噪音過大的一部分原因。因此，準(zhǔn)確計算變壓器繞組的固有振動頻率，并通過改善結(jié)構(gòu)形式使其避開軸向動態(tài)短路力的頻率對于變壓器設(shè)計和安全可靠運行具有重要的意義。

目前，變壓器軸向固有振動頻率的計算主要是基于“質(zhì)量-彈簧-阻尼”的集中參數(shù)法和有限元法，由于有限元模型可以很好地兼顧繞組結(jié)構(gòu)形式、支撐件結(jié)構(gòu)力學(xué)特性等因素，相對于傳統(tǒng)的集中參數(shù)法可以得到更加接近實際的計算結(jié)果。文獻(xiàn)[5]建立了電力變壓器固有振動頻率計算的三維有限元模型，在軸向固有振動頻率計算上取得了很好的效果；文獻(xiàn)[6]在此基礎(chǔ)上研究了軸向預(yù)緊力對固有振動頻率的影響規(guī)律，并提出了改善變壓器軸向穩(wěn)定性的措施。綜上，三維有限元模型在變壓器軸向固有振動頻率計算上具有較好的結(jié)果，但由于這種方法建模過程復(fù)雜，計算代價高，從而限制了該方法在工程實際中的普遍應(yīng)用。

本文基于彈性動力學(xué)理論[7]，根據(jù)繞組周向墊塊彈性支撐等效原則，提出了一種高效的計算電力變壓器繞組軸向固有振動頻率的二維有限元模型，其固有振動頻率的計算結(jié)果具有與三維模型一致的效果，同時相對于三維模型，其計算效率亦得到大大的提高。最后，通過對固有振動頻率計算的實例校核，驗證了本文所提二維模型方法的有效性。

2 二維有限元模型

2.1 模型的提出

大型電力變壓器的三維軸對稱模型如圖1所示。其中繞組的線餅是由絕緣墊塊隔開的線段組成的，并被壓緊在夾件之間；繞組中的墊塊在相鄰兩餅銅線中起到一個支撐的作用，是變壓器機械系統(tǒng)可靠性的一個重要的影響因素，在固有振動頻率計算模型中必須予以考慮。為了能夠兼顧到墊塊沿繞組周向布置的特點，計算模型應(yīng)屬于三維問題，但是這將會使得固有振動頻率計算的代價大大提高。

圖1 三維軸對稱模型Fig.1 3D axisymmetric model of winding

絕緣墊塊受到導(dǎo)線軸向電動力壓縮時，其受到的應(yīng)力滿足：

(1)

式中，δ為墊塊受到的應(yīng)力；F為墊塊受到的軸向壓力；∑A為兩線餅間的墊塊與線餅接觸的實際面積之和。

根據(jù)胡克定理，墊塊的應(yīng)力-應(yīng)變滿足：

(2)

式中，ε為墊塊的應(yīng)變；E為墊塊材料的彈性模量。

絕緣墊塊由于受壓導(dǎo)致線餅在軸向發(fā)生ΔZ的位移，位移與壓力之間滿足式(3)：

(3)

式中，k為墊塊的等效彈性系數(shù)，反映了墊塊對線餅的支撐作用。

由于線餅位移與墊塊的應(yīng)變滿足ΔZ∝ε，因此

(4)

聯(lián)立式(1)、式(2)和式(4)，則可以得到墊塊的等效彈性系數(shù)與墊塊材料的彈性模量，及墊塊與線餅接觸面積之間的關(guān)系：

(5)

由式(5)可以看出，只要E∑A為常數(shù)，則可以保證墊塊對線餅的支撐效果不變。由此就可以將固有振動頻率計算的三維問題簡化為二維模型進(jìn)行處理。

圖2(a)為滿足軸對稱的三維計算模型，在該計算模型中絕緣墊塊與線餅之間滿足全接觸，進(jìn)一步可以得到簡化的二維計算模型，如圖2(b)所示。為了保證墊塊對線餅的支撐效果不變，根據(jù)式(5)，墊塊材料的等效彈性模量應(yīng)滿足：

(6)

式中，Ee為等效彈性模量；r1、r2分別為線餅的內(nèi)、外半徑。

圖2 繞組軸對稱模型Fig.2 Axisymmetric models of winding

為了考慮軸對稱處理后墊塊質(zhì)量對軸向固有振動頻率的影響，還需要對墊塊的密度做等效處理，等效密度滿足：

(7)

式中，ρe為等效密度；ρ為墊塊材料的實際密度。

2.2 控制方程

假定繞組的線餅是由絕緣墊塊隔開的線段組成的，并被壓緊在夾件之間。在研究繞組中的動態(tài)過程中，可以把它視為由有彈性聯(lián)系的實體(線段)所組成的機械系統(tǒng)[8]。

基于圖2(b)提出的有限元模型，并根據(jù)彈性動力學(xué)理論，各單元運動方程為：

(8)

式中，Z為各線餅絕對位移；M為系統(tǒng)的質(zhì)量矩陣；C為系統(tǒng)的阻尼系數(shù)矩陣；K為系統(tǒng)剛度系數(shù)矩陣；F(t)為墊塊受到的軸向壓力矩陣；g為重力加速度。

忽略阻尼的影響，且沒有外力作用，可得繞組自由振動方程為：

(9)

式(9)為線性齊次常微分方程組，其解可表示為：

(10)

將式(10)代入式(9)，得：

(11)

該方程具有非零解的條件為：

(12)

解式(12)可得n個正實根ωi，則繞組軸向振動的固有振動頻率為：

(13)

3 計算結(jié)果

3.1 二維模型計算結(jié)果

本文以文獻(xiàn)[9]提出的模型為例，進(jìn)行二維計算模型有效性的驗證，模型對應(yīng)的參數(shù)如表1所示。

表1 模型參數(shù)

根據(jù)本文所提模型分析得到的固有振動頻率及對應(yīng)的振型如表2所示。

表2 固有振動頻率對應(yīng)振型

3.2 二維與三維計算比較

本文還將三維模型與二維模型的計算進(jìn)行了比較。圖3(a)為三維模型的網(wǎng)格劃分圖，圖3(b)為三維模型得到的一階軸向固有振動頻率振型。

圖3 固有振動頻率對應(yīng)振型Fig.3 Vibration models according to natural vibration frequency

表3列出了兩種模型的計算結(jié)果與計算代價的比較?？梢钥闯鰞煞N方法計算結(jié)果非常接近，但是在計算時間上，二維模型的計算效率要明顯優(yōu)于三維模型。

表3 固有振動頻率結(jié)果比較

3.3 實驗驗證

對變壓器模型的固有振動頻率進(jìn)行試驗測試，表4列出了實驗結(jié)果。可以看出，計算結(jié)果和實驗所得結(jié)果基本吻合，證明了本文所提方法的有效性。

表4 二維模型結(jié)果和實驗結(jié)果的比較

4 影響因素研究

由固有振動頻率求解控制方程式(12)可知，固有振動頻率的大小與機械系統(tǒng)等效的彈性系數(shù)及質(zhì)量有關(guān)，系統(tǒng)質(zhì)量的大小一般取決于變壓器的電磁設(shè)計和經(jīng)濟成本，調(diào)整空間不大；而系統(tǒng)的彈性系數(shù)則取決于墊塊對線餅的支撐，一方面，可以通過調(diào)整墊塊的個數(shù)及尺寸來實現(xiàn)，另一方面，可以通過改變預(yù)緊力來實現(xiàn)對軸向固有頻率的調(diào)整，是改善系統(tǒng)固有振動頻率，避免共振的有效方法[10]。

本文以表1模型參數(shù)為例研究影響變壓器繞組軸向固有振動頻率的因素。表5列出了不同墊塊的寬度對固有振動頻率的影響；圖4為對應(yīng)的一階固有振動頻率關(guān)系曲線圖?？梢钥闯?，隨著墊塊寬度的增加，墊塊對線餅有效的支撐面增大，使得系統(tǒng)等效的彈性系數(shù)k增大，因此系統(tǒng)的固有振動頻率亦成比例增大。

表5 變壓器繞組不同寬度墊塊時的固有振動頻率

圖4 墊塊寬度-頻率曲線Fig.4 Curve of pressboard’s width and frequency

表6列出了不同墊塊的數(shù)量對固有振動頻率的影響；圖5為對應(yīng)的一階固有振動頻率關(guān)系曲線圖?？梢钥闯?，增加墊塊的數(shù)量與增加墊塊的寬度對固有振動頻率具有相同的影響效果。

表6 變壓器繞組不同個數(shù)墊塊時的固有振動頻率

圖5 墊塊個數(shù)-頻率曲線Fig.5 Curve of pressboard’s number and frequency

變壓器繞組的絕緣主要是墊塊及匝絕緣，其材料為絕緣紙板[11]。對于絕緣紙板，其靜剛度為非線性，絕緣紙板在一定范圍內(nèi)(δ：0～103MPa)應(yīng)力-應(yīng)變關(guān)系可表示為：

(14)

式中，a為線性系數(shù)；b為硬化系數(shù)，可通過加載-變形試驗測試得到[12,13]。a=1.05×102MPa，b=1.75×103MPa。

則彈性模量為：

(15)

假設(shè)施加在線圈上的預(yù)緊力為T，則

(16)

因此，可以通過改變等效的墊塊的彈性模量來模擬預(yù)緊力的大小。

對應(yīng)的等效剛度系數(shù)為：

(17)

式中，h為絕緣墊塊的高度。

表7列出了變壓器繞組在不同預(yù)緊力下的固有振動頻率值。圖6為對應(yīng)的一階固有振動頻率關(guān)系曲線圖。

表7 變壓器繞組不同預(yù)緊力下的固有振動頻率

圖6 預(yù)緊力-頻率曲線Fig.6 Curve of pre-compression and frequency

可以看出，隨著軸向預(yù)緊力的增大，繞組的各階軸向振動固有頻率都在提高，這是因為繞組的固有振動頻率在很大程度上取決于裝配后的壓緊程度，絕緣墊塊的彈性系數(shù)與此有關(guān)。因此，可以通過改變等效的墊塊彈性模量來模擬預(yù)緊力的大小。

表5～表7和圖4～圖6數(shù)據(jù)表明了墊塊的個數(shù)和尺寸以及預(yù)緊力對變壓器繞組的固有振動頻率的影響。在設(shè)計變壓器時，應(yīng)當(dāng)考慮這幾點，使得變壓器繞組軸向的固有振動頻率避開軸向動態(tài)短路力的頻率(50Hz或100Hz)。

5 結(jié)論

本文基于彈性動力學(xué)理論，根據(jù)繞組周向墊塊彈性支撐等效原則，提出了一種高效的計算電力變壓器繞組軸向固有振動頻率的二維有限元模型，其固有振動頻率的計算具有與三維模型一致的結(jié)果，但其計算效率得到大大的提高。同時研究表明，增加墊塊的寬度和數(shù)量以及增大預(yù)緊力都可以提高繞組軸向的固有振動頻率，在設(shè)計時應(yīng)當(dāng)給予考慮，避免繞組產(chǎn)生共振。最后通過對固有振動頻率計算的實例校核，驗證了本文所提二維模型方法的有效性。

[ 1] 王錄亮, 劉文里, 于會鳳 (Wang Luliang, Liu Wenli, Yu Huifeng). 電力變壓器繞組軸向短路動態(tài)力的計算(Calculation of winding axial short-circuit dynamic force of power transformer) [J]. 變壓器 (Transformer), 2013,50(1): 14-17.

[ 2] Inge Johansen, Harald Riege. Numerical evaluation of natural frequencies in power transformer windings, power apparatus and systems [J]. Power Apparatus and Systems, Part III - Transactions of the American Institute of Electrical Engineers, 1962, 81 (3): 782-788.

[ 3] 李文海 (Li Wenhai). 變壓器短路強度研究成果簡介(Brief introduction of achievements in transformer short-circuit strength research) [J]. 變壓器 (Transformer), 1997, 34(10):18-22.

[ 4] S K Sahoo, S Gopalakrishna. Winding deformation analysis in power transformers using finite element method [A]. 2014 IEEE Students’ Technology Symposium (TechSym) [C]. 2014. 341-346.

[ 5] 余小輝, 李巖, 井永騰, 等 (Yu Xiaohui, Li Yan, Jing Yongteng, et al.). 變壓器繞組模型的固有振動頻率計算與分析 (Calculation and analysis of natural frequency of winding model of transformer) [J]. 變壓器 (Transformer), 2010, 47(7): 6-8.

[ 6] 邵宇鷹, 饒柱石, 謝坡岸, 等 (Shao Yuying, Rao Zhushi, Xie Poan, et al.). 預(yù)緊力對變壓器繞組固有振動頻率的影響 (The influence of pre-compression changing on the natural frequency of electric transformers windings) [J]. 噪聲與振動控制 (Noise and Vibration Control), 2006, 26(6): 51-53.

[ 7] Yan Li, Wei Zhou, Yongteng Jing, et al. Axial vibration analysis of power transformer active part under short-circuit [A]. 2011 International Conference on Electrical Machines and Systems (ICEMS) [C]. 2011.1-4.

[ 8] 汲勝昌, 王俊德, 李彥明 (Ji Shengchang, Wang Junde, Li Yanming). 穩(wěn)態(tài)條件下變壓器繞組軸向振動特性研究(Research on axial vibration characteristics of transformer windings under steady condition) [J]. 電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy), 2006, 25(1): 35-38, 48.

[ 9] Hadiya Pritesh. Coupled dynamic 3D simulation of short circuit forces [D]. Blekinge Institute of Technology, 2012.

[10] Zhang Bo, Li Yan, Fang Xiaowen. Engineering-oriented research of axial vibration in power transformer windings [A]. 2014 International Conference on Information Science, Electronics and Electrical Engineering (ISEEE) [C]. 2014. 341-345.

[11] 李洪奎,李巖 (Li Hongkui, Li Yan). 不同預(yù)緊力下變壓器繞組軸向振動模態(tài)分析(Axial vibration modal analysis of transformer windings under different levels of pre-compression) [J]. 電機與控制學(xué)報 (Electric Machines and Control), 2010, 14(8):98-101, 106.

[12] M R Patel. Dynamic response of power transformers under axial short-circuit force, Part Ⅰ- winding and clamp as individual components [J]. IEEE Transactions on Power Apparatus and Systems, 1973, PAS-92(5): 1558-1566.

[13] M R Patel. Dynamic response of power transformers under axial short-circuit force, Part II - windings and clamps as a combined system [J]. IEEE Transactions on Power Apparatus and Systems, 1973, PAS-92(5): 1567-1576.

Efficient model to calculate axial natural vibration frequency of power transformer winding

LI Kai-qi1, GUO Jian1, LIU Jun2, ZHANG An-hong2, YU Shao-jia2

(1. Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China; 2. Hangzhou Qiantang River Electric Group Co. Ltd., Hangzhou 311243, China)

In the design of transformer winding, natural vibration frequency is an important parameter. Based on the elastic dynamics theory, and according to the elastic support equivalent principle of radial pressboards, a 2D model to calculate axial vibration natural frequency of power transformer winding is presented. The 3D model to calculate natural vibration frequency can be simplified as a 2D one as the support of pressboard on the winding is same. It is tested that results of the 2D model are consistent with those of 3D model, but the former can achieve much higher calculation efficiency. It shows that increasing the width and number of pressboard can improve axial natural frequency through formula analysis and simulation, and also the relations between the changes of axial pre-compression and axial natural vibration frequencies on the windings are investigated. So these factors should be taken into consideration in the design of transformer winding, which can be an effective way to avoid the resonance of winding. Finally, the proposed 2D model’s effectiveness is proved when compared with tested ones.

power transformer; winding; efficient finite element model; axial natural vibration frequency

2015-03-08

國家自然科學(xué)基金資助項目(51307082)

李凱琪(1991-)， 女， 江蘇籍， 碩士研究生， 研究方向為電力變壓器短路強度; 郭 健(1980-)， 男， 山西籍， 副教授， 博士, 研究方向為特種電機的設(shè)計。

TM401+.1

A

1003-3076(2016)02-0048-06