The Scientific World Journal

Volume 2013 (2013), Article ID 865176, 7 pages

http://dx.doi.org/10.1155/2013/865176

Research Article

Hysteresis Modeling of Magnetic Shape Memory Alloy Actuator Based on Krasnosel'skii-Pokrovskii Model

Miaolei Zhou, Shoubin Wang, and Wei Gao

College of Communication Engineering, Jilin University, Changchun 130022, China

Received 6 March 2013; Accepted 9 April 2013

Academic Editors: A. P. Pyatakov, A. Savchuk, and H. Zhou

Copyright © 2013 Miaolei Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

As a new type of intelligent material, magnetically shape memory alloy (MSMA) has a good performance in its applications in the actuator manufacturing. Compared with traditional actuators, MSMA actuator has the advantages as fast response and large deformation; however, the hysteresis nonlinearity of the MSMA actuator restricts its further improving of control precision. In this paper, an improved Krasnosel'skii-Pokrovskii (KP) model is used to establish the hysteresis model of MSMA actuator. To identify the weighting parameters of the KP operators, an improved gradient correction algorithm and a variable step-size recursive least square estimation algorithm are proposed in this paper. In order to demonstrate the validity of the proposed modeling approach, simulation experiments are performed, simulations with improved gradient correction algorithm and variable step-size recursive least square estimation algorithm are studied, respectively. Simulation results of both identification algorithms demonstrate that the proposed modeling approach in this paper can establish an effective and accurate hysteresis model for MSMA actuator, and it provides a foundation for improving the control precision of MSMA actuator.

1. Introduction

Hysteresis phenomena are widely exhibiting between input and output response of smart material actuators; MSMA actuator is a new kind smart material actuator which has more advantages than other smart material actuators: rapid frequency response, high control precision, and high deformation rate [1–3], but hysteresis between its input and output displacement is a barrier which reduces the control precision and restricts the further application of MSMA actuator. Therefore, in order to improve the position control accuracy of MSMA actuator, it is necessary to mitigate the effects of hysteresis. Many approaches have been proposed to mitigate the effects of hysteresis, and up until now, establishing the hysteresis model of smart material actuators is the most popular solution to this problem [4–8].

Preisach-type hysteresis models, such as Preisach model [9–11], Prandtl-Ishlinskii (PI) model [12, 13], and KP model [14, 15] are the most popular ones which have been widely used in modern control techniques to mitigate the effects of unknown hysteresis. In this paper, a Preisach-type hysteresis model—modified KP model is used to describe the hysteresis of MSMA actuator. KP model can be seen as an integral of a weighting function and KP operators; compared with other Preisach-type hysteresis models, KP model has the advantage as it can describe more complex hysteresis and behave well for offline identification. Improved gradient correction algorithm and variable step-size recursive least square estimation algorithm are applied to identify the weighting parameters, and simulation experiments are investigated in this paper to verify the validity of the proposed method.

2. Structure of the Hysteresis Model of MSMA Actuator Based on KP Model

The schematic diagram of hysteresis between input magnetic induction density and output displacement of MSMA actuator is shown in Figure 1. The hysteresis makes it difficult to get high control accuracy in the MSMA actuator-based system, and in this section, the structure of the proposed hysteresis model is described.

Figure 1: Major hysteresis loop and inner hysteresis loops of MSMA actuator.

The modified KP operator is shown in Figure 2; different from the traditional KP operator, the improved KP operator has the output range , which is because the output of MSMA actuator is unidirectional [16].

Figure 2: The modified KP operator.

It can be seen from Figures 1 and 2 that the KP operator has some similar characters with the hysteresis loop of MSMA actuator. The KP model proposed in this paper can be expressed as an integral of the KP operators over a specific domain

where and are the input and output, respectively, is the operator which can transform to , and is the Preisach plane which is a triangular domain as Figure 3 shows, and , is the KP operator, is the weighting function to scale the output of the KP operator, and is the historical extreme values of the KP operator.

Figure 3: Preisach plane.

By dividing the Preisach plane into a mesh grid of as shown in Figure 3, each point in the Preisach plane represents a KP operator with its weighting value as shown in Figure 4. That means that there are KP operators in the KP model. An approximation of the integral KP model can be gotten as a parallel connection of a number of weighted KP operators

Figure 4: Schematic diagram of discrete KP model.

The function of KP operator is given in (3)

The function of is

where the ridge function for two boundaries of KP operator can be expressed by

where .

3. Identification Methods of the Hysteresis Model

From (2), it can be known that the KP model can be expressed by a parallel connection of a number of weighted KP operators like a neural network; that means that there are weighting parameters needed to be identified, and the weighting parameters identification can affect the modeling accuracy directly. In this paper, two identification methods are adopted to identify the established KP model: improved gradient correction algorithm and variable step-size recursive least square estimation algorithm.

3.1. The Improved Gradient Correction Algorithm

The recurrence formulas of gradient correction parameter estimation for deterministic problems are as follows [17]:

where is the parameters evaluation vector at time , is the weight matrix, is the vector of the KP operators’ values at time , and is the actual output at time .

Equation (6) is the traditional recurrence formula of gradient correction parameter estimation for deterministic problems; however, with this identification method, the modeling accuracy is not satisfactory, and therefore, an improved gradient correction algorithm is proposed in this paper. The algorithm of the improved gradient correction algorithm is as follows

(1) Set an expectation error . Assume that there are sets of data which are used to do the identification. (2) Assign a zero vector to the initial value , and assign a zero vector to . Set , where is the cycle number.

(i) If , go to Step (3), else set and continue, where is a maximum cycle number used to avoid the endless loop.(ii) Get the modeling error , , if , , go to Step (3), and else continue. In this step, the weight matrix , and is a unit matrix. (iii) If , set , , and and back to Step (i), else and back to Step (ii).

(3) Set , and is the weighting parameters gotten after the identification.

In the proposed improved gradient correction algorithm above, the expectation error and the maximum cycle number can be set according to the requirements.

3.2. Variable Step-Size Recursive Least Square Estimation Algorithm

In the system identification methods, least square method is the most widely used one, the recursive least square method is especially popular. In the recursive least square method, the parameters evaluation is updated at every time when a new set of observation data is gotten. In order to reduce the computations in the identification process, a variable step-size recursive least square estimation algorithm is adopted. The algorithm principle of the variable step-size recursive least square estimation algorithm is as follows [18].

Assume that the formula of least square estimation algorithm is

where is weighting parameters gotten after the identification, is the weighted factor, is the values of the KP operators, and is the actual output. Assume that , , and is the step-size which is an integer greater than 0, then there is

where and , and according to (7), there is

therefore , and assume that , then there is

According to the matrix inversion formula, there is

According to (10) and (11), the variable step-size recursive least square estimation algorithm can be derived

4. Simulation Study

In this section, simulation experiments are investigated to verify the validity of the proposed model method with different identification methods.

When the improved gradient correction algorithm is adopted and the number of KP model’s discretization lines is set to , the expectation error is set to , the maximum cycle number , and the simulation result is shown in Figures 5, 6, and 7. Figure 5 shows the comparison of actual hysteresis loop and simulation hysteresis loop, Figure 6 shows the actual output and model output of MSMA actuator, and Figure 7 shows the modeling error.

Figure 5: Comparison of actual hysteresis loop and simulation hysteresis loop (with the improved gradient correction algorithm).

Figure 6: Actual output curve and model output curve (with the improved gradient correction algorithm).

Figure 7: Modeling error (with the improved gradient correction algorithm).

From the simulation results, it can be seen that the proposed KP model with improved gradient correction algorithm can be used to establish a good hysteresis model which has the ability to describe both the major hysteresis loop and inner hysteresis loops, and the maximum modeling error is 0.0104 mm.

When the variable step-size recursive least square estimation algorithm is adopted and the number of KP model’s discretization lines is set to , the step-size of the variable step-size recursive least square estimation algorithm is set to and then the simulation result is shown in Figures 8, 9, and 10.

Figure 8: Comparison of actual hysteresis loop and simulation hysteresis loop (with the variable step-size recursive least square estimation algorithm).

Figure 9: Actual output curve and model output curve (with the variable step-size recursive least square estimation algorithm).

Figure 10: Modeling error (with the variable step-size recursive least square estimation algorithm).

The simulation results show that with the variable step-size recursive least square estimation algorithm, the established hysteresis model can also have a good performance in describing the hysteresis between input and output of MSMA actuator, and the maximum modeling error is 0.0209 mm.

In the simulation experiments above, there are some optional parameters as the number of KP model’s discretization lines , expectation error and the maximum cycle number in the improved gradient correction algorithm, and the step-size in the variable step-size recursive least square estimation algorithm. According to the mathematical principles above, the modeling accuracy can be further improved by changing these optional parameters.

With the improved gradient correction algorithm, Figure 11 shows the modeling error when the optional parameters have different values.

(a) , , and

(b) , , and

(c) , , and

(d) , , and

Figure 11: Modeling error when the optional parameters have different values (with the improved gradient correction algorithm).

From Figure 11, it can be seen that by increasing the number of discretization lines and the maximum cycle number , the modeling error can be reduced. When , , and , as shown in Figure 11(d), the maximum modeling error can be reduced to 0.0056 mm.

With the variable step-size recursive least square estimation algorithm, Figure 12 shows the modeling error when the optional parameters have different values.

(a) ,

(b) ,

(c) ,

(d) ,

Figure 12: Modeling error when the optional parameters have different values (with the variable step-size recursive least square estimation algorithm).

Just as expected, Figure 12 shows that when the number of discretization lines is increased or the step-size is reduced, the modeling error is reduced. When and , the maximum modeling error can be reduced to 0.001 mm.

The experiment results in this section can verify the validity of the proposed hysteresis model in this paper, and from the comparison of the simulation results of the two identification algorithms, it can be seen that the variable step-size recursive least square estimation algorithm has a better performance in the identification, and the modeling accuracy can be adjusted by adjusting the optional parameters in the proposed model.

5. Conclusion

A mathematic modeling method based on KP model is proposed to describe the hysteresis of the MSMA actuator in this paper. In order to identify the weighting parameters in the KP model, an improved gradient correction algorithm and a variable step-size recursive least square estimation algorithm are adopted, and the simulation results have proven the validity of the proposed hysteresis model. By adjusting the optional parameters of this hysteresis model, the modeling accuracy can be further improved, and in the simulation experiment, the modeling error can be reduced to 0.001 mm with the variable step-size recursive least square estimation algorithm. A satisfactory hysteresis model can be gotten if the optional parameters are adjusted according as the requirements of modeling accuracy and computational complexity.

Acknowledgments

This work is partially supported by the National Natural Science Foundation of China, Grant no. 51105170 and Program of Science and Technology Development Plan of Jilin province of China, Grant no. 20115015.

References

T. Krenke, S. Aksoy, E. Duman et al., “Hysteresis effects in the magnetic-field-induced reverse martensitic transition in magnetic shape-memory alloys,” Journal of Applied Physics, vol. 108, no. 4, Article ID 043914, 2010.

C. Visone, D. Davino, and A. A. Adly, “Vector preisach modeling of magnetic shape memory materials oriented to power harvesting applications,” IEEE Transactions on Magnetics, vol. 46, no. 6, pp. 1848–1851, 2010.

S. Flaga, J. Pluta, and B. Sapiński, “Pneumatic valves based on magnetic shape memory alloys: potential applications,” in Proceedings of the 12th International Carpathian Control Conference (ICCC '11), pp. 111–114, Krakow, Poland, May 2011.

M. Kamlah, “Ferroelectric and ferroelastic piezoceramics—modeling of electromechanical hysteresis phenomena,” Continuum Mechanics and Thermodynamics, vol. 13, no. 4, pp. 219–268, 2001.

M. Al Janaideh, S. Rakheja, and C. Y. Su, “An analytical generalized Prandtl-Ishlinskii model inversion for hysteresis compensation in micropositioning control,” IEEE/ASME Transactions on Mechatronics, vol. 16, no. 4, pp. 734–744, 2011.

A. Falvo, F. Furgiuele, and C. Maletta, “Hysteresis modeling of two-way shape memory effect in NiTi alloys,” Meccanica, vol. 43, no. 2, pp. 165–172, 2008.

P. B. Nguyen and S. B. Choi, “Compensator design for hysteresis of a stacked PZT actuator using a congruency-based hysteresis model,” Smart Materials and Structures, vol. 21, no. 1, Article ID 015009, 2012.

C. Natale, F. Velardi, and C. Visone, “Identification and compensation of Preisach hysteresis models for magnetostrictive actuators,” Physica B, vol. 306, no. 1–4, pp. 161–165, 2001.

Y. Bernard, E. Mendes, and F. Bouillault, “Dynamic hysteresis modeling based on Preisach model,” IEEE Transactions on Magnetics, vol. 38, no. 2, pp. 885–888, 2002.

M. Brokate, S. MacCarthy, A. Pimenov, A. Pokrovskii, and D. Rachinskii, “Modelling energy dissipation due to soil-moisture hysteresis,” Environmental Modeling and Assessment, vol. 16, no. 4, pp. 313–333, 2011.

V. Preda, M. F. Ionescu, V. Chiroiu, and T. Sireteanu, “A Preisach model for the analysis of the hysteretic phenomena,” Revue Roumaine des Sciences Techniques, vol. 55, no. 3, pp. 1–12, 2010.

L. Riccardi, D. Naso, B. Turchiano, and H. Janocha, “Adaptive modified Prandtl-Ishlinskii model for compensation of hysteretic nonlinearities in magnetic shape memory actuators,” in Proceedings of the 37th Annual Conference on IEEE Industrial Electronics Society, pp. 56–61, Bari, Italy, November 2011.

M. A. Janaideh and P. Krejčí, “Prandtl-Ishlinskii hysteresis models for complex time dependent hysteresis nonlinearities,” Physica B, vol. 407, no. 9, pp. 1365–1367, 2012.

G. V. Webb, D. C. Lagoudas, and A. J. Kurdila, “Hysteresis modeling of SMA actuators for control applications,” Journal of Intelligent Material Systems and Structures, vol. 9, no. 6, pp. 432–448, 1998.

M. C. Deng, C. G. Jiang, A. Inoue, and C. Y. Su, “Operator-based robust control for nonlinear systems with Prandtl-Ishlinskii hysteresis,” International Journal of Systems Science, vol. 42, no. 4, pp. 643–652, 2011.

H. Y. Luo, Y. J. Liao, E. Abel, Z. G. Wang, and X. Liu, “Hysteresis behaviour and modeling of SMA actuators,” in Shape Memory Alloys, C. Cismasiu, Ed., pp. 61–79, 2010.

Y. B. Hou, M. Wang, and L. Q. Wang, System Identification and Its MATLAB Simulation, Science Press, Beijing, China, 2004.

A. G. Han, Y. F. Zhou, Y. Dai, and X. D. Chen, “Improvement of the system online identification algorithm,” Journal of Huazhong University of Science and Technology, vol. 33, no. 10, pp. 4–6, 2005.