A finite element model for shape memory alloys considering thermomechanical couplings at large strains

In the present work we propose a new thermomechanically coupled material model for shape memory alloys (SMA) which describes two important phenomena typical for the material behaviour of shape memory alloys: pseudoelasticity as well as the shape memory effect. The constitutive equations are derived in the framework of large strains since the martensitic phase transformation involves inelastic deformations up to 8%, or even up to 20% if the plastic deformation after the phase transformation is taken into account. Therefore, we apply a multiplicative split of the deformation gradient into elastic and inelastic parts, the latter concerning the martensitic phase transformation. An extended phase transformation function has been considered to include the tension–compression asymmetry particularly typical for textured SMA samples. In order to apply the concept in the simulation of complex structures, it is implemented into a finite element code. This implementation is based on an innovative integration scheme for the existing evolution equations and a monolithic solution algorithm for the coupled mechanical and thermal fields. The coupling effect is accurately investigated in several numerical examples including pseudoelasticity as well as the free and the suppressed shape memory effect. Finally, the model is used to simulate the shape memory effect in a medical foot staple which interacts with a bone segment.     

形状记忆合金拟弹性行为的热力学描述

形状记忆合金是由马氏体和奥氏体组成并动态变化的两相材料 ,其拟弹性行为实质上是两相各自行为的动态组合。本文提出了形状记忆合金拟弹性行为的一种热力学描述。根据实验现象假设在感兴趣的温度和变形范围内 ,奥氏体相具有线弹性特性而马氏体相具有弹塑性特性 ,结合 Tanaka的相变描述 ,给出了小变形、初始各向同性和塑性不可压缩条件下形状记忆合金的三维本构方程。对不同温度下形状记忆合金材料的特性进行了描述 ,较好地预言了单调及循环加载下的响应和正、反相变行为及其温度影响 ,动态相变过程对应力响应的影响 ,高温相下的强度增加等。     

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.     

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.     

纳米单晶NiTi合金单程形状记忆效应的分子动力学模拟

采用基于第二近邻修正型嵌入原子势的分子动力学方法研究了纳米单晶NiTi合金的单程形状记忆效应,详细阐明了温度诱发马氏体相变和应力诱发马氏体重定向过程中纳米单晶的变形行为和微结构演化,进一步分析了加/卸载速率对NiTi合金单程形状记忆效应的影响。结果表明,NiTi纳米单晶在应力加载过程中发生马氏体重定向,卸载后存在残余应变;当加热到奥氏体转变结束温度以上时,马氏体逆相变为奥氏体相,残余应变逐渐减小,但未完全回复;随着应力加载速率的增加,重定向临界应力和模量逐渐增加;再次降温过程中不同加载速率下的原子结构演化各不相同。     

纳米单晶NiTi合金单程形状记忆效应的分子动力学模拟

采用基于第二近邻修正型嵌入原子势的分子动力学方法研究了纳米单晶NiTi合金的单程形状记忆效应,详细阐明了温度诱发马氏体相变和应力诱发马氏体重定向过程中纳米单晶的变形行为和微结构演化,进一步分析了加/卸载速率对NiTi合金单程形状记忆效应的影响。结果表明,NiTi纳米单晶在应力加载过程中发生马氏体重定向,卸载后存在残余应变;当加热到奥氏体转变结束温度以上时,马氏体逆相变为奥氏体相,残余应变逐渐减小,但未完全回复;随着应力加载速率的增加,重定向临界应力和模量逐渐增加;再次降温过程中不同加载速率下的原子结构演化各不相同。     

形状记忆合金的一种基于经典塑性理论的两相混合本构模型

形状记忆合金由马氏体相和奥氏体相动态组成,其行为实质上是两相各自行为的动态组合,根据实验现象,假设在一定的变形范围内,马氏体相为弹塑性而奥氏体相为线弹性,基于经典塑性理论和混合物理论,结合Tanaka的相变描述,得到了形状记忆合金的一种本构描述,对不同温度下形状记忆合金Au-47．5at．％Cd的铁弹性、拟弹性和形状记忆特性进行了分析,取得了与实验相吻合的结果。     

A 3D finite strain phenomenological constitutive model for shape memory alloys considering martensite reorientation

Most devices based on shape memory alloys experience both finite deformations and non-proportional loading conditions in engineering applications. This motivates the development of constitutive models considering finite strain as well as martensite variant reorientation. To this end, in the present article, based on the principles of continuum thermodynamics with internal variables, a three-dimensional finite strain phenomenological constitutive model is proposed taking its basis from the recent model in the small strain regime proposed by Panico and Brinson (J Mech Phys Solids 55:2491–2511, 2007). In the finite strain constitutive model derivation, a multiplicative decomposition of the deformation gradient into elastic and inelastic parts, together with an additive decomposition of the inelastic strain rate tensor into transformation and reorientation parts is adopted. Moreover, it is shown that, when linearized, the proposed model reduces exactly to the original small strain model.     

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.     

Crack growth resistance of shape memory alloys by means of a cohesive zone model

Crack growth resistance of shape memory alloys (SMAs) is dominated by the transformation zone in the vicinity of the crack tip. In this study, the transformation toughening behavior of a slowly propagating crack in an SMA under plane strain conditions and mode I deformation is numerically investigated. A small-scale transformation zone is assumed. A cohesive zone model is implemented to simulate crack growth within a finite element scheme. Resistance curves are obtained for a range of parameters that specify the cohesive traction-separation constitutive law. It is found that the choice of the cohesive strength t₀ has a great influence on the toughening behavior of the material. Moreover, the reversibility of the transformation can significantly reduce the toughening of the alloy. The shape of the initial transformation zone, as well as that of a growing crack is determined. The effect of the Young's moduli ratio of the martensite and austenite phases is examined.     

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.     

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.     

CONSTITUTIVE MODELING OF ROLLED SHAPE MEMORY ALLOY SHEETS TAKING INTO ACCOUNT PRE-TEXTURE AND ANISOTROPIC HARDENING

The drawing or rolling process endows polycrystal shape memory alloy with a crys- tallographic texture, which can result in macroscopic anisotropy. The main purpose of this work is to develop a constitutive model to predict the thermomechanical behavior of shape memory alloy sheets, which accounts for the crystallographic texture. The total macroscopic strain is decom- posed into elastic strain and macro-transformation strain under isothermal condition. Considering the transformation strain in local grains and the orientation distribution function of crystallo- graphic texture, the macro-transformation strain and the effective elastic modulus of textured polycrystal shape memory alloy are developed by using tensor expressions. The kinetic equation is established to calculate the volume fraction of the martensite transformation under given stress. Furthermore, the Hill＇s quadratic model is developed for anisotropic transformation hardening of textured SMA sheets. All the calculation results are in good agreement with experimental data, which show that the present model can accurately describe the macro-anisotropic behaviors of textured shape memory alloy sheets.     

??????????????????????????????

基于形状记忆合金Brinson一维热力学本构关系和von K\'{a}rm\'{a}n几何非线性薄板理论,研究了径向嵌入SMA丝复合材料加热圆板在横向均布 机械载荷作用下的弯曲响应, 获得了周边不可移简支和夹紧圆板的中心最大挠度与升温之间的关系曲线. 结果表明,形状 记忆合金丝在从马氏体向奥氏体的逆相变过程中所产生的相变回复力对板的弯曲变形具有明 显的调整作用. 通过嵌入SMA纤维丝和施加升温载荷可以主动而有效地调节受机 械载荷作用圆板的弯曲变形.     

A NEW MODEL OF SHAPE MEMORY ALLOYS

A new constitutive model of shape memory alloys ( SMAs) based on Tanaka ' s martensite fraction exponential expression is produced. This new model can present recoverable shape memory strain during different phase transformation, and reflect the action of martensite reorientation . Also it can overcome the defect of Tanaka ' s Model when the SMAs ' microstructure is fully martensite . The model is very simple and suitable for using , and the correct behavior of the model is proved by test.     

Micromechanics of plastic deformation and phase transformation in a three-phase TRIP-assisted advanced high strength steel: Experiments and modeling

The micromechanics of plastic deformation and phase transformation in a three-phase advanced high strength steel are analyzed both experimentally and by microstructure-based simulations. The steel examined is a three-phase (ferrite, martensite and retained austenite) quenched and partitioned sheet steel with a tensile strength of ~980 MPa. The macroscopic flow behavior and the volume fraction of martensite resulting from the austenite–martensite transformation during deformation were measured. In addition, micropillar compression specimens were extracted from the individual ferrite grains and the martensite particles, and using a flat-punch nanoindenter, stress–strain curves were obtained. Finite element simulations idealize the microstructure as a composite that contains ferrite, martensite and retained austenite. All three phases are discretely modeled using appropriate crystal plasticity based constitutive relations. Material parameters for ferrite and martensite are determined by fitting numerical predictions to the micropillar data. The constitutive relation for retained austenite takes into account contributions to the strain rate from the austenite–martensite transformation, as well as slip in both the untransformed austenite and product martensite. Parameters for the retained austenite are then determined by fitting the predicted flow stress and transformed austenite volume fraction in a 3D microstructure to experimental measurements. Simulations are used to probe the role of the retained austenite in controlling the strain hardening behavior as well as internal stress and strain distributions in the microstructure.     

Shakedown analysis of shape memory alloy structures

Phase transformational shakedown of a structure refers to a status that plastic strains cease developing after a finite number of loading cycles, and subsequently the structure undergoes only elastic deformation and alternating phase transformations with limited magnitudes. Due to the intrinsic complexity in the constitutive relations of shape memory alloys (SMA), there is as yet a lack of effective methods for modeling the mechanical responses of SMA structures, especially when they develop both phase transformation and plastic deformation. This paper is devoted to present an algorithm for analyzing shakedown of SMA structures subjected to cyclic or varying loads within specified domains. Based on the phase transformation and plastic yield criteria of von Mises-type and their associated flow rules, a simplified three-dimensional phenomenological constitutive model is first formulated accounting for different regimes of elastic–plastic deformation and phase transformation. Different responses possible for SMA bodies exposed to varying loads are discussed. The classical Melan shakedown theorem is extended to determine a lower bound of loads for transformational shakedown of SMA bodies without necessity of a step-by-step analysis along the loading history. Finally, a simple example is given to illustrate the application of the present theory as well as some basic features of shakedown of SMA structures. It is interesting to find that phase transformation may either increase or decrease the load-bearing capacity of a structure, depending upon its constitutive relations, geometries and the loading mode.     

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.     

Time integration of a model for martensite detwinning and reorientation under nonproportional loading using Lagrange multipliers

The paper presents a consistent integration scheme for a model of shape memory alloys that takes into account the processes of detwinning and reorientation of martensite. Lagrange multipliers are used to account for the saturation condition on the martensitic phase, which expresses the requirement on the inelastic deformation of martensite not to exceed a material-specific value. A detailed procedure for implicit time integration of the constitutive equations is given, including a comprehensive derivation of a closed-form expression for the consistent tangent moduli. The equations are implemented into a commercial finite element software, which is used to simulate the response of martensitic SMA structures to complex loading. The results are shown to agree with experimental data taken from the literature.