1 Heat Flow and Temperature Distribution in Welding
1 Heat Flow and Temperature Distribution in Welding
1 Heat Flow and Temperature Distribution in Welding
1.1 Introduction
The symbols and units used throughout this chapter are defined in Appendix 1.1.
Since heat losses from free surfaces by radiation and convection are usually negligible in
welding, the temperature distribution can generally be obtained from the fundamental differ-
ential equations for heat conduction in solids. For uniaxial heat conduction, the governing
equation can be written as:1
(i-D
where T is the temperature, t is the time, x is the heat flow direction, and a is the thermal
diffusivity. The thermal diffusivity is related to the thermal conductivity X and the volume
heat capacity pc through the following equation:
(1-2)
d-3)
and
(1-4)
The above equations must clearly be satisfied by all solutions of heat conduction problems,
but for a given set of initial and boundary conditions there will be one and only one solution.
A pre-condition for obtaining simple analytical solutions to the differential heat flow equa-
tions is that the thermal properties of the base material are constant and independent of tem-
perature. For most metals and alloys this is a rather unrealistic assumption, since both X, a,
and pc may vary significantly with temperature as illustrated in Fig. 1.1. In addition, the
thermal properties are also dependent upon the chemical composition and the thermal history
of the base material (see Fig. 1.2), which further complicates the situation.
By neglecting such effects in the heat flow models, we impose several limitations on the
application of the analytical solutions. Nevertheless, experience has shown that these prob-
lems to some extent can be overcome by the choice of reasonable average values for X, a and
pc within a specific temperature range. Table 1.1 contains a summary of relevant thermal
properties for different metals and alloys, based on a critical review of literature data. It should
be noted that the thermal data in Table 1.1 do not include a correction for heat consumed in
melting of the parent materials. Although the latent heat of melting is temporarily removed
during fusion welding, experience has shown this effect can be accounted for by calibrating
the equations against a known isotherm (e.g. the fusion boundary). In practice, such correc-
tions are done by adjusting the arc efficiency factor Tq until a good correlation is achieved
between theory and experiments.
Carbon steel
Hx-H0 = PC(T-T0 ),J/mm3
Temperature, 0C
Fig. 1.1. Enthalpy increment H7-H0 referred to an initial temperature T0 = 200C. Data from Refs.
2-4.
Table 1.1 Physical properties for some metals and alloys. Data from Refs 2 - 6 .
X9 W/mm 0C
Temperature, 0C
(b)
High alloy steel
X, W/mm 0C
Temperature, 0C
Fig. 1.2. Factors affecting the thermal conductivity X of steels; (a) Temperature level and chemical
composition, (b) Heat treatment procedure. Data from Refs. 2-4.
The concept of instantaneous heat sources is widely used in the theory of heat conduction.1 It
is seen from Fig. 1.3 that these solutions are based on the assumption that the heat is released
instantaneously at time t - 0 in an infinite medium of initial temperature T0, either across a
plane (uniaxial conduction), along a line (biaxial conduction), or in a point (triaxial conduc-
tion). The material outside the heat source is assumed to extend to x = + °° for a plane source
in a long rod, to r = °° for a line source in a wide plate, or to R = °° for a point source in a heavy
slab. The initial and boundary conditions can be summarised as follows:
T-T0 = oo for t = O and x = O (alternatively r = O or R = O)
T-J 0 = O for t = O and x * O (alternatively r > O or 7? > O)
7-T 0 = O for O < t < oo when x = ± oo (alternatively r = oo or R = oo).
It is easy to verify that the following solutions satisfy both the basic differential heat flow
equations (1-1), (1-3) and (1-4) and the initial and boundary conditions listed above:
where Q is the net heat input (energy) released at time t = O, and A is the cross section of the
rod.
(1-6)
where d is the plate thickness.
Equations (1-5), (1-6) and (1-7) provide the required basis for a comprehensive theoretical
treatment of heat flow phenomena in welding. These solutions can either be applied directly
or be used in an integral or differential form. In the next sections a few examples will be given
to illustrate the direct application of the instantaneous heat source concept to problems related
to welding.
(a) T
Fig. 1.3. Schematic representation of instantaneous heat source models; (a) Plane source in a long rod.
(b) T
X y
(C)
Fig. 1.3.Schematic representation of instantaneous heat source models (continued); (b) Line source in a
wide plate, (c) Point source in a heavy slab.
1.5 Local Fusion in Arc Strikes
The series of fused metal spots formed on arc ignition make a good case for application of
equation (1-7).
Model
The model considers a point source on a heavy slab as illustrated in Fig. 1.4. The heat is
assumed to be released instantaneously at time t = 0 on the surface of the slab. This causes a
temperature rise in the material which is exactly twice as large as that calculated from equation
(1-7):
(1-8)
(1-9)
— Dimensionless time:
d-10)
(1-11)
where qo is the net arc power (equal to Qlt(), and (Hc-Ho) is the heat content per unit volume at
the reference temperature.
(1-12)
(1-13)
Heat source
Isotherms
Fig. 1.4. Instantaneous point source model for assessment of temperatures in arc strikes.
T1
^i
Equation (1-13) has been solved numerically for different values ofCT1and T1. The results
are presented graphically in Fig. 1.5. Due to the inherent assumption of instantaneous release
of heat in a point, it is not possible to use equation (1-13) down to very small values OfCT1 and
T1. However, at some distance from the heat source and after a time not much shorter than the
real (assumed) time of heating, the calculated temperature-time pattern will be reasonably
correct. Note that the heavy broken line in Fig. 1.5 represents the locus of the peak tempera-
tures. This locus is obtained by setting 3In(OAi1VdT1 = 0:
from which
(1-14)
where Qp is the peak temperature, and e is the natural logarithm base number.
Example (1.1)
Consider a small weld crater formed in an arc strike on a thick plate of low alloy steel. Calcu-
late the cooling time from 800 to 5000C (Af875), and the total width of the fully transformed
region adjacent to the fusion boundary. The operational conditions are as follows:
where r| is the arc efficiency factor. Relevant thermal data for low alloy steel are given in
Table 1.1.
Solution
In the present case it is convenient to use the melting point of the steel as a reference tempera-
ture (i.e. 0 = 0m = 1 when Tc = TJ. The corresponding values OfZi1 and 9 (at 800 and 5000C,
respectively) are:
and
Alternatively, the same information could have been read from Fig. 1.5. Although it is
difficult to check the accuracy of these predictions, the calculated values for Ats/5 and ARlm are
considered reasonably correct. Thus, because the cooling rate is very large, in arc strikes a
hard martensitic microstructure would be expected to form within the transformed parts of the
HAZ, in agreement with general experience.
Equation (1-6) can be used for an assessment of the temperature-time pattern in spot welding
of plates.
Model
The model considers a line source which penetrates two overlapping plates of similar thermal
properties, as illustrated in Fig. 1.7. The heat is assumed to be released instantaneously at time
Heat source
Heat source d
Fig. 1.7. Idealised heat flow model for spot welding of plates.
t = 0. If transfer of heat into the electrodes is neglected, the temperature distribution is given
by equation (1-6).
This equation can be written in a dimensionless form by introducing the following group of
parameters:
— Dimensionless time:
(1-15)
(1-16)
(1-17)
(1-18)
e/n2
6/n2
T
2
T2
Figure 1.8 shows a graphical representation of equation (1-18) for a limited range of a 2 and
T2. A closer inspection of the graph reveals that the temperature-time pattern in spot welding
is similar to that observed during arc ignition (see Fig. 1.5). The locus of the peak tempera-
tures in Fig. 1.8 is obtained by setting d\n{^ln7}ldx2 - 0.
which gives
and
(1-19)
Example (1.2)
Consider spot welding of 2 mm plates of low alloy steel under the following operational con-
ditions:
Calculate the cooling time from 800 to 5000C (Af8/5) in the centre of the weld, and the cooling
rate (CR.) at the onset of the austenite to ferrite transformation. Assume in these calculations
that the total voltage drop between the electrodes is 1.6 V. The M^-temperature of the steel is
taken equal to 475°C.
Solution
If we use the melting point of the steel as a reference temperature, the parameters n2 and 6 (at
800 and 5000C, respectively) become:
and
and
Since the cooling curves in Fig. 1.8 are virtually parallel at temperatures below 8000C (i.e.
for QZn2 < 0.15), the computed values of Ar8/5 and CR. are also valid for positions outside the
weld centre-line. In the present example the centre-line solutions can be applied down to
(°"2m)2 ~ 2. According to equation (1-19), this corresponds to a lower peak temperature of:
Thermit welding is a process that uses heat from exothermic chemical reactions to produce
coalescence between metals and alloys. The thermit mixture consists of two components, i.e.
a metal oxide and a strong reducing agent. The excess heat of formation of the reaction prod-
uct provides the energy source required to form the weld.
Model
In thermit welding the time interval between the ignition of the powder mixture and the com-
pletion of the reduction process will be short because of the high reaction rates involved.
Assume that a groove of width 2L1 is filled instantaneously at time t = 0 by liquid metal of an
initial temperature Tt (see Fig. 1.10). The metal temperature outside the fusion zone is T0. If
heat losses to the surroundings are neglected, the problem can be treated as uniaxial conduc-
tion where the heat source (extending from -L 1 to +L1) is represented by a series of elementary
sources, each with a heat content of:
(1-20)
At time t this source produces a small rise of temperature at position JC, given by equation (1 -5):
(1-21)
(1-22)
Isl'srau*'*'=]
Fusion zone
Fig. 1.9. Calculated peak temperature contours in spot welding of steel plates (numerical solution). Op-
erational conditions: / = 23kA, 64 cycles. Data from Bently et al1
Fusion zone
Fig. 1.10. Idealised heat flow model for thermit welding of rails.
where erf(u) is the Gaussian error function. The error function is defined in Appendix 1.3*.
Because of the complex nature of equation (1-22), it is convenient to present the different
solutions in a dimensionless form by introducing the following groups of parameters:
*The error function is available in tables. However, in numerical calculations it is more convenient to use the
Fortran subroutine given in Appendix 1.3.
Dimensionless temperature:
(1-23)
Dimensionless time:
(1-24)
Dimensionless jc-coordinate:
(1-25)
(1-26)
Equation (1-26) has been solved numerically for different values of Q and T3. The results
are presented graphically in Fig. 1.11. As would be expected, the fusion zone itself (Q < 1)
cools in a monotonic manner, while the temperature in positions outside the fusion boundary
(Q > 1) will pass through a maximum before cooling. The locus of the HAZ peak temperatures
in Fig. 1.11 is defined by 3673T3 = 0. Referring to Appendix 1.3, we may write:
which gives
(1-27)
The peak temperature distribution is obtained by solving equation (1-27) for different com-
binations of Qm and T3m and inserting the roots into equation (1-26).
Example (1.3)
Consider thermit welding of steel rails (i.e. reduction of Fe2O3 with Al powder) under the
following operational conditions:
Calculate the cooling time from 800 to 5000C in the centre of the weld, and the total width
of the fully transformed region adjacent to the fusion boundary. The Ac3-temperature of the
steel is taken equal to 8900C.
Definition of parameters:
91
T
3
Solution
For positions along the weld centre-line (Q. = 0) equation (1-26) reduces to:
and
The locus of the 8900C isotherm in temperature-time space can be read from Fig. 1.11.
Taking the ordinate equal to 0.40, we get:
By inserting this value into equation (1-27), we obtain the corresponding coordinate of the
isotherm:
Unfortunately, measurements are not available to check the accuracy of these predictions.
Systematic errors would be expected, however, because of the assumption of instantaneous
release of heat immediately after powder ignition and the neglect of heat losses to the sur-
roundings. Nevertheless, the present example is a good illustration of the versatility of the
concept of instantaneous heat sources, since these solutions can easily be added in space as
shown here or in time for continuous heat sources (to be discussed below).
Friction welding is a solid state joining process that produces a weld under the compressive
force contact of one rotating and one stationary workpiece. The heat is generated at the weld
interface because of the continuous rubbing of the contact surfaces, which, in turn, causes a
temperature rise and subsequent softening of the material. Eventually, the material at the
interface starts to flow plastically and forms an up-set collar. When a certain amount of up-
setting has occurred, the rotation is stopped and the compressive force is maintained or slightly
increased to consolidate the weld.
(1-28)
(a)
Continuous heat source
(b) q
t
Fig. 1.12. Idealised heatflowmodel for friction welding of rods; (a) Sketch of model, (b) Subdivision of
time into a series of infinitesimal elements dt'.
At time / this heat will cause a small rise of temperature in the material, in correspondance
with equation (1-5):
(1-29)
If we substitute t"=t-1'into equation (1 -29), the total temperature rise at time t is obtained
by integrating from t"= t (t'= 0) to /"= 0 (t'= t):
(1-30)
In order to evaluate this integral, we will make use of the following mathematical transfor-
mation:
where
and
The latter integral can be expressed in terms of the complementary error function* erfc{u)
by substituting:
(1-31)
If the temperature of the contact section at the end of the heating period is taken equal to Th,
equation (1-31) can be rewritten as:
(1-32)
where t'h denotes the duration of the heating period (t < t'h). Measured contact section tem-
peratures for different metal/alloy combinations are given in Table 1.2.
Equation (1-32) may be presented in a dimensionless form by the use of the following
groups of parameters:
Dimensionless temperature:
(1-33)
Dimensionless time:
(1-34)
Based on direct readings of the voltage drop between the two work-pieces.
— Dimensionless .^-coordinate:
(1-35)
(1-36)
Equation (1-36) describes the temperature in different positions from the weld contact sec-
tion during the heating period. However, when the rotation stops, the weld will be subjected to
free cooling, since there is no generation of heat at the interface. As shown in Fig. 1.13(a) this
can be accounted for by introducing an imaginary heat source of power +qo at time t = t'h
which acts simultaneously with an imaginary heat sink of negative power -q o. It follows from
the principles of superposition (see Fig. 1.13b) that the temperature during the cooling period
is given by:9
(1-37)
where 6"(x4) and 6"(T 4 - 1) are the temperatures calculated for the heat source and the heat
sink, respectively, using equation (1-36).
Equations (1-36) and (1-37) have been solved numerically for different values of Q'and T4.
The results are presented graphically in Fig. 1.14. Considering the contact section (Q'= 0), the
temperature increases monotonically with time during the heating period, in correspondance
with the relationship:
(1-38)
q
(a)
e"
(b)
Heating
period
$ffl9
Fig. 1.13. Method for calculation of transient temperatures during friction welding; (a) Sketch of
imaginary heat source/heat sink model, (b) Principles of superposition.
(1-39)
Outside the contact section (Q / > 0), the temperature rise will be smaller and the cooling
rate lower than that calculated from equations (1-38) and (1-39).
Heating Cooling
e"
Example (1.4)
Consider friction welding of 026mm aluminium rods (Al-Cu-2Mg) under the following con-
ditions:
Calculate the peak temperature distribution across the joint. Assume in these calculations
that the thermal diffusivity of the Al-Cu-2Mg alloy is 70mm2 s"1.
Solution
Readings from Fig. 1.14 give:
Next Page
In this particular case, it is possible to check the accuracy of the calculations against in situ
thermocouple measurements carried out on friction welded components made under similar
conditions. A comparison with the data in Fig. 1.15 shows that the model is quite successful in
predicting the HAZ peak temperature distribution. In contrast, the weld heating and cooling
cycles cannot be reproduced with the same degree of precision. This has to do with the fact
that the present analytical solution omits a consideration of the plastic straining occurring
during friction welding, which displaces the coordinates and alters the heat balance for the
system.
In most fusion welding processes the heat source does not remain stationary. In the following
we shall assume that the source moves at a constant speed along a straight line, and that the net
power supply from the source is constant. Experience shows that such conditions lead to a
fused zone of constant width. This is easily verified by moving a tungsten arc across a sheet of
steel or aluminium, or by moving a soldering iron across a piece of lead or tin. Moreover,
zones of temperatures below the melting point also remain at constant width, as indicated by
the pattern of temper colours developed on welding ground or polished sheet.
It follows from the definition of pseudo-steady state that the temperature will not vary with
time when observed from a point located in the heat source. Under such conditions the tem-
perature field around the source can be described as a temperature 'mountain' moving in the
direction of welding (e.g. see Fig. 15 in Ref. 11). For points along the weld centre-line, the
temperature at different positions away from the heat source (which for a constant welding
speed becomes a time axis) may be presented in a two-dimensional plot as indicated in Fig.
1.16. Specifically, this figure shows a schematic representation of the temperature in steel
welding from the base plate ahead of the arc to well into the solidified weld metal trailing the
arc. If we consider a fixed point on the weld centre-line, the temperature will increase very
rapidly during the initial period, reaching a maximum of about 2000-22000C for positions
immediately beneath the root of the arc.11 When the arc has passed, the temperature will start
to fall, and eventually (after long times) approach that of the base plate. In contrast, an ob-
server moving along with the heat source will always see the same temperature landscape,
since this will not change with time according to the presuppositions.
It will be shown below that the assumption of pseudo-steady state largely simplifies the
mathematical treatment of heat flow during fusion welding, although it imposes certain re-
strictions on the options of the models.