This paper aims to build a constitutive model intended to describe the thermomechanical behavior of shape memory alloys. This behavior presents many facets, among them we have considered the simple way of shape memory, which is one of most important properties of shape memory alloys. Because of numerous stages of this effect, the subject was divided into three independent parts. For each part, we built the corresponding thermodynamic potential and we deduced the constitutive equations. To make this model workable, we have developed an algorithm. The simulation was performed using the NiTi as shape memory alloy.
The known behavior of conventional materials has allowed their use in many applications but the discovery of new properties, coming from a singular behavior of materials known as shape memory alloys, opened a way for other applications from medical to aerospace.
This unusual behavior has attracted a significant attention of scientists and researchers. Therefore, various models were proposed.
These models are based on thermodynamics laws and frameworks theories as generalized standard materials. Halphen and Nguyen [
These models can be classified as follows:
1) Macroscopic models: Built of the thermomechanical behavior, they are generally simpler in formulation and geared more towards engineering applications. Thus, the detailed physics of phase transformation are usually not rigorously addressed.
2) Micromechanical models are based on the micromechanics of a single crystal. Chu and James [
In this paper, we focus on the effect of simple way shape memory which is the first property discovered. The adopted approach is based on the sharing of the study on three following steps:
1) Orientation of twinned Martensite;
2) Heating for austenite;
3) Cooling for detwinned Martensite.
The thermomechanical cycle to the memory effect of a simple shape is defined the following:
1) Applying a mechanical load under Temperature T1 lower than (Temperature of transformation start of Martensite): The material is deformed first elastically, followed by an important deformation due to twinned Martensite orientation. When the load is cancelled the deformation is not fully recovered, only elastic part is recovered.
2) Heating to a temperature above (Temperature of transformation finish of Austenite): When the temperature reaches (Temperature of transformation start of Austenite) the deformation begins to recover.
3) Cooling to T1: Austenite begins to transform to Martensite and finally the material is at the origin of the cycle.
The study will concern three regions which are planes (ε,
σ); (ε, T) and axe T:
1) Plane (ε, σ): Region of mechanical loading;
2) Plane (ε, T): Region of thermal loading (heating);
3) Axe T: Region of cooling.
σs: Stress of orientation start of twinned Martensite;
σf: Stress of orientation finish of twinned Martensite;
ε0: Maximum deformation of orientation.
In each case, we consider an elementary volume
where is volume of parent phase and is volume of incipient phase.
is fraction of incipient phase.
Parent phase is twinned Martensite and incipient phase is detwinned Martensite Free energy of Gibbs:
: Young modulus of MartensiteB, C: Coefficients to be determined by tests,
: Elastic energy,
: Energy of deformation due to transformation,
: Free energy of phase change,
: Energy of interaction between two phases.
Assuming that the dissipation occurs during processing Clausius inequality can be written as following:
Let us write the term:
which is the driving force.
The dissipative force can be written:
The transformation occurs when this condition is satisfied:
It will result:
Let us use the function:
To write the evolution of the fraction we use the following expression:
The coefficients a1 and C can be determined using limit values of when:
It is denoted that parameter B is identified using pseudoelasticity test.
Parent phase is detwinned Martensite and incipient phase is Austenite Following the same approach:
Parameter H can also be determined by pseudoelasticity test; the coefficients and are defined using limit values of when:
The evolution of the fraction of Martensite can be expressed by:
Parent phase is Austenite and incipient phase is twinned Martensite
: Parameters to be defined by test using following conditions:
Using K.L. Ng. and al work [
Temperature of orientation loading T1 = 230 K; mechanical load;
3. Results (Figures 4-6)
DiscussionsThe results obtained coincide well with the experimental
data in the case of orientation test, in case of heating to make martensite changing into austenite which occurs with recovering previous deformation and finally changing of martensite into austenite during the cooling. It seems that our results are consistent with practical data of the selected material (NiTi).
In each case we see clearly that deformation is function of transformation of martensite. The cycle beginning from lower temperature reaches the final temperature with recovering the deformation caused by applied mechanical load.
This constitutive model presented in this paper was built using a simple formalism. We have divided the work into three parts, and each one has its proper particularities after simulating the model. The results appear in good agreement with experimental data. This model can be used in engineering fields.