Coupled Field Modeling of Gas Tungsten Arc Welding

Welding is used extensively in aerospace, automotive, chemical, manufacturing, electronic and power-generation industries. Thermally-induced residual stresses due to welding can significantly impair the performance and reliability of welded structures. Numerical simulation of weld pool dynamics is important as experimental measurements of velocities and temperature profiles are difficult due to the small size of the weld pool and the presence of the arc. From a structural integrity perspective of welded structures, it is necessary to have an accurate spatial and temporal thermal distribution in the welded structure before stress analysis is performed. Existing research on weld pool dynamics simulation has ignored the effect of fluid flow in the weld pool on the temperature field of the welded joint. Previous research has established that the weld pool depth/width (D/W) ratio and Heat Affected Zone (HAZ) are significantly altered by the weld pool dynamics. Hence, for a more accurate estimation of the thermally-induced stresses it is desired to incorporate the weld pool dynamics into the analysis. Moreover, the effects of microstructure evolution in the HAZ on the mechanical behavior of the structure need to be included in the analysis for better mechanical response prediction. In this study, a three-dimensional model for the thermo-mechanical analysis of Gas Tungsten Arc (GTA) welding of thin stainless steel butt-joint plates has been developed. The model incorporates the effects of thermal energy redistribution through weld pool dynamics into the structural behavior calculations. Through material modeling the effects of microstructure change/phase transformation are indirectly included in the model. The developed weld pool dynamics model includes the effects of current, arc length, and electrode angle on the heat flux and current density distributions. All the major weld pool driving forces are included, namely surface tension gradient, plasma drag force, electromagnetic force, and buoyancy. The weld D/W predictions are validated with experimental results. They agree well. The effects of welding parameters (like welding speed, current, arc length, etc.) on the weld D/W ratio are documented. The workpiece deformation and stress distributions are also highlighted. The transverse and longitudinal residual stress distribution plots across the weld bead and their variations with welding speed and current are also provided. The mathematical framework developed here serves as a robust tool for better prediction of weld D/W ratio and thermally-induced stress evolution and distribution in a welded structure by coupling the different fields in a welding process.