# Modeling Hysteresis with Inertial-Dependent Prandtl-Ishlinskii Model in Wide-Band Frequency-Operated Piezoelectric Actuator

Authors: Vahid Hassani; Tegoeh Tjahjowidodo; Albert D. Soetarto
Journal: Smart Materials Research
 (2012)

---

## Abstract

One of the major problems occurring in many technical applications is the presence of the hysteretic behavior in sensors and actuators, which causes a nonlinear relationship between input and output variables in such devices. Since the nonlinear phenomenon of hysteresis degrades the performance of the piezoelectric materials and piezoelectric drive mechanisms, for example, in positioning control framework, it has to be characterized in order to mitigate the effect of the nonlinearity in the devices. This paper is aimed to characterize and model the hysteresis in typical piezoelectric actuators under load-free and preloaded circumstances incorporating the inertial effect of the system. For this purpose, the piezoelectric actuator is modeled as a mass-spring-damper system, which is expressed in terms of a stop operator as one of the essential yet efficient hysteresis operators in the Prandtl-Ishlinskii (PI) model. The reason of utilizing the stop operator in this study is for the sake of control purposes, as the stop operator plays as the inverse of the play operator in the PI model and can be used in a feed-forward controller scheme to suppress the effect of hysteresis in general control framework. The results reveal that this model exhibits better correspondence to the measurement output compared to that of the classical PI model.


---

## Body

## 1. Introduction

Hysteresis is a nonlinear phenomenon that occurs in some types of materials such as piezoceramics, shape-memory alloys, and magnetostrictive actuators. The word “hysteresis” refers to systems that have memory such that similar loops are repeated in each cycle of operation in every dynamical system and, after removal of the force or electrical field, the system does not return to its original location. In this phenomenon, the current output depends also on the history of the input.The effect of hysteresis nonlinearity can be neglected in some systems. In contrast, ignoring this phenomenon in other types of systems that possess severe hysteresis might create undesirable consequences such as inaccuracy in open loop system, limit cycle, degradation of the tracking performance, and even instability of closed loop system.To overcome this challenge, some mathematical models, for example, the Preisach model and the PI model, have been proposed to capture the effect of hysteresis in any mechanical systems. Utilization of the PI model is addressed in many works to characterize the hysteresis due to its simplicity compared to the Preisach model. Wang et al. [1] utilized the classical PI model and neural network approximator to design a robust adaptive control, which ensures that the tracking error is bounded in the presence of unknown hysteresis in a closed-loop system. The generalized PI is proposed for modeling an asymmetric hysteresis in piezoceramic actuators, where a linear function is defined to modify the threshold of the play operator of the model [2]. Krejci and Kuhnen [3] presented a new method of constructing the inverse of the hysteresis model in an analytical form for the system with hysteresis and creep and used the proposed model in a feed-forward controller to mitigate the nonlinearity arisen from the hysteresis. Al Janaideh et al. [4] presented a method to construct the inverse of the generalized PI model to be used in a feed-forward loop cascaded with an adaptive controller to reduce the effect of hysteresis in a piezoelectric actuator, and in another paper, they utilized the generalized PI model to design a robust sliding mode controller to ensure the system stability for smart actuators [5]. A new strategy is introduced by Chen [6] to fuse the PI model with adaptive control techniques, which enhances the global stability of the system.Various rate-dependent PI models have been proposed to capture the hysteresis in different rate of input or different frequency values of input applied to the piezoelectric actuators. Al Janaideh et al. [7] proposed a new rate-dependent PI model, which utilized logarithmic threshold value and density function in the play operator to characterize the hysteresis of a piezoelectric actuator under dynamically varying input up to 200 Hz. In the paper, they also characterized the hysteresis of a piezoelectric actuator up to 500 Hz and showed that the hysteresis loop gets larger at higher frequencies as well as it generates lower displacement as the frequency is increasing [8]. Kuhnen and Janocha [9] developed a modified PI model to characterize the complex hysteresis in actuators and sensors and presented the inverse of this model in an analytical form. Ang et al. [10] proposed the hysteresis model at low frequencies up to 50 Hz and utilized a linear function relating the slope of the hysteresis loop to input rate applied to the actuator. Subsequently, an advanced model was proposed to characterize the nonlinearity in a piezoelectric actuator under ill-conditioned hysteresis [11]. An adaptive controller was designed to capture the nonlinear effects of the giant magnetostrictive actuator using the PI model to compensate for the hysteresis which was reported by Wang et al. [12]. Goldfarb and Celanovic [13] proposed a rate-independent generalized Maxwell-Slip model for force and positioning control of microrobotic systems. Quant et al. [14] proposed a nonlinear Maxwell-Slip-based model to replace the linear constitutive equations of piezoelectric actuators to make a comparison between two models. A Maxwell-Slip-based model was proposed to capture the hysteresis in piezoelectric actuator, and the parameters of this model were identified by the linear programming method [15]. Three fitness functions were proposed in terms of the Bouc-Wen model including mass-spring-damper characteristics of piezoelectric actuator, and these functions were minimized using the genetic algorithm method [16]. In other study, piezoelectric actuator was employed as Impact Drive Mechanism, and its model was constructed assuming the system as a spring-mass-damper system, where nonlinearity was modeled using the Bouc-Wen model, and finally the parameters of this system were identified by an adaptive method [17].In this study, an inertial-dependent model is proposed to overcome the large amount of errors imposed in the modeling with the classical PI model, in particular in capturing the hysteresis at high frequencies as well as high amplitudes of inputs applied to a piezoelectric actuator under load-free and preloaded conditions. As a rate-independent model, it is able to capture the hysteresis at different input magnitudes and frequencies even in the presence of preloading force.


## 2. Classical Prandtl-Ishlinskii (PI) Model

One of the well-known operator-based models, which are used to characterize the hysteresis in piezoceramics and other types of materials, is the classical PI model that utilizes two essential hysteresis operators, namely, the play operator and the stop operator.


### 2.1. Stop Operator

The rate-independent stop operator is illustrated in Figure1. Conceptually, this operator determines the linear stress-strain (w-v) relationship as in the Hook’s law, when the stress is less than the yield stress, r [7]. In analytical form, assume that Cm[0,tE] is the space of piecewise monotone continuous functions, a function er:R→R is defined by(1)er(v)=min⁡(r,max⁡(-r,v))
with r≥0. For any input v(t)∈Cm[0,tE] and any initial value ω-1∈R, the stop operator Er can be expressed by the following definition:(2)Er[v;ω-1](0)=er(v(0)-ω-1),Er[v;ω-1](t)=er(v(t)-v(ti)+Er[v;ω-1](ti)),forti<t≤ti+1,0<i≤N-1,
where 0=t0<t1<⋯<tN=tE is a part of [0,tE], such that the function v is monotone on each of the subintervals [ti,ti+1].Figure 1
Classical stop operator.


### 2.2. Play Operator

The rate-independent play operator is illustrated in Figure2. This operator implies the reciprocating motion of a piston inside a cylinder with length 2r, where the instantaneous position of center of piston is defined by v and ω as input and output subsequently [7]. In analytical form, for any input v(t)∈Cm[0,tE] and r≥0, the play operator is inductively defined by(3)Fr[v;ω-1](0)=fr(v(0)-ω-1),Fr[v;ω-1](t)=fr(v(t),Fr[v;ω-1](ti)),forti<t≤ti+1,0<i≤N-1,
where(4)fr(v,ω)=max⁡(v-r,min⁡(v,ω))ω-1 defines the initial state at time t=0.Figure 2
Classical play operator.


## 2.1. Stop Operator

The rate-independent stop operator is illustrated in Figure1. Conceptually, this operator determines the linear stress-strain (w-v) relationship as in the Hook’s law, when the stress is less than the yield stress, r [7]. In analytical form, assume that Cm[0,tE] is the space of piecewise monotone continuous functions, a function er:R→R is defined by(1)er(v)=min⁡(r,max⁡(-r,v))
with r≥0. For any input v(t)∈Cm[0,tE] and any initial value ω-1∈R, the stop operator Er can be expressed by the following definition:(2)Er[v;ω-1](0)=er(v(0)-ω-1),Er[v;ω-1](t)=er(v(t)-v(ti)+Er[v;ω-1](ti)),forti<t≤ti+1,0<i≤N-1,
where 0=t0<t1<⋯<tN=tE is a part of [0,tE], such that the function v is monotone on each of the subintervals [ti,ti+1].Figure 1
Classical stop operator.


## 2.2. Play Operator

The rate-independent play operator is illustrated in Figure2. This operator implies the reciprocating motion of a piston inside a cylinder with length 2r, where the instantaneous position of center of piston is defined by v and ω as input and output subsequently [7]. In analytical form, for any input v(t)∈Cm[0,tE] and r≥0, the play operator is inductively defined by(3)Fr[v;ω-1](0)=fr(v(0)-ω-1),Fr[v;ω-1](t)=fr(v(t),Fr[v;ω-1](ti)),forti<t≤ti+1,0<i≤N-1,
where(4)fr(v,ω)=max⁡(v-r,min⁡(v,ω))ω-1 defines the initial state at time t=0.Figure 2
Classical play operator.


## 3. Inertial-Dependent Prandtl-Ishlinskii (PI) Model

In order to determine the dynamic characteristics of the piezoelectric actuators, the linear constitutive equation has been proposed [18]. A piezoelectric element can be represented by an element with electrical and mechanical gates in linear form, which is expressed mathematically in (5) disregarding the presence of any nonlinear effects(5)T=CES-eE.
Since this equation is not involved in the nonlinear phenomenon of hysteresis, it is required to model the nonlinear characteristics of the piezoelectric actuator. For practical applications, engineering variables such as force, speed, electrical current, and voltage are required. In order to simplify the analysis, the strain, S, stress, T, and electrical field, E in (5) are substituted by measurable counterparts, namely, the displacement, x, force, F, and voltage, U, as shown in (6)F=A⋅U+Zx.
Consequently, the constitutive equation can be transformed into the Laplace domain as presented in (7), where Z represents the mechanical impedance, m is the equivalent mass of the piezoceramic, c is the damping factor, and k is the stiffness of the piezoelectric actuator(7)Z=ms2+cs+k.
Intuitively, in this model, the piezoelectric actuator is assumed as a mass-spring-damper system. Based on this assumption, basic equation of motion for such systems can be written incorporating with the PI model to take into account the both linear and nonlinear characteristics of the piezoelectric actuator. As a result, the following relationships are proposed in terms of two essential operators of the PI model, namely, the play operator:(8)X=α⋅V+β⋅FP⋅L-meqd2fP⋅O(V)dt2-ceqdfP⋅O(V)dt,
and inversely in the stop operator, it can be expressed as(9)fS⋅O(X)=β⋅FP⋅L-meqd2Xdt2-ceqdXdt-keqX,
where meq, ceq, and keq are the equivalent mass, damping, and stiffness factors, respectively. FP·L stands for the preloading force when the piezoelectric actuator is exposed to the preloading condition, which in turn involves two terms of external force and stiffness of the preloading devices which are expressed by FP·L=Fext+Keq,ext·x. In (8) and (9), V represents the input voltage applied to the piezoelectric actuator as the input of the play operator and X represents the desired displacement of the actuator as the input of the stop operator. α and β stand as correction coefficients of the input voltage in the play operator and the preloading force FP·L, respectively. fS·O(X) and fP·O(V) are subsequently the output of the stop operator and the play operator.This model takes into consideration two important terms of inertial and damping effects which are shown bymeqẌ and ceqẊ. In the sense of frequency response, these two terms are expressed by meqω2X and ceqωX. As a result, these terms are function of frequency, and they vary with variations of the input frequency. Due to this characteristic, the modeling of the hysteresis in piezoelectric actuators is expected to perform better compared to that with the classical PI model.


## 4. Experimental Set Up

A dedicated set up was built to characterize a piezoelectric actuator experimentally. In this work, the piezoceramic actuator (P-887.90, from PHYSIK Instrument Co.) was used. This actuator provides maximum displacement of36±10%μm at maximum applied voltage of 120 V. The PZT actuator is driven using a power amplifier (E.663, PHYSIK Instrument Co.), which multiplies any input voltage by gain 10 and sends it to the piezoelectric actuator. A capacitive sensor (CS1, Micro-Epsilon Co. with resolution of 0.75 nm and operating temperature range of −50°C to 200°C) was utilized to measure the displacement of PZT actuator in terms of voltage (see Figure 3).Figure 3
Configuration of the experimental setup.


## 5. Experimental Results

In order to characterize the hysteretic behavior of the piezoelectric actuator, two sets of sinusoidal inputs for low-amplitude voltage signal (20 V amplitude with 20 V bias) and high-amplitude voltage signal (40 V amplitude with 40 V bias) were applied to the piezoelectric actuator in load-free and preloading conditions. These sets of input signals were applied at varying frequency ranging from 10 Hz to 200 Hz and different preloading forces of 4 N, 8 N, and 12 N. The reason of selecting different amplitudes besides different frequencies is to analyze the hysteresis characteristic in the piezoelectric actuator in terms of amplitude and frequency variations. Different preload condition on the piezoelectric actuator is also expected to affect the hysteresis characteristics. This motivates the experiment to be conducted at different preloading forces.Applying two sets of low-amplitude and high-amplitude inputs at different frequencies besides different values of preloading force implies three important aspects in the hysteretic behavior. The first thing, which is observed from this experiment, is attributed to the rate dependency of the hysteresis curve; that is, when the rate of the input voltage is increased, the width of the hysteresis curve is also increasing. Figure4 shows that the maximum width of the hysteresis at frequency of 200 Hz with low input voltage has become 3.57 times as that of the hysteresis curve at frequency of 10 Hz with the same amount of voltage amplitude. For high amplitude voltage, the width of the hysteresis curve at frequency of 200 Hz has grown up to 4.3 times higher than that of the hysteresis width at frequency of 10 Hz. The results also show that the maximum width of hysteresis curve is expanded when the magnitude of the applied voltage is increased.Width and slope of the hysteresis curve at different frequencies and input amplitudes for load-free piezoelectric actuator.
(a)(b)Figure5 shows the hysteresis obtained from the experiments to characterize the properties of the piezoelectric actuator under preloading condition at three different preloading forces applied to the piezoceramic. From the experiment results using two input signals, the most important thing that can be concluded is that the expansion of the hysteresis curve corresponds to the preloading force. For the case with an input voltage of 40 V and frequency of 200 Hz, the width of the hysteresis is expanded 1.06 times from preloading force of 4 N to preloading force of 12 N. In the case with the input voltage of 80 V and frequency of 200 Hz, this value has shifted up to 1.11 times under preloading force of 4 N in comparison with the preloading force of 12 N. This characterization shows that the intensity of hysteresis is increasing with the preloading forces, and consequently the quantity of maximum displacement of the actuator is also decreasing with increasing preloading force.Width and slope of the hysteresis curve at different frequencies and input amplitudes for preloaded piezoelectric actuator.
(a)(b)(c)(d)(e)(f)The second finding of this experiment is regarding the relationship between the maximum stroke of the actuator and maximum input voltage applied to the actuator. For the case of 40 V input, the maximum stroke of 10.12μm is obtained at frequency of 10 Hz, and maximum stroke of 9.454 μm is obtained at frequency of 200 Hz. For the case of 80 V input, different values are acquired; meanwhile, the difference between the lowest and the highest stroke of the actuator at frequency of 10 Hz and 200 Hz is larger than that of the case in the 40 V input. These values are subsequently 21.57 μm and 13.81 μm for frequencies of 10 Hz and 200 Hz. The same phenomenon is also observed under the preloading condition.The third finding observed from this experiment is regarding a tight relationship between amplitude and frequency in one hand and asymmetrical hysteresis shape and high slope of the hysteresis curve in the other hand. As illustrated in Figures5 and 6, if the amplitude of the input voltage is increased, the asymmetric hysteresis curve will be formed at low frequencies as compared to application of the low amplitude input, and the slope of the hysteresis curve is also increasing. In the case of the input voltage with magnitude of 80 V, the asymmetric hysteresis curve is obtained at frequency of 125 Hz, while the asymmetric hysteresis is observed at 250 Hz for the lower input amplitude of 40 V.Asymmetric hysteresis curve at different frequencies and input amplitudes for load-free piezoelectric actuator.
(a)(b)


## 6. Modeling Results

In order to capture the symmetrical hysteresis shapes, which are recorded up to 200 Hz with low-amplitude and high-amplitude voltages, a stop operator was utilized. As described in Section2, a stop operator plays a role as the inverse of the play operator. The implementation of the stop operator instead of the play operator enables us to use the PI model in a feed-forward loop for positioning control purposes so that it can generate the required control input upon a desired displacement commanded to the control system.The classical PI model was implemented to characterize the hysteresis in the piezoelectric actuator by utilizing four classical stop operators, in the cases, where the symmetric hysteresis curve takes place. This modeling approach is parameterized by four threshold values ofr and multipliers of w, where r adjusts the width and height of the hysteresis and w adjusts the slope of the hysteresis curve correspondingly. These parameters are optimized upon the minimization of the following cost function:(10)J=∑i=1n(yS⋅O-ym)2,
where yS·O is the output of the stop operator in PI model and ym is the actual output which was measured from the experiment and n is the number of samples. The error minimization problem was solved using the MATLAB optimization toolbox.As can be seen in Figures7 and 8 the classical PI model is not capturing the hysteresis pattern satisfactorily, in particular at high frequencies for both cases of the unloaded and preloaded actuator.The results of modeling with the classical PI for load-free piezoelectric actuator for different input voltage and frequency: (a) 40 V; 10 Hz, (b) 80 V; 10 Hz, (c) 40 V; 25 Hz, (d) 80 V; 25 Hz, (e) 40 V; 200 Hz, (f) 80 V; 200 Hz.
(a)(b)(c)(d)(e)(f)The results of modeling with the classical PI for preloaded piezoelectric actuator under 12 N force, for different input voltage and frequency: (a) 40 V and 10 Hz, (b) 80 V and 10 Hz, (c) 40 V and 25 Hz, and (d) 80 V and 25 Hz.
(a)(b)(c)(d)In order to overcome the problem, as discussed in the previous section, the inertial-dependent model is introduced, where two important dynamic terms of inertial and damping effects, which are stated bymeqẌ and ceqẊ, are taken into account. In the sense of frequency response, these two terms are expressed by meqω2X and ceqωX. Intuitively, the physical aspect of the actuator will also manifest itself in the dynamics of the system. The classical PI model cannot be used solely in characterizing the hysteresis properties of the material. In Figures  7(a) to 7(f), it is shown that the classical PI model cannot capture the hysteresis pattern especially at the turning points, where the major lag occurs between the applied electrical field and actuator displacement. At these points, the PI model generates sharp turning points due to the absence of inertial and damping effects. In order to tackle this discrepancy, inertial-dependent model is implemented to capture the turning points as well as other points of the hysteresis curve with considerable reduction of error as illustrated in Figures  9(a)–9(d). In these figures, the results of simulation from the inertial-dependent PI model are displayed, where the amount of error between the measured output and the output of the inertial-dependent PI is significantly reduced compared to that from the classical PI model.The results of modeling with Inertial-dependent PI for load-free actuator with different input voltage and frequency: (a) 40 V and 10 Hz, (b) 80 V and 10 Hz, (c) 40 V and 200 Hz, and (d) 80 V and 200 Hz.
(a)(b)(c)(d)One of the interesting points appears from the characterization of the hysteresis at frequency of 200 Hz with high-amplitude voltage, where the hysteresis curve starts to skew resulting in an asymmetrical hysteresis. The error in this case seems to be larger compared to the cases at lower frequencies (see Figure9(d)), which implies that the asymmetric hysteresis is still not incorporated in this model, despite the considerable amount of error that has been reduced in the symmetric hysteresis cases.Figure10 shows the identification results using the inertial-dependent PI model for the piezoelectric actuator under preloading forces of 4 N, 8 N, and 12 N. The first three panels (Figures  10(a)–10(c)) represent the preloading cases for the first applied voltage signal (40 V), where the last three panels are for the case of the second voltage signal (80 V). It can be concluded that this model is also applicable to characterize the piezoelectric actuator under preloading condition with acceptable amount of error by using the inertial-dependent PI model with the consideration of the preloading force, FP·L=Fext+Keq,ext·x.The results of modeling with the inertial-dependent PI with different input voltage, frequency, and preload: (a) 40 V; 200 Hz; 4 N, (b) 40 V; 200 Hz; 8 N, (c) 40 V; 200; 12 N, (d) 80 V; 100 Hz; 4 N, (e) 80 V; 100 Hz; 8 N, and (f) 80 V; 100 Hz; 12 N.
(a)(b)(c)(d)(e)(f)The result shows that the introduction of the inertial and damping effect to the classical PI model improves the performance of the model in wide range of operation. As can be seen in Figures7 and 8, the error between the measured output (voltage) and the actual output is increased from lower frequencies to higher frequencies in the classical PI model. Also, the quantity of the error for high-amplitude voltage is larger than the error for the low-amplitude voltage case. In the inertial-dependent model, not only the amount of error remains low at different amplitudes, but this amount also remains low at different frequency values. As a result, it can be realized that the classical PI model is not capable of characterization of the hysteresis at wide range of frequency solely.On the other hand, from the simulation at different frequencies with different amplitudes in both preloading and load-free circumstances, the parameter of damping effect is varying as illustrated in Figure11. As the amplitude goes to higher values for both unloading and preloading conditions, the values of the damping factor shift to the higher values correspondingly, while the variation of the damping factor to the frequency is exponentially decreasing. On the other hand, in each set of the input voltage imposed to the actuator, the damping factor decreases with increasing the frequency value. This finding could be considered to be a guideline to develop an integrated rate and inertial-dependent model of hysteresis based on the variation of the damping factor versus frequency.Damping factor variations at different frequencies with different amplitudes (▲, low-amplitude voltage, 40 V; ●, high-amplitude voltage, 80 V).
(a)(b)(c)(d)


## 7. Conclusion

In this study, the performance of the classical PI model was evaluated using sets of data which incorporate two different magnitudes of the input voltage and varying frequency from 10 Hz to 200 Hz for different preloading conditions. It is shown that the classical Prandtl-Ishlinskii model is unable to characterize the hysteresis in piezoelectric actuators at wide range of frequency and amplitude.In order to tackle this problem, as the physical aspect of the piezoelectric actuator will manifest itself in the dynamic behavior, the inertial-dependent PI model is proposed to capture the hysteresis at different input amplitudes and frequencies. The proposed model exhibits better performance compared to that of the classical one as implied by the modeling error for wide range of applications. However, moving toward high-frequency applications, the performance of the proposed model is degrading. Incorporating the rate-dependent parameters, for example, as presented by Al-Janaideh et al. [8], to the inertial-dependent model is expected to improve the performance.



---

*Source: 164062-2012-01-24.xml*