
In this chapter, let us discuss the time domain specifications of the second order system. The step response of the second order system for the underdamped case is shown in the following figure.
All the time domain specifications are represented in this figure. The response up to the settling time is known as transient response and the response after the settling time is known as steady state response.
It is the time required for the response to reach half of its final value from the zero instant. It is denoted by td.
Consider the step response of the second order system for t ≥ 0, when ‘δ’ lies between zero and one.
c(t)=1−(e−δωnt√1−δ2)sin(ωdt+θ)
The final value of the step response is one.
Therefore, at t=td, the value of the step response will be 0.5. Substitute, these values in the above equation.
c(td)=0.5=1−(e−δωntd√1−δ2)sin(ωdtd+θ)
⇒(e−δωntd√1−δ2)sin(ωdtd+θ)=0.5
By using linear approximation, you will get the delay time td as
td=1+0.7δωn
It is the time required for the response to rise from 0% to 100% of its final value. This is applicable for the under-damped systems. For the over-damped systems, consider the duration from 10% to 90% of the final value. Rise time is denoted by tr.
At t = t1 = 0, c(t) = 0.
We know that the final value of the step response is one.
Therefore, at t=t2, the value of step response is one. Substitute, these values in the following equation.
c(t)=1−(e−δωnt√1−δ2)sin(ωdt+θ)
c(t2)=1=1−(e−δωnt2√1−δ2)sin(ωdt2+θ)
⇒(e−δωnt2√1−δ2)sin(ωdt2+θ)=0
⇒sin(ωdt2+θ)=0
⇒ωdt2+θ=π
⇒t2=π−θωd
Substitute t1 and t2 values in the following equation of rise time,
tr=t2−t1
∴
From above equation, we can conclude that the rise time t_r and the damped frequency \omega_d are inversely proportional to each other.
It is the time required for the response to reach the peak value for the first time. It is denoted by t_p. At t = t_p, the first derivate of the response is zero.
We know the step response of second order system for under-damped case is
c(t)=1-\left ( \frac{e^{-\delta \omega_nt}}{\sqrt{1-\delta^2}} \right )\sin(\omega_dt+\theta)
Differentiate c(t) with respect to ‘t’.
\frac{\text{d}c(t)}{\text{d}t}=-\left ( \frac{e^{-\delta\omega_nt}}{\sqrt{1-\delta^2}} \right )\omega_d\cos(\omega_dt+\theta)-\left ( \frac{-\delta\omega_ne^{-\delta\omega_nt}}{\sqrt{1-\delta^2}} \right )\sin(\omega_dt+\theta)
Substitute, t=t_p and \frac{\text{d}c(t)}{\text{d}t}=0 in the above equation.
0=-\left ( \frac{e^{-\delta\omega_nt_p}}{\sqrt{1-\delta^2}} \right )\left [ \omega_d\cos(\omega_dt_p+\theta)-\delta\omega_n\sin(\omega_dt_p+\theta) \right ]
\Rightarrow \omega_n\sqrt{1-\delta^2}\cos(\omega_dt_p+\theta)-\delta\omega_n\sin(\omega_dt_p+\theta)=0
\Rightarrow \sqrt{1-\delta^2}\cos(\omega_dt_p+\theta)-\delta\sin(\omega_dt_p+\theta)=0
\Rightarrow \sin(\theta)\cos(\omega_dt_p+\theta)-\cos(\theta)\sin(\omega_dt_p+\theta)=0
\Rightarrow \sin(\theta-\omega_dt_p-\theta)=0
\Rightarrow sin(-\omega_dt_p)=0\Rightarrow -\sin(\omega_dt_p)=0\Rightarrow sin(\omega_dt_p)=0
\Rightarrow \omega_dt_p=\pi
\Rightarrow t_p=\frac{\pi}{\omega_d}
From the above equation, we can conclude that the peak time t_p and the damped frequency \omega_d are inversely proportional to each other.
Peak overshoot Mp is defined as the deviation of the response at peak time from the final value of response. It is also called the maximum overshoot.
Mathematically, we can write it as
M_p=c(t_p)-c(\infty)
Where,
c(tp) is the peak value of the response.
c(∞) is the final (steady state) value of the response.
At t = t_p, the response c(t) is -
c(t_p)=1-\left ( \frac{e^{-\delta\omega_nt_p}}{\sqrt{1-\delta^2}} \right )\sin(\omega_dt_p+\theta)
Substitute, t_p=\frac{\pi}{\omega_d} in the right hand side of the above equation.
c(t_P)=1-\left ( \frac{e^{-\delta\omega_n\left ( \frac{\pi}{\omega_d} \right )}}{\sqrt{1-\delta^2}} \right )\sin\left ( \omega_d\left ( \frac{\pi}{\omega_d} \right ) +\theta\right )
\Rightarrow c(t_p)=1-\left ( \frac{e^{-\left ( \frac{\delta\pi}{\sqrt{1-\delta^2}} \right )}}{\sqrt{1-\delta^2}} \right )(-\sin(\theta))
We know that
\sin(\theta)=\sqrt{1-\delta^2}
So, we will get c(t_p) as
c(t_p)=1+e^{-\left ( \frac{\delta\pi}{\sqrt{1-\delta^2}} \right )}
Substitute the values of c(t_p) and c(\infty) in the peak overshoot equation.
M_p=1+e^{-\left ( \frac{\delta\pi}{\sqrt{1-\delta^2}} \right )}-1
\Rightarrow M_p=e^{-\left ( \frac{\delta\pi}{\sqrt{1-\delta^2}} \right )}
Percentage of peak overshoot % M_p can be calculated by using this formula.
\%M_p=\frac{M_p}{c(\infty )}\times 100\%
By substituting the values of M_p and c(\infty) in above formula, we will get the Percentage of the peak overshoot \%M_p as
\%M_p=\left ( e^ {-\left ( \frac{\delta\pi}{\sqrt{1-\delta^2}} \right )} \right )\times 100\%
From the above equation, we can conclude that the percentage of peak overshoot \% M_p will decrease if the damping ratio \delta increases.
It is the time required for the response to reach the steady state and stay within the specified tolerance bands around the final value. In general, the tolerance bands are 2% and 5%. The settling time is denoted by t_s.
The settling time for 5% tolerance band is -
t_s=\frac{3}{\delta\omega_n}=3\tau
The settling time for 2% tolerance band is -
t_s=\frac{4}{\delta\omega_n}=4\tau
Where, \tau is the time constant and is equal to \frac{1}{\delta\omega_n}.
Both the settling time t_s and the time constant \tau are inversely proportional to the damping ratio \delta.
Both the settling time t_s and the time constant \tau are independent of the system gain. That means even the system gain changes, the settling time t_s and time constant \tau will never change.
Let us now find the time domain specifications of a control system having the closed loop transfer function \frac{4}{s^2+2s+4} when the unit step signal is applied as an input to this control system.
We know that the standard form of the transfer function of the second order closed loop control system as
\frac{\omega_n^2}{s^2+2\delta\omega_ns+\omega_n^2}
By equating these two transfer functions, we will get the un-damped natural frequency \omega_n as 2 rad/sec and the damping ratio \delta as 0.5.
We know the formula for damped frequency \omega_d as
\omega_d=\omega_n\sqrt{1-\delta^2}
Substitute, \omega_n and \delta values in the above formula.
\Rightarrow \omega_d=2\sqrt{1-(0.5)^2}
\Rightarrow \omega_d=1.732 \: rad/sec
Substitute, \delta value in following relation
\theta=\cos^{-1}\delta
\Rightarrow \theta=\cos^{-1}(0.5)=\frac{\pi}{3}\:rad
Substitute the above necessary values in the formula of each time domain specification and simplify in order to get the values of time domain specifications for given transfer function.
The following table shows the formulae of time domain specifications, substitution of necessary values and the final values.
| Time domain specification | Formula | Substitution of values in Formula | Final value |
|---|---|---|---|
Delay time |
t_d=\frac{1+0.7\delta}{\omega_n} |
t_d=\frac{1+0.7(0.5)}{2} |
t_d=0.675 sec |
Rise time |
t_r=\frac{\pi-\theta}{\omega_d} |
t_r=\frac{\pi-(\frac{\pi}{3})}{1.732} |
t_r=1.207 sec |
Peak time |
t_p=\frac{\pi}{\omega_d} |
t_p=\frac{\pi}{1.732} |
t_p=1.813 sec |
% Peak overshoot |
\%M_p=\left( e^{-\left (\frac{\delta\pi}{\sqrt{1-\delta^2}} \right ) }\right )\times 100\% |
\%M_p=\left( e^{-\left (\frac{0.5\pi}{\sqrt{1-(0.5)^2}} \right ) }\right )\times 100\% |
\% \: M_p=16.32% |
Settling time for 2% tolerance band |
t_s=\frac{4}{\delta\omega_n} |
t_S=\frac{4}{(0.5)(2)} |
t_s=4 sec |