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

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

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

中圖分類號(hào)：TM 401.1

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

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

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

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

1 繞組固有振動(dòng)特性的有限元建模

1.1 理論分析

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

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

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

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

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

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

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

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

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

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

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

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

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

3 結(jié)果分析

3.1 繞組正常工況下的固有振動(dòng)特性

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

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

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

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

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

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

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

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

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

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

4 結(jié) 論

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

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

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

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

參 考 文 獻(xiàn)：

[1] 顧曉安，沈榮滋，徐基泰. 國(guó)外變壓器噪聲研究的動(dòng)向[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ù)緊力對(duì)繞組軸向振動(dòng)特性的影響[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] 王洪方，王乃慶，李同生.大型電力變壓器繞組軸向非線性振動(dòng)研究[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] 王世山，汲勝昌，李彥明.電纜繞組變壓器短路時(shí)線圈軸向穩(wěn)定性的研究[J].中國(guó)電機(jī)工程學(xué)報(bào)，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] 李巖，周偉，井永騰，等.短路條件下變壓器器身軸向振動(dòng)的分析[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ù)緊力下變壓器繞組軸向振動(dòng)模態(tài)分析[J].電機(jī)與控制學(xué)報(bào)，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ì)算[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ù)識(shí)別算法的快速實(shí)現(xiàn)[J]. 振動(dòng)與沖擊，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.

（編輯：劉素菊）