王 麗,陳振茂

(1.西安交通大學(xué) 強(qiáng)度與振動教育部重點實驗室,西安 710049;2.西安郵電學(xué)院 理學(xué)院,西安 710121)

為了確保系統(tǒng)安全運行、防止意外事故發(fā)生,對核電站等大型機(jī)械系統(tǒng)的定期檢查是必不可少的。為了根據(jù)其形狀大小決定相應(yīng)的維護(hù)對策,不僅需要檢測裂紋是否存在,還需要對裂紋進(jìn)行無損定量。合理的維護(hù)對策不但可以減少維修費用,而且可以減少因停機(jī)帶來的經(jīng)濟(jì)損失。因此裂紋無損定量對確保結(jié)構(gòu)安全及提高維護(hù)效率至關(guān)重要。

應(yīng)力腐蝕裂紋(SCC)是金屬材料在某些特定介質(zhì)和拉伸應(yīng)力共同作用下所發(fā)生的延遲破裂現(xiàn)象。隨著石油、化工、冶金和核電等工業(yè)部門的發(fā)展,SCC 的危害性已成為突出問題。因此,對SCC 的無損定量研究具有重要的意義。另一方面,SCC 主要沿晶界進(jìn)展,具有類似于樹枝分叉的復(fù)雜微觀結(jié)構(gòu),其無損定量非常困難,成了產(chǎn)業(yè)界迫切需要解決的問題。

渦流檢測(ECT)由于其方便、快速和對導(dǎo)體材料表面和近表面缺陷的高檢測靈敏度等優(yōu)勢,作為超聲檢測的一個有效補(bǔ)充,在解決實際復(fù)雜裂紋定量問題上被寄予厚望,目前已經(jīng)成功地應(yīng)用于包括SCC 在內(nèi)的一些自然裂紋的重構(gòu)[1-6]。然而,對于深度較大的SCC,即使使用導(dǎo)電性裂紋模型,在重構(gòu)信號和檢測信號相符很好時,裂紋深度的重構(gòu)精度也可能很不理想[7]。這可能是因為SCC 的形狀、電導(dǎo)率等在裂紋區(qū)域發(fā)生變化的結(jié)果。因此,對復(fù)雜形狀SCC 的重構(gòu),需要提出更符合實際的數(shù)值計算模型。

筆者通過數(shù)值模擬試驗和對人工SCC 電導(dǎo)率的檢測,在分析SCC 形態(tài)特征和電導(dǎo)率分布規(guī)律的基礎(chǔ)上,建立了分區(qū)域賦值的SCC 精細(xì)數(shù)值計算模型,改善了SCC 重構(gòu)中深度的欠評價問題。數(shù)值模擬和實際渦流檢測信號的重構(gòu)結(jié)果驗證了新模型對于SCC 重構(gòu)的有效性。

1 SCC 精細(xì)數(shù)值計算模型

渦流無損定量中SCC 等實際裂紋的理論模型尚未建立,其檢測信號的計算仍然是一個難題。電火花加工(EDM)裂紋模型假設(shè)裂紋區(qū)域電導(dǎo)率為零,但SCC 電導(dǎo)率不為零,有部分電流可流過裂紋面,因此EDM 裂紋模型不適合于SCC 的重構(gòu)。即使使用導(dǎo)電性裂紋模型,在裂紋區(qū)如果只使用統(tǒng)一的電導(dǎo)率,重構(gòu)結(jié)果誤差仍然較大。文獻(xiàn)[8]利用二維EC T 信號的分段逆策略重構(gòu)了一類復(fù)雜形狀SCC 的深度及電導(dǎo)率分布,重構(gòu)精度有所改善。在此基礎(chǔ)上,結(jié)合SCC 的切片電導(dǎo)率測量結(jié)果,筆者提出了對裂紋的電導(dǎo)率分區(qū)域賦值的SCC 精細(xì)數(shù)值計算模型。

該方法為:將被檢對象的可能裂紋區(qū)域分成層厚固定的若干層1,2, …,n,按照裂紋信號的特征,將每一層分為若干個區(qū)域1,2, …,m,每個小區(qū)域的電導(dǎo)率σij(i=1,2, …,m;j=1,2, …,n)設(shè)為一個待求值。為了減少參數(shù)個數(shù)、降低逆問題的不適定性,在深度方向設(shè)裂紋區(qū)域的電導(dǎo)率從上到下按線性規(guī)律變化,即σij=σi1+ki·σi1(i=1,2, …,m;j=1,2,…,n)。因此,需重構(gòu)的裂紋參數(shù)包括裂紋的形狀參數(shù)和電導(dǎo)率參數(shù)σi1和ki(i=1,2, …,m)。

2 基于粒子群算法和共軛梯度法的混合逆策略

由于SCC 的建模需要考慮裂紋區(qū)域電導(dǎo)率分布,所以重構(gòu)變量包括裂紋區(qū)域的電導(dǎo)率和裂紋的形狀參數(shù)。數(shù)值計算結(jié)果表明,裂紋形狀參數(shù)和電導(dǎo)率對裂紋信號的影響有很大不同,如果采用相同逆策略進(jìn)行重構(gòu),則會由于收斂速度的不同而導(dǎo)致重構(gòu)結(jié)果誤差較大。為此,筆者利用基于共軛梯度法和粒子群算法的混合逆策略,對裂紋形狀參數(shù)和電導(dǎo)率分別進(jìn)行重構(gòu)。重構(gòu)的基本思路是:首先基于裂紋初始電導(dǎo)率分布,通過共軛梯度法[9],修正裂紋的形狀參數(shù),直至殘差趨于穩(wěn)定;然后根據(jù)殘差的大小,利用粒子群算法,修正裂紋的電導(dǎo)率分布;以此反復(fù)迭代,直至達(dá)到收斂要求。

圖1 SCC 數(shù)值模擬信號對比

3 數(shù)值結(jié)果

3.1 數(shù)值模擬信號的重構(gòu)

首先,以SUS304 不銹鋼板(200 mm×100 mm×10 mm, 電導(dǎo)率1.4 MS/m)中長14 mm , 寬0.2 mm,深6 mm 的矩形裂紋信號作為檢驗信號來驗證新模型的有效性。渦流檢測線圈采用十字探頭(內(nèi)徑4.5 mm ,外徑5.5 mm,厚度0.5 mm)。采用10 kHz 激勵頻率和0.5 mm 的提離距離,當(dāng)使用新模型重構(gòu)時,信號基本符合(圖1),深度重構(gòu)結(jié)果合適(6 mm),電導(dǎo)率的重構(gòu)結(jié)果見表1。而使用均勻電導(dǎo)率模型時, 當(dāng)裂紋電導(dǎo)率分別取為4.45%,4.76%和4.94%時,重構(gòu)的深度分別是5,6,7 mm,即隨著裂紋電導(dǎo)率的增加,重構(gòu)深度會增加,重構(gòu)的信號呈整體上升趨勢,相應(yīng)信號與檢測信號相差較大。鑒于SCC 的復(fù)雜微觀結(jié)構(gòu),應(yīng)力腐蝕裂紋分區(qū)域賦值的精細(xì)數(shù)值計算模型更接近于實際,因此有望提高SCC 無損定量的精度。

表1 數(shù)值模擬電導(dǎo)率分布重構(gòu)結(jié)果與母材電導(dǎo)率之比∶%

3.2 自然缺陷試件檢測信號的重構(gòu)

圖2 和表2 所示的是采用10 kHz 激勵頻率和0.5 mm 的提離距離, 針對切斷觀察結(jié)果為長41 mm、平均深度4.7 mm 、最大深度5.1 mm 的SCC 進(jìn)行重構(gòu)的結(jié)果。使用新模型重構(gòu)的平均深度為5.0mm ,長度是34mm ,接近實際裂紋,重構(gòu)信號與實驗信號也符合較好。而常規(guī)方法使用2%的電導(dǎo)率重構(gòu)時最深處結(jié)果為4.3 mm,使用EDM裂紋模型重構(gòu)時最深處結(jié)果為3.6 mm,重構(gòu)信號與實驗信號誤差相對較大。

圖2 自然缺陷試件檢測信號重構(gòu)結(jié)果

表2 自然缺陷試件的電導(dǎo)率分布重構(gòu)結(jié)果與母材電導(dǎo)率之比∶%

以上根據(jù)數(shù)值實驗和SCC 電導(dǎo)率的測定結(jié)果,采用對裂紋區(qū)域電導(dǎo)率分區(qū)域賦值的思想,建立了SCC 精細(xì)數(shù)值計算模型,有效改善了SCC 深度反演中的欠評估問題?；跀?shù)值模擬信號和實際渦流檢測信號的重構(gòu)結(jié)果,驗證了所提模型和方法對SCC重構(gòu)的有效性。

[1] Yusa N, Chen Z, Miya K.Quantitative profile evaluation of natural crack in steam generator tube from eddy current signals[J] .Int J Appl Electromagn Mech,2000, 12(3-4):139-150.

[2] Cheng W, Kanemoto S, Komura I.Depth sizing of partial contact SCC from ECT signals[J] .NDT&E Int, 2006,39(5):374-383.

[3] Yusa N, Chen Z, Miya K.Large-scale parallel computation for the reconstruction of natural stress corrosion cracks from eddy current testing signals[J] .NDT&E Int, 2003,36(7):449-459.

[4] Yusa N, Cheng W Y, Chen Z M .Generalized neural netw ork approach to eddy current inversion to the analysis of real cracks[J] .NDT&E Int, 2002, 35(8):609-614.

[5] Yusa N, Miya K.Discussion on the equivalent conductivity and resistance of stress corrosion cracks in eddy current simulations[J] .NDT&E Int, 2009,42(1):9-15.

[6] Yusa N, Hashizume H .Evaluation of stress corrosion cracking as a function of its resistance to eddy currents[J] .Nuclear Engineering and Design, 2009, 239(12):2713-2718.

[7] Yusa N, H uang H, Miya K.Numerical evaluation of the ill-posedness of eddy current problems to size cracks[J] .NDT&E Int,2007(40):185-191.

[8] Wang L, Chen Z, Lu T J.Sizing of long stress corrosion crack from 2D EC T signals by using a multiseg ment inverse analysis strategy[J] .Int J Appl Electromagn Mech,2008, 28(1-2):155-161.

[9] Norton S J, Bow ler J R.Theory of eddy current inversion[J] .J Appl Phys,1993,73(2):501-512.