考慮彎曲剛度的懸索自由振動解析解

肖 一，卓衛(wèi)東，范立礎(chǔ)

(1.福州大學(xué) 土木工程學(xué)院，福建 福州 350108；2.同濟大學(xué) 橋梁工程系，上海 200092)

懸索的自由振動是一個古老的問題，而自1940年Tacoma橋風(fēng)毀后，其又吸引了眾多學(xué)者的關(guān)注。Pugsley，Saxon，Goodey等人[1]對懸鏈的自由振動進(jìn)行了深入研究 ，Simpson[2]，Irvine[3]，Triantafyllou[4]等人則推導(dǎo)了彈性索的自由振動解析解，尤其是Irvine的工作對其后的研究影響深遠(yuǎn)。但是，這些研究均未考慮彎曲剛度的影響。為更準(zhǔn)確地計算懸索的自振頻率，有時應(yīng)計入彎曲剛度的影響，尤其對于截面較大的拉索[5-6]，而關(guān)于振動法測索力的研究[7-9]一般均考慮了彎曲剛度的影響。文獻(xiàn)[10-11]則用半解析方法對考慮彎曲剛度的斜索固有振動問題進(jìn)行了研究。

最近，文獻(xiàn)[6]在拋物線假設(shè)下推導(dǎo)了考慮彎曲剛度的懸索面內(nèi)自由振動解析解，但其僅考慮了兩端固接邊界條件；文獻(xiàn)[12]在其已有研究基礎(chǔ)上計入彎曲剛度影響推導(dǎo)了斜索面內(nèi)固有振動解析解，但其僅考慮了兩端鉸接邊界條件。從文獻(xiàn)[13]的圖中可知，錨固端的約束條件可能為固接(直接錨固連接)，也可能為鉸接(如錨頭和錨板采用鉸接形式)，如文獻(xiàn)[5，14]就對兩種情況進(jìn)行了分析，而實際的情況往往是處于兩種狀態(tài)之間[8]。

當(dāng)不考慮彎曲剛度對靜力構(gòu)形的影響時，懸鏈線是較拋物線更接近實際情況的表達(dá)。本文以懸鏈線為索靜力構(gòu)形，分別考慮兩端固結(jié)和兩端鉸接兩種邊界條件對考慮彎曲剛度影響的懸索自由振動解析解進(jìn)行了推導(dǎo)，對已有的結(jié)果進(jìn)行了補充并糾正了文獻(xiàn)[6，12]中的不當(dāng)之處，最后對結(jié)果進(jìn)行了討論。

1 理論推導(dǎo)

如圖1所示的做微幅運動的空間懸索，其坐標(biāo)原點置于左端點。

圖1 懸索計算圖示

當(dāng)考慮彎曲剛度時，在不計阻尼和動力荷載情況下，懸索的面內(nèi)和面外線性運動方程分別為[6]:

其中:w和v分別為索z和y方向運動分量；H為索靜拉力水平分量；h為H的動力增量；EI為截面抗彎剛度；m為單位長度質(zhì)量。這里假設(shè)索截面在面內(nèi)和面外有相同的抗彎剛度。由式(1)和式(2)可知，面內(nèi)和面外運動不耦合。

由于索以承受軸向力為主，這里僅考慮由軸力產(chǎn)生的應(yīng)變而忽略彎曲應(yīng)變，此時有線性化的索方程[3]:

其中:u為x方向運動分量；s表示弧長。對其進(jìn)行無量綱化并考慮到兩端為固定的邊界條件可得:

其中，λ2=(mgl/H)2l/(HLe/EA)為 Irvine[3]引入的參數(shù)，其反映了索彈性效應(yīng)和幾何效應(yīng)的影響；Le為索有效長度[3]；EA為索軸向剛度。

1.1 面外振動

文獻(xiàn)[6，12]中并未給出面外振動解，這里首先對其做一個補充，以便于后文的討論。

易知式(4)的通解為:

其中:

兩端固接時有:

兩端鉸接時有:

可見，其解與受拉梁的結(jié)果[15-16]一致。這也體現(xiàn)了線性振動時，垂度效應(yīng)對面外運動是無影響的。

1.2 面內(nèi)振動

由于結(jié)構(gòu)和邊界條件是對稱的，分別按反對稱和對稱兩種情況討論。當(dāng)振型為反對稱形式時=0[3，6，12]，此時式(3)的通解為:

文獻(xiàn)[6]也給出了兩端固接索面內(nèi)反對稱振動解析解，通過對比發(fā)現(xiàn)，其頻率方程與本文相同，但振型函數(shù)并不滿足兩端的邊界條件。經(jīng)驗證，本文結(jié)果是正確的。

這與文獻(xiàn)[12]中的結(jié)果是不同的，由于采用了不合適的邊界條件，其頻率和振型中實際上包含對稱和反對稱部分，而式(15)和式(16)則僅為反對稱部分。

這里順便指出，由于實際上沒有考慮垂度的影響，這里求出的反對稱解實際上就是受拉梁在相應(yīng)邊界條件下的反對稱振動部分，當(dāng)把上面的中點處的邊界條件換為(1/2)=(1/2)=0時，就可得到受拉梁固有振動解的對稱部分。當(dāng)兩端固接時有:

兩端鉸接時則有:

把相應(yīng)的對稱和反對稱頻率方程相乘即可得到受拉梁固有振動的完整解，由于這里假設(shè)面內(nèi)和面外彎曲剛度相同，其就是式(8)和式(10)，這也進(jìn)一步體現(xiàn)了上面結(jié)果的正確性。

當(dāng)振型為對稱形式時，如圖1所示，為求解方便，首先將坐標(biāo)原點移至弦的中點，并沿用原來的符號表示。與文獻(xiàn)[6，12]采用的拋物線不同的是，這里采用更為準(zhǔn)確的懸鏈線索形。此時，懸鏈線表達(dá)[17]經(jīng)坐標(biāo)變換和無量綱化后得:

其中，k=mgl/H。由此可得式(3)的一個特解為:

則其通解可表示為:

而C2=C4=0，將其帶入通解中即可得到振型函數(shù)。然后，將所得振型函數(shù)帶入索方程中(此時由于坐標(biāo)系的移動，需要修改積分限)并消去≠0后，可得頻率方程為:

而C2=C4=0，將其帶入通解中即可得到振型函數(shù)，而相應(yīng)的頻率方程為:

值得指出的是，式(24)和式(26)只是一種表達(dá)形式 ，利用關(guān)系式 ξ2=α2-β2，ξ?ω=α β 和k4-ξ2k2-ξ2=(k2-α2)(k2+β2)，其還可以有其它等效的表達(dá)形式。

2 討 論

文獻(xiàn)[6]和文獻(xiàn)[12]分別對兩種邊界條件下的結(jié)果進(jìn)行了數(shù)值驗證，在λ2的實際取值范圍內(nèi)都有良好的精度。其中，文獻(xiàn)[12]的數(shù)值方法采用的索靜力構(gòu)形與本文相同。但實際上，由于他們得到的是解析解，這種驗證不但是對解的驗證，更是對近似運動方程的驗證，所以這里就不再對此進(jìn)行重復(fù)。

2.1 邊界條件的影響

圖2表示參數(shù) ξ分別取20、50、75和200[5-6]四種情況時，兩種邊界條件下前六階無量綱頻率(?ω/π)隨參數(shù) λ2的變化情況。CC表示兩端固接，HH表示兩端鉸接，Irvine表示文獻(xiàn)[3]的解。從圖中可見，考慮彎曲剛度的懸索同樣會發(fā)生Crossover現(xiàn)象[6，10]；顯然，考慮彎曲剛度的頻率要高于不考慮彎曲剛度的情況，兩端固結(jié)情況則要高于兩端鉸接情況；頻率階數(shù)越高，三者之間的差別越大；參數(shù) ξ越小，也就是彎曲剛度的比重越大時，兩種邊界條件下的頻率差別越大。

由此可見，對于需要考慮彎曲剛度的懸索的動力分析，考慮不同形式的邊界條件是必要的。

另外，正如文獻(xiàn)[8]指出的，拉索的實際約束情況往往處于固接和鉸接之間。本文的兩端固結(jié)和兩端鉸接情況相當(dāng)于實際約束情況的上下界，當(dāng)需要考慮更精確的約束條件時，由此可以對實際的拉索頻率進(jìn)行近似估計。

2.2 懸鏈線解與拋物線解對比

文獻(xiàn)[6，12]采用的靜力解實際上與Irvine[3]相同，其并未考慮彎曲剛度的影響。在忽略該影響的假定下，由于本文采用的是更準(zhǔn)確的懸鏈線索，所以從理論上講，其精度要高于文獻(xiàn)[6，12]的結(jié)果。

令文獻(xiàn)[12]中的參數(shù)ε=0時，就可得到相應(yīng)的水平索的頻率方程。這樣，在本文的無量綱表示下，考慮彎曲剛度的拋物線索的兩端固接和兩端鉸接面內(nèi)對稱頻率方程分別為[6，12]:

相比于其原始表示，式(27)、式(28)略為簡潔。

圖2 兩種邊界條件下參數(shù)ξ取不同值時前六階頻率對比

圖3 兩種邊界條件下參數(shù)k取不同值時前四階對稱頻率對比

圖3給出了 ξ=75時 ，k取 0.1、0.5、0.85和1.0時與拋物線解的前四階對稱無量綱頻率(?ω/π)隨參數(shù)λ2的變化情況。由圖3可見，兩種構(gòu)形表達(dá)的結(jié)果相差很小，僅上升段略有差別，其基本可以忽略(所以圖中未標(biāo)出各曲線的意義)，對于ξ取其它值的情況，也可得出類似的結(jié)果，這從另一個角度也驗證了式(27)和式(28)的正確。所以，相比而言，在實際應(yīng)用中可以采用式(27)和式(28)。

3 結(jié) 語

本文在已有研究基礎(chǔ)上考慮了更準(zhǔn)確的懸索靜力構(gòu)形——懸鏈線，推導(dǎo)了在兩端固接和兩端鉸接邊界條件下考慮彎曲剛度的懸索面外和面內(nèi)自由振動解析解，對已有文獻(xiàn)結(jié)論進(jìn)行了補充并修正了其不當(dāng)之處。通過討論認(rèn)為，當(dāng)分析中需要計入懸索的彎曲剛度時，應(yīng)該要考慮不同邊界條件的影響，而本文結(jié)果可以看作是實際情況的上下界。通過與已有結(jié)果對比發(fā)現(xiàn)，本文結(jié)果雖從理論上更精確，但其與已有結(jié)果的差別在實際應(yīng)用中基本可以忽略，從這個角度來講，本文結(jié)果也可以看作是對已有理論的一個補充。

[1]Rega G.Nonlinear dynamics of suspended cables，Part I.Modeling and analysis[J].Appl.Mech.Rev.，ASME，2004，(57):443-514.

[2]Simpson A.Determination of the inplane natural frequencies of multispan transmission lines by a transfer matrix method[C]//Proceedings of the IEEE，1966，113(5):870-878.

[3]Irvine H M，Caughey T K.The linear theory of free vibrations of a suspended cable[C]//Proceedings of the Royal Society of London Series A，1974，341:299-315.

[4]Triantafyllou M S.The dynamics of taut inclined cables[J].Quarterly Journal of Mechanics and Applied Mathematics，1984，37(3):421-440.

[5]Ni Y Q，Ko J M，Zheng G.Dynamic analysis of large-diameter sagged cables taking into account flexural rigidity[J].Journal of Sound and Vibration，2002 ，257(2):301-319.

[6]Ricciardi G，Saitta F.A continuous vibration analysis model for cables with sag and bending stiffness[J].Engineering Structures，2008，30(5):1459-1472.

[7]Zui H，Shinke T，Namita YH.Practical formulasfor estimation of cable tension by vibration method[J].Journal of Structural Engineering，ASCE，1996，122(6):651-656.

[8]Kim G H，ParkT.Estimation of cable tensionforce using the frequency-based system identification method[J].Journal of Sound and Vibration，2007，304(3-5):660-676.

[9]任偉新，陳 剛.由基頻計算拉索拉力的實用公式[J].土木工程學(xué)報，2005，38(11):26-31.

[10]趙躍宇，周海兵，金 波，等.彎曲剛度對斜拉索非線性固有頻率的影響[J].工程力學(xué)，2008，25(1):196-202.

[11]吳 曉，黎大志，羅佑新.斜拉索非線性固有振動特性分析[J].振動與沖擊，2003 ，22(3):37-39，11.

[12]吳慶雄，李 瀏，陳寶春.考慮彎曲剛度的拉索面內(nèi)固有振動的理論計算公式[J].工程力學(xué)，2010，27(11):9-15，27.

[13]嚴(yán)國敏.現(xiàn)代斜拉橋[M].成都:西南交通大學(xué)出版社，1996:60-68.

[14]Mehrabi A B，Tabatabai H.Unified finite difference formulation for free vibration of cables[J].Journal of Structural Engineering，ASCE，1998，124(11):1313-1322.

[15]Bokaian A.Natural frequenciesof beams under tensile axial loads[J].Journal of Sound and Vibration，1990，126:481-498.

[16]Humar J L.Dynamicsof Structures[M].AA Balkema Publishers，2002:829-876.

[17]Pugsley A.The theory of suspension bridges[M].Edward Arnold，Paris，1957:21-25.