Three-Dimensional Temperature Field Numerical Simulation of Laser Surfacing on 40Cr Steel

Literature Overview

This 2012 paper in Heat Treatment of Metals by Li Gang and colleagues from Liaoning Technical University presents a three-dimensional transient finite element analysis of the temperature field during laser surfacing on 40Cr steel plates. The simulation was conducted using the ANSYS platform with APDL (ANSYS Parametric Design Language) to implement the moving heat source, providing quantitative predictions of melt pool temperature evolution and the effects of laser power and scanning speed on the thermal cycle.

Numerical Modeling Approach

The finite element model was constructed based on the following fundamental assumptions and formulations:

Governing Equations

The transient heat conduction equation with a moving heat source is:

ρc(∂T/∂t) = ∇·(k∇T) + Q(x, y, z, t)

Where ρ is density, c is specific heat, T is temperature, k is thermal conductivity, and Q represents the volumetric heat source distribution.

Heat Source Model

The laser heat source was modeled as a Gaussian surface heat flux distribution:

q(r) = (2P/πR²) × exp(-2r²/R²)

Where P is laser power, R is the effective beam radius, and r is the radial distance from the beam center. The heat source moves at velocity v along the x-direction, creating the characteristic cometary temperature field distribution.

Material Properties

40Cr steel was modeled with temperature-dependent thermal properties:

Simulation Results Analysis

Temperature Field Distribution Characteristics

The simulation revealed that the temperature field exhibits a characteristic cometary shape during laser scanning:

• Leading edge: High temperature gradient (steep thermal gradient), indicating rapid heating of the material ahead of the beam.
• Trailing edge: Lower temperature gradient (gradual cooling), representing the heat-affected zone behind the beam.
• Peak temperature location: Slightly ahead of the beam center due to the moving heat source effect.

Effect of Laser Power

At a constant scanning speed of 6 mm/s:

The near-linear relationship between laser power and peak temperature indicates that the thermal input is the dominant factor controlling the melt pool temperature under constant scanning speed conditions.

Effect of Scanning Speed

At a constant laser power of 800 W:

The inverse relationship between scanning speed and peak temperature demonstrates that reduced dwell time allows less heat accumulation, resulting in lower maximum temperatures and potentially shallower melt pools.

Thermal Cycle and Microstructural Implications

The temperature field simulation provides critical information for predicting the microstructural evolution during laser surfacing:

1. Solidification rate: Determined by the cooling rate at the solidification temperature (approximately 1400°C for 40Cr steel), which directly influences grain size and phase formation.
2. Heat-affected zone extent: Defined by the 800-900°C isotherm contour, indicating the region where microstructural changes occur without melting.
3. Thermal stress distribution: Derived from the temperature gradient field, predicting the magnitude and direction of residual stresses.
4. Dilution estimation: The depth and width of the melt pool determine the degree of substrate dilution of the surfacing alloy.

Process Optimization Recommendations

Based on the simulation results, the following process parameter recommendations emerge for laser surfacing of 40Cr steel:

• Power range: 600-1000 W provides a practical operating window for surface hardening and thin-layer surfacing.
• Scanning speed range: 4-10 mm/s offers sufficient flexibility for different application requirements.
• Thermal input control: Maintaining thermal input between 100-150 J/mm ensures adequate melting without excessive heat-affected zone.
• Multi-pass strategy: For thicker surfacing layers, overlapping passes with 50-70% overlap and controlled inter-pass temperature are recommended.

Comparison with Conventional Surfacing Processes

Study Insights and Engineering Value

The numerical simulation approach demonstrated in this study provides a powerful tool for process development and optimization without the cost and time associated with extensive experimental trials. The ability to predict temperature fields, cooling rates, and heat-affected zone dimensions enables informed decisions regarding:

• Powder selection based on predicted solidification conditions
• Preheating requirements to minimize thermal cracking risk
• Post-weld heat treatment parameters for residual stress relief
• Maximum achievable surfacing layer thickness before distortion becomes unacceptable

The simulation also highlights the importance of the cometary temperature field shape, which has implications for the directional solidification pattern and resultant microstructural anisotropy. Understanding this anisotropy is essential for predicting the mechanical behavior of the surfacing layer under service loading conditions. For pipeline applications, where the surfacing layer must withstand internal pressure and external mechanical loading, the directional properties of the microstructure directly influence the fatigue and fracture resistance of the repaired component.