徐 帥， 程 軍， 楊繼全， 劉益劍， 裘進(jìn)浩

(1. 南京師范大學(xué)江蘇省三維打印裝備與制造重點(diǎn)實(shí)驗(yàn)室 南京，210023) (2. 南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室 南京，210016)

引 言

復(fù)合材料是由兩種或兩種以上單一材料復(fù)合而成的多相材料[1]。與單一材料相比，復(fù)合材料因其輕質(zhì)高強(qiáng)等眾多優(yōu)點(diǎn)，目前在航空航天、風(fēng)力發(fā)電等工業(yè)領(lǐng)域得到了廣泛的應(yīng)用[2-3]。盡管如此，復(fù)合材料在生產(chǎn)和使用過(guò)程中仍然會(huì)產(chǎn)生缺陷。復(fù)合材料的缺陷檢測(cè)方法有多種，例如超聲檢測(cè)、射線(xiàn)檢測(cè)、紅外熱波檢測(cè)等，但每種方法都有其特定的應(yīng)用范圍和優(yōu)缺點(diǎn)[4]。對(duì)于可導(dǎo)電的金屬或非金屬?gòu)?fù)合材料(如碳纖維復(fù)合材料)，可以采用電渦流法進(jìn)行檢測(cè)。渦流檢測(cè)是基于載流線(xiàn)圈與導(dǎo)體間的電磁感應(yīng)原理實(shí)現(xiàn)的，由于無(wú)接觸、自動(dòng)化程度高、速度快及無(wú)需表面處理等優(yōu)點(diǎn)而得到了廣泛的應(yīng)用[5]。

碳纖維復(fù)合材料由碳纖維和樹(shù)脂基體復(fù)合而成，表現(xiàn)為多相性和各向異性，因此無(wú)法從以往的經(jīng)驗(yàn)推斷出其中的渦流分布和線(xiàn)圈阻抗變化規(guī)律，必須通過(guò)測(cè)試分析結(jié)合數(shù)值計(jì)算來(lái)研究。近年來(lái)，國(guó)內(nèi)外學(xué)者在復(fù)合材料電學(xué)特性分析[6-8]、數(shù)值分析[9-12]和渦流檢測(cè)[13-16]方面開(kāi)展了大量的研究工作，取得了豐富的研究成果。但這些研究仍然存在很多問(wèn)題，表現(xiàn)為研究對(duì)象單一和不夠深入，主要集中于驗(yàn)證性的實(shí)驗(yàn)觀(guān)察以及怎樣提高信號(hào)的分辨率上，獲得的信息有限，沒(méi)有突出材料和缺陷的特點(diǎn)。在電學(xué)特性和數(shù)值分析方面只針對(duì)復(fù)合材料的渦流或?qū)щ娦赃M(jìn)行分析，沒(méi)有將二者結(jié)合，特別是層合板的方向性和界面對(duì)渦流分布和線(xiàn)圈阻抗的影響沒(méi)有得到深入研究。

基于以上背景，筆者以CFRP為對(duì)象，根據(jù)其層合結(jié)構(gòu)和電各向異性，利用有限元方法對(duì)渦流場(chǎng)與電學(xué)方向性的關(guān)系進(jìn)行了深入系統(tǒng)的分析，提出一種基于A-Φ的各向異性電磁有限元數(shù)值分析方法，以便有效地處理3階電磁場(chǎng)量的計(jì)算和層間渦流的連續(xù)性問(wèn)題。利用COMSOL軟件使方法實(shí)現(xiàn)并建立相應(yīng)的仿真模型，計(jì)算得到渦流密度在平面和厚度方向的分布規(guī)律以及材料中的裂紋缺陷對(duì)渦流分布和線(xiàn)圈阻抗的影響，并通過(guò)實(shí)驗(yàn)驗(yàn)證了模擬分析的有效性。

1 復(fù)合材料中電磁渦流場(chǎng)數(shù)值模型

1.1 各向異性電磁特性參數(shù)

碳纖維復(fù)合材料由多個(gè)碳纖維/樹(shù)脂片層按特定方向順序疊合而成，每層的電導(dǎo)率呈現(xiàn)出正交各向異性，沿纖維方向電導(dǎo)率高，垂直纖維方向電導(dǎo)率低，厚度方向最小。為了簡(jiǎn)化仿真分析，在數(shù)值模型中，應(yīng)將復(fù)合材料的單個(gè)片層等效成均勻連續(xù)材料，采用3階張量矩陣表示電導(dǎo)率σ和磁導(dǎo)率μ。電流密度J和電場(chǎng)強(qiáng)度E之間的關(guān)系[17]如下

J=σE

(1)

引入坐標(biāo)變換矩陣G

(2)

(3)

(4)

其中：σL，σT分別為CFRP沿纖維方向和垂直纖維方向的電導(dǎo)率；σCP為層間電導(dǎo)率，通常為σT的一半；θ為纖維角度。

(5)

1.2 復(fù)合材料電磁場(chǎng)求解方程

為了對(duì)碳纖維復(fù)合材料中的渦流場(chǎng)進(jìn)行有限元分析，需要列出電磁場(chǎng)控制方程。電磁場(chǎng)控制方程的基礎(chǔ)是麥克斯韋方程組，但是直接將麥克斯韋方程用于數(shù)值計(jì)算會(huì)造成未知量過(guò)多、方程組過(guò)于龐大的狀況。為了減少未知數(shù)的個(gè)數(shù)以減少計(jì)算量，引入矢量磁位A和標(biāo)量電位Φ作為控制方程中的未知數(shù)

B=×A

(6)

E=-?A/?t-Φ

(7)

其中：B和E分別為磁感應(yīng)強(qiáng)度和電場(chǎng)強(qiáng)度。

因此，筆者基于A-Φ列出電磁場(chǎng)控制方程，將研究的場(chǎng)域V分成渦流區(qū)V1和非渦流區(qū)V2兩部分，V12為V1和V2的交界面，如圖1所示。在渦流區(qū)，電場(chǎng)和磁場(chǎng)都需要描述；在非渦流區(qū)，只需要描述磁場(chǎng)。

圖1 渦流檢測(cè)仿真求解域Fig.1 Solution domain of eddy current testing simulation

基于A-Φ電磁位對(duì)法的渦流場(chǎng)計(jì)算控制方程[18]如下。

在V1內(nèi)

(8)

在V2內(nèi)

×(μ-1×A)-(μ-1·A)=Js

(9)

同時(shí)，標(biāo)量電位Φ的引入能使兩種導(dǎo)電介質(zhì)交界面處的電流密度連續(xù)性得到滿(mǎn)足，保證在復(fù)材板不同方向?qū)拥慕唤缟?，電流密度的法向分量J1n和J2n連續(xù)。

圖2 復(fù)合材料層合板結(jié)構(gòu)Fig.2 Structural outline of laminates

如圖2所示，將復(fù)合材料界面兩側(cè)電流密度的法向分量用矢量磁位和標(biāo)量電位表示

(10)

(11)

在庫(kù)侖規(guī)范條件下，為了滿(mǎn)足

J1n=J2n

(12)

必須有

(13)

由上式可見(jiàn)，Φ的梯度在交界面上有躍變，從而保證了電流的連續(xù)性。

2 基于COMSOL軟件的渦流場(chǎng)仿真

2.1 電渦流檢測(cè)仿真模型的建立

利用COMSOL仿真軟件的AC/DC模塊建立電渦流檢測(cè)(eddy current testing, 簡(jiǎn)稱(chēng)ECT)有限元仿真模型，并在軟件中設(shè)置基于A-Φ電磁位對(duì)的渦流場(chǎng)控制方程進(jìn)行計(jì)算，如圖3所示。

圖3 ECT仿真模型示意圖Fig.3 Schematic diagram of ECT simulation model

ECT仿真模型由空心的扁平線(xiàn)圈(內(nèi)徑為1.2mm、外徑為3.2mm、高為0.8mm、匝數(shù)為140)和CFRP板組成，并在板模型中設(shè)置裂紋區(qū)域，裂紋中心位于x=0，y=0處。表1列出了模型參數(shù)的詳細(xì)信息。

表1 ECT模型參數(shù)

2.2 復(fù)合材料中渦流密度計(jì)算

首先利用COMSOL軟件，分別對(duì)無(wú)裂紋單向鋪層(0o)8和正交鋪層(0o/90o)4的CFRP板中的渦流分布進(jìn)行了仿真分析。根據(jù)電導(dǎo)率測(cè)量結(jié)果[19]， CFRP板中單向片層的3階電導(dǎo)率可以按照σL=2 000S/m，σT=200S/m，σCP=100S/m進(jìn)行賦值；正交鋪層的CFRP板奇數(shù)片層的電導(dǎo)率按σL=2 000S/m，σT=200S/m，σCP=100S/m進(jìn)行賦值；偶數(shù)片層中電導(dǎo)率按照σL=200S/m，σT=20 000S/m，σCP=100S/m進(jìn)行賦值。仿真結(jié)果如圖4所示，其中1～3，4～6分別為單向復(fù)合板和正交復(fù)合板沿x，y，z軸方向的感應(yīng)電流密度分布。x方向?yàn)?°方向，y方向?yàn)?0°方向，z方向?yàn)榘搴穹较颉?/p>

圖4 單向和正交復(fù)合板沿x，y，z方向的感應(yīng)電流密度Fig.4 Induced current density along x, y, z directions in unidirectional and orthogonal composite plates

由于電流總是沿著電導(dǎo)率較大的方向傳播，觀(guān)察圖4(a～c)，單向鋪層的CFRP板上的感應(yīng)電流密度沿x軸方向(纖維方向)被拉長(zhǎng)，沿y和z軸方向(垂直纖維和厚度方向)被壓縮；觀(guān)察圖4(d～f)，正交鋪層的復(fù)合材料板上的感應(yīng)電流密度沿x軸方向(纖維方向)和沿y軸方向(垂直纖維方向)分布是相似的，比較均勻。因此，從仿真結(jié)果可以得到，復(fù)合材料鋪層方式的不同會(huì)導(dǎo)致各個(gè)方向上導(dǎo)電性的不同，呈現(xiàn)出電各向異性。

3 裂紋和頻率對(duì)線(xiàn)圈阻抗的影響

在渦流檢測(cè)中被測(cè)導(dǎo)體的導(dǎo)電性大小和損傷情況主要通過(guò)線(xiàn)圈的阻抗變化表現(xiàn)出來(lái)，線(xiàn)圈的阻抗包含了很多被測(cè)對(duì)象的信息。因此，通過(guò)不同裂紋尺寸、不同線(xiàn)圈輸入電流頻率對(duì)線(xiàn)圈阻抗造成的影響對(duì)渦流方向性檢測(cè)進(jìn)行了半定量分析。

3.1 裂紋缺陷對(duì)線(xiàn)圈阻抗的影響

首先，根據(jù)2.2節(jié)的方法對(duì)CFRP板進(jìn)行渦流檢測(cè)仿真計(jì)算，觀(guān)察裂紋對(duì)復(fù)合材料渦流場(chǎng)分布的影響，如圖5所示，圖(a)、圖(b)分別為裂紋單向板和正交板渦流密度幅值圖。

對(duì)圖5進(jìn)行分析，當(dāng)CFRP板中存在裂紋時(shí)，線(xiàn)圈產(chǎn)生的磁場(chǎng)將不能在缺陷區(qū)域感應(yīng)出電渦流，因此該區(qū)域的渦流密度值為0。同時(shí)，由于正交板的整體導(dǎo)電性?xún)?yōu)于單向板，因此正交板中的裂紋引起的渦流變化更加明顯。

圖5 裂紋單向和正交板在x=0處的渦流密度幅值圖Fig.5 Image of eddy density amplitude at x=0 at cracked unidirectional and orthogonal plates

接著，對(duì)模型中的線(xiàn)圈探頭設(shè)置運(yùn)動(dòng)屬性，分別沿x軸方向和y軸方向?qū)o(wú)裂紋和有裂紋(深度為100%，50%)單向板和正交板進(jìn)行方向性?huà)呙?。提取從x=-10～10mm、探頭每移動(dòng)1mm步長(zhǎng)時(shí)線(xiàn)圈的阻抗值，并對(duì)測(cè)量無(wú)裂紋和有裂紋板得到的阻抗值進(jìn)行差分運(yùn)算，繪制出如圖6所示的線(xiàn)圈阻抗差值模曲線(xiàn)。圖6(a，b)和(c，d)分別為沿x，y軸方向測(cè)量單向和正交復(fù)合材料板所得數(shù)據(jù)繪制的曲線(xiàn)圖。

圖6 線(xiàn)圈探頭阻抗差值模Fig.6 The absolute value of the impedance difference of the coil

對(duì)圖6進(jìn)行分析，可知不同的裂紋深度對(duì)得到的線(xiàn)圈阻抗值的影響不同，與50%深度裂紋相比，100%裂紋對(duì)線(xiàn)圈阻抗值的影響更大。同時(shí)，由于正交板的鋪層為(0o/90o)4方式，整體導(dǎo)電性要優(yōu)于單向復(fù)合材料板，對(duì)正交板測(cè)量得到的阻抗差值模要遠(yuǎn)遠(yuǎn)大于單向板。

3.2 頻率對(duì)線(xiàn)圈阻抗的影響

探頭線(xiàn)圈由銅導(dǎo)線(xiàn)繞制而成，因此有電阻R存在，向線(xiàn)圈中通入角頻率為ω的交變電流，線(xiàn)圈本身會(huì)產(chǎn)生電磁場(chǎng)，從而磁通量發(fā)生變化，這將改變線(xiàn)圈的特性并產(chǎn)生另外一種電流阻力，稱(chēng)為感抗X，所以渦流線(xiàn)圈阻抗Z為

Z=R+jX

(14)

X=ωL

(15)

其中：L為線(xiàn)圈電感。

依次取輸入電流頻率為200，500，1 000，2 000kHz，并計(jì)算CFRP板無(wú)裂紋和有裂紋(100%，50%)時(shí)線(xiàn)圈在x=0處的阻抗幅值及相位角，結(jié)果如圖7所示。圖7(a，b)和(c，d)分別為對(duì)單向板和正交板測(cè)量得到的試件板無(wú)缺陷和有缺陷(100%，50%)時(shí)線(xiàn)圈阻抗的幅值差值折線(xiàn)圖和相位折線(xiàn)圖。

圖7 電流頻率對(duì)渦流線(xiàn)圈阻抗幅值和相位的影響Fig.7 Influence of current frequency on impedance amplitude and phase of eddy current coils

由圖7(a，b)可以看出，隨著頻率的增加，線(xiàn)圈阻抗差值隨之增加，100%裂紋時(shí)線(xiàn)圈阻抗值比50%時(shí)大。對(duì)圖8(c，d)數(shù)據(jù)進(jìn)行分析，渦流線(xiàn)圈的相位角表示成

φ=arctan(ωL/R)

(16)

隨著輸入電流頻率的增加，相位角隨之減小。相位滯后是渦流檢測(cè)中一個(gè)非常重要的現(xiàn)象，表現(xiàn)為內(nèi)部損傷引起的線(xiàn)圈阻抗變化和表面損傷引起的阻抗變化相比，在時(shí)域上呈現(xiàn)一個(gè)滯后。頻率相同時(shí)，檢測(cè)100%裂紋板時(shí)線(xiàn)圈阻抗相位角比50%裂紋板的相位角大。因此，可以利用相位滯后現(xiàn)象獲得損傷的深度信息，估算出損傷的確切深度。

4 實(shí)驗(yàn)驗(yàn)證

為了對(duì)碳纖維復(fù)合材料渦流檢測(cè)方向性數(shù)值分析的有效性進(jìn)行驗(yàn)證，筆者基于ECT二維平面自動(dòng)掃查系統(tǒng)進(jìn)行了實(shí)驗(yàn)測(cè)量，如圖8所示。

圖8 ECT自動(dòng)掃查系統(tǒng)Fig.8 ECT automatic scanning system

圖9 復(fù)合板掃描成像Fig.9 Composite plate scanning imaging

首先，對(duì)CFRP板進(jìn)行二維平面掃查，選擇空心圓柱形線(xiàn)圈探頭，激勵(lì)頻率為2MHz。對(duì)單向鋪層和正交鋪層的CFRP板掃查，圖9(a，b)分別為對(duì)單向板和正交板的掃描成像結(jié)果，像素點(diǎn)強(qiáng)度用電壓幅值V表示。

從圖像中的條紋可清楚分辨出CFRP板中纖維的鋪層形式，且對(duì)正交板檢測(cè)時(shí)得到的信號(hào)幅值大于單向板。實(shí)驗(yàn)結(jié)果與仿真中板的感應(yīng)電流密度分布相對(duì)應(yīng)，符合沿纖維方向電導(dǎo)率大、正交板的整體導(dǎo)電性要優(yōu)于單向板的規(guī)律。

接著，制備出兩塊帶裂紋的正交板，裂紋長(zhǎng)度分別為10和20mm，板厚為2mm。利用圓片銑刀在板表面加工出50%和100%深度的狹縫裂紋，寬度約為0.2mm，如圖10所示。

圖10 裂紋損傷復(fù)合材料板Fig.10 Cracked composite plate

利用二維平面自動(dòng)掃查系統(tǒng)對(duì)裂紋板進(jìn)行掃查，結(jié)果如圖11所示。圖11(a，b)分別為對(duì)10和20mm裂紋成像結(jié)果。由圖11可以看出裂紋長(zhǎng)度的差異，50%裂紋位于上部，信號(hào)值較小。

圖11 10mm和20mm裂紋正交板渦流檢測(cè)Fig.11 Eddy current detection of 10 mm and 20 mm cracked orthogonal Plates

圖12 線(xiàn)性?huà)呙栊盘?hào)幅值Fig.12 Amplitude of linear scanning signal

圖13 分層板和沖擊板的渦流檢測(cè)Fig.13 Eddy current detection of layered and impact plate

最后，分別提取圖11(b)上對(duì)應(yīng)y=15mm與y=35mm處沿x軸直線(xiàn)掃描信號(hào)幅值，如圖12(a)所示；再提取圖11(b)上對(duì)應(yīng)x=30mm處沿y軸直線(xiàn)掃描信號(hào)幅值，如圖12(b)所示。將圖12(a，b)與圖6(c，d)對(duì)比可以看出，實(shí)驗(yàn)數(shù)據(jù)與仿真數(shù)據(jù)所繪曲線(xiàn)形狀和變化趨勢(shì)的吻合度很高。

為了進(jìn)一步說(shuō)明渦流法的可靠性，分別對(duì)兩塊四方向CFRP板中分層、沖擊損傷進(jìn)行了渦流檢測(cè)，測(cè)量結(jié)果如圖13所示，層合板方向角為(±45°，0°和90°，厚度為2mm。分層是在制備的時(shí)候離上表面0.5mm的深度上嵌入一張20mm×20mm的塑料薄膜，沖擊損傷是一個(gè)1kg的錘頭從2m高處自由落體墜落到板上形成的。由圖13可以看出，分層和沖擊損傷區(qū)域由于存在層間分離和纖維開(kāi)裂，與周?chē)鸁o(wú)損傷區(qū)域相比差異明顯，顏色較深，信號(hào)值較大。

以上實(shí)驗(yàn)結(jié)果充分說(shuō)明了渦流法用于碳纖維復(fù)合材料的無(wú)損探傷具有很好的效果。同時(shí)，驗(yàn)證了本研究對(duì)于各向異性碳纖維復(fù)合材料渦流檢測(cè)方向性數(shù)值分析的有效性。

5 結(jié) 論

1) 各向異性復(fù)合材料中的電渦流沿著導(dǎo)電性大的方向被拉伸，在導(dǎo)電性小的方向被壓縮，根據(jù)渦流分布可對(duì)復(fù)合材料的鋪層方式進(jìn)行分析，為復(fù)合材料中的缺陷定位提供重要的支撐數(shù)據(jù)。

2) 裂紋缺陷會(huì)引起復(fù)合材料中渦流的畸變，不同鋪層方式的復(fù)材板中渦流畸變形式不同，這對(duì)研究渦流檢測(cè)方向性具有重要的指導(dǎo)作用。

3) 復(fù)合材料的導(dǎo)電性越好，電渦流法對(duì)其檢測(cè)的精度就越高。正交鋪層的CFRP板的檢測(cè)精度優(yōu)于單向鋪層板。

4) 復(fù)材板和裂紋參數(shù)、輸入電流頻率會(huì)對(duì)缺陷檢測(cè)的精度造成不同程度的影響，這對(duì)研究渦流檢測(cè)的結(jié)果分析和技術(shù)優(yōu)化具有很大的幫助。