Numerical Simulation of Temperature and Residual Stress Fields in Self-Propagating Ceramic Electrode Hardfacing

Literature Overview

This research paper by Chen Wei and colleagues from the Department of Mechanical Engineering, Academy of Armored Force Engineering, published in Hot Working Technology in 2010 (Vol. 39, No. 5, pp. 11-13), presents a numerical simulation study of the temperature field and residual stress field during hardfacing with self-propagating high-temperature synthesis (SHS) TiC-Ni ceramic electrodes using GTAW (gas tungsten arc welding) technology. The study was supported by the National Natural Science Foundation of China (grant numbers 50675222, 50575226, 50275149). The research establishes a mathematical model for the hardfacing process and uses finite element analysis to predict and analyze the temperature distribution and residual stress evolution during the deposition of metal-ceramic hardfacing layers.

Technical Background

Self-propagating ceramic electrodes represent an advanced class of hardfacing materials that combine metallic bonding strength with the exceptional hardness and wear resistance of ceramic phases. TiC-Ni SHS electrodes produce hardfacing deposits containing titanium carbide (TiC) particles dispersed in a nickel matrix, achieving hardness levels that far exceed conventional metallic hardfacing alloys. However, the introduction of ceramic phases into the deposit introduces significant challenges related to thermal stress, cracking, and delamination.

The fundamental challenge in hardfacing with ceramic-containing electrodes is the mismatch between the thermal expansion coefficients of the ceramic phase, the metallic matrix, and the steel substrate. TiC has a thermal expansion coefficient of approximately 7.4 × 10⁻⁶ /K, while austenitic stainless steel substrates typically have coefficients in the range of 16-18 × 10⁻⁶ /K. This mismatch generates substantial residual stresses during cooling, which can lead to cracking and delamination if not properly managed.

Mathematical Model Development

The numerical simulation requires the establishment of a comprehensive mathematical model that accounts for:

Thermal Model

The thermal model is based on the heat conduction equation with moving heat source:

• Heat source: Gaussian or double-ellipsoidal distribution representing the GTAW arc
• Boundary conditions: Convective and radiative heat loss from the workpiece surface
• Initial conditions: Uniform ambient temperature
• Material properties: Temperature-dependent thermal conductivity, specific heat, and density
• Phase transformation: Latent heat effects during solidification and solid-state phase transformations

Stress Model

The residual stress model incorporates:

• Elastic-plastic constitutive law with temperature-dependent yield stress
• Thermal strain due to temperature gradients
• Plastic strain accumulation during welding
• Phase transformation strain (martensitic transformation, if applicable)
• Constraint conditions imposed by the fixture and the surrounding cold material

Key Assumptions and Simplifications

Simulation Results and Analysis

Temperature Field Distribution

The temperature field simulation reveals several important characteristics:

• The peak temperature at the weld pool center reaches values significantly above the melting point of both the electrode and substrate materials.
• The thermal gradient is steep in the direction perpendicular to the weld axis, with rapid temperature decrease from the weld center to the surrounding base material.
• The thermal cycle experienced by any given point in the workpiece is characterized by a rapid heating phase during the arc passage and a slower cooling phase as the arc moves away.
• Multi-pass hardfacing results in complex thermal histories, with each subsequent pass reheating the previously deposited layer and modifying its microstructure.

Residual Stress Field Distribution

The residual stress simulation provides critical insights into the stress state of the hardfacing deposit:

• Compressive residual stresses are generated in the weld pool region during cooling, as the contracting weld metal is constrained by the surrounding cooler material.
• Tensile residual stresses develop in the heat-affected zone and the substrate as the contracting weld metal pulls on the surrounding material.
• The stress distribution exhibits a characteristic pattern: compressive in the weld center, transitioning to tensile in the HAZ, and returning to near-zero in the far-field substrate.
• The magnitude of residual stresses is influenced by the thermal expansion coefficient mismatch, the cooling rate, and the constraint imposed by the surrounding material.

Stress Continuity and Gradation

A key finding of the simulation is the verification of stress continuity and gradation within the hardfacing layer. The study specifically examines whether the residual stress field exhibits smooth transitions or abrupt discontinuities. The results confirm that:

• The residual stress field is continuous across the deposit-substrate interface, with no abrupt jumps in stress magnitude.
• The stress gradient is moderate, indicating that the stress state does not create localized stress concentrations that would promote crack initiation.
• The stress distribution is consistent with the expected mechanical behavior of the hardfacing deposit under thermal cycling.

Engineering Implications

Process Optimization

The numerical simulation results can be used to optimize the hardfacing process parameters:

• Welding current and voltage can be adjusted to control peak temperature and thermal gradient.
• Travel speed can be modified to control cooling rate and residual stress magnitude.
• Interpass temperature can be set to minimize the thermal cycling amplitude.
• Welding sequence can be planned to distribute thermal input evenly and minimize stress accumulation.

Crack Prediction and Prevention

The residual stress analysis provides a basis for crack prediction:

• Regions of high tensile residual stress are identified as potential crack initiation sites.
• The magnitude of residual stresses can be compared with the fracture toughness of the hardfacing deposit to assess cracking risk.
• Process modifications can be designed to reduce residual stresses below the cracking threshold.
• Post-weld stress relief procedures can be specified based on the predicted stress distribution.

Integration with Experimental Validation

While numerical simulation provides valuable predictive capability, it must be validated against experimental data. The study's approach of using simulation to guide experimental design and then using experimental results to refine the simulation model represents a sound engineering methodology. Key experimental validation parameters include:

• Temperature measurements using embedded thermocouples or infrared thermography
• Residual stress measurements using X-ray diffraction or hole-drilling methods
• Microstructural examination using metallography and scanning electron microscopy
• Mechanical property testing including hardness profiling and tensile testing

Study Insights and Implications

This study demonstrates the value of numerical simulation as a tool for understanding and controlling the hardfacing process. The ability to predict temperature and residual stress distributions before conducting physical experiments saves significant time and resources, while providing a deeper understanding of the underlying physical mechanisms. For engineers working on advanced hardfacing applications involving ceramic-containing deposits, numerical simulation is an essential tool for process development and optimization.

The verification of stress continuity and gradation is particularly important for ensuring the long-term reliability of ceramic hardfacing deposits. Abrupt stress discontinuities would create localized stress concentrations that could initiate cracks and lead to premature failure. The simulation results confirm that, with appropriate process design, the residual stress field can be managed to ensure smooth transitions and acceptable stress magnitudes throughout the deposit. This provides confidence in the structural integrity of ceramic hardfacing deposits and supports their application in demanding industrial environments.