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 (exento de error de posición ).
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_2 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.
$\ \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
\$
Puesto que:
$\ \cos \beta = \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
\$
Los polos de lazo cerrado en el plano S se dibujarán como:
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\% )
\$
Nota: Este apartado no era necesario calcularlo en el examen.
$\ \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)
\$
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 + K_1 K_2 )}}{{s^2 + s(3 + K_1 K_2 ) + K_1 }} =
\$
$\ { = \lim }\limits_{x \to 0} s \cdot \frac{1}{{s^2 }} \cdot \frac{{s(s + 3 + K_1 K_2 )}}{{s^2 + s(3 + K_1 K_2 ) + K_1 }} = \frac{{3 + K_1 K_2 }}{{K_1 }} = \frac{3}{{K_1 }} + K_2
\$
Sabemos que uno de los requisitos que tiene que cumplir es que ess (vel)< 10%, por tanto:
$\ \frac{3}{{K_1 }} + K_2 < 0,1{\rm{ }} \Rightarrow {\rm{ }}K_1 > \frac{3}{{0,1 - K_2 }}{\rm{ }}{\rm{, }}K_2 < \frac{{0,1K_1 - 3}}{{K_1 }}
\$
Al igual que antes, en principio los requisitos serán δ>0,69 (por tanto β<46,37º).
A continuación escogeremos valores para K1 y K2. Para ello tendremos en cuenta que se cumplan los requisitos calculados $\ K_1 > \frac{3}{{0,1 - K_2 }}{\rm{ }}{\rm{, }}K_2 < \frac{{0,1K_1 - 3}}{{K_1 }}{\rm{ }}y{\rm{ }}\delta > 0,69
\$.
Obtenemos un valor de K2 dando un valor a K1 adecuado (p.ej., el de antes K1=35) y sustituyendo estos valores en las expresiones que hemos obtenido antes:
$\ K_2 < \frac{{0,1K_1 - 3}}{{K_1 }} = \frac{{0,1 \cdot 35 - 3}}{{35}} = 0,0142{\rm{ }} \Rightarrow {\rm{ }}K_2 < 0,0142
\$
$\ K_2 < \frac{{2\delta \sqrt {K_1 } - 3}}{{K_1 }} = \frac{{2 \cdot 0,69\sqrt {35} - 3}}{{35}} = 0,147{\rm{ }} \Rightarrow {\rm{ }}K_2 < 0,147
\$
El primero es más restrictivo, por lo que damos un valor a K2 que no esté muy lejos de este límite pero tampoco demasiado cerca (p.ej. K2=0,01).
Comprobamos los resultados obtenidos en los dos modos de hacer el ejercicio mediante simulación en Matlab.
El código que genera todas las simulaciones es:
clc, clear all, close all
K1=35;
K2=0.18;
G=tf(1,[1,3,0]);
H=tf([K2,1],1);
S=feedback(K1*G,H); %FT= Y(s)/R(s)
%Simulacion frente a entrada ESCALON (Solo para comprobar que es casi cero)
t=0:1e-3:3;
figure
[y,t]=step(S,t);
u=ones(size(t)); %Escalon
ess_pos=u(end)-y(end) %error de posicion
plot(t,y,t,u,'r--'), xlabel('Tiempo (s)'), ylabel('Escalon')
%Simulacion frente a entrada RAMPA (unitaria)
r=t; %Funcion rampa unitaria
%Modo 1: Error de realimetacion
E_REALIM=tf([1,3,0],[1,3+K1*K2,K1]);
[e_vel,t]=lsim(E_REALIM,r,t);
figure,plot(t,e_vel), xlabel('Tiempo (s)'), ylabel('Error de realimentacion')
ess_vel_realim=e_vel(end) %error de velocidad de realimentacion
ess_vel_realim_teorico=3/K1 %error de velocidad de realimentacion teorico
%Modo 1: Error real
K1=35;
K2=0.01;
G=tf(1,[1,3,0]);
H=tf([K2,1],1);
S=feedback(K1*G,H); %FT= Y(s)/R(s)
r=t; %Funcion rampa unitaria
[y2,t]=lsim(S,r,t);
ess_vel_real=r(end)-y2(end) %error de velocidad real
ess_vel_real_teorico=(3/K1)+K2 %error de velocidad real teorico
figure,plot(t,y2,t,r,'r--'), xlabel('Tiempo (s)'), ylabel('Rampa')
Como puede observarse, tanto en los resultados de la consola como en las representaciones gráficas, se cumplen las condiciones del diseño. Las pequeñas diferencias existentes entre los errores teóricos y reales (según la simulación) son producidos porque el tiempo de simulación evidentemente no es infinito. Conforme se aumenta el tiempo de simulación, los errores teóricos y reales (según la simulación) tienden a converger al mismo valor.