蒙宣伊
(湘電集團(tuán) 風(fēng)能有限公司,湘潭 411102)
激光熔覆是一種新型的材料加工與表面改性技術(shù),它可以在低成本材料上制成高性能表面,提高材料的綜合性能,降低能源消耗,該技術(shù)引起了西方各國的高度關(guān)注[1-5].
Yunchang Fu,A.Loredo等建立了熔覆過程的解析模型[6],示意圖如圖1所示.
圖1 激光熔覆模型示意圖
目前已有模擬溫度場和應(yīng)力場的文獻(xiàn),但均沒有考慮激光束與粉末的交互作用,因此,研究激光熔覆過程的應(yīng)力產(chǎn)生機(jī)理,從而提出減少和消除殘余應(yīng)力的方法,對(duì)于激光熔覆技術(shù)的完善和推廣具有重要實(shí)際意義.
基體Q235力學(xué)性能參數(shù)見表1,熔覆層材料力學(xué)性能參數(shù)如表2所示.
表1 基體Q235力學(xué)性能參數(shù)
表2 熔覆層材料力學(xué)性能參數(shù)
本文中環(huán)境溫度為293 K,工件無初始應(yīng)力,采用對(duì)稱模型,位移和對(duì)稱約束見圖2,這樣讓約束點(diǎn)盡量分散,且遠(yuǎn)離光斑加熱區(qū),避免計(jì)算產(chǎn)生偏差.
本文在模擬激光熔覆的熱應(yīng)力場時(shí),采用生死單元法實(shí)現(xiàn)熔覆層的生長過程,在熔覆開始前把熔覆單元“殺死”,并在每一步熱應(yīng)力計(jì)算時(shí),將對(duì)應(yīng)溫度場的計(jì)算結(jié)果加載到應(yīng)力場模型中.
同時(shí)采用改變單元屬性法,在開始計(jì)算前,先定義熔融金屬的材料屬性(熔融金屬包括處于熔化狀態(tài)的熔覆粉末和基體熔池金屬)材料編號(hào)為3,使它們的屈服極限和彈性模量都很低,且不隨溫度變化,熔覆材料屬性編號(hào)為2,基體材料編號(hào)為1.熔覆材料和基體的力學(xué)性能參數(shù)都隨溫度變化.
圖2 位移與對(duì)稱約束
在程序計(jì)算過程中,首先改變結(jié)構(gòu)單元為熱單元,在溫度場后處理器中讀取節(jié)點(diǎn)溫度值,將超過熔點(diǎn)的熔覆層單元和基體單元分別存入溫度表中.再轉(zhuǎn)換熱單元為結(jié)構(gòu)單元,并改變這些單元材料屬性.
圖3 應(yīng)力隨時(shí)間變化曲線圖
圖3 為熔覆層中心某點(diǎn)應(yīng)力隨熔覆時(shí)間變化曲線圖,圖中SZ表示縱向應(yīng)力,SX表示橫向應(yīng)力,SY表示厚向應(yīng)力.
當(dāng)時(shí)間t=6.667 s時(shí),光斑照射到該點(diǎn)時(shí),該點(diǎn)所在的單元被激活,并迅速達(dá)到熔點(diǎn)溫度,應(yīng)力場計(jì)算近似等于零;t=7.333 s,光斑移出該點(diǎn),該點(diǎn)很快降至熔點(diǎn)溫度,并開始凝固.
材料冷卻收縮,受到基材或自身溫度不均勻的限制,處于拉應(yīng)力狀態(tài).在塑性范圍內(nèi),拉應(yīng)力縱向最大,橫向次之,厚度方向最小.拉應(yīng)力的數(shù)值取決于材料的流動(dòng)應(yīng)力,隨溫度的下降,流動(dòng)應(yīng)力上升,拉應(yīng)力表現(xiàn)為上升趨勢(shì),當(dāng)材料進(jìn)入彈性狀態(tài)后,橫向和厚度方向上的拉應(yīng)力下降明顯,縱向應(yīng)力基本保持不變.
圖4 殘余應(yīng)力分布云圖
圖3 (a)為冷卻600 s后的縱向殘余應(yīng)力分布云圖,由圖看出,最大縱向殘余應(yīng)力σzmax=978 MPa,且靠近基體的熔覆層上,這是由于熔覆層材料的屈服強(qiáng)度、切變模量均高于基體,且熔覆層的溫度梯度大,所以導(dǎo)致殘余應(yīng)力較大.圖3(b)為冷卻600 s后橫向殘余應(yīng)力分布云圖,橫向殘余應(yīng)力最大值出現(xiàn)在熔覆層邊界與基體的交界處,σxmax=387 MPa.
圖5 至圖7分別為熔覆層中心點(diǎn)橫向、厚向和縱向的應(yīng)變變化情況,包括熱應(yīng)變(Thermal)、彈性應(yīng)變(Elastic)、塑性應(yīng)變(Plastic)和總應(yīng)變(Total)總應(yīng)變等于熱應(yīng)變、彈性應(yīng)變、塑性應(yīng)變?nèi)咧?
在時(shí)間t=6.667 s到t=7.333 s之間,該點(diǎn)在激光光斑內(nèi),處于熔化狀態(tài),因此出現(xiàn)應(yīng)變突然增大的情況.當(dāng)激光光斑移除該點(diǎn)后,該點(diǎn)溫度迅速降低,應(yīng)變都迅速減小;進(jìn)入彈性狀態(tài)后,各向總應(yīng)變均趨于平緩.
由圖5至圖6可以看出,各向應(yīng)變中,熱應(yīng)變最大,且總大于零;x、y向塑性應(yīng)變小于零,z向塑性應(yīng)變大于零.無論彈性應(yīng)變和塑性應(yīng)變符號(hào)怎么變化,總應(yīng)變(total)總大于零.說明,熱應(yīng)變是決定材料處于拉伸或壓塑塑性變形的決定因素;最終的橫向X方向總應(yīng)變小于零;厚度Y方向的應(yīng)變趨于零;縱向Z方向的應(yīng)變最大.
圖7 縱向應(yīng)變隨時(shí)間變化曲線
由于熔覆層心部組織冷卻速率要低于周圍組織,心部凝固收縮受到周圍組織的拘束作用.在三向應(yīng)力中,縱向應(yīng)力要遠(yuǎn)大于其它兩向應(yīng)力,最終表現(xiàn)為縱向產(chǎn)生拉應(yīng)變,其它兩向產(chǎn)生壓應(yīng)變,熔覆層中心的主應(yīng)變圖為第三類主應(yīng)變圖[7],如圖8所示.周圍組織對(duì)其有類似于擠壓和拉拔的作用.這進(jìn)一步說明了實(shí)際生產(chǎn)中熔覆層總是產(chǎn)生橫向裂紋這一事實(shí),熔覆層橫向裂紋如圖9所示.
圖8 第三類主應(yīng)變圖
圖9 熔覆層橫向裂紋
圖1 0和圖11為冷卻600 s后沿工件中心厚度方向從上至下縱向應(yīng)力和橫向應(yīng)力變化曲線圖.
由圖10看出,熔覆層(高1 mm)縱向應(yīng)力從上至下逐漸增大,到熔池附近達(dá)到最大值,這說明熔池附近是裂紋高發(fā)區(qū),容易產(chǎn)生橫向裂紋.過熔池后縱向應(yīng)力逐漸減小,且始終為拉應(yīng)力,對(duì)稱線上的拉應(yīng)力由兩側(cè)壓應(yīng)力補(bǔ)償.圖11表明,熔覆層上橫向應(yīng)力沿厚度方向從上至下是由壓變拉,到熔池附近拉應(yīng)力達(dá)到最大.基體則是上下兩側(cè)受拉,中間受壓.
激光熔覆過程產(chǎn)生的內(nèi)應(yīng)力是典型的由于材料溫度變化而形成的熱應(yīng)力.送粉式激光熔覆過程中,可以近似認(rèn)為所有形成熔覆層的粉末在達(dá)到基體之前全部熔化,到達(dá)基體后將熱量傳遞給基體,同時(shí)基體也直接受到激光束的照射,并在光斑內(nèi)形成熔池.處于液態(tài)的金屬流動(dòng)應(yīng)力較小,其受熱膨脹由于受到周圍材料限制所產(chǎn)生的不均勻的壓縮塑性變形和壓應(yīng)力可以忽略不計(jì).
而與熔池相鄰的高溫區(qū)(未熔化)產(chǎn)生熱膨脹,膨脹受到周圍材料的限制,產(chǎn)生不均勻的壓縮塑性變形和壓應(yīng)力;在冷卻過程中,高溫區(qū)一定程度上又被拉伸而卸載;同時(shí),熔池的冷卻凝固過程中,也受到周圍材料的限制,難以自由收縮,產(chǎn)生拉應(yīng)力與拉伸塑性變形;此外,由于金屬相變后體積發(fā)生變化,也會(huì)產(chǎn)生相應(yīng)的相變應(yīng)力;熔覆粉末與基體材料熱膨脹系數(shù)、彈性模量的差異,也會(huì)導(dǎo)致內(nèi)應(yīng)力的產(chǎn)生.
(1)熔覆層中心區(qū)域一直受到拉伸應(yīng)力,產(chǎn)生拉伸塑性變形;熔池附近的基體先受到壓應(yīng)力,產(chǎn)生壓縮塑性變形,然后隨著熔覆的進(jìn)行和冷卻再受到拉應(yīng)力,產(chǎn)生拉伸塑性變形.
(2)熔覆層縱向應(yīng)力最大,橫向應(yīng)力次之,厚度方向應(yīng)力最小.三向拉應(yīng)力當(dāng)材料在塑性狀態(tài)均表現(xiàn)為上升趨勢(shì),當(dāng)材料進(jìn)入彈性狀態(tài)后,橫向和厚度方向上的拉應(yīng)力下降明顯,縱向應(yīng)力基本保持不變.
(3)熔覆層中心靠近基體一側(cè)是裂紋敏感區(qū),該區(qū)域處于縱向拉伸、橫向和厚向壓縮的第三類主應(yīng)變狀態(tài),周圍組織對(duì)其有類似于拉拔和擠壓的作用.
總之,送粉式激光熔覆智能控制與數(shù)值模擬技術(shù)方興未艾,還有很多不確定因素需要探索,希望本文的研究工作能為該領(lǐng)域的發(fā)展起到拋磚引玉的作用.
[1]王仲任,郭殿儉,王 濤.塑性成形力學(xué)[M].哈爾濱工業(yè)大學(xué)出版社,1989:35-36.
[2]C.Lalas K.Tsirbas K.Salonitis G..Chryssolouris.Ananalytical Model of the Laser Clad Geometry[J].Int JAdv Manuf Technol,2007,32:34-41.
[3]Edson Costa Santos,Masanari Shiomi,Kozo Osakada.Rapid Manufacturing of Metal components by laser forming[J].International Journal of MachineTools&Manufacture,2006,46:1459-1468.
[4]Jichang Liu,Lijun Li.Study on Cross-section Clad Profile in Coaxial Single-pass Cladding with a lowpower laser[J].Optics&Laser Technology,2005,37:478-482.
[5]Kai Zhang,Weijun Liu,Xiaofeng Shang.Research on the Processing Experiments of Laser metal Deposition shaping[J].Optics&Laser Technology,2007,39:549-557.
[6]U.deOliveira,V.Ocelik,J.Th.M.DeHosson.Analysis of Coaxial Laser Cladding Processing Congditions[J].Surface&Coatings Technology,2005,197:130-135.
[7]Yunchang Fu,A.Loredo,B.Martin.A theoretical model for laser and powder particles interaction during laser cladding[J].Journal of Materials Processing Technology,2002,128:106-112.