吳榮興，于蘭珍，王海林，李建中

（1.寧波職業(yè)技術學院建工學院，浙江 寧波 315800；2.寧波大學機械與力學學院，浙江 寧波 315211）

利用瑞利波檢測層狀混凝土結構的材料特性

吳榮興1，2，于蘭珍1，2，王海林1，李建中1

（1.寧波職業(yè)技術學院建工學院，浙江 寧波 315800；2.寧波大學機械與力學學院，浙江 寧波 315211）

為獲得層狀混凝土結構的材料特性，利用瑞利波在半無限大混凝土和層狀混凝土中的傳播規(guī)律，分別建立兩種模型中瑞利波傳播的波速方程。基于獲得的層狀混凝土結構中瑞利波的波速方程，提出一種從實測波速來推定混凝土材料特性的反演方法。通過數(shù)值計算獲得兩種模型中波速和混凝土彈性模量的關系曲線，理論計算結果與實驗結果基本一致。研究結果表明：由于混凝土層的相互作用，結構中存在多個瑞利波振動模態(tài)以及相應的波速，同時波速和上層混凝土彈性模量之間存在對應關系。建立的研究方法和計算結果可直接應用于實際混凝土的瑞利波檢測。

瑞利波；混凝土；彈性模量；檢測；波速

0　引　言

超聲波檢測作為一種無損檢測方法，在各種混凝土結構材料性能檢測中得到了廣泛的應用[1]。但是由于混凝土內(nèi)部存在各種孔隙和裂縫以及混凝土本身粘滯性的影響，超聲脈沖波的衰減較大，因此探測深度非常有限[2-4]。同時由于存在層狀混凝土和新老混凝土的粘結使混凝土材料特性沿深度發(fā)生變化，超聲脈沖波和縱波在層狀結構之間來回反射而降低檢測精度甚至無法獲得其有效信號[5-6]。因此利用各種新波型超聲波對混凝土結構進行材料性能檢測在國內(nèi)外都得到了許多有益的嘗試[7-11]。

通過對地震波的觀察，瑞利（Rayleigh）從理論上驗證了聲表面波的存在。隨著聲表面波技術的日趨成熟，其在傳感器和無損檢測領域都得到了廣泛應用[7]。夏唐代等[8]研究了瑞利波在防滲墻中的傳播特性。馮小娟等[9]分析了飽和土中樁對瑞利波的動態(tài)響應。何旭升等[10]利用表面波并結合回彈儀來測定混凝土的質(zhì)量。利用瑞利波對混凝土結構進行無損檢測不僅能夠彌補瞬態(tài)縱波和脈沖波的不足，而且在探測深度和探測精度都有大幅提高[12-13]。在實際檢測中，各種混凝土結構大多是層狀結構。為了更加精確地檢測混凝土結構的材料特性，本文建立了瑞利波在多層混凝土結構中的波動方程，通過對波動方程的求解來推定混凝土的材料特性。

1　半無限大混凝土中的瑞利波

瑞利波在層狀混凝土結構中傳播的基本構造如圖1所示，主要由層狀混凝土結構、瑞利波發(fā)射和接收換能器等構成。從發(fā)射端發(fā)出的瑞利波通過層狀混凝土層后被接收端獲取信號，利用獲得的瑞利波信息來推定層狀混凝土的材料特性[12-13]。當上層混凝土和基體混凝土擁有一樣的材料特性時，認為瑞利波在半無限大混凝土中傳播。研究表明瑞利波的振動隨著厚度方向的增加而呈指數(shù)衰減[7]，因此假設位移的表達式為

圖1　層狀混凝土結構示意圖

式中u1（u2）、A（B）、k、β、x1（x2）、c和t分別是瑞利波的位移、振幅、波數(shù)、衰減系數(shù)、坐標、波速和時間。

這里認為半無限大混凝土結構為材質(zhì)均勻和各向同性的線彈性材料[5]，由式（1）可以得到應力為

式中λ和μ為半無限大混凝土的拉梅常數(shù)。拉梅常數(shù)與彈性模量以及泊松比的換算關系如下：

式中E和ν分別為半無限大混凝土的彈性模量和泊松比。各向同性材料的波動方程為

式中：ρ——混凝土的密度；

u咬1，u咬2——位移對時間求兩次偏導。

將式（1）和式（2）代入式（4），可以得到：

式中c2L=（λ+2μ）/ρ和c2L=μ/ρ分別為縱波和橫波波速[7]。通過對式（5）求解可以得到兩個衰減系數(shù)β為

相應的振幅比為

如圖1所示，考慮瑞利波沿x1方向傳播，那么式（5）的解為

式中A1和A2為重新設定的瑞利波振幅?？梢垣@得其應力表達式為

瑞利波在如圖1所示的半無限大混凝土中傳播時，必須滿足自由邊界條件如下：

將式（9）代入式（10），可以獲得關于振幅A1和A2的方程如下：

如果待定振幅A1和A2存在非零解，其系數(shù)行列式的值必須為零，可獲得瑞利波的波速方程如下：

在給定半無限大混凝土的彈性模量、密度和泊松比的情況下，可以對式（12）進行數(shù)值求解獲得瑞利波波速。同樣在實際檢測中可以利用檢測到的瑞利波波速來反推混凝土的材料特性[5]。為了與實驗結果進行比較，這里取ρ=2 400 kg/m3，ν=0.2。獲得的瑞利波波速和混凝土彈性模量的關系曲線如圖2所示。

圖2　半無限大混凝土中瑞利波波速和彈性模量的關系曲線

從圖中可以發(fā)現(xiàn)，在半無限大混凝土中傳播的瑞利波波速與混凝土的彈性模量存在近似二次函數(shù)關系。1、2和3點分別是半無限大混凝土C20、C25和C30混凝土彈性模量的實測值[3]。圖2表明獲得的理論曲線和實驗結果基本一致，理論曲線比較準確地反映了混凝土彈性模量變化的規(guī)律。這也為利用瑞利波檢測層狀混凝土的彈性模量奠定了基礎。

2　層狀混凝土結構中的瑞利波

當瑞利波在層狀混凝土結構中傳播時，可以認為在基體混凝土中的傳播如前面所述。而在上層混凝土層中傳播時，位移模式[7]應該修正為

同樣求出上層混凝土中的應力表達式為

式中h為上層混凝土的厚度。

將上層混凝土的位移和應力表達式以及基體混凝土的位移和應力表達式代入式（15），可以得到關于振幅A1，A2，A11，A12，A21和A22的6個方程組。若使這些振幅存在非零解，必須令這些振幅的系數(shù)行列式值為零，這樣就建立了瑞利波在層狀混凝土中傳播的波速方程。

為了運算和推導方便，對上層混凝土厚度作歸一化處理：

式中ξ是波長。

首先分析上層混凝土對瑞利波波速的影響。給定基體混凝土的材料參數(shù)E=47 GPa，ν=0.2，ρ= 2400kg/m3和上層混凝土的參數(shù)獲得的層狀混凝土中瑞利波波速與上層混凝土厚度之間的關系如圖3所示。

由于存在多層結構，可以觀察到多個瑞利波波速存在[7]。當上層混凝土厚度逐漸變小時，計算得到的瑞利波波速中最小的那個波速就是半無限大混凝土中求得的波速（2 643 m/s），這樣就驗證了推導的正確性。同時發(fā)現(xiàn)隨著上層混凝土厚度的不斷增加，所求得的瑞利波波速有減小的趨勢并最終趨向于上層混凝土的瑞利波波速（2104m/s）。這是因為隨著上層混凝土厚度增加，瑞利波作為表面波其位移和振動能量主要集中于上層混凝土，因此所獲得的瑞利波波速主要受其材料特性的影響[13]。瑞利波波速與上層混凝土厚度關系曲線還提供了一個檢測擁有不同材料特性的各混凝土層厚度的方法。

圖3　層狀混凝土中瑞利波波速與上層混凝土厚度之間的關系

在給定上層混凝土厚度的條件下，獲得了層狀混凝土中瑞利波波速與上層混凝土彈性模量的關系如圖4所示。同樣可以觀察到瑞利波的多階振動模態(tài)，在實際檢測中常用的是第1階振動模態(tài)[11，14]。實驗獲得的上層混凝土3種型號的彈性模量都略小于理論計算值，但是趨勢完全一致，主要原因是上層混凝土厚度為一個波長的情況下瑞利波主要受基體混凝土的影響[11，14]。圖4表明隨著上層混凝土彈性模量的不斷增加，所獲得的聲表面波波速都在增加，這符合實驗和前期理論的研究結果[13]。在實際檢測中，只要檢測到瑞利波波速并量取上層混凝土厚度就可以反推上層混凝土的彈性模量。

圖4　層狀混凝土中瑞利波波速與上層混凝土彈性模量的關系（上層混凝土厚度H=1ξ）

3　結束語

本文研究了瑞利波在層狀混凝土結構中的傳播特性，獲得了瑞利波波速和混凝土彈性模量之間的關系曲線。當瑞利波在半無限大混凝土中傳播時，所獲得的瑞利波波速可以推定混凝土的彈性模量，這些結果與實驗結果能良好吻合。當瑞利波在層狀混凝土中傳播時，隨著上層混凝土厚度的增加，所獲得的瑞利波波速有減小的趨勢，并且這些趨勢符合瑞利波的傳播特征。本文還分析了在給定上層混凝土厚度條件下，上層混凝土的彈性模量對瑞利波波速的影響，所獲得的曲線提供了檢測應用的方法。

[1]EVERYA G.Measurement of the near-surface elastic properties of solids and thin supported films[J].Science Technology，2002，13（5）：21-39.

[2]GOUEYGOU M，PIWAKOWSKI B，NAFFA S O，et al. Assessment of broad-band ultrasonic attenuation measurements in inhomogeneous media[J].Ultrasonics，2002，40（1-8）：77-82.

[3]趙治廣.超聲波檢測非金屬材料的研究與開發(fā)[J].實驗技術與管理，2006，23（10）：30-32.

[4]SONG W J.Nondestructive evaluation techniques using impact-generated stress waves in concrete[D].Illionis: Northwestern University，2001.

[5]張建仁，黃浩，彭輝，等.用縱波換能器測量混凝土橫波波速及動彈性模量的實驗研究 [J].中外公路，2011，31（3）：89-92.

[6]劉衛(wèi)東，張東芹，劉勝利，等.超聲波檢測混凝土結構抗凍性的方法研究[J].無損檢測，2002，24（5）：196-199.

[7]吳榮興，袁麗莉.溫度效應對聲表面波傳播的影響[J].壓電與聲光，2013，35（1）：10-15.

[8]夏唐代，石中明，楊洛.表面波法在防滲墻質(zhì)量檢測中的應用研究[J].三峽大學學報（自然科學版），2003，25（1）：56-58.

[9]馮小娟，黃義，吳炳軍.飽和土中樁對瑞利波的動力響應[J].地震研究，2007，30（1）：64-71.

[10]何旭升，張雷，李守輝，等.應用回彈儀與聲表面波儀聯(lián)合檢測隧道混凝土的質(zhì)量[J].華北水利水電學院學報，2002，23（2）：19-23.

[11]ROSEJL.Ultrasonicwavesinsolid media[D]. Cambridge：Cambridge University，1999.

[12]楊天春，何繼善，魯光銀，等.道路結構型地層瑞利波相速度頻散曲線的完整求取[J].中南大學學報（自然科學版），2013，44（2）：642-647.

[13]曲華，寧建國，李學慧.利用瑞利波頻譜分析法無損檢測混凝土結構[J].混凝土，2006（2）：76-80.

[14]AULDBA.Acousticfields and wave in solids[M]. Florida：Krieger Publishing Company，1995：112-119.

（編輯：李妮）

Detection of material properties of layered concrete structures using Rayleigh wave

WU Rongxing1，2，YU Lanzhen1，2，WANG Hailin1，LI Jianzhong1
（1.Department of Architectural Engineering，Ningbo Polytechnic，Ningbo 315800，China；2.School of Mechanical Engineering and Mechanics，Ningbo University，Ningbo 315211，China）

In order to obtain the material properties of layered concrete structures，velocity equations for Rayleigh wave of two models are established according to the propagation law of Rayleigh waves in semi-infinite concrete structures and layered concrete structures.Based on the velocity equations obtained in layered concrete structures，an inverse method is proposed to calculate the material properties with the wave velocities measured.The relation curves between wave velocities and the elastic modulus of concrete in the two models are obtained via numerical computation.And the results are similar to the experimental results.The research results show that there are many vibration modes and corresponding velocities in layered concrete structures because of concrete layer interaction，as well as a corresponding relation between wave velocities and the elastic moduli of upper concrete layers.The established research method and its computation results can be directly applied to detect the surface acoustic waves of real concrete.

Rayleigh wave；concrete；elastic modulus；detection；wave velocity

A

1674-5124（2016）03-0032-04

10.11857/j.issn.1674-5124.2016.03.008

2015-09-04；

2015-10-28

國家自然科學基金項目（11072116）；浙江省教育廳高校訪問工程師校企合作項目（FG2014041）；寧波職業(yè)技術學院青年博士創(chuàng)新項目和科研項目（2013001，NZ14001）

吳榮興（1982-），男，浙江東陽市人，講師，博士，主要從事壓電聲波檢測理論和分析。