環(huán)境溫度對連續(xù)箱梁剛構(gòu)橋相對小波熵指標(biāo)的影響研究

夏樟華， 林友勤， 張景杭， 宗周紅

(1. 福州大學(xué) 土木工程學(xué)院，福州 350108；2. 福州大學(xué) 至誠學(xué)院，福州 350002；3.東南大學(xué) 土木工程學(xué)院，南京 210096)

環(huán)境溫度對連續(xù)箱梁剛構(gòu)橋相對小波熵指標(biāo)的影響研究

夏樟華1， 林友勤1， 張景杭2， 宗周紅3

(1. 福州大學(xué) 土木工程學(xué)院，福州 350108；2. 福州大學(xué) 至誠學(xué)院，福州 350002；3.東南大學(xué) 土木工程學(xué)院，南京 210096)

以一座高速公路預(yù)應(yīng)力混凝土連續(xù)箱梁剛構(gòu)橋?yàn)楣こ瘫尘埃Y(jié)合健康監(jiān)測系統(tǒng)的長期測試數(shù)據(jù)，研究環(huán)境溫度對相對箱梁橋小波熵指標(biāo)SWT的影響規(guī)律。首先，介紹了環(huán)境激勵下橋梁結(jié)構(gòu)相對小波熵指標(biāo)的計(jì)算方法。其次，利用實(shí)時(shí)監(jiān)測數(shù)據(jù)，分析得到了箱梁橋各截面的實(shí)測相對小波熵指標(biāo)日變化和年變化規(guī)律。最后，基于逐步最優(yōu)線性回歸分析方法建立了相對小波熵指標(biāo)與溫度時(shí)、空間分布的多元線性回歸模型，并提出了相對小波熵指標(biāo)的溫度修正模型。結(jié)果表明，所建立多元線性回歸模型具有較好的擬合和預(yù)測精度，通過溫度修正可以有效地剔除環(huán)境溫度變化對預(yù)應(yīng)力混凝土箱梁橋相對小波熵指標(biāo)的影響，使其有效地應(yīng)用于橋梁的損傷識別。

預(yù)應(yīng)力混凝土連續(xù)箱梁剛構(gòu)橋；環(huán)境溫度；相對小波熵指標(biāo)；多元線性回歸模型；溫度修正模型

在結(jié)構(gòu)的損傷識別和健康診斷中，常采用動力參數(shù)作為損傷指標(biāo)。其中，小波包能量譜具有對結(jié)構(gòu)損傷敏感，抗噪聲能力強(qiáng)的特點(diǎn)。因此，利用小波包能量譜建立的損傷指標(biāo)在結(jié)構(gòu)或構(gòu)件損傷識別和健康監(jiān)測中得到了比較廣泛的應(yīng)用，包括小波熵[1-2]，相對小波熵[3]，相對能量譜[4-5]，小波包能量曲率差[6]等。在實(shí)橋應(yīng)用中，這些損傷指標(biāo)一般都需要結(jié)合環(huán)境激勵下的動力特性建立，不可避免地受到溫度的影響，可能導(dǎo)致實(shí)際結(jié)構(gòu)損傷引起的指標(biāo)變化被溫度影響所掩蓋。但目前對動力損傷指標(biāo)的溫度影響規(guī)律研究主要還是針對于橋梁結(jié)構(gòu)的模態(tài)頻率[7-12]。環(huán)境溫度對小波熵能量譜相關(guān)損傷指標(biāo)的影響研究較少，鄧揚(yáng)等[13]通過建立季節(jié)相關(guān)性模型得到了溫度對懸索橋小波包能量比影響的定量評價(jià)。對于預(yù)應(yīng)力混凝土連續(xù)箱梁橋，溫度對小波包能量譜損傷指標(biāo)的影響還研究不足，導(dǎo)致在其損傷識別和性能評估中還無法剔除溫度的影響，可能導(dǎo)致?lián)p傷誤判。因此，利用預(yù)應(yīng)力混凝土連續(xù)箱梁剛構(gòu)橋長期監(jiān)測的數(shù)據(jù)，進(jìn)行相對小波熵指標(biāo)的溫度影響規(guī)律研究，并建立相對小波熵的溫度修正模型，使相對小波熵指標(biāo)更好地用于預(yù)應(yīng)力混凝土連續(xù)箱梁剛構(gòu)橋的長期監(jiān)測、損傷判斷和性能評定中。

1 基于小波包能量譜的相對小波熵?fù)p傷指標(biāo)

采用小波包能量譜中各特征頻帶能量相對于所有頻帶能量總和的比值作為結(jié)構(gòu)損傷識別參數(shù)，能更為敏感發(fā)現(xiàn)結(jié)構(gòu)的損傷。采用各損傷特征頻帶的能量相對于所有頻帶能量總和的比值—相對小波能量pj作為結(jié)構(gòu)損傷識別參數(shù)，并按下式計(jì)算[5]：

(1)

式中:Ej為小波包頻帶能量譜和小波包時(shí)頻能量譜的第j個損傷特征頻帶的能量；E為小波包頻帶能量譜和小波包時(shí)頻能量譜所有頻帶能量總和；m為損傷特征頻帶的數(shù)量。通過相對小波能量pj的變化可以判斷結(jié)構(gòu)的損傷情況。為此將各損傷特征頻帶的相對小波能量變化組成結(jié)構(gòu)損傷特征向量Sj，即：

Sj=pjlg(pj/qj)(j=1,2,…，m)

(2)

式中:Sj為第j個特征頻帶的相對小波熵，pj和qj分別為結(jié)構(gòu)在損傷和健康狀態(tài)下第j個特征頻帶的相對小波能量。合理的損傷特征頻帶數(shù)量m的取值方法是：將前m個較大能量系數(shù)頻帶的動力響應(yīng)進(jìn)行小波包重構(gòu)，得到近似重構(gòu)響應(yīng)f′。計(jì)算原始f和重構(gòu)響應(yīng)f′的頻帶相對累積能量。頻帶相對累積能量ε由下式計(jì)算：

(3)

工程實(shí)際應(yīng)用時(shí)，可以選擇若干階頻帶計(jì)算其原始信號和重構(gòu)信號的頻帶相對累積能量，從而綜合考慮以確定特征頻帶數(shù)m的取值。

定義基于小波包能量譜的結(jié)構(gòu)損傷識別指標(biāo)-相對小波熵SWT如下：

(4)

為了有效地克服環(huán)境激勵的隨機(jī)性和不確定性，從而具有更好的環(huán)境激勵魯棒性。采用虛擬脈沖響應(yīng)函數(shù)[14]，將其與基于小波包能量譜的相對小波熵?fù)p傷識別方法相結(jié)合，應(yīng)用于環(huán)境激勵下的損傷識別指標(biāo)相對小波熵SWT的計(jì)算，基本流程如圖1所示。

圖1 相對小波熵SWT的計(jì)算流程圖Fig.1 Calculation flewt forrelative wavelet entropy index SWT

2 箱梁橋相對小波熵指標(biāo)分析

2.1 健康監(jiān)測系統(tǒng)介紹

該橋位于沈海高速福鼎至寧德段，是一座特大型的四跨預(yù)應(yīng)力混凝土連續(xù)箱梁剛構(gòu)橋。該橋上部結(jié)構(gòu)為跨徑145 m+2×260 m+l45 m，橋梁全長810 m，橋面寬為24.50 m，分為上下行兩幅，每幅橋箱梁均采用單箱單室截面。該橋設(shè)計(jì)荷載為汽超-20、掛-120，無人群荷載，設(shè)計(jì)時(shí)速為80 km/h，于2003年7月建成通車。在線實(shí)時(shí)健康監(jiān)測系統(tǒng)包括三部分：①靜態(tài)監(jiān)測系統(tǒng)。靜態(tài)監(jiān)測系統(tǒng)主要是采集和處理大橋的環(huán)境溫度與箱內(nèi)溫度以及大橋兩端的拉線位移等參數(shù)。②動態(tài)監(jiān)測系統(tǒng)。動態(tài)監(jiān)測系統(tǒng)主要是采集和處理大橋的加速度，動位移，應(yīng)變應(yīng)力等參數(shù)。③中央控制系統(tǒng)。中央控制系統(tǒng)主要負(fù)責(zé)數(shù)據(jù)的采集、數(shù)據(jù)的傳輸和系統(tǒng)的運(yùn)行控制，分為上位機(jī)顯示計(jì)算機(jī)、傳輸模塊、下位機(jī)采樣計(jì)算機(jī)以及遠(yuǎn)程控制模塊等四部分。該健康監(jiān)測系統(tǒng)監(jiān)測項(xiàng)目布置見圖2。

2.2 動力響應(yīng)監(jiān)測

在支座、邊跨1/2處、墩頂、主跨1/4處以及主跨1/2處共13個測試斷面布置了加速度傳感器，各個斷面編號從A～M，如圖2所示，每個斷面均包含豎向與橫向加速度傳感器，僅在A、B、C、E、F斷面上布置了縱向加速度傳感器，采樣頻率為200 Hz，24小時(shí)在線實(shí)時(shí)監(jiān)測。該橋的加速度傳感器采用的是中國地震局工程力學(xué)研究所的941BA型加速度傳感器，其基本參數(shù)如表1所示。加速度測點(diǎn)的橫向布置如圖3所示。

表1 加速度傳感器基本參數(shù)

圖2 大橋監(jiān)測項(xiàng)目布置圖(m)Fig.2 monitoring items arrangement of Xiabaishi bridge(m)

2.3 溫度監(jiān)測

2.3.1 測點(diǎn)布置

考慮縱、橫、豎三向溫度場分布情況，大橋縱橋向在C、E兩個斷面位置布置溫度傳感器，縱向溫度傳感器布置截面見圖2，豎、橫向布置見圖4。全橋共布設(shè)16個溫度傳感器，每個測試斷面各有8個，各溫度測點(diǎn)編號分別為C1～C8和E1～E8(分別表示C斷面1#～8#測點(diǎn)和E斷面1#～8#測點(diǎn))。通過靜態(tài)采集系統(tǒng)，每五分鐘采集一次溫度數(shù)據(jù)。該橋的環(huán)境溫度傳感器采用美國基康的BGK-3700溫度傳感器，溫度量測范圍為-30℃～70℃，精度為0.1℃，25℃阻值為3 kΩ。

圖3 加速度傳感器布置圖Fig.3 accelerators arrangement

圖4 溫度傳感器布置圖Fig.4 arrangment of temperature sensors

2.3.2 箱梁溫度分布情況

混凝土箱梁截面在受日照輻射時(shí)與外界的熱交換作用主要發(fā)生在截面的豎直方向，縱橋向溫度分布是比較類似的。在分析梁溫度場分布時(shí)，可以將三維溫度場簡化為二維溫度場來處理，這在滿足工程需要及精度方面是完全可行的。橋梁的豎向溫度測點(diǎn)可以簡化為橋面溫度、箱內(nèi)平均溫度、箱外平均溫度(箱外腹板與箱梁底板下側(cè)溫度平均值)。因此，以大橋E斷面溫度場作為分析對象來研究環(huán)境溫度與相對小波熵指標(biāo)兩者之間的關(guān)系。

2.4 相對小波熵指標(biāo)變化規(guī)律分析

2.4.1 相對小波熵指標(biāo)的計(jì)算

取A、B、C、D、E五個測試截面的豎向加速度信號來分析大橋的相對小波熵指標(biāo)SWT。由于該大橋在支座截面、邊跨1/2截面、墩頂截面、主跨1/4截面以及主跨1/2截面處對稱布置豎向加速度傳感器。因此，A、B、C、D、E五個測試截面的相對小波熵指標(biāo)SWT能夠反映全橋的相對小波熵指標(biāo)SWT。環(huán)境激勵下，對A、B、C、D、E截面的豎向加速度響應(yīng)信號進(jìn)行小波包能量譜分析時(shí)，首先要計(jì)算虛擬脈沖響應(yīng)函數(shù)，取福鼎側(cè)邊跨的L截面的加速度響應(yīng)信號作為參考點(diǎn)計(jì)算虛擬脈沖響應(yīng)函數(shù)。為了得到相對小波熵指標(biāo)，還需要確定橋梁完好結(jié)構(gòu)的基準(zhǔn)值?；鶞?zhǔn)日期的選擇要考慮兩個原則：①基準(zhǔn)日期的溫度要與該大橋合攏溫度相近；②基準(zhǔn)日期應(yīng)盡量接近該大橋的竣工日期，以確保基準(zhǔn)日期橋梁結(jié)構(gòu)為完好結(jié)構(gòu)。因?yàn)?007年5月13日的環(huán)境平均溫度為19.8℃，與該大橋合攏溫度20℃接近，而且也比較接近大橋的竣工日期，因此以2007年5月13日作為該大橋的基準(zhǔn)日期。

2.4.2 相對小波熵指標(biāo)的日變化和年變化規(guī)律

分別選取2013年11月7日(秋季)、2014年1月30日(冬季)、2014年3月6日(春季)和2014年6月27日(夏季)四天所采集的大橋A、B、C、D、E截面與福鼎側(cè)邊跨的L截面的豎向動力響應(yīng)信號進(jìn)行小波包能量譜分析，并分析損傷識別指標(biāo)。根據(jù)分析，A截面(支座截面)與C截面(墩頂截面)的相對小波熵指標(biāo)SWT晝夜變化不明顯，B截面(邊跨1/2截面)、D截面(主跨1/4截面)與E截面(主跨1/2截面)的相對小波熵指標(biāo)SWT變化為晝夜起伏的變化特征。B、D、E截面的相對小波熵指標(biāo)SWT在2014-01-30(冬季)與2014-06-27(夏季)的時(shí)候較大，在2013-11-07(秋季)與2014-03-06(春季)的時(shí)候較小，主要是因?yàn)槎?、夏季兩天的環(huán)境溫度與基準(zhǔn)溫度之差較大，一天的環(huán)境溫度變化也較大，春、秋季兩天的環(huán)境溫度與基準(zhǔn)溫度之差較小，一天的環(huán)境溫度變化較平緩。以B斷面為例，見圖5。

圖5 B截面相對小波熵指標(biāo)日變化實(shí)測曲線Fig.5 Measured dailyvariability curves of relative wavelet entropy index in section B

采用一年中不同季節(jié)的40天(共920小時(shí))的數(shù)據(jù)進(jìn)行分析(注：因?yàn)楸O(jiān)測系統(tǒng)的B截面采集箱故障了很長一段時(shí)間，因此B截面的實(shí)測樣只取了690個小時(shí))。經(jīng)過分析，A、C截面的相對小波熵指標(biāo)SWT實(shí)測曲線呈平穩(wěn)變化，B、D、E截面的相對小波熵指標(biāo)SWT實(shí)測曲線隨著季節(jié)性溫度變化出現(xiàn)比較大的起伏。B、D、E截面的相對小波熵指標(biāo)SWT實(shí)測曲線在冬、夏兩季的起伏變化較大，而在春、秋兩季的起伏變化較小。總得來說，該大橋的相對小波熵指標(biāo)SWT年變化規(guī)律基本上與日變化規(guī)律得出的結(jié)論相同。以D、E斷面為例，見圖6和圖7。

圖6 D截面相對小波熵指標(biāo)年變化實(shí)測曲線Fig.6 measured annualvariability curve of relative wavelet entropy index in section D

圖7 E截面相對小波熵指標(biāo)年變化實(shí)測曲線Fig.7 measured annualvariability curve of relative wavelet entropy index in section E

根據(jù)分析，B、D、E三個截面的相對小波熵指標(biāo)SWT的日變化與年變化規(guī)律都顯示了相對小波熵指標(biāo)SWT的變化與環(huán)境溫度存在相關(guān)關(guān)系，因此對這三個截面進(jìn)行相關(guān)關(guān)系分析。

3 相對小波熵指標(biāo)與環(huán)境溫度變化最優(yōu)線性回歸模型

3.1 最優(yōu)回歸模型

所謂“最優(yōu)線性回歸方程”，一方面是方程中包含所有對因變量有顯著影響的自變量，另一方面是方程所包含的自變量的個數(shù)要盡可能少。采用逐步回歸法(Stepwise)是向前回歸法(Forward)與向后回歸法(Backward)的結(jié)合，選擇自變量。

由于混凝土的導(dǎo)熱性能差，方向性強(qiáng)，對時(shí)間而言有明顯的滯后性，這就會導(dǎo)致橋梁動力特性變化往往滯后于環(huán)境溫度變化，基于橋梁動力特性變化的相對小波熵指標(biāo)SWT變化也往往會滯后于環(huán)境溫度的變化。因此，分析環(huán)境溫度變化對相對小波熵指標(biāo)SWT變化的影響時(shí)，不僅需要考慮環(huán)境溫度的空間分布，還應(yīng)該考慮時(shí)間分布。為了得到該大橋各個箱型截面的相對小波熵指標(biāo)SWT與各溫度測點(diǎn)與基準(zhǔn)溫度差值絕對值的溫度多元線性回歸模型，從每小時(shí)加速度信號中識別得到的各個截面的相對小波熵指標(biāo)SWT記為該時(shí)間段上的相對小波熵指標(biāo)SWT。同時(shí)，將截面各個溫度測點(diǎn)在同一時(shí)段上的平均值與基準(zhǔn)溫度差的絕對值作為該時(shí)段上的溫度記錄值，由此得到相對小波熵指標(biāo)SWT和溫度的同步時(shí)間序列：

(5)

由此可建立截面相對小波熵指標(biāo)SWT與溫度之間多元線性回歸方程為：

(6)

3.2 各個截面的最優(yōu)模型建立

逐步篩選法過程中以回歸系數(shù)顯著性檢驗(yàn)中的各自變量的F統(tǒng)計(jì)量的相伴概率值p作為自變量是否可以進(jìn)入方程或者是否可以剔除出方程的標(biāo)準(zhǔn)。在進(jìn)行逐步回歸分析前設(shè)定Entry值為0.05和Remove值為0.10。即如果一個自變量的F統(tǒng)計(jì)量的相伴概率值p小于顯著性水平0.05，則應(yīng)拒絕原假設(shè)H0，認(rèn)為該變量對因變量影響是顯著的；如果方程中一個自變量的F統(tǒng)計(jì)量的相伴概率值p大于顯著性水平0.10，則不能拒絕原假設(shè)H0，它應(yīng)該剔除出回歸方程。

由式(6)可以看出，回歸方程的自變量因子中同時(shí)含有溫度的空間分布和時(shí)間分布，為解決預(yù)應(yīng)力混凝土箱梁橋的相對小波熵指標(biāo)SWT與環(huán)境溫度變化之間時(shí)間滯后問題，考慮將時(shí)間滯后1小時(shí)、2小時(shí)和3小時(shí)的實(shí)測溫度數(shù)據(jù)分別進(jìn)行逐步最優(yōu)回歸，分別得到相應(yīng)的回歸方程和回歸結(jié)果。B、E截面的相對小波熵指標(biāo)SWT與環(huán)境溫度之間確實(shí)存在時(shí)間滯后問題，而且與時(shí)間滯后2小時(shí)的溫度數(shù)據(jù)擬合效果最佳，為最優(yōu)回歸方程。D截面的相對小波熵指標(biāo)SWT與環(huán)境溫度之間不存在時(shí)間滯后問題。最終所得的B、D、E截面的最優(yōu)回歸方程和檢驗(yàn)結(jié)果如表2～表4所示。

表2 B截面回歸方程及檢驗(yàn)

表3 D截面回歸方程及檢驗(yàn)

表4 E截面回歸方程及檢驗(yàn)

對多元回歸模型而言，經(jīng)過顯著性檢驗(yàn)，仍然不能保證方程中所有自變量都能較好反映因變量變化。因此，通過回歸系數(shù)顯著性檢驗(yàn)對每個系數(shù)進(jìn)行分析，即進(jìn)行偏相關(guān)的F檢驗(yàn)，p值用來檢驗(yàn)系數(shù)的顯著性水平。經(jīng)回歸系數(shù)檢驗(yàn)可知，各截面的各變量對其回歸方程均有顯著的線性關(guān)系。以B截面為例，見表5。

表5 B截面回歸系數(shù)及檢驗(yàn)

3.3 回代檢驗(yàn)

除了用給出的檢驗(yàn)統(tǒng)計(jì)量來檢驗(yàn)回歸方程的合理性之外，還必須對回歸方程做回代檢驗(yàn)?；卮鷻z驗(yàn)的方法就是把原始數(shù)據(jù)重新代到回歸方程來驗(yàn)證回歸結(jié)果與實(shí)測結(jié)果的擬合性能。

定義模型的擬合殘差為：

(7)

再定義均方誤差σ(x)作為誤差評判指標(biāo)：

(8)

圖8 D截面相對小波熵指標(biāo)的多元線性回歸模型的擬合殘差及其分布Fig.8 Fitting residual and its distribution for multiple linear regression model of relative wavelet entropy index in section D

誤差類型相對小波熵指標(biāo)SWT×10-2B截面D截面E截面擬合均方誤差1.481.672.43預(yù)測均方誤差1.451.312.56

3.4 外推檢驗(yàn)

圖9 E截面相對小波熵指標(biāo)的多元線性模型泛化性能分析Fig.9 Generalization performance analysis of multiple linear model of relative wavelet entropy index in section E

4 相對小波熵指標(biāo)的溫度修正模型

通過統(tǒng)計(jì)回歸分析建立起了“相對小波熵指標(biāo)SWT與溫度空間、時(shí)間分布”的多元線性回歸模型，考慮消除環(huán)境溫度變化對相對小波熵指標(biāo)的影響，提出相對小波熵指標(biāo)的溫度修正模型為：

(9)

將各截面回歸模型建立的樣本數(shù)據(jù)點(diǎn)的實(shí)測相對小波熵指標(biāo)值和由回歸模型計(jì)算得到的擬合值代入溫度修正模型中，得到剔除環(huán)境溫度變化影響后的相對小波熵指標(biāo)。消除環(huán)境溫度變化影響后，B、D、E截面的相對小波熵指標(biāo)不存在隨日溫度變化和季節(jié)性溫度變化而顯著波動的現(xiàn)象，指標(biāo)變化趨于平穩(wěn)，年均值基本都趨近于0，見表7。

表7 不同截面修正后的相對小波熵指標(biāo)的年統(tǒng)計(jì)參數(shù)

5 結(jié) 論

(1)基于預(yù)應(yīng)力混凝土連續(xù)箱梁剛構(gòu)橋健康監(jiān)測系統(tǒng)，分析了環(huán)境激勵下相對小波熵指標(biāo)SWT與溫度之間的關(guān)系，支座截面與墩頂截面的環(huán)境溫度變化對相對小波熵指標(biāo)影響不明顯。邊跨1/2截面、主跨1/4截面與主跨1/2截面的相對小波熵指標(biāo)與環(huán)境溫度變化存在顯著的正相關(guān)關(guān)系。

(2)振動信號越強(qiáng)的箱型截面對相對小波熵指標(biāo)產(chǎn)生影響的溫度測點(diǎn)越多，在回歸模型中需要考慮的自變量越多。其中，主跨1/2截面回歸方程有5個自變量，邊跨1/2截面回歸方程有4個自變量，主跨1/4截面回歸方程有2個自變量。

(3)利用逐步最優(yōu)線性回歸分析方法分析預(yù)應(yīng)力混凝土連續(xù)箱梁橋相對小波熵指標(biāo)與溫度空間分布、時(shí)間分布的關(guān)系，得到了主跨1/4截面相對小波熵指標(biāo)與溫度空間分布的多元線性最優(yōu)回歸模型，邊跨1/2截面與E主跨1/2截面的相對小波熵指標(biāo)與溫度、時(shí)間分布的多元線性最優(yōu)回歸模型。通過回推檢驗(yàn)發(fā)現(xiàn)，所建立的回歸模型具有很好的預(yù)測精度。

(4)基于最優(yōu)回歸模型建立了連續(xù)箱梁橋相對小波熵指標(biāo)溫度修正模型，該模型能較好地消除環(huán)境溫度變化引起的相對小波熵指標(biāo)變化。通過溫度修正后，邊跨1/2截面、主跨1/4截面、主跨1/2截面相對小波熵指標(biāo)已經(jīng)不存在隨季節(jié)性溫度變化而波動的現(xiàn)象，總體變化趨于平穩(wěn)。因此，修正后的相對小波熵指標(biāo)已經(jīng)基本剔除溫度的影響，可以有效地應(yīng)用于預(yù)應(yīng)力混凝土箱梁橋的損傷識別和性能評定中。

[ 1 ] REN W X, SUN Z S, XIA Y, et al. Damage identification of shear connectors with wavelet packet energy: laboratory test study[J]. Journal of Structural Engineering, 2008, 134(5): 832-841.

[ 2 ] REN W X, SUN Z S. Structural damage identification by using wavelet entropy[J]. Engineering Structures, 2008, 30(10): 2840-2849.

[ 3 ] 孫增壽, 任偉新. 基于小波熵指標(biāo)的結(jié)構(gòu)損傷檢測[J]. 振動. 測試與診斷, 2008, 28(3): 233-237. SUN Zengshou,REN Weixin. Structure damage detection based on wavelet entropy indexes[J]. Journal of Vibration, Measurement &Diagnosis, 2008, 28(3): 233-237.

[ 4 ] 劉濤, 李愛群, 丁幼亮, 等. 基于小波包能量譜的結(jié)構(gòu)損傷預(yù)警方法試驗(yàn)研究[J]. 振動與沖擊, 2009, 28(4): 4-8. LIU Tao, LI Aiqun, DING Youliang, et al. Experimental study on structural damage alarming method based on wavelet packet energy spectrum[J]. Journal of Vibration and Shock, 2009, 28(4): 4-8.

[ 5 ] 丁幼亮, 李愛群, 鄧揚(yáng). 面向結(jié)構(gòu)損傷預(yù)警的小波包能量譜識別參數(shù)[J]. 東南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 41(4): 824-828. DING Youliang, LI Aiqun, DENG Yang. Parameters for identification of wavelet packet energy spectrum for structural damage alarming[J]. Journal of southeast university (Natural Science Edition), 2011, 41(4):824-828.

[ 6 ] 余竹, 夏禾, GOICOLEA J M, 等. 基于小波包能量曲率差法的橋梁損傷識別試驗(yàn)研究[J]. 振動與沖擊, 2013, 32(5): 20-25. YU Zhu, XIA He, GOICOLEA J M, et al. Experimental study on bridge damage identification based on wavelet packet energy curvature difference method [J]. Journal of Vibration and Shock, 2013, 32(5): 20-25.

[ 7 ] SOHN H, DZWONCZYK M, STRASER E G, et al. An experimental study of temperature effect on modal parameters of the Alamosa Canyon Bridge[J]. Earthquake engineering & structural dynamics, 1999, 28(8): 879-897.

[ 8 ] ZHAO J, DEWOLF J T. Dynamic monitoring of steel girder highway bridge[J]. Journal of Bridge Engineering, 2002, 7(6): 350-356.

[ 9 ] NI Y Q, HUA X G, FAN K Q, et al. Correlating modal properties with temperature using long-term monitoring data and support vector machine technique[J]. Engineering Structures, 2005, 27(12): 1762-1773.

[10] LIU C, DEWOLF J T. Effect of temperature on modal variability of a curved concrete bridge under ambient loads [J]. Journal of Structural Engineering, 2007, 133(12): 1742-1751.

[11] 丁幼亮, 李愛群, 耿方方. 考慮環(huán)境因素影響的懸索橋整體狀態(tài)預(yù)警方法[J]. 東南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010, 40(5): 1052-1056. DING Youliang, LI Aiqun, GENG Fangfang. Monitoring and warning of health conditions for suspension bridges under varying environmental conditions[J]. Journal of southeast university (Natural Science Edition), 2010, 40(5): 1052-1056.

[12] 余印根, 宗周紅, 陳寶春, 等. 環(huán)境溫度對連續(xù)剛構(gòu)橋模態(tài)頻率的影響[J]. 振動. 測試與診斷, 2014, 34(1): 69-76. YU Yingen, ZONG Zhouhong, CHEN Baochun, et al. Effects of environmental temperature on modal frequency of continuous rigid frame bridge[J]. Journal of Vibration, Measurement &Diagnosis, 2014, 34(1): 69-76.

[13] 鄧揚(yáng), 李愛群, 丁幼亮. 環(huán)境激勵下大跨懸索橋?qū)崪y損傷預(yù)警指標(biāo)的變異性[J]. 防災(zāi)減災(zāi)工程學(xué)報(bào), 2013, 33(6): 624-630. DENG Yang, LI Aiqun, DING Youliang. Variability of measured damage alarming index of a long-span suspension bridge under ambient excitation[J]. Journal of Disaster Prevention and Mitigation Engineering, 2013, 33(6): 624-630.

[14] 丁幼亮, 李愛群, 繆長青. 基于小波包能量譜的結(jié)構(gòu)損傷預(yù)警方法研究[J]. 工程力學(xué), 2006, 23 (8): 42-48. DING Youliang, LI Aiqun, MIAO Changqing. Investigation on the structural damage alarming method based on wavelet packet energy spectrum [J]. Engineering Mechanics, 2006, 23 (8): 42-48.

Influences of environmental temperature on relative wavelet entropy index for continuous rigid-frame bridges with box sections

XIA Zhanghua1, LIN Youqin1, ZHANG Jinhang2, ZONG Zhouhong3

(1. College of Civil Engineering, Fuzhou University, Fuzhou 350108, China;2. Zhicheng College, Fuzhou University, Fuzhou 350002, China;3. College of Civil Engineering, Southeast University, Nanjing 210096, China)

In order to obtain influence laws of environmental temperature on the relative wavelet entropy index SWT, the related study based on health monitoring system of a prestressed concrete continuous rigid-frame bridge with box sections was conducted. Firstly, the calculation method of the relative wavelet entropy index for bridge structures under environment excitation was introduced. Secondly, the daily and annual varying rules of the relative wavelet entropy index of each box section of the bridge were analyzed and obtained based on the real-time monitoring data. Finally, by using the gradually optimal linear regression analysis method, the multivariate linear regression model for the relative wavelet entropy index considering time and spatial distribution of temperature was developed, and the temperature correctional model for the relative wavelet entropy index was proposed. Results showed that the multivariate linear regression model has a better fitting accuracy and a relatively high prediction precision; the effects of environmental temperature on the relative wavelet entropy index can be eliminated effectively with the temperature correctional model to make the relative wavelet entripy index be used in bridges’ damage identification effectively.

prestressed concrete continuous rigid-frame bridges with box sections；environmental temperature； relative wavelet entropy index；multivariate linear regression model；temperature correctional model

國家自然科學(xué)基金資助項(xiàng)目(51178101)

2015-09-29 修改稿收到日期：2016-02-25

夏樟華 男，博士,副研究員，1980年生

林友勤 男，博士,高級實(shí)驗(yàn)師，1973年生 E-mail:lyq@fzu.edu.cn

TU311.3

A

10.13465/j.cnki.jvs.2017.01.031