索網(wǎng)結(jié)構(gòu)中的空間兩節(jié)點(diǎn)懸鏈線索單元

劉玥君，湯愛平，劉克同

(1.東北電力大學(xué) 建筑工程學(xué)院，吉林 吉林132012； 2.哈爾濱工業(yè)大學(xué) 土木工程學(xué)院，黑龍江 哈爾濱 150090)

?

索網(wǎng)結(jié)構(gòu)中的空間兩節(jié)點(diǎn)懸鏈線索單元

劉玥君1，湯愛平2，劉克同2

(1.東北電力大學(xué) 建筑工程學(xué)院，吉林 吉林132012； 2.哈爾濱工業(yè)大學(xué) 土木工程學(xué)院，黑龍江 哈爾濱 150090)

摘要：為了分析索結(jié)構(gòu)在荷載作用下的線性和非線性反應(yīng)，根據(jù)索結(jié)構(gòu)的受力特點(diǎn)，提出了一種可以充分考慮索結(jié)構(gòu)非線性效應(yīng)的空間兩節(jié)點(diǎn)懸鏈線索單元。基于彈性懸鏈線的解析解，推導(dǎo)了索單元的切線剛度矩陣,通過迭代法求解了索結(jié)構(gòu)的剛度矩陣和內(nèi)力向量，利用該新型的索單元模型分析了靜、動(dòng)力荷載作用下結(jié)構(gòu)的非線性行為,對于非線性靜力分析，結(jié)構(gòu)平衡方程采用牛頓迭代法求解；對于非線性動(dòng)力時(shí)程分析，采用Newmark直接積分法和牛頓迭代法相結(jié)合的增量迭代法求解。數(shù)值分析結(jié)果驗(yàn)證了空間兩節(jié)點(diǎn)懸鏈線索單元模擬索結(jié)構(gòu)的精確性和可靠性。

關(guān)鍵詞：空間兩節(jié)點(diǎn)；懸鏈線索單元；靜力荷載；動(dòng)力荷載；非線性分析

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20160421.1040.030.html

近幾十年來，索結(jié)構(gòu)廣泛的應(yīng)用于輸電線、斜拉橋等大跨越空間結(jié)構(gòu)。索結(jié)構(gòu)具有很強(qiáng)的非線性特性，索結(jié)構(gòu)受力分析需要考慮擾度和大變形的影響。在索結(jié)構(gòu)的研究方面，陳太聰?shù)萚1]總結(jié)和分析了不同工況下索結(jié)構(gòu)中拉索的受力特點(diǎn)，將拉索行為歸納為已知原長和已知索力兩種典型受力情況。利用懸鏈線拉索解析解的線性展開，引入適當(dāng)?shù)募s束條件，建立兩種新型的拉索單元，提出索結(jié)構(gòu)非線性靜力分析的一般迭代算法，利用等效線自重考慮彈性變形對拉索靜力狀態(tài)的影響。魏建東等[2]基于有限元分析的基本原理，創(chuàng)建了一種新的滑移索單元，分析索結(jié)構(gòu)中索段在鞍座上的滑動(dòng)問題。唐建民等[3]根據(jù)張拉結(jié)構(gòu)的特點(diǎn)，提出了一種五節(jié)點(diǎn)等參數(shù)單元有限元模型，采用四節(jié)點(diǎn)多項(xiàng)式作為位移差值函數(shù)及單元初始形狀函數(shù)，基于修正的Lagrangian坐標(biāo)描述法，建立了非線性有限元基本方程和切線剛度矩陣，利用Newton-Raphson法進(jìn)行了實(shí)例計(jì)算。黃濤等[4]提出了具有扭轉(zhuǎn)自由度的五節(jié)點(diǎn)索單元非線性有限元法，考慮了索結(jié)構(gòu)沿軸向的轉(zhuǎn)動(dòng)自由度，這種方法不僅適用于一般的懸索結(jié)構(gòu)的分析，也適用于需要考慮扭轉(zhuǎn)的特殊懸索結(jié)構(gòu)的分析。

現(xiàn)有的索結(jié)構(gòu)數(shù)值分析技術(shù)主要采用下面兩種方法:1)基于多項(xiàng)式插值函數(shù)的有限元方法，2)基于彈性懸鏈線解析表達(dá)式的分析方法。在第1種方法中，采用插值函數(shù)來描述索結(jié)構(gòu)的非線性，這種方法常用來模擬兩節(jié)點(diǎn)和多節(jié)點(diǎn)單元及具有轉(zhuǎn)動(dòng)自由度的曲線單元，但是該方法只適合模擬具有較大預(yù)應(yīng)力的索結(jié)構(gòu)。對于長度較小，承受拉力較高的索，采用等效彈性模量來計(jì)算是足夠準(zhǔn)確的，但對于弧垂比較大的索，該方法的適用性較差[5]。對大弧垂索結(jié)構(gòu)，一般采用一系列的直線單元來模擬幾何彎曲的索結(jié)構(gòu)，多節(jié)點(diǎn)索單元是基于高階多項(xiàng)式的插值函數(shù)[6-8]來構(gòu)造，因此計(jì)算結(jié)果的精度較兩節(jié)點(diǎn)索單元有一定的提高，但多節(jié)點(diǎn)單元會(huì)破壞索結(jié)構(gòu)曲線斜率的連續(xù)性。除此之外，多節(jié)點(diǎn)單元需要使用大量的單元來建立幾何彎曲的索結(jié)構(gòu)模型，對于多索結(jié)構(gòu)體系，這種方法的計(jì)算成本巨大。第二種方法是由Jayaraman 等[9]提出的采用彈性懸鏈線的精確表達(dá)式來描述索結(jié)構(gòu)的真實(shí)力學(xué)行為，經(jīng)過Wang[10]、 Andreu[11]、 Yang[12]、Such[13]、Andreu[14]等進(jìn)行改進(jìn)。用懸鏈線單元對曲線索結(jié)構(gòu)進(jìn)行建模，與多項(xiàng)式插值函數(shù)有限元方法相比有很多優(yōu)點(diǎn)，如計(jì)算需要的自由度數(shù)少，可以完全考慮索結(jié)構(gòu)的非線性效應(yīng)。

基于彈性懸鏈線解析理論，本文開發(fā)了一種空間兩節(jié)點(diǎn)懸鏈線索單元，用來對索結(jié)構(gòu)承受的靜態(tài)和動(dòng)態(tài)荷載進(jìn)行非線性分析。考慮索結(jié)構(gòu)的預(yù)應(yīng)力效應(yīng)，推導(dǎo)并給出索單元的切線剛度矩陣和內(nèi)部力向量。采用構(gòu)造的索單元對弧垂較小的索網(wǎng)結(jié)構(gòu)和弧垂較大的輸電線索結(jié)構(gòu)進(jìn)行計(jì)算分析，并對計(jì)算結(jié)果進(jìn)行了驗(yàn)證。

1懸鏈線索單元

為了能夠模擬索結(jié)構(gòu)的真實(shí)形態(tài)，根據(jù)彈性懸鏈線單元的解析表達(dá)式，本文推導(dǎo)出來一種新的懸鏈線索單元。假設(shè)索結(jié)構(gòu)是完全柔性的，自重沿索長度均勻分布，索的橫截面積保持不變。圖1顯示

了索結(jié)構(gòu)兩懸掛點(diǎn)I、J的笛卡爾坐標(biāo)分別為(0，0，0)和(lx，ly，lz)，s和p為變形前后索上任意一點(diǎn)的曲線坐標(biāo)。索在平衡時(shí)的方程可以表示為

(1)式中：F1、F2和F3分別是索結(jié)構(gòu)在x、y和z軸的張力分量，w是索結(jié)構(gòu)的自重，索結(jié)構(gòu)的張力在拉格朗日坐標(biāo)下的表達(dá)式為

(3)式中：E是索的彈性模量，A是索的橫截面積。曲線坐標(biāo)s與笛卡爾坐標(biāo)之間的關(guān)系可以表示為

(4)

圖1　三維懸鏈線索單元Fig.1　Three-dimensional catenary cable element

(5)將方程(1)～(3)代入方程(4)，引入邊界條件(5)，索的長度在x、y和z 3個(gè)方向的分量為

(6)

式中：L0是索的原長，根據(jù)端部節(jié)點(diǎn)力(F1，F(xiàn)2，F(xiàn)3)，lx、ly和lz可以表示成如下形式：

(7)

切線剛度矩陣和相應(yīng)單元的內(nèi)力向量可以通過求解方程(6)來推導(dǎo)，對方程(6)的兩邊同時(shí)進(jìn)行微分，可以得到下面的方程：

(8)

或表示成矩陣的形式：

(9)

F是柔度矩陣，具體可以表示為

(10)

式中：Ti和Tj分別是節(jié)點(diǎn)I、J的張力，可以定義為

(11)

節(jié)點(diǎn)J的節(jié)點(diǎn)力(F4，F(xiàn)5，F(xiàn)6)可以通過平衡方程獲得

(12)

剛度矩陣通過對求柔度矩陣F的逆矩陣而得到

(13)

索單元的切線剛度矩陣和內(nèi)力向量可以表示成六自由度的形式：

(14)

(15)

如果索單元的切線剛度矩陣和內(nèi)力向量已經(jīng)確定，索長和索的弧垂如圖2可以獲得[9,14]：

(16)

(17)

式中：

(18)

圖2　斜拉索的弧垂ZsFig.2　The sag of stay cable Zs

采用彈塑性鉸模型來表示索單元的非彈性行為，假設(shè)索單元的塑性分布是集中在索的兩端，而整個(gè)索單元仍然保持彈性。如果索的軸向力大于屈服力Py=Aσy，索單元的彈性模量等于零，軸向力等于屈服力Py[15]。

2剛度矩陣及求解方法

應(yīng)用迭代算法來計(jì)算索單元的切線剛度矩陣和內(nèi)力向量。這個(gè)程序的計(jì)算需要端點(diǎn)力的初值(F1，F(xiàn)2，F(xiàn)3)，根據(jù)已知的懸鏈線方程，端點(diǎn)力的初值可以通過以下的方式獲得[9]

(19)

式中：

(20)

計(jì)算切線剛度矩陣和內(nèi)力矢量的迭代程序如下：

1)輸入w、E、A、L0，節(jié)點(diǎn)I(xi,yi,zi)和節(jié)點(diǎn)J(xj,yj,zj)。

2)計(jì)算lx0=xj-xi，ly0=yj-yi，lz0=zj-zi。

4)使用式(6)更新(lx，ly，lz)。

5)計(jì)算閉合差向量

6)如果dL小于預(yù)設(shè)的允許值，應(yīng)用式(14)計(jì)算KT，用式(15)計(jì)算Fint，如果dL大于預(yù)設(shè)的允許值，程序繼續(xù)進(jìn)行下一步。

7)計(jì)算端點(diǎn)力的修正向量dF=F-1dL。

8)更新端點(diǎn)力Fi+1=Fi+dF，然后回到4)。

如果給定索單元的張力初值為T0，索原長未知，可以采用一個(gè)類似的迭代程序來計(jì)算無應(yīng)力索的原長，具體迭代程序如下：

1)輸入w，E，A，節(jié)點(diǎn)I(xi，yi，zi)和節(jié)點(diǎn)J(xi，yi，zi)。

2)計(jì)算lx0=xj-xi，ly0=yj-yi，lz0=zj-zi。

3)對索原長L0和端點(diǎn)力(F1，F(xiàn)2，F(xiàn)3)進(jìn)行初始化：

(21)

4)使用式(6)更新(lx，ly，lz)。

5)計(jì)算閉合差向量

6)如果dL小于給定的允許值，現(xiàn)有的L0是索原長。如果dL大于預(yù)設(shè)的允許值，程序繼續(xù)進(jìn)行下一步。

7)將式(6)分別對F2，F(xiàn)3和L0進(jìn)行微分，得到矩陣D。

(22)

9)更新F2=F2+dF2,F3=F3+dF3,L0=L0+dL0，然后回到4)。

對于非線性靜力分析，每一個(gè)荷載增量的殘余力使用牛頓迭代法進(jìn)行離散，對于非線性時(shí)程分析，采用增量迭代法，以Newmark直接積分法和牛頓迭代法為基礎(chǔ)，來求解非線性運(yùn)動(dòng)方程。結(jié)構(gòu)的增量運(yùn)動(dòng)方程可以寫成如下的形式：

(23)

(24)

式中：αM和αK分別是質(zhì)量和剛度比例阻尼系數(shù)，采用將Newmark法中常平均加速度法(γ=1/2,β=1/4)，增量加速度和速度在每一時(shí)間步的第一次迭代可以寫成：

(25)

(26)

式中：n為迭代次數(shù)，將方程(25)和(26)代入到方程(23)中，增量位移可表示為

(27)

(28)

(29)

在每一時(shí)間步的第一次迭代，基于增量位移向量Δu，對時(shí)間t+Δt時(shí)的總位移，速度和加速度進(jìn)行更新。

(30)

(31)

(32)

對于每個(gè)時(shí)間步的第二次以及后續(xù)的迭代，可以求解結(jié)構(gòu)體系在殘余力ΔR作用下的解。

(33)

根據(jù)總外力F、慣性力、阻尼力和更新的內(nèi)力Fint計(jì)算殘余力ΔR為

(34)

一旦計(jì)算滿足收斂準(zhǔn)則，結(jié)構(gòu)響應(yīng)在下一步被更新為

(35)

(36)

(37)

(38)

3數(shù)值驗(yàn)證

3.1索網(wǎng)結(jié)構(gòu)

此算例來自文獻(xiàn)[4]。該問題不考慮轉(zhuǎn)動(dòng)自由度，只需將外力矩設(shè)為零即可。圖3為索網(wǎng)的離散圖，圖4為索網(wǎng)的初始幾何形狀。每根索的截面面積為1.48 cm2，彈性模量為82.68 GPa，標(biāo)號為3、4、8和11的水平位置的四根索的預(yù)應(yīng)力為37.61 MPa，其余傾斜的索的預(yù)應(yīng)力為36.69 MPa。采用12個(gè)兩節(jié)點(diǎn)非線性空間索單元?jiǎng)澐执怂骶W(wǎng)，豎直集中載荷同時(shí)作用在4、5、8和9 4個(gè)連接點(diǎn)上，大小為都為55.12 MPa，計(jì)算時(shí)取相對誤差限為0.01，為簡潔起見，表1只給出了典型連接點(diǎn)4的位移分量，同時(shí)給出了文獻(xiàn)[4]采用五節(jié)點(diǎn)索單元的計(jì)算結(jié)果。

圖3　索網(wǎng)離散圖Fig.3　The discrete figure of cable-nets

圖4　索網(wǎng)初始幾何形狀Fig.4　The initial geometry of cable-nets

分量本文結(jié)果文獻(xiàn)[4]結(jié)果X-1.4642-1.4644Y-0.1318-0.1318Z0.13160.1318

從表1可以看出，采用空間兩節(jié)點(diǎn)懸鏈線索單元與采用五節(jié)點(diǎn)索單元計(jì)算的結(jié)果非常接近，說明應(yīng)用空間兩節(jié)點(diǎn)懸鏈線索單元來模擬索網(wǎng)結(jié)構(gòu)的方法和程序是正確可靠的。

3.2輸電線動(dòng)力響應(yīng)分析算例

3.2.1算例1

索結(jié)構(gòu)在風(fēng)荷載作用下表現(xiàn)出很強(qiáng)非線性特性[3,16]，采用本文提出的兩節(jié)點(diǎn)非線性空間索單元計(jì)算索結(jié)構(gòu)的非線性時(shí)程反應(yīng)。為驗(yàn)證該方法的正確性，同時(shí)應(yīng)用ANSYS有限元分析軟件進(jìn)行計(jì)算，選取在實(shí)際工程中使用的500 kV貓頭塔，該輸電塔總高為56.2 m，呼稱高為42 m，跨距為300 m。利用空間梁單元模擬塔的桿件，將導(dǎo)線看成理想的柔性索，用桿單元來模擬導(dǎo)線，建立輸電塔-線耦合的空間有限元模型，塔底和線端處理成固結(jié)。塔的材料阻尼比取0.03，如圖5所示。導(dǎo)線選用4×LGJ-400/35型鋼芯鋁絞線，將四分裂導(dǎo)線簡化為一根導(dǎo)線建模，根據(jù)GB/T 1179-2008《圓線同心絞架空導(dǎo)線》[17]查得技術(shù)參數(shù)如表2所示。

圖5　輸電塔-線耦合的空間有限元模型Fig.5　The coupled transmission tower-line system

導(dǎo)線根數(shù)橫截面積/mm2計(jì)算抗拉斷力/N單位長度質(zhì)量/(kg·km-1)直徑/mm單根425.24103900134923.27四根1700.96415600539646.55

本文采用AR法[18]同時(shí)模擬了78點(diǎn)的空間脈動(dòng)風(fēng)速，但由于輸電塔屬于高聳結(jié)構(gòu)，所以縱向脈動(dòng)風(fēng)速譜選用沿高度變化的Kaimal譜，自然風(fēng)場對結(jié)構(gòu)離散化節(jié)點(diǎn)的脈動(dòng)作用在空間上存在位置相干性，當(dāng)輸電線路位于陡坡地段或桿塔不等高時(shí)，導(dǎo)線落差梯度較大，或由于導(dǎo)線跨距較大，導(dǎo)線垂度效應(yīng)明顯時(shí)，各離散點(diǎn)的平均風(fēng)速和風(fēng)速譜均不相同，應(yīng)同時(shí)考慮了橫向與豎向的空間相關(guān)性。

(39)

式中：SV是功率譜密度，f是頻率，V*是剪切速度，f*為無量綱的Monin坐標(biāo)。不同高度處的風(fēng)速作用的相位是不同，一般先作用于結(jié)構(gòu)的較高處，然后到達(dá)較低的結(jié)構(gòu)，作用時(shí)差為f, 可以用互功率譜密度函數(shù)來表達(dá)這些信息，互譜函數(shù)一般為復(fù)數(shù)形式：

(40)

圖6　風(fēng)速時(shí)程Fig.6　Wind velocity time-history

圖7　脈動(dòng)風(fēng)速功率譜密度Fig.7　The power spectral density of wind velocity

假設(shè)模型所處位置標(biāo)準(zhǔn)高度(10 m) 處的平均風(fēng)速為V10= 10 m/s, 地面粗糙度系數(shù)k=0.005；脈動(dòng)風(fēng)速時(shí)間間隔為0.05 s，共模擬了600 s的風(fēng)速數(shù)據(jù),風(fēng)速時(shí)程曲線如圖6所示，圖7為脈動(dòng)風(fēng)速功率譜密度，可以看出模擬譜與目標(biāo)譜吻合較好，可以用于本文中的時(shí)程分析。經(jīng)過分析，得出了子導(dǎo)線1中點(diǎn)的位移時(shí)程曲線，如圖8，圖中實(shí)線給出了由本文方法得到的子導(dǎo)線1中點(diǎn)位移時(shí)程曲線。虛線為有限元分析軟件ANSYS的計(jì)算結(jié)果，從兩條曲線的分布情況可以看

出,本文方法和ANSYS計(jì)算得到的子導(dǎo)線1中點(diǎn)處的豎向位移、橫向位移結(jié)果基本一致，但是在計(jì)算時(shí)間上，采用ANSYS計(jì)算在4 248 s時(shí)達(dá)到收斂，采用本文方法的計(jì)算達(dá)到收斂時(shí)間為700 s, 因此，應(yīng)該空間兩節(jié)點(diǎn)索單元進(jìn)行計(jì)算，在滿足計(jì)算的正確性和精度的同時(shí)，可以節(jié)省計(jì)算時(shí)間，降低計(jì)算代價(jià)。

3.2.2算例2

文獻(xiàn)[19]中采用三節(jié)點(diǎn)曲線梁法分析單跨輸電線路，導(dǎo)線為單導(dǎo)線，輸電線路主要是由輸電塔、輸電線、絕緣子、間隔棒以及其他連接金具構(gòu)成，對于一個(gè)輸電線路的截?cái)嗄Ｐ?，輸電塔、絕緣子串和相鄰跨的剛度可以用簡化的彈簧模型來表示，如圖9(a)所示。Kx、Ky、Kz分別表示輸電線兩端輸電塔、絕緣子以及相鄰子跨在x、y、z方向上的簡化邊界條件。導(dǎo)線選用LGJ-400/50型鋼芯鋁絞線，檔距為250 m，弧垂為1.3 m，跨中單元的張力為106.8 kN，通過懸鏈線理論對導(dǎo)線進(jìn)行找形分析。由于覆冰的作用，導(dǎo)線覆冰形狀為新月形，如圖9(b)是覆冰輸電導(dǎo)線的斷面圖，導(dǎo)線直徑為27.6 mm，覆冰厚度為20 mm。導(dǎo)線的軸向剛度為EA=31.1×106N，扭轉(zhuǎn)剛度為GJ=159 Nm2/rad，抗彎剛度為1 965 Nm2，導(dǎo)線的彈性模量為69 000 N/mm2，覆冰導(dǎo)線參數(shù)見表3，新月形覆冰導(dǎo)線的氣動(dòng)力系數(shù)曲線由風(fēng)洞試驗(yàn)測得，分別測試了10 m/s風(fēng)速下，風(fēng)攻角在-180°～180°范圍內(nèi)的升力系數(shù)、阻力系數(shù)和扭矩系數(shù)，如圖10。

圖8　風(fēng)荷載作用下子導(dǎo)線1中點(diǎn)位移變化曲線Fig.8　The displacement at midpoint of subconductor 1 under wind loads

參數(shù)數(shù)值覆冰厚度/mm20裸導(dǎo)線橫截面積/mm2598.28覆冰面積/mm2343.04裸導(dǎo)線單位長度質(zhì)量/(kg·m-1)1.51覆冰導(dǎo)線單位長度質(zhì)量/(kg·m-1)0.308覆冰導(dǎo)線的偏心距/mm19.17

在10 m/s風(fēng)速，風(fēng)攻角為-10°情況下，選取時(shí)間步長為0.001 s，總的仿真時(shí)間為700 s，經(jīng)過計(jì)算，得出導(dǎo)線動(dòng)力響應(yīng)位移時(shí)程，截取680～690 s的子導(dǎo)線1中點(diǎn)的位移時(shí)程曲線如圖11所示。

圖11　跨中節(jié)點(diǎn)的位移時(shí)程曲線Fig.11　The time-history curve of displacement at mid-span node

圖11中實(shí)線給出了文獻(xiàn)[19]中采用三節(jié)點(diǎn)曲線梁法來模擬的索單元得到的子導(dǎo)線1中點(diǎn)橫向位移、豎向位移時(shí)程曲線。虛線為借助本文方法得出的結(jié)果，從兩條曲線的分布情況可以看出，用本文方法和三節(jié)點(diǎn)曲線梁法計(jì)算得到的跨中節(jié)點(diǎn)處的橫向位移、豎向位移的結(jié)果達(dá)到非常好的吻合，較好的驗(yàn)證本文方法及程序的正確性，充分考慮了導(dǎo)線的非線性效應(yīng)，計(jì)算中需要的自由度數(shù)更少，有利于節(jié)省計(jì)算時(shí)間。

4結(jié)論

針對現(xiàn)有索結(jié)構(gòu)模型在模擬非線性方面的不足，本文基于彈性懸鏈線解析理論開發(fā)了一種空間兩節(jié)點(diǎn)懸鏈線索單元。經(jīng)過計(jì)算，得出的結(jié)論主要有以下幾點(diǎn)：

1)采用構(gòu)造的空間兩節(jié)點(diǎn)懸鏈線索單元編制計(jì)算機(jī)程序?qū)〈馆^小的索網(wǎng)結(jié)構(gòu)和弧垂較大的輸電線索結(jié)構(gòu)進(jìn)行計(jì)算分析?？梢跃_有效的計(jì)算懸索結(jié)構(gòu)受到的靜態(tài)和動(dòng)態(tài)荷載，并給出了相應(yīng)的切線剛度矩陣和內(nèi)力矢量，索網(wǎng)中的每個(gè)索單元可以通過一個(gè)單元來建模，通過與現(xiàn)有研究結(jié)果比較發(fā)現(xiàn)，本文提出的兩節(jié)點(diǎn)懸鏈線索單元可以很好地預(yù)測懸索結(jié)構(gòu)的線性和非線性行為。數(shù)值算例分析驗(yàn)證了本文方法的有效性和計(jì)算程序的準(zhǔn)確性。

2)采用兩節(jié)點(diǎn)索單元進(jìn)行索結(jié)構(gòu)計(jì)算無需太多的結(jié)構(gòu)單元及自由度數(shù)即可得到精度較高的計(jì)算結(jié)果，對節(jié)省計(jì)算時(shí)間降低計(jì)算代價(jià)有利。編制的計(jì)算機(jī)程序具有良好的通用性，不僅適用于輸電線路這類弧垂較大的索結(jié)構(gòu)的分析，也可以用于大跨度的懸索和索網(wǎng)結(jié)構(gòu)的分析。

參考文獻(xiàn)：

[1]陳太聰, 蘇成, 馬海濤. 索結(jié)構(gòu)靜力分析的新型拉索單元研究[J]. 工程力學(xué), 2014, 31(2): 46-52.CHEN Taicong, SU Cheng, MA Haitao. Study on the development of new-type cable elements for static analysis of cable structures[J]. Engineering mechanics, 2014, 31(2): 46-52.

[2]魏建東. 索結(jié)構(gòu)分析的滑移索單元法[J]. 工程力學(xué), 2004, 21(6): 172-176, 210.

WEI Jiandong. Sliding cable element for the analysis of cable structures[J]. Engineering mechanics, 2004, 21(6): 172-176, 210.

[3]唐建民, 董明, 錢若軍. 張拉結(jié)構(gòu)非線性分析的五節(jié)點(diǎn)等參單元[J]. 計(jì)算力學(xué)學(xué)報(bào), 1997, 14(1): 108-113.

TANG Jianmin, DONG Ming, QIAN Ruojun. A finite element method with five-node isoparametric element for nonlinear analysis of tension structures[J]. Chinese journal of computational mechanics, 1997, 14(1): 108-113.

[4]黃濤, 周金枝, 朱若燕. 具有扭轉(zhuǎn)自由度的五節(jié)點(diǎn)索單元非線性有限元法[J]. 武漢理工大學(xué)學(xué)報(bào):交通科學(xué)與工程版, 2010, 34(6): 1150-1153, 1158.

HUANG Tao, ZHOU Jinzhi, ZHU Ruoyan. A finite element method with five-node cable element with a axial rotational freedom[J]. Journal of Wuhan university of technology:transportation science & engineering, 2010, 34(6): 1150-1153, 1158.

[5]LIEW J Y R, PUNNIYAKOTTY N M, SHANMUGAM N E. Limit-state analysis and design of cable-tensioned structures[J]. International journal of space structures, 2001, 16(2): 95-110.

[6]GRECO L, IMPOLLONIA N, CUOMO M. A procedure for the static analysis of cable structures following elastic catenary theory[J]. International Journal of Solids and Structures, 2014, 51(7/8): 1521-1533.

[7]ABAD M S A, SHOOSHTARI A, ESMAEILI V, et al. Nonlinear analysis of cable structures under general loadings[J]. Finite elements in analysis and design, 2013, 73: 11-19.

[8]CHEN Z H, WU Y J, YIN Y, et al. Formulation and application of multi-node sliding cable element for the analysis of Suspen-Dome structures[J]. Finite elements in analysis and design, 2010, 46(9): 743-750.

[9]JAYARAMAN H B, KNUDSON W C. A curved element for the analysis of cable structures[J]. Computers & structures, 1981, 14(3/4): 325-333.

[10]WANG Chunjiang, WANG Renpeng, DONG Shilin, et al. A new catenary cable element[J]. International journal of space structures, 2003, 18(4): 269-275.

[11]ANDREU A, GIL L, ROCA P. A new deformable catenary element for the analysis of cable net structures[J]. Computers & structures, 2006, 84(29/30): 1882-1890.

[12]YANG Y B, TSAY J Y. Geometric nonlinear analysis of cable structures with a two-node cable element by generalized displacement control method[J]. International journal of structural stability and dynamics, 2007, 7(4): 571-588.

[13]SUCH M, JIMENEZ-OCTAVIO J R, CARNICERO A, et al. An approach based on the catenary equation to deal with static analysis of three dimensional cable structures[J]. Engineering structures, 2009, 31(9): 2162-2170.

[14]THAI H T, KIM S E. Nonlinear static and dynamic analysis of cable structures[J]. Finite elements in analysis and design, 2011, 47(3): 237-246.

[15]劉克同, 湯愛平, 劉玥君, 等. 基于MRT-LBM大渦模擬的橋梁氣動(dòng)參數(shù)數(shù)值仿真[J]. 四川大學(xué)學(xué)報(bào):工程科學(xué)版, 2013, 45(6): 87-95.

LIU Ketong, TANG Aiping, LIU Yuejun, et al. Numerical simulation for aerodynamic parameters of bridge decks using MRT-LBM-LES[J]. Journal of Sichuan university:engineering science edition, 2013, 45(6): 87-95.

[16]全國電線電纜標(biāo)準(zhǔn)化委員會(huì). GB/T 1179-2008圓線同心絞架空導(dǎo)線[S]. 北京: 中國標(biāo)準(zhǔn)出版社, 2009: 15-17.

[17]OWEN J S, ECCLES B J, CHOO B S, et al. The application of auto-regressive time series modelling for the time-frequency analysis of civil engineering structures[J]. Engineering structures, 2001, 23(5): 521-536.

[18]YAN Zhitao, LI Zhengliang, SAVORY E, et al. Galloping of a single iced conductor based on curved-beam theory[J]. Journal of wind engineering and industrial aerodynamics, 2013, 123: 77-87.

本文引用格式：

劉玥君，湯愛平，劉克同. 索網(wǎng)結(jié)構(gòu)中的空間兩節(jié)點(diǎn)懸鏈線索單元[J]. 哈爾濱工程大學(xué)學(xué)報(bào)， 2016, 37(6): 788-795.

LIU Yuejun, TANG Aiping, LIU Ketong. Spatial two-node catenary element for analyzing cable-net structures[J]. Journal of Harbin Engineering University, 2016, 37(6): 788-795.

收稿日期：2015-03-27.

基金項(xiàng)目：國家國際科技合作專項(xiàng)項(xiàng)目(2011DFA21460).

作者簡介：劉玥君(1983-)，女，講師； 通信作者：劉玥君，E-mail: liuyuejun0221@163.com.

DOI:10.11990/jheu.201503081

中圖分類號：O242.2

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

文章編號：1006-7043(2016)06-0788-08

Spatial two-node catenary element for analyzing cable-net structures

LIU Yuejun1, TANG Aiping2, LIU Ketong2

(1.School of Civil and Architecture Engineering, Northeast Dianli University, Jilin 132012, China; 2. School of Civil Engineering, Harbin Institute of Technology, Harbin 150090, China)

Abstract：To analyze the linear and nonlinear behaviors of cable-net structures based on their mechanical characteristics, this paper presents a spatial two-node catenary cable element model (STCCE), which can be used to analyze the nonlinear effect cable structures under static and dynamic loads. The tangent stiffness matrix was derived based on the analytical solutions of elastic catenary; the stiffness matrix and internal force vector were calculated using an iterative algorithm; the nonlinear behavior was analyzed using the new model STCCE under static and dynamic loading. The Newton-Raphson method was applied to sove the structural balance equation for the nonlinear static analysis. An incremental iterative method, based on the Newmark direct integration and Newton-Raphson methods, was used for nonlinear dynamic time-history analysis. The accuracy and reliability of the STCCE were verified using numerical analysis.

Keywords：spatial two-node element； catenary cable element；static load； dynamic load；nonlinear analysis

網(wǎng)絡(luò)出版日期：2016-04-21.

湯愛平(1968-)， 男，教授，博士生導(dǎo)師.