Ceic3006 Lecture 5
Ceic3006 Lecture 5
CEIC3006
Process Dynamics and Control
Week 5
Lecturer: Prof. Jie Bao
This week’s lectures
Tin, F
F,Tin Controller Valve
TT TC
e u ws T
+ C A P
IP
TR -
ws Tank
M
Steam
Condensate
F,T
Thermocouple
In the steam heated stirred tank example, the controller decides the pressure of the
compressed air p(t) to manipulate the valve opening. Deriving the transfer function for an
ideal proportional controller
Actuator:
physical (with dynamics) process triggered by the controller
directly affects process
Sensor:
monitors some property of the system and transmits signals back to
the controller
Control valve
Control valves are usually designed so that the flow rate through the
valve is a nearly linear function of the signal to the valve actuator
resulting a first-order transfer function:
Sensors/Transmitters
A sensor’ transfer function can be determined from its dynamic model,
e.g., a thermal couple
GD
Controller Valve Tank
+ P Ws T
e Gv GP
- GC
TR
Tm GM
Thermocouple Analysis of dynamics of standalone
The overall closed-loop systems systems can be applied to study:
Tin T
Gyd(s)
CEIC3006 Session 2, 2022
Closed-Loop Transfer Function
D
1. Closed-Loop Servo Response: Eliminate all variables except R(s) and Y(s) Thermocouple
Isolate Y(s)
CEIC3006 Session 2, 2022
Closed-Loop Transfer Function
D
Thermocouple
Isolating Y(s)
Servo Response:
Servo
Regulatory
CEIC3006 Session 2, 2022
Feedback Control
Steam heated tank
F,Tin TT TC
IP
Ws
Steam
Condensate
F,T
Feedback control system: Valve is manipulated to increase flow of steam to control tank
temperature
Closed-Loop
Time Constant
Larger Kc, quicker closed-loop response
CEIC3006 Session 2, 2022 18
Proportional Feedback
Final response:
Note:
for “zero offset response” we require
Tracking Disturbance rejection
Possible to eliminate offset with P-only feedback (requires infinite controller gain)
Need different control action to eliminate offset (integral)
Proportional control is often insufficient!
Proportional, Integral and Derivative Control (PID control)
CEIC3006 Session 2, 2022
19
Integral Control
Integrator is included to eliminate offset
provides reset action: by summing the historical errors
usually added to a proportional controller to produce a PI controller
• PI is the most widely used controller in industry
• optimal structure for first order processes
PI controller form
Closed-loop response
Note:
• Offset is removed, closed-loop is second order
• damping coefficient
integral time constant and controller gain can induce oscillation and change the
period of oscillation
CEIC3006 Session 2, 2022 22
PI Feedback
Effect of controller gain
Increasing controller gain:
Increases speed of response
Reduces control error quicker
May cause oscillatory response or
causes instability
PID Controllers:
• greater than 90% of all control implementations
• dates back to the 1930s
• very well studied and understood
• optimal structure for first and second order processes (given some assumptions)
• always first choice when designing a control system
Integral action
Or alternatively,
Note:
denominator provides integration to remove possibility
of steady-state errors
τ Dτ I s2 + τ I s + 1
Y ( s) = R( s) +
τ Iτ 2 1 + Kc K p
+ τ Dτ I s + τ I s + 1
K K
c p K K
c p
τ Iτ 2 τ I
s + s
K K
c p K K
c p
D( s)
τ Iτ 2 1 + Kc K p
+ τ Dτ I s + τ I s + 1
Kc K p Kc K p
In Matlab:
G = parallel (G1,G2)
Series Connection
In Matlab:
G = series (G1,G2)
Feedback Connection
In Matlab:
Gcl = feedback (G,H)
CEIC3006 Session 2, 2022
Feedback systems
General Expressions
Negative Feedback Control System: the overall controlled system, the
closed-loop system that includes the controller, actuator, process and sensor:
Where GCL(s) is the closed-loop system, Gf(s) is the feedforward system and GOL(s) is the
open-loop system
1.2
1
0.8
0.6
0.4
0.2
0 0 2 4 6 8 10 12
Time (sec.)
25
Step Response
20
15
Amplitude
10
5
0
-5
-10
-15
-20
0 10 20 30 40 50 60
Time (sec.)
12000
Step Response
10000
Amplitude
8000
6000
4000
2000
-2000
0 1 2 3 4 5 6 7 8 9 10
Time (sec.)
TR T
Gcl(s)
Stability analysis of dynamics of standalone
systems can be applied to determine the
Tin T stability of the closed-loop system!
Gyd(s)
CEIC3006 Session 2, 2022
Closed-loop Stability
D(s)
Closed-loop dynamics Gd
controller actuator process
R(s)+ Y(s)
Gc Gv Gp + +
-
sensor
Gm
GOL
Characteristic Equation:
1 + GOL(s) =0
The roots of the characteristic equation is the poles of the closed-loop system.
Unstable region is the right half plane (RHP) of the complex plane.
For nonlinear system, if its linearized system is stable, the stability holds only
locally. Moving away from the point of linearization may cause instability.
CEIC3006 Session 2, 2022
Closed-loop Stability
Root locations
Example
Roots:
Kc=2: -1.25, -0.22 + 0.44i, -0.22 - 0.44i
Kc=6: -1.48, -0.11 + 0.68i, -0.11- 0.68i
Kc=16: -1.79, 0.04 + 0.97i, 0.04 -0.97i
Step Response
3
2.5
2
Kc=16
1.5
Amplitude
0.5
Kc=6 Kc=2
0
-0.5
-1
-1.5
0 2 4 6 8 10 12 14 16 18 20
Time (sec.)
CEIC3006 Session 2, 2022
Example
To have a stable system, both roots of this characteristic equation must have negative real parts.
The roots will be negative if 40(Kc + 1) > 0.
The square root will have a value less than 7.
If 40(Kc + 1) > 0, then Kc + 1 > 0, and Kc > −1.
The closed-loop system will be stable if Kc > −1.
Analysis problems
For known closed-loop systems
1. Known process model
2. Known controller
MATLAB function ROOTS
Design problems
Unknown controller parameters
Routh array
1 an an − 2 an − 4
2 an − 1 an − 3 an − 5
3 b1 b2 b3
4 c1 c2
n +1 z1
a a −a a a a −a a
b1 = n −1 n − 2 n − 3 n , b2 = n −1 n − 4 n −5 n ,
an −1 an −1
where ba −b a ba −b a
c1 = 1 n − 3 2 n −1 , c2 = 1 n −5 3 n −1 ,
b1 b1
All elements of left column must be positive to have roots with negative real parts
Unit-step disturbance responses for the candidate controllers (First order plus time delay
(FOPTD) Model: K = 1, θ = 4, τ = 20) and G = Gd
More aggressive
Unstable
More aggressive
Best
performance
Controller settings can be adjusted to achieve the desired closed-loop system performance, a
procedure referred to as controller tuning.
CEIC3006 Session 2, 2022
Controller Design
Control design criteria
Closed-loop system is stable.
Offset is eliminated
sensor
Gm
CEIC3006 Session 2, 2022
Controller Design Methods
G
process d
R + controller actuator + C
Gc Gv Gp +
-
sensor
Gm
Isolate :
~ , controller is given by
For a desired trajectory (C/R)d and plant model G
Perfect Control
cannot be achieved, requires infinite gain
e.g.
Leads to a PI controller
Problem: the dead time makes analysis (poles and zeros) more difficult
Desired closed-loop:
Leads to a PI controller
Process parameters:
𝜏𝜏𝑐𝑐 = 12
20 2
15 1.5
y(t)
y(t)
10 1
5 0.5
0 0
0 100 200 300 0 100 200 300
Time (seconds) Time (seconds)
Consider three values of the desired closed-loop time constant: τc = 1, 3, and 10. Evaluate
the controllers for unit step changes in both the set point and the disturbance, assuming
that Gd = G. Repeat the evaluation for two cases:
a. The process model is perfect (𝐺𝐺� = 𝐺𝐺).
� = 0.9 instead of the actual value, K = 2.
b. The model gain is incorrect, 𝐾𝐾
τc = 1 τc = 3 τc = 10
� = 2)
𝐾𝐾𝑐𝑐 (𝐾𝐾 3.75 1.88 0.68
� = 0.9)
𝐾𝐾𝑐𝑐 (𝐾𝐾 8.33 4.17 1.51
τI 15 15 15
τD 3.33 3.33 3.33
τc = 1 τc = 3 τc = 10
� = 2)
𝐾𝐾𝑐𝑐 (𝐾𝐾 3.75 1.88 0.68
� = 0.9)
𝐾𝐾𝑐𝑐 (𝐾𝐾 8.33 4.17 1.51
τI 15 15 15
τD 3.33 3.33 3.33
The process information may not be accurate - field tuning may still be required to
fine tune the controller
Advantages:
Process models are not required
Disadvantage:
Tedious and time consuming
Hard to get the best controller parameters
Controller Kc τI τD
P 0.5Kcu — —
PI 0.45Kcu Pu/1.2 —
PID 0.6Kcu Pu/2 Pu/8
Remarks
Widely considered to be an "industry standard"
Advantages:
Only need one trial-and-error search
Large safety margin for Kc
Disadvantages:
Slow dynamics results the trial-and-error determination of Kcu and Pu be time consuming.
For many applications, continuous cycling is objectionable because the process is pushed to a stability limit.
CEIC3006 Session 2, 2022
Continuous Cycling Method
The Z-N controller settings are based on a 1/4 decay ratio.
The settings tend to produce oscillatory responses and large overshoots for
set-point changes.
Modification
Modified Ziegler-Nichols Settings for PID
Controller Kc τI τD
Original (1/4 decay ratio) 0.6Kcu Pu/2 Pu/8
Some overshoot 0.33Kcu Pu/2 Pu/3
No overshoot 0.2Kcu Pu/2 Pu/3
Example
Proportional
+
( θ)
KK C = 0.9 τ KK = 0.9(τ )+ 0.083
C θ
τ θ [3.33 + 0.33( τ )]
θ
= 3.33(θ )
Integral τ I
1.0 + 2.2(θ )
I
=
τ τ τ
τ
Proportional KK = 1.2(τ ) ( θ )+ 0.270
KK C = 1.35 τ
+ C θ
θ [32 + 6(θ τ )]
= 2.0(θ )
Integral τ τ =
13 + 8(θ )
I
+ τ I
Derivative
τ τ
τ
= 0.5(θ )
D τ 0.37(θ )
τ
τ τ D
τ 1.0 + 0.2(θ )
=
τ
Example:
Reset Windup
Sustained error causes the integral term to become very large and the
PID controller output eventually saturated.
Anti-reset windup
• Temporary halting the integral action whenever the controller output
saturates.
Physically unrealizable:
Needs future information
Proper systems
A transfer function whose order of numerator is no greater (or less) than its
order of denominator is a proper (or strictly proper) system.
Non-proper systems cannot be realized
Sudden changes in set point, thus e(t), will cause very large derivative term.
Eliminating the derivative kick
Take the derivative action on the measurement ym(t) rather than the error signal e(t).