支座激勵(lì)下小垂度懸索動(dòng)力響應(yīng)試驗(yàn)及數(shù)值模擬

王國(guó)威，張 鵬，盧文勝*

(1.同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室，上海 200092；2.江西理工大學(xué) 土木與測(cè)繪工程學(xué)院，江西 贛州 341000)

索網(wǎng)幕墻采光頂是近年來(lái)興起的一種建筑外圍護(hù)體系，因其美觀輕盈、通透性能優(yōu)越，被廣泛應(yīng)用于多功能商業(yè)中心、體育場(chǎng)館和機(jī)場(chǎng)航站樓等大型公共設(shè)施，如北京金成大廈[1]、聯(lián)想融科資訊中心[2]和德國(guó)海爾布隆車站[3]等。其中，懸索作為此類結(jié)構(gòu)的基本承力單元，具有剛度低、非線性強(qiáng)的特點(diǎn)，在風(fēng)、地震及環(huán)境振動(dòng)等諸多不利因素作用下，有可能發(fā)生大變形和索力突變，尤其是連接于兩主體結(jié)構(gòu)的采光頂懸索，在強(qiáng)震作用下，主體結(jié)構(gòu)若發(fā)生相對(duì)位移振動(dòng)，懸索會(huì)發(fā)生過(guò)度張拉或松弛現(xiàn)象，甚至?xí)霈F(xiàn)斷絲破壞。因此，懸索的靜動(dòng)力特性和動(dòng)力響應(yīng)研究受到廣泛關(guān)注。

國(guó)內(nèi)外學(xué)者對(duì)懸索的靜動(dòng)力特性開(kāi)展了深入的研究。Irvine等[4]根據(jù)線性理論，建立小垂跨比（小于1∶8）懸索的自由振動(dòng)方程，得到綜合反映懸索幾何形狀及材性的Irvine系數(shù)。Rega等[5]進(jìn)一步考慮懸索的強(qiáng)非線性效應(yīng)，建立了彈性懸索的大振幅自由振動(dòng)偏微分方程，研究其垂跨比在較大范圍內(nèi)的非線性現(xiàn)象，并指出懸索的軟、硬化行為取決于自身特性和振幅。沈世釗等[6]系統(tǒng)地闡述了懸索的靜力解析和有限元計(jì)算方法，并通過(guò)部分算例驗(yàn)證方法的有效性。為考慮不同邊界條件的影響，劉海濤等[7]采用多尺度法研究索–梁耦合結(jié)構(gòu)的非線性行為，分析得出主共振下兩者剛度比較大時(shí)，索從剛度硬化轉(zhuǎn)變?yōu)檐浕?，而梁的非線性行為消失；Zhang等[8–9]對(duì)比分析端部固定和增設(shè)彈簧的懸索結(jié)構(gòu)的靜動(dòng)力特性，計(jì)算表明，邊界約束降低會(huì)導(dǎo)致懸索變形增加，索力降低。

在索網(wǎng)幕墻采光頂?shù)葘?shí)際工程中，支承懸索的主結(jié)構(gòu)質(zhì)量遠(yuǎn)大于索結(jié)構(gòu)本身，故主結(jié)構(gòu)的運(yùn)動(dòng)也被視為支座激勵(lì)施加于懸索端部。Perkins[10]通過(guò)對(duì)懸索平面內(nèi)、外耦合模態(tài)的近似描述，檢驗(yàn)了所提出的兩自由度懸索非線性振動(dòng)方程的正確性。Benedettini等[11]基于連續(xù)化理論，建立支座激勵(lì)下懸索的四自由度模型，并討論了幾何性質(zhì)等參數(shù)對(duì)系統(tǒng)動(dòng)力響應(yīng)的影響。Guo等[12]考慮幾何非線性和索–支承耦合效應(yīng)建立懸索振動(dòng)簡(jiǎn)化模型，并通過(guò)數(shù)值方法進(jìn)行驗(yàn)證。此外，在索結(jié)構(gòu)振動(dòng)測(cè)試方面，胡俊等[13]實(shí)測(cè)了某懸索橋的風(fēng)振響應(yīng)，監(jiān)測(cè)數(shù)據(jù)表明，日常風(fēng)環(huán)境下吊索為高頻渦激振動(dòng)。賀文宇等[14]以兩端固定的斜索為例，通過(guò)有限差分法計(jì)算索的模態(tài)參數(shù)，探究了斜拉索傾角對(duì)索力識(shí)別精度的影響。桂成中等[15]提出一種索力增量校驗(yàn)系數(shù)的方法，用于評(píng)定拉索和吊桿的技術(shù)狀況。懸索在隨機(jī)荷載下的響應(yīng)研究則較少，Georgakis等[16]基于有限元法分析了斜拉索在豎向激勵(lì)下的面內(nèi)、外振動(dòng)位移。任淑琰等[17]采用解析方法對(duì)斜拉索在軸向隨機(jī)激勵(lì)下的參數(shù)振動(dòng)進(jìn)行分析。王波等[18]研究了斜拉索在橫橋向激勵(lì)下面內(nèi)的耦合振動(dòng)特性，結(jié)果表明，在一定的橫向激勵(lì)下，拉索由于耦合作用其面內(nèi)振動(dòng)將被激發(fā)。

懸索作為采光頂?shù)闹饕芰?gòu)件因支撐玻璃面板具有一定的垂度，主體結(jié)構(gòu)的相對(duì)位移振動(dòng)將作為支座激勵(lì)導(dǎo)致懸索構(gòu)件產(chǎn)生不利的變形。因此，本文主要研究支座激勵(lì)下具有一定垂度懸索的動(dòng)力特性和動(dòng)力響應(yīng)。綜上，現(xiàn)有研究大多集中在理論分析層面，從試驗(yàn)角度對(duì)懸索在支座激勵(lì)下的動(dòng)力響應(yīng)研究不足，且多數(shù)文獻(xiàn)尚未考慮支座隨機(jī)激勵(lì)及垂度對(duì)懸索動(dòng)力特性和響應(yīng)的影響。為進(jìn)一步開(kāi)展支座簡(jiǎn)諧激勵(lì)下懸索動(dòng)力特性、響應(yīng)的試驗(yàn)和數(shù)值分析，對(duì)比了兩者結(jié)果驗(yàn)證數(shù)值模型的正確性，并對(duì)支座隨機(jī)激勵(lì)下懸索的動(dòng)力響應(yīng)進(jìn)行了參數(shù)研究，對(duì)相鄰結(jié)構(gòu)間索網(wǎng)幕墻采光頂?shù)脑O(shè)計(jì)與應(yīng)用具有一定的工程參考價(jià)值。

1 支座簡(jiǎn)諧激勵(lì)下懸索動(dòng)力特性及響應(yīng)試驗(yàn)

1.1 試驗(yàn)設(shè)計(jì)與加載

為獲取懸索在支座激勵(lì)下的動(dòng)力特性和動(dòng)力響應(yīng)，在同濟(jì)大學(xué)土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室設(shè)計(jì)并搭建試驗(yàn)平臺(tái)，包括試驗(yàn)?zāi)Ｐ?、加載裝置及力與位移傳感器等。試件采用直徑為10 mm的1×19鍍鋅鋼絞線，跨度為7.3 m，最小破斷力為69.5 kN，其中單根鋼絲強(qiáng)度為1 320 MPa。將懸索試件一端通過(guò)耳板鉸接于反力墻上，另一端同樣與加載裝置鉸接連接，且在索上附加若干數(shù)量的配重塊，由MTS作動(dòng)器直接通過(guò)加載裝置對(duì)懸索施加水平支座激勵(lì)。

首先，建立小垂度懸索面內(nèi)的無(wú)阻尼自由振動(dòng)頻率方程[19]，其歸一化后的表達(dá)式為：

式中： ωn為懸索的正對(duì)稱模態(tài)頻率；n為頻率階數(shù)；α=E·A/H為懸索軸向剛度與索力比，H為初始索力，E和A分 別為彈性模量和截面面積；b=w/l為懸索跨中垂度（w）與跨度（l）之比。

由表達(dá)式（1）可知， αb2可以綜合考慮垂跨比、索力及材料常數(shù)對(duì)懸索振動(dòng)特性的影響。因此，下文將其作為關(guān)鍵參數(shù)研究懸索在支座激勵(lì)下的動(dòng)力特性和動(dòng)力響應(yīng)。

為研究變形特征參數(shù) αb2的影響，設(shè)計(jì)3組懸索試件參數(shù)（ αb2= 0.04、0.08、0.19），保持跨度不變，通過(guò)調(diào)整索力及配重?cái)?shù)改變垂跨比。懸索Ⅰ、Ⅱ等間隔布置8個(gè)配重，懸索Ⅲ等間隔布置17個(gè)配重，單個(gè)配重為24 kg，具體懸索試件參數(shù)見(jiàn)表1。

表1 懸索試件參數(shù)Tab.1 Parameters of the small-sagged cable

支座激勵(lì)方式為簡(jiǎn)諧位移荷載，荷載工況分變幅值和變頻率加載，考慮MTS作動(dòng)器的加載能力，兩組加載工況列于表2。試驗(yàn)的整個(gè)加載過(guò)程為：首先，通過(guò)懸索端部張緊裝置施加1 kN的初始拉力，固定限位構(gòu)件以防止安裝配重過(guò)程中懸索發(fā)生松弛；其次，釋放限位約束，檢查各傳感器及采集儀工作狀態(tài)，并通過(guò)作動(dòng)器控制索力在預(yù)定值附近；最后，采用正弦波位移激勵(lì)：Xd(t)=Dsin(2πft)，使MTS作動(dòng)器連續(xù)往復(fù)加載30周，兩組工況的總加載時(shí)間t分別為25和60 s，加載完畢后讓懸索自由振動(dòng)，直至振動(dòng)停止，測(cè)量懸索的頻率和阻尼比。

表2 試驗(yàn)加載工況Tab.2 Cases of the testing

試驗(yàn)采用力傳感器和位移計(jì)2種傳感器。在懸索兩端分別安裝一個(gè)量程為5 t的力傳感器，實(shí)時(shí)量測(cè)索力大??；在懸索跨度的1/3、1/2和2/3位置處各布置一個(gè)量程為750 mm的位移計(jì)，用來(lái)測(cè)量懸索面內(nèi)的豎向動(dòng)位移；力傳感器和位移計(jì)的采樣頻率均設(shè)定為50 Hz，并通過(guò)數(shù)據(jù)采集儀記錄各傳感器數(shù)據(jù)。在懸索下方搭設(shè)防護(hù)架以防止配重塊墜落損壞試驗(yàn)設(shè)備，試驗(yàn)測(cè)點(diǎn)及裝置示意圖和加載裝置現(xiàn)場(chǎng)照片如圖1、2所示。

圖1 試驗(yàn)測(cè)點(diǎn)及裝置示意圖Fig.1 Schematic diagram of test measuring points and device

圖2 懸索試驗(yàn)裝置Fig.2 Testing setup of the small-sagged cable specimen

1.2 試驗(yàn)結(jié)果與分析

1.2.1 懸索動(dòng)力特性

圖3為各懸索在自由振動(dòng)階段由位移計(jì)實(shí)測(cè)數(shù)據(jù)計(jì)算所得的基頻，從圖3中可以看出：隨著 αb2的增加，懸索基頻 ω1從2.07 Hz降低至1.90 Hz。在邊界條件相同的情況下，懸索的剛度隨垂跨比的增加而降低。在支座簡(jiǎn)諧激勵(lì)完畢后，以懸索Ⅰ （ αb2= 0.04）和懸索Ⅲ（ αb2= 0.19）為例，提取跨中自由振動(dòng)曲線并計(jì)算懸索等效黏滯阻尼比ξ[20]，自由衰減法計(jì)算公式如下：

圖3 不同 α b2下 懸索的基頻ω1Fig.3 Base frequencies of the small-sagged cable with differentαb2

式中，vr和vr+s分別為r和r+s時(shí)刻出現(xiàn)的相鄰正波峰。

圖4為懸索的自由振動(dòng)曲線，隨著變形特征參數(shù)αb2的增加，自由振動(dòng)曲線的指數(shù)分布規(guī)律特征逐漸明顯，通過(guò)指數(shù)擬合得到懸索Ⅰ與Ⅲ的阻尼比分別為1.2%、0.8%。

圖4 懸索Ⅰ、Ⅲ的自由振動(dòng)曲線Fig.4 Free vibration of the small-sagged Cable Ⅰ and Ⅲ

1.2.2 懸索動(dòng)力響應(yīng)

以懸索跨度1/2位置處為例，選取懸索Ⅰ和Ⅲ為研究對(duì)象，分析 αb2對(duì)支座激勵(lì)下懸索動(dòng)力響應(yīng)及動(dòng)力特性的影響，獲得動(dòng)力響應(yīng)時(shí)程曲線，如圖5所示。為便于對(duì)照，工況1和2均截取前25 s的時(shí)程曲線。由圖5（a）可知：在工況1支座激勵(lì)下，懸索面內(nèi)動(dòng)位移幅值 j?wmaxj隨 αb2的增加而增大， j ?wmaxj增加27.7 %。相反地，動(dòng)索力幅值 j?Hmaxj隨 αb2的增加而明顯減小，j?Hmaxj降 低69.1 %；當(dāng) αb2= 0.19時(shí)，除0 ～ 5 s時(shí)間段內(nèi)的索力變化較大外，5 ～ 25 s時(shí)間段內(nèi)， j?Hmaxj均在1.5 kN附近波動(dòng)，分析主要原因是：懸索的初始垂度較大，施加相同的支座激勵(lì)，索力變化較小。由圖5（b）可知，在工況2支座激勵(lì)下，面內(nèi) j?wmaxj和j?Hmaxj的變化也表現(xiàn)出相同的規(guī)律，動(dòng)位移幅值略有增大，但 j ?Hmaxj降低29.4 %。

圖5 不同α b2下懸索動(dòng)力響應(yīng)時(shí)程曲線Fig.5 Dynamic response curves of the small-sagged cable with differentαb2

2 數(shù)值模擬

2.1 數(shù)值模型建立

采用ANSYS軟件建立懸索數(shù)值模型，模型尺寸（l×A）、配重塊數(shù)目及彈性模量E均與試驗(yàn)構(gòu)件相同，泊松比μ=0.3，密度取鋼材材料密度ρ=7 800 kg/m3。懸索截面較小，忽略自身抗彎剛度對(duì)面內(nèi)振動(dòng)的影響，采用Link 10僅受拉單元模擬，通過(guò)初位移 ?L施加初始索力，配重塊采用Mass 21結(jié)構(gòu)質(zhì)量單元模擬，不考慮轉(zhuǎn)動(dòng)慣量效應(yīng)。懸索模型一端（A端）為鉸接約束，另一端（B端）釋放X方向約束，由于主要研究懸索面內(nèi)（X、Y向）的振動(dòng)，故同時(shí)約束索體Z向平動(dòng)位移。采用絕對(duì)位移加載，在支座位置處施加與試驗(yàn)相同的正弦波位移激勵(lì)Xd(t)，計(jì)算模型如圖6所示。

圖6 支座激勵(lì)下懸索振動(dòng)計(jì)算模型Fig.6 Computing model of a small-sagged cable under support excitations

具體建模過(guò)程如下：

1）根據(jù)試驗(yàn)構(gòu)件建立幾何模型，定義并賦予各部分材料屬性，將模型劃分為50個(gè)單元，并施加初始索力及設(shè)置模型邊界條件。

2）進(jìn)行懸索靜力計(jì)算，打開(kāi)“NLGEOM”命令，考慮重力和幾何非線性影響，對(duì)懸索進(jìn)行找形分析，獲得靜載下懸索 αb2的試驗(yàn)和有限元對(duì)比值。

3）修正模型坐標(biāo)，定義計(jì)算和擴(kuò)展模態(tài)數(shù)，計(jì)算懸索的1階頻率，并與試驗(yàn)實(shí)測(cè)結(jié)果對(duì)比。

4）采用完全法瞬態(tài)動(dòng)力分析，假定各階模態(tài)阻尼比相同，將實(shí)測(cè)阻尼比 ξ轉(zhuǎn)化為Rayleigh阻尼，并通過(guò)命令“ALPHAD”和“BETAD”定義阻尼系數(shù)；時(shí)間步長(zhǎng) ?t=0.05 s，采用水平支座簡(jiǎn)諧激勵(lì)進(jìn)行模擬試驗(yàn)加載；最后，提取懸索響應(yīng)結(jié)果曲線并與試驗(yàn)結(jié)果對(duì)比分析。

2.2 數(shù)值模擬與試驗(yàn)結(jié)果對(duì)比驗(yàn)證

表3給出了3組懸索試件靜載下索力H、撓度b、αb2及1階頻率的對(duì)比結(jié)果。由表3可知，有限元計(jì)算結(jié)果與試驗(yàn)結(jié)果非常接近，對(duì)于懸索Ⅰ、Ⅱ和Ⅲ，各參數(shù)誤差均在5%～12%以內(nèi)。

表3 靜力及頻率計(jì)算對(duì)比結(jié)果Tab.3 Comparison results of static and frequency calculation

以懸索Ⅱ跨度1/2位置為例，圖7給出了工況1和2下動(dòng)力響應(yīng)試驗(yàn)結(jié)果和數(shù)值模擬結(jié)果的對(duì)比，由圖7可知：1）懸索Ⅱ在工況1下的動(dòng)位移幅值 j?wmaxj和索力幅值 j?Hmaxj變 化較工況2更為明顯， j?wmaxj和 j?Hmaxj分別為44.37 mm和4.02 kN，數(shù)值模擬和試驗(yàn)結(jié)果在部分峰值處略有差異，但平均索力和平均位移變化幅值誤差均在12 %以內(nèi)，主要原因是：在高頻的支座激勵(lì)下，懸索上的配重塊剛心與質(zhì)心不完全重合，配重發(fā)生少量偏轉(zhuǎn)，對(duì)剛度有一定的影響，從而導(dǎo)致計(jì)算與試驗(yàn)存在誤差，但兩者結(jié)果的幅值變化趨勢(shì)和相位均保持同步。2）相比于工況1，工況2的位移和索力變化相對(duì)平穩(wěn)，且數(shù)值模擬和試驗(yàn)結(jié)果吻合較好，數(shù)值模擬和試驗(yàn)的 j ?wmaxj、 j ?Hmaxj的誤差分別為5.4%、6.8%。綜合上述結(jié)果可知，懸索動(dòng)力特性及響應(yīng)試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果均吻合較好，驗(yàn)證了有限元模型的準(zhǔn)確性。

圖7 懸索Ⅱ（= 0.08）試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果對(duì)比Fig.7 Comparison of test results and numerical simulations results for Cable Ⅱ(α b2 = 0.08)

在前文驗(yàn)證有限元建模方法正確的前提下，以懸索Ⅰ（ αb2=0.04）的有限元模型為例，令支座簡(jiǎn)諧激勵(lì)頻率f=2.07 Hz，研究不同阻尼比 ξ下懸索的1階共振響應(yīng)，總加載時(shí)間為15 s，獲得共振響應(yīng)時(shí)程曲線如圖8所示。從圖8可以看出：當(dāng)阻尼比 ξ=0.01時(shí)，動(dòng)位移和動(dòng)索力呈周期性的“震蕩”變化，但其幅值 j?wmaxj和 j?Hmaxj總體呈下降的趨勢(shì)；當(dāng)阻尼比 ξ= 0.20時(shí)，位移和索力變化均呈現(xiàn)平穩(wěn)變化，j?wmaxj、 j?Hmaxj分別穩(wěn)定在20 mm、4 kN附近。值得注意的是，隨著阻尼比的增加，懸索1階共振下位移和索力的響應(yīng)均受到不同程度的控制，位移控制較索力下降明顯， j?wmaxj降 低64.0 %， j?Hmaxj降低33.2 %。從側(cè)面反映出，若在工程中的索網(wǎng)采光頂懸索端部或連接件部位布置黏滯、黏彈或摩擦阻尼器等裝置，將會(huì)有效地降低懸索的位移和索力響應(yīng)。

圖8 不同阻尼比 ξ下懸索1階共振響應(yīng)時(shí)程曲線（α b2=0.04）Fig.8 First-order resonance responses of the small-sagged cable with different ξ (α b2=0.04)

3 支座隨機(jī)激勵(lì)下懸索動(dòng)力響應(yīng)分析

風(fēng)、地震及環(huán)境振動(dòng)等屬于非確定性作用，均具有一定的隨機(jī)性，因此，研究懸索在隨機(jī)激勵(lì)下的動(dòng)力響應(yīng)具有重要的現(xiàn)實(shí)意義。為討論懸索在支座隨機(jī)激勵(lì)下的動(dòng)力響應(yīng)，在前文工作的基礎(chǔ)上，以理想高斯白噪聲模擬隨機(jī)環(huán)境激勵(lì)，分析懸索控制參數(shù) αb2、阻尼比 ξ和激勵(lì)強(qiáng)度S0對(duì)懸索隨機(jī)動(dòng)力響應(yīng)的影響。

基于有限元法將懸索離散，對(duì)于阻尼較小的懸索結(jié)構(gòu)，由振型j產(chǎn)生的響應(yīng)與振型k的響應(yīng)近似滿足統(tǒng)計(jì)獨(dú)立[6]，則位移響應(yīng)譜密度Sui(ω)可表示為：

式中，SF j(ω)為 激勵(lì)的譜密度函數(shù)， ω為激勵(lì)圓頻率，ωj為 懸索第j振型圓頻率， ?ij為第i自由度的第j振型值，J為振型總數(shù)，為傳遞函數(shù)的模。對(duì)式（3）進(jìn)行積分求和，可得位移響應(yīng)均方差：

隨機(jī)激勵(lì)下懸索第e個(gè)單元的索力響應(yīng) ?Te為：

采用均值為0、激勵(lì)強(qiáng)度為S0的理想高斯白噪聲D(t)作為支座水平激勵(lì)，其自相關(guān)函數(shù)可表示為：E[D(t)D(t+τ)]=2πS0δ(τ)， 其中， δ (τ)為Dirac–δ函數(shù)，τ為時(shí)間間隔，E[·]為 [·]的期望。以試驗(yàn)中懸索Ⅰ的幾何和材料參數(shù)為例，在ANSYS中建立支座隨機(jī)激勵(lì)下懸索的數(shù)值模型，并進(jìn)行隨機(jī)振動(dòng)計(jì)算（PSD功率譜密度分析）。主要建模過(guò)程如下：

1）懸索幾何建模、單元類型、單元?jiǎng)澐旨斑吔缂s束與簡(jiǎn)諧激勵(lì)相同，將初始索力H直接施加在釋放X方向自由度的索端支座處。

2）靜力計(jì)算完畢后，定義計(jì)算模態(tài)和擴(kuò)展模態(tài)，同時(shí)令“ELCALC = YES”計(jì)算應(yīng)力結(jié)果。

3）采用位移譜作為水平支座激勵(lì)，由于軟件輸入的要求，將雙邊譜轉(zhuǎn)換為單邊譜進(jìn)行加載，定義計(jì)算頻率范圍為0.1～10.0 Hz，激勵(lì)強(qiáng)度S0為0.01～0.20 m2/(rad·s?1)。

4）在PSD分析過(guò)程中，支座激勵(lì)的縮放系數(shù)取為1，采用自動(dòng)時(shí)間步，阻尼比 ξ通過(guò)命令“DMPRAT”施加，采用“PSDCOM”命令計(jì)算懸索位移及索力響應(yīng)的均方差（1σ），即位于–1σ～1σ區(qū)間發(fā)生概率為68.3 %的結(jié)果量。

3.1 參數(shù) αb2對(duì)懸索位移及索力響應(yīng)的影響

首先，研究參數(shù) αb2對(duì)支座隨機(jī)激勵(lì)下懸索動(dòng)力響應(yīng)的影響，計(jì)算不同 αb2下的懸索位移和索力響應(yīng)，結(jié)果如圖9所示。

圖9（a）、（b）為3組 αb2的懸索跨中1/2位置處的位移、索力響應(yīng)譜，可以發(fā)現(xiàn)，兩者結(jié)果符合正態(tài)分布。由于 αb2與結(jié)構(gòu)特性直接相關(guān)，隨著 αb2的增加，懸索動(dòng)力響應(yīng)譜幅值出現(xiàn)左移，說(shuō)明結(jié)構(gòu)頻率逐漸降低，這也與圖3中試驗(yàn)的實(shí)測(cè)結(jié)果相一致。在相同的激勵(lì)強(qiáng)度S0和阻尼比 ξ下，懸索動(dòng)力響應(yīng)隨 αb2的增加而顯著加大，位移響應(yīng)譜幅值增加了3倍，索力響應(yīng)譜幅值甚至增加了20倍。圖9（c）為懸索響應(yīng)的均方差與 αb2的 關(guān)系，可知： αb2從0.04增加至1.0的過(guò)程中，位移和索力的均方差呈上升趨勢(shì)，但增加的幅度隨 αb2的增加而逐漸降低；當(dāng) αb2= 1.0時(shí)，懸索的位移均方差為0.339 m，增加了3.3倍，索力均方差為2.1 kN，增加了3.17倍。該規(guī)律反映出 αb2較大的懸索，其位移和索力響應(yīng)對(duì)支座的隨機(jī)激勵(lì)更加敏感。

圖9 不同α b2下 懸索的位移、索力響應(yīng)（S0=0.01, ξ=0.05）Fig.9 Responses of the small-sagged cable with different αb2 (S0=0.01, ξ=0.05)

3.2 阻尼比ξ 對(duì)懸索位移及索力響應(yīng)的影響

取不同的懸索阻尼比ξ，分析其對(duì)懸索位移及索力響應(yīng)的影響。不同ξ 下懸索的動(dòng)力響應(yīng)譜如圖10（a）、（b）所示，可知：當(dāng)ξ=0.01時(shí)，位移和索力響應(yīng)譜幅值最大，分別為0.58 m2/Hz和11.55 kN2/Hz；當(dāng) ξ=0.05時(shí)，位移和索力響應(yīng)譜幅值已降低至0.02 m2/Hz和0.47 kN2/Hz。對(duì)比兩組結(jié)果，可以發(fā)現(xiàn)，懸索位移和索力響應(yīng)已被顯著抑制。從圖10（c）懸索響應(yīng)的均方差曲線也能看出，位移和索力均方差隨著阻尼比ξ的增加而降低，當(dāng)ξ= 0.20時(shí)，兩者分別減小68.1 %和75.3 %；另外，兩者曲線下降趨勢(shì)逐漸變緩，表明當(dāng)ξ＜0.10時(shí)，位移和索力響應(yīng)對(duì)阻尼比的變化敏感性較高，增加 ξ帶來(lái)的抑制效果顯著。上述規(guī)律與支座簡(jiǎn)諧激勵(lì)下懸索1階共振響應(yīng)隨阻尼比 ξ的變化相同，也進(jìn)一步說(shuō)明增加懸索系統(tǒng)阻尼將明顯降低位移和索力響應(yīng)。

圖10 不同 ξ下懸索的位移、索力響應(yīng)（S0=0.01, α b2=0.04）Fig.10 Responses of the small-sagged cable with different ξ (S0=0.01, α b2=0.04)

3.3 激勵(lì)強(qiáng)度S0對(duì)懸索位移及索力響應(yīng)的影響

若支座所受的隨機(jī)環(huán)境激勵(lì)強(qiáng)度S0發(fā)生改變，系統(tǒng)的動(dòng)力響應(yīng)也將隨之變化。計(jì)算不同激勵(lì)強(qiáng)度S0下的懸索動(dòng)力響應(yīng)，如圖11所示。

圖11 不同S0下懸索的位移、索力響應(yīng)（ ξ =0.05, α b2=0.04）Fig.11 Responses of the small-sagged cable with different S0 ( ξ =0.05, α b2=0.04)

圖11表 明，在 相 同 的 αb2和 阻 尼 比 ξ 下，、及位移、索力均方差均隨S0的增加而增大；對(duì)應(yīng)的增加了7.1倍，增加了6.9倍，響應(yīng)均方差曲線斜率逐漸減小，但總體呈上升趨勢(shì)。

4 結(jié) 論

本文首先進(jìn)行了懸索在支座簡(jiǎn)諧激勵(lì)下動(dòng)力特性及響應(yīng)的試驗(yàn)和數(shù)值模擬研究，對(duì)比兩者結(jié)果驗(yàn)證模型的正確性；其次，通過(guò)數(shù)值方法計(jì)算懸索在支座隨機(jī)荷載激勵(lì)下的動(dòng)力響應(yīng)，并分析響應(yīng)譜及均方差隨 αb2、 阻尼比 ξ和激勵(lì)強(qiáng)度S0關(guān)鍵參數(shù)的變化規(guī)律，得到的主要結(jié)論如下：

1）設(shè)計(jì)并開(kāi)展懸索在支座簡(jiǎn)諧激勵(lì)下的動(dòng)力響應(yīng)試驗(yàn)，分析不同 αb2下懸索動(dòng)力響應(yīng)的變化：在兩組工況下，懸索動(dòng)位移隨 αb2的增加而小幅增加；相反地，動(dòng)索力隨 αb2的增加而明顯增大。結(jié)果表明，相同的支座簡(jiǎn)諧激勵(lì)下，相比于位移響應(yīng)，索力響應(yīng)對(duì)懸索的關(guān)鍵參數(shù) αb2更加敏感。

2）建立懸索試件有限元模型，模型結(jié)果與懸索自由振動(dòng)試驗(yàn)實(shí)測(cè)頻率、靜載下實(shí)測(cè)垂跨比及各工況下懸索試件動(dòng)力響應(yīng)試驗(yàn)結(jié)果吻合較好，從而驗(yàn)證了有限元模型的正確性。

3）對(duì)支座簡(jiǎn)諧激勵(lì)下懸索的1階共振響應(yīng)進(jìn)行有限元計(jì)算，結(jié)果表明，懸索動(dòng)位移受阻尼比 ξ的影響更大，且時(shí)程響應(yīng)變化平穩(wěn)。

4）在支座隨機(jī)激勵(lì)下，通過(guò)數(shù)值方法討論了 αb2、ξ和S0對(duì)懸索位移及索力響應(yīng)的影響。計(jì)算表明，兩種響應(yīng)譜值和均方差均隨參數(shù) αb2、激勵(lì)強(qiáng)度S0的增加而增大，隨阻尼比 ξ的增加而減小，且變化幅度逐漸降低，說(shuō)明3種參數(shù)對(duì)懸索動(dòng)力響應(yīng)的影響會(huì)逐步減弱。