PRACTICO N.
º 1
ESTUDIO DE SISTEMAS COMPLEJOS DE
TUBERÍAS EN EL ENTORNO SIMUSOL
“MODELO ANALOGICO DE UN SISTEMA DE TUBERIAS – CIRCUITO ELECTRICO”
INTRODUCCION
El modelo analógico es una herramienta fundamental en la ingeniería civil, utilizada para representar
sistemas complejos mediante modelos físicos o matemáticos con características similares, pero en
una forma diferente. Esta metodología permite analizar y predecir el comportamiento de estructuras
y materiales antes de su construcción, facilitando la toma de decisiones informadas y la optimización
de diseños. En este ensayo, se exploran las aplicaciones del modelo analógico en la ingeniería civil,
destacando su importancia en el diseño estructural y la simulación de fenómenos físicos.
El estudio de Rodríguez et al. (2002) analiza analogías entre sistemas hidráulicos y eléctricos, con el
objetivo de incorporar nuevos elementos al software de simulación de sistemas solares térmicos
SIMUSOL, basado en el programa Sceptre para circuitos eléctricos. Se demuestra que, aunque solo
existe equivalencia directa entre la ley de Ohm y el flujo laminar, las ecuaciones de régimen
turbulento pueden ajustarse para ser interpretadas mediante la misma metodología.
El trabajo resalta la complejidad del análisis de redes de tuberías (en serie, paralelo y ramificadas) y
la necesidad de métodos computacionales para optimizar los cálculos y reducir errores. Se propone
una metodología alternativa para integrar estos estudios en SIMUSOL, comparando resultados entre
circuitos hidráulicos y sus análogos eléctricos.
La analogía hidráulico-eléctrica permite modelar redes de tuberías como circuitos eléctricos, donde
la presión equivale al voltaje, el caudal a la corriente y la resistencia hidráulica a la resistencia
eléctrica. Esta aproximación facilita la resolución computacional de sistemas complejos, mejorando
la precisión y eficiencia en el análisis de flujo en conductos cerrados.
ANALOGÍAS PARA CIRCUITOS HIDRÁULICOS ESPECIALES
➢ Tuberías en serie
El análisis de este tipo de situación se realiza en función de los datos disponibles. Para el caso como
el de la figura, en el que tenemos una tubería horizontal, con estrechamiento y ensanchamiento brusco
en su parte intermedia, la determinación de los valores de energía perdida en cada uno de los puntos
indicados, puede determinarse mediante la aplicación de la ecuación de Darcy – Weisbach, las
fórmulas empíricas de Hazen –Williams, o la de Manning (Mataix, 1992). También los nomogramas,
gráficos, tablas y diagramas, proporcionan una solución rápida, aunque con menor precisión en los
resultados.
Desde el punto de vista de la hidráulica, el procedimiento seguido es el siguiente: conociendo el valor
de la energía en el extremo inicial de la tubería, la energía en la sección 3 será igual a la que tenía en
la sección 2 menos la pérdida entre los puntos 2 y 3 (𝐻3 = 𝐻2 − ℎ2−3 ). Particularizando esta ecuación
para cada tramo elegido en la tubería problema y adoptando la expresión adecuada para cada tipo de
pérdida, se puede determinar los niveles de energía totales en todas las secciones, luego restando la
porción correspondiente a energía cinética se obtienen los llamados niveles piezométricos.
Resolución mediante circuitos análogos
Este enfoque establece una analogía entre circuitos hidráulicos y eléctricos para facilitar el análisis
de sistemas de tuberías en serie mediante simulación computacional.
Componentes del Modelo Análogo:
1. Pérdidas de Carga: Representadas por resistencias eléctricas:
⎯ RF (Resistencia por fricción) → Pérdida por rozamiento en la tubería.
⎯ RC (Resistencia por contracción brusca) → Pérdida por reducción de sección.
⎯ RE (Resistencia por ensanchamiento brusco) → Pérdida por aumento de sección.
2. Energía y Flujo:
⎯ Fuente de Energía (E) → Representa la energía total inicial del sistema hidráulico.
⎯ Fuente de Corriente (JG) → Representa el caudal circulante en la tubería. Se coloca
al final del circuito para simular el flujo.
3. Instrumentación de Medición:
⎯ Caudalímetro (Amperímetro eléctrico, R7) → Se representa con una resistencia
muy pequeña y mide el flujo en una sección.
⎯ Medidor de Presión (Voltímetro eléctrico, JP) → Se representa con una fuente de
corriente de valor nulo y mide la presión en un punto.
⎯ Medidor de Nivel Piezométrico (Ev + JPE) → Representado por una fuente de
energía cinética en sentido opuesto al flujo y una fuente de corriente de valor nulo.
Cálculo y Simulación:
1. El programa de simulación calcula la caída de tensión en cada resistencia eléctrica, lo que
equivale a determinar la pérdida de carga en cada tramo de la tubería.
2. A partir de estos cálculos, se puede trazar la línea de energía y la línea piezométrica a lo largo
del sistema.
3. Los valores de las resistencias hidráulicas y fuentes de energía se determinan previamente con
datos reales, mediante cálculos manuales o en hojas de cálculo. También se pueden programar
en el software de simulación utilizando instrucciones similares a FORTRAN.
4. Como el modelo opera en el Sistema Internacional (SI), las presiones obtenidas estarán en
Pascales (Pa), por lo que, para expresarlas en metros de columna de agua, se requiere una
conversión auxiliar.
➢ Tuberías en paralelo
Cuando varias cañerías se unen de manera tal que conforman un circuito en paralelo, es necesario
verificar dos condiciones principales:
⎯ El caudal Q se reparte entre todas las tuberías.
⎯ La presión al comienzo y al final de cada rama es la misma.
Los casos más comunes son aquellos en los se conoce la pérdida de carga y se quiere determinar los
caudales en cada rama, o cuando se tiene como dato el caudal total y se desea determinar tanto la
pérdida de carga como la distribución de caudales. El problema se puede resolver empleando los
Métodos del Porcentaje y el de la Tubería Equivalente (Giles et al., 1996), utilizados especialmente
cuando el fluido circulante es agua.
Método de la Tubería Equivalente
1. Suposición Inicial:
⎯ Se adopta un valor arbitrario de pérdida de carga (𝐻) entre los puntos de bifurcación y unión.
2. Cálculo de Pérdidas de Carga por Metro:
⎯ Se determina la pérdida de carga por unidad de longitud (ℎ𝑓 ) en cada rama de la red usando
ecuaciones como Darcy-Weisbach o Hazen-Williams.
3. Determinación de los Caudales en Cada Rama:
⎯ Con la pérdida de carga por metro encontrada, se calculan los caudales que fluyen por cada
rama y se suman para obtener el caudal total (𝑄𝑇𝑂𝑇𝐴𝐿 ).
4. Definición de la Tubería Equivalente:
⎯ Se adopta un diámetro igual al mayor de las ramas y se calcula la longitud equivalente que
tendría esta tubería para producir la misma pérdida de carga (𝐻) cuando circula el caudal total.
5. Cálculo de la Nueva Pérdida de Carga por Metro:
⎯ Se determina la nueva pérdida de carga por unidad de longitud en la tubería equivalente al
hacer circular el caudal total.
6. Corrección y Ajuste de Caudales:
⎯ Con la nueva pérdida de carga calculada, se recalculan los caudales en cada rama y se ajustan
iterativamente hasta que los valores sean coherentes.
Resolución mediante circuito análogo
1. Modelo hidráulico transformado a circuito eléctrico:
El sistema hidráulico se modela como un circuito eléctrico, donde:
⎯ La presión en la cañería se representa como una fuente de tensión.
⎯ El caudal circulante por la cañería se representa como una fuente de corriente.
⎯ Las pérdidas por fricción (RF) en las cañerías se modelan como resistencias en el circuito
eléctrico.
2. Configuración de los circuitos:
Hay dos configuraciones para el sistema de cañerías:
⎯ Figura 1: Cuando se conoce la presión (PA y PB) en los puntos de la cañería. Aquí, las fuentes
de tensión (E1, E2) se introducen en el circuito eléctrico, y las resistencias hidráulicas (RF1,
RF2, RF3) se calculan a partir de las características de las cañerías y las presiones en los
puntos.
⎯ Figura 2: Cuando se conoce el caudal total (Q). En este caso, las fuentes de corriente (I1, I2,
I3) se usan, y las resistencias hidráulicas se modelan de manera similar, pero en función del
caudal.
3. Cálculo de la resistencia hidráulica (RF):
El cálculo de la resistencia hidráulica se propone como el producto de dos términos:
𝑅𝐹 = 𝐾 ⋅ 𝐼𝑅𝐹
Donde:
⎯ K es un factor relacionado con las características de la cañería (como su longitud, diámetro,
rugosidad, etc.).
⎯ IRF es el caudal a través de la resistencia hidráulica.
Cuando el caudal en las ramas del sistema es desconocido (es decir, las incógnitas son los caudales),
la resistencia RF no es una constante, sino que depende del caudal. El programa de simulación espera
que la resistencia RF sea un número ya calculado previamente, pero como ahora depende del caudal,
se necesita una forma de resolver esto.
4. Cambio de variable en la simulación:
Para resolver el problema del caudal desconocido, se realiza un cambio de variable. En lugar de usar
𝑉𝑅𝐹
directamente la fórmula de resistencia 𝑅𝐹 = 𝐾 ∙ 𝐼𝑅𝐹, se redefine IRF como , lo que permite
𝑅𝐹
escribir la resistencia en función del caudal de una forma más adecuada para el programa.
Se introduce una nueva fórmula para calcular la resistencia RF:
𝑉𝑅𝐹
𝑅𝐹 = 𝐾 ⋅ √
𝐾
Este cambio de variable permite incorporar la resistencia en el programa de simulación, que sigue
una secuencia de cálculo iterativo.
5. Problema de indeterminación y solución en el proceso iterativo:
Al empezar el proceso de simulación, el programa asigna un valor de 0 a VRF (la variable que
depende del caudal). Al calcular RF, el programa obtiene RF = 0, lo que lleva a una indeterminación
cuando intenta calcular IRF, ya que la resistencia es cero y la división entre cero no es válida, lo que
genera un error en el proceso de cálculo.
Para resolver este problema y evitar que el programa se detenga, se modifica la ecuación de la
1
resistencia para incorporar el término 1+𝑇𝐼𝑀𝐸.Esto asegura que, en la primera iteración (cuando TIME
= 0), la resistencia inicial RF sea igual a 1, evitando la indeterminación.
1
⎯ En la primera iteración: RF=1 (debido al término 1+𝑇𝐼𝑀𝐸 cuando TIME = 0).
⎯ En las iteraciones posteriores: El valor de RF aumenta en pequeñas cantidades conforme
avanza el proceso iterativo, lo que permite continuar con la simulación sin errores.
6. Ecuación final en FORTRAN:
El programa usa un lenguaje de programación (FORTRAN) para calcular las resistencias y caudales.
La ecuación implementada para calcular la resistencia es:
𝐵 1
𝑄𝑇𝑈𝐵(𝐴, 𝐵) = (𝐴 ∙ √ + )
𝐴 1 + 𝑇𝐼𝑀𝐸
Donde:
• A = K (factor relacionado con las características de la cañería),
• B = VRF (caudal a través de la resistencia),
• TIME es el tiempo de la iteración.
Resumen del proceso:
1. Se transforman los elementos hidráulicos en equivalentes eléctricos (resistencias, fuentes de
tensión y corriente).
2. Se modela la resistencia hidráulica RF como una función del caudal.
3. Se hace un cambio de variable para evitar la indeterminación en el cálculo.
4. Se modifica la ecuación de la resistencia para asegurar que la simulación avance
correctamente.
5. El programa sigue un proceso iterativo para calcular las soluciones del sistema hidráulico.
➢ Redes de tuberías
Las tres leyes principales:
⎯ Ley de la pérdida de carga: La pérdida de energía en una cañería depende del caudal que
circula a través de ella. Esta pérdida de carga se calcula usando un coeficiente β, que depende
de las características de la cañería (diámetro, longitud, material, etc.) y del caudal.
⎯ Ley de los nudos: En un "nudo" (un punto donde se unen varias cañerías), la cantidad de
fluido que entra es igual a la cantidad que sale. Es decir, no se pierde ni se gana fluido en un
nudo; solo se distribuye.
⎯ Ley de las mallas: En un circuito cerrado (una "malla"), la suma de las pérdidas de carga a
lo largo de esa malla debe ser igual a cero. Esto significa que la energía que se pierde en las
cañerías se debe equilibrar dentro de la malla.
Método de Hardy Cross:
Para resolver el sistema de cañerías, el método de Hardy Cross comienza con una aproximación
inicial (generalmente se asume un caudal unitario). Luego, se calcula el coeficiente de pérdida para
cada cañería, utilizando la ecuación de Karman-Prandtl.
Se realizan los siguientes pasos:
1. Se calcula la pérdida de carga en cada cañería.
2. Se suman las pérdidas de carga de cada malla y, si no es cero, se realiza una corrección
(cambiando los caudales en cada cañería).
3. La corrección se calcula mediante una fórmula y se ajustan los valores de los caudales.
4. Este proceso se repite varias veces, hasta que la diferencia de caudales sea muy pequeña y se
pueda considerar que el sistema ha convergido a una solución.
Aplicación en un circuito eléctrico análogo:
En este método, la red hidráulica se convierte en un circuito eléctrico equivalente, donde:
⎯ Las cañerías se representan por resistencias (que dependen del caudal y las características
de las cañerías).
⎯ Los caudales que entran o salen de los nodos se representan con fuentes de corriente en el
circuito eléctrico.
⎯ Importante: No es necesario colocar una fuente de corriente en el nodo de entrada del
fluido (el primer nodo), porque la corriente que llega a este nodo es la suma de las
corrientes que entran a través de las otras cañerías conectadas a él. Es decir, todo el caudal
total que entra al sistema se representa por las fuentes de corriente en los nodos de salida.
CONCLUSION
El análisis de redes hidráulicas complejas mediante modelos eléctricos análogos, como se demuestra
con el uso del programa SIMUSOL para simular sistemas solares, muestra que es completamente
posible representar estos sistemas de forma eficaz. Los resultados obtenidos en ejemplos
representativos y en diversas configuraciones problemáticas confirman que se pueden construir
circuitos eléctricos análogos a sistemas hidráulicos complejos, lo que facilita el diseño y la simulación
de dichos sistemas.
Si bien la simulación de tuberías en serie es relativamente sencilla, ya que los datos necesarios para
calcular los elementos del circuito están disponibles, el verdadero desafío radica en las ramificaciones
de tuberías. En estos casos, la distribución del caudal en las diferentes ramas no siempre es conocida,
lo que complica el cálculo de las resistencias hidráulicas, ya que estas son funciones del caudal. Por
ello, es crucial incorporar adecuadamente la relación entre resistencia y caudal en el modelo,
utilizando técnicas como el método de Hardy Cross y considerando la naturaleza dinámica de las
resistencias.
El desarrollo y la interpretación adecuada de este tipo de modelos son esenciales para lograr
simulaciones precisas y eficientes. Además, la correcta integración de la resistencia hidráulica como
función del caudal dentro de las rutinas de cálculo permite que el programa de simulación realice
ajustes precisos en la distribución del caudal y las pérdidas de carga, proporcionando una herramienta
útil para el diseño y análisis de redes hidráulicas complejas. En conclusión, estos enfoques permiten
mejorar la precisión y eficiencia del diseño de sistemas hidráulicos, lo cual es fundamental para el
desarrollo de infraestructuras sostenibles y eficientes.