# Bifurcation and Chaos in Integer and Fractional Order Two-Degree-of-Freedom Shape Memory Alloy Oscillators

Authors: Karthikeyan Rajagopal; Anitha Karthikeyan; Prakash Duraisamy; Riessom Weldegiorgis
Journal: Complexity
 (2018)

---

## Abstract

A two-degree-of-freedom shape memory oscillator derived using polynomial constitutive model is investigated. Periodic, quasiperiodic, chaotic, and hyperchaotic oscillations are shown by the shape memory alloy based oscillator for selected values of the operating temperatures and excitation parameters. Bifurcation plots are derived to investigate the system behavior with change in parameters. A fractional order model of the shape memory oscillator is presented and dynamical behavior of the system with fractional orders and parameters are investigated.


---

## Body

## 1. Introduction

A shape memory alloy (SMA) shows temperature and stress induced martensitic phase transformation, drastically changing the material’s mechanical properties. In other words, after an apparent plastic deformation, it will return to their original shape when heated. The same materials, in a certain temperature range, can be strained up to approx. 10% and still will return to their original shape when unloaded. These unusual effects are called thermal shape memory and superelasticity (elastic shape memory), respectively, which can use both sensor and actuator [1]; it is increasing range of potential applications, such as Biomedical [2], Robotics [3], and Aerospace. Stress-strain curve study reveals hysteretic phenomenon of SMAs. Hysteretic effect leads to nonlinearity. This property results in high damping capacity [2] which can be utilized for amplitude attenuation in vibrating system [4] so the application of these materials can be extended to control elements too.Various constitutive models of shape memory oscillators (SMO) have been developed since 1990, in order to analyze the dynamical behavior [5–9]. The study reveals that SMOs can expose to chaotic region when it is treated as single degree-of-freedom system but can extend to hyperchaos when it is treated as two-degree-of-freedom system [10, 11]. Similar property is noted down in other SMAs. According to Falk model strain and temperature are the only state variables; free energy is validated for various temperatures and shows unstable points and is not considering twinned martensite. Brinson’s model describes the reorientation process clearly compared to other models, the numerical simulations reveals the internal subloops present in thermomechanical [12] loading. Tanaka model considered martensitic volumetric fraction so it can describe compressive behavior accurately [13]. The transient response and attractor’s multistability are interesting property in nonlinear study, and the numerical analysis of SMOs nonlinear models shows these properties significantly [10]. Dynamical jumps of SMOs were studied and found that the response is not only depending on the forcing amplitude but also on the way the forcing frequency is modified [14].Experimental dynamic analysis of SMAs was carried out and found nonlinearity in phase transformation leads to complicated dynamic behavior of the system [15]. Nonlinear dynamics of SMAs are investigated experimentally and it is noted that the actuator response shows unpredictability, when low-to-moderate voltage was applied to the system and added that could affect the system design and can cause control difficulties in precise applications [3].The real systems are fractional in nature, so there is no surprise on treatment of fractional order describing behavior of the system more accurate than the integer methods. Because of holding nonlocal property, that is, the next state of a system depends not only upon its current state but also upon all of its historical states, fractional calculus can provide more realistic results. Recently fractional order calculus used to refine the results in various fields like thermodynamics, mechatronics systems, chaos theory, and biomedical system as well. When a chaotic nonlinear system is treated as integer order, it demands minimum order of 3 for chaos to appear [16]. Conversely while the system is treated as fractional order, chaos can be identified with lesser order. For example, Chau’s circuit with order 2.7 can produce chaotic attractor [17]. Recently the active and passive vibration damping capacity of SMAs are studied. Oberaigner model has been taken, and fractional order results reveal the chaotic behaviors precisely [18]. The generalized theory of electrothermoelasticity of fractional order heat transfer describes the behavior of the particles of an elastic body more realistically than the theory of generalized thermoelasticity with integer order [19]. Yu et al. developed the theory of fractional order generalized electro-magneto-thermoelasticity for anisotropic and linearly electro-magneto-thermo-elastic media and extended that the fractional order has great effect on the response when the material is imposed a sudden heating [20].


## 2. Dynamics of the Two-Degree-of-Freedom Shape Memory Alloy (2DOF SMA) Oscillator

Even though the authors of [11] have done a [0-1] test on different models of SMAs, the complete dynamical analysis of the system is less investigated. Hence we are interested in analyzing the two-degree-of-freedom (2DOF) SMA [21] shown in Figure 1 because of its higher dimension and complex behaviors.Figure 1
The two-degree-of-freedom oscillator with shape memory alloys.The dimensionless mathematical model of the 2DOF SMA polynomial constitutive model can be derived using the equation of motions and is given by [11, 21](1)x˙=y,y˙=δ1sin⁡ω1t-ε1+ε2v21μy+ε2v21μz-θ1-1+v212μθ2-1x+v212μθ2-1z+ϕ1x3-λ1x5-ϕ2v212μz-x3+λ2v212μz-x5,z˙=w,w˙=α212δ2sin⁡ω2t+ε2v21y-ε2α21+ε3v21α32w+v212θ2-1x-v212θ2-1+v212v322θ3-1z+ϕ2v212z-x3-λ2v212z-x5+ϕ3v212v322z3-λ3v212v322z5,where δ1,δ2,ω1,ω2 are the excitation parameters, ε1,ε2,ε3 are the dissipation parameters, μ is the mass relation parameter, and v21,v32,α21,α32,ϕ1,ϕ2,ϕ3,λ1,λ2,λ3 are the SMA properties. The states x and z are the dimensionless displacements, and y and w are the dimensionless velocity components of the 2DOF SMA. For the parameter values,(2)ω1=ω2=1,ε1=ε2=ε3=0.2,ϕ1=ϕ2=ϕ3=1.3e+3,λ1=λ2=λ3=4.7e+5,v21=v32=1,μ=1,δ1=0.06,δ2=0,α21=α32=1,and depending on the temperatures of the SMA given by the parameters θ1,θ2,θ3, the system changes the character as periodic, quasiperiodic, chaotic, and hyperchaotic [11]. Table 1 shows the different types of behavior of the SMA with various values of θ1,θ2,θ3 with the respective finite time Lyapunov exponents (LEs) calculated using Wolf algorithm [22] for 20000s.Table 1
Types of SMA behavior.
Parameters
Type of SMA behavior
Lyapunov exponents
θ
1
=
0.7
,
θ
2
=
3.5
,
θ
3
=
0.7
Periodic
-0.128;-0.140;-0.257;-0.275
θ
1
=
1.5
,
θ
2
=
3.5
,
θ
3
=
1.5
Quasiperiodic
0;-0.209;-0.265;-1.272
θ
1
=
1.5
,
θ
2
=
0.5
,
θ
3
=
1.5
Chaotic
0.378;0;-0.438;-0.786
θ
1
=
1.5
,
θ
2
=
1.5
,
θ
3
=
1.5
Hyperchaotic
0.453;0.233;-0.651;-0.835Figure2 shows the 2D phase portraits of the SMA system with initial conditions [0,0,0,1] for various values of θ1,θ2,θ3 as given in Table 1. It can be observed from Figure 2 that the SMA system shows periodic, quasiperiodic, chaotic, and hyperchaotic behaviors.Figure 2
2D phase portraits of the 2DOF SMA system: (a) periodic; (b) quasiperiodic; (c) chaotic; (d) hyperchaotic.
(a)
(b)
(c)
(d)Bifurcation. Bifurcation analysis is helpful in investigating the qualitative changes in the states of a system with respect to the variation in a parameter [23]. Recent techniques have shown that bifurcation analysis can be done with experiments rather than only mathematical methods. Both continuous and discrete models can be analyzed by bifurcations. To analyze the 2DOF SMA using bifurcation plots, we fix all the other parameters like in (2) and consider θ2 as the bifurcation parameter. The initial conditions for the first iteration are taken as [0,0,0,1] and are reinitialized to the end values of state trajectories at each iteration. Figure 3(a) shows the bifurcation diagram of system (1) with respect to the parameter θ2. The parameter values of the SMA system are taken as in (2) with θ1=1.5,θ3=1.5, and initial conditions for the first iteration are taken as [0,0,0,1] and are reinitialized to the end values of the state variables at the end of each iteration. It can be clearly observed from Figure 2(a) that the SMA shows periodic, quasiperiodic, chaotic, and hyperchaotic oscillations. For θ2<0.45 the SMA system shows periodic limit cycles. For 2.9≤θ2≤3.75 the SMA systems show quasiperiodic/torus attractor. Figure 3(b) shows the bifurcation of 2DOF SMA for parameter θ1 with θ2=1.5,θ3=1.5 which shows periodicity, chaos, hyperchaos, and quasiperiodicity. Period 2 limit cycles are seen for 0≤θ1≤0.75, period 3 limit cycles are exhibited for 2.2<θ1≤2.5 and 3<θ1≤3.25, and period 4 limit cycles are seen for 4.3≤θ1≤4.5. Chaotic oscillations are seen when 0.75<θ1≤1.25, 1.8<θ1≤2.2,3.25<θ1≤3.5, and 3.75<θ1≤4.3 and hyperchaos is seen for a small range of 1.25<θ1≤1.8. Quasiperiodic oscillations are shown for 3.5<θ1≤3.75. Figure 3(c) shows the bifurcation of the 2DOFSMA system with δ1 and θ1=1.5,θ2=1.5,θ3=1.5. Chaotic and hyperchaotic regions are seen in the region 0.048<δ1≤0.075 and chaotic region for 0.035<δ1≤0.045. For δ1<0.018 we could observe period 1 limit cycles and quasiperiodic oscillations for 0.018<δ1≤0.03.Figure 3
Bifurcation plots of 2DOF SMA with (a)θ2, (b) θ1, and (c) δ1 with initial conditions [0,0,0,1].
(a)
(b)
(c)


## 3. Fractional Order Two-Degree-of-Freedom Shape Memory Alloy (FO2DOFSMA) Oscillator

Fractional order calculus is as old as integer order as seen from a letter written by Leibniz to L’Hopital [24]. A more common form of general order is given by(3)Dq=dqdtqq<01q=0∫atdτ-qq>0,where q is the order. The three ways for simulating fractional order systems [25] are computational methods based on the mathematical equations, approximation through a rational system in discrete time, and approximation of the fractional system using rational function in continuous time. Most researchers prefer approximation through a rational system in discrete time in which the fractional order system is replaced by its discrete equivalent and approximated again by truncating the polynomial series, which reduces the requirement of infinite memory.Many different approaches for numerical simulations of fractional order systems have been extensively investigated [26–28]; however there is always a tradeoff between computational efficiency, complexity, and the accuracy of the approximations. Recent researchers have been working on developing fast convolution quadrature algorithms relevant to fractional differential equations because fractional calculus operators work either in continuous or discrete convolution of some form.There are three commonly used definitions of the fractional order differential operator: Grünwald-Letnikov, Riemann–Liouville, and Caputo [23, 29, 30]. We used the Grünwald-Letnikov (GL) definition to derive the fractional order model of 2DOFSMA. The Grünwald–Letnikov method is proceeding iteratively but the sum in the scheme becomes longer and longer, which reflects the memory effect. The binomial coefficients are recursively defined and show very smooth properties. The Grünwald–Letnikov derivative is used in many numerical schemes for discretizing fractional diffusion equations from a continuous Riemann–Liouville approach [31–33] and involves a discrete convolution between a “weight” or binomial coefficient function and the function of interest for differentiation. The analytical model of the binomial coefficients is well established in the literature [33].Most of the literatures have used the Caputo method for numerical solutions of the fractional order systems. But GL method has benefits over the other methods of solving fractional orders due to the smoothness of the resultant approximations [34]. Hence we use the GL method to derive the fractional order 2DOF SMA. The GL derivative can be defined as(4)Dtqaft=limh→0⁡1hq∑j=0t-a/h-1jqjft-jh=limh→0⁡1hαΔhqft,where a and t are limits of the fractional order equation, Δhqf(t) is generalized difference, h is the step size, and q is the fractional order of the differential equation.For numerical calculations the above equation is modified as(5)Dtqt-Lft=limh→0⁡h-q∑j=0Ntβjft-jh. We use the short memory principle to limit the required memory for binomial coefficients using(6)N=min⁡th,Lh.The binomial coefficients required for the numerical simulation is calculated as(7)βj=1-a+qjβj-1.Let the general form of the 3D fractional order system be defined as(8)Dqx=Ax,y,z,w,t,Dqy=Bx,y,z,w,t,Dqz=Cx,y,z,w,t,Dqw=Dx,y,z,w,t.In order to simulate system (8) using GL method we use the discretization method discussed [34, 35],(9)xtk=Axtk-1,ytk-1,ztk-1,wtk-1hqx-∑j=1Nβjqxxtk-j,ytk=Bxtk-1,ytk-1,ztk-1,wtk-1hqy-∑j=1Nβjqyytk-j,ztk=Cxtk-1,ytk-1,ztk-1,wtk-1hqz-∑j=1Nβjqzztk-j,wtk=Cxtk-1,ytk-1,ztk-1,wtk-1hqw-∑j=1Nβjqwwtk-j,where β is the binomial coefficients calculated using (11). The value of N is taken as the truncation window size L and as k when all the available memory elements are used.Let us define the FO2DOFSMA oscillator as(10)Dqxx=y,Dqyy=δ1sin⁡ω1t-ε1+ε2v21μy+ε2v21μz-θ1-1+v212μθ2-1x+v212μθ2-1z+ϕ1x3-λ1x5-ϕ2v212μz-x3+λ2v212μz-x5,Dqzz=w,Dqww=α212δ2sin⁡ω2t+ε2v21y-ε2α21+ε3v21α32w+v212θ2-1x-v212θ2-1+v212v322θ3-1z+ϕ2v212z-x3-λ2v212z-x5+ϕ3v212v322z3-λ3v212v322z5.Using (9) in (10), the discrete form of the FO2DOFSMA is given by(11)xtk=ytk-1hqx-∑j=1Nβjqxxtk-j,ytk=δ1sin⁡ω1t-ε1+ε2v21μytk-1+ε2v21μztk-1-θ1-1+v212μθ2-1xtk-1+v212μθ2-1ztk-1+ϕ1xtk-13-λ1xtk-15-ϕ2v212μztk-1-xtk-13+λ2v212μztk-1-xtk-15hqy-∑j=1Nβjqyytk-j,ztk=wtk-1hqz-∑j=1Nβjqzztk-j,wtk=α212δ2sin⁡ω2t+ε2v21ytk-1-ε2α21+ε3v21α32wtk-1+v212θ2-1xtk-1-v212θ2-1+v212v322θ3-1ztk-1+ϕ2v212ztk-1-xtk-13-λ2v212ztk-1-xtk-15+ϕ3v212v322ztk-13-λ3v212v322ztk-15hqw-∑j=1Nβjqwwtk-j.The value of N is taken as the truncation window size L and as k when all the available memory elements are used.For the parameters values given by (2) and Table 2, the FO2DOFSMA oscillator shows different behaviors as like the integer order SMA. For initial conditions [0,0,0,1], commensurate order q=0.99, and step size h=0.001, the different behaviors of the FO2DOFSMA are given in Figure 4.Table 2
Types of fractional order 2DOF SMA behavior.
Parameters
Type of SMA behavior
Lyapunov exponents
θ
1
=
1.5
,
θ
2
=
3.5
,
θ
3
=
1.5
Periodic
-0.128;-0.140;-0.257;-0.275
θ
1
=
1.5
,
θ
2
=
2
,
θ
3
=
1.5
Quasiperiodic
0;-0.209;-0.265;-1.272
θ
1
=
1.5
,
θ
2
=
1.5
,
θ
3
=
1.5
Chaotic
0.378;0;-0.438;-0.786
θ
1
=
1.5
,
θ
2
=
1
,
θ
3
=
1.5
Hyperchaotic
0.453;0.233;-0.651;-0.835Figure 4
2D phase portraits of the FO2DOFSMA oscillator: (a) periodic; (b) quasiperiodic; (c) chaotic; (d) hyperchaotic.
(a)
(b)
(c)
(d)We study the bifurcation of the FO2DOFSMA with fractional orderq and parameter θ2. The bifurcation of the FO2DOFSMA oscillator with the fractional order q is shown in Figure 5(a). The initial conditions are taken as 0,0,0,1 and are reinitialized in end of each iteration to the final values of the state variables. The FO2DOFSMA shows chaotic oscillations for q≥0.968. Figure 5(b) shows the bifurcation of the oscillator with θ2 and shows chaotic oscillations for 0.47<θ2≤0.8, 1.1<θ2≤1.6. Hyperchaotic oscillations are seen for 0.8<θ2≤1.1, quasiperiodic oscillations for 1.6<θ2≤2.2, and period 4 oscillations for 2.2<θ2≤4. For Figure 5(b), we used initial conditions without reinitialization. Figure 6 shows the bifurcation plots of the oscillator with reinitialization of the initial conditions. Figure 6(b) shows zoomed in portion of Figure 6(a) to shows the period doubling route to chaos.Figure 5
Bifurcation plots of 2DOFSMA oscillator with (a) commensurate fractional orderq and (b) parameter θ2 with q=0.99 and θ1=θ2=1.5. The initial conditions for plot (a) are taken as [0,0,0,1] and are reinitialized in end of each iteration to the final values of the state variables and for plot (b) the initial value is fixed at [0,0,0,1].
(a)
(b)Figure 6
(a) Bifurcation of the FO2DOFSMA oscillator with parameterθ2 with q=0.99 and θ1=θ2=1.5. The initial conditions for are taken as [0,0,0,1] and are reinitialized in end of each iteration to the final values of the state variables. (b) shows the zoomed in view of (a) between 4≤θ2≤6.


## 4. Conclusion

A two-degree-of-freedom polynomial constitutive shape memory alloy oscillator is investigated and it is shown that periodic, quasiperiodic, chaotic, and hyperchaotic behaviors are shown by the oscillator for various values of the excitation parameters and temperatures. Bifurcation plots are derived to show the existence of chaotic and periodic oscillations. The Grünwald–Letnikov derivative is used to derive the fractional order model of the oscillator and bifurcation plots for fractional order and parameter are derived to study the dynamical behavior of the oscillator. Compared to the integer order model, the fractional order model shows more complex chaotic behavior and the quasiperiodic region is extended as chaotic regions and the oscillator shows more extended regions of chaos with parameters. This proves that the fractional order model shows more complex behavior and hence designing controllers to suppress chaos in fractional orders is quiet complicated.



---

*Source: 8365845-2018-01-31.xml*