林?濱，杜安堯，劉再蔚，梁小虎

基于分?jǐn)?shù)階傅里葉變換的聲表面波涂層厚度檢測(cè)

林?濱，杜安堯，劉再蔚，梁小虎

(天津大學(xué)機(jī)構(gòu)理論與裝備設(shè)計(jì)教育部重點(diǎn)實(shí)驗(yàn)室，天津 300354)

聲表面波；頻散曲線；分?jǐn)?shù)階傅里葉變換

聲表面波在半無(wú)限大單一介質(zhì)表面?zhèn)鞑r(shí)，聲速僅與介質(zhì)的力學(xué)性質(zhì)有關(guān)；但當(dāng)傳播介質(zhì)為基體-涂層等非單一介質(zhì)[1-3]，或近表面存在損傷[4-6]、殘余應(yīng)力等[7-8]狀態(tài)時(shí)就會(huì)出現(xiàn)色散現(xiàn)象(不同頻率聲表面波的傳播速度不同)．通過(guò)表面波的時(shí)域信號(hào)獲取其時(shí)頻信息，進(jìn)而獲取頻散曲線來(lái)反演研究介質(zhì)的厚度是一種有效且常用的方法．傅里葉變換和小波變換都可以獲得聲表面波的頻散曲線[9-10]，但傅里葉變換在完整地保留頻率信息的基礎(chǔ)上會(huì)丟失時(shí)間信息，因此傅里葉變換不能對(duì)時(shí)變信號(hào)進(jìn)行有效處理；小波變換雖然同時(shí)在時(shí)域和頻域具有一定的分辨率，可以有效地保留信號(hào)的時(shí)變特性，但小波變換會(huì)降低頻域的分辨率，且頻率分辨率會(huì)隨著頻率的升高而大幅降低；此外小波變換的效果高度依賴于母小波的選?。虼?，本文將針對(duì)新的頻散分析方法進(jìn)行模擬仿真和實(shí)驗(yàn)分析并研究其分析效果．

1?分?jǐn)?shù)階傅里葉變換時(shí)頻分析

分?jǐn)?shù)階傅里葉變換具有多種不同的定義與解?釋[11]：線性積分變換、時(shí)頻面旋轉(zhuǎn)、特征函數(shù)與特征值、坐標(biāo)乘法和求導(dǎo)算子的變換、微分方程的解和超微分算子等．本節(jié)僅從線性積分變換對(duì)分?jǐn)?shù)階傅里葉變換進(jìn)行解釋．

1.1?分?jǐn)?shù)階傅里葉變換理論

式中為角度參數(shù)．和的關(guān)系為

因此可得

1.2?線性調(diào)頻函數(shù)時(shí)頻分析

激光激發(fā)聲表面波可以近似簡(jiǎn)化為線性調(diào)頻信號(hào)，因此可以利用加入高斯白噪聲的線性調(diào)頻函數(shù)來(lái)模擬激光激發(fā)聲表面波，然后利用分?jǐn)?shù)階傅里葉變換對(duì)模擬聲表面波進(jìn)行時(shí)頻分析，最后將分析效果與小波變換分析效果和理論時(shí)頻曲線進(jìn)行對(duì)比．

利用瑞利波包構(gòu)造線性調(diào)頻信號(hào)

式中：，為幅值參數(shù)；，為時(shí)間參數(shù)；，為調(diào)頻參數(shù)；，為瑞利參數(shù).函數(shù)的理論調(diào)頻率，理論中心頻率．為模擬實(shí)際實(shí)驗(yàn)信號(hào)，加入信噪比為－3dB的高斯白噪聲，原信號(hào)與帶噪信號(hào)時(shí)域波形如圖1所示，后續(xù)將對(duì)帶噪信號(hào)進(jìn)行分析．

圖2?帶噪信號(hào)的分?jǐn)?shù)階傅里葉變換

圖3?時(shí)頻分析與對(duì)比

2?基體-涂層材料聲表面波的分?jǐn)?shù)階傅里葉變換

圖4?激光激發(fā)聲表面波實(shí)驗(yàn)示意

表1?實(shí)驗(yàn)材料力學(xué)參數(shù)

Tab.1?Mechanical parameters of materials

為實(shí)現(xiàn)對(duì)不同傳播距離下聲表面波的拾取，將直角棱鏡和柱面透鏡固接在位移精度0.5μm的直線位移平臺(tái)上，通過(guò)驅(qū)動(dòng)器和編程控制器對(duì)五相異步電機(jī)的控制使位移平臺(tái)每次移動(dòng)1mm，進(jìn)行3次．需要注意的是，如圖6所示，實(shí)驗(yàn)過(guò)程中激發(fā)點(diǎn)需要由遠(yuǎn)及近向拾取點(diǎn)移動(dòng)，其目的是為了避免后一次激光激發(fā)聲表面波傳播經(jīng)過(guò)前一次激發(fā)損傷區(qū)域而受干擾．此外，樣品表面上激光激發(fā)點(diǎn)和振動(dòng)拾取點(diǎn)間的距離需要足夠大，使拾取點(diǎn)的振動(dòng)信號(hào)為不受熱傳導(dǎo)影響的穩(wěn)定聲表面波，本實(shí)驗(yàn)中激發(fā)點(diǎn)與拾取點(diǎn)的距離大于2cm．

圖5?壓電薄膜傳感器

圖6?不同激發(fā)位置聲表面波的拾取

圖7?不同拾取點(diǎn)的時(shí)域波形

圖8?不同拾取點(diǎn)聲表面波信號(hào)的分?jǐn)?shù)階傅里葉變換

圖9?不同拾取點(diǎn)聲表面波信號(hào)的時(shí)頻曲線

3?頻散曲線反演

因?yàn)樵嚰閽伖夤杵?，表面粗糙度較小，因此可以忽略粗糙度對(duì)聲表面波頻散曲線的影響[16]．在一般情況下，若材料的結(jié)構(gòu)和各部分的彈性常數(shù)、泊松比以及密度已知，則可以計(jì)算出材料的理論頻散曲線．凡友華等[17]提出了標(biāo)量傳遞法的相速度頻散曲線快速計(jì)算方法，但本文中的實(shí)驗(yàn)頻散曲線為群速度頻散曲線，因此利用群速度和相速度的轉(zhuǎn)換關(guān)系計(jì)算聲表面波理論群速度頻散曲線，即

式(9)為半無(wú)限大空間單一介質(zhì)中聲表面波群速度近似公式，即

式中：為泊松比；為橫波波速，可以通過(guò)彈性模量E和泊松比求得．Si(100)和中聲表面波群速度分別為4876m/s和3411m/s．圖10所示為具有4μm 氧化層的Si(100)的群速度頻散曲線，在低頻極限時(shí)，聲表面波幾乎在Si(100)基體內(nèi)傳播，其速度接近單晶硅中聲表面波速度；在高頻極限時(shí)聲表面波集中于表面的層中，其速度接近于二氧化硅中聲表面波的速度．此外，沿著具有4μm 氧化層的Si(100)的[100]晶向傳播的聲表面波，其頻散曲線的拐點(diǎn)在400MPa附近．當(dāng)頻率范圍限定在[30MHz，50MHz]時(shí)，頻散曲線可以近似為一條直線，因此可以通過(guò)聲表面波時(shí)頻信息對(duì)頻散曲線進(jìn)行線性擬合，進(jìn)而反演出層的厚度．

在[30MHz，50MHz]的頻率范圍內(nèi)，利用式(10)計(jì)算聲表面波群速度實(shí)驗(yàn)頻散曲線，即

將實(shí)驗(yàn)頻散曲線與理論頻散曲線對(duì)比，構(gòu)造頻散曲線失配函數(shù)[18]，即

圖11?激光聲表面波群速度頻散曲線

可以看出，反演計(jì)算得到的實(shí)驗(yàn)頻散曲線與理論頻散曲線在低頻段(小于35MHz)有一定的誤差．這是因?yàn)樵瓡r(shí)域信號(hào)中噪聲引起的誤差使得分?jǐn)?shù)階傅里葉變換之后形成脈沖點(diǎn)的階次產(chǎn)生一定的誤差，這個(gè)誤差會(huì)相應(yīng)地反映在信號(hào)的時(shí)頻曲線和頻散曲線上．此外，聲表面波時(shí)頻曲線和頻散曲線本質(zhì)上都是非線性曲線，對(duì)聲表面波時(shí)頻曲線進(jìn)行線性擬合時(shí)也忽略了信號(hào)的高階誤差．如要消除圖11中35MHz以下頻散曲線的誤差，就需要在對(duì)信號(hào)進(jìn)行分?jǐn)?shù)階傅里葉變換之前對(duì)信號(hào)進(jìn)行去噪處理，此外利用非線性方法對(duì)聲表面波時(shí)頻曲線和頻散曲線進(jìn)行擬合也可以降低一定的誤差．

4?結(jié)?論

聲表面波的頻散曲線包含了很多被測(cè)試件的性質(zhì)信息，如涂層厚度、彈性模量和泊松比等，因此通過(guò)控制變量法來(lái)計(jì)算頻散曲線就可以獲得基體-涂層材料的涂層厚度．而在獲取聲表面波的頻散曲線的過(guò)程中，對(duì)聲表面波信號(hào)的時(shí)頻分析是最重要的一步，分?jǐn)?shù)階傅里葉變換可以有效地對(duì)聲表面波時(shí)域信號(hào)進(jìn)行分析并獲取其時(shí)頻曲線．

(1) 利用瑞利波包和高斯白噪聲構(gòu)造線性調(diào)頻函數(shù)來(lái)模擬聲表面波信號(hào)，并用FrFT和CWT對(duì)其進(jìn)行時(shí)頻分析．雖然CWT比FrFT具有更低的中心頻率誤差，但FrFT的調(diào)頻率誤差遠(yuǎn)小于CWT，因此FrFT更能準(zhǔn)確地?cái)M合實(shí)際時(shí)頻曲線．

(3) 利用表面波時(shí)頻信息計(jì)算頻散曲線，構(gòu)造失配函數(shù)表征實(shí)驗(yàn)頻散曲線與理論頻散曲線之間的匹配度．通過(guò)對(duì)失配函數(shù)的分析得到實(shí)驗(yàn)材料的氧化層厚度為4.03μm，誤差僅為0.75%．

［1］ Liu Z，Lin B，Liang X，et al. Time-frequency analysis of laser-excited surface acoustic waves based on synchrosqueezing transform[J]. Ultrasonics，2020，106：106147.

［2］ Ollendorf H，Schneider D，Schwarz T，et al. A comparative study of the mechanical properties of TiN coatings using the non-destructive surface acoustic wave method，scratch test and four-point bending test[J]. Surface and Coatings Technology，1996，84(1/2/3)：458-464.

［3］ Hurley D C，Tewary V K，Richards A J. Surface acoustic wave methods to determine the anisotropic elastic properties of thin films[J]. Measurement Science and Technology，2001，12(9)：1486-1494.

［4］ Moumita M，Amares C，Pulkit K，et al. Effect of initial stress，heterogeneity and anisotropy on the propagation of seismic surface waves[J]. Mechanics of Advanced Materials & Structures，2020，27(3)：177-188.

［5］ Yuan M，Kang T，Kim H J，et al. A numerical model for prediction of residual stress using rayleigh waves[J]. Journal of the Korean Society for Nondestructive Testing，2011，31(6)：656-664.

［6］ Mora P，Spies M. Inversion of residual stress profiles from ultrasonic Rayleigh wave dispersion data[J]. Inverse Problems，2018，34(5)：55001.

［7］ Zerwer A，Polak M A，Santamarina J C. Rayleigh wave propagation for the detection of near surface discontinuities：Finite element modeling[J]. Journal of Nondestructive Evaluation，2003，22(2)：39-52.

［8］ Lavrentyev A I，Veronesi W A，Baaklini G Y，et al. Ultrasonic characterization of shot-peened metal surfaces

[J]. International Society for Optics and Photonics，2001，4336：88-99.

［9］ Liang X，Lin B，Liu Z. Analysis of velocity calculation methods of laser-induced surface acoustic wave[J]. Ultrasonics，2019，100：105985.

［10］ Liu X，Lin B，Liang X. Study of laser-induced surface acoustic wave propagating on materials with machined surfaces based on wavelet analysis[J]. Surface and Coatings Technology，2019，358：173-181.

［11］ Bultheel A，Martínez H. An introduction to the fractional Fourier transform and friends[J]. Cubo—A Mathematical Journal，2005，7(2)：201-221.

［12］ Bultheel A，Héctor E M S. Computation of the fractional Fourier transform[J]. Applied and Computational Harmonic Analysis，2004，16(3)：182-202.

［13］ Ozaktas H M，Arikan O，Kutay M A，et al. Digital computation of the fractional Fourier transform[J]. IEEE Trans Signal Process，1996，44(9)：2141-2150.

［14］ 趙興浩，鄧?兵，陶?然. 分?jǐn)?shù)階傅里葉變換數(shù)值計(jì)算中的量綱歸一化[J]. 北京理工大學(xué)學(xué)報(bào)，2005，25(4)：360-364.

Zhao Xinghao，Deng Bing，Tao Ran. Dimensional normalization in the digital computation of the fractional Fourier transform[J]. Transactions of Beijing Institute of Technology，2005，25(4)：360-364(in Chinese).

［15］ Bogardus H E. Third‐order elastic constants of Ge，MgO，and fused SiO2[J]. Journal of Applied Physics，1965，36(8)：2504-2513.

［16］ Paehler D，Schneider D，Herben M. Nondestructive characterization of sub-surface damage in rotational ground silicon wafers by laser acoustics[J]. Microelectronic Engineering，2007，84(2)：340-354.

［17］ 凡友華，劉家琦. 層狀介質(zhì)中瑞雷面波的頻散研究[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2001，33(5)：577-581.

Fan Youhua，Liu Jiaqi. Research on the dispersion of Rayleigh waves in multilayered media[J]. Journal of Harbin Institute of Technology，2001，33(5)：577-581(in Chinese).

［18］ Lu Y，Peng S，Du W，et al. Rayleigh wave inversion using heat-bath simulated annealing algorithm[J]. Journal of Applied Geophysics，2016，134：267-280.

Film Thickness Detection Using a Surface Acoustic Wave Based on the Fractional Fourier Transform

Lin Bin，Du Anyao，Liu Zaiwei，Liang Xiaohu

(Key Laboratory of Mechanism Theory and Equipment Design of Ministry of Education，Tianjin University，Tianjin 300354，China)

surface acoustic wave；dispersion curve；fractional Fourier transform

10.11784/tdxbz202103050

O426.9

A

0493-2137(2022)07-0714-07

2021-03-24；

2021-05-20.

林?濱（1965—??），男，博士，教授.Email：m_bigm@tju.edu.cn

林?濱，linbinph@tju.edu.cn.

國(guó)家自然科學(xué)基金資助項(xiàng)目(52075379).

the National Natural Science Foundation of China(No.52075379).

(責(zé)任編輯：王曉燕)