蔡繼興，張燁，王頔，渠旭，李賀，金光勇

（長春理工大學 固體激光技術與應用吉林省重點實驗室，長春 130022）

激光對熔石英及其表面膜層材料的損傷是一個復雜的物理過程，它主要由兩方面的性質決定：一是激光參數(shù)，如激光脈寬、能量密度、重復頻率、光斑半徑等；二是材料的自身屬性。不同的激光條件會對同一種材料產生不同的作用效果，相同激光條件下對不同的材料也會產生不同的損傷效果［1］。根據(jù)損傷效果的不同分為熱損傷、應力應變損傷、致燃損傷、爆轟波損傷以及多種損傷效應的耦合作用等。它們都在很短的時間內產生相互影響，給研究工作帶來極大的困難。1064nm增透熔石英是光學系統(tǒng)中常用的窗口材料，光路中任何窗口材料的損傷都會造成很大的損失，甚至造成整個系統(tǒng)不能穩(wěn)定運行。因此，研究1064nm增透熔石英窗口的損傷，特別是對其損傷效果及損傷機理的研究具有重要的意義和應用價值［2］。

1 理論模型

1.1 幾何模型

由于激光在空間中為高斯分布，所以計算時采用圓柱坐標系，假設熔石英及其表面膜層材料為各向同性的連續(xù)介質，激光輻照到材料表面后的能量吸收為體吸收［3］，則所建立的激光損傷1064nm增透熔石英的幾何模型如圖1所示，其中膜層厚度h=5.5×10-4mm，熔石英基底厚度d=4mm，材料半徑a=11.25mm。

圖1 幾何模型

1.2 數(shù)學模型

1.2.1 熱傳導方程

由于采用的是毫秒脈沖激光，其與1064nm增透熔石英的能量交換過程主要體現(xiàn)為熱交換，忽略材料與周圍空間的輻射和對流過程，熱傳導方程可以表示為如下形式［4］：

其中，Ti(r,z,t)為t時刻的溫度分布，ki、ρi和ci分別代表材料的熱傳導系數(shù)、密度和比熱，i=f、s分別代表膜層和熔石英基底，上述熱傳導方程的熱源項可以表示為［5］：

式中：

其中，f(r)和g(t)分別代表激光的空間和時間分布；r、z分別為柱坐標系的徑向、軸向坐標，r0為激光光斑半徑，αi為膜層和熔石英的吸收系數(shù)，τp為激光脈寬，Ri為膜層和熔石英的反射率，I0為激光的中心功率密度［6］。

1.2.2 彈性力學方程

在圓柱坐標系下，與上述熱傳導方程耦合的彈性力學方程可以表示為如下形式［7，8］：

其中，uri和uzi分別代表位移在r和z方向上的分量；μi、εi和βi分別為材料的泊松比、體應變和熱應力系數(shù)。

2 數(shù)值模擬結果和分析

圖2所示為在能量密度3.1831×103J/cm2，脈寬1ms，光斑半徑1mm的單脈沖激光作用下，1064nm增透熔石英上表面溫度隨徑向位置和時間的變化關系。由于激光束為高斯對稱分布，根據(jù)1064nm增透熔石英的熔點特性，可知1064nm增透熔石英在t=0.96ms時熔融區(qū)的半徑r=0.2mm，從而計算出激光作用后1064nm增透熔石英熔融熱損傷面積約為0.126mm2。

圖2 上表面溫度隨徑向位置和時間的變化

對應特征徑向位置和特征時刻作二維投影，給出上述激光參數(shù)條件下1064nm增透熔石英上表面特征徑向位置處，溫度隨時間的變化關系，如圖3所示；特征時刻1064nm增透熔石英上表面溫度隨徑向位置的變化關系如圖4所示。

圖3 上表面特征徑向位置溫度隨時間的變化

從圖3可知，在激光光斑半徑區(qū)域內，隨激光作用時間的增加，1064nm增透熔石英上表面各點溫度值均逐漸增加。在激光輻照中心點產生的溫升最大值約為1360K，超過了膜層材料的熔點值；在光斑半徑邊緣產生的溫升值最小，約為380K，沿材料半徑方向各點的溫升呈下降趨勢。

圖4 上表面特征時刻溫度隨徑向位置的變化

由圖4可知，激光輻照時間越長，溫升越高，材料表面產生的溫度梯度也越大。在激光輻照區(qū)域內，溫升較集中，在激光光斑邊緣附近，形成較大的溫度梯度。當激光作用結束時，形成的溫度梯度達到最大，據(jù)此可知在材料內部將產生較大的熱應力。因此，判斷1064nm增透熔石英是否發(fā)生損傷，主要依據(jù)激光作用結束時，材料內部產生的熱應力是否超過了膜層及熔石英的抗拉強度或抗壓強度。

圖5 不同能量密度條件下上表面投影中心點溫度隨脈沖寬度的變化

圖5所示為能量密度分別取2.5465×103J/cm2、2.8648×103J/cm2和 3.1831×103J/cm2激光條件下，1064nm增透熔石英上表面中心點溫度隨激光脈沖寬度的變化關系。當脈沖寬度一定時激光能量密度越高，材料表面溫度越高，由于材料表面中心點的溫度高低主要由該點的激光輻照能量密度大小決定，因此在激光光斑半徑相同的條件下，激光脈寬越大，激光能量密度值就越小，作用到材料表面后產生的溫度逐漸降低。

圖6 上表面投影中心點溫度隨時間的變化

圖6所示為激光重復頻率分別為2Hz、5Hz和10Hz條件下，1064nm增透熔石英上表面投影中心點溫度隨時間的變化關系。在激光輻照期間，中心點的溫度急劇上升，在脈沖停止作用的脈沖間隔內，由于沒有能量的聚集，材料處于冷卻階段，表現(xiàn)為激光作用中心點溫度的降低，溫升曲線呈鋸齒狀。

圖7 上表面特征徑向位置溫度隨時間的變化

圖7所示為相同激光參數(shù)條件下，1064nm增透熔石英上表面特征徑向位置溫度隨時間的變化關系，當?shù)诙€脈沖作用于材料表面后，材料上表面中心點的溫升值約為1763K，第三個脈沖作用于材料表面后，材料上表面中心點的溫升值已經達到1972K，兩者均超過了膜層及熔石英基底的熔點值。因此可以判定，當脈沖串激光作用于1064nm增透熔石英表面后，將會產生一定的溫度累積效應，其效果使材料的損傷程度加劇。

圖8 特征軸向位置溫度隨時間的變化

圖8所示為相同參數(shù)的脈沖串激光作用下1064nm增透熔石英特征軸向位置溫度隨時間的變化關系，1064nm增透熔石英上表面的溫升最大值要高于后表面的溫升最大值，且隨脈沖個數(shù)的增加，軸向不同位置處的溫升最大值依次升高，即在軸線方向上存在一定的溫度累積效應。

圖9 特征時刻溫度隨軸向位置的變化

圖9所示為前三個脈沖激光作用結束后1064nm增透熔石英溫度隨軸向位置的變化關系，可以看到隨激光脈沖個數(shù)的增加，軸向各點的溫度均上升。1064nm增透熔石英在t=101ms、201ms時，其熱熔融區(qū)的深度已經由z=3.2mm增加至z=4.0mm，熱熔融區(qū)已經由上表面至下表面貫穿整個材料。

圖10 1064nm增透熔石英溫度分布

圖10所示為脈沖串激光作用結束后1064nm增透熔石英溫度的三維分布圖，當最后一個脈沖作用結束后，其熱熔融區(qū)域半徑約為0.7mm（單脈沖作用時的熱熔融半徑約為0.2mm）?？梢姡啾扔趩蚊}沖激光，脈沖串激光將會產生一定的溫度累積效應，其效果是使得材料的損傷面積和損傷深度增加。由于激光脈沖為高斯分布，根據(jù)1064nm增透熔石英的熔點、氣化點等特性，可以得到激光脈寬為1ms，能量密度為3.1831×103J/cm2，重復頻率為10Hz條件下，3個脈沖作用截止時，其熱熔融損傷體積約為6.158mm3。

圖11 上表面徑向/環(huán)向應力隨徑向位置的變化

圖11所示為特征時刻1064nm增透熔石英上表面徑向應力、環(huán)向應力隨徑向位置的變化關系，在上表面中心處，徑向應力表現(xiàn)為拉應力，且拉應力達到最大值，隨脈沖個數(shù)的增加，拉應力逐漸增加；在半徑方向上，拉應力逐漸減小，并逐漸轉化為壓應力。同時發(fā)現(xiàn)，在激光光斑半徑邊緣附近壓應力達到了最大值。據(jù)此可以判斷，1064nm增透熔石英發(fā)生損傷時，應首先從激光作用中心點或激光光斑半徑邊緣附近開始。從圖11（b）可以看到，1064nm增透熔石英上表面環(huán)向應力始終體現(xiàn)為拉應力，且隨脈沖個數(shù)的增加，拉應力逐漸增大。在半徑方向上，拉應力先逐漸增加，并在光斑半徑邊緣附近達到最大，然后逐漸減小并最終趨向于零。根據(jù)1064nm增透熔石英的應力特性，可以得到1064nm增透熔石英在t=201ms時超過應力屈服強度區(qū)域對應的徑向半徑r=1.5mm，從而計算出激光作用后1064nm增透熔石英應力損傷面積約為7.069mm2。

圖12所示為1064nm增透熔石英徑向應力、軸向應力隨軸向位置的變化關系，在材料表面中心處，徑向應力表現(xiàn)為拉應力，且隨脈沖個數(shù)的增加，拉應力逐漸增大。沿軸線方向，拉應力逐漸減小并逐漸轉化為壓應力，壓應力在軸向上先逐漸增大，并在z=1mm附近處達到最大值，然后趨緩逐漸降低，并在距離后表面0.5mm附近處迅速減小并逐漸轉化為拉應力，拉應力在1064nm增透熔石英后表面中心位置達到最大值。

在軸線方向上，軸向應力始終表現(xiàn)為壓應力，且隨脈沖個數(shù)的增加壓應力逐漸增大。軸向應力沿軸向位置先逐漸增大并在軸向中心位置處達到最大值，然后逐漸減小，并在后表面中心位置處減小為零。根據(jù)1064nm增透熔石英的應力特性，可以得到1064nm增透熔石英在t=201ms時超過應力屈服強度區(qū)域對應的軸向深度z=4mm。

圖12 徑向/軸向應力隨軸向位置的變化關系

3 結論

本文根據(jù)熱傳導方程和彈性力學方程，在假設熔石英及其表面膜層材料為各向同性的連續(xù)介質條件下，建立了激光損傷1064nm增透熔石英隨時間變化的溫度場和應力場的理論模型。采用數(shù)值模擬方法，計算并分析了單脈沖及脈沖串條件下，激光能量密度、脈沖寬度和重復頻率的變化對1064nm增透熔石英溫度和應力的影響，并分別對其損傷效果和損傷機理進行了探討。論文研究結果可為激光損傷1064nm增透熔石英的理論和數(shù)值模擬研究提供依據(jù)，同時也為長脈沖激光對熔石英材料的損傷機理分析奠定了基礎。

［1］付耀龍，孫振皓，徐立君.長脈沖激光輻照硅材料熱應力的數(shù)值模擬［J］.長春理工大學學報：自然科學版，2013，36（3-4）：107-109.

［2］Dai G，Chen Y B，Shen Z H，et al.Analysis of laser induced thermal mechanical relationship of HfO2/SiO2high reflective optical thin film at 1064nm［J］.Chinese Optics Letters，2009，7（7）：601-604.

［3］Yoshida K，Umemura N.Influence of beam spatial distribution on the laser damage of optical material［J］.J Appl Phys，2006，49（8）：3815-3819.

［4］孫承偉.激光輻照效應［M］.北京：國防工業(yè)出版社，2002：310-314.

［5］胡漢平.熱傳導理論［M］.合肥：中國科學技術大學出版社，2010.

［6］戴罡，陸建，劉劍，等.長脈沖激光致兩種典型光學薄膜損傷閾值的測試與分析［J］.測試技術學報，2011，25（2）：122-127.

［7］宋林森，史國權，李占國.ANSYS在激光打孔溫度場仿真中的應用［J］.長春理工大學學報：自然科學版，2006，29（4）：19-21.

［8］徐秉業(yè).簡明彈塑性力學［M］.北京：高等教育出版社，2011：18-23.