景乾明，胡春梅，胡基才，張承學(xué)

(1.甘肅省電力公司，蘭州市， 730050;2.天水超高壓輸變電公司，甘肅省天水市， 741000; 3.武漢大學(xué)，武漢市， 430072)

0 引言

覆冰輸電導(dǎo)線在一定冰風(fēng)條件下所誘發(fā)的舞動是一種低頻、大振幅的自激振動［1］。導(dǎo)線舞動容易導(dǎo)致系統(tǒng)短路和連續(xù)跳閘、導(dǎo)線燒傷、斷股和斷線、各種線路金具以及桿塔的損壞［2］。隨著我國電網(wǎng)覆蓋面的不斷擴展以及電壓等級的不斷提高，穿越復(fù)雜地形及惡劣氣候條件的送電線路日益增多［3］。由于山區(qū)尤其是高海拔山區(qū)的局部特殊地形地貌不同，使得其局部氣象具有明顯的立體特征，即使在一個小范圍內(nèi)，不同地形的局部氣象指標(biāo)也存在極大的差異［4］。這些差異將直接導(dǎo)致多變的局部氣流條件和不同的導(dǎo)線覆冰［5］，使得按一般條件所設(shè)計的輸電線路在這些地區(qū)會經(jīng)常產(chǎn)生嚴(yán)重舞動，造成巨大的工程事故隱患。由于導(dǎo)線舞動的復(fù)雜性和特殊性，使用計算機仿真技術(shù)研究不同地形地貌條件下的導(dǎo)線舞動具有特殊優(yōu)勢。目前國內(nèi)外對導(dǎo)線舞動的仿真研究還處于起步階段，而運用計算機仿真技術(shù)對微地形條件下的導(dǎo)線舞動研究較少［6］，因此采用仿真技術(shù)研究微地形條件下的輸電線路舞動具有重要意義。

1 微地形及其對空間流場的影響

微地形對輸電線路舞動的影響主要體現(xiàn)在2個方面:影響局部空氣流場和導(dǎo)線覆冰。

1.1 微地形對空間流場的影響

甘肅電網(wǎng)某330 kV線路的66～67檔位于一個峽谷山坡的頂部附近，2桿塔導(dǎo)線懸掛點的高差約為40 m。特殊的地形地貌將直接改變導(dǎo)線周圍的空間流場，包括風(fēng)速、風(fēng)向和氣流性質(zhì)等。利用計算流體力學(xué)軟件可以建立相應(yīng)的流場模型。圖1為由計算機仿真計算所求得的該微地形條件下的空間流場分布。

計算結(jié)果表明，當(dāng)風(fēng)速為15 m/s的水平層流風(fēng)經(jīng)過該微地形時，風(fēng)場產(chǎn)生了嚴(yán)重的紊流現(xiàn)象。在導(dǎo)線截面所在的不同空間位置，各點的風(fēng)速和風(fēng)向均與入口處不同。覆冰導(dǎo)線中點處的實際風(fēng)速約為21 m/s，而風(fēng)向與水平面的夾角約為15°。

圖1 微地形附近流場截面Fig.1Cross section of flow field near the micro-terrain

1.2 微地形對導(dǎo)線覆冰的影響

微地形導(dǎo)致導(dǎo)線周圍風(fēng)速、風(fēng)向等發(fā)生改變，從而影響到導(dǎo)線的覆冰。一般導(dǎo)線的覆冰與導(dǎo)線接觸過冷卻水滴的有效面積相關(guān)，導(dǎo)線迎風(fēng)面凍結(jié)相對較大，而背風(fēng)面則相對較小;同時冰結(jié)厚度與風(fēng)向和導(dǎo)線之間夾角的正弦成正比，即當(dāng)風(fēng)向垂直于導(dǎo)線軸線時凍結(jié)厚度增長最快［7］。此外，導(dǎo)線的覆冰還與導(dǎo)線結(jié)構(gòu)、溫度、海拔、大氣濕度、氣壓等因素有關(guān)［8］。

根據(jù)圖1所示的流場計算結(jié)果以及對該檔線路的實際觀測，該檔線路的實際覆冰形狀為霧凇性質(zhì)的D型覆冰。

2 導(dǎo)線舞動的仿真模型

2.1 導(dǎo)線動力學(xué)模型

考慮豎直位移、水平位移和扭轉(zhuǎn)位移3個自由度，建立覆冰導(dǎo)線三自由度集中參數(shù)動力學(xué)模型［9］，如圖2所示。

圖2 三自由度集中參數(shù)系統(tǒng)模型Fig.23-D lumped-parameter system model

對圖2所示模型建立其垂直(y向)、水平(z向)及扭轉(zhuǎn)振動的綜合動力學(xué)運動方程［10］為

式中:U為風(fēng)速;D為導(dǎo)線迎風(fēng)尺寸;ρ為空氣密度; α為攻角;為扭轉(zhuǎn)振動的角速度和角加速度;y和z為豎直和水平自由度;為豎直、水平方向的速度及加速度;CL、CD、CM為升力、阻力、扭矩系數(shù);Cy為Y向空氣動力系數(shù);ξy、ξz、ξt為模型在3個方向上的結(jié)構(gòu)阻尼率;ωy、ωz、ωt為模型在3個方向上得振動固有頻率;ky、kz、kt為模型在3個方向上的剛度;mi為單位長度導(dǎo)線上的覆冰質(zhì)量;β0為初始凝冰角;r為導(dǎo)線半徑;R為特性半徑，本文可取為半徑r。

2.2 導(dǎo)線構(gòu)型

導(dǎo)線的構(gòu)型方程［11］為

式中:x為導(dǎo)線上某點與左懸掛點的水平距離;y為導(dǎo)線上某點與左懸掛點的垂直距離;l為檔距;γ為單位長度與截面積上承受的載荷;σ0為導(dǎo)線水平應(yīng)力;h為懸掛點高度差;Lh=0為兩側(cè)懸掛點等高時檔內(nèi)懸鏈線的總線長，Lh=0=(2σ0/γ)sinh(rl/2σ0)。

3 微地形下的導(dǎo)線舞動及仿真分析

二分裂330 kV線路66～67檔的結(jié)構(gòu)參數(shù)如下:檔距為203 m;兩桿塔絕緣子海拔差為39 m;導(dǎo)線型號為LGJ300/ 40;平均運行張力為10.096 kN。

對該線檔進行模態(tài)分析［12］，求得該系統(tǒng)前4階固有頻率結(jié)果如表1所示。

表1 前4階固有頻率Tab.1Natural frequency of fourth-order

在2011年2月的實際觀察中，發(fā)現(xiàn)該檔導(dǎo)線發(fā)生舞動，舞動振型為二階振型，舞動振幅的最大峰峰值約為3 m?，F(xiàn)場觀測該檔導(dǎo)線為D型覆冰，覆冰厚度為10 mm，現(xiàn)場平均風(fēng)速約為15 m/s。

3.1 考慮微地形影響的仿真計算

根據(jù)該檔線路實際微地形條件，建立導(dǎo)線系統(tǒng)的有限元模型，并利用所開發(fā)的“輸電導(dǎo)線舞動仿真計算平臺系統(tǒng)”進行動力學(xué)仿真計算，計算結(jié)果如圖3所示。

圖3 考慮微地形的仿真計算振型Fig.3Vibration mode of the simulation calculation with micro-terrain

圖4為距左塔50 m處的Y向振動時間歷程仿真曲線及頻譜圖，導(dǎo)線在600和1 200 s時的截面舞動軌跡如圖5所示。

圖4 Y向振動時間歷程及頻譜圖Fig.4Time history and spectrogram of the vibration in Y direction

從圖3、4的仿真結(jié)果可以看出，該檔導(dǎo)線舞動過程中Y向的最大振幅約為2.7 m，發(fā)生在距左塔50 m處，與圖4的實際觀測結(jié)果相近。該檔導(dǎo)線的Y向振動是主頻約為0.45 Hz的二階振動，這與表1的理論分析十分接近，也和所觀測到的二階舞動現(xiàn)象一致。

從圖5的仿真結(jié)果可以看出，導(dǎo)線在舞動過程中的截面軌跡大體上為一斜橢圓狀，與導(dǎo)線舞動的理論分析結(jié)論相一致［1］。

圖5 導(dǎo)線舞動的軌跡曲線Fig.5Locus curve of the conductor galloping

上述分析表明，該仿真結(jié)果在舞動模態(tài)、舞動振幅以及舞動軌跡方面均與理論分析和實際觀測基本相符，這表明該計算機仿真平臺系統(tǒng)所得的仿真結(jié)果可信度較高。

3.2 忽略微地形影響的仿真計算

如果忽略該處的微地形結(jié)構(gòu)，即認(rèn)為該檔導(dǎo)線只是位于海拔不同的2個山坡上，在相同的水平風(fēng)速、相同的覆冰厚度、相同的線路結(jié)構(gòu)的前提下，運用同樣的方法對該檔導(dǎo)線系統(tǒng)進行動力學(xué)仿真計算，結(jié)果如圖6所示。

圖6 忽略微地形的導(dǎo)線舞動振型Fig.6Vibration mode without micro-terrain

從圖6可以看出，不考慮微地形條件下導(dǎo)線的舞動振型為一階振型，舞動振幅最大處位于檔距中點附近位置，最大振幅約為5 m。此時檔距中點處的Y向振動時間歷程和頻譜圖如圖7所示，導(dǎo)線中點在600和1 200 s時的截面舞動軌跡如圖8所示。

圖6、7的仿真結(jié)果表明，此時該檔導(dǎo)線舞動過程中Y向的最大振幅約為5 m，發(fā)生在檔距中點處。此時該檔導(dǎo)線的Y向振動主頻約為0.22和0.45 Hz，這與表1中導(dǎo)線的一階固有頻率0.21 Hz和二階固有頻率0.42 Hz接近。表明導(dǎo)線同時存在一階和二階振動成分，同時由于一階舞動的振幅大于二階舞動的振幅，所以導(dǎo)線的振型包絡(luò)線主要呈現(xiàn)為一階的形狀。

圖8的仿真結(jié)果表明，此時導(dǎo)線在舞動過程中的截面軌跡更接近理論上的橢圓狀，這是由于忽略了微地形后，系統(tǒng)所在的流場和流場作用在導(dǎo)線上的外載荷都更加均勻所致。

圖3～8之間的對比表明，微地形的考慮與否對導(dǎo)線舞動的模態(tài)、振幅以及截面軌跡等有著顯著的影響，因而將對線路系統(tǒng)的結(jié)構(gòu)設(shè)計、防舞設(shè)計、強度計算以及疲勞壽命評估等帶來顯著影響。

4 結(jié)論

(1)在理論分析的基礎(chǔ)上，建立了覆冰導(dǎo)線的三自由度分析模型和有限元計算模型，并對某實際檔距進行了實際仿真計算。仿真計算結(jié)果不僅與理論分析相符，而且與實際舞動觀測結(jié)果相近，表明該分析模型和仿真計算可信。

(2)利用計算流體動力學(xué)仿真軟件建立了某實際線路微地形結(jié)構(gòu)的流場仿真計算模型，并進行了仿真計算。計算結(jié)果表明微地形對風(fēng)速和風(fēng)向都有不同程度的影響，而且依地形變化程度的不同相異較大，實際處理中應(yīng)根據(jù)線路實際情況分析求解。

(3)影響輸電導(dǎo)線覆冰的因素很多，應(yīng)當(dāng)綜合考慮。微地形導(dǎo)致導(dǎo)線周圍流場的變化對輸電導(dǎo)線的覆冰形狀、覆冰方向、覆冰厚度以及覆冰過程都有顯著的影響，但在實際處理中要具體情況具體分析。

(4)實際線路的仿真計算表明，微地形不僅對導(dǎo)線的舞動振幅和截面運動軌跡有較大的影響，而且還可能影響到導(dǎo)線的舞動振型。

(5)微地形對導(dǎo)線舞動影響的結(jié)論為線路設(shè)計部門和線路運行管理單位在線路的設(shè)計和防舞設(shè)計時提供了有益的指導(dǎo)。結(jié)論表明針對不同的地形地貌和氣候條件，在進行線路設(shè)計時應(yīng)慎重選擇線路的走向并分別采用不同的檔距等結(jié)構(gòu)參數(shù);對運行中有舞動隱患的線路在進行防舞改造設(shè)計時，應(yīng)針對性地采取不同的防舞方案設(shè)計和參數(shù)設(shè)計。

［1］郭應(yīng)龍，李國興，尤傳永.輸電導(dǎo)線舞動［M］.北京:中國電力出版社，2002.

［2］沈巒峰，巫世晶，王曉筍，等.基于扭振激發(fā)機理的輸電導(dǎo)線新型防舞模型［J］.電力建設(shè)，2010，31(3):34-37.

［3］王守禮，張弦.微地形微氣象對輸電線路的影響及響應(yīng)措施［J］.云南電業(yè)，2005(6):36-37.

［4］Zhou Q，Yang Z S，Chen W G． Design of conductor galloping multipoints monitoring system using ZigBee PRO wireless network technology［J］． High Voltage Engineering，2011，37( 8) : 1967-1974．

［5］蔣興良，周仿榮，王少華，等.輸電導(dǎo)線覆冰舞動機理及防治措施［J］.電力建設(shè)，2008，29(9):14-18.

［6］嚴(yán)波，胡景，周松，等.覆冰四分裂導(dǎo)線舞動數(shù)值模擬及參數(shù)分析［J］.振動工程學(xué)報，2010，23(3):310-315.

［7］魏其巍.電線覆冰機理分析及在工程設(shè)計中的應(yīng)用［J］.電力建設(shè)，2007，28(3):26-28.

［8］金西平.微地形微氣候?qū)﹄娏€路覆冰的影響［J］.供用電，2008，25(4):17-20.

［9］李欣業(yè)，張華彪，高仕趙，等.三自由度模型覆冰輸電導(dǎo)線的數(shù)值仿真分析［J］.河北工業(yè)大學(xué)學(xué)報，2010，39(3):1-5.

［10］李黎，陳元坤，夏正春，等.覆冰導(dǎo)線舞動的非線性數(shù)值仿真研究［J］.振動與沖擊，2011:30(8):107-111.

［11］石端偉.機械動力學(xué)［M］.北京:中國電力出版社，2007.

［12］Li L，Chen Y K，Xia Z C.Nonlinear numerical simulation study of iced conductor galloping［J］.Journal of Vibration and Shock，2011，30(8):107-111.

(編輯:張磊)