超弹性镍钛形状记忆合金单轴相变棘轮行为的宏观唯象本构模型

超弹性镍钛形状记忆合金因其良好的力学性能以及独特的超弹性和形状记忆效应已广泛应用于土木工程、航空航天和生物医疗等多个领域,在实际服役环境中超弹性镍钛合金元件不可避免地会承受不同应力水平的循环载荷作用,亟待建立描述相变棘轮行为(即峰值应变和谷值应变随着正相变和逆相变循环的进行不断累积)的循环本构模型.为此,基于已有的超弹性镍钛形状记忆合金在不同峰值应力下的单轴相变棘轮行为实验研究结果,在广义黏塑性框架下,对Graesser等提出的通过背应力非线性演化方程反映超弹性镍钛形状记忆合金超弹性行为的一维宏观唯像本构模型进行了拓展,考虑了正相变和逆相变过程中特征变量的差异及其随循环的演化,以非弹性应变的累积量为内变量引入了正相变开始应力、逆相变开始应力、相变应变和残余应变的演化方程,同时通过峰值应力与正相变完成应力的比值来确定演化方程中的相关系数,建立了描述超弹性镍钛合金单轴相变棘轮行为的本构模型.将模拟结果与对应的实验结果进行对比发现,建立的宏观唯像本构模型能够合理地描述超弹性镍钛形状记忆合金的单轴相变棘轮行为及其峰值应力依赖性,模型的预测结果和实验结果吻合得很好.     

ADVANCES IN TRANSFORMATION RATCHETING AND RATCHETING-FATIGUE INTERACTION OF NITI SHAPE MEMORY ALLOY

The accumulation of inelastic deformation occurring in NiTi shape memory alloy under the stress-controlled cyclic loading condition is named transformation ratcheting, since it is mainly caused by the solid-solid transformation from austenite to martensite phase and vice versa. The transformation ratcheting and its effect on the fatigue life (i.e., transformation-fatigue interaction) are key issues that should be addressed in order to assess the fatigue of NiTi shape memory alloy more accurately. In this paper, the advances in the studies on the transformation ratcheting and ratcheting-fatigue interaction of super-elastic NiTi shape memory alloy in recent years are reviewed: First, experimental observation of the uniaxial transformation ratcheting and ratcheting-fatigue interaction of super-elastic NiTi alloy under the stress-controlled cyclic loading conditions is treated, and the detrimental effect of transformation ratcheting on the fatigue life is addressed; Secondly, two types of cyclic constitutive models (i.e., a macroscopic phenomenological model and a micromechanical one based on crystal plasticity) constructed to describe the transformation ratcheting of super-elastic NiTi alloy are discussed; Furthermore, an energy-based failure model is provided and dealt with by comparing its predicted fatigue lives with experimental ones; Finally, some suggestions about future work are made.     

Study on the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy based on a new crystal plasticity constitutive model

In this paper, a crystal plasticity based constitutive model (Yu et al., 2013) is extended to describe the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy by considering the internal heat production. Two sources of internal heat productions are included in the proposed model, i.e., the mechanical dissipations of inelastic deformation and the transformation latent heat in the NiTi shape memory alloy. With an assumption of uniform temperature field in the alloy specimen, a simplified evolution law of temperature field is obtained by the first law of thermodynamics and the heat boundary conditions. An explicit scale-transition rule is adopted to extend the proposed single crystal model to the polycrystalline version. The capability of the extended polycrystalline model to describe the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy is verified by comparing the predictions with the corresponding experimental ones. The comparison demonstrates that the proposed constitutive model considering the internal heat production predicts the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy fairly well.     

A macroscopic multi-mechanism based constitutive model for the thermo-mechanical cyclic degeneration of shape memory effect of NiTi shape memory alloy

A macroscopic based multi-mechanism constitutive model is constructed in the framework of irreversible thermodynamics to describe the degeneration of shape memory effect occurring in the thermo-mechanical cyclic deformation of NiTi shape memory alloys(SMAs). Three phases,austenite A, twinned martensite Mtand detwinned martensite M~d, as well as the phase transitions occurring between each pair of phases( A → M~t, M~t→ A, A → M~d,M~d→ A, and M~t→ M~d) are considered in the proposed model. Meanwhile, two kinds of inelastic deformation mechanisms, martensite transformation-induced plasticity and reorientation-induced plasticity, are used to explain the degeneration of shape memory effects of NiTi SMAs. The evolution equations of internal variables are proposed by attributing the degeneration of shape memory effect to the interaction between the three phases(A, M~t, and M~d) and plastic deformation. Finally, the capability of the proposed model is verified by comparing the predictions with the experimental results of NiTi SMAs. It is shown that the degeneration of shape memory effect and its dependence on the loading level can be reasonably described by the proposed model.     

梯度纳米晶NiTi形状记忆合金的超弹性和形状记忆效应相场模拟

通过建立考虑两个马氏体变体的二维相场模型,对梯度纳米晶镍钛(NiTi)合金系统的超弹性、单程和应力辅助双程形状记忆过程进行了模拟和预测.模拟结果显示: 在梯度纳米晶NiTi合金的超弹性过程中,较粗晶粒的区域保留了传统粗晶的马氏体相变和逆相变特征,即局部马氏体带的形核-扩展和缩减-消失, 而随着晶粒尺寸的减小,细晶粒区域表现为均匀相变的特点, 即无局部马氏体带产生; 此外,在超弹性和形状记忆过程中,马氏体相变和重取向都首先在较粗晶粒区域开始并逐步向细晶粒区域传播,而逆相变则相反.马氏体相变和重取向的逐步扩展使梯度纳米晶NiTi合金的应力-应变和应变-温度曲线呈现出“硬化状”,其可归因于纳米多晶NiTi合金中马氏体相变对晶粒尺寸的依赖性,即随着晶粒尺寸的减小, 相变或重取向壁垒逐渐增大,马氏体相变或重取向的形核、扩展越来越困难. 可见,梯度纳米晶结构具有比传统均匀晶粒尺寸NiTi合金更宽的相变应力区间、重取向应力区间和相变温度区间,可显著提高该合金非弹性变形的可控性.      

SUPER-ELASTIC CONSTITUTIVE MODEL CONSIDERING PLASTICITY AND ITS FINITE ELEMENT IMPLEMENTATION

Based on the experimental results of super-elastic NiTi alloy, a three-dimensional super-elastic constitutive model including both of stress-induced martensite transformation and plasticity is constructed in a framework of general inelasticity. In the proposed model, transformation hardening, reverse transformation of stress-induced martensite, elastic mismatch between the austenite and martensite phases, and temperature-dependence of transformation stress and elastic modulus of each phase are considered. The plastic yielding of martensite occurred under high stress is addressed by a bilinear isotropic hardening rule. Drucker-Prager-typed transformation surfaces are employed to describe the asymmetric behavior of NiTi alloy in tension and compression. The prediction capability of the proposed model is verified by comparing the simulated results with the correspondent experimental ones. Based on backward Euler＇s integration, a new expression of consistent tangent modulus is derived. The proposed model is then implemented into a finite element package ABAQUS by user-subroutine UMAT. Finally, the validity of such implementation was verified by some numerical samples.     

Micromechanical modelling of the effect of plastic deformation on the mechanical behaviour in pseudoelastic shape memory alloys

Except for the recoverable strain induced by phase transformation, NiTi alloys are very ductile even in the martensite phase. The purpose of the present paper is to study the influence of permanent deformation, which results from plastic deformation of martensite, on the mechanical behaviour of pseudoelastic NiTi alloys. Based on phenomenological theory of martensitic transformation and crystal plasticity, a new three dimensional micromechanical model is proposed by coupling both the slip and twinning deformation mechanisms. The present model is implemented as User MATerial subroutine (UMAT) into ABAQUS/Standard to study the influences of plastic deformation on the stress and strain fields, and on the evolution of martensite transformation. Results show that with the increasing of plastic deformation the residual strain increases and the phase transformation stress–strain curves from the martensite to austenite become steeper and less obvious. Both characteristics, stabilisation of martensite and impedance of the reverse transformation, due to plastic deformation are captured.     

纳米多晶NiTi形状记忆合金超弹性的晶粒取向依赖性相场研究

本文基于Ginzburg-Landau理论,建立了一个反映纳米多晶NiTi形状记忆合金取向依赖性的二维多晶相场模型,研究了晶粒取向对其超弹性性能的影响。结果表明,纳米多晶NiTi形状记忆合金的超弹性行为依赖于晶粒取向分布,即：多晶模型中在所研究的参数变化范围内,晶粒取向分布范围越广、晶粒间取向差越大（无明显织构）,超弹性性能越差;而晶粒取向分布范围越窄、晶粒间取向差越小（织构越强）,超弹性性能越好。该晶粒取向依赖性可解释为：由于晶粒取向的不同,马氏体相变过程中相邻晶粒间的变形不匹配程度不同,因此,多晶模型中在所研究的参数变化范围内,晶粒间取向差异越大,晶界处的变形失配越严重,由此而产生的局部内应力将阻碍其附近马氏体相变的扩展,进而导致纳米多晶NiTi形状记忆合金超弹性性能下降。     

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.     

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

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

Macro versus micro-scale constitutive models in simulating proportional and nonproportional cyclic and ratcheting responses of stainless steel 304

A recent study by Hassan et al. [Hassan, T., Taleb, L., Krishna, S., 2008. Influences of nonproportional loading paths on ratcheting responses and simulations by two recent cyclic plasticity models. Int. J. Plasticity, 24, 1863–1889.] demonstrated that some of the nonproportional ratcheting responses under stress-controlled loading histories cannot be simulated reasonably by two recent cyclic plasticity models. Two major drawbacks of the models identified were: (i) the stainless steel 304 demonstrated cyclic hardening under strain-controlled loading whereas cyclic softening under stress-controlled loading, which depends on the strain-range and which the existing models cannot describe; (ii) the change in biaxial ratcheting responses due to the change in the degree of nonproportionality were not simulated well by the models. Motivated by these findings, two modified cyclic plasticity models are evaluated in predicting a broad set of cyclic and ratcheting response of stainless steel 304. The experimental responses used in evaluating the modified models included both proportional (uniaxial) and nonproportional (biaxial) loading responses from Hassan and Kyriakides [Hassan, T., Kyriakides, S., 1994a. Ratcheting of cyclically hardening and softening materials. Part I: uniaxial behavior. Int. J. Plasticity, 10, 149–184; Hassan, T., Kyriakides, S., 1994b. Ratcheting of cyclically hardening and softening materials. Part II: multiaxial behavior. Int. J. Plasticity, 10, 185–212.] and Hassan et al. [Hassan, T., Taleb, L., Krishna, S., 2008. Influences of nonproportional loading paths on ratcheting responses and simulations by two recent cyclic plasticity models. Int. J. Plasticity, 24, 1863–1889.] The first model studied is a macro-scale, phenomenological, constitutive model originally proposed by Chaboche et al. [Chaboche, J.L., Dang-Van, K., Cordier, G., 1979. Modelization of the strain memory effect on the cyclic hardening of 316 stainless steel. In: Proceedings of the Fifth International Conference on SMiRT, Div. L, Berlin, Germany, L11/3.]. This model was systematically modified for incorporating strain-range dependent cyclic hardening–softening, and proportional and nonproportional loading memory parameters. The second model evaluated is a polycrystalline model originally proposed by Cailletaud [Cailletaud, G., 1992. A micromechanical approach to inelastic behavior of metals. Int. J. Plasticity, 8, 55–73.] based on crystalline slip mechanisms. These two models are scrutinized against simulating hysteresis loop shape, cyclic hardening–softening, cross-effect, cyclic relaxation, subsequent cyclic softening and finally a broad set of ratcheting responses under uniaxial and biaxial loading histories. The modeling features which improved simulations for these responses are elaborated in the paper. In addition, a novel technique for simulating both the monotonic and cyclic responses with one set of model parameters is developed and validated.     

Transformation-induced plasticity in high-temperature shape memory alloys: a one-dimensional continuum model

A constitutive model based on isotropic plasticity consideration ispresented in this work to model the thermo-mechanical behavior ofhigh-temperature shape memory alloys. In high-temperature shapememory alloys (HTSMAs), both martensitic transformation andrate-dependent plasticity (creep) occur simultaneously at hightemperatures. Furthermore, transformation-induced plasticity isanother deformation mechanism during martensitic transformation. Allthese phenomena are considered as dissipative processes to model themechanical behavior of HTSMAs in this study. The constitutive modelwas implemented for one-dimensional cases, and the results have beencompared with experimental data from thermal cycling test foractuator applications.     

Implicit finite element formulations for multi-phase transformation in high carbon steel

An anomalous plastic deformation observed during the phase transformation of steels was implemented into the finite element modeling. The constitutive equations for the transformation plasticity originally proposed by Greenwood and Johnson [Greenwood, G.W., Johnson, R.H., 1965. The deformation of metals under small stresses during phase transformation. Proc. Roy. Soc. A 283, 403] and further extended by Leblond et al. [Leblond, J.B., Mottet, G., Devaux, J.C., 1986a. A theoretical and numerical approach to the plastic behavior of steels during phase transformations, I. Derivation of general relations. J. Mech. Phys. Solids 34, 395–409; Leblond, J.B., Mottet, G., Devaux, J.C., 1986b. A theoretical and numerical approach to the plastic behavior of steels during phase transformations, II. Study of classical plasticity for ideal-plastic phases. J. Mech. Phys. Solids 34, 411–432; Leblond, J.B., Devaux, J., Devaux, J.C., 1989a. Mathematical modeling of transformation plasticity in steels, I: case of ideal-plastic phases. Int. J. Plasticity 5, 511–572; Leblond, J.B., 1989b. Mathematical modeling of transformation plasticity in steels, II: coupling with strain hardening phenomena. Int. J. Plasticity 5, 573–591] were modified to consider the thermo-mechanical response of generalized multi-phase steel during phase transformations from austenite at high temperature. An implicit numerical solution procedure to calculate the plastic deformation of each constituent phase was newly proposed and implemented into the general purpose implicit finite element program via user material subroutine. The new algorithms include efficient calculation of consistent tangent modulus for the transformation plasticity and application of general anisotropic yield functions without limitation to the isotropic yield function. Besides the thermo-elastic–plastic constitutive equations, non-isothermal transformation kinetics was characterized by the Johnson–Mehl–Avrami–Kolmogorov (JMAK) equation and additivity relationship for the diffusional transformation, while the model proposed by Koistinen and Marburger was used for the diffusionless transformation. Numerical verifications for the continuous cooling experiments under various loading conditions were conducted to demonstrate the applicability of the developed numerical algorithms to the high carbon steel SK5.     

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.     

Stress-induced transformation behavior of a polycrystalline NiTi shape memory alloy: micro and macromechanical investigations via in situ optical microscopy

An experimental investigation of the micro and macromechanical transformation behavior of polycrystalline NiTi shape memory alloys was undertaken. Special attention was paid to macroscopic banding, variant microstructure, effects of cyclic loading, strain rate and temperature effects. Use of an interference filter on the microscope enabled observation of grain boundaries and martensitic plate formation and growth without recourse to etching or other chemical surface preparation. Key results of the experiments on the NiTi include observation of localized plastic deformation after only a few cycles, excellent temperature and stress relaxation correlation, a refined definition of “full transformation” for polycrystalline materials, and strain rate dependent effects. Several of these findings have critical implications for understanding and modeling of shape memory alloy behavior.     

应力引起的镍钛单晶形状记忆合金相变的实验研究

在单晶形状记忆合金试样中,由于没有晶粒之间的约束,它的马氏体相界面移动比多晶容易,用实验方法研究其相变的特点,对建立新的理论模型有意义,因而对它的实验分析显得重要。本文利用高分辨率的CCD系统监测到NiTi单晶形状记忆合金在拉伸时的相变伪弹性的过程;利用X射线衍射法得到了NiTi单晶试样在拉伸方向的晶向;运用高分辨率的云纹干涉技术,获得了应力引起的NiTi单晶形状记忆合金相变时的变形场;利用高分辨率、高灵敏度的红外相机记录了NiTi单晶在拉伸状态下的温度变化规律;对低温下NiTi单晶的拉伸性能做了初步的研究,得到一些有意义的现象。     

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

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

Meso-mechanical constitutive model for ratchetting of particle-reinforced metal matrix composites

Based on the Eshelby equivalent inclusion theory and Mori-Tanaka averaging method, a meso-mechanical cyclic elasto-plastic constitutive model is proposed to predict the ratchetting of particle-reinforced metal matrix composites. In the proposed model, a Hill-typed incremental formulation is used to simulate the elasto-plastic responses of the composites during cyclic loading with assumptions of elastic particle, elasto-plastic metal matrix and perfect interfacial bond between metal matrix and particles. A new nonlinear kinematic hardening rule extended from the Ohno-Abdel-Karim model [M. Abdel-Karim, N. Ohno, Kinematic hardening model suitable for ratchetting with steady-state, Int. J. Plasticity 16 (2000) 225-240] is employed to describe the ratchetting of metal matrix which dominates the ratchetting of the composites. With further assumption of spherical particles, the proposed meso-mechanical cyclic constitutive model is verified by comparing the predicted uniaxial ratchetting of SiC_P/6061Al composites with corresponding experiments obtained at room temperature [G.Z. Kang, Uniaxial time-dependent ratchetting of SiC_P/6061Al alloy composites at room and high temperature, Comp. Sci. Tech. 66 (2006) 1418-1430]. In the meantime, the effects of different tangent operators employed in the numerical implementation of the proposed model, i.e., continuum (or elasto-plastic) tangent operator C^ep and algorithmic (or consistent) one C^alg, on the predicted ratchetting are also discussed. It is concluded that the proposed model predicts the uniaxial ratchetting of SiC_P/6061Al composites at room temperature reasonably.     

NiTi 形状记忆合金热-力耦合循环变形行为宏微观实验和理论研究进展

本文对NiTi形状记忆合金热-力耦合循环变形行为研究的最新进展进行综述和评价。首先总结NiTi形状记忆合金在循环加载条件下的单轴、非比例多轴循环变形特性以及强烈的热-力耦合特性,阐述NiTi形状记忆合金在循环变形过程中出现功能性劣化的微观机理;然后,讨论在宏观和细观尺度上建立的三类NiTi形状记忆合金典型的循环本构模型,并评述代表性模型的预测能力;最后,总结已有研究存在的不足,对相关问题的进一步研究提出建议。在本构模型方面主要介绍了作者及其合作者在基于晶体塑性的热-力耦合循环本构模型方面的工作,突出了多种非弹性变形机制和强烈热-力耦合行为对形状记忆合金循环变形行为的影响。