,,
(浙江工業(yè)大學 機械工程學院,浙江 杭州 310014)
玻璃具有硬度高、透光率高和質(zhì)量輕等優(yōu)良特性,因而被廣泛應用于生產(chǎn)、生活中的各個方面,尤其是手持式便攜電子設備的顯示屏[1]。玻璃屬于脆性材料,在外力作用下僅發(fā)生很小的變形,當承受稍大于其彈性極限的力時就會發(fā)生斷裂,因此脆性材料不易加工[2]。傳統(tǒng)切割玻璃的方法是使用硬質(zhì)合金或金剛石刀具在玻璃表面進行劃線,然后采用機械外力折斷玻璃,這種方法需要復雜的后續(xù)工藝,在玻璃斷面處往往會產(chǎn)生微裂紋和殘余應力[3-4]。激光具有高能量密度、非接觸性和無污染等優(yōu)點,因而被應用于玻璃、陶瓷等脆性材料的加工中。激光切割按原理來分,主要有兩種:激光熔融法切割和激光熱裂法切割。激光熱裂法切割原理是激光使照射位置處的玻璃溫度瞬時升高,玻璃表面產(chǎn)生較大壓應力,之后玻璃溫度又迅速降低,產(chǎn)生較大的溫度梯度,使玻璃表面產(chǎn)生較大拉應力,當這個拉應力到達破壞極限時,玻璃會沿著激光掃描路徑開始斷裂[5]。
國內(nèi)外學者對激光熱裂法切割玻璃進行了大量研究,研究方向主要集中在應力場分布、斷面質(zhì)量等方面[6-7],對切割路徑偏移方面研究較少。 Nisar等[8]從激光熱源類型(體熱和面熱)、激光功率、材料厚度和切割速度等方面對熱裂法切割玻璃首尾端路徑偏移現(xiàn)象進行了研究,研究表明: CO2激光切割引起首尾端路徑偏移量比半導體激光引起的大;材料厚度越小、切割速度越慢,首尾端路徑偏移量越大。 Zhao等[9]結合 ABAQUS軟件對半導體激光非對稱切割鈉鈣玻璃進行了應力分析和裂紋動態(tài)擴展分析,得出切割路徑會向壓應力大的一側偏移。 Zhao等[10]對雙束激光熱裂法非對稱切割玻璃進行研究,研究表明雙束激光非對稱切割可以提升斷面質(zhì)量、降低路徑偏移量。筆者將從激光功率、光斑半徑、掃描速度和切割位置等方面,對激光熱裂法非對稱切割玻璃路徑偏移量的影響因素進行研究。
激光熱裂法切割過程中,玻璃的溫度隨時間變化,屬于非穩(wěn)態(tài)熱傳導問題,笛卡爾坐標下的熱傳導微分方程為
(1)
式中 :k為熱傳導系數(shù);T為溫度;ρ為密度;c為比
熱容;Q為玻璃內(nèi)部的熱源強度。
采用半導體激光切割玻璃,玻璃的吸收長度遠大于其材料厚度,因而建體熱高斯熱源模型[11]為
(2)
式中 :Ix,y,z,t為玻璃吸收的激光強度;η為吸收率;P為激光功率;r為激光光斑半徑;v為切割速度;δ(z)為吸收深度的倒數(shù)。
激光熱裂法非對稱切割玻璃裂紋類型為Ⅰ型和Ⅱ型復合裂紋[12]。當拉應力超過破壞應力時,裂紋開始擴展,表示為
(3)
式中 :E為彈性模量;γf為材料表面能;A為裂紋寬度。
本研究通過有限元軟件ABAQUS建立玻璃模型,分別用有限元方法對玻璃溫度場進行模擬,用擴展有限元方法對裂紋的動態(tài)擴展進行分析。仿真分析基于以下假設:1) 材料是均勻的,且各向同性;2) 玻璃內(nèi)部無殘余應力,且不考慮加熱過程的相變因素;3) 激光的能量完全符合高斯熱源模型;4) 外界的熱對流和熱輻射是可以忽略的。
基于以上,用FORTRAN語言編制激光移動熱源程序,用 ABAQUS建立尺寸為 30 mm×15 mm×2 mm的 DC3D8線性熱傳導玻璃模型,具體材料參數(shù)如表 1所示。為提高計算精度,對熱影響區(qū)域進行網(wǎng)格細化,如圖1所示。設定玻璃的初始溫度為 20 ℃,在玻璃模型上加載不同功率、移動速度和光斑半徑的激光移動熱源,以及在玻璃不同位置加載激光移動熱源,如圖2所示,分別得到各種條件下玻璃的溫度場。
表1 玻璃的參數(shù)Table 1 The parameters of glass
圖1 網(wǎng)格劃分Fig.1 Grid partition
圖2 不同切割位置Fig.2 Different cutting positions
Budyn等[13]在單位分解法基礎上建立了擴展有限元方法,在傳統(tǒng)有限元位移近似場中引入描述裂紋特征的擴充項,擴充后的位移近似場為
(4)
(5)
式中 :Nix為常規(guī)有限元形函數(shù);ui,aj,bkα分別為結點位移和結點加強變量;Ns為離散結構中所有結點的集合;Ncut為形函數(shù)支撐域完全被裂紋切割的結點集,用一個廣義的Heaviside函數(shù)Hx加強,位于裂紋上方取值為+1,位于裂紋下方取-1;Ntip為形函數(shù)支撐域含裂尖的結點集;r,θ分別為裂尖局部極坐標。
XFEM的基本原理是基于單位分解的思想在常規(guī)有限元單元法中位移模式中加入一些加強函數(shù)以反應不連續(xù)性, XFEM在分析不連續(xù)問題時計算網(wǎng)格和結構內(nèi)部的集合或物理界面是相互獨立的,因此能方便地分析不連續(xù)問題[14]。
本研究中使用 ABAQUS建立激光熱裂法非對稱直線切割玻璃的動態(tài)裂紋擴展模型,模型采用 C3DR8單元,初始裂紋設為 1 mm,初始裂紋與玻璃之間為硬接觸,斷裂準則選用最大周向拉應力準則,斷裂應力為 30 MPa[9],邊界條件設置為玻璃右側面的4 個端點全約束,在預定義場中使用溫度場分析的結果,以此來模擬不同條件下的裂紋擴展路徑以及應力分布情況。
不同的激光功率、移動速度和光斑半徑分別對應不同的溫度場,如圖3所示(圖3中:P為激光功率;r為光斑半徑;v為激光掃描速度; Δy為切割位置與對稱切割位置間的距離),激光功率越大、掃描速度越慢以及光斑半徑越小,對應熱流密度激光掃描位置處的熱流密度越大,從而玻璃溫度越高。
圖3 不同激光參數(shù)、不同掃描位置的溫度場分布Fig.3 Temperature distribution of different laser parameters and different scanning positions
激光非對稱切割導致切割路徑兩側玻璃的熱積累不同[9],引起切割路徑兩側應力分布不同,如圖4所示。從圖4中可以看出:裂尖位于激光光斑后方,裂尖兩側的材料表面存在不對稱的壓應力,且裂紋擴展路徑下方的壓應力比上方的壓應力更大。
圖4 玻璃表面應力分布Fig.4 Stress distribution of glass surface
圖5為激光非對稱切割路徑偏移的示意圖。圖6(a)為激光熱裂法非對稱切割玻璃的厚度方向,且垂直于切割路徑的裂尖位置的應力分布,明顯看出應力分布是不對稱的;從圖6(b)可以看出對稱切割時裂尖位置應力分布是對稱的。
圖5 路徑偏移示意圖Fig.5 Schematic of the deviation of the separation path
圖6 厚度方向裂尖應力云圖Fig.6 Stress distribution of crack tip along the thickness
圖7(a)表明路徑偏移量(指最大路徑偏移量)隨激光功率的增大而增大。因為隨著激光功率的增大,激光光斑照射位置的玻璃熱流密度變大,熱影響區(qū)域也會變大,從而引起掃描路徑兩側的應力分布差異變大,引起偏移量的增大。由圖7(a)線斜率可知:激光功率對偏移量的影響不是非常明顯。
圖7(b)表明隨著激光光斑半徑增大,路徑偏移量會先減小到一定值后,再增大。此現(xiàn)象由熱流密度的強弱和熱影響區(qū)域的大小共同決定,在一定程度上,激光光斑半徑越小,熱流密度越大,熱影響區(qū)域越小;反之,則熱流密度越小,熱影響區(qū)域越大。由圖7(b)可知:當激光光斑半徑小于 1.6 mm時,偏移量隨光斑半徑的增大而減小,即熱流密度對偏移量的影響占主導地位,熱影響區(qū)域次之;當激光光斑半徑大于 1.6 mm時,偏移量隨光斑半徑的增大而增大,即熱影響區(qū)域?qū)ζ屏康挠绊懻贾鲗У匚?,熱流密度次之?/p>
圖7(c)表明偏移量隨著激光掃描速度的增大而減小。這是因為激光掃描速度越快,熱流密度越小,從而引起掃描路徑兩側的應力分布差異變小,引起偏移量的減小。由圖7(c)線斜率可知:當掃描速度小于 7 mm/s時,偏移量的減小不明顯,當掃描速度大于 7 mm/s時,偏移量的減小比較明顯。
圖7(d)表明切割位置與對稱切割位置間的距離越大,偏移量越大。這是因為掃描路徑兩側的材料的尺寸差異越大,材料內(nèi)部間的擠壓變形量差異越大,從而引起掃描路徑兩側的應力分布差異越大,則偏移量越大。
圖7 切割參數(shù)對路徑偏移量的影響Fig.7 The effect of the cutting parameters on the deviation of separation path
激光熱裂法非對稱切割玻璃掃描路徑兩側應力呈不對稱分布,在玻璃表面上,材料尺寸大的一側的壓應力小于材料尺寸小的一側;在玻璃內(nèi)部,材料尺寸大的一側的拉應力大于材料尺寸小的一側,不對稱的應力分布導致切割路徑偏移。而在對稱切割情況下,切割路徑兩側應力分布是對稱的,不會產(chǎn)生路徑偏移現(xiàn)象。本研究結合斷裂力學的基本理論,通過有限元和擴展有限元分析,模擬激光熱裂法非對稱切割玻璃的溫度場和裂紋的動態(tài)擴展過程,研究了激光熱裂法非對稱切割玻璃路徑偏移的影響因素。研究結果表明: 1) 激光功率越大,路徑偏移量越大; 2) 激光光斑半徑小于 1.6 mm時,路徑偏移量隨光斑半徑的增大而減小,激光光斑半徑大于 1.6 mm時,路徑偏移量隨光斑半徑的減小而增大; 3) 掃描速度越大,路徑偏移量越?。?4) 與對稱切割位置間的距離越大,路徑偏移量越大。