微尺度鑄件充型過程的數(shù)值模擬

任明星, 李邦盛,, 傅恒志

微尺度鑄件充型過程的數(shù)值模擬

任明星1, 李邦盛1,2, 傅恒志2

(1. 哈爾濱工業(yè)大學(xué) 微納米技術(shù)研究中心, 哈爾濱 150001；2. 哈爾濱工業(yè)大學(xué) 材料科學(xué)與工程學(xué)院, 哈爾濱150001)

修正了傳統(tǒng)Navier-Stoke方程，并利用該方程和Fluent軟件對(duì)Zn-4Al合金充填微齒輪鑄件的充型過程進(jìn)行數(shù)值模擬。結(jié)果表明：在運(yùn)動(dòng)慣性的作用下，Zn-4Al合金在進(jìn)入微齒輪型腔后途徑齒輪盤部位時(shí)，合金并未橫向擴(kuò)展，而是保持入射狀態(tài)首先填充對(duì)面的齒輪軸，撞擊型腔壁后，產(chǎn)生二次壓頭，然后再向齒輪外圍的各齒部位反充；型腔內(nèi)氣體被高速運(yùn)動(dòng)的金屬液攪拌和切割，形成許多微氣泡，然后被帶入主流區(qū)，從排氣道排除。關(guān)鍵詞：ZnAl合金；微齒輪鑄件；數(shù)值模擬；充型過程

隨著微機(jī)械研究的深入和應(yīng)用領(lǐng)域的擴(kuò)展，三維復(fù)雜形狀金屬微構(gòu)件的制備受到越來越廣泛的重視[1]。2002年前后，借助于傳統(tǒng)鑄造工藝在成形三維復(fù)雜構(gòu)件方面的優(yōu)勢(shì)，BAUMEISTER等[2?5]提出了微熔模精密鑄造工藝，并制備了輪廓尺寸在亞毫米級(jí)的微齒輪和微渦輪。隨后，CHUNG等[6]、NOGUCHI和ABE[7]及HIROYUKI和MASAO[8]也相繼提出了具有不同特點(diǎn)的微熔模精鑄工藝。任明星等[9?11]則提出了生產(chǎn)效率更高的金屬型微精密鑄造工藝，也制備了整體尺寸在微米級(jí)的微齒輪鑄件。上述研究多側(cè)重于工藝開發(fā)與產(chǎn)品試制，對(duì)于微鑄件成形過程的理論研究稍顯滯后。由于微鑄件型腔尺寸微小，一些研究液態(tài)金屬充型流動(dòng)過程的(如水模擬、示蹤粒子及鋪設(shè)熱電偶等)傳統(tǒng)方法不再適應(yīng)。而近年來的研究[12?14]表明，用計(jì)算機(jī)數(shù)值模擬充型流動(dòng)過程的流場(chǎng)和溫度場(chǎng)分布，其方法可行，數(shù)值模擬結(jié)果準(zhǔn)確可靠。

在微鑄造成形工藝中，諸多在傳統(tǒng)鑄造過程中可忽略的因素，如表面張力和氣體反壓力等，影響非常顯著，對(duì)液態(tài)金屬的充型流動(dòng)過程產(chǎn)生較大影響。任明星等[15]曾對(duì)微尺度管道內(nèi)高溫液態(tài)金屬的微流動(dòng)規(guī)律進(jìn)行了數(shù)值模擬研究，結(jié)果表明，微流動(dòng)過程中出現(xiàn) 3個(gè)顯著微尺度效應(yīng)——附面層相對(duì)厚度的增大、“凸進(jìn)效應(yīng)”的形成和負(fù)壓力梯度區(qū)的產(chǎn)生。由此可見，微尺度條件下微鑄件的充型流動(dòng)過程與宏觀尺度鑄件的充型流動(dòng)過程不同，本文作者利用流體力學(xué)模擬軟件 Flunet，以微齒輪鑄件為模型，對(duì)液態(tài)金屬(Zn-4Al合金)在微齒輪鑄件型腔內(nèi)的充型流動(dòng)過程進(jìn)行數(shù)值模擬，為進(jìn)一步改進(jìn)微鑄件的微鑄造工藝奠定理論基礎(chǔ)。

1 數(shù)學(xué)模型

1.1 傳統(tǒng)流動(dòng)數(shù)學(xué)模型

在微鑄造成形過程中，可認(rèn)為液態(tài)金屬的流動(dòng)行為是非等溫及牛頓流體的非穩(wěn)態(tài)流動(dòng)。因此，本文作者在傳統(tǒng)宏觀流動(dòng)模型的基礎(chǔ)上，考慮微尺度空間對(duì)熔體黏度的影響，并考慮表面張力及氣體反壓力對(duì)微流體流動(dòng)行為的作用，建立能描述微尺度空間內(nèi)液態(tài)金屬微流動(dòng)行為的物理和數(shù)學(xué)模型。

根據(jù)實(shí)際情況，物理模型需進(jìn)行以下假設(shè)和簡(jiǎn)化：1) 液態(tài)金屬在整個(gè)流動(dòng)過程中為連續(xù)的不可壓縮流體；2) 流體為純黏性的牛頓流體；3) 忽略重力的影響；4) 表面張力系數(shù)不隨溫度變化。

一般地，描述宏觀尺度下液態(tài)金屬充型流動(dòng)的質(zhì)量守恒方程、Navier-Stokes (N-S)動(dòng)量守恒和能量方程分別為

為了確定自由表面的形狀，采用VOF(Volume of fluid model)模型，得到體積函數(shù)方程：

式(1)～(4)中：u、v和 w 為流體流速(m/s)；ρφM為作用在單位體積流體上的質(zhì)量力(N)；μ為流體動(dòng)力學(xué)黏度(Pa·s)；p 為流場(chǎng)中(x，y，z)點(diǎn)的壓力(Pa)；cp為流體的定壓比熱容(J/kg·K)；T為流體的熱力學(xué)溫度(K)；q˙為熱源能量(J)；φ為體積分?jǐn)?shù)。

1.2 微尺度下微流動(dòng)的附加項(xiàng)

由于微構(gòu)件尺寸微小，一些宏觀流動(dòng)中忽略的表面力(氣體反壓力和表面張力)在微尺度下的權(quán)重增加，不能忽略。此外，微尺度下黏性耗散的非線性影響增加，不能簡(jiǎn)單視為常數(shù)[15]。

1.3 微流體微流動(dòng)數(shù)學(xué)模型

綜上所述，在液態(tài)金屬微尺度充型流動(dòng)條件下，N-S方程應(yīng)修正為

式中：pf為氣體反壓力(Pa)；pσ為表面張力引起的反壓力(Pa)。

1.4 微流體流動(dòng)模型的邊界條件

邊界條件是在流體運(yùn)動(dòng)邊界上控制方程應(yīng)該滿足的條件，邊界條件的設(shè)定將對(duì)數(shù)值計(jì)算產(chǎn)生重要的影響。

運(yùn)動(dòng)邊界條件如下：1) 流動(dòng)前沿p=0；

能量邊界條件如下：

1) 入口處T=Tin；

3) 壁面邊界(z=h) T=Twall。

2 物理模型

2.1 微齒輪構(gòu)件的物理模型

圖1 微齒輪實(shí)體模型Fig.1 Solid model of micro-gear

應(yīng)用Pro/ENGINEER Wildfire2.0軟件繪出微齒輪的三維實(shí)體模型，如圖1所示。其基本數(shù)據(jù)如下：齒輪軸直徑 300 μm、軸長(zhǎng) 400 μm、齒輪盤直徑 580 μm、齒輪盤厚度300 μm。初始狀態(tài)為入口端充滿金屬液，型腔內(nèi)充滿少量空氣，排氣孔處連通大氣。

2.2 微構(gòu)件實(shí)體網(wǎng)格剖分

將三維實(shí)體模型導(dǎo)入 Gambit2.2.30，并對(duì)流動(dòng)計(jì)算區(qū)域進(jìn)行有限元網(wǎng)格劃分。如圖2所示，網(wǎng)格單元為8節(jié)點(diǎn)的六面體單元網(wǎng)格，應(yīng)用截面掃描法對(duì)微齒輪三維實(shí)體進(jìn)行幾何離散，共獲得節(jié)點(diǎn)109 413個(gè)、六面體100 350個(gè)。

將微齒輪物理模型的節(jié)點(diǎn)與網(wǎng)格劃分信息導(dǎo)入Fluent軟件，進(jìn)行數(shù)值求解。選擇分離解算器，采用一階隱式算法，壓力?速度耦合采取PISO方法，動(dòng)量方程采用 QUICK格式離散。連續(xù)性方程和動(dòng)量方程收斂殘差標(biāo)準(zhǔn)均為1×10?3；采用壓力入口邊界條件，以質(zhì)量流入口為計(jì)算起始點(diǎn)進(jìn)行計(jì)算。材料選擇常用的工業(yè)Zn-4Al合金，其物性參數(shù)見表1。初始條件如下：充型壓力為25 MPa，出口壓力為大氣壓，型腔溫度為373 K，澆注溫度為703 K。

圖2 微齒輪網(wǎng)格剖分示意圖Fig.2 Schematic diagram of mesh generation of micro-gear

表1 ZnAl合金模擬的物性參數(shù)Table 1 Physical parameters of ZnAl alloy for simulation

3 結(jié)果及討論

3.1 三維數(shù)值模擬結(jié)果分析

圖3所示為微齒輪在不同時(shí)刻充型流動(dòng)的三維模擬結(jié)果，以密度場(chǎng)表征。由圖3(a)和(b)可見，金屬液首先充填齒輪軸，然后反沖齒輪盤；由圖3(c)和(d)可見，金屬液在回填過程中，型腔內(nèi)氣體被高速金屬液攪動(dòng)、割裂，并隨著金屬液的運(yùn)動(dòng)而運(yùn)動(dòng)。三維模擬的優(yōu)點(diǎn)在于具有立體效果、直觀、可視性好，缺點(diǎn)是計(jì)算量大、耗時(shí)、且只能看到三維實(shí)體表面的流動(dòng)情況，看不到內(nèi)部金屬液的具體運(yùn)動(dòng)過程。為分析液態(tài)金屬的充型形態(tài)，采用二維剖面的方式進(jìn)行分析，即在z方向上選取一個(gè)截面，考察微齒輪鑄件內(nèi)部金屬液的微流動(dòng)過程。

圖3 不同時(shí)刻微鑄件充型流動(dòng)過程的三維密度場(chǎng)模擬結(jié)果Fig.3 3D simulation results (density ρ) of filling process of micro-casting mould at different times: (a) 6 μs; (b) 15 μs; (c) 22 μs;(d) 27 μs

3.2 微鑄型內(nèi)液態(tài)金屬充型流動(dòng)數(shù)值模擬

圖4～8所示為結(jié)合密度分布圖及速度矢量圖給出的液態(tài)金屬充填微齒輪鑄型的全過程。從微鑄件結(jié)構(gòu)看，液態(tài)金屬熔體經(jīng)過左側(cè)的齒輪軸、齒輪盤體以及右側(cè)的齒輪軸3個(gè)部位。由圖4(a)與5(a)可見，液態(tài)金屬由入口處進(jìn)入微型腔后，耗時(shí)8 μs左右充填左側(cè)齒輪軸；在充填齒輪盤時(shí)，雖然沒有鑄型管壁的約束，但并不向齒輪盤外圍的輪齒處橫向充填，而是保持初始入射狀態(tài)的運(yùn)動(dòng)方向，繼續(xù)充填右側(cè)的齒輪軸，直至充填到鑄型的后壁出口處，才開始反向充填齒輪盤空間。由圖4(b)與5(b)所示的速度矢量圖可知，液態(tài)金屬的流動(dòng)速度很大，因此，液態(tài)金屬流動(dòng)慣性很大，使液態(tài)金屬能夠保持原有的流向，越過齒輪盤繼續(xù)向前充填。此外，從速度矢量圖中可以看出，齒輪盤中的空氣在高速金屬熔體帶動(dòng)下，產(chǎn)生了速度較低的旋流。

圖6顯示流動(dòng)前沿的液態(tài)金屬撞擊微鑄型后壁，由于壁面的反作用力，壓力升高，形成二次壓頭，金屬熔體開始向齒輪外圍的各齒部位回流充填。由圖6(c)的壓力分布圖可以看到，在出口處的金屬液壓力最高。由圖6(b)和圖7(b)所示的速度矢量圖可以看出，在反向填充的金屬液前沿，氣體的速度非?？?，說明此處液態(tài)金屬對(duì)氣體的壓迫最強(qiáng)，在壓力的作用下，各輪齒內(nèi)部氣體漩渦流動(dòng)加劇，當(dāng)旋轉(zhuǎn)的氣團(tuán)與新流進(jìn)的液態(tài)金屬相遇時(shí)，在碰撞阻礙及壁面阻礙的作用下，輪齒內(nèi)的氣團(tuán)被碎裂成很多的微氣孔。同時(shí)，部分氣體在液態(tài)金屬的攜帶下不斷由排氣道排出。

隨著充填流動(dòng)過程的進(jìn)行，旋流金屬不斷壓迫、割裂氣團(tuán)。如圖8所示，在微鑄型中心部位液態(tài)金屬的流動(dòng)速度最快，不斷旋轉(zhuǎn)的氣團(tuán)次之，與壁面接觸部位的液態(tài)金屬，由于壁面無滑移及溫度下降造成黏性力增加的影響，其流動(dòng)速度顯著降低。由伯努力方程可知，流速增加，壓強(qiáng)降低，所以，快速流動(dòng)的液態(tài)金屬帶動(dòng)微氣孔向壓力較低的中心區(qū)域運(yùn)動(dòng)，在金屬液的壓迫下被拉長(zhǎng)，被主流區(qū)的高速流體帶動(dòng)，向出口端排氣口流動(dòng)，并隨著溢出排氣口的金屬液排出型腔，這一點(diǎn)在速度矢量圖中表現(xiàn)得更為明顯(見圖8(b))，圖中深色的箭頭代表高速流體，而齒根附近部分向出口方向的流線正好對(duì)應(yīng)密度圖中的空氣部分，表明這部分氣體由漩渦流動(dòng)逐漸被帶進(jìn)主流方向。

圖4 7 μs時(shí)刻微鑄件充型流動(dòng)過程的模擬結(jié)果Fig.4 Simulation results of filling process of micro-casting mould at 7 μs: (a) Distribution of density; (b) Distribution of velocity vector

圖5 10 μs時(shí)刻微鑄件充型流動(dòng)過程的模擬結(jié)果Fig.5 Simulation results of filling process of micro-casting mould at 10 μs: (a) Distribution of density; (b) Distribution of velocity vector

圖6 20 μs時(shí)刻微鑄件充型流動(dòng)過程的模擬結(jié)果Fig.6 Simulation results of filling process of micro-casting mould at 20 μs: (a) Distribution of density; (b) Distribution of velocity vector; (c) Distribution of pressure

通過以上分析可知，在微鑄件的充型過程中，液態(tài)金屬由入口進(jìn)入后，途經(jīng)齒輪盤部位時(shí)不橫向擴(kuò)散填充，而是直接射向后面的型腔壁，撞擊后產(chǎn)生二次壓頭，壓迫金屬液回流填充齒輪盤；液態(tài)金屬帶動(dòng)型腔內(nèi)的氣體旋轉(zhuǎn)，不斷割裂氣孔，并由氣孔向中心主流區(qū)域施壓，最后隨高速流體經(jīng)排氣孔排出。如不考慮金屬液的凝固，則上述過程不斷進(jìn)行，直至氣體完全排除，但實(shí)際鑄造過程，金屬液不但凝固，而且凝固非常迅速，所以，必然對(duì)充型過程及氣體的排出造成影響，這些問題有待進(jìn)一步研究。

圖7 30 μs時(shí)刻微鑄件充型流動(dòng)過程的模擬結(jié)果Fig.7 Simulation results of filling process of micro-casting mould at 30 μs: (a) Distribution of density; (b) Distribution of velocity vector

圖8 80 μs時(shí)刻微鑄件充型流動(dòng)過程的模擬結(jié)果Fig.8 Simulation results of filling process of micro-casting mould at 80 μs: (a) Distribution of density; (b) Distribution of velocity vector

4 結(jié)論

1) 通過增加氣體反壓力和毛細(xì)管力等微流動(dòng)時(shí)必須考慮的附加項(xiàng)，對(duì)描述流體宏觀流動(dòng)的N-S方程進(jìn)行修正，建立適合微尺度下液態(tài)金屬微流動(dòng)的 N-S方程以及物理和數(shù)學(xué)模型，并利用Fluent軟件實(shí)現(xiàn)微尺度鑄件充型過程的計(jì)算機(jī)數(shù)值模擬。

2) 金屬液保持入射狀態(tài)首先填充對(duì)面的齒輪軸，撞擊型腔壁后，產(chǎn)生二次壓頭，然后再向齒輪外圍的各齒部位反充，形成回流現(xiàn)象。

3) 型腔內(nèi)氣體被高速運(yùn)動(dòng)的金屬液攪拌、切割，形成許多微氣泡，然后被帶入主流區(qū)，從排氣道排除。

REFERENCES

[1] RUPRECHT R, BENZLER T, HANEMANN T, MüLLER K,KONYS J, PIOTTER V, SCHANZ G, SCHMIDT L, THIES A,WOLLMER H, HAUBELT J. Various replication techniques for manufacturing three-dimensional metal microstructures[J].Microsystem Technology, 1997, 4: 28?31.

[2] BAUMEISTER G, MUELLER K, RUPRECHT R, HAUSSELT J. Production of metallic high aspect ratio microstructures by microcasting[J]. Microsystem Technology, 2002, 8: 105?108.

[3] BAUMEISTER G, RUPRECHT R, HAUSSELT J. Microcasting of parts made of metal alloys[J]. Microsystem Technology, 2004,10: 261?264.

[4] BAUMEISTER G, RUPRECHT R, HAUSSELT J. Replication of LIGA structures using microcasting[J]. Microsystem Technology, 2004, 10: 484?488.

[5] BAUMEISTER G, HAUSSELT J, ROTH S, RUPRECHT R.Microcasting[C]//Advanced micro and nanosystems:Micro-engineering in metals and ceramics. Weinheim:Wiley-VCH, 2005: 357?393.

[6] CHUNG S, PARK S, LEE L, JEONG H. Replication techniques for a metal microcomponent having real 3D shape[J].Microsystem Technology, 2005, 11: 424?428.

[7] NOGUCHI H, ABE S. Study on microcasting[J]. International Journal of the Japan Society for Precision Engineering, 2003, 69:125?129.

[8] HIROYUKI N, MASAO M. Three-dimensional microcasting[C]//Proc of the 10th ICPE. New York, 2001:349?353.

[9] 任明星, 李邦盛, 楊 闖, 傅恒志. 納米壓痕法測(cè)定微鑄件硬度及彈性模量[J]. 中國有色金屬學(xué)報(bào), 2008, 18(1): 231?236.REN Ming-xing, LI Bang-sheng, YANG Chuang, FU Heng-zhi.Hardness and elastic modulus of microcastings by nanoindentation[J]. The Chinese Journal of Nonferrous Metals,2008, 18(1): 231?236.

[10] LI Bang-sheng, REN Ming-xing, YANG Chuang, FU Heng-zhi.Microstructure of Zn-Al4 alloy microcastings by micro precision casting based on metal mold[J]. Transactions of Nonferrous Metals Society of China, 2008, 18(2): 327?332.

[11] 任明星, 李邦盛, 楊 闖, 傅恒志. 納米壓痕法測(cè)定微鑄件室溫蠕變速率敏感指數(shù)[J]. 金屬學(xué)報(bào), 2008, 44(3): 272?276.REN Ming-xing, LI Bang-sheng, YANG Chuang, FU Heng-zhi.Measurement of creep rate sensitivity of microcastings at room temperature by using nanoindentation[J]. Acta Metallurgica Sinica, 2008, 44(3): 272?276.

[12] 周彼德, 薛 祥, 糜忠蘭, 孫小波, 張春暉, 馬 建. 鑄件三維充型過程耦合數(shù)值模擬[J]. 中國有色金屬學(xué)報(bào), 2000, 10(2):230?233.ZHOU Bi-de, XUE Xiang, MI Zhong-lan, SUN Xiaobo, ZHANG Chun-hui, MA Jian. 3D coupling numerical simulation of mold filling[J]. The Chinese Journal of Nonferrous Metals, 2000, 10(2): 230?233.

[13] WU Shi-ping, LI Chang-yun, GUO Jing-jie, SU Yan-qing, LEI Xiu-qiao, FU Heng-zhi. Numerical simulation and experimental investigation of two filling methods in vertical centrifugal casting[J]. Transactions of Nonferrous Metals Society of China,2006, 16(5): 1035?1040.

[14] 陳興潤, 張志峰, 徐 駿, 石力開. 電磁攪拌法制備半固態(tài)漿料過程電磁場(chǎng)、流場(chǎng)和溫度場(chǎng)的數(shù)值模擬[J]. 中國有色金屬學(xué)報(bào), 2010, 20(5): 937?945.CHEN Xing-run, ZHANG Zhi-feng, XU Jun, SHI Li-kai.Numerical simulation of electromagnetic field, flow field and temperature field in semi-solid slurry preparation by electromagnetic stirring[J]. The Chinese Journal of Nonferrous Metals, 2010, 20(5): 937?945.

[15] 任明星, 李邦盛, 楊 闖, 傅恒志. 微尺度型腔內(nèi)液態(tài)金屬流動(dòng)規(guī)律模擬研究[J]. 物理學(xué)報(bào), 2008, 57(8): 5063?5071.REN Ming-xing, LI Bang-sheng, YANG Chuang, FU Heng-zhi.Simulation research on the law of flow of liquid metal in micro-channels[J]. Acta Physica Sinica, 2008, 57(8):5063?5071.

Numerical simulation of filling process on casting in micro scale

REN Ming-xing1, LI Bang-sheng1,2, FU Heng-zhi2
(1. Research Center of Micro/nano Technology, Harbin Institute of Technology, Harbin 15001, China；2. School of Materials Science and Engineering, Harbin Institute of Technology, Harbin 150001, China)

The filling process of micro-gear casting with Zn-4Al alloy was simulated by applying the modified traditional Navier-Stoke equation and the Fluent software. The results show that the alloy is not horizontal expansion and keeps the initial state when the liquid metal is flowed into the cavity of micro-gear and passes through parts of gear plate by the effect of inertia. Liquid metal is anti-filled into the tooth parts outside gear by impacting posterior of cavity, resulting in the second pressure head. Many micro-bubbles are created after the gas in the cavity is mixed and cut by high-speed movement of the metal liquid, and then brought into the mainstream areas and excluded from the exhaust ports.

ZnAl alloy; micro gear casting; numerical simulation; filling process

TG21

A

1004-0609(2011)07-1675-06

國家自然科學(xué)基金資助項(xiàng)目(50475028)；黑龍江省自然科學(xué)基金資助項(xiàng)目(ZJG04-01)；微系統(tǒng)與微結(jié)構(gòu)制造教育部重點(diǎn)實(shí)驗(yàn)室開放基金資助項(xiàng)目(HIT.KLOF.2009012)

2010-07-08；

2010-10-28

任明星，講師，博士；電話：0451-86403268；E-mail：hitrmx@163.com

(編輯 陳衛(wèi)萍)