Research on Hysteresis of Piezoceramic Actuator Based on the Duhem Model

To improve the modeling accuracy of piezoceramic actuator in the precision positioning system, the Duhem hysteretic model of the piezoceramic actuator was proposed. The paper used the polynomial function to approach the piecewise continuous function and f(v) and g(v) in the Duhem model, adopted recursive least squares algorithm and gradient correction algorithm to identify parameter α, polynomial coefficients of f and g in the Duhem model, and established the nonlinear parametric model of the piezoceramic actuator. Contrasting the simulation results of recursive least squares algorithm and gradient correction algorithm, the modeling accuracy is 0.24% when adopting the recursive least squares algorithm, and the modeling accuracy is 0.11% when adopting the gradient correction method. The result showed that the gradient correction algorithm could meet the modeling accuracy better, and the structure of the algorithm is simple, adaptable, and easy to implement.


Introduction
The piezoceramic actuator is a kind of ideal drive elements in the microdisplacement technology currently and has the advantage of high positioning accuracy, large driving force, and fast response speed. Since the hysteresis, nonlinear, and creep resistance, inherent the piezoceramic actuator, the repeatability and accuracy of microdisplacement mechanism decreased, and the transient response speed becomes slower. [1][2][3][4][5][6][7][8] To decrease the impact of the nonlinearity and obtain better performance of the piezoceramic actuator, many researchers carried on an investigation in modeling and controlling of hysteresis nonlinear system. [9][10][11][12][13][14][15] Al Janaideh et al. [16] analyzed the inverse model of generalized Prandtl-Ishlinskii (PI) model to compensate the hysteresis nonlinearities of smart actuators. The experimental results verified that the inverse model of generalized PI model could be conveniently applied as a feedforward compensator and saturated hysteresis in magnetostrictive and SMA actuators. Dong and Tan [11] proposed a modified PI modeling method for rate-independent hysteresis in piezoelectric actuators and introduced a generalized backlash operator as the elementary operator into the model. By applying the proposed method, the hysteresis in the piezoelectric actuators and ultrasonic motor are illustrated, respectively. Hamimid et al. [17] proposed the minor hysteresis loops model based on parameters scaling of the modified Jiles-Atherton model by using judicious expressions. The proposed model had been applied for 3.2% Fe-Si nonoriented magnetic sheet, and the result showed the expected behavior when the flux density level was low. Chwastek [18] used the dynamic Takacs model to describe hysteresis loops in a thick nonoriented steel sheet and achieved the dynamic extension by using an additional component of the effective field. The experiment obtained a satisfactory agreement between the measured and the modeled hysteresis loops. Baghel and Kulkarni [19] proposed a hybrid technique to solve the parameter identification problem based on the Jiles-Atherton hysteresis model. And the proposed technique is flexible enough to incorporate improved GA and LM algorithms. Zirka et al. [20] analyzed the physical assumptions under the static and dynamic Jiles-Atherton (JA) hysteresis models. This led to the using in the model of a misleading entity resembling the coenergy instead of the actual energy. It is necessary to take measures to avoid this nonphysical feature; therefore the JA static model is reserved for applying in circuit simulators.
The Duhem model is a kind of differential hysteresis model, proposed by Duhem and Stefanini in 1897 [21]. The Scientific World Journal The Duhem model has the explicit function expression and is the function of input signal derivative; the output of model is related to the rate of input signal. The model is a kind of dynamic models [10], conforms the dynamic characteristic of hysteresis nonlinear in the actual intelligent materials, and could describe the hysteresis nonlinear precisely. However it is difficult to obtain the parameter , coefficients of and in the Duhem model, and it would be the obstacles for application of Duhem model.

Parameters Identification of the Hysteresis Model
The Duhem has an explicit differential expression, through adjusting the parameter , coefficients of and of the Duhem model; different hysteresis characteristics could be reflected; while identifying the parameters of the Duhem model accurately, the hysteresis model of the piezoceramic actuator could be obtained [2]. The differential function of the Duhem model is where is constant, V is the input voltage, is output displacement, (V) and (V) are piecewise continuous functions.
[ , ] represents the set of all continuous functions defined in the closed interval [ , ], to the arbitrary represents the distance between them [2]. Letting ∈ [ , ], to the arbitrary given > 0, the polynomial existed and the following equation holds: To the arbitrary given (V) ∈ [ , ] and approximation precision, there has an algebraic polynomial where is natural number and ‖ − ℎ‖ ∞ ≤ . When accuracy > 0, the order of (V) and (V) is and , respectively; the polynomials are as follows: And (5) could be transformed into As the input voltage V, output displacement , and V/ , / are measurable, while identifying the parameters , , and accurately, the parameterized model of the Duhem model could be obtained. Letting where V( ) is the input voltage of the system at time , ( ) is the output displacement at time . Letting ( ) = ( ) × , ( ) is the data vectors of the input voltage, is the identified parameter vector. That is, Let that is, The target of the parameter identification is to obtain the value of the parameter when the function is the minimum.
Applying the recursive least squares algorithm to recursive equations (11), (12), and (13), the identification parameters are Equation (11) where ( ) is the weight matrix, the effect of the weight is to control the influence of the input component.

Parameters Identification Simulation and Modeling of the Duhem Model
To verify the accuracy of the algorithms for the Duhem model identification, the paper applied recursive least squares algorithm and gradient correction algorithm to identify the parameters of the Duhem model based on MATLAB simulation software, respectively, and contrasted the influence of the two identification algorithms to the model modeling accuracy.

Identification of Recursive Least Squares Algorithm.
In the experiment, the order of the polynomial (V), (V) is = 3, = 2, respectively; that is,  Utilizing the parameter identification data, the hysteresis curve of the model is shown in Figure 2. The red and blue curves represent input-output hysteresis curve of the Duhem model and actual input-output hysteresis curve, respectively. Figure 2 showed that the output of the Duhem model and the actual output data are basically consistent. Figure 3 is the relative error curve between the model output and actual output. It can be seen that the maximum error is 0.066 m. The result of the experiment verified the validity of the recursive least squares algorithm.

Gradient Correction Algorithm.
The data of the input and output is shown in Figure 1; under the gradient correction algorithm, Λ ( ) is as follows: 4 The Scientific World Journal  The parameter identification result is shown in Figure 4. It can be seen form Figure 4 that the identification parameters tend to be stable when recursiving to = 6; the parameter identification results are as follows: Utilizing the gradient correction parameter identification results, the model hysteresis curve is shown in Figure 5, The red and blue curves represent input-output hysteresis curve of the Duhem model and real input-output hysteresis curve, respectively. The error curve between the system output and model output is shown in Figure 6. It can be seen from Figure 6 that the maximum is 0.048 m. The result also verified the validity of the gradient correction algorithm.
The identification parameters of the two algorithms are shown in Table 1. And we show part of the relative errors contrast results under the two algorithms in Table 2. It can be seen from Table 1 that the relative errors between the real output and model output under the recursive least squares algorithm could reach 0.24%, the mean square deviation of the error is 0.0263, and the maximum error is 0.066 m; in contrast, the relative errors between the actual output and model output under the gradient correction algorithm could reach 0.11%, the mean square deviation of the error is 0.0222, and the maximum error is 0.048 m.
The Scientific World Journal 5

Conclusion
The paper utilized the polynomial to approach the piecewise continuous functions and of the Duhem model, adopted the recursive least squares and gradient correction algorithm, respectively, to identify the parameter , coefficients of and of the Duhem model, and applied the identified parameters to model the Duhem model. The experiment results showed that the modeling accuracy of the recursive least squares algorithm could reach 0.24%, the mean square deviation of the error is 0.0263; the modeling accuracy of the gradient correction algorithm could reach 0.11%, the mean square deviation of the error is 0.0222. The results of the experiment certified validity of the recursive least squares algorithm and gradient correction algorithm. Contrasting with the least squares algorithm, the gradient correction algorithm is adaptable and suitable for engineering. Applying the gradient correction algorithm, the Duhem model could be established more precisely and lay the foundation for the further control research of the piezoceramic.