基于Abaqus的放射性光源激光封割技術仿真研究

李思杰，馬俊平，平杰紅，唐 顯，孫玉華

(中國原子能科學研究院，北京 102413)

放射性光源是將放射性氣體密封在涂有熒光層的玻璃管內，通過放射性氣體(通常為氣態氚)衰變發出的β粒子激發熒光物質原子而發光，具有使用壽命長、光強穩定、無需電源、無需維護等優點[1-4]。其制備流程主要包括玻璃成型、熒光粉涂覆、填充放射性物質及封割四個步驟。其中，封割是指將內充有放射性氣體的密封長玻璃管分割成合適長度的放射性光源，在分割的同時需對端口進行密封，以防止放射性氣體泄漏，因此，封割是關鍵步驟之一，封割效果決定放射性光源的成品與質量。傳統的放射性光源采用火焰封割，此方法依賴操作人員的經驗與熟練程度，放射性光源成品的美觀性和成品率都難以得到保證。和傳統的火焰封割方法相比，激光技術更適用于放射性光源的封割，其能量密度大，熱影響區域小，封割速度快，且更精確可控、安全。對放射性光源封割而言，通常采用的是CO2激光，波長為10.6 μm，比較容易被玻璃吸收，更適用于玻璃材料的加工處理[5-6]。

激光封割放射性光源是利用激光的高能量密度加熱放射性光源玻璃管，在相應的激光功率密度下，激光束的能量被玻璃管吸收，引起激光作用點的溫度急劇上升，達到熔融點后，玻璃管封割點開始軟化，隨著激光束和玻璃管的相對運動，在離心力、內外壓差和熱應力等多因素共同作用下，完成對玻璃管的封割。在進行封割技術研究時，為減少實驗次數和放射性廢物的產生，利用有限元方法進行仿真，以理論指導實驗；而Abaqus是一種功能強大的基于有限元方法的工程模擬軟件，它不僅能解決結構分析(應力/位移)問題，還能模擬和研究熱傳導、質量擴散等廣闊領域的問題[7-8]。本研究借助有限元分析軟件Abaqus，建立放射性光源激光封割的仿真模型，分析激光對放射性光源玻璃管的影響，擬為最終確立放射性光源激光封割的條件參數提供指導。

1 仿真模型

應用Abaqus有限元分析軟件對放射性光源激光封割進行仿真分析。首先建立仿真模擬的幾何模型，然后建立傳熱-應力分析物理模型，提交任務進行計算，最后對計算結果進行處理。

激光封割放射性光源的運動示意圖示于圖1。放射性光源玻璃管以一定速度旋轉，一束激光輻照到玻璃管上，玻璃管溫度上升，激光束位置的玻璃熔融，激光束停止輻照，實現放射性光源玻璃管分開并密封。

圖1 激光封割放射性光源運動示意圖Fig.1 Schematic diagram of laser sealing and cutting tritium light source

玻璃管外徑為1.5 mm，壁厚為0.25 mm，建立如圖2的計算模型及網格，激光作用區域及其附近有限元單元網格劃分較密，而其余區域網絡相對稀疏，單元類型為C3D8T。

圖2 仿真用幾何模型及其網絡Fig.2 Geometric model and its grid for simulation

激光輻照到旋轉的玻璃管上，其溫度場為無內熱源的非線性瞬態熱傳導方程[9]，表示為：

(1)

式中，T為材料結構溫度，K；t為時間，s；c為材料的定容比熱，J/(kg·K)；ρ為材料的密度，kg/m3；λ為導熱系數，W/(m·K)；r、θ、z為極坐標的三個坐標變量。

由于溫度迅速變化和溫差較大導致的熱膨脹情況不一致，因此，在激光封割過程中，放射性光源玻璃管材料內部產生較大的應力，即熱應力。熱應力與溫度和材料的膨脹系數關系為[10]：

αT(T0)(T0-Tref)

(2)

式中，αT為熱膨脹系數，1/K；Tref為參考溫度，K；T0為初始溫度，K。激光直接作用于表面薄層，對其進行加熱，然后再通過熱傳導間接作用于材料內部，因此，在應用Abaqus對激光輻照下放射性光源玻璃管的溫度和應力分布進行模擬時，激光熱源為高斯分布，并隨著進入材料的深度而線性減小，其分布為[3]：

(3)

式中，QL為激光功率，W；r0為束斑半徑，mm；d為激光熱源作用深度，mm；r為瞬態激光束位置，其中，r與x、y、z關系為：

r2=x2+y2+z2

(4)

x、y、z為t時刻激光束坐標。在封割的過程中，激光熱源固定，玻璃管轉動；但在進行仿真模擬時，由于是分析溫度場變化及相應的應力分布情況，因此不能直接施加玻璃管轉動作為邊界條件，為了模擬實際狀態，在Abaqus中將激光熱源作繞圓周運動，實現等效封割條件。

激光在沿圓柱面輻照的過程中，激光熱源作為函數型熱通量進行加載，該熱通量通過用戶子程序施加到圓周上，公式(1)～(4)構成了求解激光熱傳導方程組。在用戶子程序中，束斑直徑設置為1.5 mm，轉速為300 r/min，玻璃管壁厚即為激光熱源的滲透深度。

由于激光照射加熱時間極短，模型中熱邊界條件并未考慮對流和輻射傳熱的影響，放射性光源所用玻璃載體的泊松比為0.2，相關熱導率隨著溫度變化的相關參數設定列于表1。

表1 熱導率參數設定Table 1 Thermal conductivity parameter setting

2 仿真結果與討論

運用已建立的仿真模型，得到激光封割的仿真結果。激光功率為10 W時，放射性光源玻璃管上設置固定監測點的溫度場在玻璃管徑向上隨著時間的變化趨勢示于圖3。從圖3結果可以看出，監測點的溫度隨著封割時間的增加而規律性地迅速升高，直至達到玻璃管的熔融溫度，當t=1.15 s時，溫度達到了約1 500 K(該溫度適于放射性光源玻璃管的封割)。

圖3 監測點溫度隨著時間變化Fig.3 Temperature Vs. time of the monitoring point

在t=1.15 s時刻玻璃管的軸向溫度云圖分布結果示于圖4。從圖4結果可以看出，激光中心溫度最高，然后沿軸向兩側迅速下降，呈高斯分布；光斑中心所在的整個圓周上，溫度明顯高于周圍區域，因此，激光輻照在放射性光源玻璃管圓周上，可以實現整管同時熔融。

圖4 t=1.15 s時刻的軸向溫度場分布Fig.4 Axial temperature field distribution at time t=1.15 s

由于封割過程中激光能量高度集中，在玻璃管上的能量沉積亦是高度集中，因此在放射性光源玻璃管上會出現較大的溫度差異，造成熱沖擊，導致玻璃材料內部出現較大的熱應力。監測點應力隨著封割時間的變化趨勢示于圖5。由圖5結果可知，在t=0.03 s時便出現了最大的應力，且在此過程中可以看到，隨著激光輻照時間增加，由于熱傳遞的作用，中心區域溫差變小，應力逐漸向圓周兩側移動，但整個過程中最大熱應力均小于放射性光源玻璃的破壞應力90～110 MPa，因此，激光封割過程中玻璃管不會出現炸裂的情況。

分別對在放射性光源封割過程中的不同激光功率進行仿真分析，得到達到封割溫度時所需的最短時間及其溫度場分布，結果示于圖6。從圖6結果可以看出，隨著激光功率增大，封割放射性光源玻璃管所需的時間減小，當激光功率達到20 W以后，達到滿足封割條件所需的時間差別較小，說明此時為高斯分布的激光熱源，進入材料的滲透分布深度已接近放射性光源玻璃管的壁厚。

圖5 監測點應力隨時間變化Fig.5 Stressat the monitoring point Vs. Time

圖6 不同激光功率下的封割時間及其溫度場分布Fig.6 Sealing time and temperature distribution with different laser powers

3 實驗驗證

利用現有的激光封割系統，對仿真結果進行驗證，結果示于圖7。由圖7結果可以看出，仿真值與實驗值變化趨勢相同，在某個激光功率條件下，激光輻照封割時間計算值比實驗值略短,這是由于在仿真模擬激光封割的計算模型中，忽略了激光傳輸過程引起的氣體熱效應、周圍環境對流散熱以及放射性光源玻璃管的輻射散熱等因素?？傮w來說，仿真結果能較真實地反映激光封割放射性光源玻璃管的情況。

圖7 激光功率與激光輻照封割時間關系Fig.7 Laser power and Vs. sealing time

4 結論

本文利用有限元軟件Abaqus，對放射性光源的激光封割技術進行仿真，得到了旋轉的放射性光源玻璃管在激光輻照下隨著時間變化的溫度分布及熱應力分布，以及不同激光功率與封割時間的關系，隨著激光功率增大，所需的封割時間越小，最后通過實驗驗證，得到封割的仿真值與實驗值較為一致。本文的模擬結果和分析首次為放射性光源研制的激光封割技術提供了理論指導和實驗驗證，可為最終確立放射性光源研制的重要關鍵技術參數提供參考。