## Abstract

This study contains mathematical modelling and numerical analysis of heat transfer in laser beam welding process. The temperature field was obtained on the basis of numerical solution into unsteady heat transfer equation with convective term and volumetric heat sources taken into account. Volumetric heat source model describing laser beam power distribution in combined truncated cone–cylinder volume was developed. Due to the wide range of temperatures appearing in the process latent heat of fusion, evaporation as well as latent heat of phase transformations in solid state were taken into account in the solution algorithm. On the basis of developed numerical algorithms an analysis of heat transfer in laser butt-welded steel sheets as a three-dimensional initial-boundary problem was performed.

## Introduction

Laser beam welding has increasing application in many industries. This welding method allows joining materials previously regarded as a non-weldable or difficulty weldable [1–4]. Because of technological advantages, this technology successfully displaces conventional welding methods in many applications. However, laser welding of steel is accompanied by many phenomena not found in conventional welding [2]. Absorbed laser light vaporise material and forms a ‘keyhole’ that contains ionised vapour. The laser power is absorbed in the ionised vapour and transferred to the walls of the ‘keyhole’ forming the weld pool [3–5]. Moreover, the interaction of concentrated heat source on the material is responsible for generating high temperature gradients and high cooling rates appearing in the joint. This leads to occurring of hardening structures in the weld and heat-affected zone (HAZ) even in susceptible to welding, unalloyed and low-alloy constructional steels [6]. The quality of the joint is determined by process parameters such as: laser mode, beam radius, laser power and welding speed. Therefore, the knowledge about temperature profile of laser welded joints allows for a prediction of these parameters in terms of size and shape of the weld [7–9].

Theoretical analysis of laser-material interaction is very complex and includes various physical phenomena such as evaporation, melting, solidification, microstructure transformation, thermal and structural stress etc. [10–16]. Mathematical and numerical models usually describe chosen phenomena, due to the computational complexity of the problem. Therefore, coupled thermal phenomena and phase transformations in solid state are usually ignored in numerical analysis of welding processes [5, 8, 9, 15, 16]. However, latent heat generated during phase transformations in heat treatment of steel influences temperature distribution and has a significant effect on numerically predicted microstructure composition [17, 18].

A very important issue in modelling of laser beam welding is an appropriate selection of heat source power distribution which is primarily responsible for the melted pool shape and temperature distribution in the joint. Researchers are constantly looking for new mathematical models describing the distribution of the energy flux (taking into account the real process conditions) through theoretical study on laser induced plasma formation [19–21] or assuming Gaussian-like simplified volumetric heat source model [5, 7, 13] with decreasing energy intensity in material penetration direction, especially when developed theoretical model requires consideration of entire 3D workpiece geometry during calculations. A major problem in developing of simplified models is the determination of the size and shape of the heat source along the thickness of the workpiece corresponding to the experimental data with appropriate accuracy. Therefore, it is difficult to define a universal heat source model. From the analysis of laser beam welding process it is observed that heat source power decreases with increasing depth of penetration. Furthermore, obtained weld shape in deep penetration welding is often a mushroom-shape, mainly depending on the welding speed [4, 8].

In this study, theoretical model is presented to analyse coupled thermal phenomena and phase transformations in solid state during laser butt-welding of steel. Temperature field is calculated on the basis of numerical solution into heat transfer equation in Eulerian coordinates using finite difference method in forward Euler's time integration scheme. A heat source model is developed for the laser beam deep penetration welding, allowing better approximation of heat transfer in welded workpiece. Latent heat associated with material's state change and latent heat generated during phase transformations in solid state are taken into account in solution algorithms. Dilatometric research of S460 steel is performed to determine continuous-cooling-transformation diagram (CCT) and final fractions of structure constituents used in the model. Results of calculations are verified by the experiment. It is shown that the effect of latent heat associated with microstructure changes on temperature history in laser welding is considerable and should not be neglected in theoretical models.

Finally, it is demonstrated that developed heat source model gives a good approximation of fusion zone and HAZ geometry for mushroom-shape welds and well agrees with experimental result.

## Theoretical model

Temperature distribution in butt-welding is determined by process parameters, mainly distribution of laser beam power and welding speed. Phase transformations due to melting, solidification and evaporation as well as phase transformations in solid state are taken into account in theoretical investigations. At boundaries of the workpiece heat loss due to convection, radiation and evaporation is assumed. Sketch of considered system is illustrated in Fig. 1.

### Heat transfer

*λ*=

*λ*(

*T*) is thermal conductivity dependent on temperature,

*C*

_{ef}is effective heat capacity, is a volumetric heat source and

**v**=

**v**(x,y,z,

*t*) is a vector determined by welding speed,

*T*(x,y,z,

*t*) is a temperature and

*t*is time.

*α*is a convective coefficient,

*∊*is radiation coefficient,

*σ*is Stefan-Boltzmann constant and

*q*=

*q*(

*r*,0) is the heat flux towards the top surface of welded workpiece (

*z*= 0) in the source activity zone.

Thermal conductivity coefficient changing with temperature is assumed in theoretical model (Fig. 2). In solid state *λ* = *λ*(*T*) is defined according to data form the literature [22]. Much higher value of *λ*(*T*) is assumed at high temperatures, which corresponds to the motion of liquid material in the welding pool.

#### Thermal conductivity assumed in theoretical model

Citation: Journal of Thermal Analysis and Calorimetry J Therm Anal Calorim 110, 1; 10.1007/s10973-012-2486-0

#### Thermal conductivity assumed in theoretical model

Citation: Journal of Thermal Analysis and Calorimetry J Therm Anal Calorim 110, 1; 10.1007/s10973-012-2486-0

#### Thermal conductivity assumed in theoretical model

Citation: Journal of Thermal Analysis and Calorimetry J Therm Anal Calorim 110, 1; 10.1007/s10973-012-2486-0

### Laser beam heat source

*α*

_{z}

*d*heat source power is decreasing in a shape of a truncated cone. Remaining laser power distribution is constant in a shape of a cylinder. Heat source power distribution in truncated cone volume is expressed as follows [13]:

*Q*

_{L}is a laser beam power,

*η*is efficiency,

*d*is a total penetration depth, is a current beam radius and is a current penetration,

*V*is a volume of the heat source.

*r*

_{0}is a radius of laser beam spot at the top surface of the workpiece (

*z*= 0) and is a radius of the beam at the deep of penetration is experimentally defined coefficient depending on the welding speed.

#### Laser beam heat source power distribution **a** at the *top* and *bottom* surfaces of the joint and **b** in *cross-section* of the joint. Heat source parameters are: *Q*_{L} = 3,800 W, *η* = 75 %, *r*_{0} = 1 mm, *r*_{αz} = 0.4 mm, *α*_{z} = 0.5

#### Laser beam heat source power distribution **a** at the *top* and *bottom* surfaces of the joint and **b** in *cross-section* of the joint. Heat source parameters are: *Q*_{L} = 3,800 W, *η* = 75 %, *r*_{0} = 1 mm, *r*_{αz} = 0.4 mm, *α*_{z} = 0.5

#### Laser beam heat source power distribution **a** at the *top* and *bottom* surfaces of the joint and **b** in *cross-section* of the joint. Heat source parameters are: *Q*_{L} = 3,800 W, *η* = 75 %, *r*_{0} = 1 mm, *r*_{αz} = 0.4 mm, *α*_{z} = 0.5

### Phase transformations

*T*

_{S}) and liquidus (

*T*

_{L}) temperatures. Latent heat of fusion is defined in a body flux [23]

*ρ*

_{S}is solidus density,

*H*

_{L}is a latent heat of fusion and is a function determining the fraction of solid phase in the mushy zone.

*c*is a specific heat. The product of density and specific heat in the mushy zone depends linearly on solid fraction .

*T*

_{max}) and a full equilibrium of metal vapour pressure in the keyhole and pressure of shielding gases, effective heat capacity in temperatures exceeding boiling point is defined as follows [15]:

*i*-th phase transformation, is a volumetric fraction of

*i*-th structural constituent.

*b*= and

*n*= are determined by start (

*η*

_{s}= 0.01) and final (

*η*

_{f}= 0.99) conditions for phase transformation as well as start and final times of phase transformation (

*t*

_{s}and

*t*

_{f}) estimated using CCT diagram [29] for welded steel

*t*is a time.

*ηM*) is estimated using Koistinen–Marburger (KM) equation

*k*depends on martensite phase start and final temperatures (

*M*

_{s}and

*M*

_{f}) determined also basing on CCT diagram

## Results and discussion

The mathematical model of thermal phenomena accompanying welding processes was implemented into the computer solver. Governing equation (1) was numerically solved using finite difference method in forward Euler's time integration scheme. For the stabilization of solution algorithms Péclet number was implemented directly into central difference quotients. Numerical analysis of laser but-welding without additional material is performed as a three-dimensional problem (Fig. 4) assuming perfect contact between joined sheets. Analyzed domain with dimensions: *L* = 250 mm, *a* = 50 mm and *g* = 5 mm is considered as a half of the joint to the plane of symmetry and discretized by differential grid with constant spatial step 0.01 mm in *x* direction and linearly decreasing spatial step in *y* direction (Fig. 4b) from 0.05 mm in the heat source activity zone, where large temperature gradients are present, up to 2 mm in lower temperatures. The average time of calculations for 25 s of welding simulation lasted about 16 h. Calculations were performed on PC with Quad Core 2.40 GHz.

### Sketch of **a** analysed initial-boundary problem and **b** finite difference mesh used in calculations

### Sketch of **a** analysed initial-boundary problem and **b** finite difference mesh used in calculations

### Sketch of **a** analysed initial-boundary problem and **b** finite difference mesh used in calculations

Illustrated in Fig. 3 heat source power distribution and welding speed *v* = 0.7 m min^{−1} are used in calculations. The analysis is carried out for sheets made of S460 steel with chemical composition: 0.19 C, 1.62 Mn, 0.6 Si, 0.10 Cr, 0.09 Ni, 0.11 Cu, 0.10 V, 0.032 Al, 0.012 Nb (wt%). Thermo-physical properties used in calculations are summarized in Table 1.

Thermo-physical properties assumed in computer simulations

Thermo-physical property | Value | |
---|---|---|

Solidus temperature | T_{S}/K | 1,750 |

Liquidus temperature | T_{L}/K | 1,800 |

Boiling point | T_{b}/K | 3,010 |

Ambient temperature | T_{0}/K | 293 |

Specific heat of solid phase | c_{S}/J kg^{−1}K^{−1} | 650 |

Specific heat of liquid phase | c_{L}/J kg^{−1}K^{−1} | 840 |

Density of solid phase | ρ_{S}/kg m^{−3} | 7,800 |

Density of liquid phase | ρ_{L}/kg m^{−3} | 6,800 |

Latent heat of fusion | H_{L}/J kg^{−1} | 270 × 10^{3} |

Latent heat of evaporation | H_{b}/J kg^{−1} | 76 × 10^{5} |

Latent heat of phase transformations | H_{A→F}/J kg^{−1} | 8 × 10^{4} |

H_{A→P}/J kg^{−1} | 9 × 10^{4} | |

H_{A→B}=H_{A→M}/J kg^{−1} | 11.5 × 10^{4} | |

Convective heat transfer coefficient | α/W m^{−2}K^{−1} | 50 |

Boltzmann's constant | σ/W m^{−2}K^{−4} | 5.67 × 10^{−8} |

Surface radiation emissivity | ∊ | 0.5 |

Used in the analysis CCT diagram of S460 steel and final fractions of structure constituents were obtained on the basis of dilatometric research performed at Institute for Ferrous Metallurgy in Gliwice (Poland) using dilatometer DIL805 produced by Bähr Thermoanalyse GmbH, equipped with a LVDT measuring head, with a theoretical resolution of 0.057 μm. In dilatometric analysis a constant austenitization temperature *T*_{A} = 1,200 °C was assumed and different heating and cooling rates, simulating thermal cycles in welding. The microstructure of analyzed samples was afterwards evaluated by microstructure analysis, supported by measurements of microhardness. Obtained CCT diagram with final fractions of structure constituents is illustrated in Fig. 5.

### CCT diagram of S460 steel and final fractions of microstructure constituents

### CCT diagram of S460 steel and final fractions of microstructure constituents

### CCT diagram of S460 steel and final fractions of microstructure constituents

Figure 6 shows temperature distribution for chosen points at the top surface of the joint (*z* = 0) and in a distance *y* = 1 mm from the centre of laser beam heat source. Solidus, liquidus and boiling temperatures are marked as well as temperature range [800 °C; 500 °C], where cooling rates are determined for analysis of phase transformations in solid state. It can be observed that maximum temperatures are obtained in the centre of heat source activity zone. At this point temperature exceeds the boiling point of steel (*T*_{b} = 3,010 K), leading to the evaporation and creation of the keyhole. If maximum temperature exceeds solidus temperature (*T*_{S}) material is melted. Complete melting is achieved when temperature reaches liquidus temperature (*T*_{L}). When maximum temperature exceeds austenitization temperature (), the material is in the zone of structural transformations (HAZ). Below austenitization temperature welded material is treated as untransformed base material.

### Temperature distribution for chosen points: **a** at the *top* surface of the weld and **b** along direction of material penetration, in a distance *y* = 1 mm from the centre of laser beam heat source

### Temperature distribution for chosen points: **a** at the *top* surface of the weld and **b** along direction of material penetration, in a distance *y* = 1 mm from the centre of laser beam heat source

### Temperature distribution for chosen points: **a** at the *top* surface of the weld and **b** along direction of material penetration, in a distance *y* = 1 mm from the centre of laser beam heat source

Effective heat capacity in solid state for chosen points at the top surface of the joint is shown in Fig. 7. It can be noticed that various effective heat capacity is obtained depending on the kinetics of phase transformations in solid state and different start and final temperatures determined for each cooling rate in the weld and HAZ. It can be observed that for presented in Fig. 7 material points the effective heat capacity is determined by bainite and martensite transformation. Fig. 8 shows temperature history of a point on the top surface of the joint at a distance 1 mm from the weld line. The solid line represents temperature distribution affected by latent heat included into the model, whereas dotted line represents temperature distribution calculated without the effect of latent heat (latent heat of phase transformations in solid state was omitted in the model). It is observed that the impact of latent heat of phase transformations in solid state on temperature distribution is measurable.

### Effective heat capacity for kinetics of phase transformations calculated for two chosen thermal cycles

### Effective heat capacity for kinetics of phase transformations calculated for two chosen thermal cycles

### Effective heat capacity for kinetics of phase transformations calculated for two chosen thermal cycles

### Temperature history in a chosen point of the joint showing the influence of latent heat

### Temperature history in a chosen point of the joint showing the influence of latent heat

### Temperature history in a chosen point of the joint showing the influence of latent heat

Figure 9 illustrates temperature field at the top surface of the workpiece (from the face of the weld). The welding pool boundary is represented by a black solid line (solidus isotherm), whereas dashed line defines heat affected zone boundary (austenitization temperature *T*_{g} = 1,000 K).

### Temperature field at the *top* surface of the joint (*z* = 0, form the face of the weld)

### Temperature field at the *top* surface of the joint (*z* = 0, form the face of the weld)

### Temperature field at the *top* surface of the joint (*z* = 0, form the face of the weld)

The formation of laser butt-welded joint is presented in Fig. 10. This figure illustrates temperature distribution at different analysis time. At time *t* = 1.72 s maximum temperatures are achieved in the joint. For *t* = 1.85 s maximum size of melted zone is achieved at the top surface of the joint (*z* = 0). Maximum size of heat affected zone is reached at time *t* = 2.47 s.

### Temperature distribution in the cross-section of the weld for different simulation times

### Temperature distribution in the cross-section of the weld for different simulation times

### Temperature distribution in the cross-section of the weld for different simulation times

In order to verify the correctness of elaborated models, CO_{2} laser welding experiment was performed at Welding Institute in Gliwice. TRUMPH Lasercell 1005 with CO_{2} laser generator (*Q*_{L} = 3,800 W) was used in the experiment. Sheets made of S460 steel with dimensions: 250 × 50 × 5 mm were arranged along the longest side without a gap (perfect contact). Welding speed was set to 0.7 m min^{−1}.

Figure 11 shows the comparison of experimentally obtained weld with numerically estimated melted zone and heat affected zone geometries. As shown in this figure, the calculated geometry of the weld and HAZ agrees well with the experimental result.

### Comparison of numerically predicted melted zone and HAZ geometry with weld obtained in the experiment

### Comparison of numerically predicted melted zone and HAZ geometry with weld obtained in the experiment

### Comparison of numerically predicted melted zone and HAZ geometry with weld obtained in the experiment

## Conclusions

Three-dimensional unsteady heat transfer numerical model was developed for numerical analysis of temperature field in laser-welded sheets made of S460 steel. Volumetric heat source model was developed for a determination of laser beam power distribution in deep penetration welding. Latent heat associated with microstructure changes was included into the model. On the basis of calculated temperature field melted zone and HAZ geometry was determined. The following conclusions can be made from this study:

Estimated geometry of the weld and HAZ is very close to this obtained in the experiment, thus developed heat source model allows for a better approximation of temperature field in laser welded joints.

Theoretical model can be used not only for a prediction of joint geometry but also to analyze influence of process parameters on the size and shape of laser welded joints.

Heat generated during phase transformations in solid state included into theoretical model allows for a better approximation of temperature distribution in welded joint.

## References

1. Steen, WM 1991

*Laser material processing*Springer London.2. Dawes, C 1992

*Laser welding*Abington Publishing New York .3. Jin, X, Li, L, Zhang, Y 2002 A study of fresnel absorption and reflections in the keyhole in deep penetration laser welding.

35:2304–2310 .*J Phys D*4. Pilarczyk, J, Banasik, M, Dworak, J, Stano, S 2006 Technological applications of laser beam welding and cutting at the Instytut Spawalnictwa.

5–6:6–10.*Przeglad Spawalnictwa*5. Han, L, Liou, FW 2004 Numerical investigation of the influence of laser beam mode on melt pool.

47:4385–4402 .*Int J Heat Mass Transf*6. Bokota, A, Piekarska, W 2008 Modeling of residual stresses in laser welding.

6:19–24.*Paton Weld J*7. De, A, DebRoy, T 2005 Reliable calculations of heat and fluid flow during conduction mode laser welding through optimization of uncertain parameters.

84:101–112.*Weld J*8. Gery, D, Long, H, Maropoulos, P 2005 Effects of welding speed, energy input and heat source distribution on temperature variations in butt joint welding.

167:393–401 .*J Mater Process Technol*9. Chang, WS, Na, SJ 2002 A study on the prediction of the laser weld shape with varying heat source equations and the thermal distortion of a small structure in micro-joining.

120:208–214 .*J Mater Process Technol*10. Dowden, JM 2001

*The mathematics of thermal modeling*Taylor & Francis Group USA .11. Belhadj, A, Bessrour, J, Bouhafs, M, Barrallier, L 2009 Experimental and theoretical cooling velocity profile inside laser welded metals using keyhole approximation and Boubaker polynomials expansion.

97:911–915 .*J Therm Anal Calorim*12. Abderrazaka, K, Ben Salemb, W, Mhiria, H, Lepalecc, G, Autric, M 2008 Modelling of CO2 laser welding of magnesium alloys.

40:581–588 .*Opt Laser Technol*13. Tsirkas, SA, Papanikos, P, Kermanidis, Th 2003 Numerical simulation of the laser welding process in butt-joint specimens.

134:59–69 .*J Mater Process Technol*14. Strezov, V, Lucas, JA, Strezov, L 2003 Computer aided thermal analysis.

72:907–918 .*J Therm Anal Calorim*15. Piekarska, W, Kubiak, M 2011 Three-dimensional model for numerical analysis of thermal phenomena in laser-arc hybrid welding process.

54 23–24 4966–4974 .*Int J Heat Mass Transf*16. Rai, R, Kelly, SM, Martukanitz, RP, Debroy, T 2008 A convective heat-transfer model for partial and full penetration keyhole mode laser welding of a structural steel.

39A:98–112 .*Metall Mater Trans A*17. Lee, SJ, Lee, YK 2009 Latent heat of martensitic transformation in a medium-carbon low-alloy steel.

60:1016–1019 .*Scripta Mater*18. Woodard, PR, Chandrasekar, S, Yang, HTY 1999 Analysis of temperature and microstructure in the quenching of steel cylinders.

30B:815–822 .*Metall Mater Trans B*19. Bedenko, DV, Kovalev, OB, Krivtsun, IV 2010 Simulation of plasma dynamics in a keyhole during laser welding of metal with deep penetration.

43:105501 .*J Phys D*20. Beck, M, Berger, P, Hugel, H 1995 The effect of plasma formation on beam focusing in deep penetration welding with CO2 lasers.

28:2430–2442 .*J Phys D*21. Semak, V, Matsunawa, A 1997 The role of recoil pressure in energy balance during laser materials processing.

30:2541–2552 .*J Phys D*22. Michaleris, P, DeBiccardi, A 1997 Prediction of welding distortion.

4:172–181.*Weld J*23. Nedjar, B 2002 An enthalpy-based finite element method for nonlinear heat problems involving phase change.

80:9–21 .*Comput Struct*24. Piekarska, W, Kubiak, M, Bokota, A 2011 Numerical simulation of thermal phenomena and phase transformations in laser-arc hybrid welded joints.

56 2 409–421 .*Arch Metall Mater*25. Lee, KJ 1999 Characteristics of heat generation during transformation in carbon steels.

40:735–742 .*Scripta Mater*26. Serajzadeh, S 2004 Modelling of temperature history and phase transformations during cooling of steel.

146:311–317 .*J Mater Process Technol*27. Zhang, W, Wood, B, DebRoy, T, Elmer, JW, Palmer, TA 2003 Kinetic modeling of phase transformations occurring in the HAZ of C-Mn steel welds based on direct observations.

51 12 3333–3349 .*Acta Mater*28. Piekarska W . Numerical analysis of thermomechanical phenomena during laser welding process.

*The temperature fields, phase transformations and stresses*. Monograph series no. 135. Czestochowa; 2007.29. Orlich J , Pietrzeniuk HJ.

. vol. 3. Dusseldorf; 1973.*Atlas zur Wärmebehandlung der Stähle*