劉玥君，張新語，郭峻菘，張 猛

(1.東北電力大學(xué) 建筑工程學(xué)院，吉林 吉林 132012；2.吉林市城鄉(xiāng)規(guī)劃研究院，吉林 吉林 132001)

架空輸電線路大多建于偏遠地區(qū)，氣候較為惡劣，輸電線路在冰風荷載作用下遭受的事故屢見不鮮，因此研究輸電線路在冰風荷載作用下的可靠性問題是非常必要的.國外對輸電線路在冰風荷載作用下的可靠性研究起步比較早.A.Y.Shehata等采用有限元軟件對強風作用下輸電塔-線體系的動力響應(yīng)進行了深入分析，為改進輸電塔-線體系的抗風設(shè)計和應(yīng)用提出了有效建議[1].Yang F等人利用有限元模型在考慮覆冰剛度的情況下進行導(dǎo)線斷裂試驗，結(jié)果表明模擬值與試驗值基本一致[2].Li Qin，Wei Lin-peng，Wei Han等利用RSM(響應(yīng)面法)對大跨越輸電塔-線體系進行可靠性研究[3].S.C.Yang，T.J.Liu，H.P.Hong等利用改進的Monte-Carlo法分析了結(jié)構(gòu)體系的可靠性[4].相對于國外，國內(nèi)對輸電線路在冰風荷載下的可靠性研究起步較晚.屈成忠，張恒濤等以廣東某輸電線路為背景，建立塔線體系有限元模型，對鐵塔進行抗風能力評估和穩(wěn)定性分析[5].李妍等以某高壓輸電塔為基礎(chǔ)，針對輸電塔-線隱式非線性動力可靠度進行分析[6].侯景鵬應(yīng)用ANSYS有限元分析軟件建立了500 kV貓頭直線塔，分析比較了該模型在靜力風荷載和脈動風荷載作用下的風振響應(yīng)[7].趙桂峰與李杰等通過有限元軟件ANSYS模擬了非線性風振響應(yīng)過程，并與實測的數(shù)據(jù)進行對比分析[8].俞登科等通過引入大矩陣計算出輸電塔-線體系的整體可靠度[9].李正良，施菁華等以哈密-鄭州800 kV特高壓直流輸電工程中的一段塔-線耦合體系為研究對象，開展了一系列的風洞試驗研究[10].楊文剛，王璋奇等應(yīng)用Newmark法，研究了在脈動風荷載作用下，我國特高壓輸電塔-線體系的風振響應(yīng)[11].史天如等以220 kV貓頭型直線塔作為模擬實例，根據(jù)系統(tǒng)的分析和計算結(jié)果的可靠性及線路所在地區(qū)的氣候條件給出了對應(yīng)防災(zāi)措施的建議[12].目前我國關(guān)于輸電線路在冰風荷載作用下相關(guān)的研究都非常匱乏，因此，開展對輸電線路在冰風荷載作用下的可靠性研究具有重要意義，并在此基礎(chǔ)上分析目前覆冰輸電塔在其設(shè)計建設(shè)以及維護管理過程中的一些不足之處，以便更好地理論結(jié)合實際確保大電網(wǎng)安全可靠運行.

1 覆冰輸電塔-線耦合體系模型的建立

本文研究主塔型號為ZM2k型貓頭塔，導(dǎo)線型號為4×LGJ-630/45鋼芯鋁絞線，絕緣子采用陶瓷絕緣子，包含12片絕緣子，單片絕緣子長155 mm，絕緣子串總長1.86 m.輸電塔全高45 m，呼高42 m，水平檔距400 m，設(shè)計電壓等級220 kV.通過ANSYS有限元軟件建立輸電塔-線耦合體系模型，如圖1所示.

圖1 塔-線耦合體系有限元模型

為了驗證輸電塔-線建模合理，采用模態(tài)分析法對其進行驗證.提取輸電塔-線耦合體系部分振型如圖2所示.由各階振型圖可以看到，輸電塔-線耦合體系的自振頻率隨著階數(shù)的增加而變大.在低階時，輸電塔-線耦合體系的振型往往是一條導(dǎo)線或者地線的振動；在高階時，輸電塔-線耦合體系往往是兩條導(dǎo)線或者一條地線和一條導(dǎo)線的振動.

圖2 輸電塔-線體系振型圖

應(yīng)用改變密度法分別模擬5 mm、10 mm、15 mm、20 mm、25 mm和30 mm厚度的覆冰，分析出輸電線路的材料屬性是隨著覆冰厚度增加不斷變化的.覆冰厚度增加的同時輸電塔與導(dǎo)地線的等效直徑、等效截面積、等效體積和總質(zhì)量都是增大的，但等效密度是隨著覆冰厚度減小的.為了簡化有限元模型，假設(shè)輸電-塔線體系各部分均勻覆冰5 mm、10 mm、15 mm、20 mm、25 mm、30 mm，取冰的密度為0.9 g/cm3，由于輸電線路的檔距遠遠大于導(dǎo)線的直徑，所以在此不考慮剛度對輸電線路的影響.覆冰后輸電塔-線耦合體系的自振頻率，如圖3所示.由圖3可知，隨著覆冰厚度的增加，輸電塔線-體系的自振頻率在逐漸減小，且在覆冰厚度由10 mm增加到15 mm時變化最大.

圖3 覆冰輸電塔-線體系自振頻率圖

2 脈動風荷載的數(shù)值模擬

選取輸電塔-線體系上的風速模擬點，如圖4所示.采用諧波合成法模擬脈動風荷載，輸電塔從上到下取11個點，在導(dǎo)地線上分別選取風速模擬點，兩跨導(dǎo)線上每支導(dǎo)線選取12個風速模擬點，同樣的，地線上也要選取12個模擬點，故在輸電塔-線結(jié)構(gòu)體系上共模擬了47個風速點.

圖4 輸電塔-線結(jié)構(gòu)體系風速模擬點分布圖

根據(jù)我國規(guī)范，設(shè)輸電線路在高度10 m處的平均風速為V10=5 m/s、V10=10 m/s、V10=15 m/s三種工況，地面粗糙度系數(shù)k=0.005，脈動風速時間間隔為0.1 s，頻率步長為0.01 Hz，選擇Davenport譜模擬輸電塔-線耦合體系上選定的47個脈動風速模擬點，生成了200 s的風速時程數(shù)據(jù)與圖像.以10 m高度處平均風速V10=10 m/s為例，可以得到輸電塔-線耦合體系在該地區(qū)高度為10 m處的脈動風速圖、風速自譜圖、合成風速圖，分別如圖5、圖6、圖7所示.并將模擬的風速功率譜與Davenport目標功率譜曲線對比可知，兩條曲線趨勢一致且數(shù)值相近，論證了風速時程模擬的正確性.

圖5 脈動風速譜圖圖6 風速自譜圖

圖7 合成脈動風速圖

3 覆冰輸電塔-線耦合體系在風荷載作用下的響應(yīng)分析

3.1 塔頂位移時程響應(yīng)分析

輸電塔作為高聳結(jié)構(gòu)，其結(jié)構(gòu)的最高點在冰風荷載作用下的最大響應(yīng)位移應(yīng)該是整個結(jié)構(gòu)所有節(jié)點中最大的，于是以輸電塔最高點523號節(jié)點進行分析.通過在標準高度10 m處平均風速V10=5m/s、V10=10 m/s和V10=15 m/s，輸電-塔線體系各部分均勻覆冰0 mm、5 mm、10 mm、15 mm、20 mm、25 mm和30 mm，共計二十一種組合工況下模擬冰風荷載，進而生成位移時程曲線進行對比分析.

將模擬得到的風速時程轉(zhuǎn)化為節(jié)點的風壓時程，通過ANSYS有限元軟件進行風振響應(yīng)的模擬.在標準高度10 m處平均風速V10=5 m/s、V10=10 m/s和V10=15 m/s三種風荷載與不同覆冰厚度組合作用下，523號節(jié)點最大位移如圖8所示.

圖8 523號節(jié)點最大位移圖

由圖8可知，523號節(jié)點的位移風振響應(yīng)隨著風速的增加和覆冰厚度的增加而逐漸變大，在10 m高度處風速V10=15 m/s、覆冰厚度H=30 mm時達到位移最大值，其值為0.292 8 m.在風速相同的情況下，輸電塔-線體系在覆冰厚度從10 mm增加到15 mm時，輸電塔架位移增加值最大，最高可增加0.061 8 m.

3.2 導(dǎo)地線位移時程響應(yīng)分析

在標準高度10m處的平均風速為V10=5 m/s、V10=10 m/s和V10=15 m/s三種工況模擬的脈動風荷載作用下，不同覆冰厚度下980號節(jié)點與1056號節(jié)點最大位移，如圖9、圖10所示.

圖9 980號節(jié)點最大位移圖圖10 1506號節(jié)點最大位移圖

由圖9、圖10中也可以看到，980號節(jié)點與1506號節(jié)點的位移風振響應(yīng)隨著風速的增加和覆冰厚度的增加而逐漸變大，分別在10 m高度處風速V10=15 m/s并且30 mm覆冰時達到位移最大值，地線最大值為2.2674 m，導(dǎo)線最大值為4.6723 m.在風速相同的情況下，輸電塔-線體系在覆冰厚度從10 mm增加到15 mm時，位移增加值最大，地線位移最高可增加0.944 0 m，導(dǎo)線位移最高可增加0.925 9 m.將在各工況下的數(shù)據(jù)進行擬合得出一次功能函數(shù)為

Dmax=-0.063 24+0.0124 2v+0.005 17H.

(1)

4 覆冰輸電塔抗風動力可靠度分析

結(jié)構(gòu)的破壞主要有剛度、強度、疲勞破壞，對于輸電塔這種高聳結(jié)構(gòu)而言，剛度破壞更適用于計算結(jié)構(gòu)可靠度，因此選擇輸電塔結(jié)構(gòu)的剛度失效為失效準則.根據(jù)我國輸電線路規(guī)范，在風速為5 m/s且無覆冰情況時，輸電塔的極限位移值為3 h/1 000，為了保證輸電塔在不同地區(qū)的穩(wěn)定性，假設(shè)輸電塔的極限位移為塔高的千分之一，即輸電塔架最高點的最大位移超過輸電塔高Ht的1/1 000時，認為輸電塔剛度失效.根據(jù)極限狀態(tài)方程，可得功能函數(shù)

Z=10-3Ht-Dmax.

(2)

將公式(1)代入公式(2)，整理可得

Z=0.108 24-1.24×10-3v-5.17×10-4H.

(3)

由公式(2)可以得到，功能函數(shù)Z的值隨著是風速v和覆冰厚度H的增加單調(diào)遞減，這是由于覆冰塔架的最大位移Dmax在脈動風荷載作用下，風速越大、覆冰厚度越厚，最大位移值越大，可靠指標越小.即當Z>0時，結(jié)構(gòu)正常；當Z=0時，結(jié)構(gòu)處于極限狀態(tài)；當Z<0 時，則認為結(jié)構(gòu)失效.基于響應(yīng)面法得到輸電塔的一次功能函數(shù)的表達式后，應(yīng)用Monte Carlo法計算結(jié)構(gòu)可靠度.可計算出輸電塔的可靠指標β為1.432 94.同時也分析了當覆冰或風速為定值，改變另一變量時的輸電塔結(jié)構(gòu)的可靠度.根據(jù)輸電塔-線耦合體系的風振響應(yīng)分析得到，體系最高點的最大位移隨著風速和覆冰厚度的增加而不斷增加，可靠度指標隨之減小.為了方便分析，以10 m高度處平均風速v=10 m/s，H=30 mm時為例，輸電塔的可靠度如表2、表3所示.

表2 覆冰厚度H=30 mm時不同風速下塔架的可靠度

表3 風速v=10 m/s時不同覆冰厚度下塔架的可靠度

分析結(jié)果表明，在風速逐漸變大與輸電塔-線耦合體系上的覆冰逐漸變厚的同時，輸電塔的可靠指標和可靠概率在逐漸變?。伙L速的改變對結(jié)構(gòu)的可靠指標影響較大，可靠度指標隨風速的增加急劇減?。桓脖穸葘Y(jié)構(gòu)的可靠指標影響較小，可靠度指標隨覆冰厚度的增大而慢慢減小.在計算二次功能函數(shù)中風速與覆冰厚度兩個變量的相關(guān)性時，得到當兩個變量都為定值時，輸電塔的可靠指標與相關(guān)系數(shù)成反比，并且根據(jù)公式可以得出當覆冰厚度與風速的相關(guān)系數(shù)為1時，輸電塔可靠指標達到最?。划斚嚓P(guān)系數(shù)為-1時，可靠指標最大，兩極值相差5.3%，證明了隨機變量的相關(guān)性會影響到結(jié)構(gòu)可靠度的計算結(jié)果.

5 結(jié) 論

(1)輸電塔風振響應(yīng)的最大位移發(fā)生在輸電塔最高處，導(dǎo)地線風振響應(yīng)的最大位移發(fā)生在弧垂最低點處并且最大位移隨著風速、覆冰厚度的增大而變大;

(2)風速的改變對結(jié)構(gòu)的可靠指標影響較大，并可靠指標隨風速的增加急劇減小；覆冰厚度對結(jié)構(gòu)的可靠指標影響較小，并可靠指標隨覆冰厚度的增大而慢慢減小;

(3)在計算輸電塔結(jié)構(gòu)可靠度時覆冰厚度和風速的變量相關(guān)性起到了重要作用，在計算冰風荷載作用下的輸電塔可靠度時，不能忽略風速與覆冰厚度兩個隨機變量的相關(guān)性.