梁平華 唐 倩 余志強(qiáng) 馮琪翔

劉 威11.重慶大學(xué)機(jī)械傳動(dòng)國家重點(diǎn)實(shí)驗(yàn)室，重慶,4000442.成都飛機(jī)工業(yè)(集團(tuán))有限責(zé)任公司，成都,610092

0 引言

金屬增材制造技術(shù)可分為同步送粉/送絲成形技術(shù)和預(yù)鋪粉床成形技術(shù)兩種，預(yù)鋪粉床成形技術(shù)可分為選擇性激光熔化(selective laser melting, SLM) 和電子束選區(qū)熔化(electron beam melting, EBM)[1]，兩者都是將高能激光束或電子束通過激光振鏡或偏轉(zhuǎn)線圈投射到金屬粉床平面上，使局部金屬粉末迅速熔化而成形。在加工前建立零件的三維結(jié)構(gòu)，導(dǎo)入STL文件后，根據(jù)材料的屬性和打印機(jī)的規(guī)格對(duì)其進(jìn)行分層切片，并添加必要的支撐結(jié)構(gòu)[2]。激光熔化金屬粉末時(shí)間很短，許多基于激光的加工過程都伴隨有飛濺現(xiàn)象[3]，形成大量黑煙。為有效去除加工過程中產(chǎn)生的黑煙，同時(shí)也為了防止金屬高溫時(shí)被氧化，必須在加工面上形成有效流動(dòng)的惰性保護(hù)氣風(fēng)場[3]。FERRAR等[4]研究發(fā)現(xiàn)，在成形缸基板面上的惰性氣體風(fēng)場分布會(huì)對(duì)成形試件力學(xué)性能和可重復(fù)性產(chǎn)生顯著影響；李金梁等[5]研究了金屬選區(qū)燒結(jié)過程中氣氛保護(hù)系統(tǒng)的原理構(gòu)造；ANWAR等[6]研究了AlSi10Mg合金在SLM成形過程中惰性保護(hù)氣流速度對(duì)成形質(zhì)量的影響，發(fā)現(xiàn)相同風(fēng)速下激光掃描方向與風(fēng)場方向相反時(shí)能更有效地去除黑煙雜質(zhì)；宋長輝等[7]研究了激光掃描間距和風(fēng)場對(duì)CoCrMo合金成形質(zhì)量的綜合影響。然而，對(duì)SLM設(shè)備風(fēng)場的相關(guān)研究大都集中在風(fēng)場風(fēng)速、激光相對(duì)風(fēng)向掃描策略等對(duì)成形質(zhì)量的影響，較少有運(yùn)用CFD流體仿真技術(shù)對(duì)形成風(fēng)場的流道進(jìn)行優(yōu)化改進(jìn)以實(shí)現(xiàn)內(nèi)部風(fēng)場布局優(yōu)化。

本文基于流體仿真技術(shù)對(duì)激光選區(qū)熔化過程中的風(fēng)場布局進(jìn)行仿真，并對(duì)保護(hù)氣流道結(jié)構(gòu)進(jìn)行了多次改進(jìn)優(yōu)化和仿真檢驗(yàn)。仿真結(jié)果表明，優(yōu)化后的保護(hù)氣流道結(jié)構(gòu)能更有效地去除黑煙。

1 SLM設(shè)備成形系統(tǒng)組成

1.回風(fēng)軌道 2.成形缸 3.下進(jìn)風(fēng)口 4.激光振鏡 5.上進(jìn)風(fēng)口 6.料缸圖1 成形倉Fig.1 Forming chamber

1.與循環(huán)風(fēng)機(jī)相連的總進(jìn)風(fēng)口 2.上進(jìn)風(fēng)口 3.下進(jìn)風(fēng)口圖2 循環(huán)保護(hù)氣分支流道Fig.2 Branch runner of the circulating protective gas

以E-Plus-M250激光選區(qū)熔化設(shè)備為研究對(duì)象，測量提取得到其成形倉和循環(huán)保護(hù)氣分支流道模型，如圖1和圖2所示。成形倉整體為一個(gè)長方體空間，基本尺寸為110 mm×42 mm×48 mm，主要包括回風(fēng)軌道、成形缸、下進(jìn)風(fēng)口、上進(jìn)風(fēng)口、激光振鏡、料缸、鋪粉車等基本結(jié)構(gòu)。其中,上進(jìn)風(fēng)口布置在激光振鏡正下方，主要是為了避免減小加工過程中黑煙吸附在激光振鏡上而影響激光正常入射。下進(jìn)風(fēng)口布置在成形缸加工粉床平面上方19 mm處，主要功能是去除加工過程中產(chǎn)生的黑煙雜質(zhì)。循環(huán)保護(hù)氣分支流道由一個(gè)總進(jìn)風(fēng)口和兩個(gè)出風(fēng)口(成形倉上下進(jìn)風(fēng)口)構(gòu)成，作用是將來自循環(huán)風(fēng)機(jī)的保護(hù)氣分成兩路,分別送入成形倉對(duì)應(yīng)位置。

E-Plus-M250設(shè)備的保護(hù)氣循環(huán)系統(tǒng)由一臺(tái)帶過濾功能的循環(huán)風(fēng)機(jī)提供動(dòng)力，使成形倉在工作時(shí)保持保護(hù)氣(氮?dú)?循環(huán)流動(dòng)，從而去除激光掃描粉床面時(shí)產(chǎn)生的黑煙，同時(shí)帶走一部分激光熔化金屬粉末而產(chǎn)生的熱量。

2 SLM成形倉風(fēng)場仿真

經(jīng)由循環(huán)風(fēng)機(jī)過濾后的氮?dú)鈴某尚蝹}背后的與循環(huán)風(fēng)機(jī)相連的總進(jìn)風(fēng)口進(jìn)入，一部分進(jìn)入上進(jìn)風(fēng)口，一部分進(jìn)入下進(jìn)風(fēng)口。然而，在正常加工過程中，激光掃射粉床面時(shí)產(chǎn)生的黑煙并不能有效地被去除，仍有大量黑煙雜質(zhì)直接掉落在粉床表面，這對(duì)加工零件的致密度、力學(xué)性能都有很大的影響；掉落在粉床上的黑煙雜質(zhì)會(huì)參混在金屬粉末中而不能去除，對(duì)加工金屬粉料的純凈度也會(huì)產(chǎn)生影響。通過查看激光掃描后粉床上黑煙殘留情況，發(fā)現(xiàn)成形缸中的黑煙殘留主要集中在右側(cè)部分，甚至發(fā)現(xiàn)當(dāng)激光掃略軌跡偏右時(shí)，在料缸粉床面上也有殘留黑煙，這說明在成形缸上方的循環(huán)保護(hù)氣速度分布不均勻。

2.1 流場仿真分析理論

由于仿真模擬未涉及溫度場的模擬，故在模擬時(shí)不考慮能量方程?；贜-S方程中的質(zhì)量守恒方程和動(dòng)量守恒方程，對(duì)E-Plus-M250金屬3D打印機(jī)成形倉內(nèi)的風(fēng)場分布情況進(jìn)行速度場模擬。選擇k-ε湍流模型分別對(duì)保護(hù)氣分支流道和成形倉整體內(nèi)部的風(fēng)場進(jìn)行迭代求解[8-9]。

質(zhì)量守恒控制方程如下：

(1)

式中，ρ為流體微元體密度；t為時(shí)間；u、v、w分別為微元體在x、y、z方向的速度分量。

該方程表示流場中單位時(shí)間流體微元體中質(zhì)量的增加都等于該時(shí)間段內(nèi)流入微元體的質(zhì)量。

動(dòng)量守恒控制方程(x方向)如下：

(2)

式中，u為微元體速度矢量；p為作用于微元體上的壓力，τxx、τyx、τzx分別為微元體所受的黏性切應(yīng)力分量；Fx為微元體所受體積力分量。

式(2)表示微元體流體動(dòng)量的變化等于作用在微元體上的合力。

控制方程通式為

(3)

式中，φ為通用變量；Г為廣義擴(kuò)散系數(shù)；S為廣義源。

2.2 建立分支流道計(jì)算模型

激光加工過程中不能有效去除黑煙雜質(zhì)，是成形倉中的循環(huán)保護(hù)氣進(jìn)風(fēng)、回風(fēng)口空間布局或結(jié)構(gòu)不合理造成的。為了探究成形缸加工粉床面上方流場具體分布情況，并根據(jù)流場分布對(duì)保護(hù)氣分支流道結(jié)構(gòu)進(jìn)行優(yōu)化，首先對(duì)循環(huán)保護(hù)氣分支流道進(jìn)行了計(jì)算流體力學(xué)仿真分析，得到了成形倉上進(jìn)風(fēng)口、下進(jìn)風(fēng)口截面的風(fēng)場分布情況。

提取E-Plus-M250循環(huán)保護(hù)氣分支流道結(jié)構(gòu)尺寸參數(shù)后，在ANSYS ICEM CFD中建立等效流道模型并進(jìn)行網(wǎng)格劃分，由于流道結(jié)構(gòu)形狀不規(guī)則，因此采用非結(jié)構(gòu)化網(wǎng)格建立流道網(wǎng)格模型，建立的流道網(wǎng)格模型見圖3。

圖3 循環(huán)保護(hù)氣分支流道網(wǎng)格模型Fig.3 Grid model of gas branche

將建立好的網(wǎng)格模型導(dǎo)入ANSYS FLUENT中進(jìn)行計(jì)算。將本設(shè)備在加工時(shí)的實(shí)際工作條件作為計(jì)算仿真時(shí)的邊界條件和相關(guān)設(shè)定值，如表1所示。

表1 仿真基本參數(shù)

2.3 成形倉整體風(fēng)場仿真

為了探究成形倉工作平面上方產(chǎn)生黑煙區(qū)域的流場具體分布情況，為后續(xù)循環(huán)保護(hù)氣分支流道結(jié)構(gòu)優(yōu)化提供依據(jù)，在ICEM CFD中建立了成形倉和循環(huán)保護(hù)氣分支流道的整體流場分析模型。同樣以26 m3/h的風(fēng)場總流量進(jìn)行仿真迭代計(jì)算，得到成形倉內(nèi)風(fēng)場的速度場分布，見圖4。上進(jìn)風(fēng)口對(duì)應(yīng)的流場主要是為了保證漂浮在艙室內(nèi)的冷凝物及時(shí)從激光振鏡附近移除，減少漂浮物在激光振鏡上的吸附量；下進(jìn)風(fēng)口主要是為了去除激光掃描成形缸粉床面時(shí)產(chǎn)生的大量黑煙雜質(zhì)。

激光掃描粉床面時(shí)產(chǎn)生的黑煙主要分布在距離粉床面上方20 mm附近的區(qū)域[4]，圖5為成形缸粉床面上方20 mm(Y=4 mm)平面上的風(fēng)場速度分布云圖，云圖中左側(cè)矩形框?yàn)槌尚胃姿谄矫嬗?，即激光掃描區(qū)域，也是黑煙產(chǎn)生區(qū)域；右側(cè)矩形框?yàn)榱细姿谄矫嬗颉？擅黠@看出，在成形缸平面上方右側(cè)區(qū)域風(fēng)場的速度場明顯偏向料缸一側(cè)，沒有垂直于前后壁面，這導(dǎo)致黑煙不能被有效去除甚至有部分黑煙被吹入右側(cè)料缸的粉床面上。

1.成形缸區(qū)域(黑煙產(chǎn)生區(qū)域) 2 料缸區(qū)域圖5 成形倉黑煙區(qū)域截面速度云圖Fig.5 Velocity contour of the black smoke area

2.4 成形倉風(fēng)場仿真結(jié)果分析

將ANSYS FLUENT仿真計(jì)算得到的數(shù)據(jù)導(dǎo)入CFD POST中進(jìn)行分析。圖6所示是循環(huán)保護(hù)氣分支流道在總進(jìn)風(fēng)口流量為26 m3/h時(shí)的流場流線分布情況，可以明顯看出，經(jīng)分支流道分流后，氣體流速降低至1.5 m/s，在下進(jìn)風(fēng)口處流線分布不均勻，且在流道局部形成了渦流。進(jìn)一步分析上下進(jìn)風(fēng)口的截面流體流速，如圖7所示，下進(jìn)風(fēng)口由于橫向尺寸過大，經(jīng)分支流道分流后的流體很難在此截面上均勻分布。由圖7可以看出，靠成形倉右側(cè)(遠(yuǎn)離總進(jìn)風(fēng)口一側(cè))的流速較左側(cè)的流速偏大。

圖6 循環(huán)保護(hù)氣分支流道流場Fig.6 The flow field of gas branch

圖7 下進(jìn)風(fēng)口截面速度云圖Fig.7 The speed contour of down gas inlet

如圖1所示，下進(jìn)風(fēng)口與成形缸粉床平面中心均在垂直于成形倉前后壁面的同一直線上，因此較為理想的下進(jìn)風(fēng)口截面風(fēng)速應(yīng)當(dāng)盡可能均勻分布,才能更有效地去除激光掃描成形缸粉床面所產(chǎn)生的黑煙。為獲取下進(jìn)風(fēng)口流速分布定量數(shù)據(jù)，并為后續(xù)結(jié)構(gòu)優(yōu)化提供定量評(píng)價(jià)參考，將下進(jìn)風(fēng)口沿長度方向離散化為面積相等的6個(gè)小部分，計(jì)算各部分的平均風(fēng)速,結(jié)果如圖8所示。

圖8 下進(jìn)風(fēng)口截面流速分布Fig.8 The flow velocity distribution at the down gas inlet

將得到的下進(jìn)風(fēng)口各部分風(fēng)速進(jìn)行整理可得到表 2所示分布情況。

表2下進(jìn)風(fēng)口截面風(fēng)速分布情況

Tab.2Windspeeddistributionatdowngasinletm/s

最大風(fēng)速最小風(fēng)速平均風(fēng)速風(fēng)速分布標(biāo)準(zhǔn)差黑煙殘留情況2.158 40.976 81.532 70.516 8成形缸右側(cè)偏多

在劃分的下進(jìn)風(fēng)口的6個(gè)部分中，第5部分風(fēng)速最大，第1部分風(fēng)速最小。這說明在成形倉中成形缸粉床面上方右側(cè)部分風(fēng)場分布不均勻，實(shí)際中黑煙殘留也主要集中在這些地方,如圖9中標(biāo)注區(qū)域所示，這說明了黑煙殘留問題確實(shí)是由于下進(jìn)風(fēng)口風(fēng)速分布不均造成的。

1.下進(jìn)風(fēng)口 2.黑煙主要?dú)埩魠^(qū)域圖9 成形缸粉床面黑煙殘留情況Fig.9 Black smoke on the surface of the forming cylinder

通過對(duì)保護(hù)氣分支流道和成形倉與保護(hù)氣分支流道整體的流場仿真分析，黑煙殘留的原因可歸納為以下兩個(gè)方面：

(1)保護(hù)氣分支流道的下進(jìn)風(fēng)口不能形成流速均勻的風(fēng)場，成形缸左側(cè)區(qū)域風(fēng)速不足;

(2)保護(hù)氣下進(jìn)風(fēng)口形成的風(fēng)場在局部(右側(cè))出現(xiàn)傾斜，使得下進(jìn)風(fēng)口對(duì)面的回風(fēng)軌道很難將這部分黑煙排出。

3 保護(hù)氣分支流道結(jié)構(gòu)改進(jìn)

從上述保護(hù)氣分支流道和成形倉整體流場仿真可看出，要想有效去除激光熔化粉末過程中產(chǎn)生的黑煙，就要使成形倉下進(jìn)風(fēng)口的風(fēng)場分布更加均勻，并使其在成形缸上方產(chǎn)生的風(fēng)場盡量垂直于成形倉前后壁面。

為使得下進(jìn)風(fēng)口風(fēng)場分布均勻，考慮將成形倉的上下進(jìn)風(fēng)口分開單獨(dú)送風(fēng)，這里只針對(duì)下進(jìn)風(fēng)口結(jié)構(gòu)進(jìn)行改進(jìn)設(shè)計(jì)。將下進(jìn)風(fēng)口流道單獨(dú)設(shè)計(jì)成對(duì)稱的分叉結(jié)構(gòu),如圖10所示,內(nèi)部用5組柵格沿出風(fēng)口長度方向?qū)⑾鲁鲲L(fēng)口截面分成6個(gè)部分，從左到右依次命名為S1～S6。

圖10 改進(jìn)后的下進(jìn)風(fēng)口流道結(jié)構(gòu)Fig.10 Optimized structure of down gas inlet

改進(jìn)后的下進(jìn)風(fēng)口流道在出風(fēng)口處保持各部分截面積相同，通過調(diào)整上游柵格的分布間距進(jìn)行流場仿真迭代計(jì)算,以得到流道柵格最佳間距布置形式，使得S1～S6位置處平均風(fēng)速標(biāo)準(zhǔn)差盡可能小。原流道和優(yōu)化流道仿真迭代結(jié)果如表3所示。上述各柵格間距值如表 4所示(S1～S6對(duì)應(yīng)上游柵格間距)。

表3 下進(jìn)風(fēng)口截面風(fēng)速分布對(duì)比

表4 柵格間距系列

(a)原廠狀態(tài)

(b)有柵格(間距1)

(c)有柵格(間距2)

(d)有柵格(間距3)

(e)有柵格(間距4)

(f)有柵格(間距5)圖11 下進(jìn)風(fēng)口截面風(fēng)速分布Fig.11 Wind speed distributions at the down gas inlet

圖11所示為原廠無柵格狀態(tài)和有柵格狀態(tài)的下進(jìn)風(fēng)口截面風(fēng)場風(fēng)速分布情況，在柵格間距3情況下，下進(jìn)風(fēng)口截面風(fēng)速分布標(biāo)準(zhǔn)差最小，能達(dá)到較為理想的風(fēng)速均勻分布效果。

對(duì)下進(jìn)風(fēng)口流道結(jié)構(gòu)進(jìn)行改進(jìn)優(yōu)化，使得下進(jìn)風(fēng)口截面風(fēng)速分布更加均勻，在成形倉中激光掃描區(qū)域產(chǎn)生均勻的風(fēng)場，從而更有效地去除激光燒結(jié)金屬粉末時(shí)產(chǎn)生的黑煙。如圖 12所示，在柵格間距3時(shí)的流道柵格布置方式能在下進(jìn)風(fēng)口截面產(chǎn)生比較均勻的風(fēng)場，截面風(fēng)速分布標(biāo)準(zhǔn)差降至0.182 1 m/s，相比原無柵格狀態(tài)時(shí)的截面風(fēng)速分布有了顯著的提高。成形缸粉床面上方20 mm處風(fēng)場分布明顯較圖5有了很大改進(jìn)，風(fēng)場風(fēng)速分布更加均勻。

圖12 改進(jìn)后成形倉黑煙區(qū)域截面速度云圖Fig.12 Speed contour of the black smoke area of the optimized one

利用ANSYS FLUENT對(duì)成形倉保護(hù)氣系統(tǒng)進(jìn)行多次仿真模擬驗(yàn)證，根據(jù)仿真結(jié)果對(duì)成形倉循環(huán)保護(hù)氣供氣系統(tǒng)中的供氣流道進(jìn)行了結(jié)構(gòu)優(yōu)化。仿真結(jié)果表明，柵格間距3時(shí)(表 4)，下進(jìn)風(fēng)口截面風(fēng)速分布標(biāo)準(zhǔn)差值理論上從原廠狀態(tài)的0.472 7 m/s降低到0.182 1 m/s。對(duì)比圖5和圖 12在原廠狀態(tài)和優(yōu)化流道狀態(tài)下成形倉中黑煙產(chǎn)生區(qū)域風(fēng)場分布情況，明顯后者分布更加均勻，能更有效地去除激光燒結(jié)過程中產(chǎn)生的黑煙。

4 流道風(fēng)速分布實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證原廠進(jìn)風(fēng)流道結(jié)構(gòu)和優(yōu)化后的流道結(jié)構(gòu)流場仿真分析的正確性，說明優(yōu)化風(fēng)場分布實(shí)際的優(yōu)化效果，仍以規(guī)定的總風(fēng)速流量26 m3/h分別對(duì)原廠流道和優(yōu)化后的流道進(jìn)風(fēng)口進(jìn)行風(fēng)速測量分析。

優(yōu)化后的流道結(jié)構(gòu)采用光固化成形(SLA)的方法獲得。以進(jìn)風(fēng)口處的風(fēng)速分布為測量對(duì)象，圖13所示為進(jìn)風(fēng)流道風(fēng)速分布測量現(xiàn)場。

圖13 流道風(fēng)速分布測量現(xiàn)場Fig.13 Flow velocity distribution measurement site

多次測量S1～S6對(duì)應(yīng)進(jìn)風(fēng)口位置處的風(fēng)速分布，取其平均風(fēng)速值，可得到表5所示的風(fēng)速分布數(shù)據(jù)。

表5 優(yōu)化前后流道風(fēng)速分布

整理上述原廠流道和優(yōu)化流道的進(jìn)風(fēng)口S1～S6各個(gè)位置處的風(fēng)速分布實(shí)驗(yàn)結(jié)果，數(shù)據(jù)如表 6所示。

表6 流道風(fēng)速分布統(tǒng)計(jì)值

實(shí)驗(yàn)測得的原廠流道風(fēng)速分布與原廠流道仿真分析結(jié)果基本一致，優(yōu)化后的流道風(fēng)速分布與仿真結(jié)果也比較吻合。實(shí)驗(yàn)表明：在流道柵格間距為表4中的間距3時(shí)，下進(jìn)風(fēng)口處產(chǎn)生更加均勻的風(fēng)場。

5 結(jié)論

(1)成形倉黑煙殘留問題是由于成形倉內(nèi)風(fēng)場風(fēng)速分布不均勻造成的。下進(jìn)風(fēng)口形成的風(fēng)場左側(cè)風(fēng)速偏低且風(fēng)場向料缸側(cè)傾斜，這與粉床面上的黑煙殘留情況是相吻合的。

(2)通過優(yōu)化下進(jìn)風(fēng)口流道結(jié)構(gòu)實(shí)現(xiàn)風(fēng)場均勻化。在供氣流道中布置分流柵格。仿真和實(shí)驗(yàn)表明，在柵格間距為13.5 mm的等間距分布時(shí)，下進(jìn)風(fēng)口截面處風(fēng)速分布標(biāo)準(zhǔn)差值明顯降低。

(3)原廠流道下進(jìn)風(fēng)口截面風(fēng)速分布標(biāo)準(zhǔn)差仿真值為0.472 7 m/s，實(shí)驗(yàn)值為0.595 6 m/s，相同條件下優(yōu)化后的流道形成的風(fēng)場風(fēng)速分布的標(biāo)準(zhǔn)差仿真值為0.182 1 m/s、實(shí)驗(yàn)值為0.196 7 m/s。優(yōu)化后風(fēng)場風(fēng)速分布明顯得到改善。