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

1. Definition and Fundamental Principles

Self-Propagating High-Temperature Synthesis (SHS), also referred to as self-propagating ceramic electrode weld overlay, is a thermal-mechanical cladding technology that exploits exothermic chemical reactions between ceramic-forming elements (typically boron, silicon, chromium, and oxygen) embedded in a specialized composite electrode and the substrate material. Unlike conventional arc weld overlay processes that rely entirely on external energy input to melt filler material, SHS leverages the substantial heat generated by internal exothermic reactions to sustain a self-propagating molten pool that progressively deposits a hard, wear-resistant ceramic-metallic composite layer onto the substrate surface.

The core principle involves a composite electrode consisting of a metallic core (often low-carbon steel or alloy steel) encased in a ceramic-forming sheath. When an arc is initiated at the leading edge of the electrode, the intense localized heating triggers a vigorous exothermic reaction between the boron-silicon-chromium mixture and the molten base metal. The heat release from these reactions (typically reaching temperatures of 2,500–3,500°C in the reaction zone) exceeds the energy input from the arc, enabling the molten pool to propagate along the electrode length without continuous external energy dependence. The resulting deposit forms a gradient transition from a metallic intermetallic zone at the interface to a fully ceramic (boride, carbide, silicide) structure at the surface, providing exceptional hardness (HV 1,200–2,000+) and wear resistance.

Numerical simulation of the temperature field and residual stress field during this process is critical because the extremely rapid heating and cooling rates (on the order of 10²–10³ °C/s), coupled with the heterogeneous microstructure of the deposit, create complex thermal gradients and mechanical stresses that are difficult to characterize experimentally alone. Finite Element Analysis (FEA) provides the necessary insight into thermal cycling behavior, phase transformation kinetics, and residual stress distribution to optimize process parameters and predict service performance.

2. Category and Business Positioning

Within the cladding technology landscape, self-propagating ceramic electrode weld overlay occupies a unique position as a hybrid technology that bridges traditional arc welding with in-situ metallurgical synthesis. It belongs to the weld overlay category but is distinguished from conventional TIG/MIG weld overlay by its fundamentally different energy source and microstructural formation mechanism. In the company's three primary technology routes—TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding—SHS ceramic electrode overlay serves as a specialized, high-value-added extension of the weld overlay route, targeting applications requiring extreme wear and corrosion resistance where conventional overlay deposits are insufficient.

The numerical simulation capability associated with this technology represents an intellectual property asset that differentiates the company in the following ways:

Process Optimization: Enables virtual process development without expensive trial-and-error physical testing, accelerating time-to-market for new cladding solutions.
Quality Assurance: Provides predictive models for residual stress assessment, reducing the risk of cracking, delamination, and premature failure in service.
WPS Qualification Support: Generates the thermal cycle data and stress analysis required to develop and qualify Welding Procedure Specifications (WPS) for complex substrate geometries and thick cladding builds.
Customer Engineering Support: Offers detailed simulation reports to customers demonstrating predicted performance, service life, and failure modes, enhancing confidence in the delivered product.

3. Technical Purpose and Value

3.1 Temperature Field Simulation Objectives

The primary objective of temperature field numerical simulation is to accurately predict the thermal history experienced at every point within the deposit, interface, and substrate during and after the SHS process. Key parameters of interest include:

Peak Temperature Distribution: Determining maximum temperatures at the reaction front, deposit surface, and substrate interface to validate melting and reaction conditions.
Heating and Cooling Rates: Quantifying the rate of thermal cycling, which directly influences grain structure, phase formation, and residual stress generation.
Thermal Gradient Magnitude: Assessing the steepness of temperature gradients across the deposit-substrate interface, which is a primary driver of thermal mismatch stress.
Heat Input Distribution: Evaluating how the combined arc energy and exothermic reaction heat distribute spatially and temporally.
Molten Pool Geometry: Predicting the shape, depth, and width of the molten pool, which determines dilution rates and bond strength.

3.2 Residual Stress Field Simulation Objectives

Residual stress analysis is equally critical because the extreme thermal gradients and heterogeneous material properties in SHS deposits create significant locked-in stresses. The simulation aims to:

Quantify Stress Magnitudes: Determine the magnitude of residual stresses (tensile and compressive) in the deposit, transition zone, and substrate.
Identify Critical Stress Concentrations: Locate positions where stress concentrations may initiate cracking or delamination.
Predict Stress Evolution: Model how residual stresses develop through the sequential stages of heating, reaction, solidification, and cooling.
Evaluate Post-Weld Heat Treatment Effects: Simulate the stress relief potential of various PWHT regimes to develop effective mitigation strategies.
Assess Multi-Pass Stress Accumulation: For multi-pass builds, model how stresses from successive passes interact and accumulate.

4. Key Process Parameters and Implementation Points

4.1 Governing Process Parameters for Simulation

Parameter Category	Specific Parameter	Typical Range	Simulation Relevance
Arc Energy Input	Welding Current	80–200 A (DC)	Initial heat source intensity; affects ignition and early-stage thermal profile
Arc Energy Input	Arc Voltage	18–32 V	Determines arc power density and heat input rate
Arc Energy Input	Travel Speed	100–400 mm/min	Controls heat input per unit length and molten pool geometry
Exothermic Reaction	Heat of Reaction	1,500–3,000 J/g	Primary energy source sustaining the self-propagating molten pool
Exothermic Reaction	Reaction Zone Temperature	2,500–3,500°C	Determines melt extent, reaction kinetics, and dilution behavior
Exothermic Reaction	Composition (B, Si, Cr, O)	Varies by formulation	Influences thermophysical properties and phase formation
Thermal Properties	Thermal Conductivity (k)	15–45 W/m·K (deposit); 45–55 W/m·K (steel)	Governs heat dissipation rate and thermal gradient steepness
Thermal Properties	Specific Heat (cp)	500–900 J/kg·K	Affects thermal inertia and cooling rate
Thermal Properties	Density (ρ)	2,500–7,800 kg/m³	Varies significantly between ceramic and metallic phases
Mechanical Properties	Coefficient of Thermal Expansion (CTE)	4–12 × 10⁻⁶ /°C	Primary driver of thermal mismatch stress at interface
Mechanical Properties	Elastic Modulus (E)	200–400 GPa (deposit); 200 GPa (steel)	Determines stress response to thermal strain
Process Configuration	Electrode Diameter	10–25 mm	Affects heat input per pass and deposit thickness
Process Configuration	Substrate Thickness	5–50+ mm	Influences heat sink effect and thermal boundary conditions

4.2 Simulation Methodology and Implementation Steps

The numerical simulation of SHS ceramic electrode weld overlay follows a structured methodology:

Geometry Modeling: Create a 3D finite element model of the substrate, electrode, and expected deposit geometry. The model must account for the actual component geometry, including any curvature, thickness variations, or complex features.
Material Property Definition: Assign temperature-dependent thermophysical properties (thermal conductivity, specific heat, density, CTE, elastic modulus, yield strength) to each material region: substrate, metallic core, ceramic-forming sheath, and resulting composite deposit. These properties must reflect the phase evolution during cooling.
Heat Source Modeling: Implement a dual heat source model combining:
- An arc heat source (typically a Gaussian or double-elliptical distribution) representing the external energy input.
- A moving volumetric heat source representing the exothermic reaction zone, with heat generation rate determined by reaction kinetics and composition.
Thermal Analysis: Perform a transient thermal simulation following the actual welding sequence. The moving heat source travels along the electrode path, and the temperature field is computed at each time step. Boundary conditions include convective and radiative heat loss from exposed surfaces.
Mechanical/Thermoelastic-Plastic Analysis: Using the temperature field from the thermal analysis, perform a coupled thermoelastic-plastic analysis. The material is constrained from contraction during cooling (since it was deposited in a molten state), generating compressive stresses that relax plastically. Upon subsequent cooling, the deposit contracts against the substrate, generating tensile residual stresses in the deposit and compressive stresses in the substrate.
Post-Processing and Validation: Extract temperature histories, thermal gradients, residual stress distributions, and equivalent plastic strain from the simulation results. Validate against experimental measurements (thermocouple readings, X-ray diffraction stress measurements, or metallographic examination) to refine model accuracy.

4.3 Critical Simulation Considerations

Phase Transformation Modeling: The SHS deposit undergoes complex phase evolution during cooling, transitioning from a fully molten state through various intermetallic phases to the final ceramic structure. These phase changes involve volume changes and property transitions that significantly affect residual stress development. The simulation should incorporate phase transformation kinetics using models such as JMA (Johnson-Mehl-Avrami) or Scheil-Gulliver for solidification.
Heterogeneous Material Properties: Unlike homogeneous weld overlays, SHS deposits have a gradient microstructure with properties varying continuously from the substrate interface to the deposit surface. The simulation must use a spatially varying property field rather than a uniform material assignment.
Thermal Contact Conductance: The interface between successive passes (in multi-pass builds) or between the electrode and substrate involves imperfect thermal contact. Modeling this contact resistance is essential for accurate temperature prediction at critical interface locations.
Preheating Effects: Substrate preheating significantly alters the thermal cycle and residual stress state. The simulation should incorporate realistic preheating conditions, including temperature uniformity across the heated zone.

5. Applicable Standards and Acceptance Criteria

5.1 Weld Overlay and Cladding Standards

ASTM A240 / ASTM A213: Specifications for the substrate materials (stainless steel plates, tubes) to which SHS overlay is applied.
ASTM B136: Standard Specification for Weld Overlay (Cladding) of Steel and Cast Iron with Nickel-Cobalt-Copper Alloys—provides reference methodology for cladding qualification testing.
ASME Section IX, Part Q: Qualification of Welding, Brazing, and Fusing Procedures—governs WPS/PQR qualification for the overlay process, including thermal cycle documentation.
ASME BPVC Section II, Part D: Impact Test Requirements—relevant for thick-section substrates where residual stress may affect low-temperature toughness.
GB/T 11351: Technical conditions for clad steel plates and strips—Chinese national standard for clad plate acceptance.
GB/T 25670: Technical conditions for clad steel pipes—Chinese national standard for clad pipe acceptance.
NB/T 20339: Technical conditions for pressure vessel cladding—relevant for nuclear and pressure vessel applications.
ISO 14224: Steel plates, strips and sheets—clad steel—dimensions, mass tolerances and general technical delivery conditions.

5.2 Residual Stress and Post-Weld Heat Treatment Standards

ASME BPVC Section VIII, Div. 1, UG-113: Post-Weld Heat Treatment requirements for pressure vessels—defines PWHT parameters that must be validated through simulation for complex geometries.
NACE MR0175 / ISO 15156: Materials for use in H₂S-containing environments—residual stress levels must be controlled to prevent sulfide stress cracking in overlay deposits.
ASTM E923: Standard Practice for Determination of Residual Stresses by the Hole-Drilling Strain-Gauge Method—reference method for experimental validation of simulated residual stress values.
ASTM E765: Standard Test Method for Measurement of Residual Stress Using X-Ray Diffraction—primary validation technique for surface residual stress measurement.

5.3 Acceptance Criteria for Simulated Results

Parameter	Acceptance Threshold	Verification Method
Peak residual tensile stress in deposit	≤ 300 MPa (or ≤ 0.4 × UTS of deposit)	FEA output cross-validated with XRD measurements
Peak residual stress at interface	≤ 250 MPa tensile	FEA output cross-validated with hole-drilling method
Peak temperature at substrate surface	≤ 1,200°C (to prevent substrate microstructural degradation)	Thermocouple validation of thermal simulation
Maximum cooling rate (800°C to 600°C)	≤ 50°C/s (for ductile substrate); ≤ 200°C/s (acceptable for overlay)	Thermal simulation output vs. thermocouple data
Equivalent plastic strain	≤ 0.02 (to prevent cracking initiation)	FEA output vs. metallographic examination
PWHT effectiveness (stress reduction)	≥ 70% reduction in peak residual stress	Simulated PWHT cycle vs. post-PWHT XRD measurements

6. Common Risks and Controls

6.1 Process-Related Risks

Cracking in the Deposit: The high cooling rates and high CTE of ceramic phases create severe tensile stresses upon cooling. Control: Simulation identifies critical stress concentrations; process parameters (travel speed, electrode composition) are optimized to reduce peak stress below cracking threshold. Preheating and interpass temperature control are simulated and prescribed.
Delamination at Interface: Thermal mismatch between the hard ceramic deposit and the ductile substrate can cause interface cracking. Control: Numerical analysis quantifies interface stress; a transition layer strategy (graded composition or intermediate alloy pass) is designed and simulated to reduce CTE mismatch.

Excessive Dilution: If the molten pool penetrates too deeply into the substrate, the resulting dilution can compromise the hardness and wear resistance of the deposit. Control: Temperature field simulation predicts dilution depth; process parameters are adjusted to maintain dilution within acceptable limits (typically 10–30% substrate dilution).

Incomplete Fusion: Insufficient heat input at the trailing edge of the reaction zone can result in incomplete fusion between passes. Control: Thermal simulation verifies that the trailing edge temperature remains above the solidus temperature; overlap between passes is optimized through simulation.

6.2 Simulation-Related Risks

Model Over-Simplification: Excessive simplification of material properties (e.g., assuming linear elasticity, ignoring phase transformation) can lead to inaccurate stress predictions. Control: Use of coupled thermoelastic-plastic models with temperature-dependent properties and phase transformation kinetics; validation against experimental data.

Heat Source Model Inaccuracy: The exothermic reaction heat source is inherently difficult to model due to its dependence on composition, temperature, and local geometry. Control: Calibration of the reaction heat source model against measured temperature profiles; use of simplified but validated reaction kinetics models.

Boundary Condition Assumptions: Incorrect assumptions about convective/radiative heat loss or substrate boundary conditions can skew temperature predictions. Control: Sensitivity analysis of boundary conditions; use of measured surface temperatures to calibrate heat transfer coefficients.

Multi-Scale Property Variation: The SHS deposit has properties varying over length scales from micrometers (grain level) to millimeters (pass level). Control: Use of homogenized property fields at the macro scale; sensitivity analysis to assess the impact of property field variations on residual stress predictions.

7. Application Scenarios Across Technology Routes

7.1 Integration with TIG/MIG Weld Overlay Route

The SHS ceramic electrode overlay technology extends the conventional TIG/MIG weld overlay capability into the realm of extreme wear resistance. The numerical simulation capability developed for SHS processes is directly transferable to conventional weld overlay in the following ways:

Transition Layer Design: Simulation of the temperature and stress fields enables the design of optimal transition layer sequences (e.g., 309L → 310 → hardfacing alloy) that minimize residual stress at each interface. This is particularly important for multi-pass builds where each successive pass imposes new thermal cycles on previously deposited material.

Multi-Pass Stress Management: For thick overlay builds (common in mining equipment, paper machine rolls, and cement kiln components), simulation predicts how residual stresses accumulate across multiple passes. This enables the design of pass sequences and interpass temperatures that minimize cumulative stress.

WPS Qualification: The thermal cycle data generated from simulations supports ASME Section IX Part Q qualification by providing the documented thermal history required for procedure qualification, reducing the number of physical qualification coupons required.

Post-Weld Heat Treatment Optimization: Simulation of PWHT cycles allows the determination of optimal temperature, soak time, and cooling rate to achieve maximum stress relief without compromising the microstructural integrity of the overlay. This is particularly valuable for complex geometries where uniform PWHT is challenging.

7.2 Integration with Hydraulic Explosive Bonding Route

While hydraulic explosive bonding (hydraulic explosion cladding) is a fundamentally different technology that relies on controlled impact velocities to achieve metallurgical bonding, the numerical simulation expertise developed for SHS weld overlay contributes to hydraulic bonding in several ways:

Post-Bonding Heat Treatment Simulation: Hydraulic explosive bonds often require post-bonding heat treatment to improve interface quality and relieve residual stresses. The thermal stress simulation methodology developed for SHS overlay is directly applicable to modeling these PWHT cycles, predicting stress relief effectiveness, and identifying optimal PWHT parameters.

Interface Stress Assessment: The residual stress analysis techniques developed for SHS overlay interfaces (with their steep property gradients) are applicable to the assessment of residual stresses at the hydraulic explosive bond interface, particularly in thick-section applications where stress relaxation is incomplete.

Multi-Process Cladding Design: In hybrid cladding configurations where hydraulic explosive bonding provides the primary cladding layer and weld overlay (including SHS) provides a wear-resistant surface layer, simulation is essential for predicting the interaction between the residual stress fields from each process and ensuring overall structural integrity.

7.3 Integration with Explosion Welding Route

Explosion welding (explosive cladding) shares many characteristics with hydraulic explosive bonding and similarly benefits from the simulation expertise:

Explosion-Welded Substrate Preparation: Before applying SHS or conventional weld overlay to explosion-welded cladding, the residual stress state from the explosion welding process must be characterized. Simulation provides the analytical framework for understanding and predicting these stresses.

Overlay on Explosion-Welded Clad Plates: When SHS ceramic overlay is applied to the surface of explosion-welded clad plates (a common configuration for maximum corrosion plus wear resistance), simulation predicts the combined stress state from both processes and verifies that the total stress remains within acceptable limits.

Thermal Cycle Compatibility: Simulation verifies that the thermal cycle imposed by SHS overlay does not adversely affect the metallurgical bond quality of the underlying explosion-welded interface. This includes checking that peak temperatures remain below the bonding degradation threshold and that cooling rates do not induce cracking at the explosion weld interface.

8. Contribution to Qualification Building, Product Delivery, and Customer Value

8.1 Qualification Building

The numerical simulation capability for SHS ceramic electrode weld overlay directly contributes to the company's qualification portfolio in the following ways:

WPS Development and Qualification: Simulation provides the thermal cycle documentation required for ASME Section IX Part Q WPS qualification. By generating predicted temperature histories for various parameter combinations, the company can systematically develop qualified procedures for different substrate materials, geometries, and service conditions with reduced reliance on physical coupon testing.

Procedure Qualification Records (PQR):strong> Simulation results support PQR documentation by providing quantitative data on dilution rates, thermal cycles, and residual stress levels. This strengthens the technical basis for procedure qualification and provides regulators and customers with detailed evidence of process understanding.

Special Process Qualification: For applications requiring special qualification (nuclear, aerospace, offshore), the simulation capability demonstrates the company's technical competence and provides the detailed analytical evidence required by classification societies and regulatory bodies.

Non-Conformance Prevention: By predicting potential failure modes (cracking, delamination, excessive dilution) before physical production, the company can design processes that inherently avoid non-conformances, strengthening its quality reputation and reducing warranty exposure.

8.2 Product Delivery

Reduced Development Cycle: Simulation enables virtual optimization of process parameters, reducing the number of physical trials required to develop a production-ready process. This accelerates time-to-delivery for custom cladding solutions.

Quality Consistency: Simulation-derived process windows provide clear parameter ranges that ensure consistent quality across production batches. Operators can be trained using simulation-based process understanding, reducing human error.

Complex Geometry Capability: Simulation extends the company's capability to handle complex geometries (curved surfaces, thin-walled vessels, large-diameter pipes) where physical trial-and-error is impractical or prohibitively expensive.

Predictive Maintenance Support: Simulation models can be used to predict the remaining service life of cladded components based on accumulated thermal and mechanical cycling, enabling condition-based maintenance recommendations to customers.

8.3 Customer Value

Engineering Confidence: Providing customers with detailed simulation reports demonstrating predicted performance, residual stress levels, and service life significantly increases confidence in the delivered product. This is particularly valuable for high-value, critical-service applications where component failure carries significant safety or economic consequences.

Customized Solutions: Simulation enables the development of tailored cladding solutions for specific customer applications. Rather than offering standard products, the company can design optimized cladding configurations (material selection, pass sequence, PWHT schedule) for each customer's unique service conditions.

Failure Analysis Support: When customers experience issues with previously delivered cladded components, the simulation capability enables rapid root cause analysis by recreating the process conditions and predicting the resulting stress state and microstructure. This accelerates corrective action and maintains customer trust.

Value Engineering: Simulation can identify opportunities to reduce material usage or simplify fabrication without compromising performance. For example, simulation might reveal that a thinner overlay build achieves the same service life as a thicker build, reducing material and labor costs for the customer.

9. Summary

Numerical simulation of temperature and residual stress fields in self-propagating ceramic electrode weld overlay represents a critical technical capability that enhances the company's position in the high-performance cladding market. By providing predictive insight into the complex thermal and mechanical behavior of SHS processes, simulation enables process optimization, quality assurance, WPS qualification, and customer engineering support that would be impractical to achieve through experimental methods alone.

The integration of this simulation capability across the company's three technology routes—TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding—creates a comprehensive technical platform for addressing the full spectrum of cladding requirements. From extreme wear resistance (SHS ceramic overlay) to corrosion resistance (explosion welding) to metallurgical bonding integrity (hydraulic explosive bonding), the simulation expertise provides the analytical foundation for delivering reliable, high-performance cladding solutions.

As the company continues to expand its capabilities, the simulation infrastructure will serve as a strategic asset for qualification building, product development, and customer value creation, positioning the company as a technically differentiated provider in the global cladding technology market.