<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">MNSMS</journal-id><journal-title-group><journal-title>Modeling and Numerical Simulation of Material Science</journal-title></journal-title-group><issn pub-type="epub">2164-5345</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/mnsms.2013.34017</article-id><article-id pub-id-type="publisher-id">MNSMS-37532</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject></subj-group></article-categories><title-group><article-title>
 
 
  One Dimensional Modeling of the Shape Memory Effect
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>elkacem</surname><given-names>Meddour</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hamma</surname><given-names>Zedira</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hamid</surname><given-names>Djebaili</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Laboratory of LaSPI2A, Department of Sciences &amp;amp; Technology, University Abbas Laghrour, 
Khenchela, Algeria</addr-line></aff><aff id="aff1"><addr-line>Department of Mechanical Engineering, University of Batna, Batna, Algeria</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>samsum66@gmail.com(EM)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>04</day><month>10</month><year>2013</year></pub-date><volume>03</volume><issue>04</issue><fpage>124</fpage><lpage>128</lpage><history><date date-type="received"><day>February</day>	<month>16,</month>	<year>2013</year></date><date date-type="rev-recd"><day>March</day>	<month>16,</month>	<year>2013</year>	</date><date date-type="accepted"><day>March</day>	<month>24,</month>	<year>2013</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   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.
     
 
</p></abstract><kwd-group><kwd>Shape; Strain; Detwinned Martensite; Region</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>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.</p><p>This unusual behavior has attracted a significant attention of scientists and researchers. Therefore, various models were proposed.</p><p>These models are based on thermodynamics laws and frameworks theories as generalized standard materials. Halphen and Nguyen [<xref ref-type="bibr" rid="scirp.37532-ref1">1</xref>] used by Lexcellent and Licht [<xref ref-type="bibr" rid="scirp.37532-ref2">2</xref>], Edelen’s formalism [<xref ref-type="bibr" rid="scirp.37532-ref3">3</xref>] used by Tanaka and Nagaki [<xref ref-type="bibr" rid="scirp.37532-ref4">4</xref>].</p><p>These models can be classified as follows:</p><p>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.</p><p>2) Micromechanical models are based on the micromechanics of a single crystal. Chu and James [<xref ref-type="bibr" rid="scirp.37532-ref5">5</xref>]; James et al., [<xref ref-type="bibr" rid="scirp.37532-ref6">6</xref>]; Lexcellent et al., [<xref ref-type="bibr" rid="scirp.37532-ref7">7</xref>]; Govindjee and Hall [<xref ref-type="bibr" rid="scirp.37532-ref8">8</xref>]; Berveiller et al. [<xref ref-type="bibr" rid="scirp.37532-ref9">9</xref>]; Patoor et al. [<xref ref-type="bibr" rid="scirp.37532-ref10">10</xref>]. Generally, this class of model was derived from more fundamental thermodynamic principles and the shape strains of different martensite variants are included.</p><p>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:</p><p>1) Orientation of twinned Martensite;</p><p>2) Heating for austenite;</p><p>3) Cooling for detwinned Martensite.</p></sec><sec id="s2"><title>2. Methods</title><sec id="s2_1"><title>2.1. Presentation of the Subject (<xref ref-type="fig" rid="fig1">Figure 1</xref>)</title><p>The thermomechanical cycle to the memory effect of a simple shape is defined the following:</p><p>1) Applying a mechanical load under Temperature T<sub>1</sub> lower than <img src="5-2190040\4707946b-a151-4549-85c1-fa1d7b979dfb.jpg" /> (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.</p><p>2) Heating to a temperature above <img src="5-2190040\5758c5fd-0054-49fc-b31e-281110379a0f.jpg" /> (Temperature of transformation finish of Austenite): When the temperature reaches <img src="5-2190040\3307109b-fcf3-47b0-82cd-37d04f31a504.jpg" /> (Temperature of transformation start of Austenite) the deformation begins to recover.</p><p>3) Cooling to T<sub>1</sub>: Austenite begins to transform to Martensite and finally the material is at the origin of the cycle.</p></sec><sec id="s2_2"><title>2.2. Definition of Regions (<xref ref-type="fig" rid="fig2">Figure 2</xref>)</title><p>The study will concern three regions which are planes (ε,</p><p>σ); (ε, T) and axe T:</p><p>1) Plane (ε, σ): Region of mechanical loading;</p><p>2) Plane (ε, T): Region of thermal loading (heating);</p><p>3) Axe T: Region of cooling.</p><p>σ<sub>s</sub>: Stress of orientation start of twinned Martensite;</p><p>σ<sub>f</sub>: Stress of orientation finish of twinned Martensite;</p><p>ε<sub>0</sub>: Maximum deformation of orientation.</p></sec><sec id="s2_3"><title>2.3. Constitutive Equations (<xref ref-type="fig" rid="fig3">Figure 3</xref>)</title><p>In each case, we consider an elementary volume</p><p><img src="5-2190040\f527badf-0107-4144-aa47-ba68be6f6ce0.jpg" /></p><p>where <img src="5-2190040\22520933-afda-49e8-8d51-420ab674703c.jpg" /> is volume of parent phase and <img src="5-2190040\a21482d8-46e0-4508-8ad3-94d1afc39ef2.jpg" /> is volume of incipient phase.</p><p><img src="5-2190040\2102c44a-fc57-4b6e-bbde-efa0afa2cebc.jpg" /></p><p>is fraction of incipient phase.</p><sec id="s2_3_1"><title>2.3.1. Region 1</title><p>Parent phase is twinned Martensite and incipient phase is detwinned Martensite Free energy of Gibbs:</p><disp-formula id="scirp.37532-formula108082"><label>(1)</label><graphic position="anchor" xlink:href="5-2190040\6b98e609-af20-4c53-97c6-d342c30a8446.jpg"  xlink:type="simple"/></disp-formula><p><img src="5-2190040\5704c980-ddce-45e1-a215-cb3f23884b75.jpg" />: Young modulus of MartensiteB, C: Coefficients to be determined by tests,</p><p><img src="5-2190040\18cc9df4-7037-428c-955b-aff795283bb9.jpg" />: Elastic energy,</p><p><img src="5-2190040\95687a9d-41ee-4381-aaf6-eab2740fba48.jpg" />: Energy of deformation due to transformation,</p><p><img src="5-2190040\7a24b308-1728-4c32-b1f2-ec88b7628c14.jpg" />: Free energy of phase change,</p><p><img src="5-2190040\54399aed-e199-42a6-a555-8ff0b5735710.jpg" />: Energy of interaction between two phases.</p><p>Assuming that the dissipation occurs during processing Clausius inequality can be written as following:</p><disp-formula id="scirp.37532-formula108083"><label>(2)</label><graphic position="anchor" xlink:href="5-2190040\9938dbe5-bfa4-4710-b6bd-5f26344d2c08.jpg"  xlink:type="simple"/></disp-formula><p>Let us write the term:</p><p><img src="5-2190040\c7f4dcec-dd80-49d4-9f5a-13181cea2537.jpg" /></p><p>which is the driving force.</p><p>The dissipative force can be written:</p><disp-formula id="scirp.37532-formula108084"><label>(3)</label><graphic position="anchor" xlink:href="5-2190040\161bfd5b-9eda-4ce0-b034-1f5fd6573860.jpg"  xlink:type="simple"/></disp-formula><p>The transformation occurs when this condition is satisfied:</p><disp-formula id="scirp.37532-formula108085"><label>(4)</label><graphic position="anchor" xlink:href="5-2190040\c6c08f1b-f9ea-498e-b465-8adce0e0b4fc.jpg"  xlink:type="simple"/></disp-formula><p>It will result:</p><disp-formula id="scirp.37532-formula108086"><label>(5)</label><graphic position="anchor" xlink:href="5-2190040\af5162c4-4bdb-4234-81fb-ab51be3d8fde.jpg"  xlink:type="simple"/></disp-formula><p>Let us use the function:</p><disp-formula id="scirp.37532-formula108087"><label>(6)</label><graphic position="anchor" xlink:href="5-2190040\8a6d7259-098d-4a4c-8e86-dd198df0e912.jpg"  xlink:type="simple"/></disp-formula><p>To write the evolution of the fraction we use the following expression:</p><disp-formula id="scirp.37532-formula108088"><label>(7)</label><graphic position="anchor" xlink:href="5-2190040\685d12fe-edbd-4275-86c2-635287ca43af.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108089"><label>(8)</label><graphic position="anchor" xlink:href="5-2190040\8da38b5c-ad9b-4011-a238-48888bfe2f68.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108090"><label>(9)</label><graphic position="anchor" xlink:href="5-2190040\25da8c59-bba3-4e37-9d6a-3367177b7dae.jpg"  xlink:type="simple"/></disp-formula><p>The coefficients a<sub>1</sub> and C can be determined using limit values of <img src="5-2190040\30f9438a-7f53-4e2c-a1eb-570d994f8935.jpg" /> when:</p><disp-formula id="scirp.37532-formula108091"><label>(10)</label><graphic position="anchor" xlink:href="5-2190040\cba90ccd-6307-409d-a866-a6e74e990ffc.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108092"><label>(11)</label><graphic position="anchor" xlink:href="5-2190040\c8c351c6-d448-4bbf-8a3d-5f39b62ca013.jpg"  xlink:type="simple"/></disp-formula><p>It is denoted that parameter B is identified using pseudoelasticity test.</p></sec><sec id="s2_3_2"><title>2.3.2. Region 2</title><p>Parent phase is detwinned Martensite and incipient phase is Austenite Following the same approach:</p><disp-formula id="scirp.37532-formula108093"><label>(12)</label><graphic position="anchor" xlink:href="5-2190040\2b4090f7-2973-4691-9256-fac2c504553c.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108094"><label>(13)</label><graphic position="anchor" xlink:href="5-2190040\2c2249a2-cd2c-4f8a-85e5-76ec3ea1430e.jpg"  xlink:type="simple"/></disp-formula><p>Parameter H can also be determined by pseudoelasticity test; the coefficients <img src="5-2190040\802c6285-76db-42e3-a07d-e39f95c298db.jpg" /> and <img src="5-2190040\50f29cd1-f1e5-4eee-acce-15f02d251774.jpg" /> are defined using limit values of <img src="5-2190040\cd08520c-2d9c-4a17-b657-8b6f72e0ed48.jpg" /> when:</p><disp-formula id="scirp.37532-formula108095"><label>(14)</label><graphic position="anchor" xlink:href="5-2190040\4fed931f-9919-4260-9e44-3f1c132303f6.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108096"><label>(15)</label><graphic position="anchor" xlink:href="5-2190040\a47cb15e-0c57-4841-98ac-cf013bdb83de.jpg"  xlink:type="simple"/></disp-formula><p>The evolution of the fraction of Martensite can be expressed by:</p><disp-formula id="scirp.37532-formula108097"><label>(16)</label><graphic position="anchor" xlink:href="5-2190040\cab92eea-a8af-4e25-ac65-fd05b6049604.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s2_3_3"><title>2.3.3. Region 3</title><p>Parent phase is Austenite and incipient phase is twinned Martensite</p><disp-formula id="scirp.37532-formula108098"><label>(17)</label><graphic position="anchor" xlink:href="5-2190040\132f8201-3f8d-48db-b4b6-21f6c652ecc7.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108099"><label>(18)</label><graphic position="anchor" xlink:href="5-2190040\02c1652c-d82d-42af-9129-91abf43127dd.jpg"  xlink:type="simple"/></disp-formula><p><img src="5-2190040\b9cad321-1498-4919-a730-6b071b4f3d8a.jpg" />: Parameters to be defined by test using following conditions:</p><disp-formula id="scirp.37532-formula108100"><label>(19)</label><graphic position="anchor" xlink:href="5-2190040\ee475c66-8fd8-488f-bb12-46c46fe69610.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108101"><label>(20)</label><graphic position="anchor" xlink:href="5-2190040\1ae923ee-2923-4b36-8a61-12391be2eac2.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.37532-formula108102"><label>(21)</label><graphic position="anchor" xlink:href="5-2190040\7644dec7-2536-4f92-9704-becc611b862b.jpg"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s2_4"><title>2.4. Application of the Model</title><p>Using K.L. Ng. and al work [<xref ref-type="bibr" rid="scirp.37532-ref11">11</xref>], where a tensile test was performed on NiTi and the following parameters were determined:</p><p><img src="5-2190040\29de6c9f-2caa-492c-b44e-f1fe3f97b702.jpg" /></p><p>Temperature of orientation loading T<sub>1</sub> = 230 K; mechanical load<img src="5-2190040\bc829205-365a-4cd6-8cbf-8b113d7c5f83.jpg" />;</p><p>3. Results (Figures 4-6)</p>Discussions<p>The results obtained coincide well with the experimental</p><p>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).</p><p>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.</p></sec></sec><sec id="s3"><title>4. Conclusion</title><p>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.</p></sec><sec id="s4"><title>REFERENCES</title></sec><sec id="s5"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.37532-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">B. Halphen and Q. S. Nguyen, “Sur les Matériaux Standards Generalisés,” Journal de Mécanique, Vol. 14, 1975, pp. 39-63.</mixed-citation></ref><ref id="scirp.37532-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">C. Lexcellent and C. Licht, “Some Remarks on the Modelling of the Thermomechanical Behavior of Shape Memory Alloys,” Journal de Physique IV, Vol.1, No. C4, 1991, pp. C4-35-C4-39.</mixed-citation></ref><ref id="scirp.37532-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">D. C. Edelen, “On the Characterization of Fluxes in Nonlinear Irreversible Thermodynamics,” International Journal of Engineering Science, Vol. 12, No. 5, 1974, pp. 397-411. http://dx.doi.org/10.1016/0020-7225(74)90050-0</mixed-citation></ref><ref id="scirp.37532-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">K. Tanaka and S. Nagaki, “A Thermomechanical Description of Materials with Internal Variables in the Process of Phase Transformations,” Ingenieur-Archiv, Vol. 51, No. 5, 1982, pp. 287-299.  
http://dx.doi.org/10.1007/BF00536655</mixed-citation></ref><ref id="scirp.37532-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">C. Chu and R. D. James, “Analysis of Microstructures in Cu-14.0%Al-3.9%Ni by Energy Minimization,” Proceedings of the ICOMAT-95, Vol. 5, No. C8, 1995, pp. C8- 143-C8-149.</mixed-citation></ref><ref id="scirp.37532-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">R. D. James, R. V. Kohn and T. W. Shield, “Modeling of Branched Needle Microstructures at the Edge of a Martensite Laminate,” Proceedings of the ICOMAT-95, Vol. 5, No. C8, 1995, pp. C8-253-C8-259.</mixed-citation></ref><ref id="scirp.37532-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">C. Lexcellent, B. C. Goo, Q. P. Sun and J. Bernardint, “Characterization, Thermomechanical Behaviour and Micromechanical-Based Constitutive Model of Shape-Memory Cu-Zn-Al Single Crystals,” Acta Materialia, Vol. 44, No. 9, 1996, pp. 3773-3780.  
http://dx.doi.org/10.1016/1359-6454(95)00452-1</mixed-citation></ref><ref id="scirp.37532-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">S. Govindjee and G. J. Hall, “A Computational Model for Shape Memory Alloys,” International Journal of Solids and Structures, Vol. 37, No. 5, 2000, pp. 735-760.  
http://dx.doi.org/10.1016/S0020-7683(99)00048-7</mixed-citation></ref><ref id="scirp.37532-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">M. Berveiller, E. Patoor and M. Buisson, “Thermomechanical Constitutive Equations for Shape Memory Alloys,” Journal de Physique IV, Vol. 1, No. C4, 1991, pp. C4-387-C4-396.</mixed-citation></ref><ref id="scirp.37532-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">E. Patoor, A. Eberhardt and M. Berveiller, “On Micromechanics of Thermoelastic Phase Transition,” Proceedings of Plasticity 93: The 4th International Symposium on Plasticity and Its Applications, Baltimore, 19-23 July 1993.</mixed-citation></ref><ref id="scirp.37532-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">K. L. Ng and Q. P. Sun, “Stress-Induced Phase Transformation and Detwinning in NiTi Polycrystalline Shape Memory Alloy Tubes,” Mechanics of Materials, Vol. 38, No. 1-2, 2006, pp. 41-56.  
http://dx.doi.org/10.1016/j.mechmat.2005.05.008</mixed-citation></ref></ref-list></back></article>