唐和生　王澤宇　陳嘉緣

摘要：土木工程實際結(jié)構(gòu)損傷狀態(tài)的時間跨度通常只占總生命周期的一小部分。為解決傳統(tǒng)基于數(shù)據(jù)驅(qū)動的結(jié)構(gòu)損傷識別方法缺乏足夠多的損傷訓(xùn)練數(shù)據(jù)的問題，提出結(jié)合數(shù)字孿生和深度學(xué)習(xí)的結(jié)構(gòu)損傷識別方法，并應(yīng)用于實際工程。該方法利用數(shù)值仿真模型和在線監(jiān)測數(shù)據(jù)構(gòu)建結(jié)構(gòu)的數(shù)字孿生，以獲得不同損傷工況下結(jié)構(gòu)動力響應(yīng)的“大數(shù)據(jù)”；為了擺脫對外激勵信息的依賴，應(yīng)用經(jīng)驗?zāi)B(tài)分解法和傳遞率函數(shù)對得到的數(shù)據(jù)進行預(yù)處理；將預(yù)處理后的固有模態(tài)傳遞率函數(shù)數(shù)據(jù)作為深度學(xué)習(xí)的輸入進行訓(xùn)練，實現(xiàn)結(jié)構(gòu)的損傷識別。為驗證方法的有效性，對實際結(jié)構(gòu)未經(jīng)訓(xùn)練的監(jiān)測數(shù)據(jù)進行分析，結(jié)果表明，該方法泛化能力良好，能夠有效識別結(jié)構(gòu)損傷狀況。通過數(shù)字孿生技術(shù)解決了傳統(tǒng)方法數(shù)據(jù)匱乏的問題，不需要任何地震信息，利用固有模態(tài)傳遞率函數(shù)數(shù)據(jù)訓(xùn)練的深度神經(jīng)網(wǎng)絡(luò)仍能保持較高的損傷識別準(zhǔn)確率，二者結(jié)合可以使工程結(jié)構(gòu)健康監(jiān)測更為主動、可靠、高效。

關(guān)鍵詞：數(shù)字孿生；深度學(xué)習(xí)；固有模態(tài)傳遞率函數(shù)；損傷識別；結(jié)構(gòu)健康監(jiān)測

中圖分類號：TU317；TP183? ? ?文獻標(biāo)志碼：A? ? ?文章編號：2096-6717（2024）01-0110-12

Structural damage identification based on digital twin and deep learning

TANG Hesheng， WANG Zeyu， CHEN Jiayuan

（College of Civil Engineering， Tongji University， Shanghai 200092， P. R. China）

Abstract： The time span of the civil engineering structural damage state usually accounts for a small part of the total life cycle. In order to solve the problem that traditional data-driven structural damage identification methods lack enough damage state data for training， a structural damage identification method based on digital twins and deep learning is proposed in this paper for practical application. Firstly， the digital twin is constructed by using the numerical simulation model and online monitoring data to obtain the “big data” of the structural dynamic response under different damage conditions. Secondly， to get rid of the dependence on the external excitation， the empirical mode decomposition method and transmissibility function are used to preprocess the obtained data. Then， the damage identification is realized by using deep learning. To verify the effectiveness of this method， untrained monitoring data of structures are analyzed. The results show that the method has good generalization ability and can identify the structural damage condition effectively. The problem of data hunger is solved by digital twin technology， and the deep neural network trained by the intrinsic mode vibration transmissibility function data sets can still maintain a high accuracy of damage identification without any seismic information. The combination of the two methods can make structural health monitoring more active， reliable and efficient.

Keywords： digital twin； deep learning； intrinsic mode vibration transmissibility function； damage identification； structural health monitoring

結(jié)構(gòu)物在建成使用后，由于長期暴露于操作荷載、環(huán)境影響及各種意外事件中，不可避免地會受到一些破壞。大型工程結(jié)構(gòu)的性能退化作為大部分突發(fā)性災(zāi)難事故的導(dǎo)火索，使得采用科學(xué)的方法對大型工程結(jié)構(gòu)進行健康監(jiān)測的重要性不斷提升[1]。結(jié)構(gòu)損傷識別是結(jié)構(gòu)健康監(jiān)測系統(tǒng)的核心，是健康監(jiān)測領(lǐng)域具有挑戰(zhàn)性的課題[2]。隨著大數(shù)據(jù)技術(shù)的發(fā)展，基于數(shù)據(jù)驅(qū)動的結(jié)構(gòu)損傷識別方法逐漸成為研究熱點。在基于數(shù)據(jù)的結(jié)構(gòu)損傷識別領(lǐng)域中，機器學(xué)習(xí)由于其卓越的數(shù)據(jù)處理能力，為許多研究者所青睞[3]。如邵會辰[4]選取平面單元模態(tài)應(yīng)變能變化率作為識別指標(biāo)，應(yīng)用徑向基神經(jīng)網(wǎng)絡(luò)和廣義回歸神經(jīng)網(wǎng)絡(luò)識別4層框架模型的損傷狀態(tài)，識別誤差率在5%左右，驗證了該方法的有效性。駱勇鵬等[5]提出采用單傳感器數(shù)據(jù)結(jié)合格拉姆角場和卷積神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)損傷識別方法，對一榀鋼框架結(jié)構(gòu)完成損傷識別任務(wù)，在測試樣本上的準(zhǔn)確率達100%。Dang等[6]通過將多個加速度傳感器采集到的數(shù)據(jù)作為輸入，分別訓(xùn)練多層感知器、長短記憶神經(jīng)網(wǎng)絡(luò)和卷積神經(jīng)網(wǎng)絡(luò)，實現(xiàn)了地震動下二維鋼框架的損傷識別，其中，卷積神經(jīng)網(wǎng)絡(luò)方法相較于前兩者展示了極高的可靠性。

然而，缺乏足夠龐大的標(biāo)簽數(shù)據(jù)來訓(xùn)練和測試模型成了限制基于數(shù)據(jù)驅(qū)動機器學(xué)習(xí)的損傷識別方法應(yīng)用于實際工程的一個關(guān)鍵因素。雖然得益于傳感器的發(fā)展，來自結(jié)構(gòu)健康監(jiān)測系統(tǒng)的數(shù)據(jù)量正在增長，但也只是數(shù)據(jù)規(guī)模的一小部分。在結(jié)構(gòu)的全生命周期中，有85%～90%的時間是處正常工作運營狀態(tài)，特別缺少極端條件下出現(xiàn)損傷模式的信息特征[7]。此外，通過人工直接從實體模型中獲取數(shù)據(jù)的成本高昂、耗時，且獲取的數(shù)據(jù)集質(zhì)量直接影響著識別準(zhǔn)確率。隨著第4次工業(yè)革命的開展，人工智能、物聯(lián)網(wǎng)等技術(shù)快速發(fā)展[8]，數(shù)字孿生成為實現(xiàn)智慧建造的重要途徑，為解決上述問題提供了新的范式，為工程系統(tǒng)開發(fā)數(shù)字孿生以用于健康監(jiān)測的興趣正在迅速增長[9]。Karve等[10]開發(fā)了一個包括損傷診斷、損傷預(yù)測和任務(wù)優(yōu)化的數(shù)字孿生智能規(guī)劃方法，利用數(shù)字孿生進行GP模型訓(xùn)練后誤差為1.19%，證明了該方法可以成功用于執(zhí)行任務(wù)優(yōu)化，在保證安全的前提下實現(xiàn)預(yù)期系統(tǒng)性能的目標(biāo)。Kapteyn等[11]利用實際飛機部件和傳感器數(shù)據(jù)構(gòu)建了該飛機結(jié)構(gòu)的數(shù)字孿生，說明了基于數(shù)據(jù)驅(qū)動的數(shù)字孿生使得飛機能夠動態(tài)地重新規(guī)劃安全任務(wù)，以監(jiān)測結(jié)構(gòu)損傷和退化情況。Gardner等[12]通過本構(gòu)關(guān)系構(gòu)建3層剪切鋼框架的初始驗證模型作為數(shù)字孿生數(shù)值模型，來預(yù)測結(jié)構(gòu)每層的加速度響應(yīng)，構(gòu)建好的數(shù)字孿生模型各層加速度響應(yīng)與實際加速度響應(yīng)的均方誤差為{0.260， 2.428， 2.939}，實現(xiàn)了對加速度響應(yīng)的預(yù)測。

數(shù)字孿生的概念模型最早出現(xiàn)于2003年，由Gieves教授在美國密歇根大學(xué)的產(chǎn)品全生命周期管理課程上提出，并在2010年由美國國家航空航天局（National Aeronautics and Space Administration， NASA）在太空技術(shù)路線圖中首次引入。自此以后，數(shù)字孿生便迸發(fā)出生機與活力[13]。雖然數(shù)字孿生的定義沒有被明確下來，但研究者們一致認為，數(shù)字孿生的優(yōu)勢在于它作為物理實體的虛擬表示，能夠反映物理實體的服役情況，使得工程師能夠遠程監(jiān)控其運行狀態(tài)，提供了實時反饋和預(yù)測潛在故障的能力。不難看出，數(shù)字孿生與傳統(tǒng)基于物理實體構(gòu)建的有限元模型不同，后者在建模過程中會受到相當(dāng)大不確定因素的影響，例如，參數(shù)取值不確定、物理過程不明確等，即傳統(tǒng)有限元模型難以擬合物理實體全生命周期的表現(xiàn)，導(dǎo)致無法使用基于物理模型輸出的數(shù)據(jù)來完成結(jié)構(gòu)健康監(jiān)測任務(wù)。而數(shù)字孿生則通過從物理角度收集的實測數(shù)據(jù)來連接孿生體和物理部分，從而對基于物理模型進行校準(zhǔn)，提高模型對物理實體的表示能力。以此構(gòu)建的數(shù)字孿生模型在確保了物理可解釋性的同時，能夠輸出結(jié)構(gòu)在各種工況下的大量可靠數(shù)據(jù)，解決了基于機器學(xué)習(xí)的損傷識別方法數(shù)據(jù)匱乏的問題。目前，許多研究都矚目于將數(shù)字孿生技術(shù)與機器學(xué)習(xí)結(jié)合，Xu等[14]提出了基于深度遷移學(xué)習(xí)和數(shù)字孿生的故障診斷方法，以某汽車車身側(cè)部生產(chǎn)線為例，驗證了該方法相較于傳統(tǒng)的損傷識別方法能夠獲得更大的數(shù)據(jù)集和更可靠的識別準(zhǔn)確率。Ritto等[9]通過建立6個自由度的離散數(shù)值模型來構(gòu)建懸臂梁的數(shù)字孿生模型，并將其與支持向量機、K近鄰、二次判別器等結(jié)合應(yīng)用于懸臂梁健康監(jiān)測中，以選擇結(jié)合數(shù)字孿生后識別準(zhǔn)確率更高的機器學(xué)習(xí)算法。張勝文等[15]構(gòu)建了全生命周期的離心泵機數(shù)字孿生映射模型，應(yīng)用深度雙向長短時記憶網(wǎng)絡(luò)完成故障診斷過程的數(shù)字化、自動化和智能化。將數(shù)字孿生和機器學(xué)習(xí)共同賦能于結(jié)構(gòu)健康監(jiān)測，無疑可以推動工程領(lǐng)域向更高效的方向發(fā)展。

筆者以一棟裝有長期健康監(jiān)測系統(tǒng)的8層油阻尼器鋼框架結(jié)構(gòu)為研究對象（物理孿生體），結(jié)合數(shù)字孿生概念和深度學(xué)習(xí)技術(shù)，建立對大型土木工程結(jié)構(gòu)損傷識別的新方法。首先，根據(jù)結(jié)構(gòu)設(shè)計參數(shù)建立有限元模型，然后利用在線監(jiān)測數(shù)據(jù)進行結(jié)構(gòu)模型參數(shù)更新，構(gòu)建該結(jié)構(gòu)的數(shù)字孿生。依據(jù)該數(shù)字孿生模型構(gòu)建結(jié)構(gòu)不同損傷工況下的“大”數(shù)據(jù)?？紤]到實際工程應(yīng)用，避免依賴外激勵信息，將對此數(shù)字孿生數(shù)據(jù)和部分監(jiān)測數(shù)據(jù)進行預(yù)處理后再來訓(xùn)練深度學(xué)習(xí)，實現(xiàn)損傷分類。最后，利用結(jié)構(gòu)真實損傷情況下的數(shù)據(jù)進行驗證。

1 數(shù)字孿生框架

數(shù)字孿生理論發(fā)展至今，共整合得出了13個主要特征[16]。結(jié)合實際工程情況，將其中最為重要的5個特征列在表1中。

在這5個特征中，“狀態(tài)”彰顯了數(shù)字孿生與物理孿生之間的緊密聯(lián)系：數(shù)字孿生相較于傳統(tǒng)基于物理構(gòu)建的計算模型，能夠使用從物理角度收集的數(shù)據(jù)來校準(zhǔn)自身，即數(shù)字孿生為模型和數(shù)據(jù)相融合所形成的虛擬副本，從而度量計算模型由于參數(shù)取值的不確定、物理過程不明確等不確定因素，以提高其對物理孿生的表示能力，擬合物理孿生在物理環(huán)境下的表現(xiàn)。因此，有限元建模、模型參數(shù)更新、貝葉斯統(tǒng)計推斷等技術(shù)都可以應(yīng)用在數(shù)字孿生理念內(nèi)，并將彼此進行協(xié)調(diào)共同構(gòu)建數(shù)字孿生，以提供支持特定物理孿生目標(biāo)決策、提高預(yù)測能力。

如上所述，建立數(shù)字孿生需要一個發(fā)揮功用的計算模型并使用物理孿生測量的數(shù)據(jù)對其進行校準(zhǔn)。而后，可以利用構(gòu)建好的數(shù)字孿生進行諸多操作，這些操作可能由于成本高昂等因素難以在物理環(huán)境中完成。對于結(jié)構(gòu)損傷識別的目標(biāo)，則通過在計算機上人為破壞數(shù)字孿生的構(gòu)件模型，生成多種結(jié)構(gòu)損傷工況，以提供大量結(jié)構(gòu)損傷狀況下的可靠數(shù)據(jù)，補充健康監(jiān)測數(shù)據(jù)庫中損傷信息的缺乏，解決了傳統(tǒng)基于數(shù)據(jù)驅(qū)動的損傷識別方法受限于數(shù)據(jù)集的大小和質(zhì)量的問題。對于數(shù)據(jù)庫中“大”數(shù)據(jù)的處理，則非常適合將其與機器學(xué)習(xí)相結(jié)合，利用機器學(xué)習(xí)優(yōu)秀的特征處理能力，為結(jié)構(gòu)損傷識別賦能，補充健康監(jiān)測系統(tǒng)對結(jié)構(gòu)的監(jiān)控能力，針對性地指出損傷的發(fā)生，指導(dǎo)研究人員對兩個孿生體針對性地做出更新。圖1展示了基于數(shù)字孿生和機器學(xué)習(xí)的結(jié)構(gòu)損傷識別框架，以最大程度地發(fā)揮數(shù)字孿生和機器學(xué)習(xí)的效用。

2 數(shù)字孿生的構(gòu)建

2.1 物理孿生

以一棟油阻尼器鋼框架建筑結(jié)構(gòu)為研究對象（圖2）。該結(jié)構(gòu)東西向長48 m，南北向?qū)?.6 m，總高度34.2 m。其中，底層為挑空層，包含了層高為4 m的1層和2層，形成了結(jié)構(gòu)高度8 m的底部大空間，在后文中稱為第1層，其余各層的結(jié)構(gòu)高度均為3.8 m。為提高抗震性能，在結(jié)構(gòu)在長短邊方向均布置了阻尼器，共計56組，安裝的油阻尼器包括兩種型號，第1層采用同一型號的阻尼器，后文稱ModelⅠ，3至8層采用另一型號阻尼器，后文稱ModelⅡ。健康監(jiān)測系統(tǒng)在第1、4、8層布置了雙向加速度計（共3個），并于2016年9月完成了監(jiān)測系統(tǒng)外部裝置監(jiān)視器的安裝。圖3、圖4分別為加速度計與油阻尼器在結(jié)構(gòu)中的布置圖和油阻尼器實景圖，圖5為該建筑健康監(jiān)測系統(tǒng)的監(jiān)測界面。

在2011年3月11日日本大地震的作用下（后文簡稱311大地震），位于結(jié)構(gòu)1層的8組油阻尼器被完全破壞；位于3層和4層的油阻尼器雖然沒有完全破壞，但油液已發(fā)生了泄漏，無法再提供阻尼恢復(fù)力；5～8層的阻尼器仍處于正常工作狀態(tài)，保持良好。由于地震當(dāng)時全樓的電力系統(tǒng)失效，計算機沒有記錄下該結(jié)構(gòu)的地震數(shù)據(jù)。圖6給出了一組采集于2011年3月9日（311大地震前震）的結(jié)構(gòu)東西向?qū)崪y的加速度響應(yīng)數(shù)據(jù)及其傅里葉幅值譜。在此次地震下，結(jié)構(gòu)屬于完好狀態(tài)，使用該組數(shù)據(jù)來更新計算模型，用以構(gòu)建結(jié)構(gòu)的初始狀態(tài)數(shù)字孿生。

2.2 有限元模型構(gòu)建

首先構(gòu)建油阻尼器鋼框架結(jié)構(gòu)的計算模型，而后利用311大地震前震的實測數(shù)據(jù)進行校準(zhǔn)。圖7為利用SAP2000軟件創(chuàng)建的結(jié)構(gòu)初始有限元模型，梁柱定義為線單元中的框架單元，樓板定義為面單元中的殼單元，阻尼器定義為連接單元，各構(gòu)件的尺寸與材料性能參數(shù)均按照標(biāo)準(zhǔn)設(shè)計參數(shù)定義。

由于量測的結(jié)構(gòu)加速度響應(yīng)數(shù)據(jù)僅與平動模態(tài)相關(guān)，且初始有限元模型模態(tài)分析表明結(jié)構(gòu)的平動和轉(zhuǎn)動模態(tài)相互分離，故可將初始有限元模型簡化處理為平面模型[17]，在保證計算模型有效性的同時便于后續(xù)的更新操作。忽略構(gòu)件的豎向變形，考慮模型側(cè)向變形及節(jié)點轉(zhuǎn)動，并假定樓板為剛性[18]，以此簡化的有限元模型僅包含7個自由度，與結(jié)構(gòu)樓層相對應(yīng)，簡化模型的質(zhì)量矩陣（單位：t）和剛度矩陣（單位：MN/mm）分別為

結(jié)構(gòu)阻尼考慮經(jīng)典的阻尼模型，阻尼比為0.02。阻尼器計算模型基于文獻[17]的識別結(jié)果，采用Maxwell模型，其參數(shù)如表2所示。

采用Newmark-β法預(yù)測簡化有限元模型的各層加速度響應(yīng)，圖8為簡化有限元模型預(yù)測與實際監(jiān)測數(shù)據(jù)（311大地震前震）在強震階段（30～60 s）的加速度時域及頻域的對比?？梢园l(fā)現(xiàn)，結(jié)構(gòu)簡化有限元模型不能準(zhǔn)確模擬結(jié)構(gòu)真實的動力行為，存在較大的誤差，無法作為可操作的數(shù)字孿生體為后續(xù)損傷識別提供“大”數(shù)據(jù)，需要對簡化有限元模型進行修正以獲取更精確的響應(yīng)預(yù)測。

2.3 不確定性度量

將不確定性考量為結(jié)構(gòu)參數(shù)的不確定性?？紤]到簡化模型具有較少的自由度，采用模型修正方法識別結(jié)構(gòu)模態(tài)參數(shù)[18]，使得修正后的簡化有限元模型響應(yīng)可以很好擬合實際結(jié)構(gòu)響應(yīng)，完成對計算模型的不確定性分析。

圖9給出了在強震階段模型修正后的預(yù)測和實測加速度響應(yīng)對比及頻域?qū)Ρ?。?給出了實測、初始模型和Berman-Baruch法修正模型的前兩階模態(tài)參數(shù)?？梢园l(fā)現(xiàn)，經(jīng)修正后的模型預(yù)測響應(yīng)與實際監(jiān)測數(shù)據(jù)的擬合程度很好。由此，針對實體結(jié)構(gòu)的原始狀態(tài)數(shù)字孿生構(gòu)建完成。

對實際監(jiān)測和有限元模型預(yù)測數(shù)據(jù)的對比說明僅依靠有限元模型不能滿足數(shù)字孿生的需要。在輔以模型修正方法后，修正后的有限元模型可以很好地擬合物理孿生在物理環(huán)境下的表現(xiàn)，由此利用構(gòu)建好的數(shù)字孿生擬合結(jié)構(gòu)不同損傷工況，可以提供不同損傷狀態(tài)的“大”數(shù)據(jù)集，補充健康監(jiān)測數(shù)據(jù)庫中損傷信息的缺乏，為機器學(xué)習(xí)提供“大”數(shù)據(jù)支撐。

3 基于數(shù)字孿生和深度學(xué)習(xí)的損傷識別

3.1 損傷工況的定義

根據(jù)311大地震對實際結(jié)構(gòu)的破壞情況，基于以下假設(shè)，將實際結(jié)構(gòu)損傷情況分為如表4所示的4類。注：由于篇幅限制，僅對311大地震后結(jié)構(gòu)實際發(fā)生的損傷工況進行定義。在實際應(yīng)用中，可以對數(shù)字孿生模型進行針對性處理，以模擬可能出現(xiàn)的多種損傷情況。

假設(shè)1：由311大地震后對結(jié)構(gòu)實際的檢修結(jié)果表明，實際結(jié)構(gòu)僅發(fā)生首3層油阻尼器油液泄露的情況，未產(chǎn)生剛度損傷，因此，僅考慮阻尼器損壞的工況。

假設(shè)2：對于安裝在同一層的阻尼器，假定所有阻尼器工作時的位移、速度均相同，且所有阻尼器本身的屬性也完全一致。因此，如果某層的單個阻尼器已經(jīng)破壞，則其余的阻尼器也均考慮遭受了破壞。

假設(shè)3：對于一般的結(jié)構(gòu)，阻尼器安裝的層數(shù)越低，其發(fā)揮的效用也越大，并考慮安裝在高層的阻尼器發(fā)揮的阻滯作用有限。

假設(shè)4：基于假設(shè)3，考慮除了1、2、3層的阻尼器被完全破壞，4層及以上的阻尼器工作狀態(tài)安全無虞。

3.2 數(shù)據(jù)預(yù)處理

結(jié)構(gòu)動力響應(yīng)時頻信息包含了豐富的結(jié)構(gòu)狀態(tài)特征[19]，同時考慮到實際工程應(yīng)用，將采用經(jīng)驗?zāi)B(tài)分解時頻信息和振動傳遞率函數(shù)，構(gòu)建不依賴于外激勵的深度學(xué)習(xí)訓(xùn)練數(shù)據(jù)集。

3.2.1 經(jīng)驗?zāi)B(tài)分解

結(jié)構(gòu)的振動信號千變?nèi)f化，其中，大部分信號為非平穩(wěn)、非高斯分布和非線性的隨機信號，如地震動信號等。并且當(dāng)工程結(jié)構(gòu)出現(xiàn)損傷時上述情況顯得更為突出。而基于傅里葉變換的譜分析方法主要適用于平穩(wěn)信號（時不變信號），對非平穩(wěn)信號進行傅里葉變換只能分別給出時域或頻域的統(tǒng)計平均結(jié)果，無法滿足結(jié)構(gòu)損傷識別的要求。經(jīng)驗?zāi)B(tài)分解（Empirical Mode Decomposition，EMD）作為一種自適應(yīng)的信號處理方法，能把復(fù)雜的信號分解為有限的固有模態(tài)函數(shù)（Intrinsic Mode Function，IMF）之和，非常適合處理非線性和非平穩(wěn)信號[20]。信號x（t）經(jīng)EMD分解后可以表示為

式中：c_i為第i個固有模態(tài)函數(shù)；r_n為殘余函數(shù)。c_1，c_2…c_n分別包含了信號從高到低不同頻率段的成分，每個IMF分量（c_i）都反映了分析信號在一個時間尺度上內(nèi)在的模態(tài)特性。綜上，可以首先對數(shù)字孿生輸出的加速度信號進行EMD分解，而后基于文獻[19]和文獻[20]的研究，選取對損傷程度敏感且更能體現(xiàn)原信號趨勢走向的高階IMF分量進行后續(xù)分析。

3.2.2 傳遞率函數(shù)

傳遞率函數(shù)的概念來源于自動控制理論，其實質(zhì)上是用零初始條件下的輸入信號的拉普拉斯變換和輸出信號的額拉普拉斯變換作比，對系統(tǒng)本身的固有屬性進行描述。傳遞率函數(shù)囊括了系統(tǒng)的諸多信息，又與輸入信號本身無關(guān)，因此，對于系統(tǒng)本身的變化具有很好的描述性。

對于一般的n自由度結(jié)構(gòu)，其運動微分方程可以表達為

根據(jù)式（9）可以發(fā)現(xiàn)，振動傳遞率函數(shù)是頻響函數(shù)的函數(shù)，在反映結(jié)構(gòu)動力學(xué)特性的同時，地震激勵只作為動力源而沒有參與運算，從而避免了在未來應(yīng)用中對其進行量測。研究關(guān)注的是監(jiān)測樓層之間的傳遞率函數(shù)，因此，可以通過監(jiān)測樓層之間的頻域相應(yīng)比值來獲取結(jié)構(gòu)響應(yīng)傳遞率函數(shù)。

3.2.3 IMF與傳遞率函數(shù)的結(jié)合

基于經(jīng)驗?zāi)B(tài)分解和傳遞率函數(shù)在結(jié)構(gòu)響應(yīng)信號處理時的優(yōu)點[19]，采用將固有模態(tài)函數(shù)和傳遞率函數(shù)相結(jié)合的方法來構(gòu)建數(shù)據(jù)集。首先求解結(jié)構(gòu)監(jiān)測樓層加速度響應(yīng)信號經(jīng)EMD分解后的固有模態(tài)函數(shù)IMF，然后取每層的前3個IMF分量進行時頻變換，最后按IMF分量的順序，依次求得在該IMF分量下結(jié)構(gòu)各相鄰層的IMF傳遞率函數(shù)。具體求解流程如下。

1）在地震波的作用下，結(jié)構(gòu)監(jiān)測樓層（1、4、8層）的加速度時域響應(yīng)為

至此，完成數(shù)字孿生輸出的結(jié)構(gòu)監(jiān)測樓層加速度響應(yīng)數(shù)據(jù)的預(yù)處理。

3.3 訓(xùn)練數(shù)據(jù)集構(gòu)建

為了在天然地震波的基礎(chǔ)上得到更多組訓(xùn)練數(shù)據(jù)來充分訓(xùn)練后文的機器學(xué)習(xí)分類器，采用對天然地震波人工加噪的數(shù)據(jù)增強方法，以此法共生成的42條地震波作為不同的地震激勵進行輸入，如表5所示。對于表5中每一條地震波輸入，均模擬表4中的4種工況，計算每種情況下的IMF傳遞率函數(shù)，最終可以得到168組IMF傳遞率函數(shù)矩陣?？紤]到地震波的頻率分量較低，造成的結(jié)構(gòu)響應(yīng)也多為低頻振動，因此，僅保留IMF傳遞率函數(shù)的0～10 Hz段。圖10是在0309AccEW激勵下1、2、3層阻尼器均被完全破壞情況下的傳遞率函數(shù)矩陣。

3.4 基于卷積神經(jīng)網(wǎng)絡(luò)的損傷識別

機器學(xué)習(xí)發(fā)展至今，深度學(xué)習(xí)展現(xiàn)出了優(yōu)異的對樣本數(shù)據(jù)內(nèi)在規(guī)律的學(xué)習(xí)性能，與支持向量機、決策樹等傳統(tǒng)的機器學(xué)習(xí)方法相比，在語音、圖像識別等方面的效果得到了很大的提升。在深度學(xué)習(xí)中，卷積神經(jīng)網(wǎng)絡(luò)（Convolutional Neural Network，CNN）更是成為研究的熱點，最常用于分析視覺圖像。選用AlexNet網(wǎng)絡(luò)進行IMF傳遞率函數(shù)矩陣圖的分類，完成損傷識別的目的。

AlexNet網(wǎng)絡(luò)由8層組成，包括5層卷積層和3層全連接層（圖11）。首先，在輸入層將IMF傳遞率函數(shù)矩陣圖歸一化至227×227×3像素大小，然后在第1個卷積層使用96個卷積核提取邊緣特征，并將卷積層得到的特征圖放入激活函數(shù)層，激活函數(shù)選用ReLU函數(shù)進行非線性映射。然后將ReLU層的輸出輸入進池化層，池化層通過最大池化算法提取顯著特征。這些步驟在以下4個具有不同內(nèi)核和填充大小的卷積層中重復(fù)執(zhí)行。隨著卷積層數(shù)的增加，提取的特征數(shù)量也顯著增加。之后，將特征圖連接到全連接層。Dropout層以50%的概率將神經(jīng)元從網(wǎng)絡(luò)中移除，從而避免過擬合的發(fā)生。在輸出層使用Softmax函數(shù)對每個分類的可能性進行歸一化，得到分類結(jié)果。

3.5 訓(xùn)練過程與結(jié)果

將168組IMF傳遞率函數(shù)矩陣數(shù)據(jù)集分為4類混合輸入到AlexNet網(wǎng)絡(luò)中，設(shè)置學(xué)習(xí)速率為0.001，共設(shè)置120個epoch。在每個epoch中，圖像被隨機分成訓(xùn)練集和驗證集，其中訓(xùn)練集占總圖像的70%，驗證集占總圖像的30%（圖12）。

隨著迭代次數(shù)的增加，預(yù)測精度也在逐步提高，當(dāng)?shù)螖?shù)超過300次時，預(yù)測精度接近恒定值。圖13展示了這一訓(xùn)練過程。在驗證集的表現(xiàn)上，圖14展示了在驗證集上分類的混淆矩陣。混淆矩陣中的每一列代表了分類器預(yù)測樣本的類別，每一行代表了樣本的真實歸屬類別。其中，綠色色塊表示真實樣本被分類器正確預(yù)測，紅色色塊表示真實樣本被分類器錯誤預(yù)測?；煜仃囉蚁陆菫樽R別準(zhǔn)確率，表示所有預(yù)測正確的樣本占所有樣本的比例?；煜仃囎畹撞康男斜硎揪_率，最右側(cè)的列表示召回率。圖14展示了最終識別準(zhǔn)確率達92.3%，其中，Ⅰ類的精確率和召回率最高，其次為Ⅱ、Ⅳ類，Ⅲ類的精確率和召回率均較低。造成上述結(jié)果的可能原因為：首先，I類為未損傷工況，Ⅱ、Ⅲ、Ⅳ類均為損傷工況，因此，分類器對于I類的分類效果最好；其次，Ⅲ類損傷程度介于Ⅱ類和Ⅳ類之間，會出現(xiàn)將Ⅲ類誤分類為Ⅱ類和Ⅳ類，或?qū)ⅱ蝾惡廷纛愓`分類為Ⅲ類兩種情況。又因為對于一般結(jié)構(gòu)，阻尼器安裝的層數(shù)越低，其發(fā)揮的效用也越大，且結(jié)構(gòu)第1層為8 m大空間，其余各層均為3.8 m。綜上，分類器將Ⅲ類誤分類為Ⅱ類和Ⅳ類，造成了Ⅲ類的精確率和召回率均較低的分類結(jié)果。

3.6 實際工程結(jié)構(gòu)驗證

為評價分類器的泛化能力，驗證方法的有效性，將未訓(xùn)練過的實際結(jié)構(gòu)響應(yīng)數(shù)據(jù)經(jīng)數(shù)據(jù)預(yù)處理后構(gòu)建測試集，輸入至已訓(xùn)練好的深度學(xué)習(xí)分類器中。用于構(gòu)建測試集的實際結(jié)構(gòu)響應(yīng)數(shù)據(jù)描述如表6所示。由于篇幅限制，僅展示在地震波record20100613作用下，結(jié)構(gòu)第1、4、8層加速度響應(yīng)數(shù)據(jù)（圖15）。

圖16為根據(jù)表6構(gòu)建的測試集的訓(xùn)練結(jié)果混淆矩陣，在測試集中，深度學(xué)習(xí)分類器準(zhǔn)確識別出了實際結(jié)構(gòu)的損傷情況，測試集上識別準(zhǔn)確率為100%。展現(xiàn)該深度學(xué)習(xí)分類器良好的泛化能力，驗證了該方法在實際結(jié)構(gòu)進行實時損傷識別的有效性。

結(jié)合實際結(jié)構(gòu)發(fā)生的損傷對構(gòu)建的數(shù)字孿生定義4種損傷工況，并對4種損傷工況下輸出的數(shù)據(jù)集進行基于IMF振動傳遞率函數(shù)的數(shù)據(jù)預(yù)處理，處理好的數(shù)據(jù)集輸入到卷積神經(jīng)網(wǎng)絡(luò)中進行訓(xùn)練，最終在驗證集上識別準(zhǔn)確率達92.3%。通過對未訓(xùn)練過的數(shù)據(jù)構(gòu)建測試集，驗證了基于IMF振動傳遞率函數(shù)的數(shù)據(jù)預(yù)處理方法可以擺脫對激勵信息的依賴，深度學(xué)習(xí)分類器泛化能力良好，將二者結(jié)合可以滿足對實際工程損傷識別的需要。

4 結(jié)論

數(shù)字孿生作為健康監(jiān)測、故障診斷的新范式，使得工程結(jié)構(gòu)的監(jiān)測和管理更為主動、可靠、高效。

提出了將數(shù)字孿生和深度學(xué)習(xí)相結(jié)合的結(jié)構(gòu)損傷識別方法，并將其應(yīng)用到土木工程大型結(jié)構(gòu)中。主要結(jié)論如下：

1）數(shù)字孿生技術(shù)可以充分利用物理模型、健康監(jiān)測系統(tǒng)傳感器數(shù)據(jù)更新、運行歷史等數(shù)據(jù)，在虛擬空間中實現(xiàn)物理實體結(jié)構(gòu)的映射。利用此技術(shù)，可以通過仿真獲取物理實體結(jié)構(gòu)不同服役狀態(tài)的“大”數(shù)據(jù)。該技術(shù)突破了傳統(tǒng)基于數(shù)據(jù)驅(qū)動的結(jié)構(gòu)損傷識別方法中數(shù)據(jù)缺乏的瓶頸，使數(shù)字孿生和深度學(xué)習(xí)的結(jié)合更廣泛地在實際工程結(jié)構(gòu)損傷識別中的應(yīng)用成為可能。

2）采用基于IMF的傳遞率函數(shù)構(gòu)建的深度學(xué)習(xí)訓(xùn)練集數(shù)據(jù)可以有效擺脫對激勵信息的依賴，適合實際工程應(yīng)用。以此數(shù)據(jù)訓(xùn)練好的卷積神經(jīng)網(wǎng)絡(luò)魯棒性良好，面對未知地震激勵時仍能給出準(zhǔn)確的識別結(jié)果。

參考文獻

[1]? 李宏男， 李東升. 土木工程結(jié)構(gòu)安全性評估、健康監(jiān)測及診斷述評[J]. 地震工程與工程振動， 2002， 22（3）： 82-90.

LI H N， LI D S. Safety assessment， health monitoring and damage diagnosis for structures in civil engineering [J]. Earthquake Engineering and Engineering Vibration， 2002， 22（3）： 82-90. （in Chinese）

[2]? 朱宏平， 余璟， 張俊兵. 結(jié)構(gòu)損傷動力檢測與健康監(jiān)測研究現(xiàn)狀與展望[J]. 工程力學(xué)， 2011， 28（2）： 1-11， 17.

ZHU H P， YU J， ZHANG J B. A summary review and advantages of vibration-based damage identification methods in structural health monitoring [J]. Engineering Mechanics， 2011， 28（2）： 1-11， 17. （in Chinese）

[3]? ZHAO R， YAN R Q， CHEN Z H， et al. Deep learning and its applications to machine health monitoring [J]. Mechanical Systems and Signal Processing， 2019， 115： 213-237.

[4]? 邵會辰. 基于改進均勻設(shè)計表的損傷識別方法研究及應(yīng)用[D]. 山東青島： 青島理工大學(xué)， 2020.

SHAO H C. Research and application of damage identification method based on improved uniform design table [D]. Qingdao， Shandong： Qingdao Tehcnology University， 2020. （in Chinese）

[5]? 駱勇鵬， 王林堃， 郭旭， 等. 利用單傳感器數(shù)據(jù)基于GAF-CNN的結(jié)構(gòu)損傷識別[J]. 振動 測試與診斷， 2022， 42（1）： 169-176， 202.

LUO Y P， WANG L K， GUO X， et al. Structural damage identification using single sensor data based on GAF-CNN [J]. Journal of Vibration， Measurement ＆ Diagnosis， 2022， 42（1）： 169-176， 202. （in Chinese）

[6]? DANG H V， RAZA M， NGUYEN T V， et al. Deep learning-based detection of structural damage using time-series data [J]. Structure and Infrastructure Engineering， 2021， 17（11）： 1474-1493.

[7]? BOOYSE W， WILKE D N， HEYNS S. Deep digital twins for detection， diagnostics and prognostics [J]. Mechanical Systems and Signal Processing， 2020， 140： 106612.

[8]? WAGG D J， WORDEN K， BARTHORPE R J， et al. Digital twins： state-of-the-art and future directions for modeling and simulation in engineering dynamics applications [J]. ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems Part B Mechanical Engineering， 2020， 6（3）： 030901.

[9]? RITTO T G， ROCHINHA F A. Digital twin， physics-based model， and machine learning applied to damage detection in structures [J]. Mechanical Systems and Signal Processing， 2021， 155： 107614.

[10]? KARVE P M， GUO Y， KAPUSUZOGLU B， et al. Digital twin approach for damage-tolerant mission planning under uncertainty [J]. Engineering Fracture Mechanics， 2020， 225： 106766.

[11]? KAPTEYN M G， KNEZEVIC D J， HUYNH D B P， et al. Data-driven physics-based digital twins via a library of component-based reduced-order models [J]. International Journal for Numerical Methods in Engineering， 2022， 123（13）： 2986-3003.

[12]? GARDNER P， BORGO MDAL， RUFFINI V， et al. Towards the development of an operational digital twin [J]. Vibration， 2020， 3（3）： 235-265.

[13]? GRIEVES M W. Product lifecycle management： The new paradigm for enterprises [J]. International Journal of Product Development， 2005， 2（1/2）： 71.

[14]? XU Y， SUN Y M， LIU X L， et al. A digital-twin-assisted fault diagnosis using deep transfer learning [J]. IEEE Access， 7： 19990-19999.

[15]? 張勝文， 楊凌翮， 程德俊. 數(shù)字孿生驅(qū)動的離心泵機組故障診斷方法研究[J]. 計算機集成制造系統(tǒng)， 2023， 29（5）： 1462-1470.

ZHANG S W， YANG L H， CHENG D J. Fault diagnosis method of centrifugal pump driven by digital twin [J]. Computer Integrated Manufacturing Systems， 2023， 29（5）： 1462-1470.（in Chinese）

[16]? JONES D， SNIDER C， NASSEHI A， et al. Characterising the digital twin： A systematic literature review [J]. CIRP Journal of Manufacturing Science and Technology， 2020， 29： 36-52.

[17]? 楊朋超， 薛松濤， 謝麗宇. 消能減震建筑結(jié)構(gòu)模態(tài)參數(shù)識別的貝葉斯方法[J]. 振動工程學(xué)報， 2021， 34（4）： 671-679.

YANG P C， XUE S T， XIE L Y. Bayesian modal parameters identification of passively controlled building structures [J]. Journal of Vibration Engineering， 2021， 34（4）： 671-679. （in Chinese）

[18]? 楊朋超， 薛松濤， 謝麗宇. 結(jié)構(gòu)動力模型的改進直接修正方法及工程應(yīng)用[J]. 建筑結(jié)構(gòu)學(xué)報， 2021， 42（3）： 34-40.

YANG P C， XUE S T， XIE L Y. An improved direct method for dynamic model updating and its practical engineering applications [J]. Journal of Building Structures， 2021， 42（3）： 34-40. （in Chinese）

[19]? 顧建祖， 郝文峰， 駱英， 等. 基于固有模態(tài)函數(shù)振動傳遞率的結(jié)構(gòu)損傷識別[J]. 建筑科學(xué)與工程學(xué)報， 2011， 28（1）： 27-32.

GU J Z， HAO W F， LUO Y， et al. Structural damage identification based on intrinsic mode function vibration transmissibility [J]. Journal of Architecture and Civil Engineering， 2011， 28（1）： 27-32. （in Chinese）

[20]? HUANG N， SHEN Z， LONG S， et al. The empirical model decomposition and the hilbert spectrum for nonlinear and non-stationary time series analysis [J]. Proceedings of the Royal Society A： Mathematical， Physical and Engineering Sciences， 1998， 454（1971）： 903-995.

（編輯? 胡玲）

DOI： 10.11835/j.issn.2096-6717.2022.130

收稿日期：2022?04?27

基金項目：上海市級科技重大專項（2021SHZDZX0100）；土木工程I類高峰學(xué)科建設(shè)經(jīng)費（2022-3-YB-07）

作者簡介：唐和生（1973- ），男，博士，研究員，博士生導(dǎo)師，主要從事AI科學(xué)計算交叉研究，E-mail：thstj@#edu.cn。

Received： 2022?04?27

Foundation items： Shanghai Municipal Science and Technology Major Project （No. 2021SHZDZX0100）; Top Discipline Plan of Shanghai Universities-Class I （No. 2022-3-YB-07）

Author brief： TANG Hesheng （1973- ）， PhD， researcher， doctorial supervisor， main research interest： AI scientific computing intersection， E-mail： thstj@#edu.cn.