E R - Y-Compensadores
E R - Y-Compensadores
E R - Y-Compensadores
MATAMOROS
UNIDAD III
ESTABILIDAD RELATIVA Y
DISEÑO DE COMPENSADORES
Mayo 2012
ESTABILIDAD RELATIVA:
En el diseño de sistemas de control a menudo se
necesita ir más allá de la cuestión de estabilidad en
lazo cerrado. En particular, es deseable obtener
alguna medida cuantitativa de cuan lejos de ser
inestable está un lazo nominal; es decir, cuantificar
la estabilidad relativa del lazo.
Esta cuantificación puede lograrse definiendo medidas
que describan la distancia de la respuesta en
frecuencia nominal al punto de estabilidad crítica (-
1,0).
MÁRGENES DE GANANCIA Y FASE
El margen de ganancia (Mg) marca la ganancia adicional que
llevaría el lazo cerrado a la condición de estabilidad crítica.
Es el inverso de la magnitud de G(j) en la frecuencia
donde la fase vale =-180
El margen de fase (Mf) cuantifica el retardo de fase puro que
debería agregarse para alcanzar la misma condición de
estabilidad crítica. Es la cantidad de atraso en la frecuencia
de cruce G(j)=1 necesaria para llevar al sistema al limite
de la estabilidad
DETERMINACIÓN DE LOS MÁRGENES DE
GANANCIA Y FASE CON MATLAB
El comando “margin” calcula el margen de ganancia (Mg), el
margen de fase (Mf) y las correspondientes frecuencias de
cruce (Wcg, Wcp)
Cuando se introduce el comando “margin” a la computadora
(sin argumentos en el lado izquierdo), Matlab produce las
representaciones de Bode con los márgenes de ganancia y de
fase marcados con líneas verticales. Los argumentos
requeridos para el comando “margin” son el numerador y
denominador de la función de transferencia, es decir en la
forma de
margin(num, den)
EJEMPLO
Para un sistema cuya función de transferencia a lazo abierto es
G(s)H(s), indique si el sistema a lazo cerrado es estable y
cuales son Mf y Mg.
CASO 1 600( s 2)
G ( s) H ( s)
s ( s 2 17 s 70)
Código fuente:
num=[600 1200];
den=[1 17 70 0];
margin(num,den), grid
El resultado obtenido se muestra en la siguiente figura donde se
observa que Mg= y Mf=35.5. Se concluye que el sistema
es estable ya que ambos resultados márgenes son positivos.
CASO 2: 3s 1
G (s) H ( s)
s (5s 3 3s 2 4s 2)
Código fuente:
num=[3 1];
den=[5 3 4 2 0];
margin(num,den), grid
De acuerdo al resultado
se observa que Mg=-19.5
y Mf=-71.9
Por tanto se concluye que el sistema es inestable porque los
márgenes son negativos
CASO 3: Eliminando el polo en el origen del sistema anterior
el problema queda como:
3s 1
G (s) H (s) 3
5s 3s 2 4s 2
Código fuente:
num=[3 1];
den=[5 3 4 2];
margin(num,den), grid
En la figura se observa
que Mg= y Mf=17
Por tanto el sistema es
estable dado que los
márgenes son positivos
COMENTARIOS SOBRE MÁRGENES DE
GANANCIA Y FASE
Los márgenes de ganancia y fase de un sistema de control
son una medida de la proximidad de la traza polar al punto -1
+ j0. Por tanto pueden usarse como criterios de diseño.
Para determinar la estabilidad relativa deben considerarse
ambos.
Para un sistema de fase mínima ambos márgenes deben ser
positivos para que el sistema sea estable. Los márgenes
negativos indican inestabilidad.
Para obtener un desempeño satisfactorio del sistema el
margen de fase debe estar entre 30 y 60 y el de ganancia
debe ser mayor a 6dB.
PICO DE RESONANCIA (MR)Y
FRECUENCIA DE RESONANCIA R
La magnitud del pico de resonancia (Mr) proporciona un
indicio de la estabilidad relativa de un sistema.
Una magnitud de pico de resonancia grande indica la
presencia un par de polos dominantes en lazo cerrado con un
pequeño lo que produce una respuesta transitoria
inconveniente.
Una magnitud de pico de resonancia pequeño indica la
ausencia de un par de polos dominantes en lazo cerrado con
un pequeño, lo que indica un sistema bien amortiguado.
r es real si 0.7071. Por tanto no hay resonancia en lazo
cerrado si 0.7071
DETERMINACION DEL PICO Y FRECUENCIA DE
RESONANCIA CON MATLAB
Las órdenes de Matlab que se pueden utilizar para obtener el
pico de resonancia y la frecuencia de resonancia se muestran
en el siguiente ejemplo. En este caso se construye el
Diagrama de Bode para un sistema en lazo cerrado con una
ganancia de 15 para el controlador:
w = logspace(-1,1); C ( s) 15
R
h = tf([15],[900 420 43 13]);( s ) 900 s 3
420 s 2
43s 13
bode(h,w)
[mag,phase,w] = bode(h,w); [Mp,k] = max(mag);
pico_resonante = Mp
frecuencia_resonante = w(k)
El resultado que despliega el Matlab primero es la traza de
Bode
Sistema 1
Sistema 2
a) El sistema 1 tiene un ancho de banda de 1rad/seg y el
sistema 2 tiene un ancho de banda de 0.333rad/seg
b) El sistema 1 tiene un ancho de banda mayor al sistema 2
por lo que su tiempo de levantamiento es más pequeño y
una respuesta más rápida eso lo vemos comparando sus
respuestas al escalón.
Sistema 1
tr=2.3-0.1=2.2seg
Sistema 2
tr=6.92-0.32=6.6seg
con Kc=1
jT +1
G c (j ) = K c
jT +1
1
1-
sen( m ) = 2 ,
1 1+
2
Además:
1
m =
T
en donde G1(s)=KG(s)
El problema es determinar el valor de K que satisfaga el
requerimiento sobre la corriente estática de error
determinada.
2. Usando la ganancia K determinada se dibujan las trazas de
Bode de G1(j), el sistema con la ganancia ajustada pero sin
compensar. Calcular el valor del margen de fase.
3. Se determina el ángulo de adelanto de fase necesario para
agregar al sistema.
4. Determinar el factor de atenuación . Se establece la
frecuencia a la cual la magnitud del sistema no compensado
G1(j) es igual -20log(1/√). Se selecciona esta como la
nueva frecuencia de cruce máxima. Esta frecuencia
corresponde a m=1/(√T), y el cambio de fase máximo m
ocurre en ella.
5. Se determinan las frecuencias de esquina del compensador
en adelanto de modo siguiente:
Cero del compensador en adelanto: =1/T
Polo del compensador en adelanto: =1/T
6. Usando el valor de K del paso 1 y de establecido en el paso
4, se calcula la constante Kc a partir de:
Kc=K/
7. Se verifica el margen de ganancia para asegurar de que es
satisfactorio. De no ser así se repite el proceso de diseño
modificando la ubicación de polos y ceros del compensador
hasta obtener un resultado satisfactorio.
EJEMPLO:
Considere un sistema cuya función de transferencia de lazo
abierto es: 4
G(s)
s ( s 2)
Diseñar una compensación que proporcione un margen de fase
mayor que 50, margen de ganancia mayor que 10 dB, y
constante de velocidad Kv=20 s-¹.
SOLUCIÓN
Se usará un compensador en adelanto de la forma:
1
s
Ts +1 T
Gc (s) = K c Kc
Ts +1 s
1
T
K 10
Kc 41.67
0.24
7. El compensador queda diseñado con la siguiente función de
transferencia:
s 4.34
Gc ( s ) 41.67
s 18.11
El sistema compensado queda como:
s 4.34 4 166.68s 723.4
Gc ( s )G ( s ) 41.67
s 18 .11 s ( s 2) s ( s 18.11)( s 2)
Verifiquemos la respuesta en frecuencia:
num=[166.68 723.4];
den=[1 20.11 36.22 0];
margin(num,den), grid
Se observa que la
respuesta del sistema
compensado mejora
notablemente el
amortiguamiento y el
tiempo de establecimiento
COMPENSADOR EN ADELANTO CON
AMPLIFICADORES OPERACIONALES
El circuito
de un compensador en adelanto con A.O. es:
Análisis:
Para el primer A.O. tenemos: (inversor)
Donde:
Entonces:
El segundo A.O. es un inversor tal que:
Finalmente:
Por tanto:
1
s
Ts +1 T 1
G c (s) = K c Kc
Ts +1 s
1
T
Definimos Kc=K
K
También definimos: G1 ( s ) KG( s )
s ( s 1)(0.5s 1)
Primero ajustemos K para cumplir con la constante de error
estático de velocidad requerida.
Ts 1 sK
K v lim sGc ( s)G ( s ) lim s G1 ( s) lim K 5
s 0 s 0 Ts 1 s 0 s ( s 1)( 0.5s 1)
Se obtiene:
Mg=-4.44dB
Mf=-13
El sistema es inestable.
Esto se comprueba con la respuesta al escalón del sistema en
lazo cerrado.
>> num=[0 5];
>> den=[0.5 1.5 1 5];
>> step(num,den), grid
s 0.1
( s ) 0.53
Gccomo:
El compensador queda
s 0.01
s 0.1 1
El sistemaGc
compensado
( s )G ( s ) 0.es:
53
s 0.01 s s 1 0.5s 1
De donde:
COMPENSADOR EN ATRASO-ADELANTO
Un compensador en atraso-adelanto combina los efectos de
un compensador en adelanto con los de un compensador en
atraso. El resultado es un sistema con una mejora en la
respuesta transitoria, estabilidad y error de estado
estacionario. Para implementar un compensador en atraso-
adelanto, primero se diseña el compensador en adelanto para
lograr la respuesta transitoria y estabilidad deseadas, y
entonces se agrega un compensador en atraso para mejorar la
respuesta de estado estacionario.
CARACTERÍSTICAS DE UN COMPENSADOR
ATRASO-ADELANTO
Consideremos un compensador atraso-adelanto con la
siguiente función de transferencia:
1 1
s+ s+
T1 T2
Gc ( s ) K c
1
s T s + T
1 2
1- 1 1
sen( m ) = sen( m )
1+ 1
Tomando m=55 tenemos
1
senm senm senm 1 1 senm 1 senm
1
1 senm 1 sen (55)
10
1 senm 1 sen (55)
La frecuencia de esquina =1/T2=0.0142rad/s (polo de la
parte de atraso)
Para la parte de adelanto: Como la nueva frecuencia de cruce
ganancia es =1.42rad/s, en esa frecuencia la ganancia es
10.6 dB.