Ingenieria de Control
Ingenieria de Control
Ingenieria de Control
DISEO DE REGULADORES
DISEO DE REGULADORES PID MEDIANTE EL LUGAR DE LAS RAICES
El diseo de un regulador mediante tcnicas basadas en el lugar de las races est basado en la
capacidad de anlisis del diseador. En general, a partir de las especificaciones del controlador,
tpicamente error en rgimen permanente, velocidad de respuesta (tp, tl , ts) y sobre oscilaciones que
utiliza la accin proporcional (P) para, en la medida de lo posible, disminuir el error de rgimen
permanente consiguiendo una respuesta dinmica adecuada. Si esta no se ajusta a las
especificaciones, se introduce una accin derivativa (D) que mejora los comportamientos pero
aumenta el error, y finalmente, si es necesario, se incorpora una accin integradora tratando de no
afectar al comportamiento dinmico, disminuir o anular el error en forma permanente.
Condiciones dinmicas
Pd
Condiciones estticas
Error
SI
Pertenecen los Pd al
LDR?
Se verifica la
condicin de error?
NO
Necesitara accin derivativa R(s)= K(s*a)
SI
SI
Se verifica la
condicin de error?
NO
Fin
PD
Fin
PI
1+
K
Tangente en el punto
de inflexin
Ts
Figura 9.0-2
Mtodo
Mtodo Tipo de
Regulador
Kp
Ti
Td
1 Ta
K Tv
Mtodo I
Mtodo II
PI
0.9 Ta
K Tv
3.33 Tv
PID
1.2 Ta
K Tv
2 Tv
0.5 Kcr
PI
0.45 Kcr
0.85 Tcr
0.6 Kcr
0.5 Tcr
PID
0.5 Tcr
0.125 Tcr
a<b
De a a b se conseguir que aumente el margen de fase del sistema. No obstante, tambin aumentara
la ganancia, con lo que habr tendencia a que la frecuencia de cruce de ganancia wt se desestabilice
a frecuencias mayores y disminuya en cierta medida el margen de fase.
Efectos de la red de atraso
La funcin de transferencia de una red de atraso es:
( s+ a )
a> b
( s+ b )
Gc (s)=K
Su efecto sobre la respuesta en frecuencia del sistema es el de proporcionar una adecuacin de los
mdulos a partir de una determinada frecuencia. Este efecto puede ser aprovechado para modificar
la ganancia esttica del sistema y por lo tanto el error en rgimen permanente haga aumentar el
margen de ganancia y desplazar hacia frecuencias menores la frecuencia de cruce de ganancia, con
lo que habitualmente aumentar el margen de fase.
Ajuste de una red de adelanto o atraso
Segn los efectos explicados anteriormente, los ajustes de las redes de adelanto y de retraso pueden
hacerse segn los siguientes pasos.
Datos de partida,
Especificacin de error y
deseada
o.
Obtener
1x
x=sen m
1+ x
Obtener
m=+10 log
K
Leer wm en el Bode de K .G ( s ) como la frecuencia
en la que:
20 logKG ( s ) Km
Obtener
T=
1
Wm
1 s +1/T
Gc= .
s+1/ T
Dibujar el Bode de Gc.K G(s) y comprobar que el valor de
nuevo cumple la especificacin. Si no, repetir con mayor
Datos de partida,
Especificacin de error y
deseada
o.
20 log=B[dB]
Obtener
T=
10
gk
1 s +1/T
Gc= .
s+1 /T
Dibujar el Bode de Gc.K G(s) y comprobar que el valor de
nuevo cumple la especificacin
Figura 9.1-1
Se pide disear el regulador ideal ms sencillo que verifique:
Mp= 5%, ts=2.1 s, ep< 50%
Ante entradas escaln en X(s). Para el regulador diseado, comentar el efecto de la perturbacin
Z(s) en la salida Y(s).
tan
.100 5 =45
t ( s ) = 2.1 s =1.5
Esdecir:
w D =tan=1.5
PortantoPd=1.5j1.5
ReflejandoenlaFigura9.12
Im
1.5j
-1.5
Re
-2
-1
-1.5j
Figura 9.1-2
Como puede verse, el lugar de las races para por Pd. Por lo tanto, basta con el rgimen
proporcional. Para calcular el valor de la ganancia, aplicaremos el criterio del mdulo:
K=d1 d2=0.52 +1.52 =2.5
con K=K R .1 es decir K R =2.5 y R( s )=2.5
e p=
2.5
=1.25
( s+1)(s +2)
1
=0.44=4 <50
1+ K p
La especificacin se cumple y, por tanto, es vlido el regulador R(s)=2.5 La simulacin de salida del
sistema y(t) ante entrada escaln unitario en x(t) aparece reflejada en la Figura 9.1-3
0.55
Figura 9.1-3
Ante las entradas en escaln en Z(s), el rgimen transitorio es el mismo, ya que la funcin de
transferencia entre Z(s) e Y(s) tiene los mismos polos y ceros que la funcin entre X(s).
Y (s)
1
=
Z (s) ( s+1)( s+2)+K R
Sin embargo, el valor final en rgimen permanente vara:
y p=
1
=0.22
KR+ 2
En lugar de 1-0.44=0.55. Cuanto mayor sea el valor de Kr, menor es la influencia de transferencia
de y(t) en rgimen permanente y menor es el error ep de y con respecto a x(t).
Figura 9.2-1
Se pide disear el regulador ideal ms sencillo que verifique:
M p 5 , t s 2.1 s , e p 10
Ante entradas escaln en X(s). Para el regulador diseado, comentar el efecto de la perturbacin
Z(s) en la salida Y(s).
tan
.100 5 45
t ( s ) = 2.1 s 1.5
Esdecir:
w D =tan=1.5
PortantoPd=1.5j1.5sisedeseacumplirlasespecificacionesenellmitedado.
ReflejandoenlaFigura9.12
Im
1.5j
-1.5
Re
-2
-1
-1.5j
Figura 9.2-2
Como puede verse, el lugar de las races pasa por Pd. En principio bastara con un rgimen
proporcional. Para calcular el valor de la ganancia, aplicaremos el criterio del mdulo:
K=d1 d2=0.52 +1.52 =2.5
K p=lim
s0
e p=
2.5
=1.25
(s+1)(s +2)
1
=0.44=4
1+ K p
s+ a
s
Se elegir a suficientemente pequeo como para que el lugar de las races se vea apreciablemente.
En concreto, se tomar a 1/6 de la distancia de los polos dominantes en donde estar cerrada al eje
imaginario:
Con ello, el nuevo lugar de las races queda reflejado en la Figura 9.2-3
Figura 9.2-3
YKr=2.5.Elreguladorresulta:
R(s)=2.5
s +0.25
s
Sin embargo, el rgimen dinmico se ve afectado por la adicin de la parte integral, dado que el
diseo que se ha efectuado ha sido aproximado. Veamos a continuacin el efecto.
180arctan
][
1.1
1.1
1.1
1.1
arctan
+180arctan
+180arctan
=180
1.1b
0.9
0.1
2.1
Es decir:
=1.11
Y por tanto:
(s + )2 Wd 2=S 2 +2.84 s+ 4.27
Tambin se podra recalcular la ganancia aplicando el criterio del mdulo:
1.42
2
2
K=
2
1.5 +
Aunque el error que se comete al mantenerse 2.5 es muy pequeo. La funcin de transferencia
resulta:
M (s )=
2.5( s+ 0.25)
2.5( s+ 0.25)
=
s ( s+1 ) ( s +2 ) +2.5( s+0.25) ( s+0.16)( s 2+2.84 s+ 4.27)
3.9
(s +2.84 s+ 4.27)
2
Se obtendra
M ( p )=e
1.42
1.5
.100=5.1
t ( s ) = .0 .42=2.21 s
1
La respuesta del sistema realimentado aparece reflejada en la Figura 9.2-5
Figura 9.2-5
Figura 9.3-1
Se pide disear el regulador ideal ms sencillo que verifique:
Mp= 5%, ts=1 s, ep 25%
Ante entradas escaln en X(s).
t ( s ) = 1 s =3.14
Esdecir:
w D =tan=3.14
PortantoPd=3.14j3.14
El lugar de las races para el sistema anterior aparece reflejado en la Figura 9.3-2:
Im
3.14j
-1
-2
Re
-3.14j
Figura 9.3-2
El lugar de las races no pasa por Pd, por lo que no es suficiente con un regulador proporcional.
Hace falta un regulador de tipo PD:
R( s)=K R (s+b)
Teniendo en cuenta que interesa un valor de b grande (con el fin de que el cero introducido afecte
dinmicamente en cadena cerrada), el nuevo lugar de las races aparece reflejado en la Figura 9.3-3
Im
3.14j
Re
-b
-1
-3.14
-2
-3.14j
Figura 9.3-3
arctan
][
3.14
3.14
3.14
180arctan
+180arctan
=180
b3.14
1.14
2.14
Im
3.14j
Re
-b
-1
-3.14
-2
Figura 9.3-4
Es decir:
arctan
3.14
=54.2
b3.14
Y (t)
Figura 9.3-5
En cuanto al rgimen permanente, no se ha hecho an ninguna consideracin. Veamos si coinciden
las especificaciones:
K p=lim
s0
e p=
3.28 (s +5.4)
=8.59
(s +1)(s+ 2)
1
=0.10=10 <25
1+ K p
Figura 9.4-1
e p=
3.28 (s +5.4)
=8.59
(s +1)(s+ 2)
1
=0.10=10
1+ K p
R ( s )=K R
( s+a )( s+b )
s
El parmetro a se elegir a 1/6 de la distancia de los polos dominantes en cadena cerrada en el eje
imaginario:
a=
3.14
=0.52
6
Con ello, el nuevo lugar de las races queda reflejado en la Figura 9.4-2.
Figura 9.4-2
El regulador resulta:
R ( s )=3.28
e p=0
(s+0.52)(s +5.4)
s
El rgimen dinmico se ve afectado por la adicin de la parte integral, tal y como se ve en la Figura
9.4-3:
Figura 9.4-3
9.5 Diseo de un regulador PID real
Figura 9.5-1
Se pide disear el regulador real ms sencillo que verifique:
M p 15 , t s 1.5 s , e p 10
Ante entradas escaln en
X ( s) .
PD
t s= 1.5 s 2.1
Es decir:
W D= tan =3.5
Y
PD =2.1 j 3.5
Figura 9.5-2
PD
s +c
s+d
Figura 9.5-3
El ngulo que debe aportar este par polo-cero para que
races se obtiene aplicando el criterio del argumento.
( +1 +2 )=180
PD
Situar c en la vertical de
(c=2.1)
PD
Situar c y d tal que su bisectriz sea el ngulo buscado (Ver Figura 9.5-4):
Figura 9.5-4
Utilizando el segundo criterio, c=2 y d se obtiene aplicando el criterio del argumento tal y como se
refleja en la Figura 9.5-5 (d=3.2):
Figura 9.5-5
Para calcular el valor de la ganancia, se aplicara el criterio del mdulo:
2
s+ 2
s+ 3.2
K P =13.46
eP =
2 1
=4.21
3.2 1.2
1
=0.19=19
1+ K P
Que no es menor del 10 % pedido. As pues hace falta un regulador PID de la forma:
R ( s )=13.46
s+ a s +2
s+ b s +3.2
Con b< a , y (a, b) cercanos entre s. Se elegir a pequeo para que el lugar de las races estn
apreciablemente (a 1/6 de la distancia de los polos dominantes en cadena cerrada en el eje
imaginario):
a=
2.1
=0.35
6
1
1+ K P
K P =9=13.46
a 2 1
b 3.2 1.2
Es decir, b=0.16 y:
R ( s )=13.46
s+ 0.35 s +2
s+ 0.16 s +3.2
Figura 9.5-6
En lo referente al rgimen transitorio, la respuesta real del sistema aparece en la Figura 9.5-7:
Figura 9.5-7
9.6 Control de un sistema inestable
Se dispone de un sistema fsico cuya funcin de trasferencia es:
G ( s )=
1
(s+1)(s +2)(s1)
Se pide disear el regulador ideal ms sencillo que verifique las siguientes condiciones:
M P =5 , t s 2.86 s
Ante entradas escaln en la referencia. Comprobar si, con el regulador diseado, el sistema tendr
las condiciones de comportamiento dinmico requerido.
Solucin al problema 9.6
La posicin deseada de los polos dominantes
M p=e
tan
PD
100 5 45
t s= 2.86 s 1.1
Es decir:
W D= tan =1.1
Y
PD =1.1 j 1.1
Figura 9.6-1
Como el lugar de las races no pasa por
PD
inestabilidad el sistema, hay que modificar el lugar de las races aadiendo un regulador PD:
R ( s )=K R (s+b)
Figura 9.6-2
180arctan
][
1.1
1.1
1.1
1.1
arctan
+180arctan
+180arctan
=180
1.1b
0.9
0.1
2.1
Es decir:
180arctan
1.1
=118.25
1.1b
Valor que se encuentra entre 0 y 180, y por lo tanto b=0.5. Para calcular el valor de la ganancia se
aplica el criterio del mdulo:
K=
d 2d d 1
1
d b
K=K R1
, es decir
K R =2.97
R ( s )=2.97 ( s+0.5 )
El nuevo lugar de las races es el mostrado en la Figura 9.6-3
Figura 9.6-3
Sin embargo, si se calcula la situacin de la tercera raz:
Ahora
2 arctan
1.1
=118.25
b1.1
Valor que se encuentra entre 0 y 180, y por lo tanto b=1.75. Aplicando de nuevo el criterio de
modulo:
1.12+ 0.92 1.12+ 0.12 1.12+ 2.12
K=
=2.28
1.12 +0.6 2
Y
K R =2.28
, por lo que
R ( s )=2.28 ( s +1.75 )
Figura 9.6-4
El nuevo sistema siempre es estable. Sin embargo, los dos ceros modifican algo dinmico, tal y
como se muestra en la Figura 9.6-5:
Figura 9.6-5
9.7 Diseo de un regulador PID por el primer mtodo de Ziegler-Nichols
Se desea controlar un sistema trmico G(s) utilizando un regulador PID, como se indica en la Figura
9.7-1. Para ello se registra, mediante el transductor de temperatura de funcin de transferencia H(s),
la respuesta del sistema al escaln unitario (sistema G(s)*H(s) en lazo abierta), siendo esta la de la
Figura 9.7-2 y la Figura 9.7-3 (Detalle). Utilizando el procedimiento de Ziegler-Nichols, encontrar
el regulador adecuado.
Figura 9.7-1
Figura 9.7-2
Figura 9.7-3
Solucin al problema 9.7
Sobre las Figuras 9.7-2 y 9.7-3 se ha trazado la tangente por el punto de inflexin. En consideracin:
T v =2 segundos
T v +T a=30.6 segundos
KP=
1.2 T a
=17.16
K Tv
T i =2 T v =4
T d=0.5 T v =1
Por tanto:
La Figura 9.7-4 recoge la respuesta temporal al escaln del sistema as controlado. Observando
cmo se consigue una relacin de amortiguamiento (sobreoscilacion 1/ sobreoscilacion 2) cercana a
4:
Figura 9.7-4
9.8 Diseo de un regulador PID por el segundo mtodo de Ziegler-Nichols
Se desea controlar un sistema de funcin de transferencia desconocida. Para ello se dispone del
esquema de control de la Figura 9.8-1, pretendindose sintonizar el regulador R(s) mediante el
Ti
Td
segundo procedimiento de Ziegler-Nichols. Prescindiendo de las acciones
y
se aumenta
la ganancia
KP
Figura 9.8-1
Figura 9.8-2
Solucin al problema 9.8
En la Figura 9.8-2 puede leerse:
Se escoge un regulador PID obtenindose de la tabla de Ziegler-Nichols (segundo mtodo)
K P =0.6 K CR =0.12
T i =0.5 T CR =3.5
T d=0.12 T CR=0.84
Con lo que el regulador resulta:
R ( s )=K P (1+
0.12(1+
1
+0.84 s)
3.5 s
1
+T s )
Tis d
0.1(s+0.476)(s +0.715)
s
La Figura 9.8-3 recoge la respuesta del sistema as compensado. Como se observa excesivamente
lente y oscilatoria (aunque es estable, lo que no ocurrira si R(s)=1).
KP
la salida contra el error sera una lnea recta con una pendiente de
KP
en la banda proporcional
(figura 10.1)
100
[3]
banda proporcional
[4]
K P G P (s)
1+ K P GP ( s)
KP
Ejemplo 1
Si la planta de la figura 10.3 tiene una funcin de transferencia de
GP (s)=
1
s( s+1)
Se usa con control proporcional, Cul ser a) el tipo de sistema, b) los errores en estado estable
cuando se usa con i) una entrada escaln ii) una entrada rampa?
Respuesta
a
Go (s)=
KP
s ( s+1)
e ss
1
(s)]
1+Go ( s ) 1
1
1
}
s
KP
1+
s ( s +1 )
1
KP
Con el control integral la salida del controlador es proporcional a la integral de la seal de error e
con el tiempo, es decir.
t
SalidaK i e d ( t ) [5]
0
Donde
Ki
1
es la constante denominada ganancia integral. Esta tiene unidades de s . La figura
10.4 muestra que pasa cuando el error es de la forma de un escaln. La integral entre 0 y t es, de
hecho, el rea bajo la grfica del error entre 0 y t. as, debido a que despus de que el error
comienza, el arca se incrementa en una razn regular, la salida del controlador se debe incrementar
en una razn regular. La salida en cualquier tiempo es, entonces, proporcional a la acumulacin de
los efectos de los errores pasados.
Al tomar la transformada de Laplace de la ecuacin [5] da por resultado la funcin de transferencia,
para el controlador integral, de
Gc ( s ) =
salida ( s ) K i
= [6]
s
e (s)
As, para el sistema de la forma que se ilustre en la figura 10.5, el control integral de una funcin de
transferencia de la trayectoria directa de (
Ki
s
( )
Go ( s )=
3
,
, etcetera
nm nm
Ejemplo 2
Si la planta de la figura 10.5 tiene una funcin de transferencia de
1
s( s+1)
Y se usa con el control integral, Cul ser a) el tipo de sistema, b) los errores en estado estable
cuando se usa con i) una entrada escaln, ii) una entrada rampa y c) como se compara la estabilidad
con la que se presentara si hubiera control proporcional (como el ejemplo 1)?
Respuesta
a
Ki
2
s (s +1)
e ss
s1
(s)]
1+Go ( s ) 1
1+
e ss =lim {
s0
1
=0
s1
Ki
s 2 ( s+1 )
s
1
}
De este modo, para una entrada escaln un sistema de tipo 1 tiene un error en estado estable cero.
ii) para una entrada rampa, la entrada es 1/ s
1+
e ss =lim {
s1
Ki
s 2 ( s+1 )
1
}
s2
s0
1
=0
De esta manera, para una entrada rampa hay un error en estado estable cero, lo cual mejora la
situacin del ejemplo 1, donde se tena un control proporcional.
c
Para la situacin del control proporcional del ejemplo 1 el sistema tena una funcin de
transferencia de
KP
s ( s+ 1 )
G(s)=
KP
1+
s ( s+1 )
KP
s ( s+1 ) + K P
KP
es mayor que 0.
Para la situacin del control integral el sistema tiene una funcin de transferencia de
KP
2
G( s)=
s ( s+ 1 )
K
1+ 2 P
s ( s +1 )
Ki
2
s ( s +1 ) + K i
s +s + K i
El arreglo de Routh (vase el captulo 8) para este es
De este modo, un cambio de control proporcional a control integral, en esta instancia, da por
resultado la inestabilidad. La figura 10.6 muestra los diagramas del lugar geomtrico de las races
para las dos situaciones: diagrama a) con el control proporcional y b) con el control integral.
La reduccin en la estabilidad relativa como resultado de usar el control integral se puede resolver,
como una extensin, mediante el control proporcional integral (PI. Figura 10.7). Para tal
combinacin la salida del controlador es
t
S alidaK p e + K i e d ( t ) [8]
0
Figura 10.7
Gc ( s ) =
Ki
s
s K p+ Ki
s
K p [s+
Gc ( s ) =
Ki
Kp
( )
[ ( )]
K p s+
Gc ( s ) =
Ki
]
Kp
( )
1
i
[9]
K p s+
Gc ( s ) =
( 1 )] G ( s) [10]
p
De esta manera, mediante la funcin de transferencia el uso del control PI se adicionan un cero en -
( 1 )
i
error en estado estable para una entrada escaln. Debido a que se introducen un nuevo polo y un
nuevo cero, la diferencia entre el nmero de polos n y un numero de ceros m permanece sin cambio.
As, los ngulos de las asntotas para los lugares geomtricos de las races no cambian.
Angulos de las asintotas=
nm
Sin embargo, el punto de interseccin de las asntotas con el eje real se mueve hacia el origen y, en
consecuencia, se presenta cierta reduccin en la estabilidad relativa:
Punto de interseccion=
cambia por +
1
i
()
( 1 )
i
Ejemplo 3
Si la planta de la figura 10.7 tiene funcin de transferencia de
G p ( s )=
1
s (s +1)
[ ( )]
[ ( )]
[ ( )]
s
G ( s )=
K p s+
1+
1
G p (s )
i
s
K p [ s+0.5 ]
2
s 3 +s 2 + K p s+ 0.5 K p=0
El arreglo de Routh es:
s3
Kp
s2
0.5 K p
s1
0.5 K p
s0
0.5 K p
0.5 K p
del control proporcional al control integral resulta de la restitucin de la estabilidad. La figura 10.9
ilustra el diagrama del lugar geomtrico de las races para el sistema. Este se debera comparar con
el diagrama del lugar geomtrico de las races de la figura 10.6 cuando los controles proporcional e
integral estn separados. AL compararlo solo con el control proporcional, la estabilidad relativa se
reduce con el control PI, las asntotas estn ms cercanas al eje imaginario, pero comparado solo con
el integral el control PI ha empujado los lugares geomtricos de las races hacia la izquierda del
plano s.
Control Derivativo
Con la forma derivativa del controlador, la salida del controlador es proporcional a la razn de
cambio con el tiempo del error e, es decir
Salida=K d
de
[11]
dt
Kd
cuando hay un error de entrada rampa. Con el control derivativo, tan pronto como la seal de error
inicia puede haber una salida de controlador muy grande, puesto que esta es proporcional a la razn
de cambio de la seal de error y no a su valor. De este modo puede proporcionar una accin
correctiva grande antes de que se presente un error grande en realidad. Sin embargo, si el error es
constante entonces no hay accin correctiva, aun si el error es grande. As, el control derivativo es
insensible a las seales de error constantes o que varan con lentitud y, en consecuencia, so se usa
solo, sino con otras formas de controlador.
Al tomar la transformada de Laplace de la ecuacin [11] resulta, para el control derivativo, una
funcin de transferencia salida (s)/e(s).
Gc ( s ) =K d s [12]
Por lo tanto, para el Sistema de lazo cerrado que muestra la figura 10.11, la presencia del control
derivativo produce una funcin de transferencia en lazo abierto de:
K d S G p (s )
[13 ]
1+K d S G p ( s)
GD ( s ) =K d
([ 1 )+s ]G ( s) [14]
p
Kp
Donde d = K d
s=
1
d
Tampoco habr cambios en el tipo de sistema y, por lo tanto, en los errores en estado estable.
Gc ( s ) =
K d S G p (s )
1+ K d S G p ( s)
[( ) ]
K (s +0.5)
1
+ s G p ( s )= d
d
s( s+1)
(s)
i
b. Para una entrada rampa,
es s 2 y el error en estado estable es:
1
(s)
1+G D (s) i
1
1
K (s +0.5) s2
1+ d
s (s+1)
e ss =lim s
s0
e ss =lim s
s0
e ss =
1
0.5 K d
G ( s )=
K d [ s +0.5 ] G p (s)
1+ K d ( s+0.5) G p (s )
K d [ s +0.5 ]
G ( s )=
1
s( s+ 1)
1+ K d ( s+0.5)
G ( s )=
1
s (s +1)
K d [ s +0.5 ]
s (s +1)+ K d (s+0.5)
s +(K d )s+0.5 K d =0
El arreglo de Routh para este sistema es:
s2
0.5 K d
s1
1+ K d
s0
0.5 K d
Kd
es
positiva. La figura 10.13 muestra el diagrama del lugar geomtrico de las races para el
sistema. Esto se debera comparar con el diagrama de la figura 10.6 cuando el control
proporcional y el integral estn por separado.
Control PID
El controlador proporcional integral derivativo (PID), mejor conocido como controlador de
tres trminos, con un sistema de la forma que ilustra la figura 10.14 dar una salida, para una
entrada de error es de:
t
Salida=K p e+ K i edt + K d
0
de
[15]
dt
Ki
+ K d s [16 ]
s
, es
Kd
Kp
( )
(
G c ( s ) =K p 1+
, es
Kp
Ki
( )
Kd s
1
+
[17]
Kps Kp
G o ( s )=G c ( s ) G p ( s )K p 1+
Ki
+ s G p ( s )
is d
K p ( i s +1+ i d s2 ) G p ( s )
G o ( s )=
[18]
i s
1
s (s +1)
K p ( s 2+2 s+1 ) G p ( s ) K p ( s2 +2 s +1 )
Go ( s )=
=
2s
2 s 2( s +1)
G o ( s )=
K p ( s+1 )
2s
[ ]
[ ]
s+ 1
2 s2
(
)
G s=
s+1
1+ K p
2 s2
Kp
G ( s )=
K p ( s+1)
2 s + K p (s+1)
2
s2
Kp
s1
Kp
s0
Kp
Kp
es positiva.
La figura 10.15 ilustra el diagrama del lugar geomtrico de las races para el sistema.
Ki
proporcional,
y la ganancia
, la ganancia integral,
Ki
, y la ganancia derivativa
Kd
. La seleccin de
estas variables permite localizar polos y ceros que introduce el controlador a ser determinados y, por
lo tanto, afectan la estabilidad del sistema de control.
Para describir el proceso de seleccin de los mejores valores para el controlador se usa el trmino
sintonizacin. Existen varios mtodos para lograrlo, pero aqu solo se estudiaran dos mtodos, los de
Ziegler y Nichols. Ambos mtodos se basan en experimentacin y anlisis, estos son recetas de
cocinar tiles que se usan con mucha frecuencia. El primer mtodo a menudo se denomina me todo
de la curva de reaccin del proceso. El procedimiento con este mtodo consiste en abrir el lazo de
control de modo que no se presenten acciones de control. En general la ruptura del lazo se hace entre
el controlador y la unidad de correccin. Se aplica entonces, una seal de prueba a la unidad de
correccin y se determina la respuesta de la variable de proceso de medida, es decir la seal de error.
La seal de prueba deber ser tan pequea como sea posible. La figura 10.16 muestra la forma de la
seal de prueba y una respuesta tpica. La grafica de las seal medida se grafica contra el tiempo y
se conoce como curva de reaccin del proceso.
Kp
Ki
P/RL
0.9 P/RL
1.2P/RL
1/3.33L
1/2L
Kd
0.5L
El otro mtodo se conoce como el mtodo de la ltima ganancia . Primero, las acciones integral y
Kp
derivativa se reducen a sus valores mnimos. La constante proporcional,
, se fija en un valor
bajo y entonces, se incrementa en forma gradual. Esto es lo mismo que decir que la banda
proporcional se hace ms angosta de manera gradual. Mientras esto sucede, al sistema se le aplican
pequeas perturbaciones. El proceso continua hasta que se presentan oscilaciones. Se anota el valor
crtico de la constante proporcional, Kpc y Tc de estas. La tabla 10.2 muestra los criterios de Ziegler
y Nichols sobre cmo se seleccionan los valores de Kpc y Tc para establecer los valores del
controlador. La banda proporcional critica es 100/ Kpc .
Tabla 10.2 Criterios de Ziegler y Nichols para la ltima ganancia.
Modo de control
Proporcional solamente
Kp
0.5Kpc
Ki
Kd
Proporcional + Integral
Proporcional + Integral + Derivativo
0.45Kpc
0.6Kpc
1.2/Tc
2.0/Tc
Tc/8
Ejemplo 6
Determinar los valores que deben tener
Kp
Ki
Kd
modos, a partir de la curva de reaccin del proceso de la figura 10.17 cuando la seal de prueba fue
un 6% de cambio en la posicin de la vlvula de control.
K i=
1.2 P
1.26
=
=2.82
RL 0.017150
1
1
=
=0.0033 s1
2 L 2150
K d =0.5 L=0.5150=75 s
Ejemplo 7
Cuando se sintonizo un controlador de tres modos en un sistema de control, mediante el mtodo de
la ltima ganancia, se encontr que las oscilaciones iniciaban cuando la banda proporcional decreca
hasta un 30%. Las oscilaciones tienen un periodo de 500s. Cules son los valores apropiados de
Kp
Ki
Kd
,
,
?
Respuesta
El valor critico de
K pc
2
2
=
=0.004 s1
T c 500
Kd=
T c 500
=
=62.5 s
8
8
Realimentacin de velocidad
En muchos sistemas se involucra el posicionamiento de algn objeto, por ejemplo, el brazo de un
robot; en este caso , el requerimiento es que el sistema responda con rapidez a los errores y no
producir excesivas oscilaciones , o sobre pasos. Esto se puede lograr incorporando dentro del lazo
principal de realimentacin una trayectoria de realimentacin menor que introduzca lo que se
denomina realimentacin de velocidad, trmino que describe un lazo de realimentacin donde la
seal que se alimenta de regreso no es el valor de salida sino la razn de cambio con el tiempo de
salida. Para esta realimentacin la salida de la trayectoria se relaciona con su entrada mediante:
Salida=K v
d o
[19]
dt
G c (s) G p ( s)
[ 21]
1+G p K v s
Gc (s )G p (s)
[22]
1+ G p K v ( s ) +Gc ( s)G p (s)
G p ( s)k v s
en el
Go ( s )=
Kp
s ( s+ a)
1+
K vs
s (s +a)
Go ( s )=
Kp
s (s +a)+ K v s
G o ( s )=
Kp
[23]
s (s +a+ K v )
[ ]
100
1 2
Sobrepaso=exp
cos
100
1cos 2
cos
100
sen 2
[ ]
1
100 [24 ]
tan
Kp
, una planta
de funcin de transferencia:
G p ( s )=
1
s (s +1)
Go ( s )=
Kp
s (s +a)
1+
Kp
s (s +1)
Kp
s (s +1)+ K p
s +s + K p=0
Las races de esta ecuacin son:
1 14 K p
=0.5 K p 0.25
2
Por lo tanto la frecuencia natural angular,
wn
, es (figura 10.20)
w n20.52 + ( K p0.25 )
As, para
w n=2 rad /s
, entonces
K p=4.0
0.5
=0.25
2
Kp
s (s +a+ K v )
1K
( v )
geomtrico de las races, entonces, se convierte como muestra en la figura 10.21. Entonces:
=cos =
0.5(1+ K v )
wn
wn
0.5 ( 1+ K v ) =0.5 2
Y as,
K v =1
Compensacin
Los compensadores se pueden definir como componentes insertos en el sistema de control para
aumentar el desempeo del controlador. Si se considera que el controlador va a tener como base un
K (s+ z)
[25 ]
(s+ p)
Cuando z>p se conoce como compensador de atraso en cascada, y cuando z<p, como compensador
de adelanto en cascada. De esta manera, ambos compensadores introducen un cero y un polo; sin
embargo, entre los dos difiere la posicin relativa del polo y el cero. Un compensador de atraso en
cascada introduce el polo en lazo abierto ms cercano al origen que el cero. En muchos aspectos este
es como un control proporcional integral; sin embargo, la principal diferencia es que el polo que se
introduce no est en el origen. Del mismo modo que con el control PI ,el punto de interseccin de
las asntotas sobre el eje real se mueve a la derecha y as existe alguna reduccin en la estabilidad.
Los ngulos de las asntotas permanecen sin alteracin, puesto que no hay cambio en (n-m). Sin
embargo, puesto que no hay factor 1/s no hay cambio en el tipo de sistema. Un compensador de
adelanto en cascada introduce el cero ms cercano al origen que el polo. En muchos aspectos este es
como un control proporcional derivativo, pero la principal diferencia es que el cero no se introduce
en el origen. Como con el control PD, el punto de interseccin de las asntotas sobre el eje real se
mueve a la izquierda y, as , existe una mejora en la estabilidad. Los ngulos de las asntotas
permanecen sin alteracin, puesto que no hay cambio en (n-m). No se introduce un factor 1/s y, de
esta manera, no hay cambio en el tipo de sistema.
Ejemplo 9
Mostrar el efecto sobre el diagrama del lugar geomtrico de las races de un sistema con una funcin
de transferencia en lazo abierto de:
GO ( s )=
K
s (s +1)
Gc ( s ) =
s+ 2
s+ 8
ngulos de
consecuencia, el diagrama del lugar geomtrico de las raies es como se describe en la figura
10.22a .
K ( s+ 2)
s (s +1)( s+ 8)
De esta manera hay polos en lazo abierto en 0, -1 y -8 con un cero en -2. Las asntotas estn en
ngulos de
e intersecan el eje real en -7/2, puesto que n-m=2. Los lugares geomtricos siguen
al eje real entre 0 y -1, y entre -2 y -8. El punto de desprendimiento del eje real esta aun alrededor de
-0.5. De esta manera, el diagrama del lugar geomtrico de las races es como se ilustra en la figura
10.22b .
Puesto que la interseccin de las asntotas con el eje real se ha movido de -0.5 hasta -3.5 existe una
mejora en la estabilidad relativa.
V1= I1 R1
El amplificador operacional tiene una impedancia de entrada muy alta y de esta manera virtualmente
no fluye corriente hacia su interior a travs del punto x. Por lo tanto, la corriente, I 1, fluye solo a
travs de R2. Puesto que x es la tierra virtual, y debido a que la diferencia de potencial a travs de
R2 es Vs - V0. Por lo que:
V o=I 1 R2
Funcion de Transferencia=
V O R2
= [26 ]
V 1 R1
R2
R1
. El signo
negativo indica que la salida esta invertida, desfasada 180 grados respecto a la entrada.
Como indica la ecuacin [26], el amplificador operacional inversor tiene una ganancia de
R1
R2
de otro amplificador
operacional inversor.
R2=R1
R2
[27 ]
R1
1
idt
C
I (s)
Cs
1
[28 ]
Cs
Z 2
[29]
Z1
1
Z 2 ( s )=
Z
(
s
)
=R
1
1
Puesto que
y
Cs
, entonces:
Funcion de transferencia=
1
R1 Cs
Cuando el circuito se combina con otro circuito con amplificador operacional de funcin de
transferencia -1, como muestra la figura 10.25, da por resultado:
1
R C
Funcion de transferencia= 1 [30 ]
s
Ki=
1
R1 C
1
Cs
Funcion de transferencia=
Con
K p=
R2
R1
K i=
R 2
1
R1 C
1
Cs
R1
1
R RC
= 2 + 1 [31]
R1
s
La figura 10.27 muestra cmo se puede producir un controlador derivativo. Para el amplificador
operacional con capacitor y resistor en la lnea de entrada se tiene, segn la ecuacin [29] da por
resultado:
Funcion de transferencia=
R s Cs
R1 Cs+1
R2 Cs
[ 32]
R1 Cs+1
Kd
de
R2 Cs
Funcion de transferencia=
R s Cs
R1 Cs+1
Con
K p=
R2
R1
K i=
1
R1 C
K d =R4 C 2
arreglo tobera-aleta produce la salida correspondiente a un error cero. Cuando la presin del proceso
cambia de este valor, la aleta gira y cambia el espacio entre la aleta y la tobera. El resultado es un
cambio en la presin de salida. Esta presin cambia hasta que los fuelles de realimentacin ejercen
una fuerza para balancear aquella debida a los fuelles del proceso.
Donde
Kp
A 1 d 1=A 2 d 2
La figura 10.31 describe la forma bsica de un controlador proporcional neumtico. Una diferencia
en la presin entre los fuelles del proceso y del punto de ajuste produce un movimiento de aleta.
Esto cambia la presin de salida y la de los fuelles del proporcional. El resultado es un movimiento
de la aleta hasta que el momento de rotacin, debido a la fuerza que en los fuelles del proporcional,
se balancea a aquella que resulta de la diferencia de presin entre los fuelles de proceso y del punto
de ajuste. Mientras esto ocurre, los fuelles del integral se afectan en forma ligera debido al tiempo
de retardo que introduce la restriccin. Sin embargo, a medida que los fuelles del integral se acercan
lentamente a la presin de salida, se mueve la aleta y, en esta forma, cambia la presin de salida.
PROBLEMAS
1. Un sistema realimentado tiene una funcin de transferencia de la trayectoria directa de:
1
2
s ( s +3 s+5 )
Y realimentacin unitaria. Cul ser el error en estado estable con una entrada rampa unitaria
si
a) En la trayectoria directa, se introduce un controlador proporcional con ganancia4,
b) En lugar del controlador proporcional se usa uno integral constante de tiempo 2s?
2. Qu ceros y polos se introducen en un sistema realimentado si en la trayectoria directa se
introduce un controlador PI con una ganancia proporcional de 4 y una constante de tiempo
integral de 2s?
3. Un sistema realimentado tiene una funcin de transferencia en lazo abierto de :
G O ( s )=
s +3
s ( s+1 ) ( s+2 )
Cunto se mover el punto de interseccin de las asntotas de los lugares geomtricos de las
races con el eje si se introduce un controlador PI con una constante de tiempo integral de 5s
y una ganancia proporcional de 2?
4. Un sistema realimentado tiene un control proporcional derivativo con una planta que tiene
una funcin de transferencia en lazo abierto
1
s ( s+2 )
Seleccionar los valores de Kp y Kd de modo que el sistema tenga una frecuencia natural
angular de 0.5 rad/s y un factor de amortiguamiento relativo de 0.7.
5. Un controlador PID se usa para controlar una planta con una funcin de transferencia de
1
s ( s+3+ j1 ) ( s +3 j 1 )
Si el controlador PID tiene una constante de tiempo integral de 4s y una constante de tiempo
derivativa de 1s. Cules sern los ceros y polos en lazo abierto del sistema de control?
6. Determinar los valores idneos de Kp, Ki, Kd requeridos para un controlador de tres modos
que tiene una curva de reaccin del proceso con un retraso, L, de 2000s y un gradiente, R, de
0.010%/s cuando la seal de prueba fue un cambio de 5% en la entrada unitaria de
correccin.
7. Al sintonizar un controlador de tres modos mediante el mtodo de la ltima ganancia se
encontr que las oscilaciones iniciaron cuando la banda proporcional se increment a 20% y
estas tenan un periodo de 200s. Cules son los valores apropiados para Kp ,Ki ,Kd ?
8. Describir los efectos en un sistema de control al incluir:
a. Un control proporcional
b. Un control integral
c. Un control proporcional integral
d. Un control proporcional derivativo
e. Un control PID.
9. Explicar que significa realimentacin de velocidad y explicar la razn para usar dicho
termino.
10. Cul es el valor mximo de k, el cual dar las condiciones de estabilidad para que el sistema
que muestra la figura 10.33 con realimentacin de velocidad y posicin?
Fig10.33
K
( s +3 ) ( s+ 5 )
Cuando
a) Sin compensar
b) Con compensador de adelanto con una funcin de transferencia de (s+1)/(s+2)
c) Con un compensador de atraso con una funcin de transferencia de (s+2)/(s+1)
13. Explicar cmo se puede usar un AOP para producir las formas de control proporcional,
integra y derivativo.