Transformations in TRIP-assisted Steels: Microstructure and Properties
Transformations in TRIP-assisted Steels: Microstructure and Properties
Transformations in TRIP-assisted Steels: Microstructure and Properties
0
limit (considering the stored energy of
bainite), as illustrated in Figure 2.6.
Some elements like silicon, aluminium and phosphorous retard ce-
mentite precipitation allowing the carbon to remain dissolved in austenite.
The carbon-enriched residual austenite becomes stable against any further
martensitic transformation during subsequent cooling [Andrews, 1965], and
hence is retained. The pronounced eect of carbon in stabilising austenite is
therefore exploited in retaining austenite, thus permitting the concentration
of other solutes to be kept to a minimum.
12 2.3 Microstructural evolution
Figure 2.5: Schematic illustration of bainite reaction mechanism [Bhadeshia,
2001].
Microstructural evolution during the isothermal formation of bainite is
identical for both hot and cold rolled TRIP-assisted steels. A typical trans-
formation map is shown in Figure 2.7 [Jacques et al., 2001b]. As the bainite
reaction progresses, the residual austenite becomes enriched with carbon and
gains stability against martensitic transformation during cooling. When the
enrichment is inadequate at low holding times, the austenite may decompose
partly into martensite during cooling. On the other hand, amount of un-
13 2.3 Microstructural evolution
Figure 2.6: The incomplete reaction phenomenon: bainite reaction stops as
carbon in austenite approaches T
0
limit, x is the bulk carbon content of
the steel and A
e
3
= 1 exp[{1 }
n
] (2.1)
where f
= 1 exp(k
Z
) (2.2)
32 2.7 Strain-induced martensite formation
with empirical constants k = 28 and Z = 3.7.
A dierent parabolic type of dependence was proposed by Gerberich et
al. [1970]:
f
= A
2
1/2
(2.3)
where A
2
is an empirical constant.
As can be seen, none of these models include volume fraction of austen-
ite and therefore exclude any possible interaction between the parent and the
transformed phases. Angel [1954] and Ludwigson and Berger [1969] consid-
ered this and proposed the following model:
f
= A
1
B
V
(2.4)
where A
1
and B are constants, with B = 3 for stainless steels.
All these models are formulated on the basis of the experimental data
for fully austenitic steels. The strain-induced transformation of retained
austenite in TRIP-assisted steel has been expressed by Sugimoto et al.
[1992b] as
ln(V
0
) ln(V
) = k (2.5)
where V
0
is the
volume fraction of the untransformed austenite at a plastic strain . The
empirical constant k must be measured for each steel.
A novel way of describing strain-induced transformation of retained
austenite in these TRIP-assisted steels has been proposed by Sherif et al.
[2004], a development of the model of Sugimoto et al. [1992b]. The frac-
tion of untransformed austenite is expressed as a function of plastic strain,
33 2.8 Special properties
deformation temperature, chemical composition and the original amount of
austenite. The nal equation is given as,
ln(V
0
) ln(V
) = k
1
G
(2.6)
where k
1
takes a value of 0.002017 J mol
1
, independent of the steel compo-
sition.
All these models explain adequately the kinetics of strain-induced for-
mation of martensite. But none takes into account the possibility of mechan-
ical stabilisation at large plastic strains.
2.8 Special properties
2.8.1 Formability
Formability is dened as the ease with which an object having a complex
geometry can be manufactured using the material. This can approximately
be related to plasticity or ductility which tends to deteriorate with increasing
strength. TRIP-assisted steels are unique in this respect, they exhibit better
formability than several other steels of comparable strength [Konieczny, 2003;
Sugimoto et al., 1996].
Formability is not assessed by a single attribute but expressed as a com-
bination of properties such as deep drawability, stretch formability, stretch
angeability and bendability. The state of stress is dierent in each of these
forming operations, as illustrated schematically in Figure 2.18.
34 2.8 Special properties
Figure 2.18: State of stress existing during dierent modes of forming oper-
ations [Takahashi, 2003].
Deep drawability
This is expressed by the Lankford parameter or r value, the ratio of true
strain in the width direction to that in thickness direction. r is the average
of r values measured on specimens obtained at dierent directions, namely at
0
, 45
and 90
C, at which
iron and zinc share great anity and allow an alloy to form, whereas pure
zinc prevails at the surface. The nal product is a steel surface protected
with a zinc coating.
The cold rolled sheets are initially heated to a temperature in the (+)
phase region. This is followed by another dwelling at a lower temperature
of some 400
C to 75
C to 300
C for 48 h. Cylindrical
47 3.2 Experiments
samples of length 12 mm and diameter 8 mm were made out of this steel
and sealed in quartz tubes containing a partial pressure of argon for heat
treatment in an air furnace.
The sealed samples were austenitised at 1200
C, 1100
C, 1000
C and
950
C.
Figure 3.2: Micrograph of the sample quenched after austenitising at 1100
C.
49 3.3 Results and discussion
Figure 3.3: Micrograph of the sample quenched after austenitising at 1000
C.
Figure 3.4: Micrograph of the sample quenched after austenitising at 950
C.
Martensite plates exist with a distribution of plate lengths in each of the
samples. Growth of the initial plates is restricted only by the austenite grain
boundaries, while those formed subsequently are also stopped by the pre-
existing plates. New plates forming during the progress of transformation
therefore tend to be ner than the earlier ones, the coarsest ones in each
sample being given by the prior austenite grain size. The microstructures of
50 3.3 Results and discussion
these samples therefore get rened at lower austenitising temperatures that
lead to ner prior austenite grain sizes (Table 3.1).
Austenitising temperature Mean lineal intercept % Error
(
C) (m)
1200 410 13
1100 124 20
1000 65 18
950 25 26
Table 3.1: Variation of austenite grain size with austenitising temperature,
error represents 1.
Quenched samples in the polished but unetched conditions were exam-
ined using an optical microscope to observe any cracks. Samples austenitised
at 1200
C and 1100
C and 1000
C for
another 10 minutes before nally quenching in water. Similar microcracks as
the sample directly quenched from 1200
C: (a) etched
sample; (b) periodically cracked plate in unetched sample.
The apparent distances between two adjacent cracks in individual frag-
mented plates of martensite in the samples austenitised at 1200
C and 1100
C
were measured using optical metallography. These measurements follow a
distribution as depicted in Figure 3.8. The distances are designated l
a
to
52 3.3 Results and discussion
(a)
(b)
Figure 3.6: Microcracking in the sample austenitised at 1000
C: (a) etched
sample where cracks are dicult to detect; (b) occasional cracks visible in
unetched sample.
emphasise that they are apparent measures because stereological eects are
not taken into account. The mean values for these measurements were found
to be l
a
= 58 m and l
a
= 38 m for the samples austenitised at 1200
C and
53 3.3 Results and discussion
Figure 3.7: Microcracking in the sample austenitised at 1200
C and soaked
at 950
C before quenching.
1100
C respectively.
The minimum value observed in all the measurements, within experi-
mental error, was found to be 10 m. This is more relevant than the mean
value in the present context. If it is assumed that the critical stress transfer
length l
C
= 10 m, then an austenite grain size which is somewhat less than
this value would avoid quench cracking of the type explored in this work.
One diculty is that cracks were not detected in the T
=950
C sample with
a prior austenite grain size of 25 m, which is greater than l
C
. It is possible
that the cracks exist but have been missed because the sample contains few
plates of length greater than l
C
. Therefore, the plate length distributions
were also measured for each sample.
The apparent lengths of approximately 180 martensite plates from each
sample were measured using optical metallography. In all cases the plates
were chosen at random by arbitrarily translating the stage across a cursor
which identied the plate to be measured. The size distributions are illus-
54 3.3 Results and discussion
Figure 3.8: The frequency of observation versus the apparent distance be-
tween cracks on individual martensite plates in the samples with the coarsest
microstructures (T
C).
trated in Figure 3.9, plotted on identical scales to facilitate comparisons.
It is evident that a coarse grain size tends to have a more uniform dis-
tribution of plate sizes (Figure 3.9(a)) whereas the smallest plates dominate
when the austenite grain size is reduced (Figure 3.9(d)). This is consistent
with the work by Guimaraes and co-workers on the heterogeneity of marten-
sitic transformation when the austenite grain size is small [Guimaraes and
55 3.4 Simulating microstructural evolution
Saavedra, 1984; Guimaraes and Gomes, 1978]. Some austenite grains trans-
form via a burst of martensite, while others remain fully austenitic. Such
bursts can be expected to lead to a greater preponderance of ne plates. By
contrast, coarse grained microstructures transform uniformly with the ma-
jority of austenite grains participating, thus allowing a gradual evolution of
microstructure as the sample is cooled.
3.4 Simulating microstructural evolution
The objective in this section is to theoretically compute the evolution of
martensite plate size distribution. As mentioned earlier, formation of marten-
site plates partitions the parent austenite grain into several compartments.
Untransformed regions of austenite naturally decrease in size with more
plates of martensite forming. Thus with the progress of transformation, new
martensite plates form with reduced plate lengths generating a distribution
of plate sizes. This concept, as proposed originally by Fischer [1949], has
been used here to calculate the theoretical distribution of martensite plate
lengths:
At n
th
instant (n varying as 1, 2, 3, . . .), the number of austenite pockets
available for transformation, N, is given by
N = 2
n1
(3.1)
This is taken as the number of martensite plates forming at that instant,
assuming only one martensite plate forms from a given pocket.
The volume of each of the austenite pockets at the n
th
instant, V
n
, is
56 3.4 Simulating microstructural evolution
Figure 3.9: The distribution of apparent martensite plate lengths (T
and L
values stated). The distributions are all on the same scale to allow compar-
isons to be made. The bin size on the horizontal scale is 5 m.
57 3.4 Simulating microstructural evolution
given by
V
n
=
V
0
2
n1
(3.2)
where V
0
is the original volume of austenite grain. The volume of austenite
pockets can be converted into mean lineal intercept, invoking the relation-
ships L = 1.69a and V = 11.314a
3
for truncated octahedral grains, where L
is the mean lineal intercept, V is the volume and a is the edge length of the
unit cell [Mark, 1956].
This results in the following equation relating L
n
to V
n
at n
th
instant:
L
n
= (0.4266V
n
)
1
3
(3.3)
This can be considered to be the length of any martensite plate forming at
that instant.
There is, however, one limitation of the above formulation. Volume
of the pre-existing martensite plates has not been taken into account while
calculating the volume of austenite available for transformation, V
n
, at the
n
th
instant. Shape of the martensite plates are generally considered to be
thin circular discs with a nite volume given by r
2
c; 2r being the length of
the plate and c, the semi-thickness, which is assumed to be 0.05r. The total
martensite plate volume is thus subtracted from the instantaneous austenite
pocket volume in order to calculate the eective austenite volume.
Knowing the initial mean lineal intercept of the austenite grain size,
it is now possible to obtain a theoretical variation in the martensite plate
length using the above methodology. The theoretical and the experimental
plate length distributions obtained with dierent prior austenite grain sizes
58 3.4 Simulating microstructural evolution
are presented in Figure 3.10.
The calculated and measured curves compare well, in terms of the in-
uence of the initial austenite grain size and with respect to their general
shapes. However, the absolute reduction in length is much greater with the
calculated curves. This is because equation 3.1 implies that the number of
plates N that form in each generation n is given by N = a
n1
with a = 2.
This leads to a rapid increase in N and a correspondingly rapid decrease in
the plate length as n increases. Thus, whereas each point in the calculated
curve can represent N plates of identical length, every point in the measured
curve is a dierent plate. The gradual variation in the measured plate length
distribution implies that a < 2.
59 3.4 Simulating microstructural evolution
(a)
(b)
Figure 3.10: (a) The experimentally measured apparent lengths of martensite
plates arranged in the sequence in which they are assumed to form. (b) The
corresponding curves for the calculated length versus the order in which the
plates form.
60 3.5 Martensite in TRIP-assisted steels
3.5 Martensite in TRIP-assisted steels
In TRIP-assisted steels, austenite is retained because of its carbon enrich-
ment following the bainite reaction. Partial decomposition of the untrans-
formed austenite into martensite during cooling after bainite reaction is not
unlikely. This is because the distribution of carbon in untransformed austen-
ite after the bainite reaction is mostly heterogeneous [Matas and Hehemann,
1961]. Areas near the bainitic ferrite-austenite interface are richer in carbon
than those away from it. The austenite away from the interfaces may trans-
form into martensite on cooling to ambient temperature. Indeed there is
evidence of martensite, albeit traces, formed during heat treatment of these
materials [Jacques et al., 2001d; Zaeerer et al., 2004]. The austenite in
TRIP-assisted steels is in a nely divided state with a grain size of about
1 m [Jacques, 2003]. It therefore follows that martensite forming in these
ne regions of austenite is likely to be devoid of microcracks.
Similarly, the high-carbon martensite produced during plastic deforma-
tion of TRIP-assisted steels can also be expected to resist cracking because
of the ne plate size. However, the cracking of martensite plates discussed
in Section 3.3 results from the quench stresses. The nature of these stresses
is arbitrary and has not been characterised.
Further studies were thus carried out with the following objectives:
To investigate the cracking tendency of the martensite plates during
controlled straining.
To study the nature of the martensite formed during deformation of
austenite.
61 3.6 Experiments
3.6 Experiments
Material
The alloy used for the experiments described in Section 3.3 was modied
with nickel to depress the M
S
temperature below room temperature. A
model developed by Capdevila et al. [2003] was used to study the eect
of adding nickel on the M
S
temperature of Fe-1C-1.5Mn-1.5Si (wt%) alloy
(Figure 3.11). From this, it appears that about 8 wt% nickel can depress
the M
S
temperature of the alloy suciently below room temperature. The
alloy can therefore be expected to be fully austenitic when quenched to room
temperature after austenitising at a higher temperature.
-60
-50
-40
-30
-20
-10
0
7 8 9 10 11
M
s
t
e
m
p
e
r
a
t
u
r
e
/
o
C
Ni content / wt%
Figure 3.11: Calculated eect of Ni addition on the M
S
temperature of Fe-
1C-1.5Mn-1.5Si (wt%) alloy.
The chemical composition of the modied alloy was analysed to be Fe-
0.91C-8.42Ni-1.58Mn-1.55Si-1.15Cr-0.23Mo (wt%). Calculations using MT-
DATA [2005] with the SGTE database reveal austenite to exist as a single
62 3.6 Experiments
phase at 1100
C. The
crosshead velocity used was 0.0017 mm s
1
. The martensitic samples were
compressed at room temperature to two dierent strain levels viz. 3% and
7%. A strain of about 30% was applied in case of the austenitic sample
compressed at a temperature of -501
C.
64 3.7 Results and discussion
Metallography
The compressed samples were sectioned along the vertical direction. These
sections were then polished, etched with 2% nital and observed using optical
microscopy.
3.7 Results and discussion
Figure 3.13 depicts the microstructure of the sample quenched in liquid ni-
trogen. The micrograph reveals dark-etching lenticular plates of martensite
along with untransformed austenite. This is similar to the microstructures
presented in Figures 3.1-3.4. However, microcracking of the martensite plates
of the type illustrated in Figures 3.5-3.6, was not observed in Figure 3.13.
This is despite the comparable scales of the microstructures. It is possible
that quenching in liquid nitrogen from room temperature does not induce
enough stress to crack the martensite plates, unlike quenching in water after
soaking at high temperatures. The addition of nickel might also be thought
to have reduced the characteristic brittleness of the high-carbon martensite
formed in the present alloy.
Martensitic samples produced by quenching in liquid nitrogen were fur-
ther subjected to compression tests. Figure 3.14 presents the microstructures
of the compressed samples. The plates of martensite are seen to contain
arrays of parallel cracks across the plate length. This is similar to the mi-
crocracks observed in Figure 3.5. Furthermore, the larger plates are found
to be preferentially cracked over the smaller ones (Figure 3.14(a)). This is
consistent with the earlier quenching experiments that revealed a greater
65 3.7 Results and discussion
Figure 3.13: Plates of athermal martensite formed without any microcracks
in the present alloy.
propensity of the coarser martensite plates to crack (Section 3.4). Cracks
were more frequent following higher plastic deformation (Figure 3.14(b)).
Sometimes cracks were also visible along some portions of the martensite-
austenite interfaces in the 7% strained sample. The apparent distance be-
tween adjacent cracks was also measured in these samples, Figure 3.15. This
is found to follow a similar distribution as that observed in Figure 3.8.
Figure 3.16 presents the micrographs obtained from the sample
quenched at room temperature and then compressed by 30% at -501
C.
Figure 3.16(a) reveals clearly the dark-etching small lens shaped martensite
plates formed during straining of the austenite. Strain-induced martensite
retained the lenticular shape but the individual plates are much ner in size
than expected from the prior austenite grain size. Martensite plates cannot
traverse across strong obstacles like grain boundaries, twin boundaries or
other pre-existing martensite plates. In the present situation, growth of the
66 3.7 Results and discussion
(a)
(b)
Figure 3.14: Cracking of martensite plates after (a) 3% and (b) 7% reduction
in height of the martensitic samples.
martensite plates appears to be stied by the defects produced due to plastic
strain in austenite. Hence the plates tend to grow only upto limited sizes.
67 3.7 Results and discussion
(a)
(b)
Figure 3.15: The frequency of observing cracks versus the apparent distance
between cracks on martensite plates after (a) 3% and (b) 7% reduction in
height.
68 3.7 Results and discussion
The plates were also observed to be free of cracks.
A notable feature of the microstructure is the particular alignment of
the martensite plates along the two preferred directions (Figure 3.16). This is
due to the eect of stress on transformation [Patel and Cohen, 1953]. Marten-
site is known to form on some specic crystallographic planes of austenite
known as habit planes. Without any external stress, transformation takes
place on all the possible habit planes. This therefore does not lead to any
bias in the microstructure. By contrast, martensite tends to form on only
some specic habit planes under the inuence of stress. This explains the
particular alignment of the martensite plates along some preferred directions
in the micrographs presented in Figure 3.16.
Figure 3.16(b) depicts that martensite is formed in one of the two grains
visible in the micrograph, while there is no sign of transformation in the other
grain. This may be due to the role of stress as described above. The resolved
shear component of the external stress on the habit planes actually leads to
transformation in steels, where the volume expansion due to transformation
is small. Transformation becomes possible only above some critical magni-
tude of the resolved shear stress. The orientation of the individual austenite
grains must inuence the magnitude of the resolved shear stress on the crys-
tallographic planes. Therefore it is not unlikely that some grains may not
be favourably oriented for transformation under the inuence of stress, as
observed in Figure 3.16(b). The other possibility could be the heterogeneous
nature of the compressive strain. All the grains might not be strained to
the same extent under compression. This may also be responsible for no
transformation to occur in one of the two grains visible in Figure 3.16(b).
69 3.8 Summary
3.8 Summary
High-carbon martensite is known to be a brittle phase. Yet it does not ap-
pear to be a liability in the TRIP-assisted steels reputed for good mechanical
properties, especially formability. This apparent anomaly can be explained
considering the role of martensite plate size on the brittleness. It has been
demonstrated that the ner martensite plates are better able to resist crack-
ing under stress or strain. It is the ne grain size of austenite that allows the
high-carbon brittle martensite to form in these materials without endanger-
ing the mechanical properties.
It has been also shown that the plates of martensite formed during
straining of austenite remain ne in size. This is probably due to the defects,
created in austenite due to strain, that stie the growth of martensite. The
small plates of strain-induced martensite are observed to form only along
some preferred directions. This is due to the eect of the accompanying
stress. The ne-sized ordered martensite plates that form during the strain-
ing of austenite can be expected to be resistant to cracking. Untempered
martensite that forms in the TRIP-assisted steels is therefore not expected
to have any adverse eect on the mechanical performance, despite its high
carbon concentration.
70 3.8 Summary
(a)
(b)
Figure 3.16: Martensite plates formed during straining of the austenitic sam-
ple; (a) showing preferred alignments of the plates, (b) revealing the prefer-
ential transformation in one of the two adjacent grains of austenite.
Chapter 4
Microstructural Modelling
Allotriomorphic ferrite, carbide-free bainite and retained austenite are con-
stituents of the microstructure of TRIP-assisted steels. These phases can be
produced beginning either with hot rolled or cold rolled steel. The austen-
ite after hot rolling transforms to allotriomorphic ferrite and subsequently
to bainite. However, a two-stage annealing process is necessary after cold
rolling to produce the desired microstructure. During bainite formation, un-
transformed austenite gets enriched with carbon and is therefore retained at
ambient temperature. Processing conditions as well as the chemical compo-
sition of the steel inuences both the volume fraction and carbon content of
retained austenite. The present chapter deals with an attempt to model the
microstructure of TRIP-assisted steels using neural network techniques.
4.1 Neural network modelling
Neural network analysis represents a general method of regression that ts
a non-linear function to experimental data. This can help quantify complex
71
72 4.1 Neural network modelling
phenomena involving many variables. The subject has been reviewed by
Bhadeshia [1999a]. In linear regression, the product of each input variable
x
i
and a weight w
i
is added to a constant to obtain the output y,
y =
i
w
i
x
i
+ (4.1)
The weights and the constants can be imagined to be varied systematically
until a best-t expression of the output is reached as a function of the input
variables. The diculties with such a regression are:
1. A specic relation has to be chosen prior to the analysis.
2. The regression equation thus derived applies across the entire input
space.
These diculties can be overcome by making the network non-linear,
setting the sum of the products to the argument of a hyperbolic tangent
function,
y =
j
W
j
h
j
+
, with h
j
= tanh(
i
w
ij
x
i
+
j
) (4.2)
where W
j
is a weight and
n
(t
n
y
n
)
2
(4.3)
If all the important input variables are not considered in the analysis,
E
D
tends to increase. E
D
represents only an overall perceived level of noise
in the output and does not describe the uncertainties of the modelling pro-
cess itself. The best-t functions with the most probable weights may not
necessarily describe the uncertainties in regions of the input space with noisy
or sparse data.
Neural networks used in Bayesian framework allow the calculation of
error bars representing the uncertainty in the tting parameters [MacKay,
2003]. This invokes the possibility of tting or extrapolating many functions
into uncertain regions of the input space, still maintaining an optimum t
in regions having good quality data. A probability distribution of sets of
weights is then used to dene the tting uncertainty. It follows that the
error bars tend to increase in magnitude when data are noisy (A) or sparse
(B), as illustrated in Figure 4.2.
75 4.2 Model for retained austenite fraction
Figure 4.2: The uncertainty in dening a tting function when data are
noisy (A) or sparse (B). The thinner lines represent error bounds due to
uncertainties in determining the weights [Bhadeshia, 1999a].
4.2 Model for retained austenite fraction
4.2.1 Database
The cold rolled steel is initially heated to a temperature in the ( + )
phase region (intercritical annealing), to form ferrite and austenite. The
material is then cooled rapidly to a lower temperature at which austenite
transforms to bainite. Parameters of this two-step heat treatment along
with the chemical composition were taken as input variables to the model, as
shown in Table 4.1
1
. Figure 4.3 presents the distribution of the data points
in the input space for each variable.
1
Data used are due to Chen et al. [2002]; De Meyer et al. [1999a]; Itami et al. [1995];
Jacques et al. [2001a,b,c, 1999]; Kim et al. [2002, 2001]; Lee et al. [2002b]; Matsumura et
al. [1992, 1987b]; Nakagaito et al. [2003]; Pichler et al. [1998]; Sakuma et al. [1993, 1991a];
Traint et al. [2002, 2000]
76 4.2 Model for retained austenite fraction
Inputs Min. Max. Average S.D.
Carbon / wt% 0.0950 0.3920 0.2102 0.1020
Manganese / wt% 0.6000 1.9900 1.4057 0.2796
Silicon / wt% 0.0400 2.1000 0.9619 0.5369
Aluminium / wt% 0.0000 2.0000 0.2979 0.5236
Phosphorus / wt% 0.0000 0.2040 0.0223 0.0401
Molybdenum / wt% 0.0000 0.1400 0.0043 0.0218
Copper / wt% 0.0000 0.5100 0.0065 0.0573
IA temperature /
C 730 840 778.2653 25.7829
IA time / s 60 450 213.2143 110.8414
IT temperature /
C 300 525 404.2347 39.5554
IT time / s 0 5400 441.3673 588.7228
Output Min. Max. Average S.D.
Retained austenite / vol.% 0.1267 27.0710 10.3805 5.8054
Table 4.1: Data used for the retained austenite fraction model, IA: intercrit-
ical annealing and IT: isothermal transformation to bainite, S.D.: standard
deviation.
4.2.2 Model characteristics
The model was trained using 80% of the dataset values. The remainder
was used to test the performance of the model with unseen data. Figure 4.4
compares the normalised predicted values with the corresponding normalised
database value during training and testing steps respectively. Normalisation
of the variables is carried out as follows
x
n
=
x x
min
x
max
x
min
0.5 (4.4)
where x
n
is the normalised value of the variable x, with x
min
and x
max
as
the minimum and the maximum values respectively. The performance of the
model can be seen to be satisfactory both during the training and testing
stages.
77 4.2 Model for retained austenite fraction
0
0.1
0.2
0.3
0.4
0 10 20 30
c
a
r
b
o
n
/
w
t
%
retained austenite / vol.%
(a)
0
0.5
1
1.5
2
0 10 20 30
m
a
n
g
a
n
e
s
e
/
w
t
%
retained austenite / vol.%
(b)
0
0.5
1
1.5
2
0 10 20 30
s
i
l
i
c
o
n
/
w
t
%
retained austenite / vol.%
(c)
0
0.5
1
1.5
2
0 10 20 30
a
l
u
m
i
n
i
u
m
/
w
t
%
retained austenite / vol.%
(d)
0
0.05
0.1
0.15
0.2
0 10 20 30
p
h
o
s
p
h
o
r
o
u
s
/
w
t
%
retained austenite / vol.%
(e)
0
0.05
0.1
0.15
0 10 20 30
m
o
l
y
b
d
e
n
u
m
/
w
t
%
retained austenite / vol.%
(f)
0
0.25
0.5
0 10 20 30
c
o
p
p
e
r
/
w
t
%
retained austenite / vol.%
(g)
720
760
800
840
0 10 20 30
I
A
t
e
m
p
e
r
a
t
u
r
e
/
o
C
retained austenite / vol.%
(h)
0
100
200
300
400
500
0 10 20 30
I
A
t
i
m
e
/
s
retained austenite / vol.%
(i)
250
350
450
550
0 10 20 30
I
T
t
e
m
p
e
r
a
t
u
r
e
/
o
C
retained austenite / vol.%
(j)
0
2000
4000
6000
0 10 20 30
I
T
t
i
m
e
/
s
retained austenite / vol.%
(k)
Figure 4.3: Distribution of data in the input space for the retained austenite
fraction model
78 4.2 Model for retained austenite fraction
Figure 4.5(a) illustrates the reduction in noise with increasing number
of hidden units. The nal model is formed as a committee of models, which is
selected on the basis of the observed minimum in the test error on combining
increasing numbers of networks (Figure 4.5(b)).
Figure 4.6 presents the signicances of each of the input variables as
perceived by the models in the nal committee. Bainite reaction time appears
to be the most important factor controlling the retained austenite content in
the microstructure.
4.2.3 Model predictions
Predictions were made by varying input variables over a base steel compo-
sition of 0.2C-1.5Mn-1.0Si-0.5Al-0.02P-0Mo-0Cu (all in wt%), intercritical
annealing parameters of 780
C and 375
C.
85 4.3 Model for carbon in retained austenite
2
3
4
5
6
7
8
9
10
11
12
0 300 600 900
r
e
t
a
i
n
e
d
a
u
s
t
e
n
i
t
e
/
v
o
l
.
%
bainite reaction time / s
(a)
2
3
4
5
6
7
8
9
10
11
12
0 300 600 900
r
e
t
a
i
n
e
d
a
u
s
t
e
n
i
t
e
/
v
o
l
.
%
bainite reaction time / s
(b)
Figure 4.9: Variation in retained austenite volume fraction with bainite re-
action time at transformation temperature of (a) 350
C.
4.3 Model for carbon in retained austenite
4.3.1 Database
Data used for the model made for predicting the carbon content of retained
austenite are presented in Table 4.2
2
. Figure 4.10 presents the distribution
of each of the input variables against the output.
4.3.2 Model characteristics
As expected, increasing the number of hidden units reduces the noise level
(Figure 4.11(a)). However, the test error does not tend to decrease contin-
uously with increasing model complexity (Figure 4.11(b)). The optimum
number of models in the committee was found to be seven, when the test
error is minimum. The model appears to perform reasonably well, during
2
Data used are due to De Meyer et al. [1999a]; Itami et al. [1995]; Jacques et al.
[2001a,b,c, 1999]; Kim et al. [2002]; Matsumura et al. [1992, 1987b]; Nakagaito et al.
[2003]; Wang et al. [2004]; Zaeerer et al. [2004]
86 4.3 Model for carbon in retained austenite
Inputs Min. Max. Average S.D.
Carbon / wt% 0.1100 0.3920 0.2676 0.1150
Manganese / wt% 0.8300 1.6600 1.3019 0.1813
Silicon / wt% 0.3400 1.9900 1.2150 0.5284
Aluminium / wt% 0.0000 1.5000 0.1879 0.3982
Phosphorus / wt% 0.0000 0.0400 0.0025 0.0058
Copper / wt% 0.0000 0.5100 0.0051 0.0506
IA temperature /
C 730 830 773.2778 25.9773
IA time / s 60 360 286.9697 50.8797
IT temperature /
C 350 450 392.6263 28.4305
IT time / s 0 9725 578.7222 991.3905
Output Min. Max. Average S.D.
Carbon in retained austenite / wt% 0.49 1.8793 1.1949 0.3836
Table 4.2: Data used to model the carbon content of retained austenite, IA:
intercritical annealing and IT: isothermal transformation to bainite, S.D.:
standard deviation.
both the training and testing step (Figure 4.12(a)-(b)).
Figure 4.13 shows the model-perceived signicance of each of the input
parameters. Bainite reaction time appears to exert the strongest inuence
on the carbon content of retained austenite.
4.3.3 Model predictions
Predictions were made varying each input parameter over a base composition
of 0.2C-1.5Mn-1.0Si-0.5Al-0.02P-0Cu (wt%); intercritical annealing parame-
ters of 780
C, 500 s.
Figures 4.14(a)-(d) demonstrate that the results for carbon, aluminium and
silicon are uncertain while manganese shows an inverse trend. Manganese
shifts the T
0
curve to lower carbon contents and hence reduces the carbon
content of retained austenite.
87 4.3 Model for carbon in retained austenite
0.1
0.2
0.3
0.4
0.4 0.8 1.2 1.6 2
c
a
r
b
o
n
/
w
t
%
carbon in ret. aust. / wt%
(a)
0.9
1.2
1.5
1.8
0.4 0.8 1.2 1.6 2
m
a
n
g
a
n
e
s
e
/
w
t
%
carbon in ret. aust. / wt%
(b)
0.3
0.9
1.5
2.1
0.4 0.8 1.2 1.6 2
s
i
l
i
c
o
n
/
w
t
%
carbon in ret. aust. / wt%
(c)
0
0.4
0.8
1.2
1.6
0.4 0.8 1.2 1.6 2
a
l
u
m
i
n
i
u
m
/
w
t
%
carbon in ret. aust. / wt%
(d)
0
0.02
0.04
0.4 0.8 1.2 1.6 2
p
h
o
s
p
h
o
r
o
u
s
/
w
t
%
carbon in ret. aust. / wt%
(e)
0
0.25
0.5
0.4 0.8 1.2 1.6 2
c
o
p
p
e
r
/
w
t
%
carbon in ret. aust. / wt%
(f)
720
760
800
840
0.4 0.8 1.2 1.6 2
I
A
t
e
m
p
e
r
a
t
u
r
e
/
o
C
carbon in ret. aust. / wt%
(g)
0
100
200
300
400
0.4 0.8 1.2 1.6 2
I
A
t
i
m
e
/
s
carbon in ret. aust. / wt%
(h)
350
400
450
0.4 0.8 1.2 1.6 2
I
T
t
e
m
p
e
r
a
t
u
r
e
/
o
C
carbon in ret. aust. / wt%
(i)
0
2500
5000
7500
10000
0.4 0.8 1.2 1.6 2
I
T
t
i
m
e
/
s
carbon in ret. aust. / wt%
(j)
Figure 4.10: Distribution of data in the input space for the carbon content
of retained austenite model.
88 4.3 Model for carbon in retained austenite
N
o
i
s
e
0.04
0.08
0.12
0.16
0 4 8 12 16 20
Number of hidden units
(a)
0.4
0.8
1.2
1.6
0 4 8 12 16 20
Number of models
C
o
m
b
i
n
e
d
T
e
s
t
E
r
r
o
r
(b)
Figure 4.11: (a) Variation of noise () with increasing number of hidden
units. (b) Variation of test error with number of models in the committee.
Carbon may also be expected to have a similar eect as manganese.
But this is apparently missing in Figure 4.14(a). The concentration of car-
89 4.3 Model for carbon in retained austenite
0.6
0.4
0.2
0
0.2
0.4
0.6
0.6 0.4 0.2 0 0.2 0.4 0.6
Database values, normalised
P
r
e
d
i
c
t
e
d
v
a
l
u
e
s
,
n
o
r
m
a
l
i
s
e
d
(a)
0.6
0.4
0.2
0
0.2
0.4
0.6
0.6 0.4 0.2 0 0.2 0.4 0.6
Database values, normalised
P
r
e
d
i
c
t
e
d
v
a
l
u
e
s
,
n
o
r
m
a
l
i
s
e
d
(b)
Figure 4.12: Model predictions compared with the actual database values
during (a) training and (b) test.
bon in austenite before bainite formation is not the same as that of the
bulk steel composition. During intercritical annealing, carbon in austenite
90 4.3 Model for carbon in retained austenite
Figure 4.13: Signicances of the input parameters (A-J) to the models for
predicting carbon content of retained austenite, A: carbon, B: manganese,
C: silicon, D: aluminium, E: phosphorous, F: copper, G: intercritical anneal-
ing temperature, H: intercritical annealing time, I: bainite transformation
temperature, J: bainite transformation time.
increases due to ferrite formation. However, manganese in the austenite
formed during intercritical annealing is most likely to remain same as that
of the bulk composition. This is because partitioning of manganese requires
a longer time than partitioning of carbon. It seems that carbon enrichment
of retained austenite is aected by the bulk manganese concentration but
not bulk carbon content of the steel. It may therefore be inferred that car-
bon enrichment in retained austenite is aected by the factors that can alter
bainite reaction kinetics. This is also apparent from the eect of the heat
treatment parameters on the carbon content of retained austenite, as illus-
91 4.3 Model for carbon in retained austenite
1
1.25
1.5
1.75
2
2.25
2.5
2.75
0.1 0.2 0.3 0.4 0.5
c
a
r
b
o
n
/
w
t
%
carbon / wt%
(a)
1
1.4
1.8
2.2
2.6
0 0.5 1 1.5 2 2.5
c
a
r
b
o
n
/
w
t
%
silicon / wt%
(b)
1.2
1.4
1.6
1.8
2
2.2
2.4
0 0.5 1 1.5 2
c
a
r
b
o
n
/
w
t
%
aluminium / wt%
(c)
0.5
1
1.5
2
2.5
3
3.5
0 0.5 1 1.5 2 2.5
c
a
r
b
o
n
/
w
t
%
manganese / wt%
(d)
Figure 4.14: Eect of alloying elements on the carbon content of retained
austenite: (a) carbon, (b) silicon, (c) aluminium and (d) manganese.
trated in Figure 4.15.
Austenite formed at a higher intercritical annealing temperatures is
leaner in carbon. This can accelerate the subsequent transformation to bai-
nite, resulting in rapid carbon enrichment of the residual austenite (Fig-
ure 4.15(a)). Similarly, with a longer holding period during the intercritical
annealing, austenite gets enriched more with carbon and manganese, hence
gaining stability. Subsequent transformation of austenite to bainite becomes
sluggish, which explains the lower carbon enrichment at longer holding times
of intercritical annealing (Figure 4.15(b)).
92 4.3 Model for carbon in retained austenite
0.7
1.1
1.5
1.9
2.3
2.7
3.1
720 760 800 840 880
c
a
r
b
o
n
/
w
t
%
intercritical annealing temperature /
o
C
(a)
0.5
1
1.5
2
2.5
100 150 200 250 300 350
c
a
r
b
o
n
/
w
t
%
intercritical annealing time / s
(b)
0.2
0.6
1
1.4
1.8
2.2
300 350 400 450 500
c
a
r
b
o
n
/
w
t
%
bainite reaction temperature /
o
C
(c)
0.4
0.8
1.2
1.6
2
2.4
2.8
0 500 1000 1500 2000 2500
c
a
r
b
o
n
/
w
t
%
bainite reaction time / s
(d)
Figure 4.15: Eect of heat treatment conditions on the carbon content of
retained austenite: intercritical annealing (a) temperature and (b) time; bai-
nite transformation (c) temperature and (d) time.
Bainite transformation at a lower temperature should result in more
carbon content of the retained austenite (as expected from T
0
criterion).
This is observed up to a particular temperature, below which the carbon
content appears to fall, as can be seen in Figure 4.15(c). At low tempera-
ture, transformation kinetics can be slow. Therefore the progress of bainite
formation may not be sucient to enrich the retained austenite with car-
bon. With an increase in bainite reaction time, carbon in retained austenite
increases initially (Figure 4.15(d)). However, prolonged holding times may
93 4.4 Summary
result in some marginal loss of carbon, possibly due to cementite precipita-
tion. Predictions for longer holding times are however associated with large
uncertainties, limiting condence in the predicted behaviour. More experi-
ments are necessary in this region to conrm the predicted behaviour and its
cause.
Figure 4.16 depicts the variation of carbon in retained austenite with
progress of bainite reaction at two dierent temperatures for 0.2C-1.5Mn
(wt%) steel made with either 1.5 wt% silicon or aluminium. Carbon in re-
tained austenite appears to attain a maximum value at an intermediate time
for both the steels at both temperatures. The maximum carbon content is
higher at the lower transformation temperature for each steel. Furthermore,
the maximum carbon content is higher for the aluminium-containing steel
than that made with silicon. Both these predictions can be accounted for by
the T
0
criterion of bainite transformation theory and can therefore be judged
to be reasonable.
4.4 Summary
Neural network techniques have been adopted to create models for predict-
ing the retained austenite fraction and its carbon content as a function of
processing parameters and composition of TRIP-assisted steels.
Elements that tend to stabilise austenite over ferrite appear to increase
the nal retained austenite content. Retention of austenite is due to the car-
bon enrichment of the residual austenite after bainite formation. Therefore
the rate of bainite formation seems to have a strong inuence on retaining
higher austenite fraction. Any factor that alters the rate of bainite transfor-
94 4.4 Summary
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
0 300 600 900
c
a
r
b
o
n
c
o
n
t
e
n
t
/
w
t
%
bainite reaction time / s
(a)
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
0 300 600 900
c
a
r
b
o
n
c
o
n
t
e
n
t
/
w
t
%
bainite reaction time / s
(b)
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
0 300 600 900
c
a
r
b
o
n
c
o
n
t
e
n
t
/
w
t
%
bainite reaction time / s
(c)
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
0 300 600 900
c
a
r
b
o
n
c
o
n
t
e
n
t
/
w
t
%
bainite reaction time / s
(d)
Figure 4.16: Variation in carbon content of retained austenite with bainite re-
action time for 0.2C-1.5Mn-1.5Si (wt%) steel at transformation temperature
of (a) 375
C.
mation aects retained austenite content in these steels.
The carbon content of retained austenite is primarily controlled by the
extent of bainite transformation. However, the maximum carbon enrichment
of the retained austenite can also be inuenced by altering the T
0
point, with
suitable alloying additions or the bainite formation temperature.
The model for retained austenite fraction was further used to formulate
an optimum combination of silicon and retained austenite contents in these
steels, as explained in the following Chapter.
Chapter 5
-TRIP steel
In TRIP-assisted steels, silicon helps to retain austenite by suppressing the
precipitation of cementite. The carbon partitioned into the austenite, follow-
ing the formation of bainite, remains there and permits it to be retained at
room temperature. However, there is no quantitative theory to estimate the
exact concentration of silicon necessary to prevent cementite. The silicon
can be detrimental to the surface properties of automobile steels. So it is
important to optimise both the silicon and the retained austenite contents.
The work presented in this chapter began with the aim of nding an
optimum combination of silicon and retained austenite but the computa-
tional methods used gave surprising and interesting results. This stimulated
experiments which led to the discovery of a new concept in TRIP-assisted
steels.
To describe the story in a logical manner, the optimisation technique is
explained rst. This exploits the neural network model presented in Chapter
4. The nal part deals with the design and testing of a steel with a radically
dierent microstructure.
95
96 5.1 Optimisation
5.1 Optimisation
Numerical models can be used to optimise the chemical composition and
process parameters of a material in order to achieve the appropriate mi-
crostructure or properties [Murugananth et al., 2004]. This is done by the
repeated interrogation of the numerical model in order to attain a specic
output, in a process called optimisation (Figure 5.1).
The rst step is to identify the target, the design variables i.e., the
input variables of the numerical model that are to be varied and the param-
eters which should not be altered during the calculations. An initial set of
input variables is chosen and is used to evaluate the output parameter. The
dierence between the calculated output and the target is then used as a
tness for purpose parameter. If the tness is not satisfactory, the optimiser
algorithm decides the direction in which the design variables should be var-
ied and the whole process is then repeated. The process terminates once the
required degree of tness is attained.
Optimisation methods are classied as linear or non-linear, based on
the function relating the variables of the model [Nelder and Mead, 1965].
In practice, linear functions are seldom observed. Sequential quadratic pro-
gramming, downhill simplex methods, genetic algorithms represent some of
the non-linear methods of optimisation [Epogy, 2003].
Genetic algorithm
Genetic algorithms carry out the optimisation by mutation or recombination
and selection [Epogy, 2003]. A point is rst chosen in the design space.
A family of points is subsequently generated around the initial one. This
97 5.1 Optimisation
Figure 5.1: Flow diagram of the optimisation process.
is called mutation. By contrast, in the recombination process, values are
exchanged between a random number of points. These are done in such a
98 5.2 Optimised TRIP-assisted steel
way so that the points move in a direction allowing the objective function to
be optimised. The best of all the points is nally selected as the solution.
5.2 Optimised TRIP-assisted steel
A genetic algorithm was used in the present work to formulate an optimum
chemical composition of TRIP-assisted steel, given the objective of maximis-
ing the retained austenite fraction and keeping the silicon addition to a min-
imum. A neural network model was created relating the retained austenite
fraction to the concentrations of dierent alloying elements and the process
parameters of these materials, as described in Chapter 4. The model was
coupled with commercial software, Epogy [2003], to carry out the optimi-
sation process. The coupling and optimisation was carried out at Nanyang
Technological University, Singapore by Dr M. Murugananth.
Concentrations of all the alloying elements, except phosphorous, were
allowed to vary. The intercritical annealing temperature and the isothermal
bainite reaction temperature were also chosen as design variables. However,
the dwell time at these temperatures was kept constant. The uncertainty
in the model prediction was considered as a constraint, with a maximum
permitted value of 20%.
The algorithm converged to an alloy with the chemical composition
Fe-0.4C-2Mn-0.5Si-2Al-0.5Cu-0.02P (wt%), with a microstructure estimated
to contain 4120 vol.% of retained austenite, after intercritical annealing at
840
C.
The alloy is rich in carbon, manganese, aluminium and copper. Each of
these elements is known to help retain austenite [Girault et al., 2001; Jacques
99 5.3 Thermodynamic calculations
et al., 2001a; Kim and Lee, 1999; Sakuma et al., 1991a; Shi et al., 2002]. The
concentrations of these elements were therefore increased by the algorithm
while reducing the silicon content, as required by the objective. The result
obtained therefore appears to be physically meaningful.
5.3 Thermodynamic calculations
Amounts of various phases in equilibrium as a function of temperature were
calculated for the optimised alloy using MTDATA with the SGTE database.
Ferrite, austenite, cementite and liquid were allowed to exist. The results are
presented in Figure 5.2. The alloy starts to solidify by forming ferrite, which
reaches a maximum of about 85 wt%. This is called -ferrite. Austenite
starts to form thereafter with the amount of -ferrite decreasing. However,
the calculations indicate that austenite can never exist as a single phase in
the alloy. With cooling, the amount of austenite increases and attains a
maximum of about 80 wt% before it starts to decompose to allotriomorphic
ferrite.
Figure 5.2 also reveals the alloy to exist as a mixture of about 54 wt%
austenite and 46 wt% ferrite at 840
C and 160
C respectively. Both
the predicted intercritical annealing and the bainite reaction temperature of
300
C,
for the austenite formed during intercritical annealing at 840
C.
The initial mole fraction of carbon in a fully austenitic sample, x, can
be related to the nal volume fractions and carbon content of the bainitic
101 5.3 Thermodynamic calculations
ferrite and austenite as,
x = v
+v
(5.1)
where v
and v
and x
respectively.
The mole fraction of carbon in the austenite formed at 840
C is cal-
culated to be 0.0312. Substituting x = 0.0312, x
C for 1 h, 2 h and 4 h
respectively. Samples of the alloys A and B were also subjected to two-step
heat treatment, involving initial heating to a temperature in the (+) phase
region followed by an isothermal dwell at a lower temperature, in order to
103 5.4 Experiments
obtain bainite. Details about the heat treatment parameters are described
in the following section.
Characterisation procedures
Microstructures of the samples in both as-cast and heat-treated conditions
were characterised using optical and eld emission gun scanning electron mi-
croscopy (FEGSEM). Samples were polished using standard metallographic
techniques and etched with 2% nital solution. The electron microscope was
operated at a 10 kV accelerating voltage.
X-ray diraction tests of the samples were carried out with a scan range
of 40
-120
at 0.5
radiation was
used and the operating voltage and current were set at 40 kV and 40 mA
respectively. Samples were polished, etched with 2% nital and then once
again just ne polished before nally putting into the diractometer.
Mechanical tests
The bulk hardnesses and the microhardness of the constituent phases were
recorded with 10 kg and 50 g load respectively.
Uniaxial tensile tests were carried out with the heat-treated samples
of the alloy B, at room temperature and 100
C.
Phase Amount / wt% Composition / wt%
C Mn Si Al P Cu
Ferrite 79.5 0.2 2.1 0.4 1.7 0.2 0.4
Liquid 20.5 1.2 3.7 0.7 1 0.2 0.4
Table 5.4: Calculated amounts and the chemical composition of the phases
in equilibrium for the alloy A at 1400
C.
The liquid coexisting with the -ferrite at 1400
C for varying
time periods, mainly to conrm the high carbon content of the martensitic
109 5.5 Results and discussion
residual phase. The microhardness of the phase was found to decrease drasti-
cally with increase in tempering time, as illustrated in Table 5.5. The average
values of 20 readings with the error as one standard deviation are reported
in Table 5.5.
Figure 5.6 shows the FEGSEM images obtained from a tempered sam-
ple. Compared to the as-cast condition (Figure 5.5), the residual phase after
tempering appears with better contrast (Figure 5.6). This might be due to
carbide precipitation from the martensite during tempering, increasing the
interfacial area per unit volume thereby resulting in more intense etching.
Tempering time / h Microhardness / HV
0 80157
1 65523
2 60630
4 52936
Table 5.5: Drop in microhardness of the residual martensitic phase with
tempering time at 300
C.
The two-step heat-treatment of the as-cast samples were carried out to
produce a microstructure comprised of -ferrite, bainite and austenite. Cal-
culations using MTDATA with SGTE database revealed alloy A to consist of
60 wt% austenite in equilibrium with 40 wt% ferrite at 810
C (Figure 5.3).
Samples of the alloy were therefore heated at 810
C for 4 h.
111 5.5 Results and discussion
ond phase, the nature of which is not clear. The FEGSEM images of the
sample are presented in Figure 5.8. It appears that the second phase is made
up of bainite and retained austenite. The presence of about 14.30.7 wt%
retained austenite is revealed by X-ray diraction analysis. This is, how-
ever, seen to be less than both the amount of retained austenite predicted by
the genetic algorithm and that calculated using the bainite transformation
theory. The carbon content of the retained austenite was estimated to be
1.350.02 wt% from the lattice parameter of austenite, using the equation
proposed by Dyson and Holmes [1979].
112 5.5 Results and discussion
(a)
(b)
Figure 5.7: Optical micrographs of alloy A heat-treated to form bainite.
113 5.5 Results and discussion
(a)
(b)
Figure 5.8: FEGSEM images of the heat-treated sample of alloy A revealing
bainite.
114 5.5 Results and discussion
Alloy B
As expected from Figure 5.9, the microstructure of the alloy in the as-cast
condition reveals the presence of -ferrite, embedded in a dark-etching resid-
ual phase (Figure 5.10(a)). The microhardness of the residual phase in this
alloy was recorded as 38222 HV and that of the -ferrite as 2148 HV.
Unlike alloy A, the residual phase of this material does not therefore appear
to be martensitic. Figure 5.10(b) reveals the residual phase to be pearlite.
0
20
40
60
80
100
750 1000 1250 1500
a
m
o
u
n
t
o
f
p
h
a
s
e
/
w
t
%
Temperature /
o
C
Liquid
Ferrite
Austenite
Cementite
Figure 5.9: Amounts of phases in equilibrium at dierent temperatures for
the alloy B, calculated using MTDATA with SGTE database.
Samples of alloy B were subjected to the two-step heat-treatment. Fig-
ure 5.9 reveals the alloy to exist as a mixture of 45 wt% ferrite and 55 wt%
austenite in equilibrium at 940
C for 3 h.
Figure 5.11 depicts the microstructures of the heat-treated sample. The
distinct morphology of the -ferrite is seen to have been retained after the
heat-treatment, along with the formation of bainite. The presence of retained
austenite in the heat-treated sample was conrmed with X-ray diraction
analysis, which revealed 13.50.9 wt% of retained austenite containing about
1.170.03 wt% carbon.
The microhardnesses of -ferrite and bainite along with the bulk hard-
ness of the heat-treated sample are presented in Table 5.6. Hardness data
and point counts were used to estimate the fractions of -ferrite and bainite,
as described previously. Assuming that the densities of -ferrite and the
bainitic phase are approximately equal, the measurements eectively rep-
resent the weight fraction of each phase, as presented in Table 5.7. The
fraction of bainite is seen to be in reasonable agreement with that of the
austenite present at the intercritical annealing temperature, calculated using
MTDATA (Figure 5.9).
Phase Hardness / HV
-ferrite 2679
Bainite 51913
Bulk 38540
Table 5.6: Microhardnesses of the constituent phases and the bulk hardness
of the heat-treated sample of alloy B.
Heat-treated samples of the alloy B were mechanically tested using uni-
axial loading. The properties observed at room temperature are impressive
with 23% total elongation and an ultimate tensile strength of 1000 MPa (Fig-
117 5.5 Results and discussion
(a)
(b)
Figure 5.11: Microstructures of the heat-treated sample of alloy B, (a) optical
micrograph and (b) FEGSEM image.
118 5.5 Results and discussion
Phase Weight fraction
from hardness from point count
-ferrite 0.53 0.430.02
Bainite 0.47 0.570.03
Table 5.7: Weight fractions of -ferrite and bainite phases in the heat-treated
sample of alloy B.
ure 5.12). The 0.2% proof stress was recorded as 700 MPa. There has been
no signicant necking of the sample at fracture suggesting that the elongation
has mostly been uniform (Figure 5.13).
The true stress-strain plot in Figure 5.12 shows steady strain-hardening
continuing until fracture. This could be due to the transformation of re-
tained austenite into martensite during straining. The X-ray analysis of the
sample cut from near the fracture surface of the tensile specimen indicated
7.00.8 wt% of austenite to exist. During straining, the austenite content
was therefore reduced to half of that in the undeformed condition. The trans-
formation of austenite into martensite during straining might have played a
key role inuencing the properties.
Figure 5.14 depicts the stress-strain curve recorded for the sample
tested at 100
C.
121 5.6 Summary
Alloys C and D
Figure 5.15 indicates both the alloys C and D can be expected to have a sim-
ilar phase transformation behaviour as that of the alloys A and B. However,
it is surprising to observe that none of these materials revealed the presence
of -ferrite in the as-cast condition. Figure 5.16 depicts the light-etching
allotriomorphic ferrite formed along the prior austenite grain boundaries in
alloys C and D. The morphology of the ferrite grains in these materials is in
sharp contrast to that of the elongated -ferrite grains, as observed in alloys
A and B.
The absence of -ferrite in as-cast microstructures of alloys C and D
could be due to the cooling pattern during solidication. Figure 5.15 suggests
that -ferrite forms during solidication up to around 1400
C. At lower tem-
peratures, austenite starts gaining stability over ferrite. A fast cooling rate
during solidication can bring down the temperature rapidly below 1400
C,
thereby avoiding the chance of forming -ferrite. This might have led to the
solidication of the alloys C and D forming austenite, which has transformed
partly to allotriomorphic ferrite during subsequent cooling.
5.6 Summary
TRIP-assisted steels are made with silicon, which helps retaining austenite.
However, silicon leads to many diculties during processing and deteriorates
the surface quality of the steel. A genetic algorithm coupled with a neural
network model yielded an optimised chemical composition which is low in
silicon content but rich in aluminium, manganese and copper.
122 5.6 Summary
Experiments to conrm the calculations gave surprising results. Some
of the alloys demonstrated an enhanced tendency of -ferrite formation during
solidication, as expected from phase diagram calculations. However, there
was no evidence of the presence of -ferrite in the as-cast microstructures of
some other alloys, which could be due to the dierence in cooling pattern.
The as-cast samples of the alloys containing -ferrite were heat-treated
to form microstructures containing -ferrite, bainite and retained austen-
ite. The heat-treated samples were found to possess impressive mechanical
properties at room temperature. However, increased stability of the retained
austenite at higher temperature might have been responsible for impairing
the ductility of the specimen tested at 100
C.
123 5.6 Summary
0
20
40
60
80
100
750 1000 1250 1500
a
m
o
u
n
t
o
f
p
h
a
s
e
/
w
t
%
Temperature /
o
C
Liquid
Ferrite
Austenite
Cementite
(a)
0
20
40
60
80
100
750 1000 1250 1500
a
m
o
u
n
t
o
f
p
h
a
s
e
/
w
t
%
Temperature /
o
C
Liquid
Ferrite
Austenite
Cementite
(b)
Figure 5.15: Amounts of phases in equilibrium at dierent temperatures
for the (a) alloy C and (b) alloy D, calculated using MTDATA with SGTE
database.
124 5.6 Summary
(a)
(b)
Figure 5.16: Optical micrographs from the as-cast alloys C and D showing
allotriomorphic ferrite along the prior austenite grain boundary.
Chapter 6
Mechanical Stabilisation
Austenite can be induced to transform to martensite above the M
S
tem-
perature by stress or strain. This induced transition is widely believed to
be the cause of the excellent mechanical properties, especially the formabil-
ity, of TRIP-assisted steels. Stress adds to the thermodynamic driving force
whereas plastic strain aects the kinetics of the process. Defects created in
austenite due to strain provide additional nucleation sites. However, growth
of martensite plates can be hindered by defects in austenite and excessive
deformation can even halt the transformation, in a phenomenon known as
mechanical stabilisation. The purpose of the present chapter is to study the
role of stress and strain on martensitic transformations.
6.1 Mechanical driving force
Martensitic transformation during deformation of austenite was rst discov-
ered by Scheil [1932], the thermodynamics of which was explained by Patel
and Cohen [1953]. Stress makes transformation possible above the M
S
tem-
125
126 6.1 Mechanical driving force
perature by supplying the shortfall of the necessary thermodynamic driving
force. This is illustrated schematically in Figure 6.1. The free energy change,
G
1
, available at the M
S
temperature is the critical driving force necessary
for transformation. At a temperature T
1
above M
S
, the magnitude of the free
energy change due to transformation, G
2
, is less than the critical driving
force, G
1
. Martensite can only be formed at this temperature by supply-
ing additional energy of a magnitude U, so that U + G
2
= G
1
. This
can be done by applying stress. The work done by external stress adds to
the available chemical free energy change as a mechanical driving force for
transformation at any temperature.
Figure 6.1: A shortfall in driving force for martensitic transformation above
the M
S
temperature can be compensated for by the application of stress.
Transformation of austenite to martensite is accomplished by a perma-
nent shape deformation. This is characterised by an invariant-plane strain
deformation with and s as the dilatational and shear components of the
shape deformation respectively. The mechanical driving force due to an ex-
127 6.2 Role of strain
ternal stress, , can then be expressed as:
G
mechanical
=
N
+s (6.1)
where
N
is the normal component and is the shear component of the
applied stress resolved along the direction of the shear displacement of the
shape deformation. Both
N
and act on the crystallographic plane of
austenite on which the plate of martensite forms (Figure 6.2).
Work done by the shear component of the applied stress is always ex-
pected to be positive whereas that due to the normal component may be pos-
itive (in tension) or negative (compression). In steels, the shear component
of the shape deformation is much higher in magnitude than the dilatational
component. Thus the observed eect of stress arises mainly from the work
done by the shear component. However, hydrostatic stress interacts only
with the dilatational component and therefore is not very eective for in-
ducing transformation in steels. Hydrostatic compression, in fact, suppresses
transformation since in steels, there is a volume expansion on martensitic
transformation.
6.2 Role of strain
As is apparent from Figure 6.1, the higher the temperature above M
S
, the
greater the stress required to initiate transformation. It is possible therefore
that above a critical temperature austenite will be plastically deformed before
transformation, which is often referred as strain-induced transformation.
Plastic strain in austenite creates defects that can act as additional
128 6.2 Role of strain
Figure 6.2: Resolved normal and shear components of the external
stress [Bhadeshia, 2001].
nucleation sites for transformation. This is in contrast to the stress-assisted
transformation that starts from the same nucleation sites as athermal trans-
formation but with the thermodynamic assistance of external stress. Plastic
strain thus renders a kinetic advantage over the thermodynamic contribu-
tion of the accompanying stress [Lecroisey and Pineau, 1972; Manganon and
Thomas, 1970; Olson and Cohen, 1972; Venables, 1964].
Transformation from plastically strained austenite has been generally
observed to occur at a stress level lower than that expected from the contri-
bution of stress alone. There is a dierence in opinion regarding the cause of
such an observation:
Olson and Cohen [1979] ascribed this behaviour to the higher potency
of the nucleation sites created by plastic strain.
Onodera and Tamura [1979] explained the same behaviour in terms of
the stress concentration at defects, which may reach the critical stress
required for transformation, even when the magnitude of the applied
stress is less.
129 6.3 Mathematical formulation
In spite of the nucleation advantage, plastic strain in austenite may
actually act to retard transformation. Defects generated due to strain can
hinder the growth of martensite plates. This is generally explained in terms
of the martensite-austenite interface structure, which is made up of a set of
screw dislocations, as depicted schematically in Figure 6.3 [Christian, 2002].
Growth of a martensite plate involves translation of the glissile inter-
face, which can be made sessile by obstacles like grain boundaries. Plates of
martensite are therefore never observed to traverse austenite grain bound-
aries. Smaller defects like isolated dislocations resist the motion of the inter-
face but can often be incorporated into the lattice of martensite. However,
heavy dislocation debris in austenite generated by severe plastic deforma-
tion can halt the transformation. This is known as mechanical stabilisa-
tion [Bhadeshia, 1999b; Breedis and Robertson, 1963; Leslie and Miller, 1964;
Machlin and Cohen, 1951; Raghavan, 1992].
6.3 Mathematical formulation
The purpose here is to developed a theory to predict the critical strain nec-
essary for causing mechanical stabilisation of austenite. The problem is one
of balancing the force required to sustain the motion of the transformation
interface against the resistance due to the dislocation debris in austenite.
Dislocation theory
At low homologous temperatures, plastic deformation is caused by the trans-
lation of dislocations. The resultant plastic strain is given by [Honeycombe,
130 6.3 Mathematical formulation
Figure 6.3: Model of austenite-martensite interface (after Frank [1953]), a
screw dislocation with a lattice Burgers vector is necessary at every sixth
atomic closed-packed plane to match the atomic rearrangements [Christian,
2002].
1984]:
= bL (6.2)
where is the dislocation density, b is the Burgers vector and L is the average
distance moved by the dislocations.
For a given dislocation density, the spacing between the dislocations
is given by l =
1
2
. Under these circumstances, the mean shear stress
required to move dislocations past each other is given as [Honeycombe,
1984]:
=
Gb
8(1 )l
=
Gb
1
2
8(1 )
(6.3)
131 6.3 Mathematical formulation
where G is the shear modulus and the Poissons ratio. On combining with
equation 6.2 and noting that the force per unit length is b,
b =
Gb
3
2
8(1 )
L
(6.4)
The mean free distance L must decrease as the plastic strain increases
[Barlat et al., 2002]:
L =
D
+D
(6.5)
where D is the original grain size of austenite prior to straining and is a
coecient about equal to 1 m [Barlat et al., 2002].
Besides dislocation strengthening, solid-solution hardening adds to the
total resistance to interface motion [Olson and Cohen, 1976a-c; Ghosh and
Olson, 2001]. Solid solution hardening coecients for various solutes in
austenite have been reported for tensile strength by Irvine et al. [1969].
These values can be converted into shear stress equivalents
S
using the
Tresca criterion.
Chemical driving force
The stress
T
driving the motion of the interface originates from the chemical
free energy change G of transformation [Olson and Cohen, 1976a-c]:
T
= G (6.6)
where is a constant assumed to be equal to unity and G = G
i.e.,
the magnitude of the driving force which is a function of temperature and
the chemical composition of the alloy.
132 6.4 Results and discussion
Lattice distortions due to the accommodation of the shape deforma-
tion accompanying transformation must be accounted for in terms of stored
energy. The magnitude of this is given by Christian [1979] to be about
600 J mol
1
. This has to be subtracted from the free energy change of trans-
formation (obtained using MTDATA), in order to determine the stress driv-
ing the interface motion. The magnitude of force per unit length available
for driving the austenite-martensite interface can be obtained by multiplying
T
with the Burgers vector b.
Model
Mechanical stabilisation occurs when the force driving the interface equals
that opposing its motion:
T
= +
S
(6.7)
bG =
1
8(1 )
Gb
3
2
L
+
S
b (6.8)
This equation can now be used to calculate the critical strain for mechanical
stabilisation.
6.4 Results and discussion
The key input to the model for calculating the critical strain necessary for
mechanical stabilisation is the free energy change accompanying the transfor-
mation of austenite to martensite, G. The free energy change is a function
of temperature and the composition of the steel and can be calculated using
the software MTDATA. The software accesses the thermodynamic database
SGTE to calculate phase stabilities and free energies. The free energy values
133 6.4 Results and discussion
obtained from MTDATA were reduced by 600 J mol
1
to account for the
stored energy of martensite. The magnitudes of dierent constants used in
calculations are listed in Table 6.1.
Property Value
Shear modulus of austenite 8 10
10
Pa
Poissons ratio 0.27
Burgers vector 2.52 10
10
m
Table 6.1: Values of constants used in the model.
6.4.1 Austenitic stainless steels
A range of chemical compositions of austenitic stainless steels are listed in
Table 6.2. The reason for varying the nickel concentration was to study the
eect of the chemical driving force for transformation on the stabilisation to
martensitic transformation. All the calculations were carried out allowing
for solid-solution strengthening and for a variety of temperatures.
Alloy C Mn Si Cr Ni Mo N
Alloy 1 0.05 1.5 0.75 18.0 6.0 0.0 0.03
Alloy 2 0.05 1.5 0.75 18.0 8.0 0.0 0.03
Alloy 3 0.05 1.5 0.75 18.0 10.0 0.0 0.03
Alloy 4 0.05 1.5 0.75 18.0 12.0 0.0 0.03
316L 0.01 0.53 0.75 17.1 11.9 2.0
Table 6.2: Chemical compositions (wt%) of dierent austenitic stainless
steels.
Figure 6.4 shows the chemical driving force (before the subtraction of
the 600 J mol
1
stored energy) as calculated using MTDATA. It is seen that
these are very large driving forces, especially for 300 K which is normally the
134 6.4 Results and discussion
ambient deformation temperature. The original grain size was taken to be
40 m.
The critical strain calculated using the model for each alloy at dierent
temperatures is presented in Figure 6.5. It can be seen that the strain values
are very large, ranging from =2.94.5. The strain values increase with
decreasing nickel concentration i.e., increasing the free energy change, G,
of transformation.
Figure 6.4: Free energy change in dierent alloys as listed in Table 6.2 at
dierent temperatures.
In physical terms, this means that the mean free path L has to be very
small, less than 1 m, as illustrated in Figure 6.6. Huge plastic strains are
required to achieve the necessary defect density, so such a result may at rst
135 6.4 Results and discussion
Figure 6.5: Calculated strain to initiate mechanical stabilisation at dierent
temperatures for various alloys.
sight seem unlikely. However, experiments carried out with 316L austenitic
stainless steel deformed to a true strain of = 6.3 conrms that severe
deformation is indeed necessary to mechanically stabilise the austenite [Wang
et al., 2005].
Figure 6.7 shows the curve obtained by tting to experimental data
from the 316L experiments [Wang et al., 2005]. It was demonstrated that
thermodynamically the alloy should be able to transform completely into
martensite, were it not for the fact that the deformed austenite becomes
mechanically stabilised. In fact, the amount of strain-induced martensite is
limited to 57%, with much of this formed at low strains before the onset
of stabilisation. On the basis of the tted curve (Figure 6.7), the onset of
mechanical stabilisation, i.e. the strain at which strain-induced martensite
ceases to signicantly increase with strain, is about = 2 3.
136 6.4 Results and discussion
Figure 6.6: Variation of the mean free path with strain for mechanical sta-
bilisation, open circles for 400 K and lled circles for 300 K.
Figure 6.7: Experimental tted curve for strain-induced martensite formation
in 316L stainless steel [Wang et al., 2005].
137 6.4 Results and discussion
The calculated data for 316L are presented as points in Figure 6.5;
they show that the critical strain for mechanical stabilisation is 3.45, 2.25
and 1.45 for the temperatures 300, 400 and 500 K respectively. The large
strain data of the experiments are based on wire drawing and it is conceivable
that the temperature rises by around 200 K during the deformation [Wang
et al., 2005]. The calculations are therefore consistent with the experimental
data. Furthermore, the analysis shows that large strains needed to achieve
stabilisation is due to the large values of G for transformation in such
stainless steels at the deformation temperature.
6.4.2 TRIP-assisted steels
In TRIP-assisted steels, high-carbon retained austenite transforms to
martensite during deformation. In order to calculate the strain required
to cause mechanical stabilisation in these steels, the chemical composition of
the retained austenite grains must be known precisely.
These steels are intercritically annealed in the (+) two phase region
prior to the bainite reaction. Retained austenite gets enriched with carbon
subsequent to the bainite reaction. The carbon content of retained austenite
is generally reported for these steels [Itami et al., 1995; Jacques et al., 2001c;
Sugimoto et al., 1992b; Tomita and Morioka, 1997]. However, substitutional
solutes do not partition during bainitic transformation. Therefore, the con-
centrations of these elements in the nal retained austenite remain the same
as that before the bainite reaction. The dwell time during intercritical anneal-
ing is very short and so equilibrium partitioning of the substitutional solutes
is very unlikely. By contrast, the carbon concentration can be expected to
138 6.4 Results and discussion
attain equilibrium in austenite within the short intercritical annealing period.
The substitutional solute content of the retained austenite can be calculated
using MTDATA, assuming the existence of equilibrium and paraequilibrium
during intercritical annealing separately.
The thermodynamic driving force for transformation of austenite to
martensite was calculated using MTDATA at the temperature of deforma-
tion. Retained austenite grains in these steels are generally reported to be
very ne (about 1 m) [Jacques, 2003]. These values were used as inputs to
the model to determine the critical strain necessary for mechanical stabilisa-
tion (Figure 6.8).
0
1
2
3
4
0.5 0.7 0.9 1.1 1.3 1.5
C
r
i
t
i
c
a
l
s
t
r
a
i
n
Carbon / wt%
Equilibrium
Paraequilibrium
Figure 6.8: Critical strain necessary for onset of mechanical stabilisation in
high-carbon austenite of TRIP-assisted steels.
As seen from Figure 6.8, the critical strain for onset of mechanical sta-
bilisation in these steels decreases with increasing carbon content of retained
139 6.4 Results and discussion
austenite. This is due to the lower chemical driving force for martensitic
transformation for higher carbon containing austenite. A range of critical
strain values are predicted for each carbon level: the upper and lower limits
arise from the existence of the equilibrium or paraequilibrium partitioning
of the substitututional solutes during intercritical annealing. However, the
range of predicted strain values is small in magnitude, possibly due to the
subtle dierence between the equilibrium and the paraequilibrium concen-
trations of the substitutional solutes in the austenite.
Figure 6.8 further demonstrates the occurrence of mechanical stabili-
sation at small plastic strain in austenite with carbon content higher than
about 1.1 wt%. It appears therefore that transformation of austenite with
such a carbon level can never be strain-induced. Interestingly, high-carbon
austenite is known to be the hardest phase in the microstructures of these
steels [Furnemont et al., 2002]. Austenite might thus be thought to remain
unstrained till the late stage of the deformation. It might be possible there-
fore that high-carbon austenite in these steels transforms to martensite due
to the eect of elastic stress rather plastic strain.
Stress-assisted transformation
Stress provides the shortfall in the chemical driving force for transformation
to occur. However, the magnitude of the assistance rendered by stress de-
pends on the nature of the stress state. For a simple uniaxial loading, the
stress coecient is reported by Olson [1982] to be 0.86 J mol
1
MPa
1
. With
this value, the total driving force (chemical + mechanical) can be calculated
as a function of temperature using MTDATA.
140 6.4 Results and discussion
Figure 6.9 illustrates the total thermodynamic driving force for trans-
formation at dierent stress levels for a hypothetical composition of austenite
as Fe-1C-2Mn-2Si, all in wt%. The gure reveals the remarkable increase in
the total driving force of transformation due to stress. The diagram also
helps to understand the eect of a given stress level in terms of an equivalent
thermal undercooling necessary to achieve the same additional thermody-
namic driving force. For example, application of a stress of 400 MPa may
be considered as an equivalent undercooling of 55 K, resulting in additional
transformation.
-3100
-2800
-2500
-2200
-1900
273 293 313 333 353 373
G
/
J
/
m
o
l
Temperature / K
0 MPa
400 MPa
600 MPa
Figure 6.9: Total driving force for austenite to martensite transformation for
Fe-1C-2Mn-2Si (wt%) austenite at dierent stress levels.
It can therefore be argued that the inuence of stress on the extent of
transformation can be computed in a similar way to athermal transformation,
141 6.4 Results and discussion
using the empirical equation given by Koistinen and Marburger [1959]:
1 = exp{C(M
S
T
Q
)} (6.9)
where is the martensite volume fraction at temperature T
Q
and C is a
constant equal to 0.011 K
1
. Knowing the M
S
temperature of austenite and
the equivalent thermal undercooling produced by a given stress level, the
fraction of austenite transforming to martensite can be calculated using this
equation.
Table 6.3 presents data, reported by Jacques et al. [2001d], that have
been used to calculate the inuence of stress on transformation of austenite,
using the principle described above.
Sample Athermal
R
C in
R
M
S
(vol.%) (vol.%) (wt%) (K)
L30s 7 8.8 0.61 351 (calculated)
L60s 5 8.1 0.68 341 (calculated)
H360 0 7.9 0.85 298 (stated)
Table 6.3: Data for decomposition of residual austenite to martensite during
cooling to room temperature,
R
is retained austenite [Jacques et al., 2001d].
The M
S
temperature for sample H360 has been indicated by Jacques
et al. [2001d] to be close to room temperature i.e. 298 K, as only a small
amount of martensite was detected in this sample. On the other hand, the
M
S
temperatures for samples L30s and L60s were calculated using Koistinen-
Marburger equation, based on the reported amount of transformation at
298 K.
Stress levels corresponding to dierent strain values, at which the re-
tained austenite volume fractions were reported, were worked out from the
142 6.4 Results and discussion
relevant true stress-strain curve. These values are presented in Table 6.4.
The table also shows the expected retained austenite fractions, as calculated
considering the thermodynamic eect of stress alone.
Sample Stress / MPa Calculated
R
vol.% Reported
R
vol.%
L30s 375 5.4 6.5
463 4.9 4.8
593 4.1 3.6
659 3.8 2.7
794 3.2 1.9
L60s 379 5.0 7.9
466 4.5 5.3
572 3.9 5.3
639 3.6 3.2
698 3.3 2.2
769 3.0 1.6
H360 500 4.0 6.0
710 3.0 4.8
830 2.6 3.6
914 2.3 3.4
Table 6.4: Calculated amounts of retained austenite (
R
) and the correspond-
ing amounts as reported by Jacques et al. [2001d] at each stress level.
Table 6.4 reveals the reported extent of transformation is higher than
that calculated for samples L30s and L60s. There is no contention of mechan-
ical stabilisation at the carbon levels of these samples (Figure 6.8). Plastic
strain in austenite creates defects that act as additional nucleation sites. This
can be considered to drive the transformation to a greater extent than that
expected from the eect of stress alone.
By contrast, transformation has been reported to a lower extent than
the computed ones for sample H360. Austenite in this sample contains about
0.85 wt% carbon. The high carbon level of austenite in H360 sample implies
143 6.4 Results and discussion
a high hardness level [Furnemont et al., 2002]. Thus austenite in this sample
might not have strained plastically. It is possible therefore that transforma-
tion has been due to the eect of stress alone.
6.4.3 Bainitic steels
One of the fascinating aspects of the bainite reaction in steels is the sub-unit
mechanism of its growth. A platelet of ferrite grows to a limited size even
though there is no impingement with obstacles such as austenite grain bound-
aries. The transformation then propagates by the nucleation and growth of
another sub-unit, the collection of sub-units being known as a sheaf [Hehe-
mann, 1970]. The reason why each platelet only grows to a limited size is that
the shape deformation accompanying transformation is plastically accommo-
dated in the austenite adjacent to the plate [Bhadeshia and Edmonds, 1979].
This results in the creation of an intense dislocation debris which renders the
interface immobile, as illustrated in Figure 6.10.
Calculations were carried out for a Fe-0.4C-3Mn-2Si (wt%) steel, which
is also the material from which the micrograph in Figure 6.10 originates with
a transformation temperature of 623 K. As pointed out previously, the stored
energy has to be subtracted from the free energy change of transformation.
For bainite, the magnitude of this is given by Bhadeshia [1981] as 400 J mol
1
.
The results of calculations are presented in Table 6.5.
The shear strain associated with the growth of bainite is estimated to
be about 0.26 [Swallow and Bhadeshia, 1996]. This, in comparison with the
calculated data, can be expected to be sucient for onset of local mechanical
stabilisation, and hence a sub-unit mechanism of growth. Furthermore, the
144 6.4 Results and discussion
Figure 6.10: Intense dislocation debris at, and in the vicinity of the bainite
(
ub
)-austenite () interface, the steel has a chemical composition Fe-0.43C-
3Mn-2.12Si wt% and is transformed at 623 K [Bhadeshia and Edmonds,
1979].
Temperature / K Critical Strain
500 2.1
600 0.6
650 0.06
Table 6.5: Strain necessary for mechanical stabilisation in the bainitic steel.
length of the sub-unit is expected to become larger and this is qualitatively
consistent with experimental data on bainite, the plates become more slender
as the transformation temperature is reduced (G increased) [Garcia-Mateo
et al., 2003]. It may be further noted that martensite forms in this steel
below about 500 K [Bhadeshia and Edmonds, 1979].
145 6.4 Results and discussion
6.4.4 Athermal martensite
Austenite on cooling starts transforming to martensite only below a critical
temperature, M
S
. This temperature varies with chemical composition of the
steel. Most of the common alloying elements, except aluminium and cobalt,
are known to depress M
S
. Finer austenite grain sizes also lead to a lower
M
S
temperature [Raghavan, 1992]. This is generally believed to be due to
the hardening of austenite due to ner grain sizes, hindering the growth of
martensite. Transformation from plastically deformed austenite in a simi-
lar way can be expected to start at lower temperature, as reported in a steel
of composition Fe-0.75C-1.02Mn-0.28Si-0.11Cr-0.05Ni-0.015S-0.009P-0.08Cu
(all in wt%) [Maalekian et al., 2006]. This is becasue of the strengthening
of austenite due to a higher dislocation density than that in the undeformed
condition. A methodology is described in this section to calculate the de-
pression of the M
S
temperature caused by prior plastic strain in austenite.
The theory for mechanical stabilisation relies on balancing the force
required to move the austenite-martensite interface through a forest of dislo-
cations against the force due to the chemical free energy change (G) driving
the interface:
Gb
1
2
8(1 )
= |G| (6.10)
where b = 0.252 nm is the magnitude of the Burgers vector, G = 8 10
10
Pa
is the shear modulus of the austenite, = 0.27 is its Poissons ratio and
is the dislocation density. The extra driving force, (G), necessary for
transformation in a plastically deformed sample can be expressed as:
Gb
8(1 )
[
1
2
1
2
0
] = (G) (6.11)
146 6.4 Results and discussion
where the dislocation density has a value
0
at zero plastic strain.
G for transformation of austenite to ferrite without a composition
change was calculated using MTDATA for the Fe-0.75C-1.02Mn-0.28Si-
0.11Cr-0.05Ni-0.015S-0.009P-0.08Cu (wt%) steel as a function of tempera-
ture:
G = 7.0337T 5218 J mol
1
or (G) 7.0337T (6.12)
where T is the absolute temperature. The equation applies over the temper-
ature range 450-510 K. The above expression therefore relates the magnitude
of the change in available driving force due to the presence of dislocations in
the austenite, (G), with the corresponding reduction in the M
S
temper-
ature, T.
It follows that the depression of the transformation temperature can
be calculated as a function of the change in the dislocation density of the
austenite. However, it would be useful to be able to estimate the change
in M
S
as a function of plastic strain rather than dislocation density. It
is necessary therefore to derive a relationship between plastic strain () and
dislocation density, which will necessarily be approximate. So two alternative
methods were attempted as follows:
Data from Vodopivec et al. [1999] are used to obtain an empirical rela-
tionship as:
= 2 10
13
0
+2 10
14
m
2
(6.13)
In the alternative approach, it is assumed that only 5% of the plastic
work is stored in the material, both due to changes in the austenite grain
147 6.4 Results and discussion
surface per unit volume, S
V
, and due to the expected change in dislocation
density. By balancing the plastic work stored against the energy of defects
created, we obtain:
0.05
y
S
V
+ (6.14)
where
y
is the yield stress of austenite (measured to be 30 MPa at
1000
+ 4e
2
]
1
2
+
1
3
[e
+ 2
]
1 + 2
3
(6.15)
The value of S
V
at zero deformation is given by S
V
0
= 2/d
where d
= 26
3 m is the mean lineal intercept dening the austenite grain size [Maalekian
et al., 2006].
The stored energy due to deformation is simply calculated using the
term on the left-hand side of equation 6.14. Removing the component due
to the change in S
V
(equation 6.15) from equation 6.14 gives the change in
dislocation density.
The predicted depression in the martensite-start temperature with in-
creasing prior plastic strain in austenite is shown in Figure 6.11. The calcu-
lated trends compare well with the measured value.
148 6.5 Summary
Figure 6.11: Calculated change in M
S
as a function of plastic strain in austen-
ite prior to its transformation, the points represent experimental data.
6.5 Summary
The role of stress and strain on martensitic transformation has been studied.
A theory has been developed to calculate the critical strain necessary for the
onset of mechanical stabilisation in steels. The chemical composition of the
alloys along with the transformation temperature are used as key inputs to
the model.
The theory has been used to demonstrate why large strain values are
necessary to stabilise certain stainless steels at ambient temperature. The
theory also helps to interpret the sub-unit mechanism of bainite growth and
the depression of M
S
temperature with prior plastic strain in austenite. The
theory indicates the possibility that high-carbon austenite in TRIP-assisted
steels might actually transform to martensite due to stress, rather strain,
during deformation.
Chapter 7
Formability
TRIP-assisted steels are expected to possess enhanced formability. This is the
ability of a material to be shaped, an ability related to plasticity. Formability
is characterised by dierent test methods. The purpose of the present chapter
is to study the inuence of simple mechanical properties such as the strength
and elongation on formability.
7.1 Stretch-angeability
Several methods have been designed to characterise the forming ability, one
of which is the hole expansion test. During this test, a hole is rst made
in the steel blank. The hole is then allowed to expand under a punch until
cracks appear on the surface. The ratio of the increase in diameter to the
initial diameter is called the hole expansion ratio. The higher the ratio, the
better is the stretch-angeability of the steel.
There are examples of strong steels tested for stretch-angeability [Mat-
sumura et al., 1992; Sudo et al., 1982]. Each of these materials possesses a
149
150 7.2 Neural networks
similar ferrite-rich microstructure with a dispersion of a hard phase which
could be martensite or bainite or both. The steels may also contain retained
austenite.
There has been much research in understanding the eect of microstruc-
ture on the stretch-angeability of these steels. For example, delayed strain-
induced transformation of the retained austenite into martensite is known to
improve the properties [Sugimoto et al., 1999]. Volume fraction, hardness and
distribution of the hard phase may also inuence the performance [Hasegawa
et al., 2004]. It has even been possible to achieve superior hole-expansion
ratio over the multiphase steels by using either fully bainitic or martensitic
steels [Sugimoto et al., 2002a,b, 2000].
However, there is no proper understanding on the inuence of the sim-
ple mechanical properties, such as strength and elongation, on the stretch-
angeability. During a tensile test, the specimen rst undergoes uniform
deformation. But plastic instability sets in when the work-hardening can
no longer compensate for the increased stress due to the reduction of cross
section. The resulting non-uniform strain is generally not considered to be
useful during forming, as this might degrade the aesthetic of the formed com-
ponent. To study the relationship between the tensile properties and the hole
expansion ratio forms the aim of the present work.
7.2 Neural networks
Neural network analysis was adopted in this study. This represents a general
method of regression that ts a non-linear function to experimental data.
Details of the method can be found in Chapter 4.
151 7.3 Results and discussion
Analysis
The data presented in Table 7.1 are due to Fang et al. [2003a,b]; Hasegawa et
al. [2004]; Sudo et al. [1982]; Sugimoto et al. [1999]. They originate from steels
with microstructures of ferrite-martensite or ferrite-bainite, with or without
retained austenite. The data were used to create models for predicting the
hole expansion ratio () against various combinations of the tensile properties
viz. yield strength (YS), ultimate tensile strength (UTS), uniform elongation
(UEL), total elongation (TEL), yield ratio (YR) and the strength-elongation
product (UTS-UEL and UTS-TEL).
Each variable in the database was normalised as
x
n
=
x x
min
x
max
x
min
0.5 (7.1)
where x
n
is the normalised value of the variable x, with x
min
and x
max
as the
minimum and the maximum values respectively. The database was divided
randomly into two halves. Each model was trained using one half of the
dataset. The other half was used to test the performance of the model with
unseen data.
7.3 Results and discussion
The rst model was created with the YS, UTS, UEL, YR and UTS-UEL
as the input parameters. Figure 7.1(a) illustrates the reduction in the noise
level with increasing number of hidden units. The nal model is formed
as a committee of models based on the minimum combined test error on
increasing the number of networks (Figure 7.1(b)). The nal committee in
152 7.3 Results and discussion
Inputs Minimum Maximum Average S.D.
YS / MPa 335 783 444 92
UTS / MPa 484 1028 675 155
UEL / % 8 39 18 7
TEL / % 10 43 27 8
YR 0.49 0.82 0.67 0.09
UTS-UEL / MPa-% 6992 29715 12192 6113
UTS-TEL / MPa-% 9016 32379 18069 6894
Output Minimum Maximum Average S.D.
/ % 9 151 65 45
Table 7.1: Data used in the analysis, S.D. represents standard deviation.
this case was made of only one sub-model.
Figure 7.2 compares the measured hole expansion ratios with the cor-
responding model predictions made over the entire dataset. In general, the
predicted hole expansion ratios compare well with the corresponding mea-
sured values, but there are a few outliers. However, the uncertainty is seen to
be large. The model-perceived signicances of the input variables inuencing
the hole expansion ratio are presented in Figure 7.3. As can be seen, UTS
gained the maximum importance followed by the YS and the YR. Contrary
to general belief, UEL and UTS-UEL were perceived to be insignicant.
The apparent irrelevance of the UEL and the UTS-UEL might be real.
This is because, the attainment of UTS in the ow curve implies the end of
UEL. Therefore, UTS and UEL are actually not independent variables. To
test this hypothesis, two separate models were created. One of these had the
UTS with the UEL as the input variables while the other had the YS with
the UEL. Both these models were found to consist of only one sub-model in
the nal committee.
153 7.3 Results and discussion
N
o
i
s
e
0.08
0.12
0.16
0.2
0.24
4 8 12 16 20
Number of hidden units
(a)
0.35
0.37
0.39
4 8 12 16 20
0.41
Number of models
C
o
m
b
i
n
e
d
T
e
s
t
E
r
r
o
r
(b)
Figure 7.1: Variation of (a) noise level, , with increasing number of hidden
units and (b) combined test error with number of models in committee.
154 7.3 Results and discussion
0
40
80
120
160
0 40 80 120 160
p
r
e
d
i
c
t
e
d
/
%
measured / %
flange_1_result.dat using 4:5
Figure 7.2: Predicted hole expansion ratios compared with the measured
values using the model with YS, UTS, UEL, YR and UTS-UEL as the input
parameters.
Figure 7.3: Signicances of YS, UTS, UEL, YR and UTS-UEL as the input
variables.
155 7.3 Results and discussion
(a)
(b)
Figure 7.4: Signicances of (a) UTS and UEL and (b) YS and UEL as the
input variables.
156 7.3 Results and discussion
Figure 7.4 shows the signicances of the variables as perceived by the
model. The signicance of UEL was indeed boosted in the presence of the
YS, as compared to that with the UTS. However, its magnitude remained
small (Figure 7.4(b)). This indicates a possibility that the TEL could be a
better parameter than the UEL in describing the hole expansion ratio. A
model was therefore created with YS, UTS, TEL, YR and UTS-TEL as the
input variables. The nal committee of the model was found to have eight
sub-models.
Figure 7.5 presents the signicances of the input variables as perceived
by each sub-model in the nal committee. UTS, UTS-TEL and TEL were
found to exert considerable inuence on the hole expansion ratio, while YS
and YR did not gain any importance. However, the signicances for each
input parameters vary in dierent sub-models.
Figure 7.5: Signicances of YS, UTS, TEL, YR and UTS-TEL as the input
variables.
157 7.3 Results and discussion
Figure 7.6 depicts the predictions made using this model over the en-
tire dataset. The predicted values compare well with the measured values.
Furthermore, the magnitudes of the error bars associated with the predic-
tions using this model are much reduced as compared to the rst model
(Figure 7.2). There are no signicant outliers in the predictions. All these
suggest, the model with the TEL as an input performs appreciably better
than the one with the UEL as an input.
0
40
80
120
160
0 40 80 120 160
p
r
e
d
i
c
t
e
d
/
%
measured / %
flange_2_result.dat using 4:5
Figure 7.6: Predicted hole expansion ratios compared with the measured
values using the model with YS, UTS, TEL, YR and UTS-TEL as the input
parameters.
The above results indicate the importance of UTS, TEL and UTS-
TEL in controlling the hole expansion ratio. This appears to be physically
reasonable, given that the hole is allowed to expand until cracks appear on the
surface. This therefore tests the ability of the material to deform plastically
158 7.3 Results and discussion
until fracture, albeit locally. Two dierent models were therefore created
with the input parameters as UTS and UTS-TEL in one case while UTS and
TEL in the other.
From Figure 7.7, the signicances of both UTS-TEL and TEL can be
seen to be appreciable, although quantitatively these are less than that of the
UTS. In each case, the model consisted of only one sub-model. Figure 7.8
shows that the predictions made over the entire dataset compare well with
the corresponding measured values.
These models were used further to predict the variation in the hole
expansion ratio as function of UTS, TEL and UTS-TEL. Figure 7.9 and Fig-
ure 7.10 present the results as surface plots with the corresponding uncertain-
ties which were between 10-25. At constant UTS-TEL, the hole expansion
ratio seems to be higher at lower UTS (Figure 7.9(a)). This might be due to
the higher TEL (keeping UTS-TEL the same). However, at the same UTS,
the ratio is seen to be decreasing with increasing UTS-TEL, ultimately be-
coming constant with UTS-TEL exceeding 20000 MPa-%. This is therefore
not consistent with the benecial eect of the total elongation.
Figure 7.10(a) illustrates a more drastic reduction in the hole expansion
ratio with increasing UTS as compared to Figure 7.9(a). At constant UTS,
the ratio is however seen to remain almost constant with increasing total
elongation, although a positive correlation between the stretch-angeability
and TEL could be imagined ideally.
159 7.4 Summary
7.4 Summary
The present work identies the ultimate tensile strength to be the most
important tensile property inuencing the hole expansion ratio. The uniform
elongation loses its signicance as an input to the model, which also includes
the ultimate tensile strength as a variable. This is because both the properties
appear to be identical with respect to the ow curve.
The stretch-angeability deteriorates with increasing ultimate tensile
strength. Total elongation might be thought to be important, as hole ex-
pansion tests rely on the cracking of the surface while expanding the hole.
However, the present work suggests that increasing total elongation, keeping
the ultimate tensile strength same, might not be benecial.
160 7.4 Summary
(a)
(b)
Figure 7.7: Signicances of (a) UTS and UTS-TEL and (b) UTS and TEL
as the input variables.
161 7.4 Summary
0
40
80
120
160
0 40 80 120 160
p
r
e
d
i
c
t
e
d
/
%
measured / %
flange_6_result.dat using 4:5
(a)
0
40
80
120
160
0 40 80 120 160
p
r
e
d
i
c
t
e
d
/
%
measured / %
flange_7_result.dat using 4:5
(b)
Figure 7.8: Predicted hole expansion ratios compared with the measured
values using the models with (a) UTS and UTS-TEL and (b) UTS and TEL
as the input parameters.
162 7.4 Summary
600
800
1000
UTS / MPa
10000
20000
30000
UTS-TEL / MPa-%
40
80
120
160
predicted / %
(a)
600
800
1000
UTS / MPa
10000
20000
30000
UTS-TEL / MPa-%
10
15
20
uncertainty / %
(b)
Figure 7.9: Predicted variation in (a) the hole expansion ratio with UTS and
UTS-TEL and (b) the corresponding uncertainties.
163 7.4 Summary
600
800
1000
UTS / MPa
10
20
30
40
TEL / %
40
80
120
160
predicted / %
(a)
600
800
1000
UTS / MPa
10
20
30
40
TEL / -%
10
15
20
25
uncertainty / %
(b)
Figure 7.10: Predicted variation in (a) the hole expansion ratio with UTS
and TEL and (b) the corresponding uncertainties.
Chapter 8
Conclusions
TRIP-assisted steels have the advantage of large uniform elongation in spite
of the high strength. This makes them suitable for making cars with lighter
body weight but adequate safety. The steels contain some 10-30 vol.% high-
carbon retained austenite in ferrite-rich microstructures. The transformation
of austenite into martensite during deformation is generally believed to be
responsible for the properties of these steels. This is despite the fact that
the martensite that forms in these steels must inherit the carbon of austenite
and hence should be brittle.
It has been demonstrated here that the ner martensite plates are bet-
ter able to resist cracking. In TRIP-assisted steels, retained austenite exists
in a nely divided state with the grain size of around 1 m. The martensite
that forms from this austenite therefore remains with a ne plate size and
hence is not a liability for these steels, in spite of its high carbon content.
To maximise the properties, the stability of the retained austenite
should be such that it transforms progressively during deformation. Car-
bon content of the austenite inuences its stability. Retained austenite vol-
164
165
ume fraction and its carbon content are therefore important microstructural
parameters inuencing the properties of these steels. Using published data,
neural network models were created relating these parameters with the chem-
ical composition and the heat treatment conditions of these materials. These
models were found to be capable of illustrating the eects of various factors
on the retained austenite fraction and its carbon content, as expected from
the metallurgical principles.
The retained austenite content in steels transformed to bainite increases
with the addition of silicon. However, silicon leads to complications during
processing and deteriorates the surface condition of the steel. The neural
network model, along with a genetic algorithm, was used to formulate an al-
loy, lean in silicon concentration yet capable of retaining austenite in copious
amounts. The computational methods gave rise to a new type of material,
here designated -TRIP steel. The alloy shows an enhanced tendency for
-ferrite formation during solidication and cannot be made fully austenitic.
Usually the ferrite, that constitutes the microstructure of TRIP-assisted steel,
is generated out of the allotriomorphic transformation of austenite. By con-
trast, the present alloy gives an opportunity to substitute some of the allotri-
omorphic ferrite in these materials by -ferrite.
Experiments using alloys based on the optimisation calculation demon-
strated the formation of -ferrite during solidication. These alloys were
further heat-treated producing radically dierent microstructures containing
-ferrite, bainite and retained austenite. Specimens in heat-treated condi-
tions recorded impressive mechanical properties at room temperature. How-
ever, the ductility deteriorates when the tests are done at 100
C. This could
166
be due to the enhanced stability of the retained austenite against transfor-
mation into martensite at the high temperature.
The transformation of the retained austenite into martensite during
deformation could be either due to stress or strain. Stress generally adds to
the driving force while plastic strain retards the growth of martensite plates,
eventually causing the mechanical stabilisation of the austenite. A theory
capable of predicting the critical strain necessary for the onset of stabilisation
has been developed. Using the theory, it has been possible to explain the
requirement of the large strain for stabilising certain austenitic steels and
the sub-unit mechanism of bainite growth. Calculations involving the theory
indicate that the high-carbon retained austenite in TRIP-assisted steels may
actually transform due to stress rather strain. Further experimental work
is necessary to clarify the role of stress and strain on the transformation of
retained austenite in these steels.
There are various complex tests to assess the formability of these
materials. However, there has been no proper attempt in understanding
these results in terms of the tensile properties. Published data on stretch-
angeability of these steels were analysed using neural networks. It has been
shown that the ultimate tensile strength of the steel is the most important
criterion inuencing the stretch-angeability.
Bibliography
Andrews, K. W. Journal of The Iron and Steel Institute 203 (1965) 721
Angel, T. Journal of The Iron and Steel Institute 177 (1954) 165
Barlat, F., Glazov, M. V., Brem, J. C. and Lege, D. J. International Journal
of Plasticity 18 (2002) 919
Basuki, A. and Aernoudt, E. Journal of Materials Processing Technology
89-90 (1999) 37
Bhadeshia, H. K. D. H. and Honeycombe, R. W. K. Steels: Microstructure
and Properties, Elsevier, UK, (2006)
Bhadeshia, H. K. D. H. ISIJ International 42 (2002) 1059
Bhadeshia, H. K. D. H. Bainite in Steels, IOM Communications Ltd, UK
(2001)
Bhadeshia, H. K. D. H. ISIJ International 39 (1999) 966
Bhadeshia, H. K. D. H. Materials Science and Engineering A 273-275 (1999)
58
Bhadeshia, H. K. D. H. Acta Metallurgica 29 (1981) 1117
167
168 BIBLIOGRAPHY
Bhadeshia, H. K. D. H. and Edmonds, E. V. Metallurgical and Materials
Transactions A 10 (1979) 895
Bouet, M., Fillipone, R., Essadiqi, E., Root, J. and Yue, S. 42nd Mechani-
cal Working and Steel Processing Conference Proceedings, Iron and Steel
Society/AIME, USA, 38 (2000) 55
Bouet, M., Root, J., Es-Sadiqui, E. and Yue, S. 40th Mechanical Working and
Steel Processing Conference Proceedings, Iron and Steel Society/AIME,
USA, 36 (1998) 675
Bouet, M., Root, J., Es-Sadiqi, E. and Yue, S. Materials Science Forum
284-286 (1998) 319
Breedis, J. F. and Robertson, W. D. Acta Metallurgica 11 (1963) 547
Capdevila, C., Caballero, F. G. and Garcia de Andres, C. Journal de
Physique IV 112 (2003) 217
Chen, H. C., Era, H. and Shimizu, M., Metallurgical and Materials Transac-
tions A 20 (1989) 437
Choi, I. D., Bruce, D. M., Kim, S. J., Lee, C. G., Park, S. H., Matlock, D.
K. and Speer, J. G. ISIJ International 42 (2002) 1483
Choi, B. Y., Krauss, G and Matlock, D. K. Scripta Metallurgica 22 (1988)
1575
Christian, J. W. The Theory of Transformations in Metals and Alloys, Perg-
amon, UK (2002)
Christian, J. W. Proceedings of ICOMAT79, MIT Press, USA (1979) 230
169 BIBLIOGRAPHY
Cox, H. L. British Journal of Applied Physics 3 (1952) 72
Davies, R. G. and Magee, C. L. Metallurgical and Materials Transactions A
3 (1972) 307
De Meyer, M., Vanderschueren, D., De Cooman, B. C., ISIJ International
39 (1999) 813
De Meyer, M., Vanderschueren, D., De Cooman, B. C., 41st Mechanical
Working and Steel Processing Conference Proceedings, Iron and Steel So-
ciety/AIME, USA, (1999) 265
Dieter, G. E., Mechanical Metallurgy, McGraw-Hill, UK (1988)
Dyson, D. J. and Holmes, B Journal of The Iron and Steel Institute 208
(1979) 469
Epogy 2003 Users guide, Synaps Inc
Fang, X., Fan, Z., Ralph, B., Evans, P. and Underhill, R. Journal of Materials
Processing Technology 132 (2003) 215
Fang, X., Fan, Z., Ralph, B., Evans, P. and Underhill, R. Journal of Materials
Science 38 (2003) 3877
Fischer, J. C. Metals Transactions 185 (1949) 688
Frank, F. C. Acta Metallurgica 1 (1953) 15
Furnemont, Q., Kempf, M., Jacques, P. J., Goken, M. and Delannay, F.
Materials Science and Engineering A 328 (2002) 26
170 BIBLIOGRAPHY
Garcia-Mateo, C., Caballero, F. G. and Bhadeshia, H. K. D. H. ISIJ Inter-
national 43 (2003) 1238
Gallagher, M. F., Speer, J. G. and Matlock, D. K. Symposium on the Ther-
modynamics, Kinetics, Characterization and Modeling of: Austenite For-
mation and Decomposition, Minerals, Metals and Materials Society (TMS),
USA, (2003) 563
Gerberich, W. W., Thomas, G., Parker, E. R. and Zackay, V. F. Proceedings
of 2nd International Conference on Strength of Metals and Alloys, USA,
3 (1970)
Ghosh, G., and Olson, G. B. Journal of Phase Equilibrium 22 (2001) 199
Girault, E., Martens, A., Jacques, P., Houbaert, Y., Verlinden, B. and Van
Humbeeck, J. Scripta Materialia 44 (2001) 885
Goel, N. C., Sangal, S. and Tangri, K. Metallurgical and Materials Transac-
tions A 16 (1985) 2013
Guimaraes, J. R. C. and Saavedra, A. Materials Science and Engineering A
62 (1984) 11
Guimaraes, J. R. C. and Gomes, J. C. Acta Metallurgica 26 (1978) 1591
Guimaraes, J. R. C. Scripta Metallurgica 6 (1972) 795
Hanzaki, A. Z., Hodgson, P. D. and Yue, S. ISIJ International 35 (1995) 324
Hanzaki, A. Z., PhD Thesis, McGill University, (1994)
Hasegawa, K., Kawamura, K., Urabe, T. and Hosoya, Y. ISIJ International
44 (2004) 603
171 BIBLIOGRAPHY
Hashimoto, S., Ikeda, S., Sugimoto, K-I. and Miyake, S. ISIJ International
44 (2004) 1590
Hassani, F. and Yue, S. 41st Mechanical Working and Steel Processing Con-
ference Proceedings, Iron and Steel Society/AIME, USA, 37 (1999) 493
Hehemann, R. F. Phase Transformations, ASM, USA (1970) 397
Hiwatashi, S., Takahashi, M., Katayama, T. and Usuda, M., Journal of the
Japan Society for Technology of Plasticity 35 (1994) 1109
Honeycombe, R. W. K. The Plastic Deformation of Metals, Second Edition,
Edward Arnold, London, (1984)
Im, D. B., Lee, C. G., Kim, S. J. and Park, I. M. Journal of the Korean
Institute of Metals and Materials 40 (2002) 8
Im, D. B., Lee, C. G., Kim, S. J. and Park, I. M. Journal of the Korean
Institute of Metals and Materials 38 (2000) 447
Im, D. B., Lee, C. G., Song, B. H., Kim, S. J. and Park, I. M. Journal of the
Korean Institute of Metals and Materials 38 (2000) 1328
Imai, N., Komatsubara, N. and Kunishige, K. Japan Technical Information
Service (1992) 25
Irvine, K. J., Gladman, T. and Pickering, F. B. Journal of The Iron and
Steel Institute 207 (1969) 1017
Itami, A., Takahashi, M. and Ushioda, K. ISIJ International 35 (1995) 1121
Jacques, P. Journal de Physique IV 112 (2003) 297
172 BIBLIOGRAPHY
Jacques, P. J., Girault, E., Martens, A., Verlinden, B., Van Humbeeck, J,
and Delannay, F. ISIJ International 41 (2001) 1068
Jacques, P. J., Girault, E., Harlet, P. and Delannay, F. ISIJ International
41 (2001) 1061
Jacques, P., Furnemont, Q., Mertens, A. and Delannay, F. Philosophical
Magazine A 81 (2001) 1789
Jacques, P., Ladriere, J. and Delannay, F. Metallurgical and Materials Trans-
actions A 32 (2001) 2759
Jacques, P., Girault, E., Catlin, T., Geerlofs, N., Kop, T., Van der Zwaag,
S. and Delannay, F. Materials Science and Engineering A 273-275 (1999)
475
Jiao, S., Hassani, F., Donaberger, R. L., Essadiqi, E. and Yue, S. ISIJ Inter-
national 42 (2002) 299
Kim, S. J., Lee, C. G., Lee, T. H. and Oh, C. S. ISIJ International 42 (2002)
1452
Kim, S. J., Lee, C. G., Choi, I. and Lee, S. Metallurgical and Materials
Transactions A 32 (2001) 505
Kim, C. Metallurgical and Materials Transactions A 19 (1988) 1263
Kim, S. J. and Lee, C. G. Journal of the Korean Institute of Metals and
Materials 37 (1999) 774
Koh, H. J., Lee, S. K., Park, S. H., Choi, S. J., Kwon, S. J. and Kim, N. J.
Scripta Materialia 38 (1998) 763
173 BIBLIOGRAPHY
Koistinen, P. P. and Marburger, R. E. Acta Metallurgica 7 (1959) 59
Konieczny, A. A. Processing and Fabrication of Advanced Materials XI, ASM
International, USA, (2003) 345
Lecroisey, F. and Pineau, A. Metallurgical and Materials Transactions A 3
(1972) 387
Lee, C. G., Kim, S. J., Lee, T. H. and Lee, S. Material Science and Engi-
neering A 371 (2004) 16
Lee, C. G., Kim, S. J., Song, B. H. and Lee, S. Metals and Materials 8 (2002)
435
Lee, C. G., Kim, S. J., OH, C. S. and Lee, S. ISIJ International 42 (2002)
1162
Leslie, W. C. and Miller, R. L. ASM Transactions Quaterly 57 (1964) 972
Ludwigson, D. C. and Berger, J. A., Journal of The Iron and Steel Institute
207 (1969) 63
Maalekian, M., Lendinez, M. L., Kozeschnik, E., Brantner, H. P., Cerjak,
H., Eect of hot plastic deformation of austenite on the transformation
characteristics of eutectoid carbon steel under fast heating and cooling
conditions, (2006, unpublished work)
Machlin, E. S. and Cohen, M. Trans. AIME 191 (1951) 746
MacKay, D. J. C., Information Theory, Inference and Learning Algorithms,
Cambridge University Press, UK, (2003)
174 BIBLIOGRAPHY
Mahieu, J., De Cooman, B. C., Maki, J., and Claessens, S. Iron and Steel-
maker 29 (2002) 29
Mahieu, J., Claessens, S. and De Cooman, B. C. Verlag Stahleisen GmbH,
Galvatech 2001, Germany, (2001) 644
Manganon, P. L. and Thomas, G. Metallurgical and Materials Transactions
A 1 (1970) 1577
Manohar, P. A., Kunishige, K., Chandra, T. and Ferry, M. Materials Science
and Technology 18 (2002) 856
Manohar, P. A., Kunishige, K. and Chandra, T., Materials Science Forum
426-432 (2003) 1127
MAP STEEL MUCG73 available freely in the world wide web,
http://www.msm.cam.ac.uk/map/steel/programs/mucg73-b.html
Maki, J., Mahieu, J., De Cooman, B. C. and Claessens, S. Journal of Mate-
rials Science and Technology 19(2003) 125
Maki, J., Mahieu, J., Claessens, S. and De Cooman, B. C. Verlag Stahleisen
GmbH, Galvatech 2001, Germany, (2001) 623
Mark, C. Proceedings of Cambridge Philosophical Society 52 (1956) 216
Matas, S. J. and Hehemann, R. F. Trans. Met. Soc. AIME 221 (1961) 179
Matsumura, O., Sakuma, Y. and Takechi, H. ISIJ International 32 (1992)
1014
Matsumura, O., Sakuma, Y. and Takechi, H. Transactions of the Iron and
Steel Institute 27 (1987) 570
175 BIBLIOGRAPHY
Matsumura, O., Sakuma, Y. and Takechi, H., Scripta Metallurgica 27 (1987)
1301
Mintz, B. Materials Science Forum 426-432 (2003) 1219
Mintz, B. International Materials Reviews 46 (2001) 169
MTDATA, National Physical Laboratory, Teddington, London, 2005
Murr, L. E. Interfacial phenomena in metals and alloys, Addison-Wesley
Publication Co., USA, (1975) 131
Murugananth, M., Babu, S. S. and David, S. A. Welding Journal, Supple-
ment, October, (2004) 267
Nagasaka, A., Sugimoto, K. and Kobayashi, M. ASM International, USA,
(1996) 557
Nakagaito, T., Shimizu, T., Furukimi, O. and Sakata, K., Tetsu-to-Hagane
89 (2003) 841
Nagasaka, A., Sugimoto, K-I., Kobayashi, M., Kobayashi, Y., Hashimoto, S.
Tetsu-to-Hagane 85 (1999) 885
Nagasaka, A., Sugimoto, K-I., Kobayashi, M. and Hashimoto, S. Tetsu-to-
Hagane 85 (1999) 552
Nagasaka, A., Sugimoto, K-I., Kobayashi, M., Shirasawa, H. Tetsu-to-Hagane
84 (1998) 218
Nelder J. A. and Mead R. Computer Journal 7 (1965) 308
176 BIBLIOGRAPHY
Olson, G. B. Deformation, Processing and Structure, American Society of
Metals, USA, (1982) 391
Olson, G. B. and Cohen, M. Proceedings of US-Japan Seminar on Mechan-
ical Behaviour of Metals and Alloys Associated with Displacive Transfor-
mation, Troy, NY, (1979) 7
Olson, G. B. and Cohen, M., Metallurgical and Materials Transactions A 7
(1976) 1897
Olson, G. B. and Cohen, M. Metallurgical and Materials Transactions A 7
(1976) 1905
Olson, G. B. and Cohen, M., Metallurgical and Materials Transactions A 7
(1976) 1915
Olson, G. B. and Cohen, M. Metallurgical and Materials Transactions A 6
(1975) 791
Olson, G. B. and Cohen, M. Journal of the Less-common Metals 28 (1972)
107
Onodera, H. and Tamura, I. Proceedings of US-Japan Seminar on Mechan-
ical Behaviour of Metals and Alloys Associated with Displacive Transfor-
mation, Troy, NY (1979) 24
Patel, J. R. and Cohen, M. Acta Metallurgica 1 (1953) 531
Pereloma, E. V., Timokhina, I. B. and Hodgson, P. D. Materials Science and
Engineering A 273-275 (1999) 448
177 BIBLIOGRAPHY
Pichler, A., Traint, S., Arnolder, G., Stiaszny, P., Blaimschein, M., and
Werner, E. A. Iron and Steelmaker, USA, 30, (2003) 21
Pichler, A., Traint, S., Arnoldner, G., Stiaszny, P., Blaimschein, M. and
Werner, E. A. 44th Mechanical Working and Steel Processing Conference
Proceedings, Iron and Steel Society/AIME, USA, 40 (2002) 121
Pichler, A., Stiaszny, P., Potzinger, R., Tikal, R. and Werner, E. 40th Me-
chanical Working and Steel Processing Conference Proceedings, Iron and
Steel Society/AIME, USA 36 (1998) 259
Pickering, F. B. Physical Metallurgy and the Design of Steels, Applied Science
Publishers Ltd, England, (1978)
Pyshmintsev, I. Y., De Meyer, M., De Cooman, B. C., Savray, R. A.,
Shveykin, V. P.and Vermeulen, M., Metallurgical and Materials Trans-
actions A 33 (2002) 1659
Raghavan, V. Martensite, A tribute to Morris Cohen, ASM International,
USA, (1992) 197
Rao, B. V. N. and Rashid, M. S. Metallography 16 (1983) 19
Rashid, M. S. Formable HSLA and Dual Phase Steels, USA, (1979) 1
Rashid, M. S. Dual Phase and Cold Pressing Vanadium Steels in the Auto-
mobile Industry, Berlin, (1979) 32
Reisner, G., Warner, E. A. and Fischer, F. D. International Journal of Solid
Structures 35 (1998) 2457
178 BIBLIOGRAPHY
Rigsbee, J. M. Proceedings of the ICOMAT79 Cambridge, MA, USA (1979)
381
Sadhukhan, S., Das, K. P., Bandyopadhyay, N. R. and Banerjee, M. K.
Journal of the Institution of Engineers, India, 82 (2001) 65
Sakuma, Y., Matlock, D. K., Krauss, G. Materials Science and Technology
9 (1993) 718
Sakuma, Y., Matlock, D.K., Krauss, G., Metallurgical and Materials Trans-
actions A 23 (1992) 1221
Sakuma, Y., Matsumara, O. and Takechi, H., Metallurgical and Materials
Transactions A 22 (1991) 489
Sakuma, Y., Matsumura, O. and Akisue, O. ISIJ International 31 (1991)
1348
Sangal, S., Goel, N. C. and Tangri, K. Metallurgical and Materials Transac-
tions A 16 (1985) 2023
Scheil, E., Zeitschrift f ur anorganische and allgemeine, Chemie 207 (1932)
21
Schrader, A. and Wever, F. Arch. Eisenhittenwesen 23 (1952) 489
Sherif, M. Y., Garcia-Mateo, C., Sourmail, T. and Bhadeshia, H. K. D. H.
Materials Science and Technology 40 (2004) 319
Shi, W., Li, L., Zhou, Y., Fu, R. Y. and Wei, X. C. Heat Treatment of Metals,
China, 8 (2002) 9
179 BIBLIOGRAPHY
Song, S. M., Sugimoto, K-I., Kandaka, S., Futamura, A., Kobayashi, M. and
Masuda, S. Materials Science Research International, Japan, 9 (2003) 223
Song, S. M., Sugimoto, K-I., Kandaka, S., Futamura, A., Kobayashi, M. and
Masuda, S. Journal of the Society of Materials Science, Japan, 50 (2001)
1091
Streicher, A. M., Speer, J. G. and Matlock, D. K. Steel Research 73 (2002)
287
Sudo, M., Tsukatani, I. and Shibata, Z. Metallurgy of Continuous Annealed
Sheet Steel, AIME (1982) 310
Sugimoto, K-I., Tsunezawa, M., Hojo, T. and Ikeda, S. ISIJ International
44 (2004) 1608
Sugimoto, K-I., Nakano, K., Song, S-M. and Kashima, T. ISIJ International
42 (2002) 450
Sugimoto, K-I., Kanda, A., Kikuchi, R., Hashimoto, S. I., Kashima. T. and
Ikeda, S. ISIJ International 42 (2002) 910-915
Sugimoto, K-I., Sakaguchi, J., Iida, T. and Kashima, T. ISIJ International
40 (2000) 920
Sugimoto, K-I., Nagasaka, A., Kobayashi, M. and Hashimoto, S. I. ISIJ
International 39 (1999) 56
Sugimoto, K-I., Kobayashi, M. and Nagasaka, A. Canadian Institute of Min-
ing, Metallurgy and Petroleum (Canada) (1996) 503
180 BIBLIOGRAPHY
Sugimoto, K-I., Kobayashi, M. Nagasaka, A. and Hashimoto, S. I. ISIJ In-
ternational 35 (1995) 1407
Sugimoto, K-I. and Kobayashi, M., Iron and Steel Society/AIME, USA
(1994) 255
Sugimoto, K-I., Misu, M., Kobayashi, M. and Shirasawa, H. ISIJ Interna-
tional 33 (1993) 775
Sugimoto, K-I., Usui, N., Kobayashi, M. and Hashimoto, S. I. ISIJ Interna-
tional 32 (1992) 1311
Sugimoto, K-I., Kobayashi, M. and Hashimoto, S. I., Metallurgical and Ma-
terials Transactions A 23 (1992) 3085
Swallow, E. and Bhadeshia, H. K. D. H. Materials Science and Technology
12 (1996) 121
Takahashi, M. Nippon Steel Technical Report 88 (2003) 2-7
Takechi, H., Matsumura, O., Sakuma, Y. Japan Kokai Tokyo Koho Japan
Patent 62, 188, 729 (1987)
Tamura, I., Maki, T., Hato, H., Tomota, Y. and Okada, M. Proceedings of
2nd International Conference on Strength of Metals and Alloys, Asilomar,
CA, 3 (1970) 894
Tomita, Y. and Morioka, K. Materials Characterization 38 (1997) 243
Traint, S., Pichler, A., Stiaszny, P. and Wemer, E. A. 44th Mechanical
Working and Steel Processing Conference Proceedings, Iron and Steel So-
ciety/AIME, USA 40 (2002) 139
181 BIBLIOGRAPHY
Traint, S., Pichler, A., Tikal, R., Stiaszny, P., Wemer, E. A., 42nd Mechan-
ical Working and Steel Processing Conference Proceedings, Iron and Steel
Society/AIME, USA 38 (2000) 549
Uenishi, A., Kuriyama, Y. and Takahashi, M. Nippon Steel Technical Report
(Japan) 81 (2000) 17
ULSAB-AVC, http://www.worldautosteel.org/ulsab-avc
Venables, J. A. Philosophical Magazine 7 (1964) 35
Vodopivec, F., Reskovi, S. and Mamuzic, I. Materials Science and Technology
15 (1999) 1293
Wang, H. S., Yang, J. R., Bhadeshia, H. K. D. H. Materials Science and
Technology 21 (2005) 1323
Wang, Z. C., Kim, S. J., Lee, C. G. and Lee, T. H., Journal of Materials
Processing Technology 151 (2004) 141
Wei, X. C., Li, L., Fu, R. Y. and Shi, W. Journal of Iron and Steel Research
International, China, 10 (2003) 49
Wei, X., Fu, R., Li, L. and Shi, W., Shanghai Jinshu (Shanghai Metals),
China, 24 (2002) 32
Wei, X. C., Li, L., Fu, R. Y. and Shi. W. Acta Metallurgica Sinica (English
Letters), China 15 (2002) 285
Wiles, D. B. and Young, R. A. Journal of Applied Crystallography 14 (1981)
149
182 BIBLIOGRAPHY
Yan, B. 44th Mechanical Working and Steel Processing Conference Proceed-
ings, Iron and Steel Society, AIME, USA, 40 (2002) 509
Yan, B. and Xu, K. Iron and Steelmaker 30 (2003) 33
Yan, B. and Xu, K. 44th Mechanical Working and Steel Processing Confer-
ence Proceedings, Iron and Steel Society, AIME, USA 40 (2002) 493
Zackay, V. F., Parker, E. R., Fahr, D. and Bush, R., Trans. ASM 60 (1967)
25
Zaeerer, S., Ohlert, J. and Bleck, W. Acta Materialia 52 (2004) 2765
Zhu, C. M. S. Q., Sellars, C. M. and Bhadesiha, H. K. D. H. Materials Science
and Technology (2006) submitted