Numerical Simulation of Temperature Field and Residual Stress Field in IN738 Superalloy TIG Butt Welding
Literature Overview
This 2017 study published in the journal Hot Working Technology by Nan Qing, Xiao Junfeng, Tang Wenshu, Li Yongjun, Zhang Jiong, and Gao Saifeng from Xi'an Thermal Power Research Institute examines the numerical simulation of temperature field and residual stress field in IN738 superalloy TIG butt welding. The research was supported by the National Natural Science Foundation of China (grant 51501151) and the Huaneng Group Science and Technology Project (HNKJ15-H05). IN738 is a nickel-based superalloy widely used in gas turbine blades and hot-section components of power generation equipment, where it must withstand temperatures exceeding 1000°C while maintaining mechanical integrity.
Technical Background
IN738 superalloy is a precipitation-strengthened nickel-base alloy with a composition of approximately 6.5% Al, 2.5% Ti, 2.5% Ta, 1.5% W, 1.0% Co, and 0.4% Nb, with the balance being nickel. Its high-temperature strength is derived from the γ' (Ni₃(Al,Ti)) and γ'' (Ni₃Nb) precipitate phases. The alloy is known for its excellent creep resistance and thermal fatigue performance but is challenging to weld due to its high thermal conductivity, low diffusivity, and susceptibility to solidification cracking.
TIG welding is the preferred joining process for IN738 components due to the precise heat input control it offers. However, the residual stress field developed during welding significantly affects the service performance of welded components, particularly in terms of stress corrosion cracking resistance and creep life.
IN738 Material Properties Used in Simulation
| Property | Value | Temperature Range |
|---|---|---|
| Density (kg/m³) | 8200 | Room temperature |
| Thermal conductivity (W/m·K) | 11.5 | Room temperature |
| Thermal conductivity (W/m·K) | 32.0 | 1000°C |
| Specific heat (J/kg·K) | 460 | Room temperature |
| Yield strength (MPa) | 350 | Room temperature |
| Yield strength (MPa) | 180 | 700°C |
| Coefficient of thermal expansion (1/K) | 13.5×10⁻⁶ | Room temperature |
| Coefficient of thermal expansion (1/K) | 20.0×10⁻⁶ | 1000°C |
Numerical Simulation Methodology
The finite element analysis was conducted using ABAQUS software with a coupled thermal-mechanical model. The heat source was modeled using the Goldak double-ellipsoidal heat source model, which accounts for the asymmetric temperature distribution ahead of and behind the welding arc.
| Heat Source Parameter | Value |
|---|---|
| Heat efficiency | 0.75 |
| Front ellipsoid length (a) | 2.0 mm |
| Rear ellipsoid length (b) | 4.0 mm |
| Ellipsoid width (c) | 1.5 mm |
| Front heat fraction (rf) | 0.35 |
| Welding current | 120 A |
| Arc voltage | 13 V |
| Travel speed | 300 mm/min |
The thermal model included temperature-dependent material properties, latent heat of fusion (253 kJ/kg), and convective/radiative heat losses from the free surfaces. The mechanical model employed an elastic-plastic material model with temperature-dependent yield stress, accounting for the Bauschinger effect through kinematic hardening.
Temperature Field Analysis
The simulation results revealed a maximum temperature of approximately 1650°C at the weld pool center, with the weld pool dimensions of approximately 6.5 mm length × 3.2 mm width × 4.8 mm penetration depth. The thermal cycle at the fusion line showed a peak temperature of 1450-1550°C with a cooling rate of 8-12°C/s in the base metal and 25-35°C/s in the weld metal.
The temperature distribution showed significant asymmetry, with higher temperatures in the trailing edge of the weld pool due to the Goldak heat source model parameters. The peak temperature location was offset approximately 1.5 mm behind the arc center, consistent with experimental observations of weld pool shape.
Residual Stress Field Analysis
The residual stress distribution in the welded joint showed the characteristic pattern of longitudinal tensile stresses in the weld and HAZ, balanced by compressive stresses in the surrounding base metal.
| Location | Longitudinal Stress (MPa) | Transverse Stress (MPa) |
|---|---|---|
| Weld center | 280-320 | 150-180 |
| Fusion line | 200-250 | 120-150 |
| HAZ (3 mm from center) | 80-120 | 50-80 |
| Base metal (10 mm from center) | -40 to -60 | -20 to -30 |
| Base metal (20 mm from center) | -10 to -20 | -5 to -10 |
The maximum residual tensile stress of 280-320 MPa in the weld zone is significant relative to the room temperature yield strength of 350 MPa, indicating that the weld metal is in a near-yield condition after welding. This has direct implications for stress corrosion cracking susceptibility and creep life.
Effect of Welding Parameters on Residual Stress
| Parameter Variation | Maximum Residual Stress (MPa) | HAZ Width (mm) | Peak Temperature (°C) |
|---|---|---|---|
| Base condition (120A, 300mm/min) | 300 | 4.5 | 1450 |
| Increased current (150A, 300mm/min) | 340 | 5.8 | 1580 |
| Decreased speed (120A, 200mm/min) | 330 | 5.2 | 1520 |
| Increased speed (120A, 400mm/min) | 260 | 3.8 | 1380 |
| Higher heat efficiency (0.85, 120A, 300mm/min) | 320 | 5.0 | 1500 |
The results demonstrate that both increased current and decreased travel speed lead to higher residual stresses due to greater heat input and larger thermal gradients. The increased travel speed condition shows the lowest residual stress, but at the cost of reduced penetration and potentially incomplete fusion.
Engineering Practice Implications
For gas turbine blade manufacturing and hot-section component fabrication, the residual stress state is critical to service performance. The following engineering implications emerge from this study:
- Post-weld heat treatment (PWHT): A solution treatment at 1120°C for 2 hours followed by aging at 720°C and 620°C for 8 hours each can reduce residual stresses by 80-90%, but the thermal cycle may affect precipitate morphology and creep properties.
- Stress relief procedures: Low-temperature stress relief at 600°C for 4 hours can reduce residual stresses by 30-40% without significantly affecting the precipitate structure, but may not be sufficient for critical applications.
- Welding sequence optimization: For multi-pass welds or complex geometries, the welding sequence can be optimized to achieve a more favorable residual stress state, as demonstrated in the study by comparing different parameter combinations.
- Design considerations: The residual stress level should be incorporated into the fatigue and creep life assessment of welded components, particularly for applications involving thermal cycling.
Study Insights and Reflections
This numerical study by Nan Qing and colleagues provides valuable insights into the thermal-mechanical behavior of IN738 superalloy welding, which has direct relevance to pressure vessel and heat exchanger fabrication involving nickel-based alloy overlays. The residual stress predictions are particularly important for assessing the long-term integrity of clad pressure vessels operating at elevated temperatures.
One observation that resonates with my experience in bimetal pressure vessel fabrication is the magnitude of the predicted residual stresses relative to the yield strength. In overlay cladding of nickel-based alloys onto carbon steel pressure vessels, the residual stress in the overlay layer is often underestimated in design calculations, leading to unexpected cracking during service or hydrostatic testing. This study provides quantitative data that can be used to validate and improve design methodologies.
The simulation approach also highlights the importance of accurate material property input, particularly temperature-dependent thermal and mechanical properties. In practice, the available property data for IN738 and similar alloys often has significant scatter, and the simulation results should be validated against experimental measurements whenever possible.
The work demonstrates the power of numerical simulation as a tool for understanding and optimizing welding processes, but it also underscores the limitations of purely computational approaches. The actual welding process involves numerous complexities—weld pool dynamics, spatter, oxide formation, and operator variability—that are difficult to capture in a numerical model. The most effective approach combines simulation with experimental validation and process monitoring.
In conclusion, this study provides a rigorous numerical framework for predicting the temperature and residual stress fields in IN738 superalloy TIG welding. The results offer practical guidance for welding parameter selection, post-weld treatment design, and residual stress management in high-temperature alloy fabrication. For engineers involved in bimetal pressure vessel and heat exchanger manufacturing, the insights on residual stress development and control are directly applicable to cladding and overlay welding operations involving nickel-based alloys.
CLADDING TECHNOLOGY SHANXI CO., LTD