Analytical and Scale-Free Phase-Field Studies of $\alpha$ to $\omega$ Phase Transformation in Single Crystal Zirconium under Nonhydrostatic Loadings

TL;DR

This paper develops a scale-free phase-field model to study the alpha to omega phase transformation in single crystal zirconium under nonhydrostatic stresses, providing analytical solutions and FEM simulations that align with experimental observations.

Contribution

It introduces a novel scale-free phase-field approach for multivariant phase transformations in zirconium under complex loadings, clarifying stress effects and deriving explicit transformation conditions.

Findings

Stress effects on transformation pressures are quantified.

Analytical solutions for stress-strain and phase fractions are provided.

FEM simulations validate analytical predictions and match experimental trends.

Abstract

Zirconium (Zr) is an important engineering material with numerous practical applications. It undergoes martensitic $\alpha$ to $\omega$ phase transformation (PT) at pressures that vary from 0.67 GPa to 17 GPa under different loading conditions. Despite numerous experimental and theoretical studies, the effect of the nonhydrostatic stresses is not well understood. To separate the effect of nonhydrostatic stresses from the plastic deformation, a scale-free phase field approach (PFA) for multivariant $\alpha$ to $\omega$ PT in a single crystal Zr under general nonhydrostatic loadings is presented. Explicit conditions for the direct and reverse PTs between austenite and martensitic variants and between martensitic variants under general stress tensor are derived and analyzed. In particular, the effect of the deviatoric stresses on the PT pressures is elucidated. It is shown that their effect cannot explain much larger reduction in the transformation pressure observed during plastic flow, i.e., specific mechanisms of strain-induced phase transformations should be involved. Under assumption of the homogeneous fields in the sample, complete analytical solutions that include stress-strain curves during the PT, PT start and finish stresses (i.e., stress hysteresis), and volume fraction of the variants, are determined for different loadings. Finite element method (FEM) solutions are found for the phase field simulations of the microstructure evolution for the same loadings, as well as for two grains of the polycrystalline sample. Macroscopic averaged characteristics of the PFA solutions are well described by an analytical solution, which also simplifies their interpretations. Obtained results are in good qualitative agreement with existing experiments. In addition, some controversies of the previous approaches are analysed.