Laboratorio N°4
Laboratorio N°4
Laboratorio N°4
2. MARCO TEORICO
1
u (t ) K P e(t )
Ti
e(t )dt T
0
de(t )
dt
(1)
e(t ) r (t ) y (t )
Donde
es la seal de error, u(t) es la seal o fuerza de control, y(t) es la
salida controlada del proceso o planta, r(t) es la seal de referencia deseada o set-point,
KP es la ganancia proporcional, Ti es la constante de tiempo integral, y Td es la constante
de tiempo derivativo. El algoritmo dado en la ecuacin (1) se denomina tambin algoritmo
ISA (Instrument Society of America) o algoritmo PID no interactuante.
FIEE -UNAC
Pgina 2
1
u (s) K P 1
Td s e( s )
Ti s
(2)
e( s ) r ( s ) y ( s )
Donde
Los algoritmos PID descritos operan con una seal de error e(t), la cual es generada en
lnea (On-Line), es decir, con el controlador en plena accin, en donde tambin se
muestran las seales de disturbio o perturbaciones l (en la carga) y n (en la medicin).
+
r(t) +
e(t)
Controlador PID
l(t)
u(t)
+
Planta
n(t)
y(t)
1
u1 (t ) K P e(t )
Ti
e(t )dt Td
dy (t )
dt
(3)
1
u1 ( s) K P 1
e( s ) K P Td s y ( s )
Ti s
(4)
Otra modificacin importante es filtrar la accin derivativa del controlador PID original
mediante un filtro de primer (o segundo) orden para disminuir el ruido derivativo. Esta
FIEE -UNAC
Pgina 3
Td s
;
1 Tf s
Tf
Td
N
K p Td s
1
u1 ( s ) K P 1
y ( s) P( s ) I ( s) D( s )
e( s )
Ti s
1 Tf s
(5)
Dnde:
P ( s ) K P e( s ); I ( s )
K T s
KP
e( s ); D( s ) P d y ( s )
Ti s
1 Tf s
r(k) +
Planta
y(k)
Sea k el ndice de tiempo discreto, el cual est relacionado con el tiempo continuo t
mediante la relacin
FIEE -UNAC
t kT
Pgina 4
P K P e(s )
parte proporcional del controlador PID en tiempo continuo dada por
representada en tiempo discreto por:
vendr
P K P e(k )
(6)
I ( s)
KP
e( s )
Ti s
I (t )
KP
Ti
e(t )dt
0
La parte integral
, o lo que es lo mismo
vendr
expresada en tiempo discreto por la siguiente aproximacin trapezoidal:
I (k )
I (k )
KP
Ti
KP
Ti
k 1
i 1
T
i 1
T
2
Ti
2
I (k 1)
KP
Ti
k 1
T
i 1
I (k ) I ( k 1)
K PT
[e(k ) e(k 1)]
2Ti
(7)
D(s)
La parte derivativa
K P Td s
y (s)
1 Tf s
se puede escribir como:
(1 T f s ) D ( s) K P Td s y ( s)
D(t ) T f D (t ) K P Td y (t )
FIEE -UNAC
Pgina 5
y (t )
e
D(k ) D(k 1)
y ( k ) y (k 1)
K P Td
T
T
D(t ) T f
Tf
Despejando D(k), y considerando
como:
D(k )
Td
N
, la parte derivativa discreta queda entonces
Td
{D(k 1) K P N [ y (k ) y (k 1)]}
( NT Td )
(8)
Por consiguiente, el algoritmo PID discreto usado en las implementaciones en tiempo real
tiene la siguiente forma:
u ( k ) P ( k ) I (k ) D ( k )
(9)
3. SISTEMA DE UN HORNO TUBULAR
La figura 1 representa un horno tubular, y lo que se desea es controlar la temperatura a un
nivel de referencia de 200 C. El nivel de temperatura se sensa por medio de un sensor
de temperatura, cuyos terminales se tienen disponibles. La seal de entrada del sistema
Horno tubular
(t )
q (t )
en lazo abierto es
, (t)
la salida disponible es v(t),
es el flujo de caudal de
(t )
combustible, y
Fuel
FIEE -UNAC
q
es la temperatura
en el interior del horno.
Sensor de
temperatura
Pgina 6
v(t)
q (t ) 2q(t ) (t )
(t ) 0.5 (t ) 2q (t )
v(t ) 0.25 (t )
4. PROCEDIMIENTO DE PROBLEMA
4.1 Determinar el modelo de funcin de transferencia, considerando u = (variable de
entrada o de control), y = v (variable de salida).
4.2 Determinar por transformacin, la funcin de transferencia pulso del modelo de
funcin de transferencia obtenido en (1).
4.3 Determinar su representacin en espacio de estado de tiempo continuo, teniendo
en cuenta:
u = (variable de entrada o de control)
y = v (variable de salida)
y las siguientes variables de estado:
x1 =
FIEE -UNAC
Pgina 7
4.4 Determinar la representacin en tiempo discreto del modelo obtenido en (3), para
lo cual use el mtodo de discretizacin directa, y un periodo de muestreo de 1
segundo.
4.5 Simule la respuesta del sistema, en tiempo continuo y discreto del modelo en
espacio de estado, que permita obtener grficas muy similares. En caso que no se
logre con dicho periodo de muestreo, entonces use otro valor.
4.6 Disee un Controlador PID de temperatura, que logre obtener una respuesta
adecuada en lazo cerrado. Considere que el factor de amortiguamiento es de 0.8,
y escoja bajo su criterio el tiempo de asentamiento. Puede usar alternativamente
uno de los mtodos de Ziegler-Nichols, segn corresponda, hasta lograr sintonizar
el PID y obtener respuestas de acuerdo a lo solicitado.
r(t)
+
-
KI
s
u(t)
x Ax Bu
y Cx Du
y(t)
x
K
sistema de control
Sugerencia: ubicar la tercera raz a 10 veces la parte real de las races dominantes.
FIEE -UNAC
Pgina 8
q (t ) 2q(t ) (t )
(t ) 0.5 (t ) 2q (t )
v(t ) 0.25 (t )
G p (s)=
y (s)
0.5
=
u(s) ( s+ 0.5 ) ( s+2)
[ ]
G( z)=
z1 G(s)
Z
z
s
G( z)=
z1
0.5
z1
Z
=
Z
z
z
s ( s +0.5 ) ( s+ 2)
Donde
FIEE -UNAC
[(
A
B
C
+
+
0.5
s s+0.5 s +2
Pgina 9
G( z)=
z1 0.5
2/3 1/6
Z
+
z
s
s+0.5 s +2
G( z)=
+
; T : Periodo
z
z1 z e0.5 T ze2 T
G(z )=0.5+
(2/3)
(1/6)
+
.( z 1)
0.5T
ze
ze2T
Entonces:
FIEE -UNAC
Pgina 10
) ()
x (t )= 0.5 2 x (t )+ 0 u(t )
0
2
1
x i (k +1) x i(k)
= x i (t )
T
Matricialmente:
) ()
2T x + 0 u
x ( k+1 )= 10.5 T
0
12T (k ) T ( k )
FIEE -UNAC
Pgina 11
) ()
x ( k+1)= 0.5 2 x( k ) + 0 u( k )
0 1
1
y(k )=( 0.25 0 ) x(k)
Dnde:
G= 0.5 2
0 1
()
H= 0
1
C d=( 0.25 0 ) ;
D d =( 0 )
FIEE -UNAC
Pgina 12
Step Response
0.5
0.45
0.4
0.35
Amplitude
0.3
0.25
0.2
0.15
0.1
0.05
0
10
20
30
40
50
60
Time (sec)
t s=8 seg .
t s=8 ;
t s=
=0,8
4
0.5
rad
=8 . n=0.5 n=
=0.625
. n
0.8
seg
Gc(s) =k p 1+
FIEE -UNAC
k p (T i . T d . s2 +T i . s+1)
1
+T d s =
()
Ti s
T i. s
Pgina 13
do
s +2 . n . s+ n =0
s 2 +2. ( 0,5 ) s+ ( 0.625 )2=0 s 2+1. s+0.39=0
s1 =0.5+ j .0.37
s2=0.5 j .0 .37
( s +2 . n . s+ n ( sp 3)=0
2
GLC (s) =
k p T d . s 2+ k p . s+
FIEE -UNAC
G( s )
Y (s)
=
; G =G c(s) .G p (s)
R(s) 1+G(s) ( S)
kp
Ti
0.5
( s+0.5 ) (s +2)
Pgina 14
G ( S )=
G LC(s )=
Y (s )
G(s )
=
=
R(s ) 1+G(s)
1+
De ( I ) y ( II ) :
k pDISCRETO =k p =8,78
FIEE -UNAC
Pgina 15
k iDISCRETO =
k i T (3,9)1
=
k iDISCRETO =1,95
2
2
k d DISCRETO =
k d ( 7)
=
k dDISCRETO=7
T
1
FIEE -UNAC
Pgina 16
Sabemos de la planta:
.
x(t ) Ax(t ) Bu (t )
.. (1)
y (t ) Cx(t )
.. (2)
Ademas:
u (t ) Kx(t ) K I (t )
... (3)
.
(t ) y (t ) r (t )
.............. (4)
De (2) y (4)
.
(t ) Cx(t ) r (t )
..(5)
De (1) y (5)
x (t ) A 0 x (t ) B u (t ) 0 r (t )
1
. C 0 (t ) 0
(t )
. (6)
Entonces:
x () A 0
.
C 0
FIEE -UNAC
x () B
0
u () r ()
( ) 0
1
..(7)
Pgina 17
.
.
Si:
x(t ) x () xe (t )
(t ) () e (t )
u (t ) u () u e (t )
x e (t ) A
.
C
(
t
)
e
0 x e (t ) B
u e (t )
0 e (t ) 0
u e (t ) K
xe (t )
(
t
)
e
K I
Definimos:
A 0
C 0
K K
,
FIEE -UNAC
Pgina 18
KI
.(8)
e(t ) A xe (t ) B u e (t )
. (9)
u e (t ) K e(t )
.. (10)
K K
K I k1 k 2
kI
MP 2%
y ts (5%)=6
Adems:
e
MP=
1 2
0.78
=0.02
ts(5%)
Luego
, obteniendo
3
wn
, entonces Wn =0.641
FIEE -UNAC
Pgina 19
e(t ) ( A B K )u e (t )
sI ( A B K ) 0
Donde:
2 0
0.5
A
0
2 0
0.25 0 0
0
B 1
0
K k1 k 2
,
2 0 0
0.5
A B K 0
2 0 1 k1
0.25 0 0 0
2
0
0.5
A B K k1
2 k2 kI
0.25
0
0
FIEE -UNAC
Pgina 20
k2
kI
kI
2
1 0 0 0 .5
s 0 1 0 k1
2 k2
0 0 1 0.25
0
kI 0
0
Igualando (
) y(
FIEE -UNAC
Pgina 21
FIEE -UNAC
Pgina 22
VI.- CONCLUSIONES:
FIEE -UNAC
Pgina 23