Scribd
Cargar
15 vistas5 páginas

4ta Practica de Control

Este documento presenta un modelo matemático de un reactor químico en serie con reacciones A → B → C. Se pide linealizar el sistema de ecuaciones diferenciales que describe la variación de las concentraciones de A y B usando como variables de entrada el flujo, las concentraciones de alimentación de A y B, y como salidas las concentraciones de A y B.

Cargado por

rarocer
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
15 vistas5 páginas

4ta Practica de Control

Este documento presenta un modelo matemático de un reactor químico en serie con reacciones A → B → C. Se pide linealizar el sistema de ecuaciones diferenciales que describe la variación de las concentraciones de A y B usando como variables de entrada el flujo, las concentraciones de alimentación de A y B, y como salidas las concentraciones de A y B.

Cargado por

rarocer
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

PROBLEMA:

Considrese la operacin isotrmica de un reactor tanque agitado donde ocurre el


1 2

siguiente sistema de reacciones en serie A B C. El modelo matemtico del reactor,


el cual describe la variacin de las concentraciones de los componentes A y B, est
dado por:


= (0 ) 1 2


= (0 ) + 1 2 2

Linealizar este sistema de ecuaciones usando como variables de entrada el flujo


volumtrico Q, la concentracin del reactivo A en la corriente de alimentacin 0 , y la
concentracin de B en la misma corriente 0 . Las variables de salida, o respuesta del
sistema, son y .
Dato:

= 100 /

= 250 /

= 20/

2 = 0.5 1.

1 = 0.01 / ;

Solucin:
1. Modelo del sistema:
=
=

(0 ) 1 2

(0 ) + 1 2 2

2. Espacio de estado:
=
=

(0 1) 1. 12

(0 2) + 1. 12 2. 2

1
1
( ) = (
2
0

0 1
)( )
1 2

= 1200 ;

3. Linealizando en MatLab:
% Programa : LinPenInv.m
% Descripcion :
clc
% Variables simbolicas
syms f1 f2 x1 x2 Q Cao Cbo k1 k2 V
f1=(Q*(Cao-x1))/(V)-k1*(x1)^2
f2=(Q*(Cbo-x2))/(V)+k1*(x1)^2-k2*x2
f=[f1;f2];
% Calculo de jacobianos en Punto de Operacion
v=[x1,x2]; w=[Q,Cao,Cbo];
x1=0;x2=0;
As=subs(jacobian(f,v))
Bs=subs(jacobian(f,w))
% Dando valores a parametros del sistema
Q=100;Cao=250;Cbo=20;V=1200;k1=0.01;k2=0.5;
A=subs(jacobian(f,v))
B=subs(jacobian(f,w))
C=[1 0;0 1]
D=[0]
step(A,B,C,D)
% Fin

4. Ejecutando Programa:

f1 =
(Q*(Cao - x1))/V - k1*x1^2

f2 =
k1*x1^2 - k2*x2 + (Q*(Cbo - x2))/V

As =
[ -Q/V,
[

0]

0, - k2 - Q/V]

Bs =
[ Cao/V, Q/V, 0]
[ Cbo/V, 0, Q/V]

A=
-0.0833

0 -0.5833

B=
0.2083

0.0833

0.0167

C=
1

D=
0

0.0833

5. Sistema no Lineal en MatLab:


clc; clear all;
Cao=250;Cbo=20;V=1200;k1=0.01;k2=0.5;
t=0; h=0.001; f=[0; 0];
tacu=t; facu=f; Q=1;
for i=1:400
f=f+h*[(Q*(Cao-f(1)))/(V)-k1*(f(1))^2;(Q*(Cbo-f(2)))/(V)k1*(f(1))^2-k2*f(2)];
t=t+h;
tacu=[tacu t];
facu=[facu f]
end
%plot(tacu',yacu(1,:)','r'); grid on;
plot(tacu,facu(1,:),'r'); grid on;
axis([min(tacu) max(tacu) min(facu(1,:)) max(facu(1,:))]);
title('Sistema reactor de tanque');
xlabel('tiempo[seg]'); ylabel('Posicion DCa[mol/lt]');
plot(tacu,facu(2,:),'r'); grid on;
axis([min(tacu) max(tacu) min(facu(2,:)) max(facu(2,:))]);
title('Sistema reactor de tanque');
xlabel('tiempo[seg]'); ylabel('Posicion DCb[mol/lt]');
%-------------------------------------------------------------------------

Grafica 2: de la concentracin de DCa

Grafica 3: de la concentracin de DCb

También podría gustarte