UNIVERSIDAD ESTATAL

“PENINSULA DE SANTA ELENA”

FACULTAD CIENCIAS DE LA INGENIERIA
CARRERA DE INGENIERIA EN PETROLEO

TEMA:
CÁLCULO DE POTENCIA (O CABEZA HIDROSTÁTICA), Y
DISEÑO DE TUBERÍAS SIMPLES

MATERIA:
FACILIDADES DE SUPERFICIE I

PROFESOR:
ING. FIDEL CHUCHUCA.

INTEGRANTES:
VALERIA MENDEZ CHONILLO
ANGEL GUALE LAINEZ
FABIAN RIVERA GONZALES
FRANKLIN VARGAS SHIGUANGO

CURSO
6/1
INGENIERIA EN PETROLEO

PERIODO ACADEMICO:
2019- I
 Contenido
OBJETIVO GENERAL: ...................................................................................................................... 4
OBJETIVO ESPECÍFICOS:................................................................................................................. 4
MARCO TEÓRICO ........................................................................................................................... 4
COMPROBACION DE DISEÑO ............................................................................................. 4
CÁLCULO DE LA POTENCIA REQUERIDA....................................................................... 5
DISEÑO DE LA TUBERIA ..................................................................................................... 5
ECUACIONES PARA EL DISEÑO DE TUBERIAS SIMPLES ............................................ 5
COMPROBACIÓN DEL DISEÑO....................................................................................................... 7
CALCULO DE POTENCIA REQUERIDA............................................................................................. 8
METODO DE ITERACION DE UN PUNTO.......................................................................... 9
MÉTODO DE NEWTON-RAPHSON ................................................................................... 11
DISEÑO DE TUBERIAS SIMPLES ................................................................................................... 13
DISEÑO DE TUBERIAS SIMPLES CON ALTAS PERDIDAS MENORES .............................................. 14
LENGUAJE DE PROGRAMACION .................................................................................................. 17
COMPROBACION DE DISEÑO DARCY WEISBACH ...................................................... 17
CALCULO DE POTENCIA DARCY WEISBACH-MÉTODO DE NEWTON RAPHSON 18
COMPROBACIÓN DEL DISEÑO HAZEN WILLIAMS ..................................................... 20
CÁLCULO DE POTENCIA HAZEN WILLIAMS ............................................................... 22
DISEÑO DE TUBERIAS SIMPLES CON ALTAS PERDIDAS MENORES ...................... 23
CONCLUSIONES Y RECOMENDACIONES ...................................................................................... 26
CONCLUSIONES: ................................................................................................................. 26
RECOMENDACIONES: ........................................................................................................ 27
BIBLIOGRAFIA: ............................................................................................................................. 27
ANEXOS ....................................................................................................................................... 27
Análisis de resultados .............................................................................................................. 27
INTRODUCCION
En este proyecto se establecieron las diferentes ecuaciones de resistencia fluida que puede
ser utilizadas para llevar a cabo el diseño de tuberías. Sin embargo, este diseño puede ser
complejo debido a dos razones: en primer término, la ecuación que describe el factor de
fricción de Darcy(f) es una ecuación no explicita, hecho que dio pie al surgimiento de
ecuaciones empíricas más simples. En segundo término, los sistemas de tuberías pueden
ser complicados, dando como resultado variaciones complejas en la caída en la cabeza
piezométrica, debido a la fricción.
Aquí se plantean las metodologías para el diseño de tuberías simples utilizando las
ecuaciones ya planteadas, basadas en el estudio de Prandtl-Von Karman sobre interacción
fluido-pared sólida y en la ecuación de Darcy-Weisbach, considerada la ecuación
físicamente fundamentada para perdidas por fricción en ductos. También se utilizará para
el mismo fin la ecuación de Hazen-Williams, la más representativa de las ecuaciones
empíricas que surgieron cuando se evidencio el hecho de que la ecuación de Colebrook-
white era no explicita.
La tubería simple tiene un diámetro constante y está hecha de un solo material a lo largo
de toda su longitud. La energía que mueve el fluido dentro de ella puede ser de tipo
gravitacional (un embalse o tanque a la entrada) o mecánica (una bomba). En este último
caso, para que la tubería pueda ser considerada como simple, la bomba debe estar
localizada en uno de los extremos. La tubería simple puede tener cualquier tipo de
accesorios que produzcan perdidas menores, incluyendo válvulas para el control del
caudal, forman la base de los algoritmos de diseño de sistemas complejos de tuberías.
OBJETIVO GENERAL:
 Realizar un programa para la comprobación, cálculo de potencia (o cabeza
hidrostática), y diseño de tuberías simples.

OBJETIVO ESPECÍFICOS:
 El programa a realizarse debe permitir entrada de datos
 Permitir escoger entre comprobación, cálculo de potencia o diseño de tubería.
 Permitir escoger entre la resolución por la ecuación de Darcy-Weisbach o Hazen-
Williams (tomando en cuenta las restricciones de la ecuación)
 Utilizar el método de Newton-Raphson
 Utilizar ventanas.

MARCO TEÓRICO
TIPOS DE PROBLEMAS EN HIDRÁULICA DE DUCTOS A PRESIÓN
Los problemas en la hidráulica de tubería simples se pueden clasificar de acuerdo con la
variable desconocida en el problema. Las variables involucradas en problemas de tuberías
simples son las siguientes:
 Variables relacionadas con la tubería en sí: diámetro de la tubería(d), longitud de
la tubería(l) y rugosidad absoluta de la tubería (ks).
 Variables relacionadas con el fluido: densidad del fluido (𝜌) y viscosidad
dinámica del fluido (𝜇).
 Variables relacionadas con el esquema del sistema: coeficientes de perdidas
menores de todos los accesorios necesarios, incluyendo válvulas (Σ𝑘𝑚 ).
 Variables relacionadas con la energía impulsadora del fluido: cabeza entre el
embalse de entrada y la salida (H) o potencia de la bomba(P).
 Otras variables: aceleración de la gravedad(g) y caudal o velocidad media en la
tubería (Q o v).
Las variables que pueden ser desconocidas son el caudal, el diámetro de la tubería o la
potencia de la bomba que es necesario instalar para mover un determinado caudal.
Teniendo en cuenta la variable desconocida, los tipos de problemas son: comprobación
de diseño, cálculo de potencia requerida y diseño de la tubería.

COMPROBACION DE DISEÑO
En este tipo de problemas la tubería existe, lo cual implica que se conocen su longitud, su
diámetro y su rugosidad absoluta (se conoce el material9, al igual que todos los accesorios
y sus coeficientes de perdidas menores. Además, se conoce la energía impulsadora, ya
sea una cabeza gravitacional o una bomba y las propiedades del fluido (densidad y
viscosidad dinámica). La incógnita es el caudal que pasa por la tubería, problema típico
en el diseño de distribución de agua potable o redes de riego, en las cuales se hace un pre
dimensionamiento de los diámetros.
Variables conocidas d, ks, H (o P), Σ𝑘𝑚 , 𝜌, 𝜇, g, l
Incógnita Q (o v)

CÁLCULO DE LA POTENCIA REQUERIDA

Para estos problemas se conoce el caudal demandado, la tubería (se conocen tanto la
longitud como el diámetro, la rugosidad absoluta, los coeficientes de pérdidas menores
de los accesorios) y las propiedades del fluido (densidad y viscosidad dinámica), Se tiene
que calcular la potencia necesaria (bomba o diferencia de nivel) para mover dicho caudal
a través de la tubería. Este problema se plantea cuando se quiere utilizar una tubería
existente para mover, un cierto caudal demandado y se desea conocer la bomba que debe
ser colocada o la diferencia de nivel entre la entrada y la salida de la tubería.
Variables conocidas Incógnita
d, ks, Q (o v), Σ𝑘𝑚 , 𝜌, 𝜇, g, l H (o P)

DISEÑO DE LA TUBERIA
En este caso se conoce el caudal demandado, la potencia disponible (bomba o diferencia
de nivel generada por un embalse o tanque), algunas de las características de la tubería
(longitud y accesorios necesarios con sus correspondientes coeficientes de perdidas
menores) y las propiedades del fluido (densidad y viscosidad dinámica). Se desconoce el
diámetro necesario para permitir el paso del caudal demandado. Por lo general en cuanto
al material de la tubería se tiene dos o tres alternativas, lo cual implica que se conoce la
rugosidad absoluta de la tubería, pero es necesario hacer tantos diseños como diferentes
materiales se encuentren disponibles en el sitio del proyecto con el fin de escoger la mejor
alternativa.
Variables conocidas Incógnita
ks, Q (o v), H (o P), Σ𝑘𝑚 , 𝜌, 𝜇, g, l d

ECUACIONES PARA EL DISEÑO DE TUBERIAS SIMPLES

Mediante el uso de las ecuaciones de Colebrook-White (ecuación 1.67) conjuntamente
con la ecuación de Darcy-Weisbach (ecuación 1.36) se puede desarrollar el siguiente
procedimiento:
A partir de la figura, que corresponde a una tubería simple alimentada por un tanque de
nivel constante, se puede plantear la ecuación de Bernoulli entre los puntos 1 y 2, para
llegar a la siguiente ecuación:
𝑉22 𝑃2
ℎ1 + 𝑧1 = + 𝑧2 + + ℎ𝑓 + ∑ ℎ𝑚
2𝑔 𝛾
Si el punto 2 corresponde a la salida, de la tubería, caso muy común en este tipo de
análisis, ocurren tres cosas en la ecuación anterior. En primer lugar, el término de la
cabeza de velocidad desaparece, ya que se supone que en la salida el flujo pierde su
velocidad. En segundo lugar, y como consecuencia de lo anterior, el término de sumatoria
de pérdidas menores debe incluir un término de pérdidas por salida. Por último, la presión
en la salida es igual a la atmosférica (presión manométrica nula); luego:
𝑣2 2 𝜌2
= 0; =0
2𝑔 𝜌𝑔
Es claro que para el punto 1, ℎ1 + 𝑧1 = 𝐻, donde H es la altura del nivel de la superficie
del tanque con respecto al datum. Luego:
𝐻 = 𝑧2 + ℎ𝑓 + Σℎ𝑚

De esta última expresión se puede obtener la siguiente ecuación, la cual describe las
perdidas por fricción en función de las otras variables:
2 𝑉2
ℎ𝑓 = 𝐻 − 𝑧2 − ∑ 𝑘𝑚 2𝑔 (2.1)

Si el punto 2 no corresponde a la salida de la tubería, la ecuación (2.1) tendría que incluir

un término de cabeza de velocidad y la sumatoria de las perdidas menores no deberá
incluir un término de pérdidas por salida.
Adicionalmente debería incluir la cabeza de presión en el punto 2. Mediante la ecuación
de Darcy-Weisbach (ecuación 1.36), que también debería incluir las perdidas por fricción,
se puede despejar el factor de fricción f.
ℎ𝑓 . 𝑑. 2. 𝑔
𝑓=
𝑙𝑉 2
Lo cual implica que:

√2𝑔𝑑ℎ𝑓
√𝑓 =
√𝑙 𝑣
Y que:
1 𝑣 √𝑙
= (2.2)
√𝑓 √2𝑔𝑑ℎ𝑓

Al reemplazar la ecuación 2.2 en la ecuación 1.67 se obtiene:

1 𝑘 2.51
= −2 log10 (3.7𝑠𝑑 + ) (1.67)
√𝑓 𝑅𝑒 √𝑓

𝑣 √𝑙 𝑘𝑠 2.51𝜈√𝑙
= −2 log10 ( + )
√2𝑔𝑑ℎ𝑓 3.7 𝑑 𝑅𝑒 √2𝑔𝑑ℎ𝑓

El número de Reynolds en esta última ecuación puede reemplazarse por:

𝑣𝑑
𝑅𝑒 =
𝜈
Por consiguiente, se obtiene la siguiente expresión:

𝑣 √𝑙 𝑘𝑠 2.51𝜈𝑣√𝑙
= −2 log10 ( + )
√2𝑔𝑑ℎ𝑓 3.7 𝑑 𝑣 𝑑 √2𝑔𝑑ℎ𝑓

Finalmente, si se despeja la velocidad se encuentra una ecuación explicita para esa

variable:
−2√2𝑔𝑑ℎ𝑓 𝑘 2.51𝜇√𝑙
𝑉= 𝑙𝑜𝑔 (3.7𝑠𝑑 + 𝑑𝑔 ) (2.3)
√𝑙 √2𝑔𝑑ℎ𝑓

Esta última ecuación es la base para la solución de los tres tipos de problemas
relacionados con tuberías simples.

COMPROBACIÓN DEL DISEÑO

En este caso se conocen todas las características de la tubería: la longitud, el diámetro, la
rugosidad absoluta y los diferentes coeficientes de perdidas menores de cada uno de los
accesorios. También se conocen las propiedades del fluido, es decir, la densidad y la
viscosidad dinámica. La incógnita es la velocidad y por consiguiente, el caudal que pasa
por la tubería. Obviamente se deben conocer la cabeza disponible o la potencia y la
eficiencia de la bomba. Para resolver estos problemas se debe seguir el procedimiento
establecido en el diagrama de flujo 1.
CALCULO DE POTENCIA REQUERIDA
Todas las características de la tubería son conocidas, al igual que las del fluido. Se sabe
un caudal demandado y se pregunta por la potencia requerida, ya sea de origen
gravitacional o mecánico.
Para resolver los problemas de potencia requerida es necesario utilizar un método
numérico con el fin de averiguar el valor del factor de fricción f de Darcy en la ecuación
no explicita de Colebrook-White:
1 𝑘 2.51
= −2 log10 (3.7𝑠𝑑 + ) (1.67)
√𝑓 𝑅𝑒 √𝑓

A fin de resolver la ecuación 1.67, para f existen muchos métodos numéricos, de los
cuales dos de ellos se explican a continuación.

METODO DE ITERACION DE UN PUNTO

Para que este método pueda ser aplicado, la función no explicita debe ser de la siguiente
forma:
𝑥 = 𝑔(𝑥)
El algoritmo se desarrolla de tal manera que el valor arrojado por la función g(x) en la
iteración i se utilice como argumento x en la iteración i + 1. Si la aproximación inicial es
suficientemente cercana a la raíz α, entonces el método será convergente si α es una raíz
simple. Para el método de iteración de un punto, esto es aun cierto si la raíz no es simple,
pero la convergencia es más lenta. en ese caso. Sin embargo, en general no es posible
demostrar que una función de iteración de punto converge a una raíz múltiple aun si las
condiciones iniciales están arbitrariamente cerca de la raíz. Afortunadamente para el caso
de la ecuación de Colebrook-White, la raíz siempre es única.
En el caso de la ecuación de Colebrook-White, el método converge en 8 ó 10
aproximaciones y es muy sensible al valor inicial de f que· se suponga (semilla). En la
figura 2.5 se esquematiza el proceso de convergencia. En ésta se representa el caso de una
tubería con rugosidad relativa (k/ci) de 0.0001 y con números de Reynolds de 20000 y
200000; la línea a 45° representa la función g(x) = x. En ambos casos, la semilla para el
factor de fricción de Darcy fue f = 0.001.
El diagrama de flujo 2a esquematiza el proceso seguido por el método de Newton.
MÉTODO DE NEWTON-RAPHSON
Este método es una aceleración del método anterior, hecho que lo hace más conveniente;
por lo general se requieren sólo tres iteraciones. Sin embargo, la función g(x) debe
cumplir las siguientes condiciones especiales para que exista convergencia:
 La primera condición es que exista un intervalo l=(a,b), de modo que para todo x
perteneciente a l, la función g(x) esté definida y pertenezca a l, lo cual significa
que g(x) se aplica a sí misma. En el caso de la ecuación de Colebrook-White, para
que la función no estuviera definida se necesitaría que el logaritmo no estuviera
definido, caso imposible, ya que todos los términos dentro de la función logaritmo
son positivos.
 La segunda condición es que la función de iteración g(x) sea continua en l.
Nuevamente la función logaritmo cumple este requisito.
 La tercera condición de convergencia es que g(x) sea diferenciable en l y que la
pendiente de g(x) sea siempre menor que 1 y mayor que –l. La función
mencionada es en efecto diferenciable y su pendiente es siempre mayor que -1,
llegando a valores extremos (Re=2000, ks/d=0.00001) a ser -0.0289.
Para números de Reynolds grandes la pendiente de g(x) tiende a cero, lo cual es
menor que 1.
Teniendo en cuenta las anteriores condiciones en la iteración i+1 se calcula con base en
la aproximación de la iteración i de acuerdo con la siguiente ecuación:
𝑋 = 𝑔(𝑥)
Por consiguiente:
𝑔(𝑥) − 𝑥 = 0
El valor de la aproximación a la raíz de la ecuación en la iteración i+1 se calcula con base
en la aproximación de la iteración i de acuerdo con la siguiente ecuación:
𝑔(𝑥𝑖 ) − 𝑥𝑖
𝑥𝑖+1 = 𝑥𝑖 −
𝑔′ (𝑥𝑖 ) − 1
Donde 𝑔′ (𝑥) es la derivada de la función. Para el caso de la ecuación de Colebrook-White
la variable x es el inverso de la raíz cuadrada del factor de fricción de Darcy:
1
𝑥=
√𝑓
Por consiguiente:
𝑘𝑠 2.51𝑥𝑖
𝑔(𝑥) = −2 log10 ( + )
3.7 𝑑 𝑅𝑒
2.51
′ (𝑥)
2 𝑅𝑒
𝑔 =− ( )
ln 10 𝑘 𝑠 2.51𝑥
+ 𝑅𝑒 𝑖
3.7 𝑑
En la figura 2.6 se esquematiza el proceso de convergencia de este método. En este caso
se utilizó una tubería con rugosidad relativa (ks/d) de 0.0001 y un número de Reynolds
de 20000. El valor semilla para f fue 0.002. A pesar de que este valor estaba bastante lejos
del valor real (f = 0.0261), el método convergió muy rápidamente; sus ventajas sobre el
método anterior saltan a la vista.
DISEÑO DE TUBERIAS SIMPLES
El proceso de diseño es bastante simple porque la ecuación (2.3) es explicita para la
velocidad. Dicho proceso se esquematiza en el diagrama de flujo 4. Este diagrama puede
utilizarse tanto para tuberías de diámetros comerciales como tuberías de acero (diámetro
a la medida).
−2√2𝑔𝑑ℎ𝑓 𝑘 2.51𝜇√𝑙
𝑉= 𝑙𝑜𝑔 (3.7𝑠𝑑 + 𝑑𝑔 ) (2.3)
√𝑙 √2𝑔𝑑ℎ𝑓

Es necesario aclarar que para que el proceso de diseño converja se deben tener en cuenta
dos restricciones importantes:
 El primer diámetro supuesto tiene que ser menor que el diámetro final del diseño.
Se sugiere empezar siempre con un diámetro muy pequeño, lo cual implica un
mayor número de iteraciones, pero asegura un correcto resultado.
 La suma de las pérdidas menores debe ser inferior al 30% de las pérdidas por
fricción. Esta condición es cierta en la mayoría de diseños convencionales de
tuberías. Más adelante se tratará el caso de sistemas de tuberías en los cuales las
pérdidas menores causadas por los accesorios son mayores al 30% de las pérdidas
por fricción. Para que el diagrama de flujo 4 converja se debe cumplir la siguiente
ecuación:
Σℎ𝑚 ≤ 0.3ℎ𝑓 (2.4)

Esta última restricción en la práctica resulta irrelevante, ya que en la gran mayoría de los
sistemas de tuberías tal condición se cumple con facilidad. Para diseñar un sistema con
altas pérdidas menores, como es el caso de la tubería de succión de una bomba, se debe
seguir un algoritmo diferente.
DISEÑO DE TUBERIAS SIMPLES CON ALTAS PERDIDAS MENORES
Este proceso consiste en definir una “velocidad” de perdida que en esencia es la velocidad
que haría igualar la sumatoria de las perdidas menores y la cabeza disponible. Si se
utilizara esta definición, las ecuaciones serian:
Σℎ𝑚 = 𝐻
𝑣2
Σℎ𝑚 = (Σ𝑘𝑚 )
2𝑔
Mediante las dos ecuaciones anteriores se obtiene el siguiente resultado para la “velocidad
de pérdida":
𝑣𝑝 2 𝐻
=
2𝑔 Σ𝑘𝑚
Si se despeja Vp en esta última ecuación se llega a:
2𝑔𝐻
𝑣𝑝 2 =
Σ𝑘𝑚
2𝑔𝐻
𝑣𝑝 = √Σ𝑘 (2.5)
𝑚

Si en alguna iteración la Vi es mayor que la Vp, esto quiere decir que la velocidad Vi
implica unas pérdidas menores superiores a la cabeza disponible, lo cual es físicamente
imposible. De suceder así, se debe limitar la cabeza disponible para ser perdida por
fricción, dentro del procedimiento de diseño.
El procedimiento se esquematiza en el diagrama de flujo 5. el cual es más general que el
diagrama de flujo 4, ya que también sirve para el caso de tuberías con pérdidas menores
bajas. Una vez se ha calculado la primera velocidad de pérdida, en las demás iteraciones
esta velocidad se calcula de acuerdo con la siguiente ecuación:

2𝑔(𝐻−ℎ𝑓 )
𝑣𝑝 = √ (2.6)
Σ𝑘𝑚
LENGUAJE DE PROGRAMACION
Este algoritmo se realizó con el fin de optimizar tiempo al realizar los cálculos y obtener
mejores resultados para el diseño de tuberías simples. El lenguaje de programación
utilizado fue Python.
COMPROBACION DE DISEÑO DARCY WEISBACH
import math
###### pedir datos
#### input pide los datos
d=0.293
k_s=1.5*(10**(-6))
u=1.005*(10**(-3))
S_km=11.8
z_2=0
l=730
H=43.5
rho=998.2
h_f=H
g=9.81

##### funciones de lo que se calcula a cada rato

# math.sqrt() calcula la raiz
# math.log10() calcula el logaritmo
def calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho):
visc_cin=u/rho
a=math.sqrt(2*g*d*h_f)
b=math.sqrt(l)
c=math.log10((k_s/(3.7*d))+(2.51*visc_cin*b/(d*a)))
v=-2*a*c/b
return v

def calcular_h_f(v, H, z_2, S_km):

h_f_i=H-z_2-(S_km*(v**2)/(2*g))
return h_f_i

def calcular_h_m(H,h_i):
h_m=H-h_i
return h_m

def calcuar_Q(v_i,d):
Q=v_i*((math.pi/4)*(d**(2)))
return Q

###### hacer el calculo inicial

v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)

h_i=calcular_h_f(v_i, H, z_2, S_km )
error=abs(h_f-h_i)
print('velocidad: ', v_i)
print('h_f: ', h_f)
print('h_i: ', h_i)
print('error: ', error)
print('***********************')

##### se hacen las iteraciones

while(error>0.001):
h_f=h_i
v_i=calcular_velocidad(h_i,d,k_s, u, S_km, z_2, l, rho)
h_i=calcular_h_f(v_i, H, z_2, S_km )
error=abs(h_f-h_i)
print('velocidad: ', v_i)
print('h_f: ', h_f)
print('h_i: ', h_i)
print('error: ', error)
print('***********************')
####Calculo final con datos de ultima iteracion####
print('******Calculo Final******')
h_m=calcular_h_m(H,h_i)
print('h_m: ', h_m)

print('h_i: ', h_i)

Q=calcuar_Q(v_i,d)
print('Q: ', Q)

CALCULO DE POTENCIA DARCY WEISBACH-MÉTODO DE NEWTON

RAPHSON
import math
###### pedir datos
#### input pide los datos
d=0.1522
k_s=1.5*(10**(-6))
u=1.005*(10**(-3))
S_km=9.4
z_2=16
l=970
rho=999.1
Q=0.042
#Q=0.00000042
visc_cin=1.14*(10**(-6))
g=9.81

def calcular_Re(Q, d):

area=(math.pi*d**2)/4
vel=Q/(area)
 Re=vel*d/visc_cin
return Re

def calcular_g_de_x(k_s, d, Re, x_i):

g=(-2*math.log10((k_s/(3.7*d)+2.51*x_i/Re)))
return g

def calcular_g_prima_de_x(k_s, d, Re, x_i):

a=2.51/Re
b=k_s/(3.7*d)
c=2.51*x_i/Re
g_prima=(-2/math.log(10))*(a/(b+c))
return g_prima

def calcular_x_f(x_i, g_i, g_prima_i):

x_f=x_i-(g_i-x_i)/(g_prima_i-1)
return x_f

Re=calcular_Re(Q, d)
if Re<=2200:
f=64/Re
print('f: ', f)
else:
###### primera iteracion
f_i=0.001
x_i=1/math.sqrt(f_i)
g_i=calcular_g_de_x(k_s, d, Re, x_i)
g_prima_i=calcular_g_prima_de_x(k_s, d, Re, x_i)
x_f=calcular_x_f(x_i, g_i, g_prima_i)
f=1/(x_f**2)
error=abs(x_f-x_i)
print('g_i: ', g_i)
print('x_i: ', x_i)
print('x_f: ', x_f)
print('f: ', f)
print('error: ', error)
print('***********************')

###### iteraciones newton raphson

while(error>0.000000001):
x_i=x_f
g_i=calcular_g_de_x(k_s, d, Re, x_i)
g_prima_i=calcular_g_prima_de_x(k_s, d, Re, x_i)
x_f=calcular_x_f(x_i, g_i, g_prima_i)
f=1/(x_f**2)
error=abs(x_f-x_i)
print('g_i: ', g_i)
print('x_i: ', x_i)
 print('x_f: ', x_f)
print('f: ', f)
print('error: ', error)
print('***********************')

print('*******Calculo final****')
######### calcular h_f
area=(math.pi*d**2)/4
vel=Q/(area)
h_f=f*l*vel**2/(d*2*g)
print('h_f: ', h_f)

####### calcular suma hm

S_hm=S_km*vel**2/(2*g)
print('S_hm: ', S_hm)

######## calcular H
H=z_2+h_f+S_hm
print('H: ', H)

#### calcular potencia de la bomba

Pot=rho*Q*g*H
print('Potencia: ', Pot)

COMPROBACIÓN DEL DISEÑO HAZEN WILLIAMS

import math
###### pedir datos
#### input pide los datos
d=0.2032
S_km=6.4
l=240
H=37
h_f=H
g=9.81
C_hw=150

def velocidad(C_hw,d,h_i,l):
V=0.849*C_hw*(d**(0.63))*(h_i**(0.54)/(2.395*(l**(0.54))))
return V

def Area(d):
A=(math.pi/4)*(d**(2))
return A
def Caudal(v,A):
Q=v*A
return Q

def perdidas_menores(S_km,v,g):
S_hm=S_km*(v**(2)/(2*g))
return S_hm

def h_F_siguiente(H,S_hm):
h_ff=H-S_hm
return h_ff

def perdidas_por_friccion(l,v,C_hw,d):
h_f=6.824*l*(v**(1.851))/((C_hw**1.851)*(d**(1.167)))
return h_f

if S_km==0:
v=velocidad(C_hw,d,H,l)
print('vel: ', v)

A= Area(d)
print('A: ', A)

Q=Caudal(v,A)
print ('Q: ', Q)
else:
h_f=H
v=velocidad(C_hw,d,h_f,l)
S_hm=S_km*(v**(2)/(2*g))
h_ff=H-S_hm
Q=Caudal(v,Area(d))
error=abs(h_ff-h_f)
print('h_f:', h_f)
print('h_ff:', h_ff)
print('velocidad', v)
print('Q: ', Q)
print('*****************')

while(error>0.0001):
h_f=h_ff
v=velocidad(C_hw,d,h_f,l)
S_hm=S_km*(v**(2)/(2*g))
h_ff=H-S_hm
Q=Caudal(v,Area(d))
error=abs(h_ff-h_f)
print('h_f:', h_f)
print('h_ff:', h_ff)
print('velocidad', v)
print('Q: ', Q)
 print('S_hm: ', S_hm)
print('*****************')

CÁLCULO DE POTENCIA HAZEN WILLIAMS

import math
###### pedir datos
#### input pide los datos
d=0.1524
S_km=7.2
l=465
Q=0.102
Z_2=22
n=0.85
rho=1000
g=9.81
C_hw=120

def velocidad(Q,d):
v=(Q)/((math.pi/4)*(d**2))
return v

def perdidas_menores(S_km,v,g):
S_hm=S_km*(v**(2)/(2*g))
return S_hm

def perdidas_por_friccion(l,v,C_hw,d):
h_f=(6.824*l*(v**(1.851)))/((C_hw**1.851)*(d**(1.167)))
return h_f

def cabeza_total_requerida(Z_2,h_f,S_hm):
H=Z_2+h_f+S_hm
return H

def potencia(rho,n,Q,H,g):
P= (1/n)*rho*Q*g*H
return P

v= velocidad(Q,d)
h_f= perdidas_por_friccion(l,v,C_hw,d)
S_hm= perdidas_menores(S_km,v,g)
H= cabeza_total_requerida(Z_2,h_f,S_hm)
P= potencia(rho,n,Q,H,g)

print('S_hm: ', S_hm)

print('h_f: ', h_f)
print('H:', H)
print('P: ', P)
DISEÑO DE TUBERIAS SIMPLES CON ALTAS PERDIDAS MENORES
import math
###### pedir datos
#### input pide los datos
k_s=0.00015
u=1.14*(10**(-3))
delta_d=2*0.0254
delta_h=0.005
z_2=0
l=17
H=2.2
rho=999.1
Q_d=0.12
visc_cin=1.17*(10**(-6))
g=9.81
S_km=4.3
d=0.1016

# math.sqrt() calcula la raiz

# math.log10() calcula el logaritmo
def calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho):
visc_cin=u/rho
#print('hf', h_f)
a=math.sqrt(2*g*d*h_f)
b=math.sqrt(l)
#print('d', d)
#print('a', a)
c=math.log10((k_s/(3.7*d))+(2.51*visc_cin*b/(d*a)))
v=-2*a*c/b
return v

def Caudal(v,d):
Q_i=(v)*((math.pi/4)*(d**(2)))
return Q_i

def S_perdidas_menores(S_km,v,g):
S_hm=(S_km*((v**(2))/(2*g)))
return S_hm

def calcular_h_f(v, H, z_2, S_km):

h_f_i=H-z_2-(S_km*(v**2)/(2*g))
return h_f_i

def Velocidad_de_perdida(g,H,h_f,S_km):
v_p=math.sqrt((2*g*(H-h_f))/(S_km))
return v_p

def cabeza_disponible(h_f,S_hm):
 H_t=(h_f+S_hm)
return H_t

def ParteA_algorithmo(h_f,k_s, d, delta_d, Q_d, u, S_km, l, rho, g, H):

#h_f_1=h_f
h_f=0.5
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
while Q_i<=Q_d:
print('Parte A: Qi menor Qd')
d=d+delta_d
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
v_p=Velocidad_de_perdida(g,H,h_f,S_km)
error_v=abs(v_i-v_p)
while error_v>=0.01:
print('***** error v mayor que 0.01')
print('Q_i', Q_i)
print('v_i', v_i)
print('v_p', v_p)
print('h_f', h_f)
if v_i<v_p:
h_f=h_f+delta_h
else:
h_f=h_f-delta_h
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
while Q_i<=Q_d:
print('Parte A: Qi menor Qd 2')
d=d+delta_d
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
v_p=Velocidad_de_perdida(g,H,h_f,S_km)
error_v=abs(v_i-v_p)
return d, v_p, v_i, h_f

#inicio
h_f=H-z_2
#### aqui debe volver el loop
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
while Q_i<Q_d:
print('*******Q_i menor que Q_d')
d=d+delta_d
print('d:', d)
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
print('v_i', v_i)
 print('Q_i', Q_i)
v_p=Velocidad_de_perdida(g,H,h_f,S_km)
print('v_p', v_p)
if v_i<v_p:
h_f_1=h_f
h_f=calcular_h_f(v_i, H, z_2, S_km)
error=abs(h_f_1-h_f)
while error>0.01:
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
h_f_1=h_f
h_f=calcular_h_f(v_i, H, z_2, S_km)
error=abs(h_f_1-h_f)
print('v_i', v_i)
print('Q_i', Q_i)
print('h_f', h_f)
print('********')
while Q_i<Q_d:
d=d+delta_d
h_f=H-z_2
###### empieza de nuevo el algoritmo
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
while Q_i<Q_d:
print('Q_i menor que Q_d')
d=d+delta_d
print('d:', d)
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
print('v_i', v_i)
print('Q_i', Q_i)
print('********')
v_p=Velocidad_de_perdida(g,H,h_f,S_km)
if v_i<v_p:
h_f_1=h_f
h_f=calcular_h_f(v_i, H, z_2, S_km)
error=abs(h_f_1-h_f)
while error>0.01:
v_i=calcular_velocidad(h_f,d,k_s, u, S_km, z_2, l, rho)
Q_i=Caudal(v_i, d)
h_f_1=h_f
h_f=calcular_h_f(v_i, H, z_2, S_km)
print('v_i', v_i)
print('Q_i', Q_i)
print('********')
error=abs(h_f_1-h_f)
else:
print('v_i mayor que v_p')
 d, v_p, v_i, h_f=ParteA_algorithmo(h_f,k_s, d, delta_d, Q_d, u,
S_km, l, rho, g, H)
print('d: ', d)
print('*****final ***+*')
print('v_i', v_i)
print('v_p', v_p)
break
###
print('Fin d:', d)
print('*****final ***+*')
print('v_i', v_i)
print('v_p', v_p)

else:
print('v_i mayor que v_p')
d, v_p, v_i, h_f=ParteA_algorithmo(h_f,k_s, d, delta_d, Q_d, u, S_km,
l, rho, g, H)
print('d: ', d)
print('*****final ***+*')
print('v_i', v_i)
print('v_p', v_p)

##### fin del algoritmo

print('*****final ***')
print('v_i', v_i)
print('v_p', v_p)
S_hm=S_perdidas_menores(S_km,v_i,g)
print('S_hm', S_hm)
H_t=cabeza_disponible(h_f,S_hm)
print('H_t: ', H_t)

CONCLUSIONES Y RECOMENDACIONES
CONCLUSIONES:
 El algoritmo de Python contiene un lenguaje de programación sencillo pero eficaz
que logra determinar el caudal, la cabeza hidrostática y el diámetro de la tubería
para un buen diseño de tuberías simples.
 Este algoritmo se realizó con el fin de optimizar tiempo al realizar los cálculos de
manera inmediata.
 El algoritmo permite elegir el cálculo por las ecuaciones de Darcy-weisbach o
Hazen Williams.
 Facilita al usuario un cálculo eficaz y eficiente de forma inmediata.
RECOMENDACIONES:
 Se recomienda análisis para la toma de datos previos.
 Se debe escoger entre comprobación de diseño, cálculo de potencia, y diseño de
tubería con altas perdidas menores, asi mismo por las ecuaciones de Darcy
weisbach y Hazzen Williams, para el ingreso de los datos.
 Se recomienda ingresar los datos completos requeridos para efectuar el cálculo
del método elegido.

BIBLIOGRAFIA:
Hidráulica de tuberías- Juan G. Saldarriaga V.

ANEXOS
Análisis de resultados
1.- Ventana de Inicio del programa de Diseño de tuberías simples.
Permite escoger entre las opciones de Comprobación de diseño, cálculo de potencia, y
diseño de tuberías simples, asi mismo permite elegir entre las resoluciones por las
ecuaciones de Darcy Weisbach y Hazen Williams.

2.- Comprobación de diseño Darcy Weisbach

El programa permite la entrada de datos, una vez ingresados los datos requeridos, el
programa calcula y nos muestra los resultados con los datos de la última iteración.
3.- Cálculo de potencia Darcy Weisbach - Método de Newton Raphson.
El programa permite la entrada de datos, una vez ingresados los datos requeridos, el
programa calcula el factor de friccion de darcy y nos muestra los resultados con los datos
de la última iteración.

4.- Comprobación de diseño Hazen Williams

En este método nos muestra los resultados de comprobación de diseño, con las ecuaciones
de Hazen Williams, y calculándola velocidad, el caudal requerido para el diseño con cada
una de sus pérdidas por fricción y perdidas menores .
5.- Cálculo de potencia Hazen Williams.
Esta ventana nos muestra el calculo de un método sencillo como es el calculo de potencia
por el método de Hazen Williams. Una vez ingresando todos los datos requeridos el
programa nos muestra los resultados obtenidos.
6.-Diseño de tuberías simples
Esta ventana muestra el cálculo de diseño de tuberías simple con altas perdidas menores,
mostrando a su vez los resultados obtenidos con la última iteración, siendo estos el
diámetro de la tubería, perdidas por fricción, perdidas menores, velocidad, caudal y
velocidad de pérdida, además del cálculo de la cabeza hidrostática.