潘楚東 余嶺 劉煥林

摘要： 利用橋梁實(shí)測響應(yīng)反演橋面移動(dòng)荷載是橋梁工程的研究熱點(diǎn)之一。但既有移動(dòng)荷載正則化識(shí)別方法沒有考慮未知初始條件的影響。針對此問題，基于稀疏正則化技術(shù)提出了一種新的移動(dòng)荷載識(shí)別方法。首先分別采用冗余字典與模態(tài)空間描述移動(dòng)荷載與初始條件，計(jì)算各個(gè)荷載與初始條件分量引起的響應(yīng)，歸一化構(gòu)造用于表達(dá)響應(yīng)的像字典。然后利用像字典對響應(yīng)進(jìn)行稀疏分解，其結(jié)果即可用于估算未知移動(dòng)荷載。最后采用數(shù)值仿真算例評價(jià)所提方法的可行性與有效性，并通過實(shí)驗(yàn)室搭建的車-橋系統(tǒng)模型對新方法進(jìn)行實(shí)驗(yàn)驗(yàn)證。研究結(jié)果表明：未知初始條件下，所提方法能有效地識(shí)別移動(dòng)荷載。

關(guān)鍵詞： 移動(dòng)荷載識(shí)別； 橋梁； 初始條件； 稀疏正則化； 稀疏分解

中圖分類號： TU311.3； U441+.2 文獻(xiàn)標(biāo)志碼：A 文章編號1004-4523（2018）05-0734-10

DOI：10.16385/j.cnki.issn.1004-4523.2018.05.002

引 言

橋梁移動(dòng)荷載是橋梁主要活荷載之一，它影響著橋梁的運(yùn)營安全與使用壽命。自上世紀(jì)90年代以來，有關(guān)移動(dòng)荷載識(shí)別的研究受到廣泛關(guān)注[1-2]。早期研究主要集中在數(shù)學(xué)模型的建立上，代表性的方法有第一識(shí)別法（IMⅠ），第二識(shí)別法（IMⅡ），時(shí)域法（TDM）與頻-時(shí)域法（FTDM）。Yu和Chan對上述4種方法進(jìn)行了比較，研究表明上述4種方法均能有效地識(shí)別移動(dòng)荷載[3]。然而，由于移動(dòng)荷載識(shí)別“逆問題”特點(diǎn)，因此直接求解容易受測量噪聲的影響。

近廿年來，移動(dòng)荷載識(shí)別研究得到進(jìn)一步發(fā)展，主要表現(xiàn)為求解算法的多樣性。Law與Zhu提出基于Tikhonov正則的移動(dòng)荷載識(shí)別方法，旨在提高經(jīng)典方法的抗噪性[4-5]；Feng等結(jié)合貝葉斯正則化理論，研究移動(dòng)荷載識(shí)別與結(jié)構(gòu)參數(shù)估計(jì)[6-7]。以上工作可歸類于基于正則化技術(shù)的動(dòng)荷載識(shí)別研究[8-10]。正則化技術(shù)在求解時(shí)會(huì)引入荷載的合理約束，因此其識(shí)別結(jié)果具有較強(qiáng)的抗噪性。除了上述方法，研究人員針對不同工程模型簡化條件，又提出了許多求解方法，例如：影響線法[11-12]，將車輪接觸力簡化為移動(dòng)常力；形函數(shù)法[13]利用少數(shù)已知基函數(shù)近似移動(dòng)荷載；小波有限元法[14-15]采用小波有限元的車-橋系統(tǒng)模型。這些新方法的引入使得更準(zhǔn)確識(shí)別移動(dòng)荷載成為可能。此外，Cai所帶領(lǐng)的研究團(tuán)隊(duì)在此階段做了大量工作，提出了不少移動(dòng)荷載識(shí)別新方法，對該領(lǐng)域的發(fā)展做出了重要貢獻(xiàn)[16]。

一定程度上，既有方法能很好地識(shí)別移動(dòng)荷載，但其理論通常是建立在結(jié)構(gòu)初始條件已知的條件下，應(yīng)用要求嚴(yán)格。實(shí)際工程中，由于車流、環(huán)境、橋墩連接方式等因素的影響，在車輛上橋前，結(jié)構(gòu)通常處于振動(dòng)狀態(tài)。換言之，橋梁的初始條件是未知的。因此有必要研究未知初始條件下的移動(dòng)荷載識(shí)別問題。

基于上述原因，結(jié)合稀疏正則化，本文提出了一種未知初始條件下移動(dòng)荷載識(shí)別新方法。首先利用冗余字典展開未知移動(dòng)荷載并采用模態(tài)空間表示未知初始條件。其次求解相關(guān)像字典，并利用像字典分解實(shí)測響應(yīng)；最后利用響應(yīng)分解結(jié)果估算移動(dòng)荷載。同時(shí)，采用數(shù)值仿真與模型實(shí)驗(yàn)對所提方法的可行性與有效性進(jìn)行驗(yàn)證。

1 基本理論

1.1 移動(dòng)荷載作用下的測點(diǎn)響應(yīng) 圖1所示為移動(dòng)荷載作用下的簡支梁模型。

以10%噪聲工況為例，響應(yīng)分解結(jié)果分布如圖2所示。從圖中可見，所提方法在此工況中能有效地區(qū)分由荷載正弦分量、荷載沖擊分量和初始條件分量引起的響應(yīng)。在稀疏正則化的作用下，響應(yīng)分解結(jié)果存在大量的0值分量，說明所提方法能稀疏地提取響應(yīng)的主要特征。此外，圖2（b）中只有第一階分量（索引為801與802）為非零元素，與實(shí)際工況模擬值符合。圖3對比了3個(gè)噪聲工況下的移動(dòng)荷載識(shí)別結(jié)果，從圖中可見，移動(dòng)荷載識(shí)別結(jié)果能很好地吻合真實(shí)荷載。說明在給定的3個(gè)工況下，本文所提方法能有效地處理未知（非零）初始條件下的移動(dòng)荷載識(shí)別問題，并且具有強(qiáng)噪聲魯棒性。實(shí)際上，兼顧初始條件是大多數(shù)既有移動(dòng)荷載識(shí)別方法所欠考慮的，屬本文的主要?jiǎng)?chuàng)新點(diǎn)之一。

3 實(shí)驗(yàn)驗(yàn)證

實(shí)驗(yàn)室搭建車橋模型進(jìn)行實(shí)驗(yàn)，如圖4所示。小車在電機(jī)牽引下通過主梁，光電傳感器用于監(jiān)測小車速度。主梁為3 m長空心鋼管梁，如圖5（a）所示。桁架用于固定光電傳感器，支座如圖5（d）所示。

小車模型如圖5（c）所示。實(shí)測前后軸軸重分別為4.5955 kg與6.0290 kg，軸距為0.42 m。為研究非零初始條件的影響，在小車上橋前，采用錘擊法迫使橋梁產(chǎn)生振動(dòng)。自制錘子如圖5（b）所示，其中錘頭采用泡沫材料，以便使錘擊能量更好地集中在低頻區(qū)段。實(shí)驗(yàn)錘擊位置選取為距主梁左端支座2 m處。

在主梁底面中軸線的1/5，2/5，3/5，4/5與1/2跨處安裝有應(yīng)變片（120 Ω）。同時(shí)在1/4與1/2跨處安裝加速度傳感器（PCB，ICP 333B30）。利用LMS Test.Lab系統(tǒng)及其配套數(shù)采儀采集響應(yīng)數(shù)據(jù)，設(shè)置采用頻率為2048 Hz。

3.1 有限元模型

采用激振器激勵(lì)對主梁進(jìn)行實(shí)驗(yàn)?zāi)B(tài)分析。提取實(shí)測前4階頻率與振型。如圖6所示，橋梁均分為20個(gè)4自由度梁單元。采用豎向彈簧支撐模擬實(shí)際支座。主梁實(shí)測線密度約為ρA=7.79116 kg/m。適當(dāng)調(diào)節(jié)模型抗彎剛度EI與彈簧彈性系數(shù)kv，并最終確選為EI=1.4411×105 N·m2，kv=7.99×106 N·m-1。

表1對比了實(shí)測頻率與模型計(jì)算頻率；圖7對比了實(shí)測振型與模型計(jì)算振型。從表1與圖7可以看出，模型前4階頻率與振型均能很好地吻合主梁實(shí)測值。上述結(jié)果不能保證模型的高階模態(tài)參數(shù)接近實(shí)際橋梁的高階模態(tài)信息。但實(shí)際工程中，結(jié)構(gòu)響應(yīng)往往以低階為主。因此，該有限元模型可用于移動(dòng)荷載識(shí)別研究。分析時(shí)取前3階模態(tài)參數(shù)，利用振型疊加法計(jì)算結(jié)構(gòu)響應(yīng)。時(shí)間間隔為（1/2048） s。

3.2 靜力加載標(biāo)定應(yīng)變片

采用分級靜力加載對應(yīng)變片進(jìn)行標(biāo)定。加載所用砝碼如圖5（e）所示。表2列舉了在0.5 m處分級加載工況下的實(shí)測應(yīng)變值。其中彎矩增量的計(jì)算依據(jù)來源于結(jié)構(gòu)力學(xué)基本理論。從表2可看出，在0.5 m處加載時(shí)，跨中處的彎矩與應(yīng)變增量對為：（12.4950 N·m， 2.6983με）。依次更換加載位置（0.5 m，1 m，1.5 m，2 m，2.5 m），并對所有應(yīng)變片實(shí)行上述分析過程。提取所有工況下的彎矩與應(yīng)變增量對進(jìn)行線性擬合，結(jié)果如圖8所示。由圖8易得，應(yīng)變片靈敏度系數(shù)對應(yīng)于圖中直線斜率。

3.3 實(shí)驗(yàn)數(shù)據(jù)預(yù)分析

實(shí)驗(yàn)部分實(shí)測信號如圖9所示。采用截止頻率為200 Hz的Butterworth低通濾波器對結(jié)構(gòu)響應(yīng)進(jìn)行濾波。提取從小車上橋到小車下橋所對應(yīng)的響應(yīng)進(jìn)行后續(xù)識(shí)別分析。如表3所示，依據(jù)不同的車速與不同的錘擊力度，本文共整理9組實(shí)驗(yàn)工況用于移動(dòng)荷載識(shí)別。為了直觀地了解錘擊力度與初始條件大小的關(guān)系，圖10對比了不同錘擊力度下1/2跨測點(diǎn)的應(yīng)變響應(yīng)。從圖10中不難發(fā)現(xiàn)，在車輛上橋前，實(shí)測應(yīng)變幅值隨錘擊力度的加強(qiáng)而增大，符合實(shí)際工程經(jīng)驗(yàn)。

3.4 驗(yàn)證新方法的正確性和有效性

以工況F-3為例，考慮荷載所關(guān)心的最高頻率為200 Hz。采用1/5，2/5，3/5，4/5跨的彎矩響應(yīng)聯(lián)合1/4跨的加速度響應(yīng)識(shí)別移動(dòng)荷載。其余計(jì)算參數(shù)同數(shù)值仿真工況。

第一步：利用像字典對歸一化響應(yīng)進(jìn)行稀疏分解，結(jié)果如圖11所示。從圖11可得，在稀疏正則化的作用下，信號分解結(jié)果存在較多的0值分量。實(shí)測信號包含有由簡諧荷載、局部沖擊荷載與非零初始條件引起的分量。

第二步：計(jì)算縮放因子與信號幅值參考值，結(jié)合式（7）與（21）估算移動(dòng)荷載與結(jié)構(gòu)初始條件，結(jié)果如圖12所示。從圖12（a）與（b）可以看出，前后軸識(shí)別結(jié)果的趨勢項(xiàng)與真實(shí)軸重有較高的吻合度，說明所提方法在此工況中能有效地估算軸重。其次，前后軸識(shí)別結(jié)果的波動(dòng)相對較大。其主要原因在于，實(shí)驗(yàn)室主梁與小車車輪均采用硬質(zhì)材料，故在車速較快的情況下易產(chǎn)生較大的沖擊。此外，從圖中亦能看出除了少數(shù)有限個(gè)點(diǎn)的識(shí)別值小于0之外，絕大部分識(shí)別荷載值均在0值以上。該識(shí)別結(jié)果是比較合理的，因?yàn)楹奢d小于0意味著車輪與橋面之間存在拉力，不符合實(shí)際。圖12（c）表明在車輛上橋時(shí)，結(jié)構(gòu)初始條件不為0，即結(jié)構(gòu)處于非靜止?fàn)顟B(tài)。其主要原因是在車輛上橋前，主梁受到錘子沖擊而產(chǎn)生明顯的振動(dòng)。

（27）式中 Wtrue表示真實(shí)車重，Widentified表示估算車重。這里Widentified的計(jì)算方法為：取位于時(shí)間區(qū)間[t0.5m， t2.5m]內(nèi)前后軸識(shí)別結(jié)果和的平均值。其中t0.5m，t2.5m分別表示后軸位于主梁0.5 m時(shí)、前軸位于主梁2.5 m時(shí)所對應(yīng)的時(shí)刻。圖13給出工況F-3前后軸移動(dòng)荷載識(shí)別結(jié)果之和。從圖中可見，該工況下本文所提方法能有效地估算小車車重。圖14對比了跨中實(shí)測響應(yīng)與重構(gòu)響應(yīng)，其中重構(gòu)響應(yīng)采用上述移動(dòng)荷載識(shí)別結(jié)果進(jìn)行計(jì)算。從圖中可見，跨中實(shí)測響應(yīng)與重構(gòu)響應(yīng)具有較高的吻合度。該結(jié)果從響應(yīng)方面間接地表明了識(shí)別結(jié)果的合理性。

3.5 不同初始條件的影響

下面以快速組工況為例，分析不同初始條件對識(shí)別結(jié)果的影響。以工況F-2為例，圖15對比了考慮與不考慮初始條件下移動(dòng)荷載識(shí)別結(jié)果。其中，不考慮初始條件意味著求解時(shí)刪除系統(tǒng)矩陣中有關(guān)初始條件的列。從圖15可得，除前軸前端部分之外，考慮與不考慮初始條件對識(shí)別結(jié)果影響不大。

進(jìn)一步，以工況F-1，F(xiàn)-2與F-3為例，圖16對比了考慮與不考慮初始條件下前軸移動(dòng)荷載識(shí)別結(jié)果。從圖16（a）中可以看出，當(dāng)主梁沒受到錘擊時(shí)，識(shí)別精度不受初始條件的影響。因?yàn)橄鄬囍马憫?yīng)，初始條件引起的響應(yīng)足夠小。從圖16（b）與16（c）可以看出，當(dāng)主梁受到錘擊時(shí)，初始條件的影響不可忽略。為評估非零初始條件所帶來的影響，取不考慮初始條件時(shí)識(shí)別結(jié)果的最大幅值與車重（104.12 N）的比值進(jìn)行比較。從圖中可粗略得知，3個(gè)工況的比值約為：工況F-1，1.0倍；工況F-2，1.5倍；工況F-3，4.0倍。由此可見，不同大小的初始條件對識(shí)別結(jié)果所帶來的影響程度明顯不同。

圖17對比了初始條件識(shí)別結(jié)果。由圖17可看出，初始條件主要以一階為主，且幅值上基本呈現(xiàn)出F-1 3.6 不同車速下的計(jì)算效率

下面以工況F-3，M-3與S-3為例，討論不同車速對計(jì)算效率的影響。所采用的計(jì)算環(huán)境為：臺(tái)式計(jì)算機(jī)，4 CPU：Intel （R） Core（TM） i5-4590 @ 3.3 GHz，RAM：8G，Matlab version：R2015a。表4列出了計(jì)算詳細(xì)的CPU耗時(shí)。從表4可看出，主要的CPU耗時(shí)發(fā)生在信號分解階段。這是因?yàn)樾盘栂∈璺纸怆A段需要大量的迭代運(yùn)算，而荷載重構(gòu)階段不需要迭代運(yùn)算。其次，表4表明在上述3個(gè)工況中，CPU耗時(shí)隨著速度的下降而上升。主要是因?yàn)檐囁傧陆祵?dǎo)致車輛過橋時(shí)間變長，進(jìn)而使系統(tǒng)矩陣維數(shù)增大，影響計(jì)算效率。

3.7 所有工況識(shí)別精度對比

對所有工況進(jìn)行移動(dòng)荷載識(shí)別，并計(jì)算相應(yīng)的識(shí)別車重、1/2跨實(shí)測與重建彎矩的相關(guān)系數(shù)、1/2跨實(shí)測與重建加速度的相關(guān)系數(shù)，并將結(jié)果記錄在表5中。從表5可以看出，在所有9個(gè)工況下，所提新方法均能有效地識(shí)別小車車重；跨中實(shí)測響應(yīng)與重建響應(yīng)具有高相關(guān)性。上述結(jié)果從不同方面表明新方法能有效地處理含未知初始條件的移動(dòng)荷載識(shí)別問題，比既有方法（只處理已知初始條件問題）更具一般性，應(yīng)用條件更寬松。

4 結(jié) 論

針對既有移動(dòng)荷載識(shí)別方法較少考慮未知初始條件影響的現(xiàn)狀，本文提出了基于稀疏正則化的移動(dòng)荷載識(shí)別新方法。該方法首先利用像字典對實(shí)測信號進(jìn)行稀疏分解，然后利用稀疏分解結(jié)果估算移動(dòng)荷載。采用數(shù)值仿真與模型實(shí)驗(yàn)共同驗(yàn)證所提方法的有效性，綜合全文所得主要結(jié)論如下：

（1）基于稀疏正則化的移動(dòng)荷載識(shí)別方法能有效處理未知初始條件下的移動(dòng)荷載識(shí)別問題，更具有一般性。

（2）所提方法能有效地識(shí)別移動(dòng)荷載與結(jié)構(gòu)未知初始條件，亦能準(zhǔn)確地估算移動(dòng)小車車重，具有工程參考價(jià)值。

（3）在計(jì)算效率方面，所提方法計(jì)算耗時(shí)受小車車速影響較為顯著。在估算小車車重方面，所提方法識(shí)別精度受小車速度影響不明顯。

參考文獻(xiàn)：

[1] 陳 敏. 橋上移動(dòng)荷載識(shí)別問題的研究現(xiàn)狀[J]. 北方交通，2008， （1）： 126—127.

CHEN Min. Current situations of studying the distinguishing of moving load on bridge[J]. Northern Communications， 2008， （1）： 126—127.

[2] Zhu X Q， Law S S. Recent developments in inverse problems of vehicle-bridge interaction dynamics[J]. Journal of Civil Structural Health Monitoring， 2016， 6（1）： 107—128.

[3] Yu L， Chan T H T. Recent research on identification of moving loads on bridges[J]. Journal of Sound and Vibration， 2007， 305（1-2）： 3—21.

[4] Law S S， Chan T H T， Zhu Q X， et al. Regularization in moving force identification[J]. Journal of Engineering Mechanics， 2001， 127（2）： 136—148.

[5] Zhu X Q， Law S S. Moving loads identification through regularization[J]. Journal of Engineering Mechanics， 2002， 128（128）， 989—1000.

[6] Sun H， Feng D， Liu Y， et al. Statistical regularization for identification of structural parameters and external loadings using state space models[J]. Computer-Aided Civil and Infrastructure Engineering， 2015， 30（11）： 843—858.

[7] Feng D， Sun H， Feng M Q. Simultaneous identification of bridge structural parameters and vehicle loads[J]. Computers & Structures， 2015， 157： 76—88.

[8] 陳 震， 余 嶺. 基于截?cái)郍SVD方法的橋梁移動(dòng)荷載識(shí)別[J]. 振動(dòng)與沖擊， 2014， 33（10）：97—100.

CHEN Zhen， YU Ling. Identification of dynamic axle loads on bridge based on truncated generalized singular value decomposition[J]. Journal of Vibration and Shock， 2014， 33（10）：97—100.

[9] Bao Y Q， Li H， Chen Z C， et al. Sparse l1 optimization-based identification approach for the distribution of moving heavy vehicle loads on cable-stayed bridges[J]. Structural Control & Health Monitoring， 2016， 23（1）：144—155.

[10] Pan C D， Yu L， Liu H L. Identification of moving vehicle forces on bridge structures via moving average Tikhonov regularization[J]. Smart Materials & Structures， 2017， 26： 085041.

[11] 錢長照， 陳昌萍. 利用彎矩影響線的移動(dòng)荷載識(shí)別方法[J]. 動(dòng)力學(xué)與控制學(xué)報(bào)， 2016， 14（2）：182—185.

QIAN Changzhao， CHEN Changping. Method for moving force identification using bending moment influence line[J]. Journal of Dynamics and Control， 2016， 14（2）：182—185.

[12] 王寧波， 任偉新， 李 苗. 基于影響線的橋梁移動(dòng)荷載識(shí)別[J]. 振動(dòng)與沖擊， 2013， 32（3）：129—133.

WANG Ning-bo， REN Wei-xin， LI Miao. Moving loads identification of bridge based on influence line[J]. Journal of Vibration and Shock， 2013， 32（3）：129—133.

[13] 王 蕾， 侯吉林， 歐進(jìn)萍. 基于荷載形函數(shù)的大跨橋梁結(jié)構(gòu)移動(dòng)荷載識(shí)別[J]. 計(jì)算力學(xué)學(xué)報(bào)， 2012， 29（02）：153—158.

WANG Lei， HOU Ji-lin， OU Jin-ping. Moving force identification based on load shape function for a long-span bridge structure[J]. Chinese Journal of Computational Mechanics， 2012， 29（02）：153—158.

[14] 尤 瓊， 史治宇， 羅紹湘. 基于小波有限元的移動(dòng)荷載識(shí)別[J]. 振動(dòng)工程學(xué)報(bào)， 2010， 23（2）：188—193.

YOU Qiong， SHI Zhi-yu， LAW Siu Seong. Moving force identification based on wavelet finite element method[J]. Journal of Vibration Engineering， 2010， 23（2）：188—193.

[15] 修建權(quán)， 尤 瓊， 史治宇. 基于小波有限元法的連續(xù)梁移動(dòng)荷載識(shí)別[J]. 振動(dòng)與沖擊， 2012， 31（6）：141—144.

XIU Jian-quan， YOU Qiong， SHI Zhi-yu. Moving force identification on continuous bridges with wavelet finite element method[J]. Journal of Vibration and Shock， 2012， 31（6）：141—144.

[16] Yu Y， Cai C S， Deng L. State-of-the-art review on bridge weigh-in-motion technology[J]. Advances in Structural Engineering， 2016， 19（9）： 1—17.

[17] Qiao B， Zhang X， Wang C， et al. Sparse regularization for force identification using dictionaries[J]. Journal of Sound and Vibration， 2016， 368：71—86.

[18] Pan C D， Yu L， Liu H L， et al. Moving force identification based on redundant concatenated dictionary and weighted l1-norm regularization[J]. Mechanical Systems & Signal Processing， 2018， 98：32—49.

[19] Beck A， Teboulle M. A Fast Iterative Shrinkage Thresholding Algorithm for Linear Inverse Problems[J]. Siam Journal on Imaging Sciences， 2009， 2（1）：183—202.

Abstract： Moving force identification （MFI） from measured responses of bridges is a hot topic in bridge engineering. Current regularization-based MFI methods do not consider any influences from unknown initial conditions. In this study， a novel method is proposed for solving the MFI problem under unknown initial conditions based on sparse regularization. Firstly， the redundant dictionary and modal space are used for expressing moving forces and initial conditions， respectively. Responses induced by each component of moving forces and initial conditions are calculated and used for forming an image dictionary via normalization. The measured responses are decomposed by the image dictionary and then the decomposed results are used for estimating the moving forces. Finally， numerical simulations are used for assessing the accuracy and feasibility of the proposed method. A vehicle-bridge system model is then designed and fabricated in laboratory. A series of responses of a beam bridge subjected to moving vehicles are measured and used to verify the proposed method. The results show that the proposed method can identify the moving forces under unknown initial conditions.

Key words： moving force identification； bridge； initial conditions； sparse regularization； sparse decomposition

作者簡介： 潘楚東（1990—），男，博士研究生。電話：13427589181；E-mail：pan_cd@foxmail.com

通訊作者： 余 嶺（1963—），男，教授。電話：（020）85220476轉(zhuǎn)2012；E-mail：lyu1997@163.com