Smart Control Algorithms and Technology in Civil InfrastructuresView this Special Issue
A Constitutive Model for Superelastic Shape Memory Alloys Considering the Influence of Strain Rate
Shape memory alloys (SMAs) are a relatively new class of functional materials, exhibiting special thermomechanical behaviors, such as shape memory effect and superelasticity, which enable their applications in seismic engineering as energy dissipation devices. This paper investigates the properties of superelastic NiTi shape memory alloys, emphasizing the influence of strain rate on superelastic behavior under various strain amplitudes by cyclic tensile tests. A novel constitutive equation based on Graesser and Cozzarelli’s model is proposed to describe the strain-rate-dependent hysteretic behavior of superelastic SMAs at different strain levels. A stress variable including the influence of strain rate is introduced into Graesser and Cozzarelli’s model. To verify the effectiveness of the proposed constitutive equation, experiments on superelastic NiTi wires with different strain rates and strain levels are conducted. Numerical simulation results based on the proposed constitutive equation and experimental results are in good agreement. The findings in this paper will assist the future design of superelastic SMA-based energy dissipation devices for seismic protection of structures.
Shape memory alloys (SMAs) are a unique class of materials that have the ability to undergo large deformations, up to 8~10%, that is, at least one order of magnitude greater than common metals and alloys, and revert back to their original and undeformed shape or dimension through either applications of heat, that is, the shape memory effect (SME), or removal of stress, that is, the superelastic effect.
The particular properties of SMAs were first discovered by Chang and Read in 1951; however, it was not until after 1962 when Buechler and his colleagues found the shape memory effect in nickel-titanium (NiTi) at the Naval Ordnance Laboratory that both in-depth research and practical applications emerged. At present, SMAs have been wildly implemented in biomedical, aerospace, mechanical, and civil engineering areas [1–5], making it necessary to have a precise understanding of the special mechanical behavior of SMAs in order to fully develop and exploit their potential.
The unique mechanical behaviors of SMAs are made possible by reversible martensitic phase transformation (MPT) induced by temperature or mechanical stress between the austenitic phase () and martensitic phase (). Austenite has a body-centered cubic crystal structure (crystallographic high symmetry), stable at high temperatures and low stress values, while the martensite has a parallelogram structure (low symmetry), stable at low temperatures and high stress values, having up to 24 variations. In the stress-free state, the material is characterized by four transition temperatures, namely, martensite start temperature , martensite finish temperature , austenite start temperature , and austenite finish temperature . Below the martensite finish temperature, , residual deformations induced by the martensite reorientation due to applied stresses can be recovered by heating the material above the austenite finish temperature, , resulting in the shape memory effect (SME). Above the austenite finish temperature, , martensite is formed associated with forward MPT (from austenite to detwinned martensite) induced by stress when load exceeds the forward transformation stress; however, the material reverts to austenite at a lower stress without residual deformation associated with inverse MPT (from detwinned martensite to austenite) after unloading, resulting in the superelastic behavior.
Over the past decades, there has been a significant amount of research dedicated to the martensitic phase transformation [6, 7] and the factors that induce the MPT [8–10]. Meanwhile, some constitutive models have been proposed. A short review of the existing models is presented here:(1)phenomenological macroscopic constitutive models in terms of stress, strain, and temperature with assumed phase transformation kinetics described by preestablished simple mathematical functions proposed by Tanaka , Liang and Rogers , Brinson , Boyd and Lagoudas [14, 15], Li et al. , Tobushi et al. , and Sun and Rajapakse ; among others. (2)one-dimensional polynomial models based on Devonshire’s theory with an assumed polynomial-free energy potential, which allows superelasticity and SME description, presented by Falk et al. [19, 20];(3)thermodynamic models based on the free energy and dissipation potential developed by Patoor et al. , Sun and Hwang [22, 23], Huang and Brinson , and Boyd and Lagoudas ;(4)plastic flow models based on dislocation theories of solid state physics proposed by Graesser and Cozzarelli [26, 27], lately improved by Wilde et al. , Zhang and Zhu , and Ren et al. .
For application of SMAs in earthquake engineering, Graesser and Cozzarelli’s model  will be adopted in this paper due to its advantage of simplicity. Graesser and Cozzarelli’s model is capable of capturing both the superelastic effect and the martensitic hysteresis with their subloops due to partial transformation processes. Wilde et al.  improve the Graesser and Cozzarelli model by introducing the possibility of describing the material behavior also after phase transformation completion as well as by simulating a smooth transition between the elastic branch and the superelastic plateau. Zhang and Zhu  modified Wilde’s model to enhance the stability of numerical simulation and speed up the computation time. However, Graesser and Cozzarelli’s model and the improved editions cannot describe the strain-rate-dependent property of SMA. To consider the effect of the load path to increase its modeling accuracy, Ren et al.  improved Graesser and Cozzarelli’s model by dividing the full loop into three parts: the loading branch, unloading branch before the completion of the reverse transformation and the elastic unloading branch after the completion of the reverse transformation, where each part adopts its respective parameters. However, the improved model can only reflect one case of hysteresis for a fixed loading rate, since when the loading rate changes, the model parameters will change. Therefore, it is not convenient to perform continuous simulation of SMA devices under seismic excitations, which involve many different loading rates.
As a further research to improve Graesser and Cozzarelli’s model, this paper aims to develop a novel strain-rate-dependent constitutive model, which can simultaneously account for the effects of both strain rates and strain levels. In Section 2 of this paper, Graesser and Cozzarelli’s model is improved by adding a stress variable including the influence of strain rate. To verify the effectiveness of the proposed constitutive equation, cyclic tensile tests on superelastic NiTi wires with different strain rates and strain levels are carried out, as described in Section 3, and the experimental results and analysis are given in Section 4. In Section 5, numerical simulations are conducted to demonstrate the accuracy of the improved Graesser and Cozzarelli’s model.
2. Constitutive Model of Shape Memory Alloys
In this section, a uniaxial constitutive model, based on Graesser and Cozzarelli’s model, is developed to capture the strain-rate-dependent property of superelastic SMA materials.
2.1. Graesser and Cozzarelli’s Model
Based on Ozdemir’s model , Graesser and Cozzarelli proposed a one-dimensional hysteretic model that produced the general macroscopic stress-strain characteristics of SMAs . The equation is given as where is the one-dimensional stress, is the strain, is the elastic modulus, is the transformation stress, is a constant assumed to be any positive odd real value controlling the sharpness of the transition from the elastic state to the phase transformation, and , respectively, denote the ordinary time derivative of the stress and strain, and is the one-dimensional back stress, given by where , , and are material constants controlling the type and size of the hysteresis, the amount of elastic recovery during unloading, and the slope of the unloading stress plateau, respectively. is a constant controlling the slope of the curve in the inelastic range, given by where is the slope of the curve in the inelastic range.
is the inelastic strain, given by is the unit step function defined as is the error function defined by Graesser and Cozzarelli’s model has a relatively simple expression with parameters that can be easily acquired and it is easy to implement, yet this model can still be improved. Though the equations are written in differential form, the model is essentially strain rate independent. Studies [8–10] show that the properties are strongly dependent on the strain rate. In this paper, Graesser and Cozzarelli’s model will be improved to consider the influence of strain rate.
2.2. Improved Graesser and Cozzarelli’s Model
To simulate the strain-rate-dependent hysterestic behavior of superelastic SMAs, Graesser and Cozzarelli’s model will be improved. The following assumptions are made for this purpose.(i)The ambient temperature is constant.(ii)The effect of strain rate on the properties of SMAs can be ignored when it is below /s in this paper).
Based on the above assumptions, the stress is formulated here. The stress under dynamic loading can be separated into two parts; namely, where , , and are the stress under dynamic loading, the stress under quasistatic loading, and the stress change due to the influence of the strain rate, respectively. To keep consistency with and in symbol style, is used to represent the stress change. Actually, the value of is equal to that of which will be illustrated in Section 4.
The definition of depends on the relationship between and . The equations capturing the relationship between and , as well as can be presented from the experimental tests (presented in Sections 3 and 4) based on mathematical statistics and numerical fitting method. is given as follows.(1) For , (2) For , where is a material constant. If , then Differentiating (12) results in (3) For , Thus, its time derivative can be written as where is a material parameter related to the maximum strain amplitude, . Other unnoted parameters have already been defined.
3. Experimental Setup for Tensile Testing of Superelastic SMA Wires
In the following, cyclic tensile tests on superelastic NiTi wires with different strain rates and strain levels are carried out to verify the effectiveness of the proposed constitutive equation.
3.1. Test Specimens and Equipment
The NiTi SMA wires used for testing have a diameter of 0.5 mm. The TiNi SMA is an alloy with approximate 50.9% Ni and 49.1% Ti. Under zero external stress, the martensite start and finish temperatures and the austenite start and finish temperatures (, , , and ), measured by DSC (differential scanning calorimeter), are −73°C, −55°C, −23°C, and 5°C, respectively.
Tests were conducted using an electromechanical universal testing machine. The SMA wire specimens, with a 100 mm test length between the two custom-made grips, were subjected to triangular cyclic loading under different strain amplitudes. In order to acquire reliable hysteresis behavior of SMA, three specimens were applied in each testing. The mean values were utilized for analysis. The strains were calculated from the elongation measured by a 50 mm gage length extensometer with the stress calculated from the axial force, which was measured by a 10 KN load cell. In each test, the specimen was required to follow a triangular wave with a constant strain rate. The data was recorded automatically by a PC-based data acquisition system with a 30 Hz sampling rate. All the tests were carried out at room temperature (≈25°C). The experimental setup and strain loading curve are shown in Figures 1(a) and 1(b), respectively.
(a) Experimental setup
(b) Strain loading curve
3.2. Test Procedure
The experimental results in references [8, 9] have shown the mechanical behaviors of the SMA trend stable with the increasing of the cyclic number. Therefore, prior to testing, the NiTi SMA specimens were cycled 30 times at 6.5% strain amplitude and 1.0 × 10−3/s strain rate as a “training” process to reach a steady-state condition. The scheme of the tests is described as follows.(1)At quasistatic loading (1.0 × 10−4/s strain rate), the NiTi SMA specimens are subjected to cyclic loading with different strain levels, ranging from 2% to 6% at an increment of 2%. (2)At strain rates of 5.0 × 10−4/s, 1.0 × 10−3/s, 2.5 × 10−3/s, and 5.0 × 10−3/s, the NiTi SMA specimens are subjected to cyclic loading with different strain levels, ranging from 2% to 6% at an increment of 2%.
What needs to be pointed out is that the maximum strain rate is only 5.0 × 10−3/s due to the limitation of the experimental condition. The range is relatively low in seismic engineering. More tests with higher strain rates will be conducted in the future.
4. Results and Analysis
Many experimental studies have shown that the superelastic behavior of SMAs strongly depends on the loading rate [8–10]; however, few quantitative relationships are provided [16–18]. In this section, a stress change variable, , which is a relative quantity between the dynamic loading condition and the static loading condition at an identical strain, is introduced to evaluate the effect of the strain rate on the superelastic behavior of SMAs.
4.1. Variation of Stress during Forward Martensitic Transformation at Different Strain Rates
Figure 2(a) shows the relationship between the amount of stress change, , and the strain obtained by the tests during the forward martensitic transformation at different strain rates. The points with different shapes represent the experimental results. As we can see, increases with an increase of strain. The rate of increase of is relatively large at small strain levels; however it decreases with an increasing strain. Moreover, the stress change increases with an increasing strain rate.
On one hand, for a fixed strain rate, the relationship between stress change and strain level can be described using a power function with the index less than 1. In order to simplify the equation, the index selected is 1/2. On the other hand, for a fixed strain level, the relationship between stress change and strain rate can be captured by natural logarithms function. Therefore, based on numerical fitting tools, the stress change , as a function of the strain level and the strain rate during the forward martensitic transformation, can be expressed by the following equation: where is the strain and is the strain rate. The fitted curves calculated from (16) for different strain rates are shown in Figure 2(a). As seen from (16), the material constant in (12) is equal to 159.6.
4.2. Variation of Stress during Inverse Martensitic Transformation at Different Strain Rates
Figure 2(b) shows the relationship between the amount of stress change and the strain during the inverse martensitic transformation at different strain rates with the maximum strain amplitude at 6%. Likewise, the points with different shapes represent the experimental results. As opposed to the forward martensitic transformation, the rate of increase in is relatively small at small strain levels; however it increases with an increasing strain. Moreover, the stress change increases with an increasing strain rate. This is because the stress-strain loops of superelastic SMA translate upwards, while the branches of the curve relevant to the phase transformations harden with the increasing of strain rate, as noted in references [8, 10].
Similarly, for a given strain rate, the relationship between stress change and strain level can be described using a power function with the index more than 1. The index selected is 2. The relationship between stress change and strain rate is similar to that during the forward martensitic transformation. In the case of the maximum strain amplitude 6%, the stress change, , as a function of the strain level and the strain rate during the inverse martensitic transformation, can be expressed by the following equation: The fitted curves calculated from (17) for different strain rates are also shown in Figure 2(b). As seen from (17), the material constant in (14) is equal to 500 for the strain level at 6%.
5. Numerical Simulations
In this section, the effectiveness of the improved model to reproduce the rate-dependent superelastic behavior of SMAs wire is validated by comparing numerical results with experimental data. In order to verify suitability of the model, two different conditions, static loading and dynamic loading, are investigated.
5.1. Comparison of Model Predictions and Experimental Data under Quasistatic Loading
Figure 3 shows the comparison of the model predictions and the experimental responses of the SMA wire at 2%, 4%, and 6% strain levels under quasistatic loading with a strain rate of 1.0 × 10−4/s. The key parameters used in the models are MPa, N, , , , , and . As can be seen from Figure 3, the numerical curves predicted by the improved model agree well with the experimental data.
To make an additional quantitative observation, comparisons of the energy dissipation per cycle, tensile stress at peak strain, and equivalent viscous damping between the experimental data and numerical results are shown in Table 1. Note that the maximum differences for energy dissipation per cycle, tensile stress at peak strain, and equivalent viscous damping are 2.9%, 1.8%, and 2.6%, respectively. Therefore, the results indicate that the mechanical behavior of superelastic SMA wire under quasistatic loading was predicted well by the improved model.
5.2. Comparison of Model Predictions and Experimental Data under Dynamic Loading
Figure 4 shows the comparison of model predictions and experimental responses of the SMA wire at 4% strain under the strain rates of 5.0 × 10−4/s, 1.0 × 10−3/s, 2.5 × 10−3/s, and 5.0 × 10−3/s, respectively. The key parameters used in the models are MPa, N, , , , , and , , , and . It is obvious from Figures 4(a)–4(d) that the superelastic behavior at different strain rates is well predicted by the improved constitutive model.
Table 2 shows the comparisons of the energy dissipation per cycle, the tensile stress at peak strain, and the equivalent viscous damping between the experimental data and numerical results. It is clear that the maximum difference of energy dissipation per cycle, the tensile stress at peak strain, and the equivalent viscous damping are 6.89%, 4.01%, and 9.36%, respectively. From both Figure 4 and Table 2, it can be concluded that the improved Graesser and Cozzarelli’s model is able to capture the superelastic behavior in dynamic loading conditions within the range of rates tested in this study.
In this paper, the influence of the strain rate on the mechanical properties of SMAs is investigated. A one-dimensional strain-rate-dependent constitutive model based on Graesser and Cozzarelli’s model is proposed to predict the hysteretic behavior of superelastic SMAs. In this model, the stress is divided into two parts: the static stress and dynamic stress change. The former is based on the original Graesser and Cozzarelli’s model and describes the property under quasistatic loading. The latter one considers the effect of the strain rate. Comparisons of model predictions and experimental results at different strain levels and strain rates are performed and reveal that the improved Graesser and Cozzarelli’s model can accurately predict hysteretic behavior of superelastic SMAs under both static and dynamic loading conditions within the range of rates tested in this study. In this paper, the maximum strain rate is only 5.0 × 10−3/s due to the limitation of the experimental condition. The range is relatively low in seismic engineering. Future research will be conducted at higher strain rates to further validate the effectiveness of the model.
The research reported in this paper was supported by the National Natural Science Foundation of China (no. 51108426), China Postdoctoral Science Foundation (no. 20100471008), and Research Found for the Doctoral Program of Higher Education of China (no. 20104101120009). These supports are greatly appreciated.
G. Song, N. Ma, and H.-N. Li, “Review of applications of shape memory alloys in civil structures,” in Proceedings of the Earth and Space Conference (ASCE '04), pp. 559–566, March 2004.View at: Google Scholar
K. Tanaka, “A thermomechanical sketch of shape memory effect: one-dimensional tensile behavior,” Res Mechanica, vol. 18, no. 3, pp. 251–263, 1986.View at: Google Scholar
C. Liang and C. A. Rogers, “One-dimensional thermomechanical constitutive relations for shape memory materials,” Journal of Intelligent Material Systems and Structures, vol. 1, no. 2, pp. 207–234, 1990.View at: Google Scholar
L. C. Brinson, “One-dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable,” Journal of Intelligent Material Systems and Structures, vol. 4, no. 2, pp. 229–242, 1993.View at: Publisher Site | Google Scholar
H. Tobushi, I. Shimeno, T. Hachisuka, and K. Tanaka, “Influence of strain rate on superelastic behaviour of TiNi shape memory alloy,” Proceedings of the Institution of Mechanical Engineers Part L: Journal of Materials: Design and Applications, vol. 213, no. 2, pp. 93–102, 1999.View at: Google Scholar
E. Patoor, A. Eberhardt, and M. Berveiller, “Thermomechanical behavior of shape memory alloys,” Archives of Mechanics, vol. 40, pp. 775–794, 1988.View at: Google Scholar
Q. P. Sun and K. C. Hwang, “Micromechanics modelling for the constitutive behavior of polycrystalline shape memory alloys—II. Study of the individual phenomena,” Journal of the Mechanics and Physics of Solids, vol. 41, no. 1, pp. 19–33, 1993.View at: Publisher Site | Google Scholar | Zentralblatt MATH
J. G. Boyd and D. C. Lagoudas, “A thermodynamical constitutive model for the shape memory effect due to transformation and reorientation,” in Smart Structures and Materials, vol. 2189 of Proceedings of SPIE, pp. 276–288, 1994.View at: Google Scholar