Three-Dimensional Finite Element Simulation of TIG Welding Temperature Field in Low-Alloy Steel Thin Plates
Literature Overview and Research Background
The study authored by Guo Yanbing, Tong Yangang, and He Xiaona from the School of Materials Science and Engineering, Chongqing University, published in Hot Working Technology in 2010, addresses a critical gap in the computational understanding of gas tungsten arc welding (GTAW/TIG) processes applied to low-alloy steel thin plates. Thin-plate welding remains one of the most challenging aspects of GTAW fabrication because the combination of low thermal mass and high heat flux concentration produces steep thermal gradients, pronounced residual stresses, and significant distortion risks. The authors recognized that conventional one-dimensional or two-dimensional thermal models fail to capture the true three-dimensional heat flow behavior in thin sections, particularly near the weld root and cap transitions.
Core Technical Approach
The research employs a three-dimensional finite element method (FEM) framework to simulate the transient temperature field during TIG welding of low-alloy steel thin plates. The numerical model incorporates a moving heat source, typically represented by a Gaussian or double-ellipsoidal distribution, to account for the directional nature of the arc energy input. The following key modeling parameters are essential to achieving accurate predictions:
| Parameter | Typical Value / Range | Notes |
|---|---|---|
| Heat source type | Double-ellipsoidal (Goldak) | Captures front/back heat deposition asymmetry |
| Welding speed | 200–500 mm/min | Depends on plate thickness and joint design |
| Arc voltage | 12–20 V | Governs arc power and penetration |
| Welding current | 80–150 A | Directly affects HAZ width and distortion |
| Thermal conductivity | 45–55 W/(m·K) | Temperature-dependent for low-alloy steels |
| Specific heat | 450–700 J/(kg·K) | Must account for phase transformation |
| Latent heat | 250–300 kJ/kg | Critical near solidus/liquidus transitions |
| Element type | 8-node hexahedral / 4-node tetrahedral | Mesh refinement near weld zone |
The simulation considers the moving coordinate system technique, where the heat source translates along the weld path while the mesh remains fixed in the workpiece coordinate system. This approach avoids the computational expense of remeshing at each time step while preserving accuracy in the high-gradient weld zone.
Key Findings and Interpretation
The three-dimensional simulation reveals several phenomena that two-dimensional models cannot capture. First, the heat dissipation in the thickness direction creates a non-uniform temperature distribution through the plate, which directly influences the solidification morphology and microstructural evolution in the fusion zone. Second, the peak temperature at the weld surface exceeds 2000 K, while the temperature at the plate mid-thickness can differ by 200–400 K depending on the welding parameters. This through-thickness temperature differential is a primary driver of angular and longitudinal distortion in thin plates.
The study also highlights the importance of boundary condition accuracy. Convective heat transfer coefficients at the plate surfaces, particularly on the back side where no shielding gas is present, significantly affect the predicted cooling rates. A coefficient of 10–25 W/(m²·K) for natural convection and radiation combined is typical, but the authors emphasize that radiation becomes dominant at temperatures above 1000 K and must be modeled with the Stefan-Boltzmann law rather than a constant coefficient.
Engineering Practice Implications
From a practical standpoint, the FEM temperature field simulation provides engineers with a predictive tool for optimizing TIG welding parameters on low-alloy steel thin plates before committing to physical trials. The following engineering insights can be extracted:
- Preheat control: For low-alloy steels susceptible to cold cracking, the simulation can identify the minimum preheat temperature required to keep the HAZ above the martensite start temperature, typically in the range of 100–200°C for steels with carbon equivalents of 0.4–0.6.
- Interpass temperature management: In multi-pass welding of thin plates, the simulated temperature field allows prediction of the interpass temperature after each pass, enabling precise control to prevent excessive grain growth or hardness in the HAZ.
- Distortion prediction: The thermal strain distribution derived from the temperature field can be coupled with a mechanical FEM analysis to predict angular and longitudinal distortion, guiding the design of clamping fixtures and backing plates.
- Cooling rate optimization: The 8/5 cooling time (time from 800°C to 500°C) is a critical parameter for microstructural control. The simulation provides spatial maps of cooling rate, allowing identification of regions where the cooling rate may be too fast for thick sections or too slow for thin sections.
Critical Reflection and Limitations
While the three-dimensional FEM approach represents a significant advancement over lower-dimensional models, several limitations must be acknowledged. The model assumes a quasi-static heat source distribution that may not fully capture the dynamic behavior of the arc during actual welding, particularly in situations involving arc wandering or current pulsation. Furthermore, the thermal properties used in the simulation are typically derived from bulk material data and may not accurately represent the properties of the weld metal, which can differ substantially from the base metal. The model also does not account for the thermal effects of shielding gas flow, which can influence the back-side cooling rate and the formation of oxidation inclusions.
Despite these limitations, the study establishes a robust computational framework that, when validated against thermocouple measurements and thermographic imaging, provides reliable predictions for welding process design. The integration of such FEM tools into the welding procedure qualification process represents a paradigm shift from purely trial-and-error methods to a physics-based, predictive approach that reduces development time and material waste.
Summary
The three-dimensional FEM simulation of the TIG welding temperature field for low-alloy steel thin plates, as presented by Guo and colleagues, demonstrates that computational modeling is an indispensable tool for understanding and controlling the thermal phenomena that govern weld quality, residual stress, and distortion in thin-section GTAW fabrication. The study provides a methodological foundation that can be extended to include coupled thermo-mechanical analysis, phase transformation modeling, and microstructural prediction, ultimately enabling more efficient and reliable welding procedure development for critical low-alloy steel components.
CLADDING TECHNOLOGY SHANXI CO., LTD