A 2D finite element based on a nonlocal constitutive model describing localization and propagation of phase transformation in shape memory alloy thin structures

Here, the effects of localization and propagation of martensitic phase transformation on the response of SMA thin structures subjected to thermo-mechanical loadings are investigated using nonlocal constitutive model in conjunction with finite element method. The governing equations are derived based on variational principle considering thermo-mechanical equilibrium and the spatial distribution of the nonlocal volume fraction of martensite during transformation. The nonlocal volume fraction of martensite is defined as a weighted average of the local volume fraction of martensite over a domain characterized by an internal length parameter. The local version of the thermo-mechanical behavior model derived from micromechanics considers the local volume fraction of martensite and the mean transformation strain. A 4-noded quadrilateral plane stress element with three degrees of freedom per node accounting for in-plane displacements and the nonlocal volume fraction of martensite is developed. Numerical simulations are conducted to bring out the influence of material and geometrical heterogeneities (perturbations/defects) on the localization and propagation of phase transformation in SMA thin structures. Also, a sensitivity analysis of the material response due to the localization and the other related model parameters is carried out. The detailed investigation done here clearly shows that the localization of phase transformation has significant effect on the response of shape memory alloys.     

Martensite Strain Memory in the Shape Memory Alloy Nickel-Titanium Under Mechanical Cycling

This paper describes an experimental study of stress-induced martensitic phase transformation in the SMA Nickel-Titanium. The rich local thermo-mechanical interactions that underlie transformation are examined using three-dimensional Digital Image Correlation (strain fields) and infrared imaging (thermal fields). We quantify the complex local interactions between released/absorbed latent heat and the extent of transformation, and explore the characteristics of the phase fronts and the evolution of martensitic volume fraction. We also quantify a strong strain memory in the martensite that forms in the wake of the phase transformation front. The accommodated strain in the martensite will remain constant during loading, even as the existing phase front propagates. There also exists a remarkable strain memory in the martensite that persists from cycle to cycle, indicating that the local elastic stress fields in the martensite are driven by a dislocation structure and martensitic nuclei that largely stabilize during the first loading cycle.     

An inclusion theory for the propagation of martensite band in NiTi shape memory alloy wires under tension

In this paper, we present analytical modelling for the pure mechanical response of uniaxial tensioned NiTi wire that experiences stress-induced martensitic transformation via a propagating martensite band at the superelastic temperature regime. The model aims to predict the overall behavior of the SMA wire as a structural response containing propagating instabilities. Based on the systematic experimental investigation of Shaw and Kyriakides (Shaw, J.A., Kyriakides, S., 1995. Thermomechemical aspects of NiTi. J. Mech. Phys. Solids 43, T243–1281 and Shaw, J.A., Kyriakides, S., 1997. On the nucleation and propogation of phase transformation fronts in a NiTi alloy. Acta Mater. 45(2), 638–700), the wire is modeled as an elastic rod containing a single cylindrical transformation inclusion with a uniform axisymmetric eigenstrain. The analytical expression of the free energy of this special matrix-inclusion system is formulated and the length of the martensite band is identified as the key variable describing the transformation process of the system. Theoretical predictions on the peak stress and the subsequent steady-state propagation stress of the wire during forward and reverse transformations are provided and compared with the available experimental data. Specimen size effect on the nominal stress-strain curves and general deformation features of the wire are discussed.     

Anisotropy of martensitic transformations in modeling of shape memory alloy polycrystals

Based on the knowledge of the anisotropy associated with the martensitic transformations obtained from tension/compression experiments with oriented CuAlNi single crystals, a simple constant stress averaging approach is employed to model the SMA polycrystal deformation behaviors. Only elastic and inelastic strains due to the martensitic transformation, variant reorientations in the martensite phase and martensite to martensite transformations in thermomechanical loads are considered. The model starts from theoretical calculation of the stress-temperature transformation conditions and their orientation dependence from basic crystallographic and material attributes of the martensitic transformations. Results of the simulations of the NiTi, NiAl, and Cu-based SMA polycrystals in stress–strain tests are shown. It follows that SMA polycrystals, even with randomly oriented grains, typically exhibit tension/compression asymmetry of the shape of the pseudoelastic σ−ε curves in transformation strain, transformation stress, hysteresis widths, character of the pseudoelastic flow and in the slope of temperature dependence of the transformation stresses. It is concluded that some macroscopic features of the SMA polycrystal behaviors originate directly from the crystallography of the undergoing MT's. The model shows clearly the crystallographic origin of these phenomena by providing a link from the crystallographic and material attributes of martensitic transformations towards the macroscopic σ−ε−T behaviors of SMA polycrystals.     

A two-level micromechanical theory for a shape-memory alloy reinforced composite

A two-level micromechanical theory is developed to study the influence of the shape and volume concentration of shape-memory alloy (SMA) inclusions on the overall stress–strain behavior of a SMA-reinforced composite. The first level exists on the smaller SMA level, in which, under the action of stress, parent austenite may transform into martensite. The second level is on the larger scale consisting of the metastable SMA inclusions and an inactive polymer matrix. The evolution of martensite microstructure is evaluated from the irreversible thermodynamics, in conjunction with the micromechanics and physics of martensitic transformation. By taking martensite to exist in the form of thin plates on the micro scale and assuming SMA inclusions to be homogeneously aligned spheroids on the macro scale, the overall stress–strain behaviors of a NiTi-reinforced composite are calculated for various SMA shapes and concentrations. The results indicate that, under a tensile axial loading, martensitic transformation is easier to take place when SMA inclusions exist in the form of long fibers, but most difficult to occur when they are in the form of flat discs. In general the levels of the applied stress at which martensite transformation commences, finishes, and austenitic transformation starts, and finishes, are found to decrease with increasing aspect ratio of the SMA inclusions while the damping capacity increases with it; these properties point to the advantage of using fibrous composites for actuators or sensors under a tensile loading.     

A non-local implicit gradient-enhanced model for unstable behaviors of pseudoelastic shape memory alloys in tensile loading

In this paper, a gradient-enhanced 3-D phenomenological model for shape memory alloys using the non-local theory is developed based on a 1-D constitutive model. The method utilizes a non-local field variable in its constitutive framework with an implicit gradient formulation in order to achieve results independent of the finite element discretization. An efficient numerical approach to implement the non-local gradient-enhanced model in finite element codes is proposed. The model is used to simulate stress drop at the onset of transformation, and its performance is evaluated using different experimental data. The potential of the presented numerical approach for behavior of shape memory alloys in eliminating mesh-dependent simulations is validated by conducting various localization problems. The numerical results show that the developed model can simulate the observed unstable behaviors such as stress drop and deviation of local strain from global strain during nucleation and propagation of martensitic phase.     

超弹性形状记忆合金管单向拉伸试验的数值模拟

NiTi形状记忆合金具有很强的超弹性行为,这种超弹性行为是由于材料在应力作用下发生可逆的马氏体相变所引起。最近Sun和Lee^[4]在NiTi形状记忆合金管单向拉伸试验中观测到,应力诱导马氏体相变具有螺旋带状的形貌特征,本文对此作了数值模拟研究。采用包含应变软化效应的三线性本构关系,建立了NiTi形状记忆合金管的三维有限元模型。通过迭代计算,成功地再现了试验中所观察到的螺旋状相变带从形成到长大的全过程。数值计算结果表明,产生这一独特现象的力学机制,在于NiTi形状记忆合金管在拉伸状态下出现的局部变形失稳极其传播。     

A thermodynamical constitutive model for shape memory materials. Part I. The monolithic shape memory alloy

Pseudoelasticity and the shape memory effect (SME) due to martensitic transformation and reorientation of polycrystalline shape memory alloy (SMA) materials are modeled using a free energy function and a dissipation potential. Three different cases are considered, based on the number of internal state variables in the free energy: (1) austenite plus a variable number of martensite variants; (2) austenite plus two types of martensite; and (3) austenite and one type of martensite. Each model accounts for three-dimensional simultaneous transformation and reorientation. The single-martensite model was chosen for detailed study because of its simplicity and its ease of experimental verification. Closed form equations are derived for the damping capacity and the actuator efficiency of converting heat into work. The first law of thermodynamics is used to demonstrate that significantly more work is required to complete the adiabatic transformation than the isothermal transformation. Also, as the hardening due to the austenite/martensite misfit stresses approaches zero, the transformation approaches the isothermal, infinite specific heat conditions of a first-order transformation. In a second paper, the single-martensite model is used in a mesomechanical derivation of the constitutive equations of an active composite with an SMA phase.     

Effect of microstructure on the hardening and softening behavIors of polycrystalline shape memory alloys Part I: Micromechanics constitutive modeling

The effects of microstructure and its evolution on the macroscopic superelastic stress-strain response of polycrystalline Shape Memory Alloy (SMA) are studied by a microstructure-based constitutive model developed in this paper. The model is established on the following basis: (1) the transformation conditions of the unconstrained single crystal SMA microdomain (to be distinguished from the bulk single crystal), which serve as the local criterion for the derivation of overall transformation yield conditions of the polycrystal; (2) the micro- to macro-transition scheme by which the connection between the polycrystal aggregates and the single crystal microdomain is established and the macroscopic transformation conditions of the polycrystal SMA are derived; (3) the quantitative incorporation of three microstructure factors (i.e., nucleation, growth and orientation distribution of martensite) into the modeling. These microstructural factors are intrinsic of specific polycrystal SMA systems and the role of each factor in the macroscopic constitutive response is quantitatively modeled. It is demonstrated that the interplay of these factors will result in different macroscopic transformation kinematics and kinetics which are responsible for the observed macroscopic stress-strain hardening or softening response, the latter will lead to the localization and propagation of transformation bands in TiNi SMA. The project supported by the Research Grant Committee (RGC) of Hong Kong SAR, the National Natural Science Foundation of China and the Provincial Natural Foundation of Jiangxi Province of China     

A CONSTITUTIVE MODEL FOR TRANSFORMATION, REORIENTATION AND PLASTIC DEFORMATION OF SHAPE MEMORY ALLOYS

A constitutive model is developed for the transformation, reorientation and plastic deformation of shape memory alloys (SMAs). It is based on the concept that an SMA is a mixture composed of austenite and martensite, the volume fraction of each phase is transformable with the change of applied thermal-mechanical loading, and the constitutive behavior of the SMA is the combination of the individual behavior of its two phases. The deformation of the martensite is separated into elastic, thermal, reorientation and plastic parts, and that of the austenite is separated into elastic, thermal and plastic parts. Making use of the Tanaka’s transformation rule modified by taking into account the effect of plastic deformation, the constitutive model of the SMA is obtained. The ferroelasticity, pseudoelasticity and shape memory effect of SMA Au-47.5 at.%Cd, and the pseudoelasticity and shape memory effect as well as plastic deformation and its effect of an NiTi SMA, are analyzed and compared with experimental results.     

A thermodynamical constitutive model for shape memory materials. Part II. The SMA composite material

The phenomenological SMA equations developed in Part I are used in this second paper to derive the free energy and dissipation of a SMA composite material. The derivation consists of solving a boundary value problem formulated over a mesoscale representative volume element, followed by an averaging procedure to obtain the macroscopic composite constitutive equations. Explicit equations are derived for the transformation tensors that relate the composite transformation strain rate to the phase transformation rate in the fiber and matrix. Some key findings for the two-way SME in a SMA fiber/elastomer matrix composite are that processing-induced residual stresses alter the composite austenite start and martensite start temperatures, as well as the amount of composite strain recovered during a complete cycle of temperature and fiber martensite volume fraction. Relative to the two-way SME response of stiff-matrix composites, it was found that compliant-matrix composites: (1) complete the phase transformation over a narrower temperature range; (2) exhibit greater transformation strain during the reverse transformation; and (3) undergo an incomplete strain cycle during a complete cycle of temperature and fiber martensite volume fraction. Due to the interaction of the fiber and matrix during transformation, macroscopic proportional stressing of the composite results in non-proportional fiber stressing, which in turn causes a small amount of martensitic reorientation to occur simultaneously with the transformation.     

CHARACTERISTICS OF STRESS-INDUCED TRANSFORMATION AND MICROSTRUCTURE EVOLUTION IN Cu-BASED SMA

The mechanical behavior of shape memory alloys （SMAs） is closely related to the formation and evolution of its microstructures. Through theoretical analysis and experimental observations, it was found that the stress-induced martensitic transformation process of single crystal Cu-based SMA under uniaxial tension condition consisted of three periods： nucleation, mixed nucleation and growth, and merging due to growth. During the nucleation, the stress dropped rapidly and the number of interfaces increased very fast while the phase fraction increased slowly. In the second period, both the stress and the interface number changed slightly but the phase fraction increased dramatically. Finally, the stress and the phase fraction changed slowly while the number of interfaces decreased quickly. Moreover, it was found that the transformation could be of multi-stage： sharp stress drops at several strains and correspondingly, the nucleation and growth process occurred quasi-independently in several parts of the sample.     

Mesoscale kinetics produces martensitic microstructure

We present molecular dynamics (MD) simulations of a martensitic phase transformation studying post-transformation microstructure and moving austenite-martensite interfaces. Unlike in energy-minimisation theories, the transformation dynamics dominate the martensite morphology. We use a binary Lennard-Jones potential to describe a square-to-hexagonal transformation by shear-and-shuffle. The high-T stable square lattice and low-T hexagonal lattice represent austenite and martensite, giving four martensitic variants. Compatible twin variants have no lattice misfit and zero interfacial energies which makes our model directly comparable with the crystallographic theory of martensite. Although our dynamical interpretation is different to previous work, our MD simulations exhibit very similar martensitic morphologies to real materials. We observe the nucleation of wedge-shaped, twinned martensite plates, plate growth at narrow, travelling transformation zones, subsonic transformation waves, elastic precursors inducing secondary nucleations and the formation of martensitic domains. Martensite is produced within narrow transformation zones where atoms change their lattice sites in a co-operative manner so as to form crystallographic layers. These motions produce inertia forces on the mesoscopic length-scale which induce the formation of twin variants in the subsequent layers to transform.     

Finite Element Analysis of Shape Memory Alloy Adaptive Trusses with Geometrical Nonlinearities

This contribution deals with the nonlinear analysis of shape memory alloy (SMA) adaptive trusses employing the finite element method. Geometrical nonlinearities are incorporated into the formulation together with a constitutive model that describes different thermomechanical behaviors of SMA. It has four macroscopic phases (three variants of martensite and an austenitic phase), and considers different material properties for austenitic and martensitic phases together with thermal expansion. An iterative numerical procedure based on the operator split technique is proposed in order to deal with the nonlinearities in the constitutive formulation. This procedure is introduced into ABAQUS as a user material routine. Numerical simulations are carried out illustrating the ability of the developed model to capture the general behavior of shape memory bars. After that, it is analyzed the behavior of some adaptive trusses built with SMA actuators subjected to different thermomechanical loadings.     

Fracture in single crystal NiTi

This paper presents an investigation of the pseudoelastic transformation of two crystallographic orientations of a single crystal NiTi shape memory alloy (SMA). Both uniaxial tension and notched tension samples are considered with the tensile axes along the [1 0 0] and [1 1 1] directions. The phase transformation is observed using optical techniques in situ. For the uniaxial tension samples, martensite plates and Lüders-band structures are observed. For the notched samples, transformation structures in the [1 0 0] sample appeared predominantly on the sides of the notch and crack tip with stable crack propagation. In the [1 1 1] notched samples transformation occurs directly ahead of the notch and unstable crack propagation is observed. An available work criterion is used to predict the location of the transformation, with good agreement to the experimental observations. The different fracture behavior of the two notched sample orientations is explained utilizing the available work calculations.     

Thermomechanical behavior of shape memory alloy under complex loading conditions

The thermo-mechanical behavior of polycrystalline shape memory alloy (SMA) under multi-axial loading with varying temperature conditions has been studied by experiments. Recently the research has been extended theoretically and a mechanical model of polycrystalline SMA and the corresponding mesoscopic constitutive equations have been developed. The model presented in this paper is constructed on the basis of the crystal plasticity and the deformation mechanism of SMA. The variants in the crystal grains and the orientations of crystal grains in the polycrystal are considered in the proposed model; the constitutive equations are derived on the basis of the proposed model. The volume fraction of the martensite variants in the transformation process and the influence of the stress state on the transformation process are also considered. Some calculated results obtained by the constitutive equations are presented and compared with the experimental results. It is found that the deformation behavior of SMA under complex loading conditions can be well reproduced by the calculation of the constitutive equations.     

Simulations of localized thermo-mechanical behavior in a NiTi shape memory alloy

Previous experiments have shown that stress-induced martensitic transformation in certain polycrystalline NiTi shape memory alloys can lead to strain localization and propagation phenomena when loaded in uniaxial tension. The number of nucleation events and kinetics of transformation fronts were found to be sensitive to the nature of the ambient media and imposed loading rate due to the release/absorption of latent heat and the material's inherent temperature sensitivity of the transformation stress. A special plasticity-based constitutive model used within a 3-D finite element framework has previously been shown to capture the isothermal, purely mechanical front features seen in experiments of thin uniaxial NiTi strips. This paper extends the approach to include the thermo-mechanical coupling of the material with its environment. The simulations successfully capture the nucleation and evolution of fronts and the corresponding temperature fields seen during the experiments.     

Crystal plasticity finite element modeling of mechanically induced martensitic transformation (MIMT) in metastable austenite

A new crystal plasticity model incorporating the mechanically induced martensitic transformation in metastable austenitic steel has been formulated and implemented into the finite element analysis. The kinetics of martensite transformation is modeled by taking into consideration of a nucleation-controlled phenomenon, where each potential martensitic variant based on Kurdjumov–Sachs (KS) relationship has different nucleation probability as a function of the interaction energy between externally applied stress and lattice deformation. Therefore, the transformed volume fractions are determined following selective variants given by the crystallographic orientation of austenitic matrix and applied stress in the frame of the crystal plasticity finite element. The developed finite element program is capable of considering the effect of volume change by the Bain deformation and the lattice-invariant shear during the martensitic transformation by effectively modifying the evolution of plastic deformation gradient of the conventional rate-dependent crystal plasticity finite element. The validation of the proposed model has been carried out by comparing with the experimentally measured data under simple loading conditions. Good agreements with the measurements for the stress–strain responses, transformed martensitic volume fractions and the influence of strain rate on the deformation behavior will enable the model to be promising for the future applications to the real forming process of the TRIP aided steel.     

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.