En la construcción de edificaciones se han incorporado nuevos materiales y técnicas constructivas, lo que ha conducido al desarrollo de diversos tipos sistemas estructurales cada día más complejos y difíciles de analizar. Esto hace que las herramientas de análisis deban evolucionar en la misma medida. Los sistemas estructurales se clasifican como sistemas tradicionales y no tradicionales, así mismo las herramientas que se emplean para su cálculo. Entre las herramientas numéricas de análisis, podemos citar el Método de Elementos Finitos (MEF). Las estructuras tradicionales están conformadas con materiales homogéneos y dimensiones regulares, por esta razón el método de análisis empleado en preferencia es el MEF (1, 2, 3). Por el contrario, las estructuras No Tradicionales están conformadas por elementos heterogéneos que combinan diferentes tipos de materiales y con secciones transversales compuestas por diversos materiales, lo que dificulta definir las propiedades mecánicas de la sección transversal de los elementos que lo conforman. Usar los MEF para analizar numéricamente estructuras no tradicionales, requiere de ajustes de diversos parámetros de los modelos usados en los elementos finitos para lograr reproducir los resultados experimentales. Este proceso es complicado, y aún más complejo al tomar en cuenta la no linealidad en el comportamiento de los materiales producto del daño (perdida de rigidez) y de las deformaciones permanentes (elongaciones permanentes o pandeo local) (4, 5). Para estos casos se deben emplear métodos que puedan facilitar la simulación y reproducir estos fenómenos complejos. Otra de las desventajas del MEF, para estructuras compuestas, es que los métodos de análisis actuales requieren de muchas iteraciones para la inclusión de la no linealidad, por lo que el proceso es demándate en tiempo de cálculo, limitando su aplicabilidad.

Los sistemas tipo Sándwich son sistemas constructivos no tradicionales. Este sistema Sándwich es rápido y fácil de construir, capaz de soportar cargas externas y tener un buen comportamiento mecánico (6, 7). Los muros de corte de un sistema tipo sándwich pueden comportarse de buena forma bajo cargas laterales (8), así como también a fuerzas de impacto (9). La sección transversal de los elementos tipo sándwich está conformada por variedad de materiales, lo que dificulta realizar su modelación numérica mediante el MEF. En los elementos tipo sándwich es difícil definir las propiedades mecánicas entre los planos de contacto de los materiales que la conforman.

El avance tecnológico computacional ha permitido desarrollar algunos modelos numéricos para el análisis de estas estructuras sándwich (10). En el comportamiento de las estructuras sándwich se reflejan varios mecanismos de fallas, los cuales son analizados y explicados experimental y numéricamente en varios artículos de investigación (11). Entre algunos mecanismos de fallas se encuentran la inestabilidad local (12), el esfuerzo de corte en el núcleo, la pérdida global de la estabilidad (13). El análisis de estos mecanismos de falla requiere el uso de modelos numéricos apropiados que son difíciles de implementar, producto de la incertidumbre en la definición de los parámetros requeridos en los modelos numéricos. Producto del planteamiento anterior, se requiere desarrollar nuevos métodos que representen los mecanismos de fallas antes citados.

En este artículo se presenta el desarrollo de un modelo basado en redes neuronales artificiales (RNA), el cual consiste en un programa de recolección de datos experimentales para conocer el comportamiento inelástico de una losa de techo tipo sándwich sujeta a flexión pura. El mismo pretende emular la evolución de las deformaciones permanentes y la caída de la rigidez inicial según el daño producido en los materiales. Otro aspecto que se aborda es la disipación de energía producida por el amortiguamiento de los materiales que conforman el elemento (14). Este efecto de amortiguamiento se evidencia al observar la curva de Carga vs. Desplazamiento, donde se aprecia que el recorrido del ciclo de descarga y recarga no es el mismo, como se observará en la sección de ensayo experimental.

El uso de RNAs facilitaría el análisis, ya que emplea un conjunto de ecuaciones más simples y eficientes, para resolver el problema (15, 16), mientras que los métodos convencionales requieren del uso de un conjunto de ecuaciones más complejas. Los métodos de RNA ha sido aplicado no solo en análisis estructural, sino también en el área de geotecnia, estructuras, hidráulica y nanotecnología (17, 18, 19, 20), pero poco empleado para modelar el comportamiento de estructuras de tipo compuesto, como los paneles tipo sándwich, que se analizan en el presente trabajo.

Las variables de daños en elementos de hormigón armado miden el nivel de fisuración o agrietamiento que sufre el hormigón (4, 5), mientras que en las estructuras de acero, las variables de daño miden el grado de pandeo local del elemento de acero (1). Ambos fenómenos físicos (grietas o fisuras en el hormigón y pandeo local del perfil de acero) provocan la pérdida de rigidez del elemento y por ende de la estructura. También se propone la evaluación del daño global del panel tipo sándwich experimentalmente. Por ser un panel compuesto con diferentes materiales, el daño global que se presentan al ser sometido a acciones externas son la combinación de fisuras, grietas y pandeo del panel de techo tipo sándwich. Al final del artículo, se presenta una evaluación y comparación de esta variable de daño global obtenidas mediante simulaciones usando las RNA y los observados experimentalmente.

Se construyeron y ensayaron un grupo de paneles compuestos en el Laboratorio de Mecánica Estructural de la Universidad Centroccidental “Lisandro Alvarado”, UCLA - Venezuela. Cada panel es un elemento básico estructural no tradicional que conforma la losa de techo de una estructura más compleja. El tipo de ensayo experimental realizado a cada panel construido es el mostrado en la figura 3 y el análisis solo se realizó a flexión, ya que las variables principales a evaluar son la fuerza aplicada y la flecha que alcanza el panel hasta su falla, al someter cada panel a esfuerzos a flexión pura.

Una red neuronal artificial (RNA) es una técnica basada en el comportamiento biológico del cerebro y desarrollada para ser utilizada en aplicaciones de aprendizaje automático. La neurona es la unidad básica del cerebro y en el aprendizaje de máquinas se imita utilizando el llamado perceptrón (23). El perceptrón introducido inicialmente por Rosenblatt, se ha ido mejorando a lo largo de los años y ha evolucionado en una compleja red de perceptrones para convertirse en lo que hoy denominamos las RNA. Las RNA están compuestas principalmente por tres capas: una capa de entrada vectorial, una o más capas ocultas y una capa de salida vectorial, como se muestra en la Figura 11.

Las RNAs más conocidas y utilizadas son las de alimentación hacia delante (Feedforward), una red de este tipo suele entrenarse utilizando un algoritmo de propagación hacia atrás. La capa de entrada (Xi) está conectada a la capa oculta por medio de pesos (wij). El resultado del producto escalar, entre entradas y pesos, se pasa a una función de transferencia para activar la capa de salida (Yi). Estos resultados de salida se comparan con un conjunto de objetivos durante la fase de entrenamiento, y los errores se calculan y corrigen cambiando el valor de los pesos, esto se realiza mediante un algoritmo de optimización que minimiza los errores en forma iterativa, este proceso se conoce con el nombre de “entrenamiento” de la red.

En el presente trabajo los datos experimentales se pre-procesaron antes de la construcción del modelo RNA. Esto se efectuó con la finalidad de eliminar el ruido introducido en los datos adquiridos de los equipos electrónicos. Los valores de fuerza y desplazamiento se “suavizaron” utilizando una media móvil, con un marco de 100 valores de datos, tanto para el desplazamiento como para la fuerza. Los datos posteriores al fallo del elemento estructural se eliminaron del conjunto de entrenamiento para no afectar la generalización de la misma, con estos valores pre-procesados se entrena la RNA. En total, se utilizaron 2555 pares de valores de fuerza y deformación en el proceso de entrenamiento.

Para validar los modelos, se realizaron tres pruebas para cada modelo seleccionado con y sin amortiguación. La primera prueba tiene que ver con el proceso de entrenamiento de la RNA. El entrenamiento de los modelos usando RNAs se realizó aplicando tres conjuntos de datos como se explicó en párrafos precedentes: un conjunto para el entrenamiento, otro conjunto para validación cruzada y el último conjunto para la verificación. Luego, los resultados del modelo se confrontaron con todo el conjunto de datos observados. La Figura 16 muestra el historial de desplazamiento usado para alimentar el modelo, se puede observar que los modelos con y sin amortiguación reproducen con precisión los datos medidos en el laboratorio, ver Figuras 17 y 18. El modelo con amortiguación puede reproducir con mayor precisión los datos observados al compararlo con el modelo sin amortiguación.

La segunda prueba tiene como objetivo verificar la capacidad de los modelos de RNAs para la interpolación. Esta es una prueba complementaria para verificar si el modelo obtuvo una buena generalización. El modelo fue sometido a una nueva historia de desplazamientos; la historia de desplazamiento se seleccionó de tal manera que no coincida con los datos previamente utilizados en el entrenamiento, prueba o verificación del modelo de RNA, ver Figura 16. La última prueba realizada fue la representación de la curva de comportamiento extrapolando los datos de entrenamiento del modelo. Similar a la interpolación, la extrapolación es una forma de verificar la generalización en el entrenamiento de la red neuronal artificial.

Las figuras 19 y 20, muestran los resultados gráficos de la interpolación y extrapolación para los modelos con y sin amortiguación respectivamente. Se puede observar que ambos modelos pueden representar el comportamiento inelástico con la interpolación y la extrapolación. El modelo con amortiguación pudo extrapolar un 20% más allá del rango utilizado para el entrenamiento con buena precisión, mientras que el modelo sin amortiguación pudo extrapolar el mismo 20% pero con menos precisión. En el modelo sin amortiguación, la curva se vuelve casi horizontal cuando extrapola, mientras que el modelo con amortiguación tiene un comportamiento más estable y realista al comprar con los datos experimentales.

Los resultados muestran cómo la red neuronal artificial puede modelar el comportamiento no lineal de la curva de Fuerza vs. Desplazamiento. Una de las características importantes a destacar, es cuando el modelo está sujeto a un ciclo de carga y descarga, se puede evidenciar que al descargar no regresa a su deformación inicial, sino que tiene una deformación permanente o residual. Todas las simulaciones realizadas con los dos modelos de RNA representaron de muy buena forma la deformación permanente y la caída de la rigidez al ser comparados con la curva experimental. Demostrando que ambos modelos, con y sin amortiguación, pueden replicar las curvas de comportamiento inelástico de los paneles compuestos. Al medir la evolución de las deformaciones permanentes y las fisuras o grietas de los materiales que conforman el panel, se puede deducir que el modelo tiene “memoria” de las cargas y deformaciones anteriores a las que ha sido sometido el elemento y considera el daño global producido.

En el momento del entrenamiento de las RNAs, con y sin tomar el efecto de amortiguación del material, se puede apreciar que la evolución de las deformaciones permanentes es idéntica a la evolución durante el ensayo experimental del panel (figura 21a). También se comparó la evolución de la variable de daño global entre los resultados experimental y los resultados durante el entrenamiento de las RNAs (figura 21b), y se puede apreciar que el modelo de RNA con efecto de amortiguación se ajusta mejor al comportamiento real del sistema no tradicional evaluado en este trabajo.

Al realizar las simulaciones, con y sin efecto de amortiguación, con una historia de desplazamiento arbitraria en el rango de los datos de entrenamiento (interpolación) (figura 16 color azul) y otros datos más allá del entrenamiento (extrapolación) (figura 16 color rojo) se observa que los modelos representan el comportamiento inelástico de forma adecuada. Al analizar los resultados de las simulaciones y ver la evolución de las deformaciones permanentes se puede apreciar que se ajusta muy bien a los resultados experimentales (figura 22a). Otra variable a analizar de los resultados simulados es la degradación de la rigidez y la obtención de la evolución del daño global del panel compuesto. Observando esta variable se aprecia que el modelo con amortiguación, considerando los datos de interpolación y extrapolación, se ajusta mucho mejor a los resultados experimentales (figura 22b).

El excelente rendimiento de ambos modelos, demuestra la viabilidad de utilizar modelos RNA para el análisis y evaluación estructural conformadas por elementos compuestos. Aunque se requiere de más investigación sobre la incorporación de otras variables, como la resistencia del mortero, la malla de refuerzo y el perfil de acero, ya que estas ayudarían a mejorar el comportamiento predictivo de los modelos de RNAs. Estas variables también ayudarían a usar el modelo no solo para la simulación sino también para el diseño y evaluación de estructuras tradicionales y no tradicionales. Debido al rendimiento obtenido, se puede argumentar que es factible construir un modelo híbrido para el análisis de estructuras más complejas. Los cuales podrían combinar esquemas numéricos y redes neuronales artificiales trabajando conjuntamente. Esto podría mejorar la modelación en términos de no linealidad, así como también ahorrar tiempo de cálculo.

Dos modelos de ANN han sido exitosamente instanciados y validados. Los mismos permiten calcular la fuerza aplicada dado el desplazamiento de un elemento estructural complejo y viceversa. Los modelos instanciados son capaces de simular los efectos no lineales producidos en los materiales debido a las deformaciones plásticas en el acero y las grietas en el hormigón reforzado. Además, los modelos pueden determinar la deformación residual al descargar el elemento, lo que permite estimar el daño global producido en el elemento.

Para el elemento compuesto estudiado, el mejor rendimiento promedio (EMC) se obtuvo para las RNAs con 2 nodos ocultos y 200-250 epochs en ambos modelos. Estos valores proporcionan la mejor generalización en la fase de entrenamiento de los modelos. El modelo que considera la amortiguación interna de los materiales tiene un rendimiento ligeramente superior que el modelo sin amortiguación (EMC de 4548.85), este reproduce las ramas en bucle, lo que significa que el modelo tiene en cuenta la disipación de energía en los materiales así como el daño que se produce.

Otra de las características resaltante en los modelos entrenados, es que los mismos son capaces de interpolar y extrapolar, demostrando que las RNAs poseen una buena generalización y es posible elaborar modelos complejos. Las RNAs pueden representar el comportamiento inelástico de paneles compuestos. Representando la evolución de las deformaciones permanentes y el daño global del panel de manera adecuada. Se abre la posibilidad de emplear las RNA en combinación con modelos de tipo numérico o método de elementos finitos para la incorporación de los efectos inelásticos, creando así un modelo de tipo híbrido.