單向碳纖維增強(qiáng)復(fù)合材料的鎖相渦流熱成像檢測(cè)盲區(qū)

，，，

(西南交通大學(xué) 物理科學(xué)與技術(shù)學(xué)院 光電工程研究所與無損檢測(cè)研究中心，成都 610031)

單向碳纖維增強(qiáng)復(fù)合材料的鎖相渦流熱成像檢測(cè)盲區(qū)

郭建光，高曉蓉，郭建強(qiáng)，羅林

(西南交通大學(xué) 物理科學(xué)與技術(shù)學(xué)院 光電工程研究所與無損檢測(cè)研究中心，成都 610031)

以單向碳纖維復(fù)合材料為研究對(duì)象，使用鎖相渦流熱成像技術(shù)對(duì)試件加熱并記錄了紅外熱像圖，獲得具有較高信噪比的幅值和相位圖；對(duì)渦流加熱過程建立了有限元模型并進(jìn)行了仿真分析。試驗(yàn)結(jié)果表明：?jiǎn)蜗蛱祭w維復(fù)合材料會(huì)產(chǎn)生加熱盲區(qū)，渦流加熱存在非均勻性；仿真結(jié)果與試驗(yàn)結(jié)果有較高的一致性，渦流向量中與碳纖維垂直的分量直接導(dǎo)致了渦流加熱的非均勻性。檢測(cè)盲區(qū)的發(fā)現(xiàn)，有助于指導(dǎo)碳纖維增強(qiáng)復(fù)合材料的渦流熱成像無損檢測(cè)。

鎖相熱成像；碳纖維增強(qiáng)復(fù)合材料；渦流加熱；有限元仿真；各向異性

碳纖維增強(qiáng)復(fù)合材料是高強(qiáng)度、低密度的專用材料,被廣泛應(yīng)用于航空航天的關(guān)鍵部件中。隨著工藝改善和產(chǎn)能的提升，碳纖維增強(qiáng)復(fù)合材料(CFRP)開始在日常生活中逐漸普及，其從基礎(chǔ)材料到各類制品的制造都需要依賴特殊的焊接工藝。與金屬型材的焊接不同，碳纖維增強(qiáng)復(fù)合材料的焊接是使用渦流加熱的方法使碳纖維材料升溫，通過外加壓力使獨(dú)立的板材固化成整體的。通過在材料內(nèi)部形成渦流來加熱復(fù)合材料，因此渦流加熱是否均勻，就成為影響焊接質(zhì)量的另一個(gè)重要因素。經(jīng)過研究，渦流加熱的效果與電場(chǎng)激勵(lì)頻率、材料屬性和結(jié)構(gòu)都有關(guān)[1]，相關(guān)的加熱理論和試驗(yàn)已經(jīng)有所研究[2-3]。渦流熱成像又稱為感應(yīng)熱成像，是一種使用感應(yīng)電流來加熱被檢工件的主動(dòng)熱成像方式，被廣泛地應(yīng)用于金屬表面裂紋的無損檢測(cè)中[4]。同時(shí)，使用渦流熱成像技術(shù)研究碳纖維復(fù)合材料缺陷的工作已有開展[5-8]。筆者利用紅外相機(jī)和渦流加熱裝置對(duì)碳纖維增強(qiáng)復(fù)合材料試件進(jìn)行鎖相熱成像檢測(cè)，記錄試驗(yàn)過程中材料溫度的分布，再使用有限元仿真方法進(jìn)行論證和研究。

1 鎖相渦流熱成像的原理

鎖相渦流熱成像包含兩個(gè)基本原理，一是渦流加熱原理，二是鎖相熱成像原理。渦流加熱方法常用于金屬材料的熱處理中，其基本原理是：將待加熱物體放置在感應(yīng)線圈產(chǎn)生的交變磁場(chǎng)中，法拉第電磁感應(yīng)定律的存在使得導(dǎo)體材料中產(chǎn)生了渦流，渦流的作用直接導(dǎo)致焦耳熱的產(chǎn)生，繼而實(shí)現(xiàn)加熱的目的。

鎖相熱成像的原理：使用紅外相機(jī)采集被檢物體產(chǎn)生的輻射熱，形成二維溫度云圖進(jìn)行檢測(cè)。熱成像方法分為兩類，一類是被動(dòng)熱成像，另一類是主動(dòng)熱成像(見圖1)。被動(dòng)熱成像實(shí)施過程不需要外部激勵(lì)，被廣泛用于物體溫度的動(dòng)態(tài)監(jiān)控；主動(dòng)熱成像需要外部激勵(lì)，在工業(yè)領(lǐng)域主要用于無損檢測(cè)，熱量在檢測(cè)對(duì)象內(nèi)部擴(kuò)散的同時(shí)還會(huì)受到缺陷的阻擋，導(dǎo)致溫度分布不均勻，從而表征檢測(cè)對(duì)象的表面或表面以下的缺陷。

圖1 熱成像方法分類

主動(dòng)熱成像有兩種外部激勵(lì)方式，一種是脈沖式，即單次瞬時(shí)加熱(見圖2)，物體溫度上升和下降的規(guī)律受到結(jié)構(gòu)和材料完整性的影響，從而使得缺陷得以表征。

圖2 脈沖式加熱的溫度變化規(guī)律

另一種激勵(lì)方式是鎖相式，也是文中所使用的熱成像方法，即施加周期性強(qiáng)度變化的激勵(lì)。鎖相方法中使用低頻調(diào)制高頻交變電流信號(hào)原理示意如圖3所示。研究使用的鎖相渦流熱成像方式，就是將高頻激勵(lì)電流信號(hào)(見圖3左上)用低頻調(diào)制信號(hào)(見圖3右上)調(diào)制形成激勵(lì)信號(hào)(見圖3下)，使材料中的感應(yīng)渦流也呈現(xiàn)周期性的幅值變化。

圖3 鎖相方法中使用低頻調(diào)制高頻交變電流信號(hào)原理示意

同時(shí)，紅外相機(jī)采集形成熱圖序列，假設(shè)熱圖數(shù)量為N。采集時(shí)相機(jī)與物體保持靜止，因此不同熱圖同一坐標(biāo)位置上的像素對(duì)應(yīng)的是相同的空間位置，鎖相處理過程示意如圖4所示。隨后對(duì)整個(gè)熱圖序列進(jìn)行像素級(jí)快速傅里葉變換得到幅值和相位圖。

圖4 鎖相處理過程示意

采集到的每張熱圖里的每一個(gè)像素點(diǎn)，都是紅外相機(jī)感應(yīng)器件所能采集到的最小溫度信號(hào)單元和進(jìn)行鎖相處理的最小單元，溫度信號(hào)F(t)主要包括以下3個(gè)部分。

式中：Ds為由被檢材料輻射率波動(dòng)和探測(cè)器非一致性影響的溫度變化信號(hào)；T(t)為材料被渦流加熱而導(dǎo)致的溫度變化信號(hào)；DR為材料被線圈輻射加熱而導(dǎo)致的溫度變化信號(hào)。

其中，只有T(t)分量是對(duì)熱圖成像和鎖相運(yùn)算起積極作用的信號(hào)，其包含了渦流與材料相互作用的信息，而Ds和DR分量都是要在后期處理算法中極力消除的。在對(duì)N個(gè)熱圖的每個(gè)空間像素點(diǎn)進(jìn)行快速傅里葉變換時(shí)，相同空間位置的N個(gè)像素點(diǎn)按時(shí)間順序組成序列，分別進(jìn)行了正、余弦信號(hào)積分處理，獲得同相信號(hào)S0°和反相信號(hào)S-90°。

在F(t)分別與正余弦信號(hào)先相乘再疊加后，獲得的同、反相信號(hào)中DR分量被削弱，F(xiàn)(t)信號(hào)實(shí)現(xiàn)第一次濾波和信噪比提升。將其繪制在復(fù)平面中，可以求得同一空間位置像素點(diǎn)的幅值和相位，即由此獲得幅值和相位圖(見圖5)。

圖5 同相信號(hào)與反相信號(hào)繪制在復(fù)平面圖中獲得的幅值和相位示意

在求解每個(gè)像素點(diǎn)的相位φ時(shí)，同相和反相信號(hào)做除法，可消除Ds分量對(duì)信噪比的影響，F(xiàn)(t)信號(hào)質(zhì)量得以進(jìn)一步提升。進(jìn)行傅里葉變換后，幅值圖反映了被檢測(cè)物體的溫度強(qiáng)度，而相位圖反映了熱量輻射到相機(jī)的時(shí)間延遲情況。由于既消除了DR分量又濾掉了Ds分量，因此相位圖比幅值圖的信噪比更強(qiáng)。在相同加熱功率的情況下，鎖相熱成像法通常比脈沖熱成像法具有更高的熱圖分辨率。也是基于這一優(yōu)點(diǎn)，筆者選擇了鎖相渦流熱成像檢測(cè)方法。

2 檢測(cè)過程

2.1檢測(cè)平臺(tái)

試驗(yàn)采用鎖相渦流熱成像檢測(cè)系統(tǒng)(見圖6)，該平臺(tái)搭建使用了惠普HP6116波形產(chǎn)生器、ENI Model 2100L功率放大器、加熱頭、非制冷銅感應(yīng)線圈、Thermosensorik QWIP 384紅外相機(jī)及標(biāo)配控制單元。試驗(yàn)進(jìn)行過程中，波形產(chǎn)生器提供可選頻率和波形的感應(yīng)頻率，高頻電信號(hào)經(jīng)過功率放大器后，可以獲得最高約200 W的加熱功率。為了保證處理后的圖像擁有較高的信噪比，145 幀·s-1的紅外相機(jī)每次采集的熱圖序列包含高達(dá)1 024張甚至2 048張圖像。采集過程結(jié)束后，控制單元運(yùn)行鎖相算法獲得一張幅值圖和一張相位圖。由于采用鎖相方法，控制單元還需要保證渦流激勵(lì)過程與熱圖采集過程的同步。鎖相頻率即調(diào)制頻率，由于線圈不同，在0.4 ～36 Hz范圍內(nèi)設(shè)定調(diào)制頻率。

圖6 鎖相渦流熱成像檢測(cè)系統(tǒng)示意

根據(jù)線圈、相機(jī)和被檢材料的相對(duì)關(guān)系，有透射和反射兩種加熱模式(見圖7)。紅外成像結(jié)果均以灰度圖表示，幅值圖中顏色越淺表示溫度越高，反之則越低。

圖7 兩種加熱模式示意

圖8 單向CFRP圓形板材外觀

2.2檢測(cè)對(duì)象

圖9 單向CFRP的纖維走向?qū)е掳宀某尸F(xiàn)各向異性示意

檢測(cè)對(duì)象是一塊單向CFRP圓形板材，直徑200 mm，厚1 mm(見圖8)，被長(zhǎng)臂120 mm、短壁35 mm的矩形線圈加熱。值得注意的是，不同于金屬各方向均一的磁導(dǎo)率和電導(dǎo)率，CFRP的電磁學(xué)和熱學(xué)性能受到碳纖維走向的約束，呈現(xiàn)出各項(xiàng)異性的特點(diǎn)(見圖9)。順纖維方向的電導(dǎo)率和磁導(dǎo)率都遠(yuǎn)大于正交方向。在碳纖維中形成的感應(yīng)渦流及熱傳遞方向均會(huì)受此影響。試驗(yàn)中每根碳纖維束寬約5 mm，由若干發(fā)絲狀纖維并行排列組成。

2.3檢測(cè)結(jié)果

單向碳纖維感應(yīng)加熱幅值如圖10所示，試驗(yàn)中觀察到異常的現(xiàn)象：在矩形線圈左下拐角處，順碳纖維走向方向上有一明顯的暗色條帶，表明此區(qū)域溫度較低，未被加熱或加熱效率差，與兩側(cè)被加熱的碳纖維呈現(xiàn)鮮明的對(duì)比，該區(qū)域稱為“加熱盲區(qū)”。

圖10 單向碳纖維感應(yīng)加熱幅值圖

為了排除材料缺陷等因素對(duì)“加熱盲區(qū)”的成因造成的混淆，采取固定線圈將CFRP圓盤在同一平面內(nèi)旋轉(zhuǎn)的方法，使纖維走向與線圈長(zhǎng)臂之間的夾角按每10°為一個(gè)步長(zhǎng)旋轉(zhuǎn)，重復(fù)若干次試驗(yàn)，得出如圖11所示的結(jié)果。

圖11 旋轉(zhuǎn)角度試驗(yàn)的紅外幅值圖

每個(gè)角度(指的是碳纖維走向與線圈長(zhǎng)臂的銳角夾角)下的紅外幅值圖中，都出現(xiàn)了“加熱盲區(qū)”條紋，隨著圓盤的旋轉(zhuǎn)，這些線型條紋的角度也同步變化。在0°時(shí)，碳纖維走向與線圈長(zhǎng)臂平行，此時(shí)可以觀察到上下對(duì)稱的兩條暗紋，暗紋平行且位于長(zhǎng)臂附近。需要注意的是圖11中始終垂直于線圈長(zhǎng)臂的亮線(標(biāo)注在40°結(jié)果圖中)是試驗(yàn)中的參考線，并不是CFRP檢測(cè)對(duì)象的一部分，僅用于作為旋轉(zhuǎn)角度的參考。

3 有限元仿真

對(duì)單層單向碳纖維復(fù)合材料(單向 CFRP)的渦流加熱情景進(jìn)行了有限元仿真。將物理場(chǎng)劃分為大量的微小單元，通過求解每個(gè)單元的簡(jiǎn)單解，去逼近無限未知量的真實(shí)系統(tǒng)。使用COMSOL Multiphysics 5.1 軟件對(duì)試驗(yàn)過程進(jìn)行了合理簡(jiǎn)化(如將矩形線圈簡(jiǎn)化為圓形線圈)，仿真結(jié)果與試驗(yàn)結(jié)果有較高的一致性。

圖12 CFRP的三維感應(yīng)加熱模型

CFRP的三維感應(yīng)加熱模型如圖12所示(線圈半徑30 mm，與板材平行且懸于板材上方4 mm處)，需要建立的有限元模型共有3個(gè)域，分別是空氣立方體和位于其內(nèi)部的線圈及單向CFRP板材。單向CFRP板材設(shè)計(jì)尺寸(長(zhǎng)×寬×厚)為100 mm×100 mm×1 mm，感應(yīng)線圈設(shè)計(jì)為直徑60 mm單匝圓環(huán)，懸于板材正中央上方4 mm處。

CFRP各向異性的特點(diǎn)可以通過配置材料屬性來實(shí)現(xiàn)。根據(jù)真實(shí)的碳纖維復(fù)合材料電、磁、熱及機(jī)械屬性參數(shù)進(jìn)行模擬，單向CFRP的參數(shù)如表1所示。

表1 單向CFRP的參數(shù)

表格中，x，y，z方向的物理參數(shù)不盡相同，其中x方向?yàn)樘祭w維的延伸方向，該方向的熱導(dǎo)率和電導(dǎo)率都比其他兩個(gè)方向的大，這表征了CFRP各向異性的特點(diǎn)，較為客觀地反映了單向CFRP材料的物理性質(zhì)。模擬加熱時(shí)交變電場(chǎng)頻率為100 kHz。

3.1仿真結(jié)果與試驗(yàn)結(jié)果的對(duì)比

仿真結(jié)果展示在圖13(a)的灰度圖中，顏色由亮到暗表示溫度由高到低，有上下兩條暗色帶，互相對(duì)稱且同樣順碳纖維走向分布。將矩形線圈替換成圓線圈，重復(fù)先前板材試驗(yàn)，得到如圖13(b)，13(c)所示的幅值圖和相位圖。在相位圖中，可以比較清晰地看出試驗(yàn)結(jié)果與仿真結(jié)果有較好的一致性，該試驗(yàn)結(jié)果同樣印證了相位圖比幅值圖具有更優(yōu)的噪聲抑制功能。

圖13 圓形線圈對(duì)單向CFRP渦流加熱的仿真溫度圖與試驗(yàn)結(jié)果的對(duì)比

圖14 表面以下0.1 mm截面渦流密度分布

3.2從渦流角度分析盲區(qū)產(chǎn)生的原因

在仿真結(jié)果與試驗(yàn)結(jié)果有較好一致性的基礎(chǔ)上，從渦流分布的角度分析出現(xiàn)盲區(qū)的原因。觀察板材表層以下0.1 mm處的x,y,z三個(gè)方向分量的渦流分布。表面以下0.1 mm截面渦流密度分布如圖14所示，從圖15的仿真結(jié)果來看，x方向分量的渦流密度范圍為1 000～2 000 A·m-3，y方向分量的渦流密度范圍為20～100 A·m-3，z方向分量的渦流密度范圍為20～40 A·m-3。渦流密度的y分量與纖維走向垂直，電導(dǎo)率極低，但y分量的渦流密度分布規(guī)律與加熱盲區(qū)的圖案十分契合，兩者都顯示出線圈上下兩側(cè)是順纖維方向的。

圖15 溫度分布與渦流密度y分量的對(duì)比

這一現(xiàn)象表明與纖維方向垂直的渦流分量是熱圖圖案形成的主要貢獻(xiàn)者，該方向上纖維電導(dǎo)率低，卻生成了較高的熱量，說明渦流產(chǎn)生的焦耳熱并沒有占據(jù)主體地位。在碳纖維復(fù)合材料的感應(yīng)加熱過程中，除了焦耳熱(即電阻損耗)之外，還會(huì)產(chǎn)生電解質(zhì)損耗，感應(yīng)加熱的整體包括這兩部分。碳纖維復(fù)合材料的微觀電氣連接和電介質(zhì)關(guān)系如圖16所示，碳纖維等效為電阻，相鄰碳纖維之間相互接觸的為直接電氣連接，其他非接觸位置之間充滿高分子塑料，因此可以等效為電容，電阻損耗和電介質(zhì)損耗分別發(fā)生在碳纖維內(nèi)和高分子塑料內(nèi)。 該微觀模型由MOOK于1994年提出[9]，可以作為解釋纖維正交方向的渦流導(dǎo)致更強(qiáng)的熱量產(chǎn)生的參考。

圖16 碳纖維復(fù)合材料的微觀電氣連接和電介質(zhì)關(guān)系示意

4 結(jié)語

通過試驗(yàn)和仿真，發(fā)現(xiàn)并驗(yàn)證了單向碳纖維復(fù)合材料渦流檢測(cè)存在加熱盲區(qū)的現(xiàn)象，這種“加熱盲區(qū)”在以往文獻(xiàn)中并沒有提及。由于“盲區(qū)”的形狀和發(fā)生位置十分規(guī)律，推測(cè)應(yīng)是該處的渦流場(chǎng)產(chǎn)生了抵消，進(jìn)而導(dǎo)致渦流體密度遠(yuǎn)低于周圍的纖維束，導(dǎo)致無法加熱到相當(dāng)?shù)臏囟榷诜祱D中呈現(xiàn)暗色，因此可以推斷這是加熱原理導(dǎo)致的“盲區(qū)”而非材料本身的缺陷。該盲區(qū)的發(fā)現(xiàn)，對(duì)于有預(yù)見性地規(guī)避加熱盲區(qū)、提升復(fù)合材料加熱及焊接工藝、指導(dǎo)CFRP的渦流熱成像檢測(cè)有重要意義。

致謝本研究的試驗(yàn)和仿真過程是在德國弗朗霍夫無損檢測(cè)研究所進(jìn)行的，感謝NETZELMANN博士給予的悉心指導(dǎo)，同時(shí)還要感謝我的國內(nèi)導(dǎo)師高曉蓉教授對(duì)我研究工作的大力支持。

[1] YARLAGADDA S, KIM H J, GILLESPIE J W,et al. A study on the induction heating of conductive fiber reinforced composites[J]. Journal of Composite Materials,2002,36:401-421.

[2] MILLER A K, CHANG C, PAYNE A,et al. The nature of induction heating in graphite-fiber, polymer-matrix composite materials[J]. Sampe Journal,1990,26 (4):37-54.

[3] FINK B K, MCCULLOUGH R L, GILLESPIE J J. A local theory of heating in cross-ply carbon fiber thermoplastic composites by magnetic induction[J].Polymer Engineering and Science,2010,32(5):357-369.

[4] NETZELMANN U, WALLE G, LUGIN S,et al. Induction thermography: principle, applications and first steps towards standardization[C]//Quantitative Infrared Thermography Asia.[S.l.]:[s.n],2016.

[5] REN W, LIU J, TIAN G Y, et al. Quantitative non-destructive evaluation method for impact damage using eddy current pulsed thermography[J].Composites Part B Engineering,2013,54(1):169-179.

[6] CHENG L, TIAN G Y. Surface crack detection for carbon fiber reinforced plastic (CFRP) materials using pulsed eddy current thermography[J]. 2012,11(12):3261-3268.

[7] HE Y, TIAN G, PAN M, et al. Impact evaluation in carbon fiber reinforced plastic (CFRP) laminates using eddy current pulsed thermography[J]. Composite Structures, 2014, 109(1):1-7.

[8] HE Y, TIAN G Y, PAN M, et al. Non-destructive testing of low-energy impact in CFRP laminates and interior defects in honeycomb sandwich using scanning pulsed eddy current[J]. Composites Part B Engineering, 2014, 59(3):196-203.

[9] LANGE R, MOOK G. Structural analysis of CFRP using eddy current methods[J]. NDT & E International, 1994, 27(5): 241-248.

OntheDeadZoneofEddyCurrentHeatingofLock-inThermographyforUnidirectionalCarbonFiberReinforcedPlastic

GUO Jianguang, GAO Xiaorong, GUO Jianqiang, LUO Lin

(Photoelectric Engineering Institute and Nondestructive Testing Research Center, School of Physical Science and Technology,Southwest Jiaotong University, Chengdu 610031, China)

The unidirectional carbon fiber composite materials (UD CFRP) was taken as the research object and the lock-in eddy current thermography was utilized to heat the specimen, the high SNR of amplitude and phase diagram were calculated. Finite element model of the eddy current heating was established. Experimental results show that the UD CFRP will lead to dead zone of induction heating. Simulation and experimental results have a higher consistency and reveal that components of the eddy current vector which is perpendicular to direction of carbon fiber directly lead to non-uniformity of induction heating. This will also benefit NDT for CFRP avoiding inspection dead zone.

lock-in thermography; CFRP; eddy current heating;FEM simulation; anisotropy

2017-04-17

郭建光(1991-)，男，碩士，主要從事軌道車輛無損檢測(cè)設(shè)備和技術(shù)的研發(fā)工作

高曉蓉,gxrr@vip.163.com

10.11973/wsjc201711001

TG115.28

A

1000-6656(2017)11-0001-06