Open access peer-reviewed chapter

The Analysis of Temporary Temperature Field and Phase Transformations in One-Side Butt-Welded Steel Flats

Written By

Jerzy Winczek

Submitted: November 5th, 2015 Reviewed: April 28th, 2016 Published: September 21st, 2016

DOI: 10.5772/63994

Chapter metrics overview

1,495 Chapter Downloads

View Full Metrics

Abstract

In this chapter, the welding method applied for modelling the temperature field and phase transformations is presented. Three-dimensional and temporary temperature field for butt welding with thorough penetration was determined on the basis of analytical methods of an integral transformation and Green’s function. Structural changes of heating and cooling, proceeding in a weld (in the heat-affected zone), were described using the existing formulations of phase transformations. Considerations were illustrated by an example, for which analysis of temperature fields, developed by a moving heat source, and calculations of the distribution of particular phases (structures) were carried out. Metallographic studies of the butt joints, which were arc welded under a flux, were carried out in the empirical part of this work. Their results enabled the verification of the numerical simulation results of the phase transformations.

Keywords

  • butt-welded joint
  • temperature field
  • HAZ
  • phase transformations
  • numerical modelling
  • metallographic examination

1. Introduction

Welding is characterised by many specific features associated with variable temperatures and variable physical and mechanical properties of the welding material. The moving heat source, characteristic of welding, partial melts the joint surface and fuses an electrode. The electrode fills a joint space with liquid metal. Hence, welding elements are subjected to varying temperature ranges, that is, from ambient to that of a liquid metal.

Crystallisation and solidification, segregation of alloy elements and solutes and structural changes caused by intensive cooling occur extensively. Thermal and mechanical states and microstructure directly state about the quality of the welding joint.

Modelling the temperature field during welding was first initiated by Rosenthal [1] and Rykalin [2], who supposed the point and linear models of heat source, respectively. The adoption of a point heat source, as in the above-mentioned studies, yields results with respect to the points located near the centre of the weld, which are significantly different from the actual temperature values. Therefore, Eagar and Tsai [3] proposed a two-dimensional (2D) Gaussian-distributed heat source model and developed a solution of temperature field in a semi-infinite steel plate. Subsequently, Goldak et al. [4] introduced a double ellipsoidal three-dimensional heat source model. There are two ways of modelling the temperature field during welding: analytical [514] and numerical (the finite difference methods, infinitesimal heat balances and finite element method) [1530]. The welding methods and types of joints can be studied through these approaches [6, 20, 21, 3133]. The construction of numerical models with heightening complexity allows more essential factors for the exact description of the structural changes in the welded steel.

Advertisement

2. Temperature field in the butt-welded joint with thorough penetration

Welding is characterised by an application of the movable, concentrated heat source, which in turn makes the temperature field movable in time and space:

T=T(r,t)=T(x,y,z,t)E1

Studies are being conducted to develop models of temperature field. Such models should have a real-time shape and temperature gradients based on the geometrical dimension of the welding element and also time. Referring to the formulated problem, the solution of heat equation for isotropic medium is essential to determine a temporary temperature field:

2T(r,t)=1atT(r,t)f(r,t)λE2
where T(r,t) is temperature at rposition at ttime, a is the coefficient of temperature compensation, λ thermal conductivity and f(r,t) supplied energy per volume and time unit.

Analytical method, proposed by Geissler and Bergmann [34, 35], was chosen to solve this differential equation. A short description of the method, described in detail in the above-mentioned studies, is presented below.

The following assumptions were accepted in the calculations:

  1. quantities characterising the material properties, such as thermal conductivity, temperature compensation and thermal capacity, are constant (independent from temperature),

  2. heat waste by convection and radiation is negligible,

  3. reciprocal interaction of temperature field and phase changes is not taken into account,

  4. heat of fusion is not taken into consideration.

A sample with thickness of Dand width of B1+ B2, heated by the movable welding heat source, which is displaced at velocity valong the x-axis (Figure 1), is an illustration of the below consideration.

Figure 1.

Schematic of heating of a steel sample using welding heat source.

According to Geissler and Bergmann [34, 35], the solution of Eq. (2) can be written as a superposition of Green’s function. This leads to the following convolution of integrals as a general expression of temperature

T(r,t)=aλ0tf(r',t')G(|rr'|,tt')dx'dy'dz'dt'E3
where r(x,y,z) is a vector pointing the place on the sample, while r′= (x′,y′,z′) determines the source position. Green’s function Gdescribes the temperature field in the point of material defined by rin time t. It is caused by a point heat source acting in the r′position and at t′< ttime. Gdepends on the geometry of the sample and can be determined by transformation method, while fdefines the cross-section of the welding heat source.

A three-dimensional temperature field with the possibility of acceptance of different geometries of samples as well as the shape of the heat source can be determined from Eq. (3).

In the case of the Gauss model of heat source, we have

f(x,y)=P2πR2exp(x2+y22R2),<x<,<y<,E4

where the power of source is denoted by Pand determined for Rradius and corresponds to 1/e of its peak value (Figure 2).

Figure 2.

Gauss distribution of power of the heat source.

An infinitely long bar with the above-mentioned dimensions of cross-section was accepted in the considered example (Figure 1).

This can be written in the Cartesian coordinate system as follows:

<x<,B1yB2,0zDE5

Boundary conditions defining the surface and Green’s function were taken from the study by Carslaw and Jaeger [36]. Green’s function takes the following form:

G(|rr'|,tt')=G(x)(xx',tt')G(y)(yy',tt')G(z)(z,tt')E6
G(x)=[4πa(tt')]1/2exp((x+v(t0t')x')24a(tt'))E7
G(y)=[4πa(tt')]1/2×[n=exp([yy'2n(B2+B1)]24a(tt'))+n=1exp([yy'2nB22(n1)B1]24a(tt'))+n=1exp([yy'+2nB1+2(n1)B2]24a(tt'))]E8
G(z)=[4πa(tt')]1/2n=exp((z2nD)24a(tt'))E9
where nis the transformation number of the source.

The relationship of movement of the welding heat source to the welding element is included in G(x) function. Because a flat model of heat source was assumed, the z′coordinate is not present in G(z) function. The coordinate systems connected with the heat source and welding material coincide with t0time. Considering that the shape of the heat source is independent of time, z′is dependent on fand Gis eliminated; a modified integral to count temperature profile is thus obtained:

T(r,t)=aλ0tf(x',y')G(|rr'|,tt')dx'dy'dt'E10
z′integral is removed in comparison to Eq. (3).

Temperature distribution can be calculated after substituting Green’s function and Eq. (4) into Eq. (10). The integral over a range of variables can be evaluated, which yields the following result:

T(x,t)=ηP2πρCp0tu(t,t')n=exp((z2nD)24a(tt'))×{n=exp((x+v(t0t')))2+(y4nB)22R2+4a(tt'))F1(y)+n=exp((x+v(t0t')))2+(y4nB)22R2+4a(tt'))F2(y)+n=exp((x+v(t0t')))2+(y4nB)22R2+4a(tt'))F3(y)}dt'+T0E11
where
u(t,t')=12πa(tt')(R2+2a(tt'))E12
F1(y)=erf(R2+2a(tt')4aR2(tt')(B2R2(y4nB)R2+2a(tt')))+erf(R2+2a(tt')4aR2(tt')(B1R2(y4nB)R2+2a(tt')))E13
F2(y)=erf(R2+2a(tt')4aR2(tt')(B2R2(y2(2n1)B)R2+2a(tt')))+erf(R2+2a(tt')4aR2(tt')(B1R2(y2(2n1)B)R2+2a(tt')))E14
F3(y)=erf(R2+2a(tt')4aR2(tt')(B2R2(y+2(2n1)B)R2+2a(tt')))+erf(R2+2a(tt')4aR2(tt')(B1R2(y+2(2n1)B)R2+2a(tt')))E15

Advertisement

3. Kinetics of phase transformations in a solid state

Heating processes of steel lead to the transformation of a primary structure into austenite, while cooling leads to the transformation of austenite into ferrite, pearlite, bainite and martensite. Structural changes of a welded joint, connected with its cooling (also with hardening), develop heterogeneous image of material structure, which influences the state of stress after welding. The zone with a yield point lesser or greater than that of an indigenous material can occur in the welded joint.

Mechanical properties of the joint mostly depend on the type of welding material (its primary structure and chemical constitution of steel) and the characteristics of heat cycles accompanying welding. Temperature levels attained during heating, the hold time at a particular temperature and velocity of cooling in the 800–500°C range determine the type of structure present in the joint during and after welding.

Figure 3.

Characteristic structural areas of a welded joint, depending on the temperature and the share of the carbon in the steel.

Figure 3 shows the distribution characteristics of the weld joint zone of structural carbon steel [37] with a schematic of the fragment of an iron-carbon system and fragment of the TTT-welding diagram. Together, it is categorised into the following zones:

  1. fusion zone, which undergoes a thorough penetration and is characterised by the dendritic structure of solidification,

  2. partial joint penetration, where material is in a semi-fluid state and creates the border between the melted material and the material being converted into austenite,

  3. the course-grained structure, the so-called overheating zone,

  4. proper transformation, where perfect conversion of primary structure into austenite occurs,

  5. partial transformation between temperature A1 at the beginning of austenisation and A3 at the end of austenisation, where only a part of the structure changes into austenite,

  6. recrystallisation.

Several studies have focused on the description and numerical modelling of steel phase transformations. These studies have been reviewed by Rhode and Jeppson [38].

The type of a newly created phase depends heavily on the kinetics of heating and cooling processes. Kinetics of those processes is described by Johnson-Mehl-Avrami’s and Kolomogorov’s (JMAK) rules [39]. The amount of austenite ϕAcreated while heating the ferrite-pearlitic steel is therefore defined according to the following formula:

ϕA(T)=jϕj0(1exp(bj(T)tnj(T)))E16
where ϕj0constitutes an initial share of ferrite (jF), pearlite (jP) and bainite (jB), while constants bjand njare determined using conditions of the beginning and the end of transformation:
nj=ln(ln(0.99))ln(A1/A3),bj=0.01niA1E17

In welding processes, the volume fractions of particular phases during cooling depend on the temperature, cooling rate, and the share of austenite (in the zone of incomplete conversion 0≤ϕA≤1). In a quantitative perspective, the progress of phase transformation during cooling is estimated using additivity rule by voluminal fraction ϕjof the created phase, which can be expressed analogically in Avrami’s formula [40] by equation:

ϕj(T,t)=ϕAϕjmax{1exp[bj(T(v8/5))t(T)n(T(v8/5))]}E18
where φjmaxis the maximum volumetric fraction of phase jfor the determined cooling rate estimated on the basis of the continuous cooling diagram (Figure 4), while the integral volumetric fraction equals to:
j=1kϕj=1E19
and kdenotes the number of structural participations.

The quantitative description of dependence of the material structure and quality on temperature and transformation time of overcooled austenite during surfacing is made using the time-temperature-transformation diagram during continuous cooling, which combines the time of cooling t8/5(time when material stays within the range of temperature between 500 and 800°C, or the velocity of cooling (v8/5– (800 – 500)/t8/5) and the temperature with the progress of phase transformation (Figure 4). Those diagrams are called TTT–welding diagrams.

Figure 4.

Scheme of phase changes of overcooled austenite depending on cooling velocity within temperature range 800–500°C.

Quantitatively, the progress of phase transformation is estimated by volumetric fraction ϕjof the created phase, where idenotes ferrite (jF), pearlite (jP), bainite (jB) or martensite (jM). The volumetric fraction ϕjof the created phase can be expressed using a formula given in Eq. (27), wherein time tis replaced with a new independent variable, temperature T[41]:

ϕj=ϕAϕjmax(1exp(bjTnj))+ϕj0E20
where
nj=ln(ln(1ϕjs)/ln(1ϕjf))ln(Tjs/Tjf),bj=nj(1ϕjf)TjsE21
ϕjsϕjmax=0,01,ϕjfϕjmax=0,99E22
ϕj0is the volumetric participation of j-th structural component, which has not been converted during austenitisation; Tjs= Tjs(v8/5) and Tjf= Tjf(v8/5) are, respectively, the initial and final temperature of phase transformation of this component.

The fraction of martensite formed below the temperature Msis calculated using the Koistinen-Marburger formula [42, 43]:

ϕM(T)=ϕA/ϕMmax{1exp[μ(MsT)]},μ=ln(ϕMmin=0.1))MsMfE23
where ϕmdenotes the volumetric fraction of martensite; Msand Mfdenote the initial and final temperature of martensite transformation, respectively; Tis the current temperature of the process.

Advertisement

4. Thermal and phase transformation strains

Changes in temperature during welding cause deformations associated with the thermal expansion and deformation of the material resulting from the structural phase transformation. Deformation during the whole thermal cycle is the total deformation created during heating and cooling [44]:

ε(x,y,z,t)=εH+εCE24
where εH and εCdenote the thermal and phase transformation strains during heating and cooling, respectively.

Heating leads to an increase in the material volume, while transformation of the initial structure (ferritic, pearlitic or bainitic) in austenite causes shrinkage which is associated with different densities of the given structures. Then, the strain caused during heating is calculated as follows:

εH=εThεTrhE25
where εThis the strain caused by thermal expansion of the material:
εTh=i=A,P,F,B,M(αiϕi0(TT0)H(TA1T)+αiϕi(TTA1)H(TA3T)H(TTA1)++αA(TTA3)H(TTA3))E26
while εTrhis the phase transformation strain during heating:
εTrh=i=P,F,B,MϕiγiAE27
where γiAis the structural strain of the i-th structure in austenite, T0is the initial temperature, αiis the linear thermal expansion coefficient of the i-th structure and H(x) is the function defined as follows:
H(x)={1forx>00,5forx=00forx<0E28

During cooling, the total strain (similarly as during heating) is the sum of strains associated with thermal expansion (in this case, the shrinkage of the material) as well as structural strains. Volumetric increase can be attributed to the high density of austenite (highest among the hardening structures such as martensite, bainite, ferrite and pearlite). The strain caused during cooling can be described by the following relation [44]:

εC=εTc+εTrcE29

where εTcis the strain caused by thermal shrinkage of material:

εTc=αA(TTSOL)H(TTs)+αA(TsTSOL)H(TsT)++i=A,P,F,B,Mαiϕi(TTsi)H(TsiT)E30

while εTrcis the strain caused by phase transformation during cooling:

εTrc=i=P,F,B,MϕiγAiE31

where TSOLdenotes solidus temperature, Tsthe initial temperature of phase transformation, Tsithe initial temperature of austenite transformation in the i-th structure and γAithe structural strain of austenite in the i-th structure. In addition, due to the limitation of the existence of solid state material:

ε(x,y,z,t)=0forT>TSOLE32
Advertisement

5. The example of calculation of temperature field and phase transformations in welded flats

In the considered example, it is assumed that the welded material is steel S235 with the following material constants: specific heat Cp= 670 J/(kg·K), density ρ= 7800 kg/m3, thermal diffusivity a= 1.2 × 10−5 m2/s, the temperatures at the beginning and end of austenite transformation, respectively, A1= 996 K (723°C) and A3= 1108K (835°C), solidus temperature TS= 1763K (1490°C) and the liquidus temperature TL= 1793K (1520°C). The model of joining two butt-welded flats with a thickness of 0.012 m and a width of 0.1 m is shown in Figure 5.

Figure 5.

Scheme of the welded joint with technology chamfering.

The speed and power of the movement source are assumed to be v= 0.005 m/s and P= 18000 W, respectively. Time-varying temperature field was determined according to the formulas (Eqs. (11)(15)). The temperature field on the upper surface of the welded flats and in the longitudinal section determined by the trace of the source transition is shown in Figures 6 and 7.

Figure 6.

Temperature distribution (°C) on the surface of flats at timet= 102 s from beginning of welding.

Figure 7.

Temperature distribution (°C) in longitudinal section at timet= 102 s from beginning of welding.

The quantities B1= B2= Bare considered in calculations, which are equal to the width of one flat, thereby obtaining a temperature field symmetrical to the plane defined by vertical axis of the source and the direction of its movement. Thus, both this and the further calculation illustrations are shown on the right symmetrical cross-section which is perpendicular to a moving heat source. The image of maximum temperature isotherms in a cross-section of the joint (on the front contact of flats) is shown in Figure 8.

Figure 8.

Maximum temperatures (°C) in half of cross-section of welded flats.

On the basis of the maximum temperature field achieved in particular areas of the weld joint, specific heat-affected zones were determined (Figure 9) specified by limit temperatures TL, TS, A3and A1. The weld area is limited by a liquidus isotherm, which implies covering the part of the section wherein complete melting of the material occurred. Incomplete melting is limited by both liquidus and solidus temperatures. The heat-affected zone is determined by solidus isotherms and A1. In addition, in the region between the solidus and A3temperatures, ferrite and pearlite, the starting components of steel, have been completely transformed into austenite, whereas in the temperature range A1-A3, there is an incomplete conversion. In the area where the temperature did not exceed A1, certainly, phase transitions did not occur and the structure of the parent material (ferritic-pearlitic) was thus preserved. The participation of the heating structures was determined on the basis of the Fe-C diagram. The participation of the cooling structure was determined on the basis of knowledge of the TTT-welding diagram for welding steel S235 and the variations in temperature during cooling (Figure 10) [45].

Figure 9.

Heat-affected zones (HAZ) with selected points.

The analysis of the cooling speed v8/5showed that after complete cooling in both weld and heat-affected zone, wherein there has been complete and partial transformation of ferrite and pearlite into austenite, supercooled austenite resulted in a ferrite-pearlite structure. Changes in the temperature and structure in the considered cross-section for the selected time of the welding cycle are shown in Figures 1119. The analysis of these changes was investigated on a symmetrical half of the cross-section of the flats’ connection, which is perpendicular to the direction of source movement and x0= 0.5 m from the beginning of the flats (in half their length), that is, from the beginning of the coordinate system associated with flats. For t= 0, the beginning of the coordinate systems associated with the source and the welded object overlaps. Figures 1113 illustrate the condition where the joint area was completely filled with the welded material, there were no cooling phase transitions and the largest participation of austenite is observed. By contrast, Figures 1419 illustrate the gradual decrease in the share of austenite which forms a ferritic-pearlitic structure following the temperature drop.

Figure 10.

TTT-welding diagram for S235 steel.

Figure 11.

Volume fraction of austenite in cross-section at timet= 109 s from process beginning.

Figure 12.

Volume fraction of ferrite in cross-section at timet= 109 s from process beginning.

Figure 13.

Volume fraction of pearlite in cross-section at timet= 109 s from process beginning.

Figure 14.

Volume fraction of austenite in cross-section at timet= 134 s from process beginning.

Figure 15.

Volume fraction of ferrite in cross-section at timet= 134 s from process beginning.

Figure 16.

Volume fraction of pearlite in cross-section at timet= 134 s from process beginning.

Figure 17.

Volume fraction of austenite in cross-section at timet= 146 s from process beginning.

Figure 18.

Volume fraction of ferrite in cross-section at timet= 146 s from process beginning.

Figure 19.

Volume fraction of pearlite in cross-section at timet= 146 s from process beginning.

Changes in temperature and volume fraction of the individual structural components at the selected points in the cross-section (comp. Figure 9) are shown in Figure 20.

Point 1is located in the area of chamfering of flats; hence, the graph of temperature and austenite volume fraction starts at the moment of the joint execution. At points 2(fusion zone) and 3(full-transformation zone), a complete austenite transformation occurs. But at point 2, apart from phase transformations in the solid state, melting and solidification phase transformations occur, resulting in a diverse dendritic structure shown during metallographic testing. At point 4, a partial transformation of the primary structure into austenite occurs, which is visible on the graph.

In strain calculations, linear expansion coefficients of the particular structural elements were assumed and structural stresses (Table 1) were determined on the basis of the author’s own dilatometric research [46].

Figure 20.

Welding thermal cycles and volume fractions of phases at selected points: T—temperature, A—austenite, P—pearlite, F—ferrite.

α[1/°C]γ
Austenite2.178 × 10−5γF,P,S→A1.986 × 10−3
Ferrite1.534 × 10−5γB→A1.440 × 10−3
Pearlite1.534 × 10−5γA→F,P3.055 × 10−3
Bainite1.171 × 10−5γA→B4.0 × 10−3
Martensite1.36 × 10−5

Table 1.

Structural (γ) and thermal (α) expansion coefficients of phases.

Dilatometric graphs (thermal and structural strains as a function of temperature) for selected points of the section (comp. Figure 9) are shown in Figure 21. In Figure 21a, a dilatometric graph for point 1of the weld area is presented, where the material is deposited in the liquid phase as the molten material of the electrode, flux and partially melted edges of the flats; therefore, we observe only shrinkage “negative” strains of the cooling metal with a clear fault reflecting the strain of austenite transformation into ferrite and pearlite. Point 2(Figure 21b) is located in the area of melting of the starting material, and the dashed line in the figure reflects the change of state from solid to liquid. The solid line (bottom) is initiated by weld solidification.

Figure 21.

Dilatometric curves for selected points of cross-section.

Dilatometric graph in point 3(Figure 21c) illustrates the history of strain changes for a full thermal cycle. In the considered point, the material with a ferritic-pearlitic structure is completely transformed into austenite during heating (lower line), which is shown with a fault in the diagram in the temperature range 723–835°C. At point 4 (Figure 21d) partial conversion of ferrite and pearlite in austenite was observed, that is, at that point, the temperature during heating has exceeded the temperature of the beginning of austenitizing A1, but has not reached the final temperature of austenitizing A3. The stage of the austenite transformation is determined linearly in a numerical model by making the austenite transformation progress conditionally on the ratio (TmaxA1)/(A3A1), where Tmaxis the maximum temperature reached at a particular point of the weld joint.

Advertisement

6. Verification of the results of numerical simulation of phase transitions with the results of metallographic research

In order to verify the results of the numerical calculations, metallographic tests of the butt-welded joint were carried out. For this purpose, two flats with identical geometry as calculated, that is, two flats with a thickness of 0.012 m and width of 0.1 m, were welded. The material of the flats was steel S235. Before making the joint, sheet chamfering was conducted. Then, welding was carried out with the submerged arc welding method, under welding flux Taste-3 and with SPGΦ15 wire. Welding parameters were voltage U = 30 V, current I = 600 A and welding speed 20 m/h. The diagram of the elements prepared for joining is presented in Figure 5 (identical to the welded joint adopted for numerical analysis). Cross-section of the weld joint (the image of the sample taken for metallographic examinations) is presented in Figure 22. Metallographic analysis was performed for specific zones of a welded joint, that is, for the area of weld, heat-affected zone and parent material. Figure 23 shows an image of the middle part of the welded joint (at the junction of welded flats) with a clearly visible dendritic structure which is characteristic of solidification. Figure 24 shows an image of the structure in the right symmetrical part of the welded joint in the area from the weld to the ferrite-pearlite structure of the parent material. On the border of the weld, dendrites are visible, which change in the heat-affected zone into a structure with the Widmannstatten structure elements.

Figure 22.

Cross-section of a welded joint: sample taken for metallographic analysis, magnification 2×.

Figure 23.

Weld, magnification 140×.

Figure 24.

Transition zone, magnification 140×.

Figure 25.

Weld, Nital Etch, magnification 250×, ferritic-pearlitic dendritic structure.

Figure 26.

Heat affected zone, Nital Etch, magnification 250×, ferritic-pearlitic structure.

Next, the images of structures of individual areas are presented. In the area of the weld (Figure 25), we observe the ferritic-pearlitic structure with a small amount of supercooled pearlite and island bainite in the dendritic system characteristic of solidified structures. In the heat-affected zone (Figure 26), ferrite, pearlite and pearlite balls supercooled with a number of non-metallic inclusions are visible. Primary structure of the parent material (Figure 27) consists of ferrite and pearlite in the band system (hardly visible).

Figure 27.

Parent material, Nital Etch, magnification 250×, ferritic-pearlitic banded structure.

The results of metallographic tests show high conformity with the results of numerical simulation and testify to the correctness of the developed numerical model.

Advertisement

7. Conclusion

The proposed model for determining temperature field allows to obtain characteristic zones of the weld joint with shapes and dimensions similar to the real ones. Metallographic verification of the results of numerical simulation of phase transitions also provided satisfactory results.

Analysis of phase transformations occurring during the welding process allowed determination of the quantitative composition of the cooling structures in the weld area (in the heat-affected zone).

This approach allows for an accurate tracking of changes in phase participations and the course of thermal and phase strains, which enable to determine, within the framework of thermoplasticity theory, welding stresses temporary and residual welding stresses.

References

  1. 1. Rosenthal D. Mathematical theory of heat distribution during welding and cutting. Welding J. 1941;20:220s–234s
  2. 2. Rykalin N N. Fundamentals of heat flow in welding. AN SSSR: Moskva;1947. 272 p.
  3. 3. Eagar T W, Tsai N S. Temperature fields produced by traveling distributed heat sources. Welding J. 1983;62:346s–355s
  4. 4. Goldak J, Chakravarti A, Bibby M. A new finite element model for welding heat sources. Metall Trans. 1984;15B:299–305
  5. 5. Bo K S, Cho H S. Transient temperature distribution in arc welding of finite thickness plates. Proc Inst Mech Eng. 1990;204B3:175–183
  6. 6. Radaj D. Heat effects of welding. Temperature field, residual stress, distortion, Springer-Verlag; Berlin Heidelberg. 1992; 348 p.
  7. 7. Jeong S K, Cho H S. An analytical solution to predict the transient temperature distribution in fillet arc welds. Welding J. 1997;76:223s–232s
  8. 8. Nguyen N T, Mai Y W, Simpson S, Ohta A. Analytical approximate solution for double ellipsoidal heat source in finite thick plate. Welding J. 2004;84:82s–93s
  9. 9. Kwon Y, Weckman D C. Analytical thermal model of conduction mode double sided arc welding. Sci Technol Weld Join. 2008;13:539–549
  10. 10. Fachinotti V D, Anca A A, Cardona A. Analytical solutions of the thermal field induced by moving double-ellipsoidal and double elliptical heat sources in a semi-infinite body. Int J Num Meth Biomech Eng. 2011;27:595–607
  11. 11. Antonakakis T, Maglioni C, Vlachoudis V. Closed form solutions of the heat diffusion equation with Gaussian source. Int J Heat Mass Transf. 2013;62:314–322
  12. 12. Ghosh A, Barman N, Chattopadhyay H, Hloch S. A study of thermal behaviour during submerged arc welding. Strojniški vestnik – Journal of Mechanical Engineering. 2013;59:333–338
  13. 13. Franko A, Romoli L, Musacchio A. Modelling for predicting seam geometry in laser beam welding of stainless steel. Int J Thermal Sci. 2014;79:194–205
  14. 14. Salimi S, Bahemmat P, Haghpanahi M. A 3D transient analytical solution to the temperature field during dissimilar welding processes. Int J Mechn Sci, 2014;79:66–74
  15. 15. Na S J, Lee S Y. A study on the tree- dimensional analysis of the transient temperature distribution in gas tungsten arc welding. Proc Inst Mech Eng. 1987;201B3:149–156
  16. 16. Radaj D. Finite element analysis of welding residual stress. International Conference on Residual Stresses 2, Nancy 1988, Nov. 23–25, Elsevier Applied Science, London. 1989; p. 510–516
  17. 17. Karlsson R I, Josefson B L. Three dimensional finite element analysis of temperatures and stresses in a single – pass butt-welded pipe. Trans ASME. 1990;112:76–84
  18. 18. Mundra K, DebRoy T, Kelkar K M. Numerical prediction of fluid flow and heat transfer in welding with a moving heat source. Num Heat Transf. 1996;29A:115–129
  19. 19. Komanduri R, Hou Z B. Thermal analysis of the arc welding process: part I. General solutions. Metall Mater Trans. 2000;31B:1353–1370
  20. 20. Mahapatra M M, Datta G L, Pradhan B. Three-dimensional finite element analysis to predict the effects of shielded metal arc welding process parameters on temperature distributions and weldment zones in butt and one-sided fillet welds. Proc I Mech – J Eng Manuf. 2006;220:837–845
  21. 21. Kumar A, DebRoy T. Heat transfer and fluid flow during gas-metal-arc fillet welding for various joint configurations and welding positions. Metall Mater Trans. 2007;38A:506–519
  22. 22. Wang S, Goldak J, Zhou J, Tchernov S, Downey D. Simulation on the thermal cycle of a welding process by space-time convection-diffusion finite element analysis. Int J Thermal Sci. 2009;48:936–947
  23. 23. Lindgren L E. Finite element modeling and simulation of welding. Part 1: increased complexity. J Thermal Stresses. 2001;24:141–192
  24. 24. Hongyuan F, Qingguo M, Wenli X, Shude J. New general double ellipsoid heat source model. Sci Tech Weld Join. 2005;10:361–368
  25. 25. Deng D. FEM prediction of welding residual stress and distortion in carbon steel considering phase transformation effects. Mater Dsgn. 2009;30:359–366
  26. 26. Jiang W, Liu Z, Gong J M, Tu S T. Numerical simulation to study the effect of repair width on residual stresses of stainless steel clad plate. Int J Pres Ves Pip. 2010;87:457–463
  27. 27. Lee H T, Chen C T, Wu J L. 3D numerical study of effects of temperature field on sensitisation of Alloy 690 butt welds fabricated by gas tungsten arc welding and laser beam welding. Sci Tech Weld Join. 2010;15:605–612
  28. 28. Joshi S, Hildebrand J, Aloraier A S, Rabczuk T. Characterization of material properties and heat source parameters in welding simulation of two overlapping beads on a substrate plate. Comp Mater Sci. 2013;69:559–565
  29. 29. Fu G, Lourenco M I, Duan M, Estefen S F. Effect of boundary conditions on residual stress and distortion in T-joint welds. J Constr Steel Res. 2014;02:121–135
  30. 30. Chen Y, He Y, Chen H, Zhang H, Chen S. Effect of weave frequency and amplitude on temperature field in weaving welding process. Int J Adv Manuf Techn. 2014;75:803–813
  31. 31. Tsai N S, Eagar T W. Distribution of the heat and current fluxes in gas tungsten arc. Metall Trans. 1985;16B:841–846
  32. 32. Piekarska W, Kubiak M. Three-dimensional model for numerical analysis of thermal phenomena in laser-arc hybrid welding process. Int J Heat Mass Transf. 2011;54:4966–4974
  33. 33. Chen J, Schwenk C, Wu C S, Rethmeier M. Predicting the influence angle on heat transfer and fluid flow for new gas metal arc processes. Int J Heat Mass Transf. 2012;55:102–111
  34. 34. Geissler E., Bergmann H W. Calculation of temperature profiles heating and quenching rates during laser processing. Optoelektronik Magazin. 1987;3:430–434
  35. 35. Geissler E., Bergmann H W. 3D Temperature fields in laser transformation hardening, part I: quasi-stationary fields. Optoelektronik Magazin. 1988;4:396–403
  36. 36. Carslaw H S, Jaeger J C. Conduction of heat in solids. Clarendon Press: London 1996.
  37. 37. Ruge J. Handbook of welding technology, vol. 1: Materials. Springer Verlag;Berlin:1980
  38. 38. Rhode J, Jeppson A. Literature review of heat treatment simulations with respect to phase transformation, residual stresses and distortion. Scand J Metall. 2000;29:47–62
  39. 39. Piekarska W, Kubiak M, Bokota A. Numerical simulation of thermal phenomena and phase transformations in laser-arc hybrid welded joint. Archiv Metall Materials. 2011;56:409–421
  40. 40. Avrami M. Kinetics of phase change. I. General theory. J Chem Physics. 1939;7:1103–1112
  41. 41. Parkitny R, Winczek J. Modelling of phase transformations during multipass surfacing, In: Conf. Proc. XXXVIII Sympozjon Modelling in Mechanics, Silesian University of Technology, Gliwice. 1999; p. 219–224
  42. 42. Koistinen D P, Marburger R E. A general equation prescribing the extent of the austenite-martensite transformation in pure iron-carbon alloys and plain carbon steels. Acta Metall. 1959;7:59–60
  43. 43. Domański T, Bokota A. Numerical models of hardening phenomena of tools steel base on the TTT and CCT diagrams. Archiv Metall Mater. 2011;56:325–344
  44. 44. Winczek J. A simplified method of predicting stresses in surfaced steel rods, J Mater Proces Technol, 2012;212:1080–1088
  45. 45. Brózda J, Pilarczyk J, Zeman M. TTT-welding diagrams transformation of austenite. Śląsk; Katowice. 1983; 139 p.
  46. 46. Winczek J, Kulawik A. Dilatometric and hardness analysis of C45 steel tempering with different heating-up rates. Metalurgija. 2012;51(1):9–12

Written By

Jerzy Winczek

Submitted: November 5th, 2015 Reviewed: April 28th, 2016 Published: September 21st, 2016