Servicios Personalizados
Revista
Articulo
Indicadores
-
Citado por SciELO
-
Accesos
Links relacionados
-
Similares en SciELO
Compartir
Ingeniería sísmica
versión impresa ISSN 0185-092X
Ing. sísm no.89 Ciudad de México jul./dic. 2013
Artículos
Envolvente de resistencia lateral de piso para estructuras de mampostería confinada*
Sergio Alcocer M.1, Hugo Hernández B.2, Harry Sandoval R.3
1 Instituto de Ingeniería de la UNAM, Coordinación de Estructuras y Materiales, Apartado Postal 70-472, México DF, 041510, Tel. (+52 55) 5622-8946. salcocerm@iingen.unam.mx.
2 Universidad Michoacana de San Nicolás de Hidalgo, Edificio de Posgrado en Ingeniería Civil, CU, Morelia, Michoacán, Tel. (+52 443) 322 3500, ext. 4341, hugohbarrios@yahoo.com.mx.
3 Egresado de la Maestría en Estructuras, Programa de Maestría y Doctorado en Ingeniería, UNAM.
Artículo recibido el 28 de noviembre de 2011.
Aprobado para su publicación el 10 de junio de 2013.
Resumen
En este artículo se presenta un procedimiento para calcular la curva envolvente de resistencia lateral de estructuras de mampostería confinada con mediana o baja altura. La curva envolvente idealizada se definió con parámetros obtenidos con pruebas experimentales realizadas en México. Las curvas envolventes de resistencia lateral analíticas son comparadas con las curvas experimentales obtenidas en un modelo probado bajo carga cíclica reversible y una serie en modelos de mampostería probados en la mesa vibradora del II-UNAM. Los resultados sugieren una adecuada predicción del límite elástico, resistencia máxima y estado último. El procedimiento incluye el cálculo de las ductilidades local y global y del factor de comportamiento sísmico.
Palabras clave: curva envolvente, mampostería confinada, ductilidad, coeficiente sísmico, coeficiente sísmico de desempeño, falla por fuerza lateral, falla por carga vertical.
Abstract
This paper presents a procedure for calculating the envelope curve of lateral resistances of confined masonry structures with medium or low height. The idealized envelope curve of lateral resistance was defined by parameters obtained with experimental tests performed in Mexico. The analytical envelope curves of lateral resistance are compared with experimental curves obtained in a model tested under cyclic reversible load and a series of tests on models of masonry in shaking table of the II-UNAM. Results suggest an adequate prediction of the elastic limit, maximum resistance and ultimate state. The procedure includes the calculation of local and global ductility and seismic behavior factor.
Keywords: envelope curve, confined masonry, ductility, seismic coefficient, seismic coefficient of performance, fault lateral force, vertical load fails.
INTRODUCCIÓN
Existen varias propuestas para modelar analíticamente estructuras de mampostería confinada, éstas se pueden clasificar en: (a) Micro-modelado, (b) Macro-modelado y (c) modelado por medio de curvas envolvente de comportamiento histerético. El micro-modelado de los elementos que forman un muro de mampostería se basa principalmente en técnicas numéricas como el Elemento Finito y Elementos de Frontera (Ishibashi y Kastumata, 1994). Debido a su alto costo computacional y el uso de software especializado e incertidumbre en las propiedades de los materiales, aún no es justificable que dichas metodologías se apliquen de manera práctica.
Actualmente es posible, con cualquier software comercial, modelar estructuras de mampostería con macro-modelos estructurales. Los macro-modelos o elementos estructurales más populares son las barras y puntales a compresión y el método de la columna ancha (Holmes, 1961; Kadir, 1974; Chrysostomou et al., 1992; Mosalamet al., 1996; Madanet al., 1997; Crisafulli, 1997; Tomazevic, 1999). También se han propuesto algunas modificaciones a dichos elementos (Beyer et al., 2008; Zúñiga y Terán, 2008). En México, Bazán (1980) demostró la capacidad de la metodología de la columna ancha para representar adecuadamente la rigidez lateral elástica de muros de mampostería confinada sujetos a cargas laterales. En Zúñiga y Terán (2008) se propone una modificación a dicho modelo con el fin de estimar la respuesta no lineal de estructuras de mampostería, para lo cual calculan curvas de capacidad y umbrales de desplazamientos asociados a los estados límite. Aun cuando los autores encuentran una elevada variabilidad entre los resultados experimentales estudiados, consideran que la metodología que proponen para calcular las curvas de capacidad de los especímenes ofrece estimaciones razonablemente conservadoras y que la resistencia asociada al primer agrietamiento como la rigidez elástica se obtienen con suficiente precisión (Figura 1a).
En otro trabajo (Siete, 2011) se propone calcular curvas de fragilidad validando las curvas de capacidad calculadas con el método de la columna ancha ya que el método de puntales a compresión, usando el software IDARC-2D (Valles et al., 1996) no predice adecuadamente los resultados experimentales. El argumento del autor es que la curva calculada con el método de la columna ancha captura adecuadamente el comportamiento inicial en rigidez y resistencia, pero después del agrietamiento las curvas se separan de los resultados experimentales (Figura 1b). Alcocer et al., (2004) al comprar los valores experimentales en mesa vibradora con los obtenidos con un modelo simple, obtuvieron buenos resultados en cuanto al cortante basal pero encontraron ciertas diferencias en los ángulos de distorsión asociados, por lo que sugieren mejorar el modelado no lineal de la mampostería confinada, principalmente con el fin de capturar adecuadamente el confinamiento especialmente para grandes ángulos de distorsión. Esto también se evidencia en la Figura 2, en donde para una estructura denominada 3D (Sánchez et al., 2010) se comparan las curvas envolventes analíticas calculadas por diferentes autores y metodologías, con la curva envolvente experimental. Las gráficas indican que las diversas metodologías predicen adecuadamente la pendiente inicial que representa el comportamiento elástico del muro, pero todas difieren después de este punto.
En cuanto a los modelos analíticos basados en una curva envolvente de comportamiento histerético de muros de mampostería, Meli (1979) propuso una curva trilineal que relaciona el esfuerzo cortante con la distorsión asociada. Moroni, Astroza y Tavonatti (1994) a partir de estudios experimentales realizados en estructuras de mampostería confinada, propusieron un modelo trilineal para caracterizar el comportamiento no-lineal de la mampostería ante cargas laterales. Flores (1995) con base en los resultados experimentales obtenidos en el Centro Nacional de Prevención de Desastres (CENAPRED) propuso un modelo analítico para representar la rigidez de los lazos de histéresis en los ciclos de carga y descarga. El modelo de Flores se basa en la teoría de fricción y es similar al de Meli (1979) excepto que el tercer tramo de la curva considera degradación de resistencia. Tomazevic propuso representar el comportamiento histerético de muros de mampostería sujetos a carga lateral en combinación con carga vertical usando curvas envolventes de tipotrilineal (Tomazevic y Klemenc, 1997a) y bilineal (Tomazevic, 1999). En la Figura 3 se muestran esquemáticamente dichas curvas, los valores que las definen se pueden encontrar en las referencias antes citadas.
Las recomendaciones propuestas en Jean y Cesín (2000) corroboran que el método de la columna ancha es adecuado para el análisis y diseño elástico de estructuras de mampostería confinada y/o reforzada interiormente. De hecho, el método de la columna ancha es usado por Tena et al. (2010) para validar el procedimiento de diseño conocido como Método Simplificado. De lo comentado anteriormente se puede aceptar que el uso de la metodología de la columna ancha es adecuado para el cálculo de la rigidez elástica incluyendo deformaciones por cortante, pero para el cálculo de curvas de capacidad de estructuras de mampostería presenta inconsistencias, principalmente después de que ocurre el primer agrietamiento.
ENVOLVENTE DE RESISTENCIA LATERAL DE PISO PROPUESTA
La curva de capacidad de una estructura se obtiene a partir de un análisis estático no-lineal bajo desplazamientos laterales monótonamente crecientes. Durante el análisis se aplica una distribución de cargas laterales que mantienen un valor relativo o proporcional a la altura de la estructura pero varían en magnitud. Dicha magnitud se incrementa hasta que la estructura alcanza el desplazamiento lateral asociado a un estado máximo de utilidad que por lo general es la falla o colapso de la estructura. A cada incremento de desplazamiento lateral global de la estructura está asociada una fuerza cortante basal, que graficada versus los desplazamientos produce la curva de capacidad.
Considerando las condiciones de restricción de apoyo de los muros es posible calcular la resistencia a flexión y cortante así como la rigidez inicial efectiva de todos los muros que forman la estructura. Utilizando una curva idealizada de envolvente del comportamiento histerético se puede obtener la contribución de cada muro a la resistencia del entrepiso para calcular la curva envolvente total (Figura 4a). La curva envolvente de entrepiso es la relación entre la resistencia y el desplazamiento relativo del entrepiso, y se obtiene por superposición de las envolventes de resistencia de todos los muros en el entrepiso en consideración (Figura 4b). Esta curva proporciona información sobre el comportamiento y nivel de daño estructural, resistencia máxima y colapso de la estructura, y una vez determinada es posible asociar un índice de daño uniforme en la estructura según se muestra en la Figura 4b (Fischinger, et al., 2004).
La curva envolvente de resistencia (Figura 5) define tres estados límites característicos: (A) Límite Elástico, cuando un muro en el entrepiso alcanza el límite elástico y se agrieta; (B) Estado de Resistencia Máxima; y (C) Estado de Resistencia Última, donde la resistencia del entrepiso se deteriora sobrepasando un límite aceptable establecido. Para fines prácticos, una degradación mayor al 70% de la resistencia máxima define el estado de colapso de la estructura.
En este trabajo al desarrollar las curvas de capacidad se considera que del origen al punto (A) todos los muros están en el intervalo elástico, hasta que se alcanza el punto (B) en el que algún muro llega a su carga máxima. El punto (C) es en donde se considera que la estructura ha alcanzado la falla global.
El procedimiento aquí descrito se emplea para el cálculo de la curva evolvente de resistencia de un entrepiso que forma parte de una estructura de mampostería y que cumple con los requisitos básicos para que sea aplicable el Método de Simplificado de diseño (NTCM-RCDF, 2004). Dicho método se basa en la distribución de fuerzas laterales en estructuras simétricas con diafragmas rígidos cuando la carga lateral se aplica en una sola dirección (Tena et al., 2010).
Con base en los resultados experimentales realizados en modelos de estructuras de mampostería ensayados en mesa vibradora (Tomazevic, 1999; Vázquez, 2005; Barragán, 2005; Arias, 2005) se ha observado que el primer modo de vibrar en cortante prevalece en la respuesta dinámica. Debido a la acción de la losa de entrepiso que funciona como diafragma horizontal rígido, el daño se concentra principalmente en los elementos estructurales del piso inferior, lo anterior se ha observado tanto en pruebas experimentales como después de la ocurrencia de un sismo.
Algunos parámetros que definen el comportamiento sísmico de estructuras de mampostería pueden determinarse si se establece la relación entre la resistencia de entrepiso y su correspondiente desplazamiento, es decir, si se determina la curva envolvente de resistencia de entrepiso (Figura 6a) y además se determina el factor de comportamiento sísmico, Q , definido (Tomazevic et al., 2004) como la relación entre la fuerza lateral que produce el colapso de la estructura y la fuerza que causa la iniciación del daño, es decir, el primer cambio en la rigidez estructural,
donde Ve es la fuerza elástica y Vdu es la carga última de diseño. Para el caso de estructuras de mampostería confinada de piezas macizas se recomienda un valor de Q = 2 (NTCS-RCDF, 2004). Aunque el coeficiente sísmico no sólo depende de la ductilidad sino también de la energía histerética, para estructuras relativamente rígidas (Bazán, 1980; Newmark y Hall, 1982; Tomazevic, 1999) una forma aproximada de obtener este parámetro es:
donde µu es el factor de ductilidad global. Para estructuras de mampostería confinada considerando el valor de Q = 2 en la ecuación (2) se tiene que µu = 2.5 .Para sistemas estructurales en donde la falla se concentra en la planta baja, como las que se estudian en este trabajo, la ductilidad global y la del primer nivel se relacionan (Paulay y Priestley, 1992) con:
donde µu es la ductilidad en la planta baja y n es el número de niveles de la estructura. El término 2/3 que aparece en la ecuación (3) supone que el desplazamiento debido a la fluencia, en el primer nivel de una estructura, es aproximadamente dos terceras partes del desplazamiento total de la estructura (Paulay y Priestley, 1992). Como en este trabajo se calcularan las envolventes de resistencia de la planta baja de la estructura, será necesario calcular la ductilidad global conocida la ductilidad del primer entrepiso,
Una vez calculada la curva envolvente de resistencia del nivel 1, la ductilidad, µ1, se obtiene con el criterio de ductilidad equivalente (Priestley, Seible y Calvi, 1996), en donde la respuesta se idealiza según una curva bilineal equivalente (Figura 6b),
en donde du es el desplazamiento último equivalente y de es el desplazamiento idealizado de diseño.
Este criterio se basa en la determinación de una rigidez inicial secante al 75% de la carga última, la cual a su vez se define como la asociada al 15% de degradación de resistencia.
CURVA ENVOLVENTE IDEALIZADA
En la Figura 7 se muestra una curva típica del comportamiento histerético obtenida experimentalmente en un muro de mampostería sujeto a cargas laterales cíclicas y que relaciona la carga lateral medida con sus respectivos desplazamientos laterales. Las curvas envolventes de resistencia varían según sea el tipo de mampostería, tal que la mampostería confinada por dalas y castillos tiene un mejor comportamiento que la mampostería no confinada. El comportamiento de la mampostería por lo general se caracteriza por las deformaciones por cortante que son las que dominan. En una curva idealizada (Tomazevic y Klemenc, 1997a) que represente la envolvente de resistencia de un muro se identifican tres zonas principales: (a) Zona de comportamiento elástico, en la que se presenta agrietamiento horizontal en los castillos y agrietamiento mínimo en los elementos de mampostería, en esta zona se puede decir que se presenta un comportamiento prácticamente elástico asociado con una rigidez lateral efectiva, Ke ; (b) Zona de máxima resistencia, que está asociada a una carga lateral máxima, Vmax , y su respectivo desplazamiento, dVmáx. En esta zona la rigidez del muro, KVmáx, es menor que la rigidez anterior al primer agrietamiento, lo que le permite alcanzar la máxima resistencia; y (c) Zona de degradación de rigidez y de resistencia, definida por la carga lateral que produce el desplazamiento último o justo antes del colapso.
Con base en las tres zonas descritas se propone la curva envolvente del comportamiento idealizada mostrada en la Figura 8. La fuerza cortante resistente de diseño especificada en las NTCM-RCDF (2004) se establece como la fuerza cortante correspondiente al primer agrietamiento diagonal de la mampostería, Vagr :
donde vm* es la resistencia a compresión diagonal de la mampostería, σ es el esfuerzo de compresión y AT es el área transversal del muro, incluidos los castillos, sin transformar su área. La mampostería no se comporta elásticamente, ni siquiera para pequeñas deformaciones. Sin embargo, la rigidez elástica efectiva, Ke , se puede definir como la pendiente de la primer recta de la curva idealizada,
donde dagr es el desplazamiento asociado a la primera grieta significante en el muro. En muros sin refuerzo transversal la resistencia máxima a cortante se alcanza después de que se ha generalizado el agrietamiento inclinado. La resistencia en este punto será la suma de la resistencia al corte de los castillos y la contribución de la mampostería agrietada. En general, cuando se presenta el cortante máximo en el muro, la grieta inclinada ha penetrado los castillos. En Flores (1996) se propone que la fuerza cortante máxima se puede estimar como:
donde η = 0.30 es un factor de eficiencia y Vd representa la contribución del refuerzo vertical por acción de dovela. Dulacska, Flores y Alcocer (2001) proponen:
donde db es el diámetro de las varillas; fc' y fy son el esfuerzo a compresión del concreto y de fluencia del acero de refuerzo, respectivamente. Considerando que este valor es proporcional a la fuerza cortante de agrietamiento se tiene,
Sustituyendo la ecuación (10) en la (8) se puede escribir,
donde x vale 0.25 (Flores y Alcocer, 2001). Después del agrietamiento, el muro tiene una rigidez secante definida como la relación entre su resistencia lateral y su correspondiente desplazamiento. La rigidez asociada a la fuerza cortante máxima se puede definir como:
Sustituyendo el cortante asociado al agrietamiento obtenido con la ecuación (7) en la (12) se tiene,
El grado de degradación de la rigidez secante depende del daño, de la correlación del índice de daño y de la rigidez secante del muro. Tomazevic (1999) considera que la relación de rigideces varía entre 0.32 y 0.42. En este trabajo se considera que,
por lo que el desplazamiento asociado es:
El desplazamiento máximo se presenta para una fuerza cortante de,
tal que la relación de rigideces es:
y por tanto el desplazamiento máximo se presenta en:
La recta que une los puntos de desplazamiento máximo y cortante máximo, tiene una pendiente negativa con magnitud:
y el desplazamiento último (Figura 8) es:
es decir,
El punto en donde se presenta la resistencia última no corresponde al cortante de diseño, sino a un valor máximo experimental idealizado. Tomazevic (1999) encontró que el valor obtenido experimentalmente en más de 60 muros es:
y la rigidez secante asociada es:
para calcular el desplazamiento último, du , se considera que la relación entre rigideces tangentes para esta zona es:
El índice de degradación de rigidez lateral (Anónimo, 1999) se puede escribir en forma aproximada como,
donde R es la distorsión lateral del muro. Considerando como correcta esta ecuación, para una distorsión lateral R=0.002 , la cual controla el estado límite de resistencia (Astroza y Schmidt, 2004), se tiene:
que coincide con el valor considerado en este trabajo (ecuación 14). Y para una distorsión lateral de R = 0.0044, que controla el estado límite último (Astroza y Schmidt, 2004),
que es igual a la propuesta en este trabajo en la ecuación 24. De esta manera, el modelo de la curva envolvente de histéresis refleja matemáticamente los resultados experimentales obtenidos en México. En la Figura 9 se observa la curva envolvente propuesta, la relación entre los desplazamientos y la distorsión del muro, daño asociado, y su relación con los Estados Límite propuestos por Astroza y Schmidt (2004).
HIPÓTESIS DE CÁLCULO PARA LA CURVA ENVOLVENTE DE ENTREPISO
Se realizan las siguientes hipótesis para el cálculo de la curva envolvente:
(1) Las losas de piso trabajan como diafragmas horizontales rígidos, tal que los desplazamientos en cada uno de los muros del entrepiso i-ésimo, Di , debido a la fuerza cortante aplicada en el entrepiso, Vi , se distribuyen proporcionalmente a la rigidez lateral de cada muro, kji .
(2) Los desplazamientos debidos a los efectos por torsión son proporcionales a la rigidez del muro.
(3) Predomina el primer modo de vibrar. Se imponen desplazamientos horizontales que siguen la forma modal del primer modo (Figura 10b).
(4) La contribución individual de los muros para resistir la fuerza lateral del entrepiso dependen de los desplazamientos laterales tomados por el muro y de la forma de la curva envolvente idealizada de resistencia del muro. Los muros resisten los desplazamientos impuestos hasta agotar su capacidad de ductilidad individual, aunque pueden fallar por carga lateral, todavía son capaces de resistir carga vertical.
(5) En los muros con forma de "T", "L" o en forma de cruz, se considera que los patines están desligados verticalmente del alma.
(6) Las aberturas como puertas y ventanas mantienen un esquema uniforme en elevación.
PROCEDIMIENTO DE CÁLCULO
La metodología propuesta para calcular la curva envolvente de resistencia de entrepiso es iterativa y programable en una hoja de cálculo, el procedimiento se inicia con:
(1) Se calculan las constantes de cálculo del análisis como son: (a) el módulo de elasticidad de la mampostería y del concreto; (b) el módulo de rigidez al corte y compresión diagonal de la mampostería; (c) la masa de cada nivel considerando los entrepisos a ejes medios en la altura y (d) los esfuerzos a compresión actuantes en cada muro.
(2) Se calcula de la rigidez lateral efectiva de cada muro y la fuerza cortante resistente. En este trabajo, la rigidez lateral efectiva para cada muro se calculó con (UNIANDES, 2001):
donde h es la altura efectiva del muro, Em es el módulo de elasticidad de la mampostería, Gm es el módulo de rigidez al corte, A es el área efectiva al cortante, β es un factor que depende las condiciones de apoyo en los extremos del muro y que se considera igual a 3 para muros en voladizo y de 12 para muros doblemente empotrados, en este trabajo se consideró β=12.0 . El valor del factor correctivo, γ, es de 1.2 para muros con sección rectangular, aunque puede usarse un valor de γ=1.0 si los muros tienen castillos en sus extremos que actúan como patines; en este trabajo se consideró γ=1.2. El momento de inercia de la sección transformada, It, se puede calcular con:
donde Im es el momento de inercia de la mampostería confinada,
donde t es la espesor de la mampostería y L es la longitud del muro. La contribución del acero al momento de inercia, Ia , se puede obtener con:
y la contribución del concreto que forma los castillos, Ic , con:
donde As1 y As2 son el área de acero que forma los castillos, dsi y dci son las distancias del centroide del muro al centroide del acero y al centroide del castillo, respectivamente; bi y tc son el ancho y el espesor de los castillos, respectivamente. Las relaciones entre los módulos de elasticidad de la mampostería con el acero de refuerzo, Es , y con el concreto que forma los castillos, Ec, están definidas respectivamente por:
(3) Se forma la matriz de rigidez lateral de los muros que forman los entrepisos, la cual se puede aproximar mediante un modelo de cortante despreciando las deformaciones por flexión que presentan los muros.
(4) Para considerar los efectos por torsión (Chipol, 2001) se calculan: (a) la posición del centro de masas (CM); (b) el centro de rigidez (CR); (c) la excentricidad; (d) el momento debido a torsión; (e) la rigidez torsional y (f) los desplazamientos debidos a la torsión. El sistema estructural debe satisfacer la fórmula aproximada propuesta en el método simplificado de diseño (NTCM-RCDF, 2004) para calcular la excentricidad estática en el nivel i-ésimo, esi , que se toma como la distancia entre centroide de las áreas efectivas de los muros con respecto al centro de cortante. La excentricidad estática calculada debe ser,
Bi es la dimensión perpendicular del entrepiso medida en dirección del análisis (Figura 11).
(5) Se realiza la sumatoria de la fuerza cortante y de la rigidez de los muros en el sentido del análisis por nivel.
(6) Se calculan los desplazamientos, incluyendo los referentes a torsión.
(7) Se calcula la forma modal del primer modo de vibrar con cualquier metodología de dinámica estructural. En este trabajo se utilizó el método iterativo del vector inverso (Chopra, 1995).
(8) Se revisa el criterio de falla. Si la estructura aún no alcanza la falla se impone un desplazamiento arbitrario en el último nivel y los cálculos se repiten, regresando al paso (2).
(9) Se calcula la ductilidad y el factor de comportamiento sísmico.
(10) Finaliza el procedimiento.
El diagrama de flujo del procedimiento anteriormente descrito se muestra en forma esquemática en las Figuras 12 y 13.
CALIBRACIÓN DEL PROCEDIMIENTO
Con el fin de comprobar el procedimiento de cálculo de las curvas envolventes de resistencia lateral de piso, se calcularon las curvas envolventes de los modelos reportados en Tomazevic y Klemenc (1997b) usando la curva trilineal propuesta por dichos autores. Los modelos corresponden a una estructura de tres niveles de mampostería confinada y probados en una mesa vibradora. Al modelo de la estructura probada en la dirección longitudinal es M1 y en la dirección transversal M2. En la Figura 14 se comparan las curvas envolventes de resistencia del primer nivel calculadas analíticamente en este trabajo con las reportadas experimentalmente en dicha referencia. En las Figuras 14a y 14b se observa que los valores de la curva analítica que definen la rama ascendente, el cortante máximo, la fuerza cortante en la que se presenta el primer agrietamiento así como el desplazamiento que lo acompaña se compra bien con los valores obtenidos experimentalmente. Las curvas analíticas calculadas en este trabajo difieren a lo más en 5% con respecto a las analíticas reportadas en Tomazevic y Klemenc (1997b),atribuible a la diferencia de valores tomados como referencia en las curvas trilineales.
En la Tabla 1 se resumen los valores de ductilidad global y local, así como el del factor de comportamiento sísmico, calculados analíticamente y los obtenidos con los valores experimentales de Tomazevic y Klemenc (1997b). Además en la última columna se listan los valores del factor de comportamiento sísmico, Q ,calculado analíticamente por dichos autores. Se ve que tanto la ductilidad como los del factor de comportamiento sísmico predichos por el modelo son adecuados.
COMPARACIÓN CON RESULTADOS EXPERIMENTALES REALIZADOS EN MÉXICO
Para comparar el procedimiento y el modelo de la curva envolvente propuesta, se compararon las curvas envolventes de resistencia de piso obtenidas analíticamente con las experimentales en una serie de modelos ensayados en la mesa vibradora del Instituto de Ingeniería de la UNAM, II-UNAM (Vázquez, 2005; Barragán, 2005; Arias, 2005) y con la obtenida en el ensayo con carga cíclica reversible del modelo 3D (Sánchez, 1998). Los datos referentes a dimensiones y propiedades de los materiales empleados en los modelos se especifican en las referencias antes citadas.
El primer modelo ensayado en la mesa vibradora representa un edificio de un nivel (Vázquez, 2005). La distribución en planta de los muros del modelo inicial, M1SRSC se muestra en la Figura 15a y su vista en elevación en la Figura 15b. Debido a las características mecánicas de la mampostería usada y a la configuración estructural, el modelo resultó con rigidez y resistencia elevada, por lo que debió ser sometido a diversos ciclos de carga. En la Figura 16a se muestra la comparación de los valores experimentales con la curva envolvente calculada con el procedimiento propuesto en este trabajo. Se comprueba que el modelo se encontraba lejos de la zona del límite elástico durante la prueba.
Para continuar con las pruebas experimentales se modificó el modelo M1SRSC en su distribución lateral. La primera modificación que se realizó fue retirar dos de los muros centrales, y se llamó M1-M (Figura 15c). En la Figura 16b se muestran las curvas envolventes obtenidas analítica y experimentalmente. La prueba dinámica se suspendió cuando se presentó el primer agrietamiento, punto que coincide con el obtenido con la curva analítica.
La segunda modificación al modelo fue separar verticalmente los muros centrales (Figura 15d); a este modelo se le llamó M1-A. La separación de los muros centrales modificó ligeramente la curva experimental; la prueba se suspendió cuando ocurrió el primer agrietamiento, este punto lo predice adecuadamente la curva envolvente (Figura 17a).
La tercera modificación fue colocar en el modelo M1-A un peso de 15.7 kN (1600 kgf) y recibió el nombre de M1-B. La curva envolvente y los valores experimentales se muestran en la Figura 17b. Se puede ver que aparentemente la carga máxima en la que se suspendió la prueba resultó mayor que el valor máximo que cubre la curva envolvente; esto se explica más adelante.
La cuarta modificación al modelo consistió en colocar sobre el modelo M1-A un peso de 20.6 kN (2100 kgf) y recibió en nombre de M1-C. En la Figura 18a se muestra la curva envolvente y los valores experimentales obtenidos en este modelo. Debido a que el modelo M1-C presentaba agrietamiento importante en el momento de la prueba la curva envolvente y los valores experimentales se separan en la rama ascendente, lo cual no pasa al comparar en la misma figura los valores del modelo M1-M. También se observa que la carga máxima experimental es mucho mayor que la que se predice en forma analítica usando un módulo de cortante promedio. Sin embargo, los valores del módulo de cortante de la mampostería utilizada en el modelo presentan un coeficiente de variación (CV) de 0.39. Si se calculan las curvas envolventes con el valor medio más el CV y el valor medio menos el CV, se observa la influencia de este parámetro sobre la forma y magnitud de las cargas resistentes(Figura 18a). Para este caso la curva calculada con el valor de esfuerzo promedio más el CV, envuelve a los datos experimentales. El esfuerzo de fluencia del acero de refuerzo utilizado en los castillos durante la prueba también presentó variaciones importantes, pero como se observa en la Figura 18b, no tiene una influencia importante en la curva envolvente. En este trabajo se utilizó un valor de 440 MPa que es el valor medio del acero utilizado en la prueba.
En la segunda serie de estructuras de mampostería probadas en la mesa vibradora del II-UNAM, se modeló una estructura de dos niveles (Barragán, 2005) con la misma estructuración y características de los materiales empleados en el modelo de un nivel, M1-M, descrito en la Figura 15c. Dicho modelo recibió el nombre de M2SRCC y se muestra en la Figura 19a. La curva envolvente de resistencia del primer nivel junto con los resultados experimentales se presentan en la Figura 19b. En este caso la curva envolvente se calculó considerando:
donde fm'(kgf/cm2) es el esfuerzo a compresión de la mampostería.
No fue posible analizar el desplazamiento que se presenta en la carga lateral última ya que las pruebas experimentales se suspendieron cuando se aplicó una fuerza del 89% de la resistencia máxima. Sin embargo, se corrobora que el modelo y el procedimiento para calcular la curva envolvente de respuesta predicen adecuadamente la carga del primer agrietamiento diagonal, la carga máxima y sus respectivos desplazamientos.
Como se observa en la Figura 20 el tercer modelo de mampostería ensayado en la mesa vibradora del II-UNAM, representa una estructura de tres niveles (Arias, 2005) y tenía similares características mecánicas de los materiales empleados en los anteriores modelos: M1SRSC y M2SRCC. El modelo fue ensayado con dos configuraciones diferentes: el M3SRCC con una distribución de muros similar al modelo M1SRSC (Figura 20a) y el modelo M3-M sin los muros centrales (Figura 20b).
Las curvas envolventes de resistencia del primer nivel se calcularon con los valores considerados anteriormente. Las pruebas experimentales del modelo M3SRCC se suspendieron antes de que se presentara degradación de resistencia, por lo que los valores experimentales sólo son comparables con los analíticos en la rama ascendente (Figura 21a). Se encontró que los valores analíticos calculados en este trabajo se ajustan mejor con los experimentales que los calculados con la metodología de la columna ancha (Figura 1b). Los datos experimentales obtenidos con el modelo M3-M no son comparables con los resultados analíticos debido a que cuando se ensayó este modelo se encontraba agrietado debido a ensayos anteriores (Figura 21b).
El cuarto modelo es el ensayado en el Laboratorio de Estructuras Grandes del CENAPRED, con carga cíclica reversible y una configuración en planta ideal y simétrica (Sánchez, 1998). Este modelo es llamado 3D (Figura 22a) y consistió de dos niveles a escala 1:1 con mampostería confinada. La curva envolvente de resistencia del primer nivel y los resultados experimentales se muestran en la Figura 22b. Se observa que la curva envolvente predice muy bien la carga lateral máxima y su desplazamiento asociado; además del cortante de agrietamiento, aunque no el desplazamiento en el que se presenta. En este caso no fue posible verificar el desplazamiento a la carga última ya que el modelo no se llevó al colapso debido a que durante el trabajo experimental se planeaba repararlo.
Los resultados descritos en los párrafos precedentes corroboran que son adecuados los valores propuestos para la curva envolvente idealizada y la validez del procedimiento propuesto para el cálculo de la curva envolvente de resistencia lateral de entrepiso.
En la Tabla 2 se registran los valores de las ductilidades de entrepiso y global, así como del factor de comportamiento sísmico, calculados con la metodología desarrollada en este trabajo. En todos los casos el factor de comportamiento sísmico es mayor que el prescrito por la NTCM-RCDF (2004).
APLICACIÓN A UNA ESTRUCTURA HIPOTÉTICA
En este apartado se aplicará el procedimiento propuesto para calcular la curva envolvente de resistencia lateral de entrepiso de un edificio hipotético de cuatro niveles de mampostería confinada que satisface todos los requisitos para que sea aplicable el Método Simplificado de diseño. Tiene las siguientes características: (1) La excentricidad estática en la dirección de análisis es la máxima permitida para aplicar el Método Simplificado. (2) El cortante resistente en la planta baja, en cada sentido del análisis, está lo más próximo posible al cortante actuante de diseño, y (3) las acciones serán las máximas esperadas según la NTCM-RCDF (2004).
En la Figura 23 se observa la distribución en planta de los muros y la elevación en la dirección X del edificio por analizar. En la dirección X el edificio presenta una excentricidad de 0.83m que es el 97% del ancho en planta del edificio, siendo la excentricidad máxima permitida en esta dirección 0.85 m. La altura de los entrepisos es 2.50m y se considera que el edificio se ubica en la zona sísmica III del Valle de México. El sistema de piso está formado por una losa de 0.10m de espesor y los muros tienen un espesor de 0.12m. El diseño se realizó tal que la mayoría de los muros de planta baja están esforzados al límite del cortante último actuante. En la tabla 3 se muestran las propiedades de los materiales utilizadas en este ejemplo.
Para obtener el esfuerzo a compresión actuante en cada muro se empleó el criterio de áreas tributarias, considerando un peso de muro de 580 kgf/m y en la azotea una carga de servicio de 620 kgf/m2 y en los niveles inferiores de 633 kgf/m2. Los valores de las masa concentradas en cada entrepiso se muestran en la Tabla 4.
En el inicio del proceso del cálculo, la matriz de rigidez en la dirección de análisis X es:
En la Figura 24 se muestran las curvas envolventes de resistencia de la planta baja del edificio en consideración. Para su elaboración se utilizó:
Las curvas envolventes, en ambos sentidos, predicen un cortante de agrietamiento mayor que el cortante actuante de diseño: VuX = VuY = 454.10kN, por lo que se espera que los muros permanezcan en el intervalo elástico. Además, las deformaciones por torsión hacen que en la dirección Y se tenga una resistencia máxima ligeramente mayor que en X, a pesar de que la densidad de muros es menor. En la Figura 25 se muestra la evolución de la forma modal del primer modo de vibrar y el mecanismo de falla del llamado piso débil. Esto concuerda con lo observado experimentalmente y también con lo ocurrido después de un sismo en este tipo de estructuras.
La metodología de incrementar los desplazamientos y no las cargas laterales, resultó estable, evitándose problemas de inestabilidad numérica al invertir la matriz de rigidez en cada ciclo del proceso. Una ventaja que presenta la metodología propuesta es que se puede determinar el estado de los muros que pertenecen a un entrepiso en particular, tal y como se muestra en la Figura 26. En la Figura 26a se muestran las curvas envolventes del entrepiso 1 y 2; en la Figura 26b las curvas envolventes del entrepiso 1 y 3, y finalmente en la Figura 26c las curvas del entrepiso del nivel 1 y 4. Se puede ver que el cortante de agrietamiento disminuye en las envolventes de los pisos superiores, debido a que su capacidad de resistencia lateral es menor que la de los muros del primer nivel por tener menor esfuerzo a compresión. Cuando se cumple alguno de los criterios de falla en el nivel 1, el segundo nivel presenta una envolvente que se aproxima a la carga lateral máxima. Esto es congruente con la evidencia experimental en edificios con baja altura (Tomazevic y Klemenc, 1997).
Otro beneficio de este tipo de análisis es que permite determinar el estado de un muro en particular, para cierta fuerza cortante lateral. Como ejemplo, en las Figuras 27a a la 27c (27b), se muestran las curvas envolventes de resistencia lateral de los muros 1, 2 y 10 de la planta baja del edificio (Figura 23a), además, se indica el punto en el que el nivel 1 de la estructura alcanza el cortante de diseño. En la Figura 27d se muestran las envolventes de resistencia lateral del muro No. 1 ubicado en los diferentes niveles de la estructura. Así se pueden determinar el estado límite de comportamiento que tiene cada muro en particular y el estado de daño asociado.
En la Tabla 5 se muestran los valores de ductilidad y del factor de comportamiento sísmico calculados para este ejemplo. Se observa que el Q en ambas direcciones de análisis se mantiene dentro del límite especificado de 2 (NTCM-RCDF, 2004). Si se considera un edificio de mampostería con 5 niveles y con altura de entrepiso de 2.50m, como límite de aplicación del Método Simplificado de diseño, la ductilidad de agrietamiento, µ1, debe ser 6 para que el edificio tenga un Q = 2. En la Tabal 3 se puede ver que los valores de dicha ductilidad se acercan a dicho valor.
CONCLUSIONES
Para evaluar el comportamiento sísmico de una estructura de mampostería es necesario realizar pruebas experimentales y desarrollar procedimientos analíticos que reproduzcan los resultados experimentales para su aplicación en la práctica ingenieril. En este trabajo se propone un método para el cálculo de la curva envolvente de resistencia lateral de entrepiso de estructuras de mampostería confinada con baja o mediana altura, que cumplen con los requisitos básicos para que sea aplicable el Método Simplificado de diseño (NTCM-RCDF, 2004). El procedimiento consiste en incrementar desplazamientos laterales, suponiendo un comportamiento fundamental en el primer modo de vibrar de la estructura. La curva envolvente idealizada se basa en resultados experimentales disponibles en México. Las curvas envolventes analíticas se compararon con las obtenidas experimentalmente en modelos ensayados bajo carga cíclica reversible y en la mesa vibradora del II-UNAM. Los resultados indican una adecuada predicción de la rigidez elástica, de la carga de resistencia asociada al primer agrietamiento y al estado de resistencia máxima; en algunos casos también fue posible predecir el desplazamiento asociado a cada punto anterior. El método propuesto permite calcular la ductilidad global y la de entrepiso, así como del factor de comportamiento sísmico referido en la Norma Técnica Complementaria al Reglamento de Diseño del Distrito Federal. La metodología se aplicó a un edificio de mampostería confinada con cuatro niveles cuyos muros se encuentran esforzados al límite máximo permitido por la reglamentación vigente en el Distrito Federal. Los resultados muestran un adecuado comportamiento sísmico de este tipo de estructuras, cuantificado por el factor del comportamiento sísmico y el valor de la ductilidad calculada.
AGRADECIMIENTOS
El segundo autor desea agradecer al Instituto de Ingeniería de la UNAM y en particular al Dr. Sergio Alcocer M. por su apoyo en la realización de este trabajo mediante la estancia en el laboratorio de la Mesa Vibradora de la UNAM. De igual manera se desea agradecer a la Universidad Michoacana de San Nicolás de Hidalgo, UMSNH, y al Ing. Iván Cornejo Vásquez por el apoyo otorgado para finalizar el trabajo.
REFERENCIAS
Anónimo (1999), "Edificaciones de Mampostería para Vivienda", Fundación ICA AC, México, DF, ISBN 968-7508-65-5. [ Links ]
Alcocer S M, J Arias y L Flores (2004), "Some developments on performance-based seismic design of masonry structures", Performance Based Seismic Design Concepts and Implementation, Proceeding of an International Workshop, Bled, Slovenia, June 28, PEER 2004/05. [ Links ]
Arias, A J (2005), "Ensayos en mesa vibradora de un modelo a escala 1:2 de edificio de mampostería confinada de tres niveles", Tesis de Maestría, Universidad Nacional Autónoma de México. [ Links ]
Astroza, M y A Schmidt (2004), "Capacidad de deformación de Muros de Albañilería Confinada para Distintos Niveles de Desempeño", Revista de Ingeniería Sísmica, No. 70, pp. 59-75. [ Links ]
Barragán, T R (2005), "Ensaye de una vivienda a escala de dos niveles de mampostería confinada", Tesis de Maestría, Universidad Nacional Autónoma de México. [ Links ]
Bazán, E (1980), "Muros de mampostería ante Cargas Laterales Estudios Analíticos", Tesis Doctoral, Facultad de Ingeniería, UNAM. [ Links ]
Beyer, K, A Dazio y M Priestley (2008), "Inelastic Wide-column Models for U-Shaped Reinforced Concrete Walls", Journal of Earthquake Engineering, Taylor and Francis, pp. 1-33. [ Links ]
Chipol, A (2001), "Estudio de la respuesta sísmica de modelos tridimensionales de edificios torsionalmente acoplados", Tesis de Maestría, Programa de Maestría y Doctorado en Ingeniería, UNAM. [ Links ]
Chopra, A K (1995), Dynamics of Structures, Theory and Applications to Earthquake Engineering, Prentice Hall, Inc., ISBN 0-13-855214-2. [ Links ]
Chrysostomou, C Z, P Gergely y J F Abel (1992), "NonlinearSeismic Response of InfilledFrames", Memorias X Congreso Mundial de Ingeniería Sísmica, Madrid, España, pp. 4435-4437. [ Links ]
Coral, M (2004), "Revisión de algunas de las hipótesis del método simplificado de análisis sísmico para muros de mampostería confinada", Tesis de Maestría, Programa de Maestría y Doctorado en Ingeniería, UNAM. [ Links ]
Crisafulli, F (1997), "Seismic Behaviour of Reinforced Concrete Structures with Masonry Infills", Tesis de Doctorado, University of Canterbury. [ Links ]
Fischinger M, D Berg, T Isakovic, M Tomazevic y R Zarnic (2004), "Performance based assessment from general Methodologies to specific implementations",Performance Based Seismic Design Concepts and Implementation, Proceeding of an International Workshop, Bled, Slovenia, June 28, PEER 2004/05. [ Links ]
Flores, L (1995), "Estudio analítico de estructuras de mampostería confinada", Tesis de Licenciatura, Facultad de Ingeniería de la UNAM. [ Links ]
Flores, L y S M Alcocer (2001), "Estudio Analítico de Estructuras de Mampostería confinada", Informe de tesis de maestría, Centro Nacional de Prevención de Desastres, CENAPRED. [ Links ]
Holmes, M (1961), "Steel Frames with brickwork and concrete Infills", Proceeding of Instn.Civil Engineering, Vol. 19, pp. 473-478. [ Links ]
Ishibashi, K y H Kastumata (1994), "A Study on Nonlinear Finite Element Analysis of Confined Masonry Walls", Cuaderno de Investigación No. 15, CENAPRED. [ Links ]
Jean, R y F Cesín (2000), "Recomendaciones para el Diseño y Construcción de Estructuras de Mampostería", Primer Simposio Nacional de Estructuras de Mampostería para la Vivienda, México. [ Links ]
Kadir, M R (1974), "The Structural Behaviour of Masonry Infill Panel in Framed Structures", Tesis de Doctorado, Universidad de Edimburgo, Gran Bretaña. [ Links ]
Madan, A, A M Reinhorn, J B Mander y R E Valles (1997), "Modeling of Masonry Infill Panels for Structural Analysis", Journal of Structural Engineering, ASCE, pp. 1295-1302. [ Links ]
Meli, R (1979), "Comportamiento Sísmico de Muros de Mampostería", Reporte No. 352, Serie del Instituto de Ingeniería, Universidad Nacional Autónoma de México. [ Links ]
Moroni, M, M Astroza y S Tavonatti (1994), "Nolinear models for shear failure in confined masonry walls", The Masonry Soc. J., No. 12, Vol. 2, pp. 72-78. [ Links ]
Mosalam, K M, P Gregely, R N White y G Ayala (1996), "Effect of Geometrical configuration on the Seismic Response of Infilled Frames", Memoriasde la XI Conferencia Mundial de Ingeniería Sísmica, Acapulco, México, CDROM, artículo 1675. [ Links ]
Newmark, N M y W Y Hall (1982), Earthquake Spectra and Design, Earthquake Engineering Research Institute. [ Links ]
NTCM-RCDF (2004) "Normas Técnicas Complementarias para diseño y construcción de Estructuras de Mampostería, Reglamento de Construcciones del Distrito Federal", Gaceta Oficial del Distrito Federal, marzo. [ Links ]
NTCS-RCDF (2004) "Normas Técnicas Complementarias para Diseño por Sismo, Reglamento de Construcciones del Distrito Federal", Gaceta Oficial del Distrito Federal, Tomo II, No. 103-BIS, octubre. [ Links ]
Paulay, T, y M J N Priestley (1992), Seismic Design of Reinforced Concrete and Masonry Buildings, John Wiley and Sons, Inc. ISBN 0-471-54915. [ Links ]
Priestley, M J N, F Seible y G M Calvi (1996), Seismic Design and Retrofit of Bridges, John Wiley and Sons, Inc. ISBN 0-471-57998-X. [ Links ]
Sánchez, T A (1998), "Comportamiento de estructuras de mampostería confinada sujetas a cargas laterales", Tesis de Maestría, DEPFI-UNAM, Universidad Nacional Autónoma de México. [ Links ]
Sánchez, S, R Arroyo y S Jerez (2010), "Modelo de un grado de libertad para evaluar la curva carga lateral-distorsión en muros de mampostería confinada", Revista de Ingeniería Sísmica, No. 83, pp. 25-42. [ Links ]
Siete, B J (2011), "Fragilidad basada en desplazamiento para edificaciones de mampostería confinada", Tesis de Maestría, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán. [ Links ]
Tena, A, J Cano, A López y M A Pérez (2010), "Revisión y actualización del Método Simplificado de Análisis de Estructuras de mampostería de los Reglamentos de diseño sísmico de México", Revista de Ingeniería Sísmica, No. 82, pp. 1-34. [ Links ]
Tomazevic, M (1999), Earthquake-Resistant Design of Masonry Buildings, Series on Innovation in Structures and Construction, Vol. 1. Imperial College Press, ISBN 1-86094-066-8. [ Links ]
Tomazevic, M y I Klemenc (1997a), "Verification of Seismic Resistance of Confined Masonry Buildings", Earthquake Engineering and Structural Dynamics, Vol. 26, pp. 1073-1088. [ Links ]
Tomazevic, M y I Klemenc (1997b), "Seismic Behaviour of Confined Masonry Walls", Earthquake Engineering and Structural Dynamics, Vol. 26, pp. 1059-1071. [ Links ]
Tomazevic, M, V Bosiljkov y P Weiss (2004), "Structural Behavior Factor for masonry structures", 13th World Conference on Earthquake Engineering, Vancouver, B. C., Canada, August, paper No. 2642. [ Links ]
UNIANDES (2001), "Comportamiento sísmico de muros de mampostería confinada", Asociación Colombiana de Ingeniería Sísmica. [ Links ]
Vázquez, V A (2005), "Ensaye experimental de viviendas de mampostería confinada de un piso mediante el ensaye en mesa vibradora", Tesis de Maestría, Universidad Nacional Autónoma de México. [ Links ]
Valles, R E, A M Reinhorn, S K Kunnath, C Li y AMadan (1996), "IDARC 2D, version 4", National Center for Earthquake Engineering Research, State University of New York at Buffalo. [ Links ]
Zúñiga, O y A Terán (2008), "Evaluación basada en desplazamientos de edificaciones de mampostería confinada", Revista de Ingeniería Sísmica, No. 79, pp. 25-48. [ Links ]
* Se aceptarán comentarios y/o discusiones hasta cinco meses después de su publicación.