基于k-wave超聲場(chǎng)時(shí)域仿真研究*

李 雅，陳友興，任陽(yáng)山，王召巴

(1.中北大學(xué) 信息與通信工程學(xué)院，山西 太原 030051；2.海軍青島航保修理廠，山東 青島266071)

基于k-wave超聲場(chǎng)時(shí)域仿真研究*

李 雅1，陳友興1，任陽(yáng)山2，王召巴1

(1.中北大學(xué) 信息與通信工程學(xué)院，山西 太原 030051；2.海軍青島航保修理廠，山東 青島266071)

超聲場(chǎng)仿真研究在工業(yè)超聲探傷領(lǐng)域一直起到舉足輕重的作用。本文采用脈沖超聲波檢測(cè)原理，建立仿真平臺(tái)，利用k-wave工具箱對(duì)超聲場(chǎng)進(jìn)行了時(shí)域仿真研究，針對(duì)聲場(chǎng)在圓柱體工件內(nèi)部的傳播情況，并利用聲學(xué)傳感器接收到的回波信號(hào)，采用插值重構(gòu)模型，對(duì)工件內(nèi)部超聲場(chǎng)分布做了重構(gòu)研究和可視化實(shí)現(xiàn)，為實(shí)際工業(yè)生產(chǎn)奠定了理論基礎(chǔ)。

k-wave；超聲場(chǎng)；聲學(xué)傳感器；重構(gòu)模型

0 引言

在無損探傷過程中，探頭處發(fā)射的超聲波在缺陷中經(jīng)過發(fā)射被超聲波探傷儀接收，以確定缺陷的位置和類型。無損探傷具有檢測(cè)成本低、速度快和現(xiàn)場(chǎng)使用方便等優(yōu)點(diǎn)，多用于管道、壓力容器等圓柱體構(gòu)件的探傷檢測(cè)[1]。計(jì)算超聲學(xué)在近些年應(yīng)用于超聲波檢測(cè)的模擬，從不同尺度上入手，建立不同的超聲波模型，來模擬聲束在介質(zhì)中傳播時(shí)的聲場(chǎng)、缺陷與界面處對(duì)超聲波的作用規(guī)律，構(gòu)建超聲波檢測(cè)計(jì)算模型，整合發(fā)展超聲波檢測(cè)的技術(shù)，促進(jìn)超聲波檢測(cè)在無損檢測(cè)中的應(yīng)用。對(duì)于許多復(fù)雜情況而言，理論與實(shí)驗(yàn)手段難以進(jìn)行，但可以利用計(jì)算機(jī)環(huán)境進(jìn)行仿真計(jì)算得到數(shù)值結(jié)果[2]。本文利用MATLAB加載的k-wave工具箱模擬聲場(chǎng)傳播和建立缺陷重構(gòu)模型，使得原有聲場(chǎng)建模變得簡(jiǎn)單、快速[3]。

1 超聲場(chǎng)仿真原理

本文在均勻背景散射媒介中建立聲源模型，簡(jiǎn)化輸入?yún)?shù)在提高計(jì)算效率的同時(shí)，排除非關(guān)鍵因素，從而有效提取出超聲場(chǎng)的精確回波模型。

1.1 聲速與聲壓的計(jì)算

由流體介質(zhì)中的非線性壓力與密度關(guān)系，可以將狀態(tài)方程p?=p?(ρ?，s?)擴(kuò)展為泰勒級(jí)數(shù)中對(duì)應(yīng)的壓力項(xiàng)和密度項(xiàng)。假設(shè)非線性和壓力的變化影響(由于聲波吸收)都是二階。因此，高階壓力項(xiàng)可以被丟棄?？紤]到在很小的有限時(shí)間步長(zhǎng) δt=t1-t0內(nèi)改變流體元素的總壓強(qiáng)，泰勒級(jí)數(shù)展開可寫成：

式(1)中最后一部分代表一種能量損失也是一種聲能吸收。在熱粘性介質(zhì)中，可以結(jié)合介質(zhì)的導(dǎo)熱系數(shù)和比熱容，考慮能量守恒方程[4]，這一點(diǎn)可以表示為一種聲學(xué)損失模型：

L表示聲能損失。建立冪律吸收模型，基于分?jǐn)?shù)拉普拉斯算子定義 L的表達(dá)式為[5]：

τ和 η是吸收和色散比例系數(shù)[6]。

1.2 k-space偽譜法

偽譜法求解微分方程主要包含：(1)離散化處理，用一定的方式在整個(gè)區(qū)域上選取格點(diǎn)構(gòu)成一個(gè)網(wǎng)絡(luò)。而其它點(diǎn)上的變量值可以用插值的方法得到；(2)選取合適的基函數(shù)變量值的情況下構(gòu)造插值函數(shù)，計(jì)算格點(diǎn)上變量值的空間導(dǎo)數(shù)。偽譜法是全局方法,只能用在空間域。如果所求解的是含時(shí)微分方程，一般用簡(jiǎn)單的有限差分方法來處理時(shí)間域的計(jì)算[7]。

在k-wave中，利用快速傅里葉變換和 k-space偽譜法，將守恒方程中的粒子速度和聲學(xué)密度轉(zhuǎn)換為離散形式，計(jì)算每個(gè)時(shí)間步長(zhǎng)內(nèi)聲壓場(chǎng)的變化情況，能夠高效準(zhǔn)確地對(duì)超聲場(chǎng)進(jìn)行數(shù)值模擬。離散表達(dá)式如式(4)～式(7)所示：

聲學(xué)密度是分布在笛卡兒坐標(biāo)系中并可以引入各向異性 PML應(yīng)用的值[8]。F和 F-1表示空間傅里葉正逆變換，上標(biāo)n和n+1分別表示函數(shù)在當(dāng)前和下一時(shí)間點(diǎn)的值，i為虛數(shù)單位，Δt是時(shí)間步長(zhǎng)，kξ是 ξ處波數(shù)，κ= sinc(crefkΔt/2)，其中 cref為相關(guān)聲速。

式(7)離散方程可利用基于 Courant-Friedrichs-Lewy的 CFL數(shù)值求解時(shí)間步長(zhǎng)，Δt=CFLΔx/cmax。CFL通常取0.3能夠平衡準(zhǔn)確性和計(jì)算效率之間的關(guān)系。在每個(gè)時(shí)間步長(zhǎng)內(nèi)，質(zhì)量或力源可以通過在計(jì)算域內(nèi)添加適當(dāng)?shù)木W(wǎng)格點(diǎn)值來設(shè)置。同樣的，模擬的輸出可以通過每個(gè)時(shí)間步長(zhǎng)在特定的網(wǎng)格點(diǎn)處的聲變量記錄[9-10]。

2 聲場(chǎng)仿真模型

本文設(shè)計(jì)圓柱體工件超聲場(chǎng)模型，并通過傳感器接受回波信號(hào)。工件的設(shè)定如下：實(shí)心圓工件半徑為2 mm；人造缺陷孔工件中心圓孔工件外半徑為2 mm，內(nèi)半徑為1 mm；偏心圓孔工件外半徑2 mm，偏心孔中心為半徑的中間，半徑為0.5 mm，在距離圓柱的中心半徑為4.5 mm處設(shè)置8個(gè)等距離分布的傳感器。時(shí)間步長(zhǎng)為1 μs，由于超聲波傳播速度很快且工件尺寸很小，本文對(duì)超聲波在工件內(nèi)和水中聲速的差異不做考慮，超聲波傳播速度為1 500 m/s。

根據(jù)惠更斯原理，聲波遇到障礙物時(shí)會(huì)以障礙物作為新的聲源發(fā)射。本文仿真環(huán)境以圓柱體工件為新的聲源，并設(shè)定了8個(gè)點(diǎn)傳感器的分布位置及其編號(hào)，工件將360°接收的單脈沖信號(hào)在同一時(shí)刻發(fā)射出去，如圖1所示為不同構(gòu)件發(fā)射聲場(chǎng)模型的二維顯示圖，從中可以很直觀地觀測(cè)到聲場(chǎng)傳播過程中波形和聲壓強(qiáng)弱情況，深色到淺色區(qū)域即為聲壓從強(qiáng)到弱體現(xiàn)。

圖1 圓柱體工件及其發(fā)射聲場(chǎng)模型

圖2 為8個(gè)點(diǎn)傳感器探頭接收的超聲信號(hào)三維圖，通過與實(shí)心圓對(duì)比可以明顯看出中心圓孔和偏心圓孔發(fā)射聲場(chǎng)聲壓的異同之處：(1)由于圓柱位置及外半徑設(shè)置相同，傳感器探頭接收信號(hào)的起始和終止位置一樣；(2)由于中心圓孔的設(shè)定位置在傳感器探頭的正中心，所以8個(gè)探頭接收到的信號(hào)同一時(shí)刻達(dá)到波峰；(3)根據(jù)偏心圓孔設(shè)定的位置，5號(hào)最近點(diǎn)傳感器會(huì)首先接收到回波信號(hào)，1號(hào)最遠(yuǎn)點(diǎn)傳感器則會(huì)最后接收到超聲的二次回波信號(hào)。

3 逆時(shí)反演重構(gòu)模型

利用聲學(xué)傳感器接收到的回波信號(hào)逆時(shí)反演重構(gòu)，在圓周上每隔一度設(shè)置一個(gè)點(diǎn)傳感器，使之圍繞工件一周，充分接收到不同工件 360°的發(fā)射聲場(chǎng)，利用接收到的發(fā)射信號(hào)進(jìn)行逆時(shí)反演重構(gòu)，仿真重構(gòu)效果如圖3(b)、圖3(d)～圖3(f)所示。如圖 3(c)，圖3(e)所示，聲壓在人造缺陷孔的邊界位置將會(huì)有明顯變動(dòng)，缺陷邊界處聲壓坐標(biāo)值如表1所示。根據(jù)測(cè)量數(shù)據(jù)可知，在誤差允許范圍內(nèi)缺陷仿真邊界坐標(biāo)點(diǎn)位置與實(shí)際缺陷設(shè)定尺寸相一致。

圖2 8個(gè)傳感器接收回波信號(hào)三維顯示圖

圖3 圓柱體工件截面聲壓圖

表1 缺陷邊界位置聲壓統(tǒng)計(jì)表

4 總結(jié)

本文利用k-wave工具箱對(duì)三種典型的工件模型做了聲場(chǎng)仿真研究，設(shè)計(jì)的仿真平臺(tái)能夠達(dá)到實(shí)際檢測(cè)中難以實(shí)現(xiàn)的設(shè)定，能夠?qū)崿F(xiàn)360°同時(shí)刻采集工件的發(fā)射信號(hào)，減少人為移動(dòng)工件或探頭進(jìn)行超聲探測(cè)掃描帶來的誤差等方面的影響因素，并且能夠?qū)⒙晥?chǎng)能量具體化，直觀可視化地觀測(cè)聲場(chǎng)在工件中具體傳播過程。通過重構(gòu)結(jié)果分析，測(cè)得缺陷的位置和尺寸數(shù)據(jù)與工件的實(shí)際設(shè)定值相一致，充分說明仿真理論的正確性和仿真結(jié)果的準(zhǔn)確性。

[1]孟凡凱.水下超聲波檢測(cè)聲場(chǎng)數(shù)值模擬與實(shí)驗(yàn)分析[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2012.

[2]COX B T，TREEBY B E.Effect of Sensor Directionality on Photoacoustic Imaging:A Study Using the k-Wave Toolbox. Proc.of SPIE Vol.2010，7564，75640I1-75640I6.

[3]TREEBY B E，COX B T.K-Wave：a MATLAB toolbox for the simulation and reconstruction of photoacoustic wavefields[J].J.Biomed.Opt.，2010，15(2)，0213141-02131412.

[4]COX B T，KARA S，ARRIDGE S R，et al.K-space propagation models for acoustically heterogeneous media：Application to biomedical photoacoustics[J].J.Acoust.Soc.Am. 2007，121(40)：3453-3464.

[5]《超聲波探傷》編寫組.超聲波探傷[M].北京：電力工業(yè)出版社，1980.

[6]Tabei M，Mast T D，Waag R C.A k-space method for coupled first-order acoustic propagation equations[J].Acoust. Soc.Am，2002，111(1)：53-63.

[7]劉魯波，陳曉非，王彥賓.切比雪夫偽譜法模擬地震波場(chǎng)[J].西北地震學(xué)報(bào)，2007（01）：18-25.

[8]COX B T，KARA S，ARRIDGE S R，et al.K-space propagation models for acoustically heterogeneous media：Application to biomedical photoacoustics[J].Acoust.Soc.Am，2007，121(6)：3453-3464.

[9]Tillett J C，Daoud M I，Lacefield J C，et al.A k-space method for acoustic propagation using coupled first-order equations in three dimensions[J].Acoust.Soc.Am，2009，126 (3)：1231-1244.

[10]TREEBY B E，JAROS J，RENDELL A P.Modeling nonlinear ultrasound propagation inheterogeneo us media with power law absorption using a k-space pseudospectral method[J]. Acoust.Soc.Am，2012，131(6)：4324-4336.

Based on the k-wave ulrasonic field in time domain simulation

Li Ya1，Chen Youxing1，Ren Yangshan2，Wang Zhaoba1
(1.Information and Communication Engineering,North University of China,Taiyuan 030051，China；2.Navy Maritime Security Repair Factory,Qingdao 266071，China)

Ultrasonic field simulation research is playing a very important role in the field of industrial ultrasonic flaw detection. Based on pulse principle of ultrasound detection,this paper established the simulation platform,using the k-wave toolbox on ultrasonic field,the time-domain simulation for the transmission of sound field in the cylinder workpiece internal situation.With acoustic sensor receives the echo signal,by using the interpolation reconstruction model,the research on ultrasonic field distribution inside the workpiece made reconstruction and visualization,laying a theoretical foundation for industrial production.

k-wave；ultrasonic field；acoustic sensors；reconstruction model

TP391.9；TB553.8

A

10.16157/j.issn.0258-7998.2015.11.037

李雅，陳友興，任陽(yáng)山，等.基于k-wave超聲場(chǎng)時(shí)域仿真研究[J].電子技術(shù)應(yīng)用，2015，41(11)：132-134，139.

英文引用格式：Li Ya，Chen Youxing，Ren Yangshan，et al.Based on the k-wave ulrasonic field in time domain simulation[J]. Application of Electronic Technique，2015，41(11)：132-134，139.

2015-07-10）

李雅（1990-），通信作者，女，碩士研究生，主要研究方向：信號(hào)與信息處理，E-mail：1131680323@qq.com。

國(guó)家自然科學(xué)基金項(xiàng)目 (61201412),山西省科技攻關(guān)項(xiàng)目(20110321029),山西省青年科技研究基金(2012021011-5)