Este artículo se ocupa del diseño automático de tableros de puentes de vigas artesa pretensadas prefabricadas de coste mínimo, empleando para ello un algoritmo memético híbrido que combina la búsqueda poblacional de soluciones mediante algoritmos genéticos y una búsqueda por entornos variable. Este algoritmo se aplica a un puente formado por 2 vigas isostáticas con luces entre apoyos de entre 20 y 40m y una losa de 12m de ancho. La estructura analizada consta de 40 variables discretas. El módulo de la evaluación considera los estados límite último y de servicio que se aplican habitualmente para estas estructuras: flexión, cortante, torsor, fisuración, flechas, etc. El uso del algoritmo memético requiere previamente su calibración. Cada una de las heurísticas se procesa 12 veces, obteniéndose información estadística sobre el valor mínimo, el medio y las desviaciones. El estudio paramétrico muestra una buena correlación del coste, del número de torones y de las cuantías de acero pasivo y hormigón de la viga con la luz. Se han comprobado ahorros de entre el 8 y el 50% respecto a otras estructuras realmente ejecutadas. El procedimiento presentado permite la aplicación práctica al diseño real y su adaptación al proceso de prefabricación.
This paper deals with the minimum cost automatic design of precast bridge decks made of U-beams and an upper slab. It uses a hybrid memetic algorithm that combines the population search of solutions by genetic algorithms and a search by variable neighborhood. This algorithm is applied to a bridge made of two isostatic U-beams of 20-40m of span and a width of 12m. This example has 40 discrete variables. The evaluation module takes into account the service and ultimate limit states usually considered for these structures, i.e. flexure, shear, torsion, cracking, deflections, etc. The use of the memetic algorithm requires its previous calibration. Each of the heuristics is run 12 times, obtaining information about the minimum and average values, as well as the scatter. The parametric study showed a good correlation for the cost, the number of strands and the steel and concrete quantities with the span length. Savings have been found between 8 and 50% compared to other structures really executed. The presented procedure allows the practical application to the real design and its adaptation to the precast process.
La construcción con hormigón prefabricado presenta claras ventajas económicas cuando se fabrican en taller piezas en grandes series. El ahorro en material y en mano de obra, la elevada calidad en el producto y el rápido montaje son razones que justifican, por sí solas, el uso de la construcción prefabricada. Sin embargo, tal y como indica Yee [1], hoy en día existen motivos adicionales basados en beneficios sociales y medioambientales que justifican la adopción de la tecnología del hormigón prefabricado. A este respecto, Billington [2] proporciona una perspectiva histórica del desarrollo del hormigón pretensado (HP). Asimismo, los proyectistas han tomado buena nota de las ventajas del prefabricado cuando se trata de construir puentes con luces moderadas, de 10 a 40m [3]. En estos casos, la disminución del peso resulta fundamental para reducir los costes de elevación y transporte de las piezas. En este contexto, la optimización estructural del coste necesario para construir un puente de vigas prefabricadas constituye un área de gran interés, especialmente cuando se realizan grandes series de piezas.
Si bien los aspectos básicos relacionados con la optimización matemática se establecieron en los siglosxviii y xix con los trabajos de Lagrange o Euler, hay que esperar hasta los años cuarenta del sigloxx para que Kantorovich y Dantzing desarrollaran plenamente los principios de la programación matemática. Se puede consultar una revisión de las técnicas de optimización más importantes en la recopilación de Caballero y Grossmann [4]. Es a partir de la revolución informática de los años setenta cuando estas herramientas empezaron a ser empleadas de forma habitual en numerosas aplicaciones en las ciencias, las ingenierías y los negocios. Sin embargo, el progreso de técnicas de optimización que no requieran derivadas y que se generen a través de reglas heurísticas ha supuesto una auténtica revolución en el campo de la optimización de los problemas reales. En efecto, los métodos aproximados pueden utilizarse allí donde el elevado número de variables en juego impide la resolución de los problemas, mediante la programación matemática, en un tiempo de cálculo razonable. A estos algoritmos de optimización aproximada, cuando su uso no está restringido a un solo tipo de problemas, la comunidad científica en el ámbito de la inteligencia artificial y de la investigación operativa les ha dado el nombre de metaheurísticas [5]. Este grupo incluye una amplia variedad de procedimientos inspirados en algunos fenómenos naturales, tales como los algoritmos genéticos [6], el recocido simulado [7] o la optimización por movimientos de enjambres [8]. Liao et al. [9] presentan una revisión reciente de la aplicación de los métodos heurísticos en el campo de la gestión del proyecto y de la construcción.
En relación con la optimización de las estructuras, si bien la información más antigua hay que buscarla en el sigloxv con los trabajos de Leonardo da Vinci y de Galileo Galilei sobre la disminución del peso de estructuras de madera, hay que esperar al sigloxix, con Maxwell y Levy, y a comienzos del sigloxx, con Mitchell, para ver las primeras aportaciones en el diseño de mínimo peso de estructuras de arcos y cerchas metálicas. En 1994, Cohn y Dinovitzer [10] realizaron una amplia revisión de los métodos empleados en la optimización de estructuras, comprobando que la inmensa mayoría de las investigaciones llevadas a cabo hasta entonces se basaban en la programación matemática y en problemas más bien teóricos, con una preponderancia abrumadora de las estructuras metálicas frente a las estructuras de hormigón. Así, la aplicación de métodos heurísticos a la ingeniería estructural se remonta a los años setenta y ochenta [11–13], siendo la computación evolutiva —y en especial los algoritmos genéticos— los métodos que más se han utilizado. La revisión de Kicinger et al. [14] proporciona un completo estado de la cuestión de los métodos evolutivos aplicados al diseño estructural. Por otro lado, nuestro grupo de investigación ha presentado trabajos recientes de diseño automático y optimización de estructuras de hormigón armado con algoritmos genéticos [15] y con otras técnicas heurísticas [16–21], así como trabajos de optimización con HP [22,23].
En relación con la optimización de puentes, la revisión mencionada de Cohn y Dinovitzer [10] ya apuntaba la gran escasez de artículos publicados en esta materia. El diseño óptimo de vigas pretensadas, en especial la disposición de los tendones, es un problema clásico planteado desde hace años [24,25]. Aparicio et al. [26] presentaron un sistema de diseño asistido por ordenador de puentes de HP para carreteras, identificando cuáles eran las tipologías estructurales más eficaces. Hassanain y Loov [27] presentan una revisión del estado de la cuestión de las técnicas de optimización de puentes de hormigón. Sin embargo, tal y como apuntan Hernández et al. [28], existe cierto vacío en la investigación que se ocupe específicamente de la optimización y el diseño completo de los puentes reales. En este sentido, el trabajo de Martí y González-Vidosa [22] implicó el diseño integral y automatizado de una pasarela para peatones formada por una viga prefabricada de HP.
Siguiendo esta línea de investigación, el presente artículo se centra en el diseño automatizado de puentes de vigas artesa de HP prefabricadas empleados como pasos superiores sobre vías de comunicación. Las luces vienen impuestas por las dimensiones de la vía inferior, con rangos habituales que oscilan entre los 20 y los 40m. Estos puentes consisten en vigas de HP con forma de U con losa superior colaborante (figs. 1 y 2) y un tablero de hormigón, parcialmente prefabricado o construido in situ. Esta tipología cuenta a su favor, entre otras, con las ventajas derivadas de la prefabricación, como por ejemplo la construcción industrializada, los moldes reutilizables, los plazos reducidos de ejecución en obra y la baja interferencia con el tráfico inferior. La solución de viga en U permite eliminar completamente los poco agraciados cabezales sobre pila de los tableros de viga en doble T. La estructura que aquí es objeto de optimización está compuesta por 2 vigas de HP prefabricadas que integran en la parte superior una losa de HA para el tráfico de vehículos de 12m de anchura, siendo su tablero isostático en sentido longitudinal (fig. 2). El método para realizar la optimización se ha basado en el cálculo del coste de cada una de las soluciones evaluadas, en función de las variables geométricas, los materiales y los refuerzos activo y pasivo. A continuación se ha implementado un módulo que valora todos los estados límites relevantes, y posteriormente se ha desarrollado, de modo específico para este trabajo, un algoritmo memético híbrido capaz de encontrar soluciones optimizadas en coste. El artículo, tras realizar el planteamiento del problema de optimización, define el algoritmo citado, presenta los resultados obtenidos y señala las principales conclusiones.
2Definición del problema de optimización
El problema planteado reside en la optimización económica de un diseño estructural de HP. Se ocupa de la minimización de la función objetivo F de la expresión (1), satisfaciendo también las restricciones formuladas en la expresión (2).
Los valores x1, x2,…, xn son las variables de diseño del problema, que pueden tomar uno de los valores discretos indicados en la expresión (3). Obsérvese que la función objetivo en la expresión (1) es la suma de los precios unitarios multiplicados por las mediciones de las unidades de obra (hormigón, acero, encofrado, transporte, colocación, etc.). La expresión (2) indica las restricciones geométricas y de constructibilidad, así como todos los estados últimos y de servicio que la estructura ha de cumplir. Las soluciones que satisfacen las restricciones en la expresión (2) se denominan factibles, y las que no, soluciones no factibles. La definición del problema se basa en la tesis doctoral de Martí [29].
2.1Variables de diseño y parámetrosLa solución completa del puente de vigas artesa objeto de análisis puede definirse con 40 variables de diseño. Estas variables son discretas para facilitar la construcción efectiva de la estructura real optimizada. Se incluyen entre ellas 8 variables geométricas (figs. 2 y 3), que toman valores escalonados de centímetro en centímetro. El canto de la viga h1 oscila entre 0,50m hasta 1/17 de la luz para limitar la esbeltez mínima y para permitir el transporte de la viga por carretera. El ancho del ala inferior de la viga b1 puede variar entre 0,50 y 2,00m, mientras su espesor e1 toma valores comprendidos entre 0,15 a 0,50m. El ancho b3 de las alas superiores de la viga puede variar desde 0,15 hasta 1,00m. Tanto el espesor e2 de las almas como el espesor e3 de las alas pueden tomar valores comprendidos entre 0,10 y 0,50m. El espesor de la losa e4 varía entre 0,12 hasta 0,47m. Finalmente, la separación entre vigas Sv puede comprender valores entre 3,96 hasta 6,96m.
Las variables que definen la resistencia característica de los hormigones toman valores de entre 25MPa a 40MPa para la losa y desde 35MPa hasta 50MPa para las 2 vigas, en escalones de 5MPa. La armadura de pretensado puede definirse mediante 4 variables: (1) el número de torones en las alas superiores, con un máximo de 10; (2) el número de torones dispuesto entre las primeras, segundas y terceras capas del ala inferior, con un máximo de 98; (3) el número de tramos con fundas dispuestas en la segunda capa y (4) el número de tramos con fundas dispuestas en la tercera capa (fig. 3). El enfundado, si existe, comienza en un extremo de la viga y se prolonga por tramos consecutivos, considerándose cada tramo como la distancia entre 2 secciones de cálculo del modelo estructural. Por último, son necesarias 23 variables para definir la disposición del armado pasivo, tanto para la viga como para la losa superior (fig. 2). La armadura pasiva transversal se mantiene uniforme cada 2 tramos en que se encuentra modelizada la estructura.
El conjunto de las 40 variables de diseño y sus valores posibles permiten un número total de combinaciones realmente elevado, del orden de 3,50×1051 soluciones. Este desorbitado espacio de soluciones es el que justifica la adopción de algoritmos heurísticos para encontrar soluciones económicas en tiempos de cálculo razonables.
En cuanto a los parámetros, estos son magnitudes que se toman como datos fijos y que, por tanto, no afectan a la optimización de la estructura. Los principales parámetros se han dividido en geométricos, de carga, de coste, de armado y de exposición. Entre ellos destacan el ancho del tablero, la inclinación de las almas, la luz de las vigas, la esbeltez mínima de la viga, las cargas muertas, la distancia de transporte, la dificultad de montaje de las vigas y el tipo de aceros. Las condiciones de durabilidad son las relativas a la instrucción de hormigón EHE [30]. En la tabla 1 se recogen los principales parámetros empleados en el problema.
Parámetros del problema
| Parámetros | Notación y valores |
| Geométricos | |
| Ancho del tablero | b_total=12,00m |
| Inclinación alma (° sexagesimales) | Ia=80° |
| Pendiente cartela ala superior (1: ns3) | ns3=3 |
| División base ala superior | s3=3 |
| Pendiente cartela ala inferior (1: ni3) | ni3=3 |
| División base ala inferior | i4=4 |
| Entrega de la viga | Ent=0,47m |
| Esbeltez mínima viga | Esb=(L/17) |
| De carga | |
| Ancho de las barreras | a_bar=2×0,5m |
| Espesor nominal del pavimento | e_pav=9cm |
| Carga muerta no procedente del pavimento | Qm=2×5,0kN/m |
| De coste | |
| Distancia transporte (ida) | d_transporte=50km |
| Despunte armadura activa | 25% |
| De armado | |
| Tipo de acero pasivo (B-500-S) | fyk=500N/mm2 |
| Tipo de acero activo (Y1860-S7) | fpk=1700N/mm2 |
| Diámetro torones acero activo | Φ=0,6” |
| Armadura piel viga | Φ=8mm |
| Fundas torones | Nivel 2 y 3 |
| Esbeltez vertical cercos | 200 (longitud/Φ) |
| De exposición | |
| Ambiente de exposición externo | IIb (EHE) |
El coste del tablero es consecuencia tanto del volumen de los materiales empleados como de la mano de obra, de la maquinaria y de los medios auxiliares necesarios para su ejecución. El proceso constructivo aquí empleado contempla la fabricación de las vigas artesa en la planta de prefabricados, su transporte a obra, su colocación y la posterior ejecución de la losa in situ. La evaluación del coste, por tanto, se obtiene multiplicando los precios unitarios de cada una de las unidades de obra (tabla 2) por las mediciones correspondientes y sumando los costes fijos relacionados con el proceso constructivo.
Precios básicos de la función coste de la estructura
| Descripción | Precio unitario en viga (€) | Precio unitario en losa (€) |
| Kilogramo del acero pasivo (B-500-S) | 2,63 | 1,40 |
| Kilogramo del acero activo (Y-1860-S7) | 3,38 | NA |
| m de molde en viga | 75,11 | NA |
| m2 de encofrado en losa | NA | 30,00 |
| m3 de hormigón HA-25 | NA | 64,99 |
| m3 de hormigón HA-30 | NA | 69,95 |
| m3 de hormigón HA-35 | NA | 74,03 |
| m3 de hormigón HA-40 | NA | 79,12 |
| m3 de hormigón HP-35 | 122,25 | NA |
| m3 de hormigón HP-40 | 133,40 | NA |
| m3 de hormigón HP-45 | 142,15 | NA |
| m3 de hormigón HP-50 | 152,89 | NA |
La determinación del coste unitario del acero de la armadura pasiva depende del diámetro empleado. Para su correcta valoración, este coste se encuentra afectado por 2 coeficientes correctores: uno en función del rendimiento de su manipulación por los operarios, y otro debido al distinto precio del suministro en rollos. Dichos coeficientes se representan en la tabla 3, tomando como base el precio del acero de 12mm de diámetro.
Coeficientes correctores kilogramo de acero (B-500-S)
| Diámetro | Viga suministro | Viga manipulación | Losa suministro | Losa manipulación |
| D6 | 1,250 | 1,400 | 1,250 | 1,400 |
| D8 | 1,170 | 1,250 | 1,170 | 1,250 |
| D10 | 1,075 | 1,100 | 1,075 | 1,100 |
| D12 | 1,000 | 1,000 | 1,000 | 1,000 |
| D16 | 0,980 | 0,900 | 0,980 | 0,900 |
| D20 | 0,980 | 0,900 | 0,980 | 0,900 |
| D25 | NA | NA | 1,000 | 0,800 |
| D32 | NA | NA | 1,000 | 0,800 |
Tras el acopio temporal de la viga fabricada, esta se transporta con un coste que dependerá del peso de la viga y de la distancia a la obra. Sin que el problema pierda generalidad, en el caso analizado se han considerado distintos costes en función del peso para una distancia de transporte de 50km (tabla 4). Por otro lado, el coste de la colocación de la viga depende de factores tales como la distancia (dietas de los operarios), la longitud de la viga (grúa a emplear) y la dificultad del montaje. En nuestro caso, la dificultad del montaje se ha considerado afectada con un parámetro de valor unitario, pudiéndose variar para otros casos particulares sin que la metodología empleada pierda generalidad.
Asimismo, es importante destacar que la función de coste puede verse penalizada para el caso en que las soluciones incumplan alguna de las restricciones impuestas. Ello posibilita al algoritmo la búsqueda a través de soluciones que, aunque no sean factibles, puedan estar cerca de óptimos locales de gran calidad. La expresión (4) indica el tipo de penalización empleada:
donde F+ representa la función penalizada, F es el coste, Φ es el porcentaje máximo de incumplimiento para un estado límite, gi es la generación a la que pertenece dicha solución dentro del algoritmo, y g es el número total de generaciones, teniendo un efecto creciente la penalización conforme avanza el proceso para facilitar la convergencia final hacia soluciones factibles.2.3Modelo estructural y comprobación del tablero
El empleo de las técnicas heurísticas implica la comprobación de los estados límite a los que se encuentra sometida una estructura completamente definida. Para conocer la factibilidad de la estructura frente a los estados límite, es necesario conocer las acciones que intervienen durante el proceso constructivo de la viga prefabricada y de su puesta en obra. La viga artesa se ve sometida a las acciones que provienen de la fuerza de pretensado y del peso propio en la planta de prefabricados, del peso propio de la losa ejecutada in situ sobre la viga y del peso propio y acciones variables sobre el conjunto estructural formado por las vigas y la losa. El conocido método de la rigidez proporciona los desplazamientos en los nodos de la estructura {D}, una vez conocidos el vector de cargas nodales {F} y la matriz de rigidez [K], mediante la ecuación matricial {F}=[K] * {D}.
El primer caso analizado corresponde al estado de cargas más desfavorable hallado cuando la viga se encuentra en la planta de prefabricados y cuando se coloca en obra y recibe el hormigón fresco de la losa (fig. 4). La viga se discretiza en 20 barras y 21 nodos, usándose el método del emparrillado plano que ha sido resuelto mediante un análisis matricial lineal elástico. El segundo caso analizado se basa en un modelo estructural representado por un emparrillado plano, formado por 2 vigas paralelas discretizadas como antes, pero con conexiones entre los nudos de las barras anteriores (fig. 5). Para la representación de las características mecánicas de las barras se ha considerado un emparrillado [30] en el que se ha considerado el efecto en las tensiones longitudinales de la distorsión de la sección transversal. En el cálculo se han contemplado las combinaciones de las acciones que producen máximos flectores, cortantes y torsores para los diferentes estados límite, tanto de servicio como últimos.
![Modelización del tablero, definido como emparrillado tipo 1 según Manterola [31].](https://static.elsevier.es/multimedia/02131315/0000003000000003/v1_201407030032/S0213131513000606/v1_201407030032/es/main.assets/gr5.jpeg?xkr=ue/ImdikoIMrsJoerZ+w997EogCnBdOOD93cPFbanNdXKeWo41MsZkw/vID64cuEiPuJOyFxneHZSkizm96flfSXRuq9PKMsoqi+EW0rIT4RmdB0rXJVLjYZ01/kXBRm/vfrf54dq7cK8OTsi08o2F+1tCOSkPtWxTTtu3TNqvizVEI66cyzPF3CqNNU2lichgc/bFi0BkIN7wlN3sxhbIYX4ftTvFEpXQcrqSJI661sTJosfS/U16QnaPlLb4a4TAJCrbrYFpxbk3pYncWByYEitrkmUSzvGfgAIXd2qVQKBk31Y1iYnTFMNfOseE9pJOsud4o+734q3N8fBZlUqA== "Modelización del tablero, definido como emparrillado tipo 1 según Manterola [31].")
Modelización del tablero, definido como emparrillado tipo 1 según Manterola [31].
Las acciones permanentes consideradas de valor constante son el peso propio y las cargas muertas (pavimento y pretiles). Los pretiles constituyen una carga distribuida uniforme de 5kN/m; las cargas de valor no constante sobre la viga consideradas son el pretensado, la retracción y la fluencia. Las acciones variables aplicadas sobre el tablero se corresponden al tren de cargas en puentes de la IAP-98 [32], si bien cabe señalar que recientemente en España se ha modificado dicho tren para adaptarse al contenido de los Eurocódigos. En cada estado individual de cargas se hace actuar una sobrecarga repartida uniforme de 4kN/m2 sobre el área limitada entre 2 nodos, el eje de una viga y el extremo del tablero o el eje del mismo (fig. 6). Además, se considera una carga puntual de 600kN, que en cada estado de cargas actúa sobre un nodo y a la máxima excentricidad exterior, teniendo en cuenta la separación a los pretiles (fig. 7). La consideración del carro como carga puntual constituye una simplificación del lado de la seguridad en luces medias y pequeñas; sin embargo, en las comprobaciones realizadas en las secciones transversales del tablero se ha contemplado dicha carga puntual correspondiente a un vehículo pesado como un conjunto de 6 cargas de 100kN separadas en sentido longitudinal 1,50 y 2,00m en sentido transversal. Se ha considerado que no existe una diferencia significativa (menos de un mes) entre las edades de los hormigones de las vigas prefabricadas y de la losa, pues en caso contrario habría que contemplar los efectos de las pérdidas y deformaciones diferidas antes y después del hormigonado de la losa [33].
Por otra parte, se han considerado las restricciones habituales que contempla la normativa española para este tipo de estructuras [30,32]. Eso incluye las comprobaciones de los estados límites últimos y de servicio para flexión, cortante y torsor, así como fisuración, fatiga y las deformaciones (instantáneas y diferidas) para la envolvente de tensiones debida al peso propio y a las cargas del tráfico. El ELS de fisuración se ha comprobado para cada una de las fases constructivas de la viga (tesado y tras el peso del hormigonado de la losa) de forma que la tensión en la fibra más traccionada no supere la resistencia media a tracción. Además, se ha comprobado, siguiendo lo dispuesto en la norma EHE [30], que la fibra más traccionada de la sección no supere la resistencia media a flexotracción del hormigón. Las tensiones y las reacciones se obtienen como resultado de un programa de desarrollo propio aplicado a los 2 modelos descritos. Las flechas instantáneas se limitan a 1/250 de la luz para las cargas frecuentes y las diferidas a 1/1.000 de la luz para la combinación de cargas casi permanente. Se considera la fatiga del hormigón y del acero. Además se comprueban las rasantes entre alas/alma/losa, y las flexiones transversales locales en distintos puntos de la losa y de las alas y alma de la viga.
El uso de 2 tipos distintos de hormigón, en viga y en losa, dificulta las comprobaciones, por lo que se ha simplificado el proceso mediante la homogenización del hormigón de la losa respecto al de la viga. Para ello, se toma como factor la relación entre el módulo de elasticidad de ambos en la obtención de las características mecánicas de la secciones y en las comprobaciones en ELS, y la relación entre las resistencias características en las comprobaciones en ELU. En el primer caso, este coeficiente se aplica sobre el ancho de la losa manteniéndose el canto. En el caso de la resistencia a torsión, el factor se aplica sobre el canto para el cálculo del espesor de la pared. La acción que ejerce el pretensado en la viga se ha desglosado en 4 niveles, tres en el ala inferior y uno en las alas superiores. Se ha tomado como simplificación en este trabajo la colocación de las fundas sobre todos los torones de un mismo nivel. Cuando en una sección existen uno o varios niveles de fundas, en ellos no actúan el pretensado, no teniéndose en cuenta ni en su acción ni en las características mecánicas de la sección. Las pérdidas de la acción del pretensado sobre toda la viga se calculan por nivel y se reparten entre las secciones donde la armadura activa es adherente.
3Algoritmo memético híbridoLos algoritmos meméticos constituyen una clase de metaheurísticas estocásticas que combinan la naturaleza de búsqueda en paralelo realizada por los algoritmos evolutivos y la búsqueda local que mejora cada una de las soluciones que forman una población. El calificativo de meméticos lo aplicó por primera vez Moscato [34] en 1989 a estos algoritmos híbridos inspirándose en el término inglés meme, acuñado por Dawkins [35] para referirse al análogo del gen en el contexto de la evolución cultural. La idea de emplear la hibridación entre metaheurísticas basadas tanto en algoritmos poblacionales como en otros de búsqueda local permite mejorar la efectividad de los híbridos al combinar los esfuerzos de diversificación e intensificación [5]. Krasnogor y Smith [36] proporcionan una guía y ejemplos de aplicación para este tipo de algoritmos.
En cuanto al algoritmo de búsqueda local empleado dentro del algoritmo memético, se presenta en este trabajo una variedad de la técnica de búsqueda de entornos a gran escala denominada Very Large-Scale Neighborhood Search (VLSN). En particular, y siguiendo la clasificación propuesta por Ahuja et al. [37], la variante empleada pertenece a los métodos de búsqueda por entornos variable, Variable-Depth Neighborhood Search (VDNS). Aunque una de las primeras aplicaciones de esta estrategia puede encontrarse en la resolución de problemas de rutas [38], por lo que se ha podido averiguar, esta es la primera vez que se aplica este tipo de búsqueda local en la optimización de estructuras. VDNS se basa en realizar una búsqueda local con un movimiento hasta alcanzar un óptimo local; a continuación se cambia el movimiento a otro más amplio que permita escaparse de dicho óptimo y se sigue así hasta completar un número predefinido de movimientos, cada uno de ellos más amplio que el inmediatamente anterior. En el algoritmo propuesto en el presente artículo, el primer movimiento se define por el cambio aleatorio de una sola variable, eligiendo siempre la nueva solución si mejora a la anterior. El segundo movimiento mueve 2 variables simultáneamente, y así sucesivamente. Debe definirse, en este caso, un número de movimientos sin mejora para pasar de un movimiento al siguiente.
Por tanto, el proceso del algoritmo híbrido completo se inicia con la generación aleatoria de una población, en nuestro caso de 500 individuos (n). Cada una de estas soluciones se mejora mediante una búsqueda VDNS, hasta alcanzar un óptimo local. Para ello, comienza moviendo solo una variable, y cuando lleva 10 movimientos continuos sin mejora (no_mej), se incrementa en uno el número variables (var) simultáneas que se mueven, hasta un máximo de 5. Con esta nueva población mejorada de 500 individuos se procede a la aplicación de un algoritmo genético. Los algoritmos genéticos hacen evolucionar una población de individuos que se someten a acciones aleatorias (mutaciones y cruzamientos) de forma que se seleccionan los individuos que se consideran mejor adaptados [6]. Este algoritmo requiere la calibración de sus parámetros. Para el algoritmo memético híbrido propuesto, se han utilizado los mejores parámetros obtenidos por Martí [29] en su tesis doctoral: 500 individuos en cada población, 200 generaciones, cruzamiento igual a 0,50 y elitismo. Para la evaluación de la población se valora el coste penalizado; sin embargo, en la heurística VDNS el coste considerado está sin penalizar —es decir, solo se aceptan las soluciones que cumplen con todas las restricciones impuestas—, para evitar la divergencia prematura del algoritmo. A la nueva generación obtenida se le vuelve a aplicar el VDNS, repitiéndose el proceso hasta un máximo de 200 generaciones. En la figura 8 se muestra un diagrama básico de flujo del algoritmo memético aplicado en este estudio.
4Estudio paramétrico y discusión de resultados
Los algoritmos y el módulo de comprobación estructural se programaron en Intel® Visual Fortran Compiler Integration for Microsoft Visual Studio 2008. Cada proceso duró una media de 1.460min sobre un ordenador con un procesador INTEL® CoreTM i7 CPU X980 de 3,33GHz. Se procesaron 12 ejecuciones para cada una de las 5 luces de vano consideradas de 20, 25, 30, 35 y 40m, obteniéndose valores mínimos, medios y desviaciones de los resultados obtenidos para el coste, las variables geométricas y la armadura activa (tabla 5). Se hace notar que el valor medio de las variables no toma necesariamente los valores normalizados permitidos para cada una.
Algoritmo memético híbrido. Resultados geométricos y armadura activa
| Luz (m) | Coste (euros) | h1 (m) | e4 (m) | b1 (m) | b3 (m) | e1 (m) | e2 (m) | e3 (m) | fc,viga (MPa) | fc,losa (MPa) | p1 (n) | p2 (n) | p3 (n) | p4 (n) | Sv (m) | |
| 20 | (a) | 54.865 | 1,13 | 0,18 | 1,53 | 0,26 | 0,16 | 0,12 | 0,17 | 35 | 25 | 27 | 9 | 0 | 4 | 5,58 |
| (b) | 57.067 | 1,13 | 0,19 | 1,47 | 0,24 | 0,16 | 0,10 | 0,17 | 35,83 | 28,75 | 26 | 10,75 | 0,00 | 3,00 | 5,62 | |
| 25 | (a) | 66.591 | 1,34 | 0,19 | 1,10 | 0,23 | 0,15 | 0,10 | 0,10 | 35 | 25 | 19 | 19 | 0 | 2 | 5,62 |
| (b) | 69.480 | 1,39 | 0,19 | 1,56 | 0,25 | 0,16 | 0,10 | 0,16 | 37,08 | 27,92 | 28 | 14,67 | 0,00 | 2,80 | 5,57 | |
| 30 | (a) | 78.826 | 1,60 | 0,18 | 1,31 | 0,23 | 0,20 | 0,10 | 0,15 | 35 | 35 | 23 | 23 | 5 | 2 | 5,27 |
| (b) | 82.167 | 1,64 | 0,18 | 1,34 | 0,25 | 0,18 | 0,10 | 0,15 | 39,17 | 30,42 | 23 | 23,33 | 0,42 | 2,20 | 5,52 | |
| 35 | (a) | 99.840 | 1,82 | 0,18 | 1,63 | 0,29 | 0,19 | 0,10 | 0,16 | 40 | 30 | 29 | 29 | 0 | 2 | 5,77 |
| (b) | 101.604 | 1,93 | 0,18 | 1,61 | 0,27 | 0,19 | 0,10 | 0,15 | 38,33 | 30,00 | 29 | 26,42 | 1,17 | 2,20 | 5,56 | |
| 40 | (a) | 109.855 | 2,20 | 0,18 | 1,44 | 0,23 | 0,21 | 0,10 | 0,15 | 35 | 30 | 25 | 25 | 15 | 2 | 5,75 |
| (b) | 114.028 | 2,18 | 0,18 | 1,44 | 0,25 | 0,19 | 0,10 | 0,16 | 38,33 | 33,75 | 25 | 25,17 | 6,92 | 2,00 | 5,60 |
(a) Solución coste mínimo.
(b) Valores medios.
La figura 9 muestra una evolución típica del coste medio de la población de soluciones y el coste de la solución de menor coste respecto al tiempo de cálculo. Al principio, existe cierta dificultad en la obtención de las primeras 500 soluciones, alrededor de 11.500segundos. Ello se debe a la complejidad que encuentra el algoritmo en alcanzar una combinación de variables en la que encaje el pretensado, para elaborar una solución factible —alrededor de 23segundos—. En la figura 10 se representa la evolución de los costes medios de la población y de la mejor solución tanto para la mejor ejecución como para los valores promedio de las 12 ejecuciones realizadas del algoritmo. En todos los casos se observa que el comportamiento es similar, de forma que los costes se reducen rápidamente al principio y luego se hace de una forma más lenta hasta estabilizarse. Se aprecia para la evolución representada en la figura 10 que, en ambos casos, a partir de la generación 75 el mejor coste de ambas heurísticas mejora poco, 2,5 y 2,0%, respectivamente, y es prácticamente inapreciable a partir de la generación 100 hasta la 200. En la figura 11 se ha representado el coste en función del número de generaciones para las 12 ejecuciones realizadas.
Esta heurística híbrida permite 2 tipos de mejoras: la genética para la población de 500 individuos y la local para cada uno de ellos en las 200 generaciones. En la figura 12 se observa que aproximadamente el 90% de las mejoras que sustituyen a la mejor solución encontrada por el algoritmo hasta ese momento se produce en las primeras 120 generaciones. En la figura 13 se puede observar la efectividad del algoritmo híbrido, donde en las primeras generaciones ocurre la mejora de la mejor solución encontrada debido a la búsqueda local (VDNS), siendo necesaria la conjunción de la búsqueda local dentro de la búsqueda genética para alcanzar óptimos locales de mayor calidad. Estos resultados permiten comprobar la pertinencia de utilizar 200 generaciones del algoritmo memético con la búsqueda local VDNS definida.
La figura 14 muestra un crecimiento lineal del coste total de la estructura con la luz, de forma que el 98,23% de la variabilidad observada en el coste puede explicarse con la luz. Este incremento se debe al mayor coste del material necesario para resistir los esfuerzos crecientes y satisfacer las limitaciones impuestas a la deformación. De hecho, el coste medio se duplica cuando la luz pasa de 20 a 40m. También podemos ver en la tabla 5 que el promedio calculado para las 5 luces consideradas de la desviación del coste medio de cada una de las 12 ejecuciones respecto al valor mínimo ha sido solo del 3,50%.
En la figura 15 se observa un ajuste lineal del canto de la viga h1 respecto a luz, con un coeficiente de regresión muy elevado (R2=0,9864). De hecho, la media obtenida para las 12 ejecuciones en cada una de las 5 luces analizadas ha sido de L/18,08, valor que no supera, pero es cercano, a la limitación de esbeltez impuesta de L/17 (tabla 1). Se ha comprobado que, de todas las restricciones, la condición determinante del diseño ha sido la capacidad resistente del tablero en el ELU a flexión. En cambio, el espesor de la losa e4 no presenta una correlación apreciable con la luz, siendo el valor medio de las 60 ejecuciones realizadas en las luces analizadas de apenas 0,184m, un valor realmente bajo que tiende a reducir el peso de la estructura pero es necesario para resistir la flexión transversal de las sobrecargas variables, con una desviación típica de solo 0,008m. También resulta de interés comprobar la alta relación lineal (R2=0,7196) que existe entre el número de torones necesarios y la luz (fig. 16). El resto de valores geométricos no presentan correlaciones significativas con la luz, y sus valores medios y mínimos se pueden ver en la tabla 5. Se puede observar que los espesores de las alas y del alma de la viga son los mínimos para permitir reducir su peso pero, a la vez, poder albergar las armaduras activas en las alas y resistir a cortante y torsor en las almas. Además, la longitud de las alas superiores de las vigas b3 es reducida, con un valor medio de 0,251m, espacio suficiente para ubicar la armadura activa superior. El ala superior colabora con la losa para incrementar el canto del conjunto y así resistir mejor la solicitación a flexión transversal que produce la acción de la sobrecarga variable; sin embargo, su dimensión no se ve muy condicionada. Por otra parte, la separación media entre vigas Sv es de 5,575m, siendo esta separación inferior a la mitad de la dimensión del tablero (6m), para poder ajustarse mejor a la distribución de la sobrecarga variable, en la que el carro se encuentra desplazado del extremo.
En cuanto a las calidades del hormigón, el valor medio fck de resistencia característica de las 60 ejecuciones es de 37,75MPa para la viga y de 30,17MPa para la losa. Se comprueba que, para los costes unitarios empleados, el algoritmo no ha elegido la máxima resistencia característica disponible.
En la tabla 6 se han recogido las cuantías necesarias de armaduras pasivas y de hormigón, tanto en viga como en losa. Se observa, en las figuras 17 y 18, un crecimiento con la luz de las necesidades del acero pasivo en la viga y en la losa, con una buena correlación para el caso de la viga (R2=0,8101), pero no tanto para el caso de la losa (R2=0,3896). El valor medio del acero pasivo resultante para la viga y para la losa es de 5.375 y 9.580kg, respectivamente, lo que representa una cuantía media de 40,39kg/m2 de tablero. En la tabla 7 se ha recogido la armadura transversal en la sección de apoyo, sin inclusión de riostras, para la solución óptima obtenida por el algoritmo en cada una de las luces consideradas. En cuanto al hormigón necesario, en la figura 19 se observa un crecimiento del volumen por unidad de superficie del tablero en la viga con la luz, con un ajuste que explica el 79,76% de la variabilidad. El hormigón necesario en la losa no guarda correlación con la luz, siendo su valor medio de 0,183 m3/m2. Para todas las luces y casos analizados, es necesaria una cantidad de hormigón total de 0,288 m3/m2.
Valores medios y de las soluciones de coste mínimo de las cuantías necesarias de acero pasivo y hormigón, tanto en viga como en losa obtenidas por el algoritmo memético híbrido
| Luz (m) | Acero pasivo viga (kg) | Acero pasivo losa (kg) | Total acero pasivo (kg/m2) | Hormigón viga (m3/m2) | Hormigón losa (m3/m2) | |
| 20 | (a) | 3.725 | 7.832 | 37.042 | 0,094 | 0,180 |
| (b) | 2.917 | 7.850 | 42.725 | 0,082 | 0,186 | |
| 25 | (a) | 3.523 | 6.961 | 41.606 | 0,076 | 0,190 |
| (b) | 3.920 | 8.567 | 40.024 | 0,095 | 0,188 | |
| 30 | (a) | 4.637 | 6.900 | 31.013 | 0,100 | 0,180 |
| (b) | 4.955 | 9.627 | 39.197 | 0,100 | 0,181 | |
| 35 | (a) | 8.207 | 6.965 | 35.124 | 0,121 | 0,180 |
| (b) | 7.179 | 10.291 | 40.440 | 0,123 | 0,183 | |
| 40 | (a) | 5.634 | 10.266 | 32.317 | 0,128 | 0,180 |
| (b) | 7.904 | 11.566 | 39.575 | 0,125 | 0,180 |
(a) Solución coste mínimo.
(b) Valores medios.
Armadura transversal en la sección de apoyo sin inclusión de riostras para la solución óptima obtenida por el algoritmo memético híbrido
| Luz (m) | As6 | As7 | t1 | t2 | t3 | t4 | t5 |
| 20 | 49Ø16 | 49Ø10 | Ø10/200 | Ø12/200 | Ø8/200 | Ø20/300 | Ø20/300 |
| 25 | 61Ø10 | 49Ø8 | Ø12/200 | Ø12/200 | Ø12/200 | Ø20/200 | Ø20/300 |
| 30 | 61Ø12 | 49Ø6 | Ø16/200 | Ø12/200 | Ø8/200 | Ø16/200 | Ø16/300 |
| 35 | 73Ø12 | 49Ø6 | Ø16/200 | Ø12/200 | Ø10/200 | Ø16/250 | Ø16/300 |
| 40 | 73Ø10 | 49Ø6 | Ø10/200 | Ø10/200 | Ø6/200 | Ø20/250 | Ø20/300 |
Resulta interesante comparar la mejor solución alcanzada por el algoritmo respecto a una estructura realmente construida y calculada mediante procedimientos habituales. Se han comprobado para casos similares ahorros apreciables en torno al 7-8%. Sin embargo, en algún caso extremo, como el caso del viaducto1 del tramo Muro de Alcoy-Puerto de Albaida del proyecto de construcción de la autovía del Mediterráneo (fig. 20), el ahorro se ha estimado en un 50%, tomando los precios de la tabla 2. En este caso el puente tenía una luz de 35m y un ancho de tablero igual al de la solución optimizada, siendo el ahorro alcanzado tan importante a causa de las diferencias en la medición de las unidades de obra en materiales que pueden apreciarse en la tabla 8.
Comparación de las mediciones en las unidades de obra significativas correspondientes al viaducto 1 del tramo Muro de Alcoy-Puerto de Albaida, de luz 35m, respecto a la solución optimizada
| Unidad | Viaducto 1 | Solución optimizada |
| m3 hormigón en viga | 75,62 | 52,31 |
| m3 hormigón en losa | 107,82 | 77,63 |
| kg acero pasivo viga | 11.218,00 | 8.207,00 |
| kg acero pasivo losa | 50.484,00 | 6.966,00 |
| kg acero activo | 7.716,00 | 5.936,00 |
El trabajo presenta un algoritmo memético que realiza una hibridación entre un algoritmo genético y un algoritmo VDNS de búsqueda local aplicado con éxito al diseño automático y a la optimización del coste de tableros de puentes isostáticos de vigas artesa pretensadas prefabricadas. Este algoritmo, ejecutado varias veces, proporciona escasas variaciones respecto a la media de los valores encontrados, alcanzando buenas soluciones, incluso en el caso de una sola ejecución. Del estudio paramétrico realizado para luces comprendidas entre 20 y 40m, se ha comprobado un fuerte ajuste lineal del canto con la luz, con valor medio de L/18,08. El espesor medio de la losa resultante es de apenas 0,184m, de forma que las soluciones optimizadas tienden a reducir al máximo el peso de la estructura, pues los espesores de alas y el alma de la viga también son muy reducidos. La resistencia característica del hormigón en la viga es superior al de la losa, con valores medios de 37,75 y 30,17MPa, respectivamente, no habiendo sido necesario llegar hasta la máxima resistencia posible de 50MPa. La cuantía de armadura pasiva y el volumen por unidad de superficie necesaria en la viga crecen con la luz. Son necesarios, como valores medios, 40,39kg de armadura pasiva y 0,288 m3 de hormigón por cada metro cuadrado de tablero de puente. Por último, los resultados muestran que el algoritmo es capaz de reducir el coste de modo significativo, en torno al 8%, para este tipo de estructuras que se encuentran altamente industrializadas. Sin embargo, se han encontrado ahorros que pueden llegar al 50% en algunos casos realmente ejecutados. Como futuras líneas de investigación se apunta la consideración de otros hormigones, como los de alta resistencia o los reforzados con fibras; además, resulta de interés la optimización multiobjetivo considerando otros aspectos, como la reducción de emisiones de gases de efecto invernadero.
Los autores agradecen el aporte financiero realizado para este trabajo por parte del Ministerio de Ciencia e Innovación (Proyecto de Investigación BIA2011-23602). Los autores también agradecen los comentarios y sugerencias constructivas realizadas por los revisores anónimos.