Open access peer-reviewed chapter

Modeling of the Two-Way Shape Memory Effect

By Meddour Belkacem and Brek Samir

Submitted: January 27th 2018Reviewed: February 19th 2018Published: September 26th 2018

DOI: 10.5772/intechopen.75657

Downloaded: 319


The shape memory alloys (SMA) are distinguished from other conventional materials by a singular behavior which takes many forms depending on the thermomechanical load. The two-way shape memory effect is one of these forms. The interest that exhibits this behavior is that the material can remember two states, so this leads to many industrial applications. This thermoelastic property is driven by the temperature under residual stress of education. To model this effect in 3D, we considered stress and temperature as control variables and the fraction of martensite as internal variable; choosing Gibbs free energy expression and applying thermodynamic principles with transformation criteria have permitted to write the constitutive equations that control this behavior. The constructed model is then numerically simulated, and finally, the proposed model appears applicable in engineering.


  • two-way
  • simulation
  • hysteresis
  • transformation

1. Introduction

Shape memory alloys (SMA) as they are called are a kind of particular materials which have a singular behavior, they can be largely deformed (about 10%) under an applied mechanical, a simple heating is sufficient to recover the previous form, that is why they are called that way.

By varying controlled parameters (stress and temperature), these alloys can exhibit other properties like pseudoelasticity, two-way shape memory effect, one-way shape memory effect [1, 2], and reorientation effect [3].

These properties derive from phase transformations, i.e., higher temperature phase (austenite) to lower temperature phase (martensite).

It is observed that these phase transformations do not occur with diffusion but rather with displacement, i.e., displacement at a distance less than interatomic [4, 5].

These important properties made them requested materials in various fields such as biomedical automotive industry, fire watch devices, aeronautics, and medical devices [6].

Two-way shape memory effect is obtained after the alloy is subjected to an education [7], i.e., to a cyclic thermal load under a constant mechanical load (Figure 1).

Figure 1.

Education (50 cycles) performed on a NiTi wire under constant stress of 150 MPa [7].

It should be noted that the previous education is performed to create a field of internal stresses, which will substitute the macroscopic ones; under a cyclic thermal load, the alloy will remember two states: one is at lower temperature (martensite) and the other at higher temperature (austenite) (Figure 2).

Figure 2.

One-dimensional two-way effect. Msσ: Temperature of transformation start A- > M under stress σ.Mfσ: Temperature of transformation finish A- > M under stress σ.Asσ: Temperature of transformation start M- > A under stress σ. Afσ: Temperature of transformation finish M- > A under stress σ. ε0: Uniaxial maximum deformation.

The next steps in this paper are to build the constitutive model and simulate it using an algorithm; regarding model parameters, we will use the work of [8].

2. Methods and materials

2.1. Constitutive equations

Let us choose the following free energy expression:


SA: Fourth order tensor of complaisance; B, C: Constants to be determined, respectively, related to change of phase and interaction between austenite and martensite; f: Fraction of martensite; σ00: Tensor of stress created after education.

Assuming that dissipation is associated only with transformation (fraction of martensite) [9, 10], the second principle of thermodynamics can be written as


Let us write the driving force Fth:




Because of hysteresis, there is a dissipative force Fdi, which will oppose Fth.

We choose expression of Fdi:


Transformation occurs when


We introduce the criteria functions:


Condition of consistence gives


Doing the same with the reverse transformation


Eqs. (11) and (13) give the evolution of fraction of martensite.

The deformation resulting from transformation of austenite to martensite is denoted εT; this deformation is associated with fraction of martensite:


R is a tensor of transformation, and it can be written as the following:


2.2. Determination of constants B, C, K, and H

At the beginning of direct transformation


At the end of the direct transformation


At the beginning of the reverse transformation


At the end of the reverse transformation


2.3. Experimental data

Table 1 illustrates experimental data and material constants B, C, K, and H. The used material is CuZnAl. The test was performed under constant stress (σ=65MPa).


Table 1.

Experimental data [8, 11].

2.4. Numerical simulation

2.4.1. One-dimensional case

We considered a segment of CuZnAl submitted to a constant stress (σ=65MPa) and a thermal load (300T400K) (Figure 3).

Figure 3.

History of coupled loading.

2.4.2. Three-dimensional case

The considered specimen is a cubic element subjected to a constant mechanical load which is a triaxial traction and thermal loads (σ11=σ22=σ33=50MPa) and (300T400K) (Figures 46).

Figure 4.

Loading in the direction of σ11.

Figure 5.

Loading in the direction of σ22.

Figure 6.

Loading in the direction of σ33.

3. Results

3.1. One-dimensional case

Figure 7.

Response at (σ=65MPa).

Figure 8.

Evolution of the fraction of martensite.

3.2. Three-dimensional case

Figure 9.

Plot T−ε11.

Figure 10.

Plot T−ε22.

Figure 11.

Plot T−ε22.

Figure 12.

Evolution of fraction of martensite in direct transformation.

Figure 13.

Evolution of fraction of martensite in reverse transformation.

4. Discussion

After having written constitutive equations and criteria functions, the numerical simulation has permitted to obtain previous results.

First, we used the extracted values from the curve of the test in order to compare this curve with the response of the model in one-dimensional case, and we obtained Figure 7. This figure presents a good agreement between experimental data and the model response.

On the other hand, Figure 8, which is representing the evolution of the martensite fraction, is also in agreement with the curve in Figure 7, i.e., the direct transformation and the reverse transformation are functions of martensite fraction. We can say that the constitutive model behaves well in one-dimensional case.

For the three-dimensional case, Figures 911 show the response under the triaxial traction and thermal load; for each plot, there is a hysteresis.

The shrinking of the hysteresis in each case of the figures should be noted (Figures 911); this is due to the triaxial loading.

Despite the applied triaxial load, Figures 911 exhibit the thermomechanical cycle as it is in case of one-dimensional two-way effect (Figure 7).

Figures 12 and 13 show the evolution of fraction of martensite for each case of direct and reverse transformations, and the shapes of the plots are compatible with Figures 911 as the one-dimensional case because it was noticed previously that the deformation is proportional to the amount of martensite.

5. Conclusion

In this work, we have developed a 3D constitutive model using the principles of thermodynamics and a simple formalism, and these principles have permitted to write criteria of transformation. This macroscopic model is developed by simple formalism and assumptions.

By using an algorithm, we have implemented the model, and the response seems to be compatible with the nature of the two-way shape memory effect. In the one-dimensional case, we have observed a good agreement between the numerical and experimental plots.

It should be noted that the parameters of the model were determined by the one-dimensional test and further used in the biaxial and triaxial cases to ensure consistency of the model in different cases of loading. The implementation of the model in the algorithm is simple and practical. The obtained results testify the usability of the developed model.

At the end, we can say that this macroscopic constitutive model can be used in applications to engineering problems, in order to particularly simulate the pseudoelastic effect of shape memory alloys.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Meddour Belkacem and Brek Samir (September 26th 2018). Modeling of the Two-Way Shape Memory Effect, Shape-Memory Materials, Alicia Esther Ares, IntechOpen, DOI: 10.5772/intechopen.75657. Available from:

chapter statistics

319total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Linear Shape Memory Alloy Thermomechanical Actuators

By Velaphi Msomi and Graeme Oliver

Related Book

First chapter

Inelastic X-Ray Scattering as a Probe of the Transition Between the Hydrodynamic and the Single Particle Regimes in Simple Fluids

By Alessandro Cunsolo

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us