徐彥青 郭 彤

(東南大學(xué)土木工程學(xué)院,南京 210096)

斜拉索因其柔性大、阻尼小以及質(zhì)量輕的特點,在風(fēng)載作用下容易發(fā)生大幅振動,影響拉索的壽命以及橋梁的安全運營.因此,越來越多的學(xué)者開始關(guān)注拉索的風(fēng)致振動現(xiàn)象及其振動控制.對于傳統(tǒng)的被動阻尼器(如黏滯阻尼器、黏彈性阻尼器等),盡管目前應(yīng)用相對較多,但在阻尼可調(diào)節(jié)性或溫度穩(wěn)定性等方面仍存在一定的局限.

近年來,磁流變阻尼器作為一種智能阻尼器,在土木工程領(lǐng)域得到關(guān)注和應(yīng)用,其優(yōu)點是阻尼力可隨電流變化而改變,且阻尼可調(diào)節(jié)性及溫度穩(wěn)定性較好.Johnson等[1]針對拉索-磁流變阻尼器系統(tǒng)的半主動控制研究結(jié)果表明,與最優(yōu)被動控制相比,半主動控制下拉索的振動響應(yīng)可減小50%以上;Ni等[2]采用CLQG方法和神經(jīng)網(wǎng)絡(luò)技術(shù)設(shè)計出神經(jīng)控制器,用于研究半主動控制下拉索-磁流變阻尼器系統(tǒng)的狀態(tài)評估和反饋控制.上述研究中的主要工作集中在拉索和阻尼器的參數(shù)控制(如拉索的振動模態(tài)、長度、張力、阻尼器的安裝高度、安裝位置等)以及拉索面內(nèi)的振動控制.對于拉索的面外振動控制,Li 等[3]建立了拉索平面外風(fēng)雨激振的理論模型;Kang 等[4]對索拱結(jié)構(gòu)的面外動力特性進(jìn)行了研究;Wang 等[5]研究了拉索的面外對稱振動.然而,關(guān)于設(shè)置磁流變阻尼器的拉索面外振動研究還較少,且停留在擬半主動控制狀態(tài).鑒于此,本文根據(jù)線性自回歸模型方法得到風(fēng)速及風(fēng)荷載時程,將其導(dǎo)入ABAQUS軟件中分析因電流變化導(dǎo)致的阻尼變化及其對拉索風(fēng)振響應(yīng)的影響．編制了簡單Bang-Bang半主動控制程序,探討了拉索在半主動控制下的面外風(fēng)振響應(yīng),以考察磁流變阻尼器的減振能力.

1 風(fēng)速及風(fēng)荷載模擬

1.1 風(fēng)速模擬

拉索在風(fēng)力作用下,一般產(chǎn)生順風(fēng)向振動(面外振動)和橫風(fēng)向渦激振動,但橫風(fēng)向風(fēng)力提供的振動能量很小,因此本文僅考慮順風(fēng)向風(fēng)速的模擬.順風(fēng)向風(fēng)速由平均風(fēng)和脈動風(fēng)2部分組成.平均風(fēng)速隨高度變化的規(guī)律可用下式表示:

(1)

式中,U1和U2分別為高度Z1和Z2處的風(fēng)速;α為地面粗糙程度系數(shù).根據(jù)我國《公路橋梁抗風(fēng)設(shè)計規(guī)范》[6],取α=0.12,得到蘇通大橋南通地區(qū)的設(shè)計基本風(fēng)速為7.1 m高度處百年一遇風(fēng)速,即28.6 m/s． 南通地區(qū)平均風(fēng)速隨高度變化曲線見圖1．

本文采用線性濾波自回歸模型(AR)法模擬脈動風(fēng)速．根據(jù)文獻(xiàn)[6],順風(fēng)向脈動風(fēng)速一般可采用Kaimal譜模擬,其譜密度函數(shù)為

(2)

式中,Su為順風(fēng)向功率譜密度函數(shù);n為脈動頻率;f為莫寧坐標(biāo);u為氣流摩阻速度.

蘇通大橋距離主塔最遠(yuǎn)的斜拉索長 577 m,直徑為 20 cm.根據(jù)上述風(fēng)速模擬方法,對此斜拉索進(jìn)行風(fēng)速模擬.磁流變阻尼器的安裝高度離橋面3.5 m,約為橋塔頂部到橋面的1/66,故將該斜拉索分為66段,每段長度約為8.74 m,對應(yīng)的 67個節(jié)點從下到上順序編號為1～67.對各節(jié)點的水平脈動風(fēng)速和風(fēng)荷載時程進(jìn)行模擬,采樣時長為204.8 s,時間間隔為0.1 s,頻率分段數(shù)為2 048,截止頻率為10 rad/s.圖2給出了具有代表性的3個脈動風(fēng)速時程模擬結(jié)果,分別取自拉索底部、阻尼器安裝處以及拉索頂端,對應(yīng)于1#,2#,67#節(jié)點.

圖1 南通地區(qū)平均風(fēng)速隨高度變化曲線

圖2 脈動風(fēng)速時程模擬結(jié)果

1.2 風(fēng)荷載模擬

根據(jù)文獻(xiàn)[6],對于張拉結(jié)構(gòu)和索結(jié)構(gòu),瞬時風(fēng)載(即沿拉索長度方向的單位風(fēng)荷載)的計算公式為

(3)

式中,w(z,t)為t時刻瞬時風(fēng)載;ρ為空氣密度;V(z,t)為t時刻瞬時風(fēng)速;CH為斜拉索迎風(fēng)面阻力系數(shù);D為斜拉索等效直徑;β為斜拉索的安裝傾角．各參數(shù)的具體取值見表1.

表1 風(fēng)荷載時程模擬中的主要參數(shù)

編制程序計算斜拉索各節(jié)點的風(fēng)荷載時程．1#,2#,67#節(jié)點的脈動風(fēng)荷載時程模擬結(jié)果見圖3．

圖3 脈動風(fēng)荷載時程模擬結(jié)果

1.3 脈動風(fēng)模擬結(jié)果的驗證

為驗證模擬結(jié)果的正確性,選取1#,2#,67#節(jié)點模擬結(jié)果進(jìn)行頻域分析.圖4給出了3個節(jié)點的水平脈動風(fēng)速譜與Kaimal目標(biāo)譜的比較結(jié)果.由圖可知,模擬譜和目標(biāo)譜表現(xiàn)出良好的吻合性,從而驗證了模擬結(jié)果的正確性.

2 拉索有限元模型

2.1 有限元模型參數(shù)

根據(jù)文獻(xiàn)[7]中的Ernest公式,采用等效彈性模量法來修正拉索垂度的影響．在ABAQUS軟件中采用降溫法來實現(xiàn)拉索預(yù)應(yīng)力,其計算公式為

N=λEAΔT

(4)

式中,A為斜拉索截面積;E為彈性模量;ΔT為溫度變化量;λ為拉索溫度膨脹系數(shù).

對于蘇通大橋斜拉索上實際采用的磁流變阻尼器,根據(jù)文獻(xiàn)[8]中的試驗結(jié)果得到其隨電流變化的最大阻尼力,結(jié)果見表2.

一般而言,磁流變材料達(dá)到磁飽和狀態(tài)時,電流強(qiáng)度的增加反而會使磁場強(qiáng)度降低,磁流變阻尼器可提供的最大阻尼力減小.文獻(xiàn)[8]的試驗中僅給出了3.5 A時的數(shù)據(jù),繼續(xù)增大電流,阻尼力將會減小.

根據(jù)Weber等[8]針對蘇通大橋磁流變阻尼器的測試結(jié)果,可以得到黏滯阻尼系數(shù)為為試驗時阻尼器最大位移.文獻(xiàn)[8]采用振幅為20 mm、頻率為1 Hz的正弦波進(jìn)行阻尼力性能測試,根據(jù)得到的阻尼器滯回曲線,取fd=12.5 kN,xd=0.02 m,fr=1.29 Hz,代入式(5)進(jìn)行計算,得到黏滯阻尼系數(shù)為15.4 kN·s/m.

(5)

式中,fd為被動黏滯阻尼力;fr為拉索振動頻率;xd

(a) 1#節(jié)點

(b) 2#節(jié)點

(c) 67#節(jié)點

文獻(xiàn)[8]的試驗中未給出庫侖阻尼系數(shù),根據(jù)文獻(xiàn)[9]中磁流變阻尼器的設(shè)計方法反推此庫侖阻尼系數(shù)．設(shè)磁流變阻尼控制與最優(yōu)主動控制效果相同,其阻尼可調(diào)系數(shù)為m,采用最優(yōu)主動控制算法計算出黏滯阻尼系數(shù),根據(jù)阻尼可調(diào)倍數(shù)即可得到庫侖阻尼力.因此,在設(shè)計中認(rèn)為可調(diào)的庫侖阻尼系數(shù)一般為黏滯阻尼系數(shù)的整數(shù)倍．據(jù)此,基于文獻(xiàn)[8]中的試驗結(jié)果,其最大阻尼力為被動黏滯阻尼力的10倍,阻尼可調(diào)系數(shù)為10,最大庫侖阻尼系數(shù)為黏滯阻尼系數(shù)的9倍.在ABAQUS軟件中輸入的阻尼系數(shù)隨電流的變化情況見表3.

表3 阻尼系數(shù)隨電流變化的情況

2.2 拉索-磁流變阻尼器系統(tǒng)的有限元模型

在ABAQUS軟件中,拉索采用Euler梁單元進(jìn)行分析,劃分為544個單元．拉索兩端采用固定約束．預(yù)應(yīng)力采用文獻(xiàn)[10]中給出的數(shù)據(jù),取為12 764 kN．修正后的彈性模量為175 GPa,線膨脹系數(shù)為1.6×10-5/C°,單位質(zhì)量為102 kg/m,密度為7 800 kg/m3,截面面積為1.31×104mm2.采用降溫法模擬預(yù)應(yīng)力時,降低溫度取為348°.

建立拉索-磁流變阻尼器體系的有限元模型,其幾何關(guān)系示意圖見圖5.

圖5 拉索-磁流變阻尼器體系的幾何關(guān)系示意圖

3 拉索風(fēng)振響應(yīng)分析

本節(jié)計算不同電流情況下的結(jié)構(gòu)風(fēng)振響應(yīng),并取拉索中點作為控制點進(jìn)行分析.將1.2節(jié)模擬的風(fēng)荷載施加到結(jié)構(gòu)的各節(jié)點上,并考慮橋梁結(jié)構(gòu)設(shè)計基準(zhǔn)期為100 a,各節(jié)點間的組合系數(shù)為1.2.圖6給出了無阻尼器和電流為3.5 A狀態(tài)下的時程曲線對比,其他電流狀態(tài)下拉索各響應(yīng)的時程曲線與圖6類似.圖7給出了各響應(yīng)峰值隨電流變化的關(guān)系曲線.由圖可知,在無阻尼器以及施加電流分別為0,1.0,2.0,3.0,3.5 A六個狀態(tài)下,拉索中點的水平向位移峰值分別為15.3,14.8,13.3,12.6,12.5,12.4 mm;豎向位移峰值分別為35.0,33.9,31.4,30.2,30.0,29.9 mm;水平向加速度峰值分別為6.49,4.71,4.41,4.32,4.27,4.24 mm/s2;豎向加速度峰值分別為148.0,107.2,100.0,97.1,95.8,94.7 mm/s2.

在不出現(xiàn)磁飽和的情況下,隨著電流的增加,阻尼器的阻尼增大,結(jié)構(gòu)響應(yīng)減小.施加最大電流時的拉索風(fēng)振響應(yīng)最小,與不安裝阻尼器的情況相比,其水平向位移、豎向位移、水平向加速度、豎向加速度分別減小了18.9%,14.6%,34.7%,36.0%,表明蘇通大橋上所安裝的磁流變阻尼器可以通過電流調(diào)節(jié)來智能改變阻尼,以達(dá)到所需的減振效果.從圖7還可以發(fā)現(xiàn),電流為3.0和3.5 A時的減振效果差別不大,在實際應(yīng)用中可設(shè)置4個檔位(即電流為0,1.0,2.0,3.0 A).

(a) 水平向位移時程

(b) 水平向加速度時程

(c) 豎向位移時程

(d) 豎向加速度時程

4 拉索風(fēng)振響應(yīng)控制

半主動控制由于其良好的控制效果且耗能極少的優(yōu)點而成為近年來國內(nèi)外學(xué)者研究的熱點．本文研究的拉索-磁流變阻尼器系統(tǒng)是一種典型的半主動控制系統(tǒng).

(a) 水平向位移峰值

(b) 水平向加速度峰值

(c) 豎向位移峰值

(d) 豎向加速度峰值

在MATLAB軟件中利用經(jīng)典梁單元建立拉索的有限元模型,拉索運動方程中的質(zhì)量矩陣采用一致質(zhì)量矩陣,阻尼矩陣采用瑞利阻尼,其表達(dá)式為質(zhì)量矩陣和剛度矩陣的線性組合,即

C=α1M+β1K

(6)

式中,C,M,K分別為阻尼矩陣、質(zhì)量矩陣以及剛度矩陣;α1和β1為根據(jù)前2階振型阻尼比確定的比例系數(shù)．

4.1 拉索-磁流變阻尼器系統(tǒng)的半主動控制策略

為闡述半主動控制策略,首先介紹主動最優(yōu)控制力的計算.本文采用常見的LQR最優(yōu)主動控制算法來計算系統(tǒng)的最優(yōu)控制力,其中2個重要的控制參數(shù)權(quán)矩陣Q和R只能通過試算來確定.通過對模型進(jìn)行數(shù)次仿真,最終確定Q的控制系數(shù)為106,R的控制系數(shù)為0.1.在確定了Q和R后,采用MATLAB軟件中函數(shù)LQR得到連續(xù)狀態(tài)方程的控制力狀態(tài)反饋增益矩陣G,即

G=LQR(H,I,Q,R)

(7)

式中,H和I分別為狀態(tài)矩陣和控制矩陣.結(jié)合G和狀態(tài)向量Z,得到最優(yōu)控制力U=-GZ.

采用簡單Bang-Bang算法,編制了磁流變阻尼器-拉索系統(tǒng)的半主動控制算法,即

(8)

式中,ud(t)為控制力;fdymax為最大庫侖阻尼力;fdymin為最小庫侖阻尼力;x為結(jié)構(gòu)的位移. 該算法的控制策略為:當(dāng)結(jié)構(gòu)背離平衡位置運動時,阻尼控制器施加最大阻尼力;當(dāng)結(jié)構(gòu)靠近平衡位置運動時,施加最小阻尼力.

假定fdymin=0,令磁流變阻尼器提供的最大控制力等于LQR算法計算得出的最大主動最優(yōu)控制力Umax,并假定結(jié)構(gòu)的半主動控制與主動最優(yōu)控制的效果相同,阻尼可調(diào)系數(shù)m=10,則有

(9)

(10)

4.2 拉索-磁流變阻尼器系統(tǒng)的磁流變控制

以拉索中點作為振動控制的考察點,考察其水平向位移、豎向位移、水平向加速度和豎向加速度的減振情況,計算結(jié)果見圖8.

由圖8可知,最大主動最優(yōu)控制力為13.2 kN．在半主動控制狀態(tài)下拉索中點的水平向位移和豎向位移的最大響應(yīng)值分別為7.5和17.2 mm,水平向加速度和豎向加速度的最大響應(yīng)值分別為50.3和114.6 mm/s2.與無阻尼器狀態(tài)相比,這4類風(fēng)振響應(yīng)值分別減少了51%,51%,23%,23%.半主動控制下的拉索-磁流變阻尼器系統(tǒng)具有良好的面外減振效果,計算結(jié)果驗證了本文所編制的半主動控制程序的有效性.與面內(nèi)振動控制主要考慮控制系統(tǒng)的模態(tài)阻尼比不同,面外振動控制考慮的是風(fēng)荷載作用下的瞬態(tài)動力響應(yīng)控制,而關(guān)于拉索面外的瞬態(tài)風(fēng)振響應(yīng)控制研究目前鮮見文獻(xiàn)報道,因此本文可為今后此方面的研究提供一定的參考.

(a) 水平向位移時程

(b) 水平向加速度時程

(c) 豎向位移時程

(d) 豎向加速度時程

5 結(jié)論

1) 采用線性自回歸模型法,對蘇通大橋超長斜拉索的風(fēng)速時程以及風(fēng)荷載時程進(jìn)行了模擬．通過對比風(fēng)速功率譜的模擬值與目標(biāo)值,發(fā)現(xiàn)兩者具有良好的一致性.

2) 通過在ABAQUS軟件中建立拉索-磁流變阻尼器系統(tǒng)的有限元分析模型,計算了結(jié)構(gòu)不同電流情況下的風(fēng)振響應(yīng)．在不出現(xiàn)磁飽和的情況下,系統(tǒng)響應(yīng)隨阻尼的增加而不斷減小.

3) 在MATLAB軟件中依據(jù)簡單Bang-Bang算法編制了系統(tǒng)的半主動控制程序,對結(jié)構(gòu)進(jìn)行了振動控制計算. 仿真結(jié)果表明,半主動控制下的磁流變阻尼器可以使結(jié)構(gòu)的響應(yīng)大幅減小,產(chǎn)生良好的面外抑振效果.

)

[1] Johnson E A, Christenson R E, Spencer B F. Semiactive damping of cables with sag [J].Computer-AidedCivilandInfrastructureEngineering, 2003,18(2): 132-146. DOI:10.1111/1467-8667.00305.

[2] Ni Y Q, Chen Y, Ko J M, et al. Neuro-control of cable vibration using semi-active magneto-rheological dampers[J].EngineeringStructures, 2002,24(3): 295-307. DOI:10.1016/s0141-0296(01)00096-7.

[3] Li S, Chen Z, Wu T, et al. Rain-wind-induced in-plane and out-of-plane vibrations of stay cables[J].JournalofEngineeringMechanics, 2013,139(12): 1688-1698. DOI:10.1061/(asce)em.1943-7889.0000612.

[4] Kang H J, Zhao Y Y, Zhu H P. Out-of-plane free vibration analysis of a cable-arch structure[J].JournalofSoundandVibration, 2013,332(4): 907-921. DOI:10.1016/j.jsv.2012.10.007.

[5] Wang L, Zhao Y, Rega G. Multimode dynamics and out-of-plane drift in suspended cable using the kinematically condensed model[J].JournalofVibrationandAcoustics, 2009,131(6): 061008. DOI:10.1115/1.4000470.

[6] 中國工程建設(shè)標(biāo)準(zhǔn)化協(xié)會公路工程委員會. JTG/T D60-01—2004公路橋梁抗風(fēng)設(shè)計規(guī)范[S]. 北京:中交公路規(guī)劃設(shè)計院,2004.

[7] 楊佑發(fā), 白文軒, 魏建東. 采用Ernst公式計算斜拉橋拉索振動特性的精度分析[J]. 世界地震工程, 2008, 24(2): 50-53.

Yang Youfa, Bai Wenxuan, Wei Jiandong. Precision analysis of dynamic characteristics of cable using ernst’s modified formula in cable-stayed bridge[J].WorldEarthquakeEngineering, 2008,24(2): 50-53.(in Chinese)

[8] Weber F, Distl H. Amplitude and frequency independent cable damping of Sutong Bridge and Russky Bridge by magnetorheological dampers[J].StructuralControlandHealthMonitoring, 2014,22(2): 237-254. DOI:10.1002/stc.1671.

[9] 隋莉莉, 歐進(jìn)萍. MR減振驅(qū)動器用于結(jié)構(gòu)振動控制的算法研究[J]. 應(yīng)用力學(xué)學(xué)報, 2002, 19(3): 144-148. DOI:10.3969/j.issn.1000-4939.2002.03.032.

Sui Lili, Ou Jinping. The study on control algorithms for semi-active vibration control of structures with MR dampers[J].ChineseJournalofAppliedMechanics, 2002,19(3): 144-148. DOI:10.3969/j.issn.1000-4939.2002.03.032.(in Chinese)

[10] 張姝. 預(yù)應(yīng)力混凝土深受彎構(gòu)件力學(xué)性能研究[D]. 哈爾濱:東北林業(yè)大學(xué)土木工程學(xué)院,2015.