Aug
15
Response of Second Order System In this chapter, let us discuss the time response of second order system. Consider the following block diagram of closed loop control system. Here, an open loop transfer function, $frac{omega ^2_n}{s(s+2delta omega_n)}$ is connected with a unity negative feedback. We know that the transfer function of the closed loop control system having unity negative feedback as $$frac{C(s)}{R(s)}=frac{G(s)}{1+G(s)}$$ Substitute, $G(s)=frac{omega ^2_n}{s(s+2delta omega_n)}$ in the above equation. $$frac{C(s)}{R(s)}=frac{left (frac{omega ^2_n}{s(s+2delta omega_n)} right )}{1+ left ( frac{omega ^2_n}{s(s+2delta omega_n)} right )}=frac{omega _n^2}{s^2+2delta omega _ns+omega _n^2}$$ The power of ‘s’ is two in the denominator term. Hence, the above transfer function is of the second order and the system is said to be the second order system. The characteristic equation is – $$s^2+2delta omega _ns+omega _n^2=0$$ The roots of characteristic equation are – $$s=frac{-2omega delta _npm sqrt{(2deltaomega _n)^2-4omega _n^2}}{2}=frac{-2(deltaomega _npm omega _nsqrt{delta ^2-1})}{2}$$ $$Rightarrow s=-delta omega_n pm omega _nsqrt{delta ^2-1}$$ The two roots are imaginary when δ = 0. The two roots are real and equal when δ = 1. The two roots are real but not equal when δ > 1. The two roots are complex conjugate when 0 < δ < 1. We can write $C(s)$ equation as, $$C(s)=left ( frac{omega _n^2}{s^2+2deltaomega_ns+omega_n^2} right )R(s)$$ Where, C(s) is the Laplace transform of the output signal, c(t) R(s) is the Laplace transform of the input signal, r(t) ωn is the natural frequency δ is the damping ratio. Follow these steps to get the response (output) of the second order system in the time domain. Take Laplace transform of the input signal, $r(t)$. Consider the equation, $C(s)=left ( frac{omega _n^2}{s^2+2deltaomega_ns+omega_n^2} right )R(s)$ Substitute $R(s)$ value in the above equation. Do partial fractions of $C(s)$ if required. Apply inverse Laplace transform to $C(s)$. Step Response of Second Order System Consider the unit step signal as an input to the second order system. Laplace transform of the unit step signal is, $$R(s)=frac{1}{s}$$ We know the transfer function of the second order closed loop control system is, $$frac{C(s)}{R(s)}=frac{omega _n^2}{s^2+2deltaomega_ns+omega_n^2}$$ Case 1: δ = 0 Substitute, $delta = 0$ in the transfer function. $$frac{C(s)}{R(s)}=frac{omega_n^2}{s^2+omega_n^2}$$ $$Rightarrow C(s)=left( frac{omega_n^2}{s^2+omega_n^2} right )R(s)$$ Substitute, $R(s) = frac{1}{s}$ in the above equation. $$C(s)=left( frac{omega_n^2}{s^2+omega_n^2} right )left( frac{1}{s} right )=frac{omega_n^2}{s(s^2+omega_n^2)}$$ Apply inverse Laplace transform on both the sides. $$c(t)=left ( 1-cos(omega_n t) right )u(t)$$ So, the unit step response of the second order system when $/delta = 0$ will be a continuous time signal with constant amplitude and frequency. Case 2: δ = 1 Substitute, $/delta = 1$ in the transfer function. $$frac{C(s)}{R(s)}=frac{omega_n^2}{s^2+2omega_ns+omega_n^2}$$ $$Rightarrow C(s)=left( frac{omega_n^2}{(s+omega_n)^2} right)R(s)$$ Substitute, $R(s) = frac{1}{s}$ in the above equation. $$C(s)=left( frac{omega_n^2}{(s+omega_n)^2} right)left ( frac{1}{s} right)=frac{omega_n^2}{s(s+omega_n)^2}$$ Do partial fractions of $C(s)$. $$C(s)=frac{omega_n^2}{s(s+omega_n)^2}=frac{A}{s}+frac{B}{s+omega_n}+frac{C}{(s+omega_n)^2}$$ After simplifying, you will get the values of A, B and C as $1,: -1: and : −omega _n$ respectively. Substitute these values in the above partial fraction expansion of $C(s)$. $$C(s)=frac{1}{s}-frac{1}{s+omega_n}-frac{omega_n}{(s+omega_n)^2}$$ Apply inverse Laplace transform on both the sides. $$c(t)=(1-e^{-omega_nt}-omega _nte^{-omega_nt})u(t)$$ So, the unit step response of the second order system will try to reach the step input in steady state. Case 3: 0 < δ < 1 We can modify the denominator term of the transfer function as follows − $$s^2+2deltaomega_ns+omega_n^2=left { s^2+2(s)(delta omega_n)+(delta omega_n)^2 right }+omega_n^2-(deltaomega_n)^2$$ $$=(s+deltaomega_n)^2+omega_n^2(1-delta^2)$$ The transfer function becomes, $$frac{C(s)}{R(s)}=frac{omega_n^2}{(s+deltaomega_n)^2+omega_n^2(1-delta^2)}$$ $$Rightarrow C(s)=left( frac{omega_n^2}{(s+deltaomega_n)^2+omega_n^2(1-delta^2)} right )R(s)$$ Substitute, $R(s) = frac{1}{s}$ in the above equation. $$C(s)=left( frac{omega_n^2}{(s+deltaomega_n)^2+omega_n^2(1-delta^2)} right )left( frac{1}{s} right )=frac{omega_n^2}{sleft ((s+deltaomega_n)^2+omega_n^2(1-delta^2) right)}$$ Do partial fractions of $C(s)$. $$C(s)=frac{omega_n^2}{sleft ((s+deltaomega_n)^2+omega_n^2(1-delta^2) right)}=frac{A}{s}+frac{Bs+C}{(s+deltaomega_n)^2+omega_n^2(1-delta^2)}$$ After simplifying, you will get the values of A, B and C as $1,: -1 : and : −2delta omega _n$ respectively. Substitute these values in the above partial fraction expansion of C(s). $$C(s)=frac{1}{s}-frac{s+2deltaomega_n}{(s+deltaomega_n)^2+omega_n^2(1-delta^2)}$$ $$C(s)=frac{1}{s}-frac{s+deltaomega_n}{(s+deltaomega_n)^2+omega_n^2(1-delta^2)}-frac{deltaomega_n}{(s+deltaomega_n)^2+omega_n^2(1-delta^2)}$$ $C(s)=frac{1}{s}-frac{(s+deltaomega_n)}{(s+deltaomega_n)^2+(omega_nsqrt{1-delta^2})^2}-frac{delta}{sqrt{1-delta^2}}left ( frac{omega_nsqrt{1-delta^2}}{(s+deltaomega_n)^2+(omega_nsqrt{1-delta^2})^2} right )$ Substitute, $omega_nsqrt{1-delta^2}$ as $omega_d$ in the above equation. $$C(s)=frac{1}{s}-frac{(s+deltaomega_n)}{(s+deltaomega_n)^2+omega_d^2}-frac{delta}{sqrt{1-delta^2}}left ( frac{omega_d}{(s+deltaomega_n)^2+omega_d^2} right )$$ Apply inverse Laplace transform on both the sides. $$c(t)=left ( 1-e^{-delta omega_nt}cos(omega_dt)-frac{delta}{sqrt{1-delta^2}}e^{-deltaomega_nt}sin(omega_dt) right )u(t)$$ $$c(t)=left ( 1-frac{e^{-deltaomega_nt}}{sqrt{1-delta^2}}left ( (sqrt{1-delta^2})cos(omega_dt)+delta sin(omega_dt) right ) right )u(t)$$ If $sqrt{1-delta^2}=sin(theta)$, then ‘δ’ will be cos(θ). Substitute these values in the above equation. $$c(t)=left ( 1-frac{e^{-deltaomega_nt}}{sqrt{1-delta^2}}(sin(theta)cos(omega_dt)+cos(theta)sin(omega_dt)) right )u(t)$$ $$Rightarrow c(t)=left ( 1-left ( frac{e^{-deltaomega_nt}}{sqrt{1-delta^2}} right )sin(omega_dt+theta) right )u(t)$$ So, the unit step response of the second order system is having damped oscillations (decreasing amplitude) when ‘δ’ lies between zero and one. Case 4: δ > 1 We can modify the denominator term of the transfer function as follows − $$s^2+2deltaomega_ns+omega_n^2=left { s^2+2(s)(deltaomega_n)+(deltaomega_n)^2 right }+omega_n^2-(deltaomega_n)^2$$ $$=left ( s+deltaomega_n right )^2-omega_n^2left ( delta^2-1 right )$$ The transfer function becomes, $$frac{C(s)}{R(s)}=frac{omega_n^2}{(s+deltaomega_n)^2-omega_n^2(delta^2-1)}$$ $$Rightarrow C(s)=left ( frac{omega_n^2}{(s+deltaomega_n)^2-omega_n^2(delta^2-1)} right )R(s)$$ Substitute, $R(s) = frac{1}{s}$ in the above equation. $C(s)=left ( frac{omega_n^2}{(s+deltaomega_n)^2-(omega_nsqrt{delta^2-1})^2} right )left ( frac{1}{s} right )=frac{omega_n^2}{s(s+deltaomega_n+omega_nsqrt{delta^2-1})(s+deltaomega_n-omega_nsqrt{delta^2-1})}$ Do partial fractions of $C(s)$. $$C(s)=frac{omega_n^2}{s(s+deltaomega_n+omega_nsqrt{delta^2-1})(s+deltaomega_n-omega_nsqrt{delta^2-1})}$$ $$=frac{A}{s}+frac{B}{s+deltaomega_n+omega_nsqrt{delta^2-1}}+frac{C}{s+deltaomega_n-omega_nsqrt{delta^2-1}}$$ After simplifying, you will get the values of A, B and C as 1, $frac{1}{2(delta+sqrt{delta^2-1})(sqrt{delta^2-1})}$ and $frac{-1}{2(delta-sqrt{delta^2-1})(sqrt{delta^2-1})}$ respectively. Substitute these values in above partial fraction expansion of $C(s)$. $$C(s)=frac{1}{s}+frac{1}{2(delta+sqrt{delta^2-1})(sqrt{delta^2-1})}left ( frac{1}{s+deltaomega_n+omega_nsqrt{delta^2-1}} right )-left ( frac{1}{2(delta-sqrt{delta^2-1})(sqrt{delta^2-1})} right )left ( frac{1}{s+deltaomega_n-omega_nsqrt{delta^2-1}} right )$$ Apply inverse Laplace transform on both the sides. $c(t)=left ( 1+left ( frac{1}{2(delta+sqrt{delta^2-1})(sqrt{delta^2-1})} right )e^{-(deltaomega_n+omega_nsqrt{delta^2-1})t}-left ( frac{1}{2(delta-sqrt{delta^2-1})(sqrt{delta^2-1})} right )e^{-(deltaomega_n-omega_nsqrt{delta^2-1})t} right )u(t)$ Since it is over damped, the unit step response of the second order system when δ > 1 will never reach step input in the steady state. Impulse Response of Second Order System The impulse response of the second order system can be obtained by using any one of these two methods. Follow the procedure involved while deriving step response by considering the value of $R(s)$ as 1 instead of $frac{1}{s}$. Do the differentiation of the step response. The following table shows the impulse response of the second order system for 4 cases of the damping ratio. Condition of Damping ratio Impulse response for t ≥ 0 δ = 0 $omega_nsin(omega_nt)$ δ = 1 $omega_n^2te^{-omega_nt}$ 0 < δ < 1 $left ( frac{omega_ne^{-deltaomega_nt}}{sqrt{1-delta^2}} right )sin(omega_dt)$ δ > 1 $left ( frac{omega_n}{2sqrt{delta^2-1}} right )left ( e^{-(deltaomega_n-omega_nsqrt{delta^2-1})t}-e^{-(deltaomega_n+omega_nsqrt{delta^2-1})t} right )$ Learning working make money