Comparative Study of PID Tuning Methods For Processes With Large & Small Delay Times
Comparative Study of PID Tuning Methods For Processes With Large & Small Delay Times
Abstract— Kp
A series of simulations have been carried out to compare the PI and
PID controller tunings of different tuning methods, namely Ziegler-
Nichols(ZN), Cohen-Coon (CC), Chien-Hrones-Reswick(CHR) and
minimum error criteria ISE, ISTE, ISTSE methods for various process
R(s) E(s) Ki
C(s)
models. Using simple first order plus time delay (FOPTD) and second − s
Gp (s)
order plus time delay (SOPTD) processes in a feedback control loop,
simulations of system responses to setpoint changes were plotted and
analyzed for speed of response, stability, and robustness of these
tuning methods. It was found that minimum error criteria ISE, ISTE, Kd s
ISTSE tuning methods, although they give a relatively slow response,
were superior in stability and robustness in almost all cases. The
minimum error tuning method is easy to implement and gives the
Fig. 1. Basic control system configuration of PID
desired results using MATLAB/Simulink effortlessly. It gives faster
responses with less oscillation. This superiority is observed both for
PI and PID controllers. On the other hand, ZN and CHR gave larger
overshoots with longer settling time for PI and PID, while CC gives
very sluggish response for PI controller. The ZN and CHR tuning
methods have higher proportional gains and smaller integral time Ki
Gc = Kp + + Kd s
constants leading to very poor damping, thus they are only suitable s (1)
for processes that operate deep within the stable region while for 1
processes operating on the periphery of the stable region will be = Kp 1 + + τd s
τi s
unsafe to tune using ZN and CHR tuning methods.
Keywords—PID tuning, Ziegler-Nichols, Cohen-Coon, ISE, ISTE, where Kp , Ki and Kd are proportional gain, integral
ISTSE.
gain and derivative gain constants respectively, while τi is
integral time constant and τd is derivative time constant. The
I. I NTRODUCTION proportional controller is used when we want the controller
action to be proportional to the size of the error signal. The
Kp
Gp (s) = e−τd s (3) A. Ziegler-Nichols (ZN) step response method
(τ1 s + 1)(τ2 s + 1)
The ZN tuning method determines the value of Kp , τi
where Kp is the system gain, τd is the time delay, τ1 > 0 and τd centered on transient response of the plant. Most PID
and τ2 > 0 are the time constants. The selection of a tuning rules are based on the assumption that the plant can be
controller type (P, PI, PID) and its parameters is intimately approximated by a first-order plus time delay system, depicted
related to the model of the process to be controlled. Their in Eq. (4) whose unit-step response resemble an S-shaped
computations required that the response should have 14 decay curve with no overshoot, as depicted in Fig. 2.
ratio, minimum offset, minimum area under the load-response
curve, and other favorable properties [1] . Ideally, it is necessary k
G(s) = e−sL (4)
the closed-loop system to satisfy the following performance Ts + 1
criteria [2] : the closed-loop system must be stable; the effects The S -shaped step response is termed as the Process Re-
of disturbances are minimized, providing good disturbance action Curve in process [15,16] . Such step response curves can
rejection; rapid & smooth responses to set-point changes are be generated experimentally or from a dynamic simulation of
obtained, that is, good set-point tracking; steady-state error the plant. The S-shaped reaction curve can be characterized by
(offset) is eliminated; excessive control action is avoided; the two constants, delay time L and time constant T , which are
control system is robust, that is, insensitive to changes in determined by drawing a tangent line at the inflection point of
process conditions and to inaccuracies in the process model. the curve and finding the intersections of the tangent line with
Two FOPTD processes were selected for tuning purposes, one the time axis and the steady state level line, as shown in in Fig.
with large delay time and one with small delay time as well 2. The time constant can be identified by taking the maximum
as an SOPTD process. slope and measuring the time period between the points where
the maximum slope line crosses the initial and final response test as is more sensitive than the ZN method. The controller
lines. So this method is useful even when transfer function parameters with this method are given by formulae in Table
of the plant is not known. First a step response of the plant II, where a = kL
T and τ = L+T .
L
y(t)
C. Chien-Hrones-Reswick (CHR) PID Tuning Algorithm
The CHR PID tuning method, which is also called modified
k
ZN method, focuses on set-point regulation and disturbance
rejection. The method provides formula for 0% and 20%
overshoot. The CHR PID tuning formulas are summarized
in Table III for set point regulation. Compared with ZN
method CHR PID tuning formula uses the time constant, T , of
the process explicitly. The more heavily damped closed-loop
response, which ensures, for the ideal plant model, the quickest
response without overshoot is labeled with 0% overshoot, and
t the quickest response with 20% overshoot is labeled with 20%
a L T overshoot. Similarly, this method can also be used to design
controllers for disturbance rejection purposes, though only set
point regulation is considered in this work for comparison
purpose against other tuning methods.
TABLE III
CHR T UNING FOR SET- POINT REGULATION
Fig. 2. Open-loop process reaction curve (S-shaped response curve) for a
step change Control Mode with 0% overshoot with 20% overshoot
Kp τi τd Kp τi τd
0.3 0.7
TABLE I P a a
ZN T UNING RULE BASED ON S TEP R ESPONSE OF P LANT
0.35 0.6
PI a
1.2T a
T
Control Mode Kp τi τd
0.6 0.95
PID a
T 0.5L a
1.4T 0.47L
1
P a
0.9
PI a
3L
1.2 L
PID a
2L 2
D. Optimum PID Controller Design
The optimum settings (parameters) are usually considered
to be those which give a minimum error integral after a step
change in set point or load. Ideal reset action is assumed so
B. Cohen Coon (CC) Method that the error eventually becomes zero and the error integral
In this method control action is removed and an open loop reaches a finite limit. Several authors have used the integral
transient is introduced by a unit step change in the signal to of the square of the error as a criterion of control quality. In
the process. At the output of the measuring element the step this work, the time domain optimal tuning method of FOPID
response is recorded as process reaction curve as shown in controllers has been formulated for the control of higher
Fig. 2. Then the dynamics of process is approximated by a order processes. This technique searches for an optimal set of
first order plus transportation lag model. The CC method is controller parameters while minimizing a suitable time domain
a more complex version of the ZN method. In this method integral performance index [14] . The simple error minimization
the process reaction curve is obtained first, by an open loop criteria can be customized by a suitable choice of a time
domain performance index to have a better control action as engineering applications [19] . The proportional gain, Kp , of ZN
reported by Guzmán et al. [17] and Saptarshi et al. [18] . This and CHR methods are about twice as large as CC’s Kp . While
introduces an alternative controller design relation based on the Kp of ISE, ISTE, ISTSE are higher than that of CC’s
a performance index that considers the entire closed loop Kp but lower than ZN’s and CHR’s Kp . The time domain
response. In this work the three popular integral error criteria specifications were also used for comparing these tuning
(optimum criteria) were considered, the integral squared error methods as performance index criteria, including rise time,
(ISE) criterion, integral squared time weighted error (ISTE) setting time, overshoot, peak time, peak as depicted in Table V,
criterion, and the integral squared time-squared weighted error where less error will be considered as the best controller. The
(ISTSE) criterion. The ISE criterion penalizes large errors, shortest settling time was for the CC settings for a step change,
while the ISTE and ISTSE criteria penalize errors that persist with a value of 673.90 seconds. Similarly, ZN and CHR
for long periods. In addition the following performance crite- methods have the shortest rise time and peak times compared
rion were considered for quantitative analysis of good control, with CC as well as longest settling time and peak time too. The
the good overshoot and short settling time, settling time, peak difference in rise time after the change in setpoint was quite
value and peak time. small among all the tuning methods except the CC method,
1) Integral of Square Error (ISE): as shown in Table V. The overshoot of ZN and CHR methods
Z ∞ were found way far higher than CC method. ZN method
ISE = e2 (t)dt (5) generally aims at obtaining fast disturbance compensation [20]
0 and results in poor tuning. Despite the popularity ZN is largely
2) Integral of Time multiplied Square Error (ITSE): due to the fact that it was amongst the first tuning methods to
Z ∞ be proposed, and compared to most other tuning techniques it
ITSE = te2 (t)dt (6) is still the simplest to use. As one can expect the optimum error
0 criteria (ISE, ISTE, ISTSE) gave least oscillatory responses. It
3) Integral of Squared Time multiplied by Square Error is also found that the optimum error tuning methods gave less
(ISTSE): Z ∞ overshoot peaks (about 1.2) for a change in setpoint, while the
ISTSE = t2 e2 (t)dt (7) ZN and CHR have as high as about 1.8. The integral time of
0 ISE, ISTE, ISTSE methods on the other hand were higher than
ZN, CHR and Cohen-Coon, except ISE having lower than that
IV. R ESULTS AND D ISCUSSION of CC.
Control tuning parameters and performances were calcu-
lated/simulated for different tuning techniques considered in
Step Response
this work. The controller tuning parameters estimated by ZN, 2
CC, CHR and minimum error criteria (ISE, ISTE, ISTSE)
1.8
methods are presented below. With these tuning parameters
obtained for P, PI and PID controllers by ZN, CC, CHR 1.6
methods and the optimum error criteria (ISE, ISTE, ISTSE) for 1.4
delay time
220
Gp (s) = e−2s (8) Fig. 3. Step responses for FOPDT process with PID controller tuned by
6.5s + 1 different turning methods: ZN +, Cohen-Coon ◦, CHR 4, ISE ∗, ISTE B,
the step repones of PID closed loop systems with different ISTSE •.
tuning techniques are depicted in Fig. 3. With CHR and
ZN PID tuning methods, the resulting system exhibit a large TABLE IV
maximum overshoot in the step response method, which is PID TUNING PARAMETERS Kp , τi AND τd , FOR A FOPTD PROCESS WITH
unacceptable. In such a case, we need a series of fine-tuning LARGE DELAY TIME
methods gave large overshoots with longer settling time, while 1.4
Amplitude
overshoot similar with what was observed for PID controller 1
0.2
Step Response 0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
2 Time (seconds)
1.8
1.6
Fig. 5. Step responses for SOPDT process with PID controller tuned by
1.4 different turning methods: ZN +, Cohen-Coon ◦, CHR 4, ISE ∗, ISTE B,
ISTSE •.
1.2
Amplitude
1
TABLE VII
0.8
PID TUNING PARAMETERS Kp , τi AND τd , FOR A SOPTD PROCESS WITH
LARGE DELAY TIME
0.6
TABLE VIII
T IME DOMAIN PERFORMANCE INDEXES OF THE TUNED PID, RESPONSE
Fig. 4. Step responses for FOPDT process with PI controller tuned by
PARAMETER VALUES OF SOPDT
different turning methods: ZN +, Cohen-Coon ◦, CHR 4, ISE ∗, ISTE B,
ISTSE •.
ZN CC CHR ISE ISTE ISTSE
Rise Time 0.06 0.16 0.067 0.11 0.13 0.14
TABLE VI Settling Time 0.79 0.49 1.17 0.29 0.43 0.47
T IME DOMAIN PERFORMANCE INDEXES OF THE TUNED PI, RESPONSE
Overshoot 67.87 9.83 77.92 6.09 9.51 10.79
PARAMETER VALUES OF FOPDT
Peak 1.68 1.10 1.78 1.06 1.10 1.11
PeakTime 0.22 0.33 0.22 0.22 0.28 0.31
ZN CC CHR ISE ISTE ISTSE
Rise Time 130.32 1022.90 107.68 161.46 222.72 288.44 The tuning methods for PI controller showed the CC and
Settling Time 1636.60 2052.20 3291.30 1707.30 1069.90 820.76
Overshoot 57.21 0.00 90.78 23.12 11.63 5.41
the minimum error tuning methods have good performance
Peak 1.57 1.00 1.91 1.23 1.12 1.05 with relatively small overshoots, and settling times as depicted
Peak Time 437.98 4854.70 425.37 436.00 521.42 614.38 in Fig. 6 and Table IX. The ZN and CHR tuning methods
gave large overshoots with longer settling time. These tuning
A second process considered is a FOPDT with longer delay methods result in closed-loop systems with very poor damping.
It was observed that these, ZN and CHR tuning methods Step Response
were only suitable for open-loop stable processes and the 1.8
Amplitude
two methods also result in poor tuning. 1
0.8
0.6
Step Response
0.4
2.5
0.2
0
0 1 2 3 4 5 6 7 8 9 10
2 Time (seconds)
1.5
Fig. 7. Step responses for SOPDT process, (Gp(s) = 10s2 100 e−0.1s ),
Amplitude
+11s+1
with PID controller tuned by different turning methods: ZN +, Cohen-Coon
◦, CHR 4, ISE ∗, ISTE B, ISTSE •.
1
TABLE X R EFERENCES
T IME DOMAIN PERFORMANCE INDEXES OF THE TUNED PI, RESPONSE
PARAMETER VALUES OF SOPDT [1] S. E. LeBlanc and D. R. Coughanowr, Process Systems
Analysis and Control, 3rd ed. McGraw-Hill, 2009.
ZN CC CHR ISE ISTE ISTSE
Rise Time 0.11 0.27 0.11 0.16 0.19 0.21 [2] D. E. Seborg, D. A. Mellichamp, T. F. Edgar, and F. J.
Settling Time 7.19 1.29 8.34 1.53 1.02 1.08 Doyle III, Process Dynamics and Control, 3rd ed. Wiley,
Overshoot 103.08 12.52 105.25 34.92 26.60 22.87 2011.
Peak 2.03 1.13 2.05 1.35 1.27 1.23
PeakTime 0.34 0.54 0.34 0.39 0.43 0.45 [3] L. Sun, D. Li, K. Hu, K. Y. Lee, and F. Pan, “On tuning
and practical implementation of active disturbance
PID controller gave better performance compared to P and rejection controller: A case study from a regenerative
PI controllers with reference to settling time, rise time, offset, heater in a 1000 mw power plant,” Industrial and
ISE, IAE and IATE. Engineering Chemistry Research, vol. 55, no. 23, pp.
Similar trends were observed for SOPDT process, Gp(s) = 6686–6695, 2016, cited By 0.
100 −0.1s [4] M. Shamsuzzoha, “A unified approach for proportional-
10s2 +11s+1 e , where the ZN and CHR found to be of
higher overshoots and oscillatory responses with longer set- integral-derivative controller design for time delay
tling times. Cohen-Coon method shows better performance processes,” Korean Journal of Chemical Engineering,
close to the minimum error tuning methods ISE, IAE and vol. 32, no. 4, pp. 583–596, 2015, cited By 3.
IATE. [5] S. Zheng, X. Tang, and B. Song, “Graphical tuning
method of fopid controllers for fractional order uncertain
V. C ONCLUSION system achieving robust d -stability,” International
Different PID tuning methods are discussed in this paper. Journal of Robust and Nonlinear Control, vol. 26, no. 5,
Step response of all methods and their system characteristics pp. 1112–1142, 2016, cited By 0.
[6] S. Skogestad, “Simple analytic rules for model reduction
and PID controller tuning,” modeling, identification and
control, vol. 25, pp. 85–120, 2004.
[7] C. Grimholt and S. Skogestad, “Optimal pid-control on
first order plus time delay systems & verification of
the SIMC rules,” IFAC Proceedings Volumes, vol. 46,
no. 32, pp. 265 – 270, 2013, 10th IFAC International
Symposium on Dynamics and Control of Process
Systems.
[8] S. W. Sung and I.-B. Lee, “An improved algorithm
for automatic tuning of PID controllers,” Chemical
Engineering Science, vol. 55, no. 10, pp. 1883 – 1891,
2000.
[9] D. V. S. Strčnik and Ä. Juričić, “A magnitude optimum
multiple integration tuning method for filtered PID
controller,” Automatica, vol. 37, no. 9, pp. 1473 – 1479,
2001.
[10] J.-C. Shen, “New tuning method for PID controller,”
ISA Transactions, vol. 41, no. 4, pp. 473 – 484, 2002.
[11] G. Syrcos and I. Kookos, “Pid controller tuning using
mathematical programming,” Chemical Engineering and
Processing: Process Intensification, vol. 44, no. 1, pp.
41–49, 2005, cited By 15.
[12] M. Ramasamy and S. Sundaramoorthy, “Pid controller
tuning for desired closed-loop responses for siso systems
using impulse response,” Computers and Chemical
Engineering, vol. 32, no. 8, pp. 1773–1788, 2008, cited
By 21.
[13] N.-S. Pai, S.-C. Chang, and C.-T. Huang, “Tuning pi/pid
controllers for integrating processes with deadtime and
inverse response by simple calculations,” Journal of
Process Control, vol. 20, no. 6, pp. 726 – 733, 2010.
[14] A. O. Dwyer, Handbook of PI and PID controller tuning
rules. London, U.K.: Imperial College Press, 2006.
[15] D. Seborg, T. Edgar, and D. Mellichamp, Process Dy-
namics and Control. John Wiley and Sons, 2004.
[16] C. Smith, “Pid explained for process engineers: Part 2
- tuning coefficients,” Chemical Engineering Progress,
vol. 112, no. 2, 2016, cited By 0.
[17] J. Guzmán, T. H´’agglund, M. Veronesi, and A. Visioli,
“Performance indices for feedforward control,” Journal
of Process Control, vol. 26, pp. 26 – 34, 2015.
[18] S. Das, I. Pan, S. Das, and A. Gupta, “A novel fractional
order fuzzy PID controller and its optimal time
domain tuning based on integral performance indices,”
Engineering Applications of Artificial Intelligence,
vol. 25, no. 2, pp. 430 – 442, 2012, special Section:
Local Search Algorithms for Real-World Scheduling
and Planning.
[19] T. Olsen and N. Ito, “Implement an effective loop tuning
strategy,” Chemical Engineering Progress, pp. 42–48,
2013.
[20] F. Haugen, “The good gain method for simple experimen-
tal tuning of PI controllers,” Modeling, Identification and
Control,, vol. 33, pp. 141–152, 2012.