“Año del Fortalecimiento de la Soberanía Nacional”
AUTOMATIZACION Y CONTROL
INDUSTRIAL
LABORATORIO CALIFICADO N° 2
GRUPO 4
PROFESOR: RONALD DUEÑAS PÉREZ
ALUMNOS:
• MOQUILLAZA ROMANI JOSUE DAVID JUNIOR 1510495
• TOROBEO PARIGUANA IVAN 1020585
• CAMPOS BAYONA JUAN JOSUE 1411260
TEMA: SIMULACION DE LA DINAMICA DE TANQUE DE
CALENTAMIENTO
2022
PREGUNTA 1
SIMULACION DINAMICA:
clear;clc
global dx
% RANGO DE INTEGRACION
xi = 0;
xf = 20;
% VALORES INICIALES DE "y" Y "u"
T = 40;
T1 = 20;
Q = 0;
% PARAMETROS
F = 100;
V = 500;
rho = 1;
Cp = 1;
% SECCION DE INICIACION
% ----------------------------------------
dx = 0.0001;
XPRINT = 1.5; % INTERVALO DE IMPRESION DE DATOS cada
1.5 minutos
X = xi:dx:xf;
n = length(X);
Y = T; % INICIALIZA EL VECTOR Y
%----------------------------------------------------
for i=1:n-1
% SECCION DE ECUACIONES
%-----------------------------------------------
dT = F/V * (T1 - T) + Q/(rho*V*Cp);
% SECCION DE INTEGRACION
%-----------------------------------------------
T = T + dT*dx;
Y(i+1) = T;
end
% SECCION DE PRESENTACION DE RESULTADOS
%---------------------------------------------------
% GRAFICA DE APROXIMACION DE EULER
plot(X,Y,'b-','LineWidth',1)
xlabel('x')
ylabel('y')
axis square
ax = gca;
ax.FontName = 'ARIAL';
ax.FontSize = 12;
hold
% GRAFICA DE SOLUCION ANALITICA
% Y_A = u*(1-exp(-X/Tau));
% plot(X,Y_A,'k--','LineWidth',1)
% legend('Aprox. Euler', 'Sol. Exacta')
%hold off
% TABLA DE RESULTADOS
X_Y = PRINT_INT(X,Y,XPRINT)
RESPUESTA: Entre el minuto 7.5 y minuto 9 la temperatura será
de 24°C
PREGUNTA 2
SIMULACION DINAMICA:
clear;
clc
global dx
% RANGO DE INTEGRACION
xi = 0;
xf = 10;
% VALORES INICIALES DE "y" Y "u"
CA = 5;
CB = 1;
CA1 = 12;
% PARAMETROS
K1 = 5/3;
K2 = 5/6;
V = 8;
F = 5;
% SECCION DE INICIACION
% ----------------------------------------
dx = 0.01;
XPRINT = 0.5; % INTERVALO DE IMPRESION DE DATOS cada 30 SEGUNDOS
X = xi:dx:xf;
n = length(X);
Y = CA; % INICIALIZA EL VECTOR Y
Z = CB;
%----------------------------------------------------
for i=1:n-1
% SECCION DE ECUACIONES
%-----------------------------------------------
r1 = K1*CA;
r2 = K2*CB;
dCA = F/V *(CA1-CA) - r1;
dCB = -F/V*CB + r1 - r2;
% SECCION DE INTEGRACION
%-----------------------------------------------
CA = CA + dCA*dx;
CB = CB + dCB*dx;
Y(i+1) = CA;
Z(i+1) = CB;
end
% SECCION DE PRESENTACION DE RESULTADOS
%---------------------------------------------------
% GRAFICA DE APROXIMACION DE EULER
subplot(2,1,1),plot(X,Y,'b-','LineWidth',1)
xlabel('t(min)')
ylabel('CA(mol/l)')
axis square
axis ([0 xf 0 5])
ax = gca;
ax.FontName = 'ARIAL';
ax.FontSize = 12;
hold
subplot(2,1,2),plot(X,Z,'b-','LineWidth',1)
xlabel('t(min)')
ylabel('CB(mol/l)')
axis square
axis ([0 xf 0 10])
ax = gca;
ax.FontName = 'ARIAL';
ax.FontSize = 12;
hold
% TABLA DE RESULTADOS
X_Y = PRINT_INT(X,Y,XPRINT)
X_Z = PRINT_INT(X,Z,XPRINT)
TABLA Ca
RESPUESTA: Para Ca en el minuto 5 alcanza su valor final
TABLA Cb
RESPUESTA: Para Cb en el minuto 6.5 alcanza su valor final
PREGUNTA 3
MODELO DINAMICO
FUNCION DE TRANSFERENCIA
Funcion de transferencia para el tanque 1
Usando LAPLACE:
Funcion de transferencia del tanque 2
EJEMPLO:
Hallar F1, Fm y F2,
𝑎1 = 0.5
𝑎2 = 0.45
𝑘1 = 0.04 𝑚3 /𝑠
𝑘2 = 0.03 𝑚3 /𝑠
𝑘3 = 0.055 𝑚3 /𝑠
ℎ1 = 9 𝑚
ℎ2 = 16 𝑚
Hallando F1
𝐹1 = 𝑘1 𝑎1
𝐹1 = 0.04 ∗ 0.5 = 1200 l/min
𝐹𝑚 = 𝑘2 √ℎ1 = 0.03√9 = 0.09 = 5400 𝑙/𝑚𝑖𝑛
𝐹2 = 𝑘3 𝑎2 √ℎ2 = 0.055 ∗ 0.45 √16 = 0.099 = 5940 𝑙/𝑚𝑖𝑛
OBSERVACIONES
• Según la pregunta 1, considerando como temperatura inicial 40°C y que la cantidad de
calor disminuye al 50%, tenemos que aproximadamente a los 5 minutos el líquido tiene
una temperatura de 36°C, y que pasado 10 minutos llega a bajar hasta 34°C. Esto quiere
decir que al minuto 1 la temperatura disminuye en 1.2 y pasado otro minuto, disminuye
en 1.1, con eso se tiene que la temperatura tiene una tendencia de reducción de 0.1 en
cada minuto comenzando desde 1.2.
• En la pregunta 2 tenemos que el Ca va disminuyendo en comparación del Cb que va en
aumento.
• En la pregunta 3 se tuvo que averiguar para obtener la función de transferencia del
modelo dinámico de tanque en cascada en función del área y altura de ambos tanques,
así como también de los valores K de cada uno.
CONCLUSIONES
• Mediante el software MATLAB podemos determinar la temperatura del líquido en un
minuto determinado modificando la cantidad de calor que se requiera.
• Así también determinar los valores de Ca y Cb cada rango de minuto y determinar en
cuanto disminuye.
• Para determinar la función de transferencia de un modelo de tanque en cascada, se
procede hallando la función de transferencia de cada tanque.
BIBLIOGRAFIA
https://pdfcoffee.com/sistema-de-control-de-nivel-de-dos-tanques-en-cascada-2-pdf-
free.html
