Tesis M Bonilla

UNIVERSIDAD NACIONAL AUTNOMA
DE MXICO
FACULTAD DE INGENIERA
PROPUESTA DE INSTALACIN EXPERIMENTAL
EN EL LABORATORIO DE HIDRULICA PARA LA
ENSEANZA TERICO-PRCTICA DEL
FENMENO DE GOLPE DE ARIETE
TESIS PARA OBTENER EL TTULO DE

INGENIERO CIVIL
PRESENTA
MARIO DAVID BONILLA SNCHEZ
DIRECTOR
M.I. NIKTE NORMA OCAMPO GUERRERO

MXICO, D.F., 2014
DEDICATORIA
El esfuerzo plasmado en esta tesis es un esfuerzo compartido. Si bien, he aprendido que en la vida es
imposible avanzar solo para poder cumplir los objetivos, y que siempre ser necesario tener a las personas
adecuadas a mi lado. Es por ello que dedico todo este trabajo principalmente a mi familiaMam, Pap,
esta dedicatoria la llevaba guardada conmigo desde hace mucho tiempo y hoy que la escribo me causa
nostalgia, pero sobre todo una enorme satisfaccin. Creo fielmente que todo el camino recorrido hasta este
momento ha sido de gran aprendizaje para todos. Con el corazn en la mano les doy las gracias por
ensearme a luchar por mis ideales, a ensearme a enfrentar los retos, a ensearme a enfrentar la vida y a
disfrutarla con todos mis sentidos. Gracias por respetar mis decisiones, gracias por darme la libertad y la
confianza para explorar todos los caminos posibles, gracias por educarme para ser una mejor persona cada
da, s que nunca terminar de aprender de ustedes, pero pueden tener la certeza de que han hecho de m la
persona que yo siempre anhel ser y que a pesar de todos los cambios que pueda tener hasta el ltimo de
mis das, mi esencia siempre tratar de conservarla y siempre los amar y procurar hasta que la vida lo
permita.
De igual forma dedico esta tesis a mis hermanos, Adriana, Felipe y Alejandro. Mania, dedico mi esfuerzo
a tu apoyo en el inicio de mi carrera, por siempre brindarme t clida compaa, por procurarme siempre y
por mostrarme caminos y opciones de vida. Felipe, dedico mi esfuerzo a tu apoyo y a tus palabras siempre
bien intencionadas, y aunque nuestras vidas sean paralelas, agradezco siempre tu respeto a mis decisiones
tomadas. Alex, dedico este trabajo a tu esfuerzo por hacerme un lugar en tu vida en el inicio de este
proyecto y por compartir muchos momentos juntos, has sido siempre una persona muy importante en mi
vida.
En fines imposible escribir toda la gama de sentimientos y pensamientos en una hoja de papel, es por
ello que les digo, Gracias por disfrutar junto conmigo este xito profesional logrado!.
AGRADECIMIENTOS
Agradezco a la Facultad de Ingeniera, UNAM, por darme la oportunidad de formarme profesionalmente y
obtener el conocimiento tcnico y tico brindado por todos y cada uno de sus profesores.
Agradezco a la Mtra. Nikte Norma Ocampo Guerrero por darle la direccin adecuada a la elaboracin de
esta tesis y por toda su paciencia.
Agradezco a mis compaeros de clases, por todos los momentos vividos en esta Facultad de Ingeniera,
UNAM.
Agradezco a mis amigos y amigas, que aunque s que son pocos, pero al final s que siempre voy a poder
contar con ellos.
Agradezco a las personas que en algn momento me brindaron apoyo y palabras de aliento para continuar
adelante.
Agradezco a todas y todos los profesionistas que al inicio y a lo largo de mi vida profesional han
contribuido con compartir conmigo toda su experiencia y conocimiento y por sus ganas de ensearme a
ser un profesionista el cual tenga respeto por su trabajo y nunca olvidar que todo el conocimiento y la
experiencia adquirida debe ser siempre utilizada en beneficio de los dems.
Por ltimo quiero hacer un agradecimiento especial a mi persona. Por permitirme vivir toda experiencia
que me haga crecer, que me permita hacer algo por los dems, que me permita seguir aprendiendo y que
me permita siempre ser un mejor ingeniero.
CONTENIDO
INTRODUCCIN
1. ANTECEDENTES
4
1.1. Marco histrico
4
1.2. Prctica Golpe de ariete en el Laboratorio de Hidrulica de la Facultad de Ingeniera,
UNAM
8
2. CONCEPTOS FUNDAMENTALES DE LA TEORA DEL GOLPE DE ARIETE
2.1. Flujo permanente
2.2. Flujo transitorio
2.3. Ecuacin de continuidad para un conducto
2.3.1. Principio de conservacin de la masa
2.3.2. Ecuacin de continuidad
2.4. Ecuacin dinmica para un conducto
2.5. Descripcin del fenmeno en un conducto por gravedad
2.5.1. Celeridad de onda
2.5.2. Leyes para maniobras de cierre
2.6. Teora de la columna elstica
2.6.1. Ecuacin de continuidad
2.6.2. Ecuacin dinmica
2.7. Ecuaciones diferenciales del golpe de ariete
2.8. Interpretacin fsica de las ecuaciones integrales del golpe de ariete
2.9. Mtodos de solucin
2.9.1. Mtodo de las caractersticas
2.9.1.1.
Ecuaciones caractersticas
2.9.1.2.
Solucin de ecuaciones caractersticas y de compatibilidad
2.9.1.3.
Condiciones de frontera
2.9.1.4.
Condiciones de estabilidad y convergencia
2.9.1.5.
Ejemplo
2.9.2. Mtodo de Allievi
2.9.2.1.
Ecuaciones generales de Allievi
2.9.2.2.
Desarrollo en cadenas de Allievi
10
10
10
10
10
11
13
15
18
21
22
23
25
25
27
30
30
30
34
37
40
41
50
50
51
3. CASO DE ESTUDIO
3.1. Dispositivo experimental tanque-tubera-vlvula esfrica
3.2. Mediciones y registro de lecturas
3.3. Determinacin de condiciones iniciales
3.4. Clculo de celeridad de onda de presin en tubera
3.5. Periodo T de la onda de presin
3.6. Prdidas por friccin en la tubera
3.7. Prdidas locales
3.8. Ley para maniobra de cierre no uniforme
53
53
55
56
57
58
58
70
76
PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

ENSEANZA TERICO PRCTICA DEL FENMENO DE GOLPE DE ARIETE
3.9. Solucin aplicando el mtodo de Allievi

3.9.1. Resultados
3.10. Solucin aplicando el mtodo de las caractersticas
3.10.1. Resultados
85
87
112
122
4. PROPUESTAS DE DISPOSITIVOS EXPERIMENTALES

131
4.1. Mejoras al dispositivo experimental tanque-tubera-vlvula esfrica
131
4.1.1. Implementacin de los instrumentos de medicin y control adecuados
131
4.2. Propuesta de nuevo dispositivo
134
4.2.1. Instalacin hidrulica red-tubera en espiral-vlvula
134
4.2.2. Instrumentos de medicin y control
134
4.2.3. Implementacin de los instrumentos de medicin y control en la instalacin hidrulica
propuesta
145
5. PROPUESTA DE PRCTICAS DE LABORATORIO
149
5.1. Prctica Golpe de ariete aplicando el mtodo de las caractersticas en dispositivo
experimental existente
149
experimental propuesto.
156
6. CONCLUSIONES
183
BIBLIOGRAFA
186

INTRODUCCIN
INTRODUCCIN
En nuestro pas existen necesidades primordiales como el abastecimiento de agua potable a las
pequeas poblaciones o a las grandes ciudades. Actualmente la infraestructura existente para estos
fines carece de planeacin en algunos sectores, la cual es muy importante para evitar futuros
problemas o solucionar situaciones en el presente. Hoy en da existen proyectos en marcha ejercidos
por instituciones privadas y pblicas para rehabilitar y reordenar las redes principales de
abastecimiento, pero para ello, se deben considerar muchas variables, ya sean de tipo tcnico o
incluso social.
Por estas razones, es muy importante tener en cuenta la parte terica de las posibles soluciones, ya
que se pueden obtener simulaciones de los fenmenos fsicos que ocurren, lo que lleva a tomar
mejores decisiones.
El objetivo de esta tesis es, principalmente aplicar un mtodo numrico de solucin para el
fenmeno transitorio conocido como golpe de ariete en el dispositivo experimental instalado en el
Laboratorio de Hidrulica de la Facultad de Ingeniera de la UNAM y realizar una propuesta de
prctica de laboratorio para que los alumnos de la asignatura Hidrulica de Mquinas, conozcan y
empleen el mtodo de solucin de uso generalizado para tal efecto. Tambin se analiza la
conveniencia del dispositivo existente y se propone una nueva instalacin experimental.
Sobre el fenmeno del golpe de ariete se sabe que una de las aportaciones ms importantes sobre la
teora, fue dada a conocer por el italiano Lorenzo Allievi. Sin embargo, con la evolucin de las
computadoras en los aos 50, el investigador Gray, introdujo un procedimiento de anlisis del
fenmeno conocido como Mtodo de las Caractersticas.
Este mtodo es de uso generalizado para el anlisis y es el que se propone para la solucin de la
prctica de laboratorio, pues permite introducir algunas variables que son necesarias para simular el
fenmeno de golpe de ariete de una forma ms apegada a las circunstancias reales, ya que se pueden
incluir condiciones como el efecto de las prdidas por friccin en la tubera y la ley de cierre del
rgano de control. Tambin es posible considerar los cambios geomtricos que pudiera tener una
tubera, como sucede en la operacin real de las redes de distribucin.
En el captulo uno, se presenta un marco histrico sobre el estudio del golpe de ariete, debido que es
importante conocer el tiempo que los investigadores del tema tardaron en observar el fenmeno y
desarrollar algn mtodo para representarlo numricamente o poder dar solucin a situaciones que
se generan en las tuberas producto de este fenmeno. En este mismo captulo se presenta la
situacin actual de la prctica Golpe de Ariete que se desarrolla en el Laboratorio de Hidrulica. En
INTRODUCCIN
esta prctica se incluyen los antecedentes necesarios que el alumno debe poseer para comprender el
fenmeno fsico. As mismo se muestra el desarrollo que debe tener la prctica para dar solucin al
fenmeno transitorio con las ecuaciones de Allievi.
En el captulo dos se muestran los conceptos fundamentales que se utilizarn a lo largo de este
trabajo. Estos conceptos se refieren a la parte terica del fenmeno, as como definiciones sobre la
clasificacin de flujos. Tambin se muestra la deduccin de las ecuaciones de continuidad y la
ecuacin dinmica para un conducto mediante el principio de conservacin de la masa. Se describe
el fenmeno del golpe de ariete en un conducto por gravedad, as como las definiciones de celeridad
de onda y las leyes para maniobras de cierre. Se presenta la deduccin de la ecuacin de
continuidad y la ecuacin dinmica mediante la teora de la columna elstica. A su vez se presentan
las ecuaciones diferenciales del golpe de ariete, as como la interpretacin fsica de las ecuaciones
integrales del mismo fenmeno.
En este mismo captulo se presenta la teora del Mtodo de las Caractersticas y el Mtodo de
Allievi. Para el Mtodo de las Caractersticas se describen las ecuaciones caractersticas, y su
solucin; las ecuaciones de las condiciones de frontera aguas abajo y aguas arriba de la tubera.
Parte importante de este mtodo, son las condiciones de estabilidad y convergencia, mismas que
determinan la obtencin de resultados favorables. Al final de este captulo se presenta un ejemplo
de aplicacin, mismo que sirvi para realizar la calibracin de la hoja de clculo que posteriormente
se utilizara para la resolucin del caso experimental en el laboratorio.
Por otra parte se muestra la deduccin de las ecuaciones generales de Allievi, as como el desarrollo
en cadenas que se debe llevar a cabo para obtener una solucin.
En el captulo tres, se describe el caso de estudio, para lo cual se describe detalladamente el
dispositivo experimental utilizado. Posteriormente se muestra el procedimiento experimental para
determinar las condiciones iniciales de flujo, as como para determinar los valores de algunas
constantes necesarias en el clculo. Se presenta tambin la forma en que fueron calculadas las
prdidas totales de energa en la tubera. En seguida, se detalla la solucin del caso de estudio
mediante los dos mtodos desarrollados tericamente en el captulo dos, haciendo una comparacin
entre los resultados obtenidos con ambos mtodos.
En el captulo cuatro se hacen algunas propuestas para mejorar el dispositivo experimental actual,
as como la implementacin de instrumentos de medicin y control adecuados. Tambin se hace la
propuesta de un nuevo dispositivo y de igual forma la implementacin de instrumentos de medicin
y control.
INTRODUCCIN
Finalmente, en el captulo cinco se presenta una propuesta de prontuario de prctica de laboratorio,

utilizando como mtodo de solucin el Mtodo de las Caractersticas y empleando el dispositivo
experimental actual. Tambin se presenta una propuesta de prctica de laboratorio usando el nuevo
dispositivo experimental y como mtodo de solucin el mtodo de las caractersticas a travs de la
hoja de clculo desarrollada en el captulo dos.
As es como se desarrolla el contenido de este trabajo, el cual puede formar parte del material
didctico utilizado en el Laboratorio de Hidrulica de la Facultad de Ingeniera de la UNAM.

CAPTULO 1. ANTECEDENTES
1
1.1
ANTECEDENTES
Marco histrico
De entre los investigadores y cientficos que contribuyeron con sus estudios y experimentos
relacionados con el golpe de ariete, cabe mencionar a Michaud, quien probablemente fue el primero
en hacer investigaciones de este fenmeno en tuberas, mismas que dio a conocer en una
publicacin que apareci en Lausana, Suiza, en 1878 y cuyo ttulo fue Coups de Blier dans les
Conduites. Etude des moyens employs pour attnuer les effets (Golpe de ariete en tuberas. Estudio
de las maneras de mitigar los efectos), donde establece un anlisis matemtico del golpe de ariete
que tiene lugar en una columna de agua elstica contenida en una tubera rgida.
Posteriormente, el ruso Nicolai Egorovich Joukowsky (1847-1921) llev a cabo importantes
estudios vinculados con este fenmeno. Perfeccion sus estudios en Pars y ms tarde retorn a su
patria para pasar el resto de su vida profesional en Mosc, donde fue profesor de mecnica en el
Instituto Politcnico (1872) y posteriormente en la Universidad (1886).
Con antelacin a los estudios y resultados a que lleg el profesor Joukowsky, Helmholtz, dentro de
sus investigaciones en el campo de la Acstica haba demostrado en 1848, que la celeridad o
velocidad de propagacin del sonido en el seno de un fluido compresible contenido en una tubera,
era una funcin de la elasticidad del material de las paredes de esta ltima, y en ese mismo ao D. J.
Korteweg determin la ecuacin correspondiente.
En tales circunstancias, fue Joukowsky quien en 1897 demostr que cuando en una tubera se
presenta un incremento de presin a causa de un cierre brusco de una vlvula situada en el extremo
aguas abajo, la celeridad de la onda de presin generada por tal maniobra, al igual que el sonido,
tambin es una funcin de la elasticidad de las paredes.
En ese tiempo Joukowsky investig y experiment los efectos producidos por cierres de vlvulas,
cmaras de oscilacin, vlvulas de seguridad, etc., y en cada caso present un anlisis matemtico
riguroso, incluyendo todos los elementos necesarios, de acuerdo con las diferentes condiciones.
Los estudios y experimentos realizados por este cientfico con relacin al golpe de ariete, dieron
como resultado las primeras ecuaciones bsicas en su forma ms correcta. Tales resultados fueron
publicados en 1898 en las Memorias de la Academia Imperial de San Petersburgo y,
posteriormente, la traduccin del ruso al ingls fue hecha por la profesora O. Simin y fue editada

por la American Water Works Association (AWWA) en 1904. Las ecuaciones de Joukowsky son
consideradas como la primera gran aportacin a la teora del golpe de ariete.
El sabio italiano Lorenzo Allievi (1856-1942) es quiz la figura ms importante por sus estudios y
aportaciones a la teora del golpe de ariete. Profundiz en las investigaciones efectuadas por
Joukowsky y en 1902, partiendo de los resultados obtenidos por el cientfico ruso, los cuales
complement con ideas propias, elabor su primer trabajo importante sobre el golpe de ariete
(1903) al cual titul como Teora generale del moto perturbato del lacqua nei tubi in pressione
(Teora general del movimiento perturbado del agua en las tuberas bajo presin), que fue editado
por la sociedad italiana de ingenieros y arquitectos, y que posteriormente fue traducido al alemn
por Dubs y Bataillard y en 1909 fue publicado por la editorial Springer de Berln. El mismo Allievi
tradujo al francs su teora, que fue difundida por la Revue de Mcanique (Diario de mecnica).
Su gran inters por mejorar cada vez ms sus teoras acerca del fenmeno en cuestin, tuvo un gran
resultado con la aparicin en 1913 de su obra Teora del colpo dariete (Teora del golpe de ariete),
editada en las actas del Colegio de Ingenieros y Arquitectos de Italia. Asimismo, es necesario
destacar que su tratado de 1913 es considerado como la base para el desarrollo de los
descubrimientos realizados posteriormente en este campo de relevante importancia en la hidrulica
moderna, y que adems existen otras publicaciones del mismo Allievi, todas ellas de gran utilidad.
Lorenzo Allievi falleci en 1942.
En 1909, E. Braun aport una de las mejores contribuciones al desarrollo de la teora del golpe de
ariete, al introducir el valor de la caracterstica de la tubera () en el trabajo que present en
Stuttgart, bajo el ttulo Druckschwankungen in Rohrleitungen (Las fluctuaciones de presin en las
tuberas). En Berln, en 1910, en la revista alemana Die Turbine (La Turbina), public trabajos
similares a los expuestos por Allievi en su segundo tratado en 1912. Es por ello que en otro artculo
publicado el 16 de agosto de 1934 titulado Bermerkungen zur Theorie des Duckschwankungen in
Rohrleitungen, (Observaciones sobre la teora de las fluctuaciones de presin en las tuberas) Braun
manifiesta una prioridad sobre Allievi en la elaboracin de la teora del golpe de ariete. Sin
embargo, el nombre de Allievi destaca preponderantemente en esta teora, conocida como teora de
Allievi, porque la presentacin de los anlisis y desarrollos, en especial en su trabajo de 1913, es
una obra ejemplar en elegancia y claridad matemtica.
Por otra parte, A.H. Gibson, profesor de Ingeniera de la Universidad de Manchester, alumno y
sucesor de Osborne Reynolds, considerado como uno de los pioneros en Gran Bretaa de la
aplicacin de la mecnica de fluidos a la Ingeniera, escribi numerosas publicaciones acerca de
varias ramas de la Hidrulica y en 1925, en su libro The Mechanical Properties of Fluids (Las
propiedades mecnicas de los fluidos), proporciona una introduccin excelente a la teora del golpe
de ariete.
Posteriormente surgieron los mtodos grficos, cuyo precursor fue el suizo Othmar Schnyder, que
en 1929 desarroll un mtodo grfico de solucin del golpe de ariete en bombas, y en 1932 public
un trabajo donde expone su mtodo para resolver cualquier sistema de tuberas, especialmente para
las que alimentan turbinas.
Es necesario citar tambin al clebre cientfico francs Louis Bergeron (1876-1948), quien fue
profesor en la Escuela Central de Pars y que en 1933 trabajando independientemente de Schnyder,
desarroll el mismo mtodo grfico ideado por este ltimo y que en publicaciones en francs de los
aos 1935 y 1936 llevan su nombre. En 1937, Bergeron demostr que el mtodo grfico de anlisis
que haba desarrollado para resolver problemas de golpe de ariete poda ser aplicado a otros campos
de la teora de la elasticidad y tambin en problemas elctricos. Despus de su fallecimiento, en
1950 la editorial Dunod de Pars public su libro Du coup de blier en hydraulique au coup de
foudre en lectricit (El martillo hidrulico para aplastar a la electricidad).
Otro colaborador en el desarrollo de la teora en cuestin fue el profesor R. W. Angus de la
Universidad de Toronto, en Canad, que mejor notablemente los mtodos grfico y analtico de
solucin, mediante valiosas aportaciones. Prueba de su inters por este fenmeno son varios de sus
escritos, dentro de los cuales destacan Water Hammer in Pipes, including those supplied by
Centrifugal Pumps; Graphical Treatment y Graphical Analysis of Water Hammer Pressures in
Compound and Branched Pipes (Golpe de ariete en tuberas, incluidos los suministrados por
bombas centrfugas, tratamiento grfico y anlisis grfico de las presiones del golpe de ariete en
tuberas compuestas y ramificadas) , publicados en 1937 y 1939 respectivamente.
Cabe citar tambin al francs H. Favre, quien llev a cabo investigaciones en tuberas con
disminucin gradual de dimetro y en 1938 public sus investigaciones bajo el nombre de Theorie
des coups de blier dans les conduites caractristiques linairement variables le long de laxe
(Teora del golpe de ariete en tuberas con caractersticas linealmente variables a lo largo de su eje).
En Cours de mcanique (Cursos de mecnica), explica y desarrolla con detalle la integracin de las
ecuaciones del golpe de ariete elaborada por Riemann.
Entre 1940 y 1960 se publicaron un gran nmero de trabajos acerca del fenmeno en estudio, dentro
de los que destacan los libros publicados por Parmakian, Rich y Jaeger.
En 1953, con la aparicin de las computadoras, Gray introdujo el mtodo de las caractersticas para
poder llevar a cabo el anlisis computarizado del golpe de ariete. En 1967 se publica el libro de
Streeter, V. L. y Wylie, E. B. bajo el ttulo Hydraulic Transients (Transitorios hidrulicos), que

presenta el anlisis de los fenmenos transitorios en hidrulica mediante la aplicacin del mtodo
anteriormente citado y el de la impedancia.
En 1979 M. Hanif Chaudrhy publica su libro Applied Hidraulic Transients (Transitorios hidrulicos
aplicados), que contiene una exposicin completa de los fenmenos transitorios en hidrulica.
En 1985 el Instituto de Ingeniera de la Universidad Nacional Autnoma de Mxico, publica la obra
Diseo y Operacin Hidrulicos de Conducciones de Agua a Presin, que contiene una extensa y
detallada exposicin tanto de los fenmenos transitorios en los conductos a presin, como de los
dispositivos que existen para su control.
Sin dejar de mencionar el trabajo desarrollado por el ingeniero Francisco de Ass Tavera Escobar,
quien en su tesis profesional titulada Teora General del Golpe de Ariete y sus aplicaciones, hace
una brillante descripcin del fenmeno y presenta una exposicin detallada de los mtodos grficos
propuestos por Schnyder y Bergeron.
Finalmente, existe un software llamado Trans desarrollado por el Instituto de Ingeniera el cual
analiza el transitorio y permite establecer la adecuada ubicacin de dispositivos de alivio de acuerdo
a la magnitud del fenmeno transitorio.

1.2
Prctica Golpe de Ariete en el laboratorio de hidrulica de la Facultad de Ingeniera,

UNAM
Actualmente, la carrera de Ingeniera Civil, a travs del Departamento de Ingeniera Hidrulica,

imparte en el sexto semestre del plan de estudios, la asignatura Hidrulica de Mquinas y
Transitorios. Los objetivos de este curso, se basan principalmente en que los alumnos puedan
analizar el funcionamiento y condiciones de operacin de las turbomquinas para su seleccin en
estaciones de bombeo y en centrales hidroelctricas y dimensionar dichas instalaciones; as como
calcular las variaciones de presin debidas al golpe de ariete y comprender el funcionamiento de los
dispositivos que lo reducen.
Los temas de esta asignatura se resumen en los siguientes:
1. Conceptos generales
2. Teora general de las turbomquinas
3. Bombas
4. Turbinas
5. Flujo transitorio en conductos a presin
Especficamente el tema cinco se enfoca a describir y mencionar la importancia de los fenmenos
transitorios, en particular el golpe de ariete, sus ecuaciones bsicas y su solucin, por medio de las
ecuaciones de Allievi y el mtodo de las caractersticas. Tambin se estudian los dispositivos de
alivio que reducen el efecto del golpe de ariete.
Para una mejor comprensin de los temas y sobre todo, para que los alumnos observen la aplicacin
prctica de los conceptos tericos, el Departamento de Ingeniera Hidrulica a travs del
Laboratorio de Hidrulica imparte cuatro prcticas relacionadas con dicha asignatura.
En particular, la prctica nmero 3, titulada Golpe de Ariete, est diseada para observar la
variacin de la carga de presin a lo largo de una tubera, provocada por el cierre y apertura de una
vlvula colocada en su extremo final. Asimismo, se analiza la sobrepresin generada en una tubera
debida al cierre de la vlvula, aplicando el mtodo de las cadenas de Allievi.
Para llevar a cabo la prctica y comprender lo que sucede experimentalmente en el desarrollo de
sta, los alumnos deben poseer antecedentes tericos, tales como el concepto de flujo transitorio en
conductos a presin, as como la teora del golpe de ariete, la celeridad de una onda de presin,
conocer la clasificacin de maniobras para una vlvula y conocer las ecuaciones de Allievi. Aunque
estos conceptos fundamentales son adquiridos por los alumnos en el curso terico, se les
recomienda que, antes de comenzar la realizacin de la prctica, acudan a la bibliografa sugerida
en el temario de la asignatura o bien la que se sugiere en el prontuario de la prctica.
Cada una de las prcticas del laboratorio tiene un prontuario de prctica, donde se especifica el
objetivo de la prctica, los antecedentes tericos que debe tener el alumno, el desarrollo de la
prctica, un procedimiento para presentar la memoria de clculo, bibliografa en la cual se pueden
encontrar los conceptos necesarios, as como un cuestionario de cinco preguntas relacionadas, tanto
con la comprensin de los conceptos tericos, como con el desarrollo de la prctica.
La prctica Golpe de Ariete, utiliza como mtodo de solucin las ecuaciones de Allievi. Sin
embargo, como se ha mencionado anteriormente, en el presente trabajo se aplicar el mtodo de las
caractersticas como mtodo de solucin y se realizar la propuesta de un prontuario para realizar la
prctica Golpe de Ariete con dicho mtodo.

CAPTULO 2. CONCEPTOS FUNDAMENTALES
CONCEPTOS FUNDAMENTALES DE LA TEORA DEL GOLPE DE ARIETE
En este captulo se definen algunos conceptos fundamentales que posteriormente sern utilizados
para el anlisis del golpe de ariete en el caso experimental de estudio, tales como flujo permanente,
flujo transitorio, se establecen las ecuaciones de continuidad y dinmica para un conducto, se
describe el fenmeno de golpe de ariete en un conducto, se explica la teora de la columna elstica,
se muestran las ecuaciones diferenciales del golpe de ariete, as como la interpretacin fsica de las
ecuaciones integrales del golpe de ariete.
2.1
Flujo permanente
El flujo permanente es aquel en el que las caractersticas hidrulicas en una determinada seccin del
conducto permanecen constantes para cualquier instante y no varan con el tiempo, es decir
(
(2.1)
donde
p es la presin manomtrica en el conducto
V es la velocidad
Q es el gasto
2.2
Flujo transitorio
Se define como el flujo no permanente intermedio que tiene lugar entre el permanente inicial y el
final en un conducto. Un ejemplo particularmente importante de este tipo de flujo es el fenmeno
del golpe de ariete.
2.3
Ecuacin de continuidad para un conducto
2.3.1
Principio de conservacin de la masa
Este principio establece que de la totalidad de la masa que en la unidad de tiempo entra a un
volumen finito dentro del flujo, una parte queda almacenada dentro del mismo y el resto sale del
volumen. Lo anterior, tambin puede expresarse como

10
Cantidad neta de masa que

atraviesa la superficie de
frontera del volumen en la
unidad de tiempo.
2.3.2
Rapidez de variacin de la
masa contenida en el
volumen.
=0
Ecuacin de continuidad
En la figura 2.1 se muestra un volumen elemental limitado por la superficie interna del conducto 3,
y por las secciones transversales 1 y 2 normales a la coordenada curvilnea x, que a su vez coincide
con el eje del mismo.
La cantidad neta de masa que atraviesa la superficie A de frontera del volumen en estudio es
[
donde
A es la superficie de frontera del volumen,

V es la velocidad,
es la densidad del agua.
y la rapidez de variacin de la masa dentro del mismo
(

11
b)
CORTE A-A
3
D
1
2
dx
A
+X
a)
Figura 2.1. Volumen elemental limitado.
En consecuencia, de acuerdo con el principio de conservacin de la masa, se puede escribir
(2.2)
Sin embargo, si se acepta que la variacin longitudinal del volumen elemental es despreciable
(
), la ecuacin (2.2) se reduce a
(2.2 a)
Al desarrollar las derivadas parciales de la ecuacin anterior y ordenar trminos, se obtiene
(2.2 b)
As, tambin las derivadas totales de A y se pueden escribir como
Sustituyendo las expresiones anteriores en la ecuacin (2.2b) y dividiendo entre
resulta

12
(2.3)
que es la ecuacin de continuidad para un conducto con flujo no uniforme, no permanente y
compresible.
2.4
Ecuacin dinmica para un conducto
Para obtener la ecuacin dinmica de un conducto slo es necesario establecer el equilibrio de las
fuerzas que actan en la direccin tangencial a la coordenada curvilnea x, sobre la superficie de
frontera que limita el volumen elemental que se muestra en la figura 2.2.
Gradiente hidrulico
1
3
2
dx
Zp
P.H.C.
+x
Figura 2.2.
Estas fuerzas aplicadas en la direccin + x, son las siguientes:

a) La fuerza de superficie resultante de las que actan en las secciones 1 y 2 debida a un
gradiente de presin
(2.4)
b) La fuerza de superficie resultante debida a la resistencia al movimiento, aplicada al

volumen elemental y calculada a partir de la frmula de Darcy-Weisbach
13
(2.5)
donde
f es el factor de friccin
c) La componente tangencial del peso propio del volumen, aplicada en el centro de gravedad
del mismo
(2.6)
donde
Zp es la carga de posicin del conducto
De acuerdo con la segunda ley de Newton, la suma de estas fuerzas debe ser igual a la masa del
elemento
multiplicada por la componente tangencial de la aceleracin, lo cual permite
escribir
[
Si se recuerda que
y se divide la ecuacin anterior entre
(2.7)
, despus de ordenar
trminos resulta
Para una determinada seccin del conducto el valor
(2.8)
representa la energa total
por unidad de peso o carga total que tiene el lquido, medida con relacin a un plano horizontal de
referencia, y el trmino
puede ser definido como | | para tomar en cuenta la inversin en la
direccin del flujo cuando ste es transitorio. De acuerdo con esto, la ecuacin (2.8) adquiere la
forma
14
| |
(2.9)
que es la ecuacin dinmica para un conducto con flujo no uniforme, no permanente y compresible.
2.5
Descripcin del fenmeno en un conducto por gravedad
Con el fin de poder describir el fenmeno del golpe de ariete en un conducto por gravedad, por
facilidad, se puede considerar un conducto con seccin transversal y espesor constantes, con eje
horizontal y longitud L tal como el mostrado en la figura 2.3, donde se tiene un flujo permanente
cuando t=0, con una vlvula o compuerta (rgano de control) situada en el extremo aguas abajo que
permite regular la magnitud del gasto y que es alimentado por un depsito que contiene una gran
masa de agua cuya superficie libre tiene un nivel constante.
Depsito
V0
L
H0
P.H.C.
Figura 2.3. Conducto con seccin transversal y espesor constantes
Suponiendo que se lleva a cabo una maniobra de cierre instantneo total en el rgano de control, la
energa cintica del agua en esta seccin del conducto se transformar en energa potencial, que a su
vez, dar origen a un incremento en la carga piezomtrica original, llevndose a cabo un trabajo
elstico de deformacin en las paredes del conducto y en el agua. Esta conversin de energa
constituye lo que se conoce como onda de presin positiva o directa que se propaga hacia aguas
arriba con velocidad o celeridad a.
Esta onda modifica considerablemente las condiciones originales del conducto aguas arriba del
frente de la misma, ocasionando una deformacin en las paredes del mismo, el agua al comprimirse
aumenta su densidad y la velocidad en esta zona se reduce a cero. Aguas abajo del frente de onda, el
conducto, el agua, la carga piezomtrica y la velocidad del escurrimiento, permanecen en su estado
original (figura 2.4a).

15
Lnea de gradiente
hidrulico
Depsito
V = V0
V=0
L
H0
a
Depsito
V=0
(e) Condiciones para t = 2L / a + t
H0
Depsito
(f) Condiciones para t = 3L / a
Lnea de gradiente
hidrulico
Lnea de gradiente
hidrulico
H0
a
Depsito
V=0
H0
V = V0
V=0
(c) Condiciones para t = L / a + t
(g) Condiciones para t = 3L / a + t
Lnea de gradiente
hidrulico
Depsito
H0
V=0
L
(b) Condiciones para t = L / a
V = V0
Lnea de gradiente
hidrulico
Depsito
H0
Lnea de gradiente
hidrulico
V=0
V = V0
rgano de
control
(a) Condiciones para t + t
Depsito
Lnea de gradiente
hidrulico
Lnea de gradiente
hidrulico
V = V0
H0
Depsito
(d) Condiciones para t = 2L / a
V = V0
H0
(h) Condiciones para t = 4L / a
Figura 2.4. Propagacin de la onda de presin en un conducto causada por un cierre instantneo del rgano de
control.
Cuando el frente de onda llega al depsito en un tiempo
, la totalidad del conducto se
encuentra deformado (figura 2.4b), la velocidad de flujo es igual a cero, la carga piezomtrica es
mayor que la original y el agua se encuentra comprimida.
Como el nivel del agua en el depsito permanece constante, en esta zona no pueden presentarse
modificaciones de ninguna especie, permaneciendo la misma carga; entonces el depsito acta
como una pantalla que refleja totalmente la onda de presin y en este instante se produce la
traslacin de una onda negativa hacia el rgano de control (figura 2.4c), ocasionando que aguas
abajo del frente, el agua fluya hacia el depsito con una velocidad V0 idntica en magnitud a la
inicial, dando lugar a que el conducto y el agua adquieran sus propiedades iniciales en esta regin, y
16
cuando la onda llega al rgano en el instante
, todo el conducto y el agua contenida por ste
se encuentran en su estado original, pero con escurrimiento del rgano hacia el depsito (figura
2.4d).
A esta primera serie de eventos que consisten en la formacin de un frente de onda positivo que
viaja del rgano de control hacia el depsito y un frente de onda negativo que viaja de este ltimo
hacia el primero, se le conoce como fase directa del golpe de ariete.
Ahora bien, como consecuencia del flujo de agua hacia el depsito, se origina una cada de carga y
una contraccin de las paredes del conducto, que adopta dimensiones inferiores a las del estado
original, y una reduccin a cero de la velocidad en la zona aguas abajo del frente de onda negativa
que viaja hacia el depsito (figura 2.4e), al que llega en el instante
cuando el conducto se
encuentra contrado en su totalidad, la velocidad del agua es nula y la densidad de la misma es

mayor que la original (figura 2.4f).
A partir del instante
, el depsito funciona nuevamente como una pantalla y refleja una onda
positiva, volviendo el conducto a su estado original aguas arriba del frente de sta, con un
escurrimiento hacia el rgano y con velocidad
, de tal manera que en el instante
el
conducto vuelve a tener las caractersticas iniciales en toda su longitud (figura 2.4h).
A esta segunda serie de eventos se le conoce como fase inversa del golpe de ariete que, al igual que
la fase directa, tiene una duracin
, a la que se conoce como periodo del conducto T.
Finalmente, es necesario sealar que si no existiese el efecto de la friccin en el conducto que

transforma la energa en forma de calor, el ciclo de eventos antes descrito se repetira
indefinidamente dando lugar a una variacin de la carga piezomtrica con respecto al tiempo en el
rgano de control, tal como se indica en la figura 2.5. Sin embargo, dado que prcticamente todos
los conductos son hidrulicamente rugosos, la variacin de la carga corresponde a la mostrada en la
figura 2.6.

17
Nivel del depsito
H
H0
T=2L/a
0
2L
a
4L
a
6L
a
8L
a
12L
a
10L
a
Figura 2.5. Variacin de la carga piezomtrica en el rgano de control sin considerar prdidas por friccin.
hf
Nivel del depsito
H0
2L
a
4L
a
6L
a
8L
a
10L
a
Figura 2.6. Variacin de la carga piezomtrica en el rgano de control considerando prdidas por friccin.
2.5.1
Celeridad de onda
La celeridad de las ondas de presin en un conducto est definida por la ecuacin (2.10), en la cual
puede observarse que su valor depende tanto de las propiedades elsticas del conducto y el fluido,
como de la geometra del primero.
(2.10)

18
donde
a es la celeridad de la onda de presin
D es el dimetro del conducto
es la densidad del fluido

e es el espesor de las paredes del conducto
Ev es el mdulo de elasticidad volumtrico del fluido
Et es el mdulo de elasticidad del material del conducto
Cuando el fluido es agua dulce y en la mencionada ecuacin se aceptan valores prcticos de
y
se obtiene:
(2.11)
La ecuacin anterior permite calcular la magnitud de la celeridad de la onda de presin en un

conducto de pared delgada, cuyo espesor es menor o igual a la dcima parte del dimetro, es decir,
si
. En la tabla 2.1 se indican los valores del mdulo de elasticidad
para algunos
materiales usados en conductos, y en la tabla 2.2 se proporcionan valores del mdulo de elasticidad
volumtrico
y la densidad para algunos lquidos.

Tabla 2.1. Valores del mdulo de elasticidad Et para algunos materiales
Acero
Material
Et (Pa)
2.06 X 1011
Asbesto cemento
2.40 X 1010
P.V.C.
1.10 X 109
Fierro fundido
9.11 X 1010
Cobre
1.19 X 1011
Bronce
1.03 X 1011
Latn
1.03 X 1011
Zinc
3.63 X 1010
Plomo
1.37 X 1010
Estao
1.27 X 1011
Aluminio
7.05 X 1010
Concreto simple
1.23 X 1010
Madera
6.86 X 109
Hule
3.43 X 109
Vidrio
6.86 X 1010

19

Tabla 2.2. Valores comnmente usados del mdulo de elasticidad volumtrico Ev y de la densidad para algunos
lquidos
Lquido
Agua dulce
Agua salada
Petrleo
Gasolina
(Kg/m3)
1000
1025
850
680
Ev (Pa)
2.18 X 109
2.33 X 109
2.06 X 109
1.39 X 109
Temperatura (C)
20
15
15
15
Por otro lado, es necesario sealar que algunos autores sugieren la aplicacin de la siguiente
ecuacin para el clculo de la celeridad de la onda
(2.11 a)
donde C1 es un parmetro que depende de la relacin de Poisson del material con que est hecho
el conducto y de sus condiciones de apoyo.
Sin embargo, puede aceptarse un valor prctico de C1 igual a la unidad para la gran mayora de los
conductos, con lo que la ecuacin (2.11a) se reduce a la (2.11).
Para algunos de los materiales ms comunes en conductos de pared delgada, a partir de esta ltima
ecuacin, se puede obtener la siguiente expresin
(2.12)
donde Ka vale 0.0106 para conductos de acero, 0.091 para conductos de asbesto-cemento y 1.993
para conductos de P.V.C. Obsrvese que si el valor de Ka fuese igual a cero, para un material con
mdulo de elasticidad infinito, el valor mximo de la celeridad sera de 1,482 m/s, que es la
velocidad con la cual se propaga el sonido en el agua a una temperatura de 20C.
Tratndose de conductos de pared gruesa (figura 2.7), si se desprecia el efecto de la relacin de
Poisson , la celeridad de onda se define como
(
(
)
)
(2.13)
donde R es el radio interior del conducto.

20
Figura 2.7. Tubo de pared gruesa.
Para conductos de concreto reforzado existe alguna incertidumbre debido a la heterogeneidad del
material; sin embargo, para valuar la celeridad de onda se recurre a un conducto de acero
equivalente con un espesor virtual
dado por la siguiente ecuacin

(
Siendo
refuerzo,
2.5.2.
el espesor del tubo de concreto,

y
(2.14)
el de una placa continua de acero equivalente al
los mdulos de elasticidad del concreto y acero de refuerzo respectivamente.
Leyes para maniobras de cierre
Para obtener los valores de la carga piezomtrica y la velocidad en el rgano de control cuando ste
se somete a una maniobra de cierre o apertura, es necesario conocer la ley con la cual se efecta
dicha maniobra; para tal fin, a continuacin se presentan las ecuaciones que pueden utilizarse para
los casos ms comunes que corresponden a condiciones iniciales de apertura o cierre total.
Cierre uniforme o lineal
Para una maniobra de cierre uniforme o lineal se tiene la siguiente ecuacin
(
(2.15)
(2.15 a)
donde
para un cierre parcial, e igual a cero si ste es total.
donde
es la relacin entre el tiempo de maniobra tc y el periodo de la onda de presin T.

21
i es la ley lineal de cierre lento al inicio de la maniobra.

f es la ley de cierre lento al final de la maniobra.
Cierre no uniforme
Cuando la ley de cierre no es uniforme, se tendr una variacin de con respecto al tiempo tal
como la mostrada en la figura 2.8; en este caso resulta conveniente hacer una grfica semejante a la
que se muestra de acuerdo con las caractersticas de la maniobra e interpolar de sta el valor
deseado de , o bien, si se dispone de una computadora se puede simular la maniobra mediante
lneas rectas como las que se muestran en la misma figura.
Figura 2.8. Ley para maniobra de cierre no uniforme.
2.6
Teora de la columna elstica
Esta teora se acerca ms al comportamiento real del fenmeno y ha sido comprobada en

laboratorio. Las ecuaciones de continuidad y dinmica en este caso estn sujetas a las siguientes
hiptesis simplificatorias:
1. El conducto permanece lleno de agua todo el tiempo y la presin mnima en cualquier
seccin siempre es mayor que la de vaporizacin del fluido.
2. Las distribuciones de velocidad y presin en cualquier seccin del conducto son uniformes.
3. Las ecuaciones para el clculo de prdidas de carga cuando el flujo es permanente, tambin
son vlidas cuando ste es transitorio.

22
4. La pared del conducto y el fluido se comportan de una manera elstica lineal y tienen
pequeas deformaciones.
5. El incremento de la presin con respecto a la coordenada curvilnea x resulta pequeo
comparado con el incremento de la misma con respecto al tiempo.
6. El incremento de la carga de velocidad y la densidad del fluido resulta pequeo comparado

con el de la carga piezomtrica.
2.6.1
Ecuacin de continuidad
De acuerdo con la hiptesis 5

(2.16)
Por otra parte, de acuerdo con la fig. 2.1 el incremento de la seccin transversal del conducto se
puede expresar como
(2.17)
Segn la ley de Hooke para slidos que tienen un comportamiento elstico lineal se tiene que
(
(2.18)
Sustituyendo la ecuacin (2.12) en la (2.11), resulta
(2.19)
Observando el diagrama de cuerpo libre de las presiones sobre la pared de la tubera mostrado en la
figura 2.9 se deduce que

23
Como
O bien
(2.20)
Esta ltima ecuacin, al ser sustituida en la ecuacin (2.13) da como resultado
e
p
Figura 2.9. Diagrama de cuerpo libre de las presiones sobre la pared de la tubera
O bien, segn la hiptesis 5

(2.21)
Sustituyendo las ecuaciones (2.16) y (2.21) en la (2.3), y tomando en cuenta la hiptesis 6, esta
ltima ecuacin queda como
(
(2.22)
Recordando que Q=VA y ordenando trminos, la ecuacin anterior se puede escribir como
(2.23)
donde el parmetro a se define como la celeridad de la onda de presin en el conducto y su valor es

24
2.6.2
(2.24)
Ecuacin dinmica
La ecuacin dinmica para esta teora puede deducirse fcilmente de la ecuacin (2.9) si se
considera la hiptesis 6, la definicin de carga piezomtrica H y la relacin Q=VA, las cuales
despus de ordenar trminos permiten escribir
| |
(2.25)
Finalmente, es necesario destacar que para el clculo de las prdidas por friccin en el conducto, se
ha utilizado la ecuacin de Darcy-Weisbach; no obstante, es posible utilizar cualquier expresin de
la forma
| |
con valores de
y b que dependen de la ecuacin utilizada; as, para la
frmula de Manning, que en la prctica tambin se aplica a los conductos a presin, K1 = 8.08 gn2,
m = 1 y b = -10/3, siendo n el coeficiente de rugosidad.
2.7
Ecuaciones diferenciales del golpe de ariete
Con base en las ecuaciones de continuidad y dinmica establecidas para la teora de la columna
elstica, despreciando el efecto de la friccin y haciendo Q=VA, las ecuaciones (2.23) y (2.25) se
pueden escribir como
(2.26)
y
(2.27)
que se conocen como las ecuaciones de continuidad y dinmica del golpe de ariete. Si se recuerda
que
entonces las ecuaciones 2.26 y 2.27 pueden escribirse como

25
(2.26 a)
y
(2.27 a)
Para el caso particular de un conducto con eje horizontal y la carga piezomtrica H valuada con
respecto a un plano horizontal de comparacin que contiene a dicho eje, sta resultar igual a la
carga de presin hp, con lo cual la ecuacin (2.26a) se simplifica como
(2.28)
que es la ecuacin diferencial utilizada por Allievi para conductos con eje horizontal y seccin
transversal constante.
Sin embargo, con el fin de obtener el valor de la carga piezomtrica H en cualquier seccin del
conducto, independientemente del perfil de su eje, a continuacin se analiza la solucin e
interpretacin fsica de las ecuaciones (2.26a y 2.27a), donde la carga de presin hp se obtiene con
slo restar la carga de posicin Zp de la piezomtrica correspondiente.
Estas ecuaciones tienen la forma de la ecuacin denominada de DAlambert, cuya solucin
simultnea general fue obtenida por Riemann y para un sistema tal como el mostrado en la figura
2.10 resulta
(
[ (
)]
(2.29)
(2.29 a)
Las expresiones anteriores permiten determinar la carga piezomtrica y la velocidad en cualquier

seccin de un conducto durante el flujo transitorio en funcin de la coordenada curvilnea x con
origen en el depsito y el tiempo t.

26
H0
x1
x2
+x
x=0
P.H.C.
x=L
Figura 2.10.
2.8
Interpretacin fsica de las ecuaciones integrales del golpe de ariete
Con el fin de ayudar a la mejor comprensin del fenmeno en estudio y deducir ecuaciones de tipo
prctico que permitan la cuantificacin de sus efectos, resulta conveniente interpretar el significado
de las funciones
) y
).
Considerando que por alguna razn pudiera justificarse que la funcin

posible encontrar el efecto que la existencia de
) fuese nula, es
) traera consigo. De esta manera las
ecuaciones (2.29) y (2.29a) tomaran la forma
Si se despeja
(2.30)
(2.30 a)
) en la ecuacin (2.30a) y se sustituye su valor en la (2.30) resulta

(
(2.30 b)

27
Al realizar una maniobra de cierre en el rgano de control de la figura 2.10 cuando t = 0, se tendr
que la velocidad para el flujo permanente inicial
flujo transitorio V, es decir,
ser siempre mayor que la velocidad para el
y por lo tanto
Ahora bien, para un observador que partiera del rgano de control cuando
y viajara a lo largo
del conducto en la direccin x con una velocidad a, en un instante t1 se encontrara en la seccin

(
), donde el valor de la funcin

(
) sera (figura 2.10)
y para otro observador que partiera cuando
y viajara en las mismas condiciones, en el
mismo instante t1 se encontrara en la seccin
De acuerdo con lo anterior se deduce que
) resultando entonces que
(
(
)
) representa una onda de carga positiva que se
propaga con direccin al depsito, de tal manera que para un observador que viaja en la misma
direccin con velocidad a, su magnitud permanecer constante.
Una consideracin similar con
, aceptando que solo existiera
), conduce a la
conclusin de que esta ltima funcin representa una onda de carga negativa que se propaga del
depsito hacia el rgano de control, con un valor constante para un observador que viaja en la
direccin +x con velocidad a.
Por otra parte, como la magnitud de la carga piezomtrica H0 permanece constante en el depsito,
de la ecuacin (2.29) se obtiene para x = 0 y el instante t
( )
(2.31)
( )
Adems, si se considera una onda directa F que parte del rgano de control en el instante t, llega al
depsito cuando
, y se refleja dando origen a una onda f con la misma magnitud pero con
signo opuesto que viaja hacia el rgano al que llega en el instante
, puede afirmarse que
en la seccin del conducto correspondiente a este ltimo resulta vlida la siguiente relacin:
( )
(2.32)

28
Es decir, la magnitud de la onda f en el rgano de control para el instante t es igual a la de la onda F

con signo opuesto que parti del mismo con direccin al depsito
segundos antes.

29
2.9. MTODOS DE SOLUCIN

2.9.1
Mtodo de las caractersticas
2.9.1.1. Ecuaciones caractersticas

La base del mtodo de las caractersticas consiste en transformar las ecuaciones diferenciales
parciales de continuidad (ecuacin 2.23), y dinmica (ecuacin 2.25), en ecuaciones diferenciales
ordinarias y posteriormente resolverlas mediante un esquema explcito de diferencias finitas. Las
dos ecuaciones mencionadas se pueden escribir como
| |
(2.33)
(2.34)
Si se lleva a cabo una combinacin lineal de las ecuaciones (2.33) y (2.34), despus de multiplicar
por gA ambos miembros de la ecuacin (2.34) se tendr que
Lo que da como resultado
(
| |
(2.35)
Si las funciones H=H(x, t) y Q=Q(x, t) son soluciones de las ecuaciones (2.33) y (2.34), las
derivadas totales de estas funciones son
(2.36)
(2.37)
Igualando el primer trmino de la ecuacin (2.35) con la (2.37) y el segundo con la (2.36), se
obtiene

30
(2.38)
(2.39)
de estas dos ltimas ecuaciones se deduce que

(2.40)
(2.41)
(2.42)
Si se sustituyen las ecuaciones (2.38), (2.39), (2.41) y (2.42) en la ecuacin (2.35), resulta
| |
(2.43)
para
(2.44)
y
| |
(2.45)
para
(2.46)
Como puede verse, las ecuaciones diferenciales parciales (2.33) y (2.34) se han transformado en dos
ecuaciones diferenciales ordinarias (2.43) y (2.45), que tienen como variable dependiente al tiempo
t, y que resultan vlidas siempre que sean satisfechas las ecuaciones (2.44) y (2.46)
31
respectivamente. Matemticamente, las dos ltimas ecuaciones representan en el plano x-t dos
lneas rectas con pendiente igual a
positiva
y negativa
(figura 2.11), y se conocen como lneas caractersticas
, a lo largo de las cuales se cumplen las ecuaciones (2.43) y
(2.45), denominadas como ecuaciones de compatibilidad.

t
P(X,t 0 + t)
t0 +t
A(X1 ,t0
t0
X1
B(X2 ,t )
x
X2
Figura 2.11. Lneas caractersticas en el plano x-t.
Fsicamente representan la trayectoria recorrida por una perturbacin iniciada en cualquier seccin
de una tubera; as, por ejemplo, si se considera un conducto por gravedad, como el que se muestra
en la figura 2.12, cuyo rgano de control se somete a una determinada maniobra de cierre o apertura
, en dicho instante se originar una onda de presin que partiendo del rgano de
cuando
control, viaja hacia el depsito situado en el extremo aguas arriba, al que llegar en el instante
.

32
aguas arriba
aguas abajo
Figura 2.12. Conducto por gravedad.
Si la trayectoria de la onda mencionada se grafica en el plano x-t, se obtiene la representacin que

aparece en la figura 2.13, en esta figura se puede observar que las condiciones en la regin I
dependen solo de las del flujo permanente inicial, puesto que las correspondientes a las fronteras
situadas en los extremos aguas arriba y aguas abajo no se modifican al permanecer constante el
nivel de la superficie libre del agua en los depsitos 1 y 2 mientras que en la regin II las
condiciones dependen de las impuestas por el rgano de control. De acuerdo con esto puede
afirmarse que las lneas caractersticas CE y DE dividen al plano x-t en dos regiones, para cada una
de las cuales existe una solucin.
t
E(0,t0+T2)
t0+T2
Regin II
Regin I
t0
C(0,t0)
X=0
D(L,t0)
X=L
Figura 2.13. Conducto por gravedad con excitacin en la frontera aguas abajo.

33
2.9.1.2. Solucin de ecuaciones caractersticas y de compatibilidad

Con relacin a la figura 2.11, si las condiciones del sistema hidrulico cuando t=t0 son las del flujo
permanente inicial, y se conocen para cualquier seccin x del conducto, tal que
clculo de las condiciones para el instante
, el
puede efectuarse mediante un esquema de
diferencias finitas de primer orden, mismo que siguiendo la lnea caracterstica positiva AP permite
escribir
(2.47)
y
(2.48)
y al seguir la lnea caracterstica negativa BP resulta:
(2.49)
(2.50)
Sustituyendo las ecuaciones (2.47) y (2.48) en la ecuacin (2.43) y las ecuaciones (2.49) y (2.50) en
la ecuacin (2.45), realizando el clculo del trmino de friccin en los puntos A y B, y despus de
hacer operaciones resulta
(
(2.51)
y
(
(2.52)
La ecuacin (2.51) puede ser escrita como

(2.53)
y la ecuacin (2.52) como
(2.54)
en las cuales

34
(2.55)
(2.56)
y
(2.57)
Se observa que la ecuacin (2.53) es vlida a lo largo de la lnea caracterstica AP y la ecuacin
(2.54) a lo largo de la lnea caracterstica negativa BP. Los valores de las constantes Cp y Cn se
pueden conocer para cada instante de tiempo, y la constante Ca depende de las propiedades del
conducto. La ecuacin (2.53) ser referida como la ecuacin caracterstica positiva y la ecuacin
(2.54) como la ecuacin caracterstica negativa. En las ecuaciones (2.53) y (2.54), se tienen dos
variables, llamadas, HP y QP. Los valores de estas variables pueden ser determinados resolviendo
simultneamente estas ecuaciones.
(
(2.58)
Ahora, el valor de HP puede ser determinado con cualquiera de las ecuaciones (2.53) o (2.54).
Con el fin de explicar el uso de las ecuaciones (2.53) a (2.58), se considera el sistema mostrado en
la figura 2.12, dividiendo la tubera de descarga en un nmero de tramos iguales con longitud x,
tal como se muestra en la figura 2.14.

35

t
t
t = t0 + 2t
t
t = t0 + t
t
t = t0
x
SECCIN: 1
x
2
x
k-1
x
k
x
k+1
x
n-1
X=0
n
X=L
Puntos interiores (t>t0 )

Condiciones de Frontera aguas arriba
Condiciones de Frontera aguas abajo
Condiciones para el flujo permanente inicial (t=t0)
Figura 2.14. Malla empleada para la aplicacin del mtodo de las caractersticas.
Una vez hecho lo anterior, se determinan en cada seccin x los parmetros h0 y Q0 para las
condiciones del flujo permanente inicial que tiene lugar en el sistema cuando t = t0, y a partir de
estos ltimos se lleva a cabo el clculo de hP y QP en los puntos interiores para el instante t=t0+T,
tomando en cuenta las ecuaciones (2.53) a (2.58). Por lo que se refiere a las condiciones de frontera
del sistema, stas quedarn definidas por la ecuacin (2.53) para la ubicada aguas abajo, y por la
ecuacin (2.54) para la que se encuentra aguas arriba, adems de una ecuacin adicional para cada
frontera que permita establecer una relacin entre hP y QP en las mismas. Estas ecuaciones
adicionales dependern del tipo particular de frontera y sern establecidas en la siguiente seccin
para el caso en el que se encuentra nuestro dispositivo experimental.

36
2.9.1.3. Condiciones de frontera
Depsito con nivel constante
Independientemente de la ubicacin del depsito de almacenamiento, ya sea, aguas arriba o aguas

abajo (figura 2.15), si las prdidas locales por entrada, as como la carga de velocidad resultan
despreciables comparadas con las de friccin, entonces
(2.59)
Gradiente de Energa
2
2
Hdep
Gradiente Hidrulico
Flujo
P.H.C.
Figura 2.15. Depsito con nivel constante.
En la cual Hdep es la carga medida a partir de la superficie libre del agua hasta el plano horizontal
de comparacin. La ecuacin (2.54) para el depsito situado aguas arriba, puede ser escrita como
(2.60)
Sin embargo, si la carga de velocidad o las prdidas locales a la entrada no son despreciables,
entonces, estas pudieran ser consideradas en el anlisis como sigue.
Las prdidas a la entrada estn dadas por la ecuacin
(2.61)
Donde K es el coeficiente de prdidas por entrada. A partir de la figura 2.15 se tiene,

37
(2.62)
La solucin simultnea de las ec. 2.62 y 2.54 resulta
(2.63)
donde
(
(2.64)
Ahora, HP puede ser determinada con la ecuacin (2.62).
Vlvula situada en el extremo final aguas abajo
El gasto que pasa por la vlvula (figura 2.16) para condiciones de flujo permanente inicial est
definido por la ecuacin
(
(2.65)
En la cual el subndice 0 indica la condicin en flujo permanente inicial,

Cd = coeficiente de descarga,
Hov = carga aguas abajo de la vlvula
Av = rea de la vlvula abierta
Gradiente Hidrulico
Vlvula
P.H.C.
H
0v
Flujo
Figura 2.16. Vlvula en el extremo final aguas abajo.

38
Una ecuacin similar a la ecuacin (2.65), puede ser escrita para el estado transitorio como
(
(2.66)
Dividiendo miembro a miembro la ecuacin (2.66) entre la ecuacin (2.65), sacando raz cuadrada a
ambos lados y definiendo la relacin
) , obtenemos
)(
(
(2.67)
Sustituyendo de la ecuacin caracterstica positiva, ecuacin (2.53), por HP en la ecuacin (2.67),

resulta
(2.68)
en la cual
) (
). Resolviendo para QP y despreciando la raz negativa, tenemos

(
(2.69)
Por su parte, HP puede ser determinado con la ecuacin (2.53).

Para calcular las condiciones en estado transitorio para un cierre de vlvula, la curva -t (figura
2.17), puede ser planteada con cualquiera de las formas, ya sea tabular o con una expresin
algebraica. Notemos que = 1 corresponde a la vlvula totalmente abierta, en la cual el gasto que
pasa a travs de ella es Q0 bajo una carga de H0.
1.0
Figura 2.17. Cierre de vlvula.

39
2.9.1.4. Condiciones de estabilidad y convergencia

Con relacin al esquema de clculo que se vio en la seccin anterior, se dice que es convergente si
la solucin de las ecuaciones expresadas en trminos de diferencias finitas se aproxima a la solucin
exacta de las ecuaciones diferenciales originales conforme los valores de x y t tienden a cero. Si
el error de redondeo debido a la representacin de nmeros irracionales por un nmero finito de
dgitos significantes, aumenta en la solucin del esquema de clculo, se dice entonces que es
inestable y, en caso contrario se denomina estable; as mismo, se ha podido comprobar que la
convergencia de un esquema implica su estabilidad y viceversa.
Un criterio recomendable para la convergencia y estabilidad de las ecuaciones en diferencias finitas
para el caso de sistemas simples, tales como los que tienen conductos con caractersticas
geomtricas y elsticas constantes, es el siguiente (ver figuras 2.11 y 2.14)
(2.70)
Las soluciones al esquema de diferencias finitas con mayor precisin que se pueden obtener, se dan
si se satisface el criterio de convergencia de la ecuacin (2.70). De este modo, el criterio de
convergencia y/o estabilidad para las ecuaciones de diferencias finitas (ecuaciones 2.51 y 2.52) est
dado por la expresin
(2.71)
La expresin anterior es conocida como Condicin de estabilidad de Courant.

40
2.9.1.5. Ejemplo
El ejemplo que a continuacin se presenta consiste en determinar las cargas y los gastos a lo largo
de los 5 primeros segundos, para la tubera mostrada en la figura 2.18, considerando x = 100 m y
f = 0.016.
L = 500 m; Q = 0.6 m3/s

D = 0.61 m; a = 1000 m/s
x = 100 m
h0 = 400 m
x
0
Figura 2.18.
Solucin:
Condiciones iniciales para flujo establecido
(
)(
h0 = 400 m
h1 = 400 0.56 = 399.44 m
h2 = h1 0.56 = 398.88 m
h3 = h2 0.56 = 398.32 m
Q0 = Q1 = Q2 = Q3 = Q4 = Q5 = 0.6 m3/s
h4 = h3 0.56 = 397.76 m
h5 = h4 0.56 = 397.20 m

41
(
)(
)
)
|
(I)
(II)
Puntos interiores
(III)
(
(IV)
Condiciones de Frontera
Extremo inicial (0), h0 = 400 m = cte.
(
)(
)
(
(V)
Extremo final (5)

(
(VI)
( )
Las ecuaciones para determinar el gasto Q y la carga h, para el instante de tiempo de 0.1 segundos
son para t = 0.1 s e i = 0
Para t = 0.1 s e i = 0

42
)
(
)|
)|
)|
)|
)|
)|
))
))
)|
))

43
)
(
(
(
)|
)|
))
En la tabla 2.3 se muestran los resultados y las figuras 2.19 a 2.24. se presentan las grficas de la
variacin en la carga piezomtrica.

44
Tabla 2.3. Tabla de resultados del ejemplo de aplicacin del mtodo de las caractersticas
i
t(s)
Cn0
Cp1
Cn1
Cp2
Cn2
Cp3
Cn3
Cp4
Cn4
Cp5
h0
Q0
h1
Q1
h2
Q2
h3
Q3
h4
Q4
h5
Q5
400
0.6
399.44
0.6
398.87
0.6
398.31
0.6
397.75
0.6
397.18
0.6
0.1
-0.547
1.745
-0.545
1.744
-0.544
1.742
-0.542
1.740
-0.540
1.739
400.00
0.600
399.44
0.600
398.87
0.600
398.31
0.600
397.75
0.600
410.26
0.563
0.2
-0.547
1.745
-0.545
1.744
-0.544
1.742
-0.542
1.740
-0.615
1.739
400.00
0.600
399.44
0.600
398.87
0.600
398.31
0.600
410.79
0.563
423.34
0.525
0.3
-0.547
1.745
-0.545
1.744
-0.544
1.742
-0.617
1.740
-0.690
1.739
400.00
0.600
399.44
0.600
398.87
0.600
411.32
0.563
423.84
0.525
436.49
0.488
0.4
-0.547
1.745
-0.545
1.744
-0.618
1.742
-0.691
1.741
-0.765
1.739
400.00
0.600
399.44
0.600
411.85
0.563
424.34
0.525
436.96
0.488
449.63
0.450
0.5
-0.547
1.745
-0.619
1.744
-0.692
1.742
-0.766
1.741
-0.840
1.739
400.00
0.600
412.38
0.563
424.84
0.526
437.43
0.488
450.07
0.450
462.84
0.413
0.6
-0.621
1.745
-0.694
1.744
-0.767
1.742
-0.841
1.741
-0.915
1.740
400.00
0.526
425.33
0.526
437.89
0.488
450.52
0.451
463.26
0.413
476.04
0.375
0.7
-0.695
1.672
-0.768
1.744
-0.842
1.743
-0.916
1.741
-0.990
1.740
400.00
0.452
425.48
0.452
450.96
0.451
463.67
0.413
476.43
0.375
489.30
0.338
0.8
-0.769
1.598
-0.843
1.671
-0.917
1.743
-0.991
1.742
-1.066
1.741
400.00
0.378
425.62
0.378
451.24
0.377
476.82
0.376
489.67
0.338
502.54
0.300
0.9
-0.843
1.524
-0.917
1.597
-0.992
1.670
-1.066
1.742
-1.141
1.741
400.00
0.303
425.75
0.303
451.49
0.303
477.21
0.302
502.89
0.301
515.84
0.263
-0.918
1.450
-0.992
1.523
-1.067
1.597
-1.142
1.670
-1.217
1.742
400.00
0.229
425.87
0.229
451.73
0.228
477.56
0.228
503.37
0.226
529.12
0.225
1.1
-0.992
1.376
-1.067
1.450
-1.142
1.523
-1.217
1.596
-1.292
1.669
400.00
0.154
425.97
0.154
451.94
0.154
477.88
0.153
503.79
0.152
516.86
0.188
1.2
-1.067
1.301
-1.142
1.375
-1.217
1.449
-1.292
1.523
-1.294
1.596
400.00
0.080
426.07
0.080
452.13
0.079
478.16
0.079
491.39
0.114
504.50
0.150
1.3
-1.142
1.226
-1.217
1.301
-1.292
1.375
-1.295
1.449
-1.296
1.523
400.00
0.005
426.16
0.005
452.30
0.004
465.64
0.040
478.88
0.076
491.99
0.113
1.4
-1.217
1.152
-1.292
1.226
-1.295
1.301
-1.296
1.375
-1.298
1.449
400.00
-0.070
426.23
-0.070
439.67
-0.034
453.02
0.002
466.26
0.039
479.39
0.075
1.5
-1.292
1.077
-1.295
1.152
-1.296
1.226
-1.298
1.301
-1.299
1.375
400.00
-0.146
413.52
-0.109
426.96
-0.072
440.29
-0.036
453.53
0.001
466.66
0.037
1.6
-1.294
1.001
-1.296
1.077
-1.298
1.152
-1.299
1.226
-1.300
1.301
400.00
-0.148
400.74
-0.148
414.14
-0.111
427.47
-0.074
440.69
-0.037
453.83
0.000
1.7
-1.296
0.999
-1.298
1.001
-1.299
1.077
-1.300
1.152
-1.301
1.226
400.00
-0.150
400.63
-0.149
401.27
-0.149
414.55
-0.112
427.77
-0.075
427.80
0.000
1.8
-1.298
0.997
-1.299
0.999
-1.300
1.002
-1.301
1.077
-1.226
1.152
400.00
-0.151
400.53
-0.151
401.05
-0.151
401.58
-0.150
401.67
-0.075
401.72
0.000
1.9
-1.299
0.996
-1.300
0.997
-1.301
0.999
-1.226
1.002
-1.152
1.077
400.00
-0.152
400.42
-0.152
400.84
-0.152
388.19
-0.114
375.55
-0.075
375.56
0.000
-1.300
0.994
-1.301
0.996
-1.226
0.997
-1.152
0.999
-1.077
1.002
400.00
-0.153
400.31
-0.153
387.56
-0.115
374.82
-0.077
362.08
-0.039
349.39
0.000
2.1
-1.301
0.994
-1.226
0.995
-1.152
0.996
-1.077
0.998
-1.002
0.999
400.00
-0.154
387.15
-0.116
374.31
-0.079
361.46
-0.040
348.66
-0.002
348.61
0.000
2.2
-1.226
0.993
-1.152
0.994
-1.077
0.995
-1.002
0.996
-0.999
0.998
400.00
-0.079
374.01
-0.079
361.06
-0.042
348.15
-0.004
347.99
-0.002
347.94
0.000
2.3
-1.152
1.067
-1.077
0.993
-1.002
0.994
-0.999
0.995
-0.998
0.996
400.00
-0.005
373.92
-0.005
347.86
-0.004
347.59
-0.003
347.43
-0.001
347.38
0.000
2.4
-1.077
1.142
-1.002
1.067
-0.999
0.993
-0.998
0.994
-0.996
0.995
400.00
0.070
373.85
0.070
360.44
0.034
347.14
-0.002
346.97
-0.001
346.92
0.000
PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA ENSEANZA TERICO PRCTICA DEL FENMENO DE
GOLPE DE ARIETE
45
Tabla 2.3. Tabla de resultados del ejemplo de aplicacin del mtodo de las caractersticas (Continuacin)
i
t(s)
2.5
Cn0
-1.002
Cp1
1.217
1
Cn1
-0.999
Cp2
1.142
2
Cn2
-0.998
Cp3
1.067
3
Cn3
-0.996
Cp4
0.993
4
Cn4
-0.995
Cp5
0.994
5
h0
400.00
2.6
-0.999
1.292
-0.998
1.217
-0.996
1.142
-0.995
1.067
-0.994
0.993
400.00
2.7
-0.998
1.294
-0.996
1.292
-0.995
1.217
-0.994
1.142
-0.993
1.067
2.8
-0.996
1.296
-0.995
1.294
-0.994
1.292
-0.993
1.217
-1.067
2.9
-0.995
1.297
-0.994
1.296
-0.993
1.294
-1.067
1.291
-1.142
1
0.145
h1
386.53
0.147
399.26
400.00
0.149
1.142
400.00
1.217
400.00
2
0.109
h2
373.12
0.147
385.90
399.37
0.149
0.151
399.47
0.152
399.58
Q0
3
0.072
h3
359.83
0.110
372.61
398.74
0.149
0.151
398.95
0.152
399.16
Q1
4
0.036
h4
346.63
0.074
359.43
385.49
0.112
0.150
398.43
0.151
411.78
Q2
5
-0.001
h5
346.57
0.000
0.037
346.34
0.000
372.32
0.074
372.27
0.000
0.149
398.33
0.075
398.28
0.000
0.113
424.38
0.075
424.38
0.000
Q3
Q4
Q5
-0.994
1.299
-0.993
1.297
-1.067
1.296
-1.142
1.294
-1.217
1.291
400.00
0.153
399.69
0.153
412.40
0.115
425.10
0.077
437.81
0.039
450.46
0.000
3.1
-0.993
1.300
-1.067
1.299
-1.142
1.297
-1.217
1.296
-1.291
1.294
400.00
0.154
412.81
0.116
425.61
0.078
438.43
0.040
451.18
0.002
451.25
0.000
3.2
-1.067
1.300
-1.142
1.300
-1.217
1.298
-1.291
1.297
-1.294
1.296
400.00
0.079
425.91
0.079
438.83
0.041
451.69
0.004
451.86
0.002
451.90
0.000
3.3
-1.142
1.226
-1.217
1.300
-1.291
1.300
-1.294
1.298
-1.296
1.297
400.00
0.005
426.01
0.005
451.98
0.004
452.26
0.003
452.41
0.001
452.47
0.000
3.4
-1.217
1.152
-1.291
1.226
-1.294
1.300
-1.296
1.300
-1.297
1.298
400.00
-0.070
426.07
-0.070
439.45
-0.034
452.70
0.002
452.87
0.001
452.92
0.000
3.5
-1.291
1.077
-1.294
1.152
-1.296
1.226
-1.297
1.300
-1.298
1.300
400.00
-0.145
413.43
-0.108
426.79
-0.072
440.06
-0.036
453.21
0.001
453.27
0.000
3.6
-1.294
1.002
-1.296
1.077
-1.297
1.152
-1.298
1.226
-1.300
1.300
400.00
-0.147
400.74
-0.147
414.06
-0.110
427.31
-0.073
440.46
-0.037
453.50
0.000
3.7
-1.295
1.000
-1.297
1.002
-1.298
1.077
-1.299
1.152
-1.300
1.226
400.00
-0.149
400.63
-0.149
401.26
-0.148
414.47
-0.111
427.60
-0.074
427.65
0.000
3.8
-1.297
0.998
-1.298
1.000
-1.299
1.002
-1.300
1.077
-1.226
1.152
400.00
-0.150
400.52
-0.150
401.05
-0.150
401.57
-0.149
401.67
-0.074
401.71
0.000
3.9
-1.298
0.997
-1.299
0.998
-1.300
1.000
-1.226
1.003
-1.152
1.077
400.00
-0.151
400.42
-0.151
400.83
-0.151
388.26
-0.113
375.70
-0.075
375.69
0.000
-1.299
0.995
-1.300
0.997
-1.226
0.998
-1.152
1.000
-1.077
1.003
400.00
-0.152
400.31
-0.152
387.64
-0.115
374.97
-0.077
362.29
-0.038
349.69
0.000
4.1
-1.300
0.994
-1.226
0.995
-1.152
0.997
-1.077
0.998
-1.003
1.000
400.00
-0.153
387.23
-0.116
374.46
-0.078
361.68
-0.040
348.97
-0.002
348.90
0.000
4.2
-1.226
0.994
-1.152
0.994
-1.077
0.996
-1.003
0.997
-1.000
0.998
400.00
-0.079
374.17
-0.079
361.28
-0.041
348.47
-0.004
348.29
-0.002
348.25
0.000
4.3
-1.152
1.068
-1.077
0.994
-1.003
0.994
-1.000
0.996
-0.998
0.997
400.00
-0.005
374.06
-0.005
348.18
-0.004
347.89
-0.003
347.75
-0.001
347.68
0.000
4.4
-1.077
1.142
-1.003
1.068
-1.000
0.994
-0.998
0.994
-0.997
0.996
400.00
0.070
374.01
0.070
360.66
0.034
347.46
-0.002
347.28
-0.001
347.24
0.000
4.5
-1.003
1.216
-1.000
1.142
-0.998
1.068
-0.997
0.994
-0.996
0.994
400.00
0.144
386.61
0.108
373.29
0.072
360.05
0.035
346.95
-0.001
346.88
0.000
4.6
-1.000
1.291
-0.998
1.216
-0.997
1.142
-0.996
1.068
-0.994
0.994
400.00
0.146
399.27
0.146
385.98
0.110
372.78
0.073
359.65
0.037
346.66
0.000
4.7
-0.999
1.293
-0.997
1.291
-0.996
1.216
-0.994
1.142
-0.994
1.068
400.00
0.148
399.37
0.148
398.74
0.148
385.58
0.111
372.48
0.074
372.42
0.000
4.8
-0.997
1.295
-0.996
1.293
-0.995
1.291
-0.994
1.216
-1.068
1.142
400.00
0.150
399.48
0.150
398.96
0.149
398.44
0.148
398.34
0.074
398.29
0.000
4.9
-0.996
1.297
-0.995
1.295
-0.994
1.293
-1.068
1.291
-1.142
1.216
400.00
0.151
399.58
0.151
399.17
0.150
411.70
0.113
424.23
0.074
424.25
0.000
-0.995
1.298
-0.994
1.296
-1.068
1.295
-1.142
1.293
-1.216
1.291
400.00
0.152
399.69
0.152
412.33
0.114
424.95
0.076
437.60
0.038
450.16
0.000
PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA ENSEANZA TERICO PRCTICA DEL FENMENO DE
GOLPE DE ARIETE
46

hpz (m)
Carga Piezomtrica hpz
560
540
Seccin 1
520
Seccin 2
Seccin 3
Seccin 4
Seccin 5
500
480
460
440
420
400
380
360
340
320
300
0
0.5
1.5
2.5
3.5
4.5
Tiempo (s)
Figura 2.19.
Carga Piezomtrica hpz - Gasto Q

Seccin 1
hpz (m)
Q (m3 /s)
440
0.7
430
0.6
Carga Seccin 1
Gasto Seccin 1
420
0.5
0.4
410
0.3
400
0.2
390
0.1
380
370
-0.1
360
-0.2
0
0.5
1.5
2.5
Tiempo (s)
Figura 2.20.
3.5
4.5

47

Seccin 2
hpz (m)
Q (m3 /s)
500
0.7
450
0.6
400
0.5
350
0.4
300
Carga Seccin 2
250
Gasto Seccin 2
0.3
0.2
200
0.1
150
100
-0.1
50
0
-0.2
0
0.5
1.5
2.5
Tiempo (s)
Figura 2.21.
3.5
4.5

Seccin 3
hpz (m)
600
Q (m3 /s)
0.7
0.6
Carga Seccin 3
Gasto Seccin 3
500
0.5
400
0.4
0.3
300
0.2
200
0.1
0
100
-0.1
-0.2
0
0.5
1.5
2.5
Tiempo (s)
Figura 2.22.
3.5
4.5

48
Carga Piezomtrica hpz- Gasto Q

Seccin 4
hpz (m)
600
Q (m3 /s)
0.7
0.6
500
0.5
400
0.4
0.3
300
Carga Seccin 4
0.2
Gasto Seccin 4
200
0.1
0
100
-0.1
-0.2
0
0.5
1.5
2.5
Tiempo (s)
Figura 2.23.
3.5
4.5

Seccin 5
hpz (m)
Q (m3 /s)
600
0.7
0.6
500
0.5
400
0.4
300
0.3
Carga Seccin 5
0.2
Gasto Seccin 5
200
0.1
100
-0.1
0
0.5
1.5
2.5
Tiempo (s)
Figura 2.24.
3.5
4.5

49
2.9.1
Mtodo de Allievi
2.9.1.2. Ecuaciones generales de Allievi

El conocimiento de las funciones F y f es difcil tenerlo a la mano para resolver un determinado
problema, sin embargo, Allievi propuso un sistema de ecuaciones muy simple cuya solucin
permite calcular la variacin de la carga piezomtrica y la velocidad en la seccin adyacente
inmediatamente aguas arriba del rgano de control que se muestra en la figura 2.10. Si en las
ecuaciones (2.29) y (2.29a) se hace x=L, resulta
(
(2.72)
y
[ (
(2.73)
)]
expresiones que se conocen con el nombre de Ecuaciones de Allievi.

Por otra parte, si se hace
la ecuacin (2.32) puede escribirse como

( )
[(
) ]
donde i es un nmero adimensional entero o fraccionario, adems, si se define

[(
) ]
( )
, la ecuacin anterior queda

(2.74)
Al sustituir esta ecuacin en la (2.72) para instantes i e i-1 se obtiene

(2.75)
y
(2.76)
sumando miembro a miembro estas ltimas ecuaciones y haciendo operaciones resulta

50
(2.77)
Si se hace un razonamiento semejante con la ecuacin (2.73) se llega a
(
(2.78)
y al igualar las ecuaciones (2.77) y (2.78) resulta finalmente

(
(2.79)
La ecuacin anterior es la frmula clsica de Allievi y permite llevar a cabo un desarrollo en cadena
mediante el cual se puede obtener la carga piezomtrica en la seccin adyacente al rgano de
control para el instante i, si se conoce su valor para el instante i-1, y el incremento de velocidad
entre dichos instantes, mismo que estar determinado por la ley de cierre o apertura en el rgano de
control que se analiza ms adelante.
Es necesario subrayar que en la ecuacin original de Allievi el valor de la carga piezomtrica H que
aparece en la ecuacin (2.79), corresponde al de la carga de presin hp en un conducto de eje
horizontal; no obstante, esta ltima ecuacin es vlida para cualquier perfil del eje y se reduce a la
de Allievi, si el plano horizontal de comparacin se elige de tal manera que contenga al menos un
punto del primero en la seccin en estudio (figura 2.10), lo que da como resultado que en el rgano
de control se tenga entonces H=hp.
2.9.2.2. Desarrollo en cadenas de Allievi
Si se dividen ambos miembros de la ecuacin (2.79) entre H0 se obtiene
(
Al introducir el valor de
en el segundo miembro de la ecuacin anterior resulta

(
haciendo
, esta ltima ecuacin toma la forma

(
(2.80)

51
Ahora bien, si se aplica la ecuacin de continuidad para una seccin transversal ubicada
inmediatamente aguas arriba del rgano de control, y se toma en cuenta lo mencionado acerca del
plano horizontal de comparacin, se puede escribir
(
(
o bien
(
(2.81)
Sustituyendo la ecuacin anterior en la ecuacin (2.80), se tendr:

(
(2.82)
ecuacin que se conoce con el nombre de Ecuacin Adimensional de Allievi.

52
CAPTULO 3. CASO DE ESTUDIO
CASO DE ESTUDIO
Como se mencion en el primer captulo, se analizar el fenmeno del golpe de ariete con el
mtodo de las caractersticas, para lo cual, se cuenta con el antecedente del anlisis del fenmeno
transitorio con las ecuaciones de Allievi. En este captulo se realizar la solucin del caso
experimental con los dos mtodos citados y se har una comparacin entre los resultados obtenidos.
El objetivo principal es observar la variacin de la carga de presin a lo largo de la tubera,
provocada por el cierre de la vlvula. Con el mtodo de las cadenas de Allievi, se calcular y
analizar la sobrepresin generada en la tubera por el cierre de la vlvula y por otra parte, con el
mtodo de las caractersticas se determinar la variacin de la carga de presin a lo largo de la
tubera, utilizando las prdidas por friccin y condiciones de frontera variables.
A continuacin se presenta una descripcin del dispositivo experimental, el cual se ha utilizado para
observar el fenmeno transitorio del golpe de ariete, en la prctica nmero 3 del laboratorio de la
asignatura Hidrulica de mquinas.
3.1
Dispositivo experimental tanque tubera vlvula esfrica
En la figura 3.1 se observa un esquema del dispositivo experimental, el cual se encuentra instalado
en el Laboratorio de Hidrulica de la Facultad de Ingeniera, UNAM.
Este dispositivo est formado por un tanque o depsito con carga constante, una tubera de cobre de
4.3 metros de longitud y 2.5 pulgadas de dimetro, 5 piezmetros distribuidos a lo largo de la
tubera, en los cuales se mide la columna de agua que se genera por la presin de la onda, un
accesorio conocido como Te bridada, y una vlvula de esfera, colocada en el extremo final de la
tubera. Como dispositivo de aforo, se cuenta con un vertedor de pared delgada de seccin
triangular.

53
Piezmetros
Tubera de cobre, L=4.3 m, D=2.5"
Te con bridas, D=2.5"
Vlvula de esfera, D=2.5"
Vertedor triangular B=1.7, =45
Nsa
Nc
Figura 3.1. Dispositivo experimental
Fotografa 3.1. Perspectiva del dispositivo experimental
En la fotografa 3.1, se puede observar una perspectiva del dispositivo experimental, as como los
componentes de este.
A continuacin se presentan las mediciones que se realizaron en el dispositivo al simular el
fenmeno transitorio.

54
3.2
Mediciones y registro de lecturas
El procedimiento que se realiz para determinar las condiciones iniciales, es el que se utiliza
actualmente en el prontuario de la prctica golpe de ariete.
1. Abrir la vlvula de alimentacin del tanque y establecer una carga constante H0 y registrar su
valor, en m.
2. Medir el nivel de cresta Nc, en m, en el vertedor triangular.
3. Abrir totalmente la vlvula de esfera ubicada en el extremo final aguas abajo y purgar los
piezmetros para establecer el flujo permanente. Medir la carga de presin H, en m, en cada
uno de ellos.
4.
Medir el nivel de la superficie libre del agua Nsa, en m, en el canal de aproximacin del
vertedor de seccin triangular.
Se registraron los siguientes datos

Tabla 3.1. Registro de datos (Flujo permanente)
Depsito o tanque con

carga constante.
Vertedor triangular
Piezmetros
H0=
1.00
N c=
0.2479
Nsa=
0.3791
H1=
0.86
H2=
0.81
H3=
0.755
H4=
0.61
Para el caso de flujo no permanente se realiz un cierre brusco en la vlvula de esfera con el fin de
medir la carga mxima de presin que se present. Esta lectura nicamente se tom en el
piezmetro 4 (ver figura 3.1) y el dato obtenido fue Hmx = 1.5 m

55
3.3
Determinacin de condiciones iniciales
La ecuacin que se utiliza para el clculo del gasto en un vertedor de pared delgada de seccin
triangular es:
(3.1)
donde:
C
coeficiente de descarga del vertedor, en m1/2/s
carga en el vertedor,
aceleracin de la gravedad, en m/s2;
ngulo del vrtice del vertedor, = 45;
ancho del canal de aproximacin, B = 1.70 m;
coeficiente experimental que depende de y h segn figura 3.2
coeficiente experimental que depende de B y h segn figura 3.3
(3.2)
|, en m;
Figura 3.2. Coeficiente de gasto de vertedores triangulares en la frmula de la Universidad Catlica de Chile
Figura 3.3. Valores de K en la frmula de la Universidad Catlica de Chile para vertedores triangulares
56
Carga en el vertedor
Coeficiente experimental , con h = 0.1312 m y
Coeficiente experimental K, con
rad., de la figura 3.2
Coeficiente de descarga del vertedor C

( )
[
](
)(
Gasto Q0
(
)(
De esta forma, se determinaron las condiciones iniciales correspondientes al flujo permanente en el

dispositivo experimental, que se muestra en la tabla 3.2.
Tabla 3.2. Condiciones iniciales a flujo permanente
(m3/s)
H0
(m.c.a.)
H1
(m.c.a.)
H2
(m.c.a.)
H3
(m.c.a.)
H4
(m.c.a.)
0.00376
1.00
0.86
0.81
0.755
0.61
Q0
3.4
Clculo de celeridad de onda de presin en tubera
Se sabe que el valor de la celeridad de la onda de presin depende de las propiedades elsticas del
conducto y el fluido, as como de la geometra de la tubera. En el dispositivo experimental, todas
estas propiedades son constantes a lo largo de la tubera, por lo tanto se obtuvo un solo valor de la
celeridad a de la onda de presin, en m/s.
(3.3)

57
donde:
densidad del agua, = 1,000 Kg/m3;
Ev
mdulo de elasticidad del agua,
Et
mdulo de elasticidad de la tubera de cobre,
dimetro de la tubera, D = 0.0635 m;
espesor de la tubera, e = 0.0018 m.
;
;
Por lo tanto
)(
)(
)
)
3.5
Periodo T de la onda de presin
Se sabe que el periodo T de la onda de presin, es el tiempo que tarda la onda de perturbacin en
recorrer la distancia desde la seccin de la vlvula hasta el tanque y reflejarse hasta llegar
nuevamente a la vlvula. De este modo
Segn esta definicin con la distancia
y si la velocidad V es representada por la celeridad a
de la onda de presin, entonces, el periodo T, queda definido con la siguiente ecuacin

(3.4)
Donde, de acuerdo al caso de estudio
L = 4.3 m
a = 1,150.17 m/s
Por lo tanto
3.6
Prdidas por friccin en la tubera
Para obtener una estimacin numrica del fenmeno en estudio lo ms cercano a la realidad, fue
necesario obtener una relacin entre las mediciones realizadas a flujo permanente y los resultados
58
obtenidos en forma terica, por ello es que en este apartado, se realiza el clculo de las prdidas por
efecto de la friccin entre el fluido y la tubera, con las ecuaciones de friccin de Darcy-Weisbach,
Manning-Strickler y Hazen-Williams. Al obtener los resultados de las prdidas por friccin, adems
de calcular tambin las prdidas locales, se realiz una comparacin y se defini la ecuacin con
mayor aproximacin a las mediciones realizadas.
Tres caractersticas geomtricas de la seccin en una conduccin hidrulica muy importantes en el
clculo de las prdidas de friccin, son las siguientes.
rea hidrulica A, es decir, el rea de la seccin transversal ocupada por el lquido dentro del
conducto.
Permetro mojado P, que es el permetro de la seccin transversal del conducto en el que hay
contacto del lquido con la pared (no incluye la superficie libre si sta existe).
Radio hidrulico Rh, o sea la relacin entre el rea hidrulica y el permetro mojado de la seccin
(Rh=A/P)
Los valores de estas caractersticas geomtricas para la tubera en estudio son:
rea hidrulica A
Permetro mojado P
Radio hidrulico Rh

59
Ahora bien, para obtener las prdidas por friccin, se debe conocer el factor de friccin de la
tubera.
El factor de friccin f es funcin de la rugosidad y del nmero de Reynolds Re en el tubo.
Por otra parte, se sabe que existen distintos intentos de evaluar el efecto corrosivo del agua en
conductos, basndose en la reduccin del gasto calculado tericamente, de acuerdo con el PH del
agua y el nmero de aos de servicio de la tubera. Sin embargo, el criterio que parece ms efectivo
es el de Genijew al modificar la rugosidad absoluta del tubo nuevo, de acuerdo con el tipo de agua
que va a escurrir y el nmero de aos de servicio; esto es
(3.5)
donde
0
rugosidad del tubo (nuevo), en mm;
coeficiente que depende del grupo en el que se clasifique el agua que va a escurrir, segn la
tabla 3.3;
nmero de aos de servicio de la tubera;
rugosidad del conducto, despus de t aos de servicio, en mm.
Segn la tabla 3.3, para la tubera lisa de cobre, la rugosidad absoluta (tubera nueva de cobre) es

60

Tabla 3.3. Rugosidad absoluta en tubos comerciales
Material
, en mm
De vidrio, cobre, latn, madera (bien cepillada), acero nuevo

soldado y con una mano interior de pintura; tubos de acero de
precisin sin costura, serpentines industriales, plstico, hule
0.0015
Tubos industriales de latn
0.025
Tubos de madera
0.2 a 1
Tubos lisos
Hierro forjado
0.05
Fierro fundido nuevo
0.25
Fierro fundido, con proteccin interior de asfalto
0.12
Fierro fundido oxidado
1 a 1.5
Fierro fundido, con incrustaciones
1.5 a 3
Fierro fundido, centrifugado
0.05
Fierro fundido nuevo, con bridas o juntas de macho y campana
0.15 a 0.3
Fierro fundido usado, con bridas o juntas de macho y campana
2 a 3.5
Fierro fundido para agua potable, con bastantes incrustaciones y

dimetro de 50 a 125 mm
1a4
Fierro galvanizado
0.15
Acero rolado, nuevo
0.05
Acero laminado, nuevo

Acero laminado con proteccin interior de asfalto
0.04 a 0.1
0.05
Tubos de acero soldado de calidad normal

Nuevo
0.05 a 0.10
Limpiado despus de mucho uso
0.15 a 0.20
Moderadamente oxidado, con pocas incrustaciones

Con muchas incrustaciones
Con remaches transversales, en buen estado
0.4
3
0.1
Con costura longitudinal y una lnea transversal de remaches en

cada junta, o bien laqueado interiormente
0.3 a 0.4
Con lneas transversales de remaches, sencilla o doble; o tubos

remachados con doble hilera longitudinal de remaches e hilera
transversal sencilla, sin incrustaciones
0.6 a0.7
Acero soldado, con una hilera transversal sencilla de pernos en cada

junta, laqueado interior, sin oxidaciones, con circulacin de agua
turbia

61

Tabla 3.3. (Continuacin)
Material
, en mm
Acero soldado, con costura doble de remaches transversales, muy

oxidado. Acero remachado, de cuatro a seis filas longitudinales de
remaches, con mucho tiempo de servicio.
Tubos remachados, con filas longitudinales y transversales

a) Espesor de lmina < 5 mm
0.65
b) Espesor de lmina de 5 a 12 mm
1.95
c) Espesor de lmina > 12 mm, o entre 6 y 12 mm, si las hileras de

pernos tienen cubrejuntas
d) Espesor de lmina > 12 mm con cubrejuntas
Tubos remachados, con cuatro filas transversales y seis
longitudinales con cubrejuntas interiores
3
5.5
4
Asbesto - cemento nuevo
0.025
Asbesto - cemento, con proteccin interior de asfalto
0.0015
Concreto centrifugado, nuevo

Concreto centrifugado, con proteccin bituminosa
0.16
0.0015 a 0.125
Concreto en galeras, colado con cimbra normal de madera
1a2
Concreto en galeras, colado con cimbra rugosa de madera
10
Concreto armado en tubos y galeras, con acabado interior

cuidadosamente terminado a mano
0.01
Concreto de acabado liso
0.025
Conductos de concreto armado, con acabado liso y varios aos de

servicio
Concreto alisado interiormente con cemento
Galeras con acabado interior de cemento
0.2 a 0.3
0.25
1.5 a 1.6
Concreto con acabado normal
1a3
Concreto con acabado rugoso
10
Cemento liso
0.3 a 0.8
Cemento no pulido
1a2
Concreto presforzado Freyssinet
0.04
Concreto presforzado Bona Socoman
0.25
Mampostera de piedra, bien junteada
1.2 a 2.5
Mampostera de piedra rugosa, sin juntear
8 a 15
Mampostera de piedra, mal acabada
1.5 a 3

62
Tambin se sabe, que la tubera bajo anlisis tiene aproximadamente 15 aos en servicio.
En la tabla 3.4, se observan distintos coeficientes a de la frmula de Genijew. Se sabe que el agua
utilizada en el Laboratorio de Hidrulica es agua poco corrosiva y con un pequeo contenido de
materia orgnica, por lo tanto, se usar el valor medio del grupo I, es decir
a = 0.025.
Tabla 3.4. Coeficientes a de la frmula de Genijew
Grupo I
Agua con poco contenido mineral que no origina corrosin. Agua con un pequeo contenido de
materia orgnica y de solucin de hierro:
a vara de 0.005 a 0.055; valor medio, 0.025.
Grupo II
Agua con poco contenido mineral que origina corrosin. Agua que contiene menos de 3 mg/lt de
materia orgnica y hierro en solucin:
Grupo III
Agua que origina fuerte corrosin y con escaso contenido de cloruros y sulfatos (menos de 100 a
150 mg/lt). Agua con un contenido de hierro de ms de 3 mg/lt:
Grupo IV
Agua que origina corrosin, con un gran contenido de sulfatos y cloruros (ms de 500 a 700
mg/lt). Agua impura con una gran cantidad de materia orgnica:
Grupo V
Agua con cantidades importantes de carbonatos, pero de dureza pequea permanente, con residuo
denso de 2000 mg/lt:
a vara de 0.6 a ms que 1.
De esta forma, se puede calcular la rugosidad del conducto, de acuerdo con los aos de servicio y el
tipo de agua utilizada.
Sustituyendo los valores correspondientes en la ecuacin 3.5
(

63
El factor de friccin f es funcin del nmero de Reynolds y en el caso de un conducto cilndrico a

presin, este nmero se define como
(3.6)
Donde V es la velocidad media, D es el dimetro del conducto y la viscosidad cinemtica del
fluido.
De la ecuacin de continuidad
, la velocidad
, por lo cual
La viscosidad cinemtica del fluido , se obtiene de la figura 3.4, para una temperatura promedio de
20C.
Figura 3.4. Viscosidad cinemtica del agua y del aire a la presin atmosfrica del nivel del mar
Una vez determinada la velocidad y la viscosidad cinemtica del fluido, y sustituyendo estos
valores en la ecuacin 3.6, el nmero de Reynolds tom el siguiente valor
(

64
A partir del diagrama de Moody (figura 3.5), y con la rugosidad relativa , se determin el factor
de friccin correspondiente a la tubera del dispositivo experimental
Del diagrama de Moody con el valor de la relacin anterior y el nmero de Reynolds obtenido
anteriormente, el valor del factor de friccin
Figura 3.5. Coeficiente de friccin para cualquier tipo y tamao de tubo; diagrama universal de Moody
Una vez determinado el factor de friccin, se realiz el clculo de las prdidas por friccin. Para
ello se dividi la tubera en cinco tramos. Como se puede observar en la figura 3.6, se encuentran
colocados cinco piezmetros a lo largo de la tubera, y ello permiti establecer los tramos de
estudio.

65

pz 0
pz 1
pz 2
TRAMO
1
TRAMO
2
0,88
1,1
pz 3
TRAMO
3
pz 4
vlvula
TRAMO TRAMO
4
5
Ho
1,53
0,45
0,34
Figura 3.6. Distribucin de piezmetros
Como se mencion anteriormente, se utilizaron tres frmulas empricas para calcular las prdidas
por friccin. A continuacin se muestra el clculo realizado con la frmula de Darcy-Weisbach.
Frmula de Darcy-Weisbach
Para un flujo permanente, en un tubo de dimetro constante, la lnea de cargas piezomtricas es
paralela a la lnea de energa e inclinada en la direccin del movimiento. En 1850, Darcy, Weisbach
y otros, dedujeron experimentalmente una frmula para calcular en un tubo la prdida por friccin
(3.7)
donde
f
factor de friccin, sin dimensiones;
dimetro, en m;
longitud del tubo, en m;
hf
prdida por friccin, en m;
velocidad media, en m/s.
En la siguiente tabla, se muestran las prdidas por friccin calculadas con la frmula de DarcyWeisbach

66

Tabla 3.5. Prdidas por friccin con la frmula de Darcy-Weisbach
Longitud L
(m)
Factor de friccin f
f L/D
hf
(m)
Tramo 1
0.88
0.035
0.485
0.035
Tramo 2
1.1
0.035
0.606
0.044
Tramo 3
1.53
0.035
0.843
0.061
Tramo 4
0.45
0.035
0.248
0.018
Tramo 5
0.34
0.035
0.187
0.013
Total
4.3
0.170
Frmula de Manning-Strickler
La frmula de Manning es
(3.8)
donde
V
velocidad media, en m/s;
Rh
radio hidrulico, en m;
Sf
pendiente de friccin;
factor de rugosidad, depende del material de la tubera.
Si se despeja la pendiente de friccin de la ecuacin 3.8, resulta

(
(3.9)
Si Sf representa la relacin entre la prdida de energa y la longitud del tubo en que esta ocurre
(3.10)
De donde
(3.11)
A partir de estas definiciones, en la tabla 3.6 se presentan los resultados obtenidos con la frmula de
Manning-Strickler.

67

Tabla 3.6. Prdidas por friccin con la frmula de Manning-Strickler
L
(m)
Sf
hf
(m)
Tramo 1
0.88
0.01
0.035
0.031
Tramo 2
1.1
0.01
0.035
0.039
Tramo 3
1.53
0.01
0.035
0.054
Tramo 4
0.45
0.01
0.035
0.016
Tramo 5
0.34
0.01
0.035
0.012
Total
4.3
0.152
Frmula de Hazen-Williams
[
(3.12)
donde
hf
prdida por friccin;
depende del material del tubo, (tabla 3.7);
dimetro de la tubera, en m;
longitud de la tubera, en m.

68
Tabla 3.7. Valores de CH, , m, n, y N de acuerdo con el material de que est construido el tubo
Material
Acero corrugado
Acero con juntas lock-bar (nuevo)
Acero galvanizado (nuevo y usado)
Acero remachado (nuevo)
Acero remachado (usado)
Acero soldado o con remache avellanado y embutido (nuevo)
Acero soldado o con remache avellanado y embutido (usado)
Acero sin costura (nuevo)
Acero sin costura (usado)
Acero soldado, con revestimiento especial (nuevo y usado)
Fierro fundido limpio (nuevo)
Fierro fundido, sin incrustaciones (usado)
Fierro fundido, con incrustaciones (viejo)
Plstico
Asbesto cemento (nuevo)
Cobre y latn
Conductos con acabado interior de cemento pulido
Concreto, acabado liso
Concreto, acabado comn
Concreto monoltico, colado con cimbras deslizantes (D>1.25 m)
Concreto monoltico bien cimbrado y pulido (D>1.25 m)
Concreto monoltico bien cimbrado y sin pulir (D>1.25 m)
Concreto con acabado tosco (D>1.25 m)
Concreto con juntas de macho y campana (D>0.8 m)
Concreto con juntas toscas (D>0.5 m)
Concreto con juntas toscas (D<0.5 m)
Conductos para alcantarillado
Tubos de barro vitrificado (drenes)
Tneles perforados en roca sin revestimiento
Madera cepillada o en duelas
CH
60
135
125
110
85
120
90
130
130
110
90
150
135
130
100
130
120
110
120
0.10
0.16
0.23
0.36
0.06
0.18
-
m
0.25
0.35
0.25
0.275
0.35
0.10
0.20
0.10
N
0.014
0.015 a 0.016
0.012 a 0.013
0.013
0.010 a 0.011
0.011 a 0.0123
0.014 a 0.015
0.015 a 0.017
0.0105 a 0.012
0.0125 a 0.014
0.014 a 0.017
0.011
0.025 a 0.040
0.0105 a 0.012
N
31
28 a 26
34
31 a 27
38
36
35
30
38
27 a 26
30
28
34
-
PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA ENSEANZA TERICO PRCTICA DEL FENMENO DE GOLPE
DE ARIETE
69
En la tabla 3.8 se muestran los resultados obtenidos con la frmula de Hazen-Williams.

Tabla 3.8. Prdidas por friccin con la frmula de Hazen-Williams
L
(m)
hf
(m)
Tramo 1
0.88
130
0.0250
Tramo 2
1.1
130
0.0313
Tramo 3
1.53
130
0.0435
Tramo 4
0.45
130
0.0128
Tramo 5
0.34
130
0.0097
Total
4.3
0.1223
A continuacin se muestra una tabla resumen con los resultados obtenidos para las prdidas por
friccin calculadas con los diferentes criterios.
Tabla 3.9. Tabla resumen. Prdidas por friccin
Prdidas por friccin hf, en m

Darcy-Weisbach
Manning
Hazen-Williams
3.7
Tramo 1
0.0348
0.0311
0.0250
Tramo 2
0.0435
0.0388
0.0313
Tramo 3
0.0606
0.0540
0.0435
Tramo 4
0.0178
0.0159
0.0128
Tramo 5
0.0135
0.0120
0.0097
Total
0.1702
0.1518
0.1223
Prdidas locales
Las tuberas que se utilizan estn compuestas, generalmente por tramos rectos y curvos para
ajustarse a los accidentes topogrficos del terreno, as como a los cambios que se presentan en la
geometra de la seccin y de los distintos dispositivos para el control de las descargas (accesorios).
Estos cambios originan prdidas de energa, distintas a las de friccin, localizadas en el sitio mismo
del cambio de geometra o de la alteracin del flujo. Tal tipo de prdida se conoce como prdida
local o por accesorio. Su magnitud se expresa como una fraccin de la carga de velocidad,
70
inmediatamente aguas abajo del sitio donde se produjo la prdida; la expresin general de prdida
local es
(3.13)
donde
h
prdida de energa, en m,
coeficiente sin dimensiones que depende del tipo de accesorio que se trate, del nmero de
Reynolds y de la rugosidad del tubo,
V2/2g la carga de velocidad, aguas abajo, de la zona de alteracin del flujo (salvo aclaracin en
contrario) en m.
En los siguientes apartados, se presentan los valores del coeficiente K, de acuerdo con el tipo de
perturbacin que se presenta en el dispositivo del Laboratorio.
a) Prdida por entrada
A la entrada de las tuberas se produce una prdida por el efecto de contraccin que sufre la vena
lquida y la formacin de zonas de separacin; el coeficiente K depende, principalmente, de la
brusquedad con que se efecta la contraccin del chorro.
En la figura 3.7 se observan diferentes formas de entrada a tuberas, sin embargo, observando la
entrada a la tubera del dispositivo experimental (fotografa 3.2), se asemeja a la que se muestra en
la figura 3.7a.

71
Figura 3.7. Coeficientes de prdida -por entrada- para diferentes formas
Fotografa 3.2. Entrada a la tubera del dispositivo experimental

72
De esta forma, el valor de K por entrada, es

K = 0.5
b) Prdidas por bifurcacin
La prdida de energa en una bifurcacin de conductos depende adems del ngulo que forma la
tubera secundaria con la principal- de la relacin entre los dimetros de ambas tuberas y de la
direccin de la corriente. Dicha prdida es mayor en la unin que en la bifurcacin y se expresa
como un porcentaje de la carga de velocidad, lo que demuestra que el coeficiente K es
independiente del nmero de Reynolds.
En el dispositivo experimental, se encuentra instalada una T bridada de 2.5 pulgadas de dimetro
(elemento 7, figura 3.1), la cual origina una prdida local en el ltimo tramo de la tubera.
Para determinar el valor del coeficiente se calcul la carga total antes y despus del accesorio en las
condiciones de flujo permanente y se obtuvo experimentalmente un valor K = 1.7.
Con los coeficientes por entrada y por bifurcacin, se determinaron las prdidas locales
correspondientes, en la siguiente tabla 3.10 se muestran los resultados.
Tabla 3.10. Prdidas locales
V2/2g
(m/s)
(m)
1.1872
0.0718
Coefic.
entrada
K
0.5
Coefic.
bifurcacin
K
1.7
Prdida a la
entrada
he
Prdida por
bifurcacin
hb
(m)
(m)
0.0359
0.1221
Al sumar las prdidas por friccin y las prdidas locales se obtuvo la prdida total en la tubera. En
seguida se muestran los resultados correspondientes a cada una de las frmulas utilizadas.
Utilizando las prdidas por friccin calculadas con la frmula de Darcy-Weisbach y tomando en
cuenta la carga de velocidad, se obtuvieron las prdidas totales y las cargas piezomtricas
correspondientes, cuyos valores se muestran en las tablas 3.11 y 3.12.

73
Tabla 3.11. Prdida total con la frmula de Darcy-Weisbach
Prdida total hr, en m

0.0707
Tramo 1
Tramo 2
0.1143
Tramo 3
0.1748
Tramo 4
0.3147
Tramo 5
0.3282
Tabla 3.12. Carga piezomtrica (Frmula de Darcy-Weisbach)
Carga Piezomtrica
hpz3
hpz4
hpz1
hpz2
Hov
(m)
(m)
(m)
(m)
(m)
Calculado
0.857
0.814
0.753
0.613
0.600
Medido
0.86
0.81
0.755
0.61
Utilizando las prdidas por friccin calculadas con la frmula de Manning y de igual manera
tomando en cuenta la carga de velocidad, se obtuvieron los siguientes resultados
Tabla 3.13. Prdida total con la frmula de Manning

0.0670
Tramo 1
Tramo 2
0.1058
Tramo 3
0.1599
Tramo 4
0.2978
Tramo 5
0.3098
Tabla 3.14. Carga piezomtrica (Frmula de Manning)
hpz1
hpz2
Carga Piezomtrica
hpz3
hpz4
Hov
(m)
(m)
(m)
(m)
(m)
Calculado
0.861
0.822
0.768
0.630
0.618
Medido
0.86
0.81
0.755
0.61

74
Por ltimo, los resultados obtenidos con la frmula de Hazen-Williams, se presentan a continuacin
Tabla 3.15. Prdida total con la frmula de Hazen-Williams

0.0609
Tramo 1
Tramo 2
0.0922
Tramo 3
0.1358
Tramo 4
0.2706
Tramo 5
0.2803
Tabla 3.16. Carga piezomtrica (Frmula de Hazen-Williams)
Carga Piezomtrica
hpz3
hpz4
hpz1
hpz2
Hov
(m)
(m)
(m)
(m)
(m)
Calculado
0.867
0.836
0.792
0.658
0.648
Medido
0.86
0.81
0.755
0.61
De los resultados mostrados en las tablas 3.9 a 3.14, se puede observar que con la frmula de
Darcy-Weisbach y con el ajuste de los coeficientes de prdidas locales, se obtuvieron cargas
piezomtricas muy cercanas a las medidas en el dispositivo experimental.
En las figuras 3.8a, 3.8b y 3.8c se muestra la lnea piezomtrica medida y calculada con las
diferentes frmulas.
Ho
Hov
TRAMO 1
TRAMO 2
pz 1
TRAMO 3
pz 2
TRAMO 4
pz 3
TRAMO 5
pz 4
vlvula
Lnea de cargas piezomtricas medida

Lnea de cargas piezomtricas calculada
Figura 3.8a. Lnea de cargas piezomtricas (Frmula de Darcy-Weisbach)

75
Ho
Hov
TRAMO 1
TRAMO 2
pz 1
TRAMO 3
pz 2
TRAMO 4
pz 3
TRAMO 5
pz 4
vlvula

Figura 3.8b. Lnea de cargas piezomtricas (Frmula de Manning)
Ho
Hov
TRAMO 1
TRAMO 2
pz 1
TRAMO 3
pz 2
TRAMO 4
pz 3
TRAMO 5
pz 4
vlvula

Figura 3.8c. Lnea de cargas piezomtricas (Frmula de Hazen-Williams)
3.8
Ley para maniobra de cierre no uniforme
La ley de cierre se determin de forma experimental. Esta ley de cierre no uniforme consiste en una
relacin entre el gasto medido con la vlvula totalmente abierta y el gasto medido para un cierre
parcial de la vlvula y luego asociarla a un tiempo de cierre.
Para determinar la relacin entre gastos se implement un dispositivo de apoyo para medir el
ngulo de cierre de la vlvula. Este dispositivo de apoyo consisti en colocar sobre la vlvula una

76
base de acrlico con una graduacin de 0 a 90, ya que ste es el intervalo en el que el maneral de la
vlvula alcanza el cierre o la apertura total.
En el vertedor triangular de aforos se registraron los siguientes datos, de acuerdo al ngulo de
cierre.
Tabla 3.17. Mediciones en vertedor triangular
Cierre
()
Nc
(m)
Nsa
(m)
0.2479
0.3809
10
0.2479
0.3709
20
0.2479
0.3568
30
0.2479
0.3403
40
0.2479
0.3244
50
0.2479
0.3089
60
0.2479
0.2939
70
0.2479
0.2768
80
0.2479
0.2479
Con los datos de la tabla 3.17, se determin el gasto correspondiente a cada ngulo de cierre. El
procedimiento de clculo para determinar el gasto fue el mismo que se utilizo al calcular el gasto en
flujo permanente. Los resultados se muestran a continuacin, en la tabla 3.18.
Tabla 3.18. Gasto ngulo de cierre
Cierre
()
(m)
B
() (m)
1/2
( m /s )
( m /s )
( m3/s )
0.13
9.81
45
1.7
0.613 12.782
0.5998
0.0039
10
0.12
9.81
45
1.7
0.618 13.821
0.6047
0.0032
20
0.11
9.81
45
1.7
0.619 15.611
0.6057
0.0024
30
0.09
9.81
45
1.7
0.624 18.398
0.6106
0.0016
40
0.08
9.81
45
1.7
0.628 22.222
0.6145
0.0010
50
0.06
9.81
45
1.7
0.644 27.869
0.6302
0.0006
60
0.05
9.81
45
1.7
0.654 36.957
0.6400
0.0003
70
0.03
9.81
45
1.7
0.676 58.824
0.6615
0.0001
80
0.00
9.81
45
1.7
B/h

77
Con la medicin de estos gastos, se realiz una curva de gastos en funcin del ngulo de cierre de la
vlvula (figura 3.9). En seguida, esta curva se ajust a una ecuacin general.
(3.14)
Ley de Cierre
Vlvula de esfera
Q (m3/s)
0.0045
0.0040
0.0035
Ley de Cierre
0.0030
Ley de cierre ajuste
0.0025
0.0020
0.0015
0.0010
0.0005
0.0000
10
20
30
40
50
60
70
80
90
ngulo de cierre ( )
Figura 3.9. Curva de gastos en funcin del ngulo de cierre de la vlvula
La finalidad de la ecuacin 3.14, es conocer el gasto para cualquier ngulo de cierre, aunque cabe
destacar que el cierre total de la vlvula se produjo entre 80 y 82. Utilizando esta ecuacin de
ajuste, se recalcul el gasto en funcin del ngulo de cierre y los resultados se muestran en la tabla
3.19.
Tabla 3.19. Gastos de acuerdo a ecuacin de ajuste
()
Q
(m3/s)
0.00398
10
0.00308
20
0.00229
30
0.00162
40
0.00106
50
0.00062
60
0.00029
70
0.00008
80

78
Con los gastos de la tabla 3.19, se obtuvo la relacin
, donde Q0 es el gasto que pasa a
travs de la vlvula totalmente abierta. Esta relacin se observa en la tabla 3.20.

Tabla 3.20. Relacin
()
1.0000
10
0.7728
20
0.5746
30
0.4055
40
0.2655
50
0.1545
60
0.0727
70
0.0199
80
0.0000
Cuando la ley de cierre no es uniforme, se tiene una variacin de con respecto al tiempo de cierre
(tc). Por otra parte, se tiene la experiencia de utilizar un tiempo de cierre de un segundo cuando se
realiza la prctica Golpe de Ariete utilizando como mtodo de solucin a las cadenas de Allievi. Sin
embargo, es importante mencionar que la ley de cierre utilizada en dicha prctica es una ley de
cierre lineal. En este caso se utilizaron tiempos de cierre diferentes, con el fin de observar la
relacin entre las cargas piezomtricas cuando se presenta el flujo transitorio y el tiempo de cierre
correspondiente.
Los tiempos de cierre utilizados fueron tc = 1.0 segundos, tc = 1.2 segundos y tc = 1.5 segundos. La
ley de cierre para cada uno de estos tiempos, con respecto al ngulo de cierre, se muestra en la
siguiente tabla.

79
Tabla 3.21. Ley de cierre
()
tc = 1.0
(s)
tc = 1.2
(s)
tc = 1.5
(s)
1.0000
10
0.7728
0.1250
0.1500
0.1875
20
0.5746
0.2500
0.3000
0.3750
30
0.4055
0.3750
0.4500
0.5625
40
0.2655
0.5000
0.6000
0.7500
50
0.1545
0.6250
0.7500
0.9375
60
0.0727
0.7500
0.9000
1.1250
70
0.0199
0.8750
1.0500
1.3125
80
0.0000
1.2
1.5
Por ltimo se realiz la grfica -tc y se hizo un ajuste a la curva, obteniendo la ecuacin de la ley de
cierre no uniforme. Estos grficos y su correspondiente ecuacin de ajuste se muestran a
continuacin.

80
Tabla 3.22. Ley de cierre no uniforme
Ley de cierre no uniforme
Ecuacin
Ley de cierre
tc=1 s
t
1.00
0.80
0.60
0.40
0.20
0.00
0.00
0.20
0.40
0.60
0.80
1.00
t (s)
Ley de cierre
tc=1.2 s
t
1.00
0.80
0.60
0.40
0.20
0.00
0.00
0.20
0.40
0.60
0.80
1.00
1.20
t (s)
Ley de cierre
tc=1.5 s
t
1.00
0.80
0.60
0.40
0.20
0.00
0.00
0.30
0.60
0.90
1.20
1.50
t (s)
Para determinar la relacin entre el rea de la vlvula totalmente abierta y el rea para cada ngulo
de cierre, se hizo una medicin en un dispositivo como el mostrado en la fotografa 3.3.

81
Fotografa 3.3. Dispositivo de medicin de apertura de vlvula
El procedimiento consisti en tomar fotografas del rea de la vlvula que se iba obturando para
cada ngulo de cierre. Se realiz en intervalos de 10. En las siguientes fotografas se puede
observar ms claramente.
Aunque la relacin entre reas que finalmente se obtuvo no se utiliz directamente para determinar
la ley de cierre, puede funcionar para determinar coeficientes de prdida para vlvulas de esfera. La
tabla que a continuacin se muestra es parte de una recopilacin de informacin bibliogrfica y
contiene coeficientes de prdida para vlvulas esfricas en funcin de la relacin de reas.
Tabla 3.23. Coeficientes de prdida para vlvulas esfricas
K
0.05
A/A0
0.926
10
0.29
0.85
15
0.75
0.772
20
1.56
0.692
25
3.1
0.613
30
5.17
0.535
35
9.68
0.458
40
17.3
0.358
45
31.2
0.315
50
52.6
0.25
55
106
0.19
60
206
0.137
65
486
0.091
82

82
Fotografa 3.4. Totalmente abierta
Fotografa 3.5. Cierre 10

83
Se trasladaron las fotografas al programa Autocad y se determinaron las reas correspondientes

(figura 3.10), mismas que se muestran en la tabla 3.24.
A0
A
Q0
Q
a) Vlvula totalmente abierta (A0)
b) Vlvula con cierre (A)
Figura 3.10. Vlvula totalmente abierta y cierre parcial

84
Tabla 3.24. Relacin de reas
()
A
( m2 )
A/A0
0.0032
1.0000
10
0.0027
0.8551
20
0.0021
0.6684
30
0.0018
0.5693
40
0.0014
0.4466
50
0.0009
0.2722
60
0.0005
0.1595
70
0.0002
0.0529
80
0.0000
0.0000
Al comparar las tablas 3.21 y 3.22, se observ que varan en un pequeo porcentaje las relaciones
A/A0, lo que permite establecer que los datos experimentales que se obtuvieron para la vlvula del
caso de estudio son adecuados.
Una vez determinada la ley de cierre no uniforme, se dio solucin al caso de estudio aplicando dos
procedimientos: utilizando las ecuaciones de Allievi y el mtodo de las caractersticas. Cabe
mencionar que al utilizar las ecuaciones de Allievi, se utiliz una ley de cierre lineal y al utilizar el
mtodo de las caractersticas se us la ley de cierre no uniforme obtenida. Es por ello la importancia
de utilizar el mtodo de las caractersticas, ya que da la oportunidad de analizar el fenmeno del
golpe de ariete de manera ms cercana a la realidad.
3.9
Solucin aplicando el mtodo de Allievi
A continuacin se presenta la solucin para el caso de estudio mediante el mtodo de las cadenas de
Allievi. Anteriormente se determinaron las condiciones iniciales, el valor de la celeridad de la onda,
as como el periodo de la onda de presin. En seguida se presentan algunos clculos previos
tambin necesarios para aplicar este mtodo.
Se calcul la relacin entre el tiempo de maniobra tc y el periodo de la onda de presin T,
considerando que el tiempo de la maniobra fue de un segundo.

85
(3.18)
Lo que resulta
Se debe cerrar al nmero entero siguiente. Por lo tanto = 134

Posteriormente, el objetivo fue calcular la carga de presin Hi y la sobrepresin Hi en m, en la
seccin cercana a la vlvula, aplicando las ecuaciones de Allievi y considerando una ley de cierre
lineal, con la ecuacin
donde:
Hi
carga de presin correspondiente a cada instante relativo
del cierre en la vlvula, en
m;
H0
carga de presin medida en la seccin 0 con la vlvula cerrada, en m;
Hi
sobrepresin en cada instante relativo del cierre, en m;
relacin de cargas en cada instante relativo del cierre; se obtiene con la expresin de
Allievi.
Con estos antecedentes, se plante la solucin mediante la siguiente forma.

Como se explic en el captulo 2, la ecuacin 2.82 es la ecuacin propuesta por Allievi, misma que
reescribiremos como
(
(3.19)
Para su solucin en cadena, es ms conveniente despejar a la relacin de cargas i.

(
(3.20)
donde:
i-1
relacin de cargas en el instante inicial del intervalo de anlisis, segn la ley de cierre;
86
relacin de cargas en el instante final del periodo T, segn la ley de cierre. Para el primer
instante
Constante de Allievi
g
V0
velocidad a la salida de la tubera antes del cierre, en m/s;
ley lineal de cierre lento que se muestra en la figura 3.11 y que tiene por ecuacin
1.0
=1-i/
tc
Figura 3.11. Ley de cierre lineal
3.9.1
Resultados
Siguiendo el procedimiento anterior, se obtuvo la variacin de la carga de presin con respecto al

tiempo para los instantes i=0, i=1, i=2 e i=3, correspondientes a un tiempo de cierre de un segundo.
Para i = 0

87
Para i=1
(
Sustituyendo valores,
(
)( ))
Para i=2
(
(
)(
))

88
Para i=3
(
(
)(
))
En seguida se muestra la totalidad de los resultados de la variacin de la carga de presin para los
tiempos de cierre de 1 s, 1.2 s y 1.5 s en las tablas 3.25a, 3.25c y 3.25e y su representacin grfica
en las figuras 3.12a, 3.12b y 3.12c, respectivamente.

89
Tabla 3.25a. Variacin de H (tc = 1.0 s)
(m)
(m)
(m)
0.00748
0.99254
1.00741
1.01488
0.01488
0.01495
0.98507
1.01472
1.02965
0.02965
0.02243
0.97761
1.02192
1.04432
0.04432
0.02991
0.97015
1.02902
1.05888
0.05888
0.03739
0.96269
1.03601
1.07331
0.07331
0.04486
0.95522
1.04289
1.08762
0.08762
0.05234
0.94776
1.04966
1.10180
0.10180
0.05982
0.94030
1.05633
1.11584
0.11584
0.06729
0.93284
1.06289
1.12974
0.12974
10
0.07477
0.92537
1.06934
1.14349
0.14349
11
0.08225
0.91791
1.07568
1.15709
0.15709
12
0.08973
0.91045
1.08191
1.17053
0.17053
13
0.09720
0.90299
1.08803
1.18381
0.18381
14
0.10468
0.89552
1.09404
1.19693
0.19693
15
0.11216
0.88806
1.09995
1.20988
0.20988
16
0.11963
0.88060
1.10574
1.22265
0.22265
17
0.12711
0.87313
1.11142
1.23525
0.23525
18
0.13459
0.86567
1.11699
1.24767
0.24767
19
0.14207
0.85821
1.12245
1.25990
0.25990
20
0.14954
0.85075
1.12781
1.27194
0.27194
21
0.15702
0.84328
1.13305
1.28380
0.28380
22
0.16450
0.83582
1.13818
1.29546
0.29546
23
0.17197
0.82836
1.14321
1.30693
0.30693
24
0.17945
0.82090
1.14813
1.31821
0.31821
25
0.18693
0.81343
1.15294
1.32928
0.32928
26
0.19441
0.80597
1.15765
1.34015
0.34015
27
0.20188
0.79851
1.16225
1.35082
0.35082

90
Tabla 3.25a. Variacin de H (tc = 1.0 s) (Continuacin)
(m)
(m)
(m)
28
0.20936
0.79104
1.16674
1.36129
0.36129
29
0.21684
0.78358
1.17113
1.37156
0.37156
30
0.22431
0.77612
1.17542
1.38161
0.38161
31
0.23179
0.76866
1.17960
1.39147
0.39147
32
0.23927
0.76119
1.18369
1.40111
0.40111
33
0.24675
0.75373
1.18767
1.41055
0.41055
34
0.25422
0.74627
1.19155
1.41978
0.41978
35
0.26170
0.73881
1.19533
1.42880
0.42880
36
0.26918
0.73134
1.19901
1.43762
0.43762
37
0.27665
0.72388
1.20259
1.44623
0.44623
38
0.28413
0.71642
1.20608
1.45464
0.45464
39
0.29161
0.70896
1.20948
1.46284
0.46284
40
0.29909
0.70149
1.21278
1.47083
0.47083
41
0.30656
0.69403
1.21599
1.47863
0.47863
42
0.31404
0.68657
1.21910
1.48622
0.48622
43
0.32152
0.67910
1.22213
1.49361
0.49361
44
0.32899
0.67164
1.22507
1.50080
0.50080
45
0.33647
0.66418
1.22792
1.50779
0.50779
46
0.34395
0.65672
1.23069
1.51459
0.51459
47
0.35143
0.64925
1.23337
1.52119
0.52119
48
0.35890
0.64179
1.23596
1.52760
0.52760
49
0.36638
0.63433
1.23848
1.53382
0.53382
50
0.37386
0.62687
1.24091
1.53986
0.53986
51
0.38133
0.61940
1.24326
1.54571
0.54571
52
0.38881
0.61194
1.24554
1.55137
0.55137
53
0.39629
0.60448
1.24774
1.55686
0.55686
54
0.40377
0.59701
1.24987
1.56217
0.56217

91
(m)
(m)
(m)
55
0.41124
0.58955
1.25192
1.56730
0.56730
56
0.41872
0.58209
1.25390
1.57226
0.57226
57
0.42620
0.57463
1.25581
1.57706
0.57706
58
0.43367
0.56716
1.25765
1.58168
0.58168
59
0.44115
0.55970
1.25942
1.58614
0.58614
60
0.44863
0.55224
1.26113
1.59044
0.59044
61
0.45611
0.54478
1.26277
1.59459
0.59459
62
0.46358
0.53731
1.26435
1.59858
0.59858
63
0.47106
0.52985
1.26587
1.60241
0.60241
64
0.47854
0.52239
1.26732
1.60611
0.60611
65
0.48601
0.51493
1.26872
1.60965
0.60965
66
0.49349
0.50746
1.27006
1.61305
0.61305
67
0.50097
0.50000
1.27135
1.61632
0.61632
68
0.50845
0.49254
1.27258
1.61945
0.61945
69
0.51592
0.48507
1.27375
1.62245
0.62245
70
0.52340
0.47761
1.27488
1.62532
0.62532
71
0.53088
0.47015
1.27596
1.62806
0.62806
72
0.53835
0.46269
1.27698
1.63068
0.63068
73
0.54583
0.45522
1.27796
1.63319
0.63319
74
0.55331
0.44776
1.27890
1.63558
0.63558
75
0.56079
0.44030
1.27979
1.63785
0.63785
76
0.56826
0.43284
1.28063
1.64002
0.64002
77
0.57574
0.42537
1.28144
1.64208
0.64208
78
0.58322
0.41791
1.28220
1.64404
0.64404
79
0.59069
0.41045
1.28293
1.64590
0.64590
80
0.59817
0.40299
1.28361
1.64767
0.64767
81
0.60565
0.39552
1.28427
1.64934
0.64934

92
(m)
(m)
(m)
82
0.61313
0.38806
1.28488
1.65092
0.65092
83
0.62060
0.38060
1.28546
1.65242
0.65242
84
0.62808
0.37313
1.28601
1.65383
0.65383
85
0.63556
0.36567
1.28653
1.65516
0.65516
86
0.64303
0.35821
1.28702
1.65641
0.65641
87
0.65051
0.35075
1.28748
1.65759
0.65759
88
0.65799
0.34328
1.28791
1.65870
0.65870
91
0.68042
0.32090
1.28904
1.66162
0.66162
92
0.68790
0.31343
1.28937
1.66247
0.66247
93
0.69537
0.30597
1.28968
1.66326
0.66326
94
0.70285
0.29851
1.28996
1.66400
0.66400
95
0.71033
0.29104
1.29023
1.66468
0.66468
96
0.71781
0.28358
1.29047
1.66532
0.66532
97
0.72528
0.27612
1.29070
1.66590
0.66590
98
0.73276
0.26866
1.29091
1.66645
0.66645
99
0.74024
0.26119
1.29110
1.66694
0.66694
100
0.74771
0.25373
1.29128
1.66740
0.66740
101
0.75519
0.24627
1.29144
1.66782
0.66782
102
0.76267
0.23881
1.29159
1.66820
0.66820
103
0.77015
0.23134
1.29172
1.66855
0.66855
104
0.77762
0.22388
1.29185
1.66887
0.66887
105
0.78510
0.21642
1.29196
1.66916
0.66916
106
0.79258
0.20896
1.29206
1.66942
0.66942
107
0.80005
0.20149
1.29215
1.66965
0.66965
108
0.80753
0.19403
1.29223
1.66986
0.66986

93
(m)
(m)
(m)
109
0.81501
0.18657
1.29230
1.67005
0.67005
110
0.82249
0.17910
1.29237
1.67021
0.67021
111
0.82996
0.17164
1.29242
1.67036
0.67036
112
0.83744
0.16418
1.29247
1.67048
0.67048
113
0.84492
0.15672
1.29252
1.67060
0.67060
114
0.85239
0.14925
1.29255
1.67069
0.67069
115
0.85987
0.14179
1.29259
1.67078
0.67078
116
0.86735
0.13433
1.29261
1.67085
0.67085
117
0.87483
0.12687
1.29264
1.67091
0.67091
118
0.88230
0.11940
1.29266
1.67096
0.67096
119
0.88978
0.11194
1.29267
1.67100
0.67100
120
0.89726
0.10448
1.29269
1.67104
0.67104
121
0.90473
0.09701
1.29270
1.67107
0.67107
122
0.91221
0.08955
1.29271
1.67109
0.67109
123
0.91969
0.08209
1.29271
1.67110
0.67110
124
0.92717
0.07463
1.29272
1.67112
0.67112
125
0.93464
0.06716
1.29272
1.67113
0.67113
126
0.94212
0.05970
1.29272
1.67113
0.67113
127
0.94960
0.05224
1.29273
1.67114
0.67114
128
0.95707
0.04478
1.29273
1.67114
0.67114
129
0.96455
0.03731
1.29273
1.67114
0.67114
130
0.97203
0.02985
1.29273
1.67115
0.67115
131
0.97951
0.02239
1.29273
1.67115
0.67115
132
0.98698
0.01493
1.29273
1.67115
0.67115
133
0.99446
0.00746
1.29273
1.67115
0.67115
134
1.00194
0.00000
1.29273
1.67115
0.67115

94
t-H
tc = 1.0 s
H (m)
2.00
1.80
1.60
1.40
1.20
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.80
0.0
1.00
tc (s)
Figura 3.12a.Variacin de H (tc=1.0 s)
Finalmente, de la solucin obtenida con las ecuaciones de Allievi, observamos que la carga mxima
que se present fue
y la sobrepresin mxima
Con la medicin realizada se tiene
Comparando la carga mxima y la sobrepresin mxima calculada con la medida, obtuvimos los
siguientes porcentajes de error.
Tabla 3.25b. Porcentaje de error (tc=1.0 s)
Hmx
Allievi
1.671
Medido
1.5
% error
11%
0.671
0.89
25%

95
Resultados con tc = 1.2 segundos y tc = 1.5 segundos

Para tc = 1.2 s
Tabla 3.25c. Variacin de H (tc = 1.2 s)
(m)
(m)
(m)
0.00748
0.99379
1.00616
1.01236
0.01236
0.01495
0.98758
1.01222
1.02459
0.02459
0.02243
0.98137
1.01818
1.03669
0.03669
0.02991
0.97516
1.02403
1.04865
0.04865
0.03739
0.96894
1.02979
1.06047
0.06047
0.04486
0.96273
1.03544
1.07214
0.07214
0.05234
0.95652
1.04100
1.08367
0.08367
0.05982
0.95031
1.04645
1.09506
0.09506
0.06729
0.94410
1.05180
1.10629
0.10629
10
0.07477
0.93789
1.05705
1.11736
0.11736
11
0.08225
0.93168
1.06221
1.12828
0.12828
12
0.08973
0.92547
1.06726
1.13904
0.13904
13
0.09720
0.91925
1.07221
1.14964
0.14964
14
0.10468
0.91304
1.07707
1.16008
0.16008
15
0.11216
0.90683
1.08183
1.17035
0.17035
16
0.11963
0.90062
1.08649
1.18046
0.18046
17
0.12711
0.89441
1.09106
1.19040
0.19040
18
0.13459
0.88820
1.09553
1.20018
0.20018
19
0.14207
0.88199
1.09990
1.20978
0.20978
20
0.14954
0.87578
1.10418
1.21922
0.21922
21
0.15702
0.86957
1.10837
1.22848
0.22848
22
0.16450
0.86335
1.11246
1.23758
0.23758
23
0.17197
0.85714
1.11647
1.24650
0.24650
24
0.17945
0.85093
1.12038
1.25525
0.25525
25
0.18693
0.84472
1.12420
1.26383
0.26383
26
0.19441
0.83851
1.12794
1.27224
0.27224
27
0.20188
0.83230
1.13158
1.28048
0.28048

96
Tabla 3.25c. Variacin de H (tc = 1.2 s) (Continuacin)
(m)
(m)
(m)
28
0.20936
0.82609
1.13514
1.28854
0.28854
29
0.21684
0.81988
1.13861
1.29644
0.29644
30
0.22431
0.81366
1.14200
1.30417
0.30417
31
0.23179
0.80745
1.14530
1.31172
0.31172
32
0.23927
0.80124
1.14853
1.31911
0.31911
33
0.24675
0.79503
1.15166
1.32633
0.32633
34
0.25422
0.78882
1.15472
1.33339
0.33339
35
0.26170
0.78261
1.15770
1.34027
0.34027
36
0.26918
0.77640
1.16060
1.34700
0.34700
37
0.27665
0.77019
1.16343
1.35356
0.35356
38
0.28413
0.76398
1.16617
1.35996
0.35996
39
0.29161
0.75776
1.16885
1.36620
0.36620
40
0.29909
0.75155
1.17145
1.37229
0.37229
41
0.30656
0.74534
1.17397
1.37821
0.37821
42
0.31404
0.73913
1.17643
1.38398
0.38398
43
0.32152
0.73292
1.17881
1.38960
0.38960
44
0.32899
0.72671
1.18113
1.39507
0.39507
45
0.33647
0.72050
1.18338
1.40039
0.40039
46
0.34395
0.71429
1.18556
1.40556
0.40556
47
0.35143
0.70807
1.18768
1.41058
0.41058
48
0.35890
0.70186
1.18973
1.41546
0.41546
49
0.36638
0.69565
1.19172
1.42021
0.42021
50
0.37386
0.68944
1.19365
1.42481
0.42481
51
0.38133
0.68323
1.19552
1.42927
0.42927
52
0.38881
0.67702
1.19733
1.43360
0.43360
53
0.39629
0.67081
1.19908
1.43780
0.43780
54
0.40377
0.66460
1.20078
1.44187
0.44187

97
(m)
(m)
(m)
55
0.41124
0.65839
1.20242
1.44581
0.44581
56
0.41872
0.65217
1.20400
1.44963
0.44963
57
0.42620
0.64596
1.20554
1.45332
0.45332
58
0.43367
0.63975
1.20702
1.45689
0.45689
59
0.44115
0.63354
1.20845
1.46035
0.46035
60
0.44863
0.62733
1.20983
1.46369
0.46369
61
0.45611
0.62112
1.21116
1.46691
0.46691
62
0.46358
0.61491
1.21245
1.47003
0.47003
63
0.47106
0.60870
1.21369
1.47303
0.47303
64
0.47854
0.60248
1.21488
1.47593
0.47593
65
0.48601
0.59627
1.21603
1.47873
0.47873
66
0.49349
0.59006
1.21714
1.48143
0.48143
67
0.50097
0.58385
1.21820
1.48402
0.48402
68
0.50845
0.57764
1.21923
1.48652
0.48652
69
0.51592
0.57143
1.22022
1.48893
0.48893
70
0.52340
0.56522
1.22116
1.49124
0.49124
71
0.53088
0.55901
1.22208
1.49347
0.49347
72
0.53835
0.55280
1.22295
1.49561
0.49561
73
0.54583
0.54658
1.22379
1.49766
0.49766
74
0.55331
0.54037
1.22459
1.49963
0.49963
75
0.56079
0.53416
1.22537
1.50152
0.50152
76
0.56826
0.52795
1.22611
1.50334
0.50334
77
0.57574
0.52174
1.22681
1.50507
0.50507
78
0.58322
0.51553
1.22749
1.50674
0.50674
79
0.59069
0.50932
1.22814
1.50833
0.50833
80
0.59817
0.50311
1.22876
1.50985
0.50985
81
0.60565
0.49689
1.22935
1.51131
0.51131

98
(m)
(m)
(m)
82
0.61313
0.49068
1.22992
1.51270
0.51270
83
0.62060
0.48447
1.23046
1.51403
0.51403
84
0.62808
0.47826
1.23097
1.51530
0.51530
85
0.63556
0.47205
1.23146
1.51650
0.51650
86
0.64303
0.46584
1.23193
1.51766
0.51766
87
0.65051
0.45963
1.23238
1.51875
0.51875
88
0.65799
0.45342
1.23280
1.51980
0.51980
89
0.66547
0.44720
1.23320
1.52079
0.52079
90
0.67294
0.44099
1.23358
1.52173
0.52173
91
0.68042
0.43478
1.23395
1.52263
0.52263
92
0.68790
0.42857
1.23429
1.52347
0.52347
93
0.69537
0.42236
1.23462
1.52428
0.52428
94
0.70285
0.41615
1.23493
1.52504
0.52504
95
0.71033
0.40994
1.23522
1.52576
0.52576
96
0.71781
0.40373
1.23549
1.52644
0.52644
97
0.72528
0.39752
1.23575
1.52709
0.52709
98
0.73276
0.39130
1.23600
1.52770
0.52770
99
0.74024
0.38509
1.23623
1.52827
0.52827
100
0.74771
0.37888
1.23645
1.52881
0.52881
101
0.75519
0.37267
1.23666
1.52932
0.52932
102
0.76267
0.36646
1.23685
1.52979
0.52979
103
0.77015
0.36025
1.23703
1.53024
0.53024
104
0.77762
0.35404
1.23720
1.53066
0.53066
105
0.78510
0.34783
1.23736
1.53106
0.53106
106
0.79258
0.34161
1.23751
1.53143
0.53143

99
(m)
(m)
(m)
107
0.80005
0.33540
1.23765
1.53177
0.53177
108
0.80753
0.32919
1.23778
1.53209
0.53209
109
0.81501
0.32298
1.23790
1.53239
0.53239
110
0.82249
0.31677
1.23801
1.53267
0.53267
111
0.82996
0.31056
1.23812
1.53293
0.53293
112
0.83744
0.30435
1.23821
1.53317
0.53317
113
0.84492
0.29814
1.23830
1.53339
0.53339
114
0.85239
0.29193
1.23839
1.53360
0.53360
115
0.85987
0.28571
1.23846
1.53379
0.53379
116
0.86735
0.27950
1.23853
1.53397
0.53397
117
0.87483
0.27329
1.23860
1.53413
0.53413
118
0.88230
0.26708
1.23866
1.53428
0.53428
119
0.88978
0.26087
1.23871
1.53441
0.53441
120
0.89726
0.25466
1.23876
1.53454
0.53454
121
0.90473
0.24845
1.23881
1.53465
0.53465
122
0.91221
0.24224
1.23885
1.53475
0.53475
123
0.91969
0.23602
1.23889
1.53485
0.53485
124
0.92717
0.22981
1.23892
1.53493
0.53493
125
0.93464
0.22360
1.23896
1.53501
0.53501
126
0.94212
0.21739
1.23898
1.53508
0.53508
127
0.94960
0.21118
1.23901
1.53514
0.53514
128
0.95707
0.20497
1.23903
1.53520
0.53520
129
0.96455
0.19876
1.23905
1.53525
0.53525
130
0.97203
0.19255
1.23907
1.53529
0.53529
131
0.97951
0.18634
1.23909
1.53533
0.53533
132
0.98698
0.18012
1.23910
1.53537
0.53537

100
(m)
(m)
(m)
133
0.99446
0.17391
1.23911
1.53540
0.53540
134
1.00194
0.16770
1.23912
1.53543
0.53543
135
1.00941
0.16149
1.23913
1.53545
0.53545
136
1.01689
0.15528
1.23914
1.53547
0.53547
137
1.02437
0.14907
1.23915
1.53549
0.53549
138
1.03185
0.14286
1.23915
1.53550
0.53550
139
1.03932
0.13665
1.23916
1.53551
0.53551
140
1.04680
0.13043
1.23916
1.53553
0.53553
141
1.05428
0.12422
1.23917
1.53553
0.53553
142
1.06175
0.11801
1.23917
1.53554
0.53554
143
1.06923
0.11180
1.23917
1.53555
0.53555
144
1.07671
0.10559
1.23917
1.53555
0.53555
145
1.08419
0.09938
1.23918
1.53556
0.53556
146
1.09166
0.09317
1.23918
1.53556
0.53556
147
1.09914
0.08696
1.23918
1.53556
0.53556
148
1.10662
0.08075
1.23918
1.53556
0.53556
149
1.11409
0.07453
1.23918
1.53557
0.53557
150
1.12157
0.06832
1.23918
1.53557
0.53557
151
1.12905
0.06211
1.23918
1.53557
0.53557
152
1.13653
0.05590
1.23918
1.53557
0.53557
153
1.14400
0.04969
1.23918
1.53557
0.53557
154
1.15148
0.04348
1.23918
1.53557
0.53557
155
1.15896
0.03727
1.23918
1.53557
0.53557
156
1.16643
0.03106
1.23918
1.53557
0.53557
157
1.17391
0.02484
1.23918
1.53557
0.53557
158
1.18139
0.01863
1.23918
1.53557
0.53557

101
(m)
(m)
(m)
159
1.18887
0.01242
1.23918
1.53557
0.53557
160
1.19634
0.00621
1.23918
1.53557
0.53557
161
1.20382
0.00000
1.23918
1.53557
0.53557
t-H
tc = 1.2 s
H (m)
1.60
1.40
1.20
1.3
1.2
1.1
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.80
0.0
1.00
tc (s)
Figura 3.12b.Variacin de H (tc=1.2 s)
Se observ que para un tiempo de cierre de 1.2 s, la carga mxima fue de Hmx = 1.536 m y la
sobrepresin mxima
Con la medicin realizada, sucedi lo siguiente
Tabla 3.25d. Porcentaje de error (tc=1.2 s)
Hmx
Allievi
1.536
Medido
1.5
% error
2%
0.536
0.89
40%

102
Para tc = 1.5 s
Tabla 3.25e. Variacin de H (tc = 1.5 s)
(m)
(m)
(m)
0.00748
0.99502
1.00493
1.00988
0.00988
0.01495
0.99005
1.00976
1.01962
0.01962
0.02243
0.98507
1.01451
1.02923
0.02923
0.02991
0.98010
1.01916
1.03869
0.03869
0.03739
0.97512
1.02372
1.04800
0.04800
0.04486
0.97015
1.02819
1.05718
0.05718
0.05234
0.96517
1.03257
1.06621
0.06621
0.05982
0.96020
1.03687
1.07509
0.07509
0.06729
0.95522
1.04107
1.08383
0.08383
10
0.07477
0.95025
1.04519
1.09242
0.09242
11
0.08225
0.94527
1.04922
1.10086
0.10086
12
0.08973
0.94030
1.05317
1.10916
0.10916
13
0.09720
0.93532
1.05703
1.11731
0.11731
14
0.10468
0.93035
1.06081
1.12531
0.12531
15
0.11216
0.92537
1.06450
1.13317
0.13317
16
0.11963
0.92040
1.06812
1.14088
0.14088
17
0.12711
0.91542
1.07165
1.14844
0.14844
18
0.13459
0.91045
1.07511
1.15585
0.15585
19
0.14207
0.90547
1.07848
1.16312
0.16312
20
0.14954
0.90050
1.08178
1.17025
0.17025
21
0.15702
0.89552
1.08500
1.17723
0.17723
22
0.16450
0.89055
1.08815
1.18407
0.18407
23
0.17197
0.88557
1.09122
1.19076
0.19076
24
0.17945
0.88060
1.09422
1.19731
0.19731
25
0.18693
0.87562
1.09714
1.20373
0.20373
26
0.19441
0.87065
1.10000
1.21000
0.21000
27
0.20188
0.86567
1.10278
1.21613
0.21613

103
Tabla 3.25e. Variacin de H (tc = 1.5 s) (Continuacin)
(m)
(m)
(m)
28
0.20936
0.86070
1.10550
1.22213
0.22213
29
0.21684
0.85572
1.10815
1.22799
0.22799
30
0.22431
0.85075
1.11073
1.23372
0.23372
31
0.23179
0.84577
1.11324
1.23931
0.23931
32
0.23927
0.84080
1.11569
1.24477
0.24477
33
0.24675
0.83582
1.11808
1.25010
0.25010
34
0.25422
0.83085
1.12040
1.25531
0.25531
35
0.26170
0.82587
1.12267
1.26038
0.26038
36
0.26918
0.82090
1.12487
1.26533
0.26533
37
0.27665
0.81592
1.12701
1.27016
0.27016
38
0.28413
0.81095
1.12910
1.27486
0.27486
39
0.29161
0.80597
1.13112
1.27944
0.27944
40
0.29909
0.80100
1.13310
1.28391
0.28391
41
0.30656
0.79602
1.13501
1.28826
0.28826
42
0.31404
0.79104
1.13688
1.29249
0.29249
43
0.32152
0.78607
1.13869
1.29661
0.29661
44
0.32899
0.78109
1.14045
1.30062
0.30062
45
0.33647
0.77612
1.14215
1.30452
0.30452
46
0.34395
0.77114
1.14381
1.30831
0.30831
47
0.35143
0.76617
1.14542
1.31199
0.31199
48
0.35890
0.76119
1.14698
1.31557
0.31557
49
0.36638
0.75622
1.14850
1.31905
0.31905
50
0.37386
0.75124
1.14997
1.32243
0.32243
51
0.38133
0.74627
1.15140
1.32571
0.32571
52
0.38881
0.74129
1.15278
1.32890
0.32890
53
0.39629
0.73632
1.15412
1.33199
0.33199
54
0.40377
0.73134
1.15542
1.33499
0.33499

104
(m)
(m)
(m)
55
0.41124
0.72637
1.15667
1.33789
0.33789
56
0.41872
0.72139
1.15789
1.34071
0.34071
57
0.42620
0.71642
1.15907
1.34344
0.34344
58
0.43367
0.71144
1.16021
1.34609
0.34609
59
0.44115
0.70647
1.16131
1.34865
0.34865
60
0.44863
0.70149
1.16238
1.35113
0.35113
61
0.45611
0.69652
1.16342
1.35354
0.35354
62
0.46358
0.69154
1.16441
1.35586
0.35586
63
0.47106
0.68657
1.16538
1.35811
0.35811
64
0.47854
0.68159
1.16631
1.36029
0.36029
65
0.48601
0.67662
1.16721
1.36239
0.36239
66
0.49349
0.67164
1.16808
1.36442
0.36442
67
0.50097
0.66667
1.16892
1.36638
0.36638
68
0.50845
0.66169
1.16974
1.36828
0.36828
69
0.51592
0.65672
1.17052
1.37011
0.37011
70
0.52340
0.65174
1.17127
1.37188
0.37188
71
0.53088
0.64677
1.17200
1.37358
0.37358
72
0.53835
0.64179
1.17270
1.37523
0.37523
73
0.54583
0.63682
1.17338
1.37681
0.37681
74
0.55331
0.63184
1.17403
1.37834
0.37834
75
0.56079
0.62687
1.17466
1.37982
0.37982
76
0.56826
0.62189
1.17526
1.38123
0.38123
77
0.57574
0.61692
1.17584
1.38260
0.38260
78
0.58322
0.61194
1.17640
1.38391
0.38391
79
0.59069
0.60697
1.17694
1.38518
0.38518
80
0.59817
0.60199
1.17745
1.38640
0.38640
81
0.60565
0.59701
1.17795
1.38757
0.38757

105
(m)
(m)
(m)
82
0.61313
0.59204
1.17843
1.38869
0.38869
83
0.62060
0.58706
1.17889
1.38977
0.38977
84
0.62808
0.58209
1.17933
1.39081
0.39081
85
0.63556
0.57711
1.17975
1.39181
0.39181
86
0.64303
0.57214
1.18015
1.39276
0.39276
87
0.65051
0.56716
1.18054
1.39368
0.39368
88
0.65799
0.56219
1.18091
1.39456
0.39456
89
0.66547
0.55721
1.18127
1.39540
0.39540
90
0.67294
0.55224
1.18161
1.39621
0.39621
91
0.68042
0.54726
1.18194
1.39698
0.39698
92
0.68790
0.54229
1.18225
1.39772
0.39772
93
0.69537
0.53731
1.18255
1.39843
0.39843
94
0.70285
0.53234
1.18284
1.39911
0.39911
95
0.71033
0.52736
1.18311
1.39976
0.39976
96
0.71781
0.52239
1.18338
1.40038
0.40038
97
0.72528
0.51741
1.18363
1.40097
0.40097
98
0.73276
0.51244
1.18386
1.40153
0.40153
99
0.74024
0.50746
1.18409
1.40207
0.40207
100
0.74771
0.50249
1.18431
1.40259
0.40259
101
0.75519
0.49751
1.18451
1.40308
0.40308
102
0.76267
0.49254
1.18471
1.40354
0.40354
103
0.77015
0.48756
1.18490
1.40399
0.40399
104
0.77762
0.48259
1.18508
1.40441
0.40441
105
0.78510
0.47761
1.18525
1.40481
0.40481
106
0.79258
0.47264
1.18541
1.40520
0.40520

106
(m)
(m)
(m)
107
0.80005
0.46766
1.18556
1.40556
0.40556
108
0.80753
0.46269
1.18571
1.40591
0.40591
109
0.81501
0.45771
1.18585
1.40623
0.40623
110
0.82249
0.45274
1.18598
1.40655
0.40655
111
0.82996
0.44776
1.18610
1.40684
0.40684
112
0.83744
0.44279
1.18622
1.40712
0.40712
113
0.84492
0.43781
1.18633
1.40738
0.40738
114
0.85239
0.43284
1.18644
1.40764
0.40764
115
0.85987
0.42786
1.18654
1.40787
0.40787
116
0.86735
0.42289
1.18663
1.40810
0.40810
117
0.87483
0.41791
1.18672
1.40831
0.40831
118
0.88230
0.41294
1.18681
1.40851
0.40851
119
0.88978
0.40796
1.18689
1.40870
0.40870
120
0.89726
0.40299
1.18696
1.40887
0.40887
121
0.90473
0.39801
1.18703
1.40904
0.40904
122
0.91221
0.39303
1.18710
1.40920
0.40920
123
0.91969
0.38806
1.18716
1.40935
0.40935
124
0.92717
0.38308
1.18722
1.40949
0.40949
125
0.93464
0.37811
1.18727
1.40962
0.40962
126
0.94212
0.37313
1.18732
1.40974
0.40974
127
0.94960
0.36816
1.18737
1.40985
0.40985
128
0.95707
0.36318
1.18742
1.40996
0.40996
129
0.96455
0.35821
1.18746
1.41006
0.41006
130
0.97203
0.35323
1.18750
1.41016
0.41016
131
0.97951
0.34826
1.18754
1.41024
0.41024
132
0.98698
0.34328
1.18757
1.41033
0.41033

107
(m)
(m)
(m)
133
0.99446
0.33831
1.18760
1.41040
0.41040
134
1.00194
0.33333
1.18763
1.41047
0.41047
135
1.00941
0.32836
1.18766
1.41054
0.41054
136
1.01689
0.32338
1.18769
1.41060
0.41060
137
1.02437
0.31841
1.18771
1.41066
0.41066
138
1.03185
0.31343
1.18773
1.41071
0.41071
139
1.03932
0.30846
1.18775
1.41076
0.41076
140
1.04680
0.30348
1.18777
1.41081
0.41081
141
1.05428
0.29851
1.18779
1.41085
0.41085
142
1.06175
0.29353
1.18781
1.41089
0.41089
143
1.06923
0.28856
1.18782
1.41092
0.41092
144
1.07671
0.28358
1.18784
1.41095
0.41095
145
1.08419
0.27861
1.18785
1.41098
0.41098
146
1.09166
0.27363
1.18786
1.41101
0.41101
147
1.09914
0.26866
1.18787
1.41104
0.41104
148
1.10662
0.26368
1.18788
1.41106
0.41106
149
1.11409
0.25871
1.18789
1.41108
0.41108
150
1.12157
0.25373
1.18790
1.41110
0.41110
151
1.12905
0.24876
1.18790
1.41112
0.41112
152
1.13653
0.24378
1.18791
1.41113
0.41113
153
1.14400
0.23881
1.18792
1.41115
0.41115
154
1.15148
0.23383
1.18792
1.41116
0.41116
155
1.15896
0.22886
1.18793
1.41117
0.41117
156
1.16643
0.22388
1.18793
1.41118
0.41118
157
1.17391
0.21891
1.18794
1.41119
0.41119
158
1.18139
0.21393
1.18794
1.41120
0.41120

108
(m)
(m)
(m)
159
1.18887
0.20896
1.18794
1.41121
0.41121
160
1.19634
0.20398
1.18794
1.41121
0.41121
161
1.20382
0.19900
1.18795
1.41122
0.41122
162
1.21130
0.19403
1.18795
1.41122
0.41122
163
1.21877
0.18905
1.18795
1.41123
0.41123
164
1.22625
0.18408
1.18795
1.41123
0.41123
165
1.23373
0.17910
1.18795
1.41124
0.41124
166
1.24121
0.17413
1.18796
1.41124
0.41124
167
1.24868
0.16915
1.18796
1.41124
0.41124
168
1.25616
0.16418
1.18796
1.41124
0.41124
169
1.26364
0.15920
1.18796
1.41125
0.41125
170
1.27111
0.15423
1.18796
1.41125
0.41125
171
1.27859
0.14925
1.18796
1.41125
0.41125
172
1.28607
0.14428
1.18796
1.41125
0.41125
173
1.29354
0.13930
1.18796
1.41125
0.41125
174
1.30102
0.13433
1.18796
1.41125
0.41125
175
1.30850
0.12935
1.18796
1.41125
0.41125
176
1.31598
0.12438
1.18796
1.41125
0.41125
177
1.32345
0.11940
1.18796
1.41125
0.41125
178
1.33093
0.11443
1.18796
1.41126
0.41126
179
1.33841
0.10945
1.18796
1.41126
0.41126
180
1.34588
0.10448
1.18796
1.41126
0.41126
181
1.35336
0.09950
1.18796
1.41126
0.41126
182
1.36084
0.09453
1.18796
1.41126
0.41126
183
1.36832
0.08955
1.18796
1.41126
0.41126
184
1.37579
0.08458
1.18796
1.41126
0.41126

109
(m)
(m)
(m)
185
1.38327
0.07960
1.18796
1.41126
0.41126
186
1.39075
0.07463
1.18796
1.41126
0.41126
187
1.39822
0.06965
1.18796
1.41126
0.41126
188
1.40570
0.06468
1.18796
1.41126
0.41126
189
1.41318
0.05970
1.18796
1.41126
0.41126
190
1.42066
0.05473
1.18796
1.41126
0.41126
191
1.42813
0.04975
1.18796
1.41126
0.41126
192
1.43561
0.04478
1.18796
1.41126
0.41126
193
1.44309
0.03980
1.18796
1.41126
0.41126
194
1.45056
0.03483
1.18796
1.41126
0.41126
195
1.45804
0.02985
1.18796
1.41126
0.41126
196
1.46552
0.02488
1.18796
1.41126
0.41126
197
1.47300
0.01990
1.18796
1.41126
0.41126
198
1.48047
0.01493
1.18796
1.41126
0.41126
199
1.48795
0.00995
1.18796
1.41126
0.41126
200
1.49543
0.00498
1.18796
1.41126
0.41126
201
1.50290
0.00000
1.18796
1.41126
0.41126
t-H
tc = 1.5 s
H (m)
1.60
1.40
1.20
1.6
1.5
1.4
1.3
1.2
1.1
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.80
0.0
1.00
tc (s)
Figura 3.12c.Variacin de H (tc=1.5 s)

110
Se observ que para un tiempo de cierre de 1.5 s, la carga mxima fue de Hmx = 1.411 m y la
sobrepresin mxima
Con la medicin realizada, sucedi lo siguiente
Tabla 3.25f. Porcentaje de error (tc=1.5 s)
Hmx
Allievi
1.411
Medido
1.5
% error
6%
0.411
0.89
54%

111
3.10 Solucin aplicando el mtodo de las caractersticas

A diferencia de la solucin con las ecuaciones de Allievi, en este mtodo se consideran otras
variables como las prdidas por friccin, as como el uso de condiciones de frontera variables. Una
de ellas es la ley de cierre no uniforme determinada anteriormente, no obstante que en la solucin
de Allievi, se utiliz una ley de cierre lineal.
Con el fin de llevar cierto orden en el clculo, el primer paso fue obtener las condiciones iniciales a
flujo permanente. Posteriormente se aplic la condicin de estabilidad y convergencia del mtodo y
de esta forma se defini la malla caracterstica que se emple como apoyo para plantear el esquema
de diferencias finitas.
Se requiere conocer el intervalo de tiempo y los puntos de inters, para lo cual aplicamos la
condicin de estabilidad de Courant, ecuacin (2.71), de la cual podemos decir lo siguiente
Donde x se define de la siguiente forma.
donde
L longitud de la tubera, en m
N nmero mnimo entero de tramos en que debe dividirse el conducto
Si N = 5 tramos y L = 4.3 m, se tiene
As, t resulta:

112
De esta manera
Con esto se cumple con la condicin de estabilidad de Courant, sin embargo cabe mencionar que se
puede elegir un nmero mayor de tramos con el fin de obtener ms precisin en el clculo, pero en
este caso se considera adecuado el nmero de tramos elegido ya que la finalidad de este trabajo es
su aplicacin en la implementacin como prctica de laboratorio, y una malla ms grande implicara
mucho ms iteraciones creando dependencia con la capacidad de las computadoras personales de
los alumnos.
Una vez conocidos t y x, se define la malla caracterstica, la cual se encuentra representada en la
figura 3.13. Esta malla est compuesta de cinco secciones, de las cuales, la seccin 0 representa la
condicin de frontera aguas arriba (tanque con nivel constante), la seccin 5 representa la condicin
de frontera aguas abajo (vlvula de esfera) y las secciones 1, 2, 3 y 4, representan los puntos
interiores de la malla.

113
t=t 0 +4t
0.000747
t=t 0 +3 t
0.000747
t=t 0 +2 t
0.000747
t=t 0 + t
0.000747
t=t 0
0.8600
Seccin:
0.8600
1
0.8600
2
0.8600
3
0.8600
4
X=0
5
X=L
Puntos interiores ( t > t 0 )

Condiciones para el flujo permanente inicial ( t = t 0)
Figura. 3.13. Malla caracterstica del caso de estudio
Planteamiento de las ecuaciones caractersticas

Una vez definida la malla caracterstica, se realiz el clculo de las constantes necesarias para
realizar el planteamiento de las ecuaciones caractersticas.
De la ecuacin 2.57, obtenemos el valor de la constante Ca. Sustituyendo valores, result
As mismo, existe otra constante que est en funcin del factor de friccin obtenido, la cual
definimos con la siguiente ecuacin
Sustituyendo los valores se tiene

114
(
(
)(
)(
)
)
Con estos valores, se hizo el planteamiento de las constantes Cp y Cn, las cuales pueden ser
conocidas para cada instante de tiempo y que forman parte de las ecuaciones caractersticas. Estas
constantes quedaron definidas con las ecuaciones 2.55 y 2.56, y ahora, sustituyendo los valores se
tiene
|
(3.21)
(3.22)
Ahora, la ecuacin caracterstica positiva (ecuacin 2.53), queda de la siguiente forma

(3.23)
y la ecuacin caracterstica negativa (ecuacin 2.54), queda como
(3.24)
Los puntos interiores de la malla (figura 3.8) se definen a partir de la ecuacin 2.58
(
(3.25)
y de la ecuacin (3.23)
Sustituyendo el valor de la constante Ca en la ecuacin anterior, resulta

(
(3.26)
Las ecuaciones 3.25 y 3.26, definen el gasto Q y la carga de presin h en los puntos interiores de la
malla.
Condicin de frontera aguas arriba (Seccin 0)
En este caso, se tom en cuenta tanto la prdida local a la entrada de la tubera, como la carga de
velocidad, a fin de que el fenmeno transitorio sea modelado lo ms apegado a la realidad.
115
La condicin de frontera aguas arriba est representada por el depsito de carga constante y las
ecuaciones que determinan el gasto Q y la carga h en esta seccin son las ecuaciones 2.62 y 2.63.
Para el dispositivo experimental utilizado, se determin el valor del coeficiente de prdidas a la
entrada de la tubera con un valor de k = 0.5 (ver prdidas locales), lo que lleva a calcular el valor
de la constante k1 (ecuacin 2.64). Sustituyendo el valor de k y Ca en la ecuacin (2.64), dio como
resultado
Sustituyendo las constantes en la ecuacin (2.62), para obtener la carga h en la condicin de

frontera aguas arriba
(
)(
(3.27)
y sustituyendo constantes en la ecuacin (2.63), result
)[
(
(
)
)( )]
(3.28)
Nota: en las ecuaciones (2.62) y (2.63), Hdep est medido a partir del P.H.C. hasta el nivel de la
superficie libre de agua en el depsito. En este caso, la carga de posicin no est tomndose en
cuenta por encontrarse totalmente horizontal nuestro dispositivo experimental, por lo cual
Hdep = carga constante en el depsito, la cual en las mediciones se fij en 1.0 m.
Condicin de frontera aguas abajo (Seccin 5)
La variable Qp en la frontera aguas abajo, qued definida con la ecuacin 2.69, en la cual Cv est
dado por la siguiente relacin
(
) (
La ley de cierre no uniforme o curva -t, qued definida en la tabla 3.20 donde se observa el
comportamiento de la carga piezomtrica a lo largo de la tubera con respecto al tiempo de cierre
(tc) de la vlvula de esfera.

116
Se sabe que en flujo permanente Q0 = 0.00376 m3/s en cualquier seccin de la tubera y la carga en
la seccin aguas arriba de la vlvula H0v = 0.6 m, obtenido del clculo de las prdidas totales, por lo
que, la expresin para determinar el coeficiente Cv, queda de la siguiente forma
[ (
)] [(
)(
)]
Esquema de diferencias finitas.

Una vez definidas las condiciones iniciales, las condiciones de frontera aguas arriba de la tubera,
las condiciones de frontera aguas abajo de la tubera y los puntos intermedios de la malla
caracterstica aplicada al dispositivo experimental, se realiz el esquema explcito de diferencias
finitas. En esta seccin se presenta el clculo para el primer instante de tiempo, es decir para
y posteriormente, en el anexo A, se presentan las tablas donde se muestra el clculo
para cada instante de tiempo, as como para cada tiempo de cierre propuesto.
Se sabe que para
se presentan las condiciones iniciales a flujo permanente. Estos datos se
pueden ver en la tabla 3.26.

Seccin
Para el instante
hpz
(m)
(m3/s)
(m)
0.00
0.00376
0.8923
0.86
0.00376
0.8583
1.72
0.00376
0.8242
2.58
0.00376
0.7902
3.44
0.00376
0.7562
4.30
0.00376
0.6000
, es decir para
, se presenta la forma en que deben ser
calculados el gasto Q y la variacin de la carga h con respecto a la posicin en la tubera. Los

subndices que aparecern en las ecuaciones, por ejemplo Q1,0, representan la seccin (i) y el
instante de tiempo (t) respectivamente que se deben utilizar para realizar el clculo.

117
Condicin de frontera aguas arriba

Seccin 0: X = 0, i=0
|
(
)[
(
(
)|
)( )]
)[
)(
m.c.a.
Puntos interiores
Seccin 1: X=0.86m, i=1
|
(
|
(
)|
)|

118
Seccin 2: X=1.72 m, i=2

|
(
|
(
|
(
)|
)|
)|
Seccin 3: X=2.58 m, i=3

|
(

119
|
(
)|
)|
(
(
Seccin 4: X=3.44 m, i=4

|
(
|
(
|
(
)|

120
Condicin de frontera aguas abajo

Seccin 5: X=4.3 m, i=5
|
(
(
(
)|
Para un tiempo de cierre tc = 1.0 segundos (ecuacin 3.15):

(
(
(
)(
)
)
)(
))

121
De esta forma, se determinaron las variables Qp y hp a lo largo de la tubera para cada instante de
tiempo, as como para cada tiempo de cierre propuesto.
3.10.1 Resultados
A continuacin se presentan las grficas donde se puede observar la variacin de la carga de presin
(h) con respecto a la posicin en la tubera (X) y el instante de tiempo que se presenta (t). As
mismo, las tablas con los resultados de los clculos correspondientes, se presentan en el anexo A.

122
Tiempo de cierre (tc = 1.0 s)

X-h
tc = 1.0 s
hpz (m.c.a.)
2.0
t=0.0 s
1.8
t=0.2 s
1.6
t=0.4 s
1.4
t=0.6 s
t=0.8 s
1.2
t=1.0 s
1.0
0.8
0.6
0.4
0.2
0.0
0.00
0.86
1.72
2.58
3.44
4.30
X (m)
Figura 3.14. Carga de presin Posicin (tc=1.0 s)
Carga Piezomtrica
Seccin 1 (X=0.86 m)
hpz ( m.c.a. )
1.30
1.20
1.10
1.00
0.90
0.80
0.70
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.60
Tiempo ( s )
Figura 3.15. Seccin 1, Carga de presin Tiempo (tc=1.0 s)
Carga Piezomtrica
Seccin 2 (X=1.72 m)
hpz ( m.c.a. )
1.40
1.30
1.20
1.10
1.00
0.90
0.80
0.70
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.60
Tiempo ( s )

123
Carga Piezomtrica
Seccin 3 (X=2.58 m)
hpz ( m.c.a. )
1.60
1.40
1.20
1.00
0.80
0.60
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )
Carga Piezomtrica
Seccin 4 (X=3.44 m)
hpz ( m.c.a. )
1.80
1.60
1.40
1.20
1.00
0.80
0.60
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )
Carga Piezomtrica
Seccin 5 (X=4.3 m)
hpz ( m.c.a. )
2.00
1.80
1.60
1.40
1.20
1.00
0.80
0.60
0.40
0.20
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.00
Tiempo ( s )

124

X-h
tc = 1.2 s
hpz (m.c.a.)
2.0
t=0 s
1.8
t=0.2 s
1.6
t=0.4 s
t=0.6 s
1.4
t=0.8 s
1.2
t=1.0 s
1.0
t=1.2 s
0.8
0.6
0.4
0.2
0.0
0.00
0.86
1.72
2.58
3.44
4.30
X (m)
Carga Piezomtrica
Seccin 1 (X=0.86 m)
hpz ( m.c.a. )
1.20
1.15
1.10
1.05
1.00
0.95
0.90
0.85
0.80
0.75
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.70
Tiempo ( s )

125
Carga Piezomtrica
Seccin 2 (X=1.72 m)
hpz ( m.c.a. )
1.40
1.30
1.20
1.10
1.00
0.90
0.80
0.70
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.60
Tiempo ( s )
Carga Piezomtrica
Seccin 3 (X=2.58 m)
hpz ( m.c.a. )
1.60
1.40
1.20
1.00
0.80
0.60
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )
Carga Piezomtrica
Seccin 4 (X=3.44 m)
hpz ( m.c.a. )
1.80
1.60
1.40
1.20
1.00
0.80
0.60
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )

126
Carga Piezomtrica
Seccin 5 (X=4.3 m)
hpz ( m.c.a. )
1.80
1.60
1.40
1.20
1.00
0.80
0.60
0.40
0.20
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.00
Tiempo ( s )

127

X-h
tc=1.5 s
hpz (m.c.a.)
1.8
t=0 s
t=0.2 s
1.6
t=0.4 s
1.4
t=0.6 s
1.2
t=0.8 s
t=1.0 s
1.0
t=1.2 s
0.8
t=1.4 s
t=1.5 s
0.6
0.4
0.2
0.0
0.00
0.86
1.72
2.58
3.44
4.30
X (m)
Carga Piezomtrica
Seccin 1 (X=0.86 m)
hpz ( m.c.a. )
1.20
1.10
1.00
0.90
0.80
0.70
1.40
1.50
1.40
1.50
1.30
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.60
Tiempo ( s )
Carga Piezomtrica
Seccin 2 (X=1.72 m)
hpz ( m.c.a. )
1.30
1.20
1.10
1.00
0.90
0.80
0.70
1.30
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.60
Tiempo ( s )

128
Carga Piezomtrica
Seccin 3 (X=2.58 m)
hpz ( m.c.a. )
1.40
1.20
1.00
0.80
0.60
1.50
1.40
1.30
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )
Carga Piezomtrica
Seccin 4 (X=3.44 m)
hpz ( m.c.a. )
1.60
1.40
1.20
1.00
0.80
0.60
1.40
1.50
1.40
1.50
1.30
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )
Carga Piezomtrica
Seccin 5(X=4.3 m)
hpz ( m.c.a. )
1.60
1.40
1.20
1.00
0.80
0.60
1.30
1.20
1.10
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )

129
De los resultados obtenidos, se observa la variacin de la carga de presin respecto a la posicin en

la tubera y respecto a cada instante de tiempo. Para el tiempo de cierre propuesto de tc=1.0
segundos, la carga mxima calculada en la seccin 5, tuvo una variacin entre 1.8 y 1.9 metros.
Respecto al tiempo de cierre tc=1.2 segundos, la carga mxima calculada en la seccin 5 fue de
1.73 m y por ltimo al realizar el clculo para un tiempo de maniobra de cierre tc=1.5 segundos, la
carga mxima calculada en la seccin 5 de la tubera, fue de 1.56 m.
Recordando que la carga de presin mxima que se midi en el dispositivo experimental en la
seccin 5 fue de 1.5 metros, se puede decir que el clculo realizado con un tiempo de cierre
propuesto de tc=1.5 segundos, es el que ms se aproxima a lo medido. Sin embargo, aunque sigue
existiendo una diferencia del 5 por ciento, debido a la posicin que existe entre el piezmetro
instalado y la seccin 5 propuesta tericamente, se considera una buena aproximacin ya que
tambin intervienen factores como la exactitud al determinar la ley de cierre o al determinar las
prdidas por friccin y locales. Tambin existe otro factor importante, como los dispositivos
utilizados para medir la carga de presin, conocidos como piezmetros, los cuales se considera que
no tienen la suficiente resolucin para tomar una medida ms exacta.
Es por ello que se agregarn a este trabajo algunas propuestas para poder obtener mejores
mediciones, as como algunas modificaciones al dispositivo experimental para obtener resultados
ms representativos del fenmeno transitorio.

130
CAPTULO 4. PROPUESTAS DE DISPOSITIVOS EXPERIMENTALES
PROPUESTAS DE DISPOSITIVOS EXPERIMENTALES
En el laboratorio de hidrulica de la Facultad de Ingeniera UNAM existen dos dispositivos

experimentales para reproducir el fenmeno de golpe de ariete y realizar las mediciones asociadas.
Uno de ellos es el que se ha utilizado en el caso de estudio del captulo 3 y consta de tanque
tubera vlvula de esfera. Existe otro dispositivo que consta de instalacin hidrulica tubera en
espiral vlvula, pero no se encuentra en operacin.
Estos dispositivos podran funcionar mejor con equipo de medicin y control adecuados, por ello se
proponen las siguientes mejoras a los equipos.
4.1
Mejoras al dispositivo experimental tanque tubera vlvula esfrica
Las mejoras que se proponen al dispositivo experimental actual tienen por objetivo el poder
observar la variacin real de la onda de presin al presentarse el golpe de ariete y compararla con la
presin simulada mediante el mtodo de las caractersticas.
El dispositivo experimental actual permite observar en los piezmetros el aumento de la presin
cuando se presenta el golpe de ariete. La medicin que se hace en los piezmetros al momento de
realizar el cierre brusco en la vlvula suele ser poco precisa ya que depende de la vista del
observador, la cual puede incurrir en errores de medicin. El dato que se obtiene de la medicin es
puntual y no permite observar la variacin de la onda de presin.
Adems de los piezmetros se propone utilizar transductores de presin que tienen la capacidad de
calibrarse para medir seales elctricas y convertirlas en valores de presin.
4.1.1.
Implementacin de los instrumentos de medicin y control adecuados
Como se observa en la figura 4.1 se puede aprovechar la infraestructura del dispositivo

experimental actual. Las mejoras consisten en instalar dos sensores de presin con un intervalo de
medicin de 0 a 50 psi (0 a 3.5 kg/cm2) alimentados por una fuente de voltaje. Uno de los sensores
sera instalado en el extremo inicial de la tubera (sensor 1, figura 4.1) y el segundo sensor, sera
instalado en el extremo final de la tubera, antes de la vlvula de esfera (sensor 2, figura 4.1). Estos
sensores podrn monitorear la variacin en la onda de presin durante las maniobras de cierre o
apertura de la vlvula de esfera, as como almacenar estos datos en una tarjeta de adquisicin de
datos. Se tiene la ventaja de poder conectar una computadora a la tarjeta de adquisicin de datos y
desplegar mediante un software en tiempo real la variacin de la onda de presin.
131
COMPUTADORA
12
TARJETA DE
ADQUISICIN
DE DATOS
5
10
11
Ho = 1.0 m
Piezmetros
Tubera de cobre, L=4.3 m, D=2.5"
Te con bridas, D=2.5"
Vlvula de esfera, D=2.5"
Vertedor triangular B=1.7, 0=45
10
Sensor de presin 1
11
Sensor de presin 2
12
Tarjeta de adquisicin de datos
Nsa
Nc
Figura 4.1. Mejoras al dispositivo experimental tanque tubera vlvula de esfera
En lo que respecta al sistema electrnico del equipo (figura 4.2) la seal de voltaje que transmite
cada sensor es enviada a la tarjeta de adquisicin de datos y mediante un software, estas seales se
pueden desplegar y almacenar en una computadora personal. Tambin se pueden programar las
tareas necesarias para adquirir las seales de voltaje enviadas por los sensores de presin.

132
TARJETA
DE
ADQUISICIN
DE
DATOS
(DAQ)
SENSOR DE PRESIN 1
SENSOR DE PRESIN 2
Figura 4.2. Esquema del sistema electrnico sensor DAQ
El costo aproximado en Moneda Nacional que implican estas mejoras, se resume en la tabla 4.1.
Tabla 4.1. Presupuesto Opcin 1
CONCEPTO
Suministro de sensor de presin relativa
de 50 PSI, con salida 0-10 VDC,
alimentacin 14-30VDC, material acero
inoxidable, con rosca NPT 1/4-18
macho.
Suministro de Fuente de poder
24V/1.04A, dimensiones 104mm X
35mm X 45mm, montaje en riel DIN,
voltaje de entrada 100-240 VAC.
Suministro de Tarjeta USB con 8
entradas y 2 salidas anlogas de voltaje,
16 E/S digitales, frecuencia de
muestreo de 48 KS/s (Tarjeta de
adquisicin de datos (DAQ)).
Software TracerDAQ profesional.
UNIDAD
CANTIDAD
P.U.
IMPORTE
Pza.
$2,494.08
$4,988.16
Pza.
$622.08
$622.08
Pza.
$3,627.00
$3,627.00
Pza.
$2,786.00
$2,786.00
SUBTOTAL $12,023.24
$1,923.72
IVA 16%
$13,946.96
TOTAL

133
4.2. Propuesta de nuevo dispositivo

Otro de los dispositivos experimentales utilizados para reproducir el golpe de ariete es el que consta
de una instalacin hidrulica tubera en espiral vlvula. Con las mejoras que a continuacin se
presentan se pretende rehabilitar este dispositivo para su correcto funcionamiento ya que
actualmente no se encuentra en operacin.
4.2.1. Instalacin hidrulica red-tubera en espiral-vlvula.
El dispositivo mostrado en la figura 4.3 est formado por una tubera de cobre en espiral de
0.0127 m (1/2) de dimetro y 61 metros de longitud. Cuenta con la conexin a un manmetro de
cartula y un osciloscopio donde anteriormente se observaba la variacin de la presin en el interior
de la tubera. Tiene tambin una vlvula solenoide la cual puede cerrarse bruscamente para producir
una sobrepresin en la tubera que es abastecida por medio de su conexin a la red hidrulica.
MANMETRO
OSCILOSCOPIO
ESPIRAL DE TUBERA DE
COBRE
D = 1.27 cm
L = 61 m
e = 0.119 cm
Figura 4.3. Instalacin hidrulica tubera en espiral vlvula
Los accesorios como el osciloscopio y la vlvula solenoide no se encuentran en funcionamiento por

lo cual se propone remplazarlos por otros accesorios que permitan reproducir y observar el
fenmeno del golpe de ariete.
4.2.2. Instrumentos de medicin y control
Existen instrumentos de medicin y control de funcionamiento electrnico o electromecnico que
pueden facilitar la operacin de este tipo de sistemas, as como permitir el almacenamiento de
registros. Tales equipos son
1. Sensores de presin
2. Fuente de voltaje
3. Relevador
134
4. Vlvula solenoide
5. Tarjeta de adquisicin de datos
6. Software
A continuacin se muestra una descripcin del funcionamiento y caractersticas de los instrumentos
de medicin y control que se proponen para el nuevo dispositivo experimental.
1. Sensores de presin
Se llama sensor al instrumento que produce una seal, usualmente elctrica, que refleja el valor de
una propiedad, mediante alguna correlacin definida.
En trminos estrictos, un sensor es un instrumento que no altera la propiedad sensada; existe,
adems, el concepto estricto de transductor: un instrumento que convierte una forma de energa en
otra (o una propiedad en otra).
Para este caso los sensores de presin a implementar tanto para el dispositivo experimental que
consta de tanque tubera vlvula de esfera y el dispositivo que consta de instalacin red tubera
en espiral vlvula, entregan una seal analgica de 4 a 20 miliamperes. Cada sensor posee un
mdulo de acondicionamiento de seal capaz de entregar una seal analgica de voltaje de 0 10
volts de corriente directa, proporcional a la presin medida.
Equipo a utilizar
Los sensores de presin OEM (Fabricante de Equipos Originales) serie PMP300 recomendados para
estas opciones estn diseados con una precisin de ingeniera para adaptarse a muchas
aplicaciones de medicin de presin OEM e industriales. Su robusto diseo es resistente a la
vibracin, golpes, grandes variaciones de temperatura, RFI (Solicitud de Informacin) y muchas
otras condiciones ambientales extremas que son tpicas en aplicaciones industriales y OEM.
Su funcionamiento y fiabilidad son respaldados por la soldadura de todo el cuerpo de acero
inoxidable de la celda de medicin que elimina la necesidad de cubiertas en materiales suaves que
pueden deteriorarse con el tiempo.
Figura 4.4. Sensor de presin

135
Caractersticas
OEM
Salida amplificada
De bajo costo
Tamao compacto
Intervalo de presin de 0 a 50 PSI (0 a 3.5 Kg/cm2)
Tabla 4.2. Especificaciones de los sensores de presin
Canales
Temperatura compensada
Excitacin
Histresis
No linealidad
No repetibilidad
Temperatura en operacin
Temperatura Shift Span
Temperatura Shift Zero
Salida
Excitacin calibrada a
Altura
Peso
Dimetro
Sensor de presin
1
Min: 0 C
Max: 80 C
14 Vdc a 30 Vdc para transductores
con salida a voltaje
10 Vdc a 30 Vdc para transductores
con salida de corriente
Min: -0.16% de R.O.
Max: 0.16% de R.O.
Min: -1% de R.O.
Max: 1% de R.O.
Min: -0.1% R.O.
Max: 1% de R.O.
Min: 0 C, Max: 80 C
Min: -0.02% de carga/F
Max: 0.02% de carga/F
Min: -0.02% de R.O./F
Max: 0.02% de R.O./F
De 0 10 Vdc 4 20 mA
24 Vdc
2.4
2.8 oz
1.06

136
2. Fuente de voltaje LBCSP-142

Una fuente de voltaje es un dispositivo que convierte la tensin alterna de la red de suministro, en
una o varias tensiones, prcticamente continuas, que alimentan los distintos circuitos del aparato
elctrico al que se conecta.
Los sensores de presin a implementar requieren una fuente de voltaje para su alimentacin y para
proteger a los sensores de presin contra sobre voltajes y as evitar que puedan ser daados o arrojar
lecturas errneas.
La fuente de voltaje que se propone para alimentar a los dos sensores de presin es la siguiente:
Equipo a utilizar
Figura 4.5. Fuente de voltaje
Caractersticas
Posee una gran variedad de rangos de voltaje/corriente para diferentes tipos de aplicaciones
Diseo compacto con montaje en riel DIN
Opcin de RoHS
De bajo costo

137
Tabla 4.3. Especificaciones de la fuente voltaje
Fuente de voltaje
Modelo
25W
ENTRADA
Voltaje
100 ~ 240 Vac, 47 ~ 63 Hz
Eficiencia (a 110V de
entrada)
Corriente pico de
entrada
24V
76%
110V
20A
220V
40A
SALIDA
Indicador de salida
LED Verde
Ajuste de voltaje
-10% a 10% (Con ajuste de VR)
Ripple (Vp-p)
24V
150mV
Tolerancia de voltaje
1.0% max. (con entrada ajustada, 0 a 100% de carga)
Linealidad
Influencia de la carga
Tiempo de inicio
Mximo de 3 segundos
Tiempo en espera
13ms
Proteccin contra corto
En el modelo estndar, se reinicia automticamente
ESPECIFICACIONES AMBIENTALES
Temperatura ambiente
Operacin:-15C ~ 50C; Almacn:-20C~85C (Sin condensacin)
Humedad
Operacin: 5 a 95% RH; Almacn: 5 a 95% RH
Influencia de la temperatura
Mximo 0.05% / C ( 0 ~ 50 C)
Influencia del voltaje
Mximo 0.5% ( de 85 a 264 VAC de entrada, 100% de carga)
Vibracin
10 a 500 Hz, 2G, 10 min/ciclo en una hora cada direccin X, Y y Z
Resistencia al golpe
150 m/s2 (como 50 g) cada 3 veces en las direcciones X, Y, y Z
SEGURIDAD ELCTRICA
Estrecho dielctrico
2.0 kV AC en 1 minuto (Entre la entrada, salida, carcasa; deteccin: 20mA)
Dielctrico
Mnimo 100M (Entre entrada, salida y carcasa) a 500Vdc
EMC
EN55022 Clase B, EN61000-3-2, EN61000-3-3, EN55024, IEC61000-4-2, -3, 4, -5, -6, -8, -11
ESPECIFICACIONES MECNICAS
Dimensiones
104 x 35 x 45
Material
ABS negro. UL94V0
Terminales
Terminales a tornillo, hasta 2 x 2.5 mm, 2 cables
Montaje
Riel DIN de 35mm (EN50022)
Peso
190 g

138
3. Relevador
El rel o relevador es un dispositivo electromecnico. Funciona como un interruptor controlado
por un circuito elctrico en el que, por medio de una bobina y un electroimn, se acciona un juego
de uno o varios contactos que permiten abrir o cerrar otros circuitos elctricos independientes.
Dado que el rel es capaz de controlar un circuito de salida de mayor potencia que el de entrada,
puede considerarse, en un amplio sentido, como un amplificador elctrico. Como tal se emplearon
en telegrafa, haciendo la funcin de repetidores que generaban una nueva seal con corriente
procedente de pilas locales a partir de la seal dbil recibida por la lnea. Se les llamaba
relevadores. De ah rel.
Equipo a utilizar
En este caso se utilizar una tarjeta modelo CIO-ERB08, con su respectivo conector C37FF y su
terminal a tornillos CIO-MINI37, la cual actuar como interfaz de control del cierre y apertura de la
vlvula solenoide. Dicha tarjeta estar conectada a la tarjeta de adquisicin de datos.
Caractersticas del equipo a emplear
8 relevadores forma C
Reset en encendido
Todos los relevadores son reemplazables
Tarjeta CIO-ERB08 con 8 relevadores reemplazables
Conector C37FF
Terminal a tornillos CIOMINI37
Figura 4.6. Tarjeta CIO-ERB08 y accesorios

139
Tabla 4.4. Especificaciones del relevador
Relevadores
Capacidad de los contactos
Tiempo de funcionamiento
Tiempo de liberacin
Consumo de energa
Alimentacin requerida
Relevador
8 SPDT, forma C
10A a 28 VDC, 10A a 120 VAC
20 milisegundos
10 milisegundos
mximo de 0.95A a +5V
+ 5V, alimentacin externa o desde la PC con el
cable CPCPOWER-10 (incluido)

140
4. Vlvula solenoide
En muchas aplicaciones es necesario controlar el paso de algn tipo de flujo, desde corriente
elctrica hasta gases o lquidos. Esta tarea es realizada por vlvulas. En particular, las accionadas
por solenoides permiten su implementacin en lugares de difcil acceso y facilitan la automatizacin
del proceso al ser accionadas elctricamente.
Este tipo de vlvula se opera variando la corriente que circula a travs de un solenoide (conductor
ubicado alrededor de un mbolo, en forma de bobina). Esta corriente, al circular por el solenoide,
genera un campo magntico que atrae un mbolo mvil. Por lo general estas vlvulas operan de
forma completamente abierta o completamente cerrada, aunque existen aplicaciones en las que se
controla el flujo en forma lineal.
Al finalizar el efecto del campo magntico, el mbolo vuelve a su posicin por efecto de la
gravedad, un resorte o por presin del fluido a controlar.
Figura 4.7. Vlvula Solenoide
Para el dispositivo red tubera en espiral vlvula solenoide, se implementar esta vlvula para
efectuar el cierre brusco del flujo en la tubera y as provocar el golpe de ariete. Como ya se
coment el cierre y la apertura de esta vlvula se controlar mediante el relevador elegido.

141
Tabla 4.5. Especificaciones de la vlvula solenoide
Modelo
Vlvula solenoide
2PO160
Tipo de vlvula
2 vas normalmente abierto (NO)
Accin
Diafragma ascensor directo, Unidireccional, diferencial cero
Cv (orificio, en mm)
4.8 (16)
Presin de trabajo
0 a 90 PSI
Prueba de presin
200 PSI
Temperatura de funcionamiento
Media: -5 a 50 C; Ambiente: -5 a 45 C
Tamao de puerto (NPT)
1/2
Materiales del cuerpo
Nylon de ingeniera
Materiales de sellado
NBR (Buna N), Opciones: VITON, EPDM
Bobina de servicio
Clase H, IP65, 100% ED
Voltaje
Opciones de voltaje: 12, 24 VDC; 24, 110/120, 220/240 VAC,

50/60 Hz
Tolerancia de tensin
10% del voltaje especificado
Configuracin de los contactos
14 a 30 W
Conexiones elctricas
DIN 43650, forma A
Servicio
Lquido

142
5. Tarjeta de adquisicin de datos modelo USB 1408FS

La adquisicin de datos o adquisicin de seales, consiste en la toma de muestras del mundo real
(sistema analgico) para generar datos que puedan ser procesados por una computadora. Consiste,
en tomar un conjunto de seales fsicas, convertirlas en tensiones elctricas y digitalizarlas de
manera que se puedan procesar en una computadora. Se requiere una etapa de acondicionamiento,
que adecua la seal a niveles compatibles con el elemento que hace la transformacin a seal
digital. El elemento que hace dicha transformacin es el mdulo de digitalizacin o tarjeta de
adquisicin de datos (DAQ).
Las seales procedentes de cada sensor de presin sern enviadas a la tarjeta de adquisicin de
datos cuyas caractersticas son las siguientes:
Equipo a utilizar
Figura 4.8. Tarjeta de Adquisicin de Datos
El modelo USB 1408FS es un dispositivo basado en la tecnologa USB que fija un nuevo estndar
para la rentabilidad, basado en la PC para entradas y salidas anlogas y digitales. Proporcionando 8
entradas single ended, resolucin de 13 bits o 4 entradas diferenciales con resolucin de 14 bits,
las unidades ofrecen tasas de muestreo de hasta 48 kilomuestras/segundo con 8 rangos de entrada
seleccionables por software. La unidad tambin proporciona 2 canales de salida anloga de 12 bits,
un contador de 32 bits y 16 bits de entradas y salidas digitales.
El USB 1408FS es confiable y bastante resistente para cualquier aplicacin DAQ. Si no fuera as,
no vendra con una garanta de por vida y una garanta por daos causados.
Caractersticas del equipo a emplear
8 entradas anlogas con terminal comn single ended y 4 diferenciales

Resolucin 14 bits (Diferencial) / 13 bits (SE)
Tasa de muestreo de hasta 48 kilomuestras/segundo
2 salidas anlogas de 12 bits
16 bits de entrada y salida digitales

143
6. Software a emplear
La tarjeta de adquisicin de datos es conectada a una computadora personal. El software a emplear
Tracer DAQ profesional sirve para desplegar los valores obtenidos por la tarjeta, as como graficar y
realizar registros en la computadora para poder analizarlos posteriormente. Otra de las funciones del
software es programar las funciones que corresponden a la adquisicin de las dos seales de voltaje
[0-10 VDC] entregadas por los dos sensores de presin y programar el envo de las seales de
control [0-5 VDC] de la tarjeta de adquisicin de datos al relevador de control de la vlvula
solenoide.
TracerDAQ profesional incluye las siguientes aplicaciones:
Graficador. Aplicacin de 48 canales para registrar y graficar datos de entradas analgicas,
digitales, contadores y temperaturas de una o ms tarjetas DAQ de Measurement Computing.
Caractersticas de muestreo basadas en la velocidad de la tarjeta utilizada (hasta 1 milln de
muestras por canal), triggering usando un disipador de software interno, un disipador de software
externo (o un disipador por intervalos de tiempo), la habilidad de configurar alarmas condicionales,
agregar comentarios (anotaciones) al display de plot y volver a mostrar una adquisicin.
Osciloscopio. Aplicacin de osciloscopio con 4 canales que muestra los valores adquiridos de las
entradas anlogas. Caractersticas de muestreo basadas en la velocidad de la tarjeta utilizada y
triggering usando un disipador de software interno o un disipador de hardware externo, la habilidad
para agregar una referencia o un canal matemtico, muestra valores de medidas especficas para
cada canal y datos de canal a escala.
Generador de funciones. Genera mltiples formas de onda de 16 canales (seno, cuadrada,
triangular, constante, sierra y personalizable). La aplicacin de Generador de Funciones, genera
formas de onda para salidas anlogas a la tasa mxima del dispositivo de Measurement Computing.
El generador de funcin muestra la forma de onda, muestras por ciclo, tasa de escaneo interno y
frecuencia base de la seal de salida del canal anlogo seleccionado. El usuario puede efectuar
barrido logaritmico y lineal, cambiar la fase, ciclo de trabajo, mltiplo de pulso y amplitud de la
forma de onda.
Generador de pulsos. Aplicacin Generador de pulsos de 20 canales que muestra la forma de
onda, muestras por ciclo, tasa de escaneo interno y frecuencia base de salida del canal
contador/timer seleccionado. El usuario puede cambiar la frecuencia, el ciclo de trabajo y el estado
inicial de la forma de onda.
Incluye el CD MCC DAQ con:
Software Graficador y Registrador GRATUITO TracerDAQTM 2.0

Libreras y ejemplos para C, C++, VB, Delphi, etc.
Vis y ejemplos para LabVIEWTM
Software configurador InstaCal

144
4.2.3. Implementacin de los instrumentos de medicin y control en la instalacin hidrulica

propuesta.
Una vez elegidos los instrumentos de medicin y control, se debern implementar en la instalacin
hidrulica existente que se describi en el apartado 4.2.1, de la cual se aprovechar la tubera en
espiral de cobre, de 1/2 de dimetro y 61 m de longitud.
A esta tubera se le instalarn 2 sensores de presin (tabla 4.2), los cuales se colocarn en la parte
media de la tubera (sensor 1) y en el extremo final de la misma (sensor 2). Se ha elegido que los
sensores de presin tengan un intervalo de medicin de 0 a 50 PSI (0 A 3.5 Kg/cm2), ya que es el
rango en el que se espera se generen las sobrepresiones mximas.
Los sensores de presin elegidos necesitan ser alimentados por una fuente de voltaje (tabla 4.3).
Esta fuente de voltaje a su vez funcionar para proteger a los sensores de presin contra sobre
voltajes y a su vez para que no se enven mediciones errneas.
Al final de la tubera, despus del sensor aguas abajo, se instalar una vlvula solenoide (tabla 4.5)
con el fin de garantizar la interrupcin rpida del flujo para generar el golpe de ariete. Dicha vlvula
solenoide ser activada por una interfaz llamada relevador (tabla 4.4).
En lo que respecta al sistema electrnico y de adquisicin de datos, se utilizar una tarjeta DAQ
(Tarjeta de adquisicin de datos), la cual ser la encargada de adquirir las seales procedentes de
cada sensor.
A su vez la tarjeta de adquisicin de datos alimentar directamente al relevador el cual acta como
interfaz del control de cierre y apertura de la vlvula solenoide.
La tarjeta de adquisicin de datos estar conectada a una computadora personal.
El software a emplear servir para desplegar los valores obtenidos por la tarjeta, para graficar y
registrar los datos de las mediciones en una computadora para su anlisis posterior.
En la figura 4.9 se presenta el arreglo de la instalacin hidrulica propuesta con los instrumentos de
control y medicin adaptados.
En la figura 4.10 se muestra el esquema del arreglo electrnico de los instrumentos de control y
medicin.

145
Descarga libre
Q
Vlvula de
globo
Vlvula
solenoide
Sensor de
presin 2
Q
Espiral de
tubera de cobre
D = 1.27 cm
L = 61 m
e = 0.119 cm
Qred
Q
Vlvula
de globo
Sensor de
presin 1
Abastecimiento
de la red
hidrulica
Figura 4.9. Arreglo de la instalacin hidrulica propuesta

146
120
VAC
TARJETA
DE
ADQUISICIN
DE
DATOS
(DAQ)
RELEVADOR
ELECTROVLVULA
SENSOR DE PRESIN 1
SENSOR DE PRESIN 2
Figura 4.10. Arreglo electrnico

147
El costo aproximado en Moneda Nacional que implican estas mejoras, se resume en la tabla 4.6.
Tabla 4.6. Presupuesto Opcin 2
CONCEPTO
Suministro de sensor de presin relativa
de 50 PSI, con salida 0-10 VDC,
alimentacin 14-30VDC, material acero
inoxidable, con rosca NPT 1/4-18
macho.
Suministro de Fuente de poder
24V/1.04A, dimensiones 104mm X
35mm X 45mm, montaje en riel DIN,
voltaje de entrada 100-240 VAC.
Suministro de Tarjeta USB con 8
entradas y 2 salidas anlogas de voltaje,
16 E/S digitales, frecuencia de
muestreo de 48 KS/s (Tarjeta de
adquisicin de datos (DAQ))
Software TracerDAQ profesional
Suministro de Mdulo de 8 relevadores
tipo C de 10 A. Salidas del relevador
con terminales a tornillo. Entrada de
control con conector de 37 pines.
Entrada de voltaje de 5V (se puede
tomar desde la fuente de la PC con el
cable incluido). Incluye etapa driver de
control
Suministro de Clave de 37 conductores,
conectores hembra a hembra, 2 ft largo.
Suministro de Terminal universal a
tornillos 4 x 4
Suministro de Vlvula solenoide de
Suministro de Vlvula de globo de
UNIDAD
CANTIDAD
P.U.
IMPORTE
Pieza.
$2,494.08
$4,988.16
Pieza.
$622.08
$622.08
Pieza.
$3,627.00
$3,627.00
Pieza.
$2,786.00
$2,786.00
Pieza.
$3,144.40
$3,144.40
Pieza
$558.18
$558.18
Pieza
$1,469.86
$1,469.86
Pza.
Pza.
1
2
$1,202.50
$300.00
$1,202.50
$600.00
SUBTOTAL $18,998.18
IVA 16%
$3,039.71
TOTAL
$22,037.89

148
CAPTULO 5. PROPUESTA DE PRCTICAS DE LABORATORIO
5. PROPUESTA DE PRCTICAS DE LABORATORIO

experimental existente
Parte del objetivo de esta tesis, fue vincular las bases tericas del mtodo de las caractersticas, con
la parte experimental del fenmeno transitorio del golpe de ariete en una prctica de laboratorio. Es
por ello que, en base a los resultados adquiridos en el caso de estudio, se presenta una propuesta de
prctica de laboratorio, para que el alumno pueda conocer y tener la capacidad de aplicar el mtodo
de las caractersticas.
En esta propuesta de prctica de laboratorio se plantean los objetivos que el alumno deber alcanzar
en la realizacin de la prctica, as como los antecedentes para comprender de mejor manera el
fenmeno transitorio. En seguida se muestra el desarrollo que se debe seguir durante la prctica,
con el fin de que el alumno tome las lecturas adecuadas, las cuales sern la base para iniciar el
anlisis del golpe de ariete.
Posteriormente se indica la forma en que se deber realizar la memoria de clculo, en la cual est
incluida la solucin por medio del mtodo de las caractersticas, el cual se pretende explicar de la
manera ms sencilla y didctica posible.
Como se vio en la solucin presentada en el captulo tres, el criterio de estabilidad y convergencia
del mtodo, no permite utilizar una malla con menos de cinco tramos, ya que los resultados podran
no ser confiables, es por esta razn, que la malla propuesta para la prctica, sea igual a la presentada
en la solucin del captulo tres, con la opcin de que el alumno realice los clculos en su
computadora personal y presente nicamente los resultados por medio de grficas.
Finalmente, se muestra cmo deben ser analizados los resultados obtenidos en la memoria de
clculo, seguido de una serie de preguntas por medio de las cuales se reflejen los aspectos ms
importantes del tema en estudio, as como del mtodo empleado en la solucin.
Se presenta a continuacin la propuesta para la realizacin de la prctica de laboratorio Golpe de
ariete aplicando el mtodo de las caractersticas.

149
Objetivo
Observar la variacin de la carga de presin a lo largo de una tubera, provocada por el cierre y
apertura de una vlvula colocada en su extremo final.
Calcular y analizar la variacin de la carga de presin generada en una tubera por el cierre de una
vlvula, aplicando el mtodo de las caractersticas.
Antecedentes
Flujo permanente en conductos a presin
Flujo transitorio en conductos a presin
Teora del golpe de ariete
Celeridad de una onda de presin
Clasificacin de maniobras de vlvula
Desarrollo
Medicin de condiciones iniciales en flujo permanente.
1. Abrir la vlvula de alimentacin del tanque y establecer una carga constante H0 en el depsito y
registrar su valor, en m. Ver figura 5.1.
2. Abrir parcialmente la vlvula de esfera ubicada en el extremo final de la tubera, hasta alcanzar
el nivel de cresta Nc en el canal de aproximacin y registrar la lectura correspondiente Nc en el
vertedor de seccin triangular.
3. Abrir totalmente la vlvula de esfera y esperar a que se establezca el flujo. Cuando esto suceda,
medir el nivel Nsa, en m, en el vertedor triangular.
4. Purgar los piezmetros y medir la carga de presin H1, H2, H3 y H4, en m. Ver figura 5.1
5. Cerrar bruscamente la vlvula de esfera y en el mismo instante, medir la carga de presin
mxima Hmx, en el piezmetro no. 3 o 4. Ver figura 5.1.
6. Finalmente, realizar una maniobra de apertura y observar el comportamiento de la columna de
agua en el piezmetro.

150
5
Hmx
H0
L
0
Piezmetros
Nsa
Nc
5 Depsito con nivel constante

6 Tubera de cobre, L=4.3 m, D=2.5"
7 Te con bridas, D=2.5"

8 Vlvula de esfera, D=2.5"
9 Vertedor triangular B=1.7, =45
Registro de datos
Tabla 5.1. Registro de datos
Tanque de
carga
constante
Vertedor
triangular
Piezmetro
H0 =
Nc =
Nsa =
H1 =
H2 =
H3 =
H4 =
H4mx =
m
m
m
m
m
m
m

151
Memoria de clculo
1. Calcular el gasto Q de la tubera, en m3/s, que pasa por el vertedor triangular, segn la ecuacin
donde:
coeficiente de descarga del vertedor, en m1/2/s
carga en el vertedor,
|, en m;
|
2
aceleracin de la gravedad, en m/s ;
ngulo del vrtice del vertedor, =45;
ancho del canal de aproximacin, B=1.70 m;
coeficiente experimental que depende de y h segn figura 7.9 de la referencia 1;
coeficiente experimental que depende de B y h segn figura 7.10 de la referencia 1.
2. Calcular la celeridad de la onda a, en m/s, de la onda de presin en la tubera, con la ecuacin:

(
donde:
a
celeridad de la onda de presin, en m/s;
densidad del agua,
E0
;
Em
mdulo de elasticidad de la tubera,

;
dimetro de la tubera,
espesor de la tubera,
;
.
3. Calcular el periodo T, en segundos, de la onda de presin, como:
, donde L es la longitud
de tubera. Ver figura 5.1

152
4. Solucin aplicando el Mtodo de las Caractersticas

Ecuaciones caractersticas:
( )
( )
en las cuales:
|
( )
( )
donde:
Considerar f=0.035, t=0.000747 segundos y x=0.86 m. Ver figura 5.2

Puntos interiores:
( )
( )
Condiciones de frontera:
Condicin de frontera aguas arriba. Tanque con nivel constante
( )
donde:
(
( )
Considerar k = 0.5 (Coeficiente de prdida a la entrada de la tubera).

(
( )

153

(
donde:
(
) (
Considerar H0v=0.6 m (Carga de presin en la vlvula en flujo permanente)

Siendo la ley de cierre no uniforme definida para un tiempo de cierre tc=1.5 segundos. La ley de
cierre est dada por la siguiente ecuacin:
(
Por otra parte:

(
t
t=t0 +4t
t
t=t0 +3t
t
t=t0 +2t
t
t=t0 +t
t
t=t0
x
Seccin: 0
X=0
5
X=L
Puntos interiores (t>t0)

Condiciones para el flujo permanente
inicial (t=t0)
Figura 5.2. Malla caracterstica

154
Tabla de clculo
Tabla 5.2. Tabla de clculo
Seccin (i)
t (s)
t=t0
t=t0+t
h0
0
Q0
Cn0
Cp1
h1
Q1
Cn1
Cp2
h2
Q2
Cn2
Cp3
h3
Q3
Cn3
Cp4
h4
Q4
Cn4
Cp5
t=1.0
5. Dibujar la curva X-h, en m, para los instantes de tiempo 0, 0.4, 0.8, 1.2 y 1.5 segundos.
6. Dibujar la curva T-h, en s y m, respectivamente, en las secciones 1, 2, 3, 4 y 5.
7. Comparar la carga mxima medida en el dispositivo con la obtenida en el inciso 5, en la seccin
5.
Instrumentos de medicin
Flexmetro
Referencias bibliogrficas
1. Sotelo A., G., Hidrulica General Vol. I, Ed. Limusa, Mxico 1990.
2. Mataix, C., Mecnica de Fluidos y Mquinas Hidrulicas, Ed. Harla, Mxico, 1990.
3. Gardea V., H., Aprovechamientos Hidrulicos y de Bombeo, Facultad de Ingeniera, UNAM.
4. Mataix, C. Mecnica de Fluidos y Mquinas Hidrulicas, Ed. Harla, Mxico, 1995.
5. Mancebo, U. Teora del golpe de ariete y sus aplicaciones en ingeniera hidrulica. Ed.
Limusa, Mxico, 1997.
6. Chaudhry, M. Hanif. Applied Hidraulic Transients.
Cuestionario
El cuestionario tendr de 5 a 10 preguntas relacionadas con el desarrollo y memoria de clculo y
variarn cada semestre.

155
5
h5
Q5

experimental propuesto.
Con base en la hoja de clculo realizada en el captulo 3 para dar solucin al caso de estudio, se
realiz la solucin para el dispositivo experimental propuesto. A diferencia del dispositivo
experimental existente, el dispositivo propuesto no cuenta con un aforador de pared delgada, por lo
cual el gasto en condiciones iniciales en flujo permanente se determinar por medio de un aforo
volumtrico. As mismo, el dispositivo propuesto cuenta con una vlvula solenoide para efectuar el
cierre instantneo y provocar el golpe de ariete. Para realizar las mediciones de la variacin en la
presin, se propone la instalacin de dos sensores de presin aguas arriba y aguas abajo de la
tubera, los cuales desplegarn las mediciones por medio de una tarjeta de adquisicin de datos
hacia una computadora.
En seguida se indica el procedimiento para realizar la prctica de laboratorio, as como la solucin
por medio del mtodo de las caractersticas.

156
5.2.1. MEMORIA DE CLCULO

Medicin y registro de lecturas.
Abrir la vlvula de alimentacin general a la tubera en espiral y mantener en posicin abierta la
vlvula solenoide.
Una vez establecido el gasto en la tubera, medir con el manmetro de cartula y con el sensor de
presin 1 la carga de presin establecida (hp0). El valor esperado o supuesto de la carga de presin
inicialmente establecida es de 1.00 kg/cm2 (10.00 m de columna de agua).
Realizar tres mediciones de volumen y tiempo para calcular el gasto correspondiente con la
siguiente ecuacin.
Donde
V = volumen aforado
t = tiempo de aforo
y calcular el gasto promedio como
El valor del gasto promedio esperado es de 0.1 l/s 1 x 10-4 m3/s

Poner en posicin cerrada la vlvula solenoide y observar la evolucin de la presin en los dos
sensores de presin.
Registrar los datos obtenidos en la tabla 5.3.

157
Tabla 5.3. Registro de lecturas
Tubera de
descarga
Aforador
volumtrico
Sensores de
presin
hp0 =
10.00 m.c.a.
m3
s
m3
s
m3
s
V1 =
t1 =
V2
t2 =
V3 =
t3 =
hps1 =
hps2 =
hpmxs2 =
m.c.a.
m.c.a.
m.c.a.
Clculo de celeridad de la onda de presin en la tubera.

donde:
densidad del agua, = 1,000 Kg/m3;
mdulo de elasticidad de la tubera de cobre,
dimetro de la tubera, D = 0.0127 m;

espesor de la tubera, e = 0.00119 m.
Por lo tanto la celeridad es:
)(
)(
)
)
Periodo T de la onda de presin.
donde:
L = 61.00 m
a = 1,350.09 m/s
Por lo tanto

158
Determinacin del factor de friccin
donde
0
rugosidad del tubo (nuevo), en mm;
coeficiente que depende del grupo en el que se clasifique el agua que va a escurrir, segn la
tabla 3.4; a = 0.025

t
nmero de aos de servicio de la tubera; t = 15 aos
rugosidad del conducto, despus de t aos de servicio, en mm.
Segn la tabla 3.3, para la tubera lisa de cobre, la rugosidad absoluta (tubera nueva de cobre) es
Sustituyendo los valores en la ecuacin 3.5:

(
De la ecuacin de continuidad, la velocidad es
Si se considera una viscosidad cinemtica del fluido , se obtiene de la figura 3.4, para una
temperatura promedio de 20C.
Una vez determinada la velocidad y la viscosidad cinemtica del fluido, el nmero de Reynolds
tiene el siguiente valor
(
A partir del diagrama de Moody (figura 3.5), y con la rugosidad relativa , se determin el factor
de friccin correspondiente a la tubera del dispositivo experimental.

159
Del diagrama de Moody con el valor de la relacin anterior y el nmero de Reynolds obtenido
anteriormente, el valor del factor de friccin
Una vez determinado el factor de friccin, se realiz el clculo de las prdidas por friccin con la
ecuacin de Darcy- Weisbach.
Algunas caractersticas geomtricas e hidrulicas de la tubera son:
D (tubera) = 0.0127 m
A (tubera) = 0.00013 m2
Q0 = 0.00010 m3/s
Rh = 0.00318 m
V = 0.78941 m/s
V2/2g = 0.03176 m
Considerando 5 tramos de tubera iguales a razn de 12.2 m cada uno, la prdida por friccin
calculada con la ecuacin de Darcy-Weisbach, es
Utilizando las prdidas por friccin calculadas con la ecuacin de Darcy-Weisbach y tomando en
cuenta la carga de velocidad, se obtuvieron las siguientes cargas de presin
Tabla 5.4. Prdidas totales Darcy Weisbach

Tramo 1
1.8002
Tramo 2
3.6003
Tramo 3
5.4005
Tramo 4
7.2007
Tramo 5
9.0009

160
Tabla 5.5. Carga piezomtrica
Carga Piezomtrica, en m
Calculado
H0
hpz0
hpz1
hpz2
hpz3
hpz4
hpz5
10
9.9682
8.1681
6.3679
4.5677
2.7675
0.9674
Ley para maniobra de cierre lineal.
Una vez determinadas las condiciones iniciales a flujo permanente, se realiza la solucin mediante
el mtodo de las caractersticas.
Se requiere conocer el intervalo de tiempo y los puntos de inters, para lo cual aplicamos la
condicin de estabilidad de Courant, ecuacin (2.71), de la cual podemos decir lo siguiente
Donde x se define de la siguiente forma.
donde:
L longitud de la tubera, en m
N nmero mnimo entero de tramos en que debe dividirse el conducto
Si N = 5 tramos y L = 61.00 m, se tiene
As, t resulta:

161
De esta manera
Con esto se cumple con la condicin de estabilidad de Courant.

Una vez conocidos t y x, se define la malla caracterstica, la cual se encuentra representada en la
figura 5.3. Esta malla est compuesta de cinco secciones, de las cuales, la seccin 0 representa la
condicin de frontera aguas arriba (tubera de red), la seccin 5 representa la condicin de frontera
aguas abajo (vlvula) y las secciones 1, 2, 3 y 4, representan los puntos interiores de la malla.
t
t=t0 +4t
t
t=t0 +3t
t
t=t0 +2t
t
t=t0 +t
t
t=t0
x
Seccin:
x
1
x
2
x
3
x
4
X=0
5
X=L

inicial (t=t0)

162
Planteamiento de las ecuaciones caractersticas

Una vez definida la malla caracterstica, se realiz el clculo de las constantes necesarias para las
ecuaciones caractersticas.
De la ecuacin 2.57, obtenemos el valor de la constante Ca. Sustituyendo valores, result
Existe otra constante que est en funcin del factor de friccin obtenido, la cual definimos con la
siguiente ecuacin
Sustituyendo los valores se tiene

(
(
)(
)(
)
)
Con estos valores, se hizo el planteamiento de las constantes Cp y Cn, las cuales pueden ser
conocidas para cada instante de tiempo y que forman parte de las ecuaciones caractersticas. Estas
constantes quedaron definidas con las ecuaciones 2.55 y 2.56, y ahora, sustituyendo los valores se
tiene
|
(5.1)
(5.2)
Ahora, la ecuacin caracterstica positiva (ecuacin 2.53), queda de la siguiente forma

(5.3)
y la ecuacin caracterstica negativa (ecuacin 2.54), queda como
(5.4)
Los puntos interiores de la malla (figura ) se definen a partir de la ecuacin 2.58
(
(5.5)
y de la ecuacin (3.23)

163
Sustituyendo el valor de la constante Ca en la ecuacin anterior, resulta

(
(5.6)
Las ecuaciones 5.5 y 5.6, definen el gasto Q y la carga de presin h en los puntos interiores de la
malla.
Condicin de frontera aguas arriba (Seccin 0)
La condicin de frontera aguas arriba est representada por una tubera con una presin constante y
las ecuaciones que determinan el gasto Q y la carga h en esta seccin son las ecuaciones 2.62 y
2.63. Para el dispositivo experimental utilizado, se determin el valor del coeficiente de prdidas a
la entrada de la tubera con un valor de k = 0.5 (ver prdidas locales), lo que lleva a calcular el valor
de la constante k1 (ecuacin 2.64). Sustituyendo el valor de k y Ca en la ecuacin (2.64), dio como
resultado
Sustituyendo las constantes en la ecuacin (2.62), para obtener la carga h en la condicin de
frontera aguas arriba
(
)(
(5.7)
y sustituyendo constantes en la ecuacin (2.63), result
)[
(
(
)
)( )]
(5.8)
Condicin de frontera aguas abajo (Seccin 5)

(
(5.9)

164
Esquema de diferencias finitas

Se sabe que para
se presentan las condiciones iniciales a flujo permanente. Estos datos se
pueden ver en la tabla 5.6.

Seccin
Para el instante
hpz
(m)
(m3/s)
(m)
0.00
0.0001
9.9682
12.20
0.0001
8.1681
24.40
0.0001
6.3679
36.60
0.0001
4.5677
48.80
0.0001
2.7675
61.00
0.0001
0.9674
, es decir para
, se presenta la forma en que deben ser
calculados el gasto Q y la variacin de la carga h con respecto a la posicin en la tubera. Los

subndices que aparecern en las ecuaciones, por ejemplo Q1,0, representan la seccin (i) y el
instante de tiempo (t) respectivamente que se deben utilizar para realizar el clculo.
Condicin de frontera aguas arriba
Seccin 0: X = 0, i=0
|
(
)[
(
(
)|
)( )]
)[
(
(
)(
)
]
)
m.c.a.
165
Puntos interiores
Seccin 1: X=12.2 m, i=1
|
(
|
(
|
(
)|
)|
)|
)|
|
(
Seccin 2: X=24.4 m, i=2

|
(
|
(
|
(
(
(
|
(
)
)

166
)
)
Seccin 3: X=36.6 m, i=3

|
(
|
(
|
)
)|
)|
)|
|
(
Seccin 4: X=48.8 m, i=4

|
(
|
(

167
)|
)|
(
(
)
)
Condicin de frontera aguas abajo

Seccin 5: X=61.0 m, i=5
|
(
|
(
009
(
(
(
)
)
De esta forma, se determinan las variables Qp y hp a lo largo de la tubera para cada instante de
tiempo.

168
A continuacin se presentan las grficas donde se puede observar la variacin de la carga de presin
(h) con respecto a la posicin en la tubera (X) y el instante de tiempo que se presenta (t).

X-h
tc = 1.0s
h (mca)
18.0
16.0
14.0
12.0
t= 0 s
10.0
t = 0.2 s
8.0
t = 0.4 s
6.0
t = 0.6 s
4.0
t = 0.8 s
2.0
t = 1.0 s
0.0
0.0
10.0
20.0
30.0
40.0
50.0
60.0
70.0
X (m)
Carga Piezomtrica
Seccin 1 (X=12.2 m)
hpz ( m.c.a. )
18.60
16.60
14.60
12.60
10.60
8.60
6.60
4.60
2.60
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.60
Tiempo ( s )

169
Carga Piezomtrica
Seccin 2 (X=24.4 m)
hpz ( m.c.a. )
18.60
16.60
14.60
12.60
10.60
8.60
6.60
4.60
2.60
0.80
0.90
1.00
0.80
0.90
1.00
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.60
Tiempo ( s )
Carga Piezomtrica
Seccin 3 (X=36.6 m)
hpz ( m.c.a. )
18.40
16.40
14.40
12.40
10.40
8.40
6.40
4.40
2.40
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )

170
Carga Piezomtrica
Seccin 4 (X=48.8 m)
hpz ( m.c.a. )
18.40
16.40
14.40
12.40
10.40
8.40
6.40
4.40
2.40
0.80
0.90
1.00
0.80
0.90
1.00
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.40
Tiempo ( s )
Carga Piezomtrica
Seccin 5 (X=61.0 m)
hpz ( m.c.a. )
25.00
20.00
15.00
10.00
5.00
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.00
Tiempo ( s )

171
X-h
tc = 0.5 s
h (mca)
25.0
20.0
t= 0 s
15.0
t = 0.1 s
t = 0.2 s
10.0
t = 0.3 s
t = 0.4 s
5.0
t = 0.5 s
0.0
0.0
10.0
20.0
30.0
40.0
50.0
60.0
70.0
X (m)
Carga Piezomtrica
Seccin 1 (X=12.2 m)
hpz ( m.c.a. )
14.60
13.60
12.60
11.60
10.60
9.60
8.60
0.50
0.45
0.40
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
7.60
Tiempo ( s )

172
Carga Piezomtrica
Seccin 2 (X=24.4 m)
hpz ( m.c.a. )
18.60
16.60
14.60
12.60
10.60
8.60
6.60
0.40
0.45
0.50
0.40
0.45
0.50
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
4.60
Tiempo ( s )
Carga Piezomtrica
Seccin 3 (X=36.6 m)
hpz ( m.c.a. )
20.40
15.40
10.40
5.40
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
0.40
Tiempo ( s )

173
Carga Piezomtrica
Seccin 4 (X=48.8 m)
hpz ( m.c.a. )
20.40
15.40
10.40
5.40
0.40
0.45
0.50
0.40
0.45
0.50
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
0.40
Tiempo ( s )
Carga Piezomtrica
Seccin 5 (X=61.0 m)
hpz ( m.c.a. )
30.00
25.00
20.00
15.00
10.00
5.00
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
0.00
Tiempo ( s )

174
Las condiciones iniciales en cuanto a la presin inicial en la tubera y el gasto establecido

corresponden a valores de 1.00 Kg/cm2 y 0.0001 m3/s. Dichos valores son supuestos y adecuados a
las caractersticas de la tubera.
De los resultados obtenidos, se observa la variacin de la carga de presin respecto a la posicin en
la tubera y respecto a cada instante de tiempo. Para el tiempo de cierre propuesto de tc=1.0
segundos, la carga mxima calculada en la seccin 5, tuvo una presin mxima de 16.78 metros de
columna de agua.
Para el tiempo de cierre propuesto de tc=0.5 segundos, la carga mxima calculada en la seccin 5,
tuvo una presin mxima de 25.76 metros de columna de agua.
La variacin en la onda de presin calculada deber ser comparada con la variacin de la onda de
presin registrada por los sensores de presin.
Se presenta a continuacin la propuesta para la realizacin de la prctica de laboratorio Golpe de
ariete aplicando el mtodo de las caractersticas, para un tiempo de cierre de 1.0 segundos.

175
Objetivo
Observar la variacin de la carga de presin a lo largo de una tubera, provocada por el cierre y
apertura de una vlvula colocada en su extremo final.
Calcular y analizar la variacin de la carga de presin generada en una tubera por el cierre de una
vlvula, aplicando el Mtodo de las Caractersticas.
Antecedentes
Flujo permanente en conductos a presin
Flujo transitorio en conductos a presin
Teora del golpe de ariete
Celeridad de una onda de presin
Clasificacin de maniobras de vlvula
Desarrollo
Medicin de condiciones iniciales en flujo permanente.
1. Abrir la vlvula de alimentacin general y la vlvula solenoide hasta que se establezca el flujo.
2. Determinar el gasto mediante aforo volumtrico.
3. Conectar los sensores de presin.
4. Colocar en posicin de cerrado la vlvula solenoide y observar la evolucin de la presin en los
dos sitios de ubicacin de los sensores de presin.
5. Finalmente, realizar una maniobra de apertura y observar el comportamiento de la presin en los
dos sitios de ubicacin de los sensores de presin.

176
COMPUTADORA
TARJETA
DE
ADQUISICIN
DE DATOS
MANMETRO
Sensor de
presin
1
Sensor de
presin
2
ESPIRAL DE
TUBERA DE COBRE
D = 1.27 cm
L = 61 m
e = 0.119 cm
CONEXIN
A RED
HIDRULICA
Descarga libre
Q
Vlvula de
globo
Vlvula
solenoide
Sensor de
presin 2
Q
Espiral de
tubera de cobre
D = 1.27 cm
L = 61 m
e = 0.119 cm
Qred
Q
Vlvula
de globo
Sensor de
presin 1
Abastecimiento
de la red
hidrulica
Figura 5.11. Arreglo de la instalacin hidrulica propuesta

177
Registro de datos
Tabla 5.7. Registro de datos
Tubera de
descarga
Aforador
volumtrico
Sensores de
presin
H0 =
V1 =
t1 =
V2
t2 =
V3 =
t3 =
Hs1 =
Hs2 =
Hmxs2 =
m3
s
m3
s
m3
s
M
M
M
Memoria de clculo
1. Una vez tomadas las tres mediciones de volumen y tiempo, calcular los gastos Q1, Q2 y Q3 de la
tubera, en m3/s, segn la ecuacin
Donde
V = volumen aforado
t = tiempo de aforo
y calcular el gasto promedio como
2. Calcular la celeridad de la onda a, en m/s, de la onda de presin en la tubera, con la ecuacin:

178
donde:
a
celeridad de la onda de presin, en m/s;
densidad del agua,
Ev
;
Et
mdulo de elasticidad de la tubera,

;
dimetro de la tubera,
espesor de la tubera,
3. Calcular el periodo T, en segundos, de la onda de presin, como:
, donde L = 61.00 m es la
longitud de tubera.
4. Solucin aplicando el Mtodo de las Caractersticas
Ecuaciones caractersticas:
( )
( )
en las cuales:
|
( )
( )
donde:
Considerar f=0.059, t=0.009036 segundos y x=12.2 m. Ver figura 5.12

179
Puntos interiores:
(
( )
)
( )
Condiciones de frontera:
Condicin de frontera aguas arriba. Tanque con nivel constante

(
( )
donde:
(
( )
Considerar k = 0.5 (Coeficiente de prdida a la entrada de la tubera).

(
( )
Considerando un cierre lineal (tc = 1.0 s)

(
(
Por otra parte:

(

180
t
t=t0 +4t
t
t=t0 +3t
t
t=t0 +2t
t
t=t0 +t
t
t=t0
x
Seccin:
X=0
X=L

inicial (t=t0)
Tabla de clculo
Tabla 5.8. Tabla de clculo
Seccin (i)
t (s)
t=t0
t=t0+t
h0
0
Q0
Cn0
Cp1
h1
Q1
Cn1
Cp2
h2
Q2
Cn2
Cp3
h3
Q3
Cn3
Cp4
h4
Q4
Cn4
Cp5
t=1.0
5. Dibujar la curva X-h, en m, para los instantes de tiempo 0, 0.4, 0.8 y 1.0 segundos.
6. Dibujar la curva T-h, en s y m, respectivamente, en las secciones 1, 2, 3, 4 y 5.
7. Comparar las curvas desplegadas por los sensores de presin, con la obtenida en el inciso 5, en la
seccin 5.

181
5
h5
Q5
Instrumentos de medicin
Probeta graduada
Cronmetro
Referencias bibliogrficas
5. Mancebo, U. Teora del golpe de ariete y sus aplicaciones en ingeniera hidrulica. Ed. Limusa,
Mxico, 1997.
Cuestionario
El cuestionario tendr de 5 a 10 preguntas relacionadas con el desarrollo y memoria de clculo y
variarn cada semestre.

182
CAPTULO 6. CONCLUSIONES
CONCLUSIONES
Se definieron los conceptos principales que se utilizaron como base para comprender el
fenmeno transitorio conocido como golpe de ariete.
Se investig sobre la teora de dos mtodos de solucin para el golpe de ariete, tales como,
el mtodo de Allievi y el mtodo de las caractersticas.
Se aplic un mtodo numrico de solucin para el fenmeno transitorio golpe de ariete en el

dispositivo experimental instalado en el Laboratorio de Hidrulica de la Facultad de
Ingeniera de la UNAM. Dicho mtodo consisti en el mtodo de las caractersticas y para
fines de comparacin se utiliz tambin el mtodo de Allievi.
Para la aplicacin del mtodo de las caractersticas se obtuvo la ley para maniobra de cierre
no uniforme de la vlvula de esfera especficamente.
Se obtuvieron resultados del fenmeno transitorio golpe de ariete provocado en el

dispositivo experimental tanque tubera vlvula de esfera con el mtodo de Allievi y
con el mtodo de las caractersticas, para diferentes tiempos de cierre de la vlvula de
esfera.
Se propuso la instalacin de instrumentos de medicin en el dispositivo tanque tubera

vlvula de esfera, tales como dos sensores de presin con un rango de medicin de 0 a 50
psi, alimentados mediante una fuente de poder, una tarjeta de adquisicin de datos, y un
software para desplegar las variaciones de la presin en una computadora.
Se realiz una propuesta de prctica de laboratorio para que los alumnos de la asignatura
Hidrulica de Mquinas, conozcan y empleen el mtodo de las caractersticas como mtodo
de solucin para el golpe de ariete provocado en el dispositivo experimental tanque
tubera vlvula de esfera.

183
Se realiz una propuesta para la rehabilitacin del dispositivo experimental instalacin

hidrulica tubera en espiral vlvula, para poder representar el fenmeno golpe de ariete
en este dispositivo.
Se realiz la solucin mediante el mtodo de las caractersticas para el caso del golpe de
ariete provocado en el dispositivo instalacin hidrulica tubera en espiral vlvula a
partir de datos supuestos, ya que el dispositivo no se encuentra en funcionamiento.
Se propuso la instalacin de instrumentos de medicin y control para el dispositivo

experimental que consta de instalacin hidrulica tubera en espiral vlvula. Los
instrumentos de medicin se resumen en: dos sensores de presin con un rango de medicin
de 0 50 psi, los cuales sern alimentados por una fuente de poder. Dicho rango de presin
es apto para soportar las presiones calculadas para este caso. As mismo se propone el uso
de una tarjeta de adquisicin de datos, para captar las seales enviadas por los sensores de
presin. Dicha tarjeta viene acompaada de un software para poder desplegar en tiempo
real la variacin de la onda de presin.
En cuanto al instrumento de control del dispositivo, se uso una vlvula solenoide, la cual al
momento de cerrarse provocar el golpe de ariete. Esta vlvula ser controlada mediante la
tarjeta de adquisicin de datos utilizando como interfase un mdulo de relevadores, para
realizar la activacin de la vlvula solenoide.
El mtodo de las caractersticas resulta ser ms representativo de lo que realmente pasa en

las conducciones, ya que se pueden incluir variables como condiciones de frontera tales
como tanques de carga constante, bombas, leyes de maniobra de cierre no lineales, as como
poder tomar en cuenta los cambios de dimetro a lo largo de la conduccin, as como las
prdidas por friccin y locales que se presentan. Dichas variables no se pueden incluir en
otros mtodos como el mtodo de Allievi o el mtodo de Joukowsky.
Finalmente, el mtodo de las caractersticas resulta prctico y sencillo en su aplicacin. Se

pudieron obtener resultados muy apegados a lo medido, por lo cual se considera que la
convergencia del mtodo es aceptable.

184
Con la realizacin de la memoria de clculo del caso de estudio del captulo 3, se pudo
plantear una propuesta para prctica de laboratorio, teniendo como solucin el mtodo de
las caractersticas.
Con la realizacin de la memoria de clculo para el caso del dispositivo propuesto

instalacin hidrulica tubera en espiral vlvula, se pudo plantear una propuesta para
prctica de laboratorio, teniendo como solucin el mtodo de las caractersticas.
Por ltimo, las dos propuestas de prctica de laboratorio pretenden formar parte del material
didctico para la imparticin de la prctica Golpe de Ariete que actualmente se imparte en
el laboratorio de hidrulica de la Facultad de Ingeniera, UNAM.

185
BIBLIOGRAFA
BIBLIOGRAFA
5. Mancebo, U. Teora del golpe de ariete y sus aplicaciones en ingeniera hidrulica. Ed. Limusa,
Mxico, 1997.
7. Manual de Agua Potable, Alcantarillado y Saneamiento. Fenmenos Transitorios en Lneas de
Conduccin. Comisin Nacional de Agua, Mxico, 2007
8. Moiss Berezowsky Verduzco. Cmaras de Oscilacin. Anlisis y Diseo Hidrulico.
9. Jos Luis Snchez Bibriesca. Mtodo de las Caractersticas en Problemas de Frontera Variable.
Ingeniera Hidrulica en Mxico. Enero/Abril de 1987.

186

Tesis M Bonilla

Cargado por

Copyright:

Formatos disponibles

Tesis M Bonilla

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tesis M Bonilla

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL AUTNOMA

TESIS PARA OBTENER EL TTULO DE

MARIO DAVID BONILLA SNCHEZ

M.I. NIKTE NORMA OCAMPO GUERRERO

MARIO DAVID BONILLA SNCHEZ

MARIO DAVID BONILLA SNCHEZ

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

3.9. Solucin aplicando el mtodo de Allievi

4. PROPUESTAS DE DISPOSITIVOS EXPERIMENTALES

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

Finalmente, en el captulo cinco se presenta una propuesta de prontuario de prctica de laboratorio,

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

Streeter, V. L. y Wylie, E. B. bajo el ttulo Hydraulic Transients (Transitorios hidrulicos), que

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

Prctica Golpe de Ariete en el laboratorio de hidrulica de la Facultad de Ingeniera,

Actualmente, la carrera de Ingeniera Civil, a travs del Departamento de Ingeniera Hidrulica,

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

CAPTULO 2. CONCEPTOS FUNDAMENTALES

CONCEPTOS FUNDAMENTALES DE LA TEORA DEL GOLPE DE ARIETE

Ecuacin de continuidad para un conducto

Principio de conservacin de la masa

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

CAPTULO 2. CONCEPTOS FUNDAMENTALES

Cantidad neta de masa que

A es la superficie de frontera del volumen,

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

CAPTULO 2. CONCEPTOS FUNDAMENTALES

En consecuencia, de acuerdo con el principio de conservacin de la masa, se puede escribir

), la ecuacin (2.2) se reduce a

Al desarrollar las derivadas parciales de la ecuacin anterior y ordenar trminos, se obtiene

As, tambin las derivadas totales de A y se pueden escribir como

Sustituyendo las expresiones anteriores en la ecuacin (2.2b) y dividiendo entre

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

CAPTULO 2. CONCEPTOS FUNDAMENTALES

Ecuacin dinmica para un conducto

Estas fuerzas aplicadas en la direccin + x, son las siguientes:

b) La fuerza de superficie resultante debida a la resistencia al movimiento, aplicada al

CAPTULO 2. CONCEPTOS FUNDAMENTALES

multiplicada por la componente tangencial de la aceleracin, lo cual permite

y se divide la ecuacin anterior entre

Para una determinada seccin del conducto el valor

representa la energa total

puede ser definido como | | para tomar en cuenta la inversin en la

CAPTULO 2. CONCEPTOS FUNDAMENTALES

Descripcin del fenmeno en un conducto por gravedad

Figura 2.3. Conducto con seccin transversal y espesor constantes

PROPUESTA DE INSTALACIN EXPERIMENTAL EN EL LABORATORIO DE HIDRULICA PARA LA

CAPTULO 2. CONCEPTOS FUNDAMENTALES

(e) Condiciones para t = 2L / a + t

(f) Condiciones para t = 3L / a

(c) Condiciones para t = L / a + t

(g) Condiciones para t = 3L / a + t

(b) Condiciones para t = L / a

(a) Condiciones para t + t

(d) Condiciones para t = 2L / a

(h) Condiciones para t = 4L / a

Cuando el frente de onda llega al depsito en un tiempo

, la totalidad del conducto se