楊 倫，樓文娟，潘小濤

浙江大學(xué)結(jié)構(gòu)工程研究所，杭州310058

舞動(dòng)是輸電導(dǎo)線在惡劣風(fēng)雨雪冰凍氣候條件下產(chǎn)生的大幅低頻振動(dòng)，會(huì)對(duì)輸電塔線造成構(gòu)件及電氣方面的損壞，危害電網(wǎng)的安全和電力的正常供應(yīng)［1-2］.

自Den Hartog［3］和Nigol 等［4］分別提出豎向和扭轉(zhuǎn)舞動(dòng)機(jī)理以來(lái)，不少學(xué)者針對(duì)覆冰導(dǎo)線提出了各種二自由度［5-6］和三自由度耦合舞動(dòng)質(zhì)點(diǎn)模型［7-8］，并采用攝動(dòng)法、多尺度法等近似解析的方式求解了舞動(dòng)響應(yīng). 然而鑒于質(zhì)點(diǎn)舞動(dòng)模型無(wú)法考慮導(dǎo)線大幅振動(dòng)所致的幾何非線性效應(yīng)，Luongo［9］、Liu［10］和Yan［11］等建立了二維或三維連續(xù)體耦合非線性模型，并針對(duì)導(dǎo)線在初始條件下的舞動(dòng)穩(wěn)定性、分岔特征以及內(nèi)共振模式展開(kāi)了一系列的定性分析. 值得注意的是，以上連續(xù)體舞動(dòng)模型仍然過(guò)于簡(jiǎn)化，忽略了實(shí)際大氣環(huán)境下線路不均勻覆冰、幾何剛度時(shí)變特征、不同階振型之間的耦合效應(yīng)以及不同檔距導(dǎo)線之間的相互作用等因素的影響，不能有效反映輸電線路舞動(dòng)響應(yīng)的動(dòng)態(tài)物理特征. 為此，Desai 等［12］針對(duì)覆冰導(dǎo)線舞動(dòng)的非線性特征，提出了三節(jié)點(diǎn)拋物線索單元，并建立了用于單導(dǎo)線舞動(dòng)分析的非線性有限元法. 與商業(yè)有限元軟件提供的索單元不同的是，Desai 提出的索單元具有扭轉(zhuǎn)自由度，同時(shí)能考慮偏心覆冰對(duì)舞動(dòng)的影響. 李黎等［13］利用具有扭轉(zhuǎn)自由度的二節(jié)點(diǎn)索單元，并通過(guò)將分裂導(dǎo)線等效為單導(dǎo)線的方式，建立了連續(xù)多檔距覆冰導(dǎo)線有限元模型. 然而有關(guān)研究表明，將分裂導(dǎo)線等效為單根導(dǎo)線無(wú)法考慮子導(dǎo)線之間的尾流干擾效應(yīng)，從而導(dǎo)致得出的舞動(dòng)幅值偏小［14］. 因此，劉小會(huì)等［15］基于Desai 提出的三節(jié)點(diǎn)拋物線索單元，采用對(duì)索節(jié)點(diǎn)扭轉(zhuǎn)自由度擴(kuò)充為3 個(gè)方向上轉(zhuǎn)動(dòng)自由度的方式，建立了能夠?qū)Ω髯訉?dǎo)線施加不同氣動(dòng)荷載，且考慮間隔棒與子導(dǎo)線運(yùn)動(dòng)耦合效應(yīng)的分裂導(dǎo)線有限元模型.

值得注意的是，覆冰分裂導(dǎo)線的扭轉(zhuǎn)剛度遠(yuǎn)大于單導(dǎo)線，致使兩者在起舞機(jī)理方面有著顯著差別. 為此，本研究基于完全拉格朗日格式(Total Lagrange)，建立適用于單導(dǎo)線和分裂導(dǎo)線舞動(dòng)數(shù)值模擬的非線性有限元?jiǎng)恿Ψ治龇椒ǎ捎镁哂信まD(zhuǎn)自由度的三節(jié)點(diǎn)拋物線索單元離散覆冰單導(dǎo)線. 對(duì)于覆冰分裂導(dǎo)線，在單導(dǎo)線有限元法的基礎(chǔ)上，利用歐拉梁?jiǎn)卧M間隔棒的運(yùn)動(dòng)過(guò)程. 為盡可能地提高計(jì)算效率，提出梁?jiǎn)卧D(zhuǎn)動(dòng)自由度縮聚法實(shí)現(xiàn)間隔棒與分裂子導(dǎo)線之間的耦合，并運(yùn)用隨轉(zhuǎn)坐標(biāo)系法求解舞動(dòng)過(guò)程中的梁節(jié)點(diǎn)不平衡力. 在此基礎(chǔ)上，結(jié)合覆冰導(dǎo)線氣動(dòng)力系數(shù)的風(fēng)洞試驗(yàn)結(jié)果，分別考察覆冰單導(dǎo)線和四分裂導(dǎo)線在湍流和均勻流場(chǎng)中的起舞機(jī)理和舞動(dòng)響應(yīng)特征之異同.

1 覆冰導(dǎo)線有限元法

1.1 單導(dǎo)線的動(dòng)力平衡方程和氣動(dòng)荷載

基于完全拉格朗日格式，具有扭轉(zhuǎn)自由度的三節(jié)點(diǎn)拋物線索單元描述導(dǎo)線單元的非線性動(dòng)力平衡方程［12］為

如圖1 所示，作用于覆冰導(dǎo)線某截面的氣動(dòng)荷載可表示為

其中，F(xiàn)x、Fy和M 分別為x 向、y 向和扭轉(zhuǎn)向的節(jié)點(diǎn)氣動(dòng)荷載;ρa(bǔ)ir為空氣密度;D 為導(dǎo)線直徑;Uz為來(lái)流風(fēng)速;CL、CD和CM分別是導(dǎo)線截面的升力系數(shù)、阻力系數(shù)和扭轉(zhuǎn)系數(shù)，與導(dǎo)線覆冰形狀和風(fēng)對(duì)導(dǎo)線的攻角α 有關(guān)，

其中，˙y 為豎向運(yùn)動(dòng)速度;β 和θ 分別為初始風(fēng)攻角和t 時(shí)刻導(dǎo)線的扭轉(zhuǎn)角. 式(6)右第3 和第4 項(xiàng)分別代表導(dǎo)線豎向和扭轉(zhuǎn)運(yùn)動(dòng)速度對(duì)總風(fēng)攻角的影響.

圖1 覆冰導(dǎo)線截面風(fēng)攻角及受力Fig.1 Attack angle and aerodynamic force of iced conductor

1.2 覆冰分裂導(dǎo)線有限元模型

分裂導(dǎo)線由多根單導(dǎo)線和間隔棒共同組成. 間隔棒的作用在于保持子導(dǎo)線間距，防止子導(dǎo)線之間由于電磁吸引以及風(fēng)力而引發(fā)的相互靠近和碰撞鞭擊. 同時(shí)，受間隔棒約束作用的影響，分裂導(dǎo)線舞動(dòng)時(shí)表現(xiàn)為顯著的整體運(yùn)動(dòng). 因此從數(shù)值模擬的角度來(lái)看，只需對(duì)按照一定規(guī)律排列的多根子導(dǎo)線構(gòu)成的振動(dòng)系統(tǒng)中加入模擬間隔棒的梁?jiǎn)卧纯? 這一過(guò)程中首先需要解決的是如何高效、可靠地實(shí)現(xiàn)分裂子導(dǎo)線與間隔棒的連接. 另外，分裂導(dǎo)線發(fā)生舞動(dòng)時(shí)間隔棒的運(yùn)動(dòng)表現(xiàn)出典型的大轉(zhuǎn)動(dòng)、小應(yīng)變特征，精確求解間隔棒在運(yùn)動(dòng)過(guò)程中的單元不平衡抗力向量，是保證舞動(dòng)計(jì)算收斂性的關(guān)鍵. 本研究分別采用梁?jiǎn)卧D(zhuǎn)動(dòng)自由度縮聚和隨轉(zhuǎn)坐標(biāo)系法［16］，實(shí)現(xiàn)子導(dǎo)線與間隔棒的連接，并求解梁?jiǎn)卧?jié)點(diǎn)的不平衡抗力向量.

1.2.1 分裂子導(dǎo)線與間隔棒的連接

分裂導(dǎo)線由單導(dǎo)線和間隔棒組成，其有限元模型如圖2 所示. 各子導(dǎo)線之間通過(guò)由歐拉梁?jiǎn)卧M的間隔棒相連，其中梁?jiǎn)卧?jié)點(diǎn)有3 個(gè)平動(dòng)自由度及3 個(gè)扭轉(zhuǎn)自由度，子導(dǎo)線單元節(jié)點(diǎn)有3 個(gè)平動(dòng)自由度和1 個(gè)轉(zhuǎn)動(dòng)自由度.

圖2 分裂導(dǎo)線有限元模型Fig.2 Finite element model of bundled conductors

針對(duì)分裂子導(dǎo)線與間隔棒的連接問(wèn)題，文獻(xiàn)［15］提出一種對(duì)索單元節(jié)點(diǎn)扭轉(zhuǎn)自由度擴(kuò)張的方式來(lái)模擬子導(dǎo)線與間隔棒的連接. 即在間隔棒與子導(dǎo)線的交點(diǎn)處，將索節(jié)點(diǎn)扭轉(zhuǎn)自由度投影至梁節(jié)點(diǎn)3個(gè)方向的轉(zhuǎn)動(dòng)自由度上，交點(diǎn)處索節(jié)點(diǎn)自由度數(shù)由4 個(gè)增加至6 個(gè). 采用該連接方式會(huì)在一定程度上增加整個(gè)系統(tǒng)的自由度數(shù)，從而影響求解效率.

為盡可能地減小計(jì)算量，提高舞動(dòng)分析效率，本研究采用對(duì)梁節(jié)點(diǎn)扭轉(zhuǎn)自由度縮聚的方法實(shí)現(xiàn)梁?jiǎn)卧?jié)點(diǎn)與索單元節(jié)點(diǎn)的耦合. 在間隔棒與子導(dǎo)線的連接處，保持索單元節(jié)點(diǎn)的4 個(gè)自由度不變，將梁?jiǎn)卧?jié)點(diǎn)在整體坐標(biāo)系下的3 個(gè)轉(zhuǎn)動(dòng)自由度縮聚為1 個(gè)沿索單元軸線方向的扭轉(zhuǎn)自由度. 以導(dǎo)線初始構(gòu)型為參考構(gòu)型，t 時(shí)刻梁?jiǎn)卧谡w坐標(biāo)系下的節(jié)點(diǎn)位移向量可表示為

圖3 分裂子導(dǎo)線與間隔棒的連接Fig.3 Connection between sub-conductor and space rod

1.2.2 梁?jiǎn)卧?jié)點(diǎn)抗力的計(jì)算

分裂導(dǎo)線舞動(dòng)過(guò)程中，間隔棒的節(jié)點(diǎn)位移可分解為剛體運(yùn)動(dòng)和單元變形兩部分. 傳統(tǒng)的完全拉格朗日列式法僅適用于轉(zhuǎn)動(dòng)較小的情形，轉(zhuǎn)動(dòng)位移較大時(shí)無(wú)法有效分離單元的剛體轉(zhuǎn)動(dòng)和單元變形，導(dǎo)致在結(jié)構(gòu)非線性增量求解過(guò)程中無(wú)法精確計(jì)算由于純變形產(chǎn)生的單元抗力. 為此，本研究運(yùn)用隨轉(zhuǎn)坐標(biāo)系方法［16］解決問(wèn)題. 如圖4，選定t = 0 時(shí)刻梁?jiǎn)卧獦?gòu)型所在的局部坐標(biāo)系為初始坐標(biāo)系，則t 時(shí)刻梁?jiǎn)卧?jié)點(diǎn)在初始坐標(biāo)系下的位移向量可表示為

圖4 梁?jiǎn)卧煌瑫r(shí)刻的構(gòu)型Fig.4 Configurations of beam element at different time

在初始局部坐標(biāo)系下梁?jiǎn)卧獌啥似絼?dòng)位移增量可表示為

由于梁?jiǎn)卧谶\(yùn)動(dòng)過(guò)程中變形較小，因此在t 時(shí)刻的弦長(zhǎng)Lt可表示為

則梁?jiǎn)卧纳扉L(zhǎng)量ΔLt為

如圖5，將梁?jiǎn)卧趖 時(shí)刻局部坐標(biāo)系的坐標(biāo)原點(diǎn)移至初始時(shí)刻局部坐標(biāo)系的坐標(biāo)原點(diǎn)，則轉(zhuǎn)動(dòng)剛體位移可用梁端位移平動(dòng)位移增量和扭轉(zhuǎn)角表示為

梁節(jié)點(diǎn)在t 時(shí)刻局部坐標(biāo)系下的轉(zhuǎn)動(dòng)變形位移為

將梁?jiǎn)卧诰植孔鴺?biāo)系下的單元抗力轉(zhuǎn)換至整體坐標(biāo)系下，并進(jìn)行自由度縮聚，可得梁?jiǎn)卧谡w坐標(biāo)系下的節(jié)點(diǎn)抗力為

圖5 梁?jiǎn)卧D(zhuǎn)動(dòng)剛體位移Fig.5 Rigid rotational displacement of beam element

1.3 非線性動(dòng)力方程的求解

求得表征子導(dǎo)線和間隔棒舞動(dòng)特性的單元矩陣和荷載向量后，便可根據(jù)單元定位向量獲得單導(dǎo)線或者分裂導(dǎo)線系統(tǒng)的非線性運(yùn)動(dòng)方程. 本研究采用無(wú)條件穩(wěn)定的Newmark 法對(duì)方程直接積分求解，并運(yùn)用Newton-Raphson 法對(duì)每個(gè)時(shí)間步末尾的位移向量進(jìn)行迭代求解. Newmark 法采用平均加速度方案，即積分精度參數(shù)α 和穩(wěn)定性參數(shù)β 分別取0.25和0.5. 那么t 至t +Δt 過(guò)程中，基于完全拉格朗日格式的遞推迭代公式可表示為［18］

2 有限元算例驗(yàn)證

為驗(yàn)證本研究計(jì)算分裂導(dǎo)線方法的有效性，以某工程的單跨四分裂導(dǎo)線為例，采用本研究方法和ANSYS 有限元軟件求解導(dǎo)線舞動(dòng)響應(yīng). 分裂導(dǎo)線等間距設(shè)置5 組間隔棒，子導(dǎo)線物理參數(shù)如表1.

表1 覆冰導(dǎo)線物理參數(shù)Table 1 Physical parameters of iced conductor

ANSYS 模型中分別采用Link 10 單元和Beam 4單元離散子導(dǎo)線和間隔棒. 由于分裂導(dǎo)線的整體抗扭剛度比子導(dǎo)線繞自身軸的抗扭剛度大得多，因此在ANSYS 模型中不考慮子導(dǎo)線的抗扭剛度. 本研究計(jì)算模型中，每根子導(dǎo)線劃分20 個(gè)三節(jié)點(diǎn)拋物線索單元. 鑒于Link 10 單元是2 節(jié)點(diǎn)單元，為保證兩種方法的節(jié)點(diǎn)數(shù)量保持一致，ANSYS 計(jì)算模型中每根子導(dǎo)線劃分40 個(gè)單元.

分別在每根子導(dǎo)線跨中節(jié)點(diǎn)上同時(shí)施加400 N、400 N 和40 N·m 的垂直、水平和扭轉(zhuǎn)向突加荷載.經(jīng)計(jì)算，兩種方法所得的子導(dǎo)線跨中位移時(shí)間歷程如圖6. 不難看出，兩種方法得出的計(jì)算結(jié)果吻合很好，說(shuō)明本算法可靠有效.

圖6 子導(dǎo)線跨中位移時(shí)程Fig.6 Time histories of displacement of sub-conductor at the mid-span

3 覆冰導(dǎo)線舞動(dòng)數(shù)值模擬

3.1 覆冰導(dǎo)線氣動(dòng)力

以新月形覆冰為例，制作1∶1 四分裂導(dǎo)線節(jié)段模型，模型長(zhǎng)度為1.0 m. 鑒于舞動(dòng)發(fā)生時(shí)的風(fēng)速大多低于20 m/s，本研究試驗(yàn)風(fēng)速取為15 m/s. 利用高頻動(dòng)態(tài)測(cè)力天平，考察導(dǎo)線在均勻流和6%均勻湍流場(chǎng)下新月形覆冰四分裂導(dǎo)線整體氣動(dòng)三分力系數(shù)隨攻角的變化規(guī)律(如圖7).

圖7 覆冰導(dǎo)線氣動(dòng)三分力系數(shù)隨風(fēng)攻角的變化規(guī)律Fig.7 Aerodynamic coefficients of iced conductor versus wind attack angle

由于Den Hartog 系數(shù)和Nigol 系數(shù)分別體現(xiàn)了覆冰導(dǎo)線的豎向和扭轉(zhuǎn)穩(wěn)定性，因此結(jié)合氣動(dòng)力測(cè)試結(jié)果，給出了以上2 類系數(shù)隨初始攻角的變化規(guī)律(如圖8). 不難看出，初始攻角落在25°和175°附近時(shí)兩類系數(shù)均小于零，說(shuō)明在此情形下導(dǎo)線極有可能喪失穩(wěn)定性.

3.2 覆冰導(dǎo)線舞動(dòng)響應(yīng)分析

以表1 導(dǎo)線參數(shù)為例，采用本研究提出的有限元法，分別對(duì)單根和四分裂導(dǎo)線在均勻流和湍流條件下進(jìn)行舞動(dòng)分析. 結(jié)合圖8 給出的不穩(wěn)定攻角范圍，選取25°和175°為初始風(fēng)攻角. 鑒于舞動(dòng)屬于典型的自激振動(dòng)，主要與平均風(fēng)荷載有關(guān)，計(jì)算風(fēng)速取為13 m/s. 另外，為保證算法的收斂性，舞動(dòng)計(jì)算的時(shí)間步長(zhǎng)取為0.01 s.

3.2.1 覆冰單導(dǎo)線舞動(dòng)響應(yīng)分析

圖9 為25°攻角下覆冰單導(dǎo)線在2 類流場(chǎng)下的跨中舞動(dòng)響應(yīng). 初始風(fēng)攻角為25°時(shí)，無(wú)論在湍流還是均勻流作用下，單導(dǎo)線在垂直方向上均保持穩(wěn)定. 但在這2 種流場(chǎng)中未發(fā)生舞動(dòng)的原因有著本質(zhì)區(qū)別:在湍流場(chǎng)中，雖然Den Hartog 系數(shù)和Nigol系數(shù)均為負(fù)，但是單導(dǎo)線的抗扭剛度很小，在扭轉(zhuǎn)系不穩(wěn)定攻角范圍內(nèi)，單導(dǎo)線會(huì)發(fā)生振幅較大的扭轉(zhuǎn)舞動(dòng)，使動(dòng)態(tài)攻角極易脫離Den Hartog 系數(shù)小于零的區(qū)域，所以其豎向振動(dòng)幅值很小，可認(rèn)為保持穩(wěn)定;而流場(chǎng)為均勻流時(shí)，Den Hartog 系數(shù)和Nigol系數(shù)均為正，因此覆冰單導(dǎo)線并未發(fā)生舞動(dòng).

圖8 覆冰導(dǎo)線的Den Hartog 系數(shù)和Nigol 系數(shù)Fig.8 Nigol coefficient of iced conductor

圖9 25°風(fēng)攻角下單導(dǎo)線在兩類流場(chǎng)中的跨中舞動(dòng)響應(yīng)Fig.9 Galloping responses of mid-point of single conductor with attack angle of 25° in two types of flow field

圖10 給出當(dāng)風(fēng)向改變、使初始風(fēng)攻角變?yōu)?75°時(shí)的導(dǎo)線舞動(dòng)響應(yīng). 可見(jiàn)，湍流作用下覆冰單導(dǎo)線僅發(fā)生了大振幅的扭轉(zhuǎn)舞動(dòng)，其機(jī)理與湍流場(chǎng)中25°初始風(fēng)攻角下的舞動(dòng)機(jī)制相同. 當(dāng)流場(chǎng)為均勻流時(shí)，Nigol 系數(shù)轉(zhuǎn)變?yōu)檎鳧en Hartog 系數(shù)為負(fù)，因此單導(dǎo)線發(fā)生了振幅較大的豎向舞動(dòng). 從圖9 (b)和圖10 (b)均不難發(fā)現(xiàn)，單導(dǎo)線扭轉(zhuǎn)舞動(dòng)的特征頻率更高.

圖10 175°風(fēng)攻角下單導(dǎo)線在兩類流場(chǎng)中的跨中舞動(dòng)響應(yīng)Fig.10 Galloping responses of mid-point of single conductor with attack angle of 175° in two types of flow field

3.2.2 覆冰四分裂導(dǎo)線舞動(dòng)響應(yīng)分析

圖11 為25°初始攻角下四分裂導(dǎo)線在兩類流場(chǎng)中的跨中舞動(dòng)響應(yīng). 由于分裂導(dǎo)線的扭轉(zhuǎn)剛度較大，因此在湍流場(chǎng)中即使初始攻角落入Nigol 系數(shù)的不穩(wěn)定區(qū)域，其扭轉(zhuǎn)向僅發(fā)生了小幅振動(dòng). 在此情形下，Den Hartog 系數(shù)起主導(dǎo)作用，所以分裂導(dǎo)線發(fā)生了顯著的豎向舞動(dòng). 流場(chǎng)為均勻流時(shí)，Den Hartog 系數(shù)和Nigol 系數(shù)均大于零，因此分裂導(dǎo)線豎向和扭轉(zhuǎn)向均保持穩(wěn)定.

圖11 25°風(fēng)攻角下分裂導(dǎo)線子導(dǎo)線在兩類流場(chǎng)中的跨中舞動(dòng)響應(yīng)Fig.11 Galloping responses of mid-point of sub-conductor with attack angle of 25° in two types of flow field

圖12 175°風(fēng)攻角下分裂導(dǎo)線子導(dǎo)線在兩類流場(chǎng)中的跨中舞動(dòng)響應(yīng)Fig.12 Galloping responses of mid-point of sub-conductor with attack angle of 175° in two types of flow field

圖12 給出175°攻角下分裂導(dǎo)線的舞動(dòng)響應(yīng).結(jié)合圖8 不難看出，在湍流場(chǎng)中，雖然Nigol 系數(shù)為負(fù)，但由于分裂導(dǎo)線的抗扭剛度較大且Den Hartog 系數(shù)大于零，導(dǎo)線并未舞動(dòng).然而在均勻流場(chǎng)中，Den Hartog 系數(shù)變?yōu)樨?fù)值，因此分裂導(dǎo)線發(fā)生了振幅較大的豎向舞動(dòng). 同時(shí)，受豎向舞動(dòng)的激勵(lì)作用，分裂導(dǎo)線在扭轉(zhuǎn)向發(fā)生了振幅較小的受迫振動(dòng).

結(jié) 語(yǔ)

基于完全拉格朗日格式，采用具有扭轉(zhuǎn)自由度的三節(jié)點(diǎn)拋物線索單元和歐拉梁?jiǎn)卧謩e模擬導(dǎo)線和間隔棒，提出了梁節(jié)點(diǎn)自由度縮聚法，實(shí)現(xiàn)了子導(dǎo)線和間隔棒的連接，建立了適用于單導(dǎo)線和分裂導(dǎo)線舞動(dòng)分析的非線性有限元法. 結(jié)合覆冰導(dǎo)線在湍流場(chǎng)和均勻流場(chǎng)中的氣動(dòng)力測(cè)試結(jié)果，求解了新月形覆冰單導(dǎo)線和四分裂導(dǎo)線的舞動(dòng)響應(yīng)，并揭示了單導(dǎo)線和分裂導(dǎo)線在不同流場(chǎng)中的起舞機(jī)理.研究表明:由于單導(dǎo)線抗扭剛度較小，在Nigol 系數(shù)為負(fù)的情況下會(huì)發(fā)生大幅扭轉(zhuǎn)舞動(dòng)，所以動(dòng)態(tài)攻角極易脫離升力系數(shù)的不穩(wěn)定區(qū)域，保證了單導(dǎo)線的豎向穩(wěn)定性;對(duì)于抗扭剛度較大的分裂導(dǎo)線來(lái)說(shuō)，其氣動(dòng)穩(wěn)定性主要受升力系數(shù)的特征控制，即使在Den Hartog 系數(shù)和Nigol 系數(shù)均為負(fù)的情況下，仍易發(fā)生大幅的豎向舞動(dòng). 下一步筆者將基于D 形覆冰四、六、八分裂導(dǎo)線子導(dǎo)線氣動(dòng)力測(cè)試結(jié)果，對(duì)各子導(dǎo)線施加不同的氣動(dòng)荷載，運(yùn)用本文給出的非線性有限元算法，細(xì)致考察子導(dǎo)線間氣流干擾效應(yīng)對(duì)分裂導(dǎo)線舞動(dòng)特征的影響.

/ References:

［1］Zhu Kuanjun，F(xiàn)u Dongjie，Wang Jingchao，et al. The galloping and its preventing techniques on overhead transmission line ［J］. Electrical Equipment，2009，9(6):8-12.(in Chinese)朱寬軍，付東杰，王景朝，等. 架空輸電線路的舞動(dòng)及防治［J］. 電力設(shè)備，2009，9(6):8-12.

［2］Wei Chong，Pan Shaocheng，Tang Minggui，et al. Analysis and threatment of the galloping of 500kV transmission line ［J］. Electric Power Construction，2011，32(4):24-27.(in Chinese)魏 沖，潘少成，唐明貴，等. 500kV 輸電線路舞動(dòng)分析及治理［J］. 電力建設(shè)，2011，32(4):24-27.

［3］Den Hartog J P. Transmission line vibration due to sleet［J］. AIEE Transactions，1932，51(4):1074-1086.

［4］Nigol O，Clarke G J，Havard D G. Torsional stability of bundle conductors ［J］. IEEE Transactions on Power Apparatus and Systems，1977，96(5):1666-1674.

［5］Yu P，Popplewell N，Shah A H. Instability trends of inertially coupled galloping ［J］. Journal of Sound and Vibration，1995，183(4):663-678.

［6］Macdonald J H G，Larose G L. Two-degree-of-freedom inclined cable galloping ［J］. Journal of Wind Engineering and Industrial Aerodynamics，2008，96(3):291-326.

［7］Wang J W. Large vibrations of overhead electrical lines:a full 3-dof model for galloping studies ［D］. Province of Liège (Belgum):Universite de Liege，1996，151:1-227.

［8］Chen Xiaoming，Deng Hongzhou，Wang Zhaomin. Conductor galloping stability analysis of long-span transmission system ［J］. Engineering Mechanics，2004，21(1):56-60.(in Chinese)陳曉明，鄧洪洲，王肇民. 大跨越輸電線路舞動(dòng)穩(wěn)定性研究［J］. 工程力學(xué)，2004，21(1):56-60.

［9］Luongo A，Piccardo G. Non-linear galloping of sagged cables in 1:2 internal resonance ［J］. Journal of Sound and Vibration，1998，214(5):915-940.

［10］Liu Haiying，Zhang Qichang，Hao Shuying. A study on galloping for iced quad-bundled conductor ［J］. Journal of Vibration Engineering，2011，24(3):235-239.(in Chinese)劉海英，張琪昌，郝淑英. 覆冰四分裂導(dǎo)線舞動(dòng)研究［J］. 振動(dòng)工程學(xué)報(bào)，2011，24(3):235-239.

［11］Yan Z M，Yan Z T，Li Z，et al. Nonlinear galloping of internally resonant iced transmission lines considering eccentricity ［J］. Journal of Sound and Vibration，2012，331:3599-3616.

［12］Desai Y M，Yu P，Popplewell N，et al. Finite element modelling of transmission line galloping ［J］. Computers＆ Structures，1995，57(3):407-420.

［13］Li Li，Chen Yuankun，Xia Zhengchun，et al. Nonlinear numerical simulation study of iced conductor galloping［J］. Journal of Vibration and Shock，2011，30(8):107-111.(in Chinese)李 黎，陳元坤，夏正春，等. 覆冰導(dǎo)線舞動(dòng)的非線性數(shù)值仿真研究［J］. 振動(dòng)與沖擊，2011，30(8):107-111.

［14］Yan Bo，Li Wenyun，Zhou Song，et al. Numerical simulation on galloping of quad-bundled conductor ［J］. Journal of Vibration and Shock，2010，29(9):102-107.(in Chinese)嚴(yán) 波，李文蘊(yùn)，周 松，等. 覆冰四分裂導(dǎo)線舞動(dòng)數(shù)值模擬研究［J］. 振動(dòng)與沖擊，2010，29(9):102-107.

［15］Liu Xiaohui，Yan Bo，Zhang Hongyan，et al. Nonlinear finite element analysis for galloping of iced bundled conductors［J］. Journal of Vibration and Shock，2010，29(6):129-133.(in Chinese)劉小會(huì)，嚴(yán) 波，張宏雁，等. 分裂導(dǎo)線舞動(dòng)非線性有限元分析方法［J］. 振動(dòng)與沖擊，2010，29(6):129-133.

［16］Ma Yuquan. A new CR formulation and its application in geometrically nonlinear analysis ［J］. Science Technology and Engineering，2010，20(10):7476-7488. (in Chinese)馬玉全. 改進(jìn)CR 列式法及在幾何非線性分析中的應(yīng)用［J］. 科學(xué)技術(shù)與工程，2010，20(10):7476-7488.

［17］Yang Qi，Huang Jianyue. 3D geometric nonlinear simulation analysis of long span bridges ［J］. Central South Highway Engineering，2005，30(1):59-62. (in Chinese)楊 琪，黃建躍. 大跨度橋梁空間幾何非線性仿真分析的研究［J］. 中南公路工程，2005，30(1):59-62.

［18］Wang Xucheng，Shao Min. Basic Principles and Numerical Simulation of Finite Element Method ［M］. Beijing:Tsinghua University Press，1995:549-550.(in Chinese)王勖成，邵 敏. 有限單元法基本原理和數(shù)值方法［M］. 北京:清華大學(xué)出版社，1995:549-550.