Open access peer-reviewed chapter

Shape Memory Wires in R3

Written By

Shinya Okabe, Takashi Suzuki and Shuji Yoshikawa

Submitted: 23 November 2016 Reviewed: 12 April 2017 Published: 20 September 2017

DOI: 10.5772/intechopen.69175

From the Edited Volume

Shape Memory Alloys - Fundamentals and Applications

Edited by Farzad Ebrahim

Chapter metrics overview

1,270 Chapter Downloads

View Full Metrics


We propose a new model describing the dynamics of wire made of shape memory alloys, by combining an elastic curve theory and the Ginzburg-Landau theory. The wire is assumed to be a closed curve and is not to be stretched with deformation. The derived system of nonlinear partial differential equations consists of a thermoelastic system and a geometric evolution equation under the inextensible condition. We also show that the system has dual variation structure as well as a straight material case. The structure implies stability of infinitesimally stable stationary state in the Lyapunov sense.


  • shape memory alloys
  • elastic curve
  • thermoelastic system
  • nonlinear partial differential equations
  • Ginzburg-Landau theory
  • phase transition
  • stability
  • dual variation principle

1. Introduction

Shape memory effect arises from the phase transition of lattice structure. Although there are many models for shape memory alloys, one of the classical model is proposed by Falk, which we call the Falk model. Falk applied the Ginzburg-Landau theory for phase transition to shape memory alloys by regarding shear strain ϵ as an order parameter (see e.g., [1]). That is, the Helmholz free energy density proposed by Falk is given by


where u and θ are displacement and absolute temperature, respectively, and a positive constant θc denotes the critical temperature of the phase transition. We call the first term |x2 u|22 curvature energy density, the second term f(ϵ, θ) nonlinear elastic energy density and the third term f0(θ) thermal energy density. In other words, Falk represented the phase transition by using the form of the nonlinear elastic energy density. We also remark that for simplicity, all physical constants without the critical temperature are normalized by unity. The Falk model was proposed for straight materials. In the model, the material is built up by a stack of layers parallel to the so-called habit plane (see [2]) and assumed that the displacement u in that direction to depend only on a coordinate x perpendicular to the habit plane, that is, the variable x runs in the stacking direction. Then, the material conserves its volume.

From the point of the small deformation theory, we may use a linearized approximation relation ϵ = ∂xu. Moreover, we take f0 as the following typical form:

f0(θ) :=θθlogθ.

Then following a standard procedure of derivation of thermoelastic system (see e.g., [3]), we can derive the system of nonlinear partial differential equations called the Falk model:

{t2u +x4u=x{(xu)5(xu)3+(θθc)xu},tθx2θ=θxu txu.E2

Here, unknowns are displacement u and absolute temperature θ, and ∂x and ∂t represent partial differential operator with respect to x and t, respectively. The model is well known as one of classical models describing shape memory alloys. For the other models, we refer Fremond [4], Fremond-Miyazaki [5] and reference therein. The Falk model (Eq. (2)) has been studied actively in the mathematical literature. In the isothermal case, well-posedness, stability of solitary-wave solution, existence of steady state, travelling wave solution and invariant measure have been investigated by Fang-Grillakis [6], Falk-Laedke-Spatschek [7], Friedman-Sprekels [8], Garcke [9] and Tsutsumi-Yoshikawa [10], respectively. For the full system (Eq. (1)), the well-posedness results are found in, for example, [1113] and so on, and numerical results are found in Hoffmann-Zou [14], Niezgodka-Sprekels [15] by finite element method and in Matus-Melnik-Wang-Rybak [16] and Yoshikawa [17, 18] by the finite difference method. In particular, in Ref. [19], the stability of steady state in the Lyapunov sense was shown. More precisely, the stationary state of Eq. (2) is expressed as a nonlocal nonlinear elliptic problem. If there exists a linearized stable critical point for the functional corresponding to the elliptic problem, then for each neighbourhood U of the equilibrium, we can find a neighbourhood W of the equilibrium such that the solution of Eq. (2) with the initial data in W stays in U for any time. The proof can be shown by the dual variation principle which appears in most of the models in non-equilibrium statistical thermodynamics (see [20]).

The existence of several non-trivial steady states for low-temperature phase and low-energy case is proved in Ref. [21]. The numerical simulation given in Ref. [17] exactly indicates the properties mentioned above. The dual variation structure appears also in a multi-dimensional case [22]; however, well-posedness of the multi-dimensional model corresponding to Eq. (2) is still open in large initial data case due to the propagation of singularity. That is one of our motivations of this problem.

We mention mathematical studies on the motion of curves governed by geometric evolution equations. One of the typical objects is curve-shortening flow derived as an L2 gradient flow for the length functional of curve γ:


where s denotes the arc length parameter of γ. By Gage [23], Gage and Hamilton [24] and Grayson [25], it is well known that the curve-shortening flow shrinks simple closed curves to a point in a finite time. Since the curve-shortening flow can be regarded as a one-dimensional case of mean curvature flow for surfaces, the flow is applicable to various mathematical analysis. For example, the curve-shortening flow plays an important role in studies on phase transition. We also mention the curve-straightening flow which has been attracted a great interest and studied actively in mathematical literature. The flow is derived as an L2 gradient flow for the elastic energy

K(γ):=12γκ2 ds

where κ denotes the scalar curvature of γ. It is well known that the flow is applicable to studies on elastic curve inspired by Bernoulli and Euler. Indeed, the curve-straightening flow under the length constraint L(γ) ≡ C converges to a classical elastic curve so-called elastica. There is also an interest in the study on motion of curves governed by the L2-gradient flow for E under the inextensible condition. Under the condition, the length constraint L(γ) ≡ C is also satisfied for the condition means that the curve does not stretch. As we will state in Section 2.2, the constraint is imposed on each point of curves. Thus, a standard Lagrange multiplier theory does not work. Therefore, we have to make use of geometric properties of curves governed by the flow.

The purpose of this chapter is to derive a mathematical model describing thermoelastic deformation of shape memory wire in R3. In particular, we regard the wire as a closed space curve satisfying the inextensible condition. From the physical point of view, it may be unnatural that the wire does not stretch. However, the contribution of this chapter is to adopt a geometric analysis into a classical thermoelastic theory with phase transition inspired by Falk.


2. Setting and derivation of equations

We denote the closed curve representing shape of wire by Γ = {γ(ξ): ξ ∈ Ξ}, where the variable ξ is an arbitrary parameter not necessarily the arc length parameter. Let us define a displacement vector from a point ξ in an original shape Γ0 = {γ0(ξ): ξ ∈ Ξ} by u(ξ) (see Figure 1); namely, it holds that

Figure 1.

Original shape Γ0 and deformed curve Γ.



  • γ(ξ) := (γ1(ξ), γ2(ξ), γ3(ξ)): vector representing the shape,

  • γ0(ξ):=(γ10(ξ),γ20(ξ),γ30(ξ)): vector representing the original shape,

  • u(ξ):=(u1(ξ),u2(ξ),u3(ξ)): displacement vector.

Throughout this chapter, we denote by L the length of Γ0, and hence, the length of Γ is also L from the non-stretching assumption. To apply the idea by Falk, we need to determine the form of strain and free energy (Eq. (1)) suitable for this setting.

2.1. Definition of strain

We first consider the strain. Let γ0(ξ) be a space closed curve, where ξ is a parameter (not necessary to be the arc length parameter). For γ0 (ξ), we define the displacement vector by u(ξ), and we denote γ(ξ)=γ0(ξ)+u(ξ). Since the relation “strain ≈ line element” holds, let us first pursue line element between γ0 (ξ) and γ (ξ). From the direct calculation, we have


Here, if we assume a smallness of deformation, then we may assume


From now on, we regard the strain as


2.2. Definition of energy functional

Let γ0(ξ) be an initial closed curve and γ(ξ, t) denote a family of closed curves starting from γ0(ξ). Recalling that the arc length parameter s(ξ, t)of γ(ξ, t) is given by


we can write kinetic energy M(γ) of γ as


In a similar manner, thermal energy is defined by


and the curvature energy is expressed as


Observe that the scalar curvature κ is written as

κ={|γξ| 2γξ2(γξ2γξ2)γξ }|γξ|3.

Lastly, the nonlinear elastic energy density is given by

f(ξ γ, θ;ξ γ0) :=16(γ0ξuξ)614(γ0ξuξ)4+12(θθc)(γ0ξuξ)2

and then the nonlinear elastic energy is written as

F(ξγ,θ;ξγ0):=f(ξγ, θ;ξ γ0)|γξ|dξ

where u = γγ0. Thus, we obtain the Helmholtz energy for our setting as


From now on, let sR/LZ=: SL1 be the arc length parameter of the initial closed curve γ0 = γ0(s). It follows from the property of arc length parameter that |γ0′(s)| ≡ 1. In a similar fashion to the above equation, γ(s, t) means the closed curve deformed along evolution from γ0(s). Moreover, in what follows, we assume that γ(s, t) satisfies


which means “s is arc length parameter of γ not only the initial time but also every time t”. From the assumption, we can rewrite M, F0 and F shortly as

M(γ)=0L| tγ |2 ds,F0(θ)=0Lf0(θ)ds,F(s γ,θ;s γ0)=0Lf(sγ,θ;sγ0) ds.

Moreover, since sγs2γ= 0, from Eq. (3), it holds that

K(γ)=0L| s2γ|2ds.

Therefore, the Helmholtz free energy density is denoted by

H(t):=12||tγ(,t)||L2(SL1)2+12||s2γ(,t)||L2(SL1)2+F(sγ(,t), θ(,t);sγ0())+F0(θ(,t)).

We will explain that for the free energy under some assumptions, the following system of nonlinear partial differential equations is derived:

{t2γ+s4γ+sf,sγ(sγ, θ;sγ0)s {(v2 |s2γ|2) sγ}=0, s2v + |s2γ|2v=2 |s2γ|4|s3γ|2+|stγ|2+s2 f,sγ(sγ, θ;sγ0)sγ,tθ s2θ=θ(sγ0s(γ γ0))(tsγsγ0)E4


f,sγ(sγ, θ;s γ0)=(fsγ1, fsγ2, fsγ3)  ={( sγ0su )5( sγ0su )3+(θθc)sγ0su }sγ0.

2.3. Equation of motion

By using the Hamilton principle, we derive an equation of the motion of γ. Namely, we will derive the Euler-Lagrange equation for the functional:

H˜(γ, θ;γ0)=t2t1{12||tγ(,t)||L2(SL1)2+12||s2γ(,t)||L2(SL1)2F(sγ(,t), θ(,t);sγ0())F0(θ(,t))}dt.E5

Let us denote the variation of γ by

γ(s,t;θ):=γ(s,t)+ε ϕ(s,t)

where ε is a sufficiently small positive parameter and ϕ is sufficiently smooth and satisfies ϕ(s, t1) = ϕ(s, t2) = 0. Moreover, from the assumption Eq. (3), it is also necessary to hold that




ϕ has to satisfy


for any sSL1 and t>0. Calculating the first variation of the energy functional, we have

ddεH˜(γ, θ;γ0)|ε=0=t1t2{tγ,tϕs2γ,s2ϕf,sγ(sγ,θ;sγ0),sϕ}dt.

From the integral by parts, the right-hand side is rewritten as follows:


Then the integral Eq. (6) is equal to 0 for any ϕ satisfying ϕ(s,t1)=ϕ(s,t2)=0 and sγsϕ0. For the purpose, we define

V:={ϕ | sγsϕ0}.

The orthogonal complement V of the space V with respect to L2(SL1) inner product is given by

V={s(wsγ)w=w(s,t) is a scalar function}.E7

Here, we remark that in the case where γ is a curve embedded in three-dimensional space (not a planar curve), s2γ0 has to be satisfied for every (s,t)SL1×R+. In the end of this section, we will show the reason why V is given as above. Consequently, if for the direction t2γ+s4γsf,sγ(sγ,θ;sγ0), there exists a scalar function w = w(s, t) such that


then Eq. (6) is equal to 0.

Next, we derive the equation for w. From the assumption Eq. (3), we see that


It follows from Eq. (8) that

{s5γ+s2f,sγ(sγ,θ;sγ0)+s2(wsγ)}sγ=|stγ|2. E9

Differentiating Eq. (3), we obtain

sγs2γ=0,sγs3γ=|s2γ|2,sγs4γ=32s(|s2γ|2),sγs5γ= 2s2(|s2γ|2)+ |s3γ|2.

By the relations, we will rewrite Eq. (9). It follows from the direct calculation that

s2 (wsγ)sγ=s2w  w|s2γ|2.

Since from the definition

f,sγ(sγ, θ;sγ0)={( sγ0s(γγ0) )5( sγ0s(γγ0) )3+(θθc)(sγ0s(γγ0))}sγ0

we also obtain

s2f,sγ(sγ, θ;sγ0)sγ=s2[{(sγ0s(γγ0))5(sγ0s(γγ0))3+(θθc)(sγ0s(γγ0))}sγ0]sγ.

Therefore, substituting these into Eq. (9), we find

0=2s2(|s2γ|2)|s3γ|2+s2ww|s2γ|2+|stγ|2+s2[{(sγ0s(γγ0))5( sγ0s (γγ0))3+(θθc)(s γ0s (γγ0))}sγ0]sγ.

Here setting the new unknown v by v:=w+2|s2γ|2, we can rewrite the equation as follows:

s2v+|s2γ|2v=2 |s2γ|4|s3γ|2+|stγ|2+s2[{(sγ0s(γγ0))5(sγ0s(γγ0))3+(θθc)(sγ0s(γγ0))}sγ0]sγ.

Consequently, under given temperature θ, the equation of motion is given by

{t2γ+s4γsf,sγ(sγ,θ;sγ0)s{(v2|s2γ|2)sγ}=0, s2v+|s2γ|2v= 2|s2γ|4|s3γ|2+|stγ|2s2f,sγ(sγ, θ;sγ0)sγ.

At the rest of this section, we prove that the orthogonal complement of V is given by Eq. (7).

Lemma 1. Let γ(s,t) be a smooth curve in R3 and sSL1 be an arc length parameter of γ for any t. Suppose that s2γ0 holds for all (s,t)SL1×R+ then the orthogonal complement V of the space V={ϕ|sγsϕ0} with respect to L2(SL1) inner product is represented by

V={s(wsγ)w=w(s,t) is a scalar function}.

Proof. Observe that sγ, s2γ and the outer product sγ×s2γ are orthogonal each other. Under the assumption s2γ0, a coordinate system defined on γ consists of the vectors. Then arbitrary vector η = η(s, t) can be represented as

η(s,t)=η1(s,t)sγ(s,t)+η2(s,t) s2γ(s,t)+η3(s,t) sγ(s,t)×s2γ(s,t).

If we assume additionally ηV, then we obtain

0=sηsγ={sη1sγ+η1s2γ+sη2s2γ+η2s3γ+sη3sγ×s2γ+η3 (s2γ×s2γ+sγ×s3γ)}sγ=sη1+η2sγs3γ=sη1η2|s2γ|2

that is

sη1=|s2γ|2η2. E10

In other words, an element of V consists of η1, η2 satisfying Eq. (10) and arbitrary η3. However, we remark that we cannot take η2 freely. Indeed, in order to verify L periodicity of η1, η2, we have to show the following condition:


If ζ(s,t)=ζ1(s,t)sγ(s,t)+ζ2(s,t)s2γ(s,t)+ζ3(s,t)sγ(s,t)×s2γ(s,t)satisfies ⟨ζ, η⟩ = 0, then we see that η1 and η2 satisfy Eqs. (10) and (11) and any η3 satisfies

0L{ζ1η1+ζ2η2|s2γ|2+ζ3η3(sγ×s2γ)2}ds=0. E12

In particular, if we assume η2 ≡ 0 and η3 ≡ 0, then we infer from Eq. (10) that η1C holds true. Then, we deduce from Eq. (12) that


Now we define


Then, ϕ has the period L and satisfies sϕ=ζ1. Substituting it into Eq. (12) and using Eq. (10), we have

0=0L{sϕ η1+ζ2η2|s2γ|2+ζ3η3(sγ×s2γ)2}ds=0L{ϕsη1+ζ2η2|s2γ|2+ζ3η3(sγ×s2γ)2}ds=0L{ϕ |s2γ|2η2+ζ2η2|s2γ|2+ζ3η3(sγ×s2γ)2}ds=0L{(ϕ+ζ2)η2|s2γ|2+ζ3η3(sγ×s2γ)2}ds.

Recalling Eq. (11), we see that the vector-valued function (η2, η3) is orthogonal with (|s2γ|2,0) in the sense of L2 inner product. Therefore, there exists some function μ = μ(t) depending only on t such that

( {ϕ+ζ2} |s2γ|2, (sγ×s2γ)2ζ3)=μ(|s2γ|2,0)

that is,

ϕ+ζ2=μ, ζ30.E13



the function w(s, t) is the L periodic function and satisfies sw=ζ1. It follows from Eq. (13) that


Then ζ(s) is orthogonal with elements of V with respect to L2 inner product. Thus, we get


which completes the proof.


Remark 1. The assumption s2γ0 in Lemma 1 means that the curvature is always non-zero. If s2γ=0 at some point, we cannot determine the tangential vector ∂sγ at the point uniquely. Therefore, we need the assumption in order to give a coordinate system at every point of Γ. On the other hand, in the case of a planar curve, we do not need the assumption. Indeed, by rotating the tangential vector, we can construct a coordinate system.

Remark 2. We mention the elastic flow with the inextensible condition (Eq. (3)), more precisely, L2 gradient flow for K(γ) under the constraint (Eq. (3)). To the best of our knowledge, the problem was first considered by N. Koiso [26] for planar closed curves. With the aid of smoothing effect of the elastic energy E, the Cauchy problem on the elastic flow has a unique classical solution and the solution converges to an equilibrium state as t → ∞ in the C-topology. The result can be extended to the following case: (i) L2 gradient flow for E under the area-preserving condition and (C) [27] and (ii) L2gradient flow for Tadjbakhsh-Odeh energy functional under the constraint (C) [28]. Moreover, the result [26] was also extended to the case of space curves [29].

2.4. Derivation of heat equation

In this subsection, we study the energy law. We confirm thermal energy conservation law (the first law of thermodynamics) and the increasing law of entropy (the second law of thermodynamics). To begin with, we consider the first law of thermodynamics. According to Ref. [30], thermal energy conservation law for thermoelastic system is given by


where S, q and h are entropy, thermal velocity and external heat, respectively. In our setting, thermal transfer occurs only on wire, and the wire does not expand. Then, we may regard ∇ ⋅ q as ∂sq as the same as one-dimensional case, where s is necessary to be arc length parameter. By the same reason, the Fourier law q = ∇θ is replaced by


The Helmholtz free energy density

f˜=f˜(s2γ, sγ, θ ;γ0)=12|s2γ |2+f(sγ, θ ;γ0)+f0(θ)

and the entropy S are connected with the relation

S=f˜ θ.

Then, the conservation law, Eq. (14), is rewritten as follows:

 θ f0''(θ) tθ s2θ=θ(sγ0s(γγ0))(tsγsγ0)+h.E16

We note that the Clausius-Duhem inequality holds automatically:

tS+s(qθ)hθ .

Indeed, we observe from Eqs. (14) and (15) that

tS+s(qθ)=hsqθ+s(qθ)=hθqsθθ2=hθ+|sθθ|2hθ .

The Clausius-Duhem inequality corresponds to the second law of thermodynamics. For more precise information of the inequality, we refer to, for example, 1.11 of chapter 4 in Ref. [2].

Here, we assume external heat source h = 0 and adopt the well-known form:


Then since f0(θ)=1/θ, Eq. (16) is reduced to


We thus obtain the system of equation as Eq. (4).


3. Dual variation structure

Let us rewrite Eq. (4):

t2γ+s4γ +sf,sγ(sγ, θ;sγ0)s{(v2 |s2γ|2)sγ}=0,E17
 s2v+|s2γ|2v=2|s2γ|4|s3γ|2+|stγ|2+s2 f,sγ(sγ, θ;sγ0)sγ,E18
tθs2θ=θ(sγ0s(γγ0))(tsγsγ0), (t,s)(0,T)×SL1,E19
γ(0,s)=γ0(s), tγ(0,s)=γ1, θ(0,s)=θ0, sSL1.E20

In this section, we show that the problem in Eqs. (17)–(20) has also dual variation structure as well as the problem (Eq. (2)). The structure plays an important role to prove the dynamical stability of infinitesimally stable stationary state.

We assume that the system has sufficiently smooth solution (γ, θ) satisfying θ > 0. Then initial data also has to satisfy

|sγ0|=1, sγ0sγ1=0.


f1(sγ) :=12(γ0s(γγ0)s)2,f2(sγ) :=16(γ0s(γγ0)s)614(γ0s(γγ0)s)4θc2(γ0s(γγ0)s)2,

we have the relation f=θf1+f2. Multiplying Eq. (17) by ∂tγ and integrating it with respect to s, we obtain

ddt(12||tγ||L22+12||s2γ||L22)=f,sγ, stγ=ddt0Lf2(sγ)ds0Lθtf1(sγ)ds.

Integrating Eq. (19), we find


Then for the quantity


it holds that


Moreover, for the quantity


we deduce from Eq. (19) that


Finally, we confirm an important structure of the system (17)–(20). We denote the stationary state of θ by θ¯>0, and the corresponding equilibrium by γ satisfy the following constraint:

bE(γ0, γ1, θ0)=E(γ, 0, θ¯),s4γ=s{θ¯f1,sγ+f2,sγ}.

Eliminating θ¯ by the relation

b=Lθ¯+12||s2γ||L22+0Lf2(sγ) ds

the stationary state of this problem satisfies the following nonlinear nonlocal problem:

s4γ=s{1L(b12||s2γ||L220Lf2(sγ)ds)f1,sγ(sγ)+f2,sγ(sγ)}. E21

Eq. (21) is derived as the Euler-Lagrange equation of the functional


where y= sγH2(SL1,S2) and


We remark that the following relation between Jb and W holds true:


The relation is called semi-unfolding minimality. Thus, if (γ, θ) is a non-stationary state, b=E(γ0,γ1,θ0), y¯=sγ¯ is a linearized stable critical point of Jb=Jb(y), y H1(SL1, S2) and θ¯> 0 is a constant satisfying E(γ¯,0,θ¯) = b, then it holds that for y=sγ

Jb(y)Jb(y¯)W(sγ0, θ0)W(y¯,θ¯).

By this structure, we can infer that any infinitesimally stable stationary state is dynamically stable, that is, stable in the Lyapunov sense. A critical point y¯=sγ¯ of Jb for γ¯H2(SL1,S2) is infinitesimally stable if there exists ε0 > 0 such that any ε1(0,ε02] admits δ0 > 0 such that if ||s(γγ¯)||H1<ε0 and Jb(sγ) Jb(sγ¯)<δ0 then


The definition of infinitesimally stable is obviously weaker than the one of well-known linearized stable which means that for a critical point y¯=sγ¯ of Jb, the quadratic form


is a positive definite for any wH2(SL1,S2).

Theorem 1. Assume that θ¯>0 is a constant and that γ¯ is an infinitesimally stable critical point of Jb with constraint sγH2(SL1,S2). Then (γ¯,θ¯) is a dynamically stable in the sense that for any ε > 0, there exists δ > 0 such that if

E(γ0, γ1, θ0)= b,||s(γ0γ¯)||H1<δ, |1L0Llogθ0(s)dslogθ¯|<δE22


supt0||s(γ(t)γ¯)||H1<ε, |1L0Llogθ(t,s)dslogθ¯|<ε.

Proof. We first show the semi-unfolding minimality. From the energy conservation law, we see that

b =12||tγ||L22+12||s2γ||L22+0Lθds+0Lf2(sγ)ds.

It follows from the Jensen inequality that


Then we have

W(sγ, θ)0Lf1(sγ)dslog(1L0Lθds)0Lf1(sγ)dsLlog1L(b12||sy||L220Lf2(y)ds)=LJb(sγ).

We have thus completed to show the semi-unfolding minimality. Recall that γ¯H2(SL1,S2) is an infinitesimally stable critical point of Jb(sγ). Thus, we find ε0 > 0such that any ε1(0,ε02] admits δ0 > 0 such that if ||s(γγ¯)||H1<ε0 and Jb(sγ) Jb(sγ¯)<δ0 then


From the above properties, it holds that

Jb(sγ(t))1LW(sγ(t), θ(t))1LW(sγ0, θ0).E24

Moreover, for the constant θ¯>0, we obtain

W(sγ¯, θ¯)=L[1L0Lf1(sγ)ds log1L(b12||sy¯||L220Lf2(y¯)ds)]=LJb(s γ¯),


Jb(sγ¯)=1LW(sγ¯,θ¯). E25

Given ε > 0, setting δ(0,ε02] and satisfying Eq. (22), we have

1L|W(sγ0, θ0)W(sγ¯,θ¯)|1L(f1, sγ(sγ0)L+f1, sγ(sγ¯)L)s(γ0γ¯)L1+|1L0Llogθ0dslogθ¯|<min(δ0,ε).E26

Therefore, it follows from Eq. (24) to Eq. (26) that

Jb(sγ(t))Jb(sγ¯)=1L(W(sγ0, θ0)W(sγ¯,θ¯))<δ0.

If ||s(γ(t)γ¯)||H1=δ(ε0/2 <ε0), then we apply Eq. (23) for ε1 = δ, and hence


which is a contradiction. Thus, we have


Here, from γC([0,);H2) and s(γ0γ¯)H1<δ, it follows that


for any t ≥ 0.

From the semi-unfolding minimality (Eqs. (24) and (25)) and the linearized stability of Jb, we observe that


Then, combining Eq. (26) with Eq. (27), we have


where δ is small enough such that δ < ε / C. This completes the proof.


Remark 3. Both the existence of solution for evolution equations (Eq. (4)) and non-trivial solutions for stationary problem (Eq. (21)) are open problems. In the straight material case (i.e. the problem (2)), smooth solution for Eq. (2) is assured in [11] (we also refer to chapter 5 in [2]). The existence results of non-trivial solution for stationary problem (Eq. (21)) in low-temperature and low-energy cases can be found in [21, 31].


4. Concluding remarks

In this chapter, we propose the new mathematical model describing the movement of wire made of shape memory alloys. The derived system of nonlinear partial differential equations is a thermoelastic system with phase transition and non-stretching constraint. The Falk model (Eq. (2)) represents the dynamics for crystal as a stack of layers, whose displacement is restricted to move only on one direction. On the other hand, our model describes the dynamics of wire. We emphasize that our model allows the displacement of each direction. Thus, our model may describe a more realistic motion of wire made of shape memory alloys. Moreover, it is also interesting to regard our model as a mathematical problem on elastic curve with heat conduction. To the best of our knowledge, there is no result considering such a mathematical problem.

We mention the mathematical contribution of the present chapter. We prove the dynamical stability of an infinitesimally stable stationary state by finding the dual variation structure in our model. This property shall be applicable, for example, to assure the strength of not only a wire in an original shape but also a deformed wire. Indeed, in the straight material case, namely in the Falk model (Eq. (2)), numerical simulation shows the stability in this sense (see [17]).



This work was partially supported by JSPS KAKENHI Grant-in-Aid for Scientific Research (A) and (C), Grant Numbers 26247013 and 16K05234.


  1. 1. Falk F. Elastic phase transitions and nonconvex energy functions. In: Pitman Research Notes in Mathematics Series. Vol 185. IRSEE. Harlow: Longman Science and Technology; 1990. pp. 45–59
  2. 2. Brokate M, Sprekels J. Hysteresis and phase transitions (Applied Mathematical Sciences). 121st ed. New York: Springer-Verlag; 1996
  3. 3. Jiang S, Racke R. Evolution Equations in Thermoelasticity. Monographs and Surveys in Pure and Applied Mathematics. 112nd ed. Boca Raton: Chapman & Hall/CRC; 2000
  4. 4. Fremond M. Non-Smooth Thermomechanics. Berlin, Heidelberg: Springer-Verlag; 2002. p. 480. DOI: 10.1007/978-3-662-04800-9
  5. 5. Fremond M, Miyazaki S. Shape Memory Alloys. Wien: Springer-Verlag; 1996. p. 147. DOI: 10.1007/978-3-7091-4348-3
  6. 6. Fang Y-F, Grillakis MG. Existence and uniqueness for Boussinesq type equations on a circle. Communications in Partial Differential Equations. 1996;21:1253–1277
  7. 7. Falk F, Laedke EW, Spatschek KH. Stability of solitary-wave pulses in shape-memory alloys. Physical Review B. 1987;36:3031–3041
  8. 8. Friedman A, Sprekels J. Steady states of austenitic-martensitic domains in the Ginzburg-Landau theory of shape memory alloys. Continuum Mechanics and Thermodynamics 1990;2:199–213
  9. 9. Garcke H. Travelling wave solutions as dynamic phase transitions in shape memory alloys. Journal of Differential Equations. 1995;121:203–231
  10. 10. Tsutsumi Y, Yoshikawa S. Invariant measures for the isothermal Falk model of shape memory alloys. GAKUTO International Series. Mathematical Sciences and Applications. 2015;37:163–182
  11. 11. Sprekels J, Zheng S. Global solutions to the equations of a Ginzburg-Landau theory for structural phase transitions in shape memory alloys. Journal of Physics D. 1989;39:59–76
  12. 12. Hoffmann K-H, Zochowski A. Analysis of the thermoelastic model of a plate with nonlinear shape memory alloy reinforcements. Mathematical Methods in the Applied Sciences. 1992;15:631–645
  13. 13. Yoshikawa S. Weak solutions for the Falk model system of shape memory alloys in energy class. Mathematical Methods in the Applied Sciences. 2005;28:1423–1443
  14. 14. Hoffmann K-H, Zou J. Finite element approximations of Landau-Ginzburg’s equation model for structural phase transitions in shape memory alloys. RAIRO Mathematical Modelling and Numerical Analysis. 1995;29:629–655
  15. 15. Niezgodka M, Sprekels J. Convergent numerical approximations of the thermomechanical phase transitions in shape memory alloys. Numerische Mathematik. 1991;58:759–778
  16. 16. Matus P, Melnik RVN, Wang L, Rybak I. Applications of fully conservative schemes in nonlinear thermoelasticity: Modelling shape memory materials. Mathematics and Computers in Simulation. 2004;65:489–509
  17. 17. Yoshikawa S. A conservative finite difference scheme for the Falk model system of shape memory alloys. ZAMM Zeitschrift für Angewandte Mathematik und Mechanik. 2015;95(12):1393–1410
  18. 18. Yoshikawa S. An error estimate for structure-preserving finite difference scheme for the Falk model system of shape memory alloys. IMA Journal of Numerical Analysis. 2017;37:477–504
  19. 19. Suzuki T, Yoshikawa S. Stability of the steady state for the Falk model system of shape memory alloys. Mathematical Methods in the Applied Sciences. 2007;30:2233–2245
  20. 20. Suzuki T. Mean Field Theories and Dual Variation—Mathematical Structures of the Mesoscopic Model. 2nd ed. Amsterdam, Paris: Atlantis Press; 2016
  21. 21. Suzuki T, Tasaki S. Stationary solution to the Falk system on shape memory alloys. Mathematical Methods in the Applied Sciences. 2010;33:994–1011
  22. 22. Suzuki T, Yoshikawa S. Stability of the steady state for multi-dimensional thermoelastic systems of shape memory alloys. Discrete and Continuous Dynamical Systems—Series S. 2012;5:209–217
  23. 23. Gage M. Curve shortening makes convex curves circular. Inventiones Mathematicae. 1984;76:357–364
  24. 24. Gage M, Hamilton R. The shrinking of convex plane curves by the heat equation. Journal of Differential Geometry. 1986;23:69–96
  25. 25. Grayson MA. The heat equation shrinks embedded plane curves to round points. Journal of Differential Geometry. 1987;26:285–314
  26. 26. Koiso N. Motion of a springy circle wire. Science Reports—College of General Education, Osaka University. 1991;40:15–19
  27. 27. Okabe S. The motion of elastic planar closed curves under the area-preserving condition. Indiana University Mathematics Journal. 2007;56:1871–1912
  28. 28. Okabe S. The dynamics of elastic closed curves under uniform high pressure. Calculus of Variations and Partial Differential Equations. 2008;33:493–521
  29. 29. Koiso N. On motion of a curve towards elastica. In: Actes da la Table Ronde de Geometrie Differentielle, Semin Congr. Vol. 1, 1992; Luminy. Paris: Soc. Math. France; 1996. pp. 403–436
  30. 30. Lagnese J, Lions J-L. Modelling Analysis and Control of Thin Plates. Recherches en Mathematiques Appliquees (Research in Applied Mathematics). 6th ed. Paris, Masson; 1988
  31. 31. Suzuki T, Tasaki S. Stationary solutions to a thermoelastic system on shape memory materials. Nonlinearity. 2010;23:2023–2656

Written By

Shinya Okabe, Takashi Suzuki and Shuji Yoshikawa

Submitted: 23 November 2016 Reviewed: 12 April 2017 Published: 20 September 2017