趙國旗仇亞萍 駱英 馮侃

(江蘇大學土木工程與力學學院，江蘇鎮(zhèn)江212013)

生物、工程及交叉力學

基于細觀混凝土模型的時間逆轉(zhuǎn)損傷成像方法1)

趙國旗2)仇亞萍 駱英 馮侃

(江蘇大學土木工程與力學學院，江蘇鎮(zhèn)江212013)

提出了一種針對混凝土結(jié)構(gòu)損傷檢測的時間逆轉(zhuǎn)損傷成像方法.以檢測混凝土結(jié)構(gòu)中與骨料尺寸相近的微小損傷為目的，引入細觀混凝土隨機骨料模型，該模型將混凝土結(jié)構(gòu)視為由水泥漿基底、骨料及粘接層組成的三相復合材料，基于Monte Carlo隨機樣本原理并結(jié)合真實試件的骨料級配曲線建立.在數(shù)值模擬分析中，將生成含損傷的細觀模型導入有限元分析軟件進行超聲波場模擬，同時采用自適應(yīng)性強的時間逆轉(zhuǎn)模型(time reversedmodel,TRM)進行損傷定位.TRM分為正向檢測和逆時成像兩個部分：正向檢測過程得到包含損傷的一系列散射回波信號，從數(shù)值角度進行時間反演并作為逆時過程的輸入信號；逆時成像過程選用等效彈性參數(shù)模型，幾何尺寸與隨機骨料模型相同，時反信號在相應(yīng)幾何位置同時加載形成時反波場，時反波場在損傷位置會發(fā)生干涉疊加從而導致能量峰值的出現(xiàn),通過確定干涉峰值時刻,并獲取該時刻對應(yīng)原始波場以及小波變換能量場完成成像.與原始數(shù)據(jù)波場圖相比，小波變換處理成像結(jié)果消除了雜波干擾，成像結(jié)果更加清晰.進一步對等效彈性參數(shù)的取值進行討論，并且在骨料尺寸范圍內(nèi)調(diào)整損傷大小，結(jié)果顯示成像結(jié)果匹配度高，對于非均質(zhì)混凝土結(jié)構(gòu)的損傷檢測能很好滿足損傷定位需求.由此證明，時間逆轉(zhuǎn)成像方法對于具有復雜結(jié)構(gòu)的混凝土材料的損傷檢測具有較好的適用性.

混凝土結(jié)構(gòu)，細觀模型，時間逆轉(zhuǎn)法，超聲檢測，損傷成像

引言

當前國家經(jīng)濟發(fā)展快速，基礎(chǔ)建設(shè)規(guī)?？涨?，從樓房、橋梁到機場、海港，其數(shù)量與日俱增，而混凝土正是這些基礎(chǔ)設(shè)施的重要組成部分[1].混凝土結(jié)構(gòu)內(nèi)部在生產(chǎn)和使用過程中易產(chǎn)生損傷，損傷的積累將導致結(jié)構(gòu)失效，從而誘發(fā)災(zāi)難性事故.無損檢測技術(shù)成為目前預(yù)防災(zāi)難性事故的重要手段.

先前的損傷檢測研究將混凝土視為均質(zhì)材料，在有限元軟件模型庫中采用平均化參數(shù)處理方法；然而當結(jié)構(gòu)中損傷的尺寸與骨料相近時，平均化處理會引起檢測精度降低.基于細觀力學為背景的細觀混凝土模型已經(jīng)取得重大成果，細觀模型從細觀尺度出發(fā)將混凝土結(jié)構(gòu)視為由骨料、基底及粘接層組成的復合材料.在早期研究中以Wittmann等[2]的二維模型以及Schlangen等[3]的格構(gòu)模型模擬了混凝土不同組分的力學特性.近年來混凝土細觀分析的基礎(chǔ)是生成形狀、尺寸、位置分布都相似于真實混凝土粗骨料的隨機骨料結(jié)構(gòu)，文獻[4-5]分別對圓形骨料及橢圓形骨料的生成方法和投放技巧進行了報道；Wang等[6]較早提出隨機骨料模型并進行有限元網(wǎng)格的劃分；馬懷發(fā)等[7]依據(jù)Walaraven公式生成二維二級配和四級配混凝土隨機骨料模型；楊華等[8]對隨機骨料模型的混凝土彈性模量的確定進行了預(yù)測性研究.本文研究將細觀隨機骨料模型引入無損檢測以提高檢測精度，結(jié)合超聲檢測技術(shù)完成損傷定位成像.

近年來，超聲檢測技術(shù)由于具有能量高、穿透力強、檢測速度快以及安全方便等優(yōu)點，逐漸成為國內(nèi)外研究機構(gòu)關(guān)于混凝土無損檢測的研究重點，并得到廣泛應(yīng)用.超聲波用于混凝土檢測始于1949年，Leslie和Cheesman[9]首次采用超聲脈沖檢測混凝土結(jié)構(gòu)缺陷;然而由于實驗局限性，諸多影響因素在當時未得到驗證，因此未能用于實際工程檢測.21世紀以來，隨著電子技術(shù)的發(fā)展,超聲檢測技術(shù)進入數(shù)字時代，數(shù)值仿真模擬技術(shù)的突破促使超聲成像技術(shù)得到進一步發(fā)展，其中合成孔徑聚焦技術(shù)[10-11]、超聲相控陣技術(shù)[1213]、層析成像技術(shù)[14-15]等都取得了一定成果.時間逆轉(zhuǎn)技術(shù)由于其自適應(yīng)性特點逐漸由地球科學領(lǐng)域推廣到無損檢測領(lǐng)域(nondestructive testing,NDT).邱雷等[16-18]針對復合材料板結(jié)構(gòu)檢測中存在的信號信噪比低、頻散以及模式混疊問題提出了時間反轉(zhuǎn)聚焦成像法,有效提高了有用信號能量.Zhu等[1921]對Lamb的時逆聚焦技術(shù)提出了基于頻率--波數(shù)域內(nèi)的快速成像法.近年來,將時間逆轉(zhuǎn)技術(shù)應(yīng)用于混凝土檢測的嘗試逐漸引起學者關(guān)注.時間逆轉(zhuǎn)模型(time reversedmodeling,TRM)是探測地球物理學中一項波源定位技術(shù)[22].Saenger等[2326]將TRM技術(shù)應(yīng)用于多相高散射混凝土復合材料的檢測，TRM被用于對聲發(fā)射源進行定位和定性，與之前的檢測技術(shù)相比無需進行首時識別或預(yù)知檢測通道的相關(guān)信息.

為了彌補先前研究中將混凝土視為均質(zhì)材料而造成的仿真實驗實用性不足的缺陷,本文依據(jù)真實試件的骨料級配曲線建立更具有實際應(yīng)用價值的細觀混凝土模型，結(jié)合有限元軟件進行模擬，采用TRM技術(shù)完成損傷檢測.同時對比不同成像條件結(jié)果的精確性，提出更適用于實際混凝土結(jié)構(gòu)檢測的損傷成像方法.

1 時間逆轉(zhuǎn)成像法

時間逆轉(zhuǎn)成像法包含正向檢測和逆時成像兩個步驟：正向檢測過程完成損傷信號的采集，采用單點激勵信號得到正向波場及損傷散射波場，通過邊界傳感器接收散射回波信號；逆時成像過程是從數(shù)值角度將損傷散射信號在時域內(nèi)反轉(zhuǎn)并在相同位置加載，依據(jù)波前干涉疊加確定損傷位置.

1.1 正向檢測

確定檢測區(qū)域Ω，邊界表示為?Ω，以x∈Ω表示檢測區(qū)域中的位置.檢測時域確定為t∈[0,T]，T為終止時間.正向檢測過程用來獲取邊界一系列位移場時域信號

式中，x(k)代表傳感器位置，k為傳感器編號.

式中，N為邊界傳感器的個數(shù)[25].若加載點與損傷間的距離為L，傳感器與損傷間的距離為Lk，傳感器接收到損傷散射信號的時間為tk=(L+Lk)/vp，其中vp為縱波波速.

1.2 逆時成像

逆時過程采用等效彈性參數(shù)模型[2627]，模型參數(shù)由正向檢測過程實際測試得到，包括縱波波速vp,eff、橫波波速vs,eff以及密度ρeff.正向檢測過程中邊界?Ω記錄的位移信號為

式中,t∈[0,T].將位移信號在時域內(nèi)取反

時反信號將作為逆時成像過程的激勵信號.正向檢測過程記錄的k組信號所對應(yīng)的幾何位置將作為時反信號加載的位置，k組信號同時加載形成時反波場.信號中包含的損傷散射信號擁有同一個波源即損傷，該損傷又稱二次波源.時反處理后，傳感器加載的損傷散射信號到達損傷位置時間為

將tk=(L+Lk)/vp代入式(5)得tRk=T-L/vp，由此可見,時反過程中損傷散射信號回到損傷位置所用的時間與各個時反加載點幾何位置無關(guān)，僅與正向加載點與損傷的幾何位置以及波速有關(guān)，則k組損傷信號將同時回到波源,形成波前干涉,造成能量疊加.

1.3 成像條件

為了更好地呈現(xiàn)時反波場中的波前干涉疊加現(xiàn)象，引入成像條件進行干涉情況表征和損傷位置確定.全時域全場位移最大值成像[28]

最大位移值umax(x)為時反波場中每一點x∈Ω在t∈[0,T]內(nèi)的位移最大值.此外，另提取時反應(yīng)力場，依據(jù)全時域應(yīng)力最大值成像

式中，最大應(yīng)力值σmax(x)為時反波場中每一點x∈Ω在t∈[0,T]內(nèi)的應(yīng)力最大值.由于時反加載過程僅選取來自邊界數(shù)量有限位置上的信號，受雜波干擾會造成難以識別損傷信號形成的匯聚點，引入表征物理特性總和的能量場成像法[28]

式中，在全時域內(nèi)網(wǎng)格點處讀取應(yīng)力、應(yīng)變，文獻[29]證明了能量法對于混凝土中聲發(fā)射源的定位最為精確.

然而針對本文中混凝土損傷的主動檢測問題，由于正向檢測過程中提取的損傷散射信號比直達波、邊界反射波微弱，加上混凝土結(jié)構(gòu)對波形衰減嚴重[30]以及邊界臨近傳感器加載時會形成波形混疊等因素，全時域內(nèi)成像結(jié)果并不理想.若能獲取時反波場中損傷散射信號的干涉疊加產(chǎn)生峰值時刻對應(yīng)的波場，便能顯著提高對損傷位置的確定精度.因此，本文提出單時刻波場成像法，為了準確定位能量聚焦的位置，對波場中各檢測點采用小波變換法提取中心頻率對應(yīng)的能量幅值，并確定對應(yīng)的最大幅值出現(xiàn)時刻te,max，同時獲取該時刻時反波場以及小波變換能量場

成像定位誤差用下式進行表征

式中,?L表示定位位置與真實損傷之間的距離，Rmax為混凝土結(jié)構(gòu)中骨料的最大粒徑.

2 細觀混凝土模型的建立

從專注層面以及研究方法等側(cè)重點的不同，研究過程中將混凝土結(jié)構(gòu)分為3個尺度：宏觀尺度(macro-scale)、細觀尺度(meso-scale)以及微觀尺度(m icro-scale)[31].當檢測損傷尺寸與骨料尺寸相近時，選取細觀角度分析更為精確.

混凝土骨料的分布是隨機的，依據(jù)Monte Carlo[32]法進行隨機數(shù)的生成.生成的骨料采用經(jīng)典最大密實度曲線指導尺寸以及各尺寸對應(yīng)的數(shù)量分布

其中，D0為篩選直徑，p為經(jīng)過D0篩選的骨料所占的質(zhì)量百分比，Dmax為骨料最大粒徑.根據(jù)上式，Walraven推出二維截面內(nèi)D＜D0粒徑骨料出現(xiàn)的概率[33]計算公式

式中，pk為骨料體積占總體積百分比.據(jù)此可得出橫截面上骨料分布情況.本研究采用MATLAB軟件生成具有隨機骨料尺寸和隨機骨料分布的混凝土模型.

2.2 模型生成實例

依據(jù)實際混凝土試件,骨料粒徑為10mm～40mm,水泥富余系數(shù)為1.08，水、水泥、砂和石子之間的比例為0.42：1：1.152：2.449；混凝土試件整體尺寸大小為400mm×500mm×500mm，據(jù)此試件相關(guān)參數(shù)建立混凝土剖面--損傷--傳感器二維截面模型如圖1所示.

圖1 細觀混凝土仿真模型Fig.1 Meso-scale concretemodel

將模型導入有限元軟件進行參數(shù)設(shè)定，模型尺寸為500mm×400mm，設(shè)定損傷位置為圓形通孔損傷，中心位置在(225,265)mm處，半徑為r=30mm.假設(shè)傳感器為點源，坐標(150,400)mm至(300,400)mm每間隔12.5mm設(shè)置一個傳感器,編號為1～13.模型主要材料參數(shù)如表1所示.

表1 混凝土仿真模型主要參數(shù)Table 1 Main propertiesof the concretemodel

3 數(shù)值仿真及損傷成像

圖2 激勵信號Fig.2 Excication signal

3.1 仿真成像

(2)因為所以AE=OF=FB′．所以EO=EB′．因而點O在⊙E上．所以從而∠BAC=∠BED．所以DE∥CA．

將建立的混凝土隨機骨料模型導入COMSOL Multiphysics 5.1有限元軟件建立如圖1所示的混凝土剖面--損傷--傳感器模型，選用7號傳感器作為激勵源進行單點激發(fā)多點接收模式.激勵信號為Hanning窗調(diào)制的五周期正弦信號，如圖2(a)所示，以等效力形式加載;圖2(b)為激勵信號對應(yīng)的頻譜.由于縱波波長最大為27mm，為保證數(shù)值計算的從而出現(xiàn)匯聚點.

但由于正向檢測過程中提取的損傷散射信號較直達波、邊界反射波微弱，再加上混凝土結(jié)構(gòu)對波形衰減嚴重以及邊界臨近傳感器加載時信號會在邊界及路徑中出現(xiàn)不同形式的混疊，損傷位置之外在邊界以及路徑中也出現(xiàn)了像素亮點，致使全時域內(nèi)成像結(jié)果并不理想.因此，本文提出基于小波變換的單時刻波場成像法，選取波形干涉出現(xiàn)能量峰值時刻的全場波形來進行損傷成像，所選基函數(shù)為Gabor小波,該基函數(shù)可以同時提供時域和頻域局部化的信息[34].穩(wěn)定性與準確性，故設(shè)置模型網(wǎng)格單元尺寸不大于3×10-4m，時間步長為0.5μs.

激勵源進行信號激勵形成正向波場和損傷散射回波，全部傳感器采集信號.波形中主波包主要由直達波、損傷散射回波、邊界反射波組成，根據(jù)波速以及模型幾何尺寸關(guān)系，截取僅包含以直達波和損傷散射回波為主波形的信號(t∈[0μs,150μs]).

逆時成像過程選用等效彈性參數(shù)模型，幾何尺寸與隨機骨料模型相同，采用軟件ANSYS 14.0建立，依據(jù)文獻[26-27]進行等效參數(shù)的獲取，縱波波速vp,eff=4030.5m/s、橫波波速vs,eff=2586.3m/s以及密度ρeff=2600kg/m3.網(wǎng)格尺寸為2mm，波場傳播時間與信號時間窗值均為T=150μs.時反處理后的13組信號在相應(yīng)幾何位置同時加載形成時反波場.分別得到時反位移場、時反應(yīng)力場和時反能量場信息.

圖3(a)～圖3(c)分別為依據(jù)式(6)～式(8)的成像結(jié)果，圖中每個像素點的值為該點上對應(yīng)全時域內(nèi)的場量最大值；由于時反信號中損傷回波形成的波包在同時刻匯聚于損傷位置形成波前干涉疊加，

圖3 全時域成像Fig.3 Whole time domain imaging

基于小波變換的峰值時刻波場成像法的主要步驟如下：(1)首先對波場中各檢測點采用小波變換提取加載主頻率f0=150kHz對應(yīng)的信號幅值，基于小波變換后的全場數(shù)據(jù)進行最大幅值出現(xiàn)時刻te,max的識別.如圖4(a)所示，橫坐標表示時間，時域T=150μs，步長為0.5μs；縱坐標表示每一時刻對應(yīng)全場的能量最大值，最大峰值左側(cè)出現(xiàn)的第一個峰值是由波形在邊界處的混疊造成，全場最大能量出現(xiàn)時刻te,max=109.2μs;(2)獲取該時刻對應(yīng)原始數(shù)據(jù)波場圖.本文中為了提高計算效率采用間距2mm的網(wǎng)格點，成像過程為提高精度采用插值處理.如圖4(b)所示，圖中顯示了波形在該時刻的干涉疊加，同時也出現(xiàn)了雜波干擾，圖像聚焦點對應(yīng)實際坐標(223.9,289.5)mm；(3)對峰值時刻小波變換數(shù)據(jù)進行成像,如圖4(c)所示，圖像聚焦點對應(yīng)實際坐標(224.8,291.0)mm.

在正向檢測過程中損傷回波信號是由于損傷外表面反射引起，由結(jié)構(gòu)邊界特定位置布置的傳感器采集得到.損傷表面被認為是二次波源，時反過程中信號在波源處干涉疊加，實際損傷反射點坐標為(225,295)mm；檢測結(jié)果與實際損傷位置相差4.005mm，根據(jù)式(10)進行誤差表征，誤差為10.01%，成像結(jié)果準確地匹配了損傷上表面位置，與原始數(shù)據(jù)波場圖相比，小波變換后的數(shù)據(jù)成像結(jié)果消除了雜波干擾，成像結(jié)果更加清晰.

圖4 能量場成像圖Fig.4 Image of theenergy fiel

圖5 不同等效參數(shù)成像結(jié)果Fig.5 Damage imaging of di ff erentequivalentparameters

3.2 關(guān)于等效參數(shù)取值及成像精度的討論

原始模型為非均質(zhì)結(jié)構(gòu)，與實際結(jié)構(gòu)相比，均勻化等效彈性參數(shù)的取值可能引起損傷定位誤差.現(xiàn)依據(jù)文獻[23]的方法討論等效波速取值對損傷定位精度的影響，對波速等效分別設(shè)置為：(1)vp=vp,eff+100m/s,vs=vs,eff+100m/s；(2)vp=vp,eff-100m/s，vs=vs,eff-100m/s.

圖5(a)為設(shè)置一的模擬結(jié)果，時反波場數(shù)據(jù)經(jīng)小波變換后能量峰值出現(xiàn)時刻為te,max=114.5μs，該時刻成像聚焦點坐標為(222.7,278.9)mm，實際損傷位置間距為16.1mm，據(jù)式(10)計算誤差為40.25%；圖5(b)為設(shè)置二的模擬結(jié)果，能量峰值出現(xiàn)時刻為te,max=111.8μs，該時刻成像聚焦點坐標為(221.7,297.4)mm，實際損傷位置間距為4.08mm，誤差為10.2%.結(jié)果顯示,損傷定位點并沒有出現(xiàn)大的偏移，邊界傳感器一定程度上平衡了參數(shù)取值的差異性誤差.

本文研究的時間逆轉(zhuǎn)損傷成像方法以檢測骨料粒徑尺寸范圍內(nèi)的損傷為目的，適用于“不密實區(qū)”及“空洞”類型缺陷的檢測；此類缺陷一般是由于石子架空現(xiàn)象造成，即損傷尺寸與骨料石子尺寸相近，屬于微小型缺陷.

本文選用的混凝土結(jié)構(gòu)骨料粒徑為10mm～40mm，故縮小損傷直徑d為40mm和20mm，成像結(jié)果見圖6.圖6(a)焦點坐標為(226.1,274.1)mm，對應(yīng)損傷實際坐標為(225,285)mm，間距為10.95mm，誤差為27.3%.圖6(b)焦點坐標為(224.1,260.1)mm，對應(yīng)損傷實際坐標為(225,275)mm，間距為15.02mm，誤差為37.5%.

從圖6可看出，損傷成像結(jié)果精度隨損傷尺寸的減小而降低，當損傷尺寸小于骨料平均尺寸(27.63mm)時，成像焦點識別度降低(圖6(b))，當損傷尺寸小于骨料最小尺寸時,成像結(jié)果將無法識別.然而對于非均質(zhì)混凝土結(jié)構(gòu)中“空洞”型損傷的檢測，本文的定位結(jié)果已滿足損傷定位需求.由此證明了時間逆轉(zhuǎn)模型損傷檢測法及基于小波變換的單時刻成像方法對于具有復雜結(jié)構(gòu)的混凝土材料的損傷檢測方法依然具有較好的適用性,且操作簡便、結(jié)果清晰.

圖6 不同損傷尺寸成像結(jié)果Fig.6 Damage imaging of di ff erentdiameters

4 結(jié)論

本文基于細觀混凝土模型，結(jié)合時間逆轉(zhuǎn)成像法進行損傷檢測；通過正向檢測和逆時成像兩個過程完成損傷回波信號的提取和時反加載.時反加載在建立的與細觀模型相同尺寸的等效模型中進行，損傷回波信號在損傷處發(fā)生干涉疊加；通過提取時反過程波場數(shù)據(jù)，并進一步結(jié)合小波變換法提取信號主頻率對應(yīng)的能量幅值，最終對基于小波變換后的全場數(shù)據(jù)進行最大幅值出現(xiàn)時刻的識別并獲取該時刻對應(yīng)原始波場以及小波變換能量場，完成成像，成像結(jié)果清晰且成功匹配了損傷表面的幾何位置.數(shù)值計算結(jié)果表明，時間逆轉(zhuǎn)模型對復雜混凝土結(jié)構(gòu)的損傷定位具有強的適用性，基于小波變換的單時刻成像法能夠消除噪聲干擾提高成像精度.本文采用依據(jù)真實試件的骨料級配曲線建立細觀混凝土模型提高了研究的實際應(yīng)用價值，提出的適用于復雜混凝土結(jié)構(gòu)的損傷定位方法，對混凝土結(jié)構(gòu)內(nèi)部損傷檢測的工程實踐具有較強的指導及應(yīng)用價值.

1萬征，姚仰平，孟達.復雜加載下混凝土的彈塑性本構(gòu)模型.力學學報，2016,48(5)：1159-1171(Wan Zheng,Yao Yangping,Meng Da.An elastoplastic constitutivemodel of concrete under complicated load.Chinese JournalofTheoreticaland Applied Mechanics,2016,48(5)：1159-1171（in Chinese))

2 Wittmann FH,Roelfstra PE.Simulation and analysis of composite structures.Materials Science and Engineering,1984,68(2)：239-248

3 Schlangen E,Van M ier JGM.Simple latticemodel for numerical simulation of fracture of concretematerials and structures.Materialsand Structures,1992,25(9)：534-542

4 W riggers P,Moftah SO.Mesoscalemodels for concrete：homogenization and damage behaviour.Finite Elements in Analysis and Design,2006,42(7)：623-636

5 Leite JPB,Slow ik V,Apel J.Computationalmodel ofmesoscopic structure of concrete for simulation of fracture processes.Comput-ersand Structures,2007,85(17/18)：1293-1303

6 Wang ZM,Kwan AKH,Chan HC.Mesoscopic study of concrete I：generation of random aggregate structure and finit elementmesh.Computers&Structures,1999,70(5)：533-544

7馬懷發(fā),陳厚群,黎保琨.混凝土試件細觀結(jié)構(gòu)的數(shù)值模擬.水利學報,2004,35(10)：27-35(Ma Huaifa,Chen Houqun,Li Baokun.Meso-structure numerical simulation of concrete specimens.JournalofHydraulic Engineering,2004,35(10)：27-35(in Chinese))

8楊華,李宗利,惠弘毅.基于隨機骨料模型的混凝土彈性模量預(yù)測研究.長江科學院院報,2016,33(2)：100-105(Yang Hua,Li Zongli,HuiHongyi.Prediction ofelasticmodulusof concretebased on random aggregatemodel.JournalofYangtze River Scientifi Research Institute,2016,33(2)：100-105(in Chinese))

9 Leslie JR,CheesmanWJ.An ultrasonicmethod of studying deterioration and cracking in concrete structures.Journalofthe American Concrete Institute，1949,21(1)：17-36

10 Clayton D,Barker A,A lbright A,et al.Improved synthetic aperture focusing technique resultsof thick concrete specimens through frequency banding//Chimenti DE,Bond LJ.AIP Conference Proceedings.AIPPublishing,2016,1706(1)：020012

11 FendtKT,Mooshofer H,Rupitsch SJ,etal.Ultrasonic defect characterization in heavy rotor forgingsbymeansof the synthetic aperture focusing technique and optimizationmethods.IEEE Transactions on Ultrasonics,Ferroelectrics,and Frequency Control,2016,63(6)：874-885

12 BeniwalS,Ghosh D,GanguliA.Ultrasonic imaging of concreteusing scattered elastic wave modes.NDT&E International,2016,82(1)：26-35

13 Liu Z,Sun K,Song G,etal.Damage localization in alum inum plate w ith compact rectangularphased piezoelectric transducerarray.MechanicalSystemsand SignalProcessing,2016,70(1)：625-636

14 Vassholz M,Koberstein-Schwarz B,Ruhlandt A,etal.New X-ray tomography method based on the 3d Radon transform compatible w ith anisotropic sources.Physical Review Letters,2016,116(8)：088101-1

15李飛,余西龍,林鑫.基于TDLAS的層析成像技術(shù)TDLAT.力學學報,2014,46(1)：55-59(LiFei,Yu Xilong,Lin Xin.Parallel-beam tomography based on TDLAS.Chinese Journal ofTheoretical and Applied Mechanics,2014,46(1)：55-59(in Chinese))

16邱雷,袁慎芳,蘇永振.基于Shannon復數(shù)小波和時間反轉(zhuǎn)聚焦的復合材料結(jié)構(gòu)多源沖擊成像定位方法.航空學報,2010,31(12)：2417-2424(Qiu Lei,Yuan Shenfang.Multiple impactsource imaging and localization on composite structure based on shannon complex waveletand time reversal focusing.Acta Aeronautica EtAstronautica Sinica,2010,31(12)：2417-2424(in Chinese))

17王強,袁慎芳.主動Lamb波結(jié)構(gòu)健康監(jiān)測中信號增強與損傷成像方法.航空學報,2008,29(4)：1061-1067(Wang Qiang,Yuan Shenfang.Amplifying signal and imaging damagemethod for active lamb-wave structure health monitoring.Acta Aeronautica Et Astronautica Sinica，2008,29(4)：1061-1067(in Chinese))

18王強,袁慎芳.無參考主動Lamb波結(jié)構(gòu)損傷時反成像監(jiān)測方法.航空學報,2010,31(1)：178-183(Wang Qiang,Yuan Shenfang.No baseline time reversal imagingmethod for active lamb wave structuraldamagemonitoring.Acta Aeronautica Et Astronautica Sinica,2010,31(1)：178-183(in Chinese))

19 Zhu R,Huang GL,Yuan FG.Fastdamage imaging using the timereversal technique in the frequency-wavenumber domain.Smart Materialsand Structures,2013,22(7)：075028

20 He J,Harb MS,Yuan FG.Non-contact ultrasonic technique for Lambwave characterization in composite plates.Ultrasonics,2016,64(1)：162-169

21 Chen C,Yuan FG.Impact source identificatio in finit isotropic platesusing a time-reversalmethod：theoreticalstudy.SmartMaterialsand Structures,2010,19(10)：105028

22 Fink M,Cassereau D,Derode A,et al.Time-reversed acoustics.Reportson Progress in Physics,2000,63(12)：1933

23 Saenger EH.Time reverse characterization of sources in heterogeneousmedia.NDT&E International,2011,44(8)：751-759

24 Kocur GK,Saenger EH.Deconvolution of acoustic em issions for source localization using time reversemodeling.Journal ofSound and Vibration,2017,387(1)：66-78

25 Saenger EH,Kocur GK,Jud R,et al.Application of time reverse modeling on ultrasonic non-destructive testing of concrete.Applied MathematicalModelling,2011,35(2)：807-816

26 SaengerEH.Numericalmethods to determinee ff ectiveelastic properties.International Journal ofEngineering Science,2008,46(6)：598-605

27 Saenger EH,Shapiro SA.E ff ective velocities in fracturedmedia：a numerical study using the rotated staggered finite-d ff erence grid.Geophysical Prospecting,2002,50(2)：183-194

28 Steiner B,Saenger EH,Schmalholz SM.Time reversemodeling of low-frequency m icrotremors：application to hydrocarbon reservoir localization.GeophysicalResearch Letters,2008,35(3)：L03307

29 Kocur GK,Saenger EH,Grosse CU,et al.Time reversemodeling of acoustic em issions in a reinforced concrete beam.Ultrasonics,2016,65(1)：96-104

30蔡偉,陳文.復雜介質(zhì)中任意階頻率依賴耗散聲波的分數(shù)階導數(shù)模型.力學學報,2016,48(6)：1265-1280(CaiWei,ChenWen.Fractuionalderivativemodeling of frequency-dependentdissipative mechanism forwave propagation in complexmedia.Chinese JournalofTheoreticaland Applied Mechanics,2016,48(6)：1265-1280(in Chinese))

31唐春安,朱萬成.混凝土損傷與斷裂一數(shù)值試驗.北京：科學出版社,2003(Tang Chunan,ZhuWancheng.A Numerical Experiment of Concrete Damage and Fracture.Beijing：Science Press,2003(in Chinese))

32 Doucet A,Freitas DN,Gordon N.An introduction to sequential MonteCarlomethods.Statistics forEngineering&Information Science,2001,18(1)：3-14

33 Walraven JC,Reinhardt HW.Theory and experiments on themechanical behaviour of cracks in plain and reinforced concrete subjected to shear loading.HERON,1981,26(1A)：1-68

34 Liu Y,Li Z,ZhangW.Crack detection of fibr reinforced composite beams based on continuouswavelet transform.Nondestructive Testing&Evaluation,2010,25(25)：25-44

DAMAGEDETECTIONOFMESO-SCALECONCRETE STRUCTURESBASED ON TIME REVERSALMETHOD1)

Zhao Guoqi2)Qiu Yaping Luo Ying Feng Kan
(Faculty ofCivil Engineering and MechanicsofJiangsu University，Zhenjiang 212013，Jiangsu，China)

A specifi time reversal imaging method is proposed in this article to detect defects in concrete structures.In order to detect the damage that the scale is the same as the aggregate,ameso-scale concretemodel is introduced in thisarticle.As the concrete isa compositematerial composed by cement,aggregate,waterand concrete admixtures,the Monte Carlo random modeland the aggregate grading curve of real concrete samples are introduced for designing this finit elementmodel.Then,thedamagedmodelwasanalyzed by employing aself-adapted time reversedmodel to achieve the ultrasonicwave fiel simulation.This imagingmethod contains two steps：the firs is the forward detection.A series of reflecte echo signalsw ith damage information are obtained in this section.These

signals can be reversed in Matlab to serve as the incident signals in the next Time Reversal process;the second step is to image the damage location via interfering thewave-frontsactuated by di ff erent transducers to illustrate the peaks ofwaveform amplitudes.By determ ining the interfering wave peak time,obtaining the originalwave fiel of thatmoment,we can form ing the wavelet transform energy field and then complete the damage imaging of the concretemodelw ith defect.In the time reversalprocess,we introduced theequivalentelastic parametersas the same geometric dimension of the originalmesosacleconcretemodel to locate thedamagee ff ectively.Comparedw ith theoriginalwave field theenergy fiel modifie by wavelet transform can lower thee ff ectof theenvironmentalnoise.Finally,we discussed theequivalentelastic parameters and the damage sizes to verify the robustnessof thismethod which isapplicable inmonitoring and evaluating the damage in concrete structures.

concrete structure,meso-scalemodel,time reversalmethod,ultrasonic detection,damage imaging

TU317+.8,O347.4+1

A

10.6052/0459-1879-17-007

2017-01-04收稿，2017-03-23錄用,2017-03-24網(wǎng)絡(luò)版發(fā)表.

1)國家自然科學重點國際合作項目(11520101001)、江蘇大學高級人才啟動基金(5501480007)資助.

2)趙國旗，副教授，主要研究方向：結(jié)構(gòu)健康監(jiān)測與檢測新技術(shù)研究及應(yīng)用.E-mail：flagzhao@163.co

趙國旗,仇亞萍,駱英,馮侃.基于細觀混凝土模型的時間逆轉(zhuǎn)損傷成像方法.力學學報,2017,49(4)：953-960

Zhao Guoqi,Qiu Yaping,Luo Ying,Feng Kan.Damage detection ofmeso-scale concrete structures based on time reversalmethod.Chinese JournalofTheoreticaland Applied Mechanics,2017,49(4)：953-960