Elaboración de un modelo evolutivo híbrido de algoritmos genéticos y redes neuronales artificiales para dosificaciones de mezclas de concreto reforzadas con fibras metálicas

Development of a Hybrid Evolutionary Model of Genetic Algorithms and Artificial Neural Networks for Metal Fiber and Reinforced Concrete Mixture Dosage

Luis Octavio González Salcedo ⁽¹⁾, Aydee Patricia Guerrero Zúñiga ⁽²⁾, Silvio Delvasto Arjona ⁽³⁾, Adrián Luis Ernesto Will ⁽⁴⁾

⁽¹⁾ Doctor en Ingeniería de Materiales. Profesor Asociado, Facultad de Ingeniería y Administración, Universidad Nacional de Colombia Sede Palmira. Palmira, Colombia. Grupo de Investigación en Materiales y Medio Ambiente. logonzalezsa@unal.edu.co.

⁽²⁾ Doctor of Philosophy in Civil Engineering. Profesora Titular, Facultad de Ingeniería, Universidad del Valle. Cali, Colombia. Grupo de Investigación en Ingeniería Sísmica, Eólica, Geotécnica y Estructural. patricia.guerrero@correounivalle.edu.co.

⁽³⁾ Doctor en Nuevos Materiales y sus Tecnologías de Fabricación, Profesor Titular, Facultad de Ingeniería, Universidad del Valle, Cali, Colombia. Grupo de Materiales Compuestos. silvio.delvasto@correounivalle.edu.co.

⁽⁴⁾ Doctor en Matemáticas. Profesor Adjunto, Facultad de Ciencias Exactas y Tecnología - Universidad Nacional de Tucumán, Facultad Regional de Tucumán - Universidad Tecnológica Nacional. Tucumán, Argentina. Grupo de Investigación de Tecnologías Avanzadas de Tucumán. adrian.will.01@gmail.com.

Recibido 15 de agosto de 2015. Modificado 30 de noviembre de 2015. Aprobado 9 de diciembre de 2015.

DOI:http://dx.doi.org/10.16924/riua.v0i43.874

Algoritmos genéticos, dosificaciones de mezclas de concreto, modelo evolutivo, redes neuronales artificiales.

Resumen

En un ambiente computacional, un modelo evolutivo es desarrollado para proponer dosificaciones de mezclas de concretos reforzados con fibras metálicas, para solicitudes de resistencia a la compresión. El modelo es híbrido al incluir un sistema de dosificación basado en algoritmos genéticos y un sistema de estimación de propiedades basado en redes neuronales artificiales. Los resultados obtenidos son comparados con un sistema de dosificaciones reportado experimentalmente. Estas comparaciones muestran una aproximación en el proceso de simulación. Además, por las características del modelo, se considera de aporte en la tecnología del concreto.

Key words

Genetic Algorithms, Concrete Mixture Dosages, Evolutionary Model, Artificial Neural Networks.

Abstract

An evolutionary model is developed in a computing environment to propose metal fiber reinforced concrete mixture dosages for compressive strength applications. The model is hybrid as it includes both a dosage system based on genetic algorithm and a properties prediction system based on artificial neural networks. The results obtained are compared with experimentally reported dosages set, and the comparisons show an approximation in the simulation process. Given the characteristics of the model, it is considered a contribution to concrete technology.

Introducción

Esta investigación aborda la estimación de dosificaciones de mezclas de concreto reforzado con fibras metálicas usando Algoritmos Evolutivos, basados en la técnica de los Algoritmos Genéticos (AG). Se requiere en algunos concretos de un comportamiento de desempeño correspondiente a un ambiente aplicado, y hoy en día, el progreso de la tecnología del concreto permite hacer que el concreto encuentre estos requerimientos. Sin embargo, aunque hay manera de optimizar matemáticamente el problema en muchos criterios, los cuales están representados por funciones objetivos, el problema de decisión es multivariado y una amplia variedad de dosificaciones como soluciones son posibles (Siddique & Khan, 2011).

Los AG son algoritmos de optimización, aprendizaje y búsqueda basados en el mecanismo de la selección natural y la genética natural, y se ha comprobado su efectividad en aplicaciones en el campo de la ingeniería, especialmente para problemas de combinación. Dado que la dosificación de las mezclas de concreto no es solamente un tipo de problema de combinación, sino también de optimización combinatoria y compleja que involucra variables numéricas y cualitativas, es entonces una buena razón aplicar en este caso los AG, a la par de que la dosificación del concreto posiblemente induzca a la aparición de cálculos en funciones difíciles de evaluar o de altos errores de medición (Noguchi et al, 2003). Setién et al. (2006) ratifican el problema en la determinación de las propiedades del concreto a partir de las dosificaciones de la mezcla, y sugieren la introducción de técnicas basadas en el uso de Redes Neuronales Artificiales (RNA), dada la complejidad e incertidumbre de los factores allí involucrados.

Metodología

Para el desarrollo de la investigación se elabora un Modelo Evolutivo Híbrido, denominado así en razón de la estructura del Modelo Evolutivo basado en AG que incorpora un Modelo Neuronal basado en RNA. El AG se usa como un Sistema de Dosificación, mientras la RNA, diseñada como un Sistema de Estimación, en este caso como componente de la estructura del sistema de dosificación, se usa como Función de Estimación. El procedimiento metodológico, la definición de sus elementos y sus operadores genéticos, las bases fundamentales, y los detalles de construcción de un Modelo Evolutivo basado en AG son determinados por Goldberg (1989). En su funcionamiento típico, los AG requieren de elementos iniciales como una población posible y potencial de soluciones, la aplicación de operaciones genéticas, la determinación de una función de evaluación o de adaptación al medio, y el mecanismo de evaluación a esa adaptación o supervivencia mediante una penalización. En adición, se requiere incorporar operaciones de restricción y ajuste, para aproximar el Modelo Evolutivo al problema que se quiere simular.

Población inicial

Se plantea un conjunto que describa el sistema de dosificación en el diseño de mezclas de concreto reforzado con fibras metálicas. La dosificación del concreto consiste en la determinación de las proporciones en que han de mezclarse los componentes que lo constituyen. Dado que el concreto está constituido por diversos materiales, no existe una única dosificación posible para obtener una resistencia a la compresión f'c (se toma esta propiedad mecánica, en razón al cumplimiento en las especificaciones de construcción), sino que hay un rango de dosificaciones posibles que responden a esa misma propiedad mecánica. Por esta razón, la población inicial se define a partir de dos aspectos: las variables a considerar dentro de la dosificación y los límites del rango para cada variable. Con respecto al primer aspecto, las variables consideradas como relevantes del problema de dosificación son mostradas en la Tabla 1.

Se crea una matriz de límites para cada variable considerada en la dosificación, en la cual se definen el valor mínimo y el valor máximo de cada una de ellas, basados en trabajos realizados con anterioridad sobre optimización de componentes de dosificaciones de concreto (Will et al., 2008; Yeh, 2009; Chen & Wang, 2010), y la recopilación de 601 diseños de mezclas en un sistema de dosificación de datos experimentales reportados en la literatura especializada (González, 2014). Se muestra en la Tabla 2, la matriz de límites para la población inicial en concretos para cuatro rangos de resistencia a la compresión f'c.

Para la aplicación de los operadores genéticos se requiere la definición de una función distancia, con el fin de identificar en los posibles individuos seleccionados (soluciones seleccionadas) si existen diferencias significativas entre ellas desde el punto de vista de sus componentes genéticos, es decir, la configuración del "cromosoma". La función distancia se muestra en la Tabla 3.

Función de Estimación

Como función de estimación, se usa un sistema de predicción basado en un modelo neuronal artificial, entrenado para la estimación de la resistencia a la compresión, dada una dosificación de acuerdo con las variables de entrada consideradas en la Tabla 1. El modelo neuronal artificial consistió de una RNA, con arquitectura multicapa, que responde a un problema de decisión de dimensión 30 (resultante de las 19 variables mostradas en la Tabla 1), del cual considera 12 variables referidas a las cantidades de los materiales y 18 variables referidas a atributos o cualidades de los materiales, y cuyo desempeño se evaluó con el coeficiente de correlación lineal R. Para el entrenamiento de la RNA, se construyó un conjunto de entrenamiento conformado por 611 diseños de mezclas conducentes a una resistencia a la compresión f'c, recopilados de la literatura especializada (González, 2014). Además, se emplearon diferentes arquitecturas (capas ocultas y números de neuronas en las capas ocultas) y diferentes algoritmos de aprendizaje feedforward - blackpropagation. Como lenguaje de programación del algoritmo de la RNA se usa MATLAB® & SIMULINK® para ambiente WINDOWS® (The MathWorks Inc, 2008), y la codificación se ha basado en el Neural Network Toolbox del lenguaje mencionado (Beale et al., 2012). La RNA seleccionada como sistema de estimación correspondió a una red multicapa con las características mostradas en la Tabla 4. El Factor de Correlación R es mostrado en la Figura 1.

Restricciones y funciones de ajuste

El objetivo de los AG en esta investigación es proponer la dosificación en mezclas de concreto reforzados con fibras metálicas que satisfarán un valor de resistencia a la compresión requerida. De igual manera, se respeta un número importante de restricciones propias del proceso de diseño de la mezcla, como se hace en los procedimientos protocolarios de esta actividad. Se plantean funciones para realizar las siguientes restricciones y funciones de ajustes (ASOCRETO, 2010): Control Volumétrico (1,00 m³), Relación Agua Total / Cementante (0,40-0,60), limitación en el uso de aditivos químicos, limitación en el reemplazo del cemento por materiales cementantes y limitación en el tamaño máximo del agregado grueso debido al tamaño de la fibra metálica empleada.

Función Objetivo y Función de Penalización

Se usa una función de evaluación para múltiples criterios: la Resistencia a la Compresión, el Asentamiento del Cono de Abrams, el Control Volumétrico, la Relación Agua Total/Cementante, la limitación en el uso de los aditivos químicos, la limitación en el reemplazo del cemento por materiales cementantes, el tamaño máximo del agregado como función del diámetro de la fibra, y el Asentamiento del Cono de Abrams como función de la cantidad de agua, del tamaño máximo del agregado y del tipo del agregado, esquematizado en la representación de la Ecuación 1:

Operadores Genéticos

Como operadores de selección o reproducción (Gil, 2006) se usan Ruleta y Muestreo Universal Estocástico; como operadores de cruce (Gil, 2006) se usan Cruce Uniforme, Heurístico y Aritmético; por Recombinación Intermedia y para Permutación se usa Mapeamiento Parcial y de Ciclo, y como operador de mutación (Davis, 1991) se usa No Uniforme y Uniforme.

ARQUITECTURA DEL MODELO EVOLUTIVO

En el laboratorio computacional del Grupo de Aplicaciones de Matemática e Inteligencia Artificial de la Facultad de Ciencias Exactas y Tecnología - Universidad Nacional de Tucumán, localizada en San Miguel de Tucumán (Argentina), se elaboró un modelo básico de Algoritmos Genéticos, para el problema de dosificación en mezclas de concreto no reforzadas, cuyos resultados fueron analizados usando indicadores de comparación entre un sistema de dosificación de datos experimentales reportados y el de los resultados obtenidos por el modelo mencionado (González et al., 2009). Este modelo fue modificado, complementado y adaptado al problema de dosificaciones de mezclas en concretos reforzados con fibras, en las instalaciones de la Universidad del Valle, localizada en Santiago de Cali (Colombia); las adaptaciones correspondieron a la incorporación de las características específicas relacionadas con la generación de la población inicial, la función de estimación, las restricciones y las funciones de ajustes, la función objetivo y de penalización, y los operadores genéticos, mencionados con anterioridad. Para el modelo evolutivo basado en algoritmos genéticos se elaboraron un conjunto de programas de computador usando como lenguaje de programación para su escritura MATLAB® para ambiente Windows® (The MathWorks Inc, 2008); la codificación se basó en el Genetic Algorithm Toolbox del lenguaje mencionado (The MathWorks Inc, 2004). El conjunto de programas computacionales se estructuró en una compleja arquitectura consistente de siete módulos de primer orden, que comprende un programa principal para dar inicio a los algoritmos genéticos y seis programas auxiliares; veinticinco módulos de segundo orden, dispuestos a manera de librería para el programa principal; los módulos del primer orden, cuya funcionalidad es realizar todas las operaciones simuladas computacionalmente de un proceso de evolución.

Resultados y discusión

El modelo evolutivo basado en algoritmos genéticos se ejecutó para solucionar dosificaciones de mezclas de concreto reforzados con fibras, con requerimientos de resistencia a la compresión de 21, 30, 40, 50, 60, 70, 80, 90 y 100 MPa. A manera de ejemplo, en este documento se muestran algunos resultados para la simulación realizada, considerando en el diseño de la mezcla cemento Tipo I: superplastificante como agente retenedor de agua, arena procedente de río, grava procedente de cantera y perfil litológico I y fibra metálica tipo Hooke-end de longitud 30 mm y relación de aspecto 80. Para cada resistencia mencionada, se consideraron tres tipos de adiciones minerales a saber: en la primera, se consideró cenizas volantes, en la segunda se consideraron cenizas volantes y escoria, y en la tercera se consideró humo de sílice (Figura 2).

Para el análisis de los resultados, se presentan algunos indicadores de comparación entre el sistema de dosificación de los datos experimentales reportados y el de los resultados obtenidos por el Modelo Evolutivo (Figura 3). De la misma manera, se toman como indicadores de comparación relaciones entre los materiales cementantes y el volumen de fibra empleados en las dosificaciones (Figura 4).

Los resultados de estas comparaciones muestran que los parámetros del material mencionados tienen una similitud de tendencia para ambos sistemas, de la siguiente manera:

• Con respecto al contenido de material cementante (total de cemento, humo de sílice, cenizas volantes y escoria) (Fig. 3), se observa en ambos sistemas de dosificación que este parámetro se incrementa para mayores requerimientos de resistencia, coincidiendo con Will et al. (2008), y Yeh (2009).
• Con respecto a la relación del contenido de agua total/contenido del material cementante (W/B) (Fig. 3), se observa en ambos sistemas de dosificación que este parámetro se comporta con la relación clásica f'c vs W/B (ASOCRETO, 2010), en la cual, a medida que se incrementa el valor del parámetro, la resistencia a la compresión disminuye, coincidiendo con Will et al. (2008), y Yeh (2009).
• Con respecto a las distribuciones de frecuencia para la relación entre los contenidos de adiciones minerales (contenidos de humo de sílice, cenizas volantes y escoria) total y específicos y el contenido del material cementante (Fig. 4), se observa en ambos sistemas de dosificación que este parámetro tiene tendencias similares. En particular se hacen las siguientes apreciaciones:
• El sistema de dosificación propuesto por el AG se caracteriza por emplear reemplazos del cemento por cantidades controladas de adiciones minerales, de una manera similar a como se realiza en el sistema de dosificación reportado.
• El sistema de dosificaciones propuesto por el AG sigue una tendencia similar en cuanto a los porcentajes de reemplazo para cada adición mineral, mostrando un número importante de dosificaciones para reemplazos entre el 12-15% para humo de sílice, el 20-22% para cenizas volantes, y alrededor del 25% para escoria, lo cual concuerda con los valores máximos sugeridos en la literatura, para cada caso (Siddique & Khan, 2011).
• El sistema de dosificaciones propuesto coincide en reemplazos específicos reportados en otras investigaciones, como se detalla en la Figura 4, coincidiendo con Ujhelyi & Ibrahim, 1991, Hooton 1993, Joshi et al., 1993, Babu & Kumar, 2000, Wainwright & Rey, 2000, Wong & Razak, 2005, Koksal et al., 2008, y Yeh, 2009.

CONCLUSIONES

El modelo genera un conjunto de posibles dosificaciones, las cuales son evaluadas en la estimación de la resistencia a la compresión y categorizadas de acuerdo con el resultado obtenido, para ser escogidas como las mejores propuestas de dosificación acorde a la propiedad mecánica requerida. Las propuestas son sometidas a modificaciones, usando operaciones genéticas, para encontrar otras dosificaciones que cumplan con el mismo objetivo de alcanzar la propiedad mecánica. Los resultados muestran una alta cantidad de propuestas de dosificaciones de mezclas de concreto reforzadas con fibras, las cuales evidencian indicadores de similitud en la tendencia de diversos parámetros del material, así como una similitud en el contenido de materiales relevantes de comparación como el contenido del material cementante y del contenido de agua total, al ser comparada con un sistema de dosificaciones experimentales reportadas. Los resultados del modelo evolutivo permiten trazar las siguientes conclusiones:

• Se abre una futura agenda de investigación con las propuestas de dosificaciones de las mezclas, con la cual puede realizarse un análisis interesante para conformar clasificaciones como concretos con adiciones simples, binarias y terciarias, entre otras.
• El modelo evolutivo puede ser ajustado para usar en sus funciones objetivo y de penalización (ajustando fPenalización en la Ec.1) y la incorporación de las correspondientes funciones de estimación. Además, la estimación de otras propiedades de tal manera que la obtención de propuestas de dosificaciones que se realicen para objetivos de criterios múltiples, es decir, que satisfagan no solo valores de resistencia a la compresión, sino también otras propiedades tanto en estado fresco (asentamientos) como en estado endurecido (propiedades mecánicas y de durabilidad), abre también una futura agenda de interés para el campo de la tecnología del concreto.
• La consideración de variables tanto de tipo cuantitativa como cualitativa, como componentes de la mezcla, en similitud a las variables de entrada en el entrenamiento de las redes neuronales artificiales, posibilitan el diseño de un modelo híbrido, conformado por un modelo evolutivo cuya función de estimación es realizada por el modelo neuronal.

Se concluye que el modelo evolutivo desarrollado y expuesto en la presente investigación, en razón a sus características, constituye un aporte importante en la tecnología del concreto, específicamente en la dosificación de mezclas de concretos reforzados con fibras metálicas.

Referencias

Asociación Colombiana de Productores de Concreto -ASOCRETO. (2010). Tecnología y Propiedades. Colección Básica del Concreto. Tomo 1. Bogotá: Instituto del Concreto, ASOCRETO.         [ Links ]

Babu, K.G. & Kumar, V.S.R. (2000). Efficiency of GGBS in concrete. Cement Concrete Research, 30(7), 1031-1036. doi:10.1016/S0008-8846(00)00271-4.         [ Links ]

Beale, M.H., Hagan, M.T. & Demuth, H.B. (2012). Neural Network ToolboxTM, User's Guide. Natick: The MathWorks Inc.         [ Links ]

Chen, L. & Wang, T-S. (2010). Modeling Strength of High-Performance Concrete Using and Improved Grammatical Evolution Combined with Macrogenetic Algorithm. Journal of Computing in Civil Engineering, 24(3), 281-288. doi: 10.1061/(ASCE)CP.1943-5487.0000031.         [ Links ]

Davis, L.D. (Ed). (1991). The Handbook of Genetic Algorithms. New York: Van Nostrand Reinhold.         [ Links ]

Gil, L. N. (2006). Algoritmos Genéticos. Universidad Nacional de Colombia Sede Medellín. Medellín, Noviembre, 65p.         [ Links ]

Goldberg, D.E. (1989). Genetic Algorithms in Search: Optimization and Learning. USA: Addison-Wesley Publishing Co., Reading.         [ Links ]

González S., L.O.; (2014). Diseño de mezclas de concreto reforzado con fibras metálicas y de polipropileno, usando inteligencia artificial. (Tesis inédita de doctorado) Universidad del Valle, Santiago de Cali, Colombia.         [ Links ]

González, S., L.O., Torres A., J. & WILL, A.L.E. (2009, octubre). Evaluación de diseños de mezclas de concreto obtenidos usando Inteligencia Artificial basada en Algoritmos Genéticos. Ponencia presentada en V Congreso Internacional de Materiales, Memorias. Santiago de Cali, Colombia.         [ Links ]

Holschemacher, K., Müller, T. & Fischer, A. (2006, august). Effect of fibre type on properties of steel fibre reinforced concrete. Ponencia presentada en The Tenth East Asia-Pacific Conference on Structural Engineering & Construction (EASEC-10). Asian Institute of Technology. Bangkok, Thailand.         [ Links ]

Hooton, R.D. (1993). Influence of silica fume replacement of cement on physical properties and resistance to sulfate attack freezing and thawing, and alkali-silica reactivity. ACI Materials Journal, 90(2), 143-152.         [ Links ]

Joshi, R.C., Lohtia, R.P. & Salam, M.A. (1993). High strength concrete with high volumes of Canadian sub-bituminous coal ash. Ponencia presentada en Third International Symposium on Utilization of High Strength Concrete. Lillachammer, Norway.         [ Links ]

Köksal, F., Altun,, F., Yigit, I. & Súahin, Y. (2008). Combined effect of silica fume and Steel fibre on the mechanical properties of high strength concretes. Construction Building Materials, 22(8), 1874-1880. doi:10.1016/j.conbuildmat.2007.04.017.         [ Links ]

Lourakis, M.A. (2005). A brief description for the Levenberg-Marquardt Algorithm implemented by levmar. Heraklion: Institute of Computer Science, Foundation for Research and Technology - Hellas (FORTH).         [ Links ]

Martínez P., D. (2008). Hormigones de altas prestaciones. (Proyecto Fin de Carrera Ingeniería Técnica Civil: Arquitectura Técnica), Universidad Politécnica de Cartagena, Cartagena, España.         [ Links ]

Noguchi, T., Maruyama, I. & Kanematsu, M. (2003, may). Performance based design for concrete mixture with multi-optimizing Genetic Algorithm. Ponencia presentada en The 11th International Congress on the Chemistry of Cement (ICCC): "Cement's contribution to the development in the 21st Century". Cement & Concrete Institute. Durban, South Africa.         [ Links ]

Özgan, E. & Öztürk, A. (2007). Optimization of the Hardened Concrete Properties with GA and LP. Journal of Applied Sciences, 7(24), 3918-3926.         [ Links ]

Sánchez D.G., D. (2000). Tecnología del Concreto y del Mortero. 2ª edición. Bogotá: Biblioteca de la Construcción Bhandar Editores.         [ Links ]

Siddique, R. (2003). Effect of fine aggregate replacement with class F fly ash on the mechanical properties of concrete. Cement Concrete Research, 33(4), 539-547. doi:10.1016/S0008-8846(02)01000-1.         [ Links ]

Siddique,, R. & Khan, M.I.; (2011). Supplementary Cementing Materials. New York: Springer. Engineering Materials Collection.         [ Links ]

Sietén, J., Carrascal, I.A., Figueroa, J.F. & Polanco, J.A. (2003). Aplicación de una red neuronal artificial a la dosificación de hormigones premezclados. Materiales de Construcción, 53(270), 5-19. doi: 10.3989/mc.2003.v53.i270.270.         [ Links ]

The Mathworks Inc. (2004). Genetic Algorithm and Direct Search Toolbox for use with MATLAB®, User's Guide. Version 1. Natick: The MathWorks Inc.         [ Links ]

The Mathworks Inc. (2008). MATLAB® & SIMULINK® Release 2008a, Installation Guide for Windows. Natick: The MathWorks Inc.         [ Links ]

Ujhelyi, J.E. & Ibrahim, A.J. (1991). Hot weather concreting with hydraulic additives. Cement Concrete Research, 21(2-3), 345-354. doi:10.1016/0008-8846(91)90015-A.         [ Links ]

Wainwright, P.J. & Rey, N. (2000). The influence of ground granulated blast furnace slag (GGBS) additions and time delay on the bleeding of concrete. Cement Concrete Composites, 22(4), 253-257. doi:10.1016/S0958-9465(00)00024-X.         [ Links ]

Will, A., Folino, P., Bálzamo, H. & Etse, G. (2008). Acerca de la conveniencia del uso de algoritmos genéticos como herramienta para la dosificación de hormigones. Mecánica Computacional, 27(12), 865-880.         [ Links ]

Wong, H.S. & Razak, H.A. (2005). Efficiency of calcined kaolin and silica fume as cement replacement material for strength performance. Cement Concrete Research, 35(4), 696-702. doi:10.1016/j.cemconres.2004.05.051.         [ Links ]

Yeh, I.C. (2009). Optimization of concrete mix proportioning using a flattened simplex-centroid mixture design and neural networks. Engineering with Computers, 25(2), 179-190. doi: 10.1007/s00366-008-0113-2.         [ Links ]