Sie sind auf Seite 1von 28

# INSTRUMENTATION AND CONTROL

## TUTORIAL 1 CREATING MODELS OF ENGINEERING SYSTEMS

This tutorial is of interest to any student studying control systems and in particular the EC
module D227 Control System Engineering. The purpose of this tutorial is to introduce
students to the basic elements of engineering systems and how to create a transfer function
for them. The tutorial is mainly informative and consists of examples showing the derivation
of models for real hardware systems. The self assessment material is based on basic general
engineering knowledge.
On completion of this tutorial, you should be able to do the following.
Derive the mathematical models of basic mechanical systems.
Derive the mathematical models of basic fluid power systems.
Derive the mathematical models of basic thermal systems.
Derive the mathematical models of basic electrical systems.
Recognise the similarity between models of different systems.
Explain the standard first and second order transfer functions.
Explain the link between open and closed loop transfer functions.
If you are not familiar with instrumentation used in control engineering, you should complete
the tutorials on Instrumentation Systems.
In order to complete this tutorial, you must be familiar with basic mechanical and electrical
science. You should also be familiar with the Laplace transform and a tutorial on this may be
found in the maths section. You can also find tutorials on fluid power on this site.
Tutorial 2 in this series gives a detailed account of electric motor models and you may wish
to study this first.

D.J.DUNN

1.

INTRODUCTION

Engineering systems is a very broad area ranging from control of a power station to control of a motors
speed. The student needs to have a broad base knowledge of engineering science in order to understand the
various elements and see how many of them are mathematically the same (analogues of each other).
Different kinds of engineering systems often conform to similar laws and there are clear analogies between
electrical, mechanical, thermal and fluid systems. The basic laws which we use most often concern
Resistance R, Capacitance C, Inductance L and conservation laws. You do not need to study all these in
detail and the appropriate law will be explained as required.
Here is a table showing the main analogue components. It is useful to note that capacitance is a zero order
differential equation, resistance is a first order differential equation and Inductance/inertia/inertance is a
second order differential equation.
MECHANICAL

ELECTRICAL

THERMAL

FLUID

Spring

Electrical Capacitor

Thermal capacitor

Fluid Capacitor

x = C F = (1/k) F

Q=CV

Q = C T

M = C x p

Damper

Ohm's Law

## Heat Transfer Laws

Force = kd x velocity
F = kd dx/dt
Torque = kd x Ang.vel

V=RI
V = R dQ/dt

T = R
T = R dQ/dt

## Fluid friction Laws do

not conform to this
pattern.

No equivalence

Fluid inertance

## Newton's 2nd Law of Law of Inductors

motion
Force = Mass x acceleration
V = L d2q/dt2
F = M d2x/dt2
D'Alembert's Principles

Kirchoff's Laws

Force = 0
Moment = 0

current = 0

p = Ld2v/dt2
Law of Conservation Law of Conservation
of Mass
of Energy
Energy = constant

Lets look at the similarity of the various quantities used in these systems.

D.J.DUNN

Mass = constant

2.

SIMILARITY OF ELEMENTS

CAPACITANCE
The symbol C will be used for electrical, thermal and fluid capacitance. Mechanical capacitance is equivalent
to 1/k for mechanical systems where k is the spring stiffness.
RESISTANCE
The symbol R will be used for electrical and thermal resistance. Mechanical/hydraulic resistance is called the
damping coefficient and has various symbols.
INDUCTANCE / INERTIA / INERTANCE
The symbol L will be used for electrical inductance and fluid inertance. In mechanical systems, mass M is
the equivalent property for linear motion and moment of inertia I for angular motion.
OTHER EQUIVALENT PROPERTIES
Q is the symbol for electric charge and quantity of heat. This is equivalent to displacement in mechanical
systems, these being distance (usually x) or angle (usually ).
V is the symbol for electric voltage (potential difference or e.m.f) and is equivalent to temperature T for
thermal systems, Force F for mechanical systems and pressure p for fluid systems.
v or u is the symbol for velocity in mechanical systems and this is equivalent to electric current (I or i) and
heat flow rate .
3

## LAPLACE TRANSFORM and TRANSFER FUNCTIONS

Laplace is covered in detail in later tutorials and in the maths section. The purpose of this transform is to
allow differential equations to be converted into a normal algebraic equation in which the quantity s is just a
normal algebraic quantity. In this tutorial we should simply regard it as a shorthand method of writing
differential coefficients such that:
d
d 2
dn
2
n
becomes s
becomes
s

n becomes s
2
dt
dt
dt
4

TRANSFER FUNCTIONS

The models of systems are often written in the form of a ratio of Output/Input. If the models are turned into a
function of s it is called a transfer function and this is usually denoted as G(s).

G(s) =

Output
The output and input are functions of s.
Input

D.J.DUNN

5.

## BASIC MODELS OF MECHANICAL SYSTEMS

5.1

GENERAL PROCEDURE

## The general procedure for mechanical systems is as follows.

i.
ii.
iii.
iv.
v.
vi.

Adopt a suitable co-ordinate system with an appropriate sign convention. For linear motion, up is
positive and left to right is positive. For rotation anticlockwise is positive and clockwise is negative.
These may be ignored when convenient.
Identify any disturbing forces acting on the system (inputs to the system).
Identify displacements and/or velocities (outputs from the system).
Draw a free body diagram for each mass showing all the forces and moments acting on it.
Apply Newton's 2nd Law to each free body diagram (F = Mass x Acceleration).
Rearrange the equation(s) into a suitable form for solution by a convenient method.

## Note that unless otherwise specified, ignore gravitational effects.

Lets now examine mechanical elements in detail.

5.2

## LINEAR MECHANICAL SYSTEMS.

5.2.1 SPRING
The basic law of a mechanical spring is Force change in length. The diagram shows the model with
mechanical symbols and as a block diagram.

Figure 1
The relationship has no derivatives in it may be written as a function of t or s with no transform involved.
As a function of time we write F(t) = k x(t) where k is the spring stiffness.
As a function of s we write F(s) = kx(s)
x
(s) = 1/k = C
This can be arranged as a transfer function such that
F
C is the reciprocal of stiffness and it is called mechanical capacitance. The use of k is usually preferred in
mechanics but C is used in systems as it is directly analogous to electrical capacitance.

## 5.2.2 DAMPER or DASHPOT

Figure 2
D.J.DUNN

A damper may be idealised as a loosely fitting piston moving in a viscous fluid such that the force is directly
proportional to velocity. F v. Velocity v is the first derivative of distance so F dx/dt
dx
The basic law of a dashpot is:
kd is the damping coefficient.
F(t) = k d
dt
Changed into Laplace form.
F = kd s x
1
x
s
(s) =
Rearranged into a transfer function
kd
F

kd is the damping coefficient with units of Force/Velocity or N s/m. The diagram shows the model with
mechanical symbols and the control block.
5.2.3 MASS

When a mass is accelerated, the inertia has to be overcome and the inertia force is given by Newtons Second
Law of Motion Force = Mass x Acceleration. Acceleration is the second derivative of x with time.
d2x
Basic Law
F(t) = M 2
dt
Changed into Laplace form.
F = Ms2 x
1
x
(s) =
Rearranged into a transfer function
F
Ms 2

Figure 3
5.2.4 MASS - SPRING SYSTEM

For this spring - mass system, motion only occurs in one direction so the system has a single degree of
freedom. It is normal for the direction of motion to be expressed as the x direction regardless of the actual
direction. The free body diagram is as shown. The input is a disturbing force F which is a function of time
F(t). This could, for example, be a sinusoidal force. The output is a motion x which is a function of time x(t).
Let x be positive upwards.

Figure 4

D.J.DUNN

The input force is opposed by the spring force and the inertia force (which always opposes changes in the
motion as stated in Newtons third law of motion).
Spring force = k x
Inertia force = M d2x/dt2
D'Alembert's Principle is that all the forces and moments on the body must add up to zero. In this case it
means
F(t) - kx(t) - M d2x/dt2(t) = 0
or
F(t) = M d2x/dt2(t) + kx (t)
Changing to a function of s we have

## F(s) = Ms2 x + kx = x [Ms2 + k]

F(1/M )
Ms + k
s + k/M
x(s)
1/M
= 2
This may be shown as a transfer function. G(s) =
F(s)
s + k/M
The block diagram for use in systems is as shown.
x(s) =

Figure 5
5.2.5 SPRING DAMPER

## Force balance as a function of time.

Force balance as a function of s
Rearrange into a transfer function.

F(t) = k x + kd dx/dt
F(s) = k x + kd s x
1/k
x
(s) =
k d /k s + 1
F

The units of kd/k are seconds and this is the time constant for the system T = kd/k
1/k This is the standard first order equation which we shall study many times in these tutorials.
x
.
(s) =
Ts + 1
F

Figure 6

D.J.DUNN

## 5.2.6 MASS -SPRING - DAMPER SYSTEM

The input is the force F and the output is the movement x, both being functions of time.
Spring force
Fs = kx
Damping force
Fd = kd dx/dt
Inertia force
Fi = Md2x/dt2
The three forces oppose motion so if the total force on the system is zero then
F = Fi + Fd + Fs
2
dx
x
1/k
d x
+ kx
F(s) = Ms 2 x + k d sx + kx
G(s) =
(s) = 2
F(t) = M 2 + k d
F
dt
s (M/k) + s(k d /k) + 1
dt
If we examine the units of (M/k)1/2 we find it is seconds and this is the second order time constant also with
the symbol T. The transfer function may be written as
x
1/k
G(s) =
(s) = 2 2
F
T s + 2 Ts + 1
is the damping ratio defined as = kd/Cc and Cc is the critical damping ratio defined as (4Mk). The
k
2k d M/k
2k d M
k
=
= d and so the forgoing is correct.
term 2T is hence 2 d T =
Cc
k
4Mk
2 M k k

Figure 7
This is the standard 2nd order transfer function which will be analysed in detail later.

## WORKED EXAMPLE No.1

A mass spring system has the following parameters.
Stiffness K = 800 N/m Mass M = 3 kg Damping Coefficient kd = 20 Ns/m
i.

Calculate the time constant, critical damping coefficient and the damping ratio.

ii.

Derive the equation for the force required when the piston is accelerating.

iii.

## Use the equation to evaluate the static deflection when F = 12 N.

iv.

Use the equation to evaluate the force needed to make the mass accelerate at 4 m/s2 at the
moment when the velocity is 0.5 m/s.

D.J.DUNN

SOLUTION
i.

## T = (M/k) = (3/800) = 0.0612 seconds

cc = 4MK = (4 x 3 x 800) = 97.97 Ns/m
= kd/ cc = 20/97.97 = 0.204

ii.

## For a constant acceleration s2x = a (acceleration) and sx = v (velocity)

1/k
x
(s) = 2 2
F = kx T 2 s 2 + 2 Ts + 1
F
T s + 2 Ts + 1
F = 800x 0.0612 2 s 2 + 2 x 0.204 x 0.0612s + 1

## F = x 3 s 2 + 20 s + 800 = 0.00374s 2 x + 20 sx + 800 x

F = 3 a + 20 v + 800 x
iii.

For a constant force and a static position there will be neither velocity nor acceleration so the s
and s2 terms are zero.
F
12
= 800 x =
= 0.015 m or 15mm
x
800

iv.

## For velocity = 0.5 m/s and a = 4 m/s2

F = 3 a + 20 v + 800 x = 12 + 10 + 800x = 22 + 800 x
The deflection x would need to be evaluated from other methods x = v2/2a = 0.031 m
F = 46.8 N

1.

## A mass spring system has the following parameters.

Stiffness K = 1200 N/m Mass M = 15 kg

## Damping Coefficient kd = 120 Ns/m

i.

Calculate the time constant, critical damping coefficient and the damping ratio.
(0.112 s, 268.3 Ns/m and 0.447)

ii.

If a constant force of 22 N is applied, what will be the static position of the mass?
(18 mm)

iii.

Calculate the force needed to make the mass move with a constant acceleration of 12 m/s2 at the
point where the velocity is 1.2 m/s.
(396 N)

D.J.DUNN

## The following is the rotary equivalent of the previous work.

5.3.1 TORSION BAR

This is the equivalent of a mass and spring. A metal rod clamped at one end and twisted at the other end
produces a torque opposing the twisting directly proportional to the angle of twist. The ratio T/ is the
torsion stiffness of the torsion spring and is denoted with a k.
T is torque ( N m)
is the angle of twist (radian)
k is the torsional stiffness ( N m/rad )
Balancing the torques we have
Change to Laplace form.
Write as a transfer function.

T(t) = k (t)
T(s) = k (s)

1
(s) =
k
T

Figure 8
5.3.2 TORSION DAMPER

A torsion damper may be idealised as vanes rotating in a viscous fluid so that the torque required to rotate it
is directly proportional to the angular velocity. kd is the torsion damping coefficient in N m s/radian
T(t) = k d

d
dt

T(s) = k d s

G(s) =

1
(s) =
T
k ds

Figure 9

D.J.DUNN

## 5.3.3 MOMENT OF INERTIA

Rotating masses oppose changes to the motion and Newton's 2nd law for rotating masses is T = I d2/dt2
I is the moment of inertia in kg m2.
1
d 2

T(t) = I 2 T(s) = I s 2 G(s) = (s) = 2 . Note many text books also use J for moment of inertia.
T
Is
dt

Figure 10

## SELF ASSESSMENT EXERCISE No.2

Derive the transfer function for a mass on a torsion bar fitted with a damper and show it is another
example of the second order transfer function. T is torque and J is moment of inertia.
G(s) =

D.J.DUNN

1/k
1/k
(s) =
= 2 2
2
T
(J/k)s + (Jk d /k)s + 1
T s + 2 Ts + 1

10

## 5.3.4 GEARED SYSTEMS

When a mass is rotated through a gear system, the affect of the inertia is dramatically altered. Consider a
motor coupled to a load through a speed changing device such as a gear box. There is damping (viscous
friction) on the two bearings.

Figure 11
m is the motor rotation and o the output rotation. The gear ratio is Gr = o/m
Since this is a fixed number and is not a function of time, the speed and acceleration are also in the same
ratio.
dm/dt = m
do/dt = o
Gr = o/m is the angular velocity
d2m/dt2= m

d2o/dt2 = o

## Gr = o/m is the angular acceleration.

The power transmitted by a shaft is given by Power = T. If there is no power lost, the output and input
power must be equal so it follows that
m Tm = o To hence
Tm = o To /m = GrTo
(In reality friction significantly affects the torque)
Consider the inertia torque due the inertia on the output shaft Io.
Tm = To x Gr = Io m x Gr2
To = Ioo = Io m x Gr
Now consider the damping torque on the output shaft.
To = kdo o = kdo mGr

Tm = To x Gr = kdo m Gr2

Now consider that there is an inertia and damping torque on the motor shaft and on the output shaft. The total
torque produced on the motor shaft is
Tm = Imm + kdm m + Gr To
To = Ioo + kdo o
Tm = Imm + kdm m + Gr { Ioo + kdo o}
Tm = Imm + kdm m + Gr2 Io m + Gr2 kdo m
Tm = m (Im +Gr2 Io) + m (kdm +Gr2 kdo)
(Im +Gr2 Io) is the effective moment of inertia Ie and (kdm +Gr2 kdo) is the effective damping coefficient kde.
The equation may be written as

Tm = m (Ie) + m (kde)
Tm = d2/dt2 (Ie) + d/dt (kde)

D.J.DUNN

11

## Tm(s) = s2(Ie) + s (kde) = s{sIe + kde}

The output is the motor angle and the input is the motor torque so the geared system may be presented as a
transfer function thus. (s)/Tm (s) = (1/Ie)/s{s + Kde/Ie}

Figure 12

## WORKED EXAMPLE No.2

A DC Servo motor has a moment of inertia of 0.5 kg m2. It is coupled to an aerial rotator through a gear
reduction ratio of 10. The driven mass has a moment of inertia of 1.2 kg m2. The damping on the motor
is 0.1 N m s/rad and on the rotator bearings it is 0.05 N m s/rad.
Write down the transfer function /Tm in the simplest form. Calculate the torque required from the motor
to
i. Turn the aerial at a constant rate of 0.02 rad/s.
ii. Accelerate the rotator at 0.005 rad/s2 at the start when = 0

SOLUTION
i.

## Ie = (Im +Gr2 Io) = (0.5 + 102 x 1.2) = 120.5 kg m2.

Kde = (kdm +Gr2 kdo) = (0.1 + 102 x 0.05) = 5.1 N m s/rad.
(s)/Tm (s) = (1/Ie)/s{s + Kde/Ie}
1/I e

(s) =
s s + k de /I e
Tm

Tm = I e + K de
Tm = 120.5 + 5.1
If the rotator is moving at constant speed (acceleration) is zero. Hence:
Tm = 5.1 = 5.1 x 0.02 = 0.102 Nm
ii. When accelerating at 0.005 rad/s2 the motor acceleration is 10 times larger at 0.05 rad/s2.
Tm = 120.5 + 5.1 = 60.25 Nm when = 0

D.J.DUNN

12

## SELF ASSESSMENT EXERCISE No.3

A DC Servo motor has a moment of inertia of 12 kg m2. It is coupled to an aerial rotator through a gear
reduction ratio of 4. The driven mass has a moment of inertia of 15 kg m2. The damping on the motor is
0.2 N m s/rad and on the rotator bearings it is 0.4 N m s/rad.
Calculate the torque required from the motor to
i. Turn the aerial at a constant rate of 0.5 rad/s. (3.3 N)
ii. Accelerate the rotator at 0.02 rad/s2 at the start when = 0 (20.16 Nm)

D.J.DUNN

13

6.

6.1

## HEATING and COOLING

Consider a mass M kg at temperature 1. The mass is placed in a hot environment at temperature 2 and heat
Q is transferred into the mass causing its temperature to rise. The system could be for example, a resistance
thermometer, and we want to know how long it takes for the sensor to warm up to the same temperature as
the liquid.

Figure 13
The laws of heat transfer tell us that the temperature rise is directly proportional to the heat added so:
dQ = Mc d1 = C d1
c is the specific heat capacity. C = Mc is the thermal capacitance in Joules/Kelvin.
d
dQ
= = C 1
Divide both sides by dt and:
dt
dt
The rate of heat transfer into the mass is = C d1/dt and the rate is governed by the thermal resistance
between the liquid and the mass. This obeys a law similar to Ohms Law so that:

= (2 - 1)/R
Equating for we have

## R is the thermal resistance in Kelvin per Watt.

d1

d1

dt

=
1

1 2

d1

R
2

dt

1 2
RC

+
=
dt
RC
RC
In all systems, the product of the resistance and capacitance is the time constant T so we have
d1

+ 1 = 2
dt
T
T
Changing from a function of time into a function of s we have

1
1
s1 + 1 = 2
1 (Ts + 1) = 2
(s) =
T
T
2
(Ts + 1)

Figure 14
Note that this transfer function is the same standard first order equations derived for the spring - damper
system and thermal capacitance C is equivalent to 1/k and resistance R is equivalent to kd.
D.J.DUNN

14

6.2

## INDUSTRIAL HEATING SYSTEM

The diagram shows a schematic of an industrial process for controlling the temperature of a tank of liquid.
The pneumatic controller will not be explained here but it has an input temperature set by adjustment of the
control. The temperature of the liquid is measured with a suitable device and turned into a standard signal in
the range 0.2 1 bar. This is connected to the controller and pressure sensing devices produce another air
signal (0.2 1 bar) depending on the error. This is sent to a valve that is opened pneumatically working on
the standard range. The overall result is that if the liquid is too cool, steam is allowed through to heat the
liquid.
If the valve opened instead of closing and a cooling fluid was used instead of steam, the system control is by
cooling.
The control equipment could just as likely be all electronic. Pneumatics are used in dangerous environments
such as heating up oil tanks.

Figure 15
The model for the above system will not be derived here but it will be more complicated than simply
1
1
(s) =
because the controller has the facility to do more than proportional control. (Three term
2
(Ts + 1)

## A simple thermal heating system has a transfer function

o
1
(s) =
i
(Ts + 1)

The temperature of the system at any time is o and this is at 20oC when the set temperature i is
changed from 20 oC to 100 oC. The time constant T is 4 seconds. Deduce the formulae for how the
system temperature changes with time and sketch the graph.

D.J.DUNN

15

SOLUTION

o
1
=
i
Ts + 1

i = Ts o + o

i is a constant (100oC) at all values of time after t = 0 (the start of the change).
d
i o = Ts o = T o
dt
Let i o = x
Differentiate and - d o = dx
Rearrange and

## The equation becomes x = T

dt
dx
=
T
x

d o
dx
= T
dt
dt

t
= ln(x) + A
T
t
Substitute for x

= ln( i - o ) + A
T
When t = 0, o = starting temperature = 1 Hence
t
i - 1 = change in temperature

= 0 = ln( i - 1 ) + A A = - ln( i - 1 )
T
( - o )
t
Substitute for A and
= ln( i - o ) ln( ) = ln i
T

t
(

o
Take anti logs and e T = i

t
T

= ( i - o )

o = i e

t
T

## Put in the values T = 4

o = 100 80e-t/4

= 100 20 = 80 i = 100
Evaluating and plotting produces the result below. It is an exponential growth.

Figure 16
D.J.DUNN

16

7.

## HYDRAULIC SYSTEM MODELS

The basic theory for hydraulic and pneumatic components may be found in the tutorials on fluid power.
7.1

HYDRAULIC MOTOR

## The following is the derivation of a model for use in control

theory. The formula relating flow rate Q and speed of rotation
d
is Q = k q = k q
dt
kq is a constant known as the nominal displacement with units
of m3 per radian. is the angle of rotation in radian. Written as
a function of s this becomes Q = kq s
If we take the flow rate as the input and the angle of rotation as

1
the output the transfer function is: G(s) =
=
Q
k qs
The formula that relates system pressure p to the output torque
T is T = kq p

Figure 17

## If pressure is the input and torque the output then

T
G(s) =
= k q This is a further definition of the constant kq.
p

## WORKED EXAMPLE No.4

A hydraulic motor has a nominal displacement of 8 cm3/radian. Calculate the torque produced at a
pressure of 90 bar.
SOLUTION

## T = p kq = 90 x 105 (N/m2) x 8 x 10-6 (m3/rad) = 72 Nm

D.J.DUNN

17

7.2

HYDRAULIC CYLINDER

Figure 18
The flow rate and movement are related by the law Q = A dx/dt. Expressed as a transfer function with x
x
1
=
being the output and Q the input we have: G(s) =
Q
As
Force and pressure are related by the law F = pA. The transfer function with p as the input and F as the
F
= A
output is: G(s) =
p
7.3

## MODEL FOR A FLOW METERING VALVE AND ACTUATOR

Figure 19
The input to the system is the movement of the valve xi. This allows a flow of oil into the cylinder of Q m3/s
which makes the cylinder move a distance xo.

D.J.DUNN

18

Making a big assumption that for a constant supply pressure the flow rate is directly proportional to the valve
position we may say Q = kv xi
kv is the valve constant and examining its units we find they are m2/s
The area of the piston is A m2.
The velocity of the actuator is v = dxo/dt and this is related to the flow and the piston area by the law of
dx
Q = k vxi = A o
continuity such that
dt
Changing to a function of s this becomes kv xi = Asxo
x
1
Expressed as a transfer function we have G(s) = o (s) =
xi
(A/k v )s
The units of A/kv are seconds and we deduce this is yet another time constant T.
x
1
G(s) = o (s) =
xi
Ts
Note that this is not quite the standard first order equation 1/{Ts+1} and the difference is that the output will
keep changing for a given input, unlike the previous examples where a limit is imposed on the output.
If the actuator is a motor instead of a cylinder the equation is similar but the output is angle instead of linear
motion.

## WORKED EXAMPLE No.5

A hydraulic cylinder has bore of 90 mm and is controlled with a valve with a constant kv = 0.2 m2/s
Calculate the time constant T. Given that xi and xo are zero when t = 0, calculate the velocity of the
piston and the output position after 0.1 seconds when the input is changed suddenly to 5 mm.
SOLUTION

## A = D2/4 = 6.362 x 10-3 m2

T = A/kv = 6.362 x 10-3 / 0.2 = 0.032 seconds
G(s) =

xo
1
(s) =
xi
Ts

Tsx o = x i
T

dx o
= xi
dt

x
dx o
0.005 m
= velocity = i =
= 0.156 m/s
dt
T
0.032 s

## Velocity = distance /time

velocity is constant.

D.J.DUNN

19

## SELF ASSESSMENT EXERCISE No.4

1.

A hydraulic motor has a nominal displacement of 5 cm3/radian. Calculate the torque produced at a
pressure of 120 bar.
(60 N m)

2.

A hydraulic cylinder has bore of 50 mm and is controlled with a valve with a constant kv = 0.05 m2/s
Calculate the time constant T. Given that xi and xo are zero when t = 0, calculate the velocity of the
piston and the output position after 0.2 seconds when the input is changed suddenly to 4 mm.
(0.039 s, 0.102 m/s and 20 mm)

D.J.DUNN

20

7.4

Consider the same system but this time let the actuator move a mass M kg and have to overcome a damping
force. Further suppose that the valve now meters the pressure and not the flow rate such that the pressure
applied to the cylinder is p = kvxi. Consider the free body diagram of the actuator.

Figure 20
The applied force is due to pressure Fd and this is determined by the pressure acting on the area A such that:
Fp = pA.
The applied force is opposed by the inertia force Fi and the damping force Fd.
Fi = M d2xo/dt2 and Fd = kd dxo/dt.

d2xo

dx o
dt
dt
2
dx o
d xo
Substituting p = kvxi. we have
+ kd
k vxiA = M
2
dt
dt
2
In Laplace form we have
k v x i A = Ms x o + k d sx o
x
1
Rearranging it into a transfer function. G(s) = o (s) =
2
xi
(M/Ak v )s + (k d /Ak v )s
If we examine the units we find M/Akv = T2 where T is a time constant.
The critical damping coefficient is Cc = (4 M A kv) and the damping ratio is kd/Cc
x
1
G(s) = o (s) = 2 2
The transfer function becomes:
xi
T s + 2 Ts
Balancing forces gives

pA = M

+ kd

Figure 21
Note the similarity with the standard 2nd order equation 1/{T2s2 + 2Ts + 1}. The difference is again due to
there being no limitation on the output. If the actuator is a motor instead of a cylinder, the transfer function is
similar but the output is angle and angular quantities are used instead of linear quantities.
D.J.DUNN

21

## WORKED EXAMPLE No.6

A hydraulic cylinder has bore of 90 mm and moves a mass of 80 kg. It is controlled with a valve with a
constant kv = 20000 Pa/m. The damping coefficient is 180 Ns/m.
Calculate the time constant T, Cc and .
Given that xi and xo are zero when t = 0, calculate the initial acceleration of the mass when the input is
changed suddenly to 5 mm.
Calculate the acceleration when the velocity reaches 2 mm/s.
Calculate the velocity when the acceleration is zero.
SOLUTION

## A = D2/4 = 6.362 x 10-3 m2

T = (M/Akv)= {80/(20000 x 6.362 x 10-3) = 0.793 seconds
Cc = (4 M A kv) = 201.78 Ns/m
= kd/Cc = 0.892
x
1
G(s) = o (s) = 2 2
or in terms of time
xi
T s + 2 Ts
xi = (T2 x acceleration) + (2T x velocity)

## The initial velocity is zero.

0.005 = 0.7932 a + 0
a = 7.952 x 10-3 m/s2
When v = 0.002
0.005 = 0.7932 a + (2 x 0.892 x 0.793 x 0.002)
a = {0.005 - (2 x 0.892 x 0.793 x 0.002)}/0.7932 = 3.452 x 10-3 m/s2
The system initially accelerates and will eventually settle down to a constant velocity with no
acceleration. Put a = 0.
0.005 = 0 + (2 x 0.892 x 0.793) x velocity
velocity = 0.00353 m/s or 3.53 mm/s.

## SELF ASSESSMENT EXERCISE No.5

A hydraulic cylinder has bore of 50 mm and moves a mass of 10 kg. It is controlled with a valve with a
constant kv = 80 Pa/m. The damping coefficient is 2 Ns/m.
Calculate the time constant T, Cc and . (7.98 s, 2.5.7 Ns/m and 0.798)
Given that xi and xo are zero when t = 0, calculate the initial acceleration of the mass when the input is
changed suddenly to 10 mm. (0.157 mm/s2)
Calculate the acceleration when the velocity reaches 0.1 mm/s. (0.137 mm/s2)
Calculate the velocity when the acceleration is zero. (0.785 mm/s)

D.J.DUNN

22

8.

8.1

RESISTANCE

## Applying Ohm's Law we have V = I R

V/I = R
This may be a function of time or of s.
The equation may be expressed in terms of charge Q.
V
G(s) =
(s) = sR
Since I = dQ/dt I(s) = sQ hence
Q
This is similar to the model for the damper.
Figure 22
8.2

CAPACITANCE
The law of a capacitor is Q = C V V/Q = 1/C
This is similar to the model for spring.
Differentiating with respect to time we have dQ/dt = C dV/dt
dQ/dt is current I so the equation may be expressed as I = C dV/dt
As a function of s this becomes I (s) = C sV
1
V
The transfer function is G(s) =
(s) =
sC
I
Figure 23

8.3

INDUCTANCE

dI
d 2Q
= L 2
dt
dt
This is similar to the model for a mass and can be either a first or 2nd order equation as
required.
Expressed as a function of s we have V (s) = L sI or Ls2Q
V
G(s) =
(s) = sL or s 2 Q
I

## Faraday's Law gives us V = L

Figure 24
8.4

POTENTIOMETER

Figure 25
f the supply voltage is constant and the current is negligible, the output voltage V is directly proportional to
the position x or angle so a simple transfer function is obtained.
V
G(s) = o (s) = constant = k p (linear)
x
V
G(s) = o (s) = constant = k p (angular)

D.J.DUNN

23

8.5

R -C SERIES CIRCUIT
The input voltage Vi is the sum of the voltage over the
resistor and the capacitor so
Vi = I R + I /Cs
Vi = I (R + 1/Cs)
The output is the voltage over the capacitor so
Vo = I/Cs
Figure 26

## The transfer function is then

G(s) =

Vo
I/Cs
1
(s) =
=
Vi
I(R + 1/Cs) RCs + 1

The units of RC are seconds and this is another electrical time constant T. The transfer function may be
written as
V
1
G(s) = o (s) =
Vi
Ts + 1
This is the standard first order equation and is the same as both the spring and damper and the thermal
example.

8.6

L - C - R in SERIES
This is 3 sub-systems in series. In this case we will take the
output as the voltage on the capacitor and the input as the voltage
across the series circuit.
The input voltage is the sum of all three voltages and is found by
Vi = I R + I sL + I/sC
Figure 27

## The output voltage is Vo = I/sC The transfer function is then

V
I/Cs
1
=
G(s) = o (s) =
Vi
I(R + sL + 1/Cs) RCs + CLs 2 + 1
Vo
1
(s) = 2
Vi
s CL + sRC + 1
If we examine the units of CL we find it is seconds2 and we have yet another time constant defined as
T2 = CL and the equation may be rewritten as:
V
1
G(s) = o (s) = 2 2
Vi
T s + 2 Ts + 1
G(s) =

## The damping ratio is defined as =

and

L
is called the critical damping value.
C

L
C
Note this is the standard 2nd order equation identical to the mass-spring-damper system.
4

D.J.DUNN

24

## WORKED EXAMPLE No.7

A Capacitance of 200 F is connected in series with a resistor of 20 k as shown in figure 26. The

1
transfer function is o (s) =
i
(Ts + 1)
The voltage across the resistor is suddenly changed from 3V to 10V.
Calculate the time constant T and derive formulae for how the voltage across the capacitor varies with
time. Sketch the graph.
SOLUTION

## T = RC = 20 x 103 x 200 x 10-6 = 4 seconds

The derivation is identical to that in example 3 simply changing to V we get the result.
Vo = Vi V e
Put in the values T = 4
Vo = 10 7e

-t/4

t
T

V = 10 3 = 8

Vi = 10

Figure 28

## SELF ASSESSMENT EXERCISE No.6

1.

Calculate the time constant for an RC circuit with a resistance of 220 and capacitance of 470 nF.
(103 s)

2.

Calculate the second order time constant and the damping ratio for a R-L-C circuit with L = 5 H,
C= 60 F and R=6.8 .
(17.3 s and 11.8 )

D.J.DUNN

25

ELECTRIC MOTORS

This is covered in greater depth in tutorial 2. Here are the basic models.
9.1

## FIELD CONTROLLED MOTOR,

The main theory of electric motors is covered in another tutorial. It can be shown that for a d.c. servo motor
with field control
T = kf i f
If the motor drives an inertial load and has damping the dynamic equation becomes
d 2
d
+ kd
2
dt
dt
T = Is 2 + k d s = (Is 2 + k d s) = k f i f
T = I

G(s) =

(s) = 2 f
if
Is + sk d

This models the relationship between the shaft angle and the
control current.
Figure 29

9.2

## ARMATURE CONTROLLED MOTOR

It can be shown that the torque is related to armature voltage
d k
and resistance by the formula T = Va k

dt R a

## The torque must overcome inertia and damping as before so

T = Is 2 + k d s
Equating we get
T = (Is 2 + k d s) = (Va ks )

k 2s
k
k
Va
=
Ra
Ra
Ra

k 2s
k
=
Is 2 + k d s +
Va

Ra
Ra

Figure 30
With rearrangement we find

G(s) =

(k/R a )

(s) = 2
Va
Is + sk d k 2s/R a

This models the relationship between the angle of the shaft and the control voltage.
Worked examples and self assessment exercise for this section may be found in tutorial 2.

D.J.DUNN

26

10. CLOSED LOOP SYSTEMS TRANSFER FUNCTION WITH UNITARY FEED BACK.

Consider a simple system with an input i and output o related by the transfer function G(s). If the system is
to be a controlled system in which we require the output to change and match the value of the input (set
value), we must make the input the error e instead of the set value. The error is obtained by comparing the
output value with the input value with the signal summing device. This produces the result e = i - o and
because o is subtracted, this idea is called NEGATIVE FEED BACK. The block diagram shows that the
signal passes around a closed loop hence the name CLOSED LOOP SYSTEM.

Figure 31

G(s) = o
e
G(s) =
G(s) =

o
i - o
1
i
o

substitute e = i - o
divide the bottom line by o
rearrange
-1

i
o

-1 =

1
G(s)

1
1 + G(s)
+1 =
invert
G(s)
G(s)

o
o
G(s)
=
i
1 + G(s)

or

1
1/G(s) + 1

o
1
=
i
1/G(s) + 1

## G(s) is the transfer function of the open loop system.

Lets revisit the hydraulic open loop transfer functions derived previously.
When the hydraulic valve and actuator is turned into a closed loop system, the two transfer functions
become:
1
1
for the first order version and G(c.l) = 2 2
for the second order version.
G(c.l) =
Ts + 1
T s + 2 Ts + 1
These models are mathematically identical to the transfer function of the mass-spring- damper and the L-C-R
circuits.
Note that for any system with an open loop transfer function G(s) the closed loop transfer function with unit
1
feedback is G cl =
1/G(s) + 1

D.J.DUNN

27

## WORKED EXAMPLE No.8

An open loop system has a transfer function G(s) = 2/(s2 + 2s + 1). Derive the closed loop function when
unit feedback is used.

SOLUTION
G cl =

1
1
2
2
= 2
= 2
= 2
1/G(s) + 1
s + 2s + 1
s + 2s + 1 + 2
s + 2s + 3
+1
2

## SELF ASSESSMENT EXERCISE No.7

1.

An open loop system has a transfer function G(s) = 5/(4s2 + 2s + 2). Derive the closed loop function
when unit feedback is used.
Gcl = 5/(4s2 + 2s + 7)

2.

An open loop system has a transfer function G(s) = 10/(s3 +5s). Derive the closed loop function when
unit feedback is used.
Gcl = 10/(s3 + 5s + 10)

D.J.DUNN

28
{"found_math": false, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 0, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0}