INTRODUCCIÓN
La tendencia de diseño sismorresistente de la mayoría de los reglamentos internacionales permite que las estructuras ingresen en el intervalo inelástico durante sismos intensos, lo que proporciona dos ventajas fundamentales: un mecanismo de disipación de energía y una reducción significativa de la resistencia necesaria. Esta capacidad de las estructuras de ingresar en el intervalo de comportamiento plástico de manera estable se sustenta en la respuesta que los elementos estructurales individualmente son capaces de desarrollar, pero también en la interacción del sistema como un conjunto (incluyendo la cimentación).
En general, el estudio de la respuesta dinámica de las estructuras es un fenómeno complejo que depende de muchas variables: el comportamiento individual de los elementos que conforman la estructura, la respuesta inelástica global de la estructura considerando las posibles redistribuciones de carga en el intervalo inelástico, las posibles fuentes de excitación y sus características particulares, la respuesta dinámica del suelo y su interacción con la cimentación, etc. En este sentido, es común que el ingeniero tome decisiones durante el proceso de diseño que pretenden simplificar la complejidad y reducir el tiempo de ejecución y post-proceso de los análisis.
Entre ellas, una de las decisiones más recurrentes consiste en suponer que el sistema suelo- cimentación es lo suficientemente rígido como para despreciar sus efectos en los análisis. Sin embargo, esa simplificación no es correcta en todos los casos, puesto que en menor o mayor magnitud la base de la edificación tiene desplazamientos relativos durante un sismo intenso. Así, una condición flexible entre el sistema suelo - cimentación - estructura modificaría las propiedades dinámicas del sistema, lo que podría implicar una redistribución de cargas imprevista entre los elementos sismo-resistentes. Esa variación de la respuesta dinámica del suelo y de la estructura provocadas por la flexibilidad del sistema suelo - cimentación en un movimiento sísmico es conocida como interacción dinámica suelo - estructura (IDSE).
La mayoría de los reglamentos especializados consideran los efectos de la IDSE, a partir de una proporción que se establece en función de criterios para calcular el periodo de la estructura y amortiguamiento modificados por la flexibilidad de la base. Sin embargo, otras variables de la respuesta de la estructura, como el efecto P-Δ en las columnas, puede incrementarse con el aumento de flexibilidad de la cimentación, especialmente en edificios flexibles. De hecho, algunos estudios como Halabian y Kabiri (2004) y Barcenas y Esteva (2006) han demostrado variaciones importantes en los desplazamientos máximos que desarrollan los modelos de estructuras de concreto armado cuando se incluye la IDSE, así como variaciones en las demandas en los elementos estructurales.
Así, en esta investigación se evalúa la influencia de la interacción dinámica suelo - estructura en la respuesta inelástica en marcos de acero, que forman parte de edificios regulares de 8, 12 y 16 pisos, que fueron diseñados para las mayores demandas de aceleración propuestas para la Ciudad de México. Adicionalmente, en el estudio se ha incluido como variable:
La rigidez lateral de los modelos al considerar marcos con dos crujías contraventeadas, marcos con una crujía contraventeada y marcos sin contravientos.
La flexibilidad del sistema suelo-cimentación mediante la consideración de modelos con cimentación rígida, cimentación con pilotes y modelos desplantados en un cajón de cimentación.
CRITERIOS NORMATIVOS
Las Normas Técnicas Complementarias para el Diseño por Sismo (NTCS-04) del Reglamento de Construcciones para el Distrito Federal (RCDF-04) y el Manual de Obras Civiles de la Comisión Federal de Electricidad (MOCS-08) proponen un procedimiento para incorporar los efectos de la IDSE mediante una modificación del periodo de la estructura y su amortiguamiento; de manera que la flexibilización del sistema conlleva implícitamente la posibilidad de que la demanda de aceleración espectral sea modificada. Esta metodología establece que una vez que se determine el cortante basal correspondiente a la aceleración espectral para el periodo y el amortiguamiento modificados, el resto de las respuestas se modifiquen en la misma proporción.
En los reglamentos internacionales vigentes, las modificaciones asociadas con la IDSE más frecuentemente consideradas son: a) el alargamiento del periodo estructural y b) la modificación del amortiguamiento, debido a que son las más conocidas (Wolf 1985). Pese a esto, otras variables también deberían ser modificadas, sobre todo las que pretenden proporcionar al diseñador una estimación realista de la respuesta inelástica; por ejemplo, la ductilidad μ que desarrollará la estructura.
Específicamente, las NTCS-04 y el MOCS-08 establecen un procedimiento para determinar la modificación de la ductilidad debida a la flexibilización provocada por los efectos de la IDSE. Esta modificación, propuesta por Rosenbleth y Resendiz (1988) y reformulada posteriormente por Pérez-Rocha y Avilés (2003), se basa en que la IDSE varía la rigidez inicial Kelast (Figura 1) debido a la flexibilidad del sistema suelo - cimentación; esto ocasiona que la deformación de fluencia Uy tenga una magnitud más grande de la que se obtendría del análisis con base rígida y, por consiguiente, también se modifica la deformación última Uu. Entonces, si las deformaciones que definen la ductilidad μ se alteran, la ductilidad que puede desarrollar el sistema se reduce (por el incremento de Uy ), pese a que la estructura tenga la misma capacidad de deformación inelástica.
En ese sentido, la deformación de los análisis con base rígida están asociados únicamente a la deformación de la estructura U (Figura 1b); mientras que el desplazamiento de los análisis con base flexible depende del deformación de la estructura U y de un desplazamiento adicional provocado por la flexibilidad de la cimentación Uc que resulta de la suma del desplazamiento del giro del sistema ϕh más el movimiento de traslación de la cimentación UCIMENT.
De esta manera, la reducción teórica de la ductilidad puede determinarse con la ecuación 1, que está incluida en el Apéndice A de las NTCS-04, la cual es función de la relación entre el periodo de la estructura sobre base rígida Te y el periodo considerando la base flexible .
La ecuación anterior se desarrolló a partir de estudios del comportamiento de un oscilador de un grado de libertad (Rosenbleth y Resendiz 1988; Pérez-Rocha y Avilés 2003), que es una simplificación usada en investigaciones similares donde comúnmente la estructura se modela con un grado de libertad y el sistema de suelo - cimentación se modela con dos o tres grados de libertad (por ejemplo Ghannad y Jahankhah 2007 y Mashuli y Ghannad 2009). En contraparte, otros estudios como Ganjavi y Hao (2011) han demostrado que las demandas de resistencia y ductilidad en análisis de estructuras con múltiples grados de libertad pueden ser muy diferentes de las que se obtienen de análisis con un grado de libertad equivalente, dependiendo del número de pisos, la flexibilidad del suelo y el grado de inelasticidad. Por lo anterior, en esta investigación se pretende ampliar el conocimiento sobre el comportamiento inelástico de estructuras flexibles con múltiples grados de libertad con distintas condiciones de cimentación y evaluar los criterios de predicción de los reglamentos vigentes.
DESCRIPCIÓN DE LOS MODELOS
Los edificios estudiados se diseñaron siguiendo estrictamente los requisitos propuestos en el RCDF-04 para que los marcos de acero (trabes y columnas) en cada entrepiso sean capaces de resistir, sin considerar la contribución del sistema de contravientos, cuando menos el 50 por ciento de la fuerza sísmica actuante. En Tapia y Tena (2011) y Dávalos (2013) es posible encontrar más información de los modelos diseñados.
En la Figura 2 se muestran las dimensiones en planta y elevación de los tres tipos de marcos considerados en este estudio: marco central sin contravientos, marco intermedio con una crujía contraventeada y marco perimetral con dos crujías contraventeadas.
Los modelos de los edificios se supusieron localizados en la zona IIIb de la subzonificación sísmica del Distrito Federal (NTCS-04) para asociarlos a las mayores demandas de aceleración propuestas en el reglamento para suelo blando con poca rigidez rotacional. Las secciones transversales finales de los análisis que satisfacen estrictamente las restricciones de resistencia y rigidez del RCDF-04 se muestran en la Tabla 1; mientras que, en la Tabla 2 se presenta la información dinámica de los análisis elásticos realizados.
Modelo | Altura del edificio, HT (m) | Periodo fundamental (s) | Frecuencia circular (rad/s) | Ductilidad global μ | |
---|---|---|---|---|---|
Curva de capacidad | Curva bilineal idealizada | ||||
Edificio de 8 pisos | 28.0 | 0.40 | 15.71 | 2.64 | 1.97 |
Edificio de 12 pisos | 42.0 | 0.50 | 12.57 | 2.10 | 1.89 |
Edificio de 16 pisos | 56.0 | 0.63 | 9.97 | 1.99 | 1.86 |
Adicionalmente, se realizaron análisis no lineales ante carga estática monótona creciente en el programa Drain2DX (Prakash et al 1992) para verificar que la respuesta inelástica fuera congruente con las premisas de diseño: mecanismo de colapso columna fuerte - viga débil- contraviento más débil con una adecuada distribución en la altura de la respuesta inelástica. Para ejemplificar esto, en la Figura 3 se muestran las instantáneas de fluencia normalizadas por intervalos en código de colores de la secuencia de formación del mecanismo de colapso del marco intermedio y perimetral del modelo de 8 niveles. En la Figura 3 también se reporta la distorsión global última (
Es valioso hacer notar que la respuesta inelástica de los edificios tiene una dependencia de la altura de los modelos, pese a que fueron diseñados con la misma filosofía y condiciones de diseño (Tapia y Tena 2008). Esto es, los modelos más altos son menos dúctiles y desarrollan un mecanismo de colapso con una menor distribución de la respuesta inelástica en la altura. Con el propósito de hacer notar esta tendencia, en la Tabla 2 se muestra la ductilidad global obtenida de la curva de capacidad (cortante basal contra distorsión en la azotea) y la ductilidad global teórica que se obtuvo de curvas bilineales idealizadas. Es posible consultar otros detalles sobre la definición de estas curvas y una amplia discusión sobre la dependencia de la altura de los edificios y la ductilidad que pueden desarrollar en Tapia y Tena (2011).
En esta investigación, el desplazamiento de fluencia Uy se define como el punto de intersección entre la recta con pendiente igual a la rigidez elástica del sistema Kelast y la línea horizontal que define el cortante último. Es decir, los resultados que se discutirán a continuación están basados en la curva bilineal idealizada que se obtiene a partir de la curva de capacidad. Esta decisión pretende ajustar los resultados con el modelo elasto-plástico perfecto utilizado por Pérez-Rocha y Avilés (2003) para evaluar el criterio de las NTCS-04.
Características del suelo de desplante
En consistencia con las suposiciones de diseño, los perfiles estratigráficos para modelar la respuesta del suelo se obtuvieron con las propiedades reportadas en un Estudio de Mecánica de Suelos (ESVM 1986) de la Colonia Roma Sur, Delegación Cuauhtémoc (coordenadas: 19°24 N y 99°10 W) de la zona IIIb. El estudio forma parte de un proyecto llamado Estudios sobre sismicidad en el Valle de México, donde se realizaron pruebas in situ para obtener las velocidades de propagación de ondas de cortante por dos métodos de campo: pozo abajo (Downhole) y por sonda suspendida.
Las propiedades del suelo se definieron con base en el perfil de velocidades de propagación de ondas mostrado en la Figura 4 correspondiente a la Zona del Lago. Como se aprecia, las propiedades del suelo (Tabla 3) no presentan variaciones importantes hasta los 40 metros de profundidad, por lo que en este estudio, el suelo se idealizó como un sólo estrato.
Propiedad | Magnitud | |
---|---|---|
Estudio de Mecánica de Suelos de 1988 | Proyección al 2013 | |
Periodo del suelo estimado (Ts) | 2.34 s | 1.85 s |
Cohesión, c | 5.0 t/m2 | |
Velocidad de propagación de ondas de corte, Vs | 65.0 m/s | 77.9 m/s |
Peso volumétrico, γ | 1.30 t/m3 | |
Módulo de Poisson, v | 0.49 | |
Amortiguamiento, ξ | 0.05 |
Ahora, es conocido que el bombeo de agua en la Ciudad de México ha generado una pérdida en la presión de agua del subsuelo, induciendo a una consolidación acelerada. Esto ha provocado una reducción gradual del espesor de los depósitos y, por consiguiente, alteraciones en las propiedades dinámicas del suelo. Por esta razón, fue necesario adecuar las propiedades del suelo obtenidas en ESVM (1988) por causa del hundimiento regional del subsuelo del Valle de México con el propósito de obtener propiedades aproximadas modificadas para las condiciones cuando se desarrolló esta investigación (2013). Esta adecuación se realizó mediante el procedimiento propuesto por Aguilar et al. (2003).
Aguilar et al. (2003) estudiaron la evolución del hundimiento regional y su relación con la modificación de las propiedades del suelo, a partir de datos periódicos de bancos de nivel en la Zona del Lago y la Zona de Transición entre 1983 y 1998. Con esta información, relacionaron la velocidad de hundimiento con el espesor del suelo compresible, lo que permitió estimar la variación de esa velocidad en función del tiempo. Asimismo, Aguilar et al. (2003) establecieron el cambio del periodo fundamental del suelo Ts debido a este fenómeno, a partir de la relación entre la velocidad de hundimiento, espesor comprensible de los estratos Hs y la velocidad de propagación de las ondas de corte Vs .
De esta manera, el periodo dominante de vibración del suelo inicial se estima mediante la ecuación 2. Luego, a partir de gráficas que son función del tiempo transcurrido, es posible estimar el periodo dominante del suelo en el tiempo actual Ts 2013 (en este caso al 2013), que se utiliza para obtener la velocidad de propagación al 2013 mediante la ecuación 3 (Aguilar et al., 2003).
En las ecuaciones anteriores, Hs 1988 es el espesor del estrato cuando se realizó el Estudio de Mecánica de Suelos que se consideró igual a 38 m (Figura 4); Vs 1988 es la velocidad del estrato igual a 65 m/s (Tabla 3) ESVM (1988); Ts 2013 es el periodo dominante del suelo que se estimó mediante el procedimiento descrito en Aguilar et al., (2003) considerando que habían transcurrido 25 años y Vs 2013 es la velocidad estimada al 2013 (Tabla 3).
Como consecuencia del hundimiento regional (cuya proyección para el sitio de estudio es de 2 m en el intervalo de tiempo considerado) es natural que la cohesión c y el peso volumétrico del suelo γ también se hayan modificado; sin embargo, la estimación realizada se concentra únicamente en las propiedades dinámicas del suelo que tienen influencia en la magnitud de las funciones de impedancia, por lo que con fines prácticos se aceptaron las magnitudes publicadas en ESVM (1988). Se hace notar que el procedimiento sugiere un incremento de la velocidad de propagación de ondas de cortante de casi el 20 por ciento en 25 años (Tabla 3).
Con esta información, se calcularon dos propuestas de cimentación de los edificios independientes entre sí para considerar los efectos de flexibilidad: a) un cajón de cimentación y b) pilotes de fricción. El pre-diseño por resistencia se realizó siguiendo los lineamientos de las Normas Técnicas Complementarias para el Diseño de Cimentaciones (NTCC-04).
El cajón se diseñó como una cimentación parcialmente compensada, donde la compensación corresponde a la diferencia entre la capacidad de carga del cajón como una losa y el peso total de la edificación. Los pilotes se diseñaron primordialmente a fricción y se ajustaron sus dimensiones para determinar números de pilotes que pudieran ser distribuidos en un arreglo simétrico. En ambos casos, el pre-diseño fue regido por la combinación de cargas verticales. En la Tabla 4 se muestran las principales características de las propuestas de los tres edificios considerados en esta investigación. Es posible conocer otros detalles del procedimiento de diseño de los sistemas de cimentación en Dávalos (2013) .
Edificio | Cajón de cimentación | Pilotes de fricción | |
---|---|---|---|
Profundidad de desplante, D | Cantidad de pilotes | Longitud del pilote | |
8 niveles | 4.00 m | 121 | 20.00 m |
12 niveles | 9.00 m | 169 | 22.50 m |
16 niveles | 15.00 m | 225 | 23.00 m |
Es valioso hacer notar que la distribución de los pilotes de fricción diseñados es simétrica (Figura 5), pero no es uniforme, puesto que para hacer una propuesta de diseño eficiente, se tomó en cuenta que los marcos contraventeados producen descargas verticales más grandes a la cimentación.
CÁLCULO DE LAS FUNCIONES DE IMPEDANCIA
En los análisis, las propiedades del sistema suelo - cimentación se consideraron a través de funciones de impedancia. Las funciones de impedancia son las rigideces dinámicas de la cimentación rígida sin masa, las cuales se definen como las fuerzas y momentos necesarios para producir desplazamientos y rotaciones armónicos unitarios en la cimentación. Las funciones son dependientes de la frecuencia de la excitación ω y se modelan como resortes y amortiguadores equivalentes a las propiedades mecánicas del suelo.
El análisis simplificado de la IDSE consiste en modelar el edificio sobre resortes y amortiguadores que representan la rigidez y el amortiguamiento del sistema suelo - cimentación. Con fines de diseño, los efectos de interacción suelo - estructura suelen tenerse en cuenta únicamente para el modo fundamental de vibración de la estructura. Una función de impedancia se expresa como se muestra en la ecuación 4 (Gazetas, 1991).
En la ecuación anterior, K*(ω) es la función de impedancia del sistema suelo - cimentación, K(ω) es la rigidez dinámica del sistema suelo - cimentación, C(ω) es el coeficiente dinámico de amortiguamiento viscoso equivalente y ω es la frecuencia angular de la excitación. La rigidez dinámica K(ω) representa la rigidez e inercia del suelo en el que se encuentra la cimentación, donde la inercia es dependiente de la frecuencia de excitación, mientras que la parte imaginaria de la ecuación 4 expresa la energía disipada a través de su amortiguamiento C(ω).
Es común definir distintas magnitudes de las funciones de impedancia de acuerdo a la dirección en que se excita el sistema suelo - cimentación para representarlo como un juego de resortes y amortiguadores viscosos como se muestra en la Figura 6. En la figura, Kr es la rigidez rotacional, Kh es la rigidez horizontal y Kv es la rigidez vertical que en este estudio se modeló como infinitamente rígida (Kv →∞) ya que no se considera la influencia de la componente vertical del sismo. Además, debido a que las estructuras se estudiaron ante análisis no lineales estáticos, el amortiguamiento rotacional Cr , horizontal Ch y vertical Cv no se consideraron.
En este estudio, las funciones de impedancia se determinaron mediante el programa Dyna5 (Novak et al., 1995), que permite calcular la respuesta de varios tipos de cimentaciones (rígidas o flexibles) ante una variedad de cargas dinámicas. El programa analiza las cimentaciones con diferentes profundidades de enterramiento y proporciona las funciones de impedancia para cualquier intervalo de frecuencias.
En el programa Dyna5 el cajón de cimentación se modeló como un esquema de stratum foundation, que considera el sistema como una cimentación superficial sobre un estrato homogéneo de espesor finito como se muestra en la Figura 7a. En la figura, Hs es el espesor del estrato, D es la profundidad de desplante y Zc es la posición del centro de masa del cajón. Por otra parte, la cimentación con pilotes de fricción se utilizó el modelo pile foundation, que permite analizar grupos de pilotes desplantados en un medio estratificado, aunque en el presente estudio, debido a que el medio es homogéneo sólo se utilizó un estrato (Figura 7b). Para el cálculo de la funciones de impedancia del grupo de pilotes se consideraron los efectos de grupo.
En el programa empleado, los resultados del análisis de las cimentaciones se presentan en gráficas de la rigidez dinámica del sistema suelo-cimentación modelado y su variación con la frecuencia. El programa calcula la rigidez dinámica de traslación horizontal y vertical, rigidez de rotación, rigidez de torsión y una rigidez acoplada. En esta investigación, se tomaron en cuenta únicamente las componentes de traslación horizontal Kh y de rotación Kr para las frecuencias fundamentales obtenidas de los análisis sobre base rígida (Tabla 2). Esta consideración supone que la mayor parte del movimiento estará asociada con el modo fundamental de vibrar de la estructura, como se propone en el planteamiento de las NTCS-04 MOCS-08. Así, en la Tabla 5 se reporta la magnitud de las funciones de impedancia para cada tipo de cimentación obtenidas de los análisis del Dyna5 (Novak et al 1995).
Modelo | Frecuencia fundamental de la estructura con base rígida (rad/s) | Pilotes de fricción | Cajón de cimentación | |||
---|---|---|---|---|---|---|
Rigidez dinámica horizontal, Kh (t/cm) | Rigidez dinámica rotacional, Kr (t-cm) | Rigidez dinámica horizontal, Kh (t/cm) | Rigidez dinámica rotacional, Kr (t-cm) | |||
Edificio de 8 pisos | 15.71 | 1.203(10)3 | 7.472(10)9 | 1.043(10)3 | 9.104(10)8 | |
Edificio de 12 pisos | 12.57 | 9.533(10)2 | 3.394(10)9 | 1.451(10)3 | 2.846(10)8 | |
Edificio de 16 pisos | 9.97 | 9.053(10)2 | 3.903(10)9 | 2.036(10)3 | 5.644(10)8 |
ANÁLISIS NO LINEALES
Se realizaron análisis no lineales ante carga estática monótona creciente con el programa SAP2000 v.15 (CSI 2011) ante un patrón de cargas laterales de forma triangular invertido siguiendo los criterios del método estático de las NTCS-04, que asume que el primer modo regirá la respuesta dinámica de las estructuras.
El programa utilizado realiza un modelado de la no linealidad de los elementos con base en la propuesta del FEMA 356 (FEMA 356, 2000) mediante curvas momento - rotación no lineal para vigas y columnas y/o carga - desplazamiento para contravientos como se muestra en la Figura 8. La magnitud de los parámetros que definen la curva de capacidad de los elementos (a, b y c) se establecen para elementos de acero a partir de tablas que estipula el FEMA 356, que se obtuvieron mediante pruebas experimentales y análisis numéricos. Por otra parte, la carga crítica de pandeo por compresión Pcr de los contravientos se definió dependiendo de la sección transversal por entrepiso para cada edificio como se discute en Tapia y Tena (2013).
El comportamiento no lineal de la estructura se estudió considerando la respuesta inelástica de los elementos estructurales; esto es, considerando que las vigas y columnas pueden desarrollar articulaciones plásticas ubicadas al 5 y 95% de la longitud del elemento y que los contravientos se plastificarían por fluencia en tensión o pandeo lateral en compresión ubicando la falla al 50 por ciento de la longitud del contraviento.
Modelo de la IDSE
En los modelos desplantados sobre base flexible (cimentación con un cajón y con pilotes de fricción), la IDSE se representa por medio de un resorte rotacional de rigidez Kr y un resorte traslacional horizontal de rigidez Kh ubicados en un nodo maestro en el centro de la base, donde la base del edificio se modela con un diafragma infinitamente rígido (Figura 6).
Debido a que en este estudio se analizarán marcos planos de los edificios estudiados y a que las funciones de impedancia (Tabla 5) corresponden al sistema suelo-cimentación del edificio completo, fue necesario determinar la magnitud proporcional de la rigidez del sistema suelo - cimentación (Kr y Kh ) que le corresponde a cada marco. Para ello, las constantes de los resortes que se emplearon en los análisis con base flexible se definieron considerando la proporción a la rigidez lateral total que aporta cada marco. Los porcentajes de rigidez lateral obtenidos de los tres tipos de marcos estudiados: marco sin contravientos, marco con una crujía contraventeada y marcos con dos crujías contraventeadas (Figura 2) se muestran en la Tabla 6. Es valioso hacer notar que la rigidez total del edificio (100%) se obtiene de la suma del total de los seis marcos que conforman cada dirección global. Estos porcentajes se calcularon mediante las rigideces laterales de cada marco obtenida en los análisis y la rigidez lateral de todo el edificio.
Modelo | Marco sin contravientos | Marco con una crujía contraventeada | Marco con dos crujías contraventeadas |
---|---|---|---|
Edificio de 8 pisos | 5.42% | 16.64% | 27.94% |
Edificio de 12 pisos | 6.11% | 16.34% | 27.55% |
Edificio de 16 pisos | 7.15% | 16.43% | 26.42% |
En la Tabla 7 se muestran los periodos de vibrar de los marcos de acero considerados en este estudio. En ella, BR se refiere a los análisis que suponen un sistema suelo-cimentación infinitamente rígido (base rígida) y BF-Pilotes y BF-Cajón se refiere a los análisis con base flexible usando pilotes de fricción y el cajón de cimentación, respectivamente.
En todos los casos la IDSE incrementa el periodo fundamental de los marcos analizados con base rígida. Cuando el edificio es de baja altura (modelo de 8 pisos), la estructura más flexible se obtiene mediante el modelado de una cimentación con cajón de cimentación; sin embargo, en edificios altos (el modelo de 16 pisos), la estructura más flexible se obtiene considerando una cimentación con pilotes por fricción. Esta observación denota una ligera dependencia de la flexibilidad de los modelos entre la altura de los edificios y el sistema de cimentación, siendo que el modelo de 12 pisos es el caso intermedio.
ANÁLISIS DE RESULTADOS
En la Figura 9 se muestran las curvas de capacidad globales de los marcos considerando el cortante basal contra la distorsión total (desplazamiento de la azotea entre la altura del edificio). En ellas, se muestra el desplazamiento total Ũ provocado por la deformación de la estructura U más la deformación asociada a la rotación de la cimentación ϕh y la deformación asociada a la de la cimentación UCIMENT (Figura 1). Estos resultados se utilizaron para compararlas contra las ecuaciones que evalúan la influencia de la flexibilidad de la cimentación en la ductilidad que pueden desarrollar las estructuras (NTCS-04; Avilés y Pérez-Rocha, 2005).
En las curvas de capacidad el último paso se definió como el incremento de carga cuando ocurría una de las siguientes condiciones:
Alguno de los elementos estructurales alcanzara la deformación última que puede desarrollar físicamente (rotación en trabes y columnas y/o alargamiento ó acortamiento en contravientos), definida como el punto C de la Figura 8 en el programa de análisis empleado.
Se forme un mecanismo de piso débil que se definió como la articulación en ambos extremos de todas las columnas de un mismo entrepiso.
Se desarrolle un mecanismo plástico de nudo por la articulación de los extremos de todos los elementos que concurren en él.
En las curvas de capacidad se nota que a medida que el sistema suelo-cimentación es más flexible, la rigidez elástica del sistema es menor (Figura 9). Asimismo, se aprecia una dependencia, que se había hecho notar anteriormente, entre la altura del modelo y el sistema de cimentación que es más flexible; es decir, el edificio de 8 pisos es más flexible cuando el modelo está resuelto con un cajón de cimentación; mientras que el edificio de 16 pisos está relacionado con mayores desplazamientos con una cimentación con pilotes de fricción. El edificio de 12 pisos tiene una respuesta de los modelos con cimentación flexible muy parecida.
Además, se calcularon curvas de capacidad que consideran únicamente la deformación de la estructura U (Figura 1) para comparar su magnitud con la obtenida de los análisis con base rígida, esto con la finalidad de estudiar si la flexibilidad de la base modifica el comportamiento asociado exclusivamente con la deformación de la estructura (Figura 10). Considerando únicamente este desplazamiento, resultó que la deformación de la estructura que desarrollan los marcos es muy similar entre los análisis que consideran una base rígida y una base flexible.
Con el propósito de hacer más notoria la modificación de la respuesta inelástica asociada a la IDSE cuando se toma en cuenta únicamente la deformación de la estructura U, se calculó la respuesta inelástica para distintas magnitudes hipotéticas de la rigidez rotacional Kr . Así, en la Figura 11 se muestran algunas de las curvas de capacidad calculadas en los análisis del marco sin contravientos y marco con dos crujías contraventeadas del edificio de 16 pisos. En ella, la curva de magnitud Kr representa el comportamiento del análisis con la función de impedancia real, mientras que el resto de las curvas representan situaciones hipotéticas de suelos muy flexibles representados como una proporción de la función de impedancia.
En la Figura 11, se nota que la cimentación con pilotes de fricción es la solución más flexible en el edificio de 16 pisos y esta tendencia, comentada anteriormente, parece ser independiente de la flexibilidad del suelo, aún hasta magnitudes hipotéticas de Kr. En general, la IDSE afecta la resistencia, rigidez y capacidad de deformación que pueden desarrollar los modelos únicamente en situaciones hipotéticas, lo que implica que los efectos P-Δ asociados a la IDSE tienen poca relevancia, al menos para los marcos de esta investigación, cuando se consideran las propiedades de suelo reales incluyendo suelos con poca rigidez rotacional.
Ductilidades desarrolladas
Debido a que hay diferencias más notorias entre los resultados de los marcos menos rígidos, la siguiente discusión se centra en los análisis de marcos sin contravientos y los marcos con una crujía contraventeada; sin embargo, es posible conocer los detalles de los resultados de los marcos con dos crujías contraventeadas en Dávalos (2013).
A partir de las curvas de capacidad se construyeron curvas primarias idealizadas (Figura 12) que son las que usualmente consideran los reglamentos para definir el comportamiento inelástico del sistema. Las curvas bilineales se determinaron con el cortante de fluencia Vy , la distorsión de fluencia δy y el cortante máximo Vu . En ellas, la rigidez inicial Kelast es la pendiente de la curva en el intervalo elástico, la rigidez post-fluencia se consideró con una pendiente igual a k2 = 0 y la distorsión de fluencia teórica δy se definió como el punto de intersección entre ambas rectas.
En la Tabla 8 se muestran las magnitudes de las ductilidades μ calculadas a partir de las curvas idealizadas bilineales de los análisis de los marcos de los modelos de 8, 12 y 16 pisos. Los cálculos se desarrollaron considerando las curvas de capacidad que consideran el desplazamiento total de los marcos (es decir, Ũ = U +ϕh +UCIMENT ) y las curvas que consideran únicamente la deformación de la estructura U. Los resultados hacen notar que la IDSE tiene influencia en la ductilidad μ que los marcos desarrollan debido a que se reduce la rigidez lateral del sistema y se modifica la capacidad de deformación. En la mayoría de los casos, los análisis de modelos con cimentaciones flexibles son menos dúctiles en comparación de los resultados de los análisis con base rígida.
Igualmente, se nota que las ductilidades calculadas considerando únicamente el desplazamiento de la estructura U son bastante similares a los resultados que se calcularon considerando base rígida. Esto implica que la reducción en la ductilidad calculada está asociada a la reducción en la rigidez elástica kElastic (Figura 1), previamente identificada en Pérez-Rocha y Avilés (2003). Sin embargo, es valioso hacer notar que los modelos de 12 y 16 niveles con una crujía contraventeada no siguen esa tendencia. Estos resultados coinciden con otras investigaciones (Ganjavi y Hao 2011; Tapia y Tena 2013) que sugieren una dependencia de la ductilidad que las estructuras pueden desarrollar en función de su altura. De hecho, esta dependencia entre la ductilidad y la altura de los edificios ya se ha incluido en el reglamento de Canadá para el diseño sísmico de estructuras metálicas (Izvernari et al. 2007).
Asimismo, es posible concluir que las distorsiones obtenidas de análisis con base rígida son una buena representación de las distorsiones que se obtendrían en análisis con base flexible cuando se considera la deformación de la estructura de la estructura U únicamente. Pese a esto, cabe mencionar que el modelo con base rígida subestima el desplazamiento total que podría desarrollar una estructura con una cimentación flexible.
Ahora, las NTCS-04 proponen un criterio para estimar la ductilidad que desarrollará una estructura, a partir de su respuesta con base rígida mediante la aplicación de la ecuación 1. En la Tabla 9 se muestran las ductilidades estimadas para estructuras con base flexible considerando este procedimiento. En los cálculos se tomaron en cuenta los periodos fundamentales de la estructura con base rígida Te y los obtenidos para los modelos con base flexible mostrados en la Tabla 7 y las ductilidades con base rígida μ mostradas en la Tabla 8. En la Tabla 9, p NTC y c NTC indica la ductilidad estimada mediante el criterio de las NTCS-04 para el caso de una cimentación con pilotes por fricción y con un cajón de cimentación, respectivamente. En el cálculo sólo se utilizó la expresión para la modificación de la ductilidad incluida en las NTCS-04 (ecuación 1). La determinación de las magnitudes del periodo con base rígida (
Con base en los resultados, el criterio de estimación de las NTCS resulta ser una buena aproximación de la ductilidad que desarrollan los marcos sin contravientos. Sin embargo, ante el aumento de rigidez lateral (marcos con una crujía contraventeada), los resultados indican una tendencia desfavorable, que parece ser dependiente de la altura de los modelos.
Aunado a lo anterior, las NTCS-04 establecen un criterio para determinar la importancia de la influencia de la IDSE en el comportamiento de la estructura, a partir de la rigidez relativa de las estructuras respecto a la rigidez del sistema suelo - cimentación (ecuación 5). Según este planteamiento, los efectos de la IDSE son relevantes y no deberían de despreciarse cuando la magnitud del cociente es menor a
En la ecuación anterior, Hs es el espesor del estrato, hes la altura de la estructura para el caso de estructuras de un solo grado de libertad (Figura 1). Para el caso de estructuras de varios grados de libertad, la altura h corresponde a la altura efectiva en el primer modo. Ts y Te son el periodo dominante del estrato de suelo y de la estructura, respectivamente. En la Tabla 10 se muestran los cálculos obtenidos para evaluar la importancia de los efectos de la interacción suelo - estructura en los edificios estudiados en esta investigación. En los cálculos el periodo del suelo se consideró igual a Ts = 1.85 s (Tabla 3); el espesor del estrato igual a Hs = 40 m considerando el estudio de Mecánica de Suelos; los periodos Te de los edificios mostrados en la Tabla 2 y siguiendo el planeamiento de las NTCS la altura efectiva de la estructura se consideró igual a h= 0.7 HT , donde HT es la altura total del edificio de la Tabla 2.
Modelo | h= 0.7HT | λ |
---|---|---|
Edificio de 8 pisos | 19.6 m | 0.441 |
Edificio de 12 pisos | 29.4 m | 0.368 |
Edificio de 16 pisos | 39.2 m | 0.348 |
Con base en los resultados, los efectos de la IDSE son relevantes en los edificios estudiados localizados en suelos con poca rigidez rotacional (suelo blando). Asimismo, se nota nuevamente que los efectos son dependientes de la altura de los edificios.
CONCLUSIONES
El artículo presenta los resultados de análisis no lineales ante carga estática monótona creciente de marcos de acero de edificios regulares para estudiar la influencia de la interacción dinámica suelo - estructura (IDSE) en su comportamiento inelástico. En el estudio se eligieron marcos sin contravientos, con una crujía contraventeada y con dos crujías contraventeadas para considerar la influencia de la rigidez lateral en los resultados.
Los edificios estudiados, que tienen 8, 12 y 16 pisos, se supusieron localizados en la zona IIIb para asociarlo a las mayores demandas de aceleración que proponen las Normas Técnicas (NTCS-04) y fueron diseñados siguiendo estrictamente los criterios del Reglamentos de Construcciones para el Distrito Federal (RCDF-04). A fin de estimar los efectos de la IDSE en la respuesta inelástica se consideraron tres casos: a) análisis de los edificios considerando un sistema suelo-cimentación infinitamente rígido (base rígida), b) análisis con una cimentación con pilotes de fricción y c) análisis con un cajón de cimentación.
El diseño de la cimentación se realizó mediante los criterios de las Normas Técnicas (NTCC-04) y usando las propiedades estratigráficas de un estudio de la Colonia Roma, Delegación Cuauhtémoc en congruencia con las suposiciones del diseño elástico de los edificios. El estudio de los análisis con cimentación flexible se realizó desde dos perspectivas, considerando: a) únicamente el desplazamiento relacionado con la deformación de la estructura U y b) el desplazamiento total (Ũ= U +ϕh +UCIMENT ); es decir, considerando el desplazamiento relacionado con la deformación de la estructura U más el giro de la cimentación ϕh más el desplazamiento de la cimentación UCIMENT .
En los resultados obtenidos, se nota una ligera dependencia de la altura de los modelos He y de la rigidez de los marcos en la reducción de la capacidad de deformación inelástica asociada a la IDSE. Con base en los resultados, se nota que las ductilidades calculadas considerando únicamente el desplazamiento de la estructura U son bastante similares a los resultados que se calcularon considerando base rígida cuando se consideran marcos de acero sin contravientos; sin embargo, en edificios altos con marcos con contravientos no se sigue esa tendencia. Esta observación tiene relación con los resultados de otras investigaciones que señala que la ductilidad que desarrollan las estructuras es dependiente de la altura de los edificios.
El efecto de la IDSE modifica la ductilidad teórica que pueden desarrollar los edificios, dado que en la mayoría de los casos, los modelos con cimentaciones flexibles reportan una disminución de la ductilidad que puede desarrollar el sistema en comparación de los resultados con base rígida cuando se considera el desplazamiento total de los edificios (Ũ = U +ϕh +UCIMENT ). Sin embargo, las distorsiones obtenidas de análisis con base rígida son una buena representación de las distorsiones que se obtendrían en análisis con base flexible cuando se considera la deformación de la estructura de la estructura U únicamente. Pese a esto, se demostró que el modelo con base rígida subestima el desplazamiento total que podría desarrollar una estructura con una cimentación flexible.
Los resultados de los análisis inelásticos se compararon con la estimación de la ductilidad que desarrollará una estructura según las NTCS-04. Los resultados indican que el criterio de las NTCS-04 es una buena aproximación de la ductilidad que desarrollan los marcos de acero sin contravientos. Sin embargo, ante el aumento de rigidez lateral (marcos con contravientos) los resultados indican una tendencia desfavorable, que parece ser dependiente de la altura de los modelos.
Los resultados del estudio muestran que los efectos de la IDSE pueden ser relevantes en edificios altos estructurados con marcos de acero localizados en suelo blando dada la modificación de su desempeño inelástico. Este estudio representa una referencia racional de la variación de la rigidez inducida por la flexibilidad de la cimentación y constituyen un congruente indicio de la respuesta inelástica de la edificación en función de la IDSE; sin embargo, será necesario realizar otros estudios con una mayor cantidad de modelos antes de realizar conclusiones definitivas.