Temperature-dependent elastic anisotropy and mesoscale deformation in a nanostructured ferritic alloy

Abstract

Nanostructured ferritic alloys are a new class of ultrafine-grained oxide dispersion-strengthened steels that have promising properties for service in extreme environments in future nuclear reactors. This is due to the remarkable stability of their complex microstructures containing numerous Y-Ti-O nanoclusters within grains and along grain boundaries. Although nanoclusters account primarily for the exceptional resistance to irradiation damage and high-temperature creep, little is known about the mechanical roles of the polycrystalline grains that constitute the ferritic matrix. Here we report an in situ mesoscale characterization of anisotropic responses of ultrafine ferrite grains to stresses using state-of-the-art neutron diffraction. We show the experimental determination of single-crystal elastic constants for a 14YWT alloy, and reveal a strong temperature-dependent elastic anisotropy that leads to elastic softening and instability of the ferrite. We also demonstrate, from anisotropy-induced intergranular strains, that a deformation crossover exists from low-temperature lattice hardening to high-temperature lattice softening in response to extensive plastic deformation.

ARTICLE

Received 27 Feb 2014 | Accepted 8 Sep 2014 | Published 10 Oct 2014

Nanostructured ferritic alloys are a new class of ultrane-grained oxide dispersion-strengthened steels that have promising properties for service in extreme environments in future nuclear reactors. This is due to the remarkable stability of their complex micro-structures containing numerous Y-Ti-O nanoclusters within grains and along grain boundaries. Although nanoclusters account primarily for the exceptional resistance to irradiation damage and high-temperature creep, little is known about the mechanical roles of the polycrystalline grains that constitute the ferritic matrix. Here we report an in situ mesoscale characterization of anisotropic responses of ultrane ferrite grains to stresses using state-ofthe-art neutron diffraction. We show the experimental determination of single-crystal elastic constants for a 14YWT alloy, and reveal a strong temperature-dependent elastic anisotropy that leads to elastic softening and instability of the ferrite. We also demonstrate, from anisotropy-induced intergranular strains, that a deformation crossover exists from low-temperature lattice hardening to high-temperature lattice softening in response to extensive plastic deformation.

DOI: 10.1038/ncomms6178

Temperature-dependent elastic anisotropy and mesoscale deformation in a nanostructured ferritic alloy

G.M. Stoica1, A.D. Stoica1, M.K. Miller2 & D. Ma1

1 Chemical and Engineering Materials Division, Spallation Neutron Source, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6139, USA. 2 Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge , Tennessee 37831-6139, USA. Correspondence and requests for materials should be addressed to D.M. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6178

For high-temperature applications, polycrystalline grains can be detrimental in structural materials, such as Ni-based superalloys, but they are benecial for creep-resistant

nanostructured ferritic alloys (NFAs)13. The ultrane grains in NFAs remain unusually stable with little or no grain growth, coarsening or recrystallization up to 800900 C. Although nanoclusters (NCs) and precipitates that are embedded in the grains and along the grain boundaries4,5, together with solute segregation of W and Cr to the grain boundaries, are widely recognized to account for the excellent creep resistance2,4,68, the ultrane ferrite grains are also essential in promoting the mechanical strength, for example, via Hall-Petch hardening9. Thus, it is important to understand how the ferrite grains deform in response to stresses, particularly, at elevated temperatures.

When subjected to an applied load, each grain in a polycrystal-line NFA behaves like a single crystal. It is known that a single crystal exhibits elastic and plastic anisotropy, as its elastic stiffness and plastic strength depend on the crystallographic orientation. The mechanical response of a NFA grain, within a polycrystalline NFA aggregate, is dependent on its single-crystal properties, the crystallographic orientation of its lattice and those of its neighbouring grains10. In particular, iron and its alloys are among the most anisotropic body-centred cubic metals, for instance, o1114 is the stiffest orientation that is B2.2 times as strong as the most compliant o1004 orientation11. As such, in response to loading, differently oriented NFA grains will sustain distinctly different elastic strains and thus generate intergranular strains to maintain compatibility with its neighbouring grains10,12. Plastic deformation starts within grains whose orientation, relative to the loading axis, facilitates the easy activation of their dislocation slip systems when the critical resolved shear stress is reached; other oriented grains, which are plastically hard, will respond to the load elastically and sustain more loads leading to stress redistribution among the grains. In the case of uniaxial tensile loading, the plastically deformed grains accumulate a compressive intergranular strain, whereas the grains lying along the plastically hard orientations accommodate a tensile intergranular strain10,13.

Thus, elastic anisotropy plays a predominant role in the early stage of plastic deformation, particularly for the orientation dependence of dislocation slip activity14, while the anisotropy-induced intergranular strains contain unique information about the dislocation activities and on the deformation processes12,13,15,16. As such, it is imperative to determine the elastic anisotropy and intergranular strains of NFAs to understand their grain-level anisotropic deformation behaviours.

In this study, we have characterized in situ the mesoscale anisotropic responses of selected grains in a 14YWT NFA to uniaxial tensile loading at a broad range of temperatures, by state-of-the-art neutron diffraction at the Spallation Neutron Source (SNS), Oak Ridge National Laboratory. From the responses to elastic deformation, the single-crystal elastic constants of 14YWT were determined by considering the effect of preferred grain orientations (or texture), showing a strong temperature-dependent elastic anisotropy. From those in response to extensive plastic deformation that manifest anisotropy-induced intergranular strains, a deformation crossover was evident from a low-temperature intragranular dislocation slip mechanism to that mediated by grain boundary activities at elevated temperatures.

ResultsTensile deformation at elevated temperatures. The material used in the present study was the ultrane-grained 14YWT NFA with a nominal composition of Fe-14Cr-0.16W-0.4Ti-0.14Y-0.4O (at.%) (see alloy synthesis in Methods). Atom probe tomography

and transmission electron microscopy characterizations showed that the 14YWT NFA has a grain size of 200400 nm and contains a high number density of Y-Ti-O NCs (1-2 nm in diameter) inhomogeneously distributed in the grains and along the grain boundaries4,17. As the alloy was made by hot-extrusion of mechanically alloyed powders, the 14YWT NFA tensile specimen exhibits a strong o1104 bre texture, as seen from the neutron diffraction patterns and also consistent with the texture characterization in earlier works18.

In situ neutron diffraction characterization of the 14YWT NFA under uniaxial tensile loading at room temperature (RT), 400, 600, 750 and 800 C, respectively, was performed on VULCAN19, the Engineering Materials Diffractometer at the SNS (see Methods, Supplementary Figs 14 and Supplementary Note 1 for more details of the experimental setup). The observed true stresstrue strain curve at each temperature is shown in Fig. 1a, and selected diffraction patterns collected during tensile loading at RT are presented in Fig. 1b,c. The neutron scattering geometry (see Supplementary Fig. 1) allows for simultaneous measurements of two scattering vectors Q1 and Q2, which are parallel and perpendicular to the tensile axis, respectively. Diffraction from the crystal lattice of grains whose plane normal (that is, ohkl4)

is parallel to Q1 or Q2 and satises Braggs law will be detected, and the lattice spacing (dhkl) can thus be measured from the

diffraction peak corresponding to the (hkl) reection. If treating the crystal lattice as a natural atomic plane strain gauge, lattice strains are determined from the shifts of lattice spacing due to an applied load, that is, ehkl dhkl d

0

hkl

d0hkl where d0hkl is the stress-free reference lattice spacing. Accordingly, two lattice strain components, one parallel to the tensile axis (that is, the loading direction, LD) and the other perpendicular to the tensile axis (that is, the transverse direction, TD), are measured simultaneously under tensile loading. As the macroscopic stress is uniaxial, only one stress component value is different from zero (s3 a 0, LD -3, TD - 1 where 13 refer to the three orthogonal axes). The two components of the grain-orientation-dependent strains recorded by our experimental setup are designated as e3(hkl) for the LD and e1(hkl) for the TD. The generic hkl indices denote not only a measurement direction in the crystal space, but also refer to specic families of grain orientation. The family of grains whose plane normal ohkl4 lies along the LD or the TD is designated as ohkl4//LD or ohkl4//TD, respectively. Figure 1d shows the variation of true stress with the lattice strain for six ohkl4//LD families of grains obtained during tensile loading at RT. Thus, neutron diffraction serves as a grain selector and can provide unique insights into the grain-scale deformation mechanisms.

Crystallographic texture. Beside lattice-spacing shifts, diffraction peaks can be used to assess the preferred orientation of grains (or texture). A virtual population of randomly oriented single crystals, generated by the quasi-Monte Carlo method, facilitates the realistic calculation of weighting factors for different powder diffraction peaks detected by each of the two instrument detector banks at VULCAN, as well as the correlation between the peaks recorded in different banks. The integrated intensities of a few independent diffraction peaks, after normalization to respective intensities obtained from a reference sample with random textures, are chosen to generate inverse pole gures (IPFs) with respect to the LD and TD. For a BCC structure, nine crystal-lographic directions dene a network for texture index estimations, and the fundamental triangle is partitioned in eight triangular cells (see Supplementary Fig. 5). The texture index inside each cell is dened using a barycentric interpolation method (see more details in Supplementary Note 2). As the initial texture of the NFA specimens is expected to have the axial

2 NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6178 ARTICLE

1,200

211

110

110

RT

1,000

Intensity

Intensity

400 C

True stress (MPa)

Intensity (a.u.)

True strain d-spacing ()

35 MPa

800

1.15 1.16 1.17 1.18 1.19 1.20

411

321

310

220

d-Spacing ()

1.96 2.00 2.04 2.08

d-Spacing ()

600

600 C

420

211

200

400

640 MPa

997 MPa

200

620 MPa

310 MPa

0

800 C

750 C

3 MPa

0.0 0.1 0.2 0.3 0.4 0.5

0.5 1.0 1.5 2.0 2.5 3.0

211

110

1,100

0.002 0.003 0.004 0.005 0.006 Lattice strain

Intensity

Intensity

1,000

Intensity (a.u.)

1.15 1.16 1.17 1.18 1.19 1.20

900

332

411 321

310

True stress (MPa)

222

d-Spacing ()

1.96 2.00 2.04 2.08

d-Spacing ()

220

211

200

110

800

(110) (200) (211) (310) (321) (420)

640 MPa

700

997 MPa

620 MPa

600

310 MPa

500

3 MPa

0.5 1.0 1.5 2.0 2.5 3.0

d-Spacing ()

35 MPa

Figure 1 | Strain measurements by in situ neutron diffraction. (a) True stresstrue strain curves for tensile loading of the 14YWT NFA at RT, 400, 600, 750 and 800 C, respectively. (b,c) Selected neutron diffraction patterns collected with respect to the loading and transverse directions duringtensile loading of an NFA at RT for a loading-unloading cycle (from bottom to top). (d) True stresslattice strain curves for six ohkl4//LD families of grains obtained during tensile loading at RT.

symmetry due to hot extrusion, and the tensile deformation preserves the same symmetry, the IPFs allow for a complete texture description and can be efciently used for averaging procedures in the analysis of stresses and strains. Two IPF maps with respect to the LD and TD are depicted in Fig. 2a,b, respectively, showing clearly the presence of a strong o0114 bre texture in the NFA specimens before tensile loading (see more IPFs in Supplementary Figs 6 and 7). Under tensile loading, the o0114//LD bre texture underwent a moderate strengthening.

Stress factors and single-crystal elastic constants. To illustrate the elastic anisotropy and intergranular strain accumulations in the NFA, the evolution of the lattice strains e3(hkl) of several families of grains (ohkl4//LD) under RT tensile loading has been shown in Fig. 1d. In the elastic regime, the lattice strain is linear as a function of the applied stress (s3), where the slope represents the elastic stiffness at the grain level and manifests the effect of cubic elastic anisotropy. Upon plastic deformation, the lattice strain begins to deviate from the elastic linearity and this effect is irreversible. The difference between the current value of lattice strain and the value extrapolated from the elastic response (that is, the reversible part of the lattice strain) is termed as intergranular strain as mentioned earlier, or Type II strain eIIhkl (ref. 11) and is conceptually equivalent with the residual strain recorded after unloading elastically from the current level of stress. In the following sections, the focus is placed on analyzing

[013]

[123]

[233]

LD (max:4.2)

TD (max:2.64)

[111]

4.0

[111]

3.5

2.5

[112]

3.0

2.5

[112]

2.0

[114]

2.0

1.5

1.5

[114]

1.0

1.0

0.5

0.5

[001]

[012]

[011]

[001]

[013]

[123]

[233]

[012]

[011]

LD stress factors (TPa1)

TD stress factors (TPa1)

[111]

5.6

[111]

2.5

2.0

[112]

5.3

[112]

1.5

[114]

5.0

[114]

1.0

0.5

4.7

[001]

[013]

[123]

[233]

[012]

[011]

[001]

[013]

[123]

[233]

[012]

[011]

Figure 2 | Texture characterization and stress-factor determination. Inverse pole gures with respect to (a) LD and (b) TD for the NFA before tensile loading, showing the strong presence of a o110//LD4 bre texture. Stress-factor inverse pole gures with respect to (c) LD and (d) TD for the NFA tensile loaded at 400 C.

NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6178

the linear regime of the lattice strain evolution, in order to determine the single-crystal elastic constants that give rise to the elastic anisotropy at different temperatures.

In diffraction stress analysis for polycrystalline aggregates, there are two ways of approaching the strain response at the grain level. The rst considers the proportionality between the grain-level average normal strain components (measured by diffraction along certain single crystal directions) and the macroscopic stress components applied at the specimen level. The proportionality factors are dened as stress factors. They can be calculated starting from single-crystal elastic constants by taking into account the crystallographic texture and orientation stereology (that is, grain shape and packing) using various averaging techniques2024. This approach was also consistently used to retrieve the single-crystal elastic constants from polycrystalline diffraction data, especially when the texture inuence is absent or weak. The second approach has been developed as a tool for residual stress determination in deformed specimens after the external load is released. The macroscopic stress components are not known in this case; it is thus necessary to dene a mesoscopic average stress eld depending on the grain orientations, known as the stress orientation distribution function2528. In this case, the proportionality factors between the strain and stress components are fully determined by the single-crystal elastic constants, although additional restrictive conditions or grain interaction models are often needed to obtain physically consistent solutions for the stress eld. In the present work, as the NFA specimens exhibit a strong o1104 bre texture as demonstrated in the preceding section, the second approach has been applied to retrieve the single-crystal elastic constants from the in situ neutron diffraction data. In this approach, the orientation-dependent stress acting at the grain level is a free parameter, but, to reduce the number of parameters, the average local stress is considered to be uniaxial and parallel to the applied stress29. The FEM simulation results reported by ref. 10 seem to support this assumption. As the analysis of lattice strains was performed in the elastic region only, the experimentally determined stress factors were used as inputs in the model, while elasticity equations for single-crystals with cubic symmetry were applied for uniaxial straining30.

In the elastic regime, the longitudinal lattice strain component, e3(hkl), is related to the applied stress through the stress factor F33(hkl),

e3

hkl F33hkls3 1 The transverse lattice strain, e1(h0k0l0), can be represented as a sum of transversally isotropic bres with the ohikili4 orientation along the LD,

e1

h0k0l0 X

i

X

i

pis3 4

where pi are the average stress concentration factors attached to different texture bres and

Fn n21n22 n22n23 n23n21 5

Dn; m n21m21 n22m22 n23m23 6

S 44 S44A 1 7

A

2 S11 S12

S44

fi mj

w

i S12

1

2 S 44D ni; mj

2c44

c11 c12

8

In the above equations, ni and mi are the hkl and h0k0l0 direction cosines in the cubic lattice cell, Sij and cij the coefcients of compliance and stiffness for the single crystal, and A the anisotropy parameter.

The basic bres include o1104, o2004, o2114, o3104, o2224, o3214, o4114, o4204 and o3324 along the LD, which are given in Supplementary Note 2 and represented in Supplementary Fig. 5. As the main bre component contained almost 50% of the grains aligned with the o0114 generic direction along the LD, the pi values were treated differently for o0114 and the rest of the components. The pi corresponding to minor components (o0014, o0134, o0124, o1124 and o1234) were correlated through a rational function of F(n) and the weighted average over all the pi values was constrained to 1 (that is, the average stress is identical to the applied true stress). The values of single-crystal elastic constants were adjusted to minimize the squared difference between the experimental stress factors values and the calculated ones. The IPFs were used to determine the volume fractions wi and the correlation factors fi(mj). Although the inuence of texture is generally challenging to quantify, the o0114//LD bre texture, in this NFA, indeed facilitated the observation of the distinct difference in the transverse stress factors between the o0114 and o1124 directions, as illustrated in Fig. 3. Following equation (4), if the tensile load is applied along the o1104 direction, D(n,m) varies from to 0 and the corresponding transverse compliances change from positive to negative (see Fig. 3a). For ferritic steels, this peculiarity was recognized by ref. 31. Because of the strong o0114//LD bre texture, in our case, the transverse stress factor for the o1104 direction can effectively go positive for large elastic anisotropy (Fig. 3b). The single-crystal cubic elastic constants, c11, c12 and c44, have been deduced from the single-crystal compliances (Sij) by matching the experimental and calculated stress factors; these values, together with the shear stiffness [c0 (c11 c12)/2] and the elastic anisotropy parameter

(A), are listed as a function of temperature in Table 1 and also shown in Fig. 4. Our estimation of relative errors of S11, S12 and

S44 amounts to 2% or less for RT, 400 and 600 C, but increases to

fih0k0l0wiF13 h0k0l0; hikili

s3

~F13h0k0l0s3

2

where wi is the volume fraction of the ith bre, fi(h0k0l0) the probability of grains belonging to the ith ber that has the oh0k0l04 direction along the TD, F13(h0k0l0,hikili) the stress factor for grains with oh0k0l04 along the TD and o hikili 4 along the LD, and ~F13h0k0l04 the average transverse stress factor.

The average transverse stress factor, ~F13h0k0l0, or the

longitudinal stress factor, F33(hkl), can be determined from the slopes of the lattice strain [e1(h0k0l0) or e3(hkl)] versus true stress curves recorded during loading or unloading in the elastic regime of deformation. Because of the strong initial o1104 texture (that is, B50% of the grains are oriented as o0114//LD), only six longitudinal stress factors could be determined in the present work, and this number amounts to eight for the transverse stress factors. The hkl stress factor distribution can be represented as an

IPF map. For instance, the stress factor IPFs collected at 400 C are shown in Fig. 2c,d for LD and TD, respectively (see Supplementary Figs 810 and Supplementary Note 3 for other temperatures). The experimental values of stress factors were obtained from averaging the measured slopes recorded during a few (24) loading/unloading cycles, as a systematic drift with the degree of plastic deformation was not discernible (that is, slow texture development during loading).

Stress factors can be calculated based on the single-crystal elasticity and by considering the adjustable stress concentration factors that link the local stress eld to the applied stress. As such, equations (1) and (2) can be written as

e3 ni

S11 S 44F ni

p

is3 3

e1 mj

4 NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6178 ARTICLE

5% or more for 750 and 800 C. A strong temperature-dependent c0 and A are evident as expected from those of pure Fe, which will be discussed in greater detail below.

Intergranular strains and mesoscale deformation. As indicated in Fig. 1d, the o1104//LD grains are the rst to yield and then take on compressive intergranular strains on plastic deformation at RT. This is also the case for other higher temperatures. The

intergranular strain accumulation in the o1104//LD grains thus contains mesoscale information characteristic of the plastic deformation. In what follows, we show how the o1104 inter-granular strain accumulates as plastic deformation proceeds at various temperatures, which provide insights into the deformation mechanisms. The applied true stresses as a function of lattice strain of the o1104//LD grains obtained at RT, 400, 600, 750 and 800 C, respectively, are shown in Fig. 5. Intergranular strains (DeII) accumulated in the o1104//LD grains (see the upper left inset of Fig. 5) are calculated by subtracting the current lattice strain (e) from the projected linear elastic strain corresponding to the current applied stress (s), that is, DeII e s=ET110 where

ET110 is the slope of the s versus e curve in the elastic regime at the temperature T. Normalized intergranular strain (DeII=e110y) as a function of normalized true stress (s=s110y) for each of the ve temperatures has been plotted in the lower right inset of

Fig. 5, where s110y is dened as the yield stress necessary for the o1104//LD grains to induce a plastic strain of 50 10 6 and

e110y is the corresponding lattice strain. At low temperatures (RT and 400 C), the o1104 intergranular strain was accumulated gradually and became saturated at a value of 0:25s110y when the

applied true stress increased continuously over 1:3s110y; this

clearly depicts the initial yielding of the o1104//LD grains (xed lattice strain and taking no more loads) followed by strain hardening (increased lattice strain and taking on loads elastically).

0 30 60 90

1.5

0.5

1.0

<332> <011>

<111>

<112>

<001>

<114>

S 12+1/2S* 44D(n i,m j) (TPa1)

0.4

0.5

0.3

0.0

D(n i,m j)

0.5

0.2

1.0

0.1

1.5

2.0

0.0

Angle ()

6.0

80

4

5.5

1,000

5.0

60

800

2

<011>

<112>

600

4.5

c(GPa)

TD stress factors (TPa1)

40

400

Yield strength (MPa)

A, Anisotropy parameter

4.0

20

Yield strength

200

0

3.5

c , experiment

Fit to c

0

0

3.0

200

2

0.2 0.4 0.6 0.8 1.0

2.5

T/Tc

4

2.0

1.5

6

1.0

14YWT NFA, present work Fe-10Cr, first-principles [Ref 32] Fe-14Cr, first-principles [Ref 33]

200 400 600 800 1,000 1,200

0.5

Temperature (K)

Figure 3 | Elastic anisotropy of transverse stress factors. (a) Geometric factor D in equation (4) (left axis) and the transverse anisotropy at RT (right axis) as a function of angle towards the o0114 direction, and (b)

Temperature dependence of transverse stress factor for the o0114 and o1124 directions. Symbols and lines in b denote experimental and calculated values, respectively.

T (C)

Figure 4 | Temperature-dependent anisotropy parameter. Comparison of the temperature-dependent anisotropy parameters (A) experimentally determined in the present work with those obtained from rst-principles calculations by refs 32,33. The inset shows c0 and the 0.2% yield stress as a function of temperature.

Table 1 | Temperature-dependent mechanical properties.

Temperature (C) c11(GPa) c12 (GPa) c44 (GPa) c0 (GPa) A K(GPa) r110y (MPa) e110y (10 4) ry (MPa) 25 273 150 106 62 1.7 4.6 784 35 940 400 249 143 103 53 1.9 4.2 597 29 726 600 232 145 100 43 2.3 3.8 321 19 381 750 195 153 81 21 3.9 2.6 91 8 87 800 188 156 85 16 5.3 2.3 59 6 47

Experimentally determined temperature-dependent elastic constants (c ), anisotropy parameter (A), energy factor (K) for the 001[100] edge dislocation, stress and lattice strain corresponding to yielding of the o1104//LD grains (s110

y and e

110y ), and macroscopic yield stress (s ) for the 14YWT NFA.

0 100 200 300 400 500 600 700 800 900

NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6178

1,200

1,000

800

600

400

200

<110>

1,100

1,000

True stress (MPa)

[afii9846]

900 [afii9846]y(110)

800

[afii9830]

[afii9830]y(110)

Yielding

700

True stress (MPa)

600

0.0030 0.0035 0.0040 0.0045

Lattice strain

400 C

0.2

Normalized intergranular strain

Onset of yielding

0.1

0.0

0.1

600 C RT

400 C

0.2

800 C

0.3

0.4

0

750 C

600 C

800 C

750 C

0.5

0.8 1.0 1.2 1.4 1.6

Normalized true stress

0.000 0.001 0.002 0.003 0.004 0.005

<110> Lattice strain

RT

Figure 5 | Evolution of intergranular strains. Applied true stress as a function of the o1104 lattice strain at various temperatures. The inset at the upper-left corner shows the onset of yielding (corresponding to a plastic strain of 50 me), and the denition of type II intergranular strain (DeII). The inset at the lower-right corner represents normalized o1104 intergranular strain as a function of normalized true stress for different temperatures during plastic deformation.

At the intermediate temperature (that is, 600 C), the inter-granular strain accumulation appears to be similar to what underwent at low temperatures and also turned to be saturated but at a higher value of 0:5s110y; however, in sharp contrast,

the intergranular strain saturation continued with a decrease in the applied true stress from a value of 1:55s110y. This beha

viour indicates a mechanism of initial yielding followed by strain softening due to extensive plastic deformation. At high temperatures (750 C and 800 C), intergranular strains were accumulated immediately after the onset of yielding but, strikingly, with the decrease in the applied true stress, indicative of strain softening.

DiscussionThe experimentally determined anisotropy parameter (A) is shown in Fig. 4 as a function of temperature (T), which agrees remarkably well with the predictions obtained from rst-principles calculations by refs 32,33. This result validates the above-mentioned method that was developed to extract elastic constants from textured polycrystalline materials. The marked increase in A at elevated temperatures is attributed to a sudden drop of the shear stiffness constant (c0) as T approaches a critical value (Tc); this has been depicted in the inset of Fig. 4, together with a best t to the experimental data by the function: c0T 77:91 TTc1=2. Tc is determined to be 1,127 K (854 C)

for the NFA, in comparison with that of 1,185 K (912 C) for pure

Fe with c0 T 56:8 1 TTc

c44c0

c12 2c0c12 c44 c0 1=2 (refs 34,37). The calculated K-values, as listed in Table 1, appear to scale with the yield strength, although K shall be slightly smaller (o10%) for temperatures below 600 C where the a/2o1114 dislocations (for example, 11111 2 ) dictate the yielding34. This behaviour

suggests that, in addition to less dispersion strengthening from the NCs and possible grain-boundary-mediated softening (for example, dislocation adsorption, and so on) at increased temperatures9, the yield strength drop could be also induced by the reduction of the shear stiffness (c0) or, in other words, the large elastic anisotropy. Qualitatively, this nding may be understood in the framework of refs 34,37 which showed that dislocation pile-ups can be markedly inuenced by c0 when full anisotropic elasticity is considered.

The anisotropy-induced intergranular strains/stresses in the o1104//LD grains at RT, 400 and 600 C show a nearly identical accumulation process during the onset of yielding and in the subsequent early-stage plastic deformation regime (see the inset of Fig. 5). This intriguing similarity indicates that the same yielding mechanism is dominated by dislocation slip inside the grains, as manifested by the unaltered lattice strains and the increased intergranular strains/stresses. However, at high temperatures from 750 to 800 C, the onset of yielding is immediately followed by lattice softening as seen from decreased lattice strains; this implies that these grains sustain less and less loads, and the plastic deformation is mediated by grain-boundary activities, although the Y-Ti-O NCs and vacancy-solute interactions retard diffusion40. It is also noted that, in response to higher values of plastic deformation, a deformation crossover exists from low-temperature strain hardening to high-temperature strain

12 (ref. 34). For pure Fe, Tc is the bcc-fcc or a-g transition temperature as proved by ref. 34. The displacive a-g transition has been found to be associated with the development of the (110) soft transverse acoustic phonon mode, corresponding to shear deformation of the BCC unit cell along the phase transformation pathway3437. For the 14YWT NFA, however, the a-g transition and the presence of Tc are unexpected. Nevertheless, this transition was experimentally observed in the in situ neutron diffraction studies of a 14YWT NFA, which showed a face-centred cubic g-phase occurring at

B910 C and then disappearing at B810 C when cooling from 1,300 C to RT (see Supplementary Fig. 11). This transition is likely when considering the Fe-Cr binary phase diagram, which exhibits an (a g) two-phase region between 11.2 wt.% Cr and

13.4 wt.% Cr ranging from B846 C up to 1,200 C. Additions of other alloying elements (for example, Y and Ti, and so on.) appear to widen the (a g) region; this makes 14YWT lying,

most likely, in the vicinity of the a/(a g) boundary, as indicated

by the observed small volume fraction (that is, o1.5%) of the transformed g from neutron diffraction. Thus, the a-g transition in 14YWT appears to be marginal. However, no matter how small the transformed-g amount is, the phase transition, if present, kinetically requires the whole a phase (including every single a grain) to be elastically softening, so as to facilitate the nucleation of g from a with a specic crystallographic orientation relationship. This softening-associated phase transition provides a reasonable explanation of the critical temperature Tc, as

expected from our experimentally determined c0(T).

Also shown in the inset of Fig. 4 is the 0.2% macroscopic yield strength (sy) determined from the engineering stressstrain curves. Like c0, sy exhibits a strong temperature dependency, particularly a considerable yield strength drop above 400 C, which is consistent with previous studies in the literature9,38,39. As described earlier, the yield strength is essentially determined by the grains whose orientation, relative to the loading axis, favours easy activation of dislocation slip. For this NFA, the o1104//LD grains yield rst upon plastic deformation; the microscopic yield strength (s110y) with respect to the o1104//LD grains also decreases signicantly as the temperature increases, as shown in Table 1. It is known that the dislocation character determines the onset of crystallographic slip, and the elastic energy factor (K) of a dislocation may be considered as a strength indicator34,37. For example, the energy factor for straight edge dislocations with a Burgers vector of a/100S and a dislocation line direction of [001] (denoted here as 100[001]) can be evaluated as K

c12 c0

4p

6 NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6178 ARTICLE

softening; the crossover temperature lies between 400 and 600 C. This crossover phenomenon demonstrates that intragranular dislocation activities dictate the late-stage of plastic deformation at low temperatures, whereas intergranular grain boundary activities are predominant at high temperatures.

It should be noted that the volume fraction of NCs in a 14YWT NFA can be estimated to be o0.1%, in consideration of their sizes of B1-2 nm diameter with a number density of 1 1023 m 3

(refs 9,17,41). Moreover, these NCs have been characterized to be coherent with the ferrite matrix5,17. Such a small volume fraction, together with the nanocluster-matrix coherency, would give rise to negligibly small mist strains and stress concentrations, if any, between NCs and the ferrite matrix during elastic deformation. Thus, NCs are believed to have insignicant effect on the single-crystal elastic constants of the ferrite matrix. This may explain why there is no great difference in elastic constants between NFAs and Fe (or conventional steels), as assumed in many studies9. However, there has been no direct measurement of the single-crystal elastic constants of NFAs to prove this hypothesis. In this regard, our diffraction results provide the rst experimental evidence for validating this scenario.

In summary, anisotropic responses of the ultrane 14 YWT NFA grains when subjected to tensile stresses at a broad range of temperatures were investigated using in situ neutron diffraction. Temperature-dependent single-crystal elastic constants for the NFA were determined using a new model developed in the present work that takes into account the effect of preferred grain orientations. A strong temperature-dependent elastic anisotropy is revealed as a result of a sharp decrease in the shear stiffness constant, indicating elastic softening and instability of the ferritic matrix at elevated temperatures. The anisotropy-induced intergranular strain/stress accumulations reveal a crossover phenomenon of plastic deformation from a low-T intragranular dislocation slip mechanism to that mediated by grain boundary activities at high-T. This study provides unique insights into the grain-level mesoscale deformation mechanisms.

Methods

Alloy synthesis. The 14YWT NFA was made by mechanical alloying of the mixture of Fe-Cr-W-Ti alloy powder and Y2O3 powder, followed by canning in an evacuated jacket and hot extrusion at 850 C. The hot-extruded ingot was then annealed for 1 h at 1,000 C. One heat was used to produce this NFA, whose nominal composition is Fe-14Cr-0.16W-0.4Ti-0.14Y-0.4O (at.%). The details of alloy processing have been described elsewhere2.

Neutron diffraction experimental setup. In situ neutron diffraction experiments were performed on VULCAN19, the Engineering Materials Diffractometer at the SNS, Oak Ridge National Laboratory. Threaded dogbone-shaped NFA specimens of 2.90 mm in diameter and 20.6 mm in gauge-length were uniaxially deformed under tension using a Materials Testing System load-frame. The macroscopic engineering strain for the measurement at RT was recorded by an extensometer with a nominal gauge length of 12.7 mm. An induction heater attached to the load-frame was used for heating each specimen to a preset temperature, at which tensile loading of the specimen was conducted while being monitored in situ by neutron diffraction (see Supplementary Fig. 2). Four ne-gauge K-type thermocouples were spot weld onto each specimen to control and monitor the temperatures at locations near the specimen centre and the two ends. A uniform temperature zone of B10 mm in the centre of each specimen was established by ne-tuning the induction heater before loading. This homogeneous zone, having a temperature gradient of o3 C and being the hottest part of the specimen, was probed in situ by the neutron beam under tensile loading with a strain rate of B8 10 5 s 1. True

stresses and true strains during plastic deformation at high temperatures were evaluated from measured engineering stresses, together with engineering strains estimated from the relative cross-head displacements after a linear correction (see more details in Supplementary Note 1).

As shown in Supplementary Fig. 1, the tensile stress axis of the load-frame was mounted at 45 from the incident neutron beam direction. Diffraction patterns were recorded by the two detector banks (namely, B1 and B2) positioned at diffraction angles of / 90. As such, B1 and B2 collected diffraction from the

lattice planes perpendicular and parallel to the LD, respectively. The sampling

volume was dened by the incident beam slits (5 mm width and 3 mm height) and the radial collimators located between the specimen and the detector banks restricting the eld of view to 5 mm along the neutron beam path. The neutron diffraction data were collected continuously during the tensile loading at various temperatures. By using the data reduction software VDRIVE42, the recorded data were then sliced and binned into histograms corresponding to small temporal intervals that were synchronized with the loading parameters (force, displacement, temperature, and so on.). Singlepeak tting, employing the GSAS program in VDRIVE, was used to extract the peak position, width and intensity corresponding to each hkl reection. More details about the neutron diffraction experimental setup are available in Supplementary Note 1.

References

1. Zinkle, S. J. & Busby, J. T. Structural materials for ssion & fusion energy. Mater. Today 12, 1219 (2009).

2. Hoelzer, D. T. et al. Inuence of particle dispersions on the high-temperature strength of ferritic alloys. J. Nucl. Mater. 367, 166172 (2007).

3. Miller, M. K., Parish, C. M. & Li, Q. Advanced oxide dispersion strengthened and nanostructured ferritic alloys. Mater. Sci. Tech. 29, 11741178 (2013).

4. Yao, L., Ringer, S. P., Cairney, J. M. & Miller, M. K. The anatomy of grain boundaries: Their structure and atomic-level solute distribution. Scripta Mater. 69, 622625 (2013).

5. Hirata, A. et al. Atomic structure of nanoclusters in oxide-dispersion-strengthened steels. Nat. Mater. 10, 922926 (2011).

6. Fu, C. L., Krcmar, M., Painter, G. S. & Chen, X. Q. Vacancy mechanism of high oxygen solubility and nucleation of stable oxygen-enriched clusters in Fe. Phys. Rev. Lett. 99, 225502 (2007).

7. Brandes, M. C., Kovarik, L., Miller, M. K., Daehn, G. S. & Mills, M. J. Creep behavior and deformation mechanisms in a nanocluster strengthened ferritic steel. Acta Mater. 60, 18271839 (2012).

8. Schneibel, J. H. et al. Ultrane-grained nanocluster-strengthened alloys with unusually high creep strength. Scripta Mater. 61, 793796 (2009).

9. Schneibel, J. H., Heilmaier, M., Blum, W., Hasemann, G. & Shanmugasundaram, T. Temperature dependence of the strength of ne- and ultrane-grained materials. Acta Mater. 59, 13001308 (2011).

10. Wong, S. L. & Dawson, P. R. Inuence of directional strength-to-stiffness on the elastic-plastic transition of fcc polycrystals under uniaxial tensile loading. Acta Mater. 58, 16581678 (2010).

11. Hutchings, M. T., Withers, P. J., Holden, T. M. & Lorentzen, T. Introduction to the characterization of residual stress by neutron diffraction (Taylor&Francis, 2005).

12. Pang, J. W. L., Holden, T. M., Turner, P. A. & Mason, T. E. Intergranular stresses in Zircaloy-2 with rod texture. Acta Mater. 47, 373383 (1999).

13. Pang, J. W. L., Holden, T. M. & Mason, T. E. In situ generation of intergranular strains in an A17050 alloy. Acta Mater. 46, 15031518 (1998).

14. Clausen, B., Lorentzen, T. & Leffers, T. Self-consistent modelling of the plastic deformation of FCC polycrystals and its implications for diffraction measurements of internal stresses. Acta Mater. 46, 30873098 (1998).

15. Holden, T. M., Holt, R. A. & Pang, J. W. L. Intergranular stresses in ZIRCALOY-2. Metall. Mater. Trans. A 33, 749755 (2002).

16. Daymond, M. R., Tome, C. N. & Bourke, M. A. M. Measured and predicted intergranular strains in textured austenitic steel. Acta Mater. 48, 553564 (2000).

17. Brandes, M. C., Kovarik, L., Miller, M. K. & Mills, M. J. Morphology, structure, and chemistry of nanoclusters in a mechanically alloyed nanostructured ferritic steel. J. Mater. Sci. 47, 39133923 (2012).

18. Praud, M. et al. Study of the deformation mechanisms in a Fe-14% Cr ODS alloy. J. Nucl. Mater. 428, 9097 (2012).

19. Wang, X. L. et al. VULCANThe engineering diffractometer at the SNS. Phys. B 385, 673675 (2006).

20. Gnaeupel-Herold, T., Brand, P. C. & Prask, H. J. Calculation of single-crystal elastic constants for cubic crystal symmetry from powder diffraction data.J. Appl. Cryst. 31, 929935 (1998).21. Howard, C. J. & Kisi, E. H. Measurement of single-crystal elastic constants by neutron diffraction from polycrystals. J. Appl. Cryst. 32, 624633 (1999).

22. Bunge, H. J., Kiewel, R., Reinert, T. & Fritsche, L. Elastic properties of polycrystalsinuence of texture and stereology. J. Mech. Phys. Solids 48, 2966 (2000).

23. Matthies, S., Priesmeyer, H. G. & Daymond, M. R. On the diffractive determination of single-crystal elastic constants using polycrystalline samples.J. Appl. Cryst. 34, 585601 (2001).24. Gnaeupel-Herold, T., Creuziger, A. A. & Iadicola, M. A model for calculating diffraction elastic constants. J. Appl. Cryst. 45, 197206 (2012).

25. Behnken, H. Strain-function method for the direct evaluation of intergranular strains and stresses. Phys. Stat. Solidi A 177, 401418 (2000).

26. Popa, N. C. & Balzar, D. Elastic strain and stress determination by Rietveld renement: generalized treatment for textured polycrystals for all Laue classes.J. Appl. Cryst. 34, 187195 (2001).

NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6178

27. Wang, Y. D., Wang, X. L., Stoica, A. D., Richardson, J. W. & Peng, R. L. Determination of the stress orientation distribution function using pulsed neutron sources. J. Appl. Cryst. 36, 1422 (2003).

28. DeWit, R. Diffraction elastic constants of a cubic polycrystal. J. Appl. Cryst. 30, 510511 (1997).

29. Ungar, T., Stoica, A. D., Tichy, G. & Wang, X. L. Orientation-dependent evolution of the dislocation density in grain populations with different crystallographic orientations relative to the tensile axis in a polycrystalline aggregate of stainless steel. Acta Mater. 66, 251261 (2014).

30. Norris, A. N. Poissons ratio in cubic materials. Proc. R. Soc. A 462, 33853405 (2006).

31. Oliver, E. C., Daymond, M. R. & Withers, P. J. Interphase and intergranular stress generation in carbon steels. Acta Mater. 52, 19371951 (2004).

32. Razumovskiy, V. I., Ruban, A. V. & Korzhavyi, P. A. Effect of temperature on the elastic anisotropy of pure Fe and Fe0.9Cr0.1 random alloy. Phys. Rev. Lett. 107, 205504 (2011).

33. Zhang, H. L. et al. Elastic anomalies in Fe-Cr alloys. J. Phys.Condes. Matter 25, 195501 (2013).

34. Dudarev, S. L., Bullough, R. & Derlet, P. M. Effect of the alpha-gamma phase transition on the stability of dislocation loops in bcc iron. Phys. Rev. Lett. 100, 135503 (2008).

35. Hasegawa, H., Finnis, M. W. & Pettifor, D. G. A calculation of elastic-constants of ferromagnetic iron at nite temperatures. J. Phys. F 15, 1934 (1985).

36. Hasegawa, H., Finnis, M. W. & Pettifor, D. G. Phonon softening in ferromagnetic bcc iron. J. Phys. F 17, 20492055 (1987).

37. Fitzgerald, S. P. & Dudarev, S. L. Dislocation pile-ups in Fe at high temperature. Proc. R. Soc. A 464, 25492559 (2008).

38. Kim, J. H., Byun, T. S., Hoelzer, D. T., Kim, S. W. & Lee, B. H. Temperature dependence of strengthening mechanisms in the nanostructured ferritic alloy 14YWT: Part I-Mechanical and microstructural observations. Mater. Sci. Eng. A 559, 101110 (2013).

39. Kim, J. H. et al. Temperature dependence of strengthening mechanisms in the nanostructured ferritic alloy 14YWT: Part II-Mechanistic models and predictions. Mater. Sci. Eng. A 559, 111118 (2013).

40. Miller, M. K., Fu, C. L., Krcmar, M., Hoelzer, D. T. & Liu, C. T. Vacancies as a constitutive element for the design of nanocluster-strengthened ferritic steels. Front. Mater. Sci. China 3, 914 (2009).

41. Miller, M. K., Russell, K. F. & Hoelzer, D. T. Characterization of precipitates in MA/ODS ferritic alloys. J. Nucl. Mater. 351, 261268 (2006).

42. An, K. et al. First in situ lattice strains measurements under load at VULCAN. Metall. Mater. Trans. A 42, 9599 (2011).

Acknowledgements

This research was supported by the US Department of Energy, Division of Materials Sciences and Engineering. The neutron scattering work at the SNS was sponsored by the Scientic User Facilities Division, Ofce of Basic Energy Sciences, US Department of Energy, at Oak Ridge National Laboratory under contract DE-AC05-00OR22725 with UT-Battelle. The authors (G.M.S. and D.M.) thank Dr E.A. Payzant for useful discussions. A.D.S. and D.M. acknowledge the conceptual contribution of Professor X.-L. Wang (City University Hong Kong) to the early stage of this study. We also thank Harley Skorpenske for technical support at VULCAN.

Author contributions

G.M.S., A.D.S. and D.M. designed and performed the experiments and analysed the data. D.M. coordinated the scattering experiments. M.K.M. developed and supervised the research programme, selected and provided additional microstructural data of the 14YWT NFA. A.D.S. and D.M. wrote the manuscript after discussions and with contributions from all the authors.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Stoica, G. M. et al. Temperature-dependent elastic anisotropy and mesoscale deformation in a nanostructured ferritic alloy. Nat. Commun. 5:5178 doi: 10.1038/ncomms6178 (2014).

8 NATURE COMMUNICATIONS | 5:5178 | DOI: 10.1038/ncomms6178 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.