Scholarly article on topic 'The Influence of Grain Interactions on the Plastic Stability of Heterophase Interfaces'

The Influence of Grain Interactions on the Plastic Stability of Heterophase Interfaces Academic research paper on "Nano-technology"

Share paper
Academic journal
OECD Field of science

Academic research paper on topic "The Influence of Grain Interactions on the Plastic Stability of Heterophase Interfaces"

Materials 2014, 7, 302-322; doi:10.3390/ma7010302



ISSN 1996-1944


The Influence of Grain Interactions on the Plastic Stability of Heterophase Interfaces

Jason R. Mayeur *, Irene J. Beyerlein, Curt A. Bronkhorst and Hashem M. Mourad

Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA; E-Mails: (I.J.B.); (C.A.B.); (H.M.M.).

* Author to whom correspondence should be addressed; E-Mail:; Tel.: +1-505-665-3015.

Received: 1 November 2013; in revised form: 13 December 2013 / Accepted: 17 December 2013 / Published: 13 January 2014

Abstract: Two-phase bimetal composites contain both grain boundaries and bi-phase interfaces between dissimilar crystals. In this work, we use a crystal plasticity finite element framework to explore the effects of grain boundary interactions on the plastic stability of bi-phase interfaces. We show that neighboring grain interactions do not significantly alter interface plastic stability during plane strain compression. The important implications are that stable orientations at bimetal interfaces can be different than those within the bulk layers. This finding provides insight into bi-phase microstructural development and suggests a pathway for tuning interface properties via severe plastic deformation.

Keywords: crystal plasticity; rolling texture; nanolamellar composites; interface stability; grain interactions

1. Introduction

Bulk bimetal nanocomposites are generating much interest due to their high strength and hardness [1-10]. Recent experiments show that these bimetal nanocomposites, in contrast to bulk single phase nanocrystalline metals, are surprisingly stable in extreme conditions, such as high temperature, shock, high temperature deformation, severe plastic deformation, and radiation [1,4,11-14]. Unlike traditional composites, these nanomaterials contain a high density of interfaces, spaced nanometers apart. In addition to being strong obstacles to dislocation motion the bimetal interfaces are clearly playing a key

role in defect 'management'; however, very little is known about the physical properties of the interfaces themselves. Only very recently have atomic-scale studies revealed how deformation mechanisms, such as slip transmission, plastic slip and deformation twinning, are affected by interface structure [15-17].

Bulk bimetal nanocomposites can be fabricated by severe plastic deformation (SPD) techniques and the interface characteristics evolve as a function of the imposed deformation [1,2]. Bimetal interfaces within nanolayered composites fabricated by accumulative roll bonding (ARB) are found to stabilize to a particular interface type [18,19]. If the ARB processing pathway is changed, then another stable predominant interface prevails [20]. Both of these interfaces are distinct from those that develop in nanolayered composites fabricated by physical vapor deposition (PVD) [18]. In the case of PVD composites, the dominant interface character is driven by low formation energy. However, in SPD composites plastic stability also plays a role in the emergence of the predominant interface. Plastic stability describes the stability of the bimetal interface character (five-parameter description of a boundary) as a function of plastic deformation. Thus, to be plastically stable, both crystals comprising an interface must be able to slip without reorienting. The stability of crystal orientations for various monophase fcc and bcc single and polycrystalline metals have been well characterized experimentally and investigated using crystal plasticity models [21-26]. In contrast, much less is known about the plastic stability of bi-phase interfaces. Developing an understanding of the plastic stability of such interfaces is critical for designing processing methodologies for bulk bimetal nanocomposites.

Crystals attached to bi-phase interfaces have different microstructural 'neighborhoods' and this can affect its plastic stability. For example, ARB fabrication begins with a stack of polycrystalline sheets, each of which contain multiple grains through the thickness. As the composite is rolled, the individual layer thickness decreases as does the number of grains spanning each layer. At some critical layer thickness, only a single grain spans a given layer [10]. With further deformation, the average grain aspect ratio increases [19], such that each grain is primarily bounded by two bimetal interfaces. Thus, for composites that contain multiple grains through the thickness, interface grains are surrounded by neighboring grains of the same phase within each layer as well as grains of dissimilar phase along the interface, whereas each grain in the fine layered composites (single grain through the thickness) is primarily bounded by grains of the dissimilar phase. Therefore, it is expected that the plastic stability of the interface could vary as a function of layer thicknesses due to the potential influence or lack thereof grain interactions within individual layers.

In this work, we use a bi-phase crystal plasticity finite element model to explore the competition between grain boundary interactions and heterophase interactions on the plastic stability of a bimetal interface. To compare different interface configurations, we develop stability parameters based on the influence of interface misorientation and misorientation rates with strain. The results show that the bimetal interface constraints have a profound effect on the plastic stability. The implications of this study are that in plastically deformed bimetal composites, the grains adjoining interfaces have different orientation distributions than those within the layers. These results help to explain experimental observations of texture gradients reported in multilayered composites deformed in rolling [20]. The results of our computational study also indicate that when the layers refine to single crystalline layers via plastic deformation, the texture is primarily governed by the interface.

2. Stability of Rolling Textures

Before moving forward, it is instructive to briefly review prior work on the stability of grains as single crystals, embedded in a polycrystalline, or surrounded by bimetal interfaces. The evolution of stable orientations in rolling has been extensively studied for single phase materials. A stable orientation or orientation distribution (texture) is classified as one in which the crystal orientations do not substantially change with increasing strain. Experimentally and computationally, stable orientations can be characterized by non-decreasing change in intensity and/or volume fraction as a function of deformation [23,25,27]. These studies have underscored the importance of a number of factors that affect the texture evolution such as the imposed strain level, initial texture, grain morphology, and the slip system level hardening behavior. They found that these factors influence the inter- and intragranular slip system activity both in terms of the number of active (or dominant) systems and the relative proportions of accumulated slip on those systems, which directly drives crystal reorientation.

3. Stability of Layered Bi-Phase Composites in Rolling

Recently crystal plasticity formulations have been employed to study multilayers Cu-Nb, Cu-Ni, Zr-Nb, and Cu-Ag composites under deformation [4,28-34]. Anderson et al. [28] and Al-fadhalah et al. [29] studied Cu-Nb multilayered composites produced by PVD, which contain thermodynamically stable fcc-bcc interfaces, e.g., Kurdjumov-Sachs (KS) ND || (111)fcc || (101)bcc and II (110)fcc || (111)bcc in-plane and Nishiyama-Wasserman (NW) ND || (111)fcc || (110)bcc and (110)fcc || (001)bcc in-plane interfaces [35]. Using a composite grain model within a viscoplastic self-consistent (VPSC) scheme wherein all Cu and Nb crystals were initially joined with KS orientation relationship with (111)fcc ||( 110)bcc || ND and random in-plane orientations of the composite grains with respect to the RD, Al-fadhalah et al. [29] showed the KS interface was unstable and that the expected rolling texture characteristic of each phase developed. A separate set of simulations was performed with the constraint that the KS OR be maintained during rolling, and the results compared favorably with the experimentally measured textures of samples with 75 nm layer thickness. This result suggests that the texture evolution in these materials is intimately tied to the preservation of the KS OR, and it has been argued that the KS OR is maintained due to the symmetric activation of slip systems in the individual layers. Anderson et al. [28] used a crystal plasticity model based on the principle of minimum shear to study the stability of Cu and Nb single crystals with [111]fcc || [110]bcc || ND under rolling conditions in an effort to understand the development of the atypical rolling textures in PVD Cu-Nb. By systematically varying the initial in-plane orientations (perturbations of each crystal by a rotation about the ND), they were able to identify several mutually stable orientations as well as predict the primary texture component that was experimentally observed.

More recently, plastic stability studies have been carried out in which the interface characters were allowed to evolve [30-33,36]. In Knezevic et al. [33], the VPSC scheme was used to evaluate texture development of the Zr phase in multilayered ARB Zr-Nb composites with micron layer thickness. The number of grains through the thickness of each layer was so large that interfaces had little effect on texture evolution in the Zr phase. In [30], crystal plasticity finite element (CPFE) modeling was employed to study texture development in ARB Cu-Nb composites with layers that were several microns

thick. That study demonstrated that the bi-phase CPFE model was able to reproduce classical rolling textures for each phase that are observed for ARB Cu-Nb composites with the layer thickness > 1 fj.m. In both studies, the possibility of texture gradients due to interface effects were not studied computationally or experimentally.

In another study by Mayeur et al. [31], the two-phase CPFE method was applied to examine the intrinsic stability (i.e., in the absence of grain interaction effects) of a bicrystal during rolling. The bicrystal configuration was a periodic cell containing a single interface character, and thus models a multilayered composite containing single crystalline layers. Thus, grain interaction effects on plastic stability were not considered. That work revealed that the plastic stability of the interface is essentially dictated by the least stable of the two crystal orientations comprising the interface. These results helped explain the emergence of special interface types observed in nanolayered composites, which were associated with texture components that deviated from those seen in the composites with layer thicknesses on the order of microns. This model was also used to explain the differences observed between the strengths of ARB and PVD Cu-Nb composites in elevated temperature deformation [4].

Very recently, experiments on ARB composites with layers several microns thick pointed to a through-thickness texture gradient [20]. Understanding how such texture gradients evolved during the rolling process requires understanding the competition between interface interactions and grain-grain interactions on the plastic stability of grains at the interface. This is an issue that has not been previously studied. To this end, in the present work we use the bi-phase CPFE framework used in Hansen et al. [30] and the bicrystal model of Mayeur et al. [31] to study plastic stability of bi-phase interfaces under rolling. Herein, simulations are performed for different grain and interface configurations to systematically elucidate the effect of nearest neighbor grain interactions and its relative contribution to the plastic stability of interface character.

4. Crystal Plasticity Finite Element Modeling

The crystal plasticity constitutive model is used to study two microstructural configurations as shown in Figure 1. The first configuration is the 'multicrystal' model, which consists of two multicrystalline layers, one of Cu and the other Nb. There are multiple grains along the interface of each layer and the interface character (orientation relationship and interface plane) of one pair of joined Cu and Nb grains is given a specific relationship. The remaining grains are assigned orientations sampled from the measured texture of each phase [30]. The second configuration is the bicrystal model used by Mayeur et al. [31], which does not account for nearest neighbor grain interactions. The bicrystal model is employed to study the same set of interfaces in order to elucidate the role of neighbor grain interactions on plastic stability. The specific interfaces studied in the present work (discussed in Section 5) are chosen because they are experimentally observed in fine scale Cu-Nb multilayered composites, e.g., either the ARB Cu-Nb composites with submicron layer thickness or the PVD Cu-Nb composites.

Figure 1. Finite element models: (a) periodic bicrystal; and (b) multicrystal grain morphology.

V4 ....................

(a) (b)

4.1. Constitutive Model

The constitutive model [37] used to describe the response of the Cu and Nb phases is a classical local finite deformation formulation of single crystal plasticity based on a two-term multiplicative decomposition of the deformation gradient into elastic, Fe, and plastic, Fp, parts, i.e.,

F = FeFp (1)

Dislocation motion, which leaves lattice vectors unaltered, is encapsulated by Fp and the stretching and rotating of the lattice is given by Fe. The evolution equation for Fp is defined in terms of the plastic velocity gradient with respect to the intermediate configuration as:

Lp = FpFp-1 = ^ 7«s« 0 na (2)

where Afa is the slip system shearing rate, sa is the slip vector, and na is the slip plane normal. Since the lattice vectors are assumed to be unaltered by dislocation glide, the lattice slip vector, sa, and the slip plane normal, na, are given in the current configuration as:

s« = Fesa , na = Fe na (3)

Equation (3) provides the means to track the evolution of the crystallographic orientations of the crystal, which are embedded in the rotational part of Fe that can be obtained via the polar decomposition.

The kinetic equation for the slip system shearing rate is based on thermally-activated dislocation motion, i.e.,

Y = yoexp

. k0\ \

aj±\ P\ q

sgn(ra) (4)

where Yo is the reference shearing rate, F0 is the activation free energy, k is Boltzmann's constant, 9 is the absolute temperature, ^ is the shear modulus, is the shear modulus at 0 K, ta is the resolved

shear stress, sa is the slip resistance, s^ is the intrinsic lattice resistance, and p and q are exponents that characterize the shape of the obstacle profile. The slip resistance is governed by the following equations

^ haf Y | , sa (t = 0)

haf = [r +(1 - r) Saf ] hf

hf = h

sf - sf

where sS0 is the initial slip resistance, haf is the hardening interaction matrix, r is the latent hardening ratio, h0 is the initial hardening modulus, sf is the saturation slip resistance, and sS and A are material parameters that control the rate and temperature dependence of the saturated value of the slip resistance.

In prior work, the stability of individual fcc and bcc orientations have been found to depend on the choice of hardening parameters used in phenomenological hardening laws, such as rate sensitivity parameters, latent hardening coefficients, or in the particular case of bcc metals, the ratio of critical resolved shear stresses of the {110}(111) slip to |112}(111) slip systems [21-23,25,27]. In our simulations, we consider the {111}(110) slip systems for Cu and the {110}(111) and {112}(111) slip systems for Nb. Including both {110} and {112} slip systems is critical for capturing the plastic anisotropy of bcc metals as recently demonstrated by Wang and Beyerlein [38] in their discrete dislocation dynamics study. The constitutive parameters in the flow rule and hardening relations have been determined by fitting the single phase Cu and Nb responses to experimental data over a wide range of strain rates and temperatures. The model calibration was performed by Hansen et al. [30] and the material parameters are listed in Appendix (Table 2).

4.2. Bicrystal Model

The bicrystal model, shown in Figure 1a, consists of a Nb crystal joined to a Cu crystal with the interface normal in the x2-direction (ND), subjected to plane strain compression boundary conditions. Periodicity of the displacements is enforced [39] in both the x1 and x2 directions and the boundary conditions at the control vertices are illustrated in the schematic. In the simulations the Nb crystal is stacked on top of the Cu crystal, however, the stacking sequence is not important due to the periodicity constraints. Displacements are fully constrained at vertex 1, the x2 displacement component is constrained at vertex 2, and the xi component of displacement is fixed at vertex 3. Loading is applied at vertex 3 via displacement control at a constant strain rate of 0.08 s-1.

4.3. Multicrystal Model

The grain morphology for the multicrystal finite element model is shown in Figure 1b. The primary grains of interest that have the prescribed interface OR are located in the center of the volume element and are denoted by labels "1" and "2", where grain 1 is in the Cu layer and grain 2 is in the Nb

layer. Five different multicrystal realizations were considered for each interface. Each realization has a different orientation distribution for the grains surrounding the prescribed interface grains. The grain morphology and orientation distributions were constructed to be consistent with measurements obtained for an as-rolled ARB Cu-Nb composite material with 24 fj.m layer thickness as described by Hansen et al. [30]. The microstructural volume element has dimensions of 50 ^m x 100 ^m and contains 26 grains (14 Cu and 12 Nb). The grains are meshed using 6-noded triangular hybrid elements (CPE6H). The rolling process is modeled as plane strain compression with isochoric boundary conditions that are enforced via multipoint constraints. The rolling direction is parallel to the xi-axis and the normal direction is parallel the x2-axis. The microstructural statistical volume element was subjected to 50% rolling reduction at a strain rate of 0.1 s-1.

5. Results

The interfaces studied are summarized in Table 1. Each orientation is expressed in the rolling plane-rolling direction convention, where {hkl} is the rolling plane and {uvw) the rolling direction. The"dashed" shorthand notation for the interface lists the orientation (or component) for the Cu crystal before the dash and then the orientation for the Nb crystal after the dash. For example, the C-I interface refers to the C orientation of the Cu crystal and the I orientation of the Nb crystal. KS is short for Kurdjumov-Sachs (see Section 3), the interface characteristic of PVD Cu-Nb multilayer composites. The configurations listed in Table 1 are representative of interfaces that according to the prior calculations are intrinsically metastable (C-I, D-I), stable (Brass-I,Goss-I), and unstable (KS) in rolling.

Table 1. Interface Oreintation Relationships (ND) [RD].

Interface Cu Nb

C-I (121) [111] (121) [101]

D-I (4411) [11118] (121) [101]

Goss-I (101) [010] (121) [101]

Z-I (515) [110 1] (121) [101]

Brass-I (101) [121] (121) [101]

KS (111) [121] (101) [121]

5.1. Plastic Stability Metrics

To compare the plastic stability of the different interfaces, we consider two primary variables: the magnitude of the lattice rotation and its rate of change with respect to the deformation history. In particular, we are interested in how these quantities are influenced by grain interactions. The results are discussed in terms of what can be considered measures of the overall and instantaneous plastic stability of the interface configurations at 50% reduction. The bicrystal simulations serve to establish the baseline (intrinsic) plastic stability of the interfaces to which the multicrystal simulation results will be compared.

The total plastic stability is defined in terms of the magnitude of lattice rotation, AG, in each phase due to the applied deformation, and the angle is determined from the elastic (lattice) rotation tensor, Re, as:

/ trRe _ 1 \

AG = arccos (---j (9)

where Re is the orthogonal rotation tensor associated with the polar decomposition of Fe. The rotation angles are computed at the integration point level and then averaged over the grain to give a single value of AG. It is also informative to characterize the instantaneous stability of a crystal orientation as a function of the applied deformation. The lattice spin rate, We, is defined as:

We = skw (L - Lp) (10)

where L and Lp are total and plastic velocity gradients, respectively, defined on the current configuration. In the following, we take a simpler approach to quantifying the rate of lattice reorientation in terms of only the rotation angle. The lattice spin rate tensor, We, contains information about both the angle and axis of lattice rotation and we found metrics based on the norm of We to be too noisy for our purposes. Instead, we use the rate of change of AG with respect the height reduction, R, as a metric of instantaneous lattice reorientation. The rates, dR, are computed numerically from plots of the average value of AG as a function of R. Figure 2 shows representative plots of AG vs. R and dR vs. R.

Figure 2. Lattice rotation (a) angle; and (b) rate plotted as functions of height reduction for the Cu phase in the C-I interface. In the legend, "s" refers to the periodic bicrystal simulations and the "ml, m2,..." refer to the different realizations of the multicrystal simulations.

0.1 0.2 0.3 0.4

height reduction, R

60 50 40 30 20 10 0 -10 -20

0.0 0.1 0.2 0.3 0.4

height reduction, R

5.2. C-I and D-I Interfaces

The common interfaces in ARB composites lie between the C-I and D-I interfaces in orientation space. We first examine these interfaces using the bicrystal simulations with no grain-grain interaction effects. These simulations show that neither the ideal C-I or D-I interface is perfectly stable. Both interfaces reorient slightly under the applied rolling deformation. At 50% reduction the Cu crystal in the C-I interface rotates 5.7° about the minus TD (towards the D component) to a component the lies

between the C and D components and the Nb crystal rotates 2.04°, also about the minus TD. In contrast, the Cu crystal in the D-I interface rotates in the opposite direction (towards the C component) by an angle of 3.8° whereas the Nb crystal rotates nominally about the minus TD by 5.0°. It is interesting to note that, while the Cu crystal in each interface is rotating toward the other component, e.g., the C component is rotating toward the D component and vice versa, they do not stabilize at the same final orientation. Perhaps more surprisingly, there is some overshoot. Here overshoot means that the Cu crystal initially in the C component orientation reoriented to a position in orientation space beyond the final position of the crystal initially in the D component orientation. However, this overshoot is minor at 50% reduction and the final orientations of the Cu crystals in each case are only separated by 1.5°. The final orientations are relatively stable in that the rates of lattice rotation, as measured by ^dR^, are 0.11° and 3.99° for the Cu and Nb phases, respectively, in the C-I interface and 5.46° and 7.69° for the Cu and Nb phases, respectively, in the D-I interface. Although, the interface plane and OR are strictly lost, the C-I and D-I interfaces are classified as metastable since the crystal orientations have stabilized in close proximity to the initial configurations. Notably, the prediction that the stable interfaces lie between C-I and D-I is in agreement with experimental measurements [18,19].

Figure 3. Combined inverse pole figures at 50% reduction from the multicrystal simulations for the (a) copper phase in the C-I interface; (b) niobium phase in the C-I interface; (c) copper phase in the D-I interface; and (d) niobium phase in the D-I interface. The final orientations obtained from the periodic bicrystal simulations are indicated by the filled black circles. RD [100] /TD [010] /ND [001].

Significantly, we find that the multicrystal simulation results are consistent with the periodic bicrystal results for both the C-I and D-I interfaces. The combined inverse pole figures (IPFs) for each microstructural realization (m1-m5) for the C-I and D-I interfaces are given in Figure 3 where the final orientations from the periodic bicrystal simulations are indicated in the contour plots with filled black circles. The average misorientation angles calculated from the multicrystal simulations are actually smaller for the Cu phase in the C-I interface (2.91°) and for the Nb phase in the D-I interface (3.68°) as

compared to what is observed in the periodic bicrystal simulations. This result is somewhat surprising since it is expected that the additional heterogeneity introduced by grain interactions would serve to decrease the plastic stability. The average lattice rotations for the Nb phase in the C-I interface (3.56°) and the Cu phase in the D-I interface (7.4°) display the anticipated trend in that the total lattice rotations have been amplified due to the deformation heterogeneity. The lattice rotation rates plotted as a function of height reduction are given in Figure 4 for both interfaces. As shown, the multicrystal rotation rates display more pronounced oscillatory behavior compared to the periodic bicrystal simulation results. Notably, this plot also reveals the that the multicrystalline response curves appear to be approximating and/or oscillating about the asymptotic response of the periodic bicrystal simulations.

Figure 4. Lattice reorientation rate plotted as a function of height reduction for the (a) copper phase in the C-I interface; (b) niobium phase in the C-I interface; (c) copper phase in the D-I interface; and (d) niobium phase in the D-I interface. In the legend, "s" refers to the periodic bicrystal simulations and the "ml, m2, ..." refer to the different realizations of the multicrystal simulations.

60 50 40 30 20 10 0 -10 -20

0.0 0.1 0.2 0.3 0.4

height reduction, R

60 50 40 30 20 10 0 -10 -20

0.0 0.1 0.2 0.3 0.4

height reduction, R

Qi •w

40 30 20 10 0

0.0 0.1 0.2 0.3 0.4

height reduction, R

05 •w

0.0 0.1 0.2 0.3 0.4

height reduction, R

5.3. Goss-I and Z-I Interfaces

The Goss-I and Z-I interfaces (see Table 1) are experimentally observed in very fine layered ARB composites [40]. The bicrystal simulations reveal that the Goss-I interface is extraordinarily stable under plane strain compression as lattice rotations in each phase are less than 1°, whereas the Z-I interface is unstable. In fact, the Cu crystal in the Z-I interface rotates 10.17° about the minus TD to a near Goss orientation, while the Nb crystal is stable and rotates 1.89° about the TD. Due to the Goss-I interfaces exceptional stability, the rates of lattice reorientation in both crystals of the Goss-I interface at 50% reduction are negligible. The results also demonstrate that the final orientations of the Z-I interface, which correspond to a near Goss-I interface, are also stable as indicated by the small lattice reorientation rates of 5.37° and 3.55° for Cu and Nb, respectively. To assess the influence of grain interactions, we compare in Figure 5 the inverse pole figures for the Goss-I and Z-I interfaces in the bicrystal and multicrystal configurations. The results show that the final orientations predicted by the multicrystal calculations are consistent with the bicrystal calculations for both interfaces.

Figure 5. Combined inverse pole figures at 50% reduction from the multicrystal simulations for the (a) copper phase in the Goss-I interface; (b) niobium phase in the Goss-I interface; (c) copper phase in the Z-I interface; and (d) niobium phase in the Z-I interface. The final orientations obtained from the periodic bicrystal simulations are indicated by the filled black circles. RD [100] /TD [010] /ND [001].

In the multicrystal calculations, the Goss-I interface remains stable although magnitude of the lattice reorientation in each phase has slightly increased from <1° to average values over the five realizations of 3.21° and 3.85° for Cu and Nb, respectively. As in the bicrystal simulation, the Z-I interface reorients primarily to a Goss-I interface although the IPFs indicate there is a small amount of the {101}(111) component in the Cu phase. The average lattice rotation in the Cu phase of the Z-I interface is roughly the same as in the bicrystal simulation (9.76°), while it has increased slightly in the Nb phase to 3.89°. As in the previous set of results and shown in Figure 6, the evolution of the lattice rotation rates in the

multicrystal simulations oscillate about the bicrystal response curves. Based on these calculations, we conclude that the plastic stability of the Goss-I and Z-I interfaces is not affected by grain interactions.

Figure 6. Lattice reorientation rate plotted as a function of height reduction for the (a) copper phase in the Goss-I interface; (b) niobium phase in the Goss-I interface; (c) copper phase in the Z-I interface; and (d) niobium phase in the Z-I interface. In the legend, "s" refers to the periodic bicrystal simulations and the "m1, m2, ..." refer to the different realizations of the multicrystal simulations.

50 40 30 20 10 0 -10 -20

height reduction, R

50 40 30 20 10 0 -10 -20

0.0 0.1 0.2 0.3 0.4

height reduction, R

5.4. The Brass-I Interface

The Brass-I interface is observed at lower frequencies in Cu-Nb ARB composites as compared to the C-I and D-I interfaces; however, texture measurements suggest that it is present for all submicron to nanolayer thicknesses [10,19]. The Brass-I interface is predicted to be remarkably stable by the bicrystal simulations. Both crystals rotate less than 1° from their original orientation and the lattice reorientation rates are essentially zero throughout the duration of loading. Similarly, the multicrystal simulations reveal that the Brass-I interface is stable given that the average lattice rotations at 50% reduction over the five microstructural realizations are 3.47° and 3.72° for Cu and Nb, respectively. The only notable

difference is revealed by the inverse pole figures shown in Figure 7, which indicate that the peak direction aligned along the RD in the Cu phase has shifted from the (112) direction to a direction in the vicinity of (338). No shift is observed in the IPFs of the Nb phase. The evolution of the lattice rotation rates given in Figure 8 are qualitatively similar to those for the previously considered interfaces. As before, we find that the plastic stability of the Brass-I interface is not significantly affected by grain interactions.

Figure 7. Combined inverse pole figures at 50% reduction from the multicrystal simulations for the (a) copper phase in the Brass-I interface; and (b) niobium phase in the Brass-I interface. The final orientations obtained from the periodic bicrystal simulations are indicated by the filled black circles. RD[100]/TD[010]/ND[001].

(a) (b)

Figure 8. Lattice reorientation rate plotted as a function of height reduction for the (a) copper phase in the Brass-I interface; and (b) niobium phase in the Brass-I interface. In the legend, "s" refers to the periodic bicrystal simulations and the "ml, m2, ..." refer to the different realizations of the multicrystal simulations.

0.0 0.1 0.2 0.3 0.4

height reduction, R

0.0 0.1 0.2 0.3 0.4

height reduction, R

5.5. The KS Interface

Thus far we have examined interfaces that appear in multilayered Cu-Nb composites manufactured by ARB. As previously mentioned, the KS interface is predominant in PVD Cu-Nb composites but is not observed in ARB Cu-Nb composites [18]. In the bicrystal simulation, the KS interface is shown to be highly unstable. At 50% reduction, the Cu crystal has rotated by 25.91° about the TD and the Nb crystal

by 27.63° about the minus TD such that a 53° misorientation has developed between the two crystals with respect to the initial orientations. As might be expected, the lattice reorientation rates for both phases at 50% reduction are significantly higher than those observed for the ARB "x"-I interfaces-19.28° and 36.62°, for Cu and Nb, respectively. This result suggests that the lattice orientations have not stabilized at 50% reduction.

The multicrystal simulations, consistent with the bicrystal results, show that the KS interface is highly unstable. It is noteworthy that the average lattice rotations computed for the multicrystal calculations are approximately the same as those in the bicrystal calculations. In other words, the intrinsic instability of the interface is not further amplified by grain interaction effects, at least in terms of its average behavior. The KS interface is also unique from the "x"-I interfaces studied previously in that grain interaction effects lead to substantially different final orientations as compared to those obtained with the bicrystal model. As indicated by the IPFs, the Cu phase in the bicrystal simulations has reoriented to a texture distributed about the D component and the Nb phase to the I component. In contrast, the Cu phase develops two components in the multicrystal simulation: one is near the predicted bicrystal orientation, but closer to the C orientation, and the other is the Goss component. The Nb phase in the multicrystal calculations is not affected by this additional reorientation behavior in Cu and continues to exhibit a strong I component orientation. Similar to the other interfaces studied, the lattice rotation rate response curves oscillate about the bicrystal response curves for the respective phases. However, in contrast to the other interfaces, the lattice reorientation at 50% reduction does not appear to have stabilized, although one might expect that texture saturation is imminent since the phases have reoriented to stable rolling orientations and/or stable interface configurations, i.e., C-I and Goss-I.

Figure 9. Combined inverse pole figures at 50% reduction from the multicrystal simulations for the (a) copper phase in the KS interface; and (b) niobium phase in the KS interface. The final orientations obtained from the periodic bicrystal simulations are indicated by the filled black circles. RD[100]/TD[010]/ND[001].

6. Discussion

The interfaces studied in the present work were motivated by experimental characterization of the ARB and PVD Cu-Nb lamellar composites. This characterization has shown that single crystal-like textures develop in ARB composites when the layer thicknesses are submicron [10,18,41]. Further, two distinct classes of predominant interfaces develop: (1) the "submicron interface", which lies in the vicinity of the C-I and D-I type of interfaces and is the primary interface prior to the onset of extensive twinning in the Cu phase and (2) the "nano interface", which is near the Goss-I and Z-I interfaces and

becomes the dominant interface after substantial reorientation of the Cu phase due to twinning when the layer thicknesses are < 50 nm [40]. Prior to the onset of twinning, the Goss-I interface is also observed, but not as frequently as the near C-I/D-I interface.

Figure 10. Lattice reorientation rate plotted as a function of height reduction for the (a) copper phase in the KS interface; and (b) niobium phase in the KS interface. In the legend, "s" refers to the periodic bicrystal simulations and the "ml, m2, ..." refer to the different realizations of the multicrystal simulations.

0.0 0.1 0.2 0.3 0.4

height reduction, R

0.0 0.1 0.2 0.3 0.4

height reduction, R

Additionally, we studied the Brass-I interface, which is a secondary type of "submicron interface" in the ARB composites and the KS interface, which is the predominant interface of PVD Cu-Nb lamellar composites. In general, the periodic bicrystal and multicrystal calculations predict similar final orientations for the grains of interest at 50% reduction for the interfaces considered. The orientations from the multicrystal calculations (given by the IPF contour plots) essentially represent distributions centered about the predictions of the bicrystal results. The notable exception is the Cu crystal in the KS interface. As shown in Figure 9a, a Goss component of texture develops in the Cu phase in addition to the near |112}(111) component that is predicted by the bicrystal simulations. This result is significant because it shows that the KS interface reorients under rolling to form two interfaces, the C-I and Goss-I interfaces, which according to our calculations are stable and also prominent in the Cu-Nb lamellar composites.

Mayeur et al. [31] demonstrated that under the kinematic constraint of co-deformation the intrinsic plastic stability (no grain interactions) during rolling is determined by the less plastically stable orientation of the two making up the interface. To quantify the plastic stability of an interface, a simple metric based on the average angle of lattice rotation in each crystal was introduced. The interface stability parameter, Q, is defined as:

A6cu\ { A6ni/

Q = exp —

where A6Cu and AONb are the average angles of the lattice rotation within the Cu and Nb crystals, respectively, at a given strain level e. e is the magnitude of the true compressive strain and is given in terms of the rolling reduction, R, as

|ln (1 — R) |

Additionally, we introduce a rate form of Equation (11) as a measure of the instantaneous plastic stability of the interface at 50% reduction denoted as Q. The instantaneous stability parameter, is defined as in Equation (11) except the rotation angle and strain are replaced with their respective rates. The interface stability parameters for each of the interfaces studied are given in Figure 11 where the filled circles represent the periodic bicrystal result and the open circles represent the results for the five different orientation distributions from the multicrystal simulations.

Figure 11. Plastic stability parameters plotted for each of the considered interfaces (a) total stability; and (b) instantaneous stability at 50% reduction. Filled circles are from the periodic bicrystal simulations and open circles are for the different multicrystal realizations.

_ 1 L.J 1 - - - 1 J....Î 1 l-f 1 H 1 1 _

" , , , , I I

C-l D-l G-l Z-l B-l KS


: ! 1 Ï J r ! 1 ; o ■

O" - - -* : 8 ■

" , , , , , , 1 , ■

C-l D-l G-l Z-l B-l KS


The comparisons made in Figure 11 indicate that grain interactions reduce the total plastic stability in general. The lone exception is the C-I interface, which actually exhibits slightly improved stability for two of the microstructural realizations, although the improvement is modest. Interestingly, the metastable interfaces, i.e., those with 0.8 < Q < 0.9, appear to be less sensitive to extrinsic influences than interfaces that have high (Goss-I, Brass-I) intrinsic stability. This result is expected since one would not anticipate the near-perfect stability to be retained in the presence of heterogeneous deformations induced by grain interactions. It is also interesting to note that the plastic stability of the C-I interface is comparable to that of the Goss-I and Brass-I interfaces when grain interactions are considered. As shown in Figure 11b, all of the "x"-I interfaces have very high values of the instantaneous plastic stability parameter, which implies that orientations given by the IPFs for these interfaces have more or less stabilized at 50% reduction. Furthermore, there is much less scatter in the Q plot between the bicrystal and multicrystal calculations for these interfaces as compared to the Q plot, a result that suggests that extrinsic affects have a more pronounced effect on the magnitude of the lattice rotation than on the rate at which these interfaces establish a stable orientation. The instantaneous interface stability parameter plot also shows that the KS interface has not stabilized at 50% reduction and suggests that the C-I and Goss-I interface textures will continue to sharpen with increased deformation.

The simulation results support the experimental observation that the preferred "submicron interface" lies between the C and D components in the Cu phase and few degrees from the I component in the Nb phase. These interfaces are shown to have moderately high total intrinsic plastic stability that is

not substantially degraded by the extrinsic grain interactions. It is also shown that the Goss-I and Brass-I interfaces, which have near perfect total intrinsic stability parameter values, are quite stable in the multicrystal simulations as well; a result that is also in accord with experimental characterization of ARB Cu-Nb composites. According to the simulations, the Z-I and KS interfaces are unstable and reorient to achieve more favorable texture components. The Z-I interface moves toward a Goss-I interface and the KS interface evolves into a mix of C-I and Goss-I interfaces. The simulations indicate that the Z-I interface observed in the 10 nm ARB composite by Zheng et al. [40] is not plastically stable, wherein the Z component in Cu would diminish at the expense of the Goss component if the nanocomposite were rolled further. The KS interface results agree with the observation of Misra et al. [35] that the interface is unstable at 50% reduction for PVD Cu-Nb multilayered composites with h = 75 f^m.

Recently, it has been proposed that the emergence of predominant interfaces in ARB Cu-Nb multilayer composites can be explained through a combination of two key properties: (1) high mechanical stability under rolling deformation (plastic stability) and (2) low interface formation energy [1]. Due to the large amount of heterophase interfacial area that is created during the ARB process, the texture evolution will be biased toward interfaces that have low formation energy. Beyerlein et al. [1] have used these criteria to explain why the predominant "submicron" and "nano" interfaces emerge at the expense of other stable single phase rolling components.

7. Conclusions

This work presents a computational study of the plastic stability of Cu-Nb heterophase interfaces in rolling and the effects of neighbor grain interactions. The study focuses on interfaces that emerge in ARB Cu-Nb multilayer composites after extreme straining (the layer thickness below 1 f^m). For the six interfaces studied, the simulations show that the final orientations predicted by the multicrystal simulations with neighbor grain interactions are similar to those predicted from a bicrystal simulation with no neighbor grain interactions. In most cases, the bicrystal model yields a uniform reorientation field within each crystal, whereas the multicrystal model predicts a heterogeneous field, which is essentially distributed about the average reorientation from the bicrystal calculation. On this basis, we conclude that co-deformation constraints imposed by the bimetal interfaces are strong, while those by neighboring grain boundaries are relatively weak. As a consequence, the orientational stability conditions of grains at interfaces are altered and thus, different texture components can be expected at interfaces than within the bulk. These findings can help explain the texture gradients observed in coarse-grained multilayered composites and why special preferred orientation distributions emerge at interfaces in nano-grained multilayered composites.


Jason R. Mayeur and Irene J. Beyerlein gratefully acknowledge support by the Center for Materials at Irradiation and Mechanical Extremes, an Energy Frontier Research Center funded by the US Department of Energy, Office of Science, Office of Basic Energy Sciences under Award Number 2008LANL1026. Los Alamos National Laboratory, an affirmative action equal opportunity employer, is operated by Los Alamos National Security, LLC, for the National Nuclear Security Administration of the US

Department of Energy under Contract DE-AC52-06NA25396. Curt A. Bronkhorst and Hashem M. Mourad acknowledge the support of the Los Alamos National Laboratory Directed Research and Development (LDRD) Project DR20110029. Los Alamos National Laboratory is operated by Los Alamos National Security LLC under DOE Contract DE AC52 06NA25396.

Conflicts of Interest

The authors declare no conflicts of interest.


Calibrated constitutive parameters for the single crystal models.

Table 2. Constitutive model parameters.

Parameter Cu Nb Units

P 8960 8400 kg m-3

Cp 380 268 Jkg-i K-i

a 17.7 8.7 ^m m-i K-i

n 0 0 -

mn -36.3 -17.7 MPa K-i

Clio 179.5 242.0 GPa

mi2 -16.4 1.92 MPa K-i

Ci2o 126.4 121.2 GPa

m44 -25.7 2.18 MPa K-i

C44o 82.5 28.0 GPa

r 1.4 1.4 -

Yo 107 107 s-i

so 1 40 MPa

si 20 700 MPa

Fo 1.0 x 10-i8 2.1 x 10-i9 J

P 0.33 0.34 -

q 1.66 1.66 -

ss o 205 140 MPa

ho 200 300 MPa

A 1.5 x 10-i9 1.0 x 10-i7 J


1. Beyerlein, I.J.; Mara, N.A.; Carpenter, J.S.; Nizolek, T.; Mook, W.M.; Wynn, T.A.; McCabe, R.J.; Mayeur, J.R.; Kang, K.; Zheng, S.; et al. Interface-driven microstructure development and ultra high strength of bulk nanostructured Cu-Nb multilayers fabricated by severe plastic deformation. J. Mater. Res. 2013, 28, 1799-1812.

2. Beyerlein, I.J.; Mara, N.A.; Wang, J.; Carpenter, J.S.; Zheng, S.J.; Han, W.Z.; Zhang, R.F.; Kang, K.; Nizolek, T.; Pollock, T.M. Structure-property-functionality of bimetal interfaces. JOM 2012, 64, 1192-1207.

3. Sauvage, X.; Thilly, L.; Lecouturier, F.; Guillet, A.; Blavette, D. FIM and 3D atom probe analysis of Cu/Nb nanocomposite wires. Nanostruct. Mater. 1999, 11, 1031-1039.

4. Montclus, M.A.; Zheng, S.J.; Mayeur, J.R.; Beyerlein, I.J.; Mara, N.A.; Polcar, T.; Llorca, J.; Aldaregula, J.M.M. Optimum high temperature strength of two-dimensional nanocomposites. APL Mater. 2013,1, 052103.

5. Segal, V.M.; Hartwig, K.T.; Goforth, R.E. In situ composites processed by simple shear. Mater. Sci. Eng. A 1997, 224, 107-115.

6. Beyerlein, I.J.; Mara, N.A.; Bhattacharyya, D.; Necker, C.T.; Alexander, D.J. Texture evolution via combined slip and deformation twinning in rolled silver-copper eutectic nanocomposite. Int. J. Plastic. 2011, 27, 121-146.

7. Yasuna, K.; Terauchi, M.; Otsuki, A.; Ishihara, K.N.; Shingu, P.H. Formation of nanoscale Fe/Ag multilayer by repeated press-rolling and its layer thickness dependence of magnetoresistance. Mater. Sci. Eng. A 2000, 285, 412-417.

8. Dupouy, F.; Snoeck, E.; Casanove, M.J.; Roucau, C.; Peyrade, J.P.; Askenazy, S. Microstructural characterization of high strength and high conductivity nanocomposite wires. Scripta Mater. 1996, 34, 1067-1073.

9. Raabe, D.; Ball, J.; Gottstein, G. Rolling textures of Cu-20% Nb composite. Scripta Metall. Mater. 1992, 27,211-216.

10. Carpenter, J.S.; Vogel, S.C.; LeDonne, J.; Hammon, D.L.; Beyerlein, I.J.; Mara, N.A. Bulk texture evolution of Cu-Nb nanolamellar composites during accumulative roll bonding. Acta Mater. 2012, 60, 1576-1586.

11. Carpenter, J.S.; Zheng, S.J.; Zhang, R.F.; Vogel, S.C.; Beyerlein, I.J.; Mara, N.A. Thermal stability of Cu-Nb nanolamellar composites fabricated via accumulative roll bonding. Philos. Mag. 2013, 93, 718-735.

12. Han, W.; Demkowicz, M.J.; Mara, N.A.; Fu, E.; Sinha, S.; Rollett, A.D.; Wang, Y.; Carpenter, J.S.; Beyerlein, I.J.; Misra, A. Design of radiation tolerant materials via interface engineering. Adv. Mater. 2013, 25, 6975-6979.

13. Han, W.Z.; Cerreta, E.K.; Mara, N.A.; Beyerlein, I.J.; Carpenter, J.S.; Zheng, S.J.; Trujillo, C.P.; Dickerson, P.O.; Misra, A. Deformation and failure of shocked bulk Cu-Nb nanolaminates. Acta Mater. 2014, 63, 150-161.

14. Beyerlein, I.J.; Caro, J.A.; Demkowicz, M.J.; Mara, N.A.; Misra, A.; Uberuaga, B.P. Radiation-damage-tolerant nanomaterials. Mater. Today 2013 16, 443-449.

15. Beyerlein, I.J.; Wang, J.; Zhang, R. Mapping dislocation nucleation behavior from bimetal interfaces. Acta Mater. 2013, 61,7488-4799.

16. Beyerlein, I.J.; Wang, J.; Kang, K.; Zheng, S.J.; Mara, N.A. Twinnability of bimetal interfaces in nanostructured composites. Mater. Res. Lett. 2013, 1, 89-95.

17. Zhang, R.F.; Germann, T.C.; Wang, J.; Liu, X.Y.; Beyerlein, I.J. Role of interface structure on the plastic response of Cu/Nb nanolaminates under shock compression: Non-equilibrium molecular dynamics simulations. Scripta Mater. 2013, 68, 114-117.

18. Lee, S.B.; LeDonne, J.E.; Lim, S.C.V.; Beyerlein, I.J.; Rollet, A.D. The five-parameter Heterophase Interface Character Distribution (HICD) of physical vapor-deposited and accumulative roll-bonded Cu-Nb multilayer composites. Acta Mater. 2012, 60, 1747-1761.

19. Carpenter, J.S.; McCabe, R.J.; Beyerlein, I.J.; Wynn, T.A.; Mara, N.A. A wedge-mounting technique for nanoscale electron backscatter diffraction. J. Appl. Phys. 2013, 113, 094304.

20. Carpenter, J.S.; McCabe, R.J.; Vogel, S.C.; Rollet, A.D.; Mara, N.A.; Beyerlein, I.J. Process parameter influence on texture evolution in Cu-Nb multilayer composites fabricated via accumulative roll bonding. Metall. Mater. Trans. A 2014, in press.

21. Wagner, P.; Engler, O.; Lücke, K. Formation of Cu-type shear bands and their influence on deformation and texture of rolled f.c.c. |112}(111) single crystals. Acta Metall. Mater. 1995, 43, 3799-3812.

22. Leffers, T.; Ray, R.K. The brass-type texture and its deviation from the copper-type texture. Progr. Mater. Sci. 2009, 54, 351-396.

23. Toth, L.S.; Jonas, J.J.; Daniel, D.; Ray, R.K. Development of ferrite rolling textures in low- and extra low-carbon steels. Metall. Mater. Trans. A 1990, 21A, 2985-3000.

24. Toth, L.S.; Molinari, A.; Raabe, D. Modeling of rolling texture development in a ferritic chromium steel. Metall. Mater. Trans. A 1997, 28A, 2343-2351.

25. Raphanel, J.; Van Houtte, P. The rolling textures of b.c.c. metals. Acta Metall. 1985, 33, 1481-1488.

26. Hirsch, J.; Lücke, K. Mechanism of deformation and development of rolling textures in polycrystalline f.c.c. metals—I. Description of rolling texture development in homogeneous CuZn alloys. Acta Metall. 1988, 36, 2863-2882.

27. Zhou, Y.; Toth, L.S.; Neale, K.W. On the stability of the ideal orientations of rolling textures for f.c.c. polycrystals. Acta Metall. Mater. 1992, doi:10.1016/0956-7151(92)90481-S.

28. Anderson, P.M.; Bingert, J.F.; Misra, A.; Hirth, J.P. Rolling textures in nanoscale Cu/Nb multilayers. Acta Mater. 2003, 51, 6059-6075.

29. Al-fadhalah, K.; Tome, C.N.; Beaudoin, A.J.; Roberston, I.M.; Hirth, J.P.; Misra, A. Modeling texture evolution during rolling of a Cu-Nb multilayered system. Philos. Mag. A 2005, 85, 1419-1440.

30. Hansen, B.; Carpenter, J.; Sintay, S.; Bronkhorst, C.; McCabe, R.; Mayeur, J.; Mourad, H.; Beyerlein, I.; Mara, N.; Chen, S.; et al. Modeling the texture evolution of Cu/Nb layered composites during rolling. Int. J. Plast. 2013, 49, 71 - 84.

31. Mayeur, J.R.; Beyerlein, I.J.; Bronkhorst, C.A.; Mourad, H.M.; Hansen, B.L. A crystal plasticity study of heterophase interface character stability of Cu/Nb bicrystals. Int. J. Plast. 2013, 48, 72-91.

32. Jia, N.; Roters, F.; Eisenlohr, P.; Raabe, D.; Zhao, X. Simulation of shear banding in heterophase co-deformation: Example of plane strain compressed Cu-Ag and Cu-Nb metal matrix composites. Acta Mater. 2013, 61, 4591-4606.

33. Knezevic, M.; Beyerlein, I.J.; Nizolek, T.; Mara, N.A.; Pollock, T.M. Anomalous basal slip activity in zirconium under high-strain deformation. Mater. Res. Lett. 2013,1, 1-8.

34. Shen, Y.; Anderson, P.M. Transmission of a screw dislocation across a coherent, slipping interface. Acta Mater. 2006, 54, 3941-3951.

35. Misra, A.; Hirth, J.P.; Hoagland, R.G.; Embury, J.D.; Kung, H. Dislocation mechanisms and symmetric slip in rolled nano-scale metallic multilayers. Acta Mater. 2004, 52, 2387-2394.

36. Bronkhorst, C.A.; Mayeur, J.R.; Beyerlein, I.J.; Mourad, H.M.; Hansen, B.L.; Mara, N.A.; Carpenter, J.S.; McCabe, R.J.; Sintay, S.D. Meso-scale modeling the orientation and interface stability of Cu/Nb-layered composites by rolling. JOM 2013, 65, 431-442.

37. Bronkhorst, C.A.; Hansen, B.L.; Cerreta, E.K.; Bingert, J.F. Modeling the microstructural evolution of metallic polycrystalline materials under localization conditions. J. Mech. Phys. Solids 2007, 55, 2351-2383.

38. Wang, Z.Q.; Beyerlein, I.J. An atomistically-informed dislocation dynamics model for the plastic anisotropy and tension-compression asymmetry of BCC metals. Int. J. Plast. 2011, 27, 1471-1484.

39. Van der Sluis, O.; Schreurs, P.J.G.; Brekelmans, W.A.M.; Meijer, H.E.H. Overall behaviour of heterogeneous elastoviscoplastic materials: Effect of microstructural modelling. Mech. Mater. 2000, 32, 449-462.

40. Zheng, S.; Beyerlein, I.J.; Carpenter, J.S.; Kang, K.; Wang, J.; Han, W.; Mara, N.A. High-strength and thermally stable bulk nanolayered composites due to twin-induced interfaces. Nat. Commun. 2013, 4, 1696.

41. Carpenter, J.S.; Liu, X.; Darbal, A.; Nuhfer, N.T.; McCabe, R.J.; Vogel, S.C.; LeDonne, J.E.; Rollet, A.D.; Barmak, K.; Beyerlein, I.J.; et al. A comparison of texture results obtained using precession electron diffraction and neutron diffraction methods at diminishing length scales in ordered bimetallic nanolamellar composites. Scripta Mater. 2012, 67, 336-339.

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (

Copyright of Materials (1996-1944) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.