اسس هندسة السيطرة النظري-مرحلة ثالثة-قسم هندسة تقنيات الحاسوب
اسس هندسة السيطرة النظري-مرحلة ثالثة-قسم هندسة تقنيات الحاسوب
Control: Means measuring the value of the controlled variable of the system and
applying the manipulated variable to the system to correct or limit deviation of the
measured value from a desired value.
Plants: A plant may be a piece of equipment, perhaps just a set of machine parts
functioning together, the purpose of which is to perform a particular operation.
The output of the system has to be measured. In the figure below, we show
the system we are trying to control -the "plant"- and a sensor that measures
what the controlled system is doing.
The input to the plant is usually called the control effort, and the output of the
sensor is usually called the measured output, as shown below in the figure.
If we want to control the output, we first need to measure the output. Within the
whole system is the system we want to control - the plant - along with
a sensor that measures what the output actually is.
In our block diagram representation, we show the output signal being fed to
the sensor which produces another signal that is dependent upon the output.
We need the sensor in the system to measure what the system is doing.
To control the system we need to use the information provided by the sensor.
Usually, the output, as measured by the sensor is subtracted from the input
(which is the desired output) as shown below. That forms an error signal that
the controller can use to control the plant.
Let us take an example of simple temperature controller of the room, to clear the
concept. Suppose there is a simple heating element, which is heated up as long as
the electric power supply is switched on. As long as the power supply switch of the
heater is on the temperature of the room rises and after achieving the desired
temperature of the room, the power supply is switched off. Again due to ambient
temperature, the room temperature falls and then manually the heater element is
switched on to achieve the desired room temperature again. In this way one can
manually control the room temperature at desired level. This is an example of
manual control system.
The above system can further be improved by using timer switching arrangement
of the power supply where the supply to the heating element is switched on and off
in a predetermined interval to achieve desired temperature level of the room. There
is another improved way of controlling the temperature of the room. Here one
sensor measures the difference between actual temperature and desired
temperature. If there is any difference between them, the heating element functions
to reduce the difference and when the difference becomes lower than a
predetermined level, the heating elements stop functioning. Both forms of the
system are automatic control system.
In former one the input of the system is entirely independent of the output of the
system. Temperature of the room (output) increases as long as the power supply
switch is kept on. That means heating element produces heat as long as the power
supply is kept on and final room temperature does not have any control to the input
power supply of the system. This system is referred as open loop control system.
3 ` Dr. Mohammed Abdulla Abdulsada
But in the latter case, the heating elements of the system function, depending upon
the difference between, actual temperature and desired temperature. This
difference is called error of the system. This error signal is fed back to the system
to control the input. As the input to output path and the error feedback path create a
closed loop, this type of control system is referred as closed loop control system.
Hence, there are two main types of control system. They are as follow open loop
control system and closed loop control system.
1. Electric Hand Drier – Hot air (output) comes out as long as you keep your hand
under the machine, irrespective of how much your hand is dried.
2. Automatic Washing Machine – This machine runs according to the pre-set time
irrespective of washing is completed or not.
3. Bread Toaster - This machine runs as per adjusted time irrespective of toasting is
completed or not.
4. Automatic Tea/Coffee Maker – These machines also function for pre adjusted
time only.
5. Timer Based Clothes Drier – This machine dries wet clothes for pre – adjusted
time, it does not matter how much the clothes are dried.
2. Economical.
3. Easy to maintain.
4. Generally stable.
3. Water Level Controller– Input water is controlled by water level of the reservoir.
1. Closed loop control systems are more accurate even in the presence of non-
linearity.
4. Facilitates automation.
5. The sensitivity of system may be made small to make system more stable.
6. Stability is the major problem and more care is needed to design a stable closed
loop system.
10 Examples: Hand drier, tea maker Examples: Servo voltage stabilizer, air
conditioner
( )
( )
( )
( )
( )
( ) ( ) ( )
( ) ( )[ ( ) ( ) ( )]
( ) ( ) ( ) ( ) ( ) ( )
( )[ ( ) ( )] ( ) ( )
( ) ( )
( ) ( ) ( )
Thus the output of the closed-loop system clearly depends on both the closed-loop
transfer function and the nature of the input.
When the H(s) = 1 or unity the system is called negative unity FB control system.
-
B (S)
Black's Formula
It is a way of achieving close loop TF directly, without any equation
writing. It can be directly applied to general feedback configuration as
follows:
EX1:
( ) ( )
negative feedback)
( ) ( )
( ( ) ( ) ( ))
( ) ( )
( ) ( ) ( ))
+ +
A(S) B(S)
ii. Feed forward transfer function = A(s)
A(s) B(s)
A(s) A(s) B(s)
iii. The overall transfer function =
A(s) B(s) A(s) A(𝑆)𝐵(𝑆)𝐶(𝑆)
∙ 𝐶(𝑆)
A(s)
or by using Black's formula with 2 loop gains:
( ) ( )
( ) ( ) ( ) ( )
N(S): Disturbance
In this case, one of the inputs, D(s), is known as disturbance, while R(s) is the
reference input. Before designing a proper controller for the system, it is always
important to learn the effects of D(s) on the system.
With zero disturbance ( N(s) =0), we may calculate the response CR(S)
to the reference input only:
The fundamental law govering mechanical system is Newton’s second law. The
variables that are used to describe translational motion are acceleration, velocity,
and displacement.
Newton's law of motion states that the algebraic sum of external forces acting on a
rigid body in a given direction is equal to the product of the mass of the body and
its acceleration in the same direction.
Σ forces = Ma
where M denotes the mass, and a is the acceleration in the direction considered.
Fig. 1 illustrates the situation where a force is acting on a body with mass M. The
force equation is written as
where a(t) is the acceleration, v(t) denotes linear velocity, and y(t) is the
displacement of mass M, respectively.
𝒅𝒚
𝒇(𝒕) 𝑩
𝒅𝒕
where : B is frictional coefficient.
Any mechanical element that undergoes a change in shape (size) when subjected to
a force can be characterized as a spring.
f(t) = K. y(t)
where K is the spring constant, or simply stiffness. Eq. above implies that the force
acting on the spring is directly proportional to the displacement (deformation) of
the spring.
Eq.1
The last equation may be rearranged by equating the highest-order derivative term
to the rest of the terms:
Eq.2
Eq.3
The transfer function between Y(s) and F(s) is obtained by taking the Laplace
transform on both sides of Eq. (3) with zero initial conditions:
The voltage drop across a resistor R , inductor L and capacitor C are given as
follows:
( )
∫ ( )
To draw a block diagram for an electrical system, first write the equations that
describe the dynamic behavior of each component. Then take the Laplace
transform of equation individually in block form. Finally, assemble the elements
into a complete block diagram.
( )
∫ ( ) ( )
Assembling these two elements, we obtain the overall diagram in S-Domain for the
system.
ei eo
( )
( )
( )
( ) ( )
LS
ei eo
The thermal system shown below in Fig. a will has the block diagram shown in
Fig.b.
Where :
EX1: Consider the system shown below, simplify this diagram by using block
diagram reduction techniques?
Combining two blocks into one and finally, eliminating the unity feedback loop
results in following figure:
First, move the branch point of the path involving H1 outside the loop involving
H2, as in following figure:
By moving the summing point of the negative feedback loop containing H2 outside
the positive feedback loop containing H1, we obtain the following figure:
Eliminating the positive feedback loop of H1, we have the following figure:
Homwork2: Consider the system shown below, simplify this diagram and find the
system transfer function?
( )
( )
Most linear closed-loop systems have closed-loop transfer functions of the form:
2. If any of the coefficients are zero or negative in the presence of at least one
positive coefficient, a root or roots exist that are imaginary or that have positive
real parts. Therefore, in such a case, the system is not stable.
3. If all coefficients are positive, arrange the coefficients of the polynomial in rows
and columns according to the following pattern:
Complex Imaginary
Conjugate Axis
Poles jw Complex
Poles
Real
-2 -1 0 Axis
Real σ
Pole Real
Pole
LHP RHP
The evaluation of the b’s is continued until the remaining ones are all zero. The
same pattern of cross-multiplying the coefficients of the two previous rows is
followed in evaluating the c’s, d’s, e’s, and so on. That is,
This process is continued until the nth row has been completed. The complete
array of coefficients is triangular.
The necessary and sufficient condition that all roots of characteristic equation lie in
the left-half s plane is that all the coefficients are be positive and all terms in the
first column of the array have positive signs.
Each sign change counts as one pole in the right half of the complex plane (RHP)
and, as mentioned earlier, poles in the RHP make the system unstable. Therefore,
you want all the elements in the first column to be positive and if you get that, your
poles are safely in the LHP and your system will be stable.
Using the Routh stability criterion, determine the stability of the system?
S4 1 3 5
𝟐∗𝟑 𝟒∗𝟏
S3 2 4 0 𝟐
𝟐∗𝟓 𝟎∗𝟏
S2 1 5
𝟐
𝟏∗𝟒 𝟓∗𝟐
S1 -6
𝟏
0
S 5
𝟔∗𝟓 𝟎∗𝟏
𝟔
The number of roots in the right half-plane (RHP) is equal to the number of sign
changes in the first column of Routh array. Since there are 2 sign changes therefore
the system has 2 poles in the RHP and hence the system is unstable.
Special Cases
If a first-column term in any row is zero, but the remaining terms are not zero or
there is no remaining term, then the zero term is replaced by a very small positive
number and the rest of the array is evaluated.
If the sign of the coefficient above the zero ( ) is the same as that below it, then it
indicates that there are a pair of imaginary roots and the system is critical stable or
oscillatory.
If however, the sign of the coefficient above the zero is opposite that below it, it
indicates that there is one sign change and the system is unstable.
The system is critical stable or oscillatory since the sign of the coefficient above
the zero ( ) is the same as that below it. The system has a pair of imaginary roots
at +j and –j ( (S+2)(S-j)(S+j) ).
Using the Routh stability criterion, determine the stability of the system?
Special case: replace the zero by
very small no. ( 𝛜 )
There are two sign changes of the coefficients in the first column which indicated
of 2 poles in RHP and hence the system is unstable.
( )
S5 1 2 4
Special case : replace the zero by
S4 2 4 10
very small no. ( 𝛜 )
S3 0≈ϵ -1 0
S2 10
S1 ≈-1
S0 10
There are two sign changes of the coefficients in the first column and hence the
system is unstable.
Consider the system shown in figure below. Let us determine the range of K for
stability.
For stability, K must be positive ( K>0 ), and all coefficients in the first column
must be positive ( ). Therefore,
When K=14/9 , the system becomes oscillatory and, mathematically the oscillation
is sustained at constant amplitude.
( )( )
S3 1 2
S2 3 K
S1
S0 K
For stability all coefficients in the first column must be positive, therefore K must
be positive ( K > 0 ), and ( ). Therefore,
6 > K >0
EX6: Consider the following characteristic equation:
Using the Routh stability criterion, determine the range of K for stability?
S4 1 4+K 25
S3 2 9 0
S1
S0 25
For stability all coefficients in the first column must be positive, therefore ( K> ),
K> ( i. e. K>6.056 )
* Another method that the branches of root loci on the real axis lie to the left of
an odd number of finite poles and zeros.
3. Determine the asymptotes of the root-locus branches.
* No. of asymptotes = No. of poles (n) – No. of zeros (m)
( )
* Angles of asymptotes =
∑ ∑
* Intersection point of asymptotes with real axis σc=
( Pi and Zi : numerical values of poles and zeros respectively)
4. Find the breakaway and break-in pints.
* From the characteristic Equation ( 1 + G(S)H(S) = 0 ) find the equation of gain
parameter (K).
* The breakaway points can be simply determined from the roots of dK/ds=0.
* It should be noted that not all the solutions of Equation dK/ds=0 correspond
to actual breakaway or break-in points. If a point at which dK/ds=0 is on a root
locus, it is an actual breakaway or break-in point.
5. Find the points where the root loci may cross the imaginary axis.
* These points can be found by use of Routh’s stability criterion and find value
of K that makes the term in the first column equal zero.
32 ` Dr. Mohammed Abdulla Abdulsada
* The crossing points on the imaginary axis can then be found by solving the
auxiliary equation obtained from the one upper row.
* an alternative approach to find these points by let S=jw in the characteristic
equation, equate both the real part and the imaginary part to zero, and then solve
for w and K.
6. Taking a series of test points in the bread neighborhood of the origin of the s
plane, and finally sketch the root loci.
Note : The root-locus branches start from open-loop poles and terminate at
zeros (finite zeros or zeros at infinity).
EX.1: Draw the root locus for the unity negative feedback system has a gain and
transfer functions as follows?
( ) ( )
( )
SOL.: Imaginary
Axis
(1) Poles (S1 = 0 , S2 = -3 ) , No finite Zeros.
jw
(2) Branches on real axis : [ ∞ 0] No branch
[0 -3] branch
-3 0
[ -3 -∞ ] No branch Real Axis σ
(3) No. of asymptotes = No. of poles (n) – No. of zeros (m)
=2–0=2
( )
Angles of asymptotes =
= ± 90 (k=0)
= ± 270 (k=1) = ± 90
Therefore angles of asymptotes will be +90o and -90o .
∑ ∑
Intersection point of asymptotes with real axis σc=
( )
S(S+3) + K =0 -3 0
-90o Real Axis σ
2
K = - S(S+3) = - S - 3S
-1.5
-3 0
EX.2: Consider the negative feedback system shown in Figure below, sketch the
root-locus plot?
SOL.:
1 + G(S)H(S) = 0
Since the breakaway point must lie on a root locus between 0 and –1, it is clear that
s=-0.4226 corresponds to the actual breakaway point. Point s=-1.5774 is not on the
root locus. Hence, this point is not an actual breakaway or break-in point.
5. Determine the points where the root loci cross the imaginary axis. These points
can be found by use of Routh’s stability criterion as follows: Since the
characteristic equation for the present system is:
which yields:
The frequencies at the crossing points on the imaginary axis are thus √
.The gain value corresponding to the crossing points is K=6.
( ) ( )
( )
( )( )
(1) Poles (S1 = 2 & S2 = -1 ) , Zeros (S1= -3 )
(2) Branches on real axis : [ ∞ 2] No branch
[2 -1] branch
[ -1 -3 ] No branch
[ -3 -∞] branch
(3) No. of asymptotes = No. of poles (n) – No. of zeros (m)
=2–1=1
( )
Angles of asymptotes =
= ± 180 (k=0)
= ± 540 (k=1) = ± 180
Therefore angles of asymptotes will be +180o (or -180 o).
Since the asymptotes has 180o angle , therefore there are no intersection point
with the real axis ( i. e. parallel with real axis).
Imaginary
Axis
jw
-3 -1 0 2
Real Axis σ
1 + G(S)H(S) = 0
( )
( )
( )
( )
( )
5. Determine the points where the root loci cross the imaginary axis.
( ) ( ) ( )
S2 1 3K - 2
S1 K-1 0
S0 3K - 2
The value of K that makes the S1 term in the first column equal zero is K=1.The
crossing points on the imaginary axis can then be found by solving the auxiliary
equation obtained from the S2 row; that is,
( )
S2 = -1
-3 -1 2
EX.4: Consider the negative feedback system shown in Figure below, sketch the
root-locus plot?
1 + G(S)H(S) = 0
( )
( )
( )
( )
Now apply dK/ds=0
( )
To determine accurate root loci, several points must be found by trial and error
between the break-in point and the complex open-loop poles. For Example by
taking K= 2 , S= -2 ± j1.73 ; for K=4 , S= -3 ± j1.4 and for K= 5 , S= -3.5 ± j0.866.
-3.732
j1
-2 -1
-j1
1. Transient Response
The element of the system that vanishes with time is called transient response. It is
the response of the system from rest or equilibrium to steady state. It is denoted by
ctr(t).
2. Steady-State Response
It is the element of system response that is achieved when the time approaches to
infinity. The response of the system after transient response is basically called
steady state response. It is denoted by css(t).
Transient response doesn’t depend upon input of systems therefore can be analyzed
using step input. Steady state depends upon input and dynamics of system and can
be determined using different test signals by final value theorem.
The equation of time response will be:
c(t)= ctr(t) + css(t)
where,
1. Impulse Signal:
Impulse response in control system imitates sudden shock quality of actual input
signal. It is the reply of the system to the direct delta input.
When A=1 then the impulse signal is called Unit impulse signal.
The Laplace Transform of the impulse signal is R(s) = 1.
3. Ramp Signal:
The ramp signal tells you the constant velocity attribute of actual input signal. It is
being used to determine the behavior of system with the velocity factor.
Since the Laplace transform of the unit-step function is 1/s, substituting R(s)=1/s
into above Eq. we obtain:
The last Eq. states that initially the output c(t) is zero and finally it becomes unity.
One important characteristic of such an exponential response curve c(t) is that at
t=T the value of c(t) is 0.632, or the response c(t) has reached 63.2% of its total
change.
The unit-ramp input and the system output are shown in below Figure.
Where:
ωn is the undamped natural frequency of the system.
ζ is the damping ratio of the system.
ζ . ωn is called the attenuation ( σ).
The dynamic behavior of the second-order system can then be described in terms
of two parameters ζ and ωn.
If 0< ζ <1, the closed-loop poles are complex conjugates and lie in the left-half S
plane. The system is then called underdamped, and the transient response is
oscillatory.
If ζ =0, the transient response does not die out.
If ζ =1, the system is called critically damped.
Overdamped systems correspond to ζ >1.
A family of unit-step response curves c(t) with various values of ζ is shown below.
Transient-Response Specifications
The transient response of a practical control system often exhibits damped
oscillations before reaching steady state. In specifying the transient-response
characteristics of a control system to a unit-step input, it is common to specify the
following:
These specifications are defined in what follows and are shown graphically in
Figure.
1. Delay time (td): The delay time is the time required for the response to reach
half the final value the very first time.
2. Rise time (tr) : The rise time is the time required for the response to rise from
10% to 90%, 5% to 95%, or 0% to 100% of its final value. For underdamped
second order systems, the 0% to 100% rise time is normally used. For overdamped
systems, the 10% to 90% rise time is commonly used.
3. Peak time (tp) :The peak time is the time required for the response to reach the
first peak of the overshoot.
The amount of the maximum (percent) overshoot directly indicates the relative
stability of the system.
5. Settling time (ts) : The settling time is the time required for the response curve
to reach and stay within a range about the final value of size specified by absolute
percentage of the final value (usually 2% or 5%). The settling time is related to the
largest time constant of the control system.
The time-domain specifications just given are quite important, since most control
systems are time-domain systems; that is, they must exhibit acceptable time
responses. (This means that, the control system must be modified until the transient
response is satisfactory.)
attenuation = σ = ζ . ωn
rise time =
angle ( rad )
peak time =
( )
( ) √
maximum overshoot = Mp =
= for 5 % criterion
EX.1 : Consider the system shown in Figure below, where ζ=0.6 and ωn=5 rad-sec.
calculate the rise time tr , peak time tp, maximum overshoot Mp, and settling time
ts when the system is subjected to a unit-step input?
( ) ⁄
( )
EX.3: For the system shown in Figure (a), determine the values of gain K and
velocity-feedback constant Kh so that the maximum overshoot in the unit-step
response is 0.2 and the peak time is 1 sec?
With these values of K and Kh, obtain the rise time and settling time. Assume that
J=1kg-m2 and B=1 N-m_rad_sec?
After simplification, the block diagram of the system will be as shown in Figure
(b) below
( )
√
maximum overshoot = Mp =
The various types of controllers are used to improve the performance of control
systems.
2-Unconventional Controllers
Regardless of the size of the error, the output of the on-off controller can only be
fully ON or fully OFF, it is not proportional with the error.
►PID Controller improve the stability of the control system and to decrease
steady state error.
1) without controller?
R (S) + C (S)
𝑆 𝑆
- G(S)
1. without controller: The Simulink model and scope output is shown below
ts ≈ 13 Sec
ess ≈ 0
The system without controller has poor transient responses since the Mp and ts are
large values.
Mp= 0.75 = 75 %
Mp= 0.01 = 1 %
The system with PD controller give good transient
ts ≈ 1.75 Sec
responses since the Mp and ts are decreased
ess ≈ 0
Mp= 0.01 = 1 %
ts ≈ 1 Sec
The system with PID controller enhances transient
responses since the Mp and ts have very small values.
ess ≈ 0
80 ` Dr. Mohammed Abdulla Abdulsada