Abstract

A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

1 Introduction

Material systems comprised of zirconium alloys are used as cladding protection against the release of fission gases and other products during the operation of pressurized water reactors (PWR). This is mainly due to its favorable thermomechanical behavior, ductility, and the insignificant neutron absorption of these zirconium alloys [1]. The most commonly used alloys, for these harsh thermal environments, are Zircaloy-2, Zircaloy-4, and Zirlo (zirconium low oxidation).

However, due to thermal gradients across the cladding and the corrosion of zirconium alloys due to reaction with the coolant water, hydrogen can be transported to the zirconium alloys, and this can result in the precipitation of hydrides at the edge of zirconium-oxide layers. These hydrides generally have crystalline structures that are different from the zirconium alloy h.c.p. matrix and have an atomic structure of zirconium and hydrogen. These hydrides can initially form with circumferential orientations, but due to cycles of stress and thermal loading conditions, these hydrides can reorient to mostly radial distributions.

The orientation and distribution of hydride populations can significantly affect the ductile to brittle transition temperature (DBTT) of the parent zirconium alloy. The DBTT of hydrided zirconium alloys has been reported to range from 50 °C to 260 °C [27]. Billone et al. [8] conducted ring compression tests (RCT) with high burn-up Zircaloy-4 and Zirlo cladding specimens for temperatures ranging between 25 °C and 200 °C that indicated the transition temperatures significantly affected the hydride morphology, orientation, and crystallography. Hydride reorientation in the Zircaloy-4 can also take place during long-term storage, in the absence of a corrosive environment, and this can also significantly reduce the ductility of the cladding system [911]. Kim et al. [1215] experimentally investigated how hydrides affected the cladding mechanical behavior for different conditions, such as those thermal conditions associated with the loss of coolant accident (LOCA). Their results show that hydrides in the circumferential direction at room temperature are generally brittle with some plasticity, but can be ductile above 150 °C. In contrast, for radial hydrides, at 150 °C, micro-cracks were nucleated and as they can act as paths for crack propagation between the different crystalline structures of the matrix and the hydrides.

Chan et al. [16] conducted micro-fracture experiments to determine the local mechanical properties of multi-crystalline systems with a zirconium/zirconium hydride microstructure. Their experiments clearly indicate that different interphase crystalline structures between the δ, ZrH1.66 (f.c.c.) hydrides, the ɛ, ZrH2, and the h.c.p. matrix affect crack nucleation and propagation. How these different interfaces and interphases collectively affect fracture for the hydride population is not clearly understood due to the complex crystalline interactions between the matrix and hydride populations. For example, there are interfacial dislocations and misfit strains, due to semicoherent interfaces between the matrix and the hydrides. [17]. Due to the different lattice parameters and atomic arrangements between the parent and hydride phases, interfacial planes can form and result in heterogeneous crystalline behavior. This is because interfaces, between the closed packed planes of the matrix and the hydrides, have misorientation dislocation and resulting misfit strains. Hence, orientations relations (ORs) have to be physically represented to account for both coherent and semicoherent interfaces [18,19].

Hence, our main goal is to postulate a microstructurally based validated approach that can accurately predict and understand the thermomechanical failure behavior of the hydrided microstructure. This approach is predicated on an accurate physical representation of distributions of circumferential and radial precipitated hydride morphologies, which have distinct crystalline systems than the matrix, such that a fundamental understanding of how multiple fracture failure modes and surfaces nucleate and propagate because of hydride orientations, morphologies, and distributions and matrix interactions with a multiphase system that is comprised of h.c.p. matrix and f.c.c. and b.c.c. hydrides. The multiple-slip crystal plasticity dislocation-density-based formulation and the ORs are briefly outlined in Sec. 2, the microstructural fracture criterion is discussed in Sec. 3, the numerical implementation of the fracture approach is given in Sec. 4, the overall computational approach is given in Sec. 5, the results are presented and discussed in Sec. 6, and a summary of the results and conclusions are given in Sec. 7.

2 Dislocation-Density Evolution and Crystalline Behavior

A dislocation-density crystalline approach introduced by Shanthraj and Zikry [13] will be used to represent the evolution of dislocation densities for the different crystalline structures of the hydrided system. It is assumed in that formulation that the dislocation densities can be decomposed into mobile and immobile dislocation densities. Coupled evolution equations can then be obtained for the mobile and immobile densities as
(1)
(2)
where gsour pertains to dislocation sources, gmnter corresponds to mobile dislocation trapping and forest intersections, obstacle cross-slip, or dislocation interactions, grecov represents the immobile dislocation rearrangement and annihilation, and gimmob are coefficients related to the mobile dislocation immobilization.

The ORs are a critical microstructural characteristic that ensures different crystalline interphases are accurately represented [20]. The ORs for the f.c.c. and b.c.c. hydrides are based on the formulation developed by Ref. [21], and in this investigation, we used an OR of (0001)h.c.p.//(111)f.c.c. for the f.c.c. δ hydrides and an OR of (0001)h.c.p.//(110)f.c.c. for the b.c.c. ɛ hydrides. There are numerous experimentally observed ORs, but as discussed in Ref. [21], these are the most commonly observed experimental ORs for hydrided zirconium alloys.

3 Microstructural Fracture Criterion

The inherent fracture mode in h.c.p., b.c.c., and f.c.c. materials would be cleavage fracture. This would correspond to the (100) planes for b.c.c. and f.c.c. systems and the (1000) basal planes for the h.c.p matrix [22]. Hence, hydride distributions have to be oriented with respect to the global axis and cleavage planes have to be oriented with respect to the matrix. The first needed transformation is based on using the ORs, [T]1, and then, a second transformation from the parent matrix to the global coordinate, [T]2 is needed. The transformations are
(3)
where ncleave is normal to the cleavage plane with the highest fracture stress. It should be noted that the hydrides, both the δ-f.c.c. or ɛ-b.c.c. hydrides, will fracture on the (100) planes, which will coincide with the (0001) of the parent matrix based on the ORs. The global orientation of the cleavage planes, ni, in the current configuration, is then obtained by updating, at every time-step, the lattice rotations as
(4)
where Wij* is the elastic spin. The resulting traction stresses are the components that directly affect crack nucleation [23,24]. Hence, maximum values are monitored and are monitored with a critical fracture stress, σfracture, to initiate failure surfaces. The failure criterion is then given by
(5)
where tcleave=max(0001),planesncleaveT[σ]ncleave, t is the traction vector, such that all components of the stress tensor, σij, are ascertained for the maximum fracture stress, σfracture.

4 Numerical Implementation of the Overlapping Element Method

By splitting an element into two halves along the crack line as first proposed by Hansbo and Hansbo [25] and detailed by Wu and Zikry [26], cracks in the mesh can be obtained. The advantage of the overlap method is that it overcomes the limitations of methods, such as cohesive fracture and element enrichment. As noted by Wu and Zikry [26], the cohesive fracture can introduce compliant surfaces that can affect numerical stability. Enrichment methods are limited by elastic fracture solutions that may not be applicable for inelastic failure modes.

This fracture model is basically implemented by adding phantom nodes to overlap elements that represent discontinuous fracture surfaces [27]. A cracked element is replaced by two overlapping elements, and phantom nodes are introduced over the existing nodes. For each overlapping element, only the subdomain with area Ae1 or Ae2 (Fig. 1), corresponding to one of the two subdomains for the original cracked element, is considered as active. The internal nodal force vector of the cracked element is obtained and is given by
(6)
Fig. 1
Overlapping fracture method: top showing an element cut along the crack line to generate two overlapping elements and bottom showing a crack propagating in a mesh
Fig. 1
Overlapping fracture method: top showing an element cut along the crack line to generate two overlapping elements and bottom showing a crack propagating in a mesh
Close modal

5 Computational Approach

An approach that has been developed specifically for microstructural fracture and dislocation-density crystalline plasticity will be utilized in this investigation. Only a brief outline will be given here, the methodologies are detailed in Refs. [13,28]. An finite element (FE) method with Broyden–Fletcher–Goldfarb–Shanno (BFGS) iteration and one-point integration with hour-glass control is used to obtain the total deformation rate tensor. The dislocation-density evolution equations and the plastic deformation rate tensor is updated by an implicit–explicit fifth-order Runge–Kutta method that accounts for numerical stiffness that arises due to the field quantities associated with the different crystalline structures evolving at different rates. A Crank-Nicholson approach was used to update the local temperatures with both thermal conductivity and internal heat generation due to dissipated plastic.

6 Results and Discussion

The dislocation-density crystalline representation, the fracture nucleation framework, and the Voronoi scheme for hydride formation with OR interfaces were used to investigate the behavior of hydrided zirconium systems. The multiphase aggregate was assumed to have random low-angle misorientations, where the maximum misorientations were assumed to be less than 5 deg. Voronoi tessellations were used with different Euler angles to represent these misorientations, and plane-strain boundary conditions were used to represent cladding behavior.

Seven thousand plane-strain four-noded quadrilateral elements (Fig. 2) were used for an aggregate size of 200 μm by 200 μm, and the mesh was optimized based on a convergence analysis for different grain misorientations. For tensile plane-strain behavior, a quasi-static nominal strain rate of 200 × 10−1 s−1 was applied. The material properties that are used here were obtained for experimental measurements of Zircaloy-4, and δ and ɛ hydrides [6,7], and are listed in Table 1. The main differences between the hydrided aggregates are the b.c.c. and f.c.c. hydrides, which result in distinct fracture behavior due to the evolution of the different dislocation densities and slip-systems.

Fig. 2
200 μm × 200 μm aggregate with the grain and hydride morphologies
Fig. 2
200 μm × 200 μm aggregate with the grain and hydride morphologies
Close modal
Table 1

Properties of Zircaloy-4, δ/ɛ-hydrides

PropertyZircaloy-4δ, ZrH1.66, Hydridesɛ, ZrH2, Hydrides
Young’s modulus (GPa)90131131
Static yield stress (MPa)600700700
Poisson’s ratio0.349
Rate sensitivity coefficient45
Initial immobile dislocation density (m−2)1 × 1010
Initial mobile dislocation density (m−2)1 × 107
Burger’s vector1 × 10−10
Fraction of plastic energy to heat0.9
Thermal softening exponent0.5
Fracture stress (MPa) (N/mm2)840800800
Misorientation angle between grains and hydrides (deg)0–5
Applied tensile nominal strain, %12
Strain loading rate, s−120
PropertyZircaloy-4δ, ZrH1.66, Hydridesɛ, ZrH2, Hydrides
Young’s modulus (GPa)90131131
Static yield stress (MPa)600700700
Poisson’s ratio0.349
Rate sensitivity coefficient45
Initial immobile dislocation density (m−2)1 × 1010
Initial mobile dislocation density (m−2)1 × 107
Burger’s vector1 × 10−10
Fraction of plastic energy to heat0.9
Thermal softening exponent0.5
Fracture stress (MPa) (N/mm2)840800800
Misorientation angle between grains and hydrides (deg)0–5
Applied tensile nominal strain, %12
Strain loading rate, s−120

To fully understand and predict how δ and ɛ hydrides affect fracture nucleation and propagation, we investigated aggregates, for each crystalline hydride, with radial, circumferential, and mixed distributions with different HVF. We refer to distributions as mixed as those with high HVFs that result in aggregates that are then neither fully circumferential or radial. Furthermore, to understand the effects of temperature on deformation and local behavior, we investigated aggregate behavior with radial hydrides for three different temperatures of 20 °C, 250 °C, and 500 °C.

6.1 The Nucleation and Propagation of Multiple Cracks.

The microstructural fracture behavior of hydride zircaloy-4 for radial and circumferential hydride distributions, with different hydride crystallographic structures and different hydride volume fractions (HVFs), were investigated. The failure criterion is based on the normal components of the stress traction exceeding the experimentally obtained fracture stresses given in Table 1. The fracture formulation also accounts for the crystallographic structures of ORs that are uniquely inherent to f.c.c. and b.c.c. hydrides and the parent h.c.p. matrix microstructures [21].

We investigated three different fracture stresses to further understand the response of the Zircaloy-4 system without hydrides (Fig. 3), and these cases correspond to different fracture stresses of 720, 840, and 960 MPa, and it confirms the sensitivity of the matrix to different fracture stresses. Namely that as the fracture stress increases, the fracture is inhibited as a function of the applied strain, which is physically realistic. Based on this, we investigated the following cases to obtain a more detailed representative understanding of fracture:

  • Different morphologies with different hydride volume fractions.

  • Radial and circumferential distributions.

  • Hydrides with different crystalline structures: δ f.c.c. and ɛ b.c.c. hydrides.

Fig. 3
Global response for different applied fracture stresses
Fig. 3
Global response for different applied fracture stresses
Close modal

6.2 Different Morphologies With Different Hydrides Volume Fraction.

To better understand the mechanistic effects of HVF on the fracture behavior, we analyzed three different models with different HVFs for a fracture stress of 840 MPa and at a room temperature of 20 °C for a radial distribution of δ f.c.c. hydrides. As shown in the stress–strain response (Fig. 4), as the HVF value increases, the fracture is inhibited. For the higher HVFs, there are significant increases in plasticity, as indicated by global behavior beyond the apparent yield point, before the global unloading of the stresses. This can be further understood from the contours for the different values of the HVF, which are given for the nominal strains that correspond to the final point of unloading. The contours are for the normalized (normalized by the static yield stress) fracture stress, and the circled regions correspond to the different crack nucleation sites. All the contours correspond to nominal strains of 1.8%. As these results show, the aggregate with the higher HVF ruptures at a much higher nominal strain than the two other cases (Fig. 5) and much higher fracture stress is needed in comparison with the two other cases. This can be attributed to the significant plasticity associated with the higher nominal strains of the higher HVFs. Plasticity ahead of the crack fronts inhibits crack propagation, which requires higher fracture stresses to rupture the aggregate. This is physically consistent with fracture observations of δ radial hydrides within zirconium alloys [16].

Fig. 4
Global response for different hydrides volume fractions using fracture stress of 840 MPa
Fig. 4
Global response for different hydrides volume fractions using fracture stress of 840 MPa
Close modal
Fig. 5
Comparison between different HVFs for crack nucleation and crack propagation, circles show crack initiation sites at a nominal strain of 1.8%
Fig. 5
Comparison between different HVFs for crack nucleation and crack propagation, circles show crack initiation sites at a nominal strain of 1.8%
Close modal

6.3 Predictions for Different Hydride Distributions.

We wanted to further understand the effects of hydride orientation on this local ductile fracture behavior, so we investigated the behavior of δ-hydrides in the radial and circumferential directions with the highest fracture stress of 960 MPa. Figure 6 shows the global nominal stress–strain response for the two models at an HVF of 7%. For the radial hydrides, there is rapid unloading, in comparison with the circumferential case. For the circumferential hydride case, there is more ductility, which would delay crack rupture as the predicted results indicate. Billone et al. [8] reported that radial hydrides promote unstable crack propagation, whereas circumferential hydrides tend to promote plasticity and/or stable crack growth, so the predicted results are consistent with these experimental results. The results also showed that the hydrides act as initiation sites for crack nucleation (Fig. 7) for both the radial and circumferential hydrides. There are multiple crack nucleation sites, but not all of them coalesce to form bigger cracks. All the contours correspond to nominal strains of 2.2%. As a single crack develops, others could be blocked by potential pile-ups and accumulations of dislocation densities (Fig. 8).

Fig. 6
Nominal stress–strain response for radial and circumferential hydrides
Fig. 6
Nominal stress–strain response for radial and circumferential hydrides
Close modal
Fig. 7
Crack arrest due to extensive immobile dislocation-density accumulation for the case of δ radial f.c.c. hydrides. All the contours correspond to nominal strains of 2.2%.
Fig. 7
Crack arrest due to extensive immobile dislocation-density accumulation for the case of δ radial f.c.c. hydrides. All the contours correspond to nominal strains of 2.2%.
Close modal
Fig. 8
Comparison between radial and circumferential hydrides for crack nucleation and crack propagation. The circles denote crack initiation sites.
Fig. 8
Comparison between radial and circumferential hydrides for crack nucleation and crack propagation. The circles denote crack initiation sites.
Close modal

6.4 Comparison of Aggregates With δ and ɛ Hydrides.

We investigated two radial hydride crystallographies: δ-hydrides, which are f.c.c., and ɛ hydrides, which are b.c.c., as outlined in Table 1. We examined the effects of these two types of hydrides, which span different fracture stresses (Table 1) for an HVF of 17%. At this high HVF, the contours indicate that the hydride distribution can be considered to be mixed: neither circumferential nor radial due to the high volume fraction. Figure 9 shows the nominal stress–strain response for the two hydrides types. As these results indicate, the aggregate with the δ f.c.c. hydrides unloads gradually in comparison with the ɛ b.c.c. aggregate. This gradual unloading is due to the higher plasticity associated with the δ f.c.c. aggregate in comparison with the b.c.c. aggregate, and this can be due to the different ORs associated with the h.c.p./b.c.c. and the h.c.p/f.c.c. systems as discussed by Ref. [29]. Figure 10 shows the normalized fracture stresses for the two hydride populations. As these predictions indicate, and based on the global nominal stress strain, the gradual unloading of the aggregates with the delta hydrides results in rupture at a higher nominal strain and a higher normalized fracture stress (Fig. 10). The contours correspond to a nominal strain of 2.8%. The higher plasticity associated with the δ hydrides results in the higher fracture stresses. This is consistent with experimental observations of how the different microstructures can affect fracture [30].

Fig. 9
Nominal stress–strain response for the h.c.p. matrix with δ (f.c.c.) and ɛ (b.c.c.) hydrides.
Fig. 9
Nominal stress–strain response for the h.c.p. matrix with δ (f.c.c.) and ɛ (b.c.c.) hydrides.
Close modal
Fig. 10
Comparison between radial and circumferential hydrides for crack nucleation and crack propagation. The circles denote crack initiation sites at a nominal strain of 2.8%.
Fig. 10
Comparison between radial and circumferential hydrides for crack nucleation and crack propagation. The circles denote crack initiation sites at a nominal strain of 2.8%.
Close modal

6.5 Thermomechanical Effects of Hydrides on the Deformation of Zircaloy.

We next investigated the effects of the thermomechanical loading on the behavior of Zircaloy-4 deformation and cracking. The coupling between the heat conduction and the mechanical loading was done at each time-step. Here, we will focus on radial hydrides and present the effects of thermomechanical loading on ductility, crack nucleation, and propagation.

To understand the effects of temperature on deformation and local behavior, we subjected the aggregate with radial hydrides to three different temperatures of 20, 250, and 500 °C for an HVF of 10%, and we inhibited fracture to focus on ductility and how deformation modes evolve with increasing temperature. Figure 11 shows the global nominal stress–strain response for the three temperatures. As the temperature increases, the response is significantly more ductile. Furthermore, as the applied temperatures increase, the normalized (normalized by the static yield stress) decreases as shown in Fig. 12. This has also been observed by Kim et al. [15,31]. This ductility is due to the increased normalized immobile dislocation density as the temperature increases (Fig. 13). These dislocation-density systems are the most active normalized dislocation density for the pyramidal system, as shown.

Fig. 11
Stress–strain response for model validation
Fig. 11
Stress–strain response for model validation
Close modal
Fig. 12
Normalized normal stresses for different temperatures at 12% applied tensile strain. The stresses are normalized by the static yield stress.
Fig. 12
Normalized normal stresses for different temperatures at 12% applied tensile strain. The stresses are normalized by the static yield stress.
Close modal
Fig. 13
Hydride deformation for different temperatures for normalized dislocation-density pyramidal system (011¯1)[1¯21¯3]hcp at 12% applied tensile strain
Fig. 13
Hydride deformation for different temperatures for normalized dislocation-density pyramidal system (011¯1)[1¯21¯3]hcp at 12% applied tensile strain
Close modal

To further understand the effects of thermal softening due to high temperature on the dislocation density, a section plot is shown in Fig. 14 for a temperature of 500 °C. This line was taken along a hydride and shows the evolution of the immobile dislocation density along the x-direction. The plot shows that the dislocation-density values for both the hydrides and the parent matrix at a nominal strain of 12%. The dislocation densities for both f.c.c. hydrides and the h.c.p. matrix are almost of the same magnitude at this high temperature, which further underscores the thermal softening associated with this high temperature.

Fig. 14
Section plot of immobile dislocation density in the radial hydrides and the parent matrix at 500 °C and at 12% applied tensile strain
Fig. 14
Section plot of immobile dislocation density in the radial hydrides and the parent matrix at 500 °C and at 12% applied tensile strain
Close modal

6.6 Thermomechanical Effects on Fracture and Crack Nucleation.

To investigate the effect of temperature on fracture, we investigated two aggregates with two different temperatures of 20 °C and 500 °C. Figure 15 shows the global nominal stress–strain response for the two aggregates. As the temperature increases, the aggregates become more ductile. As noted earlier, this was also observed by Kim et al. [15,31] for δ hydrides. At room temperature, the aggregate has little ductility in comparison with the aggregate subjected to higher temperatures. Also, at room temperature, the specimens were fractured at the early stages of deformation and breached even before reaching the yield point. Similar behavior is also obtained for the circumferential hydrides (Fig. 16) at an HVF of 10%.

Fig. 15
Nominal stress–strain response for radial hydrides at room temperature and at 500 °C under the same fracture stress of 840 N/mm2 [15]
Fig. 15
Nominal stress–strain response for radial hydrides at room temperature and at 500 °C under the same fracture stress of 840 N/mm2 [15]
Close modal
Fig. 16
Nominal stress–strain response for circumferential hydrides at room temperature and at 500 °C under the same fracture stress of 840 N/mm2 [15]
Fig. 16
Nominal stress–strain response for circumferential hydrides at room temperature and at 500 °C under the same fracture stress of 840 N/mm2 [15]
Close modal

The effects of temperature on shear strain accumulation can be clearly seen in Figs. 17 and 18 for both radial and circumferential distributions. As the contours show, at room temperature, for both radial and circumferential distributions, rupture is significantly inhibited as the temperature increases. As the predictions indicate, at a nominal strain of 2.5% for both thermal cases, the shear strain for both radial and circumferential distributions, the shear strain at 500 °C is almost four times larger than what it is at 20 °C. This increased dislocation-density activity and shear strain accumulation are the mechanistic effects that inhibit fracture at high temperatures for these hydride populations and orientations. These high temperatures lead to thermal softening, shear localization, and ductile failure [2,32].

Fig. 17
Accumulated shear strain for radial hydrides for crack nucleation and crack propagation at different temperatures: 20 °C and 500 °C. Circles denote crack initiation sites and the nominal strain is 2.5%.
Fig. 17
Accumulated shear strain for radial hydrides for crack nucleation and crack propagation at different temperatures: 20 °C and 500 °C. Circles denote crack initiation sites and the nominal strain is 2.5%.
Close modal
Fig. 18
Accumulated shear strain for circumferential hydrides crack nucleation and crack propagation at different temperatures: 20 °C and 500 °C. Circles denote crack initiation sites.
Fig. 18
Accumulated shear strain for circumferential hydrides crack nucleation and crack propagation at different temperatures: 20 °C and 500 °C. Circles denote crack initiation sites.
Close modal

7 Summary

A fracture methodology specialized for simultaneous crack surface nucleation was implemented within a nonlinear finite-element framework to accurately represent a microstructural multiphase system with f.c.c. δ and b.c.c. ɛ hydrides within an h.c.p. Zircaloy-4 matrix. This validated predictive framework provides an understanding of fundamental microstructural mechanisms and characteristics that span the physical scales of crystalline structures of the hydrides, the parent matrix, and the ORs between the h.c.p. matrix and the f.c.c. and b.c.c. hydrides.

Radial, circumferential, and mixed hydride orientations were investigated to understand how inelastic mechanistic behavior of these multiphase systems evolves for radial, circumferential, and mixed orientations for different HVFs. Circumferential hydride distributions had higher strength and ductility than the radial distributions. A key aspect for both distributions is that the hydrides do contribute to ductile behavior due to dislocation-density accumulations at hydride-matrix interphases. As the HVF was increased beyond a volume fraction of 7%, the hydride distributions were mixed, in the sense that radial and circumferential hydrides could no longer be distinguished. As the HVF increased, the stresses and dislocation densities increased, which further indicates that hydrides both can strengthen the hydrided system and render it more ductile.

The hydrided ɛ b.c.c. system generally had more strength and ductility than the δ f.c.c. system, and this was mainly because of the availability of more slip-systems associated with the b.c.c. ɛ hydrides. This has resulted in the matrix with the b.c.c. ɛ hydrides having fracture modes that occurred at higher nominal strains than the matrix with the f.c.c. δ hydrides. A key aspect of the predictions was the incorporation of temperature changes that accounted for both thermal conductivity and the dissipation of plastic work. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and propagation. The present analysis indicates how underlying thermomechanical mechanisms collectively related to slip-systems, dislocation densities, local stresses, volume fractions, phase boundaries, and interfaces between crystalline systems affect and dominate local microstructural behavior. Therefore, it is essential to account for and identify these dominant microstructural mechanisms to understand how multiple fracture surfaces nucleate and propagate in these multiphase crystalline systems.

Acknowledgment

Supports from the Consortium for Advanced Simulation of Light (CASL) Water Reactors an Energy Innovation Hub for Modeling and Simulation of Nuclear Reactors, U.S. Department of Energy, Contract No. DE-AC05-00OR227 and the DOE NEUP Integrated Research Project IRP-DOE 17-13708, Development of a Mechanistic Hydride Behavior Model for Spent Fuel Cladding Storage and Transportation are gratefully acknowledged.

Conflict of Interest

There are no conflicts of interest.

References

1.
Northwood
,
D. O.
,
1985
, “
The Development and Applications of Zirconium Alloys
,”
Mater. Desi.
,
6
(
2
), pp.
58
70
.
2.
Mishima
,
Y.
, and
Okubo
,
T.
,
1972
, “
Effect of Thermal Cycling on the Stress Orientation and Circumferential Ductility in Zircaloy-2
,”
Can. Metall. Q.
,
11
(
1
), pp.
157
164
.
3.
Wallace
,
A.
,
Shek
,
G.
, and
Lepik
,
O.
,
2008
, “
Effects of Hydride Morphology on Zr-2.5Nb Fracture Toughness
,”
Zirconium in the Nuclear Industry: Eighth International Symposium
,
ASTEM Special Technical Publications, ASTM International
.
4.
Fearnehough
,
G. D.
, and
Cowan
,
A.
,
1967
, “
The Effect of Hydrogen and Strain Rate on the ‘Ductile-Brittle’ Behaviour of Zircaloy
,”
J. Nucl. Mater.
,
22
(
2
), pp.
137
147
.
5.
Mouton
,
I.
,
Breen
,
A. J.
,
Wang
,
S.
,
Chang
,
Y.
,
Szczepaniaka
,
A.
,
Kontis
,
P.
,
Stephensona
,
L. T.
,
Raabe
,
D.
,
Herbig
,
M.
,
Britton
,
T. B.
, and
Gaulta
,
B.
,
2019
, “
Quantification Challenges for Atom Probe Tomography of Hydrogen and Deuterium in Zircaloy-4
,”
Micros. Microanal.
,
25
(
2
), pp.
481
488
.
6.
Bradbrook
,
J. S.
,
Lorimer
,
G. W.
, and
Ridley
,
N.
,
1972
, “
The Precipitation of Zirconium Hydride in Zirconium and Zircaloy-2
,”
J. Nucl. Mater.
,
42
(
2
), pp.
142
160
.
7.
Hong
,
S. I.
,
Lee
,
K. W.
, and
Kim
,
K. T.
,
2002
, “
Effect of the Circumferential Hydrides on the Deformation and Fracture of Zircaloy Cladding Tubes
,”
J. Nucl. Mater.
,
303
(
2–3
), pp.
169
176
.
8.
Billone
,
M. C.
,
Burtseva
,
T. A.
, and
Einziger
,
R. E.
,
2013
, “
Ductile-to-Brittle Transition Temperature for High-Burnup Cladding Alloys Exposed to Simulated Drying-Storage Conditions
,”
J. Nucl. Mater.
,
433
(
1–3
), pp.
431
448
.
9.
Birch
,
R.
,
Wang
,
S.
,
Tong
,
V.
, and
Britton
,
B.
,
2018
, “
Short Communication: The Effect of Cooling Rate and Grain Size on Hydride Microstructure in Zircaloy-4
,”
arXiv
,
513
, pp.
221
225
.
10.
Leitch
,
B. W.
, and
Puls
,
M. P.
,
1992
, “
Finite Element Calculations of the Accommodation Energy of a Misfitting Precipitate in an Elastic-Plastic Matrix
,”
Metall. Trans. A
,
23
(
3
), pp.
797
806
.
11.
Motta
,
A. T.
,
Capolungo
,
L.
,
Chen
,
L. Q.
,
Cinbiz
,
M. N.
,
Daymond
,
M. R.
,
Koss
,
D. A.
,
Lacroix
,
E.
,
Pastore
,
G.
,
Simon
,
P. C. A.
,
Tonks
,
M. R.
,
Wirth
,
B. D.
, and
Zikry
,
M. A.
,
2019
, “
Hydrogen in Zirconium Alloys: A Review
,”
J. Nucl. Mater.
,
518
, pp.
440
460
.
12.
Ashmawi
,
W. M.
, and
Zikry
,
M. A.
,
2002
, “
Prediction of Grain-Boundary Interfacial Mechanisms in Polycrystalline Materials
,”
ASME J. Eng. Mater. Technol.
,
124
(
1
), pp.
88
96
.
13.
Shanthraj
,
P.
, and
Zikry
,
M. A.
,
2011
, “
Dislocation Density Evolution and Interactions in Crystalline Materials
,”
Acta Mater.
,
59
(
20
), pp.
7695
7702
.
14.
Kubin
,
L.
,
Devincre
,
B.
, and
Hoc
,
T.
,
2008
, “
Modeling Dislocation Storage Rates and Mean Free Paths in Face-Centered Cubic Crystals
,”
Acta Mater.
,
56
(
20
), pp.
6040
6049
.
15.
Kim
,
J.-S.-S.
,
Kim
,
T.-H.-H.
,
Kook
,
D.-H.-H.
, and
Kim
,
Y.-S.-S.
,
2015
, “
Effects of Hydride Morphology on the Embrittlement of Zircaloy-4 Cladding
,”
J. Nucl. Mater.
,
456
, pp.
235
245
.
16.
Chan
,
H.
,
Roberts
,
S. G.
, and
Gong
,
J.
,
2016
, “
Micro-scale Fracture Experiments on Zirconium Hydrides and Phase Boundaries
,”
J. Nucl. Mater.
,
475
, pp.
105
112
.
17.
Gu
,
X. F.
, and
Zhang
,
W. Z.
,
2014
, “
A Simple Method for Calculating the Possible Habit Planes Containing One Set of Dislocations and Its Applications to Fcc/Bct and Hcp/Bcc Systems
,”
Metall. Mater. Trans. A
,
45
(
4
), pp.
1855
1865
.
18.
Motta
,
A. T.
,
Yilmazbayhan
,
A.
,
da Silva
,
M. J. G.
,
Comstock
,
R. J.
,
Was
,
G. S.
,
Busby
,
J. T.
,
Gartner
,
E.
,
Peng
,
Q.
,
Jeong
,
Y. H.
, and
Park
,
J. Y.
,
2007
, “
Zirconium Alloys for Supercritical Water Reactor Applications: Challenges and Possibilities
,”
J. Nucl. Mater.
,
371
(
1–3
), pp.
61
75
.
19.
Porter
,
D. A.
, and
Easterling
,
K. E.
,
2009
,
Phase Transformations in Metals and Alloys
, (revised reprint), 3rd ed.,
CRC Press
,
Boca Raton, FL
.
20.
Mohamed
,
I.
, and
Zikry
,
M. A.
,
2021
, “
Modeling of the Microstructural Behavior of Hydrided Zirconium Alloys
,”
Comput. Mech.
21.
Ziaei
,
S.
,
Wu
,
Q.
, and
Zikry
,
M. A.
,
2015
, “
Orientation Relationships Between Coherent Interfaces in Hcp-Fcc Systems Subjected to High Strain-Rate Deformation and Fracture Modes
,”
J. Mater. Res.
,
30
(
15
), pp.
2348
2359
.
22.
Kubo
,
T.
,
Wakashima
,
Y.
,
Amano
,
K.
, and
Nagai
,
M.
,
1985
, “
Effects of Crystallographic Orientation on Plastic Deformation and SCC Initiation of Zirconium Alloys
,”
J. Nucl. Mater.
,
132
(
1
), pp.
1
9
.
23.
Musienko
,
A.
, and
Cailletaud
,
G.
,
2009
, “
Simulation of Inter- and Transgranular Crack Propagation in Polycrystalline Aggregates Due to Stress Corrosion Cracking
,”
Acta Mater.
,
57
(
13
), pp.
3840
3855
.
24.
Pouillier
,
E.
,
Gourgues
,
A. F.
,
Tanguy
,
D.
, and
Busso
,
E. P.
,
2012
, “
A Study of Intergranular Fracture in an Aluminium Alloy Due to Hydrogen Embrittlement
,”
Int. J. Plast.
,
34
, pp.
139
153
.
25.
Hansbo
,
A.
, and
Hansbo
,
P.
,
2004
, “
A Finite Element Method for the Simulation of Strong and Weak Discontinuities in Solid Mechanics
,”
Comput. Meth. Appl. Mech. Eng.
,
193
(
33–35
), pp.
3523
3540
.
26.
Wu
,
Q.
, and
Zikry
,
M. A.
,
2014
, “
Microstructural Modeling of Crack Nucleation and Propagation in High Strength Martensitic Steels
,”
Int. J. Solids Struct.
,
51
(
25–26
), pp.
4345
4356
.
27.
Morris
,
J. W.
,
2011
, “
On the Ductile-Brittle Transition in Lath Martensitic Steel
,”
ISIJ Int.
,
51
(
10
), pp.
1569
1575
.
28.
Zikry
,
M. A.
,
1994
, “
An Accurate and Stable Algorithm for High Strain-Rate Finite Strain Plasticity
,”
Comput. Struct.
,
50
(
3
), pp.
337
350
.
29.
Kubo
,
T.
,
Kobayashi
,
Y.
, and
Uchikoshi
,
H.
,
2013
, “
Determination of Fracture Strength of δ-Zirconium Hydrides Embedded in Zirconium Matrix at High Temperatures
,”
J. Nucl. Mater.
,
435
(
1–3
), pp.
222
230
.
30.
Oono
,
N.
,
Kasada
,
R.
,
Higuchi
,
T.
,
Sakamoto
,
K.
,
Nakatsuka
,
M.
,
Hasegawa
,
A.
,
Kondo
,
S.
,
Iwata
,
N. Y.
,
Matsui
,
H.
, and
Kimura
,
A.
,
2013
, “
Comparison of Irradiation Hardening and Microstructure Evolution in Ion-Irradiated Delta and Epsilon Hydrides
,”
J. Nucl. Mater.
,
442
(
1–3
SUPPL.1), pp.
S826
S829
.
31.
Hong
,
S. I.
, and
Lee
,
K. W.
,
2005
, “
Stress-Induced Reorientation of Hydrides and Mechanical Properties of Zircaloy-4 Cladding Tubes
,”
J. Nucl. Mater.
,
340
(
2–3
), pp.
203
208
.
32.
Daum
,
R. S.
,
Majumdar
,
S.
,
Liu
,
Y.
, and
And Billone
,
M. C.
,
2006
, “
Radial-Hydride Embrittlement of High-Burnup Zircaloy-4 Fuel Cladding
,”
J. Nucl. Sci. Technol.
,
43
(
9
), pp.
1054
1067
.