A Thermo-mechanical Behavior Simulation of a Niti Staple Used for the Correction of Idiopathic Scoliosis

In order to treat scoliosis, a thermo-mechanical behavior study of shape memory alloy bone staple has been proposed. A pre-stretched shape-memory-alloy bone staple, which has been heated above the transition temperature, provides greater compressive force upon insertion between two vertebrae and allows the control of idiopathic scoliosis development. Until now, the optimal design has not been reached due to the lack of appropriate design tools for shape memory alloy devices. In this paper, a shape memory alloy bone staple model is proposed by developing a user subroutine UMAT based on Boyd et al. [1] unified thermodynamic shape memory alloy constitutive law using the finite element analysis software, ABAQUS. The numerical results for superelastic and shape memory effect under the tensile and three-point bending tests are presented. Simulations of the shape memory effects and force generation of the shape memory alloy staple are also shown.


Introduction
Shape memory alloys (SMAs) show unusual but intelligent material behavior such as superelaticity (ability to recover large elastic strain) and shape memory effect (ability of a deformed SMA element to recover its initial shape when it is heated) [2].These effects are, in general, interpreted as martensitic transformations between a more ordered crystal phase, called austenite, and a less ordered crystallographic phase, called martensite [3,4].These inherent and induced characteristics of SMA have given rise to their usage in a wide array of applications.Among several other applications, SMA are used to produce shrink fit rings, actuators in the form of torque beams or wires (Bo et al., [5], Trochu et al. [6,7]), orthodontic wires and other medical devices (Auricchio et al. [8][9][10]).SMA also are useful as dampers, due to the energy dissipated in the pseudo-elastic hysteresis (Qidwai et al. [11] and Boyd et al. [1,12]).
Despite several successful applications of SMA devices, SMA is currently not of large industrial use [13].Traditionally, new products were developed by prototyping and evaluation.However, this process is very time consuming and often does not fully reveal the potential failures.Finite element analysis (FEA) can greatly reduce testing and marketing time by allowing the designer to develop and simulate products prior to prototype fab-rication.Research on modeling shape memory alloys originates from the interest expressed by several companies for new medical applications as vascular stents, osteosynthesis staples, Harrington wires for scoliosis correction, [14] and orthodontic wires [15].Current commercially available software do not support shape memory effect material properties in 3D, and have limited support on superelasticity effect.Models that were developed to represent SMA behaveior mostly are based on experimental data [16].Constitutive laws are classified under three families: thermodynamical, micromechanical, and mesomechanical models.
The thermodynamical models are based on macroscopic experimentation and identification of the principal thermomechanical parameters.The advantages of these models are simplicity, appropriateness for numerical computing, and practicality for industrial applications.However, the validation is difficult due to the lack of the multiaxial experimental data.The differences between published thermodynamical models are the variation in the internal state variables and the kinetics equation.Liang and Rogers [17] have proposed a multi dimensional constitutive model using the marten-site volume fraction as internal state variable.Lagoudas et al. [18].Auricchio et al. [10,[19][20][21][22] and Qidwai [11] have developed a unified thermodynamic constitutive law and have demonstrated the unification of models by optimizing the function of the elastic strain energy.
Micromechanical models are based on researchers following a micromechanics approach, Tanaka [16], Patoor [23]; Brinson et al. [24], Lexcellent et al. [25,26]; Huang et al. [27]; Gao et al. [28], and Sun et al. [29] have attempted to closely follow the crystallographic phenomena within the material, using thermodynamics laws to describe the transformation.These models consider the martensitic variant as a transforming inclusion and use micromechanics to calculate the interaction energy due to the phase transformation in the material.Stresses and strains are obtained as volume averages calculated over a volume in which many inclusions are considered, representing the possible variants.Models based on this approach are much more complicated than thermodynamic models and are usually computationally demanding.On the other hand, being based on an accurate definition of the material properties, these models seem to offer the most rational method to derive a highly accurate three-dimensional constitutive law.
Mesomechanical model, as proposed by Kafka [30][31][32] represent an approach, that combines the two proceedings.While taking advantage from the accuracy that offers a description of the is adopted maintaining its appropriateness for numerical computing.In this model the shape memory alloys are considered heterogeneous material, where heterogeneous implies that the material is composed of different atomic components.
The goal of this present paper is not to propose a new constitutive model for shape memory alloy but to focus on the development of a robust and efficient integration algorithm for Boyd et al. [12] unified thermodynamic model, and the integration of 3D model to simulate and to invetigate the superelasticity and the shape memory effect as proven under the tensile and three point bending test of the orthopaedic staple.

Description of the Model
The model used here consists of a mechanical law, governing the stress-strain behavior, and a kinetic law governing the crystallographic transformation.These two relationships are coupled (stress necessarily feeds into the kinetic laws and the evolving phase fraction affects the stress-strain behavior. Total strain (1) is the sum of the elastic and plastic strain due to martensite transformation A mechanical law relates the strain to the stress, temperature, and transformation strain (2) as The effective compliance tensor is defined as being the inverse of elastic rigidity tensor 1   S C where C and the thermal expansion coefficient tensor,  are given according to the martensitic volume fraction  as follows:   Based on this assumption, the following relations are induced to relate the evolution of transformation strain tensor to the evolution of the martensitic volume fraction ( 5) where the transformation tensor  is given by the following formulae: max tr  is the maximum uniaxial transformation.The activation of the transformation phase from austenite to martensite and vice versa is verified by the thermodynamic conjugate force, to : s 0 and u 0 are effective specific entropy at reference condition and effective specific internal energy at the reference condition, respectively.The value of the volume fraction  , is given by the condition = 0 through direct or inverse transformations.

π 
The hardening function   f  , is responsible for the transformation induced strain hardening in the SMA material and given by

Implementation of the Theoretical Model into a Nonlinear Finite Element Model
SMA finite element model takes into account geometric and material nonlinearities.The Newton Raphson method was used by transfoming the material and geometrical non-linear equations into a sequence of linear equations solved numerically.The return mapping algorithm for the elastic prediction and transformation correction has already been adopted by many researchers [1,11,12,18,33].The implementation of the SMA constitutive equations proposed in the framework of this project was done as a user defined material.Numerical simulations was carried out using the commercially available nonlinear finite element code ABAQUS [10,11,34] on a parallel computer.The incremental (discretized) SMA constitutive model based algorithm was implemented in the user supplied subroutine UMAT [7].This algorithm was divided into two steps: 1) A thermo elastic predictor step to preview the direction of a strain and stress variation: A transformation corrector step to compute transformation strain, correct stress value, and deduce tangent stiffness tensor: 0 The material properties described were used to simulate and validate the SMA behavior.
We have considered several examples to assess the performance of the proposed model, simulating the be-havior of SMA structure.In the following, all tests were simulated starting from an initial temperature T 0 = 294 K (where the temperature is greater than the temperature from the finish of austenite transformation A f ).

Uniaxial Response: Tension Test
An uniaxial simulation was used to show the ability of the constitutive model to properly reproduce the superelastic and shape memory effect.A three-dimensional beam with 5 mm sides and a length of 100 mm was simulated where a finite element mesh has been constructed with 80 quadratic hexahedral elements each having 20 nodes.The first extremity of the beam was fixed and a simple tensile was applied on the opposing extremity (Figure 1).

Supereleasticiy Effect
The SMA beami s subjected to a simple tensile loading-unloading with a constant temperature (T = 294 K, T < A s ).The sample has developed a recoverable deformation higher than 6% as shown in the stress-strain response in Figure 2. The result obtained from a uniaxial test induced a full transformation.It is important to note the ability of the model to reproduce the superelastic effect in a uniaxial state of stress.

Shape Memory Effect
Keeping the temperature fixed at T = 294 K (T < A s ), a complete loading and unloading cycle in terms of axial strain is achieved.Holding the axial deformation constant (ε = 0%), we initially cool the sample from T = K to T = 291 K. Keeping the temperature fixed (T = K), a complete loading-unloading cycle is also carried  out.Finally, keeping the curvature fixed, specimen is heated to back up to 320 K.The beam had remained permanently deformed due martensite's presence.With heating at T = 320 K (T > A f ), the beam recovered its initial form and the permanent deformation had disappeared.
Figure 3 shows the material response in terms of stress, strain, temperature, and the capacity simulation to model the unidirectional shape memory effect.We can observe the ability of the specimen to recover its initial state during the final heating process, that is, its ability to reproduce the shape-memory effect during bending, and reaching a temperature at above T = 320 K.

Three-Point Bending Test
The second series of tests were to investigate the response of a SMA beam undergoing pure bending.The geometry of the beam is chosen to be of 2 mm by 1 mm sides having a length of 50 mm.A finite element mesh has been constructed using quadratic hexahedral elements with 5 elements for thickness, 4 elements for width, and 30 elements for length.The extremities of the beam were supported and a force of 250 N was applied in the middle (Figure 4(a)).
An investigation of the structural response of the superelastic and the shape memory effect has been once again carried out, (Figure 4(b)).

Supereleasticiy Effect
Keeping the temperature constant (T = 294 K), the beam was subjected to a loading and unloading solicitation which is equivalent to apply a bending force.Figure 5 shows the superelastic behavior for applied force versus displacement of numerical model.The ability of the model to recover completely from the deformation during unloading at a temperature above T (A s ), is remarkable as well as its ability to reproduce the superelastic effect under bending tests.

Shape Memory Effect
Keeping the curvature fixed, the beam has been initially cooled from T = 295 K to T = 271 K.We later kept the temperature fixed (T = 271 K) to perform a complete loading-unloading cycle.Finally, keeping the curvature fixed, the specimen was heated back up to 320 K.

Three-Dimensional Orthopaedic Staple
In order to take advantage of the SMA staples' superelestic and shape memory behavior for the orthopedic scoliosis correction, it is important to perform a finite element modeling in order to determine the limitation of the staple before its validation.For this task, a mesh consisting of 4620 hexahedral elements and 516 tetrahedral elements (Figure 7) was build up.The properties of the materials used were the same as those used in previous tests (able 1) and Figure 8.  Table 1.Shape memory alloy materials parameters [11].has succeeded in reproducing the staple's shape memory effect.After rendering the two prongs normal to the staple body and removing, the applied load the martensite production is noticed due to exceeding the constraint limits at the initial transformation.The staple undergoes a permanent deformation transformation which prevented the legs from returning to its initial shape.This shape can only be recovered from heating.

Material Constants Values Model Variables
The thermomechanical loading history applied to the staple is as follows and consists of three steps: 1) Keeping the temperature fixed, T = 294 K (T < A s ), to open the prongs of the staple and applying a pressure on the extremities of the staple.
2) Unloading or removing the two previous forces.
3) Heating the staple heated up to T = 320 K (T > A f ).

Shape Memory Effect
Notice that the presence of residual constraints after unloading will shift the final transformation temperature to austenite.This indicates that increasing the temperature to A f is insufficient to recover the initial form, and therefore, must be heated to a temperature much greater until the all the martensite vanishes.
During the first test, the thermomechanic behavior of the free staple is observed.

Conclusions
The template is designed so that author affiliations are we have proposed a computation tool for the design of SMA-based devices.The tool is based on two main compounds: 1) A three dimensional constitutive model able to reproduce the basic macroscopic features of shape memory material, such as the superelasticity, the shape memory behaviour, and different response under tension (and compression); 2) The development of a robust algorithm solution and the integration into a 3D model.
Another aspect of crucial importance, is that the complex thermomechanical loading history should be taken into account, as for the case of the SMA structure.
To establish the viability of the proposed approach we have simulated the superelasticity and the shape memory effect under the tensile and three point bending test of the orthopedic staple.
The numerical investigations demonstrated that the developed formulation is a powerful tool to study the 3D model of the shape memory alloy devices in biomedical application.

Figure 1 .
Figure 1.Three-dimensional beam mesh for tensile test.Figure2.Superelastic effect for the tensile load.

Figure 2 .
Figure 1.Three-dimensional beam mesh for tensile test.Figure2.Superelastic effect for the tensile load.

Figure 4 .
Figure 4. (a) Three-dimensional beam mesh for three-point bending; (b) The structural response of the superelastic and the shape memory effect.

Figure 9 Figure 10 .
Figure 10.Stress generated after binding the legs and heating.(a) Stress distribution on the legs; (b) Variation of the total stress applied by the legs.

Figure 9 .
Figure 9. Shape memory effect of the staple with a free legs loading.(a) Loading; (b) Unloading; and (c) Heating.