Mathematical modeling of three-dimensional heat and fluid flow in a moving gas metal arc weld pool

Mathematical models capable of accurate prediction of the weld bead and weld pool geometry in gas metal arc (GMA) welding processes would be valuable for rapid development of welding procedures and empirical equations for control algorithms in automated welding applications. This article introduces a three-dimensional (3-D) model for heat and fluid flow in a moving GMA weld pool. The model takes the mass, momentum, and heat transfer of filler metal droplets into consideration and quantitatively analyzes their effects on the weld bead shape and weld pool geometry. The algorithm for calculating the weld reinforcement and weld pool surface deformation has been proved to be effective. Difficulties associated with the irregular shape of the weld bead and weld pool surface have been successfully overcome by adopting a boundary-fitted nonorthogonal coordinate system. It is found that the size and profile of the weld pool are strongly influenced by the volume of molten wire, impact of droplets, and heat content of droplets. Good agreement is demonstrated between predicted weld dimensions and experimently measured ones for bead-on-plate GMA welds on mild steel plate.     

A New Model for Keyhole Mode Laser Welding Using FLUENT

Evolution of the free surface at gas?Cliquid interface during keyhole mode welding is complex and its calculation is computationally expensive. Similarly, models based on only heat conduction without considering vapour cavity and liquid convection around it, are computationally efficient but are not effective in defining the weld pool shape especially for low conducting material like steel. In the present study a useful yet computationally efficient model has been presented for keyhole mode laser welding using commercial software FLUENT. Here instead of evolving the free surface of keyhole in a rigorous way, various possible steady keyhole shapes are assumed partially based on literature evidence and subsequently their dimensions are calculated by an overall heat balance. The estimated keyhole profile is then mapped into the thermo-fluid framework of FLUENT and steady computational fluid dynamics calculations is carried out around the keyhole that is considered rigid wall at boiling temperature. Next, an optimized keyhole shape is identified by comparing the predicted fusion lines with the experimental weld fusion lines reported in literature. Finally, using this optimized keyhole shape independent predictions are made for two materials of widely different thermal conductivities, like steel and aluminum, under different operating conditions. In all cases the results of the present simulation is found to in close agreement with experimental data and even better than the model predictions reported in literature. The present model emerges as a simple yet effective model for predicting the weld bead profile encompassing wide range of materials under different operating conditions.     

Prediction of weld pool and reinforcement dimensions of GMA welds using a finite-element model

Mathematical models of the gas metal arc (GMA) welding process may be used to study the influence of various welding parameters on weld dimensions, to assist in the development of welding procedures, and to aid in the generation of process control algorithms for automated applications. In this work, a three-dimensional (3-D), steady-state thermal model of the GMA welding process has been formulated for a moving coordinate framework and solved using the finite-element method. The model includes temperature-dependent material properties, a new finite-element formulation for the inclusion of latent heat of fusion, a Gaussian distribution of heat flux from the arc, plus the effects of mass convection into the weld pool from the melted filler wire. The influence of weld pool convection on the pool shape was approximated using anisotropically enhanced thermal conductivity for the liquid phase. Weld bead width and reinforcement height were predicted using a unique iterative technique developed for this purpose. In this paper, the numerical model is shown to be capable of predicting GMA weld dimensions for individual welds, including those with finger penetration. Also, good agreement is demonstrated between predicted weld dimensions and experimentally derived relations that describe the effects of process variables and their influence on average weld dimensions for bead-onplate GMA welds on steel plate. E. PARDO, formerly Postdoctoral Fellow, Department of Mechanical Engineering, University of Waterloo, Waterloo, ON, Canada N2L 3G1,     

Effect of Welding Speed on Gas Metal Arc Weld Pool in Commercially Pure Aluminum: Theoretically and Experimentally

Temperature and velocity fields, and weld pool geometry during gas metal arc welding (GMAW) of commercially pure aluminum were predicted by solving equations of conservation of mass, energy and momentum in a three-dimensional transient model. Influence of welding speed was studied. In order to validate the model, welding experiments were conducted under the similar conditions. The calculated geometry of the weld pool were in good agreement with the corresponding experimental results. It was found that an increase in the welding speed results in a decrease peak temperature and maximum velocity in the weld pool, weld pool dimensions and width of the heat-affected zone (HAZ). Dimensionless analyses were employed to understand the importance of heat transfer by convection and the roles of various driving forces in the weld pool. According to dimensionless analyses droplet driving force strongly affected fluid flow in the weld pool.     

A mathematical model of gas tungsten arc welding considering the cathode and the free surface of the weld pool

A two-dimensional axisymmetric numerical model, including the influence of the cathode and the free surface of the weld pool, is developed to describe the heat transfer and fluid flow in gas tungsten arc (GTA) welding. In the model, a boundary-fitted coordinate system is adopted to precisely describe the cathode shape and deformed weld-pool surface. The current continuity equation has been solved with the combined arc plasma-cathode system, independent of the assumption of current density distribution on the cathode surface, which was essential in the previous studies of arc plasma. It has been shown that the temperature profile, the current, and the heat flux to the anode show good agreement with the experimental data. Moreover, the current and the heat-flux distributions may be affected by the shape of the cathode and the free surface of the weld pool.     

Computational modeling of laser welding of Cu-Ni dissimilar couple

A three-dimensional transient model to solve heat transfer, fluid flow, and species conservation during laser welding of dissimilar metals is presented. The model is based on a control volume formulation with an enthalpy-porosity technique to handle phase change and a mixture model to simulate mixing of molten metals. Weld pool development, solidified weld pool shape, and composition profiles are presented for both stationary as well as continuous laser welding in conduction mode. Salient features of a dissimilar Cu-Ni weld are summarized and thermal transport arguments are employed to successfully explain the observations. It is found that the weld pool shape becomes asymmetric even when the heat source is symmetrically applied on the two metals forming the couple. It is also observed that convection plays an important role in the development of weld pool shape and composition profiles. As the weld pool develops, the side melting first (nickel) is found to experience more convection and better mixing. Results from the case studies of computation are compared with corresponding experimental observations, showing good qualitative agreement between the two.     

Marangoni convection in weld pool in CO<Subscript>2</Subscript>-Ar-shielded gas thermal arc welding

Small CO₂ additions of 0.092 to 10 vol pct to the Ar shielding gas dramatically change the weld shape and penetration from a shallow flat-bottomed shape, to a deep cylindrical shape, to a shallow concave-bottomed shape, and back to the shallow flat-bottomed shape again with increasing CO₂ additions in gas thermal arc (GTA) welding of a SUS304 plate. Oxygen from the decomposition of CO₂ transfers and becomes an active solute element in the weld pool and reverses the Marangoni convection mode. An inward Marangoni convection in the weld pool occurs when the oxygen content in the weld pool is over 80 ppm. Lower than 80 ppm, flow will change to the outward direction. An oxide layer forms on the weld pool in the welding process. The heavy oxide layer on the liquid-pool surface will inhibit the inward fluid flow under it and also affects the oxygen transfer to the liquid pool. A model is proposed to illustrate the interaction between the CO₂ gas and the molten pool in the welding process.     

Three-dimensional transient model for arc welding process

A direct computer simulation technique, discrete element analysis (DEA), was utilized in the development of a transient multidimensional (2-D and 3-D) mathematical model for investi-gating coupled conduction and convection heat transfer problems associated with stationary and moving arc welding processes. The mathematical formulation considers buoyancy, electro-magnetic, and surface tension driving forces in the solution of the overall heat transfer conditions in the specimen. Furthermore, the formulation of the model allows realistic consideration of the geometrical variations in the workpiece. The model treats the -weld pool surface as a truly deformable free surface, allowing for the prediction of the weld surface deformations such as the “weld crown.≓ A marked element formulation was employed to monitor the transient de-velopment of the weld pool as determined by the latent heat considerations and the calculated velocities in the weld pool. The model was utilized to simulate the heat and fluid flows in the weld pool that occur during stationary (spot) and moving (linear) gas tungsten-arc welding. Also, the present analysis considers a simple rectangular specimen and a geometrically complex specimen to demonstrate the capability of the model to simulate realistic 3-D arc welding prob-lems. The results of the present investigation clearly demonstrate the significant influence of the heat and fluid flows and the specimen geometry on the development of the weld. Comparison of the predicted and the experimentally observed fusion zone and heat-affected zone (HAZ) geometries indicate good agreement.     

On the calculation of the free surface temperature of gas-tungsten-arc weld pools from first principles: Part II. modeling the weld pool and comparison with experiments

By combining a mathematical model of the welding arc and of the weld pool, calculations are presented to describe the free surface temperature of weld pools for spot welding operations. The novel aspects of the treatment include the calculation of the heat and current fluxes falling on the free weld pool surface from first principles, a realistic allowance for heat losses due to vaporization, and a realistic allowance for the temperature dependence of the surface tension. The most important finding reported in this article is that the free surface temperature of weld pools appears to be limited by Marangoni convection, rather than heat losses due to vaporiza-tion. Furthermore, it was found that once thermocapillary flow can produce high enough surface velocities (>25 cm/s), the precise nature of the relationship between temperature and surface tension will become less important.     

Investigation on the effect of laser pulse shape during Nd:YAG laser microwelding of thin Al sheet by numerical simulation

In the welding of thin A3003 Al sheet by a Nd:YAG laser beam, the laser pulse shape plays an important role in enhancing the welding penetration stability. In order to evaluate the effect of laser pulse shape during Nd:YAG laser welding of a thin Al sheet and to predict the welding performance by numerical simulation, a three-dimensional finite differential method (FDM) analysis is presented for heating with different laser pulse shapes and related welding parameters. The simulated results give good agreement with experimental results, where a sound weld shape and crack-free weld pool are obtained. The simulation results show that the welding stability is greatly affected by the modulation of laser pulse shape for the same laser energy and welding parameters. As a rectangular laser pulse is modulated to have three stages with high, medium, and low power levels for the first, second, and third stages, respectively, more energy is absorbed in the melt pool, and the cooling rate is reduced. While a high power level at the first stage increases the laser beam absorption, the thermal energy of the third stage prevents fast cooling of the melt pool. Also, evaporation is prevented by proper modulation of the laser pulse. If the laser pulse is modulated properly, the optimum melt-pool size and cooling rate can be obtained; also, the desired weld depth and welding stability are achieved for the conduction welding mode. The numerical simulation method can be used to determine the proper conditions for good welding performance.     

A unified numerical modeling of stationary tungsten-inert-gas welding process

In order to clarify the formative mechanism of weld penetration in an arc welding process, the development of a numerical model of the process is quite useful for understanding quantitative values of the balances of mass, energy, and force in the welding phenomena because there is still lack of experimentally understanding of the quantitative values of them because of the existence of complicated interactive phenomena between the arc plasma and the weld pool. The present article is focused on a stationary tungsten-inert-gas (TIG) welding process for simplification, but the whole region of TIG arc welding, namely, tungsten cathode, arc plasma, workpiece, and weld pool is treated in a unified numerical model, taking into account the close interaction between the arc plasma and the weld pool. Calculations in a steady state are made for stationary TIG welding in an argon atmosphere at a current of 150 A. The anode is assumed to be a stainless steel, SUS304, with its negative temperature coefficient of surface tension. The two-dimensional distributions of temperature and velocity in the whole region of TIG welding process are predicted. The weld-penetration geometry is also predicted. Furthermore, quantitative values of the energy balance for the various plasma and electrode regions are given. The predicted temperatures of the arc plasma and the tungsten-cathode surface are in good agreement with the experiments. There is also approximate agreement of the weld shape with experiment, although there is a difference between the calculated and experimental volumes of the weld. The calculated convective flow in the weld pool is mainly dominated by the drag force of the cathode jet and the Marangoni force as compared with the other two driving forces, namely, the buoyancy force and the electromagnetic force.     

Modeling macro-and microstructures of Gas-Metal-Arc Welded HSLA-100 steel

Fluid flow and heat transfer during gas-metal-arc welding (GMAW) of HSLA-100 steel were studied using a transient, three-dimensional, turbulent heat transfer and fluid flow model. The temperature and velocity fields, cooling rates, and shape and size of the fusion and heat-affected zones (HAZs) were calculated. A continuous-cooling-transformation (CCT) diagram was computed to aid in the understanding of the observed weld metal microstructure. The computed results demonstrate that the dissipation of heat and momentum in the weld pool is significantly aided by turbulence, thus suggesting that previous modeling results based on laminar flow need to be re-examined. A comparison of the calculated fusion and HAZ geometries with their corresponding measured values showed good agreement. Furthermore, “finger” penetration, a unique geometric characteristic of gas-metal-arc weld pools, could be satisfactorily predicted from the model. The ability to predict these geometric variables and the agreement between the calculated and the measured cooling rates indicate the appropriateness of using a turbulence model for accurate calculations. The microstructure of the weld metal consisted mainly of acicular ferrite with small amounts of bainite. At high heat inputs, small amounts of allotriomorphic and Widmanstätten ferrite were also observed. The observed microstructures are consistent with those expected from the computed CCT diagram and the cooling rates. The results presented here demonstrate significant promise for understanding both macro-and microstructures of steel welds from the combination of the fundamental principles from both transport phenomena and phase transformation theory.     

Simulation of Heat Flow During the Welding of Thin Plates

Based on the finite difference method and the enthalpy model of Shamsundar, a computer model was developed to describe the steady state, two-dimensional heat flow during the welding of thin plates. In order to allow accurate computations of the weld pool configuration, the size of the mushy zone and the temperature distribution near the heat source, a grid mesh of variable spacings was used. The heat of fusion, the size and distribution of the heat source, the temperature dependence of thermal properties, the heat conduction in the welding direction and the surface heat loss during welding were considered. The model was first checked with Rosenthal’s analytical solution of welding heat flow using pure aluminum for examples. Experimental results of 6061 aluminum, including the width of the fusion zone and the thermal cycles at positions in both the fusion and the heat affected zones, were then compared with the calculated results of the heat flow model. The agreement was very good. Finally, in order to demonstrate systematically the quantitative effect of welding parameters such as the heat input, the welding speed and the preheating of the workpiece, a series of computations were made based upon 6061 aluminum.     

Alloying element vaporization and weld pool temperature during laser welding of AlSl 202 stainless steel

Alloying element vaporization rates, plasma composition, and the changes in weld composition during laser welding of 202 stainless steel are discussed in this paper. Iron, manganese, and chromium were the most dominant species in the plasma. During laser welding it is always a difficult task to measure the temperature of the weld pool since this region is surrounded by hot plasma. In this paper a novel technique for the determination of weld pool temperature is presented. It is demonstrated that the relative rates of vaporization of any two elements from the molten pool can serve as an indicator of weld pool temperature, irrespective of the element pair selected. The composition of the solidified region calculated from the measured values of vaporization rate, plasma composition, and the volume of the solidified region was in good agreement with the weld composition determined by electron probe microanalyis technique.     

Comparative Study on Processing Property Between CWW CO2 Gas Shielded Welding and SAW

 Cable welding wire (CWW) CO2 gas shielded welding is an innovative process arc welding with high efficiency, high quality and low consumption, in which cable wire is used as consumable electrode. CWW CO2 gas shielded welding and submerged arc welding (SAW) are used for contrast studies on processing property of high strength steel A36 used in ship structure. The results show that the shapes of weld seam, using CWW CO2 gas shielded welding and SAW, are good and no weld defect such as air hole, flaw, slag inclusion, incomplete fusion, lack of penetration and so on are found in the weld seam. Because the rotating of arc during CWW CO2 gas shielded welding process has a strong stirring effect on molten pool, the grain in the heat affected zone (HAZ) of the joints, using CWW CO2 gas shielded welding, is small. Tensile failure positions of joints by CWW CO2 gas shielded welding and SAW are all in the base metal, but tensile strength of CWW CO2 gas shielded welding joint is higher than that of SAW joint by an average of 2. 3%. The average impact energy of HAZ, using CWW CO2 gas shielded welding and SAW, is almost equal, but the average impact energy of the weld seam using CWW CO2 gas shielded welding is increased by 6%, and the average impact energy of the fusion line is increased by 7%. The 180° bending tests for the joints of CWW CO2 gas shielded welding and SAW are all qualified, and the joints hardness is all less than HV 355, but hardness of CWW CO2 gas shielded welding wire welding joint near the fusion line is obviously lower. It can be concluded that the properties of CWW CO2 gas shielded welding are better than those of the SAW joint, and CWW CO2 gas shielded welding is suitable for welding high strength steel A36 used in ship structure.     

Q345 H型钢焊接开裂原因分析

材质为Q345的热轧型钢在焊接过程中出现了开裂,导致了构件的报废。从焊缝组织成分、焊缝内应力、定位焊缝焊脚尺寸3个方面进行了分析,确认施焊过程中定位焊缝位置选择不正确、焊脚尺寸过小造成焊缝内部严重质量缺陷,从而导致焊缝开裂。     

Computational modeling of stationary gastungsten-arc weld pools and comparison to stainless steel 304 experimental results

A systematic study was carried out to verify the predictions of a transient multidimensional computational model by comparing the numerical results with the results of an experimental study. The welding parameters were chosen such that the predictions of the model could be correlated with the results of an earlier experimental investigation of the weld pool surface temperatures during spot gas-tungsten-arc (GTA) welding of Type 304 stainless steel (SS). This study represents the first time that such a comprehensive attempt has been made to experimentally verify the predictions of a numerical study of weld pool fluid flow and heat flow. The computational model considers buoyancy and electromagnetic and surface tension forces in the solution of convective heat transfer in the weld pool. In addition, the model treats the weld pool surface as a truly deformable surface. Theoretical predictions of the weld pool surface temperature distributions, the cross-sectional weld pool size and shape, and the weld pool surface topology were compared with corresponding experimental measurements. Comparison of the theoretically predicted and the experimentally obtained surface temperature profiles indicated agreement within ±8 pct for the best theoretical models. The predicted surface profiles were found to agree within ±20 pct on dome height and ±8 pct on weld pool diameter for the best theoretical models. The predicted weld cross-sectional profiles were overlaid on macrographs of the actual weld cross sections, and they were found to agree very well for the best theoretical models.     

The study of surface-active element oxygen on flow patterns and penetration in A-TIG welding

A three-dimensional mathematical model was developed to simulate the flow patterns and temperature distributions in a moving A-TIG weld pool of 304 stainless steels with different oxygen content using PHOENICS software. It is shown that the surface-active element, oxygen, is important, because it affects the weld shape by changing the flow patterns in the weld pool. The weld bead penetration and the depth/width ratio increase first sharply and then remain nearly a constant with increasing oxygen content. Depending upon the oxygen contents, three, one, or two vortexes that have different positions, strength, and directions may be found in the weld pool. Oxygen can cause significant changes in the weld shape by varying the sign of the surface tension coefficient. The situation with the maximum surface tension moves from the edge to the center with increasing oxygen content. As oxygen content exceeds a critical value, a positive surface tension coefficient dominates the flow patterns. The vortexes with opposite directions caused by positive surface tension coefficient can efficiently transfer the thermal energy from the arc, creating a deep weld pool. The critical oxygen content increases with the increase of the welding current.     

Free surface flow and heat transfer in conduction mode laser welding

Temperature profiles and fluid flow fields in a weld pool are simulated through numerical solutions of Navier-Stokes equation and the equation of conservation of energy for low power laser welding in conduction mode. Experimentally determined weld pool surface topography, peak temperature, and the secondary dendrite arm spacings are found to be in fair agreement with the corresponding theoretically calculated values.