基于空間域諧波控制的電磁超聲Lamb波模態(tài)抑制方法

翟國富 李永虔 劉玥怡 秦伊人

基于空間域諧波控制的電磁超聲Lamb波模態(tài)抑制方法

翟國富 李永虔 劉玥怡 秦伊人

（哈爾濱工業(yè)大學電器與電子可靠性研究所 哈爾濱 150001）

電磁超聲換能器（EMAT）激勵的Lamb波中經(jīng)常存在多種模式，給波形分析造成困難，甚至會干擾對缺陷信息的判斷，從而導致誤判。常規(guī)的EMAT線圈空間波數(shù)譜中含有諧波，諧波是造成多模式現(xiàn)象的重要原因。該文提出試件渦流密度分布（洛倫茲力分布）相當于對線圈電流密度分布的空間低通濾波，通過對線圈空間域進行有效設計從而抑制空間域諧波，進而實現(xiàn)激發(fā)單波數(shù)正弦形式的渦流密度分布。結(jié)合正弦脈沖寬度調(diào)制（SPWM）技術(shù)的特點，提出將時域中的SPWM技術(shù)應用于空間域以實現(xiàn)線圈參數(shù)設計。通過建立有限元仿真模型，仿真證明了正弦形式的渦流密度空間分布可以被激勵。該文以激勵S0模態(tài)Lamb波為例進行實驗，實驗結(jié)果表明，由波數(shù)譜中的3次諧波激勵的A0模態(tài)被極大地抑制，抑制比為21dB。

電磁超聲換能器 Lamb波 模態(tài)抑制 空間域諧波控制 正弦脈沖寬度調(diào)制

0 引言

電磁超聲換能器（Electromagnetic Acoustic Transducers, EMAT）通過電磁耦合的方式產(chǎn)生超聲波[1-3]，由于無需耦合劑的優(yōu)點，使得其在無損檢測和無損評估中得到廣泛研究和應用。

EMAT是產(chǎn)生超聲Lamb波的一種有效方式。Lamb波的多模式現(xiàn)象會給波形分析造成一定的困難，甚至會干擾對缺陷信息的判斷，導致誤判。因此，Lamb波模態(tài)抑制方法一直是研究的熱點。劉素貞等根據(jù)對稱模式和反對稱模式兩種模式波結(jié)構(gòu)的特點，利用雙EMAT實現(xiàn)了單一模態(tài)Lamb波的激發(fā)和接收；該團隊還研究了電磁超聲換能器激發(fā)方式對Lamb波模態(tài)的影響，并通過對加載方式的選擇實現(xiàn)了Lamb波模態(tài)控制[4-5]。M. Seher等通過有限元仿真的方式，利用遺傳算法優(yōu)化了磁鐵直徑和磁鐵提離，得到了較為純凈的A0模式，S0模式的抑制比達30dB[6]。

P. Khalili等研究了在中、高頻厚積Lamb波模式下，壓電換能器和曲折線圈EMAT對單一波模式激勵的相對能力，說明了EMAT激發(fā)單一模式的導波需要同時考慮頻率和波數(shù)[7]。李明亮等建立了電磁超聲Lamb波位移和線圈參數(shù)的關(guān)系模型[8]，根據(jù)模型可以得出，曲折線圈激發(fā)的Lamb波的聲波位移幅值與表面洛倫茲力的波數(shù)譜幅值成正比。在此基礎上，Li Mingliang和Sun Wenxiu等在線圈參數(shù)固定的情況下，通過選擇合適的工作點實現(xiàn)了模態(tài)抑制和單一模態(tài)Lamb波的激發(fā)[9-10]。Zhai Guofu等將曲折線圈等效為空間橫向濾波器，設計了可抑制諧波的濾波函數(shù)，并根據(jù)濾波函數(shù)設計了一種特殊參數(shù)的曲折線圈，實現(xiàn)了多模態(tài)現(xiàn)象的抑制[11]。

研究者們主要通過改變工作點、時域信號激勵或EMAT結(jié)構(gòu)等方式實現(xiàn)Lamb波模態(tài)抑制，設計方法大多較為復雜，實際使用時較為不便或適用范圍小。

EMAT激勵的Lamb波模態(tài)和試件表面洛倫茲力的空間波數(shù)譜直接相關(guān)。因此，可以通過控制空間域中的波數(shù)譜的方式實現(xiàn)Lamb波模態(tài)抑制。常規(guī)EMAT線圈波數(shù)譜中，除了基波外，還存在奇次諧波[11]，這是造成多模式現(xiàn)象的重要原因。因此，若能重新對EMAT線圈進行有效設計，使得其激發(fā)的洛倫茲力在空間域中僅含基波（也就是僅有單一波數(shù)），則可有效抑制因諧波而產(chǎn)生的多模態(tài)現(xiàn)象。

本文將對EMAT模型進行分析，提出新的EMAT線圈設計方法以抑制空間域諧波，最終實現(xiàn)Lamb波模態(tài)抑制。

1 EMAT模型分析

圖1 電磁超聲換能器結(jié)構(gòu)

假設曲折線圈的電流密度空間分布函數(shù)為()。對空間周期為2的線圈在空間上進行延拓，得到無限大的周期性電流密度分布。此時，函數(shù)() 可進行傅里葉級數(shù)展開，有

試件中被感應的渦流密度e大小[8, 12]為

2 線圈設計

2.1 線圈設計要求

第1節(jié)的分析中，線圈在空間域上是無限大的，實際的線圈電流密度和試件渦流密度空間分布函數(shù)應該是被窗函數(shù)調(diào)制的。傳統(tǒng)曲折線圈對應的窗函數(shù)為矩形窗。矩形窗并不會改變各諧波的次數(shù)和幅值，僅會產(chǎn)生主瓣和旁瓣。

由第1節(jié)分析可得，若想抑制因諧波而產(chǎn)生的Lamb波模態(tài)，渦流密度分布的波數(shù)譜中不應包含諧波。不包含諧波的渦流密度分布應該為單波數(shù)正弦波，其表達式為

由于試件渦流密度分布是對線圈電流密度分布的空間低通濾波，且諧波次數(shù)越高，濾波效果越好。因此，為了能抑制渦流密度分布的空間域諧波，結(jié)合EMAT線圈的特點，提出以下線圈設計要求： ①線圈空間域中必須包含基波；②線圈空間域中可以包含諧波，諧波幅值盡可能低，諧波次數(shù)盡可能大；③線圈各導線的電流密度相等或電流相等。

2.2 線圈設計方法

正弦脈沖寬度調(diào)制（Sinusoidal Pulse Width Modulation, SPWM）技術(shù)是控制領(lǐng)域的一項技術(shù)，只應用在時域中，并且要求系統(tǒng)是慣性環(huán)節(jié)[13]。EMAT的電流密度分布在空間域，而且試件渦流密度分布和線圈電流密度分布的關(guān)系并不是空間域中的慣性環(huán)節(jié)。但是，SPWM技術(shù)應用于電磁超聲領(lǐng)域有其獨特性：①SPWM脈沖序列幅值相等，可以實現(xiàn)線圈各導線的電流密度相等；②SPWM脈沖序列的幅度譜中，除基波外，還有一些高次諧波[14]，滿足2.1節(jié)提出的線圈設計要求。SPWM技術(shù)中，載波比越大，則諧波次數(shù)越高。幅值比c/m＞1時，幅值比越小，基波分量越大，諧波相對基波的幅值越低[14]。因此，只要幅值比大于1，并且盡可能小，載波比盡可能大，則滿足SPWM脈沖序列分布的線圈可以激發(fā)出近似正弦的渦流密度空間分布。因此，綜合考慮SPWM技術(shù)的特點和電磁超聲的獨特性，時域中的SPWM技術(shù)可以應用于空間域中的EMAT線圈設計，且能很好地符合本文的線圈設計要求。

SPWM基本原理如圖2所示。本文中，輸出信號就是最終的電流密度分布。輸出信號中每個的1或-1的寬度對應著線圈的寬度（正負號代表電流流向，電流正方向可自定義），每個0的寬度對應著導線間距。使用此方法可以得到線圈的參數(shù)。

圖2 SPWM基本原理

2.3 線圈參數(shù)設計

圖3 線圈參數(shù)和實際繪制的PCB線圈

表1 線圈具體參數(shù)

Tab.1 Coil parameter value

3 仿真和實驗

3.1 有限元仿真

本文首先使用COMSOL軟件建立二維頻域電磁場仿真模型，驗證SPWM線圈的有效性。建立如圖4所示的有限元法（Finite Element Analysis, FEA）仿真模型，模型參數(shù)見表2。板材選用鋁板，換能器處于鋁板中心位置。EMAT中的線圈總寬度為4個空間周期。仿真模型中線圈激勵方式設置為電壓激勵，電壓為100V。為了保證電流密度一致，每根導線設置為多匝線圈。

圖4 有限元法仿真模型

表2 模型參數(shù)

Tab.2 Model parameters

將SPWM線圈參數(shù)代入仿真模型中，提取換能器下方鋁板表面沿軸方向的洛倫茲力和渦流密度如圖5a所示。計算圖5a中洛倫茲力、渦流密度的波數(shù)譜和SPWM線圈電流密度的波數(shù)譜如圖5b所示。

由圖5a可知，洛倫茲力分布幾乎呈正弦形式，但是存在幅值很小的高次諧波，高次諧波的波數(shù)和載波的波數(shù)接近。對比渦流密度和洛倫茲力可以發(fā)現(xiàn)，洛倫茲力的邊緣幅值要大于中間部分的幅值，這是因為磁鐵邊緣的磁感應強度要大于磁鐵中心的磁感應強度。

由圖5b可知，在低波數(shù)部分，洛倫茲力波數(shù)譜的旁瓣幅值要高于SPWM線圈電流密度波數(shù)譜的旁瓣幅值，而渦流密度波數(shù)譜的旁瓣幅值和線圈電流密度波數(shù)譜的旁瓣幅值幾乎一致。因此可以確定，洛倫茲力分布的高旁瓣幅值是由于磁場分布不均勻造成的。而在高波數(shù)部分，洛倫茲力和渦流密度波數(shù)譜幾乎一致，諧波次數(shù)和線圈電流密度波數(shù)譜的諧波次數(shù)一致，但是幅值要低很多。說明實際的渦流密度分布會對高次諧波產(chǎn)生抑制效果，與理論分析一致。

圖5 洛倫茲力、渦流密度、線圈電流密度及其波數(shù)譜

由于導線內(nèi)的趨膚效應會造成各導線內(nèi)的電流密度并不嚴格一致，為了明確線圈導線內(nèi)的趨膚效應的影響，將仿真模型中各導線設置改為單導線。單導線激勵下的渦流密度分布的波數(shù)譜如圖6所示。由圖可知，導線內(nèi)的趨膚效應會導致渦流密度分布中存在3次諧波，幅值為基波幅值的5.5%。由于此幅值很小，并且接收換能器有更進一步的空間濾波作用，因此，可以認為導線內(nèi)趨膚效應造成的影響可以忽略。

圖6 考慮導線內(nèi)的趨膚效應時，渦流密度分布的波數(shù)譜

3.2 實驗

圖7 Lamb波頻散曲線和工作點

搭建實驗裝置如圖8所示。EMAT的激勵和信號處理設備為零聲科技（蘇州）有限公司研制的EMAT-2000便攜式電磁超聲探傷儀，主要技術(shù)參數(shù)見表3，參數(shù)完全滿足實驗要求。被測試件為1mm厚的6061鋁板。換能器為一發(fā)一收，發(fā)射接收EMAT的參數(shù)完全一致，中心間距為=200mm。EMAT中的磁鐵為50mm×50mm×25mm的N35型釹鐵硼磁鐵。實驗使用SPWM線圈和常規(guī)線圈作為對比，SPWM線圈參數(shù)如圖3所示，常規(guī)線圈參數(shù)為：=0.5mm，=6mm，線圈匝數(shù)為8。實驗中探傷儀和EMAT之間未添加阻抗匹配電路，這在一定程度上會降低信噪比，但不會影響實驗結(jié)果。

圖8 實驗裝置

表3 探傷儀參數(shù)

Tab.3 Instrument parameters

設置探傷儀激勵參數(shù)為周波數(shù)4，發(fā)射強度50%（約800pp），激勵頻率450kHz。設置接收參數(shù)為采樣率20MHz，采樣深度8k，平均次數(shù)64。首先使用常規(guī)線圈進行實驗，接收增益為83dB，得到的實驗結(jié)果如圖9a所示。

圖9 實驗結(jié)果

根據(jù)Lamb波群速度曲線可以計算，S0模態(tài)的群速度約為5 420m/s，A0模態(tài)的群速度約為2 880m/s。S0模態(tài)的理論到達時間約為29.2ms，A0模態(tài)的理論到達時間約為54.9ms。圖9a中第一個波包的實際到達時間為31.4ms，第二個波包的實際到達時間為56.9ms?？紤]到電路延時、計算誤差、測量誤差以及Lamb波成型需要一段距離等因素，可以判定圖中前兩個波包分別為S0模態(tài)和A0模態(tài)。A0模態(tài)與S0模態(tài)的幅值比約為0.4。

為進一步判斷S0模態(tài)和A0模態(tài)的波數(shù)，移動發(fā)射EMAT位置，在=200～300mm之間，每2mm進行一次測量，得到51組數(shù)據(jù)。對數(shù)據(jù)進行時間和空間二維傅里葉變換，得到頻率-波數(shù)如圖9b所示。由圖9b可以證實，S0模態(tài)和A0模態(tài)被激發(fā)。S0模態(tài)的中心波數(shù)約為0.52mm-1，而A0模態(tài)的中心波數(shù)約為1.53mm-1。A0模態(tài)的波數(shù)約為S0模態(tài)波數(shù)的3倍，與理論分析相一致，證明A0模態(tài)確實是由波數(shù)譜中的3次諧波激勵的。

使用SPWM線圈進行實驗，發(fā)射接收EMAT間距=200mm，接收增益為67dB（和常規(guī)線圈EMAT接收增益不一致的原因是線圈參數(shù)發(fā)生了改變，導致EMAT換能效率存在較大差異），結(jié)果如圖9c所示。A0模態(tài)得到了極大的抑制，幅值接近0。A0模態(tài)與S0模態(tài)的幅值比約為0.036。相對傳統(tǒng)線圈，抑制比達21dB。證明了該線圈有較好的Lamb波模態(tài)抑制效果。

對于仍存在的幅值很小的A0模態(tài)，本文認為，這并非是導線內(nèi)的趨膚效應引起的。由于Lamb波理論模型是基于穩(wěn)態(tài)情況下建立的，而Lamb波成型需要經(jīng)歷一個暫態(tài)過程。穩(wěn)態(tài)模型不能100%地還原這個暫態(tài)過程，導致了A0模態(tài)無法被完全消除，只能被抑制。

4 結(jié)論

諧波是造成多模式現(xiàn)象的重要原因，本文提出了基于空間域諧波控制的電磁超聲Lamb波模態(tài)抑制方法，并將時域中的SPWM技術(shù)應用于空間域線圈參數(shù)設計，仿真和實驗證明了該方法的有效性，模態(tài)抑制比達21dB。

本文分析了電磁超聲理論模型，指出試件中渦流密度分布和洛倫茲力分布一致，均是對EMAT線圈電流密度分布的空間低通濾波。渦流密度分布和電流密度分布的諧波次數(shù)和相位是一致的，但是隨著諧波次數(shù)升高、提離距離增大，渦流密度分布的諧波幅值相對電流密度分布的諧波幅值呈指數(shù)趨勢下降。

為了實現(xiàn)對渦流密度空間域諧波的控制，本文指出線圈需能夠激勵出單波數(shù)正弦形式的渦流密度分布，并提出了新的線圈設計要求：線圈空間域中必須存在基波，可以存在幅值盡可能低、次數(shù)盡可能高的諧波，且線圈各導線的電流密度或電流必須相等。結(jié)合SPWM技術(shù)的特點，指出SPWM技術(shù)符合本文的線圈設計要求，可以實現(xiàn)線圈參數(shù)設計。

被EMAT激勵的Lamb波模態(tài)是由線圈空間波數(shù)譜決定的，本文所述的基于空間域諧波控制的電磁超聲Lamb波模態(tài)抑制方法解決了波數(shù)譜中的諧波問題，從根本上解決了因諧波產(chǎn)生的多模態(tài)現(xiàn)象。該設計方法適用范圍廣，且簡單易行。

[1] 翟國富, 梁寶, 賈文斌, 等. 橫波電磁超聲相控陣換能器設計[J]. 電工技術(shù)學報, 2019, 34(7): 1441- 1448.

Zhai Guofu, Liang Bao, Jia Wenbin, et al. Design of the shear wave electromagnetic ultrasonic phased array transducer[J]. Transactions of China Electro- technical Society, 2019, 34(7): 1441-1448.

[2] 孫文秀, 劉國強, 夏慧, 等. 非鐵磁材料的電磁超聲接收過程數(shù)值模擬及實驗研究[J]. 電工技術(shù)學報, 2018, 33(19): 4443-4449.

Sun Wenxiu, Liu Guoqiang, Xia Hui, et al. Numerical simulation and experimental study of electro- magnetic acoustic transducer receiving process in nonferromagnetic material[J]. Transactions of China Electrotechnical Society, 2018, 33(19): 4443-4449.

[3] 劉素貞, 王淑娟, 張闖, 等. 鋼板電磁超聲表面波的仿真分析及缺陷定量檢測[J]. 電工技術(shù)學報, 2020, 35(1): 97-105.

Liu Suzhen, Wang Shujuan, Zhang Chuang, et al. Simulation analysis of electromagnetic acoustic surface wave of steel plate and quantitative defect detection[J]. Transactions of China Electrotechnical Society, 2020, 35(1): 97-105.

[4] 劉素貞, 張嚴偉, 張闖, 等. 基于電磁超聲雙換能器的單模態(tài)Lamb波激勵方法研究[J]. 電工技術(shù)學報, 2018, 33(1): 140-147.

Liu Suzhen, Zhang Yanwei, Zhang Chuang, et al. Research on excitation method of single-mode Lamb wave based on electromagnetic acoustic double transducer[J]. Transactions of China Electrotechnical Society, 2018, 33(1): 140-147.

[5] 劉素貞, 董碩, 張闖, 等. 電磁超聲加載方式對Lamb波模態(tài)的影響[J]. 電工技術(shù)學報, 2018, 33(19): 4426-4433.

Liu Suzhen, Dong Shuo, Zhang Chuang, et al. The influence of electromagnetic acoustic transducer loading way on the mode of Lamb wave[J]. Transa- ctions of China Electrotechnical Society, 2018, 33(19): 4426-4433.

[6] Seher M, Huthwaite P, Lowe M J S, et al. Model- based design of low frequency Lamb wave EMATs for mode selectivity[J]. Journal of Nondestructive Evaluation, 2015, 34(3): 22.

[7] Khalili P, Cawley P. Relative ability of wedge coupled piezoelectric and meander coil EMAT probes to generate single mode Lamb waves[J]. IEEE Transa- ctions on Ultrasonics Ferroelectrics & Frequency Control, 2018, 64(4): 648-656.

[8] 李明亮, 鄧明晰, 高廣健. 電磁超聲換能器激發(fā)超聲蘭姆波的模式展開分析[J]. 聲學學報, 2014, 39(3): 385-391.

Li Mingliang, Deng Mingxi, Gao Guangjian. Modal expansion analyses of ultrasonic Lamb waves generated by electromagnetic ultrasonic trans- ducers[J]. Acta Acustica, 2014, 39(3): 385-391.

[9] Li Mingliang, Deng Mingxi, Gao Guangjian. Selective generation of ultrasonic Lamb waves by electromagnetic acoustic transducers[J]. Chinese Physics B, 2016, 25(12): 124301.

[10] Sun Wenxiu, Liu Guoqiang, Xia Hui, et al. Lamb wave signal selective enhancement by an improved design of meander-coil electromagnetic acoustic trans- ducer[J]. Chinese Physics B, 2018, 27(8): 084301.

[11] Zhai Guofu, Jiang Tao, Kang Lei. Analysis of multiple wavelengths of Lamb waves generated by meander-line coil EMATs[J]. Ultrasonics, 2014, 54(2): 632-636.

[12] Hirao M, Ogi H. Electromagnetic acoustic trans- ducers: noncontacting Ultrasonic measurements using EMATs[M]. 2nd. Japan: Springer, 2018.

[13] 王兆安, 劉進軍. 電力電子技術(shù)[M]. 5版. 北京: 機械工業(yè)出版社, 2013.

[14] 邵江富. 逆變器脈寬調(diào)制輸出電壓的諧波分析[J]. 機車電傳動, 1980(4): 26-32.

Mode Suppression Method of Lamb Wave Excited by Electromagnetic Acoustic Transducers Based on Spatial Harmonic Control

（Institute of Reliability in Electrical Apparatus and Electronics Harbin Institute of Technology Harbin 150001 China）

Lamb waves excited by electromagnetic acoustic transducers (EMAT) often have multiple modes, which make waveform analysis difficult and even interfere with the judgment of defect information, leading to misjudgment. Harmonics of conventional EMAT coils in the space wavenumber spectrum are the important cause of the multi-mode phenomena. This paper proposes that the eddy current density distribution (or the Lorentz force distribution) in the specimen is equivalent to the spatial low-pass filtering of the coil current density distribution. By effectively designing the coil, the harmonics in the space domain can be suppressed, and the sinusoidal (that is, single wavenumber) eddy current density distribution can be excited. Combining the characteristics of sinusoidal pulse width modulation (SPWM) technology, it was proposed to apply the time domain SPWM technology to the space domain to design the coil parameters. A finite element simulation model was established, and the simulation proved that the sinusoidal eddy current density distribution can be excited. The experiment was conducted by taking the S0 mode Lamb wave as an example. It is shown that the A0 mode excited by the third harmonic in the wavenumber spectrum is greatly suppressed, and the suppression ratio is 21dB.

Electromagnetic acoustic transducers (EMAT), Lamb, mode suppression, spatial harmonic control, sinusoidal pulse width modulation (SPWM)

10.19595/j.cnki.1000-6753.tces.200467

TM15

翟國富 男，1964年生，教授，博士生導師，研究方向為高可靠電器設計理論與技術(shù)、產(chǎn)品質(zhì)量一致性穩(wěn)健設計理論與技術(shù)、電子系統(tǒng)及裝備可靠性和可測性設計及健康管理、電磁超聲無損檢測技術(shù)。E-mail: gfzhai@hit.edu.cn

李永虔 男，1991年生，博士研究生，研究方向為電磁超聲無損檢測技術(shù)。E-mail: 15B906029@hit.edu.cn（通信作者）

國家自然科學基金（51977044）和民用航天技術(shù)預先研究（JZKGJ20190004）資助項目。

2020-05-07

2020-06-16

（編輯 崔文靜）