'-'-'-' f,'y.4;'''-Cf;
ANALYTICAL EQUATIONS INTRODUCTION A N D S Y N O P S I S
T
here is much commonaUty in the principles underlying the different methods of thermal laser material processing. The temperature field generated during vaporization cutting is analogous to that found in keyhole welding - both involve lateral heat flow from a through-thickness energy source. The phase transformations responsible for surface hardening in ferrous alloys are the same as those that determine the heat affected zone hardness of a weld - only the thermal cycle is different. Physical processes are therefore presented here as modular components of heat flow and structural change. The modules can be used interchangeably, and applied across a range of processes. This appendix contains analytical solutions for the temperature fields around a variety of energy sources and the structural transformations that they induce in materials. Analytical equations describe expUcit relationships between the process variables. Thus the effects of changes in process variables on the structure and properties of materials can be estimated quickly, or, if a particular result is desired, combinations of processing parameters can be assessed rapidly. The notation used is defined in Table E.l.
Table E.1
Notation
Symbol
Definition
Units
A A
Absorptivity Constant of precipitate solubility product Beam area Constant of precipitate solubility product Carbon equivalent Concentration of metal species of precipitate MaCb in solution Concentration of non-metal species of precipitate MaCb in solution Diffusion coefficient, D = DQ exp —(Q/RT) Pre-exponential of diffusion coefficient Beam power density Heat flux Vickers hardness number Vickers hardness number of bainite
m^ K wt% wt%
AB
B ^eq
^m
Cc D Do E F H
m
wt% m^s-i m^s-i Js-im-2 Jm-^s"^ VPN VPN (Contd)
526
Laser Processing of Engineering Materials
Table E.1
(Contd)
Symbol
Definition
Units
Hjp Hm I
Vickers hardness number of ferrite-pearlite mixture Vickers hardness number of martensite Kinetic strength of a thermal cycle Volumetric latent heat of melting Normalized latent heat of melting, L^^=Lm/[pc{Tm-To)] Volumetric latent heat of vaporization Normalized latent heat of vaporization, ^=Ly/[pc{Tm-To)] Martensite finish temperature Martensite start temperature Activation energy Gas constant, 8.314 Temperature Initial temperature Temperature at which pearlite transforms to austenite on heating Temperature at which ferrite transforms to austenite on heating Temperature at which martensite starts to form on cooling Temperature at which martensite formation is 50% complete on cooling Temperature at which martensite formation is complete on cooling Melting temperature Dissolution temperature of a precipitate Vaporization temperature Volume fraction Volume fraction of bainite Volume fraction of ferrite-pearlite mixture Volume fraction of martensite Cooling rate at 923 K Thermal diffusivity, A/(pc) Specific heat capacity Plate thickness Base of natural logarithms, 2.718 Matrix volume fraction available for precipitate dissolution Grain size Original grain size Kinetic constant Depth of treatment Precipitate size Original precipitate size
VPN VPN -
Lffi
L*
u L* Mf Ms Q R T To TACI TAC5 TMS TM50 TM/
Tm Ts Ty
V Vb ^fp Vm
v
a c d e
f g go k
I p Po
Jm"^ —
Jm"^ — K K Jmol-i Jmol-^K-i K K K K K K K K K K Kh-i m^s-i Jkg-iK-i m — m m Various m m m
(Contd)
Appendix E: Analytical equations Table E.1
527
(Contd)
Symbol
Definition
q qmax qy r r
Js~^ J s~^ J s~^ m~^ m m
t to tp V w z ZQ Zt At2/i At8/5 At^^
Beam power Peak beam power in a Gaussian distribution Volumetric rate of heat generation Radial distance from centre of a surface heat source Lateral distance from centre of a through-thickness heat source Beam radius defined where q — qmaxl^y or beam half-width Time Time for heat to diffuse over beam radius, r|/(4a) Time taken to attain peak temperature Beam traverse rate Width (spacing of two isotherms) Depth Model displacement of workpiece surface Thermal penetration depth Time to cool from T2 to Ti Time to cool from 800 to 500°C Characteristic Ats/5 for 100% martensite formation
At^
Characteristic Ats/5 for 50% martensite formation
s
TB
Units
ni s s s ms~^ m m m m s s s
At^
Characteristic Ats/5 for 0% martensite formation
s
At?
Characteristic Ats/5 ior 0% ferrite formation
s
Atp
Characteristic Ats/s for 0% pearlite formation
s
Aty^
Characteristic Ats/5 for 50% bainite formation
s
At^ A, p r
Characteristic Ats/5 for 0% bainite formation Thermal conductivity Density Beam interaction time
s Js~^ m""^ K~^ kgm~^ s
El
EQUATIONS OF HEAT FLOW
Heat flow in laser processing can be complex. However, for many processes it may be approximated to three ftindamental conditions: steady state, transient, or quasi-steady state. Fourier's first law describes steady state conditions: F = -XVT
(ELI)
where F is the heat flux (Jm~^s~^), VT is the thermal gradient (Km~^), and k is the thermal conductivity (Js~^ m~^ K~^). In this state, the temperature field does not change with time at a location in the material. Fourier's second law describes transient conditions: 3L = ^_aW^T pc dt
(E1.2)
528
Laser Processing of Engineering Materials
-^r Figure El.l Schematic illustration of limiting geometries of heatflow:(a) radial from a point source, (b) lateral from a line source; note that a single temperature profile is used for clarity - in practice the peak temperature varies as 1/r^ and 1/r for radial and lateral heatflow,respectively where q^ is the energy generation per unit time and volume (J s ^ m ^), p is the density (kg m ^),cis the specific heat capacity (Jkg~^ K~^), T is the temperature (K), a is the thermal diffiisivity (m^ s~^), and t is the time (s). In this state, a thermal cycle is experienced at a location in the material. Quasi-steady state heatflovyrdescribes a condition in which the temperature field observed from a moving energy source remains constant. The variable ^ =x — vt is defined, where ^ is the distance of the point of interest from the source, measured along the x axis (the direction of source motion) and V is velocity. If no energy is generated, equation (El.2) becomes d^T d^T d^T _j_ _l_ a|2 ^ 3^2 ^ 9^2
—
vdT ad^
(E1.3)
In a state of quasi-steady state heat flow, transient effects can be ignored, which enables analytical expressions for temperature fields to be derived. In addition, heat flow in most laser-based processes may be considered to be limited to one of two types of source: a surface source, from which heat is conducted radially into the material (three-dimensional or *thick plate' heat flow); or a through-thickness source, from which heat is conducted laterally into the material (two-dimensional or *thin plate' heat flow). Radial and lateral heat flow are illustrated schematically for a structural steel in Figs El. 1(a) and (b), respectively. Analytical solutions to equations (E1.1)-(E1.3) can be derived by making assumptions, such as: the material is homogeneous and isotropic; heat flow occurs exclusively by conduction; and no energy is generated by material transformations. Under conditions typical of thermally induced change during laser processing, these are justifiable for the level of accuracy sought in this treatment (±5%).
u
UJ
(£
D 0 (0 >-
tr
z u 0
kJ < IL
o go a
u
I
o
S
O
n3
0 U (0 ,^ Q
Z D O < Q J ui UJ
D
H < Q.
UJ
f2
CM i
u
3 O CO
o
o pq
pq ^3
PH
in
ON in
a^
ON ON
v^
1—H
OJ
^>bJD ^ '^
C
^
^
^-^ j:) CO 00 ON
en 00 OS
CO 00 ON i-H
03
^ cd C/3
OJ
u
OJ
u
CQ
C/5
c«
C/3
— 'c^' PQ
C/5
»—' c/T cd
i-i
^ ^ 1Z CQ
n
U
oJ
'^ 00 00
^ ^ C/5
a
1 X
ON
I^
03
;3
a
c
^
•i-H
VH
S ^2
£? << i=i
U
a
\D
I
M
in
I
^ v:^ ^ '^ 1-
+
^
oq
C
'-< (:^
CO
<:?r^
)-H
s
OJ
3
^ 3 CO
cyo
'M
o3
'^M
fl ^O
03
u
^
£r OJ
u^
OJ
a ,o
b c ,o
OJ +-»
pc^
I
(N
(NOQ
:§
o
I
1-H
(N
;§
oJ
O
OJ +-* C/5
'-M
;N
+ +
H
(N
a
c o
o ^-^ ^ n M
Q>
o
8 CO
o
.2
;H
Is
O
C ;3
•4-»
M
C/5
1H Hl:
ON
ON
<
00
OS W)
1^ g-
q T^
"1 .S *§ ^
^ ^
&
o
I
<
^ ^
1^
t4^
4 ^
00
as
bJO
o
^ ^
00
C3
C/3
+
«1
3
I
K+
I
H5>
I ^
o
+ +
I
^ ^
II
+
+
+
532
Laser Processing of Engineering Materials
E3
TEMPERATURE FIELD A R O U N D A THROUGH-THICKNESS ENERGY SOURCE
Table E3.1
Analytical equations for the temperature field around a through-thickness energy source
Type
Shape
Equation
Eqn no.
Source
Stationary
Line
Aq 27tk (Ti - T2) ~d ''
E3.1
(Klemens, 1976)
E3.2
(Rosenthal, 1941; Ashby and Easterling, 1982)
•"(^) Moving
Line
Aq T(r,t) — To =
1/2
1
dT _ dt
(vd^ \Aq)
Tp-To
^Aq( \7Te) vd \2neJ
r^
1 .r A{7tkpcy/^^
T^3/2 ^^
E3.3
1/2
E3.4
pa per 2
E3.5
E3.6
47rXpce {Tp - ro)2 Atl9| e(rp - ToV
T-To
^2
la
E3.7
exp
1 1 where — = • ^""^^^el {T,-To)'
1 {T2-T0)'
and T2 > Ti 1/2
E3.8
Ine) vd \27t 1 pc {Tp,-To)
(Tp,-To)
Appendix E: Analytical equations
E4
533
KINETIC EFFECT OF A THERMAL CYCLE A N D DIFFUSION-CONTROLLED STRUCTURAL CHANGES
Table E4.1 Analytical equations for the kinetic effect of a thermal cycle and diffusion-controlled structural changes Equation Kinetic effect
I = atptxp-
—
Eqn no.
Source
E4.1
(Ion eta/., 1984)
E4.2
(Avrami, 1939; Johnson and Mehl, 1939)
E4.3
(Lifshitz and Slyozov, 1961; Wagner, 1961)
E4.4
(Ion eta/., 1984)
where Aq 1 1 Q V lirXe Tp — TQ for radial heat flow, and iTtRTty
enidtt)
a =
a=2
[WRJ^
/AqV
1 AnXpce {Tp — TQY
for lateral heat flow 9
9
Q
Grain growth
g'-g^
Precipitate coarsening
p' - PQ = -y- ^^P " ^
Precipitate dissolution
Ts =
= katpexp-
—
^ P
B A-log
/'
where I 3/2
/ = 1 — exp
oitp
_ Q
1 Tp
1 T;
and starred variables denote calibration conditions
534
E5
Laser Processing of Engineering Materials
P H A S E T R A N S F O R M A T I O N S A N D H A R D N E S S IN CARBON-MANGANESE STEELS
Because of their importance in engineering, analytical equations describing phase transformations and hardness in plain carbon-manganese steels are provided in Table E5.1. Figure E5.1 illustrates the critical cooling times for phase formation. Note that the carbon equivalent Ceq used here to calculate critical cooling times is that defined in equation (E5.7). Table E5.1 in wt%
Empirical and theoretical equations for phase transformations in steels: alloying additions
Phase transformation temperatures
Carbon equivalent Critical cooling times
Phase volume fractions Phase hardness
Average hardness
Equation
Eqn no.
Source
= 996 - 30Ni - 25Mn - 5Co + 25Si + 30A1 + 25MO + 50V TAC3(K) = 1183 --416C + 228C2 -40Cr - 50Mn - 40Ni + 800P + 60V + 130Mo + 50W + 50S T^(K) = 1810-90C TMSCK) = 812 - 423C - 30.4Mn -17.7Ni-12.lCr-7.5Mo
E5.1
(Andrews, 1965)
E5.2
(Kumar, 1968)
E5.3 E5.4
(Hansen, 1958) (Andrews, 1965)
TM50 (K) = TMS — 47
E5.5
TMf (K) = T M S - 2 1 5
E5.6
Q^ = C + Mn/12 +Si/24
E5.7
At^ = exp (17.724Ceg - 2.926) At^ = exp {l9.954Ceq - 3.944) At^ = exp(16.929Ce^ + 1.453) At^o = exp(ln(At^o • Atp/2) y^=exp{ln(0.5).(At/A40)2} n = exp{ln(0.5) • (At/At'^^f) - Vm
E5.8 E5.9
(Steven and Haynes, 1956) (Steven and Haynes, 1956) (Inagaki and Sekiguchi, 1960) (Ion eta/., 1984) (Ion eta/., 1984)
TACI (K)
Vfp = l-{Vm + Vb) Hm = 295-{-515 Ceq Hb = 223-\-U7 Ceq Hjp = UO-\-139 Ceq Hmax = VmHm + VyHy + VfpHjp
E5.10 E5.11 E5.12 E5.13 E5.14 E5.15 E5.16 E5.17 E5.18
(Ion eta/., (Ion eta/., (Ion eta/., (Ion eta/.,
1984) 1984) 1984) 1984)
(Ion eta/., (Ion eta/., (Ion eta/., (Ion eta/.,
1984) 1984) 1984) 1984)
Appendix E: Analytical equations
Cooling time 800-500°C, At (s)
At? 1
•w—•*-~....^—w
—
*^ ^ \ ^N
Bainite V ^
\
\\ \ \ \ \ ^ ^ cu
535
Martensite
\
\
\ \
^'"^''
\\\ \ \ \ \
\ \
\
^
\ v
\
\.
\ Pearlite
\
\ Hypoeutectoid steel
\
-^^=#
AfO
^^^^
^^
Hardness
-^^«
Afg
Applied energy (J mm~2)
Figure E5.1 Critical cooling times for the formation of various microstructures in hypoeutectoid carbonmanganese steels, with a schematic hardness variation superimposed (broken line)
BIBLIOGRAPHY Andrews, K.W. (1965). Empirical formulae for the calculation of some transformation temperatures Journal of the Iron and Steel Institute, 203, 721-727. Ashby, M.F. and Easterling, K.E. (1982). A first report on diagrams for grain growth in welds. Acta Metallurgical 30, (11), 1969-1978. Ashby M.R and Easterling, K.E. (1984). The transformation hardening of steel surfaces by laser beams - 1. Hypoeutectoid steels. Acta Metallurgical 32,(11), 1935-1948. Avrami, M. (1939). Kinetics of phase change I. Journal of Chemical Physics, 7,1103-1112. Bass, M. (1983a). Laser heating of solids. In: Bertolotti, M. ed. Physical Processes in Laser-material Interactions. New York: Plenum Press, p. 83. Bass, M. (1983b). Laser heating of solids. In: Bertolotti, M. ed. Physical Processes in Laser-material Interactions. New York: Plenum Press, p. 86. Bass, M. (1983c). Laser heating of solids. In: Bertolotti, M. ed. Physical Processes in Laser-material Interactions. New York: Plenum Press, p. 87. Carslaw, H.S. and Jaeger, J.C. (1959). Conduction of Heat in Solids. 2nd ed. Oxford: Clarendon Press.
Hansen, M. (1958). Constitution of Binary Alloys. New York: McGraw-HiU. Inagaki, M. and Sekiguchi, H. (1960). Continuous cooling transformation diagrams of steels for welding and their applications. Transactions of the National Institute for Metals, Japan, 2, 102-125. Ion, J.C, Easterling, K.E. and Ashby M.R (1984). A second report on diagrams of microstructure and hardness for heat-affected zones in welds. Acta Metallurgical 32, (11), 1949-1962. Johnson, W.A. and Mehl, R.R (1939). Reaction kinetics in processes of nucleation and growth. Transactions of the American Institute of Mining and Metallurgical Engineering, Iron Steel Division, 135, 416-458. Klemens, RG. (1976). Heat balance and flow conditions for electron beam and laser welding. Journal ofApplied Physics, 47, (5), 2165-2X7 A. Kumar, R. (1968). Physical Metallurgy of Iron and Steel. London: Asia. Lifshitz, I.M. and Slyozov, V.V. (1961). The kinetics of precipitation from supersaturated solid solutions. Journal of the Physics and Chemistry of Solids, 19, 35-50. Ready J.R (1971). Effects of High-Power Laser Radiation. New York: Academic Press.
536
Laser Processing of Engineering Materials
Rosenthal, D. (1941). Mathematical theory of heat Steven, W. and Haynes, A.G. (1956). The temperature distribution during welding and cutting. Welding of formation of martensite and bainite in low-alloy steels. Journal of the Iron and Steel Institute, 183, Journal 20,220s-234s. 349-359. Rosenthal, D. (1946). The theory of moving sources of heat and its application to metal treatments. Wagner, C. (1961). Theory of precipitate ageing Transactions of the American Society of Mechanical during re-solution (Ostwald ripening). Zeitschrift Engineers, 68,849-866. der ElektrochemiCy 65, 581-591. (In German.) Rykalin, N., Uglov, A. and Kokora, A. (1978). Laser Machining and Welding. Oxford: Pergamon Press.