# A Memory-Based Hysteresis Model in Piezoelectric Actuators

Authors: Guilin Zhang; Chengjin Zhang; Jason Gu
Journal: Journal of Control Science and Engineering
 (2012)

---

## Abstract

A mathematical memory-based model is proposed to capture the hysteresis behavior in piezoelectric actuators. It is observed that the ascending (descending) hysteresis curves are alike and converge to one point without memory saturation. Therefore, two, dominant curves are determined and expressed as continuous functions, and the other hysteresis curves are modeled using two dominant curves through nonlinear transforming of coordinate axis. In the event of memory saturation, a new converging point is used to compensate the model prediction error. The experimental study has been carried out and our proposed model prediction method is compared with PI model and the linear model. It shows that the proposed model prediction method is better than other two methods.


---

## Body

## 1. Introduction

Piezoelectric ceramics are widely used as actuators in nano-/micropositioning mechatronic systems due to their fast frequency response, nanometer scale resolution, and high stiffness. Since the materials of piezoelectric actuators are ferroelectric, nonlinear hysteresis behavior is commonly observed in such actuators in response to an applied electric field. Thus, it is very challenging to design high-performance servo controller for nano-/micropositioning mechatronic systems [1].In the past few years, many studies on the compensation for the positioning of the piezoelectric actuators have been reported. One simple useful method is the model-based feed-forward control scheme. This method cascades an inverse hysteresis model in series with an actuator plant to cancel out the effect of nonlinearity and achieve a relatively linear response. The more precision the inverse hysteresis model is, the better the control scheme tracking performance is. Consequently, modeling hysteresis is the first step to track control of piezoelectric actuators.Numerous research works have been done to model the hysteresis nonlinearity of piezoelectric actuators. The most well-known hysteresis model is the Preisach model [2–11]. It is the broadly used approach in modeling and compensation of hysteresis in piezoelectric actuators. However, Preisach model does not have an analytical inverse. Ge and Jouaneh improved classical Preisach model that can predict the hysteresis response of a piezoelectric actuator driven by a periodic sinusoidal or triangular input signal [5] and developed a computer-based tracking control approach for piezoelectric actuators based on linearizing the hysteresis nonlinearity [6]. Tan et al. formulated and proposed the value inversion algorithm using a class of discretized Preisach operators [9]. They proposed an adaptive identification of hysteresis in smart materials and developed an adaptive inverse control scheme to update the weight of Preisach operators [10].Another popular hysteresis model is Prandtl-Ishlinskii (PI) model, which is a subclass of the Preisach model. The advantage of PI model is that its inverse is also a PI type with different threshold and weighting values. Kuhnen developed an inverse scheme based on the inverse scheme of PI model [12, 13]. However, the PI operator in PI model has a symmetry property around the center of the loop, while hysteresis response of a piezoelectric actuator is not symmetric in practice. Bashash and Jalili investigated the nature of hysteresis in piezoelectric materials and discussed a modified PI model that describes the asymmetric and residual displacement properties simultaneously [14]. Jiang et al. proposed another modified PI model based on two asymmetric operators to model the asymmetric property of hysteresis [15]. The other disadvantage of PI model is that it cannot compensate for the saturated hysteresis. Al Janaideh Mohammad et al. proposed a generic PI model based on a generalized play operator with different loading and unloading envelop functions in order to characterize asymmetric and saturated hysteresis nonlinearities [16]. They also proposed the modified generalized PI model that has the exact analytical inverse [17].Xu and Wong applied least squares support vector machine (LS-SVM) to capture the rate-independent and rate-dependent hysteresis nonlinearities [18, 19]. Moreover, a set of models are proposed to describe the hysteresis properties including targeting turning points, curve alignment and wiping-out property. Tzen et al. applied an exponential curve to fit the hysteresis path when the piezoelectric actuator operated away from the saturation [20]. Sun et al. proposed a hysteresis model based on similarities of the hysteresis curves and the turning points [21]. Bashash and Jalili disclosed the memory-dominant nature of hysteresis in piezoelectric materials and developed a memory-based mathematical hysteresis model [22–24]. They adopted the proposed model in an inverse model-based control scheme for feedforward compensation of hysteresis nonlinearity [22].In this paper, we adopt the mathematical transformation to describe the hysteresis curve functions. The ascending and descending loading curves can be determined and expressed as two continuous functions. Then the rest of hysteresis curves adopt their shape from the dominant curves. In addition, a new converging is updated to compensate the model prediction error in the event of memory saturation. It is demonstrated that the proposed model prediction is noticeably improved compared with the model without memory operator.The remainder of this paper is organized as follows. First, a mathematical memory-based hysteresis model in piezoelectric actuators is proposed in Section2. Moreover, the approximate congruency property of the proposed model is validated in this section. Then the experimental setup and the proposed model performances are given in Section 3. In Section 4, a new converging point is updated to address the memory saturation. Finally, some concluding remarks are provided in Section 5.


## 2. Hysteresis Model

Hysteresis trajectory starts moving on the ascending loading curve as depicted in Figure1. This curve can be approximated by a monotonically increasing continuous function. As shown in Figure 1, the ascending curves are alike and converge to the upper targeting point. When the direction of the input changes, the trajectory breaks its path and moves downward on the descending loading curve that can be approximated by another monotonically increasing continuous function. All the descending curves are also alike and converge to the lower converging point. The ascending and descending dominant curves can be expressed as two monotonically continuous functions, fra(v) and frd(v), respectively. In addition, the rest of hysteresis curves can adopt their shape from these dominant curves.Figure 1
Ascending and descending dominant curves.For any ascending trajectory starting from point(v1,d1) and (v2,d2), the following linear function is obtained:(1)da(v)=kafra(v)+ba,
where da(v) represents the ascending hysteresis trajectory between points (v1,d1) and (v2,d2), when the input voltage v varies between v1 and v2, ka and ba are given by
(2)ka=du-d1fra(vu)-fra(v0),ba=d1-k*fra(v0),
where v0 and vu are the initial voltage and the upper converging voltage of the dominant curve depicted in Figure 1.Replacingka and ba from (2) into (1) yields
(3)da(v,v1,d1,v2,d2)=d1+du-d1fra(vu)-fra(v0)×(fra(v)-fra(v0)).
Once the initial voltage v1 and displacement d1 are decided, the hysteresis curves are obtained by using the dominant curve. However, when the initial voltage v1 substitute v in (3),
(4)da(v1)=d1+du-d1fra(vu)-fra(v0)(fra(v1)-fra(v0))≠d1.
Then, the ascending curve function is modified as
(5)da(v)=d1+du-d1fra(vu)-fra(v0)(fra(mv+(1-m)vu)-fra(v0)),
where m=(vu-v0)/(vu-v1) and when v=v1, da(v1)=d1 and when v=vu,da(vu)=du. It demonstrates that the model can start from the initial displacement and reach the upper converging point.Similarly, when the input voltagev varies between v1 and v2, the descending curve function is given by
(6)dd(v)=kdfrd(nv+(1-n))+bd,
where dd(v) represents the descending hysteresis trajectory between points (v1,d1) and (v2,d2), kd and bd are the same as ka and ba and given by
(7)kd=dl-d1frd(vl)-frd(vu),bd=d1-kd*frd(vu),n=vl-vuvl-v1,
where (vu,du) is the upper converging point and vu is also the initial voltage of the descending dominant curves, (vl,dl) is the lower converging point and can be identified. Then (6) can be transformed into
(8)dd(v,v1,d1,v2,d2)=d1+dl-d1frd(vl)-frd(vu)(frd(nv+(1-n)vl)-frd(vu)).
The congruency property states that any two minor hysteresis loops are identical and have the same shape if they are generated by back-and-forth variations of the input between two identical extrema. Here, we consider the case of two points (vx,dx1) and (vx,dx2)(dx1≠dx2). When the input voltage reaches vy, the displacement dy1 and dy2 are obtained by using (5). 						(9)dy1=dx1+du-dx1fra(vu)-fra(v0)(fra(mv+(1-m)vl)-fra(v0)),(10)dy2=dx2+du-dx2fra(vu)-fra(v0)(fra(mv+(1-m)vl)-fra(v0)).
Subtracting (9) from (10) yields
(11)dy2-dy1=dx2+du-dx2fra(vu)-fra(v0)×(fra(mvy+(1-m)vu)-fra(v0))-dx1+du-dx1fra(vu)-fra(v0)×(fra(mvy+(1-m)vu)-fra(v0))=(dx2-dx1)(1+fra(mvy+(1-m)vu)-fra(v0)fra(vu)-fra(v0))=(dx2-dx1)(1+σ).Whenvy≪vu, σ≈0, then dy2-dy1≈(dx2-dx1). That is to say, the congruency property of the proposed model can be approximately satisfied only under the condition that the input voltage is far away from the upper converging point.The implementation of the model needs to identify the ascending and descending dominant curves. Some structures have been proposed for these curves, including second-order polynomials, third-order polynomials, and exponential functions. Here, two-third-order polynomials are used for the approximation of the ascending and descending dominant curves. The polynomials are expressed as(12)fra(v)=∑i=03aiv3-i=a0v3+a1v2+a2v+a3,frd(v)=∑i=03div3-i=d0v3+d1v2+d2v+d3,
where ai and di are the dominant curve coefficients that can be identified by utilizing the least square identification method.The proposed model has several advantages compared to the Preisach model. Firstly, it only needs eight parameters compared with the Preisach model which needs a large number of parameters. Secondly, the numerical inversion of the proposed model can be easily obtained in real time based on the mathematical formulation.


## 3. Experimental Setup and Model Verification

To investigate the effectiveness of the proposed piezoelectric actuator hysteresis model, a set of experiment on a PST150/7/40VS12 PZT-driven is demonstrated using the PC-based control system. Figure2(a) shows the experimental block diagram. The setup includes a personal computer, a power controller, and a piezoelectric actuator, which provides maximum 40 μm displacement and an integrated high-resolution strain gauge position (SGS) sensor. The host computer generates the control codes, which are written in Visual C++ and run through the EPP (Enhanced Parallel Port). The digital signal is converted by 16-bit D/A converter and amplified by a high-voltage amplifier (PVC-150S1). The actual actuator output displacements are measured by a SGS sensor and converted to a digital signal by a 16-bit A/D converter. Then the displacement data are sent to the host computer. The experimental setup (without PC) is shown in Figure 2(b).(a) The experimental block diagram. (b) the experimental setup.
(a)(b)The dominant curves should be obtained first. A first-order reversal signal is applied to excite the piezoelectric actuator, which is shown in Figure3(a). The proposed model is applied to predict the output displacement. Figure 3(b) depicts the experimental displacement and the proposed model prediction. The maximum and root-mean-square (rms) error percentages of the proposed model are 0.96% and 0.49%, respectively.(a) The first-order reversal signal. (b) The experimental response and the proposed model response. (c) The experimental response and the linear model response, and (d) the experimental response and the PI model response.
(a)(b)(c)(d)To demonstrate the effectiveness of the proposed model, the linear model and PI models are considered for comparison. Both the linear and PI model are simulated based on the same input given in Figure3(a). The responses of these models are shown in Figures 3(c) and 3(d), respectively. The maximum and rms error percentages are 4.96% and 4.10% for linear model, and 1.32% and 0.52% for PI model, respectively. The modeling performance of the proposed model has been improved compared with the linear model and PI model.


## 4. Saturated Memory Model

The upper converging point is not a fixed point. In this section, the memory saturation is considered. The point 2 is a converging point as shown in Figure4. Path 1-2 is the dominant ascending curve. When the input signal overtakes point 2, the new hysteresis trajectory will follow the path 2–5 that is the extended path 1-2. Meanwhile, point 5 becomes the new upper converging point after the input voltage reaches point 5. Consequently, when the input signal exceeds a converging point, the new hysteresis trajectory will converge to a new point.Figure 4
Hysteresis curve movement path.We consider the case that theith ascending trajectory starts from point (vi1,di1) and (vu1,du1),vu1>vu. Because the converging point has varied in the next ascending curves, then the next ascending curves in (5) can be updated to (13) as follows:
(13)da(v)=d1+du1-d1fra(vu1)-fra(v0)×(fra(mv+(1-m)vu1)-fra(v0)).Different from the ascending curves, the descending curves converge to a fixed point(vl,dl). Thus, the descending curve functions keep unchanged in the event of memory saturation.To demonstrate the effectiveness of the saturated memory model, a simulation study is carried out. The input profile shown in Figure5(a) is applied to excite the piezoelectric actuator. The experimental and the model responses without memory operator are demonstrated in Figure 5(b). As seen from Figure 5(b), the proposed model response diverges from the correct trajectory when it hits the dominant maximum input voltage. Then the saturated memory model response is given in Figure 5(c), which demonstrates the performance enhancement of the saturated memory model.(a) The experimental input signal. (b) The experimental response and the proposed model response, and (c) the experimental response and the saturated memory model response.
(a)(b)(c)


## 5. Conclusions

Hysteresis is the main nonlinearity in piezoelectric actuators. Based on the similarities of the hysteresis trajectory, a new mathematical hysteresis model was proposed to capture the hysteresis behavior in this paper. Utilizing the mathematical transformation, a memory-based modeling framework was developed and experimentally validated on a piezoelectric actuator. Moreover, a new converging point was updated to compensate the prediction error when the hysteresis path hits the upper converging point. The experimental and simulation results demonstrate the effectiveness of the saturated memory model.



---

*Source: 498590-2012-05-30.xml*