Automation Systems
Lecture 7 - Quality of Control System and PID Controller Tuning
Jakub Mozaryn
Institute of Automatic Control and Robotics, Department of Mechatronics, WUT
Warszawa, 2019
Quality of the control system
Apart from the most important requirement of asymptotic stability, there are imposed additional requirements on control systems, concer- ning the transient (dynamic) response and steady states. They are generally referred to as quality requirements of the control system.
The requirements related to the steady state are formulated by determining the so-called static accuracy of the control system - permissible values of deviations of the system output from the set point in steady states (steady state errors) in the case of disturbances or setpoint changes.
Requirements related to the transient response in the control systems are determined by a number of indices, generally called dynamic quality in- dices of the control system.
Quality of the control system
The task of the control system is to minimize the deviation from the setpoint (when the time approaches infinity) described as an error in steady state:
e(t) = ez(t) + ew(t), (1)
where
ez(t) - error caused by disturbance,
ew(t) - error caused by change of the set point.
Quality of the control system
When rating the quality of the control LTI system, because of the superposition property, both components of the steady state error e(t) = ez(t) + ew(t), can be analyzed separately.
Steady state error caused by disturbance
Transfer function
Gz(s) = ∆ym(s)
z(s) =ez(s)
z(s) = ±Gz(s)Gob(s)
1 + Gob(s)Gr(s) (2) ez(s) = ∆ym(s) = ±Gz(s)Gob(s)
1 + Gob(s)Gr(s)· z(s) (3) Steady state error caused by disturbance
ezst.= lim
t→∞ez(t) = lim
s→0s · ez(s) (4)
ezst. = lim
s→0s · ±Gz(s)Gob(s)
1 + Gob(s)Gr(s)· z(s) (5)
Steady state error caused by change of the set point
Steady state error caused by change of the set point
Steady state error caused by change of the set point
Transfer function
Gew(s) = ew(s)
∆w (s) = −1
1 + Gob(s)Gr(s) (6) ew(s) = −1
1 + Gob(s)Gr(s)∆w (s) (7) Steady state error caused by change of the set point
ewst.= lim
t→∞ew(t) = lim
s→0s · ew(s) (8)
ewst.= lim
s→0s · −1
1 + Gob(s)Gr(s)∆w (s) (9)
Steady state error - example
Determine the steady state error of the control system shown in the figure, caused by a step change of disturbances z(t) = 2 and a step change of the setpoint ∆w (t) = 5. Assume, that in the control system there is used:
P controller, PD controller, PI controller.
Steady state error - example
Transfer function
Gob(s) = kob
(Ts + 1)4 (10)
P controller
Gr(s) = kp (11)
PD controller
Gr(s) = kp(1 + Tds) (12) PI controller
Gr(s) = kp
1 + 1
Tis
(13) Disturbance
z(t) = 2 → z(s) = 2
s (14)
Change of the set point
∆w (t) = 5 → ∆w (s) =5
s (15)
Steady state error caused by the disturbance - example, P controller
ezst.= lim
t→∞ez(t) = lim
s→0s Gob(s)
1 + Gob(s)Gr(s)z(s) (16) P controller
ezst.P = lims→0s Gob(s) 1 + Gob(s)Gr(s)
2 s =
lims→0
kob
(Ts + 1)4· 2 1 + kob
(Ts + 1)4kp
= lims→0
kob· 2 (Ts + 1)4+ kob· kp
(17)
Steady state error caused by the disturbance - P controller ezst.P = kob· 2
1 + kobkp
(18)
Steady state error caused by the disturbance - example, PD controller
ezst.= lim
t→∞ez(t) = lim
s→0s Gob(s)
1 + Gob(s)Gr(s)z(s) (19) PD controller
ezst.PD = lims→0
kob
(Ts + 1)4· 2 1 + kob
(Ts + 1)4kp(1 + Tds)
=
= lims→0
kob· 2
(Ts + 1)4+ kob· kp(1 + Tds)
(20)
Steady state error caused by the disturbance - PD controller ezst.PD= kob· 2
1 + kobkp
(21)
Steady state error caused by the disturbance - example, PI controller
ezst.= lim
t→∞ez(t) = lim
s→0s Gob(s)
1 + Gob(s)Gr(s)z(s) (22) PI controller
ezst.PI = lims→0
kob (Ts + 1)4· 2 1 + kob
(Ts + 1)4kp(1 + 1 Tis)
=
= lims→0
kob· 2
(Ts + 1)4+ kob· kp(1 + 1 Tis)
= 0
(23)
Steady state error caused by the disturbance - PI controller
ezst.PI = 0 (24)
Steady state error caused by the disturbance - summary
P controller
ezst.P = kob· 2 1 + kobkp
(25)
PD controller
ezst.PD= kob· 2 1 + kobkp
(26)
PI controller
ezst.PI = 0 (27)
Steady state error caused by change of the set point - example, P controller
ewst.= lim
t→∞ez(t) = lim
s→0s −1
1 + Gob(s)Gr(s)∆w (s) (28) P controller
ewst.P = lim
s→0s −1
1 + Gob(s)kp
5 s = lim
s→0
−5 1 + kob
(Ts + 1)4kp
= −5
1 + kobkp
(29)
Steady state error caused by change of the set point - P controller ewst.P = −5
1 + kobkp
(30)
Steady state error caused by change of the set point - example, PD controller
ewst.= lim
t→∞ez(t) = lim
s→0s −1
1 + Gob(s)Gr(s)∆w (s) (31) PD controller
ewst.PD = lims→0s −1
1 + Gob(s)kp(1 + Tds) 5 s
= lims→0
−5 1 + kob
(Ts + 1)4kp(1 + Tds)
= −5
1 + kobkp
(32)
Steady state error caused by change of the set point - PD controller ewst.PD = −5
1 + kobkp
(33)
Steady state error caused by change of the set point - example, PI controller
ewst.= lim
t→∞ez(t) = lim
s→0s −1
1 + Gob(s)Gr(s)∆w (s) (34) PI controller
ewst.PI = lims→0s −1 1 + Gob(s)kp
1 + 1
Tis
5 s
= lims→0
−5 1 + kob
(Ts + 1)4kp
1 + 1
Tis
= 0
(35)
Steady state error caused by change of the set point - PI controller
ewst.PI = 0 (36)
Steady state error caused by change of the set point - summary
P controller
ewst.P = −5 1 + kobkp
(37)
PD controller
ewst.PD = −5 1 + kobkp
(38)
PI controller
ewst.PI = 0 (39)
Conclusions about steady state errors
In a control system with a static object and P or PD control algorithm there are non-zero steady state errors in relation to the disturbances or setpoint changes respectively.
Increasing the proportional gain of the P or PD controller reduces the value of static error. Reducing the static deviation by increasing the gain kp is usually limited due to the stability of the system. (The system with PD controller reaches the border of stability at a higher gain than in the case of the regulator P).
Integral action in the controller (PI, PID) provides zero steady state errors in relation to the disturbances or setpoint changes respectively.
Dynamical quality of control system
Requirements related to the transient response in the control systems are determined by a number of indices, generally called dynamic performance quality indicies of the control system.
Groups of such indices are:
transient response indices,
indices descibing the frequency plots of the control system - magnitude and phase margins,
integral indices.
Transient response indices
To evaluate the transient response following indices are used:
Maximum error (dynamical): em - the maximum value of error after the step change of disturbance or setpoint.
Settling time: tr - it is the time between the moment of change of the set point w (t), or introduction of disturbances z(t) , and the moment when the error e(t) reaches a fixed value inside a boundary
∆e(t) (eg.∆e(t) = |0.05emax|).
Overshoot:
κ =
e2
e1
· 100% (40)
where e1 and e2are the first 2 consecutive biggest errors with opposite signs, assuming steady state value of output y (t) after transient response as the zero level (baseline).
Oscillatory transient response - disturbances
Rysunek:Oscillatory transient response of the control system to disturbances:
a) with non-zero steady state error, b) with zero steady state error
Aperiodic transient response - disturbances
Rysunek:Aperiodic transient response of the control system to disturbances: a) with non-zero steady state error, b) with zero steady state error
Oscillatory transient response - setpoint
Rysunek:Oscillatory transient response of the control system to setpoint change: a) with non-zero steady state error, b) with zero steady state error
Aperiodic transient response - setpoint
Rysunek:Aperiodic transient response of the control system to setpoint change: a) with non-zero steady state error, b) with zero steady state error
Selection of controllers
The basic premise when choosing the type of controller is dynamic cha- racteristics of the controlled process.
Rysunek:Control system
Basic equations, describing the properties of the controlled processes Gob(s) = ∆ym(s)
∆u(s) = kob
Tzs + 1e−T0s, Gob(s) = ∆ym(s)
∆u(s) = 1 Tzse−T0s
Selection of controllers
for T0
Tz < 0, 1 ÷ 0, 2 → - switch controllers (two- three- gain controllers),
for 0, 1 ¬ T0
Tz < 0, 7 ÷ 1 ÷ 0, 2 → continuous controllers, dla T0
Tz
> 1 → impulse controllers (impulse output signals).
In the case of industrial processes common ratio of T0
Tz is in the range of 0, 2 ÷ 0, 7. Therefore, in industrial control systems the most common controllers are continuous, with typical control algorithms P, PI, PD and PID.
Selection of controllers
An analysis of the controller algorithm with the process model leads to the following conclusions concerning the selection of the control algorithm:
PI algorithm provides good control only for the low frequencies of setpoint changes or disturbancs. Integral action is necessary to obtain zero error in steady state.
PI algorithm provides wider bandwidth than PID algorithm, but poorer performance for the low frequencies of setpoint changes or disturbances.
Derivative action is recommended for objects with higher order lag (such as thermal processes), because it allows the strong interaction of control even at small deviations. PD controller does not ensure the achievement of zero deviation in steady state . PID algorithm merges to a certain extent the advantages of PI and PD algorithms.
Selection of controllers
In practice, industrial controllers with the continuous algorithm are com- monly used. Their parameters (settings) can be changed (adjusted) within a wide range, so they can control properly processes with different dyna- mical properties.
Depending on the requirements of the stability and quality, the con- troller settings are selected using different selection procedures.
There are following settings of PID controller:
proportional gain kp= 0, 1 ÷ 100 integral gain Ti = 0, 1 ÷ 3600s derrivative gain Td = 0 ÷ 3600s
Selection of controllers
Methods of PID controllers tuning:
Experimental methods - usually do not allow to achieve certain quality of the control system, eg. Ziegler – Nichols, Pessen, Hassen and Offereissen, Cohen-Coon, ¨Astr¨om – Hagglund .
Tabular methods - determining the set of controller parameters based on the parameters of a mathematical model of the controlled process and the required quality criterion of the control system (like the smallest overshoot, shortest settling time. Problem: often the minimization of different quality indices base on the contrary requirements).
Autotuning, eg. relay method.
Tuning of the controllers
Ziegler-Nichols method
Type 1:
controller settings are selected on the basis of the parameters of the closed-loop control system, brought to the border of stability (by experimental excitation of the system).
It can be used to controller tuning in the control systems where processes are described by static and astatic, higher order lag elements.
Type 2:
It can be used to controller tuning in the control systems where processes are described by static higher order lag elements, controller settings are selected based on the transient response of the controlled process.
Ziegler-Nichols method
Rysunek:Functional scheme of real control system
Ziegler-Nichols method, steps 1-3 / 6
Step 1: In the manual mode (M) by changing control variable (CV), adjust the process variable ym (PV) to a state in which it is equal with the required setpoint
Step 2: Set the controller to the proportional action (switch off integral and derivative actions), set the operation point control value of the controller equal to the setting obtained in Step 1 and set the initial value of the controller gain kp> 0.
Step 3: Switch the system to automatic control (A) and if the system maintains equilibrium, by changing SP produce an impulse with some amplitude and pulse duration depending on the expected dynamics of the process; observe or record the change in the controlled variable. It is recommended to use a pulse with the amplitude of 10 % of the process value changes ym (PV) and the pulse duration of about 10 % of the estimated value of the time constant of the controlled process.
Ziegler-Nichols method, steps 4-5 / 6
Step 4: If the transient response is underdamped, set higher values of the proportional gain (Steps 1-3) until the system be on the border of stability (constant oscillations).
Step 5: From the steady oscillations read ’critical’ proportional gain kpkryt.and oscillation period Tosc.
Step 6: Set the patameters according to the table of setings developed by Ziegler-Nichols.
Rysunek:Changes of the process variable (PV) obtained during Ziegler – Nichols experiment
Ziegler-Nichols method
PID controller setting according to Ziegler-Nichols
Controller type kp Ti Td
P 0, 50kpkryt. - -
PI 0, 45kpkryt. 0, 8Tosc • PID 0, 60kpkryt. 0, 5Tosc 0, 12Tosc
Automation Systems
Lecture 7 - Quality of Control System and PID Controller Tuning
Jakub Mozaryn
Institute of Automatic Control and Robotics, Department of Mechatronics, WUT
Warszawa, 2019