Posteado por RubensaiD el 3rd Septiembre 2010 // 716 lecturas

2

MN217 Mecanica de Fluidos II. Calculo

del Factor de Friccion f mediante la

Ecuacion de Colebrook para Hp 50g

Seguramente reconocerán esta ecuación, pues se trata de la Ecuación de Colebrook, la

cual es una ecuación semi-empírica que nos permite calcular el factor de fricción f

(perdidas primerias para flujo en tuberías).

La parte bonita de esta ecuación es que nos permite “prescindir” del Diagrama de

Moody para calcular este factor. Por otro lado, la necesidad de realizar varias iteraciones

al utilizar esta ecuación es su parte negativa.

Bueno, aplicando lo aprendido en los cursos de Programación Orientada a Objetos

(Estructuras de Control de Buble) y Métodos Numéricos (Teoría de Errores) he diseñado

este pequeño programa que nos permite encontrar el valor de f con 5 cifras decimales

exactas, esto es, una cota (o tolerancia) de .

Para utilizar este programa deberemos entregar los siguientes datos: Números de Reynolds

(Re), Rugosidad Relativa y f asumido (este valor lo podemos sacar del diagrama de

Moody con una inspección poco minuciosas, pues el programa se encargara de

arreglarlo).

Ahora sí, aquí el código:

1

2

3

4

5

6

7

8

9

10

11

12

13

@Creado por Rubensaid

«

"Ec. de Colebrook"

{

{"Re= " "Numero de Reynold"}

{"“= " "Rugosidad Relativa"}

}

{ 1 0 }

{ }

{ }

INFORM

IF 1 ==

THEN

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

OBJ• DROP

0 0 0 0 0 �→ Re er f fcalc error cota fca

«

0.1 'f' STO

0.5 10 -5 ^ * 'cota' STO

cota 1 + 'error' STO

WHILE error cota >

REPEAT

er 3.71 / 2.51 Re f 0.5 ^ * / + LOG -2 * -2 ^ 'fcalc' STO

fcalc f - ABS 'error' STO

fcalc 'f' STO

END

f "f" �→TAG

5 FIX

»

END

»