Numerical Simulation of Temperature Field and Residual Stress Field in Self-Propagating Ceramic Electrode Weld Overlay

Literature Overview

This study presents a comprehensive numerical simulation approach for analyzing the thermal and mechanical behavior during self-propagating ceramic electrode weld overlay. Self-propagating ceramic welding is a specialized technique that combines the exothermic reaction of ceramic-forming elements with arc heat input to deposit ceramic-rich overlay layers with exceptional hardness and wear resistance. The challenge in this process lies in the extreme thermal gradients and complex phase transformations that occur during and after welding, making numerical simulation an essential tool for process optimization.

The research employs finite element analysis (FEA) to model the coupled thermal-mechanical behavior during multi-pass overlay welding, providing insights into residual stress distribution, distortion patterns, and potential cracking susceptibility. This predictive capability is invaluable for optimizing process parameters and preventing defects in production environments.

Thermal Field Simulation Methodology

The thermal analysis employs a moving heat source model that accounts for the arc energy input, heat conduction through the workpiece, and heat loss through convection and radiation at exposed surfaces. The ceramic electrode introduces an additional heat source through exothermic reactions between oxygen and reactive elements such as titanium, zirconium, and aluminum in the ceramic-forming compounds.

The governing heat equation incorporates:

• Joule heating from the electric arc (typically 15-25 kW)
• Exothermic heat from ceramic formation reactions (approximately 500-800 J/g of ceramic formed)
• Thermal conduction with temperature-dependent properties
• Phase transformation heat effects during solidification and cooling
• Heat loss through convection (h = 10-50 W/m²·K) and radiation (ε = 0.8-0.95)

The simulation reveals that peak temperatures during ceramic overlay welding reach 1800-2200 degrees Celsius, significantly higher than conventional weld overlay processes. The rapid heating and cooling rates (heating rates of 100-500 degrees C per second, cooling rates of 50-300 degrees C per second) create extreme thermal gradients that drive significant residual stress development.

Residual Stress Analysis and Cracking Susceptibility

The residual stress field develops through a complex sequence of thermal expansion, plastic deformation, and phase transformation. During welding, the hot zone expands plastically while the surrounding cooler material constrains this expansion. Upon cooling, the contraction of the weld zone is restrained, generating compressive residual stresses in the overlay and tensile stresses in the heat-affected zone.

The simulation results demonstrate that peak residual stresses typically reach 200-450 MPa in the overlay region, with tensile stresses concentrated at the overlay-base metal interface. These stresses are critically important because ceramic-rich overlays are inherently brittle and susceptible to cracking under tensile stress states.

Key findings from the residual stress analysis include:

• Compressive stresses in the surface layer (50-150 MPa) that are beneficial for fatigue resistance
• Tensile stresses at the interface (200-400 MPa) that represent the primary cracking risk
• Stress relaxation through plastic deformation in the heat-affected zone
• Progressive stress accumulation with multiple overlay passes
• Significant stress concentration at weld toe regions and geometric discontinuities

Process Optimization Based on Simulation Results

The numerical simulation provides a basis for optimizing process parameters to minimize residual stresses and prevent cracking. The following strategies emerge from the analysis:

1. Travel speed optimization: Higher travel speeds reduce heat input per unit length, decreasing the thermal gradient and associated residual stresses. However, excessive speed leads to incomplete fusion and reduced overlay thickness. The optimal range is typically 150-250 mm/min for most applications.
2. Multi-pass strategy: Depositing multiple thin passes (2-3 mm each) rather than a single thick pass significantly reduces peak stresses. The simulation shows that stress peaks decrease by 30-50% when using three passes instead of one thick deposit.
3. Preheating: Moderate preheating (150-250 degrees C) reduces thermal gradients and peak stresses by approximately 20-30%. However, excessive preheating can promote grain growth and reduce the beneficial compressive stress at the surface.
4. Interpass temperature control: Maintaining interpass temperatures between 150-250 degrees C allows partial stress relaxation while preventing excessive grain coarsening.
5. Backing bar design: Using copper backing bars with controlled thermal conductivity can reduce the heat input into the base metal, lowering dilution and reducing thermal stresses.

Microstructure Prediction and Validation

The simulation incorporates phase transformation kinetics using the Johnson-Mehl-Avrami-Kolmogorov (JMAK) model to predict the evolution of microstructure during cooling. The predicted microstructure includes:

• Coarse ceramic phases (TiC, ZrC, Al2O3) in the overlay
• Martensitic matrix with varying carbon content
• Decarburization zone at the interface
• Grain growth in the heat-affected zone

Experimental validation through metallographic examination and X-ray diffraction confirms the simulation predictions with good agreement. Measured residual stresses using the hole-drilling method show deviations of less than 15% from simulated values, validating the model's predictive capability.

Engineering Practice Integration

In practical applications, the simulation results guide the following engineering decisions:

• Selection of appropriate backing materials to control heat flow
• Design of stress-relief grooves or slots to accommodate thermal expansion
• Determination of optimal welding sequence for multi-pass overlays on large components
• Prediction of distortion for component fabrication planning
• Assessment of post-weld stress relief requirements

The simulation also identifies critical regions where cracking is most likely to initiate, enabling targeted process modifications. For example, when overlaying cylindrical surfaces, the simulation predicts higher stresses at the weld termination points, suggesting the use of run-off tabs and proper start/stop techniques.

Study Insights and Reflections

This research demonstrates the tremendous value of numerical simulation in understanding and optimizing complex welding processes like ceramic electrode overlay. The ability to predict thermal and stress fields before actual welding provides engineers with a powerful tool for process development and quality assurance.

The key insight is that ceramic overlay welding operates in an extreme thermal regime that cannot be adequately understood through experimental methods alone. The simulation reveals phenomena that are difficult to measure experimentally, such as transient stress states during welding and the evolution of stress through multiple passes.

However, the study also acknowledges limitations. The accuracy of predictions depends heavily on accurate material property data at extreme temperatures, which remains challenging to obtain for ceramic-rich materials. Future work should focus on developing more accurate constitutive models for ceramic-matrix composites and incorporating stochastic elements to account for material variability.

From a practical standpoint, engineers should view simulation results as guiding information rather than absolute predictions. The combination of simulation-based process planning with experimental validation provides the most reliable approach to optimizing ceramic overlay welding processes for specific applications.