Mechanical Engineering Department Middle East Technical University

ME 304 CONTROL SYSTEMS

EXERCISES

I/ Introduction & Modeling

PROBLEM I-1

A centrifuge used in training pilots is shown in the picture

and is illustrated with a conceptual block diagram below. The (t)
centrifuge consists of rotating arms, carrying a gondola,
rotated at a controlled speed to simulate the high
accelerations to which the pilots are subjected during flight.

Associate the main function of each component with one or

more of the elements of the list {reference input, comparator,
controller, amplifier, actuator, plant, sensor, transducer,
controlled output}, and write on the rectangle just above the
relevant block.

Ans. :

Comparator +
Reference Controller
Actuator Amplifier Plant
Input Output
e(t) ea(t) DC electric
Pot. Amplifier Gearbox Centrifuge
motor
(t) + _ (t)
Tm(t) T(t)
et(t) Tachometer

Transducer

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-2
Automotive cruise control systems are used to hold the vehicle speed steady at a set value
especially in long runs to prevent driver fatigue. After switched “ON” by the driver, this
system takes the control of the gas throttle using an electronically controlled stepper motor,
and keeps the vehicle speed constant. This system requires speed sensors to measure the
controlled vehicle speed.
a) Identify basic elements of the cruise control system. (Reference input, output,
actuator, controller, sensor, plant)
b) List any possible disturbances to the system during operation.
c) Draw a simple block diagram. Clearly identify each block as one of the basic
elements.
d) Classify (also explain the reason) the system as;
i. Open Loop or Closed Loop,
ii. Regulator or Servomechanism.

Ans. : Disturbances

Output
Error Detector Controller Actuator Plant
+

-
Reference
Input Sensor

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-2
Solar collectors need to be perfectly oriented towards sun in order to provide best utilization
of sun light. Therefore, solar collectors must track the sun by continuously changing
orientation during daytime. The orientation of a solar collector can be adjusted by changing
two angles, one of which is the angle between the collector axis and the horizontal (elevation,
angle φ in the figure) and the other represents the rotation of the collector about vertical
(azimuth, angle θ in the figure). These two angles can be varied smoothly using two
servomotors which are driven by a digital control card.

φ
θ

For sun tracking, two alternative methods are possible to determine the orientation (the two
angles, elevation and azimuth) of the sun in the sky;
i. A sun sensor that senses the position of the sun in the sky,
ii. A software that is loaded on the digital control card which predicts the position of
the sun according to date and location.

a) Identify basic elements of both sun tracking systems.

b) Classify (also explain the reason) both systems as;

i.Open Loop or Closed Loop,
ii.Regulator or Servomechanism.
c) List any possible disturbances to sun tracking systems and compare behavior of the
two systems under the effect of disturbances.

d) Draw simple block diagrams for the systems clearly identifying all elements.

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

Ans. :

Sensor Disturbances

Controller Actuators Plant

System 1
Output

Reference
Input

System 2

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-3:
Consider the pendulum illustrated in the figure. In
the equilibrium position shown the massless rigid O
bar connecting the mass to the pivot point O is T

vertical and the spring is undeformed. A torque T g  d

L
at the pivot point is applied to the bar, and the
angular position of the bar with respect to the
k c
equilibrium position is represented by the angle . m
Obtain the relation between the input (torque T) and the output (angle ) for small oscillations
about the origin.

Ans.: mL2  
  cd 2   kd 2  mgL   T

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-4

b k

m O

g

d
d L
L

Consider the swing axle type of independent suspension of a vehicle. The rigid half axle
pivots about the fixed body point O. A strut, which consists of a spring and a damper
coaxially placed, is pivoted on the half shaft at one end and on the body at the other. Neglect
the masses of the strut (spring and damper) and half shaft. Wheel mass is m, the spring
constant is k, and the damper damping coefficient is b. Assume the equilibrium position of
the system to be the horizontal position of the half shaft. For small motions about the
equilibrium position :
a) Identify the elements and write down the elemental equations. Hint : Give a small
rotation to the half shaft first.
b) Write down the structural equation(s) and identify them as continuity or
compatibility equation(s).
c) Obtain the input output relation (input : force F on the mass, output : angle  of the
half shaft)

 bd 2  kd 2
Ans.: c) mL   F
L L

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-5
y2
M

cs ks
y1
m

kt
y0

Although a car is a complicated mechanical system, simple models can be used with
acceptable accuracy in analyzing certain aspects of vehicle dynamics. For example, the well-
known “quarter car model” is a two degree freedom model that is used in ride comfort
studies. In this model one suspension of the car is modeled with its mass that is not supported
by the suspension springs (unsprung mass, m) and its share of the vehicle body mass (sprung
mass, M) supported by suspension springs. ks, cs, and kt represent the suspension stiffness,
suspension damping coefficient, and tire stiffness, respectively, for the suspension
considered. For the given quarter car model:

a) Identify the elements and write down the elemental equations.

b) Write down the structural equation(s) and identify them as continuity or
compatibility equation(s).
c) Obtain the differential equations relating the input (road profile yo) and the vertical
motion (y2 and y1) of the sprung and unsprung masses.

.. . . .. . .
Ans. : c) M y 2  cs y 2  k s y 2  cs y1  k s y1 , m y1  (k s  k t )y1  cs y1  k s y 2  cs y 2  k t y o
Note that this is a 2 dof system.

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-6 : (See Nise, P2.39)

A disk of inertia J and radius r is
(t)
constrained to move about the stationary J
axis A. Viscous damping with a r A
bv1
coefficient bv1 exists between the disk
and the mass M. Similarly, viscous K f(t)
M
damping with a coefficient bv2 exists
bv2
between the mass and the ground. The
motion of the mass is restrained by an x(t)

ideal spring of stiffness K.

Denote the angular position of the disk and the position of the mass by (t) and x(t). The
input to the system is the force f(t) applied to the mass and the outputs are x(t) and (t).

a) Identify the elements and write the elemental equations.

b) Write the structural equations.
c) Obtain the equation(s) of motion for the system.

Hint : This is a two-degree-of-freedom system and thus two equations are needed.

   b v1  b v2  x  Kx  rb v1  f (t) , J

Ans. : c) Mx   r 2 b v1  rb v1x  0

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-7:
A fluid system consists of two
tanks with cross-sectional areas of
Qi
A1 and A2. Valves of resistance g
coefficients R1 and R2 are fitted to
h1
the outlet of each tank. The fluid R1 Qm
flows from the first tank into the
second tank as illustrated in the A1
h2
figure. Fluid flows into the first R2 Qo
tank and out of the second tank at A2
rates of Q1 and Q2. The density of
fluid is denoted by .

Follow the procedure outlined in the first problem and obtain the input output relation (input :
Qi, output : Qo).

A1A 2 R1R 2  d 2 Qo  A R  A R  dQ 
o
Ans.: Qi    1 1 2 2
 dt   Qo
2 g 2  dt 2  g  
 

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-8

qi
g

C p1 p2 C
R R

q1 q2
Ro

qo

The fluid system shown above consists of two tanks each with a capacitance of C and three
valves of resistance coefficients R and Ro, as illustrated. The fluid flow rates from the two
tanks, and the pressures at the bottom of the tanks are denoted by q1 and q2, and p1 and p2. The
density of fluid is denoted by .

a) Identify the elements and write down the elemental equations.

b) Write down the structural equation(s) and identify them as continuity or
compatability equation(s).

c) Derive the input (qi)-output (qo) relation for the system.

dq 0
Ans. : c) C(2R o  R)  q 0  qi
dt

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM I-9:

A hydraulic accumulator is illustrated in the figure. It is y

k
used to damp out pressure pulses by storing fluid during
pressure peaks and releasing fluid during periods of low q p
pressure. qi qo

Determine the governing equation for this system in terms of the variables p (pressure at the
junction, qi (incoming volume flow rate), and qo (outgoing volume flow rate) which represent
small deviations from the equilibrium conditions.

Neglect any friction as well as any pressure drop at the junction. The piston has a mass m and
area A, the spring constant is k, and y is the displacement of the piston from the equilibrium
position.

Ans.: m   q o   k  qi  q o   A 2 p
qi  

PROBLEM I-10:

In order to harden a steel shaft, it is heated to a temperature of Tos and then it is quenched in a
small can containing mb kg of water at an initial temperature of Tob. Thus while the shaft is
cooling down, the bath temperature will be increased as a result of the quenching process.
The convective heat transfer coefficient, and the surface area and mass of the steel shaft are
denoted by h, A, and ms, respectively. The specific heats of the steel shaft and the water are
denoted by cs and cb.

Obtain the differential equations governing the variations of the shaft temperature Ts and bath
temperature Tb as functions of time.

m b cb ms cs d 2Ts dT m c m c d 2Tb dT
Ans.:  (m b cb  ms cs ) s  0 , b b s s  (m b cb  ms cs ) b  0
hA 2 dt hA 2 dt
dt dt

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

II/ Transfer Functions & Block Diagrams

PROBLEM II-1:
Consider a simple dynamic model of a vehicle travelling on a
V
rough road surface. The mass represents the mass of the body y
and the spring and damper represent the suspension springs m
and dampers. The input is the road surface profile velocity,
k c
i.e. z . If the output is the displacement y of the body,
a) obtain the input-output relation, z
b) write the transfer function,
c) draw the block diagram.

  cy  ky  cz  kz
Ans.: a) my

PROBLEM II-2:
The mechanical system V(t) (input)
shown in the figure b J, Ri, Ro
illustrates a rack and pinion
Pinion
pair. The pinion has a
Disk
companion disk as shown. Rack

The mass of the translating rack is denoted by k

m
m, b denotes the viscous damping coefficient
of the rope, and k is the stiffness of the
restraint on the rack. J, Ri and Ro represent
total inertia of the pinion and the disk, pinion x(t)
(output)
and disk radii, respectively.
For this system,
a) obtain the input-output relation,
b) write the transfer function, and
c) draw the block diagram.

J ..  R 2  .  R 
Ans.: a) (m  ) x   b i  x  kx   b i  V
R o2  R o2 
   Ro 

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-3:
The transmission system shown below consists of 2 shafts, 2 gears and a wheel. One of the
shafts is assumed to be rigid while torsional stiffness of the other is denoted by k. Inertias of
the shafts, gears and the wheel are lumped and represented by I1 to I4. R1 and R2 are the radii
of the pinion and the gear. Viscous friction exists in all bearings with viscous friction
coefficient b. If the input of the system is the torque applied on shaft 1 and the output is the
angular position of wheel,
a) obtain the input-output relation,
b) write the transfer function, and
c) draw the block diagram.

Tin R1
k
I1 II12
Wheel
b b

R2 θout
Rigid Shaft
I3 I4
b b

.. . .. .
Ans. : a) I1 1  b 1  k1  Tin  kn3 and  I2 n 2 + (I3 + I 4 )  θ3 + (n 2 + 2)b θ3 + kn 2θ3 = knθ1
 

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-4:
For the system illustrated in the figure, model the turbine as a resistive element (load) with a
resistance coefficient of RL and
a) obtain the input-output relation,
b) write the transfer function,
c) draw the block diagram.
Note that the lake is assumed to be have infinite capacitance (show that the time rate of
change of pressure becomes zero for any flow rate under this assumption !) and thus can be
modeled as a constant pressure source.
Lake
Surge Tank (Cf)

Turbine (RL)
Long pipe (Ip, Rp) Q1 P2 P3
P1 Qo
(input Valve (Rv) (output)

d 2 Qo dQo
Ans.: a) Cf I P (R v +R L ) +  Cf (R v +R L )R P +IP  +(R v +R L +R P )Qo  P1
2 dt
dt

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-5:
In the cooling system shown below, heat is removed at a rate, Qout, from the cooling room.
Although the walls of the cooler are isolated, there is still heat transfer between the cooler and
the surroundings. The thermal resistance of the cooler walls is denoted by R and the ambient
temperature is given as To. The cooler room has a thermal capacitance of Cc and the

temperature inside the cooler is denoted by Tc. Consider a mass m, with a specific heat cm,
placed in the room. Taking the convective heat transfer coefficient between the mass and the
cooler room as h, and the surface area of the mass as A,

a) obtain the differential equation relating Tm (temperature of the mass) to Qout and To,

b) write the transfer functions between Tm and Qout, and Tm and To,
c) draw the block diagram.

Cc, Tc
To

h, A
Qout
R
m, cm, Tm

d 2Tm dTm
Ans.: a) RCc R m Cm  (RCm  RCc  R m Cm )  Tm  Rq out  T0
2 dt
dt

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-6:
Consider the hydro-mechanical system illustrated in the figure below.
a) Identify and list all the elements in this system. Write down all the elemental
equations. Note that you can write two elemental equations for the cylinder&piston
as it is a transformer between hydraulic and mechanical parts of the system.
b) Write down all the structural equations; first the continuity and then the
compatibility equations.
c) Draw the detailed block diagram for the system making use of the Laplace
transforms of the elemental and structural equations.

qi (input)
g

k
kb
Cf p1 R m
qo p2
c
Ap
x y (output)
Ans. : c)

1
Fk(s)
P2(s) Ap k
_ Fm(s)
Qt(s) P1(s) _ Qo(s) X(s)
Qi(s) 1 1 1 1 sY(s) 1 Y(s)
kb +
+ _ Cf s R Aps + _ _ ms s
+ Fp(s)
Qo(s) Y(s) c
Fc(s)

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-7:
A hydraulic system is illustrated in the figure below. Assume that the fluid has a density 
and is incompressible, h2>h3, and the pipe is short. Rp, Rv, A1, and A2 denote the valve and
pipe resistances and tank areas, respectively.
a) Identify the elements and write down the corresponding elemental equations.
b) Obtain the differential equation relating the input (fluid flow rate qi into tank 1) to
the output (height of the fluid h2 in tank 2).

qi

Pipe, Rp g

h1 h2
q1
h3

Rv qo
Tank 1 A2 Tank 2
A1 Valve

 Rp A1A 2  d2 h2   Rp   dh2  ρg 
Ans. : b)   +  A1  1+  + A2  +  h2 = qi
 ρg  dt 2   R v   dt  R v 

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-8:

R1 i1 i2 R2

i3
+
Vi(t) _ C L Vo(t

For the electrical circuit illustrated above, which consists of two resistances (R1, R2), a
capacitance (C), and an inductance (L),
a) obtain the governing differential equations for the system,
I (s) V (s)
b) find the transfer functions 2 and o ,
Vi (s) Vi (s)
c) draw the block diagram considering Vi(t) as input and Vo(t) as output.

di1 i -i dV di i -i d 2i
Ans. : a) R1  1 2  i , 2 R 2  1 2  L 2  0
dt C dt dt C dt 2

PROBLEM II-9:
Complete the block diagram of a system represented by the equations below. Make sure that
correct signs at the summers are also inserted and enter the relevant variable name on every
signal.

di
eb  K b , L   b  T  Ki,   
 Ri  eb  ea , J
dt

Ea (s)
1 1 1 θ (s )
Ls Js s

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-10:
Simplify the block diagram to obtain the overall transfer function. You can choose either
block diagram algebra or manipulate the equations for the blocks.

A Q
G1

R(s) + C(s)
+ E P
G2
_ +
B T
H1
+
_ S
C(s) (1+G1 )G 2 H2
Ans. : =
R(s) 1+G 2 (H 2 -H1 )

PROBLEM II-11:

Simplify the block diagram down to a single block and obtain the overall transfer function for
the system. You can choose either block diagram algebra or manipulate the equations for the
blocks.
A Q
K1 G3
_ + C(s)
R(s) E P S
G1 G2
+ +
_
B T
K2 H2

C(s) G 1  G 2  G3 
Ans. : 
R(s) 1  K1G1  K 2 G1G 2 H 2

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-12:
Reduce the block diagram, given below, to a single block and write down the overall transfer
function for the system.

Q
H1
_
R(s) E P C(s)
G1 G2
+
_
B1 B2

+ W + X
S T V
H3 H2
+ U +

C G1G 2
Ans. : =
R 1+ G1B1 + G1G 2 H1 + G1G 2 H3B2 + G1G 2 H 2 H3

PROBLEM II-13:
C(s)
Obtain the closed loop transfer function, T(s)= , for the block diagram given below.
R(s)

H1

R(s) + + C(s)
+ + +
G1 G2 G3
_ _ _

H2

H3

(H1  G 2 )G1G 3
Ans. :
1  (H1  G 2 )G 3H3  G 2 H 2  (H1  G 2 )G1G 3

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM II-14:
C(s)
Obtain the closed loop transfer function, T(s)= , for the block diagram given below.
R(s)

G8 C(s)

R(s) + _
G1 G5 G6
_

G2 G4 G7

+ +
+ +
G3

C(s) G1G 5
Ans. : =
R(s) 1+G1G 2 +G1G 5G8 +G1G 3G 5 (G 4 +G 6G 7 )

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

III/ Control System Components

PROBLEM III-1:
A linear actuator (solenoid) consists of a mass m, and a parallel connection of a linear spring
of stiffness k1 and a linear viscous damper of coefficient c. The resistance of the solenoid coil
is R and its inductance is L. The electromagnetic coupling constant is K. The actuator excites
one end of a massless rigid bar which is pivoted at point O as illustrated in the figure. The
motion of the other end of the bar is restrained by a linear spring of stiffness k2.
a) Write down the elemental and structural equations for the system. Neglect
gravity.
b) Draw the detailed block diagram showing each element by a block. Do not
combine blocks.
c) Obtain the overall transfer function.

L1 L2
y1
O y2
i R, L
e k2 (output)
m
(input)

c k1

Ans.:
Y2 (s) K

E(s)  L1  3  L1  2  2 L1 L   L L 
 mL  s   (mR  cL)  s  (cR  k1 L  K )  k 2 L 2  s   k1 1 R  k 2 2 R 
 L2   L2   L2 L1   L2 L1 

- ..
+ 1 I(s) FP(s) - m Y1 1 1 Y1 (s) -L 2 Y2 (s)
K
E(s) Ls s L1
- -
ms
- +

R c
k
L2
k2
L1

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM III-2:

A hydraulic servomotor, controlled by Load Power Piston & Cylinder Ideal Input
Damper Lever
a rigid massless lever, drives a load as Ap
y
illustrated in the figure. m
p2 p1 c
a) Write the differential a
z q Pilot valve
equation representing the input—
(output)
output relation for the system taking x x b
(input)
as the input an z as the output. drain Ps drain

b) Obtain the overall transfer function. Note that m represents the total mass of the
load, piston, and piston rod.

You can use the linearized relation q  K1x  K 2 p .

  ab
Ans. : K 2 mz  A 2p  K 2 c z  K 2 c   x  K1A p x
 b 

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM III-3:

A centrifuge used in training pilots is shown in the picture

(t)
and is illustrated with a conceptual block diagram below. The
centrifuge consists of rotating arms, carrying a gondola,
rotated at a controlled speed to simulate the high
accelerations to which the pilots are subjected during flight.

a) Write down the elemental equations for the component

models and obtain the overall transfer function relating the
input (t) (angular position of the potentiometer knob) to the
output (t) (angular speed of the centrifuge arms).

The main specifications of the system components are listed below.

Centrifuge : Moment of inertia J.
Gearbox : Reduction ratio n, negligible moment of inertia.
DC Servomotor : Armature resistance R, Armature inductance L, Torque constant K, Back
emf constant Kb, negligible moment of inertia.

DC Amplifier : Amplification factor Ka.

Potentiometer : Constant Kp.

Tachometer : Constant KT.

b) Draw a detailed block diagram for the system.

Ans. :

 nKK a K p 
a) Ω(s) =   θ(s)
 JLs 2 + Rs + nK  K a K t + nK b  

Ea(s) I(s Tm(s) Tg(s) Ώ(s)

(s) + 1 1
Kp Ka ) K
+ _ n
_ Ls Js
_
R
Ώm(s)
Kb n

Kt
b)

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM III-4:
Consider the system illustrated in the figure. Note that the bar AB is rigid and massless
(ignore the bar between point B and the mass). Consider small displacements only. The
solenoid parameters are the mass m of the plunger, ideal spring constant k, viscous damping
coefficient c and resistance R and the inductance L of the coil. The electromagnetic coupling
constant, coil current, and the voltage applied to solenoid terminals are denoted by Kb, i, and
e. The load is modelled by a mass M and a spring of constant K, together with viscous
friction of coefficient C between the mass and the ground.
e
i x1

Solenoid A
a
x2
m, k, c,
b K
R, L, Kb
M
B C

a) Represent solenoid dynamics with two differential equations; one for the mechanical
and one for the electrical side. Note that a resistive force FL acts on the mechanical
side due to the load represented by the bar AB and the mass-spring-viscous friction
subsystem.

b) How many degrees of freedom does the system have ? Explain.

c) Write down all the relevant elemental and structural equations for the rest of the
system.

d) Draw the detailed block diagram of the system with e as the input and x2 as the
output.

Ans. :

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

k
Kb
c
_ _ s2X1(s) sX1(s) X1(s)
E(s) _
1 I(s) + 1 1 1 a X2(s)
Kb 
+ Ls _ m s s b
_
R a a K
FL(s)  M  C
b b

+ +

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM III-5:

An electro-hydro-mechanical system is shown in the figure. A dc electric motor drives a

positive displacement pump and is controlled by the armature voltage ea. Relevant parameters
are indicated on the figure. The total inertia of the motor and the pump is denoted by JT. The
load is represented by a mass-spring system. Viscous friction is assumed between the mass
and the ground.

The pump flow rate q is assumed to be proportional to to shaft speed, i.e. : q = K f 

The resistance and inertance of the pipes as well as the friction between the piston and the
cylinder are negligible.

+ ea - x
Tm Ap

q k
  
DC Positive m
Electric Disp. p
Motor Pump c
JT

Ra = armature resistance Kf = pump flow constant m = mass

La = armature inductance JT = total inertia (motor + pump) k = spring constant
Kb = back emf constant  = shaft angular position c = coefficient of viscous friction
Kt = torque constant Ap = piston area

a) Write down a complete set of equations for the dc electric motor. Clearly indicate the
terms coupling the electrical and mechanical parts. Draw a detailed block diagram of
the motor (Input : armature voltage, Ea; Output : Motor shaft angular position, m).

b) Draw a block diagram of the pump (Input : shaft speed, ; Output : pump flow rate,
q).

c) Write down a complete set of equations for the hydraulic actuator (cylinder+piston).
Draw a block diagram of the actuator (Input : pump flow rate, q; Output : load
displacement, x). Note that you need to use only one of the equations to draw the
block diagram.

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

d) Now assemble all the previously drawn individual block diagrams (Input : armature
voltage, Ea; Output : load displacement, x). Note that the angular position of the
shaft is of no interest to us and thus need not be included in the block diagram.

e) Note that the effects of the load dynamics (mass+spring+friction) do not appear in
the block diagram. This is because once the flow rate is specified, the velocity of the
mass is dictated. On the other hand, load dynamics determine the pressure at the
outlet of the pump. Thus write a differential equation to relate pump pressure p to
the displacement x of the mass. Draw a detailed block diagram of the load (Input :
load displacement, x; Output : pressure at pump outlet, p). Do not worry if you have
to enter “s” or its powers in the numerator !

Ans. :

+ 1 Ia(s) Tm(s)
Ea(s) + 1 s(s) 1 (s)
Kt
_ La s JT s s
_
Eb(s) Ra
Kb
a)

(s) Q(s)
Kf
b)

Q(s) 1 X(s)
Aps
c)

Ea(s) + 1 Ia(s) Tm(s) 1 s(s) Q(s) 1 X(s)

Kt Kf
+_ _ La s JT s Aps

Eb(s) R

Kb
d)

X(s) sX(s) s2X(s 1 P(s)

s ms
+ Ap
+ +
c

e) k

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM III-6: (See Nise, Problem 2.46)

Consider the electro-mechanical system N2=20

ct=1 Nm/rad/s
illustrated in the figure. A fixed current is
+ ea -
applied to the field circuit of the electric
motor. You can neglect the inertia of the DC R=2 m,
Electric J=1 kgm2
armature rotor and the inductance of the Motor
armature circuit. Ra=1 N1=10
Kb=1 V/rad/s
a) Write the elemental equations.
Kt=1 Nm/A
b) Draw the corresponding detailed m=1 kg x(t)
block diagram of the system. Use symbols
c=1N/m/s
only in the blocks. Do not insert numerical
values yet !
c) Insert the numerical values for the parameters and find the overall transfer function
X(s)/Ea(s) – which should not contain any symbols with the exception of Laplace variable
“s”.

Ans. : b)
1
Ct
R
1
2(s) J
R

Ea(s) 1 1(s)  1 F(s) 1 sX(s) 1
Kt n
 R a a(s) m(s)  RF(s) R  2
s X(s) Ms s X(s)
 
b(s) 4 sm(s)
c
s(s) 1
5 Kb n
c) G(s)  R
9
s2  s
4

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

IV/ Stability
PROBLEM IV-1:
Write down the characteristic equation for the 2 C(s)
R(s)
Kp
system with the block diagram shown in the + _ s1s2
figure.
s4
 
Ans. : s3 + 3s 2 + 2 1+ K p s + 8K p = 0 s

PROBLEM IV-2:
Block diagram of a position control system for a space camera is given. Show that the system
is unstable and determine the number of poles with positive real parts. All constants take
positive values only.
Integrating
Amplifier Generator Motor
Potentiometer Kg
d(s) A Km (s)
Ka Rf
+ _ s s  Tm s  1
Tm s  1

Mechanical linkage

Ans. : There are two roots with positive real parts.

PROBLEM IV-3:
Consider the control system for an inertial load, as represented by the block diagram shown.
a) Determine the value of the open loop gain (K=bKaKmKf) which results in a marginally

stable system, if the time constants for the motor (Tm) and for the load (TL) are specified
to be 0.05 and 0.20 seconds, respectively.
b) What is the undamped frequency of oscillation of the marginally stable system ?
Control Inertial
Amplifier Motor Load
d(s) Km b (s)
Ka
+ _ Tm s  1 s  TL s  1

Kf
Sensor

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

Ans. : a) 0<K<25, b) for K=25 the system poles are -25, -10j, 10j, which yield asymptotic
stability and the undamped frequency of oscillations of the marginally stable system 10 rad/s;
for K=0 the system poles are -20, -5, 0, which again result in asymptotic stability and in this
case the response is not oscillatory.

PROBLEM IV-4:

For the two characteristic equations given, use Routh’s array to determine whether there are
any roots with positive real parts, and if there are, the number of such roots.

a) D(s)= s5+7s4+19s3+25s2+16s+4

b) D(s)= s6+2s5+8s4+12s3+20s2+16s+16
Ans. : a) None, b) None.

PROBLEM IV-5:
For the characteristic equations listed below, use Matlab to define the characteristic equation
by a row vector of the coefficient of powers of s as illustrated. Then print the command
roots(p) and press enter to find the roots. Comment on the stability of the systems.

s3+12s2+46s+52=0 → p=[1, 12, 46, 52] → roots(p)

a) 3s4+10s3+5s2+8s+8=0

b) s6+2s5+8s4+12s3+20s2+16s+16=0

Ans. : a) -2.9729, 0.2533 + 0.9851i, 0.2533 - 0.9851i, -0.8670, b) -0.0000 + 2.0000i, -0.0000
- 2.0000i, -1.0000 + 1.0000i, -1.0000 - 1.0000i, 0.0000 + 1.4142i, 0.0000 - 1.4142i

PROBLEM IV-6:
Consider the block diagram for a control system. Determine the ranges of the control
parameters Kp and Ti permissable for a stable system and plot the permissable range on the
parameter plane.

R(s)  1  1 C(s)
K p 1  
+ _  Ti s   s  3
1
s 1

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

Kp
Ans. : Ti  , Ti  0, K p  0
4(3  K p )

Stable

Unstable

PROBLEM IV-7:

The characteristic polynomial of a control system is given as

D(s)  Ts3  4Ts 2  (K  3)Ts  K

Find and plot the admissable ranges of T and K on the parameter plane such that the system
has a stability margin of at least 1.

0.95

0.9

0.85

0.8

0.75
T

Stable
0.7

0.65

0.6

0.55
Unstable
0.5
2 4 6 8 10 12 14 16 18 20
K

K
Ans. : T  , T 0
2( K  1)

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM IV-8:
Consider the system with the block diagram Controller
as shown.
R(s)
Kp  s  40 C(s)
+ _ s  s 10
a) Determine the value(s) of the controller
1
gain Kp that will make the system
s  20
marginally stable.
b) What will be the frequency of oscillations of the unforced response of the marginally
stable system to an initial condition ?
c) If the value of Kp is selected as 200, determine if the stability margin of the system is less
or more than 1.
Ans. : a) Kp=600 and Kp=0, b) 28.3 rad/s, c) More than 1.

PROBLEM IV-9:
Consider the unity feedback system with the R(s) 5 s  20 C(s)
block diagram as shown. + _ s  s  2 s  3
a) Show that the system is unstable. Find the
number of unstable poles.
Assume that a proportional controller of gain K>0 is added to the system.
b) Check if it is possible to stabilize the
system by a proper choice of the R(s) 5 s  20 C(s)
K
proportional gain K. If the answer is yes,
+ _ s  s  2 s  3
determine the range of values of K for
which the system is stable.
c) For which value(s) of K will the system be marginally stable? Determine the frequency of
oscillations of the unforced response of the marginally stable system.
d) Determine the range of values of K such that the stability margin of the system is at least
0.5.
Ans. : a) Two poles with positive real parts, b) 0<K<2, c) K=2, 2.83 rad/s, d) 0.096<K<0.500.

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

V/ Transient Response

PROBLEM V-1: cs(t)

Estimate the time constant (illustrate on the figure) AK

for a first order system with the step response as

shown in the figure. Time [s]
1 2 3 4 5
0 1 2 3 4 5 6
Ans. : ~2 seconds

PROBLEM V-2:
A thermometer, reading 27oC initially, is T [oC]
tested by dipping it in water at a constant 27
temperature of 15oC. The variation of the
thermometer reading is recorded and plotted
15
as shown in the figure.
The same temperature brought back to steady
state once more such that it reads 27 oC just
0 3.6 Time [s]
above the surface of a lake.
It is then is lowered into the lake at a constant speed of 0.3 m/s. Water temperature in the lake
decreases by 1.8oC for each meter below the surface. Determine the thermometer reading
when the thermometer reaches a depth of 3 meters. What is the reading error for the
thermometer ?
Ans. : 22.46 o C , -0.86 o C

PROBLEM V-3:
Transfer functions for two different first order systems are given below. Which of these two
systems will have a faster response ? Explain.

4 0.1
G1 (s) = G 2 (s) =
3s +1 s + 0.2
Ans. : System 1

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM V-4:
The water level in the tank, shown in the figure, as well as the volume flow rates in and out of
the tank are initially constant at h, Qi , and Qo . The volume flow rate into the tank is
suddenly increased by 2 m3/s. The increase in water level in the tank is plotted versus time.
Determine the effective cross-sectional area of the tank.

H
c s[m]
(t)
K
Qi + Qi g 10

h
h p+p R1

Qo + Qo TimeTime
[s]
A 0
1 2 3 4 5

0 T
12

Hint : Write the input-output relation originally at steady state and then after the step input in
flow rate. Be careful with the units - make sure you write the units for all calculated values !
Ans. : 2.4 m 2

PROBLEM V-5: c(t)

The response of a system, to an impulse of 5
magnitude (strength) 0.1, is measured and
plotted as shown in the figure. Estimate the
response of the same system, 3 seconds after
the application of a step input of magnitude
0.2.
0 2 Time [s]
Ans. : 15.5

PROBLEM V-6:
If the overall transfer function of a third order system is given as
32
G(s) 
 s  12   s2  4s  20.25
calculate an approximate value for the undamped natural frequency.
Hint : Consider the locations of the poles

Ans. : 4.5 rad/s

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM V-7:
Consider the single degree of freedom mass-spring-
y
damper system illustrated in the figure. Assume that a M+m
load of mass m is placed on the mass M. Model the g
k c
presence of the load as an initial condition on the
displacement of the mass. Find an equivalent transfer
function, the unit step response of which will represent the
motion of the mass when the load is suddenly lifted off. M = 250 kg
m = 70 kg
Hint : Examine solved example 5-9 in Ogata. k = 12 kN/m
c = 800 N/m/s
(0.05717s 2 + 0,18294s)
Ans. : G(s) = -
s 2 + 3.2s + 48

PROBLEM V-8:
Consider a standard second order system.
a) Starting with the expression for the general unit step response, obtain an expression
to calculate the rise time.
b) Calculate the delay time for the system specified by an undamped natural
frequency of 2.5 Hz and a damping ratio of 0.5.

Ans. : 0.082 s

PROBLEM V-9:
The unit step response due to black poles only are given. Superimpose the unit step response
when all the poles are considered.
c(t) c(t)
I K K

R t t

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM V-10:
A control system designed to maintain a desired torque on a shaft is represented by the block
diagram shown. The transient response specifications to be provided are selected as follows.
Rise time < 5 s,
Maximum Overshoot < 12 %,
2 % settling time < 15 s
Determine the acceptable ranges for the parameters K1 and K2.
Control
Controller Motor
Td(s) K1 2 Ta(s)
+ _ s s4
Desired Actual
Shaft shaft
Torque K2 torque
Shaft torque sensor

Ans. : 2.16  K1  K 2  6.4

PROBLEM V-11:
Write a script m-file to determine the unit step response of a system represented by the block
diagram below with K1=1.32 and K2=1.14.

a) Plot the response.

b) Calculate the values of all the transient response specifications.
c) Compare the calculated maximum overshoot and the 5% settling time values with
- the values read from the plot, and
- the specifications of the example solved in class.

R(s) K1 1 C (s)
+ + s + 3 s
+
_
K2

Ans. : td=1.2 s, tr=2.6s, Mp=1.07 (tMp=3.58s), ts=4.25s

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM V-12:

Consider the control system represented by the

i(s) 1 o(s)
block diagram shown. 25
+ _ s  s  1
a) Calculate the maximum overshoot and 2
% settling time of the system.
b) A velocity feedback is added to the
i(s) 25 1 o(s)
system, as shown in the block diagram, to
+ _ + _  s  1 s
improve the transient response of the system.
Kt
Determine the value Kt such that the system
will have a damping ratio of 0.7.

Ans. : a) t s  2%   8s , b) Kt=0.24

PROBLEM V-13:
Consider the block diagram of a space
R(s) 0.5 3 C(s)
telescope pointing control system. Kp and + _ _ s s
Kd are positive parameters to be used to
Kd
obtain the desired performance. Kp
a) Determine the relation between the
parameters Kp and Kd such that the
damping ratio is critical.
Plot Kp versus Kd and clearly indicate the regions in which the system is underdamped
and overdamped.
b) Select ranges of values for the parameters Kp and Kd such that for a step input, the
maximum overshoot is 4% and the peak time is less than 2 seconds. In specifying the
required parameter ranges, make sure that you use the proper equality and/or inequality
signs.

Ans. : a) K p  0.042K d2 or K d2  24K p , c) Kp>3.375, Kd>6.45

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM V-14:
Transient response specifications for a unit
R(s) 1 (s)
G c (s)
step input for the system with the block + _ s s  2 
diagram shown in the figure are specified as :
Peak time = 1 s
Percent overshoot = 5 %
a) Show, without any calculations that these two specifications cannot be satisfied
simultaneously by using proportional control Gc(s)=Kp.
b) What type of controller would you use to satisfy both specifications simultaneously
without reducing the stability ? Explain!

Ans. : a) K p  2.1 or K p  18.8 required !.

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM V-15:
Unit step response of a control system is obtained and presented below.
a) Estimate (clearly indicate which values you take from the figure)

- Rise time (0 to 100 %),

- Maximum overshoot (percent),

- Peak time,

- Settling time (5 %).

b) Assume that the system is a second order system without numerator dynamics and
determine the undamped natural frequency and damping ratio.
Step Response
1.5

System: G
Peak amplitude: 1.4
Overshoot (%): 40.2
At time (sec): 4.91 System: G
System: G Settling Time (sec): 16.3
Rise Time (sec): 2.98
1
Amplitude

0.5

0
0 5 10 15 20 25 30 35
Time (sec)

Ans. : a) 3.0 s, 40 %, 4.9 s, 17 s, b)   0.28, n  0.65 rad / s .

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

VI/ Steady State Response and Error

PROBLEM VI-1:
Write down the open loop transfer function C(s)
R(s)
OLTF(s) of the closed loop system with the + _ G1 G2
block diagram as shown.
H
Ans. : G=G1G2H

PROBLEM VI-2:
A unit step input is applied to the system
R(s) K C(s)
represented by the block diagram shown.
+ _ s 2  ps  q
Determine the time domain response, i.e. c(t),
as time goes to infinity, i.e. the steady state a
sb
response, of the system.

Kb
Ans. : css = (check stability first !)
(qb + Ka)

PROBLEM VI-3:
Consider the block diagram of a system designed to control the depth of a submarine. The
actual depth is measured by a pressure transducer and is compared with the desired depth.
The error is amplified and is then used to drive the stern plane actuator.

Determine the value of the

R(s)  s  12 C(s)
Ka
 
amplifier gain such that the
+ _ s s 2  0.1
steady state error for a unit
ramp input will be less than 5
%.

Ans. : K a  2 (check stability first !)

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM VI-4:
What is the steady state error of the unity feedback system with an open loop transfer
function
 1  1
G(s)  K p  
 Ts  1  s 2
to an acceleration input of magnitude R ?
Hint : Is it defined ?

PROBLEM VI-5:
a) If the steady state error for a unit ramp
R(s)
Kp s4 2 C(s)
input of the unity feedback system, + _ s s1s2
represented by the block diagram shown, is
required to be equal to or less than 2%,
determine the minimum value of Kp.
b) Is it possible to use this value for actual operation ? If not, calculate the maximum value of
Kp that can be used and the resulting steady state error.

Ans. : a) K p  12.5 , b) Kp<3, ess  0.08 or 8%

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM VI-6:
Consider the position control of an inertial load with damping. A tachometer is built-in the
electric motor and provides a feedback proportional to the angular velocity of the load.
a) Determine the values of the control parameters Kp and Kt such that the undamped
natural frequency of the system is 1 rad/s, and the system is critically damped.
b) Calculate the values of the control parameters Kp and Kt such that the steady state error
due to a unit ramp input is less than 0.03 when the damping ratio is 0.5 and the
disturbance is zero.
c) What should be the values of Kp and Kt if the steady state error due to a unit step
disturbance is to be less than 0.01. Note that, in this case, the steady state error is defined
as the steady state value of the output itself as the desired input is now zero. Thus, you
should re-draw the block diagram with the disturbance as the input and evaluate the
output as time goes to infinity.
Disturbance
Motor + Load
J = 15 kgm2 Controller +
i(s) 1 1 o(s)
Kp
b = 3 Nm/rad/s + _ + _ +  Js  b  s
Kt
Tachometer

Ans. : a) K t =1.8 , b) Kp > 16667, Kt < 0.0298, c) K p >100 (check stability first !)

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM VI-7:
The block diagram of a unity feedback speed control system is shown in the figure.
a) Determine the type of the system. Explain how you have decided.
b) Obtain the open loop gain K of the system.
c) Decide if it is possible to obtain a steady state error less than 0.04 by a proper choice of
the controller gain Kc, if the input is specified as r(t)=0.9t.
Controller
R(s) 3 1 C (s)
Kc
+ +
_ s + 5  s + 1 0  s
_

Ans. : a) Type 1, b) K=3Kc/56, c) No.

PROBLEM VI-8:
Consider the block diagram of a nose-wheel steering system for medium size aircraft, shown
below.
a) Write down the open loop transfer function and determine the type and open loop gain K
the of the system.
b) Calculate the amplifier gain Ka, to result in the minimum possible steady state error in
response to a unit step input. Further calculate the value of the minimum achievable
steady state error.

Desired Actual
wheel heading Amplifier Clutch Actuator wheel heading
R(s) 1 15 C(s)
Ka
+ _ s5 s 2  3s  2

15K a
Ans. : a) T(s) = , type zero, open loop gain=3Ka/2,
 s +1 s + 2  s + 5
b) Ka<8.4, ess  0.106

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM VI-9:

Consider the block diagram of a unity feedback system shown in the figure.
R(s) 144 C (s)
+ s  s + 3  s + 1 2 
_

a) Determine the system type and open loop gain. Obtain the steady state error for a unit
ramp input.
b) Add a proportional controller of gain Kp>0 to the system. Determine the minimum steady
state error for a ramp input achievable with the controlled system.
c) Which controller type can be used to reduce the steady state error of this system to a
ramp input to zero. Show that you actually get zero steady state error to a ramp input
with your controller.

Ans. : a) Type 1, K=144/36=4, b) Kp<3.75, ess  0.067 .

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

VIII/ Basic Control Actions and Controllers

PROBLEM VIII-1:
Consider the system represented by the block diagram given.

a) If proportional control Gc(s)=Kp is used, determine the value of proportional gain such
that the steady state error of the system due to unit step input is equal to or less than 1%.

b) It is decided to use a proportional+

R(s) 1 C(s)
integral (PI) control. G c (s)
+ _ 4s  1
 1 
G c (s)  K p  1  
 Ti s 
instead of proportional only control. Determine the values of the proportional gain and the
integral time such that the 2% settling time is less than or equal to 4 seconds, and the steady
state error for unit ramp input is equal to or less than 1%.
Ans. : a) Kp ≥ 99, b) Kp ≥ 7, Ti≤ 0.07

PROBLEM VIII-2:
Consider the system with the block diagram as
s+2 C(s)
R(s)
 
shown. G c (s) 2
+ _ s s -3
a) Show that the uncontrolled system (with
no controller : Gc(s)=1) is always unstable.

b) Try proportional control, i.e. Gc(s)=Kp. Can the system be stabilized for some value of

the parameter Kp ?
c) Try now PI control with :
Ki
G c (s)  K p 
s
Can you find some values of the parameters Kp and Ki for which the system is stabilized ?
d) Finally use PD control with :
G c (s)  K p  K d s

Find and sketch the region for the stable operation of the system in the parameter plane.
3  2Kd
Ans. : d) K p 
2
(1  )
Kd

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM VIII-4:
A proportional control system is used for
Reaction
controlling the temperature of an Controller Heater Process
exothermal (unstable) chemical reaction. Tr(s) 2.5 0.8 T(s)
Kp
a) Determine if the system is stable for a + _ 5s 1 3s  1

proportional gain of 100. Transducer

b) To improve the stability of the Reaction
Controller Heater Process
system, derivative control is added to the
Tr(s) 2.5 0.8 T(s)
proportional control. Determine the K p 1  Td s 
+ _ 5s 1 3s  1
value of derivative gain Td for a stability Transducer
margin of at least 1, if the value of Kp is
100.
Ans. : a) Unstable, b) 0.16<Td < 1.08

PROBLEM VIII-3:
Consider the system with the block diagram as 3 C(s)
R(s)
shown. + _ 2s  1
a) Estimate the time required to reach steady
state.

b) What is the steady state error for a unit step input ?

Now insert a PI controller with the transfer function

 1 
G c (s)  K p  1  
 Ti s 

in the control system, and re-draw the block diagram. If the values of the proportional gain
and the integral time are selected as 4 and 0.15 seconds, Use Matlab to :

c) Calculate the 2% settling time.

d) Determine the steady state error for a unit step input.

e) Calculate the values of the proportional gain and the integral time to be used if the steady
state error for a unit ramp input is to be less than 1 % and the 2% settling time is to be
less than 1 second.

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM VIII-4:
Consider the position control of an inertial load with a proportional controller of gain Kp. A
tachometer is built-in the electric motor and provides a feedback proportional to the angular
velocity of the load.
a) Determine the values of the control parameters Kp and Kt such that the undamped
natural frequency of the system is at least 2 rad/s, and the system is critically damped.
b) Calculate the values of the control parameters Kp and Kt such that the steady state error
due to a unit ramp input is less than 5% when the damping ratio is 0.5 and the
disturbance is zero.
c) What should be the values of Kp and Kt if the steady state error due to a unit step
disturbance is to be less than 2.5%

Note : In specifying the required parameter values, make sure that you use the proper equality
and/or inequality signs.

Disturbance
Motor + Load
Controller +
i(s) 1 1 o(s)
Kp
+ _ + _ + 15s  3 s
Kt
Tachometer

3
Ans. : a) Kp ≥ 60 b) K p >6000 , K t  0.05- c) Kp > 40
Kp

PROBLEM VIII-5:
Consider the position control of an inertial load with the block diagram as shown in the
figure. A tachometer is built-in the electric motor. Disregard disturbance and
a) Assume that a proportional controller, Gc(s)=Kp, is used. If the undamped natural
frequency is to be at least 2.4 rad/s and the damping ratio is to be set to 0.72,
determine the (range of) values of the parameters Kp and Kt.
b) Now, consider the disturbance, D(s), as the only input to the system. Redraw the
block diagram accordingly. What are the conditions on the values of Kp and Kt if the
steady state error due to a unit step disturbance is to be less than 0.025 ?

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

 1 
c) Finally, assume that a PI controller, G c (s)  K p  1   , is used. If Kt=0.25,
 Ti s 

check the conditions on the selection of the parameters Kp and Ti for stability. If

necessary plot Kp versus Ti and clearly mark the stable and unstable regions.
Disturbance D(s)
Motor + Load
Amplifier Controller
R(s) + 3.6 1 C(s)
10 G c (s)
+ _ + _ + s s

Kt
Tachometer

10
Ans. : a) K p  0.16 , K t  6 , b) Kp  4 , c) Kp 
0.9  36Ti

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

IX/ Frequency ResponsePROBLEM IX-1:

Calculate the steady state response, x(t), of the system with the overall transfer function

X(s) s
T(s)  
Y(s) 0.2s  1

to a sinusoidal input

y(t)=0.5sin2t.


Ans. : x(t)  0.93sin 2t  68o 
PROBLEM IX-2:
Consider the unity feedback control system with the block diagram shown. Sketch the
approximate Bode plots for the open loop transfer function of the control system.

R(s) s  0.6 C(s)

12 1
+ _ s s 1 s 2  6s  36

Ans. :

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

Approximated Magnitude of Bode Plot

Approximated Phase of Bode Plot

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM IX-3:

The block diagram of a control system for the stabilization of the roll motion of a ship is
given.

a) Identify the basic factors for the open loop transfer function of the system.

b) Identify the corner frequencies separately for the magnitude and phase plots.
Remember that a corner frequency is defined as the intersection of two straight line
approximations !

c) Plot first the straight line approximations for the basic factors and then sketch the
approximate Bode plots for the open loop transfer function of the system by
graphically adding the contributions of the basic factors. Make sure that the plots are
clear and readable.

Desired Fin Ship Roll Actual

Roll Angle Controller Actuator Dynamics Roll Angle
d(s) 10 0.4 9  s +1 (s)
+ _ s Voltage s  4 Torque s 2 + 0.9s + 9

Ans. :

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM IX-4:
The block diagram of a control system for the human steering of an automobile is given.
a) Identify the basic factors for the open loop transfer function of the system.
b) Identify the corner frequencies separately for the magnitude and phase plots.
Remember that a corner frequency is defined as the intersection of two straight line
approximations !
c) Take Kp=1 and plot first the straight line approximations for the basic factors. Then
sketch the approximate Bode plots for the open loop transfer function of the system
by graphically adding the contributions of the basic factors. Mark all the lines and
make sure that the plots are clear and readable.
d) What should be the value of Kp if a magnitude of 0 [dB] is desired at a frequency of
1 [rad/s] ?

Desired Actual
Controller Position
Position Actuator Vehicle
(Driver)
R(s) 1 1 Y(s)
+ _
K p (s  0.1)
s s  10

Ans. : c)

d) 1000

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM IX-5:
Use Matlab to obtain the exact Bode plots for the open loop transfer function of the control
system in problem IX-2 and compare with the approximate plots. You can use the Matlab
command “bode(num, den)” which provides both the magnitude and phase plots. “num” and
“den” are the row vectors for the coefficients of the numerator and denominator polynomials
of the open loop transfer function. Use “help Bode” for additional information.

Ans. :
% Define s as the Laplace variable
s=tf('s');
% Define transfer function of the system
sys=12*(s+6)/(s*(s+1)*(s^2+6*s+36));
% Draw Bode diagram
BODE(sys)

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

X/ Sensitivity
PROBLEM X-1:

Consider a unity feedback control system with the feedforward transfer function

K
G(s) 
as  1

Determine the sensitivity of the closed loop transfer function T(s) with respect to the gain K
at steady state.

1
Ans. : STK ss 
1 K

PROBLEM X-2:
Consider a feedback control system with the feedforward and feedback transfer functions G(s)
and H(s), respectively.

Determine the sensitivity of the closed loop transfer function T(s) with respect to G(s) and
H(s). Can you reduce them both by choosing large controller gains, i.e. G(s)H(s)   ?

Ans. : No !

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

PROBLEM X-3:
Consider the unity feedback sytem with the block diagram shown.

a) Obtain the steady state sensitivity of the transfer function with respect to controller
gain Kp and the plant poles p1 and p2.

b) Calculate and plot, simultaneously, the magnitude of the dynamic sensitivity function
of the transfer function for this system with respect to controller gain Kp and the plant

poles p1 and p2, in the frequency range from 0 to 3 rad/s.

c) Considering the sensitivity function magnitudes at 2 rad/s, order the magnitude of the
sensitivities with respect to each parameter, from the least to the most.

Nominal values :
R(s) 10 C(s)
Kp = 10 Kp
p1 = -1
+ _ s  s  p1  s  p 2 
p2 = -5

Ans. : a) STK 
s(s  p1)(s  p2 ) , STp 
p1 s(s  p 2 ) 
,
p s(s  p1 )(s  p 2 )  10K p  1 s(s  p1 )(s  p 2 )  10K p 
   
p 2 s(s  p1 ) 
STp 
2 s(s  p1 )(s  p 2 )  10K p 
 

c)

Prof. Dr. Y. Samim Ünlüsoy

Mechanical Engineering Department Middle East Technical University
ME 304 CONTROL SYSTEMS
EXERCISES

XI/ ROOT LOCUS

PROBLEM XI-1:

The open loop transfer function of a system is given as:

K  2s  1
G(s) 
s(s  1)(s  2)(5s  1)

a) How many branches of the root locus are there ? Explain.

b) Where do they start and where do they end ?

Ans. : a) 4 branches, b) 1 finite zero at -0.5 and 3 zeroes at infinity. Root loci will start from
poles at 0, -1, -2, and -0.2. One will end at -0.5 and three at infinity.

Prof. Dr. Y. Samim Ünlüsoy