飛秒激光燒蝕血管支架材料的數值模擬

王雷雷，鄭洋洋，霍 揚，劉建英，李春霞，盧 洋

(山東理工大學機械工程學院，山東 淄博 255049)

?

飛秒激光燒蝕血管支架材料的數值模擬

王雷雷，鄭洋洋，霍 揚，劉建英，李春霞，盧 洋

(山東理工大學機械工程學院，山東 淄博 255049)

王雷雷

為了研究飛秒脈沖激光燒蝕血管支架材料的特性，利用考慮了電子之間熱傳導項的雙溫模型，采用有限差分法，對飛秒激光燒蝕NiTi合金的溫度場分布進行數值模擬，計算得到了電子溫度和晶格溫度隨時間和燒蝕深度的變化規(guī)律，進一步討論了不同激光能量密度、不同激光脈寬、不同延遲時間對電子和晶格的溫度場影響。發(fā)現血管支架材料在飛秒激光的作用下，先是電子吸收能量溫度快速升高，再通過電聲耦合作用將能量傳遞給晶格，最后兩者的溫度達到一個平衡狀態(tài)；激光能量密度主要影響電子的峰值溫度和電子與晶格的平衡溫度；脈沖寬度主要影響電子的峰值溫度和達到峰值溫度所用的時間；電子溫度隨著延遲時間的增加先升高后降低，晶格溫度隨著延遲時間的增加不斷上升。這些理論分析對實際飛秒激光加工血管支架有重要的指導意義。

飛秒激光；鎳鈦合金；血管支架；雙溫方程；電子溫度；晶格溫度

1 前 言

NiTi合金材料是形狀記憶金屬中最重要且應用最為廣泛的合金，具有優(yōu)異的形狀記憶效應和超彈性、良好的耐疲勞性能和抗腐蝕性，較強的生物相容性，因此成為一種較為理想的生物工程材料，在現代醫(yī)學介入醫(yī)療器械領域得到了廣泛應用[1-3]。其中NiTi合金血管支架的研究和制備具有重要的應用價值，它的基本原理是支架植入血管后，依靠材料的形狀記憶功能自動恢復形狀，把狹窄血管撐開，恢復遠端血管供血，在大規(guī)模的實際臨床應用中取得了很好的治療效果[4]，具有取代傳統不銹鋼血管支架的趨勢。由于支架的尺寸微小，一般尺寸小于100 μm，如果采用普通的機械加工方法，很難滿足其加工要求；光纖激光器等普通的激光加工會對支架產生熱影響區(qū)和重鑄層[5]，需要電化學拋光等后續(xù)處理，給加工工藝帶來了一定難度；而飛秒激光加工具有熱影響范圍小、切縫周圍粗糙度小、切口平行度好、定位精度高、工件的形變小、無機械應力等特點，可精確切割形狀復雜的特殊微細零件[6,7]，因此利用飛秒激光加工NiTi合金的血管支架具有獨特的優(yōu)勢。目前對飛秒激光與金屬作用的理論研究主要集中在銅、鋁、金等材料，實驗研究主要有打孔、劃線、表面周期結構等。關于飛秒激光用于血管支架加工的研究主要有：新加坡的Bryan和德國漢諾威激光中心的Ostendorf利用飛秒激光加工出的血管支架得到了相關學者的矚目[8]；程萍等人利用飛秒激光加工了以聚乳酸為材料的可降解血管支架，證明了飛秒激光在此方面的優(yōu)勢[9]；而國內關于利用飛秒激光加工金屬材料血管支架的研究卻很少，并且飛秒激光與金屬材料的作用過程受到激光器的選取參數、加工工藝環(huán)境和材料的特性等眾多因素的影響，其內在加工的機理也比較復雜，目前還沒有統一的定論來解釋飛秒激光與材料的相互作用。因此為了高效利用飛秒激光，使其在加工血管支架時達到更理想的實際效果，有必要從微觀角度理論分析飛秒激光加工工藝參數對靶材內部的溫度影響規(guī)律。

從雙溫方程入手，采用有限差分法，以NiTi合金為例，模擬了飛秒激光燒蝕血管支架材料的溫度場分布情況；同時探究了不同參數對電子和晶格溫度的影響規(guī)律。

2 計算模型

2.1 雙溫方程

前蘇聯學者Anisimov S I根據飛秒激光對金屬作用的特殊機理，在1974年提出了適用于飛秒激光和金屬作用的雙溫模型[10]。它是用來闡述飛秒激光和金屬材料相互作用的數學模型，考慮分析了燒蝕過程中光子和電子、電子之間、電子和晶格之間以及晶格之間的一系列相互作用[11]，考慮了整個過程中的能量平衡，列出了電子和晶格溫度變化的微分方程，其表達式為：

(1)

(2)

式中C表示熱容，T表示溫度，下角標e和i分別代表電子和晶格系統。方程(1)右邊第一項為電子的熱傳導項，ke為電子的熱傳導率；第二項為電子和晶格的能量耦合，g為電子晶格耦合系數；第三項S為激光熱源。方程(2)的右邊代表作用過程中電子和晶格的能量耦合，因為晶格和晶格的耦合作用極小，可以忽略，所以對于晶格而言只與電子有能量的傳遞。由方程(1)和(2)可以看出，超短激光脈沖燒蝕金屬材料時，首先是電子吸收激光光子的能量，然后通過電子與晶格之間的耦合作用，使晶格溫度升高，產生燒蝕，達到精細加工目的。

很多學者在研究雙溫方程時，經常會忽略電子之間的熱傳導項或者把電子的熱傳導率看作一個常數，而實際上當飛秒激光照射到金屬靶材表面時，電子將被加熱到很高的溫度，甚至接近或超過費米溫度。此時電子熱傳導率的表達式應為[12]：

(3)

式中μe=Te/TF,μi=Ti/TF，NiTi合金的費米溫度TF=8.12×104K,參數χ和η是和被加工材料相關的參數，分別為336 W·K-1·m-1和0.178。

飛秒激光的光源為高斯光束，因此光源項S(x,t)的表達式為[13]：

(4)

式中A代表消光系數，α為材料的吸收系數，F為激光的能量密度，τL為激光的脈沖寬度，x為距材料表面的軸向距離。

2.2 計算方法

對雙溫方程的數值求解采用有限差分法，差分采用時間向前空間中心差分法(FTCS)全隱格式，迭代法求解[14]。Δt和Δx分別為時間步長和軸向步長，i和j分別為軸向步數和時間步數，由此得雙溫模型的差分表達式為式(6)和式(7)：

(6)

(7)

(8)

(9)

exp(-αiΔx)

初始條件和邊界條件滿足：

Ti,0=T0=300K,(i=1,2,3,…,n)

Tn,j=300K,(j=0,1,2,…,m)

2.3 采用參數

模擬采用的飛秒激光器是美國Raytum Corp生產的二極管直接泵浦超快激光器。其參數為：激光能量，F=100 J/m2；脈寬，τL=300 fs；波長為1030 nm。

取時間差分步長0.1 fs，空間差分步長為1 nm，即每0.1 fs計算一次，每1 nm分為一層。NiTi合金的相關物理參數見表1。

表1 NiTi合金的參數[15]

3 計算結果與分析

圖1a和1b給出了雙溫方程模擬的電子和晶格的時間-空間溫度場?？梢钥吹剑陲w秒激光作用下，電子溫度在幾百飛秒內迅速升高，高溫電子通過碰撞進行能量傳遞，使其系統的峰值溫度達到3050 K，再通過電子熱擴散將激光的能量傳遞到靶材內部，而此時由于晶格不能直接吸收激光能量，其溫度仍比較低。然后電子通過劇烈的電聲耦合作用將能量快速地傳遞給晶格，使晶格接受能量溫度不斷升高，電子溫度不斷下降，最終兩者的溫度都在1300 K左右，達到熱平衡狀態(tài)，整個能量耦合的非平衡過程在10個皮秒左右完成[16,17]。就傳遞深度而言，電子在峰值溫度時傳遞深度約為20 nm，晶格在平衡狀態(tài)時的傳遞深度也在20 nm左右，可見飛秒單脈沖的作用只發(fā)生在材料表面，對材料內部的影響并不明顯。由此可見，飛秒激光加工材料時，在超短時間內通過極高的激光功率作用于工件，在飛秒時量級時間內，使材料從固態(tài)瞬間變?yōu)楦邷亍⒏邏旱牡入x子體狀態(tài)，以噴射的形式脫離被加工材料，實現了材料的非熱熔加工。因為脈沖的持續(xù)時間極短，溫度來不及向周圍擴散就發(fā)生了燒蝕，因此熱影響區(qū)和熱擴散的作用也很小，所以利用飛秒激光可以實現血管支架的精細加工。

圖1 電子溫度(a)和晶格溫度(b)分布Fig.1 Distribution of electrons temperature (a) and lattices temperature (b)

在其他參數一定的情況下，選擇不同的脈沖能量密度進行數值模擬，得到電子與晶格溫度場變化規(guī)律如圖2所示。圖中顯示隨著脈沖能量密度的增大，電子的峰值溫度會隨之增大，在脈沖過后電子的溫度也會增加，最終使得電子和晶格的平衡溫度也有所升高。達到耦合狀態(tài)的標志是晶格和電子的溫度達到一致，不隨時間的延長而變化，在圖上觀察到的現象就是曲線開始變得平緩。在圖2a和2b中可以看出能量為100 J/m2時，4 ps左右曲線變得平滑，也就是說此時達到了耦合狀態(tài)。800 J/m2時，直到9 ps左右曲線才開始平緩，標志著此時電子和晶格的溫度一致。隨著能量密度的增大，耦合時間延長。如果激光的能量密度繼續(xù)升高，當電子和晶格耦合時間超過晶格和晶格的熱平衡時間，熱擴散就會對燒蝕周圍區(qū)域的材料造成影響，燒蝕的過程將會是“熱”作用。由此推論對于不同的脈寬應該有一個臨界的能量密度值，當激光的能量密度小于此值為相對意義上的“冷”加工，超過這個值將產生一定熱作用，將對加工區(qū)周圍材料和加工后的工件質量產生一定影響。因此為了進一步提高加工質量和加工精度，首先應該找到這個臨界的能量密度值。這個結論對實際應用飛秒激光加工血管支架有重要的理論指導。

選取不同的脈沖寬度進行數值模擬，得到電子與晶格溫度場變化規(guī)律如圖3所示。從圖中可以看出，能量密度相同的情況下脈沖寬度越小，電子溫度上升的速度越快，達到的峰值溫度越大，所用時間也越短，而最終的電子溫度卻變化不大。對于晶格而言，脈寬越小，晶格的溫度變化越快，和電子的耦合時間越短，最終的平衡溫度基本相同，與文獻[18]的結果相近。上述結果與飛秒激光燒蝕的時間有關，激光的脈沖寬度越窄，靶材的電子在極短時間內吸收的激光能量越大，所以電子溫度快速升高，達到的峰值溫度也越大，電子通過電聲耦合作用把能量傳遞給晶格，因此晶格溫度上升的速度也會變快，達到平衡狀態(tài)時所需要的時間越短。實際的激光器脈寬為300 fs，通過以上分析可知，該參數可以實現血管支架的精密加工。

圖2 不同激光能量密度下，電子(a)和晶格(b)的溫度分布Fig.2 Distribution of electrons temperature (a) and lattices temperature (b) under different laser energy densities

圖4分析了延遲時間分別為50 fs，150 fs，l ps和6 ps時，NiTi合金材料不同深度處的電子和晶格的溫度變化。

圖3 不同激光脈沖寬度下，電子(a)和晶格(b)的溫度分布Fig.3 Distribution of electrons temperature (a) and lattices temperature (b) under different laser pulse width

圖4 不同延遲時間下，電子(a)和晶格(b)的溫度分布Fig.4 Distribution of electrons temperature (a) and lattices temperature (b) with different delay time

從總體來看，隨著深度的不斷增加，電子和晶格溫度均呈指數衰減。延遲時間在150 fs以內時，表層電子溫度在不斷上升，到150 fs溫度達到最大值約為2500 K，此后由于電子和晶格的耦合作用，電子的溫度不斷下降，6 ps左右達到最小值約為700 K。與此同時，隨著電聲耦合作用晶格溫度不斷上升，最終在6 ps左右和電子達到平衡狀態(tài)，表明此時整個耦合過程已經完畢。距表面20 nm以下溫度沒有變化，說明整個熱傳導過程只涉及到材料表面下20 nm處。因此如果合理控制飛秒激光的加工參數，可以實現NiTi合金血管支架的高精度加工。

4 結 論

為了研究飛秒激光燒蝕血管支架材料的溫度場分布，以NiTi合金為例，利用雙溫模型對飛秒激光燒蝕血管支架材料進行了數值模擬。得到了電子和晶格隨時間和燒蝕深度的溫度變化規(guī)律，并討論了激光能量密度、脈沖寬度、延遲時間對電子和晶格的溫度變化影響。計算結果表明：①在飛秒激光作用下，首先是靶材的電子溫度急劇升高，然后通過電子和晶格的耦合作用，將能量傳遞給晶格，最終兩者達到熱平衡狀態(tài)；②激光能量密度越大，電子達到的峰值溫度越大，最后電子和晶格的熱平衡溫度也越高；③激光的脈沖寬度越小，電子和晶格的溫升越快，達到平衡狀態(tài)的時間越短；④隨著燒蝕深度的增加，電子和晶格的溫度都呈指數遞減，熱傳導深度僅為幾十納米。通過以上結論分析，得出了影響飛秒激光燒蝕NiTi合金血管支架材料過程的幾個因素，對提高飛秒激光加工NiTi合金血管支架的精度和效率有重要的意義。

References

[1] Dotter C T, Buschmann R W, McKinney M K,etal.Radiology[J], 1983, 147(1): 259-260.

[2] Wu Yuanhao(吳遠浩),Zhou Xiaochen(周曉晨),Li Nan(李楠)，etal.MaterialsChina(中國材料進展)[J],2012,31(9):27-34.

[3] Meng Xianglong(孟祥龍),Cai Wei(蔡 偉).MaterialsChina(中國材料進展)[J],2011,30(9):13-20.

[4] Muhammad N, Whitehead D, Boor A,etal.AppliedPhysicsA[J], 2012, 106(3): 607-617.

[5] Huang Weidong(黃衛(wèi)東),Lv Xiaowei(呂曉衛(wèi)), Li Xin(林鑫)，etal.MaterialsChina(中國材料進展)[J],2011,30(4):1-10.

[6] Gu Xingzhong(顧興中), Ni Zhonghua(倪中華).JournalofHuazhongUniversityofScienceandTechnology(NatureScienceEdition)(華中科技大學學報(自然科學版))[J],2007,35(z1):143-146.

[7] Guo Jingzhe(郭景哲),Jia Baoxian(賈寶賢), Wang Dongsheng(王冬生),etal.MachineryDesign&Manufacture(機械設計與制造)[J], 2009 (10): 257-258.

[8] Zhou Yongheng(周永恒), Liao Jianhong(廖健宏), Meng Hongyun(蒙紅云),etal.AppliedLaser(應用激光)[J],2005,25(3):161-164.

[9] Cheng Ping(程 萍), Wei Di(位 迪), Wu Benke(吳本科),etal.OpticsandPrecisionEngineering(光學精密工程) [J],2014,22(1):63-68.

[10]Anisimov S I, Kapeliovich B L, Perelman T L.ZhEkspTeorFiz[J], 1974, 66(2): 375-377.

[11]Kim J, Na S.Optics&LaserTechnology[J], 2007, 39(7):1443-1448.

[12]Chen Anmin(陳安民), Jiang Yuanfei(姜遠飛), Liu Hang(劉航),etal.Laser&Infrared(激光與紅外)[J],2012,42(8):847-851.

[13]Elsayed-Ali H E, Norris T B, Pessot M A,etal.PhysicalReviewLetters[J], 1987, 58(12): 1212-1215.

[14]Ni Xiaochang(倪曉昌), Wang Qingyue(王清月).ChineseJournalofLasers(中國激光)[J],2004,31(3):277-280.

[15]Zhu Shijia(朱詩佳), Wang Xinlin(王新林), Zhu Weihua(朱衛(wèi)華),etal.Electro-OpticTechnologyApplication(光電技術應用) [J],2011,26(6):39-44.

[16]Yao Hongbing(姚紅兵),Yu Wenlong(于文龍),Yang Zhao(楊 昭),etal.ActaPhotonicaSinica(光子學報) [J], 2015,44(4):35-40.

[17]Wang Wenting(王文亭), Zhang Nan(張 楠), Wang Mingwei(王明偉),etal.ActaPhysicaSinica(物理學報) [J], 2013 (21):87-95.

[18]Lu Lizhong(盧立中), Shi Yunfei(石云飛), Xu Chenguang(徐晨光),etal.JournalofJiangsuUniversity(NaturalScienceEdition)(江蘇大學學報(自然科學版))[J], 2011,32(2):199-204.

(編輯 惠 瓊)

Numerical Simulation of Femtosecond Laser Ablation of Stent Material

WANG Leilei，ZHENG Yangyang，HUO Yang，LIU Jianying，LI Chunxia， LU Yang

(School of Mechanical Engineering, Shandong University of Technology, Zibo 255049, China)

In order to study the characteristics of femtosecond laser ablation of stent material, two-temperature model which considers the thermal conduction between the electronics, and finite difference method were used to simulate temperature field of NiTi alloys during femtosecond laser ablation. According to the results, the influences of laser energy density, pulse width and delay time on the temperature field of electron and lattice were discussed. The results showed that the temperature of electron and lattice finally reached a state of equilibrium. Before that, firstly, the temperature of electron increased rapidly due to energy absorption, then, the absorbed energy was transferred to lattice through electron phonon coupling effect; at last, the equilibrium would be reached. Laser energy density played a major role in the equilibrium temperature, and pulse width determined the peak temperature and the time to reach it. Consequently，the temperature of electrons first increased and then decreased with the increase of the delay time, but lattice temperature only increased with the increase of the delay time. The results may have an important guiding significance to actual processing of femtosecond laser ablation of vascular stent.

femtosecond laser; nickel titanium alloy; vascular stent; two-temperature model; electron temperature; lattice temperature

2015-09-22

國家自然科學基金項目(61303006);山東省優(yōu)秀中青年科學家科研獎勵基金項目(BS2012ZZ009)。

王雷雷，男，1990年生，碩士研究生, E-mail:354019062@qq.com

10.7502/j.issn.1674-3962.2017.03.10

TN249

A

1674-3962(2017)03-0226-05