Componente Vertical. DR Esteva
Componente Vertical. DR Esteva
Componente Vertical. DR Esteva
72 45-79 (2005)
RESUMEN
Se muestran resultados de un estudio sobre las características del componente sísmico vertical de una
pequeña muestra de registros que presentan las mayores intensidades verticales reportadas en la Base
Mexicana de Datos de Sismos Fuertes entre 1960 y 1999. Además, se muestran resultados de análisis no
lineales de vigas y marcos planos, diseñados convencionalmente, incluyendo y excluyendo el componente
sísmico vertical en los análisis dinámicos paso-a-paso. Se sintetiza un resumen de estudios previos del
movimiento vertical del terreno y de su efecto en las estructuras. Asimismo, se presenta un resumen
condensado de las consideraciones de reglamentos y normas en el mundo para tomar en cuenta los efectos
tridireccionales del movimiento del terreno. Finalmente, se presentan algunas reflexiones sobre el diseño
sísmico incluyendo el movimiento dinámico vertical, y se dan algunas sugerencias cualitativas para
entender mejor sus efectos en estructuras, particularmente en edificios de concreto reforzado.
ABSTRACT
An analysis is presented about the vertical seismic component in some ground motion records in Mexico.
These records are associated with the largest vertical intensities recorded in Mexico and reported in
Mexican Database of Ground Motion Records (Base Mexicana de Datos de Sismos Fuertes) between
1960 and 1999. Some results are shown of the nonlinear response analysis of beams and plane frames,
which were conventionally designed. The cases when the vertical seismic component is included or
ignored in the step-by-step dynamic analysis are both covered. Previous studies about vertical ground
motion and its effect on structures are summarized. Recommended provisions to include the vertical
ground motion in the seismic design of structures available in some international building codes are
discussed. Finally, some thoughts are presented about seismic design including vertical dynamic ground
motion, together with some qualitative suggestions to better understand its effects on structures,
particularly on reinforced concrete buildings.
Artículo recibido el 7 de octubre de 2004 y aprobado para su publicación el 11 de marzo de 2005. Se aceptarán comentarios y/o
discusiones hasta cinco meses después de su publicación.
(1) Departamento de Materiales, Universidad Autónoma Metropolitana Azcapotzalco, Av. San Pablo # 180, Col. Reynosa
Tamaulipas, México, DF, 02200, tperea@correo.azc.uam.mx
(2) Instituto de Ingeniería, UNAM, Ciudad Universitaria, 04510 México, DF., LEstevaM@iingen.unam.mx
45
Tiziano Perea y Luis Esteva
INTRODUCCIÓN
El componente vertical de movimientos fuertes del terreno ha sido, durante muchos años, una de
las variables menos estudiada en ingeniería sísmica. Eso se refleja en las recomendaciones
simples de algunos códigos, normas o reglamentos de construcción para la consideración del
componente sísmico vertical en la respuesta dinámica de las estructuras; asimismo lo refleja la
nula o somera documentación en el tema que presentan los comentarios de algunos reglamentos.
Lo anterior hace que el diseñador, en ocasiones, llegue a interpretaciones erróneas, y
consecuentemente, a que el efecto del componente sísmico vertical sea mal entendido o ignorado.
En este sentido es necesario estudiar, por una parte, las características del componente vertical de
los movimientos fuertes del terreno, y por otra, su efecto en las obras civiles. Con ello, será
posible tener propuestas reglamentarias transparentes y compatibles con el peligro sísmico para
considerar el componente vertical de los movimientos fuertes del terreno en la respuesta dinámica
de las estructuras.
ANTECEDENTES
Los estudios sobre movimientos del terreno relacionados con el componente sísmico vertical se
han enfocado, en general, a obtener relaciones de las intensidades sísmicas (aceleración,
desplazamiento, velocidad máximos del terreno y/o espectrales) del componente vertical entre el
horizontal (o viceversa), ya sea tomando en cuenta o ignorando su relación con otras variables
como la magnitud, la distancia epicentral, la frecuencia o periodos, entre otras. En general se
intenta predecir la intensidad vertical a través de la correspondiente intensidad horizontal.
Newmark et al. (1973) realizan uno de los primeros estudios sobre el componente vertical
de movimientos fuertes del terreno. En dicho trabajo se obtienen cocientes de las aceleraciones
máximas del terreno del componente vertical entre el horizontal de 33 registros de Estados
Unidos. El estudio señala que, en promedio, la aceleración vertical es del orden de 2/3 del
horizontal; dicho cociente ha sido adoptado en muchos códigos de diseño antisísmico en el
mundo ante la ausencia de otra información.
Otro estudio precursor del componente sísmico vertical lo realizó Rosenblueth (1975).
Este trabajo también relacionó la intensidad sísmica vertical con la horizontal, aunque fue el
primero que identificó que dicha relación es función del periodo. En este trabajo se propone un
espectro elástico para componente vertical similar al horizontal, con modificación en la escala de
ordenadas espectrales y reducción en 70% de la escala de periodos; adicionalmente propone
obtener el espectro inelástico a partir del cociente del espectro elástico entre un factor por
ductilidad vertical de tres, aunque no presenta ninguna justificación para tomar dicho factor.
Con esta metodología, aunque con una diferente base de registros, otros autores han
propuesto diferentes valores de dicho cociente aplicable a otros sitios (Kawashima et al., 1985;
Ambrasays y Simpson, 1995; Boomer y Martínez, 1996; Mohammadioun, 1996). En la tabla 1 se
muestran algunos valores del cociente V/H propuestos por algunos autores.
46
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Algunos estudios han destacado que el cociente de aceleraciones máximas del terreno
horizontal a vertical es altamente dependiente de la distancia epicentral: se observan valores altos
para distancias cortas y bajos para grandes distancias (Singh, 1995; Hudson et al., 1996;
Bozorgnia, Niazi y Campbell, 1996; Perea y Esteva, 1999). Otros más enfatizan en sus
resultados que dicho cociente es además altamente dependiente del periodo vertical, con valores
altos en periodos cortos y valores bajos en periodos altos (Singh, 1995; Hudson et al., 1996;
Ohno et al., 1996; Kusunoki, 1998; Bozorgnia, Campbell y Niazi, 2000; Perea y Esteva, 2003a).
Sharma (2000), Bozorgnia, Campbell y Niazi (2000), y Perea y Esteva (1999) han
obtenido espectros de respuesta verticales a partir de leyes de atenuación espectrales con base en
el análisis de registros sísmicos. Los dos primeros trabajos mencionados obtienen, a partir de las
expresiones de atenuación, cocientes espectrales de la intensidad vertical entre la horizontal. Los
resultados del segundo grupo de autores muestran una gran dependencia entre los componentes
horizontal y vertical con el periodo de oscilación, la distancia de la fuente al sitio, y las
condiciones locales del sitio, no así con la magnitud y el tipo de falla. Los resultados del último
grupo de autores (Perea y Esteva, 1999) muestran que la magnitud (M) del sismo y la distancia
entre el epicentro y el sitio (R) mantienen una relación similar con la amplitud espectral de ambos
componentes (β2 y β3), el horizontal y el vertical, y que son los efectos de sitio (β1) los que
principalmente controlan las diferencias en las amplitudes espectrales en dichas direcciones
(figura 1). Lo anterior obedece a que ambos componentes, el horizontal y el vertical, comparten
las mismas características de la fuente y los efectos del trayecto.
47
Tiziano Perea y Luis Esteva
Saragoni y Hart (1974), Perea, Alamilla y Esteva (2001) han propuesto generar
acelerogramas del componente vertical. Los primeros propusieron generarlos mediante el
procedimiento de pasar ruido blanco a través de un filtro que representa segmentos consecutivos
de un registro de aceleraciones del terreno por diferentes procesos estacionarios uniformemente
modulados; dichos acelerogramas simulados o generados tienen discontinuidades muestrarias en
las intensidades donde se cambia de un proceso a otro. Como alternativa, los mismos autores
proponen representar a las aceleraciones del terreno como una serie modulada de impulsos con
magnitudes aleatorias. El segundo grupo de autores genera acelerogramas del componente
vertical en suelo firme a partir de un modelo de funciones moduladoras de amplitud y frecuencia
(funciones de atenuación generalizadas). Cada historia de aceleraciones se considera como una
realización de un proceso estocástico gaussiano no estacionario, con parámetros estadísticos
dependientes de la magnitud y la distancia al sitio; para ligar estos parámetros del movimiento del
terreno con su correspondiente magnitud y distancia al sitio, se determinan funciones semi-
empíricas llamadas funciones de atenuación generalizadas. El conjunto de historias detalladas del
movimiento del terreno empleadas para obtener dichas funciones corresponden a movimientos
generados en diversas fuentes de la costa del Pacífico de México y registrados en diferentes sitios
en la vecindad de la costa. En los resultados se aprecia una dispersión significativa que
caracteriza al modelo adoptado para representar el movimiento del terreno como un proceso
gaussiano no estacionario con densidad espectral evolutiva. Esta dispersión está asociada con las
características de los acelerogramas reales incluidos en la muestra de estudio y al intervalo de los
parámetros de la fuente sísmica correspondientes a los registros.
βi
3
Componente horizontal β1 β2 β3
2 Componente vertical β1 β2 β3
0
0.0 0.2 0.4 0.6 0.8 1.0
-1
-2
-3
T (s)
48
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Comportamiento estructural
Los estudios analíticos sobre los efectos del componente vertical del movimiento sísmico en el
comportamiento estructural han consistido, fundamentalmente, en investigar la respuesta
dinámica elástica e inelástica de diferentes sistemas estructurales sujetos a la acción del
componente sísmico vertical, combinándola con los efectos del componente horizontal y de las
cargas de gravedad.
En uno de los primero trabajos analíticos (Rosenblueth, 1975) se propuso que, para fines
prácticos, se pueden considerar los efectos del componente vertical del movimiento aplicando
factores de carga para las fuerzas gravitatorias del orden de 0.8 y 1.2; se menciona que en zonas
donde la probabilidad de que ocurra un temblor intenso de epicentro cercano es pequeña, puede
omitirse el análisis sísmico por componente vertical.
49
Tiziano Perea y Luis Esteva
REFERENCIA ξ COMENTARIOS
Jennings y Husid 2% Los resultados no muestran una influencia significativa del componente sísmico vertical.
(1968)*
Iyengar y Shinozuka 1% En voladizos, se alteran considerablemente las deflexiones laterales en el extremo libre, y las fuerzas cortantes
(1972)* y momentos flectores en el extremo restringido.
Anderson y Bertero Incremento en las demandas de ductilidad local de las vigas de los pisos superiores, de las columnas de los
(1973) pisos inferiores, y de las columnas de pisos superiores.
Cheng y Ester Incremento de desplazamientos laterales en estructuras con frecuencia fundamental baja o cercana a la mitad
(1974) de la frecuencia vertical asociada con la máxima amplitud del espectro de Fourier.
Rosenblueth Los análisis realizados mostraron la aparición de aceleraciones verticales en los pisos de edificios altos, las
(1975) cuales son funciones moderadamente crecientes de la altura. Son máximas las aceleraciones en edificios de
altura intermedia y resultan poco mayores en estructuras de concreto que en metálicas.
Iyengar y Saha 1% Altera las deflexiones laterales y las envolventes de momentos y cortantes máximos en voladizos, con
(1977)* aumentos en unos casos y reducciones en otros.
Goel Incremento considerable de las demandas de ductilidad para columnas y en menor grado para vigas y
(1977) contravientos; los desplazamientos laterales permanecen relativamente inalterados.
López y Chopra Incremento de las demandas de ductilidad locales de la estructura; los desplazamientos laterales no son
(1978) afectados significativamente.
Tani y Soda Incrementos en carga axial en columnas hasta del 50%; los desplazamientos laterales máximos, cortantes en
(1980) columnas, y otras respuestas, aumentaron o disminuyeron en alrededor del 20%.
Kikuchi y Yoshimura 5% Amplificación importante de las fuerzas axiales en las columnas interiores de los pisos inferiores.
(1984)
Sovero 5% Los valores del factor de reducción por ductilidad vertical pueden oscilar entre 1 y 2.
(1986)
Papaleontiou y Roesset Para marcos de acero de 4, 10, 16 y 20 niveles, los incrementos en carga axial en columnas de entrepisos
(1993) superiores son del 72, 85, 76 y 64%; para columnas del primer entrepiso son del 56, 52, 24 y 21%.
Papazoglou y Elnashai La capacidad a momento y ductilidad en columnas de concreto reforzado se reduce debido a fluctuaciones en
(1996) la carga axial.
Yamanouchi y Hasegawa Existe influencia en las respuestas locales, pero ninguna en las respuestas globales de edificios.
(1996)
Esteva No influye significativamente en la respuesta global, sólo en la respuesta local de los miembros sometidos
(1996) principalmente a flexión.
Elnashai 2% La aportación del componente vertical en columnas es del orden de 25% al 40% del total de la fuerza axial
(1997)
Reyes-Salazar y Haldar Los incrementos de momentos en columnas son despreciables, no así con la carga axial, la cual se incrementa
(2000) considerablemente.
Kehoe y Attalla Las fuerzas inducidas en los elementos estructurales son mucho menores que las correspondientes a los
(2000) efectos de la carga de gravedad, excepto en estructuras de concreto de gran claro.
Alaghebandian, Otani y Los desplazamientos y cortantes laterales de piso son ligeramente afectados por el componente vertical del
Shiohara (2000) movimiento, pero se tienen fuertes variaciones en la carga axial en columnas.
50
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Destaca de la tabla 3 que, a la fecha, únicamente los reglamentos del Eurocódigo (CEN,
1994) y de Grecia (GSC, 1992) han propuesto un coeficiente sísmico dependiente del periodo de
vibración vertical para estimar la demanda sísmica vertical. Adicionalmente, el código de diseño
sismorresistente de Grecia (GSC, 1992) refiere al uso de la formula de Rayleigh para estimar el
periodo fundamental de vibración vertical, y a partir de él, estimar la demanda correspondiente.
Por otro lado, llama la atención que algunos reglamentos de diseño sismorresistente no
presentan ninguna consideración para tomar en cuenta el componente vertical de movimientos
del terreno para el diseño estructural, sobre todo en los países con un reconocido nivel de
desarrollo en el conocimiento de su sismicidad y riesgo sísmico; tal es el caso de Nueva Zelanda
(NZS-4203,1992), Japón (BSLEO, 1981; SSCECJ, 1980), Chile (INN, 1989) y Canadá (NBCC,
1990), todos ellos países donde para ciertas regiones es alta la probabilidad de que el componente
vertical puede alcanzar intensidades considerables. Las justificaciones para no tomar en cuenta el
componente sísmico vertical para el diseño estructural en estos países son desconocidas por los
autores, aunque no descartan que: (1) el efecto del componente sísmico vertical en las estructuras
con registros locales no ha sido estudiado, o bien, no lo suficiente para incluirlo en el reglamento
o para adaptarlo con propuestas de otros reglamentos; (2) el efecto de dicho componente en las
construcciones del sitio es pequeño, o nulo, que despreciarlo no implica ningún riesgo
51
Tiziano Perea y Luis Esteva
estructural; (3) los factores de sobrerresistencia implícitos en el diseño son suficientes para
considerar las incertidumbres relativas al componente sísmico vertical, y por lo tanto, despreciar
este componente no implica ningún riesgo estructural.
52
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
b FV = 2 ⋅ C H ⋅ I ⋅ (CM + CV )
Voladizos
(CM + CV) ⋅ C V ⋅ m i u i
Grecia GSC, 1992 b FVi = n
No Aplica
∑ mi u i
i =1
n
⎛ n
⎞ α ⋅ I ⋅ β(TV )
T = 2π ∑m u
i =1
i i
2
⎜⎜ g ∑ m i u i ⎟⎟
⎝ i =1 ⎠
Cv = 0.7
Rµ V
n
2 ⋅ CH ⋅ I
Irán IC-SRDB, 1988 b FV = (CM + CV) No Aplica
Rµ
53
Tiziano Perea y Luis Esteva
(CM i + CVi ) ⋅ H i
China GBJ-11, 1989 b FVi = FV ⋅ No Aplica
(CM + CV) ⋅ H
Edificios altos, Chimeneas, etc. Cubiertas L>24m, Volados, etc.
54
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
de tener un número menor de grados de libertad, consideren en el modelo masa puntual actuando
en un sentido horizontal. Algunos programas de análisis tienen la opción de considerar la
contribución de la masa rotacional y vertical, y modelar con masa discreta que simule el caso real
de masa distribuida. El discretizar la masa y considerar su influencia horizontal, vertical y
rotacional, implicará tener un mayor número de grados de libertad en el modelo.
Por todo lo anterior, los códigos que consideran el componente sísmico vertical presentan
recomendaciones y prácticas que, en ocasiones, se prestan a interpretaciones erróneas, lo cual
también origina que en ocasiones el efecto de tal componente sea mal entendido o ignorado.
Sobrestimar los efectos del componente vertical no necesariamente significará alcanzar en un
grado de seguridad estructural más alto, sobre todo si se piensa en un posible cambio de
mecanismo de colapso en caso de diseñar vigas más fuertes sin revisar su interacción con las
columnas. Por ello, hace falta tener propuestas reglamentarias transparentes para considerar el
componente vertical del movimiento y estudiar en el dominio de la frecuencia el efecto de diseñar
considerando el periodo vertical de la estructura, con sus combinaciones entre los demás
componentes sísmicos. Dichos análisis podrán ser comparados con análisis dinámicos paso-a-
paso que utilicen acelerogramas reales o simulados de movimientos representativos e
independientes entre sí, cuyas características de intensidad, duración y contenidos de frecuencia
sean compatibles con el peligro sísmico del sitio de interés.
Cocientes de las aceleraciones verticales máximas del terreno entre las horizontales
Para discutir la influencia del periodo de oscilación, se estudiará a continuación una muestra
específica de diez registros mexicanos en la que, según los autores, bajo ciertas circunstancias
existe una gran influencia del componente sísmico vertical debido a su cercanía con la fuente
generadora. En la figura 3 se muestra la ubicación de las estaciones asociadas a la muestra de los
registros seleccionados y enlistados en la tabla 4. En la figura 4 se muestran los respectivos
espectros de seudoaceleración (ξ=5%) y en la figura 5 se muestran los cocientes de las ordenadas
espectrales del componente vertical entre las correspondientes al componente horizontal; La
clave del registro corresponde a la que define la Base Mexicana de Datos de Sismos Fuertes.
55
Tiziano Perea y Luis Esteva
R (km)
Figura 2. Variación del cociente V/H de registros de la costa del Pacífico de México
Los cocientes espectrales de la respuesta vertical entre la horizontal que se calculan para
esta base de registros (figura 6.a) muestran dispersión considerable; sin embargo, se distingue que
la mayor dispersión se presenta en el intervalo de periodos bajos, donde además se observa que el
componente vertical incluso puede ser más intenso. En cambio, para los periodos mayores de
0.5s el componente sísmico vertical no supera el 40% de la ordenada espectral del componente
horizontal, en contradicción con el 2/3 que algunos reglamentos proponen con base en el cálculo
del cociente para periodo cero.
56
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Sa (g) Sa (g)
5.0 5.0
4.0 4.0
3.0 3.0
2.0 2.0
1.0 1.0
0.0 0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0
T (s) T (s)
a) Componente horizontal b) Componente Vertical
Figura 4. Espectros de respuesta de seudoaceleración (ξ=5%,)
57
Tiziano Perea y Luis Esteva
3.0 2.0
1.5
0.5 2.0
1.0
1.0
0.5
2.0
1.0 1.0
1.5
1.0
0.5 0.5
0.5
1.5 1.5
1.0
1.0 1.0
0.5
0.5 0.5
58
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
V/H V/H
6.0 3
5.0
4.0 2
3.0
2.0 1
1.0
0.0 0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0
T (s) T (s)
a) Cociente vertical a horizontal de registros b) µ ± 1σ y envolvente
Figura 6. Variación de la relación espectral del componente vertical a horizontal
En la figura 7, se comparan los espectros de respuesta (ξ=5%) del componente vertical de los
registros de campo libre con los estimados del escalamiento del espectro horizontal por el factor
2/3 (RCEG-90 y MDOC-93) y por la envolvente del promedio propuesta. Como se observa en
dichas gráficas, los espectros verticales obtenidos a partir del producto del espectro horizontal por
un factor de escala dependiente del periodo (envolvente) se aproximan razonablemente a los
obtenidos a partir de los registros, excepto en algunos casos para periodo corto, donde se
distinguen las mayores desviaciones. El factor de escala constante 2/3 no conduce a predicciones
correctas, particularmente para periodos verticales altos.
En dicha gráfica, se pone en evidencia que las propuestas reglamentarias vigentes para el
estado de Guerrero y zona D del MDOC son conservadoras para estructuras diseñadas con
periodos verticales mayores que 0.2s; no obstante, un gran número de estructuras puede tener en
la dirección vertical un periodo menor que 0.2s, y por tanto coincidir con el intervalo para el cual
la propuesta reglamentaria resulta ser insuficiente.
59
Tiziano Perea y Luis Esteva
Sa (g) Sa (g)
5.0 VICS800609 4.0 VCPS870207
4.5 3.5
RCEG-90, MDOC-93 RCEG-90, MDOC-93
4.0
3.0
3.5 Envolvente Envolvente
3.0 2.5
2.5 2.0
2.0 1.5
1.5
1.0
1.0
0.5 0.5
0.0 0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0
T (s) T (s)
Sa (g) Sa (g)
0.7 5.0 IAGS791015
ZACA850919
0.6 RCEG-90, MDOC-93 RCEG-90, MDOC-93
4.0
0.5 Envolvente Envolvente
3.0
0.4
0.3 2.0
0.2
1.0
0.1
0.0 0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.2 0.4 0.6 0.8 1.0
T (s) T (s)
Sa (g) Sa (g)
0.5 CALE850919 1.4 CALE970111
0.3 0.8
0.2 0.6
0.4
0.1
0.2
0.0 0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0
T (s) T (s)
Figura 7. Comparación del espectro vertical del registro de campo libre con los estimados
60
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Sa (g) Sa (g)
0.8 COPL931025 1.2 BALC941210
0.4 0.6
0.4
0.2
0.2
0.0 0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.2 0.4 0.6 0.8 1.0
T (s) T (s)
Sa (g) Sa (g)
0.6 ACAC890425 1.6 RIXC951021
1.4
0.5 RCEG-90, MDOC-93 RCEG-90, MDOC-93
1.2
Envolvente Envolvente
0.4
1.0
0.3 0.8
0.6
0.2
0.4
0.1
0.2
0.0 0.0
0.0 0.5 1.0 1.5 2.0 0.0 0.2 0.4 0.6 0.8 1.0
T (s) T (s)
Figura 7. Comparación del espectro vertical de campo libre con los estimados (cont.)
Sa (g) 0.8
RCEG-90, MDOC-93
0.7
Promedio de registros
0.6
Envolvente de registros
0.5
0.4
0.3
0.2
0.1
0.0 T (s)
0.0 0.5 1.0 1.5 2.0
Figura 8. Comparación del espectro de diseño vertical con el estimado
61
Tiziano Perea y Luis Esteva
Las expresiones generales para calcular los periodos de vibración vertical (TVN) de vigas con
masa distribuida y sus correspondientes formas modales φ(x), que se derivan a partir de la teoría
clásica de vibración de vigas, son reportadas por Young-Budynas (2002) y Chopra (1995) con
base en los resultados presentados por Huang en 1964.
2π w ⋅ L4
TVN = (2)
KN g ⋅ EI
2π w ⋅ L
β2 = (4)
T g ⋅ EI
Tabla 5. Constante KN para los primeros cinco modos de vibración vertical de vigas con
carga uniforme y distintas condiciones de apoyo.
Modo KN VE VA VV
1 K1 22.4 9.87 3.516
2 K2 61.7 39.48 22.03
3 K3 121 88.83 61.70
4 K4 200 157.91 120.90
5 K5 299 246.74 199.86
Para el caso de la viga simplemente apoyada (VA) con carga uniforme (w), la solución
cerrada para el cálculo de los periodos de vibración vertical y de las formas modales
correspondientes (Chopra, 1995), es:
62
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
2π wL4
TVN = (5)
( Nπ) 2 g ⋅ EI
Nπx
φ( x ) = sen (6)
L
El periodo fundamental en la dirección vertical también puede ser estimado con bastante
precisión y en forma sencilla con un modelo de masa discreta de acuerdo con la formula de
Rayleigh (Chopra, 1995):
n (7)
∑m u
2
i i
2π u max
T = 2π i =1
n
≈
1 .2 g
g∑ mi u i
i =1
Se modelan vigas de concreto reforzado con masa discreta en puntos espaciados a cada L/10
sobre la viga, amortiguamiento ξ=5% del crítico, y tres condiciones simples de apoyo:
doblemente empotrado (VE), simplemente apoyado (VA) y en voladizo (VV). Las vigas se
diseñaron para resistir su peso propio y una carga uniforme de 4 T/m, calibrando el claro de la
viga para tener en todos los casos un periodo fundamental vertical de 0.1s. La longitud de los
claros calculados fueron de 10 m, 6.65 m y 4 m, respectivamente.
En la tabla 6 se presentan los periodos (Tv) asociados a los modos de vibrar en la dirección
vertical obtenidos con el Drain-2DX, así como sus correspondientes factores de participación (Γ)
y porcentajes de masa modal efectiva (m*).
63
Tiziano Perea y Luis Esteva
Las rotaciones plásticas más altas se presentaron con los registros más intensos,
es decir, con los de mayor desplazamiento y seudoaceleración espectral
correspondiente al periodo fundamental vertical (0.1s); dichos registros
corresponden a las estaciones con la distancia epicentral y profundidad focal
más cercana a la fuente, como es el caso de los registros del valle de Mexicali.
Para los tres tipos de vigas estudiadas, es importante notar que al determinar los
porcentajes de refuerzo en cada sección crítica en donde pueda desarrollarse
comportamiento no lineal significativo, deben lograrse factores de seguridad
similares entre las capacidades en flexión y las correspondientes demandas para
momento positivo y negativo en cada sección. Esto tiene por objeto evitar la
ocurrencia de comportamiento asimétrico en la articulación plástica que se
forme, así como de los consiguientes efectos en las demandas excesivas de
ductilidad local y en su acumulación ante la acción de temblores excesivos.
64
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Tabla 7. Rotaciones plásticas obtenidas de los análisis con el componente sísmico vertical de
los registros seleccionados (TV=0.1s)
REGISTRO M D TIPO DE SUELO Y SITIO Sd Sa VEe VEc VAc VVe
(km) (cm) (g) R1 R2 R3 R4
VICS800609 6.1 15 Sedimento, Mexicali, Baja California 1.09 4.30 0.075 0.066 0.082 0.075
VCPS870207 5.4 8 Roca, Mexicali, Baja California 0.42 1.66 0.017 0.016 0.022 0.017
IAGS791015 6.6 10 Sedimento, Mexicali, Baja California 0.13 0.50 0.008 0.007 0.004 0.016
ZACA850919 8.1 86 Arcilla, Zacatula, Michoacán 0.13 0.51 0.003 0.003 0.005 0.002
CALE970111 6.9 40 Roca, Caleta de Campos, Michoacán 0.12 0.46 0.002 0.002 0.003 0.001
RIXC951021 6.5 112 Calizas, Tuxtla Gutiérrez, Chiapas 0.10 0.33 0.000 0.000 0.000 0.000
BALC941210 6.3 43 Roca, El Balcón, Guerrero 0.08 0.31 0.000 0.000 0.000 0.000
COPL931025 6.6 20 Roca, Copala, Guerrero 0.06 0.25 0.000 0.000 0.000 0.000
CALE850919 8.1 25 Roca, Caleta de Campos, Michoacán 0.05 0.22 0.000 0.000 0.000 0.000
Clave BMDSF: EEEEAAMMDD (EEEE: estación; AA: año; MM: mes; DD: día); M: Magnitud; D: Distancia hipocentral en km;
Sd y Sa son el desplazamiento y la seudoaceleración espectral correspondiente al periodo fundamental vertical de interés
(Tv=0.1s); y R1, R2, R3 y R4 son, respectivamente, las máximas rotaciones plásticas de: VE en los extremos empotrados (VEe),
VE en el centro del claro (VEc), VA en el centro del claro (VAc), y VV en el extremo empotrado (VVe).
Es importante señalar que para los análisis y el diseño de las vigas se depreció la
contribución de la losa en el incremento de su resistencia y rigidez, y que de haberse considerado
los efectos del componente sísmico vertical tenderían a reducirse; el estudio de este último punto
sería de interés para una investigación futura.
En el inciso anterior se presentó la ecuación (8), para estimar el periodo vertical de vigas con
condiciones simples de apoyo, la cual es función de la rigidez elástica a flexión, la carga
uniforme y el claro de la viga, y de una constante adimensional función de la condición de apoyo.
La expresión propuesta depende de la constante adimensional KN, que es función de la condición
de apoyo.
Para ocho casos, donde se conoce el periodo vertical y el cociente de la rigidez a flexión
de las vigas entre la rigidez a flexión de las columnas, se calcula la constante K1, y los datos se
ajustan a una curva asintótica (figura 9, ecuación 9), la cual presenta un adecuado ajuste con los
datos (r2=0.8) y es congruente con las condiciones de ambas fronteras (viga doblemente
empotrada: EVIV/ECIC → 0, K1=22.4; viga simplemente apoyada: EVIV/ECIC → ∞, K1=9.87).
65
Tiziano Perea y Luis Esteva
K1 K1
25 25
2
K1 = π
20 (EVIV/ECIC)+22.4 20 2
K1 = π
(EVIV/ECIC)+22.4
EVIV/ECIC+1 EVIV/ECIC+1
2
15 r = 0.8 15
10 10
5 5
1E-2 1E-1 1E+0 1E+1 1E+2 1E-2 1E-1 1E+0 1E+1 1E+2 1E+3
EvIv/∑EcIc EvIv/∑EcIc
a) Datos generadores de la función de ajuste b) Comparación con otros datos estimados
Figura 9. Comparación de la constante K1 estimada de análisis y sus curvas de ajuste
π 2 (E v I v ∑ E I ) + 22.4
K1 = c c
(9)
E vIv ∑ E I +1
c c
2π(E v I v ∑ E c I c + 1) w ⋅ L4
TV1 = 2
π (E v I v ∑ E c I c ) + 22.4 g ⋅ E v I v
(10)
Iv ⋅ H
B= (12)
Ic ⋅ L
66
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
n n
2π u max
∑ mi u i g∑ mi u i ≈
2
T = 2π (13)
i =1 i =1 1.2 g
Se estudia la influencia del componente sísmico vertical en algunos marcos de concreto reforzado
con propiedades geométricas diferentes; los marcos de estudio se muestran en la figura 10: un
nivel y una crujía (1N1C), cuatro niveles y una crujía (4N1C), cuatro niveles y dos crujías
(4N2C), y diez niveles y una crujía (10N1C). Las estructuras se diseñan conforme a las prácticas
actuales, despreciando el componente vertical. Las estructuras se modelan en el Drain-2DX
considerando masa discreta en puntos espaciados a cada L/10 sobre las vigas y un
amortiguamiento ξ=5% del crítico para los modos fundamentales de las direcciones horizontal y
vertical. Se estudia su comportamiento dinámico no lineal con análisis dinámicos paso-a-paso
utilizando algunos registros seleccionados. Las historias de aceleraciones consideradas en los
análisis son, con referencia a la clave que define la Base Mexicana de Datos de Sismos Fuertes
(BMDSF), los correspondientes a los sismos enlistados en la tabla 4.
67
Tiziano Perea y Luis Esteva
Con base en los análisis dinámicos paso-a-paso de los marcos, con registros mexicanos intensos y
aplicados independiente y simultáneamente, se presentan los siguientes casos generales:
El segundo de los casos (caso B), donde el componente vertical influyó en las respuestas
globales y locales de los marcos, fue resultado de los análisis paso-a-paso con los registros
VICS800609 y VCPS870207 (figuras 11.ii, 12.ii, 13.ii, 14, 15 y 16); con estos registros, la
seudoaceleración espectral fue mayor que 1.8g (tablas 8, 10, 12 y 13), y produjo rotaciones
plásticas en el centro y extremos de algunas vigas y en los extremos de algunas columnas, con los
valores que se enlistan en las tablas 9 y 11 para los modelos 1N1C y 4N1C, respectivamente.
Con el componente vertical del resto de los registros, cuya aceleración espectral es menor
que 0.40g, el comportamiento de los marcos resultó ser elástico (caso A), y por lo tanto, las
respuestas globales y locales de los marcos, incluyendo e ignorando dicha componente son
prácticamente idénticas (figuras 11.i, 12.i y 13.i).
68
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
a) H b) ±V c) H±V
i) Comportamiento elástico en análisis dinámicos con el componente vertical (caso A)
a) H b) ±V c) H±V
ii) Comportamiento inelástico en análisis dinámicos con el componente vertical (caso B)
Figura 11. Esquematización del mapeo de rotaciones plásticas acumuladas
Tabla 10. Demanda de la intensidad espectral para el marco 4N1C (TH=0.51s, TV=0.15s)
Desplazamiento espectral (cm) Seudoaceleración espectral (g)
REGISTRO CASO NS EW V NS EW V
VICS800609 B 4.37 4.98 1.10 0.66 0.76 1.92
VCPS870207 B 4.37 7.29 1.03 0.66 1.11 1.81
IAGS791015 A 3.49 3.96 0.16 0.53 0.60 0.29
ZACA850919 A 4.03 2.91 0.15 0.61 0.44 0.27
BALC941210 A 0.84 0.65 0.12 0.13 0.10 0.21
a) H b) ±V c) H±V
i) Comportamiento elástico en análisis dinámicos con el componente vertical (caso A)
a) H b) ±V c) H±V
ii) Comportamiento inelástico en análisis dinámicos con el componente vertical (caso B)
Figura 12. Esquematización del mapeo de rotaciones plásticas acumuladas
69
Tiziano Perea y Luis Esteva
Tabla 11. Resumen de rotaciones plásticas máximas y acumuladas obtenidas del análisis del
marco 4N1C con el registro VICS800609 y esquematizadas en la figura 12.ii
ROTACIONES MÁXIMAS EN EL CENTRO DEL CLARO DE VIGAS
Nivel HNS HEW +V -V HNS+V HNS-V HEW+V HEW-V
4 - - -0.014 -0.013 -0.017 -0.015 -0.017 -0.014
3 - - -0.012 -0.011 -0.014 -0.013 -0.013 -0.012
2 - - -0.010 -0.010 -0.012 -0.012 -0.010 -0.010
1 - - -0.007 -0.007 -0.009 -0.008 -0.007 -0.007
ROTACIONES MÁXIMAS EN LOS EXTREMOS DE VIGAS
Nivel HNS HEW +V -V HNS+V HNS-V HEW+V HEW-V
4 - - - - - - - -
3 0.006 0.009 0.013 0.012 0.028 0.026 0.025 0.028
2 0.004 0.004 0.004 0.004 0.016 0.014 0.014 0.017
1 0.012 0.016 - - 0.014 0.011 0.014 0.013
ROTACIONES MÁXIMAS EN COLUMNAS
POSITIVAS NEGATIVAS
N +V -V H H+V H-V H H+V H-V
4s 0.018 0.015 0.004 0.031 0.028 -0.002 -0.021 -0.020
4i - - - 0.001 - - - -0.001
3s - - - - - - - -
3i 0.003 0.004 0.002 0.006 0.005 - -0.007 -0.008
2s 0.005 0.005 0.002 0.015 0.015 -0.002 -0.005 -0.006
2i 0.007 0.007 0.004 0.009 0.012 -0.005 -0.022 -0.024
1s 0.005 0.005 0.003 0.011 0.016 -0.005 -0.012 -0.012
1i - - 0.005 0.005 0.005 -0.005 -0.005 -0.006
ROTACIONES ACUMULADAS EN COLUMNAS
POSITIVAS NEGATIVAS
N +V -V H H+V H-V H H+V H-V
4s 0.018 0.015 0.004 0.031 0.028 -0.002 -0.021 -0.020
4i - - - 0.001 - - - -0.001
3s - - - - - - - -
3i 0.003 0.004 0.002 0.006 0.005 - -0.007 -0.008
2s 0.005 0.005 0.002 0.015 0.015 -0.002 -0.005 -0.006
2i 0.007 0.007 0.004 0.009 0.012 -0.005 -0.022 -0.024
1s 0.005 0.005 0.003 0.011 0.016 -0.005 -0.012 -0.012
1i - - 0.011 0.007 0.008 -0.015 -0.012 -0.014
Nivel Nivel
4 4
H+V
3 3
2 2
1 1
H+V
0 0
-0.012 -0.006 0.000 0.006 0.012 -0.001 -0.001 0.000 0.001 0.001
70
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Nivel Nivel
4 4
H
H+V
H-V
3 3
2 2
1 1
H
H+V
H-V
0 0
-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 0.08 -0.010 -0.006 -0.002 0.002 0.006
Tabla 12. Demanda de la intensidad espectral para el marco 4N2C (TH=0.64s, TV=0.13s)
Desplazamiento espectral (cm) Seudoaceleración espectral (g)
REGISTRO CASO NS EW V NS EW V
VICS800609 B 6.56 8.34 0.99 0.63 0.80 2.32
VCPS870207 B 8.08 13.40 0.86 0.78 1.29 2.01
IAGS791015 A 4.50 4.29 0.15 0.43 0.41 0.34
a) H b) ±V c) H±V
i) Comportamiento con el componente vertical del registro VCPS870207 (Caso B)
a) H b) ±V c) H±V
ii) Comportamiento con el componente vertical del registro VICS800609 (Caso B)
Figura 14. Mapeo de rotaciones plásticas acumuladas del modelo 4N2C
71
Tiziano Perea y Luis Esteva
Tabla 13. Demanda de la intensidad espectral para el marco 10N1C (TH=1.0s, TV=0.12s)
Desplazamiento espectral (cm) Seudoaceleración espectral (g)
REGISTRO CASO NS EW V NS EW V
VICS800609 B 10.22 9.19 1.04 0.40 0.36 2.84
VCPS870207 B 10.90 14.52 0.74 0.43 0.57 2.01
IAGS791015 A 6.69 6.79 0.13 0.26 0.27 0.37
a) +HEW b) -V c) +HEW-V
i) Comportamiento con el componente vertical del registro VCPS870207 (Caso B)
a) +HNS b) +V c) +HNS+V
ii) Comportamiento con el componente vertical del registro VICS800609 (Caso B)
Figura 15. Mapeo de rotaciones plásticas acumuladas del modelo 10N1C
Nivel Nivel
10 10
9 9
8 8
7 7
6 6
5 5
4 4
3 3
2 2
1 1
0 0
-0.10 -0.06 -0.02 0.02 0.06 0.10 0.14 -0.005 -0.003 -0.001 0.001 0.003 0.005 0.007
δi (m) δi / H
Figura 16. Desplazamientos máximos de piso y distorsiones de los análisis unidireccional (H,
líneas continuas) y bidireccional (H+V, líneas discontinuas) del modelo 10N1C con
el registro VICS800609 (Caso B)
72
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
Las plastificaciones en los puntos medios y en los extremos de las vigas tienden a crecer para
los elementos de mayor periodo en el sentido vertical, que generalmente son las vigas de
niveles superiores, las que tienen la condición de apoyo más desfavorable. En el centro del
claro se presentan rotaciones plásticas generadas por un incremento en momentos positivos;
las articulaciones en los extremos de las vigas están asociadas a un incremento en el momento
negativo de flexión. Para proporcionar ductilidad a las secciones de mayor demanda a flexión,
será necesario dotar con una cuantía suficiente de refuerzo transversal y/o acero de
compresión, incluso para la sección central donde frecuentemente se detalla con la mayor
separación de los estribos y el mínimo refuerzo de compresión. Una inversión en los
momentos de flexión no ocurrió, y por lo tanto tampoco la plastificación de las vigas, por lo
que pareciera ser suficiente la cuantía mínima especificada en la norma para momento
negativo en el centro del claro y para momento positivo en los extremos; sin embargo, no así
para dotar de ductilidad en dichas secciones al trabajar como acero en compresión; se aclara
que en los análisis no se considera la contribución de la losa en la resistencia, rigidez y
ductilidad de las vigas.
Las plastificaciones en las columnas tienden a ser más grandes en la conexión con la viga de
mayor periodo, que generalmente corresponde a la viga de azotea, por tener la condición de
apoyo más desfavorable. Con respecto a lo anterior, la tendencia en las columnas de pisos
superiores es de un incremento en la demanda a flexión (hasta un 30%) mientras que en las
inferiores el incremento es para la demanda en carga axial (hasta un 60%). El caso del
incremento en la demanda a flexión implica un aumento en el nivel de rotaciones plásticas,
mientras que el incremento de la carga axial implica una disminución en la capacidad de
rotación; ambos efectos tienden a producir un mecanismo de falla en las columnas por
aplastamiento del concreto, lo cuál sería más grave si fuera escaso el refuerzo transversal. Se
sugiere encauzar futuras investigaciones que evalúen con detalle los efectos de la degradación
de rigidez y resistencia en los elementos.
73
Tiziano Perea y Luis Esteva
Si la estimación de los efectos del componente vertical fuera una aplicación reglamentaria, las
dimensiones al menos en las vigas tendería a incrementar, buscando disminuir sus daños. Sin
embargo, a fin de no tener un mecanismo indeseable de falla, también se debería cumplir en
principio con la propuesta reglamentaria vigente que exige que el momento resistente
factorizado en las columnas sea mayor que el de las vigas, a fin de buscar el mecanismo
columna fuerte viga débil. Con lo anterior, se buscará no sólo evitar un mecanismo
indeseable, sino también considerar las incertidumbres en las variaciones de la demanda a
flexocompresión inducida por el componente vertical de movimientos fuertes del terreno.
Aumentar la inercia de las vigas o alguna zona de ellas (vigas de sección variable). Con ello
se busca incrementar la resistencia y rigidez de la viga para disminuir las rotaciones plásticas
en la sección de mayor demanda. Lo anterior también influirá en una reducción del periodo
fundamental vertical, lo cual no implica necesariamente en una reducción de la fuerza inercial
vertical. Para el caso de vigas que formen parte de marcos, se deberá revisar la interacción
con las columnas a fin de buscar el mecanismo columna fuerte - viga débil.
Las sugerencias anteriores son aplicables también en los casos en que se desee reforzar
vigas de marcos ubicados en sitios donde el componente sísmico vertical es importante.
74
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
CONCLUSIONES
Del análisis de los registros mexicanos, se observa que la relación entre el componente vertical y
el horizontal presenta grandes dispersiones; sin embargo, es claro que su amplitud es función de
la distancia epicentral, la profundidad focal y el periodo de oscilación.
Llama la atención que los espectros de respuesta (ξ=5%) correspondientes a los registros
con el componente vertical más intenso corresponden a las estaciones Islas Agrarias (IAGS),
Victoria (VICS) y Volcán Cerro Prieto (VCPS), registros que de acuerdo con la tabla 4 presentan
las distancias epicentrales y profundidades focales más cortas, y cuyas estaciones se ubican en el
Valle de Mexicali - Baja California Norte. Futuros trabajos deberán encaminarse a estudiar con
detalle implicaciones en la respuesta dinámica del suelo y la estructura con registros similares, es
decir, registros de sismos superficiales cercanos a la fuente.
Futuros trabajos deberán encaminarse a obtener, en forma directa y para algún sitio
específico, el espectro elástico de diseño vertical, así como al estudio de algunos conceptos que
aún no quedan bien definidos: diferencias en las características dinámicas estructurales para
vibración horizontal y vertical (amortiguamiento, disipación por comportamiento inelástico), etc.
Por otro lado, también se estudió la influencia del componente sísmico vertical en la
respuesta no lineal de vigas y marcos de concreto reforzado. Para los análisis dinámicos paso-a-
paso se consideraron un conjunto de registros ligados a sismos originados en la zona de
subducción de las Costas del Pacífico de México.
75
Tiziano Perea y Luis Esteva
Con base en los análisis dinámicos paso-a-paso de las estructuras en estudio, y con
registros mexicanos de movimientos fuertes del terreno aplicados independiente y
simultáneamente, se presentan las siguientes observaciones:
De acuerdo con los resultados mostrados, las máximas rotaciones plásticas en los
elementos estructurales se obtuvieron con los registros de las mayores intensidades verticales, es
decir, con los registros que presentan las mayores aceleraciones espectrales asociadas al periodo
de vibración vertical de los modos correspondientes. Estos registros corresponden a los obtenidos
en el valle de Mexicali, y cuya característica principal es que están relacionados con sismos
superficiales y con estaciones cercanas al epicentro.
76
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
REFERENCIAS
77
Tiziano Perea y Luis Esteva
78
Componente vertical de registros sísmicos en México y su efecto en la respuesta sísmica no lineal de edificios
79