袁 野，黃曉寒，王仲剛，張 彧

(后勤工程學(xué)院 軍事土木工程系，重慶 401311)

基于頻響函數(shù)虛部奇異值熵的混凝土內(nèi)部缺陷檢測

袁 野，黃曉寒，王仲剛，張 彧

(后勤工程學(xué)院 軍事土木工程系，重慶 401311)

根據(jù)相干函數(shù)和力錘激勵的有效頻段來確定頻響函數(shù)的研究頻段，再將頻響函數(shù)虛部與信息熵理論相結(jié)合，計算頻響函數(shù)虛部構(gòu)造的Hankel矩陣的奇異值熵，最后通過概率法計算出異常情況的判斷值來實現(xiàn)對混凝土內(nèi)部缺陷的檢測。在此基礎(chǔ)上提出了一種簡單實用的混凝土缺陷成像的定性分析方法，首先利用各掃射射線的奇異譜熵值將試件分為最有可能產(chǎn)生缺陷與最不可能產(chǎn)生缺陷的區(qū)域，再計算各缺陷單元格的投影函數(shù)值，非缺陷單元格的投影函數(shù)值置0，最后用各單元格的投影函數(shù)值的等高線來反應(yīng)混凝土內(nèi)部情況，并通過試驗驗證了該方法的有效性。

奇異值熵；頻響函數(shù)；相干函數(shù)；缺陷

工程結(jié)構(gòu)發(fā)生損傷時會改變結(jié)構(gòu)的動力特性和動力響應(yīng)，因此可以根據(jù)結(jié)構(gòu)的動力特性和動力響應(yīng)的改變檢測出結(jié)構(gòu)中存在的損傷。其中基于結(jié)構(gòu)的頻率、振型等模態(tài)參數(shù)的改變進行損傷檢測的方法得到了廣泛的應(yīng)用，但這些方法屬于間接測試方法，其存在由于模態(tài)提取誤差等因素而降低檢測精度的問題，同時存在模態(tài)參數(shù)受噪聲影響較大等問題。為了克服模態(tài)參數(shù)法的這些缺點，基于振動響應(yīng)測試的結(jié)構(gòu)損傷檢測的頻響函數(shù)法受到了研究人員的關(guān)注[1]。楊彥芳[2]等提出了以頻響函數(shù)作為損傷識別的基本參量，利用主元分析和多元控制圖來識別網(wǎng)架結(jié)構(gòu)損傷的方法。高海洋，郭杏林等[3]僅利用損傷后的板結(jié)構(gòu)頻率響應(yīng)函數(shù)虛部，建立二維板結(jié)構(gòu)的損傷定位指標(biāo)。方劍青，矯桂瓊[4]用神經(jīng)網(wǎng)絡(luò)BP算法直接對所測結(jié)構(gòu)的頻響函數(shù)進行分析，用來識別結(jié)構(gòu)的狀態(tài)。頻響函數(shù)法在實際應(yīng)用中具有所使用的傳感器及設(shè)備較少，測試方便[5], 不需要專門的激振設(shè)備[6]等優(yōu)點；頻響函數(shù)的虛部更是具有良好的損傷敏感性和抗噪能力[3]。

SHANNON借鑒熱力學(xué)熵的概念，結(jié)合數(shù)理統(tǒng)計理論提出了能夠用來量化信號的復(fù)雜程度的信息熵。當(dāng)結(jié)構(gòu)出現(xiàn)損傷時，其振動狀態(tài)的復(fù)雜性會發(fā)生相應(yīng)的變化，因此可用反映振動狀態(tài)復(fù)雜程度的熵來評價結(jié)構(gòu)當(dāng)前的狀態(tài)[7]。謝中凱等[8]基于近似熵理論，利用自由振動信號對鋼筋混凝土梁結(jié)構(gòu)進行了損傷識別。劉國華等[9]基于改進的傳遞熵理論，對混凝土簡支梁進行了損傷識別。

因此筆者提出了一種利用頻響函數(shù)虛部建立Hankel矩陣，通過計算奇異值熵值來對混凝土內(nèi)部缺陷進行檢測的方法，并在此基礎(chǔ)上提出了一種定性分析缺陷的方法。

1 頻響函數(shù)的理論分析

頻響函數(shù)不僅包含了信號的幅值、頻率和相位的信息，而且還可直接獲取，所以頻響函數(shù)作為判斷結(jié)構(gòu)損傷的指標(biāo)是比較好的。對于單輸入單輸出的線性系統(tǒng)，其頻響函數(shù)定義為輸出y(t)的傅里葉變換與輸入x(t)的傅里葉變換之比，即:

(1)

然而在實際測量中，為消除響應(yīng)噪聲信號的干擾，頻響函數(shù)通過響應(yīng)和激勵間的自譜實現(xiàn)，表達式為:

(2)

式中：Syy(ω)為響應(yīng)的功率譜密度;Sxx(ω)為激勵的功率譜密度。

2 相干函數(shù)頻段的選取

圖1所示為考慮噪聲輸入的單輸入與單輸出系統(tǒng)。

圖1 考慮噪聲輸入的單輸入/單輸出系統(tǒng)

圖1中p(t)與q(t)分別表示真實輸入與輸出，輸入和輸出的測量噪聲分別為n1(t)和n2(t)，則實際輸入和輸出的測量值為x(t)=p(t)+n1(t)和y(t)=q(t)+n2(t)。假設(shè)噪聲n1(t)和n2(t)之間互不相干且噪聲與輸入輸出之間互不相關(guān)。

激勵與響應(yīng)的譜相干函數(shù)定義為：

(3)

在有噪聲的情況下，信號自譜、源信號自譜和噪聲自譜的關(guān)系為：

(4)

Syy(ω)=Sqq(ω)+Sn2n2(ω)(5)

Sxy(ω)=Spq(ω)(6)

式中：Sxx(ω)為輸入端信號自譜;Syy(ω)為輸出端信號自譜;Spp(ω)為輸入端源信號自譜;Sqq(ω)為輸出端源信號自譜;Sn1n1(ω)為輸入端噪聲信號自譜;Sn2n2(ω)為輸出端噪聲信號自譜。Sxy(ω)為輸入與輸出信號自譜;Spq(ω)為輸入與輸出源信號自譜。

相干函數(shù)化為：

(7)

相干函數(shù)反映了各頻率的信噪比，源信號在某頻段的信噪比越大，這個頻段的相干函數(shù)也越大。因此相干函數(shù)在頻域上反應(yīng)了2個信號的相關(guān)程度，根據(jù)其幅值可以估計出源信號所在的頻段[10]。選取信噪比高的頻段進行頻響函數(shù)的分析有利于得到更準確的結(jié)果。

圖2 力錘敲擊產(chǎn)生的脈沖激勵及其頻譜圖

3 激勵頻段的選取

力錘錘擊時，傳遞給結(jié)構(gòu)的沖擊力近乎為半正弦形，如圖2(a)所示，其頻譜如圖2(b)所示。fc為截止頻率，fc以下的頻段為有效激勵頻段，fc與脈沖寬度τc有一定的關(guān)系[11]：

(8)

綜合相干函數(shù)與錘擊激勵得到的各自有效頻段，試驗采用03 500 Hz頻段上的頻響函數(shù)虛部進行分析。

4 Hankel矩陣奇異譜熵

對測量得到的輸入輸出數(shù)據(jù)進行頻響分析，提取出上述研究頻段內(nèi)的頻響函數(shù)虛部的數(shù)據(jù)，將該數(shù)據(jù)組成m×n階的Hankel矩陣：

(9)

式中：m+n-1=N，m≥n。再對H矩陣進行奇異值分解，得到奇異值X0，令s=min(m,n)。奇異值熵定義為：

(10)

pi=δi(∑si=1δi)

(11)

式中：δi為第i個奇異值;pi為第i個奇異值在整個奇異值序列中所占的比重。

當(dāng)結(jié)構(gòu)發(fā)生損傷時，奇異值熵也會發(fā)生相應(yīng)的變化，損傷位置處的奇異值熵較無損傷位置處的奇異值熵大，因而可以用其來表征結(jié)構(gòu)的狀態(tài)和損傷情況[7]。

5 概率法

對于波速的判斷，南京水利科學(xué)研究院的羅騏先提出了概率法，概率法認為隨機誤差引起的混凝土的質(zhì)量波動是符合正態(tài)分布的[12]。將波速按大小順序進行排列x1≥x2≥…≥xs，然后對數(shù)據(jù)進行篩選，將位于后面明顯較小的數(shù)據(jù)視為可疑數(shù)據(jù)，再將這些可疑數(shù)據(jù)中最大的一個與排在其前面的數(shù)據(jù)組成一組新數(shù)據(jù)，對該組數(shù)據(jù)求其平均值mx和標(biāo)準差Sx，并按下式求出異常情況的判斷值xn>X0：

(12)

式中：λ為根據(jù)規(guī)范JGJ106-2003《建筑樁基檢測技術(shù)規(guī)范》查得的系數(shù)。

得到異常數(shù)據(jù)判斷臨界值后，將該值與可疑數(shù)據(jù)的最大值xn進行比較，若xnX0，則表示xn是正常值，則需對前n+1個數(shù)據(jù)進行上述的計算和判斷。基于此方法，筆者對計算得到的奇異譜熵數(shù)據(jù)進行類似的處理，根據(jù)找到的異常值來判斷混凝土內(nèi)部的缺陷。

6 定性分析方法

在樁基檢測中，可以利用聲陰影重疊法在初定缺陷區(qū)域內(nèi)仔細判定缺陷的確切位置、范圍和性質(zhì)[13]。類似于聲陰影重疊法，在2個方向進行檢測，利用奇異譜熵值分別劃出陰影區(qū)，則2個陰影區(qū)內(nèi)邊界線交叉重疊所圍成的區(qū)域，即為缺陷的范圍。為了量化陰影區(qū)，首先將試件截面離散成若干網(wǎng)格化單元，然后將射線所經(jīng)過的單元都賦予該條射線計算得到的奇異譜熵值，最后將網(wǎng)格中的奇異譜熵值進行疊加即為該單元格的投影函數(shù)值。

以一個二維模型示意圖(見圖3)來說明，模型中央設(shè)一正方形缺陷，在2個方向進行檢測，通過奇異譜熵值初判得到2個陰影區(qū)。對于經(jīng)過缺陷區(qū)域的某單元格，奇異譜熵值較大，陰影重疊越多，則穿越其途經(jīng)缺陷區(qū)域的射線越多，進而該單元投影函數(shù)值越大。由圖可知，則奇異譜熵值：A>B>C，由此判斷單元格A最有可能處于缺陷區(qū)，這與實際情況相符。

圖3 缺陷定性分析的二維模型示意

7 方形混凝土柱缺陷檢測

圖4 試件圖片

設(shè)計并制作了規(guī)格(長X寬X高)為800 mm×800 mm×150 mm的素混凝土柱試件，該缺陷試件內(nèi)部缺陷尺寸(長X寬X高)為160 mm×160 mm×150 mm，如圖4所示。試驗的儀器設(shè)備為DH5960超高速動態(tài)信號采集儀、KDL力錘和帶磁座的加速度傳感器。

7.1 模型Ⅰ

采用一發(fā)一收對測方式采集信號(均采用同一傳感器接收信號)，設(shè)置8對測點，其中第4、5測試點是經(jīng)過缺陷的情況，在每一測點敲擊6次，測點布置如圖5所示，試驗?zāi)Ｐ廷袢鐖D6所示。

利用測得的輸入與輸出數(shù)據(jù)進行頻響函數(shù)虛部奇異值熵的計算，計算結(jié)果如表1所示。

表1中數(shù)據(jù)按概率法進行計算與判斷，最終計算出來的臨界值X0=3.930 4。將其與表中數(shù)據(jù)進行對比發(fā)現(xiàn)，在4-4中只有一組數(shù)據(jù)小于臨界值，5-5中的所有數(shù)據(jù)均大于臨界值，故判斷4-4與5-5測試條件下為經(jīng)過缺陷的情況，與實際結(jié)果相吻合。

圖5 模型Ⅰ測點布置

圖6 試驗?zāi)Ｐ廷袷疽?/p>

7.2 模型Ⅱ

采用一發(fā)多收方式采集信號，設(shè)置8個采集點，在激發(fā)點敲擊6次，測點布置如圖7所示，試驗?zāi)Ｐ廷蛉鐖D8所示。

圖7 模型Ⅱ測點布置

圖8 試驗?zāi)Ｐ廷蚴疽?/p>

利用測得的輸入與輸出數(shù)據(jù)進行頻響函數(shù)虛部奇異值熵的計算，結(jié)果見表2(限于篇幅，表中只給出其中一側(cè)的計算結(jié)果)。

表1 模型Ⅰ試驗的奇異值熵計算結(jié)果

表2 模型Ⅱ試驗的奇異值熵計算結(jié)果

對表2中數(shù)據(jù)按概率法進行計算與判斷，最終計算出來的臨界值X0=3.546 2。將其與表中數(shù)據(jù)進行對比發(fā)現(xiàn)，可以初判激勵-3與激勵-6圍成的陰影區(qū)為可疑缺陷區(qū)域，再通過另一個方向的測試，得到陰影重疊區(qū)即為最有可能產(chǎn)生缺陷的區(qū)域，計算重疊區(qū)域中網(wǎng)格的疊加投影函數(shù)值(奇異譜熵值)，其余網(wǎng)格中的投影函數(shù)值置0，做等高線圖，如圖9所示。

圖9 試驗?zāi)Ｐ廷虻钠娈愔奠氐雀呔€圖

通過上述步驟，可以重建試件的截面圖像，筆者提出的定性分析方法能正確地反映出缺陷的位置，同時也能大致反映出缺陷的區(qū)域范圍，將定性分析方法與定量反演計算相結(jié)合，可以減少誤判，提高缺陷識別的準確性和可靠性。

8 結(jié)論

通過構(gòu)建合適頻段上的頻響函數(shù)虛部Hankel矩陣，并對該矩陣進行奇異值分解，計算奇異值熵，最后用概率法對缺陷進行判別。通過奇異譜熵值的計算，再結(jié)合提出的定性分析方法，能夠?qū)θ毕莸奈恢煤头秶M行比較準確的估計，可以作為定量反演的輔助判別。所提方法簡單直觀、產(chǎn)生的誤差較小、辨識度較高，而且適用于低信噪比的情況;同時也不需要結(jié)構(gòu)損傷前的數(shù)據(jù)作為參考，更利于實際工程的應(yīng)用。

[1] CARDEN E P, FANNING P. Vibration based condition monitoring:a review[J].Structural Monitoring,2004,3(4):355-377.

[2] 楊彥芳，宋玉普. 基于主元分析和頻響函數(shù)的網(wǎng)架結(jié)構(gòu)損傷識別方法[J]. 工程力學(xué)，2007，24(9):105-110.

[3] 高海洋，郭杏林，吳明勇.基于頻響函數(shù)虛部的板結(jié)構(gòu)損傷檢測方法研究[J].振動與沖擊，2012,31(12):86-91.

[4] 方劍青,矯桂瓊. 基于頻響函數(shù)的結(jié)構(gòu)狀態(tài)識別神經(jīng)網(wǎng)絡(luò)方法[J]. 振動、測試與診斷, 2007,27(1):45-47.

[5] WANG Z,LIN R M,LIM M K. Structural damage detection using measured FRF data[J]. Computer Methods in Applied Mechanics and Engineering,1997,147:187-197.

[6] YAM L H,LEUNG T P,LI D B, et al. Use of ambient response measurements to determine dynamic characteristics of slender structures[J]. Engineering Structures,1996,18(3):145-150.

[7] 李睿,于德介,曾威. 一種基于奇異譜熵和脈沖響應(yīng)的結(jié)構(gòu)損傷診斷方法[J]. 振動工程學(xué)報,2006,19(3):331-335.

[8] 謝中凱,劉國華. 近似熵在混凝土結(jié)構(gòu)損傷識別中的應(yīng)用[J]. 浙江大學(xué)學(xué)報(工學(xué)版),2013,47(3):456-463.

[9] 劉國華,謝中凱. 一種用于檢測梁結(jié)構(gòu)損傷的改進傳遞熵[J]. 振動工程學(xué)報,2014,27(1):136-144.

[10] 錢濟國. 頻響函數(shù)測量的脈沖激勵法[J]. 振動、測試與診斷,1994,14(2):32-39.

[11] 吳佳曄. 土木工程檢測與測試[M]. 北京:高等教育出版社,2015.

[12] 羅騏先.樁基工程檢測手冊[M].北京:人民交通出版社，2003:181-258.

[13] 吳慧敏.結(jié)構(gòu)混凝土現(xiàn)場檢測新技術(shù)——混凝土非破損檢測[M].長沙：湖南大學(xué)出版社，1998:161-164.

The Detection of Internal Defects of Concrete Based on the Imaginary Part of the Frequency Response Function and Singular Value Entropy

YUAN Ye, HUANG Xiao-han, WANG Zhong-gang, ZHANG Yu

(Department of Civil Engineering, Logistic Engineering University, Chongqing 401311, China)

Based on the coherence function and effective frequency band of hammer excitation,the frequency band of the frequency response function was determined. Then the imaginary part of frequency-response function and the information entropy theory were combined to calculate the singular value entropy of Hankel-matrix which was constructured by the imaginary part of the frequency-response function. Finally， probabilistic method was used to calculate the abnormal condition and examine the internal defects of concrete. Based on these, the paper offers a simple and practical method of qualitative analysis to image the concrete defect. Firstly, the specimen was divided into two zones to differ in degree of defect.Then the projection function value of the defect cell was calculated and the non-defect cell was set to zero.Lastly, the inside situation of concrete was reacted by the contour of the projection function,and the test was performed to assess availability of the method.

Singular value entropy; Frequency response function; Coherence function; Defect

2016-09-10

袁 野(1992-),男，碩士研究生，主要從事軍事特種結(jié)構(gòu)檢測及加固。

袁 野, E-mail: 39523046@qq.com。

10.11973/wsjc201704004

TG115.28

A

1000-6656(2017)04-0017-05