A constitutive model for shape memory alloys accounting for thermomechanical coupling

This paper presents a generalized Zaki-Moumni (ZM) model for shape memory alloys (SMAs) [cf. Zaki, W., Moumni, Z., 2007a. A three-dimensional model of the thermomechanical behavior of shape memory alloys. J. Mech. Phys. Solids 55, 2455-2490 accounting for thermomechanical coupling. To this end, the expression of the Helmholtz free energy is modified in order to derive the heat equation in accordance with the principles of thermodynamics. An algorithm is proposed to implement the coupled ZM model into a finite element code, which is then used to solve a thermomechanical boundary value problem involving a superelastic SMA structure. The model is validated against experimental data available in the literature. Strain rate dependence of the mechanical pseudoelastic response is taken into account with good qualitative as well as quantitative accuracy in the case of moderate strain rates and for mechanical results in the case of high strain rates. However, only qualitative agreement is achieved for thermal results at high strain rates. It is shown that this discrepancy is mainly due to localization effects which are note taken into account in our model. Analyzing the influence of the heat sources on the material response shows that the mechanical hysteresis is mainly due to intrinsic dissipation, whereas the thermal response is governed by latent heat. In addition, the variation of the area of the hysteresis loop with respect to the strain rate is discussed. It is found that this variation is not monotonic and reaches a maximum value for a certain value of strain rate.     

应力作用下NiTi形状记忆合金微结构演化的相场模拟及其本征应变率敏感性

基于Ginzburg-Landau动力学控制方程建立了NiTi形状记忆合金非等温相场模型,实现了对NiTi合金内应力诱导马氏体相变的数值模拟。同时将晶界能密度引入系统局部自由能密度,从而考虑多晶系统中晶界的重要作用。数值计算了单晶和多晶NiTi形状记忆合金在单轴机械载荷作用下微结构的动态演化过程和宏观力学行为,并重点研究了晶粒尺寸为60 nm的NiTi纳米多晶在低应变率下（0.0005～15 s?1）力学行为的本征应变率敏感性。研究结果表明,单晶NiTi合金系统高温拉伸-卸载过程中马氏体相变均匀发生,未形成奥氏体-马氏体界面。而纳米多晶系统在加载阶段出现了马氏体带的形成-扩展现象,在卸载阶段出现了马氏体带的收缩-消失现象。相同外载作用过程中,NiTi单晶系统的宏观应力-应变曲线具有更大的滞回环面积,拥有更优的超弹性变形能力。计算结果显示,在中低应变率下纳米晶NiTi形状记忆合金应力-应变关系表现出较明显的应变率相关性,应变率升高导致材料相变应力提升。这一应变率相关性主要源于相场模型中外加载荷速率与马氏体空间演化速度的相互竞争关系。     

NiTi形状记忆合金裂纹尖端热力耦合行为研究

本文对NiTi形状记忆合金I型裂纹尖端热力耦合行为进行了数值仿真分析和实验验证。建立了包含相变和热力耦合的本构模型,通过有限元计算得到了裂纹尖端附近的纵向应变、马氏体体积分数和温度场分布,依据马氏体相变情况对裂纹尖端有效应力强度因子进行了修正,揭示了加载速率对形状记忆合金裂纹尖端有效应力强度影子的影响规律。参数研究表明,随着加载频率的增加,裂纹尖端附近温度逐渐升高,马氏体相变区域逐渐缩小,有效应力强度因子呈下降趋势,形状记忆合金表现出增韧效应,有助于减缓裂纹扩展。本研究结果对于揭示热力耦合作用下超弹性形状记忆合金疲劳裂纹扩展规律具有重要参考意义。     

Hysteresis in shape-memory alloys

We present an overview of hysteresis phenomena in the martensitic transformation, and their relevance in the thermomechanical behaviour of shape-memory alloys. The first part of the paper introduces the concept of hysteresis, and the related phenomena of branching, dissipation and memory. The second part deals with revising some aspects of the thermomechanical behaviour of shape-memory alloys, emphasizing the hysteretic behaviour of single crystals and polycrystals under different driving conditions. The last part of the work is dedicated to the problem of modelling hysteresis phenomena in shape-memory alloys. Our focus is on phenomenological approaches which, as shown in the paper, account for the memory properties observed in hysteretic trajectories and open the possibility of deriving a generic energy balance for systems with hysteresis.     

Cyclic thermomechanical behavior of a polycrystalline pseudoelastic shape memory alloy

In this work, 3D finite element modeling is employed to examine the thermomechanical behavior of a polycrystalline Ni-Ti shape memory alloy in the pseudoelastic regime. It is shown that the tension-compression asymmetry during uniaxial cyclic loading is due to a preferred orientation of the crystallographic texture. In pure shear loading, the thermomechanical behavior exhibits symmetry in both senses of shear, due to the fiber texture of the specimen bar stock. It is also shown that the apparent strain rate-dependence is due to thermomechanical coupling with latent heat generation/absorption during phase transformation.     

On modelling phase propagation in SMAs by a Maxwellian thermo-viscoelastic approach

In this paper we propose a Maxwellian thermo-viscoelastic approach to the problem of phase transformation in shape memory alloys. An explicit temperature-dependent non-monotone piecewise linear stress–strain relation is considered and the corresponding free energy function is used to establish the heat propagation equation. The heat exchange between the material and its environment is also taken into account. The numerical simulations of three end-displacement rates lead to serrated hysteresis loops and result in inhomogeneous deformations and exothermic/endothermic behavior during loading/unloading tests. It is shown that the influence of the strain rate on the size and shape of the hysteresis loop is due in fact to the heat exchange between the bar and its environment. The predictions are compared qualitatively with experimental results.     

A 3-D constitutive model for shape memory alloys incorporating pseudoelasticity and detwinning of self-accommodated martensite

A 3-D constitutive model for polycrystalline shape memory alloys (SMAs), based on a modified phase transformation diagram, is presented. The model takes into account both direct conversion of austenite into detwinned martensite as well as the detwinning of self-accommodated martensite. This model is suitable for performing numerical simulations on SMA materials undergoing complex thermomechanical loading paths in stress–temperature space. The model is based on thermodynamic potentials and utilizes three internal variables to predict the phase transformation and detwinning of martensite in polycrystalline SMAs. Complementing the theoretical developments, experimental data are presented showing that the phase transformation temperatures for the self-accommodated martensite to austenite and detwinned martensite to austenite transformations are different. Determination of some of the SMA material parameters from such experimental data is also discussed. The paper concludes with several numerical examples of boundary value problems with complex thermomechanical loading paths which demonstrate the capabilities of the model.     

Thermomechanical modeling of nonlinear internal hysteresis due to incomplete phase transformation in pseudoelastic shape memory alloys

This paper presents a thermomechanical model for pseudoelastic shape memory alloys (SMAs) accounting for internal hysteresis effect due to incomplete phase transformation. The model is developed within the finite-strain framework, wherein the deformation gradient is multiplicatively decomposed into thermal dilation, rigid body rotation, elastic and transformation parts. Helmholtz free energy density comprises three components: the reversible thermodynamic process , the irreversible thermodynamic process and the physical constraints of both. In order to capture the multiple internal hysteresis loops in SMA, two internal variables representing the transition points of the forward and reverse phase transformation, \(\phi _s^f\) and \(\phi _s^r\), are introduced to describe the incomplete phase transformation process. Evolution equations of the internal variables are derived and linked to the phase transformation. Numerical implementation of the model features an Euler discretization and a cutting-plane algorithm. After validation of the model against the experimental data, numerical examples are presented, involving a SMA-based vibration system and a crack SMA specimen subjected to partial loading–unloading case. Simulation results well demonstrate the internal hysteresis and free vibration behavior of SMA.

     

Phase proportioning in CuAlBe shape memory alloys during thermomechanical loadings using electric resistance variation

The microstructure of shape memory alloys changes with the thermomechanical history of the material. During thermomechanical loading, austenite, thermally-induced martensite or stress-induced martensite can be simultaneously present in the material. In applications integrating SMA parts, utilization conditions seriously affect the microstructure and can generate macroscopic strain or stress. Consequently, during thermomechanical loadings, it is important to be able to proportion the different phases and consequently to understand the kinetic transformation. This is very useful in the development of constitutive equations. This study shows, by a series of tests, that the proposed experimental method, based on the measurement of the variation of electric resistance of CuAlBe wires, permits to determine the volume fraction of the different phases present in the material (i.e., austenite, stress-induced martensite and thermally-induced martensite). The proposed method is applied to the most common thermomechanical behavior met in engineering applications of shape memory alloys: pseudoelasticity, pseudoplasticity, recovery-stress and stress-assisted two-way shape memory effect. The proportioning method based on a mixture law integrating the resistivity of pure phases present in the SMA is first performed on different two-phase mixture cases and then applied to a three phase mixture case.     

Thermomechanical cyclic behavior modeling of Cu-Al-Be SMA materials and structures

This paper concerns the behavior of Cu-Al-Be polycrystalline shape memory alloys under cyclic thermomechanical loadings. Sometimes, as shown by many experimental observations, a permanent inelastic strain occurs and increases with the number of cycles. A series of cyclic thermomechanical tests has been carried out and the origin of the residual strain has been identified as residual martensite. These observations have been used to develop a 3D macroscopic model for the superelasticity and stress assisted memory effect of SMAs able to describe the evolution of permanent inelastic strain during cycles. The model has been implemented in a finite elements code and used to simulate the behavior of antagonistic actuators based on SMA springs under cyclic thermomechanical loading with a residual displacement appearance.     

A three-dimensional model of the thermomechanical behavior of shape memory alloys

A new macro-scale model of shape memory alloys is developed within the framework of generalized standard materials with internal constraints [Moumni, Z., 1995. Sur la modélisation du changement de phase à l’état solide. Ph.D. Thesis, École Nationale Supérieure des Ponts et Chaussées]. It is shown that the introduction of two state variables: the martensite volume fraction and the martensite orientation strain tensor, is sufficient to account for several effects exhibited by SMAs subject to thermomechanical loading, namely: self-accommodation, orientation and reorientation or martensite, as well as superelasticity and one-way shape memory. These phenomena are simulated using the same set of constitutive equations and evolution laws. A simple procedure for identifying the parameters of the model is described in detail and a validation against experimental data is conducted. The model is then used to analyze a 3D SMA structure representing a superelastic stent. Cyclic loading and other pertaining phenomena like training and two-way shape memory are considered in the second part of this paper.     

On the thermomechanical modeling of shape memory alloys

In this work, a general inelastic framework for the derivation of general three-dimensional thermomechanical constitutive laws for materials undergoing phase transformations is proposed. The proposed framework is based on the generalized plasticity theory and on some basic elements from the theory of continuum damage mechanics. More specifically, a new elaborate formulation of generalized plasticity theory capable of accommodating the multiple and interacting loading mechanisms, which occur during the phase transformations, is developed. Furthermore, the stiffness variations occurring during phase transformations are taken into account by the proposed framework. For this purpose, the free energy is decomposed into elastic and inelastic parts, not in a conventional way, but in one which resembles the elastic-damage cases. Also, a rate-dependent version of the theory is provided. The concepts presented are applied for the derivation of a three-dimensional thermomechanical constitutive model for Shape Memory Alloy materials. Numerical simulations to show qualitatively the ability of the model to capture the behavior of the shape memory alloys are also presented. Furthermore, the model has been fitted to actual experimental results from the literature.     

超弹性NiTi合金丝动力特性试验及本构模型研究

形状记忆合金(SMA)是一种兼具感知和驱动功能的功能材料,因其独特的形状记忆效应、超弹性和高阻尼等特性,成为土木工程结构振动控制的理想材料.论文研究了超弹性NiTi丝的动力特性和应变率相关的本构模型.试验测试了NiTi丝在不同应变率下的力学性能,建立了应力增量与应变率的关系方程.在试验的基础上,提出了改进的SMA本构模...     

The influence of model parameters and of the thermomechanical coupling on the behavior of shape memory devices

Shape memory materials (SMM) are receiving increasing attention for their use in applications that exploit their dynamic behavior. A thermomechanical model for devices with pseudoelastic behavior has been proposed in previous works [11] (Bernardini and Pence, 2005) [15] (Bernardini and Rega, 2005). The model takes into account several aspects of SMM behavior by means of seven model parameters.In this paper the effect of each parameter on the non-isothermal rate-dependent behavior of the device is studied, by paying particular attention to the effect of the thermomechanical coupling. Some overall synthetic indicators of the behavior of the shape memory device are defined in terms of the model parameters. By evaluating such indicators, a lot of information about the mechanical, thermal and thermomechanical effects on the device behavior can be gained before computing explicitly the response of the shape memory oscillator.The present work may provide a guide for the proper utilization of the model for the investigation of the dynamic response.     

Coupled thermomechanical dynamics of phase transitions in shape memory alloys and related hysteresis phenomena

In this paper the nonlinear dynamics of shape memory alloy phase transformations is studied with thermomechanical models based on coupled systems of partial differential equations by using computer algebra tools. The reduction procedures of the original model to a system of differential-algebraic equations and its solution are based on the general methodology developed by the authors for the analysis of phase transformations in shape memory materials with low dimensional approximations derived from center manifold theory. Results of computational experiments revealing the martensitic-austenitic phase transition mechanism in a shape-memory-alloy rod are presented. Several groups of computational experiments are reported. They include results on stress- and temperature-induced phase transformations as well as the analysis of the hysteresis phenomenon. All computational experiments are presented for Cu-based structures.     

Shape memory behaviour: modelling within continuum thermomechanics

A phenomenological material model to represent the multiaxial material behaviour of shape memory alloys is proposed. The material model is able to represent the main effects of shape memory alloys: the one-way shape memory effect, the two-way shape memory effect due to external loads, the pseudoelastic and pseudoplastic behaviour as well as the transition range between pseudoelasticity and pseudoplasticity.The material model is based on a free energy function and evolution equations for internal variables. By means of the free energy function, the energy storage during thermomechanical processes is described. Evolution equations for internal variables, e.g. the inelastic strain tensor or the fraction of martensite are formulated to represent the dissipative material behaviour. In order to distinguish between different deformation mechanisms, case distinctions are introduced into the evolution equations. Thermomechanical consistency is ensured in the sense that the constitutive model satisfies the Clausius–Duhem inequality.Finally, some numerical solutions of the constitutive equations for isothermal and non-isothermal strain and stress processes demonstrate that the various phenomena of the material behaviour are well represented. This applies for uniaxial processes and for non-proportional loadings as well.     

Differential and integrated form consistency in 1-D phenomenological models for shape memory alloy constitutive behavior

Shape memory alloys are being explored increasingly for developing smart structures and devices in aerospace, automotive and other application areas. The material behavior is highly nonlinear with coupled thermomechanical response involving temperature and/or stress induced phase transformations. Modeling the constitutive behavior of these materials poses several challenges and a few phenomenological models exist that provide a quick and reasonable approach to assess their behavior. Due to phenomenological approach, several assumptions are made in order to simplify the model and some of them introduce inconsistencies or anomalies into the model. In this paper, a frequently used approach, namely, Brinson [Brinson, L.C., 1993. One dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable. J. Intell. Mater. Syst. Struct. 4(2), 229–242.] model, is investigated. The constitutive equation is usually expressed at the outset in the differential form and the integrated form of the same is obtained. It is shown that the two forms of equations are not consistent in the Brinson [Brinson, L.C., 1993. One dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable. J. Intell. Mater. Syst. Struct. 4(2), 229–242.] model, given the assumed form of material functions. In the present work, the nature and implications of the inconsistency are highlighted. The cause of incompatibility is the inconsistent material definitions. A modified consistent constitutive model is proposed by redefining the material function which satisfies the compatibility condition. The advantages in using the proposed modified model are highlighted with numerical case studies involving hysteretic stress–strain behavior.