孫 冰,鄭緒濤,曾 晟,陳振富,郭閃閃
(南華大學(xué)城市建設(shè)學(xué)院,湖南衡陽 421001)
多點布測下錨固缺陷診斷的小波多尺度分析
孫 冰,鄭緒濤,曾 晟,陳振富,郭閃閃
(南華大學(xué)城市建設(shè)學(xué)院,湖南衡陽 421001)
為解決錨桿低應(yīng)變動測分析時,應(yīng)力波速缺乏有效的量化手段的問題,進(jìn)行多點布測方式下錨桿錨固質(zhì)量的室內(nèi)測試實驗,并采用db6小波的3尺度或4尺度分析方法對測試信號進(jìn)行分析。研究表明,多點布測的方式能準(zhǔn)確確定波速的大小和波速傳播過程中的變化大小在頂端傳感器信號中的反應(yīng);對于單一缺陷的錨桿,采用平均波速計算的缺陷位置和錨桿長度與實際結(jié)果吻合,對于多缺陷錨桿,采用分段計算法能避免各測段波速離散性的影響,提高應(yīng)力波檢測精度;波速離散性較大,采用應(yīng)力波速變化值比應(yīng)力波速更能反映錨桿質(zhì)量;錨桿的底端反射時間可在低頻信號a3或a4中反應(yīng)是由于低頻信號反應(yīng)大缺陷、高頻反應(yīng)細(xì)部缺陷的緣故。
多尺度分析;錨固缺陷;多點布測;應(yīng)力波速
隨著錨桿的大規(guī)模應(yīng)用,其質(zhì)量檢測顯得十分重要。在目前錨桿的無損檢測中,大多是利用反射波原理[1]。反射波法雖然在樁基無損檢測中應(yīng)用較為成熟,但由于錨桿與樁基應(yīng)力波的傳播差異[2],使錨桿無損檢測有別于樁基。為了提高應(yīng)力波檢測精度,較為準(zhǔn)確地對錨桿質(zhì)量進(jìn)行評價,國內(nèi)外專家學(xué)者進(jìn)行了許多研究。李義等[3-5]對應(yīng)力波在不同邊界和軸向拉伸荷載作用下的傳播規(guī)律進(jìn)行研究,得出軸向荷載與基頻呈三次冪函數(shù)關(guān)系;提出了用有效錨固長度、基頻、幅值比、衰減系數(shù)等參數(shù)來綜合確定錨桿錨固質(zhì)量優(yōu)劣。汪明武等[6-7]分析了聲頻應(yīng)力波在錨固體系中的反射相位特征和能量衰減變化規(guī)律,用幅值比來確定錨桿狀態(tài)。陳建功等[8]通過對原始信號進(jìn)行小波包分解,得到各頻帶能量特征,并建立神經(jīng)網(wǎng)絡(luò),對錨桿完整性進(jìn)行預(yù)測。夏代林等[9]依據(jù)錨固系統(tǒng)中有錨固缺陷時,不僅能量會發(fā)生變化,相位也會產(chǎn)生突變,提出了將聲波在錨固系統(tǒng)中的能量特征與相位特征相結(jié)合的方法來評定錨桿錨固質(zhì)量。M.D.Beard等[10-11]探討了超聲導(dǎo)波在錨桿無損檢測中的應(yīng)用。Dalhousie大學(xué)及太原理工大學(xué)[12-13]對超聲導(dǎo)波應(yīng)用于錨桿無損檢測進(jìn)行了研究,初步得到了超聲導(dǎo)波的傳播、衰減特性和影響因素。國內(nèi)外研究表明:在確定錨桿缺陷位置時,應(yīng)力波波速是一個十分重要的參數(shù),然而通過公式來準(zhǔn)確計算波速的大小存在難度,還未有準(zhǔn)確有效的量化手段。為此本模型試驗采用了多點布測的方法探討實際波速的大小和變化過程,并通過小波分解和多尺度分析相結(jié)合的方法來確定錨桿的缺陷位置。
小波分析基本思想是將原始信號分解成一系列具有良好的頻域定位性的基元信號,通過對基元信號的各種特征進(jìn)行分析,得到原始信號的局部特征,以達(dá)到對信號進(jìn)行局部化分析的目的。一個信號的小波變換是它在一組小波函數(shù)簇上的投影,而一組小波函數(shù)簇可以由一個小波基函數(shù)通過尺度變換和迭代運算產(chǎn)生,一般稱該小波基函數(shù)為小波函數(shù)。假設(shè)待分析信號為f(t),基小波函數(shù)φ(t)滿足允許性條件
其中,ψ(t)為φ(t)的傅里葉變換;t為時間。待分析信號的小波變換表示為
由式(2)可以看出,小波變換是把待分析信號與一組小波函數(shù)簇做內(nèi)積,將一個時間函數(shù)變換到時間-尺度相平面上,以提取原信號在尺度a(頻率)和位置b(時間)上的某些特征。a,b兩參數(shù)是連續(xù)變化的。原信號f(t)可由其小波變換Wφf(a,b)精確重構(gòu),但它也可以看成基φa,b(t)的分解。但由于a,b是連續(xù)變化的,所以基φa,b(t)之間不是線性無關(guān)的。為了避免冗余,將基離散化構(gòu)成框架,當(dāng)小波函數(shù)的伸縮平移系是正交系時,小波框架便無冗余,這就需要多尺度分析。
原信號S可以分解為低頻a1和高頻d1,低頻a1又可分解為低頻a2和高頻d2,以次類推,這就是Mallat塔式重構(gòu)算法。它說明任何信號可分解成不同頻帶的細(xì)節(jié)之和,隨著分解層數(shù)的不同,這些頻帶互不重疊且充滿整個頻率空間,也就是正交離散小波變換的時頻窗互不重疊、相互鄰接,形成對時頻平面的一種剖分。通過在時域波形上不同頻帶上能量統(tǒng)計與傅里葉變換在頻域上進(jìn)行能量統(tǒng)計不同,這恰好體現(xiàn)了小波分析具有時頻分析能力的優(yōu)勢。
2.1 模型制作
制作3根錨桿,編號分別為M-1,M-2,M-3。錨桿長度2.4 m,鋼筋外露長度為0.1 m,采用先插桿后注漿的施工方式。采用?28螺紋鋼筋,其缺陷的具體設(shè)置方式見表1,錨桿結(jié)構(gòu)簡圖如圖1所示。
表1 模型試驗錨桿設(shè)計參數(shù)簡況Table 1 Design parameters of bolts
2.2 實驗測試儀器
測試系統(tǒng)由動態(tài)信號分析儀AVANT-10、壓電式加速度傳感器、力錘、電荷放大器和電腦處理系統(tǒng)組成。為了避免在采樣時出現(xiàn)頻率折疊現(xiàn)象,采樣頻率至少是帶限信號最高頻率的2倍。
2.3 測點布置及波速計算
圖1 錨桿結(jié)構(gòu)簡圖Fig.1 Structure of bolts
由于波速計算沒有定量的公式,當(dāng)采用一個測點來計算波速時,雖然所需傳感器少,但實際測得的波速為平均波速,不能反映實際波速大小及變化情況。當(dāng)錨桿缺陷較多時,缺陷的準(zhǔn)確位置應(yīng)采用各段的波速及對應(yīng)的時間來計算,而采用平均波速計算時,其大小的選取存在很大的范圍。因此,采用多點布測的方式,達(dá)到將波速量化的目的,并討論了波速的大小和在傳播過程中的變化情況。測點布置位置Xi為第i個傳感器距鋼筋外露端錨桿端部的距離,Xi=0, 0.70,1.30,1.90,2.45 m(i=1,2,3,4,5)。
在應(yīng)力波從X1傳播到X5的過程中,相鄰兩傳感器的信號起跳點會存在時間差Δt,波速vw為
其中,ti為第i個傳感器信號起跳點時間。各測段的波速見表2,各錨桿全部測點原始信號如圖2所示。
表2 各錨桿各測段波速Table 2 W ave velocities of each measurem ent segment of anchors m/s
3.1 各錨桿原始信號
考慮到測量過程中的能量衰減,噪聲在信號中的比例越來越大,在信號尾部噪聲占主要部分,為取得較好的信噪比及缺陷信號得到清晰的反應(yīng),選取前兩個周期進(jìn)行分析,得到各錨桿頂端原始信號(圖3)。
圖2 各錨桿全部測點原始信號Fig.2 The original signal of eachmeasurement point
圖3 各錨桿頂端原始信號Fig.3 The original signals of each anchor
由圖3可以得出,M-1與M-2的底端反射較為明顯,由“峰-峰”算法可得錨桿長度為2.418 m。然而由于M-3內(nèi)部缺陷的反射波干擾,底端反射信號的相鄰時間段內(nèi)存在多個波峰,僅通過原始信號的信號特征不能準(zhǔn)確確定底端反射的時間,但M-2,M-3的原始信號中均出現(xiàn)某一段時間內(nèi)的幅值突然大于前面的相鄰段,有明顯的能量反射特征,M-1信號中個別點不能確定,有可能是受錨桿端頭露出鋼筋段的反射應(yīng)力波影響。
3.2 錨桿原始信號的小波多尺度分析
分別對3根錨桿的原始信號進(jìn)行小波變換,對于小波的選取,考慮到時間-頻率窗的大小、計算的復(fù)雜性、基小波的光滑與對稱性及逼近性等,為此選取Daubechies小波。Daubechies小波為有限緊支撐正交小波,其時域和頻域的局部化能力強(qiáng),在分析過程中能提供更具體、更實際的數(shù)字濾波器。由db N小波函數(shù)和尺度函數(shù)的波形及頻率特性知N值不同,濾波器長度L=2N不同,N越大,尺度函數(shù)和小波函數(shù)波形越光滑,濾波器頻率特性越好。然而其計算量急劇增加,由不確定原理可知,不存在完美的濾波器,所以選取db6小波,噪聲選擇未知白噪聲,使得分解后的信號中不存在處處奇異信號。通過Matlab小波工具箱,變換結(jié)果如圖4~6所示。
從圖4可以得出,由小波變換原理可知原始信號S=a3+d3+d2+d1,a3,d3,d2,d1為分解后的信號,依次對應(yīng)信號的頻率由低到高。a3為低頻信號,波形規(guī)則,表明錨桿無大尺寸缺陷;d1為以高頻低幅信號,其存在有可能對原始信號個別點判斷有影響。在多尺度分析中信號灰度的深淺反映了信號變化強(qiáng)弱,通過對理想信號的小波多尺度分析,不同頻率反映不同大小的缺陷,不同顏色對應(yīng)于不同程度的缺陷,顏色越淺,缺陷越大。因此M-1的多尺度分析表明入射信號強(qiáng),而反射信號弱,錨桿質(zhì)量良好。原始信號波形較為規(guī)則,分解后的信號也沒表現(xiàn)出缺陷特征。長度計算時,底端反射信號無論是在原始信號還是在分解后的信號中都能明顯讀出。
從圖5可知,M-2原始信號出現(xiàn)較為明顯的缺陷特征,由于從缺陷處返回的能量影響,使得首波后的一段時間內(nèi)測點振幅衰減不明顯,對2號錨桿原始信號進(jìn)行3尺度分解,分解后錨桿的底端反射也可由a3中得出,高頻信號d2中有明顯的信號突變,可由此算出缺陷的位置。在多尺度分析圖中,當(dāng)尺度為2時,有明顯的信號變化,根據(jù)突變信號的時間減去首波時間為412μs,由此可得缺陷位置為0.803 m。
圖4 M-1原始信號3尺度分析Fig.4 Wavelet decomposition and three-scale analysis on anchor 1
圖5 M-2原始信號3尺度分析Fig.5 Wavelet decomposition and three-scale analysis on anchor 2
圖6 M-3原始信號多尺度分析Fig.6 Wavelet decomposition and three-scale analysis on anchor 3
從圖6可以看出,M-3原始信號較前2個信號要復(fù)雜一些,其返回的能量較大,使得在較短的時間后振幅大于首波的振幅,表明此時的錨桿質(zhì)量存在嚴(yán)重的缺陷,對原始信號分別進(jìn)行3尺度和4尺度分解,在4尺度分解時,根據(jù)塔式重構(gòu)算法的原理有S= a4+d4+d3+d2+d1。當(dāng)尺度為3時,其表現(xiàn)出的波形特征較尺度為4時要弱。因此用4尺度分析,底端反射信號可在尺度4中得出,由此所得錨桿長度為2.313 m。在確定錨桿長度時也可以用分解后低頻信號a4來計算,缺陷引起的信號突變,在分解后高頻信號d2中有明顯的表現(xiàn),用多尺度分析中的尺度1和尺度來得出兩缺陷的位置分別為0.562,1.174 m。
由圖4~6中可以得出,各原始信號經(jīng)過小波分解后都存在一低幅高頻信號,與低頻信號相比,高頻信號衰減速度較快,傳播距離有限,但是低頻信號對于細(xì)部缺陷不敏感,所以細(xì)部缺陷在離錨桿激發(fā)端較遠(yuǎn)時,結(jié)果有可能不能反應(yīng)出細(xì)部缺陷。雖然分解后的信號對缺陷的反應(yīng)較為準(zhǔn)確,但為了減少人為選取計算點位的誤差,在確定缺陷位置時,較多是用多尺度分析中的相位突變來計算。關(guān)于缺陷的長度,在多尺度分析中由于色條的寬度較窄,選取時易引起較大的誤差,故宜從分解后的信號入手。
3.3 應(yīng)力波速大小及變化特征
由于缺陷的存在,使得錨桿中應(yīng)力波傳播過程中波速變化范圍較大,表2中,由于M-1的測段1~3無缺陷存在,各測段錨桿波速穩(wěn)定,測段4與測段3相比缺少鋼筋,從而使得波速變化較大。M-2的測段1與測段3材料相同,測段2為空漿段,應(yīng)力波從測段1~2或從測段2~3都有明顯的突變。同樣在M-3中測段3由于錨固介質(zhì)加強(qiáng),應(yīng)力波速為5 758 m/s,表明錨固介質(zhì)強(qiáng)度對波速大小有巨大影響,在文獻(xiàn)[14]中有詳細(xì)論述。空漿段段波速比密實段小,對比M-1,M-2,空漿段波速較素混凝土段波速較為接近,說明鋼筋在此段內(nèi)起到作用小,M-3空漿段波速與M-2空漿段波速不同,這是由于加強(qiáng)段的出現(xiàn)使得M-3空漿段邊界條件變化。但M-1,M-2和M-3的平均波速相對穩(wěn)定,分別為4 134, 4 070和4 320 m/s。其計算長度和缺陷位置與實際相比,結(jié)果偏小,而采用算得的平均波速,則與實際較為吻合。
從以上分析結(jié)果可知,不論是計算錨桿總長度,還是測得的缺陷距離桿端的距離,與實際的情況相比都有一定偏差,引起這種情況的原因主要有:
(1)計算時取的波速以及計算點位選取。如本文在計算時所取波速為3 900 m/s,在完好或單個缺陷時計算結(jié)果較好,而在多缺陷時誤差相對較大。通過多點布測的方式,實際傳播過程中,波速較為分散,尤其是缺陷較多時,準(zhǔn)確的計算方式為各段的波速與各段的時間相乘,而實際運用較為困難。根據(jù)文獻(xiàn)[15]中式3.76,由于平均折算剛度值與用混合物法所得的值不同,所以用該公式來計算波速存在困難。在計算點位選取時,由于人為因素,會存在一定的誤差,而波速數(shù)量級較大,由此會產(chǎn)生較大的誤差。
(2)小波分解尺度的選擇由圖6可以看出,不同的分解尺度得到的結(jié)果不同,如尺度選擇不合適,則在多尺度分析圖中不能準(zhǔn)確反應(yīng)錨桿的內(nèi)部缺陷,使分析人員得到錯誤的結(jié)論。
(1)反射波法錨桿動測分析分析的重要內(nèi)容之一就是反射信號的辨認(rèn),其方法就是獲取信號突變的位置。當(dāng)實測信號含有噪聲干擾或?qū)崪y信號比較復(fù)雜時,小波分析方法既能消除噪聲,又能保持信號的突變性。但在消噪時注意消噪層數(shù)的選擇,避免真實信號損失。
(2)小波分析技術(shù)在處理非平穩(wěn)信號時十分有效,但要注意根據(jù)實際情況選取合適的小波基及分解尺度,如相同的材料和尺寸,僅缺陷不同,分解尺度的選擇也可能不同,而它們的選取直接影響分析結(jié)果。
(3)通過多點布測的測量方式可以得出,錨桿缺陷越多時,傳播過程中的波速變化越大,但平均波速相對穩(wěn)定。在采用平均波速時,缺陷越靠近底端,結(jié)果越精確。如果用分段波速,則不存在這個問題,但需要根據(jù)缺陷時突變信號的能量、形態(tài)、長度等參數(shù)建立與波速變化之間的關(guān)系,這方面還有待深入研究。正是由于波速的復(fù)雜性,使低應(yīng)變動測法的應(yīng)用和推廣受到一定約束。
[1] 岳向紅,劉明貴,李 祺.錨桿檢測技術(shù)研究進(jìn)展[J].土工基礎(chǔ),2005,19(3):83-85.
Yue Xianghong,Liu Minggui,Li Qi.Development of bolt’s testing technique[J].Soil Eng.and Foundation,2005,19(3):83-85.
[2] 孫 冰,曾 晟,丁德馨,等.低應(yīng)變動測樁及錨桿中波的傳播規(guī)律差異研究[J].巖土力學(xué),2011,32(4):1143-1148.
Sun Bing,Zeng Sheng,Ding Dexin,et al.Research on transmit rules of stresswavewith low strain in dynamic test pile and anchorage bolt [J].Rock and Soil Mechanics,2011,32(4):1143-1148.
[3] 李 義,劉海峰,王富春.錨桿錨固狀態(tài)參數(shù)無損檢測及其應(yīng)用[J].巖石力學(xué)與工程學(xué)報,2004,23(10):1741-1744.
Li Yi,Liu Haifeng,Wang Fuchun.Nondestructive testing of parameters of bolt anchoring state and its application[J].Chinese Journalof Rock Mechanics and Engineering,2004,23(10):1741-1744.
[4] 任智敏,李 義.基于聲波測試的錨桿錨固質(zhì)量檢測信號分析與評價系統(tǒng)實現(xiàn)[J].煤炭學(xué)報,2011,36(1):191-196.
Ren Zhimin,Li Yi.Analysis of detection signal and realization on evaluation system ofboltanchoring quality based on sound wave testing[J].Journal of China Coal Society,2011,36(1):191-196.
[5] 李 義,張昌鎖,王 成.錨桿錨固質(zhì)量無損檢測幾個關(guān)鍵問題的研究[J].巖石力學(xué)與工程學(xué)報,2008,27(1):108-116.
Li Yi,Zhang Changsuo,Wang Cheng.Study on several key issues in nondestructive detection of bolt bonding integrality[J].Chinese Journal of Rock Mechanics Engineering,2008,27(1):108-116.
[6] 汪明武,王鶴齡,羅國煜,等.錨桿錨固質(zhì)量無損檢測的研究[J].工程地質(zhì)學(xué)報,1999,7(1):72-76.
Wang Mingwu,Wang Heling,Luo Guoyu,et al.Research on non-destructive test of anchoring quality of rock bolts[J].Journal of Engineering Geology,1999,7(1):72-76.
[7] 汪明武,王鶴齡.錨固質(zhì)量的無損檢測技術(shù)[J].巖石力學(xué)與工程學(xué)報,2002,21(1):126-129.
Wang Mingwu,Wang Heling.Nondestructive testing of anchoring quality[J],Chinese Journal of Rock Mechanics and Engineering, 2002,21(1):126-129.
[8] 陳建功,張永興.錨桿系統(tǒng)動測信號的特征分析[J].巖土工程學(xué)報,2008,30(7):1051-1057.
Chen Jiangong,Zhang Yongxing.Analysis on characteristics of dynamic signal for bolt anchorage system[J].Chinese Journal of Geotechnical Engineering,2008,30(7):1051-1057.
[9] 楊天春,吳燕清,夏代林.基于相位推算法的錨桿施工質(zhì)量無損檢測分析方法[J].煤炭學(xué)報,2009,34(5):629-633.
Yang Tianchun,Wu Yanqing,Xia Dailin.An analytic method for rock bolt’snon-destructive testing signalsby phase deductedmethod [J].Journal of China Coal Society,2009,34(5):629-633.
[10] Beard M D.Guided wave inspection of embedded cylindrical structures[D].London:University of London,2002.
[11] Beard M D,LoweM JS.Non-destructive testing of using guided ultrasonic waves[J].International Journal of Rock Mechanics and Mining Sciences,2003,40:527-536.
[12] Zou D H,CuiY,Madenga V,et al.Effects of frequency and grouted length on the behaviorof guided ultrasonicwaves in rock bolts[J].International Journal of Rock Mechanics&Mining Sciences,2007, 40:813-819.
[13] Zhang C S,Zou D H,Madenga V.Numerical simulation of wave propagation in grouted rock bolts and the effects ofmesh density and wave frequency[J].International Journalof Rock Mechanics& Mining Sciences,2006,43:634-639.
[14] 孫 冰,曾 晟,丁德馨,等.錨固系統(tǒng)與錨固介質(zhì)材料中應(yīng)力波波速的關(guān)系[J].防震減災(zāi)工程學(xué)報,2013,33(3):305-310.
Sun Bing,Zeng Sheng,Ding Dexin.Relationships between stress wave velocities in anchorage system and material of anchoragemedium[J].Journal of Disaster Prevention and Mitigation Engineering,2013,33(3):305-310.
[15] 周培基,霍普肯斯A K.材料對強(qiáng)沖荷載的動態(tài)響應(yīng)[M].北京:科學(xué)出版社,1985.
M ulti-scale analysis on anchoring defects diagnosis under multip lemeasuring points
SUN Bing,ZHENG Xu-tao,ZENG Sheng,CHEN Zhen-fu,GUO Shan-shan
(Institute ofUrban Construction,University of South China,Hengyang 421001,China)
In order to solve the problem of the lack of effective quantitativemethod of stresswave velocity,indoor test was carried outwith multi-point distribution measurement,and the db6 wavelet scale 3 or 4 dimension analysismethod were adopted to analyze the test signal.This study show that the size of thewave velocity can be determined accurately aswell as the signal of the change ofwave velocity can be responded in the top of the sensor by using themulti-point distribution measurement.For the anchor boltwith a single defect,the defect position is obtained by using the average wave velocity and the calculated anchor length are consistentwith the actual results,while adopting subsection calculation method for themore defects anchor to avoid the influence of discrete wave velocity of each measurement section.Compared with the stresswave velocity,its change can reflect the quality of anchor bolt better due to larger discreteness,based on that the bottom reflection time can be presented in the low frequency signal as a3or a4,since that the low frequency signals big defects while high ones represent the detail.
multi-scale analysis;anchoring defects;multi-scale analysis;stress wave velocity
TD35
A
0253-9993(2014)07-1385-06
孫 冰,鄭緒濤,曾 晟,等.多點布測下錨固缺陷診斷的小波多尺度分析[J].煤炭學(xué)報,2014,39(7):1385-1390.
10.13225/ j.cnki.jccs.2013.1122
Sun Bing,Zheng Xutao,Zeng Sheng,et al.Multi-scale analysis on anchoring defects diagnosis undermultiplemeasuring points[J].Journal of China Coal Society,2014,39(7):1385-1390.doi:10.13225/j.cnki.jccs.2013.1122
2013-08-04 責(zé)任編輯:張曉寧
國家自然科學(xué)基金資助項目(51204098);南華大學(xué)“十二五”科技創(chuàng)新團(tuán)隊項目資助項目(2012NHCXTD12);湖南省重點學(xué)科建設(shè)資助項目
孫 冰(1979—),女,河南平頂山人,副教授,博士。E-mail:sunbingzs@126.com