應(yīng)用有限元法分析變壓器繞組固有振動特性

王豐華 楊毅 何苗忠 杜勝磊 劉石 楊賢

摘 要：為進(jìn)一步理解和掌握變壓器繞組的固有振動特性，建立變壓器繞組的三維有限元模型，給出變壓器繞組正常與預(yù)緊力改變時的固有頻率及其對應(yīng)的振型。同時以某10 kV實體變壓器繞組為試驗對象，對其模態(tài)特性進(jìn)行了測試分析。計算得到的變壓器繞組固有振動特性與測試結(jié)果吻合良好，說明所建立的變壓器繞組有限元模型的正確性。此外，變壓器繞組正常狀態(tài)時，前四階固有頻率對應(yīng)的振型呈現(xiàn)不同的振動形態(tài)，且對稱性較好。繞組預(yù)緊力改變時，變壓器繞組各階固有頻率均隨預(yù)緊力增大而增大，但增幅各異，高階固有頻率對預(yù)緊力的變化更為敏感。

關(guān)鍵詞：變壓器；繞組；振動；固有頻率；振型

中圖分類號：TM 401.1

文獻(xiàn)標(biāo)志碼：A

文章編號：1007-449X（2018）04-0051-07

Abstract：To further understand and master the inherent vibration features of transformer winding， the 3-D finite element model of transformer winding is built to calculate the natural frequency and the corresponding vibration pattern. Meanwhile， the model experiment of some real 10kV power transformer was made on the transformer winding. The results of simulation and experiment match well， which verifies the correctness of the simulation model. It is shown that the vibration patterns corresponding vibration patterns were in different shapes and good symmetry when the winding was under normal condition. When the precompression force of transformer winding was changed， the first four-order natural frequencies were increased with the increasing of precompression in different amplitudes.The high order natural frequency was more sensitive to the variations of precompression.

Keywords：power transformer； winding； vibration； inherent frequency； vibration pattern

0 引 言

統(tǒng)計結(jié)果表明[1-2]，變壓器大多數(shù)故障源于出口短路引發(fā)的繞組變形或松動。突發(fā)短路時，巨大的短路電動力可能會使變壓器繞組出現(xiàn)松動或變形等故障隱患，消弱變壓器繞組的機械穩(wěn)定性，此時，即使遭受較小的短路沖擊也有可能引發(fā)變壓器繞組坍塌等大的事故。因此，有必要研究及時有效的變壓器繞組狀態(tài)檢測方法，提高變壓器及電力系統(tǒng)的安全可靠運行。

近年來，從變壓器繞組機械動力學(xué)特性出發(fā)的振動檢測法日益成為國內(nèi)外的研究熱點，鑒于對變壓器繞組振動特性進(jìn)行建模研究可為獲取繞組振動規(guī)律提供重要的理論依據(jù)，研究者從上個世紀(jì)70年代初就開始了對變壓器繞組的軸向振動特性進(jìn)行研究。主要分為兩大類：基于彈簧-質(zhì)量系統(tǒng)的變壓器軸向振動集中參數(shù)模型和有限元分析模型。第一類模型最早由Swihart[3]等提出，即將線餅等效為剛體，將絕緣墊塊等效為彈性元件，然后通過該模型得到了繞組在電磁力作用下的位移、加速度等特征。文獻(xiàn)[4-5]分析了短路電動力作用下變壓器繞組的軸向振動特性，得到了不同預(yù)緊力下線餅軸向位移隨電流和繞組分接區(qū)域的變化規(guī)律。文獻(xiàn)[6]以電纜變壓器為研究對象，計算分析了瞬態(tài)短路力激勵下的繞組振動規(guī)律，討論了預(yù)緊壓強、撐條根數(shù)等因素對繞組軸向振動的影響。文獻(xiàn)[7]定性分析了繞組軸向預(yù)緊力與器身軸向固有頻率之間的關(guān)系。但是，該等效模型無法準(zhǔn)確分析繞組幾何結(jié)構(gòu)等對其振動的影響，導(dǎo)致計算結(jié)果與繞組的實際振動特性存在一定偏差。第二類模型是隨著有限元分析法的快速發(fā)展而產(chǎn)生的。文獻(xiàn)[7-8]通過變壓器繞組有限元模型，計算分析了不同預(yù)緊力下變壓器繞組的一階固有頻率及其對應(yīng)的振型。文獻(xiàn)[9-10]根據(jù)換流變繞組的三維有限元模模型，計算分析了變壓器在穩(wěn)態(tài)和短路條件下的繞組電磁力特性及其激勵下的繞組振動及固有頻率。文獻(xiàn)[11]計算了變壓器繞組在漏磁場作用下的電動力特性，然后將短路電動力作為繞組軸向振動載荷分析了繞組的穩(wěn)定性，得出了繞組各個線餅的位移分布。但現(xiàn)有研究大都對變壓器繞組正常狀態(tài)下的振動特性進(jìn)行計算分析，對繞組本身的固有振動特性及其隨繞組狀態(tài)改變的變化規(guī)律尚不清楚。尤其是繞組松動時，變壓器的電參數(shù)如短路電抗、分布電感和分布電容等對其欠敏感，需要尋求其他有效的檢測手段，如振動檢測法等。

本文嘗試建立變壓器繞組的三維有限元模型，計算分析變壓器繞組正常及松動狀態(tài)下軸向固有振動特性及其變化規(guī)律。同時，對實體變壓器繞組為對象進(jìn)行模態(tài)特性測試，用以修正和完善仿真模型，期望獲得更為準(zhǔn)確的變壓器繞組軸向振動特性，為變壓器繞組結(jié)構(gòu)優(yōu)化設(shè)計和繞組變形振動檢測法提供依據(jù)。

1 繞組固有振動特性的有限元建模

1.1 理論分析

變壓器繞組是一個典型的復(fù)雜多自由度機械系統(tǒng)，而繞組軸向振動是電磁力激勵下的強迫振動，滿足如下動力學(xué)方程，為

式（7）中，等號左側(cè)二項分別為單元剛度矩陣和單元質(zhì)量矩陣。若將各個單元坐標(biāo)變換后的單元質(zhì)量矩陣和單元剛度矩陣組集，便可根據(jù)達(dá)朗貝爾原理得到如式（2）所示的繞組結(jié)構(gòu)動力學(xué)方程。

1.2 變壓器繞組的有限元建模

以某10kV變壓器高壓繞組為研究對象進(jìn)行有限元建模，該繞組結(jié)構(gòu)為線餅式結(jié)構(gòu)。因主要計算分析變壓器繞組的軸向振動特性，故根據(jù)變壓器繞組結(jié)構(gòu)尺寸在有限元建模時采用如下方法：

1）將繞組等效為一整塊環(huán)狀餅行實體，材質(zhì)為銅，分別設(shè)置了彈性模量、泊松比和密度等參數(shù)，由實際測試得到。繞組首末兩端各設(shè)置一圈剛性擋板，主要參數(shù)為彈性模量和泊松比，根據(jù)變壓器繞組預(yù)緊力的設(shè)置值和變壓器繞組模態(tài)測試結(jié)果等效及調(diào)整得到。

2）絕緣墊塊模型根據(jù)繞組實際結(jié)構(gòu)尺寸建立，以長方實體結(jié)構(gòu)等效。因其彈性模量隨預(yù)緊力變化而變化，故將其視為非線性材料。根據(jù)墊塊的實測應(yīng)力—應(yīng)變曲線進(jìn)行擬合得到其彈性模量變化曲線。

3）線圈、墊塊均選用solid 45單元類型；

4）對繞組上下兩端實行全約束，主要考慮軸向自由度。

根據(jù)變壓器繞組實際結(jié)構(gòu)，整個繞組模型共66層線餅，每層線餅間均勻分布了8塊墊塊。對所建立的變壓器繞組有限元模型進(jìn)行掃掠剖分后，共得到38 448個單元，如圖1所示。

2 變壓器繞組模態(tài)測試描述

實驗對象為一臺10 kV變壓器高壓繞組，其實物圖片如圖2所示，主要對其軸向模態(tài)特性進(jìn)行了測試。

圖3為模態(tài)測試分析系統(tǒng)圖。圖中，信號發(fā)生器用于產(chǎn)生20 kHz的帶寬白噪聲信號，經(jīng)過功率放大器進(jìn)行放大后，驅(qū)動垂直懸掛的激振器對變壓器繞組進(jìn)行激勵。其中，激振器端部細(xì)桿末端裝有力傳感器，用于采集激勵力的大小。數(shù)據(jù)系統(tǒng)采集激振器力信號與振動加速度傳感器信號，計算出各測點處的頻響函數(shù)。采用單點激勵多點拾振法對變壓器繞組進(jìn)行模態(tài)測試，分批次進(jìn)行。其中，加速度傳感器被固定在金屬夾件上，沿軸向布置5層，每層4個，共20個測點。

圖4為模態(tài)實驗得到的變壓器繞組振動頻響函數(shù)（vibration frequency response function，VFRF）。為清晰起見，縱坐標(biāo)采用對數(shù)坐標(biāo)表示。由圖可見，振動頻響函數(shù)的形狀隨頻率變化有著較為明顯的峰值，且均在低頻部分存在較大的噪聲干擾?？紤]到，多參考點最小二乘復(fù)頻域法（polyreference least-square complex frequency domain method， PolyMAX）特別適用于大阻尼和密集模態(tài)結(jié)構(gòu)模態(tài)參數(shù)的識別[12]，故本文在此選用該算法對變壓器繞組的模態(tài)參數(shù)進(jìn)行識別。

3 結(jié)果分析

3.1 繞組正常工況下的固有振動特性

表1分別為計算和根據(jù)繞組模態(tài)測試結(jié)果得到的變壓器繞組前四階固有頻率。由表可見，繞組前四階固有頻率的計算結(jié)果與試驗結(jié)果的偏差均在6%以內(nèi)，說明了計算結(jié)果的正確性。出現(xiàn)偏差的主要原因可能是由于實際繞組在端圈及各種夾件固定下使得仿真模型的邊界條件與實際情況存在差異。此外，變壓器繞組的前四階固有頻率均遠(yuǎn)離100 Hz的電動力激勵頻率，不易引發(fā)共振。

圖5和圖6分別為變壓器繞組前四階固有頻率對應(yīng)振型的計算結(jié)果與實驗結(jié)果。

圖6中，黑色實線部分為靜止的繞組位置，紅色虛線部分為運動中的繞組位置。由圖可見，前四階固有頻率對應(yīng)振型的仿真結(jié)果與實測結(jié)果有著較好的一致性，其中，一階振型均表現(xiàn)為繞組的整體上下振動，二階振型為繞組的左右“蹺板”式振動，三階振型為繞組兩端翹起、中間落下的形態(tài)，四階振型為繞組兩端向中間部分?jǐn)D壓的形態(tài)。顯然，這一結(jié)果進(jìn)一步說明了本文所建立的有限元仿真模型可以較為準(zhǔn)確地對變壓器繞組的機械特性進(jìn)行計算分析。

3.2 繞組預(yù)緊力變化時的固有振動特性

考慮到繞組松動（預(yù)緊力下降）是在運變壓器的常見故障之一，且如何尋找合理的變壓器繞組預(yù)緊力范圍一直是變壓器制造廠關(guān)心的焦點之一，本文在此根據(jù)所建立的有限元仿真模型，重點分析了繞組預(yù)緊力變化時的固有振動特性。

已有的研究表明，繞組線餅間的絕緣墊塊是一種剛度漸硬材料，其剛度系數(shù)隨繞組預(yù)緊力的大小有較明確的變化。由文獻(xiàn)[8]知，變壓器繞組絕緣墊塊的彈性模量E可表示為

顯然，墊塊的彈性模量隨壓強的增大而增大。故可通過改變有限元仿真模型中墊塊材料彈性模量模擬預(yù)緊力改變時的繞組固有振動特性。

圖7為變壓器繞組前四階固有頻率隨預(yù)緊力改變的變化曲線。圖中，F(xiàn)N為變壓器繞組預(yù)緊力設(shè)計值。由圖可見，當(dāng)變壓器繞組預(yù)緊力降低即繞組松動時，繞組前四階固有頻率向低頻方向偏移；預(yù)緊力增大時，繞組前四階固有頻率向高頻方向偏移。其中，繞組預(yù)緊力在80%FN～120%FN變化時，繞組前三階固有頻率的變化率約為20%；繞組預(yù)緊力在40%FN～80%FN變化時，繞組前三階固有頻率的變化率約為30%。此外，繞組第四階固有頻率較前三階固有頻率的變化更加明顯。因此，當(dāng)有松動隱患的變壓器繞組遭受頻率分量較為豐富的激勵力作用，例如突發(fā)短路的情形，若激勵力的頻率分量與繞組的固有頻率較為接近時，有可能會引發(fā)變壓器繞組共振而引起繞組坍塌，導(dǎo)致變壓器故障。同時，對變壓器制造廠來說，若已知變壓器繞組預(yù)緊力設(shè)計值時對應(yīng)的固有頻率，則可從繞組固有頻率的變化程度對繞組的壓緊程度進(jìn)行預(yù)估，從而指導(dǎo)變壓器生產(chǎn)實際。

限于論文篇幅，本文在此給出了變壓器繞組預(yù)緊力為80%FN時其前四階固有頻率對應(yīng)的振型，如圖8所示。

此時，繞組的前四階固有頻率分別為225 Hz、301 Hz、480 Hz和658 Hz。由圖8可見，繞組松動時，一階振型均表現(xiàn)為繞組的整體上下振動，二階振型為繞組的左右“蹺板”式振動，三階振型為繞組兩端翹起、中間落下的形態(tài)，四階振型為繞組兩端向中間部分?jǐn)D壓的形態(tài)，與繞組正常狀態(tài)的振型形態(tài)類似，并具有較好的對稱性。但是，繞組固有頻率的變化會使得變壓器繞組在類似激勵力下的固有振動特性發(fā)生改變。相應(yīng)地，繞組振動信號也會隨之發(fā)生改變。

4 結(jié) 論

1）變壓器繞組正常狀態(tài)下前四階固有頻率和對應(yīng)的振型與繞組激振實驗結(jié)果吻合良好，說明了所建立的有限元模型的正確性。應(yīng)用該模型，可以方便地對變壓器繞組的固有振動特性進(jìn)行分析。

2）變壓器繞組正常狀態(tài)下的前四階固有頻率所對應(yīng)的振型對稱性較好。繞組松動時，前四階固有頻率均隨預(yù)緊力的降低而降低，但變化幅度存在差異。預(yù)緊力增大時，前四階固有頻率也隨之增大。

3）繞組預(yù)緊力改變時引發(fā)的繞組固有振動特性的改變必然會使得繞組振動信號發(fā)生改變，尤其需要關(guān)注突發(fā)短路時有松動隱患的變壓器因固有頻率向低頻偏移而引發(fā)的共振問題。

本文所用仿真和模態(tài)測試手段同樣可推及更高電壓等級的變壓器繞組，考慮到變壓器繞組故障的多樣性和復(fù)雜性，對不同故障類型及故障程度下的變壓器繞組固有振動特性進(jìn)行計算分析是我們下一步的研究方向。

參 考 文 獻(xiàn)：

[1] 顧曉安，沈榮滋，徐基泰. 國外變壓器噪聲研究的動向[J]. 變壓器，2002，39（6）：33.

GU Xiaoan， SHEN Baozi， XU Jitai.Trend of transformer noise research abroad [J]. Transformer， 2002， 39（6）： 33.

[2] WANG F H， JIN Z J. Using the vibration frequency response analysis method to detect the winding deformation of power transformer[C]//IEEE PES General Meeting， Detroit， U. S. A， July 2011.

[3] SWIHART D O， WRIGHT D V. Dynamic stiffness and damping of transformer pressboard during axial short circuit vibration [J]. IEEE Transactions on Power Apparatus and Systems， 1976， 95（2）： 721.

[4] 王洪方，王乃慶，李同生.變壓器繞組軸向預(yù)緊力對繞組軸向振動特性的影響[J].電網(wǎng)技術(shù)，1999，23（9）：8.

WANG Hongfang， WANG Naiqing， LI Tongsheng. Influence of axial pre-compression level on axial vibrations in transformer windings[J]. Power System Technology， 1999， 23（9）： 8.

[5] 王洪方，王乃慶，李同生.大型電力變壓器繞組軸向非線性振動研究[J].電網(wǎng)技術(shù)，2000，24（3）：42.

WANG Hongfang， WANG Naiqing， LI Tongsheng.Axial nonlinear vibration of large power transformer winding[J]. Power System Technology， 2000， 24（3）： 42.

[6] 王世山，汲勝昌，李彥明.電纜繞組變壓器短路時線圈軸向穩(wěn)定性的研究[J].中國電機工程學(xué)報，2004， 24（2）： 166.

WANG Shishan， JI Shengchang， LI Yanming. Study on axial stability in condition of short-circuit for power transformer using XLPE insulated cable windings [J]. Proceedings of the CSEE， 2004， 24（2）： 166.

[7] 李巖，周偉，井永騰，等.短路條件下變壓器器身軸向振動的分析[J]，電工電能新技術(shù)，2012，31（3）：49.

LI Yan， ZHOU Wei， JING Yongteng， et al. Axial vibration analysis of transformer active part under short-circuit [J]. Advanced Technology of Electrical Engineering and Energy， 2012， 31（3）： 49.

[8] 李洪奎，李巖.不同預(yù)緊力下變壓器繞組軸向振動模態(tài)分析[J].電機與控制學(xué)報，2010，14（8）：98.

LI Hongkui， LI Yan. Axial vibration modal analysis of transformer windings under different level of pre-compression [J]. Electric Machines and Control， 2010， 14（8）： 98.

[9] 許加柱，羅隆福，李勇，等.新型換流變壓器繞組電磁力的分析計算[J]，高電壓技術(shù)，2007，33（6）：102.

XU Jiazhu， LUO Longfu， LI Yong， et al. Analysis and calculation of windings electromagnetic force of a novel converter transformer[J]. High Voltage Engineering， 2007， 33（6）： 102.

[10] SHAO Pengfei， LUO Longfu， LI Yong， et al.Electromagnetic vibration analysis of the winding of a new HVDC converter transformer[J]. IEEE Transactions on Power Delivery， 2012， 27（1）： 123.

[11] 郭健，林鶴云，徐子宏，等.用有限元方法分析電力變壓器繞組軸向穩(wěn)定性[J]，高電壓技術(shù)，2007，33（11）：209.

GUO Jian， LIN Heyun， XU Zihong， et al. Analysis of axial stability of power transformer windings using finite element [J]. High Voltage Engineering， 2007， 33（11）： 209.

[12] 孫鑫暉，郝木明，王淮維. PolyMAX模態(tài)參數(shù)識別算法的快速實現(xiàn)[J]. 振動與沖擊，2011，30（10）：5.

SUN Xinhui， HAO Muming， WANG Huaiwei. Fast implementation for PolyMAX modal identification algorithm[J]. Journal of Vibration and Shock， 2011， 30（10）： 5.

（編輯：劉素菊）