CLADDING TECHNOLOGY SHANXI CO., LTDMathematical Modeling of TIG Weld Penetration and Weld Pool
Literature Overview
This 1996 study by Cao Zhenning, Wu Chuansong, and Wu Lin, published in the "Welding Journal" (焊接学报) and supported by the National Natural Science Foundation of China, presents a mathematical model for predicting weld penetration and weld pool geometry in gas tungsten arc welding. The work was conducted at Shandong University of Technology and Harbin Institute of Technology, representing a significant contribution to the theoretical understanding of arc welding physics during a period when computational modeling was emerging as a critical engineering tool.
Theoretical Framework
The mathematical model developed in this study addresses the fundamental question of how arc parameters translate into weld geometry. The governing equations encompass:
- Heat transfer equation in the solid and liquid phases, incorporating conduction, convection (in the liquid pool), and phase change effects.
- Momentum equation for the fluid flow within the weld pool, accounting for buoyancy forces, electromagnetic forces, and surface tension gradients (Marangoni convection).
- Mass conservation with appropriate boundary conditions at the free surface of the weld pool.
- Arc heat flux distribution, typically modeled as a Gaussian or double-Gaussian distribution over the workpiece surface.
The key boundary condition at the weld pool surface involves the balance of heat flux from the arc, radiative heat loss, convective heat loss to the shielding gas, and the temperature-dependent surface tension gradient that drives Marangoni convection.
Key Model Assumptions and Parameters
| Parameter | Typical Value | Influence on Penetration |
|---|---|---|
| Arc power density | 10⁷–10⁸ W/m² | Higher density → deeper penetration |
| Surface tension coefficient (Al) | 0.10–0.12 N/m | Negative gradient → shallow wide pool |
| Thermal conductivity | 100–200 W/(m·K) | Higher → wider shallow pool |
| Density | 2000–2700 kg/m³ | Affects buoyancy-driven flow |
| Welding current | 50–300 A | Primary control on heat input |
| Travel speed | 1–10 mm/s | Inversely related to penetration |
| Arc length | 2–5 mm | Shorter → more concentrated heat |
Model Validation and Engineering Applications
The mathematical model was validated against experimental measurements of weld geometry for various materials including carbon steel, stainless steel, and aluminum alloys. The comparison between predicted and measured penetration depth, weld width, and reinforcement typically showed agreement within 10–15%, which was considered acceptable for process design purposes at that time.
The practical applications of this model extend across multiple domains:
- Process optimization: Determining the minimum current required for full penetration of a given plate thickness without excessive heat input.
- Weld pool monitoring: Predicting the weld pool shape in real-time based on measured arc parameters, enabling feedback control.
- Residual stress prediction: The temperature field solution provides input for thermoelastic-plastic analysis of residual stress distribution.
- Distortion prediction: Coupling the thermal model with structural analysis enables prediction of angular and longitudinal distortion.
Relevance to Cladding and Overlay Welding
For cladding and weld overlay applications, the mathematical model of weld pool behavior is particularly relevant because:
- Dilution control: The model predicts the volume fraction of base metal incorporated into the weld, which directly determines the dilution ratio. In overlay welding, dilution is a critical quality parameter that affects the corrosion resistance of the cladding layer.
- Layer thickness prediction: By understanding the weld pool geometry, engineers can predict the deposited layer thickness per pass and plan multi-pass overlay sequences accordingly.
- Interpass temperature management: The model provides the thermal field solution needed to determine the optimal interpass temperature that prevents cracking while maintaining adequate bond strength.
- Bond strength assessment: The temperature gradient at the interface between the overlay layer and the base material influences the formation of intermetallic compounds and the resulting bond strength.
Study Insights and Reflections
This mathematical modeling work represents a foundational contribution to the field of welding science. The transition from purely empirical process development to physics-based modeling enabled engineers to predict welding outcomes with greater confidence and to optimize processes systematically rather than through trial and error. For the cladding industry, this work laid the groundwork for the computational tools now routinely used in process qualification, particularly for complex overlay geometries where experimental testing alone would be prohibitively expensive. The enduring value of this research lies in its demonstration that a rigorous physical model, even with simplifying assumptions, provides engineering insights that cannot be obtained from experimental data alone.