樊程廣,潘孟春,羅飛路,羅詩途,譚項林

(國防科學技術大學機電工程與自動化學院,湖南長沙 410073）

0 引 言

纖維增強樹脂基復合材料是應用最為廣泛的復合材料,具有一系列的優(yōu)良性能,因此在航空工業(yè)領域得到了廣泛的應用,尤其是在飛機的制造上,典型應用如直升機槳葉[1-2].在纖維增強樹脂基復合材料廣泛應用的同時,人們對其性能的無損檢測和評估也在同步進行.聲學性能是評價材料的重要指標,從不同的角度來看,材料的機械性能和顯微結構等都能夠利用聲學性能來表征,因此準確測量材料聲學性能是一個關鍵的問題[3].在復合材料的無損檢測技術中,超聲檢測具有靈敏度高、穿透能力強、檢測速度快等優(yōu)點,因此發(fā)展最快、應用最廣泛[4-6].

19世紀末20世紀初,德國著名數(shù)學家David Hilbert提出了Hilbert變換方法,目前,該方法已經(jīng)廣泛應用于各工程領域.例如,Hilbert變換用于故障診斷[7]、信號解調(diào)[8]、時頻分析[9]、電壓閃變測量[10]等.Hilbert變換方法既能克服傅里葉變換方法分辨率低的缺點,又能同時克服在小波變換方法中需要選擇小波基的困難和小波變換后信號能量泄露的缺點,因此,將Hilbert變換引入復合材料聲學性能的測量是一個有意義的研究方向[11].

本文首先詳細介紹了Hilbert變換原理和聲學性能測量原理,然后利用Hilbert變換對實測的超聲回波信號進行分析,進而求得材料相關聲學性能參數(shù),并對不同種類的復合材料進行對比評估.

1 原 理

1.1 Hilbert變換

1.1.1 超聲回波信號的數(shù)學模型

超聲波入射至被測材料,當超聲波遇到非連續(xù)性結構或邊界返回時,產(chǎn)生振蕩波,即為超聲回波信號f(t),超聲回波信號具有窄帶信號的特征[12]

式中:A(t)為超聲回波信號的幅值函數(shù);ω為超聲換能器的固有激發(fā)頻率;φ(t)為相位函數(shù).

從信號分析的角度來看,超聲回波信號是對高頻載波信號的低頻調(diào)制,其中,高頻信號ν(t)為載波信號,緩慢變化的幅值函數(shù) A(t)為調(diào)幅信號,也就是超聲回波信號的包絡,體現(xiàn)了信號的整體起伏特性,包含了被測材料的關鍵信息.

1.1.2 超聲回波信號的Hilbert分析

超聲回波信號 f(t)是實信號,實信號的傅里葉變換包含負的頻率成分,不利于信號分析,因此引入解析信號的概念.Hilbert變換就是一種將時域?qū)嵭盘栕儞Q為時域解析信號的數(shù)學處理過程[13].超聲回波信號 f(r)的Hilbert變換 fh(t)定義

基于超聲回波信號 f(t)的解析函數(shù)fa(t)

式(3)表明:解析函數(shù) fa(t)的實部是超聲回波信號 f(t),虛部是超聲回波信號 f(t)的Hilbert變換 fh(t).

定義 f(ω)為 f(t)的傅里葉變換,fh(ω)為 fh(t)的傅里葉變換,fa(ω)為fa(t)的傅里葉變換,根據(jù)Hilbert變換的性質(zhì)可以得到

式中:sgn(ω)為符號函數(shù).

因此,在已有超聲回波信號 f(t)的前提下可由 f(ω)得到 fa(ω),進而求得解析函數(shù) fa(t).

圖1 模擬超聲回波信號及其Hilbert分析Fig.1 Simulation of the ultrasonic echo signalw ith its hilbert transform

根據(jù)以上分析,解析函數(shù) fa(t)又可表示成

由式(1)和式(6)可知:解析函數(shù) fa(t)和超聲回波信號 f(t)的振幅函數(shù)相同,同為 A(t),相位函數(shù)相同,同為 ωt+φ(t).因此超聲回波信號 f(t)的瞬時包絡為

1.1.3 仿真分析

圖1表示了一個模擬超聲回波信號的Hilbert分析過程,圖1(a)～(d)分別代表模擬的超聲回波信號、基于超聲回波信號的解析信號的實部、基于超聲回波信號的解析信號的虛部、超聲回波信號的包絡,從圖1中可以看出:經(jīng)過Hilbert變換后的超聲回波信號包絡與實際超聲回波信號的包絡一致.

1.2 聲學性能測量原理

聲速、聲阻抗以及衰減系數(shù)是描繪超聲波在被測材料中傳播時的聲學參數(shù),是表征材料成分和組織結構的重要參數(shù),材料的彈性模量、材料的顯微結構、形態(tài)以及相應的機械性能都能用聲學參數(shù)來表征.例如通過測量聲速可以確定材料的彈性模量,根據(jù)聲速和衰減系數(shù)可以表征材料的顯微結構[14].

1.2.1 聲 速

采用超聲脈沖回波法,根據(jù)被測材料厚度d和超聲波在材料中傳播時的回波時間間隔t,即可準確測量超聲波在被測材料中的傳播速度

1.2.2 聲阻抗

聲阻抗z的計算

式中:ρ為被測材料密度;v為聲速.

1.2.3 衰減系數(shù)

超聲波的衰減主要包括:被測材料對超聲波的衰減、耦合劑對超聲波的衰減以及超聲波擴散引起的衰減[15-16].超聲波探頭通過耦合劑和被測材料直接耦合在一起,屬于近場區(qū)域測量,因此超聲波的擴散損失忽略.耦合劑的使用對超聲波的傳播產(chǎn)生一定的能量損失,在計算時應考慮在內(nèi).因此,衰減系數(shù)

式中:A1為被測材料一次反射回波的幅值;A2為被測材料二次反射回波的幅值;α(z c,z o)為耦合劑對超聲波的衰減系數(shù);z o為耦合劑的聲阻抗;z c為被測材料的聲阻抗;d為被測材料厚度.

2 檢測系統(tǒng)

檢測系統(tǒng)原理框圖如圖2所示.檢測系統(tǒng)采用超聲脈沖回波法,超聲波探傷儀通過激發(fā)探頭激發(fā)超聲波,超聲波經(jīng)耦合劑的耦合作用到達被測試件的前表面,一部分由前表面反射回接收探頭,另一部分入射試件內(nèi)部,由后表面反射回接收探頭.接收探頭將接收到的反射回波信號通過數(shù)據(jù)采集卡進行高速采樣,并將采樣數(shù)據(jù)傳入計算機進行存儲和處理,處理結果可在顯示器上實時顯示.

圖2 檢測系統(tǒng)原理框圖F ig.2 The schematic diagram of inspection system

圖2中,界面1為耦合劑/被測試件界面,耦合劑選用機油,機油具有良好的耦合性能,保證了被測試件的前表面具有良好的耦合狀況;界面2為被測試件/空氣界面,由于空氣對超聲波有強烈的衰減作用,因此保證了后表面反射回波信號有足夠的強度返回到接收探頭.

3 實驗結果及分析

檢測系統(tǒng)中的超聲波探傷儀選用常州新區(qū)大平超聲波儀器有限公司生產(chǎn)的CUT-2007型超聲波探傷儀,選用超聲縱波探頭激發(fā)和接收超聲波,探頭工作頻率為2.5 MHz,超聲波垂直于被測試件表面入射.被測試件1選用編織玻璃纖維環(huán)氧樹脂基復合材料層壓板,尺寸大小為220mm×215 mm×7 mm;被測試件2選用單向碳纖維環(huán)氧樹脂基復合材料層壓板,尺寸大小為300mm×300mm×3mm.

圖3(a)和圖3(b)分別為被測試件1和被測試件2的實物圖.

圖3 被測試件實物圖Fig.3 Thesubstan ce d iagram of specim en

圖4和圖5分別表示被測試件1和被測試件2中實測超聲回波信號的Hilbert變換分析過程.圖4(a)和圖5(a)為實測超聲回波信號,圖4(b)和圖5(b)為去除噪聲后的超聲回波信號,圖4(c)和圖5(c)為Hilbert變換后的超聲回波包絡信號,圖4(d)和圖5(d)為歸一化的超聲回波包絡信號.圖中,S代表激發(fā)的超聲信號,R1代表一次反射回波信號,R2代表二次反射回波信號.

圖4 被測試件 1實測超聲回波信號及其Hilbert分析Fig.4 The measured ultrasonic echo signal w ith its hilber t transform of specimen 1

圖5 被測試件2實測超聲回波信號及其 Hilbert分析Fig.5 Them easured ultrason ic echo signal w ith its hilbert transform of specim en 2

由圖4和圖5可知:Hilbert變換方法能夠以較小的誤差跟蹤實測超聲回波信號的包絡,在準確跟蹤包絡信號的基礎上可進一步計算被測試件的聲學性能參數(shù).

經(jīng)過對被測試件1和被測試件2不同位置的多次測量,得到被測試件1的聲學性能參數(shù),如表1所示,被測試件2的聲學性能參數(shù)如表2所示.

對表1和表2的實驗數(shù)據(jù)分析可知:基于Hilbert變換的超聲檢測方法能夠準確測量復合材料的聲學性能,誤差控制在2.5%以內(nèi).此外,超聲波在選用的兩種復合材料中的傳播速度相近,但聲阻抗和衰減系數(shù)不同,這是由于在不考慮超聲擴散損失和使用相同耦合劑的情況下,被測試件本身的性能和組織結構對超聲波的衰減起主要作用,試件1使用玻璃纖維和環(huán)氧樹脂,纖維編織排列,試件 2使用碳纖維和環(huán)氧樹脂,纖維單向排列,不同的纖維種類、纖維體積含量以及纖維排列方式造成了聲阻抗和衰減系數(shù)的差異.

表1 被測試件1聲學性能參數(shù)(試件厚度7 mm)Tab.1 Acoustic property of specimen 1(specimen thickness:7 mm)

表2 被測試件 2聲學性能參數(shù)(試件厚度3 mm)Tab.2 Acoustic property of specimen 2(specimen thickness:3 mm)

4 結 論

1)從原理上介紹了Hilbert變換和復合材料聲學性能的測量,將Hilbert變換引入復合材料聲學性能的測量;

2)采用超聲脈沖回波法對兩種不同類型的復合材料層壓板進行檢測,利用Hilbert變換對檢測到的超聲回波信號進行分析和處理;

3)結果表明:本方法能夠準確測量復合材料聲學性能,測量誤差小于2.5%,同時,能夠有效評估復合材料本身性能.

[1] Birt EA,Smith R A.A review of NDEmethods for porositymeasurement in fibre-reinforced polymercomposites[J].Non-Destructive Testing and Condition Monitoring,2004,46(11):681-686.

[2] 吳瑞明,周曉軍,王邦凱.復合材料超聲C掃描多傳感器信息融合技術[J].傳感技術學報,2005,18(3):561-564.

Wu Ruim ing,Zhou Xiaojun,Wang Bangkai.Study onmultisensor information fusion in ultrasonic C scan system of com pos-itematerial[J].Chinese Journal of Sensors and A ctuators,2005,18(3):561-564.(in Chinese)

[3] 段劍,王從科,于波,等.纖維增強復合材料聲學性能檢測方法研究[J].測試技術學報,2004,18(1):39-42.

Duan Jian,Wang Congke,Yu Bo,et al.Research on testmethods of acoustics performance about resin-based composites[J].Journalof Testand Measurement Technology,2004,18(1):39-42.(in Chinese)

[4] 邢婭浪,徐章遂,陳鵬,等.超聲波技術在現(xiàn)代戰(zhàn)機復合材料結構損傷缺陷檢測中的應用[J].測試技術學報,2004,18(3):221-224.

Xing Yalang,Xu Zhangsui,Chen Peng,etal.App lication of ultrasonic detection in detecting flaw in compositematerial of modern battleplan[J].Journalof Test and Measurement Technology,2004,18(3):221-224.(in Chinese)

[5] Reyno ldsW N.Nondestructive testing of fibre-reinforced compositematerials[J].Materials&Design,1985,5(6):256-270.

[6] Bayraktar E,Antolovich SD,Bathias C.New developments in non-destructive controls of the compositematerials and applications in manufacturing engineering[J].Journal of Materials Processing Technology,2008,206:30-44.

[7] 李智,陳祥初,劉政波.包絡分析及其在設備故障診斷中的應用[J].測試技術學報,2002,16(2):92-95.

Li Zhi,Chen Xiangchu,Liu Zhengbo.Envelopment analysis and its application in the fau lt diagnosis[J].Journal of Test and Measurement Technology,2002,16(2):92-95.(in Chinese)

[8] 丁康,胡曉斌,謝明.解調(diào)分析的一種優(yōu)化方法[J].振動工程學報,1998,11(2):235-240.

Ding Kang,Hu Xiaobin,XieMing.An optim izedmethod of demodulation[J].Journal of Vibration Engineering,1998,11(2):235-240.(in Chinese)

[9] Soum ik Sarkar,KushalMukherjee,Asok Ray.Generalization of H ilbert transform for symbolic analysisof noisy signals[J].Signal Processing,2009,89:1245-1251.

[10] 舒泓,王毅.基于數(shù)學形態(tài)濾波和Hilbert變化的電壓閃變測量[J].中國電機工程學報,2008,28(1):111-114.

Shu Hong,W ang Yi.Voltage flicker detectionmethod based onmathematicalmorphology filterand Hilbert transform[J].Proceedings of the CSEE,2008,28(1):111-114.(in Chinese)

[11] 吳斌,鄧菲,何存富.超聲導波無損檢測中的信號處理研究進展[J].北京工業(yè)大學學報,2007,33(4):342-348.

W u Bin,Deng Fei,He Cunfu.Review of signalp rocessing in ultrasonic guided w avesnondestructive testing[J].Journal of Beijing University of Technology,2007,33(4):342-348.(in Chinese)

[12] 李瓏.超聲無損檢測和無損評價的一種新方法研究[J].聲學學報,1993,18(6):409-414.

Li Long.A new method for u ltrasonic NDT and NDE[J].Acta A custica,1993,18(6):409-414.(in Chinese)

[13] IvanW Selesnick.Hilbert Transform Pairs of W avelet Bases[J].Signal Processing Letters,2001,8(6):170-173.

[14] Graciet C,Hosten B.Simultaneousmeasurement of speed,attenuation,thickness and density with reflected ultrasonic waves in plates[C].Ultrasonic Symposium,1994,2:1219-1222.

[15] 王浩全,韓炎,曾光宇.非金屬復合材料超聲波衰減系數(shù)測定研究[J].測試技術學報,2002,16(4):249-251.

W ang Haoquan,Han Yan,ZengGuangyu.The research onmeasuring theu ltrasonic attenuation coefficientaboutnonmetal compositematerials[J].Journalof Test and Measurement Technology,2002,16(4):249-251.(in Chinese)

[16] Biwa S,Idekoba S,Ohno N.Wave attenuation in particu late polymer com posites[J].Mechanicso fMaterials,2002,34:671-682.