影響不對稱支撐圓形鋼天線風致疲勞壽命的因素

王欽華,顧 明

(1.同濟大學土木工程防災國家重點實驗室,上海 200092;2.汕頭大學 土木工程系,廣東 汕頭 515063）

對于輕質、柔性結構,風荷載是其主要的控制荷載,結構在風荷載的作用下產生較大幅度的振動,而結構在經歷了大量的應力循環后,易發生風振疲勞破壞。風振疲勞破壞的例子較多[1-7],但是對于風振疲勞研究的文獻較少。Gu等[8]對上海地區的氣象資料進行統計分析得到了該地區風速風向聯合分布函數,并對楊浦大橋進行抖振響應分析,得到每個風向角各個風速下鋼箱梁關鍵點處的應力功率譜密度函數,用蒙特卡洛方法模擬了關鍵點處的應力時程,然后用雨流計數法對應力時程進行計數得到應力范圍的分布函數,在此基礎上對楊浦大橋進行風振疲勞分析。M.Dionne和A.G.Davenport[9]建立了陣風響應因子和風振疲勞之間的關系,其文獻將結構的響應分為背景響應和共振響應兩部分,用Miner累積損傷理論計算結構的疲勞壽命,這個方法只能用于由順風向力引起的疲勞問題。M.P.Repetto和G.Solari[10]將結構的應力響應分為背景響應和共振響應,用累積損傷理論對結構進行順風向疲勞壽命分析。此后,M.P.Repetto和 G.Solari[11-13]考慮了順風向力、橫風向力以及渦激共振的影響對圓形線狀結構進行風振疲勞壽命分析,但其文獻中沒有詳細討論風向、結構表面粗糙度以及平均風速對疲勞壽命的影響。文獻[14-15]基于一鋼質天線氣動彈性模型風洞試驗結果,分別用頻域和時域的方法研究了其風致疲勞。

該文在求得結構不同風向、風速下風致響應以及結構所在位置處的風速風向聯合分布函數的基礎上,基于經典的疲勞累積損傷理論,對不對稱支撐圓形截面鋼結構進行風致疲勞壽命估計,通過實例分析的結果,詳細討論了風向、結構表面粗糙度、渦激共振以及平均風速等重要因素結構風致疲勞壽命的影響。

1 結構風致疲勞壽命分析理論

結構風致疲勞壽命的分析方法可以分為3類:第1類是時域法,其目的是從用于疲勞分析的應力時程響應中,通過雨流計數法計數得到應力范圍或者應力幅值的概率密度分布函數,第2類是基于隨機振動和概率論的頻域方法,其目的是由應力功率譜響應通過穿越分析得到應力范圍或者應力幅值的概率分布密度函數。第3類主要針對考慮背景響應不可忽略的結構風致疲勞問題,為了縮減計算量,提高計算效率,先根據隨機振動理論求得用于疲勞分析的應力響應譜,再用Monte Carlo法模擬應力時程,計數后求得應力范圍或者幅值的概率密度函數,Gu等[8]定義這種方法為頻-時域方法。以上3類方法是基于線性疲勞累積損傷理論,其主要不足在于:損傷與荷載狀態無關;累積損傷與荷載次序無關;不能考慮荷載間的相互作用。

1.1 風致疲勞壽命時域分析方法

由于結構的疲勞累積損傷與風向角和風速有關,首先要定義工況,假設每一個風向角范圍[θi,θi+1]下風速范圍[vj,vj+1]是第[i,j]個工況,每個工況出現的概率為pij,可以由風速風向分布聯合函數求得;第[i,j]個工況下,單位時間內應力范圍出現次數的數學期望為,應力范圍的概率密度函數;單位時間內,第k級應力范圍出現的次數為:

式中,應力范圍的概率密度函數 fSij是采用雨流計數法對應力時程進行計數,根據計數結果擬合得到;單位時間內應力范圍出現次數的數學期望為:

式中,nij為用雨流法計數第[i,j]個工況下的應力時程得到的應力范圍的總次數;tij為第[i,j]個工況下應力時程的時間。

第[i,j]個工況下的累積損傷為:

結構在時間T內總的累積損傷為:

1.2 風致疲勞壽命頻域分析方法[15-16]

對于隨機平穩應力過程s(t),假設其應力幅值和峰值總數為相互獨立的隨機變量。如果該應力過程的應力范圍分布用連續的概率密度函數表示,根據S-N曲線和疲勞累積損傷準則,疲勞累積損傷可以改成下列的連續形式:

式中,dn是在區間dS內應力范圍的循環次數,它等于該區間內應力范圍出現的概率 fS(S)dS和應力幅值總次數mp的乘積;m、A是與材料性質有關的常數。

第[i,j]工況下,單位時間內造成的疲勞累積損傷D0ij為:

式中,flij為第[i,j]工況下單位時間內應力范圍總數的數學期望;fSLij(s)為應力范圍概率密度函數。以上2個參數可以通過疲勞分析關鍵點處的應力功率譜密度函數而求得[13]。

T時間內由第[i,j]工況下內造成的累積損傷為:

1.3 結構風致疲勞壽命分析流程

以上簡單介紹了風致疲勞壽命分析方法,該節給出進行結構風致疲勞壽命分析流程:

1)根據結構附近氣象站處的風速風向實測資料,得到結構處的風速風向聯合分布函數,計算每個風向風速區間出現的概率以及每個風向角下不同重現期對應的最大風速。由于在計算結構的風致疲勞壽命時,需要計算每個風速風向區間造成的疲勞累計損傷,每個風向區間下風速的最大值須根據疲勞壽命作為重現期來確定。

2)假設結構的風致疲勞壽命,根據假設的疲勞壽命計算每個風向區間要考慮的最大風速。

3)計算每個風向風速區間下,用于疲勞分析的應力響應。

4)基于疲勞累計損傷理論用時域方法或者頻域方法計算結構的疲勞壽命。

5)判斷假設的疲勞壽命與計算得到的疲勞壽命的差值是否滿足要求,如果是,結束;如果否,返回到第2)步后循環計算,直到假設的疲勞壽命與計算得到的疲勞壽命的差值滿足要求為止。

為了使結構風致疲勞壽命分析流程表述更加直觀,圖1給出了流程圖。

圖1 風振疲勞壽命分析流程圖

2 實例分析

以某一位于上海的超高層建筑上的天線為實例進行分析,該高層建筑主樓高度為245 m,主樓上部天線高為87 m。天線頂端離地高度為333m。天線由下往上分段等截面,每段高度依次為34m、27 m、26 m,其直徑分別為2.4m、1.75 m、1 m,壁厚各為3 cm、2.5 cm、2 cm。在距離底部34m處有一個斜支撐其直徑為1 m,見圖2所示。

圖2 天線結構圖

文獻[14]中給出了天線結構位置處的風速風向聯合分布函數,文獻[17]介紹了求解不同風向角、風速下用于疲勞壽命分析的結構響應的方法。該文用文獻[17]中的方法對天線結構疲勞關鍵點處進行應力響應分析,結構關鍵點的位置見文獻[14]。此外,為了考慮結構表面粗糙度對結構風致疲勞壽命的影響,根據文獻[18]結構表面粗糙度ε分別選用0.002、0.0012、0.0004。以22.5度風向間隔,將0 ～ 360度風向劃分為16個風向區間。不考慮斜撐對結構風場影響時,結構關于X、Y軸對稱,所以只計算0～90度風向下的響應即可,有關風向角的定義如圖3所示 。風速范圍間隔 Δ vj=vj+1-vj為1 m/s,v00,在分析過程中以。計算了0～90度風向角、所有風速下,結構關鍵點處的名義應力響應,由于篇幅的限制這里只給出幾個典型工況下的名義應力響應。圖4給出了不同風向、風速以及結構表面粗糙度下結構的名義應力響應譜。

圖3 風向角的定義

圖4 不同風向、風速以及結構表面粗糙度下結構的名義應力響應譜

由圖4(1)可知,在90度風向角4.5 m/s風速下,應力功率譜在0.475 Hz處出現唯一峰值,這是因為在該風向角和風速下結構發生一階渦激共振。(2)圖表明,在0度風向角6.5 m/s風速下結構發生二階渦激共振。由(3)圖可知,應力功率譜出現幾個峰值,但在0.475 Hz和0.675 Hz處峰值最大,這表明在此風向角下結構的一階和二階陣型對結構響應的貢獻較大。圖4中所有的圖表明在相同的風向角和風速下,結構表面粗糙度ε越大結構的響應就越大。

圖5給出了不同表面粗糙度下結構的應力根方差隨風速的變化規律。圖5中(1)-(4)圖表明不考慮渦激共振時結構的應力響應隨著風速的增大而增大;結構的粗糙度ε對結構的應力響應影響較大。在相同風速條件下,粗糙度ε越大,結構的應力響應也越大。

圖5 不同風速風向角風速下應力根方差

在求得不同風向角、風速下結構應力響應的基礎上,對結構進行風致疲勞壽命估計。在計算結構的風振疲勞壽命時,根據《鋼結構設計規范》(GB50017-2003)規定,構件連接類別選為8類,其S-N曲線參數為m=3,A=4.1×1011。以下討論了風向、結構表面粗糙度、渦激共振以及平均風速等重要因素對結構風致疲勞壽命的影響。

2.1 平均風速對疲勞壽命的影響

根據應力響應功率譜,用Monte Carlo方法模擬了應力時程響應,并用雨流法對應力時程響應進行計數,用Goodman法則考慮了平均風荷載引起的平均應力的影響,得到了應力范圍出現的概率,對天線進行風致疲勞壽命估計。

當表面粗糙度為0.0012時,考慮平均風速的影響,用時域法計算得到的結構的風致疲勞壽命為116.37 a;不考慮平均風速影響,用時域方法計算得到的風致疲勞壽命為117.95 a。從中可以看出平均風速對結構風振疲勞壽命影響較小。圖6-圖7給出了2個典型風向角和風速下,模擬的應力時程響應,用雨流法計數應力時程并考慮平均應力影響的應力范圍及其出現的次數以及應力范圍在結構疲勞壽命期內的累積損傷。圖6中(1)圖表明90度風向角4.5 m/s風速下,結構發生一階渦激共振的應力幅值為4 MPa左右,所以(2)圖中用雨流計數法得到的接近8 MPa的應力范圍(應力幅值的2倍)出現的次數最多,(3)圖說明以某一階振型振動產生的應力時程,產生疲勞累積損傷的應力范圍為應力幅值的2倍。圖7是考慮前5階振型后得到的應力時程,其(3)圖表明出現次數最多的小應力范圍產生的疲勞累積損傷較小,而出現次數相對較少的大應力范圍造成的累積損傷較大。

圖6 圖90度風向角4.5 m/s風速下雨流計數結果和累積損傷

圖7 圖45度風向角19.5 m/s風速下雨流計數結果和累積損傷

圖8 考慮渦激共振時累積損傷隨風向角的變化

圖9 考慮渦激共振時累積損傷隨風速的變化

圖8和圖9給出了疲勞累積損傷隨風向角和風速的變化規律,從圖8中可以看出風向角對疲勞累計損傷的影響較大。從圖9中可以看出4.5 m/s、6.5 m/s以及20～30 m/s風速下產生的疲勞累計損傷較大,這是因為:在4.5和6.5 m/s風速下結構產生渦激共振,20～30 m/s風速下引起結構較大的應力。

2.2 結構表面粗糙度對疲勞壽命的影響

表面粗糙度分別為0.0004、0.0012以及0.002時,用頻域法計算得到的結構風致壽命分別為142.57 a、101.54 a以及 91.62 a。由此可見,結構表面粗糙度對圓形單體結構的風致疲勞壽命影響很大,結構的疲勞壽命隨著表面粗糙度的增加而減小。這是因為當雷諾數大于4×104時,結構表面粗糙度對結構的阻力系數影響很大,阻力系數隨著結構表面粗糙度的增加而增大,因此結構的應力響應隨著結構表面粗糙度的增加而增大。圖10給出了風速風向區間出現的概率。圖11給出了結構表面粗糙度參數為0.0004、0.0012以及0.002時,各個風速風向工況下的風致疲勞累積損傷。

圖11 每個工況下不同結構表面粗糙度對應的累積損傷

從圖11中可以看出0～25°以及150～270°風向角工況造成的風致疲勞累積損傷較大。50～120°風向角以及280～300°風向角下的疲勞累積損傷較小,這和圖10表示的結構位置處風速風向區間出現的概率一致。這說明,風向對結構的風致疲勞累積損傷影響較大,在出現概率大的風向區間內造成的風致疲勞累積損傷較大。

圖11表明:在所有風向角下,風速為 25～30 m/s區間內的累計損傷較大;風速4.5 m/s,6.5 m/s時,結構發生渦激共振,此時的累積損傷也較大;大于38 m/s的風速產生的疲勞累積損傷為零,這是因為以結構風致疲勞壽命為重現期計算得到的每個風向區間下的最大風速都小于38 m/s;不同表面粗糙度結構在各個風速風向區間內疲勞累積損傷的趨勢一致。

2.3 渦激共振對疲勞壽命的影響

該節分別用頻域和時域的方法計算了不考慮以及考慮渦激共振影響的3種表面粗糙度下結構的風致疲勞壽命,其結果如表1所示。從表1中可知結構的渦激共振對結構風致疲勞壽命有一定程度的影響,是不可忽略的。圖12給出了不同結構表面粗糙度下不考慮渦激共振時用頻域方法計算的累計損傷隨風速和風向的變化曲線。從(2)圖中可以看出:當不考慮渦激共振影響時,小于10 m/s風速造成的疲勞累計損傷很小。

表1 考慮與不考慮渦激共振時得到的風致疲勞壽命的比較 a

圖12 不考慮渦激共振時累計損傷隨風速和風向的變化

3 結論

在求得結構不同風向、風速下風致響應以及結構所在位置處的風速風向聯合分布函數的基礎上,基于經典的疲勞累積損傷理論,對單體結構進行風致疲勞壽命估計。得到的幾個重要結論為:風向對結構的風致疲勞累積損傷影響較大,在出現概率大的風向區間內造成的風致疲勞累積損傷較大,因此在進行風致疲勞壽命分析時,需根據氣象資料得到建筑物位置處的風速風向聯合分布函數;結構表面粗糙度對圓形柱狀結構的風致疲勞壽命影響很大,結構的疲勞壽命隨著表面粗糙度的增加而減小;渦激共振對結構的疲勞壽命有一定的影響,在計算圓形柱狀結構的風致疲勞壽命時,不應該忽略渦激共振的影響;平均風荷載對結構的風致疲勞壽命影響較小。

[1]B N Pritchard.Steel chimney oscillations:a comparative study of their reported performence versus predictions using existing design techniques[J].Engineering Structres,1984,6(4):315-323.

[2]C VERWIEBE,A GLOCKNER.Failure of steel chimneys due to vortex excited vibration in the second mode[C]//Proceeding of 11th international conference on wind engineering.2003.Lubbok,Tex.

[3]鄧洪洲,屠海明,王肇民.桅桿結構隨機風振疲勞研究[J].土木工程學報,2003,36(4):19-23.DENG HONG-ZHOU,TU HAI-MING,WANG ZHAO-MIN.Study on stochastic fatigue of guyed masts under wind load[J].China Civil Engineering Journal,2003,36(4):19-23.

[4]張其林,U PEIL.歐美塔桅鋼結構現狀研究[J].特種結構,1996,13(2):58-63.ZHANG QI-LIN,U PEIL.Study about the current situation ofguyed maststructure in Europe and American[J].Special Structure,1996,13(2):58-63.

[5]P T RANVIK,G ALPSTEN.Dynamic behaviour under wind loading of a 90 m steel chimney[R].Alstom Power Sweden,2002.

[6]王世村,楊曉強.高聳結構風振疲勞的隨機裂紋擴展分析[J].鋼結構,2007.8(22):43-47.WANG SHI-CUN,YANG XIAO-QIANG.Stochastic analysis of fatigue crack growth for high-rise structures under wind loading[J].Steel Construction,2007,8(22):43-47.

[7]鄧洪州,王肈民,馬星.桅桿結構風振響應及疲勞研究進展[J].特種結構,2006,23(3):14-18.DENG HONG-ZHOU,WANG ZHAO-MIN,M A XING.State of the art paper about wind-induced response and fatigue of guyed masts[J].Special Structures,2006,23(3):14-18.

[8]GU M,XU Y L,CHEN L Z,et al.Fatigue life estimation of steel girder of yangpu cable-stayed bridge due to buffeting[J].Journal of Wind Engineering and Industrial Aerodynamics,1999,80:383-400.

[9]M CORAZAO,AG DAVENPORT.Asimple relationship between the gust response factorand fatigue damage[J].Journal of Wind Engineering and Industrial Aerodynamics,1988,30:45-54.

[10]M P REPETTO,G SOLARY.Dynamic alongwind fatigue of slender vertical structure[J].Engineering Structure,2001,23:1622-1633.

[11]M P REPETTO,G SOLARY.Dynamic crosswind fatigue of slender vertical structures[J].Wind and Structures,2002,5(6):527-542.

[12]M P REPETTO.Wind-induced response and fatigue of slender vertical structures[D].Genova The University of Genova,2003

[13]M PREPETTO,G SOLARY.Directional windinduced fatigue of slender vertical structure[J].Journal of Structural Engineering,2004,130(7):1032-1040.

[14]王欽華,顧明.鋼天線風振疲勞研究[J].同濟大學學報,2008.36(8):1040-1044.WANG QIN-HUA,GU MING.Directional windinduced fatigue of a real steel antenna[J].Journal of Tongji Unversity,2008,36(8):1040-1044.

[15]王欽華,顧明.鋼質天線結構風振疲勞壽命估算方法比較研究[J].振動與沖擊,2008,27(2):78-81.WANG QIN-HUA,GU M ING.Comparison between methods of structural fatigue life estimation[J].Journal of Vibration and Shock,2008,27(2):78-81.

[16]胡毓仁,陳伯真.船舶及海洋工程結構疲勞可靠性分析[M].北京:人民交通出版社,1996.

[17]王欽華,顧明.圓形線狀結構風致響應的分析方法[J].振動工程學報,2010,23(2):151-157.WANG QIN-HUA,GU MING.Wind-induced response analysis for circular section structure[J].Journal of Vibration Engineering,2010,23(2):151-157.

[18]ESDU.Mean force,pressures and flow velocities for circular cylindrical structures:finite-length cylinders in uniform and shear flow[R].ESDU International.1981,London.