Inicie su sesión para modificar los apuntes, usar los foros, etc. Si no es usuario de la web, regístrese.
Se desea diseñar el sistema de control de la figura. Teniendo en cuenta las características propias de la aplicación, se ha optado por diseñar un sistema de control exento de error de posición y con un error de velocidad menor del 10%. El tiempo de asentamiento deberá ser inferior a 1 segundo, y el sobreimpulso máximo admitido es del 5%.
$C(s)=K_1$, $G(s)=\frac{1}{(s^2+3s)}$ y $H(s)=K_2s+1$
Resolución:
Los requerimientos del diseño son:
Error de posición --> ess (pos)= 0
Error de velocidad --> ess (vel)< 10%
Tiempo de establecimiento --> ts (2%)< 1 seg
Sobreimpulso --> Mp< 5%
En primer lugar vemos de que tipo es el sistema calculando su FT de lazo abierto, que es:
$\ FT = C(s) \cdot G(s) \cdot H(s) = K_1 \cdot \frac{1}{{(s^2 + 3s)}} \cdot (K_1 s + 1) = \frac{{K_1 K_2 s + K_1 }}{{s(s + 3)}}
\$Vemos que es de tipo 1, por lo que el eSS (pos) será cero, por lo que no hace falta calcularlo.
Para calcular el ess(vel) lo haremos de dos modos: uno usando la expresión del error de realimentación, y otro usando el error real.
Modo 1--> El error de realimentación a partir del diagrama de bloques es:
$\ \left. \begin{array}{l}
\mu (s) = Y(s) \cdot H(s) \\
E(s) = R(s) - \mu (s) \\
Y(s) = E(s) \cdot C(s) \cdot G(s) \\
\end{array} \right\}E(s) = R(s) - E(s) \cdot C(s) \cdot G(s) \cdot H(s){\rm{ }} \Rightarrow
\$
$\ \Rightarrow {\rm{ }}E(s) = R(s)\frac{1}{{1 + C(s) \cdot G(s) \cdot H(s)}} = R(s)\frac{{s(s + 3)}}{{s(s + 3) + K_1 (1 + K_2 s)}}
\$
Ahora aplicamos el teorema del valor final con una señal de entrada rampa R(s)=1/s2:
$\ e_{ss} (vel) = \ {\lim }\limits_{x \to 0} s \cdot E(s) = \ {\lim }\limits_{x \to 0} s \cdot R(s)\frac{{s(s + 3)}}{{s(s + 3) + K_1 (1 + K_2 s)}} =
\$
$\ = \ {\lim }\limits_{x \to 0} s \cdot \frac{1}{{s^2 }} \cdot \frac{{s(s + 3)}}{{s(s + 3) + K_1 (1 + K_2 s)}} = \frac{3}{{K_1 (1 + 0)}} = \frac{3}{{K_1 }}
\$
Sabemos que uno de los requisitos que tiene que cumplir es que ess (vel)< 10%, por tanto:
$\ \frac{3}{{K_1 }} < 0,1{\rm{ }} \Rightarrow {\rm{ }}K_1 > \frac{3}{{0,1}}{\rm{ }} \Rightarrow {\rm{ }}K_1 > 30
\$
Calculamos la FT de lazo cerrado para ver qué orden tiene el sistema:
$\ FTLC = \frac{{Y(s)}}{{R(s)}} = \frac{{C(s) \cdot G(s)}}{{1 + C(s) \cdot G(s) \cdot H(s)}} = \frac{{K_1 \frac{1}{{(s^2 + 3s)}}}}{{1 + K_1 \frac{1}{{(s^2 + 3s)}}(K_2 s + 1)}} =
\$
$\ = \frac{{K_1 }}{{s(s + 3) + K{}_1(1 + K_2 s)}} = \frac{{K_1 }}{{s^2 + (3 + K{}_1K_2 )s + K_1 }}
\$
Comprobamos que el sistema es de 2º orden, por lo que el sobreimpulso será:
$\ Mp = e^{\frac{{ - \delta \pi }}{{\sqrt {1 - \delta ^2 } }}}
\$
De esta expresión despejaremos el coeficiente de amortiguamiento δ, quedando:
$\ \delta = \sqrt {\frac{{\ln ^2 Mp}}{{\pi + \ln ^2 Mp}}} = \sqrt {\frac{{\ln ^2 0,05}}{{\pi + \ln ^2 0,05}}} = 0,69
\$
Como Mp debe ser menor del 5%, observamos como cambian los polos al variar δ en los dos casos, que son:
$\ p_{1,2} = - \delta \omega _n \pm j\omega _n \sqrt {1 - \delta ^2 } \left\{ \begin{array}{l}
\delta {\rm{ }}tiende{\rm{ }}a{\rm{ }}0 \Rightarrow p_{1,2} = \pm j\omega _n \sqrt {1 - \delta ^2 } \Rightarrow Oscila{\rm{ }}y{\rm{ }}aumenta{\rm{ }}Mp \\
\delta {\rm{ }}tiende{\rm{ }}a{\rm{ 1}} \Rightarrow p_{1,2} = - \omega _n \Rightarrow Mp{\rm{ }}tiende{\rm{ }}a{\rm{ }}0 \\
\end{array} \right.
\$
Por tanto, δ>0,69.
El tiempo de establecimiento en los sistemas de 2º orden para un rango del 2% es:
$\ t_s (2\% ) = \frac{4}{{\delta \omega _n }} < 1{\rm{ seg}}
\$
De donde podemos despejar la frecuencia natural no amortiguada ωn:
$\ \omega _n > \frac{4}{{\delta \cdot t_s }} = \frac{4}{{0,69 \cdot 1}} = 5,8{\rm{ }}rad/s
\$
Despreciando $\ \sqrt {1 - \delta ^2 }
\$ esto sería el valor absoluto de la parte imaginaria de los polos.
$\ \cos \beta = \frac{{\delta \omega _n }}{{\omega _n }} = \delta = 0,69{\rm{ }} \Rightarrow {\rm{ }}\beta {\rm{ = arcos 0}}{\rm{,69 = 46}}{\rm{,37}}^ \circ
\$
Como δ>0,69 esta aumenta , lo que hace que p1,2 tiendan a ser más reales, y por tanto β menor. Así que β < 46,37º.
Obtenemos la parte real de los polos a partir de la expresión de ts:
$\ \frac{4}{{\delta \omega _n }} < 1{\rm{ }} \Rightarrow {\rm{ }}\delta \omega _n > 4{\rm{ }} \Rightarrow {\rm{ - }}\delta \omega _n < - 4
\$
Por último, obtendremos la expresión para K2 comparando la FTLC con la expresión modelo de un sistema de 2º orden, debiéndose cumplir las inecuaciones anteriores.
$\ \frac{{K_1 }}{{s^2 + (3 + K{}_1K_2 )s + K_1 }} = \frac{{K_0 \omega _n^2 }}{{s^2 + 2\delta \omega _n s + \omega _n^2 }}
\$
La frecuencia natural no amortiguada será:
$\ \omega _n^2 = K_1 {\rm{ }}{\rm{, }}\omega _n^{} = \sqrt {K_1 } \Rightarrow {\rm{ }}\omega _n^{} > \sqrt {30} = 5,48{\rm{ }}rad/s
\$
Al ser este valor menor que antes es menos restrictivo, por lo que nos quedamos con el anterior para el diseño, siendo el requisito:
$\ K_1 {\rm{ = }}\omega _n^2 {\rm{ }} = 5,8^2 {\rm{ = 33}}{\rm{,64 }} \Rightarrow {\rm{ }}K_1 > 33,64
\$
Y la expresión para el coef. de amortiguamiento δ será:
$\ \delta = \frac{{3 + K_1 K_2 }}{{2\omega _n }}
\$
A continuación escogeremos valores para K1 y K2. Para ello tendremos en cuenta que se cumplan los requisitos calculados K1>33,64 y δ>0,69.
Obtenemos un valor de K2 sustituyendo estos valores en la última expresión que hemos obtenido y despejando:
$\ K_2 = \frac{{2\delta \sqrt {K_1 } - 3}}{{K_1 }} = \frac{{2 \cdot 0,69\sqrt {33,64} - 3}}{{33,64}} = 0,148{\rm{ }} \Rightarrow {\rm{ }}K_2 > 0,148
\$
Con esto, y dando un valor a K1 y a K2 que no estén muy lejos de los límites calculados pero tampoco demasiado cerca (p.ej. K1=35 y K2=0,18) obtenemos δ, que será:
$\ \delta = \frac{{3 + K_1 K_2 }}{{2\omega _n }} = \frac{{3 + 35 \cdot 0,18}}{{2\sqrt {35} }} = 0,78
\$
Comprobamos que cumple el requisito de diseño.
Y el sobreimpulso será:
$\ Mp = e^{\frac{{ - \delta \pi }}{{\sqrt {1 - \delta ^2 } }}} = e^{\frac{{ - 0,78\pi }}{{\sqrt {1 - 0,78^2 } }}} = 0,02{\rm{ }}(2\% )
\$
Modo 2--> El error real a partir del diagrama de bloques es:
$\ \left. \begin{array}{l}
E'(s) = R(s) - Y(s) \\
FTLC = \frac{{Y(s)}}{{R(s)}} = \frac{{K_1 }}{{s^2 + (3 + K{}_1K_2 )s + K_1 }} \\
\end{array} \right\}E'(s) = R(s){\rm{ }}\left( {{\rm{1 - }}\frac{{K_1 }}{{s^2 + (3 + K{}_1K_2 )s + K_1 }}} \right) \Rightarrow
\$
$\ \Rightarrow {\rm{ }}E'(s) = R(s)\left( {\frac{{s(s + 3 + K_1 K_2 )}}{{s^2 + s(3 + K_1 K_2 ) + K_1 }}} \right)
\$