Colorado School of Mines CHEN403 Ziegler-Nichols Example

Ziegler-Nichols Controller Tuning Example

The Ziegler-Nichols method uses a closed controller loop & requires the following steps:
Bring system to steady state operation.
Put on P control. Introduce a set point change and vary gain until system
oscillates continuously. This frequency is CO and M is the amplitude ratio.
Compute the following:
1
Ultimate Gain K u
M
2
Ultimate Period Pu
CO

The original Z-N tuning settings are given in the following table.
Controller Kc I D
P K u /2
PI K u /2.2 Pu /1.2
PID K u /1.7 Pu /2 Pu /8

These controller settings were developed to give a 1/4 decay ratio. However, other
settings have been recommended that are closer to critically damped control (so that
oscillations do not propagate downstream). These PID controller settings are shown in the
following table.
Controller Kc I D
Original 0.6K u Pu /2 Pu /8
Some Overshoot 0.33K u Pu /2 Pu /3
No Overshoot 0.22K u Pu /2 Pu /3

The question becomes one of how to find these parameters when we dont have a real
process but rather the idealized transfer functions for the process.

John Jechura ([email protected]) -1- Copyright 2017

April 23, 2017
Colorado School of Mines CHEN403 Ziegler-Nichols Example

General 3rd Order Overdamped System Example

Kp
R + Gc + Gp Y
1s 1 2 s 1 3 s 1
- +

As an example, let's assume we have a 3rd order overdamped system with negligible
dynamics in the final control & measurement elements. Then:

Kp
Gp s
1 s 1 2s 1 3s 1
G f s Gm s 1

The transfer function between the set point and the output under P control will be:

Kp
Kc
Gs
Gc Gp

1 s 1 2 s 1 3 s 1
1 Gc Gp Kp
1 Kc
1 s 1 2 s 1 3 s 1
KcK p

1 23 s 1 2 23 31 s 2 1 2 3 s 1 K p K c
3

The characteristic equation of this transfer function:

1 23 s3 12 23 31 s 2 1 2 3 s 1 K pK c 0

will control the stability & oscillatory nature of the systems response. Undamped
oscillatory behavior occurs at the very edge of the systems stability limit.

Routh Array. We can use a Routh array to determine this point.

John Jechura ([email protected]) -2- Copyright 2017

April 23, 2017
Colorado School of Mines CHEN403 Ziegler-Nichols Example

Row
1 1 23 1 2 3
2 12 23 31 1 K pK c
3 1 23 1 K p K c
1 2 3
1 2 23 31
4 1 K pK c

The stability limit from the 4th row is:

1
1 K pK c 0 K c
Kp

which will be satisfied for any positive K c . The stability limit from the 3rd row is:

1 23 1 K p K c 1 2 3 12 23 31 1
1 2 3 0 Kc
1 2 23 31 K p 123 Kp
1 1 1 1 1
Kc 1 2 3
Kp 1 2 3 K p

So, the ultimate gain will correspond to the stability limit:

1 1 1 1 1
K cu 1 2 3
Kp 1 2 3 K p

The period of oscillation at this ultimate gain must be calculated from the inverted Laplace
expression.

Direct Substitution. We could also use direct substitution to determine the stability limit for
this process. Substituting s j into the characteristic equation:

123 j 12 23 31 j 1 2 3 j 1 K pK c 0
3 2

1 23 3 j 12 23 31 2 1 2 3 j 1 K pK c 0
1 K K 23 31 2 1 2 3 1233 j 0 .
p c 1 2

Setting the imaginary part to zero gives us the frequency of oscillation:

John Jechura ([email protected]) -3- Copyright 2017

April 23, 2017
Colorado School of Mines CHEN403 Ziegler-Nichols Example

1 2 3
1 2 3 u 1233u 0 u
1 23

and setting the real part to zero gives us the controller gain at the stability limit:

1 2 3
12 23 31 1
1 23
1 K K
p cu 1 2 23 31 2u 0 K cu
Kp
1 1 1 1 1
1 2 3
K p 1 2 3 Kp

(which is exactly the same as that given by the Routh array analysis).

Specific 3rd Order Overdamped System Example

Let us take for example the process:

1
Gp s
s 15s 10.2s 1

The ultimate gain will be:

1 1 1 1 1
K cu 1 5 0.2 37.44
1 1 5 0.2 1

and the frequency of oscillation will be:

1 5 0.2
CO 2.48998
1 5 0.2

so the period of oscillation at the ultimate gain is:

2
Pu 2.52339
CO

The following will be the Ziegler-Nichols controller settings:

P control:

John Jechura ([email protected]) -4- Copyright 2017

April 23, 2017
Colorado School of Mines CHEN403 Ziegler-Nichols Example

Ku
Kc 18.72
2

PI control:
Ku
Kc 17.02
2.2
P
I u 2.10
1.2

PID control:
Ku
Kc 22.02
1.7
P
I u 1.26
2
P
D u 0.32
8

Lets look at a step change to the load for PI control. The closed loop transfer function will
be:

Gs
Gp

s 15s 10.2s 1
1 Gc Gp 1 1
1 Kc 1
I s s 15s 1 0.2s 1
I s

I s s 15s 1 0.2s 1 K c I s 1
I s

I s 6.2I s 6.2I s 2 I 1 K c s K c
4 3

With a step change to the load, U 1/ s and the Ziegler-Nichols settings:

1
Y s
8.1039 18.0182s 6.2s2 6.2s3 s 4
Y t 0.00573006e 5.62165t 0.0626769e 0.495759t
e 0.0412961t 0.00204812sin 1.70472t 0.0569469cos 1.70472t

John Jechura ([email protected]) -5- Copyright 2017

April 23, 2017
Colorado School of Mines CHEN403 Ziegler-Nichols Example

0.10
0.08
0.06
0.04
y(t) 0.02
0.00
-0.02
-0.04
-0.06
0 5 10 15 20
t

This figure shows the response to this disturbance. The controller settings are very
conservative there is very little damping to the response. In actual practice, these
settings would only be used as initial settings and additional tuning would be done.

John Jechura ([email protected]) -6- Copyright 2017

April 23, 2017