Dynamic Simulation of Weld Overlay Thermal Stress Based on ANSYS Platform
1. Definition and Fundamental Principles
Dynamic simulation of weld overlay thermal stress refers to the computational modeling and numerical analysis of transient thermal-mechanical phenomena occurring during the application of cladding layers via welding processes. Using the ANSYS finite element analysis (FEA) platform, engineers can predict residual stress distributions, thermal strain evolution, distortion patterns, and cracking susceptibility that develop during weld overlay operations on base materials such as carbon steel, low-alloy steel, austenitic stainless steel, and nickel-based alloys.
The fundamental physics governing this simulation encompasses coupled thermo-mechanical analysis, where the transient heat transfer equation:
ρ·cₚ·(∂T/∂t) = ∇·(k·∇T) + Q
is solved sequentially or simultaneously with the elasto-plastic mechanical equilibrium equation:
∇·σ + f_b = 0
Here, ρ is density, cₚ is specific heat capacity, k is thermal conductivity, Q is the heat source intensity (typically modeled using Goldak's double-ellipsoidal or Gaussian heat source), σ is the stress tensor, and f_b represents body forces. The coupling between thermal and mechanical fields occurs through thermal expansion strain (ε_th = α·ΔT) and temperature-dependent material properties including yield strength, elastic modulus, and plastic flow behavior.
2. Category and Business Positioning
This technical capability falls within the company's engineering design and qualification support infrastructure. It serves as a critical enabler across all three technology routes—TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding—by providing predictive insight into process outcomes before physical trials are conducted. The positioning is as follows:
- Engineering Design Support: Enables virtual prototyping of overlay schemes, reducing the number of physical WPS trials required for qualification
- Quality Assurance: Predicts residual stress fields that could lead to stress corrosion cracking (SCC), fatigue failure, or dimensional distortion in finished clad products
- Customer Value Engineering: Provides quantitative data to justify design choices, optimize welding sequences, and demonstrate compliance with acceptance criteria to end-users
- IP Development: Contributes to proprietary process knowledge databases and supports patent filings for optimized welding procedures
3. Technical Purpose and Value
3.1 Residual Stress Prediction and Mitigation
Weld overlay processes generate severe thermal gradients (typically 500–1500°C/mm near the weld pool), producing residual stresses that can reach 300–600 MPa in the heat-affected zone (HAZ) and cladding layer. These stresses, if uncontrolled, compromise:
- Resistance to stress corrosion cracking in chloride-containing environments
- Fatigue life under cyclic loading conditions
- Dimensional accuracy for precision-machined clad components
- Interfacial bonding integrity in explosion-welded laminates subjected to subsequent welding
3.2 Process Optimization
ANSYS-based simulation allows engineers to systematically evaluate:
- Optimal welding sequence (back-step, skip-weld, multi-pass strategies) to minimize distortion
- Preheat temperature requirements for thick-section components
- Interpass temperature control windows
- Post-weld stress relief (PWSR) parameters—temperature, ramp rate, hold duration
- Effect of backing plate rigidity and fixture design on constraint-induced stress levels
3.3 Qualification Efficiency
By predicting process outcomes computationally, the company can reduce the number of physical WPS qualification trials from 3–5 down to 1–2, significantly reducing qualification cost and lead time while maintaining compliance with applicable standards.
4. Key Process and Implementation Points
4.1 ANSYS Simulation Workflow
- Geometry Modeling: Create 2D or 3D models of the base material, cladding layers, and fixtures. Simplifications include symmetric boundary conditions and shell elements for thin cladding layers.
- Material Property Definition: Input temperature-dependent properties for both base material and overlay metal including thermal conductivity, specific heat, Young's modulus, yield strength, and thermal expansion coefficient.
- Heat Source Modeling: Implement Goldak's double-ellipsoidal heat source or Gaussian surface heat flux, calibrated to measured bead geometry and welding parameters.
- Boundary Conditions: Apply convection and radiation heat transfer on exposed surfaces; apply displacement constraints representing fixture rigidity.
- Sequential Coupled Analysis: First solve the transient thermal problem, then map temperature history to the mechanical analysis as a body load.
- Plastic Strain Accumulation: Implement the "birth and death" element technique to simulate progressive weld bead deposition with appropriate plastic strain relaxation.
- Post-Processing: Extract residual stress distributions (σ_x, σ_y, σ_z), distortion profiles, and thermal cycle histories at critical locations.
4.2 Key Modeling Parameters
| Parameter | Typical Range | Influence on Results |
|---|---|---|
| Welding Current (TIG) | 100–350 A | Determines heat input and penetration depth |
| Travel Speed | 50–200 mm/min | Affects thermal cycle rate and HAZ width |
| Heat Input | 0.5–8.0 kJ/mm | Primary driver of residual stress magnitude |
| Preheat Temperature | 0–350°C | Reduces thermal gradient and cracking risk |
| Interpass Temperature | 50–250°C | Controls thermal cycling between passes |
| Fixture Constraint Factor | 0.2–1.0 (relative) | Higher constraint increases residual stress |
| Element Size (near weld) | 0.5–2.0 mm | Controls stress gradient resolution |
4.3 Heat Source Calibration
Accurate simulation requires calibration of the virtual heat source against experimental data. The company employs the following calibration methodology:
- Thermocouple arrays: Place 8–16 K-type thermocouples at defined offsets from the weld centerline on coupon specimens
- Peak temperature matching: Adjust heat source efficiency (typically 60–80% for TIG, 70–85% for MIG) to match simulated peak temperatures to measured values
- Thermal cycle rate matching: Verify cooling rates at 800→600°C and 600→400°C match experimental data within ±15%
- Weld bead geometry validation: Compare simulated fusion zone dimensions with macrographically measured bead profiles
4.4 Multi-Pass Simulation Strategy
For multi-layer multi-pass overlay schemes (common in TIG/MIG cladding with 3–12 passes), the simulation must account for thermal history accumulation. Key considerations include:
- Element activation timing synchronized with actual welding sequence
- Residual stress from previous passes serves as initial condition for subsequent passes
- Thermal softening of previously deposited material during subsequent passes
- Progressive constraint evolution as the workpiece accumulates cladding layers
5. Applicable Standards and Acceptance Criteria
5.1 Standards Referenced in Simulation Validation
| Standard | Relevance to Simulation |
|---|---|
| ASME Boiler and Pressure Vessel Code, Section IX | WPS/PQR qualification requirements that simulation helps optimize |
| ASME Section VIII, Division 2 | Fracture mechanics criteria for residual stress acceptance |
| NB/T 47014—Qualification Test Procedure for Welding of Pressure Vessels | Chinese national standard for welding procedure qualification |
| GB/T 19420—Welding Procedure Specification for Clad Steel | Defines acceptable residual stress levels for clad products |
| ASTM A240/A240M | Material specifications for stainless steel overlay materials |
| ASTM E112 | Grain size measurement for HAZ characterization validation |
| NACE MR0175/ISO 15156 | H₂S service requirements—simulation predicts cracking susceptibility |
| API 579-1/ASME FFS-1 | Fitness-for-service assessment using predicted residual stress fields |
| EN ISO 17640 | Welding—Welding procedure qualification—General rules |
| GB/T 985.1 | Welding procedure qualification test methods |
5.2 Simulation Output Acceptance Criteria
- Residual stress magnitude: Predicted peak residual stress in cladding layer shall not exceed 0.5×UTS of overlay material without stress relief
- Distortion: Predicted angular distortion shall not exceed 1.5 mm/m for flat plates, 0.5° for cylindrical components
- Cracking susceptibility: Predicted cooling rate at 600°C shall not exceed 15°C/s for low-alloy steels without preheat
- Thermal cycle compliance: Simulated peak temperature in base material HAZ shall not exceed 1100°C for carbon steels (to prevent grain coarsening)
- Interfacial stress: For explosion-welded laminates subjected to subsequent weld overlay, predicted interfacial shear stress shall remain below 0.6×yield strength of the softer layer
6. Common Risks and Controls
6.1 Modeling Risks
| Risk | Description | Control Measure |
|---|---|---|
| Over-simplified geometry | Neglecting fixture effects or component thickness variation | Include representative fixtures; use 3D models for thick sections | Inaccurate material properties | Using room-temperature properties at elevated temperatures | Employ temperature-dependent property curves from literature or testing | Heat source mis-calibration | Virtual heat input not matching actual energy delivery | Systematic calibration with thermocouple arrays and bead geometry | Boundary condition errors | Incorrectly representing thermal dissipation to backing plates | Model backing plates explicitly or apply measured heat flux BCs |
| Plastic strain relaxation | Inaccurate representation of thermal softening during multi-pass | Implement proper element birth/death with stress-free reference temperature |
| Mesh sensitivity | Results dependent on element size near weld zone | Perform mesh convergence study; minimum 3 elements across weld width |
6.2 Process Risks Identified Through Simulation
- Hydrogen-induced cracking: Simulation identifies regions with high residual tensile stress combined with high cooling rates, flagging areas requiring preheat or post-weld bake-out
- Hot cracking in overlay metal: Predicted solidification cracking susceptibility in dilution-sensitive overlay alloys (e.g., Ni-base alloys on steel) through analysis of thermal gradient and strain rate
- Delamination in explosion-welded laminates: When weld overlay is applied to explosion-welded cladding, simulation predicts stress concentrations at the explosion weld interface that could cause interfacial delamination
- Dimensional distortion: Quantifies out-of-plane warpage for large flat clad plates, enabling fixture design optimization
7. Application Across the Company's Three Technology Routes
7.1 TIG/MIG Weld Overlay Applications
- Multi-pass overlay design: Simulate 3–12 pass overlay schemes on pipe, plate, and vessel components to optimize welding sequence and minimize residual stress. Example: Back-step welding sequence simulation for 304L overlay on Q345R carbon steel pipe per GB/T 19420.
- Transition layer optimization: Predict dilution effects and residual stress in multi-layer transition schemes (e.g., 309L transition → 316L overlay → 625Ni overlay) to ensure compatibility and minimize cracking risk.
- Repair weld simulation: Model stress redistribution during local repair welding on previously clad surfaces, critical for maintaining existing cladding integrity.
- Post-weld stress relief (PWSR) design: Simulate the effect of stress relief heat treatment at 620–650°C to verify stress reduction to acceptable levels while avoiding sensitization in stainless steel overlays.
7.2 Hydraulic Explosive Bonding Applications
- Post-bonding weld overlay stress prediction: When hydraulic explosive bonding produces a laminate that subsequently receives weld overlay (e.g., for pipe cap welding), simulation predicts stress concentrations at the explosion bond interface under thermal cycling.
- Thermal compatibility analysis: Evaluate whether weld overlay thermal cycles could degrade the quality of the hydraulic explosive bond interface, particularly for dissimilar material combinations (e.g., aluminum on steel).
- Fixture design for bonded components: Model the combined effect of fixture rigidity and weld thermal input on bonded laminate components to prevent interfacial separation during subsequent fabrication operations.
7.3 Explosion Welding Applications
- Post-explosion welding repair simulation: Predict residual stress interaction between explosion welding residual stresses (typically 200–400 MPa compressive in cladding) and subsequent weld overlay repair stresses.
- Welding qualification on explosion-welded clad plate: Simulate welding procedures applied to explosion-welded clad plates (e.g., SS316L/CS, Cu/CS, Ti/CS) to verify that welding does not compromise the explosion bond interface quality.
- Thermal cycle impact on explosion weld interface: Analyze whether HAZ thermal cycles from overlay welding could cause intermetallic formation or interface weakening in the explosion weld bond zone.
- Large-scale component distortion: Predict distortion of large explosion-welded clad plates during subsequent edge welding, flange welding, and overlay welding operations.
8. Contribution to Qualification Building, Product Delivery, and Customer Value
8.1 Qualification Building
ANSYS-based thermal stress simulation directly accelerates the company's qualification portfolio development:
- Reduced trial cycles: Predictive simulation eliminates unsuccessful WPS attempts, reducing qualification time by 40–60%
- Parameter range definition: Simulation identifies safe parameter windows (preheat, interpass temperature, travel speed) that meet acceptance criteria, enabling broader qualification ranges
- Standard compliance demonstration: Provides quantitative evidence that welding procedures meet requirements of ASME Section IX, NB/T 47014, and EN ISO 17640
- Novel material system qualification: Enables qualification of new overlay material combinations without extensive trial-and-error, critical for expanding the company's certified material matrix
8.2 Product Delivery Enhancement
- Distortion control: Predicted distortion profiles enable pre-compensation in fixture design and machining allowances, improving first-pass acceptance rates
- Non-destructive testing (NDT) optimization: Identification of high-stress regions guides UT/MT inspection focus, improving defect detection efficiency
- Post-weld treatment planning: Simulation-optimized PWSR parameters ensure residual stress reduction without material degradation
- Batch consistency: Simulation-validated parameters enable repeatable production across multiple units of the same design
8.3 Customer Value Demonstration
- Engineering reports: Simulation results provide customers with quantitative residual stress predictions, supporting their fitness-for-service assessments per API 579-1/ASME FFS-1
- Service life prediction: Residual stress data feeds into fatigue and corrosion-fatigue life models, enabling customers to justify inspection intervals
- Compliance documentation: Simulation-backed qualification packages satisfy regulatory requirements for pressure vessel and nuclear component applications per NB/T 47014 and GB/T 985.1
- Value engineering: Demonstrates optimal cladding thickness and material selection through stress-life analysis, potentially reducing material costs while maintaining performance
9. Advanced Simulation Capabilities and Future Development
9.1 Constitutive Model Enhancement
Advanced simulations employ sophisticated constitutive models to improve prediction accuracy:
- Chaboche kinematic hardening model: Captures cyclic plasticity and ratcheting behavior in multi-pass welding
- Temperature-dependent isotropic hardening: Accounts for thermal softening and recovery during multi-pass sequences
- Coupled damage mechanics: Predicts micro-cracking initiation in overlay metal under combined thermal and mechanical loading
- Phase transformation modeling: Incorporates solid-state phase transformations in low-alloy steel HAZ (ferrite-to-austenite and reverse) affecting volumetric strain
9.2 Integration with Digital Twin Framework
The company is developing a digital twin framework where ANSYS simulation results serve as the baseline model, updated with real-time process monitoring data (welding current, voltage, travel speed, thermocouple readings) to enable:
- In-process residual stress prediction with ±20% accuracy
- Real-time process adjustment recommendations
- Post-production stress state verification through comparison with X-ray diffraction (XRD) or neutron diffraction measurements
9.3 Machine Learning Augmentation
Surrogate models trained on ANSYS simulation databases enable rapid (seconds vs. hours) prediction of residual stress for parameter variations, facilitating real-time process optimization during production. The simulation database serves as the training ground for these AI-driven process control systems.
10. Conclusion
ANSYS-based dynamic simulation of weld overlay thermal stress represents a foundational engineering capability that permeates all aspects of the company's cladding technology operations. From accelerating WPS qualification under standards such as ASME Section IX, NB/T 47014, and GB/T 985.1, to optimizing multi-pass overlay sequences for TIG/MIG welding, predicting stress interactions in post-explosion-welding fabrication, and ensuring dimensional accuracy of hydraulic explosive bonded laminates subjected to subsequent welding, this capability delivers measurable value across the entire product lifecycle. The systematic approach to simulation development—calibrated against experimental data, validated through post-weld stress measurements, and continuously refined through production feedback—ensures that the company's simulation-based engineering decisions are reliable, repeatable, and aligned with international quality management standards including ISO 9001 and ASME NQA-1 requirements for nuclear applications.