Acessibilidade / Reportar erro

Analysis of the transformations temperatures of helicoidal Ti-Ni actuators using computational numerical methods

Abstract

The development of shape memory actuators has enabled noteworthy applications in the mechanical engineering, robotics, aerospace, and oil industries and in medicine. These applications have been targeted on miniaturization and taking full advantage of spaces. This article analyses a Ti-Ni shape memory actuator used as part of a flow control system. A Ti-Ni spring actuator is subjected to thermomechanical training and parameters such as transformation temperature, thermal hysteresis and shape memory effect performance were investigated. These parameters were important for understanding the behavior of the actuator related to martensitic phase transformation during the heating and cooling cycles which it undergoes when in service. The multiple regression methodology was used as a computational tool for analysing data in order to simulate and predict the results for stress and cycles where the experimental data was not developed. The results obtained using the training cycles enable actuators to be characterized and the numerical simulation to be validated.

shape memory effect; Ti-Ni actuators; multiple regression analysis


Analysis of the transformations temperatures of helicoidal Ti-Ni actuators using computational numerical methods

Carlos Augusto do N. OliveiraI, ** e-mail: cano.oliveira@gmail.com; Alvaro Antonio Ochoa VillaI; Cezar Henrique GonzalezI; Pablo Batista GuimarãesII; Rodrigo José FerreiraII; Severino Leopoldino Urtiga FilhoI; Niedson José da SilvaII; José Orlando Silveira RochaI

IDepartment of Mechanical Engineering, Centre for Technology and Geosciences, Federal University of Pernambuco - UFPE, Av. Acadêmico Hélio Ramos, s/n, Cidade Universitária, CEP 50740-530, Recife, PE, Brazil

IIDepartment of Mechanics, Federal Institute of Pernambuco - IFPE, Av. Luiz Freire 500, Cidade Universitária, CEP 50740-540, Recife, PE, Brazil

ABSTRACT

The development of shape memory actuators has enabled noteworthy applications in the mechanical engineering, robotics, aerospace, and oil industries and in medicine. These applications have been targeted on miniaturization and taking full advantage of spaces. This article analyses a Ti-Ni shape memory actuator used as part of a flow control system. A Ti-Ni spring actuator is subjected to thermomechanical training and parameters such as transformation temperature, thermal hysteresis and shape memory effect performance were investigated. These parameters were important for understanding the behavior of the actuator related to martensitic phase transformation during the heating and cooling cycles which it undergoes when in service. The multiple regression methodology was used as a computational tool for analysing data in order to simulate and predict the results for stress and cycles where the experimental data was not developed. The results obtained using the training cycles enable actuators to be characterized and the numerical simulation to be validated.

Keywords: shape memory effect, Ti-Ni actuators, multiple regression analysis

1. Introduction

The development of smart actuators using a shape memory effect has become very common in industrial applications1,2. The main reason for this is that these make it possible to generate very compact actuator systems. Ti-Ni alloys are an important class of alloys, due to their shape memory phenomena (one-way shape memory effect, two-way shape memory effect, and their superelasticity)3,4. This class of material has been used as a raw material for actuator devices in mechanical engineering and many other areas, including medicine, robotics, aerospace, oil and natural gas industries.

Phase transformation fatigue and mechanical cycling fatigue are important phenomena associated to smart spring actuators. Understanding shape memory degradation involving fatigue for different levels of stress and a given number of cycles is valuable information in application such as to have when studying flow control valves. The main objective of this research is to predict phase transformation parameters by using computational numerical methods monitoring the behavior of parameters for different levels of stress and number of cycles. The numerical methods proposed in this research have already been used to describe others phenomena such as drying phenomena. Research using function fit was developed depending on different conditions such as temperature, drying speed and relative humidity presented good coherence with experimental values5. Statistical parameters such as the Correlation Coefficient, Mean Bias Error and Root Mean Square Error, were used to determine the drying rate of figs6. Another paper that tackled the same subject was written by Menges and Ertekin (2006)7. In another engineering area, Hellmann et al. (1998) put forward a method to prevent the behavior of chillers by using a simple algebraic function8. The main research result was to improve the curve fit for absorption systems9.

The results obtained address the influence of the training process on the behavior of the actuator, where modifications were observed in the critical transformation temperatures and a variation in yield from shape memory due to reconfiguration of the internal defects10-12. Data analysis using numerical methods emerged as an efficient method for projecting the behavior parameters of the actuator with non-experimentally assessed tensions. The numerical method applied made the results obtained interesting and coherent as to improving understanding of the functionality of the actuators.

2. Experimental procedure

2.1. Experimental part

This research investigated the performance of the shape memory effect on commercial Ti-50.4at% Ni, 0.89 mm in diameter wire. This material was homogenized at 500 °C for a period of 24 hours and then cooled in water at 25 °C.

Differential scanning calorimetry (DSC) technique was used to identify the direct and reverse martensitic transformation temperatures (Ms, and Mf, the starting and finishing temperatures of direct martensitic transformation on cooling from austenite to martensite, As and Af, the beginning and the end temperatures of reverse martensitic transformation on heating) in the range of temperatures between - 60 °C and 100 °C, at a rate of 10 °C/min-1. This method was used to analyze the material in the as-received and heat treated conditions.

The wire was used to produce actuator with a helical spring by shaping the wire around a screw, and then submitting the set to heat treatment at 500 °C for 24 hours13,14. The actuator obtained had a 6.0 mm outer diameter, 4 active coils and was approximately 6.0 mm long.

Figure 1 shows the schematic configuration of the system used to conduct the thermomechanical training of the actuator. The device consists of an element to fix the spring, a pulley, a non-extensible wire, a rod to transmit the applied load, the load, the Linear Variation Displacement Transducer (LVDT) sensors, the temperature sensors and a data acquisition set Monitored by a computer. The load is applied at 20 °C elongating the spring (as shown latter in the results topic, at 20 °C, the wire is in the martensitic state) and then heated to 90 °C (above Af, also shown later) contracting the spring. If this cooling and heating cycle is repeated several times, it is called training process. The heating and cooling rate was estimated to be around 5°C/min and 3°C/min, respectively.

Figure 1.
Diagram of the system used to conduct on the training process of the actuator.

The evaluation of the mechanical behavior and the shape memory effect of the actuator was performed by repeating 40 times the heating and cooling cycle for each applied load (stress) of 35, 70, 105, 135, 170 and 200 MPa15,16.

Figure 2 shows the main parameters from a generic thermo-elastic deformation curve. The tangent rule identifies the critical transformation temperatures (As, Af, Ms and Mf), thermo-elastic deformations in millimeters (εt) and the linear dislocation between cycles (X)15.

Figure 2.
Tangent rule for thermo-elastic deformation as function of temperature curves15.

2.2. Numerical part

In this section the study seeks to verify the consistency of the results regarding the behavior of the actuator for non-experimentally analyzed stresses due to the limitation of the device as to producing fast heating and cooling of the cycles.

The modeling procedure consists of a multiple regression applied to the experimental data obtained in order to predict numerical values that were not evaluated during the experimental training procedure. The setting used was a non-linear method based on Levenberg-Marquardt and Gauss-Newton. Matlab software was used as a development tool by creating subroutines to calculate the general settings17.

Initially, non-linear regressions were performed to calculate the coefficients for known stress. Thereafter, the results were evaluated by linear regression for the total adjustment, i.e., the verification of the simulated and experimental data, for each value of stress and applied number of cycles.

To determine the coefficients of the equations for each critical transformation temperature (As, Af, Ms, Mf), a non-linear regression was performed for each of the forty cycles. Then calculation of the general coefficients as a function of the each applied stress was made in order to obtain a general numerical model that would predict the temperatures evolution as function of number of cycles, i.e., As, Af, Ms and Mf.

Equations 1 to 4 were defined to simulate the behavior of the critical temperatures of the transformation as function of the cycles evaluated. The general coefficients are shown in Equations 5 to 8. The adjustment consists of coefficients that vary with the stress and the number of cycles.

where a, b, c and d, were the temperature coefficients of the equations, which were implicit functions of the tension applied to the material by means of training. A linear form was used due to the simplicity of the analysis5,6,7,18.

In order to validate the equations and the data obtained by multiple regressions method, the following statistical parameters were used:

  • Correlation coefficient (r

    2): defines the accuracy of the simulated values in relation to the experimental ones as a function of the straight line, i.e., the proximity of all values being the same for each stress condition during the evolution of the cycles

    19;

  • Mean bias error (MBE): a total residue, between the experimental and simulated values

    19;

  • Root mean square error (RMSE), the total squared residue between the experimental and simulated data

    19.

Equations 9, 10 and 11 were used to calculate the statistical parameters mentioned above.

Where:

Valuesexperimental; Values from experimental data.

Valuessimulated; Values from simulated data.

N; number of experimental data.

3. Results and discussions

3.1. Calorimetric analysis

Figure 3 shows curves representing heat flow as a function of temperature for the as-received condition. As far as no transformation peaks were observed, it is supposed that the wire is in as drawn condition, that is, work hardened condition inhibiting the reversibility of direct and reverse martensitic transformation4,20.

Figure 3.
DSC curves of the as received material.

Figure 4 shows the calorimetric curves that designate phase transformation after heat treatment at 500 °C for 24 hours. On heating, a reverse martensitic transformation peak is observed with As and Af temperatures of 56 and 73 °C, respectively. On cooling, a direct martensitic transformation temperature with Ms and Mf of 39 and 26 °C are observed. The respectives transformations enthalpies on heating and on cooling were ∆HH = 25.0 J/g and ∆HC = 26.0 J/g, a tipical values for NiTi shape memory alloys.

Figure 4.
Calorimetric curve for the sample when heat treated at 500 °C for 24 hours.

3.2. Training process of the actuators

Thermomechanical cycles were performed in a silicone oil thermal bath as shown elsewhere15. A constant shear stress was applied to the actuator wire and this enabled the recovery of the shape of the activator to be assessed as a function of the applied stress and the number of cycles. The training process was developed in this research to induce martensite variants reconfiguration in accordance to the axis along the stress is applied11,21-24.

The heating and cooling cycles brought about by training enable the critical temperature, thermal hysteresis and the thermo-elastic strain to be studied. The thermo-elastic strain is able to indicate the efficiency of the shape memory effect by means of linear displacement.

Figure 5 shows the strain curves as a function of temperature for nominal shear stress of 35, 70, 105, 135, 170 and 200 MPa.


Figure 5. Strain as a function of temperature for. a) 35 MPa e b) 70 MPa, c) 105 MPa, d) 135 MPa, e) 170 MPa and f) 200 MPa.

The analysis of the strain curves identifies that 35 MPa stress exhibited approximately 13.0 mm of thermo-elastic deformation after 40 training cycles and that stresses of 70 and 105 MPa showed a continuous trend for thermo-elastic deformation to increase. The 70 MPa applied stress exhibited a thermo-elastic increase from 26.0 to 28.0 mm and 105 MPa stress showed an increase from 35.0 to 41.0 mm. Thermo-elastic deformation behavior showed an increase for 35, 70 and 105 MPa stress. Above these parameters the tendency starts to decrease and thermal hysteresis starts to increase. Narrow thermo-elastic results were observed for 35 and 70 MPa. The main reason for this behavior was the small stress applied that was insufficient to totally induce martensitic variants orientation.

Figure 6 shows the thermo-elastic behavior as a function of the number of cycles for the actuators and Figure 7 shows the Ms and As temperatures as that for applied stress.


Figure 7. Temperature behavior as function of number of cycles to 35, 70, 105, 135, 170 and 200 MPa stress.a) Ms and b) As.

Figure 6.
Thermo-elastic deformation behavior as function of number of cycles for 35, 70, 105, 135, 170 and 200 MPa stress.

The start temperature of the martensitic transformation (Ms) tends to increase with stress in accordance with the Clausius-Clapeyron relation4,22,23. During the training process, Ms temperature decreases with number of cycles for small stress value, while for stress above 135 MPa the evolution is the inverse. In the first case, the Ms temperature decreases because of dislocation reorientation which modifies the internal stress and points the martensitic variant in the preferential direction16,24. In the second case, the Ms temperature increases for high stress due to fast martensitic reconfiguration and exhaustion (practically all martensitic variants were immediately activated by stress)23. The results for As temperature showed that this parameter decreases for all stresses studied. However, for a small stress the variation shown is almost 5 °C and for higher stresses, the variation is 12 °C.

The increase behavior observed on the Ms temperature and the decrease behavior on the As temperature results in the hysteresis reduction. This phenomenon is a result of an efficient martensitic orientation process due to the thermomechanical training process that induces more martensitic variants activation which needs less energy to initiate transformation10,11.

3.3. Modeling the temperature curves (As, Af, Ms, Mf)

The development of the curve fit in Matlab enabled the linear correlation of the statistical parameters to be obtained. The Mean Bias Error and Root Mean Square Error were obtained to validate the regressions applied to the critical temperatures of the phase transformation. Table 1 presents these parameters as a function of each of the critical transformation temperatures studied for the actuators developed.

Table 1. Results of the statistical parameters used for validating the multiple regression applied to the adjustment of the temperature curves.

As can be seen in Table 1, for all temperature data, the values of the statistical parameters enabled it to be verified that the regressions applied provided satisfactory results. The values shown in Table 1 certify the validity of this model for the non-linear regression, since it has a large r2 value, almost unity, and small values for both the mean bias error and the root mean square error.

3.3.1. Temperature (As)

Figure 8 shows the temperature curves adjusted by the regressions of the experimental values with the simulated ones for each applied stress and the results of the linear regressions as a function of the number of cycles. The experimental values are represented by solid lines and the simulated ones by symbols. In this case, the adjustments were made individually. It was observed that the behavior of the curves was similar, since the simulated values are close to the experimental ones, with a maximum relative error of 4% at 200 MPa stress. This error is due to the behavior of the equations for this load which led to the final correlations obtaining this value, but even that it is considered a very good fit.


Figure 8. (a) Comparison of the experimental values with those adjusted by regression for temperature As and (b) Linear regression of the experimental and simulated values of the temperature As for different applied tensions.

On comparing the simulated values with the general equation for temperature As (Equation 1), and the experimental values for this temperature (Figure 8) it is observed that the multiple regression methodology applied reveals that the simulated values met a 95% probability within the range limited between maximum and minimum values, thus ensuring an optimal setting according to the function of a straight line (x = y)6.

Figure 9 shows the data for the residual As temperature as a function of the simulated values based on the multiple regression of the experimental data and the comparison of these curves after the multiple regression of the data for temperature As.


Figure 9. (a) Comparison of the residual values versus the simulated values of temperature As by means of the multiple regression and (b) final comparison of the experimental and simulated values of temperature As as function of the evolution of cycles and tensions.

The data presented in Figure 9 a enable the observation to be made that the residue confirms a good fit was made - the residues from a fitted curve are defined as the differences between the experimental data and the simulated data at each load value - since this behavior does not show a systematic pattern, thus showing the efficiency of the model applied. On comparing the simulated values with the general equation for As temperature (Equation 1), and the experimental values for this temperature (Figure 9b) it is observed that the multiple regression methodology applied reveals that the simulated values met a 95% probability within the range limited between maximum and minimum values, thus ensuring an optimal setting according to the function of a straight line (x = y)6. The comparison of the residues with the simulated values provides a maximum relative error of 3% and a minimum one of 0.2%, which once again validates the model used to simulate the parameter of As temperature for the tension range between 35 and 200 MPa.

3.3.2. Temperature (Ms)

Figure 10 shows the curves of the temperature Ms adjusted by the regressions of the experimental values with the simulated ones for each stress applied and the results of the linear regressions as function of the number of cycles. The experimental values are represented by solid lines and the simulated ones by symbols. In this case, the adjustment is made individually, and it is observed that the behavior of the curves is generally similar, since the simulated values are close to the experimental ones, with a maximum relative error of 3% for 105 MPa stress. This error is due to the behavior of the equations for this load which led to the final correlations obtaining this value, but even that it is considered a very good fit.


Figure 10. (a) Comparison of the experimental values with the adjusted ones by regression for temperature Ms and (b) Linear regression of the experimental and simulated values of temperature Ms for different applied tensions.

On comparing the simulated values with the general equation of temperature Ms (Equation 2), and the experimental values of the temperature As, (Figure 10), it is observed that the multiple regression methodology applied reveals that the simulated values met a 95% probability within the range limited between maximum and minimum values, thus ensuring an optimal setting according to the function of a straight line (x = y)6.

Figure 11 shows the data of the residual temperature Ms as a function of the simulated values based on the multiple regression of the experimental data and the comparison of these curves after the multiple regression of the data for temperature Ms.


Figure 11. (a) Comparison of the residue versus the simulated values of temperature Ms by the multiple regression and (b) final comparison of the experimental and simulated values of temperature Ms as a function of the evolution of the cycles and tensions.

The data presented in Figure 11.a enable it to be observed that the residue confirms the good fit achieved, since this behavior does not show a systematic pattern, thus showing the efficiency of the applied model. The comparison of the simulated residues with the simulated values provides a maximum relative error of 6% and a minimum one of 0.6 %, which once again validates the model used for simulating of the temperature parameter AS for the tension range between 35 and 200 MPa.

3.3.3. Temperatures (Af and Mf)

The data on fitting the curves to the temperature parameters of Af and Mf showed a similar behavior to that obtained for the parameters As and Ms. However, by analyzing the comparison of the residues results, relative errors of 2% maximum and 0.13% minimum were obtained for the parameter Af. Mf had maximum and minimum errors of 7% and 0.7% respectively. These errors are due to the behavior of the equations for these loads that led to the final correlations obtained for these values, but even that it is considered a very good fit.

Figure 12 shows the result of the final comparison between the experimental and simulated values of the temperature parameters Af and Mf.


Figure 12. Final comparison of experimental and simulated values of the parameters of temperature. (a) Af and (b) Mf.

3.4. Simulation

Figures 13 and 14 show the behavior of temperatures (As, Ms, Af, Mf) for stress ranging from 35 MPa to 200 MPa and the evolution of the cycles, as per the equations provided by the multiple regression of the experimental values.


Figure 13. Behavior of the temperature parameters for stress between 35 MPa and 200 MPa. (a) As and (b) Ms.


Figure 14. Behavior of the temperature parameters for stress between 35 MPa and 200 MPa. (a) Af and (b) Mf.

Using the results of the computational simulation, Figures 13 and 14 show that the behavior of all temperatures (As, Ms, Af, Mf) follows a set pattern, thus obeying the reaction of the material to external stimuli that impose real changes on transformation temperatures. Therefore, we have as the default (standard) behavior, reducing the As temperatures with the evolution of the cycles and their increasing when using stress of a higher intensity. Thus, the Ms temperatures remain showing a mild trend of increasing with the evolution of the cycles and when tensions of higher intensity are used. This behavior is also seen at Af and Mf temperatures. Therefore, multiple regression using temperature data, cycles and tensions of thermomechanical training in a Ti-Ni actuator was able to produce results by generating information suitable for simulation testing that would take too long to be studied by means of conventional experimental analysis.

4. Conclusions

Based on the results, we conclude that:

  • This paper showed the experimental results for the training process. Transformation temperatures and thermal hysteresis were modified due to the applied stress and the evolution of the cycles. Stress fields and martensitic reconfiguration induced changes in temperature transformation due to the applied stress and number of cycle's evolution. These phenomena were observed because of the variation in temperature and thermo-elastic deformation during the training procedure;

  • The curve fitting methodology applied is valid to simulate the conditions, where the experimental data was not developed on a conventional way, of training in a Ti-Ni actuator by thermomechanical analysis, thus facilitating the characterization of phase transformation of the actuator according to the condition that was imposed; and

  • The numerical results showed that the behavior of all temperatures (A

    s, M

    s, A

    f and M

    f) follows a definite pattern and obey the material reaction to external stimuli that impose modifications on the transformation temperature.

Acknowledgements

The authors are grateful to the Brazilian agencies: Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Fundação de Amparo à Ciência e Tecnologia do Estado de Pernambuco (FACEPE) for their financial support and the scholarships they awarded to this research project.

Received: January 19, 2012

Revised: January 29, 2013

  • 1. Andersen D, Pedersen A and Sangesland SAS. Detailed study of shape memory alloys in oil well applications Trondheim: SINTEF Petroleum Research; 1999.
  • 2. Gore J, Bowles A, Maylin M, Chandrasekaran L, Forsyth D and Buyers M. High temperature shape memory alloy actuators through mechanical treatments for an oil & gas down-hole valve. Industrial and Commercial Applications of Smart Structures Technologies, SPIE Digital Library, 2008. v. 6930. http://dx.doi.org/10.1117/12.776306
  • 3. Delaey L, Krishnan RV, Tas H and Warlimont H. Thermoelasticity, pseudoelasticity and the memory effects associated with martensitic transformations. Journal of Materials Science. 1974; 9:14. http://dx.doi.org/10.1007/BF00552939
  • 4. Otsuka K and Ren X. Physical metallurgy of Ti-Ni-based shape memory alloys. Progress in Materials Science 2005; 50(5):511-678. http://dx.doi.org/10.1016/j.pmatsci.2004.10.001
  • 5. Akpinar EK and Bicer Y. Mathematical modelling of thin layer drying process of long green pepper in solar dryer and under open Sun. Energy Conversion and Management 2008; 49(6):1367-1375. http://dx.doi.org/10.1016/j.enconman.2008.01.004
  • 6. Xanthopoulos G, Oikonomou, N and Lambrinos G. Applicability of a single-layer drying model to predict the drying rate of whole figs. Journal of Food Engineering 2007. 81(3):553-559. http://dx.doi.org/10.1016/j.jfoodeng.2006.11.033
  • 7. Menges HO and Ertekin C. Mathematical modeling of thin layer drying of Golden apples. Journal of Food Engineering 2006; 77(1):119-125. http://dx.doi.org/10.1016/j.jfoodeng.2005.06.049
  • 8. Hellmann HM, Schweigler C and Ziegler F. The characteristic equations of absorption chillers. In: Proceedings of the International Sorption Heat Pump Conference; 1999; Munich. Munich; 1999. p. 4.
  • 9. Gutiérrez-Urueta G, Rodríguez P, Ziegler F, Lecuona A and Rodríguez-Hidalgo MC. Extension of the characteristic equation to absorption chillers with adiabatic absorbers. International Jounal of Refrigeration 2012; 35(3):709-718. http://dx.doi.org/10.1016/j.ijrefrig.2011.10.010
  • 10. Gonzalez CH, Oliveira CAN, Pina EAC, Urtiga Filho SL, Araújo Filho OO and Araújo CJ. Heat Treatments and Thermomechanical Cycling Influences on the R-Phase in Ti-Ni Shape Memory Alloys. Materials Research-Ibero-American Journal of Materials 2010; 13(3):325-331. http://dx.doi.org/10.1590/S1516-14392010000300008
  • 11. Paula ASC, Canejo JPHG, Martins RMS and Braz Fernandes FM. Effect of thermal cyclin on the transformation temperature ranges of a Ni-Ti shape memory alloy. Materials Science and Engineering: A. 2004; 378(1-2):92-96. http://dx.doi.org/10.1016/j.msea.2003.11.057
  • 12. Somsen C, Zähres H, Kastner J, Kakeshita T and Saburi T. Influence of thermal annealing on the martensitic transitions in Ni-Ti shape memory alloys. Materials Science and Engineering: A. 1999; 273:310-314. http://dx.doi.org/10.1016/S0921-5093(99)00362-7
  • 13. De Araújo CJ, Silva E and Gonzalez CH. Thermal alarm using a shape memory alloy helical spring. In: Proceedings of the 16th Brazilian Congress of Mechanical Engineering Uberlândia; 2005; Uberlândia. Uberlândia; 2005. v. 15, p. 8.
  • 14. Oliveira CAN, Gonzalez CH, De Araújo CJ, Pina EAC, Urtiga Filho SL and Araújo Filho OO. Caracterização do Efeito Memória de forma Reversível de Molas de Cu-Zn-Al. Revista Eletrônica de Materiais e Processos 2009; 4(3):76-86.
  • 15. Oliveira CAN, Gonzalez CH, Araújo CJ, Araújo Filho OO and Urtiga Filho SL. Thermoelastic Properties on Cu-Zn-Al Shape Memory Springs. Materials Research 2011; 13(2):219-223. http://dx.doi.org/10.1590/S1516-14392010000200016
  • 16. De Araújo CJ, Gonzalez CH, Morin M and Guénin G. Influence of the Mechanical Loading History on the Stress Assisted Two Way Memory Effect in a Ti-Ni-Cu Alloy. In: Anais do 61ş Congresso Anual da Associação Brasileira de Metalurgia e Materiais; 2006; Rio de Janeiro. Rio de Janeiro; 2006. p. 8.
  • 17. Steven C and Raymond C. Numerical Methods for Engineers McGraw-Hill College; 2009.
  • 18. Togrul IT and Pehlivan D. Mathematical modelling of solar drying of apricots in thin layers. Journal of Food Engineering 2002; 55(3):209-216. http://dx.doi.org/10.1016/S0260-8774(02)00065-1
  • 19. Chaudhuri S and Drton M. On the Bias and Mean-square error of the Sample Minimum and the Maximum Likelihood Estimator for Two Ordered Normal Means Washington: Department of Statistics, University of Washington; 2003.
  • 20. Oliveira CAN, Gonzalez CH, Pina EAC, Urtiga Filho SL, Araújo Filho OO and De Araújo CJ. Precipitates formation in Ti-Ni equiatomic alloys due to annealing heat treatment. Materials Science Forum 2010; 643:49-54. http://dx.doi.org/10.4028/www.scientific.net/MSF.643.49
  • 21. Wang ZG, Zu XT, Dai JY, Fu P and Feng XD. Effect of Thermomechanical Training Temperature on the Two-Way Shape Memory Effect of Ti-Ni and Ti-Ni-Cu Shape Memory Alloys Springs. Materials Letters 2003; 57(9-10):1501-1507. http://dx.doi.org/10.1016/S0167-577X(02)01014-5
  • 22. Huang W and Goh HB. On the long-term stability of two-way shape memory alloy trained by reheat treatment. Journal of Materials Science Letters 2001; 20(19):1795-1797. http://dx.doi.org/10.1023/A:1012543402147
  • 23. Lahoz R, Gracia-Villa L and Puertolas JA. Training of the two-way shape memory effect by bending in NiTi alloys. Journal of Engineering Materials and Technology-Transactions of the Asme 2002; 124(4):397-401. http://dx.doi.org/10.1115/1.1495001
  • 24. Villafana A, Masse M, Portier R and Pons J. Thermomechanical testing machine conceived for the study of shape memory alloys. Application to the training and testing of the two-way memory effect. Journal de Physique IV France 1997; 7(C5):655-660. http://dx.doi.org/10.1051/jp4:19975104
  • *
    e-mail:
  • Publication Dates

    • Publication in this collection
      10 May 2013
    • Date of issue
      Aug 2013

    History

    • Received
      19 Jan 2012
    • Reviewed
      29 Jan 2013
    ABM, ABC, ABPol UFSCar - Dep. de Engenharia de Materiais, Rod. Washington Luiz, km 235, 13565-905 - São Carlos - SP- Brasil. Tel (55 16) 3351-9487 - São Carlos - SP - Brazil
    E-mail: pessan@ufscar.br