UNIVERSIDAD DEL BÍO-BÍO

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA

Ejercicios
410 377 Modelación y Sim. de Proc. y Mec.

Ejercicio 1.
Para el tanque de calentamiento de la figura, se considera que no hay perdidas de calor y el
revolvedor garantiza temperatura homogénea.

F1(t), T1(t)
ρ, cp

T2(t)

h(t)
Q(t)

F2(t), T2(t)
A
vapor

Se realiza el siguiente análisis para determinar su modelo dinámico.

Análisis de variables.
• Flujo volumétrico de entrada F1(t): puede ser manipulado por la válvula de flujo. Se
mide en [m3/h]
• Flujo volumétrico de salida F2(t): sus variaciones alrededor del punto de operación
se considerarán proporcionales a las variables de la presión hidrostática en el fondo
del tanque. Se mide en [m3/h]
• Temperatura de entrada T1(t): depende del proceso anterior. Se mide en [°C].
• Temperatura en el Tanque T2(t): se considera homogénea en el tanque por el
revolvedor e igual a la de la salida. Variará en función de T1(t), Q(t) y la masa del
fluido dentro del tanque. Se mide en [°C].
• Calor suministrado por unidad de tiempo Q(t): puede ser manipulado por la válvula
de vapor y se mide en [kcal/h].
• Nivel del liquido en el tanque h(t): variará según F1(t) y F2(t). Se mide en [m].

Profesor Daniel Quezada., Septiembre 2023

UNIVERSIDAD DEL BÍO-BÍO
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA

Clasificación de variables.
• Entradas manipulables: Flujo volumétrico de entrada F1(t) y Calor suministrado Q(t). Se
pueden modificar a total voluntad.
• Entrada perturbadora: Temperatura de entrada T1(t). Depende del medio exterior.
• Salidas (objetivos de control): En términos de producción, Flujo volumétrico de salida
F2(t). En términos de calidad, Temperatura en el Tanque T2(t). En términos de seguridad,
Nivel del líquido en el tanque h(t).
• Parámetros: Densidad del fluido ρ (kg/m3), Calor especifico del fluido cp (kcal/kg°C),
Área de la sección transversal del tanque (m2).
Modelación.
• Balance de masa.
d ( m (t ))
G1 ( t ) − G2 ( t ) =
dt
d (  Ah ( t ) )
 F1 ( t ) −  F2 ( t ) =
dt
d ( h (t ))
F1 ( t ) − F2 ( t ) = A
dt

• Balance de energía.

d ( Q2 ( t ) )
Q1 ( t ) + Q ( t ) − Q2 ( t ) =
dt
d (  Ah ( t ) c pT2 ( t ) )
 F1 ( t ) c pT1 ( t ) + Q ( t ) −  F2 ( t ) c pT2 ( t ) =
dt
1 d ( h ( t ) T2 ( t ) )
F1 ( t ) T1 ( t ) + Q ( t ) − F2 ( t ) T2 ( t ) = A
cp dt

• Relación algebraica.

F2 ( t ) = k p ( t ) − patm
F2 ( t ) = k  h ( t ) g

Puntos de operación.
Conocidos/deseados: F1o, T1o, T2o, ho
F1o − F2 o = 0  F1o = F2 o

Profesor Daniel Quezada., Septiembre 2023

UNIVERSIDAD DEL BÍO-BÍO
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA

1
F1oT1o + Q − F2 oT2 o = 0  Qo =  c p F1o (T2o − T1o )
cp o

F2o
F2o = k  ho g  k =
 ho g

Diagrama en bloques.
El diagrama en bloques que resume el proceso a partir de sus ecuaciones dinámicas se
presenta a continuación.
F1(t) h(t)
1 1
+
A s
_

F2(t)
k ρg10-5

_
T1(t) 1 1
× + A s
×
÷
T2(t)
Q(t) 1
c p

Profesor Daniel Quezada., Septiembre 2023

 UNIVERSIDAD DEL BÍO-BÍO
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA

Simulación.
a) Analizar el siguiente código y replicar en un Script de Matlab.
1 % Fluido agua
2 p = 1000; % densidad (kg/m3)
3 g = 10; % gravedad (m/s2)
4 cp = 1; % calor específico (kcal/kgC)
5
6 A = 1; % área del tanque (m2)
7
8 % punto de operación
9 F1o = 1; % flujo de entrada (m3/min) (variable de entrada manipulable)
10 T1o = 50; % temperatura de entrada (C) (variable de entrada perturbadora)
11 T2o = 70; % temperatura de salida deseada (C) (variable de salida)
12 ho = 1; % nivel en el tanque (m) (variable de salida)
13
14 % del balance de masa en estado estacionario (d(...)/dt = 0) se llega a:
15 F2o = F1o;
16 % de la ecuación de relación, con la presión en bar, para que se cumpla el
17 % punto de operación, la tubería de salida se debe diseñar para que:
18 k = F2o/sqrt(p*g*1e-5*ho);
19 % del balance de energía en estado estacionario (d(...)/dt = 0) se llega a:
20 Qo = p*cp*F1o*(T2o-T1o);
21
22 % analizar lo que ocurre en las salidas (T2 y h) cuando varían las entradas:
23
24 % variación en el flujo de entrada
25 % en t=tdF1 el F1 va de F1o a F1o+dF1
26 dF1 = 0.5;
27 tdF1 = 10;
28
29 % variación en la temperatura de entrada
30 % en t=tdT1 la T1 va de T1o a T1o+dT1
31 dT1 = 5;
32 tdT1 = 30;
33
34 % variación en el calor suministrado en t=tdQ el Q va de Qo a Qo+dQ
35 dQ = 5000;
36 tdQ = 50;

Profesor Daniel Quezada., Septiembre 2023

UNIVERSIDAD DEL BÍO-BÍO
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA

b) Implementar el modelo para la salida de nivel.

El bloque del Flujo volumétrico de entrada y el Integrador deben configurarse de la siguiente

manera:

Profesor Daniel Quezada., Septiembre 2023

UNIVERSIDAD DEL BÍO-BÍO
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA

La simulación entregará las siguientes curvas:

c) Completar el modelo agregando la parte térmica del proceso.

Profesor Daniel Quezada., Septiembre 2023

UNIVERSIDAD DEL BÍO-BÍO
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA

Las entradas deben configurarse de la siguiente manera:

La simulación entregará las siguientes curvas:

Profesor Daniel Quezada., Septiembre 2023