不同材料儲(chǔ)罐底板中Lamb波頻散曲線(xiàn)的繪制與對(duì)比分析

程 莉，焦 陽(yáng)，李光海

(1.河北科技大學(xué)信息科學(xué)與工程學(xué)院，河北石家莊 050018；2.中國(guó)特種設(shè)備檢測(cè)研究院，北京 100013)

?

不同材料儲(chǔ)罐底板中Lamb波頻散曲線(xiàn)的繪制與對(duì)比分析

程 莉1，焦 陽(yáng)1，李光海2

(1.河北科技大學(xué)信息科學(xué)與工程學(xué)院，河北石家莊 050018；2.中國(guó)特種設(shè)備檢測(cè)研究院，北京 100013)

Lamb波頻散特性是進(jìn)行Lamb波無(wú)損檢測(cè)的前提，為了給儲(chǔ)罐底板的Lamb波檢測(cè)提供參考依據(jù)，根據(jù)Rayleigh-Lamb方程，采用二分法求解無(wú)限大板中Lamb波的頻散曲線(xiàn)，再通過(guò)Matlab軟件對(duì)求解過(guò)程進(jìn)行編程，繪制出Lamb波在不同儲(chǔ)罐底板材料中傳播的相速度和群速度曲線(xiàn)。為便于對(duì)比分析，繪制出各種材料S0模態(tài)相速度頻散特性曲線(xiàn)及A0模態(tài)相速度頻散特性曲線(xiàn)。通過(guò)對(duì)不同材料頻散特性的分析與比較，獲取對(duì)各種材料Lamb波檢測(cè)可選的模態(tài)為S0和A0模態(tài)，針對(duì)S0模態(tài)和A0模態(tài)分別給出了每種材料在實(shí)際檢測(cè)中可選的頻厚積范圍。

超聲學(xué)；Lamb波；頻散曲線(xiàn)；儲(chǔ)罐底板；無(wú)損檢測(cè)

程 莉，焦 陽(yáng)，李光海.不同材料儲(chǔ)罐底板中Lamb波頻散曲線(xiàn)的繪制與對(duì)比分析[J].河北工業(yè)科技,2016,33(5):385-391.

CHENG Li, JIAO Yang, LI Guanghai.Drawing and comparative analysis of Lamb wave dispersion curves in tank bottom plate of different materials [J].Hebei Journal of Industrial Science and Technology,2016,33(5):385-391.

隨著科學(xué)技術(shù)的發(fā)展，大型儲(chǔ)罐在工業(yè)生產(chǎn)和國(guó)家戰(zhàn)略?xún)?chǔ)備中的應(yīng)用日益廣泛且種類(lèi)復(fù)雜多樣。儲(chǔ)罐底板受所在地環(huán)境的影響容易產(chǎn)生腐蝕，導(dǎo)致結(jié)構(gòu)失效，所以人們對(duì)儲(chǔ)罐的安全性能越來(lái)越重視。為了更好地對(duì)儲(chǔ)罐底板這樣的大型板狀材料進(jìn)行方便、快捷的檢測(cè)，彌補(bǔ)傳統(tǒng)檢測(cè)技術(shù)的不足，超聲導(dǎo)波技術(shù)應(yīng)運(yùn)而生[1-9]。Lamb波是超聲無(wú)損檢測(cè)中最常見(jiàn)的一種導(dǎo)波形式，20世紀(jì)初由LAMB先生研究無(wú)限大板中正弦波問(wèn)題而得名。Lamb波檢測(cè)具有快速便捷的特點(diǎn)，非常適合板形結(jié)構(gòu)的大面積無(wú)損檢測(cè)，而Lamb波的頻散曲線(xiàn)就是為不同材料的板結(jié)構(gòu)選擇合適的模態(tài)和頻率進(jìn)行檢測(cè)的主要依據(jù)[10-13]。本文對(duì)超聲 Lamb 波的傳播及其頻散曲線(xiàn)的繪制做了詳盡的介紹，繪制出不同儲(chǔ)罐底板材料的頻散曲線(xiàn)，得到了超聲Lamb波在不同材料結(jié)構(gòu)中的傳播特性，并分別對(duì)不同材料儲(chǔ)罐底板的超聲Lamb波傳播特性和適宜檢測(cè)的單個(gè)模態(tài)特性進(jìn)行了分析，并得出結(jié)論。

1 Lamb波的傳播特性及Rayleigh-Lamb方程

Lamb波是在具有2個(gè)平行表面的結(jié)構(gòu)中傳播的，是由橫波和縱波相互耦合形成的一種應(yīng)力波。Lamb波傳播的實(shí)質(zhì)是平板中質(zhì)點(diǎn)的振動(dòng)，因其是板中的導(dǎo)波，又名為板波，根據(jù)平板上下表面質(zhì)點(diǎn)振動(dòng)方向不同可以把Lamb波分成對(duì)稱(chēng)S和反對(duì)稱(chēng)A兩種模態(tài)[10-18]，如圖1、圖2所示。

圖1 對(duì)稱(chēng)型Lamb波Fig.1 Symmetric Lamb wave

圖2 反對(duì)稱(chēng)型Lamb波Fig.2 Anti symmetric Lamb wave

每種模式又可以進(jìn)一步分成相速度不同的若干種模態(tài)，分別用S0,S1,S2,…,Sn和A0,A1,A2,…,Am表示。描述Lamb波頻散特性的是Rayleigh-Lamb方程[11-18]：

對(duì)稱(chēng)模態(tài)(S模態(tài))為

(1)

反對(duì)稱(chēng)模態(tài)(A模態(tài))為

(2)

其中：

(3)

式中：cL為縱波速度；cT為橫波速度；k為波數(shù)；f為L(zhǎng)amb波的頻率；cP為相速度；d為被測(cè)結(jié)構(gòu)厚度的1/2。

cL和cT可由式(4)、式(5)計(jì)算：

(4)

(5)

式中：E為被測(cè)結(jié)構(gòu)的楊氏模量；ρ和v分別表示結(jié)構(gòu)所用材料的密度和泊松比。

群速度cg可以通過(guò)相速度cP來(lái)求解，cg與cP有如下關(guān)系[11-12,15-17]：

(6)

(7)

因?yàn)樵谔囟ǖ牟牧舷?，其厚度d是一定的，所以式(7)可以表示為群速度cg與頻厚積fd的關(guān)系：

(8)

2 Lamb波頻散曲線(xiàn)繪制

由于對(duì)稱(chēng)模態(tài)與反對(duì)稱(chēng)模態(tài)的求解方式相同，所以下面以反對(duì)稱(chēng)模態(tài)為例闡述Lamb波頻散曲線(xiàn)的求解、Matlab編程及頻散曲線(xiàn)繪制[11]。

2.1 頻散曲線(xiàn)求解

由于Lamb波是多模態(tài)和頻散的，且板狀材料的選擇決定了橫波波速和縱波波速，在厚度一定時(shí)一個(gè)頻率下方程可能有多個(gè)解，即有多個(gè)模態(tài)的相速度cP，因此，將頻散方程轉(zhuǎn)變?yōu)檫m當(dāng)?shù)男问?，可以方便得到相速度cP與頻厚積fd的關(guān)系曲線(xiàn)，繼而求得群速度cg與頻厚積fd的關(guān)系曲線(xiàn)，這些曲線(xiàn)就是Lamb波的頻散曲線(xiàn)[16-17]。

以反對(duì)稱(chēng)模態(tài)為例，將式(3)代入式(2)中可得，

(9)

化簡(jiǎn)可得：

(10)

再約分提取相速度cP后，化簡(jiǎn)可得：

(11)

為了方便求相速度cP與頻厚積fd的關(guān)系，采用二分法計(jì)算，把式(11)轉(zhuǎn)化為各因式相乘的形式，令

(12)

使各個(gè)因式都為實(shí)數(shù)后，由y=0即可用二分法在Matlab里求Rayleigh-Lamb方程在實(shí)數(shù)域范圍內(nèi)的解，從而得到頻散曲線(xiàn)。

2.2Matlab編程與頻散曲線(xiàn)的繪制

使用迭代和二分的方法對(duì)Lamb方程進(jìn)行求解，可以先確定fd在所需范圍內(nèi)的掃描步長(zhǎng)，代入方程函數(shù)，然后在某一fd下掃描cP，分別代入cP1和cP2的值， 當(dāng)方程函數(shù)值y有為零的或者符號(hào)相反時(shí)，則有一合適的cP值對(duì)應(yīng)當(dāng)前的fd，為使方程結(jié)果滿(mǎn)足所需精度，可以在cP1和cP2之間進(jìn)行進(jìn)一步的二分計(jì)算[17-22]。

同時(shí)，在代入方程時(shí)應(yīng)當(dāng)考慮cP,cL和cT的關(guān)系，確定方程中根號(hào)下部分的虛實(shí)，從而保證得到的結(jié)果是實(shí)數(shù)解[16]。例如，在程序中為使方程簡(jiǎn)化，可以令

p=abs(sqrt((cp/ct)^2-1))；

q=abs(sqrt((cp/cl)^2-1))；

則cp<=ct時(shí)，p和q都為復(fù)數(shù)，又由-i*sin(ix)=sinhx和cos(ix)=coshx可將方程化為y=-4*p*q*sinh(pi*fd/cp*q)*cosh(pi*fd/cp*p)+(-p^2-1)^2*sinh(pi*fd/cp*p)*cosh(pi*fd/cp*q)。

同理，當(dāng)cp>ct且cp<=cl時(shí)，y=-4*p*q*sinh(pi*fd/cp*q)*cos(pi*fd/cp*p)+(p^2-1)^2*sin(pi*fd/cp*p)*cosh(pi*fd/cp*q)。

當(dāng)cp>cl時(shí)，y=4*p*q*sin(pi*fd/cp*q)*cos(pi*fd/cp*p)+(p^2-1)^2*sin(pi*fd/cp*p)*cos(pi*fd/cp*q)。

將上述式子在程序中寫(xiě)入1個(gè)子程序中，那么在主程序中也可方便調(diào)用。

計(jì)算完相速度cP與頻厚積fd的對(duì)應(yīng)關(guān)系后，畫(huà)出相速度的頻散曲線(xiàn)，之后群速度cg則可以直接調(diào)用已得數(shù)據(jù)進(jìn)行計(jì)算，從而畫(huà)出群速度曲線(xiàn)。

圖3即為整個(gè)頻散曲線(xiàn)繪制的流程圖。

根據(jù)上述求解方法對(duì)Lamb波在儲(chǔ)罐底板常用的a3鋼板中傳播的頻散曲線(xiàn)進(jìn)行了計(jì)算，將a3鋼板的材料參數(shù)E=210GPa，v=0.3，ρ=7.85g/cm3輸入程序中，繪制得到相應(yīng)的相速度頻散曲線(xiàn)和群速度頻散曲線(xiàn)如圖4所示[20-22]。

由圖4可得Lamb波頻散曲線(xiàn)的特點(diǎn)[11]如下：

1) 在低頻情況下，模態(tài)較為簡(jiǎn)單，只有S0和A0 2種模態(tài)，而頻率越高模態(tài)種類(lèi)越復(fù)雜，會(huì)同時(shí)產(chǎn)生2種以上的模態(tài)；

2) 除了S0和A0模態(tài)，其他模態(tài)都有截止頻率，且截止頻率以下不存在該模態(tài)；

3) 各個(gè)模態(tài)的群速度不固定，而是隨頻率增加而上下浮動(dòng)并逐漸趨向某一固定值。

圖3 Matlab頻散曲線(xiàn)繪制流程圖 Fig.3 Flow chart of drawing frequency dispersion curve on Matlab

3 不同材料儲(chǔ)罐底板的頻散特性分析

為了獲取Lamb波對(duì)不同儲(chǔ)罐底板材料檢測(cè)可選的模態(tài)，選取3種常用的儲(chǔ)罐底板材料，即鋼板、鋁板及玻璃鋼板，分別繪制出頻散曲線(xiàn)并進(jìn)行分析比較。算例中已給出鋼板的頻散曲線(xiàn)，利用上述的計(jì)算方法所繪制的另外2種材料的頻散曲線(xiàn)如圖5、圖6所示。

由圖4—圖6可得不同材料的Lamb波頻散曲線(xiàn)的特點(diǎn)：

1) 不同材料的儲(chǔ)罐底板的頻散曲線(xiàn)大概分布和復(fù)雜程度類(lèi)似，且模態(tài)分布簡(jiǎn)單的情況都在低頻厚積時(shí)，為使檢測(cè)中盡可能地克服多模態(tài)的問(wèn)題，最好選用A0和S0模態(tài)進(jìn)行檢測(cè)；

2) 不同材料的群速度隨頻厚積的變化而變化的過(guò)程相似，而其大小卻不太相同，鋼板和鋁板群速度總體大小接近，鋼板的群速度在某些模態(tài)下要比鋁板高，而玻璃鋼的群速度最小。

分析得到，為了方便檢測(cè)，選擇頻率和模態(tài)時(shí)應(yīng)該盡量選擇低頻厚積，適合的模態(tài)有S0和A0[20-23]，接下來(lái)分別對(duì)不同材料下這2個(gè)模態(tài)的相速度進(jìn)行比較。

3.1 不同材料S0模態(tài)相速度頻散特性分析

將3種材料的S0模態(tài)畫(huà)在一張坐標(biāo)圖上，如圖7所示，3種材料的S0模態(tài)走勢(shì)相同，都是在低頻厚積下有稍高的相速度，隨著頻厚積的加大快速轉(zhuǎn)折到低相速度的位置，并逐漸趨向固定值。

對(duì)比后發(fā)現(xiàn)，雖然3種材料的曲線(xiàn)走勢(shì)類(lèi)似，但是某一頻厚積下的相速度大小不同，鋼板相速度最大，鋁板次之，玻璃鋼相速度最小。另外轉(zhuǎn)折點(diǎn)位置不同，玻璃板轉(zhuǎn)折處的頻厚積最低，之后是鋁板，鋼板的轉(zhuǎn)折處頻厚積最大。

圖4 Matlab計(jì)算并繪制a3鋼板的頻散曲線(xiàn) Fig.4 Frequency dispersion curve of a3 steel plate calculate and drawn by matlab

圖5 鋁板的頻散曲線(xiàn)(E=70 GPa，v=0.33，ρ=2.7 g/cm3) Fig.5 Frequency dispersion curve of aluminum plate (E=70 GPa, v=0.33, ρ=2.7 g/cm3)

圖6 玻璃鋼板的頻散曲線(xiàn)(E=21 GPa，v=0.3，ρ=1.85 g/cm3)Fig.6 Frequency dispersion curve of glass fibre reinforced plastic plate (E=21 GPa, v=0.3, ρ=1.85 g/cm3)

圖7 3種材料S0模態(tài)相速度頻散曲線(xiàn) Fig.7 Phase velocity curves of S0 modes of three materials

由此可知，不同材料適合檢測(cè)頻厚積的范圍不同，鋼板范圍最寬，為0～2 MHz·mm;鋁板次之，為0～1.95 MHz·mm;而玻璃鋼可檢測(cè)范圍最窄，為0～1.1 MHz·mm。

3.2 不同材料A0模態(tài)相速度頻散特性分析

將3種材料的A0模態(tài)畫(huà)在一張坐標(biāo)圖上，如圖8所示，3種材料的A0模態(tài)走勢(shì)相同，都是在低頻厚積下相速度較低，隨著頻厚積的加大而相速度逐漸上升，并慢慢趨向1個(gè)固定值。

圖8 3種材料A0模態(tài)相速度頻散曲線(xiàn) Fig.8 Phase velocity curves of A0 modes of three materials

對(duì)比后發(fā)現(xiàn)，雖然3種材料的曲線(xiàn)走勢(shì)類(lèi)似，但是和對(duì)稱(chēng)模態(tài)相同的是某一頻厚積下的相速度大小不同，鋼板相速度最大，鋁板次之，玻璃鋼相速度最小。另外相速度上升的位置不同，玻璃鋼上升處的頻厚積最低，之后是鋁板，鋼板的上升處頻厚積最大。由此可知，不同材料適合檢測(cè)頻厚積的范圍不同，鋼板范圍最寬，為0.8～2 MHz·mm；鋁板次之，為0.8～1.95 MHz·mm；而玻璃鋼可檢測(cè)范圍最窄，為0.7～1.1 MHz·mm。

總得來(lái)看，在不同材料適合檢測(cè)的頻厚積范圍下，S0模態(tài)在稍微低頻處的相速度相對(duì)較大，而A0模態(tài)在合適的范圍下，則是在相對(duì)較高的頻厚積處相速度稍大，由此可以進(jìn)一步縮短范圍，便于選取合適頻率的激勵(lì)信號(hào)。但由于不同材料的頻散曲線(xiàn)不同，可選范圍也不相同，還需在檢測(cè)前對(duì)各個(gè)材料進(jìn)行進(jìn)一步的分析并選取合適頻率合適模態(tài)的激勵(lì)源。

4 結(jié) 論

本文根據(jù)Rayleigh-Lamb方程，采用二分法并通過(guò)Matlab編程計(jì)算獲取儲(chǔ)罐底板常用3種材料的頻散特性曲線(xiàn)。對(duì)3種材料的頻散特性曲線(xiàn)進(jìn)行整體分析，并對(duì)S0模態(tài)下和A0模態(tài)下3種材料的頻散特性進(jìn)行對(duì)比分析，得出以下結(jié)論：

1)3種材料的相速度和群速度隨頻厚積的變化趨勢(shì)相似，但對(duì)應(yīng)同一頻厚積其大小不同。鋼板與鋁板接近，玻璃鋼最小。但都在低頻厚積時(shí)模態(tài)分布簡(jiǎn)單，為使檢測(cè)中盡可能地克服多模態(tài)的問(wèn)題，應(yīng)選用A0和S0模態(tài)。

2)對(duì)于S0模態(tài)，3種材料適合檢測(cè)頻厚積的范圍不同，鋼板范圍最寬，為0～2 MHz·mm；鋁板次之，為0～1.95 MHz·mm；而玻璃鋼可檢測(cè)范圍最窄，為0～1.1 MHz·mm。

3) 對(duì)于A0模態(tài)，3種材料適合檢測(cè)頻厚積的范圍不同，鋼板范圍最寬，為0.8～2 MHz·mm；鋁板次之，為0.8～1.95 MHz·mm；而玻璃鋼可檢測(cè)范圍最窄，為0.7～1.1 MHz·mm。

[1] 李楊. 儲(chǔ)罐底板超聲導(dǎo)波檢測(cè)傳感器研制及應(yīng)用研究[D]. 北京: 北京工業(yè)大學(xué), 2012.

LI Yang. Transducers Development and Application for Ultrasound Guided Wave Inspection on Oil Tank Bottom Plate[D]. Beijing: Beijing University of Technology, 2012.

[2] 王勇，沈功田，閆河，等. 壓力容器無(wú)損檢測(cè)——大型常壓儲(chǔ)罐的無(wú)損檢測(cè)技術(shù)[J]. 無(wú)損檢測(cè), 2005,27(9): 487-490.

WANG Yong, SHEN Gongtian, YAN He, et al. Nondestructive testing of pressure vessels: Nondestructive testing technique for large atmospheric storage tanks[J]. Nondestructive Testing, 2005,27(9): 487-490.

[3] 孫曉云，楊陽(yáng)，劉東輝，等. 基于ARM & Windows CE的錨桿無(wú)損檢測(cè)系統(tǒng)軟件設(shè)計(jì)[J].河北科技大學(xué)學(xué)報(bào)，2012，33(3)：244-247.

[4] 劉東輝，米博，趙東松，等. 基于HHT的錨桿錨固質(zhì)量聲波無(wú)損檢測(cè)[J].河北科技大學(xué)學(xué)報(bào)，2013,34(4)：208-212.

[5] 帥健，許學(xué)瑞，韓克江. 原油儲(chǔ)罐檢修周期[J]. 石油學(xué)報(bào), 2012,33(1): 157-163.

SHUAI Jian, XU Xuerui, HAN Kejiang. Overhaul interval of crude oil tanks[J].Acta Petrolei Sinica, 2012,33(1): 157-163.

[6] 楊志軍. 儲(chǔ)罐底板磁檢測(cè)技術(shù)研究[D]. 大慶: 大慶石油學(xué)院, 2003.

YANG Zhijun. The Study of Magnetic Flux Leakage Technology for the Inspection of Tank Floors[D]. Daqing: Northeast Petroleum University, 2003.

[7] 劉鎮(zhèn)清，劉曉. 超聲無(wú)損檢測(cè)的若干新進(jìn)展[J]. 無(wú)損檢測(cè), 2000,22(9): 403-405.

LIU Zhenqing, LIU Xiao. Development of ultrasonic nondestructive testing techniques[J]. Nondestructive Testing, 2000,22(9): 403-405.

[8] 張清華. 基于超聲和渦流的復(fù)合式無(wú)損檢測(cè)技術(shù)[D]. 廣州: 華南理工大學(xué), 2010.

ZHANG Qinghua. Technique of Complex NDT based on Ultrasonic & Eddy Current Testing[D]. Guangzhou:South China University of Technology, 2010.

[9] 宋志東. 超聲導(dǎo)波技術(shù)在管道缺陷檢測(cè)中的研究[D]. 天津: 天津大學(xué), 2006.

SONG Zhidong. Study on Pipeline Defects Detection by Ultrasonic Guided Waves[D]. Tianjin: Tianjin University, 2006.

[10] MANSFIELD T L. Lamb wave inspection of aluminum sheet[J]. Materials Evaluation,1975, 33(4):96-100.

[11] 閻石，張海鳳，蒙彥宇. Lamb波頻散曲線(xiàn)的數(shù)值計(jì)算及試驗(yàn)驗(yàn)證[J]. 華中科技大學(xué)學(xué)報(bào)(城市科學(xué)版), 2010,27(1): 1-4.

YAN Shi, ZHANG Haifeng, MENG Yanyu. Numerical calculation and experimental validation for Lamb wave dispersion curves[J]. Journal of HUST(Urban Science Edition), 2010,27(1): 1-4.

[12]徐可北. 金屬薄板蘭姆波檢驗(yàn)技術(shù)[J]. 無(wú)損檢測(cè), 1999,21(10):461-465.

XU Kebei. Nondestructive testing of metallic plates by Lamb waves[J]. Nondestructive Testing, 1999, 21(10):461-465.

[13]黃志強(qiáng)，王明波，王維斌，等. 基于超聲Lamb波技術(shù)的儲(chǔ)罐底板缺陷檢測(cè)[J]. 科技導(dǎo)報(bào), 2013, 31(16): 47-50.

HUANG Zhiqiang, WANG Mingbo, WANG Weibin, et al. Defect detection for the bottom plate of oil tank based on ultrasonic Lamb waves[J]. Science & Technology Review, 2013, 31(16): 47-50.

[14]陳強(qiáng). 基于Lamb波的薄板鋼結(jié)構(gòu)損傷識(shí)別研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2013.

CHEN Qiang. Damage Identification Research for Steel Plate Based on Lamb Wave[D]. Harbin: Harbin Engineering University, 2013.

[15]李博成. Lamb波結(jié)構(gòu)損傷識(shí)別中信號(hào)處理與采集的研究與實(shí)驗(yàn)[D]. 成都: 電子科技大學(xué), 2013.

LI Bocheng. Research and Experiment on Signal Processing and Acquisition of Structural Damage Identification based on Lamb Wave[D]. Chengdu: University of Electronic Science and Technology of China, 2013.

[16]陳亮，王恩報(bào)，馬豪，等. 基于有限元方法的薄板超聲Lamb波頻散曲線(xiàn)計(jì)算[J]. 實(shí)驗(yàn)室研究與探索, 2014,33(11): 28-31.

CHEN Liang, WANG Enbao, MA Hao, et al. Calculation of ultrasonic Lamb wave dispersion curve of the thin plate by the finite element method[J]. Research and Exploration in Laboratory, 2014,33(11): 28-31.

[17]許西寧，余祖俊，朱力強(qiáng). 圖解法求解Lamb波頻散方程[J]. 電子測(cè)量與儀器學(xué)報(bào), 2012,26(11): 966-971.

XU Xining, YU Zujun, ZHU Liqiang. A graphical method to solve a dispersion equation of Lamb wave[J]. Journal of Electronic Measurement and Instrument, 2012,26(11): 966-971.

[18]曹正敏，林莉，李喜孟. 蘭姆波頻散曲線(xiàn)的繪制與試驗(yàn)驗(yàn)證[J]. 理化檢驗(yàn)(物理分冊(cè)), 2008,44(9): 482-500.

CAO Zhengmin, LIN Li, LI Ximeng. Plotting and experimental validation of Lamb wave dispersion curve[J]. Physical Testing and Chemical Analysis. Part: A Physical Testing, 2008,44(9): 482-500.

[19]劉增華，何存富，吳斌. 單層板頻散曲線(xiàn)的矩陣表示及MATLAB實(shí)現(xiàn)[J]. 無(wú)損檢測(cè), 2005,27(5): 225-227.

LIU Zenghua, HE Cunfu, WU Bin. Matrix expression for dispersion curve of single plate and its implementation by MATLAB[J]. Nondestructive Testing, 2005,27(5): 225-227.

[20]李迎. 復(fù)合材料Lamb波檢測(cè)中S0單模式信號(hào)的提取[J]. 無(wú)損檢測(cè), 2012,34(8): 30-33.

LI Ying. Extracting the S0 single mode of Lamb wave in composite material[J]. Nondestructive Testing, 2012,34(8): 30-33.

[21]陳軍，李志浩，林莉，等. 鋁板中Lamb波檢測(cè)的實(shí)驗(yàn)研究[J]. 應(yīng)用聲學(xué), 2011,30(2): 98-104.

CHEN Jun, LI Zhihao, LIN Li, et al. Experiment investigations of Lamb waves in an aluminum plate[J]. Applied Acoustics, 2011,30(2): 98-104.

[22]嚴(yán)剛，周麗. 基于Lamb波的復(fù)合材料結(jié)構(gòu)損傷成像研究[J]. 儀器儀表學(xué)報(bào), 2007,28(4): 583-589.

YAN Gang, ZHOU Li. Study on damage imaging for composite structure using Lamb wave[J]. Chinese Journal of Scientific Instrument, 2007,28(4): 583-589.

[23]禹化民，王維斌，呂小青，等. 儲(chǔ)罐底板缺陷蘭姆波檢測(cè)的影響因素[J]. 無(wú)損檢測(cè), 2015,37(7): 19-24.

YU Huamin, WANG Weibin, LYU Xiaoqing, et al. Influencing factors for the Lamb wave detection of tank bottom defects[J]. Nondestructive Testing, 2015,37(7): 19-24.

Drawing and comparative analysis of Lamb wave dispersion curves in tank bottom plate of different materials

CHENG Li1, JIAO Yang1, LI Guanghai2

(1. School of Information Science and Engineering, Hebei University of Science and Technology, Shijiazhuang, Hebei 050018, China; 2. China Special Equipment Inspection and Research Institute, Beijing 100013, China)

Lamb wave dispersion characteristic is the precondition of lamb wave nondestructive detection. The article provides an important reference for the Lamb wave detection of the tank bottom plate. Based on the Rayleigh lamb equation, Lamb wave dispersion curves in the infinite plate are obtained according to dichotomy in this paper. Through programming the Matlab software for the solving process, phase velocity curves and group velocity curves that Lamb wave propagate in different tank bottom plat materials are drawn. In order to facilitate the comparison and analysis, the S0 mode phase velocity dispersion characteristic curve and the A0 mode phase velocity dispersion characteristic curve of various materials are also plotted. Through the analysis and comparison of the dispersion characteristics of different materials, the optional mode for Lamb wave detection for different materials are S0 mode and A0 mode. For S0 mode and A0 mode, the optional frequency-thickness product ranges for each material in the actual detection are given, respectively.

ultrasonics; Lamb wave; frequency dispersion curve; tank bottom plate; nondestructive testing

1008-1534(2016)05-0385-07

2016-04-27;

2016-06-17；責(zé)任編輯：李 穆

國(guó)家科技支撐計(jì)劃項(xiàng)目(211BAK06B03)

程 莉(1991—)，女，河南周口人，碩士研究生，主要從事超聲檢測(cè)和傳感器方面的研究。

焦 陽(yáng)教授。E-mail:jyang939@163.com

TK172

A

10.7535/hbgykj.2016yx05005