Control Systems: Fundamentals & Modelling
Control Systems: Fundamentals & Modelling
Handout Part 1:
Summary:
This handout sits alongside the slides (it has the key info in and some spaces for you to
complete as we go). It introduces the control system (along with the Lecture 1 handout) and
sets the scene for why we design closed-loop controllers. It reminds you of some of the
basic math (section2) – we will not cover this in the lectures, just use it later. Then we’ll be
discussing how to build mathematical models of engineering systems, together with
examples. These are the foundation for analysis upon which control systems design is built.
Control Systems Part 1: Control Fundamentals & System Modelling
• The typical closed-loop control system and its benefits (over open-loop)
• The basic aims of the control engineer in designing a control system
• The basic (Time Domain) methods of assessment of performance
A list of commonly used control engineering terms is also included at the end of the section.
Control Aim:
To make the system output “track” or follow the command as closely as possible
In open-loop control we have a feed-forward controller which is typically the inverse of the
plant. The problem is it does not cope well with disturbances or uncertainty:
Disturbances –
Unknown/unmeasured inputs from outside the system being controlled which affect its
operation .
Uncertainties
E.g. variations in the weight of a car affecting the cruise control system
R Dixon Page 2 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Disturbances
Uncertainty
Controller Input
signal δθ
Command Error u
+ Output
C(s) System y
yin
- Compensator
“Feedback”
Control Aims:
To make the system output “track” or follow the command as closely as possible
In closed-loop control, the output (y) is fed-back to the controller which acts to minimise the
error (i.e. the difference) between the command and the output.
Exercise
Can we imagine the difference in an overall system response in open-loop and closed loop
scenarios, when a change is made to the system?
For a cruise control application, assume a 1 degree change in throttle position gives a
5km/hr change in speed. Also a 1% change in slope (disturbance) gives a -2mph change in
speed.
Disturbance,
0%
Car
Controller -2
(km/hr) / %
+
1/5 8deg 5 +
40km/hr deg/(km/hr)
40km/hr
(km/hr)/deg
R Dixon Page 3 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Disturbance,
0%
Car
Controller -2
(km/hr) / %
1 7.8deg
40km/hr 7.8 5 39 km/hr
+
deg/(km/hr) (km/hr)/deg
-
In the closed-loop you will see that we start with an error of 1km/hr. But after the disturbance
it settles to only 1.25km/hr. Compare with the open-loop (10km/hr).
Note: We have not yet covered the block-diagram math for you to work out (easily) the
steady state of the closed-loop system. But hopefully you get the idea – you can see what
happens for various values of the error.
Sensors
Controller System
Command Output
The system is the engineering system to be controlled. In general the aim is for the output
of the system to track (follow) the command of the system.
R Dixon Page 4 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
1.2.1 Tracking
Here, the primary requirement is for the system to follow a continuously-varying command
input (although it will also have to react to disturbances). E.g. position control of a robot. In
that case we seek to make the output y follow command yin with a particular speed of
response and a particular steady state error, whilst not using too much control action, u.
1.2.2 Regulation
Now, the primary requirement is for the system to maintain its output in response to external
disturbances. Although there may also be occasional variations in the “setpoint”/command.
E.g. temperature control of a central heating system.
In practice many control systems are a combination of the regulation and the tracking
problem – this involves a trade-off in the design of the controller (between tracking and
regulation responses).
R Dixon Page 5 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
The Laplace Transforms of these most commonly used inputs can be found in the table of
Laplace Transforms.
A response to a ‘step-input’ can be seen in the graph below. The output ‘y’ is plotted
alongside the input command ‘yin’.
R Dixon Page 6 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
• The Transient part where the output is changing from its start position to its new
state of equilibrium
• The Steady State where the new state of equilibrium has been reached.
Understanding key attributes within each of these conditions allows for careful design of
controllers. One of the key items that we can affect by our design decisions as control
engineers is the control gain within C(s). Adjusting this will change the way the system
responds to changing inputs and the control compensator is designed until a satisfactory
response is achieved in the transient and steady state parts of the response. A typical effect
of changing the control gain is shown below:
Question: what are the main differences in the time response between high and low
controller gains? Note your thoughts here. Think about the transient and steady state parts
of the response.
R Dixon Page 7 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
R Dixon Page 8 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
From a time response it is possible to assess the absolute stability as in it is possible to see
that either the response:
Exercise
Draw a block diagram of a closed-loop system which has all these elements included and
labelled
R Dixon Page 9 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
A complex number is a number that can be put in the form a+bj, where a and b are real
numbers and j is the imaginary unit (the square root of −1).
In this expression, a is called the real part and b the imaginary part of the complex number.
Complex numbers extend the one-dimensional number to a two-dimensional complex-plane
by using the horizontal axis for the real part and the vertical axis for the imaginary part.
In the figure below, the complex number a + bj can be identified with the point (a, b).
Alternatively a complex number can also be expressed in terms of polar co-ordinates.
Cartesian Co-ordinates:
x = a + bj Im a x
Polar Co-ordinates:
R
x = 𝑅𝑅𝑒𝑒 𝑗𝑗𝑗𝑗
ϕ b
= 𝑅𝑅 (𝑐𝑐𝑐𝑐𝑐𝑐(𝜑𝜑) + 𝑗𝑗. 𝑠𝑠𝑠𝑠𝑠𝑠(𝜑𝜑)) Re
Gain (magnitude): R =
Angle (argument): j =
B B
The Gain (R) is given by: =
A A
B
And the Angle (𝜙𝜙) by: ∠ = ∠B − ∠A
A
R Dixon Page 10 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
The formal mathematical definitions, properties and theorems are quite involved (see your
math notes!). In this module we will focus on key aspects relevant for control systems
analysis & design – these are quite straightforward with a little practice. First a reminder of
some of the math….
∞
F ( s ) = L{ f (t )} = ∫
0
f (t )e − st dt
This transform is also known as the one-sided Laplace Transform (LT). The Laplace
operator, s, is a complex number:
𝑠𝑠 = 𝜎𝜎 + 𝑗𝑗𝑗𝑗
f (t ) = L-1{F ( s )}
The commonly useful transforms and their inverses are listed in LT tables available in control
text books. There are also a number of Theorems which are useful to control engineers.
R Dixon Page 11 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
n!
3 tn
s n +1
1
4 e − at
s+a
ω
5 sin(ωt )
s +ω2
2
s
6 cos(ωt )
s +ω2
2
ω
7 e − at sin(ωt )
(s + a )2 + ω 2
(s + a )
8 e − at cos(ωt )
(s + a )2 + ω 2
In the table above, all the functions f(t) are defined for t ≥ 0. For example, u(t) is unit step
function which has amplitude 0 for t < 0 and unity (1) for t ≥ 0.
The table above can be used for Laplace Transforms and for Inverse Laplace Transforms. In
other words, we can go from f(t) to F(s) and vice-versa.
R Dixon Page 12 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Some of the key theorems and properties associated with Laplace transforms can be found
below…
R Dixon Page 13 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
We can use the LT to re-express time domain expressions using the Laplace operator ‘s’.
This allows dynamic systems described by differential equations to be represented in a form
which can be manipulated and handled algebraically to generate transfer functions.
The use of the Laplace operator is also helpful in forming block diagram representations of
dynamic systems (we will see soon how useful this can be).
There are important theoretical techniques in the frequency domain which arise from
considering ‘s’ as a complex number. These and allow us to use our models to analyse
systems in the frequency and the time domains – and through doing so design control
systems.
Time Domain
Differential Solution in t
Equations (time response)
Algebraic
Solution in s
Equations
Frequency Domain
(Frequency response)
Both methods can produce transfer functions relating outputs to inputs and both are
mathematically equivalent so they will give the same answer for a given system.
R Dixon Page 14 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
𝑌𝑌(𝑠𝑠) 1
𝐺𝐺(𝑠𝑠) = = 2
𝑈𝑈(𝑠𝑠) 𝑠𝑠 + 4𝑠𝑠 + 2
Where ‘s’ is the Laplace Operator (as discussed earlier), and G(s) is used to represent the
ratio of Y(s) to U(s).
3.1.1 Stability
The behaviour of the overall time response of the transfer function can be identified by
calculating the position of the poles on the complex s-plane (below). The poles and zeros of
a transfer function can be found by equating the polynomial expressions for the denominator
and numerator (respectively) to zero, and solving to find the values of ‘s’ for which this
occurs.
Typically we plot these values of s on the complex s-plane (as per below) with its imaginary
and real parts. Looking at the diagram (and the skull and cross-bone area)….
For a system to be stable, the real part of the pole should -ve, and the pole should therefore
sit on the left hand side of the s-plane.
R Dixon Page 15 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Question: Take a look at line 4 of Table 1-1 (a first order system with a single, non-complex
pole that occurs when s = -a). Considering the inverse Laplace transform of the transfer
function, can you see why ‘a’ must be a positive value for stability?
• Stable - if the natural output response decays to zero as time tends to infinity. I.e. all
the poles have negative real parts.
• Marginally stable - if the natural response oscillates with constant bounds (peaks)
as time tends to infinity. Or if it increases indefinitely at a constant rate. I.e. one or
more of the poles has zero real part(s) and all others have negative real parts.
• Unstable - if the natural response grows without bounds as time approaches infinity.
I.e. one or more of the poles has positive real parts.
The following graphs show examples of stable and unstable systems for ‘first order’ systems
(that have a single pole) and second order systems (that have 2 poles).
Unstable
R Dixon Page 16 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
For example, given the arbitrary second order system described by the differential equation:
𝑌𝑌(𝑠𝑠) 1
= 2
𝑈𝑈(𝑠𝑠) 𝑠𝑠 + 4𝑠𝑠 + 2
Transfer functions are useful to relate single outputs to single inputs, but intermediary
information can be hard to obtain, and more complex systems can become algebraically
difficult to work with. Most of the systems that will be used in this module will be single-input-
single-output systems, but it is always worth having other approaches on hand if required.
A comment on Laplace Transforms for control engineering (especially for this module):
When taking Laplace transforms during this course, it is possible (and desirable) to assume zero
initial conditions – which means taking the Laplace transform is as simple as replacing d/dt
with ‘s’. So
R Dixon Page 17 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Block diagrams are often constructed by considering the differential system equations and
building the system up using component parts. The basic building blocks will first be
introduced, then an example system will be constructed.
G G
Gain: x(t) y(t) X(s) Y(s)
Differentiator d/dt s
x(t) y(t)
X(s) Y(s)
Integrator ∫ (.) dt
x(t) y(t)
These blocks can be connected in numerous different ways and combinations. For example,
this is a feedback network:
+
X(s) A(s) Y(s)
-
B(s)
Note that A(s) and B(s) are transfer functions, and a ‘sum block’ is used to combine signals.
R Dixon Page 18 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Blocks in series:
U (s ) Y (s ) U (s ) Y (s )
B(s) A(s) B( s ) A( s )
U (s ) + Y (s )
U (s ) Y (s )
A(s ) A( s )
- 1 + A( s ) B( s )
B (s )
The last transformation is probably the most useful one to remember in this series as the first
is relatively straight forward, and the second isn’t too common.
The best way to start building a block diagram is to re-arrange the equation so that the most
differentiated output variable is the subject of the equation is on the left hand side.
Start drawing the diagram by considering 𝑦𝑦̈ (𝑡𝑡) and what the equation above represents. To
calculate 𝑦𝑦̈ (𝑡𝑡), we need the input 𝑢𝑢(𝑡𝑡) summed with 𝑦𝑦̇ (𝑡𝑡) (and some coefficient) and 𝑦𝑦(𝑡𝑡)
(again with a coefficient). How is 𝑦𝑦̇ (𝑡𝑡) and 𝑦𝑦(𝑡𝑡) created from the starting point 𝑦𝑦̈ (𝑡𝑡)using
blocks? How do we include a coefficient? Start drawing the required elements over the
page...
R Dixon Page 19 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
𝑦𝑦̈ (𝑡𝑡)
Now the diagram can be converted to a transfer function using the ‘simplify feedback loop
transformation in Section 3.2.2. (hint: Start with the ‘inner loop’ then do the outer loop’)
R Dixon Page 20 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
u y y y
1 1
+ +
- - s s
4
2
u 1 y
s + 4s + 2
2
...which gives the same answer as the algebraic approach used in Section 3.1.2.
An advantage to this approach is that a block diagram exposes more information than a
single transfer function. These diagrams are often used in conjunction with modelling
software (such as Simulink). When these simulations are executed in software, access to not
only input and output variables is possible (𝑢𝑢(𝑡𝑡) and 𝑦𝑦(𝑡𝑡)) but to other internal variables,
such as 𝑦𝑦̇ (𝑡𝑡) and 𝑦𝑦̈ (𝑡𝑡) or even other items such as ‘error’ signals or ‘drive’ signals. In the
algebraic world, separate transfer functions would be required for each of these.
Furthermore, outside of the classroom, if you are asked to model a dynamic system as part
of your job doing it two different ways gives (some) confidence that you got the right answer.
- Always assuming you did get the same result from both approaches!!
R Dixon Page 21 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
4 Modelling Examples
Use these to practise creating systems models. Don’t forget that there are plenty more
exercises to try in the recommended text books.
x (m)
Damper constant
b=2(N/m/s)
Mass, u , Applied
M =10kg Force (N)
Spring constant
k=5(N/m)
2. Equations of Motion
R Dixon Page 22 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
R Dixon Page 23 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
4.2 A DC Motor
The direct current (DC) motor was one of the first machines developed to convert electrical
power into mechanical power. Its history can be traced back to Michael Faraday who
experimented with a number of rotating electromechanical machines in the first part of the
19th Century. The physical principle which produces a mechanical force from the
interactions of an electric current and a magnetic field is known as “Faraday's law of
induction”, which was discovered in 1831.
Today, DC Motors come in all shapes and sizes and are used in many applications: e.g.,
fans/blowers and pumps, machine tools, conveyor belts, household appliances, power tools,
disk drives, electric/hybrid vehicles, ship propulsion etc.
The motor is expressed by two equations, one describing the electrical part of the system,
and another for the mechanical.
R Dixon Page 24 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Variables Parameters
V Voltage across (Volts) R Armature resistance (Ohms)
windings
I Armature Current (Amperes) L Armature inductance (Henries)
Ve Back Emf (Volts) Ke Back emf Constant V/rads-1
𝜃𝜃̇ Rotational velocity rad/s
Variables Parameters
𝜃𝜃̈ Rotational rad/s2 J Armature inertia kgm2
acceleration
Te Induced torque from Nm Tl Motor Torque load Nm
motor (assumed zero)
Tf Torque due to Nm D Friction Coefficient Nm/rads-1
friction
The equation is formed considering the sum of torques (the rotational form of the Newton’s
second law):
R Dixon Page 25 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
The first step is to take the laplace transforms of the equations of interest:
𝜃𝜃(𝑠𝑠) 𝐾𝐾𝑡𝑡
Should give….. =
𝑉𝑉(𝑠𝑠) 𝑠𝑠(𝐿𝐿𝐿𝐿𝑠𝑠 2 +(𝑅𝑅𝑅𝑅+𝐿𝐿𝐿𝐿)𝑠𝑠+𝐾𝐾𝑡𝑡 𝐾𝐾𝑒𝑒 +𝐷𝐷𝐷𝐷)
Furthermore, we can make use of the fact that 𝜃𝜃̇ (𝑠𝑠) = 𝑠𝑠𝑠𝑠(𝑠𝑠). This can give us a transfer
function relating output velocity to applied voltage quite easily by a quick substitution:
R Dixon Page 26 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Suggested procedure:
R Dixon Page 27 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Answer:
D
R
Ke
There are some common assumptions made to simplify our DC motor model.
θ( s ) Kt
=
V (s) LJs + ( RJ + LD ) s + K t K e + DR
2
Often the dynamics of the electrical circuit are much faster than those of the mechanical
system. This can lead to an assumption that L is very small and negligible. I.e. L=0. In this
situation we end up with a first order model of the motor. Also it is quite common to assume
that D is approximately zero.
θ( s ) Kt
=
V (s) RJs + K t K e
1/ K e
=
( RJ / K t K e ) s + 1
K
Namely, G (s) = where K = 1 / K e and τ = RJ / K t K e
τs + 1
In this case, K is the steady-state gain of the system and is the time taken for the system to
go from 0 to 63% of the final value. The mathematical reasons behind this will be clear later
in the module. However, it is useful in the case where an experimental step response of the
system can be obtained because we can then directly measure K and τ .
Motor Exercise – build the model in Simulink, run a step response and check K and τ .
Use parameters:
-1
Back e.m.f. constant Ke = 0.024 V/rads
Torque constant KT = 0.024 Nm/A
Armature resistance R = 4Ω
-5 2
Movement of inertia J = 1.1 x 10 kgm
R Dixon Page 28 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
4.3 Car Model (for Cruise Control) - this section is for you to look at in your time.
In order to model our car for the purposes of designing a cruise control system we need to
make some assumptions:
1. That the engine has a torque controller which can arrange for the car to transmit a
given (accelerating) force to the road.
2. That the drag due to road friction and air resistance can be assumed proportional to
velocity of the vehicle.
3. That we can ignore the (relatively fast) dynamics in the engine and drive-train.
With these assumptions in mind we can draw a free body diagram as follows:
Where:
And:
Applying Newton’s second law to the car free body diagram, we have:
𝑢𝑢 − 𝑏𝑏𝑏𝑏 = 𝑚𝑚𝑣𝑣̇
R Dixon Page 29 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Assume zero initial conditions, take the Laplace transform and show that the output/input
transfer function is given by
V(s) 1
=
U(s) ms + b
R Dixon Page 30 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
2. Numerical Solution
Now we can put in the numbers to obtain the actual numeric transfer function:
V(s) 1
=
U(s) 1080s + 108
As a TF as block diagram
Then, we could do some scaling to get the input and output in more normal units given the
size of a car → [input force in kN and output speed in km/hr]
A straightforward way to do this is using the block diagram. So we define a new input u’(kN)
and a new output v’(km/hr).
Or,
V' (s) 10
=
U' (s) 3s + 0.3
R Dixon Page 31 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Quick Exercise
Now compare with the standard first order system introduced in DC motor example (A4,
Ex2),
K
G (s) =
τs + 1
What is the steady state gain K for the motor car and what is the time constant τ?
K =………..
τ =………..
See how these compare to the response to a step input as produced by Matlab/Simulink
(below).
R Dixon Page 32 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
Description:
The exercise is concerned with the tank system shown in Figure 1. The requirement is to
calculate its transfer function. Then to add proportional level controller (with a gain of G) and
calculate the closed-loop transfer function (this would allow us to begin to characterise its
dynamic performance).
h
qout
Variable Parameters
s
q the control input (a m3s-1 c the (linearised) 0.1 m2s-1
flow of fluid) outlet coefficient
h the level of the m A Cross sectional 0.5 m2
liquid area of the tank
qout is the flow out of the 𝑞𝑞𝑜𝑜𝑜𝑜𝑜𝑜 = ℎ𝑐𝑐
tank
2. Using the algebraic approach take the Laplace Transform assuming zero initial
conditions and write down the TF
3. Use the block diagram approach to model the tank level – then simplify to get the TF.
4. Draw a block diagram showing the system model with a proportional controller with
gain, G, applied to the system. Write down the closed-loop transfer function.
R Dixon Page 33 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
R Dixon Page 34 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
R Dixon Page 35 of 36
Control Systems Part 1: Control Fundamentals & System Modelling
n!
3 tn
s n +1
1
4 e − at
s+a
ω
5 sin(ωt )
s +ω2
2
s
6 cos(ωt )
s +ω2
2
ω
7 e − at sin(ωt )
(s + a )2 + ω 2
(s + a )
8 e − at cos(ωt )
(s + a )2 + ω 2
R Dixon Page 36 of 36