60 37 PB
60 37 PB
60 37 PB
eISSN: 2477-8990
Volumen 54, Nro. 1, Agosto · Octubre 2024
revistapolitecnica.epn.edu.ec
www.epn.edu.ec
ISSN: 1390-0129
eISSN: 2477-8990
Volumen 54, Nro. 1, Agosto · Octubre 2024
revistapolitecnica.epn.edu.ec
www.epn.edu.ec
TEMÁTICA Y ALCANCE
La Revista Politécnica es una publicación periódica trimestral editada por la Escuela Politécnica Nacional del Ecuador, creada en el año
1961, siendo la primera revista científica ecuatoriana, cuyo objetivo es contribuir al conocimiento científico y tecnológico, mediante la
publicación de estudios científicos relacionados con las áreas de ciencias básicas (Física, Química, Biología y Matemática) e ingenierías
(Química y Agroindustria, Civil y Ambiental, Eléctrica y Electrónica, Geología y Petróleos, Mecánica, y Sistemas). La Revista
Politécnica está dirigida a profesionales e investigadores que trabajan en estos campos del conocimiento.
EDITORA CO-EDITOR
Jenny Gabriela Torres, Ph.D. Luis Alberto Morales, Ph.D.
Escuela Politécnica Nacional Escuela Politécnica Nacional
[email protected] [email protected]
CONSEJO EDITORIAL
Ph.D. José Aguilar Ph.D. Víctor Hugo Hidalgo Ph.D. José Luis Paz
Universidad de los Andes, Venezuela Escuela Politécnica Nacional, Ecuador Universidad Nacional Mayor de San Marcos,
Perú
Ph.D. Hernán Álvarez Ph.D. Majid Khorami (C) Ph.D. Nelson Pérez
Universidad Nacional Colombia, Colombia Universidad Tecnológica Equinoccial, Universidad de los Andes, Venezuela
Ecuador
Ph.D. Santiago Arellano Chalmers Ph.D. Hugo Leiva Ph.D. Franklin Rivas
University of Technology, Suecia Yachay Tech University, Ecuador Universidad Técnica Federico Santamaría,
Chile
Ph.D. Carlos Ávila Ph.D. Francisco León Ph.D. Andrés Rosales
Escuela Politécnica Nacional, Ecuador Universidad de los Andes, Venezuela Escuela Politécnica Nacional, Ecuador
Ph.D. Juan Carlos De los Reyes Ph.D. Rui Pedro Paiva Ph.D. Sebastián Taco
Escuela Politécnica Nacional, Ecuador University of Coimbra, Portugal Escuela Politécnica Nacional, Ecuador
La Revista Politécnica está incluida en diferentes índices, siendo los principales SCOPUS, Scientific Electronic Library Online
(SciELO), Red de Revistas Científicas de América Latina y el Caribe, España y Portugal (Redalyc), Catálogo 2.0 de Latindex,
Directory of Open Access Journals (DOAJ).
Se autoriza la reproducción total o parcial de su contenido siempre y cuando se cite la fuente. Los conceptos expresados son de
responsabilidad exclusiva de sus autores.
Ing. Ricardo Villarroel MSc. Karina Játiva MSc. María Eufemia Cristian Basurto
[email protected] [email protected] Torres Diseñador Gráfico
AUTORIDADES
ESCUELA POLITÉCNICA NACIONAL
Vicerrector de Investigación, Rector Vicerrectora de Docencia
Innovación y Vinculación Tarquino Sánchez, Ph.D. Valentina Ramos, Ph.D.
Marco Santórum, Ph.D.
Cambio climático y ciencia
El cambio climático, impulsado en gran parte por todas las actividades humanas, está
alterando los patrones climáticos globales y afectando la biodiversidad, los ecosistemas y las
comunidades en todo el mundo. Los informes científicos y las evidencias empíricas son
claros: la temperatura global sigue aumentando, los fenómenos meteorológicos extremos se
están volviendo más frecuentes, y el nivel del mar está subiendo a un ritmo alarmante.
Ante esta problemática surge la necesidad de una acción coordinada. No basta con reconocer
la gravedad del problema; es esencial que gobiernos, universidades, empresas, comunidades
y ciudadanos trabajen juntos para implementar soluciones efectivas. La transición hacia
fuentes de energía renovables, la adopción de prácticas sostenibles en la agricultura y la
promoción de tecnologías limpias representan pilares fundamentales en la lucha contra el
cambio climático y el desarrollo de un futuro más sostenible. Estas acciones no solo ayudan
a reducir las emisiones de gases de efecto invernadero, sino que también fomentan un uso
más eficiente de los recursos naturales, contribuyendo a la preservación del medio ambiente
y al bienestar de las generaciones futuras.
En el marco del Día Internacional del Cambio Climático, la Revista Politécnica presenta un
nuevo volumen que incluye a cinco de sus diez artículos que están relacionados con
alternativas para reducir el impacto ambiental del ser humano en la naturaleza, como el uso
de materiales reciclados para la elaboración de concreto, generación de modelos matemáticos
que optimizan la energía eólica, desarrollo de fuentes alternativas de energía renovable para
dispositivos electrónicos portátiles y el estudio de herramientas de optimización de
parámetros y condiciones operativas de sistemas fotovoltaicos.
October 24 marks the International Day of Climate Change. In this context, it is important to
highlight that the climate crisis has reached global dimensions, becoming one of the most
urgent challenges of our time. This day not only marks an opportunity to reflect on the
impacts of climate change, but also to reaffirm our commitment, as academics and
researchers, towards raising, investigating and proposing effective and sustainable solutions.
Climate change, driven in large part by all human activities, is altering global weather
patterns and affecting biodiversity, ecosystems and communities around the world. Scientific
reports and empirical evidence are clear: global temperatures continue to rise, extreme
weather events are becoming more frequent, and sea levels are rising at an alarming rate.
Faced with this problem, there is a need for coordinated actions. It is not enough to recognize
the seriousness of the problem; it is essential that governments, universities, businesses,
communities and citizens work together to implement effective solutions. The transition to
renewable energy sources, the adoption of sustainable practices in agriculture and the
promotion of clean technologies represent fundamental pillars in the fight against climate
change and the development of a more sustainable future. These actions not only help to
reduce greenhouse gas emissions, but also promote a more efficient use of natural resources,
contributing to the preservation of the environment and the well-being of future generations.
Innovation and knowledge are our best allies in this fight. It is imperative that scientific
research and engineering contribute to developing new technologies and approaches that
reduce our carbon emissions. It is also essential that education, dissemination of research in
these aspects and public awareness be strengthened to drive a change in the way we live.
In the framework of the International Day of Climate Change, the Revista Politécnica
presents a new volume that includes five of its ten articles related to alternatives to reduce
the environmental impact of human beings on nature, such as the use of recycled materials
for the production of concrete, the generation of mathematical models that optimize wind
energy, the development of alternative sources of renewable energy for portable electronic
devices and the study of tools for the optimization of parameters and operating conditions of
photovoltaic systems.
The Editorial Board of Revista Politécnica calls on readers and society in general to take an
active role against climate change. Every action counts, and every decision taken can make
a significant difference. Only a joint and sustained effort will protect our planet for future
generations. We invite you to explore the articles in this volume and reflect on how you can
contribute to the necessary change. The climate crisis is not a distant problem; it is an
immediate challenge that requires our attention and determined action.
Contenido
Vol. 54, No. 1
Agosto – Octubre 2024
7
Quispe, Reynaldo; Huamani, Renzon
Modelos Autorregresivos para Predecir la Velocidad del Viento en la Zona Rural
Socabaya de Perú, 2022
Autoregressive Models for Forecasting Wind Speed in the Rural Area of Socabaya, Peru, 2022
15
Gallegos, Álvaro; Ulquiango, Oscar; Lima, Paulina
Comportamiento Hidráulico del Flujo en el Modelo Numérico de un Disipador de Energía
Tipo Vórtice
Hydraulic Flow Behavior in the Numerical Model of a Vortex Energy Sink
25
Rodríguez, Martín; Baes, Mario; Jara, Ricardo; López, Yadira; Orozco, Santos
Closed-loop Speed Control for a Three-Phase Alternating-Current Motor using a Modbus
Network
Control de Velocidad en Malla Cerrada para un Motor de Corriente Alterna Trifásico
utilizando una Red Modbus
33
Flores Berenguer, Isaida; González Haramboure, Yoermes; García Tristá, Jenny; Rosete Suárez,
Alejandro
Artificial Neural Network Model to Predict the Factor of Safety in Earth Dams Subjected
to Rapid Drawdown
Modelo de Redes Neuronales Artificiales para Predecir el Factor de Seguridad en Presas de
Tierra Sometidas a Desembalse Rápido
45
Montenegro, Valery; Hernández, Delvis; Domínguez, Anmary; Vergara, Fidedigna;
Arrocha, Jonatha; González-Valoys, Ana
Estabilización de Suelos Expansivos en Panamá con Mezclas de Residuos de Concha de
Mar
Stabilization of Expansive Soils in Panama with Seashell Residue Mixtures
53
Portilla, Miguel; Ludeña, Raúl; Asanza, Víctor; Dávila, Miguel; Nevarez, Manuel
Diseño e Implementación de un Prototipo Portable de Generación Piezoeléctrica
Design and Implementation of a Portable Piezoelectric Power Generation Prototype
65
Valles, José M.; González-Longatt, Francisco; Riquelme-Rodríguez, José Miguel, Angeles-
Camacho, César
Optimización Simultanea de los Parámetros de Esquemas de Desconexión Automática de
Carga y Convertidores Grid-forming
Simultaneous Optimization of the Parameters of Under-Frequency Load Shedding Schemes and
Grid-forming Converters
77
Monje Yovera, Paul Mauricio; Muñoz Perez, Socrates Pedro; Villena Zapata, Luigi Italo;
García Chumacero, Juan Martín; Diaz Ortiz, Edwin Adolfo; Quispe Osorio, José Luis
Effect of Total Replacement of Quarry Aggregate by River Aggregate on the Strength of
Concrete
Efecto de la Sustitución Total del Agregado de Cantera por Agregado de Río en la Resistencia
del Concreto
87
Quispe, Reynaldo; Juli, Verónica
Uso de Minimización Irrestricta de una Función Espectral para Estimar la Zona Visible
en Matlab 19.0
Use of Unrestricted Minimization of a Spectral Function to Estimate the Visible Zone in Matlab
19.0
97
Carpio, Bryan; Muñiz, Abel
Efectos del Cartón y Cerámica Reciclada en las Propiedades del Concreto Hidráulico
Effects of cardboard and recycled ceramics on the properties of hydraulic concrete
Modelos Autorregresivos para Predecir la Velocidad del Viento en la Zona Rural Socabaya de Perú, 2022 7
1
Universidad Nacional José María Arguedas, Departamento Académico de Ciencias Básicas, Andahuaylas, Perú
2
Universidad Nacional de San Agustín de Arequipa, Departamento Académico de Física, Arequipa, Perú
Resumen: La zona rural del Perú presenta una gran variabilidad de vientos, su desconocimiento debido a la ausencia
de modelos de predicción tiene un efecto desfavorable en la agricultura, infraestructura, seguridad, transporte y
optimización de la energía eólica. Este trabajo obtiene un modelo de media móvil integrada autorregresiva (ARIMA)
de predicción de la rapidez de viento en el lenguaje de programación R para la zona rural del Perú, distrito de Socabaya
de Arequipa. La metodología consiste en el método cuantitativo y técnica documental, con una muestra de 334 datos
de viento del año 2022 de la estación meteorológica Administración Nacional de Aeronáutica y Espacio (NASA)
(POWER, 2023) para el distrito de Socabaya. Utilizando R, se aplicó el método paramétrico: Dicker-Fuller, Levene,
D'Agostino, Kwiatkowski–Phillips–Schmidt–Shin (KPSS) y diferenciación (d=1) para lograr la normalidad y
estacionariedad de los datos. Se analiza la función de autocorrelación simple (ACF) y función de autocorrelación
parcial (Parcial ACF) mediante un proceso de ajuste recursivo, criterio de información de Akaike (AIC) para elegir
el mejor modelo ARIMA de predicción. El resultado obtenido es el modelo de predicción de viento ARIMA (1, 1,
2), con una precisión de error de escala absoluta de media (MASE) de 0,849. Se concluye que el modelo ARIMA
obtenido puede usarse para predecir la velocidad de viento en Socabaya a corto plazo, desde el 29 de noviembre hasta
8 de diciembre del 2022, y su aleatoriedad estaría influenciado por la variabilidad climática y cantidad de datos del
año 2022.
Palabras clave: Simulación, modelos estocásticos, velocidad del viento, modelo ARIMA, predicciones de viento
Keywords: Simulation, stochastic models, wind speed, ARIMA model, wind forecasts
*[email protected] CC 4.0
Recibido: 03/10/2023
Aceptado:24/03/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.01
(Lesage et al., 2022). Los modelos ARIMA fueron planteados midió cada una de las técnicas empleadas (Pilco & Acurio,
por Box y Jenkins en 1970 en su obra "Análisis de Series 2019).
Temporales: Pronóstico y Control", la cual muestra que la
predicción de viento depende de las características La velocidad de viento no solo se utiliza para el diseño e
topográficas de la superficie terrestre, donde una serie instalación de grandes parques eólicos, sino también para
temporal puede ser inferida a partir de un proceso de mantener segura la red eléctrica mediante la predicción de
generación de datos (Uriel, 1985). velocidad de viento utilizando un modelo estocástico ARIMA
y un modelo basado en redes neuronales Feed-Forward
En algunos lugares de África como Tanzania, laderas de Artificial Neural Network (FFANN) o perceptrón multicapa
Kilimanjaro, centro de Sudán los fuertes vientos transportan la (MLP) (Dumitru & Gligor, 2019).
arena a los campos de cultivo secándolas y poniendo en peligro
su producción agrícola (Stigter et al., 2002). El uso de datos climáticos de temperatura, precipitación,
humedad, horas de sol y velocidad de viento permitieron
El viento es la energía renovable (ER) no contaminante más estimar las próximas 24 horas en un sistema de distribución de
utilizada en la Unión Europea y constituye una fuente de agua para riego, a través del uso de métodos ARIMA y redes
energía eólica de bajo costo para almacenar energía en los neuronales (RNCs) (Pulido et al., 2002).
parques eólicos marinos, sin embargo, es vulnerable al cambio
climático global, lo cual afecta la vida silvestre y la Las predicciones de velocidad de viento tienen profundo
supervivencia de las especies que dependen de su condición interés para el desarrollo de turbinas eólicas y las estrategias
reproductiva personal (Ziemba et al., 2017). Las fluctuaciones para construir grandes parques eólicos, pero estas predicciones
de la velocidad del viento son uno de los principales obstáculos se ven afectadas por la turbulencia atmosférica. Por eso, se
para su explotación en la red eléctrica. A una altura de 40 m, utilizaron modelos estocásticos de series temporales ARIMA,
la velocidad media del viento en muchas regiones de Libia está con datos de 20 meses y una variable gaussiana para modelizar
entre 6 m/s y 7,5 m/s. A pesar de esto, hay varias épocas del la variable meteorológica de viento (Sim et al., 2018).
año cuando soplan vientos cálidos y secos en cambio en la
región oriental durante un año la velocidad promedio del La predicción de velocidad de viento a corto plazo, 6-8 horas,
viento fue 8,21 m/s (Badi et al., 2023). para un sistema eléctrico exige cada vez una mejor predicción
en el Báltico. Los métodos de análisis de series de tiempo
A partir del siglo XIX, las investigaciones de predicción han utilizaron modelos ARIMA y sus predicciones se validaron a
aumentado, debido a la disponibilidad de datos meteorológicos través del error cuadrático medio de la raíz (RMSE), error
online y propuestas de modelos de pronóstico ARIMA (p,d,q) porcentual absoluto medio (MAPE) y error absoluto medio
con coeficientes p, d, q (Bórawski et al., 2020; Salman & (MAE) (Grigonytė & Butkevičiūtė, 2016).
Kanigoro, 2021).
En las islas de Baja California, se pronosticó la velocidad de
Las mejoras en la previsión de energía eólica tienen potencial viento mediante un modelo ARIMA y una red neuronal
de incrementar la cantidad de reservas necesarias en sistemas artificial (ANN), con errores obtenidos para pronosticar la
con cantidades significativas de energía eólica, y finalmente tendencia no lineal que el modelo ARIMA no pudo identificar.
reducir el costo de la electricidad en dichos sistemas mediante Además, se utilizaron los errores estadísticos: MAE, error
un análisis de series temporales con modelos ARIMA para medio cuadrado (MSE) y RMSE (Cadenas & Rivera, 2010).
pronosticar la futura producción de energía eólica con datos
históricos (Hodge et al., 2011). En Oaxaca y Metepec, se realizó una comparación de
predicción de velocidad de viento entre el modelo ARIMA y
En China, los fuertes vientos produjeron descarrilamiento de modelos autorregresivos no lineales con entrada exógena
ferrocarriles, por ello, se centró en predicción de alta precisión (NARX), tomando una base de datos de presión, temperatura,
de velocidad de viento para sistemas de alerta de viento de los viento y radiación solar. Se determinó la influencia de las
ferrocarriles. Con ese fin, se utilizó un modelo híbrido que variables meteorológicas sobre la velocidad del viento. Se
combina la descomposición modal (EMD) y ARIMA concluyó que la red neuronal artificial fue 2,3% mejor que el
mediante el modelo RARIMA (Hui et al., 2015). modelo ARIMA, debido a su carácter no lineal (Cadenas et al.,
Nuevas tecnologías para obtener energía eléctrica a partir de 2016).
energías renovables llevaron a buscar metodologías que
permitan identificar zonas con mejores condiciones de Debido a la ausencia de modelos de pronóstico de velocidad
aprovechamiento de energía eólica mediante la utilización de de viento en Arequipa, se han producido daños en la
modelos ARIMA aplicados a Tunja (Colombia), y se logró agricultura, materiales y pérdidas humanas, por esta razón es
identificar 4 lugares óptimos para ubicar aerogeneradores necesario hacer un pronóstico mediante un modelo ARIMA.
(Benia et al., 2022; Gómez, 2016). Por ello, la ciencia ha Con este fin, se cuenta con datos de velocidad de viento del
avanzado en los estudios realizados con el fin de conocer las Satélite ambiental operativo geoestacionario (GOES), con el
condiciones meteorológicas (temperatura, velocidad de viento, propósito de beneficiar al sector de agricultura (aumento de
presión atmosférica, humedad, radiación solar, etc.) a través cosechas), porque se encontrará una manera de anticipar las
del uso de técnicas como ARIMA, así como el coeficiente U fuertes velocidades de viento. Así también promoverá la
de Theil y el test de Diebold-Mariano (DM) con el que se elaboración de más modelos estocásticos para futuras
Prueba de Dicker-Fuller
Estimate Std. t-value Pr(>|t|)
a) 0,05233 -25,18 2 x10-16
Valores críticos tau1 para estadísticas de prueba:
1pct 5pct 10pct
Velocidad de viento (m/s)
100 Según la Figura 3a, se observa que no hay una tendencia a que
0 100 200 300
los puntos de las distribuciones sigan una línea recta. Por lo
Día tanto, la serie temporal no tiene un comportamiento normal.
b)
Figura 2. a) Diagrama de cajas b) Serie temporal de velocidad de viento del
2022 en R
-0.2
0.0
2
1
-0.3
0
0 5 10 15 20 25 0 5 10 15 20 25
Retrasos Retrasos
-1
móviles (MA), MA(q): MA (1), MA (2). Est. Std. 0,47059 0,1472 0,46845 0,4412
• Para la representación gráfica Parcial ACF, se observa z value 0,0856 -0,168 -1,4713 -0,514
2.5 % -0,88208 -0,313 -1,60738 -1,092
más de tres retardos significativos, se tomarán los dos 97.5 % 0,96261 0,2639 0,22870 0,6378
primeros retardos debido a que sobresalen más de los p-value = 0.8726
límites de confianza obteniendo los modelos KPSS Trend Truncation p-value
autorregresivos (AR), AR (p): AR (1), AR (2). lag parameter
0,078452 5 0,1
Para el criterio de elección de los retardos, se tomó en * MLE (log likelihood)
cuenta también los valores menores AIC, generados
por la función auto.arima de R. A partir de la Tabla 3, se pueden observar los modelos con
• Se obtiene el orden de las partes medias diferentes valores (p,d,q), con su error estimado estándar (Est.
autorregresivas, integradas y móviles del modelo, Std.), z value obtenidos, varianza de los residuos (σ2), intervalo
ARIMA (p d q): (1, 1, 1), (2, 1, 1), (1, 1, 2) (2, 1, 2). de confianza entre 2.5 % y 97.5 %, respecto a los coeficientes
de los modelos. Se observa que los modelos (1, 1, 2) y (2, 1,
2) tienen valores MLE: -1 638,33, -1 638,32, respectivamente.
partir de la Prueba KPSS el valor p=0,1 >0,05, entonces la Mejor modelo ARIMA
serie de tiempo diferenciada es estacionaria en tendencia. Comando Modelo ARIMA ar1 ma1 ma2
(1,1,2) -0,0265 -0,6315 -0,2912
s.e. 0,1894 0,1800 0,1643
A continuación, analizaremos el comportamiento residual del auto.arima
σ2 = 1227: log likelihood = -1641,18
modelo ARIMA (1, 1, 2). AIC=3 290,36 AICc*=3 290,48 BIC=3 305,56
Precisión del pronóstico en R
A partir de la Figura 5, se puede identificar que los ACF de los MAE MAPE MASE
residuos muestran autocorrelaciones no significativas (vea
28,140 18,713 0,849
Figura 5b) y tienen un comportamiento normal
*AICc: Criterio de información de Akaike de segundo orden
respectivamente (ver Figura 5c), y según la gráfica de densidad
3.4 Predicción mediante el uso de un modelo ARIMA
espectral se verifica que no hay un patrón de residuales (ver
Figura 5a), por lo tanto, se puede calcular el pronóstico.
La Figura 6 representa el modelo predictivo ARIMA (1,1,2).
1000
200
100
10 100
0.00 2.25
2.20
-0.05
2.15
-0.10
0 5 10 15 20 25 2.10
Retrasos 2.5 5.0 7.5 10.0
Enero
b)
12.5 b)
Figura 6. a) Predicción del modelo ARIMA (1, 1, 2), b) Comparación Log
10.0 [Velocidad de viento] desde 29 de noviembre hasta 8 de diciembre del 2022
en R
Frecuencias
Debido a la disponibilidad de datos meteorológicos online Badi, I., Pamučar, D., Stević, Ž., & Muhammad, L. J. (2023).
(POWER, 2023), se pudo proponer el modelo ARIMA (1,1,2) Wind farm site selection using BWM-AHP-MARCOS
para Socabaya. Esto contribuye al aumento de investigación method: A case study of Libya. Scientific African, 19,
de predicción como Bórawski et al. (2020) y Salman & e01511. https://doi.org/10.1016/J.SCIAF.2022.E01511
Kanigoro (2021). Benia, T. K., Goswami, S., & Banerjee, A. (2022). Application
of SPSS for Forecasting of Renewable Energy as Future
Los fuertes vientos de Socabaya produjeron accidentes de Energy in India. In Advanced Machine Intelligence and
tránsito, desbordamiento de tierras en las chacras por ello el Signal Processing (pp. 379-391). Singapore: Springer
trabajo se centra en un modelo de predicción ARIMA de buena Nature Singapore.https://doi.org/10.1007/978-981-19-
precisión MASE de 0.849, a diferencia de China que utilizaron 0840-8_28
un modelo RARIMA (Hui et al., 2015; Pilco & Acurio, 2019). Bórawski, P., Bełdycka-Bórawska, A., Jankowski, K. J.,
Dubis, B., & Dunn, J. W. (2020). Development of wind
La predicción de velocidad de viento para Socabaya a corto energy market in the European Union. Renewable
plazo desde el 29 de noviembre hasta 8 de diciembre del 2022 Energy, 161, 691-700.
a diferencia de la predicción, 6-8 horas en el Báltico hecha por https://doi.org/10.1016/j.renene.2020.07.081
Grigonytė & Butkevičiūtė (2016). Cadenas, E., & Rivera, W. (2010). Wind speed forecasting in
three different regions of Mexico, using a hybrid
En la zona rural de Socabaya, se pronosticó la velocidad del ARIMA–ANN model. Renewable Energy, 35(12),
viento mediante un modelo ARIMA (1,1,2) a diferencia de las 2732-2738.
islas de Baja California, en las cuales, el modelo ARIMA no https://doi.org/10.1016/j.renene.2010.04.022
pudo pronosticar recurriendo a una ANN (Cadenas & Rivera, Cadenas, E., Rivera, W., Campos-Amezcua, R., & Heard, C.
2010). (2016). Wind speed prediction using a univariate
ARIMA model and a multivariate NARX model.
Una limitación del trabajo fue el tamaño de la muestra que Energies, 9(2), 109. https://doi.org/10.3390/en9020109
delimitó la aplicación de los métodos paramétricos y Dönmez, K., Çetek, C., & Kaya, O. (2022). Air traffic
determinación del mejor modelo ARIMA de la serie temporal. management in parallel-point merge systems under wind
uncertainties. Journal of Air Transport Management,
4. CONCLUSIONES 104, 102268.
https://doi.org/10.1016/j.jairtraman.2022.102268
Con base en los resultados, se obtuvo el modelo de predicción Dumitru, C. D., & Gligor, A. (2019). Wind energy forecasting:
óptima ARIMA (1, 1, 2) de la serie temporal de velocidad de a comparative study between a stochastic model
viento 2022 en R para el Distrito de Socabaya-Arequipa, la (ARIMA) and a model based on neural network
cual cumple con las condiciones de estacionariedad en (FFANN). Procedia manufacturing, 32, 410-417.
varianza, normalidad, invertibilidad, independencia de https://doi.org/10.1016/j.promfg.2019.02.234
residuos y posee menor error. Gómez Arévalo, G. A. (2016). Desarrollo de una metodología
para obtener un atlas de viento de una microzona con
El modelo ARIMA (1, 1, 2) es viable por su precisión MASE datos disponibles de corto plazo en Tunja, Colombia.
mayor al 0,8 entre datos de velocidad de viento de NASA y los https://hdl.handle.net/10901/9816
datos pronosticados del modelo. Por lo tanto, el modelo sirve Grigonytė, E., & Butkevičiūtė, E. (2016). Short-term wind
para predecir la velocidad del viento y evitar desastres en la speed forecasting using ARIMA model. Energetika,
agricultura, accidentes de transportabilidad, así como para 62(1-2). https://doi.org/10.6001/energetika.v62i1-
mejorar la producción de productos agrícolas, salud, desarrollo 2.3313
económico y seguridad pública del distrito y será Hodge, B.S., Zeiler, A., Brooks, D., Blau, G.E., Pekny, J.F., &
preponderante para estudios del cambio climático. Reklatis, G. (2011). Improved Wind Power Forecasting
with ARIMA Models. Computer-aided chemical
Los procedimientos utilizados en la metodología pueden ser engineering, 29, 1789-1793.
útiles para determinar modelos de predicción de vientos en https://doi.org/10.1016/B978-0-444-54298-4.50136-7
zonas con altos niveles de velocidades de viento. Hui, L., Tian, H. Q., & Li, Y. F. (2015). An EMD-recursive
ARIMA method to predict wind speed. J Wind End Ind
No se consideraron los modelos ARIMA (1, 1, 1), (2, 1, 1), (2, Aerodyn, 141, 27e38.
1, 2) para la predicción de velocidad de viento porque sus https://doi.org/10.1016/j.jweia.2015.02.004
residuos eran correlacionados, y no se logró obtener el valor Jastrzebska, A., Hernández, A. M., Nápoles, G., Salgueiro, Y.,
más pequeño AIC, además que su intervalo de confianza no es & Vanhoof, K. (2022). Measuring wind turbine health
muy alto. using fuzzy-concept-based drifting models. Renewable
Energy, 190, 730-740.
Se ha logrado obtener un modelo de predicción de velocidad https://doi.org/10.1016/j.renene.2022.03.116
de viento que podría estar influenciado por la variabilidad Shewhart, W. A., Wilks, S. S., Balding, D. J., Cressie, N. A.
climática (presión, temperatura, radiación solar, etc.) C., Fitzmaurice, G. M., Givens, G. H., Goldstein, H.,
Molenberghs, G., Scott, D. W., Smith, A. F. M., Tsay, estudios de Postgrado de Maestría en Ciencias, Física
R. S., Hunter, E. J. S., Johnstone, I. M., Kadane, J. B., & aplicada a la Medicina y Radioprotección en la UNSA
Teugels, J. L. (2019). WILEY SERIES IN PROBABILITY (2022), docente seleccionado por PMESUT del programa
AND STATISTICS. de capacitación en investigación científica (2022),
http://www.wiley.com/go/permissions docente seleccionado para la presentación de POSTER de
Li, X., Cao, Y., Yu, X., Xu, Y., Yang, Y., Liu, S., ... & Wang, un trabajo de investigación en la UNAMAD (2022),
Z. L. (2022). Breeze-driven triboelectric nanogenerator Diplomado de Especialización en Docencia de Educación
for wind energy harvesting and application in smart Superior (2021).
agriculture. Applied Energy, 306, 117977.
https://doi.org/10.1016/j.apenergy.2021.117977 Renzon, Huamani, Egresado de
Pilco Sánchez, V. K., & Acurio Martinez, W. D. (2019). la Escuela Profesional de Física
Técnicas estadísticas para la modelación y predicción de la UNSA (2022), realiza
de la temperatura y velocidad de viento en la provincia trabajos de investigación en el
de Chimborazo (Bachelor's thesis, Escuela Superior Área de Energías Renovables,
Politécnica de Chimborazo). seleccionado para la
http://dspace.espoch.edu.ec/handle/123456789/10955 presentación de POSTER de un
POWER Data Access Viewer. (2023). Recuperado de trabajo de investigación en la
https://power.larc.nasa.gov/data-access-viewer/ UNAMAD (2022).
(consultado el 20 de diciembre de 2022)
Pulido Calvo, I., Roldán Cañas, J., López Luque, R., &
Gutiérrez Estrada, J. C. (2002). Técnicas de predicción
a corto plazo de la demanda de agua. Aplicación al uso
agrícola. Ingeniería del agua, 9(3), 319-331.
https://doi.org/10.4995/ia.2002.2622
Salman, A. G., & Kanigoro, B. (2021). Visibility forecasting
using autoregressive integrated moving average
(ARIMA) models. Procedia Computer Science, 179,
252-259. https://doi.org/10.1016/j.procs.2021.01.004
Sim, S. K., Maass, P., & Lind, P. G. (2018). Wind speed
modeling by nested ARIMA processes. Energies, 12(1),
69. https://doi.org/10.3390/en12010069
Stigter, C. J., Mohammed, A. E., Nasr Al-Amin, N. K.,
Onyewotu, L. O. Z., Oteng’i, S. B. B., & Kainkwa, R.
M. R. (2002). Agroforestry solutions to some African
wind problems. Journal of Wind Engineering and
Industrial Aerodynamics, 90(10), 1101–1114.
https://doi.org/10.1016/S0167-6105(02)00224-6
Uriel, E. (1985). Análisis de series temporales modelos Arima
(No. 517.2/U77a).
Zahedi, R., Ghorbani, M., Daneshgar, S., Gitifar, S., &
Qezelbigloo, S. (2022). Potential measurement of Iran's
western regional wind energy using GIS. Journal of
Cleaner Production, 330, 129883.
https://doi.org/10.1016/j.jclepro.2021.129883
Ziemba, P., Wa̧tróbski, J., Zioło, M., & Karczmarczyk, A.
(2017). Using the PROSA method in offshore wind farm
location problems. Energies, 10(11).
https://doi.org/10.3390/en10111755
BIOGRAFÍAS
1
Universidad Central del Ecuador, Facultad de Ingeniería y Ciencias Aplicadas, Quito, Ecuador
Resumen: Para analizar y comprender el comportamiento del flujo en estructuras hidráulicas, se construyen modelos
a escala, sin embargo, en el caso del disipador de energía tipo vórtice, la naturaleza helicoidal del flujo dificulta su
análisis y la medición de datos, por lo tanto, se ha creado un modelo numérico a escala real que permitiera un análisis
más eficiente del flujo en este tipo de disipador. Se aplicó la metodología de la dinámica de fluidos computacional
mediante el software adecuado. Primero fue creada la geometría del disipador, luego el mallado, posteriormente se
establecieron las condiciones de borde y parámetros físicos, con lo cual se resolvió el modelo obteniendo resultados
que fueron validados con datos experimentales del modelo físico, para finalmente analizar el prototipo. Los resultados
obtenidos no solo demostraron la aplicabilidad de la dinámica de fluidos computacional (CFD) en el diseño hidráulico
y su calibración, sino que también proporcionaron una distribución espacial de las características hidráulicas a lo
largo del prototipo, aspecto esencial en el diseño de este tipo de estructuras.
*[email protected]
Recibido: 15/03/2024
Aceptado:21/06/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.02
CC 4.0
Tomando en cuenta que en un modelo numérico se pueden Figura 1. Diagrama de flujo. Metodología CFD
simular diferentes estructuras hidráulicas y del material que se
requiriese, es imperativo considerar valores de rugosidad 2.1. Creación de la geometría
acorde al material de la estructura. En esta investigación, se
han desarrollado simulaciones con material acrílico para el La simulación numérica se ha llevado a cabo partiendo de la
modelo a escala, así como hormigón para el prototipo y tras creación de la geometría del disipador tipo vórtice mediante
realizar un análisis de sensibilidad se ha podido concluir que Solidworks Student 2023 (SolidWorks, 2023), para esto ha
este parámetro puede obviarse tras una correcta configuración sido necesario remitirse a los planos realizados por Chiluisa y
de datos iniciales y validar los resultados en función de datos Guanoluisa (2015).
experimentales, con el fin de optimizar el modelo y tiempos de
simulación. 2.2. Mallado de la geometría
El análisis de los resultados obtenidos en la simulación La discretización del dominio computacional del disipador
numérica del prototipo permitirá una comprensión más tipo vórtice se ha generado a través de un mallado mixto
profunda del comportamiento hidráulico del disipador de combinando elementos hexaédricos (malla estructurada) y
energía en un flujo estacionario y su aplicación en el diseño tetraédricos (malla no estructurada) según se recomienda en
hidráulico. Esto debido a la versatilidad que ofrece CFD con Yang et al. (2021). Las zonas que poseen cambios de dirección
la generación de diversos escenarios, por lo que este estudio y secciones variadas en donde la malla no estructurada se
servirá como una base para futuros diseños de disipadores de encuentra en la cámara espiral, conducto vertical y cámara de
energía (Ogunsesan et al., 2021). disipación. Además, se ha realizado el refinamiento cuatro
veces más pequeño de la malla en las aristas y cambios de
2. METODOLOGÍA sección del disipador.
El diagrama de flujo de la Figura 1 describe el proceso del En relación con la discretización del dominio computacional,
trabajo. también se llevó a cabo un análisis de independencia mediante
la realización de cinco escenarios con diferentes tamaños de
celdas. El escenario que cumple con la calidad de mallado,
parametrización de calidad de ANSYS CFX (Ansys, 2023) y
con la calibración de resultados tiene las siguientes – ε, que tardó 4h 16min de simulación, debido a su menor gasto
características: tamaño de celda de 0.01m, un total de 307 766 computacional en comparación con el modelo RNG k – ε, que
nodos y 4h16minutos de simulación. En relación con los tardó 5h 10min de simulación.
parámetros del mallado óptimo, se han obtenido los siguientes
valores: 0.25 para la oblicuidad, 0.83 para la ortogonalidad, Los datos climáticos usados en el modelo han sido: presión de
2.03 para el aspecto del elemento y 0.94 para la calidad del referencia de 0.71 atm, gravedad de 9.81 m/s2 y temperatura
elemento. del aire a 20°C. A pesar de que la simulación numérica podría
realizarse para cualquier caudal, esta investigación se centra
Se ha considerado la limitante de nodos que se pueden utilizar en la simulación de cinco caudales específicos, como lo indica
en el programa en su versión estudiantil (512 000 nodos). la Tabla 1, utilizados con el propósito de realizar
Además, las características del computador utilizado fueron: comparaciones con los resultados obtenidos del modelo físico
16 GB de memoria RAM, SSD de 512 GB y un procesador y validar la simulación. Esta comparación se lleva a cabo en
Core i7 de 12ª generación. cuatro secciones, las cuales se muestran en la Figura 2.
diferencia entre el resultado del modelo numérico y el valor por lo que existen errores porcentuales de más del 100%. Esta
experimental dividido por el valor experimental, son discrepancia se atribuye a la turbulencia en el disipador, un
generalmente aceptables, siendo menores al 10% (Yang et al., fenómeno complejo que puede ser desafiante de modelar con
2021). precisión.
Los errores que exceden el margen establecido se encuentran Con respecto a la sección 4-4, tanto para calados como
en las secciones 1 y 2 de velocidades, para un caudal de 40.78 velocidades, errores mayores al 10% son atribuibles a su
l/s. Este error refleja la incertidumbre inherente en la proximidad a la cámara de disipación. En esta zona, se produce
recopilación de datos físicos, pues para caudales anteriores, el un flujo mixto de agua y aire, fenómeno que el software puede
error es inferior al 10%. tener dificultades para simular con precisión, resultando en
estas variaciones. A pesar de estas anomalías, los resultados
Es crucial destacar que en la sección 3-3, correspondiente a la presentados en la Tabla 1 proporcionan una base sólida para
cámara de disipación, el comportamiento del flujo es variable, avanzar con la simulación y el análisis del prototipo.
3. DISCUSIÓN DE RESULTADOS
5
Qp =Qm ∙ λ2 (1)
Dado que la parametrización del modelo numérico ha sido
validada previamente con los resultados del modelo físico, se
dispone de la configuración adecuada para simular el disipador Siendo lambda (λ) la relación entre las dimensiones
de energía tipo vórtice en un escenario distinto. En este caso, longitudinales del prototipo y el modelo físico.
se trataría de una geometría de mayor tamaño, equivalente al
prototipo. Este enfoque permite la proyección de los resultados Es importante destacar que los resultados obtenidos en este
obtenidos en la simulación del modelo a escala a un sistema de estudio están sujetos a cierto grado de incertidumbre,
tamaño real, proporcionando una herramienta valiosa para el atribuible a diversas fuentes. En primer lugar, el método de
diseño y optimización de disipadores de energía tipo vórtice Volumen de Fluido (VOF) implementado en el software
en aplicaciones prácticas. ANSYS que presenta limitaciones en la simulación precisa del
flujo bifásico agua-aire (Yang et al., 2021). Específicamente,
Para la transposición de las condiciones iniciales de un modelo puede surgir dificultades al resolver con exactitud las
físico a un prototipo se hace uso de la teoría de la similitud, interfaces entre ambos fluidos, un desafío que se intensifica en
propuesta por Kline, que menciona: “Si dos sistemas obedecen flujos con alta turbulencia o en presencia de fenómenos
al mismo grupo de ecuaciones y condiciones gobernantes, y si complejos como la atomización.
los valores de todos los parámetros y las condiciones se hacen
idénticas, los sistemas deben exhibir comportamientos Además, las condiciones ambientales pueden introducir
similares con tal de que exista una solución única para el grupo incertidumbre adicional. Un conocimiento insuficiente de los
de ecuaciones y condiciones” (Sotelo, 1974). efectos de estas condiciones sobre la medición, o una medición
imperfecta de las mismas, puede generar variaciones en los
Es por ello que esta simulación se la realizará con un caudal resultados.
de 25.46 m3/s y 35.54 m3/s, que por similitud hidráulica está
en correspondencia con el caudal de 29.22 l/s y 40.78 l/s, ya Por último, los datos de entrada iniciales, obtenidos a través de
que la relación entre el caudal del modelo y prototipo está mediciones en el disipador, están sujetos a un margen de error
establecida con la siguiente fórmula: inherente. Estos errores pueden propagarse a través del modelo
numérico y afectar la precisión de los resultados.
4,0
A pesar de estas fuentes de incertidumbre, es relevante señalar
que no han impactado de manera significativa los resultados 3,8
de este estudio, ya que la validación de los resultados ha
demostrado ser óptima. Sin embargo, es crucial tener en cuenta 3,6
estas fuentes de incertidumbre al interpretar los resultados y al 3,4
Calado (m)
aplicar el modelo computacional a prototipos a escala real.
3
3,2 Q = 25.46 m /s
3
3.1 Análisis de sensibilidad Q = 35.54 m /s
3,0
-4
Cota (m)
-6
-8
-10
7 8 9 10 11 12 13 14
Velocidad (m/s)
Figura 5. Velocidad en el tubo vertical
-2
-4
Cota (m)
-6
Figura 7. Longitud máxima del tubo vertical. Caudal: 25.46 m3/s
-8 3
Q = 25.46 m /s Una de las limitantes que presenta el disipador tipo vórtice en
3
Q = 35.54 m /s su funcionamiento, es la velocidad máxima con la que puede
-10 fluir el agua en el tubo vertical, tal como se puede apreciar en
la Figura 8. A medida que el tubo vertical se extiende, se
-12 observa que existe una velocidad límite que aumenta en
15 20 25 30 35 40 45 función del caudal que se está desarrollando. Esta velocidad
Ángulo (°)
límite es atribuida a la presencia del núcleo de aire en el
Figura 6. Ángulo de inclinación del flujo en el tubo vertical interior del tubo vertical, el cual genera resistencia al
movimiento del flujo (Rodal et al., 2011).
Una observación acertada de este comportamiento es que el
vector velocidad coincide con la dirección de las líneas de Para el diseño del tubo vertical es necesario considerar una
flujo, como es de esperarse, el cual se puede medir en ANSYS correcta aireación, así como su funcionalidad, para esto se ha
obteniéndose valores importantes para realizar cálculos hecho uso del valor lambda (λ), mismo que representa la
analíticos. relación entre el área del núcleo de aire y el área de la sección
del tubo vertical. Considerando que se recomienda un valor
Con respecto a la longitud máxima que el tubo vertical puede por encima de 0.25, se ha examinado este número a través de
poseer, es importante mencionar el trabajo de Hager (1985), la simulación numérica y se ha encontrado que la variación de
quién presentó una fórmula para estimar esta longitud, en la lambda a lo largo del tubo vertical alcanza un valor mínimo de
cual la velocidad tangencial disminuye su magnitud, dejando 0.41, esto garantiza el correcto funcionamiento del prototipo
de tener un impacto significativo en la resultante de velocidad (Chan et al., 2018).
lo cual provoca que el flujo adquiera un comportamiento
vertical. La expresión matemática que describe esta relación es Finalmente, se observaron presiones negativas en la unión del
la siguiente: tubo vertical con la cámara de disipación. Para evaluar estas
6 4
presiones es necesario recurrir al índice de cavitación el cual
3K5 Q 5 responde a la siguiente fórmula.
Zd = (π∙D) (2)
g
P -P
IC= 0.5abs∙ v2 ∙vρ (3)
Donde, K es el inverso del coeficiente de Manning, g es la
gravedad, Q es el caudal y D es el diámetro del tubo vertical.
Donde, Pabs es la presión absoluta, PV es la presión de vapor,
Aplicando la fórmula (1) se han obtenido longitudes máximas v es la velocidad y ρ es la densidad del agua.
en función de los caudales simulados, con lo cual se ha
modificado la geometría del tubo vertical y se consiguió
visualizar el comportamiento de las líneas de flujo como se
limitaciones. En particular, los modelos RANS proporcionan energética adecuada, dado que la eficiencia de la disipación
una descripción promediada en el tiempo del flujo turbulento, superó el 80 %.
pudiéndose no capturar adecuadamente los detalles finos de la
turbulencia, como se observó en la cámara de disipación, con Al aumentar la altura del tubo vertical, se pudo observar que
errores mayores al 10%. las velocidades adquieren una tendencia a crecer con forme el
flujo circula hacia la cámara de disipación. Sin embargo, esta
Para obtener un análisis más detallado y específico del flujo velocidad llega a un límite donde permanece prácticamente
dentro de la cámara, sería beneficioso utilizar un modelo LES constante debido a la interacción que existe entre el núcleo de
(Large eddy simulation), el cual se centra en resolver las aire y las paredes del disipador.
estructuras turbulentas más grandes en el flujo, mientras que
los efectos de las estructuras más pequeñas se modelan Es crucial garantizar que la longitud del conducto vertical se
(Capote et al., 2008). Sin embargo, estos modelos requieren encuentre dentro del intervalo de la altura máxima. Esto se
una resolución de malla mucho más fina y, por lo tanto, debe a que longitudes superiores resultarían en la pérdida del
consumen significativamente más recursos computacionales. comportamiento helicoidal, lo que generaría un movimiento en
Aunque el tiempo de simulación con LES puede ser de caída libre. Como consecuencia, se perdería la función original
semanas, este enfoque no se descarta si se dispone de para la cual se diseñó la estructura.
suficientes recursos y tiempo. De hecho, el uso de LES puede
ser muy beneficioso para obtener una comprensión más REFERENCIAS
profunda de los fenómenos de flujo complejos.
Ansys. (2023). ANSYS CFX [versión estudiantil para PC]
Es importante tener en cuenta estas consideraciones al (2023 R1).
seleccionar el modelo de turbulencia más apropiado para una Capote, J. A., Alvear, D., Abreu, O. V, Lázaro, M., & Espina,
determinada simulación. P. (2008). Influencia del modelo de turbulencia y del
refinamiento de la discretización espacial en la exactitud
4. CONCLUSIONES de las simulaciones computacionales de incendios.
Revista Internacional de Métodos Numéricos Para
El modelo de turbulencia k-ε logró una mejor simulación del Cálculo y Diseño En Ingeniería, 24(3).
comportamiento del flujo en el disipador de energía. Los datos Chan, S. N., Qiao, Q. S., & Lee, J. H. W. (2018). On the three-
experimentales concordaron de manera adecuada con los dimensional flow of a stable tangential vortex intake.
resultados del modelo numérico. Además, este modelo de Journal of Hydro-Environment Research, 21(2018), 29-
turbulencia demostró ser más eficiente en términos de tiempo 42. https://doi.org/10.1016/j.jher.2018.07.001
de resolución en comparación con el modelo RNG k-ε. Estos Chan, S. N., Qiao, Q. S., & Lee, J. H. W. (2019). Flow Features
hallazgos respaldan la posibilidad de emplear el modelo k-ε en of an Unstable Tangential Vortex Intake. Green Energy
futuros proyectos que involucren flujos con naturaleza and Technology. https://doi.org/10.1007/978-3-319-
helicoidal, así amplían las posibilidades de aplicación en 99867-1_55
diversos escenarios. Chiluisa C., & Guanoluisa C. (2015). Verificación del diseño
hidráulico de un disipador de energía tipo vórtice
El análisis de sensibilidad realizado en este estudio reveló que mediante modelo físico construido en laboratorio.
la rugosidad de las paredes del disipador de energía tiene un Universidad Central del Ecuador.
impacto mínimo en los resultados de la simulación. Esta Crispino, G., Contestabile, P., Vicinanza, D., & Gisonni, C.
observación crucial sugiere que la exclusión de este parámetro (2021). Energy Head Dissipation and Flow Pressures in
podría simplificar significativamente el modelo. Tal Vortex Drop Shafts. Water, 13(2), 165.
simplificación no solo acelera el proceso de simulación, https://doi.org/10.3390/w13020165
especialmente en proyectos de gran escala, sino que también Del Giudice, G., Gisonni, C., & Rasulo, G. (2009). Vortex
reduce notablemente el tiempo requerido para obtener Drop Shaft for Supercritical Flow. In Advances in Water
resultados confiables. Lo más importante es que esta Resources and Hydraulic Engineering (pp. 1515–1520).
optimización se logra sin sacrificar la precisión y la fiabilidad Springer Berlin Heidelberg.
del modelo. Por lo tanto, la exclusión del parámetro de https://doi.org/10.1007/978-3-540-89465-0_262
rugosidad podría ser una estrategia valiosa para mejorar la Fernández J. (2011). Técnicas numéricas en ingeniería de
eficiencia de las simulaciones de disipadores de energía tipo fluidos (Editorial Reverté, Ed.; 1st ed.).
vórtice, lo que tiene implicaciones significativas para la Hager, W. H. (1985). Head‐Discharge Relation for Vortex
modelación y el diseño en la ingeniería hidráulica. Shaft. Journal of Hydraulic Engineering, 111(6).
https://doi.org/10.1061/(asce)0733-
En el proceso de validación del modelo numérico se pudo 9429(1985)111:6(1015)
inferir que la sección 3 incluida en la cámara de disipación, fue Mrope, H. A., Chande Jande, Y. A., & Kivevele, T. T. (2021).
la única zona donde no se logró una concordancia con los A Review on Computational Fluid Dynamics
resultados experimentales, ya que su error excedió el 10 %. Applications in the Design and Optimization of
Esto trasladado al prototipo solo resalta que los valores en esta Crossflow Hydro Turbines. Journal of Renewable
región pueden carecer de fiabilidad debido a su Energy, 2021, 1–13.
comportamiento impredecible. Pero en un enfoque más https://doi.org/10.1155/2021/5570848
general, estas zonas continúan demostrando una disipación
Ogunsesan, O. A., Hossain, M., & Droubi, M. G. (2021). Oscar, Ulquiango, Ingeniero
Computational fluid dynamics modelling of multiphase civil graduado de la Universidad
flows in double elbow geometries. Proceedings of the Central del Ecuador, sobresale
Institution of Mechanical Engineers, Part E: Journal of por su experiencia en
Process Mechanical Engineering, 235(6), 1835–1846. presupuestos, ejecución de
https://doi.org/10.1177/09544089211021744 obras civiles y la aplicación de
Rodal E., Echávez G., & Ruiz G. (2011). Estudio experimental la CFD en proyectos
en modelo hidráulico del pozo de alta caída. Instituto hidráulicos. Apasionado por el
Argentino de Recursos Hídricos. crecimiento personal y
Rodríguez F. (2016). Modelación numérica del flujo en un profesional, busca
descargador a vórtice aplicando mecánica de fluidos constantemente aprender y
computacional (CFD). Escuela Politécnica Nacional. enriquecer sus conocimientos.
Ruiz Solorio, G. (2014). MODELOS DE TURBULENCIA
EN MECANICA DE FLUIDOS COMPUTACIONAL.
XXII CONGRESO NACIONAL DE HIDRAULICA. Paulina, Lima, Profesora en el área
Sentyabov, A. V., Timoshevskiy, M. V., Pervunin, K. S., de hidráulica de la carrera de
Gavrilov, A. A., Markovich, D. M., & Dekterev, A. A. Ingeniería Civil en la Universidad
(2016). Numerical and experimental investigation of Central del Ecuador desde 2014,
cavitation flow around naca0015 hydrofoil. Bulletin of ostenta un Doctorado en Ingeniería
the Tomsk Polytechnic University, Geo Assets Civil de la Universidad de Nuevo
Engineering, 327(8). México, Estados Unidos, obtenido
SolidWorks. (2023). SolidWorks [versión estudiantil para en 2019. Tras graduarse en
PC] (2023). Ingeniería Civil en la UCE en 2004,
Sotelo, G. (1974). Hidráulica General (Limusa, Ed.). Lima se especializó en Recursos
Usman, M., Shahid, S., Ali, S., & Ullah, M. K. (2023). Hídricos y Ciencias del Agua con
Numerical simulations of turbulent and flow un Magíster de la Escuela Politécnica Nacional en 2007 y
characteristics of complex river reach in Pakistan. un Máster Complementario en Gestión de Riesgos
Environmental Engineering Research, 28(1). Naturales en Bélgica en 2009. Destaca por su trabajo en
https://doi.org/10.4491/eer.2021.369 estudios de ingeniería hidráulica desde 2004, con
Yang, Z., Yin, J., Lu, Y., Liu, Z., Yang, H., & Xu, G. (2021). contribuciones significativas en desarrollo de proyectos
Three-dimensional flow of a vortex drop shaft spillway hidráulicos relacionados con la Ingeniería Civil en el
with an elliptical tangential inlet. Water (Switzerland), Ecuador.
13(4). https://doi.org/10.3390/w13040504
Zhang, W., Wang, J., Zhou, C., Dong, Z., & Zhou, Z. (2018).
Numerical simulation of hydraulic characteristics in a
vortex drop shaft. Water (Switzerland), 10(10).
https://doi.org/10.3390/w10101393
Zhao, C.-H., Zhu, D. Z., Sun, S.-K., & Liu, Z.-P. (2006).
Experimental Study of Flow in a Vortex Drop Shaft.
Journal of Hydraulic Engineering, 132(1).
https://doi.org/10.1061/(asce)0733-
9429(2006)132:1(61)
BIOGRAFÍAS
1 Universidad Tecnológica del Norte de Aguascalientes, Departamento de Ingeniería Mecatrónica, Rincón de Romos, México
2 Universidadde Guadalajara, Departamento de Ciencias Exactas y Tecnología, Lagos de Moreno, México
3 Universidad Autónoma de la Ciudad de México, Colegio de Ciencia y Tecnología, Ciudad de México, México
Abstract: This paper presents a closed-loop controller for speed regulation of a three-phase induction motor by means of
a Modbus network, which is compound only by industrial-purpose hardware, namely two programmable logic controllers,
a human-machine interface and a variable-frequency drive. This represents a step beyond the state of the art because the
results in similar works reported in the literature were obtained using personal computers, academic-purpose hardware
or a combination of both academic-industrial technologies. The speed regulation was achieved through a proportional-
integral controller executed in a programmable logic controller, whose feedback was supplied by another similar device
completely dedicated to read and decode an industrial optical encoder; furthermore, the control signals were commanded
to the alternating current motor using a variable-frequency drive. The proposed control network was tested through diverse
regulation experiments, under different conditions, in order to assess its effectiveness and robustness. The obtained results
were successful in terms of accuracy of the speed set-point tracking in comparison to the most similar works of the
state of the art, namely zero steady-state error was achieved, 12 seconds earlier for 75 % larger references and 22 %
less RMSE than other works, which are the best results reported so far. Such successful achievements encourage the
later implementation of more advanced control techniques on three-phase-motor based industrial machinery using only
industrial-purpose hardware.
Keywords: Three-phase motor control, PLC based control, Modbus networks, AC motor speed regulation
Palabras claves: Control de motores trifásicos, Control basado en PLC, Redes Modbus, Regulación de velocidad de
motores de CA
1. I NTRODUCTION kare and Jagtap, 2017). Another relevant technology for indus-
trial processes is the programmable logic controller (PLC), whose
reliability, operation simplicity, electrical and mechanical robust-
Several machines that compound diverse manufacturing proces- ness and high speed turn it into the preferred device in industry for
ses are actuated by three-phase alternating current (AC) induction the execution of automation and control algorithms (Dukare and
motors because of its simplicity, robustness and cost, therefore, Jagtap, 2017). Additionally, the possibility of developing PLC net-
controlling its operation is a highly relevant task in industry (Du-
*[email protected]
Recibido: 23/11/2023
Aceptado: 11/03/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.03
CC BY 4.0 Revista Politécnica, Agosto - Octubre 2024, Vol. 54, No. 1
Martín Rodríguez; Mario Baes; Ricardo Jara; Yadira López; Santos Orozco
26
works that might include other industrial-purpose hardware, such environments is still a drawback for such approaches since the ro-
as human-machine interfaces (HMI), increases the aforesaid pre- bustness of the system might be compromised. Recently, a very
ference of PLC and stand out their versatility. For such a reason, interesting fully industrial setup for induction motors control was
several manufacturing companies in Mexico are reluctant to the reported by Vadi et al. (2022), these authors use a PLC Profibus
use of computers, embedded systems or laboratory-purpose hard- network to communicate the PLC with the VDF, but the closed-
ware and software to automate or control their machinery, even loop control action is still missing.
though these devices might execute cutting-edge automation and Considering the information found in the literature, the problem
control algorithms that could be beneficial to optimize or improve of implementing closed-loop control systems using just industrial-
industrial processes. purpose hardware is not fully overcome. In this context, the contri-
Considering the aforementioned problem, an appropriate alternati- bution of this paper is twofold: in one hand, up to the best authors’
ve for the industry could be the implementation of automation and knowledge, this is the first time that a PI controller is analyzed,
control algorithms using just industrial-purpose hardware. In this clearly explained and fully programmed within a PLC using an
context, and recalling the idea of the wide usage of three-phase industrial encoder as feedback. This represents a step towards the
AC motors, several researchers have reported PLC-based control implementation of more complex control systems on industrial de-
systems for such electrical machines (Verma and Potdar, 2017) in vices. On the other hand, the control system uses a Modbus net-
the literature. One of the most remarkable results was presented work in order to communicate the feedback-analyzing PLC with
by Saad and Arrofiq (2012), in which the authors programmed a the control PLC, the VFD and also an HMI that is useful to com-
fuzzy controller for an AC motor within the PLC, but they used a mand the target speed and also to monitor the measured velocity.
laboratory-type sensor, namely a DC motor adapted to work as a In order to assess the correct operation and the robustness of the
generator. Similarly, Howimanporn et al. (2016) reported an opti- full setup, the following experiments were put into effect: firstly,
mized PID controller for a conveyor belt actuated by an AC motor, the controller was commanded with large target velocities to pro-
nonetheless, the authors also used laboratory-purpose hardware to ve its stability in the whole error space; later on, the importance
send the control signals to the motor. Some examples of the usage of the integral action is evaluated and highlighted through the dis-
of laboratory-purpose software can be found in the works of Al- turbance rejection capability assessment. The successful results
Manfi (2019) and Awdaa et al. (2020); in both cases the closed- demonstrated that the proposed control architecture is suitable to
loop control for the motor is implemented by means of high-cost effectively control three-phase AC motors by means of only indus-
software executed on a personal computer (PC) in addition to the trial hardware, eliminating the necessity of high-cost additional
PLC. software or laboratory-purpose hardware that is more likely to fail
On the other hand, some works that only use industrial hardwa- in industrial environments. Furthermore, the developed control al-
re have been reported as well, for instance, the speed control of gorithm yielded better results than those reported in the literature
a three-phase motor for a centrifugal machine of a sugar manu- in terms of time of response and quantitative performance indices.
facturing process was reported by Pramudijanto et al. (2015). The Finally, the implementation of a Modbus network stands out the
authors implemented a predictive PID controller on a PLC, ho- versatility of the industrial hardware and opens the possibility of
wever, they do not specify the technology used for feedback, the involving more devices and adding functionalities depending on
execution of the control algorithm is not explained and their results the application requirements.
are not clearly presented in plots, but only a 10.79 root median The organization of this paper is as follows: firstly, the used
squared error (RMSE) is presented. In addition, Mahesh and Ra- methods are described in Section 2, namely the industrial hard-
machandra (2017) presented a PLC-based P controller using just ware network is described in addition to the mathematical model
industrial-purpose hardware, obtaining successful results with an of the plant as well as the control system design and its implemen-
average of 2 % steady-state error. Nevertheless, the authors do not tation. Then, the Section 3 addresses the experiments carried out
show plots of the behavior of the speed and they do not present and the obtained results are discussed and contrasted against the
theoretical foundations. Another PLC based closed-loop contro- state of the art. Finally, the Section 4 concludes the paper.
ller without laboratory-purpose hardware is reported in the work
2. M ETHODOLOGY
of Demir et al. (2019), achieving a zero steady-state error in 40
seconds for their largest step reference from 0 to 1000 [RPM]. Below, the full system setup used for this research is described,
Despite these successful results, it is worth mentioning that the as well as the development and implementation of the proposed
authors did not program the control algorithm but they only imple- control technique using the industrial hardware network.
mented and configured the proportional-integral-derivative (PID)
controller functional block of the PLC programming environment.
2.1 System Description
This represents a limitation if it is pretended to implement more
advanced control techniques. The implemented experimental system is displayed in Figure 1. It
Besides, some researchers have tried to exploit the versatility of is composed of two control and power panels, respectively, as well
the PLC networks to control the speed of AC motors. For example, as an independent station dedicated for the speed data acquisition
the use of open-platform communications (OPC) servers using from the three-phase AC motor. The first panel is integrated by a
Ethernet and other field buses to connect the PLC with variable- Micrologix 1500 PLC (hereinafter PLC1) and a HMI Panel View
frequency drives (VFD) were presented by Velagic et al. (2011) Plus 600 touch screen. Furthermore, the 1761-NET-ENI and 1761-
and Al-Manfi (2019), in order to implement PID controllers. Ho- NET-AIC modules have been adapted to the PLC1 to communi-
wever, the use of a non-industrial-purpose PC within industrial cate it with other devices through the Ethernet and RS845 com-
is not applicable. Besides, the gain tuning procedure was conduc- panel with a range from 0 to 1,800 [RPM] and later on is sent to
ted by analyzing the pole of the closed-loop transfer function with the PLC1 using the Modbus network. The latter is obtained from
only the proportional controller, which is given by the encoder attached to the shaft of the AC motor; the generated
pulses are counted at a 20 [kHz] frequency with the PLC2 and la-
Ω (z) Ka
= . (4) ter transferred to the PLC1 also through the Modbus network; the
Ωd (z) z + (Ka − b) use of the PLC2 is justified since only this device was capable to
accurately read the encoder pulses while spinning at the nominal
Considering the stability condition, −1 < z < 1 (Ogata, 2015), the
shaft speed. Next, the PLC1 computes the control signal u(k), as it
proportional gain can be computed with
was mentioned before, and transforms it into frequency units with
b−1 1+b a range from 0 to 600, which is also sent via the Modbus network
<K< . (5) to the VFD. Finally, such a signal is translated to the range from
a a
0.0 to 60.0 [Hz] and then applied to the AC motor. The block dia-
Once the stability is ensured with the proportional controller, the gram of the full control system is depicted in Figure 3, where all
integral action was added yielding the following closed-loop trans- the signals and the devices processing them are illustrated. In ad-
fer function (Landau and Zito, 2010) dition, the flow diagram of the control system algorithm executed
by the PLC1 is shown in Figure 4, where the interaction between
Ts Ts
Ka z+1− the PLC and other devices can be observed meanwhile it executes
Ω (z) Ti Ti
= , (6) the proposed control system.
Ωd (z)
2
Ts Ts
z + K a − b z + Ka 1 −
Ti Ti
by commanding diverse large steps as the assigned reference va- In order to provide a quantitative result to show which performan-
lues, which are displayed in Table 1. Likewise, it is possible to ce is better among both controllers (Balestrino et al., 2006), the
appreciate the graphic behavior of such changes in Figure 5, and Integral of the Absolute Error (IAE) and the Integral of the Abso-
how the response is established in the assigned value for each case lute value multiplied by Time of the Error (ITAE) were computed
analyzed, nullifying in all cases the steady-state error and yielding as follows
an overshoot of less than 2 %. N
1
IAE = ∑ |e(k) + e(k + 1)|, (8)
Table 1. Large Step References in [RPM] 2 k =0
Stages Desired [RPM] N
1
1 1,800 ITAE = ∑ k|e(k) + e(k + 1)|, (9)
2 50 2 k =0
3 1,100
4 40 where N is the number of samples of the experiment. The resulting
5 1,200 values are listed in Table 2, where it can be observed that the PI
6 20
controller significantly reduces both indices.
1250
1000 140
750 120
500 100
Shaft speed [RPM]
250
80
0
0 20 40 60
80 100 120 140 60
Time [s] 40
Figure 5. Time evolution of the controlled shaft speed with large
set-point changes. 20 d
0
0
30 10
40 20
50 60
Time [s]
3.2 Disturbances Rejection Assessment Figure 6. Time evolution of the controlled shaft speed with P control
and application of disturbances
After validating the proper behavior of the PI controller while fa-
cing sudden changes in the speed reference value, some other ex-
periments were carried out to visualize and compare the responses
40
obtained using both the PI controller and a proportional (P) con-
Control signal [Hz×10]
400 d
140
350 120
250 80
200
60
150
100 40
50 20
0 0
0 10 20 30 40 50 0 10
30 4020 50
Time [s] Time [s]
Figure 8. Time evolution of the controlled shaft speed with PI control Figure 9. Time evolution of the PI control signal with application of
and application of disturbances disturbances
use of only industrial-purpose networked hardware to implement Landau, I. D. and Zito, G. (2010). Digital control systems: design,
more complex control systems rather than traditional high-cost identification and implementation. Springer.
academic-purpose hardware and software, which are less interes-
ting for manufacturing industry in Mexico. Mahesh, K. and Ramachandra, B. (2017). Speed control of three
phase induction motor using PLC under open and closed loop
ACKNOWLEDGEMENTS condition. Int. J. Engineering Research Applications, 7(1): 34–
39.
The authors would like to thank the Universidad Autónoma de https://doi.org/10.9790/9622–0701043439.
la Ciudad de México for supporting the project UACM CCYT-
2023-IMP-02. This contribution is part of the results of such a Moudgalya, K. M. (2007). Digital control. John Wiley & Sons.
project. In addition, the National Researcher Candidate with CVU
Nicola, M. and Nicola, C.-I. (2021). Tuning of PI speed controller
No. 319148 thanks the National Humanities, Science and Techno-
for PMSM control system using computational intelligence. In
logy Council (CONAHCYT) for the support given.
Proc. of IEEE 21st Int. Sym. on Power Electronics (EE). pp: 1–
REFERENCES 6.
https://doi.org/10.1109/Ee53374.2021.9628297 .
Al-Manfi, A. A. (2019). Implementation of remote self-tuning
Ogata, K. (2015). Discrete-time control systems. Pearson.
fuzzy PID controller for induction motor through ethernet. J. of
Electrical Systems, 15(1): 1059–1065. Pramudijanto, J., Ashfahani, A., Fatoni, A., and Nugroho, S. A.
http://journal.esrgroups.org/jes/papers/15_1_2.pdf. (2015). PLC-based PID-predictive controller design for 3-phase
induction motor on centrifugal machine for sugar manufactu-
Avdeev, B. A. and Vyngra, A. V. (2020). Simulation of varia-
ring process. In Proc. of IEEE Int. Con. on Advanced Mechatro-
ble frequency controlled AC induction motor operating on non-
nics, Intelligent Manufacture, and Industrial Automation (ICA-
linear load. In Proc. of IEEE Conf. of Russian Young Resear-
MIMIA). pp: 49–52.
chers in Electrical and Electronic Engineering (EIConRus).
https://doi.org/10.1109/ICAMIMIA.2015.7508001.
pp: 2346–2349.
https://doi.org/10.1109/EIConRus49466.2020.9038963. Saad, N. and Arrofiq, M. (2012). A PLC-based modified-fuzzy
Awdaa, M., Obed, A., and Saleh, A. (2020). Three-phase induc- controller for PWM-driven induction motor drive with constant
tion motor SVPWM-FOC control based on PLC Matlab trans- V/Hz ratio control. Robotics and Computer-Integrated Manu-
lation approach. In IOP Conf. Series: Materials Science and facturing, 28(2): 95–112.
Engineering. pp: 1–8. https://doi.org/10.1016/j.rcim.2011.07.001.
https://doi.org/10.1088/1757-899X/745/1/012028. Ustun, S. V. and Demirtas, M. (2009). Modeling and control of
Balestrino, A., Landi, A., Medaglia, M., and Satler, M. (2006). V/f controlled induction motor using genetic-ANFIS algorithm.
Performance indices and tuning in process control. In Proc. Energy conversion and management, 50(3): 786–791.
of IEEE 14th Mediterranean Conf. on Control and Automation https://doi.org/10.1016/j.enconman.2008.09.032.
(MED). pp: 1–6. https://doi.org/10.1109/MED.2006.328852. Vadi, S., Bayindir, R., Toplar, Y., and Colak, I. (2022). Induc-
Bhatti, S. A., Malik, S. A., and Daraz, A. (2016). Compari- tion motor control system with a programmable logic controller
son of PI and IP controller by using Ziegler-Nichols tuning (PLC) and profibus communication for industrial plants—An
method for speed control of DC motor. In Proc. of IEEE Int. experimental setup. ISA Transactions, 122: 459–471.
Conf. on Intelligent Systems Engineering (ICISE). pp: 330– https://doi.org/10.1016/j.isatra.2021.04.019.
334. https://doi.org/10.1109/INTELSE.2016.7475144. Velagic, J., Kaknjo, A., Osmic, N., and Dzananovic, T. (2011).
Bolton, W. (2021). Instrumentation and Control systems. Newnes. Networked based control and supervision of induction motor
using OPC server and PLC. In Proc. of IEEE Int. Symp. on
Demir, M. C., Cengiz, T., and Eminoğlu, İ. (2019). Ex- Electronics in Marine (ELMAR). pp: 251–255.
perimental setup for AC motor speed control by https://ieeexplore.ieee.org/abstract/document/6044283.
using PLC. J. of Scientific Perspectives, 3(3): 177–
188. https://doi.org/10.26900/jsp.3.018. Verma, P. and Potdar, R. M. (2017). A review on AC drive
controller using PLC. Int. J. of Innovative Research in Elec-
Dukare, P. D. and Jagtap, P. (2017). Induction motor con- trical, Electronics, Instrumentation and Control Engineering,
trol by stator voltage control method using PLC: A review. 5(9):pp: 53–55.
Int. Research J. of Engineering and Technology, 4(4): 92– https://doi.org/10.17148/IJIREEICE.2017.5913.
95. https://www.irjet.net/archives/V4/i4/IRJET–V4I422.pdf.
Howimanporn, S., Thanok, S., Chookaew, S., and Sootkaneung,
W. (2016). Speed control technique for conveyor using PSO
based PID with programmable logic controller. In Proc. of
IEEE/SICE Int. Symp. on System Integration (SII). pp: 670–675.
https://doi.org/10.1109/SII.2016.7844076.
Abstract: Rapid drawdown has been identified as one of the most frequent causes of slope failures due to the effects
associated with drought and operational changes when incorporating hydroelectric plants, which influence the filling
level of earth dams. The main goal of this research is to obtain predictive models based on Artificial Neural Networks
that return the factor of safety of the upstream slope in homogeneous earth dams in the face of the effect of rapid
drawdown. Three geometries and 40 soils were defined to form the embankment, from which hybrid numerical
models of transient water flow with unsaturated soils were built, considering three discharge speeds. From these
results, a database was built to develop the predictive models, by means of the KNIME program and an algorithm
based on Artificial Neural Networks. The behavior of the factor of safety as a function of time is also analyzed to
establish its recovery intervals. Main results show that the minimum factor of safety is obtained between 52 % and
88 % of the total drawdown time. Regarding the predictive models, the adjusted R2 determination coefficients were
greater than 95 % in all cases and the errors remained below 10 %. This demonstrates a high effectiveness of this
algorithm and the validity of its application to geotechnical problems.
Keywords: earth dams, slope stability, factor of safety, rapid drawdown, artificial neural networks
Palabras clave: presas de tierra, estabilidad de taludes, factor de seguridad, desembalse rápido, redes neuronales
artificiales
1. INTRODUCTION1 water resources. The dams built around the world represent
fundamental solutions to meet this objective (Toapaxi et al.,
Since ancient times, the availability of water has been a 2015).
limiting factor for the development of society, which has
forced man to build structures for the use and management of
*[email protected] CC 4.0
Recibido: 12/12/2023
Aceptado:22/03/2024
Publicado en línea: 01/08/2024
10.33333/rp.vol54n1.04
The complex behavior of soil and rock as construction moment); as shown in Equation 1. Normalization has shown
materials limits the understanding of their applications in its usefulness in preprocessing to apply mining techniques,
engineering projects. For this reason, engineering designs are including clustering algorithms, artificial neural networks, and
based on fundamental principles of simplification, where part classification (Kim et al., 2014; Patra & Basudhar, 2003; Singh
of their characteristics is assumed. Sometimes this means that et al., 2023).
the behavior of the soil is underestimated or overestimated in
design work (Shahin, 2016). 𝑋−𝑋𝑚𝑖𝑛
𝑋𝑁 = 𝑋 (1)
𝑚𝑎𝑥−𝑋𝑚𝑖𝑛
embankment. They all have drainage prisms and a foundation simplified from a mathematical and computational point of
considered waterproof and non-deformable, since it is the view (Fredlund et al., 2019).
object of study in the stability analysis.
In this research, the computational simulation proposed is a
The behavior of the embankment was studied in 40 different hybrid type. All the analyzes of the transient leaks were
soil variants (Flores et al., 2023), all classified as high or low obtained using the Finite Element Method (FEM) and the
compressibility clays (CL or CH) according to the Unified Soil stability studies were obtained by the use of the Morgenstern-
Classification System (SUCS), and whose minimum, average Price Method, as one of the most precise among the Limit
and maximum values of the main geotechnical parameters are Equilibrium Methods ( LEM) (El-Hazek et al., 2020; Jasim et
shown in Table 1. al., 2017; Vahedifard et al., 2020). All of the above is done in
the GeoStudio (2018) program.
Table 1. General properties of embankment soils
Specific The result of combining three heights for the embankment (15,
Cohesion Friction Permeability
weight
(kN/m3)
(kPa) angle (°) (m/s) 30 and 40 meters), with their associated geometric
Minimum 16.6 13.8 12.6 4.3 × 10−8
characteristics, as well as 40 combinations of soils and three
discharge speeds (0.1, 0.15 and 0.3 m/day) yield a total of 360
Half 18.3 31.2 17.6 7.7 × 10−7 models. It is taken into account that the FS were obtained daily
Maximum 19.4 67.5 29.5 7.7 × 10−6 and that the relationships between speed and height establish
diversity in terms of duration for each case. Finally, the
The characteristic curves of all modeled soils are included resulting database will have a total of 1.46 × 103 sets of input
(Figure 1), obtained using the method of Aubertin et al. (2003). variables with their respective FS as the response variable.
Figure 2. Implemented work methodology for the use of machine learning in the prediction of FS under rapid drawdown
Figure 3. Tested configurations for calibrating Artificial Neural Network for rapid drawdown
Figure 4. Artificial Neural Network structures implemented for rapid drawdown (a) with 13 input variables and (b) with six input variables
In Figure 5 the decrease in FS during rapid drawdown is neuron adds to the values of the corresponding variable. In the
abrupt. When up to 20 % of the total time has passed the FS first layer there are no high addition relationships between any
has decreased between 50 % to 60 % compared to its initial of the variables. However, in the second layer, these addition
value (always greater than 1.5; therefore, it is considered stable relations are presented in variables 0 and 8. As for the third
and safe). Then between 20 % and 50 % of the total time the layer (output), the red values of neuron 8 imply that the growth
decrease in FS reaches values between 60 % and up to 100 % of this variable is inversely proportional to the value of the FS.
potentially endangering the safety of the slope (if values lower Once the training of each of the applied ANN structures has
than 1 are reached during the process) with the own been carried out the test results are shown between the
repercussions related to its possible failures. However, the numerical model and the prediction obtained for the
maximum recovery of FS does not exceed 30 % of the initial configuration with 13 and six input variables (Figure 6).
value oscillating mostly below 15 %.
From the results shown in Figure 6, a good fit between the
This behavior is sustained for the geometries of 15 and 40 original values and those predicted by the tool for these
meters of embankment height which is reflected in Table 3, analyses can be considered. The robustness of the predictions
which shows the extreme values of FS recovery as a function concerning the results obtained with the numerical model
of time for all the cases studied. guarantees high precision when applying tools based on
Artificial Intelligence such as ANN to solve complex
Table 3. FS recovery intervals as a function of time for the three geometries geotechnical problems such as the stability of slopes against
studied
the effects of rapid drawdown. Furthermore, the results
15 meters 30 meters 40 meters presented in this research guarantee two approaches: a general
52 % 56 % 56 % one, which includes all the soil parameters involved in the
Recovery time interval analysis and a simplified one, which is of vital importance,
88 % 85 % 85 %
since it is known that it is not always possible to have all the
variables geotechnical techniques resulting from soil studies.
Additionally, it was observed that for the same geometry, there The nodes trained from this analysis can be used to make other
is a variation in these intervals when the speed of occurrence predictions for cases not considered in this research.
of the discharge is modified. However, this variation is less
than 5 % in all the cases analyzed and for the three speeds 4. CONCLUSIONS
considered so it is considered irrelevant for these analyses.
From this study, it is recommended to take into account that Rapid drawdown can be caused by phenomena such as drought
the recovery interval of the FS during drawdown is wide or operational changes. Hybrid modeling of this phenomenon
(around 30 % of the time) a behavior that can be associated to obtain FS on earthen dam slopes allows modeling these
with the fluctuation of the physical-mechanical properties of transient water flow problems, considering unsaturated soils. In
the component soils of the embankment, in particular these cases, the reduction of the FS is observed with the
cohesion, a parameter that presents a high variability in the consequent loss of stability and safety of the upstream slope.
study despite the fact that all the soils studied are clays with
similar classifications (CL or CH). A meticulous examination of results derived from numerical
models enables a discerning evaluation of the fluctuation in FS
3.1 Predictive models with ANN during rapid drawdown processes, as outlined in the cited
literature. This underscores the imperative to delve into the
Understanding the internal functioning of ANNs and how the study of this phenomenon and its potential ramifications for the
model is obtained is extremely complex and currently receives optimal operation of reservoir.
a lot of attention (Barredo-Arrieta et al., 2020). In fact, the field
of Explainable AI (XAI) is dedicated to trying to make Thorough analyses and numerical normalization within the 0 to
understandable these models, which achieve high precision but 1 scale reveal distinct trends. The rapid decline in FS is notably
are difficult to understand. For all these reasons, Table 4 shows abrupt, especially within the initial 20 % of the total elapsed
the results corresponding to the biases and the weights time, resulting in reduction of approximately 60 % FS. The
corresponding to the simplified model. pivotal points, marking the minimum FS, occurs within the
50 % to 90 % timeframe of the overall drawdown duration. This
In Table 4, the colors show how the weights and bias are phase is followed by a recuperative period, associated with the
distributed between the neurons of the different layers of the stabilization of pore pressures within the embankment.
ANN. The intermediate layers of ANNs do not have a Substantially, the FS recovery typically does not surpass 15 %
predefined meaning, but rather function as hidden of the previously attained minimum FS.
characteristic variables that the model detects from the data
(Barredo-Arrieta et al., 2020). It is observed (Table 4) that the The intervals of FS recovery exhibit slight variation (less than
weights with small values, close to 0 (in yellow) are the 5 %) for identical geometries, contingent on the drawdown
predominant ones in the model. In this case, this behavior is velocity. Consequently, these intervals are not deemed
observed in the first layer between variable 4 (velocity) and representative for the analyses. The notable fluctuations in FS
neuron 7. Extremely negative values (in red) mean that the recovery values are linked to the variation in the physical-
neuron opposes the values taken by the corresponding mechanical properties of the constituent soils within the
variable. Extremely positive values (in green) mean that the embankment.
Figure 6. Relationship between the numerical model and the prediction with the Artificial Neural Network with 13 and six input variables for the Factor of
Safety with rapid drawdown
El-Hazek, A. N., Abdel-Mageed, N., & Hadid, M. (2020). network). Journal of Engineering Science and
Numerical and experimental modelling of slope stability Technology, 16(5), 3712–3725.
and seepage water of earthfill dam. Journal of Water and Jasim, F. H., Vahedifard, F., Ragno, E., AghaKouchak, A., &
Land Development, 44(I_III), 55–64. Ellithy, G. (2017). Effects of Climate Change on
https://doi.org/10.24425/jwld.2019.127046.Abstract Fragility Curves of Earthen Levees Subjected to
Flores, I., Haramboure, Y. G., Días, J., & Tristá, J. G. (2022). Extreme Precipitations. Geo-Risk, 498–507.
Estudio de la estabilidad de taludes ante procesos de Jong, S. C., Ong, D. E. L., & Oh, E. (2021). State-of-the-art
desembalse rápido en presas. Ingeniería Hidráulica y review of geotechnical-driven artificial intelligence
Ambiental, XLIII(03), 64–72. Available in: techniques in underground soil-structure interaction.
https://search.ebscohost.com/login.aspx?direct=true&d Tunnelling and Underground Space Technology, 113.
b=edsdoj&AN=edsdoj.5af06550472049faafda0b4e8ff9 https://doi.org/10.1016/j.tust.2021.103946
49ff&lang=es&site=eds-live&scope=site. Kim, J., Kwon, O.-I., Baek, Y., & Kim, C.-Y. (2014). A Study
Flores, I., Hernández, F. M., Haramboure, Y., & Tristá, J. on the Variation of Ground Safety Factor by Earthworks.
(2023). Relationship between soil physical-mechanical The Journal of Engineering Geology, 24(3), 333–341.
parameters and earth dam slope stability using stochastic https://doi.org/10.9720/kseg.2014.3.333
and numerical modeling. Geociencias, 42(3), 387–402. Komasi, M., Mohammadzadeh, A., & Beiranvand, B. (2019).
https://doi.org/ 10.5016/geociencias.v42i3.17683 Optimization of horizontal drain dimensions in
Flores, I., Tristá, J. G., & Haramboure, Y. G. (2020). heterogeneous earth dams using Artificial Neural
Estabilidad de taludes durante un desembalse rápido en Network (ANN) (Case study : Marvak dam). Journal of
presas de tierra con suelos parcialmente saturados. Applied Research in Water and Wastewater, 12, 109–
Ingeniería y Desarrollo, 38(1), 13–31. 116. https://doi.org/10.22126/arww.2019.3817.1101
http://www.scielo.org.co/scielo.php?script=sci_arttext Maneta, M., & Shnabel, S. (2003). Aplicación de Redes
&pid=S0122- Neuronales Artificiales para determinar la distribución
34612020000100013&lng=en&nrm=iso&tlng=es%0A espacial de la humedad del suelo en una pequeña cuenca
http://www.scielo.org.co/scielo.php?script=sci_abstract de drenaje. Estudios preliminares. Estudios de La Zona
&pid=S0122- No Saturada Del Suelo, VI, 295–304.
34612020000100013&lng=en&nrm=iso&tlng=es Montoya, E. (2018). Metodología para la aplicación de Redes
Flores, I., Tristá, J. G., & Haramboure, Y. G. (2021). Neuronales Artificales para sistemas de alerta de
Evaluación de la estabilidad de taludes en presas de deslizamientos provocados por lluvias en regiones
tierra empleando Redes Neuronales Artificiales. Revista montañosas. PhD Thesis. Universidad de Brasilia.
de La Universidad de Zulia, 32, 261–283. Nanehkaran, Y. A., Licai, Z., Chengyong, J., Chen, J., Anwar,
https://doi.org/http://dx.doi.org/10.46925//rdluz.32.17 S., Azarafza, M., & Derakhshani, R. (2023).
Fredlund, D., Xing, A., & Huang, S. (1994). Predicting the Comparative Analysis for Slope Stability by Using
permeability function for unsaturated soils using the Machine Learning Methods. Applied Sciences
soil-water characteristic curve. Canada Geotechnical (Switzerland), 13(3), 1–14.
Journal, 31, 533–546. https://doi.org/10.3390/app13031555
Fredlund, D. G., & Xing, A. (1994). Equations for the soil- Patra, C. R., & Basudhar, P. K. (2003). Generalized solution
water characteristic curve. Canadian Geotechnical procedure for automated slope stability analysis using
Journal, 31, 533–546. inclined slices. Geotechnical and Geological
Fredlund, M., Mckeown, R., Lu, H., & Gitirana, G. (2019). Engineering, 21(3), 259–281.
Comparison of 2D and 3D anchor analysis https://doi.org/10.1023/A:1024923329646
methodologies. 1999, 359–367. Pinyol, M., Alonso, E. E., & Olivella, S. (2008). Rapid
https://www.researchgate.net/publication/340438381_C drawdown in slopes and embankments. Water
OMPARISON_OF_2D_AND_3D_ANCHOR_ANAL Resources Researh, 44, 1–22.
YSIS_METHODOLOGIES https://doi.org/10.1029/2007WR006525
Ganasan, E., John, A., & Meng, L. (2021). Settlement Pinyol, N., Alonso, E. E., Corominas, J., & Moya, J. (2013).
Prediction Model in Consideration of Static Loading on Influencia del desembalse en la estabilidad de una
Soft Clay by Utilising Machine Learning Method. ladera. Caso Real. 37–39. VIII Simposio Nacional sobre
Recent Trends in Civil Engineering and Built Taludes y Laderas Inestables. Barcelona. Available in:
Environment, 2(1), 551–560. https://search.ebscohost.com/login.aspx?direct=true&d
https://doi.org/https://doi.org/10.30880/rtcebe.2021.02. b=edsair&AN=edsair.dedup.wf.001..57c83103345e083
01.060 529765aeac339b1f9&lang=es&site=eds-
Haramboure, Y. G., Flores B., I., & Tristá, J. G. (2021). Efecto live&scope=site
de desembalse en la estabilidad de presas de tierra : dos Riedmiller, M., & Braun, H. (1993). A Direct Adaptive
casos de estudio en Cuba . Ingeniería Hidráulica y Method for Faster Backpropagation Learning : The
Ambiental, XLII(1), 42–53. Available in: RPROP Algorithm. IEEE Access, 16, 586–591.
https://search.ebscohost.com/login.aspx?direct=true&d https://doi.org/O-7803-0999-5/93/$03
b=edsdoj&AN=edsdoj.bc970eb6c00a4b9ca277830ee2d Romer, G., Ayala, L., Chagas, F., & Silva, D. A. (2019).
a395c&lang=es&site=eds-live&scope=site Estudio de Estabilidad de Taludes en Presas de Tierra en
Jabbar, A., & Muataz, I. (2021). Stability and Seepage of Earth Condiciones No Saturadas Cuando Ocurre Descenso
Dams with Toe Filter (Calibrated with artificial neural Rápido de Nivel de Agua. Geotechnical Engineering in
the XXI Century: Lessons Learned and Future Faculty of Civil. He has more than 12 publications and has
Challenges, 1804–1814. participated in 8 events. both national and international. She is
https://doi.org/10.3233/STAL190237 a reviewer for several national and international journals. She
Salazar, F., Toledo, M. A., Oñate, E., & Morán, R. (2015). An has published a monograph that has an ISBN.
empirical comparison of machine learning techniques
for a dam behaviour modelling. Structural Safety, 59,
118. González Haramboure, Yoermes,
https://doi.org/https://doi.org/10.1016/j.strusafe.2015.0 Graduated in Hydraulic
5.001 Engineering (2005) from the
Santillán, D., Fraile-Ardanuy, J., & Toledo, M. Á. (2014). Technological University of
Predicción de lecturas de aforos de filtraciones de presas Havana “José Antonio Echeverría”
bóveda mediante redes neuronales artificiales. (Cujae). Master in Civil
Tecnología y Ciencias Del Agua, V(3), 81–96. Available Engineering (Geotechnics
in: Specialty) (2010) and Doctor in
https://search.ebscohost.com/login.aspx?direct=true&d Technical Sciences (Since 2014,
b=edsair&AN=edsair.doajarticles..987aeb382aafcbb0b Civil Specialty). He has been a Full
a932817dc676ed1&lang=es&site=eds-live&scope=site Professor since 2022. He has received the “Alma Mater”
Shahin, M. A. (2016). State-of-the-art review of some artificial award, “Golden Chalk” and Distinction for Cuban Education,
intelligence applications in pile foundations. Geoscience among others. He has 19 publications in referenced magazines
Frontiers, 7, 33–34. and has participated in 20 events, both national and
https://doi.org/https://doi.org/10.1016/j. international. He is a reviewer in two international journals and
gsf.2014.10.002. has published 6 monographs that have ISBNs.
Singh, P., Bardhan, A., Han, F., Samui, P., & Zhang, W.
(2023). A critical review of conventional and soft
computing methods for slope stability analysis. García Tristá, Jenny, Graduated
Modeling Earth Systems and Environment, 9(1), 1–17. in Civil Engineering (2006) at the
https://doi.org/10.1007/s40808-022-01489-1 Technological University of
Toapaxi, J., Hidalgo, X., Valencia, N., & Castro, M. (2015). Havana “José Antonio Echeverría”
Bases y Criterios para el Dimensionamiento Hidráulico (Cujae). Doctor in Technical
de Formas Especiales de Saltos en Esquí en Presas. Sciences (2016, Civil Specialty).
Revista Politécnica, 35(3), 1–6. Available in: She has been a Full Professor since
https://search.ebscohost.com/login.aspx?direct=true&d 2022. She has received the
b=edsair&AN=edsair.doajarticles..898fe082a53ba82ae Forjadores Stamps of the Future
7c3ec4762a8c503&lang=es&site=eds-live&scope=site Young Researcher award and the
Vahedifard, F., Jasim, F. H., Tracy, F. T., Abdollahi, M., distinction “For Cuban education” awarded by the Ministry of
Alborzi, A., & Aghakouchak, A. (2020). Levee Fragility Higher Education (MES). He has more than 20 publications
Behavior under Projected Future Flooding in a Warming and has participated in 8 events, both national and
Climate. Journal of Geotechnical and international. She is a reviewer for several national and
Geoenvironmental Engineering, 146(12), 1–12. international journals.
https://doi.org/10.1061/(ASCE)GT.1943-5606.0002399 Rosete Suárez, Alejandro,
Wengang, Z., Hanlong, L., Wang, L., Xing, Z., & Yanmei, Z. Graduated as an Engineer in
(2023). Application of Machine Learning in Slope Automated Management Systems
Stability Assessment (Springer (ed.)). (1993) at the Technological
https://doi.org/https://doi.org/10.1007/978-981-99- University of Havana “José
2756-2 Antonio Echeverría” (Cujae).
Master in Computer Science
Applied to Engineering and
BIOGRAPHIES Architecture (1995) and Doctor in
Technical Sciences (2000,
Flores Berenguer, Isaida, Computer Science Specialty). He has received several awards,
Graduated in Civil Engineering including three Awards from the Cuban Academy of Sciences,
(2016) at the Technological two Awards from the Ministry of Higher Education (MES),
University of Havana “José and three Special Distinctions from the MES, among others.
Antonio Echeverría” (Cujae). President of the Permanent Court (Doctorates) of Automation
Master in Civil Engineering and Computing from 2019 to 2022 and member of that court
(Communication Routes Specialty) since 2005. Full Academician of the Cuban Academy of
(2020) and Aspiring for a Doctorate Sciences since 2018.
in Technical Sciences (Since 2022,
Civil Specialty). She has been an
Assistant Professor since 2022. She
has received the Young Researcher award awarded by the
1
Universidad Tecnológica de Panamá, Facultad de Ingeniería Civil, Panamá, Panamá
2
Universidad Tecnológica de Panamá, Centro Experimental de Ingeniería, Panamá, Panamá
3
Secretaría Nacional de Ciencia Tecnología e Innovación de Panamá, Sistema Nacional de Investigación, Panamá,
Panamá
Resumen: Los suelos expansivos son característicos de regiones tropicales, presentando grandes cambios
volumétricos cuando su contenido de agua varía; estos cambios suelen presentar problemas estructurales en las
cimentaciones provocados directamente por el comportamiento del material sobre el cual han sido construidas. Para
controlar la expansión en suelos se emplean materiales cementantes, como el carbonato de calcio (CaCO3), que es un
mineral que se encuentra en muchas partes del planeta. Los exoesqueletos de animales marinos, residuos de
actividades industriales como la pesca y la extracción de arena marina, que pudieran reutilizarse como una fuente
rica en CaCO3. El objetivo de este estudio es evaluar la viabilidad para la estabilización de suelos expansivos en
Panamá, a través de la compactación controlada y la adición de estabilizantes como la concha marina y el cemento,
y así determinar la proporción óptima en el suelo expansivo del estudio, atendiendo al ODS 12 “Producción y
consumo responsable”. Se caracterizó el suelo expansivo, la composición química del suelo y de la concha de mar,
y se determinó el esfuerzo máximo a soportar a través del ensayo de la compresión no confinada del suelo expansivo
con las adiciones realizadas, para determinar la mezcla y proporción óptima de estabilización. La mezcla de
cemento/concha marina al 35 % obtuvo los mejores resultados, esto muestra viabilidad para el empleo de los residuos
de concha marina como material cementante para la estabilización de suelos expansivos en Panamá.
Palabras clave: Conchas marinas, suelo expansivo, reutilización, estabilización, compresión no confinada del suelo
*[email protected] 10.33333/rp.vol54n1.05
Recibido: 05/01/2024 CC 4.0
Aceptado:10/05/2024
Publicado en línea: 31/08/2024
Figura 3. Concha de mar pulverizada sin lavar Figura 4. Prueba de expansión del suelo de Metetí, Darién
mayor porcentaje de CaO en comparación a la muestra de pesquera. El reemplazo del uso de materiales como las calizas
concha sin lavar. permite caminar a una nueva era de construcción sostenible sin
la explotación de recursos limitados a través de un ciclo
Al adicionar un 35 % de concha pulverizada sin lavar a la económico verde (Chilakala et al., 2019).
muestra de suelo expansivo, se logró mejorar la capacidad de
carga intrínseca del mismo en un 13,7 %; mientras que se Considerando que los suelos expansivos predominan en
dieron mejoras de hasta un 27,8 % en la mezcla de suelo al Panamá, el uso de la concha marina pulverizada es una
agregarle cemento al 10 %. alternativa mucho más factible como estabilizador de
expansión contrapuesto, al remplazo del suelo en su totalidad
Además, tomando en cuenta los resultados presentados en la por suelos no expansivos, ya que no sólo disminuirá costos de
Figura 6 de esfuerzo deformación, se puede inferir que la transporte y mano de obra porque es un recurso de fácil acceso,
mezcla de S + 35CS tuvo el comportamiento óptimo entre sino que permitirá aumentar la resistencia a la compresión y la
todas las mezclas analizadas, 4,5 veces más el esfuerzo durabilidad del material.
máximo del suelo expansivo seco. Esto se debe principalmente
a la elevada presencia de CaO en la muestra, como resultado
de la combinación de la concha pulverizada y el cemento.
Con base en los resultados de compresión no confinada, se Conforme en los resultados de compresión no confinada, se
puede inferir que la mejora con respecto al aumento de la puede observar que la mezcla del 35 % de cemento junto con
capacidad de carga en los cilindros ensayados aumenta concha pulverizada sin lavar (en partes iguales) tuvo una
proporcionalmente con el porcentaje adicionado de concha capacidad de carga mucho mayor de 15,43 kg/cm2 a
marina pulverizada, lo que indica una posible neutralización comparación de la capacidad de soporte del suelo expansivo
de la expansión propia de este tipo de suelos. con cemento que fue de 4,37 kg/cm2. Lo que supone una
mejora a la capacidad máxima del suelo expansivo natural.
Cabe destacar que el aumento de proporción de concha
pulverizada posee un límite de efectividad en la mezcla y El agregado fino natural puede ser sustituido por conchas
variará dependiendo de no solo las propiedades del suelo marinas trituradas hasta en un 35 % en peso para producir
ensayado, sino de la composición química en lo que respecta concreto ordinario estándar con la resistencia adecuada y
al porcentaje de carbonato de calcio. densidad. Esto demuestra la viabilidad del estudio y su
potencial de aplicación en el día a día de la industria estructural
Es posible que la efectividad del concreto esté ligada a la (Anggraini et al., 2023; Chilakala et al., 2019; Ruíz & Farfán,
presencia de aditivos externos, sin embargo, las propiedades 2016; Vinod et al., 2020).
intrínsecas del CaCO3 propias de la concha como agente
cementante contribuyen a la durabilidad de este y a la Para el entendimiento del comportamiento a futuro de esta
disminución del contenido de vacíos. mezcla, se recomienda la evaluación de todas las propiedades
de este agregado como componente estabilizante en los suelos
5. CONCLUSIÓN expansivos, el estudio de la cohesión, densidad conjunta y
ensayos de capacidad de carga a escala. Además, se
Las conchas marinas son un recurso natural renovable que pasa recomienda utilizar la concha lavada para aumentar el
desapercibido, el uso de éste como materia prima de carbonato porcentaje de óxido de calcio, y con un secado a mayor
de calcio en los suelos y cementos permite tener un nuevo temperatura, para aumentar la capacidad de cementación, y por
enfoque en el uso de residuos de la industria arenera y lo tanto la capacidad de carga máxima.
Jonatha, Arrocha, es
Delvis, Hernández, Estudiante Licenciado en Edificaciones.
de la carrera de Ingeniería Técnico del Laboratorio de
Geológica de la Universidad Geotecnia (LABGEO) del
Tecnológica de Panamá, Centro Experimental de
Panamá. Fue olímpico y Ingeniería de la Universidad
finalista representante de la Tecnológica de Panamá (UTP).
República de Panamá en las II Tiene experiencia en
Olimpiadas Centroamericanas y perforaciones de pozos para
del Caribe de Biología, 2018. investigación geotécnica,
Trabajó en la División de ensayos de campo y laboratorio
Dragado del Canal de Panamá, de mecánica de suelos y mecánica de rocas. Ha
2019. Es investigador del Centro de Investigaciones contribuido en la elaboración de publicaciones científicas
Arqueológicas del Istmo Fundación El Caño, Panamá. en revistas indexadas en SCOPUS. Es Investigador de
Además, ha contribuido en una publicación científica en Proyectos de Investigación financiados por la SENACYT,
el área de geoquímica en la revista I+D Tecnológico de la y asesor de proyectos de Geotecnia a estudiantes de la
Universidad Tecnológica de Panamá. Áreas de interés: UTP.
geología, petrología, mineralogía y minería.
1
Pontificia Universidad Católica del Ecuador - Sede Esmeraldas PUCESE, Programa de Maestría en Electricidad,
Esmeraldas, Ecuador
2
SDAS Research Group, Ben Guerir 43150, Moroco
Resumen: Este documento presenta una investigación sobre el diseño y evaluación de un sistema de recolección de
energía piezoeléctrica utilizando una alfombra con múltiples elementos piezoeléctricos. La investigación inicia con
una introducción sobre la importancia de desarrollar fuentes alternativas de energía renovable para dispositivos
electrónicos portátiles. Posteriormente, se presenta una revisión de literatura sobre tecnologías de recolección de
energía del movimiento humano, enfocándose en materiales piezoeléctricos. La metodología describe el diseño de
un prototipo de alfombra piezoeléctrica compuesta por una capa de goma, una lámina de acrílico y múltiples
elementos piezoeléctricos conectados en paralelo a un circuito recolector de energía. Se realizaron pruebas para
caracterizar la respuesta de los sensores piezoeléctricos y evaluar el sistema de recolección de energía bajo
diferentes configuraciones. Los resultados indican que un solo elemento generó 18,59 μJ, mientras que 10
elementos conectados produjeron 297,4 μJ con un voltaje de 0,57 V. Además, el sistema permite la adquisición
remota de datos a través de la plataforma IoT Thingspeak.
Palabras clave: Energía renovable; Recolección de energía; Movimiento humano; Alfombra piezoeléctrica; IoT
Keywords: Enewable energy; Energy harvesting; Human movement; Piezoelectric mat; IoT
1
1. INTRODUCCIÓN Sin embargo, estos no son los únicos problemas energéticos
actuales. Aunque los hidrocarburos han sido la principal
El uso y abuso ilimitado de recursos energéticos finitos como fuente de energía, se ha demostrado que no son infinitos y su
el carbón, los hidrocarburos y el gas natural, tanto en hogares uso excesivo está dañando el medio ambiente a escala global.
como en grandes industrias, están teniendo un impacto Por lo tanto, el sistema energético actual no es sostenible y
negativo en los ciclos vitales del planeta como el ciclo del varias naciones buscan cambiar su matriz energética a fuentes
agua. Esto conduce a fenómenos como lluvias ácidas, limpias para lograr mayor eficiencia energética mediante
desertificación, sequías, contaminación de ríos y océanos, y alternativas renovables (Niasar et al., 2020).
la degradación de la capa de ozono (Abdal & Leong, 2019).
*[email protected]
Recibido: 23/11/2023
Aceptado:19/07/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.06
CC 4.0
Existen métodos para minimizar los efectos de la actividad las tecnologías más comunes para recolectar energía
humana sobre el medio ambiente. En este contexto, las biomecánica del movimiento humano son el Piezoeléctrico,
fuentes de energía limpia emergen como una alternativa el electromagnético y el triboeléctrico, respectivamente (Lin
importante. Estas fuentes se producen de manera constante y et al., 2017; Moure et al., 2016; Saha et al., 2008).
no se agotan a una escala significativa para los seres
humanos. Se originan a partir de fenómenos naturales que no Los materiales piezoeléctricos funcionan según el principio
alteran el equilibrio de la Tierra, no generan residuos y su del efecto piezoeléctrico, es decir, la propiedad del material
tasa de consumo es menor que la velocidad de regeneración para generar electricidad tras la deformación mecánica (Kang
de los recursos utilizados (Huang et al., 2022). et al., 2016). Estos materiales generalmente caen en los
modos de trabajo transversal, longitudinal y piezotrónico
Las fuentes renovables de energía sobresalen debido a que (Narita & Fox, 2018). Varios investigadores han utilizado
son inagotables, renovables y no contaminantes. Además, estos modos de trabajo para recolectar energía del
según algunos estudios (Vizzari et al., 2021), el uso a largo movimiento humano in vitro e in vivo para alimentar
plazo de energías renovables conserva la vida humana al dispositivos portátiles inteligentes como zapatos inteligentes,
tiempo que promueve sociedades más sostenibles y textiles inteligentes, dispositivos implantables y dispositivos
equitativas. de monitoreo biomédico (Zhang et al., 2019). Del mismo
modo, los materiales electromagnéticos basados en la
Con el rápido desarrollo de la microelectrónica, los captación de energía del movimiento humano son bien
dispositivos electrónicos inteligentes, se ha minimizado el conocidos por su mayor potencia de salida. Los dispositivos
uso de fuentes de energía externas, como baterías (Fan et al., de captación de energía electromagnética funcionan según la
2016) . No obstante, existe un interés considerable para ley de inducción electromagnética de Faraday, por lo tanto, la
mejorar el almacenamiento de la batería, ya que los sistemas propiedad del material de inducir corriente eléctrica cuando
antes descritos aún requieren una carga regular y un proceso atraviesa el campo magnético (Zhang et al., 2019).
de reemplazo que reduce significativamente sus aplicaciones
(Song et al., 2019). Por lo tanto, para alimentar aparatos Dependiendo de su mecanismo, los materiales piezoeléctricos
electrónicos de baja potencia se necesita un suministro de se pueden clasificar en recolectores de energía rotativos,
energía constante y sostenible. La solución más factible es oscilatorios e híbridos (Dinulovic et al., 2015). Los colectores
extraer energía del medio ambiente y transformarla en rotativos suelen utilizar un tren de engranajes y turbinas
electricidad (Wang et al., 2018). debido a sus requisitos de alta potencia (Niroomand &
Foroughi, 2016). Los recolectores oscilatorios usan el sistema
El estudio presentado por Owusu & Asumadu-Sarkodie masa-resorte-amortiguador que generalmente funciona en
(2016) menciona a la energía solar y eólica como fuentes de excitaciones de baja frecuencia con salidas de baja potencia,
energía renovable y limpia. Sin embargo, estas no y los recolectores electromagnéticos híbridos usan rotores
proporcionan energía constante debido a su dependencia de excéntricos para acoplar la energía cinética producida por el
las condiciones climáticas. Por el contrario, el movimiento movimiento humano (Wu et al., 2018).
humano está presente en todas partes, exhibiendo un gran
potencial para generar energía biomecánica ilimitada, limpia En el contexto de dispositivos electrónicos portátiles
y sostenible para la electrónica inteligente (Xie & Cai, 2014). autoalimentados, los recolectores de energía electromagnética
que aprovechan el movimiento humano, suelen tener
Los movimientos biomecánicos humanos, como el dimensiones más amplias y carecer de diseños
movimiento de los dedos, caminar, correr, escribir a máquina miniaturizados. Como resultado, a menudo se combinan con
e incluso los desplazamientos diminutos dentro del cuerpo otros materiales, como materiales piezoeléctricos y
humano, como el movimiento de los pulmones, los latidos triboeléctricos, para su integración en la electrónica portátil
del corazón y la contracción muscular, poseen una cantidad (Quan et al., 2015). Finalmente, la tecnología más reciente y
abundante de energía residual que se puede aprovechar avanzada para el aprovechamiento de energía del movimiento
(Roundy & Trolier-McKinstry, 2018). Pequeños humano es la recopilación de energía triboeléctrica, que
desplazamientos como el movimiento de las extremidades opera en el conocido fenómeno llamado triboelectrificación
inferiores, como el balanceo de las piernas, puede generar la en el que los materiales se cargan eléctricamente al contacto
máxima potencia debido a sus pares más altos que los por fricción (Wang et al., 2015). El desplazamiento relativo
movimientos de otras partes del cuerpo (Fan et al., 2017). Sin entre materiales cargados eléctricamente da como resultado
embargo, aprovechar la energía del movimiento humano es la generación de electricidad. Se han propuesto cuatro tipos
bastante complejo y desafiante debido a los movimientos de modos de trabajo, en función del tipo de movimiento
multidimensionales y la frecuencia ultra baja del movimiento relativo entre los dos materiales triboeléctricos y la posición
humano (Kim et al., 2018). del electrodo (Randriantsoa et al., 2023).
Bajo este contexto, se necesitan tecnologías de recolección de Los materiales piezoeléctricos, parte de la categoría más
energía que puedan recopilar de manera efectiva el potencial amplia de ferroeléctricos, poseen una característica distintiva:
del movimiento humano. El progreso en el desarrollo de su estructura molecular permite la existencia de dipolos
dispositivos electrónicos portátiles autoalimentados depende eléctricos locales (Zhao & You, 2014). Estos materiales
considerablemente del rendimiento de las tecnologías de tienen la capacidad de convertir el estrés mecánico en carga
recolección de energía (Alhawari et al., 2018). Actualmente, eléctrica -efecto piezoeléctrico directo- y viceversa -efecto
piezoeléctrico inverso- (Cao et al., 2019). Trascendentales en de energía piezoeléctrica. Los autores destacan la importancia
la generación de energía de forma sostenible, especialmente de las soluciones de captación de energía piezoeléctrica en el
mediante la utilización de cerámicas piezoeléctricas. Sin desarrollo de redes de sensores inalámbricas autoalimentadas.
embargo, el almacenamiento de energía convertida es un Además, el estudio compara diferentes modelos para
desafío significativo (Wen & Xu, 2018). Aunque se puede soluciones de captación de energía piezoeléctrica, incluyendo
utilizar directamente, la energía generada por un dispositivo modelos analíticos, modelos numéricos y modelos
piezoeléctrico no resonante, que no produce una señal experimentales. Asimismo, se discute las ventajas y
constante, debe almacenarse, generalmente, en una batería. desventajas de cada modelo y proporciona un análisis
En este sentido, los captadores de energía que aprovechan las detallado de su rendimiento. La conclusión de este estudio
fuentes de energía mecánicas y vibratorias, incluyendo el dicta que la selección del modelo apropiado para las
movimiento humano, presentan un futuro prometedor para la soluciones de captación de energía piezoeléctrica depende de
creación de sistemas energéticos autoalimentados en diversos la aplicación específica y del nivel de precisión deseado. Por
entornos, que incluyen el cuerpo humano (Qian et al., 2018). lo que, los autores sugieren que la investigación futura
debería centrarse en el desarrollo de modelos más precisos y
El presente trabajo de investigación tiene como objetivo eficientes para las soluciones de captación de energía
principal desarrollar un dispositivo que no solo sea funcional piezoeléctrica con el fin de mejorar su rendimiento y
y eficiente en la conversión de energía, sino también fiabilidad.
fácilmente transportable y económico, satisfaciendo así las
necesidades y demandas de un público más amplio. Para lo Chen et al. (2020) exploran la posibilidad de utilizar
cual, se emplea una metodología secuencial que permite el materiales piezoeléctricos para recolectar energía de las
diseño tanto a nivel de hardware y software de la alfombra vibraciones generadas por el viento. Este trabajo también
piezoeléctrica compuesta por una capa de goma, una lámina examina la factibilidad de emplear la energía obtenida
de acrílico y múltiples elementos piezoeléctricos conectados mediante piezoelectricidad para alimentar nodos de sensores
en paralelo a un circuito recolector de energía. inalámbricos. La investigación revela las conclusiones de
Posteriormente, se realizaron pruebas para caracterizar la experimentos realizados con el fin de establecer la
respuesta de los sensores piezoeléctricos y evaluar el sistema configuración más eficiente de los transductores
de recolección de energía bajo diferentes configuraciones. piezoeléctricos para la recolección de energía. Los autores
Donde se obtuvo, que un solo elemento generó 18,59 μJ, del estudio tratan también los desafíos que se presentan en el
mientras que 10 elementos conectados produjeron 297,4 μJ diseño de sistemas de recolección de energía, como la
con un voltaje de 0,57 V. necesidad de una administración eficaz de la potencia y el
almacenamiento de energía. En última instancia, se concluye
El artículo se encuentra estructurado en cuatro secciones que la recolección de energía mediante piezoelectricidad
distintas. En la segunda sección, se describen los trabajos puede ser una tecnología con gran potencial para alimentar
relacionados, mientras que la tercera sección presenta la nodos de sensores inalámbricos, si bien se requieren más
metodología empleada en la investigación. En la cuarta investigaciones para optimizar el rendimiento de estos
sección, se exponen los resultados de las predicciones sistemas.
realizadas, acompañados de su correspondiente análisis.
Finalmente, dentro de la sección 5 se abordan las discusiones Qian et al. (2019) exploran el uso de la hebilla de correa
y conclusiones del artículo. amplificada mecánicamente para una mochila utilizando un
actuador piezoeléctrico apilado para la extracción de energía,
2. TRABAJOS RELACIONADOS y se logra una potencia de salida máxima de 0,4 mW.
Existen varios estudios que abordan el problema desde Guido et al. (2016) abordan la captación de energía
diferentes perspectivas, por ejemplo, Ali et al. (2023) revisan biocompatible de los movimientos humanos mediante el
los avances recientes en materiales piezoeléctricos apilamiento de películas delgadas de nitruro de aluminio con
biodegradables, incluidos los polímeros naturales y sintéticos, un voltaje pico de 0,7 V y una potencia eléctrica de 0,2 uW.
sus principios, aplicaciones avanzadas y retos para usos Asimismo, Kim et al. (2018) realizan la mejora de un
médicos. El artículo destaca la evolución de los materiales recolector de energía flexible para fuerzas de entrada
convencionales no reciclables a materiales funcionales, multidireccionales, fabricado con películas de PDMS y
biodegradables, ecológicos y biocompatibles. También se PVDF con un pico de tensión de salida generado de 1,75 V.
analiza los métodos de degradación de estos materiales
piezoeléctricos mediante estudios in vitro e in vivo. El Zhao & You (2014) recomiendan el uso de una multicapa de
artículo concluye que estas mejoras en los materiales PVDF para diseñar PEH a partir de energía mecánica en
piezoeléctricos biodegradables y los microsistemas podrían zapatos, la cual se origina por el movimiento humano con
permitir nuevas aplicaciones en el campo biomédico, como la una potencia de salida de 1 mW. De forma similar, Cao et al.
detección de señales fisiológicas, el diagnóstico de problemas (2019) desarrollan un colector insertado en el zapato con un
médicos, la apertura de la barrera hematoencefálica y la diseño estructural en zigzag compuesto por dos capas de
estimulación de la cicatrización o el crecimiento de tejidos. acero inoxidable y una película de PVDF. Este diseño se
implementó en un humano corriendo en una cinta a 6 km/h,
Por otro lado, Raj & Ranganathan (2021) abordan el análisis lo cual produjo un pico de tensión de salida de 5,032 V.
comparativo de distintos modelos de soluciones de captación
Wen & Xu (2018) diseñan un recolector de energía para los 3. METODOLOGÍA PARA LA GENERACIÓN DE
pasos humanos durante la marcha basado en el mecanismo de ENERGÍA PIEZOELÉCTRICA
amplificación de fuerza, con el fin de amplificar la fuerza
aplicada en la pila piezoeléctrica y se examina la alta Para el desarrollo de un sistema de recolección de energía, se
potencia de salida en comparación con los diseños existentes. ideó un prototipo detallado cuya representación se puede
observar en la Figura 1.
Qian et al. (2018) utilizan un mecanismo de amplificación de
fuerza y se fabrica un colector de energía para el calzado.
Este consiste en varias pilas piezoeléctricas y dos placas de
aluminio para caminar a diferentes velocidades. En
consecuencia, se obtiene una potencia máxima de salida de
28 mW. Del mismo modo, se han realizado diferentes
estudios sobre el mecanismo de amplificación de la fuerza
con el fin de maximizar la potencia de salida en los
dispositivos de recolección de energía piezoeléctrica.
𝑔33𝑙𝐹
𝑉= (1)
𝜋𝑟 2 Figura 4. Placa PCB de alfombra piezoeléctrica
4. RESULTADOS
Figura 5. Prueba de comportamiento de carga de la batería
Luego de diseñar el prototipo de la alfombra piezoeléctrica,
3.2 Sistema de Adquisición de datos se realizaron diversas pruebas de rendimiento en conjunto
con la plataforma de Medición de Internet de las Cosas (IoT,
Se utilizaron herramientas electrónicas y de conectividad a por sus siglas en inglés) ThingSpeak. Este análisis se enfocó
Internet para adquirir las mediciones de la alfombra en dos áreas principales: en primer lugar, se examinaron los
piezoeléctrica. El dispositivo empleado en este proceso está resultados de las mediciones de las propiedades del voltaje de
compuesto por varios componentes clave, que incluyen un salida del material piezoeléctrico; en segundo lugar, se
módulo Node MCU ESP 8266, un sensor de voltaje INA 219 evaluaron los resultados del sistema de recolección de
y un display IC2 16X2, como se detalla en la Tabla 1. energía. A través de la evaluación de estas métricas clave, se
busca obtener una caracterización completa de las
Tabla 1. Características técnicas de circuito de adquisición de datos
Componente Voltaje de Corriente de Forma de
capacidades y las posibles aplicaciones de la alfombra
Operación Operación comunicación piezoeléctrica para la recolección de energía.
Generador 1,2 – 5 V Varía según Análogo
Piezoeléctrico modelo 4.1 Mediciones de las características de la tensión de
INA 219 3,3V o 5 V Varía según Análogo salida piezoeléctrica
aplicación
I2C 16X2 5V Varía según Serial Los resultados obtenidos en estas mediciones muestran un
uso voltaje mínimo de 0,1 V y un máximo de 2,8 V luego de
efectuar pruebas con 10 muestras, el voltaje se generaba
cuando el piezoeléctrico estaba sometido a una carga y
La adquisición de datos de voltaje y corriente se lleva a cabo
desaparecía o llegaba a cero en ausencia de cambios de
mediante la codificación de un sistema que involucra el
sensor INA 219 y la placa de desarrollo Node MCU ESP presión en el mismo, lo que se puede clasificar como un
impulso. Sin embargo, los resultados de las pruebas con
8266. En este proceso, la placa se configura con credenciales
diferentes pesos mostraron una inconsistencia en el voltaje de
de red específicas, incluyendo Identificación de Conjunto de
salida para un elemento en particular. La alfombra
Servicios (SSID, por sus siglas en inglés) y contraseña, para
piezoeléctrica de 30 x 30 cm, como se muestra en la Figura 7,
permitir la transferencia de los datos recopilados a la
fue empleada para ajustar la trayectoria de caída sobre el
plataforma ThingSpeak. El diagrama de este procedimiento
piezoeléctrico, la cual no logró asegurar una respuesta
se encuentra detallado en el algoritmo 1 (Figura 6). Además,
los datos obtenidos de las pruebas de experimentación se piezoeléctrica consistente.
pueden acceder mediante el siguiente enlace:
𝐼 = 𝐹∆𝑡 (2)
https://thingspeak.com/channels/2285032.
𝐼 = 𝑚𝑔∆𝑡 (3)
El Algoritmo 1 describe el procedimiento lógico secuencial 𝑚
para realizar la adquisición de datos provenientes desde la 𝐼 = 70 𝑘𝑔 ∗ 9,8 2 ∗ 1𝑠 (4)
𝑠
alfombra piezoeléctrica. Una vez obtenidos los datos, se
realiza la conexión de la placa de desarrollo NODE MCU 𝐼 = 686 𝑘𝑔 ∗ 𝑚/𝑠 (5)
ESP 8266 mediante la ID y la API KEY del canal de
ThingSpeak. De manera posterior, mediante el sensor INA La Ecuación 2 indica que el impulso (I) es igual a la fuerza
2019 se recopila el voltaje y la corriente. En este sentido, el (F) por el cambio en el tiempo (Δt). Es decir, la corriente
voltaje se multiplica por un factor de ajuste para obtener una generada es proporcional a la fuerza aplicada y al tiempo en
buena precisión. Este factor viene dado por el fabricante del que se aplica dicha fuerza. Por otro lado, la Ecuación 3
sensor. Los datos se recopilan cada 2 ms y se envían a la especifica que la fuerza (F) es equivalente a la masa (m) por
plataforma IoT de Thing Speak. la aceleración gravitacional (g). Por lo tanto, la corriente es
proporcional a la masa que ejerce la fuerza y al intervalo de rango de frecuencias de operación abarcaba desde 400 Hz
tiempo. hasta 9 kHz, lo que significa que tenía una alta sensibilidad a
las frecuencias elevadas o cuando se aplicaba la misma carga.
Los datos característicos indicaron que este piezoeléctrico era
más reactivo a las altas frecuencias o cuando se aplicaba el
mismo peso de carga. Incluso una ligera inclinación de la
presión resultaba en un corto período de presión y una mayor
amplitud de la tensión de impulso. La medición de las
características de la tensión de salida también podía utilizarse
para determinar las especificaciones del convertidor de
potencia.
Figura 7. Piezoeléctricos dispuestos sobre azulejo (Izquierda). Vista Los resultados de las mediciones del sistema de captación de
superior de alfombra piezoeléctrica (Derecha) energía muestran que la tensión a través del condensador o
La respuesta del impulso piezoeléctrico, según la Ecuación 2, batería aumenta con la adición de cada elemento, aunque el
se ve influenciada por la magnitud de la fuerza aplicada y la aumento no es lineal absoluto. Además, cada elemento
duración del intervalo de tiempo durante el cual se ejerce la adicional aumentaba la tensión del condensador en 0,1 V.
fuerza. En caso de que la fuerza sea mayor y el intervalo de Los resultados más altos se obtuvieron con una configuración
tiempo sea más extenso, el impulso generado será más en paralelo de 25 elementos conectados, donde la tensión del
significativo. Esto se ejemplifica en la Figura 8, que condensador alcanzó 0,57 V.
representa el ángulo de caída del peso sobre la superficie de
la alfombra. En este contexto, la Figura 8(a) muestra una
superficie plana, mientras que la Figura 8(b) representa una
superficie inclinada.
𝑝 = 𝐹/𝐴 (6)
Figura 9. Comportamiento de carga y nivel de voltaje de la batería en base a
distintos arreglos de sensores piezoeléctricos
De acuerdo con la Ecuación (6), las variaciones en el ángulo
o en la configuración superficial de un objeto en relación con La Figura 9 muestra las diferencias de cálculo para diferente
la fuerza ejercida tendrán un impacto en la presión que se número de sensores piezoeléctricos conectados en serie y
transmite a la superficie de otro objeto. La relación entre la paralelo. En este marco, el arreglo 1 consta de un solo sensor
fuerza (F) con el área de la sección transversal o superficie piezoeléctrico. El arreglo 2 consta de tres sensores
(A) es inversamente proporcional. Cuanto mayor sea la piezoeléctricos conectados en serie y paralelo, de forma
sección transversal o la superficie afectada por la presión, mixta. El arreglo 3 consiste en 6 elementos piezoeléctricos
menor será la fuerza generada. Por el contrario, cuanto más conectados de forma mixta. Por último, el arreglo 4
estrecha sea la sección transversal, mayor será la fuerza corresponde a 25 sensores piezoeléctricos conectados en una
generada. Por lo tanto, existe una correlación entre el ángulo configuración mixta. Según el cálculo, los 25 elementos
de caída del peso y la respuesta de la tensión piezoeléctrica. deberían haber producido una tensión de 1 V con un aumento
Cuando el peso cae precisamente en un ángulo de 90 grados de 0,1 V por cada elemento adicional. Sin embargo, en
perpendiculares a la superficie de la alfombra, como se condiciones reales, la tensión de la batería solo pudo alcanzar
muestra en la Figura 8(a), el periodo es más prolongado que 0,57 V.
cuando las pesas caen y forman un ángulo, como se muestra
en la Figura 8(b). Cuanto más largo es el periodo, más baja es La prueba de carga también puede ser calculada a partir de la
la frecuencia; cuanto más corto es el periodo, más alta es la información sobre el condensador o la batería, la energía que
frecuencia. se convierte y la energía generada por la alfombra
piezoeléctrica. La energía convertida sigue la relación
El piezoeléctrico empleado en la alfombra piezoeléctrica fue definida en la Ecuación 7, la cual está vinculada a la energía
un zumbador piezoeléctrico con una dimensión de 35 mm. Su
Resumen: Al reemplazar la generación síncrona con generación basada en convertidores, se reduce la inercia
rotacional disponible en el sistema eléctrico. Esta disminución de inercia modifica el comportamiento dinámico del
sistema, lo que provoca cambios más pronunciados en la frecuencia del mismo y aumenta la posibilidad de activar
esquemas de disparo automático de carga (DAC). Con el fin de evitar la disminución de la inercia rotacional, en la
literatura se han elaborado esquemas de control para los convertidores que emulan el comportamiento de las máquinas
síncronas, conocidos como controles grid-forming. Este artículo presenta una metodología general que permite
optimizar simultáneamente los parámetros de los esquemas DAC y convertidores grid-forming en sistemas de baja
inercia rotacional. El objetivo es minimizar la carga disparada durante eventos de baja frecuencia. Se presentan dos
aportaciones: la primera es el procedimiento para evaluar la función objetivo, que consiste en evaluarla mediante el
uso de resultados de una simulación en el dominio del tiempo; la segunda es la integración de parámetros de control
grid-forming al problema de optimización, con el fin de aprovechar la respuesta rápida de los convertidores y
proporcionar una respuesta más eficiente ante eventos de baja frecuencia. La metodología ha sido validada mediante
diversos casos de estudio desarrollados usando el sistema de potencia tradicional de 9 barras, presentado
originalmente en el libro de P.M. Anderson. Los resultados obtenidos corroboran que la metodología propuesta reduce
la activación del esquema DAC durante un evento de baja frecuencia.
*[email protected]
Recibido: 15/08/2023
Aceptado:19/07/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.07
CC 4.0
La resistencia R y la reactancia X de la impedancia en serie Z pendiente establecida por una ganancia de regulación primaria
se calcula mediante las siguientes ecuaciones: (Kundur, 2017). Este es un control que ha demostrado ser útil
para repartir la potencia activa entre múltiples unidades de
Control grid- generación considerando una desviación de la frecuencia. De
forming manera similar, los ingenieros han tomado las ventajas
I aportadas por la simplicidad de este control para efectuar el
control de la potencia reactiva y el voltaje. A continuación, se
muestran las ecuaciones usadas para modelar el lazo de control
R jX
de potencia activa control droop:
VT θT 𝑑𝑃𝑎𝑐𝑡 𝑃𝑚𝑒𝑎 − 𝑃𝑎𝑐𝑡
Vm θm = (4)
𝑑𝑡 𝑇𝑝
∆𝜔𝑟 = (𝑃𝑠𝑒𝑡 − 𝑃𝑎𝑐𝑡 )𝑚𝑝 (5)
Modelo del convertidor 𝑑𝜃𝑚
= ∆𝜔𝑟 𝜔𝑏 (6)
Figura 1. Circuito que representa el modelo RMS de un convertidor del tipo 𝑑𝑡
VSC mostrando el lazo de control para el tipo grid-forming donde Pset es el punto de consigna de potencia en por unidad
(pu), Pmea es la potencia activa medida en pu, Pact es la potencia
𝑢𝑘 activa medida con un retraso de tiempo Tp, mp es el droop de
𝑍= (1)
100 potencia activa, ∆ωr es la desviación de la velocidad angular
𝑃𝑐𝑢 en pu, ωb es velocidad angular base en rads/s, y θm es el ángulo
𝑅= (2)
1000 ∗ 𝑆𝐵 del vector de tensión de secuencia positiva.
𝑋 = √𝑍 2 − 𝑅2 (3) A continuación, se muestran las ecuaciones correspondientes
donde uk es la impedancia de cortocircuito en %, el cual se a las ecuaciones del lazo de control de potencia reactiva:
puede derivar de las impedancias en de los filtros del 𝑑𝑄𝑎𝑐𝑡 𝑄𝑚𝑒𝑎 − 𝑄𝑎𝑐𝑡
convertidor y del transformador de interconexión. El valor Pcu = (7)
𝑑𝑡 𝑇𝑞
representa las pérdidas en el cobre en kW en los elementos del ∆𝑣 = (𝑄𝑠𝑒𝑡 − 𝑄𝑎𝑐𝑡 )𝑚𝑞 (8)
convertidor, y SB es la potencia base del sistema en MVA. 𝑉𝑚 = 𝑉𝑠𝑒𝑡 + ∆𝑣 (9)
donde Qset es el punto de ajuste de potencia reactiva en pu, Qmea
2.2 Controles grid-forming para convertidores es la potencia reactiva medida en pu, Qact es la potencia activa
medida con un retraso de tiempo Tq, mq es el droop de potencia
Los modelos de control denominados grid-forming se han reactiva, ∆v es la desviación de voltaje, Vset es el punto de
planteado en la literatura científica como una opción válida ajuste de voltaje y Vm es la magnitud del voltaje de secuencia
para superar algunos de los retos de los sistemas de potencia positiva vector.
con alta penetración de convertidores y baja participación de
generadores sincrónicos, con un gran potencial para mejorar la 2.2.2 Máquina síncrona virtual (VSM)
respuesta de los convertidores ante disturbios del sistema
(Habibullah et al., 2021). La técnica de control denominada la máquina virtual síncrona
El control del tipo grid-forming se puede realizar usando lazos (virtual synchronous machine, VSM) está basada en el intento
de control de frecuencia y voltaje de los VSC. En la literatura de emular el comportamiento de una máquina síncrona,
científica moderna, se pueden encontrar múltiples propuestas modelando la dinámica electromecánica en una versión
para este tipo de control. En este trabajo, se ha decidido reducida de la ecuación de oscilación de la máquina sincrónica
considerar los modelos de control de tipo grid-forming para (Zhong, 2016). A continuación, se muestran las ecuaciones
VSC, los cuales se consideran como soluciones dominantes en usadas el control del tipo VSM:
los futuros sistemas de potencia: 𝑑∆𝜔𝑟 𝑃𝑠𝑒𝑡 − 𝑃𝑎𝑐𝑡 − 𝐷𝑝 ∆𝜔𝑟
Control droop (Rosso et al., 2021) = (10)
𝑑𝑡 𝑇𝑎
Máquina Síncrona Virtual (VSM) (Musca et al., 2022) 𝑑𝜃𝑚
Synchronverter (Sync) (Deepak et al., 2019; Zhong & = ∆𝜔𝑟 𝜔𝑏 (11)
𝑑𝑡
Weiss, 2011) La constante de tiempo Ta es la denominada constante de
Oscilador virtual despachable (dVOC) (Ajala et al., 2021; aceleración en la máquina sincrónica, y es conocido por la
Lu et al., 2021) teoría de máquinas eléctricas que la constante de aceleración
En las siguientes sub-secciones, se da una breve descripción está relacionada con la denominada constante de inercia
de estos modelos de grid-forming, que incluye detalles de las rotacional de la máquina H, (Ta = 2H). En algunos documentos
ecuaciones necesarias para el modelado. científicos, se considera a H como la inercia del generador
síncrono y puede denominarse constante de inercia virtual.
2.2.1 Control de droop La ecuación de oscilación expandida incluye un término
relacionado con el amortiguamiento de desviaciones de la
La base del control droop es emular el comportamiento del velocidad. El amortiguamiento Dp se usa generalmente en el
control primario de frecuencia que se utiliza en generador modelo clásico de la máquina síncrona y representa el cambio
síncrono tradicional. El principio de control primario es que de carga proporcional al cambio de frecuencia. En el modelo
cualquier cambio en la potencia del sistema produce un cambio VSM, la amortiguación Dp se puede emplear de manera similar
en la frecuencia, este cambio generalmente está dado por una al droop de frecuencia.
Algunos autores proponen que el control de voltaje del modelo (Dhople et al., 2013). Este control controla el VSC para emular
VSM sea equivalente a un modelo de excitador proporcional la dinámica de amortiguamiento no lineal, como el oscilador
simple. A continuación, se muestran las ecuaciones utilizadas Van der Pol (Lu et al., 2019).
para representar el control de voltaje de VSM: Estos osciladores pueden generar oscilaciones periódicas,
𝑑𝑉𝑎𝑐𝑡 𝑉𝑚𝑒𝑎 − 𝑉𝑎𝑐𝑡 autosostenidas y estables, y cuando se aprovechan como
= (12) controladores para convertidores, estos ofrecen sincronización
𝑑𝑡 𝑇𝑣
𝑉𝑥 = (𝑉𝑟𝑒𝑓 − 𝑉𝑎𝑐𝑡 )𝐾𝑣 (13) sin comunicación y uso compartido de la potencia, así como
𝑉𝑚𝑎𝑥 , 𝑉𝑥 ≥ 𝑉𝑚𝑎𝑥 regulación de voltaje y frecuencia. Sin embargo, es importante
𝑉𝑚′ = { 𝑉𝑥 , 𝑉𝑚𝑖𝑛 < 𝑉𝑥 < 𝑉𝑚𝑎𝑥 (14) mencionar que el VOC no es capaz de controlar la potencia
𝑉𝑚𝑖𝑛 , 𝑉𝑥 ≤ 𝑉𝑚𝑖𝑛 activa y reactiva inyectada por el VSC a la red, por lo que
𝑑𝑉𝑚 𝑉𝑚′ − 𝑉𝑚 recientemente, se ha desarrollado el denominado VOC
= (15) despachable (dVOC) (Ajala et al., 2021).
𝑑𝑡 𝑇𝑙𝑝𝑓
El controlador está compuesto por un circuito tanque LC
donde Vmea es el voltaje medido en terminales del convertidor, resonante y fuentes de corriente y voltaje no lineales
Vact es el voltaje medido con un retraso de tiempo Tv, Kv es la dependientes del estado. Los autores en Lu et al. (2021)
ganancia del control de voltaje, Vx es una variable auxiliar derivan un modelo dinámico para controlar el voltaje y la
limitada por Vmin y Vmax que son el valor máximo y mínimo del frecuencia del modelo dVOC. Estas ecuaciones se presentan a
voltaje de salida V’m, Vm es el voltaje de salida retrasado por continuación:
Tlpf que es la constante de tiempo de respuesta del control. 𝑑𝑉𝑚 𝜉 𝑘𝑣 𝑘𝑖
2
= 2 𝑉(2𝑉𝑛𝑜𝑚 − 2𝑉 2 ) − (𝑄 − 𝑄∗ ) (26)
𝑑𝑡 𝑘𝑣 3𝐶𝑉
2.2.3 Synchronverter 𝑑𝜃𝑚 𝑘𝑣 𝑘𝑖
= 𝜔𝑛𝑜𝑚 − (𝑃 − 𝑃∗ ) (27)
El control, denominado Synchronverter (SynC) (Zhong & 𝑑𝑡 3𝐶𝑉 2
Weiss, 2011), es considerado una técnica de control para
Se elige los factores de escala kv y ki como:
habilitar a los VSC cómo grid-forming. Este control es similar
al modelo VSM, ya que busca emular el comportamiento de
𝑘𝑣 = 𝑉𝑛𝑜𝑚 (28)
un generador síncrono. El controlador Synchronverter también 3𝑉𝑛𝑜𝑚
utiliza las ecuaciones de oscilación para controlar la 𝑘𝑖 = (29)
frecuencia; sin embargo, en lugar de la potencia activa, el 𝑆𝑟𝑎𝑡𝑒𝑑
SynC utiliza los pares Tset y Te. El control SynC incluye una
caída de control de potencia reactiva utilizada para regular el Al establecer la derivada 𝑑𝜃𝑚 /𝑑𝑡 = 𝜔, el modelo por unidad
voltaje a través de la excitación de la máquina virtual. A del dVOC queda definido como:
continuación, se muestran las ecuaciones pertinentes al SynC:
𝑃𝑠𝑒𝑡 𝑑𝑉𝑚 2)
1 ∗ )
𝑇𝑠𝑒𝑡 = (16) = 2𝜉𝑉𝑝𝑢 (1 − 𝑉𝑝𝑢 − (𝑄𝑝𝑢 − 𝑄𝑝𝑢 (30)
𝜔𝑟 𝑑𝑡 √2𝐶𝑉𝑝𝑢
𝑑∆𝜔𝑟 𝑇𝑠𝑒𝑡 − 𝑇𝑒 − ∆𝜔𝑟 𝐷𝑝 1 ∗ )
= (17) 𝜔𝑝𝑢 = 𝜔𝑟𝑒𝑓 − (𝑃𝑝𝑢 − 𝑃𝑝𝑢 (31)
𝑑𝑡 𝑇𝑎 2
√2𝐶𝜔𝑏 𝑉𝑝𝑢
𝜔𝑟 = 𝜔𝑟𝑒𝑓 + ∆𝜔𝑟 (18)
𝑑𝜃𝑟 Sin embargo, al inspeccionar la ecuación para el control de
= ∆𝜔𝑟 𝜔𝑏 (19)
𝑑𝑡 voltaje, se puede ver que el voltaje en estado estable solo puede
𝜋
𝜃𝑚 = 𝜃𝑟 + (20) ser V = 1.00 pu debido a la relación cuadrática. Para solucionar
2 esto, el voltaje base debe configurarse diferente al voltaje
𝑄𝑒𝑟𝑟 = 𝑄𝑠𝑒𝑡 − (𝑉𝑚𝑒𝑎 − 𝑉𝑟𝑒𝑓 )𝐷𝑞 − 𝑄𝑒 (21) nominal, Vb ≠ Vnom. A partir de esto, las ecuaciones dinámicas
𝑑𝑀𝑓 𝑖𝑓 𝑄𝑒𝑟𝑟 del dVOC son las siguientes:
= (22)
𝑑𝑡 𝐾𝑣
𝑇𝑒 = 𝑀𝑓 𝑖𝑓 ∗ 𝐼𝑎𝑏𝑠 sin(𝜃𝑟 − 𝜃𝑖 ) (23) 𝑑𝑉𝑚 2 2
𝑉𝑝𝑢 𝑉𝑏 2
𝑉𝑛𝑜𝑚 ∗ )
𝑉𝑚 = 𝜔𝑟 𝑀𝑓 𝑖𝑓 (24) = 2𝜉𝑉𝑝𝑢 (1 − 2 ) − (𝑄𝑝𝑢 − 𝑄𝑝𝑢 (32)
𝑑𝑡 𝑉𝑛𝑜𝑚 √2𝐶𝑉𝑝𝑢 𝑉𝑏2
𝑄𝑒 = −𝑉𝑚 𝐼𝑎𝑏𝑠 cos(𝜃𝑟 − 𝜃𝑖 ) (25) 2
𝑉𝑛𝑜𝑚 ∗ )
𝜔𝑝𝑢 = 𝜔𝑟𝑒𝑓 − (𝑃𝑝𝑢 − 𝑃𝑝𝑢 (33)
2 𝑉2
donde Pset es el punto de consigna de potencia en pu, Te es el √2𝐶𝜔𝑏 𝑉𝑝𝑢 𝑏
par eléctrico, Qe es la potencia reactiva, Iabs es la magnitud del
fasor de corriente medido en el convertidor, θi es el ángulo del 3. METODOLOGÍA DE OPTIMIZACIÓN DE
fasor de corriente, Mfif es la excitación simplificada de la PARÁMETROS
máquina virtual, θr el ángulo del rotor virtual que retrasa el
ángulo del fasor de voltaje por θm por π/2. 3.1 Métricas de frecuencia
2.2.4 Oscilador virtual despachable (dVOC) La respuesta de la frecuencia estable de un sistema de potencia
tradicional siguiente a un evento de pérdida de generación (o
El modelo de control denominado dispatchable virtual incremento de la demanda) se muestra en la Figura 2. La
oscillator (dVOC) es un control del tipo grid-forming para respuesta en frecuencia puede ser caracterizada por varios
VSC inspirado en el control de oscilador virtual (VOC) indicadores o métricas que definen la calidad de la respuesta.
Las métricas más comúnmente utilizadas son: inspirada en el trabajo presentado originalmente en Acosta et
𝒅𝒇 al. (2020) y Sanchez Gorostiza et al. (2020).
(1) La derivada 𝒅𝒕 de la frecuencia a lo largo del tiempo
El objetivo del problema de optimización es reducir el total de
conocida como rate of change of frequency (RoCoF),
(2) La frecuencia mínima fmin durante el disturbio, y el tiempo la desconexión de carga en el sistema (Ptot), manteniendo la
frecuencia de estado estable (fss) posterior al evento, dentro de
en el que ocurre la frecuencia minina.
los límites operativos. En este estudio, se considera un
(3) La frecuencia fss a la que se estabiliza el sistema.
esquema DAC que desconecta la carga por bloques
escalonados. Dichos esquemas, actúan desconectando un valor
de carga Pshed,i, utilizando una lógica de disparo cuando la
frecuencia cae por debajo de un valor ajustado (f<fthr,i). La
lógica de disparo va a actuar con un retraso de tiempo Td,i, que
puede ser diferente para cada etapa del esquema DAC.
La optimización de la configuración de un esquema DAC
permite tener una menor desconexión de carga ante grandes
disturbios que causen un evento de baja frecuencia. Partiendo
Figura 2. Típico evento de baja frecuencia estable, indicando tres métricas de esta idea, la función objetivo f(x) del problema de
comúnmente usadas. fmin, fss, df/dt
optimización se puede definir cómo el total de la desconexión
3.2 Planteamiento general de un problema de optimización de carga del sistema durante un evento de baja frecuencia:
𝑛𝑟 𝑛𝑠
En esta investigación, se busca minimizar la desconexión de
carga de un sistema de potencia durante un evento de baja 𝑃𝑡𝑜𝑡 = 𝑓(𝐱) = ∑ ∑ 𝑃𝑠ℎ𝑒𝑑𝑖,𝑗 (35)
frecuencia. Para lograr eso, se toma ventaja de la respuesta 𝑗=1 𝑖=1
rápida y capacidad de regulación de frecuencia que ofrecen los
convertidores electrónicos de potencia con control de tipo donde Pshed,ij es el deslastre de la etapa i del esquema DAC j,
grid-forming. En tal sentido, es necesario formular ns es el número de etapas activadas del esquema DAC, y nr el
apropiadamente el problema de investigación en forma de un número de esquemas de DAC.
problema de optimización matemática. El problema de optimización busca minimizar el deslastre de
Un problema de optimización puede ser visto como maximizar carga f(x), mediante la apropiada definición de las variables de
o minimizar el valor de una función objetivo f(x). Esto se logra control x.
calculando los valores del vector de variables de control x, para min[𝑓(𝐱)] (36)
𝐱
el cual la función objetivo evalúa al valor máximo o mínimo. Teniendo ya la función objetivo y el objetivo del problema de
Sin embargo, el problema típicamente incluye ciertas optimización, el siguiente paso es definir las variables de
restricciones a fin de incluir condiciones adicionales que deben control.
ser evaluadas. Las restricciones pueden ser formuladas Este estudio considera los límites operativos definidos por la
matemáticamente en varias formas, pueden ser límites norma PRC-006-SERC-02 (SERC, 2008). Siguiendo los
máximos y mínimos del vector x, o condiciones que se deben lineamientos de esta norma, la primera etapa del esquema
cumplir. DAC debe actuar dentro del rango de frecuencia 59,3 Hz
La forma general del problema de optimización se puede ≤f≤59,5 Hz. Para ajustar las etapas siguientes, estas deben
generalizar de la siguiente manera: tener actuar en un valor de fthr que tenga una diferencia
respecto al ajuste anterior de al menos ∆fthr = 0,2 Hz, pero no
min[𝑓(𝐱)] (34) más de ∆fthr = 0,5 Hz. En tal sentido, se tienen cuatro posibles
𝐱
sujeto a: variables de control para el problema de optimización que son
𝑐(𝐱) ≤ 0 Pshed, fthr, ∆fthr y Td. En este caso, se toman como variables de
𝑐𝑒𝑞(𝐱) = 0 control Pshed, fthr, y ∆fthr, la variable Td no se considera, ya que
𝐀𝐱 ≤ 𝐛 tiene poco efecto sobre la función objetivo. Entonces, el vector
𝑨𝒆𝒒𝐱 = 𝒃𝒆𝒒 de variables de control es definido de la siguiente manera:
𝑙𝑏 ≤ 𝒙 ≤ 𝑢𝑏
donde b y beq son vectores, A y Aeq son matrices, c(x) y ceq 𝐱 = [𝐏𝐬𝐡𝐞𝐝 𝐟𝐭𝐡𝐫 ∆𝐟𝐭𝐡𝐫 ]𝑇 (37)
(x) son funciones no lineales que devuelven vectores, por 𝐏𝐬𝐡𝐞𝐝 = [𝑃𝑠ℎ𝑒𝑑1,1 … … … 𝑃𝑠ℎ𝑒𝑑𝑁𝑟,𝑛𝑟 ] (38)
último, lb y ub son los límites inferiores y superiores de las 𝐟𝐭𝐡𝐫 = [𝑓𝑡ℎ𝑟1,1 … … … 𝑓𝑡ℎ𝑟1,𝑛𝑟 ] (39)
variables que conforman al vector x. ∆𝐟𝐭𝐡𝐫 = [∆𝑓𝑡ℎ𝑟1,1 … … … ∆𝑓𝑡ℎ𝑟𝑁𝑟−1,𝑛𝑟 ] (40)
3.3 Metodología Propuesta 1: Optimización de los ajustes
Donde Pshed es un vector que contiene cada valor de deslastre
DAC
de carga, fthr es un vector con el primer valor de activación de
cada esquema DAC, ∆fthr es un vector de deltas de frecuencia
La hipótesis de esta metodología es que los ajustes de los
correspondientes a los valores de activación de las etapas del
esquemas de DAC pueden ser óptimamente ajustados de
DAC, Nr es el número de etapas de los DAC.
manera tal que el valor total de la carga a ser desconectada sea
Cómo siguiente paso, se utilizaron las recomendaciones de la
minimizado considerando las restricciones operativas del
norma PRC-006-SERC-02, para definir los límites superiores
sistema. La metodología propuesta en esta sección está
e inferiores, así como las funciones de desigualdad. La norma
también da indicaciones para los retrasos de tiempo Td, inercia rotacional. En esta metodología, no se toman en cuenta
mencionando que deben ser de al menos seis ciclos o 0,1s. los controles de voltaje de los convertidores grid-forming
Teniendo esto en cuenta, se definieron los límites superiores e debido a que se están considerando como desacoplados los
inferiores de la siguiente manera: comportamientos de la frecuencia y el voltaje.
Entonces, el vector de variables de control se define de la
59,3 𝐻𝑧 ≤ 𝑓𝑡ℎ𝑟,1 ≤ 59,5 𝐻𝑧 (40) siguiente manera:
𝑇
0,2 𝐻𝑧 ≤ ∆𝑓𝑡ℎ𝑟 ≤ 0,5 𝐻𝑧 (41) 𝐱 = [𝐏𝐬𝐡𝐞𝐝 𝐟𝐭𝐡𝐫 ∆𝐟𝐭𝐡𝐫 𝐑 𝐠𝐟𝐦 𝐓𝐚 ] (46)
0 ≤ 𝑃𝑠ℎ𝑒𝑑 ≤ 1 (42)
𝐑 𝐠𝐟𝐦 = [𝑅𝑔𝑓𝑚1 … … … 𝑅𝑔𝑓𝑚𝑛𝑔 ] (47)
Para finalizar se deben definir las desigualdades no lineales 𝐓𝐚 = [𝑇𝑎1 … … … 𝑇𝑎𝑛𝑔 ] (48)
que debe respetar el algoritmo de optimización. Estas se Dónde Rgfm y Ta son vectores conformados por las ganancias
definen considerando que la frecuencia de estado estable fss de regulación primaria y de inercia virtual de los convertidores
debe alcanzar un valor que se mantenga entre el rango de
grid-forming. Para definir los límites inferiores y superiores
operación fssmin ≤ fss ≤ fssmax:
del valor de Rgfm se ha considerado lo establecido en el código
de red mexicano (Energía, 2016), donde se recomiendan los
𝑐1(𝐱) = 𝑓𝑠𝑠 − 𝑓𝑠𝑠𝑚𝑎𝑥 ≤ 0 (43) siguientes valores:
𝑐2(𝐱) = 𝑓𝑠𝑠𝑚𝑖𝑛 − 𝑓𝑠𝑠 ≤ 0 (44)
3% ≤ 𝐑 𝐠𝐟𝐦 ≤ 8% (49)
Por último, se agrega una desigualdad no lineal extra, esta se
basa en que la norma PRC-006-SERC-02, no permite que el Con esta modificación del vector de variables de control (x), y
parámetro de la etapa final de un esquema DAC sea menor a agregando los respectivos límites, se puede considerar este
58,4 Hz.
problema de optimización como una extensión del problema
de optimización de parámetros del esquema DAC.
𝑐3(𝐱) = 𝑓𝑡ℎ𝑟,𝑁𝑠 ≤ 58,4 𝐻𝑧 (45)
3.5 Implementación de la metodología optimización
Dónde Ns es el número de etapas del relevador.
La metodología propuesta 1 tiene la capacidad de optimizar los Para formular y resolver el problema de optimización expuesto
ajustes del esquema DAC para minimizar el deslastre de carga, en las metodologías anteriormente propuestas, se ha
pero no toma en consideración el posible uso de los implementado un esquema centrado en el uso del lenguaje de
convertidores equipados con control del tipo grid-forming. En programación de MATLAB® El uso de MATLAB® se
la siguiente subsección se presenta una metodología donde los fundamenta, debido a que un objetivo a futuro es probar y
parámetros del control grid-forming son incluidos en el demostrar las metodologías propuesta en tiempo real, usando
problema de optimización. la plataforma de modelación y simulación RT-LAB® de
OPAL-RT®; la cual utiliza MATLAB® y SIMULINK®.
3.4 Metodología Propuesta 2: Optimización de los ajustes En la Figura 3, se muestra un diagrama general de la
UFLS y parámetros de control grid-forming implementación para resolver el problema de optimización. La
interfaz utiliza diferente software de modelado y simulación,
En esta sección, los autores proponen una metodología permitiendo que a futuro pueda ser expandida fácilmente. Para
original para incluir en el problema de optimización modelar los convertidores grid-forming, se ha utilizado el
simultanea de los parámetros del esquema DAC y el control software OpenModelica®, el cual cuenta con la habilidad de
del tipo convertidores grid-forming, con el objetivo de la exportar los modelos desarrollados en formato ‘fmu’.
reducción de deslastre de carga en eventos de baja frecuencia. Posteriormente, estos modelos pueden ser integrados de
En esta propuesta la función objetivo es nuevamente definida manera sencilla al simulador de transitorios electromecánicos
como el total del deslastre de carga del sistema que fue ePHASORSIM® y así ser resueltos en régimen dinámico. A su
definido previamente en la Ecuación (35) y se busca minimizar vez, ePHASORSIM® es una librería para el software
la función f(x) según lo establecido en la Ecuación (36). SIMULINK® que está incluida en la plataforma de simulación
Las variables de control empleadas en la optimización deben RT-LAB®. Al operar dentro de SIMULINK®, se pueden
incluir los parámetros de interés del esquema DAC, según la también utilizar bloques de SIMULINK® para hacer control de
Ecuación (37) y del control grid-forming. A manera de elementos del sistema de potencia. En este caso, el esquema
ejemplo se definen las variables controlables del modelo VSM DAC ha sido modelado dentro de SIMULINK®, y para
presentado en la Sección 2,2,2. calcular la frecuencia de los buses del sistema, se está
El control VSM tiene como parámetros: la constante de inercia utilizando un bloque de SIMULINK® que utiliza la derivada
virtual Ta, la constante de amortiguamiento virtual Dp, y la del ángulo de fase para estimar la frecuencia de cada bus. Esta
constante de regulación de voltaje Kv. En esta investigación, se señal de frecuencia es utilizada para evaluar los disparos de los
eligió utilizar las constantes Dp y Ta como variables de control. esquemas de DAC.
La razón tras esta decisión, es que se puede sintonizar el valor
de Dp al valor de las constantes de regulación primaria si se
define como Dp = 1/Rgfm. La constante de aceleración Ta del
control VSM se puede incluir en el planteamiento del
problema de optimización, especialmente para mejorar la
respuesta en frecuencia del sistema en escenarios de baja
Figura 3. Esquema general de la implementación del algoritmo de optimización. Se muestra el enlace entre las diferentes herramientas de software
utilizadas
Inicialmente, se valida la metodología propuesta usando el Escenario 1: Optimización de esquemas DAC en un sistema
sistema WSCC 9-barras. En este studio, se utiliza una versión con generación síncrona. Se optimizan los esquemas DAC
modificada de la originalmente propuesta por el Prof. Fouad, incluyendo solo la dinámica electromecánica de las máquinas.
las modificaciones consisten en utilizar el modelo de máquina La idea de este escenario es verificar la validez de la
síncrona llamado GENROU para modelar los generadores; metodología y la implementación, demostrar que el concepto
dicho modelo es equivalente al Modelo 2,2 en Machinery funciona. Este escenario tiene como objetivo mostrar la
Committee of the IEEE Power (2019), cuando se consideran evolución de la frecuencia en el tiempo al ocurrir un disturbio
las reactancias subtransitorias xd'' y xq'' tienen el mismo valor. en la red. La simulación consiste en la desconexión repentina
Se agregan controles a los generadores, en este caso se utilizan de la máquina G2 en el tiempo t = 0,50s, la simulación termina
EXST1 y TGOV1 para el excitador y el gobernador, en t=20,00s y el paso de tiempo es de ∆t=0,001s. Esto se
traduce en una pérdida de generación P = 163,00 MW, y Q = de tiempo es de ∆t=0,001s. A esto se le agrega que el G3 es
6,65 Mvar. La reserva de potencia activa entre las máquinas cambiado por un convertidor grid-forming llamado GFM3.
G1 y G3 es de P = 293,53 MW. En este caso, se realiza una Se aplicarán diferentes controles grid-forming a GFM3 con el
optimización de la configuración del esquema de DAC para objetivo de analizar cómo difiere la solución óptima al incluir
sistema de 9 barras operando con máquinas síncronas. Este diferentes dinámicas de control. Los controles grid-forming se
sistema de prueba cuenta con tres cargas, las cargas se ubican sintonizaron para que se tenga una respuesta similar al
en los buses 5, 6, y 8. Cada carga contará con un esquema DAC generador síncrono G3. Las condiciones iniciales del esquema
a ser optimizado, y se llamarán L5, L6, y L8, respectivamente. DAC son las mismas que en la Tabla 1.
Para este caso de optimización, las condiciones iniciales de los En la Tabla 3, se muestran los resultados de las optimizaciones
esquemas DAC se muestran en la Tabla 1. antes mencionadas. En algunos casos, los resultados de la
metodología de optimización se ajustan parámetros de la
Tabla 1. Condiciones iniciales del esquema DAC tercera etapa de los esquemas de DAC. Esto debido a que
fthr [Hz] Td [s] Pshed [pu] durante las primeras iteraciones se llegaba a superar la
frecuencia de operación de la tercera etapa, pero en los casos
L5 59,45 0,15 0,30
finales no es necesaria su activación. Esto debido a que se
L6 59,20 0,15 0,15 están optimizando los parámetros de la frecuencia de
activación de las etapas en conjunto con la desconexión de
L8 58,95 0,15 0,05 carga por etapa.
Se puede observar que la desconexión de carga final es similar
En la Figura 4, se muestra la evolución del algoritmo de al Escenario 1. Esto valida dos puntos: 1) La metodología es
optimización en cada iteración. El deslastre de carga final es aplicable a casos con máquinas síncronas y convertidores grid-
de Pshed = 1,26495 pu de la carga, lo que representa un forming. 2) Los convertidores grid-forming son capaces de
deslastre de 126,495 MW de potencia activa, para mantener la emular satisfactoriamente el comportamiento de una máquina
frecuencia de estado estacionario dentro de los límites síncrona. La única excepción a este caso es el control dVOC,
preestablecidos de 59,5 Hz < fss < 60,5 Hz, y se respeta los el cual está diseñado para proveer una respuesta más rápida
límites de las variables de optimización presentadas en las
ecuaciones (40-42). Tabla 3. Resultados de optimización del caso 1, escenario 2. Se
muestran los parámetros óptimos de los DAC cuando el sistema opera
con un convertidor grid-forming y se pruebas diferentes técnicas de
control. Se resalta en gris las etapas no activadas en la solución final.
Pshd1 Pshed2 Pshed3 fthr1 fthr2 fthr3
VSM Total
[pu] [pu] [pu] [Hz] [Hz] [Hz]
L5 0,025 0,000 0,044 0,070 59,489 59,239 58,989
L6 0,400 0,269 0,044 0,669 59,328 59,079 58,829
L8 0,337 0,188 0,048 0,526 59,323 59,073 58,823
Total 1,264
Pshd1 Pshed2 Pshed3 fthr1 fthr2 fthr3
Droop Total
Figura 5. Curva de convergencia del algoritmo de optimización de DAC [pu] [pu] [pu] [Hz] [Hz] [Hz]
para el Caso 1 L5 0,149 0,042 0,063 0,191 59,434 59,184 58,934
L6 0,364 0,231 0,044 0,595 59,480 59,230 58,980
En la Tabla 2, se muestran los resultados finales de los
parámetros óptimos de los esquemas DAC. L8 0,332 0,155 0,048 0,487 59,485 59,235 58,985
Total 1,273
Tabla 2. Resultados optimización del Escenario 1. Sistema tradicional con Pshd1 Pshed2 Pshed3 fthr1 fthr2 fthr3
GS. SynC Total
[pu] [pu] [pu] [Hz] [Hz] [Hz]
Pshed1 Pshed2 Pshed3 Total fthr1 fthr2 fthr3
[pu] [pu] [pu] [pu] [Hz] [Hz] [Hz] L5 0,000 0,000 0,063 0,000 59,398 59,148 58,898
L5 0,218 0,000 0,000 0,219 59,369 59,119 58,869 L6 0,487 0,362 0,026 0,848 59,355 59,105 58,855
L6 0,289 0,184 0,106 0,581 59,300 59,050 58,800 L8 0,281 0,113 0,048 0,395 59,332 59,082 58,832
Tabla 4. Resultados de optimización del Caso 2. Se muestran los parámetros óptimos de los DAC cuando el sistema opera con un convertidor grid-
forming y se pruebas diferentes técnicas de control. Se resalta en gris las etapas no activadas en la solución final.
Mattsson, S. E., & Elmqvist, H. (1997). Modelica - An sistemas de potencia con penetración de convertidores.
International Effort to Design the Next Generation https://orcid.org/0009-0002-9563-1560
Modeling Language. IFAC Proceedings Volumes, 30(4),
151–155. https://doi.org/10.1016/s1474-
6670(17)43628-7 Francisco M., Gonzalez-
Musca, R., Gonzalez-Longatt, F., & Sánchez, C. A. G. (2022). Longatt, es fundador y líder de
Power System Oscillations with Different Prevalence of DIgEnSys-Lab (Laboratorio de
Grid-Following and Grid-Forming Converters. Sistemas de Energía Digital) en
Energies, 15(12), 1–19. el Departamento de Ingeniería
https://doi.org/10.3390/en15124273 Eléctrica, Tecnología de la
OPAL-RT Technologies. (2023). ePHASORSIM Información y Cibernética de la
Documentation. https://opal- Universidad del Sudeste de
rt.atlassian.net/wiki/spaces/PEUD/pages/144505576/Na Noruega, Noruega; Actualmente
tive+Library trabaja en el Centro de
Ortega, Á., & Milano, F. (2019). Combined Frequency and Tecnología de Sistemas de
RoCoF Control of Converter-Interfaced Energy Storage Energía Renovable (CREST) en Loughborough,
Systems. IFAC-PapersOnLine, 52(4), 240–245. Universidad en el Reino Unido. Su principal área de
https://doi.org/10.1016/j.ifacol.2019.08.198 interés es potenciar la operación segura y resiliente de los
PowerWorld Corporation. (n.d.). PowerWorld Block sistemas de energía aprovechando las ventajas de las
Diagrams. https://www.powerworld.com/files/Block- tecnologías digitales. https://orcid.org/0000-0002-7157-
Diagrams-18.pdf 9844
Rahmann, C., & Castillo, A. (2014). Fast frequency response
capability of photovoltaic power plants: The necessity of
new grid requirements and definitions. Energies, 7(10), José Miguel, Riquelme-
6306–6322. https://doi.org/10.3390/en7106306 Dominguez, nació en España en
Rosso, R., Wang, X., Liserre, M., Lu, X., & Engelken, S. 1993. Recibió el Grado en
(2021). Grid-Forming Converters: Control Approaches, Ingeniería Industrial y el Máster en
Grid-Synchronization, and Future Trends—A Review. Ingeniería Industrial por la
IEEE Open Journal of Industry Applications, 2(May), Universidad de Sevilla. Obtuvo el
93–109. https://doi.org/10.1109/ojia.2021.3074028 título de doctor en Ingeniería
Sanchez Gorostiza, F., Gonzalez‐Longatt, F., & Rueda, J. L. Eléctrica y Electrónica por la
(2020). Multi‐objective optimal provision of fast Universidad Politécnica de Madrid.
frequency response from EV clusters. IET Generation, Su experiencia laboral incluye el
Transmission & Distribution, 14(23), 5580–5587. operador del sistema eléctrico
https://doi.org/10.1049/iet-gtd.2020.0717 español, Red Eléctrica de España. Desde 2018 forma parte del
SERC. (2008). PRC-006-SERC-02. 1–13. Departamento de Ingeniería Eléctrica de la Universidad
Zhong, Q. C., & Weiss, G. (2011). Synchronverters: Inverters Politécnica de Madrid, donde actualmente ejerce de Profesor
That Mimic Synchronous Generators. IEEE Ayudante. Sus principales áreas de interés son: sistemas
Transactions on Industrial Electronics, 58(4), 1259– fotovoltaicos conectados a la red, y la estabilidad y el control
1267. https://doi.org/10.1109/TIE.2010.2048839 de sistemas eléctricos de potencia. https://orcid.org/0000-
Zhong, Q. C. (2016). Virtual Synchronous Machines: A 0001-9741-0388
unified interface for grid integration. IEEE Power
Electronics Magazine, 3(4), 18–27.
https://doi.org/10.1109/MPEL.2016.2614906 César, Angeles-Camacho,
recibió la licenciatura y maestría
BIOGRAFÍAS en ingeniería eléctrica del
Instituto Tecnológico de
José Martín, Valles Canales, Morelia, Morelia, Michoacán,
es un estudiante de Doctorado México, en 1992 y 2000,
en el Instituto de Ingeniería de respectivamente, y el doctorado
la UNAM. Su especialidad es en de la Universidad de Glasgow,
Ingeniería Eléctrica, Sistemas Glasgow, Escocia, Reino Unido,
Eléctricos de Potencia. Tiene en 2005. Es Profesor Asociado de
una Maestría en Ingeniería la Universidad Nacional
Eléctrica en Sistemas Eléctricos Autónoma de México, México. Es fundador y líder del
de potencia por la UNAM, y es Laboratorio de Investigación de Redes Eléctricas
Ingeniero en Energías Inteligentes en el Instituto de Ingeniería, UNAM. Sus
Renovables por el Tecnológico Nacional de México. intereses de investigación incluyen el desarrollo de
Actualmente está realizando una estancia en la aplicaciones para redes eléctricas inteligentes, modelado
Universidad de Loughborough en el Reino Unido y su y simulación de sistemas eléctricos de potencia.
trabajo se centra en la estabilidad de frecuencia en https://orcid.org/0000-0002-9499-2842
Abstract: Currently, the collection and crushing process of fine and coarse aggregate has become costly and
transportation time has increased due to the complexity of certain areas in the Peruvian highlands. As a consequence,
river aggregate is extracted as a total substitution of fine and coarse aggregate for the preparation of structural
concrete, but it is not known for sure if it is totally viable to use this aggregate for the benefit of construction interests.
Therefore, the objective of this study is to evaluate the total substitution of fine and coarse aggregate by river
aggregate obtained from the city of Cajamarca, Peru, on the mechanical properties of concrete. Concrete specimens
were manufactured with fine and coarse aggregate and others with river aggregate from the Huaquillo and Portachuelo
quarries, respectively. It is considered the designs of 175 kg/cm2 and 210 kg/cm2, in addition to tests such as slump,
unit weight, compressive strength, flexural strength and an analysis of variance (ANOVA) in Tukey's block. An
experimental campaign was carried out with cementitious pastes to evaluate the effect on physical and mechanical
properties due to the use of river aggregate; the samples containing fine aggregate and coarse aggregate revealed good
results in terms of mechanical properties. However, concrete made with river aggregate meets the minimum required
theoretical design strength and is suitable to be used as a total substitute, which is a new and important scientific issue
to highlight.
Keywords: coarse aggregate; concrete; fine aggregate; mechanical properties; river aggregate
Palabras clave: árido grueso; hormigón; árido fino; propiedades mecánicas; árido de río
*[email protected]
Recibido: 20/09/2024
Aceptado:23/06/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.08
CC 4.0
the different properties it can contain, such as resistance, sand in an optimum amount was 5 to 10% by mass has a
corrosion, durability, etc. (Flores Fernandez et al., 2019; favorable effect on the workability and strength of concrete.
Grinys et al., 2020), being the material that is used in large Likewise, the researchers, Limantara et al. (2017), mention in
quantities worldwide (Da Silva et al., 2017). Likewise, fine their findings on the compressive strength test that employing
and coarse aggregates form the skeleton of concrete and the natural river material, they obtained a value of 19.47 MPa
generally occupy 60-75 % of the volume (Yaragal et al., 2019). while with a crushed mixture, they obtained a value of 21.12
One of the constant and increasingly asked questions is MPa, respectively at 28 days of curing.
whether the deposits are suitable for preparing concrete
(Oliveira et al., 2020). Aggregates can greatly influence the Other research studies according to Eziefula et al. (2020)
fresh and hardened properties of concrete, in particular by their showed that the compressive strength increased with granite
shape and size characteristics (Güçlüer, 2020). The different by 38.13 N/mm2, river stone by 34.57 N/mm2 and local stone
mechanical properties mainly depend on the material attributes by 31.96 N/mm2 with the river aggregate being higher than
to make concrete strong and durable (Qing-Xiang et al., 2020). local stone. On the other hand, according to Laserna and
Such that, each professional controls the quality according to Montero (2016), they obtained a 15 % increase in compressive
his experience (Orozco et al., 2018). strength for natural river aggregates with a substitution ratio of
100 %, while a decrease of approximately 10-25 % was
However, in the construction industry (García et al., 2023), the observed for crushed natural aggregate mixtures. With respect
environmental issue should not only be seen by using natural to flexural strength, according to Hachani et al. (2017), they
ingredients, but also the recycling of materials from those mentioned that the flexural strength increased by 46 % in 28
found on riverbanks or beaches, the use of these materials can days compared to concrete with crushed aggregate. Finally, it
reduce costs (Fadzilla Sari et al., 2020). River sediments as was concluded that a concrete with pebbles will be better
aggregate material for concrete production have been employed in concretes with medium strength or non-structural
investigated by various researchers (Yan et al., 2022). In concrete.
addition, 30 % of the sediment volume is mainly sand and is
easily reusable, for concrete (Beddaa et al., 2020). Few studies have been carried out on the total substitution of
Nevertheless, in many cases, in order to use this material, it is fine and coarse aggregate from quarries by river aggregate in
necessary to transport it from long distances to the source of the concrete mix, presenting a current knowledge gap on the
origin, where prices increase considerably, one of the solutions certainty of its adequate use or not in structural concrete. This
is that the processed coarse aggregate is replaced by natural work contributes on the use of the total substitution of fine and
granular material (Zega et al., 2006). coarse aggregate by river aggregate, in order to reduce the
exploitation and overexploitation of aggregate quarries in the
This is of increasing concern as this material presents different highlands of Peru, and on the use that can be given to river
durability problems, as a result of the passage of time and the aggregate specifically depending on the characteristics of this
interaction with the environment (Solís-Carcaño and Alcocer- material. This generates the feasibility of its use for the
Fraga, 2019). Thus, the reuse of river sediments as aggregates construction of structural or non-structural elements to the
becomes a sustainable management technique since it offers surrounding communities.
the potential to obtain economic benefits by reducing up to 41
% the cost of concrete production (Beddaa et al., 2021). Recap, 2. MATERIALS AND EXPERIMENTAL PROCESS
that the use of natural stone aggregates is especially resilient,
as their procurement is net from river beds for concrete 3.1 Ordinary Portland Cement
production (Tugrul Tunc, 2018). Emphasizing a problem of
scarcity of natural aggregate materials, priority is given to the Ordinary Portland cement Type I (OPC) - Weight 42.50
use of river aggregate, which is a mixed material between kilograms was used to obtain a conventional hydraulic binder
aggregates and sands. In relation to this, the search for a new without any type of additive. The properties of cement are
source of fine aggregate substitute for the production of listed in Table 1.
concrete, which can be substituted naturally by river granular
material is important (Ararsa et al., 2018). Currently, in Peru Table 1. Properties of ordinary Portland cement
there are few scientific studies on the use of river aggregate as Reference
Testing Requeriments
standard
a total substitute for quarry aggregate, that is a gap existing Specific gravity 3.15 ASTM C150
until today. Initial stting 45 minutes ASTM C191
Final setting 375 minutes ASTM C191
While research demonstrates relevant data using river granular
material, the findings of Aïssoun et al. (2015) show that 3.2 Processed aggregate and river aggregate
concrete produced with crushed stone exhibited a 22-42 %
lower surface slump compared to similar concrete made with The aggregates used in the preparation of the concrete were:
river aggregate. Additionally, the increase in fines content fine aggregate and coarse aggregate from quarries. The quarry
from 8 to 18 % led to a significant reduction in slump due to aggregate was obtained from two locations, the first was the
the fine particles in the sand aggregate. Mixes made with Huaquillo quarry and the second the Portachuelo quarry,
crushed aggregate had 10 % higher compressive strength at 56 extracted from the city of San Ignacio, Cajamarca, Peru. The
days than those with river aggregate. However, other characteristics of the fine and coarse aggregates of the quarries
researchers, Duc-Trong et al. (2022), mentioned that crushed analyzed are shown in Table 2, the granulometric curves of the
The results do not vary with the quarry aggregates because the
river aggregate maintains a heterogeneous mixture, since it
was obtained from two different natural quarries, the first A
(Huaquillo) located at coordinates E 729 061.48 m; S 9 311
862.87 m, and the second B (Portachuelo) located at
coordinates E 727 478.96 m; S 9 435 647.70 m.
The concrete mix design was carried out under the ACI 211.1
standard, for a design strength of 175 and 210 kg/cm2 which
used fine aggregate, coarse aggregate and for the mix with
treatment the total replacement of 100 % of fine and coarse
aggregate by river aggregate considering each study quarry, a
water-cement ratio of 0.74 and 0.68 was used, being constant
for each design of 175 and 210 kg/cm2, respectively, no
Figure 2. Grain size distribution curve of quarry aggregate plasticizer was used, the proportions are shown in Table 4.
Figure 4. Compressive strength results Portachuelo Quarry (a and b); Huaquillo Quarry (c and d)
Figure 5. Flexural strength results Portachuelo Quarry (a and b); Huaquillo Quarry (c and d)
Figure 6(b), in reference to the variable design compressive Figure 6(c), in reference to the flexural strength variable of
strength 210, the p-value of significance of the ANOVA test design 175, the p-value of significance of the ANOVA test in
in blocks, presented a value less than 0.05 (p = 0.00213 < blocks, resulted to be less than 0.05 (p = 0.000903 < 0.05).
0.05), i.e., we can affirm that there is a significant difference That is to say, there is a significant difference between
between treatments T6 and T8, this affirmation can be treatments T5 and T7, also the Tukey Post Hoc test, allowed
supported with the Tukey multiple comparisons test, where us to support establishing different letters on the box plots,
both box plots show different letters, being treatment T6 the finding that treatment T5, was where the significantly higher
one that presented the significantly higher compressive flexural strength was recorded, finding a sample average
strength, reaching a sample mean compressive strength of flexural strength equal to 37.18 kg/cm².
243.33 kg/cm².
T5 T7
T6 T8
Treatment
Treatment
T5 T7
T6 T8
Treatment
Treatment
T5 T7
Treatment T6 T8
Treatment
Figure 6. Boxplot with Tukey's block ANOVA and multiple comparisons of the study: according to the variable’s compressive strength (kg/cm²), flexural
strength (kg/cm²) and slump test (in) for each corresponding design.
the ANOVA test in blocks presented a value greater than 0.05 2), 597 - 609. https://doi.org/10.1617/s11527-015-0522-
(p = 0.634 > 0.05), and the p-value of significance of the 9
ANOVA test in blocks presented a value greater than 0.05 (p Ararsa, W., Tucay Quezon, E., & Aboneh, A. (2018).
= 0.634 > 0.05), as well as equal letters are observed above the Suitability of Ambo Sandstone Fine Aggregate as an
box and whisker diagrams supporting the Tukey's multiple Alternative River Sand Replacement in Normal
comparisons test. Concrete Production. American Journal of Civil
Engineering and Architecture, 6(4), 140-146.
Finally, in Figure 6(f), no significant difference was visualized https://doi.org/10.12691/ajcea-6-4-2
between treatments T6 and T8 in reference to the variable Beddaa, H., Fraj, A. B., & Ducléroir, S. (2021). Experimental
settling test of design 210, where a p-value of significance study on river sediment incorporation in concrete as a
greater than 0.05 (p = 0.177 > 0.05) was found, as well as equal full aggregate replacement: Technical feasibility and
letters in the box and whiskers diagrams, which were assigned economic viability. Construction and Building
by Tukey's multiple comparisons test. Materials, 313, 125425.
https://doi.org/10.1016/j.conbuildmat.2021.125425
4. CONCLUSIONS Beddaa, H., Ouazi, I., Ben Fraj, A., Lavergne, F., & Torrenti,
J.-M. (2020). Reuse potential of dredged river sediments
The study shows that river aggregate significantly influences in concrete: Effect of sediment variability. Journal of
the mechanical properties of concrete by totally replacing the Cleaner Production, 265, 121665.
fine and coarse aggregate in the preparation of concrete for https://doi.org/10.1016/j.jclepro.2020.121665
structural purposes, which leads to the following conclusions Da Silva, M., Pepe, D. S., Mende de Andrade, R., Shubert
to be considered: Pfeil, M., & Toledo Filho, R. (2017). Rheological and
mechanical behavior of High Strength Steel Fiber-River
This study on the effect of river aggregate in total replacement Gravel Self Compacting Concrete. Construction and
of the processed aggregate (fine and coarse aggregate), Building Materials, 150, 606-618.
showed reliable results on the effect of river aggregate on the https://doi.org/10.1016/j.conbuildmat.2017.06.030
design strengths of 175 and 210 kg/cm², respectively. Duc-Trong, N., Duy-Liem, N., & MyNgoc-Tra, L. (2022). An
experimental investigation on the utilization of crushed
It was observed that with river aggregate from the Huaquillo sand in improving workability and mechanical
and Portachuelo quarries, the compressive strength was lower resistance of concrete. Construction and Building
between 3.7 % and 4.39 % and between 6.38 % and 7.62 %, Materials, 326, 126766.
corresponding to 175 kg/cm2 and 210 kg/cm2 respectively. https://doi.org/10.1016/j.conbuildmat.2022.126766
Eziefula, U. G., Opara, H. E., & Eziefula, B. I. (2020). Strength
It was observed that with river aggregate from the Huaquillo of concrete produced with different sources of
and Portachuelo quarries, the flexural strength was reduced aggregates from selected parts of Abia and Imo States of
between 8.47 % and 9.51 % and between 7.47% and 9.66 %, Nigeria. Journal of Engineering, Design and
corresponding to the strength of 175 kg/cm2 and 210 kg/cm2 Technology, 18(5), 1053 - 1061.
respectively. https://doi.org/10.1108/JEDT-08-2019-0220
Fadzilla Sari, F., Limantara, A. D., Ridwan, A., Gardjito, E.,
It is concluded that it is feasible to use river aggregate for the Subiyanto, B., Sudarmanto, H., . . . . Wiwoho
elaboration of structural concrete, totally substituting the Mudjanarko, S. (2020). Laboratory Testing on The
independent processed aggregates over the fine and coarse Standard Mixed Designed Paving with Bamboo
aggregates, since it presented values above the theoretical Material as Smooth and Rough Aggregate. IOP
design resistances. Conference Series: Earth and Environmental Science.
498, pp. 1-8. Jawa Timur, Indonesia: IOP Publishing
Consider that the results obtained include certain Ltd. https://doi.org/10.1088/1755-1315/498/1/012032
characteristics of aggregates belonging to the highlands of Flores Fernández, A. M., Villafranca Castillo, A. J., &
Peru, which is not an indicator that can be used in different Reconco Amaya, J. A. (2019). Concrete with recycled
areas of the country itself or abroad, as it should be analyzed aggregate: a construction material option with
in the respective tests per location. sustainability criteria. Unitec, 8, 120-124.
https://doi.org/10.5377/innovare.v8i2.9023
In addition, further studies on aggregates such as X-ray García, J., Arriola, G., Villena, L., Muñoz, S. (2023). Strength
disphraction, to determine the chemical and mineralogical of Concrete Using Partial Addition of Residual Wood
composition and its direct or indirect influence on the strength Ash with Respect to Cement. Revista Politecnica, 52(1).
of concrete, as well as chemical tests, are recommended. 45-54. https://doi.org/10.33333/rp.vol52n1.05
Grinys, A., Augonis, A., Daukšys, M., & Pupeiki, D. (2020).
REFERENCES Mechanical properties and durability of rubberized and
SBR latex modified rubberized concrete. Construction
Aïssoun, B. M., Soo-Duck, H., & Kamal, K. H. (2015). and Building Materials, 248, 118584.
Influence of aggregate characteristics on workability of https://doi.org/10.1016/j.conbuildmat.2020.118584
superworkable concrete. Materials and Structures, 49(1- Güçlüer, K. (2020). Investigation of the effects of aggregate
textural properties on compressive. Journal of Building
Engineering, 27, 9.
https://doi.org/10.1016/j.jobe.2019.100949
Hachani, M. I., Kriker, A., & Seghiri, M. (2017). Experimental BIOGRAPHIES
study and comparison between the use of natural and
artificial coarse aggregate in concrete mixture. Energy Paul Mauricio, Monje
Procedia, 119, 182-191. Yovera, Civil Engineer from
https://doi.org/10.1016/j.egypro.2017.07.067 Universidad Señor de Sipán,
Laserna, S., & Montero, J. (2016). Influence of natural with experience in the
aggregates typology on recycled concrete strength development of road
properties. Construction and Building Materials, 115, infrastructure projects,
78-86. reinforced concrete and others.
https://doi.org/10.1016/j.conbuildmat.2016.04.037 Studies of reinforced concrete
Limantara, A. D., Widodo , A., Winarto , S., Krisnawati , L. specialization programs with
D., & Mudjanarko , S. W. (2017). Optimizing the use of the use of specialized software
natural gravel Brantas river as normal concrete mixed and studies on the quality of
with quality fc = 19.3 Mpa. 4th International Conference concrete.
on Civil and Environmental Engineering for
Sustainability, IConCEES 2017. Langkawi.
https://doi.org/10.1088/1755-1315/140/1/012104 Socrates Pedro, Muñoz
Oliveira, M. A., Scotto, M. G., Barbosa, S., Freire de Andrade, Perez, Civil Engineer
C., & Freitas, M. d. (2020). Morphological controls and graduated from the
statistical modelling of boulder transport by extreme Universidad Nacional Pedro
storms. Marine Geology, 426, 16. Ruiz Gallo - Lambayeque in
https://doi.org/10.1016/j.margeo.2020.106216 2006, with 11 years of
Orozco, M., Avila, Y., Restrepo, S., & Parody, A. (2018). experience in execution of civil
Influencing factors in concrete quality: a survey of works and production of
relevant actors in the concrete industry. Revista precast and prestressed and 16
Ingeniería de Construcción, 33, 161-172. years of collegiate, with a
https://doi.org/10.4067/S0718-50732018000200161 Master's degree in Earth
Qing-Xiang, M., Dandan, L., & Yang, L. (2020). Mesoscale Sciences with mention in Geotechnics from the
computational modeling of concrete-like particle- Universidad Nacional de San Agustin - Arequipa, Doctor
reinforced composites with non-convex aggregates. in Civil Engineering from the Universidad Nacional del
Computers & Structures, 240, 106349. Santa. My main values are loyalty, responsibility,
https://doi.org/10.1016/j.compstruc.2020.106349 companionship, punctuality and the desire to excel.
Solís-Carcaño, R. G., & Alcocer-Fraga, M. A. (2019). Former undergraduate professor in the courses of Soil
Durability of concrete with high absorption aggregates. Mechanics, Fluid Mechanics, Geology at Universidad
Engineering, Research and Technology, 20, 1-13. Católica San Pablo, Universidad Nestor Caceres
https://doi.org/10.22201/fi.25940732e.2019.20n4.039 Velasquez, Universidad Alas Peruanas, Former
Tugrul Tunc, E. (2018). Strength Properties of Hardened Postgraduate Professor at Universidad Nestor Caceres
Concrete Produced with Natural Aggregates for Velasquez in the course of Soil and Rock Mechanics,
Different Water/Cement Ratios. European Journal of Professor at Universidad Nacional Toribio Rodriguez de
Science and Technology, 14, 280 - 287. Mendoza.
https://doi.org/10.31590/ejosat.486093
Yan, J., Zhong, S., Chen, S., Lv, Y., Yang, L., Peng, G., & Luigi Italo, Villena Zapata,
Deng, A. (2022). Study on the Application of Sediment- Statistical Engineer from the
Based Embankment Building and Ultra-High- National University of Trujillo,
Performance Concrete (UHPC) Preparation in the Master in Education with
Resource Utilization of Yellow River Sediment. mention in higher education
Materials, 15, 5668. didactics from the Universidad
https://doi.org/10.3390/ma15165668 Privada Antenor Orrego and
Yaragal, S. C., Basavana, G., & Rajasekaran, C. (2019). doctoral candidate in statistics
Characterization and performance of processed lateritic and mathematics at the
fine aggregates in cement mortars and concretes. Universidad Nacional del Santa.
Construction and Building Materials, 200, 10-25. As work experience I have been head of the Statistics Unit
https://doi.org/10.1016/j.conbuildmat.2018.12.072 of the Bagua Support Hospital "Gustavo Lanatta Lujan", I
Zega, C. J., Taus, V. L., & Di Maio, A. A. (2006). Physico- have worked in the Statistics and Informatics Unit and
mechanical behavior of recycled concretes made with Epidemiology Unit - Regional Health Management La
pebbles. Technical Bulletin, 44, 17-26. Retrieved Libertad. Teaching experience in universities and author
September 10, 2023, of books and articles in indexed journals Scopus and Web
http://ve.scielo.org/scielo.php?script=sci_arttext&pid= of Science.
S0376-723X2006000300003&lng=es&tlng=.
Resumen: Debido a la simplicidad y amplia aplicabilidad, la minimización irrestricta en la zona visible (Vis) es una
herramienta importante para resolver muchos problemas de optimización de parámetros y condiciones operativas de
sistemas fotovoltaicos. El propósito de este trabajo es usar la minimización irrestricta de una función objetivo
simplificada ξ para estimar la zona Vis. Utilizamos el método cuantitativo y técnica documental, con una muestra de
34 datos experimentales Vis de la estación de la Universidad Heredia. Se utilizó el Software Matlab 19.0 y se aplicó
el método paramétrico: minimización de un modelo matemático mediante el algoritmo básico con corrección de
Armijo mediante el backtracking. El resultado obtenido es un optimizador factible de ξ en 33 iteraciones, la cual
determinó un modelo de transferencia de la zona Vis de turbidez por aerosoles () 3,69x10-2 y capa de ozono (l)
57,40x10-2 cm con parámetros estadísticos de incertidumbre 0,132%, 2,066% para el error de sesgo medio relativo
(rMBE) y error cuadrático medio relativo (rRMSE) respectivamente. Se concluyó que la atmósfera de la Universidad
Heredia presenta un cielo blanco alternado turbio sin agujero de ozono.
Palabras clave: Turbidez por aerosoles, función espectral, corrección de Armijo, zona Vis, minimización irrestricta
Keywords: Aerosol turbidity, spectral function, Armijo correction, Vis zone, unrestricted minimization
*[email protected]
Recibido: 22/11/2023
Aceptado:24/07/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.09
CC 4.0
2021), el problema de ubicar instalaciones desagradables 0,610 µm atenuada que llega directamente a la superficie de la
(Alamatsaz et al., 2021), etc. Tierra, considera rayos monocromáticos extraterrestres (io) en
(Wm-2µm-1) incidentes sobre la atmósfera alta, que varían
1.3 Relación de la minimización irrestricta y la Zona Vis: debido a que esta capa está compuesta de aerosoles, ozono y
Aplicaciones y Relevancia vapor de agua distribuidos irregularmente en todo su volumen,
produciendo efectos simultáneos combinados de
Los métodos de minimización irrestricta permiten hacer transmitancias independientes (vea la Ecuación 1) de
visible información oculta en los datos de condiciones dispersión Rayleigh (tr), absorción por ozono (to), aerosoles
óptimas, lo que posibilita así una mejor comprensión y (ta) y vapor de agua (tw) (Lisenko, 2018).
manipulación de los sistemas fotovoltaicos (PV) (Duan et al.,
2023). La zona Vis representa el rango de parámetros y i(,l,,w)=(r/ro)2io()trtotatw (1)
condiciones en el que se puede operar eficazmente un PV, Donde:
mientras que la minimización irrestricta permite explorar y ro=1.496x108 km, es la distancia media Sol-Tierra, r es el radio
optimizar este espacio de manera eficiente. de órbita elíptica de la Tierra alrededor del Sol, en uno de sus
focos (Stutzmann & Csoklich, 2022) y d es el número del día,
Se propone la aplicación de un algoritmo novedoso y eficiente, (r/ro)2=[1+0.033cos(2d/365)]2. (2)
basado en el flujo turbulento de optimización a base de agua La masa de aire relativa mr, en función del ángulo cenital solar
para resolver problemas de flujo de potencia óptimo en redes z en grados (Bai & Zong, 2021) es:
de energía con energía solar fotovoltaica y unidades de turbina
eólica. La minimización irrestricta juega un papel importante mr=[cosz+0.15(93,885-z)-1.253]-1 (Bosca Berga, 1995). (3)
en la determinación de las configuraciones óptimas de la red La masa de ozono mo, en función de la altura ho=22 km donde
de energía, mientras que la Zona Vis se refiere a las se da la máxima concentración de ozono es:
restricciones y consideraciones sobre la generación de energía
mo=[(1+ho6370-1)/(cos2z+2ho6370-1)0,5]. (4)
solar y eólica en ubicaciones específicas (Alghamdi, 2023).
La masa de aire ma, en función de la presión estándar Po=1,013
mb y la presión local P en mb,
La extracción de parámetros de modelos fotovoltaicos
presenta algoritmos de optimización avanzados. La ma=mrP/Po (Lefèvre et al., 2013). (5)
minimización irrestricta se aplica en la búsqueda de los Además,
parámetros óptimos del modelo, mientras que la Zona Vis
influye en la generación de datos de entrada del modelo, como tr(,ma)=exp(-ma/[4(115,641-1,335-2)]), (6)
la irradiancia solar y la velocidad del viento (Eslami et al., es llamada dispersión Rayleigh (Gueymard & Kambezidis,
2022). 2004),
to(ko,l,mo)=exp(-kolmo) (Zo et al., 2014), (7)
El desarrollo de un modelo de programación para optimizar la
siendo ko el coeficiente de absorción de ozono en cm-1, l el
energía de una concentración de energía solar de un sistema espesor de ozono en cm,
fotovoltaico (PV). La minimización irrestricta se utiliza para
maximizar las ganancias y determinar un calendario de energía ta(,,ma)=exp(--,ma) (Costa et al., 2012), (8)
rentable, mientras que la Zona Vis se refiere a las fluctuaciones siendo el coeficiente de turbidez por aerosoles sin unidades,
en la disponibilidad de energía solar y eólica a lo largo del =102,74x10-2 para <0,5 µm y =1,206 para >0,5 µm, ka
tiempo (Hamilton et al., 2020). el coeficiente de absorción por aerosoles dada en cm-1,
0
ξ(a)=Σ[ℋ(a)-1]2 (10)
1 2 3 4
a3.
Longitud de onda (µm) Donde ℋ(a)=i(,a)/ie sin unidades, siendo a=[a1 a2 a3]t=[l
Figura 2. Banda espectral solar extraterrestre io e incidente directa w]t.
atenuada ie sobre la superficie terrestre (Wright, 2003) La Ecuación 10 es la función a minimizar y corresponde al
problema de minimización sin restricciones.
Tabla 1. Datos experimentales de ie, io (Wright, 2003), ko y kw A partir de la Tabla 1, se puede observar que kw=0 para la
(Iqbal, 1983) zona Vis, se sustituye ese valor. En la Ecuación 9 se obtiene:
𝜆 io() i e ko
kw
(µm) (Wm-2µm-1) (Wm-2µm-1) (cm-1) tw(w,(0),mr)=exp[-0,2385w(0)mr/(1+20,07w(0)mr)0,45], (11)
0,450 208,852x101 103,261x101
0,003 0 entonces,
0,450 203,469x101 114,130x101 0,003 0
0,455 201,316x101 115,217x101 0,004 0
tw(w,(0),mr)=1, (12)
0,460 200,239x101 117,391x101 0,006 0 reemplazando la Ecuación 12 en la Ecuación 1 obtenemos,
0,465 199,163x101 118,478x101 0,008 0
0,470 198,086x101 119,565x101 0,009 0 i(,l,)=(r/ro)2io()tr(,ma)to(ko,l,mo)ta(,,ma),
(13)
0,475 194,856x101 119,565x101 0,012 0 esta ecuación reduce un grado de libertad en la Ecuación 10,
0,480 192,703x101 120,652x101 0,014 0 se obtiene:
0,485 190,550x101 122,826x101 0,021 0
0,490 192,703x101 122,826x101 0,021 0 ξ(a)=Σ[ℋ(a)-1]2 (14)
0,495 191,627x101 123,913x101 0,025 0 a2.
0,500 192,703x101 125,000x101 0,030 0 Con a=[a1 a2]t=[l ]t, queda la función objetivo simplificada a
0,505 189,474x101 126,087x101 0,035 0 ξ: 2→ con región de factibilidad ϑ={ a2} (Dwail &
0,510 187,321x101 127,174x101 0,040 0
Shiker, 2020).
0,515 190,550x101 127,174x101 0,045 0
0,520 186,244x101 128,261x101 0,048 0 Cuando se obtiene un minimizador local a* de la Ecuación 14,
0,525 181,938x101 129,348x101 0,057 0 la diferencia entre los valores de i(,a) y ie se minimiza.
0,530 188,397x101 130,435x101 0,063 0 Una estrategia para resolver el problema de optimización sin
0,535 190,550x101 130,435x101 0,070 0 restricciones es:
0,540 189,474x101 130,435x101 0,075 0 Paso 1: Resuelva el sistema de ecuaciones
0,545 188,397x101 129,348x101 0,080 0
0,550 187,321x101 129,348x101 0,085 0 grad ξ(a)=0, (15)
0,555 188,397x101 129,348x101 0,095 0 y obtenga una solución a*.
0,560 185,167x101 128,261x101 0,103 0 Paso 2: Verificar si la matriz hessiana es a* es definida
0,565 183,014x101 128,261x101 0,110 0
positiva, es decir hessian ξ(a)>0. Si fuera el caso, entonces a*
0,570 184,091x101 128,261x101 0,120 0
0,575 181,938x101 127,174x101 0,122 0
es un minimizador local estricto del problema de optimización.
0,580 179,785x101 127,174x101 0,120 0 Esta estrategia puede ser implementada en un programa en
0,585 178,708x101 128,261x101 0,118 0 Matlab 19.0 mediante el Algoritmo básico con corrección de
0,590 177,632x101 127,174x101 0,115 0 Armijo usando backtracking, dado de la forma:
0,595 175,478x101 126,087x101 0,120 0 Sean 0,1, a(k) y d(k), tal que
0,600 174,402x101 125,000x101 0,125 0
visualizar que la mínima altura correspondería a l=57,40x10-2 También, mostramos los resultados de los parámetros de
cm y =3,69x10-2. incertidumbre del modelo de transferencia de la zona Vis en la
Tabla 4. Igualmente, se observa que las incertidumbres son
pequeñas, siendo un buen indicativo que el modelo de
transferencia de radiación da una buena estimación de la zona
1.2
Vis.
Función objetivo simplificado
0.6
1250
0.4
1200
0.2 1150
0
0 1100
0.04 r de
eso (cm)
0.03 0.02 0.5
Coeficiente de turbidez por 0.01 0 sp
E ono
aerosoles 1050
oz
Figura 4. Ilustración de ξ en relación a l y utilizando Matlab 19.0 1000
440 460 480 500 520 540 560 580 600 620
Longitud de onda (nm)
Para corroborar esta afirmación mediante la Figura 5, ξ en
función de mantiene constante el dato optimizado de l. a)
También, se visualiza que la función ξ tiene su menor altura en
Intensidad de radiación Vis (W/m2µm)
10
1600
9
1400
Función objetivo simplificado
7 1200
6
1000
440 460 480 500 520 540 560 580 600 620
5 Longitud de onda (nm)
4
b)
3
Figura 6. Visualización comparativa a) entre ie y i, b) entre io y i
2 utilizando Matlab 19.0
función objetivo, lo que limita a encontrar minimizadores Alamatsaz, K., Fatemi Ghomi, S. M. T., & Iranpoor, M.
locales. (2021). Minimal covering unrestricted location of
obnoxious facilities: bi-objective formulation and a case
Los resultados obtenidos a través del enfoque de minimización study. OPSEARCH, 58(2), 351–373.
irrestricta en la optimización de parámetros l, del modelo de https://doi.org/10.1007/s12597-020-00487-0
la función objetivo simplificada en la zona Vis han demostrado Alghamdi, A. S. (2023). Optimal Power Flow of Hybrid
ser válidos y robustos en varios aspectos. La validez de los Wind/Solar/Thermal Energy Integrated Power Systems
resultados se sustenta en la mejora de la precisión y robustez Considering Costs and Emissions via a Novel and
en la estimación de algunos parámetros del modelo Efficient Search Optimization Algorithm. Applied
Sciences (Switzerland), 13(8).
fotovoltaico de Dan (Duan et al., 2023). Esta mejora en la
https://doi.org/10.3390/app13084760
precisión contribuye a una extracción más confiable de los
Aliyeva, Y. N., Mammadova, K. A., & Huseynova, A. N.
parámetros tanto la zona Vis como de módulos fotovoltaicos,
(2022). A Theoretical Model of RGB Attenuation of
lo que sugiere que el enfoque de minimización irrestricta es Solar Radiation Components Under Strong Aerosol
efectivo en mejorar la calidad de los resultados obtenidos. Pollution of the Atmosphere. Aerosol Science and
Debido a que nuestra función objetivo simplificada tuvo baja Engineering, 1-4. https://doi.org/10.1007/s41810-021-
dimensionalidad, su convergencia fue rápida en comparación 00120-7
con la de los modelos de Premkumar y Li (Li et al., 2023; Bai, J., & Zong, X. (2021). Global solar radiation transfer and
Premkumar et al., 2023). its loss in the atmosphere. Applied Sciences
(Switzerland), 11(6).
Una posible fuente de error fue la sensibilidad a la escala de https://doi.org/10.3390/app11062651
las variables l, y los criterios del algoritmo de Armijo Basith, A., Ulin Nuha, M., Prastyani, R., & Winarso, G.
respecto a la terminación de la iteración, que afectó la (2019). COMMUNICATIONS IN SCIENCE AND
convergencia del algoritmo y por ende la precisión de los TECHNOLOGY Aerosol optical depth (AOD) retrieval
resultados. El algoritmo de minimización irrestricta pudo ser for atmospheric correction in Landsat-8 imagery using
sensible a las variaciones de las condiciones iniciales second simulation of a satellite signal in the solar
provocando soluciones muy diferentes, lo que podría afectar la spectrum-vector (6SV). In Communications in Science
precisión de la estimación de los parámetros del modelo de la and Technology (Vol. 4, Issue 2).
zona Vis. https://earthexplorer.usgs.gov/
Bosca Berga, J. V. (1995). Contribucion al estudio de la
Se sugiere establecer la técnica de algoritmo de Armijo para el radiacion solar y de la determinacion de la turbiedad
diseño de procesamiento de optimización de datos de la zona atmosferica. Aplicacion a Valencia y Sevilla.
Vis, esto puede implicar el uso de sensores, controladores, [Universitat Politècnica de València]. In Riunet.
software y hardware especializado con el fin de lograr https://doi.org/10.4995/Thesis/10251/4502
resultados más favorables. Borah, P., & Gupta, D. (2020). Unconstrained convex
minimization based implicit Lagrangian twin extreme
CONCLUSIONES learning machine for classification (ULTELMC).
Applied Intelligence, 50(4), 1327–1344.
Conforme con los hallazgos obtenidos, se deduce que se logró https://doi.org/10.1007/s10489-019-01596-0
estimar la Zona Visible mediante el cálculo de modelos: Costa, R. S., Martins, F. R., & Pereira, E. B. (2012). Aerossóis
transferencia de la zona Vis, modelo matemático “función atmosféricos e a quantificação do recurso energético
objetivo simplificada” y uso de la técnica “algoritmo básico solar: experimentos em modelo de transferência
con corrección de Armijo usando backtracking”. Mediante la radiativa. Revista Brasileira de Energia Solar, 3(1), 62-
implementación del algoritmo de Armijo en un programa en el 62. https://doi.org/10.59627/rbens.2012v3i1.79
Software Matlab 19.0, se obtuvo un minimizador local a* de Cox, J. L., Hamilton, W. T., Newman, A. M., Wagner, M. J.,
una función objetivo sin restricciones. A partir de los datos de & Zolan, A. J. (2023). Real-time dispatch optimization
minimizador, concluimos que los datos empíricos de radiación for concentrating solar power with thermal energy
de la zona Vis pertenecen a una atmósfera sin agujero de ozono storage. Optimization and Engineering, 24(2), 847–884.
con cielo blanco alternado turbio para la estación del https://doi.org/10.1007/s11081-022-09711-w
Departamento Académico de la Universidad de Heredia. de Wijn, A. G., Casini, R., Carlile, A., Lecinski, A. R., Sewell,
La ilustración de ξ en relación con permite determinar la S., Zmarzly, P., Eigenbrot, A. D., Beck, C., Wöger, F.,
restricción del valor del coeficiente de turbidez por aerosoles, & Knölker, M. (2022). The Visible Spectro-Polarimeter
el cual podría ser utilizado para el uso del método de of the Daniel K. Inouye Solar Telescope. Solar Physics,
optimización con restricciones de desigualdad. 297(2). https://doi.org/10.1007/s11207-022-01954-1
Di Pietro, S. (2022). Processes of urban transition to
REFERENCIAS autonomous decentralized systems of renewable energy.
Estudios Demograficos y Urbanos, 37(3), 807–837.
Abunima, H., Teh, J., & Jabir, H. J. (2019). A new solar https://doi.org/10.24201/edu.v37i3.2073
radiation model for a power system reliability study. Drozdowski, M., & Shakhlevich, N. V. (2021). Scheduling
IEEE Access, 7, 64758–64766. divisible loads with time and cost constraints. Journal of
https://doi.org/10.1109/ACCESS.2019.2916168
Scheduling, 24(5), 507–521. Based on Tri Thin Layer Material in Visible Spectrum
https://doi.org/10.1007/s10951-019-00626-6 for Solar Cell Energy Harvesting. Plasmonics, 15(4),
Duan, Z., Yu, H., Zhang, Q., & Tian, L. (2023). Parameter 1061–1069. https://doi.org/10.1007/s11468-020-01132-
Extraction of Solar Photovoltaic Model Based on 8
Nutcracker Optimization Algorithm. Applied Sciences Ilango, R., Loff, B., & Oliveira, I. C. (2020). NP-Hardness of
(Switzerland), 13(11). circuit minimization for multi-output functions. Leibniz
https://doi.org/10.3390/app13116710 International Proceedings in Informatics, LIPIcs, 169.
Dwail, H. H., & Shiker, M. A. K. (2020). Using a trust region https://doi.org/10.4230/LIPIcs.CCC.2020.22
method with nonmonotone technique to solve Iqbal, M. (1983). A CLOUDLESS-SKY ATMOSPHERE
unrestricted optimization problem. Journal of Physics: AND ITS OPTICS. In An Introduction to Solar
Conference Series, 1664(1). Radiation (pp. 85–105). Elsevier.
https://doi.org/10.1088/1742-6596/1664/1/012128 https://doi.org/10.1016/b978-0-12-373750-2.50010-0
Emiola, I., Adem, R. (Junio, 2021). Comparison of Jiang, H., & Dong, Y. (2016). A nonlinear support vector
Minimization Methods for Rosenbrock Functions. Work machine model with hard penalty function based on
presented at 2021 29th Mediterranean Conference on glowworm swarm optimization for forecasting daily
Control and Automation (MED), PUGLIA, pp. 837-842. global solar radiation. Energy Conversion and
https://doi.org/10.1109/MED51440.2021.9480200 Management, 126, 991–1002.
Eslami, M., Akbari, E., Seyed Sadr, S. T., & Ibrahim, B. F. https://doi.org/10.1016/j.enconman.2016.08.069
(2022). A novel hybrid algorithm based on rat swarm Kairouz, P., Ribero, M., Rush, K., & Thakurta, A. (2020). Fast
optimization and pattern search for parameter extraction dimension independent private adagrad on publicly
of solar photovoltaic models. Energy Science and estimated subspaces. arXiv preprint arXiv:2008.06570
Engineering, 10(8), 2689–2713. Lefèvre, M., Oumbe, A., Blanc, P., Espinar, B., Gschwind, B.,
https://doi.org/10.1002/ese3.1160 Qu, Z., Wald, L., Schroedter-Homscheidt, M., Hoyer-
Farag, M. M., Patel, N., Hamid, A. K., Adam, A. A., Bansal, Klick, C., Arola, A., Benedetti, A., Kaiser, J. W., &
R. C., Bettayeb, M., & Mehiri, A. (2023). An Optimized Morcrette, J. J. (2013). McClear: A new model
Fractional Nonlinear Synergic Controller for Maximum estimating downwelling solar radiation at ground level
Power Point Tracking of Photovoltaic Array under in clear-sky conditions. Atmospheric Measurement
Abrupt Irradiance Change. IEEE Journal of Techniques, 6(9), 2403–2418.
Photovoltaics, 13(2), 305–314. https://doi.org/10.5194/amt-6-2403-2013
https://doi.org/10.1109/JPHOTOV.2023.3236808 Li, J., Qin, C., Yang, C., Ai, B., & Zhou, Y. (2023). Extraction
Fei, L., Yin, Y., Yang, M., Zhang, S., & Wang, C. (2021). of Single Diode Model Parameters of Solar Cells and PV
Wearable solar energy management based on visible Modules by Combining an Intelligent Optimization
solar thermal energy storage for full solar spectrum Algorithm with Simplified Explicit Equation Based on
utilization. Energy Storage Materials, 42, 636–644. Lambert W Function. Energies, 16(14).
https://doi.org/10.1016/j.ensm.2021.07.049 https://doi.org/10.3390/en16145425
Ferry, M. W., Gill, P. E., Wong, E., & Zhang, M. (2021). Li, Y., Liu, Z., & Liu, H. (2019). A subspace minimization
Projected-Search Methods for Bound-Constrained conjugate gradient method based on conic model for
Optimization. http://arxiv.org/abs/2110.08359 unconstrained optimization. Computational and Applied
Filos-Ratsikas, A., Kanellopoulos, P., Voudouris, A. A., & Mathematics, 38(1). https://doi.org/10.1007/s40314-
Zhang, R. (2023). Settling the Distortion of Distributed 019-0779-7
Facility Location. http://arxiv.org/abs/2301.01604 Lisenko, S. A. (2018). Atmospheric correction of multispectral
Franc, V., Yermakov, A., Balasubramanian, V. N., & Tsang, satellite images based on the solar radiation transfer
I. (2021). Learning Maximum Margin Markov approximation model. Atmospheric and Oceanic Optics,
Networks from examples with missing labels. In 31(1), 72–85.
Proceedings of Machine Learning Research (Vol. 157) https://doi.org/10.1134/S1024856018010116
Ghosh, S., Acharya, T., & Maity, S. P. (2019). On outage Manoharan, H., Teekaraman, Y., Kuppusamy, R., &
minimization in RF energy harvesting relay assisted Radhakrishnan, A. (2021). A Novel Optimal Robotized
bidirectional communication. Wireless Networks, 25(7), Parking System Using Advanced Wireless Sensor
3867–3881. https://doi.org/10.1007/s11276-018-01924- Network. Journal of Sensors, 2021.
1 https://doi.org/10.1155/2021/2889504
Gueymard, C. A., & Kambezidis, H. D. (2004). SOLAR Masoom, A., Kosmopoulos, P., Bansal, A., & Kazadzis, S.
RADIATION AND DAYLIGHT MODELS. In Solar (2020). Solar energy estimations in india using remote
Spectral Radiation (pp. 221-301). Elsevier. sensing technologies and validation with sun
https://doi.org/10.1016/B978-075065974-1/50013-9 photometers in urban areas. Remote Sensing, 12(2).
Hamilton, W. T., Husted, M. A., Newman, A. M., Braun, R. https://doi.org/10.3390/rs12020254
J., & Wagner, M. J. (2020). Dispatch optimization of Mims, F. M. (2022). A 30-Year Climatology (1990-2020) of
concentrating solar power with utility-scale Aerosol Optical Depth and Total Column Water Vapor
photovoltaics. Optimization and Engineering, 21(1), and Ozone over Texas. Bulletin of the American
335–369. https://doi.org/10.1007/s11081-019-09449-y Meteorological Society, 103(1), E101–E109.
Hoque, A., & Islam, M. T. (2020). Numerical Analysis of https://doi.org/10.1175/BAMS-D-21-0010.1
Single Negative Broadband Metamaterial Absorber
1
Universidad César Vallejo, Escuela Profesional de Ingeniería Civil, Lima, Perú
2
Universidad César Vallejo, Escuela Profesional de Ingeniería Civil, Lima, Perú
Resumen: El problema del concreto hidráulico en canales abiertos se deriva de factores de diseño, construcción
deficiente y condiciones climatológicas adversas que hacen que el concreto sea menos duradero y propenso a
desarrollar defectos superficiales que debilitan su resistencia física y mecánica. El estudio analizó el efecto de la
ceniza de cartón (CC) y cerámica reciclada (RC) en las propiedades físicas y mecánicas del concreto hidráulico. La
investigación fue de diseño experimental, que incluyó ensayos de absorción de agua, permeabilidad, resistencia a la
compresión y resistencia a la flexión. Se utilizaron diferentes combinaciones para sustituir el cemento en proporciones
de CC 1,5 % + RC 1 %, CC 3 % + RC 2 %, CC 4,5 % + RC 3 %, y CC 6 % + RC 4 %, que luego se compararon con
la norma ASTM y una muestra control. Los resultados indicaron una absorción de agua de 0,061 g/s y una
permeabilidad de 297,09 g/m2, ambas inferiores. Las mayores resistencias mecánicas a los 7, 14 y 28 días alcanzaron
resistencias a compresión de 191,75 ± 1,22 kg/cm2, 274,38 ± 3,06 kg/cm2, 323,34 ± 6,85 kg/cm2, y resistencias a
flexión de 31,32 ± 0,04 kg/cm2, 33,23 ± 0,16 kg/cm2 y 35,76 ± 0,31 kg/cm2. Se concluye que la combinación óptima
para los ensayos físicos y mecánicos se consigue con la combinación de CC 3 % y RC 2 %; para una mayor cantidad
de CC y RC respecto a la muestra óptima, el concreto tiende a deteriorarse y perder resistencia.
Palabras clave: Concreto hidráulico, ceniza de cartón, cerámico reciclado, resistencia física y mecánica
Keywords: Hydraulic concrete, cardboard ash, recycled ceramic, physical and mechanical resistance
*[email protected]
Recibido:23/11/2023
Aceptado:19/07/2024
Publicado en línea: 31/08/2024
10.33333/rp.vol54n1.10
CC 4.0
Los aditivos que se utilizan en los materiales cementosos resistencia a las acciones externas debido a la dureza del
modernos para mejorar sus características mecánicas y de material.
durabilidad, su uso generalizado en la fabricación del concreto
los hace notoriamente muy costosos (Paul et al., 2022). El Por su parte, Fernández-Torrez et al. (2022) compararon
concreto ha encontrado un amplio uso en aplicaciones de concretos fabricados con residuos de caucho de neumáticos
construcción, por lo que analizar el reciclaje y la reutilización (RCN); utilizaron porcentajes de 5 %, 10 % y 20 % a las
del polvo de cerámico procedente de los escombros de edades de 7 y 14 días. Los resultados muestran que los RCN
construcción como material cementante suplementario en la superiores al 5 % reducen considerablemente la mala
fabricación del concreto es una muy buena opción (Zhang et absorción en un 70 %, pero pierde resistencia a compresión y
al., 2023). En comparación con el uso de agregados naturales, flexión en relación con la muestra sin contenido de RCN.
la utilización de residuos de construcción y demolición en la
producción de concreto y mortero tiene un impacto Salih et al. (2022), en su investigación, sustituyeron
medioambiental reducido (Barrios et al., 2021). parcialmente el Cemento Portland Ordinario hasta en un 71 %,
utilizando cenizas volantes, escoria granulada de alto horno
Las actividades industriales, construcción y demolición (GGBS) y microsílice para producir concreto autocompactante
producen grandes cantidades de residuos, como los residuos duraderos y de alta resistencia. Se diseñaron dos grupos de
de ladrillos, cerámicas, el papel y el cartón (de Matos et al., mezclas. El primero contenía 14 %, 23,4 % y 32,77 % de
2021). Los restos de papel y cartón contribuyen en gran cenizas volantes y 6,4 % de microsílice. El segundo contenía
medida a la contaminación de los vertederos y los centros de 32,77 %, 46,81 % y 65,5 % de GGBS y 6,4 % de microsílice.
producción. El sector manufacturero produce cada año más de La mayoría de las pruebas reveló una alta resistencia a la
estos residuos que contribuyen a la contaminación y la compresión y durabilidad, la permeabilidad y la absorción al
demanda de ampliación de vertederos van en aumento, al igual agua presentó un efecto significativo.
que el consumo de energía y recursos naturales que son de
preocupación mundial (Solahuddin & Yahaya, 2021). Saad et al. (2020), en su estudio evaluaron el comportamiento
del uso de nuevos materiales en pavimentos de concreto rígido
A nivel mundial existe un movimiento creciente para reciclar a las edades de 7 y 14 días. Ellos utilizaron dosificaciones de
y reutilizar los residuos de construcción y demolición, se humo de sílice (0 % + 20 %) y caucho (0 % + 25 %). Los
calcula un aproximado del 35 % de la cantidad producida se resultados mostraron que la resistencia a la compresión simple
trasladan a vertederos sin ningún tratamiento posterior disminuía hasta un 20 % cuando se añadía caucho granulado
(Menegaki & Damigos, 2018). En consecuencia, la al 25 %. La adición de humo de sílice aumentó la resistencia a
sostenibilidad de la construcción y el reciclado de los residuos la compresión y a la flexión.
se están convirtiendo en preocupaciones medioambientales de
primer orden (Amin et al., 2023). Según Sadagopan et al. (2023), en su investigación se
utilizaron agregados de concreto triturado (ACT) y GGBS en
Los materiales reciclados procedentes de demoliciones de reemplazo de parcial del cemento en relaciones a/c de 0,42 y
construcciones con bajas emisiones de carbono podrían 0,48. Los agregados resultantes muestran una mayor densidad,
producirse para sustituir al Clinker Portland como polvo de lo que contribuye a un aumento en la resistencia a la
cemento endurecido y pueden ofrecer una forma de reducir el compresión del concreto. Por tanto, el preprocesamiento
CO2 y la carbonatación mineral en los materiales cementosos. mecánico se puede mejorar para producir ACT, lo que
Su aplicación como tecnología de producción de cemento favorece a un comportamiento del concreto más dúctil.
contribuiría a contrarrestar las consecuencias del aumento de
las emisiones de gases de efecto invernadero (Li et al., 2022). El estudio unió los agregados para que actúen como adhesivo
Por ello, varios investigadores han empezado a explorar y den la resistencia deseada al concreto, al sustituir el cemento
materiales de construcción cementosos alternativos que por combinaciones de CC y RC en proporciones específicas
puedan utilizarse como sustitutos parciales del cemento mediante el análisis de los efectos al utilizan materiales
Portland (Demissew, 2022). reciclados. De este modo, se reducen los daños causados por
el cemento al medio ambiente y, al mismo tiempo, los costos
A medida que aumenta la necesidad del concreto de alta excesivos generados por los proyectos de construcción. Las
resistencia en diversos contextos, ganan popularidad combinaciones utilizadas en el estudio se prepararon en
productos como los compuestos de concreto hidráulicos proporciones variables de CC y RC, concretamente CC 1,5 %
fabricados a partir de residuos industriales o reciclados. + RC 1 %, CC 3 % + RC 2 %, CC 4,5 % + RC 3 %, y CC 6 %
Además de su resistencia, los compuestos tienen otras + RC 4 %. Asimismo, se compararon con la norma ASTM y
propiedades deseables, como baja contracción por secado, baja con una muestra control sin contenido de CC y RC.
permeabilidad y resistencia a productos químicos y altas
temperaturas (Alhazmi et al., 2021). El objetivo principal de este estudio es analizar el efecto de la
ceniza de cartón (CC) y cerámica reciclada (RC) en las
Brahimi et al. (2022) estudiaron el uso de residuos con material propiedades físicas y mecánicas del concreto hidráulico. Los
de Polietileno Tereftalato para sustituir al cemento, utilizaron materiales proceden de vertederos y demoliciones de obras de
porcentajes de 3,0 g, 1,75 g, 1,0 g y 2,25 g. Los resultados construcción a escala mundial, materiales que contribuyen a la
revelaron que los geo polímeros preparados contienen un gran descontaminación medioambiental y mejoras de las
porcentaje de cuarzo en forma de sílice, presentando una mejor propiedades físicas y mecánicas del material compuesto.
La importancia del estudio es especialmente relevante, porque Fase tres: consistió en la recogida de agregados pétreos. Se
con frecuencia se desechan cantidades sustanciales de estos eligió la cantera de Vicho ubicada en la región Cusco como
materiales, que acaban en vertederos o escombreras. Esta fuente de agregados gruesos (AG), mientras que la cantera de
investigación abordó la importancia de utilizar materiales Huambutio fue seleccionada para los agregados finos (AF).
reciclables en el concreto hidráulico mediante el diseño de Estas canteras se consideraron óptimas para la producción de
mezclas y la realización de ensayos de resistencia física y concreto por ofrecer materias primas de alta calidad, cumplir
mecánica que aportaron en la mitigación de los efectos las normas técnicas aplicables y ser fácilmente accesibles por
perjudiciales del deterioro en el concreto hidráulico. carretera. Este factor de accesibilidad contribuye a reducir el
costo por m3 de material entregado en la obra.
2. METODOLOGÍA
Fase cuatro: consistió en la producción de probetas cilíndricas
La investigación se clasificó como aplicada, porque resolvió y rectangulares, que posteriormente, se sometieron a curado de
los problemas asociados al concreto hidráulico. Esta disciplina intervalos de tiempo específicos de 7, 14 y 28 días naturales
académica se caracterizó por su énfasis en la aplicación para la resistencia mecánica. Asimismo, las pruebas de
práctica de la información adquirida. Además, el estudio absorción de agua y de permeabilidad se realizaron con
utilizó un diseño experimental, porque varió la variable probetas en edad de 28 días, debido a que este es el tiempo
independiente mediante una manipulación deliberada para ver necesario para determinar su máxima capilaridad y resistencia.
su impacto en la variable dependiente.
2.2. Estudios de laboratorio:
La población la conformaron 240 probetas cilíndricas de 300
mm x 15 mm de diámetro y 60 probetas rectangulares de 42 En esta sección, se realizaron pruebas de resistencia a la
cm x 15 cm x 12 cm, las mismas que se dividieron en una compresión, a la flexión, absorción de agua, y permeabilidad
muestra de 180 probetas para ensayos de resistencia mecánica mediante ensayos de laboratorio estandarizados por las normas
(compresión y flexión) y una muestra de 60 probetas para internacionales y nacionales como son:
ensayos de resistencia física (absorción y permeabilidad).
• ASTM C33/CMMM-18 (2018)
La investigación utilizó la técnica de observación directa, • ASTM D854-14 (2023)
porque los investigadores participaron activamente en todas • ASTM C1585-20 (2020)
las fases de elaboración de las muestras y en las pruebas de • ACI PRC-522 R-10 (2011)
laboratorio, incluidas las evaluaciones físicas y mecánicas. El • IRAM1871 (2021)
método de muestreo empleado fue no probabilístico, es decir, • ASTM C33/C39M-21 (2021)
basado en el juicio, tras considerar detenidamente las
• ASTM C78/C78M (2022)
características tanto de la población como de la muestra,
• NTP 339.034 (2021)
teniendo en cuenta las directrices de la Norma Técnica Peruana
• NTP 339.078 (2022)
(NTP) y la norma ASTM (American Society for Testing and
Materials). • MTC E-107-200 (2016)
2.1. Estudios de Campo: Se realizó el ensayo tanto de los agregados gruesos como de
los finos, y los valores resultantes se presentan en la Tabla 1.
En el presente estudio, los investigadores supervisaron
Tabla 1. Tipologías de los materiales naturales y reciclados
minuciosamente cada uno de los procedimientos de la Material AG AF CC RC
recolección de material reciclado, agregados pétreos naturales Peso específico (g/cm3) 2,49 2,51 2,22 2,62
y los ensayos de laboratorio, que se llevaron a cabo en cuatro Absorción (%) 1,16 0,76 - -
fases, las cuales se detallan a continuación. Contenido de humedad (%) 2,03 1,37 1,01 1,45
Módulo de fineza 6,90 2,37 0,16 0,26
Peso compactado seco (kg/m3) 1 596,45 1 630,02 1 327,21 1 439,33
Fase uno: consistió en la puesta en marcha del procedimiento Peso suelto seco (kg/m3) 1 580,87 1 495,56 1 203,27 1 356,92
de recogida de cartón por parte de la organización “Pro- TMN (pulg) 1/2" N° 4 - -
Ecología”. Esta empresa se encarga de la selección y Nota: Agregado Grueso (AG), Agregado Fino (AF), Ceniza de Cartón (CC),
Cerámico Reciclado (RC), Tamaño Máximo Nominal (TMN)
recuperación sistemática de materiales de papel y cartón en
diversos puntos de la ciudad del Cusco. El material de cartón Se formuló, de acuerdo con la metodología ACI-211.1. (2002),
se sometió a combustión para crear cenizas, que luego se un diseño de mezcla plástica de 280 kg/cm2. En el diseño de la
refinaron hasta convertirlas en polvo. El peso específico del mezcla se utilizó cemento portland IP, conocido por su
polvo resultante fue evaluado mediante ensayos de laboratorio, gravedad específica de 2,85 g/cm3. Los resultados obtenidos
lo cual arrojó un valor de 2,22 g/cm3. del diseño de la mezcla plástica se presentan en la Tabla 2.
Fase dos: consistió en la recuperación de los residuos de Tabla 2. Diseño de mezcla plástica
cerámica generada durante la construcción del Pabellón "B" en Dosificación 280 kg/cm2
Material
el campus Cusco de la Universidad Continental. Volumen Tandas Valores de diseño
Posteriormente, la sustancia fue sometida a pulverización y Cemento 1,0 pie3/pie3 42,5 kg/bolsa 528 kg/m3
Agua de diseño 16,0 l/pie3 16,0 l/bolsa 205 kg/m3
posteriormente a medición de densidad, alcanzando una Agregado fino 0,9 pie3/pie3 39,4 kg/bolsa 490 kg/m3
gravedad específica de 2,62 g/cm3. Agregado grueso 1,8 pie3/pie3 80,0 kg/bolsa 994 kg/m3
A partir de los datos obtenidos en el diseño de la mezcla, fue Tabla 3 muestra, a continuación, los valores obtenidos
posible determinar las cantidades de materiales empleadas en relacionados con las cantidades necesarias para los cuatro
la preparación de las probetas cilíndricas y rectangulares. La diseños y una muestra de control sin contenido de CC y RC.
3. RESULTADOS Y DISCUSIÓN
La Tabla 4 presenta los resultados de la variación de la El estudio realizado por Fernández-Torrez et al. (2022),
velocidad de absorción de agua en el concreto hidráulico con demostró que la utilización de caucho de neumáticos dio lugar
ceniza de cartón (CC) y cerámico reciclado (RC), denotando a una reducción de la absorción de agua del 2,72 %, 0,67 % y
para la muestra control un valor de 0,101 g/s ± 0,047 g/s, la 0,55 %. Los resultados indican que el uso de materiales
combinación CC 1,5 % + RC 1,0 % alcanzó una absorción de reciclados mejora la absorción de agua en diferentes
0,072 g/s ± 0,048 g/s, la combinación CC 3,0 % + RC 2,0 % porcentajes. El presente estudio logró con éxito una reducción
obtuvo un absorción de 0,061 g/s ± 0,045 g/s, la combinación de la porosidad en la absorción de agua de 0,061 g/s utilizando
CC 4,5 % + RC 3,0 % logró una absorción de 0,085 g/s ± 0,043 una combinación de CC 3,0 % y RC 2,0 %. Está mezcla en
g/s, y la combinación CC 6,0 % + 4,0 % obtuvo un absorción particular se considera la más favorable en términos de
de 0,090 g/s ± 0,045 g/s, respectivamente. minimización de la velocidad de absorción de agua. Para
incrementos mayores de CC y RC, la velocidad de absorción
mostró una tendencia a aumentar.
3.2. Permeabilidad:
La Tabla 5 muestra los resultados de la permeabilidad en la curva (S1/2) que representan las masas en g/m2, que es el
mediante el ensayo de succión capilar. Los valores variables valor del incremento de masa por unidad de área en dos
momentos diferentes: 60 y 294 horas. La muestra control observaron valores de permeabilidad del 33,3 %, 40 %, 47 %
alcanzó un valor promedio de 867,09 g/m2, la combinación CC y 53,3 %, lo que demostró una influencia sustancial en la
1,5 % + RC 1,0 % obtuvo un valor de 366,13 g/m2, la mejora de la capacidad de absorción de agua. En este estudio,
combinación CC 3,0 % + RC 2,0 % alcanzó un valor de 297,09 la utilización de materiales reciclados demuestra una notable
g/m2, la combinación CC 4,5 % + RC 3,0 % presentó un valor reducción de la permeabilidad cuando se usan materiales
de 634,98 g/m2, y la combinación CC 6,0 % + RC 4,0 % logró reciclados hasta una combinación de CC 3,0 % y RC 2,0 %.
un valor de 669,69 g/cm2, respectivamente. Asimismo, cuando se usan por encima de la mezcla óptima
estos valores adoptan la forma de una línea creciente que se
Según Salih et al. (2022), se realizó una investigación del vuelve casi horizontal, lo que demuestra que la permeabilidad
impacto de las cenizas volantes y microsílice en la mejora de del concreto aumenta con el tiempo.
la capacidad de forma sofisticada. Los investigadores
La Tabla 6 muestra los datos obtenidos de las variaciones en la compresión para la muestra control de 330 y 445 Mpa y con
la resistencia a la compresión a las edades de 7, 14 y 28 días. las combinaciones obtuvieron 228 y 245 MPa. Según la norma
La muestra control alcanzó valores de resistencia de 191,75 ± ASTM C39 (2021), la resistencia mínima requerida para un
1,22 kg/cm2, 240,67 ± 0,53 kg/cm2 y 282,72 ± 1,91 kg/cm2. La diseño de 280 kg/cm2 a los 7 días es de 191,48 kg/cm2, a los
combinación CC 1,5 % + RC 1,0 % alcanzó resistencias de 14 días 239,22 kg/cm2 y a los 28 días 280 kg/cm2. Sin
207,48 ± 2,87 kg/cm2, 258,60 ± 1,90 kg/cm2 y 305,90 ± 4,35 embargo, la muestra control supera a la norma en un 0,14 % a
kg/cm2, la combinación CC 3,0 % + RC 2,0 % obtuvo los 7 días, 0,61 % a los 14 días y 0,97 % a los 28 días, la
resistencias de 218,88 ± 5,16 kg/cm2, 274,38 ± 3,06 kg/cm2 y combinación de CC 1,5 % + RC 1,0 % supera en un 8,36 % a
323,34 ± 6,85 kg/cm2, la combinación CC 4,5 % + RC 3,0 % los 7 días, 8,10 % a los 14 días y 9,25 % a los 28 días, la
logró resistencias de 190,21 ± 0,67 kg/cm2, 234,19 ± 0,99 combinación de CC 3,0 % + RC 2,0 % superan en un 14,31 %,
kg/cm2 y 276,57 ± 1,17 kg/cm2, la combinación CC 6,0 % + 14,70 % y 15,48 % en los mismos intervalos de tiempo. El
RC 4,0 % alcanzó resistencias de 180,91 ± 1,38 kg/cm2, 223,12 presente estudio alcanzó una mayor resistencia con el uso de
± 2,27 kg/cm2 y 262,87 ± 2,84 kg/cm2, respectivamente. la combinación de CC 3,0 % + 2,0 %, superando a la muestra
control y a la norma ASTM C-39. De esta manera, se conduce
Saad et al. (2020) afirmaron que el uso de un 20 % de humo a una disminución de la resistencia a la compresión cuando se
de sílice y un 25 % de caucho permitió alcanzar resistencias a utilizan cantidades superiores de la combinación óptima.
La Tabla 7 y la Tabla 8 presentan los datos de las variaciones Los resultados del estadístico de la prueba Anova, que se
de la resistencia a la flexión en tres intervalos de tiempo utilizó para contrastar los ensayos físicos (absorción y
diferentes: 7, 14 y 28 días. La muestra control alcanzó una permeabilidad) que se hicieron en el concreto hidráulico,
resistencia de 28,89 ± 0,17 kg/cm2, 30,68 ± 0,12 kg/cm2 y pueden verse en la Tabla 9. Se observó un nivel de
32,93 ± 0,29 kg/cm2. La combinación CC 1,5 % + RC 1,0 % significación de 0,006 para la prueba de absorción, mientras
obtuvo una resistencia de 30,19 ± 0,12 kg/cm2, 32,43 ± 0,16 que se observó un nivel de significación de 0,045 para la
kg/cm2, y 35,02 ± 0,12 kg/cm2. La combinación CC 3,0 % + prueba de permeabilidad. Estas cifras indican que los valores
RC 2,0 % logró resistencias de 31,32 ± 0,04 kg/cm2, 33,23 ± p-valué son inferiores al nivel alfa predeterminado de 0,05.
0,16 kg/cm2 y 35,76 ± 0,31 kg/cm2. La combinación CC 4,5 %
Tabla 10. Anova – Ensayos mecánicos
+ RC 3,0 % alcanzó resistencias de 26,58 ± 0,25 kg/cm2, 27,22
Suma de Media
± 0,42 kg/cm2 y 27,45 ± 0,32 kg/cm2, y la combinación CC 6,0 cuadrados
gl
cuadrática
F Sig.
% + RC 4,0 % lograron resistencias de 24,43 ± 0,35 kg/cm2, Resistencia a la compresión
25,16 ± 0,36 kg/cm2 y 25,56 ± 0,36 kg/cm2, respectivamente. Entre
3395,516 4 848,879 73,104 0,000
grupos
Dentro
En su reciente estudio, Brahimi et al. (2022) examinaron la 7 días
de 116,12 10 11,612
utilización de residuos que contienen material PET como grupos
sustituto del cemento. Los autores observaron que la Total 3511,636 14
incorporación de este material de desecho dio lugar a una Entre
5239,249 4 1309,812 226,06 0,000
mayor resistencia mecánica, atribuida a la dureza del material. grupos
Dentro
Según la norma ASTM C78 (2022), el módulo de rotura de la 14 días
de 57,941 10 5,794
resistencia de diseño 280 kg/cm2 el 10 % es de 28 kg/cm2 como grupos
mínimo estándar para las edades de 7, 14 y 28 días, la muestra Total 5297,19 14
control logro a superar con porcentajes de 3,18 %, 9,57 % y Entre
8177,67 4 2044,417 86,317 0,000
grupos
17,59 %, la combinación de CC 1,5 % + RC 1,0 % superó en
Dentro
porcentajes de 7,84 %, 15,81 % y 25,07 %, la combinación de 28 días
de 236,849 10 23,685
CC 3,0 % + RC 2,0 % superó en porcentajes de 11,85 %, 18,69 grupos
% y 27,70 % en los mismos intervalos de tiempo. La Total 8414,519 14
investigación de los autores y el presente estudio demuestran Resistencia a la flexión
Entre
que el uso de materiales reciclados en cantidades grupos
85,232 4 21,308 301,527 0,000
cuidadosamente controladas tiene el potencial de hacer que el Dentro
7 días
concreto hidráulico sea más resistente a la flexión, por lo que de 0,707 10 0,071
se consigue la mayor resistencia con la mezcla de CC 3,0 % y grupos
RC 2,0 %, a mayor cantidad de esta combinación la resistencia Total 85,939 14
Entre
a la flexión del concreto tiende a disminuir. grupos
133,643 4 33,411 297,921 0,000
Dentro
14 días
3.5. Prueba estadística Anova: de 1,121 10 0,112
grupos
Tabla 9. Anova – Ensayos físicos Total 134,764 14
Suma de Media Entre
gl F Sig. 246,213 4 61,553 483,124 0,000
cuadrados cuadrática grupos
Absorción Dentro
28 días
Entre grupos 0,021 4 0,005 3,844 0,006 de 1,274 10 0,127
Dentro de grupos 0,124 90 0,001 grupos
Total 0,145 94 Total 247,487 14
Permeabilidad
Entre grupos 2 877 079,82 4 719 269,954 2,593 0,045 El análisis estadístico, es decir, la prueba Anova, se llevó a
Dentro de grupos 18 030 505,11 65 277 392,386 cabo para examinar los resultados mostrados en la Tabla 10.
Total 20 907 584,93 69
Los valores significativos obtenidos tanto para la prueba de
resistencia a la compresión como para la prueba de flexión de CC y RC, la resistencia a la flexión tiende a reducirse.
resultaron ser 0,000 para las edades de 7, 14 y 28 días. Los Además, esta conclusión y resultados alcanzados se respaldan
resultados obtenidos son estadísticamente significativos, ya con los datos mostrados en la Tabla 10 del análisis estadístico
que se sitúan por debajo del nivel alfa de significación Anova-Varianza, que alcanzó un valor p de 0,000, inferior al
especificado (p < 0,05). nivel de significación predeterminado de alfa 0,05 para grupos
a las edades de 7, 14 y 28 días. En este contexto, se infiere que
4. CONCLUSIONES la CC y la RC tienen efectos significativos sobre las
propiedades del concreto hidráulico.
El estudio concluye que la velocidad de absorción de agua con
el uso de cenizas de cartón (CC) y cerámica reciclada (RC) en 5. REFERENCIAS
concentraciones de CC 3,0 % y RC 2,0 % disminuye la
porosidad en un promedio de 0,061 g/s ± 0,045 g/s en Abubakar Tadda, M., Ahsan, A., Imteaz, M., Shitu, A.,
comparación con la muestra control que alcanzo 0,101 ± 0,047 Abdulbaki Danhassan, U., & Idris Muhammad, A.
g/s, lo que se debe principalmente por la forma natural de los (2020). Operation and Maintenance of Hydraulic
materiales. Igualmente, estas propiedades se pierden cuando se Structures. In Hydraulic Structures - Theory and
supera la mezcla óptima de CC 3,0 % y RC 2,0%, provocando Applications. IntechOpen.
mayor absorción capilar en el concreto. Esta conclusión y los https://doi.org/10.5772/intechopen.91949
resultados se apoyan en los datos mostrados en la Tabla 9 del ACI PRC-522 R-10 American Concrete Institute. (2011).
análisis estadístico Anova-Varianza, que alcanzó un p-valor de Report on Pervious Concrete (Reapproved 2011).
0,006, inferior al nivel de significación predeterminado de alfa https://www.concrete.org/store/productdetail.aspx?Item
0,05. En este contexto, la CC y la RC tienen efectos positivos ID=52210&Format=PROTECTED_PDF&Language=E
sobre las propiedades físicas del concreto hidráulico, nglish&Units=US_AND_METRIC
recomendándose el uso de estos materiales reciclados. ACI-211.1 American Concreto Institute®. (2002, August 9).
La investigación concluye que la permeabilidad con el uso de Diseño de mezclas Método ACI.
ceniza de cartón (CC) y cerámico reciclado (RC) de una https://www.concrete.org/topicsinconcrete/topicdetail/
combinación de CC 3,0 % y RC 2,0 % alcanzó un valor metodo%20aci%20211?search=metodo%20aci%20211
promedio de 297,09 g/m2, presentando así una capacidad de Alhazmi, H., Shah, S. A. R., Anwar, M. K., Raza, A., Ullah,
permeabilidad reducida, lo que indica una permeabilidad M. K., & Iqbal, F. (2021). Utilization of Polymer
relativamente baja en comparación con la muestra control que Concrete Composites for a Circular Economy: A
obtuvo 867,09 g/m2. Del mismo modo, la entrada de los Comparative Review for Assessment of Recycling and
agentes que aceleran la velocidad de deterioro en el concreto Waste Utilization. Polymers, 13(13), 2135.
hidráulico. La confirmación de esta conclusión se respalda con https://doi.org/10.3390/polym13132135
los datos de los resultados mostrados en la Tabla 9 del análisis Amin, M., Agwa, I. S., Mashaan, N., Mahmood, S., & Abd-
estadístico Anova-Varianza, que alcanzó un p-valor de 0,045, Elrahman, M. H. (2023). Investigation of the Physical
inferior al nivel de significación predeterminado de alfa 0,05. Mechanical Properties and Durability of Sustainable
Así, se establece que el efecto sobre las propiedades del Ultra-High Performance Concrete with Recycled Waste
concreto hidráulico con el uso de CC y RC es favorable. Glass. Sustainability, 15(4), 3085.
https://doi.org/10.3390/su15043085
Se establece que la resistencia a la compresión con ceniza de ASTM C33/CMMM-18. (2018, April 20). Standard
cartón (CC) y Cerámico Reciclado (RC) de CC 3,0 % + RC Specification for Concrete Aggregates.
2,0 % alcanzó la mayor resistencia con un valor de 323,34 https://www.astm.org/c0033_c0033m-18.html
kg/cm2 ± 6,85 kg/cm2 después de un periodo de curado de 28 ASTM C39/C39M-21. (2021, March 9). Standard Test
días, superando a la muestra control que obtuvo 282,73 kg/cm2 Method for Compressive Strength of Cylindrical
± 1,91 kg/cm2 y a la norma ASTM C-39 con diseño establecido Concrete Specimens.
de 280 kg/cm2. Asimismo, cuando se utiliza mayor cantidad https://www.astm.org/c0039_c0039m-21.html
de materiales de CC y RC, la resistencia a la compresión tiende ASTM C78/C78M. (2022, March 30). Standard Test for
a disminuir. Estos resultados se respaldan con los datos Flexural Strength of Concrete (Using Simple Beam with
mostrados en la Tabla 10 del análisis estadístico Anova- Third-Point Loading).
Varianza, que alcanzó un p-valor de 0,000, inferior al nivel de https://www.astm.org/c0078_c0078m-22.html
significación predeterminado de alfa 0,05 para grupos a las ASTM C1585-20. (2020, September 22). Standard Test
edades de 7, 14 y 28 días. En consecuencia, se determina que Method for Measurement of Rate of Absorption of Water
el uso de materiales reciclados de CC y RC, tiene efectos by Hydraulic-Cement Concretes.
significativos en las propiedades del concreto. https://www.astm.org/c1585-20.html
ASTM D854-14. (2023, May 15). Standard Test Methods for
Se concluye que la resistencia a la flexión con ceniza de cartón Specific Gravity of Soil Solids by Water Pycnometer
(CC) y cerámico reciclado (RC) de CC 3,0 % y RC 2,0 % (Withdrawn 2023). https://www.astm.org/d0854-
alcanzó la mayor resistencia de 35,76 kg/cm2 ± 0,31 kg/cm2 14.html
después de un periodo de curado de 28 días, superando la Barrios, A. M., Vega, D. F., Martínez, P. S., Atanes-Sánchez,
muestra control que obtuvo 32,93 kg/cm2 ± 0,29 kg/cm2 y a la E., & Fernández, C. M. (2021). Study of the properties
norma ASTM C78 con una resistencia requerida de 28 kg/cm2. of lime and cement mortars made from recycled ceramic
De esta forma, cuando se utiliza mayor cantidad de materiales aggregate and reinforced with fibers. Journal of Building
revistapolitecnica.epn.edu.ec
ISSN: 1390-0129
eISSN: 2477-8990
Volumen 54, Nro. 1, Agosto · Octubre 2024
revistapolitecnica.epn.edu.ec
www.epn.edu.ec