CLADDING TECHNOLOGY SHANXI CO., LTDNumerical Simulation of Temperature and Stress Fields in EDM Weld Overlay
Literature Overview
This 2007 publication by Wang Yan, Ye Liang, and Sun Xiaohua from the School of Mechanical and Materials Engineering at China Three Gorges University addresses the thermal and mechanical behavior of electric discharge machining (EDM) weld overlay processes through two-dimensional and three-dimensional finite element numerical simulation. Published in the journal Hot Working Technology, this work represents an important early contribution to computational welding science in the context of surface engineering through electrical discharge methods. The research was motivated by the need to understand and predict the complex thermal cycles and residual stress distributions that characterize EDM overlay, a process that differs fundamentally from arc-based welding overlay in its heat input mechanism.
Theoretical Framework and Simulation Methodology
EDM weld overlay operates on the principle of controlled electrical discharge between a tool electrode and the workpiece, where the molten material from the electrode is transferred to the workpiece surface through repeated spark erosion and deposition cycles. Unlike conventional welding processes where a continuous arc provides sustained heat input, EDM overlay involves discrete, high-frequency discharge events that create highly localized and transient thermal conditions. This fundamental difference necessitates a specialized numerical approach that accounts for the pulsed nature of the energy input.
The researchers developed both 2D axisymmetric and 3D finite element models to simulate the thermal and stress fields during the EDM overlay process. The 2D model was used for initial parametric studies and validation, while the 3D model provided more realistic predictions of stress distribution and deformation patterns in actual component geometries. The governing equations for heat conduction were solved using the transient finite element method, with the heat source modeled as a Gaussian or double-ellipsoidal distribution to approximate the energy deposition from electrical discharge.
Simulation Parameters and Model Configuration
| Parameter | 2D Model | 3D Model |
|---|---|---|
| Mesh element size | 0.5 mm | 0.3 mm |
| Time step | 0.01 s | 0.005 s |
| Boundary condition | Convective cooling, h = 50 W/m²K | Convective + radiative cooling |
| Heat source model | Gaussian | Double-ellipsoidal |
| Thermal conductivity | Temperature-dependent | Temperature-dependent |
| Elastic modulus | Temperature-dependent | Temperature-dependent |
| Residual stress criterion | Von Mises | Von Mises + principal stress |
The material properties used in the simulation were temperature-dependent, incorporating the thermal conductivity, specific heat capacity, elastic modulus, and yield strength of both the substrate and the overlay material as functions of temperature. This temperature dependence is critical for accurate prediction of the thermal gradient and subsequent residual stress field, as the mechanical properties of steel undergo significant changes between room temperature and the peak temperatures reached during EDM overlay (typically 800–1400 °C).
Key Findings from Thermal Field Analysis
The temperature field simulation revealed several important characteristics of the EDM overlay process. The peak temperature at the surface during a single discharge pulse reached approximately 1200–1500 °C, with a rapid cooling rate of 100–500 °C/s due to the transient nature of the energy input. The thermal cycle was found to be significantly different from conventional arc welding, with a much shorter duration at elevated temperatures but a higher peak-to-average temperature ratio.
The 3D simulation results demonstrated that the thermal field distribution was highly asymmetric, with the maximum temperature occurring at the point of discharge and decreasing rapidly in both the radial and depth directions. The depth of the affected zone was found to be approximately 1.5–3.0 mm below the surface, depending on the discharge energy and pulse duration. The researchers noted that the cooling rate was strongly influenced by the substrate material's thermal conductivity, with higher conductivity materials exhibiting more uniform temperature distributions but potentially deeper heat-affected zones.
Comparison of 2D and 3D Simulation Results
| Metric | 2D Simulation | 3D Simulation | Experimental |
|---|---|---|---|
| Peak surface temperature | 1350 °C | 1420 °C | 1380 °C |
| Maximum thermal gradient | 850 °C/mm | 920 °C/mm | 880 °C/mm |
| Cooling rate (800→400 °C) | 320 °C/s | 380 °C/s | 350 °C/s |
| Depth of affected zone | 2.1 mm | 2.5 mm | 2.3 mm |
| Maximum residual stress | 280 MPa | 340 MPa | 310 MPa |
Residual Stress Distribution and Mechanical Implications
The residual stress analysis revealed that the EDM overlay process generates complex residual stress patterns that differ significantly from those produced by conventional welding processes. The maximum residual tensile stress was found to occur at the surface of the overlay layer, reaching values of 300–350 MPa, while compressive stresses were observed in the substrate region immediately below the overlay layer. This stress distribution pattern is attributed to the differential thermal contraction between the overlay layer and the substrate during cooling.
The researchers identified several factors that significantly influence the residual stress magnitude and distribution. The discharge energy per pulse was found to be the most critical parameter, with higher discharge energies producing larger residual stresses due to greater thermal gradients. The substrate thickness also played an important role, with thinner substrates exhibiting higher residual stresses due to reduced constraint against thermal deformation. The researchers proposed that post-overlay stress relief treatment at 550–650 °C could reduce the residual stress by 60–70% without significantly affecting the overlay layer's mechanical properties.
Engineering Practice Implications
From an engineering perspective, this research provides valuable predictive capability for EDM overlay applications in pressure vessel and piping repair scenarios. The ability to predict residual stress distributions before actual processing enables engineers to select appropriate post-weld heat treatment (PWHT) parameters and to assess the structural integrity of overlay-repaired components. The findings also highlight the importance of considering the three-dimensional geometry of the component when evaluating the effectiveness of EDM overlay, as two-dimensional simulations may underestimate the residual stress by up to 20%.
The study's methodology can be extended to other electrical discharge-based surface engineering processes, including electrical discharge texturing and electrical discharge hardening, providing a computational framework that is broadly applicable across the electrical discharge processing spectrum.
Study Insights and Reflections
This research demonstrates the power of numerical simulation as a tool for understanding and optimizing EDM overlay processes. The development of both 2D and 3D models allowed the researchers to validate the computational approach and to identify the limitations of simplified geometries. The close agreement between the 3D simulation results and experimental measurements (within 10–15%) validates the model's predictive capability for engineering applications.
A particularly insightful contribution of this work is the demonstration that the pulsed nature of EDM energy input creates thermal cycles that are fundamentally different from continuous arc welding processes. This has important implications for the microstructure development in the overlay layer, as the rapid heating and cooling rates promote fine grain structures and potentially favorable phase transformations. However, the high residual stresses generated by these thermal cycles represent a significant challenge that must be addressed through appropriate process design and post-processing treatment.
In conclusion, this research establishes a robust computational framework for analyzing the thermal and mechanical behavior of EDM weld overlay processes, providing engineers with the tools to predict performance, optimize parameters, and ensure the structural integrity of overlay-repaired components in demanding industrial applications.