Chapter 6 - Transient Response Analysis
Chapter 6 - Transient Response Analysis
1. Introduction
It was stated previously in lecture #1 that the first step in analyzing a control system
was to derive a mathematical model of the system. Once such a model is obtained,
various methods are available for the analysis of system performance.
Typical Test Signals: The commonly used test input signals are those of step
functions, ramp functions, acceleration functions, impulse functions, sinusoidal
functions, and the like. With these test signals, mathematical and experimental
analyses of control systems can be carried out easily since the signals are very simple
functions of time.
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
Dr. AHMED MUSTAFA HUSSEIN
1 Instructor: Dr. Ahmed Kareem Al-Bakri
Al-Furat Al-Awsat Technical University Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Musiab F !"א#"א$%&
Technical College א
4th class Room
Jא
No. 1
Electrical./
Power0123&
Engineering K, Techniques אא
Email: [email protected]
אא
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.
The time response of a control system consists of two parts as shown in Fig. 1;
a) Transient response
b) Steady-state response.
where Ctr(t) is the transient response and Css(t) is the steady-state response.
The transient response of a practical control system often exhibits damped
oscillations before reaching a steady state. If the output of a system at steady state
does not exactly agree with the input, the system is said to have steady state error.
This error is indicative of the accuracy of the system. In analyzing a control system,
we must examine transient-response behavior and steady-state behavior.
2. Transient Response
2.1 First-Order system
Consider the first-order system shown in Fig. 2.
For a unit step input whose Laplace transform is 1/S, the output C(S) is given by
The above equation indicates that initially (at t = 0) the output c(t) is zero and finally
(at t = ∞) it becomes unity as shown in Fig. 3.
By the same way, in two time constants (t = 2T), the response reaches 86.5% of the
final value. At t = 3T, the response reaches 95% of its final value. At t = 4T, the
system response reaches 98.2% of its final value. Finally at t = 5T, the response
reaches 99.3% of the final value. Thus, for t ≥ 4T, the response remains within 2% of
the final value. As seen from the equation of c(t), the steady state value (c(t) = 1) is
reached mathematically only after an infinite time. In practice, however, a reasonable
estimate of the response time is the length of time the response curve needs to reach
and stay within the 2% line of the final value, or four time constants.
2.2 Second-Order Systems
Consider the 2nd order control system shown in Fig. 4, whose T.F. is given as:
This form is called the standard form of the second-order system, where ζ and ωn are
the damping ratio and undamped natural frequency, respectively.
This result can be obtained directly by using a table of Laplace transforms tables.
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. 5.
As the parameters ζ changes, the location of the system poles S1 and S2 are change.
Therefore, the dynamic behavior of the second-order system is also changes. The
nature of the roots s1 and s2 of the characteristic equation with varying values of
damping ratio ζ can be shown in the complex plane as shown in Fig. 6.
Delay Time: The delay time td is the time needed for the response to reach half (50%)
of its final value.
Delay time can be calculated from this formula;
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
Since ≠ 0, therefore
The cosine terms in the above equation cancel each other. Therefore, dc(t)/dt,
evaluated at t = tp, can be simplified to
This means
Since the peak time corresponds to the first peak overshoot, =
=
Settling Time: The settling time ts is the time required for the response curve to reach
and stay within ± 2% of the final value. In some cases, 5% instead of 2%, is used as
the percentage of the final value. The settling time is the largest time constant of the
system.
The settling time corresponding to ± 2% or ± 5% tolerance band may be measured in
terms of the time constant {T = l/ (ζ ωn)}
Dr. AHMED
Dr. AhmedMUSTAFA HUSSEIN
9 Instructor: Kareem Al-Bakri
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א J א
4th class Room No. 1
Al-Musiab Technical College
./
Electrical 0123&
Power K, Techniques
Engineering אא
Email: [email protected]
Summary:
Matlab software package can be used to calculate and plot the step response of the
second order system. The following m-file can be used.
Dr. AHMED
Instructor: MUSTAFA
Dr. Ahmed HUSSEIN
Kareem Al-Bakri
10
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א J א
4th class Room No. 1
Al-Musiab Technical College
./
Electrical Power0123& K, Techniques
Engineering אא
Email: [email protected]
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?
Dr. AHMED MUSTAFA HUSSEIN
11 Instructor: Dr. Ahmed Kareem Al-Bakri
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א J א
4th class Room No. 1
Al-Musiab Technical College
./
Electrical 0123&
Power K, Techniques
Engineering אא
Email: [email protected]
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.
As shown in the feedback control system given above, the type of controller used is
PID controller. The PID terms are stand for:
P: Proportional,
I: Integral,
D: Derivative
These correlations may not be exactly accurate, because Kp, Ki, and Kd are
dependent on each other. In fact, changing one of these variables can change the
effect of the other two.
Consider the Mass (m), spring (k), and damper (b) problem given in the above figure.
Dr. AHMED
Instructor: MUSTAFA
Dr. Ahmed HUSSEIN
Kareem Al-Bakri
14
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א J א
4th class Room No. 1
Al-Musiab Technical College
Electrical./
Power0123&
EngineeringK,Techniques אא
Email: [email protected]
From the system response shown above, the Mass-spring and damper system, is
suffering from the following problems:
Rise time is improved (Tr=0.1) and steady-state error is improved (Ess=0.95) but the
system overshoot is deteriorated (Mp~1.1). Settling time (Ts=1.2)
Dr. AHMED
Instructor: MUSTAFA
Dr. Ahmed Kareem HUSSEIN
Al-Bakri
16
Al-Furat Al-Awsat Technical University Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١F !"א#"א$%&
Al-Musiab Technical College
א
4th class Room
Jא
No. 1
Electrical./
Power0123& K, Techniques
Engineering אא
Email: [email protected]
Rise time and steady-state error are not affected. But the system overshoot is
improved (Mp~1.05) and settling time is improved (Ts~0.5)
Dr. AHMED MUSTAFA HUSSEIN
17 Instructor: Dr. Ahmed Kareem Al-Bakri
Al-Furat Al-Awsat Technical University Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١F !"א#"א$%&
Al-Musiab Technical College
א
4th class Room
Jא
No. 1
./
Electrical 0123&
Power K, Techniques
Engineering אא
Email: [email protected]
Example #1
Consider the system shown in Fig. 9, where ζ = 0.6 and ωn = 5 rad/sec. Let us obtain
the rise time tr, peak time tp, maximum overshoot Mp, and settling time tp when the
system is subjected to a unit-step input.
Fig. 9.
Dr. AHMED
Instructor: Dr.MUSTAFA HUSSEIN
Ahmed Kareem Al-Bakri
20
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א
4th class Room
Jא
No. 1
Al-Musiab Technical College
Electrical./
Power0123& K, Techniques
Engineering אא
Email: [email protected]
Example #2
Consider the control system whose closed loop poles are given in Fig. 10.
Fig. 10.
Find
Example #3
Determine the values of Td, Tr, Tp and Ts for the control system shown in Fig. 11.
Fig. 11
Dr. AHMED MUSTAFA HUSSEIN
21 Instructor: Dr. Ahmed Kareem Al-Bakri
Al-Furat Al-Awsat Technical University Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١F !"א#"א$%&
Al-Musiab Technical College
א
4th class JRoom
א
No. 1
./
Electrical 0123&
Power K, Techniques
Engineering אא
Email: [email protected]
Fig. 12
Dr. AHMED
Instructor: MUSTAFA
Dr. Ahmed HUSSEIN
Kareem Al-Bakri
22
Al-Furat Al-Awsat Technical University Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١F !"א#"א$%&
Al-Musiab Technical College
א
4th class Room
Jא
No. 1
Electrical./
Power 0123& K,Techniques
Engineering אא
Email: [email protected]
Example #4
For the system shown in Fig. 13, 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 =
1 kg-m2 and B = 1 N-m/rad/sec.
Since Mp = 0.2;
where
Example #5
When the system shown in Fig. 14 (a) is subjected to a unit-step input, the system
output responds as shown in Fig.14 (b). Determine the values of K and T from the
response curve.
Report:
Determine the values, of K and k of the closed-loop system shown in Fig. so that the
maximum overshoot in unit-step response is 25% and the peak time is 2 sec. Assume
that J = 1 kg-m2.
Example #6
The T.F. of a closed-loop, unity feedback control system is
C(S) K
=
R(S) S + 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) ω(
=
R(S) S + 2ξω( S + ω(
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,
Dr. AHMED
Instructor: MUSTAFA
Dr. Ahmed HUSSEIN
Kareem Al-Bakri
26
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א J א
4th class Room No. 1
Al-Musiab Technical College
./
Electrical 0123&
Power K, Techniques
Engineering אא
Email: [email protected]
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
− 0 − 1.24904577
Rise Time (Tr) = = = 0.63085 7 8
3
9:;<=< >? @7ℎ>> B = C/ C
= 0.35085 = 35.085%
4 4
FG = = =4
H 0.31623 × 3.1623
At K = 36
ωn = √36 = 6
ξ = 1 / 6 = 0.16667
Β = cos-1 (1 / 6) = 80.40593177 ̊ = 1.4334825 rad
ωd = √35 = 5.9160798
− 0 − 1.4334825
Rise Time (Tr) = = = 0.2938 7 8
√35
4 4
FG = = =47 8
H 0.16667 × 6
At K = 100
ωn = √100 = 10
ξ = 1 / 10 = 0.1
Β = cos-1 (0.1) = 84.261 ̊ = 1.47063 rad
ωd = √99 = 9.94987
− 0 − 1.47063
Rise Time (Tr) = = = 0.167938 7 8
√99
4 4
FG = = =47 8
H 0.1 × 10
Rise Time Maximum Overshoot Settling Time
10 0.63085 7 8 35.085% 47 8
36 0.2938 7 8 58.8% 47 8
100 0.167938 7 8 72.925% 47 8
Dr. AHMED MUSTAFA
Kareem HUSSEIN
27 Instructor: Dr. Ahmed Al-Bakri
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א J א
4th class Room No. 1
Al-Musiab Technical College
./
Electrical Power0123& K, Techniques
Engineering אא
Email: [email protected]
Based on information given in the table, by increasing the system gain from 10 to 100, the
rise time and steady-state error are decreased (improved) which is V.Good. On the other
hand, the Maximum overshoot is increased (deteriorated).
Example #7
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?
d) If both derivative and integral terms are ignored, calculate the damping ratio,
maximum overshoot, rise time, peak time and settling time, then draw a free-hand
sketch for the system output c(t).
H(S) = 1,
1
80 L1 + + F N
FM
K() =
+ 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 )
K() =
+ 8 + 80
80
O = lim K() = =1
Q→S 80
1 1
TGG = = = 0.5
1 + O 1+1
Dr. AHMED
Instructor: Dr. AhmedMUSTAFA HUSSEIN
Kareem Al-Bakri
28
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א
4th class Room
Jא
No. 1
Al-Musiab Technical College
Electrical./
Power0123&
EngineeringK,Techniques אא
Email: [email protected]
Process
R(S) + E(S) C(S)
80
_
+ 8 + 80
V() 80
=
W() + 8 + 160
The system characteristic equation is
+ 8 + 160 = 0
The standard form of second order system characteristic equation is
+ 2 + = 0
By comparing the coefficients
ωn = √160 = 12.649 rad/sec
2ξ ωn = 8 → ξ = 0.3162
XY
Z[\ Y
Maximum overshoot = = = 0.35096 = 35.096 %
β = cos-1 0.3162 = 71.5667° = 1.2491 rad
Dr. AHMED
Dr. AhmedMUSTAFA HUSSEIN
29 Instructor: Kareem Al-Bakri
Instructor: Dr. Ahmed Kareem Al-Bakri
E٣٥١١
Al-Furat F !"א#"א$%&
Al-Awsat Technical University א J א
4th class Room No. 1
Al-Musiab Technical College
./
Electrical Power0123& K, Techniques
Engineering אא
Email: [email protected]
Rise Time Tr
−0 − 1.2491
F] = = = 0.1577 7 8.
12.649√1 − 0.3162
Peak Time Tp
F = = = 0.2618 7 8.
12.649√1 − 0.3162
Settling Time Ts
3 3
FG = = = 0.75 7 8. (^97 _ >` ± 5% >b @9`8 )
0.3162 × 12.649
4 4
FG = = = 1.0 7 8. (^97 _ >` ± 2% >b @9`8 )
0.3162 × 12.649
Dr. AHMED
Instructor: MUSTAFA
Dr. Ahmed HUSSEIN
Kareem Al-Bakri
30