Next Article in Journal
Correction: Wang et al. Parametric Formula for Stress Concentration Factor of Fillet Weld Joints with Spline Bead Profile. Materials 2020, 13, 4639
Next Article in Special Issue
Sintering Behavior, Thermal Expansion, and Environmental Impacts Accompanying Materials of the Al2O3/ZrO2 System Fabricated via Slip Casting
Previous Article in Journal
Ablation of Bone Tissue by Femtosecond Laser: A Path to High-Resolution Bone Surgery
Previous Article in Special Issue
Al2O3/ZrO2 Materials as an Environmentally Friendly Solution for Linear Infrastructure Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of Particulate Migration Behavior within Molten Pool during TIG-Assisted Droplet Deposition Manufacturing of SiC Particle-Reinforced Aluminum Matrix Composites

College of Digital Technology and Engineering, Ningbo University of Finance & Economics, Ningbo 315175, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(9), 2430; https://doi.org/10.3390/ma14092430
Submission received: 26 March 2021 / Revised: 24 April 2021 / Accepted: 4 May 2021 / Published: 7 May 2021
(This article belongs to the Special Issue Advances in Metal and Ceramic Matrix Composites)

Abstract

:
A transient three-dimensional (3D) numerical model was established to illustrate the heat transfer, fluid flow and particle migration behaviors in the molten pool during TIG-assisted droplet deposition manufacturing (DDM) of SiC particle-reinforced aluminum matrix composites (AMCs). The effect of temperature-dependent physical properties and the interaction between the SiC reinforcement and the liquid metal matrix were considered. A double-ellipsoidal volumetric heat source model was adopted to simulate the energy interactions between the pulse square-wave variable polarity TIG welding arc and the moving substrate. Free surface fluctuations of molten pool due to arc force and sequential droplet impact are calculated with volume of fluid (VOF) method in a fixed Eulerian structured mesh. The numerical model, capable of capturing the impact, simultaneous spread, and phase change of the droplets as well as the motion trajectory and terminate distribution state of the reinforcement particles, is key tool to understand the formation mechanism of the TIG-assisted DDM of SiC particle-reinforced AMCs. The numerical model was validated by the metallographic observations, and the calculated particle distribution and solidification morphology of deposited layer agree well with the experimental measurements.

1. Introduction

Compared to unreinforced metal matrix, ceramic-reinforced metal matrix composites (MMCs) can be utilized in extreme working conditions and some high-end engineering industries (i.e., automotive, aerospace, etc.) due to their higher mechanical, thermal, and fatigue properties at elevated temperatures [1]. The research and application of these materials such as particle-reinforced aluminum matrix composites (AMCs) will not only significantly expand the scope of potential applications but also achieve their economic potential [2]. Particle-reinforced AMCs are generally prepared by powder metallurgical technique [3] and melt-metallurgical production method [4,5,6]. Reinforced AMCs have the improved mechanical properties, such as the tensile strength [7], the yield strength [8], and the wear resistance [9]. However, it is difficult for the abovementioned preparation methods of particle-reinforced AMCs to fabricate complex-shaped components [10], but also, they have problems of high cost and energy consumptions.
Additive manufacturing (AM) is the manufacturing process of components through layer-based material consolidation (ASTM, 2012). Recently, the emerging manufacturing technology with tailorable composition and properties, lowered cost and energy usage, and high design flexibility has gained popularity in the fabrication of high-performance materials. Different research groups investigated the feasibility to use direct metal laser sintering (DMLS) to achieve Al-based metal matrix composites [11,12]. SiC particles with different size and volume fraction were employed to analyze the behaviors of the composite, demonstrating that cracking occurs during the DMLS of these composites. It is found that obvious reaction occurs between aluminum melt and SiC particles, leading to the formation of Al4SiC4 [13]. It is also proved that in any case, the presence of the ceramic reinforcement on hardness is evident. Meanwhile, there is a growing research to develop MMCs via selective laser melting (SLM) process [14,15]. However, due to the low absorption of the laser and the high thermal conductivity of the aluminum alloys, the deposition rate of aluminum alloys and Al-based composites via laser AM is very low [16]. Although wire arc additive manufacturing (WAAM) has also proven its capability to produce medium-to-large-scale aluminum components, aluminum-based filler wire containing ceramic particles are rarely seen. Therefore, it is necessary to develop a novel and highly efficient forming method to realize the additive manufacturing of complex AMCs components.
In this study, a novel hybrid additive manufacturing process, termed as TIG-assisted droplet deposition manufacturing, is proposed to fabricate SiC particle-reinforced AMCs components, the process principle of the additive manufacturing technique will be described in detail in Section 2. A 3D numerical simulation model based on finite volume method (FVM) was proposed to understand and analyze the molten pool dynamics and particle-melt interactions within the molten pool during the TIG-assisted DDM of SiC particle-reinforced AMCs. The resultant distribution state of SiC particles in the solidified aluminum matrix was experimentally examined, which is in a good agreement with the results calculated by simulation.

2. Principle of TIG-Assisted DDM Process

Figure 1 shows the schematic illustration of the TIG-assisted droplet deposition manufacturing process. The equipment of TIG-assisted DDM includes an induction-heated crucible, a piezoelectric ceramic rod, a gas protection device, and a 3D mobile platform. The TIG-assisted process can produce ceramic-reinforced AMCs components by synchronized controlling the 3D mobile platform and the injected molten aluminum droplets with ceramic reinforcements. A programmable multi-axis controller (PMAC) is used to control the motion of the mobile platform. A high-dynamic welding camera was applied to monitor the forming process. Under the combined action of the high frequency vibration by the piezoelectric ceramic rod, surface tension, and hydrostatic pressure, the molten aluminum will flow through the inner flow channel of graphite nozzle and to be injected in form of successive droplets. When the injected molten aluminum droplets impact and contact with the molten pool formed on the substrate by the pulsed variable polarity TIG welding arc, coalescence, spontaneous spreading and non-equilibrium solidification would take place. Thus, a metallurgical bonding can be achieved. The solidified deposition layers of reinforced AMCs always keep synchronous with the substrate. Variable polarity TIG welding arc has a two-fold effect in the TIG-assisted DDM: First, the oxides on the surface of aluminum alloy substrate could be removed timely. Fortunately, the VP-TIG welding arc can clean oxide film timely. Second, sound metallurgical bonding between the deposited layers could be achieved by applying the preceding VP-TIG welding arc, which substantially improve the high-temperature wettability and spreadability of molten aluminum droplets on the surface of previous deposited layers or substrate.
TIG-assisted DDM process has high deposition efficiency, low cost, great freedom in material selection, etc. It has a wide application prospects in aerospace, shipping, rail transportation and other fields. Compared to a single TIG welding and droplet deposition process, TIG-assisted DDM involves more process parameters, and thus, relying on experience or trial-and-error methods to determine optimized process parameters will be a long cycle and high-cost problem. With the development of computer technology and numerical calculation methods, using computer simulation technology to study the TIG-assisted DDM process of particle-reinforced AMCs has become possible, which is helpful to reveal the forming mechanism and optimize process parameters.
In this work, the transient thermal fluid fields and particles migration behavior during the deposition process were calculated and analyzed using a 3D finite volume method model. The molten pool dynamics and the kinetic characteristics of SiC particles within the molten pool were well understood by visualizing the thermo-fluid field. This will provide theoretical basis for a deeper understanding of the mechanism and engineering application of TIG-assisted DDM.

3. Description of the Mathematical Model

3.1. Physical Description of TIG-Assisted DDM Process

TIG-assisted DDM of the particle-reinforced AMCs involves complex non-equilibrium metallurgical processes, which exhibits impacting and spreading of molten aluminum droplets with SiC particles, free-surface fluctuation, and unsteady heat transfer behaviors accompanied with solidification phase change. Arc force, buoyant force, and surface tension gradient force produce vortex flows in the melt, which directly affect the flow pattern and the associated physical metallurgy process. After cooling and solidification of the melt, the deposited layer has a definite distribution of SiC particles. The calculations of the morphology of deposited layer and the resultant distribution state of reinforcing particles are the important goals of this work.

3.2. Basic Assumptions

The melt motion induced by the impacting of molten aluminum droplets on the molten pool by the preceding TIG welding arc can be assumed to be laminar with a viscous incompressible heat-conducting fluid. In fluid dynamics, the Boussinesq approximation in a 3D Cartesian coordinate system is used in the field of buoyancy-driven flow. The influence of shielding gas on the flow pattern in the molten pool is ignored. The volume thermal expansion of the SiC particles is neglected. The fraction of particles in the liquid metal matrix is small, and thus, the particles collision with each other, and the influence of particles on the melt flow and surface tension can also be ignored. The thermal physical properties of the SiC particles are assumed to be constant.

3.3. Governing Equations

The governing equations describing the heat-mass transport of molten pool- impacting droplets with SiC particles include the mass, momentum, and energy conservation equations, as shown in Equations (1), (3), and (5). The morphology of SiC particle-reinforced AMCs deposited layers was calculated by adopting the volume of fluid (VOF) method, and the corresponding equation is shown in Equation (6).
Mass conservation equation,
V = m ˙ s ρ
where V is the velocity vector of melt, m ˙ s is the rate of the mass source of droplets, and ρ is density.
In the TIG-assisted DDM, molten aluminum droplets will be periodically added to the calculation model as a mass source term. It can be expressed as Equation (2).
m s = ρ d Δ t
where ρd is the droplet density, and ∆t is the time interval between adjacent droplets.
Momentum conservation equation,
V t + V V = 1 ρ P + μ k 2 V K V + m ˙ s ρ ( V s V ) + G
where P is the hydrodynamic pressure, μk is the dynamic viscosity, V s is the velocity vectors for the mass sources, G is the body acceleration due to body force (e.g., gravity, Lorentz force, and buoyancy forces), K is the deceleration coefficient in the semisolid zone, and K can be written as
K = C [ F s 2 ( 1 F s ) 3 + B ]
where C is a constant representing semisolid zone morphology, Fs is the solid fraction, B* is the positive zero used to avoid division by zero. The Carman–Kozeny equation is used to model the flow in the semisolid zone, which is based on the derivation from the Darcy model and the assumption that this region is able to be treated as a porous media.
Energy conservation equation,
h t + V h = 1 ρ ( k T ) + h ˙ s + S q
where k is the thermal conductivity of the material, h is the specific enthalpy, T is the temperature, Sq is the other energy source terms, and h ˙ s is the specific enthalpy source, which is associated with the mass source in Equation (1).
Volume of fluid (VOF) equation,
F t + ( V F ) = F s
where F is the volume fraction and Fs is the volume source from the molten droplets.
Solidification model.
The enthalpy method is used to solve the melting and solidification problem. The internal energy equation of the system is
h t + v h = 1 ρ ( k T )
where k is the thermal conductivity of the material, ρ is density, T is time, and h is enthalpy; it can be expressed as
h = { ρ s C s T ( T T s ) h ( T s ) + h s l T T s T l T s ( T s T T l ) h ( T l ) + ρ l C l ( T T l ) ( T l T )
where ρs and ρl are solid and liquid densities; Cs and Cl are solid and liquid specific heats, respectively; Ts and Tl are solid and liquid temperatures of materials, and hsl is latent heat of fusion.

3.4. Driving Forces

The main driving forces acting on the molten pool are Lorentz force, arc pressure, surface tension, gravity, and buoyancy force. In this study, the arc pressure is approximated as Gaussian density distribution, as given below
P a r c ( x , y ) = μ 0 I 2 4 π 2 σ r 2 exp ( ( r 2 2 σ r 2 ) )
where μ0 is the magnetic permeability, and σr is the distribution parameter of arc pressure.
The arc drag force acting on the molten pool surface can be represented as follows.
P D r a g ( r ) = P M a x r r S h e a r exp ( r r S h e a r ) 2
where rShear is the distribution parameter of arc drag force.
The Lorentz force is considered by adopting the elliptically symmetric welding current density in the EMF model [17]. The equations relating EMF are listed as below.
F x = J z × B θ x r a
F z = J r × B θ
J z = I 2 π 0 λ J 0 ( λ r a ) exp ( λ 2 σ x 2 ) sinh [ λ ( c z ) ] sinh ( λ c ) d λ
J r = I 2 π 0 λ J 1 ( λ r a ) exp ( λ 2 σ x 2 ) cosh [ λ ( c z ) ] sinh ( λ c ) d λ
B θ = μ m I 2 π 0 J 1 ( λ r a ) exp ( λ 2 σ x 2 ) sinh [ λ ( c z ) ] sinh ( λ c ) d λ
r a = ( x x 0 ) 2 + [ σ x σ y ( y y 0 ) ] 2
where I is the welding current, Fi are the components of Lorentz force in i-direction (i = x, y, z), Jz and Jr are the axial and radial current density, Bθ is the angular component of magnetic field, J0 and J1 are the zero order and one order Bessel function, respectively.
The surface tension of the liquid AMCs can be expressed as a linear function of temperature:
γ = γ 0 + k sur ( T T 0 )
where γ0 is the surface tension coefficient at T0, ksur is the temperature coefficient of surface tension.

3.5. Droplet Model with SiC Particles

In the TIG-assisted DDM process, the reinforcing phase in the deposited layer is introduced as the form of SiC particle within molten aluminum droplets. Droplet transition is achieved through the coupling drive of piezoelectric ceramic excitation and the flow channel of graphite nozzle. The relevant parameters of the droplets and the SiC particles are shown in Table 1.
The motion control equation of the particle-reinforced phase is:
d u p d t = 1 ρ p P + g + β ( u u ) | u u | ρ ρ p
where u = u p + u d i f f u s i o n , up is the average velocity of SiC particles, ρp is the average density of SiC particles, g is the total volume force to which SiC particles are subjected, u and P are respectively the flow velocity and pressure of the melt around the particle phase.
The motion trajectories and thermodynamics properties of SiC particles are described by Lagrangian method. It has been proved that the trajectory approach can be effectively used to simulate the transport properties of dispersed phase.
For an individual particle p with a radius rp, the velocity vector V p = ( u p , v p , w p ) of the mass center and the volume-averaged temperature of the particle are considered. The related equations of motion and heat transfer for the individual spherical particle p are written in the form
d V p d t = 3 8 ρ C d ρ p r p | V V p | ( V V p ) + g ρ p ρ ρ p
d X p d t = V p
ρ p 4 3 π r p 3 d H p d t = 4 π r p 2 [ k m 2 r p N u p ( T p T ) ]
T p = { H p / c s , H p < c s T m p T m p , c s T m p H p c s T m p + H m p T m p + ( H p c s T m p H m p ) / c m , H p > c s T m p + H m p
C d ( Re p ) = 24 Re p ( 1 + 0.179 Re p 0.5 + 0.013 Re p )
Re p 10 3 , Re p = 2 r p ρ m | V V p | μ
N u p = 2 + 0.459 Re p 0.55 Pr 0.33 ; N u p = ( 2 r p α T p ) / k m ; Pr = μ c m / k m
where Cd is the drag coefficient, H p = T 0 T p c p ( T ) d T is the enthalpy of particles, rp is the particle radius, V p is the particle velocity, Tpm is the liquidus temperature of the SiC particles, and Hpm is the latent heat of fusion of SiC particles. The heat flux to an individual particle includes heat exchange with the liquid metal matrix.

3.6. Boundary Conditions

3.6.1. Energy Boundary Conditions

The energy boundary condition on the top surface of the substrate can be described by the heat input from the arc heat flux and the heat loss from convection and radiation, as shown in Equation (25).
k T n = q ˙ a r c q ˙ c o n v q ˙ r a d
where q a r c , q c o n v , q r a d are the arc heat input and convective and radiative heat loss, respectively, and n is the surface normal.
The molten pool created by the preceding variable polarity TIG welding arc is a prerequisite to ensure that a complete metallurgical bonding layer at the droplets/molten pool interface can be formed. The establishment of a reasonable heat flux field model of the variable polarity TIG welding arc is the key to obtaining accurate numerical calculation results. Based on an in-depth analysis of the TIG-assisted DDM, a dual ellipsoidal heat source model was used as the heat source model of TIG welding molten pool. The numerical simulation process considers the effect of the variable polarity pulse welding current on the flow pattern of the molten pool.
The arc heat input qarc is given as
q a r c = η U I
where η is the thermal efficiency of the arc, U is the arc voltage, I is the arc current.
The heat flux density qf on the front half axis of the double ellipsoid at (x, y, z) is
q f = 6 3 q a r c f f π a f b c exp ( 3 [ ( x v 0 t ) 2 a f 2 + y 2 b 2 + z 2 c 2 ] )
The heat flux density qr on the second half axis of the double ellipsoid at (x, y, z) can be calculated as
q r = 6 3 q a r c f r π a r b c exp ( 3 [ ( x v 0 t ) 2 a r 2 + y 2 b 2 + z 2 c 2 ] )
where af, ar, b and c are the distribution parameters of the arc heat source; ff and fr are the heat distribution coefficients of the two parts of the heat source, and v0 is the arc moving speed, t is time.
For single-track deposition SiC-particle reinforced AMCs (the welding torch is at an angle of 65° with the substrate surface), the distribution parameters of the dual ellipsoidal heat source model are DCEN stage: af = 4 mm, ar = 6 mm, b = 4 mm, c = 1.2 mm, DCEP stage: af = 3 mm, ar = 4 mm, b = 3 mm, c = 0.9 mm. The model distribution parameters under other process conditions are also determined according to the above method.
For other surfaces:
k T n = h A ( T T 0 ) σ ε r ( T 4 T 0 4 )
where hA is convection coefficient, σ and εr are Stefan–Boltzmann constant and radiation emissivity, respectively.

3.6.2. Momentum Boundary Conditions

For Marangoni convection, surface tension gradient on the free surface of the molten pool is equal to the shear stress:
μ u z = γ T T x μ v z = γ T T y
The normal pressure boundary conditions at the molten pool surface can be expressed as
p + 2 μ V n n = P a r c γ R c
where V n is the normal velocity, and Rc is radius of the surface curvature.

3.7. Numerical Method

Finite volume method was adopted to discrete the related governing equations, the structural finite difference grid was adopted to define the computational domain, and the grid size is 200 μm. Time step is set to be 5 × 10−6 s. The explicit method was applied to update the flow fields. The pressure equation (Poisson equation) is solved by Successive Over Relaxation (SOR) method to satisfy the continuity equation, and then, the energy equation was solved by the implicit method. Finally, the VOF equation was used to update the configuration of the free surface. These steps will be repeated at every time step until the specified simulation time is reached.

3.8. Properties of As-Used Materials and TIG-Assisted DDM Processing Parameters

The metal matrix feedstock material was a commercial ZL101 aluminum alloy; its chemical composition (mass fraction, %) is Si 6.6%, Mg 0.28%, Fe < 0.16%, Mn < 0.10%, Zn < 0.10%, and Al balanced. The liquidus and solidus temperature of SiC/ZL101 composites were tested by differential scanning calorimetric (DSC) method. The thermophysical properties of the as-used materials in TIG-assisted DDM are listed in Table 1.
The material of substrate was 2024 aluminum alloy (Al-Cu-Mg) in initial T6 condition with a thickness of 6 mm. The chemical composition of the substrate is 4.9Cu-1.8Mg-0.9Mn-0.5Si-0.5Fe-bal.Al (wt%). The material was cut into several pieces with 200 × 200 × 6 mm3 dimensions. The initial temperature of the substrate was controlled at 200 °C. A Fronius Magic Wave 3000 welding power source (Fronius International GmbH, Pettenbach, Austria) was adopted. The VP-TIG welding process used was variable-balance square-wave AC mode. Before deposition, oil and other impurities that sit on the substrate surface were removed using acetone. Pure argon (99%) was used as a shielding gas with flow rates of 15~18 L/min for VP-TIG, the working distance is about 6 mm. Initial conditions and the process parameters of TIG-assisted DDM are shown in Table 2. Solidification contact angle is assumed to be constant.

4. Experiment Details

4.1. TIG-Assisted DDM PROCESS of SiC-Particle Reinforced AMCs

In this study, the experimental material was 4 vol pct SiC reinforced ZL101 (Al-7Si-0.3Mg) composite, a customized TIG-assisted DDM machine was used. The particle size range is between 30 and 72 μm, with an average particle size of 45 μm, obtained using a Malvern Mastersizer 3000E (Malvern, UK). Before TIG-assisted DDM experiment, as-cast SiC/ZL101 composites were machined into feedstocks with mass of 1 kg. The Al-based MMCs was melted in an induction-heated crucible at 720~730 °C and then incubated for 15 min. The AMCs melt was ejected in the form of droplets from the graphite nozzle via piezoelectric ceramic rod. The falling droplets of SiC/ZL101 composites will continuously impact on the molten pool surface, and thus, droplets are arrested, and their final solid shape determined. A CMOS Near Infrared (NIR) camera (XIRIS XVC-1000e) was used to record high-resolution images of the morphology evolution of molten pool surface during droplet impact and the bright welding arc with an exposure time of 80 μs, a frame rate of 30 fps, and a spectral range up to 1500 nm. To reduce the risk of forming adhesions between the impacting AMCs droplets and the tip of the electrode, the location of the droplet impact and the center of molten pool in the horizontal direction was kept at a distance of 2 mm.

4.2. Microstructural Characterization

Metallographic examinations were conducted, samples for metallographic observation were cut, ground, and polished according to standard procedures. A solution consisting of HF (2 mL), HCl (3 mL), HNO3 (5 mL), and distilled water (180 mL) was taken as an etching agent, with an etching time of 15 s. Microstructural analysis of the dispersions of SiC reinforcing-particles in the ZL101 matrix was performed using a light optical microscope (LOM) (Olympus-BX60M, Tokyo, Japan).

5. Results and Discussion

5.1. Experiment Verification

The surface morphology and cross-section of the as-deposited sample are employed to validate the simulation results. Figure 2a shows the surface morphology of two single-track deposited layers with SiC reinforcing particles, they are formed under the same process conditions as listed in Table 2. The deposition direction is from right to left because the droplet impact position is fixed; the substrate tends to move from right to left.
The distributions of the reinforcement particles in the matrix were characterized by using scanning electron microscopy (SEM) equipped with energy dispersive spectrometry (EDS). SiC particles with irregular shapes is no-uniformly distributed in the Al alloy matrix (Figure 3a). Meanwhile, it was observed that there are no agglomerations of SiC particles and free regions in the Al matrix. The interface between the SiC particles and the Al matrix is relatively clear. According to the EDS line scanning results in Figure 3c, light signs of interface reaction were detected.
Table 3 gives the detailed feature sizes (i.e., width, height and the depth of penetration) of the calculated and experimental results from metallographic samples. From Figure 2b, it is clear that the simulated feature sizes of the single-track deposition geometry are in good agreement with experimental results. The relative errors are within 9.5%. A qualitative comparison between the experimental and simulation results for the final distribution state of SiC reinforcing particles on the cross section of the deposited layer is also presented in this figure. Clearly, both the quantitative and qualitative comparisons show that the established numerical model exhibits high accuracy and can be used to reveal the complex physics behind the droplet impact effects on molten pool dynamics and particulate migration behaviors.

5.2. Impact, Spreading, and Solidification Behaviors in the TIG-Assisted DDM Process

Figure 4 shows the solid fraction distribution and the surface deformation of a single-track deposited layer in the early stage of TIG-assisted DDM, which are associated with the impact, spreading, and solidification process.
As shown in Figure 4a, a shallow molten pool on the top surface of the moving substrate was firstly created by the fixed preceding TIG welding arc. The surface of molten pool was deformed under the action of arc forces and surface tension force, an obvious depression zone at the front of the molten pool can be seen, the melted metal will flow backward to form a swelling at the rear region of molten pool. As time increased, the successive molten aluminum alloy droplets with SiC reinforcing particles were released from specified heights, H, with a vertical velocity of 1.2 m/s. When a molten aluminum alloy droplet impacts on the rear region of molten pool, it coalesces firstly with the receiving liquid metal with particle reinforcement, as shown in Figure 4b. Then, most of the melt will rapidly spread on high-temperature solid surface around the molten pool. As shown in Figure 4c, a liquid crown will be visible above the surface. As the deposition process continues, the solidified microstructure in the deposition area is continuously moved from right to left under the action of the horizontally moving substrate (Figure 4d). The introduction of droplet energy leads to an increase in local heat accumulation, which in turn increases the size of the molten pool, the corresponding results are shown in Figure 4e–h.

5.3. Migration Behaviors of the SiC Particles in the Molten Pool

During the deposition process, there are complex two-way coupling drag forces between the SiC particle-reinforced phase and the molten aluminum matrix. At the same time, the annexation and repulsion of the interface between the reinforced phase and the solidified phase of SiC particles directly affect the distribution characteristics of the enhanced phase in the solidified tissue. Figure 5 shows the distribution position of the SiC particle-reinforced phase in the deposited layer.
It can be seen from Figure 5 that in the solidified microstructure of the deposited layer, the SiC particle-reinforced phase is mainly concentrated on both sides of the sedimentary layer due to the flow of the molten pool and the drag force of the melt. It is consistent with the experimental results. It is proved that the molten pool thermal flow field calculation model established in this paper can provide a theoretical basis for the quantitative study of heat and mass transport and solidification behavior in the TIG-assisted DDM of SiC particle-reinforced AMCs components.
Figure 6 demonstrates the temperature evolution of four fixed points during one deposition period. Here, the vertical distances from P1, P2, P3, P4 to the top surface of substrate are all kept at 2 mm, the coordinates of the four points are P1(−6, 0, 2), P2(−4.5, 0, 2), P3(−3, 0, 2) and P4(0, 0, 2). Point P4 is located below the arc action position.
It was clear that the magnitude of temperature fluctuation at the specified points was found to be closely related to the molten pool behavior induced by the droplet impact. During coalescence, after droplet impact the force is more towards the crater center than the crater sidewalls, and thus, the thermal energy of the impacting droplet is absorbed and dissipated by the molten. On the other hand, the temperature at the four points tend to decrease when the incoming droplet start to contact the molten pool. Therefore, a decrease in the temperature at the points can be seen. However, the magnitude of temperature fluctuation at point P3 are the more obvious, which can be attributed to the flow oscillation and the thermal energy variation induced by the impacting droplet and the variable polarity gas tungsten arc welding power.

6. Conclusions

In this study, a 3D numerical model based on the finite volume method was established to clarify the mechanism of the particulate migration and the molten pool behaviors induced by successive droplet impact and variable-polarity TIG welding arc. The resultant surface morphology of the AMCs deposited layer and the distribution state of SiC particle reinforcement were experimentally acquired, which are in a good agreement with the results predicted by simulation. The main conclusions can be summarized as follows:
(1)
The feasibility to fabricate ceramic-reinforced aluminum matrix composite components by TIG-assisted DDM has been proven.
(2)
The cross-sectional morphology of a single-track deposited layer and the distribution characteristics of the SiC particle-reinforced phase reasonably agree with the experimental results, and the relative errors are within 9.5%, which indicated that the model can reveal the formation mechanism of TIG-assisted DDM.
(3)
The non-uniformly distribution of SiC reinforcing particles on the cross-section of deposited layer could be found, which can be attributed to the viscous dissipation effect and the time-delayed velocity difference during the droplet spreading and recoiling stage.

Author Contributions

Conceptualization, J.D. and Y.Z.; methodology, J.D. and Y.Z.; validation, J.D.; Resources, Y.Z.; formal analysis, H.W., H.L. and X.Z.; data curation, Y.Z.; writing—original draft, J.D.; writing—review and editing, Y.Z.; supervision, J.D.; project administration, J.D. and Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The research is financially supported by the Ningbo Natural Science Foundation (202003N4071) and Zhejiang Provincial Natural Science Foundation (LQ21E010004).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the Corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chawla, N.; Chawla, K.K. Metal Matrix Composites; Spinger: New York, NY, USA, 2013. [Google Scholar]
  2. Schmidt, A.; Siebeck, S.; Götze, U.; Wagner, G.; Nestler, D. Particle-Reinforced Aluminum Matrix Composites (AMCs)—Selected Results of an Integrated Technology, User, and Market Analysis and Forecast. Metals 2018, 8, 143. [Google Scholar] [CrossRef] [Green Version]
  3. Behera, M.P.; Dougherty, T.; Singamneni, S. Conventional and Additive Manufacturing with Metal Matrix Composites: A Perspective. Procedia Manuf. 2019, 30, 159–166. [Google Scholar] [CrossRef]
  4. Islam, M.; Khalid, Y.; Ahmad, I.; Almajid, A.A.; Achour, A.; Dunn, T.J.; Akram, A.; Anwar, S. Microstructural Evaluation of Inductively Sintered Aluminum Matrix Nanocomposites Reinforced with Silicon Carbide and/or Graphene Nanoplatelets for Tribological Applications. Metall. Mater Trans. A 2018, 49, 2963–2976. [Google Scholar] [CrossRef]
  5. Kanth, U.R.; Rao, P.S.; Krishna, M.G. Mechanical behaviour of fly ash/SiC particles reinforced Al-Zn alloy-based metal matrix composites fabricated by stir casting method. J. Mater. Res. Technol. 2019, 8, 737–744. [Google Scholar] [CrossRef]
  6. Srinivasan, R.; Hemanth Shrinivasan, B.; Jeevan Prasath, K.; Johnson Saleth, R.; Anandhan, R.D. Experimental investigation of aluminium hybrid metal matrix composites processed through squeeze casting process. Mater. Today Proc. 2020, 27, 1821–1826. [Google Scholar] [CrossRef]
  7. Abdizadeh, H.; Baghchesara, M.A. Investigation on mechanical properties and fracture behavior of A356 aluminum alloy based ZrO2 particle reinforced metal-matrix composites. Ceram. Int. 2013, 39, 2045–2050. [Google Scholar] [CrossRef]
  8. Chawla, N.; Shen, Y.L. Mechanical Behavior of Particle Reinforced Metal Matrix Composites. Adv. Eng. Mater. 2001, 3, 357–370. [Google Scholar]
  9. Sharifi, E.M.; Karimzadeh, F. Wear behavior of aluminum matrix hybrid nanocomposites fabricated by powder metallurgy. Wear 2011, 271, 1072–1079. [Google Scholar] [CrossRef]
  10. Manfredi, D.; Calignano, F.; Krishnan, M.; Canali, R.; Ambrosio, E.P.; Biamino, S.; Ugues, D.; Pavese, M.; Fino, P. Additive Manufacturing of Al Alloys and Aluminium Matrix Composites (AMCs); Monteiro, W.A., Ed.; Light Metal Alloys Applications; IntechOpen: London, UK, 2014. [Google Scholar]
  11. Simchi, A.; Godlinski, D. Densification and microstructural evolution during laser sintering of A356/SiC composite powders. J. Mater. Sci. 2011, 46, 1446–1454. [Google Scholar] [CrossRef]
  12. Ghosh, S.K.; Saha, P. Crack and wear behavior of SiC particulate reinforced aluminium based metal matrix composite fabricated by direct metal laser sintering process. Mater. Design. 2011, 32, 139–145. [Google Scholar] [CrossRef]
  13. Simchi, A.; Godlinski, D. Effect of SiC particles on the laser sintering of Al–7Si–0.3Mg alloy. Scr. Mater. 2008, 59, 199–202. [Google Scholar] [CrossRef]
  14. Dadbakhsh, S.; Mertens, R.; Hao, L.; Van Humbeeck, J.; Kruth, J.P. Selective Laser Melting to Manufacture “In Situ” Metal Matrix Composites: A Review. Adv. Eng. Mater. 2019, 21, 1801244. [Google Scholar] [CrossRef] [Green Version]
  15. Chang, F.; Gu, D.; Dai, D.; Yuan, P. Selective laser melting of in-situ Al4SiC4+SiC hybrid reinforced Al matrix composites: Influence of starting SiC particle size. Surf. Coat. Technol. 2015, 272, 15–24. [Google Scholar] [CrossRef]
  16. Louvis, E.; Fox, P.; Sutcliffe, C.J. Selective laser melting of aluminium components. J. Mater. Process. Tech. 2011, 211, 275–284. [Google Scholar] [CrossRef]
  17. Kumar, A.; DebRoy, T. Calculation of three-dimensional electromagnetic force field during arc welding. J. Appl. Phys. 2003, 94, 1267–1277. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the TIG-assisted DDM process of particulate reinforced AMCs.
Figure 1. Schematic diagram of the TIG-assisted DDM process of particulate reinforced AMCs.
Materials 14 02430 g001
Figure 2. Top surface morphology and cross-section of a single-track deposited layer: (a,c) experimental results, (b,d) simulated results.
Figure 2. Top surface morphology and cross-section of a single-track deposited layer: (a,c) experimental results, (b,d) simulated results.
Materials 14 02430 g002
Figure 3. SEM microstructures of ZL101 Al matrix composites in the specified region (Figure 2c): (a) AMC-SiC; (b) interface between SiC particles and the Al alloy matrix; (c) EDS line scanning of the SiC particle and (d) EDS spectrums corresponding to the locations in (c).
Figure 3. SEM microstructures of ZL101 Al matrix composites in the specified region (Figure 2c): (a) AMC-SiC; (b) interface between SiC particles and the Al alloy matrix; (c) EDS line scanning of the SiC particle and (d) EDS spectrums corresponding to the locations in (c).
Materials 14 02430 g003
Figure 4. Solid fraction distribution and surface deformation of a single-track deposited layer in the early stage of TIG-assisted DDM.
Figure 4. Solid fraction distribution and surface deformation of a single-track deposited layer in the early stage of TIG-assisted DDM.
Materials 14 02430 g004
Figure 5. Distribution position of SiC particle-reinforced phase in the deposited layer.
Figure 5. Distribution position of SiC particle-reinforced phase in the deposited layer.
Materials 14 02430 g005
Figure 6. Temperature evolution at specified points (P1, P2, P3 and P4) during one deposition period.
Figure 6. Temperature evolution at specified points (P1, P2, P3 and P4) during one deposition period.
Materials 14 02430 g006
Table 1. The relevant thermophysical properties of SiC particle and ZL101 aluminum alloy.
Table 1. The relevant thermophysical properties of SiC particle and ZL101 aluminum alloy.
Thermal PropertiesValue
Solid density of ZL101ρs = 2500 (kg·m−3)
Liquid density of ZL101ρl = 2420 (kg·m−3)
Solidus temperatureTs = 829 (K)
Liquidus temperatureTl = 905 (K)
Latent heat of fusionL = 4.29 × 105 (J·kg−1)
Latent heat of evaporationLev = 2.8 × 107 (J·kg−1)
Thermal conductivity of solid phaseks = 167 (W·m−1·K−1)
Thermal conductivity of liquid phasekl = 79.5 (W·m−1·K−1)
Specific heat capacity of solid phaseCs = 1194 (J·kg−1·K−1)
Specific heat capacity of liquid phaseCl = 1265 (J·kg−1·K−1)
Dynamic viscosityν = 0.008 (Pa·s)
Surface tensionγ0= 0.871 (N·m−1)
Temperature coefficient of surface tensionksur = −3.5 × 10−4 (N·m−1·K−1)
Thermal expansion coefficientρT = 1.5 × 10−4 (K−1)
Particle densityρp = 3100 (kg·m−3)
Particle diameterrp = 40 (m−6)
Particle fractionf = 0.04
Thermal conductivity of particlekp = 120 (W·m−1·K−1)
Specific heat capacity of particleCp = 750 (J·kg−1·K−1)
Liquidus temperature of particleTpm = 1955 (K)
Latent heat of fusion of particleHpm = 3.7 × 105 (J·kg−1)
Ambient temperatureT0 = 298 (K)
Magnetic permeabilityμ0 = 1.26 × 10−6 (H·m−1)
Convection heat transfer coefficienthc = 40 (W·m−2·K−1)
Radiation emissivityε = 0.4
Stefan-Boltzmann constantσ = 5.67 × 108 (W·m−2·K−4)
Table 2. Initial conditions and the process parameters of TIG-assisted DDM.
Table 2. Initial conditions and the process parameters of TIG-assisted DDM.
Process ParametersValue
EP ratio50 (%)
Current cyclet0 = 20 (ms)
Arc current in EN phaseIen = 100 (A)
Arc current in EP stageIep = 200 (A)
Arc voltage at EN stageUen = 14.3 (V)
Arc voltage at EP stageUep = 12.1 (V)
Arc heat absorptionηen = 0.7
ηep = 0.4
Substrate moving speedv = 7 (mm s−1)
Torch angleθ = 65 (°)
Droplet temperatureTd = 1020 (K)
Droplet diameterd0 = 3 (mm)
Droplet frequencyf = 10 (Hz)
Impact velocity of dropletsv0 = 1.2 (mm)
Table 3. Comparison of calculated and experimental feature sizes of a single-track deposited layer.
Table 3. Comparison of calculated and experimental feature sizes of a single-track deposited layer.
SimulationExperimentError
Deposition width10.6 mm11.2 ± 0.3 mm5.4%
Deposition height3.1 mm3.3 ± 0.2 mm6.1%
Depth of penetration1.2 mm1.1 ± 0.1 mm9.1%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Du, J.; Wang, H.; Li, H.; Zhao, X. Numerical Analysis of Particulate Migration Behavior within Molten Pool during TIG-Assisted Droplet Deposition Manufacturing of SiC Particle-Reinforced Aluminum Matrix Composites. Materials 2021, 14, 2430. https://doi.org/10.3390/ma14092430

AMA Style

Zhang Y, Du J, Wang H, Li H, Zhao X. Numerical Analysis of Particulate Migration Behavior within Molten Pool during TIG-Assisted Droplet Deposition Manufacturing of SiC Particle-Reinforced Aluminum Matrix Composites. Materials. 2021; 14(9):2430. https://doi.org/10.3390/ma14092430

Chicago/Turabian Style

Zhang, Yubin, Jun Du, Huafeng Wang, Hua Li, and Xiaoyun Zhao. 2021. "Numerical Analysis of Particulate Migration Behavior within Molten Pool during TIG-Assisted Droplet Deposition Manufacturing of SiC Particle-Reinforced Aluminum Matrix Composites" Materials 14, no. 9: 2430. https://doi.org/10.3390/ma14092430

APA Style

Zhang, Y., Du, J., Wang, H., Li, H., & Zhao, X. (2021). Numerical Analysis of Particulate Migration Behavior within Molten Pool during TIG-Assisted Droplet Deposition Manufacturing of SiC Particle-Reinforced Aluminum Matrix Composites. Materials, 14(9), 2430. https://doi.org/10.3390/ma14092430

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop