Next Article in Journal
Fabrication of Optical Fibers with Multiple Coatings for Swelling-Based Chemical Sensing
Next Article in Special Issue
A Review of Microrobot’s System: Towards System Integration for Autonomous Actuation In Vivo
Previous Article in Journal
Visualization of Local Concentration and Viscosity Distribution during Glycerol-Water Mixing in a Y-Shape Minichannel: A Proof-of-Concept-Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Digitized Representation of the Modified Prandtl–Ishlinskii Hysteresis Model for Modeling and Compensating Piezoelectric Actuator Hysteresis

1
School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 637459, Singapore
2
Invenco Group Limited, Auckland 1023, New Zealand
*
Author to whom correspondence should be addressed.
Submission received: 21 July 2021 / Revised: 4 August 2021 / Accepted: 5 August 2021 / Published: 10 August 2021
(This article belongs to the Special Issue Imaging-Guided Intelligent Micromachines)

Abstract

:
Piezoelectric actuators are widely used in micromanipulation and miniature robots due to their rapid response and high repeatability. The piezoelectric actuators often have undesired hysteresis. The Prandtl–Ishlinskii (PI) hysteresis model is one of the most popular models for modeling and compensating the hysteresis behaviour. This paper presents an alternative digitized representation of the modified Prandtl–Ishlinskii with the dead-zone operators (MPI) hysteresis model to describe the asymmetric hysteresis behavior of piezoelectric actuators. Using a binary number with n digits to represent the classical Prandtl–Ishlinskii hysteresis model with n elementary operators, the inverse model can be easily constructed. A similar representation of the dead-zone operators is also described. With the proposed digitized representation, the model is more intuitive and the inversion calculation is avoided. An experiment with a piezoelectric stacked linear actuator is conducted to validate the proposed digitized MPI hysteresis model and it is shown that it has almost the same performance as compared to the classical representation.

1. Introduction

Piezoelectric actuators are widely used in micromanipulation applications due to their rapid response and high repeatability, for example, microsurgical robots [1,2], autofocus optical systems [3,4,5], precise fabrication [6,7,8] and other applications [9,10,11]. One of the biggest challenges while using the piezoelectric actuators in dynamic applications is to model and compensate for the undesired complex hysteresis.
Various methods have been proposed for modeling and compensating the hysteresis behavior. The existing methods can be classified into physics-based models and phenomenology-based models [12]. Physics-based models [13] are often derived on the basis of physical principles of certain material or system properties [14,15,16]. Physics-based models often require a deep understanding of the causes of hysteresis and are often specific to the related properties.
The phenomenology-based model can also be roughly divided into three main groups [17]: differential-based models, such as the Bouc-wen [18] model, the Duhem [19] model, the Dahl, and the LuGre [20] model; neural-network models, such as the back propagation neural network based model [21], the gated recurrent unit based model [22], the neural network adaptive control method [23] and etc; operator-based models, such as the Preisach [24] model, the Krasnosel’skii–Pokrovskii(KP) [25] model, the Maxwell-slip [26] model, the Prandtl–Ishlinskii(PI) [27] model, and their variations.
The PI hysteresis model is one of the most popular models because of its simple implementation and analytical inversion [28]. Ang [29] concluded the relationship between the rate and the operators and modified the PI hysteresis model to model the rate-dependent hysteresis behavior. Tan [30] extended the PI hysteresis model to compensate for the ill-conditioned hysteresis behavior with a negative gradient. Mohammad [31] proposed a nonlinear play operator for modeling the asymmetric hysteresis of the Shape Memory Alloy (SMA) operators. Kuhnen [32] proposed to use the dead-zone operators for modeling the memory-free asymmetric hysteresis behavior.
In the above literature, the design procedure generally consists of parameters identification of the description model and the construction of the inverse model as a desired compensator. The mathematical complexity of the identification and inversion problem depends on the selected modeling approaches [33]. The amount of calculation may increase significantly with an increase in the number of elementary operators.
In this paper, we present an alternative digitized representation with which the inverse model can be easily constructed to compensate for the hysteresis behavior. The inversion calculation is avoided. The proposed representation is also more intuitive due to the geometric meanings of the parameters of the operators.
The paper is organized as follows: Section 2 describes the classical representation from which the proposed representation can be derived. Section 3 gives the experimental results with discussion. Section 4 covers the conclusion.

2. Hysteresis Mathematical Model

This section describes the classical representation and an alternative digitized representation of the classical PI hysteresis model. A similar representation of the dead-zone operators is also described, the construction of the inverse model with the proposed representation is also presented.

2.1. Classical PI Hysteresis Model

The classical PI hysteresis model has a series of elementary operators, the system output can be obtained by a weighted sum of n elementary operators with different thresholds. As shown in Figure 1, the working principles of the classical representation of the PI hysteresis model can be represented using the following equations:
y ( 0 ) = y 0
y ( t ) = i = 1 n w i · m a x { u ( t ) r i , m i n { u ( t ) + r i , y e i ( t T ) } }
where w i is the weight corresponding to the ith elementary operator, r i is the threshold of the ith elementary operator, n is the number of the elementary operators and y e i is the output of the ith elementary operator.
When the parameter r i and w i are identified, the hysteresis curve model in the ( u , y ) plane is shifted so that it is in the positive section of the plane. An example with three elementary operators is shown in Figure 2, b w i = 2 r i is the input value x i of turning point where the slope of the hysteresis curve is changed. The k t h output value y k can be formulated as follows [34]:
y k = i = 1 k ( b w k + 1 b w i ) · w i
The Figure 2 can also be regarded as a series of 4 segments. Let ( Δ x i , Δ y i ) by the input and output value of each of the n segments, the value of the Δ x i is the change of b w i while the value of the Δ y i is the change of y k . The slope g r a d k of the kth segment can be calculated using the following equation:
g r a d 1 = 0 ;
g r a d k = i = 1 k w i ;
As b w i = 2 r i , the value of Δ x i can be calculated using the following equations:
Δ x 1 = 2 · r 1 ;
Δ x i = 2 · ( r i + 1 r i ) ;
Δ x n = 2 · ( u A r n 1 ) ;
The value of Δ y k can then be obtained using the following equations:
Δ y 1 = 0 ;
Δ y k = ( i = 1 i w ( i ) ) · Δ x k
The paired property ( Δ x i , Δ y i ) represents the change of input and output value of each segment of the hysteresis curve and is thus more intuitive and the parameters are easier to be identified. The value of the paired property can be derived from the parameters of the classical representation. With this paired property, an alternative digitized representation of the classical PI hysteresis model is described in the next section.

2.2. Digitized Classical PI Hysteresis Model

The state of the classical PI hysteresis model with n operators is represented with a binary number B = b n b 2 b 1 with n digits. Each digit of the binary number is a elementary operator with a binary state that can be either 0 or 1. Each elementary operator b i also has a paired property ( Δ x i , Δ y i ) , where Δ x i represents a change in the system input and Δ y i represents a change in the system output. An example with two operators is shown in Figure 3. With an increase in the number of operators, the system output and input can be approximated to be a linear relationship within each operator. The paired properties H of a PI hysteresis model with n elementary operators can be expressed in one equation:
H = Δ x 1 Δ x n Δ y 1 Δ y n
The model also has a paired property ( X B , Y B ) that corresponds to its state. The state of the system always changes to either one of its two neighbors: upper state or lower state. Let b 1 be the least significant digit and b n be the most significant digit, the change of the state is always from the least significant digit. Let the least significant digit with a binary state of 0 at the current state be the A digit, let the least significant digit with a binary state of 1 at the current state be the Z digit.
When the state of the system changes to its upper neighbor, the binary state of the b A digit is changed from 0 to 1, the paired property of the system is increased by ( Δ x A , Δ y A ) . Denote X B t , Y B t be the system input and output of the current state and X B t + 1 , Y B t + 1 be the system input and output after a change, the change to its upper neighbor can be shown in the following equations:
X B t + 1 = X B t + Δ x A
Y B t + 1 = Y B t + Δ y A
Similarly, when the state of the system changes to its lower neighbor, the binary state of the b Z digit is changed from 1 to 0, the paired property of the system is decreased by ( Δ x Z , Δ y Z ) . Denote X B t , Y B t be the system input and output of the current state, and X B t + 1 , Y B t + 1 be the system input and output after a change, the change to its lower neighbor can be shown in the following equations:
X B t + 1 = X B t Δ x Z
Y B t + 1 = Y B t Δ y Z
The traditional PI hysteresis model with n elementary operators has 2 n possible states. A simple example with 3 operators can be shown in Figure 4. The modeling and control accuracy may be improved with an increase in the number of operators.
When the system input is changed to x with the current state B T , the system output can be calculated with the following procedure:
1. If the system input value, x, is greater than or equal to X B T and less than X B T + Δ x A , the value of the system output y is calculated using (16).
2. If the system input value, x, is greater than X B T + Δ x A , the state of the system is changed to its upper neighbor by setting the b A digit from 0 to 1, the ( X B T , Y B T ) is updated using (12) and (13). The procedure is then repeated.
3. If the system input value, x, is less than X B T and greater than X B T Δ x Z , the value of the system output y is calculated using (17).
4. If the system input value, x, is is less than X B T Δ x Z , the state of the system is changed to its lower neighbor by setting the b Z digit, the ( X B T , Y B T ) is updated using (14) and (15). The procedure is then repeated.
y = Y B T + Δ y A Δ x A × ( x Δ X B T )
y = Y B T + Δ y Z Δ x Z × ( x Δ X B T )

2.3. Digitized Dead-Zone Operators

With the dead-zone operators, the modified PI hysteresis model can model the memory-free asymmetric hysteresis behavior. The state of the dead-zone operators can be represented with a number C. The n paired properties can be used to help describe the values of the n dead-zone operators as shown in (18). An example with two operators is shown in Figure 5.
S = Δ x 1 Δ x m Δ y 1 Δ y m
The model also has a paired property ( X C , Y C ) that corresponds to its state, the relationship between the paired property and the state of the model is described in the following equations:
X C = X 0 + i = 1 C Δ X i ; C = 1 , , m
Y C = Y 0 + i = 1 C Δ Y i ; C = 1 , , m
When the system input is changed to x with the current state C, the system output can be calculated with the following strategy:
1.
If the system input value, x, is greater than or equal to X C and less than X C + Δ x C , the value of the system output y is calculated using (21).
2.
If the system input value, x, is greater than X C + Δ x C and C is less than m, the state of the system is changed to its upper neighbor by increasing C, the ( X C , Y C ) is updated by adding ( Δ x C , Δ y C ) . The procedure is then repeated.
3.
If the system input value, x, is less than X C and C is greater than 1, the state of the system is changed to its lower neighbor by decreasing C, the ( X C , Y C ) is updated by subtracting ( Δ x C , Δ y C ) . The procedure is then repeated.
y = Y C + Δ y C Δ x C × ( x Δ X C )

2.4. Inverse Model

The key idea of the parameters identification of the inverse model is to find the reflection of the resultant hysteresis curves about the 45 line. A linear response is obtained by cascading the inverse hysteresis operator Γ 1 as a feedfoward controller with the actual hysteresis. The proposed model for the inverse feedfoward controller is illustrated in Figure 6.
Unlike previous work, the parameters of the inverse model using the proposed representation can be found by simply exchanging Δ x with Δ y within the paired property, as shown in the following equations:
H 1 = Δ y 1 Δ y n Δ x 1 Δ x n
S 1 = Δ y 1 Δ y m Δ x 1 Δ x m
Due to its recursive nature, the calculation of the parameters identification of the inverse model does not change much with an increase in the number of the operators, while that of the modified PI hysteresis model is proportional to the number of the elementary operators. With the proposed digitized representation, the inversion calculation is avoided.

3. Experimental Results

In this section, the asymmetric hysteresis behavior of a piezoelectric stacked linear actuator, MPO-050100 (Nanofaktur, Stuttgart, Germany), is modeled using the classical representation as well as the proposed representation. Two experiments have been performed. The first experiment is to identify the parameters and to model the hysteresis behavior of the piezoelectric actuator. The second experiment is to test the performance of the proposed representation deploying the inverse model as a feedfoward compensator.

3.1. Experimental Setup

As seen from Figure 7, the input for the driver is a 10 V peak to peak, 1 Hz sinusoidal wave, is produced by a 16-bit D/A card. The input is passed to an amplifier, EMO-050100 (Nanofaktur, Stuttgart, Germany). The piezoelectric stacked linear actuator will deform, and the output displacement will be measured by a laser sensor, ZW-7010 (Omron, Kyoto, Japan) and is then converted to an analog signal which is received by a 16-bit A/D card. The experimental setup is shown in Figure 8.

3.2. Experimental Results

The hysteresis behavior of the piezoelectric stacked linear actuator under periodic control inputs is recorded. The parameters of the digitized representation and the classical representation of the modified PI model with the dead-zone operators are then identified with n = 25 , the identified parameters are used to model the hysteresis behavior of the piezoelectric stacked linear actuator. With the least square fitted models, the modeling results using the classical and the digitized representation with dead-zone operators are superimposed to the measured hysteresis behavior as shown in Figure 9. It is shown that the modeling results are almost the same using classical representation and the digitized representation. The root mean square error (RMSE) between the modeled hysteresis behavior and the measured hysteresis behavior are 0.50 μ m using either of the two representation with a stroke of 66.14 μ m. The modeling error is around 0.76 % of its stroke. This proves that the proposed representation has almost the same performance compared to the classical representation while modeling the asymmetric actuator hysteresis.
The inverse model is constructed and is then applied to compensate for the hysteresis behavior of the system at a changing amplitude sinusoidal wave with x d ( t ) = A × s i n ( 2 × π × t ) + 68.15 ( μ m). The compensation results are shown in Figure 10. The RMSE between the measured compensated position and the desired position is 0.36 μ m, with a maximum amplitude of 60 μ m. The error is around 0.6 % .

3.3. Discussion

Experimental results show that the proposed digitized representation performs almost the same as the classical representation while modeling the hysteresis behavior of the piezoelectric stacked linear actuator. This is because the digitized representation is derived from the classical representation and the two representation are equal mathematically. With n = 25 , the RMSE for modeling the piezoelectric actuator hysteresis is around 0.76 % . The compensation experimental results with the RMSE of around 0.6 % is even slightly better than the modeling results. This also proves that the proposed representation can compensate for the asymmetric actuator hysteresis. The description and compensation experimental results prove that the proposed digitized representation and its inverse model can model and compensate for the asymmetric asymmetric actuator hysteresis.
With the proposed digitized representation, the construction of the inverse model is less mathematically complicated and this improvement may be significant with a large value of n. The parameters of the proposed representation is also more intuitive with their geometric meanings; this may help with the parameters identification. The proposed representation is more intuitive and may be easier to be modified further.

4. Conclusions

The proposed digitized representation is an alternative to the classical representation. The inverse model can be easily constructed with the proposed representation and the inversion calculation is avoided. The proposed representation is also more intuitive. The proposed representation and its inverse model are validated with description and compensation experimental results.

Author Contributions

Conceptualization, C.Z., Y.N.A. and W.T.A.; Software, C.F.; Validation, C.Z.; Writing—original draft, C.Z.; Supervision, W.T.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Latt, W.; Tan, U.X.; Shee, C.; Ang, W. A compact hand-held active physiological tremor compensation instrument. In Proceedings of the 2009 IEEE/ASME International Conference on Advanced Intelligent Mechatronics, Singapore, 14–17 July 2009; pp. 711–716. [Google Scholar]
  2. MacLachlan, R.A.; Becker, B.C.; Tabarés, J.C.; Podnar, G.W.; Lobes, L.A.; Riviere, C.N. Micron: An actively stabilized handheld tool for microsurgery. IEEE Trans. Robot. 2011, 28, 195–212. [Google Scholar] [CrossRef] [Green Version]
  3. Bravo-Zanoguera, M.E.; Laris, C.A.; Nguyen, L.K.; Oliva, M.; Price, J.H. Dynamic autofocus for continuous-scanning time-delay-and-integration image acquisition in automated microscopy. J. Biomed. Opt. 2007, 12, 034011. [Google Scholar] [CrossRef] [PubMed]
  4. Aguirre, A.D.; Sawinski, J.; Huang, S.W.; Zhou, C.; Denk, W.; Fujimoto, J.G. High speed optical coherence microscopy with autofocus adjustment and a miniaturized endoscopic imaging probe. Opt. Express 2010, 18, 4222–4239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Chen, X.; Chen, Z.; Li, X.; Shan, L.; Sun, W.; Wang, X.; Xie, T.; Dong, S. A spiral motion piezoelectric micromotor for autofocus and auto zoom in a medical endoscope. Appl. Phys. Lett. 2016, 108, 052902. [Google Scholar] [CrossRef]
  6. Woronko, A.; Huang, J.; Altintas, Y. Piezoelectric tool actuator for precision machining on conventional CNC turning centers. Precis. Eng. 2003, 27, 335–345. [Google Scholar] [CrossRef] [Green Version]
  7. Okazaki, Y. A micro-positioning tool post using a piezoelectric actuator for diamond turning machines. Precis. Eng. 1990, 12, 151–156. [Google Scholar] [CrossRef]
  8. Abis, C.; Unal, F.; Mugan, A. Active vibration control with piezoelectric actuator on a lathe machine with a gain controller. In Proceedings of the 2011 IEEE International Conference on Mechatronics, Beijing, China, 7–10 August 2011; pp. 19–22. [Google Scholar]
  9. Junwu, K.; Zhigang, Y.; Taijiang, P.; Guangming, C.; Boda, W. Design and test of a high-performance piezoelectric micropump for drug delivery. Sens. Actuators A Phys. 2005, 121, 156–161. [Google Scholar] [CrossRef]
  10. Mori, K.; Munemoto, T.; Otsuki, H.; Yamaguchi, Y.; Akagi, K. A dual-stage magnetic disk drive actuator using a piezoelectric device for a high track density. IEEE Trans. Magn. 1991, 27, 5298–5300. [Google Scholar] [CrossRef]
  11. Chu, C.L.; Fan, S.H. A novel long-travel piezoelectric-driven linear nanopositioning stage. Precis. Eng. 2006, 30, 85–95. [Google Scholar] [CrossRef]
  12. Xie, S.L.; Liu, H.T.; Mei, J.P.; Gu, G.Y. Modeling and compensation of asymmetric hysteresis for pneumatic artificial muscles with a modified generalized Prandtl–Ishlinskii model. Mechatronics 2018, 52, 49–57. [Google Scholar] [CrossRef]
  13. Harrison, R.G. A physical model of spin ferromagnetism. IEEE Trans. Magn. 2003, 39, 950–960. [Google Scholar] [CrossRef]
  14. Haverkamp, R.; Reggiani, P.; Ross, P.J.; Parlange, J.Y. Soil water hysteresis prediction model based on theory and geometric scaling. Environ. Mech. Water Mass Energy Transf. Biosph. 2002, 129, 213–246. [Google Scholar]
  15. Zheng, J.; Li, Y.; Li, Z.; Wang, J. Transient multi-physics analysis of a magnetorheological shock absorber with the inverse Jiles–Atherton hysteresis model. Smart Mater. Struct. 2015, 24, 105024. [Google Scholar] [CrossRef]
  16. Knobloch, T.; Rzepa, G.; Illarionov, Y.Y.; Waltl, M.; Schanovsky, F.; Stampfer, B.; Furchi, M.M.; Mueller, T.; Grasser, T. A physical model for the hysteresis in MoS 2 transistors. IEEE J. Electron Devices Soc. 2018, 6, 972–978. [Google Scholar] [CrossRef]
  17. Yu, Y.; Zhang, C.; Wang, Y.; Zhou, M. Neural network-based iterative learning control for hysteresis in magnetic shape memory alloy actuator. IEEE/ASME Trans. Mechatron. 2021, 1. [Google Scholar] [CrossRef]
  18. Rakotondrabe, M. Bouc–Wen modeling and inverse multiplicative structure to compensate hysteresis nonlinearity in piezoelectric actuators. IEEE Trans. Autom. Sci. Eng. 2010, 8, 428–431. [Google Scholar] [CrossRef] [Green Version]
  19. Si, Z.Y.; Bai, X.X.; Qian, L.J.; Zhong, W.M. An enhanced Duhem model of magnetostrictive material-based actuators. In Behavior and Mechanics of Multifunctional Materials XIII; International Society for Optics and Photonics: Bellingham, WA, USA, 2019; Volume 10968, p. 1096818. [Google Scholar]
  20. Piatkowski, T. Dahl and LuGre dynamic friction models—The analysis of selected properties. Mech. Mach. Theory 2014, 73, 91–100. [Google Scholar] [CrossRef]
  21. Wu, Y.; Fang, Y.; Ren, X. A high-efficiency Kalman filtering imaging mode for an atomic force microscopy with hysteresis modeling and compensation. Mechatronics 2018, 50, 69–77. [Google Scholar] [CrossRef]
  22. Wu, Y.; Fang, Y.; Liu, C.; Fan, Z.; Wang, C. Gated recurrent unit based frequency-dependent hysteresis modeling and end-to-end compensation. Mech. Syst. Signal Process. 2020, 136, 106501. [Google Scholar] [CrossRef]
  23. Zhao, X.; Su, Q.; Chen, S.; Tan, Y. Neural network adaptive control of nonlinear systems preceded by hysteresis. J. Intell. Mater. Syst. Struct. 2021, 32, 104–112. [Google Scholar] [CrossRef]
  24. Song, G.; Zhao, J.; Zhou, X.; De Abreu-García, J.A. Tracking control of a piezoceramic actuator with hysteresis compensation using inverse Preisach model. IEEE/ASME Trans. Mechatron. 2005, 10, 198–209. [Google Scholar] [CrossRef]
  25. Li, Z.; Shan, J.; Gabbert, U. Inverse compensation of hysteresis using Krasnoselskii-Pokrovskii model. IEEE/ASME Trans. Mechatron. 2018, 23, 966–971. [Google Scholar] [CrossRef]
  26. Liu, Y.; Xie, S.; Du, D.; Qi, N. A Finite-Memory Discretization Algorithm for the Distributed Parameter Maxwell-Slip Model. IEEE/ASME Trans. Mechatron. 2020, 25, 1138–1142. [Google Scholar] [CrossRef]
  27. Al Janaideh, M.; Rakheja, S.; Su, C.Y. An analytical generalized Prandtl–Ishlinskii model inversion for hysteresis compensation in micropositioning control. IEEE/ASME Trans. Mechatron. 2010, 16, 734–744. [Google Scholar] [CrossRef]
  28. Yang, M.J.; Gu, G.Y.; Zhu, L.M. Parameter identification of the generalized Prandtl–Ishlinskii model for piezoelectric actuators using modified particle swarm optimization. Sens. Actuators A Phys. 2013, 189, 254–265. [Google Scholar] [CrossRef]
  29. Ang, W.T.; Khosla, P.K.; Riviere, C.N. Feedforward controller with inverse rate-dependent model for piezoelectric actuators in trajectory-tracking applications. IEEE/ASME Trans. Mechatron. 2007, 12, 134–142. [Google Scholar] [CrossRef] [Green Version]
  30. Tan, U.X.; Latt, W.T.; Shee, C.Y.; Riviere, C.N.; Ang, W.T. Feedforward controller of ill-conditioned hysteresis using singularity-free Prandtl–Ishlinskii model. IEEE/ASME Trans. Mechatron. 2009, 14, 598–605. [Google Scholar]
  31. Al Janaideh, M.; Mao, J.; Rakheja, S.; Xie, W.; Su, C.Y. Generalized Prandtl–Ishlinskii hysteresis model: Hysteresis modeling and its inverse for compensation in smart actuators. In Proceedings of the 2008 47th IEEE Conference on Decision and Control, Cancun, Mexico, 9–11 December 2008; pp. 5182–5187. [Google Scholar]
  32. Kuhnen, K. Modeling, identification and compensation of complex hysteretic nonlinearities: A modified Prandtl–Ishlinskii approach. Eur. J. Control 2003, 9, 407–418. [Google Scholar] [CrossRef]
  33. Gu, G.Y.; Yang, M.J.; Zhu, L.M. Real-time inverse hysteresis compensation of piezoelectric actuators with a modified Prandtl–Ishlinskii model. Rev. Sci. Instrum. 2012, 83, 065106. [Google Scholar] [CrossRef]
  34. Rakotondrabe, M. Classical Prandtl–Ishlinskii modeling and inverse multiplicative structure to compensate hysteresis in piezoactuators. In Proceedings of the 2012 American Control Conference (ACC), Montreal, QC, Canada, 27–29 June 2012; pp. 1646–1651. [Google Scholar]
Figure 1. Working principles of the PI hysteresis model.
Figure 1. Working principles of the PI hysteresis model.
Micromachines 12 00942 g001
Figure 2. Example of (shifted) hysteresis obtained with three elementary operators.
Figure 2. Example of (shifted) hysteresis obtained with three elementary operators.
Micromachines 12 00942 g002
Figure 3. The paired properties of the digitized representation with two operators.
Figure 3. The paired properties of the digitized representation with two operators.
Micromachines 12 00942 g003
Figure 4. A digitized representation of the classical PI hysteresis model with 3 operators.
Figure 4. A digitized representation of the classical PI hysteresis model with 3 operators.
Micromachines 12 00942 g004
Figure 5. Paired properties of the dead-zone operators for extending the classical PI hysteresis model for asymmetric hysteresis behavior.
Figure 5. Paired properties of the dead-zone operators for extending the classical PI hysteresis model for asymmetric hysteresis behavior.
Micromachines 12 00942 g005
Figure 6. Block diagram of the inverse model as a feedforward controller.
Figure 6. Block diagram of the inverse model as a feedforward controller.
Micromachines 12 00942 g006
Figure 7. Experimental architecture.
Figure 7. Experimental architecture.
Micromachines 12 00942 g007
Figure 8. Experimental setup.
Figure 8. Experimental setup.
Micromachines 12 00942 g008
Figure 9. Modeling results using the classical representation and the digitized representation with dead-zone operators at 1 Hz, sinusoidal wave.
Figure 9. Modeling results using the classical representation and the digitized representation with dead-zone operators at 1 Hz, sinusoidal wave.
Micromachines 12 00942 g009
Figure 10. Compensation results at a changing amplitude sinusoidal wave.
Figure 10. Compensation results at a changing amplitude sinusoidal wave.
Micromachines 12 00942 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, C.; Feng, C.; Aye, Y.N.; Ang, W.T. A Digitized Representation of the Modified Prandtl–Ishlinskii Hysteresis Model for Modeling and Compensating Piezoelectric Actuator Hysteresis. Micromachines 2021, 12, 942. https://0-doi-org.brum.beds.ac.uk/10.3390/mi12080942

AMA Style

Zhou C, Feng C, Aye YN, Ang WT. A Digitized Representation of the Modified Prandtl–Ishlinskii Hysteresis Model for Modeling and Compensating Piezoelectric Actuator Hysteresis. Micromachines. 2021; 12(8):942. https://0-doi-org.brum.beds.ac.uk/10.3390/mi12080942

Chicago/Turabian Style

Zhou, Chao, Chen Feng, Yan Naing Aye, and Wei Tech Ang. 2021. "A Digitized Representation of the Modified Prandtl–Ishlinskii Hysteresis Model for Modeling and Compensating Piezoelectric Actuator Hysteresis" Micromachines 12, no. 8: 942. https://0-doi-org.brum.beds.ac.uk/10.3390/mi12080942

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop