Research Article
Design and performance evaluation of an exponentially weighted moving average–based adaptive control for piezo-driven motion platform
Advances in Mechanical Engineering 2018, Vol. 10(6) 1–14 Ó The Author(s) 2018 DOI: 10.1177/1687814018767194 journals.sagepub.com/home/ade
Yung Ting, Tho Van Nguyen and Jia-Ci Chen
Abstract In this article, building a controlled system with velocity feedback in the inner loop for a platform driven by piezoelectric motors is investigated. Such a motion control system is subject to disturbance such as friction, preload, and temperature rise in operation. Especially, temperature rise is an essential problem of using piezoelectric motor, but very few research works address this topic in depth. Exponentially weighted moving average method has been widely used in process control to deal with systematic change and drift disturbance. It is attempted to map the exponentially weighted moving average method and the predictor corrector control with two exponentially weighted moving average formulas into a run-to-run model reference adaptive system for velocity control. Using a predictive friction model, a dead-zone compensator is built that can reduce the friction effect and provide an approximately linear relation of the input voltage and the output velocity for the subsequent exponentially weighted moving average or predictor corrector control control design. Comparison of the exponentially weighted moving average, predictor corrector control, and proportional–integral–derivative controllers is carried out in experiment with different speed patterns on a single-axis and a bi-axial platform. The results indicate that the proposed run-to-run-model reference adaptive system predictor corrector control is superior to the other methods. Keywords Piezoelectric motor, exponentially weighted moving average, predictor corrector control, run-to-run model reference adaptive system
Date received: 2 November 2016; accepted: 20 February 2018 Handling Editor: Mario L Ferrari
Introduction Because of the inherent characteristics of good resolution, small size, and less electro-mechanical impedance (EMI), piezoelectric ceramic motors are popular and are used for high-precision positioning application purposes. Various kinds of piezoelectric motors have been developed in decades and classified into standing wave and traveling wave. Besides the problems of inaccurate assembly or improper preload, these piezoelectric motors are sensitive to temperature rise and their performance becomes worse during long-time operation.1,2
These negative factors are the disadvantages and serious concerns of using piezoelectric motors. Miscellaneous wedge-type piezoelectric motors usually driven by standing wave to propel a carriage via the
Department of Mechanical Engineering, Chung Yuan Christian University, Taoyuan, Taiwan, R.O.C. Corresponding author: Yung Ting, Department of Mechanical Engineering, Chung Yuan Christian University, 32023 Taoyuan, Taiwan, R.O.C. Email:
[email protected]
Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).
2 contact tip were designed to improve the positioning capability and accuracy.3–5 In this article, a wedge-type piezoelectric motor is selected as an example to investigate the disturbance problem of temperature effect in particular. Friction usually causes a dead-zone area which is another problem. For instance, dead zone is apparent when motion starts as well as manipulates forward and reverse repetitively.6 Also, small applied voltage is usually assigned for slow speed requirement neighboring to the set-point change, which may render dead zone in particular. Minh et al. compared with three methods of reducing friction effect while using piezoelectric actuator devices to drive a linear platform. From experimental results, they concluded that the friction compensator based on the LuGre model outperforms the bang-bang control based on a static friction model and the friction state observer based on the LuGre model.7 Lee and Tomizuka designed a control structure coupled with friction compensation for robust high-speed high-accuracy motion control systems. The friction compensator could cope with the friction using a cancellation term in the friction model.8 Following this method, a dead-zone compensator is also designed in this study. On the basis of a predictive friction model included in the control scheme, a linear input/output relation is available, which is important to the subsequent exponentially weighted moving average (EWMA) and predictor corrector control (PCC) design. Miscellaneous approaches utilized self-tuning discrete-time control, sliding mode control, genetic algorithm, neural network methods, and so on to deal with the time-varying parameters from nonlinearity and disturbance so as to preserve high-precision positioning while operating with piezoelectric motors.9–11 Therefore, building a suitable robust controller for such a timevarying system is demanding. EWMA algorithm is one of the most basic run-to-run (RtR) control schemes widely used for monitoring and tracking small changes in the volatility. It has been applied to miscellaneous applications like economics, inventory, and quality control. Process drift and environment disturbance are common factors that often occur in many manufacturing processes. The EWMA method has been extensively studied to overcome such problem in recent decades, especially in semiconductor manufacturing processes. Stability conditions for the process encountered drifting or noise for a first-order MISO (multiple-input-singleoutput) process were derived.12,13 A self-tuning EWMA controller was established to deal with linear drift and random noise, which could improve the semiconductor operating processes.14 A PCC method on the basis of internal model control (IMC) structure was proposed for etching process to remove the impact of process drift and reduce the offset.15
Advances in Mechanical Engineering As known, the relation between the input voltage applied to the piezoelectric motor and the output velocity of the platform is similar to chemical mechanical polishing (CMP) process with linear relation between the input recipe and the output target. Therefore, an attempt is to use EWMA for velocity control objective. Nonetheless, disturbance due to improper assembly, friction, and temperature rise in operation is similar to the CMP process encountered with systematic change, which would cause the predictive linear relation unlikely to keep and appear bias. Increasing the control effort or even deteriorating the system performance is inevitable. Further investigation of solving this problem is thus necessary. Osburn et al.,16 probably the first ones, proposed the model reference system using the performance index minimization method. This adaptive control scheme later has been proved to be one of the most popular methods in the growing field of adaptive control.17 A robust adaptive control scheme based upon the model reference for continuous-time system was designed with an unknown input dead zone to achieve both stabilization and tracking within a desired precision.18 Through adapting system parameters, the problem of model inaccuracy and uncertainty can be overcome so as to gain better robustness. In our previous work, mapping the EWMA into a control system under the IMC configuration was developed.19 Also, the analytic and experimental results proved that the EWMA with a cascade proportional–integral (PI) controller could efficiently deal with the temperature rise problem of stack-type piezoelectric actuator. In this study, the equivalent transfer function of statistical EWMA and PCC is derived. The PCC with two EWMA formulas is proved to be a combination of an integrator and a double integrator. A new approach of implementing these control methods in a model reference adaptive system (MRAS) to form a velocity feedback control scheme is designed. Various velocitytracking trajectories of square and trapezoidal types for single-axis platform as well as circular type for bi-axial platform are assigned for these controllers. The experimental results show that the PCC performs better.
Problem statements Using wedge-type piezoelectric motors to drive a single-axis platform as well as a bi-axial platform is investigated and shown in Figure 1. A nearly linear relation with a constant slope between the applied input voltage and the output velocity of the platform is measured and presented in Figure 2. However, there exists a typical dead zone due to unavoidable friction effect. Besides, a focus in this article is as operating piezoelectric ceramic actuator for a while, so heat is generated to increase temperature around the ceramics. This intrinsic
Ting et al.
3
Figure 3. Temperature rise in operation versus time.
Dead-zone compensator
Figure 1. Wedge-type piezoelectric motor structure: (a) single platform and (b) bi-axial platform.
Figure 2. Output (velocity) versus input (command voltage).
phenomenon deteriorates its performance significantly. In Figure 3, temperature of the piezoelectric motor gradually increases as operation time increases. Using cooling system or intermittent stop during operation may cool down the actuator; however, design of a suitable control method by fully understanding the system behavior subject to temperature rise in operation is the primary goal.
In Figure 2, dead zone occurs in the range of small operating voltage. An equivalent model of the input/ output relation is illustrated in Figure 4.7 An approach to obtain a linear relation depicted in Figure 4(a) is by subtracting Figure 4(b) from Figure 4(c). In this method, the dead zone can be determined based on the voltage-to-velocity profile of the platform in practice. Offset occurs in the region of dead zone and its value may not necessarily be the same in both directions. Depending on the specific driver, motor, and mechanical structure, dead zone can be slightly different.20 The use of a dead-zone compensator requires the command threshold to be established, which is the point when motion starts.6 Using the dead-zone compensator, the relation between input voltage and output velocity corresponding to various operating temperatures (20°, 30°, 40°, 50°, 60°) is measured respectively. Figure 5 shows the selected example of operating temperature at 20°C, 40°C, and 60°C. Approximate slopes of the linear relation for various operating temperatures are found and shown in Figure 6. In this study, normal operating condition of temperature about 20°C is chosen as the reference. Based upon the slopes presented in Figure 6, offset is defined as the difference from the zero point and drift is defined as the difference of linear relation, which are both calculated and listed in Table 1. As seen, drift for each operating condition is found with a linearly increasing trend. Also, it is seen that the linear relation is symmetric for positive and negative applied voltage, which represents the forward and reverse motion of the platform.
Controller design To deal with a time-varying system as addressed above, a modified adaptive control scheme is developed, which applies the RtR control method using EWMA onto an MRAS structure. Based on such a control scheme,
4
Advances in Mechanical Engineering
Table 1. Slope versus operating temperature. Temperature (° C)
Offset a
Slope b
Drift s
20 30 40 50 60
0 2:35 3:248 5:154 7:438
5:945 6:335 6:593 6:736 6:976
0:22 2:25 3:705 5:479 7:884
additional controller added in the feedback path can enhance the robustness to deal with the deterministic drift. Besides the general EWMA, PCC is also considered for its capability of dealing with the drift problem. With the MRAS velocity control loop, an outer loop is simply constructed for positioning control purpose so as to complete a dual control loop.
Modified MRAS
Figure 4. Equivalent output/input model.
The structure of MRAS for velocity feedback control is illustrated in Figure 7, where rn is the assigned velocity command after the dead-zone compensator and also the target value T on the aspect of EWMA; GC is the transfer function of the equivalent EWMA controller ^ P are the or the equivalent PCC controller; GP and G transfer function of the actual and estimated model, respectively; yn and ^yn are the actual and estimated output velocity of the platform, respectively; d is the disturbance. As illustrated in Figure 7, it gives yn = (rn GC en )
GP +d ^P G
ð1Þ
where e is the error between the actual and estimated output velocity, that is, en = yn ^yn . To rewrite (equation (1)), it yields ^ P )un + d (1 + GC )en = (GP G
ð2Þ
Equation (2) is rewritten as 1+
Figure 5. Velocity versus applied voltage (20°C, 40°C, and 60°C).
^ P) GC GP n (GP G rn + d e = ^P ^P G G
ð3Þ
It is assumed that the input rn and the disturbance d are not correlated. By taking discrete-time transformation, the transfer function of output en to independent input rn and d in the MRAS structure is derived respectively as ^P en (z) GP G = n ^ P + GC GP r (z) d = 0 G ^P en (z) G = ^ P + GC GP d(z) rn = 0 G
ð4Þ
ð5Þ
EWMA controller While the dead zone is eliminated, an approximate linear relation of the input (applied voltage) and output (velocity) can be described by21 Figure 6. Slope change versus operating temperature.
Yn = dn + rn Xn + dn + en
ð6Þ
Ting et al.
5
Figure 7. EWMA or PCC controller embedded on an MRAS structure.
where n(= 1, 2, . . . ) is the sample order; Y is the output; X is the input variable (process recipe); d is the system bias or offset; r is the system gain or slope; e is the process white noise (;N(0, s2 )); d is the deterministic disturbance. The predictive model is defined by Guo et al.21 Y^n + 1 = ^dn + r^n + 1 Xn + 1
ð7Þ
where Y^n + 1 is the predictive output whose value is expected to nearly reach the target T after n runs; parameters ^ d and r^ represent the estimated values of the actual d and r in equation (6), respectively. Following equation (7), the input variable is defined as13,21 Xn + 1 =
T ^dn r^n + 1
ð8Þ
and the estimated parameter can be adapted by ^ dn = v(Yn r^n Xn ) + (1 v)^dn1
ð9Þ
where v is the weighting factor in the range of ½0, 1. By writing equation (9) for two successive runs and substituting equation (8) into it, the estimated parameter is rewritten as ^ dn = v(en1 + Y^ n r^n Xn ) + (1 v)^dn1 = v(en1 + ^dn1 ) + (1 v)^dn1 = ven1 + ^dn1
the EWMA contains an integrator of 1=(1 z1 ) for the error to adjust the parameter d. Corresponding to the MRAS block diagram shown in Figure 7, the variables are defined as X [un , Y [yn , and Y^ [^yn . Equations (4) and (5) are given as ^ P )(z 1) en (z) (GP G = n ^ r (z) d = 0 GP (z 1) + vGP ^ P (z 1) en (z) G = ^ d(z) rn = 0 GP (z 1) + vGP
z1 e(z) 1 z1
Dynamic model. It is noted that using the static model GP with constant slope r may have error in computing the output. The control command is defined 6 10 V for input voltage 6 100 V in this study. As shown in Figure 8, a number of step input control commands (0:7, 3:0, 7:0) V are applied to the piezoelectric motor.
ð10Þ
ð11Þ
According to Figure 7, the EWMA controller is defined as GC = vz1 =(1 z1 ). Hence, it is seen that
ð13Þ
Further investigation of the system model either static or dynamic is also concerned.
where en1 is defined as the error between the actual output Yn and the expected output Y^n . A discrete-time expression for equation (10) is derived as ^d(z) = v
ð12Þ
Figure 8. Step input response.
6
Advances in Mechanical Engineering By the Jury stability test,22 the following bounds for the stability region of the RtR-EWMA discrete-time controller embedded on an MRAS structure are given by 0\
vr \2 r^
ð18Þ
For weighting factor v is in the range of ½0, 1, it concludes r 0\ \2 r^
This derived stability condition is same as the conditional convergence for the EWMA controller addressed by Ingolfsson and Sachs.12 The stability region for r^ to guarantee convergence of the sequence fYn g is ½vr=2, + ‘. For step input, the tracking error
Figure 9. Performance of using static or dynamic model.
The settling time is about 50 ms. If the control updating time is designed to be equal to the settling time, it is similar to use a static model but without concerning the transient state. However, the sampling period can be designed much shorter to consider the dynamic effect. An example of using EWMA employed with sampling period Ts = 5 ms for a step input of 10 mm/s is investigated and shown in Figure 9. It is seen that large overshoot occurs while using static model. Obviously, using dynamic model to estimate the output velocity ^yn during the transient state for computing an appropriate new control input is necessary. With sampling period Ts = 5 ms, system behavior is identified as a secondorder form given by 2 ^ P (z) = b2 z + b1 z + b0 G 2 z + a1 z + a0
ð14Þ
a0 = 0:3653, where coefficients a1 = 1:352, b2 = 1240, b1 = 2140, and b0 = 966:2. Stability condition for EWMA. According to Figure 6, the linear relation between the input voltage and output speed is obtained with a constant slope r by means of dead-zone compensation. Thus, the actual and the estimated models are assigned with linear relations as ^ P = r^, respectively. Equations (12) and GP = r and G (13) can be rewritten as en (z) (r r^)(z 1) = rn (z)d = 0 r^(z 1) + vr en (z) r^(z 1) = d(z) rn = 0 r^(z 1) + vr
(r r^)(z 1) 1 enr jss = lim (1 z1 ) z!1 r^(z 1) + vr 1 z1
ð20Þ
=0 As seen, the steady-state error of the system converges to zero; in other words, the output Yn converges to the target T. Disturbance and steady-state error. As described above, the suitable operating condition of temperature is assigned as 20°C, which is chosen as the reference index. As seen in Figure 6, a linear drift exists, which is treated as a deterministic ramp disturbance in the sense of control design. In equation (15), the error due to a ramp disturbance becomes end jss = lim½(1 z1 )end (z) z!1 " # r^(z 1) gz1 1 = lim (1 z ) z!1 r^(z 1) + vr (1 z1 )2 ð21Þ =
g^ r vr
where g is a constant. The steady-state error found in equation (21) is g^ r=vr is same as the value computed by statistical method.12 That is, like the condition of linear drift noise (d = ns) happens in the aging problem in CMP process or subject to wear-out, EWMA method described in equations (6)–(8) is proved unable to obtain zero steady-state error.12,23
ð15Þ
PCC controller ð16Þ
As seen, the characteristic equation is D(z) = r^z (^ r vr)
ð19Þ
ð17Þ
As described above, the EWMA controller is not able to deal with the problem of ramp-type disturbance theoretically. Instead of using EWMA with one-stage estimation, a PCC is proposed with two-stage parameter estimation15 and expressed by
Ting et al.
7
Xn + 1 =
^n T ^dn f r^n + 1
ð22Þ
Then, the transfer functions for en to rn and en to d in the MRAS structure are derived as
ð23Þ
en (z) (z 1)(z 1)2 ð31Þ = z2 + (zvS 2)z + z(vP vS ) + 1 rn (z)d = 0 en (z) (z 1)2 ð32Þ = z2 + (zvS 2)z + z(vP vS ) + 1 d(z) rn = 0
where ^ dn = v1 (Yn r^n Xn ) + (1 v1 )^dn1 ^ n = v2 (Yn r^n Xn ^dn1 ) + (1 v2 )f ^ n1 f (0\v1 , v2 \1)
ð24Þ
v1 and v2 are the weights for equations (23) and (24), ^ n is used to compensate for the error respectively, and f incurred by ^ dn . Equation (23) can be rewritten as ^ dn = v1 (en1 + Y^ n r^n Xn ) + (1 v1 )^dn1 ^ n1 + ^dn1 = v1 en1 + v1 f
ð25Þ
where z denotes the ratio of the slopes of the actual plant and the estimated model, r=^ r. The characteristic equation of (31) and (32) can be expressed as D(z) = z2 + (zvS 2)z + z(vP vS ) + 1
Using the Jury stability test, the stability conditions are determined as 8 .0 < zvP 4 + z(vP 2vS ) .0 : jz(vP vS ) + 1j \1
A discrete-time expression for equation (25) is derived as ^ d(z) =
z1 ^ ½e(z) + f(z)v 1 1 z1
ð26Þ
Similarly, equation (24) can be rewritten as ^n f
= v2 (en1 + Y^ n r^n Xn ^dn1 ) ^ n1 + (1 v2 )f ^ = v2 en1 + fn1
ð33Þ
ð34Þ
Substituting vS = v1 + v2 and vP = v1 v2 into equation (34) yields (i) (ii) (iii)
ð27Þ
zv1 v2 .0 4 + zv1 v2 2z(v1 + v2 ).0 0\z(v1 + v2 v1 v2 )\2
ð35Þ
A discrete-time expression for equation (27) is derived as
From (i) in equation (35), it can be concluded that z.0. Solving (ii) and (iii), one can have
1 ^ = v2 z e(z) f(z) 1 z1
r 4 \ r^ 2v1 + 2v2 v1 v2
ð36Þ
r 2 \ r^ v1 + v2 v1 v2
ð37Þ
ð28Þ
From equations (26) and (28), it yields ^ ^ = d(z) + f(z)
1 z (v1 + v2 ) 1z 1 1 2 z e(z) + v1 v2 1z 1
and
h
ð29Þ
From equations (36) and (37), it is clear to conclude
Hence, it is seen that the PCC contains an integrator as well as a double integrator as shown in equation (29). Stability condition for PCC controller. According to Figure 7 and equation (29), the transfer function of the equivalent PCC controller in the MRAS structure can be derived as GC = v S
1 2 z1 z + v P 1 1z 1 z1
r 4 2 \ \ r^ 2v1 + 2v2 v1 v2 v1 + v2 v1 v2
Thus, the stability bound of the PCC controller in the MRAS structure is chosen as r 4 0\ \ r^ 2v1 + 2v2 v1 v2
ð39Þ
An attempt is to narrow the region of choosing weighting factors v1 and v2 , so equation (39) is rewritten as
ð30Þ
where vS = v1 + v2 and vP = v1 v2 ; both weighting factors v1 and v2 are in a range of [0,1].
ð38Þ
v1 (2 v2 ) + 2v2 \4 then
r^ r
ð40Þ
8
Advances in Mechanical Engineering
v1 \
4 rr^ 2v2
ð41Þ
2 v2
Table 2. Velocity-tracking error (mm/s)—square. Controls
Thus, the feasible region of (v1 , v2 ) is determined as ( Q = (v1 , v2 )j0\v1 min 1,
4 rr^ 2v2 2 v2
!
) , 0\v2 \1 ð42Þ
Appropriate weighting factors v1 and v2 can be chosen following the above constraints for stability.
PID EWMA PCC
Conditions At 20°C
At 40°C
At 60°C
0:688114 0:691712 0:678228
0:865387 0:878629 0:671148
0:959172 0:935570 0:685127
PID: proportional–integral–derivative; EWMA: exponentially weighted moving average; PCC: predictor corrector control.
Table 3. Velocity-tracking error overshoot—square.
Disturbance and steady-state error of PCC controller. PCC is proved as an integral-double-integral controller (I-II) in equation (30) and in previous research.15,21 As seen, PCC is quite similar to EWMA but with two-stage parameter estimation. Besides adjustment of the estimate parameter d, the additional parameter f is useful to deal with the error addressed in equation (21). Using PCC, the steady-state error due to a ramp disturbance is given by
Controls
PID EWMA PCC
Conditions At 20°C
At 40°C
At 60°C
0:174 0:123 0:082
0:235 0:159 0:091
0:203 0:183 0:116
PID: proportional–integral–derivative; EWMA: exponentially weighted moving average; PCC: predictor corrector control.
1 n end j= ss lim½(1 z )ed (z) z!1
= lim½(1 z1 )
Table 4. Velocity-tracking error (mm/s)—trapezoidal.
(z 1)2 gz1 2 z + (zvS 2)z + z(vP vS ) + 1 (1 z1 )2 ð43Þ (z 1)g = lim 2 z!1 z + (zvS 2)z + z(vP vS ) + 1 =0
Controls
z!1
PID EWMA PCC
Conditions At 20°C
At 40°C
At 60°C
0:313228 0:249927 0:237904
0:439307 0:357493 0:274411
0:662491 0:520866 0:297553
PID: proportional–integral–derivative; EWMA: exponentially weighted moving average; PCC: predictor corrector control.
Therefore, PCC is able to overcome the problem of process drift and converge to the desired target, which is also proved by statistical control method in terms of the expected value and variance.21,23
Table 5. Velocity-tracking error overshoot—trapezoidal. Controls
Experiment and discussion A single-axis platform and a bi-axial platform driven by wedge-type piezoelectric motors were used for testing. An optical quadrature encoders (MicroE Systems II 5600 with 1.22 nm resolution per count) connected to NI PCI-DAQ 6229 via DB-15 was used to provide feedback of the instantaneous position and velocity for the platform. Sampling period of 5 ms is chosen for using a regular PC-based control system. Based on the feasible region defined in equation (38), the weighting factors for EWMA and PCC are selected to be (v1 = 0:252) and (v1 = 0:252 and v2 = 0:011), respectively. Proportional–integral–derivative (PID) control is also implemented for comparison with the proposed control method. The PID gains are selected as Kp = 0:00018, Ki = 0:00000326, and Kd = 0:00000045
PID EWMA PCC
Conditions At 20°C
At 40°C
At 60°C
0:073 0:065 0:044
0:086 0:078 0:052
0:151 0:118 0:058
PID: proportional–integral–derivative; EWMA: exponentially weighted moving average; PCC: predictor corrector control.
by means of the minimization of the integral-of-timemultiplied by absolute error (ITAE).24 For a single-axis platform, two types of tests were chosen with arbitrarily chosen trajectories (square command standing for constant velocity and trapezoidal command standing for constant velocity and acceleration). The platform is driven continuously by piezoelectric motor from temperature 20°C up to 60°C. Selected
Ting et al.
9
Figure 10. Velocity-tracking and error performance (20°C): (a) tracking performance, (b) tracking error, and (c) zoom of (b).
Figure 11. Velocity-tracking and error performance (40°C): (a) tracking performance, (b) tracking error, and (c) zoom of (b).
examples of operating temperature (20°C, 40°C, 60°C) are presented with implementation of PID, EWMA, and PCC control methods associated with the deadzone compensator. Velocity-tracking performances for
square command and trapezoidal command are shown in Figures 10–15. The root-mean-square error (RMSE) and error overshoot computed for each case are presented in Tables 2–5. Since the trapezoidal command is
10
Advances in Mechanical Engineering
Figure 12. Velocity-tracking and error performance (60°C): (a) tracking performance, (b) tracking error, and (c) zoom of (b).
Figure 13. Velocity-tracking and error performance (20°C): (a) tracking performance, (b) tracking error, and (c) zoom of (b).
a smoother trajectory, it is apparently seen that the RMSE of the trapezoidal command is smaller than that of the square command. From the results, it indicates that PCC outperforms EWMA controller as expected.
Both EWMA and PCC perform better than the PID controller. It is obvious to see the PCC controller has the advantages of minimum overshoot while moving direction (forward and reverse) is changing as shown in
Ting et al.
11
Figure 14. Velocity-tracking and error performance (40°C): (a) tracking performance, (b) tracking error, and (c) zoom of (b).
Figure 15. Velocity-tracking and error performance (60°C): (a) tracking performance, (b) tracking error, and (c) zoom of (b).
(b) and (c) of Figures 10–12. Also, velocity-tracking performance is focused on the critical condition of setpoint change during the switching from acceleration to
constant velocity as shown in (b) and (c) of Figures 13– 15. It is found that the PCC controller has the smallest error on an average. As the temperature increases, the
12
Advances in Mechanical Engineering
Figure 19. Velocity-tracking plot of y-axis (20 mm).
Figure 16. Compilation of contour tracking plot (20 mm).
Figure 20. Velocity-tracking error of y-axis (20 mm).
Table 6. Norm of error (mm/s) on 20-mm circular trajectory.
Figure 17. Velocity-tracking plot of x-axis (20 mm).
Controls
Norm of error
PID EWMA PCC
0:63337 0:44332 0:34752
PID: proportional–integral–derivative; EWMA: exponentially weighted moving average; PCC: predictor corrector control.
Figure 18. Velocity-tracking error of x-axis (20 mm).
tracking error and overshoot are more serious. Hence, a conclusion is drawn that the proposed PCC controller preserves better performance.
A contour test is also carried out on a bi-axial platform to evaluate the performance of the proposed control method with sinusoidal command. As presented in Figure 16, a circular trajectory with a stroke of diameter 20 mm is tested under operating temperature 60°C. A simple position feedback is designed in the outer loop for tracking position. In the zoom area, it is seen the PCC preserves better positioning accuracy. The tracking performance on x-axis and y-axis are shown in Figures 17–20 as well as listed in Table 6. The tracking performance on both axes shows that norm of the error of the proposed PCC control method is (0.3475) mm/s is superior to the EWMA and PID with (0.4433) mm/s and (0.6333) mm/s, respectively.
Ting et al.
Conclusion With a dead-zone compensator, friction is effectively reduced to give an approximate linear input/output relation for the following control design using EWMA or PCC. Temperature rise in operation affects the system performance significantly, which is one of the main focuses in this article. From experiment, the effect due to temperature rise is found as a linear drift bias and treated as a deterministic ramp disturbance. Based upon the configuration of an MRAS, a velocity feedback control employed with EWMA or PCC method is developed. The error due to ramp disturbance is calculated and verified to be the same as that using statistical control. As expected, the PCC is able to cope with the ramp disturbance caused by temperature rise while operating piezoelectric motors. Mapping onto a control block diagram and under a configuration of MRAS also gives a convenient way in design and analysis. Different speed patterns for a single-axis platform and a circular planar motion for a bi-axial platform are carried out to verify the effectiveness of the proposed method. All the testing results indicate that the EWMA and PCC controllers are verified to outperform the PID with optimal gains and PCC is superior to the other methods. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by the National Science Council under the grant NSC101-2622-E-033-001-CC2 and Soundwide Technology Corp.
References 1. Senjyu T, Miyazato H, Yokoda S, et al. Speed control of ultrasonic motors using neural network. IEEE Trans Power Electr 1998; 13: 381–387. 2. Yung Ting, Chun-Chung Li and Cheng-Min Lin, ‘‘Controller Design for High-frequency Cutting Using A Piezo-driven Microstage’’, Precision Engineering 2011; 35: 455–463. 3. Schmoeckel F and Worn H. Remotely controllable mobile microrobots acting as nano positioners and intelligent tweezers in scanning electron microscopes (SEMs). In: 2001 proceedings of the IEEE international conference on robotics and automation, Vol. 4, Seoul, 21–26 May 2001, pp.3909–3913. New York: IEEE. 4. Urban S and Stefan J. Fabrication of monolithic piezoelectric drive units for a miniature robot. J Micromech Microeng 2002; 12: 582–589.
13 5. Kawano H, Ando H, Hirahara T, et al. Application of a multi-DOF ultrasonic servomotor in an auditory teleexistence robot. Trans Rob 2005; 21: 790–800. 6. Xu J-X, Xu J and Lee TH. Iterative learning control for a linear piezoelectric motor with a nonlinear input deadzone. In: Proceedings of the 2004 IEEE international conference on control applications, Vol. 2, Taipei, 2–4 September 2004, pp.1001–1006. New York: IEEE. 7. Minh TN, Ohishi K, Takata M, et al. Accurate friction compensation for precision stage using synchronous piezoelectric device driver. In: 2006 9th IEEE international workshop on advanced motion control, Istanbul, 27– 29 March 2006, pp.109–114. New York: IEEE. 8. Lee HS and Tomizuka M. Robust motion controller design for high-accuracy positioning systems. IEEE Trans Ind Electr 1996; 43: 48–55. 9. Snitka V, Zukauskas D and Mizariene V. Positioning to nanometer resolution with ultrasonic actuators. In: Transducers ’95. The 8th international conference on solidstate sensors and actuators and Eurosensors IX, Stockholm, 25–29 June 1995, pp.409–412. New York: IEEE. 10. Wai R-J and Lee J-D. Intelligent motion control for linear piezoelectric ceramic motor drive. IEEE T Syst Man Cy B 2004; 34: 2100–2111. 11. Wai R-J and Tu C-H. Design of total sliding-mode-based genetic algorithm control for hybrid resonant-driven linear piezoelectric ceramic motor. IEEE Trans Power Electr 2007; 22: 563–575. 12. Ingolfsson A and Sachs E. Stability and sensitivity of an EWMA controller. J Qual Tech 1993; 25: 271–287. 13. Sachs E, Hu A and Ingolfsson A. Run by run process control: combining SPC and feedback control. IEEE Trans Semicond Manufac 1995; 8: 26–43. 14. Smith TH and Boning DS. A self-tuning EWMA controller utilizing artificial neural network function approximation techniques. IEEE T Compon Pack A 1997; 20: 121–132. 15. Butler SW and Stefani JA. Supervisory run-to-run control of polysilicon gate etch using in situ ellipsometry. IEEE Trans Semicond Manufac 1994; 7: 193–201. 16. Osburn PV, Whitaker HP and Kezer A. New developments in the design of model reference adaptive control systems. IAS 29th annual meeting, Institute of Aeronautical Sciences, New York, New York, January 1961, paper No. 61-39. 17. Parks P. Liapunov redesign of model reference adaptive control systems. IEEE Trans Automat Cont 1966; 11: 362–367. 18. Wang XS, Hong H and Su CY. Model reference adaptive control of continuous-time systems with an unknown input dead-zone. IEEE Proc Cont Theory Appl 2003; 150: 261–266. 19. Ting Y, Li CC and Lin CM. Controller design for highfrequency cutting using a piezo-driven microstage. Precision Eng 2011; 35: 455–463. 20. Nanomotion. The piezoelectric effect, 2008, http:// www.nanomotion.com/piezo-ceramic-motor-technology/ piezoelectric-effect/ 21. Guo R-Y, Chen A and Che J-J. Run-to-run control schemes for CMP process subject to deterministic drifts. In:
14 Semiconductor manufacturing technology workshop, Hsinchu, Taiwan, 14–15 June 2000, pp.251–258. New York: IEEE. 22. Ogata K. Discrete-time control systems. Upper Saddle River, NJ: Prentice Hall, 1987.
Advances in Mechanical Engineering 23. Castillo ED. Statistical process adjustment for quality control. New York: Wiley, 2002. 24. Astrom KJ and Hagglund T. PID controllers—theory, design, and tuning. 2nd ed. Paris: ISA, 1995.