CAPITULO IV-CONTROL I Analisis de La Respuesta Transitoria y Estacionaria
CAPITULO IV-CONTROL I Analisis de La Respuesta Transitoria y Estacionaria
CAPITULO IV-CONTROL I Analisis de La Respuesta Transitoria y Estacionaria
CAPITULO IV
ANALISIS DE LA RESPUESTA TRANSITORIA Y ESTACIONARIA
4.1. INTRODUCCIÓN:
Primer paso para analizar un sistema de control era obtener un modelo
matemático del mismo.
Existen varios métodos para el análisis del desempeño del sistema.
En la práctica la señal de entrada para un sistema de control no se
conoce con anticipación, pero es de naturaleza aleatoria, y la entrada
instantánea no puede expresarse en forma analítica.
Sólo en algunos casos especiales se conoce con anticipación la señal de
entrada y se puede expresar en forma analítica o mediante curvas; tal es
el caso del control automático de herramientas de corte.
En el análisis y diseño de sistemas de control, debemos tener una base
de comparación del desempeño de diversos sistemas de control.
Esta base se configura especificando las señales de entrada de prueba
particulares y comparando las respuestas de varios sistemas a estas
señales de entrada.
Muchos criterios de diseño se basan en tales señales o en la respuesta
del sistema a los cambios en las condiciones iniciales (sin señales de
prueba).
El uso de señales de prueba se justifica porque existe una correlación
entre las características de respuesta de un sistema para una señal de
entrada de prueba común y la capacidad del sistema de manejar las
señales de entrada reales.
4.2. SEÑALES DE PRUEBA TÍPICAS: Las señales de prueba que se usan
regularmente son las funciones:
Escalón
Rampa
Parábola
Impulso
Senoidales, etc.
Con estas señales de prueba, es posible realizar con facilidad análisis de
sistemas de control:
Matemáticos y Experimentales
Ya que las señales son funciones del tiempo muy simples.
R(s) 1 C(s)
Ts + 1
La relación entrada-salida se obtiene mediante:
C (s) 1
=
R ( s ) Ts + 1
En lo sucesivo:
Analizaremos las respuestas del sistema a entradas tales como la función:
“escalón unitario, rampa unitaria e impulso unitario”. Se supone que las
condiciones iniciales son cero.
Observe que todos los sistemas que tienen la misma función de
transferencia exhibirán la misma salida en respuesta a la misma entrada.
Para cualquier sistema físico dado, la respuesta matemática recibe una
interpretación física.
4.5.1 RESPUESTA ESCALÓN UNITARIO DE SISTEMAS DE PRIMER ORDEN:
Como la transformada de Laplace de la función escalón unitario
es 1/s, sustituyendo R(s) = 1/s en la ecuación anterior,
obtenemos
1 1
C (s) =
Ts + 1 s
Si se desarrolla C(s) en fracciones parciales simples se obtiene:
1 T 1 1
C (s) = - = -
s Ts + 1 s s + 1
T
Si se toma la transformada inversa de Laplace de la ecuación
anterior, se obtiene:
t
-
c (t ) = 1 - e T , para t �0
La ecuación anterior plantea que la salida c(t) es inicialmente
cero y al final se vuelve unitaria.
Una característica importante de tal curva de respuesta
exponencial c(t) es que:
Para t = T, el valor de c(t) es 0.632, o
Que la respuesta c(t) alcanzó 63.2% de su cambio total.
Esto se aprecia con facilidad sustituyendo t = T en c(t). Es decir:
c(T ) = 1 - e -1 = 0, 632
Observe que, conforme más pequeña es la constante de tiempo
T, más rápida es la respuesta del sistema.
R(s) + k C(s)
- s( Js + B )
C (s) k
= 2
t ( s) Js + Bs + k
Que puede escribirse como:
k
C ( s) J
=
t ( s) B B k B B k
[s + + ( ) 2 - ][ s + - ( )2 - ]
2J 2J J 2J 2J J
Los polos en lazo cerrado son complejos si:
B 2 - 4 Jk p 0 , y son reales si B 2 - 4 Jk �0
En el análisis de la respuesta transitoria, es conveniente escribir:
k B
= wn2 , = 2xwn = 2s
J J
Donde:
σ: Atenuación
wn: frecuencia natural no amortiguada
ζ: factor de amortiguamiento relativo del sistema.
B B
Bc = 2 Jk x= =
Bc 2 Jk
B: amortiguamiento real
Bc: amortiguamiento critico.
s + xwn
l -1 [ ] = e -xwnt coswd t ;
( s + xwn ) + wd
2 2
wd
l -1 [ ] = e -xwnt senwd t
( s + xwn ) + wd
2 2
1 x 2 -1)wn t 1 x 2 -1)wn t
c(t ) = 1 + e - (x + - e - (x -
2 x - 1(x + x - 1)
2 2
2 x - 1(x - x - 1)
2 2
wn e - s1t e - s2t
c (t ) = 1 + - ( ) , para t≥0
2 x 2 - 1 s1 s2
Donde:
s1 = (x + x 2 - 1)wn ; s1 = (x - x 2 - 1)wn
Por tanto, la respuesta c(t) incluye dos términos exponenciales que
decaen.
La figura siguiente contiene una familia de curvas c(t) con diversos
valores de ζ, en donde la abscisa es la variable adimensional w nt. Las
curvas solo son funciones de ζ y se obtienen
tiempo del sistema de control. Los objetivos del diseño del sistema en
cuestión determinan cuál criterio de error en porcentaje usar.
Las especificaciones en el dominio del tiempo que se proporcionaron son
muy importantes, dado que casi todos los sistemas de control son
sistemas en el dominio del tiempo; es decir, deben presentar respuestas
de tiempo aceptables. (Esto significa que el sistema de control debe
modificarse hasta que la respuesta transitoria sea satisfactoria.) Observe
que, si especificamos los valores de td, t,, tp, ts y Mp, la forma de la curva
de respuesta queda prácticamente determinada. Esto se aprecia con
claridad en la figura.
1- x 2 w
tanwd tr = - =- d
x s
Por tanto, el tiempo de subida tr es:
1 w p -b
tr = tan -1 ( d ) =
wd -s wd
Donde β se define en la fig. sgte. Es evidente que para un valor pequeño de
tr, wd debe ser grande:
senwd t p = 0
O bien:
wd t p = 0, p , 2p ,3p ,.....
Como el tiempo pico corresponde al primer pico sobreelongación máximo,
wd t p = p . Por tanto,
p
tp =
wd
El tiempo pico tp corresponde a medio ciclo de la frecuencia de oscilación
amortiguada.
Sobreelongación máxima Mp: la sobreelongación máxima se presenta en el
p
tiempo pico o en t = t p = . Por tanto, a partir de la Ecuación, Mp se
wd
obtiene como:
M p = c (t p ) - 1
c(t p ) - c(�)
Mp =
c(�)
1-x 2 )p
= e - ( s / wd ) p = e - ( x /
El porcentaje de sobrelongación máxima es e - (s / wd )p x100%
Si el valor final c(∞) de la salida no es la unidad, entonces se necesita utilizar
la ecuación siguiente:
c(t p ) - c(�)
Mp =
c(�)
Tiempo de asentamiento ts: para un sistema subamortiguado de segundo
orden, la respuesta transitoria se obtiene a partir de la ec.
e -xwnt 1- x 2
c(t ) = 1 - sen(wd t + tan -1 ), para t �0
1- x 2 x
Las curvas 1 �(e -xwnt / 1 - x 2 ) son las curvas envolventes de la respuesta
transitoria para una entrada escalón unitario. La curva de respuesta c(t)
Par de curvas envolventes para la curva de respuesta a escalón unitario del sistema.
La constante de tiempo de estas curvas envolventes es 1/ xwn
La velocidad de decaimiento de la respuesta transitoria depende del valor
de la constante de tiempo 1/ xwn . Para wn un determinado, el tiempo de
asentamiento ts es una función del factor de amortiguamiento relativo x .
p - b 3.14 - b
tr = =
wd 4
Donde β se obtiene mediante:
wd 4 p
b = tan -1 = tan -1 = 53.13�= 53.13 x = 0,93 rad
s 3 180
Por lo tanto el tiempo de subida ts es:
3.14 - 0.93
tr = = 0,55 seg
4
Tiempo pico tp: es:
p 3.14
tp = = = 0, 785 seg
wd 4
Sobre elongación máxima Mp: es:
s 3
-( )p - ( )3.14
wd
Mp =e =e 4
= 0, 095
Por tanto, el porcentaje de sobre elongación máxima es 9.5%.
Tiempo de asentamiento ts: es:
Para el criterio de 2%:
4 4
ts = = = 1.33 seg
s 3
Para el criterio de 5%:
3 3
ts = = = 1 seg
s 3
Ejemplo:
Para el sistema de la figura.
R(s) + + k 1 C(s)
- - Js + B s
kh
R(s) + k C(s)
- s ( Js + B + kkh )
p -b
tr =
wd
Donde:
wd
b = tan -1 = tan -1 1.95 = 1.10
s
Por tanto, tr es:
tr = 0.65 seg .
Tiempo de asentamiento ts: Para el criterio del 2%
4
ts = = 2.48 seg.
s
Para el criterio del 5%
3
ts = = 1.86 seg.
s
Análisis de la respuesta transitoria con MATLAB
EJEMPLO 3
Considere el sistema siguiente:
MATLAB Programa 1:
>> A=[-1 -1;6.5 0];
>> B=[1 1;1 0];
>> C=[1 0;0 1];
>> D=[0 0;0 0];
>> step(A,B,C,D)
MATLAB Programa 2
>> % ***** En este programa se representan las Curvas de
>> % respuesta escalón de un sistema que tiene dos entradas
>> %(u1 y u2) y dos salidas (y1 y y2) *****
>> % ***** Se representa en primer lugar la curva de respuesta
>> % escalón cuando la entrada es u1. A continuación se
>> % representa la curva de respuesta escalón cuando la entrada
>> % es u2 *****
>> % ***** Introducir las matrices A, B, C, y D *****
>> A=[-1 -1;6.5 0];
>> B=[1 1;1 0];
>> C=[1 0;0 1];
>> D=[0 0;0 0];
>> D=[0 0;0 0];
>> % ***** Para representar la curva de respuesta escalón
>> % cuando la entrada es u1, introduzca la orden
>> % 'step(A,B,C,D,1)' *****
>> step(A,B,C,D,1)
>> grid
>> title ('Respuesta a un escalón: Entrada=u1 (u2=0)')
>> text(3.4, -0.06,'Y1')
>> text(3.4, 1.4,'Y2')
>> step(A,B,C,D,2)
>> grid
>> title ('Respuesta a un escalón: Entrada=u2 (u1=0)')
La evaluación de las b continúa hasta que todas las restantes son cero.
Se sigue el mismo patrón de multiplicación cruzadas de los coeficientes
de las dos filas anteriores al evaluar las c, las d, las e, etc. Es decir,
b1a3 - a1b2 b1a5 - a1b3 b1a7 - a1b4
c1 = , c2 = , c3 =
b1 b1 b1
c b -bc c b -b c
d1 = 1 2 1 2 , d2 = 1 3 1 3
c1 c1
4.12. EJEMPLO:
Se va aplicar el criterio de estabilidad de Routh al siguiente polinomio de tercer
orden:
a0 s 3 + a1 s 2 + a2 s + a3 = 0
Donde todos los coeficientes son números positivos. El arreglo de coeficientes se
convierte en:
s3 a0 a2
s2 a1 a3
a1a2 - a0 a3
s1
a1
s0 a3
La condición de que todas las raíces tengan partes iguales negativas se obtiene
mediante
a1a2 �
a0 a3
EJEMPLO:
Se va aplicar el criterio de
estabilidad de Routh a la ecuación
de sexto orden:
Donde todos los coeficientes son
números positivos. El arreglo de
coeficientes se convierte en:
EJEMPLO:
Determinar el rango de estabilidad para el valor de un parámetro.
s 4 + 3s 3 + 3s 2 + 2s + k = 0
El arreglo de coeficientes se convierte en:
s4 1 3 k
s3 3 2 0
7
s2 k
3
9
s1 2- k
7
s0 k
Para la estabilidad, k debe ser positivo, y todos los coeficientes de la primera
columna deben serlo también. Por tanto,
9 9 14
2- k =0 , k =2, k=
7 7 9
14
��
k 0
9
14
Cuando: k=
9
El sistema se vuelve oscilatorio y, matemáticamente, la oscilación se mantiene
en una amplitud constante.
Obsérvese que los rangos de los parámetros de diseño que conduzcan a la
estabilidad se pueden determinar utilizando el criterio de estabilidad de
Routh Hurwizt.
EJEMPLO:
Consideremos el sistema de la Figura.
s +1
G(s) =
Con s ( s - 1)( s + 6)
s3 1 k -6
s2 5 k
4k - 30
s1
5
s0 5
Para que el sistema permanezca estable, necesitamos que
4k - 30
f0 y kf0
5
30
4k - 30 = 0 , k= = 7,5
4
Porque si K > 7.5 el sistema permanecerá estable
Figura (a) Diagrama de bloques de un controlador de encendido y apagado; (b) diagrama de bloques
de un controlador de encendido y apagado con una brecha diferencial.
U ( s ) ki
=
E ( s) s
Si se duplica el valor de e(t), el valor de u(t) varía dos veces más rápido.
Para un error de cero, el valor de u(t) permanece estacionario. En
ocasiones, la acción de control integral se denomina control de reajuste
(reset). La figura muestra un diagrama de bloques de tal controlador.
Figura (a) Diagrama de bloques de un controlador proporcional-derivativo; (b) y (c) diagramas que
muestran una entrada rampa unitaria y la salida del controlador.
TAREA
Efectos de las Acciones de Control Integral y Derivativo sobre el Desempeño de
un Sistema:
Acción de control integral.
Control integral de los sistemas de control del nivel de liquido
Respuesta a perturbaciones de par (control proporcional).
Respuesta a perturbaciones de par (control proporcional-integral).
Acción de control derivativa.
Control proporcional de sistemas con carga de inercia.