利用信號反向傳播法移除Lamb頻散特性

1、利用頻散補償法消除蘭姆波頻散效應倪園，程建政，張德俊，周玉祿（中科院武漢物理與數學研究所 武漢 430071）摘 要 蘭姆波的多模式及頻散特性使得其在激勵、傳播以及信號處理等方面較為復雜。本文在對蘭姆波頻散特性進行了分析的基礎上，用頻散補償方法對在鋼板中傳播10cm距離的蘭姆波信號作模擬計算，所得信號明顯地減小了因頻散而引起脈沖信號拉長。在實驗中，對厚度為3mm鋁板中傳播了10cm的頻散信號進行處理，也得到較好的結果。關鍵詞 蘭姆波；頻散補償；群速度；相速度Using Dispersion Compensation Method to Remove the Dispersion of Lamb

2、 Wave NI yuan, CHENG jianzheng, ZHANG dejun, ZHOU yulu（Wuhan Institute of Physics and Mathematics the Chinese Academy of Sciences Wuhan 430071）Abstract Multi-modes and dispersion of Lamb wave make it very complicated to be used in NDE. This paper analyzes the dispersion of Lamb wave, and proposes a

3、dispersion compensation method to remove the dispersion of the received signal. Using this method, S0 mode of Lamb wave propagating 10cm in steel plate has been simulated, then attained a signal with high spatial resolution. In the experiment, the signal propagating 10cm in the 3mm thickness Al plat

4、e is processed, and achieve a good result.Keywords Lamb wave; Dispersion compensation; Group velocity; Phase velocity1 引言隨著現代工業(yè)的發(fā)展，對板狀構件內部的缺陷進行檢測、定征的需求越來越迫切，特別是近年來發(fā)展起來的復合材料以其密度小、強度高、耐高溫等優(yōu)點而被認為是航空航天工業(yè)的理想新材料。用體聲波檢測板狀材料時，盲區(qū)問題是一個難以克服的困難，而蘭姆波自然成為在線檢測板結構件的有力手段1。蘭姆波本質上是二維的，衰減較三維體波小，所以傳播距離長，工程技術人員因此希望應用它對大型

5、的板狀材料進行快速的無損檢測。但是由于蘭姆波的頻散特性使得其在激勵、傳播、接收以及信號處理等方面比較復雜，這使得蘭姆波在無損檢測技術領域里并未得到充分應用。本文在對蘭姆波的頻散特性作定性分析的基礎上，針對傳播了一段距離的S0模式蘭姆波進行了模擬計算，經頻散補償后，將信號從時域變換到空間域，從而精確地表示出蘭姆波包絡傳播的距離及缺陷的位置。通過實驗，對在鋁板中傳播了10cm距離的A0模式蘭姆波進行了處理，得到較好結果。該研究對實現蘭姆波在工業(yè)中的應用具有重要意義。 2 蘭姆波頻散分析蘭姆波是在一種厚度與激勵聲波波長為相同數量級的聲波導中由縱波和橫波合成的特殊形式的應力波，它是20世紀初期英國力學

6、家蘭姆解平板自由邊界條件下波動方程而發(fā)現的。按照板內振動位移的分布形態(tài)不同，可分成對稱模式與反對稱模式，它們分別有一系列不同頻散特性的各階模式。在厚度為2b的無限大固體自由平板中，蘭姆波的頻散方程2為：對稱模式： (1a) 反對稱模式： (1b) ，式中k0為沿板水平方向的波數，b為1/2板厚，=2f為角頻率，cl為縱波速度（常數），cs為橫波速度(常數）。由于k0與的關系不是線性的，因而蘭姆波的相速度不是常數，而是隨頻率變化的。根據式(1a)(1b)，并參照文獻3的方法，我們運用matlab軟件，計算出鋼（cl=5790m/s, cs =3200m/s,）的相速度頻散曲線如圖1所示.。 圖1

7、鋼的相速度曲線 圖2 鋼的群速度曲線由于頻散效應，蘭姆波脈沖信號在板中傳播包絡將逐漸拉長，其傳播速度為群速度cg。cg定義為，據此容易計算出鋼的群速度頻散曲線如圖2所示。3 蘭姆波頻散補償在一個適當的薄板材結構上，用發(fā)射換能器激發(fā)某種特定的蘭姆波模式。發(fā)射換能器所在的位置設定為原點，換能器的輸出信號為f(t),所激發(fā)的聲場用u(x,t)來表示，t為傳播時間，x為傳播距離，可推導出板中任意一點的聲信號為： （2）k()為角頻率時的波數，F()是f(t)的傅里葉變換。信號f(t)傳播一段距離后，遇到缺陷反射回來的信號設為g(t)。g(t)通常由不同缺陷反射信號的疊加組成，設各反射信號的傳播距離分別

8、為dj，反射處的反射系數為Aj()，則g(t)可表示為： （3）式(3)中，假設在缺陷反射處不存在模式轉換,且只考慮一個缺陷的情況，即j=1，那么g(t)可視為接收換能器在距離發(fā)射換能器為兩倍缺陷距離處接收到的信號，如圖3所示。 圖3 反向傳播示意圖通過下式： x=Vgr t （4）將信號從時域轉換到空域。Vgr為F()中心頻率所對應的群速度。由式(4)得到的空域信號將在空間上拉得較長，不能精確地計算發(fā)射信號的傳播距離，因此需要采用頻散補償方法對g(t)作處理。 （5）式(5)在頻域對接收信號在每個頻率成份上進行補償。其中G()為g(t)的傅里葉變換，k()為角頻率對應的波數，可由得到。由群速

9、度計算公式，代入式(5)得到： （6）其中。由式(6)可計算出頻散補償后的空域信號。該式為反傅里葉變換式3，易于數值求解。式(5)到式(6)的轉換是頻散補償的關鍵，需以相速度和群速度頻散曲線為基礎。4 模擬計算 采用鋼板進行模擬計算。圖4a為發(fā)射換能器輸出信號f(t)，其中心頻率為2MHz，帶寬為1MHz。假設所激勵起的只有S0一種模式，在距離發(fā)射換能器10cm處用另一個換能器接收，不考慮換能器性能的影響，根據式(2)可模擬計算出頻散信號g(t)如圖4b。若不經頻散補償處理，由式(4)，可將g(t)轉換為空域信號h(x)，如圖4c。顯見未經頻散補償處理得到的信號空間持續(xù)相當長，無法確定傳播的距

10、離。而由式(6)經頻散補償處理后得到的信號如圖4d，可以看出信號壓縮效果較好，能精確的表明傳播距離。 圖4 模擬仿真示意圖模擬仿真步驟為：（1） 將頻散信號g(t)至少補8倍以上的零點，以減小頻域的取樣間隔，避免信號在頻域失真。設補零后的信號總長度為n，間隔為，將補零后的g(t)作FFT變換，得到G()，頻域間隔為。（2） 計算以為變量的板材的相速度cp和群速度cg，插值使其與G()的頻率范圍及步長一致。（3） 計算。用式，將轉換為k。但此時所得到的k是不等間距的。（4） 將H(k)以進行插值，得到新的H(k)，利用式6即可得到h(x)。（5） 設信號h(x)為m個點，步長為，因頻散補償后的信

11、號總長度不應小于初始的發(fā)射信號，應滿足式,而=1/（m），可知<1/（）。當時基信號轉換到空域信號時，空間奈奎斯特波數不應小于奈奎斯特頻率處的波數，由此可得出：，而fNyq=1/（2），=1/（2），得到。此過程可得出、的最大值，所以m的取值應。5 實驗研究 51 實驗裝置 實驗裝置如圖5所示，發(fā)射與接收換能器均為2.5MHz寬帶超聲縱波直探頭，采用猝發(fā)音發(fā)生器激勵發(fā)射換能器產生脈沖超聲波4。panametrics5052UA型超聲波分析儀將接收信號放大。被測試樣為3mm厚的鋁板，兩換能器表面均涂上少量耦合劑，直接緊壓在鋁板的同一側，超聲波分析儀將接收信號放大后送入TDS1012數字示波

12、器，通過示波器數據存儲卡將數據讀入PC機進行處理。 圖5 實驗裝置示意圖52 實驗結果實驗選用的材料是3mm厚的鋁板，將接收換能器置于距發(fā)射換能器10cm處。微調猝發(fā)音的發(fā)射頻率，發(fā)現在1.1MHz時，產生的蘭姆波模式較為單一，經波形數據分析可知，所記錄下來的信號為A0模式。由群速度及相速度曲線可知，S0模式應該在A0模式之前到達，但是由于S0模式信號太弱，無法提取有用信號進行分析，因此,實驗時針對A0模式進行分析處理。不經頻散補償，接收時域信號直接轉換成空間域信號，如圖6a所示，而頻散補償后的信號如圖6b所示。 圖6 實驗結果 從圖6a可以看出，未經補償的頻散信號在空間的脈沖寬度約為150m

13、m，而經頻散補償后的脈沖寬度約為50mm，壓縮效果較為明顯。但與模擬的結果相比，信號壓縮得還不夠理想。分析其原因可能是實驗中板材的參數不精確，即板材的實際參數與計算相速度、群速度曲線所用參數有偏差所致。另外，在實驗的過程中，因信噪比較低，存在噪聲干擾，這些都可能造成實驗結果的誤差。6結論本文對蘭姆波的頻散特性及蘭姆波在板中的傳播作了定性分析。在理論分析的基礎上，采用頻散補償方法，通過模擬和實驗對在板中傳播了一段距離后的頻散信號作處理，得到壓縮的空域信號，較好地提高了蘭姆波檢測的分辨率及缺陷定位精確，同時也驗證了該方法的有效性。此外，實驗研究可在提高信噪比方面進行改進，以得到更理想的結果，本文作者將就此開展進一步的研究。參考文獻1 杜功煥,等.聲學基礎M.南京:南京大學出版社,2001,537-540.2 應崇福,等.超聲在固體中的散射M.北京:國防工業(yè)出版社,1994,16-22.3 鄭祥明,趙玉珍,史耀武.蘭姆波頻散曲線的計算J無損檢測,2003,25(2)4 Paul D.Wilcox.A rapid signal processing technique to remove the effect of dispersion from guided wavw