陳友興，吳其洲，趙建輝，王召巴，石兵華，劉洋帆

(1.中北大學(xué)儀器科學(xué)與動(dòng)態(tài)測(cè)試教育部重點(diǎn)實(shí)驗(yàn)室，太原 030051;2.中北大學(xué) 電子測(cè)試技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，太原 030051;3.空軍西安飛行學(xué)院，西安 710306)

0 引言

合金火箭彈彈體具有比強(qiáng)度高、質(zhì)量輕的優(yōu)點(diǎn)，在現(xiàn)代武器中得到廣泛應(yīng)用。但其在生產(chǎn)過程中容易產(chǎn)生裂紋、氣孔、夾雜等缺陷，影響其使用性能，甚至帶來安全隱患，因此需要在其使用前進(jìn)行無(wú)損檢測(cè)。文獻(xiàn)［1－2］采用超聲檢測(cè)技術(shù)完成了小口徑火箭彈彈頭和合金彈體的缺陷檢測(cè)，并與X射線檢測(cè)結(jié)果比較，說明超聲檢測(cè)方法的有效性。在實(shí)際應(yīng)用中，由于檢測(cè)工件的特殊性，不適合用水浸方式耦合，一般采用噴水方式耦合，噴水裝置的結(jié)構(gòu)和水流會(huì)影響缺陷回波特征的提取，如果要從根本上去除噪聲，就要改進(jìn)噴水裝置的結(jié)構(gòu)。對(duì)于非通用的檢測(cè)設(shè)備，會(huì)大大增加研制成本和周期［3］。

本文在研究過程中避開了噴水裝置機(jī)械設(shè)計(jì)的問題，從信號(hào)軟處理的角度，利用自適應(yīng)濾波抵消原理抑制了噴水裝置帶來的噪聲，彌補(bǔ)了機(jī)械設(shè)計(jì)復(fù)雜的問題，方法具有通用性。此外，為了提高缺陷的準(zhǔn)確定位精度，將系統(tǒng)分析的理論引入到缺陷的超聲檢測(cè)中，利用反褶積方法準(zhǔn)確提取缺陷特征信息，實(shí)現(xiàn)了缺陷的準(zhǔn)確定位和顯示。

1 合金火箭彈彈體的超聲自動(dòng)檢測(cè)方法

圖1是一合金火箭彈彈體的剖面示意圖，由實(shí)心圓柱體、空心圓錐體和空心圓柱體3段組成。圖2為超聲自動(dòng)檢測(cè)原理示意圖。

圖1 合金火箭彈彈體剖面圖Fig.1 Cross-section drawing of the alloy rocket body

圖2 合金火箭彈彈體剖面圖Fig.2 Cross-section drawing of the alloy rocket body

圖2(a)是實(shí)現(xiàn)超聲自動(dòng)檢測(cè)的示意圖，通過旋轉(zhuǎn)工件和軸向移動(dòng)探頭實(shí)現(xiàn)整個(gè)工件的全面探測(cè)，形成C掃描圖。在工件檢測(cè)過程中，需要根據(jù)3段不同的結(jié)構(gòu)檢測(cè)方法有所區(qū)別，超聲探頭的軌跡也有所調(diào)整。當(dāng)檢測(cè)實(shí)心圓柱體段時(shí)，采用圖2(b)的方式檢測(cè)，探頭垂直于工件正下方;當(dāng)檢測(cè)空心圓柱體時(shí)，采用圖2(c)的方式檢測(cè)，探頭在圓周方向偏離工件正下方，在軸向上如圖2(a)的C所示;當(dāng)檢測(cè)圓錐體段時(shí)，探頭在圓周方向如圖2(c)所示，在軸向上如圖2(a)的B所示。本文采用噴水耦合方式，即解決了超聲檢測(cè)的耦合問題，也方便了超聲探頭的靈活調(diào)整。

本系統(tǒng)采用聚焦超聲探頭，其中心頻率為2.5 MHz，探頭晶片直徑為10 mm，探頭聚焦深度為50 mm。超聲采集卡的采樣頻率為50 MHz，8位量化采樣，量程為－1～1 V。檢測(cè)時(shí)圓周方向和軸向方向每間隔2 mm采集1次。根據(jù)文獻(xiàn)［2］的計(jì)算方法可計(jì)算出，檢測(cè)過程中不同位置的探頭姿態(tài)，即圖2(b)、(c)中的D、X和探頭旋轉(zhuǎn)角度。

圖3是采集得到的一些典型信號(hào)。由圖3可見，采用上述方法可很好地檢測(cè)出缺陷，但回波信號(hào)中還存在干擾噪聲，這些噪聲會(huì)影響到缺陷的特征提取，甚至?xí)斐扇毕莸恼`判。通過分析得出，圖3中的干擾噪聲是由于水流、水盒邊界造成，其出現(xiàn)的位置和形狀基本相似，可預(yù)先設(shè)定干擾固有參考信號(hào)，采用自適應(yīng)濾波技術(shù)進(jìn)行去噪。

圖3 檢測(cè)過程中的典型回波信號(hào)Fig.3 Typical echo signal in the ultransonic testing

2 自適應(yīng)濾波的噪聲抑制

自適應(yīng)濾波器的原理如圖4所示。其中，x(n)為輸入信號(hào);d(n)為噪聲參考信號(hào);y(n)濾波器得到近似噪聲信號(hào);r(n)是自適應(yīng)濾波去噪后的結(jié)果;e(n)是d(n)與y(n)的誤差信號(hào)。自適應(yīng)濾波器的濾波系數(shù)受到誤差信號(hào)e(n)控制，根據(jù)e(n)的值和自適應(yīng)算法自動(dòng)調(diào)整［4－5］。

根據(jù)圖4的自適應(yīng)濾波算法，以圖5的信號(hào)為參考信號(hào)，對(duì)采集的超聲原始信號(hào)進(jìn)行去噪處理。圖6列出了圖3相應(yīng)數(shù)據(jù)的結(jié)果，本文采用matlab下的adaptnlms命令實(shí)現(xiàn)自適應(yīng)濾波。

圖4 自適應(yīng)濾波原理圖Fig.4 Principle of adaptive filtering

圖5 高斯包絡(luò)噪聲參考信號(hào)模型Fig.5 Gaussian envelope noise reference signal model

圖6 自適應(yīng)濾波器處理結(jié)果Fig.6 Processing results of adaptive filtering

對(duì)比圖6和圖3可看出，經(jīng)過自適應(yīng)濾波去噪后，噪聲干擾得到很好的抑制，缺陷信號(hào)并沒有受到影響，基本保持著原有的能量和特征，說明了上述去噪方法的有效性。

3 基于反褶積技術(shù)的缺陷特征提取

要實(shí)現(xiàn)缺陷精確檢測(cè)與定位，還要根據(jù)缺陷的特征量(幅度、時(shí)間延時(shí))進(jìn)行重建得出。自適應(yīng)濾波方法雖然很好地去除了噪聲干擾，但還不能直接得出缺陷特征(幅度、時(shí)間延時(shí))。反褶積技術(shù)將超聲回波等效成激勵(lì)信號(hào)與反射系數(shù)函數(shù)的褶積，那么反射系數(shù)函數(shù)即可由超聲回波與激勵(lì)信號(hào)的反褶積得到，這樣就可由反射系數(shù)函數(shù)得出缺陷的特征量。下面介紹其原理［6］。

采用超聲反射法探測(cè)一工件時(shí)，超聲探頭發(fā)射一脈沖信號(hào)經(jīng)耦合液和工件反射后由超聲探頭接收，如果將每次反射當(dāng)作沖激函數(shù)組成，那么超聲回波信號(hào)可表示成:

式中 x(t)為超聲信號(hào);b(t)為激勵(lì)信號(hào)(子波);ξ(t)為反射系數(shù)函數(shù)。

在頻域內(nèi)的表達(dá)式可寫成:

式中 X(ω)、Bω()和 ξ(ω)分別為超聲信號(hào)頻譜、子波頻譜和反射系數(shù)函數(shù)頻譜。

如果令:

則

在時(shí)域內(nèi):

由此可知，可根據(jù)超聲激勵(lì)信號(hào)，求出a(t)，再利用式(5)反射系數(shù)函數(shù)ξ(t)。

圖7是圖6(b)、(d)經(jīng)反褶積得到的結(jié)果。從圖7可看出，缺陷特征脈沖明顯，可很容易得出代表缺陷的信息(幅度、時(shí)間延時(shí))。

圖7 基于反褶積的缺陷特征提取結(jié)果Fig.7 Feature extraction results of defect based on deconvolution method

4 缺陷的C掃描顯示與定位

4.1 缺陷檢測(cè)與C掃描顯示

對(duì)合金火箭彈按第2章描述的方法進(jìn)行檢測(cè)，將檢測(cè)的數(shù)據(jù)采用自適應(yīng)濾波噪聲抑制方法和反褶積技術(shù)進(jìn)行處理，得出有缺陷的幅值和時(shí)間延時(shí)特征量，圖8是根據(jù)缺陷幅值特征量進(jìn)行C掃描重建得到的檢測(cè)結(jié)果。從圖8可看出，工件在實(shí)心圓柱體、空心圓錐體和空心圓柱體部分都存在缺陷。

為了證明檢測(cè)方法的有效性，采用X射線對(duì)工件進(jìn)行檢測(cè)。根據(jù)工件的結(jié)構(gòu)，采用2次不同劑量透射得到2張照片，圖9是2張照片拼接的結(jié)果。從圖9可看出，工件中共存在4處缺陷:實(shí)心圓柱體部分2個(gè)孔狀缺陷(1.5 mm人工孔)、空心圓錐體部分1個(gè)內(nèi)表面劃痕缺陷(人工)和空心圓柱體部分1個(gè)孔狀缺陷(1.5 mm 人工孔)。

圖8 C掃描重建缺陷結(jié)果Fig.8 C-scan image of the defect

圖9 X射線檢測(cè)結(jié)果Fig.9 Result of X-ray testing

對(duì)比圖8和圖9可看出，超聲C掃描能夠檢測(cè)出缺陷的有無(wú)和在軸向上的位置，但不能對(duì)缺陷在圓周方向進(jìn)行定位，也沒辦法區(qū)分在同一圓周截面上的2處缺陷。因此，僅采用C掃描圖不能很好地體現(xiàn)缺陷的完整信息。下面根據(jù)缺陷的時(shí)間延時(shí)特征量進(jìn)行缺陷定位。

4.2 缺陷定位的計(jì)算方法

根據(jù)對(duì)缺陷特征信號(hào)的時(shí)間延時(shí)，可推導(dǎo)出缺陷在工件圓周截面的位置(ρ，θ)。其中，ρ為缺陷距離工件中心的距離;θ為缺陷偏移水平軸的角度。由于實(shí)心體部分和空心體部分的檢測(cè)方式有所不同，缺陷的定位方法也有所不同。

(1)實(shí)心體部分

根據(jù)圖1(b)，探頭垂直入射工件下表面，故

式中 R為彈體橫截面半徑;l為超聲波入射點(diǎn)位置與缺陷處的直線距離，l=CΔt/2;C為合金材料聲速;Δt為下表面回波特征與缺陷特征的時(shí)間間隔(如圖7(a)所示)。

由于彈體轉(zhuǎn)動(dòng)的是由電機(jī)控制，所以彈體1圈檢測(cè)的總點(diǎn)數(shù)N是已知的，用n表示對(duì)彈體進(jìn)行1圈檢測(cè)中的第n個(gè)檢測(cè)數(shù)據(jù)，則

(2)空心體部分

根據(jù)圖1(c)所示，超聲波偏離中心軸X入射工件表面，此時(shí)屬于斜入射情況，那么根據(jù)超聲波的傳播規(guī)律遵循Snell定律和幾何關(guān)系，可得

式中 α為入射角;β為折射角;R為工件檢測(cè)位置的外徑;C為合金材料聲速;C0為水聲速。

根據(jù)三角形余弦定理，可得到缺陷所在位置到圓心的距離:

與式(6)一樣，這里l=CΔt/2。

該缺陷位置所對(duì)應(yīng)極坐標(biāo)圓心角為

實(shí)心體部分根據(jù)式(6)和式(7)，空心體部分根據(jù)式(9)和式(10)可計(jì)算出缺陷的位置，圖10是3個(gè)不同部位的缺陷定位圖。

圖10 缺陷定位圖Fig.10 Defect location map

比較圖10的定位結(jié)果、圖9的X射線檢測(cè)結(jié)果以及工件的實(shí)際的缺陷情況，發(fā)現(xiàn)圖10的缺陷定位很準(zhǔn)確，基本能夠確定缺陷所在的位置。說明采用的檢測(cè)方法和信號(hào)處理技術(shù)是有效的。

5 結(jié)論

(1)通過與X射線檢測(cè)結(jié)果比較得出，采用的超聲自動(dòng)檢測(cè)方法在提高檢測(cè)效率時(shí)，可有效地檢測(cè)出合金火箭彈彈體內(nèi)部的缺陷。

(2)針對(duì)噴水耦合帶來的噪聲干擾，采用自適應(yīng)濾波技術(shù)進(jìn)行噪聲抑制，避免了噴水裝置機(jī)械設(shè)計(jì)復(fù)雜的問題。結(jié)果表明，在噪聲得到很好抑制的同時(shí)，確保缺陷信息能量和特征不受影響。

(3)將系統(tǒng)分析的理論引入到缺陷的超聲檢測(cè)中，采用反褶積技術(shù)進(jìn)行缺陷特征量的提取，通過缺陷定位算法驗(yàn)證特征量提取的準(zhǔn)確性。

［1］金永，王召巴，丁戰(zhàn)陽(yáng)，等.一種小口徑火箭彈彈頭超聲檢測(cè)方法［J］.固體火箭技術(shù)，2010，33(1):115-118.

［2］戚勵(lì)文，金永，王召巴，等.合金彈體棒狀坯料超聲檢測(cè)方法研究［J］.彈箭與制導(dǎo)學(xué)報(bào)，2009，29(4):269-272.

［3］原可義，韓贊東，王柄方等.復(fù)合材料噴水耦合超聲C掃查檢測(cè)系統(tǒng)的研制［J］.航空制造技術(shù)，2009，(15):101-103.

［4］Paulo S R Diniz .Adaptive filtering［M］.Springer，2013:209-213.

［5］高鷹，謝勝利.一種變步長(zhǎng)LMS自適應(yīng)濾波算法及分析［J］.電子學(xué)報(bào)，2001，29(8):1-3.

［6］Wei Liang，Huang Zuo-ying，Que Pei-wenb.Sparse deconvolution method for improving the time-resolution of ultrasonic NDE signals［J］.NDT＆E International，2009(42):430-434.