分析熔于一冷双组分金属粉末层与恒热流外文文献翻译、中英文翻译
分析熔于一冷双组分金属粉末层与恒热流外文文献翻译、中英文翻译,分析,一冷双,组分,金属粉末,热流,外文,文献,翻译,中英文
subcooled Yuwen 11 is investigated points. solid approximate solution, respectively. The eects of porosity, Stefan number, and subcooling on the surface temperature and solid liquid interface are also investigated. The present work provides a strong foundation upon which the investigation of complex selectively fusing a thin layer of the powders with scan- imate diameter of the laser beam. The particular melting and resolidification induced by a directed laser sity change of the powder bed accompanies the melting process. Melting and solidification in 1-D semi-infinite body with density change under the boundary condition of the first kind have been investigated by Zckert and Drake 9, Crank 10, Carslaw and Jaeger 11 and * Corresponding author. Tel.: +1 573 884 6939; fax: +1 573 884 5090. E-mail address: zhangyumissouri.edu (Y. Zhang). Applied Thermal Engineering 26 1359-4311/$ - see front matter C211 2005 Elsevier Ltd. All rights reserved. ning laser beam. After sintering of a layer, a new layer of the powder is deposited in the same manner and a 3-D part can be built in a layer-by-layer process. A mixed metal powder bed, which contains two types of the metal powders possessing significantly dierent melting points, is used extensively in direct SLS of metal powders 2,3. The high melting point powder never melt in the sintering process and plays a significant role as the support structure necessary to avoid boiling phenom- enon, which is the formation of spheres with the approx- beam. It is a good starting point to investigate a simpli- fied 1-D model to get a better understanding of the melting process in direct SLS before a much more com- plicated 3-D model is investigated. Fundamentals of melting and solidification have been investigated extensively and detailed reviews are avail- able in Refs. 7,8. Melting in SLS of the metal powders is significantly dierent from the normal melting since the volume fraction of the gas in the powders decreases significantly after melting. Therefore, a significant den- three-dimensional selective laser sintering (SLS) process can be based. C211 2005 Elsevier Ltd. All rights reserved. Keywords: Melting; Metal; Powder layer 1. Introduction Direct Selective Laser Sintering (SLS) is an emerging technology of Solid Freeform Fabrication (SFF) via which 3-D parts are built from the metal-based powder bed with CAD data 1. A fabricated layer is created by material properties and methods of material analysis of the metal-based powder system for SLS applications are addressed by Storch et al. 4 and Tolochko et al. 5. Fundamental issues on direct SLS are thoroughly re- viewed by Lu et al. 6. In fabrication of near full density objects from metal powder, direct SLS is realized via Analysis of melting in a metal powder layer with Tiebing Chen, Department of Mechanical and Aerospace Engineering, University Received 1 February 2005; Available online Abstract Melting of a subcooled two-component metal powder layer mixture of two metal powders with significantly dierent melting physical model. The temperature distributions in the liquid and doi:10.1016/j.applthermaleng.2005.07.034 two-component constant heat flux Zhang * of Missouri-Columbia, Columbia, MO 65211, United States accepted 18 July 2005 October 2005 analytically. The powder bed considered consists of a Shrinkage induced by melting is taken into account in the phases are obtained using an exact solution and an integral (2006) 751765 Nomenclature c p specific heat (J kg C01 K C01 ) h sl latent heat of melting or solidification (J kg C01 ) k thermal conductivity (W m C01 K C01 ) K g dimensionless thermal conductivity of gas K s dimensionless eective thermal conductivity of unsintered powder q 00 heat flux (W m C02 ) s solidliquid interface location (m) S dimensionless solidliquid interface location s 0 location of liquid surface (m) S 0 dimensionless location of liquid surface Sc subcooling parameter Ste Stefan number t time (s) T temperature (K) w velocity of liquid phase (m s C01 ) W dimensionless velocity of the liquid phase z coordinate (m) Z dimensionless coordinate Greek symbols a thermal diusivity (m 2 s C01 ) 752 T. Chen, Y. Zhang / Applied Thermal Charach and Zarmi 12. It should be noted that melting during SLS occurs under the boundary condition of specified heat flux instead of specified temperature. Goodman and Shea 13 studied melting and solidifica- tion in the finite slab under a specified heat flux by using the heat balance integral method. Zhang et al. 14 investigated the melting problem in a subcooled semi- infinite region subjected to constant heat flux heating. Zhang et al. 15 solved melting in a finite slab with the boundary condition of the second kind by using a semi-exact method. Shrinkage formation due to density change during the solidification process in 2-D cavity was investigated numerically by Kim and Ro 16, who concluded that the density change played a more impor- tant role than convection in the solidification process. Zhang and Faghri 17 analytically solved a one- dimensional melting problem in a semi-infinite two- component metal powder bed subjected to a constant heating heat flux. Eects of the porosity of the solid phase, initial subcooling parameter and dimensionless thermal conductivity of the gas were investigated. Since SLS of the metal powder is actually a layer-by-layer pro- cess, it is necessary to investigate melting in a mixed me- tal powder bed with the finite thickness during the SLS process. In this paper, melting of the mixed powder bed with finite thickness subjected to constant heating heat flux will be investigated. C22a dimensionless thermal diusivity b parameter to distinguish between two melting cases d thermal penetration depth (m) D dimensionless thermal penetration depth e volume fraction of gas(es) (porosity for unsintered powder) h dimensionless temperature q density (kg m C03 ) s dimensionless time / volume fraction of the low melting point powder in the powder mixture Subscripts g gas i initial l liquid phase m melting point p sintered part s unsintered solid (mixture of two solid pow- ders) Engineering 26 (2006) 751765 2. Physical model The physical model of the melting problem is shown in Fig. 1. A powder bed with finite thickness contains two metal powders with significantly dierent melting points. The initial temperature of the powder bed is below the melting point of the low melting point pow- der. At time t = 0, a constant heat flux, q 00 , is suddenly applied to the top surface of the powder bed, and the bottom surface of the powder bed is assumed to be adiabatic. Since the initial temperature of the powder bed is below the melting point of the low melting point powder, its melting does not start simultaneously with the addition of heat heating. Only after a finite period of time of preheating, in which the surface tem- perature of the powder reaches the melting point of the low melting point powder, will the melting start. The powder with the high melting point will never melt during the entire process. Therefore, the problem can be subdivided into two problems: one being heat conduction during preheating and the other being melting. The physical model is considered as a conduc- tion-controlled problem. The eect of natural convec- tion in the liquid region due to the temperature dierence is not considered since the temperature is highest at the liquid surface and decreases with increas- ing z. T. Chen, Y. Zhang / Applied Thermal Engineering 26 (2006) 751765 753 2.1. Duration of preheating During preheating, pure conduction heat transfer oc- curs in the powder mixture. The governing equation and the corresponding initial and boundary conditions for the preheating problem are z Fig. 1. Physical q s s 0 H 0 a s o 2 T s oz 2 oT s ot ; 0 z H 0 ; t t m 1 T T i ; 0 z H 0 ; t 0 2 C0 k s oT s oz q 00 ; z 0; t t m 3 oT s oz 0; z H 0 ; t t m 4 2.2. Melting After melting starts, the governing equation in the liquid phase is a l o 2 T l oz 2 oT l ot w oT l oz ; s 0 z t m 5 where w is the velocity of liquid surface induced by the shrinkage. Since the liquid is incompressible, the shrink- age velocity w is w ds 0 dt ; s 0 t m 7 The governing equation for the solid phase and its cor- responding boundary conditions are a s o 2 T s oz 2 oT s ot ; st t m 9 Liquid-solid interface Low melting point powder High melting point powder model. Original surface Liquid surface oz The temperature at the solidliquid interface satisfies T l z;tT s z;tT m ; z st; t t m 10 The energy balance at the solidliquid interface is k s oT s oz C0 k l oT l oz 1 C0e s /q l h sl ds dt ; z st; t t m 11 Based on the conservation of mass at the solidliquid interface, the shrinkage velocity, w, and the solidliquid interface velocity, ds/dt, have the following relationship 17: w e s C0e l 1 C0e l ds dt 12 2.3. Non-dimensional governing equations By defining the following dimensionless variables: h l qc p p T l C0 T m Uq l h sl h s qc p p T s C0 T m Uq l h sl Sc qc p p T m C0 T i Uq l h sl ; s a p t H 2 ; Z z H S s H ; S 0 s 0 H ; D d H ; W w C1 H a p K s k s k p 1 C0e s ; K g k g k p ; C22a s a s a p Ste q 00 H Uq l h sl a p 13 The non-dimensional governing equation and the corre- sponding initial and boundary conditions for the pre- heating problem become o 2 h s oZ 2 1 C22a s C1 oh s os ; 0 Z 1; s s m 14 h C0Sc; 0 Z 1; s 0 15 oh s oZ C0 Ste K s 1 C0e s ; Z 0; s s m 24 K s oh s oZ C0 1 C0e l 1 C0e s oh l oZ dS ds ; Z Ss; s s m 25 W e s C0e l 1 C0e l dS ds ; S 0 Z t m 26 3. Approximate solutions When the top surface of the mixed metal powder bed is subjected to constant flux heating, the heat flux will penetrate through the top surface and conduct down- ward the bottom surface. The depth to which the heat flux penetrates at an instant in time is defined as the thermal penetration depth, beyond which there is no heat conduction. Goodman and Shea 13 introduced a 00 754 T. Chen, Y. Zhang / Applied Thermal Engineering 26 (2006) 751765 oh s oZ 0; Z D; s s m 18 For melting, the non-dimensional equation and corre- sponding boundary conditions are o 2 h l oZ 2 oh l os W oh l oZ ; S 0 Z s m 19 W dS 0 ds ; S 0 s m 21 o 2 h s oZ 2 1 C22a s C1 oh s os ; Ss Z s m 22 Fig. 2. Validation of analytical parameter, b = q H/2k s (T m C0 T i ), to classify two cases of melting in a finite slab. When b is greater than 1, the top surface temperature reaches the melting point in a shorter time than the thermal penetration depth reaches the bottom surface, indicating that a shorter preheating time is needed. If b is less than 1, the surface tempera- ture is still below the melting point when the thermal penetration depth has reached the bottom surface. Pre- heating continues until the top surface temperature reaches the melting point of low melting point powder. The parameter b can also be expressed using non- dimensional parameters defined in Eq. (13), i.e., b = Ste/2K s Sc(1 C0 e s ). It can be seen that the value of b is determined by four basic non-dimensional parame- ters: Stefan number Ste, subcooling parameter Sc, eec- tive thermal conductivity of the solid phase K s and solutions. T. Chen, Y. Zhang / Applied Thermal Engineering 26 (2006) 751765 755 volume fraction of gas e s in the solid phase. Preheating and melting for both b 1 will be discussed. 3.1. Preheating 3.1.1. b 1 The heat-balance integral method 18,19 is employed here. Integrating the heat-conduction Eq. (14) with re- spect to Z from 0 to D, the integral equation is obtained. oh s oZ D;sC0 oh s oZ 0;s C20C21 1 C22a s d ds H ScD27 where H R D 0 h s Z;sdZ. h s (Z, s) is assumed to be a second degree polynomial function which satisfies boundary conditions specified by Eqs. (16)(18). Then h s (Z, s) can be determined (b) Fig. 3. Eect of porosity in the liquid phase l l l l h s Z;sC0Sc Q 2K s D1 C0e s DC0 Z 2 28 The Eqs. (16)(18) and (28) can be substituted into Eq. (27) and then an ordinary dierential equation for the thermal penetration depth, D, is obtained which can be solved easily. D 6 C1 C22a s C1s 1=2 29 When the thermal penetration depth reaches the bottom surface, i.e., D = 1, the temperature distribution in the powder bed is h s Z;sC0Sc Ste 2K s 1 C0e s 1 C0 Z 2 ; 0 Z 1; s s D1 s m 30 (a) l l l l on surface temperature (Ste = 0.02). 756 T. Chen, Y. Zhang / Applied Thermal Engineering 26 (2006) 751765 which becomes the initial condition of the next stage of preheating. After the thermal penetration depth reaches the bottom, the problem becomes a conduction problem in a finite slab. In a manner analogous to that described previously, the temperature of the powder is h s Z;sC0Sc Ste 2K s 1 C0e s 1 C0 Z 2 Ste C1 C22a s K s 1 C0e s C2sC0s D1 ; 0 Z 1; s D1 0; s s m 36 h s 0;sC0Sc Ste C1 6a s s p =2 C1 K s C11 C0e s C138; Z 0; 0 s s D=1 , the problem becomes melting in a finite slab. The temperature distribution in the solid, h s (Z, s), and the liquidsolid interface location, S can be ob- tained by solving Eqs. (22)(24) using the integral approximate method identical to the case of b 1. 4. Results and discussion The validation of the analytical solution was conducted by comparing the results with the numerical results obtained from Chen and Zhang 20, who inves- tigated the two-dimensional melting and resolidification of a two-component metal powder layer in SLS process subjected to a moving laser beam. In order to use the two-dimensional code in Ref. 20 to solve melting in (a) (b) Fig. 8. Eect of subcooling on surface l a powder layer subjected to constant heat flux, the Gaussian laser beam was replaced by a constant heat- ing heat flux on the top of the entire powder bed and the laser scanning velocity was set to zero in numerical solution. The parameters used in the present paper were converted into corresponding parameters in Ref. 20 for purpose of code validation. The comparisons of instantaneous locations of liquid surface and liquid solid interface obtained by analytical and numerical solutions are shown in Fig. 2. It can be seen that the preheating time obtained by the analytical and numer- ical solutions are almost the same. The locations of liquid surface and liquidsolid interface obtained by analytical and numerical solutions move at very similar trends. The time it takes to completely melt the entire powder layer obtained from analytical solution is about l temperature (Ste = 0.15). T. Chen, Y. Zhang / Applied Thermal Engineering 26 (2006) 751765 761 l 4% longer than that obtained from the numerical solution. The eects of porosity, subcooling, dimensionless thermal conductivity and Stefan number on the surface temperature, location of the liquid surface, and the loca- tion of the solidliquid interface of the powder bed will be investigated. Fig. 3 shows how the surface tempera- ture is influenced by the porosity in the liquid phase for Ste = 0.02 and several dierent subcooling parame- ters. The eect of shrinkage is isolated by fixing the sub- cooling parameter, porosity of the solid phase, and the dimensionless thermal conductivity. It can be seen that the surface temperature increases as porosity in the liquid phase increases. This is because the eective ther- mal conductivity decreases with increasing volume frac- tion of the gas. When Sc = 0.1, the preheating time is (a) (b) Fig. 9. Eect of subcooling on the location of the liquid much shorter compared to when Sc = 3.0. The eect of shrinkage on the surface temperature for Ste = 0.15 is shown in Fig. 4. As we can see, the increase of poros- ity in the liquid phase results in higher surface tempera- tures and that a higher Sc requires a longer preheating time. When Sc = 3.0, one can observe that the duration of the melting process is shortened significantly when Ste increases from 0.02 to 0.15. Fig. 5 shows the loca- tions of solidliquid interface and liquid surface corre- sponding to the conditions of Fig. 3. The solidliquid interface moves faster when more gas is driven out from the liquid. It follows that the corresponding location of the liquid surface moves downward significantly due to the shrinkage of the mixed metal powder bed. The locations of solidliquid interface and liquid surface corresponding to the conditions of Fig. 4 are shown in l surface and the liquidsolid interface (Ste = 0.02). 762 T. Chen, Y. Zhang / Applied Thermal Engineering 26 (2006) 751765 l Fig. 6. The decrease of porosity in the liquid phase also expedites the motion of the solidliquid interface and liquid surface downward. Fig. 7 shows the eect of the initial subcooling on the surface temperature
收藏