Design Coefficients and Cracking Paths of Circular, Elliptic and Triangular Isolated Slabs

Juárez-Luna Gelacio, Caballero-Garatachea Omar

Departamento de Materiales, Área de Estructuras. Universidad Autónoma Metropolitana. Correo: gjl@correo.azc.uam.mx.

Asesores en Ingeniería Estructural. Correo: caballerog12@gmail.com.

Recibido: octubre de 2012, ]]> Reevaluado: diciembre de 2012,
Aceptado: febrero de 2013

Resumen

Se determinan, mediante la mecánica del daño continuo, coeficientes de diseño y trayectorias de agrietamiento en losas circulares, elípticas y triangulares. El comportamiento no lineal del concreto reforzado se simuló con un modelo de daño, el cual considera ablandamiento, mientras que el comportamiento del acero de refuerzo con un modelo de plasticidad bilineal con superficie de fluencia de Von Misses. Los modelos constitutivos y el método de los elementos finitos se validan comparando los resultados numéricos con los experimentales de una losa rectangular reportada en la literatura. La distribución del acero de refuerzo se propone con base en análisis elásticos lineales con elementos placa. Los coeficientes se calcularon con análisis elásticos y no lineales, observándose que en estos últimos existe una variación debido a la redistribución de esfuerzos. Se muestran las trayectorias del agrietamiento, las cuales son dependientes de las condiciones de apoyo. Los coeficientes y la distribución del acero de refuerzo para losas con estas geometrías son una base en la reglamentación para el análisis y diseño de losas con este tipo de geometrías.

Descriptores: losa, modelo de daño, coeficientes de diseño, elementos finitos, agrietamiento.

Abstract

Design coefficients and cracking paths for circular, elliptical and triangular labsare obtained by damage mechanics were determined. The nonlinear behavior of the reinforced concrete was simulated with a damage model which considers softening, whereas the behavior of the reinforcing steel was modeling with a bilinear plasticity model which has a Von Mises yield surface. The constitutive models and the finite element method were validated by comparing the computed numerical results with the experimental results of a rectangular slab reported in the scientific literature. The distribution of the reinforcing steel was proposed based on linear elastic analysis with plate elements. The coefficients were computed with elastic and nonlinear analyses, noticing that the last has a variation because of the stress redistribution. The crack paths are shown, which are dependent on the boundary conditions. The coefficients and the distribution of the reinforcing steel are a base on the regulation for the analysis and design of slabs with this kind of geometry.

Keywords: slab, damage model, design coefficients, finite elements, cracking.

Introducción

Los sistemas de piso son subsistemas estructurales de los más antiguos y utilizados por la humanidad, pues se les han dado diversos usos en la mayoría de las construcciones. La gran aceptación que tienen estos elementos como parte de un sistema estructural, se debe a las ventajas que presentan en cuanto a la transferencia de cargas verticales y horizontales.

Los sistemas de piso en estructuras de concreto y mampostería generalmente contemplan las geometrías rectangulares; sin embargo, algunos diseños arquitectónicos requieren de geometrías irregulares en planta. Para el diseño de losas de geometría rectangular existen métodos de análisis elásticos, o bien, se utilizan coeficientes tabulados en normas de diseño, mientras que las losas no rectangulares, como las mostradas en la figura 1, se caracterizan por no presentar formas comunes de estructuración y por carecer de metodologías de análisis.

Para el análisis de losas rectangulares, Westergaard y Slater (1921) propusieron una metodología con base en la teoría de la elasticidad, particularmente de losas apoyadas en dos direcciones y losas planas apoyadas en capiteles circulares mediante una redistribución de esfuerzos dependientes de la redistribución de rigideces, estos autores consideraron dos variables elementales: las condiciones de apoyo y el tipo de carga; algunos de estos resultados se compararon con pruebas experimentales reportadas en la literatura. Tiempo después, Marcus (1929) desarrolló un procedimiento aproximado, mediante el cual los momentos máximos elásticos en una placa apoyada sobre vigas rígidas se pueden determinar tomando como base el análisis de una franja de losa y obteniendo factores de reducción para corregir los momentos obtenidos. Además, realizó cuatro tablas de coeficientes de momentos, los cuales proporcionan valores para la distribución del panel de carga en dos direcciones, la evaluación de momentos negativos en bordes continuos y positivos para la mitad del centro del claro.

Fue hasta 1940 que el Comité del American Concrete Institute, ACI, realizó modificaciones a las propuestas realizadas por Westergaard y Slater (1921), extendiendo el método de análisis para cubrir todas las posibles combinaciones de paneles en un sistema de piso. Esto llevó a una tabla de coeficientes para momento positivo en el centro del claro y momento negativo en los bordes, ya sean continuos o discontinuos, considerando las condiciones de forma y de apoyo adyacentes al panel. Posteriormente, Newmark y Siess (1950) propusieron un método de diseño para losas en dos direcciones, resultado de estudios analíticos mediante un proceso de distribución de momentos; estos autores estudiaron diferentes condiciones de placas continuas con el fin de verificar los efectos relativos de las variables que influyen en el desarrollo de momentos en losas soportadas por vigas rígidas. El desarrollo en detalle del procedimiento aproximado de distribución de momentos fue reportado por estos mismos autores en Siess y Newmark (1950).

Para determinar los coeficientes de momentos en losas rectangulares, como los que se tienen tabulados en los reglamentos actuales de diseño, se realizaron calibraciones con pruebas de laboratorio como las reportadas por Bach y Graf (1915), Mayes et al. (1959), Hatcher et al. ((1960), Hatcher et al. (1961),Vanderbilt et al. (1961), Gamble et al. (1961), Jirsa et al. (1962), Girolami et al. (1970), Foster et al. (2004), Casadei et al. (2005) y Galati et al. (2008), entre otros. Es de interés señalar que las referencias mencionadas son la base de las investigaciones y aplicaciones que se realizan en el análisis y diseño de losas rectangulares actualmente. Sin embargo, al revisar normatividades vigentes de diseño de losas como el ACI 318-11, Eurocódigo 2 (EC2), y el Reglamento de Construcciones para el Distrito Federal (RCDF-04) y sus Normas Técnicas Complementarias para Diseño y Construcción de Estructuras de Concreto (NTCC-04), no se encontraron recomendaciones para el análisis y diseño de losas con geometría diferente a la rectangular.

Para determinar los momentos de diseño en losas rectangulares, el ACI 318-11 sugiere los métodos directo y el equivalente, las NTCC-04 proporcionan una tabla de coeficientes y el EC2, aunque sólo recomienda utilizar algún método válido, también proporciona como opción el uso de una tabla de coeficientes reportada en IStructE (2006). Es de interés mencionar que la tabla de coeficientes que proporcionan las NTCC-04 fue tomada del ACI 318-63, los cuales, como lo demuestran Gutiérrez-Morgado y Juárez-Luna (2012), en algunos casos subestiman los momentos de diseño.

El resto del artículo se describe a continuación. En la siguiente sección se presentan los modelos constitutivos del concreto y el acero utilizados en los análisis no lineales. La tercera sección presenta la validación del método de los elementos finitos y de los modelos constitutivos con el modelado numérico de un experimento de laboratorio reportado en la literatura. Después se presentan los modelos numéricos de las losas en estudio, determinándose la distribución del acero con análisis elásticos y las curvas carga contra desplazamiento al centro, resultado de los análisis no lineales del material. En la quinta sección se determinan los coeficientes de diseño para las losas con las geometrías estudiadas. La sexta sección muestra el inicio y la evolución de las trayectorias del agrietamiento de losas empotradas y simplemente apoyadas. En la última sección se presentan las conclusiones derivadas de esta investigación.

Modelos constitutivos

Para modelar el comportamiento del concreto se utilizó un modelo de daño con diferente umbral en tensión y en compresión, el cual fue propuesto por William y Warnke (1975), consistente con el comportamiento constitutivo experimental reportado en la literatura, como se muestra en la figura 2a. El concreto considera el ablandamiento por deformación después de alcanzar el esfuerzo último; aunque existen otras funciones de ablandamiento como la exponencial, bilineal, trilineal, etcétera, en los análisis no lineales, este ablandamiento se idealizó con una función lineal (figura 2b). Para modelar el acero de refuerzo se utilizó un modelo de plasticidad con superficie de fluencia de Von Mises, en el cual se supone el mismo umbral en tensión y compresión, como se muestra en la figura 3a; el endurecimiento se idealizó como una función bilineal mostrada en la figura 3b.

Los modelos no lineales se realizaron en el programa de elementos finitos Ansys 13.0. El concreto se discretizó con el elemento finito sólido hexaedro de ocho nodos del tipo solid 65 mostrado en la figura 4a, el cual tiene la capacidad de simular de manera distribuida el agrietamiento en tensión y aplastamiento en compresión, en tanto que el acero se discretizó con el elemento finito unidimensional lineal link 180 mostrado en la figura 4b, al que se le asignó el modelo de plasticidad de Von Mises. Ambos elementos finitos tienen tres grados de libertad traslacional por nodo, así como la capacidad de representar deformaciones grandes.

Validación

En este artículo no se realizaron modelos experimentales de laboratorio debido a que está fuera del objetivo del trabajo, por lo que, para validar los modelos

constitutivos del concreto y el acero, se modeló uno de los seis especímenes probados por Girolami et al. (1970). El espécimen consistió de una losa de 1.829 * 1.829 m y 0.044 m de espesor, en la que se aplicaron únicamente cargas verticales (figura 5a). Cada carga vertical se aplicó sobre cuatro gatos, distribuida sobre 16 placas, mediante cuatro árboles de carga como se muestra en la figura 5b. Cada placa de distribución se centró sobre los huecos del espécimen de la losa. Las placas de distribución de cargas eran cuadradas de 0.2038 m por lado y de 0.0254 m de ancho; cada placa tenía un asiento esférico cóncavo para una junta esférica cóncava. El concreto tenía un módulo elástico E_c =19.90 GPa y su resistencia a compresión era de σ_uc =31.026 MPa, el acero de refuerzo tenía un módulo elástico E_a = 206 GPa y un esfuerzo de fluencia σ_y = 330.95 MPa.

El sistema de losa y vigas perimetrales tenía acero en la parte superior e inferior, diseñado para soportar una carga uniformemente distribuida de 7.182 kPa. La distribución del acero de refuerzo de la losa se muestra en la figura 6a y el de las vigas en la figura 6b. El hecho que las vigas perimetrales tengan bastante acero por cortante, estribos del #7, se debe a que a cinco de los seis especímenes probados por Girolami et al. (1970) se les aplicó cargas en su plano, además de las cargas perpendiculares a su superficie; con este acero por cortante se buscó evitar el daño en las vigas y transmitirlo al plano de las losas.

El análisis no lineal de esta losa se realizó considerando los dos ejes de simetría del espécimen, por lo que sólo se modeló un cuarto del espécimen, restringiendo respectivamente los grados de libertad perpendiculares a cada eje de simetría. La malla de acero de refuerzo de la losa se muestra en la figura 7a, la cual se estructuró de tal forma que los elementos lineales coincidan con los bordes de los elementos sólidos, mientras que la malla del concreto de la losa y las vigas con la malla del acero de refuerzo embebido se muestra en la figura 7b, los nodos del acero de refuerzo coinciden con los nodos de los elementos sólidos. En éste y en el resto de los ejemplos numéricos del artículo se consideró el acoplamiento perfecto entre los grados de libertad del acero y del concreto.

Ejemplos numéricos

Para el estudio de losas, inicialmente se realizaron análisis elásticos con elementos finitos placa, los cuales tienen un grado de libertad traslacional y dos angulares en cada nodo; posteriormente se realizaron análisis no lineales con elementos sólidos con ocho nodos y barras con dos nodos, ambos elementos con tres grados de libertad traslacionales por nodo. Con los análisis elásticos se determinan las distribuciones de momento, las cuales son de utilidad para proponer el acero de refuerzo superior e inferior; mientras que con los análisis no lineales se determinaron las trayectorias de agrietamiento, así como la evolución de los coeficientes de diseño al existir una redistribución de esfuerzos debido al agrietamiento.

En los ejemplos elásticos se utilizó un mallado con elementos finitos placa para obtener la distribución de momentos de las losas circular, elíptica y triangular mostradas respectivamente en las figuras 10, 12 y 14. En estos modelos se utilizó un módulo elástico equivalente al concreto reforzado de E = 21.673 GPa y relación de Poisson de v = 0.25 bajo condiciones de apoyo empotrado, pues las condiciones de apoyo simple proporcionan distribuciones de momentos que sólo requieren acero en la parte inferior de la losa.

En los modelos no lineales se utilizaron elementos sólidos en 3D para el concreto simple con un módulo elástico E_c = 17GPa, relación de Poisson v = 0.20, esfuerzo último a tensión σ_ut = 5MPa y a compresión σ_uc = 25 MPa. Para el acero de refuerzo se utilizó un módulo elástico E_a =2.05 * 10⁵ MPa, relación de Poisson v = 0.30 y esfuerzo de fluencia σ= 410 MPa. En estos modelos el espesor de la losa es 0.10 m con un recubrimiento de 0.025 m para ambos lechos de acero. Los modelos se analizaron en condiciones de apoyo simple y empotrado, en los cuales se aplicaron incrementos de carga distribuida sobre la superficie. En la generación del mallado se consideraron los dos ejes de simetría de las losas, modelando sólo una cuarta parte de su geometría para reducir el costo computacional.

Distribución del acero de refuerzo

La obtención de momentos elásticos en losas circulares se realizó mediante un modelo numérico con radio r = 2 m, en el que el mallado se estructuró con elementos cuadriláteros placa, evitando tener elementos triangulares en el centro con relaciones de aspecto relativamente grandes que pudieran producir inestabilidades numéricas. A partir de la distribución de momentos mostrada en la figura 10, se identificaron las zonas de momento negativo en la losa circular, la cual está delimitada por la línea punteada y los bordes, obteniéndose una longitud del acero superior en el borde de L_neg = r/2; es de interés mencionar que esta longitud es análoga a una losa cuadrada en la que se recomienda una longitud de acero negativo de un cuarto del claro, que para el caso de losa circular es un cuarto del diámetro. Se propuso una distribución típica del acero de refuerzo como se muestra en la figura 11, en la que el acero positivo se coloca en la parte inferior de la losa y el acero negativo ubicado en la parte superior en los bordes, el cual baja al emparrillado inferior en la zona central mediante un columpio o bayoneta realizado en la línea punteada.

La distribución de momentos en la sección elíptica se determinó mediante modelos en los que se varió la relación del eje menor b respecto al eje mayor a, i.e., b/a=0.5,0.6,0.7,0.8 y 0.9, respectivamente, considerando constante el valor del eje mayor a = 4 m; al igual que en las losas circulares. En el mallado de las losas elípticas se utilizaron elementos cuadriláteros placa, de tal forma que los nodos de los elementos coincidan con los nodos del acero de refuerzo. A partir de la distribución de momentos mostrada en la figura 12 se identificaron las zonas de momento negativo en la losa elíptica, teniendo que la longitud del acero superior es en el eje mayor l_neg= a/3 y en el eje menor l_neg= b/3. Con base en esta distribución de momentos, se propuso el armado típico mostrado en la figura 13, en el que en la zona central se tiene sólo acero en el lecho inferior, y en el borde, delimitado por la línea punteada, se tiene acero en las partes superior e inferior.

La distribución de momentos de la sección triangular se determinó mediante un modelo numérico, el cual tiene una altura a = 4 m y de base 2a/√3 y un mallado estructurado con elementos cuadriláteros placa, que evita singularidades numéricas. En la figura 14 se muestran las longitudes correspondientes al momento negativo en ambas direcciones: en dirección horizontal se tiene una longitud de a/3, mientras que en dirección vertical se tiene una longitud de a/7 en la base y a/2 en parte superior. Con base en estas longitudes, en la figura 15 se muestra el armado propuesto con acero sólo en la parte inferior de la zona central, delimitada por la línea punteada y con acero en el lecho superior e inferior en los bordes.

Análisis no lineal

El análisis no lineal de las losas elípticas se realizó considerando relaciones eje menor entre eje mayor b/a = 0.5, 0.6, 0.7, 0.8 y 0.9. El acero de refuerzo está formado por barras de acero del número 3 separadas a cada 0.20 m, que incluye el acero superior, inferior y bastones, la malla del refuerzo discretizada con elementos lineales se muestra en la figura 18a. La malla del concreto de la losa con el acero embebido se muestra en la figura 18b, observe que la malla de concreto se estructuró para que los nodos de los elementos sólidos coincidan con los del acero de refuerzo y para tener sólo elementos hexaedros. Las curvas de la variación de la carga distribuida con respecto al desplazamiento en el centro de la losa se muestran en la figura 19, donde se observa que la carga necesaria para desarrollar un desplazamiento al centro de la losa es inversamente proporcional a la relación b/a en ambas condiciones de apoyo; además se observa que esta carga en las losas simplemente apoyadas es aproximadamente 20% de la carga de las losas empotradas para desarrollar el mismo desplazamiento en el centro.

Respecto al análisis de la losa triangular, el acero de refuerzo consiste de barras de acero del número 3 a cada 0.20 m, su mallado se realizó con elementos finitos lineales como se muestra en la figura 20a, la malla del concreto de la losa con el acero embebido se muestra en la figura 20b, se realizó el mallado de la mitad del modelo por sólo tener un eje de simetría en esta losa triangular. Las curvas de la variación de la carga distribuida con respecto al desplazamiento en el centro de la losa se muestran en la figura 21, donde nuevamente se observa que las condiciones de apoyo influyen en el comportamiento global, de tal forma que la losa triangular simplemente apoyada requiere de una carga aproximada de 20% en comparación con la losa empotrada para desarrollar el mismo desplazamiento en el centro.

Determinación de coeficientes de diseño

donde σ_x es la función que describe el esfuerzo en la dirección x, z es el eje de referencia y t es el peralte de la losa. En los análisis elásticos presentados en la sección de distribución del acero de refuerzo no fue necesario utilizar la ecuación (1) para determinar los momentos, puesto que el elemento placa los proporcionó directamente.

En los análisis no lineales, el cálculo de los momentos flexionantes con la ecuación (1) se realizó integrando en ocho tramos la función de esfuerzos, la cual se describe por funciones de interpolación lineales, como se muestra en la figura 22 en el ancho de la losa; el orden de las funciones de interpolación de esfuerzos es consistente con el de la derivada de las funciones de interpolación de los desplazamientos nodales.

Una vez obtenido el momento, los coeficientes se calcularon con la siguiente ecuación:

donde α es el coeficiente de diseño, w carga por unidad de área, l longitud de referencia (r para la losa circular, bpara la elíptica y a para la triangular) y M el momento calculado con la ecuación (1).

Los coeficientes calculados con análisis no lineales presentan una variación, decremento o incremento, con respecto a los coeficientes elásticos, tomándose los valores máximos de la variación respecto al incremento de carga; esto se debe a la influencia del agrietamiento en el concreto y el endurecimiento del acero de refuerzo en las zonas más esforzadas de la losa. Por lo tanto, se propone utilizar aquellos coeficientes con valores máximos con el fin de garantizar la seguridad estructural de las losas, para que tanto el concreto como el acero de refuerzo brinden una resistencia adecuada en el intervalo elástico.

Los valores de los coeficientes calculados en esta sección deben multiplicarse por 10^-4wl² para obtener momentos de diseño por unidad de ancho. El índice I indica la condición de colado monolíticamente, mientras que el índice II indica la condición de colado no monolíticamente. Es en esta última condición que los coeficientes negativos en los bordes tienen un valor de cero, pues sólo existe momento positivo en toda la losa debido a la condición de apoyo.

La obtención de los coeficientes en losas circulares se realizó en el centro, donde se presentan los momentos máximos positivos, y en los extremos A y B, donde se presentan los momentos máximos negativos, los cuales se muestran en la figura 23. Los valores de los coeficientes son idénticos en A, B y en cualquier punto sobre el borde, por tratarse de una estructura que se puede considerar como axisimétrica.

La variación de los coeficientes de diseño calculados para la losa circular en la condición de apoyo simple se muestran en la figura 24a, la cual es idéntica en cualquier dirección por ser una estructura axisimétrica; mientras que en el caso de la losa empotrada se presenta el momento máximo negativo en los extremos y positivo en el centro como se muestra en la figura 24b. En ambas condiciones de apoyo se observa que al iniciar el agrietamiento los coeficientes decrecen en el centro y en los bordes debido al agrietamiento en el concreto, posteriormente los coeficientes se incrementan debido a la redistribución de esfuerzo y plastificación del acero, sólo los coeficientes en los extremos superan su valor en el intervalo elástico y los del centro quedan debajo de este valor. De los coeficientes mostrados en la figura 24 se tomaron los valores de los coeficientes de diseño elásticos, los cuales se muestran en la tabla 2, para el caso de los coeficientes en el intervalo inelástico se tomaron los correspondientes a los valores máximos, que para el centro corresponden a los mismos valores elásticos y para los bordes se tiene un valor mayor a este intervalo como se muestra en la tabla 3.

Losa elíptica

La obtención de los coeficientes en las losas elípticas se realizó en los puntos donde se obtienen los valores de momento máximo negativo, puntos A y B, y momento máximo positivo, centro (figura 25). En losas con esta geometría se tienen diferentes magnitudes de los coeficientes en las dos direcciones en el centro; en el punto A, los coeficientes tienen valores para determinar el momento en la dirección horizontal y para el punto B, el momento en la dirección vertical.

Las figuras 26 a 30 (27, 28, 29) muestran la variación del coeficiente para las relaciones 0.5, 0.6, 0.7, 0.8 y 0.9, respectivamente. La variación de los coeficientes en el centro de las losas es mayor en la dirección y para la condición de apoyo simple, este mismo comportamiento se observa en el intervalo elástico y en el no lineal, aunque en el segundo intervalo, la magnitud de los coeficientes decrece y posteriormente vuelven a incrementarse sin alcanzar la magnitud del intervalo elástico. Para el caso de losas con empotramiento, la magnitud máxima de los coeficientes se presentan en el extremo B en el intervalo elástico, posteriormente, en el intervalo no lineal, la magnitud de estos coeficientes decrece y después muestran un incremento, particularmente para las relaciones 0.5 a 0.7; en el resto de las relaciones el momento máximo se presentó en el centro en la dirección x en el intervalo no lineal, esto se debe a la influencia de la redistribución de esfuerzos por el agrietamiento del concreto y la plastificación del acero de refuerzo. Como resultado de este análisis, los coeficientes elásticos e inelásticos se muestran respectivamente en las tablas 4 y 5.

Losas triangulares

En losas triangulares, las zonas A, B y C presentan momento negativo, por el contrario, en el centro presentan momento positivo, como se muestran en la figura 31. Los momentos en C se calcularon en la dirección θ = 30° respecto al eje horizontal, puesto que los esfuerzos principales en C se presentaron en esta dirección. Las zonas de interés para la condición simplemente apoyada son en el centro, porque sólo se presenta momento positivo, mientras que para la condición empotrada son de interés todos los puntos mostrados en la figura 31 por presentarse momentos con ambos signos.

En la losa triangular simplemente apoyada, la magnitud de los coeficientes son mayores en el centro en la dirección horizontal en el intervalo elástico, condición que cambia cuando se está en el intervalo no lineal, puesto que los momentos máximos se presentan en el centro en dirección y, como se muestra en la figura 32a. Por otra parte, en la condición de empotramiento, el momento máximo se presenta en la zona C en ambos intervalos, aunque en principio existe una reducción del coeficiente, posteriormente éste se incrementa superando su magnitud respecto al intervalo elástico (figura 32b). Con base en lo anterior, los coeficientes para los intervalos elásticos e inelásticos obtenidos para losa triangular se muestran, respectivamente, en las tablas 6 y 7.

Trayectorias de agrietamiento

Las trayectorias de agrietamiento son de utilidad para identificar zonas de mayor concentración de esfuerzos, con las que se determinan cualitativamente la ubicación del acero de refuerzo. Dentro de la mecánica computacional existen principalmente tres modelos para representar el agrietamiento:

]]> 1. La mecánica de fractura, que localiza en una línea o superficie las grietas, pero tiene la desventaja de necesitar grietas iniciales, proceso de remallado costoso y se aplica sólo a materiales frágiles;

2. El agrietamiento distribuido, que las grietas aparecen en cualquier parte del modelo, por lo que no necesita grietas iniciales. Una de sus ventajas principales en que no necesita remallado, pero presenta problemas de atoramiento de esfuerzos; y

3. El modelo de discontinuidades interiores, que toma las ventajas de los dos anteriores, al no necesitar el remallado e introducir las discontinuidades (grietas) en cualquier lugar y dirección dentro de los elementos finitos. En estos dos últimos no se observa como tal una discontinuidad física (grieta), sino que se presentan zonas con desplazamientos mayores y, en consecuencia, deformaciones con gradientes grandes.

Es de interés mencionar que en este artículo el agrietamiento se modeló como un agrietamiento distribuido, asociado al estado de esfuerzos en los puntos de integración Gaussiana del elemento, donde el material alcanza la superficie de falla, por lo que no se muestran discontinuidades físicas (grietas), sino que se presentan zonas con grandes desplazamientos y, en consecuencia, concentración de gradientes grandes de deformaciones donde se degrada el material. El software ANSYS, donde se realizaron estos análisis, despliega círculos en los puntos de integración que se agrietan o aplastan; el agrietamiento se muestra con un círculo en dirección del plano de agrietamiento, mientras que el aplastamiento se muestra con un octaedro. Si la grieta se abre y posteriormente se cierra, se despliega una X dentro del círculo. En las figuras mostradas en esta sección se delimitaron las zonas que contienen los puntos de integración dañados, las cuales se muestran en color rojo donde los estados de esfuerzo están dominados por la tensión y en azul los que están en compresión.

Losas circulares

En losas circulares empotradas en sus bordes, el agrietamiento inicia en todo el borde por la parte superior de la losa debido a que es axisimétrica, como se muestra en la figura 33a, posteriormente, el agrietamiento se presenta al centro del claro en la parte inferior (figura 33b). Se observó que el agrietamiento evoluciona, de manera que en la parte superior éste se propaga de los bordes hacia la parte central, no obstante, en el centro se presenta un aplastamiento y en la parte inferior, el agrietamiento se propaga del centro hacia los bordes.

Cuando las losas circulares están simplemente apoyadas en sus bordes, el agrietamiento inicia en el centro de la parte inferior, el cual, posteriormente, se propaga hacia los bordes, como se muestra en la figura 34b, y en la parte superior se presenta un aplastamiento incipiente como se muestra en la figura 34a. El que nada más exista momento positivo en la losa demuestra que sólo existen esfuerzos a compresión en la parte superior de la losa, por lo que aparecen aplastamientos, y esfuerzos de tensión en la parte inferior, con lo que aparecen agrietamientos.

Losas elípticas

En las losas elípticas empotradas, el agrietamiento inicia en la parte superior de los bordes aproximadamente perpendiculares al eje mayor, como se muestra en la figura 35a, posteriormente, el agrietamiento aparece en el centro de la parte inferior de la losa, como se muestra en la figura 35b. Al incrementar la carga distribuida, el agrietamiento se propaga análogamente a la losa circular, pues en la parte superior de la losa, el agrietamiento se propaga de los bordes hacia la parte central generándose aplastamiento en el centro; mientras que en la parte inferior de la losa el agrietamiento se propaga del centro hacia el borde.

Es en la condición de apoyo simple, que las losas elípticas comienzan a agrietarse en su parte inferior, como se observa en la figura 36b, posteriormente, éste se propaga del centro hacia los bordes. Al incrementar la carga distribuida, un incipiente aplastamiento aparece en la parte superior de la losa como se observa en la figura 36a.

Losas triangulares

El agrietamiento en la losa triangular empotrada inicia en los tres bordes de la parte superior de la losa, como se muestra en la figura 37a, donde se tienen estados de esfuerzos negativos; posteriormente, el agrietamiento aparece en la zona central de la parte inferior de la losa, como se muestra en la figura 37b. Al incrementar la carga distribuida, el agrietamiento se propaga del centro de los bordes hacia los vértices en la parte superior, en cuanto que éste se propaga del centro a los vértices en la cara inferior.

Cuando la losa triangular está simplemente apoyada en sus bordes, el agrietamiento inicia en el centro de la parte inferior de la losa, propagándose hacia los vértices conforme se incrementa la carga distribuida, como se muestra en la figura 38b. Por otra parte, debido a los estados de esfuerzos a compresión en la parte superior se presenta aplastamiento, como se muestra en la figura 38a.

En general, cuando las losas se encuentran empotradas, el agrietamiento inicia en los bordes de la parte superior, posteriormente en el centro de la parte inferior. En las losas simplemente apoyadas, el agrietamiento inicia en la cara inferior de la losa, al centro, propagándose en todo el sistema; posteriormente se presenta un aplastamiento incipiente al centro de la parte superior. El que el agrietamiento sea mayor en la cara inferior de las losas respecto al aplastamiento que aparece en la cara superior se debe directamente a la redistribución de esfuerzos y a la diferencia de resistencia que tiene el concreto al someterse a tensión que es de 5% a 20% de la resistencia a compresión.

El comparar la curva carga contra desplazamiento de los resultados numéricos con los experimentales reportados en la literatura se válida el MEF y los modelos constitutivos utilizados para el estudio de losas con geometría triangular, circular y elíptica, pues a pesar de que la curva experimental muestra una recuperación de desplazamiento, se considera que éste posiblemente se debió a un deslizamiento de los instrumentos de medición, pues ante una carga vertical incremental como la del experimento no puede presentarse una recuperación de desplazamientos al centro.

Se descarta el método de líneas de fluencia para determinar las trayectorias de agrietamiento y los coeficientes de diseño o la carga de colapso, pues a pesar de que este método proporciona trayectorias de agrietamiento, no lo proporciona en la parte superior o inferior de la losa; además, estos métodos sólo proporcionan un valor de la carga de colapso, que puede ser mayor o menor a la carga real de colapso. Por lo anterior, se utilizó la mecánica numérica que proporciona el inicio y evolución del agrietamiento, así como la variación de los momentos en el intervalo no lineal.

El armado propuesto para cada geometría de losa, con base en análisis elásticos de placas homogéneas de concreto reforzado resulta adecuado, pues se obtiene una buena aproximación de la distribución de momentos y de la localización e inicio del agrietamiento en el sistema. Estas longitudes de acero negativo y positivo se recomiendan para el diseño y construcción de losas con estas geometrías. De las gráficas carga distribuida contra desplazamiento al centro de la losa con las tres geometrías estudiadas, se observó que la carga distribuida sobre las losas con apoyo simple es aproximadamente 20% de la magnitud de la carga sobre las losas empotradas en sus apoyos para desarrollar un mismo desplazamiento en el centro de la losa.

Los coeficientes elásticos obtenidos con los modelos de elementos finitos sólidos y barra son consistentes en el intervalo elástico con aquellos obtenidos con elementos finitos placa. Después, cuando se agrieta el concreto y/o fluye el acero de refuerzo, el valor de estos coeficientes cambia debido a la redistribución de esfuerzos.

En las tablas de coeficientes, resultado de los análisis no lineales, se proporcionan los valores máximos como los de diseño, con la finalidad de garantizar un comportamiento estructural adecuado. Sin embargo, se recomienda utilizar los coeficientes del análisis elástico, pues se puede observar que en la mayoría de los casos los coeficientes decrecen debido al agrietamiento y luego se incrementan por la plastificación del acero, teniendo en la mayoría de los casos el valor de los coeficientes elásticos como máximo.

En general, el agrietamiento en las losas empotradas inicia en los bordes de la cara superior, posteriormente en el centro de la parte inferior, propagándose al incrementar la carga; mientras que en las losas simplemente apoyadas inicia en el centro del claro en la cara inferior, propagándose hacia los bordes.

Agradecimientos

A la Universidad Autónoma Metropolitana por las facilidades proporcionadas en la realización de este trabajo y al patrocinio proporcionado por el Programa de Mejoramiento del Profesorado. El primer autor agradece al proyecto "182736 Análisis y diseño de losas de concreto" auspiciado por la Secretaría de Educación Pública (SEP) y al Consejo Nacional de Ciencia y Tecnología (CONACyT). El segundo autor agradece la beca al CONACyT por la beca para la realización de estudios de maestría.

Referencias

Ansys 13.0.1, Ansys Inc., Estados Unidos, 2010.         [ Links ]

Bach C., Graf O. Tests with Simply Supported, Quadratic Reinforced Concrete Plates (en Alemán Versuchemitallseitigaufliegenden, quadratischen andrechteckigeneisenbetonplatten), Deutscher Ausschuss für Elsenbeton, 30, Berlin, 1915.         [ Links ]

ACI 318-63. Building Code Requirements for Structural Concrete (ACI-318-63) and Commentary, Farmington Hills, MI, USA, American Concrete Institute, 1963.         [ Links ]

ACI 318-11. Building Code Requirements for Structural Concrete (ACI-318-11) and Commentary, Farmington Hills, MI, USA, American Concrete Institute, 2011.         [ Links ]

Casadei P., Parretti R., Nanni A., Heinze T. In Situ Load Testing of Parking Garage Reinforced Concrete Slabs: Comparison Between 24h and Cyclic Load Testing. Practice Periodical on Structural Design and Construction. ASCE, volumen 10 (número 1), febrero de 2005: 40-48 [en línea]. Disponible en: doi: http://dx.doi.org/10.1061/(ASCE)1084-0680(2005)10:1(40)        [ Links ]

EC2. Eurocode 2: Design of Concrete Structures, British Standards Institution, 1992.         [ Links ]

Foster S.J., Bauley D.G., Burgess I.W., Plank R.J. Experimental Behaviour of Concrete Floor Slabs at Large Displacements. Engineering Structures, volumen 26, 2004:1231-1247. [en línea]. Disponible en: doi: http://dx.doi.org/10.1016/j.engstruct.2004.04.002        [ Links ]

Galati N., Nanni A., Tumialan J.G., Ziehl P.H. In-Situ Evaluation of Two Concrete Slab Systems, I: Load Determination and Loading Procedure. Journal of Performance of Constructed Facilities, ASCE, volumen 22 (número 4), 2008: 207-216. [en línea]. Disponible en: doi: http://dx.doi.org/10.1061/(ASCE)0887-3828(2008)22:4(207)        [ Links ]

Gamble W.L., Sozen M.A., Siess C.P. An Experimental Study of a Two-Way Floor Slab, Structural Research Series Núm. 211, Departamento de Ingeniería Civil, Universidad de Illinois, 1961, p. 326.         [ Links ]

Girolami A.G., Sozen M.A., Gamble W.L. Flexural Strength of Reinforced Concrete Slabs with Externally Applied in-Plane Forces, Report to the Department of Defense, Universidad de Illinois, Urbana, Illinois, 1970, p. 176.         [ Links ]

Gutiérrez-Morgado P., Juárez-Luna G. Coeficientes numéricos en losas rectangulares aisladas, Memorias, XVIII, en: Congreso Nacional de Ingeniería Estructural, Acapulco, Guerrero, CDROM, Artículo Núm. 15-10, 1-12, octubre-noviembre, 2012.         [ Links ]

Hatcher D.S., Sozen M.A, Siess C.P. An Experimental Study of a Quarter-Scale Reinforced Concrete Flat Slab Floor, Structural Research Series Núm. 200, Departamento de Ingeniería Civil, Universidad de Illinois, 1960, p. 288.         [ Links ]

Hatcher D.S., Sozen M.A., Siess C.P. A Study of Tests on a Flat Plate and a Flat Slab, Structural Research Series Núm. 217, Departamento de Ingeniería Civil, Universidad de Illinois, 1961, p. 504.         [ Links ]

IStructE. Manual for the Design of Concrete Building Structures to Eurocode 2, Institution of Structural Engineers, The Concrete Center, 2006, p. 141.         [ Links ]

Jirsa J.O., Sozen M.A., Siess C.P. An Experimental Study of a Flat Slab Floor Reinforced with Welded Wire Fabric, Structural Research Series Núm. 248, Departamento de Ingeniería Civil, Universidad de Illinois, 1962, p. 188.         [ Links ]

Marcus H. Die Vereinfachteberechnungbiegsamerplatten, Julius Springer, Berlín, 1929, p. 126.         [ Links ]

Mayes G.T., Sozen M.A., Siess C.P. Tests on a Quarter-Scale Model of a Multiple-Panel Reinforced Concrete Flat Plate Floor, Structural Research Series Núm. 181, Departamento de Ingeniería Civil, Universidad de Illinois, 1959.         [ Links ]

Newmark N.M., Siess C.P. Proposed Design Specifications for Two-Way Floor Slabs. Journal Proceedings, American Concrete Institute, volumen 46 (número 4), abril de 1950: 597-607.         [ Links ]

NTCC-04. Normas Técnicas Complementarias para Diseño de Estructuras de Concreto, Gaceta Oficial del Distrito Federal, Tomo II, No. 103-BIS, octubre de 2004, p. 107.         [ Links ]

RCDF-04. Reglamento de Construcciones del Departamento del Distrito Federal, Gaceta Oficial del Departamento del Distrito Federal, enero de 2004, p. 40.         [ Links ]

Siess C.P., Newmark N.M., Moments in Two-Way Concrete Floor Slabs, Bulletin Series Núm. 385, Universidad de Illinois, Urbana, Illinois, febrero de 1950, p. 124.         [ Links ]

Vanderbilt M.D., Sozen M.A., Siess C.P. An Experimental Study of a Reinforced Concrete Two-Way Floor Slab with Flexible Beams, Structural Research Series Núm. 228, Departamento de Ingeniería Civil, Universidad de Illinois, 1961, p. 188.         [ Links ]

William K.J., Warnke E.D. Constitutive Model for the Triaxial Behavior of Concrete, Proceedings of the International Association for Bridge and Structural Engineering, ISMES, Bergamo, Italy, 19, 1975, p. 174.         [ Links ]

Westergaard H.M., Slater W.A. Moments and Stresses in Slabs. Proceedings of the American Concrete Institute, volumen 17 (número 2), febrero de 1921: 415-538.         [ Links ]

Semblanza de los autores

Gelacio Juárez-Luna. Es ingeniero civil y maestro en ciencias con especialidad en estructuras, ambos por el Instituto Politécnico Nacional. Obtuvo el grado de doctor en ingeniería por la Universidad Nacional Autónoma de México en 2006. Es profesor-investigador en la Universidad Autónoma Metropolitana desde 2009. Su trabajo de investigación está orientado al desarrollo y aplicación de elementos finitos para resolver problemas estructuras en su evolución al colapso.