非連續(xù)阻抗粘接結(jié)構(gòu)脫粘缺陷的稀布陣列超聲成像*

張輝 朱文發(fā)? 范國(guó)鵬 張海燕

1) (上海工程技術(shù)大學(xué)城市軌道交通學(xué)院,上海 201620)

2) (上海大學(xué)通信與信息工程學(xué)院,上海 200444)

無砟軌道是典型的非連續(xù)阻抗粘接結(jié)構(gòu),在重載、環(huán)境惡化等因素的影響下,脫粘現(xiàn)象頻發(fā),嚴(yán)重危及列車運(yùn)行安全.本文提出一種高精度陣列超聲快速成像方法,建立考慮層間聲速差異的非連續(xù)阻抗粘接結(jié)構(gòu)聲速理論模型,采用射線追蹤方法獲取聲波在介質(zhì)中的傳播路徑和時(shí)間.基于實(shí)數(shù)編碼設(shè)計(jì)高自由度稀布陣列,構(gòu)建稀布陣列合成孔徑聚焦技術(shù)(synthetic aperture focusing technique,SAFT)成像,提高檢測(cè)效率.在無砟軌道結(jié)構(gòu)上的實(shí)驗(yàn)結(jié)果表明: 射線追蹤方法能夠準(zhǔn)確計(jì)算超聲波的傳播路徑和傳播時(shí)間,提高檢測(cè)精度;優(yōu)化設(shè)計(jì)的稀布陣列方向圖主瓣寬度窄,旁瓣增益低,可提高檢測(cè)效率和聲場(chǎng)指向性;脫粘缺陷成像誤差在±5%以內(nèi)時(shí),稀布陣列SAFT 成像方法效率提高了30.9%,可為該類缺陷檢測(cè)提供理論支撐.

1 引言

多層復(fù)合粘接結(jié)構(gòu)能夠均勻化應(yīng)力分布,避免應(yīng)力集中問題,廣泛應(yīng)用于航空航天、汽車、軌道交通等工業(yè)領(lǐng)域.如高鐵基礎(chǔ)結(jié)構(gòu)中大量使用的無砟軌道結(jié)構(gòu)就是一種由軌道板、砂漿層、底座板和基床組成的典型非連續(xù)阻抗多層混凝土復(fù)合粘接結(jié)構(gòu),各層間的材料、厚度、密度均不相同.作為列車運(yùn)行的主要承載結(jié)構(gòu),在高速列車不斷沖擊以及外部環(huán)境的長(zhǎng)期作用下,脫粘現(xiàn)象頻發(fā)[1].脫粘缺陷降低了無砟軌道承載能力,危及列車運(yùn)行安全,嚴(yán)重時(shí)將導(dǎo)致列車傾覆等重大安全事故[2?4].目前,針對(duì)脫粘缺陷常用的現(xiàn)代無損檢測(cè)方法主要有沖擊回波法、紅外熱成像法、探地雷達(dá)法、聲發(fā)射法、超聲波等.沖擊回波法[5,6]只能定性分析被檢測(cè)物體的內(nèi)部是否存在缺陷,難以定量分析.紅外熱成像法[7,8]通過被檢試樣表面紅外輻射性能的改變實(shí)現(xiàn)檢測(cè),難以實(shí)現(xiàn)多層復(fù)合粘接結(jié)構(gòu)脫粘缺陷的檢測(cè).探地雷達(dá)法[9,10]的電磁波在遇到鋼筋時(shí)無法繼續(xù)傳播,該方法不適用于分布有大量鋼筋的無砟軌道脫粘缺陷檢測(cè).聲發(fā)射法[11]受到環(huán)境噪聲影響較大,檢測(cè)準(zhǔn)確度較低,并且該方法不可逆,信號(hào)無法重復(fù)獲得.

超聲波具有傳播能量大,指向性好等優(yōu)勢(shì),在對(duì)鋼筋混凝土結(jié)構(gòu)檢測(cè)時(shí),能夠穿過內(nèi)部鋼筋不受其影響,越來越多的學(xué)者使用超聲波檢測(cè)混凝土結(jié)構(gòu)內(nèi)部損傷.尤其是利用線性周期超聲陣列通過全矩陣捕獲方式獲取檢測(cè)信號(hào),再對(duì)全矩陣數(shù)據(jù)進(jìn)行后處理和虛擬聚焦,就能夠?qū)崿F(xiàn)缺陷成像,精確地可視化分析缺陷的位置和形狀,如全聚焦成像方法.然而,全聚焦成像方法建立在單一均勻介質(zhì)假設(shè)基礎(chǔ)上的超聲成像理論,難以直接運(yùn)用于非連續(xù)阻抗結(jié)構(gòu)中脫粘缺陷的成像.無砟軌道結(jié)構(gòu)中不同層之間聲阻抗不同,超聲波在每一層傳播的速度都不同[12].超聲波在無砟軌道結(jié)構(gòu)中傳播路徑更加復(fù)雜,傳播路徑和傳播時(shí)間難以通過常規(guī)的計(jì)算方法獲得.射線追蹤可以解決聲波在介質(zhì)中不同層之間的傳播時(shí)延問題,但常規(guī)的射線軌跡追蹤法[13?15]是基于反向映射原理,在計(jì)算聲波傳播時(shí)間的過程中均存在反復(fù)迭代計(jì)算,成像效率低.

此外,利用全矩陣數(shù)據(jù)進(jìn)行超聲成像計(jì)算耗時(shí)多,而無砟軌道結(jié)構(gòu)的檢修時(shí)間一般在凌晨0—4 時(shí)進(jìn)行,檢測(cè)時(shí)間短暫.為了提高超聲成像效率,很多學(xué)者在滿足時(shí)效性要求方面進(jìn)行了研究.一方面是通過提高計(jì)算硬件的性能,比如,圖形處理器(graphics processing unit,GPU)加速技術(shù)[16]以及多場(chǎng)可編程門陣列(multi-field-programmable gate array,multi-FPGA)并行計(jì)算[17]等.然而,利用高性能計(jì)算硬件在提高計(jì)算效率的同時(shí)增加了檢測(cè)系統(tǒng)的成本,難以在軌道交通運(yùn)維市場(chǎng)進(jìn)行大面積推廣與工程應(yīng)用.另一方面,采用陣元稀疏算法,減少陣元使用數(shù)量,構(gòu)建稀疏矩陣數(shù)據(jù),在保證成像質(zhì)量的條件下,提高計(jì)算效率,保證成像的實(shí)時(shí)性[18].李明飛等[19]對(duì)比研究了遺傳算法和粒子群算法優(yōu)化光學(xué)相控陣天線配置,分別實(shí)現(xiàn)了不同種群數(shù)量遺傳算法和粒子群算法迭代優(yōu)化.Bray 等[20]利用遺傳算法進(jìn)行線性陣列稀疏設(shè)計(jì),獲得了較好的陣列聲學(xué)特性.Peng 等[21]使用稀疏的方法對(duì)單層結(jié)構(gòu)進(jìn)行檢測(cè),同時(shí)增大陣列孔徑尺寸,提高了成像效率,由于陣列孔徑尺寸發(fā)生變化,檢測(cè)時(shí)聲場(chǎng)也會(huì)隨之發(fā)生變化,導(dǎo)致成像的質(zhì)量不高.以上各種方法在解決超聲陣列稀疏優(yōu)化問題時(shí),均采用二進(jìn)制編碼方式,陣元位置一般是周期性等間距排列,陣元自由度低,聲場(chǎng)指向性和聲束聚焦性能差,難以獲得最佳的陣列聲學(xué)性能.

因此,針對(duì)非連續(xù)阻抗粘接因聲阻抗差異大,無法精確計(jì)算聲波在各層之間的傳播時(shí)間,以及全矩陣數(shù)據(jù)進(jìn)行超聲成像計(jì)算耗時(shí),難以實(shí)時(shí)成像、檢測(cè)效率低的問題,本文采用射線追蹤方法,利用直線掃描轉(zhuǎn)換技術(shù)Bresenham 算法,Snell 定律和Fermat 原理,建立射線追蹤的映射過程,準(zhǔn)確計(jì)算超聲波在介質(zhì)中不同層之間的傳播時(shí)間,提高檢測(cè)精度.構(gòu)建基于實(shí)數(shù)編碼的稀布陣列合成孔徑聚焦技術(shù)(synthetic aperture focusing technique,SAFT)成像方法,以陣列方向圖的最小化峰值旁瓣電平為目標(biāo),采用改進(jìn)的遺傳算法優(yōu)化陣元位置,陣元位置呈非周期性排列,陣元自由度顯著提高,增強(qiáng)了聲場(chǎng)指向性和聲束聚焦性能,提高了檢測(cè)效率.

2 非連續(xù)阻抗粘接結(jié)構(gòu)的超聲SAFT成像方法

2.1 多層復(fù)合結(jié)構(gòu)的聲速模型

在多層結(jié)構(gòu)中,當(dāng)聲波以一定的入射角到達(dá)多層固體介質(zhì)分層界面時(shí),會(huì)發(fā)生反射、折射等現(xiàn)象,多層介質(zhì)每層的聲場(chǎng)都可能是多個(gè)單波聲場(chǎng)的疊加,層內(nèi)波的類別與聲波入射角度有關(guān).多層介質(zhì)中聲波的傳播路徑如圖1 所示.多層結(jié)構(gòu)中存在多個(gè)水平層,編號(hào)由上向下依次為1,2,···,n.α,θ分別表示聲波的入射角和折射角,v表示聲波在各層的波速.Z軸為垂直方向,X軸為層狀平行方向.

圖1 多層結(jié)構(gòu)中聲波的傳播Fig.1.Sound wave propagation in multilayer structure.

多層結(jié)構(gòu)中,在不考慮Snell 定理的情況下,超聲波從S點(diǎn)可以通過多條傳播路徑傳播至多層結(jié)構(gòu)中的第二層中的缺陷點(diǎn)F點(diǎn),如圖2 所示.像素點(diǎn)在第二層中時(shí),超聲波在多層結(jié)構(gòu)內(nèi)的傳播路徑為超聲波從激勵(lì)點(diǎn)S以α角入射到第一層中,入射波在第一層中傳播至第一層與第二層界面處的R點(diǎn).超聲波在R點(diǎn)處發(fā)生角度為θ的折射后傳播至第二層中的缺陷F點(diǎn).已知入射點(diǎn)S和F點(diǎn)坐標(biāo)時(shí),Snell 定律為

圖2 聲波在多層結(jié)構(gòu)內(nèi)的傳播路徑Fig.2.Sound wave propagation path in the multi-layer structure.

其中,v1和v2分別表示第一層和第二層中超聲波的波速.根據(jù)Snell 定律,結(jié)合圖3 可得

圖3 聲波在多層結(jié)構(gòu)中的傳播示意圖 (a) F 在第二層內(nèi);(b) F 在第n 層內(nèi)Fig.3.Schematic diagram of sound wave propagation in multi-layer structure: (a) F is in the second layer;(b) F is in the nth layer.

聯(lián)立(1)式,得出聲波在多層結(jié)構(gòu)內(nèi)的傳播聲時(shí)T F:這里,v1表示超聲波折射前的波速,v2表示超聲波折射后的波速,(xS,zS) 表示入射點(diǎn)S的坐標(biāo),(xR,zR)表示折射點(diǎn)R的坐標(biāo),(x,z) 表示F點(diǎn)坐標(biāo).

(3)式是關(guān)于x R的四次方程,當(dāng)像素點(diǎn)F在第二層內(nèi)時(shí),求解(3)式即可得到聲波從第一層傳播至第二層時(shí)折射點(diǎn)位置坐標(biāo).當(dāng)像素點(diǎn)F在第n層內(nèi)時(shí),由A點(diǎn)傳播至第n層要穿過n– 1 個(gè)界面,如圖3(b)所示,發(fā)生n– 1 次折射,此時(shí)需要求解關(guān)于折射點(diǎn)坐標(biāo)的多元四次方程,才能計(jì)算出像素點(diǎn)在第n層時(shí)的折射點(diǎn)坐標(biāo).求解多層結(jié)構(gòu)聲速解析計(jì)算模型可準(zhǔn)確獲得超聲波在多層結(jié)構(gòu)中的傳播時(shí)間.由于無砟軌道結(jié)構(gòu)中軌道板、砂漿層和底座板的聲阻抗差異較大,因此可以將其等效為復(fù)合彈性板結(jié)構(gòu),在該多層結(jié)構(gòu)上使用超聲線性陣列的橫波探頭建立激發(fā)和接收的聲傳播模型.將砂漿層與底座板之間的脫粘缺陷視為局部脫空,超聲陣列的理論模型如圖4 所示.由于超聲波在界面處會(huì)發(fā)生折射,超聲波在無砟軌道結(jié)構(gòu)內(nèi)的傳播路徑不再是兩點(diǎn)之間的直線距離.因此,準(zhǔn)確計(jì)算出超聲波在無砟軌道中傳播時(shí)間的關(guān)鍵在于如何找到真正的超聲波折射路徑.

圖4 無砟軌道脫粘缺陷模型Fig.4.Debonding defect model of ballastless track.

2.2 多層結(jié)構(gòu)射線追蹤SAFT 成像

根據(jù)Snell 定律可知,超聲波在從一種介質(zhì)傳播到另一種介質(zhì)時(shí),在兩層不同介質(zhì)的界面處勢(shì)必產(chǎn)生折射,超聲波在各層內(nèi)的折射線段就組成了其實(shí)際傳播路徑.脫粘缺陷精確表征的難點(diǎn)在于精確求得聲波的傳播時(shí)延,即準(zhǔn)確求得聲波的傳播路徑,而射線追蹤技術(shù)正是解決此問題的有效方法之一.基于Snell 定律和Fermat 原理的射線追蹤方法的原理是通過迭代運(yùn)算尋找超聲波傳播耗時(shí)最短的路徑.該方法在路徑選擇時(shí),需要計(jì)算全部目標(biāo)像素點(diǎn)所對(duì)應(yīng)的所有可能路徑之后再進(jìn)行迭代選擇.

無砟軌道是三層結(jié)構(gòu),需要將目標(biāo)像素點(diǎn)處在第一層(軌道板)、第二層(砂漿層)和第三層(底座板)進(jìn)行分塊計(jì)算,如圖5 所示.

圖5 目標(biāo)像素點(diǎn)示意圖Fig.5.Schematic diagram of target pixels diagram.

目標(biāo)像素點(diǎn)F(x,z)在軌道板內(nèi),此時(shí)不存在折射現(xiàn)象,只需利用軌道板中的聲波速度即可求得超聲波的傳播時(shí)間T(x,z) :

目標(biāo)像素點(diǎn)F(x,z)在砂漿層內(nèi),超聲波從激勵(lì)點(diǎn)S(xS,zS)傳播至軌道板與砂漿層界面發(fā)生折射后,傳播至目標(biāo)像素點(diǎn)F(x,z).將軌道板與砂漿層的界面劃分成M個(gè)折射點(diǎn) (R1,R2,···Ri,···RM) .目標(biāo)像素點(diǎn)F(x,z)在砂漿層時(shí),總的傳播時(shí)間可以由超聲波在折射線上的傳播時(shí)間相加求得.根據(jù)Fermat原理,最短的傳播時(shí)間為超聲波由點(diǎn)S(xS,zS)傳播至F(x,z)的時(shí)間,即

式中,TSR和T RF分別表示超聲波從發(fā)射點(diǎn)S(xS,zS)傳播至軌道板-砂漿層界面的傳播時(shí)間以及超聲波在軌道板-砂漿層界面折射后傳播至砂漿層中目標(biāo)像素點(diǎn)F(x,z)的傳播時(shí)間.

目標(biāo)像素點(diǎn)F(x,z)在底座板內(nèi),超聲波從發(fā)射點(diǎn)S(xS,zS)傳播至軌道板與砂漿層界面時(shí)發(fā)生折射后,傳播至砂漿層中,最后在砂漿層與底座板界面處發(fā)生折射后進(jìn)入底座板中,在這個(gè)過程中,超聲波發(fā)生了兩次折射,傳播路徑分為三段.軌道板-砂漿層界面可以離散成M個(gè)折射點(diǎn),砂漿層與底座板界面可以離散成P個(gè)折射點(diǎn)(K1, K2,···Kj,··· KP).將超聲波在三段路徑上的時(shí)間相加,即可得到目標(biāo)像素點(diǎn)F(x,z)在底座板的情況下超聲波的傳播時(shí)間.超聲波由發(fā)射點(diǎn)至目標(biāo)像素點(diǎn)的傳播時(shí)間可以根據(jù)Fermat 原理求得:

無砟軌道層間聲速差異較大,采用恒定聲速對(duì)脫粘缺陷進(jìn)行檢測(cè)的誤差大,通過射線追蹤能有效求得聲波的傳播路徑和聲波的傳播時(shí)間,檢測(cè)精度得到明顯提升.但是,在成像時(shí)每個(gè)像素點(diǎn)都需要單獨(dú)求解,射線追蹤計(jì)算過程復(fù)雜、運(yùn)算時(shí)長(zhǎng),難以實(shí)現(xiàn)實(shí)時(shí)成像,檢測(cè)效率低.為了提高檢測(cè)效率,采用智能優(yōu)化算法,減少陣元使用數(shù)量,在保證成像質(zhì)量的條件下,提高計(jì)算效率,保證成像的實(shí)時(shí)性.

3 稀布陣列優(yōu)化方法

使用超聲陣列進(jìn)行檢測(cè)時(shí),采用線性陣元周期等間距排列且傳統(tǒng)的智能優(yōu)化算法對(duì)周期陣列進(jìn)行稀疏設(shè)計(jì),本質(zhì)上是在周期陣列中抽取若干陣元進(jìn)行激發(fā)/接收,陣元自由度低,聲場(chǎng)指向性和聲束聚焦性能差,無法做到陣列中特定陣元激發(fā)/接收.稀布陣列的陣元在滿足最小陣元間距條件下,呈非均勻分布于陣列孔徑范圍內(nèi),陣元排布自由度高,聲場(chǎng)指向性和聲束聚焦性能好,能有效減少數(shù)據(jù)量,降低計(jì)算復(fù)雜度,提高檢測(cè)效率.

3.1 稀布陣列優(yōu)化算法

遺傳算法是一種具有高效并且能夠在全局范圍內(nèi)并行的搜索方法,在搜索過程中首先要對(duì)具體問題的解空間進(jìn)行規(guī)劃,優(yōu)化問題搜索空間采用二進(jìn)制編碼方式轉(zhuǎn)換生成一組可行的解.在遺傳算法中盡管二進(jìn)制編碼在實(shí)施的過程中操作簡(jiǎn)單,但是無法滿足一些連續(xù)性或維度高的問題的求解,這是因?yàn)閷?duì)于連續(xù)性問題離散成二進(jìn)制化問題時(shí)本身就有誤差,從而影響了求解的精度.另外,二進(jìn)制編碼的稀疏陣列的陣元間距通常為半波長(zhǎng)的整數(shù)倍,該布置方式的自由度還有待提高.對(duì)于二進(jìn)制編碼的弊端可以采用實(shí)數(shù)編碼進(jìn)行線性陣列的稀布布陣來解決.稀布陣列的陣元在滿足最小陣元間距條件下,陣元位置隨機(jī)分布,陣元排布自由度高.

獲得稀布陣列中陣元位置的具體過程為: 首先在進(jìn)行遺傳操作之前要對(duì)需要優(yōu)化的問題進(jìn)行實(shí)數(shù)編碼,通過實(shí)數(shù)編碼的形式把需要優(yōu)化的問題轉(zhuǎn)化成為計(jì)算機(jī)能夠識(shí)別的方式,變量參數(shù)的范圍為[0,L ?(N ?1)dc]內(nèi)的任意實(shí)數(shù),其中dc表示相鄰兩陣元最小距離間隔,通常選取為半個(gè)波長(zhǎng).設(shè)置初始化種群中的個(gè)體數(shù),對(duì)種群適應(yīng)度函數(shù)進(jìn)行初始化并計(jì)算初始化種群中所有個(gè)體的適應(yīng)度值.然后進(jìn)行遺傳算子操作,得到進(jìn)化生成新的種群.最后,當(dāng)達(dá)到給出的最大迭代次數(shù)后,最終得到一個(gè)具有最低旁瓣水平的稀布陣列.

3.2 稀布陣列設(shè)計(jì)模型

稀布陣列的陣元位置在孔徑范圍內(nèi)隨機(jī)分布,由于方向圖中的峰值旁瓣電平與線性陣列中傳感器分布位置的非線性關(guān)系較強(qiáng),很難采用解析法求解出最優(yōu)解.鑒于此,建立稀布陣列優(yōu)化模型,在相同的陣列孔徑、陣元數(shù)目、最小陣元間距條件下,顯著降低峰值旁瓣電平.假設(shè)陣列的孔徑為L(zhǎng),陣元的個(gè)數(shù)為N,第m個(gè)陣元相對(duì)于第一個(gè)陣元的位置為d m,入射波與線性陣列法向的夾角為θ,主波束的指向角度為θ0,陣列的方向圖函數(shù)為[22]

為了確保稀布陣列的孔徑與滿陣的孔徑相等,陣列左右兩端必須要布置陣元,由此可得:

為了更好地理解稀布陣列陣元之間的間距,將陣元間距拆分為d=x+kdc的形式,

聯(lián)立(10)式—(12)式可以推導(dǎo)出:

通過上述變換,實(shí)現(xiàn)了將陣元間距的計(jì)算轉(zhuǎn)化為對(duì)xi的求解,將搜索空間的范圍由[0,L]壓縮為[0,L ?(N ?1)dc].綜上所述,優(yōu)化模型可以表述為: 尋找滿足(13)式解的集合,使最大旁瓣水平盡可能小.

3.3 多層結(jié)構(gòu)稀布陣列SAFT 成像

在維持相同的陣列孔徑、陣元數(shù)目以及最小陣元間距的前提下,基于遺傳算法可以設(shè)計(jì)出合適的稀布陣列,由于優(yōu)化前后的孔徑大小不變,因而優(yōu)化前后的點(diǎn)散射函數(shù)也保持不變.基于射線追蹤的多層結(jié)構(gòu)稀布SAFT 成像過程為: 首先,將目標(biāo)檢測(cè)區(qū)域劃分成若干個(gè)像素點(diǎn),根據(jù)像素點(diǎn)坐標(biāo)值范圍計(jì)算該點(diǎn)位于多層結(jié)構(gòu)中的某一結(jié)構(gòu)層;然后,采用射線追蹤方法精確計(jì)算每?jī)蓚€(gè)激發(fā)和接收換能器下所有像素點(diǎn)對(duì)應(yīng)的超聲波傳播時(shí)間,并提取對(duì)應(yīng)的采樣數(shù)據(jù),對(duì)所有采樣數(shù)據(jù)進(jìn)行累加;最后,得到圖像范圍內(nèi)所有像素點(diǎn)的像素值為

式中,NT表示稀布發(fā)射陣元數(shù),NR表示稀布接收陣元數(shù);WT表示稀布發(fā)射陣元的權(quán)重,WR表示稀布接收陣元的權(quán)重;Srt表示稀布陣元NT激發(fā)的信號(hào)被稀布陣元NR所接收的信號(hào);Trt表示聲波由稀布發(fā)射陣元出發(fā),經(jīng)缺陷散射后到達(dá)稀布接收陣元的總傳輸時(shí)間.

4 實(shí) 驗(yàn)

實(shí)驗(yàn)采用計(jì)算機(jī)配置為中央處理器Intel(R)Core(TM)i5-4590 3.30 GHz,內(nèi)存容量為8 GB,硬盤容量為1 TB.實(shí)驗(yàn)地點(diǎn)為滬杭高鐵金山北站,實(shí)驗(yàn)對(duì)象為運(yùn)營(yíng)中的CRTS II (China Railway TrackSystem II) 型無砟軌道,檢測(cè)其中的砂漿層與底座板之間的脫粘缺陷,該結(jié)構(gòu)具體材料參數(shù)詳見表1.使用A1040 MIRA 混凝土超聲斷層掃描儀進(jìn)行數(shù)據(jù)采集,該儀器的主要參數(shù)在表2 列出.圖6 為A1040 MIRA 混凝土超聲斷層掃描儀對(duì)CRTS II型無砟軌道砂漿層與底座板之間的脫粘缺陷進(jìn)行探測(cè)的示意圖.脫粘缺陷的厚度約為3 mm,脫粘缺陷的長(zhǎng)度約為260 mm.

表1 CRTS II 型無砟軌道工程材料參數(shù)Table 1.Engineering material parameters of the ballastless track of CRTS II.

表2 A1040 MIRA 混凝土超聲斷層掃描儀主要參數(shù)Table 2.Main parameters of A1040 MIRA concrete ultrasonic tomography scanner.

圖6 脫粘缺陷示意圖Fig.6.Schematic diagram of debonding defects diagram.

5 結(jié)果與分析

以實(shí)數(shù)編碼的稀布陣列優(yōu)化方法為基礎(chǔ),進(jìn)行脫粘缺陷的檢測(cè).根據(jù)優(yōu)化后保留陣元總數(shù)的不同分別做以下分析.當(dāng)稀布陣列中陣元數(shù)量為11 個(gè)時(shí),陣元相對(duì)位置為[0,38,77,107,137,167,197,227,257,299,330] mm,最優(yōu)峰值旁瓣電平Spsl=–15.99 dB,陣列優(yōu)化過程如圖7 所示.成像時(shí)間為4.88 s,成像結(jié)果及剖面圖分別如圖8 和圖9 所示.

圖7 11 個(gè)陣元稀布陣列優(yōu)化過程 (a)方向圖;(b)迭代曲線;(c)陣元相對(duì)位置Fig.7.Optimization process of 11 elements thinned array: (a) Direction diagram;(b) iterative curve;(c) array element position.

圖8 11 個(gè)陣元稀布陣列成像結(jié)果Fig.8.Imaging results of the thinned array of 11 elements.

圖9 11 個(gè)陣元稀布陣列成像結(jié)果剖面圖 (a)縱向剖面圖;(b)橫向剖面圖Fig.9.Cross section of imaging results of thinned array of 11 elements: (a) Longitudinal cross section;(b) transverse cross section.

通過對(duì)成像圖8 的像素進(jìn)行分析可知,11 個(gè)陣元稀布成像的信噪比為19.04 dB,縱向剖面如圖9(a)所示的第一個(gè)波峰為Z1=228 mm,第二個(gè)波峰Z2=533 mm,因此,脫粘缺陷定位誤差為2 mm,底座板底部偏差為3 mm;橫向剖面圖9(b)的等高線與剖面曲線交點(diǎn)的取值范圍為(–126 mm,128 mm),即表征的脫粘長(zhǎng)度為254 mm,長(zhǎng)度誤差為6 mm.

當(dāng)稀布陣列中陣元數(shù)量為9 個(gè)時(shí),陣元相對(duì)位置為[0,41,91,132,168,202,240,281,330] mm,最優(yōu)峰值旁瓣電平SPSL=–15.35 dB,陣列優(yōu)化過程如圖10 所示.成像時(shí)間為3.37 s,成像結(jié)果及剖面圖分別如圖11 和圖12 所示.

圖10 9 個(gè)陣元稀布陣列優(yōu)化過程 (a)方向圖;(b)迭代曲線;(c)陣元相對(duì)位置Fig.10.Optimization process of 9 elements thinned array:(a) Direction diagram;(b) iterative curve;(c) array element position.

通過對(duì)成像圖11 的像素進(jìn)行分析可知,9 個(gè)陣元稀布成像的信噪比為17.46 dB,縱向剖面如圖12(a)所示的第一個(gè)波峰為Z1=227 mm,第二個(gè)波峰Z2=534 mm,因此,脫粘缺陷定位誤差為3 mm,底座板底部偏差為4 mm.橫向剖面如圖12(b)所示的等高線與剖面曲線交點(diǎn)的取值范圍為(–135 mm,113 mm),即表征的脫粘長(zhǎng)度為248 mm,長(zhǎng)度誤差為12 mm.

圖11 9 個(gè)陣元稀布陣列成像結(jié)果Fig.11.Imaging results of thinned array of 9 elements.

圖12 9 個(gè)陣元稀布陣列成像結(jié)果剖面圖 (a)縱向剖面圖;(b)橫向剖面圖Fig.12.Cross section of imaging results of thinned array of 9 elements: (a) Longitudinal cross section;(b) transverse cross section.

6 結(jié)論

無砟軌道是由軌道板、砂漿層、底座板和基床組成的非連續(xù)阻抗混凝土粘接結(jié)構(gòu),層間聲速差異顯著.傳統(tǒng)的超聲檢測(cè)采用恒定聲速進(jìn)行成像,檢測(cè)誤差大、效率低.建立多層結(jié)構(gòu)聲速模型,采用射線追蹤方法準(zhǔn)確計(jì)算超聲波的傳播路徑和傳播時(shí)間.提出實(shí)數(shù)編碼的稀布陣列方法,利用多層結(jié)構(gòu)合成孔徑方法采集超聲信號(hào),進(jìn)行成像,分析結(jié)果得出以下結(jié)論.

1)針對(duì)無砟軌道層間聲阻抗差異較大,將其視為非均勻聲速?gòu)?fù)合彈性板結(jié)構(gòu),建立多層聲速理論模型,采用射線追蹤方法精確求得聲波在層間聲速差異大的介質(zhì)中的傳播路徑和走時(shí),提高了檢測(cè)精度.

2)采用實(shí)數(shù)編碼的稀布陣列方法,進(jìn)行SAFT成像.實(shí)驗(yàn)結(jié)果表明: 稀布陣列具有更窄的主瓣寬度,更低的旁瓣增益,陣元自由度得到顯著提高.當(dāng)成像誤差在±5%以內(nèi)時(shí),稀布陣列SAFT 成像方法的效率提高了30.9%.

3)隨著有效陣元個(gè)數(shù)的減少,稀布陣列下的旁瓣增益都在逐漸增大,脫粘缺陷的檢測(cè)精度也隨之減低,為保證成像質(zhì)量和檢測(cè)精度,有效陣元數(shù)量必須位于合理的區(qū)間范圍內(nèi).