Tarea 3 C

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 9

UNIVERSIDAD DE G U A D A L A J A R A . CUCEI.

DIVISION DE INGENIERIAS. DEPARTAMENTO DE INGENIERIA QUIMICA


CONTROL DE P R O C E S O S . Tarea 3. Entregar 13 de marzo de 2018.
..

1. Considere el ejemplo de un tanque de nivel con flujo de salida proporcional a la raíz


cuadrada del nivel.
a) Utilizando el balance de masa, desarrolle el modelo del sistema.
b) Encuentre la aproximación lineal del modelo anterior y expréselo en variables de
desviación.
c) Compare l a respuesta de ambos modelos a un cambio del flujo F1 de 1 lt=seg,
gra…cando ambas respuestas en la misma Figura.
Datos: Estado Estacionario: F1s = 30; F2s = 10; = 11:6190; A = 10; tiempo de
intergracion de 0 a 10,
Anexar a la solucion de este problema, los comandos necesarios de MATLAB para
generar la ecuacion diferencial, integrarlas y graficarlas.

function p21
clear
clc

t0=[0 50];
h0=0;

[t,hs]=ode45(@p2,t0,h0);

plot(t,hs)

end

%estado estacionario
function dhsdt=p2(~,hs)

As=10;
F1s=30;
F2s=10;
Bs=11.619;

dhsdt=(F1s+F2s-Bs*sqrt(hs))/As;
hold on
end

1
function p2v
clear
clc

t0=[36 50];
h0=11.84;

[t,h]=ode45(@p5,t0,h0);

plot(t,h)

end

function dhdt=p5(~,h)

A=10;
F1=31;
F2=10;
B=11.619;

dhdt=(F1+F2-B*sqrt(h))/A;
end

2
2. Considere el caso de tres tanques de nivel colocados en serie.
a) Desarrolle el modelo matematico de este sistema dinamico sabiendo que se ali-
mentan dos corrientes al primer tanque (F y F ) y que las corrientes de salida son
proporcionales al nivel del tanque.
b) Integre las ecuaciones diferenciales utilizando MATLAB usando como condiciones
inciales un valor 50% del valor en el estado estacionario.
function P2T3

tspan=[0,100];
h10=0;
h20=2;
h30=1;

%Solucion Ecuacion de diseño


[h,t]=ode45(@P2T,tspan,[h10,h20,h30]);

subplot(1,1,1),plot(h(:,1),t)
subplot(1,1,2),plot(h(:,2),t)
subplot(1,1,3),plot(h(:,3),t)

end

function dhdt=P2T(~,h)

x=h(1);
T=h(2);
P=h(3);

%Datos calculados
%Datos
a1=5;
a2=10;
a3=6;
3
F1=10;
F2=10;
A1=10;
A2=15;
A3=12;

%Ecuacion de Diseño
dhdt(1)=(F1+F2-(a1*x))/A1;
dhdt(2)=((a1*x)-(a2*T))/A2;
dhdt(3)=((a2*T)-(a3*P))/A3;
dhdt=dhdt';

end

c) Gra…que su solucion en la misma Figura.


Anexar a la solucion de este problema, los comandos necesarios de MATLAB para
generar la ecuacion diferencial, integrarlas y gra…carlas
DATOS: F = 10; F = 10; 1 = 5; 2 = 10; 3 = 6; A1 = 10; A2 = 15; A3 = 12

3. Para l os sistemas dinámicos lineales siguientes:


a) Determine los polos.
b) Clasi…que el estado estacionario de acuerdo al tipo y carácter de su estabilidad.
c) Dibuje su plano de fase utilizando al menos 10 condiciones iniciales.
.

4
dx 1 2 dx 6 8
a) = x d) = x
dt 3 4 dt 2 2
dx dx 1 1
1 2 = x
b) = x e) dt 5 3
dt 3 4

dx 1 3 dx 2 6
c) = x f) = x
dt 4 2 dt 1 2
.

function p3t3
clear
clc

t0=[0 10];
x10=-30;
x20=0;

x0=[x10 x20] ;

[t,x]=ode45(@p3t,t0,x0);

plot(t,x(:,1),t,x(:,2))
plot(x(:,1),x(:,2))

end

function dxdt=p3t(~,x)
global a11 a12 a21 a22
a11=1;%input('Valor a11= ');
a12=-2;%input('Valor a12= ');
a21=3;%input('Valor a21= ');
a22=-4;%input('Valor a22= ');

dxdt(1,:)=a11*x(1,:)+a12*x(2,:);
dxdt(2,:)=a21*x(1,:)+a22*x(2,:);

end

a)

5
_____________________________________

b)

6
_____________________________________

c)

__________________________________ 7
d)
_______________________________

e)

_______________________________________

8
f)

______________________________________

4. Considere de nuevo el problema 5 de la Tarea 1 que es un reactor químico continuo


de tanque agitado (RCTA) en donde se está llevando a cabo una reacción de primer
orden no isotérmica, solo que ahora se encuentra operando en estado no estacionario.
.
a) Cual es ahora el modelo que describe este proceso dinamico (expreselo en variables
adimensionales)
b) Integre las ecuaciones diferenciales para al menos 5 condiciones iniciales con los
datos siguientes
b.1 = 8:7 = 0:15 = 30 y b.2 = 25 = 0:25 = 30
c) Gra…que los resultados en el denominado plano de fase ( x(t) vs y(t) ) para todas
las condiciones iniciales.
.
Anexar a la solucion de este problema, los comandos necesarios de MATLAB para
generar la ecuacion diferencial, integrarlas y gra…carlas

También podría gustarte