Análisis Básico de una Red Simple de Tuberías: Teoría y Resolución Iterativa

1. Configuración Básica de una Red Simple y Condiciones a Verificar

Una red cerrada (malla) se compone topológicamente de tres elementos geométricos:

• Nudos (o Nodos): Puntos donde convergen dos o más tuberías, o donde el fluido entra o sale del sistema.
• Tramos: Los segmentos de tubería que conectan dos nudos.
• Mallas: Circuitos cerrados formados por un conjunto de tramos conectados en serie que regresan al punto de partida.

Para que una red hidráulica esté correctamente resuelta y físicamente equilibrada, debe cumplir obligatoriamente dos leyes fundamentales (análogas a las Leyes de Kirchhoff en circuitos eléctricos):

Condición 1: Conservación de la Masa (Ecuación de Continuidad en los Nudos)

La suma algebraica de los caudales que entran y salen de cualquier nudo debe ser exactamente cero. Físicamente, el agua no puede comprimirse ni desaparecer en una intersección.

$$ \sum Q_{entran} = \sum Q_{salen} \implies \sum Q = 0 $$

Condición 2: Conservación de la Energía (Ecuación de Pérdidas en las Mallas)

La suma algebraica de las pérdidas de carga ($h_f$) a lo largo de cualquier circuito cerrado (malla) debe ser cero.

Si partimos de un nudo con una energía $H$, y recorremos la malla sumando las pérdidas a favor del flujo y restando las pérdidas en contra del flujo, al regresar al mismo nudo debemos tener exactamente la misma energía $H$.

$$ \sum h_f = 0 $$

2. Ecuación de Hazen-Williams (Fórmula Empírica)

Para evaluar las pérdidas de carga en los tramos ($\sum h_f$), necesitamos una ecuación que relacione el caudal ($Q$) con la geometría de la tubería. Aunque la ecuación de Darcy-Weisbach es la más rigurosa, en redes de distribución de agua a temperatura ambiente, los ingenieros utilizan casi universalmente la Ecuación de Hazen-Williams por su facilidad matemática (no requiere iterar el factor de fricción).

La pérdida de carga por fricción ($h_f$) en el Sistema Internacional se expresa como:

$$ h_f = \frac{10.67 \cdot L \cdot Q^{1.852}}{C^{1.852} \cdot D^{4.87}} $$

Donde:

• $L$ = Longitud de la tubería ($m$)
• $D$ = Diámetro interior ($m$)
• $Q$ = Caudal volumétrico ($m^3/s$)
• $C$ = Coeficiente de rugosidad de Hazen-Williams (adimensional, depende del material. Ej. PVC $C \approx 150$, Acero usado $C \approx 100$).

Para simplificar el cálculo de redes, agrupamos todos los términos geométricos en una constante de resistencia física del tubo ($K$):

$$ K = \frac{10.67 \cdot L}{C^{1.852} \cdot D^{4.87}} $$

De este modo, la pérdida de carga se reduce a una ecuación general:

$$ h_f = K \cdot Q^{n} $$

(Donde $n = 1.852$ para Hazen-Williams).

3. Consideraciones Básicas y Método de Hardy-Cross

Resolver una red de tuberías implica encontrar el caudal $Q$ en cada tramo. Como las pérdidas dependen de $Q^{1.852}$, el sistema genera ecuaciones algebraicas no lineales que son imposibles de resolver directamente.

En 1936, el ingeniero Hardy Cross desarrolló un método de "relajación iterativa" que permite resolver estas ecuaciones a mano (hoy programado en softwares como EPANET).

Demostración del Método de Hardy-Cross

1. Suposición Inicial: Se asume un caudal inicial ($Q_a$) para cada tramo de la red de tal manera que se cumpla estrictamente la Condición 1 (Continuidad en cada nudo: lo que entra igual a lo que sale).
2. El Error: Como los caudales fueron "adivinados", al calcular las pérdidas en una malla, la suma no será cero, existirá un error.
3. $$ \sum h_f = \sum (K \cdot Q_a^n) \neq 0 $$
4. El Factor de Corrección ($\Delta Q$): El caudal real ($Q_{real}$) en cualquier tubería es igual al caudal asumido ($Q_a$) más una pequeña corrección matemática ($\Delta Q$).
5. $$ Q_{real} = Q_a + \Delta Q $$
6. Para que la Condición 2 se cumpla con los caudales reales:
7. $$ \sum K (Q_a + \Delta Q)^n = 0 $$
8. Expansión de Taylor: Como $\Delta Q$ es pequeño en comparación con $Q_a$, expandimos el binomio usando la serie de Taylor de primer orden y despreciamos los términos de orden superior ($\Delta Q^2$, $\Delta Q^3$):
9. $$ (Q_a + \Delta Q)^n \approx Q_a^n + n \cdot Q_a^{n-1} \cdot \Delta Q $$
10. Sustituyendo esto en la sumatoria:
11. $$ \sum (K \cdot Q_a^n) + \sum (K \cdot n \cdot Q_a^{n-1} \cdot \Delta Q) = 0 $$
12. Como $\Delta Q$ es un valor único de corrección para toda la malla en esa iteración, lo sacamos de la segunda sumatoria y lo despejamos:
13. $$ \Delta Q \sum (n \cdot K \cdot Q_a^{n-1}) = - \sum (K \cdot Q_a^n) $$

Llegamos a la Ecuación Fundamental de Hardy-Cross:

$$ \Delta Q = - \frac{\sum K Q_a^n}{\sum |n \cdot K \cdot Q_a^{n-1}|} = - \frac{\sum h_f}{n \sum \left| \frac{h_f}{Q_a} \right|} $$

(Nota: En el denominador se utiliza el valor absoluto para asegurar que la tasa de cambio de la derivada siempre sume de manera positiva y el método no diverja).

4. Método Iterativo Simplificado (Algoritmo de Cálculo)

Para aplicar este método de manera estructurada (típicamente en una tabla de Excel), se siguen los siguientes pasos:

1. Topología y Signos: Dibuja la red. Define una convención de signos geométrica para cada malla. Por ejemplo: Sentido horario (+) positivo, Sentido antihorario (-) negativo.
2. Distribución Inicial: Asigna caudales $Q_a$ arbitrarios a cada tramo, asegurando que en cada nudo $\sum Q = 0$. Si el flujo en un tramo va en sentido horario, su $Q_a$ es positivo.
3. Cálculo de Parámetros: Calcula la constante de resistencia $K$ para cada tubería.
4. Evaluación de la Malla: Para cada tramo de la malla, calcula:

• La pérdida de carga: $h_f = K \cdot Q_a \cdot |Q_a|^{0.852}$ (Esta forma conserva el signo algebraico del flujo).
• La derivada de la pérdida: $\frac{h_f}{Q_a} = K \cdot |Q_a|^{0.852}$ (Siempre positivo).

1. Cálculo de la Corrección: Suma algebraicamente todos los $h_f$ de la malla. Suma todos los términos $h_f/Q_a$. Aplica la fórmula de $\Delta Q$:
2. $$ \Delta Q = - \frac{\sum h_f}{1.852 \cdot \sum \left| \frac{h_f}{Q_a} \right|} $$
3. Actualización: Corrige los caudales de la malla: $Q_{nuevo} = Q_{anterior} + \Delta Q$.

• Regla de Tramos Compartidos: Si una tubería pertenece a dos mallas contiguas, su caudal se corrige dos veces. Recibe su propia $\Delta Q_{malla1}$ y resta la corrección de la malla vecina $\Delta Q_{malla2}$.

1. Iteración: Repite los pasos 4 a 6 utilizando los $Q_{nuevo}$ hasta que el valor de $\Delta Q$ sea despreciable (ej. $\Delta Q < 0.001 m^3/s$). Cuando esto ocurre, la red está equilibrada y se han encontrado los caudales reales.