付章建， 岳祖潤， 梁文彥

(1.石家莊鐵道大學(xué) 土木工程學(xué)院，河北 石家莊 050043； 2.河北省建筑科學(xué)研究院有限公司，河北 石家莊 050021; 3.哈爾濱工程大學(xué) 航天與建筑工程學(xué)院，黑龍江 哈爾濱 150001)

隨著人類需求的多樣化，人類對于建筑結(jié)構(gòu)的要求已經(jīng)遠(yuǎn)遠(yuǎn)不限于滿足安全性、適用性、耐久性，而對于結(jié)構(gòu)舒適性提出了越來越多的新要求，尤其是近20年，樓蓋結(jié)構(gòu)的舒適度問題已經(jīng)引起了世界各國的廣泛關(guān)注，如果樓蓋結(jié)構(gòu)振動(dòng)問題控制不當(dāng)，不僅會(huì)影響建筑居住者的身體健康，也會(huì)影響精密設(shè)備的正常工作和壽命，甚至?xí)斐山Y(jié)構(gòu)的破壞，目前，英美等國已經(jīng)進(jìn)行了大量的實(shí)測研究，頒布了多項(xiàng)版本的規(guī)程、指南[1-9]。為了滿足結(jié)構(gòu)對于舒適度的要求，我國也頒布了一系列的相關(guān)規(guī)程、指南[10-16]，采用結(jié)構(gòu)豎向自振頻率、豎向加速度峰值等評價(jià)指標(biāo)進(jìn)行結(jié)構(gòu)舒適度評價(jià)。

本文通過實(shí)踐和理論相結(jié)合的方式，首先通過模態(tài)測試獲取所測試局部樓蓋的動(dòng)力特性，然后通過FEM進(jìn)行模擬分析。通過FEM分析與現(xiàn)有結(jié)果比較，驗(yàn)證了本文測試方法和模擬結(jié)果的正確性與可靠性，該測試方法不僅可以應(yīng)用于復(fù)雜結(jié)構(gòu)的結(jié)構(gòu)動(dòng)力特性預(yù)估及優(yōu)化設(shè)計(jì)，而且可以用于診斷和預(yù)報(bào)結(jié)構(gòu)系統(tǒng)的故障，識(shí)別結(jié)構(gòu)系統(tǒng)的載荷。

1 樓蓋模態(tài)測試

1.1 工程概況

本次測試工程為河北省衡水市桃城區(qū)某小學(xué)體育館，該結(jié)構(gòu)為地上2層的混凝土框架結(jié)構(gòu)，一層為梁板式樓蓋，平面布置圖如圖1所示。

圖1 一層樓蓋平面布置Fig.1 The plane layout of the first floor

本次測試區(qū)域?yàn)殛幱安糠?，測試區(qū)域橫向跨度7 500 mm，縱向跨度6 900 mm，橫向主梁截面尺寸為300 mm×650 mm，縱向主梁、次梁截面尺寸均為250 mm×600 mm，柱子截面尺寸均為400 mm×400 mm，樓板為現(xiàn)澆樓板，厚度為110 mm，結(jié)構(gòu)梁、板、柱混凝土強(qiáng)度等級(jí)均為C30。

結(jié)構(gòu)模態(tài)作為結(jié)構(gòu)的一個(gè)重要特性，在處理結(jié)構(gòu)振動(dòng)問題時(shí)，必須對其動(dòng)力特性有全面的了解，結(jié)構(gòu)的動(dòng)力特性通常用各階模態(tài)參數(shù)(模態(tài)頻率、模態(tài)振型和模態(tài)阻尼)描述，通過對結(jié)構(gòu)的模態(tài)測試可以得到準(zhǔn)確的動(dòng)力特性參數(shù)，進(jìn)而與理論計(jì)算結(jié)果進(jìn)行校驗(yàn)。

本模態(tài)測試用到的儀器設(shè)備主要有：16通道東華(DHDAS-5921)動(dòng)態(tài)信號(hào)采集分析系統(tǒng)一臺(tái)、100 m的2×0.3RVVP 信號(hào)屏蔽線、東華2D001VDH610V加速度傳感器及信號(hào)接收軟件等。其中，加速度傳感器靈敏度指標(biāo)在(300～370 mV/(m·s-2))，可以滿足對樓蓋結(jié)構(gòu)測試信號(hào)要求。現(xiàn)場測試圖如圖2所示。

圖2 動(dòng)力特性現(xiàn)場測試Fig.2 Field test of dynamic characteristics

1.2 測試原理

本測試采用自然激勵(lì)技術(shù)法(natural excitation technique，NEXT)，通過計(jì)算出自然激勵(lì)下實(shí)測信號(hào)的互功率譜，經(jīng)過傅里葉變換得到結(jié)構(gòu)測試各點(diǎn)與參考點(diǎn)之間的互相關(guān)函數(shù)，然后與時(shí)域識(shí)別方法特征系統(tǒng)實(shí)現(xiàn)算法(eigensystem realization algorithm，ERA)相結(jié)合進(jìn)行模態(tài)參數(shù)識(shí)別。

設(shè)一個(gè)n個(gè)自由度的線性體系，其一般運(yùn)動(dòng)方程為：

(1)

式中：M、C、K分別為質(zhì)量矩陣、阻尼矩陣、剛度矩陣；x為位移；f(t)為激勵(lì)力的列向量。

(2)

(3)

對于n自由度系統(tǒng)，A為2n×2n階矩陣，由于A矩陣的上半部分由零和單位矩陣組成，故A矩陣的未知元素為n×2n=2n2個(gè)，方程適用于任何時(shí)刻ti。

因此，如果在n個(gè)坐標(biāo)上對位移、速度、加速度進(jìn)行測量，每一組采樣y(i)可得到n個(gè)線性非奇次方程，如果在2n個(gè)不同瞬時(shí)進(jìn)行采樣，便可得到2n2個(gè)方程，即可確定n個(gè)未知量，則式(3)可以寫成：

(4)

(5)

令方程的解為：

y=ψexp(λit)

(6)

式中：λ為A的特征根，代入式(5)可得：

λψ=Aψ

(7)

從式(7)可以看出，對模態(tài)參數(shù)的識(shí)別問題可轉(zhuǎn)化為廣義特征值問題。A矩陣的特征值可給出系統(tǒng)的自然頻率、阻尼信息及其特征矢量，即模態(tài)矢量。

然而，欲進(jìn)行時(shí)域識(shí)別，則需要對加速度、速度、位移信息同時(shí)進(jìn)行測量，但在工程實(shí)踐中這樣做是不現(xiàn)實(shí)的，也是不需要的，對于n自由度的線性系統(tǒng)位移響應(yīng)可表示為：

(8)

對式(8)進(jìn)行2次微分，得到：

(9)

(10)

(11)

(12)

因此，只取得其中一種響應(yīng)數(shù)據(jù)，便可得到系統(tǒng)特征矩陣A，本次測試選用加速度作為響應(yīng)數(shù)據(jù)。

1.3 測試系統(tǒng)

本文采用的測試系統(tǒng)主要由4個(gè)部分構(gòu)成：

1)激振部分，由于該結(jié)構(gòu)不易獲得激勵(lì)力，且環(huán)境激勵(lì)下不會(huì)影響結(jié)構(gòu)的正常工作，更加符合實(shí)際工作狀態(tài)，因此，本測試采用環(huán)境激勵(lì)作為激振振源。

2)拾振部分，本測試采用拾振部分為接觸式的測量設(shè)備——壓電晶體式加速度傳感器，二次儀表采用電壓放大器。

3)數(shù)據(jù)采集及譜分析部分，將測試所得到的響應(yīng)信號(hào)傳遞到計(jì)算機(jī)上，對其進(jìn)行FFT變換的譜分析及傳遞函數(shù)分析。

4)曲線擬合及模態(tài)參數(shù)識(shí)別部分，采用計(jì)算機(jī)通過自互功率譜的方式進(jìn)行模態(tài)識(shí)別。

這4個(gè)部分協(xié)調(diào)工作，從而完成測試。

1.4 測點(diǎn)布置

根據(jù)《建筑工程容許振動(dòng)標(biāo)準(zhǔn)(GB 50868-2013)[16]的要求，結(jié)構(gòu)板構(gòu)件可視為連續(xù)彈性體，理論上具有無窮振型。在振動(dòng)過程中各點(diǎn)的振幅不同，具有隨機(jī)性，因此測試中采用多點(diǎn)測試統(tǒng)計(jì)平均方法。

此外，振動(dòng)測試方向應(yīng)與結(jié)構(gòu)樓板的垂直方向(法向)一致，同一構(gòu)件上的測試點(diǎn)應(yīng)等間隔均勻布置。對于板構(gòu)件的振動(dòng)測試，測點(diǎn)數(shù)量不應(yīng)少于5個(gè)。本測試對象為一層局部樓蓋3～4×B～C，采用35個(gè)測點(diǎn)進(jìn)行測試，以期得到良好的測試效果，拾振器參數(shù)設(shè)置如表1所示，測點(diǎn)布置如圖3所示。

圖3 測點(diǎn)布置Fig.3 The plane layout of the test point

表1 拾振器參數(shù)設(shè)置(采樣頻率100 Hz,電壓測量)Table 1 Parameter setting of vibration picker (sampling frequency 100 Hz, voltage measurement)

1.5 實(shí)測過程

由于本測試采用環(huán)境激勵(lì)的方式進(jìn)行，為減少信號(hào)干擾，所以測試時(shí)間選擇在晚上進(jìn)行。通過單點(diǎn)拾振獲取各測點(diǎn)的振動(dòng)加速度，為了提高測試效率與準(zhǔn)確性，本次測試共采用13個(gè)拾振器，整個(gè)測試過程中，A18傳感器不動(dòng)，其余傳感器測試完一組數(shù)據(jù)后進(jìn)行移動(dòng)，傳感器布置圖如圖4所示。通過3組測試完成了所測樓蓋的模態(tài)測試，所測得的加速度時(shí)程曲線如圖5所示。每組測試時(shí)間為25 min，數(shù)據(jù)采集頻率為100 Hz，采集完成后通過自互功率譜的方式進(jìn)行模態(tài)識(shí)別，確定模態(tài)參數(shù)。

圖4 傳感器布置Fig.4 Layout of sensors

圖5 7號(hào)測點(diǎn)加速度時(shí)程曲線Fig.5 Acceleration time history curve for measuring point 7

1.6 測試結(jié)果

采用自互功率譜對現(xiàn)場實(shí)測數(shù)據(jù)進(jìn)行模態(tài)參數(shù)的識(shí)別以及振型的識(shí)別，計(jì)算分析得到該結(jié)構(gòu)一層局部樓蓋3～4×B～C的一階、二階模態(tài)振型如圖6所示，振動(dòng)頻率和阻尼比如表2所示。

2.3.2 疾病經(jīng)濟(jì)風(fēng)險(xiǎn)度:由表5可知，門診自付費(fèi)用給低收入組和中低收入組的患者帶來了很大的風(fēng)險(xiǎn)(矯正RR值>1)，比如，低收入組門診就診經(jīng)濟(jì)風(fēng)險(xiǎn)度約是高收入組的18倍。隨著收入水平的增加，患者門診就診的經(jīng)濟(jì)風(fēng)險(xiǎn)也逐漸緩解。

表2 樓蓋頻率和阻尼比測試結(jié)果Table 2 Test results of floor frequency and damping ratio

圖6 實(shí)測樓蓋振型Fig.6 Vibration mode of measured floor

2 樓蓋模態(tài)FEM模擬

通過采用有限元方法對樓蓋模態(tài)進(jìn)行模擬計(jì)算，與試驗(yàn)結(jié)果進(jìn)行校核，從而判定建立有限元模型時(shí)所引入的一系列人為假設(shè)與實(shí)際情況是否相符，結(jié)構(gòu)的動(dòng)力特性是否符合要求，進(jìn)而應(yīng)用于進(jìn)一步的動(dòng)力特性研究。

2.1 有限元模型的建立

為了準(zhǔn)確模擬樓蓋結(jié)構(gòu)的振動(dòng)特性，本次分析采用整體建模的方式，盡可能的模擬實(shí)際工作狀態(tài)，以充分考慮局部樓蓋周邊構(gòu)件的約束作用。采用ANSYS軟件進(jìn)行三維建模，BEAM4 單元模擬柱單元，BEAM189單元模擬梁單元，SHELL63單元模擬板單元，采用映射網(wǎng)格劃分六面體單元，得到有限元模型如圖7所示。

圖7 ANSYS模型Fig.7 The model of ANSYS

2.2 參數(shù)選取

為了準(zhǔn)確模擬結(jié)構(gòu)的動(dòng)力特性，取混凝土密度為2 500 kg/m3，彈性模量為3.0×1010Pa，泊松比為0.2，混凝土板的阻尼取為0.05，本次有限元模態(tài)分析方法采用 Block Lanczos 模態(tài)分析方法，該法適用于大型有限元模型模態(tài)分析，并采用PCG求解器，得到了如圖8所示的一階、二階模態(tài)振型。

圖8 樓蓋振型模擬結(jié)果Fig.8 Simulation results of floor vibration mode

3 結(jié)果對比分析

本文對數(shù)值分析結(jié)果與實(shí)測結(jié)果進(jìn)行對比，如表3所示，可以看出有限元數(shù)值分析得到的模態(tài)參數(shù)與實(shí)測數(shù)據(jù)識(shí)別的模態(tài)參數(shù)結(jié)果擬合良好，最大偏差僅為5.9%，足以滿足工程分析的精度要求，從而驗(yàn)證了本文模態(tài)測試方法的正確性。

表3 2種方法結(jié)果對比Table 3 Comparison of the results of two methods

通過對有限元結(jié)果與實(shí)測結(jié)果的進(jìn)一步分析，可以得出以下結(jié)論：

1)由于所測樓蓋為局部樓蓋，只有當(dāng)結(jié)構(gòu)的模態(tài)較高時(shí)，才會(huì)表現(xiàn)出局部樓板的振動(dòng)，因此，需要在高階模態(tài)中準(zhǔn)確尋找所測樓板的振型及頻率。

2)2種方法得到的一階、二階振型圖基本吻合，且振動(dòng)頻率均在14.5 Hz左右，一層樓蓋3～4×B～C會(huì)出現(xiàn)顯著的振動(dòng)。最大振幅達(dá)到8 mm，說明 14.5 Hz 為該樓板的敏感振動(dòng)頻率，當(dāng)外界振動(dòng)接近該頻率附近時(shí)極易引發(fā)結(jié)構(gòu)的共振。

3) 本次動(dòng)力測試的實(shí)測結(jié)果與所建立有限元模型動(dòng)力特性較為相近，足以說明該模型可以用于進(jìn)一步的結(jié)構(gòu)動(dòng)力分析與研究。

4 結(jié)論

1)通過現(xiàn)場實(shí)測和有限元模擬結(jié)果得到的局部樓蓋振型圖及振動(dòng)頻率對比，驗(yàn)證了本文測試方法正確性及模擬結(jié)果可靠性。

2)采用環(huán)境激勵(lì)下的振動(dòng)響應(yīng)數(shù)據(jù)進(jìn)行模態(tài)參數(shù)識(shí)別對于復(fù)雜環(huán)境下大型結(jié)構(gòu)模態(tài)測試是一種行之有效的方法，而且更加符合實(shí)際情況和邊界條件。

為準(zhǔn)確得到局部樓蓋振動(dòng)特性，建議采用整體建模方式，以充分考慮樓蓋結(jié)構(gòu)周邊構(gòu)件約束作用，獲得可靠的模擬結(jié)果，利于開展進(jìn)一步的動(dòng)力分析和研究。

在環(huán)境激勵(lì)下，由于高階模態(tài)無法被激發(fā)出，因此，對于高階模態(tài)的識(shí)別需要進(jìn)一步研究。同時(shí)，如何將子結(jié)構(gòu)綜合成整體結(jié)構(gòu)進(jìn)行分析也是下一步研究的重點(diǎn)。