Métodos numéricos aplicados al cálculo hidráulico en canales de regadío de Bagua
Numerical methods applied to the hydraulic calculation of irrigation canals in Bagua
Métodos numéricos aplicados para o cálculo hidráulico dos canais de irrigação em Bagua
Ronald Omar Estela Urbina1*, Jhon Jairo Danducho Paati1 Segundo Marcial Chiclote Alcalde1, Fernando Alain Incio Flores1 Nemesio Santamaría Baldera1, Roger Álvaro Fernández Villarroel1 Carlos Carcausto Quispe2, Giovanna Guzmán Cáceres3 María Elena Cárdenas León4, Daniel Jesús Castro Vargas4
El presente artículo evidencia la investigación realizada, ponderando la importancia de la hidráulica en la formación académica del futuro ingeniero civil, así mismo teniendo en cuenta la escasa infraestructura de riego en la provincia de Bagua, que es eminentemente agrícola. El objetivo fue la revisión de los factores hidráulicos de los distintos conductos abiertos existentes, ya sean naturales o artificiales, así como la elaboración de modelos matemáticos que permitan simular numéricamente el comportamiento hidráulico. Esta investigación se realizó en el curso de Métodos Numéricos de la carrera de Ingeniería Civil de la Universidad Nacional Intercultural “Fabiola Salazar Leguía” de Bagua –UNIFSLB durante el año 2020 y 2021. La metodología consistió en tres fases: La primera fase consistió en la obtención de datos en campo, en la segunda fase se realizó cálculos informáticos utilizando el lenguaje de programación Python, y en la última fase metodológica se procedió a elaborar e implementar modelos o métodos matemáticos. Los resultados obtenidos para el caudal y la velocidad difieren aproximadamente en 3% del valor teórico respecto al valor medido en campo, de la misma forma se implementó programas para el cálculo hidráulico. Finalmente se codificó los métodos numéricos para la simulación numérica y gráfica.
Palabras claves: Método numérico, canal hidráulico, Manning, Python.
This article evidences the research carried out, considering the importance of hydraulics in the academic training of the future civil engineer, also taking into account the scarce irrigation infrastructure in the province of Bagua, which is eminently agricultural. The objective was the review of the hydraulic factors of the different existing open conduits, whether natural or artificial, as well as the elaboration of mathematical models that allow numerically simulating the hydraulic behavior. This research was carried out in the Numerical Methods course of the Civil Engineering career of the National Intercultural University "Fabiola Salazar Leguía" of Bagua -UNIFSLB during the year 2020 and 2021. The methodology consisted of three phases: The first phase consisted of field data collection, in the second phase computer calculations were performed using the Python programming language, and in the last methodological phase we proceeded to develop and implement models or mathematical methods. The results obtained for flow and velocity differ by approximately 3% from the theoretical value with respect to the value measured in the field, and programs for hydraulic calculation were also implemented. Finally, the numerical methods for numerical and graphical simulation were codified.
Keywords: Numerical method, hydraulic channel, Manning, Python.
DOI: https://doi.org/10.55996/dekamuagropec.v3i1.70
1 Universidad Nacional Intercultural Fabiola Salazar Leguía de Bagua. Perú, correo: restela@unibagua.edu.pe; jdanduchop@unibagua.edu.pe; schiclotea@unibagua.edu.pe; fincio@unibagua.edu.pe; nsantamaria@unibagua.edu.pe; rfernandez@unibagua.edu.pe
2Universidad Nacional del Altiplano. Perú, correo: ccarcausto@unap.edu.pe
3Universidad Nacional San Antonio Abad del Cuzco. Perú, correo: giovanna.guzman@unsaac.edu.pe
4Universidad Nacional Autónoma de Chota. Perú, correo: mcardenas@unach.edu.pe; djcastrov@unach.edu.pe
Este artigo mostra a pesquisa realizada, considerando a importância da hidráulica na formação acadêmica do futuro engenheiro civil, bem como levando em conta a escassa infra-estrutura de irrigação na província de Bagua, que é eminentemente agrícola. O objetivo era rever os fatores hidráulicos dos diferentes condutos abertos existentes, sejam naturais ou artificiais, bem como a elaboração de modelos matemáticos para simular numericamente o comportamento hidráulico. Esta pesquisa foi realizada no curso de Métodos Numéricos do curso de Engenharia Civil da Universidade Nacional Intercultural "Fabiola Salazar Leguía" de Bagua -UNIFSLB durante os anos de 2020 e 2021. A metodologia consistia em três fases: A primeira fase consistiu na coleta de dados de campo, na segunda fase foram realizados cálculos computadorizados utilizando a linguagem de programação Python, e na última fase metodológica procedemos ao desenvolvimento e implementação de modelos ou métodos matemáticos. Os resultados obtidos para o fluxo e velocidade diferem em aproximadamente 3% do valor teórico em relação ao valor medido no campo, e da mesma forma, foram implementados programas para o cálculo hidráulico. Finalmente, os métodos numéricos para a simulação numérica e gráfica foram codificados.
Palavras-chave: Método numérico, calha hidráulica, Manning, Python.
Donde V representa la velocidad del fluido, Rh es el radio hidráulico, S pendiente del canal, n es el coeficiente de rugosidad.
Haciendo uso del principio de continuidad podríamos definir el caudal (Q) que trastada una determinada sección hidráulica (A):
El volumen de agua que transporta una determinada sección transversal en el tiempo se le denomina caudal como podemos apreciar en la expresión anterior, de la cual se deduce que no solo depende de la máxima sección hidráulica, sino también de las condiciones de revestimiento y conservación de la infraestructura, esto es de la rugosidad que presente el perímetro mojado del canal (Asalde, 2020).
Donde para .
para .
Integración numérica: Regla de Simpson
Los métodos numéricos existentes para estimar el área debajo de una curva son variados, sin embargo, los más usuales y elementales como el trapezoidal y rectangular generan márgenes de error considerables, por lo tanto, planteamos como algoritmo la Regla de Simpson (Figura 1), cuya segmentación de la figura no es lineal, lo que permite un mejor ajuste y cercanía a la figura original (Peña & Rivas, 2019).
Figura 1. Integración numérica: Regla de Simpson
Para esta regla dividiremos el intervalo [a,b] en subintervalos de longitud
Siendo este el valor de la integral donde el valor de sea siempre un número par.
Interpolación numérica: Polinomio de Lagrange
En algunos casos no se dispone de una ecuación o función que permita la manipulación de las variables para su respectivo cálculo, todo lo contrario, se cuenta con pares ordenados
El polinomio de Lagrange se expresa(Carnero, 2018):
Diferenciación numérica
En hidráulica, una de las mediciones directas que podemos realizar es el volumen hídrico por de una determinada sección transversal, tomando como referencia el tiempo, ambas magnitudes físicas que se relacionan con el caudal, teniendo en cuenta que este es la razón temporal de cambio del voumen (Duran, Páez, & García, 2018).
Si conocemos los puntos , los cuales se hallan equi espaciados:
.
.
.
Para el primer intervalo [], se evalúo numéricamente la primera derivada de la función en el punto medio del intervalo asumido:
Forma del canal y ubicación |
Ancho de solera(b) en metros |
Espejo de agua(T) en metros |
Talud(Z) |
Pendiente(S) |
Tirante de agua(y) en metros |
Material de construcción de canal(n) |
Rectangular(Achaguay bajo)*** |
0.95 |
0.95 |
0 |
0.0048 |
0.23 |
0.014 |
Rectangular(Achaguay bajo)** |
0.95 |
0.95 |
0 |
0.0045 |
0.03 |
0.014 |
Rectangular(Achaguay bajo)* |
0.95 |
0.95 |
0 |
0.0056 |
0.1352 |
0.014 |
Parabólico de tierra(Peca Palacios)** |
0 |
2.08 |
0 |
0.0123 |
0.261 |
0.050 |
Parabólico de tierra(Peca Palacios)** |
0 |
3.62 |
0 |
0.0084 |
0.18 |
0.050 |
Trapezoidal(Amoja-Aramango)*** |
1.62 |
2.47 |
1.88 |
0.0032 |
0.80 |
0.015 |
Trapezoidal (San Francisco)** |
0.39 |
0.69 |
1.6
|
0.0024 |
0.24 |
0.013 |
Trapezoidal(San Francisco)* |
0.382 |
0.56 |
0.75 |
0.02619 |
0.12 |
0.015 |
Trapezoidal(La Peca)** |
0.5 |
1.6 |
0.37 |
0.0005 |
0.2 |
0.014 |
Rectangular(La Peca)** |
0.65 |
0.065 |
0 |
0.0008 |
0.5 |
0.014 |
Rectangular(La puntilla – Valencia)** |
0.90 |
0.90 |
0 |
0.005 |
0.06 |
0.013 |
Trapezoidal (La Peca)*** |
0.40 |
0.94 |
1.48 |
0.087 |
0.40 |
0.013 |
Rectangular(La Peca)* |
0.80 |
0.80 |
0 |
0.087 |
0.662 |
0.013 |
Rectangular (San Francisco)** |
0.61 |
0.61 |
0 |
0.0026 |
0.23 |
0.014 |
Rectangular(Toma comunal de Bagua) |
0.95 |
0.95 |
0 |
0.0067 |
0.17 |
0.013 |
Rectangular(La puntilla)*** |
0.94 |
0.94 |
0 |
0.005 |
0.168 |
0.013 |
Factores hidráulicos como espejo de agua(T), tirante de agua(y), ancho de solera(b), Talud(Z) o pendiente(S) se obtenidas mediante mediciones directas, haciendo uso de instrumentos de medición de longitud como cintas métricas, cuya resolución o lectura mínima es el milímetro (0.001 m), así también niveles de burbuja (Figura 3).
Forma del canal |
Longitud del tramo(L) en metros |
Tiempo promedio(t) en segundos |
Estimación del área hidráulica(A) en metros cuadrados |
Velocidad del fluido(v) en mestros por segundo |
Estimación del caudal(Q) en metros cúbicos por segundo |
Parabólico de tierra(Peca Palacios)* |
10.00 |
19.800 |
|
0.52 |
|
Rectangular(La Peca)* |
10.00 |
14.130 |
0.3250 |
0.7077 |
0.2300 |
Trapezoidal(La Peca)* |
7.00 |
6.3915 |
0.2100 |
1.0952 |
0.2300 |
Rectangular(Achaguay bajo)*** |
8.00 |
7.2900 |
0.2185 |
1.0970 |
0.2808 |
Rectangular(Achaguay bajo)** |
5.00 |
10.760 |
0.1254 |
0.4647 |
0.0597 |
Rectangular(Achaguay bajo)* |
9.10 |
7.8790 |
0.1284 |
1.1549 |
0.1483 |
Parabólico de tierra( Peca Palacios)* |
15.00 |
21.440 |
3.2151 |
0.6996 |
2.2527 |
Rectangular (San Francisco)*** |
4.40 |
6.2520 |
0.1403 |
0.7038 |
0.0987 |
Trapezoidal (San Francisco)*** |
5.00 |
6.3080 |
0.2220 |
0.7926 |
0.1760 |
Rectangular(Toma comunal de Bagua)** |
9.00 |
8.2800 |
0.1650 |
1.0869 |
0.1793 |
Rectangular(La puntilla)*** |
6.00 |
4.1250 |
0.1579 |
|
|
Trapezoidal(Amojao-Aramango)*** |
12.00 |
9.7400 |
1.6360 |
1.2320 |
2.0156 |
Figura 6. Código fuente en Python |
Figura 10. Solución numérica para
el tirante de un canal trapezoidal
Una de las fortalezas del lenguaje Python fue la realización de gráficos con mucha facilidad, así como el ingreso y manejo de funciones, lo que permitió codificar los métodos numéricos y la búsqueda de raíces de forma gráfica.
def ejecutar():
def f(x):
f=float(ba.get())*x-float(Q.get())/float(V.get())+float(z.get())*x**2
return f
def d(x):
d=float(ba.get())+2*float(z.get())*x
return d
y=float(yp.get())
print (y)
cont=0
x_x=y
while abs(f(x_x))>0.00001:
Lista_x.insert(str(cont),str(round(x_x, 5)))
Lista_fx.insert(str(cont),str(round(f(x_x), 8)))
Lista_dx.insert(str(cont),str(round(d(x_x), 3)))
cont=cont+1
x_x=x_x-f(x_x)/d(x_x)
Lista_x.place(x=30,y=220)
Lista_fx.place(x=160,y=220)
Lista_dx.place(x=300,y=220)
Alfaro, M., Guerra, R., & Olivares, A. (2020). Evaluación de la profundidad de recursión de la solución analítica de la ecuación de Colebrook-White en la exactitud de la predicción del factor de fricción. Ingeniería, Investigación y Tecnología, 21(4), 1–15. https://doi.org/10.22201/fi.25940732e.2020.21.4.036
Aponte, C. (2019). Diseño y construcción de una canaleta Parshall para el Laboratorio de Hidráulica de la Universidad Santo Tómas (UNIVERSIDAD SANTO TOMÁS). Retrieved from https://repository.usta.edu.co/bitstream/handle/11634/30431/2020carlosaponte.pdf?sequence=6
Arévalo, B., & González, A. (2020). Enseñanza de integrales mediante métodos numéricos pormedio del uso de las TIC. Sinergias Educativas, 5(4), 21–34. Retrieved from http://sinergiaseducativas.mx/index.php/revista/article/view/155/442
Asalde, J. (2020). Diseño del canal de riego con máxima eficiencia hidráulica mediante el software Hec-Ras, en el distrito de Salas
(Universidad San Martín de Porres). Retrieved from https://hdl.handle.net/20.500.12727/6942
Cadena, C., Saltos, M., & Villalta, M. (2018). Comparación de tipos de flujos, para diferentes secciones de canales. (Escuela Superior Politécnica del Litoral). Retrieved from https://www.researchgate.net/publication/322555700_Comparacion_de_tipos_de_flujos_para_diferentes_secciones_de_canales
Caparicona, J. (2020). Elaboración de un programa informático de aplicación para diseño agronómico e hidráulico en el método de riego por goteo. Revista Apthapi, 6(2), 1935–1953. Retrieved from http://ojs.agro.umsa.bo/index.php/ATP/article/view/410
Carnero, E. (2018). Análisis numérico de infiltración de agua en un suelo residual no saturado utilizando la curva de retención de agua. Revista de Investigaciones Altoandinas, 20(4), 439–450. Retrieved from http://www.scielo.org.pe/scielo.php?script=sci_arttext&pid=S2313-29572018000400007
Deza, J., & Castañeda, H. (2020). Hierofonía y matemática del canal Cumbemayo. Cajamarca. Perú. Revista Ciencia y Desarrollo, 23(4), 67–102. https://doi.org/http://revistas.uap.edu.pe/ojs/index.php/CYD/index
Duran, M., Páez, J., & García, P. (2018). Modelado numérico y análisis experimental para flujos en un medio poroso homogéneo a través de suelos. Iteckne, 15(1), 24. https://doi.org/10.15332/iteckne.v15i1.1961
Fatehi, B., Hajikandi, H., Hassanzadeh, Y., & Jamali, S. (2019). Investigación experimental y analítica de los efectos de las células secundarias actuales en las características del salto hidráulico en canales trapezoidales. Tecnologia y Ciencias Del Agua, 10(3), 190–218. https://doi.org/10.24850/j-tyca-2019-03-08
Fernández de Córdova, C., León, A., Rodríguez, Y., Martínez, P., & Meneses, D. (2018). Influencia del método de estimación en el coeficiente de Manning para cauces naturales. Ingeniería Hidráulica y Ambiental, 39(1), 17–31.
Gonzáles, D. (2020). Modelo hidráulico virreinal utilizado en el centro histórico de Lima para abstecimiento de agua. Devenir, 7(14), 91–106. https://doi.org/https://doi.org/10.21754/devenir.v7i14.811
Gutierrez, E. (2020). Consideraciones mínimas para el diseño hidráulico de proyectos de infraestructura menor en canales abiertos (Universidad Peruana Unión). Retrieved from https://drive.google.com/file/d/1_lR8G2xOKhcG-2YVMyc1XXEeJ99GQQ2s/view
Hernández, J., & Martínez, S. (2019). Tránsito de avenidas en vasos: ¿hidrológico o hidráulico? Tecnologia y Ciencias Del Agua, 10(6), 147–177. https://doi.org/10.24850/j-tyca-2019-06-06
Hernández, Y., Rivas, R., & Feliu, V. (2020). Control automático de la distribución de agua en sistemas de riego: revisión y retos. Revista Ingeniería Electrónica, Automática y Comunicaciones I, 41(2), 80–97. Retrieved
Ladino, E., García, C., & García, M. (2020a). Flujo crítico en canales abiertos, solución numérica mediante el método de Newton-Raphson para aplicación Android 4.0. Tecnura, 24(63), 99–114. https://doi.org/10.14483/22487638.16196
Ladino, E., García, C., & García, M. (2020b). LA IMPLICANCIA ECONÓMICA MEDIANTE NEWTON RAPSHON PARA EL DESARROLLO DE UNA APLICACIÓN ANDROID PARA EL DISEÑO DEL DIÁMETRO DE TUBERÍAS A PRESIÓN. Aglala, 11(1), 29–50. Retrieved from http://revistas.curnvirtual.edu.co/index.php/aglala/article/view/1577
Lima, R., Cedeño, J., & Padilla, M. (2020). Aplicación de los métodos numéricos en la enseñanza superior. Revista Sinapsis, 1(16), 1–7. Retrieved from https://www.itsup.edu.ec/sinapsis
Pantaleón, J., Pérez, P., Cordero, R., Torres, A., García, B., Castillo, H., … Rivas, G. (2019). Determinación del coeficiente de rugosidad para encachado usado en el revestimiento de cunetas en la República Dominicana. Ciencia y Sociedad, 38(3), 551–568.
Peña, C., Ramirez, M., & Rivas, E. (2019). La integración numérica en la integral definida : caso de estudio. Revista Espacios, 40(19), 23.
Pineda, J., & Chica, E. (2020). Métodos numéricos para el desarrollo de una turbina hidrocinética tipo Gorlov. Revista UIS Ingenierías, 19(3), 187–205. https://doi.org/10.18273/revuin.v19n3-2020018
Sablón, L., & De Arma, A. (2019). Implementación de métodos numéricos para resolución de sistemas de ecuaciones no lineales. Ciencia & Futuro, 9(1), 74–100. Retrieved from http://revista.ismm.edu.cu/index.php/revista_estudiantil/article/view/1772
Sánchez, E., Gómez, M., & Bladé, E. (2020). Análisis numérico 3D de las características del flujo en un canal curvo. Ingeniería Del Agua, 24(3), 157–169. https://doi.org/10.4995/ia.2020.12276
Segarra, J. (2020). Análisis De Los Métodos Numéricos En Ecuaciones Diferenciales Ordinarias Utilizando Mathematica. Revista Ingeniería, Matemáticas y Ciencias de La Información, 7(13), 13–23. https://doi.org/10.21017/rimci.2020.v7.n13.a72
Terán, J., & Rúa, C. (2018). El método de Newton para raices complejas. Revista EIA, 15(29), 97–108. https://doi.org/https://doi.org/10.24050/reia.v15i29.1131