Saltar al contenido
Control Automático Educación

Modelo de 2 tanques esféricos

Hola controleros y controleras sean bienvenidos una vez más a otra entrada de la pagina de Control Automático Educación, donde vamos a continuar con nuestro curso de Análisis de Sistemas. Hoy vamos a aprender como obtener el modelo matemático de un sistema de dos tanques esféricos en cascada a través de ecuaciones diferenciales, vamos a linealizar el sistema en un punto de operación para obtener sus representaciones en función de transferencia y en espacio de estados y finalmente simularemos la planta en MATLAB/SIMULINK.

Antes que nada te invito a VER EL CURSO GRATUITO DE ANALISIS DE SISTEMAS.

Modelo de 2 tanques esféricos

A continuación se muestra un sistema de 2 tanques esféricos, el cual conforma un sistema MIMO (multi-Input Multi-Output) con dos variables controladas h_1 y h_2 y dos variables manipuladas F_{i_1} y F_{i_2}.

Modelo tanques esfericos

La dinámica de los tanques puede ser modelada usando el principio de conservación de Masa:

\dfrac{dm}{dt}=\dfrac{dm_{in}}{dt}-\dfrac{dm_{out}}{dt}

Se asume que el fluido dentro de los tanques es incompresible y que el fluido posee una densidad constante, de esa forma las ecuaciones diferenciales de los dos tanques son formulados como la variación del volumen en ambos tanques

(1) \dfrac{dV_1}{dt}=F_{i_1}-F_{o_1}

(2) \dfrac{dV_2}{dt}=F_{i_2}-F_{o_2}

Donde V_i es el volumen acumulado en cada tanque y F_{i} y F_{o} son los caudales volumetricos de entrada y salida de cada tanque respectivamente. Dabiendo que dV=Adh:

(3) \dfrac{A_1(h_1)dh_1}{dt}=F_{i_1}-F_{o_2}

(4) \dfrac{A_2(h_2)dh_2}{dt}=F_{i_2}-F_{o_2}

Además, el área de la superficie libre de fluido en los tanques 1 y 2 podrían expresarse en función de las de alturas de los líquidos h_1 y h_2, asi tomando en cuenta que dV=Adh procedemos a derivar el volumen de un tanque esferico con relacion a la altura, sin embargo para esto deberemos conocer la ecuación matemática del volumen de un casquete esférico:

V=\dfrac{\pi h^2}{3}(3R-h)

Si no sabes obtener esta ecuación, te voy a dejar el procedimiento para el cálculo del casquete esférico a través del calculo integral:

Volumen del Casquete Esférico

El casquete esférico es el sólido generado a partir de una esfera cuando está seccionada por un plano:

Casquete esférico Tanques


Por lo tanto, los límites de integración serán de r - h a r.

Partiendo de la ecuación de la circunferencia con centro en el origen:

x=\sqrt{r^2+y^2}

Supongamos que el casquete esférico de altura h es formado por cilindros infinitos de alturas infinitesimales dy y radios x, donde x es variable para cada punto de h:

modelado casquete esférico


Sabemos que el volumen de un cilindro es:

V_c=A_bh=\pi r^2 h

Y en este caso:

V_c=\pi x^2 dy

La suma de estos cilindros infinitos de alturas infinitesimales forma el volumen del casquete esférico en los límites r - h y r. Entonces su volumen estará dado por:

V=\int_{r-h}^{r}\pi x^2 dy

V=\pi \int_{r-h}^{r} x^2 dy

Sustituyendo la ecuación del circulo tenemos


V=\pi \int_{r-h}^{r} (\sqrt{r^2+y^2})^2 dy

V=\pi \int_{r-h}^{r} (r^2+y^2) dy

Resolviendo

V=\pi\left[r^2y - \dfrac{y^3}{3} \right]_{r-h}^{r}

aplicamos limites

V=\pi\left[\left( r^2r - \dfrac{r^3}{3}\right) - \left( r^2(r-h) - \dfrac{(r-h)^3}{3} \right) \right]

V=\pi\left[\left( \dfrac{3r^3-r^3}{3}\right) - \left( \dfrac{3r^3-3hr^2-r^3-3hr^2+3rh^2-h^3}{3} \right) \right]

V=\dfrac{\pi h^2}{3}(3r-h)

Continuación del Modelo de los 2 Tanques Esféricos

Volviendo al ecuacionamiento dinámico de los dos tanques esféricos acoplados:

(5) A_1=\pi(2R_1h_1-h_1^2)

(6) A_2=\pi(2R_2h_2-h_2^2)

Los caudales de salida que pasan por una sección transversal y a través de una válvula pueden ser modelados como

(7) F_{o_1}=k_1a_1\sqrt{h_1}

(8) F_{o_2}=k_2a_2\sqrt{h_2}

De esa forma la ecuación diferencial que representa los tanques acoplados viene dado por:

(9) \pi(2R_1h_1-h_1^2)\dfrac{dh_1}{dt}=F_{i_1}-k_1a_1\sqrt{h_1}

(10) \pi(2R_2h_2-h_2^2)\dfrac{dh_2}{dt}=F_{i_2}+k_1a_1\sqrt{h_1}-k_2a_2\sqrt{h_2}

despejando el sistema tenemos:

(11) \dfrac{dh_1}{dt}=\dfrac{1}{\pi(2R_1h_1-h_1^2)}\left(F_{i_1}-k_1a_1\sqrt{h_1}\right)

(12) \dfrac{dh_2}{dt}=\dfrac{1}{\pi(2R_2h_2-h_2^2)}\left(F_{i_2}+k_1a_1\sqrt{h_1}-k_2a_2\sqrt{h_2}\right)

En el caso de los flujos de entrada, que provienen de las bombas, vamos a dejarlos indicados como F_{i_1} y F_{i_2}, sin embargo, por causa de usar bombas aquí, podemos representar esa entrada como:

F_{i_1}=K_pu_1(t) y F_{i_2}=K_pu_2(t) donde u_1(t) y u_2(t) son el voltaje aplicado a las bombas y el K_p es la constante de la bomba dada en unidades de caudal por voltaje, por ejemplo (L/v).

Estado Estacionario

Basta con cerar las derivadas y encontrar los estacionarios para ambas alturas.

Comenzamos con la altura 1

0=\dfrac{F_{i_1}^*}{\pi(2R_1h_1^*-(h_1^*)^2)}-\dfrac{k_1a_1\sqrt{h_1^*}}{ \pi(2R_1h_1^*-(h_1^*)^2)}

(13) h_1^*=\left(\dfrac{F_{i_1}^*}{k_1a_1}\right)^2

Ahora la altura 2

0=\dfrac{F_{i_2}^*}{\pi(2R_2h_2^*-(h_2^*)^2)}+\dfrac{k_1a_1\sqrt{h_1^*}}{\pi(2R_2h_2^*-(h_2^*)^2)}-\dfrac{k_2a_2\sqrt{h_2^*}}{\pi(2R_2h_2^*-(h_2^*)^2)}

h_2^*=\left(\dfrac{F_{i_2}^*}{k_2a_2}+\dfrac{F_{i_1}^*}{k_2a_2}\right)^2

(14) h_2^*=\left(\dfrac{F_{i_2}^*}{k_2a_2}\right)^2+\dfrac{2F_{i_1}^*F_{i_2}^*}{(k_2a_2)^2}+\left(\dfrac{F_{i_1}^*}{k_2a_1}\right)^2

Espacio de Estados

A partir de las ecuaciones (11) y (12) es posible obtener la representacón en espacio de estados del sistema linealizado usando la Jacobiana y empleando los puntos de equilibrio de las ecuaciones (13) y (14). Para eso reescribimos las ecuaciones (11) y (12) como:

(15) G_1(h_1,h_2,F_{i_1},F_{i_2})=\dfrac{1}{\pi(2R_1h_1-h_1^2)}\left(F_{i_1}-k_1a_1\sqrt{h_1}\right)

(16) G_2(h_1,h_2,F_{i_1},F_{i_2})=\dfrac{1}{\pi(2R_2h_2-h_2^2)}\left(F_{i_2}+k_1a_1\sqrt{h_1}-k_2a_2\sqrt{h_2}\right)

\begin{bmatrix}\dot{h_1}\\ \\\dot{h_2}\end{bmatrix}=\begin{bmatrix}\dfrac{\partial G_1}{\partial h_1} & \dfrac{\partial G_1}{\partial h_2}\\ \\\dfrac{\partial G_2}{\partial h_1} & \dfrac{\partial G_2}{\partial h_2}\end{bmatrix}\begin{bmatrix}{h_1}\\ \\ \\h_2\end{bmatrix}+\begin{bmatrix}\dfrac{\partial G_1}{\partial F_{i_1}} & \dfrac{\partial G_2}{\partial F_{i_2}}\\ \\\dfrac{\partial G_2}{\partial F_{i_1}} & \dfrac{\partial G_2}{\partial F_{i_2}}\end{bmatrix}\begin{bmatrix}{F_{i_1}}\\ \\ \\F_{i_2}\end{bmatrix}

donde:

\dfrac{\partial G_1}{\partial h_1}=\dfrac{-F_{i_1}(2\pi R_1-2\pi h_1^*)}{(2\pi R_1 h_1^2-\pi (h_1^*)^2)^2}+ \dfrac{ \dfrac{-k_1a_1}{2\sqrt{h1^*}}(2\pi R_1h_1^*-\pi (h_1^*)^2)+k_1a_1\sqrt{h_1^*}(2\pi R_1-2\pi h_1^*) } { (2\pi R_1 h_1^2-\pi (h_1^*)^2)^2 }

\dfrac{\partial G_2}{\partial F_{i_2}}=0

\dfrac{\partial G_2}{\partial F_{i_1}}= \dfrac{ \dfrac{k_1a_1}{2\sqrt{h_1^*}} }{ (2\pi R_1 h_2^2-\pi (h_2^*)^2) }

    \begin{align*}\dfrac{\partial G_2}{\partial F_{i_2}}= & \dfrac{1}{(2\pi R_2 h_2^2-\pi (h_2^*)^2)^2} \\& \left(-F_{i_2}(2\pi R_2-2\pi h_2^*)+ \dfrac{-k_2a_2}{2\sqrt{h_2^*}}(2\pi R_2h_2^*-\pi (h_2^*)^2)+k_2a_2\sqrt{h_2^*}(2\pi R_2-2\pi h_2^*) \right.\\& \left. -k_1a_1\sqrt{h_1^*}(2\pi R_2-2\pi h_2^*)\right)\end{align*}

\dfrac{\partial G_1}{\partial F_{i_1}}=\dfrac{1}{\pi(2R_1h_1^*-h_1^2)}

\dfrac{\partial G_2}{\partial F_{i_2}} = \dfrac{\partial G_2}{\partial F_{i_1}} = 0

\dfrac{\partial G_2}{\partial F_{i_2}}=\dfrac{1}{\pi(2R_2h_2^*-h_2^2)}

Comparación Modelo No Lineal y Modelo Lineal

La siguiente grafica muestra la comparación entre los dos modelos, en un punto muy cercano al punto de equilibrio donde fue linealizado el sistema, con una pequeña desviación de la entrada del tanque 1 de magnitud 0.1

Tanques Esfericos
Linealización Tanques Esfericos

Simulación de 2 Tanques Esféricos en Matlab

A continuación esta disponible el código que te permitirá modelar y simular los tanques esfericos acoplados que vimos en esta entrada. Este código tiene la capacidad de obtener el perfil dinámico de los tanques, obtener el modelo lineal en un punto de operación y finalmente hacer una comparación entre el modelo lineal y no lineal.

Para tener acceso al código basta con que compartas esta información para ayudar que este sitio web siga llegando a mas personas y continue compartiendo información de interés.

Los dos Archivos están disponibles aquí

>>>> BAJAR ARCHIVOS <<<<<

Código Principal

Modelo de los 2 Tanques Esféricos

Bibliografia

  • Tavakolpour-Saleh, A. ; Setoodeh, A. ; Ansari, E. (2016), ‘Iterative Learning Control of Two Coupled Nonlinear Spherical Tanks ‘, World Academy of Science, Engineering and Technology, Open Science Index 119, International Journal of Mechanical and Mechatronics Engineering, 10(11), 1862 – 1869.

Eso es todo por la entrada del dia de hoy, espero les haya gustado y hayan aprendido algo nuevo. Si te ha servido el contenido de esta entrada, de los videos y los códigos de implementación y deseas apoyar mi trabajo invitandome a un café super barato, puedes hacerlo en el siguiente link:

👉 Invitar a Sergio a un Café ☕️

Que esten muy bien, nos vemos en la siguiente entrada.