Revis ta Ingenieras Universidad de Medelln 5(9): 59-74 julio-diciembre de 2006

Deficiencias, limitaciones, ventajas y desventajas

de las metodologas de anlisis ssmico no lineal
Mauricio Alexander Mora*, Jess Daniel Villalba** y Esperanza Maldonado***

Recibido: 27/07/2006
Aceptado: 27/09/2006

RESUMEN
El presente trabajo est dirigido a la comparacin de las metodologas de anlisis ssmico: dinmico
cronolgico no lineal, pushover tradicional, y dos procedimientos modales. La comparacin se realiz
a partir del estudio de un edificio de mediana altura, aporticado, de concreto reforzado, representativo
de las edificaciones tpicamente construidas en Bucaramanga. La estructura fue sometida a uno de los
sismos definidos por el estudio de microzonificacin de la ciudad, escalado a diferentes valores de ace-
leracin pico efectiva. Los resultados obtenidos fueron comparados en trminos de desplazamientos y
derivas, con lo cual se pudieron detectar algunas deficiencias, limitaciones, ventajas y desventajas de las
metodologas utilizadas, concluyendo que las metodologas modales son una buena aproximacin para
la determinacin del comportamiento ssmico de estructuras de concreto reforzado.
Palabras clave: anlisis ssmico inelstico, anlisis dinmico cronolgico no lineal, anlisis esttico no
lineal, pushover tradicional, pushover modal.

ABSTRACT
This research presents the comparison of different non linear analysis methodologies: the non linear
response history analysis, the traditional pushover methodology, and two modal pushover methodolo-
gies. A medium rise reinforced concrete building, with moment resisting space frame, representative of
the typical buildings in Bucaramanga, is examined using one of the earthquakes in the microzonation
study of the city, which was scaled to different values of effective peak acceleration. The comparison of
the results (floor displacements and interstory drifts) permitted to find deficiencies, limitations, advan-
tages and disadvantages of the methodologies used. It was found that Modal Methodologies are a good
approach to assess the seismic performance of reinforced concrete buildings.

* Ingeniero Civil. Investigador grupo INME, Universidad Industrial de Santander. E-mail: [email protected]

** Ingeniero Civil, Investigador grupo INME, Universidad Industrial de Santander. E-mail: [email protected]

*** Doctora en Ingeniera de Caminos, Canales y Puertos. Profesora Titular Escuela de Ingeniera Civil, Universidad Industrial de
Santander. E-mail: [email protected]
60 Mora et al.

INTRODUCCIN cacin de la metodologa pushover modal, en la cual

se considera la participacin de varios modos.
El modelamiento del comportamiento de las
estructuras ante eventos ssmicos es complejo y El principal objetivo de esta investigacin es pre-
la comunidad cientfica an no ha definido un sentar la metodologa de anlisis modal, como una
nico procedimiento para ello. alternativa al riguroso anlisis dinmico cronolgi-
co no lineal, conservando la facilidad de aplicacin
En la actualidad este fenmeno puede ser modela- de un anlisis pushover tradicional y superando
do de dos formas. La primera como un evento de algunas de sus deficiencias y limitaciones. Lo ante-
naturaleza dinmica, para lo que se utiliza Anlisis rior se lograr mediante la aplicacin de cada una
Dinmico Cronolgico No Lineal (Non Linear Res- de las anteriores metodologas y la comparacin
ponse History Analysis, NLRHA), el cual permite de los resultados obtenidos en trminos de des-
conocer la variacin de cualquier respuesta de la plazamientos y derivas de piso.
edificacin (desplazamientos de piso, derivas, fuerza
cortante, etc.) en el tiempo. Asimismo, posee una
DESCRIPCIN
base fsica rigurosa y es complicado de realizar. De
otra parte, el efecto del sismo en las edificaciones DE LAS METODOLOGAS
puede modelarse como un evento esttico, median-
Anlisis dinmico cronolgico no lineal
te la utilizacin del procedimiento conocido como
Anlisis Esttico No Lineal o Pushover, descrito en
(NLRHA, Nonlinear Response History
algunos documentos como FEMA-273 o ATC-40. A Analysis)
pesar de tener algunas limitaciones desde el punto
La esencia de la accin ssmica es la naturaleza di-
de vista metodolgico y del alcance y exactitud de
nmica de la carga, lo que conlleva a la variabilidad
sus resultados, este procedimiento ha tenido una
en el tiempo de las fuerzas internas del sistema y
mayor aceptacin.
su estado (definido en trminos de desplazamien-
El Anlisis Pushover Tradicional (FEMA-273, ATC- tos y deformaciones). Por lo tanto, mientras un
40) slo tiene en cuenta el aporte del primer modo problema esttico tiene una nica solucin inde-
de vibracin, considerado imperante por poseer, pendiente del tiempo, la solucin de un problema
generalmente, el mayor factor de participacin dinmico requiere la descripcin del sistema en
modal. Un aspecto importante en la aplicacin de todos los instantes dentro del perodo en estudio,
esta metodologa es la construccin de la curva de lo cual conlleva a grandes esfuerzos analticos y
capacidad de la edificacin. Esta curva se obtiene numricos que pueden, a su vez, implicar el uso
empujando el edificio horizontalmente mediante de una herramienta computacional.
una distribucin de carga predefinida, la cual se in-
De otra parte, el anlisis debe estar basado en una
crementa paso a paso hasta alcanzar un valor mximo
apropiada representacin del movimiento del terreno
de desplazamiento en el ltimo piso, o el colapso de
y en principios aceptados por la dinmica estructu-
la estructura. Los documentos gua no especifican
ral. Para sistemas de mltiples grados de libertad la
un tipo de distribucin de carga a utilizar, sino que
ecuacin dinmica se puede expresar como:
permiten cierta libertad. Lo anterior se convierte en
una de las principales deficiencias de la metodologa, Ec. 1
dado que la curva de capacidad del edificio es sensible
a la distribucin de carga utilizada. En donde M, C y K son las matrices de masa, amor-
tiguamiento y rigidez, respectivamente, U es el vector
Una posible solucin a las limitaciones de la tcnica de desplazamientos, U es el vector de velocidad,
pushover tradicional, antes mencionadas, es la apli- es el vector de aceleraciones, g(t) la historia de ace-

Revista Ingenieras Universidad de Medelln

 Deficiencias, limitaciones, ventajas y desventajas ... 61
leraciones en el tiempo en la base del edificio, y cada estos servirn para compararlos con los obtenidos
elemento del vector de influencia es igual a 1. por los mtodos de anlisis no lineal.
El vector U est relacionado con el nmero de grados Anlisis esttico no lineal-pushover tradicional
de libertad o posibilidades de movimiento de cada
una de las masas concentradas. Generalmente, las El Anlisis esttico no lineal o anlisis pushover
masas se concentran en cada piso y se supone que no es otro tipo de metodologa para determinar la
se pueden desplazar ni vertical ni rotacionalmente, respuesta no lineal de la estructura. A diferencia
permitiendo solo que se desplacen horizontalmente. del anlisis dinmico cronolgico no lineal, en
Con estas suposiciones el sistema se reduce a dos la metodologa pushover se calculan los valores
posibilidades de movimiento o grados de libertad mximos de cualquier respuesta de la edificacin,
por masa: uno en el eje x y otro en el eje y. y no su historia en el tiempo.
Si el anterior sistema de ecuaciones (Ec. 1) se plan- El anlisis tiene como base dos conceptos pri-
tea en una direccin, se tendr slo un grado de mordiales: la capacidad y la demanda. La primera
libertad para cada masa, y como stas se concentran puede entenderse como una caracterstica propia
en cada piso, el total de desplazamientos es igual a de la edificacin que depende de factores como la
N, donde N es el nmero de pisos; por consiguien- geometra de los elementos, la cantidad de refuer-
te, el vector U tiene tantos elementos como pisos zo, las propiedades de los materiales como rigidez
tiene la edificacin. Igual sucede con las matrices y ductilidad, entre otras. La demanda depende de
M, C y K, las cuales tienen N x N elementos. un sismo en particular, representado mediante
La solucin de este sistema de ecuaciones es parti- un acelerograma o un espectro de diseo, y hace
cularmente complicada porque, si bien, la matriz referencia a las fuerzas y deformaciones impuestas
de masa de la edificacin, M, es constante en el por ste. De manera que la demanda, a diferencia
tiempo, las matrices de amortiguamiento y rigidez, de la capacidad, no es en teora un valor constante,
C y K, respectivamente, cambian cuando se supera pues depende del conjunto de fuerzas externas o
el lmite de comportamiento elstico. aceleraciones a las que se somete la edificacin.

Sin embargo, reconstruir las matrices de rigidez y Para evaluar la capacidad de la edificacin se debe
amortiguamiento cada vez que un elemento (viga, construir la curva de capacidad, la cual se obtiene
columna, muro) supera la fluencia no es el nico graficando los valores mximos de cortante basal y
problema. La parte derecha de la ecuacin (1) re- desplazamiento en el ltimo piso. El procedimiento
presenta las fuerzas externas a las que se somete la para construir la curva de capacidad consiste en
edificacin, siendo g(t) la aceleracin en la base aplicar cargas laterales monotnicas, las cuales se van
del edificio en el tiempo, es decir, el acelerograma, incrementando hasta llevar la estructura al colapso.
de manera que la ecuacin de movimiento es dife- En la figura 1 se muestra un esquema de las fuerzas
rente para cada valor de aceleracin. que se aplican por piso a la estructura; la sumatoria
de todas las fuerzas Fi es igual al cortante basal V.
En el presente estudio, el Anlisis Dinmico Crono-
lgico No Lineal (NLRHA) se lleva a cabo utilizando La edificacin responde a la aplicacin de las cargas,
el software IDARC 2D versin 6 (Inelastic Damage desplazndose lateralmente hasta que alcanza un mxi-
Analysis of Reinforced Concrete Structures; State mo umax. El par de valores u y V representa un punto
University of New York at Buffalo, Department of en la curva de capacidad de la edificacin, de manera
Civil Engineering, Julio 2004) y se asume que los que para determinar otros puntos se carga la estructura
resultados obtenidos mediante este tipo de anlisis nuevamente con fuerzas mayores y, as, hasta generar
corresponden al punto de referencia, es decir, que un mecanismo de colapso en la estructura.

5 (9):59-74 Julio - diciembre de 2006

62 Mora et al.

Figura 1. Procedimiento pushover para determinar la curva de capacidad de la edificacin

La demanda vara dependiendo del sismo al que modo; posteriormente, todo el sistema de N grados
se somete la edificacin y se evala mediante de libertad (MDF) se reduce a N sistemas de un
espectros de respuesta (o diseo). El concepto solo grado de libertad (SDF), que se someten a la
de espectro de respuesta aparece como un medio accin ssmica, para, finalmente, combinar todos
conveniente para resumir la respuesta mxima de los aportes modales y encontrar una aproximacin
sistemas lineales y no lineales de un grado de liber- a la respuesta total.
tad (SDF) a un movimiento ssmico en particular.
Igualmente, permite aplicar el conocimiento de El procedimiento tambin se presenta como una
la dinmica estructural al diseo de estructuras y solucin a las deficiencias encontradas en el An-
al desarrollo de los requerimientos de resistencia lisis pushover tradicional, ya que se define una
a fuerzas ssmicas consignados en los cdigos sola distribucin de carga horizontal y toma en
(Chopra, 2001a). cuenta el aporte de todos los modos de vibracin
a la respuesta total de la estructura.
En la literatura pueden encontrarse artculos de
varios expertos, relacionados con este tipo de an- Inicialmente, para aplicar este mtodo, es nece-
lisis, tales como Krawinkler y Seneviratna (1998), sario desarrollar un procedimiento de Anlisis
Mwafy y Elnashai (2001), Inel et al (2003), Lpez modal desacoplado de la historia de respuesta
(2004), Sierra (2004), entre otros. (UMRHA), que es una buena aproximacin al
Anlisis dinmico cronolgico no lineal (NLRHA)
o Anlisis dinmico, para sistemas de mltiples
Anlisis no lineal modal grados de libertad (MDOF).
En el ao 2001, Anil K. Chopra y Rakesh K. Goel
publicaron el trabajo titulado A modal pusho- Anlisis modal desacoplado de la historia
ver procedure to estimate seismic demands for de respuesta (UMRHA, Uncoupled Modal
buildings: Theory and preliminary evaluation, Response History Analysis)
el cual es una combinacin de las metodologas
presentadas anteriormente (Anlisis dinmico La respuesta inelstica de un edificio de N-pisos
cronolgico no lineal y esttico no lineal). Ini- con simetra en planta en los dos ejes ortogonales x
cialmente, se implementa un Anlisis esttico no y y, cuando se somete a un movimiento ssmico en
lineal (pushover) para evaluar la capacidad por la direccin de uno de sus ejes de simetra (x o y),

Revista Ingenieras Universidad de Medelln

 Deficiencias, limitaciones, ventajas y desventajas ... 63
se puede determinar como funcin del tiempo me- Finalmente, se combinan las respuestas modales
diante el anlisis UMRHA (Chopra y Goel, 2001). para determinar la total utilizando las siguientes
El procedimiento a seguir parte de determinar la ecuaciones:
frecuencia natural n y los modos de vibracin n
bajo comportamiento lineal; para el i-simo modo N N
se elabora la curva de capacidad utilizando una sola u (t ) = u n (t ) = n n Dn (t )
forma de distribucin de carga sn*, que depende de n =1 n =1
la matriz de masa y de la forma modal. La curva
obtenida es idealizada como una relacin bilineal y
luego convertida a la relacin Fsn/Ln Dn, utilizando
las siguientes ecuaciones: Ec. 4

Ec. 2 A partir del resultado anterior se puede determinar
el valor pico, ro, de la respuesta total
Siendo Vbny el cortante en el punto de fluencia;
M n* la masa modal efectiva; Dny es el desplaza- Anlisis pushover modal
miento de fluencia del sistema de un grado de (MPA, Modal Pushover Analysis)
libertad equivalente del MDF para el modo n;
u rny el desplazamiento de fluencia en el ltimo Chopra y Goel (2001), propusieron el Anlisis
piso; Fsn /Ln es la fuerza resistente; n el factor de pushover modal (MPA) como una herramienta
participacin modal. para determinar la respuesta inelstica pico de un
edificio sometido a un sismo.
A su vez, se calcula la historia de deformacin, Dn(t),
y la seudo-aceleracin, An(t), para el i-simo modo En este mtodo se desarrollan los primeros pasos
de vibracin del sistema inelstico de un grado del anlisis modal desacoplado (UMRHA) hasta
de libertad asumiendo que la relacin fuerza-de- convertir la curva de capacidad bilineal en la
formacin es la relacin Fsn /Ln Dn. Con ello, se relacin Fsn/Ln Dn. Luego se determina la defor-
determina la historia de varios tipos de respuesta macin pico, Dn, del i-simo modo de vibracin del
como derivas de piso, fuerza cortante en la base, sistema inelstico de un grado de libertad (SDF)
etctera, mediante las ecuaciones: utilizando la relacin fuerza-deformacin. Luego,
se calcula el desplazamiento mximo en el ltimo
Ec. 3 piso urno asociado al i-simo modo de vibracin del
En donde un(t) representa la historia en el tiempo de sistema inelstico SDF utilizando la ecuacin:
los desplazamientos de piso y rn(t) es la variacin en u rno = n rn D n Ec. 5
el tiempo de cualquier tipo de respuesta, teniendo
en cuenta que r stn corresponde a la respuesta modal Cuando se alcance la deformacin urno, se extrae
esttica, es decir, el valor esttico de r (respuesta) de la base de datos pushover los valores de otras
debido a las fuerzas externas s*n. respuestas, rno. Se repite lo anterior para tantos
modos como sean necesarios, hasta obtener sufi-
El procedimiento anterior se realiza para tantos ciente exactitud. Al igual que para la metodologa
modos como sean necesario para obtener una UMRHA la utilizacin de los dos o tres primeros
buena exactitud. Por lo general, los primeros dos modos es suficiente, por lo general.
o tres modos son suficientes.

5 (9):59-74 Julio - diciembre de 2006

64 Mora et al.

Finalmente, se determina la respuesta total utili- por prticos ortogonales entre s. Cuenta con una
zando la regla de combinacin SRSS (raz cuadrada planta de 11.6 m 10.5 m, y una altura de 17.7 m.
de la suma de cuadrados): En la figura 2 se presenta una vista en planta de los
niveles 1-6 y la elevacin de los prticos en cada
direccin. Los prticos son iguales en cada direc-
cin: 3 luces de 4.4, 2.45 y 4.75 m en la direccin
Ec. 6 este-oeste (E-W: direccin x) y 2 luces y 2 voladizos
exteriores de 4.15 y 1.1 m respectivamente en la
direccin norte-sur (N-S: direccin y).
APLICACIN
DE LAS METODOLOGAS El sistema de entrepiso est conformado por placas
rgidas aligeradas apoyadas en la direccin este-oeste
Seleccin de la estructura (E-W). La planta del nivel 7 es igual a la de los dems
niveles, excepto por los voladizos. Las columnas son
La estructura utilizada corresponde a un edificio de seccin cuadrada, de 0.50.5 m y las vigas rectan-
de concreto reforzado, de 7 pisos, regular en gulares de 0.350.45 m. La resistencia del concreto
planta y altura, representativo de las edificaciones utilizado en placas, vigas y columnas es de 21 MPa,
construidas en el municipio de Bucaramanga. El y la del acero de refuerzo tanto longitudinal como
sistema de resistencia ssmica est conformado transversal es de 420 MPa (Mora y Villalba, 2005).

Figura 2. (a) Vista en planta de los pisos 1-6; (b) Prticos A, B, C; (c) Prticos 1, 2, 3 y 4; y Refuerzo
longitudinal y transversal: (d) columnas y (e) vigas. Todas las medidas en metros

Revista Ingenieras Universidad de Medelln

 Deficiencias, limitaciones, ventajas y desventajas ... 65
Definicin de la accin ssmica
Para la aplicacin de las metodologas se utiliz el
sismo de Coalinga, el cual garantiza que la estruc-
tura se comporte inelsticamente para el intervalo
de aceleraciones que se va a trabajar. El registro
fue tomado en la estacin USGS 1608 Oil Fields
Fire Station a las 02:39 del 22 de julio de 1983,
present una magnitud de 5.77, profundidad de 7.4
km y aceleracin pico efectiva de 0.213g. Este sismo
es uno de los definido por INGEOMINAS para
la realizacin del Estudio de microzonificacin
ssmica del rea metropolitana de Bucaramanga
(Ingeominas, 2001). En la figura 3 se presenta el
registro de los 16 primeros segundos del sismo,
dentro de los cuales se present la aceleracin
pico efectiva.
0.3

0.2
Aceleracin del terreno (g)

0.1

0
0 2 4 6 8 10 12 14 16

-0.1

-0.2

-0.3
Tiempo (seg)

Figura 4. Distribucin lateral de fuerzas en

Figura 3. Sismo Coalinga la construccin de la curva de capacidad. (a)
Distribucin uniforme; (b) Distribucin ELF; (c)
Anlisis esttico no lineal- pushover Distribucin triangular invertida o lineal
tradicional
Para la construccin de las curvas de capacidad
Con el fin de observar la influencia de la distri- se utiliz el software IDARC 2D; los resultados
bucin lateral de fuerzas en la construccin de la se presentan en la figura 5. Estos muestran que
curva de capacidad, sta se calcul para la direccin existe una considerable variacin en funcin de
x del edificio en estudio, tomando una distribucin la distribucin de carga horizontal utilizada. Las
uniforme, una distribucin ELF (equivalent lateral distribuciones ELF y triangular invertida son muy
force) y una triangular invertida, con los valores similares, pero difieren de la distribucin uniforme,
que se muestran en la figura 4. incluso en la rigidez inicial.

5 (9):59-74 Julio - diciembre de 2006

66 Mora et al.

2400 se presenta el punto de comportamiento para cada

2000
distribucin de carga y el espectro elstico para el
Distribucin Triangular sismo de Coalinga.
CORTANTE BASAL [KN] .

Distribucin Uniforme
1600

Distribucin ELF Una vez conocido el punto de comportamiento es

1200
posible determinar casi cualquier valor de respuesta
800
mxima de la edificacin. En la figura 7(a) puede ver-
se el perfil de desplazamientos de piso mximos como
400
porcentaje de la altura de la edificacin (17.7 m),
0
obtenido mediante cada una de las tres distribuciones
0 2 4 6
DESPLAZAMIENTO LTIMO PISO [cm]
8 10 12
horizontales de carga y el perfil asumido verdadero
mediante la aplicacin del Anlisis dinmico crono-
Figura 5. Curvas de capacidad utilizando lgico no lineal. De igual forma, en la figura 7(b) se
diferentes distribuciones de carga horizontal presenta el perfil de mximas derivas de entrepiso
El siguiente paso consiste en determinar la demanda como porcentaje de la altura de cada piso.
ssmica. En primer lugar, es necesario definir num-
ricamente el movimiento ssmico, primero como
acelerograma y despus como espectro de demanda.
A manera de clculo tipo se utiliz el sismo Coalinga,
escalado uniformemente a 0.5g (figura 3).
El punto mximo hasta el que es llevada la edificacin
se conoce como punto de comportamiento, y para en-
contrarlo se han propuesto diferentes metodologas.
En esta investigacin se utiliz el procedimiento A del
ATC-40, que es bsicamente un proceso de prueba y
error: una vez se tienen en una misma grfica (forma-
to ADRS) el espectro de capacidad de la edificacin
y el espectro elstico de demanda (amortiguamiento
= 5 %), se supone un punto de comportamiento
cualquiera, luego se construye una representacin
bilineal de la capacidad, que debe tener la misma
pendiente inicial que el espectro, y la misma rea
bajo la curva hasta el punto que se asume es el punto
de comportamiento. Posteriormente, se determina
el amortiguamiento desarrollado por la estructura y
se construye el espectro de demanda para este valor
de amortiguamiento. Si el espectro de demanda
calculado intercepta al espectro de capacidad en el
punto previamente definido, entonces la suposicin
fue correcta y ese es el punto de comportamiento.
Si los puntos no coinciden, es necesario suponer un
Figura 6. Punto de comportamiento: (a)
nuevo punto de comportamiento y repetir el proceso
Distribucin uniforme (eq = 28.74);
tantas veces como sea necesario hasta que los dos
(b) Distribucin ELF (eq = 30.51) y (c)
espectros se intercepten en ese punto. En la figura 6
Distribucin triangular invertida (eq = 30.49)

Revista Ingenieras Universidad de Medelln

 Deficiencias, limitaciones, ventajas y desventajas ... 67

Figura 7. Perfiles de (a) desplazamiento de piso mximos y (b) mximas derivas de entrepiso,
normalizados por la altura total de la edificacin

Con el fin de observar la influencia de la intensidad piso mximos y en los valores pico de las derivas de
del movimiento de tierra, el sismo Coalinga se esca- entrepiso, obtenidos para la distribucin uniforme,
l para que alcanzara una aceleracin pico efectiva ELF y triangular, respectivamente. Se asume que
igual a 0.2g, 0.3g, 0.4g, 0.6g y 0.7g. Para cada valor los resultados obtenidos por el NLHRA correspon-
se construy el espectro de respuesta y se aplic el den a los ms cercanos a la realidad, por lo cual, la
mismo procedimiento descrito anteriormente. validacin y comparacin de los resultados se hizo
en funcin de ellos.
En las figuras 8 y 9 se presentan resmenes con
los porcentajes de error en los desplazamientos de

Figura 8. Error en los mximos desplazamientos de piso: (a) Distribucin uniforme; (b) Distribucin
ELF; (c) Distribucin triangular

5 (9):59-74 Julio - diciembre de 2006

68 Mora et al.

Figura 9. Error en las derivas de entrepiso mximas: (a) Distribucin uniforme; (b) Distribucin ELF;
(c) Distribucin triangular

En el caso de los desplazamientos, puede observar- representacin bilineal y, luego, al igualar el

se que no siempre la misma distribucin horizontal punto de comportamiento supuesto con el
de carga genera los mayores o menores porcentajes punto de corte de los espectros de capacidad
de error. Por ejemplo, para 0.5g, los porcentajes de (real) y demanda.
error ms grandes se obtienen para la distribucin 5. El procedimiento slo tiene en cuenta el
uniforme, pero para 0.3g se obtienen los menores. aporte del primer modo de vibracin.
Con las derivas pasa lo contrario los porcentajes de 6. La distribucin de carga horizontal seleccio-
error ms grandes, sin importar el factor de escala, nada.
se presentan para la distribucin uniforme. Las
distribuciones ELF y triangular generan errores
Anlisis esttico no lineal modal
muy similares, debido a la similitud de la curva
de capacidad obtenida, ver Figura 5. La respuesta del edificio de siete pisos se deter-
Los errores, en algunos casos mayores que el 50%, min por medio de los dos mtodos modales
se pueden relacionar con: inelsticos aproximados UMRHA (Anlisis
1. Las suposiciones inherentes a la teora utilizada modal desacoplado de la historia de respuesta)
para obtener la expresin que define el clculo y MPA (Anlisis pushover modal) y, se compar
del amortiguamiento equivalente (eq). con los resultados del riguroso procedimiento
2. La idealizacin del espectro de capacidad NLRHA utilizando el software IDARC-2D. Para
como una relacin bilineal. garantizar que la edificacin respondiera dentro
3. El espectro de respuesta se obtiene para del rango inelstico, inicialmente se utiliz el
un sistema de un grado de libertad (SDF), acelerograma del sismo Coalinga escalado a 0.6g.
y luego se aplica a un sistema de mltiples En la figura 10 se presentan las propiedades di-
grados de libertad (MDF). nmicas de los tres primeros modos de vibracin
4. La suma de errores, primero al igualar las del edificio en la direccin x y la distribucin
reas bajo el espectro de capacidad real y la lateral de fuerzas:

Revista Ingenieras Universidad de Medelln

 Deficiencias, limitaciones, ventajas y desventajas ... 69

Figura 10. (a) Forma y periodo de los tres primeros modos

y (b) Distribucin de fuerzas horizontales, en KN

Las curvas de capacidad, para cada uno de los tres piso, o cualquier otra respuesta de la edificacin.
primeros modos de vibracin obtenidas utilizando En la figura 12 se presenta la variacin del despla-
la distribucin lateral de fuerzas de la figura 10, se zamiento en el primer y ltimo piso, obtenida me-
presentan en la figura 11. diante el anlisis UMRHA y mediante el NLRHA.
Puede verse que el anlisis UMRHA constituye una
Como se mencion anteriormente, el anlisis buena aproximacin y, en algunos casos, puede
UMRHA permite determinar una aproximacin a utilizarse como alternativa, pues es ms fcil de
la historia en el tiempo de los desplazamientos de implementar que el anlisis NLRHA.

Figura 11. Curvas de capacidad para los tres primeros modos de vibracin (direccin x)

5 (9):59-74 Julio - diciembre de 2006

70 Mora et al.

U(t) (cm)
U(t) (cm)

Figura 12. Mtodo UMRHA vs NLRHA. (a) Historia de desplazamientos del piso 7 y (b) Historia de
desplazamientos del piso 1.
En la figura 13(a) se muestra el perfil de desplazamientos de piso mximos, obtenidos mediante UMR-
HA y NLRHA. El porcentaje de error ms grande ocurre en el piso 5, en donde se alcanza un 14.5% y
el menor en el piso 2, con 3.27%.

Figura 13. Mtodo UMRHA vs NLRHA. (a) Perfil de desplazamientos mximos (Sismo Coalinga
escalado a 0.6g) y (b) Error en los desplazamientos de piso para diferentes factores de escala

Revista Ingenieras Universidad de Medelln

 Deficiencias, limitaciones, ventajas y desventajas ... 71
Con el fin de observar el cambio de los porcentajes Finalmente, los mximos desplazamientos de piso y
de error en la medida en que se vara la magnitud del mximas derivas de entrepiso del edificio aplicando
sismo, el mismo procedimiento se llev a cabo para el sismo Coalinga escalado a 0.6g en la direccin x,
el sismo Coalinga escalado a 0.1g, 0.2g, 0.3g, 0.4g, se evaluaron mediante el Anlisis pushover modal
0.5g y 0.7g. Los resultados se presentan para cada (MPA) tomando como base las curvas de capacidad
piso en la figura 13(b). Como era de esperarse, los por modo de la figura 11. En las figuras 14 y 15
errores ms pequeos se presentan para los menores se presenta el perfil de desplazamientos y derivas
factores de escala (0.1g y 0.2g) pues la estructura an mximas para cada piso, teniendo en cuenta uno,
no ha incursionado en el rango inelstico. Para fac- dos y tres modos, y el perfil calculado mediante
tores de escala intermedios (0.3g, 0.4g y 0.5g) se pre- NLRHA. Tambin se presenta el error teniendo en
sentan los porcentajes de error ms grandes debido cuenta el primer modo, los dos primeros modos y,
a que corresponden a la zona de menor correlacin finalmente, los tres primeros modos de vibracin
entre la curva de capacidad real y la idealizada, que del edificio. El anlisis de las grficas nos permite
es la que en realidad se utiliza. El error tiende a concluir que el primer modo de vibracin no es
disminuir cuando la estructura ha incursionado en suficiente cuando se quiere implementar un an-
gran medida en la fluencia (0.6g y 0.7g). lisis modal.

7 7
(a) (b)
6 1 Modo 6
5 2 Modos 5
NLRHA
Piso No
Piso No

4 3 Modos 4
3 3 3 Modos 1 Modo

2 2 2 Modos

1 1
0 0
0.00 0.50 1.00 1.50 0 5 10 15 20

Desplazamiento (% Altura) Error (% )

Figura 14. Mtodo MPA vs. NLRHA. (a) Perfil de desplazamientos mximos;
(b) % de error considerando uno, dos y tres modos
7 7
(a) 2 Modos (b)
6 6

1 Modo NLRHA
5 5
Piso No

Piso No

4 4
3 Modos 3 Modos
3 3

2 2 Modos 2
1 Modo
1 1

0 0
0.0 0.5 1.0 1.5 2.0 0 10 20 30 40 50
Deriva (% Altura de piso) Error (% )

Figura 15. Mtodo MPA vs. NLRHA. (a) Perfil de derivas mximas;
(b) % de error considerando uno, dos y tres modos

5 (9):59-74 Julio - diciembre de 2006

72 Mora et al.

El anlisis MPA se implement adicionalmente desplazamientos de piso y derivas de entrepiso.

para el sismo Coalinga escalado a 0.1g, 0.2g, Los errores relativos al anlisis NLRHA se pre-
0.3g, 0.4g, 0.5g y 0.7g. Nuevamente se tuvo en sentan en la figura 16, encontrndose conclu-
cuenta el aporte de los tres primeros modos de siones similares a las halladas para el mtodo
vibracin y se determinaron los valores pico de UMRHA.

100 100
90 90 7 6 5 4 3 2 1
Piso 7 6 5 4 3 2 1
(a) Piso (b)
80 80
70 70
Error (%)

60

Error(%)
60
50 50
40 40
30 30
20 20
10 10
0 0
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7

Factor de escala del sismo Factor de escala del sismo

Figura 16. Anlisis Pushover Modal (MPA). Error como funcin del factor de escala del sismo; (a)
Desplazamientos y (b) Derivas de entrepiso

CONCLUSIONES ral force) y triangular las curvas son prcticamente

iguales, pero existe gran diferencia, incluso en la
El Anlisis dinmico no lineal (NLRHA) es la rigidez antes de la fluencia, con la curva obtenida
metodologa de anlisis ssmico ms completa que mediante la distribucin uniforme.
existe en la actualidad, pues permite conocer la
variacin en el tiempo de cualquier respuesta de Tambin se puede concluir que la distribucin
la edificacin. Sin embargo, existen metodologas de carga que genera los errores ms grandes es la
alternativas ms utilizadas, debido principalmente distribucin uniforme. Los errores disminuyen
a que son ms sencillas de implementar, aunque utilizando las otras dos distribuciones, sin embar-
menos exactas. Para determinar el grado de go, siguen siendo mayores que aquellos obtenidos
exactitud de stas, los desplazamientos de piso y mediante los procedimientos UMRHA y MPA.
derivas de entrepiso de un edificio de concreto Las explicaciones son bsicamente dos; primero,
reforzado, aporticado, de siete pisos, regular en el procedimiento pushover tradicional slo tiene
planta y altura, sometido al sismo Coalinga, fueron en cuenta el primer modo de vibracin, y segun-
comparados contra aquellos obtenidos mediante do, la incertidumbre en torno a la seleccin de
el Anlisis dinmico no lineal, encontrndose que la distribucin horizontal de carga con la que se
los mayores porcentajes de error se presentaron construye la curva de capacidad. Puede verse que
para el procedimiento pushover tradicional. este factor tiene gran influencia en los resultados
obtenidos, y no existe un criterio unificado para
Se encontr que la construccin de la curva de seleccionar solamente una.
capacidad, necesaria para implementar un an-
lisis pushover, es muy sensible a la distribucin Los dos problemas mencionados en el prrafo an-
horizontal de carga, lo que en teora no debera terior se pueden solucionar con los procedimien-
suceder, ya que la capacidad de la edificacin es tos UMRHA y MPA. En los dos procedimientos
nica. Para las distribuciones ELF (equivalent late- modales se utiliza solamente una distribucin

Revista Ingenieras Universidad de Medelln

 Deficiencias, limitaciones, ventajas y desventajas ... 73
horizontal de carga, proporcional a la masa y forma menos correlacin entre las curvas de capacidad
del modo, para construir la curva de capacidad; de real e idealizada, esto es, cerca de la fluencia de
esta manera se obtiene slo una curva para cada la estructura. Cuando el lmite de fluencia se
modo y se elimina la incertidumbre sobre cul dis- sobrepasa significativamente, los errores tienden
tribucin utilizar, presentada en el procedimiento a disminuir.
pushover tradicional.
Se puede observar que los porcentajes de error se
Los errores obtenidos en las dos metodologas mo- reducen significativamente, especialmente en las
dales (UMRHA y MPA) aumentan cuando la edifi- derivas de entrepiso, cuando se tiene en cuenta
cacin ingresa en el rango inelstico, reflejando de el segundo modo y, en menor medida, cuando se
esta forma la dificultad e incertidumbres existentes considera el tercer modo. Sin embargo, el aporte de
para tratar de predecir el comportamiento no lineal los dems modos es insignificante; por tal razn se
de las estructuras. Las causas de dichos errores son, puede concluir que los dos o tres primeros modos
de vibracin son suficientes.
en primer lugar, la falta de una teora rigurosa,
principalmente en el anlisis UMRHA. Los dos De las metodologas alternativas de anlisis ssmico
procedimientos se basan en algunas suposiciones estudiadas, la que produce los errores ms peque-
y extensiones del anlisis elstico de edificaciones os es el anlisis MPA, por tal razn, se concluye
que no son del todo ciertas cuando el anlisis es que es la ms adecuada; no obstante, debe mencio-
inelstico. En segundo lugar, se induce al error al narse que en algunas oportunidades los porcentajes
idealizar la curva de capacidad como una relacin de error fueron grandes. Esto permite plantear la
bilineal. Los resultados permitieron determinar necesidad de mejorar los procedimientos actuales
que los mayores errores se cometen en la parte de de anlisis y diseo ssmico.

BIBLIOGRAFIA
APPLIED TECHNOLOGY COUNCIL. 1996. Seismic evaluation and retrofit of concrete buildings, Vol 1. California.

CHOPRA, K. 2001a. Dynamics of structures: Theory and applications to earthquake engineering. New Jersey. University of
California at Berkeley. Prentice Hall.

CHOPRA, K. & GOEL, R. (2001). A modal pushover procedure to estimate seismic demands for buildings: Theory and
preliminary evaluation. PEER Report, Pacific Earthquake Engineering Research Center, College of Engineering, Uni-
versity of California Berkeley.

FEMA. 1996. NEHRP Guidelines for the seismic rehabilitation of buildings. FEMA 273, Federal Emergency Management
Agency.

INEL, M., TJHIN, T. & ASCHHEIM, M (2003). The significance of lateral load pattern in pushover analysis. Proceedings
Fifth National Conference on Earthquake Engineering Istanbul, Turkey. Paper No: AE-009.

KRAWINKLER, H. & SENEVIRATNA, G. (1998). Pros and cons of a pushover analysis of seismic performance evaluation.
Engineering Structures, Vol. 20 pp. 452-464.

LOPEZ, A. 2004. A review of existing pushover methods for 2-d reinforced concrete buildings. PhD Dissertation, Rose
School.

MORA, M. & VILLALBA, J. 2005. Anlisis ssmico inelstico y clculo de ndices de dao en una edificacin de mediana
altura utilizando diferentes metodologas Tesis de Grado. Universidad Industrial de Santander.

5 (9):59-74 Julio - diciembre de 2006

74 Mora et al.

MWAFY, A. & ELNASHAI, A. (2001). Static pushover versus dynamic collapse analysis of RC buildings. Engineering Struc-
tures, Vol. 23 pp. 407424.

SIERRA, N. 2004. Evaluacin comparativa de metodologas de anlisis inelstico para edificios de concreto reforzado. Tesis
de Maestra, Universidad de los Andes.

SUBDIRECCIN DE AMENAZAS GEOAMBIENTALES, INGEOMINAS, 2001, Zonificacin sismogeotcnica indicativa

del rea metropolitana de Bucaramanga. Bogot.

Revista Ingenieras Universidad de Medelln