劉海濤，魏明海，肖儀清，林 坤

(1.哈爾濱工業(yè)大學(xué) 深圳研究生院，廣東 深圳 518055;2.沈陽建筑大學(xué) 營造與工程管理系，沈陽 110168)

索－梁耦合結(jié)構(gòu)應(yīng)用廣泛。對其進行的傳統(tǒng)分析往往局限于單個結(jié)構(gòu)，即單個索或梁構(gòu)件［1－5］。索 －梁耦合結(jié)構(gòu)作為整體，系統(tǒng)的非線性行為不僅由索的幾何非線性引起，且亦因索、梁間模態(tài)耦合而產(chǎn)生。將索、梁作為整體用于斜拉結(jié)構(gòu)必存在多種內(nèi)、外共振形式及聯(lián)合形式。因此，對該耦合結(jié)構(gòu)非線性響應(yīng)進行深入研究具有重要意義。

索－梁耦合結(jié)構(gòu)源于Fujino等［6］對斜拉索參數(shù)振動的研究，其對模型進行一定簡化，僅考慮梁面內(nèi)、外及索面內(nèi)位移3個自由度，將索面外位移通過梁端約束耦合一起。Xia等［7］用解析方法考察拉索在端部隨機位移激勵的響應(yīng)，采用文獻［6］的3自由度索橋耦合振動簡化模型，由等效線性化方法獲得在主梁上作用豎向白噪聲荷載的系統(tǒng)響應(yīng)。Fung等［8］通過Hamilton原理建立索－梁耦合時變系統(tǒng)知，拉索張力、長度為隨時間變化的函數(shù)。Gattulli等［9－11］研究索 －梁組合結(jié)構(gòu)整體、局部與耦合模態(tài)的存在及相互間影響并用試驗及有限元進行驗證。王濤等［12］用拉索非線性解析振動方程與有限元非線性動力時程積分結(jié)合方法研究索、梁發(fā)生大幅振動的特性表明，拉索局部振動與整體結(jié)構(gòu)相關(guān)效應(yīng)較明顯。趙躍宇等［13－14］研究索－梁耦合結(jié)構(gòu)可能出現(xiàn)的內(nèi)共振模式，利用數(shù)值模擬方法研究該系統(tǒng)因內(nèi)共振引起的振動行為，并將索－梁耦合結(jié)構(gòu)擴展到索－曲梁耦合結(jié)構(gòu)，研究其面內(nèi)振動特性。馮維明等［15－16］將索－梁耦合結(jié)構(gòu)作為整體研究其耦合系統(tǒng)振動特性。而對索－梁耦合結(jié)構(gòu)在內(nèi)外共振聯(lián)合激勵下的非線性分析較少，尤其考慮模態(tài)耦合影響。為此，本文建立考慮模態(tài)相互耦合的索－梁耦合結(jié)構(gòu)面內(nèi)振動方程;用多尺度法對耦合運動方程進行解耦;并數(shù)值分析索－梁耦合結(jié)構(gòu)在內(nèi)外共振聯(lián)合激勵下的非線性特性及系統(tǒng)參數(shù)對非線性影響。

1 理論模型建立

索－梁耦合結(jié)構(gòu)簡化模型見圖1。由于研究耦合結(jié)構(gòu)面內(nèi)非線性特性，故假設(shè)條件為:① 不考慮梁的扭矩及大變形;② 不考慮索的抗剪剛度及抗彎剛度;③ 將索的重力垂度曲線近似為為拋物線;④ 用Lagrangian應(yīng)變描述索軸向伸長;⑤ 索質(zhì)量遠小于梁質(zhì)量;⑥ 索、梁變形本構(gòu)關(guān)系服從虎克定律且各點受力均勻。

圖1 索－梁耦合結(jié)構(gòu)模型Fig.1 The model of a cable-beam structure

在以上假設(shè)條件下，索－梁耦合結(jié)構(gòu)振動方程組及邊界條件無量綱化后表達形式［17］為

式中:下標(biāo)1，2分別表示梁、索;ρ，χ，θ分別為索與梁的質(zhì)量比、剛度比及傾角。

考慮梁、索橫向位移關(guān)系，表達式為

式中:φ1(x)，φ2(x)分別為梁、索模態(tài)，形式為

利用Galerkin方法對索－梁耦合結(jié)構(gòu)運動微分方程進行一階模態(tài)處理，獲得內(nèi)、外共振聯(lián)合激勵下系統(tǒng)二自由度非線性常微分方程為

式中:fij為激勵幅值;ai，bi，ci，aij，bij，aijk，bijk為梁、索Galerkin截斷系數(shù)。

由式(5)知，由于梁、索模態(tài)相互耦合，即使考慮梁模型是線性的，仍有非線性項(平方項、立方項)存在于梁的振動方程中;而索的振動方程中非線性項較僅考慮索幾何非線性時多。

2 攝動分析

2.1 內(nèi)外共振模式分析

引入多尺度參數(shù)ε，式(5)變換為

設(shè)梁、索解的表達式為

將式(6)代入式(7)，合并ε同次項且令各項系數(shù)為 0，則有

ε0階:

ε1階:

式(8)解的復(fù)數(shù)形式為

式中:A1，A2為待定函數(shù)分別為A1，A2的復(fù)共軛形式。

將式(10)代入式(9)，得

式中:ω1，ω2為梁、索頻率;cc為式(11)中函數(shù)共扼項;NST為式(11)中函數(shù)非長期項。

由式(11)知，索－梁耦合結(jié)構(gòu)存在多種內(nèi)、外共振模式，如 ω1=2ω2，ω1= ω2，ω1=ω2/2 的內(nèi)共振模式;Ω= ω1，Ω =2ω1，…，Ω = ω2，Ω =2ω2，…等外共振模式。通常結(jié)構(gòu)第一階模態(tài)占據(jù)振動的主要響應(yīng)，且模態(tài)階數(shù)越高對結(jié)構(gòu)振動響應(yīng)貢獻越小。因此，本文僅考慮梁、索各自第一階模態(tài)研究索－梁耦合非線性特性。

2.2 振動方程近似解與運動穩(wěn)定性

2.2.1 外激勵與梁主共振

考慮外激勵與梁主共振時，內(nèi)、外共振模式關(guān)系為

式中:σ1為調(diào)頻參數(shù);ε為遠小于1的參數(shù)。

將式(12)代入式(11)，并令長期項等于零，得

為求解式(13)，將A1，A2表示成極函數(shù)形式，即

將式(14)代入式(13)，并分離實、虛部，整理得

由式(15)知，考慮模態(tài)耦合影響的索－梁耦合結(jié)構(gòu)在不同σ1時表現(xiàn)出軟、硬化行為，而σ1值取決于系統(tǒng)的平方及立方非線性項，致耦合系統(tǒng)非線性特性更復(fù)雜。

2.2.2 外激勵與索主參數(shù)共振

考慮外激勵與索主參數(shù)共振時，梁由于1∶2內(nèi)共振機制會產(chǎn)生亞諧波共振。因此，內(nèi)、外共振模式關(guān)系式可表示為

將式(16)代入式(11)，并令長期項等于零，得復(fù)數(shù)形式的平均方程為

考慮式(17)，并分離實、虛部，整理得

由式(18)知，外激勵雖作用在梁上，但通過梁、索間內(nèi)共振機制作用已轉(zhuǎn)移到索的振動方程，梁此時相當(dāng)于被動激勵。因此，索的響應(yīng)較外激勵與梁主共振時大。而在內(nèi)外共振下，由于索的振動方程中出現(xiàn)阻尼項，若想激發(fā)外共振則需較大外界激勵幅值。

3 數(shù)值算例分析

由索－梁耦合結(jié)構(gòu)模型圖1，選具有5組不同系統(tǒng)參數(shù)耦合模型，模型參數(shù)見表1。通過每兩組模型對比分析系統(tǒng)參數(shù)對索－梁耦合結(jié)構(gòu)的非線性特性影響。

表1 不同系統(tǒng)參數(shù)的索－梁耦合結(jié)構(gòu)Tab.1 The different system parameters of the cable-beam coupled system

3.1 外激勵與梁主共振

考慮梁與索的內(nèi)共振關(guān)系為1∶2，研究外激勵直接作用于梁上且與梁產(chǎn)生1∶1主共振時，系統(tǒng)不同參數(shù)對幅頻響應(yīng)曲線影響。激勵幅值f11=0.005時，系統(tǒng)不同參數(shù)對梁、索各自幅頻響應(yīng)曲線影響見圖2。由圖2看出，梁、索響應(yīng)在σ1的某些范圍內(nèi)有多個解，不同系統(tǒng)初始條件均可獲得，但只能是其中的兩個(圖2(a))，雖本文認(rèn)為梁模型為線性，但其幅頻響應(yīng)曲線卻表現(xiàn)出非線性特征，即隨激勵頻率變化不同參數(shù)系統(tǒng)分別呈剛度軟(系統(tǒng)A、系統(tǒng)B、系統(tǒng)C)、硬化特征(系統(tǒng)D、系統(tǒng)E);而梁的剛度存在從軟化到硬化的轉(zhuǎn)變現(xiàn)象，盡管該現(xiàn)象僅在系統(tǒng)剛度比參數(shù)增加時發(fā)生(如系統(tǒng)D、系統(tǒng)E)。值得注意的是，系統(tǒng)剛度比參數(shù)增加使梁的幅頻曲線呈硬化特征時，剛度比繼續(xù)增加僅對振幅響應(yīng)有輕微影響。

圖2 索－梁耦合結(jié)構(gòu)不同系統(tǒng)參數(shù)對幅頻響應(yīng)曲線影響Fig.2 Effect of the different system parameters on the frequency-responses of the cable-beam coupled system

外激勵幅值對索－梁耦合結(jié)構(gòu)系統(tǒng)幅頻響應(yīng)曲線影響見圖3、圖4。其中，圖3以剛度軟化系統(tǒng)B為研究對象，圖4則以剛度硬化系統(tǒng)D為研究對象。與線性結(jié)構(gòu)不同，外激勵幅值大小只改變系統(tǒng)振動幅值，不會引起跳躍;而對非線性結(jié)構(gòu)，當(dāng)激勵幅值達到一定程度時，系統(tǒng)非線性被激發(fā)，呈明顯跳躍現(xiàn)象。由圖3、圖4可知，無論系統(tǒng)處于剛度軟化狀態(tài)或剛度硬化狀態(tài)，外激勵幅值增大時，系統(tǒng)非線性行為均更顯著，且振幅響應(yīng)均有明顯增加。

圖3 激勵幅值對系統(tǒng)B幅頻響應(yīng)曲線影響Fig.3 Effect of the amplitude of excitation on the frequency-responses for the system B

圖4 激勵幅值對系統(tǒng)D的幅頻響應(yīng)曲線影響Fig.4 Effect of the amplitude of excitation on the frequency-responses for the system D

圖5 索－梁耦合結(jié)構(gòu)不同系統(tǒng)參數(shù)對幅頻響應(yīng)曲線影響Fig.5 Effect of the different system parameters on the frequency-responses of the cable-beam coupled system

3.2 外激勵與索主參數(shù)共振

外激勵頻率Ω=2即外激勵與索產(chǎn)生1∶1主參數(shù)共振時，系統(tǒng)不同參數(shù)對幅頻響應(yīng)曲線影響見圖5～圖7。

圖6 外激勵幅值對系統(tǒng)B幅頻響應(yīng)曲線影響Fig.6 Effect of the amplitude of excitation on the frequency-responses for the system B

圖7 外激勵幅值對系統(tǒng)D激勵頻率與響應(yīng)幅值關(guān)系影響Fig.7 Effect of the amplitude of excitation on the frequency-responses for the system D

由圖5(a)知，梁的幅頻響應(yīng)曲線總處于剛度軟化狀態(tài)，且系統(tǒng)質(zhì)量比及索的軸壓比與等參數(shù)對梁幅頻響應(yīng)曲線均影響明顯，尤其質(zhì)量比參數(shù)增加最顯著:質(zhì)量比參數(shù)增加不僅使響應(yīng)振幅增大，且使幅頻曲線整體向σ2軸負方向移動(如系統(tǒng)B與系統(tǒng)C對比)。尤其當(dāng)系統(tǒng)剛度比較大時(如系統(tǒng)D、系統(tǒng)E)，梁的幅頻響應(yīng)完全消失，說明梁、索作為整體結(jié)構(gòu)應(yīng)用時，外激勵與索作主參數(shù)共振情況僅在系統(tǒng)剛度比參數(shù)較小時發(fā)生。由圖5(b)知，索的幅頻響應(yīng)曲線仍有兩個尖峰，但不同于圖2(b)內(nèi)容，無論系統(tǒng)參數(shù)如何變化兩個尖峰均保持同一趨勢，說明索的振動行為仍受梁振動影響，但其影響程度較弱，使索的振動行為支配系統(tǒng)振動行為。由圖5看出，系統(tǒng)的質(zhì)量比、索的垂跨比及軸壓比等參數(shù)對索的幅頻響應(yīng)曲線影響類似梁的影響，區(qū)別在于，增加質(zhì)量比參數(shù)會使索能在更大共振頻率范圍內(nèi)產(chǎn)生非線性行為，如系統(tǒng)B與系統(tǒng)C的對比。此外，對系統(tǒng)剛度比參數(shù)影響，索的幅頻曲線隨剛度比增加出現(xiàn)從剛度硬化到軟化的轉(zhuǎn)變現(xiàn)象，較圖2(b)內(nèi)容，可清楚發(fā)現(xiàn)兩者的剛度轉(zhuǎn)變現(xiàn)象完全相反，后者隨剛度比增加較大尖峰發(fā)生從剛度軟化到硬化狀態(tài)轉(zhuǎn)變，如系統(tǒng)B與系統(tǒng)D的比較。

基于上述現(xiàn)象，以索的剛度變化趨勢為指向，外激勵幅值對索－梁耦合結(jié)構(gòu)系統(tǒng)幅頻響應(yīng)曲線影響見圖6、圖7。其中，圖6以剛度硬化系統(tǒng)B為研究對象，圖7以剛度軟化系統(tǒng)D為研究對象。由圖6、圖7知，無論系統(tǒng)中索的剛度處于軟化狀態(tài)或硬化狀態(tài)，當(dāng)外激勵幅值增大時，系統(tǒng)中梁與索的非線性行為均趨向于更顯著，且振幅響應(yīng)均有明顯增加。對比圖4、圖6知，索－梁耦合結(jié)構(gòu)在外激勵與索作主參數(shù)共振時，欲使振幅響應(yīng)達到與外激勵與梁主共振級別，激勵幅值需增大104倍。

4 結(jié)論

研究內(nèi)、外共振聯(lián)合激勵下索－梁耦合結(jié)構(gòu)的非線性特性，利用多尺度法分析系統(tǒng)可能存在的多種內(nèi)、外共振模式。考慮索、梁間模態(tài)具有1:2內(nèi)共振關(guān)系下，對外激勵與梁發(fā)生主共振及與索發(fā)生主參數(shù)振動時的對應(yīng)工況分別進行攝動分析，獲得系統(tǒng)四維極坐標(biāo)形式的平均方程，結(jié)論如下:

(1)荷載僅作用在梁上時，索－梁耦合結(jié)構(gòu)由于模態(tài)耦合影響，存在兩種外共振機制，即荷載與梁外共振及荷載與索外共振，兩種外共振均能使系統(tǒng)振動呈非線性行為，后者所需激勵幅值較前者大104倍。

(2)無論荷載與梁發(fā)生主共振或與索發(fā)生主參數(shù)共振，即使考慮梁為線性模型，由于梁與索相互耦合振動影響，梁仍表現(xiàn)出多解、不穩(wěn)定及跳躍等非線性行為;盡管只考慮索的一個模態(tài)，但其非線性行為呈現(xiàn)雙模態(tài)特征，即兩個幅頻曲線。對前者外共振，索的最大響應(yīng)總與梁的響應(yīng)保持一致;而對后者外共振，該一致性消失，且隨激勵幅值增加非線性行為愈顯著。

(3)荷載與梁作主共振時，索的垂跨比及系統(tǒng)質(zhì)量比對梁的非線性影響微小，但卻能顯著影響索的非線性行為;系統(tǒng)剛度比增加會使梁發(fā)生從剛度軟化到硬化狀態(tài)轉(zhuǎn)變。此時索的非線性行為將產(chǎn)生兩次跳躍現(xiàn)象;荷載與索發(fā)生主參數(shù)共振時系統(tǒng)質(zhì)量比對索的影響較索的垂跨比影響更顯著;系統(tǒng)剛度比較大時梁的非線性完全消失，索發(fā)生從剛度硬化到軟化轉(zhuǎn)變。

(4)鑒于剛度比參數(shù)對索－梁耦合結(jié)構(gòu)有重要影響，建議對索－梁耦合結(jié)構(gòu)進行振動控制時應(yīng)避免用增加剛度策略或應(yīng)盡量減小因控制措施造成的系統(tǒng)剛度變化，否則將有可能引起結(jié)構(gòu)更復(fù)雜的振動行為。

［1］Nayfeh A H，Arafat H N，Chin C M，et al.Multi-mode interactions in suspended cables［J］.Journal of Sound and Vibration，2002，8:337 －387.

［2］王波，徐豐，張海龍.端部激勵下空間傾斜拉索非線性振動特性研究［J］.振動與沖擊，2009，28(5):172－175.WANG Bo， XU Feng， ZHANG Hai-long. Non-linear vibration characteristics of spatial inclined cables under periodic support excitation［J］.Journal of Vibration and Shock，2009，28(5):172 －175.

［3］吳曉，黎大志，羅佑新.斜拉索非線性固有振動特性分析［J］.振動與沖擊，2003，22(3):34－39.WU Xiao， LI Da-zhi， LUO You-xin. Nonlinear natural vibration character analysis of stay cables［J］.Journal of Vibration and Shock，2003，22(3):34 －39.

［4］Ghayesh M H. Subharmonic dynamics of an axially accelerating beam［J］.Archive of Applied Mechanics，2012，82:1169－1181.

［5］羅帥，劉紅軍，王剛.考慮橋面運動的斜拉索減振模型［J］.深圳大學(xué)學(xué)報:理工版，2010，27(4):471 －473.LUO Shuai，LIU Hong-jun，WANG Gang.The characteristic of stay-cable damping system in consideration of bridge deck vibration［J］.Journal of Shenzhen University Science and Engineering，2010，27(4):471 －473.

［6］Fujino Y，Warnitchai P，Pacheco B M.An experimental and analytical study of autoparametric resonance in a 3DOF model of cable-stayed-beam［J］.Nonlinear Dynamics，1993，4:111－138.

［7］Xia Y，F(xiàn)ujino Y.Auto-parametric vibration of a cable-stayedbeam structure under random excitation［J］.Journal of Engineering Mechanics-Asce，2006，132:279 －286.

［8］Fung R F，Lu L Y，Huang S C.Dynamic modelling and vibration analysis of a flexible cable-stayed beam structure［J］.Journal of Sound and Vibration，2002，254(4):717－726.

［9］Gattulli V，Morandini M，Paolone A.A parametric analytical model for non-linear dynamics in cable-stayed beam［J］.Earthquake Engineering ＆ Structural Dynamics，2002，31:1281－1300.

［10］Gattulli V，Lepidi M.Nonlinear interactions in the planar dynamics of cable-stayed beam［J］.International Journal of Solids and Structures，2003，40:4729 －4748.

［11］Gattulli V，Lepidi M，Mjh G，et al.One-to-two global-local interaction in a cable-stayed beam observed through analytical，finite element and experimental models［J］.International Journal of Non-Linear Mechanics，2005，40:571－588.

［12］王濤，沈銳利，李洪.斜拉橋索－梁相關(guān)振動概念及其研究方法初探［J］.振動與沖擊，2013，32(20):29 －34.WANG Tao，SHEN Rui-li，LI Hong.Primary exploration for concept and studying method of cable-beam vibration in a cable-stayed bridge［J］.Journal of Vibration and Shock，2013，32(20):29－34.

［13］趙躍宇，蔣麗忠，王連華，等.索－梁組合結(jié)構(gòu)的動力學(xué)建模理論及其內(nèi)共振分析［J］.土木工程學(xué)報，2004，37:69－72.ZHAO Yue-yu，JIANG Li-zhong，WANG Lian-hua，et al.The dynamical modeling theory and internal resonance of cable-beam composite structure［J］.China Civil Engineering Journal，2004，37:69 －72.

［14］趙躍宇，楊相展，劉偉長，等.索－梁組合結(jié)構(gòu)中拉索的非線性響應(yīng)［J］.工程力學(xué)，2006，23(11):153－158.ZHAO Yue-yu，YANG Xiang-zhan，LIU Wei-chang，et al.Nonlinear response of cables in cable-stayed beam structure［J］.Engineering Mechanics，2006，23(11):153 －158.

［15］馮維明，高黎黎，金棟平.索－梁耦合系統(tǒng)解的穩(wěn)定性分析［J］.應(yīng)用力學(xué)學(xué)報，2008，25(2):284－288.FENG Wei-ming， GAO Li-li， JIN Dong-ping. Nonlinear dynamic analysis for coupled structure of cable-stayed beam［J］.Chinese Journal of Applied Mechanics，2008，25(2):284－288.

［16］馮維明，高黎黎.索－梁耦合系統(tǒng)非線性振動分析［J］.振動工程學(xué)報，2008，21(2):115－119.FENG Wei-ming，GAO Li-li.Nonlinear vibration analysis for coupled structure of cable-stayed beam［J］.Journal of Vibration Engineering，2008，21(2):115 －119.

［17］Wei M H，Xiao Y Q，Liu H T.Bifurcation and chaos of a cable-beam coupled system under simultaneous internal and external resonances［J］.Nonlinear Dynamics，2012，67(3):1969－1984.