MITUCONTROL - Transient Response - LEC 9
MITUCONTROL - Transient Response - LEC 9
(TP 08 )
TRANSIENT RESPONSE
American University in
Cairo School of Science
and Engineering
Transient Response:
If the inputs to a control system are gradually changing functions of time, then a ramp
function of time may be a good test signal. Similarly, if a system is subjected to
sudden disturbances, a step function of time may be a good test signal; and for a
system subjected to shock inputs, an impulse function may be best. Once a control
system is designed on the basis of test signals, the performance of the system in
response to actual inputs is generally satisfactory. The use of such test signals enables
one to compare the performance of all systems on the same basis.
Effect of Pole & Zero on the System Response:
Consider the first-order system that has one real zero at -2 and one real pole at -5, as
shown below, therefore the open-loop T.F. is:
𝐶(𝑠) 𝑆 + 2
=
𝑅(𝑠) 𝑆 + 5
Considering a step input, i.e. R(s) = 1/S;
The first part of c(t) is called the forced response (it is independent on time), but the
second term is called natural response (it is dependent on time)
Another example: Consider the third-order system with open-loop T.F. is:
𝐶(𝑠) 𝑆+3
=
𝑅(𝑠) (𝑆 + 2)(𝑆 + 4)(𝑆 + 5)
2. Transient Response
1. First-Order system
Consider the first-order system shown in Fig. 2.
If we plot the output c(t) versus time, such kind of plot is dependent on the two
parameters ζ and ωn. A family of curves at different values of ζ is shown in Fig. 7.
𝑆2 + 2𝜁𝜔𝑛𝑆 + 𝜔𝑛2 = 0
Complete square of the above equation we get;
2
(𝑆 + 𝜁𝜔𝑛 )2 + (𝜔𝑛 √1 − 𝜁2) =0
To summarize our observations. The natural response has four different shapes:
1) Undamped Response: poles are pure imaginary at jn, Natural response:
Undamped sinusoid with radian frequency equal to the imaginary part of the poles,
𝐶(𝑡) = 𝐴𝑐𝑜𝑠(𝜔𝑛𝑡 − 𝜑)
2) Underdamped Response: poles are complex at d jd, Natural response: Damped
sinusoid with an exponential envelope whose time constant is equal to the reciprocal
of the pole’s real part. The radian frequency of the sinusoid, the damped frequency of
oscillation, is equal to the imaginary part of the poles,
𝐶(𝑡) = 𝐴 𝑒−𝜎𝑑𝑡 𝑐𝑜𝑠(𝜔𝑑𝑡 − 𝜑)
3) Critically-damped Response: poles are pure real and equal at −1, Natural
response: One term is an exponential whose time constant is equal to the reciprocal of
2. Rise Time: The rise time tr is the time required for the response to rise from 10%
to 90%. Or the time required to rise from 0% to 100% of its final value.
We obtain the rise time tr by letting c(tr) = 1
√1 − 𝜁2
tan(𝜔𝑑𝑡𝑟) = −
𝜁
1 −1
√1 − 𝜁2 𝜋−𝛽
𝑇𝑟 = 𝑡𝑎𝑛 (− ) =
𝜔𝑑 𝜁 𝜔𝑑
𝐶(𝑡 𝑝 ) − 𝐶(∞)
𝑀𝑎𝑥. %𝑜𝑣𝑒𝑟𝑠ℎ𝑜𝑜𝑡 (𝑀𝑝) = × 100%
𝐶(∞)
The maximum overshoot occurs at the peak time. Therefore
𝜋
−𝜁𝜔𝑛 𝜁
𝑀𝑝 = 𝐶(𝑡𝑝) − 1 = −𝑒 𝜔𝑑 (cos(𝜋) + sin(𝜋))
√1 − 𝜁2
𝜋𝜁
−
𝑀𝑝 = 𝑒 √1− 𝜁 2
The effect of damping ratio on the maximum overshoot is shown in Fig. 13.
5. Settling Time: The settling time ts is the time required for the response curve to
reach and stay within ± 2% or ± 5% of the final value. The settling time is the largest
time constant of the system.
In the previous Matlab code we consider some Matlab functions such as tf and step.
What is tf and step? and how can we use them?
Steady-State Error
The difference between the input and output of a system in the limit as time goes to
infinity, and it will be discussed in more details in next chapter.
Feedback PID controller – How does it work?
Consider the Mass (M), spring (K), and damper (B) problem given above;
𝑈(𝑠) 𝐾𝑝𝑆 + 𝐾𝑖
=
𝐸(𝑠) 𝑆
Find
So, we must calculate the angle β first based on Fig. below, as follows:
where
𝐾
𝜔𝑛 = √ (1)
𝐽
𝐵
2𝜁𝜔𝑛 = (2)
𝐽
But it is given that Ts = 2 for 2% criterion,
4
𝑇𝑠 = 2 = → 𝜁𝜔𝑛 = 2
𝜁𝜔 𝑛
From (2), B/J = 4
It is given that Mp = 0.2 → = 0.456
From (1), = 2J/K
Example #7
The T.F. of a closed-loop, unity feedback control system is
C(S) K
=
R(S) S2 + 2S + K
If the system gain (K) is set at three different values of 10, 36 and 100
- Calculate the rise time, maximum overshoot, and settling time at each value of K,
- At which value of K the system response is superior.
C(S) ωn2
= 2 2
R(S) S + 2ξω Sn + ω n
By Comparing,
n = K
= 1 / K
Β = cos-1 (1 / K)
d = n 1- 2 = √K(1-1/k) = √(K-1)
Based on 2% criteria, it is found that Ts = 4T,
At K = 10
n = 10 = 3.1623
= 1 / 10 = 0.31623
Β = cos-1 (1 / 10) = 71.56505 ̊ = 1.24904577 rad
d = K-1 = 3.0
𝜋− 𝛽 𝜋 − 1.24904577
Rise Time (Tr) = = = 0.63085 𝑠𝑒𝑐
𝜔𝑑 3
𝑀𝑎𝑥𝑖𝑚𝑢𝑚 𝑜𝑣𝑒𝑟𝑠ℎ𝑜𝑜𝑡 𝑀𝑝 = 𝑒−𝜋𝜉/√1−𝜉2 = 0.35085 = 35.085%
4 = 4 =4
𝑇𝑠 =
𝜉 𝜔𝑛 0.31623 × 3.1623
At K = 36
n = 36 = 6
= 1 / 6 = 0.16667
𝜋−𝛽 𝜋−
Rise Time (Tr) = 1.4334825
= = 0.2938 𝑠𝑒𝑐
𝜔𝑑 √35
𝜋− 𝛽 𝜋−
Rise Time (Tr) = = 1.47063 = 0.167938 𝑠𝑒𝑐
𝜔𝑑 √99
Example #8
A 3-term (PID) controller is used to control a process with unity feedback as shown
in Fig. 3, where Ti and Td are the integral and derivative time constant, respectively.
For unit step input,
a) If Td = 3.5, and the integral term is ignored, calculate the steady-state error,
b) If Ti = 2.0, and Td as given in (a), calculate the steady-state error,
c) Which steady-state error obtained from (a) and (b) is better. Why?
H(S) = 1,
80 (1 + 1 + 𝑇𝑑 𝑆)
𝑇𝑖𝑆
𝐺(𝑆) =
𝑆2 + 8𝑆 + 80
Since unit step input, we calculate the position error coefficient Kp
a) Ti is set to ∞ to ignore the integral term
Td = 3.5
80(1 + 3.5 𝑆)
𝐺(𝑆) =
𝑆2 + 8𝑆 + 80
80
𝐾𝑝 = lim 𝐺(𝑆) = =1
𝑆→0 80
1 = 1 = 0.5
𝐸𝑠𝑠 =
1 + 𝐾𝑝 1+1
b) Ti =2.0 and Td = 3.5
80 (1 + 1 + 3.5 𝑆)
𝐺(𝑆) = 2𝑆
2
𝑆 + 8𝑆 + 80
∞
𝐾𝑝 = lim 𝐺(𝑆) = =∞
𝑆→0 80
1 = 1 =0
𝐸𝑠𝑠 =
1 + 𝐾𝑝 1+∞
c) the steady-state error in case (b) is better than that of (a) because the integral term
is employed, therefore the system type is increased by one, so that the error is
reduced to 0.
d) Ti is set to ∞ to ignore the integral term
32 Dr. Samar Akef
Td is set to 0 to ignore the derivative term
The overall system is shown in the figure below
Process
R(S) + E(S) C(S)
80
_ 2
𝑆 + 8𝑆 + 80
𝐶(𝑆) 80
=
𝑅(𝑆) 𝑆2 + 8𝑆 + 160
The system characteristic equation is
𝑆2 + 8𝑆 + 160 = 0
The standard form of second order system characteristic equation is
𝑆2 + 2𝜁𝜔𝑛 𝑆 + 𝜔𝑛2 = 0
By comparing the coefficients
ωn = √160 = 12.649 rad/sec
2ξ ωn = 8 → ξ = 0.3162
𝜋𝜁
−
√1− 𝜁 2
Maximum overshoot = 𝑀𝑝 = 𝑒 = 0.35096 = 35.096 %
β = cos-1 0.3162 = 71.5667 = 1.2491 rad
Rise Time Tr
𝜋−𝛽 𝜋 − 1.2491
𝑇𝑟 = = = 0.1577 𝑠𝑒𝑐.
𝜔𝑑 12.649√1 − 0.31622
Peak Time Tp
𝜋 𝜋
𝑇𝑝 = = = 0.2618 𝑠𝑒𝑐.
𝜔𝑑 12.649√1 − 0.31622
Settling Time Ts
3 3
𝑇𝑠 = = 0.75 𝑠𝑒𝑐. (𝐵𝑎𝑠𝑒𝑑 𝑜𝑛 ± 5% 𝑡𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑒)
𝜁 𝜔𝑛 0.3162 × 12.649
=
4 4
𝑇𝑠 = = 1.0 𝑠𝑒𝑐. (𝐵𝑎𝑠𝑒𝑑 𝑜𝑛 ± 2% 𝑡𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑒)
𝜁 𝜔𝑛 0.3162 × 12.649
=