Sociedad Mexicana de Ingeniería Estructural

CONFIABILIDAD SÍSMICA DE EDIFICIOS CON EXCENTRICIDADES

TORSIONALES DE RIGIDEZ Y DE RESISTENCIA

1 2 3
Yasser Picazo Gama , Orlando Javier Díaz López y Luis Esteva Maraboto

RESUMEN

Se estudia la influencia de las excentricidades de rigidez y de resistencia en la respuesta sísmica de edificios

de comportamiento no lineal ante la acción de dos componentes horizontales ortogonales simultáneas del
movimiento del terreno. Se desarrollan los criterios y la metodología necesaria para la obtención de la
confiabilidad de los sistemas estructurales, así como criterios para definir sistemas simplificados de referencia
que pueden ser utilizados para estudios paramétricos. Se analiza la respuesta sísmica de un conjunto de
edificios con distintas configuraciones tridimensionales (simétricos y asimétricos). Los resultados de estos
análisis se utilizan para estimar funciones de confiabilidad de los sistemas asimétricos y se comparan con la
confiabilidad que presentan los sistemas simétricos. Estos resultados pueden ser utilizados para obtener
factores de reducción de ordenadas espectrales que deban aplicarse para el diseño sísmico de los sistemas
asimétricos a fin de lograr los mismos niveles de confiabilidad que los sistemas simétricos.

ABSTRACT

A study is presented of the influence of stiffness and strength eccentricities on the seismic response of
building systems with nonlinear behavior under the action of two simultaneous orthogonal horizontal ground
motion components. Attention is given to the criteria and methods necessary to estimate the seismic
reliability functions of the structural systems considered, as well as to define simplified reference systems that
can be applied in parametric studies. The seismic response of a set of buildings is analyzed, considering
different configurations (symmetric and asymmetric). The results are applied to the estimation and
comparison of the reliability functions of symmetric and asymmetric systems. These results can be applied to
the determination of reduction factors of spectral ordinates to be used for the design of asymmetric systems,
in order to attain the same reliability levels as those of the symmetric ones.

INTRODUCCIÓN

Los criterios de diseño sísmico vigentes para edificios con planta asimétrica se basan en análisis dinámicos de
sistemas lineales. Para tomar en cuenta la acción simultánea de las dos componentes horizontales ortogonales
del movimiento del terreno, junto con la influencia de las excentricidades de rigidez y de resistencia, se parte
de un análisis de respuesta para cada componente por separado, incluyendo una excentricidad accidental, y se
aplica un criterio de superposición que toma en cuenta que los máximos de las respuestas asociadas a cada
componente individual no ocurren simultáneamente.

Buena parte de los edificios que se construyen presentan asimetrías en planta, tanto en rigideces como en
resistencias laterales; estos edificios están expuestos a movimientos sísmicos que incluyen varias

1
Instituto de Ingeniería, UNAM, Av. Universidad No. 3000, Col. Universidad Nacional Autónoma de
México, CU, Delegación Coyoacán, 04510, México, D.F. Teléfono, (55) 5623-3600 ext. 8479; fax: (55)
5623-3600 ext. 8050; YPicazoG@iingen.unam.mx
2
Instituto de Ingeniería, UNAM, Av. Universidad No. 3000, Col. Universidad Nacional Autónoma de
México, CU, Delegación Coyoacán, 04510, México, D.F. Teléfono, (55) 5623-3600 ext. 8472; fax: (55)
5623-3600 ext. 8050; odil@pumas.iingen.unam.mx
3
Instituto de Ingeniería, UNAM, Av. Universidad No. 3000, Col. Universidad Nacional Autónoma de
México, CU, Delegación Coyoacán, 04510, México, D.F. Teléfono, (55) 5623-3653; fax: (55) 5623-3600 ext.
8050; LEstevaM@iingen.unam.mx

1
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

componentes simultáneas; cuando menos, en dos direcciones horizontales ortogonales. La influencia del
comportamiento no lineal en la respuesta dinámica tri-dimensional de uno de estos sistemas no puede tomarse
en cuenta mediante la aplicación directa de los factores de reducción de la capacidad lateral requerida
determinados para sistemas planos o simétricos, libres de excentricidades en planta. En los sistemas
asimétricos sujetos a dos componentes horizontales simultáneas del movimiento del terreno, la respuesta
dinámica no lineal depende, además de las rigideces y resistencias laterales en cada una de dos direcciones
ortogonales, de las excentricidades en rigidez y resistencia en cada dirección, así como de la rigidez torsional
del sistema con respecto a un eje vertical, y de su capacidad para resistir momentos de torsión con respecto al
mismo eje. Además, al estimar la capacidad lateral de cada elemento vertical, como una columna o un núcleo
de elevadores, debe tomarse en cuenta la interacción entre las fuerzas internas en cada uno de estos elementos
ante la acción de las dos componentes simultáneas del movimiento del terreno.

Los factores de reducción de ordenadas espectrales que se proponen en los criterios convencionales de diseño
sísmico se basan en los resultados de estudios sobre la respuesta dinámica no lineal de sistemas planos, o de
sistemas simétricos con respecto a un plano vertical, sometidos a una componente horizontal del movimiento
del terreno alojada en el plano de simetría. Es por lo tanto necesario realizar estudios sistemáticos sobre la
respuesta sísmica no lineal de edificios asimétricos en planta, sujetos a dos componentes horizontales
ortogonales simultáneas, con el fin de contar con información sobre los valores de los factores de reducción,
como los mencionados arriba, que deberían aplicarse para el diseño sísmico de tales edificios.

El estudio que se presenta en este trabajo tiene como objetivo evaluar la influencia de las excentricidades de
rigidez y de resistencia en la respuesta sísmica de edificios de comportamiento no lineal ante la acción de dos
componentes horizontales ortogonales simultáneas del movimiento del terreno. Para ello se desarrollan los
criterios y la metodología necesaria para la obtención de la confiabilidad de los sistemas estructurales. Se
estudia la respuesta sísmica de un conjunto de edificios con distintas condiciones de efectos tridimensionales.
Los resultados del estudio pueden emplearse para: obtener factores de transformación para estimar las
distorsiones laterales máximas de los elementos de borde de los sistemas con planta asimétrica a partir de los
correspondientes sistemas simétricos, funciones que relacionen los niveles de confiabilidad sísmica con el
valor de la intensidad, en términos de los parámetros que definan las propiedades globales del sistema, así
como factores de reducción de ordenadas espectrales que deben aplicarse para el diseño sísmico de los
sistemas asimétricos, a fin de lograr los mismos niveles de confiabilidad que corresponden a los sistemas
simétricos.

Además, se definen sistemas simples de referencia (SSR) los cuales pueden utilizarse para realizar estudios
paramétricos, cuyos resultados permitan establecer criterios generales para estimar valores de la respuesta
máxima de sistemas estructurales con excentricidades torsionales de rigidez y de resistencia. Las propiedades
de estos sistemas se obtienen a partir de las propiedades globales del sistema detallado de interés.

Con los resultados de este estudio se pretende desarrollar criterios aplicables en condiciones prácticas para
diseño sísmico con niveles uniformes de confiabilidad.

PLANTEAMIENTO

SISTEMAS DETALLADOS

El estudio de la influencia de las excentricidades de rigidez y de resistencia sobre la respuesta dinámica global
y las funciones de confiabilidad sísmica de edificios de comportamiento no lineal, se realiza sobre un
conjunto de edificios detallados que presentan diferentes características de excentricidades de rigidez y de
resistencia. Estos sistemas se someten a dos componentes horizontales ortogonales simultáneas del
movimiento del terreno. Los resultados de estos análisis se utilizan para estimar funciones de valores
máximos de la respuesta de los marcos (de borde e internos) que conforman la estructura, así como de su
respuesta global.

2
 Sociedad Mexicana de Ingeniería Estructural
Índices de confiabilidad basados en IRRS

A partir de estudios de respuesta dinámica paso-a-paso, se determinan valores del índice de reducción de
rigidez secante, IRRS (Esteva e Ismael, 2004)
k0  k
I RRS  (1)
k0
Donde K0 es la rigidez tangente inicial asociada con la curva cortante basal-desplazamiento de azotea, que
resulta de un análisis de empuje lateral y K es la rigidez secante (cortante basal entre desplazamiento de
azotea lateral) cuando el desplazamiento de azotea lateral alcanza su valor máximo absoluto durante la
respuesta sísmica del sistema. La condición de falla se alcanza cuando IRRS=1.0. Cuando la estructura se
comporta en forma elástica IRRS =0.0.

Picazo (2012) ha estudiado la aplicación del criterio basado en la determinación del índice de reducción de
respuesta sísmica (IRRS) para la estimación de las funciones de confiabilidad sísmica, en el caso de sistemas
asimétricos en planta sujetos a dos componentes sísmicas horizontales ortogonales simultáneas. Con el fin de
evitar ciertas incongruencias en su aplicación propone el siguiente criterio.

Sean k0x, k0y, kx, ky, las rigideces iniciales y las rigideces secantes para la distorsión máxima en las direcciones
(x, y), respectivamente, referidas a la fuerza cortante en la base para el sistema completo y al movimiento del
centro de gravedad del nivel superior; k0y k la rigidez rotacional inicial (momento torsionante en la base,
dividido entre el giro del extremo superior con respecto a un eje vertical) y la rigidez rotacional secante para
el instante en que ocurre la máxima rotación del extremo superior del sistema.

Definimos ahora IRRS como el máximo de los siguientes valores:

k0 x  k x
(2a)
k0 x

k0 y  k y
(2b)
k0 y

k 0  k
(2c)
k 0

En teoría este criterio es adecuado; sin embargo, para el cálculo del valor del índice IRRS para el grado de
libertad rotacional, se tiene el inconveniente de que cuando el sistema se encuentra en el intervalo inelástico,
el centro de rigidez del sistema va cambiando de acuerdo a como fluyen los diferentes elementos que forman
los diferentes marcos que componen la estructura. Debido a esto, es necesario recalcular en cada momento en
que suceda un cambio de rigidez en los elementos estructurales el centro de rigidez del sistema, para obtener
el momento torsionante que se presenta y con ello calcular la rigidez rotacional instantánea (secante). Dado lo
anterior, en este estudio no se toma en cuenta el valor de IRRS en rotación y solo se considera el valor máximo
de tal índice en las direcciones x y y. Se piensa que tomarlo de esa forma no es inadecuado, ya que en general
las estructuras difícilmente fallan por momento torsionante, sin que ocurra simultáneamente el colapso por
deformaciones excesivas en traslación.

Relación entre IRRS y margen de seguridad

Para la estimación de la función de confiabilidad sísmica, (y), de sistemas estructurales, se utiliza el índice
de confiabilidad  de Cornell (1969), el cual se define como la relación del valor medio a la desviación
estándar del margen de seguridad. Aquí y es la intensidad del movimiento del suelo. La función de
confiabilidad se obtiene mediante la determinación de los indicadores de segundos momentos de la función de
distribución de probabilidades de ZF = lnYF, donde YF es el valor mínimo de la intensidad que produce el
colapso de la estructura (IRRS=1.0). El margen de seguridad, ZM, se define como el logaritmo natural de la

3
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

relación de la capacidad del sistema a la amplitud de su respuesta, dada una intensidad y. En este caso dicho
margen de seguridad se puede tomar como
Y 
Z M  ln  F  (3)
 y 
Por lo tanto, la función de confiabilidad se puede expresar como
E Z F   ln( y )
 ( y)  (4)
Z F

E(∙) y (∙) son la media y la desviación estándar, respectivamente.

La función de confiabilidad, (y), se pude obtener a partir de una muestra de pares de valores de Z y IRRS, los
cuales se usan para la estimación de la media y la desviación estándar de Z(u), definida como el logaritmo
natural de la intensidad aleatoria Y que corresponde al valor de IRRS=u. Los valores de Z en la muestra que
corresponden a valores de IRRS=1.0 son límites superiores de ZF=lnYF. Por lo tanto, los puntos
correspondientes a estos valores de ZF no pueden incorporarse en la muestra para la estimación de E(Z(u)) y
(Z(u)) utilizando un análisis de regresión de mínimos cuadrados convencional. Sin embargo, ellos pueden
ser incluidos en el proceso de estimación, mediante un análisis de máxima verosimilitud (Díaz y Esteva, 2006
y Díaz, Ismael y Esteva, 2008). En caso de que se utilice este último procedimiento, el valor esperado y la
desviación estándar de Z(u), se representan ahora por las funciones E(Z(u)|1) y (Z(u)|2), de forma dada y
parámetros 1 y 2 desconocidos; éstos se estiman usando una muestra que incluye m valores de IRRS menores
que 1.0 y n valores iguales a 1.0. Para la muestra observada, la función de verosimilitud de ( 1, 2) puede
expresarse como
m n
L(1 ,  2 )   f Z ( zi | ui , 1 ,  2 ) FZ ( z j | u j , 1 ,  2 ) (5)
i 1 j 1

SISTEMAS SIMPLIFICADOS DE REFERENCIA

Con el fin de establecer criterios generales para estimar valores de la respuesta máxima de sistemas
estructurales con excentricidades torsionales de rigidez y de resistencia, se utilizan sistemas simples de
referencia (SSR). Las propiedades de estos sistemas se obtienen a partir de las propiedades globales del
sistema de interés con propiedades medias (Picazo, 2012).

C22

C12

C00

C11
C21

Figura 1 Sistema Simplificado de Referencia “SSR”

El SSR se representa mediante un sistema de un solo nivel, con una cubierta horizontal, con tres grados de
libertad en el plano de la cubierta: dos de traslación y uno de rotación con respecto a un eje vertical (figura 1).
Las propiedades mecánicas de los miembros verticales que soportan la cubierta se representan mediante
curvas de fuerza cortante en la base vs desplazamientos laterales, y de momentos de torsión vs rotación de la

4
 Sociedad Mexicana de Ingeniería Estructural
cubierta con respecto a un eje vertical ante carga cíclica (funciones V vs δ y MT vs θ), de manera de
representar las propiedades globales de modelos detallados de múltiples niveles, tomando en cuenta en forma
aproximada la interacción entre las fuerzas cortantes en las columnas ante la acción de las dos componentes
horizontales ortogonales. Mediante el empleo de estos modelos se obtienen relaciones entre las distorsiones
laterales máximas en los elementos resistentes de borde de los sistemas. Por razones de simplicidad, la
degradación de rigideces y de resistencias ante cargas cíclicas se representa en todos los casos mediante el
mismo modelo, tomando los parámetros que definen dicha degradación de manera que conduzcan a funciones
representativas de las propiedades globales de sistemas típicos.

Las propiedades globales de los sistemas se representan mediante los siguientes parámetros:

a) Frecuencias naturales en traslación en las dos direcciones ortogonales principales, x, y.

b) Relación entre la frecuencia natural fundamental en traslación y en torsión pura, con respecto a un
eje vertical que pase por el centro de masa
c) Función V vs δ en cada dirección (x, y); para los casos de estructuras con excentricidades de rigidez y
de resistencia: al determinar estas curvas se tomarán en cuenta los tres grados de libertad del
desplazamiento de la cubierta.
d) Aportación de los elementos resistentes en cada dirección a la función MT vs θ.
e) Excentricidad de rigideces en cada una de las direcciones horizontales consideradas.
f) Excentricidad de resistencias en cada una de las direcciones horizontales consideradas.
Criterios para determinar propiedades de modelos simplificados (Picazo, 2012)

SSR para edificios con planta rectangular con asimetrías en rigidez y resistencia

El sistema simplificado de referencia (SSR) está integrado por cinco columnas, una en cada esquina y una en
el centro, designadas como Cij, con j = 0, 1, 2 (figura. 1). Cada columna tiene rigidez y resistencia en dos
direcciones ortogonales. Tales propiedades se designan como sigue:

Kijx = rigidez lateral de la columna Cij en la dirección paralela al eje x.

Rijx = resistencia lateral de la columna Cij en la dirección paralela al eje x.
Kijy = rigidez lateral de la columna Cij en la dirección paralela al eje y.
Rijy = resistencia lateral de la columna Cij en la dirección paralela al eje y.

Las propiedades del sistema simple de referencia se obtienen a partir de las propiedades globales del sistema
de interés. Estas propiedades se determinan mediante dos análisis ELSE (empuje lateral seudo-estático) de
traslación, en cada una de las dos direcciones ortogonales de interés, imponiendo restricciones laterales que
impidan el giro y el desplazamiento de las plantas del edificio en una dirección diferente a la del empuje
lateral. De estos análisis se obtienen las funciones V0 vs δN de fuerza cortante en la base vs desplazamiento en
el extremo superior (azotea) de cada subsistema (marco, sub-sistema dual) resistente en la dirección de
interés. Las pendientes de estas curvas en el intervalo de comportamiento lineal representarán las rigideces
laterales de cada uno de los sub-sistemas 2D (marcos o sub-sistemas duales) contenidos en planos paralelos a
la dirección de interés. La línea de acción de las resultantes de tales rigideces lineales en la base del sistema
estructural considerado determina una de las coordenadas del centro de torsión. Realizando un análisis
semejante en la dirección ortogonal, se determina el centro de torsión en la base del sistema.

Determinación de funciones V() para cada uno de los elementos resistentes del SSR

Las funciones V(δ) de los elementos resistentes del SSR en cada dirección se determinan a partir de los
resultados de un análisis ELSE (empuje lateral seudo-estático) del modelo detallado, mediante las curvas
obtenidas de la fuerza cortante resistente en la base contra el desplazamiento lateral de azotea, Vivsδi, para
cada marco resistente, así como también de la curva de capacidad global que relaciona la fuerza cortante basal
contra el desplazamiento de azotea Vxvsδ0x, o Vyvsδ0y, según el caso.

De acuerdo con esto, se pueden establecer las siguientes relaciones:

5
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

Dirección x:

M x1  0 x   Vix ( 0 x )  yi ; para yi  0 (6a)

M x 2  0 x   Vix ( 0 x )  yi ; para yi  0 (6b)

Para los elementos del SSR se tiene:

M x1  0 x   Vx1 ( 0 x )  b / 2 (7a)
M x 2  0 x   Vx 2 ( 0 x )  b / 2 (7b)
Vx 0  0 x   Vx ( 0 x )  Vx1 ( 0 x )  Vx 2 ( 0 x ) (7c)

Dirección y:

M y1  0 y   Viy ( 0 y )  xi ; para xi  0 (8a)

M y 2  0 y   Viy ( 0 y )  xi ; para xi  0 (8b)

Para los elementos del SSR se tiene:

M y1  0 y   Vy1 ( 0 y )  a / 2 (9a)
M y 2  0 y   Vy 2 ( 0 y )  a / 2 (9b)
Vy 0  0 y   Vy ( 0 y )  Vy1 ( 0 y )  Vy 2 ( 0 y ) (9c)

En las expresiones anteriores:

xi es la abscisa del plano vertical que contiene al i-ésimo marco paralelo a la dirección y, con respecto al
origen definido por el eje vertical que contiene a los centros de masa de las plantas del sistema detallado.
yi es la ordenada del plano vertical que contiene al i-ésimo marco paralelo a la dirección x, con respecto al
origen definido por el eje vertical que contiene a los centros de masa de las plantas del sistema detallado.
Mx1(δ0x), Mx2(δ0x) son las aportaciones al momento torsionante en la base, correspondientes a los dos marcos
extremos del SSR en la dirección x, expresadas como función del desplazamiento lateral de azotea referido
al centro de masa.
My1(δ0y), My2(δ0y) son las aportaciones al momento torsionante en la base, correspondientes a los dos marcos
extremos del SSR en la dirección y, expresadas como función del desplazamiento lateral de azotea referido
al centro de masa.
Vx1(δ0x), Vx2(δ0x) son las funciones de la fuerza cortante resistente en la base contra el desplazamiento lateral
de azotea referido al centro de masa, correspondientes a los dos marcos extremos del SSR en la dirección
x.
Vy1(δ0y), Vy2(δ0y) son las funciones de la fuerza cortante resistente en la base contra el desplazamiento lateral
de azotea referido al centro de masa, correspondientes a los dos marcos extremos del SSR en la dirección
y.
Vx0(δ0x) es la función de la fuerza cortante resistente en la base contra el desplazamiento lateral de azotea,
correspondiente al elemento central del SSR, en la dirección x.
Vy0(δ0y) es la función de la fuerza cortante resistente en la base contra el desplazamiento lateral de azotea,
correspondiente al elemento central del SSR, en la dirección y.

Dada la función Vx1(δ0x), pueden determinarse las funciones que ligan el momento flexionante en un plano
paralelo a la dirección x, en cada extremo de las columnas C11 y C21, con los giros de las articulaciones
plásticas correspondientes. A cada una de estas funciones puede ajustarse una expresión bilineal, de donde se

6
 Sociedad Mexicana de Ingeniería Estructural
obtendrán las propiedades K11x, K21x, R11x y R21x. De igual manera pueden obtenerse las propiedades Kijx, Rijx,
Kijy y Rijy en cada columna mostrada en la figura 1.

MOVIMIENTOS SÍSMICOS

Para los análisis de respuesta sísmica de los sistemas analizados, se utiliza una familia de movimientos
sísmicos simulados. El método de simulación se puede encontrar en Esteva et al (2010). El procedimiento
descrito en dicho trabajo permite simular las dos componentes ortogonales con intensidades reales aleatorias,
simuladas de manera de satisfacer un valor especificado del coeficiente de correlación entre ellas.

Intensidad sísmica

En este trabajo se toma como medida de intensidad de los movimientos sísmicos a la media cuadrática de las
ordenadas espectrales en dos direcciones ortogonales, Sa.

2
S aEO  S aNS
2
Sa  (10)
2

Donde SaEO y SaNS son la ordenada del espectro de seudo-aceleraciones para un amortiguamiento de 5% del
movimiento sísmico en la dirección este-oeste y norte-sur, respectivamente, para un periodo de vibración
dado.

Para la presentación de los resultados de la respuesta sísmica de los sistemas estudiados, se utiliza como
medida de la intensidad sísmica a una intensidad normalizada, q. Ésta se define como la media cuadrática de
las ordenadas espectrales en dos direcciones ortogonales, Sa, normalizada con respecto a la aceleración de
fluencia del sistema en la dirección en que se tiene la resistencia mínima

Sa
q (11)
Vy m

Donde y es el cortante basal de fluencia mínimo del sistema con propiedades medias (de entre las dos
direcciones ortogonales) y es la masa traslacional total media del sistema.

APLICACIONES

Los modelos detallados tienen una configuración geométrica en planta de forma rectangular, de dimensiones
ancho y largo de 16.5m y 22.0m, en todos sus niveles. Cada uno de ellos tiene una altura total de 17m (6
niveles); la altura del primer nivel es de 3.0m y la de los niveles superiores es de 2.8m. Los modelos poseen
una estructuración a base de marcos de concreto reforzado, con tres y cuatro crujías en las direcciones x e y,
respectivamente (figura 2).

En este estudio, se efectúo el análisis y diseño estructural de cuatro sistemas detallados con base en el
Reglamento de Construcciones del Distrito Federal (RCDF2004) mediante el programa ECOgcW (Corona,
2007), considerando para ello, los valores medios tanto de las cargas gravitacionales como de las propiedades
mecánicas de los elementos del sistema. El primero de ellos fue diseñado de tal forma que no tuviera
excentricidad de rigidez en todos sus niveles, en ambas direcciones ortogonales, incluyendo la excentricidad
accidental; los otros tres, fueron obtenidos a partir de este sistema simétrico de partida, los cuales fueron
diseñados con una excentricidad torsional de rigideces de esx=0.2 y esy=0.1 en cada uno de sus niveles, y
también, con un valor de la excentricidad accidental igual a cero. Lo anterior se logra a partir de aplicar en el
sistema simétrico de partida ciertos factores que permitan cambiar la rigidez lateral de los marcos de concreto
reforzado, al cambiar las dimensiones de columnas y trabes, de tal manera, de satisfacer simultáneamente en
ambas direcciones ortogonales valores especificados de las excentricidades de rigidez.

7
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

Dichos sistemas se designaron como SS6-EDS, SA6-EDS, SA6-EDSR y SA6-EDSA respectivamente. Para el
diseño de los dos primeros, se utilizó el espectro de diseño de seudo aceleraciones (EDS) que se establece en
las normas, mientras que, para el caso de los dos últimos, se usó el EDS correspondiente, pero con ordenadas
reducidas a la mitad y aumentadas al doble, respectivamente. En la tabla 1 se muestran los valores obtenidos
de la excentricidad estática de entrepiso de los sistemas estudiados, y en la tabla 2 se proporcionan los
parámetros empleados para su análisis dinámico modal espectral.

x CT

x
Figura 2 Vista en planta de los sistemas detallados

Tabla 1 Valores obtenidos de la excentricidad estática de entrepiso de los sistemas de estudio

SA6 diseñado con EDS SA6 diseñado con EDS reducido SA6 diseñado con EDS
Sistema
(SA6-EDS) a la mitad (SA6-EDSR) aumentado al doble (SA6-EDSA)
Nivel esx/a esy/b esx/a esy/b esx/a esy/b
1 0.209 0.106 0.205 0.097 0.202 0.101
2 0.208 0.104 0.197 0.098 0.205 0.113
3 0.212 0.100 0.197 0.101 0.204 0.108
4 0.217 0.097 0.210 0.099 0.202 0.110
5 0.218 0.098 0.210 0.101 0.206 0.113
6 0.201 0.096 0.188 0.096 0.198 0.089

Tabla 2 Parámetros utilizados para realizar el análisis modal espectral de los sistemas

Sistema Zona Grupo Q k e Ts (s) Cs (Sa/g) a0 (Sa/g) Ta (s) Tb (s)

SS6-EDS IIIb B 2 1 0.0*b 2.0 1.2 0.25 1.175 2.4
SA6-EDS IIIb B 2 0.9 0.0*b 2.0 1.2 0.25 1.175 2.4
SA6-EDSR IIIb B 2 0.9 0.0*b 2.0 0.6 0.25 1.175 2.4
SA6-EDSA IIIb B 2 0.9 0.0*b 2.0 2.4 0.25 1.175 2.4
En esta tabla, b es la dimensión en planta perpendicular a la dirección de análisis

Donde:

Q Factor de comportamiento sísmico

k Factor reductivo de Q en función de la regularidad
e Excentricidad accidental
Ts Periodo dominante más largo del terreno
C s Coeficiente sísmico correspondiente a la meseta del espectro de seudo aceleraciones
a 0 Aceleración máxima del terreno
Ta Periodo característico del espectro de diseño
Tb Periodo característico del espectro de diseño

8
 Sociedad Mexicana de Ingeniería Estructural
Para efecto de estimar las funciones de confiabilidad de los sistemas detallados, es preciso realizar el análisis
de empuje lateral seudo estático “ELSE” y el análisis dinámico no lineal paso a paso de dichos sistemas con
propiedades medias. Por tal motivo, se hace uso del programa de análisis no lineal de estructuras Ruaumoko
3D (Carr, 2007), en donde, para el análisis y el modelado de los sistemas se consideró lo siguiente: el sistema
de piso se modeló como diafragma rígido, los elementos estructurales se modelaron como elementos viga
(modelo viga de un componente de Giberson); también, se consideró la reducción de la capacidad resistente
de las columnas debido a la interacción de los momentos actuantes y resistentes en los extremos de las
mismas, en ambas direcciones ortogonales, producidos por la acción de dos componentes horizontales
ortogonales simultáneas, mediante la ecuación dada por Bresler (1960):
 
 Mx   My 
      1.0 (12)
 
 x0   y0 
M M

En esta ecuación, el valor adoptado del parámetro α es de 1.5. Asimismo, para representar la matriz de
amortiguamiento de la estructura se seleccionó el modelo de amortiguamiento de Rayleigh, el cual se define
como
C  M   K (13)

Aquí, M y K son las matrices de masa y rigidez de la estructura, y α y β son parámetros que permiten definir
el nivel de amortiguamiento viscoso en dos frecuencias de vibrar diferentes, que generalmente son los dos
primeros modos de vibrar; para estos dos modos, se seleccionó el cinco por ciento del amortiguamiento
crítico. Además, M se representó como matriz de masa diagonal. Por otro lado, el paso de tiempo
seleccionado para la integración de la ecuación del movimiento fue de 0.001s. En esta etapa del estudio no se
consideró el efecto de la interacción suelo-estructura en los análisis de respuesta de los sistemas estructurales.
Estos se consideraron con base rígida.

Con base en lo anterior, se llevó a cabo el análisis ELSE de los sistemas detallados de manera puramente
traslacional e independiente, en las direcciones x e y, restringiendo los grados de libertad de giro y de
desplazamiento lateral en el sentido transversal a la dirección de análisis, de cada nivel de piso. De esta forma,
se obtuvieron las curvas de capacidad lateral global del sistema, que relacionan la fuerza cortante en la base
contra el desplazamiento lateral de azotea referido al centro de masa V(δ0). En la figura 3 se grafican estas
curvas para cada dirección ortogonal de los sistemas SS6-EDS y SA6-EDS.

Dado que interesa definir el valor de la rigidez inicial y los valores de la fuerza cortante en la base y el
desplazamiento lateral del sistema, cuando éste comienza a experimentar un comportamiento no lineal, a
dichas curvas ELSE se realiza un ajuste bilineal por mínimos cuadrados, igualando las áreas por debajo de la
curva real y la curva de ajuste. Así, los valores obtenidos se presentan en la tabla 3; también, en esta misma
tabla, se presentan los valores totales correspondientes de la masa traslacional , la masa rotacional , y el
periodo fundamental de vibrar T, de los sistemas mencionados en el párrafo anterior (Picazo, 2012). Los
valores de y presentados en la tabla 3 se calcularon sumando únicamente las masas e inercias rotacionales
de todos los niveles de la estructura. Una forma más adecuada de calcular estos valores se da en Picazo
(2012). En dicha referencia la masa equivalente en traslación se puede calcular como

M i i
M * i 1
N
(14)
i 1
i

y el momento polar de inercia equivalente con respecto al eje vertical que pasa por su centro de torsión
“global”

9
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

J i i
J 
* i 1
N
(15)

i 1
i

Figura 3 Curvas de la fuerza cortante en la base vs Desplazamiento lateral de azotea, en la

dirección x e y, respectivamente

Tabla 3 Propiedades dinámicas de los sistemas detallados

Kx Vyx yx Ky Vyy yy m J T

Sistema
(kg/cm) (kg) (cm) (kg/cm) (kg) (cm) (kg*s2/cm) (kg*cm*s2) (s)
SS6-EDS 5.90x104 7.98x105 13.54 5.90x104 8.22x105 13.94 1941.0 1.4191x109 0.76
SA6-EDS 7.50x104 9.82x105 13.09 8.80x104 1.01x106 11.47 2048.7 1.4979x109 0.85

Análogamente, se llevó a cabo el análisis sísmico no lineal paso a paso de cincuenta estructuras simuladas.
Estas estructuras se simularon para tomar en cuenta la incertidumbre asociada a las propiedades geométricas y
mecánicas del sistema, así como las cargas que actúan sobre ellas. De estos análisis fue posible obtener las
curvas de respuesta sísmica de historia-tiempo; en la figura 4 se muestran dichas curvas para uno de los
sistemas simulados, SA6-EDS, para cada dirección ortogonal.

Figura 4 Curvas de la fuerza cortante en la base vs Desplazamiento lateral de azotea

10
 Sociedad Mexicana de Ingeniería Estructural
MODELOS SIMPLIFICADOS DE LOS SISTEMAS DETALLADOS

Los sistemas simplificados de referencia se definen de acuerdo al procedimiento descrito en la subsección

“SISTEMAS SIMPLIFICADOS DE REFERENCIA”, considerando al sistema detallado con propiedades
medias. Para caracterizar las propiedades elástica e inelástica de los elementos columna de este sistema, se
emplean las curvas de capacidad lateral global del sistema detallado en consideración, Vxvsδ0x y Vyvsδ0y, como
las mostradas en la figura 3, así como también sus correspondientes curvas de la fuerza cortante en la base
contra el desplazamiento lateral en el extremo superior, Vivsδi, de cada uno de sus marcos, en las dos
direcciones ortogonales; a partir de esta información y con base en el planteamiento descrito anteriormente, se
determinan las funciones (δ) de los elementos del SSR, para cada dirección x e y.

En la figura 5 se presentan las funciones (δ) de los elementos del SSR obtenidas para el sistema detallado
SA6-EDS. De aquí en adelante se presentan los resultados obtenidos para los SSR correspondientes a los
sistemas detallados SS6 y SA6, diseñados con EDS, los cuales serán designados como sistemas simétrico y
asimétrico de un nivel, “SS1-EDS y SA1-EDS”, respectivamente. Las dimensiones en planta y la altura del
SSR son las mismas que las de su correspondiente sistema detallado. Las propiedades elásticas de los
elementos de los sistemas SS1-EDS y SA1-EDS, en las dos direcciones ortogonales, se obtienen a partir de
efectuar el ajuste bilineal por mínimos cuadrados de las curvas V(δ) (figura 5).

Con base en la rigidez elástica calculada, se determinan las dimensiones de los elementos mediante las
siguientes relaciones:

b  8 144 I y 3 I x (16)

h  8 144 I x 3 I y (17)

Para determinar las propiedades inelásticas de dichos elementos, las curvas (δ) se transforman a los
diagramas equivalentes de momento contra curvatura M vs φ, con base en el procedimiento descrito en el
apéndice.

En la tabla 4 se presentan las propiedades mecánicas obtenidas de los elementos del sistema SA1-EDS. En
ella se muestran las dimensiones de las columnas, b y h, los momentos resistentes de fluencia en cada
dirección ortogonal, Myx y Myy, y las ductilidades de curvatura y resistencia residual Duct1, Duct2 y Rduct
respectivamente, que definen el modelo de degradación de resistencia usado en el programa Ruaumoko (Carr,
2007). Debido a que el programa Ruaumoko únicamente permite definir un modelo de degradación de
resistencia en una dirección dada del elemento, aquí se tomó el criterio de obtener el promedio de las
ductilidades de curvatura que se tienen en ambas direcciones ortogonales de una sección en particular. Este
mismo criterio fue tomado para realizar el análisis sísmico de los sistemas detallados.

Tabla 4 Valores de las propiedades mecánicas de los elementos columna del SA1-EDS

b h Myx Myy
Columna Duct1 Duct2 Rduct
(cm) (cm) (kg*cm) (kg*cm)
C12 100.7 196.8 1.95E+08 7.39E+07 15.0 21.1 0.3
C22 170.8 165.0 1.95E+08 1.80E+08 12.6 23.0 0.3
C00 179.0 172.2 2.49E+08 2.83E+08 11.8 17.6 0.3
C11 112.8 140.2 1.08E+08 7.39E+07 14.7 17.6 0.3
C21 191.3 117.6 1.08E+08 1.80E+08 12.4 19.5 0.3

FUNCIONES DE CONFIABILIDAD

Para determinar el índice de confiabilidad o índice “β” propuesto por Cornell (1969), de los sistemas
detallados y simplificados, se aplica lo descrito en la sección “ PLANTEAMIENTO”. De acuerdo con ello,
el índice de reducción de rigidez secante IRRS de los sistemas simulados en las dos direcciones ortogonales se

11
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

calcula a partir de los resultados obtenidos del análisis ELSE y el análisis dinámico no lineal paso a paso,
aplicando las ecuaciones 2a y 2b. Por otro lado, el índice de Cornell se calcula aplicando la ecuación 4, en
donde la medida de la intensidad sísmica normalizada se toma igual a q (ecuación 11).

A continuación se presentan las curvas de vulnerabilidad sísmica de los sistemas detallados y las
correspondientes curvas de vulnerabilidad de los sistemas simplificados, para los casos SS1-EDS, SS6-EDS,
SA1-EDS y SA6-EDS.

Figura 5 Funciones V ( ) de los elementos del SA1-EDS para ambas direcciones ortogonales,
obtenidas a partir del sistema SA6-EDS

12
 Sociedad Mexicana de Ingeniería Estructural

Figura 6 Funciones de confiabilidad sísmica de los sistemas detallados y simplificados, en

términos de la intensidad normalizada q

De la figura 6 se puede observar que al comparar las funciones de confiabilidad de los sistemas SS1-EDS y
SS6-EDS, el índice de confiabilidad β para una misma intensidad normalizada es mayor para el primer
sistema que para el segundo mencionado; similarmente, para los sistemas SA1-EDS y SA6-EDS se tiene la
misma tendencia, pero ahora la diferencia entre los índices β para un mismo valor de la intensidad
normalizada en ambos sistemas, es significativamente mayor de lo que se tenía para el caso del sistema
simétrico.

CONCLUSIONES

En el caso de los sistemas detallados presentados en este trabajo (casos EDS), ambos sistemas, simétrico y
asimétrico, fueron diseñados de tal forma que cumplieran con los criterios de diseño establecidos en el
RCDF2004. Esto llevó a diseños que no son completamente comparables en el sentido de no mantener iguales
o parecidos ciertos valores de los parámetros que definen al sistema (p.ej. periodo de vibración, resistencia
lateral, capacidad lateral de deformación, etc.). En este caso el sistema SA6-EDS, con propiedades medias,
presentó valores del periodo de vibración y de la capacidad lateral resistente y de deformación mayores que
las del caso simétrico. Aun considerando estas diferencias se encuentra que la confiabilidad del sistema
simétrico (SS1-EDS) es menor que la del sistema asimétrico (SA1-EDS).

Las diferencias que se tienen en el valor del índice β entre los modelos detallados y simplificados se deben
principalmente a las simplificaciones que se realizaron al caracterizar los elementos del modelo SSR, como
son el estimar su rigidez inicial a partir de trazar una secante a la rama ascendente a la curva de capacidad
lateral del sistema, lo cual modificó los periodos de vibrar; además, no se consideró en el sistema simplificado
un modelo histerético que representara, de manera equivalente, las características de degradación de rigidez y
resistencia de su respectivo sistema detallado, a partir de estimar, con suficiente precisión, la formación de los
lazos histeréticos obtenidos de dichos sistemas, al excitar su base ante un movimiento sísmico. Para el caso de
los sistemas asimétricos, estas diferencias del índice β son mayores, ya que dentro del comportamiento no
lineal de un sistema complejo asimétrico su respuesta sísmica puede variar más fuertemente con respecto a
aquélla que se tendría en un sistema simplificado de un solo nivel, debido a la influencia que tiene el
fenómeno de torsión a lo largo de la altura del sistema complejo, principalmente en los elementos de la
periferia, lo cual ocurre de forma diferente en un sistema simple. Todo lo anterior puede explicar porque el
índice β para un mismo nivel de intensidad q, es menor para el sistema SS1-EDS, que aquél correspondiente
para el sistema SA1-EDS.

13
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

El estudio de los sistemas SA6-EDSR y SA6-EDSA y otros sistemas con distintos niveles y configuraciones
estructurales permitirá entender mejor los resultados presentados aquí sobre la influencia de la asimetría en la
respuesta de las estructuras. Estos mismos resultados podrán utilizarse para poder establecer factores de
reducción de ordenadas espectrales que deban aplicarse para el diseño sísmico de los sistemas asimétricos a
fin de lograr los mismos niveles de confiabilidad que los sistemas simétricos.

Dados los resultados que se obtienen con los SSR, es necesario realizar más estudios sobre la forma de definir
tales sistemas, de tal forma que puedan representar en forma más adecuada el comportamiento de los sistemas
detallados simétricos o asimétricos.

REFERENCIAS

Bresler B, (1960), “Design Criteria for Reinforced Columns under Axial Load and Biaxial Bending”, J.
ACI, Vol. 32, No. 5, pp 481-490.

Carr A J, (2007), “Ruaumoko 3D Users Manual”, University of Canterbury, Christchurch, New Zealand.

Cornell C A, (1969), “A probability based structural code”, Journal of the American Concrete Institute. 66
(12).

Corona G, (2007), “Programa ECOgcW de análisis y diseño de edificios de concreto”, g.c. ingeniería y
diseño, s.c.

Díaz-López O, and Esteva L, (2006), “Seismic reliability analysis of complex nonlinear systems, using
secant-stiffness-reduction indicators”, Institute of Engineering, National University of Mexico.

Díaz-López O, Ismael E, Esteva L, (2008), “About efficient algorithms for the determination of seismic
reliability functions of multistory buildings”, IFIP WG7.5 Reliability and optimization of structural
systems, Toluca, México.

Esteva, L, Díaz, O, García, J, (2010) “Estimación de funciones de vulnerabilidad y riesgo ante sismos en
construcciones en el Distrito Federal: Manual de lineamientos para identificar sistemas que ameritan
evaluación, para determinar sus niveles de vulnerabilidad y riesgo y para decidir sobre acciones
pertinentes de demolición o rehabilitación”, Informe de Investigación elaborado para la Secretaría de
Protección Civil del Gobierno del Distrito Federal, Instituto de Ingeniería, UNAM.

Esteva L and Ismael E, (2004), “A maximum likelihood approach to system reliability with respect to
seismic collapse”, Proc. 11th IFIP WG7.5 Working Conference, Banff, Canada.

Parme A L, Nieves J M, Gouwens A J, (1966), “Capacity of Reinforced Rectangular Columns Subject to

Biaxial Bending”, Journal of the American Concrete Institute. Vol. 63, pp 912-923.

Paulay T, Priestley M, (1992), “Seismic design of reinforced concrete and masonry buildings”, John Wiley
& Sons, USA.

Picazo Y, (2012), “Modelo simplificado para estimar las funciones de vulnerabilidad y riesgo sísmico de
edificios asimétricos en planta”, Tesis de doctorado (en proceso), División de Posgrado, Universidad
Nacional Autónoma de México.

RCDF (2004), “Reglamento de Construcciones del Distrito Federal”, Gaceta Oficial del Distrito Federal,
Mexico.

14
 Sociedad Mexicana de Ingeniería Estructural
APÉNDICE

DETERMINACIÓN DE LAS PROPIEDADES MECÁNICAS DE LAS ARTICULACIONES

PLÁSTICAS CORRESPONDIENTES A LOS ELEMENTOS DEL SSR

Para modelar los elementos columna del SSR se utilizó el modelo viga de un componente de Giberson que se
emplea en el programa Ruaumoko, el cual está formado por un elemento elástico y por una zona de
articulación plástica en los extremos del mismo (figura 7). Para este modelo viga sujeto a doble curvatura en
sus extremos, se consideró una longitud de articulación plástica definida de la siguiente manera (Paulay &
Priestley, 1992):

L p  0.08  L 2  0.022 f yl  d bl (A-1)

Donde:

L p Longitud de articulación plástica

L Longitud del elemento
f yl Esfuerzo de fluencia del acero de refuerzo longitudinal
d bl Diámetro del acero de refuerzo longitudinal

Figura 7 Elemento viga de un componente de Giberson

Dado que se desea establecer un modelo bilineal de degradación de resistencia en función de la demanda de
ductilidad para caracterizar las zonas de articulación plástica de las columnas del SSR de acuerdo con el
programa Ruaumoko, se emplean las curvas de capacidad lateral de la fuerza cortante en la base contra el
desplazamiento lateral de azotea V vs δ de los elementos del sistema simplificado en cuestión, y a partir de
estas curvas se obtienen sus respectivos diagramas equivalentes del momento resistente contra curvatura M vs
φ para cada una de las columnas de dicho sistema. Lo anterior se logra a partir de realizar un ajuste bilineal
por mínimos cuadrados en dichos diagramas y con base en las relaciones que se mencionan a continuación:

M I  VI  h / 2 (A-2)

t   y   p (A-3)

 p   p / Lp (A-4)

 I  T   E (A-5)

 y  M ycol / EI (A-6)

En las expresiones anteriores

15
XVIII Congreso Nacional de Ingeniería Estructural Acapulco, Guerrero 2012

MI, VI son el momento resistente en la base y la fuerza cortante resistente en el elemento del SSR, dentro del
intervalo inelástico, respectivamente.

φp, Lp son la curvatura y la rotación plástica en los extremos de las columnas del SSR, respectivamente. Ésta
última, se define por la relación dada entre el desplazamiento lateral inelástico en el elemento del SSR y la
altura total del sistema detallado, θp= I/h.

φy, φt son las curvaturas de fluencia y total en los extremos de las columnas del SSR, respectivamente.

δT, δE son los desplazamientos laterales total y elástico en el elemento del SSR, respectivamente; aquí,
δE=V/K.

Para determinar la curvatura de fluencia de las columnas de esquina correspondientes a un extremo y a una
dirección dada en el SSR, por ejemplo, las columnas C11 y C12 en la dirección x (ver figura 5), se toma el
momento de inercia Iy con respecto a los ejes locales de uno de dichos elementos, y el momento de fluencia de
una columna se calcula como Mycol=My/2; donde My=(Vy*h)/2. La fuerza cortante de fluencia Vy se define a
partir de realizar el ajuste bilineal de la correspondiente curva V vs δ para el marco correspondiente. De esta
manera, se calcula el valor de la curvatura de fluencia para las columnas mencionadas anteriormente. Cabe
mencionar aquí, que para el elemento central del SSR (la columna C00), la curvatura de fluencia para una
dirección dada se obtiene de manera similar, pero ahora, el momento de fluencia es igual a Mycol=My.
Finalmente, los diagramas Mvsφ de cada elemento del SSR se normalizan con sus respectivos valores del
momento y la curvatura de fluencia, es decir: MI/Mycol y φt/φy respectivamente. Con base en lo anterior, se
obtiene el modelo bilineal de degradación de resistencia que se emplea en el programa Ruaumoko “Factor de
Carga vs Ductilidad de Curvatura”.

16