Interfacial energy and dissipation in martensitic phase transformations. Part I: Theory

A model of evolving martensitic microstructures is formulated that incorporates the interfacial energy and dissipation on three different scales corresponding to the grain boundaries attained by martensite plates, the interfaces between austenite and martensite plates, and the twin interfaces within martensite plates. Three different time scales are also considered in order to clarify the meaning of rate-independent dissipation related to instabilities at more refined temporal and spatial scales. On the slowest time scale, the process of deformation and martensitic phase transformation is investigated as being composed of segments of smooth quasi-static evolution separated by sudden jumps associated with creation or annihilation of interfaces. A general evolution rule is used in the form of minimization of the incremental energy supply to the whole compound thermodynamic system, including the rate-independent dissipation. Close relationship is shown between the evolution rule and the thermodynamic condition for stability of equilibrium, with no a priori assumption on convexity of the dissipation function. It is demonstrated that the extension of the minimum principle from the first-order rates to small but finite increments requires a separate symmetry restriction imposed on the state derivative of the dissipation function. Formulae for the dissipation associated with annihilation of interfaces are proposed that exhibit limited path-independence and satisfy that symmetry requirement. By exploiting the incremental energy minimization rule with the help of the transport theorems, local propagation conditions are derived for both planar and curved phase transformation fronts. The theory serves as a basis for the algorithm for calculation of the stress-induced evolution of martensitic microstructures along with their characteristic dimensions and related hysteresis loops in shape memory alloys; the examples are given in Part II of the paper.     

Size effects on the martensitic phase transformation of NiTi nanograins

The analysis of nanocrystalline NiTi by transmission electron microscopy (TEM) shows that the martensitic transformation proceeds by the formation of atomic-scale twins. Grains of a size less than about 50 nm do not transform to martensite even upon large undercooling. A systematic investigation of these phenomena was carried out elucidating the influence of the grain size on the energy barrier of the transformation. Based on the experiment, nanograins were modeled as spherical inclusions containing (0 0 1) compound twinned martensite. Decomposition of the transformation strains of the inclusions into a shear eigenstrain and a normal eigenstrain facilitates the analytical calculation of shear and normal strain energies in dependence of grain size, twin layer width and elastic properties. Stresses were computed analytically for special cases, otherwise numerically. The shear stresses that alternate from twin layer to twin layer are concentrated at the grain boundaries causing a contribution to the strain energy scaling with the surface area of the inclusion, whereas the strain energy induced by the normal components of the transformation strain and the temperature dependent chemical free energy scale with the volume of the inclusion. In the nanograins these different energy contributions were calculated which allow to predict a critical grain size below which the martensitic transformation becomes unlikely. Finally, the experimental result of the atomic-scale twinning can be explained by analytical calculations that account for the transformation-opposing contributions of the shear strain and the twin boundary energy of the twin-banded morphology of martensitic nanograins.     

Size effects in martensitic microstructures: Finite-strain phase field model versus sharp-interface approach

A finite-strain phase field model for martensitic phase transformation and twinning in shape memory alloys is developed and confronted with the corresponding sharp-interface approach extended to interfacial energy effects. The model is set in the energy framework so that the kinetic equations and conditions of mechanical equilibrium are fully defined by specifying the free energy and dissipation potentials. The free energy density involves the bulk and interfacial energy contributions, the latter describing the energy of diffuse interfaces in a manner typical for phase-field approaches. To ensure volume preservation during martensite reorientation at finite deformation within a diffuse interface, it is proposed to apply linear mixing of the logarithmic transformation strains. The physically different nature of phase interfaces and twin boundaries in the martensitic phase is reflected by introducing two order-parameters in a hierarchical manner, one as the reference volume fraction of austenite, and thus of the whole martensite, and the second as the volume fraction of one variant of martensite in the martensitic phase only. The microstructure evolution problem is given a variational formulation in terms of incremental fields of displacement and order parameters, with unilateral constraints on volume fractions explicitly enforced by applying the augmented Lagrangian method. As an application, size-dependent microstructures with diffuse interfaces are calculated for the cubic-to-orthorhombic transformation in a CuAlNi shape memory alloy and compared with the sharp-interface microstructures with interfacial energy effects.     

Modelling of laminated microstructures in stress-induced martensitic transformations

This paper is concerned with micromechanical modelling of stress-induced martensitic transformations in crystalline solids, with the focus on distinct elastic anisotropy of the phases and the associated redistribution of internal stresses. Micro-macro transition in stresses and strains is analysed for a laminated microstructure of austenite and martensite phases. Propagation of a phase transformation front is governed by a time-independent thermodynamic criterion. Plasticity-like macroscopic constitutive rate equations are derived in which the transformed volume fraction is incrementally related to the overall strain or stress. As an application, numerical simulations are performed for cubic β₁ (austenite) to orthorhombic γ₁′ (martensite) phase transformation in a single crystal of Cu-Al-Ni shape memory alloy. The pseudoelasticity effect in tension and compression is investigated along with the corresponding evolution of internal stresses and microstructure.     

An upper bound to the free energy of n-variant polycrystalline shape memory alloys

In the last two decades, the problem of computing the elastic energy of phase transforming materials has been studied by a variety of research groups. Due to the non-quasiconvexity of the underlying multi-well landscape, different relaxation methods have been used in order to estimate the quasiconvex envelope of the energy density, for which no explicit expression is known at present.This paper combines a recently developed lamination bound for monocrystalline shape memory alloys which relies on martensitic twinned microstructures with the work of Smyshlyaev and Willis [1998a. A ‘non-local’ variational approach to the elastic energy minimization of martensitic polycrystals. Proc. R. Soc. London A 454, 1573–1613]. As a result, a lamination upper bound for n-variant polycrystalline martensitic materials is obtained.The lamination bound is then compared with Reuß- and Taylor-type estimates. While, for given volume fractions, good agreement of lamination upper and convexification lower bounds is obtained, a comparison using energy-minimizing volume fractions computed from the various bounds yields larger differences. Finally, we also investigate the influence of the polycrystal's texture. For a strong ellipsoidal texture, we observe even better agreement of upper and lower bounds than for the case of isotropic statistics.     

Coupling between martensitic phase transformations and plasticity: A microstructure-based finite element model

A microstructural finite element (MFE) model is developed to capture the interaction between martensitic transformations and plasticity in NiTi shape memory alloys (SMAs). The interaction is modeled through the grain-to-grain redistribution of stress caused by both plasticity and phase transformation, so that each mechanism affects the driving force of the other. A unique feature is that both processes are modeled at a crystallographic level and are allowed to operate simultaneously. The model is calibrated to pseudoelastic data for select single crystals of Ti–50.9at.%Ni. For polycrystals, plasticity is predicted to enhance the overall martensite volume fraction at a given applied stress. Upon unloading, residual stress can induce remnant (retained) martensite. For thermal cycling under load bias, plasticity is observed to limit the net transformation strain/cycle and increase the hysteretic width. Deformation processing, via plastic pre-straining at elevated temperature, is shown to dramatically alter subsequent pseudoelastic response, as well as induce two-way shape memory behavior during no-load thermal cycling. Overall, the model is suitable at smaller imposed strains, where martensite detwinning is not expected to dominate.     

The influence of surface energy on stress-free microstructures in shape memory alloys

Certain alloys such as In-Tl, Ni-Ti, Ag-Cd or Cu-Al-Ni display remarkable mechanical properties such as the shape memory effect or pseudo-elasticity. This behaviour is related to a solid-solid phase transformation which leads to a complicated microscopic arrangement of different phases. In recent studies such microstructures have been analyzed by the minimization of elastic energy. We discuss the influence of additional surface energy terms on the existence of stress-free microstructures both in the nonlinear and a geometrically linear setting.

---

Sommario Certe leghe come quelle di In-Tl, Ni-Ti, Ag-Cd o Cu-Al-Ni mostrano proprietà meccaniche notevoli quali la memoria di forma o la pseudoelasticità. Questo comportamento è determinato da una trasformazione di fase solido-solido che conduce a complicati arrangiamenti a livelo microscopico. In studi recenti tali microstrutture sono state analizzate attraverso la minimizzazione dell'energia elastica. Noi discutiamo l'influenza di termini addizionali di energia superficiale sull'esistenza di microstrutture in uno stato naturale sia in un contesto lineare che non lineare.

     

Application of a partially relaxed shape memory free energy function to estimate the phase diagram and predict global microstructure evolution

In this work recent results regarding bounds of the relaxed free energy functions of a broad class of shape memory materials of arbitrary symmetry are leveraged to develop a simple and efficient numerical method to analyze several aspects of the thermomechanical response of such materials. This approach is shown to be useful for construction of the austenitic phase diagram used in many early phenomenological models. It is also demonstrated that the resulting implementation is suitable for finite element analysis, has desirable numerical properties, and makes realistic quantitative predictions regarding the evolution of the lattice correspondence variants, transformation stress and strain levels, and orientation dependence.     

Micromechanical modeling of stress-induced phase transformations. Part 2. Computational algorithms and examples

Based on the theory developed in Part 1 of this paper [Levitas, V.I., Ozsoy, I.B., 2008. Micromechanical modeling of stress-induced phase transformations. Part 1. Thermodynamics and kinetics of coupled interface propagation and reorientation. Int. J. Plasticity. doi:10.1016/j.ijplas.2008.02.004], various non-trivial examples of microstructure evolution under complex multiaxial loading are presented. For the case without interface rotation, the effect of the athermal thresholds for austenite (A)–martensite (M) and martensitic variant M_I–variant M_II interfaces and loading paths on stress–strain curves and phase transformations was studied. For coupled interface propagation and rotation, two types of numerical simulations were carried out. The tetragonal–orthorhombic transformation has been studied under general three-dimensional interface orientation and zero athermal threshold. The cubic–tetragonal transformation was treated with allowing for an athermal threshold and interface reorientation within a plane. The effect of the athermal threshold, the number of martensitic variants and an interface orientation in the embryo was studied in detail. It was found that an instability in the interface normal leads to a jump-like interface reorientation that has the following features of the energetics of a first-order transformation: there are multiple energy minima versus interface orientation that are separated by an energy barrier; positions of minima do not change during loading but their depth varies; when the barrier disappears (i.e. one of the minima transforms to the local saddle or maximum points), the system rapidly evolves toward another stable orientation. Depending on the loading and material parameters, we observed a large continuous change in interface orientation, a jump in interface reorientation, a jump in volume fractions and stresses, an expected stress relaxation during the phase transition and unexpected stress growth during the transition because of large change in elastic moduli.     

Interaction between phase transformations and dislocations at the nanoscale. Part 1. General phase field approach

Thermodynamically consistent, three-dimensional (3D) phase field approach (PFA) for coupled multivariant martensitic transformations (PTs), including cyclic PTs, variant–variant transformations (i.e., twinning), and dislocation evolution is developed at large strains. One of our key points is in the justification of the multiplicative decomposition of the deformation gradient into elastic, transformational, and plastic parts. The plastic part includes four mechanisms: dislocation motion in martensite along slip systems of martensite and slip systems of austenite inherited during PT and dislocation motion in austenite along slip systems of austenite and slip systems of martensite inherited during reverse PT. The plastic part of the velocity gradient for all these mechanisms is defined in the crystal lattice of the austenite utilizing just slip systems of austenite and inherited slip systems of martensite, and just two corresponding types of order parameters. The explicit expressions for the Helmholtz free energy and the transformation and plastic deformation gradients are presented to satisfy the formulated conditions related to homogeneous thermodynamic equilibrium states of crystal lattice and their instabilities. In particular, they result in a constant (i.e., stress- and temperature-independent) transformation deformation gradient and Burgers vectors. Thermodynamic treatment resulted in the determination of the driving forces for change of the order parameters for PTs and dislocations. It also determined the boundary conditions for the order parameters that include a variation of the surface energy during PT and exit of dislocations. Ginzburg–Landau equations for dislocations include variation of properties during PTs, which in turn produces additional contributions from dislocations to the Ginzburg–Landau equations for PTs. A complete system of coupled PFA and mechanics equations is presented. A similar theory can be developed for PFA to dislocations and other PTs, like reconstructive PTs and diffusive PTs described by the Cahn–Hilliard equation, as well as twinning and grain boundaries evolution.     

Interface propagation and microstructure evolution in phase field models of stress-induced martensitic phase transformations

Analytical solutions for diffuse interface propagation are found for two recently developed Landau potentials that account for the phenomenology of stress-induced martensitic phase transformations. The solutions include the interface profile and velocity as a function of temperature and stress tensor. An instability in the interface propagation near lattice instability conditions is studied numerically. The effect of material inertia is approximately included. Two methods for introducing an athermal interface friction in phase field models are discussed. In the first method an analytic expression defines the location of the diffuse interface, and the rate of change of the order parameters is required to vanish if the driving force is below a threshold. As an alternative and more physical approach, we demonstrate that the introduction of spatially oscillatory stress fields due to crystal defects and the Peierls barrier, or to a jump in chemical energy, reproduces the effect of an athermal threshold. Finite element simulations of microstructure evolution with and without an athermal threshold are performed. In the presence of spatially oscillatory fields the evolution self-arrests in realistic stationary microstructures, thus the system does not converge to an unphysical single-phase final state, and rate-independent temperature- and stress-induced phase transformation hysteresis are exhibited.     

Phase transitions induced by extension in a slender SMA cylinder: Analytical solutions for the hysteresis loop based on a quasi-3D continuum model

In this paper, we propose a quasi-3D continuum model to study the rate-independent hysteresis phenomenon in phase transitions of a slender shape memory alloy (SMA) cylinder subject to the uniaxial tension. Based on the three-dimensional field equations and the traction-free boundary conditions, by using a coupled series-asymptotic expansion method, we manage to express the total elastic potential energy of the cylinder in terms of the leading order term of the axial strain. We further consider the rate-independent dissipation effect in a purely one-dimensional setting. The mechanical dissipation functions are also expressed in terms of the axial strain. The equilibrium configuration of the cylinder is then determined by using the principle of maximizing the total energy dissipation. An illustrative example with some special chosen material constants is further considered. Free end boundary conditions are proposed at the two ends of the cylinder. By conducting a phase plane analysis and through some calculations, we obtain the analytical solutions of the equilibrium equation. We find that the engineering stress–strain curves corresponding to the obtained solutions can capture some important features of the experimental results. It appears that the analytical results obtained in this paper reveal the multiple solutions nature of the problem and shed certain light on the instability phenomena during the phase transition process.     

Micromechanical modeling of stress-induced phase transformations. Part 1. Thermodynamics and kinetics of coupled interface propagation and reorientation

The universal (i.e. independent of the constitutive equations) thermodynamic driving force for coherent interface reorientation during first-order phase transformations in solids is derived for small and finite strains. The derivation is performed for a representative volume with plane interfaces, homogeneous stresses and strains in phases and macroscopically homogeneous boundary conditions. Dissipation function for coupled interface (or multiple parallel interfaces) reorientation and propagation is derived for combined athermal and drag interface friction. The relation between the rates of single and multiple interface reorientation and propagation and the corresponding driving forces are derived using extremum principles of irreversible thermodynamics. They are used to derive complete system of equations for evolution of martensitic microstructure (consisting of austenite and a fine mixture of two martensitic variants) in a representative volume under complex thermomechanical loading. Viscous dissipation at the interface level introduces size dependence in the kinetic equation for the rate of volume fraction. General relationships for a representative volume with moving interfaces under piece-wise homogeneous boundary conditions are derived. It was found that the driving force for interface reorientation appears when macroscopically homogeneous stress or strain are prescribed, which corresponds to experiments. Boundary conditions are satisfied in an averaged way. In Part 2 of the paper [Levitas, V.I., Ozsoy, I.B., 2008. Micromechanical modeling of stress-induced phase transformations. Part 2. Computational algorithms and examples. Int. J. Plasticity (2008)], the developed theory is applied to the numerical modeling of the evolution of martensitic microstructure under three-dimensional thermomechanical loading during cubic-tetragonal and tetragonal-orthorhombic phase transformations.     

Internal loops in pseudoelastic behaviour of Ti-Ni shape memory alloys: Experiment and modelling

The stress-strain isothermal hysteresis loops due to the incomplete martensitic transformation are analysed for Ti-Ni shape memory alloys. Experiments show the existence of two distinct yield lines for phase transition; one for the forward transformation austenitemartensite (AM), the other for the reverse transformation MA. The tensile behaviour of single crystals with only one yield line (AM) [1] can be considered as an ideal case. An extension of a thermodynamic model for pseudoelasticity [2] allows these two yield lines to be taken into account.

---

Sommario Per leghe Ti-Ni con memoria di forma vengono analizzati i cicli di isteresi isotermici tensione-deformazione prodotti da una incompleta trasformazione martensitica. Gli esperimenti mostrano l'esistenza di due distinte linee di snervamento per la transizione di fase, una verso la trasformazione austenitemartensite (AM), l'altra per la trasformazione inversa MA. Il comportamento a trazione di un singolo cristallo con una sola linea di snervamento (AM) [1], può essere considerato un caso ideale. L'estensione ad un modello termodinamico pseudo-elastico [2] consente di analizzare queste due linee di snervamento.

     

Interfacial excess energy, excess stress and excess strain in elastic solids: Planar interfaces

In this paper, interfacial excess energy and interfacial excess stress for coherent interfaces in an elastic solid are reformulated within the framework of continuum mechanics. It is shown that the well-known Shuttleworth relationship between the interfacial excess energy and interfacial excess stress is valid only when the interface is free of transverse stresses. To account for the transverse stress, a new relationship is derived between the interfacial excess energy and interfacial excess stress. Dually, the concept of transverse interfacial excess strain is also introduced, and the complementary Shuttleworth equation is derived that relates the interfacial excess energy to the newly introduced transverse interfacial excess strain. This new formulation of interfacial excess stress and excess strain naturally leads to the definition of an in-plane interfacial stiffness tensor, a transverse interfacial compliance tensor and a coupling tensor that accounts for the Poisson's effect of the interface. These tensors fully describe the elastic behavior of a coherent interface upon deformation.     

Transition waves in bistable structures. II. Analytical solution: wave speed and energy dissipation

We consider dynamics of chains of rigid masses connected by links described by irreversible, piecewise linear constitutive relation: the force-elongation diagram consists of two stable branches with a jump discontinuity at the transition point. The transition from one stable state to the other propagates along the chain and excites a complex system of waves. In the first part of the paper (Cherkaev et al., 2004, Transition waves in bistable structures. I. Delocalization of damage), the branches could be separated by a gap where the tensile force is zero, the transition wave was treated as a wave of partial damage. Here we assume that there is no zero-force gap between the branches. This allows us to obtain steady-state analytical solutions for a general piecewise linear trimeric diagram with parallel and nonparallel branches and an arbitrary jump at the transition. We derive necessary conditions for the existence of the transition waves and compute the speed of the wave. We also determine the energy of dissipation which can be significantly increased in a structure characterized by a nonlinear discontinuous constitutive relation. The considered chain model reveals some phenomena typical for waves of failure or crushing in constructions and materials under collision, waves in a structure specially designed as a dynamic energy absorber and waves of phase transitions in artificial and natural passive and active systems.     

The free energy of mixing for n-variant martensitic phase transformations using quasi-convex analysis

The construction of effective models for materials that undergo martensitic phase transformations requires usable and accurate functional representations for the free energy density. The general representation of this energy is known to be highly non-convex; it even lacks the property of quasi-convexity. A quasi-convex relaxation, however, does permit one to make certain estimates and powerful conclusions regarding phase transformation. The general expression for the relaxed free energy is however not known in the n-variant case. Analytic solutions are known only for up to 3 variants, whereas cases of practical interests involve 7-13 variants. In this study we examine the n-variant case utilizing relaxation theory and produce a seemingly obvious but very powerful observation regarding a lower bound to the quasi-convex relaxation that makes practical evolutionary computations possible. We also examine in detail the 4-variant case where we explicitly show the relation between three different forms of the free energy of mixing: upper bound by lamination, the Reuß lower bound, and a lower estimate of the -measure bound. A discussion of the bounds and their utility is provided; sample computations are presented for illustrative purposes.     

Magnetic field-induced martensitic phase transformation in magnetic shape memory alloys: Modeling and experiments

In this work, a continuum based model of the magnetic field induced phase transformation (FIPT) for magnetic shape memory alloys (MSMA) is developed. Hysteretic material behaviors are considered through the introduction of internal state variables. A Gibbs free energy is proposed using group invariant theory and the coupled constitutive equations are derived in a thermodynamically consistent way. An experimental procedure of FIPT in NiMnCoIn MSMA single crystals, which can operate under high blocking stress, is described. The model is then reduced to a 1-D form and the material parameter identification from the experimental results is discussed. Model predictions of magneto-thermo-mechanical loading conditions are presented and compared to experiments.