不銹鋼增材制造件表面缺陷激光超聲檢測

房海基　葉國良　呂波　張艷喜　高向東

摘要：利用激光超聲技術(shù)檢測增材制造件表面缺陷，采用有限元方法模擬熱彈機制下激光激發(fā)超聲波的傳播過程，對不同探測位置接收的表面波反射信號進行了分析，研究缺陷的深度和寬度對表面波反射信號的影響。對316L不銹鋼增材制造件進行了激光超聲檢測試驗，驗證了該模型的正確性。利用小波軟閾值去噪算法處理采集到的激光超聲信號。結(jié)果表明，數(shù)值模擬與試驗結(jié)果基本一致，表面波與缺陷作用產(chǎn)生的RS波和RR波的到達時間差可以檢測缺陷的深度，缺陷的寬度對檢測結(jié)果幾乎無影響。

關(guān)鍵詞：激光超聲檢測;數(shù)值模擬;表面缺陷;增材制造 中圖分類號：TG454

前言

金屬增材制造是通過材料逐點或逐層疊加的方法制造實體零件的技術(shù)。相較于傳統(tǒng)去除成形的減材制造方法，金屬增材制造有特殊的成形工藝，可快速而精確地制造出任意復雜形狀的零件＼[1–3＼]。但該工藝的特殊性也使增材制造件表現(xiàn)出孔隙性、開裂和分層等致命的弱點，會降低制件的合金性能＼[4-5＼]。特別是表面裂紋的孕育擴展而造成瞬間斷裂，嚴重影響著生產(chǎn)安全。因此，需要對金屬結(jié)構(gòu)件進行早期診斷、評估，及時修復裂紋或更換結(jié)構(gòu)件。

金屬增材制造件的結(jié)構(gòu)特點要求無損檢測手段能夠減小檢測盲區(qū)，并且實現(xiàn)精密測量，而目前國內(nèi)外的研究主要集中在超聲、電磁、和計算機斷層掃描等檢測方法＼[6＼]。激光超聲檢測技術(shù)利用脈沖激光在材料表面激發(fā)出縱波、橫波和聲表面波等多模式超聲波來檢測缺陷，具有非接觸、靈敏度高和空間分辨率高等優(yōu)點＼[7-8＼]，被廣泛用于表面、亞表面和內(nèi)部缺陷的檢測。針對金屬材料的表面裂紋缺陷，國內(nèi)外很多學者進行了研究。Zhou等人＼[9＼]、Guan等人＼[10＼]從試驗和數(shù)值模擬角度，解釋了激光激發(fā)表面波與材料表面開口缺陷的作用機理?；诓煌疃取挾鹊娜毕菰跁r間、頻率、能量反射率和透射率上的不同，實現(xiàn)缺陷的定位和定量。王威等人＼[11＼]發(fā)現(xiàn)激光遠場激發(fā)時有利于各種模式轉(zhuǎn)換波的分離，散射回波特征點到達時間差與缺陷深度、寬度有線性關(guān)系。Kromine等人＼[12＼]提出了移動掃描光源技術(shù)，當脈沖激光掃描至缺陷的近場區(qū)域時，信號幅度不斷增大，從而確定裂紋的位置。Ni等人＼[13＼]提出了一種雙光源檢測系統(tǒng)，利用散射超聲模式傳播時間計算裂紋取向角，對數(shù)值模擬和試驗結(jié)果進行了比較，證明了該方法的準確性。

關(guān)于增材制造件缺陷檢測已有許多研究，包括在線檢測和無損檢測。由增材制造過程和零件使用過程引起的表面微缺陷對傳統(tǒng)無損檢測方法提出了挑戰(zhàn)。目前，利用激光超聲檢測技術(shù)檢測增材制造件質(zhì)量被證實是可行的，且具有一定的優(yōu)勢，得到研究者的廣泛關(guān)注。以316L不銹鋼增材制造件為研究對象，建立激光激發(fā)超聲波的數(shù)值模型，分析表面波與缺陷相互作用的傳播過程，研究不同探測位置、不同深度和寬度缺陷對表面波反射波信號的影響。最后結(jié)合激光超聲試驗進行驗證，根據(jù)表面波與缺陷作用產(chǎn)生的RS波與RR波的到達時間，實現(xiàn)對缺陷的位置、深度的檢測。

1理論與數(shù)值模型

1.1激光超聲理論

根據(jù)脈沖激光的光功率密度強度和材料特性的不同，激光超聲的激發(fā)機制主要包括熱彈機制和燒蝕機制。為了滿足對金屬材料無損檢測的要求，采用激光熱彈機制進行仿真與試驗。在熱彈機制下，當較低光功率密度強度的脈沖激光入射到材料表面時，部分激光能量被吸收并轉(zhuǎn)換為熱能，在輻照區(qū)域附近迅速地升溫，導致快速地熱膨脹，從而產(chǎn)生超聲波。

激光輻照材料激發(fā)表面波的控制方程＼[14＼]為：

式中：T（x，y，t）為t時刻的溫度分布;U（x，y，t）為t時刻的瞬態(tài)位移;T0為環(huán)境溫度;β=（3λ+2μ）α為材料熱彈耦合系數(shù)，λ和μ為Lame常數(shù);α為材料的熱膨脹系數(shù);ρ，c和k分別為材料的密度、比熱容和熱傳導系數(shù)。忽略由應力場產(chǎn)生的微小應變對材料表面的影響，T0βSymbolQC@

·U·項可去除，式（1）進一步化簡為：

作用于材料上表面的激光能量為：

式中：I0為激光中心功率密度;A（T）為材料表面的吸收率;f（x），g（t）為激光的空間分布和時間分布，分別表示為：

式中：r0，t0分別為激光光斑半徑和激光脈沖上升時間。

1.2數(shù)值模型建立

圖1中的數(shù)值模型材料是均勻、各向同性的316L不銹鋼，其材料參數(shù)見表1。該簡化二維模型的整體尺寸為35mm×10mm，表面缺陷的寬度為w，深度為h，缺陷距離模型左邊界23mm。激光作用于距離模型左邊界10mm處，在缺陷左側(cè)3mm處設(shè)置信號接收點Det1，用于采集表面波反射信號。數(shù)值模擬過程中將材料的上表面設(shè)置為自由邊界，模型的側(cè)邊和底面設(shè)置為吸收層，具有低反射特點，能夠減少其它邊界反射波帶來的影響。

激光的脈沖上升時間為10ns，激光光源的能量為1mJ，激光光斑半徑為0.3mm。由于激光作用在金屬材料上時，只有表面和亞表面能夠吸收激光的能量產(chǎn)生熱膨脹。為了精確模擬出激光超聲產(chǎn)生彈性波在材料中的傳播，這里選用自由四邊形網(wǎng)格，在激光輻照附近區(qū)域的網(wǎng)格大小為20μm，其余區(qū)域網(wǎng)格大小為100μm。另外求解器設(shè)置最重要的一個是時間步長的設(shè)置，小的時間步長將提高精度，但過于小的時間步長又將減緩模擬效率。這里取時間步長為10ns，計算總時長10μs。

2數(shù)值模擬結(jié)果

2.1瞬時聲場分布

圖2所示為激光作用在316L不銹鋼材料時，激發(fā)超聲波在不同時刻的聲場分布圖。從圖2a可以看出材料表面吸收激光的瞬時高能量而激發(fā)出表面波（Rwave）、橫波（Swave）和縱波（Lwave）等，表面波只集中在表面附近傳播，并且在深度上存在振動幅度的降低。圖2b和圖2c分別是在4.66μs和5.65μs的聲場分布，表面波在遇到缺陷時會發(fā)生模式轉(zhuǎn)換，產(chǎn)生反射波（Rr）和透射波（Rt），其中大部分能量被反射。

試驗發(fā)現(xiàn)表面反射波中含有大量表面缺陷相關(guān)信息，可以作為檢測表面缺陷的有效數(shù)據(jù)，進而通過信號特征有效地提取缺陷特征。由于表面波的質(zhì)點振動位移有2個分量：1個垂直于表面，1個沿著表面內(nèi)波的傳播方向，所以可以在模型上表面測量位移確定反射波的到達時間和振動幅度。

2.2不同探測位置接收的信號

在數(shù)值模型表面距離激光中心8mm，9mm，10mm，11mm和12mm點處探測表面波信號，保持表面缺陷的寬度為0.2mm，深度為1mm，與激光中心的距離13mm不變。獲取表面波信號，如圖3所示。可以看到不同探測點接收的表面波信號時域波形基本相似，依次出現(xiàn)掠面縱波、掠面橫波、表面波和各類反射波，有兩個明顯的波峰，分別是直達表面波（R）和反射表面波（RR）。直達表面波就是從激光作用中心傳播到信號探測點處引起的，反射表面波則是表面波與缺陷相互作用后形成的反射波。進一步觀察該圖，隨著探測點遠離激光作用中心而靠近缺陷時，直達表面波和反射表面波的到達時間都呈線性變化，容易混疊在一起，但其幅值在該范圍內(nèi)受到的影響較小。

2.3缺陷深度、寬度對檢測的影響

上述分析只針對一種缺陷深度，為了進一步分析表面波與不同尺寸表面缺陷的相互作用規(guī)律，在圖1模型的基礎(chǔ)上，在Det1處設(shè)置信號接收點，測量垂直表面位移。如圖4所示是測量深度分別為0.5mm，1.0mm，1.5mm和2.0mm，寬度均為0.2mm的矩形表面缺陷的位移波形圖。圖4a中，可以看到表面波（R）的能量較強，Rr放大圖中有2個明顯的信號RR和RS，RS波為缺陷底端向樣品表面?zhèn)鞑サ哪Ｊ睫D(zhuǎn)換橫波。RR波的振動幅度隨著深度的增加而有顯著的增大，RS波的振動幅度變化則相反，說明隨著深度的增加，更多的能量沿原路徑直接反射。R波和RR波的到達時間基本相同，而RS波的到達時間則隨著深度的增加而延后，并呈線性關(guān)系，說明RR波的到達時間僅僅與缺陷的位置有關(guān)，而缺陷的深度測量則依賴于RS波的到達時間。

在圖1所示的模型中，保持缺陷的深度為1.0mm不變，寬度分別為0.2mm，0.5mm和1.0mm時，在相同探測點處接收到的波形發(fā)生改變，如圖5所示。由前面的分析可知，反射表面波RR的到達時間與缺陷的尺寸無關(guān)，從圖5來看，不同寬度缺陷的RR波到達時間基本一致，而RS波的到達時間隨著寬度的增加有所延遲，但變化不明顯，說明缺陷的寬度變化對反射波的傳播并無主導性的影響，RS波的到達時間只與缺陷深度有關(guān)。但隨著寬度的增大，RS波的幅值大小有輕微的下降。

3激光超聲試驗與結(jié)果對比

3.1激光超聲試驗平臺與樣品

如圖6所示為激光超聲檢測試驗平臺。主要包括激光器（Nimma400型）、激光測振儀（TWM532型）、示波器（TekMDO3024型）、光電探測器、夾具和光學配件（三棱鏡、聚焦鏡）等。試驗中主要對316L不銹鋼材料的激光選區(qū)熔化（SLM）增材制造樣品進行檢測，如圖7所示的2塊100mm×15mm×10mm的樣品表面切割出深度為0.5mm，1.0mm，1.5mm和2.0mm，寬度為0.5mm的槽，表面經(jīng)過加工處理。316L不銹鋼粉末粒度為15～53μm，打印參數(shù)為：激光功率250W，掃描速度1000mm/s，打印線寬100μm，打印層厚30μm，打印精度0.05mm。試驗中激光激發(fā)點和信號探測點的位置與仿真保持一致。每個信號采集點的采樣率設(shè)置為500M/s，每組數(shù)據(jù)的采樣長度為10k。

3.2有限元模型驗證

為了驗證激光激發(fā)表面波有限元模型的正確性，首先對無缺陷模型的仿真信號和試驗信號進行了對比，在圖1的激光激發(fā)點位置右側(cè)9mm的A點處接收表面波信號，結(jié)果如圖8所示。2個信號中特征較明顯的是直接從激發(fā)源產(chǎn)生首達表面波，其到達時間分別為3.3μs和3.203μs，計算得到表面波在316L不銹鋼材料中的速度分別為2730m/s和2810m/s，兩者誤差大約為2.84%。進一步，對有缺陷模型的仿真信號和試驗信號進行了對比，信號接收點距離激光源10mm（B點），結(jié)果如圖9所示。表面波遇到缺陷反射回到信號接收點的總路程為16mm，測得2個信號的反射表面波到達時間分別為5.89μs和5.68μs，計算得到傳播速度分別為2716m/s和2817m/s，誤差為3.58%。對比無缺陷模型和有缺陷模型的表面波波形和計算得到的傳播速度，可以看到仿真模型和試驗結(jié)果基本一致，說明仿真模型的正確性。

3.3信號處理與結(jié)果對比

在實際表面缺陷檢測過程中，由于外界因素，激光超聲信號中除了表現(xiàn)缺陷特征的有用信號，還摻雜著噪聲。小波閾值去噪法作為應用最廣泛的去噪方法，其基本原理是對經(jīng)過小波分解的小波系數(shù)與選定的閾值進行比較，去除噪聲部分的系數(shù)，并將有用信號部分的系數(shù)閾值化處理，最后重構(gòu)成降噪后的信號。這里選用小波軟閾值函數(shù)，其表達式為：

式中：閾值T=σ2lnN;其中N為采樣點數(shù);σ為第k層噪聲的標準差，按下式估算獲得：

式中：median（x）表示取第k層分解小波系數(shù)絕對值的中值。

圖10為試驗中檢測缺陷的表面波形。圖11為檢測0.5mm和1.0mm深度缺陷時經(jīng)過小波降噪的信號，這里采用db4作為小波基函數(shù)，進行5層小波分解?？梢钥闯觯边_表面波的幅值較高，經(jīng)過表面缺陷反射作用后，接收到較為微弱的回波及各類模態(tài)轉(zhuǎn)換后的復雜波形，經(jīng)過小波去噪的表面波信號不但去除了噪聲，并且在最大程度上保留了很多和缺陷相關(guān)的特征。經(jīng)過分析，在檢測0.5mm深缺陷時，RR波和RS波的到達時間為5.68μs和5.94μs，時間差為0.26μs，檢測1.0mm深缺陷時，RR波和RS波的到達時間為5.56μs和6.02μs，時間差為0.46μs。為了避免單次試驗的偶然性，對每個深度的缺陷進行3次數(shù)據(jù)采集并計算RS波與RR波的到達時間差，其它試驗條件保持一致。增材制造缺陷檢測試驗結(jié)果與仿真結(jié)果對比如圖12所示。從對比可以看出，試驗與仿真結(jié)果基本一致，造成一定偏差的原因主要是：①實際線切割樣品缺陷帶來的尺寸偏差;②每次試驗移動樣品位置時存在人為誤差;③表面處理導致樣品表面粗糙度與吸收的激光能量差異。另外，增材制造件的成形工藝表明在鋪粉方向和樣品沉積方向的表面波傳播存在一定的差異，這也是試驗結(jié)果偏差的原因。但從擬合曲線看，線性趨勢基本一致，利用RS波與RR波的到達時間差能夠有效地檢測增材制造件表面缺陷的深度。

4結(jié)論

以316L不銹鋼增材制造件為研究對象，基于激光超聲的熱彈機制，通過有限元方法建立含有表面缺陷的數(shù)值模型，利用激光激發(fā)表面波與表面缺陷作用的反射回波信號，研究表面波在不同深度、寬度缺陷的傳播特征。通過激光超聲試驗進行波速計算驗證了該模型的正確性，深度檢測試驗結(jié)果表明：數(shù)值模擬與試驗結(jié)果基本一致。激光激發(fā)表面波與缺陷作用產(chǎn)生的RS波和RR波的到達時間差可以檢測缺陷深度，并且到達時間差隨缺陷深度增加呈線性變化。缺陷的寬度變化對缺陷位置、深度的檢測幾乎無影響。同時表明激光超聲檢測增材制造件表面缺陷的有效性，為增材制造在線檢測提供參考。

參考文獻

[1]楊義成，黃瑞生，方乃文，等.光粉交互對同軸送粉增材制造能量傳輸?shù)挠绊懀躘J＼].焊接學報，2020，41（6）：19-23.

＼[2＼]王中豪，韓紅彪，陳俊潮，等.一種逐條預置焊絲和逐層重熔的增材制造方法＼[J＼].焊接，2020（3）：35-39.

＼[3＼]趙慧慧，尹晨豪，封小松，等.旁路等離子在航天用軸結(jié)構(gòu)增材修復中的應用＼[J＼].焊接，2020（1）：47-50.

＼[4＼]ZhaoC，ParabND，LiXX，etal.Criticalinstabilityatmovingkeyholetipgeneratesporosityinlasermelting＼[J＼].Science，2020，370：1080-1086.

＼[5＼]YinYan，ZhangYuan，DongKaiji，etal.Thedevelopmentof3DprintingtechnologyandthecurrentsituationofcontrollingdefectsinSLMtechnology＼[J＼].ChinaWelding，2020，29（3）：9-19.

＼[6＼]孫長進，趙宇輝，王志國，等.增材新概念結(jié)構(gòu)無損檢測技術(shù)發(fā)展現(xiàn)狀及趨勢研究＼[J＼].真空，2019，56（4）：65-70.

＼[7＼]DixonS，BurrowsSE，DuttonB，etal.DetectionofcracksinmetalsheetsusingpulsedlasergeneratedultrasoundandEMATdetection＼[J＼].Ultrasonics，2011，51（1）：7-16.

＼[8＼]CernigliaD，MontinaroN.Defectdetectioninadditivelymanufacturedcomponents：laserultrasoundandlaserthermographycomparison＼[J＼].ProcediaStructuralIntegrity，2018，8：154-162.

＼[9＼]ZhouZG，ZhangKS，ZhouJH，etal.Applicationoflaserultrasonictechniquefornoncontactdetectionofstructuralsurfacebreakingcracks＼[J＼].Optics&LaserTechnology，2015，73：173-178.

＼[10＼]GuanJF，ShenZH，NiXW，etal.Numericalsimulationofthereflectedacousticwavecomponentsinthenearfieldofsurfacedefects＼[J＼].JournalofPhysicsDAppliedPhysics，2006，39（6）：1237-1243.

＼[11＼]王威，仲政，潘永東.激光遠場激發(fā)表面波在開口缺陷處的散射回波＼[J＼].激光技術(shù)，2015，39（2）：157-165.

＼[12＼]KromineAK，F(xiàn)omitchovPA，KrishnaswamyS，etal.Applicationsofscanninglasersourcetechniquefordetectionofsurfacebreakingdefects＼[J＼].SymposiumonAppliedPhotonics，2000，4076：252-259.

＼[13＼]NiCY，DongLM，ShenZH，etal.Theexperimentalstudyoffatiguecrackdetectionusingscanninglaserpointsourcetechnique＼[J＼].Optics&LaserTechnology，2011，43（8）：1391-1397.

＼[14＼]AchenbachJD.Laserexcitationofsurfacewavemotion＼[J＼].JournaloftheMechanicsandPhysicsofSolids，2003，51（11-12）：1885-1902.

房海基簡介：碩士;主要從事焊接自動化方向的研究。

基金項目：廣州市科技計劃（202002020068，202002030147）

源文獻：房?；?，葉國良，呂波，等.不銹鋼增材制造件表面缺陷激光超聲檢測＼[J＼].焊接，2021（9）：16-21.