(四川大學(xué) 制造科學(xué)與工程學(xué)院，成都 610065)

工業(yè)零部件在生產(chǎn)、加工及使用過(guò)程中，會(huì)受到外部載荷的作用，從而在材料的局部區(qū)域產(chǎn)生應(yīng)力集中，致使材料的力學(xué)性能降低，甚至導(dǎo)致材料失效。因此，對(duì)材料或結(jié)構(gòu)的應(yīng)力狀態(tài)進(jìn)行檢測(cè)與監(jiān)測(cè)，特別是采用無(wú)損檢測(cè)方法進(jìn)行檢測(cè)與監(jiān)測(cè)是很有必要的。在眾多的無(wú)損檢測(cè)方法中，渦流熱成像檢測(cè)(ECT)技術(shù)可以快速地檢測(cè)試件表面或近表面的缺陷，如表面裂紋[1-3]、腐蝕缺陷[4-6]、焊接缺陷[7-9]、疲勞損傷[10-12]等。同超聲法、射線法等相比，ECT方法具有非接觸、單次檢測(cè)面積大、對(duì)復(fù)雜形狀零件的易適應(yīng)性等優(yōu)點(diǎn)，可用于大范圍內(nèi)的應(yīng)力先驗(yàn)檢測(cè)。近幾年，不少專(zhuān)家學(xué)者將渦流熱成像檢測(cè)技術(shù)應(yīng)用于金屬材料沖擊損傷應(yīng)力、拉伸應(yīng)力的檢測(cè)中，并取得了一定的研究成果。通過(guò)檢測(cè)多次沖擊后的鋼板，劉錄葉等證明了應(yīng)用ECT檢測(cè)沖擊殘余應(yīng)力是可行的[13]；利用ECT檢測(cè)拉伸應(yīng)力，有學(xué)者給出了材料縱向和橫向電磁特性的變化規(guī)律[14]以及應(yīng)力和熱導(dǎo)率之間的線性耦合關(guān)系[15]；雷慶等有效定量了單向拉伸載荷下鋼結(jié)構(gòu)的表面應(yīng)力[16]等。但在實(shí)際應(yīng)力檢測(cè)中，預(yù)制一定量大小及分布厚度的應(yīng)力樣本試件比較困難，結(jié)果量化也比較復(fù)雜。因此，筆者采用數(shù)值仿真分析手段，通過(guò)設(shè)置材料局部區(qū)域的電導(dǎo)率、熱導(dǎo)率等參數(shù)的大小及其厚度來(lái)模擬實(shí)際應(yīng)力層，建立渦流熱成像仿真模型并分析試件表面溫度響應(yīng)的變化情況，試圖為渦流熱成像定量檢測(cè)材料或結(jié)構(gòu)的應(yīng)力提供理論參考。

1 ECT應(yīng)力檢測(cè)原理

應(yīng)用ECT技術(shù)檢測(cè)導(dǎo)體材料的應(yīng)力分布時(shí)，含應(yīng)力區(qū)域的感應(yīng)渦流密度大小是異于正常區(qū)域的，且與應(yīng)力幅值相關(guān)，應(yīng)力幅值的改變會(huì)導(dǎo)致材料電磁特性參數(shù)和傳熱特性參數(shù)的改變，從而使得感應(yīng)焦耳熱在材料內(nèi)部傳導(dǎo)的速度發(fā)生改變，最終表現(xiàn)為熱成像圖像溫度場(chǎng)的不均勻分布和應(yīng)力區(qū)與非應(yīng)力區(qū)的溫升變化規(guī)律的不一致。通過(guò)對(duì)ECT熱成像圖像/序列的后期處理分析，可以獲得材料的應(yīng)力狀態(tài)，其檢測(cè)理論模型如圖1所示。電導(dǎo)率變化主要影響檢測(cè)過(guò)程中的感應(yīng)加熱，熱導(dǎo)率主要影響冷卻過(guò)程中的熱傳導(dǎo)。

圖1 ECT應(yīng)力檢測(cè)理論模型

1.1 檢測(cè)原理

渦流熱成像技術(shù)在檢測(cè)過(guò)程中所涉及到的物理過(guò)程有：感應(yīng)渦流加熱、熱傳導(dǎo)、紅外輻射[17]。當(dāng)激勵(lì)線圈中通有一定頻率的交變電流時(shí)，會(huì)產(chǎn)生電磁感應(yīng)現(xiàn)象，即在被測(cè)導(dǎo)體內(nèi)會(huì)感應(yīng)出相同頻率的電渦流，由麥克斯韋爾方程推導(dǎo)出該渦流場(chǎng)的控制方程如式(1)所示。

(1)

式中：μ為材料的磁導(dǎo)率；A為磁矢勢(shì)；σe為材料的電導(dǎo)率；ε為介電常數(shù)；J為外部電流密度;t為電磁激勵(lì)作用時(shí)間。

產(chǎn)生的渦流會(huì)產(chǎn)生集膚效應(yīng)，積聚在被測(cè)導(dǎo)體表面，其集膚深度dJ的計(jì)算公式如式(2)所示。

(2)

式中：f為激勵(lì)電流頻率。

由焦耳定律可知，部分渦流會(huì)在導(dǎo)體內(nèi)部由電能轉(zhuǎn)換為熱能，產(chǎn)生的焦耳熱Q正比于渦流密度Js和電場(chǎng)強(qiáng)度E，其計(jì)算公式如式(3)所示。

(3)

由式(3)可知，影響焦耳熱的因素主要有導(dǎo)體材料的電導(dǎo)率和渦流密度，但因?yàn)殡妼?dǎo)率也會(huì)影響渦流的密度值，因此焦耳熱與電導(dǎo)率不是單純的負(fù)冪次關(guān)系，其變化規(guī)律還需要進(jìn)一步探討。

在感應(yīng)加熱階段產(chǎn)生的焦耳熱Q會(huì)在被測(cè)試件內(nèi)部進(jìn)行傳遞，其傳遞規(guī)律及熱擴(kuò)散深度dT的計(jì)算公式如式(4)，(5)所示。

(4)

(5)

式中：ρ為材料密度；Cp為材料的比熱容；T為材料熱力學(xué)溫度；λ為材料的熱導(dǎo)率；α=λ/(ρCp)為導(dǎo)體熱擴(kuò)散系數(shù);t為熱擴(kuò)散的時(shí)間。

1.2 應(yīng)力對(duì)導(dǎo)體材料屬性的影響

1.2.1 金屬電阻應(yīng)變效應(yīng)

導(dǎo)體材料內(nèi)部應(yīng)力會(huì)使得材料內(nèi)部原子間的距離和點(diǎn)陣的動(dòng)畸變發(fā)生變化，導(dǎo)致材料的電阻率發(fā)生變化，從而影響材料的電導(dǎo)率[15，18]。一般情況下，電導(dǎo)率隨應(yīng)力的變化呈負(fù)指數(shù)關(guān)系，當(dāng)指數(shù)的絕對(duì)值遠(yuǎn)遠(yuǎn)小于1時(shí)，則略去高階小項(xiàng)，近似成為線性關(guān)系。其影響規(guī)律遵循式(6)[19]。

(6)

式中：σe1為導(dǎo)體存在應(yīng)力時(shí)的電導(dǎo)率；σe0為導(dǎo)體無(wú)應(yīng)力時(shí)的電導(dǎo)率；σ為應(yīng)力值(拉應(yīng)力為正、壓應(yīng)力為負(fù))；αr為應(yīng)力系數(shù)。

即當(dāng)試件存在拉應(yīng)力時(shí)，其應(yīng)力區(qū)域的電導(dǎo)率將減小相反，當(dāng)試件存在壓應(yīng)力時(shí)，其應(yīng)力區(qū)域的電導(dǎo)率將增大。

在檢測(cè)過(guò)程中，由于被測(cè)試件產(chǎn)生的焦耳熱會(huì)影響材料本身的溫度，且電導(dǎo)率會(huì)隨溫度的變化而發(fā)生變化，其關(guān)系遵循式(7)。

(7)

式中：σeT，σeT0分別為材料在溫度T，T0時(shí)的電導(dǎo)率；β為導(dǎo)體電阻溫度系數(shù)。

1.2.2 魏德曼-弗蘭茲定律

由魏德曼-弗蘭茲定律可知，在室溫下很多金屬的熱導(dǎo)率與電導(dǎo)率之比幾乎相同，如式(8)所示。因此當(dāng)被測(cè)試件內(nèi)部存在應(yīng)力時(shí)，應(yīng)力區(qū)域的熱導(dǎo)率也會(huì)隨著電導(dǎo)率的改變而發(fā)生相應(yīng)改變。

(8)

設(shè)λ0為無(wú)應(yīng)力時(shí)的熱導(dǎo)率，λ1為存在應(yīng)力時(shí)的熱導(dǎo)率。設(shè)置仿真模型參數(shù)時(shí)，同電導(dǎo)率一樣，當(dāng)試件存在拉應(yīng)力時(shí)，將試件應(yīng)力區(qū)域的熱導(dǎo)率λ1同σe1等比例減??；相反，當(dāng)試件存在壓應(yīng)力時(shí)，將λ1同σe1等比例增大。

圖2 渦流熱成像檢測(cè)的仿真模型及其網(wǎng)格劃分

2 ECT應(yīng)力檢測(cè)仿真模型

渦流熱成像檢測(cè)過(guò)程中，被測(cè)試件溫升的熱源主要來(lái)自導(dǎo)體材料的電磁損耗，在進(jìn)行數(shù)值仿真分析時(shí)，需要耦合電磁場(chǎng)和溫度場(chǎng)。文章用COMSOL MUTIPHYSICS多物理仿真軟件建立了脈沖渦流熱成像檢測(cè)的三維模型，模型及其網(wǎng)格劃分如圖2所示。選取電磁場(chǎng)與熱傳導(dǎo)兩個(gè)物理場(chǎng)進(jìn)行耦合，并將熱傳導(dǎo)模式設(shè)置為固體傳熱，通過(guò)瞬態(tài)分析獲取激勵(lì)加熱階段和自然冷卻階段的溫度場(chǎng)變化情況。激勵(lì)線圈直徑為6 mm，長(zhǎng)150 mm，距檢測(cè)表面的提離為1 mm；導(dǎo)體材料的尺寸(長(zhǎng)×寬×厚)為200 mm×120 mm×10 mm，中部存在(長(zhǎng)×寬)20 mm×20 mm的區(qū)域?yàn)殡妼?dǎo)率及熱導(dǎo)率變化層，厚度可變，模擬被測(cè)試件存在應(yīng)力的區(qū)域；初始溫度設(shè)置為293 K，激勵(lì)電流設(shè)置為500 A，加載頻率為250 kHz，加熱時(shí)間為500 ms，冷卻時(shí)間為1 000 ms，即總的檢測(cè)時(shí)間為1 500 ms。

仿真所用導(dǎo)體材料的參數(shù)見(jiàn)表1，由式(5)計(jì)算出材料正常區(qū)域在檢測(cè)結(jié)束后的熱擴(kuò)散深度可達(dá)到4.231 mm，即在研究應(yīng)力層厚度對(duì)溫度響應(yīng)曲線的影響時(shí)，最大厚度不能超過(guò)4.231 mm。

表1 仿真所用材料的參數(shù)

3 COMSOL仿真結(jié)果討論

在分析仿真結(jié)果的過(guò)程中，分別考查被測(cè)導(dǎo)體表面點(diǎn)1(位于試件應(yīng)力層中部)和點(diǎn)2(位于遠(yuǎn)離應(yīng)力層的正常區(qū)域處)的溫度響應(yīng)變化情況。如圖3所示，點(diǎn)1,點(diǎn)2的檢測(cè)值分別代表應(yīng)力層及遠(yuǎn)離應(yīng)力層的正常區(qū)域的表面溫度值，以下曲線原始數(shù)據(jù)點(diǎn)均是基于點(diǎn)1,點(diǎn)2處的溫度值所繪制。

圖3 數(shù)據(jù)取點(diǎn)位置示意

3.1 電導(dǎo)率、熱導(dǎo)率變化與檢測(cè)結(jié)果的關(guān)系

電導(dǎo)率、熱導(dǎo)率變化層即所模擬的應(yīng)力層，其厚度設(shè)置為1 mm，導(dǎo)體正常區(qū)域的電導(dǎo)率σe0設(shè)置為5.0×106S·m-1，為得到較光滑的仿真曲線，應(yīng)力區(qū)域的電導(dǎo)率σe1設(shè)置如下：從1.0×105S·m-1以1.0×105S·m-1的增量逐次增加至5.0×106S·m-1，模擬拉應(yīng)力層；再繼續(xù)從5.0×106S·m-1以1.0×105S·m-1的增量逐次增加至1.0×107S·m-1，模擬壓應(yīng)力層。其中，熱導(dǎo)率與電導(dǎo)率呈等比例變化，因此后續(xù)結(jié)果分析中只列出了電導(dǎo)率與溫度值的關(guān)系。結(jié)合仿真結(jié)果，探討上述變化過(guò)程中被測(cè)試件表面溫升曲線的響應(yīng)情況。改變層電導(dǎo)率對(duì)導(dǎo)體表面溫度分布的影響如圖4所示。

圖4 改變層電導(dǎo)率對(duì)導(dǎo)體表面溫度分布的影響

圖5 不同電導(dǎo)率對(duì)應(yīng)點(diǎn)1,點(diǎn)2處的溫度響應(yīng)曲線

由圖4可見(jiàn)，被測(cè)導(dǎo)體拉應(yīng)力層區(qū)域的溫度明顯高于標(biāo)準(zhǔn)件的，且拉應(yīng)力越大，應(yīng)力層感應(yīng)渦流越集中，加熱結(jié)束及冷卻結(jié)束時(shí)的高溫區(qū)越明顯；相反，壓應(yīng)力層區(qū)域的溫度低于標(biāo)準(zhǔn)件的，壓應(yīng)力越大，應(yīng)力層感應(yīng)渦流越分散，低溫區(qū)范圍越大。不同電導(dǎo)率對(duì)應(yīng)點(diǎn)1，點(diǎn)2處的溫度響應(yīng)曲線如圖5所示，可以看出電導(dǎo)率的改變對(duì)點(diǎn)1處的溫度影響較大，基于此可評(píng)估應(yīng)力大小的變化情況，而遠(yuǎn)離應(yīng)力區(qū)域的點(diǎn)2處溫度幾乎不受電導(dǎo)率變化影響，可以其作為拉、壓應(yīng)力層溫度的參照值。

圖6 加熱結(jié)束時(shí)電導(dǎo)率變化對(duì)點(diǎn)1、點(diǎn)2處溫度值的影響

圖6是應(yīng)力區(qū)域點(diǎn)1處與正常區(qū)域點(diǎn)2處在加熱結(jié)束時(shí)的溫度T1,T2隨電導(dǎo)率、熱導(dǎo)率變化的曲線圖，從圖中可以看出，拉應(yīng)力層電導(dǎo)率越小，應(yīng)力區(qū)溫度高于正常區(qū)且溫差越大；壓應(yīng)力層電導(dǎo)率越大，應(yīng)力區(qū)溫度低于正常區(qū)且溫差也越大。通過(guò)式(6)反演導(dǎo)體電導(dǎo)率與其應(yīng)力絕對(duì)值(不考慮拉、壓應(yīng)力的正負(fù)問(wèn)題)之間的關(guān)系，如式(9)所示，并繪制應(yīng)力大小變化與點(diǎn)1、點(diǎn)2處溫度差值(T1-T2)之間的關(guān)系曲線圖，如圖7所示。

(9)

圖7 應(yīng)力變化對(duì)溫度的影響

結(jié)果表明，在一定范圍內(nèi)，導(dǎo)體應(yīng)力值與表面應(yīng)力區(qū)域和正常區(qū)域的溫度差近似為線性關(guān)系。拉應(yīng)力越大，應(yīng)力區(qū)與正常區(qū)的正溫差越大；壓應(yīng)力越大，應(yīng)力區(qū)與正常區(qū)的負(fù)溫差越大。

3.2 應(yīng)力層厚度變化與檢測(cè)結(jié)果的關(guān)系

圖8 應(yīng)力層厚度變化對(duì)導(dǎo)體表面溫度分布的影響

類(lèi)比裂紋缺陷深度的影響，研究導(dǎo)體應(yīng)力層的厚度改變時(shí)，其表面溫度響應(yīng)曲線的變化規(guī)律,正常區(qū)域、拉應(yīng)力層以及壓應(yīng)力層的電導(dǎo)率分別設(shè)置為5.0×106,2.0×106,8.0×106S·m-1。則由式(5)可計(jì)算得到冷卻階段結(jié)束時(shí)模擬拉、壓應(yīng)力層的熱擴(kuò)散深度分別為2.68,5.35 mm。因此，模擬應(yīng)力層的厚度不能超過(guò)2.68 mm。在仿真分析時(shí)，將最大應(yīng)力層厚度設(shè)置為2.5 mm，應(yīng)力層厚度從0.5 mm開(kāi)始,以0.1 mm為增量逐漸疊加。應(yīng)力層厚度變化對(duì)導(dǎo)體表面溫度分布的影響如圖8所示。

由圖8可見(jiàn)，拉應(yīng)力層厚度越大時(shí)，靠近應(yīng)力區(qū)的感應(yīng)渦流越多、越集中，因此應(yīng)力層表面溫度分布越集中，最大溫度值越大；相反，壓應(yīng)力層厚度增大時(shí)，感應(yīng)渦流逐漸向正常區(qū)域靠近，應(yīng)力層表面溫度分布越分散，最大溫度值越小。對(duì)應(yīng)應(yīng)力層厚度變化時(shí)，在加熱結(jié)束及冷卻結(jié)束時(shí)點(diǎn)1處的溫度值均相應(yīng)發(fā)生變化，如圖9所示。因此可以通過(guò)導(dǎo)體應(yīng)力層表面溫度的變化來(lái)反演應(yīng)力層的厚度值。

圖9 不同應(yīng)力層厚度對(duì)應(yīng)點(diǎn)1處的溫度響應(yīng)曲線

圖10 應(yīng)力層厚度對(duì)點(diǎn)1處溫度值的影響

由圖8,9可以看出，拉應(yīng)力層在冷卻階段結(jié)束時(shí)的溫度隨厚度的變化明顯，而壓應(yīng)力層在加熱結(jié)束時(shí)的溫度隨厚度的變化明顯，在加熱結(jié)束時(shí)模擬壓應(yīng)力層的熱擴(kuò)散深度已達(dá)到3.1 mm，大于最大厚度2.5 mm，因此后續(xù)探究應(yīng)力層厚度與溫度變化之間的關(guān)系時(shí)，拉應(yīng)力選取冷卻階段結(jié)束時(shí)的數(shù)據(jù)，壓應(yīng)力層選取加熱結(jié)束時(shí)的數(shù)據(jù)。應(yīng)力層厚度對(duì)點(diǎn)1處溫度值的影響如圖10所示，導(dǎo)體應(yīng)力區(qū)域應(yīng)力層厚度與表面溫度近似為線性函數(shù)關(guān)系。拉應(yīng)力層厚度越大，其表面溫度就越高；壓應(yīng)力層厚度越大，其表面溫度就越低。

4 結(jié)論

討論了渦流熱成像法檢測(cè)應(yīng)力的原理，建立了數(shù)值仿真模型，利用多物理場(chǎng)耦合的有限元方法，分析了導(dǎo)體應(yīng)力層厚度及電導(dǎo)率、熱導(dǎo)率大小變化對(duì)表面不同區(qū)域點(diǎn)的溫度響應(yīng)曲線的影響，根據(jù)仿真結(jié)果，給出了應(yīng)力大小及應(yīng)力層厚度變化與溫度之間的關(guān)系。

(1) 在一定范圍內(nèi)，導(dǎo)體應(yīng)力值與其表面檢測(cè)溫度近似為線性關(guān)系，即可以通過(guò)導(dǎo)體表面各點(diǎn)的溫度變化情況反演得到導(dǎo)體表面應(yīng)力的大小。若導(dǎo)體表面存在拉應(yīng)力，則應(yīng)力區(qū)的溫度將高于正常區(qū)的，且拉應(yīng)力越大，應(yīng)力層的溫度越高；若導(dǎo)體表面存在壓應(yīng)力，則應(yīng)力區(qū)的溫度將低于正常區(qū)的，且壓應(yīng)力越大，應(yīng)力層的溫度越低。

(2) 在已知應(yīng)力值大小的條件下，導(dǎo)體應(yīng)力層厚度與其表面檢測(cè)溫度近似為線性函數(shù)關(guān)系，拉應(yīng)力層越厚，對(duì)應(yīng)表面溫度越高；壓應(yīng)力層越厚，對(duì)應(yīng)表面溫度越低。即可以通過(guò)溫度值的大小來(lái)分別確定拉、壓應(yīng)力層厚度。