MATLAB Métodos Numéricos 3 Nota
MATLAB Métodos Numéricos 3 Nota
MATLAB Métodos Numéricos 3 Nota
DE CAJAMARCA
FACULTAD DE INGENIERA
ESCUELA ACADMICO PROFECIONAL DE INGENIERA
CIVIL
APLICACIN DE MATLAB
TEMA: A LA
INGENIERA
TUBERAS
SIMPLES
Integrantes:
CHILN JULCA, Royer Alfredo.
HUAMN VILLANUEVA, Thania Celfa.
NARVAJO GOICOCHEA, Katerine.
TACILLA ROJAS, Vctor Manuel.
VILLANUEVA TERN, Harold Edwin Ronald.
DE ENERO DE 2017
MTODOS NUMRICOS
NDICE
NDICE................................................................................................................ 2
INTRODUCCIN................................................................................................... 3
OBJETIVOS.......................................................................................................... 4
MARCO TERICO.......................................................................................... 5
1.
1.
2.
COMPROBACIN DE DISEO......................................................................7
5.
6.
DISEO DE IA TUBERA............................................................................... 8
7.
COMPROBACIN DE DISEO......................................................................9
2.
While:....................................................................................................... 10
For:........................................................................................................... 10
Ejercicios..................................................................................................... 11
3.
1.
Ejercicio de caudal................................................................................... 15
Ejercicio de dimetro................................................................................ 18
CONCLUSIONES................................................................................................ 21
BIBLIOGRAFA................................................................................................... 22
ANEXOS............................................................................................................ 23
MTODOS NUMRICOS
INTRODUCCIN
MTODOS NUMRICOS
OBJETIVOS
Objetivos principales:
Objetivos especficos:
Los resultados del caudal deben ser idnticos a los resultados de un desarrollo
normal por iteraciones.
MTODOS NUMRICOS
1.
MARCO TERICO
Variables relacionadas con la energa impulsora del fluido: Cabeza entre el embalse de
entrada y la salida (H) o potencia de la bomba (P),
Las variables que pueden ser desconocidas son el caudal, el dimetro de la tubera 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: Comprobacin de diseo, clculo de potencia
requerida y diseo en s de la tubera.
1.
Figura 1.1 Lneas de energa total y de gradiente hidrulico (o lnea piezomtrica) para un sistema
de tuberas a presin conectando dos tanques y utilizando una bomba que produce una cabeza HB
hmi = prdida menor i.
hf= prdida por friccin en la tubera i,
2.
Figura 1.2 Lneas de energa total y de gradiente hidrulico (o lnea piezomtrica) para un sistema
de tubera a presin que une dos tanques. En este caso, la potencia requerida para mover el fluido
es suministrada por la gravedad.
3.
Figura 1.3 Lneas de energa total y de gradiente hidrulico (o lnea piezomtrica) para un sistema
de tubera a presin en una central de generacin hidroelctrica. HT representa la cabeza
absorbida por la turbina.
4.
COMPROBACIN DE DISEO
Se conocen su longitud, su dimetro y su rugosidad absoluta (se conoce el material), al igual que
todos los accesorios y sus coeficientes de prdidas menores. Adems, se conoce la energa
impulsora, ya sea una cabeza gravitacional o una bomba y las propiedades del fluido (densidad y
viscosidad dinmica). La incgnita es el caudal que pasa por la tubera, problema tpico en el
diseo de redes de distribucin de agua potable o redes de riego, en las cuales se hace un predimensionamiento de los dimetros.
5.
Para estos problemas se conoce el caudal demandado, la tubera (se conocen tanto la longitud
como el dimetro, la rugosidad absoluta, los coeficientes de prdidas menores de los accesorios) y
las propiedades del fluido (densidad y viscosidad dinmica), Se tiene que calcular la potencia
necesaria (bomba o diferencia de nivel) para mover dicho caudal a travs de la tubera. Este
problema se plantea cuando se quiere utilizar una tubera 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 tubera.
6.
DISEO DE IA TUBERA
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 caractersticas de la tubera (longitud y
accesorios necesarios con sus correspondientes coeficientes de prdidas menores) y las
propiedades del fluido (densidad y viscosidad dinmica). Se desconoce el dimetro necesario para
permitir el paso del caudal demandado. Por lo general, en cuanto al material de la tubera se tienen
dos o tres alternativas, lo cual implica que se conoce la rugosidad absoluta de la tubera, pero es
necesario hacer tantos diseos como diferentes materiales se encuentren disponibles en el sitio del
proyecto con el fin de escoger la mejor alternativa. ECUACIONES PARA EL DISEO DE
TUBERAS SIMPLES
Mediante el uso de las ecuaciones de Colebrook-White conjuntamente con la ecuacin de
Darcy-Weisbach tenemos:
Esta ltima ecuacin es la base para la solucin de los tres tipos de problemas relacionados con
tuberas simples mencionados anteriormente.
MTODOS NUMRICOS
7.
COMPROBACIN DE DISEO
2.
While:
while expression, statements, endevala una expresin , y repite la ejecucin de un grupo de
instrucciones en un bucle mientras que la expresin es verdadera. Una expresin es verdadera
cuando su resultado es no vaco y slo contiene elementos no nulos (lgicos o numricos
reales). De lo contrario, la expresin es falsa.
Ejemplo:
n = 10;
f = n;
mientras que n> 1
n = n-1;
f = f * n;
final
disp ([ 'n! =' num2str (f)])
norte! = 3628800
For:
for index = values, statements, endejecuta un grupo de instrucciones en un bucle para un nmero
determinado de veces. valuestiene una de las siguientes formas:
initVal: endVal- Incrementar la indexvariable a partir initValde endValpor 1, y repetir la ejecucin
de statementshasta indexes mayor que endVal.
initVal: step: endVal- Incremento indexpor el valor stepde cada iteracin, o
decrementos indexcuando stepes negativo.
valArray- Crear un vector columna, indexdesde las columnas posteriores de serie valArrayen cada
iteracin. Por ejemplo, en la primera iteracin, . El bucle se ejecuta un mximo de veces, donde es
el nmero de columnas de , dados por . La entrada puede ser de cualquier MATLAB tipo de
datos, incluyendo un vector de caracteres, serie de clulas, o
estructura.index = valArray(:,1)nnvalArraynumel(valArray(1,:))valArray
MTODOS NUMRICOS
ejemplo:
for v = 1.0:-0.2:0.0
disp(v)
end
1
0.8000
0.6000
0.4000
0.2000
0
3.
EJERCICIOS
Para hallar el caudal de una tubera.
Se desea calcular el caudal de agua que puede ser movido a travs de una tubera de
PVC, de
12 pulgadas de dimetro nominal y 730 m de longitud, que conecta dos tanques de
abastecimiento de agua potable con una diferencia de nivel de 43.5 m, El dimetro real
-6
de la tubera es de 293 mm y su rugosidad absoluta es de 1.5 x10 m. Todos los
accesorios que forman parte del sistema, incluyendo la entrada y Ia salida, implican un
coeficiente global de prdidas menores km de 11.8. El agua se encuentra a 20C.
Para el agua a 20"C se tienen las siguientes
caractersticas:
Con los datos anteriores se puede seguir el procedimiento del diagrama de flujo; los
resultados de las iteraciones se muestran en la siguiente tabla. Este procedimiento es
fcilmente programable.
Los resultados de este ejemplo se pueden resumir as:
Ht=43.5;
Z2=0.0;
E=0.0001;
g=9.81;
hf=Ht-Z2;
Hf=Ht-Z2;
HFf=Ht-Z2;
L=[730,730];
Ks=[0.0000015,0.0000015];
SKm=[11.8,11.8];
D=[0.293,0.293];
i=1;
while i>=0
for i=1
t=(2.0*sqrt(2.0*g*D(i)*hf))/(sqrt(L(i)));
xx=(2.51*visco*sqrt(L(i)))/(D(i)*sqrt(2.0*g*D(i)*hf));
v1=-t*log10((Ks(i)/(3.7*D(i)))+xx);
hf1=Ht-Z2-SKm(i)*(v1^2/(2*g))
end
for i=2
t=(2.0*sqrt(2.0*g*D(i)*Hf))/(sqrt(L(i)));
xx=(2.51*visco*sqrt(L(i)))/(D(i)*sqrt(2.0*g*D(i)*Hf));
v2=-t*log10((Ks(i)/(3.7*D(i)))+xx);
hf2=Ht-Z2-SKm(i)*(v2^2/(2*g))
end
if abs(hf1-hf)<=E
break
end
hf=hf1;
i=i+1;
if abs(hf2-Hf)<=E
break
end
Hf=hf2;
i=i+1;
end
%Calculo de rea hidrulica
A1=pi()*D(1)^2/4;
A2=pi()*D(2)^2/4;
%Calculo del caudal en m3/s
Q1=v1*A1/2;
Q2=v2*A2/2;
%Perdida de carga por friccin
HF1=hf;
HF2=Hf;
%Clculo de la perdida menores hm1=Ht-HF1
hm=Ht-HF2
%Clculo de la velocidad
v1;
v2;
%Clculo del coeficiente de friccin
f1=(2*g*D(1)*HF1)/(2*L(1)*v1^2);
f2=(2*g*D(2)*HF2)/(2*L(2)*v2^2);
%Clculo del coeficiente de friccin total
f=f1+f2
%Clculo de la velocidad total
MTODOS NUMRICOS
v=v1/2+v2/2
%clculo del caudaltotal%
Qt=Q1+Q2
10
end
fprintf('el diametro es %39.3f\n',d(i));
1.
Ejercicio de caudal.
Trial>> MATLAB_CAUDAL
hf1 =
24.571697169005077
hf2 =
24.571697169005077
hf1 =
33.313549987369839
hf2 =
33.313549987369839
hf1 =
29.324524836729310
hf2 =
29.324524836729310
MTODOS NUMRICOS
11
hf1 =
31.156641257243557
hf2 =
31.156641257243557
hf1 =
30.317460832310250
hf2 =
30.317460832310250
hf1 =
30.702337012828075
hf2 =
30.702337012828075
hf1 =
30.525923044147724
hf2 =
30.525923044147724
hf1 =
30.606806960722260
hf2 =
30.606806960722260
hf1 =
30.569727124063846
MTODOS NUMRICOS
12
hf2 =
30.569727124063846
hf1 =
30.586726698497287
hf2 =
30.586726698497287
hf1 =
30.578933297575411
hf2 =
30.578933297575411
hf1 =
30.582506199995617
hf2 =
30.582506199995617
hf1 =
30.580868203669397
hf2 =
30.580868203669397
hf1 =
30.581619144601582
hf2 =
30.581619144601582
MTODOS NUMRICOS
13
hf1 =
30.581274875430676
hf2 =
30.581274875430676
hf1 =
30.581432705852805
hf2 =
30.581432705852805
hf1 =
30.581360348412026
hf2 =
30.581360348412026
hm =
12.918567294147195
f=
0.011211605578142
v=
4.634649438813166
Qt =
0.312494236703457
Trial>>
Ejercicio de dimetro.
Trial>> MATLAB_DIAMETRO
caudal =
MTODOS NUMRICOS
14
0.0040
A=
0.0044
v1 =
0.9088
caudal =
0.0086
A=
0.0079
v1 =
1.0984
caudal =
0.0252
A=
0.0177
v1 =
1.4287
caudal =
0.0539
A=
0.0314
v1 =
MTODOS NUMRICOS
15
1.7170
caudal =
0.0971
A=
0.0491
v1 =
1.9774
caudal =
0.1567
A=
0.0707
v1 =
2.2174
HF =
1.3730
el dimetro es
Trial>>
0.300
MTODOS NUMRICOS
16
CONCLUSIONES
Se obtienen los mismos resultados, usando MATLAB y usando una tabla de Excel.
MTODOS NUMRICOS
17
BIBLIOGRAFA
MTODOS NUMRICOS
18
ANEXOS
MTODOS NUMRICOS
19
MTODOS NUMRICOS
20