許明三，江堯峰，姚耀伍，陳相檔，曾壽金，葉建華

SLM成形薄壁件尺寸偏差預(yù)測與控制研究

許明三a,b，江堯峰a,b，姚耀伍a,b，陳相檔a,b，曾壽金a,b，葉建華a,b

（福建工程學(xué)院 a.福建省智能加工技術(shù)及裝備重點(diǎn)實(shí)驗(yàn)室 b.機(jī)械與汽車工程學(xué)院，福州 350118）

針對選區(qū)激光熔化成形薄壁件過程中存在的變形較大、精度低等問題，通過獲得最優(yōu)工藝參數(shù)區(qū)間來減小薄壁件的變形。利用有限元軟件分析薄壁件成形過程中溫度場和應(yīng)力場的演化規(guī)律；建立形變量預(yù)測模型并進(jìn)行試驗(yàn)驗(yàn)證，研究工藝參數(shù)對薄壁件尺寸偏差的影響，得到激光功率、掃描速度與形變量之間的關(guān)系，實(shí)現(xiàn)對形變量的預(yù)測和控制。隨著掃描層數(shù)的增加，熔池的最高溫度和熱影響區(qū)也隨之增大，等溫線越密集，溫度梯度越大，最終趨于穩(wěn)定；薄壁件成形過程中，出現(xiàn)兩側(cè)壁邊緣向內(nèi)傾斜、上側(cè)邊緣出現(xiàn)內(nèi)凹的現(xiàn)象，薄壁件的最大應(yīng)力隨層數(shù)的增加而減小，最大熱應(yīng)力主要分布在薄壁件底層的兩端；形變量隨激光功率的增大而增大，隨掃描速度的增大而減小，薄壁件的形變量最小約為0.02 mm；試驗(yàn)驗(yàn)證所建立的數(shù)學(xué)模型誤差在10%左右，誤差較小，可以對形變量進(jìn)行良好的預(yù)測和控制。激光功率100～200 W、掃描速度800～1 000 mm/s為最優(yōu)參數(shù)區(qū)間；降低能量密度可以有效降低薄壁件形變量，提高其精度。

選區(qū)激光熔化；仿真分析；薄壁件；形變；數(shù)學(xué)模型

選區(qū)激光熔化（Selective Laser Melting，SLM）技術(shù)是金屬增材制造技術(shù)的一種，采用高能激光束對特定區(qū)域的金屬粉末進(jìn)行熔化后層層堆積成形。選區(qū)激光熔化技術(shù)理論上可以加工成形任意復(fù)雜的金屬零部件[1-2]，且具有加工時間短、成形精度高、成形件性能穩(wěn)定的優(yōu)點(diǎn)，在航空航天、醫(yī)療、軍工等方面擁有良好的應(yīng)用前景[3-4]。

由于SLM成形薄壁件困難，國內(nèi)外學(xué)者在相關(guān)技術(shù)領(lǐng)域進(jìn)行了大量研究[5]，陳帥等[6]研究了工藝參數(shù)對TiB2/AlSi10Mg密度與硬度的影響，利用不同參數(shù)獲得的SLM成形件優(yōu)于鑄造件。楊永強(qiáng)等[7]研究了工藝參數(shù)對SLM成形不銹鋼薄壁件的影響，得出了壁厚絕對誤差極限值為20 μm。陳俠宇等[8]通過對18NI300進(jìn)行多目標(biāo)的優(yōu)化，獲得了一組最佳的成形參數(shù)和優(yōu)異的成形件性能。在仿真研究中，Zhang等[9]通過ANSYS數(shù)值模擬和試驗(yàn)研究了不同跨度的AlSi10Mg拱橋結(jié)構(gòu)制造過程中的變形，跨度越大，拱彎曲變形越大，降低激光能量密度有助于減少拱變形。Li等[10]分析了在SLM過程中雙懸臂溫度場與殘余應(yīng)力場的演變，通過優(yōu)化使得懸臂上層的拉伸殘余應(yīng)力降低了70%。李保強(qiáng)等[11]利用熱-結(jié)構(gòu)耦合點(diǎn)曝光掃描成形應(yīng)力場，分析了點(diǎn)曝光時間和點(diǎn)間距對殘余應(yīng)力的影響。可知，在SLM成形薄壁件過程中，金屬粉末在熱源作用下極易產(chǎn)生較高的溫度梯度，導(dǎo)致較大的熱應(yīng)力，致使零件翹曲變形[12-13]。不同工藝參數(shù)決定著溫度梯度、熱應(yīng)力的變化，進(jìn)而影響零件的變形。因此，針對工藝參數(shù)對SLM成形無支撐薄壁件的多層熱力耦合行為及形變的影響還有待補(bǔ)充[14]。

文中通過對SLM成形薄壁件溫度場、應(yīng)力場的仿真研究，分析翹曲變形機(jī)理，探尋工藝參數(shù)對選區(qū)激光熔化成形薄壁件發(fā)生翹曲變形的影響機(jī)制，建立預(yù)測模型，獲得最佳工藝參數(shù)組合，為成形薄壁零件提供指導(dǎo)。

1 有限元模型

1.1 模型假設(shè)

采用ANSYS軟件建立薄壁零件的仿真模型，進(jìn)行有限元仿真，共定義基板、粉末、實(shí)體的3種材料類型，通過激光熱源對3種材料作用模擬整個成形過程。SLM成形薄壁件過程中，由于激光光斑直徑細(xì)小，作用在粉末表面的能量集中，同時熱源移動速度快，熔池形貌難以穩(wěn)定，周圍溫度場分布不均勻等原因，導(dǎo)致成形過程非平衡物理冶金過程和熱物理過程十分復(fù)雜[8]。為簡化有限元對SLM成型過程的模擬仿真，提出以下假設(shè)：

1）假設(shè)粉層為連續(xù)均勻介質(zhì)，不考慮粉末之間的傳熱。

2）不考慮熔池內(nèi)部的蒸發(fā)和熱毛細(xì)力作用。

3）假設(shè)形變遵循Von Mises屈服準(zhǔn)則、塑性流動性準(zhǔn)則和等向強(qiáng)化準(zhǔn)則。

4）當(dāng)時間間隔較小時，認(rèn)為材料的熱力學(xué)性能、熱應(yīng)力應(yīng)變線性變化。

1.2 幾何模型與網(wǎng)格劃分

模型進(jìn)行網(wǎng)格劃分時，應(yīng)當(dāng)綜合考慮模擬精度要求與計(jì)算效率，整體分為基板和成形件兩個部分?？拷勰┑牟糠志W(wǎng)格密集，采用映射的網(wǎng)格劃分方式，尺寸為0.02 mm；遠(yuǎn)離粉末的部分網(wǎng)格稀疏，選用網(wǎng)格尺寸為0.05 mm×0.05 mm×0.05 mm，采用solid90號單元，成形件總共分為4層，每層網(wǎng)格尺寸為0.01 mm×0.01 mm×0.03 mm，計(jì)算溫度場時，采用solid70號單元，計(jì)算應(yīng)力場時則轉(zhuǎn)變單元類型為solid45。為便于觀察溫度場變化情況，將模型從中點(diǎn)劃分，取其中一側(cè)進(jìn)行觀察，如圖1所示。粉末以單道的方式向上堆積成形，激光掃描方向從左往右，總共4層，箭頭表示激光掃描方向，同時在每層的中點(diǎn)選取4個觀測點(diǎn)。SLM成形薄壁件工藝參數(shù)設(shè)置如表1所示。

圖1 三維有限元模型

表1 有限元參數(shù)設(shè)置

Tab.1 Finite element parameter setting

1.3 控制方程

1.3.1 溫度場有限元模型

SLM成形薄壁件是指金屬粉末受激光熱源的影響局部受熱重熔的過程，粉末床表面溫度由于熱源的移動而不斷改變，因此，成形過程是一個典型的非線性瞬態(tài)熱傳導(dǎo)過程，其內(nèi)部的熱量變化可用式（1）表示[15]。

式中：為材料密度，g/cm3；為材料比熱容，J/(kg·K)；為粉床系統(tǒng)的溫度，℃；為激光作用時間，s；為內(nèi)熱源，W/m3；為粉末材料的導(dǎo)熱系數(shù)。

在SLM中，粉末床是通過吸收激光能量來熔化的。在這一過程中，材料內(nèi)部傳導(dǎo)的熱能約為熱對流和熱輻射[16-18]。

1.3.2 熱源模型

在有限元仿真成形過程中，激光熱源對溫度場和應(yīng)力場的計(jì)算精度有很大影響?？紤]到激光熔化金屬粉末過程中激光熱源存在一定的穿透效果，因此，熱源模型選取Gauss體熱源模型，熱源高度方向上指數(shù)衰減，激光光學(xué)穿透深度（）設(shè)定為鋪粉厚度，激光能量以熱流密度的形式施加在單元節(jié)點(diǎn)上。Gauss分布移動體熱源的數(shù)學(xué)表達(dá)式見式（3）[19]。

式中：為激光功率，W；為粉體材料對激光的吸收率，取0.3[12]；0為光斑半徑，mm；為激光掃描速度，mm/s。

1.3.3 應(yīng)力場有限元模型

通過模型計(jì)算應(yīng)力時，采用熱-結(jié)構(gòu)耦合（間接耦合）的方式進(jìn)行研究。對基板進(jìn)行約束，將計(jì)算好的溫度場作為前提施加給對應(yīng)位置以完成應(yīng)力場的計(jì)算，其中，熱-結(jié)構(gòu)耦合分析流程示意圖如圖2所示。粉末在激光熱源作用下熔化的區(qū)域存在兩種變形，即彈性變形與塑性變形。彈性變形由于溫差過大致使材料受熱膨脹而產(chǎn)生；塑性變形則由于熱應(yīng)力值大于材料屈服極限而產(chǎn)生。在整個模擬過程中，材料應(yīng)力應(yīng)變行為遵循熱-彈-塑性原理。在塑性變形過程中，應(yīng)力應(yīng)變行為應(yīng)遵循流動準(zhǔn)則和等向強(qiáng)化準(zhǔn)則[9]。

圖2 熱-應(yīng)力耦合分析流程

1.4 材料物性

在對SLM加工過程進(jìn)行模擬分析時，著重考慮不同溫度下材料的熱物性參數(shù)和力學(xué)性能參數(shù)，二者在仿真模型計(jì)算精度方面有著重要的影響。在溫度場模擬過程中，熱物性參數(shù)為主要考慮對象，如熔點(diǎn)、密度、比熱容和熱導(dǎo)率等；在應(yīng)力場模擬過程中，力學(xué)性能參數(shù)為主要考慮點(diǎn)，如泊松比、線膨脹系數(shù)、屈服強(qiáng)度和彈性模量等，具體如圖3所示。該仿真試驗(yàn)中，基板材料選取45鋼，加工粉末選取316L不銹鋼粉末，其固相線為1 399 ℃，液相線為1 450 ℃，泊松比為0.3。

2 仿真結(jié)果與分析

2.1 溫度場分析

溫度場分析由冷卻速度和成形過程中熔池的尺寸與熱影響區(qū)的變化展開。圖4為4個觀測點(diǎn)溫度隨時間的變化曲線。當(dāng)溫度曲線超過熔點(diǎn)的紅色虛線時，粉末為熔融態(tài)?？梢钥闯觯麄€曲線出現(xiàn)類似于阻尼的振蕩“鋸齒波”，除了最后一層外，前3層的溫度曲線均超過熔點(diǎn)兩次。下一層中點(diǎn)的第2個溫度峰值出現(xiàn)在上層中點(diǎn)第1個峰值的時間點(diǎn)上，表明當(dāng)激光移動到上層的中點(diǎn)處，粉末快速熔化形成熔池時，多余的激光能量使得相鄰下層出現(xiàn)重熔的現(xiàn)象，下層中點(diǎn)出現(xiàn)第2個溫度峰值。隨著薄壁高度的增加，傳熱路徑增長，隨后幾次傳熱過程中下層的峰值逐次減少并趨于穩(wěn)定。從圖4可知，當(dāng)激光靠近又遠(yuǎn)離中點(diǎn)時，中點(diǎn)升溫速度明顯快于降溫速度，符合選區(qū)激光熔化技術(shù)快速熔化和快速冷卻的特性。由于熱量的積累和預(yù)熱效應(yīng)，4段曲線的最高溫度呈上升趨勢，其峰值對應(yīng)每個掃描矢量末端的激光位置，說明激光經(jīng)過每層的最高溫度呈升高趨勢。熔池的最高溫度分別為2 064、2 192、2 299、2 396 ℃。由于選區(qū)激光熔化快速冷卻的特性，熔池的溫度并非隨著層數(shù)的增加而無限制地升高，而是隨著層數(shù)的增加逐漸趨于穩(wěn)定。

圖3 316L熱物性和力學(xué)性能參數(shù)

圖4 溫度隨時間的變化曲線

圖5為4層中心的溫度場。溫度云圖中不同顏色代表著不同的溫度范圍，其中深紅色代表1 435 ℃以上的溫度區(qū)域，而該溫度區(qū)域?yàn)?16L不銹鋼的液相溫度[20]，因此，該區(qū)域?yàn)槿鄢販囟取D5顯示了從第1層到第4層溫度場云圖分布情況，可以看出，每層的剖面熔池形狀類似于“彗星”[21]，后面熱影響區(qū)相當(dāng)于彗星的尾巴。后面熱影響區(qū)范圍遠(yuǎn)大于前面熱影響區(qū)，但是等溫線密集程度卻小于熔池前面熱影響區(qū)。產(chǎn)生這一現(xiàn)象的原因是材料導(dǎo)熱系數(shù)不同，熔池前面等溫線密集是因?yàn)槿鄢厍岸螢槲赐耆刍姆勰?，而后面是?jīng)熔化后凝固的固體狀態(tài)，隨著各層參數(shù)的增大，層與層之間的熱積累效應(yīng)使得熔池和熱影響區(qū)有緩慢變大的趨勢。第1層基板與加工料的熱導(dǎo)率不同，使得粉末熔化過程中部分熱量被基板傳遞出去，當(dāng)激光掃描第2層時，第1層殘留的余熱對第2層起到預(yù)熱效果，同時上一層受熱熔化的能量也傳遞給了下一層，使得下層產(chǎn)生了重熔現(xiàn)象，因此，層與層之間產(chǎn)生良好的冶金結(jié)合，最后得到的成形樣件有著良好的致密性和機(jī)械性能。熔池的尺寸影響成形件的尺寸精度和力學(xué)性能，熔道的搭接最終影響成形件的整體。

圖5 每層的熱影響區(qū)

2.2 應(yīng)力分析

為了探究薄壁零件應(yīng)力場的演變過程，通過4個觀測點(diǎn)獲得激光功率為100 W時等效應(yīng)力隨時間的變化曲線，如圖6所示。從圖6中可以看出，各層中點(diǎn)的等效應(yīng)力循環(huán)曲線和溫度曲線相似，經(jīng)歷多個波峰波谷。當(dāng)激光熱源經(jīng)過第1層中點(diǎn)時，粉末熔化形成熔池，此時液態(tài)熔池里熱應(yīng)力為零，隨著激光的快速離去，熔池冷卻，金屬凝固，從而產(chǎn)生應(yīng)力。等到金屬完全凝固，與基板接觸部分產(chǎn)生相互作用，應(yīng)力迅速升高至250 MPa。隨后移動熱源分別經(jīng)過第3、4層中點(diǎn)時，第1層中點(diǎn)經(jīng)歷了3次溫度回升，但最高應(yīng)力值卻在下降，這是因?yàn)楫?dāng)激光經(jīng)過第2層中點(diǎn)時，第1層中點(diǎn)處出現(xiàn)重熔的現(xiàn)象[22]，使得這一層的熱應(yīng)力得到釋放。同時，第2、3、4層對第1層也有著應(yīng)力釋放的作用[23]。由圖6可以看出，第1層產(chǎn)生的應(yīng)力最大，因?yàn)榧す鈱Ψ勰┻M(jìn)行熔化時，基板與粉末層的熱導(dǎo)率和熱膨脹系數(shù)不一樣，同時，由于選區(qū)激光熔化技術(shù)速熱速冷的特點(diǎn)，從而產(chǎn)生巨大的溫度梯度，使得第1層應(yīng)力峰值達(dá)到最大。第3、4層波峰數(shù)較第1、2層逐次減少，表明前幾層殘留的余溫對后幾層起到保溫作用，溫度梯度隨之降低，從而產(chǎn)生的最大應(yīng)力隨之減小，最終穩(wěn)定為150 MPa。在成形薄壁零件時會出現(xiàn)零件與基板脫離或斷裂的現(xiàn)象，尤其出現(xiàn)最底部的兩端，圖6虛框內(nèi)為第1層沿路徑方向等效應(yīng)力分布的幾何形狀，可以明顯看出，最大的拉應(yīng)力出現(xiàn)在零件的兩端，約為300 MPa。其中，右端應(yīng)力最大的原因是該位置為激光最后掃描的位置，激光停留在此處造成大的溫度梯度和熱應(yīng)力，從而產(chǎn)生大的形變。當(dāng)整體打印完成，薄壁零件會出現(xiàn)兩側(cè)壁的邊緣向內(nèi)傾斜，上側(cè)邊緣出現(xiàn)內(nèi)凹的現(xiàn)象。

圖6 各層應(yīng)力曲線

2.3 工藝參數(shù)對形變量的影響

對不同工藝參數(shù)成形的薄壁零件進(jìn)行仿真分析，獲得成形件形變數(shù)值，如表2所示。所獲得數(shù)據(jù)采用MATLAB進(jìn)行擬合，分析激光功率和掃描速度對形變量的影響，獲得工藝參數(shù)對形變量影響的關(guān)系圖，如圖7所示。其中，綠色實(shí)線代表影響趨勢，形變量隨著激光功率的增大而增大，隨著掃描速度的增大而減小。圖8為工藝參數(shù)對形變量影響的三維關(guān)系圖。如圖8a所示，在不同功率下，掃描速度不同，對薄壁件變形影響不同。當(dāng)激光功率為150 W時，掃描速度為400～600 mm/s，薄壁件的變形變化趨勢較為平緩；掃描速度為600～1 000 mm/s，薄壁件的變形量呈線性遞減。而當(dāng)激光功率為300 W時，掃描速度為400～800 mm/s，薄壁件的變形較大，呈線性增長；掃描速度控制在800～1 000 mm/s時，薄壁件的變形受掃描速度的影響較小。結(jié)合圖8b可知，掃描速度在較低的激光功率下對薄壁件變形的影響小于功率較大時。因?yàn)楫?dāng)激光功率過高時，掃描熔道吸收的能量過多，使得熔池溫度升高，導(dǎo)致黏度增大出現(xiàn)球化，以及熔池內(nèi)部的炸裂產(chǎn)生孔洞，造成嚴(yán)重形變。隨著掃描速度的增大或激光功率的減小，能量密度隨之減小，可以有效地降低上述現(xiàn)象的發(fā)生幾率，所以變形較小。結(jié)合圖8的分析，為獲得變形量較小的薄壁件，當(dāng)激光功率選擇100～200 W、掃描速度在800～ 1 000 mm/s時為最優(yōu)區(qū)間，所得到薄壁零件的形變量最小約為0.02 mm。

表2 仿真參數(shù)與結(jié)果

Tab.2 Simulation parameters and results

2.4 試驗(yàn)驗(yàn)證

為了驗(yàn)證仿真的準(zhǔn)確性，根據(jù)表2的工藝參數(shù)進(jìn)行試驗(yàn)成形，使用德國SLM-solutions公司開發(fā)的SLM-125打印機(jī)打印成形，采用商業(yè)化316L不銹鋼粉末，平均粒徑為15～45 μm，試驗(yàn)件是仿真模型尺寸的15倍，其試驗(yàn)設(shè)備與粉末形貌如圖9所示。

2.4.1 試驗(yàn)結(jié)果分析

首先對仿真模型進(jìn)行成形試驗(yàn)，圖10顯示了單道4層樣件的截面形貌，從圖中可以看出，隨著層數(shù)的增加，熔化獲得的截面寬度增加，球化、掛渣等現(xiàn)象減少。這是因?yàn)樵诳焖偃刍湍踢^程中，通過傳導(dǎo)的熱損失高于對流輻射[24]。在單道多層的模型中，成形樣件相對較小，周圍粉末具有較低的熱導(dǎo)率，因此熱量主要通過底層基板傳遞。然而，隨著層數(shù)的增加，熱傳導(dǎo)路徑增長，熔池逐漸增大。在實(shí)際成形過程中，隨著樣件的逐層堆積，熔池范圍逐漸趨于穩(wěn)定，與仿真的結(jié)果相一致。

圖7 不同參數(shù)下的趨勢圖

圖8 參數(shù)對形變量影響的三維圖和等高線圖

圖9 試驗(yàn)設(shè)備及316L不銹鋼粉末SEM圖像

圖10 單道4層樣件截面形貌

通過打印成形獲得的薄壁樣件如圖11所示，圖11a從左往右分別代表的是不同的掃描速度，從上往下代表著不同的激光功率，其中，線框圈出的部分分別表示激光功率為250 W(1)、300 W(2)，掃描速度為400 mm/s時的成形樣件，相對于其他樣件，其顏色呈暗黑色，高能量密度使得成形過程中產(chǎn)生過燒現(xiàn)象，造成較大的變形，同時，圖11b中每個薄壁零部件底部與基板相接觸的地方均出現(xiàn)黑色區(qū)域，由于基板與成形件的材質(zhì)不一樣而使得導(dǎo)熱性和膨脹系數(shù)不一樣，導(dǎo)致在初始成形時，成形件與基板的接觸端底部出現(xiàn)較大的溫度梯度和熱應(yīng)力，因此，在成形薄壁零件時，底部常出現(xiàn)翹曲和分離等現(xiàn)象。在試驗(yàn)前期，可通過采用基板高溫預(yù)熱減小溫度梯度，從而避免此現(xiàn)象的發(fā)生。

圖11 打印成形獲得的薄壁樣件

2.4.2 數(shù)學(xué)模型的建立與驗(yàn)證

采用三坐標(biāo)測量儀對成形薄壁的寬度進(jìn)行測量，薄壁件的兩側(cè)從上到下總共測量12個點(diǎn)，測量5次以減少測量誤差，以平面度作為變形。獲得的樣件測量結(jié)果如圖12所示，采用最大絕對誤差值來評判最大的形變量值。為了量化上述結(jié)果，根據(jù)試驗(yàn)數(shù)據(jù)建立激光功率、掃描速度對形變量的一階回歸方程，如式（4）所示。

擬合得到的公式如式（5）所示，該方程的相關(guān)系數(shù)=0.86，根據(jù)系數(shù)大小可知，激光功率對形變的影響大于掃描速度的影響。

為了更加直觀地反映線性關(guān)系，對該方程進(jìn)行改進(jìn)和轉(zhuǎn)化，如式（6）所示。

式中：1、2、3為利用數(shù)據(jù)得到的模型參數(shù)。在本研究中，運(yùn)用線性假設(shè)對不同工藝下的形變量參數(shù)進(jìn)行預(yù)測。

得到相關(guān)系數(shù)R=0.944，發(fā)現(xiàn)式（7）比式（4）擬合程度更好，更為準(zhǔn)確。其中，形變量與預(yù)測值如圖13所示，所得到的點(diǎn)在斜率約為0.5的直線上下浮動，證明所得模型擬合良好。根據(jù)圖12實(shí)際試驗(yàn)測量的相對誤差建立工藝參數(shù)與薄壁形變量的關(guān)系，如圖14所示，所得趨勢圖與仿真結(jié)果一致。結(jié)合式（7）可得出激光功率對形變量有顯著的正效應(yīng)，而掃描速度對形變量有負(fù)效應(yīng)，通過合理選擇工藝參數(shù)能夠?qū)Ρ”诩淖冃芜M(jìn)行有效控制。

圖13 改進(jìn)后的擬合曲線

圖14 不同工藝參數(shù)下的形變趨勢圖

對上述數(shù)學(xué)模型進(jìn)行驗(yàn)證，取不同于表2的工藝參數(shù)進(jìn)行成形試驗(yàn)，測得的結(jié)果如表3所示，相對誤差約為12.2%、9.59%，與Zhang等[25]建立的數(shù)學(xué)模型相比，其最優(yōu)工藝參數(shù)下的預(yù)測值與實(shí)際值之差相近，驗(yàn)證了方程的準(zhǔn)確性。同時，通過對圖14中工藝參數(shù)對形變量影響的分析可知，激光功率在100～200 W，掃描速度在800～1 000 mm/s范圍內(nèi)的形變量最小，成形最優(yōu)。

表3 數(shù)學(xué)模型驗(yàn)證

Tab.3 Mathematical model verification

3 結(jié)論

通過有限元仿真分析SLM成形薄壁零件過程中熔池的熱行為與應(yīng)力的演變，建立了激光功率和掃描速度與成形件形變量的關(guān)系模型并通過試驗(yàn)驗(yàn)證，得出以下主要結(jié)論。

1）隨著層數(shù)的增加，熔池的熱影響區(qū)大小有緩慢增大的趨勢；熔池溫度越高，等溫線越密集，溫度梯度越大，最終趨于穩(wěn)定；熱應(yīng)力從頂層向基板逐漸增大，最大熱應(yīng)力主要分布在薄壁零件底層的兩端。

2）形變量隨著激光功率的增大而增大，隨著掃描速度的增大而減??；降低能量密度可以有效降低形變量。

3）建立SLM成形薄壁件形變量模型能夠有效地預(yù)測變量；形變量最小的工藝參數(shù)區(qū)間是激光功率為100～200 W，掃描速度為800～1 000 mm/s。

[1] LI Ying-li, ZHOU Kun, TAN Peng-fei, et al. Modeling Temperature and Residual Stress Fields in Selective Laser Melting[J]. International Journal of Mechanical Sciences, 2018, 136: 24-35.

[2] 曾壽金, 吳啟銳, 葉建華. 選區(qū)激光熔化成型316L不銹鋼多孔結(jié)構(gòu)的力學(xué)性能[J]. 紅外與激光工程, 2020, 49(8): 67-75.

ZENG Shou-jin, WU Qi-rui, YE Jian-hua. Mechanical Properties of 316L Stainless Steel Porous Structure Formed by Selective Laser Melting[J]. Infrared and Laser Engineering, 2020, 49(8): 67-75.

[3] 陳旭斌, 葛翔, 祝毅, 等. 選擇性激光熔化零件微觀結(jié)構(gòu)及摩擦學(xué)性能研究[J]. 機(jī)械工程學(xué)報, 2018, 54(3): 63-72.

CHEN Xu-bin, GE Xiang, ZHU Yi, et al. A Study on Microstructure and Tribology Performance of Samples Processed by Selective Laser Melting (SLM)[J]. Journal of Mechanical Engineering, 2018, 54(3): 63-72.

[4] 楊志成. 鈦合金薄壁結(jié)構(gòu)選區(qū)激光熔化成形機(jī)理及力學(xué)性能研究[D]. 太原: 中北大學(xué), 2022.

YANG Zhi-cheng. Study on Selective Laser Melting Forming Mechanism and Mechanical Properties of Titanium Alloy Thin-Walled Structure[D]. Taiyuan: North University of China, 2022.

[5] 陳帥, 劉建光, 王衛(wèi)東, 等. 激光選區(qū)熔化成形薄壁件研究進(jìn)展[J]. 精密成形工程, 2020, 12(5): 122-131.

CHEN Shuai, LIU Jian-guang, WANG Wei-dong, et al. Research Progress in Thin-Walled Parts Formed by Selective Laser Melting[J]. Journal of Netshape Forming Engineering, 2020, 12(5): 122-131.

[6] 陳帥, 劉建光, 王衛(wèi)東, 等. TiB2/AlSi10Mg激光選區(qū)熔化成形工藝研究[J]. 精密成形工程, 2021, 13(3): 154-161.

CHEN Shuai, LIU Jian-guang, WANG Wei-dong, et al. Research on the Selective Laser Melting Process of TiB2/AlSi10Mg[J]. Journal of Netshape Forming Engineering, 2021, 13(3): 154-161.

[7] 楊永強(qiáng), 羅子藝, 蘇旭彬, 等. 不銹鋼薄壁零件選區(qū)激光熔化制造及影響因素研究[J]. 中國激光, 2011, 38(1): 60-67.

YANG Yong-qiang, LUO Zi-yi, SU Xu-bin, et al. Study on Process and Effective Factors of Stainless Steel Thin-Wall Parts Manufactured by Selective Laser Melting[J]. Chinese Journal of Lasers, 2011, 38(1): 60-67.

[8] 陳俠宇, 黃衛(wèi)東, 張偉杰, 等. 基于灰色關(guān)聯(lián)分析的選區(qū)激光熔化成形18Ni300模具鋼多目標(biāo)工藝優(yōu)化[J]. 中國激光, 2020, 47(5): 341-351.

CHEN Xia-yu, HUANG Wei-dong, ZHANG Wei-jie, et al. Multiple Targets Technology Optimization Based Grey Relative Analysis of 18Ni300 Die Steel Formed by Selective Laser Melting[J]. Chinese Journal of Lasers, 2020, 47(5): 341-351.

[9] ZHANG Xiao-chuan, WANG Ji-wu, KANG Jin-wu, et al. The Dynamic Arch Bending Mechanism of Flat Bridge Structure of AlSi10Mg during SLM Process[J]. Materials & Design, 2020, 188: 108469.

[10] LI C, LIU J F, FANG X Y, et al. Efficient Predictive Model of Part Distortion and Residual Stress in Selective Laser Melting[J]. Additive Manufacturing, 2017, 17: 157-168.

[11] 李保強(qiáng), 李忠華, 白培康, 等. 選區(qū)激光熔化AlSi10Mg應(yīng)力場數(shù)值模擬研究[J]. 應(yīng)用激光, 2019, 39(2): 211-216.

LI Bao-qiang, LI Zhong-hua, BAI Pei-kang, et al. Numerical Simulation of Stress Field for AlSi10Mg Fabricated by Selective Laser Melting[J]. Applied Laser, 2019, 39(2): 211-216.

[12] HUSSEIN A, HAO Liang, YAN Chun-ze, et al. Finite Element Simulation of the Temperature and Stress Fields in Single Layers Built Without-Support in Selective Laser Melting[J]. Materials & Design (1980-2015), 2013, 52: 638-647.

[13] 姚耀伍. SLM成型AlSi10Mg傾斜薄壁件變形研究[D]. 福州: 福建工程學(xué)院, 2022.

YAO Yao-wu. Deformation of SLM-formed AlSi10Mg Inclined Thin-Walled Parts[D]. Fuzhou: Fujian University of Technology, 2022.

[14] 柯林達(dá), 殷杰, 朱海紅, 等. 鈦合金薄壁件選區(qū)激光熔化應(yīng)力演變的數(shù)值模擬[J]. 金屬學(xué)報, 2020, 56(3): 374-384.

KE Lin-da, YIN Jie, ZHU Hai-hong, et al. Numerical Simulation of Stress Evolution of Thin-Wall Titanium Parts Fabricated by Selective Laser Melting[J]. Acta Metallurgica Sinica, 2020, 56(3): 374-384.

[15] 李雅莉. 選區(qū)激光熔化AlSi10Mg溫度場及應(yīng)力場數(shù)值模擬研究[D]. 南京: 南京航空航天大學(xué), 2015.

LI Ya-li. Numerical Investigation on Temperature Field and Stress Field during Selective Laser Melting of AlS10Mg[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2015.

[16] STREEK A, REGENFUSS P, EXNER H. Fundamentals of Energy Conversion and Dissipation in Powder Layers during Laser Micro Sintering[J]. Physics Procedia, 2013, 41: 858-869.

[17] HUANG Wei-bo, ZHANG Yi-min. Finite Element Simulation of Thermal Behavior in Single-Track Multiple-Layers Thin Wall Without-Support during Selective Laser Melting[J]. Journal of Manufacturing Processes, 2019, 42: 139-148.

[18] MISHRA A K, KUMAR A. Numerical and Experimental Analysis of the Effect of Volumetric Energy Absorption in Powder Layer on Thermal-Fluidic Transport in Selective Laser Melting of Ti6Al4V[J]. Optics & Laser Technology, 2019, 111: 227-239.

[19] 艾百運(yùn), 姜勇, 柯林達(dá), 等. 基于ANSYS的SLM成形熱彈塑性有限元分析[J]. 航天制造技術(shù), 2016(3): 10-15.

AI Bai-yun, JIANG Yong, KE Lin-da, et al. Thermal Elastoplastic FEM Analysis in SLM Process Based on ANSYS[J]. Aerospace Manufacturing Technology, 2016(3): 10-15.

[20] 段成紅, 郝曉杰, 羅翔鵬. 選區(qū)激光熔化316L溫度場研究[J]. 應(yīng)用激光, 2018, 38(5): 748-753.

DUAN Cheng-hong, HAO Xiao-jie, LUO Xiang-peng. Study on Temperature Field of Selective Laser Melting 316L[J]. Applied Laser, 2018, 38(5): 748-753.

[21] DEMIR A G, MAZZOLENI L, CAPRIO L, et al. Complementary Use of Pulsed and Continuous Wave Emission Modes to Stabilize Melt Pool Geometry in Laser Powder Bed Fusion[J]. Optics & Laser Technology, 2019, 113: 15-26.

[22] LIU Shi-wen, ZHU Hai-hong, PENG Gang-yong, et al. Microstructure Prediction of Selective Laser Melting AlSi10Mg Using Finite Element Analysis[J]. Materials & Design, 2018, 142: 319-328.

[23] 丁林山. 基于瞬態(tài)有限元分析的選擇性激光熔化多層掃描路徑優(yōu)化[D]. 武漢: 華中科技大學(xué), 2018.

DING Lin-shan. Optimization of Multi-Layer Scanning Path for Selective Laser Melting Based on Transient Finite Elements Analysis[D]. Wuhan: Huazhong University of Science and Technology, 2018.

[24] ?IRIN T B, KAYNAK Y. Prediction of Residual Stress and Distortion in Laser Powder Bed Fusion Additive Manufacturing Process of Inconel 718 Alloy[J]. Procedia CIRP, 2021, 99: 330-335.

[25] ZHANG Luo, ZHANG Sha-sha, ZHU Hai-hong, et al. Horizontal Dimensional Accuracy Prediction of Selective Laser Melting[J]. Materials & Design, 2018, 160: 9-20.

Size Deviation Prediction and Control of Thin-walled Parts Formed by SLM

XU Ming-sana,b, JIANG Yao-fenga,b, YAO Yao-wua,b, CHEN Xiang-danga,b, ZENG Shou-jina,b, YE Jian-huaa,b

(a. Fujian Key Laboratory of Intelligent Machining Technology and Equipment, b. School of Mechanical and Automotive Engineering, Fujian University of Technology, Fuzhou 350118, China)

The work aims to reduce the deformation of thin-walled parts by obtaining the optimal process parameter range, so as to solve the problems of large deformation and low precision in the process of selective laser melting forming. The evolution law of temperature field and stress field during the forming process of thin-walled parts was analyzed by finite element software. The shape variable prediction model was established and verified by experiments. Then, the effect of process parameters on the size deviation of thin-walled parts was studied and the relationship between laser power, scanning speed and shape variable was obtained to realize the prediction and control of shape variable. With the increase of the number of scanning layers, the maximum temperature and heat affected zone of the molten pool also increased. The denser the isotherms were, the greater the temperature gradient was, which finally tended to be stable. During the forming process of thin-walled parts, the edges of both sides of the wall inclined inward, and the upper edge was concave. The maximum stress of thin-walled parts decreased with the increase of the number of layers, and the maximum thermal stress was mainly distributed at both ends of the bottom layer of thin-walled parts. The shape variable increased with the increase of laser power and decreased with the increase of scanning speed. The minimum shape variable of thin-walled parts was about 0.02 mm. The experiments verified that the error of the established mathematical model was minor, about 10%, which could predict and control the shape variables well. The optimal parameter range covers the laser power of 100-200 W and the scanning speed of 800-1 000 mm/s. Reducing the energy density can effectively reduce the shape variable of thin-walled parts and improve the precision.

selective laser melting; simulation analysis; thin-walled parts; shape variable; mathematical model

10.3969/j.issn.1674-6457.2023.02.011

TN249

A

1674-6457(2023)02-0086-09

2022?04?22

2022-04-22

國家自然科學(xué)基金（51575110）；福建省自然科學(xué)基金（2020J01827）

National Natural Science Foundation of China (51575110); Natural Science Foundation of Fujian Province (2020J01827)

許明三（1974—），男，碩士，教授，主要研究方向?yàn)榧す庠霾闹圃?、硬脆材料加工技術(shù)及智能制造裝備設(shè)計(jì)。

XU Ming-san (1974-), Male, Master, Professor, Research focus: laser additive manufacturing, hard and brittle material processing technology, and intelligent manufacturing equipment design.

許明三, 江堯峰, 姚耀伍, 等. SLM成形薄壁件尺寸偏差預(yù)測與控制研究[J]. 精密成形工程, 2023, 15(2): 86-94.

XU Ming-san, JIANG Yao-feng, YAO Yao-wu, et al. Size Deviation Prediction and Control of Thin-walled Parts Formed by SLM[J]. Journal of Netshape Forming Engineering, 2023, 15(2): 86-94.