汪慶桃，錢秋冬，李 霄，王崇旭

(1. 國(guó)防科技大學(xué) 軍事基礎(chǔ)教育學(xué)院， 湖南 長(zhǎng)沙 410072； 2. 國(guó)防科技大學(xué) 空天科學(xué)學(xué)院， 湖南 長(zhǎng)沙 410073； 3. 國(guó)防科技大學(xué) 前沿交叉學(xué)科學(xué)院， 湖南 長(zhǎng)沙 410073)

激光輻照下金屬板的動(dòng)態(tài)響應(yīng)對(duì)于激光加工、未爆彈藥銷毀等具有重要的意義，Chun等[1]指出，這一過(guò)程中同時(shí)包含著多種復(fù)雜的物理過(guò)程。在激光輻照下，大量的激光能量被材料表面吸收并轉(zhuǎn)化為熱能，熱量向內(nèi)部傳遞使得材料迅速升溫，形成相應(yīng)的溫度場(chǎng)。若材料溫度達(dá)到其熔點(diǎn)或沸點(diǎn)，則會(huì)發(fā)生熔融甚至氣化現(xiàn)象，由此產(chǎn)生的粒子繼續(xù)被加熱形成蒸汽等離子體，使得材料發(fā)生一系列變化，產(chǎn)生損傷效應(yīng)[2]。相關(guān)實(shí)驗(yàn)結(jié)果也表明，激光與材料之間相互作用最顯而易見(jiàn)的效果是材料熔融等燒蝕破壞現(xiàn)象[3]。

Sami[4]、Anisimov[5]、Valette[6]等采用數(shù)值計(jì)算方法研究了激光輻照下金屬材料的溫度變化，并用物理方法估算了材料表面徑向熱擴(kuò)散區(qū)的范圍。Le Harzic等[7]利用透射電子顯微鏡觀察了金屬鋁在激光輻照作用下的熱影響區(qū)，發(fā)現(xiàn)表面燒蝕損傷明顯。王偉平[8-9]、Chimier[10]等建立了激光輻照下的熱傳導(dǎo)數(shù)值模型，由此得到了金屬靶的溫度場(chǎng)分布。在此基礎(chǔ)上，梁業(yè)廣[11-12]、趙鳳艷[13]等考慮金屬材料熔融相變的影響，分析了模型中材料參數(shù)和激光參數(shù)對(duì)計(jì)算結(jié)果的影響。通過(guò)研究，王譯那等[14]進(jìn)一步發(fā)現(xiàn)激光輻照的整個(gè)破壞過(guò)程主要包括升溫、穿孔和孔徑增大三個(gè)階段。趙方東[15]、何雅靜[16]等對(duì)激光輻照金屬靶過(guò)程進(jìn)行了數(shù)值模擬，得到了金屬靶熔化深度、速率以及不同深度處溫度等隨激光作用時(shí)間的變化規(guī)律。

本文以45#鋼為研究對(duì)象，采取實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的方法，對(duì)激光輻照下靶板的熱響應(yīng)特性展開(kāi)研究，為激光加工金屬材料或者破壞金屬結(jié)構(gòu)提供理論參考。

1 激光輻照下45#鋼靶的熱響應(yīng)實(shí)驗(yàn)

1.1 實(shí)驗(yàn)條件

實(shí)驗(yàn)采用連續(xù)光纖激光器系統(tǒng)，該系統(tǒng)由電源、激光器以及冷卻裝置三個(gè)部分組成，波長(zhǎng)為1 080 nm。測(cè)溫裝置采用接觸式與非接觸式兩種方式測(cè)溫。接觸式測(cè)溫采用24通道MT-X溫度記錄儀、SH-K型鎳鉻-銅鎳熱電偶實(shí)現(xiàn)數(shù)據(jù)采集，測(cè)溫范圍為-100～1 768 ℃，精度為±0.5 ℃，分辨率可達(dá)0.1 ℃。非接觸式測(cè)溫采用FLUKE手持式Tis75紅外熱像儀，其測(cè)溫范圍為-20～1 500 ℃。鋼靶為直徑100 mm圓柱形，厚度分別為3 mm、4 mm、5 mm、8 mm四種情況，實(shí)驗(yàn)裝置示意如圖1所示。實(shí)驗(yàn)中激光發(fā)射功率分別為546 W、699 W兩種功率，光路長(zhǎng)度為4.8 m，光斑直徑為10 mm。

圖1 實(shí)驗(yàn)裝置示意Fig.1 Schematic diagram of the experimental device

實(shí)驗(yàn)前借助激光觀測(cè)儀調(diào)節(jié)激光器光路與金屬靶板相對(duì)位置，使激光入射方向垂直于靶板正面且激光輻照光斑中心點(diǎn)位于金屬靶板圓心處，選取靶板受輻照面圓心為正面中心點(diǎn)。在靶板輻照的背面以圓心為起點(diǎn)沿徑向設(shè)置4個(gè)SH-K型熱電偶，接觸應(yīng)牢固可靠，熱電偶的另一端與MT-X多路溫度記錄儀相連接，完成數(shù)據(jù)的記錄與儲(chǔ)存，每?jī)蓚€(gè)熱電偶之間的間隔為10 mm，如圖2所示。

圖2 靶板背面熱電偶測(cè)溫點(diǎn)Fig.2 Thermocouple temperature measurement point on the back of the target

鋼靶受到激光輻照作用后吸收熱能，溫度上升。在靶板正面，入射激光光斑中心處溫度最高，研究中記作靶板峰值溫度，峰值溫度超出金屬熔點(diǎn)范圍后金屬靶發(fā)生相變，呈現(xiàn)為液態(tài)。由于激光能量在光斑中心最高并沿半徑衰減，輻照靶材后會(huì)在表面形成近似圓形的液相區(qū)域，將該圓的直徑記作熔池半徑，而液相區(qū)域縱向發(fā)展在中心處最大，該值記作熔池深度。沿徑向溫度逐漸降低，從而導(dǎo)致金屬表面發(fā)生氧化反應(yīng)的程度不同，形成了明暗層次分明的光暈痕跡，光暈近似呈圓形，且尺寸大于熔池。本文以該光暈痕跡直徑來(lái)表示溫度傳導(dǎo)的范圍，實(shí)驗(yàn)結(jié)束后，利用直尺對(duì)每一塊靶板表面的溫度傳導(dǎo)范圍進(jìn)行測(cè)量記錄，如圖3所示。

圖3 溫度傳導(dǎo)范圍的測(cè)定Fig.3 Determination of the temperature conduction range

1.2 實(shí)驗(yàn)結(jié)果

表1為激光輻照作用下的實(shí)驗(yàn)結(jié)果。其中，T3次實(shí)驗(yàn)在輻照過(guò)程中未觀察到燒蝕孔貫穿靶板，輻照結(jié)束后對(duì)熔池深度進(jìn)行測(cè)量時(shí)發(fā)現(xiàn)整個(gè)熔池區(qū)域金屬已全部熔融，液態(tài)金屬雖然尚未完全滴落，但探針插入熔池后可以完全貫穿，故可認(rèn)為靶板已經(jīng)實(shí)現(xiàn)燒蝕穿孔。

由表1可知：熔池尺寸與激光功率直接相關(guān)，當(dāng)激光功率為699 W時(shí)，對(duì)于厚度小于4 mm的靶板均可以成功穿透；同時(shí)，當(dāng)靶板厚度一定時(shí)，隨著激光功率的增大，熔池的直徑和深度均增大；當(dāng)激光功率相同時(shí)，隨著靶板厚度的增加，熔池的直徑與深度均相應(yīng)減小。
圖4給出了幾個(gè)典型工況的溫度場(chǎng)響應(yīng)曲線。

表1 激光輻照實(shí)驗(yàn)結(jié)果

(a) T5

(b) T6圖4 溫度場(chǎng)隨時(shí)間變化曲線Fig.4 Variation curves of temperature with time

由圖4可知，在激光輻照作用下，金屬靶板正背面溫度均有所升高但二者溫升速度和幅值存在較大差異，究其原因，靶板正面能量來(lái)源為入射激光能量，其轉(zhuǎn)化為熱量并直接被吸收，而背面能量來(lái)源為正面所吸收能量的熱傳導(dǎo)；同時(shí)可以發(fā)現(xiàn)，溫度傳導(dǎo)范圍與溫度峰值存在一定的正向關(guān)系，當(dāng)溫度峰值升高時(shí)溫度傳導(dǎo)范圍也隨之增大，溫度場(chǎng)峰值溫度和溫度傳導(dǎo)范圍隨著靶板厚度的增大而減小。

2 激光輻照下45#鋼板熱響應(yīng)數(shù)值模擬

2.1 數(shù)值模型及驗(yàn)證

基于有限元軟件ANSYS建立三維數(shù)值模型，靶板離散成八節(jié)點(diǎn)實(shí)體單元Solid70，采用映射法劃分網(wǎng)格。為了統(tǒng)籌計(jì)算規(guī)模和計(jì)算效率之間的關(guān)系，采用變網(wǎng)格技術(shù)，即在激光輻照及影響區(qū)域采用較為精細(xì)的網(wǎng)格?？紤]到模型的對(duì)稱性，僅建立四分之一模型，這樣在總體網(wǎng)格數(shù)量保持不變的情況下可以盡可能地使網(wǎng)格尺寸更加精細(xì)。靶板為圓柱形結(jié)構(gòu)，直徑為100 mm，材料參數(shù)參考文獻(xiàn)[17]。當(dāng)材料的溫度超過(guò)熔點(diǎn)范圍最高溫度時(shí)，判定該材料被侵蝕，程序自動(dòng)刪除該網(wǎng)格單元的質(zhì)量，用以模擬燒蝕穿孔現(xiàn)象。

圖5為激光功率密度為2 000 W/cm2、光斑直徑為10 mm的連續(xù)激光輻照10 mm厚度金屬板時(shí)熱響應(yīng)云圖。

從圖5(a)可以看出，在激光輻照初期，靶板迅速升溫，溫度傳導(dǎo)到一定的范圍。靠近光斑中心區(qū)域溫度較高，沿著徑向方向往外溫度值逐漸降低。此時(shí)，金屬靶板吸收能量?jī)H用于溫度增長(zhǎng)，尚未形成熔池。當(dāng)t=4.55 s時(shí)，溫度升高至超出材料熔點(diǎn)范圍上限，單元被刪除，形成孔洞，如圖5(b)所示。隨著激光輻照過(guò)程的持續(xù)，溫度傳導(dǎo)的范圍不斷擴(kuò)大，溫度持續(xù)上升，孔洞不斷向縱深發(fā)展，如圖5(c)～(d)所示。當(dāng)t=10.55 s時(shí)，靶板背面中心區(qū)域單元溫度超出熔點(diǎn)上限，此時(shí)實(shí)現(xiàn)靶板的貫穿，如圖5(e)所示。在此之后，隨著溫度的進(jìn)一步升高，穿孔沿徑向發(fā)展，孔形由倒圓錐狀逐漸向倒圓臺(tái)狀過(guò)渡，并呈現(xiàn)出向圓柱狀轉(zhuǎn)變的趨勢(shì)，如圖5(f)所示。

(a) t=3.0 s (b) t=4.55 s (c) t=7.05 s

(d) t=9.05 s (e) t=10.55 s (f) t=30.0 s圖5 金屬靶板燒蝕穿孔過(guò)程Fig.5 Process of ablation and perforation of metal target

圖6為工況T6(激光功率699 W、光斑直徑10 mm、金屬板厚度8 mm)的數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果的溫度場(chǎng)分布比較。圖6(a)～(b)中靶板表面中心區(qū)域溫度值最高，向四周傳播的同時(shí)溫度值下降，由于徑向溫度梯度的存在，表面呈現(xiàn)規(guī)則分布的年輪狀溫度場(chǎng)，層次分明，與實(shí)驗(yàn)結(jié)果中靶板表面產(chǎn)生的光暈現(xiàn)象吻合較好。圖6(c)～(d)中各特征點(diǎn)的數(shù)值模擬溫升曲線與實(shí)驗(yàn)結(jié)果中所測(cè)得的曲線趨勢(shì)基本一致。

(a) 數(shù)值模擬溫度場(chǎng)分布(a) Temperature field distribution of numerical simulation

(b) 實(shí)驗(yàn)溫度場(chǎng)分布(b) Temperature field distribution of experimental results

(c) 正面中心點(diǎn)溫升曲線(c) Temperature rise curve of front center point

(d) 背面特征點(diǎn)溫升曲線(d) Temperature rise curve of characteristic point on the back圖6 數(shù)值模擬與實(shí)驗(yàn)結(jié)果對(duì)比Fig.6 Comparison of numerical simulation and experimental results

表2將幾個(gè)典型工況下靶板燒蝕穿孔特性的數(shù)值模擬與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，可以發(fā)現(xiàn)，對(duì)于熔池直徑與深度尺寸的數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合也較好，由此證明了數(shù)值模擬結(jié)果的準(zhǔn)確性。

2.2 影響因素分析

采用數(shù)值模擬方法，研究靶板厚度、激光功率以及光斑半徑等參數(shù)對(duì)激光輻照下鋼靶板熱響應(yīng)特性的影響規(guī)律，數(shù)值模擬工況及模擬結(jié)果如表3所示。

2.2.1 靶板厚度的影響

靶板厚度對(duì)熱響應(yīng)特性影響很大，在輻照功率相同的情況下，若靶板過(guò)厚，則熱耗散增大，傳導(dǎo)到靶板背面使其升溫變緩，同時(shí)燒蝕穿孔也相對(duì)較為困難。數(shù)值計(jì)算結(jié)果表明，當(dāng)靶板厚度較小(1 mm和2 mm)時(shí)，激光輻照下靶板溫度迅速升高，靶板被灼穿。圖7給出了在激光功率為1 000 W、光斑直徑為10 mm時(shí)靶板峰值溫度、溫度傳導(dǎo)范圍、熔池尺寸等參量隨靶板厚度的變化曲線。

表2 數(shù)值模擬與實(shí)驗(yàn)結(jié)果對(duì)比

表3 數(shù)值模擬結(jié)果

由圖7可知，靶板厚度對(duì)金屬靶板熱響應(yīng)特性影響顯著，靶板峰值溫度、溫度傳導(dǎo)范圍以及熔池尺寸等均隨靶板厚度的增大而減小，且近似呈指數(shù)變化關(guān)系。

(a) 金屬靶板峰值溫度(a) Peak temperature of metal target plate

(b) 溫度傳導(dǎo)范圍(b) Temperature conduction range

(c) 熔池直徑(c) Melting pool diameter

(d) 熔池深度(d) Melting pool depth圖7 金屬靶板熱響應(yīng)特性與靶板厚度關(guān)系Fig.7 Relationship between the thermal response characteristics of the metal target and the target thickness

2.2.2 激光功率的影響

激光功率是金屬靶板激光輻照過(guò)程中的關(guān)鍵參數(shù)，直接決定了整個(gè)輻照過(guò)程的能量輸入大小。圖8給出了當(dāng)光斑直徑為10 mm時(shí)，不同激光功率情況下10 mm厚靶板的峰值溫度、溫度傳導(dǎo)范圍和熔池尺寸等參量的變化曲線。

由圖8可知，靶板峰值溫度、溫度傳導(dǎo)范圍以及熔池尺寸均與激光功率近似呈指數(shù)變化關(guān)系，各參數(shù)均隨激光功率的增大而增大。其中，當(dāng)激光功率為500 W和750 W時(shí)，金屬靶板表面溫升有限，未形成熔池。分析原因：當(dāng)靶板由8 mm進(jìn)一步增加到10 mm時(shí)，由于熱傳導(dǎo)所導(dǎo)致的能量耗散進(jìn)一步增大。

(a) 金屬靶板峰值溫度(a) Peak temperature of metal target plate

(b) 溫度傳導(dǎo)范圍(b) Temperature conduction range

(c) 熔池直徑(c) Melting pool diameter

(d) 熔池深度(d) Melting pool depth圖8 金屬靶板熱響應(yīng)特性與激光功率關(guān)系Fig.8 Relationship between the thermal response characteristics of the metal target and the laser power

對(duì)表3中數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比可以發(fā)現(xiàn)，燒蝕穿孔熔池直徑的數(shù)值模擬值普遍偏小，原因主要有：①數(shù)值模擬中以單元溫度為失效準(zhǔn)則，設(shè)置單元溫度超出熔點(diǎn)范圍后對(duì)其進(jìn)行“殺死”操作，但ANSYS中單元溫度為八節(jié)點(diǎn)溫度取平均值，與單元位置處真實(shí)溫度存在差異。由于單元較遠(yuǎn)處節(jié)點(diǎn)溫度值較低，當(dāng)單元尺寸較大時(shí)數(shù)值模擬結(jié)果與實(shí)際溫度相比偏低，增大了誤差。②靶板表面有雜質(zhì)，并且氧化后吸熱較強(qiáng)，數(shù)值模擬中未考慮由此帶來(lái)的吸收率增長(zhǎng)，導(dǎo)致表面溫升較低，熔池生成受限。

2.2.3 光斑直徑的影響

光斑直徑是激光輻照過(guò)程中十分重要的參數(shù)，決定了入射激光的能量分布范圍。圖9給出了當(dāng)激光功率為2 000 W時(shí)，不同光斑直徑情況下10 mm厚鋼靶的峰值溫度、溫度傳導(dǎo)范圍和熔池尺寸。由圖9可知，峰值溫度、溫度傳導(dǎo)范圍以及熔池深度與光斑直徑近似呈指數(shù)變化關(guān)系，熔池直徑與光斑直徑近似呈二次多項(xiàng)式變化關(guān)系，鋼靶峰值溫度以及熔池尺寸隨光斑直徑的增大而減小，溫度傳導(dǎo)范圍隨光斑直徑的增大而增大。其中，當(dāng)光斑直徑較小時(shí)(10 mm和15 mm)，靶板表面溫升較高，但能量相對(duì)集中限制了熔池向外擴(kuò)展，導(dǎo)致熔池直徑較小。

(a) 金屬靶板峰值溫度(a) Peak temperature of metal target plate

(b) 溫度傳導(dǎo)范圍(b) Temperature conduction range

(c) 熔池直徑(c) Melting pool diameter

(d) 熔池深度(d) Melting pool depth圖9 金屬靶板熱響應(yīng)特性與光斑直徑關(guān)系Fig.9 Relationship between the thermal response characteristics of the metal target and the spot diameter

2.3 特征參數(shù)隨初始條件的關(guān)系

由數(shù)值模擬結(jié)果可以得出靶板峰值溫度、溫度傳導(dǎo)范圍和熔池尺寸等特征參數(shù)隨靶板厚度、激光功率及光斑直徑等初始條件的變化關(guān)系，假定特征參數(shù)F隨初始條件的變化形式為：

F=k0·Ik1·Lk2·dk3

(1)

式中，ki(i=0，1，2，3)為待定參數(shù)。

根據(jù)表3中數(shù)值模擬結(jié)果，借助函數(shù)繪圖軟件Origin對(duì)靶板峰值溫度、溫度傳導(dǎo)范圍、熔池直徑、熔池深度隨靶板厚度、激光功率及光斑直徑的變化關(guān)系進(jìn)行數(shù)據(jù)擬合，可以得到式(1)中相關(guān)參數(shù)，如表4所示。

表4 參數(shù)擬合結(jié)果

3 結(jié)論

采用實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的方法，對(duì)激光輻照下45#鋼板的響應(yīng)特性進(jìn)行了研究，主要結(jié)論如下：

1)在入射激光功率和光斑直徑相同的條件下，靶板峰值溫度、溫度傳導(dǎo)范圍以及熔池尺寸均隨靶板厚度的增大而減小，在已有結(jié)果范圍內(nèi)呈指數(shù)關(guān)系；

2)在靶板厚度和光斑直徑相同的條件下，靶板峰值溫度、溫度傳導(dǎo)范圍以及熔池尺寸均隨著激光功率的增大而增大，在已有結(jié)果范圍內(nèi)呈指數(shù)關(guān)系；

3)在入射激光功率和靶板厚度相同的條件下，隨著光斑直徑的增大，靶板峰值溫度以及熔池尺寸呈現(xiàn)減小的趨勢(shì)，而溫度傳導(dǎo)范圍呈現(xiàn)增大的趨勢(shì)。