激光輻照光學(xué)窗口的溫度場(chǎng)分析

易 瑔,戴顏斌,楊建昌,張宇航

(陸軍裝甲兵學(xué)院兵器與控制系,北京 100072)

1 引 言

在軍事上,如何干擾和反制敵方的觀察設(shè)備或制導(dǎo)導(dǎo)彈,一直是地面部隊(duì)關(guān)注的難題。隨著激光技術(shù)的進(jìn)步,小型化、高功率的激光武器系統(tǒng)逐漸成為可能,世界各國都在加緊研制用于反導(dǎo)、干擾的激光武器[1]。輻照激光在對(duì)抗光電偵查系統(tǒng)時(shí),往往首先作用到系統(tǒng)的窗口材料上,隨著能量的注入使目標(biāo)溫度上升,最后達(dá)到高溫熔融或者應(yīng)力破裂的效果。因此,研究激光輻照窗口材料的溫度場(chǎng),有利于厘清激光破壞窗口材料的機(jī)理,進(jìn)一步為激光武器的指標(biāo)設(shè)計(jì)、性能優(yōu)化提供科學(xué)理論。

對(duì)于玻璃材料的溫度場(chǎng)分析,國內(nèi)外有許多研究成果[2-7]。但有些研究將模型等效為無限大或半無限大平板,采用一維的熱傳導(dǎo)模型研究,忽略熱量在平面方向上的擴(kuò)散；有些采用二維熱傳導(dǎo)模型研究靶材表面的溫度分布,雖然溫度分布規(guī)律在熱平衡之前是相同的,但是因?yàn)閷⑦吔鐥l件設(shè)置為絕熱,與實(shí)際情況差異較大。因此,為彌補(bǔ)上述不足,并為軍用激光武器的設(shè)計(jì)選型提供理論依據(jù),本文基于傅里葉定律的三維的熱傳導(dǎo)模型,設(shè)置換熱的邊界條件,計(jì)算得到其解析解,并對(duì)不同窗口材料的輻照情況進(jìn)行計(jì)算和仿真,分析溫度場(chǎng)變化的規(guī)律。

2 三維溫度場(chǎng)理論模型

本文研究的對(duì)象是形狀為平行平板玻璃的光學(xué)窗口,在溫度場(chǎng)分析中假設(shè)長(zhǎng)寬都為A,厚度為B,如圖1所示。

圖1 物理模型示意圖

考慮到窗口玻璃的厚度一般遠(yuǎn)小于長(zhǎng)、寬,玻璃正面和背面的熱交換遠(yuǎn)大于四條窄邊,因此在模型中假設(shè)玻璃的四條窄邊絕熱,兩個(gè)正方形平面與空氣換熱。在XOY平面接受激光的輻照,沒有其他熱源。因此,熱傳導(dǎo)方程如下:

(1)

(2)

式中,a是光吸收系數(shù);R是玻璃反射系數(shù)；P為激光總功率；激光束腰半徑ω0。假設(shè)光束幾乎不擴(kuò)散,在材料表面的光斑半徑也為ω0。

邊界條件為:

(3)

因?yàn)槟繕?biāo)關(guān)于z軸中心對(duì)稱,因此只取空間中第一象限研究,式(3)中前兩式表示目標(biāo)在窄邊絕熱,并且在XOZ,YOZ面因?yàn)閷?duì)稱不會(huì)發(fā)生熱交換；后兩式表示目標(biāo)在激光入射面和玻璃背面與空氣發(fā)生對(duì)流換熱,換熱系數(shù)為H,假設(shè)兩處換熱系數(shù)相等。設(shè)初始溫度T0為與時(shí)間無關(guān)的常數(shù)。

定解方程是關(guān)于空間和時(shí)間的二階偏微分方程,利用積分變換法[8]可以將關(guān)于空間的二階偏微分項(xiàng)去掉,變成容易求解的一階偏微分方程。建立正、逆變換:

(4)

(5)

Z(ηp,z)=ηpcos(ηpz)+Hsin(ηpz),N(ηp)

(6)

(7)

(8)

式中：

(9)

(10)

3 計(jì)算結(jié)果與規(guī)律分析

在計(jì)算和仿真中,對(duì)三種常用的光學(xué)窗口材料進(jìn)行計(jì)算,其熱物性參數(shù)如表1所示。光源選取波長(zhǎng)1064 nm、總功率100 W的連續(xù)激光,光斑半徑為1 mm,能量呈高斯分布(標(biāo)準(zhǔn)差σ=0.5),室溫20 ℃。窗口材料設(shè)定為長(zhǎng)50 mm、高50 mm、寬5 mm的平行平板。

表1 幾種窗口材料的熱物性參數(shù)

以K9玻璃為例,經(jīng)過500 s激光輻照后,表面溫度場(chǎng)趨于穩(wěn)定,其分布如圖2(a)所示。人由于光源是高斯光束,熱平衡后的溫度場(chǎng)也近似于高斯分布,如圖2(b)所示。

圖2 K9玻璃受幅照表面溫度場(chǎng)分布及等溫線

3.1 激光參數(shù)對(duì)溫度場(chǎng)影響

該模型中,將激光的作用等效為一個(gè)能量密度隨位置變化的體積熱源,假設(shè)使用的是高斯激光,并且光在通過玻璃時(shí)逐漸被吸收,則體積熱源的能量分布是沿Z軸逐漸衰減、在垂直于Z軸的任意截面呈正態(tài)分布。因此,研究激光參數(shù)變化對(duì)溫度場(chǎng)的影響,只需要考慮激光參數(shù)對(duì)等效熱源的影響即可。

根據(jù)解析式(8),某點(diǎn)溫度正比于該點(diǎn)熱源功率密度對(duì)于g時(shí)間t的積分,這意味著越大的激光功率P和越小的光斑半徑r,能產(chǎn)生更高的溫升,符合我們對(duì)激光輻射的常規(guī)認(rèn)識(shí),如圖3所示。

圖3 不同參數(shù)下等效熱源功率密度表面分布

對(duì)于激光光斑抖動(dòng)、束散或漂移的情況,可以視作增大輻照光斑、降低能量密度進(jìn)行處理。例如,采用蒙文等提出的光斑抖動(dòng)模型[9],將半徑2 mm、總功率100 W的激光光斑等效為抖動(dòng)光斑,其半徑和功率變?yōu)?

R′=0.6827Ltanθ+R

(11)

(12)

式中,R,R′指原光斑和等效光斑的半徑；P,P′指原光束總功率和等效光束功率；L為光源到材料表面距離；θ為誤差偏離角。假設(shè)距離為1 m,偏離角0.2 mrad,對(duì)等效熱源密度的影響如圖4所示。

圖4 引入光斑抖動(dòng)前后的等效熱源功率密度對(duì)比

3.2 材料參數(shù)對(duì)溫度場(chǎng)影響

在激光參數(shù)不變的情況下,材料性質(zhì)的差異顯然是溫度場(chǎng)差異的決定因素。由公式中分析可知,導(dǎo)熱系數(shù)k,密度ρ,比熱C,反射率R和體積尺寸等參數(shù)決定了某時(shí)刻的溫度場(chǎng)數(shù)值。對(duì)于常見的三種光學(xué)窗口材料,整理了其熱學(xué)性質(zhì),并在同種光源和不同光源情況下對(duì)其溫度場(chǎng)進(jìn)行了分析。

下面是幾種常見的可見光和紅外波段窗口材料激光輻照仿真的結(jié)果。可見長(zhǎng)時(shí)間輻照后,材料整體溫度分布與等效熱源功率密度分布輪廓相似,且因?yàn)檫_(dá)到了熱平衡,此時(shí)三種材料都沒有達(dá)到熔融或軟化點(diǎn)如圖5所示。

圖5 三種材料溫度場(chǎng)仿真結(jié)果

(13)

當(dāng)材料尺寸相同、輻照時(shí)間相同時(shí),公式(13)中的特征值和積分時(shí)間是相同的；c是與時(shí)間無關(guān)的常數(shù)項(xiàng)。因此,可以看出目標(biāo)材料的溫升大致上與密度ρ、比熱C和反射率R成負(fù)相關(guān)；與吸收系數(shù)a成正相關(guān)。

考慮到激光毀傷光學(xué)窗口材料的需求,理解吸收系數(shù)a對(duì)于溫度場(chǎng)的影響十分必要,但是通過固定光源參數(shù)、改變材料的方法,不能考察吸收系數(shù)對(duì)溫度場(chǎng)的獨(dú)立影響,因?yàn)椴牧系钠渌麉?shù)例如密度、比熱也改變了。吸收系數(shù)是相關(guān)于材料本身和光波長(zhǎng)的屬性,因此這里考察同種材料對(duì)應(yīng)不同波長(zhǎng)激光的溫度場(chǎng)變化情況,就能分析吸收系數(shù)對(duì)溫度場(chǎng)的影響。在反射率R接近于0,厚度為d,透射率為T時(shí),可以認(rèn)為:

(14)

因此,對(duì)于不同波段激光輻照的同種材料,其透過率和反射率越小,吸收系數(shù)越大,對(duì)應(yīng)單位時(shí)間內(nèi)輻照溫升就越大。以K9玻璃在1.064 μm、1.3 μm、1.5 μm這三種光源下的溫升情況為例,如圖6所示。

圖6 不同吸收系數(shù)對(duì)K9玻璃輻照中心點(diǎn)溫升的影響

由圖可見,隨著K9玻璃的光吸收系數(shù)增大,溫度場(chǎng)上升速度大幅增加,并且換熱情況下的熱平衡溫度也上升了。

3.3 邊界換熱系數(shù)對(duì)溫度場(chǎng)影響

一些文獻(xiàn)中,計(jì)算時(shí)將邊界條件設(shè)置為絕熱,雖然能得到形式簡(jiǎn)單的解,但是會(huì)造成一下錯(cuò)誤:第一,材料溫度隨著輻照時(shí)間延長(zhǎng),可以無限上升；第二,沿光軸方向的溫度基本相同。這是因?yàn)榻^熱條件下的模型不考慮與空氣對(duì)流換熱,只有熱量輸入沒有輸出,所以不會(huì)達(dá)到熱平衡點(diǎn)；同時(shí)不考慮散熱也影響了沿光軸方向的溫升計(jì)算,使得溫度在此方向變化很小。使用這樣的模型會(huì)導(dǎo)致忽略溫度上升限制,與實(shí)際情況不符的后果,并且計(jì)算熱應(yīng)力時(shí)也會(huì)受到影響。圖7和圖8是將邊界條件設(shè)置為與空氣換熱,計(jì)算輻照1064 nm激光輻照500 s時(shí),與同條件的絕熱情況進(jìn)行對(duì)比。

圖7 換熱與絕熱條件下溫升情況對(duì)比

圖8 K9玻璃主光軸溫度分布

由圖可見,換熱條件下溫度升高到一定數(shù)值后達(dá)到穩(wěn)定,不隨輻照時(shí)間增加而上升,但是絕熱條件下溫度可以一直上升；換熱條件下沿光軸方向的溫度分布是中部高、兩側(cè)低,并且峰值靠近輻照面,這是由于兩個(gè)邊界正進(jìn)行對(duì)流換熱,并且等效熱源在前表面能量密度大造成的。絕熱條件下的光軸方向,溫度相差不超過0.04 ℃,差異不大。

換熱模型中,考慮模型達(dá)到熱平衡時(shí),材料溫度不再上升,激光的能量輸入與空氣換熱能量輸出持平:

(15)

式中,S指換熱表面面積,該模型中考慮入射面和出射面與空氣換熱,左式表示熱平衡時(shí)的流出能量,右式表示流入能量。顯然當(dāng)激光參數(shù)不變時(shí),右式為常數(shù),因此平衡溫度T與換熱系數(shù)H成反比關(guān)系,即換熱系數(shù)越大,熱平衡時(shí)的溫度越低,如圖9所示。

圖9 不同換熱系數(shù)下K9玻璃輻照中心溫升情況

4 結(jié) 論

利用積分變換的方法,求解了包含換熱條件的三維溫度場(chǎng)模型,并對(duì)三種常用光學(xué)窗口材料進(jìn)行計(jì)算和仿真。將激光光源參數(shù)轉(zhuǎn)換為等效熱源的屬性,通過改變等效熱源的參數(shù)可以計(jì)算溫度場(chǎng)隨波長(zhǎng)、功率、光斑半徑以及光斑抖動(dòng)等影響的變化情況,分析了大激光功率和小光斑半徑有利于產(chǎn)生更高的中心溫度；對(duì)于不同材料,其密度、比熱和反射率越小越利于溫升,對(duì)應(yīng)的光吸收系數(shù)越高,溫升越快,熱平衡溫度也更高；對(duì)比了常見的絕熱模型和本文的換熱模型在熱平衡情況和主光軸溫度分布上的差異,材料實(shí)際上受輻照最終達(dá)到熱平衡,且平衡溫度與換熱系數(shù)成反比、與光源能量成正比,本文建立的換熱模型更符合實(shí)際。