Fuzzy PID Feedback Control of Piezoelectric Actuator with Feedforward Compensation

Piezoelectric actuator is widely used in the field of micro/nanopositioning. However, piezoelectric hysteresis introduces nonlinearity to the system, which is the major obstacle to achieve a precise positioning. In this paper, the Preisach model is employed to describe the hysteresis characteristic of piezoelectric actuator and an inverse Preisachmodel is developed to construct a feedforward controller. Considering that the analytical expression of inverse Preisach model is difficult to derive and not suitable for practical application, a digital inverse model is established based on the input and output data of a piezoelectric actuator. Moreover, to mitigate the compensation error of the feedforward control, a feedback control scheme is implemented using different types of control algorithms in terms of PID control, fuzzy control, and fuzzy PID control. Extensive simulation studies are carried out using the three kinds of control systems. Comparative investigation reveals that the fuzzy PID control system with feedforward compensation is capable of providing quicker response and better control accuracy than the other two ones. It provides a promising way of precision control for piezoelectric actuator.


Introduction
As is known, human enters the world of micro/nanolevel with the inventions of scanning tunneling microscopy (STM) [1] and atomic force microscope (AFM) [2].One key technology in STM and AFM is micro/nanopositioning.Actually, micro/nanopositioning has been applied in more and more fields nowadays.Regarding the drive principle in micro/ nanopositioning system, piezoelectric actuator is popular because of its high stiffness, fast response, and several other outstanding features.However, piezoelectric actuator introduces some obvious limitations, such as hysteresis, creep, and vibration characteristics.How to realize the precise control of piezoelectric actuator is a hot research topic in recent years.
Generally, under open-loop voltage control, the piezoelectric actuator produces 10%-15% error with respect to full range [3].With the increase of the input signal frequency, the error will even reach to 35% [4].So, the hysteresis characteristics of piezoelectric actuator are the main problem to be overcome.In the literature, a physical explanation for the hysteresis phenomenon from a macroscopic viewpoint was given by Chen and Montgomery [5].Yet, piezoelectric actuator exhibits more complex hysteresis nonlinearity [6].In particular, the output signal not only depends on the input signal, but also relates to the history of the system state.Thus, for the same input signal under different states, the output signal will be different.In addition, previous studies have shown that the frequency of the input signal also affects the output signal and error.
To realize the control of the piezoelectric actuator to cater for the requirement of micro/nano positioning, appropriate mathematical models can be established to characterize the piezoelectric hysteresis.Researchers have established various models from different perspectives to describe the hysteresis nonlinearity of piezoelectric actuator.As shown in Figure 1, the hysteresis models can be mainly classified into two types: physics-based models and phenomenological models.Physics-based models are used to describe the basic physical principle of material and the hysteresis models are obtained in view of the relations of energy, displacement, and so on [7][8][9].Alternatively, phenomenological models start from the characteristic of hysteresis curve.They are employed to describe the hysteresis curve by using the effective math model directly, without paying attention to the physical meaning [10][11][12][13].Specifically, Preisach model is a popular phenomenological hysteresis model, and it has been widely used.It is able to give an accurate description of the characteristics of the hysteresis nonlinearity.In this research, the Preisach model is used to describe the hysteresis nonlinearity of a piezoelectric actuator.
Concerning the control scheme of hysteresis nonlinear system, the inverse model of hysteresis nonlinearity is usually obtained first.Then, the inverse model is used to construct a feedforward controller to compensate for the hysteresis effect of the system.Using the inverse model as the feedforward compensation directly is a simple and effective method.The inverse compensation model and the piezoelectric actuator, which are connected in series, can be considered as a linear system.In order to achieve this goal, a lot of previous works have been conducted in the literature.To name a few, Leang and Devasia [14] adopted iterative learning control strategy to solve the Preisach model's inverse compensation control issue.Krejci and Kuhnen [15] derived the inverse analytical expression of traditional Prandtl-Ishlinskii (P-I) model and reduced the tracking error by one order of magnitude.Xu and Wong [16] built an inverse hysteresis model using support vector machines for compensating the hysteresis nonlinearity of piezoelectric actuator and then demonstrated that it is more effective than Bouc-Wen model and P-I model via experimental studies.
Control method based on inverse hysteresis model is simple and intuitive, but it has many drawbacks, such as heavy computational burden and complicated system structure.In particular, because of the complexity of the hysteresis model, finding out the analytic inverse model directly is difficult.Most of the time, numerical inverse model is used to approximate the exact model; thus it appears that the inverse model is not unique.Moreover, a standalone feedforward control is not sufficient to totally cancel out the positioning error because there always exist certain degrees of model error.Hence, a closed-loop feedback control can be designed to combine with the feedforward control in order to further mitigate the control error.Particularly, PID control is widely used because of its simple construction [17].In the literature, Tan et al. [18] proposed a learning type of PID controller and tried to enhance the robustness of the system.Additionally, intelligent controllers based on fuzzy logic and neural networks have been applied extensively in the control of piezoelectric actuator [19,20].The inverse Preisach model can also be used as a feedforward compensation which is added to PID feedback control [21].Moreover, this compound control method has also been applied in the joint angle control of a manipulator driven by pneumatic artificial muscles [22].In addition, Chen et al. [23] proposed a control method which combines the inverse Preisach compensation model with the indirect adaptive controller.An adaptive inverse model has also been proposed which is updated by least mean square algorithm [24].Recently, more different control methods for piezoelectric actuator have been proposed [25][26][27].The whole purpose of these control approaches is to achieve a precise and stable control.
As an important branch of intelligent control, fuzzy control is a control method on the basis of fuzzy set theory, fuzzy language variables, and fuzzy logic reasoning.It spans a wide application in various fields of control and automation.As a combination of fuzzy control and PID control, the fuzzy PID control is a popular control approach.Although both fuzzy control and fuzzy PID control have been widely used, it is unclear how fuzzy control performs in comparison with fuzzy PID control in piezoelectric actuator control.In this research, a comparison study of fuzzy control and fuzzy PID control with feedforward compensation is conducted for precision motion control of a piezoelectric actuator.Through a series of simulation comparative studies, some useful conclusions are derived.
The following parts of the paper are organized as follows.Section 2 gives a brief review of the Preisach model.Three kinds of controllers are then constructed in Section 3. Section 4 performs simulation studies of the three controllers.Some conclusions are drawn in Section 5.

Preisach Model
Preisach model was originally used to study the physical principle of magnetic hysteresis characteristics in phenomenon of  magnetization [10,28].After forty years of its development, the mathematicians Krasnosel'skii and Pokrovskii [29] separated the physical meaning of Preisach model in the 70s of the 20th century, gave a kind of pure mathematics characteristic model definition, and expanded the application field of the Preisach model.Nowadays, Preisach model has become one of the most widely used hysteresis models.

Model Expression. The mathematical description of classic Preisach model is shown as follows:
where () is the input of system; () is the output of system; (, ) is weighting function; ,  are the "rise threshold" and "fall threshold, " respectively; and   [()] is hysteresis operator with the value of +1 or −1.
Generally, Preisach model solves the current input response through the integration of historical input operations, and it has the characteristics of the global memory.In order to characterize the hysteresis of piezoelectric actuators, the Preisach model can be established as follows.
As shown in Figure 2, based on the rule of Preisach model  ⩾ , the integral area of (1) constructs the right triangle .The right triangle vertex is ( 0 ,  0 ), and the hypotenuse is the straight line  = .Any point (, ) within S corresponds to a hysteresis operator   .When   = 1, its  + part lies in area ; when   = 0, its  − part lies in area .(, ) is a distribution function which is defined within the triangle  and its value obeys the statistical law.In addition, (, ) = 0 lies in the area outside .
From Figure 3, the corresponding output of piezoelectric actuators can be calculated as (2)

Model Discretization.
It is found that although (2) can be used to calculate the output displacement, it is very difficult to solve.So, it is necessary to discretize this equation in order to facilitate its usage.
When the input () starts from 0 and increases to   , the output is (  ,   ).Then, () monotonically decreases to   , which produces an output (  ,   ).The change of the output is defined as (  ,   ): As shown, Figure 3(a) is the trajectory of input signal and Figure 3(b) is the domain of integration. + and  − are departed into  1 ,  2 , and  3 .From this, the final output can be calculated: Combining ( 4) with the definition in (3) yields When () is increasing, When () is decreasing, Through ( 6) and ( 7), the response of the output signal can be found out at any time.It is notable that only the nonmemory part needs to be considered to obtain an expression for the input signal based on this algorithm [30].
For illustration, a simulation result is shown in Figure 4.This figure clearly shows the Preisach curve after discretization.The upper one's simulation calculation time is  = 10; the lower one is  = 100.We can see that the discretization curve can be used to describe the hysteresis loop.Yes displacement as input, the compensator gives a compensation control signal for the piezoelectric actuator.This can reduce the effects of the hysteresis phenomenon and make the controlled model close to linear.The effectiveness of this control design has been demonstrated by Ge and Jouaneh [3].

Controller Design
To represent the nonlinear plant of the piezoelectric actuator, a Preisach model is connected in series with the transfer function () to describe the dynamics system with hysteresis characteristics.The combination of these two components is taken as the controlled plant.

PID Control with Feedforward Compensation.
A PID controller in the continuous time domain can be described as follows: A popular formula for the digital PID control realization is where   ,   , and   are the proportional coefficient, integral coefficient, and differential coefficient, respectively.Δ () is the corresponding increment value.In addition,   =   (/  ),   =   (  /), and  is sampling period.The block diagram of PID control with feedforward compensation is given in Figure 6.
Considering the system stability, response speed, overshoot, and steady-state precision, the tuning roles of   ,   , and   are given as follows.
(a) If   is too small, it will reduce the accuracy.The response speed is slow too.And it will extend the settling time and degrade the system performance.
(b) The role of   is to eliminate the steady-state error of the system.The static error in the system will be reduced faster when   is increased.But if   is too high, it will produce larger overshoot amount.If   is too low, it is difficult to eliminate steady-state error; this will reduce the precision of the system.
(c) The effect of   is to improve the system's dynamic characteristics.It could suppress the change of the error.But if   is too high, it will extend the settling time and reduce the robustness of the system.

Fuzzy Control with Feedforward Compensation.
Fuzzy control is a computer control method on the basis of fuzzy set theory, fuzzy language variables, and fuzzy logic reasoning.Fuzzy controller is the core of fuzzy control, and the key issue of fuzzy controller design is the determination of fuzzy control rules.Fuzzy control rule table is a series of control rules summed up by the expert or the operator according to their manual control experience.The error  and error change rate  are relatively easy to obtain in the control process.Hence, they are employed as the input language variables of the fuzzy controller.In addition,  is output linguistic variable.Thus,  = (, ).As shown in Figure 7, the designed fuzzy controller consists of the steps of fuzzification of inputs, making fuzzy control rules, and defuzzification.After adding the feedforward compensation, the fuzzy tracking control with feedforward compensation is depicted in Figure 8.

Fuzzy PID Control with Feedforward Compensation.
The later simulation results show that fuzzy control with feedforward compensation is not sufficient to produce a satisfactory result.In order to improve the control performance while ensuring the dynamic performance of system, the PID control and fuzzy control are combined together to reduce the shortcomings of each controller.The block diagram of the control scheme is shown in Figure 9.
In this system, the error  and error change rate  are input signals, and the correction values of PID (Δ  , Δ  , and Δ  ) are the outputs.Based on the change of  and , Δ  , Δ  , and Δ  are modified at every time instant to enable the system good dynamic and static characteristics.At last, the values of PID control parameters are obtained.
In the following section, the three types of controllers are implemented and a comparison investigation is carried out through simulation studies.It is notable that, from herein until the end of this paper, () in each following figure is used to express the series connection of the second-order model and the Preisach hysteresis model.

Results of PID Control with Feedforward Compensation.
A PID control with feedforward compensation scheme is realized in MATLAB Simulink, as shown in Figure 10.The PID control parameters are adjusted according to the tuning rules as described in Section 3.2.1.Because of the manual adjustment, the tuning efficiency is low.
Figures 11 and 12 show the system responses to a step input and a sinusoidal input, respectively, where the dashed lines represent the reference inputs and the solid lines represent  the output responses.It can be seen that the system dynamic performance is not good enough, but the steady-state error is small as shown in both cases.

Results of Fuzzy Control with Feedforward Compensation.
Based on the control block diagram, the MATLAB simulation model of the system is developed as shown in Figure 13.
To implement the fuzzy control, the practical values of input variables  and  need to be converted into language variable values.This kind of translation is termed fuzzification, which relies on the membership functions as shown in Figure 14.Then, the language variable values are taken as input.By defining certain control rules, the output fuzzy sets are obtained.This process is called fuzzy inference, as shown in Figure 15.The control rule is derived from expert's experience of operation and control of the system, and they can be edited in the form of fuzzy conditional statement, as shown in Figure 16.At last, the fuzzy output is treated through defuzzification process, which makes the control decision of the system and completes the process of fuzzy control.
Figures 17 and 18 show the system responses to step input and sinusoidal input, respectively.It is seen that the step response of the system is very nice with a fast response, no overshoot, and almost no steady-state error.However, although the dynamic performance is improved, the error of the sinusoidal tracking is large and is not reduced much as compared with PID control.

Results of Fuzzy PID Control with Feedforward Compensation.
In this subsection, a fuzzy PID control is realized to further reduce the steady-state error of the system.The MATLAB simulation model of the fuzzy PID control with feedforward compensation is shown in Figure 19, where the embedded Simulink modules of PID controller and fuzzy controller are shown in Figures 20 and 21, respectively.Generally, different values of  and  require different PID parameter settings.The design objective of the fuzzy PID control is to greatly improve the steady-state control precision without losing too much dynamic performance.For these reasons, the fuzzy control rules are designed.For instance, the setting rules and language description of the parameter   are shown in Table 1.Similar rules are designed for parameters   and   .Figure 22 illustrates the MATLAB settings of the fuzzy control rules.
Moreover, Figures 23 and 24 illustrate the system response to step input and sinusoidal input, respectively.It is found that the system dynamic performance is very good.Most importantly, the sinusoidal signal tracking error is significantly reduced close to zero.In order to further test the fuzzy PID control system with feedforward compensation, more simulation studies have been conducted to examine its performance under different frequencies (5x and 20x) of the input signal.The results are shown in Figure 25.In general, with the improvement of the input frequency, the response of the system will be degraded.But it can be seen from the diagram that the response of the fuzzy PID control system with feedforward compensation does not change much; it keeps a good control result.

Discussion on Control Results
. Preliminary testing shows that the feedforward compensator based on inverse Preisach model is able to mitigate the influence of hysteresis greatly.Thus, the feedforward compensator is employed in the three types of feedback control systems in simulation testing.For a clear comparison, the simulation results of the three kinds of control systems are shown in Table 2.
It is found that the PID control with feedforward compensation delivers a small sinusoidal tracking error, but its dynamic performance is the worst as reflected by the step response results.In addition, the major problem is that the adjustment of PID parameters is a complicated process with low efficiency.
Besides, the fuzzy control with feedforward compensation has great dynamic response.However, its sinusoidal tracking ability is poor, and the control result is not accurate enough for the majority of applications.
Alternatively, the fuzzy PID control with feedforward compensation not only can realize the accurate control similar to PID controller, but also can improve the dynamic performance of system greatly.This is enabled by the designed fuzzy control rules, which are used to modify the PID parameters online, making the system have good learning ability and adaptability.The only problem with fuzzy PID control is that it demands a heavier computation than traditional way.
Additionally, in the aforementioned simulations, the fuzzy control rules are finely tuned to produce the overshoot as small as possible.This implies that the challenge of fuzzy control design lies in the tuning of these inference rules.To meet higher control requirements, more experiences on operation are needed to design more appropriate fuzzy rules.

Conclusions
This paper presents the design and simulation study of fuzzy PID control with feedforward compensation for precision motion control of a piezoelectric actuator.An inverse Preisach model is developed to construct a feedforward compensator.Based on the feedforward compensation, three kinds of feedback controller are designed and realized.Comparative investigations reveal that the fuzzy PID control is superior over PID control and fuzzy control in terms of both step response and sinusoidal response performance.Future work will be conducted to tune the fuzzy rules automatically to reduce the work load of fuzzy control design.

Figure 3 :
Figure 3: Discretization of the model of the Preisach model.

Figure 11 :
Figure 11: Step response of PID control with feedforward compensation.

Figure 12 :
Figure 12: Sinusoidal tracking results of PID control with feedforward compensation.

Figure 13 :Figure 14 :
Figure 13: MATLAB simulation model of fuzzy control with feedforward compensation.

Figure 15 :
Figure 15: The fuzzy inference system editor.

Figure 24 :
Figure 24: Sinusoidal tracking results of fuzzy PID control with feedforward compensation.

Figure 25 :
Figure 25: Responses under different frequencies of input sinusoidal signal.
By cascading the inverse hysteresis compensator and piezoelectric actuator, a linear model is obtained approximately.Furthermore, feedback control can be employed to improve the control precision and enhance the robustness of the system.There are many popular feedback control methods in the literature.This paper employs PID control, fuzzy control, and fuzzy PID control.Moreover, the feedforward control based on the inverse Preisach model and feedback control are combined together to improve the control performance.Without loss of generality, the transfer function () of the plant is represented by a second-order mode.Its expression is shown below: 3.2.Closed-Loop Feedback Control.

Table 1 :
Fuzzy control rule table of parameter   .

Table 2 :
Comparison of the simulation results of the three controllers.