About this Journal Submit a Manuscript Table of Contents
Journal of Applied Mathematics
Volume 2013 (2013), Article ID 427213, 8 pages
Research Article

Modeling and Control for Giant Magnetostrictive Actuators with Rate-Dependent Hysteresis

School of Automation Science and Electrical Engineering, Beijing University of Aeronautics & Astronautics, Beijing 100191, China

Received 21 February 2013; Accepted 29 July 2013

Academic Editor: Constantinos Siettos

Copyright © 2013 Ping Liu 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.


The rate-dependent hysteresis in giant magnetostrictive materials is a major impediment to the application of such material in actuators. In this paper, a relevance vector machine (RVM) model is proposed for describing the hysteresis nonlinearity under varying input current. It is possible to construct a unique dynamic model in a given rate range for a rate-dependent hysteresis system using the sinusoidal scanning signals as the training set input signal. Subsequently, a proportional integral derivative (PID) control scheme combined with a feedforward compensation is implemented on a giant magnetostrictive actuator (GMA) for real-time precise trajectory tracking. Simulations and experiments both verify the effectiveness and the practicality of the proposed modeling and control methods.

1. Introduction

Magnetostrictive actuator has broad applications in the super-precision tracking system and microvibration control systems due to its prominent capabilities of high energy densities, large stroke, and fast response. Figure 1 illustrates the components of a GMA. The stroke and output force is provided by the Terfenol-D rod in response to a varying magnetic field generated by the surrounding solenoid coils. A biased magnetic field generated by the permanent magnet and the prestress are introduced to produce bidirectional actuation and to improve performance of the Terfenol-D rod, respectively [1].

Figure 1: Section view of GMA.

However, the GMA’s wider application is restricted due to an inherent property of the magnetostrictive materials, the hysteresis nonlinearity, which can cause undesirable inaccuracy, oscillation, or even instability of the systems. The hysteresis is embodied as a nonlinear relationship between the input voltage and the output displacement and is dependent not only on the amplitude but also on the frequencies of the input signals; hence, the hysteresis system is rate-dependent [2]. As shown in Figure 2, the shape of the hysteresis loops changes significantly with the frequencies of the input signals. Consequently, it is a major challenge to capture the complicated rate-dependent hysteretic behavior precisely, and the modeling, identification, and control of the rate-dependent hysteresis have attracted much attention for its sophisticated nature and extensive applications.

Figure 2: Rate-dependent hysteresis of GMA.

The existing techniques for the rate-dependent hysteresis modeling can be classified into three categories. The first category is based on phenomenological models. Recent propositions include the modified Preisach model [3], the modified Prandtl-Ishlinskii model [4], and the generalized Prandtl-Ishlinskii model [5]. The second category contains the physical models. For instance, JinHyoung et al. developed the semilinear Duhem model to describe the rate-dependent hysteresis [6, 7]. The third category is the intelligent methods, such as artificial neural network (ANN) [810], fuzz tree (FT) [11], and all providing efficient ways to model the rate-dependent hysteresis. Control of smart structures with consideration of hysteresis nonlinearity has also received much attention. A very effective approach in dealing with the hysteresis in control system is to find an accurate model of the hysteresis and then to use its inverse compensation to cancel out the hysteresis [12, 13]. In addition, some feedback control scheme including adaptive control [14], robust control [15], and optimal control [16] were also adopted in the controller synthesis.

In this paper, we define the ability for a model to correctly describe the data of some other frequencies not in the training set as the generalization ability of frequency. A universal model with generalization ability of frequency and high precision is always required in a practical application for convenience and efficiency.

In this paper, we study the modeling and trajectory tracking control of a GMA with rate-dependent hysteresis nonlinearity under variable input signal frequencies. First of all, an RVM model, which is an inherent online machine learning technique using a more flexible and sparser function without additional regularization parameters, is employed as a universal model to capture the rate-dependent hysteresis. Then, the compensation with a PID feedback algorithm based on the proposed hysteresis model is applied to the GMA to cancel out hysteresis for real-time trajectory tracking control. Experiments are carried out to validate the feasibility and effectiveness of the proposed model and the control schemes.

2. Rate-Dependent Hysteresis Model

2.1. RVM Model

The relevance vector machine (RVM) introduced by Tipping [17, 18] is a probabilistic model similar to the support vector machine (SVM), but where the training takes place in a Bayesian framework, and predictive distributions of the outputs instead of point estimates are obtained. RVM evades the complexity by producing models that have both structure and parameterization processes that are appropriate to the information content of the data.

Given a data set of input-target pairs considering scalar-valued target functions only, we follow the standard probabilistic formulation and assume that the targets are samples from the model with additive noise as follows: where is independent sample from some noise process which is further assumed to be mean-zero Gaussian with variance , is the design matrix with , wherein , is kernel function, and is the weight vector. In this work, the radial basis function (RBF) kernel function is used. Due to the assumption of independence of the , the likelihood of the complete data set can be written as

With as many parameters in the model as training examples, we would expect maximum likelihood estimation of and from (2) to lead to severe overfitting. To avoid this, a common approach is to impose some additional constraints on the parameters, for example, through adopting a Bayesian perspective, and constrain the parameters by defining an explicit prior probability distribution over them.

We encode a preference for smoother (less complex) functions by making the popular choice of a zero-mean Gaussian prior distribution over : with the hyperparameters . Importantly, there is an individual hyperparameter associated independently with every weight, moderating the strength of the prior thereon. According to Bayes’ rule, the posterior distribution over the weights is thus given by where the posterior covariance and mean are, respectively; with . From the marginal likelihood, the probability distribution of is given by with . Values of and cannot be obtained in closed form, and here we summarize formulae for their iterative re-estimation. For , differentiation of (6), equating to zero and rearranging, gives where is the th posterior mean weight from (5), and we have defined the quantities by with the th diagonal element of the posterior weight covariance from (5) computed with the current and values. For the noise variance , differentiation leads to re-estimation

Note that the in the denominator refers to the number of data examples and not the number of basis functions.

The learning algorithm thus proceeds by repeated application of (7) and (9), concurrent with updating of the posterior statistics and from (5) until some suitable convergence criteria have been satisfied.

In practice, during re-estimation, we generally find that many of the tend to infinity. From (4), this implies that becomes highly peaked at zero; that is, we are a posteriori certain that those are zero. The corresponding basis functions can thus be pruned, and sparsity is realized.

2.2. Modeling Performance

In this paper, the multivalued mapping is transformed into a single-valued one by employing the current and previous inputs and previous outputs as exogenous inputs, the nonlinear regression model of the form which is called feedforward model, NARX model, or series-parallel model [19], where and are the input and the true output of the system at discrete time instant , is the estimated output of the model at time instant , and and are the dynamic orders.

The effects of the order ( and ) selection on the modeling results are investigated based on multisimulation. It is found that the larger the model order, the lower the training error. However, the testing error is not monotonic with respect to the order selection. We choose the order ( and ) to make a compromise between the training and testing errors. Then, we choose , .

We first apply the sinusoidal scanning signals data from 1 Hz to 100 Hz to obtain the nonlinear system based on RVM, then we test the rate-dependent model using respectively the single-signal data from 1 Hz to 100 Hz and some sets of compound-signal data. Figure 3 gives the sinusoidal scanning signals data generated by the MATLAB order and the modeling result of the model data.

Figure 3: The input signal and the hysteresis loop of the sinusoidal scanning signal.

Our experimental data are obtained at a sampling frequency of 10 KHz from GMA manufactured by the Beijing University of Aeronautics and Astronautics [20]. At every frequency, we take 500 pairs of data, , from all experimental data and from the model simulation, respectively. The modeling results can be evaluated through the root mean square error () and the relative error (), which are defined, respectively, as where is the number of data, is the experimental output, and is the output calculated from the model.

Figure 4 and Table 1 show the comparison between the rate-dependent hysteresis loop measured by experiments and those simulated based on the RVM modeling method at several single-signal and compound-signal inputs.

Table 1: Modeling results.
Figure 4: Results of RVM modeling.

Remarks (1) The number of the modeling data must be enough big, then the frequency range of the sinusoidal scanning signal can cover the needed frequency.
(2) The experimental data are compared with the modeling results in Figure 4 and Table 1 to reveal the outstanding generalization performance of the universal rate-dependent RVM hysteresis model.

2.3. Trajectory Tracking Control

The experimental equipment was constructed to achieve 1-DOF trajectory tracking control as shown in Figure 5. The GMA was manufactured by Zhang et al. [20], measuring Φ50 mm × 200 mm with Φ7 mm × 80 mm Terfenol-D rod, and was driven by GF-20 amplifier operating in current mode, which was controlled by a computer with dSPACE (DS1103) control board. The displacement was measured by an eddy current sensor with a 0.1 resolution. The input signal was generated by the computer with dSPACE and transported to the GMA through amplifier; then, the output displacement generated by the GMA was transported by the eddy current sensor to the computer, and then we got the experiment data .

Figure 5: The diagram of experiment equipment.

The basic idea for controller synthesis for hysteresis system is to design a right inverse model to cancel the hysteresis nonlinearity. However, there is no any exact analytical inverse solution to the RVM model, and there are deviations with the variable input frequency in the direct inverse compensation control. In order to improve the accuracy of tracking control, a PID feedback controller combined with compensation in the feed-forward loop is used for real-time tracking control. Figure 6 shows the block diagram of the tracking control system. The hysteresis nonlinearity could be eliminated by the inverse compensation in the feed-forward loop, and the PID feedback controller is adopted to deal with remaining nonlinear uncertainties generated by the inverse compensation with the change on the input frequency, which can be seen as the disturbance.

Figure 6: The diagram of the combined control scheme.

The nonlinear system (10) is said to be invertible at if there is a subset of , such that for , for any distinct inputs , [21]. So we construct the inverse compensation of the model based on

Although the inverse compensator is not the analytical inverse solution of the original model, the compensation effect is also good, which can be seen in Figure 7. And the small part uncompensated can be treated by the PID controller.

Figure 7: The compensation effect of the inverse model.

The combined feedback controller, suitable for both single signal and compound signal, is fixed with the different input signal. Figure 8 and Table 2 show the trajectory tracking control results of the combined control scheme. For convenience, the time intervals shown in the figures are shorter with the increasing frequency. Note that the sampling frequency of tracking control is the same as that used for modeling.

Table 2: Tracking control results.
Figure 8: Results of combined tracking control.

Remarks (1) The reference signal should be constrained in the range of the output of the GMA between and ; otherwise, the overtopped reference signal may not be attained. To solve the problem, one method is to predict the range of GMA output during the modeling process.
(2) The RVM model we derived is a numerical model, and an analytical form is not available, so the PID controller’s parameters cannot be directly computed. Thus, we adjust the parameters of PID controller through offline simulation.
(3) The tracking control results in Figure 8 and Table 2 reveal that the combined controller is effective at both single-input and compound-input frequencies, and the RE all less than 13%, the allowable range of the project.

3. Conclusion

The main contribution of this paper is to propose a systematic approach for applications of a rate-dependent hysteresis nonlinear model to precise tracking control of the GMA. We model the rate-dependent hysteresis based on the RVM in a wide frequency range. Based on the model, an inverse compensation in the feed-forward loop combined with a PID feedback controller has been developed. Experimental results demonstrate that the RVM model can capture the rate-dependent hysteresis effects, and our identification and control schemes are effective.

Conflict of Interests

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work; there is no professional or other personal interest of any nature or kind in any product, service and/or company that could be construed as influencing the position presented in, or the review of, this paper.


This work is supported by National Nature Science Foundation of China (nos. 91016006, and 91116002) and the Fundamental Research Funds for the Central Universities (nos. 30420111109, and 30420120305).


  1. Y. Ma and J. Mao, “Application of a fuzzy tree based dynamic hysteresis model to precise tracking control for a smart structure,” in Proceedings of the 7th Asian Control Conference (ASCC '09), pp. 489–494, August 2009. View at Scopus
  2. J. Oh and D. S. Bernstein, “Semilinear Duhem model for rate-independent and rate-dependent hysteresis,” IEEE Transactions on Automatic Control, vol. 50, no. 5, pp. 631–645, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  3. Y. Yu, Z. Xiao, N. G. Naganathan, and R. V. Dukkipati, “Dynamic Preisach modelling of hysteresis for the piezoceramic actuator system,” Mechanism and Machine Theory, vol. 37, no. 1, pp. 75–89, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  4. U. X. Tan, T. L. Win, C. Y. Shee, and W. T. Ang, “Rate-dependent hysteresis model of piezoelectric using singularity free Prandtl-Ishlinskii model,” in Proceedings of the IEEE International Symposium on Computational Intelligence in Robotics and Automation (CIRA '07), pp. 356–361, June 2007. View at Publisher · View at Google Scholar · View at Scopus
  5. M. A. Janaideh, S. Rakheja, and C.-Y. Su, “Experimental characterization and modeling of rate-dependent hysteresis of a piezoceramic actuator,” Mechatronics, vol. 19, no. 5, pp. 656–670, 2009. View at Publisher · View at Google Scholar · View at Scopus
  6. J. Oh and D. S. Bernstein, “Identification of rate-dependent hysteresis using the semilinear Duhem model,” in Proceedings of the American Control Conference (AAC '04), pp. 4776–4781, July 2004. View at Publisher · View at Google Scholar · View at Scopus
  7. J. Oh and D. S. Bernstein, “Semilinear Duhem model for rate-independent and rate-dependent hysteresis,” IEEE Transactions on Automatic Control, vol. 52, no. 3, pp. 576–582, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  8. S. Yu, B. Shirinzadeh, G. Alici, and J. Smith, “Sliding mode control of a piezoelectric actuator with neural network compensating rate-dependent hysteresis,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 3641–3645, April 2005. View at Publisher · View at Google Scholar · View at Scopus
  9. R. Dong, Y. Tan, H. Chen, and Y. Xie, “A neural networks based model for rate-dependent hysteresis for piezoceramic actuators,” Sensors and Actuators A, vol. 143, no. 2, pp. 370–376, 2008. View at Publisher · View at Google Scholar · View at Scopus
  10. J. P. Lien, A. York, T. Fang, and G. D. Buckner, “Modeling piezoelectric actuators with Hysteretic Recurrent Neural Networks,” Sensors and Actuators A, vol. 163, no. 2, pp. 516–525, 2010. View at Publisher · View at Google Scholar · View at Scopus
  11. J. Mao and H. Ding, “Intelligent modeling and control for nonlinear systems with rate-dependent hysteresis,” Science in China, vol. 52, no. 4, pp. 656–673, 2009. View at Publisher · View at Google Scholar · View at Scopus
  12. R. C. Smith, “Inverse compensation for hysteresis in magnetostrictive transducers,” Mathematical and Computer Modelling, vol. 33, no. 1–3, pp. 285–298, 2001. View at Publisher · View at Google Scholar · View at Scopus
  13. X. Zhang, Y. Tan, M. Su, and Y. Xie, “Neural networks based identification and compensation of rate-dependent hysteresis in piezoelectric actuators,” Physica B, vol. 405, no. 12, pp. 2687–2693, 2010. View at Publisher · View at Google Scholar · View at Scopus
  14. X. Chen, T. Hisayama, and C.-Y. Su, “Pseudo-inverse-based adaptive control for uncertain discrete time systems preceded by hysteresis,” Automatica, vol. 45, no. 2, pp. 469–476, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  15. Y. Feng, C. A. Rabbath, H. Hong, M. Al Janaideh, and C.-Y. Su, “Robust control for shape memory alloy micro-actuators based flap positioning system,” in Proceedings of the American Control Conference (ACC '10), pp. 4181–4186, July 2010. View at Scopus
  16. W. S. Oates and R. C. Smith, “Nonlinear optimal control techniques for vibration attenuation using magnetostrictive actuators,” Journal of Intelligent Material Systems and Structures, vol. 19, no. 2, pp. 193–209, 2008. View at Publisher · View at Google Scholar · View at Scopus
  17. C. M. Bishop and M. E. Tipping, “Variational relevance vector machines,” in Proceedings of the 16th Conference on Uncertainty in Artificial Intelligence (UAI '00), pp. 46–53, San Francisco, CA, USA, 2000.
  18. M. E. Tipping, “Sparse Bayesian learning and the relevance vector machine,” Journal of Machine Learning Research, vol. 1, no. 3, pp. 211–244, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  19. J. A. K. Suykens and J. Vandewalle, “Recurrent least squares support vector machines,” IEEE Transactions on Circuits System, vol. 47, no. 7, pp. 1109–1114, 2000.
  20. H. B. Zhang, C. B. Jiang, Z. B. Wang, and H. B. Xu, “Effect of compressive stress on magnetostriction hysteresis of 110 oriented Tb0.29Dy0.48Ho0.23Fe2 crystal,” Journal of Alloys and Compunds, vol. 475, no. 1-2, pp. 35–37, 2009.
  21. K. J. Hunt and D. Sbarbaro, “Neural networks for nonlinear internal model control,” IEE Proceedings D of Control Theory and Applications, vol. 138, no. 5, pp. 431–438, 1991.