周志紅，吳邵慶，費(fèi)慶國(guó)

(1.東南大學(xué)土木工程學(xué)院，江蘇南京 210096;2.江蘇省工程力學(xué)分析重點(diǎn)實(shí)驗(yàn)室，江蘇南京 210096)

隨著我國(guó)公路交通事業(yè)的發(fā)展，超重超速車(chē)輛逐漸成為人們關(guān)注的主要問(wèn)題，直接影響橋梁的工作狀態(tài)和使用壽命［1］。同時(shí)，隨著許多大跨度、輕質(zhì)化、大柔度新型橋梁的大量涌現(xiàn)，車(chē)載特性直接影響橋梁動(dòng)力響應(yīng)特性，這使得移動(dòng)車(chē)輛荷載作用下橋梁的動(dòng)力響應(yīng)問(wèn)題成為有關(guān)研究者關(guān)注的熱點(diǎn)［2］。

在車(chē)橋耦合振動(dòng)問(wèn)題中，車(chē)輛模型可分為整車(chē)模型、半車(chē)模型和單輪車(chē)輛模型等［3-4］。求解車(chē)橋系統(tǒng)動(dòng)力學(xué)方程的方法又分為直接積分法和振型疊加法。趙永剛等［5］考慮橋梁路邊不平整度，建立了大跨度斜拉橋空間有限元模型，并利用隨機(jī)分析方法獲得橋梁動(dòng)響應(yīng)。肖新標(biāo)等［6］建立了勻速移動(dòng)質(zhì)量模型作用下的簡(jiǎn)支梁動(dòng)態(tài)響應(yīng)模型，實(shí)現(xiàn)了車(chē)橋耦合振動(dòng)系統(tǒng)的仿真，其仿真結(jié)果具有快速、簡(jiǎn)單和靈活的特點(diǎn)。Moghimi等［7］利用仿真分析了車(chē)速等因素對(duì)車(chē)橋耦合系統(tǒng)的影響。王解軍等［8］建立了適用于大跨橋梁車(chē)輛振動(dòng)計(jì)算的車(chē)橋耦合單元模型，分析比較了重型汽車(chē)與輕型汽車(chē)作用下的動(dòng)力沖擊系數(shù)。然而以上方法均側(cè)重于理論分析，未應(yīng)用于實(shí)際工程。劉文淑等［9］利用小波的多尺度分析，對(duì)實(shí)測(cè)鋼箱梁上的應(yīng)力數(shù)據(jù)去噪，使得處理后的數(shù)據(jù)更能反映真實(shí)應(yīng)力水平，但未應(yīng)用實(shí)測(cè)響應(yīng)數(shù)據(jù)進(jìn)行車(chē)橋系統(tǒng)參數(shù)的識(shí)別。

本文將小波變換技術(shù)引入車(chē)軸輪間距和車(chē)速監(jiān)測(cè)中，提出一種針對(duì)大跨度橋梁上車(chē)輪軸間距和車(chē)速監(jiān)測(cè)的新型數(shù)據(jù)處理方法。根據(jù)車(chē)輪過(guò)橋時(shí)橋面的動(dòng)響應(yīng)，實(shí)現(xiàn)對(duì)大跨度橋梁結(jié)構(gòu)上車(chē)輛行駛狀況的監(jiān)測(cè)，所提方法能夠較準(zhǔn)確地得到車(chē)輪軸間距、車(chē)速和車(chē)輛進(jìn)入/離開(kāi)橋梁的時(shí)間間隔。

1 基于小波變換的車(chē)輪軸間距和車(chē)速狀態(tài)監(jiān)測(cè)原理

小波變換作為一種時(shí)頻分析工具，與傳統(tǒng)傅里葉變換相比，能夠很好地捕捉不同頻率下信號(hào)分量的局部信息，并準(zhǔn)確解構(gòu)和重建有限的非周期和非平穩(wěn)信號(hào)［10-12］，在信號(hào)分析和參數(shù)識(shí)別中具有明顯的優(yōu)勢(shì)，并有較多的工程應(yīng)用。小波變換包括連續(xù)小波變換(CWT)和離散小波變換(DWT)等，本文利用基于DWT的方法對(duì)橋上車(chē)輛狀態(tài)進(jìn)行檢測(cè)識(shí)別分析。具體思路為:利用小波的多尺度分析將車(chē)輛經(jīng)過(guò)時(shí)實(shí)測(cè)橋面動(dòng)應(yīng)變數(shù)據(jù)進(jìn)行多尺度分解，利用每個(gè)尺度上的小波系數(shù)的峰值來(lái)檢測(cè)車(chē)輪入橋和出橋的時(shí)刻，通過(guò)車(chē)輛過(guò)橋的時(shí)間以及前后輪在橋上停留時(shí)間來(lái)計(jì)算車(chē)的實(shí)際速度，然后根據(jù)前后輪入橋和出橋引起的動(dòng)應(yīng)力峰值來(lái)識(shí)別車(chē)輪軸間距。

2 青衣南大橋?qū)嵗?yàn)證

為驗(yàn)證本文方法在大跨度橋梁結(jié)構(gòu)車(chē)輛車(chē)輪軸間距和車(chē)速監(jiān)測(cè)中的適用性，將該方法應(yīng)用于香港青衣南大橋監(jiān)測(cè)中。香港青衣南大橋?yàn)槿邕B續(xù)混凝土橋梁，總長(zhǎng)73 m，寬10.58 m，水平方向傾斜角為27°，即該橋的縱向斜率為6.75%。采用雙箱梁型截面，鋼筋混凝土結(jié)構(gòu)的等效彈性模量為26 GPa，密度為24.5 kN/m2。具體步驟為:(a)對(duì)香港青衣南大橋建立有限元模型，利用動(dòng)力學(xué)分析得到行車(chē)狀態(tài)下橋面各處的動(dòng)態(tài)應(yīng)變信號(hào)，利用基于DWT的方法進(jìn)行車(chē)速、車(chē)輪軸間距的識(shí)別;(b)根據(jù)橋上實(shí)測(cè)數(shù)據(jù)，同樣利用基于DWT的方法進(jìn)行車(chē)速、車(chē)輪軸間距的識(shí)別，驗(yàn)證本文提出的方法的適用性。

2.1 有限元數(shù)值仿真

圖1為建立的有限元模型以及測(cè)點(diǎn)排布示意圖。假設(shè)車(chē)輪軸間距為4.3 m，車(chē)速為20 m/s情況下，利用基于有限元模型動(dòng)力學(xué)分析結(jié)果得到各處的動(dòng)應(yīng)變結(jié)果，動(dòng)應(yīng)變信號(hào)的采樣頻率為100 Hz。

圖1 有限元模型中應(yīng)變片排布位置Fig.1 Strain gauge installation on finite element model

基于DWT的方法，采用‘db4’小波對(duì)單一車(chē)輛入橋時(shí)橋上動(dòng)應(yīng)變信號(hào)進(jìn)行10個(gè)尺度的分解，圖2為某動(dòng)應(yīng)變信號(hào)在尺度1的細(xì)節(jié)信號(hào)。從圖2可以發(fā)現(xiàn)，信號(hào)中出現(xiàn)了與兩車(chē)輪軸入橋相對(duì)應(yīng)的2個(gè)峰值，峰值之間包含22個(gè)時(shí)間步長(zhǎng)，即后軸入橋比前軸早0.22 s。如假設(shè)車(chē)速為20 m/s且事先已知，則可以計(jì)算得到車(chē)輪軸間距為4.4 m，與實(shí)際車(chē)輪軸間距4.3 m非常接近。識(shí)別橋上車(chē)輛的車(chē)速可以根據(jù)已知距離的2個(gè)應(yīng)變片的信號(hào)或其DWT分解結(jié)果計(jì)算得到。車(chē)載作用下的橋梁，由于車(chē)輪作用點(diǎn)越靠近應(yīng)變片，對(duì)應(yīng)的動(dòng)應(yīng)變信號(hào)幅值就越大，因此，可以將兩應(yīng)變片的距離除以?xún)蓱?yīng)變信號(hào)達(dá)到峰值的時(shí)間差來(lái)獲得車(chē)輪的行駛速度。

2.1.1 噪聲水平對(duì)識(shí)別結(jié)果的影響

由于試驗(yàn)數(shù)據(jù)往往包含噪聲，在此考慮信噪比(SNR)分別為10dB和20dB的車(chē)輪軸間距識(shí)別效果。依然采用基于DWT的方法將帶噪聲的信號(hào)進(jìn)行10個(gè)尺度的分解，取出第4個(gè)尺度的細(xì)節(jié)來(lái)識(shí)別車(chē)輪軸間距，識(shí)別結(jié)果見(jiàn)表1。由表1可知，包含噪聲的原始信號(hào)的識(shí)別精度低于無(wú)噪聲情況，隨著動(dòng)應(yīng)變信號(hào)信噪比的降低，車(chē)輪軸間距的識(shí)別誤差會(huì)增加。

圖2 某動(dòng)應(yīng)變信號(hào)在尺度1的細(xì)節(jié)信號(hào)Fig.2 Detail strain signal at level 1

2.1.2 小波類(lèi)型和車(chē)速對(duì)識(shí)別結(jié)果的影響

采用不同小波類(lèi)型以及車(chē)速對(duì)識(shí)別結(jié)果的影響見(jiàn)表2。由表2可以看出，不同的小波類(lèi)型對(duì)識(shí)別結(jié)果有一定的影響，除基于‘db4’小波的識(shí)別精度隨著車(chē)速的增大而降低以外，其他類(lèi)型的小波識(shí)別精度都隨著車(chē)速的增大而減小。其原因可以歸結(jié)于各小波的特點(diǎn):‘db4’小波從本質(zhì)上屬于離散小波，能夠用小波系數(shù)精確地描述信號(hào)細(xì)節(jié);而其他小波如‘haar’或‘sym4’，其本身就有明確的表達(dá)式。當(dāng)車(chē)速較快，動(dòng)應(yīng)變信號(hào)包含較多高頻分量時(shí)，細(xì)節(jié)分量的峰值將會(huì)出現(xiàn)偏移，導(dǎo)致識(shí)別精度下降?？傮w來(lái)說(shuō)，隨著車(chē)速的降低，車(chē)輪軸間距的識(shí)別精度有所提高，因此在車(chē)速較高時(shí)，建議采用‘db4’小波。

表1 考慮噪聲時(shí)的車(chē)輪軸間距識(shí)別Table 1 Vehicle axle spacing identification with consideration of noise

表2 不同小波類(lèi)型以及車(chē)速下車(chē)輪軸間距的識(shí)別結(jié)果Table 2 Vehicle axle spacing identification with different wavelets and vehicle speeds

2.2 基于實(shí)測(cè)數(shù)據(jù)監(jiān)測(cè)和識(shí)別驗(yàn)證

為了應(yīng)用本文提出的方法進(jìn)行車(chē)輪軸間距的識(shí)別和車(chē)速的檢測(cè)，在香港青衣南大橋上進(jìn)行振動(dòng)測(cè)試。在三跨橋的邊跨布置傳感器，標(biāo)準(zhǔn)卡車(chē)經(jīng)過(guò)橋面時(shí)，獲得如圖3所示傳感器布置方案下的2，5，8，11和14位置的動(dòng)應(yīng)變數(shù)據(jù)?，F(xiàn)場(chǎng)測(cè)試采樣頻率為500 Hz，所用的標(biāo)準(zhǔn)車(chē)輛為三軸車(chē)輛，其中第1和第2輪軸間距為4.43 m，第2和第3輪軸間距為1.32 m。利用布置的兩已知距離的光柵傳感器測(cè)得經(jīng)過(guò)時(shí)刻，將傳感器距離除以經(jīng)過(guò)的總時(shí)間計(jì)算得到(在試驗(yàn)中盡量保持車(chē)輛勻速行駛)車(chē)速為15.28 m/s。

2.2.1 車(chē)輪軸間距和車(chē)速監(jiān)測(cè)

由于實(shí)測(cè)數(shù)據(jù)中橋面動(dòng)應(yīng)變數(shù)據(jù)的信噪比較低，利用基于DWT的方法得到各細(xì)節(jié)部分信號(hào)由于車(chē)輪經(jīng)過(guò)應(yīng)變片產(chǎn)生的峰值較難判斷，因此還需引入信號(hào)的近似部分來(lái)識(shí)別車(chē)輪軸間距和車(chē)速。信號(hào)的近似部分受噪聲影響較小。圖4為布置在第1跨跨中兩邊的典型應(yīng)變片采集到的信號(hào)的近似部分，由圖4可知，近似信號(hào)中的突變由車(chē)輪軸入橋引起，可以根據(jù)3個(gè)突變位置的時(shí)刻來(lái)計(jì)算軸間距。由此識(shí)別得到的車(chē)輪軸間距和車(chē)速見(jiàn)表3(采用第4尺度的近似信號(hào))。

圖3 現(xiàn)場(chǎng)測(cè)試中應(yīng)變片排布位置Fig.3 Strain gauge installation in field test

由于試驗(yàn)測(cè)得的信號(hào)信噪比較低，利用信號(hào)細(xì)節(jié)部分識(shí)別會(huì)難以判斷峰值，通過(guò)近似信號(hào)可以較好地識(shí)別間距相對(duì)較大的第1和第2輪軸間距，而對(duì)間距很小(兩輪胎緊挨)的車(chē)輪軸間距識(shí)別誤差較大。相比之下，車(chē)速監(jiān)測(cè)的準(zhǔn)確度相對(duì)較高。識(shí)別精度也與實(shí)測(cè)信號(hào)的質(zhì)量有關(guān)，不同傳感器其信號(hào)識(shí)別效果不同，建議選擇質(zhì)量較高的傳感器分別識(shí)別之后進(jìn)行平均，以獲得更好的識(shí)別結(jié)果。

圖4 第一跨跨中位置前/后應(yīng)變片信號(hào)近似部分Fig.4 Approximations for two strain signals before/after mid-span of first span

3 結(jié) 語(yǔ)

本文利用小波變換，提出一種基于橋面實(shí)測(cè)動(dòng)應(yīng)變數(shù)據(jù)來(lái)獲得橋梁車(chē)輛信息的識(shí)別方法，以實(shí)現(xiàn)對(duì)大跨度橋梁結(jié)構(gòu)上車(chē)輛行駛狀況的監(jiān)測(cè)。數(shù)值仿真結(jié)果表明:采用基于DWT的方法能夠較為準(zhǔn)確地得到橋上車(chē)輪軸間距和車(chē)速的信息，車(chē)速大小、噪聲水平以及所采用的小波種類(lèi)都會(huì)對(duì)識(shí)別結(jié)果產(chǎn)生影響。香港青衣南大橋現(xiàn)場(chǎng)測(cè)試分析表明，基于DWT的方法在實(shí)際大跨混凝土橋梁上行車(chē)監(jiān)測(cè)是可行的，識(shí)別得到的車(chē)速精度略高于車(chē)輪軸間距的精度，同時(shí)，車(chē)輪軸間距的大小對(duì)識(shí)別精度的影響也較大，從動(dòng)應(yīng)變信號(hào)中能夠較為準(zhǔn)確地識(shí)別前后輪距較大的車(chē)輪軸間距，對(duì)于標(biāo)準(zhǔn)車(chē)第2軸和第3軸車(chē)輪緊挨、軸間距很短的情況，識(shí)別精度較低，證明了采用基于DWT方法的有效性。

表3 由近似部分信號(hào)識(shí)別得到的車(chē)輪軸間距和車(chē)速Table 3 Vehicle axle spacing and speed identification from approximation signals

［1］朱勁松，肖汝誠(chéng)，何立志.大跨度斜拉橋智能可靠度評(píng)估方法研究［J］.土木工程學(xué)報(bào)，2007，40(5):41-48.(ZHU Jinsong，XIAO Rucheng，HE Lizhi.Reliability assessment of large-span cable-stayed bridges based on artificial intelligence［J］.China Civil Engineering Journal，2007，40(5):41-48.(in Chinese))

［2］MICHALTSOS G T.Dynamic behaviour of a single-span beam subjected to loads moving with variable speeds［J］.Journal of Sound and Vibration，2002，258(2):359-372.

［3］李小珍，張黎明，張潔.公路橋梁與車(chē)輛耦合振動(dòng)研究現(xiàn)狀與發(fā)展趨勢(shì)［J］.工程力學(xué)，2008，25(3):230-240.(LI Xiaozhen，ZHANG Liming，ZHANGJie.State-of-the-art review and trend of studies on coupling vibration for vehicle and highway bridge system［J］.Engineering Mechanics，2008，25(3):230-240.(in Chinese))

［4］張潔.公路車(chē)輛與橋梁耦合振動(dòng)分析研究［D］.成都:西南交通大學(xué)，2007.

［5］趙永剛，盛洪飛，陳明非.大跨度斜拉橋車(chē)-橋耦合振動(dòng)分析的數(shù)值解法［J］.公路工程，2009，34(1):143-146.(ZHAO Yonggang，SHENG Hongfei，CHEN Mingfei.Numerical method based on vehicle-bridge coupling vibration analysis of large span cable-stayed bridges［J］.Highway Engineering，2009，34(1):143-146.(in Chinese))

［6］肖新標(biāo)，沈火明.移動(dòng)荷載作用下橋梁的系統(tǒng)仿真［J］.振動(dòng)與沖擊，2005，24(1):121-123.(XIAO Xinbiao，SHEN Huoming.Dynamic simulation of bridge subjected to moving load［J］.Journal of Vibration and Shock，2005，24(1):121-123.(in Chinese))

［7］MOGHIMI H，RONAGH H R.Impact factors for a composite steel bridge non-linear dynamic simulation［J］.International Journal of Impact Engineering，2008，35:1228-1243.

［8］王解軍，張偉，吳衛(wèi)祥.重型汽車(chē)荷載作用下簡(jiǎn)支梁橋的動(dòng)力反應(yīng)分析［J］.中南公路工程，2005，30(2):55-62.(WANGJiejun，ZHANG Wei，WU Weixiang.Analysis of dynamic responses of simply supported girder bridge under heavy moving vehicles［J］.Central South Highway Engineering，2005，30(2):55-62.(in Chinese))

［9］劉文淑，林咸志，薛濤，等.基于小波變換的鋼箱梁應(yīng)力監(jiān)測(cè)多尺度分析［J］.河海大學(xué)學(xué)報(bào):自然科學(xué)版，2010，38(4):407-410.(LIU Wenshu，LIN Xianzhi，XUE Tao，et al.Multi-scale analysis of stress monitoring of steel box girders based on wavelet transform［J］.Journal of Hohai University:Nature Science，2010，38(4):407-410.(in Chinese))

［10］LIU Zejia，JIANG Bin，TANG Liqun.Features of long-term health monitored strains of a bridge with wavelet analysis［J］.Theoretical and Applied Mechanics Letters，2011，1(5):051006-051009.

［11］CAPRIOI L，CIGADA A，RAVEGLIA D.Rail inspection in track maintenance:a benchmark between the wavelet approach and the more conventional fourier analysis［J］.Mechanical Systems and Signal Processing，2007，21:631-652.

［12］魏明果.實(shí)用小波分析［M］.北京:北京理工大學(xué)出版社，2005.