On the Fremond’s constitutive model for shape memory alloys

The remarkable properties of shape memory alloys have increasing the interest in applications in different areas varying from biomedical to aerospace hardware. Despite the large number of applications, the modeling of SMA is the objective of many researches developed in order to describe all details of its thermomechanical behavior. The present contribution revisits a constitutive model presented by Savi et al. (2002), which is built up on the classical Fremond’s model, in order to contemplate the horizontal enlargement of the stress–strain hysteresis loop. Numerical simulations present qualitative agreement with experimental data, showing pseudoelastic, one-way and two-way shape memory effects.     

A thermodynamically-consistent microplane model for shape memory alloys

In microplane theory, it is assumed that a macroscopic stress tensor is projected to the microplane stresses. It is also assumed that 1D constitutive laws are defined for associated stress and strain components on all microplanes passing through a material point. The macroscopic strain tensor is obtained by strain integration on microplanes of all orientations at a point by using a homogenization process. Traditionally, microplane formulation has been based on the Volumetric–Deviatoric–Tangential split and macroscopic strain tensor was derived using the principle of complementary virtual work. It has been shown that this formulation could violate the second law of thermodynamics in some loading conditions. The present paper focuses on modeling of shape memory alloys using microplane formulation in a thermodynamically-consistent framework. To this end, a free energy potential is defined at the microplane level. Integrating this potential over all orientations provides the macroscopic free energy. Based on this free energy, a new formulation based on Volumetric–Deviatoric split is proposed. This formulation in a thermodynamic-consistent framework captures the behavior of shape memory alloys. Using experimental results for various loading conditions, the validity of the model has been verified.     

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.     

A three-dimensional phenomenological model for martensite reorientation in shape memory alloys

In this work, we propose a macroscopic phenomenological model that is based on the classical framework of thermodynamics of irreversible processes and accounts for the effect of multiaxial stress states and non-proportional loading histories. The model is able to account for the evolution of both twinned and detwinned martensite. Moreover, reorientation of the product phase according to loading direction is specifically accounted for. Towards this purpose the inelastic strain is split into two contributions deriving, respectively, from creation of detwinned martensite and reorientation of previously existing martensite variants. Computational tests demonstrate the ability of the model to simulate the main aspects of the shape memory response in a one-dimensional setting and some of the features that have been experimentally found in the case of multiaxial non-proportional loading histories. Experimental non-proportional loading paths have also been simulated and a good qualitative agreement between numerical and experimental response is observed.     

A 3-D phenomenological constitutive model for shape memory alloys under multiaxial loadings

This paper presents a new phenomenological constitutive model for shape memory alloys, developed within the framework of irreversible thermodynamics and based on a scalar and a tensorial internal variable. In particular, the model uses a measure of the amount of stress-induced martensite as scalar internal variable and the preferred direction of variants as independent tensorial internal variable. Using this approach, it is possible to account for variant reorientation and for the effects of multiaxial non-proportional loadings in a more accurate form than previously done. In particular, we propose a model that has the property of completely decoupling the pure reorientation mechanism from the pure transformation mechanism. Numerical tests show the ability to reproduce main features of shape memory alloys in proportional loadings and also to improve prediction capabilities under non-proportional loadings, as proven by the comparison with several experimental results available in the literature.     

Three-dimensional modeling and numerical analysis of rate-dependent irrecoverable deformation in shape memory alloys

Shape memory alloys (SMAs) provide an attractive solid-state actuation alternative to engineers in various fields due to their ability to exhibit recoverable deformations while under substantial loads. Many constitutive models describing this repeatable phenomenon have been proposed, where some models also capture the effects of rate-independent irrecoverable deformations (i.e., plasticity) in SMAs. In this work, we consider a topic not addressed to date: the generation and evolution of irrecoverable viscoplastic strains in an SMA material. Such strains appear in metals subjected to sufficiently high temperatures. The need to account for these effects in SMAs arises when considering one of two situations: the exposure of a conventional SMA material (e.g., NiTi) to high temperatures for a non-negligible amount of time, as occurs during shape-setting, or the utilization of new high temperature shape memory alloys (HTSMAs), where the elevated transformation temperatures induce transformation and viscoplastic behaviors simultaneously. A new three-dimensional constitutive model based on established SMA and viscoplastic modeling techniques is derived that accounts for these behaviors. The numerical implementation of the model is described in detail. Several finite element analysis (FEA) examples are provided, demonstrating the utility of the new model and its implementation in assessing the effects of viscoplastic behaviors in shape memory alloys.     

A finite strain kinematic hardening constitutive model based on Hencky strain: General framework, solution algorithm and application to shape memory alloys

The logarithmic or Hencky strain measure is a favored measure of strain due to its remarkable properties in large deformation problems. Compared with other strain measures, e.g., the commonly used Green-Lagrange measure, logarithmic strain is a more physical measure of strain. In this paper, we present a Hencky-based phenomenological finite strain kinematic hardening, non-associated constitutive model, developed within the framework of irreversible thermodynamics with internal variables. The derivation is based on the multiplicative decomposition of the deformation gradient into elastic and inelastic parts, and on the use of the isotropic property of the Helmholtz strain energy function. We also use the fact that the corotational rate of the Eulerian Hencky strain associated with the so-called logarithmic spin is equal to the strain rate tensor (symmetric part of the velocity gradient tensor). Satisfying the second law of thermodynamics in the Clausius-Duhem inequality form, we derive a thermodynamically-consistent constitutive model in a Lagrangian form. In comparison with the available finite strain models in which the unsymmetric Mandel stress appears in the equations, the proposed constitutive model includes only symmetric variables. Introducing a logarithmic mapping, we also present an appropriate form of the proposed constitutive equations in the time-discrete frame. We then apply the developed constitutive model to shape memory alloys and propose a well-defined, non-singular definition for model variables. In addition, we present a nucleation-completion condition in constructing the solution algorithm. We finally solve several boundary value problems to demonstrate the proposed model features as well as the numerical counterpart capabilities.     

Optimization of hysteretic characteristics of damping devices based on pseudoelastic shape memory alloys

In order to develop a fundamental understanding and the feasibility of SMA devices for passive vibration control, an undamped SDOF system with a pseudoelastic SMA restoring force is investigated to find the basic relationship between the shape of the hysteresis loop of SMA elements and their performance as a damping device. The dynamic characteristics of the device are evaluated by the steady-state response at the resonance point in order to focus on the damping effect. Dynamic analysis utilizing the equivalent linearization approach results in two major findings that, to the best of the authors’ knowledge, have not yet been reported in the literature. These results which characterize the unique behavior of the SMA hysteresis include: (a) for a given excitation amplitude, the “scale” of the hysteresis loop, which is a measure of displacement and restoring force, needs to be adjusted so that the response sweeps the maximum loop but does not exceed it; (b) the ratio of the area confined within the hysteresis loop to the area of a corresponding envelope of triangular shape should be as large as possible. The results of this study would be quite useful not only as a guideline for the design of actual SMA devices, but also as a basis for the development of new autoadaptive materials in future.     

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.     

A constitutive relation for pseudo-elastic behaviour in shape memory alloys

A constitutive relation to describe pseudo-elastic deformation in shape memory alloys is presented in this paper. It is capable of describing deformation behaviour of polycrystalline materials under triaxial stress state as well as of monocrystalline materials under one-dimensional condition. Total strain rate is supposed to be composed of elastic strain rate and transformation strain rate. Deformation behaviour of Cu−Zn−Sn alloy and Ti−ni alloy is simulated by use of the proposed constitutive relation. it is shown that simulated results are in a good agreement with experimental data. The project supported by National Natural Science Foundation of China.     

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.     

Three-dimensional phenomenological thermodynamic model of pseudoelasticity of shape memory alloys at finite strains

The familiar small strain thermodynamic 3D theory of isotropic pseudoelasticity proposed by Raniecki and Lexcellent is generalized to account for geometrical effects. The Mandel concept of mobile isoclinic, natural reference configurations is used in order to accomplish multiplicative decomposition of total deformation gradient into elastic and phase transformation (p.t.) parts, and resulting from it the additive decomposition of Eulerian strain rate tensor. The hypoelastic rate relations of elasticity involving elastic strain rate are derived consistent with hyperelastic relations resulting from free energy potential. It is shown that use of Jaumann corotational rate of stress tensor in rate constitutive equations formulation proves to be convenient. The formal equation for p.t. strain rate , describing p.t. deformation effects is proposed, based on experimental evidence. Phase transformation kinetics relations are presented in objective form. The field, coupled problem of thermomechanics is specified in rate weak form (rate principle of virtual work, and rate principle of heat transport). It is shown how information on the material behavior and motion inseparably enters the rate virtual work principle through the familiar bridging equation involving Eulerian rate of nominal stress tensor.

---

     

Thermomechanical coupling in shape memory alloys under cyclic loadings: Experimental analysis and constitutive modeling

In this paper, we examine the influence of thermomechanical coupling on the behavior of superelastic shape memory alloys subjected to cyclic loading at different loading rates. Special focus is given to the determination of the area of the stress-strain hysteresis loop once the material has achieved a stabilized state. It is found that this area does not evolve monotonically with the loading rate for either transient or asymptotic states. In order to reproduce this observation analytically, a new model is developed based on the ZM model for shape memory alloys which was modified to account for thermomechanical coupling. The model is shown to predict the non-monotonic variation in hysteresis area to good accord. Experimentally observed variations in the temperature of SMA test samples are also correctly reproduced for lower strain rates.     

Dislocation density based constitutive model for ultrasonic assisted deformation

The mechanical behaviour of metallic materials subjected to plastic deformation is altered with the superposition of ultrasonic vibrations. A significant effect is the reduction of flow stress or acoustic softening. This phenomenon is utilized in metal forming and other deformation based manufacturing processes. Experimental investigations on ultrasonic assisted tensile tests focus on the effect of ultrasonic vibrations along the longitudinal axis of the specimen, whereas the manufacturing processes employs in transverse directions. In the present work, transverse ultrasonic vibrations are imposed during a uniaxial tensile test using an aluminium alloy. The trend of acoustic softening due to transverse direction vibrations is similar to that along longitudinal direction. A dislocation density based constitutive model is extended to model the softening due to ultrasonic effect. The predicted results agree well with the experimental observations.     

Recoverable strains in composite shape memory alloys

New upper bounds are proposed for a generic problem of geometric compatibility, which covers the problem of bounding the effective recoverable strains in composite shape memory alloys (SMAs), such as polycrystalline SMAs or rigidly reinforced SMAs. Both the finite deformation and infinitesimal strain frameworks are considered. The methodology employed is a generalization of a homogenization approach introduced by Milton and Serkov [2000. Bounding the current in nonlinear conducting composites. J. Mech. Phys. Solids 48, 1295-1324] for nonlinear composites in infinitesimal strains. Some analytical and numerical examples are given to illustrate the method.     

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.     

Modeling the dynamic behavior of shape memory alloys

The paper studies the single degree of freedom vibration of a rigid mass suspended by a thin-walled shape memory alloy tube under torsional loading. The behavior is analyzed for the cases of quasiplasticity (low temperatures) and pseudoelasticity (high temperatures) on the basis of an improved version of the Müller–Achenbach model. To illustrate the strong hysteresis-induced damping capacity and the non-linear vibration characteristics, both, free and forced vibrations are considered in the first part of the paper. This is done on the basis of an isothermal version of the model, while the second part of the paper focuses on the effect of non-constant temperature caused by the rate-dependent release and absorption of latent heats.     

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.     

A micromechanics-inspired constitutive model for shape-memory alloys that accounts for initiation and saturation of phase transformation

A constitutive model to describe macroscopic elastic and transformation behaviors of polycrystalline shape-memory alloys is formulated using an internal variable thermodynamic framework. In a departure from prior phenomenological models, the proposed model treats initiation, growth kinetics, and saturation of transformation distinctly, consistent with physics revealed by recent multi-scale experiments and theoretical studies. Specifically, the proposed approach captures the macroscopic manifestations of three micromechanial facts, even though microstructures are not explicitly modeled: (1) Individual grains with favorable orientations and stresses for transformation are the first to nucleate martensite, and the local nucleation strain is relatively large. (2) Then, transformation interfaces propagate according to growth kinetics to traverse networks of grains, while previously formed martensite may reorient. (3) Ultimately, transformation saturates prior to 100% completion as some unfavorably-oriented grains do not transform; thus the total transformation strain of a polycrystal is modest relative to the initial, local nucleation strain. The proposed formulation also accounts for tension–compression asymmetry, processing anisotropy, and the distinction between stress-induced and temperature-induced transformations. Consequently, the model describes thermoelastic responses of shape-memory alloys subject to complex, multi-axial thermo-mechanical loadings. These abilities are demonstrated through detailed comparisons of simulations with experiments.