Smart Control Algorithms and Technology in Civil Infrastructures
View this Special IssueResearch Article  Open Access
Hui Qian, Hongnan Li, Gangbing Song, Wei Guo, "A Constitutive Model for Superelastic Shape Memory Alloys Considering the Influence of Strain Rate", Mathematical Problems in Engineering, vol. 2013, Article ID 248671, 8 pages, 2013. https://doi.org/10.1155/2013/248671
A Constitutive Model for Superelastic Shape Memory Alloys Considering the Influence of Strain Rate
Abstract
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 strainratedependent 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 SMAbased energy dissipation devices for seismic protection of structures.
1. Introduction
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 nickeltitanium (NiTi) at the Naval Ordnance Laboratory that both indepth 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 bodycentered 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 stressfree 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 [11], Liang and Rogers [12], Brinson [13], Boyd and Lagoudas [14, 15], Li et al. [16], Tobushi et al. [17], and Sun and Rajapakse [18]; among others. (2)onedimensional polynomial models based on Devonshire’s theory with an assumed polynomialfree 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. [21], Sun and Hwang [22, 23], Huang and Brinson [24], and Boyd and Lagoudas [25];(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. [28], Zhang and Zhu [29], and Ren et al. [30].
For application of SMAs in earthquake engineering, Graesser and Cozzarelli’s model [26] 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. [28] 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 [29] 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 strainratedependent property of SMA. To consider the effect of the load path to increase its modeling accuracy, Ren et al. [30] 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 strainratedependent 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 strainratedependent property of superelastic SMA materials.
2.1. Graesser and Cozzarelli’s Model
Based on Ozdemir’s model [30], Graesser and Cozzarelli proposed a onedimensional hysteretic model that produced the general macroscopic stressstrain characteristics of SMAs [26]. The equation is given as where is the onedimensional 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 onedimensional 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 strainratedependent 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.
Differentiating (7) results in where The parameters in (9) have already been defined in Section 2.1 and can be obtained experimentally by the stressstain curve under quasistatic loading.
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 custommade 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 PCbased 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 steadystate 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.
(a)
(b)
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 stressstrain 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 ratedependent 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.
(a)
(b)
(c)
(d)
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.

6. Conclusions
In this paper, the influence of the strain rate on the mechanical properties of SMAs is investigated. A onedimensional strainratedependent 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.
Acknowledgments
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.
References
 S. Saadat, J. Salichs, M. Noori et al., “An overview of vibration and seismic applications of NiTi shape memory alloy,” Smart Materials and Structures, vol. 11, no. 2, pp. 218–229, 2002. View at: Publisher Site  Google Scholar
 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
 L. Janke, C. Czaderski, M. Motavalli, and J. Ruth, “Applications of shape memory alloys in civil engineering structures—overview, limits and new ideas,” Materials and Structures, vol. 38, no. 279, pp. 578–592, 2005. View at: Publisher Site  Google Scholar
 H. N. Li and L. S. Huo, “Advances in structural control in civil engineering in China,” Mathematical Problems in Engineering, vol. 2010, Article ID 936081, 23 pages, 2010. View at: Publisher Site  Google Scholar
 R. Johnson, J. E. Padgett, M. E. Maragakis, R. Desroches, and M. S. Saiidi, “Large scale testing of nitinol shape memory alloy devices for retrofitting of bridges,” Smart Materials and Structures, vol. 17, no. 3, Article ID 035018, 13 pages, 2008. View at: Publisher Site  Google Scholar
 M. S. Alam, M. A. Youssef, and M. Nehdi, “Utilizing shape memory alloys to enhance the performance and safety of civil infrastructure: A review,” Canadian Journal of Civil Engineering, vol. 34, no. 9, pp. 1075–1086, 2007. View at: Publisher Site  Google Scholar
 K. L. Ng and Q. P. Sun, “Stressinduced phase transformation and detwinning in NiTi polycrystalline shape memory alloy tubes,” Mechanics of Materials, vol. 38, no. 12, pp. 41–56, 2006. View at: Publisher Site  Google Scholar
 M. Dolce and D. Cardone, “Mechanical behaviour of shape memory alloys for seismic applications 2. Austenite NiTi wires subjected to tension,” International Journal of Mechanical Sciences, vol. 43, no. 11, pp. 2657–2677, 2001. View at: Publisher Site  Google Scholar
 R. Desroches, J. McCormick, and M. Delemont, “Cyclic properties of superelastic shape memory alloy wires and bars,” Journal of Structural Engineering, vol. 130, no. 1, pp. 38–46, 2004. View at: Publisher Site  Google Scholar
 D. Wolons, F. Gandhi, and B. Malovrh, “Experimental investigation of the pseudoelastic hysteresis damping characteristics of shape memory alloy wires,” Journal of Intelligent Material Systems and Structures, vol. 9, no. 2, pp. 116–126, 1998. View at: Publisher Site  Google Scholar
 K. Tanaka, “A thermomechanical sketch of shape memory effect: onedimensional tensile behavior,” Res Mechanica, vol. 18, no. 3, pp. 251–263, 1986. View at: Google Scholar
 C. Liang and C. A. Rogers, “Onedimensional 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, “Onedimensional constitutive behavior of shape memory alloys: thermomechanical derivation with nonconstant 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
 J. G. Boyd and D. C. Lagoudas, “A thermodynamical constitutive model for shape memory materials. Part I. The monolithic shape memory alloy,” International Journal of Plasticity, vol. 12, no. 6, pp. 805–842, 1996. View at: Publisher Site  Google Scholar
 J. G. Boyd and D. C. Lagoudas, “A thermodynamical constitutive model for shape memory materials. Part II. The SMA composite material,” International Journal of Plasticity, vol. 12, no. 7, pp. 843–873, 1996. View at: Publisher Site  Google Scholar
 H. Li, C.X. Mao, and J.P. Ou, “Strain selfsensing property and strain rate dependent constitutive model of austenitic shape memory alloy: experiment and theory,” Journal of Materials in Civil Engineering, vol. 17, no. 6, pp. 676–685, 2005. 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
 S. Sun and R. K. N. D. Rajapakse, “Simulation of pseudoelastic behaviour of SMA under cyclic loading,” Computational Materials Science, vol. 28, no. 34, pp. 663–674, 2003. View at: Publisher Site  Google Scholar
 F. Falk, “Model free energy, mechanics, and thermodynamics of shape memory alloys,” Acta Metallurgica, vol. 28, no. 12, pp. 1773–1780, 1980. View at: Publisher Site  Google Scholar
 F. Falk and P. Konopka, “Threedimensional Landau theory describing the martensitic phase transformation of shapememory alloys,” Journal of Physics, vol. 2, no. 1, pp. 61–77, 1990. View at: Publisher Site  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—I. Derivation of general relations,” Journal of the Mechanics and Physics of Solids, vol. 41, no. 1, pp. 1–17, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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
 M. S. Huang and L. C. Brinson, “A multivariant model for single crystal shape memory alloy behavior,” Journal of the Mechanics and Physics of Solids, vol. 46, no. 8, pp. 1379–1409, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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
 E. J. Graesser and F. A. Cozzarelli, “Shapememory alloys as new materials for aseismic isolation,” Journal of Engineering Mechanics, vol. 117, no. 11, pp. 2590–2608, 1991. View at: Publisher Site  Google Scholar
 E. J. Graesser and F. A. Cozzarelli, “Proposed threedimensional constitutive model for shape memory alloys,” Journal of Intelligent Material Systems and Structures, vol. 5, no. 1, pp. 78–89, 1994. View at: Publisher Site  Google Scholar
 K. Wilde, P. Gardoni, and Y. Fujino, “Base isolation system with shape memory alloy device for elevated highway bridges,” Engineering Structures, vol. 22, no. 3, pp. 222–229, 2000. View at: Publisher Site  Google Scholar
 Y. Zhang and S. Zhu, “A shape memory alloybased reusable hysteretic damper for seismic hazard mitigation,” Smart Materials and Structures, vol. 16, no. 5, pp. 1603–1613, 2007. View at: Publisher Site  Google Scholar
 W. J. Ren, H. N. Li, and G. Song, “A onedimensional strainratedependent constitutive model for superelastic shape memory alloys,” Smart Materials and Structures, vol. 16, no. 1, pp. 191–197, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Hui Qian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.