何存富，劉青青，焦敬品，劉 飛，吳 斌

（北京工業(yè)大學(xué)機(jī)械工程與應(yīng)用電子技術(shù)學(xué)院，北京 100124）

其中：級(jí)數(shù)中的每一項(xiàng)代表了一個(gè)在空間具有固定模式并按一定頻率振動(dòng)的駐波?？梢钥闯?，振動(dòng)解利用

基于振動(dòng)模態(tài)分析的鋼軌中超聲導(dǎo)波傳播特性數(shù)值計(jì)算方法

何存富，劉青青，焦敬品，劉 飛，吳 斌

（北京工業(yè)大學(xué)機(jī)械工程與應(yīng)用電子技術(shù)學(xué)院，北京 100124）

針對(duì)鋼軌等復(fù)雜異型波導(dǎo)結(jié)構(gòu)中超聲導(dǎo)波頻散特性難于求解的問題，通過振動(dòng)模態(tài)分析提取出鋼軌中超聲導(dǎo)波傳播特性。基于振動(dòng)特征頻率和波動(dòng)解互相轉(zhuǎn)化的原理，將振動(dòng)模態(tài)分析獲得的特征頻率解，轉(zhuǎn)換成波動(dòng)解，進(jìn)而提取出鋼軌的頻散特性。同時(shí)，還可以由結(jié)構(gòu)振動(dòng)分析中獲得的變形信息，計(jì)算出超聲導(dǎo)波的波結(jié)構(gòu)。通過對(duì)計(jì)算出的鋼軌中低階超聲導(dǎo)波模態(tài)頻散特性和波結(jié)構(gòu)分析，選擇出適合鋼軌無損檢測(cè)的超聲導(dǎo)波模態(tài)類型和頻率范圍，為后續(xù)專用超聲導(dǎo)波傳感器設(shè)計(jì)和檢測(cè)實(shí)驗(yàn)提供理論支持。

模態(tài)分析；特征頻率法；鋼軌；超聲導(dǎo)波；頻散特性；波結(jié)構(gòu)

高速鐵路由于具有高速、大容量等優(yōu)勢(shì)，已經(jīng)成為現(xiàn)代社會(huì)最重要的交通運(yùn)輸方式之一。隨著我國高速鐵路運(yùn)營里程的跨越式發(fā)展以及客運(yùn)高速化、貨運(yùn)重載化的發(fā)展，對(duì)高速鐵路運(yùn)輸?shù)陌踩砸笕找嫣岣?。特別是“7·23”動(dòng)車追尾特大交通事故后，高速鐵路運(yùn)輸?shù)陌踩愿莻涫苋藗冴P(guān)注。

作為高速鐵路運(yùn)輸系統(tǒng)的重要組成部件，軌道具有引導(dǎo)機(jī)車前進(jìn)、承受機(jī)車荷載以及將其傳遞至軌枕上的功能。由于自身材質(zhì)缺陷以及長期受到各種外力及環(huán)境腐蝕綜合作用，在軌道上會(huì)產(chǎn)生局部損傷甚至宏觀斷裂，從而導(dǎo)致列車出軌、傾覆等重大行車安全事故，造成人員傷亡和財(cái)產(chǎn)損失。因此，需要利用各種無損檢測(cè)技術(shù)對(duì)鋼軌的損傷狀態(tài)進(jìn)行檢測(cè)。目前，常見的鋼軌檢測(cè)方法有超聲波法、軌道電路、應(yīng)變片檢測(cè)法、光纖斷軌檢測(cè)法等。這些常規(guī)方法還不能滿足現(xiàn)代軌道安全評(píng)價(jià)的需要，迫切需要發(fā)展高效、快捷的軌道無損檢測(cè)新技術(shù)。作為一種新型無損檢測(cè)技術(shù)，超聲導(dǎo)波由于具有傳播距離遠(yuǎn)、檢測(cè)效率高等優(yōu)點(diǎn)，對(duì)于鋼軌大范圍缺陷檢測(cè)具有特殊優(yōu)勢(shì)，因此，鋼軌超聲導(dǎo)波檢測(cè)技術(shù)也是目前國內(nèi)外無損檢測(cè)領(lǐng)域研究的熱點(diǎn)問題。

由于超聲導(dǎo)波的檢測(cè)效果很大程度上取決于所選超聲導(dǎo)波模態(tài)的傳播特性。因此，對(duì)結(jié)構(gòu)進(jìn)行超聲導(dǎo)波檢測(cè)時(shí)，首先需要根據(jù)結(jié)構(gòu)中超聲導(dǎo)波的傳播特性，如頻散、波結(jié)構(gòu)、衰減及可激勵(lì)性，選擇合適的超聲導(dǎo)波模態(tài)類型和頻率范圍。對(duì)于規(guī)則截面波導(dǎo)（如板、桿、管）中超聲導(dǎo)波的傳播特性，其數(shù)值計(jì)算方法已經(jīng)很成熟［1－5］，且已有商用的數(shù)值分析軟件［6］。但對(duì)于鋼軌等異型截面波導(dǎo)中超聲導(dǎo)波傳播特性的研究，商用軟件無法完成其分析計(jì)算。因此，異型截面波導(dǎo)中超聲導(dǎo)波傳播特性研究目前仍是超聲導(dǎo)波研究領(lǐng)域的難點(diǎn)和熱點(diǎn)問題。

國外很多學(xué)者［7－8］開展了異型截面波導(dǎo)中超聲導(dǎo)波傳播特性研究。例如，Rose等［9］利用半解析有限元的方法對(duì)任意截面波導(dǎo)中超聲導(dǎo)波的傳播問題進(jìn)行了研究，計(jì)算出100 kHz頻率范圍鋼軌中主要導(dǎo)波模態(tài)的相速度、群速度頻散曲線，并在軌頭踏面、軌頭側(cè)面進(jìn)行了超聲導(dǎo)波試驗(yàn)驗(yàn)證。同樣采用半解析有限元方法，Lee等［10］采用半解析有限元的方法對(duì)鋼軌中瑞利波及高頻蘭姆波傳播特性進(jìn)行了研究，并通過試驗(yàn)證明其對(duì)軌頭下方的橫向裂紋檢測(cè)具有很好的敏感性。Coccia等［11］采用半解析有限元對(duì)不同激勵(lì)條件下鋼軌軌頭中導(dǎo)波傳播特性進(jìn)行了研究，結(jié)果發(fā)現(xiàn)，對(duì)稱激勵(lì)時(shí)，能量主要集中在軌頭踏面；非對(duì)稱激勵(lì)時(shí)，能量集中于軌頭側(cè)面。Hesse等［12］采用半解析有限元法對(duì)鋼軌中表面波傳播特性進(jìn)行了理論分析。

綜上所述，對(duì)于任意截面波導(dǎo)（包括鋼軌）中超聲導(dǎo)波傳播特性的計(jì)算方法較單一，多采用半解析有限元計(jì)算頻散特性，然后結(jié)合其他方法仿真或?qū)嶒?yàn)進(jìn)行驗(yàn)證。針對(duì)以上研究現(xiàn)狀，基于結(jié)構(gòu)振動(dòng)解與波動(dòng)解的互化原理，本文利用振動(dòng)模態(tài)分析的方法進(jìn)行鋼軌中超聲導(dǎo)波傳播特性的計(jì)算。

1 波動(dòng)特性的振動(dòng)求解原理

在給定擾動(dòng)源及邊界條件、初始條件下，彈性體的動(dòng)力響應(yīng)有波動(dòng)解和振動(dòng)解兩種形式。介質(zhì)中擾動(dòng)質(zhì)點(diǎn)間的彈性力逐漸傳播的過程，稱為彈性波。對(duì)于有界介質(zhì)，由于擾動(dòng)在其邊界上來回反射，使得整個(gè)物體在其平衡位置附近呈現(xiàn)出一種周期性振蕩現(xiàn)象，稱為彈性體的振動(dòng)。彈性波和彈性體的振動(dòng)之間存在著本質(zhì)的內(nèi)在聯(lián)系，可以看作是同一物理問題的不同表現(xiàn)形式。即擾動(dòng)一開始總是以行波的方式將能量傳播出去，而當(dāng)物體有界時(shí)，由于行波的來回反射，最終使物體處于定常的運(yùn)動(dòng)狀態(tài)，則表現(xiàn)為振動(dòng)現(xiàn)象。彈性體的振動(dòng)是波動(dòng)過程的一種特殊形式。振動(dòng)解可表示為無窮級(jí)數(shù)形式：無窮多個(gè)駐波的迭加描述了行波的傳播，而駐波是由相同頻率的簡諧行波迭加得到。因此，通過傅里葉級(jí)數(shù)可以將波動(dòng)解與振動(dòng)解聯(lián)系起來［13］。波動(dòng)特性可以通過以下振動(dòng)分析求解得到。

其中：級(jí)數(shù)中的每一項(xiàng)代表了一個(gè)在空間具有固定模式并按一定頻率振動(dòng)的駐波。可以看出，振動(dòng)解利用

彈性體振動(dòng)方程可以表示為：

其中：M為質(zhì)量，C為阻尼系數(shù)，K為彈性系數(shù)。

利用特征頻率法對(duì)彈性體振動(dòng)方程進(jìn)行動(dòng)力學(xué)求解，得到不同特征模態(tài)的振型Φ。通過對(duì)不同頻率下不同特征模態(tài)的振型分析，可以計(jì)算出不同頻率下該模式波的波數(shù)，從而得到該模式波的相速度

其中：cp為相速度，L為模型長度，k為波數(shù)。

進(jìn)而求得其群速度

cg為群速度，Δf為頻率差，Δk為波數(shù)差。

因此，利用振動(dòng)解與波動(dòng)解互化原理，通過對(duì)結(jié)構(gòu)的振動(dòng)模態(tài)分析，可以計(jì)算出任意截面波導(dǎo)的相速度及群速度頻散特性，同時(shí)，還可以由結(jié)構(gòu)振動(dòng)分析中獲得的變形信息，計(jì)算出超聲導(dǎo)波的波結(jié)構(gòu)。下面以鋼軌異型截面波導(dǎo)為例，說明基于振動(dòng)模態(tài)分析對(duì)任意截面波導(dǎo)中超聲導(dǎo)波傳播特性計(jì)算的有效性。

2 鋼軌振動(dòng)模態(tài)的有限元模型

基于COMSOLMultiphysics有限元軟件平臺(tái)，進(jìn)行鋼軌振動(dòng)模態(tài)分析。參考GB 2585－2007標(biāo)準(zhǔn)［14］，對(duì)60 kg／m鋼軌進(jìn)行有限元分析，軌底寬為150 mm，鋼軌高度為176mm，軌頭寬為73mm。由于結(jié)構(gòu)的對(duì)稱性，在計(jì)算鋼軌中對(duì)稱模態(tài)導(dǎo)波傳播特性時(shí)，可以對(duì)其有限元分析模型進(jìn)行簡化，取鋼軌的一半進(jìn)行分析，模型如圖1所示。其中，x方向?yàn)殇撥墮M截面的寬度方向，y方向?yàn)殇撥墮M截面的高度方向，z方向?yàn)殇撥夐L度方向。

圖1 鋼軌有限元分析模態(tài)Fig.1 Rail finite element analysismode

建立鋼軌模型后，需要對(duì)其進(jìn)行參數(shù)設(shè)置，主要包括材料屬性、邊界條件、網(wǎng)格劃分、求解方式等。其中的材料性能參數(shù)設(shè)置主要包括：鋼軌密度為7 824 kg／m3，楊氏模量為2.09×1011Pa，泊松系數(shù)0.29，熱膨脹系數(shù)為1.21×10－51／K。有限元分析網(wǎng)格劃分中設(shè)定的最大單元尺寸為0.000 3，網(wǎng)格單元層數(shù)為250層。求解方式選擇特征頻率求解，設(shè)定分析的特征頻率范圍為0～50 kHz。

3 結(jié)果分析

圖2和圖3給出50 kHz低頻段范圍內(nèi)，幾個(gè)典型模態(tài)的振型云圖。從圖2可以看出，隨著頻率的增大，振動(dòng)從遍及整個(gè)截面逐漸過渡為集中在截面的局部區(qū)域，振動(dòng)趨向復(fù)雜。

從圖2（c）和圖3中I～Ⅴ幾種不同模態(tài)在46 kHz附近的振動(dòng)云圖可以看出，模態(tài)Ⅰ、Ⅱ、Ⅳ、Ⅴ振動(dòng)主要集中在軌底部分，模態(tài)Ⅲ振動(dòng)主要集中于軌頭部分，且模態(tài)振動(dòng)幅度有差別。

圖2 模態(tài)Ⅰ振型圖Fig.2 Shapes ofmodalⅠ

圖3 模態(tài)振型圖Fig.3 Modal shapes

3.1 頻散特性

在模態(tài)分析基礎(chǔ)上，根據(jù)不同頻率下不同特征模態(tài)的振型，計(jì)算出不同頻率下該模式波的波數(shù)。圖4給出通過模態(tài)分析獲得50 kHz頻率范圍內(nèi)對(duì)稱模態(tài)導(dǎo)波的波數(shù)－頻散曲線。

根據(jù)式（2）及式（3），可進(jìn)一步計(jì)算得到50 kHz頻率范圍內(nèi)鋼軌的相速度及群速度頻散曲線，結(jié)果如圖5和圖6所示。

可以看出，隨著頻率的增加，鋼軌中超聲導(dǎo)波模態(tài)的數(shù)量急劇增加，且不同模態(tài)的導(dǎo)波均存在頻散現(xiàn)象。因此，在進(jìn)行鋼軌超聲導(dǎo)波檢測(cè)參數(shù)選擇時(shí)，需要考慮超聲導(dǎo)波的這些傳播特性，通過探頭設(shè)計(jì)及檢測(cè)參數(shù)優(yōu)化避免多模態(tài)激勵(lì)，并避開頻散嚴(yán)重的檢測(cè)頻帶。

為清晰反映圖2～3中五種模態(tài)導(dǎo)波的頻散特性，圖7給出該五種模態(tài)導(dǎo)波的群速度頻散曲線。其中，模態(tài)Ⅰ為粉色圓點(diǎn)曲線，模態(tài)Ⅱ?yàn)樘焖{(lán)色圓點(diǎn)曲線，模態(tài)Ⅲ為紅色圓點(diǎn)曲線，模態(tài)Ⅳ為綠色圓點(diǎn)曲線，模態(tài)Ⅴ為藍(lán)色圓點(diǎn)曲線。

圖4 鋼軌中超聲導(dǎo)波的波數(shù)頻散曲線Fig.4 Wavenumber dispersion curves of ultrasonic guided wave in rail

圖5 鋼軌相速度頻散曲線Fig.5 Phase velocity dispersion curves of rail

圖6 鋼軌群速度頻散曲線Fig6.Group dispersion curves in rail

3.2 波結(jié)構(gòu)

根據(jù)不同頻率下各模態(tài)的位移分析，可以得到不同頻率下各模態(tài)導(dǎo)波的波結(jié)構(gòu)。圖8～圖12給出46 kHz附近Ⅰ～Ⅴ五種模態(tài)導(dǎo)波的波結(jié)構(gòu)。對(duì)比可以發(fā)現(xiàn)，模態(tài)Ⅰ、Ⅱ、Ⅳ、Ⅴ的位移主要集中在軌底部位，而模態(tài)Ⅲ的位移集中在軌頭部位。對(duì)比模態(tài)Ⅰ、Ⅱ、Ⅳ、Ⅴ沿不同方向位移分布可知，在軌底部位，模態(tài)Ⅳ沿導(dǎo)波傳播方向位移較大，且在軌底部分分布較均勻，因此，該頻率的模態(tài)Ⅳ在軌底中傳播能力強(qiáng)，對(duì)整個(gè)軌底部分的缺陷具有較高靈敏度，適合軌底無損檢測(cè)使用。

圖8 在46.4 kHz處模態(tài)Ⅰ的波結(jié)構(gòu)Fig.8Wave structure of modalⅠat 46.4 kHz

圖9 模態(tài)Ⅱ在46.6 kHz時(shí)的波結(jié)構(gòu)Fig.9 Wave structure of modalⅡat 46.6 kHz

圖10 模態(tài)Ⅲ在46.3 kHz時(shí)的波結(jié)構(gòu)Fig.10 Wave structure of modalⅢat46.3 kHz

圖11 模態(tài)Ⅳ在46.2 kHz時(shí)的波結(jié)構(gòu)Fig.11 Wave structure of modalⅣat46.2 kHz

圖12 模態(tài)Ⅴ在46.5 kHz時(shí)的波結(jié)構(gòu)Fig.12 Wave structure of modalⅤat46.5 kHz

4 結(jié) 論

利用振動(dòng)模態(tài)分析法，本文實(shí)現(xiàn)鋼軌異型截面波導(dǎo)中超聲導(dǎo)波頻散和波結(jié)構(gòu)等波動(dòng)特性計(jì)算。并在低頻50 kHz范圍內(nèi)，對(duì)五種典型模態(tài)導(dǎo)波的頻散特性和波結(jié)構(gòu)進(jìn)行了對(duì)比分析。研究發(fā)現(xiàn)，在46 kHz附近，各模態(tài)導(dǎo)波頻散相對(duì)較小，且Ⅰ、Ⅱ、Ⅳ、Ⅴ模態(tài)位移主要集中在軌底，而其中的中模態(tài)Ⅳ在軌底中傳播能力最強(qiáng)，且在軌底部位具有較均勻的位移分布，適合于軌底部位的無損檢測(cè)。模態(tài)Ⅲ的位移集中在軌頭部位，適合于軌頭部位的無損檢測(cè)。后續(xù)將基于理論分析結(jié)果，利用選定的導(dǎo)波模態(tài)類型和頻率范圍，開展軌底和軌頭部位的無損檢測(cè)實(shí)驗(yàn)研究。

［1］鄭祥明，趙玉珍，史耀武.蘭姆波頻散曲線的計(jì)算［J］.無損檢測(cè)，2003，25（2）：66－68.

ZHENG Xiang-ming，ZHAO Yu-zhen，SHI Yao-wu.Calculation for Lamb wave dispersion curves［J］.NDT，2003，25（2）：66－68.

［2］劉增華，何存富，吳 斌.單層板頻散曲線的矩陣表示及Matlab實(shí)現(xiàn)［J］.無損檢測(cè)，2005，27（5）：225－227.

LIU Zeng-hua，HE Cun-fu，WU Bin.Matrix expression for dispersion curve of single plate and its implementation by Matlab［J］.NDT，2005，27（5）：225－227.

［3］閻 石，張海鳳，蒙彥宇.Lamb波頻散曲線的數(shù)值計(jì)算及試驗(yàn)驗(yàn)證［J］.華中科技大學(xué)學(xué)報(bào)（城市科學(xué)版），2010，27（1）：1－4.

YAN Shi，ZHANG Hai-feng，MENG Yan-yu.Numerical calculation and experimental validation for Lamb wave dispersion curves［J］.J.of HUST，（Urban Science Edition）2010，27（1）：1－4.

［4］何存富，李隆濤，吳 斌.超聲導(dǎo)波在管道中傳播的數(shù)值模擬［J］.北京工業(yè)大學(xué)學(xué)報(bào)，2004，30（2）：129－133.

HE Cun-fu，LI Long-tao，WU Bin.Numerical simulation of guided ultrasonic wave transmitted in pipe［J］.Journal of Beijing University of Technology，2004，30（2）：129－133.

［5］何存富，李隆濤，吳 斌.空心圓柱體中的周向超聲導(dǎo)波［J］.機(jī)械工程學(xué)報(bào)，2004，40（8）：7－12.

HE Cun-fu，LI Long-tao，WU Bin.Guided circum ferential waves in hollow cylinders［J］.Chinese Journal ofMechanical Engineering，2004，40（8）：7－12.

［6］Lowe M.Matrix techniques formodeling ultrasonic waves in multilayeredmedia［J］.IEEE Transactions on Ultrasonics，F(xiàn)erroelectrics，and Frequency Control，1995，42（4）：525－542.

［7］Bray D E，Egle D M.Rayleigh wave dispersion in the coldworked layer of used railroad rail［J］.Journal of the Acoustical Society of America，1987，64（3）：845－851.

［8］Hirao M，Kyukawa M，Sotani Y，et al.Rayleigh wave propagation in a solid with a cold-worked surface layer［J］.Journal of Nondestructive Evaluation，1981，2（1）：43－49.

［9］Hayashi T，Song W J，Rose J L.Guided wave dispersion curves for a barwith an arbitrary cross-section，a rod and rail example［J］.Ultrasonices，2003，41（3）：175－183.

［10］lee CM，Rose JL，Cho Y.A guided wave approach to defect detection under shelling in rail［J］.NDT＆E International，2009，42（3）：174－180.

［11］Coccia S，Bartoli I，Marzani A，et al.Numerical and experimental study of guided waves for detection of defects in the rail head［J］.NDT＆E International，2011，44（1）：93－100.

［12］Hesse D，Cawley P.Surfacewavemodes in rails［J］.Journal of the Acoustical Society of America，2006，120（2）：733－740.

［13］楊桂通，張善元.彈性動(dòng)力學(xué)［M］.北京：中國鐵道出版社，1988：1－15.

［14］GB 2585－2007.鐵路用熱軋鋼軌［S］.中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局，中國國家標(biāo)準(zhǔn)化管理委員會(huì)，2008.

Propagation characteristics of ultrasonic guided wave in rails based on vibration modal analysis

HE Cun-fu，LIU Qing-qing，JIAO Jing-pin，LIU Fei，WU Bin
（College of Mechanical Engineering and Applied Electronic Technology，Beijing University of Technology，Beijing 100124，China）

Aiming at difficult problems to solve the dispersion characteristics of ultrasonic guided wave in complex waveguide structures，such as，rails，here the vibration modal analysis method was used to extract the propagation characteristics of ultrasonic guided wave.Based on the principle that vibration characteristic frequencies and wave solutions can convert each other，the solutions of characteristic frequencies obtained with vibration modal analysis were converted into wave solutions，and then the dispersion characteristics of railswere extracted.Deformation information was gained from vibration analysis of a structure，and with it the wave structure of ultrasonic guided wave was calculated.The dispersion information and wave structure information of lower ordermodes of ultrasonic guided wave were analyzed，the modal type of ultrasonic guided wave and its frequency range were selected，all of them provided a theoretical support for designing special ultrasonic guided wave sensors andmeasuring tests.

modal analysis；method of characteristic frequencies；rail；ultrasonic guided wave；dispersion characteristics；wave structure

TB553

A

國家自然科學(xué)基金項(xiàng)目（51235001，11272017）；北京市自然科學(xué)基金項(xiàng)目（1122005）

2012－12－12 修改稿收到日期：2013－03－19

何存富男，博士，教授，1958年7月生