Next Article in Journal
Synergistic Effects of Pr6O11 and Co3O4 on Electrical and Microstructure Features of ZnO-BaTiO3 Varistor Ceramics
Previous Article in Journal
Effect of Wear Conditions, Parameters and Sliding Motions on Tribological Characteristics of Basalt and Glass Fibre Reinforced Epoxy Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of Keyhole and Weld Pool Behaviors in Ultrasonic-Assisted Plasma Arc Welding Process

MOE Key Lab for Liquid-Solid Structure Evolution and Materials Processing, Institute of Materials Joining, Shandong University, Jinan 250061, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(3), 703; https://doi.org/10.3390/ma14030703
Submission received: 31 December 2020 / Revised: 24 January 2021 / Accepted: 28 January 2021 / Published: 2 February 2021

Abstract

:
The acoustic radiation force driving the plasma jet and the ultrasound reflection at the plasma arc-weld pool interface are considered to modify the formulas of gas shear stress and plasma arc pressure on the anode surface in ultrasonic-assisted plasma arc welding (U-PAW). A transient model taking into account the dynamic changes of heat flux, gas shear stress, and arc pressure on the keyhole wall is developed. The keyhole and weld pool behaviors are numerically simulated to predict the heat transfer and fluid flow in the weld pool and dynamic keyhole evolution process. The model is experimentally validated. The simulation results show that the acoustic radiation force increases the plasma arc velocity, and then increases both the plasma arc pressure and the gas shear stress on the keyhole wall, so that the keyholing capability is enhanced in U-PAW.

1. Introduction

Plasma arc welding (PAW) can form a keyhole channel inside the weld pool along the plate thickness, and join materials of mid-thickness in a single pass [1]. The stability of the weld pool and keyhole is the prerequisite for obtaining a large penetration and good weld quality [2]. However, in the PAW process, the dynamic stability of the weld pool and keyhole is poor, and the welding process-parameter window is narrow, which restricts its wider applications in the industry [3,4]. To further improve the weld pool dynamic stability and the keyholing capability in PAW, Wu et al. developed a process variant of PAW, i.e., ultrasonic-assisted plasma arc welding (U-PAW) [5], which exerts ultrasonic vibration on the tungsten electrode, as schematically illustrated in Figure 1. Different from ultrasonic-assisted TIG welding and ultrasonic-assisted MIG welding [6,7], due to the special structure of plasma arc welding torch and the setback of the tungsten electrode, a U-PAW torch was designed. The ultrasonic vibration at the end of the tungsten electrode is directly transmitted into the plasma arc, leading to the enhancement of heat-pressure characteristics, therefore the keyholing capability of plasma arc can be improved [8,9].
To study the process mechanism of U-PAW, Wu et al. [5] conducted experiments under different welding process parameters, and found that by the application of ultrasonic vibration, the plasma arc pressure is increased, so that an open keyhole can be formed at a higher welding velocity or lower welding current in U-PAW. The current density on the anode surface in U-PAW is increased due to the further constricted arc column [8]. Wang et al. [9] used the controlled pulse current waveform in U-PAW and carried out the closed-loop experiments. They found that with the application of ultrasonic vibration, the peak level of the welding current for a specific steel plate thickness is decreased. The above experimental results show that the ultrasonic vibration improves the keyholing ability in U-PAW and broadens the welding process window. An experimental investigation can only provide the test data for specific welding conditions, while the analytical method and numerical simulation can reveal the underlying physical mechanism in the welding process. However, the analytical formulae for thermal conduction in welding was derived based on some unreasonable assumptions for simplification, and large errors exist when they are used to predict the temperature distribution near and in the weld pool [10]. On the other hand, although the analytical method can save the calculation time [11,12], it can only be used to calculate the temperature field but not the flow field in the molten pool. Therefore, it is essential to conduct numerical modeling and simulation of fluid flow and heat transfer in the molten pool in both PAW and U-PAW for a complete understanding of the U-PAW process mechanism.
Recently, many researchers have conducted a numerical analysis of weld pool and keyhole behaviors in PAW [13,14,15,16,17,18,19]. However, there is little such work on the newly developed U-PAW. Li et al. [20,21] developed a model to study the interaction between the ultrasonic vibration and plasma arc from macroscopic and microscopic aspects, but did not consider the influence of plasma arc on the weld pool. When the plasma jet collides the workpiece surface, it will produce a stagnation pressure, called the arc pressure. According to the interaction between the ultrasonic and plasma arcs, the author deduced the arc pressure model [22]. However, the arc pressure was determined by assuming that the acoustic energy was totally absorbed by the plasma arc.
In this study, a mathematical model for weld pool and keyhole behaviors in U-PAW is developed. The arc pressure and gas shear stress formulae are modified by considering the acoustic radiation force and the ultrasonic reflection at the plasma arc-weld pool interface. The effects of the arc pressure and gas shear stress on the fluid flow and keyhole behavior in the weld pool were analyzed quantitatively, and the numerical simulation results are verified experimentally.

2. Formula of Plasma Arc Pressure in U-PAW

An ultrasonic action on the plasma arc medium will produce the corresponding acoustic radiation force, which is directly proportional to the sound intensity [23]. Based on the acoustic impedance theory, the interface between the plasma arc and weld pool can be taken as a hard boundary. The total ultrasound reflection occurs at the interface between the plasma arc and weld pool. The difference of the acoustic radiation force at the end of the tungsten electrode and the anode surface, F ¯ , is written as follows [24]:
F ¯ = 2 I 0 S t c ( 1 e 2 α L )
where I 0 is the sound intensity, S t is the area of tungsten end, c is the sound speed, L is the arc length, and α is the attenuation coefficient of sound pressure in the plasma arc.
In the plasma arc, during time interval Δ t , the acoustic radiation force will do some work equal to F ¯ c Δ t , while the electric field will do some work equal to η i I U Δ t , where η i is the conversion efficiency [15], I is the welding current, and U is the arc voltage. According to the kinetic energy theorem, we obtain the following:
F ¯ c Δ t + η i I U Δ t = 1 2 m v 2 2 1 2 m v 1 2
m = ρ Ar Q plas Δ t
where m is equal to the mass of the plasma gas ejected from the nozzle exit within the time interval Δ t , v 1 is the velocity of the plasma gas at the nozzle exit, v 2 is the velocity of the plasma arc near the anode surface, ρ Ar is the plasma gas density, and Q plas is the plasma gas flow rate.
Since the vibration of the tungsten end is harmonic vibration, the amplitude of the sound pressure is p u = ρ Ar c v u , and the amplitude of the vibration velocity is v u = 2 A π f , where A is the vibration amplitude, and f is the frequency. Therefore,
I 0 = p u 2 v u 2 = 1 2 ( ρ A r c v u ) v u = 1 2 ρ A r c ( 2 π A f ) 2 = 2 π 2 A 2 f 2 c ρ A r
Combined with Equations (1)–(4), we get the following:
v 2 2 v 1 2 = 2 F ¯ c + 2 η i I U ρ A r Q p l a s
Based on Equations (1) and (4), we get:
F ¯ c = 2 I 0 S t ( 1 e 2 α L ) = 4 π 2 A 2 f 2 c ρ A r S t ( 1 e 2 α L )
v 1 = Q p l a s S n o z
where S n o z is the nozzle exit area.
Manipulating Equations (5)–(7), we can get:
v 2 = 8 c A 2 π 2 f 2 S t ( 1 e 2 α L ) Q p l a s + 2 η i I U ρ A r Q p l a s + ( Q p l a s S n o z ) 2
The plasma arc velocity v 2 produces pressure on the surface of the anode. The arc pressure is assumed in the Gaussian distribution on the anode surface. In addition, considering the action of its own magnetic field inside the plasma arc, a coefficient β M is introduced as follows:
β M = μ 0 I 2 4 π 2 1 A I
where μ 0 is the permeability in a vacuum, and A I is a constant.
With the modified plasma arc velocity v 2 described by Equation (8), the similar derivation procedure in [22] may be used for the next step, and not iterated here. The modified arc pressure formula is written as follows:
P ( x , y ) = 3 β M r n o z 2 r p 2 ρ A r exp ( 3 ( x 2 + y 2 ) r p 2 ) 8 ρ A r c A 2 π 2 f 2 S t ( 1 e 2 α L ) Q p l a s + 2 η i I U Q p l a s ρ A r S n o z 2 + ( Q p l a s S n o z ) 4
where r P is the radius of the plasma arc pressure, r n o z is the nozzle exit radius, and (x,y) is the coordinate away from the torch axis.
Inside the keyhole, the keyhole radius is gradually smaller from the top to the bottom surface of the workpiece, which means that the keyhole channel further shrinks the plasma arc [25]. Therefore, the plasma arc pressure inside the keyhole increases with the keyhole depth, and a coefficient χ p is introduced to describe this effect as follows:
χ p = 1 + Γ ( x , y ) Z L + H n w
where Γ ( x , y ) = z is the function of the keyhole wall, H n w is the distance between the nozzle exit and the top of the workpiece surface, and Z L is the thickness of the workpiece. Finally, the plasma arc pressure is written as follows:
P ( x , y , z ) = 3 χ p β M r n o z 2 r p 2 ρ A r exp ( 3 ( x 2 + y 2 ) r p 2 ) 8 ρ A r c A 2 π 2 f 2 S t ( 1 e 2 α L ) Q p l a s + 2 η i I U Q p l a s ρ A r S n o z 2 + ( Q p l a s S n o z ) 4
Equation (12) includes the contribution components of the plasma arc pressure in U-PAW, including the acoustic effect term, electric field term, and plasma gas flow rate term. The early studies only correlated the arc pressure with the square of the welding current [26,27]. Though Li et al. [28] and Lang et al. [29] proposed other factors such as density, the nozzle exit area, and plasma gas flow rate, their formulae were suitable for the conventional PAW. Equation (12) is for the U-PAW process. The plasma arc pressure was measured by the water-cooled copper anode method. The detailed measurement method may be found in [5]. According to the calculation results in the literature [30], the arc pressure drops to a very small value at r = 3 mm, so the plasma arc pressure action radius r P is taken as 3 mm. For the plasma arc pressure on the anode surface, χ p = 1 . Under the welding conditions listed in Table 1, the distribution of PAW and U-PAW arc pressure was calculated by Equation (12), as shown in Figure 2. Compared with the previous formula [22], the arc pressure distribution predicted by the modified formula is more consistent with the experimental measurements. It can be seen that the arc pressure in U-PAW is greater than that in PAW. Therefore, Equation (12) is more reasonable.

3. Gas Shear Stress Formula

In the literature, the maximum value of the wall gas shear stress produced by the vertical jet is expressed as follows [31]:
τ = ( ρ A r μ D ) 1 2 ( D L ) 2 v 2 1.5
where μ is the plasma arc viscosity, D is the nozzle diameter, L is the distance between the nozzle exit and the top of the workpiece surface, and v 2 is the plasma arc velocity.
The plasma arc velocity v 2 at the workpiece surface is described by Equation (8). The maximum value of the plasma arc gas shear stress is written as follows:
τ u = ( ρ A r μ D ) 1 2 ( D L ) 2 χ p 1.5 ( 8 ρ A r c A 2 π 2 f 2 S t ( 1 e 2 α L ) + 2 η i I U ρ A r Q p l a s + ( Q p l a s S n o z ) 2 ) 0.75
For PAW without the ultrasonic vibration, the values of A and F are zero, so that Equation (14) can also be used to calculate the gas shear stress in PAW. Meng et al. [32] showed that the plasma shear stress increases linearly near the center of the weld pool surface, and then decreases rapidly. According to the results of the gas shear stress described by Wu et al. [18], its dimensionless distribution function is as follows:
g ( r ) = { 0.18 r 4 0.95 r 3 + 1.37 r 2 + 0.0957 ( r 0.208   m m ) 1.896 exp ( 0.315 r ) + 0.022 ( r > 0.208   m m )
where r = x 2 + y 2 + z 2 .
Therefore, the plasma arc gas shear stress formula is as follows:
τ u ( x , y , z ) = τ u g ( r )
Under the conditions listed in Table 1, the calculated gas shear stress on the anode surface are shown in Figure 3. It can be seen that after the ultrasound application, the plasma arc velocity increases, and then the peak value of the gas shear stress on the workpiece surface increases, and the peak value of the gas shear stress is 2 mm away from the torch axis, which is different from the result in the literature [18,19]. In this paper, the nozzle exit radius is 3.2 mm, so the arc pressure and gas shear stress are small. In the variable polarity PAW, with a peak welding current of 150 A, the peak value of the gas shear stress on the workpiece surface was about 85 Pa [33]. Bai et al. [34] gave a peak gas shear stress of 80 Pa when the welding current was 169 A. Therefore, the calculated value of the gas shear stress is reasonable.
Equation (16) gives the shear stress on the workpiece surface. In PAW, the weld pool surface is concave and deformed, so that the gas shear stress distribution will change with the keyhole shape. According to the description in the literature [32], when the weld pool surface is concave and deformed, there is a following relation between the plasma gas shear stress vector τ , the normal vector n of the deformed weld pool surface, and the gas shear stress acting direction vector o = ( x , y , 0 ) before the deformation of the weld pool surface, as follows:
{ τ n = 0 | τ | = τ ( x , y , z ) τ = A o + B n
where A and B are arbitrary real numbers.

4. Model of the Weld Pool and Keyhole

Owing to the symmetry of the workpiece, half of the calculation domain is used in the molten pool simulation. As shown in Figure 4, the size of the calculation domain is 26 mm × 14 mm × 10 mm. The thickness of the workpiece is 4 mm, and the air layer is 2 mm above and below the plate. The grid is evenly divided with a size of 0.2 mm. In this model, the welding torch is stationary while the workpiece is moving. The Cartesian coordinate system is established by taking the intersection of the center axis line of the welding torch and the upper surface of the workpiece as the coordinate origin.
Due to the interaction between the ultrasound, weld pool, and plasma arc in the U-PAW process, the physical phenomena involved are very complicated. The following simplifications and assumptions have to be made in the current research: (1) In the U-PAW process, since the density and sound velocity in liquid metal are relatively high, the acoustic streaming driving force caused by the difference of the acoustic energy density can be neglected. Therefore, the influence of sound pressure on the molten metal flow and heat transfer is temporarily not considered; (2) the weld pool flow is incompressible and laminar. Thereby, the heat transfer and fluid flow in the molten pool can be described by following governing equations.
Continuity equation:
ρ t + ( ρ v ) = 0
where ρ is the density, t is the time, and v is the fluid velocity vector.
Momentum equation:
t ( ρ v ) + ( ρ v v ) = p + μ 2 v + F v
where p is the pressure in fluid, μ is viscosity, F v is the momentum source term.
Energy equation:
t ( ρ H ) + ( v ρ H ) = ( k T ) + Q v
where H is the enthalpy, k is the thermal conductivity, T is the temperature, and Q v is the source term of the energy equation.
VOF equation:
t ( ϕ ) + ( ϕ v s ) = 0
0 < ϕ < 1
According to the description in [15], Equation (21) cannot accurately describe the interface of the plasma arc and the molten pool, so that a function δ is introduced to describe the gas-liquid interface and describe the physical interface definition function δ 1 of the upper surface of the workpiece:
δ 1 = { 1 if   ϕ ε 1   and   | ϕ | ε 2   and   ϕ z ε 3   0 else  
The physical interface defining function δ 2 describing the upper and lower surfaces of the workpiece is as follows:
δ 2 = { 1 if   ϕ ε 1   and   | ϕ | ε 2 0 else  
where ε 1 , ε 2 , ε 3 are very small constants.
The heat source, surface tension, arc pressure, electromagnetic force, gas shear stress, buoyancy, and the resistance force in a mushy zone are taken as the boundary conditions of the weld pool. The treatment method of the heat source model and arc pressure at the plasma arc-weld pool interface has been described in detail in our previous paper [22].
Table 2 lists two study cases. The thermophysical properties of a 304 stainless steel can be obtained from our previous work [22]. The boundary conditions are written by the user defined function (UDF), and are brought into the source terms of energy and momentum conservation equations to be solved. The equations are discretized by the finite volume method (FVM). The governing equations and boundary conditions were solved by the numerical method.

5. Experimental Process

Figure 5 is a schematic diagram of the U-PAW welding process setup. A vision system was developed to observe the keyhole exit and keyholing time from the backside of the workpiece [35]. As shown in Figure 6a, no open keyhole was detected in PAW. For U-PAW, Figure 6b shows that an open keyhole was firstly monitored at an instant of 4.6 s.

6. Results and Discussion

Figure 7 and Figure 8 show the calculated keyhole dynamic evolution process in PAW and U-PAW, respectively. At the initial time, the peak heat flux is mainly distributed at the top part of the keyhole wall. With the increase of the keyhole depth, the consumed plasma arc heat gets larger, the heat flux at the bottom part of the keyhole is lower, and the heat flux at the bottom part of the keyhole is smaller than that at the top part. As shown in Figure 7e, at t = 4.01 s, the gas shear force, Marangoni shear force, and arc pressure on the keyhole wall reach equilibrium. Therefore, the keyhole depth does not change any more in PAW. However, in U-PAW, the larger arc pressure promotes the keyhole formation, and the bottom of the keyhole continuously transmits the plasma arc heat flow and arc pressure to the weld pool. As shown in Figure 8f, at t = 4.2 s, an open keyhole is formed.
Figure 9 and Figure 10 show the fluid flow and heat transfer in the weld pool during the dynamic keyhole evolution in PAW and U-PAW, respectively. The main driving forces of the molten metal flow include the Marangoni shear force caused by the surface tension gradient, gas shear stress, arc pressure, electromagnetic force, and buoyancy [18,19]. Figure 9 shows the temperature and fluid flow fields at the longitudinal cross-section in PAW. At an instant of 1.21 s, the workpiece is melted to form a shallow and little weld pool under the thermal action of the plasma arc. At this time, the surface of the weld pool is a little bit depressed. At 2.02 s, a bigger depression of the weld pool surface is caused by the plasma arc pressure, forming a shallow blind keyhole, and the molten metal flows to the rear of the weld pool. At t = 4.01 s, since the plasma arc heat and pressure on the keyhole wall reaches equilibrium, the keyhole depth reaches a value of about 2.5 mm, and does not change any more. Under this PAW condition, an open keyhole cannot be formed, but the bottom of the weld pool reaches the melting point. In Figure 9e, there are two flow circulations in the weld pool. Under the action of gas shear stress and Marangoni force, liquid metal flows from the central high temperature zone near the pool center to the low temperature zone at the rear part of the pool, forming a clockwise circulation at the rear-upper part of the weld pool. Under the action of the arc pressure and electromagnetic force, the liquid metal is driven to flow to the bottom of the weld pool, and a counterclockwise circulation is formed in most of the weld pool parts, which is consistent with the X-ray observation results in [36].
In U-PAW, due to the exerted ultrasonic vibration, the plasma jet velocity increases, and both the gas shear stress and arc pressure increase. Therefore, the ultrasonic-assisted plasma arc has a stronger keyholing capability. As shown in Figure 10, although the fluid flow trend in U-PAW is almost the same as that in PAW, an open keyhole is formed at an instant of 4.2 s. Figure 10f demonstrates that after the keyhole is formed, the liquid metal at the front keyhole wall is very thin, while there is more liquid metal on the rear wall of the keyhole. The computed results about whether the open keyhole is formed or not and the establishment time of the open keyhole are in agreement with the experimental results.
Figure 11 shows the comparison of the experimental and calculated weld cross section in U-PAW. It can be seen that the calculated weld profile agree with the experimentally measured one. Table 3 lists the detailed comparison between the numerical calculation and experimental results of the weld cross section and keyhole exit size in U-PAW. The calculated weld widths at the top and bottom surfaces are in good agreement with the measured ones, but there is still a deviation of the fusion line track along the workpiece thickness direction. Similarly, the calculated keyhole exit size is smaller than the measured value. In our previous research, the influence of the gas shear force on the weld pool flow was neglected [22]. Therefore, the weld width at the top side is 8 mm, which is narrower than the experimental results. Considering the gas shear force, the top side weld width is closer to the experimental results. In this paper, we consider the distribution of gas shear stress on the keyhole wall through a simplified mathematical model. In the future work, we will get the distribution of the gas shear stress on the anode surface by modeling the plasma arc. In a next step, we will further improve the model and improve the numerical calculation accuracy.

7. Conclusions

(1)
The experimental results show that under the same welding process parameters, an open keyhole can be formed after applying the ultrasonic vibration, which further shows that U-PAW improves the keyholing ability.
(2)
The acoustic radiation force is considered to modify the formula of the plasma arc pressure in the ultrasonic-assisted PAW, and the prediction accuracy of the plasma arc pressure on the anode surface in U-PAW is improved.
(3)
With the modified velocity of the plasma arc near the anode surface, a formula is proposed to calculate the gas shear stress on the workpiece surface.
(4)
The effects of ultrasonic vibration on the plasma arc pressure and gas shear stress are considered in modeling the weld pool and keyhole behaviors in U-PAW. The dynamic keyhole evolution behavior in the weld pool is quantitatively analyzed.
(5)
The numerical simulation results of the weld pool widths, establishment time of open keyhole, and keyhole exit sizes are in agreement with the experimental ones.

Author Contributions

J.Q., numerical simulation, formal analysis, validation, original draft; C.W., conceptualization, review and editing, supervision; Y.L., formal analysis, review and editing, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant no. 51775312).

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. Wu, C.S.; Wang, L.; Ren, W.J.; Zhang, X.Y.J. Plasma arc welding: Process, sensing, control and modeling. Manuf. Process. 2014, 16, 74–85. [Google Scholar] [CrossRef]
  2. Liu, Z.M.; Cui, S.L.; Luo, Z.; Zhang, C.Z.; Wang, Z.M.; Zhang, Y.C. Plasma arc welding: Process variants and its recent developments of sensing, controlling and modeling. J. Manuf. Process. 2016, 23, 315–327. [Google Scholar] [CrossRef]
  3. Jia, C.B.; Liu, X.F.; Wu, C.S.; Lin, S.B. Stereo analysis on the keyhole and weld pool behaviorsin K-PAW with triple CCD cameras. J. Manuf. Process. 2018, 32, 754–762. [Google Scholar] [CrossRef]
  4. Lu, Z.Y.; Zhang, W.; Jiang, F.; Chen, S.J.; Yan, Z.Y. A primary study of variable polarity plasma arc welding using a pulsed plasma gas. Materials 2019, 12, 1666. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Wu, C.S.; Zhao, C.Y.; Zhang, C.; Li, Y.F. Ultrasonic vibration assisted keyholing plasma arc welding. Weld J. 2017, 96, 279–286. [Google Scholar]
  6. Sun, Q.; Lin, S.; Yang, C.; Zhao, G. Penetration increase of AISI 304 using ultrasonic assisted tungsten inert gas welding. Sci. Technol. Weld. Join. 2009, 14, 765–767. [Google Scholar] [CrossRef]
  7. Fan, Y.; Yang, C.; Lin, S.; Fan, C.; Liu, W. Ultrasonic wave assisted GMAW. Weld J. 2012, 91, 91s–99s. [Google Scholar]
  8. Zhang, C.; Wu, C.S.; Tian, S.S. Effect of ultrasonic vibration on current density and keyholing capability of plasma arc. Sci. Technol. Weld. Join. 2020, 25, 422. [Google Scholar] [CrossRef]
  9. Wang, L.; Zhang, C.; Wu, C.S. Experimental study on controlled pulse keyholing plasma arc welding assisted by ultrasonic vibration. Int. J. Adv. Manuf. Technol. 2020, 107, 4995–5009. [Google Scholar] [CrossRef]
  10. Wu, C.S. Welding Thermal Processes and Weld Pool Behaviors; CRC Press/Taylor & Francis Group: Roca Raton, FL, USA, 2011. [Google Scholar]
  11. Mahmood, M.A.; Popescu, A.C.; Hapenciuc, C.L.; Ristoscu, C.; Visan, A.; Oane, M.; Mihailescu, I.N. Estimation of clad geometry and corresponding residual stress distribution in laser melting deposition: Analytical modeling and experimental correlations. Int. J. Adv. Manuf. Technol. 2020, 111, 77–91. [Google Scholar] [CrossRef]
  12. Mahmood, M.A.; Popescu, A.C.; Oane, M.; Ristoscu, C.; Chioibasu, D.; Mihai, S.; Mihailescu, I.N. Three-Jet Powder Flow and Laser–Powder Interaction in Laser Melting Deposition: Modelling Versus Experimental Correlations. Metals 2020, 10, 1113. [Google Scholar] [CrossRef]
  13. Li, Y.; Feng, Y.H.; Zhang, X.X.; Wu, C.S. An improved simulation of heat transfer and fluid flow in plasma arc welding with modified heat source model. Int. J. Therm. Sci. 2013, 64, 93–104. [Google Scholar] [CrossRef]
  14. Xu, G.; Hu, J.; Tsai, H.L. Three-dimensional modeling of the plasma arc in arc welding. J. Appl. Phys. 2008, 104, 103301. [Google Scholar] [CrossRef]
  15. Li, T.Q.; Wu, C.S. Numerical simulation of plasma arc welding with keyhole-dependent heat source and arc pressure distribution. Int. J. Adv. Manuf. Technol. 2015, 78, 593–602. [Google Scholar] [CrossRef]
  16. Li, T.Q.; Chen, L.; Zhang, Y.; Yang, X.M.; Lei, Y.C. Metal flow of weld pool and keyhole evolution in gas focusing plasma arc welding. Int. J. Heat Mass Transf. 2020, 150, 119296. [Google Scholar] [CrossRef]
  17. Li, Y.; Feng, Y.H.; Zhang, X.X.; Wu, C.S. An evolutionary keyhole-mode heat transfer model in continuous plasma arc welding. Int. J. Heat Mass Transf. 2018, 117, 1188–1198. [Google Scholar] [CrossRef]
  18. Wu, D.S.; Nguyen, A.V.; Tashiro, S.; Hua, X.M.; Tanaka, M. Elucidation of the weld pool convection and keyhole formation mechanism in the keyhole plasma arc welding. Int. J. Heat Mass Transf. 2019, 131, 920–931. [Google Scholar] [CrossRef]
  19. Wu, D.S.; Tashiro, S.; Hua, X.M.; Tanaka, M. A novel electrode-arc-weld pool model for studying the keyhole formation in the keyhole plasma arc welding process. J. Phys. D Appl. Phys. 2019, 52, 165203. [Google Scholar] [CrossRef]
  20. Li, Y.F.; Wu, C.S.; Chen, M. Effects of ultrasonic vibration on the transport coefficients in plasma arc welding. Metals 2020, 10, 312. [Google Scholar] [CrossRef] [Green Version]
  21. Li, Y.F.; Wu, C.S.; Chen, M.A. Numerical analysis of the ultrasound induced arc pressure increment in plasma arc welding. J. Phys. D Appl. Phys. 2019, 52, 1–14. [Google Scholar] [CrossRef]
  22. Qiao, J.N.; Wu, C.S.; Li, Y.F. Numerical and experimental investigation of keyholing process in ultrasonic vibration assisted plasma arc welding. J. Manuf. Process. 2020, 50, 603–613. [Google Scholar] [CrossRef]
  23. Francisco, J.T.K.; Kai, K. A computational modeling approach of the jet-like acoustic streaming and heat generation induced by low frequency high power ultrasonic horn reactors. Ultraso. Sonochem. 2011, 18, 1263–1273. [Google Scholar]
  24. David, C. Fundamentals and Applications of Ultrasonic Waves; CRC Press: Montrel, QC, Canada, 2002. [Google Scholar]
  25. Xu, B.; Jiang, F.; Chen, S.J.; Tanaka, M.; Tashiro, S.; Nguyen, A.V. Numerical analysis of plasma arc physical characteristics under additional constraint of keyhole. Chin. Phys. B 2018, 27, 034701. [Google Scholar] [CrossRef]
  26. Allum, C.J. Gas flow in the column of a TIG welding arc. J. Phys. D Appl. Phys. 1981, 14, 1041–1059. [Google Scholar] [CrossRef]
  27. Lin, M.L.; Eagar, T.W. Pressures produced by gas tungsten arcs. Metall. Mater. Trans. B 1986, 17, 601–607. [Google Scholar] [CrossRef]
  28. Li, T.Q.; Wu, C.S. An analytic formula describing the plasma arc pressure distribution. China Weld. 2014, 23, 7–11. [Google Scholar]
  29. Lang, R.Q.; Han, Y.Q.; Bai, X.Y.; Hong, H.T. Prediction of the weld pool stability by material flow behavior of the perforated weld pool. Materials 2020, 13, 303. [Google Scholar] [CrossRef] [Green Version]
  30. Li, Y.F.; Wu, C.S.; Chen, M.A. Numerical analysis of the heat-pressure characteristics in ultrasonic vibration assisted plasma arc. J. Appl. Phys. 2020, 128, 114903. [Google Scholar] [CrossRef]
  31. Phares, D.J.; Smedley, G.T.; Flagan, R.C. The wall shear stress produced by the normal impingement of a jet on a flat surface. J. Fluid. Mech. 2000, 418, 351–375. [Google Scholar] [CrossRef] [Green Version]
  32. Meng, X.M.; Qin, G.L.; Bai, X.Y.; Zou, Z.D. Numerical analysis of undercut defect mechanism in high speed gas tungsten arc welding. J. Mater. Process. Technol. 2016, 236, 225–234. [Google Scholar] [CrossRef]
  33. Pan, J.J.; Hu, S.S.; Yang, L.J.; Chen, S.J. Numerical analysis of the heat transfer and material flow during keyhole plasma arc welding using a fully coupled tungsten–plasma–anode model. Acta Mater. 2016, 118, 221–229. [Google Scholar] [CrossRef]
  34. Bai, X.W.; Colegrove, P.; Ding, J.L.; Zhou, X.M.; Diao, C.L.; Bridgeman, P.; Hönnige, J.R.; Zhang, H.O.; Williams, S. Numerical analysis of heat transfer and fluid flow in multilayer deposition of PAW-based wire and arc additive manufacturing. Int. J. Heat Mass Transf. 2018, 124, 504–516. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, G.K.; Wu, C.S.; Liu, X.F. Single vision system for simultaneous observation of keyhole and weld pool in plasma arc welding. J. Mater. Process. Technol. 2015, 215, 71–78. [Google Scholar] [CrossRef]
  36. Nguyen, A.V.; Tashiro, S.; Tanaka, M.; Hanh, B.V. Experimental investigation on weld pool formation process in plasma keyhole arc welding. J. Phys. D Appl. Phys. 2018, 51, 1–14. [Google Scholar]
Figure 1. Schematic illustration of the ultrasonic-assisted plasma arc welding process.
Figure 1. Schematic illustration of the ultrasonic-assisted plasma arc welding process.
Materials 14 00703 g001
Figure 2. Plasma arc pressure distribution.
Figure 2. Plasma arc pressure distribution.
Materials 14 00703 g002
Figure 3. Gas shear stress distribution (welding current 80 A).
Figure 3. Gas shear stress distribution (welding current 80 A).
Materials 14 00703 g003
Figure 4. Calculation domain.
Figure 4. Calculation domain.
Materials 14 00703 g004
Figure 5. Schematic illustration of the ultrasonic-assisted plasma arc welding (U-PAW) process setup.
Figure 5. Schematic illustration of the ultrasonic-assisted plasma arc welding (U-PAW) process setup.
Materials 14 00703 g005
Figure 6. The captured images of the keyhole exit (a) PAW (b) U-PAW.
Figure 6. The captured images of the keyhole exit (a) PAW (b) U-PAW.
Materials 14 00703 g006
Figure 7. The dynamic keyhole evolution process in PAW (a) 1.21, (b) 2.02, (c) 2.5, (d) 3.2, (e) 4.01, (f) 4.6 s.
Figure 7. The dynamic keyhole evolution process in PAW (a) 1.21, (b) 2.02, (c) 2.5, (d) 3.2, (e) 4.01, (f) 4.6 s.
Materials 14 00703 g007
Figure 8. The dynamic keyhole evolution process in U-PAW (a) 1, (b) 2.01, (c) 2.6, (d) 3, (e) 3.31, (f) 4.2 s.
Figure 8. The dynamic keyhole evolution process in U-PAW (a) 1, (b) 2.01, (c) 2.6, (d) 3, (e) 3.31, (f) 4.2 s.
Materials 14 00703 g008
Figure 9. The calculated temperature field, fluid flow, and keyholing time in the PAW case.
Figure 9. The calculated temperature field, fluid flow, and keyholing time in the PAW case.
Materials 14 00703 g009
Figure 10. The calculated temperature field, fluid flow, and keyholing time in the U-PAW case.
Figure 10. The calculated temperature field, fluid flow, and keyholing time in the U-PAW case.
Materials 14 00703 g010
Figure 11. Weld cross section of U-PAW (a) measured (b) calculated.
Figure 11. Weld cross section of U-PAW (a) measured (b) calculated.
Materials 14 00703 g011
Table 1. The used welding parameters.
Table 1. The used welding parameters.
ParametersValue
Density of Ar (kg/m3)1.784
Adjusting constant (N)1.256 × 10−6
Plasma gas flow rate (L/min)2.8
Conversion coefficient2.04 × 10−4
Nozzle diameter (mm)3.2
Distance from electrode to workpiece (m)0.005
Shielding gas flow rate (L/min)20
Sound velocity of argon (m/s)341
Adjusting constant (N)4.58 × 10−4
Vibration frequency (KHz)25
Vibration amplitude (μm)20
Nozzle exit area (m2)8.0384 × 10−6
Table 2. Simulation cases.
Table 2. Simulation cases.
Welding Current (A)Welding Speed (mm/min)Ultrasonic VibrationWorkpiece Thickness (mm)
PAW80100No4
U-PAW80100Yes4
Table 3. Comparison of the experimental and calculated results.
Table 3. Comparison of the experimental and calculated results.
ItemMeasuredCalculated
Keyholing time (s)4.64.2
Keyhole exit length (mm)0.980.9
Keyhole exit width (mm)1.20.8
Weld width (top) (mm)8.88.9
Weld width (bottom) (mm)2.83
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qiao, J.; Wu, C.; Li, Y. Numerical Analysis of Keyhole and Weld Pool Behaviors in Ultrasonic-Assisted Plasma Arc Welding Process. Materials 2021, 14, 703. https://doi.org/10.3390/ma14030703

AMA Style

Qiao J, Wu C, Li Y. Numerical Analysis of Keyhole and Weld Pool Behaviors in Ultrasonic-Assisted Plasma Arc Welding Process. Materials. 2021; 14(3):703. https://doi.org/10.3390/ma14030703

Chicago/Turabian Style

Qiao, Junnan, Chuansong Wu, and Yongfeng Li. 2021. "Numerical Analysis of Keyhole and Weld Pool Behaviors in Ultrasonic-Assisted Plasma Arc Welding Process" Materials 14, no. 3: 703. https://doi.org/10.3390/ma14030703

APA Style

Qiao, J., Wu, C., & Li, Y. (2021). Numerical Analysis of Keyhole and Weld Pool Behaviors in Ultrasonic-Assisted Plasma Arc Welding Process. Materials, 14(3), 703. https://doi.org/10.3390/ma14030703

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