張大英 王錄民 王樹明

1 引 言

為了獲得工作狀態(tài)下結(jié)構(gòu)的動力特性，基于環(huán)境激勵的時域模態(tài)參數(shù)識別方法[1-5]得到了廣泛的應(yīng)用，尤其在識別大型橋梁結(jié)構(gòu)的動力特性、測定拉索索力、進(jìn)行健康檢測和損傷評估等方面取得了重要的研究成果[6-10]。鑒于該方法具有無須激勵設(shè)備、測試方便快捷、無須打斷結(jié)構(gòu)工作狀態(tài)、識別結(jié)果可靠等優(yōu)點，筆者將該方法應(yīng)用于鋼筋混凝土立筒倉結(jié)構(gòu)中展開了動力特性研究。

鋼筋混凝土立筒倉結(jié)構(gòu)剛度較大，不同于普通的建筑結(jié)構(gòu)，更不同于橋梁結(jié)構(gòu)。從外形上看，單倉多為高徑比較大的獨立薄壁筒體結(jié)構(gòu)[11]，很難將其簡化為平面結(jié)構(gòu)而獲取所需的真實模態(tài)；從支承形式上看，立筒倉有柱支承、筒壁支承和外筒內(nèi)柱支承等方式，其模態(tài)各不相同；從組成方式看，立筒倉分單倉和群倉，單倉為中心對稱結(jié)構(gòu)，群倉由單倉按照一定排列方式澆筑成一個整體得到，為軸對稱結(jié)構(gòu)，單倉和群倉的模態(tài)必然不同。因此，將基于環(huán)境激勵的時域模態(tài)參數(shù)識別方法應(yīng)用于立筒倉結(jié)構(gòu)中可以綜合考慮以上各個因素對其模態(tài)的影響。

然而，立筒倉結(jié)構(gòu)高度較大，傳感器需要固定在倉壁表面，試驗難度較大。為此，先在立筒倉結(jié)構(gòu)模型上開展環(huán)境激勵試驗，探索切實可行的測試方案，為測試在役立筒倉奠定基礎(chǔ)。文中以某單倉模型為研究對象，研究了其環(huán)境激勵測試方案，識別了其模態(tài)參數(shù)，研究結(jié)果可以用于指導(dǎo)在役立筒倉結(jié)構(gòu)的環(huán)境激勵測試和模態(tài)參數(shù)識別。

2 單倉的環(huán)境激勵測試

2.1 模型簡介

根據(jù)相似理論設(shè)計了某柱支承糧倉的縮尺模型，如圖1所示，分別為2×3的群倉模型和一個獨立的單倉模型。群倉模型利用鋼化玻璃做頂蓋，單倉模型沒有設(shè)計頂蓋。制作模型的材料為微?；炷?，彈性模量E=8.983×103MPa，泊松比μ=0.2，密度ρ=2 000 kg/m3。模型制作好后，先進(jìn)行了地震模擬振動臺試驗，研究了其在各級地震作用下的動力反應(yīng)及倉內(nèi)貯料的動力反應(yīng)特征[12，13]。因此，模型局部出現(xiàn)了損傷，但結(jié)構(gòu)整體依然牢固，可以對其進(jìn)行環(huán)境激勵測試獲得有效的動力響應(yīng)數(shù)據(jù)，從而識別結(jié)構(gòu)的模態(tài)參數(shù)。

圖1 立筒倉結(jié)構(gòu)模型Fig.1 Silos test model

文中對圖1中的單倉模型進(jìn)行環(huán)境激勵測試和模態(tài)參數(shù)識別研究。模型尺寸為：總高度2.5 m，柱子高度0.5 m，柱截面50 mm×50 mm，柱子數(shù)量12根，環(huán)梁中心到單倉頂部高度2 m，環(huán)梁高度50 mm，單倉外半徑389 mm，倉壁和漏斗壁厚度均為14 mm。測試時單倉內(nèi)沒有裝貯料。

2.2 測點布置方案

單倉為一筒體結(jié)構(gòu)，從幾何形狀分析，為一曲面體，要獲得它的模態(tài)參數(shù)尤其是振型，不能像測試橋梁結(jié)構(gòu)、高層房屋建筑結(jié)構(gòu)等一樣沿著x,y和z三個互相垂直的方向布置傳感器，應(yīng)該結(jié)合其自身的結(jié)構(gòu)特點沿著倉壁圓周的徑向(定義為R向)和切向(定義為θ向)布置傳感器。從結(jié)構(gòu)剛度分析，單倉結(jié)構(gòu)剛度大，環(huán)境激勵能量較弱，因此測試時需要選擇動力反應(yīng)較明顯的部位布置測點，如倉頂部和柱子部位；單倉各部分剛度分布不均，柱子是整個結(jié)構(gòu)的剛度薄弱部分，而柱子以上的筒體結(jié)構(gòu)剛度大，因此，環(huán)梁為結(jié)構(gòu)剛度突變的分界線，在環(huán)梁上需要布置測點。此外，要觀察筒體部分的動力反應(yīng)，需要在筒體范圍內(nèi)再布置一定數(shù)量的測點。圖2給出了單倉測點布置方案的立面圖和平面圖。

圖2 單倉測點布置圖Fig.2 Arrangement diagram of measuring points of single silo

分析圖2(a)僅在部分柱子的某些部位布置了測點，主要原因是：①地震模擬振動臺試驗對柱子的局部有損傷，這些部位不能牢固安放傳感器；②與群倉模型相鄰的柱子，由于操作空間小無法布設(shè)傳感器。分析圖2(b)僅在圓周270°范圍內(nèi)布置了測點，其余90°范圍內(nèi)沒有布置測點，原因是單倉的這一部位與群倉模型相鄰，空隙較小，沒有足夠的工作空間。

2.3 測試方案

圖2共設(shè)置了36個測點位置，每個測點處傳感器的測試方向可能為R向、θ向或既有R向又有θ向，因此自制了薄壁角鋼，如圖3所示，可以在角鋼的兩個肢背上同時布置沿著R向和θ向的兩個傳感器。測試時先將角鋼用膨脹螺栓固定在倉壁表面，保持橫平豎直，然后再將傳感器固定于角鋼上，如圖4所示。

圖3 角鋼固定圖Fig.3 Steel angles fixed to silo wall

圖4 傳感器固定圖Fig.4 Sensors fixed to steel angles

不考慮R向和θ向的耦合作用對單倉模態(tài)的影響，首先根據(jù)傳感器的測試方向?qū)y點分為兩組：第一組R向和第二組θ向。第一組分四批進(jìn)行測試(表1)，第二組分兩批進(jìn)行測試(表2)，兩組測試的參考點為C37號，位置統(tǒng)一，測試方向不同。以單倉中心指向外為R+，指向內(nèi)為R-，定義測試方向繞單倉中心逆時針旋轉(zhuǎn)為θ+，順時針旋轉(zhuǎn)為θ-。

表1單倉R+向測點分批

Table1Measuringpointsinbatchesinpositiveradialdirection

批次測點一C2C3C4C5C6C7C8C9C10二C11C12C13C14C15C16C17C18C19C20三C21C22C23C24C25C26C27C28C29C30四C31C32C33C34C35C36C1

表2單倉θ+向測點分批

Table2Measuringpointsinbatchesinpositivehoopdirection

批次測點一C1C31C32C33C34C35C36 二C21C22C23C24C25C26C27C28C29C30

2.4 測試數(shù)據(jù)及預(yù)處理

設(shè)置采樣頻率200 Hz，上限頻率100 Hz，對各批測點進(jìn)行多次環(huán)境激勵測試。測試時實時觀察各測點的加速度響應(yīng)曲線、頻譜圖、概率密度函數(shù)、概率分布函數(shù)等，發(fā)現(xiàn)異常及時停止采集數(shù)據(jù)，找出問題并解決后再開始采集數(shù)據(jù)。

對采集得到的原始數(shù)據(jù)需要進(jìn)行預(yù)處理[14]后再用來分析系統(tǒng)的模態(tài)參數(shù)，這里采用最小二乘法和五點三次平滑法消除原始信號的趨勢項，使振動曲線更加光滑，但是幅值隨平滑次數(shù)m的增大而減小，因此需要根據(jù)實際情況，結(jié)合多次試算，確定合適的平滑次數(shù)。設(shè)置上下限頻率分別為45 Hz和5 Hz，對原始信號采用帶通窗函數(shù)法進(jìn)行數(shù)字濾波，僅保留所需范圍的頻率信號。圖5給出了R向第一批測點C3的加速度原始數(shù)據(jù)和預(yù)處理數(shù)據(jù)曲線。

圖5 R+向第一批測點C3的加速度信號Fig.5 Acceleration signal of C3 in first batch (R+)

3 單倉模態(tài)參數(shù)識別

3.1 頻率和阻尼比識別

對所有測試得到的測點加速度信號進(jìn)行預(yù)處理后，采用改進(jìn)的數(shù)據(jù)驅(qū)動隨機子空間方法(Improved-DD-SSI方法)[15]識別單倉的頻率f、振型ψ和阻尼比ξ。在計算過程中截取一段數(shù)據(jù)后，確定Improved-DD-SSI方法中建立H矩陣的行數(shù)2i，然后進(jìn)行數(shù)值計算。根據(jù)試驗測試時布置的參考點數(shù)目rE，設(shè)定計算時用到的參考點個數(shù)r≤rE，設(shè)定頻率容差ew<2%、阻尼比容差eξ<5%和模態(tài)置信因子MAC>90%，作出穩(wěn)定圖，用以確定系統(tǒng)的階數(shù)N。從而計算出對應(yīng)N值的系統(tǒng)頻率f、振型ψ和阻尼比ξ。

圖6給出了單倉R向第一批、第二批、第三批和第四批計算得到的穩(wěn)定圖，圖7給出了單倉θ向第一批和第二批計算得到的穩(wěn)定圖。穩(wěn)定圖中符號“·f” 代表頻率穩(wěn)定，“·fv”代表頻率和振型穩(wěn)定，“·fz”代表頻率和阻尼比穩(wěn)定，“*”代表穩(wěn)定極點。穩(wěn)定圖中同時作出了各批次測點的自功率譜疊加曲線。

分析穩(wěn)定圖圖6和圖7，每批測點的自功率譜疊加曲線的峰值個數(shù)不同，需要找到各批數(shù)據(jù)中相一致的峰值才能確定結(jié)構(gòu)模態(tài)，可以發(fā)現(xiàn)利用峰值拾取法很難獲得結(jié)構(gòu)的高階模態(tài)。再分析圖中的穩(wěn)定極點，對應(yīng)結(jié)構(gòu)的低階模態(tài)， 當(dāng)系統(tǒng)階數(shù)N<10時穩(wěn)定極點就出現(xiàn)了，而對應(yīng)結(jié)構(gòu)的高階模態(tài)，當(dāng)系統(tǒng)階數(shù)N較大時，穩(wěn)定極點才出現(xiàn)。另外，在計算中還發(fā)現(xiàn)，穩(wěn)定極點的出現(xiàn)亦與塊H矩陣的行數(shù)2i有關(guān)，若i較小，則很難找到穩(wěn)定極點，因此需要經(jīng)過多次試算找到合適的i值，然后利用Improved-DD-SSI方法識別結(jié)構(gòu)的模態(tài)參數(shù)。

圖6 單倉R向各批次數(shù)據(jù)穩(wěn)定圖Fig.6 Stabilization diagrams for different batches in R direction

圖7 單倉θ向各批次數(shù)據(jù)穩(wěn)定圖Fig.7 Stabilization diagrams for different batches in θ direction

表3列出了利用時域Improved-DD-SSI方法和頻域峰值拾取法識別得到的單倉頻率和阻尼比。分析兩種方法的識別結(jié)果，各階頻率值非常接近，相互驗證了識別結(jié)果的準(zhǔn)確性。但各階阻尼比存在一定的誤差，Improved-DD-SSI方法識別的一階第一個阻尼比和二階第二個阻尼比大于峰值拾取法識別結(jié)果，其余阻尼比前者均小于后者。由于阻尼比識別結(jié)果的準(zhǔn)確性和精度目前仍然是個難題，所以在此不以阻尼比作為判斷識別結(jié)果準(zhǔn)確性的依據(jù)。

3.2 模態(tài)振型識別

圖8給出了利用Improved-DD-SSI方法識別得到的單倉歸一化振型圖，圖中一階、二階和四階振型都是用兩視圖表示的。根據(jù)測點布置方案，識別得到了單倉前四階模態(tài)振型。第一階頻率f=11.98Hz對應(yīng)的振型為彎曲振型，如圖8(a)所示；第二階頻率f=21.70Hz對應(yīng)的振型為筒壁翹曲，振型平面形狀近似為橢圓，如圖8(b)所示；第三階頻率f=30.69Hz對應(yīng)的振型為立筒倉倉壁隨柱子的扭轉(zhuǎn)振型，如圖8(c)所示；第四階頻率f=38.41Hz對應(yīng)的振型為彎剪振型，振型平面上有三個外凸和三個內(nèi)凹的波形，如圖8(d)所示。表3中各階模態(tài)其他頻率對應(yīng)的振型是沿與圖8中同階振型方向垂直的相同振型，沒有在圖中列出。

4 結(jié) 論

通過對立筒單倉模型的環(huán)境激勵測試、模態(tài)參數(shù)識別和結(jié)果分析，得到以下結(jié)論：

(1) 立筒單倉為薄壁筒體結(jié)構(gòu)，環(huán)境激勵測試時采用自制角鋼固定傳感器的位置，為R向和θ向同時采集數(shù)據(jù)提供了最佳方案，試驗測點數(shù)目較多，分批進(jìn)行測試，并設(shè)置了共同的參考點實現(xiàn)振型數(shù)據(jù)歸一化。

(2) 環(huán)境激勵法測試結(jié)果和時域Improved-DD-SSI方法識別結(jié)果表明，采用自然環(huán)境振動的方法，足以識別出立筒單倉所感興趣的主要模態(tài)參數(shù)，是立筒倉結(jié)構(gòu)動力實驗的實用方法。

圖8 單倉模態(tài)振型Fig.8 Vibration mode shapes of single silo

表3單倉頻率和阻尼比

Table3Frequenciesanddampingratiosofsinglesilo

模態(tài)階數(shù)Improved-DD-SSI方法頻率/Hz阻尼比峰值拾取法頻率/Hz阻尼比111.982.67%11.931.57%13.850.54%13.691.72%218.820.70%18.771.02%21.700.89%21.70.6%330.690.03%29.720.44%438.410.51%38.711.01%40.750.44%40.861.31%

(3) 對比分析Improved-DD-SSI方法和傳統(tǒng)峰值拾取法識別結(jié)果，實現(xiàn)了互補和驗證，識別低階模態(tài)時兩者差異不大，但識別高階模態(tài)時前者優(yōu)于后者。

立筒倉屬于剛度較大的高聳構(gòu)筑物，它的阻尼機制非常復(fù)雜，尤其是工作狀態(tài)下倉內(nèi)裝有糧食、煤炭等貯料，存在使用過程中裝料、卸料的問題，而且為了保持倉內(nèi)適宜的溫度和濕度，必然需要通風(fēng)等機械設(shè)備的持續(xù)運轉(zhuǎn)，產(chǎn)生的噪聲對識別結(jié)果的影響及影響程度需要作進(jìn)一步的研究。此外，通過環(huán)境激勵測試識別的立筒倉阻尼的實用性亦是需要進(jìn)一步研究的課題。

[ 1 ] van Overschee P, de Moor BL. Subspace identification for linear systems: theory-implementation-application [M]. Dordrecht, Netherlands: Kluwer Academic Publishers, 1996.

[ 2 ] 常虹，竇立軍，殷琨.基于環(huán)境激勵的模態(tài)參數(shù)識別方法比較研究[J].長春工程學(xué)院學(xué)報(自然科學(xué)版)，2009，10(2)：10-12．

Chang Hong, Dou Lijun, Yin Kun. Comparative study of modal parameter identification methods on ambient excitation[J]. Journal of Changchun Institute of Technology(Natural Science Edition), 2009, 10(2): 10-12. (in Chinese)

[ 3 ] Rainieri C, Fabbrocino G. Automated output-only dynamic identification of civil engineering structures [J]. Mechanical Systems and Signal Processing, 2010, 24: 678-695.

[ 4 ] Broownjohn J M W, Magalhaes F, Caetano E, et al. Ambient vibration re-testing and operational modal analysis of the Humber Bridge [J]. Engineering Structures, 2010, 32: 2003-2018.

[ 5 ] 付春，姜紹飛，牟海東.基于改進(jìn)HHT的結(jié)構(gòu)模態(tài)參數(shù)識別方法[J]．應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報，2011，19(4)：583-590．

Fu Chun, Jiang Shaofei, Mu Haidong. Structural modal parameters identification method based on improved Hilbert-Huang transform[J]. Journal of Basic Science and Engineering, 2011, 19(4): 583-590. (in Chinese)

[ 6 ] 樊可清，倪一清，高贊明．改進(jìn)隨機子空間系統(tǒng)識別方法及其在橋梁狀態(tài)檢測中的應(yīng)用[J]．中國公路學(xué)報，2004，17(4)：70-73．

Fan Keqing, Ni Yiqing, Gao Zanming. Improved stochastic system identification approach with its application in bridge condition monitoring[J]. China Journal of Highway and Transport, 2004, 17(4):70-73.(in Chinese)

[ 7 ] 高 敏，劉劍鋒，張啟偉.超高索塔自立狀態(tài)下環(huán)境振動試驗與分析[J].結(jié)構(gòu)工程師，2008，24(3)：129-134.

Gao Min, Liu Jianfeng, Zhang Qiwei. Ambient vibration test of super-high cable-tower in its self-erect state[J]. Structural Engineers, 2008, 24(3):129-134. (in Chinese)

[ 8 ] 蒲黔輝，秦世強，施洲,等．環(huán)境激勵下鋼筋混凝土拱橋模態(tài)參數(shù)識別[J]．西南交通大學(xué)學(xué)報，2012，47(4)：539-545．

Pu Qianhui, Qin Shiqiang, Shi Zhou, et al. Modal parameter identification of reinforced concrete arch bridge under ambient excitation[J]．Journal of Southwest Jiaotong University, 2012, 47(4): 539-545. (in Chinese)

[ 9 ] 葉錫鈞，顏全勝，李 健 等．基于環(huán)境激勵的大跨度斜拉橋模態(tài)參數(shù)和索力識別[J]．振動與沖擊，2012，31(16)：157-163．

Ye Xijun, Yan Quansheng, Li Jian et al. Modal identifi cation and cable tension estimation of long span cable-stayed bridge based on ambient excitation[J]．Journal of Vibration and Shock, 2012, 31(16): 157-163. (in Chinese)

[10] 曹碩，張啟偉.基于實測的蘇通大橋疲勞壽命評估[J].結(jié)構(gòu)工程師，2011，27(2)：84-89.

Cao Shuo, Zhang Qiwei. Fatigue evaluation of Sutong bridge based on actual measurement[J]. Structural Engineers, 2011, 27(2):84-89. (in Chinese)

[11] 龔耀清，張偉杰．高層建筑結(jié)構(gòu)筒體的共振響應(yīng)[J]．廣西大學(xué)學(xué)報(自然科學(xué)版)，2010，35(4)：549-554．

Gong Yaoqing, Zhang Weijie. Resonant response of tubular structure of high-rise buildings[J]. Journal of Guangxi University(Natural Science Edition), 2010, 35(4): 549-554. (in Chinese)

[12] 王錄民，張華.鋼筋混凝土立筒群倉結(jié)構(gòu)抗震模型的設(shè)計與制作[J].結(jié)構(gòu)工程師，2008，24(3)：126-128.

Wang Lumin, Zhang Hua. Design and construction for a seismic model of reinforced concrete group silos[J]. Structural Engineers, 2008, 24(3):126-128. (in Chinese)

[13] 王錄民，張華，盧文勝,等.柱承式群倉結(jié)構(gòu)模型模擬地震振動臺試驗研究[J].建筑結(jié)構(gòu),2010,(10)：41-43．

Wang Lumin, Zhang Hua, Lu Wensheng, et al. Study of shaking table tests on the model of group silo structures [J]．Building Structure, 2010, (10): 41-43. (in Chinese)

[14] 王濟(jì)，胡曉．MATLAB在振動信號處理中的應(yīng)用[M]．北京：中國水利水電出版社，2006.

Wang Ji, Hu Xiao. MATLAB applications in vibration signal processing[M]. Beijing: China WaterPower Press, 2006. (in Chinese)

[15] 張大英．基于改進(jìn)的數(shù)據(jù)驅(qū)動隨機子空間方法的鋼筋混凝土立筒倉動力參數(shù)研究[D]．合肥：合肥工業(yè)大學(xué)，2011．

Zhang Daying. Dynamic parameters research for reinforced concrete silos based on improved data-driven stochastic subspace method [D]. Hefei: Hefei University of Technology, 2011. (in Chinese)