Analisis de Respuesta Sismica en Direcciones Ortogonales
Analisis de Respuesta Sismica en Direcciones Ortogonales
Analisis de Respuesta Sismica en Direcciones Ortogonales
Cmo citar?
Fascculo completo
Pgina de la revista
www.redalyc.org
Proyecto acadmico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
RESUMEN
En este artculo se estudia el problema que representa la combinacin de los efectos ssmicos en las
estructuras causados por la accin simultnea de los dos componentes ortogonales horizontales del
movimiento del suelo. Mediante el uso de la teora de vibraciones aleatorias se desarrollan algunas
expresiones analticas que permiten combinar los efectos ssmicos causados por la accin de ambos
componentes ortogonales horizontales del movimiento del suelo con el propsito de estimar la
mxima respuesta elstica, bi-direccional de la estructura. La principal hiptesis que se hace en el
desarrollo de estas expresiones, se refiere a la representacin del espectro de amplitudes de Fourier
de ambos componentes del movimiento del suelo mediante funciones tipo delta de Dirac. Por medio
del anlisis elstico bi-direccional, paso a paso de distintos modelos estructurales, se verifica la
precisin del planteamiento de combinacin propuesto. El procedimiento desarrollado es aplicable
exclusivamente al caso de suelo blando y considera explcitamente el ngulo de incidencia del
temblor y el tipo de respuesta que se analiza (ortogonal o colineal).
ABSTRACT
This work studies the problem that represents the combination of the seismic effects on structures
caused by the action of two horizontal orthogonal components of the ground motion. By using
random vibrations some analytical expressions were developed to combine the two horizontal
orthogonal seismic effects in order to estimate the elastic bi-directional peak response. The main
hypothesis made in the development of these expressions was to represent the Fouriers spectrum
for the two orthogonal components of the ground motion by a Diracs delta. Through the fully bidirectional elastic step by step analysis of some structural models, it is verified the accuracy of
the proposed method. The exposed procedure is exclusively applicable to soft soil. It considers
explicitly the angle of earthquake incidence and the kind of response in terms of the direction of its
components (orthogonal or collinear).
INTRODUCCIN
Durante la ocurrencia de un temblor, las estructuras se ven sometidas a un campo complejo del
movimiento del suelo. Sin embargo, en la prctica es comn analizarlas, nicamente, bajo la accin de dos
componentes ortogonales horizontales del movimiento del suelo. En algunos casos (regiones cerca del
epicentro del temblor, por ejemplo) se incluye un tercer componente ortogonal, el cual acta en la
direccin vertical.
Artculo recibido el 14 de noviembre de 2007 y aprobado para su publicacin el 20 de noviembre de 2008. Se aceptarn
comentarios y/o discusiones hasta cinco meses despus de su publicacin
(1)
Facultad de Ingeniera, Universidad Autnoma del Estado de Mxico, Cerro de Coatepec S/N, Ciudad Universitaria, Toluca,
Mxico, CP. 50130, jvaldes@uaemex.mx.
(2)
Instituto de Ingeniera, Universidad Nacional Autnoma de Mxico, Mxico, D.F, 04510, mors@pumas.iingen.unam.mx
91
Cuando el anlisis de la estructura se lleva a cabo mediante tcnicas espectrales, la mayora de los
cdigos recomiendan analizar la estructura en forma independiente en dos direcciones ortogonales para
obtener las mximas respuestas unidireccionales que interesen (mximas respuestas unidireccionales).
Bajo este criterio, el principal problema es cmo combinar las mximas respuestas unidireccionales (Rx max
y Ry max en la figura 1) para estimar la mxima respuesta bi-direccional (Rxy max en la figura 1).
Rxy max
R
t1
Sx
x
y
Sy
Rx max
Ry max
t2
Sx
t3
x
x
y
Sy
c) Respuesta ssmica considerando nicamente
la accin del movimiento del suelo en la direccin y
92
93
R xy ( t ) = R x2 ( t ) + R y2 ( t )
(1)
R xy ( t ) = R x ( t ) + R y ( t )
(2)
R x y (t)
R y (t)
R x (t)
R x (t)
R y (t)
a ) R e s p u e s ta o rto g o n a l
b ) R e s p u e s ta c o lin e a l
R xy ( t )max R xy ( t )max
,
R x ( t )max R y ( t )max
= Min
(3)
94
R y ( t )max
Rx ( t )max
= Min
(4)
VIBRACIONES ALEATORIAS
En el dominio de la frecuencia, cualquier respuesta de inters de una estructura sometida a la accin
simultnea de ambos componentes ortogonales horizontales del movimiento del suelo puede calcularse de
la siguiente forma: Si los componentes de la respuesta son colineales, entonces
R xy ( ) = g1 ( )&x&s ( ) + g 2 ( ) &y&s ( )
(5)
R xy ( ) = g 12 ( )&x&s2 ( ) + g 22 ( ) &y&s2 ( )
(6)
donde g1() y g2() son funciones deterministas que dependen de las caractersticas de la estructura
(rigidez, amortiguamiento, geometra, distribucin de masa, etc.), as como de la respuesta que en
particular se calcule. &x&s () y &y&s () son las transformadas de Fourier de los componentes ortogonales
horizontales de la aceleracin del suelo.
En general, Rxy() es un nmero complejo cuya amplitud al cuadrado est dada por:
2
R xy ( ) = R xy ( )R *xy ( )
(7)
(8)
(9)
95
k=
&y&s ( )
&x&s ( )
(10)
donde k vara de 0.4 a 1.0 aproximadamente y depende de la orientacin de los ejes que se utilicen para
descomponer el movimiento del suelo (Valds, 2004).
lA()l
8.00
6.00
4.00
2.00
0.00
0.00
0.50
1.00
1.50
2.00
f(hz)
Figura 3. Espectro de amplitudes de Fourier correspondiente al registro obtenido durante el temblor del 19
de septiembre de 1985 en la estacin SCT de terreno blando en la ciudad de Mxico
Bajo estas consideraciones, las ecs. 8 y 9 pueden rescribirse de la siguiente forma. Para el caso de
componentes colineales (ec. 8):
R xy ( s ) = g 12 ( s ) + k 2 g 22 ( s ) + 2 k g 1 ( s )g 2 ( s )real [ ( s )]
2
(11)
(12)
96
E R xy2 max =
Fp
2Ts
xy
( ) d
(13)
donde Fp es el factor pico y Ts la duracin de la fase intensa del movimiento del suelo.
Evaluando la ec. 13 para el caso donde los componentes de la respuesta son colineales (ec. 11)
resulta:
E R xy2 max =
F p g12 ( s ) + k 2 g 22 ( s ) + 2 k g1 ( s )g 2 ( s )real [ ( s )]
2Ts
(14)
E Rxy4 max =
(15)
Ntese que la integral de la ec. 13 ha desaparecido de las ecs. 14 y 15 debido a que los espectros de
amplitudes de Fourier del movimiento del suelo son deltas de Dirac.
S slo acta un componente del movimiento del suelo, entonces la respuesta puede calcularse
como sigue. Para la accin del temblor en la direccin x:
R x ( ) = g 1 ( )&x&s ( )
(16)
R y ( ) = g 2 ( ) k &y&s ( )
(17)
F px g12 ( s )
F py k 2 g 22 ( s )
E R x2 max =
E R y2 max =
(18)
2Tsx
(19)
2Tsy
donde Tsx y Tsy son las duraciones de las fases intensas del movimiento del suelo en cada direccin. Bajo
estas condiciones, el cociente (ec. 3) puede calcularse como sigue:
Si los componentes de la respuesta son colineales
97
c =
F1 E Rxy2 max
F2 E( Rx2max )
E( Rxy max )
E( Rx max )
(20)
c =
2Ts
F1
F2
Fpx g12 (s )
2Tsx
o =
F1 4 E R xy4 max
F2 4 E ( R x4max )
Fp
o =
F1
F2
(g
4
1
E ( R xy max )
E ( R x max )
(21)
( s ) + k 4 g 24 ( s ) + 2 k 2 g 12 ( s ) g 22 ( s )real 2 [ ( s )]
2Ts
F px g 14 ( s )
2Tsx
Las ecs. 20 y 21 fueron formuladas considerando que la respuesta individual debida a la accin del
sismo en x es mayor que la correspondiente a y (Rx mx > Ry mx).
F1 y F2 son funciones que dependen de las distribuciones de probabilidad de las variables aleatorias
que se estudian. El propsito de estas funciones es transformar la raz cuadrada de los valores esperados
de las mximas respuestas cuadrticas al valor esperado de la mxima respuesta. Para simplificar las ecs.
20 y 21, se asume que F1/F2=1. Adems, se considera que los factores pico (Fp, Fpx) y las duraciones de
las fases intensas (Ts, Tsx) son los mismos, tanto para los efectos bi-direccionales como para los
unidireccionales. De esta forma, la variable puede calcularse de acuerdo a las siguientes expresiones. Si
los componentes de la respuesta son colineales
c = 1 + 2 + 2 real [ ( s )]
(22)
o = 4 1 + 4 + 2 2 real 2 [ ( s )]
(23)
donde
k g 2 ( s )
g1( s )
(24)
98
99
c =
1 + 2 + 2 real [ ( s )] 1
(25)
o =
1 + 4 + 2 2 real 2 [ ( s )] 1
(26)
Ecuacin 22
2.0
1.8
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
1.00
0.80
0.60
0.40
0.20
0.00
0.0
0.2
0.4
0.6
0.8
1.0
1.00
0.80
1.30
0.60
0.40
1.20
0.20
0.00
1.10
1.00
0.0
0.2
0.4
0.6
0.8
1.0
100
1.00
0.6
0.80
0.4
0.60
0.2
0.40
0.0
0.20
-0.2
0.00
-0.4
-0.6
-1.0
-0.8
-0.6
-0.4
-0.2
0.0
0.2
0.4
0.6
0.8
1.0
0.80
0.60
0.6
0.40
0.20
0.4
0.00
0.2
0.0
-1.0
-0.8
-0.6
-0.4
-0.2
0.0
0.2
0.4
0.6
0.8
1.0
101
(
(
2
Fx E R x max
=
) Fy E R y2 max
E( R x max )
E( R y max
) F
=
) F
F px g12 ( s )
x
y
F py
2Tsx
k 2 g 22 ( s )
(27)
2Tsy
donde Fx y Fy son funciones que transforman la raz cuadrada del valor esperado de la mxima respuesta
cuadrtica en el valor esperado de la mxima respuesta. Si se considera que los factores pico y la duracin
de la fase intensa del movimiento del suelo son los mismos en el numerador y en el denominador, adems
de que Fx/Fy = 1, entonces
E( R x max )
E( R y max )
g1 ( s )
=
g2(s ) k
(28)
De esta forma, es el cociente de las respuestas unidireccionales pico debidas a la accin individual
del sismo en la direccin x o y. Esta conclusin es importante, ya que las respuestas individuales pico (x o
y) pueden calcularse por medio de las tcnicas usuales de espectros de respuesta..
EJEMPLOS DE APLICACIN
Se resuelven dos ejemplos para verificar la precisin del procedimiento expuesto en este trabajo. En
el primer ejemplo se analiza un modelo de tres grados de libertad (dos desplazamientos ortogonales
horizontales y un giro). En el segundo ejemplo se analiza un edificio de concreto reforzado de cuatro
niveles. El principal objetivo de estos ejemplos es comparar los resultados obtenidos utilizando el
procedimiento expuesto en este trabajo y las reglas de combinacin existentes, respecto a la respuesta bidireccional exacta, la cual corresponde a un anlisis elstico, bi-direccional, paso a paso. Los
acelerogramas que se utilizan corresponden a distintas estaciones de terreno blando de la ciudad de
Mxico.
Ejemplo 1
La figura 6 muestra el modelo de tres grados de libertad que se analiza en este ejemplo. Las
variables que definen el modelo son: a, m, Kx, Bx = ex /a, By = ey /ar, r = b/a y =Ky /Kx. Donde Kx y Ky son
las rigideces de traslacin en las direcciones ortogonales, ex y ey la excentricidad en cada direccin, a y b
las dimensiones en planta del modelo y m la masa traslacional.
102
ey K4
b/2
Ky
CR
CM
K1
Kx
K2
ex
b/2
K3
a/2
a/2
Primero, se muestra la aplicacin del procedimiento propuesto en este artculo mediante el anlisis
de un caso especial para el modelo B, en el cual By = 0.0527. Como resultado del anlisis modal
unidireccional para el acelerograma del temblor del 19 de septiembre de 1985 (Ms=8.1), obtenido en suelo
blando de la ciudad de Mxico (estacin SCT, Ts = 2.0 s), se obtienen los siguientes resultados
correspondientes a la fuerza cortante que acta sobre el marco 1: V1 = 0.604 Ton (sismo en la direccin
x) y V1 = 0.602 Ton (sismo en la direccin y). Por otra parte, los resultados obtenidos del anlisis
elstico, bi-direccional, paso a paso, son los siguientes: V1=1.02 Ton (sismo actuando en la direcciones x y
y, simultneamente), V1 = 0.757 Ton (sismo actuando en x y y), V1 = -0.757 Ton (sismo en x y y) y V1 =
1.02 Ton (sismo en x y y).
El objetivo es combinar los resultados del anlisis modal unidireccional para obtener una
estimacin de los resultados obtenidos en el anlisis elstico, bi-direccional, paso a paso. Para fines de
diseo el resultado bi-direccional que interesa es V1 = 1.02 Ton.
De acuerdo al mtodo SRSS, la mxima fuerza cortante bi-direccional que se obtiene es
2
VMAX ,SRSS = (0.604) + ( 0.602 )2 = 0.8527 Ton . Si se utiliza una regla de combinacin para = 0.30 (RDF
1993), entonces VMAX,=0.30 = 0.604 + (0.3)(0.602) = 0.7846 Ton..
La regla CQC3 es otra regla que puede utilizarse para calcular los efectos ssmicos bi-direccionales
(Menun y Der Kiureghian, 1998; Smeby y Der Kiureghian, 1985). Esta regla toma en consideracin,
explcitamente, la correlacin entre las respuestas modales, adems de la correlacin entre los
componentes del movimiento del suelo. Esta regla permite calcular la respuesta bi-direccional mxima,
como funcin del ngulo de incidencia del temblor. Para simplificar su aplicacin, se ha propuesto una
frmula simplificada que permite calcular en forma directa la mxima respuesta bi-direccional (rcr), sin
necesidad de calcular explcitamente el ngulo crtico de incidencia del temblor (Lpez y Torres, 1997;
Lpez et al., 2000). Para fines de este artculo, la regla CQC3 se define mediante la expresin simplificada
que se comenta.
103
rcr = rx 1 + k 2
) 1 + 2
2
CQC
+ 1 k 2
2
1 CQC
rxy
+
rx ry
2
CQC
(29)
donde k es el cociente del espectro de respuesta que se utilice para el anlisis en la direccin x, dividido
entre el espectro de respuesta para la direccin y, CQC = ry /rx, rx y ry son las respuestas unidireccionales
en las direcciones indicadas y rxy la covarianza entre las respuestas unidireccionales rx y ry, la cual se
calcula de la siguiente forma
(30)
donde
1
2
rk = ij rki rkj ; k = x , y
i j
(31)
Analizando el intervalo que contiene a los distintos valores de la respuesta crtica (rcr), que se
pueden calcular con la ec. 29, para todos los posibles valores de CQC, k y rxy /rx ry, se observa que:
rx rcr 2 rx , lo cual en trminos de representa una variacin de 1 2 . Los mximos valores para
la respuesta crtica bi-direccional que se pueden calcular con la ec. 29, son los mismos que se pueden
calcular mediante la regla SRSS. La regla CQC3 (ec. 29) nunca dar mayores valores para la respuesta
crtica bi-direccional que los que se obtengan con la regla SRSS. Para este caso particular, la mayor fuerza
cortante que puede calcularse aplicando la ec. 29, se obtiene para los siguientes valores de las variables
involucradas en el clculo, k = 1, CQC = 0.602/0.604 = 0.996 y rxy /rx ry = 1, los cuales hacen que:
VCQC3 = 1 + (0.9966) 2 (0.604) = 0.8527 Ton .
104
105
Error (%)
0.00
gama
alfa=0.3
-10.00
SRSS
-20.00
CQC3
-30.00
-40.00
0.00
0.10
0.20
By
0.30
0.40
0.50
a) Modelo A, marco 4
Error (%)
Error en la estimacin
15.00
10.00
5.00
0.00
-5.00
-10.00
-15.00
-20.00
-25.00
-30.00
-35.00
0.00
gama
alfa=0.3
SRSS
CQC3
0.10
0.20
By
0.30
0.40
0.50
b) Modelo B, marco 1
Error en la estimacin
10.00
5.00
Error (%)
0.00
gama
-5.00
alfa=0.3
-10.00
SRSS
-15.00
CQC3
-20.00
-25.00
0.00
0.10
0.20
0.30
0.40
0.50
By
c) Modelo C, marco 2
Figura 7. Errores en la estimacin de la mxima fuerza cortante bi-direccional que acta en algn
marco estructural de los modelos analizados en el ejemplo 1, para diferentes reglas de combinacin
106
Error en la estimacin
Modelo C, Tx=4.00 s
10.00
5.00
Error (%)
0.00
gama
alfa=0.3
-5.00
SRSS
-10.00
CQC3
-15.00
-20.00
-25.00
0.00
0.10
0.20
By
0.30
0.40
0.50
15.00
10.00
gama
Error (%)
5.00
alfa=0.3
0.00
SRSS
CQC3
-5.00
-10.00
-15.00
0.00
0.10
0.20
By
0.30
0.40
0.50
Ejemplo 2
En la figura 9 se muestra la estructura que se analiza en este ejemplo, la cual corresponde a un
edificio de concreto reforzado de cuatro niveles. Se analizan dos respuestas, el esfuerzo normal en la
columna de esquina identificada como columna A, y el cortante basal. Se utilizan los acelerogramas
registrados en la estacin SCT de la ciudad de Mxico (Ts = 2.0 s) registrados durante el temblor del 19 de
septiembre de 1985 (Ms = 8.1), los cuales se rotaron para analizar la direccin de mxima correlacin
cruzada entre los componentes del movimiento del suelo (k = 1). Se analizaron todas las posibles
combinaciones de signo para la excitacin y se consider una parte real de la coherencia de 0.6.
107
CONCLUSIONES
Se desarroll un procedimiento para estimar la mxima respuesta bi-direccional de las estructuras
desplantadas en terreno blando. Las estimaciones hechas con este procedimiento se compararon respecto a
108
Regla 100%+30%
- 8.30 %
-32.50 %
Regla
SRSS
-27.50 %
-5.00 %
-18.50 %
5.00 %
-1.00 %
- 30.00 %
-20.00 %
109
REFERENCIAS
American Society of Civil Engineers (ASCE), 1986. Seismic Analysis of Safety Nuclear Structures and
Commentary on Standard for Seismic Analysis of Safety Related Nuclear Structures, New York.
Building Seismic Safety Council. NEHRP, 1997. Recommended provisions for seismic regulations for
new buildings. FEMA302, Washington, DC.
California Department of Transportation (Caltrans), 1990. Bridge Design Specifications, Sacramento,
CA.
COVENIN, 1998. Edificaciones Sismorresistentes. Norma Venezolana COVENIN 1756.
Heredia-Zavoni, E y R Machicao-Barrionuevo (2004), Response to orthogonal components of ground
motion and assessment of percentage combination rules. Earthquake Engineering and Structural
Dynamics, Vol. 33, pp. 271-284.
International Conference of Building Officials (ICBO), 1997. Uniform Building Code, Vol. 2, Structural
Engineering Design Provisions, Whittier, CA., 492 pp
Lpez, O A y R Torres (1997), The critical angle of seismic incidence and the maximum structural
response. Earthquake Engineering and Structural Dynamics, Vol. 26, pp. 881-894.
Lpez, O A, A K Chopra, y J J Hernndez (2000), Critical response of structures to multicomponent
earthquake excitation. Earthquake Engineering and Structural Dynamics, Vol. 29(12), pp. 17591778.
Menun, C y A Der Kiureghian (1998), A replacement for the 30%, 40% and SRSS rules for
multicomponent seismic analysis. Earthquake Spectra, Vol. 14, No. 1, pp. 153-156.
Newmark, N M (1975), Seismic design criteria for structures and facilities, trans-Alaska pipeline
system. Proceedings of the U.S. National Conference Earthquake Engineering, EERI, pp. 94-103.
Nutt, R V (1996), ATC-32: Provisional recommendations for the seismic design of bridges.
Proceedings of the 4th Caltrans Seismic Workshop, California Department of Transportation,
Sacramento, CA.
Reed, J W y R P Kennedy (1996), Discussion of A Clarification of Orthogonal Effects in a TreeDimensional Seismic Analysis by E.L. Wilson, I. Suharwardy, and A. Habibullah. Earthquake
Spectra Vol. 12, No. 2, pp. 353-356.
Reglamento de Construcciones para el Distrito Federal (RDF), 1993. Departamento del Distrito Federal,
Mxico.
Reyes-Salazar, A, J A Jurez-Duarte, A Lpez-Barraza, y J I Velzquez-Dimas (2004), Combined effect
of the horizontal components of earthquakes for moment resisting steel frames. Steel and
Composite Structure An International Journal, Vol. 4, No 3, pp. 189-209.
Rosenblueth, E y H Contreras (1977). Approximate design for multicomponent earthquakes. Journal of
the Engineering Mechanics, Division (ASCE), Vol. 103, pp. 881-893.
Santa Cruz, S, E Heredia-Zavoni, y R S Harichandran (1999), Variacin de la coherencia del movimiento
ssmico del terreno en frecuencias bajas y modelos para el Valle de Mxico, Memorias del XII
Congreso Nacional de Ingeniera Ssmica, Vol. 1, Morelia- Michoacn, Mxico, pp. 197-206.
110
111